一种基于核磁共振曲线的饱和及非饱和土渗透系数预测方法技术

本发明专利技术提供一种基于核磁共振曲线的饱和土和非饱和土渗透系数预测方法。结合渗透试验和核磁共振实验对饱和土渗透比例常数kc进行标定，实现了利用饱和土的T2分布曲线预测饱和土渗透系数；对于不同含水量状态土体的非饱和渗透系数的预测，需要对这些含水量状态下的土体分别进行核磁共振实验，利用某一状态下土样的T2分布曲线预测该状态下的相对渗透系数，进而得到非饱和渗透系数的预测值。此外，本发明专利技术还提出了T2分布曲线出现缺级现象时的处理方法，提高了预测精度。本发明专利技术所述方案能够克服试验直接测定饱和土体及非饱和土体渗透系数时产生的人力和物力以及大量时间的消耗，克服了现有的渗透系数预测方法以经验法为主且预测精度较低的缺点。

【技术实现步骤摘要】
一种基于核磁共振曲线的饱和及非饱和土渗透系数预测方法
本专利技术属于岩土工程中饱和土及非饱和土渗透系数研究的
，具体涉及一种基于核磁共振曲线的饱和及非饱和土渗透系数预测方法。
技术介绍
非饱和土体中存在封闭气泡，这使得水在土中渗透过程中需要先绕过气泡，增大了土中水流程的绕曲度，这使得非饱和土中的渗流非常缓慢，不易测量。并且随着孔隙比和饱和度的不同，非饱和土渗透系数变化剧烈。由于试验测量非饱和土渗透系数的方法费时费力，目前，研究者们从微观孔隙的角度或通过分析土水特征曲线与相对渗透系数的关系提出了许多间接预测的模型，但这些模型的预测精度有待提高。通过定量研究建立更为准确的饱和土及非饱和土渗透系数的模型具有重要的研究意义。
技术实现思路
本专利技术的目的就是针对上述技术的不足，提供一种基于核磁共振曲线的饱和及非饱和土渗透系数预测方法，克服了现有的渗透系数预测方法以经验法为主且预测精度较低的缺点，实现利用核磁共振试验对任意干密度土样的饱和土渗透系数和任意含水量的非饱和土渗透系数进行预测。为实现上述目的，本专利技术所设计的基于核磁共振曲线的饱和及非饱和土渗透系数预测方法，结合达西定律以及核磁共振曲线提出一种基于核磁共振曲线的饱和及非饱和土体渗透系数预测模型，饱和土体渗透系数和非饱和土体相对渗透系数的预测模型分别为：公式(1)为饱和土体渗透系数预测模型、公式(2)为非饱和土体相对渗透系数预测模型；式中：kc为综合渗透比例常数，对于同一土样其值为常数；T2i为核磁共振试验中测得的第i个弛豫时间，利用核磁共振试验测得的弛豫时间将土体内部孔隙按孔径大小划分为n级，θi为第i级孔隙通道充满水分的体积含水量，m为第m级孔隙通道且m＜n；由饱和土体渗透系数和非饱和土体相对渗透系数的关系，得到非饱和土体渗透系数预测模型为：kw＝kr×ks(19)；利用饱和土体渗透系数试验测得某一干密度土体的饱和土体渗透系数ks实测,将实测的饱和土体渗透系数ks实测代入公式(1)即得到该土样在该干密度下渗透比例常数kc0；对于不同干密度的同一土样，其渗透比例常数kc近似相等，因此，将kc0替换公式(1)中的kc，从而计算得到任意干密度土样的核磁共振曲线预测其饱和土体渗透系数；对非饱和土体进行核磁共振试验测得其核磁曲线(T2i和θi的关系曲线)，然后将T2i和θi代入公式(2)中对其非饱和土体相对渗透系数进行预测。进一步地，饱和土体渗透系数预测模型公式(1)具体推导步骤为：1)利用微观孔隙通道，水分在土中的渗流满足达西定律公式(3)：公式(3)中，q为单位渗水量；A为试样截面面积；v为渗透速度；k为土的渗透系数；J为水力梯度，J＝hw/l，l为试样的长度，hw为水头损失，hw＝hf+hj，hf为沿程损失，hj为局部水头损失；根据达西提出的沿程水头损失的通用表达式(4)：表达式(4)中，λ为沿程阻力系数，d为连通孔隙的等效孔径，一般土中水的流速较慢，常处于层流状态，其沿程阻力系数表示为公式(5)：公式(5)中，Re为雷诺数，其表达式(6)为：表达式(6)中，ρ为流体的密度，μ为粘度；由公式(4)、公式(5)及表达式(6)可得公式(7)：局部水头损失按公式(8)计算：公式(8)中，ζ为局部阻力系数，由实验确定；比较公式(7)和公式(8)，公式(7)表明沿程水头损失与渗透速度的一次方成正比，公式(8)表明局部水头损失与渗透速度的二次方成正比，故在计算中可忽略局部水头损失带来的影响，即可认为hw＝hf，则J＝hf/l，即公式(3)可得：由公式(7)、公式(9)两式得：将组成土体的所有连通孔隙通道的渗透系数k叠加起来即可得到理论饱和渗透系数为：公式(11)中，A为试样的横截面积，Ai为孔隙通道的横截面积；假设第i级孔隙通道实际长度与土样长度l比值为pi，则实际长度为pil，若相应第i级孔隙通道总体积为θiVT，其中θi为该级孔隙通道充满水分的体积含水量；由核磁共振试验测得θi的值；VT为测试土样的总体积，则相应通道的横截面面积为：假设pi为常数，将公式(12)代入公式(11)得到饱和土体渗透系数预测模型ks：2)由Korringa、Seevers和Torrey于1962年提出的KST模型的理论可知土体孔隙结构与弛豫时间之间的关系表示为：将公式(14)整理得：d＝λρ2T2(15)结合公式(13)、公式(15)两式可得：令kc＝γλ2ρ22/32μpi，由公式(16)可得：进一步地，用核磁共振试验对非饱和土体相对渗透系数进行预测，非饱和土体渗透系数预测模型公式(2)具体推导步骤为：假设总孔隙通道有n级，现只有第1～m级孔隙充满水(m<n)。由《基于微观孔隙通道的饱和/非饱和土渗透系数模型及其应用》(水利学报,2017,48(06):702-709)非饱和土体相对渗透系数表示为：公式(17)中，θi为第i级孔隙充满水的体积含水量，ψi为相应于第i级孔隙孔径di的基质吸力，ψi＝4Tscosα/di，其中，温度不变时Tscosα为常数；公式(17)可简化为：联立公式(15)和公式(18)知：由饱和土体渗透系数和非饱和土体渗透系数的关系，得到非饱和土体渗透系数：kw＝kr×ks(19)。值得注意的是：针对土体中出现孔隙缺级的情况，利用本专利技术中提出的预测模型计算非饱和土相对渗透系数时应将缺级段后的孔隙舍弃。其原因为：通常情况下，孔径相近的孔隙(第i级及邻级孔隙)相互连接形成等效孔径为di的孔隙通道，其渗透特性符合达西定律。对于缺级段左右两侧的孔隙，左侧小孔隙和右侧大孔隙的孔径相差巨大(dl>>ds)。当土体中大孔隙左侧出现缺级时，大孔隙数量也极少，不足以形成等效孔径为dl的孔隙通道，通常大孔隙与小孔隙连接形成孔隙通道，其等效孔径受小孔隙孔径ds的限制。除此之外，由于形成通道的大、小孔隙的孔径相差很大，通道内墨水瓶效应明显，阻碍土中水力传导的进行。因此，将这一类大孔隙不能形成等效孔径为dl的通道，对渗透系数不起控制作用，可近似舍去。本专利技术基于核磁共振曲线的饱和及非饱和土渗透系数预测方法，一是针对某一干密度不同含水量的试样，本专利技术方法需要对各含水量状态下的试样进行核磁共振试验，然后利用T2分布预测其相对渗透系数，而已有的方法只需对这一干密度下的饱和土体进行核磁共振试验，然后利用归一化模型预测不同含水量时土体的相对渗透系数；二是本专利技术提出了土体中出现的缺级现象时的预测方法，考虑了孔隙通道中的墨水瓶效应对土体渗透特性的影响，相比之下，本专利技术提出的非饱和土相对渗透系数的预测方法更能反映土中渗流的真实状态，因此预测效果更准确。与现有技术相比，本专利技术具有以下优点：通过本专利技术所述方案能够克服通过试验直接测定饱和土体及非饱和土体渗透系数时产生的人力和物力以及大量时间的消耗。通过对所需状态下土体进行核磁共振试验，即可预测该状态下土样的饱和土体渗透系数或该状态下非饱和土体的渗透系数，进而可以大幅度减小试验工作量，克服了现有的渗透系数预测方法以经验法为主且预测精度较低的缺点，实现利用核磁共振试验对任意干密度土样的饱和土渗透系数和任意含水量的非饱和土渗透系数进行预测。附图说明图1为通过核磁曲线划分土体孔隙的示意图；图2为核磁共振曲线中不连续孔隙的区分；图3为本专利技术具体实本文档来自技高网...

【技术保护点】
1.一种基于核磁共振曲线的饱和及非饱和土渗透系数预测方法，其特征在于：结合达西定律以及核磁共振曲线提出一种基于核磁共振曲线的饱和及非饱和土体渗透系数预测模型，饱和土体渗透系数和非饱和土体相对渗透系数的预测模型分别为：

【技术特征摘要】
1.一种基于核磁共振曲线的饱和及非饱和土渗透系数预测方法，其特征在于：结合达西定律以及核磁共振曲线提出一种基于核磁共振曲线的饱和及非饱和土体渗透系数预测模型，饱和土体渗透系数和非饱和土体相对渗透系数的预测模型分别为：公式(1)为饱和土体渗透系数预测模型、公式(2)为非饱和土体相对渗透系数预测模型；式中：kc为综合渗透比例常数，对于同一土样其值为常数；T2i为核磁共振试验中测得的第i个弛豫时间，利用核磁共振试验测得的弛豫时间将土体内部孔隙按孔径大小划分为n级，θi为第i级孔隙通道充满水分的体积含水量，m为第m级孔隙通道且m＜n；由饱和土体渗透系数和非饱和土体相对渗透系数的关系，得到非饱和土体渗透系数预测模型为：kw＝kr×ks(19)；利用饱和土体渗透系数试验测得某一干密度土体的饱和土体渗透系数ks实测,将实测的饱和土体渗透系数ks实测代入公式(1)即得到该土样在该干密度下渗透比例系数kc0；对于不同干密度的同一土样，其渗透比例系数kc近似相等，因此，将kc0替换公式(1)中的kc，从而计算得到任意干密度土样的核磁共振曲线预测其饱和土体渗透系数；对非饱和土体进行核磁共振试验测得其核磁曲线，然后将T2i和θi代入公式(2)中对其非饱和土体相对渗透系数进行预测。2.根据权利要求1所述基于核磁共振曲线的饱和及非饱和土渗透系数预测方法，其特征在于：饱和土体渗透系数预测模型公式(1)具体推导步骤为：1)利用微观孔隙通道，水分在土中的渗流满足达西定律公式(3)：公式(3)中，q为单位渗水量；A为试样截面面积；v为渗透速度；k为土的渗透系数；J为水力梯度，J＝hw/l，l为试样的长度，hw为水头损失，hw＝hf+hj，hf为沿程损失，hj为局部水头损失；根据达西提出的沿程水头损失的通用表达式(4)：表达式(4)中，λ为沿程阻力系数，d为连通孔隙的等效孔径，一般土中水的流速较慢，常处于层流状态，其沿程阻力系数表示为公式(5)：公式(5)中，Re为雷诺数，其表达式(6)为：表达式(6)中，ρ为流体...

【专利技术属性】
技术研发人员：陶高梁，雷达，庄心善，胡其志，梅利芳，徐维生，赵蔚，李梓月，李奕，匡航，
申请(专利权)人：湖北工业大学，
类型：发明
国别省市：湖北,42