《Earth's Future》：Small Reservoirs Offer a New Perspective on Flood Reduction in Large Basins
编辑推荐：
摘要
大型河流网络连续体内单个水库的洪水削减潜力，因水库特征与网络属性之间的复杂相互作用而仍未得到充分理解。本研究探讨了一系列相对较小的水库是否能在调节下游洪水中发挥重要作用，并评估相较于传统已知的水库特征，水库网络属性如何影响这一作用。本研究的独特贡献在于
摘要
大型河流网络连续体内单个水库的洪水削减潜力，因水库特征与网络属性之间的复杂相互作用而仍未得到充分理解。本研究探讨了一系列相对较小的水库是否能在调节下游洪水中发挥重要作用，并评估相较于传统已知的水库特征，水库网络属性如何影响这一作用。本研究的独特贡献在于：在覆盖美国德克萨斯州墨西哥湾沿岸415,000 km2区域的基于过程的流域水文模型中，集成了81座主要水库（库容介于6 × 10?至8 × 10? m3），并模拟了跨越15,000个河流河段的洪峰流量及对应的下游洪水淹没图。本研究的三项关键结论如下：（a）较小的水库可显著削减下游洪水，这表明仅纳入大型水库的传统洪水风险管理研究方法可能只呈现了部分图景。（b）较小水库的洪水削减效应在上游区域更为有效，不过这一现象可能与干旱程度及整体水资源可用性有关。（c）尽管水库规模至关重要，但它并非决定其下游洪水削减潜力的首要因素；网络属性的影响，如集水区面积、上游水库数量、上游水库的累积最大储水容量以及河流等级（即位置），与水库规模同等重要，且往往更为重要。本研究发现的更广泛影响不仅限于洪水领域，还为应对老化水坝与河流连通性等新兴挑战提供了基础性见解。
小型水库网络在大流域洪水削减中的作用：基于德克萨斯州墨西哥湾沿岸流域的集成模拟研究
1 研究背景与科学问题
随着全球气候变化与人类活动的加剧，河流洪水风险正面临前所未有的严峻挑战。据估算，到本世纪末，全球洪水风险较2010年基线将增长二十倍，潜在年损失可达13.7万亿美元。在此背景下，上游蓄滞洪设施在流域综合管理中的作用日益受到重视，水库作为重要的径流调节基础设施，其削减下游洪水的能力已成为水文水资源领域的核心议题。然而，现有洪水风险模型与规划工作大多聚焦于少数大型水库，对数量众多、分布广泛的小型水库在流域尺度上的洪水调节功能关注不足。这种“重大型、轻小型”的研究范式，主要源于传统观念认为大型水库因其巨大的蓄水容量对区域水安全具有更显著的影响，从而在水资源管理与监管机构中获得了更多优先关注。
然而，从流域系统的整体视角来看，大型河流网络中的水库并非孤立存在，而是通过河网相互连接、级联运行。多座水库的联合行为往往偏离单座水库效应的简单叠加：其相互作用既可能表现为协同增效，也可能表现为相互拮抗。因此，仅关注少数大型水库的研究路径无法完整揭示水库网络在流域尺度上的真实洪水调节能力。与此同时，决定水库洪水削减潜力的因素也远不止传统认知中的蓄水容量和调度运行策略，水库在河网中的空间位置、上下游连接关系等网络属性同样可能发挥关键作用，但这一维度在以往研究中常常被忽略。
新一带卫星遥感技术（如冰、云与陆地高程卫星-2（ICESat-2）和地表水与海洋地形卫星）的发展，使得获取小型湖泊与水库的高精度数据成为可能。然而，尽管数据可用性大幅提升，许多洪水风险研究仍延续着聚焦大型水库的传统路径。在此背景下，一个关键的科学问题亟待回答：在大型流域河网连续体中，相对较小规模的水库群能否在削减下游洪水中发挥显著作用？水库的网络属性相较于其蓄水能力等传统特征，又将对下游洪水削减产生什么样的影响？为回答上述问题，研究选择美国德克萨斯州墨西哥湾沿岸流域作为研究区域，开展了系统的集成模拟研究。
2 研究区概况与数据方法
研究选取美国德克萨斯州墨西哥湾沿岸流域作为研究区域，覆盖面积约415,000 km2。该区域具有三个典型特征：（1）拥有完善的公共水库属性数据，便于集成到过程模型中；（2）是美国受水库调控程度最高的河流系统之一，共有超过7,400座已测绘的水坝与水库，其中122座受到定期监测，每年储存约930万英亩-英尺的地表水；（3）德克萨斯州的突发性洪水和洪水相关死亡人数位居全美前列，使该区域水库的防洪潜力具有重要的区域和全国意义。研究区域内约31%的水库以洪水控制为主要设计用途，其余水库分别用于供水灌溉（25%）、休闲娱乐（20%）以及水力发电、航运和鱼类与野生动物保护等（24%）。
研究采用了基于过程的土壤和水评估工具（SWAT）模型，构建了两类模拟情景：（1）基准情景，模型集成全部81座水库。根据全球湖泊与湿地数据库（GLWD）定义的5 × 10? m3库容阈值，这81座水库包括45座库容范围在6 × 10?至4 × 10? m3之间的小型水库，以及36座库容超过该阈值、最高达8 × 10? m3的大型水库。各水库的蓄水能力（即最大容积）、正常水面面积及出流率分别取自国家水坝清单（USACE）和德克萨斯州水务发展委员会（TWDB）等权威数据库。洪水削减效应以100年一遇设计暴雨为驱动条件进行模拟。（2）水库情景：采用逐座移除法，即将某一座水库的相关参数从模型中移除，模拟无该水库条件下的天然来流状态（即非调节流态），以量化单座水库对下游洪水的边际贡献。该方法的优势在于从现实全水库运行状态出发，避免了多水库同时移除时非线性系统交互对归因分析造成干扰，确保了河-库网络连通性的完整性，并使每一座水库的独立效应均可系统评估。
水文模型构建采用了30 m地形数据重采样后的90 m数字高程模型（DEM），将河网划分为14,829个河段。径流模拟采用SCS曲线数法进行产流估算，Penman-Monteith方程计算蒸散发，变容蓄量法进行河道汇流演算。模型以日为时间步长运行，水库蓄水量平衡遵循特定的水量平衡方程。模型采用2010–2019年1 km分辨率的日尺度气象观测数据进行率定与验证，其中39个水文站点用于率定（2013–2016年），31个站点用于验证（2016–2019年），并采用加权克林-古普塔效率（KGE）作为目标函数衡量模拟精度。此外，采用ResOpsUS水库运行时序数据集对36座水库的出流量模拟结果进行了补充评估。
在洪水淹没制图方面，研究采用高度临近排水（HAND）方法，基于合成水位-流量关系曲线（SRC），以模型生成的100年一遇洪峰流量为上游边界条件，生成90 m分辨率的洪水淹没范围。与传统的二维水动力模型相比，HAND方法计算效率高、可复现性强，且对多水库大区域模拟具有良好的适应性。
为系统解析水库洪水削减潜力的影响因素，研究定义了七大潜在决定因素，分别归属于水库特征和网络属性两个类别：水库特征包括最大储水容量、正常条件下面面积以及平均出流率（反映调度运行策略）；网络属性包括相对集水面积、上游水库数量、上游水库累积最大储水容量以及目标水库所处的河序。由于决定因素间存在高度多重共线性，研究采用岭回归模型配合沙普利加法解释（SHAP）方法，量化了各因素对洪水削减的相对贡献。
3 研究结果与关键发现
3.1 径流模拟精度评估
在70个水文站点中，52个站点的KGE值为正，表明模拟结果总体达到中等至良好精度。大型河流（河序≥4）中29/38个站点取得了合理KGE值，小型上游河流（河序<4）中23/32个站点表现良好。气候湿润的东北部流域（如内切斯-萨宾河与特里尼蒂河流域）模拟精度相对较高，而半干旱的西部和西南部流域（如努埃塞斯河流域）精度较低，这一空间差异与研究区域干旱梯度分布相符。在水库出流验证方面，36个验证位置中有24个KGE为正，表明模拟出流与历史运维数据总体吻合。综合分析表明，模型能够有效模拟湿润或洪水条件下的径流过程，基本再现了洪峰及中高流量条件，但在干旱环境或枯水期表现欠佳。
3.2 小型水库的洪水削减显著性
以位于科罗拉多河流域的特拉维斯湖（该区域第二大水库）和位于布拉索斯河流域的帕洛平托湖为例进行对比分析。特拉维斯湖可将其紧邻下游流域出口的洪峰流量削减约4,900 m3/s，相对于无水库情景削减接近100%，下游洪水淹没面积减少近300 km2。虽然帕洛平托湖的库容仅为特拉维斯湖的约1/17，但其洪峰削减比例仍高达94%（约290 m3/s），淹没面积减少近50 km2。这一结果直观表明，小型水库同样能够在大洪水事件中发挥显著削减作用。
从全流域尺度来看，81座水库中仅有一小部分大型水库展现出显著的洪峰和淹没削减能力，而相当大比例的小型水库也表现出可比的削减潜力。在大型水库中约56%可将洪峰削减50%以上，而小型水库中约60%实现了类似削减水平。深入分析发现，部分具有较大蓄水容量的水库仅表现出边际性的洪水和淹没削减量，说明水库规模（蓄水体积）并不恒定地决定其洪水削减潜力。具有显著洪水削减效果的小型水库主要分布于降水强度较高的北/东北部上游流域，而那些降水强度较低的西部和西南部水库仅表现出名义上的下游削减作用，这体现了大规模降水格局对水库效能的重要影响。
3.3 水库特征与网络属性的差异性影响
对于81座水库中的13座，水库特征组中的高蓄水容量和面积产生了较高的SHAP值，显示出其在削减下游洪水中的较强影响。然而，绝大多数小型水库的SHAP值聚集于负值区，进一步证实：在德克萨斯州墨西哥湾海岸河网背景下，水库规模总体上并非决定其下游洪水削减潜力的强有力因素。另一方面，水库特征组中的出流率变量几乎总是产生较低的SHAP值，表明单座水库的调度运行对其下游洪水削减的影响相对有限。
与之形成鲜明对比的是，网络属性组中75%的水库展现出正的SHAP值，而这些水库主要位于流域上游。不过，不同河序间的SHAP值变化表明，位于主河道上的水库可能比支流上的水库对下游洪水具有更强的影响。这些发现共同强调了在大河网水库梯级评估中纳入网络属性的关键意义：如果仅关注水库规模这一单一因素，则可能过度简化复杂系统，并低估水库群的联合洪水削减能力。
4 讨论与启示
4.1 分级衔接与互补效应
本研究结果对“大型水库主导洪水削减”的传统认知提出了挑战。尽管小型水库蓄水容量有限，但其数量和分布广度可累积形成显著的蓄滞效应。在100年一遇暴雨条件下，超过一半的水库（无论大小）将洪峰削减了50%以上。大型水库在系统中扮演着“闸门”角色，而数量众多的小型水库则有助于削减上游洪水，防止大型下游水库被过度淹没，体现了大小水库梯级间的互补作用。研究人员提出的逐库“移除”情景方法具有创新性，可广泛应用于坝体拆除决策等更广泛的系统分析领域。
4.2 淹没制图的新视角
当前水库效应量化研究主要基于流量指标的统计分析，而以淹没范围为风险指标的评估相对有限。尽管研究结果表明洪峰削减与淹没削减之间存在强相关性，但仅关注单一流量指标可能高估或低估水库的真实调节功能。流量指标反映了洪水规模，却无法捕捉洪水的空间分布格局。在针对内切斯-萨宾流域包含14座水库的深度变化实验中，小型水库削减洪量1.40 × 10? m3，大型水库削减1.08 × 10? m3，小型水库贡献了约13%的全流域削减量，该效应虽主要集中在局部上游区域，但具有可量化且水文意义上的显著性。流量与淹没双重指标的整合框架为水库管理、老化坝体修复与拆除评估以及流域韧性评估提供了更全面的信息基础。
4.3 洪水强度敏感性
在内切斯-萨宾流域开展的多重现期（10年、50年与100年一遇）模拟实验中，小型水库对洪水强度表现出显著的敏感性：10年一遇低强度洪水下其洪峰削减率不足1%；50年一遇中等洪水时增至约16%，表明蓄水空间得到了有效利用；但100年一遇极端洪水时降至约12%，此时蓄水空间更易被超越。这表明小型水库的洪峰削减并非与洪水强度成线性比例关系。相比之下，大型水库在所有洪水强度下均表现出相对稳定的削减效果。
4.4 局限性与展望
研究人员也坦诚指出模型存在的不确定性来源：包括半干旱地区低流量模拟偏差、校准-验证期相对较短、全流域统一假设4天100年一遇降雨事件的空间均匀性、出流率采用日均值而对实际调度规则的表征有限，以及90 m分辨率DEM对HAND淹没制图精度的影响等。值得指出的是，经面积折减系数（ARF）修正降水的敏感性试验表明，修正后面平均的洪峰削减量约降低10%，中位数减少6.3%，说明研究结论具有稳健性。研究进一步指出，由于不同水文气候区（如南亚季风区或美国西部融雪主导流域）的水库相互作用机制可能存在差异，本研究的结论向全球其他水库梯级系统推广时需持谨慎态度。
5 核心结论与研究意义
本研究的核心结论可归纳为以下三个层面：（1）小型水库在下游洪水削减中确实发挥着实质性作用，仅聚焦于少数大型水库的研究框架可能遗漏关键的流域系统动态效应。（2）小型水库的洪水削减效应在上游区域更为明显，但这一现象与研究区域的整体水资源可用性密切相关。（3）尽管水库规模具有影响，但网络相关因素——包括集水面积、上游水库数量、上游水库累积最大储水容量以及河序（即位置）——往往同等甚至更为关键地决定了水库的洪水削减潜力。
这项研究发表在《Earth's Future》期刊，其科学贡献不仅在于提供了小型水库在大流域中发挥显著洪水削减功能的量化证据，更在于方法论层面的创新——将洪峰流量与淹没制图双重指标相结合，系统区分了水库特征与网络属性的相对作用。研究结果为更具适应性和可持续性的水库管理策略制定提供了科学依据，并对水坝老化风险评估、河流连通性维持以及新水库选址优化等新兴挑战具有重要的参考价值。未来研究可将洪水风险概念更紧密地嵌入水库梯级建模，纳入下游人口暴露和基础设施影响评估，以推动流域洪水风险管理从传统流量调控范式向综合性风险减缓范式的转变。
