Seismic modeling plays an important role in geophysics and seismology for estimating the response of seismic sources in a given medium. In this work, we present a MATLAB-based package, FDwave3D, for synthetic wavefield and seismogram modeling in 3D anisotropic media. The seismic simulation is carried out using the finite-difference method over the staggered grid, and it is applicable to both active and passive surveys. The code package allows the incorporation of arbitrary source mechanisms and offers spatial derivative operators of accuracy up to tenth-order along with different types of boundary conditions. First, the methodological aspects of finite-difference method are briefly introduced. Then, the code has been tested and verified against the analytical solutions obtained for the homogeneous model. Further, the numerical examples of layered and overthrust models are presented to demonstrate its reliability.

地震建模在地球物理学和地震学中起着重要作用，用于估计给定介质中地震源的响应。在这项工作中，我们提出了一个基于 MATLAB 的软件包FDwave3D，用于 3D 各向异性介质中的合成波场和地震图建模。地震模拟是在交错网格上使用有限差分方法进行的，它适用于主动和被动勘测。该代码包允许结合任意源机制，并提供精度高达十阶的空间导数算子以及不同类型的边界条件。首先，简要介绍有限差分法的方法论方面。然后，针对齐次模型获得的解析解对代码进行了测试和验证。此外，还给出了分层模型和逆推模型的数值例子，以证明其可靠性。