固体地球物理研究生学术报告第四期（2）

The finite element method (FEM) and its variants such as spectral element method (SEM) have emerged as a powerful tool in seismic wave simulations when considering the model with topography and internal discontinuities. In these methods, the free boundary conditions can be naturally satisfied, and highly irregular topography can be treated directly, so they are suitable for high precision surface wave simulations. We present the mathematical background of the FEM and SEM in solving the seismic wave equation coupled with the Perfectly Matched Layer (PML) in the time domain. The workflows of FEM and SEM from generating the mesh to obtaining the waveforms will be discussed in both regional and global scales. A novel Multi-Axial Perfectly Matched Layer (M-PML) developed in our recent research will be discussed. The proposed M-PML successfully achieved stable simulations in some anisotropic medium where the conventional PML fails, and our method saves more memory and remains computational efficient than the M-PML in literatures.