Application of the multiaxial perfectly matched layer (M-PML) to near-surface seismic modeling with Rayleigh waves

Perfectly matched layer (PML) absorbing boundaries are widely used to suppress spurious edge reflections in seismic modeling. When modeling Rayleigh waves with the existence of the free surface, the classical PML algorithm becomes unstable when the Poisson’s ratio of the medium is high. Numerical errors can accumulate exponentially and terminate the simulation due to computational overflows. Numerical tests show that the divergence speed of the classical PML has a nonlinear relationship with the Poisson’s ratio. Generally, the higher the Poisson’s ratio, the faster the classical PML diverges. The multiaxial PML (M-PML) attenuates the waves in PMLs using different damping profiles that are proportional to each other in orthogonal directions. The proportion coefficients of the damping profiles usually vary with the specific model settings. If they are set appropriately, the M-PML algorithm is stable for high Poisson’s ratio earth models. Through numerical tests of 40 models with Poisson’s ratios that varied from 0.10 to 0.49, we found that a constant proportion coefficient of 1.0 for the x- and z-directional damping profiles is sufficient to stabilize the M-PML for all 2D isotropic elastic cases. Wavefield simulations indicate that the instability of the classical PML is strongly related to the wave phenomena near the free surface. When applying the multiaxial technique only in the corners of the PML near the free surface, the original M-PML technique can be simplified without losing its stability. The simplified M-PML works efficiently for homogeneous and heterogeneous earth models with high Poisson’s ratios. The analysis in this paper is based on 2D finite difference modeling in the time domain that can easily be extended into the 3D domain with other numerical methods.

[1]  Yinhe Luo,et al.  Rayleigh-Wave Dispersive Energy Imaging Using a High-Resolution Linear Radon Transform , 2008 .

[2]  Jianghai Xia,et al.  Inversion of high frequency surface waves with fundamental and higher modes , 2003 .

[3]  Robert W. Graves,et al.  Simulating seismic wave propagation in 3D elastic media using staggered-grid finite differences , 1996, Bulletin of the Seismological Society of America.

[4]  Moshe Reshef,et al.  A nonreflecting boundary condition for discrete acoustic and elastic wave equations , 1985 .

[5]  B. Luke,et al.  Improved parameterization to invert Rayleigh-wave data for shallow profiles containing stiff inclusions , 2007 .

[6]  A. Levander.  Fourth-order finite-difference P-SV seismograms , 1988 .

[7]  Qing Huo Liu,et al.  PERFECTLY MATCHED LAYERS FOR ELASTODYNAMICS: A NEW ABSORBING BOUNDARY CONDITION , 1996 .

[8]  Jianghai Xia,et al.  Generating an Image of Dispersive Energy by Frequency Decomposition and Slant Stacking , 2007 .

[9]  Jianghai Xia,et al.  Numerical investigation of implementation of air-earth boundary by acoustic-elastic boundary approach , 2007 .

[10]  Jean-Pierre Berenger,et al.  A perfectly matched layer for the absorption of electromagnetic waves , 1994 .

[11]  Patrick Joly,et al.  Stability of perfectly matched layers, group velocities and anisotropic waves , 2003 .

[12]  Thorne Lay,et al.  Quantitative Seismology, Second Edition , 2003 .

[13]  Chrysoula Tsogka,et al.  Application of the PML absorbing layer model to the linear elastodynamic problem in anisotropic hete , 1998 .

[14]  Jianghai Xia,et al.  Estimation of Elastic Moduli in a Compressible Gibson Half-space by Inverting Rayleigh-wave Phase Velocity , 2006 .

[15]  Sebastiano Foti,et al.  Surface-wave analysis for building near-surface velocity models — Established approaches and new perspectives , 2010 .

[16]  J. Kristek,et al.  3D Heterogeneous Staggered-grid Finite-difference Modeling of Seismic Motion with Volume Harmonic and Arithmetic Averaging of Elastic Moduli and Densities , 2002 .

[17]  Jianghai Xia,et al.  Comparing shear-wave velocity profiles inverted from multichannel surface wave with borehole measurements , 2002 .

[18]  C. Tsogka,et al.  Application of the perfectly matched absorbing layer model to the linear elastodynamic problem in anisotropic heterogeneous media , 2001 .

[19]  J. Carcione,et al.  Seismic modelingSeismic modeling , 2002 .

[20]  Peter Moczo,et al.  Efficient Methods to Simulate Planar Free Surface in the 3D 4th-Order Staggered-Grid Finite-Difference Schemes , 2002 .

[21]  L. Knopoff,et al.  Fast Surface Wave and Free Mode Computations , 1972 .

[22]  R. Madariaga.  Dynamics of an expanding circular fault , 1976, Bulletin of the Seismological Society of America.

[23]  D. Komatitsch,et al.  An unsplit convolutional perfectly matched layer improved at grazing incidence for the seismic wave equation , 2007 .

[24]  Rune Mittet,et al.  Free-surface boundary conditions for elastic staggered-grid modeling schemes , 2002 .

[25]  W R Hudson,et al.  DETERMINATION OF IN SITU SHEAR WAVE VELOCITIES FROM SPECTRAL ANALYSIS OF SURFACE WAVES , 1982 .

[26]  Jianghai Xia,et al.  Estimation of near‐surface shear‐wave velocity by inversion of Rayleigh waves , 1999 .

[27]  Kristel C. Meza-Fajardo,et al.  A Nonconvolutional, Split-Field, Perfectly Matched Layer for Wave Propagation in Isotropic and Anisotropic Elastic Media: Stability Analysis , 2008 .

[28]  J. Virieux.  P-SV wave propagation in heterogeneous media: Velocity‐stress finite‐difference method , 1986 .

[29]  Erik H. Saenger,et al.  Finite-difference modeling of viscoelastic and anisotropic wave propagation using the rotated staggered grid , 2004 .

[30]  É. Delavaud,et al.  Interaction between surface waves and absorbing boundaries for wave propagation in geological basins: 2D numerical simulations , 2005 .

[31]  Gaetano Festa,et al.  PML Absorbing Boundaries , 2003 .