CN112666048B - 一种利用动态核极化增强测试材料润湿性及孔径大小的方法 - Google Patents

一种利用动态核极化增强测试材料润湿性及孔径大小的方法 Download PDF

Info

Publication number
CN112666048B
CN112666048B CN201910983179.XA CN201910983179A CN112666048B CN 112666048 B CN112666048 B CN 112666048B CN 201910983179 A CN201910983179 A CN 201910983179A CN 112666048 B CN112666048 B CN 112666048B
Authority
CN
China
Prior art keywords
detected
porous material
model
dnp
nmr
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
Application number
CN201910983179.XA
Other languages
English (en)
Other versions
CN112666048A (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.)
Wuhan Institute of Physics and Mathematics of CAS
Original Assignee
Wuhan Institute of Physics and Mathematics of CAS
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 Wuhan Institute of Physics and Mathematics of CAS filed Critical Wuhan Institute of Physics and Mathematics of CAS
Priority to CN201910983179.XA priority Critical patent/CN112666048B/zh
Publication of CN112666048A publication Critical patent/CN112666048A/zh
Application granted granted Critical
Publication of CN112666048B publication Critical patent/CN112666048B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02ATECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/30Assessment of water resources

Landscapes

  • Testing Resistance To Weather, Investigating Materials By Mechanical Methods (AREA)

Abstract

本发明公开了一种利用动态核极化增强测试材料润湿性及孔径大小的方法,其步骤为:1、向浸润液中加入能DNP增强其NMR信号强度的自由基;2、取N个不同的待测多孔材料,将各待测多孔材料分别用浸润液分别浸润;3、将各各待测多孔材料进行NMR实验,计录各待测多孔材料的NMR信号强度,记为S1、S2……SN;4、将各待测多孔材料进行DNP‑NMR实验,计录各待测多孔材料的DNP增强的NMR信号强度,记为Sd1、Sd2……SdN;5、计算各待测多孔材料的DNP增强倍数M1、M2……MN;6、将M1、M2……MN进行比较,分析M1、M2……MN的变化趋势来判断浸润液对于各多孔材料是否具有润湿性以及各待测多孔材料的孔径大小。该方法简单,操作方便,测试时间短,测试效率高,且准确可靠。

Description

一种利用动态核极化增强测试材料润湿性及孔径大小的方法
技术领域
本发明属于磁共振应用技术领域,具体涉及一种利用动态核极化增强测试材料润湿性及孔径大小的方法。
背景技术
储层的润湿性是直接影响储层残余油饱和度、束缚水饱和度、油水两相相对渗透率、毛细管压力、油水分布等参数的关键因素之一。了解储层的润湿性及孔径分布有利于研究储层残余流体相的类型、评估储层含量,从而制定合理的开采方案,提高石油开采率。
NMR技术为测试、分析岩心等多孔材料提供了一种有效的研究方法,依据材料中两相流体饱和度、离心后残余流体相的含量或弛豫分布情况,研究材料的润湿性,从而为分析储层特性提供参考。上述方法主要依赖于离心处理后样品内束缚水饱和度、残余油饱和度等信息。然而,该状态下样品内的流体含量少,NMR信号弱,加之岩心样品常在低场条件下测试,NMR信号强度会进一步减小,为保证测试结果的可靠性,需要长时间的累加实验,样品测试分析周期较长。
润湿相流体与材料表面相互作用,使流体的扩散、弛豫特性不同于自由状态,将流体的弛豫-扩散特性或弛豫-弛豫特性相结合可以获取流体在材料中的润湿性。一般地2DNMR方法不需要对样品进行离心处理,因而样品内流体含量相对较大,但仍受限于低场NMR信噪比的限制,从而使测得的弛豫、扩散分布宽化,T1/T2或T2-D并不能有效地反映润湿相的改变。此外,为确保2D方法的准确性,需要远大于1DNMR方法所需的累加时间。由此,为提高低场条件下NMR方法在分析样品润湿性的可靠性与分析测试效率,需要提高NMR信噪比或者发展新的方法。
发明内容
基于上述现有技术,本发明提供了一种利用动态核极化增强测试材料润湿性及孔径大小的方法,该方法简单,操作方便,测试时间短,测试效率高,且准确可靠。
实现本发明上述目的所采用的技术方案为:
一种利用动态核极化增强测试材料润湿性及孔径大小的方法,包括如下步骤:
1、向浸润液中加入能DNP增强其NMR信号强度的自由基;
2、取N个不同的待测多孔材料,所有的多孔材料的材质均为同一类,将各待测多孔材料分别用溶解有自由基的浸润液分别浸润,浸润完成后擦除各待测多孔材料表面的残余液体;
3、将各待测多孔材料进行NMR实验,采集各待测多孔材料的NMR信号,计录各待测多孔材料的NMR信号强度,记为S1、S2……SN
4、在与步骤1.3同等的实验条件下,施加微波激励,将各待测多孔材料进行DNP-NMR实验,采集各待测多孔材料DNP增强的NMR信号,计录各待测多孔材料的DNP增强的NMR信号强度,记为Sd1、Sd2……SdN
5、计算各待测多孔材料的DNP增强倍数M1、M2……MN,M1=|Sd1/S1|、M2=|Sd2/S2|……MN=|SdN/SN|;
6、将M1、M2……MN进行比较,分析M1、M2……MN的变化趋势,若无明显变化,则浸润液对于该类多孔材料具有非润湿性,若变化趋势明显,则浸润液对于该类多孔材料具有润湿性,而且待测多孔材料的增强倍数越大,则待测多孔材料的孔径越大。
进一步,所述的浸润液为水或油。
进一步,所述的NMR信号强度可用时域信号幅度、信号频域峰值或者谱峰面积进行表征。
与现有技术相比,本发明的有益效果和优点在于:
1、该方法利用润湿相流体与待测多孔材料接触时,浸润液扩散系数及弛豫时间减小,多孔材料孔径越小,两者减小幅度越大,从而使与动态核极化相关的泄露因子、耦合因子减小,相应地多孔材料内浸润液的动态核极化增强随之变小,依据多孔材料的动态核极化增强的变化趋势可直观反映多孔材料的润湿性、孔径大小。
2、该方法可以提高NMR信号灵敏度高,测试时间短,测试效率高,且准确可靠。
3、该方法采用动态核极化增强在材料中的变化趋势表征材料润湿性、孔径大小,可实现批量测试,且数据处理简单,无需对NMR数据进行反演。
具体实施方式
下面结合具体实施例对本发明进行详细说明。
多孔材料模拟模型的构建:
选取粒径分别为77μm、68μm、40μm、27μm、20μm、10μm的玻璃珠,将同一粒径的玻璃珠随机堆叠于10mm样品管中,构建6种多孔材料模拟模型,分别命名为模型A1、模型A2、模型A3、模型A4、模型A5、模型A6;
由式1)计算每一种多孔材料模拟模型中孔隙的平均孔径:
Figure BDA0002235865690000031
式1),
Figure BDA0002235865690000032
为孔隙度,相同粒径构建模型的孔隙度平均值为0.45,d粒径为玻璃珠的粒径;
根据式1)进行计算,模型A1、模型A2、模型A3、模型A4、模型A5、模型A6的孔隙的平均孔径为72.4um、63.9um、37.6um、25.4um、18.8um、9.4um。
以下所有实施例的DNP-NMR实验均在0.06T DNP谱仪上进行,测试序列为单脉冲序列,测试时微波激励施加时间为1s。
实施例1
1、选取水作为浸润液,选择TEMPO(四甲基哌啶氮氧化物)作为自由基,将TEMPO溶解于水中,配制成浓度为10mM的TEMPO水溶液;
2、向模型A1、模型A2、模型A3、模型A4、模型A5、模型A6的样品管中分别注入TEMPO水溶液,对其中堆叠的玻璃珠进行浸润,浸润完成后吸除样品管内上层自由状态的水溶液;
3、将模型A1、模型A2、模型A3、模型A4、模型A5、模型A6进行NMR实验,采集模型A1、模型A2、模型A3、模型A4、模型A5、模型A6的FID信号,计算并记录模型A1、模型A2、模型A3、模型A4、模型A5、模型A6的NMR谱图的谱峰面积;
4、在与步骤3同等的实验条件下,施加微波激励,微波功率为5W,将模型A1、模型A2、模型A3、模型A4、模型A5、模型A6分别进行DNP-NMR实验,采集模型A1、模型A2、模型A3、模型A4、模型A5、模型A6的DNP增强的FID信号,计算并记录模型A1、模型A2、模型A3、模型A4、模型A5、模型A6的DNP-NMR谱图的谱峰面积;
5、计算模型A1、模型A2、模型A3、模型A4、模型A5、模型A6的DNP增强倍数;
6、重复步骤4-5两次,两次施加微波的微波功率分别为10W和20W;
7、步骤3-6测试的结果如表1所示:
Figure BDA0002235865690000033
Figure BDA0002235865690000041
由表1可知,水中1H的DNP增强在6种模型中呈现出明显的减小趋势,表明水与多孔材料表面充分接触,随多孔材料孔径、渗透率的变化,表面弛豫在多孔材料内水相弛豫中所占比重随之而变,从而改变材料内的DNP增强,由此可知,水对于玻璃珠具有润湿性。
此外,在相同微波功率条件下6种模拟模型中1H的DNP增强依次减小,与模型内孔隙的孔径成正相关。
实施例2
1、选取5#矿物油作为浸润液,选择TEMPO(四甲基哌啶氮氧化物)作为自由基,将TEMPO溶解于5#矿物油中,配制成浓度为10mM的TEMPO矿物油溶液;
2、向模型A1、模型A2、模型A3、模型A4、模型A5、模型A6中的样品管分别注入TEMPO矿物油溶液,对其中堆叠的玻璃珠进行浸润,浸润完成后吸除样品管内上层自由状态的矿物油溶液;
3、将模型A1、模型A2、模型A3、模型A4、模型A5、模型A6进行NMR实验,采集模型A1、模型A2、模型A3、模型A4、模型A5、模型A6的FID信号,计算并记录模型A1、模型A2、模型A3、模型A4、模型A5、模型A6的NMR谱图的谱峰面积;
4、在与步骤3同等的实验条件下,施加微波激励,微波功率为5W,将模型A1、模型A2、模型A3、模型A4、模型A5、模型A6分别进行DNP-NMR实验,采集模型A1、模型A2、模型A3、模型A4、模型A5、模型A6的DNP增强的FID信号,计算并记录模型A1、模型A2、模型A3、模型A4、模型A5、模型A6的DNP-NMR谱图的谱峰面积;
5、计算模型A1、模型A2、模型A3、模型A4、模型A5、模型A6的DNP增强倍数;
6、重复步骤4-5两次,两次施加微波的微波功率分别为10W和20W;
7、步骤3-6测试的结果如表2所示:
Figure BDA0002235865690000042
Figure BDA0002235865690000051
由表2可知,油中1H的DNP增强在6种模型中无明显变化,表明随多孔材料孔径或渗透率的变化,油相与多孔材料表面相互作用对油相弛豫影响不大,由此可知,5#矿物油对于玻璃珠没有润湿性。
而且在相同微波功率条件下6种模拟模型中1H的DNP增强无明显变化,表明对于没有润湿性的浸润液浸润的多孔材料,其DNP增强与多孔材料的孔隙的孔径、渗透率等无明显相关性。
由实施例1和实施例2的测试结果可知,玻璃珠表现为水润湿性,且增强越大的模拟模型,对应的孔径越大,与已知条件相一致。
实施例3
1、选取水作为浸润液,选择TEMPO(四甲基哌啶氮氧化物)作为自由基,将TEMPO溶解于水中,配制成浓度为10mM的TEMPO水溶液;
2、选取渗透率分别为950md、100md、50md、3md的砂岩进行测试,四种砂岩按照渗透率从大到小分别命名为砂岩B1、砂岩B2、砂岩B3和砂岩B4,压汞法测得砂岩B1、砂岩B2、砂岩B3、砂岩B的平均孔径依次为22.1um、9.7um、7.5um、0.8um,将砂岩B1、砂岩B2、砂岩B3和砂岩B4分别采用TEMPO水溶液进行浸润,浸润完成后取出各砂岩并擦除各砂岩表面残余的水分;
3、将砂岩B1、砂岩B2、砂岩B3和砂岩B4进行NMR实验,采集砂岩B1、砂岩B2、砂岩B3和砂岩B4的FID信号,计算并记录砂岩B1、砂岩B2、砂岩B3和砂岩B4的NMR谱图的谱峰面积;
4、在与步骤3同等的实验条件下,施加微波激励,微波功率为5W,将砂岩B1、砂岩B2、砂岩B3和砂岩B4分别进行DNP-NMR实验,采集砂岩B1、砂岩B2、砂岩B3和砂岩B4的DNP增强的FID信号,计算并记录砂岩B1、砂岩B2、砂岩B3和砂岩B4的DNP-NMR图谱的谱峰面积;
5、计算砂岩B1、砂岩B2、砂岩B3和砂岩B4的DNP增强倍数;
6、重复步骤4-5两次,两次施加微波的微波功率分别为10W和20W;
7、步骤3-6测试的结果如表3所示:
Figure BDA0002235865690000061
由表3可知,砂岩样品水中1H的DNP增强随砂岩不同而明显变化,表明岩心表面与水之间存在明显的相互作用,水对于砂岩具有润湿性,相应地可以得出四组砂岩的平均孔径大小依次为B1>B2>B3>B4,四种砂岩的这两种特性与已知砂岩的特性相一致,进一步表明采用DNP增强可以快速、批量测试样品的润湿性及孔径大小分布。

Claims (3)

1.一种利用动态核极化增强测试材料润湿性及孔径大小的方法,其特征在于包括如下步骤:
1.1、向浸润液中加入能DNP增强其NMR信号强度的自由基;
1.2、取N个不同的待测多孔材料,所有的多孔材料的材质均为同一类,将各待测多孔材料分别用溶解有自由基的浸润液浸润,浸润完成后擦除各待测多孔材料表面的残余液体;
1.3、将各待测多孔材料进行NMR实验,采集各待测多孔材料的NMR信号,计录各待测多孔材料的NMR信号强度,记为S1、S2……SN
1.4、在与步骤1.3同等的实验条件下,施加微波激励,将各待测多孔材料进行DNP-NMR实验,采集各待测多孔材料DNP增强的NMR信号,计录各待测多孔材料的DNP增强的NMR信号强度,记为Sd1、Sd2……SdN
1.5、计算各待测多孔材料的DNP增强倍数M1、M2……MN,M1=|Sd1/S1|、M2=|Sd2/S2|……MN=|SdN/SN|;
1.6、将M1、M2……MN进行比较,分析M1、M2……MN的变化趋势,若无明显变化,则浸润液对于该类多孔材料具有非润湿性,若变化趋势明显,则浸润液对于该类多孔材料具有润湿性,而且待测多孔材料的增强倍数越大,则待测多孔材料的孔径越大。
2.根据权利要求1所述的利用动态核极化增强测试材料润湿性及孔径大小的方法,其特征在于:所述的浸润液为水或油。
3.根据权利要求1所述的利用动态核极化增强测试材料润湿性及孔径大小的方法,其特征在于:所述的NMR信号强度可用时域信号幅度、信号频域峰值或者谱峰面积进行表征。
CN201910983179.XA 2019-10-16 2019-10-16 一种利用动态核极化增强测试材料润湿性及孔径大小的方法 Active CN112666048B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910983179.XA CN112666048B (zh) 2019-10-16 2019-10-16 一种利用动态核极化增强测试材料润湿性及孔径大小的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910983179.XA CN112666048B (zh) 2019-10-16 2019-10-16 一种利用动态核极化增强测试材料润湿性及孔径大小的方法

Publications (2)

Publication Number Publication Date
CN112666048A CN112666048A (zh) 2021-04-16
CN112666048B true CN112666048B (zh) 2022-06-07

Family

ID=75400315

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910983179.XA Active CN112666048B (zh) 2019-10-16 2019-10-16 一种利用动态核极化增强测试材料润湿性及孔径大小的方法

Country Status (1)

Country Link
CN (1) CN112666048B (zh)

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101316616A (zh) * 2005-12-01 2008-12-03 通用电气医疗集团股份有限公司 应用三苯甲基自由基和顺磁金属离子的动态核极化(dnp)的方法
EP2146215A1 (de) * 2008-07-18 2010-01-20 Bruker BioSpin AG Apparatur zur Durchführung von DNP-NMR Messungen mit Kompensationsanordnung, sowie Verfahren zur Auslegung der Kompensationsanordnung
CN101915716A (zh) * 2010-07-06 2010-12-15 中国石油天然气股份有限公司 一种判断储层岩石的润湿性的方法
WO2016038131A1 (en) * 2014-09-11 2016-03-17 Universite Claude Bernard Lyon I Structured and porous materials particularly suitable as polarizing agents in dissolution dynamic nuclear polarization
CN109030292A (zh) * 2018-09-26 2018-12-18 西南石油大学 一种致密岩石润湿性确定的新方法
CN110082382A (zh) * 2019-05-08 2019-08-02 中国科学院武汉物理与数学研究所 一种利用动态核极化进行油水两相nmr信号分离的方法
CN110082381A (zh) * 2019-05-08 2019-08-02 中国科学院武汉物理与数学研究所 一种利用dnp-mri测量油水分布方法

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101316616A (zh) * 2005-12-01 2008-12-03 通用电气医疗集团股份有限公司 应用三苯甲基自由基和顺磁金属离子的动态核极化(dnp)的方法
EP2146215A1 (de) * 2008-07-18 2010-01-20 Bruker BioSpin AG Apparatur zur Durchführung von DNP-NMR Messungen mit Kompensationsanordnung, sowie Verfahren zur Auslegung der Kompensationsanordnung
CN101915716A (zh) * 2010-07-06 2010-12-15 中国石油天然气股份有限公司 一种判断储层岩石的润湿性的方法
WO2016038131A1 (en) * 2014-09-11 2016-03-17 Universite Claude Bernard Lyon I Structured and porous materials particularly suitable as polarizing agents in dissolution dynamic nuclear polarization
CN109030292A (zh) * 2018-09-26 2018-12-18 西南石油大学 一种致密岩石润湿性确定的新方法
CN110082382A (zh) * 2019-05-08 2019-08-02 中国科学院武汉物理与数学研究所 一种利用动态核极化进行油水两相nmr信号分离的方法
CN110082381A (zh) * 2019-05-08 2019-08-02 中国科学院武汉物理与数学研究所 一种利用dnp-mri测量油水分布方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
Dynamic Nuclear Polarization Enhanced Nuclear Magnetic Resonance and Electron Spin Resonance Studies of Hydration and Local Water Dynamics in Micelle and Vesicle Assemblies;Evan R.McCarney等;《Langmuir》;20081231;第10062-10072页 *
Overhauser dynamic nuclear polarization-enhanced NMR relaxometry;john M.franck等;《Microporous and Mesoporous Materials》;20131231;第113-118页 *

Also Published As

Publication number Publication date
CN112666048A (zh) 2021-04-16

Similar Documents

Publication Publication Date Title
US11378509B2 (en) Evaluation method for hydrogen-bearing components, porosity and pore size distribution of organic-rich shale
CN109030311B (zh) 基于核磁共振t2谱敏感参数的孔隙结构分类与识别方法
US7397240B2 (en) Method of measuring rock wettability by means of nuclear magnetic resonance
Yan et al. A robust NMR method to measure porosity of low porosity rocks
Al-Muthana et al. Wettability as a function of pore size by NMR
WO2013037093A1 (zh) 恒定梯度场核磁共振岩样分析方法及仪器
WO2020224420A1 (zh) 一种利用动态核极化进行油水两相nmr信号分离的方法
WO2020224419A1 (zh) 一种利用dnp-mri测量油水分布方法
WO2015053952A1 (en) Nuclear magnetic resonance (nmr) distributions and pore information
US20240027379A1 (en) Method for quantitative evaluation on sensitivity of shale oil and gas reservoir to injected fluids
Fleury Characterization of shales with low field NMR
Meng et al. Rock fabric of tight sandstone and its influence on irreducible water saturation in Eastern Ordos Basin
Branco et al. NMR study of carbonates wettability
Zhang et al. Influence of pore throat structure and the multiphases fluid seepage on mobility of tight oil reservoir
CN112666048B (zh) 一种利用动态核极化增强测试材料润湿性及孔径大小的方法
CN115078210A (zh) 页岩孔隙结构测试方法
CN110031496A (zh) 一种评价致密储层可动流体分布特征的方法
CN114235641A (zh) 一种核磁共振测量致密岩石润湿性的方法
CN115389387A (zh) 一种评价岩心伤害的实验方法
Zhang et al. Direct determination of surface relaxivity in isolated kerogen by pulsed-field gradient NMR
CN114486977B (zh) 一种不同成因孔隙空间特征的定量化评价方法
Haber-Pohlmeier et al. Relaxation in a Natural Soil: Comparison of Relaxometric Imaging, T 1-T 2 Correlation and Fast-Field Cycling NMR
CN115753864A (zh) 低渗透油藏油水微观分布定量表征方法及装置
CN117007627A (zh) 一种利用核磁共振技术定量表征页岩有机孔孔隙度的分析方法
WO2022204297A1 (en) The combined discrete gaussian analysis of micp and nmr t2 distributions of multi-modal carbonate rocks

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