在当今世界各地,尤其是在中国的西部地区,大量与石油、天然气及矿产等资源相关的地震勘探工作在起伏地表地区进行。与平原地区的地震勘探相比,在山地开展地震勘探要面临一些特殊问题,例如：数据采集困难且采集到的数据质量差、散射干扰严重、静校正不准、成像效果不好甚至不成像等。造成这些问题的主要原因有：①起伏地表地区复杂的地质条件造成地震波场的结构也很复杂,而长期以来人们对这一特定区域地震波的传播规律认识不够深入,很多理论和实际问题还没有得到很好的解决；②由于地质条件的复杂和缺乏大量的正演模拟数据的支持,起伏地表地区地震数据采集的有效性和针对性相对较差；③传统的地震勘探理论和处理技术都是基于地表为水平的假设,对于起伏地表问题常采用静校正来处理,然而在地表剧烈起伏地区即使采用了细致的静校正技术也很难取得好的成像效果。解决上述问题的一个基本途径是增加对地震波场的认识,其主要方法有两种,即波动方程数值模拟和射线追踪。本文选择后者来研究该问题,主要研究对象是起伏地表条件下的地震波走时与射线路径的计算问题,主要研究内容包括：①为了给起伏地表条件下的地震波走时与射线路径计算寻找合适的算法,首先对常规走时与射线路径的计算方法进行了分析及对比：分别阐述了两点射线追踪法、最短路径射线追踪法、有限差分法、线性插值法、波前构建法以及其它方法的基本原理、发展历程等问题,并基于相关文献对这些算法进行了综合对比,最后基于对比的结果选择快速推进法和线性插值法作为本文研究问题的算法基础。②针对二维起伏地表条件下地震波走时的计算问题,对常规快速推进法进行改进,提出了阶梯网格迎风差分法和不等距差分法：阶梯网格迎风差分法运用阶梯网格建立起伏地表模型,采用改进后的迎风差分格式进行局部走时计算。不等距差分法运用不等距网格建立起伏地表模型,通过在迎风差分格式中引入不等距差分格式并综合应用Huygens原理和Fermat原理进行局部走时计算。两种方法在整体实现时均采用起伏地表条件下的窄带技术作为波前扩展方式。计算精度分析表明两种算法均能达到很好的计算精度。③为了计算二维起伏地表条件下的地震波走时,对常规线性插值法进行改进,提出了一种混合网格线性插值法：对常规线性插值法的局部实现策略和波前扩展方式均进行一定程度的改进,然后采用混合网格建立起伏地表模型,采用三角网格和正方形网格中的插值公式进行局部走时计算,并利用起伏地表条件下的窄带技术作为波前扩展方式。计算精度分析表明混合网格线性插值法能够达到很好的计算精度。④对二维起伏地表条件下的曲网格法进行了初步的研究：推导了坐标变换法和正交贴体网格法在计算空间中的程函方程,给出了求解这两个方程的数值实现策略,并对坐标变换法进行了初步的计算精度分析,最后讨论了曲网格法存在的一些技术难点。⑤对②、③、④中提出的几种二维起伏地表条件下的地震波走时算法进行了对比分析,并给出了相应的计算实例：首先对常规快速推进法和线性插值法进行精度和效率的对比,然后对笛卡尔坐标系下的三种算法进行对比,最后对曲网格法和笛卡尔坐标系下的算法进行了对比。⑥提出二维起伏地表条件下的地震波射线路径的计算方法,并对算法进行了相应的计算精度、稳定性及模型的适应性分析：起伏地表条件下的地震波射线路径计算是从已计算出的起伏地表条件下的地震波走时分布信息出发利用线性插值法来完成的。通过算法分析和具体的计算实例表明算法能达到很好的计算精度,同时算法也能稳定地适应任意起伏地形和复杂介质模型。⑦提出三维起伏地表条件下的地震波走时计算方法,并给出了算法的稳定性及计算精度分析和计算实例：采用不等距网格建立三维起伏地表模型,推导了三维不等距网格中的不等距差分公式。算法的稳定性、计算精度及计算实例分析表明基于不等距差分法的三维起伏地表条件下的地震波走时计算方法简单易行且能达到很好的计算精度,同时也能稳定地适应任意的三维起伏地形和复杂介质模型。⑧利用上述③、⑥中提出的起伏地表条件下的地震波走时与射线路径的计算方法,有针对性地计算了起伏地表条件下各种波型的走时与射线路径,并研究了起伏地表条件下的地震波的一些传播规律：利用③中提出的混合网格线性插值法和⑥中提出的射线路径计算方法分别计算了起伏地表条件下初至波、首波、透射波、绕射波、反射波以及多次反射波的走时与射线路径,通过对这些波的走时与射线路径分布情况的分析得出了一些起伏地表条件下地震波的传播规律。基于上述研究内容,本论文主要取得了如下研究成果：①提出了一套以阶梯网格迎风差分法、不等距差分法以及混合网格线性插值法为主体的二维起伏地表条件下的地震波走时计算方法,这些方法均能适应二维任意起伏地形和任意复杂介质且兼顾计算精度及效率；②引入曲网格法解决起伏地表问题,建立了曲网格中用于走时计算的基本方程,并提出了相应的数值实现策略,同时还分析了算法在实现时还有待解决的一些技术难点；③提出了一种基于线性插值的起伏地表条件下的射线路径计算方法,该方法能在保证一定计算精度的条件下适应二维任意起伏地形和任意复杂介质；④提出了一种基于不等距差分法的三维起伏地表条件下的地震波走时计算方法,该方法能在保证一定计算精度的条件下适应三维任意起伏地形和任意复杂介质；⑤提出了起伏地表条件下的各种波型的走时与射线路径的计算方案,并得出了一些起伏地表条件下地震波的传播规律。综上所述,本论文所做的研究工作和所得出的一些结论及研究成果对于起伏地表地区地震勘探效果的提高有着一定的理论意义和实际价值,同时也有着很广阔的应用前景。