基于预条件共轭梯度法的宽带约束波阻抗反演

发布时间：2024-03-09 05:46

　　以地质信息、测井资料为约束条件,构建宽频带初始波阻抗模型,研究参数置换法建立拟线性矩阵方程,对反演解加入硬约束条件,采用非线性预条件共轭梯度数值算法,改进搜索方向,降低系数矩阵条件数,提高收敛速度,实现从带限的地震数据资料中恢复出宽频带波阻抗的空间分布特征。通过砂泥岩互层复杂速度模型数值反演,在5%的随机噪声干扰下,反演得到速度模型也能够描述出储层和逆断层的位置,由于噪音强度分布不同,造成各道反演得到的绝对速度值存在偏差,但最大偏差值仅为真实值的±5%。进一步验证算法可行有效,反演精度高,具有一定的抗噪能力。

【文章页数】：7 页

【部分图文】：

图1宽带约束波阻抗反演流程图

(3)建立反演目标函数，采用参数置换法对原始波阻抗进行参数变换，得到式(6)近似线性系统;采用更改的预处理共轭梯度法求解，对初始波阻抗模型进行迭代修改，直到求解的模型波阻抗合成的地震记录与原始地震数据最佳匹配为止，进而获得最优的高分辨率波阻抗模型，宽带约束波阻抗反演流程见图1。4....

图2十二层单道理论速度模型及其初始模型

设计一个12层单道理论速度模型及一个6层初始速度模型，相关系数仅为0.785(图2)，采样时间为2ms，密度常数1.0，采样点数235个。选取地震子波为零相位阻尼余弦子波，采样率2ms，子波的长度为40ms，主频为60Hz，地震记录是由地震子波与反射系数相褶积形成，加入....

图3地震记录中加入5%随机噪声反演结果对比

利用fortran语言编制计算机程序进行波阻抗反演，其反演效果如图3所示:由图3可知，在原始地震记录中加入5%的随机噪声，反演结果与理论模型相关系数可达0.997，剖面与绝对波阻抗值均与理论值大体相同;对比采用稀疏约束反褶积法[13-14]，由迭代求出的反射系数存在误差，小的反射....

图4砂泥岩互层岩性油气藏模型

以一个61道砂泥岩互层理论速度模型模拟实际油气储层。该模型大致可分为5层:地表为速度较低的第四系覆盖层;第二层为泥岩层;第三层为楔状砂岩体储层，局部含有油气资源;而第四层泥岩层局部存在逆断层，使得下部第五层含水砂岩向上侵入，各层的岩性及速度值见表1。选取的地震子波为阻尼余弦子波，....

本文编号：3922988