任青峰，謝輝鈺，張東平，翁 玲

(貴州省赤水市氣象局，貴州 赤水 564700)

赤水河、習(xí)水河流域降水對(duì)水位影響的研究

任青峰，謝輝鈺，張東平，翁 玲

(貴州省赤水市氣象局，貴州 赤水 564700)

該文利用赤水河、習(xí)水河降水、流量、徑流量、水位等資料，建立了赤水河、習(xí)水河暴雨—水位計(jì)算公式，為赤水市小流域暴雨洪澇災(zāi)害防御工作提供決策依據(jù)：赤水9月出現(xiàn)暴雨過程對(duì)赤水河水位影響最大；赤水24 h降水達(dá)到50 mm時(shí)，赤水河水位正常，達(dá)到100 mm時(shí)，赤水河水位將超過警戒水位，其中9月超過危急水位；赤水24 h降水達(dá)到150 mm時(shí)，赤水河水位超過危急水位。習(xí)水河由于缺乏水文數(shù)據(jù)，建立公式偏差較大，預(yù)警準(zhǔn)確度偏低，提升習(xí)水河暴雨—洪水預(yù)警準(zhǔn)確率需要對(duì)習(xí)水河水文資料開展收集，從而修正偏差值。

赤水河；習(xí)水河；暴雨；水位

1 引言

1.1 研究背景

赤水地處于云貴高原向湖南丘陵和四川盆地過渡地帶的斜坡地帶，境內(nèi)地形復(fù)雜多樣，由于暴雨引發(fā)的洪澇等自然災(zāi)害時(shí)有發(fā)生[1]。2014年8月10日晚至11日凌晨, 赤水市普降中到大雨，局地暴雨，但地處習(xí)水河、赤水河流域上游的習(xí)水縣境內(nèi)發(fā)生罕見特大暴雨，全縣有15個(gè)鄉(xiāng)鎮(zhèn)出現(xiàn)暴雨天氣，其中有3個(gè)鄉(xiāng)鎮(zhèn)大暴雨 、5個(gè)鄉(xiāng)鎮(zhèn)特大暴雨，最大降水量為良村鎮(zhèn)334.3 mm，最大小時(shí)雨強(qiáng)達(dá)到84.1 mm，均突破當(dāng)?shù)赜袣庀笥涗浺詠淼臍v史極值。貴州省氣候中心對(duì)這次特大暴雨定性為極端天氣氣候事件。

特大暴雨導(dǎo)致習(xí)水河短時(shí)間迅速形成特大洪水，使地處習(xí)水河下游的赤水市東部5鄉(xiāng)鎮(zhèn)相繼出現(xiàn)超歷史警戒水位，洪峰持續(xù)約2 h。沿岸的石堡鄉(xiāng)、官渡鎮(zhèn)、長期鎮(zhèn)、長沙鎮(zhèn)政府所在地集鎮(zhèn)分別超警戒水位4.7 m、7.4 m、9.9 m、4.8 m，均為有記錄以來的最高水位。由于習(xí)水河上游鄉(xiāng)鎮(zhèn)累計(jì)雨量大、雨勢(shì)猛，加之特殊的地質(zhì)結(jié)構(gòu)和800 m的地形落差，迅猛的洪水夾雜著泥土、砂石導(dǎo)致赤水市東部5鄉(xiāng)鎮(zhèn)發(fā)生了嚴(yán)重的塌方、滑坡、泥石流和山洪災(zāi)害，造成大規(guī)模的房屋垮塌，交通、水利、電力、通訊等基礎(chǔ)設(shè)施嚴(yán)重?fù)p毀。

由于中小河流站網(wǎng)偏少，缺乏必要的應(yīng)急監(jiān)測(cè)手段[2]，為更好的預(yù)防洪澇災(zāi)害，構(gòu)建赤水市主要河流流域防洪減災(zāi)預(yù)警體系，有必要科學(xué)分析赤水河、習(xí)水河降水量和水位關(guān)系,為赤水市的小流域暴雨洪澇災(zāi)害防御工作提供決策依據(jù)。

1.2 自然特征

1.2.1 赤水河 赤水河位于川江南岸，流域處在云南、貴州、四川三省的接壤地帶，發(fā)源于云南省鎮(zhèn)雄縣北部的雨河鎮(zhèn)花果頂梁子，海拔2 000 m，流域集水面積約為20 440 km2，長度約520 km，落差1 588 m，平均比降約3‰。赤水河流域地勢(shì)東、南、西三面高，北部低，大部分為山區(qū)，局部為丘陵和沖積盆地。西與永寧河、南廣河以及橫江的上游分界，南、東隔大婁山脈與烏江流域接壤，東北部緊鄰綦江，北接長江干流及其支流塘河。

1.2.2 習(xí)水河 習(xí)水河，為赤水河下游河口段匯入的東岸大支流。發(fā)源于習(xí)水縣寨壩，東南流至良村下折西北流經(jīng)官渡、長沙入四川境，至合江城附近注入赤水河。流域面積1 600余km2，貴州省境內(nèi)為1 513 km2。包括川、黔兩省的習(xí)水、赤水、合江3縣。河道全長116 km，總落差725 m，平均比降約6‰；貴州境內(nèi)河段長86 km，總落差700 m，平均比降約8‰。流域地勢(shì)東北高而西南低，海拔400～1 200 m，山高坡陡，傍河很少，臺(tái)地多為岸坡上梯田。

2 數(shù)據(jù)和方法

本文的資料來自于赤水河、習(xí)水河多年每月平均降水、流量、徑流量、水位等資料。對(duì)降水與水位關(guān)系進(jìn)行分析，通過參證對(duì)比，獲得降水與水位聯(lián)系公式[3]，并通過歷史暴雨—水位資料檢驗(yàn)，從而建立有一定參考價(jià)值的赤水河、習(xí)水河暴雨—水位計(jì)算模型，為赤水市小流域暴雨洪澇災(zāi)害防御工作提供決策依據(jù)。

3 赤水河降水與水位關(guān)系分析

3.1 降水量、徑流量分析

對(duì)赤水氣象站多年月平均降水與赤水水文站多年月平均徑流量相關(guān)性分析，從圖1可以看出，降水在5—8月最強(qiáng)，而徑流量也在5—8月最大，時(shí)間分布基本一致。

3.2 暴雨、大洪水年內(nèi)發(fā)生機(jī)率分析

從圖2可以看出，赤水河發(fā)生暴雨幾率最高的是7—8月，而赤水河發(fā)生大洪水幾率最大的是6—7月，對(duì)比圖1結(jié)論可以看出在7月赤水出現(xiàn)的強(qiáng)降水將有極大幾率導(dǎo)致赤水河出現(xiàn)洪水。

圖2 赤水河暴雨、大洪水年內(nèi)發(fā)生幾率月分布Fig.2 Monthly distribution of probability of rainstorm and flood in Chishui River

3.3 降水、流量、徑流量關(guān)系分析

通過赤水河多年月平均降雨量與多年月平均徑流量對(duì)應(yīng)分析，應(yīng)用統(tǒng)計(jì)回歸方法建立降水與徑流量關(guān)系圖[4](圖3)，其線性相關(guān)系數(shù)R2=0.934 2，誤差不超過7%。

流量Q，指單位時(shí)間內(nèi)通過某一過水?dāng)嗝娴乃俊搅髁縒，是指時(shí)段內(nèi)通過河流某一斷面的總水量。以所計(jì)算時(shí)段的時(shí)間乘以該時(shí)段內(nèi)的平均流量，就得徑流總量W，即：

W=Q×Δt，直接可用降水與流量開展對(duì)比分析，從而建立降水—流量關(guān)系方程[5]。

由圖4分析，建立的降水—流量線性相關(guān)系數(shù)R2=0.926 0，誤差不超過8%。

圖3 赤水河降水、徑流量線性關(guān)系Fig.3 Linear relation of precipitation and runoff in Chishui River

圖4 赤水河降水、流量線性關(guān)系Fig.4 Linear relation of precipitation and discharge in Chishui River

3.4 水位、流量分析

查閱赤水河水文資料可知，2000年后由于赤水河赤水城區(qū)—復(fù)興段修筑防洪堤，河床變窄，赤水水文站測(cè)流斷面發(fā)生沖淤，流量有所增加。為使水位流量關(guān)系主要反映現(xiàn)狀情況，堤防標(biāo)準(zhǔn)以下的洪水采用2000年以來的資料定線，超過堤防標(biāo)準(zhǔn)的洪水采用建堤前的H-Q關(guān)系。通過查閱赤水河水文資料，建立水位—流量多項(xiàng)式(圖5)，其相關(guān)系數(shù)R2=0.9766，誤差不超過3%。其中水位采用吳淞高程系統(tǒng)，與黃海高程系統(tǒng)的換算關(guān)系為：黃海高程=凍結(jié)高程-187.452 m。

圖5 赤水河水位流量關(guān)系Fig.5 Water level discharge relation of Chishui River

3.5 暴雨、水位分析

由于圖4降水流量關(guān)系是基于多年月平均雨量、流量開展計(jì)算，故本文討論發(fā)生暴雨、洪水時(shí)的情況采用多年月平均雨量、流量為基數(shù)，以赤水河多年月平均降雨量與水位變幅之間建立有一定機(jī)制原理的對(duì)應(yīng)關(guān)系式[6]。以月平均流量通過圖5關(guān)系計(jì)算月平均水位作為基準(zhǔn)水位，根據(jù)歷史暴雨—洪水發(fā)生情況，以5—9月為研究范圍，將降水分為50 mm、100 mm、150 mm、200 mm 4個(gè)量級(jí)，計(jì)算出現(xiàn)降水時(shí)，水位在基準(zhǔn)水位上的增加量，從而確立暴雨—水位關(guān)系：

H=(R-R0/D0)×C1+H0

(1)

其中H為計(jì)算水位，R為設(shè)定降雨量，R0為月平均降雨量，D0為月平均降雨日數(shù)，C1為降雨—水位關(guān)系計(jì)算結(jié)果相關(guān)系數(shù)，H0為月基準(zhǔn)水位。

H0=Q0×C2+C3

(2)

其中Q0為月平均流量，C2、C3為水位—流量關(guān)系計(jì)算結(jié)果相關(guān)系數(shù)。

雖然公式(1)、(2)及式中常數(shù)是通過歷史平均數(shù)據(jù)計(jì)算得出，但是數(shù)據(jù)來源可靠，跨度較長，可以一定程度上反映真實(shí)暴雨—水位關(guān)系，后期通過延長數(shù)據(jù)年限，增加氣候變化變量，可進(jìn)一步完善公式。通過公式(1)計(jì)算出暴雨水位值，結(jié)果如表1：

表1 赤水河暴雨—水位計(jì)算結(jié)果Tab.1 Calculation results of storm water level in Chishui River

表1計(jì)算結(jié)果表明，赤水降水達(dá)到50 mm時(shí)，赤水河水位最高為9月226.77 m，達(dá)不到警戒水位229.0 m；赤水降水達(dá)到100 mm時(shí)，赤水河水位最高為9月231.68 m，超過危急水位231.0 m，且其余月份均超過警戒水位；赤水降水達(dá)到150 mm時(shí)，赤水河水位最高為9月239.31 m，超過危急水位8.31 m，其余月份均超危急水位；赤水降水達(dá)到200 mm時(shí)，赤水河水位最高為9月241.50 m，全部月份都超過有數(shù)據(jù)記錄的歷史極值。

由于上文所述結(jié)果是按照24 h降水計(jì)算，故實(shí)際情況中應(yīng)將降水持續(xù)時(shí)間作為強(qiáng)度修正量計(jì)算。

3.6 結(jié)論檢驗(yàn)

以歷史暴雨—洪水?dāng)?shù)據(jù)作為檢驗(yàn)，查閱歷史赤水城區(qū)洪災(zāi)資料、氣象資料可知[7]，1998年8月赤水降雨量297.1 mm，赤水河市中段洪水水位231.84 m，受淹780戶，房屋2 126間，受災(zāi)人口2 850人，受災(zāi)面積4 650 m2，15個(gè)企業(yè)遭受損失，帶入公式(1)計(jì)算可得如下結(jié)果。

表2 歷史數(shù)據(jù)檢驗(yàn)Tab.2 Historical data check

由于表2計(jì)算時(shí)采用歷史數(shù)據(jù)差值，非實(shí)際過程降雨量，故水位計(jì)算結(jié)果比實(shí)際記錄偏高1.21m，較為接近。

4 習(xí)水河降水與水位關(guān)系分析

根據(jù)前文赤水河降水與水位關(guān)系分析方法，開展習(xí)水河降水與水位關(guān)系分析。

由于習(xí)水河干流原石筍水文站實(shí)測(cè)流量資料系列為1959—1961年，時(shí)間短,代表性差,故直接引用石筍水文站資料不能滿足要求。選用流域面積與習(xí)水河流域相近的桐梓水文站為參證水文站。根據(jù)桐梓水文站的徑流資料，按習(xí)水河流域面積與桐梓站控制面積的倍比，將控制站的徑流系列放大得到習(xí)水河年月徑流系列[8]。

按照上文赤水河暴雨洪水計(jì)算方法，計(jì)算出暴雨水位值，結(jié)果如表3：

表3 習(xí)水河暴雨—水位計(jì)算結(jié)果Tab.3 Calculation results of storm water level in Xishui River

表3計(jì)算結(jié)果表明，習(xí)水降水達(dá)到50 mm時(shí)，習(xí)水河水位最高為5月395.59 m，由于二郎壩水文站歷史水位極差為6.5 m，故習(xí)水河水位變化率為17.32%，各月平均變化率16.53%；習(xí)水降水達(dá)到100 mm時(shí)，習(xí)水河水位變化率最高為7月51.23%，各月平均變化率40.03%；習(xí)水降水達(dá)到150 mm時(shí)，習(xí)水河水位變化率最高為5月75.48%，各月平均變化率67.23%；習(xí)水降水達(dá)到200mm時(shí)，習(xí)水河水位變化率最高為5月95.90%，各月平均變化率87.65%。

5 結(jié)論

①赤水河水文、氣象數(shù)據(jù)較完整，建立公式偏差小，可作為實(shí)際暴雨—洪水預(yù)測(cè)預(yù)報(bào)的參考。

②赤水9月出現(xiàn)暴雨過程對(duì)赤水河水位影響最大；赤水24 h降水達(dá)到50 mm時(shí)，赤水河水位正常；赤水24 h降水達(dá)到100 mm時(shí)，赤水河水位將超過警戒水位，其中9月超過危急水位；赤水24 h降水達(dá)到150 mm時(shí)，赤水河水位超過危急水位。

③習(xí)水河由于缺乏水文數(shù)據(jù)，建立公式偏差較大，預(yù)警準(zhǔn)確度偏低，提升習(xí)水河暴雨—洪水預(yù)警準(zhǔn)確率需要對(duì)習(xí)水河水文資料開展收集，從而修正偏差值。

④由于本文采用定量計(jì)算，故在實(shí)際應(yīng)用過程中需根據(jù)降雨持續(xù)時(shí)間及強(qiáng)度變化情況開展修訂，以接近實(shí)際情況。

[1] 姚正蘭,謝和林,宋芳.遵義市暴雨氣候特征分析[J].貴州氣象,2017,41(2):14-19.

[2] 劉志雨,楊大文,胡健偉.基于動(dòng)態(tài)臨界雨量的中小河流山洪預(yù)警方法及其應(yīng)用[J].北京師范大學(xué)學(xué)報(bào)(自然科學(xué)版),2010,46(3):317-321.

[3] 何艷虎，林凱榮.近期氣候變化對(duì)東江流域水資源的影響研究[J].中國農(nóng)村水利水電，2011,(6):7-10.

[4] 羅喜平,周明飛,汪超.貴州省冬季地表(0cm)溫度預(yù)報(bào)探討[J].貴州氣象,2016,40(4):1-5.

[5] 李勇,劉壽東,李登文，等.氣象水文藕合預(yù)報(bào)模式預(yù)警洪澇一以湄潭縣河流為例[J].氣象科學(xué)，2007,27(12):115-120.

[6] 周筱蘭,張禮平,王仁喬,等.應(yīng)用最優(yōu)化訂正法制作長江上游面雨量預(yù)報(bào)[J].氣象,2003,29(3):31-33.

[7] 朱定真,陳玉石,王清樓等.汛期長江中下游降雨量與下游水位增幅的預(yù)測(cè)[J].氣象科學(xué),2002,22(3):352-356.

[8] 郭良,唐學(xué)哲,孔凡哲.基于分布式水文模型的山洪災(zāi)害預(yù)警預(yù)報(bào)系統(tǒng)研究及其應(yīng)用[J].中國水利,2007,(14):23.

StudyontheinfluenceoftheChishuiRiver,XishuiRiverBasinPrecipitationonWaterLevel

REN Qingfeng，XIE Huiyu，ZHANG Dongping，WENG Ling

(Chishui Meteorological Bureau of Guizhou Province， Chishui 564700， China)

This paper uses the Chishui River, Xishui river flow, rainfall, runoff, water level data, establishes the calculation formula of the Chishui River, Xishui River storm - water level, provides the decision-making basis for the small watershed of Chishui city flood disaster prevention work: rainstorm occurred in Chishui in September, the greatest impact on the Chishui River; Chishui 24h precipitation at 50mm the Chishui river was normal; Chishui 24h precipitation reaches 100mm, the Chishui river water level will exceed the warning level, which in September exceeded the critical water level; Chishui 24h precipitation reaches 150mm, the Chishui river water level exceeds the critical level. Xishui river due to the lack of hydrological data, establish the formula of large deviations, early warning accuracy is low, Xishui River storm flood warning and improve accuracy of Xishui river hydrological data need to be collected to carry out, to correct the deviation.

Chishui River; Xishui River; rainstorm; water level

1003-6598(2017)05-0058-04

2017-02-21

任青峰(1989-)，男，助工，主要從事天氣預(yù)報(bào)服務(wù)工作，Email:jesusr@163.com。

貴州省氣象局青年基金項(xiàng)目(QN[2015]12號(hào))資助：赤水河、習(xí)水河上游降水對(duì)下游水位影響的預(yù)警機(jī)制。

P426.61+3

B