有机质生烃成因的地层超压预测方法技术

一种有机质生烃成因的地层超压预测方法，包括，优选参考井；对钻井实测样品总有机碳值标定并计算单井总有机碳预测曲线；优选与总有机质含量参数；建立相关常数与总有机碳含量函数关系；分析确定区域生烃门限；针对生烃门限以下地层，分析地震叠前道集质量及井震；优选目标区已钻井开展时深标定与子波提取；求取待预测井点或区域岩石纵横波速度、密度信息及常数与总有机碳含量函数关系式所需弹性参数；利用弹性参数及公式求取待预测井点或区域内总有机碳含量；利用单井预测总有机碳含量与实测地层孔隙压力数据得回归关系式；利用总有机碳和回归关系式计算待预测井点或区域内的超压预测成果；可较高精度完成烃源岩发育区域地层超压定量化预测。超压定量化预测。超压定量化预测。

【技术实现步骤摘要】
有机质生烃成因的地层超压预测方法


[0001]本专利技术属于油气田勘探
，尤其涉及一种有机质生烃成因的地层超压预测方法。

技术介绍

[0002]异常高压是含油气盆地普遍存在的现象。渗透性地层中存在的异常高压对钻井安全造成严重威胁。因此，含油气盆地钻前超压预测对钻井安全起到至关重要的作用。异常高压的形成机制多达数十种，而普遍存在且在地层超压形成过程中占据主导地位的是泥岩欠压实作用和有机质生烃作用形成的孔隙超压。
[0003]目前，针对泥岩欠压实成因地层孔隙超压定量化预测手段有Eaton法、平衡深度法等，并且这些方法在实践过程中取得了较好的应用效果。但是有机质生烃成因的超压机理与泥岩欠压实成因的超压机理完全不同，这些成熟的欠压实地层超压预测方法不再适用于有机质生烃成因的地层超压预测。目前基于Bowers(卸载曲线)方法虽然可以开展钻后生烃超压分析，但尚没有见到其在钻前应用的实例，其主要限制因素在于应用该方法时涉及参数众多，在钻前取值难度较大。而随着含油气盆地勘探程度的逐渐加深，越来越多的探井被部署在沉积盆地的深凹区域并钻探深部的烃源岩地层，这些都对有机质生烃成因的地层超压定量化预测技术提出迫切需求。

技术实现思路

[0004]本专利技术目的在于提供一种有机质生烃成因的地层超压预测方法,以解决有机质生烃成因的地层超压定量化预测的技术问题。
[0005]为实现上述目的，本专利技术的有机质生烃成因的地层超压预测方法的具体技术方案如下：
[0006]第一步，优选目标区内可用于地层生烃分析的多口参考井；参考井需测得有较多的总有机碳(TOC)及镜质体反射率(Ro)等实验信息，同时包含地层纵横波速度、电阻率、密度等测井信息及实测压力信息；
[0007]第二步，基于第一步取得的已钻井实测样品的总有机碳(TOC)值进行标定，利用Passey经验公式计算单井总有机碳(TOC)预测曲线；计算公式为：
[0008][0009][0010]式中：ΔlgR为声波时差曲线和电阻率曲线叠合幅度差；R
t
为电阻率测井曲线；R
基线
为基线电阻率；Δt为声波时差测井曲线；w(TOC)为有机碳含量；L
OM
为热变指数，反映有机质成熟度，与TOC质量分数有对应关系，也可用镜质体反射率R
O
代替；
[0011]第三步，将第一步取得的测井信息根据阵列声波及密度计算岩石物理弹性参数，
采用交会图方法优选与总有机质含量(TOC)相关的参数；
[0012]第四步，采用逐步回归或神经网络非线性映射方法建立相关常数与总有机碳含量的函数关系；该关系式为：TOC＝Fun(λ，K，V
p
，...)
ꢀꢀ
(1.3)；
[0013]第五步，基于目标区内各已钻井镜质体反射率(Ro)与深度交会分析确定区域生烃门限；
[0014]第六步，针对生烃门限以下地层，开展地震叠前道集质量及井震AVO分析，并进行道集目标处理进一步提高道集质量，为叠前弹性参数反演提供高质量的道集数据；
[0015]第七步，优选目标区已钻井开展时深标定与子波提取，以原始VSP(垂直地震剖面)数据为初始时深关系，优选地震子波频率、相位等关键参数，通过整体平移和局部拉伸保证合成记录与原始地震的波组关系具有较好一致性。已钻井的选择依据是钻至或钻穿可能发生超压的层段，井径质量较好，同时具备速度、密度且平面分布相对均匀的已钻井；
[0016]第八步，利用高精度初始模型，开展基于共反射点(CMP)道集的叠前弹性参数反演，求取待预测井点或区域的岩石纵横波速度、密度信息并同时计算公式(1.3)所需弹性参数；
[0017]第九步，利用第八步反演获得弹性参数以及第四步建立的公式(1.3)求取待预测井点或区域内总有机碳含量(TOC)；
[0018]第十步，利用第二步计算的单井预测总有机碳含量与实测地层孔隙压力数据进行交会获得回归关系式；
[0019]第十一步，利用第九步获得的总有机碳(TOC)和第十步获得的交会关系式计算待预测井点或区域内的超压预测成果。
[0020]本专利技术的有机质生烃成因的地层超压预测方法具有以下优点：该方法是通过综合运用地球化学分析、岩石物理分析、地震反演等手段，建立一套烃源岩发育区生烃成因机制的地层孔隙超压定量化预测的方法，该方法避免了现有有机质生烃成因地层超压钻前预测模型参数众多、确定困难的问题，为烃源岩发育区深部地层钻探面临的钻井液密度配比和井身结构设计问题提供科学指导，保证钻井安全高效施工。
附图说明
[0021]图1为本专利技术有机质生烃成因的地层超压预测方法流程示意图。
[0022]图2为本专利技术有机质生烃成因的地层超压预测方法单井总有机碳含量(TOC)预测效果图。
[0023]图3A为本专利技术有机质生烃成因的地层超压预测方法弹性参数中拉梅系数与总有机碳含量(TOC)交会图。
[0024]图3B为本专利技术有机质生烃成因的地层超压预测方法弹性参数中体积模量与总有机碳含量(TOC)交会图。
[0025]图3C为本专利技术有机质生烃成因的地层超压预测方法弹性参数中纵波速度与总有机碳含量(TOC)交会图。
[0026]图3D为本专利技术有机质生烃成因的地层超压预测方法弹性参数中密度与总有机碳含量(TOC)交会图。
[0027]图3E为本专利技术有机质生烃成因的地层超压预测方法弹性参数中泊松比与总有机
碳含量(TOC)交会图。
[0028]图3F为本专利技术有机质生烃成因的地层超压预测方法弹性参数中纵横波速比与总有机碳含量(TOC)交会图。
[0029]图4为本专利技术有机质生烃成因的地层超压预测方法镜质体反射率与深度交会确定区域生烃门限图。
[0030]图5为本专利技术有机质生烃成因的地层超压预测方法已钻井合成地震记录标定及岩性柱状图。
[0031]图6A为本专利技术有机质生烃成因的地层超压预测方法密度模型构建图。
[0032]图6B为本专利技术有机质生烃成因的地层超压预测方法的横波速度模型构建图。
[0033]图7为本专利技术有机质生烃成因的地层超压预测方法总有机碳与实测压力系数交会图。
[0034]图8为本专利技术有机质生烃成因的地层超压预测方法预测区域内生烃成因超压预测成果。
具体实施方式
[0035]为了更好地了解本专利技术的目的、结构及功能，下面结合附图，对本专利技术一种有机质生烃成因的地层超压预测方法做进一步详细描述。
[0036]实施例：
[0037]如图1所示，本专利技术有机质生烃成因的地层超压预测方法包括以下步骤：
[0038]第一步，优选目标区内可用于地层生烃分析的多口参考井；参考井需测得有较多的总有机碳(TOC)及镜质体反射率(Ro)等实验信息，同时包含地层纵横波速度、电阻率、密度等测井信息及实测压力信息；
[0039]第二步，基于第一步取得的已钻井实测样品的总有机碳(TOC本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种有机质生烃成因的地层超压预测方法，其特征在于，所述方法包括以下步骤：第一步，优选目标区内可用于地层生烃分析的多口参考井；参考井需拥有丰富的总有机碳(TOC)及镜质体反射率(Ro)等实验信息，同时包含地层纵横波速度、电阻率、密度等测井信息及实测压力信息；第二步，基于第一步取得的已钻井实测样品的总有机碳(TOC)值进行标定，利用Passey经验公式计算单井总有机碳(TOC)预测曲线；计算公式为：计算公式为：式中：ΔlgR为声波时差曲线和电阻率曲线叠合幅度差；R
t
为电阻率测井曲线；R
基线
为基线电阻率；Δt为声波时差测井曲线；w(TOC)为有机碳含量；L
OM
为热变指数，反映有机质成熟度，与TOC质量分数有对应关系，也可用镜质体反射率R
O
代替；第三步，将第一步取得的测井信息根据阵列声波及密度计算岩石物理弹性参数，采用交会图方法优选与总有机质含量(TOC)相关的参数；第四步，采用逐步回归或神经网络非线性映射方法建立相关常数与总有机碳含量的函数关系；该关系式为：TOC＝Fun(λ，K，V
p
，...)

【专利技术属性】
技术研发人员：郭军，周东红，张志军，何玉，周星，李英，秦童，李少轩，姜本厚，谭辉煌，陈平，汤国民，
申请(专利权)人：中海石油中国有限公司天津分公司，
类型：发明
国别省市：