一种预测饱和黏性土非线性渗透的分形方法

技术领域

本发明涉及岩土工程领域，特别是一种预测饱和黏性土非线性渗透的分形方法。

背景技术

饱和黏性土的非线性渗透现象是一种低渗透介质渗流现象，它的推导预测与海洋工程、石油勘探与开发、岩土工程等领域紧密联系，由于低渗透介质非均质性、渗透性较小等特性导致其渗流过程的机理相较复杂难以进行预测，对隧道工程的设计和堤坝防渗等重大项目的开展都带来了巨大的考验。因此，如何准确且简便地预测低渗透介质非线性渗透是对隧道工程的设计和堤坝防渗具有科学的指导作用。

然而，目前非线性渗透多以经验公式为主，而已有理论预测模型，如陶高梁等人所提出的非线性渗透模型(Tao G,Huang Z,Xiao H,et al.A new nonlinear seepage modelfor clay soil considering the initial hydraulic gradient of microscopicseepage channels[J].Computers and Geotechnics,2023,154:105179.)，对饱和黏性土非线性渗透的预测计算较为复杂，部分理论预测模型无法进行不同水力梯度下非线性渗透的连续性预测，预测效果也有待提高。故需要提出一种更加有效的预测手段用于解决现有预测方法所存在的不足。

发明内容

本发明的目的在于，提供一种预测饱和黏性土非线性渗透的分形方法，用于解决现有非线性渗透模型无法进行不同水力梯度下非线性渗透的连续性预测，且预测精度较低的问题。

为解决上述技术问题，本发明提供了一种预测饱和黏性土非线性渗透的分形方法，包括如下步骤：

S1.基于微观渗透通道的起始水力梯度，构建饱和黏性土非线性渗透分形模型；

S2.利用所述饱和黏性土非线性渗透分形模型，计算得到不同水力梯度下饱和黏性土的渗流速度预测值，通过所述渗流速度预测值对不同水力梯度下饱和黏性土的渗流状态进行预测；

所述饱和黏性土非线性渗透分形模型的表达式为：

其中，所述的v

本发明的有益效果是：区别于现有技术的情况，本发明提供了一种预测饱和黏性土非线性渗透的分形方法，以陶高梁提出的非线性渗透模型为基础结合分形理论，提出了考虑微观渗透通道起始水力梯度的新的饱和黏性土非线性渗透分形模型，解决了陶高梁的非线性渗透模型的局限性；该模型能对非线性渗透进行连续性预测，简化现有非线性渗透模型计算步骤，弥补了现有非线性渗透模型只能进行散点值预测而无法进行连续预测的不足，提高了预测精度，为饱和黏性土渗透的理论体系的系统建立提供了一种新的理论方法。

附图说明

图1是本发明中基于土壤水分特征曲线SWRC的孔隙渗流通道的分布特征模型图；

图2是本发明中土体孔隙渗透通道分形剖面图；

图3a是本发明实施例1中湖南红黏土的SWCC曲线的分形维数求解图；

图3b是本发明实施例2中湖南红黏土的SWCC曲线的分形维数求解图；

图3c是本发明实施例3中湖南红黏土的SWCC曲线的分形维数求解图；

图4a是本发明实施例1中饱和湖南黏性土渗流速度实测值与预测值对比图；

图4b是本发明实施例2中饱和湖南黏性土渗流速度实测值与预测值对比图；

图4c是本发明实施例3中饱和湖南黏性土渗流速度实测值与预测值对比图。

具体实施方式

下面将结合本发明实施例，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，均属于本发明保护的范围。

对于本发明中所提出的预测饱和黏性土非线性渗透的分形方法，具体步骤如下：

S1.基于微观渗透通道的起始水力梯度，构建饱和黏性土非线性渗透分形模型。S1的具体步骤如下：

S11.构建基于土壤水分特征曲线SWRC的孔隙渗流通道的分布特征模型。该步骤中，SWRC表示基质吸力-饱和度表征的土水保持曲线，图2中该模型的土壤颗粒中具有若干孔径不等的孔隙通道，渗流的孔隙通道记为渗流孔，而未渗流的孔隙通道记为无渗流孔。具体地，当基质吸力达到土体进气值ψ

S12.基于孔隙渗流通道的分布特征模型，得到土体的饱和渗透系数k

(1)式中：γ为层流阻力系数，u为黏度，r

将不同级别孔隙通道的渗透系数进行叠加得到土体的饱和土渗透系数，具体表达式如下：

(2)式中：A

基于S1步骤中基于孔隙渗流通道的分布特征模型，相应第i级孔隙通道的横截面面积为：

(3)式中：L为实际土体的长度，p

简化之后相应第i级孔隙通道的渗透系数表示为：

(4)式中：k

基于(4)式叠加求和，土体饱和渗透系数最终表示为：

(5)式中：k

S13.构建考虑微观渗透通道起始水力梯度的黏性土非线性渗透简化模型。该步骤中，将饱和土体中不同尺寸孔隙通道的毛细水开始发生迁移时，饱和孔隙通道的上下压力差ΔP

ΔP

(6)式中：ψ

结合微观孔隙通道的水压差ΔP与基质吸力ψ之间的关系式(6)和伯努利方程可得到不同尺寸第i级孔隙通道的起始水力梯度I

(7)式中：I

土体的渗流速度可以视为所有发生渗流的土体孔隙渗透通道的渗流速度叠加，如式(8)所示：

v

(8)式中：v

通过式(4)所计算的第i级孔隙通道渗透系数方程与饱和黏性土非线性渗透模型中式(7)所计算的第i级孔隙通道对应的起始水力梯度方程代入式(8)中，得到不同水力梯度下的饱和黏性土整体的渗流速度预测模型，如式(9)所示：

对式(10)进行合并可得渗流速度预测模型简化形式：

S14.基于土体孔隙的分形特征模型，得到不同水力梯度下饱和黏性土非线性渗透分形模型。

由于土壤水分特征曲线SWRC具有良好的分形特性，本质上反映了土体孔隙分布具有较强的分形行为；从而SWRC可以看作反映土体不同尺寸孔隙渗透通道的间接指标，土体渗透通道由SWRC所反映的许多大小不同的孔隙通道组成，这种孔隙渗透通道分布形式同样具有良好的分形特性。图2所示土体孔隙渗透通道分形剖面图，土体孔隙分形渗透通道中，当基质吸力达到ψ

该步骤中，假设最大级别渗透通道孔径r

f(r)＝Cr

(11)式中：常数C＝(3-D)/(k

根据函数(12)可得到土体中尺寸小于或者等于r的孔隙的总体积V(≤r)，表达式如下：

假设当土体孔隙中尺寸小于或者等于r的孔隙通道完全充满水，乘以水的密度，则单位土体质量含水率w表达式为：

当土体中最小尺寸r

根据Young-Laplace方程可得基质吸力ψ与有效孔径r之间的关系，表示为：

(15)式中：T

不同尺寸的孔隙通道所对应的基质吸力大小不同，其中最大尺寸孔径的孔隙通道所对应的基质吸力可以视为进气值ψ

将式(15)、式(16)分别代入式(13)和式(14)中，可分别得到如下关系式：

通过分形理论和毛细理论结合式(17)和式(18)提出了新的土壤水分特征曲线的分形模型：

土体整体处于饱和条件时，质量含水率w

(20)式中：G

将(20)式带入(19)式中可得土壤水分特征曲线质量含水率w表征的分形模型具体表达式：

通过质量含水率w与体积含水率θ和饱和度S

(21)～(23)式中：饱和度S

将式(10)的渗流速度预测模型进行简化得到积分形式：

(24)式中：θ

对式(22)进行求导，可得分形模型体积含水率θ与基质吸力ψ之间的微分关系式：

dθ＝θ

将(25)式代入(24)式中，可得基质吸力范围表征的非线性渗透分形模型，表达式如下：

(26)式中根据积分定义可得，等式建立之后原非线性渗流模型中ψ

对式(26)求解积分之后，可得饱和黏性土非线性渗透分形模型的表达式：

S2.利用所述饱和黏性土非线性渗透分形模型，计算得到不同水力梯度下饱和黏性土的渗流速度预测值。S2的具体步骤如下：

S21.标定不同尺寸的孔隙通道中发生饱和渗流时的起始水压差与利用轴平移技术测得的基质吸力之间的关联系数c。该步骤中，测定土体进气值ψ

S22.计算不同加压下的实际水力梯度I。该步骤中，根据SWRC所反映分段形式，通过伯努利方程计算渗流时施加的不同大小水压差ΔP

S23.计算分维数D。该步骤中，选取基质吸力ψ>ψ

S24.标定k

S25.基于所述饱和黏性土非线性渗透分形模型，计算得到不同水力梯度下饱和黏性土的渗流速度预测值V

下面通过具体实施例对上述预测饱和黏性土非线性渗透的分形方法的预测效果进行表征分析。实施例1～3的土样均选自湖南红黏土，具体土样干密度情况如表1所示。

表1

实施例1～3的土样均采用前述S2步骤所述方法进行预测，利用已测的非线性渗透试验结果，采用前述步骤S23进行分维数的求解，结果如附图3a～3c所示；采用前述步骤S25计算得到不同水力梯度下饱和黏性土的渗流速度预测值并进行拟合，同时与现有技术中陶高梁采用的非线性渗流模型以及Halex模型所得预测结果进行对比，所得模型预测精度的验证如附图4a～4c所示。

由附图3a～3c结果可知，除干密度为1.5g/cm

由附图4a～4c示三种渗透模型的预测值与实测值的对比可知，在三种干密度饱和湖南黏性土低水力梯度下，本文非线性渗透分形模型更接近于实测值。因此本文提出的渗透分形新模型的预测值与实测值数据更为吻合，预测效果更加准确。

与Halex模型相比，本发明研究提出的模型对黏土的非线性渗流趋势提供了全面的描述，预测值与实测值之间显示出良好的一致性。相反，Halex模型无法捕捉黏土的非线性渗流趋势，预测性能较差。此外，本研究提出的模型和陶高梁的非线性渗流模型都能够描述黏土的非线性渗流行为。然而，陶高梁的模型无法预测非线性连续渗流，并且在干密度增加时预测性能降低。相反，本研究中提出的饱和粘土非线性渗流分形模型能够成功预测非线性连续渗流。它简化了陶高梁模型的复杂计算过程，同时提高了预测精度。重要的是，该模型在预测干密度较高的饱和粘土的渗流速度方面仍然具有适用性，为后续的工程运用更好的提供理论指导。

区别于现有技术的情况，本发明提供了一种预测饱和黏性土非线性渗透的分形方法，以陶高梁提出的非线性渗透模型为基础结合分形理论，提出了考虑微观渗透通道起始水力梯度的新的饱和黏性土非线性渗透分形模型，解决了陶高梁的非线性渗透模型的局限性；该模型能对非线性渗透进行连续性预测，简化现有非线性渗透模型计算步骤，弥补了现有非线性渗透模型只能进行散点值预测而无法进行连续预测的不足，提高了预测精度，为饱和黏性土渗透的理论体系的系统建立提供了一种新的理论方法。

需要说明的是，以上各实施例均属于同一发明构思，各实施例的描述各有侧重，在个别实施例中描述未详尽之处，可参考其他实施例中的描述。

以上所述实施例仅表达了本发明的实施方式，其描述较为具体和详细，但并不能因此而理解为对发明专利范围的限制。应当指出的是，对于本领域的普通技术人员来说，在不脱离本发明构思的前提下，还可以做出若干变形和改进，这些都属于本发明的保护范围。因此，本发明专利的保护范围应以所附权利要求为准。