断面法水库库容计算的算法细节

3．3 第一层截锥的体积计算
6
底部锥体以上的计算层是按照分层截锥来计算的， 但是从下面数 的第一层形状比较特殊。见图 6，构成第一个计算层的两个底面中有 一个（对应于上断面）只有一条直线上边，下边则是由断面的最低河 底高程附近的部分河底曲线围成的。 如果上断面的最低河底高程不是 恰好等于某个分层高程时，在图 6 中就是 z 2 不等于 z j-1 时，如果对该 层采用截锥公式计算体积， 则其与以底部锥体计算的部分之间就会有 一块儿体积需要另外考虑。
j
Vzj
3b ( A
i i 1
1
1i

A1i A2 i A2 i )
j=1,2,3……
(3)
采用分层——累加的方法有利于在计算中处理间距随高程不同 的问题，而且有助于削弱上下断面形状的相似性差异所带来的误差。 由于 2 次关系不具有可加性， 分层水平面的间隔距离——分层高 程步长的采用值对（3）式的结果是有影响的，不同的分层步长导致 不同的断间库容计算结果。但是数学上可以证明，由于参与计算的函 数都是连续或分段连续的，所以分层高程步长趋近于零时， （3 ）式中 体积和的极限是存在的。工程实际中，不求这个极限，而是对分层步
127.5531 127.5520 0.0474 0.0464
127.0058 126.9946 0.0532 0.0420
123.1561 123.1569 0.0555 0.0563
119.6931 119源自文库7028 0.0449 0.0546
118.0917 118.0896 0.0552 0.0530
断面法水库库容计算的算法细节
刘炜
（黄河水利委员会 水文局，河南 郑州 450004 ）
摘要：本文论述了断面法水库库容计算的基本算法模型及间距采用、底部锥体和 回水末端处理等细节问题。 关键词：断面法 库容计算 算法
断面法是水库库容及冲淤量测算的常规方法之一， 断面法分为加 密断面法和基本断面法。 前者是通过在水库水系各干支流上布设足够 密集的测量断面（称为加密断面） ，实现对水库库容接近于地形法测 图精度的精确测算。通过减少参与计算的断面数量，经过反复对比计 算，并依照水库河道测量的有关规范要求，从加密断面中选取出一定 数量和足够代表性的断面，固定下来进行历年的常规测量和库容计 算，就构成了基本断面法，基本断面法又称为固定断面法， “断面法” 在一般情况下所指的也就是基本断面法。 基本断面是在对比计算基础 上确定的，数量少且对于库区地形变动的代表性强。因此，基本断面 法可以在相当长的时期（基本断面代表期）内，以较低的成本和较短 的测量周期实现对水库库容的准确测算。直到水库经过多年运行， 河 床形态和冲淤规律发生了显著变迁时，基本断面需要从新确定。 在基本断面代表期内， 影响库容及冲淤量成果准确性的主要因素 有两方面，一是外业测量的质量控制，二是数据处理与计算方法。本 文就后者的若干细节问题进行讨论。 1．基本算法模型及公式 水库断面法计算通常采用截锥体概化， 即假设将上下两个断面间
上断面
河底倾角θ
Z2
底部锥体 底部锥体顶面面积A2 z2
下断 面
Z1
下断面
上断面
下断面
上断面
间距bz2 Z2 Z2
间距 bz3
间距bz2
θ
Z1
Z3 Z1
(a)
(b)
图 5 底部锥体模型
回水末端锥体和底部锥体的情况很类似， 本文仅以后者为例进行 叙述。 当底部柱体完全处在一个间距 ——高程分段区间内时，例如
300 280 260 240 220 高程(m) 200 180 160 140 120 100 800 连续函数间距 分段函数间距
900
1000
1100
1200 间距(m)
1300
1400
1500
1600
图 3 断面间距高程函数
3．2 回水末端及底部锥体 在水库支流上游的最后一个实设断面以上， 河底高程到达水库设 计水位的位置虚设一个断面，称为回水末端断面。回水末端断面纯粹
定义高程范围在 z j-1 ,z j 之间的第 j 层体积为两个锥体体积之差
V j V底锥 , z j V底锥 , z j1
1 bz 2 ( z j A2 z j z j 1 A2 z j 1 ) 3 z2
(6)
当底部锥体的高程范围跨过若干个间距——高程分段区间时， 情 况就比较复杂，以图 5(b) 为例，(z 1,z 3]对应间距值 b(z3) ，(z3,z2] 对应 间距值 b(z2 )，就是说在整个锥体内，间距值不唯一。这时有两种处 理方法，一种办法是忽略间距在(z 1,z 2]区间内的变化，只采用最高高 程 z 2 所对应的间距，则第 j 层的体积仍然用（6）式计算。 另一种办法是以间距——高程函数的跳跃点 z 3 将底部锥体分成 上下两个部分， 间距分别采用。 在实际计算中， 这并不需要专门处理 ， 因为间距——高程函数是在整米数处分段的， 而分层高程步长常采用 1 米。只需要按照计算层所属的区间采用相应的间距值来定义分层体 积即可。
3
长事先作出规定，以该规定值的计算结果为准。例如对于百米级高差 的峡谷型水库（如黄河小浪底水库）规定标准计算流程中采用的分层 高程步长为 1 米。 3．算法细节 3．1 断面间距的采用 依据等高线量算得到的断面间距构成一张表格， 其内容是若干个 高程上的相邻断面间距值。 通过构造一个间距——高程的连续插值函 数， 就可以求出任意高程上的断面间距， 这相当于图 3 中的连续曲线 。 但工程上通常采用更简单的做法： 以上下两条等高线高程上的间距平 均值作为这两个高程范围内的间距来使用。例如，当计算中需要使用 175m ～200m 米之间任一高程（>175,<=200）的间距值时，都采用高 程 175m 和 200m 的实际量算间距 1270 和 1136 的平均值 1203 。这实 际上是把间距——高程作为一个分段常数函数来处理的办法。 相当于 图 3 中的分段曲线，下文中用 b(z) 或 bz 表示这个函数。
1
的河道按概化间距拉直后，其容积立体构成一个截锥体：上下断面分 别对应该截锥的两个底面，概化间距对应截锥的高。在截锥体假设的 情况下， 计算区段内的任意河道横断面在宽和深两个方向上都被认为 是沿河长线性变化的。因而其面积在上下断面间以 2 次关系变化。
A2 b
A1
b+a
图 1 断面间容积立体的截锥体假设
上断面
下断面
Zj Z2 Zj- 1
四面体
Zj Z j- 1
Zj
截锥计算层
Zj
Z2
四面体
截锥计算层 锥体计算层
Z2
Zj-1
Zj-1
(a)
图 6 第一层截锥的体积分划
(b)
一种处理方法是按照其确定的边界，直接补充一个四面体体积， 见图 6(a) 。 第二种方法见图 6(b)， 以高程 z 2 增加一个水平分层， 把 z j-1,z j 之间认定的体积分成上下两部分，下面按照一个底部锥体层计算， 上 面按照一个截锥层计算，求和得到整层体积。 通常情况下断面面积是按照 1 米分层事先计算的， 第二种方法涉 及到一个额外的非整米处的断面面积（即 A 2z2 ）计算，相对来说略为 复杂一些。 3．4 支流河口和干流坝前
4
是为计算模型的需要而设置的，它被定义为一个点，以该点为顶点， 最后一个实断面为底面，构成一个锥体（见图 4） 。回水末端锥体描 述了我们对支流最上游部分容积体的形状认定。
设计水位
回水末端
河流上游最后一个实设断面
图 4 回水末端锥体模型
在一般的相邻断面区间，上下断面的最低河底高程不相等，在水 平分层时， 其中较高的最低河底高程所对应的水平面以下的容积体被 认为是一个锥体，称为底部锥体（见图 5）。
（2）
上式即为水库库容计算的基本公式，一般称为截锥（体体积）公 式。给定一个高程，计算出水库所有基本断面的断面面积，对相邻断 面采用（2）式计算断间库容，加起来即可实现该高程下水库库容的 粗略计算。 2．体积的分层计算 截锥体概化假设计算区段上下断面形状是概化相似的， 在河道原
2
型中，这种条件不可能准确达到。另一方面，由于实际的断面间河长 是与高程有关的， 用等高线法确定的断面间距也是一个随高程变化的 量。因而，在实际的库容计算中，通常不是直接采用（2）式，而是 通过分层——累加的方法来计算逐级的断间库容。
113.2614 113.2715 0.0501 0.0603
9
zi zi-1 zi-2
断面间距 bi
面积A2i
面积A 1i 一个截锥计算层，对应于（ 3 ）式中的一个求和项
图 2 断间库容的分层计算模型
如图 2 所示，以等间距的水平面将容积立体分成若干层，对每层 的体积应用截锥公式计算，再累计求和得到各个高程下的断间库容。 设 A1i ， A2i ，b i 分别为上下断面在第 i 层边界内的面积和概化间距， 则第 j 层层顶高程对应的断间库容为：
8
表1
第一层截 锥的处理 (1) (2) (2)-(1) 补充锥体 (3) (4) (3)-(1) (4)-(1) 计算层 补充四面体 最高间距 最高间距 不处理 不处理 底部锥体 间距采用 最高间距 相应间距
黄河小浪底水库历年总库容算法结果对比
1999 127.5056 127.5088 0.0032 2000 126.9526 126.9559 0.0033 2001 123.1006 123.1036 0.0029 2002 119.6482 119.6503 0.0021 2003 118.0366 118.0383 0.0017 2004 113.2113 113.2128 0.0015
j
Vzj
3b ( A
i i 1
1
1i
j j 1 A1 i A2 i A2 i ) bi A1i bi ( A1 i A2 i ) i 1 i 1 2
(8)
j=1,2,3……
也就是说，支流河口区段的容积立体可以按照分层柱体来认定。 同理， 干流坝前区段的容积也是通过在大坝位置虚设断面按照分 层柱体来对待，同样采用（8）式来计算。 4 实际算例 针对底部截锥的间距采用和第一截锥的体积计算这两个细节问 题， 对黄河小浪底水库从 1999 年到 2004 年的 6 个汛前测次数据采用 不同的处理方法进行了库容实例计算，对比结果见表 1。
z1 , z2 175,200 ，见图 5(a) ，间距值是唯一的，其整体体积为：
1 V底锥 bz 2 A2 z 2 3
5
任意高程 z z1 , z 2 下的体积也直接按锥体计算
1 1 bz2 V底锥 , z ctg z A2 z z A2 z 3 3 z2
如图 1 所示，断间容积立体按照截锥体假设，其体积为两个锥体 体积之差：
1 1 V ( a b) A2 aA1 3 3
（1）
根据锥体的性质，有比例关系：
a2 A 1 2 ( a b) A2
解出
a
bA1 b A1 A2 A2 A1
代入（1）式简化后得到：
1 V b( A1 A1 A2 A2 ) 3
7
支流河口
干流坝前
图 7 支流河口区段及 干流坝前区段位置示意图
对于某些支流，由于地形条件及其他因素的限制，无法恰好在河 口位置设置基本断面。于是，支流最下游的实设断面与河口之间的一 块儿容积就需要另做考虑。目前的做法是在河口位置虚设河口断面， 计算时完全借用上游实设断面的数据。因而在计算公式上有如下关 系：
V j
1 bz j ( z j A2 z j z j 1 A2 z j 1 ) 3 z2
（7）
其中 bzj 的取值，显然有： 当 z j1 , z j z1 , z 3 时 bzj =b(z3)， z j1 , z j z 3 , z2 时 bzj =b(z2) 对 V j 进行累加后得到的体积对应于图 5(b)中实线所代表的部分。