Application of the Multiaxial Perfectly Matched Layer to Near‐Surface Seismic Modeling with Rayleigh Waves

Chong Zeng,Jianghai Xia,Richard D. Miller,G. P. Tsoflias
DOI: https://doi.org/10.4133/1.3614264
2011-01-01
Abstract:PreviousNext You have accessSymposium on the Application of Geophysics to Engineering and Environmental Problems 2011Application of the Multiaxial Perfectly Matched Layer to Near‐Surface Seismic Modeling with Rayleigh WavesAuthors: Chong ZengJianghai XiaRichard MillerGeorgios TsofliasChong Zeng, Jianghai Xia, Richard Miller, and Georgios Tsofliashttps://doi.org/10.4133/1.3614264 SectionsAboutPDF/ePub ToolsAdd to favoritesDownload CitationsTrack CitationsPermissions ShareFacebookTwitterLinked InRedditEmail Abstract 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 takes values greater than about 0.38. 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 non‐linear 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 in orthogonal directions. If the proportion coefficients of the damping profiles 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 find that a constant proportion coefficient of 1.0 is sufficient to stabilize the M‐PML for all 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 loosing its stability. The simplified M‐PML works efficiently for both homogeneous and heterogeneous earth models with high Poisson's ratios. The analysis reported on here was based on 2‐D finite difference modeling in the time domain that can easily be extended into the 3‐D domain with other numerical methods.Permalink: https://doi.org/10.4133/1.3614264FiguresReferencesRelatedDetails Symposium on the Application of Geophysics to Engineering and Environmental Problems 2011ISSN (online):1554-8015Copyright: 2011 Pages: 642 publication data© 2011 Copyright © 2011 The Environmental and Engineering Geophysical SocietyPublisher:Environmental & Engineering Geophysical Society HistoryPublished Online: 06 Jul 2011 CITATION INFORMATION Chong Zeng, Jianghai Xia, Richard Miller, and Georgios Tsoflias, (2011), "Application of the Multiaxial Perfectly Matched Layer to Near‐Surface Seismic Modeling with Rayleigh Waves," Symposium on the Application of Geophysics to Engineering and Environmental Problems Proceedings : 131-131. https://doi.org/10.4133/1.3614264 Plain-Language Summary PDF DownloadLoading ...
What problem does this paper attempt to address?