| 413 | 0 | 352 |
| 下载次数 | 被引频次 | 阅读次数 |
地下介质往往表现为粘弹特性,研究基于粘弹假设的数值模拟方法对于正确认识地震波的传播规律和提高地震勘探精度具有重要意义。波动方程正演过程中纵横波的保幅解耦是研究准确的粘弹介质中地震波传播机理的前提,基于散度和旋度算子的纵横波解耦方法会使波场的相位和振幅产生畸变,且解耦后的波场在极性反转位置上无法与分离前混合波场各分量对应。在散度和旋度算子上再做一次梯度和旋度运算的波场分离方法虽然能够克服上述缺陷,但存在保幅性差等问题。本文从粘弹介质中的一阶速度-应力方程出发,推导了矢量纵横波分离的波数域表达式,结合有限差分思路给出了其在空间域的求解方法。本文方法利用纵横波的传播速度对现有的矢量波场分离方法进行振幅校正,并将校正结果分别作为纵波与横波对时间的二阶偏导,实现了粘弹介质中的纵横波分离。模型试算结果表明,本文方法能够克服现有方法的缺陷,获得更具保幅性的波场分离结果。
Abstract：Due to the viscoelastic properties of underground rocks, it is important to study the numerical simulation method based on viscoelastic assumption for correctly understanding the propagation mechanism of seismic waves and improving seismic exploration accuracy. Amplitude-preserving decoupling of P-wave and S-wave is the prerequisite for accurately understand propagation mechanism of seismic wave in viscoelastic medium. The divergence and curl operators may cause the phase and amplitude distortion of P-wave and S-wave, moreover the physical meaning and polarity inversion position of decoupled wave field is not inconsistent with the original mixed wave field. Another wave field separation method which separate the P-wave and S-wave with gradient and curl operations on divergence and curl operators can overcome the above shortcomings, but the algorithm has the disadvantage of poor amplitude preservation. In this paper, a wavenumber domain expression of vectorial P-and S-wave separation formulas is derived which was based on one-order velocity-stress viscoelastic equations, and its solution method in spatial domain is also according to the finite difference method. In our method, the velocity of P-and S-waves are used to correct the amplitude of the existing wave field separation method, and the correction results are used as the second-order partial derivatives of P-wave and S-wave for time, respectively. Based on above methods, the P-and S-wave separation in viscoelastic media is realized. The experimental results show that the proposed method can overcome the shortcomings of the existing methods and obtain more amplitude-preserving wave field separation results.
[1] Stockes G G.On the theories of internal friction of fluids in motion and of the equilibrium and motion of elastic solids[J].Trans Cambridge Philos Soc,1845,8(Part Ⅲ):287-319.
[2] Emmerich H,Korn M.Incorporation of attenuation into time-domain computations of seismic wave fields[J].Geophysics,1987,52(9):1252-1264.
[3] Carcione J M.Viscoacoustic wave propagation simulation in the earth[J].Geophysics,1988,53(6):769-777.
[4] Carcione J M,Dan K,Kosloff R.Wave propagation simulation in a linear viscoacoustic medium[J].Geophysical Journal International,1988,93(2):393-401.
[5] Robertsson J O A,Blanch J O,Levander A,et al.3-D viscoelastic finite-difference modeling[J].Geophysics,1994,59(9):1444-1456.
[6] 刘财,张智,邵志刚,等.线性粘弹体中地震波场伪谱法模拟技术[J].地球物理学进展,2005,20(3):640-644.Liu C,Zhang Z,Shao Z G.Pseudo spectral forward modeling of seismic wave in linear viscoelasic solid[J].Progress in Geophysics,2005,20(3):640-644.
[7] 孙成禹,印兴耀.三参数常Q 粘弹性模型构造方法研究[J].地震学报,2007,29(4):348-357.Sun Y C,Yin X Y.Construction of constant Q viscoelastic model with three parameters[J].Acta Seismologica Sinica,2007,29(4):348-357.
[8] 孟凡顺,郭海燕,等.复杂地质体粘滞弹性波正演模拟的有限差分法[J].青岛海洋大学学报(自然科学版),2000,30(2):315-320.Meng F S,Guo H Y.Viscoelastic wave simulating in complex medium by finite difference method[J].Journal of Ocean University of Qingdao,2000,30(2):315-320.
[9] 王美霞.双相及粘弹性介质中波传播方程的保辛算法及其波场模拟[D].北京:清华大学,2012.Wang M X.Symplectic Stereo Modelling Method for Wave Equations in Two-Phase and Viscoelastic Media and Its Numerical Simulations[D].Beijing:Tsinghua University,2012.
[10] 张明,王润秋.有限元解三维弹性波方程的并行算法[J].计算数学,1995,2:127-135.Zhang M,Wang R Q.a parallel algorithm with finite element for solving 3-D elastic wave equations[J].Mathematica Numerica Sinica,1995,2:127-135.
[11] 赵童鑫.粘弹性介质中的声波方程及其有限差分数值模拟[D].成都:西南石油大学,2016.Zhao T X.Acoustic Wave Equation in Viscoelastic Media and Finite Difference Numerical Simulation[D].Chengdu:Southwest Petroleum University,2016.
[12] 董良国,马在田,曹景忠.一阶弹性波方程交错网格高阶差分解法稳定性研究[J].地球物理学报,2000,439(6):856-864.Ma L G,Ma Z T,Cao J Z.A study on stability of the staggered grid high order difference method of first order elastic wave equation[J].Chinese Journal of Geophysics,2000:43(6):856-864.
[13] Virieux J.P-SV wave propagation in heterogeneous media:Velocity-stress finite-difference method[J].Geophysics,1986,51(4):889-901.
[14] 胡楠,何兵寿.三维各向同性介质矢量波场保幅分离方法[J].煤炭学报,2017,42(9):2420-2426.Hu L,He B S.Amplitude-preserving separation method of 3D isotropic medium vector wave field[J].Journal of the China Coal Society,2017,42(9):2420-2426.
[15] Berenger J P.A perfectly matched layer for the absorption of electromagnetic waves[J].Journal of Computational Physics,1994,114(2):185-200.
[16] Weng C C,Weedon W H.A 3D perfectly matched medium from modified maxwell's equations with stretched coordinates[J].Microwave & Optical Technology Letters,1994,7(13):599-604.
[17] Aki K,Richards P G.Quantitative Seismology[M].2nd Ed.[s.l.]:W.H.Freeman,1980:5-43.
[18] Sun R,Chow J,Chen K J.Phase correction in separating P-and S-waves in elastic data[J].Geophysics,2001,66(5):1515-1518.
[19] Sun R,Mcmechan G A,Chuang H H.Amplitude balancing in separating P-and S-waves in 2D and 3D elastic seismic data[J].Geophysics,2011,76(3):S103.
[20] 李振春,雍鹏,黄建平,等.基于矢量波场分离弹性波逆时偏移成像[J].中国石油大学学报(自然科学版),2016,40(1):42-48.Li Z C,Yong P,Huang J P.Elastic wave reverse time migration based on vector wave field separation[J].Journal of China University of Petroleum(Edition of Natural Science),2016,40(1):42-48.
[21] 王鹏飞,何兵寿.基于行波分离的三维弹性波矢量场点积互相关成像条件[J].石油地球物理勘探,2017,52(3):477-483.Wang P F,He B S.Vector field dot product cross-correlation imaging based on 3D elastic wave separation[J].OGP,2017,52(3):477-483.
[22] 何兵寿,张会星.多分量波场的矢量法叠前深度偏移技术[J].石油地球物理勘探,2006,41(4):368-374.He B S,Zhang H X.Vector prestack depth migration of multi-component wavefield[J].OGP,2006,41(4):368-374.
[23] Chang W F.Elastic reverse-time migration[J].Geophysics,1987,52(10):1365-1375.
[24] Sun R,McMechan G,Hsiao H,et al.Separating P-and S-waves in prestack 3D elastic seismograms using divergence and curl[J].Geophysics,2004,69(1),286-297.
[25] Tang H G,He B S,Mou H B.P-and S-wave energy flux density vectors[J].Geophysics,2016,81(6):T357-T368.
[26] Chen T,He B S.A normalized wave field separation cross-correlation imaging condition for reverse-time migration based on Poynting Vector[J].Applied Geophysics,2014,11(2):158-188.
基本信息:
DOI：10.16441/j.cnki.hdxb.20180397
中图分类号:P631.4
引用信息:
[1]侯志强,尹文笋,李键,等.基于纵横波保幅分离的粘弹介质弹性波正演模拟[J],2020,50(01):82-92.DOI:10.16441/j.cnki.hdxb.20180397.
基金信息:
国家重大科技专项项目(2016ZX05027-002)资助~~