一阶速度—压力常分数阶黏滞声波方程及其数值模拟

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

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

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

摘要与传统的整数阶黏滞波动方程相比，分数阶拉普拉斯算子黏滞方程能更准确地匹配目前广泛使用的常Q模型，而且分数阶黏滞波动方程中控制振幅衰减和相位变化的算子是显式分离的，这对于发展稳定的衰减补偿逆时偏移算法至关重要。首先基于时间域二阶位移形式的常分数阶拉普拉斯算子黏滞声波方程，推导了一阶速度-压力形式常分数阶拉普拉斯算子黏滞声波方程；为了模拟更加真实的振幅变化信息，在新的黏滞声波方程中考虑了密度空变的影响；为了避免由傅里叶变换的周期性而引入的虚假反射，提出了一种适用于分数阶黏滞声波方程的卷积型完全匹配层（CPML）吸收边界加载方法；最后采用交错网格伪谱法进行数值模拟。均匀介质中数值解与解析解的对比证实了该一阶速度-压力常分数阶黏滞声波方程能准确描述常Q模型，BP盐丘模型的地震波场模拟结果证实了其对复杂介质的适用性。

	服务

	把本文推荐给朋友
	加入我的书架
	加入引用管理器
	E-mail Alert
	RSS
	作者相关文章
	陈汉明
	汪燚林
	周辉

关键词 ：黏滞声波方程, 数值模拟, 分数阶拉普拉斯算子, 交错网格, 伪谱法, 完全匹配层

Abstract：We develop a viscoacoustic wave equation with fractional-order Laplacians,which is better than the traditional integral-order viscoacoustic equation because the new equation more accurately describes the widely used constant-Q model.The operators related to amplitude attenuation and phase change are explicitly independent of each other,which is important for the robust reverse time migration with attenuation compensation.We formulate the first-order velocity-pressure viscoacoustic wave equation with the constant fractional-order Laplacians based on the second-order displacement viscoacoustic equation with the constant fractional-order Laplacians in time domain.To better model amplitude variation,space-varying density is involved in the new equation.To avoid spurious reflections caused by the periodicity of the Fourier transform,a convolutional perfectly matched layer (CPML) is employed as the absorbing boundary for the fractional-order Laplacian viscoacoustic equation.Numerical simulations are fulfilled using staggered-mesh pseudo-spectral method.We compare the numerical solution with the analytic solution for the homogeneous medium,and we find the new equation accurately describes the constant-Q model.We also verify its feasibility for complicated media through seismic wave field simulation using the BP salt dome model.

Key words： viscoacoustic wave equation numerical simulation fractional-order Laplacians staggered mesh pseudo-spectral method perfectly matched layer

收稿日期: 2019-08-13

PACS:

P631

基金资助:本项研究受国家重点研发计划项目“多信息相容约束高效全波形反演方法研究”（2018YFA0702502）和国家自然科学基金项目“基于高维最优传输距离目标函数的衰减介质全波形反演研究”（41804111）、“变分数阶拉普拉斯算子粘滞声波方程正演、逆时偏移和全波形反演研究”（41630314）、“弹簧网络模型和格子玻尔兹曼模型耦合的含流体孔隙介质波场模拟方法研究”（41874130）联合资助。

通讯作者: 周辉,北京市昌平区府学路18号中国石油大学(北京)地球物理学院地球物理系,102249。Email:huizhou@cup.edu.cn E-mail: huizhou@cup.edu.cn

作者简介: 陈汉明,博士,1987年生;2010年获长江大学勘查技术与工程专业学士学位,2013年获中国石油大学(北京)地球探测与信息技术专业硕士学位,2017年获中国石油大学(北京)地质资源与地质工程专业博士学位,并留校工作。现为中国石油大学(北京)地球物理学院地球物理系讲师、SEG/EAGE会员。主要从事复杂波动方程数值模拟、复杂介质波场传播、波动方程偏移成像和全波形反演方面的研究。目前已在国内外SCI/EI期刊以第一作者发表论文十余篇。

引用本文:

陈汉明, 汪燚林, 周辉. 一阶速度—压力常分数阶黏滞声波方程及其数值模拟[J]. 石油地球物理勘探, 2020, 55(2): 302-310. CHEN Hanming, WANG Yilin, ZHOU Hui. A novel constant fractional-order Laplacians viscoacoustic wave equation and its numerical simulation method. Oil Geophysical Prospecting, 2020, 55(2): 302-310.

链接本文:

http://www.ogp-cn.com/CN/10.13810/j.cnki.issn.1000-7210.2020.02.009 或 http://www.ogp-cn.com/CN/Y2020/V55/I2/302

[1]	Chai X,Wang S,Yuan S,et al.Sparse reflectivity inversion for nonstationary seismic data[J].Geophy-sics,2014,79(3):V93-V105.
[2]	Yuan S,Wang S,Tian N,et al.Stable inversion-based multitrace deabsorption method for spatial continuity preservation and weak signal compensation[J].Geophysics,2016,81(3):V199-V212.
[3]	李海山,杨午阳,雍学善.二维时间域黏声波全波形反演[J].石油地球物理勘探,2018,53(1):87-94.LI Haishan,YANG Wuyang,YONG Xueshan.Visco-acoustic wave full waveform inversion in the 2D time domain[J].Oil Geophysical Prospecting,2018,53(1):87-94.
[4]	Kjartansson E.Constant Q-wave propagation and attenuation[J].Journal of Geophysical Research,1979,84(B9):4737-4748.
[5]	孙成禹,印兴耀.三参数常Q黏弹性模型构造方法研究[J].地震学报,2007,29(4):348-357.SUN Chengyu,YIN Xinyao.Construction of constant-Q viscoelastic model with three parameters[J].Acta Seismologica Sinica,2007,29(4):348-357.
[6]	杜启振.各向异性黏弹性介质伪谱法波场模拟[J].物理学报,2004,53(12):4428-4434.DU Qizhen.Wavefield forward modeling with the pseudo-spectral method in viscoelastic and azimuthally anisotropic media[J].Acta Physica Sinica,2004,53(12):4428-4434.
[7]	蔡瑞乾,孙成禹,伍敦仕,等.黏声波动方程变机制数有限差分正演[J].石油地球物理勘探,2019,54(3):529-538.CAI Ruiqian,SUN Chengyu,WU Dunshi,et al.Finite-difference numerical modeling with variable mechanisms for viscoacoustic wave equation[J].Oil Geophysical Prospecting,2019,54(3):529-538.
[8]	苑春方,彭苏萍,张中杰,等.Kelvin-Voigt均匀黏弹性介质中传播的地震波[J].中国科学(D辑),2006,35(10):957-962.YUAN Chunfang,PENG Suping,ZHANG Zhongjie,et al.Seismic wave propagation in Kelvin-Voigt viscoelastic media[J].Science in China (Series D),2006,35(10):957-962.
[9]	李晓波,董良国.黏弹介质中可变网格地震波传播数值模拟[J].石油物探,2012,51(1):1-10.LI Xiaobo,DONG Liangguo.Variable-grid numerical simulation of seismic wave propagation in visco-elastic media[J].Geophysical Prospecting for Petroleum,2012,51(1):1-10.
[10]	孙成禹,肖云飞,印兴耀,等.黏弹介质波动方程有限差分解的稳定性研[J].地震学报,2010,32(2):147-156.SUN Chengyu,XIAO Yunfei,YIN Xingyao,et al.Study on the stability of finite difference solution of visco-elastic wave equations[J].Acta Seismologica Sinica,2010,32(2):147-156.
[11]	Chen W,Holm S.Fractional Laplacian time-space mo-dels 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.
[12]	Carcione J M,Cavalint F,Mainardi F,et al.Time-domain modeling of constant-Q seismic waves using fractional derivative[J].Pure and Applied Geophy-sics,2002,159(7):1719-1736.
[13]	孙成禹,乔志浩,伍敦仕,等.常Q衰减介质分数阶波动方程优化有限差分模拟[J].地震学报,2017,39(3):343-355.SUN Chengyu,QIAO Zhihao,WU Dunshi,et al.Modeling of wave equation with fractional derivative using optimal finite-difference method in constant-Q attenuation media[J].Acta Seismologica Sinica,2017,39(3):343-355.
[14]	刘财,胡宁,郭智奇,等.基于分数阶时间导数常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 viscoelastic constitutive relation with a fractional temporal derivative[J].Chinese Journal of Geophysics,2018,61(6):2446-2458.
[15]	Carcione J M.A generalization of the Fourier pseudospectral method[J].Geophysics,2010,75(6):A53-A56.
[16]	Zhang Y,Zhang P,Zhang H.Compensating for visco-acoustic effects in reverse-time migration[C].SEG Technical Program Expanded Abstracts,2010,29:3160-3164.
[17]	Zhu T,Harris J M.Modeling acoustic wave propagation in heterogeneous attenuating media using decoupled fractional Laplacians[J].Geophysics,2014,79(3):T105-T116.
[18]	吴玉,符力耘,陈高祥.基于分数阶拉普拉斯算子解耦的黏声介质地震正演模拟与逆时偏移[J].地球物理学报,2017,60(4):1527-1537.WU Yu,FU Liyun,CHEN Gaoxiang.Forward mode-ling and reverse time migration of viscoacoustic media using decoupled fractional Laplacians[J].Chinese Journal of Geophysics,2017,60(4):1527-1537.
[19]	Treeby B E,Cox B T.Modeling power law absorption and dispersion for acoustic propagation using the fractional Laplacian[J].Journal of the Acoustical Society of America,2010,127(5):2741-2748.
[20]	Zhu T,Harris J M,Biondi B.Q-compensated reverse-time migration[J].Geophysics,2014,79(3):S77-S87.
[21]	Chen H,Zhou H,Li Q,et al.Two efficient modeling schemes for fractional Laplacian viscoacoustic wave equation[J].Geophysics,2016,81(5):T233-T249.
[22]	Chen H,Zhou H,Jiang S,et al.Fractional Laplacian viscoacoustic wave equation low-rank temporal extrapolation[J].IEEE Access,2019,7(1):93187-93197.
[23]	Li Q,Zhou H,Zhang Q,et al.Efficient reverse time migration based on fractional Laplacian viscoacoustic wave equation[J].Geophysical Journal International,2016,204(1):488-504.
[24]	Wang N,Zhou H,Chen H,et al.A constant fractional-order viscoelastic wave equation and its numerical simulation scheme[J].Geophysics,2018,83(1):T39-T48.
[25]	Wang Y,Zhou H,Chen H,et al.Adaptive stabilization for Q-compensated reverse time migration[J].Geophysics,2018,83(1):S15-S32.
[26]	Chen H,Zhou H,Zu S.Simultaneous inversion of velocity and Q using a fractional Laplacian constant-Q wave equation[C].Extended Abstracts of 79th EAGE Conference and Exhibition,2017,ThP409.
[27]	Chen H,Zhou H,Rao Y,et al.A matrix-transform numerical solver for fractional Laplacian viscoacoustic wave equation[J].Geophysics,2019,84(4):T283-T297.
[28]	Chen H,Zhou H,Li Y.Application of unsplit convo-lutional perfectly matched layer for scalar arbitrarily wide-angle wave equation[J].Geophysics,2014,79(6):T313-T321.
[29]	Komatitsch D,Martin R.An unsplit convolutional perfectly matched layer improved at grazing incidence for the seismic wave equation[J].Geophysics,2007,72(5):SM155-SM167.
[30]	Carcione J M,Kosloff D,Kosloff R.Wave propagation simulation in a linear viscoelastic medium[J].Geophysical Journal International,1988,95(3):597-611.