在现代地球物理勘探中,目标体越来越复杂,所需的精度和分辨率越来越高。普通地面地震因为地表低降速带的影响,产生的记录和井中地震相比,精度较低。此外,井中的地震频率更高,可以实现更好的纵向分辨率。为了理解在复杂地质条件下形成的波场信息,需要一种精确的正演方法。对于复杂的地质体来说,非规则网格剖分十分适用,它可以根据模型的地质参数来调整剖分的大小,使之更加符合地质条件。在适用于不规则网格划分的正演方法中,最简单的方法是有限元方法。本论文从实际出发,为了解决复杂地质条件和起伏地表条件下的正演模拟问题,运用Delaunay三角网剖分技术,将地下复杂的构造和带有起伏界面的地质模型剖分为Delaunay三角形,并运用有限单元法求解在该区域内的弹性波波动方程,得出在相同网格情况下,非规则三角网能够更好的模拟起伏地层和起伏界面,求取的波场更为精细。论文在有限元求解过程中,运用逐个单元法代替整体组装矩阵,有效的降低了有限元方法的内存占用率。由于逐个单元法容易实现并行化处理,论文将程序用Openmp加速,实现计算速度的提升。在压制频散时采用集中质量矩阵和分散质量矩阵组合的方法,有效压制了频散,但是这导致质量矩阵不是对角矩阵,因此,论文用前置共轭梯度迭代加速求解线性方程组。同时在边界的处理上,论文改进了Sarma吸收边界,用二次递增分段函数函数代替一次函数,降低了边界人工反射。论文在模型正演部分,从简单到复杂,包括了起伏地层界面、起伏地表、各向同性介质、各向异性介质等,同时正演了VSP记录和井间记录,并且对记录中的波场信息进行了详细的说明。