1. Shale Gas Research Institute,PetroChina Southwest Oil & Gas field Company,Chengdu 610051,China 2. School of Geosciences,China University of Petroleum(East China),Qingdao 266580,China 3. Chendu Engineering Corporation Limited,Power China,Chengdu 610072,China
The seismic forward modeling technique is critical to seismic exploration.Moreover,it shows a faster rate and higher calculation efficiency in the frequency domain than in the time domain.Presently,there is a need to complete the forward calculation in the frequency domain efficiently and accurately.The specific problems include the numerical dispersion and the high memory consumption for calculating and decomposing impedance,which should be reduced by improving the calculation efficiency.Different from the conventional direct method,this study adopted the least-squares conjugate gradient (LSCG) method used to determine the impedance matrix for the frequency-domain forward modeling and proposed an expression for wavenumber compensation to suppress the numerical dispersion.The numerical tests of simple and complex models show that the LSCG method can effectively reduce the calculation time and that the frequency-domain forward modeling method based on wavenumber compensation can effectively suppress the numerical dispersion and thus improve the precision of wave field simulation.
Cheng K Y. Numerical simulation of forward modeling of wavefield separation of directional traveling waves in anisotropic elastic media[J]. Lithologic Reservoirs, 2014, 26(5):91-96.
[2]
苏巍. 河流相砂岩储层地震响应特征与识别技术研究[D]. 长春: 吉林大学, 2007.
[2]
Su W. A study on seismic response characteristic and identification technology of fluvial sandstone reservoir[D]. Changchun: Jilin University, 2007.
[3]
Lysmer J, Drake L A. A finite element method for seismology[J]. Methods in Computational Physics:Advances in Research and Applications, 1972, 11:181-216.
[4]
Shin C S. Nonlinear elastic wave Inversion by Blocky Parameterization[D]. Tulsa: University of Tulsa, 1988.
[5]
Jo C, Shin C, Suh J. An optimal 9-point finite-difference frequency space 2-D scalar wave extrapolator[J]. Geophysics, 1996, 61(2):529-537.
doi: 10.1190/1.1443979
[6]
Min J, Shin C, Kwon B, et al. Improved frequency-domain elastic wave modeing using weighted-averaging difference operators[J]. Geophysics, 2000, 65(3):884-895.
doi: 10.1190/1.1444785
[7]
Stekl I, Pratt K G. Frequency-domain finite accurate difference visco-elastic modeing by using rotated operators[J]. Geophysics, 1998, 63(5):1779-1794.
doi: 10.1190/1.1444472
[8]
Hustedt B, Opert S, Virieux J. Mixed-grid and staggered grid finite-difference methods for frequency-domain acoustic wave modeling[J]. Geophysics, 2004, 157(3):1269-1296.
[9]
Paige C C, Saunders M A. LSQR:An algorithm for sparse liner equation and sparse least squares[J]. ACM Transaction Mathematical Software, 1982, 8(1):43-71.
doi: 10.1145/355984.355989
[10]
Bjorck A, Duff I S. A direct method for the solution of sparse liner least squares problems[J]. Linear Algebra and its Applications, 1980, 34(1):43-67.
doi: 10.1016/0024-3795(80)90158-5
[11]
张京思, 揣媛媛, 边立恩. 正演模拟技术在渤海油田 X 井区砂体连通性研究中的应用[J]. 岩性油气藏, 2016, 28(3):127-132.
[11]
Zhang J S, Tuan Y Y, Bian L E. Application of forward modeling to study of sand body connectivityin X well field of Bohai Oilfield[J]. Lithologic Reservoirs, 2016, 28(3):127-132.
Feng D S, Wang X Y. Elastic wave propagation simulation in anisotropic media and random media using high-order difference method of rotation staggered grids based on convolutional perfectly matched layer[J]. Geophysical and Geochemical Exploration, 2018, 42( 4) :766-776.
Wang G W, Wang H Y, Li H Q, et al. Application of seismic forward simulation technology in deep reflection seismic profile detection[J]. Geophysical and Geochemical Exploration, 2021, 45( 4) :970-980.
Cao S H, Chen J B. A 17-point Scheme and its numerical implementation for high-accuracy modeling of frequency-domain acoustic equation[J]. Chinese Journal of Geophysics, 2012, 55(10):3440-3449.