定线方法 (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 的函数。这意味着模型能够根据河道断面的几何形状，自动调整幂律指数，从而更灵活、更准确地描述复杂的水位-流量关系

目录