CN110306977A - 核磁共振t2谱含烃校正方法及系统 - Google Patents
核磁共振t2谱含烃校正方法及系统 Download PDFInfo
- Publication number
- CN110306977A CN110306977A CN201810260213.6A CN201810260213A CN110306977A CN 110306977 A CN110306977 A CN 110306977A CN 201810260213 A CN201810260213 A CN 201810260213A CN 110306977 A CN110306977 A CN 110306977A
- Authority
- CN
- China
- Prior art keywords
- magnetic resonance
- nuclear magnetic
- spectrum
- wave spectrum
- water wave
- 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
- 238000005481 NMR spectroscopy Methods 0.000 title claims abstract description 98
- 238000001228 spectrum Methods 0.000 claims abstract description 176
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 118
- 238000005457 optimization Methods 0.000 claims abstract description 35
- 239000011148 porous material Substances 0.000 claims abstract description 26
- 230000035699 permeability Effects 0.000 claims abstract description 22
- 238000002474 experimental method Methods 0.000 claims abstract description 14
- 238000012937 correction Methods 0.000 claims description 12
- 238000000034 method Methods 0.000 abstract description 14
- 230000000694 effects Effects 0.000 abstract description 7
- 230000008569 process Effects 0.000 abstract description 5
- 238000011158 quantitative evaluation Methods 0.000 abstract description 3
- 239000003921 oil Substances 0.000 description 15
- 239000011435 rock Substances 0.000 description 13
- 230000008901 benefit Effects 0.000 description 4
- 238000009792 diffusion process Methods 0.000 description 4
- 239000004215 Carbon black (E152) Substances 0.000 description 3
- 239000012530 fluid Substances 0.000 description 3
- 229930195733 hydrocarbon Natural products 0.000 description 3
- 150000002430 hydrocarbons Chemical class 0.000 description 3
- 230000005311 nuclear magnetism Effects 0.000 description 3
- 239000010779 crude oil Substances 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 239000003208 petroleum Substances 0.000 description 2
- 230000003595 spectral effect Effects 0.000 description 2
- 238000009736 wetting Methods 0.000 description 2
- 241000772415 Neovison vison Species 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000011157 data evaluation Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 239000004744 fabric Substances 0.000 description 1
- 238000001764 infiltration Methods 0.000 description 1
- 230000008595 infiltration 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
- 238000011160 research Methods 0.000 description 1
- 229920006395 saturated elastomer Polymers 0.000 description 1
- 239000004575 stone Substances 0.000 description 1
- 238000012916 structural analysis Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
Classifications
-
- 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
- E21B49/00—Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
- G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
- G01R33/46—NMR spectroscopy
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V3/00—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
- G01V3/14—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation operating with electron or nuclear magnetic resonance
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/02—Agriculture; Fishing; Forestry; Mining
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/30—Assessment of water resources
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Mining & Mineral Resources (AREA)
- Geology (AREA)
- General Physics & Mathematics (AREA)
- High Energy & Nuclear Physics (AREA)
- Environmental & Geological Engineering (AREA)
- Business, Economics & Management (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Condensed Matter Physics & Semiconductors (AREA)
- General Health & Medical Sciences (AREA)
- Geochemistry & Mineralogy (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Geophysics (AREA)
- Agronomy & Crop Science (AREA)
- Animal Husbandry (AREA)
- Marine Sciences & Fisheries (AREA)
- Remote Sensing (AREA)
- Health & Medical Sciences (AREA)
- Economics (AREA)
- Fluid Mechanics (AREA)
- Human Resources & Organizations (AREA)
- Marketing (AREA)
- Primary Health Care (AREA)
- Strategic Management (AREA)
- Tourism & Hospitality (AREA)
- General Business, Economics & Management (AREA)
- Theoretical Computer Science (AREA)
- Geophysics And Detection Of Objects (AREA)
- Other Investigation Or Analysis Of Materials By Electrical Means (AREA)
Abstract
公开了一种核磁共振T2谱含烃校正方法及系统。该方法可以包括:采集目标储层的岩心样品进行核磁共振实验,获取核磁共振T2谱以及相应深度的孔隙度φ、渗透率K、几何平均值T2g;根据孔隙度φ、渗透率K、几何平均值T2g计算特征参数;将核磁共振T2谱划分为束缚水波谱和可动水波谱,通过特征参数重构可动水波谱,获得优化可动水波谱;将优化可动水波谱与束缚水波谱进行组合,即获得饱含水状态下的核磁共振T2谱。本发明将含烃的核磁共振T2谱校正到100%饱含水的核磁共振T2谱,以实现利用其连续定量评价储层孔隙结构的目的,流程简单,易于操作,实际应用效果好。
Description
技术领域
本发明涉及石油地球物理勘探领域,更具体地,涉及一种核磁共振T2谱含烃校正方法及系统。
背景技术
核磁共振测井自从二十世纪九十年代被引入测井行业以来,由于其在获取储层束缚水饱和度、总孔隙度、有效孔隙度和渗透率等参数以及定量评价储集层孔隙结构方面的作用,受到岩石物理学家及测井分析家的普遍重视。尤其是在储集层孔隙结构定量评价方面,具有其他测井方法所无法比拟的独特优势。
据核磁共振的基本原理,核磁共振T2弛豫时间主要由三部分来组成,包括体积弛豫,表面弛豫和扩散弛豫,即:
式中,T2是孔隙流体横向弛豫时间,T2B是流体体弛豫,T2S是流体表明弛豫,T2D是流体扩散弛豫。
利用核磁共振测井评价孔隙结构的理论基础认为,对于水润湿和岩石而言,对于100%饱含水的岩石,在磁场均匀,扩散系数不大且岩石孔隙具有规则几何形状的情况下,其体积弛豫和扩散弛豫均可以忽略。此时,T2弛豫时间主要是表面弛豫,此种情况下,T2弛豫时间表示为:
式中,ρ2为岩石的表面弛豫率,S为岩石孔隙表面积,V为岩石孔隙体积,rC为岩石孔喉半径。由此可以看出,对于100%饱含水的水润湿相岩石而言,核磁共振T2弛豫时间与岩石的孔隙半径成正比,即孔隙结构较小的岩石,对应的T2弛豫时间也很短,在核磁共振T2谱上谱峰的位置相对较左。反之,孔隙半径较大的岩石,对应的T2弛豫时间较长,在核磁共振T2谱上谱峰的位置相对较右。因此,根据核磁共振T2谱的形态和相对位置来评价储集层的孔隙半径大小及分布。
然而,当储层岩石的孔隙空间含有非润湿相的烃之后,由于烃的体积弛豫不可忽略,导致核磁共振T2谱的形态和位置会发生变化。此时,核磁共振T2谱不能很好地反映岩石的孔隙分布。前人的研究成果主要集中在基于100%含水的核磁共振T2谱表征孔隙结构上(刘堂宴,马在田,付容珊.核磁共振谱的岩石孔喉结构分析[J].石油地球物理勘探,2003,38(3):328-333。邵维志,丁娱娇,刘亚。核磁共振测井在储层孔隙结构评价中的应用[J].测井技术,2009,33(1):52-56。),而对于非润湿相烃对核磁共振T2谱的影响,鲜有人提出。肖亮(肖亮.利用核磁共振测井资料评价储集层孔隙结构的讨论.新疆石油地质,2008,29(2):260-263)分析了孔隙含烃对核磁共振T2谱的影响,利用“三孔隙组分百分比法”来评价孔隙结构,但是并未提供具体的含烃校正方法。因此,有必要开发一种核磁共振T2谱含烃校正方法及系统。
公开于本发明背景技术部分的信息仅仅旨在加深对本发明的一般背景技术的理解,而不应当被视为承认或以任何形式暗示该信息构成已为本领域技术人员所公知的现有技术。
发明内容
本发明提出了一种核磁共振T2谱含烃校正方法及系统,其能够将含烃的核磁共振T2谱校正到100%饱含水的核磁共振T2谱,以实现利用其连续定量评价储层孔隙结构的目的,流程简单,易于操作,实际应用效果好。
根据本发明的一方面,提出了一种核磁共振T2谱含烃校正方法。所述方法可以包括:采集目标储层的岩心样品进行核磁共振实验,获取核磁共振T2谱以及相应深度的孔隙度φ、渗透率K、几何平均值T2g;根据所述孔隙度φ、渗透率K、几何平均值T2g计算特征参数;将所述核磁共振T2谱划分为束缚水波谱和可动水波谱,通过所述特征参数重构所述可动水波谱,获得优化可动水波谱;将所述优化可动水波谱与所述束缚水波谱进行组合,即获得饱含水状态下的核磁共振T2谱。
优选地,所述特征参数包括所述核磁共振T2谱的截止值T2C、峰值T2f、右边界T2R。
优选地,根据孔隙度φ、渗透率K、几何平均值T2g通过公式(1)、公式(2)、公式(3)计算峰值T2f、右边界T2R:
T2f=a1*(T2g)2+a2*(T2g)+a3 (1)
T2R=b1*T2s+b2 (2)
其中,T2s为孔隙结构特征参数,a1、a2、a3、b1、b2为公式拟合系数。
优选地,所述通过所述特征参数重构所述可动水波谱,获得优化可动水波谱包括:根据所述峰值T2f与截止值T2C的相对位置确定优化可动水波谱的波峰位置和/或波谱形态。
优选地,根据所述峰值T2f与截止值T2C的相对位置确定优化可动水波谱的波峰位置和/或波谱形态包括:当峰值T2f位于截止值T2C左侧时,所述可动水波谱的形态为单斜分布;当峰值T2f在截止值T2C右侧时,所述可动水波谱的波峰位置为峰值T2f,形态为正态分布。
优选地,还包括:采集储层油样进行原油性质实验,获得原油粘度,当原油粘度为轻质油时,对核磁共振T2谱进行校正。
根据本发明的另一方面,提出了一种核磁共振T2谱含烃校正系统,可以包括:参数获取单元,用于采集目标储层的岩心样品进行核磁共振实验,获取核磁共振T2谱以及相应深度的孔隙度φ、渗透率K、几何平均值T2g;计算单元,用于根据所述孔隙度φ、渗透率K、几何平均值T2g计算特征参数;优化可动水波谱单元,用于将所述核磁共振T2谱划分为束缚水波谱和可动水波谱,通过所述特征参数重构所述可动水波谱,获得优化可动水波谱;组合单元,用于将所述优化可动水波谱与所述束缚水波谱进行组合,即获得饱含水状态下的核磁共振T2谱。
优选地,所述特征参数包括所述核磁共振T2谱的截止值T2C、峰值T2f、右边界T2R。
优选地,根据孔隙度φ、渗透率K、几何平均值T2g通过公式(1)、公式(2)、公式(3)计算峰值T2f、右边界T2R:
T2f=a1*(T2g)2+a2*(T2g)+a3 (1)
T2R=b1*T2s+b2 (2)
其中,T2s为孔隙结构特征参数,a1、a2、a3、b1、b2为公式拟合系数。
优选地,所述通过所述特征参数重构所述可动水波谱,获得优化可动水波谱包括:根据所述峰值T2f与截止值T2C的相对位置确定优化可动水波谱的波峰位置和/或波谱形态:所述优化可动水波谱的波峰位置和波谱形态由峰值T2f与截止值T2C的相对位置确定包括:当峰值T2f位于截止值T2C左侧时,所述可动水波谱的形态为单斜分布;当峰值T2f在截止值T2C右侧时,所述可动水波谱的波峰位置为峰值T2f,形态为正态分布。
本发明具有其它的特性和优点,这些特性和优点从并入本文中的附图和随后的具体实施方式中将是显而易见的,或者将在并入本文中的附图和随后的具体实施方式中进行详细陈述,这些附图和具体实施方式共同用于解释本发明的特定原理。
附图说明
通过结合附图对本发明示例性实施例进行更详细的描述,本发明的上述以及其它目的、特征和优势将变得更加明显,其中,在本发明示例性实施例中,相同的参考标号通常代表相同部件。
图1示出了根据本发明的核磁共振T2谱含烃校正方法的步骤的流程图。
图2示出了根据本发明的一个实施例的峰值与几何平均值的关系图。
图3示出了根据本发明的一个实施例的右边界与孔隙结构特征参数的关系图。
图4示出了根据本发明的一个实施例的饱含水状态下的核磁共振T2谱的示意图。
具体实施方式
下面将参照附图更详细地描述本发明。虽然附图中显示了本发明的优选实施例,然而应该理解,可以以各种形式实现本发明而不应被这里阐述的实施例所限制。相反,提供这些实施例是为了使本发明更加透彻和完整,并且能够将本发明的范围完整地传达给本领域的技术人员。
图1示出了根据本发明的核磁共振T2谱含烃校正方法的步骤的流程图。
在该实施例中,根据本发明的核磁共振T2谱含烃校正方法可以包括:步骤101,采集目标储层的岩心样品进行核磁共振实验,获取核磁共振T2谱以及相应深度的孔隙度φ、渗透率K、几何平均值T2g;步骤102,根据孔隙度φ、渗透率K、几何平均值T2g计算特征参数;步骤103,将核磁共振T2谱划分为束缚水波谱和可动水波谱,通过特征参数重构可动水波谱,获得优化可动水波谱;步骤104,将优化可动水波谱与束缚水波谱进行组合,即获得饱含水状态下的核磁共振T2谱。
在一个示例中,特征参数包括核磁共振T2谱的截止值T2C、峰值T2f、右边界T2R。
在一个示例中,根据孔隙度φ、渗透率K、几何平均值T2g通过公式(1)、公式(2)、公式(3)计算峰值T2f、右边界T2R:
T2f=a1*(T2g)2+a2*(T2g)+a3 (1)
T2R=b1*T2s+b2 (2)
其中,T2s为孔隙结构特征参数,a1、a2、a3、b1、b2为公式拟合系数,根据核磁共振实验确定。
在一个示例中,通过特征参数重构可动水波谱,获得优化可动水波谱包括:根据峰值T2f与截止值T2C的相对位置确定优化可动水波谱的波峰位置和/或波谱形态。
在一个示例中,根据峰值T2f与截止值T2C的相对位置确定优化可动水波谱的波峰位置和/或波谱形态包括:当峰值T2f位于截止值T2C左侧时,可动水波谱的形态为单斜分布;当峰值T2f在截止值T2C右侧时,可动水波谱的波峰位置为峰值T2f,形态为正态分布。
在一个示例中,还包括:采集储层油样进行原油性质实验,获得原油粘度,当原油粘度为轻质油时,对核磁共振T2谱进行校正。
具体地,采集目标储层的岩心样品进行核磁共振实验,采集储层油样进行原油性质实验,获得原油粘度,当原油粘度为中等粘度时,即原油密度在0.852-0.930g/cm3之间,孔隙含烃对核磁共振T2谱形态不会造成影响,无需对核磁共振T2谱进行校正;当原油粘度为重质油时,即原油密度在0.931-0.998g/cm3之间,饱和油岩样的核磁共振T2谱形态已经完全失真,已经不能用于反映储层的孔隙结构;对于原油粘度为轻质油时,即原油密度小于0.852g/cm3时,获取核磁共振T2谱以及相应深度的孔隙度φ、渗透率K、几何平均值T2g,采用直方图方法,选取频率最高值作为所有岩心统一的截止值T2C,进而根据孔隙度φ、渗透率K、几何平均值T2g通过公式(1)、公式(2)、公式(3)计算峰值T2f、右边界T2R;根据截止值T2C,将核磁共振T2谱划分为束缚水波谱和可动水波谱,对于束缚水波谱,孔隙含烃对其形态不会造成影响,无需做核磁共振T2谱含烃校正,通过特征参数重构可动水波谱,优化可动水波谱的波峰位置和波谱形态由峰值T2f与截止值T2C的相对位置确定:当峰值T2f位于截止值T2C左侧时,可动水波谱的形态为单斜分布;当峰值T2f在截止值T2C右侧时,可动水波谱的波峰位置为峰值T2f,形态为正态分布;将优化可动水波谱替换可动水波谱,与束缚水波谱进行组合,即获得饱含水状态下的核磁共振T2谱。
当可动水波谱的形态为正态分布时,通过公式(4)、(5)进行重建:
其中,μ为峰值T2f,T2max为T2谱最大值,S(T2)为T2谱幅度,T2为核磁横向弛豫。
本发明将含烃的核磁共振T2谱校正到100%饱含水的核磁共振T2谱,以实现利用其连续定量评价储层孔隙结构的目的,流程简单,易于操作,实际应用效果好。
应用示例
为便于理解本发明实施例的方案及其效果,以下给出一个具体应用示例。本领域技术人员应理解,该示例仅为了便于理解本发明,其任何具体细节并非意在以任何方式限制本发明。
采集目标储层的岩心样品进行核磁共振实验,采集储层油样进行原油性质实验,获得原油粘度为15cp,当储层含烃时,核磁共振T2谱需要做含烃校正。获取核磁共振T2谱以及相应深度的截止值T2C为36ms,利用核磁实验资料确定T2峰值T2f、T2右边界T2R,并建立与峰值T2f、右边界T2R的相关关系,如图2和图3所示,拟合公式如下:
T2f=0.03027(T2g)2+3.7821(T2g)-8.8057 (6)
T2R=59.599*(T2s)+108.8 (7);
通过截止值T2C将核磁共振T2谱划分为束缚水波谱和可动水波谱,对可动水波谱进行重构,根据峰值T2f与截止值T2C确定优化可动水波谱的波谱形态为正态分布,通过公式(4)、(5)进行重建。将优化可动水波谱替换可动水波谱,与束缚水波谱进行组合,即获得饱含水状态下的核磁共振T2谱,如图4所示,第4列为原始的核磁共振T2谱,根据截止值进行划分,左侧为束缚水波谱,右侧为可动水波谱,将优化可动水波谱替换可动水波谱,得到第5列,即为含烃校正后的核磁共振T2谱,即饱含水状态下的核磁共振T2谱,第6列为校正前后对比。
综上所述,本发明将含烃的核磁共振T2谱校正到100%饱含水的核磁共振T2谱,以实现利用其连续定量评价储层孔隙结构的目的,流程简单,易于操作,实际应用效果好。
根据本发明的核磁共振T2谱含烃校正系统,可以包括:参数获取单元,用于采集目标储层的岩心样品进行核磁共振实验,获取核磁共振T2谱以及相应深度的孔隙度φ、渗透率K、几何平均值T2g;计算单元,用于根据孔隙度φ、渗透率K、几何平均值T2g计算特征参数;优化可动水波谱单元,用于将核磁共振T2谱划分为束缚水波谱和可动水波谱,通过特征参数重构可动水波谱,获得优化可动水波谱;组合单元,用于将优化可动水波谱与束缚水波谱进行组合,即获得饱含水状态下的核磁共振T2谱。
优选地,特征参数包括核磁共振T2谱的截止值T2C、峰值T2f、右边界T2R。
优选地,根据孔隙度φ、渗透率K、几何平均值T2g通过公式(1)、公式(2)、公式(3)计算峰值T2f、右边界T2R:
T2f=a1*(T2g)2+a2*(T2g)+a3 (1)
T2R=b1*T2s+b2 (2)
其中,T2s为孔隙结构特征参数,a1、a2、a3、b1、b2为公式拟合系数。
优选地,通过特征参数重构可动水波谱,获得优化可动水波谱包括:根据峰值T2f与截止值T2C的相对位置确定优化可动水波谱的波峰位置和/或波谱形态:优化可动水波谱的波峰位置和波谱形态由峰值T2f与截止值T2C的相对位置确定包括:当峰值T2f位于截止值T2C左侧时,可动水波谱的形态为单斜分布;当峰值T2f在截止值T2C右侧时,可动水波谱的波峰位置为峰值T2f,形态为正态分布。
本系统将含烃的核磁共振T2谱校正到100%饱含水的核磁共振T2谱,以实现利用其连续定量评价储层孔隙结构的目的,流程简单,易于操作,实际应用效果好。
本领域技术人员应理解,上面对本发明的实施例的描述的目的仅为了示例性地说明本发明的实施例的有益效果,并不意在将本发明的实施例限制于所给出的任何示例。
以上已经描述了本发明的各实施例,上述说明是示例性的,并非穷尽性的,并且也不限于所披露的各实施例。在不偏离所说明的各实施例的范围和精神的情况下,对于本技术领域的普通技术人员来说许多修改和变更都是显而易见的。
Claims (10)
1.一种核磁共振T2谱含烃校正方法,包括:
采集目标储层的岩心样品进行核磁共振实验,获取核磁共振T2谱以及相应深度的孔隙度φ、渗透率K、几何平均值T2g;
根据所述孔隙度φ、渗透率K、几何平均值T2g计算特征参数;
将所述核磁共振T2谱划分为束缚水波谱和可动水波谱,通过所述特征参数重构所述可动水波谱,获得优化可动水波谱;
将所述优化可动水波谱与所述束缚水波谱进行组合,即获得饱含水状态下的核磁共振T2谱。
2.根据权利要求1所述的核磁共振T2谱含烃校正方法,其中,所述特征参数包括所述核磁共振T2谱的截止值T2C、峰值T2f、右边界T2R。
3.根据权利要求2所述的核磁共振T2谱含烃校正方法,其中,根据孔隙度φ、渗透率K、几何平均值T2g通过公式(1)、公式(2)、公式(3)计算峰值T2f、右边界T2R:
T2f=a1*(T2g)2+a2*(T2g)+a3 (1)
T2R=b1*T2s+b2 (2)
其中,T2s为孔隙结构特征参数,a1、a2、a3、b1、b2为公式拟合系数。
4.根据权利要求2所述的核磁共振T2谱含烃校正方法,其中,所述通过所述特征参数重构所述可动水波谱,获得优化可动水波谱包括:
根据所述峰值T2f与截止值T2C的相对位置确定优化可动水波谱的波峰位置和/或波谱形态。
5.根据权利要求4所述的核磁共振T2谱含烃校正方法,其中,根据所述峰值T2f与截止值T2C的相对位置确定优化可动水波谱的波峰位置和/或波谱形态包括:
当峰值T2f位于截止值T2C左侧时,所述可动水波谱的形态为单斜分布;当峰值T2f在截止值T2C右侧时,所述可动水波谱的波峰位置为峰值T2f,形态为正态分布。
6.根据权利要求1所述的核磁共振T2谱含烃校正方法,其中,还包括:
采集储层油样进行原油性质实验,获得原油粘度,当原油粘度为轻质油时,对核磁共振T2谱进行校正。
7.一种核磁共振T2谱含烃校正系统,包括:
参数获取单元,用于采集目标储层的岩心样品进行核磁共振实验,获取核磁共振T2谱以及相应深度的孔隙度φ、渗透率K、几何平均值T2g;
计算单元,用于根据所述孔隙度φ、渗透率K、几何平均值T2g计算特征参数;
优化可动水波谱单元,用于将所述核磁共振T2谱划分为束缚水波谱和可动水波谱,通过所述特征参数重构所述可动水波谱,获得优化可动水波谱;
组合单元,用于将所述优化可动水波谱与所述束缚水波谱进行组合,即获得饱含水状态下的核磁共振T2谱。
8.根据权利要求7所述的核磁共振T2谱含烃校正系统,其中,所述特征参数包括所述核磁共振T2谱的截止值T2C、峰值T2f、右边界T2R。
9.根据权利要求8所述的核磁共振T2谱含烃校正系统,其中,根据孔隙度φ、渗透率K、几何平均值T2g通过公式(1)、公式(2)、公式(3)计算峰值T2f、右边界T2R:
T2f=a1*(T2g)2+a2*(T2g)+a3 (1)
T2R=b1*T2s+b2 (2)
其中,T2s为孔隙结构特征参数,a1、a2、a3、b1、b2为公式拟合系数。
10.根据权利要求8所述的核磁共振T2谱含烃校正系统,其中,所述通过所述特征参数重构所述可动水波谱,获得优化可动水波谱包括:
根据所述峰值T2f与截止值T2C的相对位置确定优化可动水波谱的波峰位置和/或波谱形态:
当峰值T2f位于截止值T2C左侧时,所述可动水波谱的形态为单斜分布;当峰值T2f在截止值T2C右侧时,所述可动水波谱的波峰位置为峰值T2f,形态为正态分布。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810260213.6A CN110306977B (zh) | 2018-03-27 | 2018-03-27 | 核磁共振t2谱含烃校正方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810260213.6A CN110306977B (zh) | 2018-03-27 | 2018-03-27 | 核磁共振t2谱含烃校正方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110306977A true CN110306977A (zh) | 2019-10-08 |
CN110306977B CN110306977B (zh) | 2021-10-29 |
Family
ID=68073733
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810260213.6A Active CN110306977B (zh) | 2018-03-27 | 2018-03-27 | 核磁共振t2谱含烃校正方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110306977B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112326710A (zh) * | 2020-11-10 | 2021-02-05 | 中国石油大学(华东) | 一种泥页岩饱和水后各离心状态核磁共振t2谱校正方法 |
CN113720869A (zh) * | 2021-08-26 | 2021-11-30 | 中国地质大学(北京) | 油基泥浆侵入低孔隙度水层核磁共振t2几何均值校正方法与系统 |
Citations (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2000026696A1 (en) * | 1998-10-30 | 2000-05-11 | Numar Corporation | Nmr logging apparatus and methods for fluid typing |
WO2001013142A1 (en) * | 1999-08-13 | 2001-02-22 | Halliburton Energy Services, Inc. | Nmr apparatus and methods for measuring volumes of hydrocarbon gas and oil |
WO2003025625A1 (en) * | 2001-09-19 | 2003-03-27 | Halliburton Energy Services, Inc. | Method and system for using conventional core data to calibrate bound water volumes derived from nmr logs |
CN1570670A (zh) * | 2004-05-08 | 2005-01-26 | 大庆石油管理局 | 利用岩石激发极化谱确定孔隙结构和地层渗透率的方法 |
WO2005119300A1 (en) * | 2004-05-27 | 2005-12-15 | Baker Hughes Incorporated | Method of detecting, quantifying and correcting borehole contaminations from multi-frequency, multi-sensitive volume nmr logging data |
CN1869733A (zh) * | 2005-05-27 | 2006-11-29 | 中国石油天然气股份有限公司 | 核磁共振测井t2谱t2截止值的确定方法 |
CN102141637A (zh) * | 2010-01-28 | 2011-08-03 | 中国石油天然气股份有限公司 | 一种利用核磁共振测井资料连续定量评价储集层孔隙结构的方法 |
CN103487837A (zh) * | 2013-09-13 | 2014-01-01 | 同济大学 | 拟饱含水核磁共振自旋回波信号的分解与合成方法 |
CN104698020A (zh) * | 2013-12-06 | 2015-06-10 | 中国石油天然气股份有限公司 | 一种松散岩心微观孔隙结构特征参数的采集处理方法 |
CN104863574A (zh) * | 2014-02-20 | 2015-08-26 | 中国石油化工股份有限公司 | 一种适用于致密砂岩储层的流体识别方法 |
CN106050225A (zh) * | 2016-06-06 | 2016-10-26 | 中国石油天然气集团公司 | 一种核磁共振测井t2谱100%纯水谱的确定方法 |
US20160370492A1 (en) * | 2015-02-26 | 2016-12-22 | Halliburton Energy Services, Inc. | Methods and systems employing nmr-based prediction of pore throat size distributions |
CN106351652A (zh) * | 2016-11-25 | 2017-01-25 | 中国地质大学(北京) | 一种含烃储集层核磁共振测井t2谱形态校正方法 |
CN106501144A (zh) * | 2016-09-13 | 2017-03-15 | 中国石油大学(华东) | 一种基于核磁共振双截止值的致密砂岩渗透率计算方法 |
-
2018
- 2018-03-27 CN CN201810260213.6A patent/CN110306977B/zh active Active
Patent Citations (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2000026696A1 (en) * | 1998-10-30 | 2000-05-11 | Numar Corporation | Nmr logging apparatus and methods for fluid typing |
WO2001013142A1 (en) * | 1999-08-13 | 2001-02-22 | Halliburton Energy Services, Inc. | Nmr apparatus and methods for measuring volumes of hydrocarbon gas and oil |
WO2003025625A1 (en) * | 2001-09-19 | 2003-03-27 | Halliburton Energy Services, Inc. | Method and system for using conventional core data to calibrate bound water volumes derived from nmr logs |
US20030094946A1 (en) * | 2001-09-19 | 2003-05-22 | Halliburton Energy Services, Inc. | Method and system for using conventional core data to calibrate bound water volumes derived from NMR logs |
CN1570670A (zh) * | 2004-05-08 | 2005-01-26 | 大庆石油管理局 | 利用岩石激发极化谱确定孔隙结构和地层渗透率的方法 |
WO2005119300A1 (en) * | 2004-05-27 | 2005-12-15 | Baker Hughes Incorporated | Method of detecting, quantifying and correcting borehole contaminations from multi-frequency, multi-sensitive volume nmr logging data |
CN1869733A (zh) * | 2005-05-27 | 2006-11-29 | 中国石油天然气股份有限公司 | 核磁共振测井t2谱t2截止值的确定方法 |
CN102141637A (zh) * | 2010-01-28 | 2011-08-03 | 中国石油天然气股份有限公司 | 一种利用核磁共振测井资料连续定量评价储集层孔隙结构的方法 |
CN103487837A (zh) * | 2013-09-13 | 2014-01-01 | 同济大学 | 拟饱含水核磁共振自旋回波信号的分解与合成方法 |
CN104698020A (zh) * | 2013-12-06 | 2015-06-10 | 中国石油天然气股份有限公司 | 一种松散岩心微观孔隙结构特征参数的采集处理方法 |
CN104863574A (zh) * | 2014-02-20 | 2015-08-26 | 中国石油化工股份有限公司 | 一种适用于致密砂岩储层的流体识别方法 |
US20160370492A1 (en) * | 2015-02-26 | 2016-12-22 | Halliburton Energy Services, Inc. | Methods and systems employing nmr-based prediction of pore throat size distributions |
CN106050225A (zh) * | 2016-06-06 | 2016-10-26 | 中国石油天然气集团公司 | 一种核磁共振测井t2谱100%纯水谱的确定方法 |
CN106501144A (zh) * | 2016-09-13 | 2017-03-15 | 中国石油大学(华东) | 一种基于核磁共振双截止值的致密砂岩渗透率计算方法 |
CN106351652A (zh) * | 2016-11-25 | 2017-01-25 | 中国地质大学(北京) | 一种含烃储集层核磁共振测井t2谱形态校正方法 |
Non-Patent Citations (15)
Title |
---|
LAI, J ET AL.: "Fractal analysis of tight shaly sandstones using nuclear magnetic resonance measurements", 《AAPG BULLETIN》 * |
ZOU, YL ET AL.: "A Rapid NMR T-2 Inversion Method Based on Norm Smoothing", 《APPLIED MAGNETIC RESONANCE》 * |
刘堂宴等: "核磁共振谱的岩石孔喉结构分析", 《石油地球物理勘探》 * |
刘美杰等: "江苏油田核磁共振录井综合应用研究", 《中国优秀硕士学位论文全文数据库基础科学辑》 * |
张春雷等: "核磁共振测井岩石实验变T2截止值及储层评价应用研究", 《中国优秀硕士学位论文全文数据库基础科学辑》 * |
张晋言等: "核磁共振测井在深层砂砾岩孔隙结构及有效性评价中的应用", 《测井技术》 * |
肖亮等: "利用核磁共振测井资料评价储集层孔隙结构的讨论", 《新疆石油地质》 * |
肖飞等: "核磁共振测井应用研究——碎屑岩储层孔隙结构、动态T_2截止值、稠油识别研究", 《中国优秀硕士学位论文全文数据库基础科学辑》 * |
肖飞等: "核磁共振测井连续表征储层孔隙结构方法研究", 《石油天然气学报》 * |
苏俊磊等: "应用核磁共振测井资料评价储层孔隙结构的改进方法", 《吉林大学学报(地球科学版)》 * |
许明标等: "《现代储层保护技术》", 31 December 2016, 中国地质大学出版社 * |
谭茂金等: "有机页岩核磁共振实验分析与解释模型建立", 《2014年中国地球科学联合学术年会》 * |
邵维志等: "利用T2谱形态确定T2截止值的方法探索", 《测井技术》 * |
邵维志等: "核磁共振测井在储层孔隙结构评价中的应用", 《测井技术》 * |
魏佳音等: "P型核磁共振储层参数计算及流体识别评价方法", 《中国优秀硕士学位论文全文数据库基础科学辑》 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112326710A (zh) * | 2020-11-10 | 2021-02-05 | 中国石油大学(华东) | 一种泥页岩饱和水后各离心状态核磁共振t2谱校正方法 |
CN112326710B (zh) * | 2020-11-10 | 2022-06-03 | 中国石油大学(华东) | 一种泥页岩饱和水后各离心状态核磁共振t2谱校正方法 |
CN113720869A (zh) * | 2021-08-26 | 2021-11-30 | 中国地质大学(北京) | 油基泥浆侵入低孔隙度水层核磁共振t2几何均值校正方法与系统 |
Also Published As
Publication number | Publication date |
---|---|
CN110306977B (zh) | 2021-10-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Wang et al. | A comprehensive method for studying pore structure and seepage characteristics of coal mass based on 3D CT reconstruction and NMR | |
CN103353462B (zh) | 一种基于核磁共振成像的岩石非均质性定量评价方法 | |
CN106050225B (zh) | 一种核磁共振测井t2谱100%纯水谱的确定方法 | |
WO2017024700A1 (zh) | 一种计算烃源岩中有机碳含量的装置 | |
CN102998322A (zh) | 恒定梯度场核磁共振岩样分析方法及仪器 | |
CN104819923A (zh) | 基于核磁共振的低渗透砂岩储层孔隙结构定量反演方法 | |
Wu et al. | Characterization of pore structure of tight sandstone reservoirs based on fractal analysis of NMR echo data | |
CN102262041A (zh) | 一种基于多谱孔隙结构分析的饱和度确定方法 | |
CN110231272A (zh) | 致密砂岩孔径与核磁共振t2值转换关系的确定方法及系统 | |
CN104453874A (zh) | 一种基于核磁共振的砂砾岩储层含油饱和度的计算方法 | |
Zou et al. | Inversion of nuclear magnetic resonance echo data based on maximum entropy | |
CN106930754A (zh) | 一种储层自适应的核磁共振测井数据采集方法 | |
CN105986813A (zh) | 准致密储层快速评价方法及准致密储层多指标评价方法 | |
CN110306977A (zh) | 核磁共振t2谱含烃校正方法及系统 | |
CN106644879B (zh) | 一种确定岩心不同孔隙组分渗透率贡献值的方法及装置 | |
Chen et al. | Contributions of pore-throat size distribution to reservoir quality and fluid distribution from NMR and MIP in tight sandy conglomerate reservoirs | |
CN106096081B (zh) | 缝洞型底水油藏的储量估算方法 | |
CN106351652B (zh) | 一种含烃储集层核磁共振测井t2谱形态校正方法 | |
CN106290103A (zh) | 一种页岩气储层中粘土微孔孔隙度的测定方法 | |
Deng et al. | Integrated petrophysical log characterization for tight carbonate reservoir effectiveness: A case study from the Longgang area, Sichuan Basin, China | |
CN109541689B (zh) | 一种基于反射波能量特征的介质密实度评价方法 | |
CN116027453A (zh) | 一种水合物混合层饱和度定量评价方法及装置 | |
CN106324008B (zh) | 核磁共振测井波谱的校正方法 | |
CN109856688A (zh) | 基于核磁测井双tw极化增强法的流体性质识别方法 | |
CN106156456A (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 |