快速Fourier变换波动方程基准面校正方法研究

详细信息   

• 责任者：卢回忆^c中国科学院地质与地球物理研究所岩石圈演化国家重点实验室^d1980^e博士研究生 ; 符力耘 ; 蒋韬
• 刊名：地球物理学进展
• 出版年：2010
• 卷期：25^b4总102
• 关键词：波动方程 ; 静校正20^a波场延拓 ; 速度分析 ; 模型 ; 算法 ; 傅里叶分析30^a基准面校正 ; 复杂地表 ; 快速Fourier变换
• 页码：1313～1322
• 索书号：P260.6/208-1
• 图表：^c表1参35插图
• 分类号：2100
• 文摘期号：2011/04

摘要

当地表起伏剧烈、近地表速度横向变化较大时，基于地表一致性假设的常规静校正方法存在着较大误差。波动方程基准面静校正方法能很好地解决起伏地表和复杂近地表结构问题，但计算量巨大，特别是三维波动方程基准面校正，适应横向任意速度变化、计算精度较高的有限差分或其混合的方法波动方程基准面校正涉及海量的计算和存储操作。为了提高波动方程基准面校正的计算效率，该文研究一类只用快速Fourier变换(FFT)实施波动方程基准面校正的方法，采用相移(PS)、分裂步(SSF)和一阶退化(DP1)三种具有相同算法结构、但不同计算效率、适应不同地表复杂程度的Fourier变换延拓算子。PS和SSF算子只适应于速度横向变化较弱的起伏地表;DP1通过在两个分裂步之间作波数域线性插值来实现波场延拓，将常规的SSF算法推广适应强速度横向变化介质和大角度传播波场。该文着重比较了基于这三种延拓算子的逐层延拓累加波动方程基准面校正方法对地表起伏和近地表速度横向变化的适应能力和计算效率，给出了一个相对定量的评估，以便针对不同的地表复杂程度合理选择合适的FFT波动方程基准面校正方法，既满足了精度又提高了计算效率。