考虑裂隙面积率的裂隙性黏土优势流双域入渗规律研究_周峙.pdf

引用格式：周峙，罗易，张家铭，等考虑裂隙面积率的裂隙性黏土优势流双域入渗规律研究安全与环境工程，（）：，y y d ，（）：考虑裂隙面积率的裂隙性黏土优势流双域入渗规律研究周峙，罗易，张家铭，孙狂飙（湖北经济学院工程管理系，湖北 武汉 ；中国地质大学（武汉）岩土钻掘与防护教育部工程研究中心，湖北 武汉 ；中国地质大学（武汉）工程学院，湖北 武汉 ；安徽省交通控股集团有限公司，安徽 合肥 ）摘要：为模拟降雨条件下裂隙性黏土中的雨水入渗过程，基于双孔隙域入渗理论和 入渗模型，引入土体裂隙平面几何参数，提出考虑裂隙面积率的裂隙性黏土优势流双域入渗模型，探究降雨强度、土体裂隙面积率和裂隙域饱和渗透系数对土体基质域与裂隙域积水时间、入渗量和入渗深度的影响，揭示土体干缩开裂后的优势流入渗规律。结果表明：提出的考虑裂隙面积率的裂隙性黏土优势流双域入渗模型形式简单、计算方便、物理意义明确且较好地体现了土体裂隙对雨水入渗过程的影响；降雨强度增大，土体基质域与裂隙域两域积水时间缩短，入渗量和入渗深度增幅较小；土体裂隙面积率增大，土体基质域积水时间不变、裂隙域积水时间延长，且土体基质域与裂隙域的入渗量分别与土体裂隙面积率呈反比和正比；土体裂隙域饱和渗透系数增大对土体基质域入渗的影响较小，仅导致土体裂隙域积水时间延长，入渗量和入渗深度增大；土体裂隙域优势流最终入渗量和入渗深度对降雨强度变化的响应不敏感，主要受土体裂隙面积率和裂隙域饱和渗透系数的控制。关键词：裂隙性黏土；土体裂隙；裂隙面积率；优势流；基质域与裂隙域；双域入渗模型中图分类号：；文章编号：（）收稿日期：开放科学（资源服务）标识码（）：基金项目：国家自然科学基金面上项目（）；岩土钻掘与防护教育部工程研究中心开放基金项目（）；湖北省水利重点科研项目（）；安徽省交控建设管理有限公司科技攻关项目（）；河南省地质矿产开发局地质科研项目（）作者简介：周峙（），男，博士，讲师，主要从事黏土优势流及致灾机理方面的研究。：通讯作者：张家铭（），男，博士，副教授，博士生导师，主要从事特殊性岩土与工程地质方面的研究。：o o o o o o ，（.，y ，；.d ，y d ，y （），；.y ，y （），；.d .，d.，）：y ，y ，y，y ，y 第 卷第期 年月安 全 与 环 境 工 程 y ，y ，y，y ，y y ，y y ，y ，y y y y o ：y；安徽沿江地区广泛存在的裂隙性黏土是导致路基边坡失稳、路基病害的主要致灾土体。土体开裂形成的雨水入渗优势通道，会导致大量的地表径流绕 过 土 体 基 质 直 接 进 入 裂 隙 底 部，形 成 优 势流。优势流使雨水快速入渗至土体深部，是造成地下水污染、水资源流失、路基边坡失稳破坏的重要影响因素。有关土体裂隙产生的优势流，目前多采用染色示踪方法对其最终的 入 渗 深 度 及 分 布 范 围 进 行 研究。已有研究结果表明，土体裂隙产生的优势流与裂隙平面几何参数高度相关。然而，由于优势流入渗的复杂性，上述土体裂隙优势流研究仍处于定性研究阶段。为了模拟雨水在具有土体裂隙等大孔隙土体中的运移过程，等、和 基于双孔隙域理论将土体孔隙分为大孔隙域和小孔隙域，认为雨水入渗只发生在大孔隙域内（如土体裂隙域），据此提出了裂隙性黏土优势流入渗模型，并利用 方程、方程、运动波方程等描述大孔隙域内的水体流动特征。然而，雨水在裂隙性黏土中的入渗同时发生在土体基质域与裂隙域内，上述假设显然与真实情况不符。为此，等 基于 方程提出了双孔隙域入渗模型，可模拟雨水在大孔隙域和小孔隙域内的入渗过程，该模型因具有明确的物理意义已得到了广泛使用。上述学者提出的裂隙性黏土优势流入渗模型包含了大量难以从试验中获取的参数，同时还需要复杂的建模及数值求解方法，不仅削弱了模型的预测能力，还限制了模型的进一步推广应用。尽管已有很多学者曾提出简化的裂隙性黏土优势流入渗模型，但仅简单地采用降雨阈值来区分基质流和优势流，物理意义不明确，且模型大多未考虑土体的初始条件（如土体含水率）及其入渗能力的时变特征，难以准确反映裂隙性黏土优势流入渗的变化规律。入渗模型因其物理意义明确、形式简单 而 被 广 泛 应 用 于 一 维 均 质 入 渗 问 题 的 研 究中 。该模型可考虑土体的初始条件（如土体含水率）及其入渗能力的时变特征，同时可预测积水时间、累积入渗量和湿润锋深度等水文参数，对研究地表径流入渗、灌溉效率和边坡稳定性评价具有重要的意义。目前仅有少量学者将 入渗模型用于模拟裂隙性黏土中的优势流，但同样存在参数过多、物理意义不明确等问题。图像处理技术在土体裂隙发展演化规律的研究中已成为一种常用的重要研究手段，诸多学者 基于该技术提取土体表面裂隙几何参数（如土体裂隙面积、宽度、长度等），研究了土体裂隙的形成机理及其影响因素，但结合土体裂隙平面几何参数定量分析雨水入渗过程的模型尚不多见。为了解决上述问题，本文基于双孔隙域入渗理论，将裂隙性黏土分为土体基质域与裂隙域，并结合 入渗模型，引入土体裂隙平面几何参数，提出考虑裂隙面积率的裂隙性黏土优势流双域入渗模型，并探讨降雨强度、土体裂隙面积率和裂隙域饱和渗透系数对土体基质域与裂隙域两域积水时间、入渗量和入渗深度的影响规律。裂隙性黏土优势流概念模型平面上将土体中多条裂隙按照等效面积与体积视为多个规则的裂隙域 图（）、（），其中为土体裂隙域所占体积权重系数。为了更好地研究土安全与环境工程 ：第 卷体裂隙优势流入渗规律，模型做了如下假设：图裂隙性黏土优势流入渗概念模型 （）降雨强度为一定值，假设雨水仅在竖直方向运移，不考虑通过裂隙壁在基质域与裂隙域之间的水分交换。（）假设土体基质域的饱和渗透系数远小于土体裂隙域的饱和渗透系数。将土体基质域和裂隙域的积水时间分别表示为，和，将降雨历时表示为。当降雨强度大于土体两域饱和渗透系数时，考虑上述假设可将雨水入渗过程分为个阶段：土体基质域积水前，。该阶段土体基质域与裂隙域的入渗能力均大于降雨强度，因此土体不产生积水，土体各域的雨水入渗速率和流量分别等于按面积占比分配的降雨强度和降雨量。土体裂隙域积水前、基质域积水后，。如图（）所示，该阶段土体基质域的入渗量将小于按其面积占比分配的降雨量，其入渗量可采用文献 提出的弱积水条件下土体的入渗速率预估模型计算得到，并假设土体基质域表面积水后，积水将直接进入裂隙域，即降雨量与土体基质域入渗量的差值为优势流入 渗 量。土 体 裂 隙 域 积 水 后，。如图（）所示，该阶段土体裂隙域达到饱和并积水，其总入渗速率将小于降雨强度，同时土体裂隙域最大入渗速率将由其饱和渗透系数控制（假设土体裂隙域内雨水入渗模式为单位水力梯度下的活塞流），土体基质域的入渗速率仍由文献 提出的弱积水条件下土体的入渗速率预估模型进行计算。基于上述假设及方法，可求得各域雨水入渗量，在已知土体基质域和裂隙域的容水孔隙空间（土体饱和体积含水率和初始体积含水率）的条件下，可进一步求得各域雨水入渗的湿润锋深度。模型控制方程推导2.1双孔隙域 Green-Ampt 入渗模型裂隙性黏土的总入渗速率可以定义如下：（）（）式中：为土体总入渗速率（），为土体水通量（），为土体表面积（），为土体入渗速率（），其下标和分别表示土体基质域与裂隙域；为土体裂隙域所占体积权重系数，本文中单位厚度为，该体积权重系数与裂隙面积率相等。将土体 基 质域 内 的 入渗视为活塞流（入渗模型中的基本假定），并对式（）应用 y定律，可得：（）（）（）（）（）式中：为 裂 隙 性 黏 土 的 整 体 饱 和 渗 透 系 数（）；为 土 体 基 质 域 的 饱 和 渗 透 系 数（）；为 土 体 裂 隙 域 的 饱 和 渗 透 系 数（）；和分别为土体基质域和土体裂隙域的压力水头（）；为竖向水力梯度。根据 入渗方程解：（）式中：为土体表面积水深度（）；为湿润锋处土体基质吸力水头（）；为湿润锋深度（），可表示为 其中，为土体累积入渗量（）；为土体容水孔隙空间（），可表示为，和分别为土体的饱和含水率和初始体积含水率。弱积水条件下，式（）可写为（）根据土体累积入渗量与入渗速率的导数关系，联立式（）与式（），可得湿润锋深度与时间的关系式为（）（）（）将、代入式（），可得：（）（）第期周峙等：考虑裂隙面积率的裂隙性黏土优势流双域入渗规律研究2.2分阶段土体基质域与裂隙域入渗方程（）当土体（基质域与裂隙域）入渗能力均大于降雨强度时，土体不产生积水。即当，时，土体各域的入渗速率等于按面积占比分配的降雨强度，有：（）（）（）（）（）则该阶段土体各域的累积入渗量为（）（）（）（）当土体裂隙域积水前、基质域积水后，即当入渗时间，时，应首先确定土体基质域的积水时间。由 入渗模型可知，当土体基质域的入渗速率等于其降雨强度时，积水便会产生。根据文献 所提出的弱积水条件下土体的入渗速率预估模型，土体基质域的入渗速率可表示为（），（）式中：为土体基质域的饱和渗透系数（）；为常数，常取值为；，为土体基质域中湿润锋处土体基质吸力水头（）；为土体基质域容水孔隙空间（）。令，联立式（）与式（），可得土体基质域积水时间，的隐式表达式为（，），（），确定后，土体基质域内的累积入渗量可表示为（），（）将式（）代入式（）并积分，可得：（），（，），（，）（）对于土体裂隙域，在其未被雨水填满之前，土体裂隙域入渗速率可由降雨强度与土体基质域的入渗速率之差求得，即：（）其相应的土体裂隙域累积入渗量为（），（，），（，）（）（）土体裂隙域积水后，当，时，首先确定土体裂隙域积水时间。由于假设土体裂隙域中的雨水入渗模式为单位梯度下的活塞流，则土体裂隙域积水饱和时的入渗速率为（）式中：为土体裂隙域的饱和渗透系数（）。当满足以下条件时，土体裂隙域产生积水：（）将式（）代入式（），可得：（）（，），（）于是，土体裂隙域积水时间，的隐式表达式为（）（，），（）该阶段土体基质域累积入渗量同式（），土体裂隙域累积入渗量可表示为，（，）（），（），（，），（）式中：，为土体基质域在裂隙域积水时的累积入渗量（）。土体裂隙域积水后，当土体缓倾斜时，整个土体表面开始产生径流，此时土体基质域与裂隙域均处于弱积水状态，两者入渗速率和入渗量仍可分别按式（）、（）和式（）、（）计算。由上述推导过程可见，通过式（）、（）、（）、（）和（），可求得降雨期间土体基质域与裂隙域任一时刻的累积入渗量。最后，土体各域的湿润锋入渗深度按照其容水孔隙空间大小可表示为（）（）（）式中：、分别为土体基质域与裂隙域的湿润锋入渗深度（）；为土体裂隙域孔隙率；为土体基质域容水孔隙空间（），土体裂隙域容水孔隙空间可设置为 。为了分析裂隙性黏土优势流入渗规律，本文定义优势流的相对入渗量为、入渗量占比为和相对入渗深度为，各指标的计算公式如下：安全与环境工程 ：第 卷（）（）（）（）（）（）总体说来，本文所提出的裂隙性黏土优势流双域入 渗 模 型 主 要 包 括 以 下 参 数：（）、（）、（）、（）、，（）、。其中，可通过土体渗透试验求得；可通过经验赋值反算获取；可通过将裂隙图像进行二值化后提取裂隙面积求得；，可由文献 提出的预估模型求得，即：，（）.（）（）式中：和为由 模型得到的土水特征曲线拟合参数。裂隙性黏土优势流入渗过程模拟及参数分析为了探讨降雨条件和裂隙发育程度对裂隙性黏土优势流入渗过程的影响，本文开展了降雨强度、土体裂隙面积率和裂隙域饱和渗透系数对裂隙性黏土优势流入渗过程的数值模拟，影响因子参数均按照倍、倍递增。是土壤非饱和渗流领域的一款成熟的计算软件，其采用有限元求解 方程，内部嵌有双重渗透模块，并有较为全面的土壤参数数据库。本文土体基质域参数选择 软件提供的黏土，具体参数见表，采用迭代方法求解隐函数根。土体裂隙域参数及不同工况模拟参数见表，其中工况至工况表示不同降雨强度对土表土体基本参数表 y y 土体类型 残余体积含水率，饱和体积含水率，初始体积含水率，（）基质域饱和渗透系数（）黏土 表不同工况模拟参数表 工况模型参数（）降雨历时 ，（）（）工况 工况 工况 工况 工况 工况 工况 工况 工况 体基质域与裂隙域累积入渗量的影响，工况至工况表示不同土体裂隙面积率对土体基质域与裂隙域累积入渗量的影响，工况至工况表示不同土体裂隙域饱和渗透系数对土体基质域与裂隙域累积入渗量的影响。不同降雨强度对土体基质域与裂隙域累积入渗量的影响分析图为不同降雨强度条件下土体基质域与裂隙域入渗过程的模拟结果。图（）（）为不同降雨强度条件下土体基质 域 与 裂 隙 域 的 累 积 入 渗 量 变 化 曲 线。由图（）（）可知：不同降雨强度下土体基质域累积入渗量均大于裂隙域，且随降雨强度的增大土体两域累积入渗量均有所增加，但后期增幅较小；土体基质域与裂隙域积水时间均随降雨强度的增大而显著减小，且两者时间间隔越来越小。图（）（）为不同降雨强度下土体基质域与裂隙域的入渗深度变化曲线。由图（）（）可知：不同降雨强度下土体基质域入渗深度曲线先陡增后缓增，土体裂隙域入渗深度前期增长缓慢小于基质域，随后以近似线性趋势快速增长并超过土体基质域入渗深度；随降雨强度的增大，土体裂隙域入渗深度缓增段时间缩短，其超过基质域的时间也逐渐减小。图（）（）为不同降雨强度下优势流相对入渗量、入渗量占比和入渗量相对入渗深度的变化曲线。由图（）（）可知：总体而言，不同降雨强度下各优势流入渗指标的变化趋势较为一第期周峙等：考虑裂隙面积率的裂隙性黏土优势流双域入渗规律研究致，均在土体基质域积水前保持为一常数，随后随降雨时间先陡增后趋缓；优势流入渗指标在降雨初期随降雨强度的增大而增加，但随降雨时间的增加各指标间差距逐渐减小，至降雨末期，最大降雨强度（）下的优势流入渗指标小于降雨强度为 时的优势流入渗指标。tp c,tp m,tp m,tp m,=0.87 mintp c,=0.9 mintp c,1.20.80.40累积入渗量/cm1.51.00.50累积入渗量/cm1.51.00.50累积入渗量/cm时间/min时间/min020406080100 120020406080100 120020406080100 120时间/min时间/min时间/min020406080100120020406080100120020406080100120时间/min121086420入渗深度/cm121086420入渗深度/cm121086420入渗深度/cm0 400.320.240.160.08.优势,相对入渗量=/f IImc优势,入渗量/比=/F IIctotal0.320.240.160.083.63 02 41 81 20 6.优势,相对入渗深度=/D zzcm020406080100 120时间/min020406080100 120时间/min020406080100 120时间/mintp c,tp m,rrr=0.018 cm/min0.036 cm/min0.072 cm/mintp c,tp m,rrr=0.018 cm/min0.036 cm/min0.072 cm/mintp c,tp m,rrr=0.018 cm/min0.036 cm/min0.072 cm/min=,=0.050.096 cm/minKc=,=0.050.096 cm/minKc=,=0.050.096 cm/minKc(a)0.018 cm min降雨强度=/时累积入渗量变化曲线r(b)0.036 cm min降雨强度=/时累积入渗量变化曲线r(c)0.072 cm min降雨强度=/累积入渗量变化曲线r(d)0.018 cm min降雨强度=/时r入渗深度变化曲线(e)0.036 cm min降雨强度=/r时入渗深度变化曲线(f)0.072 cm min降雨强度=/r时入渗深度变化曲线深 度不同降雨强度下优势,入渗量/比变化曲线h不同降雨强度下优势,相对入渗量变化曲线(g)深 度不同降雨强度下优势,相对入渗深度变化曲线i0.050.096=,=Kccm/min基质=裂隙=时0.050.096=,=Kccm/min基质=裂隙=0.050.096=,=Kccm/min基质=裂隙=0.050.096=,=Kccm/min基质=裂隙=0.050.096=,=Kccm/min基质=裂隙=0.050.096=,=Kccm/min基质=裂隙=t=30.6 mint=10.9 mint=8.9 min图不同降雨强度条件下土体基质域与裂隙域入渗过程的模拟结果 不同裂隙面积率对土体基质域与裂隙域累积入渗量的影响分析图为不同裂隙面积率条件下土体基质域与裂隙域入渗过程的模拟结果。图（）（）为不同裂隙面积率条件下土体基质域与裂隙域的累积入渗量变化曲线。由图（）（）可知：当降雨强度和土体裂隙域饱和渗透系数保持不变时，随裂隙域所占体积权重系数的增大，土体基质域与裂隙域的累积入渗量曲线逐步靠近直至相交。其中，土体基质域最终累积入渗量减小，土体隙域累积入渗量增长速度及增幅显著增大；当.时，土体裂隙域累积入渗量曲线超过基质域。此外，随增大，土体基质域积水时间不变，裂隙域积水时间显著增大。图（）（）为不同裂隙面积率条件下土体基质域与裂隙域的入渗深度变化曲线。由图（）（）可知：随增大，土体基质域入渗深度曲线保持不变，土体裂隙域入渗深度曲线缓增段时间延长，其超过基质域的时间也逐渐增长，最终其入渗深度减小。图（）（）为不同土体裂隙面积率条件下优势流相对入渗量、入渗量占比和相对入渗深度的变化曲线。由图（）（）可知：优势流相对入渗量和入渗量占比均随裂隙面积率的增加而增加，当 .时，优势流入渗量占比达 ；相比之下，优势流入渗量相对入渗深度与相对入渗量的变化趋势相反，其随裂隙面积率的增加而减小。安全与环境工程 ：第 卷rK=0.018 cm/mincm/minc,=0.096基质域裂隙域(c)0.2裂隙面积率=时累积入,量变化曲线tp c,tp m,tp c,tp m,tp c,tp m,1.51.00.50累积入,量化 cm1.51.20.90.60.302.42.11.81.51.20.90.60.30累积入,量化 cm累积入,量化 cm020406080100120时间化 min020406080100120时间化 min020406080100120时间化 min020406080100120时间化 min020406080100120时间化 min020406080100120时间化 min020406080100120时间化 min020406080100120时间化 min020406080100120时间化 min121086420入,深度化 cm121086420入,深度化 cm121086420入,深度化 cm2.01.61.20.80.40优势流相对入,量=化f IImc优势流相对入,深度=化D zzcm优势流入,量占比=化F IIctotal0.90.60.303.22.82.42.01.61.20.80.4(a)0.05裂隙面积率=时累积入,量变化曲线(b)0.1裂隙面积率=时累积入,量变化曲线(d)0.05裂隙面积率=时入,深度变化曲线(e)0.1裂隙面积率=时入,深度变化曲线(f)0.2裂隙面积率=时入,深度变化曲线(g)不=裂隙面积率 下优势流相对入,量变化曲线(h)不=裂隙面积率 下优势流入,量占比变化曲线(I)不=裂隙面积率 下优势流相对入,深度变化曲线rK=,=0.036 cm/min0.096 cm/minc=0.050.10.2tp c,tp m,=0.050.10.2=0.050.150.25rK=,=0.036 cm/min0.096 cm/mincrK=,=0.036 cm/min0.096 cm/minctp c,tp m,tp c,tp m,rK=0.036 cm/mincm/minc,=0.096基质域裂隙域rK=0.036 cm/mincm/minc,=0.096基质域裂隙域rK=0.018 cm/mincm/minc,=0.096基质域裂隙域rK=0.036 cm/mincm/minc,=0.096基质域裂隙域rK=0.036 cm/mincm/minc,=0.096基质域裂隙域t=10.9 mint=11.7 mint=14.6 min图不同裂隙面积率条件下土体基质域与裂隙域入渗过程的模拟结果 不同裂隙域饱和渗透系数对土体基质域与裂隙域累积入渗量的影响分析图为不同裂隙域饱和渗透系数条件下土体基质域与裂隙域入渗过程的模拟结果。图（）（）为不同裂隙域饱和渗透系数条件下土体基质域与裂隙域的累积入渗量变化曲线。由图（）（）可知：当土体裂隙域所占体积权重系数和降雨强度保持不变时，随裂隙域饱和渗透系数的增大，土体基质域与裂隙域的累积入渗量曲线逐步靠近。其中，土体基质域累积入渗量不变，土体裂隙域累积入渗量增长速度及增幅显著增大，可以推测，当足够大时，土体裂隙域累积入渗量曲线终将超过基质域。同时，随着增大，土体基质域积水时间不变，裂隙域积水时间显著增大。图（）（）为不同裂隙域饱和渗透系数条件下土体基质域与裂隙域的入渗深度变化曲线。由图（）（）可知：随增大，土体基质域入渗深度曲线保持不变，土体裂隙域入渗深度增长速率显著增加，超过基质域入渗深度的时间也显著缩短，其最终入渗深度增幅分别达到 和 。图（）（）为不同裂隙域饱和渗透系数下优势流相对入渗量、入渗量占比和相对入渗深度的变化曲线。由图（）（）可知：不同裂隙域饱和渗透系数条件下各优势流入渗指标均在土体裂隙域积水前一段时间基本保持不变，且各入渗指标值相同（如均为 ，均为 ，均为 ），随后随降雨时间先陡增后缓增；优势流入渗指标在降雨初期随增大而减小，但随降雨时间的增加，值较大的各优势流入渗指标快速增加，各指标间的差距逐渐增大。第期周峙等：考虑裂隙面积率的裂隙性黏土优势流双域入渗规律研究(g)不同裂隙域饱和渗透系数下优势流相对入渗量变化曲线(h)不同裂隙域饱和渗透系数下优势流入渗量占比变化曲线(I)不同裂隙域饱和渗透系数下优势流相对入渗深度变化曲线(c)0.192 cm min裂隙域饱和渗透系数=流时=积入渗量变化曲线Kc(a)048 c裂隙域饱和渗透系数时=积入渗量变化曲线Kc=流0.m min(b)0.096 cm min裂隙域饱和渗透系数=流时=积入渗量变化曲线Kc(d)0.048 cm min裂隙域饱和渗透系数=流时入渗深度变化曲线Kc(e)0.096 cm min裂隙域饱和渗透系数=流时入渗深度变化曲线Kc(f)0.192 cm min裂隙域饱和渗透系数=流时入渗深度变化曲线Kc020406080100 120时间流 min020406080100 120时间流 min020406080100 120时间流 min020406080100120时间流 min020406080100120时间流 min020406080100120时间流 min020406080100 120时间流 min020406080100120时间流 min020406080100120时间流 mintp m下tp c下tp c下tp m下tp c下=3.79tp m下=3.651.51.20.90.60.301.51.00.50=积入渗量流 cm=积入渗量流 cm=积入渗量流 cm1.51.20.90.60.30121086420入渗深度流 cm24201612840入渗深度流 cm6543210入渗深度流 cm8765432100.60.50.40.30.20.11.00.80.60.40.20优势流相对入渗量=流f IImc优势流入渗量占比=流F IIctotal优势流相对入渗深度=流D zzcmtp c下tp m下tp c下tp m下tp c下tp m下=下=0.050.036 cm/minrKKKccc=0.048 cm/min0.096 cm/min0.192 cm/minKKKccc=0.048 cm/min0.096 cm/min0.192 cm/minKKKccc=0.048 cm/min0.096 cm/min0.192 cm/min=下=0.050.036 cm/minr=下=0.050.036 cm/minr0.05=下 r=0.036 cm/min基质域裂隙域0.05=下 r=0.036 cm/min基质域裂隙域0.05=下 r=0.036 cm/min基质域裂隙域0.05=下 r=0.036 cm/min基质域裂隙域0.05=下 r=0.036 cm/min基质域裂隙域0.05=下 r=0.036 cm/min基质域裂隙域t=40.2 mint=10.9 mint=7.13 min图不同裂隙域饱和渗透系数条件下土体基质域与裂隙域入渗过程的模拟结果 y 讨论4.1降雨强度对裂隙性黏土雨水入渗过程的影响降雨强度增大不仅使土体基质域与裂隙域两域的累积入渗量增加，而且还显著缩短了土体两域的积水时间。但两域积水后土体入渗速率将不再由降雨强度控制，不同降雨强度下土体两域的最终累积入渗量主要取决于积水后的降雨时长。当降雨强度增大时，土体两域积水时间与间隔越来越短，积水后的降雨时长差值也相应减小。这是由于土体裂隙域面积占比较小，雨水下渗所要充填的孔隙空间体积较小，尽管土体裂隙域的累积入渗量小于基质域，但其湿润锋下渗速度快于基质域。当土体裂隙域湿润锋入渗深度超过基质域时，可视为优势流产生的标志时间。由计算结果可知，随降雨强度的增大，优势流现象产生速度加快，这与前人试验观察得到暴雨期间的优势流产生频率更高的结果一致，。裂隙面积率和裂隙域饱和渗透系数对裂隙性黏土雨水入渗过程的影响土体裂隙面积率通过影响各域的受水面积及容水孔隙空间进而影响入渗量和入渗深度。当土体裂隙面积率增大时，土体基质域受雨面积减小，其累积入渗量与容水空间体积同比例减小，因此土体基质域内最终入渗深度不变。土体裂隙面积率降低使优势流入渗量减小，但同时土体裂隙域内相同深度的容水空间减小幅度更大，因此其优势流相对入渗深度反而增大。当土体裂隙面积率分别为原先的、（从.分别减小至.、.）时，土体裂隙域容水空间与之等比例缩小，优势流最终入渗量减小安全与环境工程 ：第 卷 和 。该模拟结果或可解释诸多前人试验中的普遍现象 ，即在裂隙或者大孔隙不易被观察到的土体中，优势流仍使雨水快速下渗至土体深部。土体裂隙域饱和渗透系数变化仅影响土体裂隙域的入渗过程，这是由于土体裂隙域饱和渗透系数的增大加快了优势流现象的产生，尤其使其入渗深度明显增大。实际过程中，较大的土体裂隙域渗透系数往往对应着土体裂隙的发育程度较高、连通性较好，能使雨水更快入渗至土体深部。不同因素对裂隙性黏土参数指标的敏感性分析为了评价上述不同因素对裂隙性黏土雨水入渗过程的影响程度，对各影响因素增加至倍、倍时的土体基质域和裂隙域的两域入渗量、入渗深度、优势流入渗量占比和优势流产生时间的变化规律进行了分析，其分析结果见表。表裂隙性黏土双域入渗指标随影响因素的变化情况 影响因素因素增加倍数倍双域入渗指标变化倍数倍土体基质域入渗量土体基质域入渗深度土体裂隙域入渗量土体裂隙域入渗深度优势流入渗量占比优势流产生时间降雨强度 土体裂隙面积率 土体裂隙域饱和渗透系数 由表可知：对于土体基质域入渗过程，主要受降雨强度和裂隙面积率的影响，但影响程度较小；对于土体裂隙域入渗过程，土体裂隙面积率和土体裂隙域饱和渗透系数的增加均会显著增加其入渗量和入渗量占比，但优势流入渗深度更易受土体裂隙域饱和渗透系数的影响；裂隙性黏土优势流现象产生时间受降雨强度和土体裂隙域饱和渗透系数的影响较大。4.4计算模型讨论本文所提出的裂隙性黏土优势流双域入渗模型采用 入渗方程描述土体基质域与裂隙域的雨水入渗过程，并假设土体裂隙域中的水力梯度为单位梯度。土体中的毛细水作用力（如湿润锋处的基质势）取决于土体的进气值与孔隙分布特征，由于土体裂隙的孔隙较大，进气值极低，因此单位水力梯度的假设是合理的。然而，当裂隙开度较大时，土体基质域积水将会以薄膜流（幕流）的形式沿裂隙壁快速渗入裂隙底部 ，本文基于单位水力梯度下的土体裂隙域活塞流入渗模型可能会低估优势流的入渗深度。但是，该过程或可通过调整土体裂隙域的容水孔隙率项（）的表达形式来进行校准。此外，雨水入渗时，土体基质域遇水膨胀使土体裂隙逐渐闭合，有研究 认为土体裂隙闭合会极大地降低优势流的入渗速率和入渗深度，甚至导致优势流现象的消失。但也有学者 研究发现，即使土体裂隙闭合后，优势流现象将继续存在且其入渗速率和入渗深度仍显著大于土体基质域。本文所提出的裂隙性黏土优势流入渗模型暂未考虑裂隙闭合过程，可能高估了优势流的入渗量，后续研究可通过改变土体裂隙面积率的表达形式来描述这一动态过程。尽管本文所提出的裂隙性黏土优势流入渗模型存在上述问题，但其形式简单、计算方便、物理意义明确且较好地体现了土体裂隙对雨水入渗过程的影响，为研究降雨条件下裂隙性黏土优势流入渗规律提供了新思路。结论与建议（）对于土体基质域入渗，降雨强度和裂隙面积率对其入渗的影响较小，两域积水时间会随降雨强度的增大而缩短，其入渗量和入渗深度增幅较小。（）对于土体裂隙域入渗，降雨强度增大会加快优势流现象的产生速度，土体裂隙面积率和裂隙域饱和渗透系数增加均会显著增加其入渗量和入渗量占比，尤其是土体裂隙域饱和渗透系数的增加会大幅增大优势流入渗深度和延长积水时间。（）土体裂隙优势流最终入渗量和入渗深度对降雨强度变化的响应不显著，主要受土体裂隙面积率和裂隙域饱和渗透系数的控制。（）本文所提出的裂隙性黏土优势流入渗模型适用于土体裂隙开度较小的情况，未来针对该模型的改进应进一步考虑裂隙侧壁幕流的影响、土体裂第期周峙等：考虑裂隙面积率的裂隙性黏土优势流双域入渗规律研究隙动态变化过程以及两域之间的水平入渗过程。参考文献：周峙，张家铭，宁伏龙，等降雨入渗下裂土边坡水分运移时空特征与失稳机理 交通运输工程学报，（）：何岱洵，张家明，陈茂，等 极端干湿循环作用下饱和泥浆红黏土的干缩裂隙发育特征安全与环境工程，（）：，y ，y d y，（）：y 张文杰，杨金康 填埋场垃圾中优先流特征的染色示踪试验岩土力学，（）：，y ，（）：（），（）：梁冰，张柴，刘磊基于优势流效应的填埋场渗沥液回灌过程预测科学技术与工程，（）：吕捷，樊秀峰，吴振祥 基于降雨染色示踪试验的大孔隙流特性研究长江科学院院报，（）：，：y y 赵宽耀，许强，刘方洲，等黄土中优势通道渗流特征研究岩土工程学报，（）：，y d y，（）：（）d ，（）：，y ，：，y ，（）：，（）：许旭堂，简文彬，吴能森，等 降雨诱发残积土坡失稳的模型试验中 国 公 路 学 报，（）：常金源，包含，伍法权，等 降雨条件下浅层滑坡稳定性探讨岩土力学，（）：孙狂飙，袁超，周峙，等人工模拟降雨条件下坡面侵蚀特性的模型试验研究安全与环境工程，（）：y ：，y y y d y，（）：y ，y y ，（）：（）张家铭，罗易，周峙，等 基于足尺模型试验的边坡裂隙发展演化规律中南大学学报（自然科学版），（）：倪可，蔡奕，魏乾坤，等 柴油污染粉质黏土干裂隙的特征 中南大学学报（自 然 科 学 版），（）：罗易，张家铭，周峙，等降雨蒸发条件下土体开裂临界含水率演变规律研究岩土力学，（）：，（）：y ，y ，y ：y ，（）：梁越，陈鹏飞，林加定，等基于透明土技术的多孔介质孔隙流动特性研究岩土工程学报，（）：y d ，（）：，y，d ，（）：（）程世涛，张维，任鹏 非饱和裂隙膨胀土边坡降雨入渗和变形的数值模拟 四川建筑，（）：魏凌傲，郜迪，王鹏，等降雨条件下膨胀土裂隙边坡深层浅层滑坡渗透稳定性分析水利科学与寒区工程，（）：，y d ，：安全与环境工程 ：第 卷