基于光流矢量的黏滞声波方程Q补偿逆时偏移Viscoacoustic Wave Equation Q-Compensated Reverse Time Migration Based on Optical Flow Vectors
单闻博,孟建盛,何兵寿
摘要(Abstract):
逆时偏移是解决复杂介质成像问题的有效工具。但由于地下介质的黏滞性,地震波在实际传播过程中会受到地层衰减效应的影响,导致其能量减弱并产生相位频散。在此情况下,基于完全弹性介质假设的逆时偏移会出现深部照明能量不足、偏移剖面能量偏弱及分辨率过低等问题。此外,由于常规互相关成像条件无法避免同路径成像,逆时偏移成像结果浅部会产生强烈的低波数噪声,进而难以获得高精度、高分辨率的成像剖面,因此需要在逆时偏移过程中对地震波的振幅和相位进行校正与补偿,并同时实现低波数噪声的压制。本文从解耦的拉普拉斯黏声方程出发,采用Q补偿算法,并结合光流矢量波场分离技术,研究了一种基于光流矢量的黏滞声波方程Q补偿逆时偏移成像方法,在实现黏声补偿校正的基础上提升了逆时偏移噪声的压制效果。
关键词(KeyWords): 黏滞声波方程;Q补偿逆时偏移;行波分离;光流矢量
基金项目(Foundation): 国家科技重大专项(2024ZD1004201);; 国家自然科学基金项目(42274149)
作者(Author): 单闻博,孟建盛,何兵寿
参考文献(References):
- [1]王胜,何兵寿,陈婷.任意广角方程逆时偏移的脉冲响应及模型试算[J].中国煤炭地质,2014, 26(2):49-54.
- [2]豆辉,徐逸鹤.黏性逆时偏移成像研究进展[J].地球物理学进展,2019, 34(5):1866-1877.
- [3]KJARTANSSON E. Constant Q-wave propagation and attenuation[J]. Journal of Geophysical Research:Solid Earth, 1979, 84(B9):4737-4748.
- [4]DENG F, MCMECHAN G A. True-amplitude prestack depth migration[J]. Geophysics, 2007, 72(3):S155-S166.
- [5]胡江涛,王华忠.基于解析时间波场外推与波场分解的逆时偏移方法研究[J].地球物理学报,2015, 58(8):2886-2895.
- [6]CARCIONE J M, CAVALLINI F, MAINARDI F, et al. Timedomain modeling of constant-Q seismic waves using fractional derivatives[J]. Pure and Applied Geophysics, 2002, 159(7):1719-1736.
- [7]CHEN W, HOLM S. Fractional Laplacian time-space models for linear and nonlinear lossy media exhibiting arbitrary frequency powerlaw dependency[J]. The Journal of the Acoustical Society of America,2004, 115(4):1424-1430.
- [8]TREEBY B E, COX B T. Modeling power law absorption and dispersion for acoustic propagation using the fractional Laplacian[J].The Journal of the Acoustical Society of America, 2010, 127(5):2741-2748.
- [9]ZHU T, HARRIS J M. Modeling acoustic wave propagation in heterogeneous attenuating media using decoupled fractional Laplacians[J]. Geophysics, 2014, 79(3):T105-T116.
- [10]王喜鹏.采用上下行波分离的Q补偿逆时偏移[D].北京:中国科学院大学(中国科学院精密测量科学与技术创新研究院),2023.
- [11]项守龙,鲍五堂,刘丕哲,等.风积沙漠地区三维地震勘探资料处理中Q补偿应用实例[J].中国煤炭地质,2014, 26(9):61-64.
- [12]ZHU T, HARRIS J M, BIONDI B. Q-compensated reverse-time migration[J]. Geophysics, 2014, 79(3):S77-S87.
- [13]吴玉,符力耘,陈高祥.基于分数阶拉普拉斯算子解耦的黏声介质地震正演模拟与逆时偏移[J].地球物理学报,2017, 60(4):1527-1537.
- [14]常海明,何兵寿,杨佳佳,等.逆时偏移中的归一化加权互相关成像条件及应用[J].中国煤炭地质,2014, 26(1):51-55.
- [15]杜启振,朱钇同,张明强,等.叠前逆时深度偏移低波数噪声压制策略研究[J].地球物理学报,2013, 56(7):2391-2401.
- [16]牟海波,何兵寿.地震波逆时偏移中的层位校正与去噪方法[J].中国煤炭地质,2016, 28(4):71-75.
- [17]LIU F, MORTON S A, LEVEILLE J P, et al. Reverse-time migration using one-way wavefield imaging condition[C]//SEG International Exposition and Annual Meeting, 2007:2170-2174.
- [18]XUE H, LIU Y. Reverse-time migration using multidirectional wavefield decomposition method[J]. Applied Geophysics, 2018, 15(2):222-233.
- [19]FEI T W, LUO Y, YANG J, et al. Removing false images in reverse time migration:The concept of de-primary[J]. Geophysics,2015, 80(6):S237-S244.
- [20]王一博,郑忆康,薛清峰,等.基于Hilbert变换的全波场分离逆时偏移成像[J].地球物理学报,2016, 59(11):4200-4211.
- [21]YOON K, MARFURT K J. Reverse-time migration using the Poynting vector[J]. Exploration Geophysics, 2006, 37(1):102-107.
- [22]DU Q, ZHU Y, BA J. Polarity reversal correction for elastic reverse time migration[J]. Geophysics, 2012, 77(2):S31-S41.
- [23]DUAN Y, SAVA P. Scalar imaging condition for elastic reverse time migration[J]. Geophysics, 2015, 80(4):S127-S136.
- [24]HORN B K P, SCHUNCK B G. Determining optical flow[J].Artificial intelligence, 1981, 17(1-3):185-203.
- [25]ZHANG Q. RTM angle gathers and specular filter(SF)RTM using optical flow[C]//SEG International Exposition and Annual Meeting,2014:3816-3820.
- [26]吴成梁,王华忠,冯波,等.基于CLG光学流和波场分解的逆时偏移角度道集提取方法研究[J].地球物理学报,2021, 64(4):1375-1388.
- [27]解闯,宋鹏,邹志辉,等.基于光流矢量行波分离的相关加权逆时偏移成像[J].地球物理学报,2022, 65(6):2260-2275.
- [28]冯彦苍,刘文卿,张会星,等.基于光流矢量的TI介质逆时偏移方法[J].岩性油气藏,2026, 38(1):115-125.
- [29]GAZDAG J. Modeling of the acoustic wave equation with transform methods[J]. Geophysics, 1981, 46(6):854-859.
- [30]谢志南,廖振鹏.多次透射公式在波动数值模拟中的一种实现方案[J].固体力学学报,2010, 31(4):422-426.