王文化^1,2,,
文晓涛^1,2,,,
张波³,
张懿疆^1,2,
王威⁴
1. 成都理工大学油气藏地质及开发工程国家重点实验室, 成都 610059
2. 成都理工大学地球物理学院, 成都 610059
3. 美国阿拉巴马大学地质科学系, 塔斯卡卢萨 35401
4. 中国石油化工股份有限公司勘探分公司, 成都 610041

基金项目: 国家自然科学基金(41774142), 国家自然科学基金联合基金项目(U1562111), 国家科技重大专项(2016ZX05002-004)联合资助


详细信息

作者简介: 王文化, 男, 1990年生, 博士研究生, 主要从事地震波场数值模拟方面的研究工作.E-mail:wwh_cdut@outlook.com

通讯作者: 文晓涛, 男, 1976年生, 教授、博士生导师, 主要从事储层地球物理综合预测、烃类检测方面的研究工作.E-mail:wenxiaotao@cdut.cn

中图分类号: P631

收稿日期:2019-05-05
修回日期:2019-09-06
上线日期:2020-06-05

Adaptive finite difference forward modeling of two-phase media based on the Remez iterative algorithm for solving differential coefficients

WANG WenHua^1,2,,
WEN XiaoTao^1,2,,,
ZHANG Bo³,
ZHANG YiJiang^1,2,
WANG Wei⁴
1. State Key Laboratory of Oil and Gas Reservoir Geology and Exploitation, Chengdu University of Technology, Chengdu 610059, China
2. College of Geophysics, Chengdu University of Technology, Chengdu 610059, China
3. Department of Geological Sciences, The University of Alabama, Tuscaloosa 35401, USA
4. Exploration Company of SINOPEC, Chengdu 610041, China


More Information

Corresponding author: WEN XiaoTao,E-mail:wenxiaotao@cdut.cn

MSC: P631

--> Received Date: 05 May 2019
Revised Date: 06 September 2019
Available Online: 05 June 2020


摘要
摘要:利用传统有限差分方法对基于Biot理论的双相介质波动方程进行数值求解时，由于慢纵波的存在，数值频散效应较为明显，影响模拟精度.相对于声学近似方程及普通弹性波方程，Biot双相介质波动方程在同等数值求解算法和精度要求条件下，其地震波场正演模拟需要更多的计算时间.本文针对Biot一阶速度-应力方程组发展了一种变阶数优化有限差分数值模拟方法，旨在同时提高其正演模拟的精度和效率.首先结合交错网格差分格式推导Biot方程的数值频散关系式.然后基于Remez迭代算法求取一阶空间偏导数的优化差分系数，并用于Biot方程的交错网格有限差分数值模拟.在此基础上把三类波的平均频散误差参数限制在给定的频散误差阈值和频率范围内，此时优化有限差分算子的长度就能自适应非均匀双相介质模型中的不同速度区间.数值频散曲线分析表明：基于Remez迭代算法的优化有限差分方法相较传统泰勒级数展开方法在大波数范围对频散误差的压制效果更明显；可变阶数的优化有限差分方法能取得与固定阶数优化有限差分方法相近的模拟精度.在均匀介质和河道模型的数值模拟实验中将本文变阶数优化有限差分算法与传统泰勒展开算法、最小二乘优化算法进行比较，进一步证明其在复杂地下介质中的有效性和适用性.
关键词: 正演模拟/
Remez迭代算法/
最小二乘/
双相介质/
自适应有限差分

Abstract:When the finite difference method (FDM) is used to numerically solve the two-phase medium wave equation based on Biot theory, the numerical dispersion effect is obvious due to the existence of slow compressional waves (P-wave), which affects the accuracy of numerical simulation. Compared with the acoustic approximation equation and the ordinary elastic wave equation, the Biot two-phase medium wave equation requires more calculation time for the seismic wave-field forward modeling under the same numerical solution algorithm and accuracy requirements. To solve this problem, this work developed a variable-order optimal FDM (VOFDM) for Biot first-order velocity-stress equations to improve the accuracy and efficiency of the forward modeling. Firstly, the numerical dispersion relation of Biot equation is derived by combining the staggered-grid FD (SFD) scheme. Then, based on the Remez iterative (RI) algorithm, the optimal differential coefficients of the first-order spatial partial derivatives are obtained and used for the SFD of the Biot equation. On this basis, the average dispersion error parameters of three kinds of waves are limited to a given dispersion error threshold and frequency range, so that the length of the FD operator can be adaptive to different velocity intervals in the inhomogeneous two-phase media model. The numerical dispersion curve analysis shows that the optimized FDM based on the RI algorithm is more effective in suppressing the dispersion error in the large wave number range than the traditional Taylor-series expansion (TE) method. The VOFDM can achieve similar simulation accuracy with the fixed-order optimized FDM (FOFDM). The numerical simulation experiments on homogeneous medium and river models prove VOFDM's effectiveness and applicability in complex underground media compared with the TE algorithm and least squares (LS) optimization algorithm.
Key words:Forward modeling/
Remez iterative algorithm/
Least squares/
Two-phase medium/
Adaptive finite difference



PDF全文下载地址:

http://www.geophy.cn/data/article/export-pdf?id=dqwlxb_15495