地震波场数值模拟是地震学领域中一项非常重要的技术方法,其对特定地质模型做适当的抽象并将其简化为数学模型,采用数值计算方法获取地震响应。地震波场数值模拟的目标就是在地下地质体特征已知的情况下,通过计算得到理论时深关系和时间切片图。其对于地震数据的采集、地震资料的处理和解释具有极大的辅助作用。地震波场数值模拟的方法主要包括两大类,即波动方程法和射线追踪法。波动方程数值模拟方法实质上是求解地震波波动方程,其模拟的地震波场中包含了地震波传播的所有信息。射线追踪法是将地震波波动理论简化为射线理论,它主要考虑的是地震波传播的运动学特征。由于波动方程法是通过对计算区域网格化来对方程进行求解,因而其计算量比较大,运算速度较慢,此外该方法对计算机内存要求也非常高。故本文使用射线追踪法对地震波场进行数值模拟。传统的射线追踪法主要解决的是P波的射线路径问题以及其对应的旅行时问题,而对于整个地震波场的地震波衰减以及多次波、转换波等问题都没有进行考虑。相较于波动方程数值模拟方法,射线追踪法在获取较高运算速度的同时损失掉了地震波场的动力学信息。本文在使用射线追踪法进行波场模拟的同时拟加入衰减信息以及追踪转换波以弥补射线追踪法在动力学上的缺陷。论文研究层状弹性介质模型,主要考虑了两方面的衰减：透射损失以及球面扩散。在透射损失方面,主要引用了Zoeppritz方程计算反射系数以及透射系数;在球面扩散方面,论文则结合水平层状介质对球面扩散因子进行了推导。在射线追踪方面,本文结合试射法以及弯曲迭代法完成传统的P波射线追踪外,考虑到了P波以及SV波在非垂直入射的情况下会产生转换波。因而对每条非垂直入射的射线同时计算其透射纵波、反射纵波、透射横波以及反射横波,进行全波射线追踪。同时还利用该算法进行了多次波的追踪。全波射线追踪改进了传统的射线追踪的算法,糅合使用试射法与弯曲迭代法确定射线入射点,同时追踪同类波与转换波,一次波与多次波,引入了“四叉树”对各类波进行数据管理,丰富了所追踪地震波的类型。在全波射线追踪过程中,于完成射线路径追踪以及旅行时计算的同时,对各地层间的射线匹配其反射系数、透射系数以及球面扩散因子,完成该段射线的能量衰减的计算。在得到的合成地震记录中涵盖了纵横波信息、同类波与转换波信息。这样在取得较高计算速度的同时也获取了较为完整的合成地震记录,在模型计算中取得了较好的效果。