Nonsplit complex-frequency shifted perfectly matched layer combined with symplectic methods for solving second-order seismic wave equations — Part 1: Method

Geophysics ◽  
2018 ◽  
Vol 83 (6) ◽  
pp. T301-T311 ◽  
Author(s):  
Xiao Ma ◽  
Dinghui Yang ◽  
Xueyuan Huang ◽  
Yanjie Zhou
Keyword(s):  
Boundary Condition ◽  
Seismic Wave ◽  
Wave Equations ◽  
Perfectly Matched Layer ◽  
Absorbing Boundary Conditions ◽  
Second Order ◽  
Complex Frequency ◽  
Absorbing Boundary ◽  
Time Layer ◽  
The Time Domain

The absorbing boundary condition plays an important role in seismic wave modeling. The perfectly matched layer (PML) boundary condition has been established as one of the most effective and prevalent absorbing boundary conditions. Among the existing PML-type conditions, the complex frequency shift (CFS) PML attracts considerable attention because it can handle the evanescent and grazing waves better. For solving the resultant CFS-PML equation in the time domain, one effective technique is to apply convolution operations, which forms the so-called convolutional PML (CPML). We have developed the corresponding CPML conditions with nonconstant grid compression parameter, and used its combination algorithms specifically with the symplectic partitioned Runge-Kutta and the nearly analytic SPRK methods for solving second-order seismic wave equations. This involves evaluating second-order spatial derivatives with respect to the complex stretching coordinates at the noninteger time layer. Meanwhile, two kinds of simplification algorithms are proposed to compute the composite convolutions terms contained therein.

Nonsplit complex-frequency-shifted perfectly matched layer combined with symplectic methods for solving second-order seismic wave equations — Part 2: Wavefield simulations

Geophysics ◽  
2019 ◽  
Vol 84 (3) ◽  
pp. T167-T179 ◽  
Author(s):  
Xiao Ma ◽  
Dinghui Yang ◽  
Xijun He ◽  
Xueyuan Huang ◽  
Jiaxing Song
Keyword(s):  
Boundary Condition ◽  
Wave Equation ◽  
Free Surface ◽  
Seismic Wave ◽  
Wave Equations ◽  
Perfectly Matched Layer ◽  
Second Order ◽  
Surface Boundary ◽  
Complex Frequency ◽  
Symplectic Methods

The perfectly matched layer (PML) is an efficient artificial boundary condition that has been routinely implemented in seismic wave modeling. However, the effective combination of PML with symplectic numerical schemes for solving seismic wave equations has rarely been studied. In a companion paper, we have developed a complex-frequency-shifted convolutional PML (CPML) with a nonconstant compression grid parameter for solving the time-domain second-order seismic wave equation. Subsequently, we combine this CPML with two classes of symplectic methods to formulate symplectic partitioned Runge-Kutta (SPRK) + CPML and nearly analytic SPRK (NSPRK) + CPML, both of which are properly synchronized. To further investigate their validity, the two algorithms are then applied to acoustic and elastic wave simulations in typical geologic models, including a heterogeneous acoustic model, several isotropic and orthotropic elastic models, and an isotropic elastic model with a free-surface boundary. Relevant numerical results demonstrate the effectiveness of our CPML and combination algorithms. Specifically, the numerical accuracy and stability of the CPML that we develop are greatly improved compared with the classic split-field PML. Moreover, the final model with the free-surface boundary condition indicates that the nonconstant grid-compression parameter can eliminate the unstable modes at the free surface in the PML domain. The (N)SPRK + CPML that we propose is prospective for future application in other complex models and wave-equation-based migration and inversion.

A stable auxiliary differential equation perfectly matched layer condition combined with low-dispersive symplectic methods for solving second-order elastic wave equations

Geophysics ◽  
2019 ◽  
Vol 84 (4) ◽  
pp. T193-T206 ◽  
Author(s):  
Xiao Ma ◽  
Yangjia Li ◽  
Jiaxing Song
Keyword(s):  
Differential Equation ◽  
Elastic Wave ◽  
Wave Equations ◽  
Perfectly Matched Layer ◽  
Second Order ◽  
Elastic Model ◽  
Complex Frequency ◽  
Symplectic Method ◽  
Symplectic Methods ◽  
Compression Parameter

The stable implementation of the perfectly matched layer (PML), one of the most effective and popular artificial boundary conditions, has attracted much attention these years. As a type of low-dispersive and symplectic method for solving seismic wave equations, the nearly-analytic symplectic partitioned Runge-Kutta (NSPRK) method has been combined with split-field PML (SPML) and convolutional complex-frequency shifted PML (C-CFS-PML) previously to model acoustic and short-time elastic wave modelings, not yet successfully applied to long-time elastic wave propagation. In order to broaden the application of NSPRK and more general symplectic methods for second-order seismic models, we formulate an auxiliary differential equation (ADE)-CFS-PML with a stabilizing grid compression parameter. This includes deriving the ADE-CFS-PML equations and formulating an adequate time integrator to properly embed their numerical discretizations in the main symplectic numerical methods. The resulting (N)SPRK+ADE-CFS-PML algorithm can help break through the constraint of at most second-order temporal accuracy that used to be imposed on SPML and C-CFS-PML. Especially for NSPRK, we implement the strategy of neglecting the treatment of third-order spatial derivatives in the PML domain and obtain an efficient absorption effect. Related acoustic and elastic wave simulations illustrate the enhanced numerical accuracy of our ADE-CFS-PML compared with SPML and C-CFS-PML. The elastic wave simulation in a homogeneous isotropic medium shows that compared to NSPRK+C-CFS-PML, the NSPRK+ADE-CFS-PML is numerically stable throughout a simulation time of 2 s. The synthetic seismograms of the 2D acoustic SEG salt model and the two-layer elastic model demonstrate the effectiveness of NSPRK+ADE-CFS-PML for complex elastic models. The stabilization effect of the grid compression parameter is verified in the final homogeneous isotropic elastic model with free-surface boundary.

scholarly journals ABSORBING BOUNDARY CONDITIONS FOR A WAVE EQUATION WITH A TEMPERATURE-DEPENDENT SPEED OF SOUND

2013 ◽  
Vol 21 (02) ◽  
pp. 1250028 ◽  
Author(s):  
IGOR SHEVCHENKO ◽  
MANFRED KALTENBACHER ◽  
BARBARA WOHLMUTH
Keyword(s):  
Wave Equation ◽  
Boundary Conditions ◽  
Speed Of Sound ◽  
Wave Equations ◽  
Focused Ultrasound ◽  
Absorbing Boundary Conditions ◽  
Second Order ◽  
High Intensity Focused Ultrasound ◽  
Temperature Dependent ◽  
Absorbing Boundary

In this work, new absorbing boundary conditions (ABCs) for a wave equation with a temperature-dependent speed of sound are proposed. Based on the theory of pseudo-differential calculus, first- and second-order ABCs for the one- and two-dimensional wave equations are derived. Both boundary conditions are local in space and time. The well-posedness of the wave equation with the developed ABCs is shown through the reduction of the original problem to an equivalent one for which the uniqueness and existence of the solution has already been established. Although the second-order ABC is more accurate, the numerical realization is more challenging. Here we use a Lagrange multiplier approach which fits into the abstract framework of saddle point formulations and yields stable results. Numerical examples illustrating stability, accuracy and flexibility of the ABCs are given. As a test setting, we perform computations for a high-intensity focused ultrasound (HIFU) application, which is a typical thermo-acoustic multi-physics problem.

scholarly journals A Well-Posed and Discretely Stable Perfectly Matched Layer for Elastic Wave Equations in Second Order Formulation

2012 ◽  
Vol 11 (5) ◽  
pp. 1643-1672 ◽  
Author(s):  
Kenneth Duru ◽  
Gunilla Kreiss
Keyword(s):  
Elastic Wave ◽  
Wave Equations ◽  
Perfectly Matched Layer ◽  
Second Order ◽  
Order System ◽  
Stable Node ◽  
Complex Frequency ◽  
Elastic Wave Equations ◽  
Anisotropic Elastic ◽  
Well Posed

AbstractWe present a well-posed and discretely stable perfectly matched layer for the anisotropic (and isotropic) elastic wave equations without first re-writing the governing equations as a first order system. The new model is derived by the complex coordinate stretching technique. Using standard perturbation methods we show that complex frequency shift together with a chosen real scaling factor ensures the decay of eigen-modes for all relevant frequencies. To buttress the stability properties and the robustness of the proposed model, numerical experiments are presented for anisotropic elastic wave equations. The model is approximated with a stable node-centered finite difference scheme that is second order accurate both in time and space.

scholarly journals Frequency-domain elastic wavefield simulation with hybrid absorbing boundary conditions

2019 ◽  
Vol 16 (4) ◽  
pp. 690-706
Author(s):  
Zhencong Zhao ◽  
Jingyi Chen ◽  
Xiaobo Liu ◽  
Baorui Chen
Keyword(s):  
Boundary Condition ◽  
Free Surface ◽  
Boundary Conditions ◽  
Frequency Domain ◽  
Time Domain ◽  
Absorbing Boundary Conditions ◽  
Elastic Media ◽  
Absorbing Boundary ◽  
Wavefield Simulation ◽  
The Time Domain

Abstract The frequency-domain seismic modeling has advantages over the time-domain modeling, including the efficient implementation of multiple sources and straightforward extension for adding attenuation factors. One of the most persistent challenges in the frequency domain as well as in the time domain is how to effectively suppress the unwanted seismic reflections from the truncated boundaries of the model. Here, we propose a 2D frequency-domain finite-difference wavefield simulation in elastic media with hybrid absorbing boundary conditions, which combine the perfectly matched layer (PML) boundary condition with the Clayton absorbing boundary conditions (first and second orders). The PML boundary condition is implemented in the damping zones of the model, while the Clayton absorbing boundary conditions are applied to the outer boundaries of the damping zones. To improve the absorbing performance of the hybrid absorbing boundary conditions in the frequency domain, we apply the complex coordinate stretching method to the spatial partial derivatives in the Clayton absorbing boundary conditions. To testify the validity of our proposed algorithm, we compare the calculated seismograms with an analytical solution. Numerical tests show the hybrid absorbing boundary condition (PML plus the stretched second-order Clayton absorbing condition) has the best absorbing performance over the other absorbing boundary conditions. In the model tests, we also successfully apply the complex coordinate stretching method to the free surface boundary condition when simulating seismic wave propagation in elastic media with a free surface.

A nonsplit complex frequency-shifted PML based on recursive integration for FDTD modeling of elastic waves

Geophysics ◽  
2007 ◽  
Vol 72 (2) ◽  
pp. T9-T17 ◽  
Author(s):  
Francis H. Drossaert ◽  
Antonios Giannopoulos
Keyword(s):  
Boundary Condition ◽  
Elastic Waves ◽  
Absorbing Boundary Conditions ◽  
Attractive Alternative ◽  
Complex Frequency ◽  
Electromagnetic Modeling ◽  
Absorbing Boundary ◽  
Split Field ◽  
Recursive Integration ◽  
Stretching Function

In finite-difference time-domain (FDTD) modeling of elastic waves, absorbing boundary conditions are used to mitigate undesired reflections that can arise at the model’s truncation boundaries. The perfectly matched layer (PML) is generally considered to be the best available absorbing boundary condition. An important but rarely addressed limitation of current PML implementations is that their performance is severely reduced when waves are incident on the PML interface at near-grazing angles. In addition, very low frequency waves as well as evanescent waves could cause spurious reflections at the PML interface. In electromagnetic modeling, similar problems are circumvented by using a complex frequency-shifted stretching function in the PML formulas. However, in elastic-wave modeling using the conventional PML formulation — based on splitting the velocityand stressfields — it is difficult to adopt a complex frequency--shifted stretching function. We present an alternative implemen-tation of a PML that is based on recursive integration and does not require splitting of the velocity and stress fields. Modeling re-sults show that the performance of our implementation using a standard stretching function is identical to that of the convention-al split-field PML. Then we show that the new PML can be modi-fied easily to include the complex frequency-shifted stretching function. Results of models with an elongated domain show that this modification can substantially improve the performance of the PML boundary condition. An efficient implementation of the new PML requires less memory than the conventional split-field PML, and, therefore, is a very attractive alternative to the con-ventional PML. By adopting the complex frequency-shifted stretching function, the PML can accommodate a wide variety of model problems, and hence it is more generic.

Instabilities in applying absorbing boundary conditions to high‐order seismic modeling algorithms

Geophysics ◽  
1998 ◽  
Vol 63 (3) ◽  
pp. 1017-1023 ◽  
Author(s):  
Antonio Simone ◽  
Stig Hestholm
Keyword(s):  
Boundary Condition ◽  
Acoustic Wave ◽  
Elastic Wave ◽  
Wave Equations ◽  
Absorbing Boundary Conditions ◽  
Absorbing Boundary Condition ◽  
Seismic Modeling ◽  
Absorbing Boundary ◽  
Numerical Grid ◽  
Numerical Discretization

The problem of artificial reflections from grid boundaries in the numerical discretization of elastic and acoustic wave equations has long plagued geophysicists. Even if modern computers have made it possible to extend the synthetics over more wavelengths (equivalent to larger propagation distances), efficient absorption methods are still needed to minimize interference from unwanted reflections from the numerical grid boundaries. In this study, we examine applicabilities and stabilities of the optimal absorbing boundary condition (OABC) of Peng and Toksöz (1994, 1995) for 2-D and 3-D acoustic and elastic wave modeling. As a basis for comparison, we use exponential damping (ED) (Cerjan et. al., 1985), in which velocities and stresses are multiplied by progressively decreasing terms when approaching the boundaries of the numerical grid.

Application of perfectly matched layer for scalar arbitrarily wide-angle wave equations

Geophysics ◽  
2013 ◽  
Vol 78 (1) ◽  
pp. T29-T39 ◽  
Author(s):  
Hanming Chen ◽  
Hui Zhou ◽  
He Lin ◽  
Shangxu Wang
Keyword(s):  
Wave Equation ◽  
Wave Equations ◽  
Perfectly Matched Layer ◽  
Evanescent Waves ◽  
Absorbing Boundary Conditions ◽  
Wide Angle ◽  
Absorbing Boundary ◽  
Computation Cost ◽  
Propagating Waves ◽  
Derivation Procedure

Arbitrarily wide-angle wave equation (AWWE) is a space domain, high-order one-way wave equation (OWWE). Its accuracy can be arbitrarily increased, and it is amenable to easy numerical implementation. Those properties make it outstanding among the existing OWWEs and further enable it to be a desirable tool for migration. We extend the perfectly matched layer (PML) to 3D scalar AWWE to provide a good approach to suppress artifacts arising at truncation boundaries. We follow the concept of complex coordinate stretching, and the derivation procedure of PML for AWWE is straightforward. An existing finite-difference scheme is adopted to fit the split PML formulation and its stability is observed through numerical examples. The performance of the developed PML condition is compared with two different wave-equation based absorbing boundary conditions. Numerical results illustrate that the PML condition used in AWWE propagator can effectively absorb the propagating waves and evanescent waves at a price of limited additional computation cost.

Sign in / Sign up

Export Citation Format

Share Document