基于水量平衡原理的支流调蓄作用分析.docxVIP

基于水量平衡原理的支流调蓄作用分析 在洪水的过程中，由于支流的存在，当主流的洪水达到上升阶段时，主流的洪水将流入支流，阻塞支流的水位，即主流的洪水将淹没支流。另一个。 方面,当干流洪水处于落水阶段时,支流中被壅高的水体又将回灌干流,称为支流回灌。支流倒灌和回灌都引起干、支流水力要素如水位、流量等的变化。以前的洪水演进模型,就支流对干流洪水的影响处理只是仅仅把支流的流量叠加到干流同时刻的流量上,即忽略了支流对干流洪水的调蓄作用,这对于一般洪水问题作这样的近似处理是可以的,但对于溃坝产生的特大洪水,这种考虑显然过分粗糙。 当支流汇入干流时,一般情况下交汇角α90°。实验表明:当α70°时,支流的倒灌及回灌受交汇角的影响不大,而主要受交汇处支流水面宽及支流底坡坡度影响,综合体现为对应于某一水深支流所形成的库容,这为建立支流倒回灌模型提供了科学依据。 1 全球思想及其数学模型的研究 1.1 干流回灌 当干流洪水处于起涨阶段时,干流洪水将会向支流倒灌,因此支流在Δt时段内库容增加,水位上涨;当干流洪水落水阶段时,引起支流回灌干流,因此支流库容减少,水位下降。所以,支流的倒回灌可以设想为一个一定库容的水库的蓄泄问题,并由此建立支流水库的蓄泄关系。 依据干支流交汇处水流的连续性建立干支流联结条件,联解支流水库的蓄泄关系和干流非恒定流方程组,运用迭代法即可同时求解支流倒回灌问题和干流洪水波向下游演进问题。 1.2 非恒流的基本特征 设支流在断面b-b处与干流交汇,如图1所示,将支流视作一水库,将支流的汇入视作干流的集中旁侧入流,其流量称为干支流交换流量,记为qe(t),并约定干流注入支流时为正,反之为负;q0(t)为支流上游原有来流量。此两流量的汇入将使支流的库容和水位发生变化,依据水量平衡,可得如下支流水库蓄泄关系。 ΔWb=∫tt0tt0[q0(t)+qe(t)]dt=∫hbh0hbh0F(h)dh(1) 式中,ΔWb为Δt时段内支流库容变化;hb和h0分别对应于t和t0时刻干支流交汇处的水位;F(h)为对应于水位h时的支流库水面面积。 (1)式与梯形公式近似,得 另一方面,交换流量qe的存在又要影响干流非恒定流向下游的演进。 干流洪水波向下游的演进由圣维南方程组控制,具体形式为: ?Q?x+?(A+A0)?t=0(3)?Q?t+?(Q2/A)?x+gA(?h?x+Sf+Se)=0(4)?Q?x+?(A+A0)?t=0(3)?Q?t+?(Q2/A)?x+gA(?h?x+Sf+Se)=0(4) 式中,Q为流量;A为干流有效过水断面面积;A0为非河槽蓄水的断面面积(滩地面积);x为顺水流方向的距离;t为时间;Sf为摩阻比降,Sf=(n2Q2)/(A2R4/3);Se为局部损失,Se=kΔ(Q/A)2/(2gΔx);g为重力加速度;h为水位;k为干流河道断面收缩扩张系数,收缩时取正值,否则为0;n为满宁糙率系数;R为水力半径。 对于洪水波向下游演进的数学模型,有关资料很多。在此采用隐式加权四点差分格式进行求解,其差分方程如下: θ[Qj+1i+1-Qj+1iΔxi]+(1-θ)[Qji+1-QjiΔxi]+1Δtj[(A+A0)j+1i+1+(A+A0)j+1i-(A+A0)ji+1-(A+A0)ji]=0(5)(Qj+1i+1+Qj+1i-Qji+1-Qji2Δtj)+θ[(Q2/A)j+1i+1-(Q2/A)j+1iΔxi+gˉAj+1)?(hj+1i+1-hj+1iΔxi+ˉSj+1f+Sj+1e)]+(1-θ)?[(Q2/A)ji+1-(Q2/A)jiΔxi+gˉAj(hji+1-hjiΔxi+ˉSjf+Sje)]=0(6)θ[Qj+1i+1?Qj+1iΔxi]+(1?θ)[Qji+1?QjiΔxi]+1Δtj[(A+A0)j+1i+1+(A+A0)j+1i?(A+A0)ji+1?(A+A0)ji]=0(5)(Qj+1i+1+Qj+1i?Qji+1?Qji2Δtj)+θ[(Q2/A)j+1i+1?(Q2/A)j+1iΔxi+gAˉˉˉj+1)?(hj+1i+1?hj+1iΔxi+Sˉˉj+1f+Sj+1e)]+(1?θ)?[(Q2/A)ji+1?(Q2/A)jiΔxi+gAˉˉˉj(hji+1?hjiΔxi+Sˉˉjf+Sje)]=0(6) 其中: ˉA=(Ai+Ai+1)/2ˉSf=n2ˉQ2/(ˉA2ˉR3/4)ˉQ=(ˉQi+ˉQi+1)/2ˉR=(ˉA/ˉB)ˉB=(Bi+Bi+1)/2。 式中,B为干流河道过水断面水面宽;i为沿河道或分的空间网格;j为时间网格;θ为加权因子,0θ1.0。 在用差分方程求解非恒定流时,必须首先知道初始条件和边界条件。初始条件为初始时刻(t=0)各断面的水位和流量(h和Q),本模型假定初始时刻之