地球介质的各向异性是普遍存在的,地震各向异性研究已经成为地震学研究领域中的前沿课题之一。地震波正演模拟是模拟地震波在地球介质中的传播过程,并研究地震波的传播特性与地球介质参数的关系,通过正演模拟以达到对实际观测地震记录的最优逼近。开展地震波的正演模拟研究,对人们正确认识地震波的传播规律,验证所求地球模型的正确性;进行实际地震资料的地质解释与储层预测以及地球资源开发等,均具有重要的理论和实际意义。本文对具有垂直对称轴的三维横向各向同性(VTI)介质qP波在频率域和时间域进行了正演模拟。在频率域正演中,从三维VTI介质频率空间域qP波波动方程出发,在常规差分算子的基础上构造了适合三维VTI介质的频率空间域加权平均有限差分算子,利用Gauss-Newton优化方法算出VTI介质加权平均差分算子系数,并对常规差分算子和加权平均差分算子进行了速度频散精度对比分析。本文把螺旋边界和吸收边界应用在频率域正演上,在求取大型带状稀疏差分矩阵方面,作者提出了改进算法;并给出了三维频率域正演的数值实例。在三维VTI介质qP波时间域正演模拟方面,首先对三维VTI介质qP波时间域波动方程进行降阶处理,接着利用提出的加权平均差分算子对波动方程进行差分,并对加权平均差分算子进行了速度频散精度对比分析,还提出了一种新的分裂的完全匹配层来处理时间域边界,并进行了边界处理效果对比分析,最后进行了三维时间域正演模拟。频率域和时间域数值模拟表明:本文提出的算法有效地压制了数值频散,成功地减少了人为边界反射,具有较高的模拟精度,为三维VTI介质qP波的偏移和反演奠定了理论基础。