基于Chebyshev自褶积组合窗的有限差分算子优化方法
详细信息 本馆镜像全文    |  推荐本文 | | 获取馆网全文
摘要
有限差分法广泛应用于地震波数值模拟、成像和波形反演中,差分数值解的精度直接影响着地震成像和反演的效果.因为有限差分算子可以通过截断伪谱法的空间褶积序列得到,而截断窗函数的属性影响有限差分算子逼近微分算子的精度.具体地讲,窗函数的幅值响应的主瓣和旁瓣决定了有限差分算子逼近的精度,主瓣越窄,旁瓣衰减越大,则有限差分算子逼近微分算子的精度越高,更好地压制数值频散.基于此认识,本文提出了一种基于Chebyshev自褶积组合窗截断逼近的有限差分算子优化方法.Chebyshev自褶积组合窗的主瓣较窄,且旁瓣衰减大,其可通过只调节三个参数,更直观和可视化地控制主瓣和旁瓣的形状,改变有限差分算子逼近微分算子的精度;该窗函数截断逼近的有限差分算子不仅有较大的谱覆盖范围,而且精度误差波动较小,这表明低阶的差分算子可以达到高阶算子的精度,且逼近误差更稳定;从经济上来讲,将有效地减少模拟计算花费,提高计算效率.
The finite-difference method has been widely utilized in seismic wave numerical modeling,seismic imaging and full waveform inversion.The accuracy of finite-difference numerical solutions directly affects results of seismic imaging and inversion.Using an autoconvolution combined window function to truncate spatial convolutional counterpart of the pseudospectral method,optimized explicit finite-difference operators are derived.The truncated window function method is used to get optimized finite-difference operators.Firstly,we analyze the influence on the accuracy of finite-difference operators caused by the properties of main lobe and side lobe in the amplitude response of truncated window functions.Secondly,based on the methods of auto-convolution and weighted combination,a window function which has narrower main lobe and larger attenuation of side lobe is designed,correspondingly bringing higher accuracy of finite-difference approximation.Finally,we use the window function to get optimized finite-difference operators.From the analysis of window functions,the factors that affect the accuracy of finitedifference approximation can be summarized in two aspects:(1)The width of main lobe of window functions.(2)The attenuation of side lobe of window functions.The window functions which have narrower main lobe and larger attenuation of side lobe can yield higher accuracy of finite-difference approximation.The Chebyshev window function has appropriate width of main lobe and attenuation of side lobe to get better stability of accuracy error on the premise of maintaining the appropriate wave number coverage range.Furthermore,the auto-convolution method will increase the attenuation of side lobe,however,widen the main lobe.Weighted combination can remedy the defect,choose different weight coefficients to reduce the width of main lobe.Combining the two methods,a specific window function named Chebyshev autoconvolution combined window is designed.Compared with the accuracy curves of approximation between the finite-difference operators truncated by Chebyshev auto-convolution combined window and the conventional operators,the former lead to great accuracy in a bigger frequency region.Tests on a homogeneous model and the Marmousi model show that the dispersion caused by the operators based on the Chebyshev auto-convolution combined window is quite weak under the same order of the conventional operators.Further comparison with the finite-difference operators truncated by improved binomial window,our operators still have a bigger frequency coverage range and smaller fluctuation of accuracy error.The finite-different operators based on the Chebyshev auto-convolution combined window have higher accuracy than that of the conventional operators and operators based on the improved binomial window.Our optimized eighth-order and twelfth-order operators can respectively reach,even exceed the accuracy of the conventional twelfth-order and twenty-fourth-order finitedifference operators,and the maximum deviation of absolute error is within [0,0.0004].For higher-order operators,the accuracy increase becomes more obvious.The results of elastic wavefield numerical modeling demonstrate that the operators based on the Chebyshev autoconvolution combined window can efficiently suppress the numerical dispersion and has greater modeling accuracy under the same discretizations without extra computing costs.In addition,through adjusting parameters of the auto-convolution combined window,we can adjust the accuracy of finite-difference approximation as required,visually and intuitively.
引文
Alterman Z,Karal F C Jr.1968.Propagation of elastic waves inlayered media by finite difference methods.Bulletin ofSeismological Society of America,58(1):367-398.
    Boris J P,Book D L.1973.Flux-corrected transport.I.SHASTA,a fluid transport algorithm that works.Journal of ComputationalPhysics,11(1):38-69.
    Carcione J M,Kosloff D,Behle A,et al.1992.A spectral schemefor wave propagation simulation in 3-D elastic-anisotropic media.Geophysics,57(12):1593-1607,doi:10.1190/1.1443227.
    Cheng B J,Li X F,Long G H.2008.Seismic waves modeling byconvolutional Forsyte polynomial differentiator method.ChineseJ.Geophys.(in Chinese),51(2):531-537.
    Chu C L,Stoffa P L.2012.Determination of finite-differenceweights using scaled binomial windows.Geophysics,77(3):W17-W26,doi:10.1190/GEO2011-0336.1.
    Dablain M A.1986.The application of high-order differencing to thescalar wave equation.Geophysics,51(1):54-66,doi:10.1190/1.1442040.
    Diniz P S R,da Silva E A B,Netto S L.2012.Digital SignalProcessing System Analysis and Design.Beijing:China MachinePress.
    Dong L G,Ma Z T,Cao J Z,et al.2000.A staggered-grid highorder difference method of one-order elastic wave equation.Chinese J.Geophys.(in Chinese),43(3):411-419.
    Fornberg B.1987.The pseudospectral method:Comparisons withfinite differences for the elastic wave equation.Geophysics,52(4):483-501,doi:10.1190/1.1442319.
    Gazdag J.1981.Modeling of the acoustic wave equation withtransform methods.Geophysics,46(6):854-859,doi:10.1190/1.1441223.
    Holberg O.1987.Computational aspects of the choice of operatorand sampling interval for numerical differentiation in large-scalesimulation of wave phenomena.Geophysical Prospecting,35(6):629-655,doi:10.1111/j.1365-2478.1987.tb00841.x.
    Igel H,Mora P,Riollet B.1995.Anisotropic wave propagationthrough finite-difference grids.Geophysics,60(4):1203-1216,doi:10.1190/1.1443849.
    Kelly K R,Ward R W,Treitel S,et al.1976.Syntheticseismograms:A finite-difference approach.Geophysics,41(1):2-27,doi:10.1190/1.1440605.
    Kosloff D D,Baysal E.1982.Forward modeling by a Fouriermethod.Geophysics,47(10):1402-1412,doi:10.1190/1.1441288.
    Lee C,Seo Y.2002.A new compact spectral scheme for turbulencesimulations.Journal of Computational Physics,183(2):438-469,doi:10.1006/jcph.2002.7201.
    Liu Y,Li C C,Mou Y G.1998.Finite-difference numericalmodeling of any even-order accuracy.OGP(in Chinese),33(1):1-10.
    Liu Y,Sen M K.2009a.A new time-space domain high-order finitedifferent method for the acoustic wave equation.Journal ofComputational Physics,228(23):8779-8806,doi:10.1016/j.jcp.2009.08.027.
    Liu Y,Sen M K.2009b.Numerical modeling of wave equation by atruncated high-order finite-difference method.EarthquakeScience,22(2):205-213,doi:10.1007/s11589-009-0205-0.
    Liu Y,Sen M K.2010.Acoustic VTI modeling with a time-space domain dispersion-relation-based finite-difference scheme.Geophysics,75(3):A11-A17,doi:10.1190/1.3374477.
    Liu Y.2013.Globally optimal finite-difference schemes based onleast squares.Geophysics,78(4):T113-T132,doi:10.1190/geo2012-0480.1.
    Pei Z L,Mou L G.2003.A staggered-grid high-order differencemethod for modeling seismic wave propagation ininhomogeneous media.Journal of China University ofPetroleum(Edition of Natural Science)(in Chinese),27(6):17-27.
    Robertsson J O A,Blanch J O,Symes WW,et al.1994.Galerkinwavelet modeling of wave propagation:Optimal finite-differencestencil design.Mathematical and Computer Modelling,19(1):31-38,doi:10.1016/0895-7177(94)90113-9.
    Yan H Y,Liu Y.2013a.Visco-acoustic prestack reverse-timemigration based on the time-space domain adaptive high-orderfinite-difference method.Geophysical Prospecting,61(5):941-954,doi:10.1111/1365-2478.12046.
    Yan H Y,Liu Y.2013b.Pre-stack reverse-time migration based onthe time-space domain adaptive high-order finite-differencemethod in acoustic VTI medium.Journal of Geophysics andEngineering,10(1):015010,doi:10.1088/1742-2132/10/1/015010.
    Yang D H,Teng J W.1997.FCT finite difference modeling ofthree-component seismic records in anisotropic medium.OGP(in Chinese),32(2):181-190.
    Yang L,Yan H Y,Liu H.2014.Least squares staggered-gridfinite-difference for elastic wave modelling.Exploration Geophysics,45(4):255-260,doi:10.1071/EG13087.
    Zhang J H,Yao Z X.2013.Optimized finite-difference operator forbroadband seismic wave modeling.Geophysics,78(1):A13-A18,doi:10.1190/GEO2012-0277.1.
    Zhou B,Greenhalgh S A.1992.Seismic scalar wave equationmodeling by a convolutional differentiator.Bulletin of theSeismological Society of America,82(1):289-303.
    程冰洁,李小凡,龙桂华.2008.基于广义正交多项式褶积微分算子的地震波场数值模拟方法.地球物理学报,51(2):531-537.
    董良国,马在田,曹景忠等.2000.一阶弹性波方程交错网格高阶差分解法.地球物理学报,43(3):411-419.
    刘洋,李承楚,牟永光.1998.任意偶数阶精度有限差分法数值模拟.石油地球物理勘探,33(1):1-10.
    裴正林,牟永光.2003.非均匀介质地震波传播交错网格高阶有限差分法模拟.石油大学学报(自然科学版),27(6):17-27.
    杨顶辉,腾吉文.1997.各向异性介质中三分量地震记录的FCT有限差分模拟.石油地球物理勘探,32(2):181-190.

版权所有:© 2023 中国地质图书馆 中国地质调查局地学文献中心