随着勘探对象的日益复杂化,使得研究的模型介质也越来越复杂,为了更详细地了解和认识复杂介质波场,进行地震波场数值模拟是非常必要的。波动方程有限差分法是近几年最流行和最受欢迎的数值模拟方法之一,其优点就是模拟波场信息丰富、实现简单和运算速度较快等。波动方程有限差分法虽然有很多优点,但是也存在一些不足之处,其中稳定性条件就是一个很重要且急需解决的问题。稳定性条件直接决定正演模拟结果的成败,在粘弹介质波动方程的有限差分解中,介质的粘弹特性会增加数值解的不稳定性。本文从粘弹性理论入手,研究了Kelvin-Voigt和Maxwell在不同精度差分格式下的稳定性条件,而对于标准线性粘弹模型,只给出了稳定性条件与速度和品质因子的数值解。当空间网格步长一定时正演模拟中的时间步长是由模型中的最大速度决定的,速度越大时间步长越小,所需的运算量也越大,为了在不影响结果的前提下提高波动方程有限差分法模拟的效率,本文提出了波动方程变时间步长有限差分数值模拟。近地表的主要特点就是低速度、低密度和低品质因子。随着山地、沙漠等地表复杂区地震勘探的发展,近地表地震数值模拟技术受到了地球物理学家的广泛关注。结合上述得到的结果进行地震波场近地表效应的模拟及分析,对这一部分的模拟采用自由边界条件,对近地表的瑞雷面波进行了模拟和定量分析,确定模拟出的是瑞雷面波。粘弹性介质是更符合实际介质,进行了含有低速层的粘弹性介质在自由边界条件下的数值模拟。还从Q值大小、震源子波频率等方面定量讨论粘弹介质地震波衰减影响因素。实际介质(和油气藏有关的介质)往往是非均匀的,非均匀介质中存在很多微小异常且分布极不规则(通过测井数据和岩心样品就可以发现)。对于上述情况实现了随机介质的构造和模拟,得到非均匀介质中的波场。最后对几个实际速度模型进行数值模拟,得到了复杂储层(含有孔洞、裂缝等)下的全波场和储层在不同含水饱和度下的波场,并且通过频谱和时频分析对结果进行了分析,还研究了低速层和近地表对波场模拟的影响。