CN110823780B - 一种饱和冻结岩石未冻水含量的计算方法 - Google Patents

一种饱和冻结岩石未冻水含量的计算方法 Download PDF

Info

Publication number
CN110823780B
CN110823780B CN201911133940.7A CN201911133940A CN110823780B CN 110823780 B CN110823780 B CN 110823780B CN 201911133940 A CN201911133940 A CN 201911133940A CN 110823780 B CN110823780 B CN 110823780B
Authority
CN
China
Prior art keywords
pore
rock
rock sample
porosity
radius
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
CN201911133940.7A
Other languages
English (en)
Other versions
CN110823780A (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.)
Central South University
Original Assignee
Central South University
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 Central South University filed Critical Central South University
Priority to CN201911133940.7A priority Critical patent/CN110823780B/zh
Publication of CN110823780A publication Critical patent/CN110823780A/zh
Application granted granted Critical
Publication of CN110823780B publication Critical patent/CN110823780B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume, or surface-area of porous materials
    • G01N15/08—Investigating permeability, pore-volume, or surface area of porous materials
    • G01N15/088—Investigating volume, surface area, size or distribution of pores; Porosimetry
    • 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

Abstract

本发明公开了一种饱和冻结岩石未冻水含量的计算方法。该方法首先利用核磁共振技术获得岩石试样的孔隙度和T2谱图，根据岩石T2谱图与孔隙尺寸分布的关系，计算得出孔隙度累积分布曲线的解析表达式；将孔隙累积分布曲线的解析表达式联立冻结过程孔隙尺寸与相变温度的数学模型，从而得出饱和岩石未冻水含量的计算公式。该方法只需通过核磁共振技术获得岩石试样的孔隙度和孔隙度累积分布曲线，即可计算出不同温度下饱和岩石的未冻水含量。该方法适用于岩石类脆性材料的未冻水含量计算，具有计算过程简单、便捷，且计算结果直观、可靠的特点。

Description

一种饱和冻结岩石未冻水含量的计算方法
技术领域
本发明属于岩土工程技术领域，具体涉及一种饱和冻结岩石未冻水含量的计算方法。
背景技术
在寒区矿山开采及岩土工程的建设过程中，冻融循环作用对岩石力学性能劣化、岩体稳定性及地质灾害防治产生了重要影响。在低温环境下，岩石孔隙中的水发生相变，水结冰后体积膨胀约为9％，从而产生冻胀力，进而对岩石产生冻胀变形和损伤。而在冻结过程中，未冻水是液态水迁移的主要来源，对冰的形成、冻胀力及其造成的变形程度等都有着重要的影响。可以说，孔隙中的未冻水含量是影响岩石低温冻胀变形的主要因素，它对冻胀力学、传热、水分迁移等问题都有直接的关联，对于进行寒区岩体工程低温多场耦合计算参数取值至关重要。
已公开的未冻水含量计算模型主要是针对冻土介质，考虑到岩石与土体的本质区别，这些数学计算模型应用于冻岩中的未冻水含量计算存在一定的局限性。
发明内容
本发明的目的是提供一种饱和冻结岩石未冻水含量的计算方法。该方法在基于热力学理论推导出冻结过程孔隙尺寸与相变温度的数学关系的基础上，结合核磁共振技术获得的岩石累积孔隙度分布曲线，建立了计算冻结过程未冻水含量的数学模型。
为解决上述技术问题，本发明采用如下技术方案：
一种饱和冻结岩石未冻水含量的计算方法，包括如下步骤：
步骤1：选取并加工岩石试样；
步骤2：采用真空饱和装置对岩石试样进行真空饱水处理；
步骤3：采用低场核磁共振测试设备，测试出岩石试样的孔隙度、T2谱图和孔隙度累积分布曲线；
步骤4：获得岩石试样的孔隙度累积分布曲线的解析表达式；
步骤5：计算不同温度下未冻水膜的厚度。
步骤6：计算不同温度下孔隙水结冰的临界半径。
步骤7：将上述步骤中的参数代入式(1)计算出不同温度T下饱和冻结岩石未冻水含量与岩石试样的孔隙度的比值Sw，
Figure GDA0002979506070000021
其中，Sw为饱和冻结岩石中的未冻水含量占孔隙度的比例；
Vw为饱和冻结岩石未冻水含量；
n为岩石试样的孔隙度；
r为岩石试样的孔隙半径；
φ(r)为岩石试样的孔隙度累积曲线的解析表达式；
rc为孔隙水结冰的临界半径；
φ(rc)为岩石试样中小于孔隙水结冰的临界半径rc的孔隙含量；
rwf为未冻水膜的厚度。
进一步的，所述步骤3中岩石试样的孔隙度、T2谱图和孔隙度累积分布曲线由核磁共振测试设备测得；根据式(2)中岩石试样的T2弛豫时间与岩石试样孔隙尺寸的关系式，将T2谱图转化为不同孔隙尺寸的孔隙度累积分布曲线。
r＝ρFsT2 (2)
其中，T2为岩石试样的横向弛豫时间；
Fs为孔隙几何形状因子；
ρ为表面弛豫率；
r为孔隙半径。
进一步的，所述步骤3中孔隙度累积分布曲线即是表征岩石试样中不同尺寸的孔隙含量的曲线。
进一步的，所述步骤4中岩石试样不同孔隙尺寸的孔隙度累积曲线的解析表达式如式(3)所示
Figure GDA0002979506070000022
其中，
r为孔隙半径；
r1为岩石试样孔隙半径范围的最小值，取T2弛豫时间的最小截至值代入式(2)计算得到；
r2为岩石试样孔隙半径范围的最大值，取T2弛豫时间的最大截至值代入式(2)计算得到；
Figure GDA0002979506070000031
为孔隙半径尺寸在区间[r1,r2]内的孔隙含量表达式，由孔隙度累积分布曲线拟合而得；
n1为岩石试样总的孔隙含量，可将孔隙半径r＝r2代入孔隙含量表达式中计算得出。
进一步的，所述步骤5中未冻水膜的厚度rwf可由式(4)计算，
Figure GDA0002979506070000032
其中，T0为0℃的开氏温度；
T为岩石试样冻结时的温度，一般取0℃至实验最低温度的开氏温度；
未冻水膜的厚度rwf的单位为nm。
进一步的，所述步骤6中孔隙水结冰的临界半径rc可由式(5)计算，
Figure GDA0002979506070000033
其中，孔隙水结冰的临界半径rc的单位为nm。
进一步的，本发明适用于饱和岩石类脆性材料的未冻水含量计算。
与现有技术相比，本发明的有益效果在于：利用核磁共振技术，通过岩石试样的T2分布曲线与孔隙尺寸分布的关系，联立冻结过程孔隙尺寸与相变温度的数学公式，得出饱和岩石未冻水含量的计算方式。由于核磁共振技术的高效、便捷、无损等优点，能够清晰直观的获取岩石试样的孔隙尺寸分布，且操作难度较低。
本发明具有计算方法简单、便捷，计算结果直观、可靠的特点，可为岩石低温冻胀过程中液态水迁移机制、多场耦合计算参数选取等研究提供理论参考。
附图说明
图1为本发明的流程图；
图2为砂岩试样的T2谱图和孔隙累积分布曲线。
具体实施方式
以下结合附图，对本发明实施例中的技术方案进行清楚、完整地描述。显然所描述的实施例仅为本发明的优选例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
参见图1，一种饱和冻结岩石未冻水含量的计算方法，包括如下步骤：
(1)选取砂岩为岩石试样，根据常规岩石力学试验标准，将岩样加工成高径比2:1(100mm×Φ50mm)的圆柱体。这里需要解释说明的是，岩石试样采用砂岩仅仅为示例性说明，本发明中的岩石试样也可以采用其他种类的岩石。
(2)采用真空饱和装置，设备设置抽气时间和压力后，对砂岩试样进行真空饱水处理。
(3)利用低场核磁共振设备，对真空饱水处理后的砂岩试样进行测试，得到了饱和砂岩试样的孔隙度、T2谱图和孔隙度累积分布曲线，如图2所示。进而可由式(2)计算得到砂岩试样不同孔隙尺寸的孔隙度累积分布曲线，式(2)为
r＝ρFsT2 (2)
其中，T2为岩石试样的横向弛豫时间；
Fs为孔隙几何形状因子；
ρ为表面弛豫率；
r为孔隙半径。
由图2可知砂岩试样的弛豫时间范围为0.099ms～413.56ms，则T2的最小截止值为0.099ms，最大截止值为413.56ms；假设砂岩试样的孔隙为理想圆柱体，取孔隙几何形状因子Fs＝2，砂岩的表面弛豫率ρ＝10μm/s，代入式(2)可得到砂岩试样的T2与孔隙半径的线性关系式，如式(6)所示
r＝20T2 (6)
再将T2的最小截止值和最大截止值代入式(3)，即可计算出砂岩试样的孔隙尺寸的最小孔隙半径为r1＝20×0.099＝1.98nm，最大孔隙半径为r2＝20×413.56＝8271.20nm。
根据式(3)，可得到砂岩试样不同孔隙尺寸的孔隙度累积分布曲线(试验累计孔隙度)。
对砂岩试样不同孔隙尺寸的孔隙累积分布曲线进行拟合。根据拟合结果可得到孔隙半径尺寸在区间[r1,r2]内砂岩试样的孔隙含量表达式，如式(7)所示；
Figure GDA0002979506070000041
其中，
Figure GDA0002979506070000042
为孔隙半径尺寸在区间[r1,r2]内的孔隙含量表达式；r为岩石试样的孔隙半径。
将r2＝8271.20nm代入式(7)即可计算出砂岩试样总的孔隙含量n1＝54.18％。
将式(7)和n1＝54.18％、r1＝1.98nm以及r2＝8271.20nm代入式(3)即可得到砂岩试样不同孔隙尺寸的孔隙度累积分布曲线的解析式，即式(8)：
Figure GDA0002979506070000051
其中，φ(r)为岩石试样的孔隙度累积曲线的解析式；
r为孔隙半径；
(4)计算不同温度下未冻水膜的厚度rwf和临界半径rc
分别取温度为-0.5℃、-1℃、-2℃、-5℃、-10℃、-20℃、-30℃，代入式(4)中，计算出不同温度下未冻水膜厚度rwf，式(4)为：
Figure GDA0002979506070000052
其中，T0为0℃的开氏温度；
T为岩石试样冻结时的温度，一般取0℃至实验最低温度的开氏温度；
未冻水膜的厚度rwf的单位为nm；
未冻水膜厚度rwf的计算结果如表1所示。
表1不同温度条件下未冻水膜厚度rwf
Figure GDA0002979506070000053
分别取温度为-0.5℃、-1℃、-2℃、-5℃、-10℃、-20℃、-30℃，代入式(5)中计算不同温度下孔隙水结冰的临界半径rc，式(5)为：
Figure GDA0002979506070000054
其中，孔隙水结冰的临界半径rc的单位为nm。
孔隙水结冰的临界半径rc的计算结果如表2所示。
表2不同温度条件下孔隙水结冰的临界半径rc
Figure GDA0002979506070000055
(5)将表1和表2中的参数带入式(1)，式(1)为
Figure GDA0002979506070000061
其中，Sw为饱和冻结岩石未冻水含量占总孔隙体积的比例；
Vw为饱和冻结岩石未冻水含量；
n为岩石试样的孔隙总体积；
r为岩石试样的孔隙半径；
φ(r)为岩石试样的孔隙度累积曲线公式；
rc为孔隙水结冰的临界半径；
φ(rc)为岩石试样中小于孔隙水结冰的临界半径rc的孔隙含量；
rwf为未冻水膜的厚度。
计算出不同温度下砂岩试样的未冻水含量占孔隙总体积的比值，计算结果如表3所示。
表3不同温度条件下未冻水含量占孔隙中体积比例
Figure GDA0002979506070000062
上述实施例仅仅是清楚地说明本发明所作的举例，而非对实施方式的限定。对于所属领域的普通技术人员来说，在上述说明的基础上还可以做出其它不同形式的变化或变动。这里也无需也无法对所有的实施例予以穷举。而由此所引申出的显而易见的变化或变动仍处于本发明的保护范围之中。

Claims (3)

1.一种饱和冻结岩石未冻水含量的计算方法，其特征在于：包括如下步骤：
步骤1：选取并加工岩石试样；
步骤2：对岩石试样进行真空饱水处理；
步骤3：采用低场核磁共振测试设备，测试出岩石试样的孔隙度、T2谱图和孔隙度累积分布曲线；
步骤4：获得岩石试样的孔隙度累积分布曲线的解析表达式；
步骤5：计算不同温度下未冻水膜的厚度；
步骤6：计算不同温度下孔隙水结冰的临界半径；
步骤7：将上述步骤中的参数代入式(1)中，计算出不同温度T下饱和冻结岩石未冻水含量与岩石试样的孔隙度的比值Sw，
Figure FDA0002979506060000011
其中，Sw为饱和冻结岩石未冻水含量占孔隙度的比例；
Vw为饱和冻结岩石未冻水含量；
n为岩石试样的孔隙度；
r为岩石试样的孔隙半径；
φ(r)为岩石试样的孔隙度累积曲线的解析表达式；
rc为孔隙水结冰的临界半径；
φ(rc)为岩石试样中小于孔隙水结冰的临界半径rc的孔隙含量；
rwf为未冻水膜的厚度；
孔隙半径依据步骤3中核磁共振测试设备获得的T2谱图计算得到，计算公式如下：
r＝ρFsT2 (2)
其中，T2为岩石试样的横向弛豫时间；
Fs为孔隙几何形状因子；
ρ为表面弛豫率；
r为孔隙半径；
所述步骤4中岩石试样的孔隙度累积分布曲线的解析表达式如式(3)所示：
Figure FDA0002979506060000021
其中，r为孔隙半径；
r1为岩石试样孔隙半径范围的最小值，取T2弛豫时间的最小截至值代入式(2)计算得到；
r2为岩石试样孔隙半径范围的最大值，取T2弛豫时间的最大截至值代入式(2)计算得到；
Figure FDA0002979506060000022
为孔隙半径尺寸在区间[r1,r2]内的孔隙含量表达式，由孔隙度累积分布曲线拟合而得；
n1为岩石试样总的孔隙含量，
Figure FDA0002979506060000023
所述步骤5中未冻水膜的厚度rwf由式(4)计算，
Figure FDA0002979506060000024
其中，T0为0℃的开氏温度；
T为岩石试样冻结时的温度，取0℃至实验最低温度的开氏温度；
所述步骤6中孔隙水结冰的临界半径rc由式(5)计算，
Figure FDA0002979506060000025
2.根据权利要求1所述的计算方法，其特征是：所述岩石试样为脆性岩石材料。
3.根据权利要求1所述的计算方法，其特征是：孔隙几何形状因子Fs＝2，表面弛豫率ρ＝10μm/s。
CN201911133940.7A 2019-11-19 2019-11-19 一种饱和冻结岩石未冻水含量的计算方法 Active CN110823780B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911133940.7A CN110823780B (zh) 2019-11-19 2019-11-19 一种饱和冻结岩石未冻水含量的计算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911133940.7A CN110823780B (zh) 2019-11-19 2019-11-19 一种饱和冻结岩石未冻水含量的计算方法

Publications (2)

Publication Number Publication Date
CN110823780A CN110823780A (zh) 2020-02-21
CN110823780B true CN110823780B (zh) 2021-05-18

Family

ID=69556783

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911133940.7A Active CN110823780B (zh) 2019-11-19 2019-11-19 一种饱和冻结岩石未冻水含量的计算方法

Country Status (1)

Country Link
CN (1) CN110823780B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113237914B (zh) * 2021-05-17 2021-12-28 中国科学院西北生态环境资源研究院 一种基于核磁信号强度测定水合物中未水合水含量的方法
CN113433292A (zh) * 2021-06-25 2021-09-24 中国矿业大学 冻融循环作用下散体试样宏微观结构演化试验系统及方法
CN114384232B (zh) * 2022-01-20 2022-10-28 中国科学院西北生态环境资源研究院 一种基于吸附和毛细耦合效应的未冻水模型

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104634718A (zh) * 2015-03-05 2015-05-20 中国石油大学(华东) 应用核磁共振表征致密砂岩孔径分布的标定方法
CN104777181A (zh) * 2014-12-31 2015-07-15 中国石油天然气股份有限公司 致密油核磁共振t2截止值及流体饱和度确定方法、装置
CN105241920A (zh) * 2015-11-06 2016-01-13 天津城建大学 利用比热计算确定土冻结过程中未冻水含量的方法
CN105466830A (zh) * 2015-11-19 2016-04-06 中国石油天然气股份有限公司 储层砂岩孔隙喉道尺寸分布识别方法
CN106872529A (zh) * 2017-04-07 2017-06-20 天津城建大学 通过测量电阻率确定冻土未冻水含量的方法
RU2654832C1 (ru) * 2017-06-21 2018-05-22 Автономная некоммерческая образовательная организация высшего образования "Сколковский институт науки и технологий" Способ определения содержания незамерзшей воды в мерзлых грунтах
CN108107068A (zh) * 2017-12-11 2018-06-01 内蒙古工业大学 在冻融循环下轻骨料混凝土孔隙结冰速率的测试方法
CN110231272A (zh) * 2019-07-09 2019-09-13 中国地质大学（北京） 致密砂岩孔径与核磁共振t2值转换关系的确定方法及系统

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10353104B2 (en) * 2012-03-29 2019-07-16 Baker Hughes, A Ge Company, Llc Carbonate permeability by pore typing

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104777181A (zh) * 2014-12-31 2015-07-15 中国石油天然气股份有限公司 致密油核磁共振t2截止值及流体饱和度确定方法、装置
CN104634718A (zh) * 2015-03-05 2015-05-20 中国石油大学(华东) 应用核磁共振表征致密砂岩孔径分布的标定方法
CN105241920A (zh) * 2015-11-06 2016-01-13 天津城建大学 利用比热计算确定土冻结过程中未冻水含量的方法
CN105466830A (zh) * 2015-11-19 2016-04-06 中国石油天然气股份有限公司 储层砂岩孔隙喉道尺寸分布识别方法
CN106872529A (zh) * 2017-04-07 2017-06-20 天津城建大学 通过测量电阻率确定冻土未冻水含量的方法
RU2654832C1 (ru) * 2017-06-21 2018-05-22 Автономная некоммерческая образовательная организация высшего образования "Сколковский институт науки и технологий" Способ определения содержания незамерзшей воды в мерзлых грунтах
CN108107068A (zh) * 2017-12-11 2018-06-01 内蒙古工业大学 在冻融循环下轻骨料混凝土孔隙结冰速率的测试方法
CN110231272A (zh) * 2019-07-09 2019-09-13 中国地质大学（北京） 致密砂岩孔径与核磁共振t2值转换关系的确定方法及系统

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
Experimental investigation on the effects of ambient freeze–thaw cycling on creep mechanical properties of sandstone under step loading;Jielin Li,et al.;《IEEE ACCESS》;20181231;第1-9页 *
冻融作用下岩石细观结构损伤的低场核磁共振研究;李杰林 等;《西安科技大学学报》;20180331;第38卷(第2期);第266-272页 *
冻融作用下砂岩孔隙结构损伤特征研究;李杰林 等;《岩土力学》;20190930;第40卷(第9期);第3524-3532页 *

Also Published As

Publication number Publication date
CN110823780A (zh) 2020-02-21

Similar Documents

Publication Publication Date Title
CN110823780B (zh) 一种饱和冻结岩石未冻水含量的计算方法
Choukroun et al. Thermodynamic data and modeling of the water and ammonia-water phase diagrams up to 2.2 GPa for planetary geophysics
Li et al. Experimental measurement and mathematical model of permeability with methane hydrate in quartz sands
Staykova et al. Formation of porous gas hydrates from ice powders: diffraction experiments and multistage model
CN104196524B (zh) 一种欠饱和煤储层开发的气水产出动态相渗曲线测量方法
Sabil et al. Estimations of enthalpies of dissociation of simple and mixed carbon dioxide hydrates from phase equilibrium data
CN106481332B (zh) 用于确定页岩气多段压裂水平井内外区动态储量的方法
Li et al. Permeability measurements of quartz sands with methane hydrate
Shiozaki et al. Measuring the heat capacity in a Bose-Einstein condensation using global variables
Knoebel et al. Second virial coefficients of binary mixtures of benzene with methanol, ethanol, acetone, and diethyl ether
Sawamura et al. Solubility of sodium chloride in water under high pressure
Tahavvor et al. Experimental and numerical study of frost formation by natural convection over a cold horizontal circular cylinder
Ahn et al. Experimental characterization of production behaviour accompanying the hydrate reformation in methane-hydrate-bearing sediments
Gupta et al. Shock-induced vaporization of anhydrite and global cooling from the K/T impact
Wilson et al. The influence of strain rate and presence of dispersed second phases on the deformation behaviour of polycrystalline D2O ice
Sakoda et al. Transient temperature and pressure behavior of high-pressure 100 MPa hydrogen during discharge through orifices
Jaupart Archean thermal regime and stabilization of the cratons
Detian et al. Very low gas flow measurements for UHV/XHV and leak calibration
Legatski et al. Dispersion coefficients for gases flowing in consolidated porous media
Wei et al. Mechanisms of the memory effect of clathrate hydrates
Diller The specific heats (Cv) of dense simple fluids
CN106404610B (zh) 一种水合物动力学抑制剂的评价方法及其在筛选中的应用
Nakagawa et al. Dissociation and specific heats of gas hydrates under submarine and sublacustrine environments
Shen et al. Novel method for calculating the effective stress coefficient in a tight sandstone reservoir
Li et al. Systematic study of steam–water capillary pressure

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