CN111241460A - 一种复杂致密储层孔隙度计算方法 - Google Patents
一种复杂致密储层孔隙度计算方法 Download PDFInfo
- Publication number
- CN111241460A CN111241460A CN202010058727.0A CN202010058727A CN111241460A CN 111241460 A CN111241460 A CN 111241460A CN 202010058727 A CN202010058727 A CN 202010058727A CN 111241460 A CN111241460 A CN 111241460A
- Authority
- CN
- China
- Prior art keywords
- porosity
- logging
- por
- cnl
- correction
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Physics (AREA)
- Data Mining & Analysis (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Computational Mathematics (AREA)
- Pure & Applied Mathematics (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Algebra (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
一种复杂致密储层孔隙度计算方法,第一步,岩心归位,建立传统AC计算POR模型;第二步,计算△SP值;第三步,基于e(1‑△SP)对AC和CNL进行泥质校正;第四步,建立复杂致密储层孔隙度解释模型;第六步,对复杂致密储层进行孔隙度POR计算;经实例计算结果表明,本发明能准确计算复杂致密储层孔隙度,具有条件简单、适用范围广、解释精度高及易于推广的优点,为储层储量解释打下重要的基础,具有重要的应用前景和经济价值。
Description
技术领域
本发明涉及油气田勘探技术领域,具体涉及一种复杂致密储层孔隙度计算方法。
背景技术
储层孔隙度POR是进行储层储量计算的重要基础。关于国内外孔隙度计算有多种方法,主要有以下五种:一、AC与CNL、CNL与DEN等交会图解释孔隙度;二、采用wylli公式(孙波,2016),采用声波时差Δt、流体声波时差Δtf与骨架声波时差Δtma进行孔隙度Φ计算:Δt=Φ×Δtf+(1-Φ)×Δtma;三、采用孔隙度POR与声波时差AC、补偿密度DEN、补偿中子CNL、泥质含量Vsh、GR进行回归(郭晓磊,2018;司马立强等,2008;王自亮,2016;张鹏,2014,陈硕等,2019);四、基于核磁共振或ECS测井计算孔隙度或(侯克均等,2019;佘刚等,2019);五、基于自然电位计算泥质砂岩储层孔隙度(任杰等,2017;梁霄等,2017)。第一种方法:采用交会图计算孔隙度,这个方法存在问题是其粘土骨架参数不好确定,不同人有不同的选择,不同井的参数也会有差别(杨黎明等,2019),应用较困难,且人为性较大,其解释精度较低;第二种方法采用wylli公式计算储层孔隙度,但wylli公式对于孔隙度低(5-15%)和高孔隙度(>30%)不适用(王晓光等,2015),致密储层孔隙度一般小于15%,所以wylli公式计算孔隙度误差较大。第三种方法,主要建立各种回归分析,要计算准确孔隙度的前提条件需要提前计算好泥质含量Vsh值,然后用Vsh去校正各种计算模型的孔隙度,而恰恰是由于复杂储层含有高伽马砂岩,给储层泥质含量准确计算带来很大困难,而泥质含量解释准确与否,关系到孔隙度解释的精度,正因为如此,有些研究者直接采用AC等测井参数计算孔隙度,方法简便,但对于一套砂体内由于泥质含量有波动增加的情况,会造成其解释孔隙度偏高问题,计算结果与实际情况相差较大。总体来说,该种计算孔隙度方法对于复杂致密储层孔隙度计算条件较苛刻、适用性较低;第四种方法,主要采用特殊测井,比如核磁共振、成像测井或ECS测井,尽管特殊测井在复杂储层计算孔隙度有较高精度,但由于这些测井数据成本很高,所以数量很少,不宜推广;第五种方法,直接用自然电位计算计算储层孔隙度,该方法在缺少三孔隙度测井情况下有一定效果,但对于致密储层,由于SP异常幅度较低,孔隙度解释会偏低,此外SP与泥质含量、孔隙度、和地层水矿物度有关,所以单纯用SP计算高伽马储层孔隙度,其精度较低。
常规测井在高伽马储层泥质Vsh计算方面困难很大,尽管自然能谱曲线计算泥质含量Vsh比较准,但其因成本高,所测得井数量少,无法大面积使用。现有技术中针对含有高伽马复杂致密储层主要采用常规测井计算孔隙度,存在的共同缺陷即若没有在先计算泥质含量Vsh情况下,无法对其解释孔隙度进行泥质含量校正,从而导致其解释孔隙偏高问题。
发明内容
为了克服上述现有技术的缺陷,本发明的目的在于提供一种复杂致密储层孔隙度计算方法,采用常规储层的测井参数,不需要提前计算泥质含量Vsh,就能实现对声波时差AC和补偿中子CNL进行泥质含量动态校正,构建不受放射性(高伽马储层)影响的孔隙度计算模型,能较准确计算含有高伽马复杂致密储层孔隙度,为储层储量的解释模型打下重要基础,具有重要的应用前景和经济价值。
为了达到上述目的,本发明的技术方案为:
一种复杂致密储层孔隙度计算方法,包括以下步骤:
步骤一:岩心归位,建立传统AC计算POR模型
将岩心分析物性数据加载在测井wis数体中,根据综合声波时差和补偿密度DEN测井资料,对岩心分析数据进行岩心归位,建立岩心分析孔隙度与声波时差AC模型:
POR=a*AC+b (1)
式中:POR是岩心分析孔隙度,AC为测井声波时差,a和b为方程系数;
步骤二:计算△SP值:
式中:△SP为测井SP异常幅度差,幅度差越大,储层泥质含量越低,反之,储层泥质含量越高,SPmin为砂岩测井SP值,SPmax为泥岩处SP值,SP为某个深度的测井自然电位值;
步骤三:基于e(1-△SP)对AC和CNL进行泥质校正
△SP为测井SP异常幅度差,1-△SP越大,说明△SP越小,△SP越小有两种情况,一是泥质含量增加,二是砂岩本身泥质含量低,但物性致密,孔隙度物性值很小,其声波时差较小,AC<=215us/m;
若AC>215us/m,则利用(e(1-△SP))2动态校正因子对测井声波时差AC和补偿中子CNL进行泥质校正:
AC校正=AC-5*(e(1-△SP))2 (3)
CNL校正=CNL/(e(1-△SP))2 (4)
步骤四:建立复杂致密储层孔隙度解释模型
建立岩心分析孔隙度POR与AC、补偿密度DEN、测井光电截面PE、CNL、CNL校正、AC校正、GR相关性分析,优选POR与AC校正和CNL校正建立回归模型:
POR=b0*e(b1*T) (6)
式中:AC校正为泥质校正后测井声波时差AC,CNL校正为泥质校正测井补偿中子CNL,AC为测井声波时差,CNL为测井补偿中子,△SP为SP异常幅度差,b0、b1是方程系数,为常数;
步骤五:对复杂致密储层进行孔隙度POR计算
为了避免第三步骤中第二种情况AC<=215us/m时被过度泥质校正,第二种情况AC<=215us/m时采用公式(1)计算孔隙度POR;
AC>215us/m时,利用步骤四中公式(5)和公式(6)进行致密储层孔隙度POR计算。
本发明的有益效果是:
本发明是一种复杂致密储层孔隙度计算方法,不需要提前计算泥质含量Vsh,就能实现对声波时差AC和补偿中子CNL进行泥质含量动态校正,有效避免现有技术中存在若不先计算泥质含量Vsh,就无法实现对其解释孔隙度进行泥质校正的技术缺陷。本发明建立的基于e(1-△SP)动态因子泥质校正的孔隙度解释模型,不仅适用于常规储层,也可适用于含有高伽马砂岩的储层。因此,具有条件简单、适用范围广、方法简便、解释精度高及易于推广的特点,为致密复杂储层储量参数的进一步解释打下重要的基础,具有重要的应用前景和经济价值。
附图说明
图1为复杂致密储层孔隙度计算的流程框图。
图2为实施例中的孔隙度计算模型。
图3为实施例中的传统AC计算孔隙度与岩心分析孔隙度交会图。
图4为实施例中的AC校正/CNL校正计算孔隙度与岩心分析孔隙度交会图。
具体实施方式
下面结合实例对本发明的具体应用作详细说明。
参照图1,一种复杂致密储层孔隙度计算方法,包括以下步骤:
步骤一:岩心归位,建立传统AC计算POR模型
将岩心分析物性数据加载在测井wis数体中,根据综合声波时差和补偿密度DEN等测井资料,对岩心分析数据进行岩心归位。建立岩心分析孔隙度与声波时差AC模型(表1),其相关性R=0.772:
POR=0.296*AC-56.955 (1)
式中:POR是岩心分析孔隙度,AC为测井声波时差,a=0.296、b=-56.955。
表1 AC计算POR模型参数
模型 | R | R平方 | 调整后R平方 | 标准偏斜度错误 |
1 | .772<sup>a</sup> | .597 | .595 | 1.955326 |
a.预测值:(常数),AC
步骤二:计算△SP值
式中:△SP为测井SP异常幅度差,幅度差越大,储层泥质含量越低,反之,储层泥质含量越高,SPmin为砂岩测井SP值,SPmax为泥岩处SP值,SP为某个深度的测井自然电位值。
步骤三:基于e(1-△SP)对AC和CNL进行泥质校正
△SP为测井SP异常幅度差,1-△SP越大,说明△SP越小,△SP越小有两种情况,一是泥质含量增加,二是砂岩本身泥质含量低,但物性致密,孔隙度等物性值很小,其声波时差AC<=215us/m。若AC>215us/m,则利用(e(1-△SP))2动态校正因子对测井声波时差AC和补偿中子CNL进行泥质校正:
AC校正=AC-5*(e(1-△SP))2 (3)
CNL校正=CNL/(e(1-△SP))2 (4)
步骤四:建立复杂致密储层孔隙度解释模型
建立岩心分析孔隙度POR与AC、DEN(补偿密度)、PE(测井光电截面)、CNL、CNL校正、AC校正、GR相关性分析(表2),优选POR与AC校正和CNL校正建立回归模型。
表2岩心分析孔隙度POR与其它参数相关性分析
从表2中,没有进行泥质含量校正,AC与POR相关性0.772,CNL与POR相关性0.534,AC经过5*(e1-△SP)2泥质含量校正后,其相关性升至0.877,相关性增加0.105,CNL经过(e1 -△SP)2泥质含量校正后,其相关性达0.813,相关性提高0.279。上述两个参数相关性增加一个数量级,泥质校正效果很明显。因此,优选POR与AC校正、CNL校正建立回归模型(表3):
表3本发明模型总计及参数评估
因变量:POR
自变数为AC校正/CNL校正。
POR=16.542*e(-0.025*T) (6)
式中:POR为孔隙度,AC校正为泥质校正后测井声波时差AC,CNL校正为泥质校正测井补偿中子CNL,AC为测井声波时差,CNL为测井补偿中子,△SP为SP异常幅度差,b0=16.542、b1=-0.025;
表4模型相关性表
模型 | R | R平方 |
1 | .894 | .800 |
模型的相关性高达0.894,经泥质校正后的孔隙度POR解释模型可靠(表4)。
步骤五:对复杂致密储层进行孔隙度POR计算
为了避免第三步骤中第二种情况(AC<=215us/m)被过度泥质校正,第二种情况(AC<=215us/m)采用公式(1)计算孔隙度POR。AC>215us/m时,利用步骤四中公式(5)和公式(6)进行致密储层孔隙度POR计算。
以A267井为例,通过与岩心分析分析POR对比分(附图3、附图4、表5)析,传统方法计算复杂致密储层孔隙度平均绝对误差1.46%,相对误差21.29%,本发明平均绝对误差0.55%,相对误差6.87%,绝对误差相比统统方法减少62.3%,相对误差减少14.42%。传统AC计算孔隙度在常规储层(GR<100)地层精度相对较好,但是其GR>=100,在含有放射性储层,其相对误差较大,且误差范围较宽,相对误差10-66.73%,平均相对误差22.86%,而本发明在GR>=100,其相对误差3-18.63%,绝大部分相对误差小于10%,平均相对误差8.14%。在高GR储层(GR>=100),本发明计算孔隙度比传统计算的相对误差减少14.72%。
计算孔隙度POR与岩心分析POR对比来看,传统AC计算孔隙度与岩心分析POR孔隙度模型相关性0.77(图3),而AC校正/CNL校正计算孔隙度相关性高达0.89(图4),本发明相关性比传统提高0.12,解释精度提高一个数量级别。综上所述,本发明的复杂致密储层孔隙度不依赖必须先计算泥质含量Vsh条件下,就能实现对声波时差AC和补偿中子CNL进行泥质含量动态校正,其孔隙度解释精度显著比传统方法高,该计算方法不仅适用于常规储层,也适用含有高伽马储层,计算效果良好。
表5孔隙度POR计算误差分析(A267井)
Claims (2)
1.一种复杂致密储层孔隙度计算方法,其特征在于,包括以下步骤:
步骤一:岩心归位,建立传统AC计算POR模型;
步骤二:计算△SP值;
步骤三:基于e(1-△SP)对AC和CNL进行泥质校正;
步骤四:建立复杂致密储层孔隙度解释模型;
步骤五:对复杂致密储层进行孔隙度POR计算。
2.根据权利要求1所述的一种复杂致密储层孔隙度计算方法,其特征在于,具体为:
步骤一:岩心归位,建立传统AC计算POR模型;
将岩心分析物性数据加载在测井wis数体中,根据综合声波时差和补偿密度DEN测井资料,对岩心分析数据进行岩心归位,建立岩心分析孔隙度与声波时差AC模型:
POR=a*AC+b (1)
式中:POR是岩心分析孔隙度,AC为测井声波时差,a和b为方程系数;
步骤二:计算△SP值
式中:△SP为测井SP异常幅度差,幅度差越大,储层泥质含量越低,反之,储层泥质含量越高,SPmin为砂岩测井SP值,SPmax为泥岩处SP值,SP为某个深度的测井自然电位值;
步骤三:基于e(1-△SP)对AC和CNL进行泥质校正
△SP为测井SP异常幅度差,1-△SP越大,说明△SP越小,△SP越小有两种情况,一是泥质含量增加,二是砂岩本身泥质含量低,但物性致密,孔隙度物性值很小,其声波时差较小,AC<=215us/m;
若AC>215us/m时,则利用(e(1-△SP))2动态校正因子对测井声波时差AC和补偿中子CNL进行泥质校正:
AC校正=AC-5*(e(1-△SP))2 (3)
CNL校正=CNL/(e(1-△SP))2 (4)
步骤四:建立复杂致密储层孔隙度解释模型
建立岩心分析孔隙度POR与AC、补偿密度DEN、测井光电截面PE、CNL、CNL校正、AC校正、GR相关性分析,优选POR与AC校正和CNL校正建立回归模型:
POR=b0*e(b1*T) (6)
式中:AC校正为泥质校正后测井声波时差AC,CNL校正为泥质校正测井补偿中子CNL,AC为测井声波时差,CNL为测井补偿中子,△SP为SP异常幅度差,b0、b1是方程系数,为常数;
步骤五:对复杂致密储层进行孔隙度POR计算
为了避免第三步骤中第二种情况AC<=215us/m时被过度泥质校正,第二种情况AC<=215us/m时采用公式(1)计算孔隙度POR;
AC>215us/m时,利用步骤四中公式(5)和公式(6)进行致密储层孔隙度POR计算。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010058727.0A CN111241460B (zh) | 2020-01-19 | 2020-01-19 | 一种复杂致密储层孔隙度计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010058727.0A CN111241460B (zh) | 2020-01-19 | 2020-01-19 | 一种复杂致密储层孔隙度计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111241460A true CN111241460A (zh) | 2020-06-05 |
CN111241460B CN111241460B (zh) | 2021-08-24 |
Family
ID=70864924
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010058727.0A Active CN111241460B (zh) | 2020-01-19 | 2020-01-19 | 一种复杂致密储层孔隙度计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111241460B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111812736A (zh) * | 2020-07-22 | 2020-10-23 | 西南石油大学 | 一种致密砂岩无水气藏含气性评价方法 |
CN112228037A (zh) * | 2020-09-04 | 2021-01-15 | 长江大学 | 一种基于声波传播理论的孔隙度解释方法和装置 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1381734A (zh) * | 2002-03-28 | 2002-11-27 | 王振信 | 碳氢比地层流体饱和度测井方法 |
US20050281792A1 (en) * | 1997-08-13 | 2005-12-22 | Short Jay M | Phytases, nucleic acids encoding them and methods of making and using them |
CN101832133A (zh) * | 2010-01-28 | 2010-09-15 | 中国石油集团川庆钻探工程有限公司 | 密度孔隙度和中子孔隙度差值储层流体类型判别方法 |
CN102052070A (zh) * | 2010-12-13 | 2011-05-11 | 中国石油集团川庆钻探工程有限公司 | 复杂岩性碎屑岩地层组分测井处理方法 |
CN103821505A (zh) * | 2014-02-28 | 2014-05-28 | 中国石油大学(北京) | 砂岩油气输导层地球物理-地质-地球化学检测方法及装置 |
CN104636819A (zh) * | 2014-12-31 | 2015-05-20 | 中国石油天然气集团公司 | 一种利用储层有效厚度加权系数进行储层定量产能预测的方法 |
US20150260034A1 (en) * | 2014-03-14 | 2015-09-17 | Schlumberger Technology Corporation | Determining mineralogy of an earth formation using linear regressions |
CN106468172A (zh) * | 2016-09-30 | 2017-03-01 | 西安石油大学 | 一种超低渗砂岩油藏低阻储层测井解释方法 |
-
2020
- 2020-01-19 CN CN202010058727.0A patent/CN111241460B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20050281792A1 (en) * | 1997-08-13 | 2005-12-22 | Short Jay M | Phytases, nucleic acids encoding them and methods of making and using them |
CN1381734A (zh) * | 2002-03-28 | 2002-11-27 | 王振信 | 碳氢比地层流体饱和度测井方法 |
CN101832133A (zh) * | 2010-01-28 | 2010-09-15 | 中国石油集团川庆钻探工程有限公司 | 密度孔隙度和中子孔隙度差值储层流体类型判别方法 |
CN102052070A (zh) * | 2010-12-13 | 2011-05-11 | 中国石油集团川庆钻探工程有限公司 | 复杂岩性碎屑岩地层组分测井处理方法 |
CN103821505A (zh) * | 2014-02-28 | 2014-05-28 | 中国石油大学(北京) | 砂岩油气输导层地球物理-地质-地球化学检测方法及装置 |
US20150260034A1 (en) * | 2014-03-14 | 2015-09-17 | Schlumberger Technology Corporation | Determining mineralogy of an earth formation using linear regressions |
CN104636819A (zh) * | 2014-12-31 | 2015-05-20 | 中国石油天然气集团公司 | 一种利用储层有效厚度加权系数进行储层定量产能预测的方法 |
CN106468172A (zh) * | 2016-09-30 | 2017-03-01 | 西安石油大学 | 一种超低渗砂岩油藏低阻储层测井解释方法 |
Non-Patent Citations (2)
Title |
---|
李博克 等: "白豹油田长6 油藏储层岩矿及孔隙结构特征", 《西安石油大学学报( 自然科学版)》 * |
石先达 等: "Z油田阜宁组低—特低孔储层孔隙度测井计算改进方法", 《复杂油气藏》 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111812736A (zh) * | 2020-07-22 | 2020-10-23 | 西南石油大学 | 一种致密砂岩无水气藏含气性评价方法 |
CN111812736B (zh) * | 2020-07-22 | 2021-03-16 | 西南石油大学 | 一种致密砂岩无水气藏含气性评价方法 |
CN112228037A (zh) * | 2020-09-04 | 2021-01-15 | 长江大学 | 一种基于声波传播理论的孔隙度解释方法和装置 |
Also Published As
Publication number | Publication date |
---|---|
CN111241460B (zh) | 2021-08-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Li et al. | Uncertainty quantification of reservoir performance predictions using a stochastic optimization algorithm | |
CN108376188B (zh) | 黄土湿陷性的评价计算方法 | |
CN107783187B (zh) | 一种将测井速度和地震速度结合建立三维速度场的方法 | |
CN111255435B (zh) | 一种复杂储层泥质含量计算方法 | |
CN111241460B (zh) | 一种复杂致密储层孔隙度计算方法 | |
CN105893679B (zh) | 低产水平井续流修正试井解释方法 | |
CN111502652B (zh) | 一种三孔介质气藏水平井产量递减及生产动态预测方法 | |
CN110955982B (zh) | 一种变质岩储层渗透率计算方法、装置及计算机存储介质 | |
CN110727035A (zh) | 一种低渗强非均质气藏气水层识别方法 | |
CN111706318B (zh) | 一种确定低渗储层剩余油分布状况的方法 | |
CN107130955B (zh) | 井底流压的确定方法及储集体天然能量的确定方法 | |
CN111476448B (zh) | 一种多层合采直井动态产量劈分方法 | |
CN109667576B (zh) | 一种高矿化度成因低阻油层测井识别方法 | |
CN107526117B (zh) | 基于自动编码和超限学习联合网络的声波速度预测方法 | |
CN109655394B (zh) | 一种孔喉结构参数约束下的核磁共振t2谱渗透率计算方法 | |
CN112818501B (zh) | 基于动态监测数据校正碳酸盐岩油藏静态渗透率的方法 | |
CN113464129B (zh) | 煤层气藏平均地层压力计算方法和装置 | |
CN106593423A (zh) | 一种油藏流体类型的识别方法及装置 | |
CN113123783B (zh) | 一种精细评价地层体积含量的方法 | |
CN105239998A (zh) | 一种基于海上油田储层细分的部井方法 | |
CN113622908B (zh) | 一种水侵气藏废弃地层压力的确定方法 | |
CN116930023A (zh) | 致密砂岩相控分类孔隙度测井精细解释方法和装置 | |
CN112561223B (zh) | 油气井增产潜力的定量评价方法 | |
CN106526694A (zh) | 一种致密油藏识别方法及装置 | |
CN114109349B (zh) | 一种致密砂岩储层孔隙度指数/饱和度指数的确定方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |