CN104991280B - 一种基于微观孔喉结构指数的射孔层段选择方法及装置 - Google Patents
一种基于微观孔喉结构指数的射孔层段选择方法及装置 Download PDFInfo
- Publication number
- CN104991280B CN104991280B CN201510382782.4A CN201510382782A CN104991280B CN 104991280 B CN104991280 B CN 104991280B CN 201510382782 A CN201510382782 A CN 201510382782A CN 104991280 B CN104991280 B CN 104991280B
- Authority
- CN
- China
- Prior art keywords
- pore throat
- throat radius
- mrow
- distributed data
- magnetic resonance
- 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
- 239000011148 porous material Substances 0.000 title claims abstract description 179
- 238000010187 selection method Methods 0.000 title abstract 2
- 238000005481 NMR spectroscopy Methods 0.000 claims abstract description 61
- 238000009826 distribution Methods 0.000 claims abstract description 61
- 230000015572 biosynthetic process Effects 0.000 claims abstract description 29
- 238000000034 method Methods 0.000 claims abstract description 27
- 239000011435 rock Substances 0.000 claims abstract description 25
- 238000011160 research Methods 0.000 claims description 18
- 230000008569 process Effects 0.000 claims description 9
- 239000004744 fabric Substances 0.000 claims description 8
- 238000012545 processing Methods 0.000 claims description 6
- 238000009533 lab test Methods 0.000 claims description 5
- 239000012530 fluid Substances 0.000 claims description 3
- 238000009825 accumulation Methods 0.000 claims 2
- 230000005311 nuclear magnetism Effects 0.000 claims 1
- 238000011156 evaluation Methods 0.000 description 10
- 239000011159 matrix material Substances 0.000 description 7
- 230000008901 benefit Effects 0.000 description 5
- 238000012360 testing method Methods 0.000 description 5
- 230000000694 effects Effects 0.000 description 4
- 239000003208 petroleum Substances 0.000 description 4
- 238000011161 development Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 238000003860 storage Methods 0.000 description 3
- 238000009736 wetting Methods 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 238000005266 casting Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- QSHDDOUJBYECFT-UHFFFAOYSA-N mercury Chemical compound [Hg] QSHDDOUJBYECFT-UHFFFAOYSA-N 0.000 description 2
- 229910052753 mercury Inorganic materials 0.000 description 2
- 235000016838 Pomo dAdamo Nutrition 0.000 description 1
- 244000003138 Pomo dAdamo Species 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000012512 characterization method Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 238000000280 densification Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 230000002349 favourable effect Effects 0.000 description 1
- 210000000867 larynx Anatomy 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000035699 permeability Effects 0.000 description 1
- 230000000704 physical effect Effects 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 238000004451 qualitative analysis Methods 0.000 description 1
- 238000013441 quality evaluation Methods 0.000 description 1
- 238000004445 quantitative analysis Methods 0.000 description 1
- 238000011158 quantitative evaluation Methods 0.000 description 1
- 238000013139 quantization Methods 0.000 description 1
- 238000004335 scaling law Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明涉及一种基于微观孔喉结构指数的射孔层段选择方法及装置,其中,方法包括:建立核磁共振横向弛豫时间分布数据和岩心样本孔喉半径分布数据之间的关系式;选取目标层段的核磁共振横向弛豫时间分布数据;根据目标层段的核磁共振横向弛豫时间分布数据获得地层孔隙度;利用目标层段的核磁共振横向弛豫时间分布数据获得目标层段的孔喉半径分布数据;利用目标层段的孔喉半径分布数据确定目标层段的最大孔喉半径和孔喉半径中值;根据地层孔隙度、最大孔喉半径和孔喉半径中值建立孔喉品质指数;根据孔喉品质指数值对目标层段进行分类,对比不同层段的孔喉品质指数差异,选择值大的处理层段作为射孔层段。
Description
技术领域
本发明涉及油气勘探技术领域,特别涉及一种基于微观孔喉结构指数的射孔层段选择方法及装置。
背景技术
随着世界能源需求的增加,非常规油气资源越来越受到重视。致密油气是非常规油气资源中比较现实的部分,也是未来油气增储上产的重要接替资源。我国致密油气资源分布广泛,潜力巨大,是目前最具现实意义和勘探开发价值的油气资源。致密油气储层孔渗低、非均质性强,致密油气储层的微观孔喉结构与产能密切相关,储层品质测井定量评价以及储层测井精细分类已经成为提高勘探开发效益的关键。因此,开展致密储层微观孔喉结构研究,提出合理有效的致密储层微观孔喉结构评价与分类方法,对于发现有利勘探目标,有效开展储层品质评价,实现致密油气规模效益开发具有重大意义。
致密储层孔隙结构复杂,储层性质差异大,依靠传统的岩石物理分类方法很难有效地进行储层分类。目前通常依据孔隙度、渗透率等物性特征以及铸体薄片、压汞毛管压力曲线得到的孔隙类型和孔隙结构信息,如面孔率、平均孔径、分选系数、孔隙组合类型、排驱压力可动流体饱和度等来建立储层分类标准,但在没有岩心分析资料时,则难以应用该方法对储层进行有效分类。在实际应用中,通常需要采用测井资料进行储层分类。考虑到岩心压汞实验和铸体薄片周期长且不能应用于地层连续评价,常规的以孔隙度为主的测井评价方法在致密油气储层分类中应用效果不理想,当前缺乏应用测井资料开展致密油气储层微观孔喉结构定量评价及在此基础上优选射孔层段的有效方法。核磁共振测井能连续反映储层孔径分布,如何通过核磁共振测井提取储层微观孔喉信息的表征参数并构建储层微观孔喉结构定量评价参数是致密油气储层品质测井评价的重点,以此为基础,实现致密油气储层分类以及射孔层段优选。
发明内容
为解决现有技术的问题,本发明提出一种基于微观孔喉结构指数的射孔层段选择方法及装置,利用核磁共振测井资料计算储层孔喉半径分布,提取微观孔喉信息,构建储层微观孔喉品质指数PTI,实现对致密油气储层微观孔喉结构的定量评价、储层分类与射孔层段优选。
为实现上述目的,本发明提供了一种基于微观孔喉结构指数的射孔层段选择方法,包括:
确定岩心样本的核磁共振横向弛豫时间分布数据T2和孔喉半径分布数据,并建立核磁共振横向弛豫时间分布数据T2和岩心样本孔喉半径分布数据之间的关系式;
选取目标层段的核磁共振横向弛豫时间分布数据T2;
根据目标层段的核磁共振横向弛豫时间分布数据T2获得地层孔隙度φ;
根据核磁共振横向弛豫时间分布数据T2和岩心样本孔喉半径分布数据之间的关系式,利用目标层段的核磁共振横向弛豫时间分布数据T2获得目标层段的孔喉半径分布数据;
利用目标层段的孔喉半径分布数据确定目标层段的最大孔喉半径和孔喉半径中值;
根据地层孔隙度φ、最大孔喉半径和孔喉半径中值建立孔喉品质指数PTI;
根据孔喉品质指数PTI值对目标层段进行分类,对比不同层段的孔喉品质指数PTI差异,选择PTI值大的处理层段作为射孔层段。
优选地,所述核磁共振横向弛豫时间分布数据T2和岩心样本孔喉半径分布数据之间的关系式为:
其中,A、B、C、D、E为实验室配套实验资料刻度确定;T2为核磁共振横向弛豫时间分布数据,Rpt为孔喉半径分布数据。
优选地,所述地层孔隙度φ的表达式为:
其中,BIN为孔隙度分量个数,即核磁共振T2分布的布点数,Ami为T2分布第i个分量的幅度值。
优选地,所述孔喉品质指数PTI的表达式为:
PTI=ω1*f1(Rp max)+ω2*f2(Rp50)+ω3*f3(φ)
其中,w1、w2、w3均为权重系数,通过层次分析法确定;Rpmax表示最大孔喉半径,为非润湿相流体进入孔隙网络时最先突入的喉道值;Rp50表示孔喉半径中值,为孔喉半径分布累积频率为50%处所对应的孔喉半径值;f1、f2、f3均为孔喉参数的归一化函数;
min(Rp max)表示研究区最大孔喉半径分布的极小值;max(Rp max)表示研究区最大孔喉半径分布的极大值;min(Rp50)表示研究区孔喉半径中值分布的极小值;max(Rp50)表示研究区孔喉半径中值分布的极大值;min(φ)表示研究区孔隙度分布的极小值;max(φ)表示研究区孔隙度分布的极大值。
为实现上述目的,本发明还提供了一种基于微观孔喉结构指数的射孔层段选择装置,包括:
关系式建立单元,用于确定岩心样本的核磁共振横向弛豫时间分布数据T2和孔喉半径分布数据,并建立核磁共振横向弛豫时间分布数据T2和岩心样本孔喉半径分布数据之间的关系式;
核磁共振横向弛豫时间分布数据选取单元,用于选取目标层段的核磁共振横向弛豫时间分布数据T2;
地层孔隙度获取单元,用于根据目标层段的核磁共振横向弛豫时间分布数据T2获得地层孔隙度φ;
孔喉半径分布数据获取单元,用于根据核磁共振横向弛豫时间分布数据T2和岩心样本孔喉半径分布数据之间的关系式,利用目标层段的核磁共振横向弛豫时间分布数据T2获得目标层段的孔喉半径分布数据;
最大孔喉半径和孔喉半径中值获取单元,用于利用目标层段的孔喉半径分布数据确定目标层段的最大孔喉半径和孔喉半径中值;
孔喉品质指数建立单元,用于根据地层孔隙度φ、最大孔喉半径和孔喉半径中值建立孔喉品质指数PTI;
射孔层段选择单元,用于根据孔喉品质指数PTI值对目标层段进行分类,对比不同层段的孔喉品质指数差异,选择孔喉品质指数差异最大的处理层段作为射孔层段。
优选地,所述关系式建立单元建立的核磁共振横向弛豫时间分布数据T2和岩心样本孔喉半径分布数据之间的关系式为:
其中,A、B、C、D、E为实验室配套实验资料刻度确定;T2为核磁共振横向弛豫时间分布数据,Rpt为孔喉半径分布数据。
优选地,所述地层孔隙度获取单元得到的地层孔隙度φ的表达式为:
其中,BIN为孔隙度分量个数,即核磁共振T2分布的布点数,Ami为T2分布第i个分量的幅度值。
优选地,所述孔喉品质指数建立单元获得的孔喉品质指数PTI的表达式为:
PTI=ω1*f1(Rp max)+ω2*f2(Rp50)+ω3*f3(φ)
其中,w1、w2、w3均为权重系数,通过层次分析法确定;Rpmax表示最大孔喉半径,为非润湿相流体进入孔隙网络时最先突入的喉道值;Rp50表示孔喉半径中值,为孔喉半径分布累积频率为50%处所对应的孔喉半径值;f1、f2、f3均为孔喉参数的归一化函数;
min(Rp max)表示研究区最大孔喉半径分布的极小值;max(Rp max)表示研究区最大孔喉半径分布的极大值;min(Rp50)表示研究区孔喉半径中值分布的极小值;max(Rp50)表示研究区孔喉半径中值分布的极大值;min(φ)表示研究区孔隙度分布的极小值;max(φ)表示研究区孔隙度分布的极大值。
上述技术方案具有如下有益效果:本技术方案根据核磁共振横向弛豫时间分布数据计算孔喉半径分布以及地层孔隙度,构建表征储层微观孔喉结构的孔喉品质指数,利用孔喉品质指数实现储层微观孔喉结构参数的定量计算、储层孔喉品质指数定量计算和射孔层段优选,对优选层段进行射孔试油。具有简单、直观、区分度高、可靠性好等优点,具有明显的实际应用效果。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明提出的一种基于微观孔喉结构指数的射孔层段选择方法流程图;
图2是利用本实施例1分别对两口井的储层微观孔喉品质指数计算结果对比图与射孔层段选择实例;
图3是本发明实施例1提供的一种基于微观孔喉结构指数的射孔层段选择装置框图;
图4是本发明实施例2提供的储层微观孔喉品质指数计算模块结构示意图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
如图1所示,为本发明提出的一种基于微观孔喉结构指数的射孔层段选择方法流程图。包括:
步骤101):确定岩心样本的核磁共振横向弛豫时间分布数据T2和孔喉半径分布数据,并建立核磁共振横向弛豫时间分布数据T2和岩心样本孔喉半径分布数据之间的关系式;
步骤102):选取目标层段的核磁共振横向弛豫时间分布数据T2;
步骤103):根据目标层段的核磁共振横向弛豫时间分布数据T2获得地层孔隙度φ;
在步骤103中,所述地层孔隙度φ的表达式为:
其中,BIN为孔隙度分量个数,即核磁共振T2分布的布点数,Ami为T2分布第i个分量的幅度值。
步骤104):根据核磁共振横向弛豫时间分布数据T2和岩心样本孔喉半径分布数据之间的关系式,利用目标层段的核磁共振横向弛豫时间分布数据T2获得目标层段的孔喉半径分布数据;
在步骤104中,所述核磁共振横向弛豫时间分布数据T2和岩心样本孔喉半径分布数据之间的关系式为:
其中,A、B、C、D、E为实验室配套实验资料刻度确定;T2为核磁共振横向弛豫时间分布数据,Rpt为孔喉半径分布数据。
步骤105):利用目标层段的孔喉半径分布数据确定目标层段的最大孔喉半径和孔喉半径中值;
步骤106):根据地层孔隙度φ、最大孔喉半径和孔喉半径中值建立孔喉品质指数PTI;
在步骤106中,所述孔喉品质指数PTI的表达式为:
PTI=ω1*f1(Rp max)+ω2*f2(Rp50)+ω3*f3(φ)
其中,w1、w2、w3均为权重系数,通过层次分析法确定;Rpmax表示最大孔喉半径,为非润湿相流体进入孔隙网络时最先突入的喉道值;Rp50表示孔喉半径中值,为孔喉半径分布累积频率为50%处所对应的孔喉半径值;f1、f2、f3均为孔喉参数的归一化函数;
min(Rp max)表示研究区最大孔喉半径分布的极小值;max(Rp max)表示研究区最大孔喉半径分布的极大值;min(Rp50)表示研究区孔喉半径中值分布的极小值;max(Rp50)表示研究区孔喉半径中值分布的极大值;min(φ)表示研究区孔隙度分布的极小值;max(φ)表示研究区孔隙度分布的极大值。
对于本技术方案来说,步骤106具体包括:
S1:计算权重系数w1、w2、w3。
采用层次分析法,确定权重系数w1、w2、w3。
层次分析法又称AHP构权法(Analytic hierarchy process,简写为AHP),是将复杂的评价对象排列为一个有序的递阶层次结构的整体,然后在各个评价项目之间进行两两的比较、判断,计算各个评价项目的相对重要性系数,即权重。AHP构权法又分为单准则构权法和多准则构权法,在此介绍单准则构权法及具体步骤。
(1)确定指标的量化标准
层次分析法的核心问题是建立一个构造合理且一致的判断矩阵,判断矩阵的合理性受到标度的合理性的影响。所谓标度是指评价者对各个评价指标(或者项目)重要性等级差异的量化概念。确定指标重要性的量化标准常用的方法有:比例标度法和指数标度法。比例标度法是以对事物质的差别的评判标准为基础,一般以5种判别等级表示事物质的差别。当评价分析需要更高的精确度时,可以使用9种判别等级来评价,见表1。
表1 比例标度值体系表(重要性分数xij)
(2)确定初始权数
初始权数的确定常常采用定性分析和定量分析相结合的方法。具体操作为:根据分析研究的目的、已经建立的评价指标体系和初步确定的指标重要性的量化标准给出各个指标的权重,分别计算各个指标权重的平均数和标准差,通过不断优化得到初始权数。
(3)对初始权数进行处理
第一步,建立判断矩阵A。通过对评价指标的评价,进行两两比较,其初始权数形成判断矩阵A,见下表。判断矩阵A中第i行和第j列的元素xij表示指标xi与xj比较后所得的标度系数。
表2 参数的判断矩阵A
指标 | Rpmax | Rp50 | φ |
Rpmax | 1 | 1/5 | 1/3 |
Rp50 | 5 | 1 | 3 |
φ | 3 | 1/3 | 1 |
第二步,计算判断矩阵A中的每一行各标度数据的几何平均数,记作wi。
第三步,进行归一化处理。归一化处理是利用公式计算,依据计算结果确定各个指标的权重系数。
Σwi=w1+w2+w3=3.8717
w1=0.1047
w2=0.6370
w3=0.2583
S2:计算孔喉参数的归一化函数f1、f2、f3。
选取工区储层最小和最大的Rpmax、Rp50、φ,即:min(Rp max)、max(Rp max),min(Rp50)、max(Rp50),min(φ)、max(φ),采用下面的公式得到归一化函数f1、f2、f3:
其中:max(Rp max)=0.75、min(Rp max)=0.3、max(Rp50)=0.25、min(Rp50)=0.1、max(φ)=12、min(φ)=4。
S3:计算孔喉品质指数PTI。
根据所述最大孔喉半径Rpmax、孔喉半径中值Rp50和孔隙度φ的权重系数以及归一化函数建立孔喉品质指数PTI,其表达式为:
PTI=ω1*f1(Rp max)+ω2*f2(Rp50)+ω3*f3(φ)
步骤107):根据孔喉品质指数PTI值对目标层段进行分类,对比不同层段的孔喉品质指数PTI差异,选择PTI值大的处理层段作为射孔层段。
如图4所示,为本发明提出的一种基于微观孔喉结构指数的射孔层段选择装置框图。包括:
关系式建立单元401,用于确定岩心样本的核磁共振横向弛豫时间分布数据T2和孔喉半径分布数据,并建立核磁共振横向弛豫时间分布数据T2和岩心样本孔喉半径分布数据之间的关系式;
核磁共振横向弛豫时间分布数据选取单元402,用于选取目标层段的核磁共振横向弛豫时间分布数据T2;
地层孔隙度获取单元403,用于根据目标层段的核磁共振横向弛豫时间分布数据T2获得地层孔隙度φ;
孔喉半径分布数据获取单元404,用于根据核磁共振横向弛豫时间分布数据T2和岩心样本孔喉半径分布数据之间的关系式,利用目标层段的核磁共振横向弛豫时间分布数据T2获得目标层段的孔喉半径分布数据;
最大孔喉半径和孔喉半径中值获取单元405,用于利用目标层段的孔喉半径分布数据确定目标层段的最大孔喉半径和孔喉半径中值;
孔喉品质指数建立单元406,用于根据地层孔隙度φ、最大孔喉半径和孔喉半径中值建立孔喉品质指数PTI;
射孔层段选择单元407,用于根据孔喉品质指数PTI值对目标层段进行分类,对比不同层段的孔喉品质指数差异,选择孔喉品质指数差异最大的处理层段作为射孔层段。
对于本实施例来说,利用本技术方案获得的射孔层段进行射孔试油来验证本技术方案的优越性。
通常情况下,所选层段计算的PTI值越大,说明储层孔喉品质越好,储层品质越好。因此,根据PTI值,对储层进行分类(表3),并优选Ⅰ、Ⅱ类储层作为射孔层段试油开发。
表3 基于储层微观孔喉品质指数PTI的储层分类表
如图2所示,为利用本技术方案获得的A1井的不同层段的孔喉品质指数结果图。如图3所示,为利用本技术方案获得的A2井的不同层段的孔喉品质指数结果图。根据计算结果,A1井45号层、46号层的PTI值分别为0.24~0.57(平均值0.39)、0.14~0.46(平均值0.32),测井评价结果为差储层(Ⅳ类)。A2井104号层、105号层、106号层的PTI值分布为0.57~0.87(平均值0.73)、0.29~0.66(平均值0.51)、0.51~1.00(平均值0.75),测井评价104号层、106号层为较好储层(Ⅱ类),说明A2井104号层、106号层储层品质较好,故优选这两层段进行射孔试油,试油结果为日产油13.35吨,产能较好。
应用本发明方法实现了对储层进行有效分类和射孔层段优选,具有简单、直观、区分度高、可靠性好等优点,具有明显的实际应用效果。
以上所述的具体实施方式,对本发明的目的、技术方案和有益效果进行了进一步详细说明,所应理解的是,以上所述仅为本发明的具体实施方式而已,并不用于限定本发明的保护范围,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (6)
1.一种基于微观孔喉结构指数的射孔层段选择方法,其特征在于,包括:
确定岩心样本的核磁共振横向弛豫时间分布数据T2和孔喉半径分布数据,并建立核磁共振横向弛豫时间分布数据T2和岩心样本孔喉半径分布数据之间的关系式;
选取目标层段的核磁共振横向弛豫时间分布数据T2;
根据目标层段的核磁共振横向弛豫时间分布数据T2获得地层孔隙度φ;
根据核磁共振横向弛豫时间分布数据T2和岩心样本孔喉半径分布数据之间的关系式,利用目标层段的核磁共振横向弛豫时间分布数据T2获得目标层段的孔喉半径分布数据;
利用目标层段的孔喉半径分布数据确定目标层段的最大孔喉半径和孔喉半径中值;
根据地层孔隙度φ、最大孔喉半径和孔喉半径中值建立孔喉品质指数PTI;
根据孔喉品质指数PTI值对目标层段进行分类,对比不同层段的孔喉品质指数PTI差异,选择PTI值大的处理层段作为射孔层段;
其中,所述孔喉品质指数PTI的表达式为:
PTI=w1*f1(Rpmax)+w2*f2(Rp50)+w3*f3(φ)
其中,w1、w2、w3均为权重系数,通过层次分析法确定;Rpmax表示最大孔喉半径,为非润湿相流体进入孔隙网络时最先突入的喉道值;Rp50表示孔喉半径中值,为孔喉半径分布累积频率为50%处所对应的孔喉半径值;f1、f2、f3均为孔喉参数的归一化函数;
min(Rpmax)表示研究区最大孔喉半径分布的极小值;max(Rpmax)表示研究区最大孔喉半径分布的极大值;min(Rp50)表示研究区孔喉半径中值分布的极小值;max(Rp50)表示研究区孔喉半径中值分布的极大值;min(φ)表示研究区孔隙度分布的极小值;max(φ)表示研究区孔隙度分布的极大值。
2.如权利要求1所述的方法,其特征在于,所述核磁共振横向弛豫时间分布数据T2和岩心样本孔喉半径分布数据之间的关系式为:
<mrow>
<msub>
<mi>R</mi>
<mrow>
<mi>p</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msubsup>
<mi>ET</mi>
<mn>2</mn>
<mi>D</mi>
</msubsup>
</mrow>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>+</mo>
<mfrac>
<mi>A</mi>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>BT</mi>
<mn>2</mn>
</msub>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mi>C</mi>
</msup>
</mfrac>
<mo>)</mo>
</mrow>
</mfrac>
</mrow>
其中,A、B、C、D、E为实验室配套实验资料刻度确定;T2为核磁共振横向弛豫时间分布数据,Rpt为孔喉半径分布数据。
3.如权利要求1所述的方法,其特征在于,所述地层孔隙度φ的表达式为:
<mrow>
<mi>&phi;</mi>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mrow>
<mi>B</mi>
<mi>I</mi>
<mi>N</mi>
</mrow>
</munderover>
<msub>
<mi>Am</mi>
<mi>i</mi>
</msub>
</mrow>
其中,BIN为孔隙度分量个数,即核磁共振T2分布的布点数,Ami为T2分布第i个分量的幅度值。
4.一种基于微观孔喉结构指数的射孔层段选择装置,其特征在于,包括:
关系式建立单元,用于确定岩心样本的核磁共振横向弛豫时间分布数据T2和孔喉半径分布数据,并建立核磁共振横向弛豫时间分布数据T2和岩心样本孔喉半径分布数据之间的关系式;
核磁共振横向弛豫时间分布数据选取单元,用于选取目标层段的核磁共振横向弛豫时间分布数据T2;
地层孔隙度获取单元,用于根据目标层段的核磁共振横向弛豫时间分布数据T2获得地层孔隙度φ;
孔喉半径分布数据获取单元,用于根据核磁共振横向弛豫时间分布数据T2和岩心样本孔喉半径分布数据之间的关系式,利用目标层段的核磁共振横向弛豫时间分布数据T2获得目标层段的孔喉半径分布数据;
最大孔喉半径和孔喉半径中值获取单元,用于利用目标层段的孔喉半径分布数据确定目标层段的最大孔喉半径和孔喉半径中值;
孔喉品质指数建立单元,用于根据地层孔隙度φ、最大孔喉半径和孔喉半径中值建立孔喉品质指数PTI;
射孔层段选择单元,用于根据孔喉品质指数PTI值对目标层段进行分类,对比不同层段的孔喉品质指数差异,选择孔喉品质指数差异最大的处理层段作为射孔层段;
其中,所述孔喉品质指数建立单元获得的孔喉品质指数PTI的表达式为:
PTI=w1*f1(Rpmax)+w2*f2(Rp50)+w3*f3(φ)
w1、w2、w3均为权重系数,通过层次分析法确定;Rpmax表示最大孔喉半径,为非润湿相流体进入孔隙网络时最先突入的喉道值;Rp50表示孔喉半径中值,为孔喉半径分布累积频率为50%处所对应的孔喉半径值;f1、f2、f3均为孔喉参数的归一化函数;
min(Rpmax)表示研究区最大孔喉半径分布的极小值;max(Rpmax)表示研究区最大孔喉半径分布的极大值;min(Rp50)表示研究区孔喉半径中值分布的极小值;max(Rp50)表示研究区孔喉半径中值分布的极大值;min(φ)表示研究区孔隙度分布的极小值;max(φ)表示研究区孔隙度分布的极大值。
5.如权利要求4所述的装置,其特征在于,所述关系式建立单元建立的核磁共振横向弛豫时间分布数据T2和岩心样本孔喉半径分布数据之间的关系式为:
<mrow>
<msub>
<mi>R</mi>
<mrow>
<mi>p</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msubsup>
<mi>ET</mi>
<mn>2</mn>
<mi>D</mi>
</msubsup>
</mrow>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>+</mo>
<mfrac>
<mi>A</mi>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>BT</mi>
<mn>2</mn>
</msub>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mi>C</mi>
</msup>
</mfrac>
<mo>)</mo>
</mrow>
</mfrac>
</mrow>
其中,A、B、C、D、E为实验室配套实验资料刻度确定;T2为核磁共振横向弛豫时间分布数据,Rpt为孔喉半径分布数据。
6.如权利要求4所述的装置,其特征在于,所述地层孔隙度获取单元得到的地层孔隙度φ的表达式为:
<mrow>
<mi>&phi;</mi>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mrow>
<mi>B</mi>
<mi>I</mi>
<mi>N</mi>
</mrow>
</munderover>
<msub>
<mi>Am</mi>
<mi>i</mi>
</msub>
</mrow>
2
其中,BIN为孔隙度分量个数,即核磁共振T2分布的布点数,Ami为T2分布第i个分量的幅度值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510382782.4A CN104991280B (zh) | 2015-07-02 | 2015-07-02 | 一种基于微观孔喉结构指数的射孔层段选择方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510382782.4A CN104991280B (zh) | 2015-07-02 | 2015-07-02 | 一种基于微观孔喉结构指数的射孔层段选择方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104991280A CN104991280A (zh) | 2015-10-21 |
CN104991280B true CN104991280B (zh) | 2017-09-01 |
Family
ID=54303116
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510382782.4A Active CN104991280B (zh) | 2015-07-02 | 2015-07-02 | 一种基于微观孔喉结构指数的射孔层段选择方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104991280B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105785457B (zh) * | 2016-05-13 | 2018-04-06 | 吉林大学 | 基于双极性脉冲的地面核磁共振横向弛豫时间的测量方法 |
CN106526079B (zh) * | 2016-10-27 | 2019-04-16 | 中国石油大学(北京) | 一种研究致密砂岩孔喉结构动态变化的方法 |
CN110244369B (zh) * | 2019-06-28 | 2020-10-13 | 中国石油大学(北京) | 储层束缚和可动流体分布确定方法、装置及*** |
CN110987751B (zh) * | 2019-11-15 | 2021-05-18 | 东北石油大学 | 一种致密储层三维空间孔喉定量分级评价方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102141637A (zh) * | 2010-01-28 | 2011-08-03 | 中国石油天然气股份有限公司 | 一种利用核磁共振测井资料连续定量评价储集层孔隙结构的方法 |
CN104634718A (zh) * | 2015-03-05 | 2015-05-20 | 中国石油大学(华东) | 应用核磁共振表征致密砂岩孔径分布的标定方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6714009B2 (en) * | 2002-05-16 | 2004-03-30 | Schlumberger Technology Corporation | Method for the inversion of CPMG measurements enhanced by often repeated short wait time measurements |
US20140132259A1 (en) * | 2012-11-13 | 2014-05-15 | Schlumberger Technology Corporation | Nmr method to determine grain size distribution in mixed saturation |
-
2015
- 2015-07-02 CN CN201510382782.4A patent/CN104991280B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102141637A (zh) * | 2010-01-28 | 2011-08-03 | 中国石油天然气股份有限公司 | 一种利用核磁共振测井资料连续定量评价储集层孔隙结构的方法 |
CN104634718A (zh) * | 2015-03-05 | 2015-05-20 | 中国石油大学(华东) | 应用核磁共振表征致密砂岩孔径分布的标定方法 |
Non-Patent Citations (1)
Title |
---|
利用核磁共振测井资料计算平均孔喉半径;赵文杰;《油气地质与采收率》;20090331;第16卷(第2期);第43-45页 * |
Also Published As
Publication number | Publication date |
---|---|
CN104991280A (zh) | 2015-10-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104819923B (zh) | 基于核磁共振的低渗透砂岩储层孔隙结构定量反演方法 | |
CN108222925A (zh) | 页岩气储层分级综合评价方法 | |
CN106526693A (zh) | 裂缝识别方法和装置 | |
CN104991280B (zh) | 一种基于微观孔喉结构指数的射孔层段选择方法及装置 | |
CN104278989B (zh) | 一种获取低孔低渗储层饱和度指数的方法 | |
Zhang et al. | Experimental study on the pore structure characteristics of tight sandstone reservoirs in Upper Triassic Ordos Basin China | |
CN105134195A (zh) | 一种基于测井资料的页岩气储层品质评价方法 | |
CN105275456B (zh) | 一种利用测井资料识别优质泥页岩的方法 | |
CN105447762B (zh) | 一种流体替换的低渗透油藏水淹信息的计算方法 | |
CN104806232B (zh) | 一种确定碳酸盐岩储层孔隙度下限的方法 | |
CN104750884B (zh) | 基于多因素非线性回归的页岩油气富集指数定量评价方法 | |
CN106050225A (zh) | 一种核磁共振测井t2谱100%纯水谱的确定方法 | |
CN104564042B (zh) | 一种页岩储层脆性的评价方法 | |
CN104863574B (zh) | 一种适用于致密砂岩储层的流体识别方法 | |
CN104564041A (zh) | 基于开发渗透率下限的低渗透碎屑岩储层有效性评价方法 | |
CN109025982A (zh) | 页岩气开发层段的分级评价方法、装置及终端设备 | |
CN102913233A (zh) | 一种基于无因次对比图版的优势流动通道识别方法 | |
CN107346455A (zh) | 一种用于识别页岩气产能的方法 | |
CN104977613A (zh) | 基于多信息的碳酸盐岩岩相古地理重构方法及装置 | |
Zhu et al. | Rapid identification of high-quality marine shale gas reservoirs based on the oversampling method and random forest algorithm | |
CN104047598A (zh) | 非均质古岩溶碳酸盐岩储层产能预测方法 | |
CN105986813A (zh) | 准致密储层快速评价方法及准致密储层多指标评价方法 | |
CN107102377B (zh) | 定量预测致密砂岩油气有利勘探区的方法 | |
CN103376468A (zh) | 基于神经网络函数逼近算法的储层参数定量表征方法 | |
CN112363242A (zh) | 基于测录井融合的储层流体识别方法与装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |