基于水量平衡方程的水庫調(diào)洪演算基礎(chǔ)數(shù)據(jù)誤差傳播規(guī)律研究

周 斌

(汕尾市水利水電規(guī)劃設(shè)計(jì)院，廣東 汕尾 516600)

為確定水庫洪水位，需要根據(jù)水量平衡原理進(jìn)行水庫調(diào)洪演算。由于方程各項(xiàng)積分難以直接計(jì)算，需要時(shí)間分段、逐時(shí)段連續(xù)求解，工程中常采用梯形法進(jìn)行單步時(shí)段積分。有諸多學(xué)者嘗試了Runge-Kutta法、Lagrange插值[1]、三次樣條插值函數(shù)[2]等諸多數(shù)值方法以提高單步時(shí)段積分的精度。筆者曾對梯形法的積分誤差和傳播進(jìn)行了討論，提出了最大可能誤差的計(jì)算方法[3]以控制計(jì)算誤差風(fēng)險(xiǎn)。入庫洪水、泄洪能力和庫容的基礎(chǔ)數(shù)據(jù)如存在誤差，也將直接影響到計(jì)算的調(diào)洪成果精度，曾有學(xué)者對以上三者隨機(jī)過程產(chǎn)生的蓄洪量的特性進(jìn)行了相關(guān)研究[4]。近年來隨著洪水預(yù)報(bào)技術(shù)的發(fā)展，不少學(xué)者也開展了基于水文預(yù)報(bào)的調(diào)洪演算研究[5]。因此，分析基礎(chǔ)數(shù)據(jù)誤差對調(diào)洪成果精度的影響是必要的。

1 調(diào)洪演算的水量平衡基本方程

對任一時(shí)段(如t1→t2)，入庫水量減去出庫水量，等于該時(shí)段內(nèi)水庫增加或減少的蓄水量[5]，此時(shí)水庫水位由z1升降至z2，需滿足水量平衡方程：

(1)

式中Q(t)——入庫流量過程；q(z(t))——出庫流量過程；F(z)——水庫庫面面積；z(t)——庫水位隨時(shí)間的變化過程。

采用梯形數(shù)值積分法[6]，式(1)可改寫為常見的實(shí)用公式(以下稱實(shí)用水量平衡方程)：

(2)

式中Qt1、Qt2——時(shí)段始、末的入庫流量；qz1、qz2——時(shí)段始、末的出庫流量；Vz1、Vz2——時(shí)段始、末的庫容；Δt——計(jì)算時(shí)段，Δt=t2-t1。

2 時(shí)段內(nèi)的單步誤差方程及單步誤差傳遞規(guī)律

2.1 單一時(shí)段內(nèi)的水位變化量

式(2)為非線性方程，已知時(shí)段初的庫水位z1可以直接解得時(shí)段末的庫水位z2。為方便后續(xù)討論，先討論一個(gè)近似解答。

(3)

2.2 入庫洪水誤差產(chǎn)生的時(shí)段末水位單步誤差方程

在計(jì)算時(shí)段(如t1→t2)內(nèi)入庫洪水流量存在相對誤差ξQ，即入庫洪水流量可表述為(1+ξQ)Q(t)。由于誤差的存在，將使時(shí)段末水位產(chǎn)生誤差εQ-z 2(相對誤差為ξQ-z2)，即時(shí)段末的水位為z2+εQ-z2。則式(2)變?yōu)椋?/p>

(4)

將式(4)采用一階泰勒公式在z2處展開，并略去高階微量，可寫為：

(5)

式(5)中消去式(2)，整理后有：

(6)

(7)

2.3 泄洪能力誤差產(chǎn)生的時(shí)段末水位單步誤差方程

如出庫流量(泄洪能力)存在相對誤差ξq，即出庫流量可表述為(1+ξq)q(z)。由于誤差的存在，將使時(shí)段末水位產(chǎn)生誤差εq-z2(相對誤差為ξq-z2)，即時(shí)段末的水位為z2+εq-z2。式(2)可寫為：

(8)

將式(8)采用一階泰勒公式在z2處展開，并略去高階微量，可寫為：

(9)

式(9)中消去式(2)，整理后有：

(10)

(11)

2.4 庫容誤差產(chǎn)生的時(shí)段末水位單步誤差方程

如庫容曲線存在相對誤差ξV，即庫容可表述為(1+ξV)V(z)。由于誤差的存在，將使時(shí)段末水位產(chǎn)生誤差εV-z2 2(相對誤差為ξV-z2)，將使得時(shí)段末水位變?yōu)閦2+εV-z2，式(2)可寫為：

(12)

將式(12)采用一階泰勒公式在z2處展開，并略去高階微量，可寫為：

(13)

式(13)中消去式(2)，整理后有：

(14)

(15)

2.5 水位誤差敏感性排序

(16)

(17)

(18)

2.6 單步誤差傳遞規(guī)律

根據(jù)前述方程，基礎(chǔ)數(shù)據(jù)的相對誤差會產(chǎn)生單步計(jì)算時(shí)段末的水位誤差，匯總得從基礎(chǔ)數(shù)據(jù)的相對誤差傳遞至?xí)r段末的庫水位單步誤差的規(guī)律見表1。

表1 單步誤差傳遞規(guī)律

接近庫水位極大值(滿足Qt1+Qt2=2qz1)處ξQ-z2和ξq-z2存在奇點(diǎn)，且ξQ-z2和ξq-z2在奇點(diǎn)兩側(cè)正負(fù)號相反。

εV-z2在庫水位極大值兩側(cè)正負(fù)號相反。

εq-z2和εV-z2在庫水位極大值點(diǎn)前有交點(diǎn)，交點(diǎn)前|εq-z2|<|εV-z2|。

3 庫水位誤差累積方程及傳播規(guī)律

3.1 時(shí)段初、末的誤差傳播基本方程

對任一時(shí)段(如t1→t2)，如時(shí)段初存在水位誤差εz1(初始水位為z1+εz1)，將使時(shí)段末的水位變?yōu)閦2+εz2，有[3]：

(19)

3.2 洪水相對誤差產(chǎn)生的水位絕對誤差累積方程

水庫調(diào)洪演算時(shí)，第n時(shí)段末的庫水位總存在相對誤差ξQ，可通過式(6)和式(19)計(jì)算累積誤差，有：

(20)

3.3 出庫流量相對誤差產(chǎn)生的水位絕對誤差累積方程

水庫調(diào)洪演算時(shí)，第n時(shí)段末的庫水位總存在相對誤差ξq，可通過式(10)和式(19)計(jì)算累積誤差，有：

(21)

3.4 庫容相對誤差產(chǎn)生的水位絕對誤差累積方程

水庫調(diào)洪演算時(shí)，第n時(shí)段末的庫水位總存在相對誤差ξV，可通過式(14)和式(19)計(jì)算累積誤差，有：

(22)

3.5 累積誤差的傳播規(guī)律

4 算例

某中型水庫位于廣東省海豐縣境內(nèi)的黃江支流龍船溪上，壩址集雨面積23 km2，主河長9.2 km，河道加權(quán)平均坡降1%，樞紐建筑物有1座主壩、5座副壩、溢洪道、2座灌溉輸水涵管、2座灌溉進(jìn)水閘。

溢洪道進(jìn)口控制段為無閘控制的開敞式寬頂堰， 堰頂高程為12 m，凈寬35 m，控制段后接坡度為1∶8的陡槽段，陡槽后接消力池。以堰頂作為起調(diào)水位，遭遇洪水時(shí),即從溢洪道自動按最大渲泄能力下泄洪水的調(diào)度方式計(jì)算遭遇百年一遇洪水的水庫調(diào)洪過程。為揭示基礎(chǔ)資料誤差對調(diào)洪成果的影響，在調(diào)洪成果的基礎(chǔ)上直接分析入庫洪水偏小10%(ξQ=-0.1)、溢洪道泄流能力偏大10%(ξq=0.1)、庫容偏大10%(ξV=0.1)產(chǎn)生的庫水位誤差情況。調(diào)洪過程和水位誤差變化過程見圖1，最高洪水位附近的誤差計(jì)算成果摘錄見表2。

從誤差圖可見，水位單步絕對誤差在庫水位上升期總體上以入庫流量的誤差影響最顯著，出庫流量的誤差影響最??；庫水位降落期則以出庫流量誤差影響最大，庫容誤差產(chǎn)生了制約庫水位下降的作用；累積水位誤差以入庫流量誤差的影響最為顯著，出庫流量誤差產(chǎn)生了持續(xù)了降低庫水位的影響；時(shí)段相對誤差的絕對值仍以入庫流量誤差影響最大，入庫流量和出庫流量的相對誤差在庫水位上升至最高點(diǎn)附近有奇點(diǎn)。入庫洪水、泄洪能力和庫容三要素的相對誤差對該水庫最高庫水位的影響程度，入庫洪水最大、泄洪能力最小。

表2 誤差計(jì)算成果摘錄

5 結(jié)語

采用水量平衡方程調(diào)洪演算時(shí)，洪水、泄洪能力和庫容等原始數(shù)據(jù)可能存在誤差，將影響到計(jì)算的調(diào)洪成果精度，分析評價(jià)這些基礎(chǔ)數(shù)據(jù)誤差對調(diào)洪成果的影響對工程安全是十分有必要的。各基礎(chǔ)數(shù)據(jù)的相對誤差對時(shí)段內(nèi)的水位單步誤差有著不同的影響，庫水位上升期以入庫洪水流量最為敏感，且對庫水位相對誤差存在放大效應(yīng)；庫水位升、降期，庫容誤差產(chǎn)生的水位單步誤差正負(fù)號相反，且相對誤差大于0時(shí)對庫水位相對誤差存在縮減效應(yīng)；最高庫水位處入庫洪水流量和泄洪能力的水位單步相對誤差存在奇點(diǎn)。對于水位傳播累積誤差，有著不同的影響，在庫水位上升期，入庫洪水流量最為敏感；庫容相對誤差大于0時(shí)對庫水位相對誤差仍保持了縮減效應(yīng)。