Frequency-Domain Finite Element Modeling of Seismic Wave Propagation Under Different Boundary Conditions

Seismic wave propagation in complex terrains, especially in the presence of air layers, plays a crucial role in accurate subsurface imaging. However, the influence of different boundary conditions on seismic wave propagation characteristics has not been fully explored. This study employs the finite element method (FEM) to simulate and analyze seismic wavefields under different boundary conditions, including perfectly matched layer (PML), Neumann free boundary conditions, and air layer conditions. First, the finite element solution for the 2D frequency-domain acoustic wave equation is introduced, and the correctness of the algorithm is validated using a homogeneous model. Then, both horizontal and undulating terrain interfaces are designed to investigate the kinematic and dynamic characteristics of the wavefields under different boundary conditions. The results show that PML boundaries effectively absorb seismic waves, prevent reflections, and ensure stable wave propagation, making them an ideal choice for simulating open boundaries. In contrast, Neumann boundaries generate significant reflected waves, particularly in undulating terrains, complicating the wavefield characteristics. Introducing an air layer alters the dynamics of the wavefield, leading to energy leakage and multi-path effects, which are more consistent with real-world seismic-geophysical models. Finally, the computational results using the Overthrust model under different boundary conditions further demonstrate that different boundary conditions significantly affect wavefield morphology. It is essential to select appropriate boundary conditions based on the specific simulation requirements, and boundary conditions with an air layer are most consistent with real seismic geological models. This study provides new insights into the role of boundary conditions in seismic numerical simulations and offers theoretical guidance for improving the accuracy of wavefield simulations in realistic geological scenarios.