Weighted-averaging finite-element method for 2D elastic wave equations in the frequency domain SCIE SCOPUS

Cited 24 time in WEB OF SCIENCE Cited 27 time in Scopus
Title
Weighted-averaging finite-element method for 2D elastic wave equations in the frequency domain
Author(s)
Min, DJ; Shin, C; Pratt, RG; Yoo, HS
Alternative Author(s)
민동주; 유해수
Publication Year
2003-04
Abstract
we present a weighted-averaging frequency-domain finite-element method for an accurate and efficient 2D elastic wave modeling technique. Our method introduces three kinds of supplementary element sets in addition to a basic element set that is used in the standard finite-element method. By constructing global stiffness and mass matrices for four kinds of element sets and then averaging them with weighting coefficients, we obtain a new global stiffness and mass matrix. With optimal weighting coefficients determined by a Marquardt-Levenberg method to minimize grid dispersion and grid anisotropy, we can reduce the number of nodal points per shear wavelength from 33.3 (using the standard finite-element method) and 20 (using the eclectic method) to 5, with the errors of group velocities no larger than 1%. By reducing the number of grid points per wavelength, we achieve a 97% and 75% reduction of computer memory required to store the complex impedance matrix for a band-type matrix solver and a nested dissection method, respectively, compared with those of the eclectic method. Our method gives approximate solutions compatible with exact solutions for an infinite homogeneous, a semi-infinite homogeneous (Lamb's problem), and a horizontal two-layer model with fewer grid points than the standard and the eclectic method. A major advantage of the weighted-averaging finite-element method for the elastic wave equation is that it provides solutions very close to correct solutions, for Lamb's problem economically, unlike most of the displacement approaches. In addition, our scheme makes the complex impedance matrix symmetric, which satisfies reciprocity. Seismic forward modeling techniques that satisfy reciprocity are of critical importance in seismic imaging and inversion because we can economically calculate a Jacobian matrix using the reciprocity. Successful simulation of a large-size model shows that our method can be used for the simulation of wave propagation in the geological model needed in the reverse-time migration or seismic inversion.
ISSN
0037-1106
URI
https://sciwatch.kiost.ac.kr/handle/2020.kiost/5524
DOI
10.1785/0120020066
Bibliographic Citation
BULLETIN OF THE SEISMOLOGICAL SOCIETY OF AMERICA, v.93, no.2, pp.904 - 921, 2003
Publisher
SEISMOLOGICAL SOC AMER
Subject
SYNTHETIC SEISMOGRAMS; HETEROGENEOUS MEDIA; DIFFERENCE METHOD; SCALAR; PROPAGATION; EXTRAPOLATOR; COMPUTATION; INVERSION; OPERATORS; GRIDS
Type
Article
Language
English
Document Type
Article
Files in This Item:
There are no files associated with this item.

qrcode

Items in ScienceWatch@KIOST are protected by copyright, with all rights reserved, unless otherwise indicated.

Browse