Nonuniform 3D finite difference elastic wave simulation on staggered grids
KAUST DepartmentApplied Mathematics and Computational Science Program
Extreme Computing Research Center
Computer, Electrical and Mathematical Science and Engineering (CEMSE) Division
Office of the President
KAUST Grant NumberOSR-2019-CCF-3666.4
Permanent link to this recordhttp://hdl.handle.net/10754/676318
MetadataShow full item record
AbstractWe present an approach to simulate the 3D isotropic elastic wave propagation using nonuniform finite difference discretization on staggered grids. Specifically, we consider simulation domains composed of layers of uniform grids with different grid spacings, separated by nonconforming interfaces. We demonstrate that this layer-wise finite difference discretization has the potential to significantly reduce the simulation cost, compared to its fully uniform counterpart, by more than an order of magnitude on geophysically representative 3D problems. Stability of such a discretization is achieved by using specially designed operators, which are variants of the standard finite difference operators with adaptations near boundaries or interfaces, and penalty terms, which are appended to the discretized wave system to weakly impose boundary or interface conditions. Combined with specially designed interpolation operators, the discretized wave system is shown to preserve the energy conserving property of the continuous elastic wave equation, and a fortiori ensure the stability of the simulation. Numerical examples are presented to demonstrate the efficacy of the proposed simulation approach
CitationGao, L., Ghattas, O., & Keyes, D. (2022). Nonuniform 3D finite difference elastic wave simulation on staggered grids. GEOPHYSICS, 1–79. https://doi.org/10.1190/geo2021-0512.1
SponsorsThe authors thank the editors and referees for carefully reviewing this article and providing valuable suggestions. This research used resources of the Core Labs of King Abdullah University of Science and Technology (KAUST). This research was funded by KAUST award OSR-2019-CCF-3666.4 and U.S. National Science Foundation Frontera award 2033468.
PublisherSociety of Exploration Geophysicists
Showing items related by title, author, creator and subject.
Growth-induced axial buckling of a slender elastic filament embedded in an isotropic elastic matrixO'Keeffe, Stephen G.; Moulton, Derek E.; Waters, Sarah L.; Goriely, Alain (International Journal of Non-Linear Mechanics, Elsevier BV, 2013-11) [Article]We investigate the problem of an axially loaded, isotropic, slender cylinder embedded in a soft, isotropic, outer elastic matrix. The cylinder undergoes uniform axial growth, whilst both the cylinder and the surrounding elastic matrix are confined between two rigid plates, so that this growth results in axial compression of the cylinder. We use two different modelling approaches to estimate the critical axial growth (that is, the amount of axial growth the cylinder is able to sustain before it buckles) and buckling wavelength of the cylinder. The first approach treats the filament and surrounding matrix as a single 3-dimensional elastic body undergoing large deformations, whilst the second approach treats the filament as a planar, elastic rod embedded in an infinite elastic foundation. By comparing the results of these two approaches, we obtain an estimate of the foundation modulus parameter, which characterises the strength of the foundation, in terms of the geometric and material properties of the system. © 2013 Elsevier Ltd. All rights reserved.
Target-oriented high-resolution elastic full waveform inversion with an elastic redatuming methodLi, Yuanyuan; Alkhalifah, Tariq Ali (GEOPHYSICS, Society of Exploration Geophysicists, 2022-05-25) [Article]Characterizing the elastic properties in deep-buried reservoirs beneath complex overburden structures remains challenging for seismic inversion. Elastic full waveform inversion (FWI) is capable of quantitatively estimating the subsurface elastic properties with reasonably high resolution. However, elastic FWI using high frequencies is computationally expensive because fine discretization is required to stabilize the wavefield simulation. Additionally, it is challenging for elastic FWI to obtain high-resolution inversion results of the deep targets due to complex overburden structures and limited energy illumination of the target of interest. To overcome these limitations, we propose a target-oriented high-resolution elastic FWI scheme by using estimated elastic data for a virtual survey deployed just above a zone of interest. Specifically, we present an elastic redatuming approach to retrieve virtual elastic data by solving an iterative least-squares data fitting problem. In addition to the elastic multicomponent data recorded at the acquisition surface, an estimate of the overburden model is required. A synthetic data example shows that an overburden model estimated by elastic FWI with a low-frequency band enables our elastic redatuming approach to reconstruct the virtual elastic data representing the seismic reflection response from the target zone. We then perform regularized elastic FWI by using the redatumed elastic data to estimate the elastic properties in the target zone with reduced computational cost. In the numerical examples, the Marmousi2 synthetic data and 2D Volve field data are used to demonstrate the feasibility and practicality of the proposed method.
geodynamics/specfem3d: SPECFEM3D_Cartesian simulates acoustic (fluid), elastic (solid), coupled acoustic/elastic, poroelastic or seismic wave propagation in any type of conforming mesh of hexahedra (structured or not). It can for instance model seismic waves propagating in sedimentary basins or any other regional geological model following earthquakes. It can also be used for non destructive testing or for ocean acoustics.Abreu, Rafael; Peter, Daniel; Thomas, Christine (Github, 2013-11-11) [Software]SPECFEM3D_Cartesian simulates acoustic (fluid), elastic (solid), coupled acoustic/elastic, poroelastic or seismic wave propagation in any type of conforming mesh of hexahedra (structured or not). It can for instance model seismic waves propagating in sedimentary basins or any other regional geological model following earthquakes. It can also be used for non destructive testing or for ocean acoustics.