• Login
    View Item 
    •   Home
    • Research
    • Articles
    • View Item
    •   Home
    • Research
    • Articles
    • View Item
    JavaScript is disabled for your browser. Some features of this site may not work without it.

    Browse

    All of KAUSTCommunitiesIssue DateSubmit DateThis CollectionIssue DateSubmit Date

    My Account

    Login

    Quick Links

    Open Access PolicyORCID LibguideTheses and Dissertations LibguideSubmit an Item

    Statistics

    Display statistics

    Nonuniform 3D finite difference elastic wave simulation on staggered grids

    • CSV
    • RefMan
    • EndNote
    • BibTex
    • RefWorks
    Thumbnail
    Name:
    geo2021-0512.1 (1).pdf
    Size:
    3.732Mb
    Format:
    PDF
    Description:
    Accepted Manuscript
    Download
    Type
    Article
    Authors
    Gao, Longfei
    Ghattas, Omar
    Keyes, David E. cc
    KAUST Department
    Applied Mathematics and Computational Science Program
    Extreme Computing Research Center
    Computer, Electrical and Mathematical Science and Engineering (CEMSE) Division
    Office of the President
    KAUST Grant Number
    OSR-2019-CCF-3666.4
    Date
    2022-03-29
    Permanent link to this record
    http://hdl.handle.net/10754/676318
    
    Metadata
    Show full item record
    Abstract
    We 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
    Citation
    Gao, 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
    Sponsors
    The 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.
    Publisher
    Society of Exploration Geophysicists
    Journal
    GEOPHYSICS
    DOI
    10.1190/geo2021-0512.1
    Additional Links
    https://library.seg.org/doi/10.1190/geo2021-0512.1
    ae974a485f413a2113503eed53cd6c53
    10.1190/geo2021-0512.1
    Scopus Count
    Collections
    Articles; Applied Mathematics and Computational Science Program; Extreme Computing Research Center; Computer, Electrical and Mathematical Science and Engineering (CEMSE) Division

    entitlement

     

    Related items

    Showing items related by title, author, creator and subject.

    • Thumbnail

      Growth-induced axial buckling of a slender elastic filament embedded in an isotropic elastic matrix

      O'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.
    • Thumbnail

      Elastic orthorhombic anisotropic parameter inversion: An analysis of parameterization

      Oh, Juwon; Alkhalifah, Tariq Ali (GEOPHYSICS, Society of Exploration Geophysicists, 2016-09-15) [Article]
      The resolution of a multiparameter full-waveform inversion (FWI) is highly influenced by the parameterization used in the inversion algorithm, as well as the data quality and the sensitivity of the data to the elastic parameters because the scattering patterns of the partial derivative wavefields (PDWs) vary with parameterization. For this reason, it is important to identify an optimal parameterization for elastic orthorhombic FWI by analyzing the radiation patterns of the PDWs for many reasonable model parameterizations. We have promoted a parameterization that allows for the separation of the anisotropic properties in the radiation patterns. The central parameter of this parameterization is the horizontal P-wave velocity, with an isotropic scattering potential, influencing the data at all scales and directions. This parameterization decouples the influence of the scattering potential given by the P-wave velocity perturbation fromthe polar changes described by two dimensionless parameter perturbations and from the azimuthal variation given by three additional dimensionless parameters perturbations. In addition, the scattering potentials of the P-wave velocity perturbation are also decoupled from the elastic influences given by one S-wave velocity and two additional dimensionless parameter perturbations. The vertical S-wave velocity is chosen with the best resolution obtained from S-wave reflections and converted waves, little influence on P-waves in conventional surface seismic acquisition. The influence of the density on observed data can be absorbed by one anisotropic parameter that has a similar radiation pattern. The additional seven dimensionless parameters describe the polar and azimuth variations in the P- and S-waves that we may acquire, with some of the parameters having distinct influences on the recorded data on the earth's surface. These characteristics of the new parameterization offer the potential for a multistage inversion from high symmetry anisotropy to lower symmetry ones. © 2016 Society of Exploration Geophysicists.
    • Thumbnail

      Normalized nonzero-lag crosscorrelation elastic full-waveform inversion

      Zhang, Zhendong; Alkhalifah, Tariq Ali; Wu, Zedong; Liu, Yike; He, Bin; Oh, Juwon (GEOPHYSICS, Society of Exploration Geophysicists, 2018-11-23) [Article]
      Full-waveform inversion (FWI) is an attractive technique due to its ability to build high-resolution velocity models. Conventional amplitude-matching FWI approaches remain challenging because the simplified computational physics used does not fully represent all wave phenomena in the earth. Because the earth is attenuating, a sample-by-sample fitting of the amplitude may not be feasible in practice. We have developed a normalized nonzero-lag crosscorrelataion-based elastic FWI algorithm to maximize the similarity of the calculated and observed data. We use the first-order elastic-wave equation to simulate the propagation of seismic waves in the earth. Our proposed objective function emphasizes the matching of the phases of the events in the calculated and observed data, and thus, it is more immune to inaccuracies in the initial model and the difference between the true and modeled physics. The normalization term can compensate the energy loss in the far offsets because of geometric spreading and avoid a bias in estimation toward extreme values in the observed data. We develop a polynomial-type weighting function and evaluate an approach to determine the optimal time lag. We use a synthetic elastic Marmousi model and the BigSky field data set to verify the effectiveness of the proposed method. To suppress the short-wavelength artifacts in the estimated S-wave velocity and noise in the field data, we apply a Laplacian regularization and a total variation constraint on the synthetic and field data examples, respectively.
    DSpace software copyright © 2002-2022  DuraSpace
    Quick Guide | Contact Us | KAUST University Library
    Open Repository is a service hosted by 
    Atmire NV
     

    Export search results

    The export option will allow you to export the current search results of the entered query to a file. Different formats are available for download. To export the items, click on the button corresponding with the preferred download format.

    By default, clicking on the export buttons will result in a download of the allowed maximum amount of items. For anonymous users the allowed maximum amount is 50 search results.

    To select a subset of the search results, click "Selective Export" button and make a selection of the items you want to export. The amount of items that can be exported at once is similarly restricted as the full export.

    After making a selection, click one of the export format buttons. The amount of items that will be exported is indicated in the bubble next to export format.