基于上下行波分解的分数阶拉普拉斯算子黏滞声波逆时偏移成像方法

doi:10.13810/j.cnki.issn.1000-7210.2023.04.012

摘要
图/表
参考文献
相关文章 (15)

全文: PDF (4493 KB) HTML []
输出: BibTeX | EndNote (RIS)

摘要地下介质通常存在黏滞性，使地震波在传播过程中产生能量衰减及（相位）频散效应，降低了偏移成像照明能量。采用互相关成像条件的双程波动方程逆时偏移技术，当速度模型精度不高或存在较强速度梯度时，可能会产生低频噪声和假象。为此，提出一种基于上下行波分解的解耦分数阶拉普拉斯算子黏滞声波逆时偏移成像（Q补偿逆时偏移）方法，采用解耦的分数阶拉普拉斯算子黏滞声波方程分别对震源波场和检波点波场进行外推，在此过程中校正地下介质对地震波产生的衰减效应；引入自适应稳定因子应对地震波衰减补偿过程中数值不稳定问题；采用时间解析波场外推方法，在时间切片上分别对震源、检波点波场的传播方向进行分解；最终运用互相关成像条件对下行的震源波场和上行的检波点波场进行成像。数值测试结果表明，所采用的Q补偿逆时偏移方法可校正由地下介质黏性特性引起的地震波振幅衰减和（相位）频散效应，该成像方法能显著减少传统逆时偏移成像中的低频噪声与假象。

	服务

	把本文推荐给朋友
	加入我的书架
	加入引用管理器
	E-mail Alert
	RSS
	作者相关文章
	王喜鹏
	张庆臣
	毛伟建

关键词 ： Q补偿逆时偏移, 波场分解, 低频噪声, 分数阶

Abstract：The subsurface media typically exhibit viscosity, which causes energy attenuation and(phase) dispersion effects during seismic wave propagation and in turn reduces the illumination energy of migration imaging. The twoway wave equation reverse-time migration with cross-correlation imaging conditions may produce low-frequency noise and false images when there is no highly accurate velocity models or strong velocity gradients. To this end, this paper proposes a viscoacoustic reverse-time migration method based on decoupled fractional Laplacian with up and down-going wave decomposition(Q-compensated reverse-time migration method). The viscoacoustic wave equation based on decoupled fractional Laplacian is employed to separately extrapolate the source wavefield and receiver wavefield and correct the attenuation effects generated by subsurface media during the process. An adaptive stability factor is introduced to handle numerical instability arising from the attenuation compensation of seismic waves. The time-domain analytical wavefield extrapolation method is adopted to separately decompose the propagation directions of the source and receiver wavefields on each time slice. Finally, the cross-correlation imaging condition is utilized for imaging of the down-going source wavefield and up-going receiver wavefield. Numerical tests show that the Q-compensated reverse-time migration method can correct the amplitude attenuation and(phase) dispersion effects caused by the viscosity characteristics of subsurface media. This imaging method significantly reduces the low-frequency noise and false images in traditional reverse-time migration imaging.

Key words： Q-compensated reverse-time migration wavefield decomposition low-frequency noise fractional order

收稿日期: 2022-09-17

PACS:

P631

基金资助:本项研究受国家自然科学基金重点项目“海底四分量高精度地震偏移成像方法研究”（42130808）资助。

通讯作者: 毛伟建,湖北省武汉市洪山区徐东大街340号中国科学院精密测量科学与技术创新研究院计算与勘探地球物理研究中心,430077。Email:wjmao@whigg.ac.cn E-mail: wjmao@whigg.ac.cn

作者简介: 王喜鹏 1996年生;2020年毕业于东北石油大学,获地球物理学专业学士学位;2023年获中国科学院大学固体地球物理学专业硕士学位;现入职中国石油集团东方地球物理勘探有限责任公司,主要从事黏声介质正演模拟、逆时偏移成像与噪声压制方面的研究。

引用本文:

王喜鹏, 张庆臣, 毛伟建. 基于上下行波分解的分数阶拉普拉斯算子黏滞声波逆时偏移成像方法[J]. 石油地球物理勘探, 2023, 58(4): 872-882. WANG Xipeng, ZHANG Qingchen, MAO Weijian. Viscoacoustic reverse-time migration imaging based on fractional Laplacian with up-going and down-going wave decomposition. Oil Geophysical Prospecting, 2023, 58(4): 872-882.

链接本文:

http://www.ogp-cn.com/CN/10.13810/j.cnki.issn.1000-7210.2023.04.012 或 http://www.ogp-cn.com/CN/Y2023/V58/I4/872

[1]	陈生昌, 刘亚楠, 李代光. 基于波动方程偏移的地震数据地层成像[J]. 石油地球物理勘探, 2022, 57(5):1105-1113.CHEN Shengchang, LIU Yanan, LI Daiguang. Stratigraphic imaging with seismic data based on wave equation migration[J]. Oil Geophysical Prospecting, 2022, 57(5):1105-1113.
[2]	吴丹, 吴海莉, 李群, 等. 应用自适应矩估计的快速最小二乘逆时偏移[J]. 石油地球物理勘探, 2022, 57(2):386-394.WU Dan, WU Haili, LI Qun, et al. Fast least-squares reverse time migration with adaptive moment estimation[J]. Oil Geophysical Prospecting, 2022, 57(2):386-394.
[3]	ZHOU J, BIRDUS S, HUNG B, et al. Compensating attenuation due to shallow gas through Q tomography and Q-PSDM, a case study in Brazil[C]. SEG Technical Program Expanded Abstracts, 2011, 30:3332-3336.
[4]	KJARTANSSON E. Constant Q-wave propagation and attenuation[J]. Journal of Geophysical Research:Solid Earth, 1979, 84(B9):4737-4748.
[5]	LIU H P, ANDERSON D L, KANAMORI H. Velocity dispersion due to anelasticity;implications for seismology and mantle composition[J]. Geophysical Journal International, 1976, 47(1):41-58.
[6]	EMMERICH H, KORN M. Incorporation of attenuation into time-domain computations of seismic wave fields[J]. Geophysics, 1987, 52(9):1252-1264.
[7]	PAPOULIA K D, PANOSKALTSIS V P, KURUP N V, et al. Rheological representation of fractional order viscoelastic material models[J]. Rheologica Acta, 2010, 49(4):381-400.
[8]	MOCZO P, KRISTEK J. On the rheological models used for time-domain methods of seismic wave propagation[J]. Geophysical Research Letters, 2005, 32(1):L01306.
[9]	BLANCH J O, ROBERTSSON J O A, SYMES W W. Modeling of a constant Q:Methodology and algorithm for an efficient and optimally inexpensive viscoelastic technique[J]. Geophysics, 1995, 60(1):176-184.
[10]	CHEN W, HOLM S. Fractional Laplacian time-space models for linear and nonlinear lossy media exhibiting arbitrary frequency power-law dependency[J]. The Journal of the Acoustical Society of America, 2004, 115(4):1424-1430.
[11]	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.
[12]	刘文卿, 王宇超, 雍学善, 等. 基于GPU/CPU叠前逆时偏移研究及应用[J]. 石油地球物理勘探, 2012, 47(5):712-716.LIU Wenqing, WANG Yuchao, YONG Xueshan, et al. Prestack reverse time migration on GPU/CPU co-parallel computation[J]. Oil Geophysical Prospecting, 2012, 47(5):712-716.
[13]	ZHU T, HARRIS J M. Modeling acoustic wave propagation in heterogeneous attenuating media using decoupled fractional Laplacians[J]. Geophysics, 2014, 79(3):T105-T116.
[14]	ZHU T, HARRIS J M, BIONDI B. Q-compensated reverse-time migration[J]. Geophysics, 2014, 79(3):S77-S87.
[15]	LI Q, ZHOU H, ZHANG Q, et al. Efficient reverse time migration based on fractional Laplacian viscoacoustic wave equation[J]. Geophysical Journal International, 2015, 204(1):488-504.
[16]	TREEBY B E. Acoustic attenuation compensation in photoacoustic tomography using time-variant filtering[J]. Journal of Biomedical Optics, 2013, 18(3):036008.
[17]	WANG Y. Inverse Q-filter for seismic resolution enhancement[J]. Geophysics, 2006, 71(3):V51-V60.
[18]	WANG Y, ZHOU H, CHEN H, et al. Adaptive stabilization for Q-compensated reverse time migration[J]. Geophysics, 2018, 83(1):S15-S32.
[19]	陈汉明, 周辉, 田玉昆. 分数阶拉普拉斯算子黏滞声波方程的最小二乘逆时偏移[J]. 石油地球物理勘探, 2020, 55(3):616-626.CHEN Hanming, ZHOU Hui, TIAN Yukun. Least-squares reverse-time migration based on a fractional Laplacian viscoacoustic wave equation[J]. Oil Geophysical Prospecting, 2020, 55(3):616-626.
[20]	CHEN H, ZHOU H, RAO Y. An implicit stabilization strategy for Q-compensated reverse time migration[J]. Geophysics, 2020, 85(3):S169-S183.
[21]	YOON K, MARFURT K J. Reverse-time migration using the Poynting vector[J]. Exploration Geophysics, 2006, 37(1):102-107.
[22]	康玮, 程玖兵. 叠前逆时偏移假象去除方法[J]. 地球物理学进展, 2012, 27(3):1163-1172.KANG Wei, CHENG Jiubing. Methods of suppressing artifacts in prestack reverse time migration[J]. Progress in Geophysics, 2012, 27(3):1163-1172.
[23]	GUITTON A, KAELIN B, BIONDI B. Least-squares attenuation of reverse-time-migration artifacts[J]. Geophysics, 2007, 72(1):S19-S23.
[24]	ZHANG Y, SUN J. Practical issues in reverse time migration:true amplitude gathers, noise removal and harmonic source encoding[J]. First Break, 2009, 27(1):53-60.
[25]	LIU F, ZHANG G, MORTON S A, et al. Reverse-time migration using one-way wavefield imaging condition[C]. SEG Technical Program Expanded Abstracts, 2007, 26:2170-2174.
[26]	LIU F, ZHANG G, MORTON S A, et al. An effective imaging condition for reverse-time migration using wavefield decomposition[J]. Geophysics, 2011, 76(1):S29-S39.
[27]	杜启振, 朱钇同, 张明强, 等. 叠前逆时深度偏移低频噪声压制策略研究[J]. 地球物理学报, 2013, 56(7):2391-2401.DU Qizhen, ZHU Yitong, ZHANG Mingqiang, et al. A study on the strategy of low wavenumber noise suppression for prestack reverse-time depth migration[J]. Chinese Journal of Geophysics, 2013, 56(7):2391-2401.
[28]	FEI T W, LUO Y, YANG J, et al. Removing false images in reverse time migration:The concept of deprimary[J]. Geophysics, 2015, 80(6):S237-S244.
[29]	胡江涛, 王华忠. 基于解析时间波场外推与波场分解的逆时偏移方法研究[J]. 地球物理学报, 2015, 58(8):2886-2895.HU Jiangtao, WANG Huazhong. Reverse time migration using analytical time wavefield extrapolation and separation[J]. Chinese Journal of Geophysics, 2015, 58(8):2886-2895.
[30]	ZHOU H, LIU Y, WANG J. Viscoacoustic reverse time migration by attenuation compensation and wavefield decomposition[C]//SEG 2018 Workshop:SEG Seismic Imaging Workshop, Beijing, China, 2019, 90-93.
[31]	BLAND D R. The Theory of Linear Viscoelasticity[M]. Dover Publications Inc., Mineola, 2016.
[32]	CLAERBOUT J F. Toward a unified theory of reflector mapping[J]. Geophysics, 1971, 36(3):467-481.
[33]	HU L Z, MCMECHAN G A. Wave-field transformations of vertical seismic profiles[J]. Geophysics, 1987, 52(3):307-321.