CN106769771B - 一种基于低场核磁共振技术的非饱和土渗透系数的测量方法 - Google Patents
一种基于低场核磁共振技术的非饱和土渗透系数的测量方法 Download PDFInfo
- Publication number
- CN106769771B CN106769771B CN201710012453.XA CN201710012453A CN106769771B CN 106769771 B CN106769771 B CN 106769771B CN 201710012453 A CN201710012453 A CN 201710012453A CN 106769771 B CN106769771 B CN 106769771B
- Authority
- CN
- China
- Prior art keywords
- sample
- nuclear magnetic
- section
- moisture content
- signal amount
- 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.)
- Expired - Fee Related
Links
- 238000001764 infiltration Methods 0.000 title claims abstract description 17
- 230000008595 infiltration Effects 0.000 title claims abstract description 17
- 239000002689 soil Substances 0.000 title claims abstract description 11
- 238000005481 NMR spectroscopy Methods 0.000 title claims abstract description 10
- 238000005516 engineering process Methods 0.000 title claims abstract description 6
- 238000000691 measurement method Methods 0.000 title claims abstract description 5
- 238000000034 method Methods 0.000 claims abstract description 22
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 10
- 238000012360 testing method Methods 0.000 claims description 8
- 238000004364 calculation method Methods 0.000 claims description 4
- 230000035699 permeability Effects 0.000 claims description 3
- 229920006395 saturated elastomer Polymers 0.000 claims 1
- 238000011156 evaluation Methods 0.000 abstract 1
- 239000000523 sample Substances 0.000 description 41
- 238000001514 detection method Methods 0.000 description 4
- 238000002474 experimental method Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 230000005311 nuclear magnetism Effects 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 239000004744 fabric Substances 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000009659 non-destructive testing Methods 0.000 description 1
- 230000003204 osmotic effect Effects 0.000 description 1
- 230000035515 penetration Effects 0.000 description 1
- 238000009738 saturating Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
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
-
- 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
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- General Physics & Mathematics (AREA)
- Pathology (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- High Energy & Nuclear Physics (AREA)
- Immunology (AREA)
- Health & Medical Sciences (AREA)
- Engineering & Computer Science (AREA)
- Environmental & Geological Engineering (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geochemistry & Mineralogy (AREA)
- Geology (AREA)
- Dispersion Chemistry (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明提供一种基于低场核磁共振技术的非饱和土渗透系数的测量方法。该方法是以低场核磁共振技术为基础,结合瞬时截面法,进行非饱和渗透系数评价的方法。本发明以核磁信号量的大小为依据,监测试样非饱和渗流状态下多个截面的核磁信号量的变化,依据核磁信号量的大小对试样的质量含水率分布进行推算,结合瞬时截面法,通过Darcy定律计算得到非饱和渗透系数。
Description
技术领域
本发明属于岩土工程领域,涉及一种基于低场核磁共振技术的非饱和土渗透系数的测量方法。
背景技术
非饱和土在自然界广泛分布,许多工程问题都涉及到非饱和土的渗流分析。实验室常常采用瞬时截面法量测非饱和土的渗透系数。瞬时截面法需实时监测试样渗流状态下多个截面的质量含水率(吸力)的变化,最后通过Darcy定律计算得到渗透系数,即:
式中:A为所测试样的横截面积,Q为Δt时间内流过试样断面面积A的水的流量,;it为t时刻试样内一点的水力梯度。
显然,根据公式(1)并通过相关实验可以得到渗透系数kw。瞬时截面法的困难之处在于截面质量含水率(吸力)的测定。目前的截面质量含水率的检测手段都需要检测探头与土样直接贴合,不但使土样截面积发生变化,还会因接触面存在缝隙而影响渗流的均匀性,导致水的流速计算产生误差。另外,传统测量探头尺寸较大,导致各个截面间隔较大,降低了试验的精度。因此,现有渗透试验装置不是对实现自然界中的渗透特性的真实模拟。
发明内容
本发明所要解决的问题是提供一种能无损检测土样渗透系数的方法,克服现有技术存在的上述技术问题。
为实现上述目的,本发明采用以下技术方案予以解决:
1)基于核磁共振信号量分布曲线建立质量含水率分布曲线;
2)基于持水特征曲线和质量含水率分布曲线建立吸力分布曲线;
3)以质量含水率分布曲线和吸力分布曲线为基础,通过Darcy定律计算非饱和渗透系数。
建立质量含水率分布曲线的具体方法如下:
取试样第一个被测截面即进水截面在试验过程中的核磁信号量的均值S′作为标准核磁信号量,可以表示为:
式中,S0,i为试样进水截面、在第i次检测时的核磁信号量。
根据试样的饱和质量含水率,将S′标定为试样饱和质量含水率所对应的核磁信号量。试样饱和质量含水率与其对应的核磁信号量的比例关系可以表示为:
式中,ω为试样的饱和质量含水率。
试样其他截面的在t时刻的质量含水率可以表示为:
ωj,t=Sj,t·X (3)
式中,ωj,t为试样j截面在t时刻的质量含水率,Sj,t为试样j截面在t时刻的核磁信号量。
建立吸力分布曲线的具体方法如下:
根据试样的持水特征曲线(通过其他相关试验获得或参考相关材料的文献资料),将所得的质量含水率分布曲线转换为相应的吸力分布曲线。
计算渗透系数的具体方法如下:
(1)计算瞬时水力梯度
在时刻t,试样内一点的水力梯度可以表示为水头剖面的坡度:
iw=dhw/dx (4)
这里:iw——t时刻试样内点w的水力梯度;
dhw/dx——水头剖面中点w的坡度;
x——考虑点的位置。
(2)计算瞬时流量
某点的瞬时流量Q等于Δt时间内流过试样断面面积A的水的体积。在Δt时间内,通过试样某点的水的总体积可以通过体积质量含水率剖面在Δt时间间隔内的积分求得:
这里:θ——体积质量含水率;
x——考虑点的位置;
L——试样长度。
(3)平均渗透系数计算
有了瞬时流量和瞬时梯度,根据Darcy定律,可以求出Δt时间内的平均渗透系数:
这里:A——试样截面面积;
与现有技术相比,本发明具有以下有益效果:
本发明以核磁信号量的大小为依据,监测试样非饱和渗流状态下多个截面的核磁信号量的变化,依据核磁信号量的大小对试样的质量含水率分布进行推算,结合瞬时截面法,通过Darcy定律计算得到非饱和渗透系数。由于本发明采用核磁共振技术间接获取试样截面质量含水率分布,达到了不影响渗流形态、缩小试样尺寸、缩短试验周期的目的。
附图说明
图1为本发明核磁信号量分布曲线;
图2为本发明质量含水率分布曲线;
图3为本发明持水特征曲线;
图4为本发明吸力分布曲线;
图5为本发明水力梯度计算示意图;
图6为本发明流量计算示意图;
图7为本发明渗透系数与吸力的关系效果图;
具体实施方式
下面结合附图对本发明做进一步详细说明:
参见图1至图7,本发明包括以下步骤:1、建立质量含水率分布曲线;2、建立吸力分布曲线;3、通过Darcy定律计算非饱和渗透系数。
1、建立质量含水率分布曲线
试样的核磁信号分布曲线如图1所示。取试样第一个被测截面即进水截面在试验过程中的核磁信号量的均值S′作为标准核磁信号量,可以表示为:
式中,S0,i为试样进水截面、在第i次检测时的核磁信号量。
根据试样的饱和质量含水率,将S′标定为试样饱和质量含水率所对应的核磁信号量。试样饱和质量含水率与其对应的核磁信号量的比例关系可以表示为:
式中,ω为试样的饱和质量含水率。
试样其他截面的在t时刻的质量含水率如图2所示,可以表示为:
ωj,t=Sj,t·X (3)
式中,ωj,t为试样j截面在t时刻的质量含水率,Sj,t为试样j截面在t时刻的核磁信号量。
2、建立吸力分布曲线的具体方法如下:
根据试样的持水特征曲线如图3所示,将所得的质量含水率分布曲线转换为相应的吸力分布曲线,如图4所示。
3、计算渗透系数的具体方法如下:
(1)计算瞬时水力梯度
在时刻t,试样内一点的水力梯度可以表示为水头剖面的坡度,如图5所示:
iw=dhw/dx (4)
这里:iw——t时刻试样内点w的水力梯度;
dhw/dx——水头剖面中点w的坡度;
x——考虑点的位置。
(2)计算瞬时流量
某点的瞬时流量Q等于Δt时间内流过试样断面面积A的水的体积。在Δt时间内,通过试样某点的水的总体积可以通过体积质量含水率剖面在Δt时间间隔内的积分求得,如图6所示:
这里:θ——体积质量含水率;
x——考虑点的位置;
L——试样长度。
(3)平均渗透系数计算
有了瞬时流量和瞬时梯度,根据Darcy定律,可以求出Δt时间内的平均渗透系数,如图7所示:
这里:A——试样截面面积;
本发明具有以下优点:
本发明以核磁信号量的大小为依据,监测试样非饱和渗流状态下多个截面的核磁信号量的变化,依据核磁信号量的大小对试样的质量含水率分布进行推算,结合瞬时截面法,通过Darcy定律计算得到非饱和渗透系数。由于本发明采用核磁共振技术间接获取试样截面质量含水率分布,达到了不影响渗流形态、缩小试样尺寸、缩短试验周期的目的。
以上内容仅为描述本发明的思想,不能以此限定本发明的保护范围,凡是在本发明提出的技术思想的基础上所做的任何具体实施的技术方案的改动,均落入本发明的权利要求书的保护范围之内。
Claims (1)
1.一种基于低场核磁共振技术的非饱和土渗透系数的测量方法,其特征在于,包括以下步骤:
1)基于核磁共振信号量分布曲线建立质量含水率分布曲线;
2)基于持水特征曲线和质量含水率分布曲线建立吸力分布曲线;
3)以质量含水率分布曲线和吸力分布曲线为基础,通过Darcy定律计算非饱和渗透系数;
所述步骤1)中,建立质量含水率分布曲线的具体方法如下:
取试样第一个被测截面即进水截面在试验过程中的核磁信号量的均值S′作为标准核磁信号量,表示为:
式中,S0,i为试样进水截面在第i次检测时的核磁信号量,
根据试样的饱和质量含水率,将S′标定为试样饱和质量含水率所对应的核磁信号量;
试样饱和质量含水率与其对应的核磁信号量的比例关系表示为:
式中,ω为试样的饱和质量含水率;
试样其他截面的在t时刻的质量含水率表示为:
ωj,t=Sj,t·X (3)
式中,ωj,t为试样j截面在t时刻的质量含水率,Sj,t为试样j截面在t时刻的核磁信号量所述步骤2)中,建立吸力分布曲线的具体方法如下:
根据试样的持水特征曲线,将根据步骤1)所得的质量含水率分布曲线转换为相应的吸力分布曲线;
所述步骤3)中,计算渗透系数的具体方法如下:
(1)计算瞬时水力梯度
在时刻t,试样内一点的水力梯度表示为水头剖面的坡度:
iw=dhw/dx (4)
这里:iw——t时刻试样内点w的水力梯度;
dhw/dx——水头剖面中点w的坡度;
x——点的位置;
(2)计算瞬时流量
某点的瞬时流量Q等于Δt时间内流过试样断面面积A的水的体积;在Δt时间内,通过试样某点的水的总体积通过体积质量含水率剖面在Δt时间间隔内的积分求得:
这里:θ——体积质量含水率;
x——点的位置;
L——试样长度;
(3)平均渗透系数计算
有了瞬时流量和瞬时梯度,根据Darcy定律,求出Δt时间内的平均渗透系数:
这里:A——试样截面面积。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710012453.XA CN106769771B (zh) | 2017-01-09 | 2017-01-09 | 一种基于低场核磁共振技术的非饱和土渗透系数的测量方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710012453.XA CN106769771B (zh) | 2017-01-09 | 2017-01-09 | 一种基于低场核磁共振技术的非饱和土渗透系数的测量方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106769771A CN106769771A (zh) | 2017-05-31 |
CN106769771B true CN106769771B (zh) | 2019-07-05 |
Family
ID=58950370
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710012453.XA Expired - Fee Related CN106769771B (zh) | 2017-01-09 | 2017-01-09 | 一种基于低场核磁共振技术的非饱和土渗透系数的测量方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106769771B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111678853A (zh) * | 2020-06-19 | 2020-09-18 | 中国科学院武汉岩土力学研究所 | 一种测量高压实膨胀土饱和渗透系数的方法 |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108918379B (zh) * | 2018-06-22 | 2021-09-03 | 同济大学 | 高压实膨润土非饱和渗透系数测量用渗透仪及测量系统 |
CN109613041B (zh) * | 2018-12-11 | 2022-04-22 | 西安理工大学 | 一种基于核磁共振的冻土渗透系数测定系统及方法 |
CN110231359B (zh) * | 2019-07-02 | 2020-07-07 | 吉林大学 | 一种基于地面磁共振弛豫信号的水力特征参数估计方法 |
CN111795990B (zh) * | 2020-01-09 | 2021-12-21 | 中国矿业大学(北京) | 一种基于核磁共振技术测定岩土体中水分迁移过程的方法 |
CN111443027A (zh) * | 2020-06-02 | 2020-07-24 | 东北大学 | 基于磁共振成像技术进行岩石非饱和渗流定量计算的方法 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7567079B2 (en) * | 2007-06-08 | 2009-07-28 | University Of New Brunswick | Methods suitable for measuring capillary pressure and relative permeability curves of porous rocks |
CN105467098B (zh) * | 2015-11-16 | 2017-03-22 | 天津城建大学 | 对基于轴平移技术测试的非饱和土土水特征曲线的修正方法 |
CN105910941B (zh) * | 2016-04-12 | 2018-10-23 | 天津城建大学 | 基于压力板仪的冻土未冻水含量的测试方法 |
CN106353357B (zh) * | 2016-11-08 | 2018-08-03 | 西安理工大学 | 一种渗流作用下砂土介质细观结构变化的监测装置及方法 |
-
2017
- 2017-01-09 CN CN201710012453.XA patent/CN106769771B/zh not_active Expired - Fee Related
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111678853A (zh) * | 2020-06-19 | 2020-09-18 | 中国科学院武汉岩土力学研究所 | 一种测量高压实膨胀土饱和渗透系数的方法 |
CN111678853B (zh) * | 2020-06-19 | 2021-09-14 | 中国科学院武汉岩土力学研究所 | 一种测量高压实膨胀土饱和渗透系数的方法 |
Also Published As
Publication number | Publication date |
---|---|
CN106769771A (zh) | 2017-05-31 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106769771B (zh) | 一种基于低场核磁共振技术的非饱和土渗透系数的测量方法 | |
CN107727679B (zh) | 一种表征深层碳酸盐岩岩石物理学特征方法 | |
CN105588796B (zh) | 一种精确快速测定土壤渗透系数的装置 | |
Vanapalli et al. | Axis translation and negative water column techniques for suction control | |
CN105115873B (zh) | 混凝土抗渗性的无损检测方法 | |
CN109444201A (zh) | 一种测定致密岩心多孔介质中多相流体流动特征的核磁共振实验装置及方法 | |
CN107504924B (zh) | 一种宽域值接触热阻三元测试方法及装置 | |
CN106353357B (zh) | 一种渗流作用下砂土介质细观结构变化的监测装置及方法 | |
Jian-Ping et al. | A QUANTITATIVE EVALUATION METHOD OF LOW PERMEABLE SANDSTONE PORE STRUCTURE BASED ON NUCLEAR MAGNETIC RESONANCE (NMR) LOGGING: A CASE STUDY OF Es4 FORMATION IN THE SOUTH SLOPE OF DONGYING SAG | |
CN110595953B (zh) | 一种页岩混合润湿性的实验测试装置及方法 | |
Lu et al. | Constant flow method for concurrently measuring soil-water characteristic curve and hydraulic conductivity function | |
CN110133035B (zh) | 核磁共振耦合恒速压汞定量表征砂岩储层孔喉结构的方法 | |
CN109781765A (zh) | 一种计算致密储层束缚水液膜厚度的新方法 | |
CN109444020A (zh) | 一种基于3d打印技术的岩体裂隙渗流微观特征观测设备及观测方法 | |
Dong et al. | Measurement of thermal conductivity function of unsaturated soil using a transient water release and imbibition method | |
Alsherif et al. | Measuring the soil water–retention curve under positive and negative matric suction regimes | |
Stagno et al. | Microstructural features assessment of different waterlogged wood species by NMR diffusion validated with complementary techniques | |
CN109856688A (zh) | 基于核磁测井双tw极化增强法的流体性质识别方法 | |
Chen et al. | Investigation on unsaturated hydraulic conductivity with nuclear magnetic resonance technology and instantaneous profile method | |
Malaya et al. | A stud on Wetting-soil water characteristic curve of a sand soil | |
Lin et al. | Soil-water characteristics and permeability of compacted loess considering effects of dry density | |
CN105004639A (zh) | 一种凝胶在岩石微通道中的扩散性的可视化评价方法 | |
PRESSIONS | On the control of low negative water pressures in laboratory tests on unsaturated sand | |
CN110487835A (zh) | 一种计算致密油气藏储层饱和度指数的新方法 | |
Chen et al. | Hydraulic path dependence of shear strength for compacted loess |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20190705 |