Deep borehole log evidence for fractal distribution of fractures in crystalline rock

Summary Sonic velocity and electrical resistivity logs run to a depth of 3.5 km in crystalline rock near the San Andreas fault at Cajon Pass in southern California correlate over scale-lengths both small (sub-metre) and large (tens to hundreds of metres). No such correlations are seen with the more lithologically sensitive natural gamma intensity log. The correlation between the sonic velocity and electrical resistivity logs suggests that a non-lithologic property of the crystalline rock controls fluctuations. In situ fracture intensity is a logical candidate for the controlling rock property. The fluctuations of the individual sonic velocity and electrical resistivity logs are examined with the Hurst rescaled range parameter over borehole log intervals 1.5 m < L < 1500 m. For log fluctuations arising from a scale-invariant physical process the Hurst rescaled range scales with data interval as LH, 0 < H < 1. A purely random sequence of in situ fractures produces a scaling exponent H= 0.50. Fluctuations in the Cajon Pass sonic velocity and electrical resistivity logs yield H˜ 0.70 evidence that in situ fracture sets tend to occur in clusters rather than at purely random intervals. The tendency for fracture clustering over log intervals 1.5 m < L < 1500 m suggests that fracture formation is a fractal process independent of length-scale in which larger fracture intervals form from clustering of numerous smaller fracture intervals. Seismic reflectivity derived from the borehole sonic velocity log is also scale independent over the range of data intervals 1.5m < L < 1500 m with a Hurst exponent H= 0.21. If we associate fracture clustering with crustal fault formation, the Cajon Pass borehole sonic velocity and electrical resistivity logs predict that crustal faults scale fractally with fractal dimension D˜ 2.30. The equivalent b-value for earthquake size distribution is b˜D/2˜ 1.15. On this hypothesis observed b-values < 1.15 indicate a tendency for earthquakes to cluster on existing (weak?) faults. A scale-independent mechanics for crystalline rock fracture formation in which larger scale fractures form as clusters of smaller scale fractures provides a mechanical link between the small pervasive stress aligned flaws, cracks and microfractures which can impart anisotropy to crystalline rock and the larger scale fractures associated with finite strain and faulting. Thus in crustal regions of active but low strain, fracture anisotropy is observed to be aligned with the inferred maximum principal stress, while in actively faulted crust, fracture anisotropy is observed to be more nearly fault parallel as if the fractures are aligned by the finite strain faulting process.

[1]  Ray series modeling of seismic wave travel times and amplitudes in three‐dimensional heterogeneous anisotropic crystalline rock: Borehole vertical seismic profiling seismograms from the Mojave Desert, California , 1990 .

[2]  M. Zoback,et al.  Self‐similar distribution and properties of macroscopic fractures at depth in crystalline rock in the Cajon Pass Scientific Drill Hole , 1992 .

[3]  X. R. Shih,et al.  Observation of shear wave splitting from natural events: South Moat of Long Valley Caldera, California, June 29 to August 12, 1982 , 1990 .

[4]  T. Blenkinsop Correlation of paleotectonic fracture and microfracture orientations in cores with seismic anisotropy at Cajon Pass Drill Hole, southern California , 1990 .

[5]  H. Kanamori,et al.  A moment magnitude scale , 1979 .

[6]  Keiiti Aki,et al.  Magnitude‐frequency relation for small earthquakes: A clue to the origin of ƒmax of large earthquakes , 1987 .

[7]  M. Savage,et al.  Shear wave anisotropy and stress direction in and near Long Valley Caldera, California, 1979–1988 , 1990 .

[8]  Benoit B. Mandelbrot,et al.  Fractal Geometry of Nature , 1984 .

[9]  Stuart Crampin,et al.  Seismic-wave propagation through a cracked solid: polarization as a possible dilatancy diagnostic , 1978 .

[10]  L. Frazer Dynamic elasticity of microbedded and fractured rocks , 1990 .

[11]  H. E. Hurst,et al.  Long-Term Storage Capacity of Reservoirs , 1951 .

[12]  B. Mandelbrot A Fast Fractional Gaussian Noise Generator , 1971 .

[13]  Mark W. Vincent,et al.  Laumontite mineralization in rocks exposed north of San Andreas Fault at Cajon Pass, southern California , 1988 .

[14]  B. Gutenberg,et al.  Seismicity of the Earth , 1970, Nature.

[15]  Robert W. Clayton,et al.  Finite difference simulations of seismic scattering: Implications for the propagation of short‐period seismic waves in the crust and models of crustal heterogeneity , 1986 .

[16]  P. Bak Self-organized criticality , 1990 .

[17]  K. Aki,et al.  Observation and modelling of fault-zone fracture seismic anisotropy — II. P-wave polarization anomalies , 1987 .

[18]  K. Aki,et al.  Observation and modelling of fault‐zone fracture seismic anisotropy ‐ I. P, SV and SH travel times , 1987 .

[19]  G. King,et al.  The mechanisms of finite brittle strain , 1992 .

[20]  L. Jones Focal mechanisms and the state of stress on the San Andreas Fault in southern California , 1988 .

[21]  J. B. Fletcher,et al.  Shear wave splitting in the Anza Seismic Gap, southern california: Temporal Variations as possible precursors , 1988 .

[22]  D. L. Anderson,et al.  Theoretical Basis of Some Empirical Relations in Seismology by Hiroo Kanamori And , 1975 .

[23]  P. Bak,et al.  Self-organized criticality , 1991 .

[24]  Jer-Ming Chiu,et al.  Temporal changes in shear wave splitting during an earthquake swarm in Arkansas , 1990 .

[25]  T. Hewett Fractal Distributions of Reservoir Heterogeneity and Their Influence on Fluid Transport , 1986 .

[26]  Jon B. Fletcher,et al.  Changes in shear wave splitting at Anza near the time of the North Palm Springs Earthquake , 1990 .

[27]  P. Burton,et al.  LONG-TERM EARTHQUAKE RECURRENCE CONSTRAINED BY TECTONIC SEISMIC MOMENT RELEASE RATES , 1986 .

[28]  B. Atkinson Subcritical crack propagation in rocks: theory, experimental results and applications , 1982 .

[29]  S. Crampin,et al.  Shear-wave polarizations near the North Anatolian Fault – I. Evidence for anisotropy-induced shear-wave splitting , 1985 .

[30]  Yan Sun,et al.  Oriented microfractures in Cajon Pass Drill Cores: Stress field near the San Andreas Fault , 1990 .

[31]  K. Aki,et al.  The fractal nature of the inhomogeneities in the lithosphere evidenced from seismic wave scattering , 1985 .

[32]  Geoffrey C. P. King,et al.  The accommodation of large strains in the upper lithosphere of the earth and other solids by self-similar fault systems: the geometrical origin of b-Value , 1983 .

[33]  T. Henyey,et al.  Stress orientation inferred from shear wave splitting in basement rock at Cajon Pass , 1988 .

[34]  T. Henyey,et al.  Fracture related reflectors in basement rock from vertical seismic profiling at Cajon Pass , 1988 .

[35]  J. Hudson Wave speeds and attenuation of elastic waves in material containing cracks , 1981 .

[36]  Takayuki Hirata,et al.  A correlation between the b value and the fractal dimension of earthquakes , 1989 .

[37]  R. Evans,et al.  Earthquake prediction: a new physical basis , 1984 .

[38]  P. Bak,et al.  Earthquakes as a self‐organized critical phenomenon , 1989 .

[39]  B. Mandelbrot,et al.  Fractional Brownian Motions, Fractional Noises and Applications , 1968 .

[40]  S. Crampin,et al.  Seismic fracture anisotropy in the Earth's crust: An overview , 1990 .

[41]  S. Kaneshima Origin of crustal anisotropy: Shear Wave splitting studies in Japan , 1990 .

[42]  Larry R. Myer,et al.  Anisotropy in seismic velocities and amplitudes from multiple parallel fractures , 1990 .

[43]  E. W. James,et al.  Implications of zeolites and their zonation in the Cajon Pass Deep Drillhole , 1988 .