Moment tensors and micromechanical models

Abstract A numerical modelling approach that simulates cracking and failure in rock and the associated seismicity is presented and a technique is described for quantifying the seismic source mechanisms of the modelled events. The modelling approach represents rock as an assemblage of circular particles bonded together at points of contact. The connecting bonds can break under applied stress forming cracks and fractures in the modelled rock. If numerical damping is set to reproduce realistic levels of attenuation, then energy is released when the bonds break and seismic source information can be obtained as damage occurs. A technique is described by which moment tensors and moment magnitudes can be calculated for these simulated seismic events. The technique basically involves integrating around the source and summing the components of force change at the surrounding particle contacts to obtain the elements of the moment tensor matrix. The moment magnitude is then calculated from the eigenvalues of the moment tensor. The modelling approach is tested by simulating a well-controlled experiment in which a tunnel is excavated in highly stressed granite while microseismicity is recorded. The seismicity produced by the model is compared to the actual recorded seismicity underground. The model reproduces the spatial and temporal distribution of seismicity observed around the tunnel and also the magnitudes of the events. A direct comparison between the actual and simulated moment tensors is not performed due to the two-dimensional nature of the model, however, qualitative comparisons are presented and it is shown that the model produces intuitively realistic source mechanisms. The ability to obtain seismic source information from the models provides a unique means for model validation through comparison with actual recorded seismicity. Once it is established that the model is performing in a realistic manner, it can then be used to examine the micromechanics of cracking, failure and the associated seismicity and to help resolve the non-uniqueness of the geophysical interpretation. This is demonstrated by examining in detail the mechanics of one of the modelled seismic events by observation of the time dependence of the moment tensor and by direct examination of the particle motions at the simulated source.

[1]  P. Cundall,et al.  A discrete numerical model for granular assemblies , 1979 .

[2]  R. Paul Young,et al.  Moment tensor inversion of induced microseisnmic events: Evidence of non-shear failures in the -4 < M < -2 moment magnitude range , 1992 .

[3]  R. P. Young,et al.  Simulating acoustic emissions in bonded-particle models of rock , 2000 .

[4]  C. Martin,et al.  The strength of massive Lac du Bonnet granite around underground openings , 1993 .

[5]  T. Jordan,et al.  Optimal estimation of scalar seismic moment , 1982 .

[6]  R. P. Young,et al.  Acoustic emission and ultrasonic velocity study of excavation-induced microcrack damage at the underground research laboratory , 1993 .

[7]  David D. Pollard,et al.  Distinct element modeling of deformation bands in sandstone , 1995 .

[8]  R. Paul Young,et al.  Evidence for extensile crack initiation in point source time-dependent moment tensor solutions , 1997 .

[9]  E. Wainwright,et al.  Rockbursts and Seismicity in Mines , 1984 .

[10]  R. Madariaga.  Dynamics of an expanding circular fault , 1976, Bulletin of the Seismological Society of America.

[11]  F. Gilbert.  Excitation of the Normal Modes of the Earth by Earthquake Sources , 1971 .

[12]  R. P. Young,et al.  Point-source inversion neglecting a nearby free surface: simulation of the Underground Research Laboratory, Canada , 2001 .

[13]  A. McGarr,et al.  Some comparisons between mining-induced and laboratory earthquakes , 1994 .

[14]  P. A. Cundall,et al.  Modeling notch-formation mechanisms in the URL Mine-by Test Tunnel using bonded assemblies of circular particles , 1998 .

[15]  R. Paul Young,et al.  Lithological Controls on Seismicity in Granitic Rocks , 2000 .

[16]  C. Scholz.  The Mechanics of Earthquakes and Faulting , 1990 .

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

[18]  G. Backus.  Interpreting the seismic glut moments of total degree two or less , 1977 .

[19]  Peter Mora,et al.  Numerical simulation of earthquake faults with gouge: Toward a comprehensive explanation for the heat flow paradox , 1998 .

[20]  B. Gutenberg,et al.  Seismicity of the Earth and associated phenomena , 1950, MAUSAM.

[21]  J. Morgan,et al.  Numerical simulations of granular shear zones using the distinct element method: 1. Shear zone kinematics and the micromechanics of localization , 1999 .

[22]  John A. Hudson,et al.  Defining the scalar moment of a seismic source with a general moment tensor , 1999 .

[23]  Emmanuel M Detournay,et al.  Discrete element modelling of rock cutting , 1999 .

[24]  Erik Eberhardt,et al.  Numerical modelling of three-dimension stress rotation ahead of an advancing tunnel face , 2001 .

[25]  D. Giardini.  Systematic analysis of deep seismicity: 200 centroid-moment tensor solutions for earthquakes between 1977 and 1980 , 1984 .

[26]  D. Potyondy,et al.  Bonded-Particle Simulations of the In-Situ Failure Test At Olkiluoto , 2001 .

[27]  Brian McPherson,et al.  Simulation of sedimentary rock deformation: Lab‐scale model calibration and parameterization , 2002 .

[28]  P. Cundall,et al.  Modelling Rock Using Bonded Assemblies of Circular Particles , 1996 .

[29]  R. Paul Young,et al.  Micromechanical modeling of cracking and failure in brittle rocks , 2000 .