土壤水分非均匀运动的分形模型与优先流通道表征研究综述

土壤水分非均匀运动的分形模型与优先流通道表征研究综述土壤是由大小不一的颗粒、孔隙和团聚体共同构成的复杂多孔介质，其内部结构的空间异质性直接决定了水分入渗、溶质运移等过程的非均匀特征。传统的基于达西定律的均匀流假设往往低估了水分在土壤中的实际运移速度，也无法解释田间观测中普遍存在的溶质快速穿透、水分分布不连续等现象。优先流作为土壤水分非均匀运动的典型表现形式，是指水分绕过部分土壤基质，沿着大孔隙、裂缝、根孔等优势通道快速运移的过程，其形成机制与土壤孔隙结构的分形特征高度相关。分形理论作为描述不规则、自相似结构的有效工具，为量化土壤孔隙系统的复杂特征、构建非均匀水流运动模型提供了新的研究路径。一、土壤水分非均匀运动的核心特征与驱动机制土壤水分非均匀运动的核心特征表现为入渗过程中湿润锋的不规则跳跃、水分分布的斑块化以及运移速率的多尺度差异性。在田间尺度下，即使是质地看似均一的壤土，水分入渗深度的变异系数也可达到30%以上，而在存在大量生物孔隙的林地土壤中，这一数值甚至超过60%。这种非均匀性直接导致传统的理查兹方程模拟结果与实测值出现显著偏差，尤其是在短历时强降雨、灌溉等快速入渗场景下，均匀流模型的预测误差可达2~3个数量级。驱动土壤水分非均匀运动的因素可分为内在结构因素和外在环境因素两类。内在因素中，土壤孔隙系统的多尺度分布是最核心的控制条件：直径大于0.1mm的大孔隙虽然通常仅占土壤总孔隙度的5%~20%，却承担了60%以上的快速水分通量；团聚体之间的裂隙、动物活动形成的通道、植物根系腐烂后留下的孔道等结构单元，共同构成了优先流的优势路径，其形态和连通性具有明显的自相似特征。外在因素方面，降雨强度超过土壤基质入渗能力会引发大孔隙优先流启动，土壤初始含水率的差异会改变水分的流动路径选择，而土壤干湿交替、冻融循环等过程会进一步改造孔隙结构，加剧水分运动的非均匀程度。二、分形理论在土壤孔隙结构表征中的应用进展分形理论的核心是识别复杂系统在不同尺度下的自相似规律，通过分形维数量化结构的不规则程度。土壤作为自然风化和生物活动共同作用的产物，其颗粒组成、孔隙分布均表现出典型的分形特征，这为分形模型的应用提供了物理基础。早期的分形研究主要聚焦于土壤颗粒的粒径分布特征，通过质量分形维数描述颗粒大小的异质性。研究表明，砂质土壤的颗粒分形维数通常在2.2~2.5之间，而黏质土壤可达到2.7~2.9，分形维数越高意味着细颗粒含量越高，孔隙结构越复杂。随着成像技术的发展，CT扫描、核磁共振等技术开始被用于获取土壤内部三维孔隙结构，学者们进一步提出了孔隙表面分形维数、孔隙体积分形维数等参数，实现了对孔隙形态、连通性的量化表征。其中，孔隙分形维数与土壤饱和导水率存在显著的幂函数关系，相较于传统的经验公式，基于分形维数的导水率预测模型在不同质地土壤中的适用性提高了40%以上。针对孔隙结构的多重分形特征，后续研究逐步突破了单一分形维数的局限。多重分形谱能够刻画不同尺度孔隙的局部异质性，其谱宽参数可以反映孔隙大小的分布范围，谱的不对称性则对应了大孔隙与小孔隙的相对占比。在优先流发育的土壤中，多重分形谱的宽度通常比均匀流土壤大20%以上，且谱形呈现明显的右偏特征，表明大孔隙对整体孔隙系统的贡献更为突出。这些研究成果为从结构角度解释水分非均匀运动提供了关键参数。三、基于分形特征的非均匀水流运动模型构建传统的水流运动模型假设土壤水力参数是均质的，无法反映优先流的非线性特征，而分形模型通过将孔隙结构的分形规律嵌入水力参数表达式，实现了对非均匀运动过程的定量描述。目前主流的分形水流模型可分为三类：孔隙尺度的分形网络模型、宏观尺度的分形理查兹方程，以及考虑多尺度耦合的双重分形模型。孔隙尺度的分形网络模型将土壤孔隙系统抽象为具有自相似特征的分叉网络，通过设定不同级别孔隙的直径、长度和连通性规则，模拟水分在孔隙网络中的流动路径。这类模型能够复现湿润锋的指状突进现象，并且可以解释为什么当入渗通量超过临界值时，优先流的发生概率会呈指数增长。研究显示，当孔隙网络的分形维数在2.6~2.8之间时，模拟得到的水分穿透曲线与林地土壤的实测结果拟合度最高，其决定系数R²普遍大于0.85。宏观尺度的分形理查兹方程对传统模型的改进主要体现在水力参数的分形表达上。基于孔隙分形特征推导得到的土壤水分特征曲线和非饱和导水率公式，不再依赖于经验性的vanGenuchten参数，而是直接将分形维数作为核心输入参数。与传统模型相比，分形导水率模型在低含水率段的预测精度提升尤为明显，能够更好地描述干旱条件下水分沿细小孔隙的缓慢运移过程。针对优先流的快速运动特征，部分学者进一步提出了活动流场模型，通过引入分形动态的活动流场占比参数，实现了对非均匀流场随时间演化过程的动态模拟。双重分形模型则同时考虑了基质域和大孔隙域的分形特征，将两个域的水分运动过程通过质量交换项耦合。其中基质域的孔隙分形维数较低，水流运动符合基本的达西规律，而大孔隙域的分形维数较高，水流以紊流为主。这类模型在农业面源污染模拟中得到了广泛应用，相较于传统的双重孔隙度模型，其对农药、氮肥等溶质快速穿透的预测误差降低了35%以上，能够更准确地评估污染物向地下水的淋失风险。四、优先流通道的定量表征方法体系优先流通道的形态、数量和连通性直接决定了非均匀水流的运移效率，其定量表征是分形模型应用的前提。目前的表征方法可分为原位观测法、染色示踪法和数值模拟法三类，不同方法适用于不同的研究尺度。原位观测法主要包括土壤切片技术、CT扫描和探地雷达探测。土壤切片技术能够获取厘米尺度的孔隙结构参数，通过图像识别可以直接统计优先流通道的直径、密度和弯曲度，但其缺点是需要破坏原状土壤，难以实现连续观测。工业CT扫描的分辨率可达到微米级别，能够重建土壤三维孔隙网络，直接计算通道的连通性分形维数，目前已成为实验室尺度优先流研究的标准手段。田间尺度的优先流通道探测则多采用探地雷达，通过解析雷达信号的振幅和频谱特征，可以反演2米深度范围内的大孔隙分布，其分形维数计算结果与开挖实测值的误差通常小于10%。染色示踪法是识别优先流路径最直观的方法，通过将亮蓝、罗丹明等染色剂随入渗溶液施加到土壤中，然后开挖剖面记录染色区域的分布特征。基于染色图像计算的分形维数可以量化优先流的发育程度：染色区域分形维数小于1.5时表明水流以均匀基质流为主，介于1.5~1.8之间时为基质流与优先流共存，大于1.8时则表明优先流占据主导地位。为了提高示踪的定量化水平，近年来学者们开发了数字图像处理技术，能够自动识别染色区域的连通斑块，计算优先流通道的数量、平均直径和连通性指数，实现了对染色模式的快速批量分析。数值模拟法主要通过构建虚拟孔隙结构来表征优先流通道的统计特征。基于分形理论的孔隙生长模型可以根据实测的分形维数生成与真实土壤结构高度相似的虚拟介质，在此基础上开展水流运动模拟，能够定量分析不同通道特征对水分运移的影响。模拟结果表明，优先流通道的连通性对水分通量的贡献远大于通道直径，当连通性分形维数提高0.2时，相同水头下的入渗通量可提升1~2倍，这一结论为农田耕作措施改良提供了理论依据——通过深松、免耕等措施增加大孔隙的连通性，能够有效提高土壤的入渗能力，减少地表径流。五、现有研究的不足与未来发展方向尽管分形模型在土壤水分非均匀运动研究中已经取得了显著进展，但仍存在三个关键问题有待解决：一是多重分形参数的物理意义尚不明确，不同研究中得到的分形维数与水力参数的关系差异较大，缺乏统一的换算标准；二是现有的分形模型大多假设孔隙结构是静态的，没有考虑生物活动、干湿交替等过程对孔隙结构的动态改造，难以模拟长期耕作条件下的优先流演化规律；三是尺度转换问题尚未突破，实验室厘米尺度得到的分形参数无法直接应用于田间百米尺度的水文模拟，不同尺度下分形特征的耦合关系仍不清晰。未来研究的重点方向包括三个方面：首先，开展多方法融合的优先流通道表征研究，将CT扫描的高精度结构信息、染色示踪的动态流动信息和田间原位监测的通量信息相结合，建立结构特征与运动过程之间的定量关系，明确分形参数的物理意义和适用范围。其次，开发动态分形演化模型，将土壤动物活动、根系生长、干湿循环等过程纳入模型框架，实现对优先流通道形成、发育和消亡过程的动态模拟。最后，构建多尺度分形耦合模型，