[发明专利]一种煤层气储层识别与丰度预测的方法

摘要

本发明是非常规气藏气勘探综合测井数据对煤层气储层识别与丰度预测的方法，在岩芯上读取特殊岩性地层深度和厚度，在测井曲线上将相应深度的地层层段删除，对工区内的测井数据进行计算得到煤层的物性数据的平均值和方差，标准化处理后建立多参数煤层气丰度预测公式，得到地层含气量变化曲线，根据地层含气量变化曲线，进行煤层气储层识别和丰度预测。本发明方法简便易行。成本低，运行效率高，可准确识别薄的煤层，煤层气储层丰度预测结果可靠。

主权项

一种煤层气储层识别与丰度预测的方法，特点是采用以下步骤：1)钻井并在井中采集常规地球物理测井数据；2)消除储层埋深和压力对测井数据的影响；3)在岩芯剖面上，读取碳酸盐岩、火成岩特殊岩性地层深度位置和地层厚度，在测井曲线上将相应深度的地层层段数据删除，使其不参与统计；4)根据煤层的物性分布特点，对工区内的测井数据进行计算得到煤层的物性数据的平均值和方差，利用均值—方差法对测井资料进行标准化处理；5)在工区内已有实验室含气量测试结果参考井测井资料基础上，建立多参数煤层气丰度预测公式；6)对工区内没有进行实验室含气量测定的每口井中的测井曲线进行上述步骤2)至4)的处理，利用步骤5)建立的丰度预测公式对工区内没有进行实验室含气量测定的井进行含气量预测，得到每口井的地层含气量变化曲线；7)根据地层含气量变化曲线，进行煤层气储层识别，确定煤层气储层的厚度和埋藏深度，以及储层中煤层气丰度变化规律，步骤5)所述的建立多参数煤层气丰度预测公式过程是：(1)选取工区中同时具有测井数据和实验室含气量测试结果的井作为参考井，用于丰度预测公式的建立；对这些井中的测井数据进行步骤1)至4)的处理；(2)建立多参数煤层气丰度预测公式：Vgas＝A0+A1X1+A2X2+A3X3+A4X4+A5X5+A6X6   (1)式中：Vgas是煤层气含气量数据，A0是预测公式中的常数项；X1＝{x11，x21，…xn1}是地层的密度测井数据，其中，xn1是密度测井数据在目的层段中第n个点的数值；A1是密度变量的加权系数；X2＝{x12，x22，…xn2}是地层的纵波速度，其中，xn2是纵波速度数据在目的层段中第n个点的数值；A2是纵波速度变量的加权系数；X3＝{x13，x23，…xn3}是地层的孔隙度数据，其中，xn3是孔隙度测井数据在目的层段中第n个点的数值，A3是孔隙度变量的加权系数；X4＝{x14，x24，…xn4}是地层的伽马测井数据，其中，xn4是伽马测井数据在目的层段中第n个点的数值，A4是伽马变量的加权系数；X5＝{x15，x25，…xn5}是地层的电阻率测井数据，其中，xn5是电阻率测井数据在目的层段中第n个点的数值，A5是电阻率变量的加权系数；X6＝{x16，x26，…xn6}是井径变化数据，其中，xn6是井径数据在目的层段中第n个点的数值，A6是井径变量的加权系数；(3)利用参考井中已有的实验室测定的煤层气含量和测井数据作为已知数据，对公式(1)进行多元回归分析，得到下列方程组：(CCT)A0A1A2A3A4A5A6=CVgas1Vgas2Vgas3...Vgasn---(2)]]>其中：公式(2)中的A0、A1、A2、A3、A4、A5、A6为公式(1)中的加权系数，{Vgas1，Vgas2，…Vgasn}是参考井在目的层段内的实验室测定的煤层气含量数据值，其中的Vgasn是目的层段内第n个点的含气量数据值；公式(3)中的{x11，x21，…xn1}、{x12，x22，…xn2}、{x13，x23，…xn3}、{x14，x24，…xn4}、{x15，x25，…xn5}、{x16，x26，…xn6}的定义与公式(1)中的定义相同，是参考井中已知的不同类型的测井数据在目的层段内的离散采样值；(4)用乔里斯基(Cholesky)分解算法解公式(2)，就得出公式(1)的加权系数A0、A1、A2、A3、A4、A5、A6，完成丰度预测公式。