鄭長陵 周沛勇
(新沂水文水資源監(jiān)測中心 221400)(邳州水文水資源監(jiān)測中心 221300)
沂河是沂沭泗水系中最大的山洪河道,發(fā)源于山東省沂源縣沂山,自江蘇邳州的齊村流入江蘇省徐州市境內(nèi),在新沂的苗圩附近匯入駱馬湖。河道總長353km(其中徐州境內(nèi)45.5km),流域總集水面積11820km2,在彭家道口由分沂入沭水道分水入沭河,在江風口由邳蒼分洪道分水入中運河。
20世紀90年代中后期,由于河床淤積、新建徐海公路壩、河灘栽種高稈植物等綜合因素影響,沂河徐州下游段(華沂—駱馬湖口段)的行洪能力下降到歷史最低水平。各級防汛部門高度重視,在東調(diào)南下沂河20年一遇治理工程(Ⅰ期工程,設(shè)計流量7000m3/s)、50年一遇治理工程(Ⅱ期工程,設(shè)計流量8000m3/s)中,先后采取拆壩建橋、中泓開挖等整治措施,使該河段行洪能力明顯增強。
為合理進行洪水調(diào)度,科學進行防洪決策,本文對東調(diào)南下沂河Ⅰ期工程、Ⅱ期工程的治理效果及沂河徐州下游段現(xiàn)狀行洪能力進行分析。
華沂(壩上)水位站位于沂河徐州境內(nèi)段中部,是沂河徐州下游段的重要水情控制站,該站水位—流量關(guān)系的演變直接反映沂河徐州下游段行洪能力的變化。因此,對該站不同時期水位—流量關(guān)系進行分析可以達到對該河段不同時期行洪能力進行分析的目的。
2000年初,曾就沂河徐州境內(nèi)段行洪能力進行過分析,得到1972~1998年間華沂(壩上)站的4種水位—流量關(guān)系。這4種水位—流量關(guān)系表明:在不到30年的時間里,該河段的行洪能力呈現(xiàn)明顯的三個階梯性降低[1]。
1999年汛后至2000年汛前,沂河徐州境內(nèi)段Ⅰ期治理工程實施了徐海公路壩拆除、堰頭以下段6.6km中泓開挖等項目。項目實施后,華沂(壩上)站呈現(xiàn)出2000~2008年水位—流量關(guān)系。2008年汛后至2009年汛前,沂河徐州境內(nèi)段Ⅱ期治理工程實施了華沂漫水壩閘改橋、華沂壩至鐵路橋段4.8km中泓開挖、沂河入湖口段4.2km中泓開挖等項目。項目實施后,華沂(壩上)站呈現(xiàn)出2009~2011年水位—流量關(guān)系。分析上述兩種水位—流量關(guān)系,即可得到這兩個時期該河段的河道行洪能力。將這兩種關(guān)系和20世紀中后期的水位—流量關(guān)系進行對比分析,即看出沂河階段性工程治理的防洪效果。
華沂(壩上)站1951~1971年為國家水文站,1972年流量測驗設(shè)施上遷20km至港上后降為水位站,從此該站僅有水位資料,要建立該站洪峰水位—洪峰流量關(guān)系,首先要解決其洪峰流量資料的來源問題。根據(jù)沂沭泗地區(qū)暴雨洪水往往具有較高同步性的水文特性,經(jīng)過歷史水文資料對比分析,結(jié)合防汛部門工作經(jīng)驗,綜合考慮洪水傳播特性和區(qū)間集水面積比重較小的實際,用上游港上水文站洪峰流量的90%推求華沂(壩上)站洪峰流量。
搜集2000~2011年沂河36場洪水資料,將2000~2008年的26場洪水組成第一組資料,將2009~2011年的10場洪水組成第二組資料,分別摘錄各場洪水華沂(壩上)站洪峰水位和港上站相應洪峰流量,用Qm華=0.90Qm港計算華沂站洪峰流量,分別繪制兩組資料 Gm~Qm相關(guān)圖[2](2000~2008年相關(guān)圖略,2009~2011年相關(guān)圖見圖1),選配的兩個相關(guān)線數(shù)學模型公式分別為:Qm=82.8187G2m-3549.3082Gm+36263.3399(2000~2008年)和 Qm=162.3336G2m-7843.0153Gm+94925.6439(2009~2011年),兩關(guān)系線均通過符號檢驗、適線檢驗、偏離數(shù)值檢驗,系統(tǒng)誤差和隨機不確定度均滿足一類精度巡測站單一曲線法水位—流量關(guān)系定線精度指標要求[3]。
圖1 華沂(壩上)站2009年至2011年水位—流量相關(guān)曲線
根據(jù)上述兩種關(guān)系,用得到的兩個數(shù)學模型公式,分別計算2000~2008年、2009~2011年各水位級對應的流量值。為便于比較分析,同時引用2000年分析成果中1994~1998年的水位流量關(guān)系[1],計算與引用成果匯總見下表。根據(jù)表中數(shù)據(jù)繪制華沂(壩上)站不同時期水位—流量關(guān)系對比分析圖,見圖2。
沂河華沂(壩上)站不同時期水位—流量關(guān)系表 單位:m3/s
圖2 華沂(壩上)站不同時期水位—流量關(guān)系對比分析
上表所列華沂(壩上)站三個時期的水位—流量關(guān)系,即代表著三個時期沂河徐州下游段的行洪能力,從上表和圖2可看出三種關(guān)系存在明顯的兩次突變,這兩次變化就是沂河20年一遇工程治理、50年一遇工程治理的明顯效果。
根據(jù)上述分析,沂河東調(diào)南下Ⅰ期、Ⅱ期工程治理階段性增加了沂河徐州下游段的河道行洪能力,大大減輕了沂河防汛壓力。
a.沂河20年一遇治理工程使徐州下游段河道行洪能力提高了15%左右,當發(fā)生4000m3/s左右的中等洪水時可使同級洪峰流量下的洪峰水位降低0.3m左右。
b.沂河50年一遇治理工程使徐州下游段河道行洪能力又提高了30%左右,當發(fā)生4000m3/s左右的中等洪水時可使同級洪峰流量下的洪峰水位再降低0.8m左右。
c.當沂河徐州下游段今后發(fā)生3500m3/s以下洪水時,2009年以后的水位—流量關(guān)系線可直接用作該河段洪水預報方案;當發(fā)生3500m3/s以上洪水時,可將新出現(xiàn)的水位—流量關(guān)系加入重新綜合定線,進行洪水預報方案修訂。也可根據(jù)沂沭泗洪水調(diào)度方案直接建立利用臨沂站洪峰流量預報華沂(壩上)站洪峰水位的預報方案。
d.建議流域機構(gòu)和地方政府加大河道管理力度,鞏固來之不易的豐碩成果。
1 徐州市青年科技工作者協(xié)會.徐州市第四屆青年學術(shù)年會論文集[C].徐州:[出版者不詳],2002.
2 蔡文祥,許大明.水文計算[M].南京:河海大學出版社,1989.
3 水利部長江水利委員會水文局.SL 247—1999水文資料整編規(guī)范[S].北京.中國水利水電出版社,2000.