水位趨勢法在西苕溪橫塘村站的應(yīng)用

元曉華

(浙江省水文局，浙江 杭州 310001)

?

水位趨勢法在西苕溪橫塘村站的應(yīng)用

元曉華

(浙江省水文局，浙江 杭州 310001)

指出了西苕溪流域橫塘村站水位受上游洪水及下游回水頂托的雙重影響，水位流量關(guān)系曲線呈不規(guī)則繩套曲線，水位預(yù)報(bào)難度較大。從實(shí)際情況出發(fā)，采用簡單且效果明顯的水位趨勢法預(yù)報(bào)橫塘村站水位過程，取得了較為滿意的效果，可以作為其他預(yù)報(bào)方法的重要補(bǔ)充和參考。

頂托；繩套；水位預(yù)報(bào)；水位趨勢法

1　西苕溪流域概況

西苕溪位于湖州市的西南部，杭嘉湖地區(qū)西部，上游為崇山峻嶺，孝豐至安城一帶為狹長河谷平原，安城以下地勢逐漸展開，至梅溪與長興平原相連。整個流域西南高、東北低，因上游為地形復(fù)雜的山區(qū)，坡陡流急，下游河道蜿蜓曲折，加上受東苕溪導(dǎo)流及太湖頂托影響，造成洪水期泄水不暢[1]。橫塘村站為西苕溪流域的重要控制站，其水位情況在防汛工作中有著重要的警示作用，洪峰水位的預(yù)報(bào)更是下游各控制站洪峰水位預(yù)報(bào)的依據(jù)。因此，研究橫塘村站的水位預(yù)報(bào)方法，做好水位預(yù)報(bào)工作是十分必要的(圖1)。

圖1橫塘村站流域

2　預(yù)報(bào)方法

橫塘村河道站位于安城下游約3 km處，以上為山區(qū)性河流，以下河床平緩，其水位流量關(guān)系曲線受上游洪水漲落、下游回水頂托等多種因素的影響，呈復(fù)式繩套模式：同樣的流量對應(yīng)的水位相差很大。該站流量預(yù)報(bào)工作一直能取得令人滿意的效果，但因其水位流量關(guān)系紊亂，水位預(yù)報(bào)工作難度非常大。水力學(xué)模型能解決復(fù)雜的河道水位預(yù)報(bào)問題，但這種方法對資料條件要求很高，西苕溪流域不具備使用該方法的資料條件。將復(fù)雜的問題簡單化，對于水文上難以解決的問題來說不失為一種簡單易行的方法，常能取得意想不到的效果。筆者考慮主要影響因子，采用水位趨勢法研究橫塘村站水位預(yù)報(bào)問題。 受洪水漲落影響時，流量公式為[2]：

(1)

(2)

受變動回水影響時，流量公式為[2]：

(3)

(3)式中，n為糙率，A為斷面面積，R為水力半徑，Se為能面比降，在斷面穩(wěn)定、河道順直時，可用水面比降近似代替。其中，n、A、R為水位的函數(shù)，Se用河流上、下游兩站的水位差即落差ΔZ，除以兩斷面的距離計(jì)算出來，因此，

(4)

由于Z～f(Q,A)，對于主要受洪水漲落和變動回水雙重影響的斷面來說，則有：

Z～f(Q1,Q2,A)，其中，A為水位的函數(shù)，由公式(2)、(4)，可得如下關(guān)系：

(5)

針對橫塘村站水位和其主要影響因素建立回歸方程為：

(6)

(7)

采用最小二乘法擬合方程曲線，以歷史數(shù)據(jù)點(diǎn)群與回歸曲線的離差平方和最小為目標(biāo)函數(shù)，即可求出回歸方程(6)式中的系數(shù)。

3　回歸分析及參數(shù)確定

取橫塘村站1980～2007年洪峰水位接近或超過警戒水位6.67 m的洪水作為歷史洪水，共25場，其中22場參加參數(shù)率定，3場用于檢驗(yàn)。

預(yù)見期和精度是矛盾的：一方面，預(yù)見期越短、精度越高，即越靠近當(dāng)前時刻回歸方程越能反應(yīng)當(dāng)前時刻的水位趨勢；另一方面，預(yù)見期越長，精度越差，但對預(yù)報(bào)工作有利。因此不能為了追求精度而犧牲預(yù)見期，也不能為了追求預(yù)見期而犧牲精度，兩者要權(quán)衡。另外，回歸階數(shù)代表當(dāng)前時刻水位與(預(yù)見期)個小時之前個時段的各要素相關(guān)，但并不是越大相關(guān)性越好，因?yàn)槊繄龊樗加懈髯缘奶匦?，統(tǒng)計(jì)學(xué)方法從大容量樣本中尋找規(guī)律，如果樣本一致性不高，則其尋找的規(guī)律就是“大眾化”的。例如，緩漲和陡漲的兩場洪水，回歸階數(shù)會相差很大，對于緩漲洪水來說，可能過去6、7個時段的要素都能反應(yīng)未來的水位趨勢，而對于陡漲洪水來說，可能只有過去2、3個時段的要素才能反應(yīng)未來的水位趨勢。

給出幾組預(yù)見期和回歸階數(shù)的值進(jìn)行回歸分析，目標(biāo)函數(shù)值對比如表1所示。

表1　回歸分析結(jié)果　m2

綜合考慮預(yù)見期和模擬精度，當(dāng)預(yù)見期取3 h，回歸階數(shù)n取為2 h，建模效果較好，此時，求得回歸方程如下：

(8)

4　精度分析及模型檢驗(yàn)

根據(jù)《水文情報(bào)預(yù)報(bào)規(guī)范》GB/T 22482-2008，洪峰水位預(yù)報(bào)以預(yù)見期內(nèi)實(shí)測變幅的20%為許可誤差，過程預(yù)報(bào)誤差以確定性系數(shù)作為指標(biāo)。精度評定結(jié)果如表2所示(取其中的6場洪水為例)。

表2　精度評定

注：ZO,P為實(shí)測洪峰水位，TO,P為實(shí)測峰現(xiàn)時間；ZC,P為計(jì)算洪峰水位，TC,P為計(jì)算峰現(xiàn)時間；dz為洪峰水位誤差；dt為峰現(xiàn)時間誤差；dc為確定性系數(shù)

由表2可知， 22場洪水確定性系數(shù)均較高(在0.97以上)，峰現(xiàn)時間也在允許誤差范圍之內(nèi)，按照規(guī)范，洪峰水位也只有一場洪水模擬不合格。但總體來看，該方法對洪峰的模擬不夠理想，誤差較大。模擬效果圖舉例如圖2所示。

由圖2可見，模擬計(jì)算的水位在漲、落洪段模擬效果較好，而在洪峰附近波動較大，究其原因，是該方法自身的缺陷造成的。在峰段，水位漲落趨勢有突變，首先由原來的上漲變?yōu)闈q平，之后，由原來的漲平變?yōu)橄陆?，而這種實(shí)際發(fā)生的突變在發(fā)生之前無法被該方法識別，舉例來說，當(dāng)實(shí)際洪峰發(fā)生時，該方法卻認(rèn)為水位還會像之前那樣上漲。等到洪峰過后，水位呈下降趨勢，且超過回歸階數(shù)(n)個時段之后，該方法就會熟練“記憶”水位下降的趨勢，從而取得良好的模擬效果。另外一個問題是，該法只解決了水位過程預(yù)報(bào)的問題，卻不能回答洪峰具體何時出現(xiàn)。根據(jù)洪水漲落洪規(guī)律，洪峰水位出現(xiàn)于洪峰流量之后，對歷史洪水統(tǒng)計(jì)分析，橫塘村站洪峰水位約在洪峰流量之后3 h左右出現(xiàn)。因此，結(jié)合流量預(yù)報(bào)，可大致清楚洪峰水位的出現(xiàn)時間。

以2007年以后的三場洪水檢驗(yàn)?zāi)Ｐ?，模擬精度如表3所示。

表3結(jié)果表明，模型參數(shù)通過檢驗(yàn)，方案合理，可用于實(shí)時洪水過程作業(yè)預(yù)報(bào)。洪峰水位預(yù)報(bào)因方法自身缺陷誤差較大，還需要進(jìn)一步研究和改善。

圖2　歷史洪水模擬效果

洪號ZO,P/mTO,PZC,P/mTC,Pdz/cm洪峰許可誤差/cmdt峰時許可誤差/hdc200908117.622009-8-116:007.852009-8-115:002325-130.998201208097.702012-8-92:008.002012-8-823:003088-34.20.996201310088.542013-10-83:009.102013-10-83:005562030.997合格率/%100100100

5　結(jié)語

根據(jù)西苕溪流域橫塘村站水位預(yù)報(bào)困難的現(xiàn)狀，采用簡單的水位趨勢法研究水位預(yù)報(bào)問題。該方法對洪峰水位的過渡段模擬不良，容易造成大的誤差，有待改進(jìn)，而對于整個洪水過程而言，尤其是漲、落洪段，模擬效果良好。此外，該方法只解決了過程預(yù)報(bào)的問題，卻不能預(yù)報(bào)峰現(xiàn)時間，筆者利用水力學(xué)知識和統(tǒng)計(jì)學(xué)方法，從歷史洪水找規(guī)律，分析出橫塘村站洪峰水位約在洪峰流量之后3 h左右出現(xiàn)。由于該方法能夠較好地把握整個洪水過程，可以結(jié)合流量預(yù)報(bào)方法一起使用，作為其他水位預(yù)報(bào)方法的重要補(bǔ)充。

[1]席桂平. 落差指數(shù)法在西苕溪上游站洪水預(yù)報(bào)中的應(yīng)用[J]. 浙江水利水電?？茖W(xué)校學(xué)報(bào)，2002，14(4)：11～12.

[2]林傳真，周忠遠(yuǎn). 水文測驗(yàn)與查勘[M]. 南京：河海大學(xué)出版社，1987.

Application of Water Level Trend Method in Hengtang Village Station of Xitiao River

Yuan Xiaohua

(ZhejiangProvincialHydrographicOffice,Hangzhou,Zhejiang310001,China)

This paper pointed out that since the water level of Hengtang Village station in Xitiao River was influenced both by upstream flood and downstream backwater, stage-discharge relation curve showed an irregular loop curve and stage forecast was very difficult. From this point, we used simple and effective water level trend method to forecast water level of Hengtang Village, and achieved satisfactory results. Therefore the method could be used as an important supplement and reference.

back water; loop; stage forecast; water level trend method

2016-06-25

元曉華(1981—)，女，朝鮮族，工程師碩士，主要從事水文學(xué)及水資源相關(guān)科研工作。

X83

A

1674-9944(2016)16-0188-04