本文发表于《水科学进展》期刊，网络首发时间为2026年1月30日。研究由河海大学的盛奕华、李致家、刘志雨与水利部海委漳卫南运河管理局的徐宁合作完成。该研究针对水库入库洪水预报中的一个关键难题——库区干支流面积悬殊导致建模时存在尺度效应与计算成本矛盾，提出并验证了一种基于TOKASIDE模型的空间多尺度嵌套预报方法，旨在提升无实测资料库区的洪水预报精度。

研究的学术背景聚焦于水文模拟与洪水预报领域。水库的防洪安全与优化调度高度依赖于高精度的入库洪水预报。入库洪水通常由上游干流来水、库区区间支流与坡面来水以及库面直接降雨径流三部分构成。其中，库区区间往往是缺乏水文观测资料的“无资料地区”，这使得传统水文模型难以通过率定来准确确定其产汇流参数，成为预报的核心难点。此外，无资料库区通常地形复杂，干支流汇水面积差异巨大。若采用统一的高分辨率网格进行全域模拟，虽然能精细刻画地形和河网，但计算成本高昂，难以满足实时预报需求；若采用较低的粗分辨率，又会平滑掉关键的地形细节，导致产汇流过程失真，降低模拟精度。目前，虽有研究利用分布式水文模型对无资料库区进行整体参数率定，但所得参数是混合了干流演算与区间产汇流的“综合参数”，物理意义不清，且在不同水文条件下的适应性可能降低。多尺度嵌套模拟技术被认为是解决上述矛盾的可行途径，该技术在水动力学和城市洪涝模拟中已有应用，但在流域产汇流过程模拟，特别是应用于无资料库区洪水预报的研究和实践相对较少。

本研究详细的工作流程可以概括为理论解析、方法构建、验证与应用三个核心阶段。

第一阶段：理论解析与尺度效应量化。 此阶段旨在深入理解TOKASIDE模型模拟结果对空间分辨率的依赖关系，为后续的参数尺度转换提供理论基础。研究选择了一个有实测资料的流域——漳河观台区间（面积1204平方公里）作为“试验场”。首先，研究者准备了数字高程（DEM）、土地覆盖和土壤类型等分布式输入数据，并将DEM重采样为12个不同的空间分辨率（从100米到2000米）。接着，进行了一系列数值实验：1) 分析输入数据的尺度依赖性：借鉴分形理论和地貌瞬时单位线（GIUH）理论，计算了不同分辨率下河网的分形维数，并统计了所有网格单元的流路长度与坡度分布。结果显示，河网分形维数在100-500米的高分辨率区间稳定在1.17，在600-2000米的较低分辨率区间则稳定在1.05，揭示了流域地貌结构存在两个明显的“标度不变性区间”。流路长度和坡度的概率密度分布也随分辨率变化呈现规律性偏移。2) 量化模型输出的尺度效应：以1000米分辨率下率定好的模型参数为基准，直接将其移植到其他11个分辨率下进行洪水模拟。结果发现，模拟流量对分辨率变化呈现出系统性偏差：在100-500米分辨率区间，模拟流量系统性偏低，洪峰平均减小21.3%，过程线坦化；而在600-2000米分辨率区间，模拟流量系统性偏高，洪峰平均增大2.6%。这种效应在洪峰流量大于1500立方米/秒的大洪水事件中尤为显著，其纳什效率系数波动剧烈，而在小洪水事件中则不敏感。进一步分析表明，大洪水中地表径流占比高，对汇流路径（受分辨率影响）极其敏感，从而放大了尺度效应。

第二阶段：构建参数尺度转换方案。 基于第一阶段的发现，研究指出分辨率变化主要改变了地形特征（流路长度、坡度）的表达，从而影响汇流速度。为了在不同分辨率下保持水动力过程的相似性，需要对表征水流阻力的关键参数——地表曼宁糙率（ns）和河道曼宁糙率（nch）——进行补偿性调整。研究者基于GIUH理论，提出了一个参数尺度转换方法。其核心是定义一个表征汇流能力的尺度因子（K），该因子与流域内所有网格单元的流路长度之和成正比，与流域平均坡度成反比。转换流程为：首先确定一个参考分辨率（本研究选1000米）及其对应的已率定参数；然后计算目标分辨率下的尺度因子K_target与参考分辨率下的K_ref的比值；最后，根据此比值对曼宁糙率参数进行缩放（公式：n_target = n_ref × (K_target / K_ref)）。该方案在分形维数标度不变的区间内分别进行。

第三阶段：方法验证与无资料库区应用。 此阶段旨在验证尺度转换方案的有效性，并将其应用于真正的无资料库区。首先，在观台区间验证尺度转换方案。选取100、200、400、800米四个目标分辨率（覆盖了后续多尺度建模的关键跨度），应用上述转换方案调整其糙率参数。结果显示，转换后，所有分辨率下的模拟洪水过程线都紧密簇拥在1000米基准线周围，洪峰形状和位相高度一致。定量评估表明，尺度转换使平均洪峰相对误差绝对值降低了13.9%，平均纳什效率系数提高了8.5%，特别是在大洪水的洪峰校正上效果突出。验证成功后，研究将整套空间多尺度嵌套模拟技术应用于无资料的岳城水库库区（面积290.5平方公里）。该库区包含多条无水文站的支沟，在暴雨期间产流显著。嵌套建模的策略是：依据各支流流域的面积和地形复杂度，为其分配不同的网格分辨率（例如，面积小但重要的观台沟用100米，面积较大的都党沟用800米），并通过一个统一的400米分辨率模型进行汇流衔接。将经过观台区间验证的尺度转换关系应用于确定各分辨率子区域的参数，从而构建了库区的多尺度嵌套预报模型。利用2016年和2021年的两场历史洪水进行检验，并与传统的全域1000米单一分辨率模拟结果进行对比。

研究取得了一系列明确的结果。在尺度效应分析方面，证实了流域地貌数字化表达和TOKASIDE模型模拟结果均具有强烈的尺度依赖性，且这种依赖性在大洪水事件中更为显著，这为针对性地开展尺度转换研究提供了依据。在参数尺度转换验证方面，基于GIUH理论构建的曼宁糙率转换方法被证明是有效的，它能显著降低不同分辨率模型间的系统性偏差，使得移植自参考分辨率的参数能在目标分辨率下产生一致的水文响应。最重要的成果体现在岳城水库库区的应用上。多尺度嵌套模拟方案成功复现了两场历史洪水的入库过程。与单一分辨率（1000米）模拟相比，多尺度方案在保持整体过程拟合优度（两场洪水纳什效率系数均超过0.7）的同时，显著提升了洪峰预报精度：对于“2016071906”号洪水，洪峰误差从-45%改善至-3.57%；对于“2021072104”号洪水，误差从-42.98%改善至-30.5%。两场洪水平均，洪峰模拟误差降低了约27%。这一改进对于依赖准确洪峰进行调度的水库防洪安全至关重要。

本研究的结论是，针对无资料库区入库洪水预报中因尺度效应和参数确定困难导致的精度难题，所发展的基于TOKASIDE模型的空间多尺度嵌套模拟技术是有效且可行的。该方法通过解析尺度效应、建立参数跨尺度转换关系，实现了对库区不同子流域“因地制宜”地采用不同分辨率建模，在不过度增加计算成本的前提下，显著提高了洪水预报，特别是洪峰预报的精度与可靠性。

本研究的亮点和创新之处在于：第一，科学问题聚焦且有深度：不仅指出了无资料库区洪水预报的普遍难题，更深入剖析了其内在矛盾——尺度效应与计算效率的平衡，并提出了创新的解决方案。第二，理论扎实，方法严谨：研究并非简单提出嵌套思路，而是从分形理论和地貌学理论出发，系统解析了模型尺度效应的成因，并在此基础上建立了具有物理基础的参数尺度转换关系，使整个方法体系有据可依。第三，验证与应用逻辑清晰：研究设计巧妙，先在有资料流域（观台区间）进行理论探索和方法验证，再将成熟的技术体系应用于无资料库区（岳城水库），形成了从理论到方法再到实践的完整闭环。第四，应用价值显著：所提出的多尺度嵌套方案无需在无资料区重新率定参数，大幅节省了建模时间和成本，同时显著提升了预报精度，为水库实际的防洪调度决策提供了直接、可靠的技术支撑。

此外，研究过程中还有一些有价值的发现。例如，研究明确指出尺度效应的显著程度与产流机制密切相关：以地表径流为主导的大洪水对汇流路径变化（由分辨率引起）极其敏感，因而尺度效应明显；而以壤中流为主导的小洪水则对尺度变化不敏感。这一认识有助于指导研究者将尺度研究的重点放在大洪水事件上，提升了研究的针对性和效率。同时，研究中对河网分形维数在不同分辨率区间呈现稳定值的发现，也为分尺度建立参数转换关系提供了重要的理论依据。总的来说，这项工作为分布式水文模型在复杂无资料地区的应用提供了一个可复制、可推广的高效精准框架。