[发明专利]一种碎屑岩不整合结构体测井综合识别方法

摘要

本发明公开了一种碎屑岩不整合结构体测井综合识别方法，利用测井曲线划分碎屑岩不整合内部结构，通过多参数测井曲线拟合，构建不整合孔隙综合判别参数、岩性判别曲线、泥质含量曲线，结合不整合内部结构的特点，定量的划分不整合内部结构，为准确评价不整合的控藏作用，提高油气勘探与预测的有效性提供技术支撑。

主权项

1.一种碎屑岩不整合结构体测井综合识别方法，其特征在于，包括：(1)分析工区构造背景，确定不整合结构发育类型，完成目标钻井选取；(2)遴选不整合结构敏感测井曲线，排除曲线异常值，选取相应岩石骨架参数及流体参数，完成不整合结构敏感参数计算；具体包括：S1：计算不整合孔隙综合判别参数U，具体包括：S11、选取对不整合孔隙敏感的补偿中子CNL、密度DEN、电阻率R和自然伽马GR进行曲线变化幅度差计算，计算公式为ΔY_i＝(|Y_i‑1‑Y_i|+|Y_i+1‑Y_i|)/Y_i；其中，ΔY_i为ΔCNL、ΔDEN、ΔR或ΔGR，Y_i为当前深度点的测井曲线值，Y_i‑1和Y_i+1为当前深度点邻近的两点的测井曲线值；S12、对各个参数进行归一化计算其中，对于补偿中子CNL、电阻率R、自然伽马GR的变化幅度差采用如下公式进行归一化计算：Z_i＝(ΔY_i‑ΔY_min)/(ΔY_max‑ΔY_min)；其中，Z_i表示当前深度点上述变化幅度差曲线归一化数值，分别用Z_N、Z_R、Z_G表示，ΔY_i当前深度点该种测井曲线实际变化幅度差值，ΔY_max表示该种测井曲线的最大变化幅度差，ΔY_min表示该种测井曲线的最小变化幅度差；密度DEN变化幅度差采用如下公式进行归一化计算：Z_D＝1‑(ΔY′‑ΔY′_min)/(ΔY′_max‑ΔY′_min)；其中，Z_D表示当前深度点密度测井变化幅度差归一化数值，ΔY′表示当前深度点密度测井变化幅度差值，ΔY′_max表示密度测井的最大变化幅度差，ΔY′_min表示密度测井的最小变化幅度差；S13、构建不整合孔隙综合判别参数U利用归一化计算得到的变化幅度差Z_N、Z_D、Z_R及Z_G构建不整合孔隙综合判别参数U：U＝[P(Z_N)+P(Z_D)+P(Z_R)+P(Z_G)]/4；其中，P为相对各种曲线的函数，满足以下条件：式中，a、b为异常临界值，且有a<b；S2：构建岩性判别曲线密度测井计算总孔隙度的公式为：其中，ρ_b为密度测井曲线读出的测量值，ρ_f为孔隙中流体的密度值，ρ_ma为岩石骨架密度值；补偿中子测井计算总孔隙度的公式为：其中，为补偿中子曲线上读出的测量值，为岩石骨架中子值，为孔隙中流体的中子值；基于密度测井计算的总孔隙度和补偿中子测井计算的总孔隙度构建岩性判别曲线：其中，为补偿中子测井计算的总孔隙度，为密度测井计算的总孔隙度；S3：计算泥质含量曲线泥质含量的变化会引起GR曲线幅度的变化，可用于区分不整合三层结构,并可借此对比风化泥岩与正常泥岩，泥质含量指数I_GR计算公式为：其中，GR_log为实际测量值，GR_max为该曲线的最大测量值，GR_min为该测井曲线的最小测量值；则泥质含量曲线V_sh计算公式为：其中，GCUR为希尔奇指数，根据实验室取心分析资料确定，以上完成了不整合内部结构敏感参数计算；(3)不整合内部结构划分根据不整合内部结构的特点，通过不整合孔隙综合判别参数U、岩性判别曲线泥质含量曲线V_sh三者的结合，综合分析，从而完成不整合内部结构的定量划分。