劉田田,楊明非,孫 瑞,房 磊,俞琳琳
江蘇省水文水資源勘測(cè)局 徐州分局,江蘇徐州 221000
警戒水位、保證水位等防汛特征水位是各級(jí)政府防汛決策指揮的重要依據(jù)。運(yùn)河站是中運(yùn)河重要水情控制站,其現(xiàn)有防汛特征水位是徐州市防汛抗旱指揮部1974年大洪水后,根據(jù)當(dāng)時(shí)的水情、工情及防洪能力等確定的。由于目前工情水情的改變,該水位已多年未超警,失去了原有的警戒作用。
中運(yùn)河運(yùn)河站預(yù)報(bào)區(qū)間為微山湖、會(huì)寶嶺水庫(kù)、江風(fēng)口閘以下、運(yùn)河站以上區(qū)域,集水面積6 102 km2,其中包括微山湖湖西地區(qū)326 km2,藺家壩閘至解臺(tái)閘區(qū)間262 km2,解臺(tái)閘至劉山閘區(qū)間870 km2,韓莊閘至臺(tái)兒莊閘區(qū)間1 345 km2,林子站以上流域面積2 085 km2。本流域內(nèi)有江蘇省的邳州市和山東省的蘭陵縣,故俗稱邳蒼地區(qū)[1]。
韓莊閘至運(yùn)河站河長(zhǎng)72.4 km,其中,韓莊閘至魯蘇邊界河段稱韓莊運(yùn)河,河長(zhǎng)42.5 km;魯蘇邊界至駱馬湖為中運(yùn)河,魯蘇邊界至運(yùn)河站河長(zhǎng)29.9 km。左岸有嶧城大沙河、陶溝河、邳蒼分洪道等較大支流匯入,右岸有伊家河于臺(tái)兒莊鎮(zhèn)西南匯入。藺家壩閘至中運(yùn)河河口稱京杭運(yùn)河,是不牢河取直而成,河長(zhǎng)73 km,其中藺家壩閘至解臺(tái)閘河長(zhǎng)26 km,解臺(tái)閘至劉山閘河長(zhǎng)41 km。林子站至運(yùn)河站26.3 km,臺(tái)兒莊閘至運(yùn)河站37.4 km。
運(yùn)河站以上流域北部為山區(qū),坡度大,來(lái)水水流湍急。邳州市北部地勢(shì)低洼,運(yùn)西地區(qū)為平原,來(lái)水水勢(shì)較為平緩。韓莊運(yùn)河河寬400~750 m,蘇魯省界至運(yùn)河鎮(zhèn)河寬1 300~71 800 m。
運(yùn)河站位于徐州市邳州市運(yùn)河鎮(zhèn)前索家村,距邳州市主城區(qū)較近,設(shè)立于1950年。掌握中運(yùn)河進(jìn)入駱馬湖水量、沙量,以及調(diào)水水量和水質(zhì)情況,為水文預(yù)報(bào)、分析計(jì)算,以及防汛、防旱、規(guī)劃治理、水資源管理收集資料。
運(yùn)河站為國(guó)家基本水位站,水位等各項(xiàng)水文資料的測(cè)驗(yàn)和整編滿足各時(shí)期的相關(guān)測(cè)驗(yàn)規(guī)范和水文資料整編規(guī)范的要求,歷年來(lái),基本資料均通過(guò)主管部門的審查。因此,資料成果可靠、質(zhì)量較高,可用于防汛特征水位復(fù)核調(diào)整計(jì)算。
運(yùn)河水文站每年汛前對(duì)所有水尺進(jìn)行校測(cè),汛期使用過(guò)的水尺汛后再校測(cè)一次;基本水準(zhǔn)點(diǎn)每年檢查一次,發(fā)現(xiàn)異常情況時(shí)及時(shí)校測(cè),逢0、5年份必須校測(cè)一次。因此,資料系列具有一致性,可用于防汛特征水位復(fù)核分析計(jì)算。
運(yùn)河水文站觀測(cè)時(shí)間較長(zhǎng),已包含完整的豐平枯周期,故可認(rèn)為運(yùn)河水文站1974—2020年以來(lái)的水位系列具有較好的代表性,可用于分析防汛特征水位。
3.1.1 原警戒水位存在問(wèn)題 警戒水位是指江河漫灘行洪,堤防可能發(fā)生險(xiǎn)情,需要開(kāi)始加強(qiáng)防守的水位。運(yùn)河站現(xiàn)有防汛特征水位是徐州市防汛抗旱指揮部在1974年根據(jù)當(dāng)時(shí)的水情、工情確定的,中運(yùn)河運(yùn)河站保證水位26.5 m,相應(yīng)行洪流量不超過(guò)6 500 m3/s,警戒水位25.5 m。
由于1974年大洪水的洪峰水位、洪峰流量迄今為止仍為新中國(guó)成立以來(lái)的歷史最高。近年來(lái),隨著水情、工情的變化,區(qū)域防洪能力已有較大提高,該水位明顯偏高。運(yùn)河站近年來(lái)水位偏低,汛期經(jīng)常低水位行洪,近十幾年水位沒(méi)有一次超警,該水位已失去應(yīng)有的警戒作用。
3.1.2 距平法對(duì)年最高水位的分析 距平法樣本為年極值取樣,一組數(shù)據(jù)某一要素與多年平均值的差即該要素Hi的距平值。
對(duì)運(yùn)河水文站1974—2020年資料系列進(jìn)行距平分析,年最高水位距平、累計(jì)距平分析見(jiàn)圖1。根據(jù)最高水位距平分析圖可知,運(yùn)河水文站歷年最高水位在0軸附近上下波動(dòng),波動(dòng)符合正常規(guī)律。從累積距平圖看,1974—1989、2008—2020年的年最高水位處于回落階段,2002—2008年的年最高水位處于抬升階段。
圖1 運(yùn)河站年最高水位距平、累計(jì)距平分析圖
運(yùn)河站水位、流量1989年以前總體偏低,波動(dòng)幅度較?。?989以年后為控制駱馬湖水資源,人工干預(yù)增多,1989—2002年的水位變幅增大,流量變化幅度相應(yīng)增加;2002年中運(yùn)河沿線配套建筑及影響工程完工,致使2002—2020年運(yùn)河水文站流量明顯變高,變幅也隨之增加;2008年河道擴(kuò)挖,堤防加固,同時(shí)處理了多個(gè)穿堤涵洞,運(yùn)河沿線防洪保護(hù)區(qū)的防洪標(biāo)準(zhǔn)由20年一遇提高到50年一遇,2008—2020年的水位較前一時(shí)段低。
考慮到運(yùn)河站近年來(lái)工程現(xiàn)狀,分別采用1974—2020、1974—2002和2003—2020年3個(gè)資料系列進(jìn)行分析。
3.1.3 超不同等級(jí)水位的天數(shù)分析 以0.10 m為統(tǒng)計(jì)水位段,統(tǒng)計(jì)日平均水位超不同等級(jí)的天數(shù),運(yùn)河站1974—2020年不同系列不同水位級(jí)出現(xiàn)天數(shù)分布見(jiàn)圖2。
圖2 超不同等級(jí)水位出現(xiàn)天數(shù)分布圖
分析運(yùn)河水文站逐日平均水位可 知,1974—2020年的日平均水位多年平均值為22.56 m,歷史最高水位26.42 m(1974年8月15日),年水位超25.50 m(現(xiàn)警戒水位)的情況僅在3年中發(fā)生過(guò),分別為1974、1993、2005年。
從日平均水位超不同等級(jí)天數(shù)的統(tǒng)計(jì)分析可知,在1974—2020年,超過(guò)25.50 m的天數(shù)僅為0.11 d/年;在1974—2002年中,超過(guò)25.50 m的天數(shù)為0.17 d/年;在2003—2020年,超 過(guò)25.50 m的天數(shù)為0。
根據(jù)運(yùn)河站1974—2020年的日平均水位資料系列,以0.10 m為統(tǒng)計(jì)水位段進(jìn)行日平均水位超不同等級(jí)的概率分析,統(tǒng)計(jì)累積概率見(jiàn)圖3。
圖3 運(yùn)河站不同系列累計(jì)概率分布圖
分析日平均水位超不同等級(jí)的累積概率,在1974—2020年,運(yùn)河站日平均水位超過(guò)26.00 m的累積概率為0.65%,超過(guò)25.5 m的累積概率為1.08%。
分階段統(tǒng)計(jì)看,在1974—2002年,超過(guò)25.50 m的累積概率為2.16%;在2003—2020年,超過(guò)25.50 m的累積概率為0。
3.1.4 年最高水位頻率分析 年最高水位頻率分析采用P-Ⅲ型頻率曲線法,通過(guò)計(jì)算年頻率對(duì)應(yīng)的年最高水位,在頻率曲線上查到相應(yīng)的水位值[2]。
對(duì)1974—2020年的逐年最高水位由大到小進(jìn)行排頻,由頻率計(jì)算程序進(jìn)行適線分析,頻率曲線見(jiàn)圖4,分析結(jié)果見(jiàn)表1。
表1 運(yùn)河站年最高水位頻率分析成果表
圖4 頻率曲線圖
根據(jù)頻率分析成果,2年一遇、5年一遇、10年一遇水位為分別23.8 m,24.5 m,25.0 m,現(xiàn)狀警戒水位25.5 m,相當(dāng)于22年一遇。
3.1.5 運(yùn)河站與駱馬湖防汛特征水位互相關(guān)系 中運(yùn)河為復(fù)式河床斷面,小水時(shí)經(jīng)河槽排水,大水時(shí)由灘地行洪。其中,自省界至大王廟長(zhǎng)11.8 km,河底高程17.0 m,底寬50~80 m;大王廟至運(yùn)河鎮(zhèn)長(zhǎng)17.1 km,河底高程18.0~16.6 m,河 底 寬80~100 m;運(yùn)河鎮(zhèn)至房亭河口長(zhǎng)6.3 km,河底高程16.5 m,底寬100~140 m;房亭河口至窯灣長(zhǎng)12.4 km,河底高程16.5 m,河底寬150~200 m;窯灣至二灣,河底高程16.0~16.5 m,寬160~130 m。為行洪需要,房亭河口以下還開(kāi)挖了淺槽,至窯灣長(zhǎng)12.4 km,河底高程19.5 m,底寬170~180 m。中運(yùn)河運(yùn)河鎮(zhèn)以下主槽糙率為0.022 5,淺槽糙率為0.025,灘面糙率為0.033。
中運(yùn)河主要節(jié)點(diǎn)流量水位成果見(jiàn)表2。
表2 中運(yùn)河主要節(jié)點(diǎn)流量水位成果表
考慮運(yùn)河鎮(zhèn)行洪5 500 m3/s,河床糙率采用已有成果,參考運(yùn)河主要水位節(jié)點(diǎn)成果表,駱馬湖按現(xiàn)警戒水位為23.50 m起算,運(yùn)河站水位約為25.50 m。
3.1.6 警戒水位確定 根據(jù)頻率分析結(jié)果,運(yùn)河站原警戒水位25.5 m的重現(xiàn)期約為22年一遇。以該水位計(jì)算,自1974年以來(lái),平均每年有0.11 d超警戒,超25.5 m的累積概率為1.08%。據(jù)統(tǒng)計(jì),1974—2018年,僅 有1974、1993、2005年最高水位超過(guò)警戒水位。
根據(jù)運(yùn)河水文站1974—2020年運(yùn)河日平均水位超不同等級(jí)天數(shù)分析,2003—2020年運(yùn)河站水位下降趨勢(shì)明顯,超過(guò)25.5 m的水位于1974—2002年。
警戒水位應(yīng)根據(jù)歷史洪水重現(xiàn)期,選擇適宜的超警天數(shù)。根據(jù)頻率分析結(jié)果,運(yùn)河站水位按2~5年一遇考慮較為合適,警戒水位應(yīng)設(shè)置在23.8~24.5 m之間。
警戒水位要考慮河段普遍漫灘和重要堤段臨水,并達(dá)到一定高度,結(jié)合工程現(xiàn)狀,堤防工程歷史出險(xiǎn)情況等予以綜合研究確定。2008年韓莊運(yùn)河、中運(yùn)河及駱馬湖堤防工程基本實(shí)施完成后,運(yùn)河沿線防洪保護(hù)區(qū)的防洪標(biāo)準(zhǔn)由20年一遇提高到50年一遇,防洪能力有了較大的提高,但仍有局部堤防堤段受運(yùn)河河勢(shì)影響,水流頂沖,歷年是防汛的重點(diǎn)地段,運(yùn)河水位過(guò)高會(huì)對(duì)邳蒼分洪道行洪產(chǎn)生很大壓力,而且根據(jù)運(yùn)河站最新大斷面,當(dāng)運(yùn)河站水位高于24.53 m時(shí),左岸開(kāi)始漫灘行洪,考慮到運(yùn)河站部分堤防相對(duì)薄弱,警戒水位不宜超過(guò)該水位。綜合以上分析,運(yùn)河站警戒水位設(shè)定為24.5 m。
保證水位是指保證堤防及其附屬工程安全擋水的上限水位。
堤防的高度、寬度、坡度及堤身、堤基質(zhì)量已達(dá)到規(guī)劃設(shè)計(jì)標(biāo)準(zhǔn)的河段,其設(shè)計(jì)洪水位即為保證水位。堤防工程尚未達(dá)到規(guī)劃設(shè)計(jì)標(biāo)準(zhǔn)的河段,可按安全防御相應(yīng)的洪水位確定。例如,堤防工程防洪能力能夠防御超規(guī)劃設(shè)計(jì)標(biāo)準(zhǔn)的歷史最高洪水位,保證水位可按歷史最高洪水位確定[3]。
新中國(guó)成立以來(lái),中運(yùn)河多次進(jìn)行治理,2008年實(shí)施了沂沭泗河洪水東調(diào)南下韓莊運(yùn)河、中運(yùn)河及駱馬湖堤防加固工程。工程完成后,南四湖和運(yùn)河沿線防洪保護(hù)區(qū)的防洪標(biāo)準(zhǔn)由20年一遇提高到50年一遇。因此,根據(jù)中運(yùn)河的堤防實(shí)際現(xiàn)狀,確定運(yùn)河水文站的保證水位仍為26.5 m,與《淮河流域綜合規(guī)劃》中運(yùn)河的50年一遇防洪水位26.5 m相適應(yīng)。
綜合中運(yùn)河現(xiàn)狀灘地高程調(diào)查成果,結(jié)合超不同等級(jí)水位的天數(shù)及最高水位頻率分析,警戒水位此次復(fù)核將運(yùn)河站警戒水位由25.5 m下降至24.5 m,保證水位仍為26.5 m不變。