土壤水分特征曲线中的滞后现象对非饱和渗流模拟的影响研究综述

土壤水分特征曲线中的滞后现象对非饱和渗流模拟的影响研究综述土壤水分特征曲线（SoilWaterCharacteristicCurve,SWCC）作为描述非饱和土壤中基质势与体积含水量之间非线性关系的核心本构曲线，是水文过程模拟、农业水土工程及地质灾害防控等领域的基础参数。在土壤吸湿（含水量由低到高）与脱湿（含水量由高到低）过程中，相同基质势对应的含水量存在显著差异，该现象被定义为SWCC的滞后效应。长期以来，多数非饱和渗流模型为简化计算，默认采用单一脱湿曲线表征土壤水力特性，但大量野外试验与数值模拟结果表明，滞后现象的忽略会导致渗流场演化规律的预测偏差可达30%以上，直接影响边坡稳定性评估、农田灌溉效率测算及地下水补给量估算等应用场景的结果可靠性。滞后现象的物理机制与定量表征方法滞后现象的产生是土壤孔隙结构、流体特性与界面作用共同驱动的结果，其物理机制可归纳为三个核心维度。首先是“墨水瓶效应”：土壤孔隙普遍存在“狭颈宽腹”的非规则形态，脱湿过程中宽腹孔隙中的水分需克服狭颈处的更大毛细张力才能排出，导致相同基质势下脱湿过程的含水量高于吸湿过程，该机制在砂性土中对滞后效应的贡献占比可达60%以上。其次是接触角滞后：土壤颗粒表面的接触角在吸湿过程中（液体推进）通常大于脱湿过程中（液体后退），相同毛细压力下接触角越大，孔隙所能容纳的水量越高，这一效应在黏粒含量较高的土壤中表现尤为显著，当黏粒含量超过30%时，接触角滞后对总滞后量的贡献可提升至40%。最后是孔隙封闭效应：吸湿过程中气体被快速入渗的水分包裹形成封闭气泡，占据部分孔隙空间，导致相同基质势下吸湿含水量低于脱湿过程，封闭气泡的体积占比通常在2%~15%区间波动，与土壤初始含水率、入渗速率直接相关。针对滞后现象的定量表征，目前主流的建模方法可分为三类。经验域模型以Kool-Parker模型为代表，通过引入形状参数分别拟合吸湿与脱湿曲线，再通过插值函数计算干湿循环路径下的扫描曲线，该类模型参数少、计算效率高，在农田水文模拟中应用广泛，但缺乏物理机制支撑，对复杂干湿循环路径的适应性较差。物理域模型以孔隙网络模型为核心，通过重构土壤三维孔隙结构，模拟不同干湿路径下孔隙的充排水过程，能够从微观层面解释滞后现象的演化规律，但计算成本极高，仅适用于厘米级尺度的理想土样模拟。半经验半物理模型以Scott等提出的独立域模型为代表，将土壤孔隙划分为若干独立的孔隙组，每个孔隙组对应固定的吸湿与脱湿临界基质势，通过统计孔隙组的分布特征实现扫描曲线的预测，该类模型兼顾计算效率与物理合理性，是当前工程应用领域的主流选择。不同表征方法的适用场景差异显著：经验模型适合长时序大尺度的流域水文模拟，物理模型适合微观机制研究，半经验半物理模型则更适配边坡渗流、污染物迁移等中等尺度的工程问题。滞后效应对非饱和渗流场演化规律的影响非饱和渗流模拟的核心是求解Richards方程，而SWCC作为方程的关键输入项，其滞后特性首先会改变基质势与含水量的映射关系，进而影响渗流场的时空分布特征。在降雨入渗场景中，若忽略滞后效应，模拟得到的湿润锋推进速度会明显高于实测值：以砂壤土为例，当降雨强度为10mm/h、持时24h时，采用单一脱湿曲线模拟的湿润锋深度可达72cm，而考虑滞后效应的模拟结果仅为58cm，偏差达24%。其原因在于，吸湿过程中相同基质势对应的含水量更低，土壤的实际持水能力弱于脱湿曲线的预测值，更多水分以地表径流形式流失而非入渗到深层土壤，该偏差在间歇性降雨场景中会进一步放大——当存在多轮干湿循环时，忽略滞后效应的入渗量累积偏差可超过40%。滞后效应还会显著影响土壤孔隙水压力的动态响应过程。在蒸发条件下，表层土壤含水量逐渐降低进入脱湿路径，若采用单一吸湿曲线进行模拟，会高估相同含水量下的基质势绝对值，导致模拟的表层孔隙水压力偏高，与实际监测数据的偏差可达20~50kPa。对于层状土壤剖面，滞后效应的影响存在明显的层间差异：砂质土层的滞后环宽度较小但变化速率快，黏质土层的滞后环更宽且对长期干湿循环的记忆效应更强，两种土层交界面处的水力梯度模拟结果受滞后效应的影响最为显著，部分场景下的模拟偏差甚至会改变渗流方向，将原本的向下入渗错误模拟为向上毛管水上升。此外，滞后现象会改变非饱和渗流的动态平衡特征。在周期性干湿交替条件下，考虑滞后效应的土壤含水量会呈现出明显的“滞回环”特征，即经过多轮循环后含水量不会回到初始值，而是在一个固定区间内波动，且波动幅度与滞后环的宽度正相关。若忽略滞后效应，模拟结果会错误地认为含水量最终会收敛到一个稳定值，导致对土壤水库调蓄能力的评估出现偏差。在地下水浅埋区，滞后效应还会影响降水入渗对地下水的补给过程：研究表明，忽略滞后效应会使得地下水补给量的模拟结果偏高15%~30%，且土壤质地越细，偏差越大。滞后效应在不同应用场景中的影响差异在农业水土工程领域，滞后效应直接影响农田灌溉制度的制定与水分利用效率的评估。当采用喷灌等间歇性灌溉方式时，土壤处于反复的干湿循环过程中，若忽略滞后效应，会高估土壤的保水能力，导致灌溉定额的设计值偏高，不仅造成水资源浪费，还可能引发深层渗漏与肥料淋失风险。华北平原冬小麦种植区的试验数据表明，考虑滞后效应的灌溉制度可减少12%的灌溉水量，同时作物产量不受影响，而基于单一脱湿曲线制定的灌溉方案会导致20%以上的水分通过深层渗漏损失。对于盐碱地改良场景，滞后效应会影响盐分的运移路径：脱湿过程中水分向上运动携带盐分至表层，吸湿过程中水分向下入渗淋洗盐分，忽略滞后效应会导致盐分累积量的模拟偏差可达40%，直接影响改良方案的效果评估。在地质灾害防控领域，滞后效应对边坡稳定性的影响尤为突出。降雨入渗过程中，坡体表层土壤处于吸湿路径，基质势降低速率慢于单一脱湿曲线的预测结果，使得浅层土体的有效应力维持在较高水平，若忽略滞后效应，会过早预测边坡失稳，导致预警准确率降低。三峡库区某堆积层边坡的监测数据显示，2023年雨季期间，采用考虑滞后效应的渗流模型模拟的浅层孔隙水压力与实测值的平均误差仅为5.8kPa，而忽略滞后效应的模拟误差达31.2kPa，对应的边坡安全系数计算偏差超过15%。在极端强降雨场景下，滞后效应还会延缓孔隙水压力的消散过程，当降雨停止后，考虑滞后效应的边坡安全系数回升速度比不考虑滞后的结果慢2~3天，这一差异对于雨后边坡的风险管控具有重要指导意义。在水文水资源领域，滞后效应是流域水文模型参数不确定性的重要来源。针对江淮丘陵区小流域的模拟研究表明，将滞后效应纳入SWCC表征后，流域径流深的模拟效率系数（NSE）从0.67提升至0.82，洪峰流量的模拟偏差从21%降至8%。对于干旱半干旱地区的地下水补给评估，滞后效应的影响更为显著：由于该地区的干湿循环周期长、幅度大，土壤的滞后记忆效应可延续数月甚至数年，忽略滞后效应会导致地下水补给量的年际波动特征被严重低估，无法准确评估气候变化对地下水资源的影响。内蒙古草原区的试验结果显示，考虑滞后效应的地下水补给量模拟结果与同位素示踪法的实测结果偏差仅为7%，而传统模型的偏差达34%。滞后效应纳入渗流模拟的挑战与研究趋势尽管滞后效应的重要性已得到广泛认可，但将其全面纳入工程尺度的非饱和渗流模拟仍面临三大挑战。首先是参数获取难度大：传统的SWCC试验仅测定脱湿曲线，而获取完整的滞后曲线需要进行多轮干湿循环试验，试验周期从一周到数月不等，且对试验设备的控制精度要求极高，工程勘察阶段通常难以提供足够的实测参数。其次是计算复杂度提升：考虑滞后效应后，Richards方程的求解需要额外追踪每个计算单元的干湿路径与历史状态，三维大尺度模型的计算量会增加2~5倍，对于流域尺度的长时序模拟而言，计算成本难以承受。最后是尺度效应问题：目前滞后效应的表征参数多基于室内厘米级土样试验获得，而野外米级至千米级尺度的土壤存在显著的空间异质性，微观滞后效应在宏观尺度上的等效表征机制尚不明确，直接套用室内参数会带来新的误差。针对上述挑战，当前相关研究主要沿着三个方向推进。一是滞后参数的快速测定方法研究：通过时域反射仪（TDR）、张力计等原位监测设备的组合应用，结合机器学习算法反演SWCC的滞后参数，将参数获取周期从数月缩短至数天，目前该方法在均质土壤中的反演精度已达85%以上。二是高效数值算法的开发：通过预处理技术、自适应时间步长算法与并行计算框架的结合，将考虑滞后效应的渗流模型计算效率提升3~4倍，部分商业水文软件（如HYDRUS、MODFLOW）已集成简化的滞后模块，可支撑中等尺度的工程应用。三是多尺度滞后表征框架的构建：采用孔隙尺度-土柱尺度-坡面尺度的跨尺度试验，结合尺度提升方法，建立考虑空间异质性的宏观滞后参数计算模型，目前已在砂质土壤区域取得阶段性成果，未来将逐步拓展到黏质土与复合土壤类型。从长期发展来看，滞后效应与其他非饱和土水力过程的耦合机制研究是未来的重点方向。一方面，需