基于多倍角公式的一步波场外推法  

One-step extrapolation method based on the multiple-angle formula

作　　者：柯璇[1,2] 石颖 王银凤[5,6] KE Xuan;SHI Ying;WANG YinFeng(School of Earth Sciences,Northeast Petroleum University,Daqing 163318,China;Key Laboratory of Oil and Gas Reservoir Formation Mechanism and Resource Evaluation in Heilongjiang Province,Daqing 163318,China;Institute of Unconventional Oil and Gas,Northeast Petroleum University,Daqing 163318,China;Key Laboratory of Continental Shale Hydrocarbon Accumulation and Efficient Development,Ministry of Education,Northeast Petroleum University,Daqing 163318,China;College of Geo-Exploration Science and Technology,Jilin University,Changchun 130021,China;School of Mathematics and Statistics,Northeast Petroleum University,Daqing 163318,China)

机构地区：[1]东北石油大学地球科学学院,大庆163318 [2]黑龙江省油气藏形成机理与资源评价重点实验室,大庆163318 [3]东北石油大学非常规油气研究院,大庆163318 [4]东北石油大学陆相页岩油气成藏与高效开发教育部重点实验室,大庆163318 [5]吉林大学地球探测科学与技术学院,长春130021 [6]东北石油大学数学与统计学院,大庆163318

出　　处：《地球物理学报》2021年第7期2480-2493,共14页Chinese Journal of Geophysics

基　　金：国家自然科学基金项目(41804133,41930431);东北石油大学优秀科研人才培育基金(GLJHB201601);中央支持地方高校改革发展资金人才培养支持计划项目(140119001);中国博士后科研项目(2019M651253);黑龙江省博士后科研项目(LBH-Z18040);黑龙江省油气藏形成机理与资源评价重点实验室开放基金项目(KL20190105);大庆市指导性科技计划项目(zd-2019-02);东北石油大学人才引进科研启动经费项目(2019KQ05)联合资助.

摘　　要：为了提高地震波场正演模拟的准确性和稳定性,针对一步波场外推法地震波场正演,本文提出了基于多倍角公式的耦合方程组解法.借助欧拉公式,将一步波场外推法的复数波场延拓方程转化为两个实数波场耦合的方程组,结合多倍角公式和泰勒展开式精确逼近包含拟微分算子的简谐函数算子,利用谱方法求解拟微分算子,进而推导了一种基于多倍角公式的一步波场外推法的耦合方程组.相比于常规一步波场外推法中复数方程的矩阵解法,本文方法能够显著减少傅里叶变换次数,降低计算成本.此外,本文推导了稳定性条件,为正确选取地震波场模拟参数提供了理论依据.基于二维匀速模型和复杂构造模型的数值测试表明,本文方法能够在大时间步长情况下保持外推波场稳定,计算效率较高.We propose a coupled-equations solution for the one-step extrapolation method to solve the acoustic wave equation.We use the Euler formula to transform the one-step extrapolation equation for a complex wavefield into coupled equations of two real wavefields.Then this multiple-angle formula is applied to more accurate and stable approximations of the harmonic function operators of pseudodifferential operators.The polynomial Taylor expansions are adopted in the coupled equations and then the Fourier method is used to solve this operator and we derive a kind of coupled equations of the one-step extrapolation method.The method workflow requires fewer iterations of the Fourier transform than the existing one-step wave extrapolation matrix method.Besides,we also derive the stability conditions,which could provide a quantitative basis for choosing the multiple-angle parameters and the orders of polynomials to ensure the proposed algorithm stable.Numerical experiments for a two-dimensional constant-velocity model and complex structural model show that our method can keep stable of the extrapolation of the wavefield with a large time step and is more effective than the one-step wave extrapolation matrix method.

关 键 词：一步外推法 耦合方程组 多倍角公式 

分 类 号：P631[天文地球—地质矿产勘探]