洪水调节课程设计：从水量平衡到水库调洪演算全流程解析

简介：一份面向水利水电工程专业学生的洪水调节课程设计参考文档，完整展示三峡大学该课程设计的任务要求与计算思路。内容涵盖设计目的、工程基本资料、洪水标准确定，以及列表试算法、半图解法推求下泄流量、库容与水位变化过程的详细流程，并配有水位库容关系、溢洪道泄流能力等关键数据。文档以某以发电为主、兼顾防洪供水的中型水利枢纽为案例，依据工程规模确定其等别为Ⅲ等，并详细说明设计洪水标准采用100年一遇、校核洪水标准采用1000年一遇的查表过程，对理解洪水标准选取和调洪演算方法有直接帮助。资源为单个Word文档，共1个doc文件，压缩包大小约816KB，便于下载后按需修改与排版。目前已有87人学习下载，适合正在完成洪水调节课程设计、需要对照计算步骤与成果整理方法的水利类本科生参考。

1. 洪水调节课程设计：一场围绕水量平衡的水库推演

洪水调节课程设计最容易让人误判的地方，是觉得它只是一张张填满数据的表格。真正动手算过一遍就会发现，整份设计反复在回答一个问题：一场百年或千年洪水入库之后，水库到底拦多少、放多少、在什么时刻放，才能保证大坝和下游都不出问题。三峡大学这门课程设计把这个问题压缩成了一个可以完整复现的工程案例，从工程等别确定、洪水标准选取，到溢洪道泄流能力公式、列表试算法、半图解法，最后落到设计洪水位和调洪库容。适合正在做课程设计的水利水电专业学生，也适合刚入行做水库调度、防洪计算的工程师作为调洪演算的对照样本。

2. 计算前置：工程等别、洪水标准与溢洪道泄流能力

洪水调节不是拿到洪水过程线就直接开算，第一步是先回答“这个工程够不够格、按什么标准算”。等别定错了，后面所有设计洪水和校核洪水参数都会跟着错。

2.1 工程等别怎么定：先分等再取标准

水利水电枢纽工程按照水库总库容、防洪、治涝、灌溉、供水、发电等指标综合分等。该工程以发电为主，装机容量 5000 kW，年发电量 1372 × 10^4 kW·h，挡水建筑物为混凝土面板坝，溢洪道堰顶高程 519 m，正常蓄水位 525 m，电站发电引用流量 10 m³/s。

课程设计里通常需要同时用“水库总库容”和“装机容量”两个指标查表，然后取其中较高的工程等别。这里按库容指标可定为中型，按装机容量指标属于小（2）型范围，规范要求取两者中的最高等级，最终工程等别定为 III 等中型。这一步看着简单，却是最容易丢分的地方：只按装机容量定等，会导致设计洪水标准偏低；只按库容定等，又可能让校核洪水标准偏严，后面泄洪建筑物尺寸就得跟着放大。

工程等别确定后，永久性水工建筑物级别也跟着定下来。III 等工程的主要建筑物为 3 级，洪水标准需要查《水利水电工程洪水标准》中“山区、丘陵区永久性水工建筑物洪水标准”一栏。

2.2 设计洪水标准：100 年与 1000 年各管什么

混凝土面板坝本质上属于土石坝类型，洪水校核标准要按土石坝一列选取。III 等 3 级建筑物，设计洪水标准为 100—50 年重现期，校核洪水标准为 1000—500 年重现期。该课程设计取设计洪水标准 100 年（频率 P=1%），校核洪水标准 1000 年（频率 P=0.1%）。

这里的“100 年一遇”不是一百年只来一次，而是每年发生的概率为 1%。设计洪水用于确定正常运用情况下泄洪建筑物要能安全通过的最大洪水，校核洪水则对应非常运用情况，校核时的允许风险更低、对坝体安全裕度要求更高。计算工况的区别直接反映在溢洪道开启孔数上：设计洪水来临时用左右 2 孔泄洪，校核洪水来临时 3 孔全部参与泄洪。

2.3 泄流能力公式：下泄流量 q 的骨架

溢洪道采用 3 孔 8 m × 6 m 弧形门控制，堰顶高程 519 m。无闸门控制时，溢洪道泄流按堰流公式计算：

q = n · ε · m · b · √(2g) · H0^1.5 + 10

各参数含义如下表：

H0 是溢洪道堰顶以上的水头。水库水位 Z 越高，堰顶水头越大，下泄流量按 1.5 次方增长。要注意的是，当 Z 低于 519 m 时，溢洪道不出流，但发电洞仍按 10 m³/s 计入下泄。这个 10 m³/s 在调洪过程中始终存在，属于“底流”，不能漏掉。

import math def spillway_q(z, n_holes=2): """ 计算溢洪道总泄流量（含发电引用流量） 参数： z : 水库水位，m n_holes: 开启孔数，设计洪水取 2，校核洪水取 3 """ g = 9.81 eps = 0.92 m_coef = 0.48 b = 8.0 h0 = z - 519.0 if h0 <= 0: return 10.0 q_spill = n_holes * eps * m_coef * b * math.sqrt(2 * g) * h0 ** 1.5 return q_spill + 10.0

这段代码直接对应上面的公式，h0 小于等于 0 时只返回发电引用流量。实际演算中，H0 会随水位变化，所以 spillway_q 在每一轮迭代里都会被重新调用，而不是一次性算死。系数 ε 和 m 在初步设计阶段可以按规范推荐值取，但若墩头形状、堰面型式与假定差异较大，需要重新查表修正。

3. 列表试算法：迭代不出来的 q 就是不该出现的 q

列表试算法是洪水调节课程设计的核心方法之一，它的思路很直白：在每一个时段内，水库的蓄水量变化等于来水量减去出水量，然后根据蓄水量反推水位，再根据水位反推下泄流量，最后把下泄流量代回水量平衡方程检查是否一致。

3.1 水量平衡方程与蓄泄方程

时段水量平衡写成：

V2 = V1 + (Q1 + Q2) / 2 × Δt - (q1 + q2) / 2 × Δt

其中 Q1、Q2 是时段初、末的入库洪水流量，q1、q2 是时段初、末的下泄流量，V1、V2 是时段初、末的库容。这里用时段平均流量近似代替瞬时过程，Δt 的选择直接影响计算精度。该课程设计取 Δt = 2 h，换算成秒是 7200 s。

把水库容积曲线 V = f(Z) 和溢洪道泄流公式 q = f(Z) 联立，消去 Z 后得到蓄泄方程 q = f(V)。水量平衡方程中，Q1、Q2、q1、V1 都是已知量，未知的只有 q2 和 V2，而它们之间又通过 q = f(V) 相互约束。直接解方程比较繁琐，所以采用试算法反复逼近。

3.2 试算流程与收敛方式

试算的标准流程是：

1. 假设一个时段末下泄流量 q2。
2. 代入水量平衡方程求 V2。
3. 由库容曲线查得 Z2。
4. 由泄流公式计算新的 q2'。
5. 比较 q2 与 q2'，若差值小于允许误差，则本时段结束；否则用 q2' 代替 q2 重新计算。

这个迭代过程本身是收敛的，课程设计中可以手算，也可以写循环程序。手算时建议把每轮试算记成一张表，方便老师检查迭代路径；写程序时则设一个最大迭代次数，防止因为初始值给得太离谱导致不收敛。

手算试算过程通常记录成这样：

实际课程设计中每个时段要试算 3 到 5 轮，整个设计洪水过程约 30 h、15 个时段，手算工作量不小。用程序迭代可以大幅减少重复劳动。

3.3 用 Python 实现列表试算法

将水位库容曲线离散成列表后，用 numpy 的插值函数可以快速完成“水位查库容、库容查水位”的双向转换。

import numpy as np # 水位和对应库容，单位：m / 万m³ z_list = np.array([520, 522, 524, 526, 528, 530, 532, 534, 536, 538, 540]) v_list = np.array([4604.99, 4750, 4900, 5060, 5230, 5410, 5600, 5810, 6030, 6300, 6670]) def v_of_z(z): return np.interp(z, z_list, v_list) * 1e4 # 万m³ -> m³ def z_of_v(v): return np.interp(v, v_list, z_list) def routing_step(v1, q1, q_in_avg, dt_s, n_holes=2): """ 单个时段列表试算 参数： v1 : 时段初库容，m³ q1 : 时段初下泄流量，m³/s q_in_avg : 时段平均入库流量，m³/s dt_s : 计算时段长度，s """ v2 = v1 + (q_in_avg - q1) * dt_s # 初始用 q1 代替 q2 for _ in range(30): z2 = z_of_v(v2) q2 = spillway_q(z2, n_holes=n_holes) v2_new = v1 + (q_in_avg - (q1 + q2) / 2) * dt_s if abs(v2_new - v2) < 1e-3: return q2, v2_new, z2 v2 = v2_new raise RuntimeError("试算未收敛，请检查初始假设或计算时段")

这里 v_of_z 返回的是国际单位 m³，z_of_v 则相反，两者都用 np.interp 做线性插值。水位库容曲线在 520—540 m 区间内接近线性，用分段直线近似误差可控。若实际工程的库容曲线在这个区间内弯曲明显，需要加密插值点，而不是简单拉直线。

时段平均入库流量 q_in_avg 由相邻两个时刻的 Q 取平均得到。这段函数没把 “闸门控制阶段” 单独建模，原因是课程设计里闸门全开前的库水位被控制在防洪限制水位，下泄流量始终等于入库流量，不会产生蓄量变化；真正需要试算的只是闸门全开、溢洪道进入自由流态之后的时段。

3.4 从试算表里读结果

整条过程线算完后，把入库洪水过程线 Q-t 和下泄流量过程线 q-t 画在一起，两条曲线的交点就是最大下泄流量 qmax 出现的时刻。这个时刻也是水库水位最高的时刻，因为在此之前入库大于出库，水库持续蓄水；在此之后出库大于入库，水位开始回落。

该工程设计洪水工况下，qmax 出现在 t=8 h 附近，对应最高水位就是设计洪水位；校核洪水工况下，qmax 出现在 t=7 h 附近，对应校核洪水位。如果你的课程设计算出来的峰值时刻偏差超过 2 h，优先检查 Δt 是否取错、初始库容是否对上了防洪限制水位。

4. 半图解法：一条辅助线省掉一半试算

列表试算法虽然直观，但每个时段都要迭代。半图解法的思路是把 q=f(V) 曲线改画成两组辅助线，让每一时段的计算变成一次简单的查表加法，不需要反复逼近。

4.1 V/Δt ± q/2 辅助线为什么能成立

把水量平衡方程两边同时除以 Δt：

V2 / Δt + q2 / 2 = V1 / Δt - q1 / 2 + (Q1 + Q2) / 2

左边是时段末库容项 V2/Δt 加上时段末出流的一半 q2/2，右边第一项是时段初的 V1/Δt 减去 q1/2，再加上时段平均入库流量。对固定的水库来说，V/Δt - q/2 和 V/Δt + q/2 都只是水位 Z 的单值函数，于是可以预先画出两条辅助线：

• f1(Z) = V/Δt - q/2
• f2(Z) = V/Δt + q/2

这样调洪计算就从“试算 q2”变成了“查线”：

1. 由时段初 Z1 查 f1(Z1)。
2. 加上时段平均入库流量，得到 f2(Z2)。
3. 由 f2(Z2) 在辅助线上反查 q2 和 Z2。
4. 继续下一时段。

每一步都是直接查表，不再需要假设和校核，这是半图解法最大的优势。

4.2 辅助线的计算与绘制

辅助线计算仍然基于泄流公式和库容曲线。以 Z=540 m、3 孔全开为例：

• V = 6670 万 m³ = 6.67 × 10^7 m³
• Δt = 2 h = 7200 s
• V/Δt = 9264 m³/s
• q = 3 孔泄流量 + 10 ≈ 4530 m³/s
• f1 = V/Δt - q/2 ≈ 6999 m³/s
• f2 = V/Δt + q/2 ≈ 11529 m³/s

辅助线表的标准形式如下：

绘图时横轴取 f1 或 f2，纵轴取水位 Z，q 曲线单独画在另一张图上。实际课程设计中由于 q 随 Z 变化较缓，辅助线可以近似成折线，计算分段线性方程即可，不必拟合高次曲线。

4.3 半图解法调洪计算与代码

半图解法每一步都只需要算一个累加值。下面的代码把辅助线构建和调洪过程合在一起：

def build_aux_curves(z_list, v_list, dt_s, n_holes=2): v_m3 = np.array(v_list) * 1e4 q_arr = np.array([spillway_q(z, n_holes) for z in z_list]) f_minus = v_m3 / dt_s - q_arr / 2 f_plus = v_m3 / dt_s + q_arr / 2 return f_minus, f_plus, q_arr def semi_graphical_routing(z1, v1, q1, q_in_avg, f_minus, f_plus, q_arr, z_list, dt_s): f1 = v1 / dt_s - q1 / 2 f2_target = f1 + q_in_avg # 由 f2 反查 q idx = np.searchsorted(f_plus, f2_target) if idx >= len(f_plus): idx = len(f_plus) - 1 q2 = np.interp(f2_target, f_plus[idx-1:idx+1], q_arr[idx-1:idx+1]) z2 = np.interp(q2, q_arr, z_list) v2 = (f2_target - q2 / 2) * dt_s return q2, v2, z2

np.searchsorted 先定位 f2 所在的区间，再用 np.interp 做局部线性插值，避免整条曲线非单调导致的反查错误。半图解法要求 f2 辅助线随水位单调递增，实际工程中泄流能力随水位上升而增大，通常满足单调条件。若某个水位区间出现异常，先检查库容曲线是否有跳变。

5. 设计洪水与校核洪水：两种工况的边界条件差异

同一个水库，设计洪水 P=1% 和校核洪水 P=0.1% 的调洪结果完全不同。差异不仅来自入库洪水过程线，还来自泄洪建筑物运行方式。

5.1 孔数、起调水位与计算时段

设计洪水工况下溢洪道只开 2 孔，校核洪水工况下 3 孔全开。孔数变化直接改变 q=f(Z) 曲线，所以两个工况要分别建立泄流曲线，不能共用一张辅助线。

起调水位都是防洪限制水位，起调库容 4604.99 万 m³。汛期洪水来临前，水库先通过闸门控制使下泄流量等于入库流量，维持水位不变；当入库流量继续增大、闸门全开后，水库才开始真正蓄水。因此调洪计算第一时段初的下泄流量 q1 就是当时的入库流量，不一定是溢洪道按自由流算出来的值。

计算时段均取 Δt=2 h，校核洪水过程线比设计洪水更陡，峰值更大，但时段步长保持一致才能比较两种方法的结果。

5.2 结果怎么判定：qmax、Zmax、Vmax

调洪结果的三个关键参数：

• qmax：最大下泄流量，出现在 Q-t 与 q-t 曲线交点；
• Zmax：最高水位，与 qmax 同时发生；
• Vmax：最大库容，对应 Zmax，调洪库容 = Vmax - V起调。

按课程设计文档中的数据，半图解法设计洪水工况得到最大库容 V=5386.579 万 m³，校核洪水工况对应库容约 6090.004 万 m³。这两个库容值可以直接用来反推最高水位，也可以作为两种算法互相校核的基准。列表试算法和半图解法理论上结果应基本一致，如果误差超过 1% 到 2%，说明辅助线绘制或插值环节出了问题。

5.3 常见误用与排错

排错时先把第一时段手算一遍，对比水量平衡方程左右两侧。水位 520 m 附近溢洪道不出流，q 只有 10 m³/s，最容易验证单位是否统一。

6. 把调洪成果落成工程结论：洪水位、调洪库容与曲线复核

课程设计做到最后，不能只交几张过程表，要把调洪成果转化成工程结论：设计洪水位多高、校核洪水位多高、调洪库容多大、最大下泄流量能否被下游安全接纳。

6.1 从演算过程提炼设计参数

设计洪水工况对应的是“水库正常运用时能承受的最高水位”，校核洪水工况对应“非常运用情况下水库的极限水位”。这两个水位决定了大坝坝顶高程、防浪墙高度和溢洪道边墙超高。调洪库容则直接参与水库调度图的制定：汛限水位以下用于防洪，汛限水位以上才是动态调洪空间。

课程设计报告里建议把每个工况的 qmax、Zmax、Vmax、调洪库容列成一张汇总表，下面附 Q-t、q-t、Z-t 和 V-t 四条曲线。曲线用同一时间轴对齐，峰现时刻标出来，老师一眼就能看出你的调洪结果是否合理。

6.2 用曲线交点复核最大下泄流量

最大下泄流量不必从表格里肉眼找，可以直接用插值求两条过程线的交点：

def find_intersection(t_list, q_in_arr, q_out_arr): for i in range(1, len(t_list)): if (q_in_arr[i] - q_out_arr[i]) * (q_in_arr[i-1] - q_out_arr[i-1]) < 0: t0 = t_list[i-1] t1 = t_list[i] d0 = q_in_arr[i-1] - q_out_arr[i-1] d1 = q_in_arr[i] - q_out_arr[i] t_qmax = t0 + (t1 - t0) * (-d0) / (d1 - d0) qmax = np.interp(t_qmax, np.array(t_list), np.array(q_out_arr)) return t_qmax, qmax return None, None

这个函数用线性插值寻找 Q-t 和 q-t 的交点。两个曲线在交点处符号发生变化，说明这一点之前入库大于出库、之后出库大于入库，正好对应 qmax 和 Zmax 出现的时刻。把这个交点的数值与表格里手算得到的峰值对比，能有效发现手工试算中某一个时段是否抄错了数字。

课程设计中还有一种常见检查方法：将列表试算法和半图解法的 q-t 曲线画在同一张图里，若两条曲线重合得不好，优先检查辅助线绘制时 V/Δt 的单位，以及插值点是否足够密。两条曲线重合良好后，设计成果才算真正闭合。

本文还有配套的精品资源，点击获取