TTI介质qP波方程频率—空间域加权平均有限差分算子
详细信息 本馆镜像全文    |  推荐本文 | | 获取馆网全文
摘要
波动方程有限差分方法能够较精确地模拟任意非均匀介质中的地震波场,但它本身存在着数值频散问题。在具有倾斜对称轴的横向各向同性介质(TTI介质)地震波正演模拟中,为了解决常规有限差分算子的数值频散问题,本文构造了频率—空间域qP波方程加权平均有限差分算子,求取了归一化相速度,并根据最优化理论中的高斯牛顿法确定了加权平均差分算子的最优加权系数。利用常规差分算子和加权平均差分算子对归一化相速度进行了频散分析,并对均匀TTI介质(包括各向同性介质和椭圆各向异性介质)中的qP波地震波场进行了有限差分数值模拟。结果表明:加权平均有限差分算子具有较高的数值精度,能有效地压制常规有限差分算子的数值频散,为TTI介质频率—空间域qP波正演模拟奠定了基础。
Wave-equation finite-difference algorithm can more preciously simulate seismic wavefield for any non-uniform medium,but have the issue of numeric dispersion. In a seismic wave forward simulation in titled transversely isotropic medium (TTI medium) with titled symmetric axis, in order to solve the issue of numeric dispersion of ordinary finite-difference operator, the paper constructed the weighted mean finite-difference operator of qP-wave equation in frequency-space domain,computed normalized phase velocity and determined the optimal weighted coefficient of weighted mean difference operator according to Gauss-Newton approach of optimization theory. Using ordinary difference operator and weighted mean difference operator to analyze the dispersion of normalized phase velocity and carry out numeric finite-difference simulation of qP-wave seismic wavefield in u-niform TTI medium (including isotropic medium and ellipsoid anisotropic medium). The simulated results showed that the weighted mean finite-difference operator is characterized by higher numeric precision and capable to effectively suppress the numeric dispersion by ordinary finite-difference operator,laying the foundation of qP-wave forward simulation of TTI medium in frequency-space domain.
引文
[1]牟永光,裴正林.三维复杂介质地震数值模拟.北京:石油工业出版社,2005
    [2]Jo C H,Shin C and Suh J H.An opti mal9-point,fi-nite-difference,frequency-space,2-D scalar wave ex-trapolator.Geophysics,1996,61(2):529~537
    [3]Shin C and Sohn H.Afrequency-space2-Dscalar wave extrapolator using extend25-point finite-differ-ence operator.Geophysics,1998,63(1):289~296
    [4]ˇStekl I and Pratt R G.Accurate viscoelastic modeling by frequency-domain finite differences using rotated operators.Geophysics,1998,63(5):1779~1794
    [5]Min DJ,Shin C,Kwon B Det al.I mproved frequen-cy-domain elastic wave modeling using weighted-av-eraging difference operators.Geophysics,2000,65(3):884~895
    [6]吴国忱,梁锴.VTI介质频率—空间域准P波正演模拟.石油地球物理勘探,2005,40(5):535~545
    [7]Min DJ,Yoo HS,Shin C,Hyun HJ et al.Weigh-ted-averaging finite-element method for scalar wave e-quation in the frequency domain.J Seism Explor,2002,11:197~222.
    [8]Min D J,Shin C,Pratt R G,et al.Weighted-avera-ging finite-element method for2D elastic wave equa-tions in the frequency domain.BSSA,2003,93(2):904~921.
    [9]牛滨华,何樵登,孙春岩.六方各向异性介质方位矢量波动方程及其相速度.石油物探,1994,33(1):19~29
    [10]Winterstein D F.Velocity anisotropy terminology for geophysicists.Geophysics,55(8):1070~1088
    [11]吴国忱.各向异性介质地震波传播与成像.山东东营:中国石油大学出版社,2005
    [12]Alkhalifah T.Acoustic approxi mation for processing in transversely isotropic media.Geophysics,1998,63(2):623~631
    [13]Alkhalifah T.An acoustic wave equation for aniso-tropic media.Geophysics,2000,65(4):1239~1250
    [14]Zhang L B,Rector J WⅢ and Hoversten G M.An a-coustic wave equation for modeling intilted TI media.Expanded Abstracts of 73rd SEG Mtg,2003,153~156.
    [15]Thomsen L.Weak elastic anisotropy.Geophysics,1986,51(10):1954~1966
    [16]吴国忱.波场模拟中的数值频散分析与校正策略.地球物理学进展,2005,20(1):58~65

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