A 3D discrete numerical elastic lattice method for seismic wave propagation in heterogeneous media with topography

[1] A three-dimensional elastic lattice method for the simulation of seismic waves is presented. The model consists of particles arranged on a cubic lattice which interact through a central force term and a bond-bending force. Particle disturbances are followed through space by numerically solving their equations of motion. A vacuum free-surface boundary condition is implicit in the method. We demonstrate that a numerical implementation of the method is capable of modelling seismic wave propagation with complex topography. This is achieved by comparing the scheme against a finite-difference solution to the elastodynamic wave equation. The results indicate that the scheme offers an alternative 3D method for modelling wave propagation in the presence of strong topography and subsurface heterogeneity. We apply the method to seismic wave propagation on Mount Etna to illustrate its applicability in modelling a physical system.

[1]  F. Sánchez-Sesma,et al.  Rayleigh waves modeling using an elastic lattice model , 2003 .

[2]  Arbabi,et al.  Elastic properties of three-dimensional percolation networks with stretching and bond-bending forces. , 1988, Physical review. B, Condensed matter.

[3]  H. Igel,et al.  Seismic wave simulation in the presence of real volcano topography , 2003 .

[4]  M. P. Anderson,et al.  Elastic and fracture properties of the two-dimensional triangular and square lattices , 1994 .

[5]  E. M. Lifshitz,et al.  Theory of Elasticity: Vol. 7 of Course of Theoretical Physics , 1960 .

[6]  Christopher J. Bean,et al.  Numerical simulation of seismic waves using a discrete particle scheme , 2000 .

[7]  Christopher J. Bean,et al.  Fracture properties from seismic data - a numerical investigation , 2002 .

[8]  Bernard A. Chouet,et al.  A free-surface boundary condition for including 3D topography in the finite-difference method , 1997, Bulletin of the Seismological Society of America.

[9]  S. Timoshenko,et al.  Theory of elasticity , 1975 .

[10]  R. Olness,et al.  Two‐dimensional computer studies of crystal stability and fluid viscosity , 1974 .

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

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

[13]  David D. Pollard,et al.  Distinct element modeling of structures formed in sedimentary overburden by extensional reactivation of basement normal faults , 1992 .