南城子水庫洪水序列展延和調(diào)洪演算

盧金江，寧 聰，張凱馨

(1. 鐵嶺市南城子水庫管理局，遼寧 鐵嶺 112304；2.河海大學(xué) 水文水資源學(xué)院，江蘇 南京 210098；3.南京水利科學(xué)研究院，江蘇 南京 210029)

設(shè)計洪水是確定水利水電工程建設(shè)規(guī)模及制定運行管理策略的重要依據(jù)[1]。推求設(shè)計洪水的常用方法主要是洪水頻率分析和 PMP/PMF 估算[2-4]。洪水頻率分析需要較長的連續(xù)洪水序列，若洪水序列較短或?qū)崪y期內(nèi)存在缺測年份，經(jīng)常通過相關(guān)分析對洪水序列進行插補延長。常用方法有利用上下游站資料插補、利用本站峰量關(guān)系插補和利用流域暴雨資料插補[5]。在南城子水庫設(shè)計洪水復(fù)核中，實測期洪水資料存在缺測年份，必須對洪水序列加以延長。本文通過建立南城子站天然洪水與實測洪水的相關(guān)關(guān)系推求缺測洪水，實現(xiàn)了實測期洪水序列的延長，以供設(shè)計洪水復(fù)核和調(diào)洪演算使用。

1 研究資料

1.1 水庫概況

南城子水庫位于遼寧省開原市威遠堡鎮(zhèn)南城子村，遼河三級支流葉赫河下游。水庫始建于1958年，主體工程于1965年竣工，是一座大(2)型水庫。1994年，南城子水庫進行了除險加固設(shè)計，加固后水庫防洪標(biāo)準(zhǔn)為100 a一遇設(shè)計,5000 a一遇校核。

轉(zhuǎn)山湖水庫位于南城子水庫上游41.58 km處,吉林省梨樹縣梨樹園子村附近。水庫建于1978年，控制面積121 km2，河長18.32 km，河道平均比降3.98‰，防洪標(biāo)準(zhǔn)為100 a一遇洪水設(shè)計，500 a一遇洪水校核。

1.2 水文資料

南城子水文站始建于1952年，目前擁有1952—1958年、1962—2012年系列的實測流量資料。此外，遼寧省水利水電勘測設(shè)計院分別在1966年及1987年用清河的東石人溝站資料對南城子站缺測年份洪水進行了插補，插補年份為1935—1945、1950、1951、1959—1961年。

轉(zhuǎn)山湖水庫對南城子水庫運行有很大影響，在設(shè)計洪水復(fù)核中不可忽視?，F(xiàn)有1994年南城子水庫除險加固設(shè)計時收集的轉(zhuǎn)山湖水庫1981—1990年水庫觀測資料，1991年以后該水庫資料缺失。該水庫調(diào)度運行方式未發(fā)生變化。

2 設(shè)計洪水復(fù)核

2.1 歷次洪水復(fù)核成果

1966年的《南城子水庫溢洪道技術(shù)設(shè)計》，水庫標(biāo)準(zhǔn)為100 a一遇洪水設(shè)計，1000 a一遇洪水校核。設(shè)計洪峰由手冊查算，七日洪量頻率計算采用系列為1935—1945年、1950—1960年22 a系列，1951年洪水按重現(xiàn)期120 a的第二位處理，100 a一遇洪峰流量為2810 m3/s，七日洪水總量1.69億m3；1000 a一遇洪峰流量為4430 m3/s，七日洪水總量2.52億m3[6]。

1994年的《南城子除險加固初步設(shè)計》，剔除了1951年之前插補的年份系列，采用洪峰洪量系列為1951—1990年系列，其中1951—1980年系列采用《修訂遼河流域規(guī)劃》的成果；同時考慮轉(zhuǎn)山湖水庫的影響,對1981—1990年的洪水進行了還原計算，1951年洪水按重現(xiàn)期147 a(1884—1990年)的第二位處理[7]，成果見表1。

表1 除險加固初步設(shè)計階段采用設(shè)計洪水成果

2.2 南城子洪水序列延長

本次設(shè)計洪水復(fù)核中1935—1945年，1950—1980年洪水系列采用《修訂遼河流域規(guī)劃》成果，1981—1990年洪水系列采用南城子水庫除險加固初步設(shè)計成果。

同時對南城子站1991—2012年洪水系列進行了展延，由南城子水庫觀測資料推求實測洪峰、洪量。由于缺乏1991—2012年轉(zhuǎn)山湖水庫觀測資料，1991—2012年南城子水庫天然洪水按1981—1990年南城子站天然洪水與實測洪水相關(guān)關(guān)系推求。通過點繪相關(guān)點距發(fā)現(xiàn)天然洪水與實測洪水間存在線性相關(guān)關(guān)系，建立一元線性回歸模型：

yi=a+bxi

(1)

(2)

(3)

(4)

式中:yi為天然洪水，m3/s;xi為實測洪水，m3/s;n為樣本數(shù)，各相關(guān)關(guān)系見圖1～圖4。

圖1 南城子天然洪水與實測洪水相關(guān)圖

圖2 南城子站天然1 d與實測1 d洪量相關(guān)圖

圖3 南城子站天然3 d與實測3 d洪量相關(guān)圖

圖4 南城子站天然3 d與實測3 d洪量相關(guān)圖

經(jīng)本次延長后南城子洪水系列為1935—1945、1950—2012年，其中1951年洪水按重現(xiàn)期170年(1844—2012年)的第二位處理。南城子洪水頻率計算考慮兩種方案：

(1)洪水系列選用1935—1945、1950—2012年共計74 a系列。

(2)洪水系列選用1951—2012年連續(xù)62 a系列，剔除1951年之前的插補年份。

經(jīng)頻率適線兩種情況計算的洪水水文參數(shù)相同。

比較表1和表2可知，由于近年無特大洪水，本次系列延長后，南城子水庫設(shè)計洪水小于1994年除險加固初步設(shè)計成果。為安全計，本次設(shè)計洪水維持水庫除險加固初步設(shè)計成果，見表2。

表2 本次復(fù)核南城子設(shè)計洪水成果表

2.3 設(shè)計洪水合理性分析

繪制南城子水庫及臨近各水庫的設(shè)計洪峰、洪量與流域面積、設(shè)計三日洪量與流域面積的雙對數(shù)相關(guān)圖，分析點據(jù)分布是否與流域暴雨及地形等因素相適應(yīng)[8]，結(jié)果見圖5～圖6。公式見式(5)～式(6)：

Q=aS2/3+b

(5)

W3=cS+d

(6)

式中：Q為設(shè)計洪峰，m3/s;W3為設(shè)計三日洪量，106m3;S為流域面積，km2;a,b,c,d為擬合系數(shù)。

可見各水庫設(shè)計洪峰基本與面積呈2/3次方關(guān)系，設(shè)計三日洪量與面積呈線性關(guān)系，符合洪水地區(qū)分布規(guī)律，表明南城子水庫設(shè)計洪水成果基本合理。

圖5 南城子及附近各水庫設(shè)計洪峰與面積關(guān)系圖

圖6 南城子及附近各水庫設(shè)計3 d洪量與面積關(guān)系圖

2.4 南城子設(shè)計洪水

現(xiàn)有1951年南城子洪水過程線和1953年南城子實測洪水過程。1951年洪水過程線由暴雨徑流關(guān)系推求，并與清河干流新立屯和區(qū)間洪水過程線組合到開原。其洪峰與開原站調(diào)查值相接近，過程線形狀也較相似。1994年水庫除險加固初設(shè)階段用1953年實測洪水過程為典型進行了比較，1953年洪水過程較矮胖，加倍放大以后峰型修正太多，過程線變形過大[7]。故最終采用1951年洪水過程線為典型，主峰按洪水調(diào)查值3090 m3/s修正。南城子100 a一遇及5000 a一遇設(shè)計洪水過程線見圖7。

圖7 南城子100 a一遇及5000 a一遇設(shè)計洪水過程線圖

3 調(diào)洪演算

汛限水位為149.00 m，調(diào)洪原則為采用溢洪道泄流，以水位控制為主，水位達到149.00 m，當(dāng)入庫流量小于溢洪道下泄流量時，以入庫流量為下泄流量；隨著入庫流量增加逐步加大溢洪道下泄流量直至達到其最大下泄流量。

表3 南城子水庫洪水調(diào)節(jié)成果表

由表3可見，本次復(fù)核成果與除險加固初步設(shè)計成果差別不大，本次復(fù)核的100 a一遇設(shè)計洪水位152.29 m，略高于原設(shè)計洪水位152.09 m(20 cm)，本次復(fù)核的5000 a一遇校核洪水位156.04 m，略低于原校核洪水位156.31 m(27 cm)。鑒于兩次復(fù)核結(jié)果相差較小，特征水位仍沿用除險加固成果，即100 a一遇設(shè)計洪水位為152.09 m，5000 a一遇校核洪水位為156.31 m[9]。

4 結(jié) 論

在缺乏上游轉(zhuǎn)山湖水庫觀測資料的條件下，本文提出建立1989—1990年南城子天然洪水和實測洪水的相關(guān)關(guān)系，利用南城子實測洪水插補1991—2012年南城子天然洪水序列的方法，并通過點繪洪峰、洪量與流域面積關(guān)系圖驗證設(shè)計洪水的合理性。利用延長后的洪水序列推求設(shè)計洪水和調(diào)洪演算，取得了良好的效果，證明該方法在洪水序列插補延長上行之有效。