王 平
(山西省忻州市水文水資源勘測分局,山西 忻州 034000)
大芳河地處代縣、原平交界,屬海河流域子牙河水系,為滹沱河一級支流。河流發(fā)源于山西省原平市沿溝鄉(xiāng)炭峪村西北,向東南經(jīng)代縣牛村、官莊、海子村、邵家寨,于原平市沿溝鄉(xiāng)河頭村匯入滹沱河。河源高程1859.6 m,河口高程829.3 m,域面積69.3 km2,河長21 km,河道平均比降25‰。本次劃界涉及大芳河段位于代縣境內(nèi),河段長7.17 km,上游流域面積17.9 km2,河段平均比降46.9‰。代縣屬典型的溫帶大陸性季風(fēng)氣候,四季分明,變化顯著,冬春較長,夏秋較短;冬季寒冷少雪,春季干燥多風(fēng),夏季炎熱,秋季涼爽。年平均氣溫7.2℃~9.3℃,平川為7.8℃~8.4℃,丘陵區(qū)7.2℃~9.0℃,土石山區(qū)6.4℃。平原區(qū)極端最低氣溫-24.5℃,極端最高氣溫35.9℃;山區(qū)極端最低氣溫-30℃,極端最高氣溫32℃。多年平均降雨量平川區(qū)442 mm,丘陵區(qū)500 mm,山區(qū)580 mm。雨量隨海拔高度的增高而增加。降雨量季節(jié)變化很大,年內(nèi)分配很不均勻,主要集中在7 月、8 月,6 月~9 月份占67%~77%,10 月~次年2月份占21%,3 月~5 月份占10%~17%。降水年際變化也極不均勻。最大年降水量為772.8 mm(1967 年),最小年降水量為244.9 mm(1965 年)。多年平均蒸發(fā)量1865.7 mm,干旱指數(shù)為3.2。忻州市大芳河河道炭峪村下游二廣高速- 邵家村東,規(guī)劃河長7.17 km,均位于忻州市代具境內(nèi)。
由于大芳河流域內(nèi)無水庫,屬于無水文觀測資料地區(qū)。依據(jù)《水利水電工程設(shè)計洪水計算規(guī)范》(SL 44-2006)和《山西省水文計算手冊》的有關(guān)規(guī)定,本次擬采用流域模型法[1]、推理公式法和經(jīng)驗公式法[2]進行設(shè)計洪水分析計算,并對計算成果進行比較,合理確定設(shè)計流域洪峰流量。
(1)設(shè)計點暴雨計算
代縣大芳河位于《山西省水文計算手冊》的水文分區(qū)中區(qū),計算10 min、60 min、6 h、24 h、3 d 點暴雨及其參數(shù)。由于該流域計算面積特小,因此,假定洪水與暴雨同頻率。本次規(guī)劃采用的設(shè)計洪水頻率為20 a 一遇,所以設(shè)計暴雨頻率采用20 a 一遇。由《山西省水文計算手冊》附圖查得控制斷面設(shè)計暴雨參數(shù)查圖成果見表1。
表1 大芳河控制斷面流域設(shè)計暴雨參數(shù)查圖成果表
通過合理性檢查滿足要求,由《山西省水文計算手冊》附錄表1-2 查得相對應(yīng)的模比系數(shù)KP,并按式Hp=KP·H 計算各歷時設(shè)計點暴雨量,并按下式求得暴雨公式的參數(shù)。
式中:Hp(t)為同頻率各標(biāo)準(zhǔn)歷時設(shè)計雨量;Sp為設(shè)計雨力,即1 h 設(shè)計雨量,mm/h;λ 為經(jīng)驗參數(shù);t 為暴雨歷時,h;n、ns為雙對數(shù)坐標(biāo)系中設(shè)計暴雨時~強關(guān)系曲線的坡度及t=1 h 時的斜率。
(2)設(shè)計面暴雨計算
由上式計算設(shè)計面暴雨量初值,然后求解暴雨參數(shù)λ,其值應(yīng)滿足0≤λ≤12。根據(jù)各歷時設(shè)計面暴雨量初值求出暴雨公式的三個參數(shù)Sp、λ、ns,求出不同歷時面暴雨量。設(shè)計暴雨計算成果見表2。
表2 P=5%設(shè)計暴雨計算成果
(1)設(shè)計凈雨深計算
依據(jù)《山西省水文計算手冊》計算設(shè)計凈雨深,計算結(jié)果見表3。
表3 大芳河控制斷面以上流域設(shè)計凈雨深計算成果表
(2)凈雨過程計算依據(jù)《山西省水文計算手冊》計算時段凈雨。把計算出的時段凈雨按序位編號安排在設(shè)計時雨型“模板”中的相應(yīng)序位位置,即得主雨日的凈雨過程。計算結(jié)果見表4。
表4 設(shè)計斷面P=5%主雨日凈雨過程計算結(jié)果 單位:mm
推理公式法由設(shè)計暴雨、推理產(chǎn)流、推理匯流三個子系統(tǒng)構(gòu)成,其中設(shè)計暴雨、推理產(chǎn)流與降雨徑流模型法相近。推理匯流計算包括求解洪峰流量Qm、流域匯流時間τ 及推理洪水過程線。在求解最大流量Qm時采用數(shù)值法。
根據(jù)下式聯(lián)立方程組求解。
式中:A 為流域面積,km2;L 為為河長,km;J 為河流縱比降,‰;m 為匯流參數(shù),根據(jù)流域下墊面的實際情況,從《山西省水文計算手冊》表7.3.1.2 中合理選用;hR,p為設(shè)計洪水凈雨深,mm。
推理公式法計算成果見表5。
表5 P=5%推理公式法計算成果
由設(shè)計暴雨推求相應(yīng)頻率設(shè)計洪峰流量的地區(qū)經(jīng)驗公式法,是一種集產(chǎn)流與匯流于一體的洪峰流量計算方法,多用于計算小匯水面積涉水工程的設(shè)計洪峰流量。
式中:Qp為頻率為P 的設(shè)計洪峰流量,m3/s;A 為涉水工程控制的流域面積,km2;Cp為與頻率p 和地類有關(guān)的經(jīng)驗參數(shù),查《山西省水文計算手冊》表7.3-6 取0.172;N1、β 為經(jīng)驗參數(shù),取N1=0.92,取β=0.050;Sp0為涉水工程控制流域內(nèi)定點概率雨力的面平均值,即設(shè)計定點雨力,mm/h。
P=5%推理公式法計算得Cp=30 min,Sp0=20,設(shè)計洪峰流量54.7 m3/s。
代縣大芳河流域河道治導(dǎo)線規(guī)劃及河道劃界設(shè)計洪水計算,根據(jù)涉水工程規(guī)模、流域水文條件、工程重要性、水文資料等按《山西省水文計算手冊》所提供的設(shè)計洪水計算方法進行選擇。流域模型法、推理公式法與經(jīng)驗公式法計算成果存在較大差別,見表6。流域水文模型法是由流域產(chǎn)流模型、流域匯流模型耦合而成,流域的產(chǎn)、匯流地類分別劃分,計算成果相對可靠。因此,代縣大芳河劃界采用流域水文模型法設(shè)計洪水計算成果,即設(shè)計洪峰流量20.1 m3/s。
表6 P=5%設(shè)計洪水成果比較
本文通過采用流域模型法、推理公式法和經(jīng)驗公式法進行代縣大芳河流域設(shè)計洪水分析計算成果比較的結(jié)果表明,流域水面線計算結(jié)果和實測大芳河控制斷面歷史洪水位調(diào)查成果吻合較好,所以認(rèn)為大芳河流域水文模型法設(shè)計洪水成果合理,能作為代縣大芳河流域河道治導(dǎo)線規(guī)劃、河道劃界及防洪排澇工程建設(shè)洪水設(shè)計的依據(jù)。