频率域全波形反演是重要的地震成像方法,而频率域波动方程数值模拟是频率域全波形反演的基础.对于大规模的问题,由于受存储和计算量的限制,基于LU分解的直接方法一般不再适用,而是采用迭代方法.基于多重网格预条件的双共轭梯度稳定化方法是一种重要的迭代方法.本文重点讨论了多重网格预条件求解过程中的松弛因子选择方法,研究结果表明,(1)对于一般选取的松弛因子,随模型复杂性的增加,所能计算的重数逐渐下降,方法的实用性也随之下降;(2)对于复杂模型,采用局部模式分析方法选取松弛因子,提高了所能计算的重数,保证了多重网格方法的收敛性和实用性.这些研究成果对基于多重网格预条件的迭代算法的实际应用具有重要意义.
Frequency-domain full waveform inversion(FWI)is an important method for seismic imaging,and frequency-domain modeling is the basis of frequency-domain FWI.For large-scale problems,the LU-decomposition-based direct method is no longer applicable due to limitations of storage and computational time.Instead,iterative methods are employed.The bi-conjugategradient-stabilized method,which is based on multi-grid precondition,is an important iterative method.The preconditioner based on a heavily damped wave equation is approximately solved with one multi-grid cycle.The multi-grid method is implemented with a weighted Jacobi smoothing,a standard full-weighting coarsening,a linear interpolation,and a matrix-free implementation.To make it convergent for complex model,local model analysis is used to obtain the relaxation factor in the implementation of weighted Jacobi smoothing.Numerical experiments reveal:(1)for a generally-chosen relaxation factor,the levels of themulti-grid decrease as the complexity of the model increases,and,accordingly,the method becomes less practical;(2)for complex models,the relaxation factor obtained by local mode analysis increases the levels and reduces the number of iterations for each single frequency.The bi-conjugate-gradient-stabilized method based on multi-grid precondition obtains its efficiency and precision by using one multi-grid cycle for the approximate inversion of the preconditioner.To obtain reasonably fast convergence of the multi-grid method for complex model,local model analysis is applied in the relaxation factor selection.Compared with a generally-chosen relaxation factor,the relaxation factor obtained by local mode analysis increases the levels and guarantees the convergence and practicality of the multi-grid method.These conclusions are of great significance for application of the multi-grid-precondition-based iterative algorithm.