CN102998322A - 恒定梯度场核磁共振岩样分析方法及仪器 - Google Patents
恒定梯度场核磁共振岩样分析方法及仪器 Download PDFInfo
- Publication number
- CN102998322A CN102998322A CN2011102724771A CN201110272477A CN102998322A CN 102998322 A CN102998322 A CN 102998322A CN 2011102724771 A CN2011102724771 A CN 2011102724771A CN 201110272477 A CN201110272477 A CN 201110272477A CN 102998322 A CN102998322 A CN 102998322A
- Authority
- CN
- China
- Prior art keywords
- magnetic resonance
- nuclear magnetic
- oil
- sample
- rock
- 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 175
- 239000011435 rock Substances 0.000 title claims abstract description 167
- 238000012284 sample analysis method Methods 0.000 title abstract 2
- 239000012530 fluid Substances 0.000 claims abstract description 114
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 79
- 238000000034 method Methods 0.000 claims abstract description 45
- 238000005259 measurement Methods 0.000 claims abstract description 27
- 238000009826 distribution Methods 0.000 claims abstract description 23
- 239000000523 sample Substances 0.000 claims description 164
- 239000003921 oil Substances 0.000 claims description 110
- 238000005084 2D-nuclear magnetic resonance Methods 0.000 claims description 78
- 238000009792 diffusion process Methods 0.000 claims description 75
- 239000010779 crude oil Substances 0.000 claims description 46
- 230000035699 permeability Effects 0.000 claims description 37
- 229910000831 Steel Inorganic materials 0.000 claims description 18
- 239000010959 steel Substances 0.000 claims description 18
- 238000000685 Carr-Purcell-Meiboom-Gill pulse sequence Methods 0.000 claims description 15
- 230000008569 process Effects 0.000 claims description 13
- 229920006395 saturated elastomer Polymers 0.000 claims description 13
- 230000000306 recurrent effect Effects 0.000 claims description 11
- 238000012937 correction Methods 0.000 claims description 10
- 238000011156 evaluation Methods 0.000 claims description 9
- XEEYBQQBJWHFJM-UHFFFAOYSA-N Iron Chemical compound [Fe] XEEYBQQBJWHFJM-UHFFFAOYSA-N 0.000 claims description 8
- 238000012360 testing method Methods 0.000 claims description 8
- 238000000354 decomposition reaction Methods 0.000 claims description 7
- 230000007246 mechanism Effects 0.000 claims description 7
- 230000008859 change Effects 0.000 claims description 6
- 230000010354 integration Effects 0.000 claims description 6
- 238000005086 pumping Methods 0.000 claims description 6
- 238000013461 design Methods 0.000 claims description 5
- 229910052742 iron Inorganic materials 0.000 claims description 4
- 238000004026 adhesive bonding Methods 0.000 claims description 3
- 230000003321 amplification Effects 0.000 claims description 3
- 230000005540 biological transmission Effects 0.000 claims description 3
- 239000003990 capacitor Substances 0.000 claims description 3
- 239000000696 magnetic material Substances 0.000 claims description 3
- 238000003913 materials processing Methods 0.000 claims description 3
- 239000000203 mixture Substances 0.000 claims description 3
- 238000003199 nucleic acid amplification method Methods 0.000 claims description 3
- 238000012545 processing Methods 0.000 claims description 3
- 238000007789 sealing Methods 0.000 claims description 3
- 238000012546 transfer Methods 0.000 claims description 3
- 238000004804 winding Methods 0.000 claims description 3
- 241000826860 Trapezium Species 0.000 claims description 2
- 238000013459 approach Methods 0.000 claims description 2
- 238000001228 spectrum Methods 0.000 abstract description 26
- 238000002495 two-dimensional nuclear magnetic resonance spectrum Methods 0.000 abstract 1
- DCAYPVUWAIABOU-UHFFFAOYSA-N hexadecane Chemical compound CCCCCCCCCCCCCCCC DCAYPVUWAIABOU-UHFFFAOYSA-N 0.000 description 16
- 239000011148 porous material Substances 0.000 description 15
- 238000011161 development Methods 0.000 description 11
- 230000000694 effects Effects 0.000 description 11
- 239000007789 gas Substances 0.000 description 10
- 239000011159 matrix material Substances 0.000 description 10
- 238000005516 engineering process Methods 0.000 description 9
- 238000002474 experimental method Methods 0.000 description 6
- 238000005520 cutting process Methods 0.000 description 5
- 238000006073 displacement reaction Methods 0.000 description 5
- 239000003208 petroleum Substances 0.000 description 5
- 238000011160 research Methods 0.000 description 4
- 239000004215 Carbon black (E152) Substances 0.000 description 3
- 208000002925 dental caries Diseases 0.000 description 3
- 229930195733 hydrocarbon Natural products 0.000 description 3
- 150000002430 hydrocarbons Chemical class 0.000 description 3
- 229910052739 hydrogen Inorganic materials 0.000 description 3
- 239000001257 hydrogen Substances 0.000 description 3
- -1 G&W Substances 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 2
- 239000010426 asphalt Substances 0.000 description 2
- 230000015572 biosynthetic process Effects 0.000 description 2
- 238000011088 calibration curve Methods 0.000 description 2
- 230000001808 coupling effect Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000005284 excitation Effects 0.000 description 2
- 239000012467 final product Substances 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- 230000003993 interaction Effects 0.000 description 2
- 230000005415 magnetization Effects 0.000 description 2
- 230000009467 reduction Effects 0.000 description 2
- 239000007787 solid Substances 0.000 description 2
- UFHFLCQGNIYNRP-UHFFFAOYSA-N Hydrogen Chemical compound [H][H] UFHFLCQGNIYNRP-UHFFFAOYSA-N 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000012153 distilled water Substances 0.000 description 1
- 230000005672 electromagnetic field Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000000704 physical effect Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
Images
Classifications
-
- 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
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N24/00—Investigating or analyzing materials by the use of nuclear magnetic resonance, electron paramagnetic resonance or other spin effects
- G01N24/08—Investigating or analyzing materials by the use of nuclear magnetic resonance, electron paramagnetic resonance or other spin effects by using nuclear magnetic resonance
- G01N24/081—Making measurements of geologic samples, e.g. measurements of moisture, pH, porosity, permeability, tortuosity or viscosity
-
- 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/48—NMR imaging systems
- G01R33/54—Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
-
- 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/28—Details of apparatus provided for in groups G01R33/44 - G01R33/64
- G01R33/38—Systems for generation, homogenisation or stabilisation of the main or gradient magnetic field
- G01R33/3806—Open magnet assemblies for improved access to the sample, e.g. C-type or U-type magnets
-
- 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/28—Details of apparatus provided for in groups G01R33/44 - G01R33/64
- G01R33/38—Systems for generation, homogenisation or stabilisation of the main or gradient magnetic field
- G01R33/383—Systems for generation, homogenisation or stabilisation of the main or gradient magnetic field using permanent magnets
-
- 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/445—MR involving a non-standard magnetic field B0, e.g. of low magnitude as in the earth's magnetic field or in nanoTesla spectroscopy, comprising a polarizing magnetic field for pre-polarisation, B0 with a temporal variation of its magnitude or direction such as field cycling of B0 or rotation of the direction of B0, or spatially inhomogeneous B0 like in fringe-field MR or in stray-field imaging
-
- 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/48—NMR imaging systems
- G01R33/483—NMR imaging systems with selection of signals or spectra from particular regions of the volume, e.g. in vivo spectroscopy
- G01R33/4833—NMR imaging systems with selection of signals or spectra from particular regions of the volume, e.g. in vivo spectroscopy using spatially selective excitation of the volume of interest, e.g. selecting non-orthogonal or inclined slices
- G01R33/4835—NMR imaging systems with selection of signals or spectra from particular regions of the volume, e.g. in vivo spectroscopy using spatially selective excitation of the volume of interest, e.g. selecting non-orthogonal or inclined slices of multiple slices
Landscapes
- Physics & Mathematics (AREA)
- High Energy & Nuclear Physics (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Analytical Chemistry (AREA)
- Geochemistry & Mineralogy (AREA)
- Health & Medical Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Geophysics (AREA)
- Signal Processing (AREA)
- Condensed Matter Physics & Semiconductors (AREA)
- Remote Sensing (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
本发明涉及一种恒定梯度场核磁共振岩样分析方法及仪器;在由磁体产生的恒定梯度场下进行核磁共振测量获得测量数据;将测量的核磁共振数据转化成核磁共振二维谱D-T2;对恒定梯度场的标准样进行测量、反演得到标准样核磁共振二维谱D-T2;测量岩样,获得岩样中流体核磁共振二维谱D-T2,根据实际测量岩样的核磁共振二维谱D-T2进行流体类型识别;根据岩样中流体的核磁共振二维谱D-T2计算岩样的流体性质和岩石物性参数;对岩样进行单切片扫描,获得岩样的局部油、水饱和度;对岩样进行连续的切片扫描,得岩样轴向油、水饱和度分布,可动流体饱和度分布;本方法可更好的评价储层和识别流体。
Description
技术领域
本发明涉及一种恒定梯度场核磁共振岩样分析仪器和利用核磁共振二维谱技术识别流体类型和测量岩石物性参数的方法。
背景技术
近些年,核磁共振技术在石油勘探开发领域得到高度重视和发展,已广泛应用于核磁共振录井、核磁共振测井和低渗透储层评价等方面,从而推动了核磁共振仪器及其测量方法在石油勘探开发中快速发展。
中国专利200410083878.2中描述了用岩石核磁共振弛豫信号测量地层岩石物性的设备及测量方法。该设备包括磁体、探头、前置放大器、功率放大器、核磁共振控制器和控制计算机等部分。该设备通过测量获得一维的弛豫时间T2谱,可以快速给出孔隙度、渗透率、可动流体饱和度和含油饱和度等岩石物性参数。
核磁共振仪器除了提供弛豫时间T2以外,还可以提供流体的扩散系数D。目前所有的可以测量流体扩散系数的实验室核磁共振仪器采用的是电梯度,利用电流的方式来获得梯度磁场,利用电流来控制梯度的大小和方向,在医学的核磁共振仪器方面也广泛采用该种方式。但是,由于在石油勘探开发中原油的扩散系数较小,所以需要很大的梯度场,才能测量到原油的扩散系数,对油和水进行识别。梯度场的大小是由瞬间电流的大小决定的,但是,瞬间电流是不可能无限放大的,它的大小将受到仪器的限制。特别是产生瞬间电流越大,将会导致***越复杂。
目前,核磁共振仪器动态驱替实验都是驱替后,将岩心取出,然 后进行核磁共振测量。而核磁共振高温高压核磁共振实验就根本无法进行。主要是由于核磁共振仪器的电磁场无法穿越金属壳,在探头内部不能形成梯度磁场,从而无法进行扩散编码,测量流体的扩散系数。这些都抑制了核磁共振实验室技术的发展。
因此,急需可以测量扩散系数,能够进行高温高压实验的恒定梯度场核磁共振岩样分析仪以及核磁共振二维谱识别流体类型和评价岩石物性参数的方法。
发明内容
本发明的目的是提供一种恒定梯度场核磁共振岩样分析仪及测量方法,以测量扩散系数、弛豫时间、原油粘度、孔隙度、渗透率、含油、水饱和度及分布、可动流体饱和度等岩石物性参数和流体性质参数,解决石油勘探开发中的流体识别和储层评价中的难点问题,提高室内核磁共振装备与技术水平,满足油气田勘探开发需求,促进核磁共振岩样技术在石油勘探开发领域快速发展。
本发明采用的技术方案是:恒定梯度场核磁共振岩样分析仪由测量装置1和控制***4构成;测量装置1由梯度磁体2、探头3和温控***构成;梯度磁体2是一个封闭的腔体,探头2位于梯度磁体腔体的中心位置,探头2与控制***4通过电缆与控制***4的T/R匹配电路10连接;控制***4通过USB口与计算机连接。
恒定梯度场核磁共振岩样分析仪的梯度磁体2由轭板14、磁钢15、极板16、侧轭板17构成;轭板14和侧轭板17组成梯度磁体腔体,磁钢15与轭板14连接,极板16截面为梯形结构,与磁钢15连接;轭板、极板、侧轭板采用电工纯铁DT4C材料加工制成,磁钢采用2:17SmCo磁性材料制成。磁钢15与轭板14,极板16与磁钢采用BJ-39胶粘接。梯度磁体2会产生纵向均匀而横向存在梯度的磁场B0形成样 品探测区,磁场方向沿着z轴,梯度方向沿x轴。两极板的角度在10~30°之间,其两极板间形成的气隙在150~364mm之间。(详见专利200910092839.1公布的一种用于恒定梯度场岩样分析仪的梯度永磁体)。梯度磁体2外表面是一组温控***6,会对梯度磁体进行10~60度恒温控制,减少温度对仪器测量结果的影响。
恒定梯度场核磁共振岩样分析仪探头3由一个不导电的圆形样品管18,回路线圈19构成;回路线圈19缠绕在样品管18上;回路线圈19为螺旋管线圈、鞍形线圈或其它适用的线圈。探头可作为发射机将激励信号以射频的方式发射到待测样品中和作为接收机接收核磁共振信息。探头3内部是样品探测区20。特别是恒定梯度场核磁共振岩样分析仪可以放置核磁共振高温高压探头。
恒定梯度场核磁共振岩样分析仪控制***4由主控制器7,频率发生器8,功率放大器9,T/R匹配电路10,前置放大器11、接收器12和带有缓存A/D转换器13构成。主控制器7主要有3个功能,产生并发射用以激发核磁共振信号的脉冲序列,接收核磁共振回波信号以及将该信号计算处理上传至计算机,主控制器7通过USB接口与计算机5进行互联及数据传输。频率发生器8与主控制器7进行采用电缆互联,该频率发生器将主控制器产生的核磁共振激励信号进行电流驱动,所产生的信号送至功率放大器9进行功率放大。功率放大器9与频率发生器8采用电缆互联,该功率放大器可将信号放大到几百瓦,用来激发待测样品的核磁共振信息。探头3以及发射器/接收器(T/R)匹配电路10通常包括一个谐振电容器、一个T/R转换开关和阻抗匹配电路,与功率放大器9和前置放大器11采用电缆连接;为获得更准确的回波信号数据,要求在射频激励信号发射完之后尽量短的时间内接收回波信号,通过T/R发射接收转换开关可以快速有效的完成这一转 换过程;阻抗匹配电路用来对射频电路进行阻抗匹配。经过前置放大器将核磁共振回波信号放大约70DB,其幅值经放大后可以接收器12采集接收,然后由一个带有缓存器的A/D转换器13输出到主控制器7上,以提供用于进一步使用和分析的所需输出数据。前置放大器11与接收器12,接收器12与A/D转换器13之间采用电缆连接。
同时,该***还包括辅助设备计算机5,对控制***进行控制,从而来控制整套仪器。计算机5与控制***4中的主控制器7之间连接,采用USB口的方式。
恒定梯度场核磁共振岩样分析仪的测量方法分为:A步骤:恒定梯度场脉冲序列及参数设计;B步骤:获得核磁共振二维谱;C步骤:恒定梯度场核磁共振定标;D步骤:识别流体类型;E步骤:计算岩石物性参数;F步骤:获得岩样轴向饱和度分布。
A步骤:恒定梯度场脉冲序列及参数设计
脉冲序列是核磁共振技术的灵魂,决定了核磁共振的应用领域。本发明涉及一个CGMF-CPMG脉冲序列。90度脉冲和180度脉冲之间的时间间隔为脉冲间隔tau,单位:us,Nk为tau的个数,单位:个。当k=1时,脉冲间隔tau的值为 ,一个90度脉冲后面跟着P1个180度脉冲,此时CGMF-CPMG脉冲序列就简化成CPMG脉冲序列;当k=2时,一个90度脉冲后面跟着P2个180度脉冲……一直到等于k时,一个90度脉冲后面跟着Pk个180度脉冲。该技术的原理是由于分子的布朗运动造成核磁共振信号强度的衰减,而这种衰减与流体分子的扩散系数关联。90度脉冲的作用是使宏观磁化矢量旋转到横向平面上,连续的180度脉冲的作用是重新聚焦除了扩散衰减导致的相散以外的信号,获得自旋回波衰减T2的其它机制的信息。通过在恒定梯度场下改变CGMF-CPMG脉冲序列中tau的值,来产生扩散效应。恒定梯度场是由 磁体产生的,一直作用于样品探测区。
根据被测样品中饱和流体设计CGMF-CPMG脉冲序列中回波串k,脉冲间隔tau,180度脉冲个数Ne,等待时间RD参数,进行核磁共振测量获得测量数据;k值为2~8,tau的值分布在150us~20000us之间,Ne的值分布在128~30720之间,等待时间分布在2000ms~12000ms之间;
B步骤:获得核磁共振二维谱;
核磁共振技术在石油工业中的应用,主要是利用核磁共振仪来测量岩石孔隙中饱和流体氢核的信号。这些应用不仅与流体本身有关,同时还受到岩石的孔隙结构的影响。流体中氢核的弛豫时间的影响因素包括体弛豫,表面弛豫和扩散弛豫。体弛豫T2B是自旋核与晶格之间能量转移的过程,弛豫时间受自旋-自旋偶合作用和共振频率的影响。表面弛豫T2S是多孔介质中孔隙流体分子间的相互作用,弛豫时间和表面弛豫率以及孔隙的比表面有关。扩散弛豫和分子的扩散系数D,磁场梯度G,脉冲间隔tau有关。因此,在岩石孔隙中T2弛豫时间的表达式为:
式中:T2app为表观弛豫时间,T2B是体弛豫时间,T2S是表面弛豫时间,T2D是扩散弛豫时间。
对T2S有:
对T2D有:
其中:D是流体的扩散系数,γ是旋磁比常数,G是磁场梯度,tau是脉冲时间间隔。
根据步骤A中的恒定梯度场的脉冲序列测量的信号b(t,tau)的一般表达式为:
b(t,tau)=∫∫f(D,T2)k2(t,TE,D,T2)dDdT2+ε (4)
离散形式为:
弛豫时间和扩散系数在对数坐标系轴上均匀分布,方程(5)就简化成:
bik=fljEik,lj+εik (6)
其中:i=1,……nk,k=1,……q,l=1,……p,j=1,……m
i表示第k个回波串的第i个回波,无量纲;
k表示第k个回波串,无量纲;
l表示提前选择的第l个扩散系数,无量纲;
j表示提前选择的第j个弛豫时间,无量纲;
nk是第k个回波串的回波的个数,无量纲;
q是不同tauk的回波串的个数,无量纲;
p表示提前选择的扩散系数的个数,单位:个;
m表示提前选择的弛豫时间的个数,单位:个;
bik表示脉冲间隔为tauk的第k个回波串的第i个回波的振幅,无量纲;
flj表示在扩散系数为Dl和弛豫时间为T2j时的振幅,无量纲;
γ是旋磁比,单位:MHz/T;
G是磁场梯度,单位:Gauss/cm;
tauk是第k个回波串的脉冲间隔,单位:us。
方程(6)的矩阵形式为:
在求解出的核磁共振二维谱的幅度中很可能存在小于零的分量,这在物理意义上来说是不合逻辑的。在解方程7的时候,还必须加上非负约束,即flj≥0。方程7中的矩阵E高度奇异,条件数高,并且实际测的信号伴随着噪声干扰,给反演带来非常大的难题。
本发明中是利用改进的奇异值分解法进行反演,获得核磁共振二维谱。奇异值分解法给出的解是||Ef-b||2最小意义下的一个最优解,但不一定满足非负约束条件。传统的SVD算法通过缩减E迭代求解,等于降低了解的维数,丢掉了振荡最厉害的坐标分量。正是由于缺失了这部分坐标分量,才导致谱图的不连续。为了获得连续的谱图,改进的奇异值分解法将求解Δf的最小二乘解。假设已知一个初始解f0,令b0=Ef0,这样原来的方程可以变为:E(f-f0)=b-b0,即EΔf=Δb。若求得||EΔf-Δb||2最小意义下的最优解Δf,则f0+Δf是||Ef-b||2最小意义下的最优解f。由于求解的是Δf,所以在实现非负约束时,就不再缩小矩 阵E,而只须将f小于零的分量改为零,再重新迭代计算Δf即可,直到f所有坐标分量满足非负约束。由于E在迭代求解的过程中不发生变化,所以求||EΔf-Δb||2最小意义下的最优解无须再进行奇异值分解。与原算法相比较,新算法以计算Δf和Δb代替了矩阵E的奇异值分解过程,不必每次循环都要对矩阵E进行奇异值分解,这样就大大减少了计算量和计算时间。
C步骤:恒定梯度场核磁共振定标
定标是核磁共振技术测量岩样孔隙度不可缺少的步骤。恒定梯度场核磁共振岩样分析仪定标,是根据步骤A和B利用恒定梯度场核磁共振岩样分析仪对恒定梯度场的标准样进行测量、反演得到标准样的核磁共振二维谱D-T2,获得单位体积信号与孔隙度之间的关系线。恒定梯度场核磁共振标准样一共有12个,孔隙度分别为0.5%,1%,2%,3%,6%,9%,12%,15%,18%,21%,24%,27%。每次定标至少选择其中的5个进行测量,主要步骤如下:
(1)把标准样放入恒定梯度场核磁共振岩样分析仪探头中利用提前设计的脉冲序列参数进行测量,得到一套多指数衰减数据;
(2)输入标准样体积、孔隙度,对样品核磁共振信号进行归一化处理,得到样品的单位体积信号量;
(3)重复步骤(1)和(2)测量不同孔隙度的标准样;
(4)软件***将自动进行线性拟合得到核磁共振二维谱定标曲线。单位体积核磁共振信号与孔隙度之间的关系式:
y=ax+b (8)
其中,y代表单位体积核磁共振信号量,x代表核磁共振孔隙度(%),a代表斜率,b代表截距。
D步骤:识别流体类型;
根据步骤A和B,获得实际岩样的核磁共振二维谱D-T2进行流体类型识别。核磁共振二维谱一维是扩散系数D,另一维是弛豫时间T2,所以通过核磁共振二维谱可同时获得被测样品的扩散系数和弛豫时间。由于体相流体,例如原油、天然气和水,不仅弛豫时间T2不同,而且扩散系数D也不同,所以可以利用核磁共振二维谱快速识别流体类型。体相流体水的扩散系数是一常数,与温度有关。气的扩散系数与温度和压力有关。原油的的扩散系数和弛豫时间之间存在线性关系。
Dw(T2)=Dw(T) (9)
Dg(T2)=Dg(T,P) (10)
Do=αT2 (11)
恒定梯度场核磁共振岩样分析仪根据油气水扩散系数和弛豫时间的这种特性,建立核磁共振二维谱解释模板,可以快速识别流体类型。
孔隙中的流体受孔隙表面的影响,弛豫时间和扩散系数都会受到不同程度的影响,导致在一维的弛豫时间T2谱和一维的扩散系数D谱上,部分油和水的信号会互相重叠,无法直接区分。核磁共振二维谱D-T2在一维弛豫时间T2谱的基础上增加了一维扩散系数,从立体空间上可以更好的识别孔隙中流体类型。
E步骤:计算岩石物性参数
恒定梯度场核磁共振岩样分析仪可用于测量流体弛豫和扩散的性质,特别是岩石孔隙中的流体性质。由于这些性质对于油气水是不同的,所以恒定梯度场核磁共振岩样分析仪能够提供油、气和水的比例,以及岩石孔隙度、渗透率、含油饱和度,可动流体饱和度、原油系数等参数,尤其是可以获得岩样的轴向含油饱和度分布和含水饱和度分布。
(1)孔隙度
核磁共振技术直接对被测样品中的氢核(1H)进行探测,可以准确测量岩样孔隙内流体的量,当岩样孔隙内充满流体时,流体量与孔隙体积相等。因此,采用核磁共振技术可以准确测量岩样孔隙度。
由公式(8)得:
利用核磁共振二维谱技术测量岩样孔隙度时,首先,必须保证岩样充分饱和;其次,利用恒定梯度场核磁共振标准样进行定标,获得核磁共振二维谱定标线,即单位体积核磁共振信号与孔隙度之间的关系;最后,将测量饱和的岩样中流体的核磁共振信号带入公式(12),即可获得岩样的核磁共振孔隙度。
(2)含油饱和度
含油饱和度是评价油田储层好坏的重要参数之一,测量含油饱和度的关键就是区分油和水。核磁共振二维谱技术可以测量油和水的弛豫时间和扩散系数,根据核磁共振二维谱解释模板快速识别流体类型,区分油和水。恒定梯度场核磁共振岩样分析仪软件可以利用点击选择式计算含油饱和度。利用鼠标选择核磁共振二维谱中油的区域,电脑将自动计算选择区域体积分与总体积分的比值,获得含油饱和度。以此类推,可以分别得到含油、气、水的饱和度。
(3)可动流体饱和度
弛豫时间谱T2反映了岩石孔径分布情况,而当孔径小到某一程度后,孔隙中的流体将被毛管力所束缚,而无法流动。因此,在T2弛豫时间谱上存在一个界限,当孔隙流体的弛豫时间大于某一弛豫时间时,流体为可动流体,反之为束缚流体。这个弛豫时间界限,称为可动流体截止值。可动流体T2截止值通常需要通过核磁共振和离心法相结合 来确定。一维弛豫时间T2谱上的截止值,在核磁共振二维谱上就是一条截止线。核磁共振二维谱可以根据截止线,扩散系数线,油相关系线等,综合判断、选择可动流体区域,***将自动计算选择区域体积分与总体积分的比值,获得可动流体饱和度。同理,可获得可动水饱和度,束缚水饱和度。可动流体饱和度已经是评价储层的一个重要指标,研究可动流体饱和度对低渗透油田勘探和开发具有重要的意义。
(4)渗透率
核磁共振二维谱反映了地层孔径分布,而地层岩石渗透率又与孔径(孔喉)有一定的关系,因此,可以根据核磁共振二维谱计算出岩样渗透率,实验表明可以利用经验公式11来计算岩样渗透率。
其中,Knmr是核磁共振渗透率,φnmr是核磁共振孔隙度,Swi是束缚水饱和度,C是待定系数,具有地区经验性,需要通过室内核磁共振岩心分析来精确确定。根据经验公式(13)计算的核磁共振渗透率,很好的解决了根据Coates模型计算时,含油饱和度高、含油密度大的岩样的核磁共振渗透率偏小的问题。
(5)原油系数
不同性质的原油,单位体积核磁共振信号量不同,这是由于原油的成分复杂,各种成分含量不同引起的。由于密度高的原油含有大量的重组份,特别是稠油,含有较高的胶质沥青质,其碳氢比较大,单位体积核磁共振信号量较小。因此,当岩样中含有粘度较高的原油时,测量岩样的单位体积核磁共振信号偏小,导致岩样孔隙度、渗透率、可动流体饱和度和含油饱和度等参数有不同程度的偏差,必须进行校正。从而引入了原油系数η,即单位体积恒定梯度场核磁共振标准样 的信号量与单位体积原油的核磁共振信号量的比值:
式中:η为原油系数,A标准样为标准样单位体积核磁共振信号量,φ标准样为标准样孔隙度,A原油为被测原油单位体积核磁共振信号量。
引入原油系数后,由于原油的信号和总信号都发生变化,所以可动流体饱和度、束缚水饱和度、可动水饱和度等参数需要根据解释模板,重新选择区域,进行计算。孔隙度、渗透率、含油饱和度等参数需要对公式进行修正。
原油系数对含油饱和度的修正公式:
原油系数对孔隙度的修正公式:
原油系数对渗透率的修正公式:
其中,φ修正后是修正后的核磁共振孔隙度,K修正后是修正后的核磁共振渗透率,η是原油系数,S油修正后是修正后的含油饱和度,A油是含油信号量,A水是含水信号量,Swi修正后是修正后的束缚水饱和度。
F步骤:获得岩样轴向饱和度分布
恒定梯度场核磁共振岩样分析仪切片厚度是0.3cm,根据步骤A~F,可以对岩样进行切片扫描,获得岩样的局部油、水饱和度;同时重复步骤A~F,可以对岩样进行连续的切片扫描,获得岩样轴向油、水饱和度分布,可动流体饱和度分布等,可以更好的实时的评价储层和识别流体,特别是低孔、低渗以及复杂储集空间的油气藏评价。
本发明的有益效果:本发明提供了一种恒定梯度场核磁共振岩样分析仪。该仪器采用了梯度场的设计和二维的反演,不仅可以获得一维弛豫时间T2谱,提供孔隙结构、孔隙度、渗透率、含油饱和度、可动流体饱和度等岩石物性参数,而且恒定梯度场核磁共振岩样分析仪可以进行切片扫描,提供油、气、水扩散系数、原油粘度、孔隙度分布,油、水饱和度分布等,可以更好的实时的评价储层和识别流体,特别是低孔、低渗以及复杂储集空间的油气藏评价。
恒定梯度场核磁共振岩样分析仪可以放置高温高压探头,模拟核磁共振测井仪井下作业,为实验室内研究地层状态下流体的渗流机理研究提供了新的技术手段,为认识和解决核磁共振测井仪研制过程中的技术瓶颈提供了新方法。
附图说明
图1是一方块图,表示一种恒定梯度场核磁共振岩样分析仪
图2恒定梯度场核磁共振岩样分析仪测量装置的结构示意图,主要包括梯度磁体和探头。
图3是恒定梯度场核磁共振岩样分析仪测量岩样物性参数、识别流体类型和获得油、水饱和度分布等方法的流程图。
图4是CGMF-CPMG脉冲序列,用来进行恒定梯度场下的核磁共振测量。
图5是核磁共振二维谱定标线,单位体积核磁共振信号与孔隙度 之间关系图
图6利用CGMF-CPMG脉冲序列测量的一组正十六烷和蒸馏水的多指数衰减曲线。
图7是利用核磁共振二维反演得到的正十六烷和水的核磁共振二维谱三维立体图
图8是利用核磁共振二维反演得到的正十六烷和水的核磁共振二维谱等高线图
图9a是岩心中饱和油和水的核磁共振二维谱D-T2,
图9b是核磁共振二维谱投影图-弛豫时间T2谱,
图9c是核磁共振二维谱投影图-扩散系数D谱
图10是对岩样进行多次切片扫描,获得的岩心内部轴向含油饱和度分布
图11是对岩样进行多次切片扫描,岩心内部轴向含水饱和度分布
具体实施方式
恒定梯度场核磁共振岩样分析仪由梯度磁体2、探头3、控制***4以及辅助设备计算机5构成。梯度磁体2是一个独立的封闭的腔体。探头2在梯度磁体腔体的中心位置,与梯度磁体构成测量装置1。探头2与控制***3采用电缆连接。控制***3与计算机采用USB口连接。
恒定梯度场核磁共振岩样分析仪的梯度磁体2由轭板14、磁钢15、极板16、侧轭板17构成。其特征是磁钢15与轭板14连接,极板16与磁钢15连接,磁体极板16采用截面为梯形的软铁,其中轭板、极板、侧轭板采用电工纯铁DT4C材料加工制成,磁钢采用2:17SmCo磁性材料制成。磁钢15与轭板14,极板16与磁钢采用BJ-39胶粘接。 梯度磁体2会产生纵向均匀而横向存在梯度的磁场B0形成样品探测区,磁场方向沿着z轴,梯度方向沿x轴。两极板的角度在10~30°之间,其两极板间形成的气隙在150~364mm之间。(详见专利200910092839.1公布的一种用于恒定梯度场岩样分析仪的梯度永磁体)。梯度磁体2外表面是一组温控***6,会对梯度磁体进行10~60度恒温控制,减少温度对仪器测量结果的影响。
恒定梯度场核磁共振岩样分析仪探头3由一个不导电的圆形样品管18,回路线圈19构成。回路线圈19缠绕在样品管18上。回路线圈19包括螺旋管线圈、鞍形线圈等任何其它适用的线圈。探头可作为发射机将激励信号以射频的方式发射到待测样本中和作为接收机接收核磁共振信息。探头3内部是样品探测区20。特别是恒定梯度场核磁共振岩样分析仪可以放置核磁共振高温高压探头。
恒定梯度场核磁共振岩样分析仪控制***4由主控制器7,频率发生器8,功率放大器9,T/R匹配电路10,前置放大器11和接收器12构成。主控制器7主要有3个功能,产生并发射用以激发核磁共振信号的脉冲序列,接收核磁共振回波信号以及将该信号计算处理上传至计算机,主控制器7通过USB接口与计算机5进行互联及数据传输。频率发生器8与主控制器7进行采用电缆互联,该频率发生器将主控制器产生的核磁共振激励信号进行电流驱动,所产生的信号送至功率放大器9进行功率放大。功率放大器9与频率发生器8采用电缆互联,该功率放大器可将信号放大到几百瓦,用来激发待测样品的核磁共振信息。探头3以及发射器/接收器(T/R)匹配电路10通常包括一个谐振电容器、一个T/R转换开关和阻抗匹配电路,与功率放大器9和前置放大器11采用电缆连接;为获得更准确的回波信号数据,要求在射频激励信号发射完之后尽量短的时间内接收回波信号,通过T/R发射 接收转换开关可以快速有效的完成这一转换过程;阻抗匹配电路用来对射频电路进行阻抗匹配。经过前置放大器将核磁共振回波信号放大约70DB,其幅值经放大后可以接收器12采集接收,然后由一个带有缓存器的A/D转换器13输出到主控制器7上,以提供用于进一步使用和分析的所需输出数据。前置放大器11与接收器12,接收器12与A/D转换器13之间采用电缆连接。
同时,该***还包括辅助设备计算机5,对控制***进行控制,从而来控制整套仪器。计算机5与控制***4中的主控制器7之间连接,采用USB口的方式。
恒定梯度场核磁共振岩样分析仪的测量方法分为:A步骤:恒定梯度场脉冲序列及参数设计;B步骤:获得核磁共振二维谱;C步骤:恒定梯度场核磁共振定标;D步骤:识别流体类型;E步骤:计算岩石物性参数;F步骤:获得岩样轴向饱和度分布。如图3所示,是恒定梯度场核磁共振岩样分析仪测量岩样物性参数、识别流体类型和获得油、水饱和度分布等方法的流程图。
A步骤:恒定梯度场脉冲序列及参数设计
脉冲序列是核磁共振技术的灵魂,决定了核磁共振的应用领域。本发明涉及一个CGMF-CPMG脉冲序列。90度脉冲和180度脉冲之间的时间间隔为脉冲间隔tau,单位:us,Nk为tau的个数,单位:个。当k=1时,脉冲间隔tau的值为 ,一个90度脉冲后面跟着P1个180度脉冲,此时CGMF-CPMG脉冲序列就简化成CPMG脉冲序列;当k=2时,一个90度脉冲后面跟着P2个180度脉冲……一直到等于k时,一个90度脉冲后面跟着Pk个180度脉冲。表1是正十六烷和水的一套CGMF-CPMG的脉冲参数,其中,k=4。分别有4个不同的tau值和4个Ne的值,磁场梯度是恒定值,约为20gauss/cm,RD都是12s。由于所 采用的等待时间RD足够长,氢原子核全被充分极化,就不存在对T1的依赖关系。所以,在这种条件下获得的是扩散系数D和弛豫时间T2的信息。该技术的原理是由于分子的布朗运动造成核磁共振信号强度的衰减,而这种衰减与流体分子的扩散系数关联,通过在恒定梯度场下改变CGMF-CPMG脉冲序列中tau的值,来产生扩散效应。90度脉冲的作用是使宏观磁化矢量选择到横向平面上,连续的180度脉冲的作用是重新聚焦除了扩散衰减导致的相散以外的信号,获得自旋回波衰减T2的其它机制的信息。
表1正十六烷和水的一套CGMF-CPMG的脉冲参数
k | tau(us) | Ne(个) | G(gauss/cm) | RD(s) |
1 | 150 | 30720 | 20 | 12 |
2 | 900 | 1536 | 20 | 12 |
3 | 1800 | 512 | 20 | 12 |
4 | 3000 | 256 | 20 | 12 |
B步骤:获得核磁共振二维谱;
核磁共振技术在石油工业中的应用,主要是利用核磁共振仪来测量岩石孔隙中饱和流体氢核的信号。这些应用不仅与流体本身有关,同时还受到岩石的孔隙结构的影响。流体中氢核的弛豫时间的影响因素包括体弛豫,表面弛豫和扩散弛豫。体弛豫T2B是自旋核与晶格之间能量转移的过程,弛豫时间受自旋-自旋偶合作用和共振频率的影响。表面弛豫T2S是多孔介质中孔隙流体分子间的相互作用,弛豫时间和表面弛豫率以及孔隙的比表面有关。扩散弛豫和分子的扩散系数D,磁场梯度G,脉冲间隔tau有关。因此,在岩石孔隙中T2弛豫时间的表达式为:
式中:T2app为表观弛豫时间,T2B是体弛豫时间,T2S是表面弛豫时间,T2D是扩散弛豫时间。
对T2S有:
对T2D有:
其中:D是流体的扩散系数,γ是旋磁比常数,G是磁场梯度,tau是脉冲时间间隔。
根据步骤A中的恒定梯度场的脉冲序列测量的信号b(t,tau)的一般表达式为:
b(t,tau)=∫∫f(D,T2)k2(t,TE,D,T2)dDdT2+ε (4)
离散形式为:
如图6所示,为测量正十六烷和水利用表1中的脉冲序列参数测量得到的一套多指数衰减曲线,该曲线遵循公式5的规律。
弛豫时间和扩散系数在对数坐标系轴上均匀分布,方程(5)就简化成:
bik=fljEik,lj+εik (6)
其中:i=1,……nk,k=1,……q,l=1,……p,j=1,……m
i表示第k个回波串的第i个回波,无量纲;
k表示第k个回波串,无量纲;
l表示提前选择的第l个扩散系数,无量纲;
j表示提前选择的第j个弛豫时间,无量纲;
nk是第k个回波串的回波的个数,无量纲;
q是不同tauk的回波串的个数,无量纲;
p表示提前选择的扩散系数的个数,单位:个;
m表示提前选择的弛豫时间的个数,单位:个;
bk表示脉冲间隔为tauk的第k个回波串的第i个回波的振幅,无量纲;
flj表示在扩散系数为Dl和弛豫时间为T2j时的振幅,无量纲;
γ是旋磁比,单位:MHz/T;
G是磁场梯度,单位:Gauss/cm;
tauk是第k个回波串的脉冲间隔,单位:us。
方程(6)的矩阵形式为:
在求解出的核磁共振二维谱的幅度中很可能存在小于零的分量,这在物理意义上来说是不合逻辑的。在解方程7的时候,还必须加上非负约束,即flj≥0。方程7中的矩阵E高度奇异,条件数高,并且实际测的信号伴随着噪声干扰,给反演带来非常大的难题。
本发明中是利用改进的奇异值分解法进行反演,获得核磁共振二维谱。奇异值分解法给(SVD)出的解是||Ef-b||2最小意义下的一个最优 解,但不一定满足非负约束条件。传统的SVD算法通过缩减E迭代求解,等于降低了解的维数,丢掉了振荡最厉害的坐标分量。正是由于缺失了这部分坐标分量,才导致谱图的不连续。为了获得连续的谱图,改进的奇异值分解法将求解Δf的最小二乘解。假设已知一个初始解f0,令b0=Ef0,这样原来的方程可以变为:E(f-f0)=b-b0,即EΔf=Δb。若求得||EΔf-Δb||2最小意义下的最优解Δf,则f0+Δf是||Ef-b||2最小意义下的最优解f。由于求解的是Δf,所以在实现非负约束时,就不再缩小矩阵E,而只须将f小于零的分量改为零,再重新迭代计算Δf即可,直到f所有坐标分量满足非负约束。由于E在迭代求解的过程中不发生变化,所以求||EΔf-Δb||2最小意义下的最优解无须再进行奇异值分解。与原算法相比较,新算法以计算Δf和Δb代替了矩阵E的奇异值分解过程,不必每次循环都要对矩阵E进行奇异值分解,这样就大大减少了计算量和计算时间。
图7是对图6的测量数据进行核磁共振二维反演得到了正十六烷和水的核磁共振二维谱D-T2三维立体图,图8是与图7对应的核磁共振二维谱等高线图。由于常温下,水的扩散系数约为2.5*10-5cm2/s,所以,从图7和图8中的核磁共振二维谱上可以判断,上边的谱峰是水的信号,下边的谱峰是正十六烷的信号。水的扩散系数
C步骤:恒定梯度场核磁共振定标
定标是核磁共振技术测量岩样孔隙度不可缺少的步骤。恒定梯度场核磁共振岩样分析仪定标,是根据步骤A和B利用恒定梯度场核磁共振岩样分析仪对恒定梯度场的标准样进行测量、反演得到标准样的核磁共振二维谱D-T2,获得单位体积信号与孔隙度之间的关系线。恒定梯度场核磁共振标准样一共有12个,孔隙度分别为0.5%,1%,2%,3%,6%,9%,12%,18%,18%,21%,24%,27%,每次定标至少选择 其中的5个进行测量。图5是选择了其中的9个进行定标,获得的核磁共振二维谱定标线。主要步骤如下:
(1)把孔隙度为27%的标准样放入恒定梯度场核磁共振岩样分析仪探头中利用提前设计的脉冲序列参数进行测量,得到一套多指数衰减数据;
(2)输入标准样为27%的标准样的体积和孔隙度,恒定梯度场核磁共振二维谱反演软件将对样品的测量数据进行反演以及归一化处理,得到样品的单位体积信号量;
(3)重复步骤(1)和(2),分别测量孔隙度为24%,21%,15%,12%,9%,6%,3%,2%的标准样;
(4)软件***将自动对测量的9个标准样进行线性拟合,获得核磁共振二维谱定标曲线,线性度为0.9985。
y=0.3703x-0.0703
D步骤:识别流体类型;
根据步骤A和B,获得实际岩样的核磁共振二维谱D-T2进行流体类型识别。核磁共振二维谱一维是扩散系数D,另一维是弛豫时间T2,所以通过核磁共振二维谱可同时获得被测样品的扩散系数和弛豫时间。由于体相流体,例如原油、天然气和水,不仅扩散系数不同,而且弛豫时间也不同,所以可以利用核磁共振二维谱快速识别流体类型。水的扩散系数与温度有关,常温条件下,水的扩散系数约为2.5*10-5cm2/s。由于原油成分复杂,扩散系数分布广,约为10-7~10-5cm2/s。原油存在体弛豫和扩散弛豫,体弛豫时间T2B与原油粘度成反比,即T2B∝η-1;扩散系数和原油粘度也成反比,即D∝T/300η。因此,原油体弛豫时间T2B和扩散系数D之间存在线性关系,即Do=αT2,其中α≈1.045*10-5cm2/s2。气体的扩散比油或者水的扩散快的多,且与 温度、压力以及气体的密度有关,分布范围约为10-4~10-3cm2/s。如图7和图8所示,在核磁共振二维谱上水峰和油峰互相独立,根据油和水扩散系数的特性可以清晰的分辨油和水。
恒定梯度场核磁共振岩样分析仪软件通过油气水扩散系数和弛豫时间的性质,建立了核磁共振二维谱识别流体类型的模板。图9a是岩心中饱和油和水的核磁共振二维谱D-T2,图9b是图9a核磁共振二维谱投影图-弛豫时间T2谱,图9c是图9a核磁共振二维谱投影图-扩散系数D谱。图9a中两条水平线,上边的是气体的扩散系数线,下边是水的扩散系数线,右边的斜线是油的D-T2关系线。图9a中处在D-T2关系线上的谱峰是油的信号,其它的是水的信号。因为岩样中完全饱和了油和水,而且岩样是强亲水的,所以油相是以体弛豫和扩散弛豫为主,而表面弛豫为辅。因此,油相受孔隙表面作用较小,弛豫时间和扩散系数与体相状态接近。由于岩样是强亲水的,所以水相受到体弛豫、扩散弛豫和表面弛豫三种弛豫共同作用。因此,水相受孔隙表面作用力较大,弛豫时间和扩散系数都受到不同程度的影响。从而导致在一维的弛豫时间T2谱和一维的扩散系数谱D上,部分油和水的信号会互相重叠,而无法直接区分,如图9b和图9c所示。核磁共振二维谱D-T2在一维弛豫时间T2谱的基础上增加了一维扩散系数,从立体空间上可以更好的识别孔隙中流体类型。
E步骤:计算岩石物性参数
根据步骤D获得的实际岩样的核磁共振二维谱D-T2计算岩石物性参数。恒定梯度场核磁共振岩样分析仪可用于测量流体弛豫和扩散的性质,特别是岩石孔隙中的流体性质。由于这些性质对于油气水是不同的,所以恒定梯度场核磁共振岩样分析仪能够提供油、气和水的比例,以及岩石孔隙度、渗透率、含油饱和度,可动流体饱和度、原油 系数等参数,下面以图9a中所示岩样为例,分别计算上述参数,该岩样常规孔隙度为14.28%,渗透率为50.8mD,含油饱和度为33.98%。
(1)孔隙度
核磁共振技术直接对被测样品中的氢核(1H)进行探测,可以准确测量岩样孔隙内流体的量,当岩样孔隙内充满流体时,流体量与孔隙体积相等。因此,采用核磁共振技术可以准确测量岩样孔隙度。利用核磁共振二维谱技术测量岩样孔隙度时,首先,必须保证岩样充分饱和;其次,利用恒定梯度场核磁共振标准样进行定标,获得核磁共振二维谱定标线,即单位体积核磁共振信号与孔隙度之间的关系线;最后,将测量饱和的岩样中流体的核磁共振信号利用定标线进行刻度,即可获得岩样的核磁共振孔隙度。根据公式17得,图9a中所测岩样孔隙度x=2.7*(y+0.0703)=2.7*(5.2186+0.0703)=14.28。
(2)可动流体饱和度
弛豫时间谱T2反映了岩石孔径分布情况,而当孔径小到某一程度后,孔隙中的流体将被毛管力所束缚,而无法流动。因此,在T2弛豫时间谱上存在一个界限,当孔隙流体的弛豫时间大于某一弛豫时间时,流体为可动流体,反之为束缚流体。这个弛豫时间界限,称为可动流体截止值。可动流体T2截止值通常需要通过核磁共振和离心法相结合来确定。一维弛豫时间T2谱上的截止值,在核磁共振二维谱上就是一条截止线。图9a中的竖线就是该岩样的核磁共振二维谱可动流体截止线。根据核磁共振二维谱解释模板,即可动流体截止线,扩散系数线,油相关系线等,综合判断、选择可动流体截止线右侧部分为可动流体区域,***将自动计算选择区域体积分与总体积分的比值,可动流体饱和度为80.86%。束缚流体饱和度为100%-80.86%=19.14%。可动流体饱和度已经是评价储层的一个重要指标,研究可动流体饱和度对低渗 透油田勘探和开发具有重要的意义。
(3)含油饱和度
含油饱和度是评价油田储层好坏的重要参数之一,测量含油饱和度的关键就是区分油和水。核磁共振二维谱技术可以测量油和水的弛豫时间和扩散系数,根据核磁共振二维谱解释模板快速识别流体类型,区分油和水。恒定梯度场核磁共振岩样分析仪软件利用点击选择式计算含油饱和度。例如,在图9a中利用鼠标选择处在核磁共振二维谱中D-T2关系线上的谱峰为油的区域,电脑将自动计算选择区域体积分与总体积分的比值,含油饱和度为,31.19%。由于岩样中只饱和了油和水,所以,含水饱和度为68.81%。
(4)渗透率
核磁共振二维谱反映了地层孔径分布,而地层岩石渗透率又与孔径(孔喉)有一定的关系,因此,可以根据核磁共振二维谱计算出岩样渗透率,实验表明可以利用经验公式11来计算岩样渗透率。
其中,Knmr是核磁共振渗透率,φnmr是核磁共振孔隙度,Swi是束缚水饱和度,C是待定系数,具有地区经验性,需要通过室内核磁共振岩心分析来精确确定,国内主要分布在3~12。根据经验公式11计算的核磁共振渗透率,很好的解决了根据Coates模型计算时,含油饱和度高、含油密度大的岩样的核磁共振渗透率偏小的问题。
(5)原油系数
不同性质的原油,单位体积核磁共振信号量不同,这是由于原油的成分复杂,各种成分含量不同引起的。由于密度高的原油含有大量的重组份,特别是稠油,含有较高的胶质沥青质,其碳氢比较大,单位体积核磁共振信号量较小。因此,当岩样中含有粘度较高的原油时, 测量岩样的单位体积核磁共振信号偏小,导致岩样孔隙度、渗透率、可动流体饱和度和含油饱和度等参数有不同程度的偏差,必须进行校正。例如:孔隙度等于24%的标准样,单位体积核磁共振信号量为1258,某原油单位体积核磁共振信号量为4680,则公式12得,原油系数 无量纲。输入原油系数后,恒定梯度场核磁共振岩样分析仪将根据公式15~17自动修正岩样的孔隙度、渗透率、含油饱和度、可动流体饱和度等参数。
F步骤:获得岩样轴向饱和度分布
恒定梯度场核磁共振岩样分析仪切片厚度是0.3cm,根据步骤A~F,可以对岩样进行切片扫描,获得岩样的局部油、水饱和度;同时重复步骤A~F,可以对岩样进行连续的切片扫描,获得岩样轴向油、水饱和度分布,可动流体饱和度分布等,可以更好的实时的评价储层和识别流体,特别是低孔、低渗以及复杂储集空间的油气藏评价。
实验选取的岩心孔隙度为23.82%,渗透率为12.44mD。利用驱替实验使岩心达到束缚水饱和油状态,最终出水量为2.89ml,经计算含油饱和度为77.62%。如图10、图11所示,分别是在饱和油束缚水状态下利用核磁共振二维谱技术测量岩心轴向含油、水饱和度分布。入口端含油饱和度高达88.42%,高出常规含油饱和度10.8%,入口端向出口端含油饱和度逐渐减小,其中从入口端到2.4cm处,含油饱和度逐渐较小,在2.4cm到出口端,含油饱和度降低幅度大幅增加,特别是出口端3.0cm处,含油饱和度降低至48.63%,低于常规含油饱和度28.99%,这是由于出口端的末端效应,导致含油饱和度较低,无法达到束缚水状态。总体而言核磁共振二维谱技术可以测量岩心内部的轴向含油饱和度分布,整块岩心的含油饱和度平均值为75.41%,与常规含油饱和度绝对误差是2.21%。岩心的轴向含水饱和度正好与含油饱 和度相反,如图11所示。从入口端到出口端,岩心驱替后的含油饱和度总的趋势是降低的,含水饱和度是升高的,但局部含油、水饱和度高低不同,主要是由于驱替过程中的的毛管力作用、贾敏效应、指进作用等引起的,导致岩心中油和水分布不均匀。
尽管已经结合优选实例介绍了本发明,但是其目的并不是将本发明限制于所提出的具体形式,而是应当使本发明覆盖在本发明的构思和范围内能够合理做出的各种改进、替换以及等价物,正如在所附的权利要求书所述的那样。
Claims (6)
1.一种恒定梯度场核磁共振岩样分析方法,其特征在于
1)在由磁体产生的恒定梯度场下,根据被测样品中饱和流体设计的CGMF-CPMG脉冲序列中回波串k,脉冲间隔tau,180度脉冲个数Ne,等待时间RD参数,进行核磁共振测量获得测量数据;k值为2~8,tau的值分布在150us~20000us之间,Ne的值分布在128~30720之间,等待时间分布在2000ms~12000ms之间;
2)将步骤1)测量的核磁共振数据转化成核磁共振二维谱D-T2;采用CGMF-CPMG脉冲序列测量的数据遵循多指数衰减规律:
其中:其中:i=1,……nk,k=1,……q,l=1,……p,j=1,……m
i表示第k个回波串的第i个回波,无量纲;
k表示第k个回波串,无量纲;
l表示提前选择的第l个扩散系数,无量纲;
j表示提前选择的第j个弛豫时间,无量纲;
nk是第k个回波串的回波的个数,无量纲;
q是不同tauk的回波串的个数,无量纲;
p表示提前选择的扩散系数的个数,单位:个;
m表示提前选择的弛豫时间的个数,单位:个;
bik表示脉冲间隔为tauk的第k个回波串的第i个回波的振幅,无量纲;
flj表示在扩散系数为Dl和弛豫时间为T2j时的振幅,无量纲;
γ是旋磁比,单位:MHz/T;
G是磁场梯度,单位:Gauss/cm;
tauk是第k个回波串的脉冲间隔,单位:us。
利用改进的奇异值分解法对公式(5)进行反演,获得核磁共振二维谱D-T2。
3)用步骤1)和步骤2)对恒定梯度场的标准样进行测量、反演得到标准样核磁共振二维谱D-T2;取恒定梯度场核磁共振标准样12个,孔隙度分别为0.5%,1%,2%,3%,6%,9%,12%,15%,18%,21%,24%,27%,每次定标至少选择其中的5个进行测量、反演得到的不同孔隙度的标准样的核磁共振二维谱,然后对核磁共振二维谱进行体积分,得到该标准样的核磁共振信号;核磁共振信号与标准样体积的比值就是单位体积核磁共振信号,对单位体积核磁共振信号和孔隙度进行线性拟合,得到它们的关系线:
y=ax+b (8)
其中,y代表单位体积核磁共振信号量,x代表核磁共振孔隙度(%),a代表斜率,b代表截距;
当测量岩样时,测量出单位体积岩样的核磁共振信号,就可以计算出岩样的孔隙度;
4)用步骤1)和步骤2)测量岩样,获得岩样中流体核磁共振二维谱D-T2,根据实际测量岩样的核磁共振二维谱D-T2进行流体类型识别;水的扩散系数是一常数,与温度有关;气的扩散系数与温度和压力有关;原油的的扩散系数和弛豫时间之间存在线性关系;
Dw(T2)=Dw(T) (9)
Dg(T2)=Dg(T,P) (10)
Do=αT2 (11)
恒定梯度场核磁共振岩样分析仪根据油气水扩散系数和弛豫时间的这种核磁共振特性,建立核磁共振二维谱解释模板,根据实际测量 岩样的核磁共振二维谱D-T2在解释模板中的位置来划分油气水,快速识别流体类型;
5)根据步骤4)获得的岩样中流体的核磁共振二维谱D-T2计算岩样的孔隙度、渗透率、含油饱和度、可动流体饱和度、原油系数;
a)通过对岩样核磁共振二维谱进行体积分,根据步骤3)带入恒定梯度场核磁共振定标线,计算岩样核磁共振孔隙度;
b)根据核磁共振孔隙度和束缚水饱和度计算岩样渗透率:
其中,Knmr是核磁共振渗透率,φnmr是核磁共振孔隙度,Swi是束缚水饱和度,C是待定系数;
c)根据步骤4)建立的核磁共振二维谱解释模板快速识别流体类型,区分油和水;利用鼠标选择核磁共振二维谱中油的区域,***将自动计算选择区域体积分与核磁共振二维谱总体积分的比值,获得含油饱和度;以此类推,可以分别得到含油、气、水的饱和度;
d)根据步骤4)建立的核磁共振二维谱解释模板中的可动流体T2截止线,扩散系数线,油相关系线等,综合判断、选择可动流体区域,***将自动计算选择区域体积分与总体积分的比值,获得可动流体饱和度;
e)当岩样中含有粘度较高的原油时,测量岩样的单位体积核磁共振信号偏小,导致岩样孔隙度、渗透率、可动流体饱和度和含油饱和度等参数有不同程度的偏差,必须进行校正;原油系数为单位体积恒定梯度场核磁共振标准样的信号与单位体积原油的核磁共振信号的比值:
式中:η为原油系数,A标准样为标准样单位体积核磁共振信号,φ标准样 为标准样孔隙度,A原油为被测原油单位体积核磁共振信号;
6)恒定梯度场核磁共振岩样分析仪切片厚度是0.3cm,根据步骤1)~5)岩样进行单切片扫描,获得岩样的局部油、水饱和度;同时,重复步骤1)~5),对岩样进行连续的切片扫描,可以获得岩样轴向油、水饱和度分布,可动流体饱和度分布,能更好的实时的评价储层和识别流体。
3.一种权利要求1所述的恒定梯度场核磁共振岩样分析仪,其特征在于:分析仪由测量装置(1)和控制***(4)构成;测量装置由梯度磁体(2)、探头(3)和温控***构成;梯度磁体是一个封闭的腔体,探头位于梯度磁体腔体的中心位置,探头与控制***通过电缆与控制***的T/R匹配电路(10)连接;控制***通过USB口与计算 机(5)连接。
4.根据权利要求3所述的恒定梯度场核磁共振岩样分析仪,其特征在于:恒定梯度场核磁共振岩样分析仪的梯度磁体(2)由轭板(14)、磁钢(15)、极板(16)、侧轭板(17)构成;轭板和侧轭板组成梯度磁体腔体,磁钢与轭板连接,极板截面为梯形结构,与磁钢连接;轭板、极板、侧轭板采用电工纯铁DT4C材料加工制成,磁钢采用2:17SmCo磁性材料制成;磁钢与轭板,极板与磁钢采用BJ-39胶粘接;梯度磁体会产生纵向均匀而横向存在梯度的磁场B0形成样品探测区,磁场方向沿着z轴,梯度方向沿x轴;两极板的角度在10~30°之间,其两极板间形成的气隙在150~364mm之间;梯度磁体外表面是一组温控***,会对梯度磁体进行10~60度恒温控制,减少温度对仪器测量结果的影响。
5.根据权利要求3所述的恒定梯度场核磁共振岩样分析仪,其特征在于:恒定梯度场核磁共振岩样分析仪探头(3)由一个不导电的圆形样品管(18),回路线圈(19)构成;回路线圈缠绕在样品管上;回路线圈为螺旋管线圈、鞍形线圈或其它适用的线圈;探头能作为发射机将激励信号以射频的方式发射到待测样品中和作为接收机接收核磁共振信息;探头内部是样品探测区(20)。
6.根据权利要求3所述的恒定梯度场核磁共振岩样分析仪,其特征在于:恒定梯度场核磁共振岩样分析仪控制***(4)由主控制器(7),频率发生器(8),功率放大器(9),T/R匹配电路(10),前置放大器(11)、接收器(12)和带有缓存A/D转换器(13)构成;主控制器产生并发射用以激发核磁共振信号的脉冲序列,接收核磁共振回波信号和将该信号计算处理上传至计算机,主控制器通过USB接口与计算机(5)进行互联及数据传输;频率发生器与主控制器采用电缆互 联,频率发生器将主控制器产生的核磁共振激励信号进行电流驱动,所产生的信号送至功率放大器进行功率放大;功率放大器与频率发生器采用电缆互联,功率放大器将信号放大到几百瓦,用来激发待测样品的核磁共振信息;探头和发射器/接收器T/R匹配电路包括一个谐振电容器、一个T/R转换开关和阻抗匹配电路,与功率放大器和前置放大器采用电缆连接;通过T/R发射接收转换开关完成发射接收转换过程;阻抗匹配电路用来对射频电路进行阻抗匹配;经过前置放大器将核磁共振回波信号放大70DB,其幅值经放大后接收器采集接收,然后由带有缓存器的A/D转换器输出到主控制器上,提供用于下一步使用和分析的所需输出数据;前置放大器与接收器,接收器与A/D转换器之间采用电缆连接。
Priority Applications (4)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201110272477.1A CN102998322B (zh) | 2011-09-14 | 2011-09-14 | 恒定梯度场核磁共振岩样分析方法及仪器 |
EP11872407.9A EP2765409B1 (en) | 2011-09-14 | 2011-09-28 | Nuclear magnetic resonance rock sample analysis method and instrument with constant gradient field |
PCT/CN2011/001641 WO2013037093A1 (zh) | 2011-09-14 | 2011-09-28 | 恒定梯度场核磁共振岩样分析方法及仪器 |
US14/344,285 US9625601B2 (en) | 2011-09-14 | 2011-09-28 | Nuclear magnetic resonance rock sample analysis method and instrument with constant gradient field |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201110272477.1A CN102998322B (zh) | 2011-09-14 | 2011-09-14 | 恒定梯度场核磁共振岩样分析方法及仪器 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102998322A true CN102998322A (zh) | 2013-03-27 |
CN102998322B CN102998322B (zh) | 2014-08-06 |
Family
ID=47882516
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201110272477.1A Active CN102998322B (zh) | 2011-09-14 | 2011-09-14 | 恒定梯度场核磁共振岩样分析方法及仪器 |
Country Status (4)
Country | Link |
---|---|
US (1) | US9625601B2 (zh) |
EP (1) | EP2765409B1 (zh) |
CN (1) | CN102998322B (zh) |
WO (1) | WO2013037093A1 (zh) |
Cited By (26)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104215921A (zh) * | 2013-05-30 | 2014-12-17 | 北京大学 | 一种带有固定梯度的nmr磁体 |
CN104316554A (zh) * | 2014-10-27 | 2015-01-28 | 中国石油天然气股份有限公司 | 储层孔隙水可动性测试方法 |
CN104634804A (zh) * | 2013-11-08 | 2015-05-20 | 中国石油天然气股份有限公司 | 一种利用核磁共振t2谱确定储层相对渗透率的方法 |
CN104819923A (zh) * | 2015-05-17 | 2015-08-05 | 西南石油大学 | 基于核磁共振的低渗透砂岩储层孔隙结构定量反演方法 |
CN104215652B (zh) * | 2014-08-21 | 2016-08-31 | 中国石油天然气股份有限公司 | 确定油气饱和度的方法和装置 |
CN106066494A (zh) * | 2016-05-24 | 2016-11-02 | 中国地质大学(北京) | 一种火成岩核磁共振孔隙度校正方法和t2分布校正方法 |
CN106383365A (zh) * | 2016-10-28 | 2017-02-08 | 中国地质大学(北京) | 一种利用图版校正火成岩核磁共振孔隙度的方法 |
CN107676082A (zh) * | 2016-08-01 | 2018-02-09 | 中国石油化工股份有限公司 | 用测井合成曲线修正岩心含油饱和度的方法 |
CN107807078A (zh) * | 2017-09-13 | 2018-03-16 | 中国石油天然气股份有限公司 | 地层条件下岩心可动水饱和度在线检测方法 |
CN107941838A (zh) * | 2017-11-14 | 2018-04-20 | 西安石油大学 | 一种二氧化碳驱油过程中沥青质沉淀对孔喉分布影响的定量评价方法 |
CN108005633A (zh) * | 2017-11-27 | 2018-05-08 | 中石化石油工程技术服务有限公司 | 致密气藏的二维核磁共振测井观测模式参数确定方法 |
CN108459042A (zh) * | 2018-06-04 | 2018-08-28 | 中国科学院电工研究所无锡分所 | 一种核磁共振岩心孔隙流体成像方法 |
CN108468536A (zh) * | 2018-05-08 | 2018-08-31 | 中国海洋石油集团有限公司 | 基于紫外光谱微观驱油模拟实验的含油饱和度监测方法 |
CN108827853A (zh) * | 2018-05-22 | 2018-11-16 | 西南石油大学 | 基于核磁共振的致密储层岩电测量装置及测量方法 |
CN109025924A (zh) * | 2018-05-08 | 2018-12-18 | 中国海洋石油集团有限公司 | 基于微观岩石薄片的含油饱和度动态监测平台 |
CN109870734A (zh) * | 2019-03-08 | 2019-06-11 | 中国石油天然气集团有限公司 | 三分量重力梯度场、三分量磁力梯度场获取装置和方法 |
CN111351813A (zh) * | 2020-03-17 | 2020-06-30 | 无锡鸣石峻致医疗科技有限公司 | 一种基于非均匀场磁共振***的表观扩散系数测量方法 |
CN111366601A (zh) * | 2020-03-05 | 2020-07-03 | 中国石油天然气股份有限公司 | 基于核磁共振t2谱的乳化液成分相似程度定量评价方法 |
CN111380890A (zh) * | 2018-12-29 | 2020-07-07 | 中国石油天然气股份有限公司 | 核磁共振井下二维谱流体识别探头及装置 |
CN111380891A (zh) * | 2018-12-29 | 2020-07-07 | 中国石油化工股份有限公司 | 构建钻井液地层油判识图版的方法及钻遇油层判识方法 |
CN111379560A (zh) * | 2018-12-29 | 2020-07-07 | 中国石油化工股份有限公司 | 一种构建储层油水性质判识图版的方法 |
CN111595887A (zh) * | 2019-11-29 | 2020-08-28 | 苏州纽迈分析仪器股份有限公司 | 样品含氢流体量一维空间分布的检测方法 |
CN112129801A (zh) * | 2020-09-14 | 2020-12-25 | 北京大学 | 一种多孔介质的选层t2弛豫谱快速测试方法 |
CN113624797A (zh) * | 2020-05-07 | 2021-11-09 | 中国石油天然气股份有限公司 | 致密岩心内的流体检测方法及装置 |
CN113781452A (zh) * | 2021-09-15 | 2021-12-10 | 大庆油田有限责任公司 | 页岩含油饱和度的确定方法及装置、电子设备和存储介质 |
CN113933333A (zh) * | 2021-09-09 | 2022-01-14 | 西安石油大学 | 致密储层co2驱沥青质沉积对岩石润湿性变化特征评价方法 |
Families Citing this family (30)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2015193142A1 (en) * | 2014-06-19 | 2015-12-23 | Total Sa | Nmr method for determining non-oil volume of a rock sample |
US11300531B2 (en) * | 2014-06-25 | 2022-04-12 | Aspect Ai Ltd. | Accurate water cut measurement |
CN104297278B (zh) * | 2014-09-28 | 2016-06-01 | 中国石油大学(华东) | 改良式CPMG序列T2-G采集参数t0自适应方法 |
CN105891248A (zh) * | 2015-04-17 | 2016-08-24 | 北京大学 | 一种高温高压岩石物性及渗流机理核磁共振在线测试装置 |
CN106324009B (zh) * | 2015-07-06 | 2018-05-11 | 中国石油化工股份有限公司 | 岩心驱替实验设备和实验方法 |
MX2018000120A (es) * | 2015-07-31 | 2018-03-22 | Halliburton Energy Services Inc | Aparato y metodo para procesar e interpretar datos de adquisicion de registros de resonancia magnetica nuclear (rmn). |
US10156531B2 (en) * | 2015-10-27 | 2018-12-18 | Saudi Arabian Oil Company | Measuring connectivity between different pore types in porous media |
CN105606517B (zh) * | 2016-03-25 | 2017-05-03 | 中国地质大学(北京) | 一种结合核磁共振测量低渗透储层相对渗透率的仪器 |
CN105954170B (zh) * | 2016-05-04 | 2017-02-22 | 中国石油大学(华东) | 一种考虑背景噪声信号量的核磁孔隙度计算方法 |
US10488352B2 (en) * | 2017-01-27 | 2019-11-26 | Saudi Arabian Oil Company | High spatial resolution nuclear magnetic resonance logging |
US10345251B2 (en) | 2017-02-23 | 2019-07-09 | Aspect Imaging Ltd. | Portable NMR device for detecting an oil concentration in water |
EP3593159B1 (en) * | 2017-03-07 | 2022-06-01 | Saudi Arabian Oil Company | Absolute porosity and pore size determination of pore types in media with varying pore sizes using nmr |
US10338014B2 (en) * | 2017-07-27 | 2019-07-02 | Saudi Arabian Oil Company | Estimating formation properties using saturation profiles |
US10545105B2 (en) | 2017-07-27 | 2020-01-28 | Saudi Arabian Oil Company | Estimating formation properties using saturation profiles |
CN111721684A (zh) * | 2019-03-22 | 2020-09-29 | 中国石油天然气股份有限公司 | 一种砾岩含油饱和度测定装置与方法 |
CN109856176B (zh) * | 2019-03-27 | 2022-03-08 | 中国石油大学(华东) | 基于nmr和ltna的致密储层全尺寸孔隙定量表征方法 |
BR112021020493A2 (pt) * | 2019-04-26 | 2021-12-07 | Hyperfine Inc | Técnicas para o controle dinâmico de um sistema de imagem de ressonância magnética |
US11460602B2 (en) | 2019-09-05 | 2022-10-04 | Saudi Arabian Oil Company | Systems and methods for saturation logging of hydrocarbon wells |
US11112373B1 (en) | 2020-02-26 | 2021-09-07 | Saudi Arabian Oil Company | Systems and methods for slice selective nuclear magnetic resonance testing of fractured core plugs to determine in-situ pore volume |
US11143607B2 (en) * | 2020-03-13 | 2021-10-12 | King Fahd University Of Petroleum And Minerals | Method for evaluation of permeability anisotropy using NMR diffusion measurements for oil and gas wells |
CN112834543B (zh) * | 2020-04-28 | 2024-05-14 | 苏州纽迈分析仪器股份有限公司 | 基于脉冲梯度硬件结构的一维空间选层t2谱测试方法 |
CN112816516A (zh) * | 2021-02-07 | 2021-05-18 | 中国石油天然气集团有限公司 | 一种饱含甲烷岩样高温高压核磁共振t2谱实验室测量方法 |
CN113433157B (zh) * | 2021-06-24 | 2023-12-15 | 西南石油大学 | 基于核磁共振t2谱建立随机单元等效岩心模型的方法 |
US11614417B2 (en) | 2021-07-06 | 2023-03-28 | Saudi Arabian Oil Company | Determining saturation in low resistivity pay zones |
CN113933260B (zh) * | 2021-09-15 | 2024-06-11 | 核工业北京地质研究院 | 一种热液铀矿床流体活动中心的识别方法 |
US11788978B2 (en) | 2021-11-12 | 2023-10-17 | Saudi Arabian Oil Company | Method and system for determining geological models using spatial properties based on nuclear magnetic resonance data |
US11821862B2 (en) | 2021-11-12 | 2023-11-21 | Saudi Arabian Oil Company | Method for measuring the spatial water permeability profile of porous media by using non-destructive nuclear magnetic resonance technique |
CN115788373B (zh) * | 2021-11-26 | 2024-04-30 | 中国石油天然气集团有限公司 | 一种co2与页岩油藏油气置换规律评价方法 |
CN114880755B (zh) * | 2022-07-07 | 2022-09-20 | 西南交通大学 | 铁路隧道的围岩分级方法、装置、设备及可读存储介质 |
CN118243591A (zh) * | 2024-05-29 | 2024-06-25 | 苏州纽迈分析仪器股份有限公司 | 一种非规则形状岩心的孔隙度测试方法及装置 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1193387A (zh) * | 1995-06-26 | 1998-09-16 | 纽马公司 | 采用扩散和弛豫测井测量技术评估地岩的核磁共振***和方法 |
US20030169040A1 (en) * | 2002-03-11 | 2003-09-11 | Schlumberger Technology Corporation | Method and apparatus for NMR measurement of wettability |
US20030214287A1 (en) * | 2002-05-15 | 2003-11-20 | Boqin Sun | Methods of decoupling diffusion effects from relaxation times to determine properties of porous media containing fluids and multi-dimensional representation of those properties |
CN1707290A (zh) * | 2004-06-04 | 2005-12-14 | 施卢默格海外有限公司 | 使用脉冲场梯度nmr测量确定流体取样测井工具中流体性质的方法和装置 |
CN102024546A (zh) * | 2009-09-09 | 2011-04-20 | 中国石油天然气股份有限公司 | 一种用于恒定梯度场岩样分析仪的梯度永磁体 |
Family Cites Families (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4710713A (en) * | 1986-03-11 | 1987-12-01 | Numar Corporation | Nuclear magnetic resonance sensing apparatus and techniques |
US6111408A (en) * | 1997-12-23 | 2000-08-29 | Numar Corporation | Nuclear magnetic resonance sensing apparatus and techniques for downhole measurements |
US7034528B2 (en) | 2003-02-27 | 2006-04-25 | Schlumberger Technology Corporation | Methods for formation evaluation based on multi-dimensional representation of nuclear magnetic resonance data |
US6937014B2 (en) * | 2003-03-24 | 2005-08-30 | Chevron U.S.A. Inc. | Method for obtaining multi-dimensional proton density distributions from a system of nuclear spins |
CN100373172C (zh) | 2004-10-21 | 2008-03-05 | 北京大学 | 用岩石核磁共振弛豫信号测量地层岩石物性的设备及测量方法 |
US7893692B2 (en) | 2008-11-03 | 2011-02-22 | Schlumberger Technology Corporation | Method for estimating the formation productivity from nuclear magnetic resonance measurements |
US8427145B2 (en) * | 2010-03-24 | 2013-04-23 | Schlumberger Technology Corporation | System and method for emulating nuclear magnetic resonance well logging tool diffusion editing measurements on a bench-top nuclear magnetic resonance spectrometer for laboratory-scale rock core analysis |
-
2011
- 2011-09-14 CN CN201110272477.1A patent/CN102998322B/zh active Active
- 2011-09-28 WO PCT/CN2011/001641 patent/WO2013037093A1/zh active Application Filing
- 2011-09-28 US US14/344,285 patent/US9625601B2/en active Active
- 2011-09-28 EP EP11872407.9A patent/EP2765409B1/en active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1193387A (zh) * | 1995-06-26 | 1998-09-16 | 纽马公司 | 采用扩散和弛豫测井测量技术评估地岩的核磁共振***和方法 |
US20030169040A1 (en) * | 2002-03-11 | 2003-09-11 | Schlumberger Technology Corporation | Method and apparatus for NMR measurement of wettability |
US20030214287A1 (en) * | 2002-05-15 | 2003-11-20 | Boqin Sun | Methods of decoupling diffusion effects from relaxation times to determine properties of porous media containing fluids and multi-dimensional representation of those properties |
CN1707290A (zh) * | 2004-06-04 | 2005-12-14 | 施卢默格海外有限公司 | 使用脉冲场梯度nmr测量确定流体取样测井工具中流体性质的方法和装置 |
CN102024546A (zh) * | 2009-09-09 | 2011-04-20 | 中国石油天然气股份有限公司 | 一种用于恒定梯度场岩样分析仪的梯度永磁体 |
Non-Patent Citations (4)
Title |
---|
刘美杰: "江苏油田核磁共振录井综合应用研究", 《中国优秀硕士学位论文全文数据库 基础科学辑》 * |
顾兆斌 等: "2D NMR技术在石油测井中的应用", 《波谱学杂志》 * |
顾兆斌 等: "基于核磁共振二维谱技术识别储层流体类型", 《西南石油大学学报(自然科学版)》 * |
顾兆斌 等: "核磁共振二维谱反演", 《波谱学杂志》 * |
Cited By (44)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104215921A (zh) * | 2013-05-30 | 2014-12-17 | 北京大学 | 一种带有固定梯度的nmr磁体 |
CN104634804B (zh) * | 2013-11-08 | 2016-10-26 | 中国石油天然气股份有限公司 | 一种利用核磁共振t2谱确定储层相对渗透率的方法 |
CN104634804A (zh) * | 2013-11-08 | 2015-05-20 | 中国石油天然气股份有限公司 | 一种利用核磁共振t2谱确定储层相对渗透率的方法 |
CN104215652B (zh) * | 2014-08-21 | 2016-08-31 | 中国石油天然气股份有限公司 | 确定油气饱和度的方法和装置 |
CN104316554A (zh) * | 2014-10-27 | 2015-01-28 | 中国石油天然气股份有限公司 | 储层孔隙水可动性测试方法 |
CN104316554B (zh) * | 2014-10-27 | 2016-07-06 | 中国石油天然气股份有限公司 | 储层孔隙水可动性测试方法 |
CN104819923A (zh) * | 2015-05-17 | 2015-08-05 | 西南石油大学 | 基于核磁共振的低渗透砂岩储层孔隙结构定量反演方法 |
CN104819923B (zh) * | 2015-05-17 | 2017-10-10 | 西南石油大学 | 基于核磁共振的低渗透砂岩储层孔隙结构定量反演方法 |
CN106066494A (zh) * | 2016-05-24 | 2016-11-02 | 中国地质大学(北京) | 一种火成岩核磁共振孔隙度校正方法和t2分布校正方法 |
CN106066494B (zh) * | 2016-05-24 | 2018-05-18 | 中国地质大学(北京) | 一种火成岩核磁共振孔隙度校正方法和t2分布校正方法 |
CN107676082A (zh) * | 2016-08-01 | 2018-02-09 | 中国石油化工股份有限公司 | 用测井合成曲线修正岩心含油饱和度的方法 |
CN107676082B (zh) * | 2016-08-01 | 2020-07-14 | 中国石油化工股份有限公司 | 用测井合成曲线修正岩心含油饱和度的方法 |
CN106383365A (zh) * | 2016-10-28 | 2017-02-08 | 中国地质大学(北京) | 一种利用图版校正火成岩核磁共振孔隙度的方法 |
CN106383365B (zh) * | 2016-10-28 | 2019-03-19 | 中国地质大学(北京) | 一种利用图版校正火成岩核磁共振孔隙度的方法 |
CN107807078B (zh) * | 2017-09-13 | 2019-11-08 | 中国石油天然气股份有限公司 | 地层条件下岩心可动水饱和度在线检测方法 |
CN107807078A (zh) * | 2017-09-13 | 2018-03-16 | 中国石油天然气股份有限公司 | 地层条件下岩心可动水饱和度在线检测方法 |
CN107941838A (zh) * | 2017-11-14 | 2018-04-20 | 西安石油大学 | 一种二氧化碳驱油过程中沥青质沉淀对孔喉分布影响的定量评价方法 |
CN107941838B (zh) * | 2017-11-14 | 2019-08-23 | 西安石油大学 | 一种二氧化碳驱油过程中沥青质沉淀对孔喉分布影响的定量评价方法 |
CN108005633A (zh) * | 2017-11-27 | 2018-05-08 | 中石化石油工程技术服务有限公司 | 致密气藏的二维核磁共振测井观测模式参数确定方法 |
CN108468536B (zh) * | 2018-05-08 | 2023-04-25 | 中国海洋石油集团有限公司 | 基于紫外光谱微观驱油模拟实验的含油饱和度监测方法 |
CN108468536A (zh) * | 2018-05-08 | 2018-08-31 | 中国海洋石油集团有限公司 | 基于紫外光谱微观驱油模拟实验的含油饱和度监测方法 |
CN109025924A (zh) * | 2018-05-08 | 2018-12-18 | 中国海洋石油集团有限公司 | 基于微观岩石薄片的含油饱和度动态监测平台 |
CN109025924B (zh) * | 2018-05-08 | 2024-04-05 | 中国海洋石油集团有限公司 | 基于微观岩石薄片的含油饱和度动态监测平台 |
CN108827853A (zh) * | 2018-05-22 | 2018-11-16 | 西南石油大学 | 基于核磁共振的致密储层岩电测量装置及测量方法 |
CN108827853B (zh) * | 2018-05-22 | 2024-02-06 | 西南石油大学 | 基于核磁共振的致密储层岩电测量装置及测量方法 |
CN108459042A (zh) * | 2018-06-04 | 2018-08-28 | 中国科学院电工研究所无锡分所 | 一种核磁共振岩心孔隙流体成像方法 |
CN111380891A (zh) * | 2018-12-29 | 2020-07-07 | 中国石油化工股份有限公司 | 构建钻井液地层油判识图版的方法及钻遇油层判识方法 |
CN111380890B (zh) * | 2018-12-29 | 2024-03-01 | 中国石油天然气股份有限公司 | 核磁共振井下二维谱流体识别探头及装置 |
CN111379560A (zh) * | 2018-12-29 | 2020-07-07 | 中国石油化工股份有限公司 | 一种构建储层油水性质判识图版的方法 |
CN111380890A (zh) * | 2018-12-29 | 2020-07-07 | 中国石油天然气股份有限公司 | 核磁共振井下二维谱流体识别探头及装置 |
CN111379560B (zh) * | 2018-12-29 | 2023-04-07 | 中国石油化工股份有限公司 | 一种构建储层油水性质判识图版的方法 |
CN109870734B (zh) * | 2019-03-08 | 2024-05-28 | 中国石油天然气集团有限公司 | 三分量重力梯度场、三分量磁力梯度场获取装置和方法 |
CN109870734A (zh) * | 2019-03-08 | 2019-06-11 | 中国石油天然气集团有限公司 | 三分量重力梯度场、三分量磁力梯度场获取装置和方法 |
CN111595887A (zh) * | 2019-11-29 | 2020-08-28 | 苏州纽迈分析仪器股份有限公司 | 样品含氢流体量一维空间分布的检测方法 |
CN111366601A (zh) * | 2020-03-05 | 2020-07-03 | 中国石油天然气股份有限公司 | 基于核磁共振t2谱的乳化液成分相似程度定量评价方法 |
CN111366601B (zh) * | 2020-03-05 | 2023-05-26 | 中国石油天然气股份有限公司 | 基于核磁共振t2谱的乳化液成分相似程度定量评价方法 |
CN111351813A (zh) * | 2020-03-17 | 2020-06-30 | 无锡鸣石峻致医疗科技有限公司 | 一种基于非均匀场磁共振***的表观扩散系数测量方法 |
CN111351813B (zh) * | 2020-03-17 | 2021-09-24 | 无锡鸣石峻致医疗科技有限公司 | 一种基于非均匀场磁共振***的表观扩散系数测量方法 |
CN113624797B (zh) * | 2020-05-07 | 2024-03-26 | 中国石油天然气股份有限公司 | 致密岩心内的流体检测方法及装置 |
CN113624797A (zh) * | 2020-05-07 | 2021-11-09 | 中国石油天然气股份有限公司 | 致密岩心内的流体检测方法及装置 |
CN112129801A (zh) * | 2020-09-14 | 2020-12-25 | 北京大学 | 一种多孔介质的选层t2弛豫谱快速测试方法 |
CN113933333A (zh) * | 2021-09-09 | 2022-01-14 | 西安石油大学 | 致密储层co2驱沥青质沉积对岩石润湿性变化特征评价方法 |
CN113933333B (zh) * | 2021-09-09 | 2024-06-07 | 西安石油大学 | 致密储层co2驱沥青质沉积对岩石润湿性变化特征评价方法 |
CN113781452A (zh) * | 2021-09-15 | 2021-12-10 | 大庆油田有限责任公司 | 页岩含油饱和度的确定方法及装置、电子设备和存储介质 |
Also Published As
Publication number | Publication date |
---|---|
EP2765409A4 (en) | 2015-10-14 |
US9625601B2 (en) | 2017-04-18 |
EP2765409B1 (en) | 2019-11-13 |
CN102998322B (zh) | 2014-08-06 |
EP2765409A1 (en) | 2014-08-13 |
WO2013037093A1 (zh) | 2013-03-21 |
US20140285196A1 (en) | 2014-09-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102998322B (zh) | 恒定梯度场核磁共振岩样分析方法及仪器 | |
CN100373172C (zh) | 用岩石核磁共振弛豫信号测量地层岩石物性的设备及测量方法 | |
CN1132019C (zh) | 采用扩散和弛豫测井测量技术评估地岩的核磁共振***和方法 | |
CN103353462B (zh) | 一种基于核磁共振成像的岩石非均质性定量评价方法 | |
Hfirlimann et al. | Diffusion-editing: New NMR measurement of saturation and pore geometry | |
Yan et al. | A robust NMR method to measure porosity of low porosity rocks | |
CN100504446C (zh) | 通过比较不同勘测深度处的核磁共振响应来探测烃的方法 | |
CN105849509A (zh) | 解释nmr信号以给出气/液***的多相流体流动测量的方法 | |
Newgord et al. | Simultaneous assessment of wettability and water saturation using 2D NMR measurements | |
CN104453874A (zh) | 一种基于核磁共振的砂砾岩储层含油饱和度的计算方法 | |
CN107748126A (zh) | 一种获取岩石孔隙尺寸和孔隙表面弛豫率的核磁共振方法 | |
Deng et al. | Effects and corrections for mobile NMR measurement | |
Tian et al. | Rock-type definition and pore characterization of tight carbonate rocks based on thin sections and MICP and NMR experiments | |
CN103675722A (zh) | 岩石t2-g实验采集参数自动匹配方法 | |
Carneiro et al. | Evaluating pore space connectivity by NMR diffusive coupling | |
Shi et al. | A proposed NMR solution for multi-phase flow fluid detection | |
McPhee et al. | Nuclear magnetic resonance (NMR) | |
Shen et al. | Optimized inside-out magnetic resonance probe for soil moisture measuring in situ | |
Liao et al. | NMR fluid analyzer applying to petroleum industry | |
Garcia et al. | Integrated workflow to estimate permeability through quantification of rock fabric using joint interpretation of nuclear magnetic resonance and electric measurements | |
Wang et al. | Evaluation of the pore structure of reservoirs based on NMR T 2 Spectrum decomposition | |
US11598735B2 (en) | Combined discrete gaussian analysis of MICP and NMR T2 distributions of multi-modal carbonate rocks | |
XIE et al. | The (T2, D) NMR logging method for fluids characterization | |
Sun et al. | NMR isotherm studies of gas shales | |
Xu et al. | A new permeability predictive model based on NMR data for sandstone reservoirs |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |