基于不同類型指標的水庫汛期分期研究

孫 新

（撫順大伙房水庫發(fā)電有限責任公司，遼寧 撫順 113007）

大伙房水庫位于渾河中上游，座落在遼寧撫順境內(nèi)，控制流域面積5437 km2，是一座以防洪、灌溉和城市供水為主，兼顧發(fā)電和養(yǎng)殖的綜合利用大型水利工程[1]。

大伙房水庫設計汛限水位126.40 m，興利蓄水位131.50 m，設計洪水位136.00 m，可能最大洪水校核洪水位138.83 m?？値烊?2.68 億m3，興利庫容12.96 億m3，下游要求300年一遇洪水保證沈陽安全，100年一遇洪水保證撫順市安全，50年一遇洪水保證農(nóng)田安全[2-3]。水庫農(nóng)灌面積達150 萬畝，工業(yè)需水量5.9 億m3。防洪與興利矛盾十分突出，1960年棄水14 億m3，汛后蓄水卻比正常蓄水量少2 億m3，完全是因為汛限水位全期控制太低造成的；1964年連續(xù)發(fā)生四次較大洪水，又因過早抬高汛限水位使下游受淹。吸取歷史經(jīng)驗教訓，科學地進行汛期分期是必要的。

1 暴雨洪水的季節(jié)性變化規(guī)律

1.1 暴雨的時間分布特征

（1）逐旬降水量分析

圖1 多年平均逐旬降水量分布趨勢

結合水庫1962年～2019年共58年逐旬降水量統(tǒng)計資料，7 月下旬降水量最大，達到了90.7 mm，其次為7 月中旬和8 月上旬。在汛期中降水量主要集中在7 月上旬到8 月中旬，平均旬降水量67.5 mm，遠大于汛期平均降水量48.2 mm。

（2）25mm 以上降水量級分析

根據(jù)1962年～2019年共58年逐日降水量資料，統(tǒng)計日降水量在25 mm 以上的降水，并按照旬統(tǒng)計該降水量下的次數(shù)。由圖2 可知，大于25 mm 降水主要集中在7 月上旬～8 月中旬，達到了216 次，占整個汛期的70%。其中7 月下旬最大，達到了64 次，占整個汛期的21%。

圖2 大于25mm 以上降水次數(shù)分布趨勢

（3）1～3 日最大降水量分析

根據(jù)1962年～2019年共58年逐日降水量資料，統(tǒng)計6月～9 月范圍內(nèi)1～3 日最大降水量次數(shù)。由圖3 可知，1～3 日最大降水量集中分布在7 月下旬到8 月中旬，7 月下旬1～3 日最大降水次數(shù)平均在20 次左右，8 月上旬1～3 日最大降水次數(shù)平均近15 次，8 月中旬1～3 日最大降水次數(shù)平均近10 次，其中7 月下旬中2 日內(nèi)最大降水量次數(shù)達到了近25 次。通過對58年的1～3 日最大降水量統(tǒng)計分析，1 日、2 日、3 日的最大降水量分別為159.2 mm、222 mm、270 mm。

圖3 1～3 日最大降水量次數(shù)分布趨勢

根據(jù)降水時空分布的分析，降水具有明顯的季節(jié)性分布規(guī)律，整個汛期分為前汛期、主汛期和后汛期三個階段。逐旬降水量指標的主汛期發(fā)生在7 月上旬～8 月中旬；25mm 以上降水量級和1～3 日最大降水量兩項指標的主汛期一致，發(fā)生在7 月下旬到8 月中旬。

1.2 洪水的時間分布特征

（1）逐旬流量時間分布

根據(jù)1972年～2019年6 月～9 月逐旬入庫流量資料，繪制了逐旬入庫流量變化趨勢線。同時與逐旬降水量進行對比分析。由圖5 可知，在整個汛期階段7 月下旬～8 月下旬水庫的入庫流量達到了最大值，均在逐旬平均流量以上，其中7 月下旬～8 月中旬的逐旬流量大于120 m3/s（汛期旬降水量×徑流系數(shù)）。從發(fā)生的次數(shù)來看，旬流量大于220 m3/s 次數(shù)發(fā)生在7 月下旬到8 月上旬。從降水量與流量兩項指標的分析結果來看，基本是一致的，降水量、入庫流量的最大值均是發(fā)生在整個汛期的7 月下旬～8 月下旬，但是入庫流量略有滯后。

圖4 旬降雨量與旬入庫流量對比圖

圖5 大于1000 m3/s 洪峰流量分布趨勢

（2）大于1000 m3/s 的洪峰流量時間分布

根據(jù)洪水流量資料，統(tǒng)計了大于1000 m3/s 入庫流量的頻次、峰值與日期。由圖6 可知，大于1000 m3/s 洪水主要集中在7月～8 月，其中近90%的洪水分布在7 月下旬～8 月中旬，大于2000 m3/s 洪水主要分布在7 月下旬～8 月中旬，最大的峰值發(fā)生在1995年7 月30 日，流量達到了10700 m3/s。

（3）最大洪峰流量分析

統(tǒng)計1972年～2019年共48年大伙房水庫峰值及時間，結果見表1。由表1 可知，最大峰值出集中在7 月下旬～8 月中旬，所占比重達到75.5%，其中最大的三場洪水分別為1995年7 月30 日，峰值10700 m3/s；1960年8 月4 日，峰值7630 m3/s；2005年8 月13 日，峰值6040 m3/s。

表1 1972年～2019年最大洪峰發(fā)生時間統(tǒng)計

根據(jù)洪水時空分布的分析，洪水具有明顯的季節(jié)性分布規(guī)律，整個汛期分為前汛期、主汛期和后汛期三個階段。逐旬流量指標的主汛期發(fā)生在7 月中旬～8 月下旬；洪峰流量與最大洪峰兩項指標的主汛期一致，發(fā)生在7 月下旬～8月中旬。

2 大伙房水庫模糊統(tǒng)計法檢驗汛期

為確定汛期分期的合理性，分別采用日降雨量、日入庫流量為指標模糊統(tǒng)計法進行檢驗。模糊統(tǒng)計法是一種通過描述汛期的模糊性，即非汛期與汛期之間的中介過渡性來進行汛期分期的方法。通過利用主汛期隸屬函數(shù)來刻畫時間對主汛期的隸屬程度，以實現(xiàn)時間序列的聚類劃分。在給定進入汛期的指標降雨量、徑流量閾值后，依據(jù)多年實測資料，統(tǒng)計出時間隸屬于汛期的隸屬度，最后對擬合后的隸屬函數(shù)取一定的閥值，將超過閥值的部分作為主汛期，之前作為前汛期，之后作為后汛期。

從氣象降雨分析，9 月底以后、6 月1 日前的降大雨（大于25 mm）可能性及日入庫流量（汛期平均降雨量乘以徑流深折合成的流量）均很小，為此以6 月～9 月為汛期變化范圍，確定汛期與非汛期界線的指標為YT徑流=120 m3/s、YT降雨=25 mm/d。根據(jù)樣本數(shù)據(jù)，當大于等于YT的起止時間來確定汛期范圍。經(jīng)計算分析，日降水量指標的主汛期發(fā)生在7 月10 日～8 月19 日，日入庫流量指標的主汛期發(fā)生在7 月19 日～8 月17日，兩項指標隸屬度變化曲線見圖6。

圖6 降雨量和入庫流量隸屬度曲線對比

由圖6 可知，日降水量、日流量作為分析指標，隸屬度均為大于等于0.7 作為劃入主汛期，結果顯示日降水量作為指標，7 月10 日～8 月19 日間進入主汛期；日流量作為指標，7月19 日～8 月17 日間進入主汛期。

綜合以上的分析，兩項指標的主汛期時間是基本一致的，只是起止時間略有差異，但是考慮到主汛期來臨時洪水流量較大以及結束時安全等因素，最終確定主汛期為7 月11 日～8 月20 日。

3 大伙房水庫汛期分期結果與分析

大伙房水庫采用數(shù)理統(tǒng)計、模糊統(tǒng)計兩種方法分期結果見表2。

表2 水庫汛期分期成果統(tǒng)計

由表2 可知，在數(shù)理統(tǒng)計法中采用旬降水量、25 mm 以上降水量級、1～3 日最大降水量、旬流量和洪峰流量為指標研究汛期分期時，主汛期結束的時間都是8 月中旬。

當旬降水大于汛期旬降水量、25 mm 以上降水量級次數(shù)占汛期的大于11%、1～3 日最大降水量次數(shù)占汛期的大于15%、年最大洪峰流量次數(shù)占汛期的大于15%，均作為入主汛期考慮。

在模糊統(tǒng)計法中采用日降水量和日流量為指標，當日降水量與日流量大于等于旬降水量與徑流系數(shù)之積的流量值時，隸屬度值大于等于0.7，進入主汛期，也說明了數(shù)理統(tǒng)計法以各項指標的汛期劃分。

根據(jù)水庫的多年運行實際狀況，在每年的8 月20 日之后抬升汛限水位后，普遍在接下來的月份來水量較少，有效庫容很難蓄滿，達不到正常蓄水位，因此設置過渡階段，根據(jù)降水量的多少逐漸抬升水位。采用洪峰流量指標，由于洪水的滯后性，主汛期為7 月下旬。

根據(jù)兩種方法及實際運行情況的分析，確定水庫主汛期為7 月21 日～8 月20 日。

4 結論

通過大伙房水庫汛期分期研究，汛期分期成果對大伙房水庫汛期防洪、調蓄供水具有一定的指導和依據(jù)，使水庫發(fā)揮了最佳效益。由于水庫汛期分期受到的影響因素眾多，汛期分期結果與水庫實際運行時，主汛期可能會出現(xiàn)提前或滯后，還需要根據(jù)當時的各種條件做相應的微調，但從總的趨勢與范圍是基本沒有變化的，進一步說明了汛期分期成果是合理的。隨著氣候、降雨及徑流等因素的變化，汛期分期要做定期的復核，同時管理部門要加強實際運行時水庫的調度。為大伙房水庫防汛、調蓄調度相協(xié)調且發(fā)揮最大效益提供有意義的指導。