CN105891089A - 储层渗透率确定方法及装置 - Google Patents

储层渗透率确定方法及装置 Download PDF

Info

Publication number
CN105891089A
CN105891089A CN201610330618.3A CN201610330618A CN105891089A CN 105891089 A CN105891089 A CN 105891089A CN 201610330618 A CN201610330618 A CN 201610330618A CN 105891089 A CN105891089 A CN 105891089A
Authority
CN
China
Prior art keywords
reservoir
core sample
magnetic resonance
spectrum
nuclear magnetic
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
Application number
CN201610330618.3A
Other languages
English (en)
Other versions
CN105891089B (zh
Inventor
谢然红
刘秘
段野
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China University of Petroleum Beijing
Original Assignee
China University of Petroleum Beijing
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by China University of Petroleum Beijing filed Critical China University of Petroleum Beijing
Priority to CN201610330618.3A priority Critical patent/CN105891089B/zh
Publication of CN105891089A publication Critical patent/CN105891089A/zh
Application granted granted Critical
Publication of CN105891089B publication Critical patent/CN105891089B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N15/00Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
    • G01N15/08Investigating permeability, pore-volume, or surface area of porous materials
    • EFIXED CONSTRUCTIONS
    • E21EARTH OR ROCK DRILLING; MINING
    • E21BEARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B47/00Survey of boreholes or wells
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N15/00Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
    • G01N15/08Investigating permeability, pore-volume, or surface area of porous materials
    • G01N15/088Investigating volume, surface area, size or distribution of pores; Porosimetry
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V5/00Prospecting or detecting by the use of ionising radiation, e.g. of natural or induced radioactivity
    • G01V5/04Prospecting 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 = aT 2 , S m a x b φ c ;
其中,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 = aT 2 , S m a x b φ c ;
其中,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 = aT 2 , S m a x b φ c ; - - - ( 1 )
其中,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,从而得到如下多元统计关系式:
k = 0.3793 × T 2 , S m a x 0.88578 φ 1.407 . - - - ( 3 )
具体的实施例中,根据所述多元统计关系式,以及储层的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 = aT 2 , S m a x b φ c ;
其中,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 (10)

1.一种储层渗透率确定方法,其特征在于,包括:
选取储层岩心样品,确定储层岩心样品的渗透率及核磁共振回波串;
对储层岩心样品的核磁共振回波串进行反演,获得储层岩心样品的核磁共振T2谱,计算储层岩心样品的核磁共振T2谱的特征参数及储层岩心样品的孔隙度;
建立储层岩心样品的渗透率与储层岩心样品的核磁共振T2谱的特征参数和储层岩心样品的孔隙度的多元统计关系式;
根据所述多元统计关系式,以及储层的T2谱的特征参数及储层的孔隙度,确定储层的渗透率。
2.如权利要求1所述的方法,其特征在于,所述计算储层岩心样品的核磁共振T2谱的特征参数,包括:
将储层岩心样品的核磁共振T2谱从长弛豫组分到短弛豫组分反向累加并归一化,获得储层岩心样品的核磁共振T2谱反向累积归一化分布S;
将S和T2的乘积S×T2与S关系曲线中的最大值max(S×T2)作为储层岩心样品的核磁共振T2谱的特征参数T2,Smax
3.如权利要求1所述的方法,其特征在于,所述建立储层岩心样品的渗透率与储层岩心样品的核磁共振T2谱的特征参数和储层岩心样品的孔隙度的多元统计关系式,包括建立如下多元统计关系式:
k = aT 2 , S m a x b φ c ;
其中,k为储层岩心样品的渗透率,mD;φ为储层岩心样品的孔隙度,小数;T2,Smax为储层岩心样品的核磁共振T2谱的特征参数;a、b、c为待求的参数,通过拟合得到。
4.如权利要求3所述的方法,其特征在于,所述建立储层岩心样品的渗透率与储层岩心样品的核磁共振T2谱的特征参数和储层岩心样品的孔隙度的多元统计关系式,进一步包括:
将多元统计关系式两边取对数,转化为线性表达式:
lg(k)=lg(a)+blg(T2,Smax)+clg(φ);
用最小二乘法拟合所述线性表达式,得到参数a、b、c。
5.如权利要求3或4所述的方法,其特征在于,所述根据所述多元统计关系式,以及储层的T2谱的特征参数及储层的孔隙度,确定储层的渗透率,包括:
确定储层的核磁共振回波串;
对储层的核磁共振回波串进行反演,获得储层的核磁共振T2谱,计算储层的核磁共振T2谱的特征参数及储层的孔隙度;
根据多元统计关系式以及储层的T2谱的特征参数及储层的孔隙度,确定储层的渗透率。
6.一种储层渗透率确定装置,其特征在于,包括:
样品参数确定模块,用于选取储层岩心样品,确定储层岩心样品的渗透率及核磁共振回波串;
样品参数计算模块,用于对储层岩心样品的核磁共振回波串进行反演,获得储层岩心样品的核磁共振T2谱,计算储层岩心样品的核磁共振T2谱的特征参数及储层岩心样品的孔隙度;
关系式建立模块,用于建立储层岩心样品的渗透率与储层岩心样品的核磁共振T2谱的特征参数和储层岩心样品的孔隙度的多元统计关系式;
储层渗透率确定模块,用于根据所述多元统计关系式,以及储层的T2谱的特征参数及储层的孔隙度,确定储层的渗透率。
7.如权利要求6所述的装置,其特征在于,所述样品参数计算模块具体用于:
将储层岩心样品的核磁共振T2谱从长弛豫组分到短弛豫组分反向累加并归一化,获得储层岩心样品的核磁共振T2谱反向累积归一化分布S;
将S和T2的乘积S×T2与S关系曲线中的最大值max(S×T2)作为储层岩心样品的核磁共振T2谱的特征参数T2,Smax
8.如权利要求6所述的装置,其特征在于,所述关系式建立模块具体用于建立如下多元统计关系式:
k = aT 2 , S m a x b φ c ;
其中,k为储层岩心样品的渗透率,mD;φ为储层岩心样品的孔隙度,小数;T2,Smax为储层岩心样品的核磁共振T2谱的特征参数;a、b、c为待求的参数,通过拟合得到。
9.如权利要求8所述的装置,其特征在于,所述关系式建立模块进一步用于:
将多元统计关系式两边取对数,转化为线性表达式:
lg(k)=lg(a)+blg(T2,Smax)+clg(φ);
用最小二乘法拟合所述线性表达式,得到参数a、b、c。
10.如权利要求8或9所述的装置,其特征在于,所述储层渗透率确定模块具体用于:
确定储层的核磁共振回波串;
对储层的核磁共振回波串进行反演,获得储层的核磁共振T2谱,计算储层的核磁共振T2谱的特征参数及储层的孔隙度;
根据多元统计关系式以及储层的T2谱的特征参数及储层的孔隙度,确定储层的渗透率。
CN201610330618.3A 2016-05-18 2016-05-18 储层渗透率确定方法及装置 Active CN105891089B (zh)

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 true CN105891089A (zh) 2016-08-24
CN105891089B 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 (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106383078A (zh) * 2016-09-20 2017-02-08 中国石油天然气股份有限公司 岩石水驱效率的确定方法及装置
CN106644879A (zh) * 2016-11-17 2017-05-10 中国石油天然气股份有限公司 一种确定岩心不同孔隙组分渗透率贡献值的方法及装置
CN107329182A (zh) * 2017-07-26 2017-11-07 中国石油天然气股份有限公司 一种确定储层渗透率的方法及装置
CN107843937A (zh) * 2016-09-20 2018-03-27 中国石油化工股份有限公司 低渗砂岩储层有效渗流能力的核磁共振表征方法
CN108918396A (zh) * 2018-09-13 2018-11-30 中国地质大学(北京) 一种用三轴应力测试致密储层岩石渗透率的方法
CN109612900A (zh) * 2018-11-28 2019-04-12 中国石油大学(北京) 一种储层岩石渗透率预测方法、装置及存储介质
CN110837132A (zh) * 2018-08-15 2020-02-25 中国石油化工股份有限公司 碳酸盐岩测井渗透率预测方法
CN111707695A (zh) * 2020-06-16 2020-09-25 中国石油大学(北京) 基于核磁共振数据确定储层渗透率的方法、装置及设备
CN113495046A (zh) * 2020-04-03 2021-10-12 中国石油天然气股份有限公司 储层类型的确定方法、装置及存储介质

Families Citing this family (1)

* Cited by examiner, † Cited by third party
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

Citations (7)

* Cited by examiner, † Cited by third party
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
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谱确定储层相对渗透率的方法

Patent Citations (7)

* Cited by examiner, † Cited by third party
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
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)

* Cited by examiner, † Cited by third party
Title
常文会等: ""用核磁共振测井技术评价储层渗透率特性"", 《测井技术》 *
黄乔松等: ""核磁共振渗透率模型研究与应用"", 《青岛大学学报》 *

Cited By (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106383078B (zh) * 2016-09-20 2019-05-07 中国石油天然气股份有限公司 岩石水驱效率的确定方法及装置
CN107843937A (zh) * 2016-09-20 2018-03-27 中国石油化工股份有限公司 低渗砂岩储层有效渗流能力的核磁共振表征方法
CN106383078A (zh) * 2016-09-20 2017-02-08 中国石油天然气股份有限公司 岩石水驱效率的确定方法及装置
CN106644879A (zh) * 2016-11-17 2017-05-10 中国石油天然气股份有限公司 一种确定岩心不同孔隙组分渗透率贡献值的方法及装置
CN106644879B (zh) * 2016-11-17 2019-09-10 中国石油天然气股份有限公司 一种确定岩心不同孔隙组分渗透率贡献值的方法及装置
CN107329182A (zh) * 2017-07-26 2017-11-07 中国石油天然气股份有限公司 一种确定储层渗透率的方法及装置
CN110837132A (zh) * 2018-08-15 2020-02-25 中国石油化工股份有限公司 碳酸盐岩测井渗透率预测方法
CN110837132B (zh) * 2018-08-15 2022-03-08 中国石油化工股份有限公司 碳酸盐岩测井渗透率预测方法
CN108918396A (zh) * 2018-09-13 2018-11-30 中国地质大学(北京) 一种用三轴应力测试致密储层岩石渗透率的方法
CN109612900A (zh) * 2018-11-28 2019-04-12 中国石油大学(北京) 一种储层岩石渗透率预测方法、装置及存储介质
CN109612900B (zh) * 2018-11-28 2020-05-12 中国石油大学(北京) 一种储层岩石渗透率预测方法、装置及存储介质
CN113495046A (zh) * 2020-04-03 2021-10-12 中国石油天然气股份有限公司 储层类型的确定方法、装置及存储介质
CN111707695A (zh) * 2020-06-16 2020-09-25 中国石油大学(北京) 基于核磁共振数据确定储层渗透率的方法、装置及设备
CN111707695B (zh) * 2020-06-16 2021-09-10 中国石油大学(北京) 基于核磁共振数据确定储层渗透率的方法、装置及设备

Also Published As

Publication number Publication date
CN105891089B (zh) 2018-10-12

Similar Documents

Publication Publication Date Title
CN105891089B (zh) 储层渗透率确定方法及装置
CN103884633B (zh) 一种确定岩石渗透率的方法及装置
EP3333360B1 (en) Method of acquiring rock component content of stratum
Timur Effective porosity and permeability of sandstones investigated through nuclear magnetic resonance principles
CN104278989B (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
CN103225500B (zh) 一种应用三参数自洽迭代算法的新型水淹层测井评价方法
CN104634804B (zh) 一种利用核磁共振t2谱确定储层相对渗透率的方法
CN103674811B (zh) 一种核磁共振孔隙度测量的校正方法、装置及系统
CN104215652B (zh) 确定油气饱和度的方法和装置
CN104237957B (zh) 一种核磁共振测井t2截止值的确定方法
CN107525553B (zh) 一种确定多相流体组分流量的方法及装置
CN109100812B (zh) 基于核磁共振的岩石孔隙分形维数评价方法及装置
CN103760082A (zh) 页岩层系致密储层原油有效可动空间确定方法和装置
CN106525888A (zh) 一种测试致密油藏润湿性的方法及装置
CN106355571A (zh) 一种白云岩储层质量的确定方法及装置
CN107679358A (zh) 一种确定储层渗透率的方法及装置
CN103513270A (zh) 一种基于岩石声学特性的气层识别评价方法及装置
CN105158796A (zh) 确定toc含量的方法和装置
CN106644879A (zh) 一种确定岩心不同孔隙组分渗透率贡献值的方法及装置
Saboorian-Jooybari et al. Fracture identification and comprehensive evaluation of the parameters by dual laterolog data
NO20120654A1 (no) Fremgangsmåte og anordning for inkorporering av intern gradient og begrenset diffusjon i NMR-inversjon ved brønnlogging
CN104374679B (zh) 一种确定储层渗透率贡献值的方法
CN105844011B (zh) 一种基于毛管模型的渗透率计算方法
CN110320575A (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