1. Xi'an Centre of Geological Survey, China Geological Survey, Xi'an, Shaanxi 710054, China; 2. Institute of Geophysics, School of Geological Engineering and Geomatics, Chang'an University, Xi'an, Shaanxi 710054, China; 3. Department of Earth Sciences, Khalifa University of Science and Technology, Abu Dhabi 2533, UAE
Abstract:The 2.5D seismic wavefield numerical stimulation employs the point source in 2D geological models to calculate 3D seismic wavefields. In this paper, we present a generalized 2.5D first-order time-domain wave equation that can be applied to different media (acoustic isotropic, elastic isotropic, and elastic anisotropic) and various boundary conditions (acoustic free-surface, solid free-surface, and solid-liquid boundary). The wave equation is solved by a curvilinear-grid finite-difference method. A comparison of 2.5D numerical solutions, 3D analytic solutions, and 3D numerical solutions in different homogeneous medium models (acoustic isotropic, elastic isotropic, and elastic anisotropic) verifies the correctness of the derived equation and the numerical solution method. It also demonstrates that compared with the 3D nume-rical method, the 2.5D numerical method has great advantages in calculation efficiency and memory footprint. The 2D numerical solutions cannot be applied directly in that they suffer significant amplitude distortion and phase shifts due to an artificial line source applied in this method. The results of numerical experiments show that the proposed 2.5D numerical simulation method can be applied to geological models with different boundary conditions (acoustic free-surface, solid free-surface, and solid-liquid boundary). In addition, unlike the 2D wavefield numerical simulation method, the 2.5D method can be directly employed to process actual point source observation data such as 2.5D reverse-time migration.
基金资助:本项研究受国家重点研发计划“典型覆盖区航空地球物理技术示范与处理解释软件平台开发”所属课题“北秦岭华阳川地区隐伏铀矿空—地—井协同勘查技术示范研究”(2017YFC0602205)和Reward of Khalifa University of Science and Technology (CIRA-2018-48)联合资助。
李青阳, 吴国忱, 段沛然. 准规则网格高阶有限差分法非均质弹性波波场模拟[J]. 石油地球物理勘探, 2019, 54(3):539-550.LI Qingyang, WU Guochen, DUAN Peiran. Elastic wavefield forward modeling in heterogeneous media based on the quasi-regular grid high-order finite difference[J]. Oil Geophysical Prospecting, 2019, 54(3):539-550.
[2]
Zhu J, Dorman J. Two-dimensional, three-component wave propagation in a transversely isotropic medium with arbitrary-orientation-finite-element modelling[J]. Geophysics, 2000, 65(3):934-942.
[3]
Ge Z, Chen X. An efficient approach for simulating wave propagation with the boundary element method in multilayered media with irregular interfaces[J]. Bulletin of the Seismological Society of America, 2008, 98(6):3007-3016.
[4]
谭文卓, 吴帮玉, 李博, 等. 梯形网格伪谱法地震波场模拟[J]. 石油地球物理勘探, 2020, 55(6):1282-1291.TAN Wenzhuo, WU Bangyu, LI Bo, et al. Seismic wave simulation using a trapezoid grid pseudo-spectral method[J]. Oil Geophysical Prospecting, 2020, 55(6):1282-1291.
[5]
张庆朝, 朱国维, 周俊杰, 等. TTI介质qP波伪谱法正演模拟[J]. 石油地球物理勘探, 2019, 54(2):302-311.ZHANG Qingchao, ZHU Guowei, ZHOU Junjie, et al. qP-wave numerical simulation in TTI media with pseudo-spectral method[J]. Oil Geophysical Prospecting, 2019, 54(2):302-311.
[6]
Zhou B, Greenhalgh S A. 3-D frequency-domain seismic wave modelling in heterogeneous, anisotropic media using a Gaussian quadrature grid approach[J]. Geophysical Journal International, 2011, 184(1):507-526.
[7]
Dormy E, Tarantola A. Numerical simulation of elastic wave propagation using a finite volume method[J]. Journal of Geophysical Research:Solid Earth, 1995, 100(B2):2123-2133.
[8]
Zhang Y, Duan L, Xie Y. A stable and practical implementation of least-squares reverse time migration[J]. Geophysics, 2015, 80(1):V23-V31.
[9]
Li Y E, Demanet L. Full-waveform inversion with extrapolated low-frequency data[J]. Geophysics, 2016, 81(6):R339-R348.
[10]
Zhou B, Greenhalgh S A, Hansruedi M. 2.5-D frequency-domain seismic wave modelling in heteroge-neous, anisotropic media using a Gaussian quadrature grid technique[J]. Computer and Geosciences, 2012, 39(1):18-33.
[11]
Auer L, Nuber A M, Greenhalgh S A, et al. A critical appraisal of asymptotic 3D-to-2D data transformation in full-waveform seismic crosshole tomography[J]. Geophysics, 2013, 78(6):R235-R247.
[12]
Williamson P R, Pratt R G. A critical review of the acoustic wave modelling procedure in 2.5 dimensions[J]. Geophysics, 1995, 60(2):591-595.
[13]
Song Z, Williamson P R. Frequency-domain acoustic-wave modeling and inversion of crosshole data:Part I, 2.5-D modeling method[J]. Geophysics, 1995, 60(3):784-795.
[14]
Zhou B, Greenhalgh S A. An adaptive wavenumber sampling strategy for 2.5D seismic-wave modelling in the frequency-domain[J]. Pure and Applied Geophysics, 2006, 163(7):1399-1416.
[15]
Zhou B, Greenhalgh S A. A damping method for the computation of the 2.5-D Green's function for arbitrary acoustic media[J]. Geophysical Journal International, 1998, 133(1):111-120.
[16]
Maokun L, Vladimir D, Aria A, et al. A 2.5D finite-difference algorithm for elastic wave modeling using near-optimal quadratures[J]. Geophysics, 2016, 81(4):T155-T162.
[17]
Takenaka H, Kennett B L N. A 2.5-D time-domain elastodynamic equation for plane-wave incidence[J]. Geophysical Journal International, 1996, 125(2):F5-F9.
[18]
Novais A, Santos L T. 2.5D finite-difference solution of the acoustic wave equation[J]. Geophysical Prospecting, 2005, 53(4):523-531.
[19]
Xiong J L, Lin Y, Abubakar A, et al. 2.5-D forward and inverse modelling of full-waveform elastic seismic survey[J]. Geophysical Journal International, 2013, 193(2):938-948.
[20]
Zhang W, Chen X F. Traction image method for irregular free surface boundaries in finite difference seismic wave simulation[J]. Geophysical Journal International, 2006, 167(1):337-353.
[21]
李芳, 曹思远, 孙建国, 等. 2.5D傅氏变换法声波方程数值模拟及精度分析[J]. 石油地球物理勘探, 2009, 44(3):276-281.LI Fang, CAO Siyuan, SUN Jianguo, et al. Numeric simulation and precision analysis of acoustic equation by 2.5D Fourier transform[J]. Oil Geophysical Prospecting, 2009, 44(3):276-281.
[22]
Aki K, Richards P G. Quantitative Seismolog(2nd Edition)[M]. University Science Books, Sausalito, CA, 2002.
[23]
蒋丽丽, 孙建国. 基于Poisson方程的曲网格生成技术[J]. 世界地质, 2008, 27(3):298-305.JIANG Lili, SUN Jianguo. Source terms of elliptic system in grid generation[J]. Global Geology, 2008, 27(3):298-305.
[24]
Hixon R, Turkel E. Compact implicit MacCormack-type schemes with high accuracy[J]. Journal of Computational Physics, 2000, 158(1):51-70.
[25]
Tam C K, Webb J C. Dispersion-relation-preserving finite difference schemes for computational acoustics[J]. Journal of Computational Physics, 1993, 107(2):262-281.
[26]
Hixon R. Evaluation of a high-accuracy MacCormack-type scheme using benchmark problems[J]. Journal of Computational Acoustics, 1997, 6(3):291-305.
[27]
Zhou B, Moosoo W, Greenhalgh S, et al. Generalized stiffness reduction method to remove the artificial edge-effects for seismic wave modelling in elastic anisotropic media[J]. Geophysical Journal Internatio-nal, 2020, 220(2):1394-1408.
[28]
Vavryuk V. Asymptotic green's function in homoge-neous anisotropic viscoelastic media[J]. Proceedings of the Royal Society A:Mathematical, Physical and Engineering Sciences, 2007, 463(2):2689-2707.