^Doherty,2001;张永刚,2007;LanandZhang,2011a,0口2012;MaandZhang,2014a,2014b)与矿产资源地球的表面常常是起伏或是不规则的,且实际的地球物(Farquharson以a/.,2008;Lelievreetai.,2009,2012)地理观测绝大部分是在这些地表崎岖的地方进行的球物理勘探中,地球物理学家也遇到同样的问题.剧烈的地(RobertssonandHolliger,1997).例如,深地震测深中为了形起伏给地震探测工作提出了严峻的挑战.理解造山带、盆地等的形成机制,测线常常穿过这些地表剧起伏地表下的地震波传播数值模拟,不仅是对起伏地表烈起伏的地区(Teng过,1987,2003;Zengetal.,1995;观测的地震资料进行偏移和反演的基础,而且对后续处理结LiandMooney,1998;Gaoetal.,2000,2005;Zhangand果的解释具有重要意义.地震波场和地震波走时场作为描述Klemperer,2010).在油气(阎世信等,2000;Yilmazand地震波传播的两种重要途径,常常通过对它们的模拟来研究
地震波的传播规律和传播特征.下面我们就从地震波场和地提高,提出了一些有针对性的地震波场数值模拟方法.包括震波走时场两方面来阐述起伏地表下的地震波传播数值模有限元法,伪谱法,边界元法,谱元法和有限差分法等.这些拟方法的研究进展.方法在模拟起伏地表的地震波场时各有优劣,下面对这些方,扣体诎志"CiVtW?键她忆粉估进H、丨七吐m由、ttH法逐一介绍,表1简单的概括了这些数值方法在进行复杂地1起舰表下的地震波场随模拟細开究进展表的波场觀时的优点和缺点.近年来,随着计算机技术的‘决速发展和计算能力的大幅表1起伏地表下的地震波场数值模拟方法Table1Methodsforsimulatingseismicwavefieldsinmediawithanirregularsurface数值模拟方法^^有限元法处理不规则地表方便耗内存大,计算速度慢,成本高谱元法处理不规则地表方便成本高德在怠上田出左,丨、处理地表结构复杂或地表剧烈起伏以及地下结构复杂的情况时雜歡的误差边界元法处理不规则地表方便不适用于地表速度变化较大情况有限差分法处理复杂地形比細难1.1有限元法(FE)维横向各向同性介质波动方程的正演模拟,特别是对边界吸有限元法基于变分原理,采用分段近似,基于三角形网收问题作了有效的处理.Funmura等(1995)首次针对伪谱法格来剖分模型,该剖分保证了复杂地层形态模拟的逼真性提出了反周期扩展边界方法;Takenaka和王彦宾等(Rialetal.,1992;ToshinawaandOhmachi,1992).但有限(TakenakaeZaZ.,1999,Wangetal.,2000)利用伪谱法分元法算法复杂,计算速度慢,占用内存和运算量均较大.为别计算了球对称全球模型和具有垂向速度梯度的沉积盆地此,发展了一些有限元和其他方法结合的混合方法.例如,模型中地震波的传播问题.之后不久,WangandTakenakaMoczo等用离散波数方法模拟震源激发和下部介质中地震(2001)利用不连续网格傅立叶伪谱多域方法模拟了区域地波的传播,而通过有限元方法来模拟沿起伏地表传播的波球模型中弹性波的传播,接下来Wang等(2001)利用伪谱法(Moczoetal.,1997;Galisetal.,2008);张美根等(2002)模拟了二维柱坐标下全地球模型中地震波的传播.赵志新等研究了各向