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

摘要

本发明公开了一种黏土颗粒间孔隙等效孔径的方法，拟合得到完整的土-水特征曲线，选取基质吸力于144000kPa～0kPa之间的数据，将基质吸力转化为相对湿度RH，将质量含水量转化为每单位质量土体内水填充孔隙的体积，计算开尔文半径，然后计算实际孔隙半径和水膜的厚度，再计算平均孔隙半径，计算土体的等效毛细直径，最后计算最大毛细上升高度。以完整的土-水特征曲线为载体，建立在已有研究的基础之上，通过大量实验并结合统计理论，推导得到不同黏粒含量土样孔径的等效方法，为最大毛细上升高度的计算提供了一种有利的途径。

主权项

一种黏土颗粒间孔隙等效孔径的方法，其特征在于，具体按照以下步骤进行：步骤1，首先制备含黏粒土样进行土‑水特征实验，并通过MATLAB用(1)式拟合得到完整的土‑水特征曲线，以质量含水量表示，拟合公式如下：$\mathrm{\θ}=({\mathrm{\θ}}_s-{\mathrm{\θ}}_r){\[\frac{1}{1+a^n{\left(u_a-u_w\right)}^n}\]}^m+{\mathrm{\θ}}_r;---\left(1\right)$式中，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 \left(-\frac{(u_a-u_w)v_w}{RT}\right);---\left(4\right)$由以上推导可得，毛细半径r以相对湿度表示形式如下：$r=-\frac{2T_s{v}_{w}cos\mathrm{\α}}{RTln\left(RH\right)};---\left(5a\right)$或表示为基质吸力形式：$r=\frac{2T_s\cos \mathrm{\α}}{u_a-u_w};---\left(5b\right)$式中，v_w为水蒸气的偏摩尔体积；R为通用气体常数；T为热力学温度；T_s为毛细管内液相表面张力；α为土颗粒与孔隙水之间的接触角；步骤3，当RH或u_a‑u_w为第i步增量时，单位质量土体被气体或水填充的孔隙体积用式(6)表示：$V_p^i=\frac{w^i}{{\mathrm{\ρ}}_w};---\left(6\right)$式中，wⁱ为含水量；假设固液相互作用只对土颗粒表层附近水膜密度产生显著影响，则毛细吸附区ρ_w本质上是常数，ρ_w＝1g/cm³；通过(7)式计算单位质量土体内水所占的孔隙体积减小量及累积孔隙体积变化量。则每单位质量固体内水的排出量：${\mathrm{\ΔV}}_p^i=V_p^{i-1}-V_p^i;---\left(7\right)$式中，i＝2，3，4……；步骤4，假设接触角α＝0，通过式(5a)或式(5b)中的任何一式计算开尔文半径；步骤5，计算实际孔隙半径和吸附水膜厚度表达如下：$r_p^i=r_k^i+t^i;---\left(8\right)$其中吸附水膜厚度tⁱ用Halsey方程表示如下：$t^i=\mathrm{\τ}{\[-\frac{5}{ln\left({\mathrm{RH}}^i\right)}\]}^{1/3};---\left(9\right)$式中，tⁱ为相对湿度在第i增量步时所对应的吸附水膜厚度；τ为吸附水分子的有效直径；假设液体水分子所占区域的横截面面积约为水的摩尔体积v_w＝1.8×10^‑5m³/mol，Avogadro常数N_A＝6.02×10²³/mol，则吸附水分子的有效直径τ为：步骤6，用(10)式计算平均孔隙半径${\left(r_p^i\right)}_{avg}=\frac{r_p^i+r_p^{i-1}}{2};---\left(10\right)$式中，i＝2，3，4……；步骤7，选择孔隙体积大于0.01cm³所对应的孔径，根据下式计算黏土颗粒间孔隙等效孔径：$d=2r=\frac{2}{\mathrm{\β}}\frac{\mathrm{\Σ}\[{\mathrm{\ΔV}}_p^i\·{\left(r_p^i\right)}_{avg}\]}{{\mathrm{\Σ\ΔV}}_p^i};---\left(11\right)$式中，d为黏土颗粒间孔隙等效孔径；为孔隙体积在0.01cm³以上时所对应的孔径；β为无量纲修正系数，取值0.02；步骤8，将计算得到的黏土颗粒间孔隙等效孔径代入公式(12)求得最大毛细上升高度：$h_{max}=\frac{4T_{s}cos\mathrm{\α}}{{\mathrm{d\ρ}}_{w}g};---\left(12\right)$式中，g为重力加速度。