Show simple item record

dc.contributor.authorZeng, Chong
dc.contributor.authorXia, Jianghai
dc.contributor.authorMiller, Richard D.
dc.contributor.authorTsoflias, Georgios P.
dc.date.accessioned2015-03-18T20:34:09Z
dc.date.available2015-03-18T20:34:09Z
dc.date.issued2011-05
dc.identifier.citationChong Zeng, Jianghai Xia, Richard D. Miller, and Georgios P. Tsoflias (2011). ”Application of the multiaxial perfectly matched layer (M-PML) to near-surface seismic modeling with Rayleigh waves.” Application of the multiaxial perfectly matched layer (M-PML) to near-surface seismic modeling with Rayleigh waves, 76(3), T43-T52. http://dx.doi.org/10.1190/1.3560019.en_US
dc.identifier.urihttp://hdl.handle.net/1808/17116
dc.descriptionThis is the published version. Reuse is subject to Society of Exploration Geophysicists terms of use and conditions.en_US
dc.description.abstractPerfectly 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.en_US
dc.publisherSociety of Exploration Geophysicistsen_US
dc.titleApplication of the multiaxial perfectly matched layer (M-PML) to near-surface seismic modeling with Rayleigh wavesen_US
dc.typeArticle
kusw.kuauthorXia, Jianghai
kusw.kuauthorMiller, Richard D.
kusw.kuauthorTsoflias, Georgios P
kusw.kudepartmentGeologyen_US
dc.identifier.doi10.1190/1.3560019
kusw.oaversionScholarly/refereed, publisher version
kusw.oapolicyThis item meets KU Open Access policy criteria.
dc.rights.accessrightsopenAccess


Files in this item

Thumbnail

This item appears in the following Collection(s)

Show simple item record