Seismic forward modeling in fluid-solid media based on equivalent staggered grid scheme
Wu Jianlu1,2, Wu Guochen1,2, Wang Wei3, He Jing3
1. China University of Petroleum(East China), Qingdao, Shandong 266580, China; 2. Laboratory for Marine Mineral Resources, Qingdao National Laboratory for Marine Science and Technology, Qingdao, Shandong 266071, China; 3. Bohai Petroleum Research Institute, Tianjin Branch, CNOOC, Tianjin 300459, China
Abstract:In order to improve the stability and accuracy of seismic simulation for fluid-solid coupling media, the effect of density variation is taken into account based on the equivalent staggered grid.For fluid-solid coupling interface, the first-order stress-velocity equations are selected, in which the fluid domain is represented by second-order pressure equations and solid domain is computed by second-order displacement equations, and the method for coupling the fluid and solid domain by interface is carefully illustrated.The method can save the memory of the modeling parameters and reduce the operation time of seismic modeling.Comparing with the results in fine size model condition, the accuracy of the proposed method is verified.Then the approach is applied on complex Marmousi 2 model.The results show that the proposed method is very stable and precise, in which the accuracy is proved to be the same compared with the results of standard staggered grid method.
Hou G,Wang J and Layton A.Numerical methods for f-structure interaction-a review.Computer Physics Communications,2012,12(2):337-377.
[2]
Virieux J.P-SV wave propagation in heterogeneous media:Velocity-stress finite-difference method.Geophysics,1986,51(4):889-901.
[3]
Burns D R and Stephen R A.Three-dimensional numerical modeling of geoacoustic scattering from seafloor topography.The Journal of the Acoustical Society of America,1990,88(5):2338-2345.
[4]
Dougherty M E and Stephen R A.Seismo/acoustic propagation through rough seafloors.The Journal of the Acoustical Society of America,1991,90(5):2637-2651.
[5]
Stephen R A.A review of finite difference methods for seismo-acoustics problems at the seafloor.Reviews of Geophysics,1988,26(3):445-458.
[6]
Van Vossen R,Robertsson J and Chapman C.Finite-difference modeling of wave propagation in a f-solid configuration.Geophysics,2002,67(2):618-624.
[7]
De Basabe J D and Sen M K.A comparison of finite-difference and spectral-element methods for elastic wave propagation in media with a f-solid interface.Geophysical Journal International,2015,200(1):278-298.
[8]
Stephen R A.A comparison of finite difference and reflectivity seismograms for marine models.Geophy-sical Journal International,1983,72(1):39-57.
[9]
Sochacki J S.Interface conditions for acoustic and elastic wave propagation.Geophysics,1991,56(2):168-181.
[10]
Lee H,Lim S and Min D.2D time-domain acoustic-elastic coupled modeling:a cell-based finite-difference method.Geosciences Journal,2010,13(4):407-414.
[11]
曲英铭,黄建平,李振春等.基于单元交错网格的变坐标系正演模拟方法在声-弹介质中的应用.石油物探,2015,54(5):582-591.Qu Yingming,Huang Jianping and 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.
[12]
Qu Y,Huang J and Li Z et al.A hybrid grid method in an auxiliary coordinate system for irregular f-solid interface modeling.Geophysical Journal International,2017,208(3):1540-1556.
[13]
Komatitsch D,Barnes C and Tromp J.Wave propagation near a fluid-solid interface:A spectral-element approach.Geophysics,2000,65(2):623-631.
[14]
Carcione J M and Helle H B.The physics and simulation of wave propagation at the ocean bottom.Geophysics,2004,69(3):825-839.
宣领宽,张文平,明平剑等.基于不同时间步长时域非结构有限体积法模拟声-弹性耦合问题.固体力学学报,2013,34(2):158-168.Xuan Lingkuan,Zhang Wenping and Ming Pingjian et al.Numerical simulation of acoustic-elastic interaction based on muti-time-step unstructured finite volume time domain method.Chinese Journal of Solid Mechanics.2013,34(2):158-168.
[17]
Choi Y,Min D J and Shin C.Two-dimensional waveform inversion of multi-component data in acoustic-elastic coupled media.Geophysical Prospecting,2008,56(6):863-881.
[18]
Bae H S,Shin C and Cha Y et al.2D acoustic-elastic coupled waveform inversion in the Laplace domain.Geophysical Prospecting,2010,58(6):997-1010.
[19]
Bae H S,Pyun S and Chuang W et al.Frequency-domain acoustic-elastic coupled waveform inversion using the Gauss-Newton conjugate gradient method.Geophysical Prospecting,2011,60(3):413-432.
[20]
Yu P,Geng J and Li X et al.Acoustic-elastic coupled equation for ocean bottom seismic data elastic reverse time migration.Geophysics,2016,81(5):S333-S345.
[21]
Yu P,Geng J and Li X et al.Separating quasi-P-wave in transversely isotropic media with a vertical symmetry axis by synthesized pressure applied to ocean-bottom seismic data elastic reverse time migration.Geophysics,2016,81(6):C295-C307.
[22]
黄金强,李振春和黄建平等.TTI介质Lebedev网格高阶有限差分正演模拟及波型分离.石油地球物理勘探,2017,52(5):915-927.Huang Jinqiang,Li Zhenchun and Huang Jianping et al.Lebedev grid high-order finite-difference modeling and elastic wave-mode separation for TTI media.OGP,2017,52(5):915-927.
[23]
高正辉,孙建国,孙章庆等.基于完全匹配层构建新方法的2.5维声波近似方程数值模拟.石油地球物理勘探,2016,51(6):1128-1133.Gao Zhenghui,Sun Jianguo and Sun Zhangqing et al.2.5D acoustic approximate equation numerical simulation with a new construction method of the perfectly matched layer.OGP,2016,51(6):1128-1133.
[24]
梁文全,王彦飞,杨长春.声波方程数值模拟矩形网格有限差分系数确定法.石油地球物理勘探,2017,52(1):56-62.Liang Wenquan,Wang Yanfei and Yang Changchun.Acoustic wave equation modeling with rectangle grid finite difference operator and its linear time space domain solution.OGP,2017,52(1):56-62.
[25]
Xu K and McMechan G.2D frequency-domain elastic full-waveform inversion using time-domain modeling and a multi step-length gradient approach.Geophy-sics,2014,79(2):R41-R53.
[26]
赵慧,刘颖,李予国.自适应有限元海洋大地电磁场二维正演模拟.石油地球物理勘探,2014,49(3):578-585.Zhao Hui,Liu Ying and Li Yuguo.Adaptive finite element forward modeling for two-dimensional marine magnetotelluric fields.OGP,2014,49(3):578-585.
[27]
黄一凡,胡祥云,韩波.自适应有限元法的大地电磁二维各向异性正反演.石油地球物理勘探,2016,51(4):809-820.Huang Yifan,Hu Xiangyun and Han Bo.Forward and inverse modeling of the magnetotelluric field in 2D anisotropic media with an adaptive finite element method.OGP,2016,51(4):809-820.
[28]
Bartolo L D,Dors C and Mansur W.A new family of finite-difference schemes to solve the heterogeneous acoustic wave equation.Geophysics,2012,77(5):T187-T199.
[29]
Bartolo L D,Dors C and Mansur W.Theory of equi-valent staggered-grid schemes:application to rotated and standard grids in anisotropic media.Geophysical Prospecting,2015,63(5):1097-1125.