1.一种基于核磁共振曲线的饱和及非饱和土渗透系数预测方法,其特征在于:结合达西定律以及核磁共振曲线提出一种基于核磁共振曲线的饱和及非饱和土体渗透系数预测模型,饱和土体渗透系数和非饱和土体相对渗透系数的预测模型分别为:公式(1)为饱和土体渗透系数预测模型、公式(2)为非饱和土体相对渗透系数预测模
2 2
型;式中:kc=γλρ2 /32μpi,γ为流体重度、λ为沿程阻力系数、ρ2表示表面弛豫率、μ为粘度、pi为第i级孔隙通道实际长度与土样长度l的比值;对于同一土样其值为常数;T2i为核磁共振试验中测得的第i个弛豫时间,利用核磁共振试验测得的弛豫时间将土体内部孔隙按孔径大小划分为n级,θi为第i级孔隙通道充满水分的体积含水量,假设总孔隙通道有n级,现只有第1~m级孔隙充满水,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)中,g为重力加速度,λ为沿程阻力系数,d为连通孔隙的等效孔径,一般土中水的流速较慢,常处于层流状态,其沿程阻力系数表示为公式(5):公式(5)中,Re为雷诺数,其表达式(6)为:
表达式(6)中,ρ为流体的密度,μ为粘度;
由公式(4)、公式(5)及表达式(6)可得公式(7):公式(7)中,γ为流体重度;
局部水头损失按公式(8)计算:
公式(8)中,g为重力加速度,ζ为局部阻力系数,由实验确定;
比较公式(7)和公式(8),公式(7)表明沿程水头损失与渗透速度的一次方成正比,公式(8)表明局部水头损失与渗透速度的二次方成正比,故在计算中可忽略局部水头损失带来的影响,即可认为hw=hf,则J=hf/l,即公式(3)可得:由公式(7)、公式(9)两式得:
将组成土体的所有连通孔隙通道的渗透系数k叠加起来即可得到理论饱和渗透系数为:公式(11)中,di为第i级孔隙通道的等效直径,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)
公式(14)、公式(15)中,T2表示横向弛豫时间,ρ2表示表面弛豫率,S/V表示单位孔隙的表面积S与单位岩石中孔隙体积V的比值;
结合公式(13)、公式(15)两式可得:
2 2
令kc=γλρ2/32μpi,由公式(16)可得:
3.根据权利要求2所述基于核磁共振曲线的饱和及非饱和土渗透系数预测方法,其特征在于:用核磁共振试验对非饱和土体相对渗透系数进行预测,非饱和土体渗透系数预测模型公式(2)具体推导步骤为:假设总孔隙通道有n级,现只有第1~m级孔隙充满水,m
公式(17)可简化为:
联立公式(15)和公式(18)知:
由饱和土体渗透系数和非饱和土体渗透系数的关系,得到非饱和土体渗透系数:kw=kr×ks (19)。