Application of the finite difference method in seismic source and wave diffraction simulation

Hong, Moonki; Bond, Leonard J.

doi:10.1111/j.1365-246x.1986.tb01969.x

Cited by 11 publications

(8 citation statements)

References 11 publications

(6 reference statements)

Supporting

Mentioning

Contrasting

Order By: Relevance

“…The decomposition of the elastic wavefield into separate P-and S-wave potential fields opens up a convenient way of representing seismic sources numerically, in particular explosive sources. Hong & Bond (1986) used an area source to simulate an explosive source. With the displacement-potential formulation, a source with only P-wave potential simulates an extensional/compressional explosive source whereas a source with only S-wave potential simulates a unidirectional shear source.…”

Section: Source Excitationmentioning

confidence: 99%

Finite-difference simulation of P--SV-wave propagation: a displacement-potential approach

Cao

Greenhalgh

1992

Geophysical Journal International

View full text Add to dashboard Cite

S U M M A R YA new finite-difference formulation is presented for the simulation of P-SV-wave propagation in heterogeneous, isotropic media. In this approach, the wave equations are separated as two sets of equations: one for the displacement fields and the other for potential fields. These two sets of equations involve only first-order spatial derivatives. In terms of potentials, t h e P-SV-wavefield can be split into P-and S-wave potential fields, which may open an opportunity of simulating P-and S-wave propagation in different ways. By assuming constant density and shear modulus, both P-and S-wave potentials can be expressed in the form of scalar wave equations. Thus, the Lindman 'free space' boundary condition (Lindman 1975) for scalar waves can be used. An improvement can be made by adding a dissipation zone to absorbing boundaries. F o r a 2-D model, four equations are used, which is one less than Virieux's (1986) velocity-stress approach. Therefore, the new algorithm is more efficient and requires less computer memory. T h e calculation is naturally performed in a staggered-grid manner which gives an excellent result for both internal discontinuities and model edges.

show abstract

Section: Source Excitationmentioning

confidence: 99%

Finite-difference simulation of P--SV-wave propagation: a displacement-potential approach

Cao

Greenhalgh

1992

Geophysical Journal International

View full text Add to dashboard Cite

show abstract

“…So far most applications of finite differences have been limited to simple geometric shapes. For instance, much earlier work deals with restrictive models in which the boundaries are either parallel to the coordinates or permit only a 45 degree ramp (Alterman and Karal, 1968;Loewenthal, 1970, 1972;Munasinghe and Farnell, 1973;Ilan et al, 1975Ilan et al, , 1979Ilan, 1978;Fuyuki and Matsumoto, 1980;Fuyuki and Nakano, 1984;Hong and Bond, 1986). Boore et al (1981) simulated vertically incident body waves impinging on a 45 degree ramp by using the same formulation as in Ilan et al (1975).…”

Section: Introductionmentioning

confidence: 97%

Free‐boundary conditions of arbitrary polygonal topography in a two‐dimensional explicit elastic finite‐difference scheme

1988

View full text Add to dashboard Cite

We present a simple method for simulating 2-D elastic waves in a model with free-surface topography of polygonal shape, i.e., a continuous but irregular surface composed of line segments. Our method requires special treatment for each of the six specific cases involving line segments of various slopes as well as transition points between the sloping segments. For brevity, only nonnegatively sloping segments are specifically included.On an inclined free surface, vanishing stress conditions are implemented using a rotated coordinate system parallel to the inclined boundary. At transition points on the topography between line segments, we use a first-order approximate boundary condition in a locally rotated coordinate system aligned with the bisector of the corner. As in the classical one-sided explicit approximation scheme widely used for the flat freesurface case, these extrapolation formulas are accurate to first order in spatial increment. Numerical tests indicate that the present scheme is stable over a range of Poisson's ratios of practical interest (v > 0.3) for fairly complicated geometric shapes consisting of ridges and valleys with both steep and gentle slopes. Stability for complicated shapes enables us to study realistic problems for which the topography plays a significant role in shaping the wave field and for which analytical solutions are not generally available.

show abstract

“…Source functions for both point and finite elastic dislocations can be modelled through specification of equivalent body forces (cf. Burridge, Lapwood & Knopoff 1964) implemented in finite-difference form (Alterman & Aboudi 1970;Aboudi 1971;Hong & Bond 1986). A propagating rupture is easily simulated by sequentially initiating a series of point sources oriented parallel to the fault, with displacements parallel to the slip direction.…”

Section: Theorymentioning

confidence: 99%

“…The latter include full waveform modelling (cf. Langston & Helmberger 1975;Archuleta 1976;Hong & Bond 1986;Benz & Smith 1987) and inversion (cf. Stump & Johnson 1977;Spudich 1980).…”

Section: Introductionmentioning

confidence: 99%

Wavefield extrapolation of body waves for 3-D imaging of earthquake sources

Chang

McMechan

1991

Geophysical Journal International

View full text Add to dashboard Cite

S U M M A R Y Reverse-time imaging of earthquake source parameters is extended from two to three spatial dimensions and from two-to three-component recordings. Provided that the recording aperture is sufficiently large, and the data are not spatially aliased, source radiation patterns, in both time and space, can be reconstructed by elastic reverse-time propagation of body waves; such reconstructions are necessarily partial as only the energy that was recorded is available for reconstruction. For a point source, the origin time and the (3-D) spatial location of an event can be reconstructed by extracting the time and position, of the best focused energy from the backward propagating wavefield. For a spatially and temporally extended source, biased estimates of fault position and rupture/slip time history can be estimated if the recording aperture is sufficient; for smaller recording apertures, the time history can still be estimated if the fault location and geometry are known a priori. The latter is viable for a reasonable number (<200) of three-component recordings, and so is potentially applicable to real data. All these processes assume that a smoothed representation of the 3-D velocity distribution in the volume containing the source and receivers is available.

show abstract

Application of the finite difference method in seismic source and wave diffraction simulation

Cited by 11 publications

References 11 publications

Finite-difference simulation of P--SV-wave propagation: a displacement-potential approach

Finite-difference simulation of P--SV-wave propagation: a displacement-potential approach

Free‐boundary conditions of arbitrary polygonal topography in a two‐dimensional explicit elastic finite‐difference scheme

Wavefield extrapolation of body waves for 3-D imaging of earthquake sources

Contact Info

Product

Resources

About