《2.5维地震波场褶积微分算子法数值模拟》.pdfVIP

第23卷 第4期 地球物理学进展 V01．23No．4 2008年8月(页码：1099～1105) PRoGRESSINGEOPHYSICS Aug．2008 2．5维地震波场褶积微分算子法数值模拟 程冰洁1， 李小凡2 (1．中石化西南分公司博士后科研工作站．成都610081}2．中国科学院地质与地球物理研究所，北京100029) 摘要早期的褶积微分算子都是基于正反傅立叶变换而实现的，其精度比四阶有限差分的精度稍高，本文将计算 数学中的Forsyte广叉正交多项式微分算予与褶积算子相结合，构建了一个新的快速、高精度褶积微分算子，其计算 结果非常接近实验函数微分的精确值，精度与16阶有限差分的精度相当，远优于错格伪谱法的精确度。另外，2。5维 数值模拟比二维模拟可以更真实地模拟三雏介质的某个剖面的波场，并且2．5维地震波模拟的计算量比三雏模拟的 计算量及计算耗时要大大减少．本文利用基于Forsyte广义正交多项式褶积微分算子法计算2．5维非均匀介质地震 波场，模拟结果表明，该算法的计算速度快，计算精度高，能够直观、高效地反映复杂介质中波场的传播规律，并且2．5 维波场数值模拟具有更高的计算效率，是一种jF常值得深入研究并广泛应用的方法． 关键词褶积微分算子，Forsyte广义正交多项式，地震波场，2．5维数值模拟 中图分类号P631 文献标识码A 文章编号 1004—2903(2008)04—1099—07 ofelastic the 2．5D waves convolutional modelling using differentiatormethod CHENG Xiao—fan2 Bing—jiel，LI PostdoctoralCenter，SouthwestBranch (1．Mobile Co．，SINOPEC，Chengdu610081，China’ 2．Institute and 100029，China) ofGeologyGeophysics，ChineseAcademyofScience，Beijing Abstraet convolutionaldifferentiatorswereallbasedontheFouriertransformation，theirwasalittle Early precision thanthatofthefour-orderfinitedifference．Forthe and oIseismic higer improvingprecisionefficiency modeling，the authors anew Differentiator developmodelingapproach(ConvolutionalForsytePolynomial Method)byusingopti— for in mizedconvolutionalfor differentiationand finite-differencetimedifferentiation