CN104879126B - 一种基于阵列感应测井的碎屑岩储层流体识别方法 - Google Patents
一种基于阵列感应测井的碎屑岩储层流体识别方法 Download PDFInfo
- Publication number
- CN104879126B CN104879126B CN201510246298.9A CN201510246298A CN104879126B CN 104879126 B CN104879126 B CN 104879126B CN 201510246298 A CN201510246298 A CN 201510246298A CN 104879126 B CN104879126 B CN 104879126B
- Authority
- CN
- China
- Prior art keywords
- mrow
- array induction
- induction logging
- msub
- mover
- 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.)
- Active
Links
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明公开了一种基于阵列感应测井的碎屑岩储层流体识别方法，首先利用区块已试油的储层，读取层位阵列感应测井电阻率值，然后计算阵列感应斜率特征参数与相关特征参数，依据试油结论分系列建立阵列感应斜率特征参数与相关特征参数交会图，确定油、气、水分界区域；最后对一口新井目的层位的流体进行准确判识。本发明实现了在不考虑电阻率绝对值大小情况下，利用基于阵列感应测井资料确定储层流体性质的方法，能够准确、可靠的对储层流体进行判识。
Description
技术领域
本发明属于石油勘探中的复杂油气藏储层流体性质测井评价技术领域，涉及一种针对低阻油气层快速判别流体性质的方法，具体涉及一种基于阵列感应测井的碎屑岩储层流体识别方法。
背景技术
储层所含流体性质是测井解释评价的主要内容之一，通常储层流体的识别主要依赖于储层的电学性质，当储层岩性较纯、物性较好、孔隙单一、孔喉匹配关系好的情况下，利用储层的电阻率值可以对储层所含流体做出准确判定，即：
式中:Rt为储层的电阻率，R1为区域油层电阻率下限值，R2为区域水层电阻率上限值。
当储层类型复杂时,岩性和流体性质同时对储层电性产生影响,储层的电阻率绝对值无法正确反映储层流体性质；特别是当区域存在低电阻率油气层时，依靠电阻率值的大小去判断储层流体性质往往得到相反的结论；这种情况下计算的储层含油饱和度严重失真，油气层漏失现象普遍存在，致使测井解释的可信度严重下降。
发明内容
本发明的目的在于提供一种基于阵列感应测井的碎屑岩储层流体识别方法，能够在不考虑储层电阻率绝对值大小的情况下，快速、准确、连续的判断储层流体的性质。
为达到上述目的，本发明采用的技术方案为：
一种基于阵列感应测井的碎屑岩储层流体识别方法，包括以下步骤：
1)选择区域内经过试油验证的储层，获取不同试油层位的阵列感应测井电阻率值；
2)根据步骤1)中得到的阵列感应测井电阻率值，计算得到阵列感应斜率特征参数及相关特征参数；
3)根据步骤2)得到的阵列感应斜率特征参数及相关特征参数，建立阵列感应斜率特征参数与相关特征参数的交会图，并在交会图上确定出油、气、水的界限区域；
4)对于待识别储层流体性质的新井，先计算其阵列感应斜率特征参数及相关特征参数，并将计算结果投射到交会图上，再按照步骤3)得到的油、气、水的界限区域识别该井中储层流体的性质。
在进行步骤1)时，收集区域内具有阵列感应测井资料并且经过试油验证的储层，获取具有代表性的不同层位的阵列感应测井电阻率值。
所述的具有代表性的不同层位是指包含有油层、气层、水层的若干层位。
所述阵列感应斜率特征参数通过以下公式计算得到：
其中，k为阵列感应斜率特征参数；
n为阵列感应测井径向上的不同探测深度曲线数；
yi为不同探测深度的电阻率；
xi为不同探测深度电阻率对应的横坐标；
为xi的平均值，为yi的平均值。
所述相关特征参数通过以下公式计算得到：
其中，r为相关特征参数；
n为阵列感应测井径向上的不同探测深度曲线数；
yi为不同探测深度的电阻率；
xi为不同探测深度电阻率对应的横坐标；
为xi的平均值，为yi的平均值。
所述n＝5，yi和xi中的i为1～5的正整数，y1～y5分别为10英寸、20英寸、30英寸、60英寸和90英寸探测深度的电阻率值，x1～x5分别为1、2、3、4、5。
所述步骤3)在建立交会图时，以不同的试油结论为系列，以步骤2)计算的阵列感应斜率特征参数为纵坐标，以相关特征参数为横坐标，绘制得到交会图。
所述的不同的试油结论包括含油水层、水层、油层、气层和油气同层。
相对于现有技术，本发明的有益效果为：
本发明提供的基于阵列感应测井的碎屑岩储层流体识别方法，先获取经试油验证的储层中不同层位的阵列感应测井电阻率值，然后计算出阵列感应斜率特征参数及相关特征参数，再建立交会图并确定出油、气、水的界限区域，最后计算新井的阵列感应斜率特征参数及相关特征参数并投射到交会图上，根据油、气、水的界限区域判断出该新井中的储层流体性质。该方法能够针对低阻油气层快速判别流体性质，在不考虑储层电阻率绝对值大小的情况下，利用阵列感应测井电阻率值得到交会图，以反映井内泥浆在含有不同流体的储层的径向侵入所导致的电阻率变化的不同特征，从而能够快速、准确、连续、可靠的判断出目标储层流体的性质，为现场测井快速解释油气判识提供了一种简单有效的方法。
附图说明
为了更清楚地说明本发明实施例中的技术方案，下面将对实施例描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1是本发明提供的基于阵列感应测井的碎屑岩储层流体识别方法的流程图；
图2是本发明实施例中提供的阵列感应斜率特征值与相关特征值的交会图；
图3是本发明实施例中提供的某一口井的流体识别测井解释图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚，下面将结合附图对本发明作进一步地详细描述。
本发明提供的基于阵列感应测井的碎屑岩储层流体识别方法，实现了在不考虑储层电阻率绝对值大小的情况下，通过现有的阵列感应测井资料，反映泥浆在含有不同流体的储层径向上的侵入致使电阻率变化特征不同，来快速、准确、连续、可靠的判断储层流体的性质(要求泥浆电阻率大于等于两倍地层水电阻率)。
本发明提供的基于阵列感应测井的碎屑岩储层流体识别方法，具体包括以下步骤：
1)利用区域已试油的储层，选择具有试油结论与阵列感应测井资料的井，读取不同试油层位的阵列感应测井电阻率值；
其中对于一口井，其阵列感应测井电阻率值和阵列感应测井资料可由现有的阵列感应测井仪器直接获得。
2)利用步骤1)中的阵列感应测井电阻率值计算阵列感应斜率特征参数k(无量纲)与相关特征参数r(无量纲)；
公式(2)、(3)中：
n为阵列感应测井径向上不同探测深度曲线数；
yi为不同探测深度的电阻率，
xi为不同探测深度电阻率对应的横坐标，取自然数，
为xi的平均值，为yi的平均值，
n、yi、xi可通过阵列感应测井仪器直接获得。
其中根据不同的阵列感应测井仪器得到的阵列感应测井径向上不同探测深度曲线数n为5条或6条。
以本发明实施例中具体使用的阵列感应测井仪器为例，其测得的n＝5，此时yi和xi中的i为1～5的正整数，y1～y5分别为10英寸、20英寸、30英寸、60英寸、90英寸探测深度的电阻率值，x1～x5分别为1、2、3、4、5。
3)以不同的试油结论(包括含油水层、水层、油层、气层和油气同层)为系列，以计算的k为纵坐标，r为横坐标，绘制阵列感应斜率特征参数与相关特征参数的交会图；
4)在步骤3)所建立的交会图上确定出油、气、水的各自区域以及油、气、水的相互界限；
5)新井处理，对于具有阵列感应测井资料的新井，利用公式(2)和(3)计算目的储层的斜率特征参数k与相关特征参数r，然后将新井目的层位k值与r值数据点投到交会图上，利用在交会图上所确定的油、气、水界限区域快速判识目的层位的储层流体的性质。
下面结合具体实施例对本发明作进一步详细说明。
参见图1，本发明提供的基于阵列感应测井的碎屑岩储层流体识别方法，按照以下步骤进行操作：
步骤一，选取区域内具有阵列感应测井资料的井，在已经试油验证的层位上读取层内具有代表性的阵列感应测井电阻率值(包括10英寸、20英寸、30英寸、60英寸、90英寸探测深度的电阻率值)，要求选择的层位内包含有油层、气层、水层。
步骤二，利用步骤一所获取的不同探测深度的阵列感应电阻率数据，按照公式(2)和(3)分别求取每一个层位的阵列感应斜率特征参数k与阵列感应相关特征参数r，具体计算结果如表1所示。
其中AT10为10英寸探测深度的电阻率值、AT20为20英寸探测深度的电阻率值、AT30为30英寸探测深度的电阻率值、AT60为60英寸探测深度的电阻率值、AT90为90英寸探测深度的电阻率值，即AT10、AT20、AT30、AT60和AT90分别为y1～y5的值。
表1
步骤三，利用步骤二计算出来的阵列感应斜率特征参数k与相关特征参数r，根据试油结论分系列绘制出交会图，如图2所示，根据交会图确定出油、气、水各自的区域以及油、气、水界限。
步骤四，对于一口具有阵列感应测井资料的新井，利用公式(2)和(3)计算储层剖面连续的k与r值，将其投射到交会图中，结合交会图中划分出的油、气、水界限区域，快速得出目的层位流体的性质。
如图3中，1号层落入油层区域，2、3、4号层落入水层区域，与表1中的结果相符合。说明本发明提供的储层流体识别方法能够准确、可靠的判识储层流体的特性，在低阻油气储层中有很好的应用效果。
上述实施例序号仅仅为了描述，不代表实施例的优劣。
以上所述仅为本发明的较佳实施例，并不用以限制本发明，凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
Claims (7)
1.一种基于阵列感应测井的碎屑岩储层流体识别方法，其特征在于，包括以下步骤：
1)选择区域内经过试油验证的储层，获取不同试油层位的阵列感应测井电阻率值；
2)根据步骤1)中得到的阵列感应测井电阻率值，计算得到阵列感应斜率特征参数及相关特征参数；
所述阵列感应斜率特征参数通过以下公式计算得到：
<mrow>
<mi>k</mi>
<mo>=</mo>
<mfrac>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mn>1</mn>
<mi>n</mi>
</munderover>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<mover>
<mi>x</mi>
<mo>&OverBar;</mo>
</mover>
<mo>)</mo>
</mrow>
<mrow>
<mo>(</mo>
<msub>
<mi>y</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<mover>
<mi>y</mi>
<mo>&OverBar;</mo>
</mover>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mn>1</mn>
<mi>n</mi>
</munderover>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<mover>
<mi>x</mi>
<mo>&OverBar;</mo>
</mover>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
</mrow>
其中，k为阵列感应斜率特征参数；
n为阵列感应测井径向上的不同探测深度曲线数；
yi为不同探测深度的电阻率；
xi为不同探测深度电阻率对应的横坐标；
为xi的平均值，为yi的平均值；
3)根据步骤2)得到的阵列感应斜率特征参数及相关特征参数，建立阵列感应斜率特征参数与相关特征参数的交会图，并在交会图上确定出油、气、水的界限区域；
4)对于待识别储层流体性质的新井，先计算其阵列感应斜率特征参数及相关特征参数，并将计算结果投射到交会图上，再按照步骤3)得到的油、气、水的界限区域识别该井中储层流体的性质。
2.根据权利要求1所述的基于阵列感应测井的碎屑岩储层流体识别方法，其特征在于：在进行步骤1)时，收集区域内具有阵列感应测井资料并且经过试油验证的储层，获取具有代表性的不同层位的阵列感应测井电阻率值。
3.根据权利要求2所述的基于阵列感应测井的碎屑岩储层流体识别方法，其特征在于：所述的具有代表性的不同层位是指包含有油层、气层、水层的若干层位。
4.根据权利要求1所述的基于阵列感应测井的碎屑岩储层流体识别方法，其特征在于：所述相关特征参数通过以下公式计算得到：
<mrow>
<mi>r</mi>
<mo>=</mo>
<mfrac>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mn>1</mn>
<mi>n</mi>
</munderover>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>x</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<mover>
<mi>x</mi>
<mo>&OverBar;</mo>
</mover>
</mrow>
<mo>)</mo>
</mrow>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>y</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<mover>
<mi>y</mi>
<mo>&OverBar;</mo>
</mover>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
<msqrt>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mn>1</mn>
<mi>n</mi>
</munderover>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>x</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<mover>
<mi>x</mi>
<mo>&OverBar;</mo>
</mover>
</mrow>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<munderover>
<mo>&Sigma;</mo>
<mn>1</mn>
<mi>n</mi>
</munderover>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>y</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<mover>
<mi>y</mi>
<mo>&OverBar;</mo>
</mover>
</mrow>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
</mfrac>
</mrow>
其中，r为相关特征参数；
n为阵列感应测井径向上的不同探测深度曲线数；
yi为不同探测深度的电阻率；
xi为不同探测深度电阻率对应的横坐标；
为xi的平均值，为yi的平均值。
5.根据权利要求1或4所述的基于阵列感应测井的碎屑岩储层流体识别方法，其特征在于：所述n＝5，yi和xi中的i为1～5的正整数，y1～y5分别为10英寸、20英寸、30英寸、60英寸和90英寸探测深度的电阻率值，x1～x5分别为1、2、3、4、5。
6.根据权利要求1所述的基于阵列感应测井的碎屑岩储层流体识别方法，其特征在于：所述步骤3)在建立交会图时，以不同的试油结论为系列，以步骤2)计算的阵列感应斜率特征参数为纵坐标，以相关特征参数为横坐标，绘制得到交会图。
7.根据权利要求6所述的基于阵列感应测井的碎屑岩储层流体识别方法，其特征在于：所述的不同的试油结论包括含油水层、水层、油层、气层和油气同层。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510246298.9A CN104879126B (zh) | 2015-05-14 | 2015-05-14 | 一种基于阵列感应测井的碎屑岩储层流体识别方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510246298.9A CN104879126B (zh) | 2015-05-14 | 2015-05-14 | 一种基于阵列感应测井的碎屑岩储层流体识别方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN104879126A CN104879126A (zh) | 2015-09-02 |
| CN104879126B true CN104879126B (zh) | 2017-09-26 |
Family
ID=53946777
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510246298.9A Active CN104879126B (zh) | 2015-05-14 | 2015-05-14 | 一种基于阵列感应测井的碎屑岩储层流体识别方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104879126B (zh) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107832574B (zh) * | 2017-09-19 | 2020-07-10 | 中国石油天然气集团公司 | 一种基于随钻测井的水平井水淹层解释方法 |
| CN110133724B (zh) * | 2019-05-06 | 2021-01-29 | 中国石油天然气集团有限公司 | 一种应用阵列声波测井数据进行气层识别的方法 |
Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2004083898A1 (en) * | 2003-03-17 | 2004-09-30 | Electromagnetic Geoservices As | Method and apparatus for determining the nature of submarine reservoirs |
| CN101196460A (zh) * | 2007-10-26 | 2008-06-11 | 辽河石油勘探局 | 一种岩石润湿性的评价方法 |
| CN102789003A (zh) * | 2012-08-15 | 2012-11-21 | 中国石油天然气股份有限公司 | 利用阵列感应测井资料进行地层参数反演的方法及装置 |
| CN102830441A (zh) * | 2012-05-16 | 2012-12-19 | 陶冶 | 包含数据挖掘的油气层综合解释评价方法 |
| CN102979517A (zh) * | 2012-12-04 | 2013-03-20 | 中国海洋石油总公司 | 一种复杂油气储层饱和度定量评价方法 |
| CN103132993A (zh) * | 2013-03-26 | 2013-06-05 | 西北大学 | 一种逐步识别低渗复杂岩性油藏中的油层和水层的方法 |
| CN103775057A (zh) * | 2013-12-27 | 2014-05-07 | 中国石油天然气股份有限公司 | 一种致密油气藏有效储层的判识方法及装置 |
| CN103883318A (zh) * | 2014-03-05 | 2014-06-25 | 中国石油天然气股份有限公司 | 一种中低矿化度地层水储层的参数反演方法及装置 |
-
2015
- 2015-05-14 CN CN201510246298.9A patent/CN104879126B/zh active Active
Patent Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2004083898A1 (en) * | 2003-03-17 | 2004-09-30 | Electromagnetic Geoservices As | Method and apparatus for determining the nature of submarine reservoirs |
| CN101196460A (zh) * | 2007-10-26 | 2008-06-11 | 辽河石油勘探局 | 一种岩石润湿性的评价方法 |
| CN102830441A (zh) * | 2012-05-16 | 2012-12-19 | 陶冶 | 包含数据挖掘的油气层综合解释评价方法 |
| CN102789003A (zh) * | 2012-08-15 | 2012-11-21 | 中国石油天然气股份有限公司 | 利用阵列感应测井资料进行地层参数反演的方法及装置 |
| CN102979517A (zh) * | 2012-12-04 | 2013-03-20 | 中国海洋石油总公司 | 一种复杂油气储层饱和度定量评价方法 |
| CN103132993A (zh) * | 2013-03-26 | 2013-06-05 | 西北大学 | 一种逐步识别低渗复杂岩性油藏中的油层和水层的方法 |
| CN103775057A (zh) * | 2013-12-27 | 2014-05-07 | 中国石油天然气股份有限公司 | 一种致密油气藏有效储层的判识方法及装置 |
| CN103883318A (zh) * | 2014-03-05 | 2014-06-25 | 中国石油天然气股份有限公司 | 一种中低矿化度地层水储层的参数反演方法及装置 |
Non-Patent Citations (1)
| Title |
|---|
| 川西E区须家河组致密碎屑岩储层流体识别方法及优选——以E101井为例;蒋健美等;《国外测井技术》;20130430(第2期);34-36 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN104879126A (zh) | 2015-09-02 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104278991B (zh) | 盐湖相烃源岩有机碳和生烃潜量的多元测井计算方法 | |
| CN109838230B (zh) | 油藏水淹层的定量评价方法 | |
| CN108252709B (zh) | 一种致密砂岩油藏的油水性质识别方法及系统 | |
| CN105426612A (zh) | 一种地层组分最优化确定方法及装置 | |
| WO2019232224A1 (en) | Systems and methods for special core analysis sample selection and assessment | |
| CN105447762B (zh) | 一种流体替换的低渗透油藏水淹信息的计算方法 | |
| CN107832574B (zh) | 一种基于随钻测井的水平井水淹层解释方法 | |
| CN109117505A (zh) | 一种基于介电实验的孔隙结构储层冲洗带含水饱和度计算方法 | |
| CN110782187B (zh) | 基于含水率的枯竭油气藏型储气库动态评价方法及系统 | |
| Warren et al. | An evaluation of the significance of permeability measurements | |
| CN104879126B (zh) | 一种基于阵列感应测井的碎屑岩储层流体识别方法 | |
| CN107092995A (zh) | 地热资源评价方法和装置 | |
| CN104514553B (zh) | 一种岩屑砂岩层中气层的识别方法及其应用 | |
| Neves et al. | Multi-attribute seismic volume facies classification for predicting fractures in carbonate reservoirs | |
| AU2017279838B1 (en) | Method for classifying deep rock geofacies based on data mining | |
| CN104636592A (zh) | 通过常规测井获取油层指示标志的方法 | |
| CN111350497A (zh) | 一种基于含水率的低饱和度油层动态评价方法 | |
| CN113775326A (zh) | 可动水饱和度的评价方法、装置、电子设备及介质 | |
| CN104834934B (zh) | 一种用于识别储层流体的核体俘获法 | |
| CN110873904B (zh) | 流体识别方法及装置 | |
| Rafiee et al. | A new cementation factor correlation in carbonate parts of oil fields in south-west Iran | |
| CN109973068B (zh) | 油藏注水诱导裂缝的识别方法及装置 | |
| Mogbolu | Identification and quantification of thinly bedded low resistivity pay in the Niger Delta | |
| Ding et al. | Identification and characterization of high-permeability zones in waterflooding reservoirs with an ensemble of methodologies | |
| Holmes et al. | Mixed reservoir wetting in unconventional reservoirs and interpretation of porosity/resistivity cross plots, derived from triple-combo log data |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| EXSB | Decision made by sipo to initiate substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |