二维频率空间域粘声波正演模拟研究
|
摘要
地球介质本身的非弹性导致大地吸收效应,黏滞性影响波场的所有频率成分且对高频成分的影响更大,从而降低了地震记录和地震成像的垂向分辨率.通常使用经验公式描述吸收衰减.在时间域模拟吸收衰减需大量的计算机内存、计算效率较低,模拟Q值非常困难;而频率空间域可以通过引入复速度和Q值,使用任意的经验公式而进行粘滞介质的数值模拟.本文采用非常节省计算机内存需求的嵌套剖分网格排序法,扩展频率空间域声波多炮模拟快速计算流程到粘声波介质而方便高效地描述吸收衰减效应.数值试验证明频率空间域很容易模拟吸收与衰减效应,为后续的保振幅处理提供了很好的基础.
The earth is anelastic which leads to earth absorption phenomena,viscosity affects all frequencies component of the wavefield,the higher frequency component,the larger it affects,decrease the seismic recording and imaging's vertical resolution.Usually we use empirical formulation to describe absorption and attenuation.Modeling absorption and attenuation need huge computer's storage in time-space domain,low computational efficiency,it is very difficult in modeling absorption and attenuation.While we can use any empirical formulation in frequency-space domain to model absorption and attenuation numerical simulation,through introducing complex velocity and attenuation factor,adopt the very saving computer's storage's nested dissection,extend the multishot modeling fast computational procedure to viscosity,easily describe absorption and attenuation with a high efficiency.Numerical simulations on Marmousi demonstrate that it is very easy to model absorption and attenuation effect in frequencyspace domain,providing good tools for the subsequent amplitude preservation treatment.
The earth is anelastic which leads to earth absorption phenomena,viscosity affects all frequencies component of the wavefield,the higher frequency component,the larger it affects,decrease the seismic recording and imaging's vertical resolution.Usually we use empirical formulation to describe absorption and attenuation.Modeling absorption and attenuation need huge computer's storage in time-space domain,low computational efficiency,it is very difficult in modeling absorption and attenuation.While we can use any empirical formulation in frequency-space domain to model absorption and attenuation numerical simulation,through introducing complex velocity and attenuation factor,adopt the very saving computer's storage's nested dissection,extend the multishot modeling fast computational procedure to viscosity,easily describe absorption and attenuation with a high efficiency.Numerical simulations on Marmousi demonstrate that it is very easy to model absorption and attenuation effect in frequencyspace domain,providing good tools for the subsequent amplitude preservation treatment.
引文
[1]Parra J O.The transversely isotropic poroelastic wave equation including the Biot and the squirt mechanisms theory and application[J]:Geophysics,62:309-318.1997.
[2]Pride,S.R.,and J.G.Berryman,Linear dynamics of double-porosity dual-permeability materials.I.Governing equations and acoustic attenuation:Physical Review E 68, 036603(1-10).2003.
[3]张立彬.地震波非弹性吸收与反Q补偿稳定性研究[博士论文].上海:同济大学海洋地质国家重点实验室,2009. Zhang Li Bin.Research on The Seismic Anelatic Absorption and The Stability of Inverse Q Compensation[D].Shanghai: State Key Laboratory of Marine Geology,Tongji University, 2009.
[4]马昭军,刘洋.地震波衰减反演研究综述[J].地球物理学进展, 2005,20(4):1074-1082. MA Zhao-jun,LIU Yang.A summary of research on seismic attenuation[J].Progress in Geophysics(in Chinese),2005, 20(4):1074-1082.
[5]徐平,王宝善,张尉,林建民,陈顒.利用互相关函数求地震波衰减.地球物理学报,2006,49(6):1738-1744. XU Ping,WANG Bao-Shan,ZHANG Wei,LIN Jian-Min, CHEN Yong.Estimating seismic attenuation using crosscorrelation function.Chinese J.Geophys.(in Chinese), 2006,49(6):1738-1744.
[6]李振春,王清振.地震波衰减机理及能量补偿研究综述[J].地球物理学进展,2007,22(4):1147-1152. LI Zhen-chun,WANG Qing-zhen.A review of research on mechanism of seismic attenuation and energy compensation [J],Progress in Geophysics(in Chinese),2007,22(4): 1147-1152.
[7]李子顺.地震波衰减规律及其恢复方法[J].地球物理学进展, 2007,22(5):1545-1551. LI Zi-shun.Seismic wave attenuation and recovering method [J].Progress in Geophysics(in Chinese),2007,22(5): 1545 - 1551.
[8]陈文超,高静怀.基于改进的最佳匹配地震子波的地震资料衰减特性分析.地球物理学报,2007,50(3):837-843. CHEN Wen-Chao,GAO Jing-Huai.Characteristic analysis of seismic attenuation using MBMSW wavelets.Chinese J. Geophys.(in Chinese),2007,50(3):837-843.
[9]魏修成,卢明辉,巴晶,杨慧珠.含黏滞流体各向异性孔隙介质中弹性波的频散和衰减[J].地球物理学报,2008,51(1): 213-220. Wei X C,Lu M H,Ba J,Yang H Z.Dispersion and attenuation of elastic waves in a viscous fluid-saturated anisotropic porous solid[J].Chinese J.Geophys.(in Chinese),2008,50(1): 213-220.
[10]Versteeg,R.,The Marmousi experience:Velocity model determination on a synthetic complex data set.The Leading Edge,13:927-936,1994.
[11]Bjorn Ursin,Tommy Toverud,Comparison of seismic dispersion and attenuation models.CWP408-References, 2002.
[12]Churl-Hyun Jo.An optimal 9-point,finite-difference frequencyspace 2-D scalar wave extrapolator,Geophysics,61(2): 529-537,1996.
[13]George A,Liu J W.Computer Solution of Large Sparse Positive Definite Systems,Prentice Hall,New York,1981.
[14]任浩然,王华忠,龚婷.标量地震波频率-空间域有限差分法数值模拟.石油物探,2009,48(1):20-26. Ren Hao ran,Wang Hua zhong,Gong ting.Scalar seismic wave finite-difference numerical simulation in frequencyspace domain.GEOPHYSICAL PROSPECTING FOR PETROLEUM,2009,48(1):20-26.
[15]廖建平.频率空间域地震波数值模拟方法研究[博士论文]. 上海:同济大学,2009. Liao Jianping.Research on Seismic Wave Numerical Modeling Method in Frequency-space Domain[D]. Shanghai:Tongji University,2009.
[16]廖建平,刘和秀,王华忠,等.快速高精度的频率空间域声波数值模拟方法研究[J].地球物理学进展.2011,26(4):1359 -1363. Liao J P,Liu H X,Wang H,et al.Study on rapid highly accurate acoustic wave numerical simulation in frequency space domain[J].Progress in Geophys.(in Chinese),2011,26 (4):1359-1363.
[17]Slekl and R.Gerhard Pratt.Accurate viscoelastic modeling by frequency-domain finite differences using rotated operators. Geophysics,63(5):1779-1794,1998.
[18]李庆忠.走向精确勘探的道路—高分辨率地震勘探系统工程剖析[M].北京:石油工业出版社,1994. Li Qing-Zhong.The Way to Obtain a Better Resolution in Seismic Prospecting—A Systematical Analysis of High Resolution Seismic Exploration[M].Beijing:PETROLEUM INDUSTRY PRESS,1994.
[2]Pride,S.R.,and J.G.Berryman,Linear dynamics of double-porosity dual-permeability materials.I.Governing equations and acoustic attenuation:Physical Review E 68, 036603(1-10).2003.
[3]张立彬.地震波非弹性吸收与反Q补偿稳定性研究[博士论文].上海:同济大学海洋地质国家重点实验室,2009. Zhang Li Bin.Research on The Seismic Anelatic Absorption and The Stability of Inverse Q Compensation[D].Shanghai: State Key Laboratory of Marine Geology,Tongji University, 2009.
[4]马昭军,刘洋.地震波衰减反演研究综述[J].地球物理学进展, 2005,20(4):1074-1082. MA Zhao-jun,LIU Yang.A summary of research on seismic attenuation[J].Progress in Geophysics(in Chinese),2005, 20(4):1074-1082.
[5]徐平,王宝善,张尉,林建民,陈顒.利用互相关函数求地震波衰减.地球物理学报,2006,49(6):1738-1744. XU Ping,WANG Bao-Shan,ZHANG Wei,LIN Jian-Min, CHEN Yong.Estimating seismic attenuation using crosscorrelation function.Chinese J.Geophys.(in Chinese), 2006,49(6):1738-1744.
[6]李振春,王清振.地震波衰减机理及能量补偿研究综述[J].地球物理学进展,2007,22(4):1147-1152. LI Zhen-chun,WANG Qing-zhen.A review of research on mechanism of seismic attenuation and energy compensation [J],Progress in Geophysics(in Chinese),2007,22(4): 1147-1152.
[7]李子顺.地震波衰减规律及其恢复方法[J].地球物理学进展, 2007,22(5):1545-1551. LI Zi-shun.Seismic wave attenuation and recovering method [J].Progress in Geophysics(in Chinese),2007,22(5): 1545 - 1551.
[8]陈文超,高静怀.基于改进的最佳匹配地震子波的地震资料衰减特性分析.地球物理学报,2007,50(3):837-843. CHEN Wen-Chao,GAO Jing-Huai.Characteristic analysis of seismic attenuation using MBMSW wavelets.Chinese J. Geophys.(in Chinese),2007,50(3):837-843.
[9]魏修成,卢明辉,巴晶,杨慧珠.含黏滞流体各向异性孔隙介质中弹性波的频散和衰减[J].地球物理学报,2008,51(1): 213-220. Wei X C,Lu M H,Ba J,Yang H Z.Dispersion and attenuation of elastic waves in a viscous fluid-saturated anisotropic porous solid[J].Chinese J.Geophys.(in Chinese),2008,50(1): 213-220.
[10]Versteeg,R.,The Marmousi experience:Velocity model determination on a synthetic complex data set.The Leading Edge,13:927-936,1994.
[11]Bjorn Ursin,Tommy Toverud,Comparison of seismic dispersion and attenuation models.CWP408-References, 2002.
[12]Churl-Hyun Jo.An optimal 9-point,finite-difference frequencyspace 2-D scalar wave extrapolator,Geophysics,61(2): 529-537,1996.
[13]George A,Liu J W.Computer Solution of Large Sparse Positive Definite Systems,Prentice Hall,New York,1981.
[14]任浩然,王华忠,龚婷.标量地震波频率-空间域有限差分法数值模拟.石油物探,2009,48(1):20-26. Ren Hao ran,Wang Hua zhong,Gong ting.Scalar seismic wave finite-difference numerical simulation in frequencyspace domain.GEOPHYSICAL PROSPECTING FOR PETROLEUM,2009,48(1):20-26.
[15]廖建平.频率空间域地震波数值模拟方法研究[博士论文]. 上海:同济大学,2009. Liao Jianping.Research on Seismic Wave Numerical Modeling Method in Frequency-space Domain[D]. Shanghai:Tongji University,2009.
[16]廖建平,刘和秀,王华忠,等.快速高精度的频率空间域声波数值模拟方法研究[J].地球物理学进展.2011,26(4):1359 -1363. Liao J P,Liu H X,Wang H,et al.Study on rapid highly accurate acoustic wave numerical simulation in frequency space domain[J].Progress in Geophys.(in Chinese),2011,26 (4):1359-1363.
[17]Slekl and R.Gerhard Pratt.Accurate viscoelastic modeling by frequency-domain finite differences using rotated operators. Geophysics,63(5):1779-1794,1998.
[18]李庆忠.走向精确勘探的道路—高分辨率地震勘探系统工程剖析[M].北京:石油工业出版社,1994. Li Qing-Zhong.The Way to Obtain a Better Resolution in Seismic Prospecting—A Systematical Analysis of High Resolution Seismic Exploration[M].Beijing:PETROLEUM INDUSTRY PRESS,1994.
版权所有：© 2023 中国地质图书馆 中国地质调查局地学文献中心