深地油气精准导航随钻伽马测井曲线调优方法与系统技术方案

本发明专利技术属于地质勘探领域，具体涉及了一种深地油气精准导航随钻伽马测井曲线调优方法与系统，旨在解决的现有的钻进技术中对于钻井的深度的测量精确性不足的问题

【技术实现步骤摘要】
深地油气精准导航随钻伽马测井曲线调优方法与系统


[0001]本专利技术属于地质勘探，具体涉及了一种深地油气精准导航随钻伽马测井曲线调优方法与系统
。

技术介绍

[0002]深地精准导航技术可用于深部油气勘探，指导钻井人员优化井轨迹到油气藏指定位置以获得最大的泄油面积和最佳的采收率，成为提高单井油气产量和油田开发效益的前沿技术
。
在实时钻进过程中，由于仪器
、
泥浆钻井液等非地层因素对测井曲线产生不可忽视的影响，干扰了后续对地层岩性反演的工作
。
因此，需要一种适合于随钻测井资料的调优算法，用于消除不同井次测井曲线的量纲差异
、
非地层因素的噪声影响
、
测量深度误差，提升储层随钻预测模型的准确性与效率
。

技术实现思路

[0003]为了解决现有技术中的上述问题，即现有的钻进技术中对于钻井的深度的测量精确性不足的问题，本专利技术提供了一种深地油气精准导航随钻伽马测井曲线调优方法，所述方法包括：步骤
S100
，实时获取随钻
GR
曲线数据并进行异常值处理获得剔除异常值的
GR
曲线数据；所述随钻
GR
曲线数据为以
batch
的形式获得；获取已钻井的历史
GR
曲线数据；步骤
S200
，对当前
batch
的剔除异常值的
GR
曲线数据根据所述历史
GR
曲线数据进行拉伸
、
平移变换，获得标准化后的随钻
GR
参数数据；步骤
S300
，判断所述标准化后的随钻
GR
参数数据和历史
GR
曲线数据的合理性，得到合理性判断结果；若所述合理性判断结果为通过，则进入步骤
S400
；若所述合理性判断结果为不通过则进入步骤
S100
并调整异常值处理的判定条件并重新进行合理性判断；步骤
S400
，基于已钻井的历史
GR
曲线数据进行波形指示反演，获得已钻井的波形指示反演数据体；步骤
S500
，基于所述已钻井的波形指示反演数据体和预先设计的钻井轨迹，获取钻前预测井曲线；步骤
S600
，基于所述标准化后的随钻
GR
参数数据，进行方差属性特征提取，获得实时更新的随钻方差属性曲线；基于所述钻前预测井曲线，进行方差属性特征提取，获得钻前预测井曲线的方差属性曲线；步骤
S700
，基于所述实时更新的随钻方差属性曲线进行求导获取实时更新的地层突变点检测结果；基于所述钻前预测井曲线的方差属性曲线进行求导获取钻前预测井曲线的地层突变点检测结果；
步骤
S800
，基于所述实时更新的地层突变点检测结果和钻前预测井曲线的地层突变点检测结果，获取设定的4个控制点，基于所述设定的4个控制点进行深度变换获得随钻精确深度；步骤
S900
，基于所述随钻精确深度完成当前
batch
的精确深度矫正，回到步骤
S100
基于下一
batch
的随钻
GR
曲线数据进行随钻精确深度的确定
。
[0004]在一些优选的实施方式中，所述异常值处理，其方法为：预设以样本点为圆心的判定半径
r
和判别比例
s
；从当前
batch
的随钻
GR
曲线数据进行无放回地抽取出第一样本点
p
，第一样本点
p
，通过异常值判定公式进行判定，所述异常值判定公式为：；；表示在预设的判定半径
r
和判别比例
s
的情况下第一样本点
p
为异常值的判定值，表示计数函数，表示一个中的计数算子，表示在一个中无放回地抽取样本点
q
时判断样本点
q
为样本点
p
的临近点的个数，指样本点
p
与样本点
q
之间的欧几里得距离；当得值大于预设得阈值时判定为异常值，剔除对应得第一样本点
。
[0005]在一些优选的实施方式中，所述步骤
S200
，具体包括：步骤
S210
，计算当前
batch
的剔除异常值的
GR
曲线数据的均值估计值和方差估计值：；；表示一个
batch
中的计数算子，
j
表示样本序号，
i
表示样本序号，
N
与的值相同；步骤
S220
，计算历史
GR
曲线数据的标准差，对当前
batch
的剔除异常值的
GR
曲线数据进行拉伸压缩变换：；表示拉伸压缩变换后得到的
GR
数据，表示当前
batch
的剔除异常值的
GR
曲线数据的第
k
个样本点的
GR
值；步骤
S230
，对所述拉伸压缩变换后得到的
GR
数据，基于历史
GR
曲线数据的均值标
准值进行平移变换：；表示平移变换后得到的标准化后的随钻
GR
参数数据
。
[0006]在一些优选的实施方式中，所述步骤
S300
，具体包括：步骤
S310
，分别计算变换后的随钻
GR
参数数据的概率分布曲线和已钻井的
GR
参数数据的概率分布曲线；步骤
S320
，基于和计算概率分布曲线面积交集和概率分布曲线面积并集；步骤
S330
，基于所述和计算质量控制参数；步骤
S340
，比对所述质量控制参数
U
与设定的质量控制阈值的大小；当所述质量控制参数
U
大于或等于设定的质量控制阈值时，合理性判断结果为通过，进入步骤
S400
；当所述质量控制参数
U
小于设定的质量控制阈值时，合理性判断结果为不通过，进入步骤
S100
并调整圆心的判定半径
r
和判别比例
s。
[0007]在一些优选的实施方式中，所述步骤
S400
，具体包括：利用白噪声满足高斯分布的规律，将已钻井的历史
GR
曲线数据表示为：；其中，表示已钻井的历史
GR
曲线数据，表示待求解的已钻井的波形指示反演数据体，表示随机噪声；根据中心极限定理，也满足高斯分布，确定初始目标函数为：；其中，表示与后验信息有关的函数，表示已钻井的历本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.
一种深地油气精准导航随钻伽马测井曲线调优方法，其特征在于，所述方法包括：步骤
S100
，实时获取随钻
GR
曲线数据并进行异常值处理获得剔除异常值的
GR
曲线数据；所述随钻
GR
曲线数据为以
batch
的形式获得；获取已钻井的历史
GR
曲线数据；步骤
S200
，对当前
batch
的剔除异常值的
GR
曲线数据根据所述历史
GR
曲线数据进行拉伸
、
平移变换，获得标准化后的随钻
GR
参数数据；步骤
S300
，判断所述标准化后的随钻
GR
参数数据和历史
GR
曲线数据的合理性，得到合理性判断结果；若所述合理性判断结果为通过，则进入步骤
S400
；若所述合理性判断结果为不通过则进入步骤
S100
并调整异常值处理的判定条件并重新进行合理性判断；步骤
S400
，基于已钻井的历史
GR
曲线数据进行波形指示反演，获得已钻井的波形指示反演数据体；步骤
S500
，基于所述已钻井的波形指示反演数据体和预先设计的钻井轨迹，获取钻前预测井曲线；步骤
S600
，基于所述标准化后的随钻
GR
参数数据，进行方差属性特征提取，获得实时更新的随钻方差属性曲线；基于所述钻前预测井曲线，进行方差属性特征提取，获得钻前预测井曲线的方差属性曲线；步骤
S700
，基于所述实时更新的随钻方差属性曲线进行求导获取实时更新的地层突变点检测结果；基于所述钻前预测井曲线的方差属性曲线进行求导获取钻前预测井曲线的地层突变点检测结果；步骤
S800
，基于所述实时更新的地层突变点检测结果和钻前预测井曲线的地层突变点检测结果，获取设定的4个控制点，基于所述设定的4个控制点进行深度变换获得随钻精确深度；步骤
S900
，基于所述随钻精确深度完成当前
batch
的精确深度矫正，回到步骤
S100
基于下一
batch
的随钻
GR
曲线数据进行随钻精确深度的确定
。2.
根据权利要求1所述的深地油气精准导航随钻伽马测井曲线调优方法，其特征在于，所述异常值处理，其方法为：预设以样本点为圆心的判定半径
r
和判别比例
s
；从当前的随钻
GR
曲线数据进行无放回地抽取出第一样本点
p
，第一样本点
p
，通过异常值判定公式进行判定，所述异常值判定公式为：；；表示在预设的判定半径
r
和判别比例
s
的情况下第一样本点
p
为异常值的判定值，表示计数函数，表示一个中的样本个数，表示在一个中无放回地抽取样本点
q
时判断样本点
q
为样本点
p
的临近点的个数，指样本点
p
与样本点
q
之间的欧几里得距离；
当得值大于预设得阈值时判定为异常值，剔除对应得第一样本点
。3.
根据权利要求2所述的深地油气精准导航随钻伽马测井曲线调优方法，其特征在于，所述步骤
S200
，具体包括：步骤
S210
，计算当前
batch
的剔除异常值的
GR
曲线数据的均值估计值和方差估计值：；；表示一个
batch
中的计数算子，
j
表示样本序号，
i
表示样本序号，
N
与的值相同；步骤
S220
，计算历史
GR
曲线数据的标准差，对当前
batch
的剔除异常值的
GR
曲线数据进行拉伸压缩变换：；表示拉伸压缩变换后得到的
GR
数据，表示当前
batch
的剔除异常值的
GR
曲线数据的第
k
个样本点的
GR
值；步骤
S230
，对所述拉伸压缩变换后得到的
GR
数据，基于历史
GR
曲线数据的均值标准值进行平移变换：；表示平移变换后得到的标准化后的随钻
GR
参数数据
。4.
根据权利要求2所述的深地油气精准导航随钻伽马测井曲线调优方法，其特征在于，所述步骤
S300
，具体包括：步骤
S310
，分别计算变换后的随钻
GR
参数数据的概率分布曲线和已钻井的
GR
参数数据的概率分布曲线；步骤
S320
，基于和计算概率分布曲线面积交集和概率分布曲线面积并集；步骤
S330
，基于所述和计算质量控制参数；步骤
S340
，比对所述质量控制参数
U
与设定的质量控制阈值的大小；当所述质量控制参数
U
大于或等于设定的质量控制阈值时，合理性判断结果为通过，进入步骤
S400
；当所述质量控制参数
U
小于设定的质量控制阈值时，合理性判断结果为不通过，进入...

【专利技术属性】
技术研发人员：田飞，张江云，底青云，郑文浩，杨永友，郑健，曹文静，
申请(专利权)人：中国科学院地质与地球物理研究所，
类型：发明
国别省市：