湘潭縣主汛期涓水流域洪峰與降水量關系分析

胡 堅， 黃昕怡， 何澤良， 林明麗， 敬鋼球

(1.湖南省湘潭市氣象局，湖南湘潭 411100；2.湖南省湘潭縣氣象局，湖南湘潭 411200)

?

湘潭縣主汛期涓水流域洪峰與降水量關系分析

胡 堅1， 黃昕怡2， 何澤良1， 林明麗1， 敬鋼球2

(1.湖南省湘潭市氣象局，湖南湘潭 411100；2.湖南省湘潭縣氣象局，湖南湘潭 411200)

摘要利用2011~2014年涓水流域上游10個區域氣象站的日降水量及射埠水文站洪峰水位高度與流量資料，用算術平均法求得流域日面雨量，與洪峰過程進行分析，結果表明,流域累積面雨量達≥60.0 mm時，其后1~2 d將有洪峰產生；若累積面雨量≥90.0 mm，且最大日面雨量達100.0 mm以上，其后1 d將有流量達1 000 m3/s以上洪峰產生；洪水過程中，洪峰水位和流量與洪峰出現前1 d日面雨量的關系最為顯著。利用逐步回歸和多元回歸方法建立預報方程模型，發現在站點資料齊全的情況下，采用多元回歸方法具有更好的預測準確率；若站點降水量資料不全，則采用最優因子逐步回歸模型。

關鍵詞汛期降水量；洪峰；關系;湘潭縣；涓水流域

涓水位于湖南省中部，是湘江一級支流，發源于婁底市雙峰縣昌山，流經衡山縣新橋、白果、湘潭縣花石、射埠，從河口注入湘江。射埠水文站設立于1972年1月，是涓水河下游控制站，屬省級重要水文站。涓水在湘潭縣內流程67 km，流域面積1 069 km2，占全縣總面積的42.6%，是境內面積最大的一條河流。據統計，湘潭縣近20年共發生12次大的洪水，其中9次出現在涓水流域，出現時間均在主汛期4~9月。如2011年6月10~14日，涓水流域出現強降水，涓水超警戒水位，造成兩岸農田受災808.8 hm2，房屋倒塌481 間等，經濟損失0.4億。因此掌握涓水流域主汛期洪峰與降水量的關系，能夠主動為防洪搶險工作提供科學的決策依據，有利于防汛工作的順利開展。肖立炎等[1]利用萍水河上游觀測站各時段降水量與萍鄉市城區河內漲水關系，建立了二者之間的預測模型；王歡等[2]利用逐步回歸方程，建立了信江河上游雨量站點與梅港水文站洪峰期之間的預測方法；吳璐等[3]研究了沙澧河上游關鍵區面雨量與致洪之間的關系。筆者將利用涓水流域出現洪峰的數據資料，建立遇強降水時洪峰與降水量的關系及預報模型，以期更好地為當地政府及各部門提供服務，提升防災減災能力和水平。

1資料和方法

1.1資料來源利用涓水流域荷葉站、新橋站、青山橋站、龍口站、花石水庫站、分水站、排頭站、嚴沖水庫站、錦石站及射埠站10個區域自動氣象站2011~2014年4~9月08：00~次日08：00≥25.0 mm的降水量及出現日期，射埠水文站2011~2014年出現洪峰時的水位、流量及日期。若洪水過程中站點雨量缺測太多，無法準確計算流域面雨量則以空缺方式顯示。

1.2分析方法

1.2.2逐步回歸法。利用逐步回歸方法對洪峰期射埠水文站的水位高度、最大流量與前期各站點日降水量的關系進行研究，建立洪峰水位、最大流量的預測模型。

2結果與分析

2.1洪水過程雨洪關系分析洪峰的出現離不開強降水過程[8]，從2011~2014年涓水流域出現11次洪水過程的洪峰(超過警戒水位46.5m)流量與流域站點當天及前1~2d面雨量(表1)可看出,涓水流域發生1 000m3/s以上的洪峰，當天或前1d一般有較強的降水過程，一般累積面雨量達90.0mm以上，或最大日面雨量≥100.0mm;低于1 000m3/s的洪峰中，僅一次累積面雨量<50.0mm，其余過程累積面雨量均大于60.0mm;洪峰出現時間一般滯后強降水1d，11次洪水過程中有9次過程最強降水出現在洪峰前1d;涓水流域洪峰出現與前1d降水量密切相關。綜上所述，當涓水流域出現累積面雨量達90.0mm以上，且最大日面雨量≥100.0mm時，其后1d流域將有流量達1 000m3/s以上洪峰產生；若當累積面雨量達60.0mm以上時，其后1~2d流域將有洪峰產生。因此，可將流域日面雨量≥100.0mm的降水過程視為大洪水預報臨界值，將過程累積面雨量≥60.0mm的降水過程視為涓水流域可致洪的預報臨界值。

2.2涓水流域面雨量與射埠站洪峰的關系分析11次洪水過程射埠水文站每日最大流量、洪峰水位與當天及前1~2d涓水流域逐日面雨量，分別計算相關關系，結果發現(表2)，射埠水文站洪峰和流量與涓水流域逐日面雨量的相應關系不盡相同。其中，前1d涓水流域日面雨量與洪峰及流量的相關關系均最好，相關系數分別為0.83、0.84，當天次之；前2d面雨量與洪峰及流量的相關關系最弱，且相關系數均為負值，原因可能是由于涓水流域為小流域，各站點降水匯聚到射埠水文站的過程中沿程流量不斷增大，發生坦化變形，再加上流域匯流是高阻尼系統，具有高度衰減性，導致匯流中流量衰減量過大，匯聚到射埠站的流量反而減小造成的[9]。由此表明射埠水文站洪峰預報關鍵為涓水流域前1d面雨量。

表1　涓水流域流量與面雨量統計關系　mm

表2　射埠站洪水洪峰、流量與涓水流域日面雨量的相關關系

2.3洪峰水位及流量預報方程研究前期雨量對洪峰的影響，僅選用洪峰形成時期的樣本，剔除其余樣本，最后共得到11個有效樣本。涓水流域內，由于各區域氣象站到射埠水文站的距離不同，同一時間的降水對射埠站洪峰的影響程度和時間不一樣，但由于涓水流域面積較小，因此僅考慮前2 d逐日降水量對射埠站洪峰的影響。選取涓水流域10個氣象站點當天及前1~2 d的面雨量，資料樣本為11個，采用逐步回歸和多元回歸2種方法，建立射埠站洪峰過程水位以及最大日流量的3種統計預報回歸方程。設x0、x1、x2分別代表洪峰出現日、前1 d及前2 d涓水流域日面雨量，應用一元、二元回歸得出的洪峰水位高度逐步回歸方程分別為y1=45.84+0.03x1、y2=45.865+0.023x1+0.051x0，最大日流量逐步回歸方程分別為Y1=191.302+8.452x1、Y2=198.014+6.473x1+13.695x0，而運用多元回歸得到的洪峰水位高度多元回歸方程為y3=45.41+0.052x0+0.027x1+0.013x2、最大日流量多元回歸方程為Y3=80.681+14.055x0+7.340x1+3.300x2。

用以上方程，在2011~2014年涓水流域10次(2013年6月15~16日洪水過程時間相隔過近，得出方程效果較差，在此視為1次)洪水來臨前，用相應的日雨量實況對射埠站洪水洪峰高度和流量進行預測計算，結果發現計算結果與實況比較較為接近(表3)，可見，統計回歸方法對涓水流域洪峰、流量預報具有一定的指導意義，在面雨量預報具有一定精度的情況下，可提前2 d預報涓水射埠站洪峰流量。

表3　射埠站洪水流量、洪峰方程預報結果

這3個模型均能很好地預報出射埠站洪峰值，為甄選出最適合射埠站洪峰流量預報方程，采用Pearson相關系數來確定最優預報方程，由預報結果與實況值相關系數(表4)可見，對于涓水流域射埠水文站洪峰流量預報，多元回歸方法得出的結果與實際結果最為接近，洪峰流量的預報、實測數據的相關系數分別高達0.85、0.86；若由于數據缺失，無法采取多元回歸模型時，應采用最優因子逐步回歸模型(即步入1個因子，在此步入因子為前1 d面雨量)來預報，即y1、Y1的預報方程，其預報與實測數據相關性也較高。

表4　方程預報結果與實測值相關系數

3結論

通過分析射埠水文站洪峰過程與涓水流域降水量的關系，得出以下結論：

(1) 涓水流域出現累積面雨量達90.0 mm以上，且最大日面雨量≥100.0 mm時，其后1 d流域將有流量達1 000 m3/s以上洪峰產生；若當累積面雨量達60.0 mm以上時，其后1~2 d流域將有洪峰產生。

(2) 涓水流域為小型流域，預測洪峰時采用其前2 d的降水量，面雨量采用算術平均法計算。洪峰與其出現前1 d面雨量的關系最為顯著。

(3) 利用逐步回歸和多元回歸方法對洪峰期內射埠水文站水位、流量與各氣象站日降水量關系進行了研究，建立了洪峰預報方程模型。表明在站點資料齊全的情況下，采用多元回歸方法具有更好的預測準確率；若站點降水量資料不全，則采用最優因子逐步回歸模型。

參考文獻

[1] 肖立炎, 鄔清元, 鄧學友, 等.萍水河上游流域降水與城區河內洪峰關系的研究[J].江西氣象科技, 2006, 27(2):17-19.

[2] 王歡,單九生,鄧虹霞.信江流域洪峰期水位預測方法初探[J].暴雨災害,2007,26(4):352-354.

[3] 吳璐,靖春悅,張曉歐,等.沙澧河上游關鍵區汛期強降水面雨量對洪峰的影響[J].環境與科學,2010,33(3):29-32.

[4] 王銘才.大氣科學常用公式[M].北京：氣象出版社, 1994：518-519.

[5] 董官臣, 冶林茂.面雨量在天氣預報中的應用[J].氣象, 2000,26(1):9-13.

[6] 孟遂珍, 彭治班, 趙秀英，等.流域平均降水量的一種算法[J].北京氣象學院學報，2001(2)：64-68.

[7] 畢寶貴, 徐晶, 林建.面雨量計算方法及其在海河流域的應用[J].氣象, 2003, 29(8):39-42.

[8] 徐晶, 林建, 姚學祥, 等.七大江河流域面雨量計算方法及應用[J].氣象, 2001, 27(11):13-16.

[9] 何明瓊, 熊傳輝.清江流域、長江上游面雨量與致洪關系分析[J].四川氣象, 2000, 20(3):49-53.

收稿日期2015-12-09

作者簡介胡堅(1979- )，男，湖南永州人，工程師，碩士，從事氣象服務和應用氣象研究

基金項目2015年湖南省氣象局短平快課題(XQKJ15B141)。

中圖分類號S 164

文獻標識碼A

文章編號0517-6611(2016)02-214-02