一种基于分频目标迭代体的调谐振幅反演方法技术

本发明专利技术公开了一种基于分频目标迭代体的调谐振幅反演方法，属于石油地震勘探开发技术领域。首先对地震数据进行频带预处理，保证地震资料具有满足目标储层厚度要求的频带宽度；其次通过广义S算法将原始地震数据体分成一系列调谐振幅体；再次利用GR曲线作为目标曲线对每个单频体进行神经网络反演，进一步强化井震匹配效果；最后利用加权调谐振幅反演方法重组优选的数据体，并重复上述两步骤直到得到理想的井震匹配效果的数据体。利用该方法能有效改善厚、薄储层井震不匹配问题，非常有利于厚、薄储层的精细刻画。储层的精细刻画。储层的精细刻画。

【技术实现步骤摘要】
一种基于分频目标迭代体的调谐振幅反演方法


[0001]本专利技术属于石油勘探地震数据处理与解释的
，尤其涉及一种基于分频目标迭代体的调谐振幅反演方法。

技术介绍

[0002]厚、薄储层超过地震分辨率，在地震剖面上表现为多个反射轴或无响应，影响储层刻画。改善厚储层的显示效果，增强井震匹配程度一直是地球物理研究的难点问题之一。随着油气勘探开发工作的持续推进，厚、薄储层的准确刻画倍受重视：勘探阶段，厚、薄储层的准确识别与精细刻画，有利于快速研究圈闭潜力储量，为探井决策提供依据。开发阶段，厚、薄储层的准确识别与精细刻画，对于准确研究砂体沉积微相，剩余油分布及调整井实施有巨大帮助。因此厚、薄储层的准确刻画具有重要意义。
[0003]井震结合的反演数据体能有效改善厚、薄储层显示效果，利于储层平面刻画。现有的波阻抗反演方法在很大程度上改善了厚、薄储层刻画，提高了井震匹配效果，但井间外推存在一定假象；阻抗剖面与子波褶积过程中，受子波分辨率影响，超分辨率储层仍然无法有效刻画，为厚、薄储层精细研究带来困难。

技术实现思路

[0004]本专利技术的目的在于提供一种基于分频目标迭代体的调谐振幅反演方法，以解决
技术介绍
的问题。
[0005]为实现上述目的，本专利技术的一种基于分频目标迭代体的调谐振幅反演方法的具体技术方案如下：
[0006]一种基于分频目标迭代体的调谐振幅反演方法，包括以下步骤：
[0007]分析地震资料频带宽度，若地震资料频宽不满足目标储层厚度要求，则对地震资料的频宽做拓宽频带宽度处理；
[0008]利用匹配追踪算法对地震资料进行分频，得到系列连续离散频率体；
[0009]分析目标储层厚度范围，根据目标储层厚度范围，利用频率与厚度的关系式计算出对应的频率范围，从系列连续离散频率体中，选出合适的调谐频率数据体组合；
[0010]选出合适的调谐频率数据体组合作为输入，以GR曲线为目标曲线，利用神经网络算法进行目标反演，得到对应的系列伽马数据体；
[0011]所得到的GR数据体进行加权调谐振幅反演，得到与井匹配效果理想的数据体。
[0012]进一步的，对地震资料频宽做处理方法选用反Q滤波方法。
[0013]进一步的，若得到的数据体与井匹配效果不理想，修改加权参数重新反演，重复往下运算，经过迭代反演，直到得到理想的数据体。
[0014]进一步的，若得到的数据体与井匹配效果不理想，修改所述合适的调谐频率数据体组合，重复往下运算，经过迭代反演，直到得到理想的数据体。
[0015]进一步的，所述利用匹配追踪算法对地震资料进行分频的计算方法为：
[0016][0017]进一步的，所述神经网络算法进行目标反演的计算方法为：
[0018]信号向前传播算法：
[0019]假设训练样本为X＝(X1,X2,
…
X
n
)
[0020]隐含层第i个节点值为net
i
:
[0021][0022]隐含层第i个节点的输出yi
[0023][0024]输出层第k个节点的输入net
k
为
[0025][0026]其输为出o
k
为
[0027][0028]信号向后传播以求取权重参数，定义误差函数为E
[0029][0030]反向计算隐藏层，误差为
[0031][0032]反向计算输入层，误差为
[0033][0034]根据梯度下降法，校正权值参数
[0035][0036][0037]进一步的，所述加权调谐振幅反演，所用计算方法如下：
[0038][0039]式中s
i
为被选中合适的调谐频率数据体组合经过下一步骤中的神经网络算法处理后得到GR体集，a
i
为赋予的权重系数；引入调整因子F∈[-2，-1]∪[1，2]，该参数用于对权重系数a
i
的调整；
[0040]当F∈[1，2]时，有利于刻画薄层砂体或厚层中受隔夹层影响的砂体；
[0041]当F∈[-2，-1]时，有利于刻画厚层砂体；
[0042]a
i
的表达式为：
[0043][0044][0045]相比较现有技术而言，本专利技术具有以下有益效果：
[0046]1、本专利技术根据反演砂体厚度需要，灵活控制参与运算数据体权重参数的算法。
[0047]2、该算法最大的优势就是根据反演结果的反馈，可直接修改调整因子F，从而控制数据体权重，使最终反演数据体达到理想效果。
[0048]3、该算法直接有效，实现简单，并且每一次迭代结果直观可视，并直接作为下一次参数调整的依据。
附图说明
[0049]图1为本专利技术的流程图；
[0050]图2为原始地震剖面一；
[0051]图3为对应图2的本专利技术求取的地震剖面一；
[0052]图4为原始地震剖面二；
[0053]图5为对应图4的本专利技术求取的地震剖面二。
具体实施方式
[0054]为了更好地了解本专利技术的目的、结构及功能，下面结合附图1-5，对本专利技术的理解。
[0055]如图1，具体实施方式包括以下步骤：
[0056](1)分析地震资料频带宽度，若地震资料频宽不满足目标储层厚度要求，需要对地震资料频宽做处理，如反Q滤波等方法，拓宽频带宽度；
[0057](2)利用广义S算法对地震资料进行分频，得到一系列连续离散频率体；
[0058](3)分析目标储层厚度范围，根据目标储层厚度范围，利用频率与厚度的关系式计算出对应的频率范围，在步骤(2)中优选合适的调谐频率数据体组合；
[0059](4)用步骤(3)优选的每个数据体作为输入，以GR曲线为目标曲线，利用神经网络算法进行目标反演，得到对应的系列伽马数据体；
[0060](5)利用步骤(4)中所得到的GR数据体进行加权调谐振幅反演，得到新的数据体，若得到的数据体与井匹配效果不理想，可以在步骤(4)中修改加权参数重新反演，或者修改
步骤(3)中所选的数据体组合，依步骤往下运算，直到得到理想的数据体。
[0061]步骤(2)中所用计算方法如下：
[0062][0063]步骤(4)中所用BP神经网络计算方法如下：
[0064]信号向前传播算法：
[0065]假设训练样本为X＝(X1,X2,
…
X
n
)
[0066]隐含层第i个节点值为net
i
:
[0067][0068]隐含层第i个节点的输出yi
[0069][0070]输出层第k个节点的输入net
k
为
[0071][0072]其输为出o
k
为
[0073][0074]信号向后传播以求取权重参数，定义误差函数为E
[0075][00本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于分频目标迭代体的调谐振幅反演方法，其特征在于，包括以下步骤：分析地震资料频带宽度，若地震资料频宽不满足目标储层厚度要求，则对地震资料的频宽做拓宽频带宽度处理；利用匹配追踪算法对地震资料进行分频，得到系列连续离散频率体；分析目标储层厚度范围，根据目标储层厚度范围，利用频率与厚度的关系式计算出对应的频率范围，从系列连续离散频率体中，选出合适的调谐频率数据体组合；选出合适的调谐频率数据体组合作为输入，以GR曲线为目标曲线，利用神经网络算法进行目标反演，得到对应的系列伽马数据体；所得到的GR数据体进行加权调谐振幅反演，得到与井匹配效果理想的数据体。2.根据权利要求1所述的基于分频目标迭代体的调谐振幅反演方法，其特征在于，对地震资料频宽做处理方法选用反Q滤波方法。3.根据权利要求1所述的基于分频目标迭代体的调谐振幅反演方法，其特征在于，若得到的数据体与井匹配效果不理想，修改加权参数重新反演，重复往下运算，经过迭代反演，直到得到理想的数据体。4.根据权利要求1所述的基于分频目标迭代体的调谐振幅反演方法，其特征在于，若得到的数据体与井匹配效果不理想，修改所述合适的调谐频率数据体组合，重复往下运算，经过迭代反演，直到得到理想的数据体。5.根据权利要求1所述的基于分频目标迭代体的调谐振幅反演方法，其特征在于，所述利用匹配追踪算法对地震资料进行...

【专利技术属性】
技术研发人员：马佳国，明君，夏同星，王建立，周卿，李文滨，郄莹，于茜，王腾，王波，
申请(专利权)人：中海石油中国有限公司天津分公司，
类型：发明
国别省市：