| 频率域波动方程正演基于多重网格预条件的迭代算法研究 |
| |
| 引用本文： | 曹健,陈景波,曹书红.频率域波动方程正演基于多重网格预条件的迭代算法研究[J].地球物理学报,2015,58(3):1002-1012. |
| |
| 作者姓名： | 曹健 陈景波 曹书红 |
| |
| 作者单位： | 1. 中国科学院地质与地球物理研究所, 中国科学院油气资源研究重点实验室, 北京 100029;2. 中国科学院大学, 北京 100049 |
| |
| 基金项目： | 国家自然科学基金(41474104,41274134)资助 |
| |
| 摘 要： | 频率域全波形反演是重要的地震成像方法,而频率域波动方程数值模拟是频率域全波形反演的基础.对于大规模的问题,由于受存储和计算量的限制,基于LU分解的直接方法一般不再适用,而是采用迭代方法.基于多重网格预条件的双共轭梯度稳定化方法是一种重要的迭代方法.本文重点讨论了多重网格预条件求解过程中的松弛因子选择方法,研究结果表明,(1)对于一般选取的松弛因子,随模型复杂性的增加,所能计算的重数逐渐下降,方法的实用性也随之下降;(2)对于复杂模型,采用局部模式分析方法选取松弛因子,提高了所能计算的重数,保证了多重网格方法的收敛性和实用性.这些研究成果对基于多重网格预条件的迭代算法的实际应用具有重要意义.
|
| 关 键 词： | 多重网格方法 松弛因子 局部模式分析 平均导数优化格式 |
| 收稿时间： | 2014-07-07 |
Studies on iterative algorithms for modeling of frequency-domain wave equation based on multi-grid precondition |
| |
| CAO Jian;CHEN Jing-Bo;CAO Shu-Hong.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. |
| |
| Authors: | CAO Jian;CHEN Jing-Bo;CAO Shu-Hong |
| |
| Institution: | 1. Key Laboratory of Petroleum Resources Research, Institute of Geology and Geophysics, Chinese Academy of Sciences, Beijing 100029, China;2. University of Chinese Academy of Sciences, Beijing 100049, China |
| |
| Abstract: | 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-conjugate-gradient-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.#br#Numerical experiments reveal:(1) for a generally-chosen relaxation factor, the levels of the multi-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. |
| |
| Keywords: | Multi-grid method Relaxation factor Local mode analysis Average-derivative optimal scheme |
| 本文献已被 CNKI 等数据库收录！ |
| 点击此处可从《地球物理学报》浏览原始摘要信息 |
| 点击此处可从《地球物理学报》下载免费的PDF全文 |
|
