CN110702722B - 基于核磁共振t2谱计算含水合物多孔介质分形维数的方法 - Google Patents
基于核磁共振t2谱计算含水合物多孔介质分形维数的方法 Download PDFInfo
- Publication number
- CN110702722B CN110702722B CN201911076405.2A CN201911076405A CN110702722B CN 110702722 B CN110702722 B CN 110702722B CN 201911076405 A CN201911076405 A CN 201911076405A CN 110702722 B CN110702722 B CN 110702722B
- Authority
- CN
- China
- Prior art keywords
- porous medium
- pore
- hydrate
- water
- magnetic resonance
- 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
- 238000001228 spectrum Methods 0.000 title claims abstract description 43
- 238000005481 NMR spectroscopy Methods 0.000 title claims abstract description 32
- 238000000034 method Methods 0.000 claims abstract description 51
- 238000004458 analytical method Methods 0.000 claims abstract description 8
- 238000012417 linear regression Methods 0.000 claims abstract description 8
- 238000004364 calculation method Methods 0.000 claims abstract description 6
- 229920006395 saturated elastomer Polymers 0.000 claims abstract description 5
- 239000011148 porous material Substances 0.000 claims description 67
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 53
- 239000012530 fluid Substances 0.000 claims description 14
- CURLTUGMZLYLDI-UHFFFAOYSA-N Carbon dioxide Chemical compound O=C=O CURLTUGMZLYLDI-UHFFFAOYSA-N 0.000 claims description 8
- VYPSYNLAJGMNEJ-UHFFFAOYSA-N Silicium dioxide Chemical compound O=[Si]=O VYPSYNLAJGMNEJ-UHFFFAOYSA-N 0.000 claims description 6
- 238000002347 injection Methods 0.000 claims description 6
- 239000007924 injection Substances 0.000 claims description 6
- 239000006004 Quartz sand Substances 0.000 claims description 4
- 229910002092 carbon dioxide Inorganic materials 0.000 claims description 4
- 239000001569 carbon dioxide Substances 0.000 claims description 4
- 230000001186 cumulative effect Effects 0.000 claims description 4
- 238000006243 chemical reaction Methods 0.000 claims description 2
- 230000015572 biosynthetic process Effects 0.000 description 8
- 150000004677 hydrates Chemical class 0.000 description 4
- 238000000354 decomposition reaction Methods 0.000 description 3
- 230000000052 comparative effect Effects 0.000 description 2
- 239000002245 particle Substances 0.000 description 2
- 239000004576 sand Substances 0.000 description 2
- 230000003595 spectral effect Effects 0.000 description 2
- 238000004611 spectroscopical analysis Methods 0.000 description 2
- 238000003889 chemical engineering Methods 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 239000000243 solution Substances 0.000 description 1
Images
Classifications
-
- 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
-
- 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)
- High Energy & Nuclear Physics (AREA)
- Health & Medical Sciences (AREA)
- Analytical Chemistry (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Environmental & Geological Engineering (AREA)
- Engineering & Computer Science (AREA)
- Chemical & Material Sciences (AREA)
- Geochemistry & Mineralogy (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
本发明属于能源与化工领域,公开了基于核磁共振T2谱计算含水合物多孔介质分形维数的方法。这种低场核磁共振实现水合物生成过程中多孔介质分形维数的计算的方法,首先利用低场核磁共振获得容器中完全水饱和状态下的T2分布谱,然后获得水合物生成过程中的T2分布谱,最后基于T2数据进行分形维数计算。该方法可以较准确地基于线性回归分析直接计算出多孔介质内生成水合物过程中的分形维数,预测多孔介质内水合物动态形状变化。
Description
技术领域
本发明属于能源与化工领域,涉及一种基于核磁共振T2谱计算含水合物多孔介质的分形维数的方法。
背景技术
水合物作为一种清洁型新型能源,主要赋存于多孔介质储层中,因而对含水合物多孔介质以及其中水合物的生成和分解特性进行研究具有重大的现实意义。然而,水合物的生成和分解都会改变多孔介质的孔隙结构,从而影响多孔介质的物性参数,因此用传统的欧氏几何方法无法对多孔介质结构的几何特性进行准确的表征。分形维数的提出提供了一种有效的方法,它是一种可以预测多孔介质内水合物动态形状变化的关键性参数。此外,用低场核磁共振(场强≦0.3T)获取T2分布谱,可以有效地表征多孔介质内部孔隙流体的含量以及分布的变化。因此用低场核磁共振T2谱计算分形维数,可以对含水合物的多孔介质的几何特性进行动态而又准确的研究。
发明内容
本发明提供了一种基于核磁共振T2谱计算含水合物多孔介质的分形维数的方法。该方法利用低场核磁共振对多孔介质内水合物的生成过程进行监测,通过T2分布谱获取水合物生成过程中多孔介质的孔隙流体含量以及分布变化,对数据进行分析,定量得出水合物生成过程中的分形维数。
本发明提供了一种基于核磁共振T2谱计算含水合物多孔介质的分形维数的方法,该方法可以准确地确定水合物生成过程中多孔介质内流体含量以及分布,从而基于线性回归分析直接计算出分形维数。
本发明的基于核磁共振T2谱计算含水合物多孔介质分形维数的方法,采用低场核磁共振实现多孔介质内孔隙流体含量以及分布的确定,方法首先利用低场核磁共振获得容器中完全水饱和状态下的T2分布谱一,然后驱出部分水后获得T2分布谱二,最后获得水合物生成过程中的T2分布谱三,基于T2数据进行分形维数的计算。
具体步骤如下:
步骤一:获得完全水饱和状态下的T2分布谱一和部分水饱和状态下的T2分布谱二
将石英砂作为多孔介质填入低场核磁共振容器中并压实,注水实现完全水饱和后,控制温度在1℃并维持稳定,用低场核磁共振获取多孔介质完全水饱和状态下的T2分布谱一,标定此时孔隙水体积为V0,该体积也是多孔介质的孔隙体积;驱出部分水,用低场核磁共振获取此时部分水饱和状态下的T2分布谱二,标定此时孔隙水体积为V0’,则初始水饱和度为W=V0’/V0。
步骤二:获得水合物生成过程中的T2分布谱三
注入二氧化碳气体至系统压力高于相平衡压力1350kPa,并保持压力恒定,进行水合物的生成过程,用低场核磁共振获取水合物生成过程中的T2分布谱三,标定各个时刻的孔隙水体积为Vt;
判断T2分布谱三是否存在不连续阶段。若存在,则进行步骤三;若不存在,则用低场核磁共振获取水合物生成结束后的T2分布谱四,将此时的累积含水量分布曲线的最大值投射到水合物生成过程中某一时刻的累积含水量分布曲线上,交点所对应的T2值即为该时刻含水合物多孔介质多对应的T2截止值(T2cutoff)。根据T2截止值划分出两种孔隙:一种是在水合物生成过程中不能转化为水合物的束缚水孔隙(小孔隙,T2<T2cutoff),另一种是在水合物生成过程中可以完全转化为水合物的自由水孔隙(大孔隙,T2>T2cutoff)。计算出孔隙水体积Vt所占孔隙体积V0的百分比St,并排除束缚水空间的影响,对大孔隙部分所对应的lgT2与lgSt数据进行线性回归分析,得到水合物生成过程中的分形维数D。
步骤三:基于T2数据进行分形维数计算
定义不连续阶段的T2最大值为T2c。对比步骤一和步骤二的多孔介质完全水饱和状态下的T2分布谱一以及水合物生成过程中的T2分布谱三,用公式(1) 计算出相应时刻孔隙流体体积所占孔隙体积的百分比St:
依据多孔介质孔隙流体的核磁共振横向弛豫机制和多孔介质毛细管压力与孔径之间的关系,用公式(2)、(3)推导出多孔介质毛细管压力Pc与T2之间的关系(4):
其中Fs是几何形状因子,当孔隙为球状时Fs=3,孔隙为柱状时Fs=2;ρ是孔隙表面的T2弛豫强度(μm/ms);r是孔径(cm);Pc是当孔径为r时的毛细管压力(Mpa);σ为孔隙流体的表面张力(mN/m);θ为孔隙流体与孔隙表面间的接触角;为转换系数。
其中Pcmin为多孔介质中的最大孔径所对应的毛管压力,T2max为多孔介质中最大孔隙所对应的横向弛豫时间。
根据多孔介质分形理论,可以推导出多孔介质毛管压力曲线的分形表达式 (5)和多孔介质基于低场核磁共振T2谱的分形表达式(6):
对式(6)两边取对数可得:
lgSt=(3-D)lgT2+(D-3)lgT2max (7)
用上式(7)直接对T2>T2c的孔隙部分所对应的lgT2与lgSt数据进行线性回归分析,得到水合物生成过程中的分形维数D。
本发明的有益效果是,应用该方法可以对多孔介质内部结构的几何特性进行准确的表征,动态地监测水合物生成/分解过程对多孔介质的内部形状的影响。这种方法可以有效简化计算分形维数的步骤。
附图说明
图1是T2谱分布曲线连续的T2截止值确定方法示意图。
图2是T2谱分布曲线连续时的分形维数确定方法。
图3是T2谱分布曲线不连续的T2谱处理方法示意图。
图4是T2谱分布曲线不连续时的分形维数确定方法。
具体实施方式
下面结合技术方案和附图详细叙述本发明的对比例和实施例。
对比例
(1)获得完全水饱和状态下的T2分布谱一和部分水饱和状态下的T2分布谱二
将石英砂(粒径范围为80-120目)作为多孔介质填入容器中并压实,填砂模型长度为5cm,直径为2.5cm,体积为24.544cm3,孔隙度为42.03%。注水实现完全水饱和后,控制温度在1℃并维持稳定,用低场核磁共振获取多孔介质完全水饱和状态下的T2分布谱一。此时孔隙水体积V0=10.315ml,即多孔介质的孔隙体积为10.315cm3;驱出部分水,测定此时孔隙水体积V0’=4.722ml,则初始水饱和度为W=V0’/V0=45.78%。
(2)获得水合物生成过程中的T2分布谱三
用注气柱塞泵注入二氧化碳气体至系统压力达到3500kPa,并用泵保持压力恒定,进行水合物的生成过程,用低场核磁共振获取水合物生成过程中30min 时刻的T2分布谱三,如图1所示。判断该T2分布谱不存在不连续阶段,且此时孔隙水体积Vt=1.633ml。
(3)基于T2数据进行分形维数计算
计算出此时孔隙流体所占孔隙体积的百分比St=15.83%、T2cutoff=7.08ms、 T2max=10000ms,并直接对T2cutoff<T2<T2max所对应的lgT2与lgSv数据进行线性回归分析,得到水合物生成过程中的分形维数D=2.606,如图2所示。
实施例1
(1)获得完全水饱和状态下的T2分布谱一和部分水饱和状态下的T2分布谱二
将石英砂(粒径范围为80-120目)作为多孔介质填入容器中并压实,填砂模型长度为5cm,直径为2.5cm,体积为24.544cm3,孔隙度为42.03%。注水实现完全水饱和后,控制温度在1℃并维持稳定,用低场核磁共振获取多孔介质完全水饱和状态下的T2分布谱一。此时孔隙水体积V0=10.315ml,即多孔介质的孔隙体积为10.315cm3;驱出部分水,测定此时孔隙水体积V0’=4.722ml,则初始水饱和度为W=V0’/V0=45.78%。
(2)获得水合物生成过程中的T2分布谱三
用注气柱塞泵注入二氧化碳气体至系统压力达到3500kPa,并用泵保持压力恒定,进行水合物的生成过程,用低场核磁共振获取水合物生成过程中30min 时刻的T2分布谱二,如图3所示。判断该T2分布谱存在不连续阶段,且此时孔隙水体积Vt=1.494ml。
(3)基于T2数据进行分形维数计算
计算出此时孔隙流体所占孔隙体积的百分比St=14.48%、T2c=22.39ms、 T2max=10000ms,并直接对T2c<T2<T2max所对应的lgT2与lgSv数据进行线性回归分析,得到水合物生成过程中的分形维数D=2.665,如图4所示。
Claims (1)
1.基于核磁共振T2谱计算含水合物多孔介质分形维数的方法,其特征在于,包括以下步骤:
步骤一:获得完全水饱和状态下的T2分布谱一和部分水饱和状态下的T2分布谱二;
将石英砂作为多孔介质填入低场核磁共振容器中并压实,注水实现完全水饱和后,控制温度在1℃并维持稳定,用低场核磁共振获取多孔介质完全水饱和状态下的T2分布谱一,标定此时孔隙水体积为V0,该体积也是多孔介质的孔隙体积;驱出部分水,用低场核磁共振获取此时部分水饱和状态下的T2分布谱二,标定此时孔隙水体积为V0’,则初始水饱和度为W=V0’/V0;
步骤二:获得水合物生成过程中的T2分布谱三;
注入二氧化碳气体至系统压力高于相平衡压力1350kPa,并保持压力恒定,进行水合物的生成过程,用低场核磁共振获取水合物生成过程中的T2分布谱三,标定各个时刻的孔隙水体积为Vt;
判断T2分布谱三是否存在不连续阶段;若存在,则进行步骤三;若不存在,则用低场核磁共振获取水合物生成结束后的T2分布谱四,将此时的累积含水量分布曲线的最大值投射到水合物生成过程中任一时刻的累积含水量分布曲线上,交点所对应的T2值即为该时刻含水合物多孔介质多对应的T2截止值T2cutoff;根据T2截止值划分出两种孔隙:一种是在水合物生成过程中不能转化为水合物的束缚水孔隙,即小孔隙,T2<T2cutoff;另一种是在水合物生成过程中可以完全转化为水合物的自由水孔隙,即大孔隙,T2>T2cutoff;计算出孔隙水体积Vt所占孔隙体积V0的百分比St,并排除束缚水空间的影响,对大孔隙所对应的lgT2与lgSt数据进行线性回归分析,得到水合物生成过程中的分形维数D;
步骤三:基于T2数据进行分形维数计算;
定义不连续阶段的T2最大值为T2c;对比步骤一和步骤二的多孔介质完全水饱和状态下的T2分布谱一以及水合物生成过程中的T2分布谱三,用公式(1)计算出相应时刻孔隙流体体积所占孔隙体积的百分比St:
依据多孔介质孔隙流体的核磁共振横向弛豫机制和多孔介质毛细管压力与孔径之间的关系,用公式(2)、(3)推导出多孔介质毛细管压力Pc与T2之间的关系(4):
其中Fs是几何形状因子,当孔隙为球状时Fs=3,孔隙为柱状时Fs=2;ρ是孔隙表面的T2弛豫强度μm/ms;r是孔径cm;Pc是当孔径为r时的毛细管压力Mpa;σ为孔隙流体的表面张力mN/m;θ为孔隙流体与孔隙表面间的接触角;为转换系数;
其中Pcmin为多孔介质中的最大孔径所对应的毛管压力,T2max为多孔介质中最大孔隙所对应的横向弛豫时间;
根据多孔介质分形理论,推导出多孔介质毛管压力曲线的分形表达式(5)和多孔介质基于低场核磁共振T2谱的分形表达式(6):
对式(6)两边取对数可得:
lgSt=(3-D)lgT2+(D-3)lgT2max(7)
用上式(7)直接对T2>T2c的孔隙部分所对应的lgT2与lgSt数据进行线性回归分析,得到水合物生成过程中的分形维数D。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911076405.2A CN110702722B (zh) | 2019-11-06 | 2019-11-06 | 基于核磁共振t2谱计算含水合物多孔介质分形维数的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911076405.2A CN110702722B (zh) | 2019-11-06 | 2019-11-06 | 基于核磁共振t2谱计算含水合物多孔介质分形维数的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110702722A CN110702722A (zh) | 2020-01-17 |
CN110702722B true CN110702722B (zh) | 2021-01-19 |
Family
ID=69205370
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911076405.2A Active CN110702722B (zh) | 2019-11-06 | 2019-11-06 | 基于核磁共振t2谱计算含水合物多孔介质分形维数的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110702722B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111707695B (zh) * | 2020-06-16 | 2021-09-10 | 中国石油大学(北京) | 基于核磁共振数据确定储层渗透率的方法、装置及设备 |
CN112284999B (zh) * | 2020-10-26 | 2022-06-17 | 中国石油大学(华东) | 砂岩孔径分布测定方法及其应用 |
CN113237914B (zh) * | 2021-05-17 | 2021-12-28 | 中国科学院西北生态环境资源研究院 | 一种基于核磁信号强度测定水合物中未水合水含量的方法 |
Family Cites Families (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7771609B2 (en) * | 2002-08-16 | 2010-08-10 | Aerogel Technologies, Llc | Methods and compositions for preparing silica aerogels |
FR3033645B1 (fr) * | 2015-03-12 | 2018-10-12 | Institut National De Recherche En Sciences Et Technologies Pour L'environnement Et L'agriculture (Irstea) | Procede d'identification et de quantification des proportions relatives des differentes fractions cristallisees presentes dans un coulis d'hydrate et son utilisation pour les fluides frigoporteurs |
US10578568B2 (en) * | 2015-04-29 | 2020-03-03 | Colorado School Of Mines | Water/oil/gas emulsions/foams characterization using low field nuclear magnetic resonance |
WO2018024748A1 (en) * | 2016-08-01 | 2018-02-08 | Centre Hospitalier Universitaire D'angers | Multi-targeted fibrosis tests |
CN108802087A (zh) * | 2018-05-08 | 2018-11-13 | 中国石油大学(华东) | 一种基于低场核磁共振的岩石中天然气水合物生成效率的定量评价方法 |
CN109100812B (zh) * | 2018-07-25 | 2020-06-05 | 中国石油大学(北京) | 基于核磁共振的岩石孔隙分形维数评价方法及装置 |
CN109374489B (zh) * | 2018-08-21 | 2020-07-07 | 中国地质大学(武汉) | 联合x-ct技术的水合物沉积物nmr弛豫信号量标定装置和方法 |
-
2019
- 2019-11-06 CN CN201911076405.2A patent/CN110702722B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN110702722A (zh) | 2020-01-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110702722B (zh) | 基于核磁共振t2谱计算含水合物多孔介质分形维数的方法 | |
CN104196524B (zh) | 一种欠饱和煤储层开发的气水产出动态相渗曲线测量方法 | |
Liu et al. | Study on the seepage characteristics of coal based on the Kozeny-Carman equation and nuclear magnetic resonance experiment | |
Clarke et al. | Determination of the activation energy and intrinsic rate constant of methane gas hydrate decomposition | |
CN105507870A (zh) | 一种砂岩储层无填砂水力裂缝导流能力确定方法 | |
CN112360422B (zh) | 一种页岩气藏压裂水平井产量预测方法及系统 | |
AU2020103953A4 (en) | Method and system for predicting production of fractured horizontal well in shale gas reservoir | |
CN106706496A (zh) | 一种致密油/页岩油的纳米级含油孔径分布的测定方法 | |
CN107271314A (zh) | 一种测量煤岩吸附膨胀系数的方法 | |
CN109342297A (zh) | 基于压汞实验的煤中孔隙校正方法 | |
CN107703275A (zh) | 一种甲烷水合物相平衡研究的高压实验装置及方法 | |
CN103995961A (zh) | 一种测试井的生产数据分析方法 | |
CN111537544B (zh) | 一种提高核磁共振t2谱表征致密储层孔径分布精度的转换方法 | |
CN106124378A (zh) | 一种基于压汞法毛管压力曲线的复杂岩石分形维数分布谱构建方法 | |
CN108268712B (zh) | 核磁共振确定孔隙介质毛管压力的方法和装置 | |
Yang et al. | An integrated method of measuring gas permeability and diffusion coefficient simultaneously via pressure decay tests in shale | |
Wang et al. | Methane hydrate formation and decomposition properties during gas migration in porous medium | |
CN114741932A (zh) | 基于粒子群算法的磁纳米粒子体积分数优化方法 | |
Wu et al. | Elucidation of evolutionary characteristics of pore structure in anthracite through microwave-LN2 freeze–thaw cycling | |
CN115860266B (zh) | 一种页岩气/煤层气井产能评价方法、系统及电子设备 | |
CN111006985B (zh) | 地质条件下陆相湖盆致密储层孔喉有效性定量评价方法 | |
CN109799164B (zh) | 吸附气量的确定方法及装置 | |
Li et al. | Pore size distribution of a tight sandstone reservoir and its effect on micro pore‐throat structure: a case study of the Chang 7 member of the Xin'anbian Block, Ordos Basin, China | |
CN112485281B (zh) | 动态测定多孔介质中气体水合物饱和度和渗透率的方法 | |
Zheng et al. | Pore structure reconstruction and moisture migration in porous media |
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 |