复杂介质地震波场模拟的单程波方法
|
摘要
单程波方法是常用的地震正演模拟方法之一。与双程波方法相比,单程波方法的计算量成倍减少,且能保持较高的精度。但由于单程波方法不能获得准确的入射角信息,因而多采用垂直入射的反射系数,使得模拟波场的振幅特征不是实际地层的真实反映。为此,提出了结合频率-波数域反射系数近似公式和相移加插值的改进单程波方法。分别采用双程波虚谱法、常规单程波方法和改进单程波方法,基于Marmousi模型进行了数值模拟,对比模拟结果发现,改进的单程波方法可以准确地模拟振幅随入射角的变化,且具备与双程波方法相同的精度。
The one-way wave method is commonly used for simulating seismic wave,compares with the two-way wave simulation method, the computation of one-way method becomes time of reduction,also maintains higher precision,but the one-way wave method is unable to obtain the accurate incidence angle information,thus the normal incidence reflectivity is used,but the amplitude characteristics of simulated wave field is not the real reflection of actual formations. Therefore,we proposed an improved one-way wave method by combining reflectivity approximate formula in frequency-wave number domain and phase-shift plus interpolation method.Two-way wave pseudo spectral method,conventional one-way wave method and the improved one-way wave method were separately adopted to simulate single shot recorder based on Marmousi model.Comparison on simulated results show that the improved one-way wave method can accurately simulate the changes of amplitude via incident angle and has the same precision with two-way wave method.
The one-way wave method is commonly used for simulating seismic wave,compares with the two-way wave simulation method, the computation of one-way method becomes time of reduction,also maintains higher precision,but the one-way wave method is unable to obtain the accurate incidence angle information,thus the normal incidence reflectivity is used,but the amplitude characteristics of simulated wave field is not the real reflection of actual formations. Therefore,we proposed an improved one-way wave method by combining reflectivity approximate formula in frequency-wave number domain and phase-shift plus interpolation method.Two-way wave pseudo spectral method,conventional one-way wave method and the improved one-way wave method were separately adopted to simulate single shot recorder based on Marmousi model.Comparison on simulated results show that the improved one-way wave method can accurately simulate the changes of amplitude via incident angle and has the same precision with two-way wave method.
引文
1 Cerveny V,Molotkov I A,Psencik I.Ray method in seismology[M].Praha,Czech:Karlova University, 1977.63-72
2 Chanpman C H,Drummond R.Body wave seismograms in inhomogeneous media using Maslov asymptotic theory[J].Bulletin of the Seismological Society of America,1982,72(6B):S277- S317
3 Cerveny V.Synthetic body wave seismograms for laterally varying layered structures by the Gaussian beam method[J].Geophysical Journal International,1983, 73(2):389-426
4 Boore D M.Finite difference methods for seismic waves:methods in computational physics[M].New York,USA:Academic Press,1972.21-22
5 Keny K R,Ward R W,Treitel S,et al.Synthetic seismograms: a finite-difference approach[J].Geophysics, 1976,41(2):2-27
6 Smith W D.A nonreflecting plane boundary for wave propagation problem[J].Journal of Computational Physics,1974,15(4):492-503
7 Wong H L,Jennings P C.Effect of canyon topography on strong ground motion[J].Bulletin of the Seismological Society of America,1975,65(5):1 239-1 257
8 Orzag S A.Spectral methods for problems in complex geometries[J].Journal of Computational Physics, 1980,37(1):70-92
9熊高军,张琳.共炮记录正演模拟检波点下延记录原理[J].石油物探,1999,38(2):43-49
10熊小兵,贺振华.相位移加有限差分法波动方程正演模拟[J].石油物探,1998,37(3):22-28
11熊晓军,贺振华,黄德济.三维波动方程正演及模型应用研究[J].石油物探,2005,44(6):554-556
12刘礼农,崔凤林,张剑峰.三维复杂构造中地震波模拟的单程波方法[J].地球物理学报,2004,47(3):514- 520
13吴永国,贺振华,黄德济.串珠状溶洞模型介质波动方程正演与偏移[J].地球物理学进展,2008,23(2): 539-544
14 Gazdag J,Sguazzero P.Migration of seismic data by phase shift plus interpolation[J].Geophysics,1984,49 (2):124-131
15 de Bruin C G M,Wapenaar C P A,Berkhout A J.Angle -dependent reflectivity by shot record migration[J]. Geophysics,1990,55(9):1 093-1 096
16 Wapenaar C P A,Berkhout A J.Elastic wave field extrapolation [M].Amsterdam,NetherLands:Elsevier Science,1989.102-105
2 Chanpman C H,Drummond R.Body wave seismograms in inhomogeneous media using Maslov asymptotic theory[J].Bulletin of the Seismological Society of America,1982,72(6B):S277- S317
3 Cerveny V.Synthetic body wave seismograms for laterally varying layered structures by the Gaussian beam method[J].Geophysical Journal International,1983, 73(2):389-426
4 Boore D M.Finite difference methods for seismic waves:methods in computational physics[M].New York,USA:Academic Press,1972.21-22
5 Keny K R,Ward R W,Treitel S,et al.Synthetic seismograms: a finite-difference approach[J].Geophysics, 1976,41(2):2-27
6 Smith W D.A nonreflecting plane boundary for wave propagation problem[J].Journal of Computational Physics,1974,15(4):492-503
7 Wong H L,Jennings P C.Effect of canyon topography on strong ground motion[J].Bulletin of the Seismological Society of America,1975,65(5):1 239-1 257
8 Orzag S A.Spectral methods for problems in complex geometries[J].Journal of Computational Physics, 1980,37(1):70-92
9熊高军,张琳.共炮记录正演模拟检波点下延记录原理[J].石油物探,1999,38(2):43-49
10熊小兵,贺振华.相位移加有限差分法波动方程正演模拟[J].石油物探,1998,37(3):22-28
11熊晓军,贺振华,黄德济.三维波动方程正演及模型应用研究[J].石油物探,2005,44(6):554-556
12刘礼农,崔凤林,张剑峰.三维复杂构造中地震波模拟的单程波方法[J].地球物理学报,2004,47(3):514- 520
13吴永国,贺振华,黄德济.串珠状溶洞模型介质波动方程正演与偏移[J].地球物理学进展,2008,23(2): 539-544
14 Gazdag J,Sguazzero P.Migration of seismic data by phase shift plus interpolation[J].Geophysics,1984,49 (2):124-131
15 de Bruin C G M,Wapenaar C P A,Berkhout A J.Angle -dependent reflectivity by shot record migration[J]. Geophysics,1990,55(9):1 093-1 096
16 Wapenaar C P A,Berkhout A J.Elastic wave field extrapolation [M].Amsterdam,NetherLands:Elsevier Science,1989.102-105
版权所有：© 2023 中国地质图书馆 中国地质调查局地学文献中心