基于滞回效应的非恒定渐变流条件下径流监测方法评估与优选框架

《Water Resources Research》：A Hysteresis-Based Framework for Evaluating and Selecting Streamflow Monitoring Approaches Under Unsteady Gradually Varied Flow Conditions

　　

摘要专业翻译 由水位-流量关系曲线（HQRC）估算的径流无法充分表征非恒定渐变流（UGVF）的时空复杂性，且通常在不附带不确定度的情况下对外发布，致使其被当作确定性量值对待。随着数据用户对跨机构比较及科学应用场景中径流估算的可信度提出更高要求，上述局限性日益突出。HQRC误差中一个占主导地位但长期缺乏定量刻画的主要来源，是由UGVF所特有的普遍非恒定与非均匀流条件引发的滞回效应（hysteresis）。水文测量机构试图通过替代方案来缓解这些影响，例如采用HQRC修正方法，或引入以流速或自由水面比降（free-surface slope）补充水位的独立多变量方法。然而，评估这些替代方案的有效性仍具挑战性：HQRC修正方法依赖于简化的洪水波表征，而多变量方法则保留了传统HQRC公式的半经验基础。为弥合这一评估缺口，本研究提出了一种以滞回效应为核心的性能度量指标，用于量化HQRC修正方法与多变量方法所获得的二元关系相对于常规HQRC定线之间的偏离程度。该指标整合了三个互补性指示量：（a）流量偏差幅度（flow deviation magnitude）；（b）水文过程线洪峰相位（hydrograph peak phasing）；以及（c）累积水量差异（cumulative volumetric differentiation），共同刻画滞回严重程度、洪峰时序变化及滞回对径流估算的累积影响。该框架被应用于美国地质调查局（USGS）三个涵盖不同水文地貌条件的流量站。研究结果为在主导洪水波行为（运动波型vs.非运动波型）背景下评估与优选替代监测方法提供了实践指导，并综合形成了一套面向UGVF条件下监测系统优选的结构化决策工作流程。

论文解读

一、研究背景与问题

径流记录是全球应用最广泛的环境数据集之一，支撑水资源评估、洪水预报、生态管理、水利工程及水文模型开发等众多领域。然而，全球超过95%的流量站仍依赖一个多世纪前发展起来的水位-流量关系曲线（HQRC）方法进行径流监测。该方法基于恒定均匀流假设，通过单一水位观测结合半经验定线推求流量。尽管HQRC具有简单、低成本和长历史记录等优势，其隐含的运动波（kinematic wave）假设在天然河流条件下经常被违反。当水流行为偏离这一理想化状态时，基于HQRC的径流记录可能包含显著误差。

非恒定渐变流（UGVF）是天然河流中普遍存在的水力条件，涵盖非恒定浅水流动及本质上非均匀的洪水波。在UGVF条件下，水力变量呈现路径依赖性，产生滞回关系（hysteretic relationships），即相同水位在洪水涨水段与退水段对应不同流量。这种滞回现象并非特例，而是天然河流的普遍特征，尤其在低比降系统中，回水、蓄水效应和波扩散对水力控制作用显著。然而，现有研究大多聚焦于评估单一监测方法，缺乏一个共同框架来系统评估替代方法再现UGVF水力特征的能力。滞回通常仅被视为定线误差来源，而非可测量的洪水波动力学表征，导致从业者缺乏判断常规HQRC监测何时足够、何时需要更先进方法的客观标准。

二、研究内容与结论

本研究旨在填补这一评估缺口，提出一个以滞回效应为核心的监测方法评估框架。该框架不再仅通过流量量值比较监测方法，而是评估各方法再现洪水波动力学基本特征的能力。研究提出三个互补性指示量：流量偏差幅度（flow deviation magnitude）、水文过程线洪峰相位（hydrograph peak phasing）和累积水量差异（cumulative volumetric differentiation），共同构成统一评估基础。研究在美国地质调查局（USGS）三个不同规模的流量站开展应用验证，包括佐治亚州亚特兰大市查塔胡奇河站（USGS 02336000，最小）、西弗吉尼亚州查尔斯顿市卡诺瓦河站（USGS 03198000，中等）和俄亥俄州艾恩顿市俄亥俄河站（USGS 03216070，最大）。三个站点分别代表湿润亚热带和湿润大陆性气候条件，河床比降从0.0017递减至0.00001，覆盖了从运动波主导到强滞回效应的不同水力条件。

研究得出以下主要结论：HQRC估算的径流在洪水驱动的UGVF条件下经常超过通常接受的±5%不确定度阈值；HQRC的使用会系统性地引入流量估算偏差，影响洪峰流量、洪水时间、事件水量及污染物负荷等水文过程线衍生量；纳入额外实测变量的监测方法普遍优于仅依赖水位的单一变量方法。研究进一步将滞回分析整合为监测策略优选的决策工作流程，帮助水文测量机构判断现有HQRC定线是否足够，或是否需要增强监测配置。

三、关键技术方法

研究人员采用了三类主要技术方法。第一类是基于圣维南方程（Saint-Venant equations）的物理框架，将监测方法锚定在浅水明渠流动的完整动量方程中，区分运动波、扩散波（diffusive wave）和动力波（dynamic wave）等不同洪水波类型。第二类是HQRC修正方法，包括Jones（1915）的运动波修正、Boyer（1939）的扩散波修正、Rátky（2000）和Fenton（2001）的准稳态动力波修正，以及Fread（1975）的全动力波修正，后者通过数值迭代求解并纳入详细断面几何，被用作参考流量（Q_FRE）。第三类是多变量监测方法，包括指数流速定线法（IVRC）和连续比降面积法（CSA），前者需要连续指数流速测量，后者利用多个水位传感器的自由水面比降（FSS）估算流量。所有分析基于USGS公开的15分钟间隔实测数据，数据来源于上述三个流量站及邻近辅助水位站。

四、研究结果

3.1 研究站点与仪器布局

研究人员选取三个不同规模的USGS流量站进行验证。亚特兰大站位于湿润亚热带气候区，为比降站（slope-gaging station）；查尔斯顿站和艾恩顿站位于湿润大陆性气候区，为独立IVRC站。三个站点的河床比降依次递减（0.0017、0.00005、0.00001），河宽从50米增至440米，流量范围从35–1,158 m3/s扩展至350–13,200 m3/s，覆盖了从中小河流到大河的广泛水力条件。IVRC站点的FSS通过邻近流量站的水位数据估算获得。

3.2 实验证据

亚特兰大站：分析了2009–2020年间三场持续3–5天的洪水事件。水位与FSS峰值之间存在随事件量级增大而增加的相位滞后，导致更大的滞回环。WSE-FSS关系对UGVF效应尤为敏感，其环向方向与水位-流量关系相反。五种HQRC修正方法产生相似的动态流量模式，尤其在退水段高度一致，涨水段仅有微小差异。独立比降法（CSA）的流量与Fread修正流量吻合良好，即使在漫滩条件下滞回仍然明显。

查尔斯顿站：分析了2020–2024年间三场持续6–8天的洪水事件。与亚特兰大相比，查尔斯顿表现出较小的相位滞后和滞回环，表明即使在大事件中滞回效应也较弱，其强度受河道比降和暴雨急变性等局地因素控制。FSS与流速峰值的相位滞后序列与亚特兰大相似但不显著。对于最大事件，所有HQRC修正方法均成功再现了流量动态，但IVRC和CSA估算与Fread参考流量存在明显差异。

艾恩顿站：分析了2017–2025年间三场持续5–6天的洪水事件。作为三个站点中河床比降最低者，艾恩顿表现出最强的滞回行为。随着事件量级增大，FSS和流速峰值相对于水位峰值的滞后增大，各水力变量间的滞回环也相应扩大。五种HQRC修正方法均令人满意地再现了Fread参考流量。与独立多变量方法（IVRC和CSA）的比较显示，涨水段流量略有低估，但退水期吻合良好。

五、讨论与结论

讨论部分的核心贡献在于将滞回分析整合为监测策略优选的结构化决策框架。该框架首先基于历史观测的滞回特征评估站点的主导洪水波类型（运动波、扩散波或动力波），然后通过成本效益分析权衡精度提升与仪器、数据处理和运行成本之间的关系。弱滞回且以运动波行为主导的站点可继续依赖常规HQRC定线；而持续呈现扩散波或动力波行为的站点则应考虑采用能够表征UGVF过程的增强监测方案。

研究结论指出，基于HQRC的估算在UGVF条件下的局限性凸显了建立可扩展不确定度量化框架的必要性。滞回不仅是非恒定流动力学的副产品，更是有价值的诊断与评估工具。通过流量偏差、洪峰相位和水量指标量化滞回行为，所提出的方法能够客观评估径流监测性能，为非均匀流条件下优选适配监测策略提供物理基础。研究同时指出，当前结论基于有限站点、事件和监测方法，需要在更广泛的水文气候条件、河道形态和滞回强度范围内开展进一步验证，以建立通用性能阈值并提高可迁移性。未来研究应扩展滞回方法至更多监测技术，建立与滞回严重程度关联的标准化不确定度指标，并发展可迁移的分类方案，推动从技术驱动监测向水力信息驱动的监测设计转变。