本文主要是通过波场数值正演模拟技术来研究粘弹性各向同性介质中地震波的传播问题。文中涉及粘弹性力学理论、地震波传播理论、有限差分理论、波场正演模拟技术等相关知识。所以首先介绍了粘弹性力学理论的基础知识，它是研究波在粘弹性介质中传播的理论基础，充分考虑了物质的粘弹性性质：蠕变、松弛、粘性、阻尼、内摩擦、滞弹性和弹性后效等等；其次根据粘弹性力学的基本理论介绍了粘弹性介质中地震波的传播问题，从而了解了地震波在介质中传播时能量变化同衰减系数和频率之间的关系，合理的解释了经典弹性波理论所不能解释的预期地震记录与实际接收的地震记录之间的差异；再次介绍了正演模拟用到的数值方法——有限差分方法的理论基础，包括常规有限差分法和多网格有限差分法。虽然数值模拟的方法有许多种，但是当精度要求一般时，有限差分技术是各种方法中实现最简单和最有效的方法，其中时间二阶和空间四阶的FD算法是一种好的选择(Carcione，2002)，尤其随着高阶交错网格有限差分法的发展与应用，我们可以根据需要对时间和空间做任意阶的有限差分，从而可以大大的提高计算精度、节约大量的计算时间；最后以前面各部分介绍的理论为基础，结合给出的地质模型进行波场正演数值模拟，并根据波场合成记录进行频谱分析，研究振幅谱随频率的变化、频率随偏移距的变化、频率随深度的变化。地震模拟技术是模拟地震波在地下传播的技术，其目的就是在假定地下结构模型的情况下，预测在给定的检波点生成的地震记录，该技术作为一项重要的检验手段，一方面可用于检验解释结果和各种反演算法的正确性，另一方面可用于设计和评估野外观测系统。本文采用的是全波方程法模拟，即由方程的解中给出了全波场的值。根据论文前部分所介绍的理论基础结合前人的研究成果推导了适合粘弹性介质中地震波传播的波动方程。又根据需要进一步简化速度－应力公式，使其便于用交错网格有限差分法进行数值模拟，推导了交错网格中时间和空间上的高阶差分近似公<WP=90>式，文中使用时间上二阶空间上四阶的交错网格有限差分数值方法。针对给出的地质模型(陆地地质模型和海洋地质模型)的特殊性，做了自由表面边界条件和流体－固体边界条件等特殊处理，成功的对两种地质模型进行了波场数值模拟。通过模拟得到了下面几点认识和结论：粘弹性介质充分了考虑了实际介质对地震波的吸收衰减作用，合理的解释了经典弹性波理论不能解释的预期地震记录同实际接收的地震记录之间的差异，使理想介质理论研究向实际介质理论研究前进了一步。交错网格有限差分法是一种高效的数值模拟方法，在提高计算速度和计算精度方面有着常规有限差分方法不可比拟的优势。随着高阶交错网格有限差分技术的应用，使得该方法的应用前景更为广阔。正演模拟技术为我们直观的了解地震波在复杂介质中的传播规律提供了方便。我们可以通过模拟得到的地震波传播快照和地震波场合成记录来获取地震波在介质中传播的信息。从而为勘探地震数据反演提供了必要的依据。从模拟得到的波场合成记录上可以清楚的看到粘弹性介质对地震波的吸收作用，但是对其具体影响，尚需通过AVO反射系数来分析，是需要以后进一步做的工作。通过海洋介质模型模拟得到的波场合成记录可以看出，当海水层较深时地震勘探受海洋多次波影响较小；当海水层较浅时受海洋多次波影响较严重。所以海洋多次波消除技术对海洋地震勘探有着深远的影响。振幅衰减特征：P波的衰减特征：零偏移距（即法向入射）时P波的振幅较强，随着偏移距的增大P波振幅减小。转换波的衰减特征：零偏移距时，没有转换波，随着偏移距的增大，出现转换波，当达到一定角度时，转换波的能量大于其衰减量，所以振幅增强，但随着偏移距的继续增大，其振幅减小。频率的衰减特征：振幅谱随频率的增大而减小，但不是单调递减；频率随深度的增大而减小，因为地层对频率有吸收作用，并且对高频的吸收作用比对低频的吸收作用强。