黃志強,陳思翌,閆團進,季亞輝
(1.湖北省水利水電科學(xué)研究院,武漢 430070;2.陜西省引漢濟渭工程建設(shè)有限公司,西安 710100)
中國山洪災(zāi)害每年都造成重大人員傷亡和基礎(chǔ)設(shè)施、生態(tài)環(huán)境的毀滅性破壞,在活動強度、爆發(fā)規(guī)模、經(jīng)濟損失、人員傷亡等方面均居世界前列[1-3]。分布廣泛、發(fā)生頻繁,突發(fā)性強、預(yù)測預(yù)防難度大,成災(zāi)快、破壞性強、季節(jié)性強、區(qū)域性明顯是中國山洪災(zāi)害的主要特點[4-6]。為防治及減少山洪災(zāi)害的發(fā)生,亟需研究山洪災(zāi)害易發(fā)地區(qū)小流域洪水發(fā)生的規(guī)律及預(yù)測預(yù)警方法,作為山洪災(zāi)害防治措施的依據(jù)[7]。
宣恩縣地處鄂西南,地跨東經(jīng)109°11′~109°55′,北緯29°33′~30°12′之間,東南與湖南省龍山、桑植等縣接壤,行政區(qū)總面積2 730 km2。屬云貴高原延伸部分,地處武陵山和齊躍山的交接部位,縣境東南部、中部和西北邊緣,橫旦著幾條東北至西南走向的大山嶺,形成許多臺地、崗地、小型盆地、平壩、橫狀坡地和山谷、峽等地貌。
宣恩縣屬中亞熱帶季風(fēng)濕潤型山地氣候,隨海拔高程的變化,呈明顯的垂直差異。海拔800.00 m以下的低山帶,四季分明,冬暖夏熱,雨熱同步,年均氣溫15.8 ℃,無霜期294 d,年降水量1 491.3 mm,年日照時數(shù)1 136.2 h;海拔800.00~1 200.00 m的二高山地帶,春遲秋早,濕潤多雨,光溫不足,年均氣溫13.7 ℃,無霜期263 d,年降水量1 635.3 mm,年日照時數(shù)1 212.4 h;海拔1 200.00 m以上的高山地帶,氣候冷涼,冬長夏短,易澇少旱,年均氣溫8.9 ℃,無霜203 d,年降水量1 867 mm,年日照時數(shù)1 519.9 h。
宣恩縣境內(nèi)河網(wǎng)密布,縱橫交錯,共有大小河溪121條,總長度537.65 km,其中流域面積10 km2以上的河流34條,河長在10 km以上的21條。以中部的龍崩山為分山嶺,形成相對獨立的南北兩大水系4條大河,北部貢水、雞籠洞河流歸清江,南部酉水、頭坪河流進沅江,匯入洞庭湖。宣恩縣河流水系情況見圖1。
圖1 宣恩縣河流水系分布圖
本次分析的51個小流域,主要位于宣恩縣李家河鎮(zhèn)(16條)、高羅鎮(zhèn)(11條)及沙道溝鎮(zhèn)(8條)等南部地區(qū),在貢水、高羅河、冉大河、白水河、長潭河等流域均有分布,具體分布情況見圖2。其中流域面積小于5 km2的有1條,介于5~20 km2之間的有13條,介于20~100 km2之間的有24條,介于100~200 km2的有11條,大于200 km2的有2條。最小流域匯流面積為3.85 km2,最大的小流域匯流面積為262.78 km2。小流域坡降介于7.7‰~138.9‰之間。小流域匯流面積匯總見表1。
表1 小流域匯流面積匯總表
圖2 宣恩縣計算小流域分布圖
宣恩縣小流域缺少短歷時實測暴雨資料,結(jié)合山洪災(zāi)害調(diào)查評價項目的要求[8-9],本次小流域設(shè)計暴雨計算采用查算《湖北省暴雨參數(shù)等值線圖集》(2008年,以下簡稱《圖集》)及《湖北省暴雨徑流查算圖表》[10](1985年,以下簡稱《圖表》)方法。小流域暴雨統(tǒng)計參數(shù)(點雨量均值、Cv值)根據(jù)流域中心位置直接查《圖集》;暴雨歷時選取常規(guī)標(biāo)準(zhǔn)歷時(10 min、1、6、24 h)和流域匯流時間2種情況;暴雨頻率為5年一遇(20%)、10年一遇(10%)、20年一遇(5%)、50年一遇(2%)、100年一遇(1%)5種。Kp值及點面換算系數(shù)均由《圖表》查算得出;設(shè)計暴雨統(tǒng)一采用《圖表》中所推薦的5種典型歷時的設(shè)計雨型,為避免漏峰、又兼顧設(shè)計雨型的獲取,計算統(tǒng)一以0.5 h為計算時段。據(jù)此,計算得出各個小流域不同頻率下的暴雨過程。
設(shè)計洪水常用的計算方法有流量推求法與暴雨推求法[11-14],其中后者又包括瞬時單位線法( 降雨徑流相關(guān)圖) 、瞬時單位線法( 初損后損) 、推理公式法、經(jīng)驗公式法和流域模型法[15]等。根據(jù)方法比選,本次設(shè)計洪水計算采用暴雨途徑的瞬時單位線法,分產(chǎn)流計算和匯流計算2個步驟。
(1) 產(chǎn)流計算
產(chǎn)流計算采用降雨徑流相關(guān)圖法。在收集到宣恩縣宣恩(二)站降雨徑流相關(guān)圖后,根據(jù)降雨過程及降雨開始時的前期影響雨量Pa在圖上求出凈雨過程。按如下步驟進行計算。
假定n個時段降雨,第k(k≤n)時段雨量為Pk,降雨開始時的前期影響雨量為Pa(本次取Pa=0.8Wm)。在P-Pa-R圖上Pa對應(yīng)的線上,由1-k時段累積降雨∑Pk查得1-k時段累積凈雨∑Rk,則第k時段凈雨Rk=∑Rk-∑Rk-1,依次可求得凈雨過程。
(2) 匯流計算
匯流計算采用瞬時單位線法。
宣恩縣處于第10水文分區(qū),其瞬時單位線地區(qū)綜合公式為:
m1=0.8F0.3L0.1j-0.06;n=0.69F0.224j0.092
巖溶地區(qū),流域內(nèi)天坑面積超過10%時:
m1=0.76F0.33L0.11j-0.07f-0.12;
式中:m1、n為瞬時單位線法匯流參數(shù);F為流域面積,km2;F天為天坑面積,km2;L為干流河長,km;j為干流平均比降,‰。
由于宣恩縣為有巖溶地質(zhì)和天坑的流域,故各頻率設(shè)計洪水均不考慮非線性改正。
根據(jù)瞬時單位線參數(shù),將設(shè)計頻率的瞬時單位線轉(zhuǎn)換為時段單位線。由設(shè)計頻率的凈雨過程及時段單位線作匯流計算,便得到設(shè)計地表徑流過程。
由穩(wěn)損產(chǎn)生的地下徑流計算公式如下。
當(dāng)t≤T時:
當(dāng)t>T時:
Qt=Qg·e-β(t-T)
其中:
β=0.133F-0.28
式中:Q0為起漲流量,m3/s;Qg為地下徑流洪峰,m3/s;β為退水指數(shù);T為地面徑流過程線底寬,h;fc為穩(wěn)損值,mm;tc為凈雨歷時,h。
將地表徑流過程與地下徑流過程疊加,即得設(shè)計洪水過程線及設(shè)計洪水。
根據(jù)以上計算方法,計算得出宣恩縣51個小流域的設(shè)計洪水成果,為分析成果的合理性,選取整理各個小流域的過程洪量W,繪制與流域面積的關(guān)系圖,見圖3。
圖3 匯流面積F與過程洪量W關(guān)系圖
由圖3可知,過程洪量隨流域面積的增大而增大,隨頻率的減小而增大,并且過程洪量與流域面積呈明顯一次線性相關(guān)關(guān)系(相關(guān)指數(shù)0.98以上),符合地區(qū)總體規(guī)律。
整理宣恩縣51個小流域匯流面積、不同頻率洪峰流量數(shù)據(jù),并進行小流域匯流面積與各頻率洪峰流量進行回歸分析,經(jīng)對比指數(shù)、線性、對數(shù)、多項式、冪等函數(shù)模型后最終采用冪函數(shù)進行擬合分析。分析結(jié)果見表2,小流域匯流面積與各頻率洪峰流量關(guān)系見圖4。
表2 匯流面積F與各頻率洪峰流量Q回歸分析表
圖4 匯流面積F與洪峰流量Q關(guān)系圖
圖5 匯流面積ln F與洪峰流量ln Q關(guān)系圖
為驗證小流域匯流面積F與洪峰流量Q冪函數(shù)關(guān)系是否正確,對51個小流域匯流面積及不同頻率洪峰流量分別取自然對數(shù),在此基礎(chǔ)上進行線性回歸分析及顯著性檢驗。分析檢驗結(jié)果見表2,小流域匯流面積與各頻率洪峰流量關(guān)系見圖5。
根據(jù)分析檢驗結(jié)果,不同頻率洪峰流量lnQ與小流域匯流面積lnF回歸方程顯著性檢驗P值均<0.01,即擬定的冪函數(shù)回歸方程擬合效果很顯著。
根據(jù)以上宣恩縣51條小流域匯流面積、洪峰流量數(shù)據(jù)整理及相關(guān)關(guān)系分析可知:
(1) 隨著匯流面積的增加,各頻率洪峰流量均不斷增大,但增長速率不斷減??;
(2) 隨著匯流面積的增加,頻率越小,洪峰流量增大的速度越快;
(3) 洪峰流量與匯流面積呈冪函數(shù)關(guān)系,回歸方程相關(guān)指數(shù)達(dá)到0.99以上,顯著性檢驗P值<0.01。
本文對宣恩縣51條小流域進行了5種頻率的設(shè)計洪水計算,根據(jù)小流域匯流面積與洪峰流量的相關(guān)分析可知,宣恩縣小流域洪峰流量與匯流面積呈冪函數(shù)關(guān)系,回歸方程相關(guān)指數(shù)達(dá)到0.99以上,顯著性檢驗P值<0.01,對同類地區(qū)小流域設(shè)計洪水成果選取具有重要參考意義。