Forward modeling and illumination analysis in visco-acoustic VTI media
Li Jinli1,2, Liu Jianxun1,2, Jiang Chunxiang1,2, Wang Xiaojiang1,2, Yue Hangyu1,2, Li Guangcai1,2
1. The National Center for Geological Exploration Technology, Langfang, Hebei 065000, China; 2. Institute of Geophysical and Geochemical Exploration, Chinese Academy of Geological Sciences, Langfang, Hebei 065000, China
Abstract:The visco-elasticity and anisotropy widely exist in the earth.Therefore,the attenuation term of visco-acoustic quasi-differential equation is extended to VTI quasi-differential equations to propose visco-acoustic VTI quasi-differential equations,which care for both visco-elastic and VTI properties of media.The finite difference method is used to solve the equations,in which attenuation terms are solved in the wave-number domain,while other terms are solved in the time-space domain.On this basis,we develop a two-way wave illumination analysis method in visco-acoustic VTI media.A homogeneous model forward modeling proves that the proposed visco-acoustic VTI quasi-differential equations can describe dynamics and kinematics characteristics of seismic wave propagating in the visco-acoustic VTI media.The forward modeling for a sag model shows that the VTI property causes traveltime and waveform difference,and the visco-elasticity leads to amplitude attenuation.Finally,the illumination analysis illustrates that the visco-elastic and VTI properties affect the energy distribution of the wave field.
Qu Yingming,Li Zhenchun,Huang Jianping et al.Viscoacoustic anisotropic full waveform inversion.Journal of Applied Geophysics,2017,136:484-497.
[2]
Day S M,Minster J B.Numerical simulation of wavefields using a Pade approximant method.Geophysics Journal Royal Astronomical Society,1984,78(1):105-118.
[3]
Emmerich H,Korn M.Incorporation of attenuation into time-domain computations of seismic wave fields.Geophysics,1987,52(9):1252-1264.
[4]
Carcione J M.Viscoacoustic wave propagation simulation in the earth.Geophysics,1988,53(6):769-777.
[5]
Tal-Ezer H,Carcione J M,Kosloff D.An accurate and efficient scheme for wave propagation in linear visco-elastic media.Geophysics,1990,55(10):1366-1379.
[6]
Robertsson J O A,Blanch J O,Symes W W.Visco-elastic finite-difference modeling.Geophysics,1994,59(9):1444-1456.
[7]
奚先,姚姚.二维粘弹性随机介质中的波场特征分析.地球物理学进展,2004,19(3):608-615.Xi Xian,Yao Yao.The analysis of wave field characteristics in 2-D viscoelastic random medium.Progress in Geophysics,2004,19(3):608-615.
[8]
邵志刚,傅容珊,黄建华等.线性粘弹性介质中地震波场数值模拟.地球物理学进展,2007,22(4):1135-1141.Shao Zhigang,Fu Rongshan,Huang Jianhua et al.Modeling seismic wave in linear viscoelastic media.Progress in Geophysics,2007,22(4):1135-1141.
[9]
郭振波,田坤,李振春等.利用非线性最优化方法构建近似常Q黏弹模型.中国石油大学学报(自然科学版),2014,38(2):52-58.Guo Zhenbo,Tian Kun,Li Zhenchun et al.Constructing nearly constant-Q viscoelastic model using the nonlinear optimization method.Journal of China University of Petroleum,2014,38(2):52-58.
[10]
王小杰,栾锡武,王延光等.黏弹性介质叠前地震反演方法.石油地球物理勘探,2016,51(3):544-555.Wang Xiaojie,Luan Xiwu,Wang Yanguang et al.Prestack seismic inversion of viscoelastic medium.OGP,2016,51(3):544-555.
[11]
Helbig K.Elliptical anisotropy-Its significance and meaning.Geophysics,1983,48(7):825-832.
[12]
Dellinger J and Muir F.Imaging reflections in elliptically anisotropic media.Geophysics,1988,53(12):1616-1618.
Cphen J K.Analytic study of the effective parameters for determination of the NMO velocity function in transversely isotropic media.Geophysics,1997,62(6):1855-1866.
[15]
Alkhalifah T.Acoustic approximations for processing in transversely isotropic media.Geophysics,1998,63(2):623-631.
[16]
Alkhalifah T.An acoustic wave equation for anisotropic media.Geophysics,2000,65(4):1239-1250.
[17]
Grechka V,Zhang L and Rector J W.Shear waves in acoustic anisotropic media.Geophysics,2004,69(2):576-582.
[18]
Duveneck E,Bakker P M.Stable P-wave modeling for reverse-time migration in tilted TI media.Geophysics,2011,76(2):S65-S75.
[19]
Zhang Y,Zhang H Z and Zhang G Q.A stable TTI reverse time migration and its implementation.Geo-physics,2011,76(3):WA3-WA11.
[20]
Chu C,Macy B K and Anno P D.Pure acoustic wave propagation in transversely isotropic media by the pseudospectral method.Geophysical Prospecting,2013,61(3):556-567.
[21]
Zhan G,Pestana R C and Stoffa P L.Decoupled equations for reverse time migration in tilted transversely isotropic media.Geophysics,2012,77(2):T37-T45.
[22]
Zhan G,Pestana R C and Stoffa P L.An efficient hybrid pseudo spectral/finite-difference scheme for solving the TTI pure P-wave equation.Journal of Geophysics and Engineering,2013,10(2):025004,doi:10.1088/1742-2132/1012/025004.
[23]
李娜,李振春,黄建平等.Lebedev网格与标准交错网格耦合机制下的复杂各向异性正演模拟.石油地球物理勘探,2014,49(1):121-131.Li Na,Li Zhenchun,Huang Jianping et al.Numerical simulation with coupling Lebedev and standard staggered grid schemes for complex anisotropic media.OGP,2014,49(1):121-131.
[24]
曲英铭,黄建平,李振春等.一种基于非规则网格的地震波射线追踪方法.石油物探,2014,53(6):627-632.Qu Yingming,Huang Jianping,Li Zhenchun et al.A seismic ray tracing method basd on irregular grids.GPP,2014,53(6):627-632.
[25]
吴国忱,梁锴.VTI介质频率-空间域准P波正演模拟.石油地球物理勘探,2005,40(5):63-73.Wu Guochen,Liang Kai.Quasi P-wave forward mode-ling in frequency-space domain in VTI media.OGP,2005,40(5):63-73.
[26]
Huang Jianping,Qu Yingming,Li Zhenchun et al.A hybrid grid method in an auxiliary coordinate system for irregular fluid-solid interface modeling.Geophysical Journal International,2017,208(3):1540-1556.
[27]
曲英铭,黄建平,李振春等.分层坐标变换法起伏自由地表弹性波叠前逆时偏移.地球物理学报,2015,58(8):544-560.Qu Yingming,Huang Jianping,Li Zhenchun et al.Elastic wave forward modeling and pre-stack reverse time migration of irregular free-surface based on la-yered mapping method.Chinese Journal of Geophysics,2015,58(8):544-560.
[28]
Huang J,Qu Y,Li Z et al.Variable coordinate system forward modeling method for irregular surface based on dual-variable grid.Applied Geophysics,2015,12(1):101-110.
[29]
曲英铭,黄建平,李振春等.分层映射法起伏自由地表弹性波正演模拟与波场分离.石油地球物理勘探,2016,51(2):261-271.Qu Yingming,Huang Jianping,Li Zhenchun et al.Elastic wave forward modeling and wave-field separation of irregular free-surface based on the layered mapping method.OGP,2016,51(2):261-271.
[30]
曲英铭,黄建平,李振春等.基于单元交错网格的变坐标系正演模拟方法在声-弹介质中的应用.石油物探,2015,54(5):582-591.Qu Yingming,Huang Jianping,Li Zhenchun et al.The application of variable coordinate system forward modeling using cell-based staggered grids in acoustic-elastic coupled medium.GPP,2015,54(5):582-591.
[31]
刘定进,蒋波,李博等.起伏地表逆时偏移在复杂山前带地震成像中的应用.石油地球物理勘探,2016,51(2):315-324.Liu Dingjin,Jiang Bo,Li Bo et al.Rugged-topography reverse time migration in complex piedmont zone.OGP,2016,51(2):315-324.
[32]
Bear G,Lu C,Lu R et al.The construction of subsurface illumination and amplitude maps via ray tracing.The Leading Edge,2000,19(7):726-728.
[33]
Muerdter D and Ratcliff D.Understanding subsalt illumination through ray-trace modeling,Part 1:Simple 2-D salt models.The Leading Edge,2001,20(6):578-594.
[34]
Muerdter D,Kelly M and Ratcliff D.Understanding subsalt illumination through ray-trace modeling,Part 2:Dipping salt bodies,salt peaks,and nonreciprocity of subsalt amplitude response.The Leading Edge,2001,20(7):688-697.
[35]
Gelius L J,Lecomte I and Tabti H.Analysis of the resolution function in seismic prestack depth imaging.Geophysical Prospecting,2002,50(5):505-515.
[36]
Hoffmann J.Illumination,resolution,and image quality of PP-and PS-waves for survey planning.The Leading Edge,2000,20(9):1008-1014.
[37]
Wu R S,Chen L.Directional illumination analysis using beamlet decomposition and propagation.Geo-physics,2006,71(4):S147-S159.
[38]
Xie X B,Jin S W,Wu R S et al.Wave-equation-based seismic illumination analysis.Geophysics,2006,71(5):S169-S177.
[39]
陈雪,韩立国,杨贺龙等.粘弹VTI介质单程波正演模拟与照明分析.石油物探,2015,54(6):643-651.Chen Xue,Han Liguo,Yang Helong et al.One-way wave equation forward modeling and illumination analysis in viscoelastic VTI media.GPP,2015,54(6):643-651.
[40]
Cao J,Wu R S.Full-wave directional illumination analysis in the frequency domain.Geophysics,2009,74(4):165-165.
[41]
雷涛,顾汉明,李列等.面向目标靶区的双程波动方程地震定向照明分析.石油物探,2014,53(4):437-443.Lei Tao,Gu Hanming,Li Lie et al.Target-oriented seismic dirctional illumination analysis based on two-way acoustic wave equation.GPP,2014,53(4):437-443.
[42]
单联瑜,刘洪,匡斌等.基于胜利典型地质模型的波动方程地震波照明分析研究.石油地球物理勘探,2009,44(1):1-6.Shan Lianyu,Liu Hong,Kuang Bin et al.Analysis on seismic wave illumination of wave equation based on Shengli typical geologic model.OGP,2009,44(1):1-6.
[43]
赵贤正,刘国昌,史原鹏等.基于局部相似属性的数据依赖地震照明分析方法.石油地球物理勘探,2016,51(6):1062-1068.Zhao Xianzheng,Liu Guochang,Shi Yuanpeng et al.Seismic illumination analysis with local similarity.OGP,2016,51(6):1062-1068.
[44]
Carcione J M.A generalization of the Fourier pseudospectral method.Geophysics,2010,75(6):A53-A56.
[45]
Bai J Y,Chen G Q,Yingst D et al.Attenuation compensation in viscoacoustic reserve-time migration.SEG Technical Program Expanded Abstracts,2013,32:3825-3830.
[46]
Suh S,Yoon K,Cai J et al.Compensating visco-acoustic effects in anisotropic reverse-time migration.SEG Technical Program Expanded Abstracts,2012,31:1-5.
[47]
Duvenck E and Bakker P M.Stable P-wave modeling for reverse-time migration in tilted TI media.Geo-physics,2011,76(2):S65-S75.
[48]
徐文才,杨国权,李振春等.基于GSLS模型TI介质衰减拟声波方程.地球物理学报,59(6):2232-2244.Xu Wencai,Yang Guoquan,Li Zhenchun et al.Pseudo acoustic equation for TI medium attenuation based on GSLS model.Chinese Journal of Geophysics,59(6):2232-2244.
[49]
Collino F and Tsogka C.Application of the perfectly matched absorbing layer model to the linear elasto-dynamic problem in anisotropic heterogeneous media.Geophysics,2001,66(1):294-307.