Application of a perfectly matched layer in seismic wavefield simulation with an irregular free surface

Abstract

A B S T R A C TRecently, an effective and powerful approach for simulating seismic wave propagation in elastic media with an irregular free surface was proposed. However, in previous studies, researchers used the periodic condition and/or sponge boundary condition to attenuate artificial reflections at boundaries of a computational domain. As demonstrated in many literatures, either the periodic condition or sponge boundary condition is simple but much less effective than the well-known perfectly matched layer boundary condition. In view of this, we intend to introduce a perfectly matched layer to simulate seismic wavefields in unbounded models with an irregular free surface. We first incorporate a perfectly matched layer into wave equations formulated in a frequency domain in Cartesian coordinates. We then transform them back into a time domain through inverse Fourier transformation. Afterwards, we use a boundary-conforming grid and map a rectangular grid onto a curved one, which allows us to transform the equations and free surface boundary conditions from Cartesian coordinates to curvilinear coordinates. As numerical examples show, if free surface boundary conditions are imposed at the top border of a model, then it should also be incorporated into the perfectly matched layer imposed at the top-left and top-right corners of a 2D model where the free surface boundary conditions and perfectly matched layer encounter; otherwise, reflections will occur at the intersections of the free surface and the perfectly matched layer, which is confirmed in this paper. So, by replacing normal second derivatives in wave equations in curvilinear coordinates with free surface boundary conditions, we successfully implement the free surface boundary conditions into the perfectly matched layer at the top-left and top-right corners of a 2D model at the surface. A number of numerical examples show that the perfectly matched layer constructed in this study is effective in simulating wave propagation in unbounded media and the algorithm for implementation of the perfectly matched layer and free surface boundary conditions is stable for long-time wavefield simulation on models with an irregular free surface.

An upwind fast sweeping scheme for calculating seismic wave first‐arrival travel times for models with an irregular free surface

“…In recent numerical simulation of wave propagation, effective attempts have been made to deal with surface topography taking the advantage of transformation from Cartesian to curvilinear coordinates (Hestholm and Rudd ; Tessmer and Kosloff ; Appelo and Petersson ; Lan and Zhang ; Lan et al . ). We extended this idea to travel‐time computation and developed a Lax–Friedrichs (LF) fast sweeping scheme to calculate seismic wave first‐arrival travel times for models with rugged surfaces (Lan and Zhang ,b).…”

A novel hybrid method based on discontinuous Galerkin method and staggered‐grid method for scalar wavefield modelling with rough topography

“…The finite‐difference method (FDM) is most popular (Alford et al ., 1974; Lisitsa et al ., 2010; Zhang & Chen, 2010; Virieux et al ., 2011; Lan et al ., 2015). Due to programming simplicity and low storage capacity, FDM has been widely used in numerical modelling (Wang & Liu, 2007; Liu & Sen, 2009; Lisitsa & Vishnevskiy, 2010; Lisitsa et al ., 2012), reverse time migration (Herrmann & Li, 2012; Zhao & Sen, 2019) and inversion problem (Ren & Liu, 2015).…”