CN110306977B - 核磁共振t2谱含烃校正方法及系统 - Google Patents
核磁共振t2谱含烃校正方法及系统 Download PDFInfo
- Publication number
- CN110306977B CN110306977B CN201810260213.6A CN201810260213A CN110306977B CN 110306977 B CN110306977 B CN 110306977B CN 201810260213 A CN201810260213 A CN 201810260213A CN 110306977 B CN110306977 B CN 110306977B
- Authority
- CN
- China
- Prior art keywords
- spectrum
- magnetic resonance
- nuclear magnetic
- value
- movable water
- 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
- 238000005481 NMR spectroscopy Methods 0.000 title claims abstract description 94
- 239000004215 Carbon black (E152) Substances 0.000 title claims abstract description 27
- 238000000034 method Methods 0.000 title claims abstract description 27
- 229930195733 hydrocarbon Natural products 0.000 title claims abstract description 26
- 238000012937 correction Methods 0.000 title claims abstract description 22
- 150000002430 hydrocarbons Chemical class 0.000 title claims abstract description 20
- 238000001228 spectrum Methods 0.000 claims abstract description 179
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 119
- 239000011148 porous material Substances 0.000 claims abstract description 34
- 230000035699 permeability Effects 0.000 claims abstract description 22
- 239000011435 rock Substances 0.000 claims abstract description 19
- 238000002474 experimental method Methods 0.000 claims abstract description 16
- 229920006395 saturated elastomer Polymers 0.000 claims abstract description 14
- 239000010779 crude oil Substances 0.000 claims description 22
- 230000003595 spectral effect Effects 0.000 claims description 19
- 239000003921 oil Substances 0.000 claims description 12
- 125000001183 hydrocarbyl group Chemical group 0.000 claims description 6
- 230000000694 effects Effects 0.000 abstract description 7
- 238000001225 nuclear magnetic resonance method Methods 0.000 abstract description 2
- 238000009792 diffusion process Methods 0.000 description 4
- 239000012530 fluid Substances 0.000 description 4
- 238000009736 wetting Methods 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 235000021167 banquet Nutrition 0.000 description 1
- 238000012512 characterization method Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 239000000295 fuel oil Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000000655 nuclear magnetic resonance spectrum Methods 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 238000011158 quantitative evaluation Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
Images
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)
- Engineering & Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Geology (AREA)
- Mining & Mineral Resources (AREA)
- General Physics & Mathematics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Environmental & Geological Engineering (AREA)
- Business, Economics & Management (AREA)
- High Energy & Nuclear Physics (AREA)
- Health & Medical Sciences (AREA)
- Animal Husbandry (AREA)
- Primary Health Care (AREA)
- Strategic Management (AREA)
- Tourism & Hospitality (AREA)
- Human Resources & Organizations (AREA)
- General Business, Economics & Management (AREA)
- General Health & Medical Sciences (AREA)
- Theoretical Computer Science (AREA)
- Economics (AREA)
- Marine Sciences & Fisheries (AREA)
- Marketing (AREA)
- Remote Sensing (AREA)
- Agronomy & Crop Science (AREA)
- Geophysics (AREA)
- Fluid Mechanics (AREA)
- Geochemistry & Mineralogy (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Condensed Matter Physics & Semiconductors (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)进行重建:
本发明将含烃的核磁共振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 (6)
1.一种核磁共振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为公式拟合系数。
2.根据权利要求1所述的核磁共振T2谱含烃校正方法,其中,所述通过所述特征参数重构所述可动水波谱,获得优化可动水波谱包括:
根据所述峰值T2f与截止值T2C的相对位置确定优化可动水波谱的波峰位置和/或波谱形态。
3.根据权利要求2所述的核磁共振T2谱含烃校正方法,其中,根据所述峰值T2f与截止值T2C的相对位置确定优化可动水波谱的波峰位置和/或波谱形态包括:
当峰值T2f位于截止值T2C左侧时,所述可动水波谱的形态为单斜分布;当峰值T2f在截止值T2C右侧时,所述可动水波谱的波峰位置为峰值T2f,形态为正态分布。
4.根据权利要求1所述的核磁共振T2谱含烃校正方法,其中,还包括:
采集储层油样进行原油性质实验,获得原油粘度,当原油为轻质油时,对核磁共振T2谱进行校正。
5.一种核磁共振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为公式拟合系数。
6.根据权利要求5所述的核磁共振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 CN110306977A (zh) | 2019-10-08 |
CN110306977B true 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) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112326710B (zh) * | 2020-11-10 | 2022-06-03 | 中国石油大学(华东) | 一种泥页岩饱和水后各离心状态核磁共振t2谱校正方法 |
CN113720869B (zh) * | 2021-08-26 | 2023-03-31 | 中国地质大学(北京) | 油基泥浆侵入低孔隙度水层核磁共振t2几何均值校正方法与系统 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1570670A (zh) * | 2004-05-08 | 2005-01-26 | 大庆石油管理局 | 利用岩石激发极化谱确定孔隙结构和地层渗透率的方法 |
CN102141637A (zh) * | 2010-01-28 | 2011-08-03 | 中国石油天然气股份有限公司 | 一种利用核磁共振测井资料连续定量评价储集层孔隙结构的方法 |
CN104698020A (zh) * | 2013-12-06 | 2015-06-10 | 中国石油天然气股份有限公司 | 一种松散岩心微观孔隙结构特征参数的采集处理方法 |
Family Cites Families (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6366087B1 (en) * | 1998-10-30 | 2002-04-02 | George Richard Coates | NMR logging apparatus and methods for fluid typing |
US6661226B1 (en) * | 1999-08-13 | 2003-12-09 | Halliburton Energy Services, Inc. | NMR apparatus and methods for measuring volumes of hydrocarbon gas and oil |
AU2002341788B2 (en) * | 2001-09-19 | 2008-01-10 | Halliburton Energy Services, Inc. | Method and system for using conventional core data to calibrate bound water volumes derived from NMR logs |
US7227355B2 (en) * | 2004-05-27 | 2007-06-05 | Baker Hughes Incorporated | Method of detecting, quantifying and correcting borehole contaminations from multi-frequency, multi-sensitive-volume NMR logging data |
CN100349013C (zh) * | 2005-05-27 | 2007-11-14 | 中国石油天然气股份有限公司 | 核磁共振测井t2谱t2截止值的确定方法 |
CN103487837B (zh) * | 2013-09-13 | 2016-11-23 | 同济大学 | 拟饱含水核磁共振自旋回波信号的分解与合成方法 |
CN104863574B (zh) * | 2014-02-20 | 2018-03-13 | 中国石油化工股份有限公司 | 一种适用于致密砂岩储层的流体识别方法 |
US10324222B2 (en) * | 2015-02-26 | 2019-06-18 | Halliburton Energy Services, Inc. | Methods and systems employing NMR-based prediction of pore throat size distributions |
CN106050225B (zh) * | 2016-06-06 | 2019-03-12 | 中国石油天然气集团公司 | 一种核磁共振测井t2谱100%纯水谱的确定方法 |
CN106501144B (zh) * | 2016-09-13 | 2018-10-26 | 中国石油大学(华东) | 一种基于核磁共振双截止值的致密砂岩渗透率计算方法 |
CN106351652B (zh) * | 2016-11-25 | 2019-03-29 | 中国地质大学(北京) | 一种含烃储集层核磁共振测井t2谱形态校正方法 |
-
2018
- 2018-03-27 CN CN201810260213.6A patent/CN110306977B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1570670A (zh) * | 2004-05-08 | 2005-01-26 | 大庆石油管理局 | 利用岩石激发极化谱确定孔隙结构和地层渗透率的方法 |
CN102141637A (zh) * | 2010-01-28 | 2011-08-03 | 中国石油天然气股份有限公司 | 一种利用核磁共振测井资料连续定量评价储集层孔隙结构的方法 |
CN104698020A (zh) * | 2013-12-06 | 2015-06-10 | 中国石油天然气股份有限公司 | 一种松散岩心微观孔隙结构特征参数的采集处理方法 |
Non-Patent Citations (4)
Title |
---|
P型核磁共振储层参数计算及流体识别评价方法;魏佳音等;《中国优秀硕士学位论文全文数据库基础科学辑》;20131231(第12期);第A011-215页 * |
核磁共振测井岩石实验变T2截止值及储层评价应用研究;张春雷等;《中国优秀硕士学位论文全文数据库基础科学辑》;20161231(第12期);第A011-123页 * |
核磁共振测井应用研究——碎屑岩储层孔隙结构、动态T_2截止值、稠油识别研究;肖飞等;《中国优秀硕士学位论文全文数据库基础科学辑》;20130131(第1期);第A011-215页 * |
江苏油田核磁共振录井综合应用研究;刘美杰等;《中国优秀硕士学位论文全文数据库基础科学辑》;20110430(第4期);第A011-183页 * |
Also Published As
Publication number | Publication date |
---|---|
CN110306977A (zh) | 2019-10-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Kuchuk et al. | Pressure behavior of laterally composite reservoirs | |
CN110306977B (zh) | 核磁共振t2谱含烃校正方法及系统 | |
Wu et al. | Characterization of pore structure of tight sandstone reservoirs based on fractal analysis of NMR echo data | |
CN104849763B (zh) | 趋于低频趋势的多井测井曲线标准化方法 | |
CN106897531B (zh) | 一种低渗透石灰岩储层渗透率的定量评价方法 | |
Chen et al. | Efficient incorporation of global effects in upscaled models of two-phase flow and transport in heterogeneous formations | |
CN110687597B (zh) | 一种基于联合字典的波阻抗反演方法 | |
Wang et al. | Anchoring multi‐scale models to micron‐scale imaging of multiphase flow in rocks | |
Khalili et al. | Formation factor for heterogeneous carbonate rocks using multi-scale Xray-CT images | |
Saidian | Effect of rock composition and texture on pore size distributions in shales: Applications in low field nuclear magnetic resonance | |
Tripathi et al. | Exporting petrophysical properties of sandstones from thin section image analysis | |
Meng et al. | An NMR log echo data de-noising method based on the wavelet packet threshold algorithm | |
CN111827995B (zh) | 一种基于核磁共振的超压砾岩储层渗透率的计算方法 | |
Gómez et al. | Spectral structure-oriented filtering of seismic data with self-adaptive paths | |
CN114674722A (zh) | 一种不同岩性储层的分类与孔隙结构定量表征方法及装置 | |
Meridji et al. | Fluid identification in complex clastic reservoirs using 2d NMR Maps: A case study from Saudi Arabia | |
Rosa et al. | Use of the NMR diffusivity log to identify and quantify oil and water in carbonate formations | |
Madi et al. | Unconventional Shale Play in Oman: Preliminary Assessment of the Shale Oil/Shale Gas Potential of the Silurian Hot Shale of the Southern Rub al-Khali Basin | |
CN112444864A (zh) | 薄层高分辨率反演方法及存储介质 | |
CN114154283A (zh) | 储层大孔径孔隙度计算方法、装置、电子设备及介质 | |
Hinai et al. | Pore geometry in gas shale reservoirs | |
Woodhouse et al. | Sw and Hydrocarbon Pore Volume Estimates in Shaly Sands—Routine Oil-Based-Mud Core Measurements Compared With Several Log Analysis Models | |
Manchuk et al. | Advances in micromodeling using resistivity borehole images | |
Prince | Textural and Diagenetic Controlson Sandstone Permeability | |
CN111472750B (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 |