Efficient Global Matrix Approach to the Computation of Synthetic Seismograms,
SACLANT ASW RESEARCH CENTRE LA SPEZIA (ITALY)
Pagination or Media Count:
A numerically efficient global matrix approach to the solution of the wave equation in horizontally stratified environments is presented. The field in each layer is expressed as a superposition of the field produced by the sources within the layer and an unknown field satisfying the homogeneous wave equations, both expressed as integral representations in the horizontal wavenumber. The boundary conditions to be satisfied at each interface then yield a linear system of equations in the unknown wavefield amplitudes, to be satisfied at each horizontal wavenumber. As an alternative to the traditional propagator matrix approaches, this solution technique yields both improved efficiency and versatility. Its global nature makes it well suited to problems involving many receivers in range as well as depth and to calculations of both stresses and particle velocities. The global solution technique is developed in close analogy to the finite element method, thereby reducing the number of arithmetic operation to a minimum and making the resulting computer code very efficient in terms of computation time.
- Theoretical Mathematics