黃碧波,鄧 懿,鐘 波
(四川省成都水文水資源勘測(cè)局,四川 成都 611130)
干五里河屬長(zhǎng)江流域岷江水系支流,總河長(zhǎng)22 km,落差1874 m,河道平均比降5.8‰,流域面積108 km2。干五里河流經(jīng)四川省崇州市三郎鎮(zhèn),有和平水文站控制流域面積78.7 km2。干五里河流域是典型的山區(qū)型暴漲暴落洪水特性,以暴雨洪水為主,河道調(diào)蓄能力小,坡降較陡,沿途接納來(lái)自眾多山溪洪流,匯流快,洪水持續(xù)時(shí)間短,但漲幅大,洪峰高,洪水過(guò)程呈多峰尖瘦型。災(zāi)害特征以汛期多發(fā)山洪災(zāi)害、地質(zhì)災(zāi)害為主,對(duì)沿河低凹地帶的鄉(xiāng)鎮(zhèn)及居民生命財(cái)產(chǎn)安全造成重大威脅。
暴雨歷時(shí)分析根據(jù)流域大小和產(chǎn)匯流特性,確定小流域設(shè)計(jì)暴雨所需要考慮的最長(zhǎng)暴雨歷時(shí)及其典型歷時(shí)。常規(guī)設(shè)計(jì)暴雨時(shí)段有10分鐘、1小時(shí)、6小時(shí)、24小時(shí)4種標(biāo)準(zhǔn)歷時(shí),根據(jù)干五里河流域暴雨特征,本次暴雨歷時(shí)采用前3種標(biāo)準(zhǔn)歷時(shí)計(jì)算,計(jì)算暴雨的頻率為5年一遇、10年一遇、20年一遇、50年一遇、100年一遇5種。根據(jù)《四川省暴雨等值線圖集》中干五里河流域等值線分布情況,可查得小流域 4個(gè)時(shí)段(10min、1h、6 h和 24 h)暴雨量和變差系數(shù)[1]。
根據(jù)流域設(shè)計(jì)暴雨參數(shù),分別求出10 min、1 h、6 h、24 h的面暴雨,計(jì)算結(jié)果見(jiàn)表1。
表1 干五里河流域設(shè)計(jì)暴雨成果表(mm)
2.1.2 P-Ⅲ頻率曲線法
根據(jù)干五里河實(shí)測(cè)降雨資料,采用P-Ⅲ頻率曲線進(jìn)行適線得到結(jié)果,與干五里河流域和平站實(shí)測(cè)降水資料排頻對(duì)比,設(shè)計(jì)頻率差值見(jiàn)表2。
表2 干五里河流域設(shè)計(jì)暴雨成果表(mm)
由表2可知,10分鐘設(shè)計(jì)暴雨的成果誤差最大,依次到24 h的設(shè)計(jì)暴雨成果誤差最小。
設(shè)計(jì)洪水計(jì)算包括設(shè)計(jì)暴雨計(jì)算、產(chǎn)流(凈雨)計(jì)算、匯流計(jì)算3個(gè)環(huán)節(jié)。中小流域設(shè)計(jì)洪水計(jì)算方法采用以下兩種計(jì)算方法[2]:
(2)在有實(shí)測(cè)洪水資料的地區(qū),可采用頻率適線法推求出洪水特征值(均值、變差系數(shù)CV、偏差系數(shù)CS后,計(jì)算出不同設(shè)計(jì)標(biāo)準(zhǔn)的洪峰流量、時(shí)段流量,其洪水設(shè)計(jì)過(guò)程線采用典型洪水過(guò)程放大法。
2.2.1 推理公式法
采用《四川省中小流域暴雨洪水手冊(cè)》(下簡(jiǎn)稱“手冊(cè)”)推薦的推理公式法計(jì)算小流域的設(shè)計(jì)洪水,其公式如下:
式中:Q為最大流量(m3/s);Ψ為洪峰徑流系數(shù);S為暴雨雨力,即最大1 h暴雨量(mm/h);τ為流域匯流時(shí)間(h);n為暴雨公式指數(shù);F 為集水面積(km2)。
2.2.2 瞬時(shí)單位線法
根據(jù)《四川省中小流域暴雨洪水計(jì)算手冊(cè)》中有關(guān)綜合瞬時(shí)單位線的分區(qū),結(jié)合各小流域的地理分布情況,分別確定各小流域的各類分區(qū),按各分區(qū)的參數(shù)進(jìn)行洪水分析計(jì)算。設(shè)計(jì)洪水流量過(guò)程由用綜合瞬時(shí)單位線推求的設(shè)計(jì)地面徑流過(guò)程再加上由穩(wěn)定入滲量形成的地下徑流過(guò)程和深層地下水補(bǔ)給的基流量組成。
2.2.3 實(shí)測(cè)資料計(jì)算法
和平水文站共有35年(1981-2015年)實(shí)測(cè)系列洪水資料,加入調(diào)查到的歷史洪水資料組成一個(gè)不連續(xù)洪峰流量系列采用P—Ⅲ型理論頻率曲線進(jìn)行適線,目估適線法確定統(tǒng)計(jì)參數(shù)并推求設(shè)計(jì)值。和平水文站設(shè)計(jì)洪水成果詳見(jiàn)表3。
夏日里,蔣海峰回老家,帶水仙芝一起游雙峰山,放飛心情,向山谷大聲呼喚,向藍(lán)天白云傾吐了愛(ài)意。他倆手牽著手,穿過(guò)茂密的竹林,并肩躺在青山之上。
表3 和平水文站設(shè)計(jì)洪水成果表
通過(guò)3種方法計(jì)算結(jié)果比較分析:推理公式法和綜合瞬時(shí)單位線計(jì)算成果與實(shí)測(cè)均有一定的差異,其中推理公式法成果較實(shí)測(cè)值成果偏小11.9%~18.9%;綜合瞬時(shí)單位線計(jì)算成果較實(shí)測(cè)值成果偏大26.7%~33.9%。以實(shí)測(cè)資料為真值,綜合分析來(lái)看,小流域分析計(jì)算設(shè)計(jì)洪水時(shí),推理公式相比瞬時(shí)單位線法計(jì)算結(jié)果更接近實(shí)測(cè)值,因此在無(wú)資料小流域建議采用暴雨洪水推理公式計(jì)算方法及成果。
API(前期雨量指數(shù))是一種傳統(tǒng)的用于流域產(chǎn)流預(yù)報(bào)的方法。由于方法主要針對(duì)獨(dú)立的次洪水產(chǎn)流量計(jì)算,故稱其為次洪API模型(配合單位線法預(yù)報(bào)匯流)。1969年,美國(guó)天氣局西納(Sittner)等在“用API型水文模型綜合連續(xù)過(guò)程線法”一文中,將API方法與單位線法結(jié)合,發(fā)展成為一個(gè)完整地模擬降雨徑流過(guò)程的流域模型,稱其為連續(xù)的API流域模型。
根據(jù)干五里河流域水情站點(diǎn)資料系列情況?;贏PI預(yù)報(bào)模型原理主要選取鉆洞水文站和和平水文站36年(1981~2016年)的水文資料。
API模型主要有4個(gè)基礎(chǔ)參數(shù):基流、土壤最大含水量Im、Em和K值、P-Pa-R曲線。
(1)基流:取和平水文站年最小流量多年均值作為干五里河基流值,值為0.278 m3/s。
(2)土壤最大含水量Im:選取資料時(shí)按照前期久旱不雨,突然發(fā)生大暴雨的資料,計(jì)算該流域平均雨量P及其所產(chǎn)生的徑流深R,依據(jù)流域水量平衡方程式求得最大初損Im=70 mm。
(3)Em和K值詳見(jiàn)表4。
表4 Em和K值
(4)P-Pa-R曲線如圖1。
圖1 P-Pa-R曲線
和平水文站控制流域面積78.7 km2,單位線時(shí)段長(zhǎng)取1小時(shí),單位線凈流深 R=10 mm,單位線為:0、0.5、28.5、7、0.14、0.15、0.1、0。
根據(jù)干五里河小流域前期資料統(tǒng)計(jì)情況,利用API模型,采用和平站1989~2009年的資料進(jìn)行率定,2010年至2016年共7年的資料用于評(píng)定檢驗(yàn)。采用API模型對(duì)2010~2016年共20次洪水樣本進(jìn)行評(píng)定,結(jié)果為:洪峰流量的合格率為60%,洪峰出現(xiàn)時(shí)間的合格率為80.0%。
利用和平水文站以上小流域設(shè)計(jì)洪水資料,采用推理公式法、綜合瞬時(shí)單位線法與和平站實(shí)測(cè)洪水法進(jìn)行計(jì)算,對(duì)所得成果進(jìn)行比對(duì)可知,小流域分析計(jì)算設(shè)計(jì)洪水時(shí),推理公式相比瞬時(shí)單位線法計(jì)算結(jié)果更接近實(shí)測(cè)值。干五里河流域產(chǎn)流預(yù)報(bào)采用API流域模型,主要參數(shù)是基流、土壤最大含水量Im、Em和K值、P-Pa-R曲線,對(duì)其進(jìn)行精度評(píng)定,通過(guò)20次洪水樣本確定洪峰流量的合格率為60%,洪峰出現(xiàn)時(shí)間的合格率為80%。
隨著中小河流站點(diǎn)建設(shè),建議優(yōu)化合理的布置上游的雨量站點(diǎn),逐步補(bǔ)充資料,完善預(yù)報(bào)方案,有助于實(shí)現(xiàn)中小河流流域的水情自動(dòng)預(yù)報(bào),為產(chǎn)流預(yù)報(bào)提供有效的數(shù)據(jù)監(jiān)測(cè)依據(jù)。