CN105891089B - 储层渗透率确定方法及装置 - Google Patents
储层渗透率确定方法及装置 Download PDFInfo
- Publication number
- CN105891089B CN105891089B CN201610330618.3A CN201610330618A CN105891089B CN 105891089 B CN105891089 B CN 105891089B CN 201610330618 A CN201610330618 A CN 201610330618A CN 105891089 B CN105891089 B CN 105891089B
- Authority
- CN
- China
- Prior art keywords
- reservoir
- core sample
- magnetic resonance
- nuclear magnetic
- permeability
- 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
- 230000035699 permeability Effects 0.000 title claims abstract description 107
- 238000000034 method Methods 0.000 title claims abstract description 27
- 238000005481 NMR spectroscopy Methods 0.000 claims abstract description 117
- 238000001228 spectrum Methods 0.000 claims abstract description 106
- 238000000491 multivariate analysis Methods 0.000 claims abstract description 47
- 230000005311 nuclear magnetism Effects 0.000 claims description 10
- 238000009826 distribution Methods 0.000 claims description 9
- 238000010606 normalization Methods 0.000 claims description 8
- 230000000694 effects Effects 0.000 abstract description 6
- 230000008901 benefit Effects 0.000 abstract description 5
- 238000010586 diagram Methods 0.000 description 14
- 238000004458 analytical method Methods 0.000 description 8
- 238000004590 computer program Methods 0.000 description 7
- 238000003860 storage Methods 0.000 description 7
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 6
- 238000012545 processing Methods 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 4
- 239000011148 porous material Substances 0.000 description 4
- 238000009738 saturating Methods 0.000 description 4
- 238000011156 evaluation Methods 0.000 description 3
- 230000006870 function Effects 0.000 description 3
- 238000005259 measurement Methods 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 2
- 230000001186 cumulative effect Effects 0.000 description 2
- 238000002592 echocardiography Methods 0.000 description 2
- 239000012530 fluid Substances 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 208000019585 progressive encephalomyelitis with rigidity and myoclonus Diseases 0.000 description 2
- 239000011435 rock Substances 0.000 description 2
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 230000005251 gamma ray Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000001764 infiltration Methods 0.000 description 1
- 230000008595 infiltration Effects 0.000 description 1
- 238000011545 laboratory measurement Methods 0.000 description 1
- 238000012821 model calculation Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
- G01N15/08—Investigating permeability, pore-volume, or surface area of porous materials
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B47/00—Survey of boreholes or wells
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
- G01N15/08—Investigating permeability, pore-volume, or surface area of porous materials
- G01N15/088—Investigating volume, surface area, size or distribution of pores; Porosimetry
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V5/00—Prospecting or detecting by the use of ionising radiation, e.g. of natural or induced radioactivity
- G01V5/04—Prospecting or detecting by the use of ionising radiation, e.g. of natural or induced radioactivity specially adapted for well-logging
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- General Physics & Mathematics (AREA)
- Pathology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- Health & Medical Sciences (AREA)
- Immunology (AREA)
- Dispersion Chemistry (AREA)
- Geophysics (AREA)
- Engineering & Computer Science (AREA)
- Analytical Chemistry (AREA)
- Mining & Mineral Resources (AREA)
- Geology (AREA)
- Fluid Mechanics (AREA)
- Environmental & Geological Engineering (AREA)
- Geochemistry & Mineralogy (AREA)
- High Energy & Nuclear Physics (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
本发明公开了一种储层渗透率确定方法及装置,其中方法包括:选取储层岩心样品,确定储层岩心样品的渗透率及核磁共振回波串;对储层岩心样品的核磁共振回波串进行反演,获得储层岩心样品的核磁共振T2谱,计算储层岩心样品的核磁共振T2谱的特征参数及储层岩心样品的孔隙度;建立储层岩心样品的渗透率与储层岩心样品的核磁共振T2谱的特征参数和储层岩心样品的孔隙度的多元统计关系式;根据所述多元统计关系式,以及储层的T2谱的特征参数及储层的孔隙度,确定储层的渗透率。本发明具有简单、精确、可靠性好等优点,具有明显的实际应用效果,为油藏分类、产能预测和储层建模提供了有力的技术支持。
Description
技术领域
本发明涉及油气勘探中的测井评价技术领域,尤其涉及储层渗透率确定方法及装置。
背景技术
渗透率评价对于油藏分类、产能预测和储层建模具有重要意义。储层渗透率主要受孔隙度、束缚水饱和度和孔隙结构等因素影响。目前通常利用岩心分析渗透率和孔隙度等参数的统计关系建立计算渗透率的地区经验公式。在中高孔渗的常规储层,孔隙结构相对简单,这种经验公式一般能够满足储层渗透率计算精度要求。但该方法并不适用于具有复杂孔隙结构储层渗透率的计算。
核磁共振测井作为一种新的测井技术,所测信号与岩石孔隙流体中氢核的数量有关,能得到孔隙度、孔径分布等信息,对渗透率评价具有独特的优势。目前利用核磁共振测井求取储层渗透率主要依据Coates模型和SDR模型,分别为:
(Coates模型)
(SDR模型)
其中,FFI为自由流体孔隙度;BVI为束缚水孔隙度;T2LM为T2谱的几何均值;φ为孔隙度;c1、c2、m1、m2、n1、n2为模型参数。
这两个模型在求取储层渗透率时均存在一定的局限性。Coates模型确定的渗透率对束缚水饱和度的值很敏感。SDR模型主要用T2谱的几何均值来反映整个T2谱的信息,对水层有比较好的计算效果,而对于油气层,计算结果不可靠。
发明内容
本发明实施例提供一种储层渗透率确定方法,用以获得准确的储层渗透率,该方法包括:
选取储层岩心样品,确定储层岩心样品的渗透率及核磁共振回波串;
对储层岩心样品的核磁共振回波串进行反演,获得储层岩心样品的核磁共振T2(横向弛豫时间)谱,计算储层岩心样品的核磁共振T2谱的特征参数及储层岩心样品的孔隙度;
建立储层岩心样品的渗透率与储层岩心样品的核磁共振T2谱的特征参数和储层岩心样品的孔隙度的多元统计关系式;
根据所述多元统计关系式,以及储层的T2谱的特征参数及储层的孔隙度,确定储层的渗透率;
所述计算储层岩心样品的核磁共振T2谱的特征参数,包括:
将储层岩心样品的核磁共振T2谱从长弛豫组分到短弛豫组分反向累加并归一化,获得储层岩心样品的核磁共振T2谱反向累积归一化分布S;
将S和T2的乘积S×T2与S关系曲线中的最大值max(S×T2)作为储层岩心样品的核磁共振T2谱的特征参数T2,Smax;
所述建立储层岩心样品的渗透率与储层岩心样品的核磁共振T2谱的特征参数和储层岩心样品的孔隙度的多元统计关系式,包括建立如下多元统计关系式:
其中,k为储层岩心样品的渗透率,mD;φ为储层岩心样品的孔隙度,小数;T2,Smax为储层岩心样品的核磁共振T2谱的特征参数;a、b、c为待求的参数,通过拟合得到。
一个实施例中,所述建立储层岩心样品的渗透率与储层岩心样品的核磁共振T2谱的特征参数和储层岩心样品的孔隙度的多元统计关系式,进一步包括:
将多元统计关系式两边取对数,转化为线性表达式:
lg(k)=lg(a)+blg(T2,Smax)+clg(φ);
用最小二乘法拟合所述线性表达式,得到参数a、b、c。
一个实施例中,所述根据所述多元统计关系式,以及储层的T2谱的特征参数及储层的孔隙度,确定储层的渗透率,包括:
确定储层的核磁共振回波串;
对储层的核磁共振回波串进行反演,获得储层的核磁共振T2谱,计算储层的核磁共振T2谱的特征参数及储层的孔隙度;
根据多元统计关系式以及储层的T2谱的特征参数及储层的孔隙度,确定储层的渗透率。
本发明实施例还提供一种储层渗透率确定装置,用以获得准确的储层渗透率,该装置包括:
样品参数确定模块,用于选取储层岩心样品,确定储层岩心样品的渗透率及核磁共振回波串;
样品参数计算模块,用于对储层岩心样品的核磁共振回波串进行反演,获得储层岩心样品的核磁共振T2谱,计算储层岩心样品的核磁共振T2谱的特征参数及储层岩心样品的孔隙度;
关系式建立模块,用于建立储层岩心样品的渗透率与储层岩心样品的核磁共振T2谱的特征参数和储层岩心样品的孔隙度的多元统计关系式;
储层渗透率确定模块,用于根据所述多元统计关系式,以及储层的T2谱的特征参数及储层的孔隙度,确定储层的渗透率;
所述样品参数计算模块具体用于:
将储层岩心样品的核磁共振T2谱从长弛豫组分到短弛豫组分反向累加并归一化,获得储层岩心样品的核磁共振T2谱反向累积归一化分布S;
将S和T2的乘积S×T2与S关系曲线中的最大值max(S×T2)作为储层岩心样品的核磁共振T2谱的特征参数T2,Smax;
所述关系式建立模块具体用于建立如下多元统计关系式:
其中,k为储层岩心样品的渗透率,mD;φ为储层岩心样品的孔隙度,小数;T2,Smax为储层岩心样品的核磁共振T2谱的特征参数;a、b、c为待求的参数,通过拟合得到。
一个实施例中,所述关系式建立模块进一步用于:
将多元统计关系式两边取对数,转化为线性表达式:
lg(k)=lg(a)+blg(T2,Smax)+clg(φ);
用最小二乘法拟合所述线性表达式,得到参数a、b、c。
一个实施例中,所述储层渗透率确定模块具体用于:
确定储层的核磁共振回波串;
对储层的核磁共振回波串进行反演,获得储层的核磁共振T2谱,计算储层的核磁共振T2谱的特征参数及储层的孔隙度;
根据多元统计关系式以及储层的T2谱的特征参数及储层的孔隙度,确定储层的渗透率。
本发明实施例中,选取储层岩心样品,确定储层岩心样品的渗透率及核磁共振回波串;对储层岩心样品的核磁共振回波串进行反演,获得储层岩心样品的核磁共振T2谱,计算储层岩心样品的核磁共振T2谱的特征参数及储层岩心样品的孔隙度;建立储层岩心样品的渗透率与储层岩心样品的核磁共振T2谱的特征参数和储层岩心样品的孔隙度的多元统计关系式;根据所述多元统计关系式,以及储层的T2谱的特征参数及储层的孔隙度,确定储层的渗透率;具有简单、精确、可靠性好等优点,具有明显的实际应用效果,为油藏分类、产能预测和储层建模提供了有力的技术支持。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。在附图中:
图1为本发明实施例中储层渗透率确定方法的示意图;
图2为本发明实施例中核磁共振T2谱及其反向累加归一化分布S的示意图;
图3为本发明实施例中计算核磁共振T2谱的特征参数的示意图;
图4为本发明实施例中本发明实施例方法确定的渗透率与岩心分析渗透率的交会图;
图5为本发明实施例中实际储层渗透率计算结果与岩心分析渗透率的具体实例对比图;
图6为本发明实施例中储层渗透率确定装置的示意图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚明白,下面结合附图对本发明实施例做进一步详细说明。在此,本发明的示意性实施例及其说明用于解释本发明,但并不作为对本发明的限定。
为了获得准确的储层渗透率,本发明实施例提供一种储层渗透率确定方法。图1为本发明实施例中储层渗透率确定方法的示意图,如图1所示,该方法包括:
步骤101、选取储层岩心样品,确定储层岩心样品的渗透率及核磁共振回波串;
步骤102、对储层岩心样品的核磁共振回波串进行反演,获得储层岩心样品的核磁共振T2谱,计算储层岩心样品的核磁共振T2谱的特征参数及储层岩心样品的孔隙度;
步骤103、建立储层岩心样品的渗透率与储层岩心样品的核磁共振T2谱的特征参数和储层岩心样品的孔隙度的多元统计关系式;
步骤104、根据所述多元统计关系式,以及储层的T2谱的特征参数及储层的孔隙度,确定储层的渗透率。
由图1所示流程可以得知,本发明实施例中是基于核磁共振T2谱进行储层渗透率计算,其中通过确定储层岩心样品的渗透率、计算储层岩心样品的核磁共振T2谱的特征参数及储层岩心样品的孔隙度,建立储层岩心样品的渗透率与储层岩心样品的核磁共振T2谱的特征参数和储层岩心样品的孔隙度的多元统计关系式,利用该多元统计关系式,以及储层的T2谱的特征参数及储层的孔隙度,确定出储层的渗透率,相对于现有技术而言,没有对束缚水饱和度的值敏感的缺点,无论是对水层还是对油气层都能得到准确、可靠的储层渗透率,并且计算过程简单,具有明显的实际应用效果,可以为油藏分类、产能预测和储层建模提供有力的技术支持。
具体实施时,先选取储层岩心样品,确定储层岩心样品的渗透率及核磁共振回波串。实施例中,可以选取研究区块岩心样品,获得测量得到的每块岩心样品的渗透率和核磁共振回波串,具体的,例如可以获得实验室岩心渗透率分析仪器测量得到的岩心样品的渗透率,获得低场核磁共振岩心分析仪测量得到的岩心样品的核磁共振回波串。实施例中,可以按照《岩心分析方法(SY/T5336-2006)》测量岩心样品的渗透率,按照《岩样核磁共振参数实验室测量规范(SY/T6490-2007)》测量岩心样品的核磁共振回波串。
实施时对储层岩心样品的核磁共振回波串进行反演,可以获得储层岩心样品的核磁共振T2谱,进而可以计算出储层岩心样品的核磁共振T2谱的特征参数及储层岩心样品的孔隙度。将核磁共振回波串反演可以得到每块岩心样品的核磁共振T2谱,根据核磁共振测井或者常规测井可以计算得到岩心样品的孔隙度,或者根据实验室岩心孔隙度分析仪器测量得到岩心样品的孔隙度。
具体的实例中,在计算储层岩心样品的核磁共振T2谱的特征参数时,可以包括:将储层岩心样品的核磁共振T2谱从长弛豫组分到短弛豫组分反向累加并归一化,获得储层岩心样品的核磁共振T2谱反向累积归一化分布S;将S和T2的乘积S×T2与S关系曲线中的最大值max(S×T2)作为储层岩心样品的核磁共振T2谱的特征参数T2,Smax。
图2为本发明实施例中核磁共振T2谱及其反向累加归一化分布S的示意图,如图2所示,实线为核磁共振T2谱,虚线为T2与S关系曲线。图3为本发明实施例中计算核磁共振T2谱的特征参数的示意图,如图3所示,其中空心三角形即为T2谱的特征参数。
实施例中,建立储层岩心样品的渗透率与储层岩心样品的核磁共振T2谱的特征参数和储层岩心样品的孔隙度的多元统计关系式,可以包括建立如下多元统计关系式:
其中,k为储层岩心样品的渗透率,mD;φ为储层岩心样品的孔隙度,小数;T2,Smax为储层岩心样品的核磁共振T2谱的特征参数;a、b、c为待求的参数,通过拟合得到。
实施例中,建立储层岩心样品的渗透率与储层岩心样品的核磁共振T2谱的特征参数和储层岩心样品的孔隙度的多元统计关系式,可以进一步包括:
将多元统计关系式两边取对数,转化为线性表达式:
lg(k)=lg(a)+blg(T2,Smax)+clg(φ); (2)
用最小二乘法拟合所述线性表达式,得到参数a、b、c。
表1给出了储层岩心样品的渗透率、储层岩心样品的核磁共振T2谱的特征参数和储层岩心样品的孔隙度的具体实例数据。
表1岩心样品的孔隙度、渗透率和T2谱的特征参数
以表1中所示实例数据为例,用最小二乘法拟合上述线性表达式,可以得到参数a=0.3793、b=0.88578、c=1.407,从而得到如下多元统计关系式:
具体的实施例中,根据所述多元统计关系式,以及储层的T2谱的特征参数及储层的孔隙度,确定储层的渗透率,可以包括:确定储层的核磁共振回波串;对储层的核磁共振回波串进行反演,获得储层的核磁共振T2谱,计算储层的核磁共振T2谱的特征参数及储层的孔隙度;根据多元统计关系式以及储层的T2谱的特征参数及储层的孔隙度,确定储层的渗透率。
仍以表1所示实例数据为例,实施时对于实际地层,获取待求储层的核磁共振回波串,对核磁共振回波串反演得到核磁共振T2谱,并且计算孔隙度,计算核磁共振T2谱的特征参数,利用公式(3)计算储层渗透率。图4为本发明实施例中本发明实施例方法确定的渗透率与岩心分析渗透率的交会图,如图4所示,本发明实施例方法确定的渗透率与岩心分析渗透率(实验测量得到的岩心样品渗透率)的相关性好,本发明实施例方法的计算精度高。
图5为本发明实施例中实际储层渗透率计算结果与岩心分析渗透率的具体实例对比图。如图5所示,本例中对某地区某井计算渗透率,图5中,第一道为自然伽马曲线(GR)和井径曲线(CAL),第二道为深度道,第三道为实测核磁共振回波串,第四道为核磁共振T2谱(T2-DIST)和T2几何均值(T2LM),第五道为核磁共振孔隙度(TCMR)和岩心分析孔隙度(PHIT-CORE),第六道为核磁共振T2谱的特征参数(T2SMAX),第七道为Coates模型计算的渗透率(PERM-COATES),第八道为SDR模型计算的渗透率(PERM-SDR),第九道为本发明实施例方法计算的渗透率(PERM-NEW),第十道为测井解释结论。从图5可以看出,本发明实施例方法计算的储层渗透率比Coates模型和SDR模型计算结果更精确。
基于同一发明构思,本发明实施例中还提供了一种储层渗透率确定装置,如下面的实施例所述。由于该装置解决问题的原理与储层渗透率确定方法相似,因此该装置的实施可以参见储层渗透率确定方法的实施,重复之处不再赘述。
图6为本发明实施例中储层渗透率确定装置的示意图,如图6所示,该装置可以包括:
样品参数确定模块601,用于选取储层岩心样品,确定储层岩心样品的渗透率及核磁共振回波串;
样品参数计算模块602,用于对储层岩心样品的核磁共振回波串进行反演,获得储层岩心样品的核磁共振T2谱,计算储层岩心样品的核磁共振T2谱的特征参数及储层岩心样品的孔隙度;
关系式建立模块603,用于建立储层岩心样品的渗透率与储层岩心样品的核磁共振T2谱的特征参数和储层岩心样品的孔隙度的多元统计关系式;
储层渗透率确定模块604,用于根据所述多元统计关系式,以及储层的T2谱的特征参数及储层的孔隙度,确定储层的渗透率。
具体实施时,样品参数计算模块602具体可以用于:
将储层岩心样品的核磁共振T2谱从长弛豫组分到短弛豫组分反向累加并归一化,获得储层岩心样品的核磁共振T2谱反向累积归一化分布S;
将S和T2的乘积S×T2与S关系曲线中的最大值max(S×T2)作为储层岩心样品的核磁共振T2谱的特征参数T2,Smax。
具体实施时,关系式建立模块603具体可以用于建立如下多元统计关系式:
其中,k为储层岩心样品的渗透率,mD;φ为储层岩心样品的孔隙度,小数;T2,Smax为储层岩心样品的核磁共振T2谱的特征参数;a、b、c为待求的参数,通过拟合得到。
具体实施时,关系式建立模块603可以进一步用于:
将多元统计关系式两边取对数,转化为线性表达式:
lg(k)=lg(a)+blg(T2,Smax)+clg(φ);
用最小二乘法拟合所述线性表达式,得到参数a、b、c。
具体实施时,储层渗透率确定模块604具体可以用于:
确定储层的核磁共振回波串;
对储层的核磁共振回波串进行反演,获得储层的核磁共振T2谱,计算储层的核磁共振T2谱的特征参数及储层的孔隙度;
根据多元统计关系式以及储层的T2谱的特征参数及储层的孔隙度,确定储层的渗透率。
综上所述,本发明实施例中,选取储层岩心样品,确定储层岩心样品的渗透率及核磁共振回波串;对储层岩心样品的核磁共振回波串进行反演,获得储层岩心样品的核磁共振T2谱,计算储层岩心样品的核磁共振T2谱的特征参数及储层岩心样品的孔隙度;建立储层岩心样品的渗透率与储层岩心样品的核磁共振T2谱的特征参数和储层岩心样品的孔隙度的多元统计关系式;根据所述多元统计关系式,以及储层的T2谱的特征参数及储层的孔隙度,确定储层的渗透率;具有简单、精确、可靠性好等优点,具有明显的实际应用效果,为油藏分类、产能预测和储层建模提供了有力的技术支持。
本领域内的技术人员应明白,本发明的实施例可提供为方法、系统、或计算机程序产品。因此,本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且,本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器,使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品,该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上,使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理,从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
以上所述的具体实施例,对本发明的目的、技术方案和有益效果进行了进一步详细说明,所应理解的是,以上所述仅为本发明的具体实施例而已,并不用于限定本发明的保护范围,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (6)
1.一种储层渗透率确定方法,其特征在于,包括:
选取储层岩心样品,确定储层岩心样品的渗透率及核磁共振回波串;
对储层岩心样品的核磁共振回波串进行反演,获得储层岩心样品的核磁共振T2谱,计算储层岩心样品的核磁共振T2谱的特征参数及储层岩心样品的孔隙度;
建立储层岩心样品的渗透率与储层岩心样品的核磁共振T2谱的特征参数和储层岩心样品的孔隙度的多元统计关系式;
根据所述多元统计关系式,以及储层的T2谱的特征参数及储层的孔隙度,确定储层的渗透率;
所述计算储层岩心样品的核磁共振T2谱的特征参数,包括:
将储层岩心样品的核磁共振T2谱从长弛豫组分到短弛豫组分反向累加并归一化,获得储层岩心样品的核磁共振T2谱反向累积归一化分布S;
将S和T2的乘积S×T2与S关系曲线中的最大值max(S×T2)作为储层岩心样品的核磁共振T2谱的特征参数T2,Smax;
所述建立储层岩心样品的渗透率与储层岩心样品的核磁共振T2谱的特征参数和储层岩心样品的孔隙度的多元统计关系式,包括建立如下多元统计关系式:
其中,k为储层岩心样品的渗透率,mD;φ为储层岩心样品的孔隙度,小数;T2,Smax为储层岩心样品的核磁共振T2谱的特征参数;a、b、c为待求的参数,通过拟合得到。
2.如权利要求1所述的方法,其特征在于,所述建立储层岩心样品的渗透率与储层岩心样品的核磁共振T2谱的特征参数和储层岩心样品的孔隙度的多元统计关系式,进一步包括:
将多元统计关系式两边取对数,转化为线性表达式:
lg(k)=lg(a)+blg(T2,Smax)+clg(φ);
用最小二乘法拟合所述线性表达式,得到参数a、b、c。
3.如权利要求1或2所述的方法,其特征在于,所述根据所述多元统计关系式,以及储层的T2谱的特征参数及储层的孔隙度,确定储层的渗透率,包括:
确定储层的核磁共振回波串;
对储层的核磁共振回波串进行反演,获得储层的核磁共振T2谱,计算储层的核磁共振T2谱的特征参数及储层的孔隙度;
根据多元统计关系式以及储层的T2谱的特征参数及储层的孔隙度,确定储层的渗透率。
4.一种储层渗透率确定装置,其特征在于,包括:
样品参数确定模块,用于选取储层岩心样品,确定储层岩心样品的渗透率及核磁共振回波串;
样品参数计算模块,用于对储层岩心样品的核磁共振回波串进行反演,获得储层岩心样品的核磁共振T2谱,计算储层岩心样品的核磁共振T2谱的特征参数及储层岩心样品的孔隙度;
关系式建立模块,用于建立储层岩心样品的渗透率与储层岩心样品的核磁共振T2谱的特征参数和储层岩心样品的孔隙度的多元统计关系式;
储层渗透率确定模块,用于根据所述多元统计关系式,以及储层的T2谱的特征参数及储层的孔隙度,确定储层的渗透率;
所述样品参数计算模块具体用于:
将储层岩心样品的核磁共振T2谱从长弛豫组分到短弛豫组分反向累加并归一化,获得储层岩心样品的核磁共振T2谱反向累积归一化分布S;
将S和T2的乘积S×T2与S关系曲线中的最大值max(S×T2)作为储层岩心样品的核磁共振T2谱的特征参数T2,Smax;
所述关系式建立模块具体用于建立如下多元统计关系式:
其中,k为储层岩心样品的渗透率,mD;φ为储层岩心样品的孔隙度,小数;T2,Smax为储层岩心样品的核磁共振T2谱的特征参数;a、b、c为待求的参数,通过拟合得到。
5.如权利要求4所述的装置,其特征在于,所述关系式建立模块进一步用于:
将多元统计关系式两边取对数,转化为线性表达式:
lg(k)=lg(a)+blg(T2,Smax)+clg(φ);
用最小二乘法拟合所述线性表达式,得到参数a、b、c。
6.如权利要求4或5所述的装置,其特征在于,所述储层渗透率确定模块具体用于:
确定储层的核磁共振回波串;
对储层的核磁共振回波串进行反演,获得储层的核磁共振T2谱,计算储层的核磁共振T2谱的特征参数及储层的孔隙度;
根据多元统计关系式以及储层的T2谱的特征参数及储层的孔隙度,确定储层的渗透率。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610330618.3A CN105891089B (zh) | 2016-05-18 | 2016-05-18 | 储层渗透率确定方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610330618.3A CN105891089B (zh) | 2016-05-18 | 2016-05-18 | 储层渗透率确定方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105891089A CN105891089A (zh) | 2016-08-24 |
CN105891089B true CN105891089B (zh) | 2018-10-12 |
Family
ID=56716344
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610330618.3A Active CN105891089B (zh) | 2016-05-18 | 2016-05-18 | 储层渗透率确定方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105891089B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2022053991A1 (en) * | 2020-09-10 | 2022-03-17 | Khalifa University of Science and Technology | Morphology decoder to predict heterogeneous rock permeability with machine learning guided 3d vision |
Families Citing this family (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107843937B (zh) * | 2016-09-20 | 2019-07-12 | 中国石油化工股份有限公司 | 低渗砂岩储层有效渗流能力的核磁共振表征方法 |
CN106383078B (zh) * | 2016-09-20 | 2019-05-07 | 中国石油天然气股份有限公司 | 岩石水驱效率的确定方法及装置 |
CN106644879B (zh) * | 2016-11-17 | 2019-09-10 | 中国石油天然气股份有限公司 | 一种确定岩心不同孔隙组分渗透率贡献值的方法及装置 |
CN107329182B (zh) * | 2017-07-26 | 2019-08-06 | 中国石油天然气股份有限公司 | 一种确定储层渗透率的方法及装置 |
CN110837132B (zh) * | 2018-08-15 | 2022-03-08 | 中国石油化工股份有限公司 | 碳酸盐岩测井渗透率预测方法 |
CN108918396A (zh) * | 2018-09-13 | 2018-11-30 | 中国地质大学(北京) | 一种用三轴应力测试致密储层岩石渗透率的方法 |
CN109612900B (zh) * | 2018-11-28 | 2020-05-12 | 中国石油大学(北京) | 一种储层岩石渗透率预测方法、装置及存储介质 |
CN113495046B (zh) * | 2020-04-03 | 2024-06-25 | 中国石油天然气股份有限公司 | 储层类型的确定方法、装置及存储介质 |
CN111707695B (zh) * | 2020-06-16 | 2021-09-10 | 中国石油大学(北京) | 基于核磁共振数据确定储层渗透率的方法、装置及设备 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102253069A (zh) * | 2011-04-20 | 2011-11-23 | 中国石油天然气股份有限公司 | 一种根据核磁共振t2谱确定渗透率的方法和装置 |
CN103335928A (zh) * | 2013-05-30 | 2013-10-02 | 中国石油天然气集团公司 | 一种测量孔隙岩石渗透率的方法和装置 |
CN103674803A (zh) * | 2013-11-18 | 2014-03-26 | 中国石油天然气集团公司 | 一种计算煤岩割理渗透率的方法及系统 |
CN104101562A (zh) * | 2013-04-15 | 2014-10-15 | 中国石油天然气集团公司 | 一种确定储层渗透率的方法 |
CN104374679A (zh) * | 2014-10-23 | 2015-02-25 | 中国石油天然气集团公司 | 一种确定储层渗透率贡献值的方法 |
CN104634804A (zh) * | 2013-11-08 | 2015-05-20 | 中国石油天然气股份有限公司 | 一种利用核磁共振t2谱确定储层相对渗透率的方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20060272812A1 (en) * | 2005-06-04 | 2006-12-07 | Gang Yu | Method for analyzing drill cuttings using nuclear magnetic resonance techniques |
-
2016
- 2016-05-18 CN CN201610330618.3A patent/CN105891089B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102253069A (zh) * | 2011-04-20 | 2011-11-23 | 中国石油天然气股份有限公司 | 一种根据核磁共振t2谱确定渗透率的方法和装置 |
CN104101562A (zh) * | 2013-04-15 | 2014-10-15 | 中国石油天然气集团公司 | 一种确定储层渗透率的方法 |
CN103335928A (zh) * | 2013-05-30 | 2013-10-02 | 中国石油天然气集团公司 | 一种测量孔隙岩石渗透率的方法和装置 |
CN104634804A (zh) * | 2013-11-08 | 2015-05-20 | 中国石油天然气股份有限公司 | 一种利用核磁共振t2谱确定储层相对渗透率的方法 |
CN103674803A (zh) * | 2013-11-18 | 2014-03-26 | 中国石油天然气集团公司 | 一种计算煤岩割理渗透率的方法及系统 |
CN104374679A (zh) * | 2014-10-23 | 2015-02-25 | 中国石油天然气集团公司 | 一种确定储层渗透率贡献值的方法 |
Non-Patent Citations (2)
Title |
---|
"核磁共振渗透率模型研究与应用";黄乔松等;《青岛大学学报》;20041231;第17卷(第4期);第1.3部分 * |
"用核磁共振测井技术评价储层渗透率特性";常文会等;《测井技术》;20051231;第29卷(第6期);全文 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2022053991A1 (en) * | 2020-09-10 | 2022-03-17 | Khalifa University of Science and Technology | Morphology decoder to predict heterogeneous rock permeability with machine learning guided 3d vision |
Also Published As
Publication number | Publication date |
---|---|
CN105891089A (zh) | 2016-08-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105891089B (zh) | 储层渗透率确定方法及装置 | |
CN103884633B (zh) | 一种确定岩石渗透率的方法及装置 | |
CN105464654B (zh) | 确定致密砂岩饱和度指数的方法及装置 | |
Mao et al. | Estimation of permeability by integrating nuclear magnetic resonance (NMR) logs with mercury injection capillary pressure (MICP) data in tight gas sands | |
Peters | Simple consistent models for water retention and hydraulic conductivity in the complete moisture range | |
CN101892837B (zh) | 地层因数确定方法及含油饱和度确定方法 | |
CN105223116B (zh) | 一种基于核磁共振谱系数法计算束缚水饱和度的方法 | |
CN104634804B (zh) | 一种利用核磁共振t2谱确定储层相对渗透率的方法 | |
CN104278989B (zh) | 一种获取低孔低渗储层饱和度指数的方法 | |
US10996364B2 (en) | System and method for determining water-filled porosity and water salinity in well cores and logs | |
CN105804734B (zh) | 一种利用核磁共振测井识别稠油储层的方法 | |
CN109100812B (zh) | 基于核磁共振的岩石孔隙分形维数评价方法及装置 | |
CN107688037B (zh) | 一种利用核磁测井t2分布确定井下岩石粒度曲线的方法 | |
CN109667573A (zh) | 三维页岩储层孔隙压力预测方法、装置和电子设备 | |
CN106355571A (zh) | 一种白云岩储层质量的确定方法及装置 | |
CN106644879B (zh) | 一种确定岩心不同孔隙组分渗透率贡献值的方法及装置 | |
CN107818379A (zh) | 吸附气量确定方法和装置 | |
CN109025961A (zh) | 页岩储层含气饱和度的计算方法、装置和电子设备 | |
CN104374679B (zh) | 一种确定储层渗透率贡献值的方法 | |
CN105844011B (zh) | 一种基于毛管模型的渗透率计算方法 | |
CN106570235B (zh) | 一种致密砂岩岩电参数获取方法及装置 | |
CN111208566B (zh) | 基于sca模型的孔缝参数反演方法、装置及存储介质 | |
CN107991710A (zh) | 一种储层孔径分布获取方法及装置 | |
Qobi et al. | Permeability determination from Stoneley waves in the Ara group carbonates, Oman | |
CN103675945B (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 |