Acoustic wave numerical simulation in pseudo-depth domain
Li Qingyang1, Huang Jianping1, Li Zhenchun1, Li Na2, Qu Yingming1, Su Yun2
1. School of Geoscience, China University of Petroleum (East China), Qingdao, Shandong 266580, China;
2. Geophysical Exploration Research Institute, Zhongyuan Oilfield Branch Co., SINOPEC, Puyang, Henan 457001, China
Abstract:Surface seismic exploration usually meets challenges like low-velocity bodies, surface topography, steep structures, fractures and so on. Seismic wavefield simulation with conventional Cartesian finite difference forward modeling method in these areas often needs very small grid spacing to ensure the accuracy and stability. So it leads local oversampling problem. To overcome this problem, we introduce the idea of pseudo-depth to the forward modeling. Using the gradient and divergence formula in curvilinear coordinates, we derive the first-order velocity-stress equation in pseudo-depth domain and implement the pseudo-depth domain seismic modeling algorithm. Meanwhile, considering uneven horizontal sampling, we introduce adaptive variable-length spatial operators to calculate the horizontal spatial derivatives and propose seismic forward modeling method with adaptive variable-length spatial operators in the pseudo-depth domain. Finally, by extending the forward operator to reverse time migration, we implement pseudo-depth domain reverse time migration with adaptive variable-length spatial operators. Model test results show that our pseudo-depth domain method with adaptive variable-length spatial operators can reduce about 25% of computer time consuming and nearly 30% of the storage without influencing the accuracy of the results.
Pei Zhenglin,Mu Yongguang.Numerical simulation of seismic wave propagation.Progress in Geophysics, 2004,19(4):933-941.
[3]
张剑锋.弹性波数值模拟的非规则网格差分法.地球物理学报,1998,41(增刊):357-366.
Zhang Jianfeng.Non-orthogonal grid finite-difference method for numerical simulation of elastic wave propagation.Chinese Journal of Geophysic,1998,41(Sl):357-366.
[4]
Mcozo P.Finite-difference technique for SH waves in 2-D media using irregular grids-application to the seismic response problem.Geophys J Int,1989,99(2):321-329.
[5]
Kristek J,Moczo P.Seismic wave propagation in viscoelastic media with material discontinuities-a 3D 4th-order staggered-grid finite-difference modeling.Bull Seism Soc Am,2003,93(5):2273-2280.
Huang Chao,Dong Liangguo.Staggered-grid high-order finite-difference method in elastic wave simulation with variable grids and local time-steps.Chinese Journal of Geophysic,2009,52(11):2870-2878.
[9]
张慧,李振春. 基于双变网格算法的地震波正演模拟.地球物理学报,2011,54(1):77-86.
Zhang Hui,Li Zhenchun.Seismic wave simulation method based on dual-variable grid.Chinese Journal of Geophysic,2011,54(1): 77-86.
[10]
Mi Tieliang,Ma Jianwei,Hervé Chauris et al.Multilevel adaptive mesh modeling for wave propagation in layer media.Journal of Seismic Exploration,2010,19(2):121-139.
[11]
Liu Y and Sen M K.A new time-space domain high-order finite-difference method for the acoustic wave equation.Journal of Computational Physics,2009,228(23):8779-8806.
[12]
Alkhalifah T,Fomel S and Biondi B.The space-time domain:Theory and modelling for anisotropic media. Geophysical Journal International,2001,144(1):105 -113.
[13]
Ma X and Alkhalifah T.Wavefield extrapolation in pseudo-depth domain.73rd EAGE Conference & Exhibition Extended Abstracts, 2011,A014.
[14]
Ma X and Alkhalifah T.Wavefield extrapolation in pseudo-depth domain.74rd EAGE Conference & Exhibition Extended Abstracts, 2012,P212.
[15]
Liu Y, Sen M K.Finite-difference modeling with adaptive variable-length spatial operators.Geophysics, 2011, 76(4):T79-T89.
[16]
Thompson J F,Warsi Z U A,Mastin C W.Numerical Grid Generation:Foundations and Applications. Amsterdam: North-Holland,1985.