Characteristics of seismic wave field in frequency-space domain in strong attenuation media
ZHANG Yi1, WANG Yun1, CHEN Benchi2, WANG Xiang-chun1
1. School of Geophysics and Information Technology, China University of Geosciences(Beijing), Beijing 100083, China; 2. Sinopec Science & Technology Department, Beijing 100728, China
Abstract:As a viscoelastic medium model,the strong attenuation model applies Biot’s basic idea to modify the viscous dissipation coefficient,and describes the absorption and attenuation characteristics of the medium through the equation of motion.Compared with some conventional viscoeastic medium models and viscoelastic-porous medium models,the strong attenuation model can more easily and accurately describe the strong attenuation properties of media such as heavy oil reservoirs and loose sediments near the surface.In this paper,we implement a method of 25-point frequency-space domain finite difference to simulate the wave field characteristics of a strong attenuation model medium,and study the attenuation mechanism of seismic waves caused by three physical factors including porosity,fluid viscosity and viscoelasticity of medium.From the numerical simulation results,we can find that porosity,fluid viscosity and viscoelasticity are all important factors on inducing strong attenuation and high-velocity dispersion of seismic waves,and their influence on attenuation of S-wave is more obvious than that of P-wave.Of which,viscoelasticity is the key factor on affecting the attenuation of high-frequency component of seismic waves,and both porosity and fluid viscosity can cause the attenuation of energy in the effective frequency range of seismic waves,especially porosity has the strongest attenuation effect.We also study the influence of the shallow medium with strong attenuation on the deep wave field,and find that the strong attenuation model is more practical in describing the low-velocity medium by comparing the wave field simulation results under the complete elasticity theory and general viscoelasticity.These results provide a reference to the study of strong attenuation media involving different attenuation mechanisms,and also lay a foundation for the establishment of the comprehensive strong attenuation compensation theory.
Carcione J M.Seismic modeling in viscoelastic media[J].Geophysics,1933,58(1):110-120.
[2]
牛滨华,孙春岩.地震波理论研究进展——介质模型与地震波传播[J].地球物理学进展,2004,19(2):255-263.NIU Binhua,SUN Chunyan.Developing theory of propagation of seismic waves:medium model and propagation of seismic waves[J].Progress in Geophysics,2004,19(2):255-263.
[3]
奚先,姚姚.二维黏弹性随机介质中的波场特征[J].石油地球物理勘探,2004,39(4):381-387.XI Xian,YAO Yao.Characteristics of wavefield in 2-D viscoeiastic random medium[J].Oil Geophysical Prospecting,2004,39(4):381-387.
[4]
李合群,孟小红,赵波,等.塔里木沙漠区地震数据品质与沙层Q吸收[J].石油地球物理勘探,2010,45(1):28-34.LI Hequn,MENG Xiaohong,ZHAO Bo,et al.Seismic data quality and sand layer Q absorption in Tarim desert area[J].Oil Geophysical Prospecting,2010,45(1):28-34.
[5]
宋吉杰,禹金营,王成,等.近地表介质Q估计及其在塔河北部油田的应用[J].石油物探,2018,57(3):436-442.SONG Jijie,YU Jinying,WANG Cheng,et al.Q estimation for near-surface media and its application in the Northern Tahe Oilfield,China[J].Geophysical Prospecting for Petroleum,2018,57(3):436-442.
[6]
牛滨华,孙春岩.黏弹性介质与地震波传播[M].北京:地质出版社,2007.
[7]
Stokes G G.On the aberration of light[J].The London,Edinburgh,and Dublin Philosophical Magazine and Journal of Science,1845,27(177):9-15.
[8]
孙成禹.地震波理论与方法[M].山东东营:中国石油大学出版社,2007.
[9]
Carcione J M,Kosloff D,Kosloff R.Wave propagation simulation in a linear viscoacoustic medium[J].Geophysical Journal International,1988,93(2):393-401.
[10]
Carcione J M,Kosloff D,Kosloff R.Wave propagation simulation in a linear viscoacoustic medium[J].Geophysical Journal International,1988,95(3):597-611.
[11]
张金波,杨顶辉,贺茜君,等.求解双相和黏弹性介质波传播方程的间断有限元方法及其波场模拟[J].地球物理学报,2018,61(3):926-937.ZHANG Jinbo,YANG Dinghui,HE Qianjun,et al.Discontinuous Galerkin method for solving wave equations in two-phase and viscoelastic media[J].Chinese Journal of Geophysics,2018,61(3):926-937.
[12]
刘财,胡宁,郭智奇,等.基于分数阶时间导数常Q黏弹本构关系的含黏滞流体双相VTI介质中波场数值模拟[J].地球物理学报,2018,61(6):2446-2458.LIU Cai,HU Ning,GUO Zhiqi,et al.Numerical simulation of the wavefield in a viscous fluid-saturated two-phase VTI medium based on the constant-Q viscoela-stic constitutive relation with a fractional temporal derivative[J].Chinese Journal of Geophysics,2018,61(6):2446-2458.
[13]
Biot M A.Theory of propagation of elastic waves in a fluid-saturated porous solid,I:Low-frequency range[J].The Journal of the Acoustical Society of America,1956,28(2):168-178.
[14]
Mavko G,Nur A.Melt squirt in the asthenosphere[J].Journal of Geophysical Research,1975,80(11):1444-1448.
[15]
Dvorkin J,Nur A.Dynamic poroelasticity:A unified model with the squirt and the Biot mechanisms[J].Geophysics,1993,58(4):524-533.
[16]
杜启振,刘莲莲,孙晶波.各向异性粘弹性孔隙介质伪谱法波场正演模拟[J].油气地球物理,2007,5(1):22-26.DU Qizhen,LIU Lianlian,SUN Jingbo.Wavefield numerical modeling with the pseudo-spectral method in anisotropic viscoelastic porous media[J].Petroleum Geophysics,2007,5(1):22-26.
[17]
Nie J X,Ba J,Yang D H,et al.BISQ model based on a Kelvin-Voigt viscoelastic frame in a partially satura-ted porous medium[J].Applied Geophysics,2012,9(2):213-222.
[18]
凌云,杜向东,曹思远.基于Zener线性体的黏弹孔隙介质衰减频散特征分析[J].地球物理学进展,2017,32(1):205-209.LING Yun,DU Xiangdong,CAO Siyuan.Attenuation and dispersion characteristics analysis in visco-poroelastic medium based on Zener linear solid[J].Progress in Geophysics,2017,32(1):205-209.
[19]
谢佩瑜,杨顶辉.近地表强衰减介质中的地震波传播模型[J].地球物理学报,2018,61(3):917-925.XIE Peiyu,YANG Dinghui.Seismic wave propagation model in near-surface strong attehuation media[J].Chinese Journal of Geophysics,2018,61(3):917-925.
[20]
Lysmer J,Lawrence A D.A finite element method for seismology[J].Methods in Computational Physics,1972,11:181-216.
[21]
Min D J,Shin C,Changsoo S,et al.Improved frequency-domain elastic wave modeling using weighted-ave-raging difference operators[J].Geophysics,2000,65(3):884-895.
[22]
李桂花,冯建国,朱光明.黏弹性VTI介质频率空间域准P波正演模拟[J].地球物理学报,2011,54(1):200-207.LI Guihua,FENG Jianguo,ZHU Guangming.Quasi-P wave forward modeling in viscoelastic VTI media in frequency-space domain[J].Chinese Journal of Geophysics,2011,54(1):200-207.
[23]
Min D J,Shin C,Yoo H S.Free-surface boundary condition in finite-difference elastic wave modeling[J].Bulletin of the Seismological Society of America,2004,94(1):237-250.
[24]
刘财,杨庆节,鹿琪,等.双相介质中地震波的频率-空间域数值模拟[J].地球物理学报,2014,57(9):2885-2899.LIU Cai,YANG Qingjie,LU Qi,et al.Simulation of wave propagation in two-phase porous media using a frequency-space domain method[J].Chinese Journal of Geophysics,2014,57(9):2885-2899.
[25]
周聪,刘江平,罗银河,等.二维频率域全波场有限差分数值模拟方法[J].石油地球物理勘探,2014,49(2):278-287.ZHOU Cong,LIU Jiangping,LUO Yinhe,et al.2D full-wavefield modeling in frequency domain using finite-difference[J].Oil Geophysical Prospecting,2014,49(2):278-287.
[26]
Yang L,Yang D H,Nie J X.Wave dispersion and a-ttenuation in viscoelastic isotropic media containing multiphase flow and its application[J].Science China:Physics,Mechanics & Astronomy,2014,57(6):1068-1077.
[27]
Lines L R,Treitel S.A review of least-squares inversion and its application to geophysical problems[J].Geophysical Prospecting,1984,32(2):159-186.
[28]
杨庆节,刘财,郭智奇,等.基于BISQ模型双相各向同性介质弹性波传播的频率-空间域有限差分模拟[J].地球物理进展,2015,30(1):249-260.YANG Qingjie,LIU Cai,GUO Zhiqi,et al.Wave propagation in two-phase isotropic medium based on BISQ model in frequency-space domain[J].Progress in Geophysics,2015,30(1):249-260.
[29]
Collino F,Tsogka C.Application of the perfectly matched absorbing layer model to the linear elastodynamic problem in anisotropic heterogeneous media[J].Geophysics,2001,66(1):294-307.
[30]
罗玉钦,刘财.多轴复频移近似完全匹配层弹性波模拟[J].石油地球物理勘探,2019,54(5):1024-1033.LUO Yuqin,LIU Cai.Multi-axial compex-frequency shifting nearly perfectly matched layer for seismic forward modeling in elastic media[J].Oil Geophysical Prospecting,2019,54(5):1024-1033.
[31]
Berenger J P.A perfectly matched layer for the ab-sorption of electromagnetic waves[J].Journal of Computational Physics,1994,114(2):185-200.
[32]
高凤霞.频率域波动方程多参数全波形反演方法研究[D].吉林长春:吉林大学,2014.
[33]
殷文,印兴耀,吴国忱,等.高精度频率域弹性波方程有限差分方法及波场模拟[J].地球物理学报,2006,49(2):561-568.YIN Wen,YIN Xingyao,WU Guochen,et al.The method of finite difference of high precision elastic wave equations in the frequency domain and wavefield simulation[J].Chinese Journal of Geophysics,2006,49(2):561-568.
[34]
Cui Q H,Rui Y J,Shang X M.Near-surface absorption compensation technique and its application[C].Near Surface Geophysics Asia Pacific Conference,2013,517-520.
[35]
邹冠贵.孔隙介质地震波传播及衰减特征评价研究[D].北京:中国矿业大学,2010.ZOU Guangui.Propagation of Elastic Wave in a Fluid-saturated Porous Media and Attenuation Characteristic Evaluation[D].China University of Mining,Beijing,2010.
[36]
赵秋芳,云美厚,朱丽波,等.近地表Q值测试方法研究进展与展望[J].石油地球物理勘探,2019,54(6):1397-1418.ZHAO Qiufang,YUN Meihou,ZHU Libo,et al.Progress and outlook of near-surface quality factor Q measurement and inversion[J].Oil Geophysical Prospecting,2019,54(6):1397-1418.