矿井储层有效厚度下限的获取方法技术

本发明专利技术公开了一种矿井储层有效厚度下限的获取方法。本发明专利技术的基本原理为在待测地层内取岩心，对该岩心分别进行实验室化验分析得到该岩心的岩心孔隙度数据和岩心渗透率数据，然后向所述岩心注入汞，得到注入汞饱和度达到50％时对应的毛管压力，即汞饱和度中值压力，将汞饱和度中值压力与岩心孔隙度、岩心渗透率分别拟合，得到二者的回归关系式。通过一阶导数、二阶导数等方式求取该幂函数曲线的曲率最大值，此处对应的孔隙度、渗透率即为下限值。本发明专利技术通过采用能更好反应岩石孔喉结构的汞饱和度中值压力参数计算储层物性下限值，弥补了常规方法只能通过经验统计方法得到储层物性下限的不足，提高了获取的储层物性下限的准确性，同时提高了储层物性下限的效率。

【技术实现步骤摘要】
【专利摘要】本专利技术公开了一种。本专利技术的基本原理为在待测地层内取岩心，对该岩心分别进行实验室化验分析得到该岩心的岩心孔隙度数据和岩心渗透率数据，然后向所述岩心注入汞，得到注入汞饱和度达到50％时对应的毛管压力，即汞饱和度中值压力，将汞饱和度中值压力与岩心孔隙度、岩心渗透率分别拟合，得到二者的回归关系式。通过一阶导数、二阶导数等方式求取该幂函数曲线的曲率最大值，此处对应的孔隙度、渗透率即为下限值。本专利技术通过采用能更好反应岩石孔喉结构的汞饱和度中值压力参数计算储层物性下限值，弥补了常规方法只能通过经验统计方法得到储层物性下限的不足，提高了获取的储层物性下限的准确性，同时提高了储层物性下限的效率。【专利说明】
本专利技术涉及地球物理测井
，具体地指一种矿井储层有效厚度下限的获取 方法。
技术介绍
在油气田勘探开发过程中，需要测井技术人员对获得的井下测井资料进行及时的 测井解释并划分出有效储层（油层、气层），工作的关键在于确定划分有效储层的标准（也 即确定孔隙度和渗透率下限值）。现有的为含油产状法， 该方法详细记载与中国海洋石油南海西部公司企业标准Q/HXJ3026-93中，然而在长时间 的生产实践中技术人员发现上述方法获取的矿井储层有效厚度下限的准确性较低，究其原 因主要是该方法的主观性较大，用到了很多估算的过程。另外，该方法中一个重要的测量参 数为含油性，对于气田这种方法并不适用。
技术实现思路
本专利技术的目的就是要提供一种，利用该方法可 显著提高获取的矿井储层有效厚度下限的准确性，并且该方法对气田同样适用，具有较强 的通用性。 为实现此目的，本专利技术所设计的，其特征在于， 它包括如下步骤： 步骤1 :在待测地层内取岩心，对该岩心分别进行实验室化验分析得到该岩心的 岩心孔隙度数据和岩心渗透率数据，然后向所述岩心注入汞，得到注入汞饱和度达到50% 时对应的毛管压力，即汞饱和度中值压力； 步骤2:将岩心孔隙度数据与汞饱和度中值压力建立第一幂函数关系式yi = alXlb1，其中，Xl为岩心孔隙度数据，yi为汞饱和度中值压力， ai和h均为幂函数常数；将岩 心渗透率数据与汞饱和度中值压力建立第二幂函数关系式y2 = a2x2b2,其中，x2为岩心渗透 率数据，y 2为汞饱和度中值压力，a2和b2均为幂函数常数； 步骤3 :通过非线性拟合的方法对第一幂函数关系式和第二幂函数关系式进行幂 函数拟合求取幂函数常数bp a2和b2 ; 步骤4 :在第一幂函数关系式中，曲率最大处表示汞饱和度中值压力变化率最大， 此处即为幂函数曲线突变点，该幂函数曲线突变点对应的X轴数值为岩心孔隙度下限值； 在第二幂函数关系式中，曲率最大处表示汞饱和度中值压力变化率最大，此处即 为幂函数曲线突变点，该幂函数曲线突变点对应的X轴数值为岩心渗透率下限值； 步骤5 :将第一幂函数关系式yi = alXlbl代入第一幂函数关系式的曲率方程； 【权利要求】1. 一种，其特征在于，它包括如下步骤： 步骤1 :在待测地层内取岩心，对该岩心分别进行实验室化验分析得到该岩心的岩心 孔隙度数据和岩心渗透率数据，然后向所述岩心注入汞，得到注入汞饱和度达到50%时对 应的毛管压力，即汞饱和度中值压力； 步骤2 :将岩心孔隙度数据与汞饱和度中值压力建立第一幂函数关系式yi = alXlb1，其 中，Xl为岩心孔隙度数据，yi为汞饱和度中值压力，％和h均为幂函数常数；将岩心渗透率 数据与汞饱和度中值压力建立第二幂函数关系式y 2 = a2x2b2,其中，x2为岩心渗透率数据， y2为汞饱和度中值压力，a2和b2均为幂函数常数； 步骤3 :通过非线性拟合的方法对第一幂函数关系式和第二幂函数关系式进行幂函数 拟合求取幂函数常数ap bp a2和b2 ; 步骤4 :在第一幂函数关系式中，曲率最大处表示汞饱和度中值压力变化率最大，此处 即为幂函数曲线突变点，该幂函数曲线突变点对应的X轴数值为岩心孔隙度下限值； 在第二幂函数关系式中，曲率最大处表示汞饱和度中值压力变化率最大，此处即为幂 函数曲线突变点，该幂函数曲线突变点对应的X轴数值为岩心渗透率下限值； 步骤5 :将第一幂函数关系式yi = alXlbl代入第一幂函数关系式的曲率方程；其中，y/为第一幂函数关系式的一阶导，yi"为第一幂函数关系式的二阶导； 得到第一幂函数关系式的曲率将第二幂函数关系式y2 = a2x2b2代入第二幂函数关系式的曲率方程；其中，y，为第二幂函数关系式的一阶导，y2"为第二幂函数关系式的二阶导； 得到第二幂函数关系式的曲率：步骤6 :计算第一幂函数关系式的曲率的最大值，即对第一幂函数关系式的曲率根据 如下公式求一阶导数：在第一幂函数关系式的曲率的一阶导数为零时，得到的X1的值为第一幂函数关系式的 曲率的最大值； 计算第二幂函数关系式的曲率的最大值，即对第二幂函数关系式的曲率根据如下公式 求一阶导数：在第二幂函数关系式的曲率的一阶导数为零时，得到的^的值为第二幂函数关系式的 曲率的最大值； 步骤S7 :将上述第一幂函数关系式的曲率一阶导数的关系式简化为：将上述第二幂函数关系式的曲率一阶导数的关系式简化为：步骤S8 :对于第一幂函数关系式的曲率一阶导数的关系式，由于第一幂函数关系式中 ai和h均不为零，若和h有一个为零，则第一幂函数关系式变为yi = 0,此时第一幂函 数关系式的曲率没有最大点；且h不为1，若h为1则第一幂函数关系式变为yi = alXl，没 有曲率最大点；同时，Xl大于零，由于Xl代表岩心孔隙度数据，所以 Xl必然大于零，因此第 一幂函数关系式的曲率的一阶导数为零时的条件为：则有：将％和h代入上式并求出Xl即可得到岩心孔隙度数据下限； 对于第二幂函数关系式的曲率一阶导数的关系式，由于第二幂函数关系式中％和132均 不为零，若a2和b2有一个为零，则第二幂函数关系式变为y2 = 0,此时第二幂函数关系式的 曲率没有最大点；且b2不为1，若b2为1则第二幂函数关系式变为y 2 = a2x2，没有曲率最大 点；同时，x2大于零，由于x 2代表岩心渗透率数据，所以x2必然大于零，因此第二幂函数关 系式的曲率的一阶导数为零时的条件为： 则有：将a2和b2代入上式并求出x2即可得到岩心渗透率数据下限。【文档编号】G01N15/08GK104122182SQ201410307584【公开日】2014年10月29日 申请日期:2014年6月30日 优先权日:2014年6月30日 【专利技术者】谢玉洪, 周家雄, 何胜林, 胡向阳, 张恒荣 申请人:中国海洋石油总公司, 中海石油(中国)有限公司湛江分公司, 武汉时代地智科技股份有限公司本文档来自技高网...

【技术保护点】
一种矿井储层有效厚度下限的获取方法，其特征在于，它包括如下步骤：步骤1：在待测地层内取岩心，对该岩心分别进行实验室化验分析得到该岩心的岩心孔隙度数据和岩心渗透率数据，然后向所述岩心注入汞，得到注入汞饱和度达到50％时对应的毛管压力，即汞饱和度中值压力；步骤2：将岩心孔隙度数据与汞饱和度中值压力建立第一幂函数关系式y1＝a1x1b1，其中，x1为岩心孔隙度数据，y1为汞饱和度中值压力，a1和b1均为幂函数常数；将岩心渗透率数据与汞饱和度中值压力建立第二幂函数关系式y2＝a2x2b2，其中，x2为岩心渗透率数据，y2为汞饱和度中值压力，a2和b2均为幂函数常数；步骤3：通过非线性拟合的方法对第一幂函数关系式和第二幂函数关系式进行幂函数拟合求取幂函数常数a1、b1、a2和b2；步骤4：在第一幂函数关系式中，曲率最大处表示汞饱和度中值压力变化率最大，此处即为幂函数曲线突变点，该幂函数曲线突变点对应的X轴数值为岩心孔隙度下限值；在第二幂函数关系式中，曲率最大处表示汞饱和度中值压力变化率最大，此处即为幂函数曲线突变点，该幂函数曲线突变点对应的X轴数值为岩心渗透率下限值；步骤5：将第一幂函数关系式y1＝a1x1b1代入第一幂函数关系式的曲率方程；k1=|y1′′|(1+y1′2)32]]>其中，y1′为第一幂函数关系式的一阶导，y1″为第一幂函数关系式的二阶导；得到第一幂函数关系式的曲率k1=|a1b1(b1-1)x1(b1-2)|(1+a12b12x1(2b1-2))32]]>将第二幂函数关系式y2＝a2x2b2代入第二幂函数关系式的曲率方程；k2=|y2′′|(1+y2′2)32]]>其中，y2′为第二幂函数关系式的一阶导，y2″为第二幂函数关系式的二阶导；得到第二幂函数关系式的曲率：k2=|a2b2(b2-1)x2(b2-2)|(1+a22b22x2(2b2-2))32]]>步骤6：计算第一幂函数关系式的曲率的最大值，即对第一幂函数关系式的曲率根据如下公式求一阶导数：k1′=a1b1(b1-1)(b1-2)x1(b1-3)(1+a12b12x1(2b1-2))32-a1b1(b1-1)x1(b1-2)32(1+a12b12x1(2b1-2))12a12b122(b1-1)x1(2b1-3)(1+a12b12x1(2b1-2))3]]>在第一幂函数关系式的曲率的一阶导数为零时，得到的x1的值为第一幂函数关系式的曲率的最大值；计算第二幂函数关系式的曲率的最大值，即对第二幂函数关系式的曲率根据如下公式求一阶导数：k2′=a2b2(b2-1)(b2-2)x2(b2-3)(1+a22b22x2(2b2-2))32-a2b2(b2-1)x2(b2-2)32(1+a22b22x2(2b2-2))12a22b222(b2-1)x2(2b2-3)(1+a22b22x2(2b2-2))3]]>在第二幂函数关系式的曲率的一阶导数为零时，得到的x2的值为第二幂函数关系式的曲率的最大值；步骤S7：将上述第一幂函数关系式的曲率一阶导数的关系式简化为：k1′=a1b1(b1-1)x1(b1-3)[(b1-2)(1+a12b12x1(2b1-2))32-3a12b12(b1-1)x1(2b1-2)(1+a12b12x1(2b1-2))12](1+a12b12x1(2b1-2))3]]>将上述第二幂函数关系式的曲率一阶导数的关系式简化为：k2′=a2b2(b2-1)x2(b2-3)[(b2-2)(1+a22b22x2(2b2-2))32-3a22b22(b2-1)x2(2b2-2)(1+a22b22x2(2b2-2))12](1+a22b22x2(2b2-2))3]]>步骤S8：对于第一幂函数关系式的曲率一阶导数的关系式，由于第一幂函数关系式中a1和b1均不为零，若a1和b1有一个为零，则第一幂函数关系式变为y1＝0，此时第一幂函数关系式的曲率没有最大点；且b1不为1，若b1为1则第一幂函数关系式变为y1＝a1x1，没有曲率最大点；同时，x1大于零，由于x1代表岩心孔隙度数据，所以x1必然大于零，因此第一幂函数关系式的曲率的一阶导数为零时的条件为：(b1-2)+(b1-2)a12b12x1(2b1-2)-3(b1-1)a12b12x1(2b1-2)=0]]>则有：x1(2b1-2)=b1-2(2b1-1)a12b12]]>将a1和b1代入上式并求出x1即可得到岩心孔隙度数据下限；对于第二幂函数关系式的曲率一阶导数的关系式，由于第二幂函数关系式中a2和b2均不为零，若a2和b2有一个为零，则第二幂函数关系式变为y2＝0，此时第二幂函数关系式的曲率没有最...

【技术特征摘要】

【专利技术属性】
技术研发人员：谢玉洪，周家雄，何胜林，胡向阳，张恒荣，
申请(专利权)人：中国海洋石油总公司，中海石油中国有限公司湛江分公司，武汉时代地智科技股份有限公司，
类型：发明
国别省市：北京;11