定线方法 (Rating Curve Methods)
定线方法的选择完全由水位流量关系的类型决定:关系是什么样,方法就是什么样。实测点据的分布形态分两类——单值关系(一个水位对应唯一流量,点据密集成线)和多值关系(同一水位对应多个流量,点据散开成绳套或一簇簇)。多值关系的成因主要有洪水涨落影响、回水顶托、断面冲淤变化。
从曼宁公式看定线方法的来源
明渠流量的基本公式(曼宁公式):
\begin{equation} Q = (\frac{1}{n})·A(H)·R(H)^{\frac{2}{3}}·S^{\frac{1}{2}} \end{equation}
- $n$:糙率系数(河床粗糙程度,近似常数)
- $A$:过水断面面积——水位的函数 A(H)
- $R$:水力半径(A/湿周)——水位的函数 R(H)
- $S$:能坡(水面比降)
由此可见, 在稳定流条件下,假定坡降 $S$ 和粗糙系数 $n$ 基本不变,流量$Q$主要取决于$A$ 和 $R^{\frac{2}{3}}$ 。
幂律方程
对于宽浅的天然河流,水深 $h$ 远小于河宽 $B$,此时湿周 $P \approx B$,水力半径 $R = \frac{A}{P} \approx \frac{B \cdot h}{B} = h$。过水面积 $A$ 也可以近似表示为水深 $h$ 的函数:$A \propto h^m$。将这些几何关系代回曼宁公式: \begin{equation} Q \propto h^m \cdot (h)^{2/3} = h^{m + 2/3} \end{equation} 如果考虑到河床最低点的高程(即零流量水位 $h_0$,当水位 $h \le h_0$ 时,流量 $Q = 0$),将实际有效水深写为 $(h - h_0)$,最终即可导出如下形式的幂律方程(Power Law Equation): \begin{equation} Q = C (h - h_0)^b \end{equation} 其中:
- $h$:实测水文站水位。
- $h_0$:零流量水位(Zero-flow Stage),代表水流停止时的底高程,是强物理约束参数。
- $C$:比例系数(包含了河道坡降 $S$、阻力 $n$ 及断面宽度的综合效应)。
- $b$:幂指数(通常在 $1.5 \sim 2.5$ 之间)。对于标准矩形断面,$b \approx 1.67$(即 $5/3$);对于三角形或抛物线断面,$b$ 会更大。
分段幂律方程
如果天然河道是一个标准的规则几何体(如平整的矩形混凝土渠道),用单一的幂律方程 $Q = C(h - h_0)^b$ 就能完美拟合全水位的流量。但天然河道是非规则的,最典型的是复式断面(Compound Channel):
简图如下:
当水位发生变化时,过水断面的几何形态和水力特性会发生非连续突变.
在低水期(水位仅在主槽内):水流完全限制在主槽中。此时主槽河床较光滑($n$ 较小),断面较窄,对应的几何指数为 $b_1$,比例系数为 $C_1$。 \begin{equation} Q = C_1 (h - h_0)^{b_1} \quad (h_0 < h \le h_1) \end{equation}
在高水期(水位漫过主槽,漫上滩地):一旦水位超过突破点 $h_1$(主槽台口水位/转折点 Breakpoint),水流开始漫上滩地:过水面积 $A$ 剧增:水深稍微增加一点,过水面积就会成倍扩大。阻力 $n$ 突变:滩地上通常长满了植被、树木或有建筑物,粗糙度远大于主槽。水力半径 $R$ 发生突变:湿周 $P$ 突然大幅增加。此时,旧的参数 $C_1, b_1$ 完全失效。必须以 $h_1$ 为界,启用第二段甚至第三段参数($C_2, b_2$),以适应滩地上的水力学规律: \begin{equation} Q = C_2 (h - h_0)^{b_2} \quad (h > h_1) \end{equation}
所谓“分段幂律方程”,就是用分段函数的数学形式,去表达不同水深阶段下由曼宁公式决定的水力学规律: \begin{equation} Q(h) = \begin{cases} C_1 (h - h_0)^{b_1}, & h_0 < h \le h_1 \quad \text{(主槽段/低水)} \\ C_2 (h - h_0)^{b_2}, & h_1 < h \le h_2 \quad \text{(过渡段/中水)} \\ C_3 (h - h_0)^{b_3}, & h > h_2 \quad \text{(漫滩段/特大洪水)} \end{cases} \end{equation} 一些开源软件开源软件,本质上就是在寻找最佳的转折点 $h_1, h_2$ 和零流量水位 $h_0$,使得每个区段内的拟合结果既严格符合曼宁公式的物理约束,又能顺畅地连接成一条光滑且连续的水位流量关系曲线。
所有定线方法,都是对“S 听不听话、A 稳不稳定”两个问题不同假设下的曼宁公式退化。
- 单一曲线法 = 假设 S 恒定、断面稳定:曼宁公式退化成单值函数 $Q = f(H)$。幂函数形式 $Q = a(H − H_0)^b$ 中的指数 b 有真实水力含义:宽浅矩形断面 A ∝ h、R ∝ h,Q ∝ h^(5/3),b ≈ 1.67;三角形(V 形)断面 A ∝ h²,b ≈ 2.67。b 的范围(约 1.5~2.7)对应断面形状——定线时看拟合出的 b 值就能反推断面形态是否合理。
- 落差指数法 = 把变化的 S 显式写进公式:涨水落水、回水顶托时 S 不再恒定——绳套的物理根源就在 S。用实测落差 ΔH 代表 S(S ≈ ΔH/L),代入曼宁公式整理得 $Q = a(H − H_0)^b·(ΔH/ΔH_c)^c$,其中 c 的理论值就是 1/2——因为 S 在曼宁公式里开平方根。实际定线时 c 在 0.5 附近试算修正,起点不是拍脑袋定的。
- 抵偿河长法 = 用选断面吸收 S 的变化:不修正 S,而是选一段“抵偿河长”,使中断面水位本身能抵偿上下游落差的变化,Q 与中断面水位的关系重新变回单值。
- 改正水位法 / 临时曲线法 = A 变了(断面冲淤):冲淤改变 A(和 R),同一水位下流量自然不同。改正水位法把水位“改正”到同一断面面积的当量水位,本质是把被冲淤改写的 A(H) 修正回来;临时曲线法则按冲淤阶段各定一条曲线,分段逼近变化的 A(H)。
- 高水延长 = 直接搬曼宁公式:实测断面几何(A、R 随 H 延伸),按曼宁公式把曲线外推到高水位——规范里称“曼宁公式法延长”。
单一曲线法
适用于单值关系:河床稳定、不受回水/冲淤影响,点据密集成线。
- 图解法:点绘实测点,目估或透明尺适线,三点法检验——传统手艺
- 分析法:回归拟合(幂函数 $Q = a(H − H_o)^b$ 取对数成直线,或多项式),计算机定线的常规做法
特点:简单、精度高,是定线方法的“主航道”。
水力因素法
适用于复式关系(绳套曲线)。核心思想:流量不只是水位的函数,还受水力因素(主要是落差)影响:$Q = f(H, \Delta H)$。把落差作为校正参数,把绳套“拉直”成单值关系。
- 落差指数法:用上/下游水位差 ΔH 校正,$Q = a(H − H_o)^b·(\Delta H/\Delta H_c)^c$,先试算落差再迭代定线
- 抵偿河长法:在上游选一段“抵偿河长”,用中断面水位定线——选得好时,中断面水位自动抵消涨落影响
特点:能处理绳套,但需要落差观测资料、计算较复杂。
特殊情况处理
断面变形、变动回水等情形下,用补充手段:
- 临时曲线法:断面冲淤使关系改变时,按冲淤时段分别定线(每段一条临时曲线),用测点时序判断当前用哪条——冲淤频繁的沙质河床常用
- 改正水位法:把冲淤影响折算成水位改正数,将实际水位改正到同一断面条件的当量水位,再按单值定线
- 连时序法:受变动回水影响、点据散乱时,按实测点时间先后顺序连成折线,推流时按时间查线——最“笨”但最稳
方法选择
点据是否密集单值?
├─ 是 → 单一曲线法(图解 / 分析)
└─ 否 → 看散开的原因:
├─ 洪水涨落、回水 → 水力因素法(落差指数 / 抵偿河长)
├─ 断面冲淤 → 临时曲线法 / 改正水位法
└─ 变动回水严重、点据乱 → 连时序法
检验与延长
无论用哪种方法定线,之后都要做两件事:
- 检验:符号检验、适线检验、偏离数值检验(GB 50179 有明确精度指标),不合格重定
- 延长:实测覆盖不到的高水位(曼宁公式外推)和低水位(延伸到断流水位 H₀)
软件与库
针对定线需求,一些组织推出了库和软件
- ratingcurve是美国地质调查局 USGS 开源定线库,采用优化算法拟合分段幂函数(Segmented Power Law Model),解决了多段定线中断点(Control Points)和过零点 $h_0$ 极难收敛的优化难题
- HydRA (nr-hydrometric-rating-curve)是加拿大不列颠哥伦比亚省(BC 省)政府开发的一个水文评级曲线应用程序, 旨在帮助用户基于实测的水位(stage) 和流量(discharge) 数据,自动或手动拟合方程,并生成可用于将连续水位数据转换为流量的表或方程
- bdrc是一个用 R 语言开发的包,使用贝叶斯分层模型来拟合水位-流量关系曲线。它最大的亮点是实现了 Hrafnkelsson 等人(2022)提出的“广义幂律模型”。传统幂律模型的形式是 $Q(h) = a(h-c)^b$,其中指数 b 是一个固定常数。而广义幂律模型将其改进为 $Q(h) = a(h-c)^{f(h)}$,这里的指数 f(h) 变成了水位 h 的函数。这意味着模型能够根据河道断面的几何形状,自动调整幂律指数,从而更灵活、更准确地描述复杂的水位-流量关系