[发明专利]黏土颗粒间孔隙等效孔径的方法

权利要求书

1.一种黏土颗粒间孔隙等效孔径的方法，其特征在于，具体按照以下步骤进行：

步骤1，

首先制备含黏粒土样进行土-水特征实验，并通过MATLAB用(1)式拟合得到完整的土- 水特征曲线，以质量含水量表示，拟合公式如下：

θ=(θs-θr)[11+an(ua-uw)n]m+θr;---(1)]]>

式中，u_a-u_w为基质吸力；θ_r为残余含水率；θ_s饱和含水率；a、m及n均为拟合参数，参数n与 土体孔径分布有关，参数m与土-水特征曲线的整体对称性有关；

步骤2，

以(u_a-u_w)_i＝(u_a-u_w)_i-1+50(i-1)，其中i＝2，3，4……，(u_a-u_w)₁＝100kPa的方式选取基 质吸力于144000kPa～0kPa之间的数据，通过(4)式将基质吸力转化为相对湿度RH：

RH=exp(-(ua-uw)vwRT);---(4)]]>

由以上推导可得，毛细半径r以相对湿度表示形式如下：

r=-2TsvwcosαRT ln(RH);---(5a)]]>

或表示为基质吸力形式：

r=2Tscosαua-uw;---(5b)]]>

式中，v_w为水蒸气的偏摩尔体积；R为通用气体常数；T为热力学温度；T_s为毛细管内液相 表面张力；α为土颗粒与孔隙水之间的接触角；

步骤3，

当RH或u_a-u_w为第i步增量时，单位质量土体被气体或水填充的孔隙体积用式(6)表 示：

Vpi=wiρw;---(6)]]>

式中，wⁱ为含水量；

假设固液相互作用只对土颗粒表层附近水膜密度产生显著影响，则毛细吸附区ρ_w本质 上是常数，ρ_w＝1g/cm³；

通过(7)式计算单位质量土体内水所占的孔隙体积减小量及累积孔隙体积变化量。则 每单位质量固体内水的排出量：

ΔVpi=Vpi-1-Vpi;---(7)]]>

式中，i＝2，3，4……；

步骤4，

假设接触角α＝0，通过式(5a)或式(5b)中的任何一式计算开尔文半径；

步骤5，

计算实际孔隙半径和吸附水膜厚度表达如下：

rpi=rki+ti;---(8)]]>

其中吸附水膜厚度tⁱ用Halsey方程表示如下：

ti=τ[-5ln(RHi)]1/3;---(9)]]>

式中，tⁱ为相对湿度在第i增量步时所对应的吸附水膜厚度；τ为吸附水分子的有效直 径；

假设液体水分子所占区域的横截面面积约为水的摩尔体积v_w＝1.8×10^-5m³/mol，Avogadro常数N_A＝6.02×10²³/mol，则吸附水分子的有效直径τ为：

步骤6，

用(10)式计算平均孔隙半径

(rpi)avg=rpi+rpi-12;---(10)]]>

式中，i＝2，3，4……；

步骤7，

选择孔隙体积大于0.01cm³所对应的孔径，根据下式计算黏土颗粒间孔隙等效孔 径：

d=2r=2βΣ[ΔVpi·(rpi)avg]ΣΔVpi;---(11)]]>

式中，d为黏土颗粒间孔隙等效孔径；为孔隙体积在0.01cm³以上时所对应的孔 径；β为无量纲修正系数，取值0.02；

步骤8，

将计算得到的黏土颗粒间孔隙等效孔径代入公式(12)求得最大毛细上升高度：

hmax=4Tscosαdρwg;---(12)]]>

式中，g为重力加速度。