论文摘要
地震数值模拟在地震勘探中具有重要作用。按照原理可分为基于射线理论的射线追踪方法和基于波动理论的波动方程正演。射线追踪在复杂地表、复杂构造中的反射波追踪是一个难点。本文通过对试射法、全局最短路径法、常速度梯度法的研究,取得了一定的效果。首先,本文对试射法推导了透射线和反射线的传播方向,并将其应用到起伏地表和复杂构造模型中。模型试算表明,算法稳定,可适应复杂介质模型。然而,试射法无法适应强速度变化、难于求出多值走时中的全局最小走时,为此,研究了全局最短路径法。根据图形理论和惠更斯原理,求得任意点的全局最短路径,同时实现了反射波追踪,并对误差进行了分析,算例分析表明,全局最短路径法能够适应强速度变化,可以求出全局最小走时,但只适合两点射线追踪,且存在射线分叉问题。通过研究常速度梯度法基本原理,采用网格剖分的方法,求解一组射线路径的解析表达式,得到了射线追踪路径,结合斯奈尔定律,实现了反射波射线追踪。计算反射波路径时需要有层位信息,常规的层位结构模型,无法适应界面起伏较大或存在许多逆断层的复杂模型,本文提出一种线性连接法和分段结构模型,使得层位模型的建立简洁方便,并将其应用到常速度梯度法中取得了良好的效果。当模型为复杂介质时,采用普通射线法制作合成地震记录,会存在焦散区、空白区等奇异性区域。本文采用高斯射线束代替普通射线进行正演模拟,克服了普通射线法的缺陷,通过模型试算,证明了方法的正确性,且有较高的精度。射线追踪计算速度快,但无法很好的反映波场的动力学特征。因此,本文对波动方程正演中的非均匀介质二维声波、粘声波、弹性波以及粘弹性波交错网格高阶有限差分正演模拟进行了研究,并以复杂介质模型进行了试算。探讨了三维声波、弹性波高阶交错网格有限差分正演模拟,并以简单层状模型进行了试验。对有限差分中涉及的边界条件、稳定性及频散进行了分析。给出了透明边界、旁轴近似边界和PML边界。对前人得到的几种稳定性条件进行了分析和试算,得到了适于本文的稳定性条件。通过理论推导,对声波有限差分引起频散的几个参数进行了定量分析。从CRP叠加入手,将射线追踪方法应用于观测系统设计,分析每道覆盖次数,优化观测系统参数。同时,采用照明分析对波动方程在观测系统中的应用进行了探讨。最后对正演软件的研发进行了介绍和应用。