流体体积模量参数提取方法技术

本申请提供一种流体体积模量参数提取方法

【技术实现步骤摘要】
流体体积模量参数提取方法、流体识别方法、设备及介质


[0001]本申请涉及油气地球物理勘探
，具体而言，涉及一种流体体积模量参数提取方法
、
流体识别方法
、
设备及介质
。

技术介绍

[0002]弹性阻抗反演最早在
1999
年提出来的，是一类基于
Zoeppritz
方程及其近似表达式的叠前反演方法，结合了不同入射角的地震资料和测井资料，可以定量求取多种流体敏感弹性参数，减少了误差的累计，反演精度显著提高，逐渐成为地球物理储层预测的主流技术
。
[0003]随着油田精细高效开发地不断深入，薄层
、
薄互层成为勘探重点，对反演技术提出了愈来愈高的要求，常规反演方法的分辨率较低，无法满足流体精细刻画的要求
。
[0004]同时，现有技术中提取的常规流体敏感因子数据体，假象明显，导致流体识别的可靠性下降
。

技术实现思路

[0005]本申请实施例的目的在于提供一种流体体积模量参数提取方法
、
流体识别方法
、
设备及介质，其能够提高弹性阻抗数据体分辨率，并且降低流体识别假象，提高流体反演的钻遇率
。
[0006]第一方面，提供了一种流体体积模量参数提取方法，包括以下步骤：
[0007]S1、
准备目的工区的地震资料和测井资料；
[0008]S2、
确定出地震资料对应的有效的入射角范围，将地震道集对应的有效入射角范围由小到大划分为至少三个角度区间，并在每个角度区间中选定位于区间中的一典型角度；
[0009]S3、
将每个角度区间中的所有地震道集各自进行叠加，得到每个角度区间对应的部分叠加体；
[0010]S4、
利用测井资料确定出流体体积模量曲线
、
固体刚性参数曲线
、
密度曲线和每个典型角度对应弹性阻抗曲线；
[0011]S5、
以测井资料得到的弹性阻抗曲线作约束，结合每个角度区间的部分叠加体，进行弹性阻抗反演，得到每个角度区间的典型角度对应的弹性阻抗数据体，并由测井位置对应的弹性阻抗数据体抽取出井旁道弹性阻抗数据；
[0012]S6、
以流体体积模量曲线
、
固体刚性参数曲线
、
密度曲线和每个典型角度对应的弹性阻抗曲线构建计算关系，并将流体体积模量曲线
、
固体刚性参数曲线
、
密度曲线和井旁道弹性阻抗数据带入前述计算关系中，确定出流体体积模量
K
f
、
固体刚性参数
f
m
和密度
ρ
各自的数据体
。
[0013]在一种可实施的方案中，所述步骤
S5
包括以下步骤：
[0014]S51、
对每个典型角度对应的弹性阻抗曲线的小波域进行不同频段多尺度滤波，并
针对弹性阻抗曲线的中低频部分建立起一个超出地震频带范围的高频初始模型；
[0015]S52、
利用部分叠加体求取相对阻抗，结合每个典型角度对应的弹性阻抗曲线，求取绝对阻抗，建立似然函数；
[0016]S53、
在贝叶斯理论框架下
,
并利用马尔可夫链蒙特卡罗随机模拟对高频初始模型进行随机扰动，反复迭代干扰，直至得到的反演结果同时符合中频地震信息和弹性阻抗曲线的结构特征，并计算出与每个典型角度对应的且频带范围超出地震频带范围的弹性阻抗数据，并由测井位置对应的弹性阻抗数据体抽取出井旁道弹性阻抗数据体
。
[0017]在一种可实施的方案中，所述步骤
S6
包括以下步骤：
[0018]S61、
以标准化弹性阻抗方程为基础，并以流体体积模量曲线
、
固体刚性参数曲线
、
密度曲线和每个典型角度对应的弹性阻抗曲线构建三项标准化弹性阻抗方程，并且方程中包含有分别与流体体积模量曲线
、
固体刚性参数曲线
、
密度曲线相关的三个方程系数；
[0019]S62、
将步骤
S61
中的弹性阻抗方程两边取自然对数形成对数方程，然后将每个部分叠加体的数据带入对数方程，形成每个部分叠加体对应的矩阵；
[0020]S63、
将流体体积模量曲线
、
固体刚性参数曲线
、
密度曲线和井旁道弹性阻抗数据带入步骤
S52
中对应的矩阵中，分别计算出每个典型角度对应的不同的方程系数；
[0021]S64、
再将每个典型角度对应的方程系数带入步骤
S62
的对数方程中，形成计算流体体积模量
K
f
、
固体刚性参数
f
m
和密度
ρ
的方程组，求解方程组得到流体体积模量
K
f
、
固体刚性参数
f
m
和密度
ρ
的数据体
。
[0022]在一种可实施的方案中，其特征在于，所述角度区间的数量被设置为大于等于三个
。
[0023]根据本申请的第二方面，还提供了一种流体识别方法，包括前述方案中的流体体积模量参数提取方法，还包括：以提取到的流体体积模量
K
f
的数据体进行流体识别
。
[0024]根据本申请的第三方面，还提供了一种计算机存储介质，其存储有计算机程序，该程序被处理器执行时实现前述方案中的流体体积模量参数提取方法
。
[0025]根据本申请的第四方面，还提供了一种计算机设备，包括存储器及处理器，所述存储器存储有计算机程序，该程序被处理器执行时实现前述方案中的流体体积模量参数提取方法
。
[0026]与现有技术相比，本申请的有益效果至少包括：
[0027]本实施例的方案中，由于利用测井资料得到的弹性阻抗曲线作约束，并结合每个角度区间的部分叠加体进而得到不同角度区间中的典型角度对应的弹性阻抗数据体，有效提高了弹性阻抗数据体的分辨率
。
同时，由于流体体积模量
K
f
不受孔隙度的影响，因此，利用本实施例流体体积模量参数提取方法得到的流体体积模量
K
f
进行流体识别时，可以降低假象，尤其提高了对薄层
、
薄互层的识别准确度和可靠性
。
附图说明
[0028]为了更清楚地说明本申请实施例的技术方案，下面将对实施例中所需要使用的附图作简单地介绍，应当理解，以下附图仅示出了本申请的某些实施例，因此不应被看作是对范围的限定，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.
一种流体体积模量参数提取方法，其特征在于，包括以下步骤：
S1、
准备目的工区的地震资料和测井资料；
S2、
确定出地震资料对应的有效入射角范围，将地震道集对应的有效入射角范围由小到大划分为至少三个角度区间，并在每个角度区间中选定一典型角度；
S3、
将每个角度区间中的所有地震道集各自进行叠加，得到每个角度区间对应的部分叠加体；
S4、
利用测井资料确定出流体体积模量曲线
、
固体刚性参数曲线
、
密度曲线和每个典型角度对应弹性阻抗曲线；
S5、
以测井资料得到的弹性阻抗曲线作约束，结合每个角度区间的部分叠加体，进行弹性阻抗反演，得到每个角度区间的典型角度对应的弹性阻抗数据体，并由测井位置对应的弹性阻抗数据体抽取出井旁道弹性阻抗数据；
S6、
以流体体积模量曲线
、
固体刚性参数曲线
、
密度曲线和每个典型角度对应的弹性阻抗曲线构建计算关系，并将流体体积模量曲线
、
固体刚性参数曲线
、
密度曲线和井旁道弹性阻抗数据带入前述计算关系中，确定出流体体积模量
K
f
、
固体刚性参数
f
m
和密度
ρ
各自的数据体
。2.
根据权利要求1所述的流体体积模量参数提取方法，其特征在于，所述步骤
S5
包括以下步骤：
S51、
对每个典型角度对应的弹性阻抗曲线的小波域进行不同频段多尺度滤波，并针对弹性阻抗曲线的中低频部分建立起一个超出地震频带范围的高频初始模型；
S52、
利用部分叠加体求取相对阻抗，结合每个典型角度对应的弹性阻抗曲线，求取绝对阻抗，建立似然函数；
S53、
在贝叶斯理论框架下
,
并利用马尔可夫链蒙特卡罗随机模拟对高频初始模型进行随机扰动，反复迭代干扰，直至得到的反演结果同时符合中频地震信息和弹性阻抗曲线的结构特征，并计算出与每个典型角度对应的且频带范围超出地震频带范围的弹性阻抗数据体，并由测井位置对应的弹性阻抗数据体抽取出井旁道弹性阻抗数据
。3.
根据...

【专利技术属性】
技术研发人员：王慧欣，张尚虎，刘苗，蒋涔，晏玉环，张大海，
申请(专利权)人：中国石油化工股份有限公司上海海洋油气分公司，
类型：发明
国别省市：