2.1 基本原理

2.1.1 圣维南方程组在河道流量演算中的简化应用

圣维南方程组的水力学解法以瞬态法和特征线法为代表,可以绘出两断面间的水面曲线,但计算工作十分繁琐,且要有详尽的河道地形资料。对河道流量演算而言,只需由上断面入流推算下断面出流,至于两断面间的水文情势是无关紧要的。为此,以水量平衡方程代替连续方程(见图2.1),以槽蓄方程代替动力方程,则式(1.1)及式(1.3)分别改为

式中:Qr——河段入流量,m3/s;

Qc——河段出流量,m3/s;

S——河段槽蓄量,m3

t——时间,h。

图2.1 水量平衡示意图

式(2.1)和式(2.2)便是圣维南方程组的简化形式,是流量演算法直接推求出流过程的基本方程组。式(2.2)是非线性的,一般的解法是用线性的槽蓄方程取代式(2.2)。有时为了便于求解,将式(2.1)化为有限差形式

式中:Qr,1, Qr,2——河段初、末入流量;

Qc,1, Qc,2——河段初、末出流量;

S1, S2——河段初、末槽蓄量;

Δt——时段长。

在实际的流量演算中,一次洪水的起始河槽蓄量并不为零,设为S0,如果以W表示滞蓄量,则S=S0+W,故式(2.3)改为

2.1.2 槽蓄曲线的特性

河段内槽蓄量的大小,是由河段内各处断面面积所决定的,若河段为棱柱形,则槽蓄量的大小由各处水位所决定。当水流呈稳定流状态时,槽蓄量是任一断面处水位的函数,而此时的水位与稳定流流量存在着单值关系。如用河段下断面(出口断面),则有

式中:Zx——下断面(出口断面)水位,m;

Qw——稳定流流量,m3/s。

式中:Sw——下断面水流处于稳定流流量时相应的槽蓄量,m3

一般情况下,水流是处在不稳定流状态,则上述关系不复存在,这是由于河床的冲淤变化使水面不是直线和附加比降的作用所致。在一般情况下水面近似直线,故用直线表示水面线其误差可以忽略(流量演算诸法均以水面线为直线的假定为前提)。

在不考虑河床冲淤变化时,影响水位流量关系和槽蓄曲线的只有附加比降。由图2.2可知,在涨洪段的水面比降i1大于稳定流水面比降i0;在落洪段水面比降i2小于稳定流水面比降 i0。由曼宁(Manning)公式Q=为糙率,i为水面比降,h为水深)知,在同一水深h(即同一水位Z)情况下,稳定流流量,涨洪段流量,落洪段流量,其中h=Z-Z0, Z0为河底高程。因为i1i0i2,所以QzQwQl。故水位流量关系曲线呈反时针方向绳套,见图2.3。

图2.2 河段不稳定流示意图

图2.3 不稳定流ZxQc关系曲线

同样,在同一Zx情况下,由图2.2可知,涨洪段SzSw,落洪段SlSw,则水位蓄量关系曲线也呈反时针方向绳套,见图2.4。至于槽蓄曲线,即SQc关系曲线,情况就不那么简单。若将ZxQc关系曲线放在第四象限(此时ZxQc关系曲线呈顺时针方向),将ZxS关系曲线放在第三象限,见图2.5。在同一Zx情况下,涨洪段时SSw, QcQw。落洪段时SSw, QcQw。即槽蓄量和流量均随附加比降增加而增加,亦随附加比降减少而减少。因此,SQc关系曲线与ZxS曲线方向相反且绳套—定比ZxS曲线为小,甚至会出现相反的情况,这就取决于ZxQc曲线和ZxS曲线的绳套大小,若ZxQc曲线的绳套甚小,而ZxS曲线绳套较大,我们近似地以ZxQw曲线代替ZxQc曲线,由于QcQw,即在同一Zx时流量几乎不变,由ZxS曲线可知,涨洪段SSw,故此时S落在SwQw曲线左侧的a点。而在落洪段SSw,此时S落在SwQw曲线右侧的b点。所以SQc曲线呈顺时针方向绳套,见图2.5(a)。若ZxS曲线的绳套甚小,而ZxQc曲线的绳套较大;同理,近似地以ZxSw曲线代替ZxS曲线,由于SSw,即在同一Zx时槽蓄量几乎不变,由ZxQc可知,涨洪段QcQw,故此时Qc落在SwQw曲线右侧的c点。落洪段QcQw,此时Qc落在SwQw曲线左侧的d点。因此SQc曲线呈逆时针方向绳套,见图2.5(b)。综上所述,槽蓄量不仅是流量的函数,而且还是水面比降的函数,在下游出流量已定时,上游入流量便可反映水面比降,故一般用式(2.2)表示之。槽蓄曲线的坡度,其量纲为时间,K为反映调蓄(滞蓄)能力大小的参数,称之为调蓄滞时。当K大时表明调蓄能力大,当K小时表明调蓄能力小。

图2.4 不稳定流ZxS关系曲线

图2.5 附加比降对槽蓄曲线的影响示意图

(a)ZxQc绳套小而ZxS绳套大(b)ZxQc绳套大而ZxS绳套小

通过上述分析可知,流量演算的中心问题就是如何处理槽蓄曲线。