土壤水分特征曲线中的孔隙连通性对非饱和渗透系数预测偏差的影响研究综述

土壤水分特征曲线中的孔隙连通性对非饱和渗透系数预测偏差的影响研究综述土壤非饱和渗透系数是描述水分在非饱和土壤中运移能力的核心参数，广泛应用于农田灌溉管理、地下水补给评估、污染物迁移模拟及边坡稳定性分析等领域。当前非饱和渗透系数的预测模型大多以土壤水分特征曲线（SoilWaterRetentionCurve,SWRC）为基础，通过假定孔隙尺寸分布与水力传导特性的对应关系实现间接估算。但大量野外实测与室内试验结果表明，这类模型的预测结果与实测值常存在1~3个数量级的偏差，其核心原因之一是传统SWRC仅能表征孔隙尺寸与含水率、基质势的静态关系，无法反映孔隙之间的连通性特征，而连通性恰恰是控制水分在非饱和状态下流动路径连续性的关键因素。一、土壤孔隙连通性的内涵与表征方法土壤孔隙连通性是指孔隙空间在三维空间中相互连接的程度，可分为几何连通性与水力连通性两个层面：几何连通性描述孔隙之间的物理连接状态，仅关注孔隙结构本身的拓扑特征；水力连通性则进一步考虑水分占据孔隙的情况，只有被连续水相填充的孔隙才会对水流产生贡献。在非饱和状态下，随着基质势降低（土壤逐渐脱水），大孔隙首先排空，水分仅赋存于较小的孔隙中，此时即使几何上连通的孔隙，若被空气阻断也会失去水力传导功能，因此水力连通性是影响非饱和渗透系数的直接因素。目前孔隙连通性的表征方法主要分为三类：拓扑学参数法：通过孔隙网络模型提取连通性指标，常见参数包括孔隙配位数（每个孔隙连接的喉道数量）、欧拉数（表征孔隙空间的连通性整体水平，数值越低连通性越好）、连通孔隙体积占比等。这类参数能够直接反映孔隙结构的拓扑特征，但需要借助CT扫描、核磁共振等成像技术获取三维孔隙结构，试验成本较高，难以应用于大尺度野外样品分析。基于SWRC的间接推导法：利用不同吸力段SWRC的形态特征反演连通性，典型代表是孔隙连通性指数l（即Mualem模型中的曲率参数）。Mualem假设l的默认值为0.5，但大量试验表明该值在不同土壤中变异范围可达-5到10，其数值差异本质上反映了孔隙连通性的差异。此外，通过SWRC的分形特征也可计算连通性分形维数，表征孔隙空间的自相似连通特征。这类方法无需复杂的成像设备，成本较低，但依赖于SWRC的测量精度，且参数的物理意义仍存在一定争议。**tracer试验法**：通过示踪剂的穿透曲线计算有效孔隙度与可动水比例，间接表征水力连通性。当孔隙连通性较好时，示踪剂穿透速度快，穿透曲线峰值高且拖尾短；若连通性差，示踪剂易被滞留于死端孔隙中，穿透曲线呈现明显的拖尾特征。这类方法能够直接反映实际水流过程中的连通性，但无法区分几何连通性与水力连通性的贡献，且试验周期较长。二、传统非饱和渗透系数预测模型的连通性假设缺陷当前主流的非饱和渗透系数预测模型均建立在“等效毛细管束”假设基础上，将土壤孔隙简化为一组互不连通、平行排列的毛细管，认为水力传导能力仅由孔隙尺寸分布决定，完全忽略了孔隙之间的连接关系。这类假设在饱和状态下偏差较小，因为几乎所有连通孔隙都被水填充，流动路径的连续性主要受大孔隙控制；但在非饱和状态下，该假设的局限性被显著放大，直接导致预测偏差。以应用最广泛的vanGenuchten-Mualem（VGM）模型为例，其非饱和渗透系数计算公式为：$$K(\theta)=K_s\left(\frac{\theta-\theta_r}{\theta_s-\theta_r}\right)^l\left[1-\left(1-\left(\frac{\theta-\theta_r}{\theta_s-\theta_r}\right)^{1/m}\right)^m\right]^2$$其中$K_s$为饱和渗透系数，$\theta$为体积含水率，$\theta_s$为饱和含水率，$\theta_r$为残余含水率，$m$为vanGenuchten模型的形状参数，$l$为孔隙连通性指数。传统应用中通常将$l$取为固定值0.5，该值来源于对砂质土壤的拟合结果，但在黏粒含量较高、团聚体发育或存在大孔隙的土壤中，这一假设会带来显著偏差：对于团聚化黑土，由于团聚体内孔隙与团聚体间孔隙形成两级连通网络，当含水率处于田间持水量附近时，大孔隙排空但团聚体内孔隙仍保持连通，实测渗透系数比VGM模型预测值高1~2个数量级，若将$l$调整为2可显著降低偏差，此时参数$l$的增大本质上反映了小孔隙网络的高连通性。对于存在生物孔隙或裂隙的林地土壤，SWRC在低吸力段（0~100kPa）呈现明显的台阶状特征，表明大孔隙占比高但连通性存在异质性，传统模型假设孔隙均匀分布，无法反映大孔隙之间的连通性差异，预测结果通常比实测值低2~3个数量级。对于压实的农业土壤，机械压实导致孔隙被挤压，连通性显著降低，死端孔隙占比增加，此时实测渗透系数远低于模型预测值，$l$值甚至会出现负值，反映出水分流动路径的高度不连续性。另一类常用的Brooks-Corey模型同样存在类似问题，其假设孔隙尺寸分布服从幂函数关系，且流动路径的弯曲度为常数，没有考虑随着含水率降低，连通性下降导致的弯曲度增大效应。研究表明，当土壤含水率从饱和降低到残余含水率时，实际流动路径的弯曲度可增大2~5倍，而传统模型假设弯曲度不变，直接导致低含水率段渗透系数的预测值普遍偏高。三、孔隙连通性影响预测偏差的作用机制孔隙连通性对非饱和渗透系数预测偏差的影响本质上是通过改变水流路径的有效传导面积和弯曲度实现的，其作用机制可分为三个层面：有效传导孔隙的截断效应：传统模型假设只要孔隙尺寸大于对应吸力下的阈值，就会参与导水，但实际上只有属于连续水相网络的孔隙才具有传导能力。当基质势升高，部分原本连通的孔隙被空气填充，导致连续水相网络被截断，有效传导孔隙的占比远低于模型假设的尺寸分布占比。例如在砂壤土中，当基质势为30kPa时，传统模型假设有30%的孔隙参与导水，但实际通过核磁共振观测到的连续水相孔隙占比仅为15%，直接导致模型预测值是实测值的2倍左右。这种截断效应在低含水率段尤为显著，当含水率接近残余含水率时，连续水相网络基本破碎，水分仅通过孤立的孔隙膜流动，此时模型预测值与实测值的偏差可达3个数量级以上。流动路径的弯曲度动态变化：传统模型通常将弯曲度设为常数，或仅考虑与含水率的线性关系，但实际上弯曲度随孔隙连通性的变化呈现非线性特征。当孔隙连通性较好时，水流可以沿着最短路径流动，弯曲度较低；随着含水率降低，连通性下降，水流需要绕过被空气占据的孔隙，流动路径变得更加曲折，弯曲度急剧增大。研究表明，对于壤质黏土，当含水率从饱和降低到0.2cm³/cm³时，弯曲度从1.3增大到4.7，而传统模型假设弯曲度为2，导致低含水率下渗透系数预测值比实测值高1~2个数量级。不同质地土壤的弯曲度变化规律存在显著差异：砂质土壤孔隙分布均匀，连通性变化平缓，弯曲度随含水率的变化幅度较小；黏质土壤孔隙结构复杂，连通性随含水率变化剧烈，弯曲度的非线性特征更加明显。孔隙连通性的空间异质性效应：自然土壤的孔隙连通性存在显著的空间异质性，尤其是在剖面尺度上，耕作层、犁底层、淀积层的孔隙连通性差异巨大，而传统模型通常假设土壤是均质的，仅采用单一的SWRC代表整个土层的特征，忽略了连通性的空间变异。例如在华北平原的农田剖面中，耕作层由于受农事活动影响，大孔隙发育且连通性好，而犁底层经过长期压实，孔隙连通性仅为耕作层的1/5，若采用耕作层的SWRC预测整个剖面的渗透系数，会导致犁底层的渗透系数被高估2个数量级，最终模拟的地下水补给量比实际观测值高40%以上。此外，团聚体尺度的连通性异质性也会带来偏差：团聚体内的小孔隙连通性好，而团聚体之间的大孔隙在低吸力下排空后会阻断水流路径，传统模型将两类孔隙视为统一的尺寸分布，无法反映这种两级网络的连通性差异，导致预测结果在田间持水量附近出现显著偏差。四、考虑孔隙连通性的预测模型改进方向为降低非饱和渗透系数的预测偏差，近年来学者们围绕孔隙连通性的表征与融入开展了大量研究，主要改进方向包括以下三类：连通性参数的本地化率定：针对传统模型中连通性指数$l$取固定值的缺陷，大量研究通过实测数据建立$l$值与土壤基本属性的定量关系，实现参数的本地化估算。例如基于全球1000多组土壤数据的统计分析表明，$l$值与黏粒含量呈显著负相关，与砂粒含量呈显著正相关，据此建立的经验公式可将VGM模型的预测偏差降低30%~50%。此外，也有研究通过SWRC的形状参数估算$l$值，例如利用vanGenuchten模型的进气值和斜率参数反演孔隙连通性，无需额外的渗透试验即可实现参数估算，适用于缺乏实测渗透数据的区域。双孔隙/多孔隙模型的构建：针对具有两级孔隙结构的团聚化土壤或裂隙土壤，将孔隙系统分为两个或多个相互独立又存在联系的子系统，分别表征不同尺度孔隙的连通性特征。典型的双孔隙模型将孔隙分为基质孔隙（团聚体内的小孔隙）和结构孔隙（团聚体间的大孔隙或裂隙），分别构建两类孔隙的SWRC和渗透系数公式，考虑两者之间的水分交换和连通性差异。研究表明，这类模型在黏粒含量大于30%的团聚化土壤中，预测精度比传统单孔隙模型提高1~2个数量级，尤其在中低含水率段的偏差降低最为显著。基于孔隙网络模型的直接模拟：借助CT扫描技术获取土壤三维孔隙结构，构建真实的孔隙网络模型，直接模拟不同基质势下的水分分布和流动过程，无需依赖等效毛细管束假设。这类模型能够explicit考虑孔隙的几何形状、配位数、连通性等拓扑特征，预测结果与实测值的偏差通常小于50%。但目前该方法的应用仍受限于成像分辨率和计算成本，难以应用于大尺度样品，未来随着计算机性能的提升和成像技术的普及，有望成为非饱和渗透系数预测的重要手段。五、当前研究存在的不足与未来展望尽管围绕孔隙连通性的研究已取得显著进展，但当前仍存在几个关键问题亟待解决：一是孔隙连通性的定量表征仍缺乏统一标准，不同方法获得的连通性参数之间缺乏可比性，如何建立几何连通性与水力连通性之间的定量转换关系，仍是未来研究的核心难点。二是现有改进模型大多针对特定土壤类型或特定含水率段，普适性不足，如何将孔隙连通性的影响融入到统一的理论框架中，构建适用于不同质地、不同结构土壤的通用预测模型，仍是领域内的重要挑战。三是大尺度上孔隙连通性的空间变异特征仍不明确，如何将点尺度的连通性观测结果扩展到田间或流域尺度，解决非饱和渗流模拟的尺度效应问题，是实际应用中