The numerical simulation of visco-acoustic media using the high-order Fourier finite-difference method can reflect the seismic response of high-dip strata with geo-absorption effects more accurately.It can adapt to any lateral speed changes and suppress the dispersion and background noise at large dip angles caused by the finite-difference method.The migration precision of the high-dip strata depends on the determination of the constant coefficient of the difference operator and the calculation of the order.In this study,the gradient descent method was used to optimize the high-order finite-difference correction item in the Fourier finite-difference operator.According to the optimization results of the relative errors and constraint coefficients,the approximation effects of higher-order equations were achieved without increasing the order of the equation and then were expanded to viscoelastic media.Using the designed model,it can be concluded that the proposed method is applicable to the forward simulation of the strong spatial variable speed media with absorption and attenuation effects and has high calculation accuracy and efficiency.Accurate seismic numerical simulation of complex geological structures further confirmed the effectiveness of this method.
Xu K, Sun Z D. Least-squares reverse time migration based on visco-acoustic attenuation compensation[J]. Geophysical Prospecting for Petroleum, 2018, 57(3):419-427.
Liao J P, Wang H Z, Liu H X, et al. Accurate visco-acoustic wave finite difference numerical simulation in frequency space domain[J]. Geophysical and Geochemical Exploration, 2011, 35(4):541-545.
[8]
Stolt R H. Migration by fourier transform[J]. Geophyscis, 2012, 43(1): 23-48.
Li J L, Li Z C, Guan L P, et al. The method of seismic attenuation and energy compensation[J]. Geophysical and Geochemical Exploration, 2015, 39(3):456-465.
Deng W Z, Li Z C, Wang Y G, et al. The least-squares reverse time migration for visco-acoustic medium based on a stable reverse-time propagator[J]. Geophysical and Geochemical Exploration, 2015, 39(4):791-796.
Li J L, Qu Y M, Liu J X, et al. A model study of three-dimensional viscoacoustic least-squares reverse time migration[J]. Geophysical and Geochemical Exploration, 2018, 42(5):1013-1025.
[13]
赵连锋. 井间地震波速与衰减联合层析成像方法研究[D]. 成都: 成都理工大学, 2002.
[13]
Zhao L F. Study on crosswell seismic tomography combing velocity and attenuation[D]. Chengdu: Chengdu University of Technology, 2002.
Zhang J H, Wang W M, Zhao L F, et al. Fourier finite-different forward modeling in viscoacoustic media[J]. Oil Geophysical Prospecting, 2008, 43(2):174-178.
[16]
Stoffa P L, Fokkema J T, de Luna Freire R M, et al. Split-step Fourier migration[J]. Geophysics, 1990, 55(4):410-421.
doi: 10.1190/1.1442850
Wang H Z, Ma Z T, Cao J Z. Three dimensional one-pass migration using paraxial approximate equation with optimized coefficients[J]. Oil Geophysical Prospecting, 1998, 33(2):170-184.
[22]
Liu L N, Zhang J F. 3D wavefield extrapolation with optimum split-step Fourier method[J]. Society of Exploration Geophysicists, 2006, 71(3):95-108.
[23]
Lee M W, Suh S Y. Optimization of one-way wave equations[J]. Geophysics, 1985, 50(10):1634-1637.
doi: 10.1190/1.1441853
[24]
Kadalbajoo M K, Awasthi A. Crank-Nicolson finite difference method based on a midpoint upwind scheme on a non-uniform mesh for time-dependent singularly perturbed convection-diffusion equations[J]. International Journal of Computer Mathematics, 2008, 85(5):771-790.
doi: 10.1080/00207160701459672
[25]
Claerbout J F. Imaging the earth’s interior[J]. Geophysical Journal International, 1986, 86(1):217.
[26]
Li Z M, Liu C L. An ideal depth extrapolation of two-dimensional seismic wave field[J]. Oil Geophysical Prospecting, 1990, 25(5):517-528.