基于余弦调制Chebyshev窗的弹性波高精度正演  

High precision elastic wave equation forward modeling based on cosine modulated Chebyshev window function

作　　者：郑婉秋[1] 孟小红[1] 刘建红[2] 王建[1] 

机构地区：[1]中国地质大学(北京)地球物理与信息技术学院,北京100083 [2]中国石油集团东方地球物理勘探有限责任公司物探技术研究中心,河北涿州072751

出　　处：《地球物理学报》2016年第7期2650-2662,共13页Chinese Journal of Geophysics

基　　金：国家重大科研装备研制项目(ZDYZ2012-1-02-04);国家自然科学基金(41474106)联合资助

摘　　要：有限差分时间域正演是弹性波逆时偏移和全波形反演的基础,正演的计算精度也控制着偏移结果的准确性,若精度不高,则在偏移、反演后会带来假象.为了有效提高正演精度,本文结合窗函数优化方法,在窗函数截断伪谱法空间褶积序列以逼近有限差分算子的基础上,提出了一种基于Chebyshev窗的余弦调制模型,在原始Chebyshev窗的基础上引入了调制次数和调制范围,通过调节这两个参数可以人工可视化的调节截断误差,新的窗函数继承了Chebyshev窗的特点,在不明显降低截断谱范围的基础上明显降低了截断误差.本文针对不同正演阶数N,给出了一组经验调制系数,并通过数值模拟方法,对比了新方法、改进二项式窗和基于最小二乘优化方法的正演效果.结果表明,基于余弦调制的Chebyshev窗控制数值频散的能力更强,在大网格下可以得到更精确的正演结果.从经济角度分析,该方法减小了计算花费,提高了计算效率.The finite difference forward modeling is the basis of elastic wave reverse-time migration and full waveform inversion in the time domain.The accuracy of forward modeling also controls the accuracy of seismic imaging and inversion.The migration or inversion will bring illusion if the accuracy is not high.We can get optimized explicit finite difference operators by using the window function to truncate spatial convolution counterpart of the pseudo-spectral method.Based on this,a cosine modulated Chebyshev window is designed.On the basis of the original Chebyshev window,the modulation times and modulation domain are introduced,and we can adjust truncation error visually by controlling these two parameters.As the new window function inherits the character of Chebyshev window,we observe that the spectral range using the modulated window function for truncation is significantly broader than using the conventional window function with stable error.For different forward modeling orders N,we give a set of empirical modulation factors and compare the forward modeling effect of the new method andimproved binomial window by the numerical simulation method.The results demonstrate that the operators based on the cosine modulated Chebyshev window can efficiently suppress the numerical dispersion and get more accurate forward modeling results on the large grid.From economic perspective,this method reduces the computational cost and improves efficiency.

关 键 词：有限差分 数值频散 窗函数 弹性波 余弦调制 

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