Modelling the compliance of crustal rock—II. Response to temporal changes before earthquakes

SUMMARY There have been several claims that seismic shear waves respond to changes in stress before earthquakes. The companion paper develops a stress-sensitive model (APE) for the behaviour of low-porosity low-permeability crystalline rocks containing pervasive distributions of fluid-filled intergranular microcracks, and this paper uses APE to model the behaviour before earthquakes. Modelling with APE shows that the microgeometry and statistics of distributions of such fluid-filled microcracks respond almost immediately to changes in stress, and that the behaviour can be monitored by analysing seismic shear-wave splitting. The physical reasons for the coupling between shear-wave splitting and differential stress are discussed. In this paper, we extend the model by using percolation theory to show that large crack densities are limited at the grain-scale level by the percolation threshold at which interacting crack clusters lead to pronounced increases in rock-matrix permeability. In the simplest formulation, the modelling is dimensionless and almost entirely constrained without free parameters. Nevertheless, APE modelling of the evolution of fluid-saturated rocks under stress reproduces the observed fracture criticality and the narrow range of shear-wave azimuthal anisotropy in crustal rocks. It also reproduces the behaviour of temporal variations in shear-wave splitting observed before and after the 1986, M= 6, North Palm Springs earthquake, Southern California, and several other smaller earthquakes. The agreement of APE modelling with a wide range of observations confirms that fluid-saturated crystalline rocks are stress-sensitive and respond to changes in stress by critical fluid-rock interactions at the microscale level. This means that the effects of changes in stress and other parameters can be numerically modelled and monitored by appropriate observations of seismic shear waves.

[1]  M. Biot.  General Theory of Three‐Dimensional Consolidation , 1941 .

[2]  M. Biot.  Theory of Propagation of Elastic Waves in a Fluid-Saturated Porous Solid. II. Higher Frequency Range , 1956 .

[3]  M. Biot.  Theory of Propagation of Elastic Waves in a Fluid‐Saturated Porous Solid. I. Low‐Frequency Range , 1956 .

[4]  Gene Simmons,et al.  Stress‐induced velocity anisotropy in rock: An experimental study , 1969 .

[5]  Amos Nur,et al.  Effects of stress on velocity anisotropy in rocks with cracks , 1971 .

[6]  Amos Nur,et al.  Dilatancy, pore fluids, and premonitory variations of ts/tp travel times , 1972, Bulletin of the Seismological Society of America.

[7]  D. L. Anderson,et al.  Earthquake Prediction: Variation of Seismic Velocities before the San Francisco Earthquake , 1973, Science.

[8]  M. L. Sbar,et al.  Premonitory Changes in Seismic Velocities and Prediction of Earthquakes , 1973, Nature.

[9]  I. Gupta.  Premonitory Variations in S-Wave Velocity Anisotropy before Earthquakes in Nevada , 1973, Science.

[10]  S. Kirkpatrick.  Percolation and Conduction , 1973 .

[11]  P. Richards,et al.  Spatial and temporal variations in ts /tp and in P wave residuals at Blue Mountain Lake, New York: Application to earthquake prediction , 1975 .

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

[13]  S. Hartzell,et al.  The horse canyon earthquake of August 2, 1975—Two-stage stress-release process in a strike-slip earthquake , 1979, Bulletin of the Seismological Society of America.

[14]  S. Crampin,et al.  Observations of dilatancy-induced polarization anomalies and earthquake prediction , 1980, Nature.

[15]  B. Wilkins.  Slow crack growth and delayed failure of granite , 1980 .

[16]  Comments on papers about shear-wave splitting in dilatancy-induced anisotropy by I. N. Gupta and by A. Ryall and W. U. Savage , 1981 .

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

[18]  Stuart Crampin,et al.  A review of wave motion in anisotropic and cracked elastic-media , 1981 .

[19]  A. Provost,et al.  Scaling rules in rock fracture and possible implications for earthquake prediction , 1982, Nature.

[20]  Tamaz Chelidze,et al.  Percolation and fracture , 1982 .

[21]  J. Dienes.  PERMEABILITY, PERCOLATION AND STATISTICAL CRACK MECHANICS , 1982 .

[22]  R. Kranz.  Microcracks in rocks: a review , 1983 .

[23]  S. Crampin.  Anisotropy in exploration seismics , 1984 .

[24]  J. B. Fletcher,et al.  Analysis of analog and digital records of the 1982 Arkansas earthquake swarm , 1984 .

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

[26]  B. Atkinson.  Subcritical crack growth in geological materials , 1984 .

[27]  Stuart Crampin.  Effective anisotropic elastic constants for wave propagation through cracked solids , 1984 .

[28]  Stuart Crampin,et al.  Analysis of records of local earthquakes: the Turkish Dilatancy Projects (TDP1 and TDP2) , 1985 .

[29]  E. Rutter,et al.  On the relationship between deformation and metamorphism, with special reference to the behavior of basic rocks , 1985 .

[30]  David C. Booth,et al.  Shear-wave polarizations on a curved wavefront at an isotropic free surface , 1985 .

[31]  Iain Bush,et al.  Estimating the internal structure of reservoirs with shear‐wave VSPs , 1986 .

[32]  R. M. Alford,et al.  Shear data in the presence of azimuthal anisotropy: Dilley, Texas , 1986 .

[33]  T. Chelidze.  Percolation theory as a tool for imitation of fracture process in rocks , 1986 .

[34]  H. A. Willis,et al.  Azimuthal Anisotropy: Occurrence And Effect On Shear-Wave Data Quality , 1986 .

[35]  S. Crampin.  Geological and industrial implications of extensive-dilatancy anisotropy , 1987, Nature.

[36]  David Beamish,et al.  The Turkish Dilatancy Project (TDP3): multidisciplinary studies of a potential earthquake source region , 1987 .

[37]  C. Rai,et al.  Shear‐wave velocity anisotropy in sedimentary rocks: A laboratory study , 1988 .

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

[39]  J. Sochacki.  Absorbing boundary conditions for the elastic wave equations , 1988 .

[40]  Colin M. Sayers,et al.  Stress-induced ultrasonic wave velocity anisotropy in fractured rock , 1988 .

[41]  S. Crampin,et al.  Shear-wave splitting showing hydraulic dilation of pre-existing joints in granite , 1989 .

[42]  S. Crampin.  SUGGESTIONS FOR A CONSISTENT TERMINOLOGY FOR SEISMIC ANISOTROPY , 1989 .

[43]  J. Dienes,et al.  Transport properties of rocks from statistics and percolation , 1989 .

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

[45]  David Vere-Jones,et al.  Percolation Theory: A Model For Rock Fracture? , 1990 .

[46]  K. Heffer,et al.  Scaling Relationships in Natural Fractures: Data, Theory, and Application , 1990 .

[47]  Peter M. Shearer,et al.  Quantitative measurements of shear wave polarizations at the Anza Seismic Network, southern California: Implications for shear wave splitting and earthquake prediction , 1990 .

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

[49]  S. Crampin.  An alternative scenario for earthquake prediction experiments , 1991 .

[50]  J. B. Fletcher,et al.  Comment on “Quantitative measurements of shear wave polarizations at the Anza Seismic Network, southern California: Implications for shear wave splitting and earthquake prediction” by Richard C. Aster, Peter M. Shearer, and Jon Berger , 1991 .

[51]  P. Shearer,et al.  Reply [to “Comment on ‘Quantitative measurements of shear wave polarizations at the Anza Seismic Network, southern California: Implications for shear wave splitting and earthquake prediction’ by Richard C. Aster, Peter M. Shearer, and Jon Berger”] , 1991 .

[52]  J. Munster,et al.  Microcrack-induced seismic anisotropy of sedimentary rocks , 1991 .

[53]  S. Siegesmund,et al.  The effect of oriented microcracks on seismic velocities in an ultramylonite , 1991 .

[54]  D. Turcotte.  Fractals, chaos, self-organized criticality and tectonics , 1992 .

[55]  I. W. Farmer,et al.  Fluid Flow in Discontinuous Rocks , 1993 .

[56]  Amos Nur,et al.  Dynamic poroelasticity: A unified model with the squirt and the Biot mechanisms , 1993 .

[57]  S. Crampin.  Do you know of an isolated swarm of small earthquakes , 1993 .

[58]  Mark A. Meadows,et al.  Seismic detection of a hydraulic fracture from shear-wave VSP data at Lost Hills Field, California , 1994 .

[59]  Stuart Crampin,et al.  The fracture criticality of crustal rocks , 1994 .

[60]  Y. Géraud.  Variations of connected porosity and inferred permeability in a thermally cracked granite , 1994 .

[61]  S. I. Tkachenko,et al.  Phase equilibria in fluid systems at high pressures and temperatures , 1994 .

[62]  C. Neuzil.  Abnormal pressures as hydrodynamic phenomena , 1995 .

[63]  K. Shmulovich,et al.  Melting of albite and dehydration of brucite in H2O–NaCl fluids to 9 kbars and 700–900°C: implications for partial melting and water activities during high pressure metamorphism , 1996 .

[64]  J. Hudson,et al.  The mechanical properties of materials with interconnected cracks and pores , 1996 .

[65]  J. Hudson,et al.  Diffraction of seismic waves by cracks with application to hydraulic fracturing , 1997 .

[66]  Stuart Crampin,et al.  Modelling the compliance of crustal rock—I. Response of shear‐wave splitting to differential stress , 1997 .

[67]  Yuan Gao,et al.  Temporal changes in shear-wave splitting at an isolated swarm of small earthquakes in 1992 near Dongfang, Hainan Island, southern China , 1998 .