摘要
数值模拟方法在研究地震波波场传播和地下构造情况中起着重要的作用。在实际应用中,三维地震波场交错网格数值模拟方法应用广泛,但是模拟计算需要消耗大量的时间,计算成本很高。本文将变阶数方法应用到三维交错网格有限差分中,通过实际误差适配差分计算阶数,结合应力-速度方程实现了三维地震波场的正演模拟。利用GPU并行技术改进三维网格计算方法,并结合实际模型修改并行计算的模式实现了高性能运算,数倍提高了计算的效率。本文通过对三维实际模型中波场和地震记录的分析对比,并结合理论推导,验证了基于GPU的变阶数正演模拟方法应用于复杂地层三维地震波场数值模拟中的实用性和有效性。
引文
[1]Chen,H.,Wang,X.,and Zhao H..Arotated staggered grid finite-difference with the adsorbing boundary condition of a perfectly matched layer.Chinese Science Bulletin,2006,51(19):2304-2314.
[2]Liu,Y,and Sen,M.K..An imphcit staggered-grid finite-difference method for seismic modeling.Geophysical Journal International,2009,179:459-474.
[3]Wang,W Z.,Hu,T.Y,Lu,X.M.,Qin,Z.,Li,Y D.,and Zhang,Y.Variable-order rotated staggered-grid method for elastic-wave forward modeling.Applied Geophysics,2015,12(3):389-400.
[4]Sanders,J.,and Kandrot,E..CUDA by Example:An Introduction to General-purpose GPU Programming.U.S.A:Addison-Wesley Professional,2010.
[5]Farber,R..CUDA Application Design and Development.U.S.A:Elsevier Inc.,2011.
[6]Michea,D.,Komatitsch D..Acceleration a three-dimensional finite-difference wave propagation code using GPU graphics cards.Geophysical Journal International,2010,182:389-402.