CN112098295A - 一种基于核磁共振技术测定无粘性土渗透系数的方法 - Google Patents
一种基于核磁共振技术测定无粘性土渗透系数的方法 Download PDFInfo
- Publication number
- CN112098295A CN112098295A CN202010956548.9A CN202010956548A CN112098295A CN 112098295 A CN112098295 A CN 112098295A CN 202010956548 A CN202010956548 A CN 202010956548A CN 112098295 A CN112098295 A CN 112098295A
- Authority
- CN
- China
- Prior art keywords
- soil
- permeability coefficient
- sample
- magnetic resonance
- nuclear magnetic
- 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.)
- Withdrawn
Links
- 239000002689 soil Substances 0.000 title claims abstract description 105
- 230000035699 permeability Effects 0.000 title claims abstract description 78
- 238000005481 NMR spectroscopy Methods 0.000 title claims abstract description 47
- 238000000034 method Methods 0.000 title claims abstract description 37
- 238000005516 engineering process Methods 0.000 title claims abstract description 22
- 239000002245 particle Substances 0.000 claims abstract description 63
- 238000012360 testing method Methods 0.000 claims abstract description 48
- 239000011148 porous material Substances 0.000 claims abstract description 47
- 238000002360 preparation method Methods 0.000 claims description 27
- 239000004575 stone Substances 0.000 claims description 24
- 229920006395 saturated elastomer Polymers 0.000 claims description 18
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 17
- 230000005415 magnetization Effects 0.000 claims description 9
- 238000004364 calculation method Methods 0.000 claims description 6
- 229920001971 elastomer Polymers 0.000 claims description 4
- 238000012512 characterization method Methods 0.000 claims description 3
- 230000005484 gravity Effects 0.000 claims description 3
- 239000011800 void material Substances 0.000 claims description 3
- 230000014509 gene expression Effects 0.000 claims description 2
- VYPSYNLAJGMNEJ-UHFFFAOYSA-N Silicium dioxide Chemical compound O=[Si]=O VYPSYNLAJGMNEJ-UHFFFAOYSA-N 0.000 description 10
- 239000011324 bead Substances 0.000 description 10
- 239000011521 glass Substances 0.000 description 10
- 239000006004 Quartz sand Substances 0.000 description 9
- 239000012530 fluid Substances 0.000 description 5
- 238000004458 analytical method Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 239000000203 mixture Substances 0.000 description 3
- 238000012986 modification Methods 0.000 description 3
- 230000004048 modification Effects 0.000 description 3
- 230000000903 blocking effect Effects 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 239000004927 clay Substances 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011005 laboratory method Methods 0.000 description 1
- 238000002156 mixing Methods 0.000 description 1
- 230000035515 penetration Effects 0.000 description 1
- 230000000704 physical effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 239000004576 sand Substances 0.000 description 1
- 238000009738 saturating Methods 0.000 description 1
- 238000007790 scraping Methods 0.000 description 1
- 238000007789 sealing Methods 0.000 description 1
- 238000003756 stirring Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000010998 test method Methods 0.000 description 1
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
-
- 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 Health & Medical Sciences (AREA)
- Analytical Chemistry (AREA)
- Pathology (AREA)
- Immunology (AREA)
- High Energy & Nuclear Physics (AREA)
- Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Biochemistry (AREA)
- Geochemistry & Mineralogy (AREA)
- Engineering & Computer Science (AREA)
- Environmental & Geological Engineering (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geology (AREA)
- Dispersion Chemistry (AREA)
- Sampling And Sample Adjustment (AREA)
Abstract
本发明公开了一种基于核磁共振技术测定无粘性土渗透系数的方法，通过核磁共振试验，获得不同颗粒形状的无粘性土在同一孔隙率下的横向弛豫时间T2分布曲线，首先利用T2分布曲线表征土体的内部孔隙结构，再将横向弛豫时间T2分布曲线上的峰值点对应的横向弛豫时间T2peak与土颗粒的平均粒径进行线性拟合，从而基于经验公式推导出横向弛豫时间T2分布曲线上的峰值点对应的横向弛豫时间T2peak与渗透系数的拟合公式。本发明利用核磁共振技术能够快捷、高效、科学地表征无粘性土的渗透系数，对试验土样的需求量小，大幅度减小了试验成本。
Description
技术领域
本发明涉及一种基于核磁共振技术测定无粘性土渗透系数的方法，属于地下工程技术领域。
背景技术
近年来，随着科学技术快速地发展，研究土体内部孔隙结构的方法也日益丰富。土体的内部孔隙结构影响着其物理力学性质，特别地，土体的内部孔隙结构也是影响土体渗透率的重要因素之一，土颗粒的直径与形状又同时影响着土体内部的孔隙结构。学者们也一直致力于探索修正经验公式中颗粒形状对土体渗透率的影响系数，但利用土体内部的孔隙等效半径表征土体渗透率的研究尚不明确。
室内试验需用重塑土测定土体的渗透系数，但重塑土的传统制样要求较高，操作难度较大，最终可能导致实验结果与实际值有较大的偏差。室内常用的测试土体渗透系数的方法为：“常水头法”和“变水头法”。“常水头法”适用于测定渗透性较大的砂性土的渗透系数，但试验用土量较大，在制样过程中难以控制试样的密实度和孔隙比，试验过程的繁琐步骤也会带来较大的试验误差；“变水头法”则适用于测定渗透系数较小的粘性土。借助“常水头法”和“变水头法”所得到的试验结果是一个趋近值，且精确度不高，具有不稳定性。利用三轴仪也可测试土体的渗透系数，但试验过程需要花费的时间较长，试验设备成本较高。上述常规实验室测定土体渗透系数的方法，均无法直观表征土体颗粒形状、孔隙特征和渗透系数之间的关系。
因此，亟需一种既能建立颗粒形状、孔隙等效半径和渗透系数之间的关系，又能快捷、高效、科学地测定无粘性土渗透系数的渗透试验方法。
发明内容
针对上述现有技术的不足，本发明提供了一种基于核磁共振技术测定无粘性土渗透系数的方法，利用核磁共振的横向弛豫时间T2建立颗粒形状、孔隙等效半径与渗透系数之间的联系，进一步地用T2分布曲线峰值点所对应的横向弛豫时间T2peak表征无粘性土试样的渗透系数。
本发明的目的是通过以下技术方案来实现的：一种基于核磁共振技术测定无粘性土渗透系数的方法，该方法包括以下步骤：
1)按预设的孔隙比制备无粘性土试样，试样的饱和度需达到97％以上，试样制备完成后进行核磁共振试验，得到横向弛豫时间T2分布曲线；
2)根据试验用土建立土样孔隙模型(土样孔隙模型根据具体试验的土样颗粒特征建立)，利用土样孔隙模型结合核磁共振原理得到孔隙等效半径R与横向弛豫时间T2的关系式R＝aρ2T2，其中，a为孔隙形状相关系数，ρ2为土颗粒表面的T2表面弛豫强度；
3)对步骤1)得到的T2分布曲线进行线性拟合得到拟合曲线，将拟合曲线和关系式R＝aρ2T2相结合，得到平均粒径D与T2分布曲线峰值点所对应的横向弛豫时间T2peak(峰值横向弛豫时间T2peak)的关系式D＝f(T2peak)；
4)根据Kozeny-Carman渗透系数计算公式得到平均粒径D与渗透系数k的关系式k＝g(D)，再将D＝f(T2peak)代入k＝g(D)，即可得到T2分布曲线峰值点所对应的横向弛豫时间T2peak与渗透系数k的关系式k＝h(T2peak)；
5)将待测无粘性土进行核磁共振试验得到T2分布曲线，带入关系式k＝h(T2peak)中，从而得到待测无粘性土的渗透系数k。
进一步地，所述的步骤1)具体为：
(1.1)制备试样时采用湿法制样，根据预设孔隙比(例如可以取孔隙比e＝0.7)，取适当质量(由选用饱和器的容积计算得到)且粒径均匀的试验土体在饱和器中分若干层逐层击实(例如可以设定分5层)；
(1.2)试样制备完成后，需要在真空饱和缸中对试样进行真空饱和，饱和度需达到97％以上，未达到饱和度要求则继续饱和；
(1.3)试样达到饱和要求后取出，擦干表面水分，放入核磁共振分析仪中进行扫描试验，导出试验数据，记录横向弛豫时间T2分布曲线。
进一步地，所述的步骤1)中，试样的饱和器为单独设计的一种适用于核磁共振试验的一体化饱和试样制备装置，包括装置主体结构、顶部透水石、顶盖、底部透水石和底盖；所述装置主体结构呈圆柱筒状，在两端设有搭接口，所述顶部透水石、底部透水石分别通过搭接口设置在装置主体结构内部的两端处，所述顶盖、底盖分别固定在装置主体结构的两端，且盖面在固定后分别与两端的透水石紧密接触；所述顶盖、底盖的盖壁上的内部圈角处设有橡胶圈，与透水石的外边沿接触，起到阻隔水分向外渗流的作用。
进一步地，所述的步骤2)具体为：
(2.1)基于核磁共振原理，横向弛豫时间T2与横向磁化矢量的表征关系式为：
式中：Mxy(t)为在时间t的回波串的幅度，即横向磁化矢量的幅度；M0为最大横向磁化矢量，与温度、主磁场强和质子含量有关；
(2.2)通过对核磁共振中FID信号的傅里叶转换，土样中孔隙水的T2分布可以表示为：
(2.3)将建立的土样孔隙模型代入式(2)，得到：
R＝aρ2T2 (3)
进一步地，所述的步骤2)中，当孔隙为柱形时，a＝2，当孔隙为球形时，a＝3；试验选取表面弛豫强度ρ2为10μm/s。
进一步地，所述的步骤4)具体为：
(4.1)引入水温在20℃时的Kozeny-Carman渗透系数计算公式：
式中：k为渗透系数，cm/s；γ为水的重度，N/m3；μ为水的黏度，Pa·s；Cf为形状系数；
为比表面积，mm-1；Deff为有效粒径；e为孔隙比；
(4.2)由式(4)推算出不同颗粒形状的无粘性土的平均粒径D与渗透系数k的关系式为：
k＝g(D) (5)
式中：g是关于平均粒径D与渗透系数k的函数；
(4.3)将D＝f(T2peak)代入式(5)可得到横向弛豫时间T2与渗透系数k的关系式：
k＝h(T2peak) (6)
式中：h是关于T2分布曲线峰值点所对应的横向弛豫时间T2peak与渗透系数k的函数
进一步地，所述的步骤(4.1)中，对于圆形土颗粒，Cf＝5；对于角粒状土颗粒，Cf＝14；γ/μ＝9.93×104(cm·s)-1；Deff为有效粒径，采用同一粒径级颗粒时可取平均粒径。
本发明的有益效果是：
1、本发明建立在明确的理论基础之上，采用核磁共振技术对无粘性土试样进行扫描试验，利用横向弛豫时间T2建立起颗粒形状、孔隙等效半径和渗透系数之间的联系，进一步地用T2分布曲线峰值点所对应的横向弛豫时间T2peak表征无粘性土试样的渗透系数。
2、本发明采用核磁共振技术通过对试样孔隙内的流体量、流体性质，以及流体在试样孔隙中的核磁共振特征等测试分析，基于已有理论和经验公式，快速获得孔隙度、可动流体饱和度等物性和流体参数。在对土体内部微观孔隙测定和分析研究的同时，还可以根据峰值横向弛豫时间T2peak来表征土体的渗透特性，测试速度快，可快速完成多个试样的测试工作，测试结果的准确度相对较高。
3、本发明在制样过程中采用了一种适用于核磁共振试验的一体化饱和试样制备装置，此装置密封性好且具有透水饱和度高的特点，饱和后擦干表面水分即可直接放入核磁共振分析仪中进行扫描试验。另外，此装置简易便捷，不仅能够提高试验效率，还能够在很大程度上减少人为操作对试样的扰动。
4、相比传统测定渗透系数的方法，本发明利用核磁共振技术能够快捷、高效、科学地表征无粘性土的渗透系数，对试验土样的需求量小，大幅度减小了试验成本。另外，简易的操作和较高的准确度也使得核磁共振技术在测定土体渗透特性方面具有广阔的发展前景。
5、本发明与现有方法进行了比较，结果显示：本发明提出的方法与现有常用的方法所测定的数据结果处在一个合理的误差阈值内，可运用于实际工程中。
附图说明
图1为本发明的实例流程图；
图2为土样孔隙模型图；
图3为玻璃珠(0.800～1.000mm)的T2分布曲线；
图4为石英砂(0.850～1.250mm)的T2分布曲线；
图5为T2峰顶点时间与玻璃珠的平均粒径的关系曲线；
图6为T2峰顶点时间与石英砂的平均粒径的关系曲线；
图7为饱和试样制备装置主体部分剖面图；
图8为饱和试样制备装置顶盖和底盖剖面图；
图9为饱和试样制备装置顶盖俯视图；
图10为饱和试样制备装置整体侧视图；
其中：图3～图6为核磁共振试验数据整理图；图7～图10为饱和试样制备装置示意图；图中：装置主体结构1、搭接口1-1、透水石2、土样3、顶盖4-1、底盖4-2、橡胶圈5、盖面6、盖壁7。
具体实施方式
下面结合附图和具体实施例对本发明作进一步详细说明。
如图1所示的实施例流程图，本发明提出的一种基于核磁共振技术测定无粘性土渗透系数的方法，包括如下步骤：
1)按预设的孔隙比制备无粘性土试样，试样的饱和度需达到97％以上，试样制备完成后进行核磁共振试验，得到横向弛豫时间T2分布曲线；
2)根据试验用土建立土样孔隙模型(土样孔隙模型根据具体试验的土样颗粒特征建立)，利用土样孔隙模型结合核磁共振原理得到孔隙等效半径R与横向弛豫时间T2的关系式R＝aρ2T2，其中，a为孔隙形状相关系数，ρ2为土颗粒表面的T2表面弛豫强度；
3)对步骤1)得到的T2分布曲线进行线性拟合得到拟合曲线，将拟合曲线和关系式R＝aρ2T2相结合，得到平均粒径D与T2分布曲线峰值点所对应的横向弛豫时间T2peak的关系式D＝f(T2peak)，其中，f是关于平均粒径D与T2分布曲线峰值点所对应的横向弛豫时间T2peak的函数；
4)根据Kozeny-Carman渗透系数计算公式得到平均粒径D与渗透系数k的关系式k＝g(D)，再将D＝f(T2peak)代入k＝g(D)，即可得到T2分布曲线峰值点所对应的横向弛豫时间T2peak与渗透系数k的关系式k＝h(T2peak)；
5)将待测无粘性土进行核磁共振试验得到T2分布曲线，带入关系式k＝h(T2peak)中，从而得到待测无粘性土的渗透系数k。
关于步骤1)的具体实施说明如下：
(1.1)将颗粒均匀的土样(本实例所用试验土样如表1所示)加适量的无气水充分混合后，平均分五次放入如图7～图10所示的饱和试样制备装置内，且逐层击实至相应的高度，前4层在击实后在土样的表面用刮刀刮花，以防试样内部结构分层；
表1试验材料物理特征及预设孔隙比
(1.2)制样完成后，将饱和试样制备装置的顶盖、底盖旋紧，放入真空饱和缸中饱和至97％饱和度，未达到饱和度要求时则继续饱和；
(1.3)将饱和完成的试样取出后擦干表面水分，放入调试好的核磁共振分析仪中进行扫描试验，可得到如图3、图4所示的横向弛豫时间T2分布曲线图，且对试验数据进行记录和整理。
关于步骤2)的具体实施说明如下：
(2.1)本实例采用的试验土样为玻璃珠和石英砂，故根据玻璃珠和石英砂的特点建立如图2所示的土样孔隙模型；
(2.2)基于核磁共振原理，横向弛豫时间T2与横向磁化矢量的表征关系式为：
式中：Mxy(t)为在时间t的回波串的幅度，即横向磁化矢量的幅度；M0为最大横向磁化矢量，与温度、主磁场强和质子含量有关；
(2.3)通过对核磁共振中FID信号的傅里叶转换，土样中孔隙水的T2分布可以表示为：
(2.4)将建立的土样孔隙模型代入式(2)，得到：
R＝aρ2T2 (3)
式中：a为孔隙形状相关系数，当孔隙为柱形时，a＝2；当孔隙为球形时，a＝3。则有：
Rb＝3ρ2T2 (4)
Rs＝2ρ2T2 (5)
式中：Rb为玻璃珠的孔隙等效半径；Rs为石英砂的孔隙等效半径。
关于步骤3)的具体实施说明如下：
(3.1)对步骤1)得到的T2分布曲线进行线性拟合得到拟合曲线，如图5、图6所示，通过线性分析并结合上述式(3)，将平均粒径D与T2分布曲线峰值点所对应的横向弛豫时间T2peak分别拟合在公式中：
D＝f(T2peak) (6)
式中：f是关于平均粒径D与T2分布曲线峰值点所对应的横向弛豫时间T2peak的函数。
(3.2)玻璃珠和石英砂的平均粒径D与T2分布曲线峰值点所对应的横向弛豫时间T2peak的具体拟合结果如下：
Db＝8×10-7T2 2-4.0×10-4T2+2.03×10-1 (7)
Ds＝5×10-7T2 2+5.0×10-4T2-8.61×10-2 (8)
式中：Db、Ds分别为玻璃珠和石英砂的平均粒径，单位mm。
关于步骤4)的具体实施说明如下：
(4.1)引入水温在20℃时的Kozeny-Carman渗透系数计算公式：
式中：k为渗透系数，cm/s；γ为水的重度，N/m3；μ为水的黏度，Pa·s；γ/μ＝9.93×104(cm·s)-1；Cf为形状系数，对于圆形土颗粒，Cf＝5，对于角粒状土颗粒，Cf＝14；
为比表面积，mm-1；其中，Deff为有效粒径，采用同一粒径级颗粒时可取平均粒径；e为孔隙比，设孔隙比函数为f(e)，令f(e)＝e3/(1+e)，本实例中预设孔隙比e＝0.7，则f(e)＝0.2018。
(4.2)由式(9)推算出不同颗粒形状的无粘性土的平均粒径D与渗透系数k的关系式为：
k＝g(D) (10)
式中：g是关于平均粒径D与渗透系数k的函数。
结合本实例中的具体参数，可推算玻璃珠和石英砂这两种颗粒形状的砂土的平均粒径与渗透系数的关系式：
当土颗粒为圆形颗粒时：kb＝1.1133Db 2 (11)
当土颗粒为角粒状颗粒时：ks＝0.3976Ds 2 (12)
式中：kb为玻璃珠的渗透系数；ks为石英砂的渗透系数。
(4.3)将上述式(6)代入式(10)中可得到T2分布曲线峰值点所对应的横向弛豫时间T2peak与渗透系数k的关系式：
k＝h(T2peak) (13)
式中：h是关于T2分布曲线峰值点所对应的横向弛豫时间T2peak与渗透系数k的函数。
将上述式(7)、式(8)代入式(11)、式(12)中，分别得到玻璃珠和石英砂的渗透系数k与T2分布曲线峰值点所对应的横向弛豫时间T2peak的具体关系式：
关于步骤5)的实施说明如下：
通过上述各具体步骤得到的式(14)和式(15)即为基于核磁共振技术测定无粘性土的渗透系数的最终关系式，从表2中可以看到本发明提出的方法与Terzaghi等、朱崇辉等提出的经验公式计算的渗透系数值的对比。通过对比和分析后发现：Terzaghi公式并未考虑颗粒形状与孔隙结构对渗透系数的影响；朱崇辉公式反映了渗透系数与颗粒组成的相互关系，适用于颗粒形状不一和孔隙结构复杂的砂土；而核磁共振技术在一定程度上弥补了经验公式的不足，同时也已通过试验验证了该技术的合理性和准确性，扩展了该技术的适用土样的范围。
表2各粒径级的渗透系数
注：表中kTerzaghi/k和k朱崇辉/k为比值关系，与“1”作比较。
关于步骤1)中的饱和试样制备装置具体实施说明如下：
饱和试样制备装置为单独设计的一种适用于核磁共振试验的一体化饱和试样制备装置，如图7-10所示，包括装置主体结构1、顶部透水石、顶盖4-1、底部透水石和底盖4-2；装置主体结构1呈圆柱筒状，在两端设有搭接口1-1，顶部和底部的透水石2分别通过搭接口1-1设置在装置主体结构1内部的两端处，顶盖4-1、底盖4-2分别固定在装置主体结构1的两端，且盖面6在固定后分别与两端的透水石2紧密接触；顶盖4-1、底盖4-2的盖壁7上的内部圈角处设有橡胶圈5，与透水石2的外边沿接触，起到阻隔水分向外渗流的作用。
装置主体结构1采用PVC材料，内外表面均平整光滑；顶部和底部的透水石2的形状大小完全一致，其直径大于装置主体结构内两端搭接口之间的圆柱状空间的直径，其外径与搭接口处所形成的圆截面直径一致。
饱和试样制备装置的具体操作过程如下：
1)将底部透水石装入装置主体结构的相应位置，制样开始；
2)将添加适量无气水的土样3搅拌均匀后分5次加入装置内，并逐层击实至相应的高度；
3)制样完成后，加盖顶部透水石，再分别在顶部和底部透水石的两端施加夹具进行固定；
4)固定好后，将整体结构放入真空饱和缸中进行饱和；
5)饱和工作完成后，将整体结构取出，擦拭表面水分，小心去掉夹具后，立即加盖顶盖和底盖，即完成试样的制备；
6)试样制备完成后，放入核磁共振分析仪内，随即开展核磁共振试验。
以上所述仅是本发明的优选实施方式，虽然本发明已以较佳实施例披露如上，然而并非用以限定本发明。任何熟悉本领域的技术人员在不脱离本发明技术方案范围情况下，都可利用上述揭示的方法和技术内容对本发明技术方案做出许多可能的变动和修饰，或修改为等同变化的等效实施例。因此，凡是未脱离本发明技术方案的内容，依据本发明的技术实质对以上实施例所做的任何的简单修改、等同变化及修饰，均仍属于本发明技术方案保护的范围内。
Claims (7)
1.一种基于核磁共振技术测定无粘性土渗透系数的方法，其特征在于，包括如下步骤：
1)按预设的孔隙比制备无粘性土试样，试样的饱和度需达到97％以上，试样制备完成后进行核磁共振试验，得到横向弛豫时间T2分布曲线；
2)根据试验用土建立土样孔隙模型，利用土样孔隙模型结合核磁共振原理得到孔隙等效半径R与横向弛豫时间T2的关系式R＝aρ2T2，其中，a为孔隙形状相关系数，ρ2为土颗粒表面的T2表面弛豫强度；
3)对步骤1)得到的T2分布曲线进行线性拟合得到拟合曲线，将拟合曲线和关系式R＝aρ2T2相结合，得到平均粒径D与横向弛豫时间T2分布曲线上的峰值点对应的横向弛豫时间T2peak的关系式D＝f(T2peak)；
4)根据Kozeny-Carman渗透系数计算公式得到平均粒径D与渗透系数k的关系式k＝g(D)，再将D＝f(T2peak)代入k＝g(D)，即可得到横向弛豫时间T2分布曲线上的峰值点对应的横向弛豫时间T2peak与渗透系数k的关系式k＝h(T2peak)；
5)将待测无粘性土进行核磁共振试验得到T2分布曲线，带入关系式k＝h(T2peak)中，从而得到待测无粘性土的渗透系数k。
2.根据权利要求1所述的一种基于核磁共振技术测定无粘性土渗透系数的方法，其特征在于，所述的步骤1)具体为：
(1.1)制备试样时采用湿法制样，根据预设孔隙比，取粒径均匀的试验土体在饱和器中分若干层逐层击实；
(1.2)试样制备完成后，需要在真空饱和缸中对试样进行真空饱和，饱和度需达到97％以上，未达到饱和度要求则继续饱和；
(1.3)试样达到饱和要求后取出，擦干表面水分，放入核磁共振分析仪中进行扫描试验，导出试验数据，记录横向弛豫时间T2分布曲线。
3.根据权利要求1所述的一种基于核磁共振技术测定无粘性土渗透系数的方法，其特征在于，所述的步骤1)中，试样的饱和通过一体化饱和试样制备装置实现，包括装置主体结构、顶部透水石、顶盖、底部透水石和底盖；所述装置主体结构呈圆柱筒状，在两端设有搭接口，所述顶部透水石、底部透水石分别通过搭接口设置在装置主体结构内部的两端处，所述顶盖、底盖分别固定在装置主体结构的两端，且盖面在固定后分别与两端的透水石紧密接触；所述顶盖、底盖的盖壁上的内部圈角处设有橡胶圈，与透水石的外边沿接触。
4.根据权利要求1所述的一种基于核磁共振技术测定无粘性土渗透系数的方法，其特征在于，所述的步骤2)具体为：
(2.1)基于核磁共振原理，横向弛豫时间T2与横向磁化矢量的表征关系式为：
式中：Mxy(t)为在时间t的回波串的幅度，即横向磁化矢量的幅度；M0为最大横向磁化矢量；
(2.2)通过对核磁共振中FID信号的傅里叶转换，土样中孔隙水的T2分布可以表示为：
(2.3)将建立的土样孔隙模型代入式(2)，得到：
R＝aρ2T2 (3)。
5.根据权利要求1所述的一种基于核磁共振技术测定无粘性土渗透系数的方法，其特征在于，所述的步骤2)中，当孔隙为柱形时，a＝2，当孔隙为球形时，a＝3；试验选取表面弛豫强度ρ2为10μm/s。
6.根据权利要求1所述的一种基于核磁共振技术测定无粘性土渗透系数的方法，其特征在于，所述的步骤4)具体为：
(4.1)引入水温在20℃时的Kozeny-Carman渗透系数计算公式：
式中：k为渗透系数，cm/s；γ为水的重度，N/m3；μ为水的黏度，Pa·s；Cf为形状系数；
为比表面积，mm-1；Deff为有效粒径；e为孔隙比；
(4.2)由式(4)推算出不同颗粒形状的无粘性土的平均粒径D与渗透系数k的关系式为：
k＝g(D) (5)
式中：g是关于平均粒径D与渗透系数k的函数；
(4.3)将D＝f(T2peak)代入式(5)可得到横向弛豫时间T2与渗透系数k的关系式：
k＝h(T2peak) (6)
式中：h是关于横向弛豫时间T2分布曲线上的峰值点对应的横向弛豫时间T2peak与渗透系数k的函数。
7.根据权利要求6所述的一种基于核磁共振技术测定无粘性土渗透系数的方法，其特征在于，所述的步骤(4.1)中，对于圆形土颗粒，Cf＝5；对于角粒状土颗粒，Cf＝14；γ/μ＝9.93×104(cm·s)-1；Deff为有效粒径，采用同一粒径级颗粒时可取平均粒径。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010956548.9A CN112098295A (zh) | 2020-09-11 | 2020-09-11 | 一种基于核磁共振技术测定无粘性土渗透系数的方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010956548.9A CN112098295A (zh) | 2020-09-11 | 2020-09-11 | 一种基于核磁共振技术测定无粘性土渗透系数的方法 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN112098295A true CN112098295A (zh) | 2020-12-18 |
Family
ID=73752048
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010956548.9A Withdrawn CN112098295A (zh) | 2020-09-11 | 2020-09-11 | 一种基于核磁共振技术测定无粘性土渗透系数的方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112098295A (zh) |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112730102A (zh) * | 2020-12-29 | 2021-04-30 | 浙大城市学院 | 一种基于标准差椭圆的颗粒类材料剪切带演化测定方法 |
| CN113075106A (zh) * | 2021-03-26 | 2021-07-06 | 桂林理工大学 | 一种砾石土侵蚀阶段的水力参数及涌水量的确定方法 |
| CN114034621A (zh) * | 2021-11-09 | 2022-02-11 | 中国电建集团成都勘测设计研究院有限公司 | 土体渗透系数的获取方法及装置、模型构建方法 |
| CN114324103A (zh) * | 2021-12-15 | 2022-04-12 | 青岛理工大学 | 一种饱和黏土体渗透系数测定方法 |
| CN114577673A (zh) * | 2022-01-18 | 2022-06-03 | 重庆大学 | 一种等粒径煤粉的孔隙润湿性nmr表征方法 |
-
2020
- 2020-09-11 CN CN202010956548.9A patent/CN112098295A/zh not_active Withdrawn
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112730102A (zh) * | 2020-12-29 | 2021-04-30 | 浙大城市学院 | 一种基于标准差椭圆的颗粒类材料剪切带演化测定方法 |
| CN113075106A (zh) * | 2021-03-26 | 2021-07-06 | 桂林理工大学 | 一种砾石土侵蚀阶段的水力参数及涌水量的确定方法 |
| CN114034621A (zh) * | 2021-11-09 | 2022-02-11 | 中国电建集团成都勘测设计研究院有限公司 | 土体渗透系数的获取方法及装置、模型构建方法 |
| CN114324103A (zh) * | 2021-12-15 | 2022-04-12 | 青岛理工大学 | 一种饱和黏土体渗透系数测定方法 |
| CN114324103B (zh) * | 2021-12-15 | 2023-12-12 | 宿迁学院 | 一种饱和黏土体渗透系数测定方法 |
| CN114577673A (zh) * | 2022-01-18 | 2022-06-03 | 重庆大学 | 一种等粒径煤粉的孔隙润湿性nmr表征方法 |
| CN114577673B (zh) * | 2022-01-18 | 2024-02-02 | 重庆大学 | 一种等粒径煤粉的孔隙润湿性nmr表征方法 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN112098295A (zh) | 一种基于核磁共振技术测定无粘性土渗透系数的方法 | |
| Estabragh et al. | Laboratory investigation of the effect of cyclic wetting and drying on the behaviour of an expansive soil | |
| Puppala et al. | Swell and shrinkage characterizations of unsaturated expansive clays from Texas | |
| Wang et al. | A slurry consolidation approach to reconstitute low-plasticity silt specimens for laboratory triaxial testing | |
| Wang et al. | Characterization of freeze–thaw effects within clay by 3D X-ray computed tomography | |
| CN108645676B (zh) | 基于超声波技术的季冻区路基土材料损伤检测与评价方法 | |
| CN109470580B (zh) | 一种评价我国不同海域粘性土强度参数的方法 | |
| CN110672494B (zh) | 一种多孔混凝土不同孔隙率的快速测定方法 | |
| CN110702541B (zh) | 一种湿化作用下路基土多级动力加载累积变形试验方法 | |
| Toll et al. | Tensiometer techniques for determining soil water retention curves | |
| CN107102066A (zh) | 一种室内超声检测气泡混合轻质土强度的装置及方法 | |
| WO2004013613A1 (en) | Determination of envelope volume and density of porous samples | |
| CN104931392A (zh) | 一种采用吸附测量固体表面自由能的装置及方法 | |
| CN104792974A (zh) | 一种采用全自动真密度仪测试水泥浆体化学收缩的方法 | |
| Milatz et al. | Settlements in unsaturated granular soils induced by changes in saturation and suction | |
| Wang et al. | Soil-water characteristic surface model of soil-rock mixture | |
| CN113654888B (zh) | 一种炭质泥岩永久变形快速预测方法 | |
| Dominguez-Quintans et al. | Quality assessment of a new in-mould slurry deposition method for triaxial specimen reconstitution of clean and silty sands | |
| CN109520905A (zh) | 一种透水混凝土透水非均匀性测定装置及其测试方法 | |
| CN112182872A (zh) | 一种黄土浸水孔隙变化影响下的渗流特性的分析方法 | |
| Musso et al. | Monitoring drying and wetting of a cement bentonite mixture with Electrical Resistivity Tomography | |
| Chen et al. | Effect of compactness degree on the hydraulic properties for coarse soils of high-speed railway embankment | |
| Goodman et al. | Effects of sample thickness and shape on cracking patterns and crack depth in a compacted clay | |
| Malaya et al. | A stud on Wetting-soil water characteristic curve of a sand soil | |
| CN111504874A (zh) | 一种饱和非饱和土水特征参数测试装置及实验方法 |
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 | ||
| WW01 | Invention patent application withdrawn after publication | ||
| WW01 | Invention patent application withdrawn after publication |
Application publication date: 20201218 |