Undulating surface body-fitted grid seismic modeling based on fully staggered-grid mimetic finite difference
Li Qingyang1, Huang Jianping1, Li Zhenchun1, Li Na2, Wang Chengquan3, Zhang Yanyan2
1. School of Geoscience, China University of Petroleum (East China), Qingdao, Shandong 266580, China;
2. Geophysical Research Institute, Zhongyuan Oilfield Branch Co., SINOPEC, Puyang, Henan 457001, China;
3. Geophysical Research Institule, Huabei Oilfield company, PetroChina, Renqiu, Hebei 062552, China
Abstract：In order to accurately simulate wave field characteristics for undulating surface, we present a novel finite-difference method. The method is based upon the fully staggered grid combined with body-fitted grid to honor undulating surface and a mimetic solution for the free-surface condition. This algorithm has three main advantages: ①The use of body-fitted grid to discrete topography can not only eliminate the ladder of discrete problems, but also adapt to very complex topography. ②Introducing the fully-staggered grid to curvilinear coordinates, this method can avoid not only interpolation error when standard staggered grid used, but also achieve high-frequency oscillations with the collocated-grid. Therefore, the method improves the simulation accuracy and reduces the complexity of the algorithm degrees. ③The mimetic finite-difference has the same differential accuracy in boundary and internal to implement free boundary conditions in curvilinear coordinates. Thereby this method improves the simulation accuracy further and has a more reliable theoretical basis compared with traction image method. Model tests confirm that the proposed method matches almost to the analytical solution with a high precision and adapt to various complex models. Owing to a relaxed dispersion condition, the proposed method can obtain results with very few points per wavelength when compared to other similar approaches.
Zhang Hua, Li Zhenchun, Han Wengong.Review of seismic wave numerical simulation from irregular topography. Progress in Exploration Geophysics,2007,30(5): 334-339.
[2]
孙建国.复杂地表条件下地球物理场数值模拟方法评述.世界地质,2007,26(3):345-362] Sun Jianguo. Methods for numerical modeling of geophysical fields under complex topographical conditions: a critical review.Global Geology,2007,26(3):345-362.
Wang Xueqiu,Sun Jianguo.The state-of-the-art in numerical modeling including surface topography with finite-difference method.Progress in Geophysics, 2008, 23(1): 40-48.
[5]
Virieux J. P-SV wave propagation in heterogeneous media: velocity-stress finite-difference method. Geophysics,1986, 51(4):889-901.
Xue Dongchuan, Wang Shangxu, Jiao Shujing. Wave equation finite-element modeling including rugged topography and complicated medium. Progress in Geophysics, 2007, 22(2):522-529.
[7]
Moczo P, Bystricky E, Kristek J et al. Hybrid modelling of P-SV seismic motion at inhomogeneous viscoelastic topographic structures. Bulletin of the Seismological Society of America,1997,87(5):1305-1323.
[8]
Komatitsch D and Tromp J.Spectral-element simulations of global seismic wave propagation I:Validation.Geophysical Journal International,2002,149(2):390-412.
[9]
Tessmer E and Kosloff D.3-D elastic modeling with surface topography by a Chebychev spectral method. Geophysics, 1994, 59(3): 464-473.
[10]
Ge Z X,Wang C Y,Lei T et al. 2.5D boundary element simulation of wave propagation in a vertical fluid-filled borehole with irregular-shape.Applied Geophysics, 2007, 4(3):1-9.
Wang Xiuming, Zhang Hailan. Modeling the seismic wave in the media with irregular free interface by the finite-difference method. Science in China (Series G), 2004, 34(5): 481-493.
[12]
Hayashi K, Bums D and Toksoz M. Discontinuous-grid finite-difference seismic modeling including surface topography. Bulletin of the Seismological Society of America, 2001, 91(6):1750-1764.
[13]
Hestholm S and Ruud B.2-D finite-difference elastic wave modeling including surface topography.Geophysical Prospecting, 1994, 42(5): 371-390.
[14]
董良国,郭晓玲,吴晓丰等.起伏地表弹性波传播有限差分法数值模拟.天然气工业,2007, 27(10):38-41] Dong Liangguo, Guo Xiaoling, Wu Xiaofeng et al. Finite difference numerical simulation for the elastic wave propagation in rugged topography.Natural Gas Industry, 2007, 27(10): 38-41.
[15]
Tarrass I,Giraud L,Thore P. New curvilinear scheme for elastic wave propagating in presence of curved topography.Geophysical Prospecting,2011,59(5):889-906.
Chu Chunlei,Wang Xiutian. Seismic modeling with a finite-difference method on irregular triangular grids. Periodical of Ocean University of China(Natural Science Edition),2005,35(1): 43-48.
Shi Ruiqi,Wang Shangxu, Guo Rui et al. Finite-element numerical modeling of elastic wave on unstructured meshes. OGP, 2013, 48(6): 915-923.
[18]
Fornberg B. The pseudo-spectral method: Accurate representation in elastic wave calculations. Geophy-sics, 1988, 53(5): 625-637.
[19]
Komatitsch D, Coute F, Mora P. Tensorial formulation of the wave equation for modelling curved interfaces.Geophysical Journal International,1996,127(1):156-168.
[20]
Zhang W and Chen X F. Traction image method for irregular free surface boundaries in finite difference seismic wave simulation. Geophysical Journal International, 2006, 167(1): 337-353.
[21]
Zhang W, Shen Y, Zhao L.Three-dimensional anisotropic seismic wave modelling in spherical coordinates by a collocated-grid finite-difference method. Geophysical Journal International, 2012,188(3): 1359-1381.
Zhu Hejun, Zhang Wei, Chen Xiaofei. Two-dimensional seismic wave simulation in anisotropic media by non-staggered finite-difference method. Chinese Journal of Geophysics, 2009, 52(6):1536-1546.
[23]
Appelo D and Petersson N A.A stable finite difference method for the elastic wave equation on complex geometries with free surfaces.communications in computational Physics,2009,5(1): 84-107.
Lan Hanqiang,Liu Jia and Bai Zhiming.Wave-field simulation in VTI media with irregular free surface. Chinese Journal of Geophysics, 2011,54(8):2072-2084.
Tang Wen,Wang Shangxu,Yuan Sanyi.Stability analysis of differential strategy for rugged topography by second-order elastic wave equation based on coordinate transformation. GPP, 2013,52(5): 457-463.
Qiu Lei,Tian Gang,Shi Zhanjie et al.Finite-difference method for seismic wave numerical simulation in presence of topography—In generally orthogonal curvilinear coordinate system. Journal of Zhejiang University (Engineering Science), 2012, 46(10):1923-1930.
Sun Jianguo,Jiang Lili.Orthogonal curvilinear grid generation technique used for numeric simulation of geophysical fields in undulating surface condition. OGP, 2009, 44(4): 494-500.
Jiang Lili,Sun Jianguo.Source terms of elliptic system in grid generation.Global Geology,2008, 27(3):298-305.
[29]
Lisitsa V,Vishnevskiy D. Lebedev scheme for the numerical simulation of wave propagation in 3D anisotropic elasticity double dagger. Geophysical Prospecting, 2010, 58(4): 619-635.
[30]
Castillo J, Hyman J, Shashkov M et al. Fourth- and sixth-order conservative finite difference approximations of the divergence and gradient. Applied Numerical Mathematics, 2001, 37(1-2):171-187.
[31]
Rojas O, Day S, Castillo J et al.Modelling of rupture propagation using high-order mimetic finite differences.Geophysical Journal International,2008,172(2):631-650.
