Frequency-domain seismic numerical modeling by applying hybrid conjugate gradient iteration
LIU Wen'ge1, ZHOU Milu1, PENG Haotian2, LIU Fulie1, MOU Qisong1
1. School of Geoscience and Technology, Southwest Petroleum University, Chengdu, Sichuan 610500, China; 2. Exploration and Development Research Institute, PetroChina Southwest Oil & Gasfield Company, Chengdu, Sichuan 610000, China
Abstract：At present, LU decomposition method or the Krylov subspace iteration method is usually used to solve the time-harmonic wave equation in frequency-domain seismic modeling. The direct decomposition method takes up a lot of memory and is time-consuming, and it can hardly simulate high-dimension, large-density seismic acquisition. The gene-ral Krylov subspace iteration method, however, converges slowly and may not converge while dea-ling with complex models. On the basis of the e-xisting conjugate gradient algorithms, an optimized Krylov subspace method-the hybrid conjugate gradient algorithm, is developed to solve time-harmonic wave equations. The numerical simulations of the layered medium model and the standard model indicate that compared with LU decomposition method, the proposed method can effectively reduce the memory requirement and calculation time on the premise of ensuring accuracy. Compared with the stable bi-conjugate gradient iteration method, it has better computational stability in dealing with complex models.
LYSMER J，DRAKE L A. A finite element method for seismology[J]. Methods in Computational Physi-cs Advances in Research and Applications，1972，11(2)：181-216.
[3]
MARFURT K J. Accuracy of finite-difference and finite-element modeling of the scalar and elastic wave equations[J]. Geophysics，1984，49(4)：533-549.
[4]
JO C H，SHIN C，SUH J H. An optimal 9-point，finite difference，frequency-space，2D scalar wave extrapolator[J]. Geophysics，1996，61(2)：529-537.
[5]
SHIN C，SOHN H. A frequency-space 2-D scalar wave extrapolator using extended 25-point finite-difference operator[J]. Geophysics，1998，63(1)：289-296.
[6]
吴国忱，梁锴. VTI介质频率—空间域准P波正演模拟[J]. 石油地球物理勘探，2005，40(5)：535-545.WU Guochen，LIANG Kai. Quasi P-wave forward modeling in frequency-space domain in VTI media[J]. Oil Geophysical Prospecting，2005，40(5)：535-545.
[7]
刘璐，刘洪，刘红伟. 优化15点频率—空间域有限差分正演模拟[J]. 地球物理学报，2013，56(2)：644-652.LIU Lu，LIU Hong，LIU Hongwei. Optimal 15-point finite difference forward modeling in frequency-space domain[J]. Chinese Journal of Geophysics，2013，56(2)：644-652.
[8]
范娜，成景旺，秦雷，等. 一种优化的频率域三维声波有限差分模拟方法[J]. 地球物理学报，2018，61(3)：1095-1108.FAN Na，CHENG Jingwang，QIN Lei，et al. An optimal method for frequency-domain finite-difference solution of 3D scalar wave equation[J]. Chinese Journal of Geophysics，2018，61(3)：1095-1108.
[9]
唐超，文晓涛，王文化. 基于最小范数优化交错网格有限差分系数的波动方程数值模拟[J]. 石油地球物理勘探，2021，56(5)：1039-1047.TANG Chao，WEN Xiaotao，WANG Wenhua. Numerical simulation of wave equations based on minimum-norm optimization of staggered-grid finite difference coefficients[J]. Oil Geophysical Prospecting，2021，56(5)：1039-1047.
[10]
TANG X D，LIU H，ZHANG H，et al. An adaptable 17-point scheme for high-accuracy frequency-domain acoustic wave modeling in 2D constant density media[J]. Geophysics，2015，80(6)：T211-T221.
[11]
XU W H，GAO J H. Adaptive 9-point frequency-domain finite difference scheme for wavefield modeling of 2D acoustic wave equation[J]. Journal of Geophy-sics and Engineering，2018，15(4)：1432-1445.
[12]
殷文. 基于频率域高阶有限差分法的正演模拟及并行算法[J]. 吉林大学学报(地球科学版)，2008，38(1)：144-151.YIN Wen. Forward modeling and parallel algorithm based on high-order finite-difference method in frequency domain[J]. Journal of Jilin University(Earth Science Edition)，2008，38(1)：144-151.
[13]
马召贵，王尚旭，宋建勇. 频率域波动方程正演中的多网格迭代算法[J]. 石油地球物理勘探，2010，45(1)：1-5.MA Zhaogui，WANG Shangxu，SONG Jianyong. Multigrid iterative algorithm in frequency domain wave equation forward modeling[J]. Oil Geophysical Prospecting，2010，45(1)：1-5.
[14]
曹健，陈景波，曹书红. 频率域波动方程正演基于多重网格预条件的迭代算法研究[J]. 地球物理学报，2015，58(3)：1002-1012.CAO Jian，CHEN Jingbo，CAO Shuhong. Studies on iterative algorithms for modeling of frequency-domain wave equation based on multi-grid precondition[J]. Chinese Journal of Geophysics，2015，58(3)：1002-1012.
[15]
DU Z L，LIU J J，LIU W G，et al. Frequency-space domain acoustic wave simulation with the BiCGstab(1) iterative method[J]. Journal of Geophysics and Engineering，2016，13(1)：70-77.
[16]
BELONOSOV M，KOSTIN B，NEKLYUDOV D，et al. 3D numerical simulation of elastic waves with a frequency-domain iterative solver[J]. Geophysics，2018，83(6)：T333-T344.
[17]
熊治涛，唐新功，李丹丹. 二维电性各向异性极化体的频率域响应[J]. 石油物探，2020，59(3)：481-490.XIONG Zhitao，TANG Xingong，LI Dandan. Frequency responses of a two-dimensional anisotropic polarization body[J]. Geophysical Prospecting for Petroleum，2020，59(3)：481-490.
[18]
JIANG W，CHEN X H，JIANG S S，et al. Multiscale finite-difference method for frequency-domain acoustic wave modeling[J]. Geophysics，2021，86(4)：T193-T210.
[19]
JAYSAVAL P，DATTA D，SEN M，et al. A schur complement based fast 2D finite-difference multimodel modeling of acoustic wavefield in the frequency domain[C]. SEG Technical Program Expanded Abstracts，2017，36：4086-4090.
[20]
颜红芹. 基于声波散射理论的零井源距VSP波场数值模拟方法[J]. 石油地球物理勘探，2020，55(3)：567-574.YAN Hongqin. Numerical simulation of zero-offset VSP wave field based on the acoustic scattering theory[J]. Oil Geophysical Prospecting，2020，55(3)：567-574.
[21]
韩利，韩立国，李翔，等. 二阶声波方程频域PML边界条件及频域变网格步长并行计算[J]. 吉林大学学报(地球科学版)，2011，41(4)：1226-1232.HAN Li，HAN Liguo，LI Xiang，et al. PML boundary conditions for second-order acoustic wave equations and variable grid parallel computation in frequency-domain modeling[J]. Journal of Jilin University(Earth Science Edition)，2011，41(4)：1226-1232.
[22]
郭振波，李振春. 起伏地表条件下频率空间域声波介质正演模拟[J]. 吉林大学学报(地球科学版)，2014，44(2)：683-693.GUO Zhenbo，LI Zhenchun. Acoustic wave modeling in frequency-spatial domain with surface topography[J]. Journal of Jilin University(Earth Science Edition)，2014，44(2)：683-693.
[23]
吴悠，吴国忱，李青阳，等. 频率—空间域非均质声波有限差分模拟[J]. 石油地球物理勘探，2022，57(2)：342-356.WU You，WU Guochen，LI Qingyang，et al. A finite difference scheme in frequency-space domain to solve heterogeneous acoustic wave equation[J]. Oil Geophysical Prospecting，2022，57(2)：342-356.
[24]
PLESSIX R E，MULDER W A. Separation-of-variables as a preconditioner for an iterative Helmholtz solver[J]. Applied Numerical Mathematics，2003，44(3)：385-400.
[25]
HESTENES M R，STIEFEL E. Methods of conjugate gradients for solving linear systems[J]. Journal of Research of the National Bureau of Standards，1952，49(6)：409-436.
[26]
DAI Y H，YUAN Y. A nonlinear conjugate gradient method with a strong global convergence property[J]. SIAM Journal on Optimization，1999，10(1)：177-182.
[27]
BABAIE-KAFAKI S，GHANBARI R. A hybridization of the Hestenes-Stiefel and Dai-Yuan conjugate gradient methods based on a least-squares approach[J]. Optimization Methods and Software，2015，30(4)：673-681.
[28]
ZHANG L，ZHOU W J，LI D H. Some descent three-term conjugate gradient methods and their global convergence[J]. Optimization Methods and Software，2007，22(4)：697-711.
