三道河子荒溝水文站產匯流方案成果分析

徐敏亮

(黑龍江省牡丹江水文局，黑龍江 牡丹江 157000)

?

三道河子荒溝水文站產匯流方案成果分析

徐敏亮

(黑龍江省牡丹江水文局，黑龍江 牡丹江 157000)

摘要：牡丹江市位于長白山脈以北、張廣才嶺以東，地貌特征以山區(qū)、丘陵為主，河流眾多，雨量大，暴雨頻繁，水資源豐富，水文狀況十分復雜。為做好測區(qū)內洪水預報工作，牡丹江水文局編制了測區(qū)內的實用洪水預報方案，文章以三道河子荒溝站的產匯流預報方案為例，對三道河子自然地理特征，水文氣象條件，荒溝水文站預報斷面情況進行了詳細闡述，對三道河子荒溝水文站以上流域產、匯流方案進行了分析，建立了三道河子荒溝站的預報方案模型，進行了方案評定，達到了為該站的短期洪水預報提供科學可靠依據(jù)的目的。

關鍵詞：產匯流；三道河子；徑流深；單位線；成果評定

0前言

牡丹江市境內河流眾多，水文站網布設密度不足，絕大多數(shù)河流沒有水文站，2015年中小河流水文監(jiān)測系統(tǒng)建立了18處水文站，一定程度上彌補了中小河流水文監(jiān)測的空白，但新建站沒有足夠的資料系列進行洪水預報。此外，2 000 km2以下流域面積的山區(qū)河流，河流比降大，匯流時間短，如果沒有成型的預報方案或模型，很難進行洪水預報。為此，我們選取了三道河子荒溝站作為代表，編制產匯流預報方案，以期能為其它200多條河流的洪水預報提供依據(jù)。

2流域概況

2.1自然地理特征

三道河子是牡丹江左岸支流，發(fā)源于張廣才嶺黑龍江省海林市東南部,河長79.5 km，流域面積1 455 km2，坡降4.55‰。自上而下主要支流有芭藁砬子河、腰營河、干明益河、新房子河等。

三道河子流域地處山區(qū)，張廣才嶺主脈綿亙于西部，山頂多為半渾圓狀，群山起伏，坡陡流急。在荒溝站斷面下游12 km處匯入牡丹江(現(xiàn)蓮花水庫)。本區(qū)域為侵蝕剖蝕低山地形，土層較松，透水性良好，樹木茂密，尤其是流域上游有相當一部分較厚的枯枝落葉層，土壤以暗棕壤為主，其余為草甸土、沼澤土、水稻土等。本站測驗河段順直部分約700 m，左岸漫灘到山腳，寬約50 m，右岸有土路。

流域內地形多山起伏、坡陡流急、土層較、透水性良好、樹木茂密，尤其流域上游有相當一部份較厚的枯枝落葉層，對本流域長期無雨或少雨的情況下，保持一定數(shù)量的地下徑流補給，是該流域基本特點之一。

2.2水文氣象特征

本站屬東亞季風區(qū)，即中溫帶大陸性季風氣候，具有冬長嚴寒干燥，夏季降水集中，氣候溫熱、濕潤；春季多大風，降水少，易發(fā)生干旱；秋季氣溫下降快常有霜凍發(fā)生。年平均氣溫為3.2℃，1月平均氣溫-18.1℃，7月平均氣溫21.8℃。無霜期85-130天。由于受季風影響，本站以偏西風為主，年平均風速2.5m/s。

多年平均降水量568.2 mm，降水的年際變化較大，歷年最大降水量為1994年804.6 mm，歷年最低(1975年)降水量僅292.8 mm。

本站降水主要集中在汛期6—9月，多年平均6—9月降水量為412.8 mm，約占年降水量的72.7%。歷年最大6—9月降水量624.9 mm(1960年)，歷年最小6—9月降水量220.9 mm(1975年)，

歷年最大月降水量291.8 mm(1960年8月)，歷年最大日降水量114.9 mm(1960年8月23日)，歷年最大1 h降水量45.3 mm(1977 年7月29日02時)。

本站洪水為典型的雨洪徑流，年最大洪峰集中出現(xiàn)在汛期7—10月，尤以7、8月份最為集中，其出現(xiàn)最大洪峰的概率為75%。洪水來源是三道河子干流。洪水過程大多是單式洪峰，洪峰停留時間一般在2 h左右，漲落水歷時都比較短，峰形較瘦,洪水漲落過程歷時為1～2 d。

3預報站基本情況

荒溝水文站位于黑龍江省海林市柴河林業(yè)局板橋林場。集水面積1 327 km2，是牡丹江左岸支流三道河子控制站，蓮花水電站的入庫站之一。

該站測驗河段位于三道河子下游唯一的峽谷處，此峽谷長約1 km?；舅邤嗝嫖挥跍y驗河段中間，流速儀兼浮標測流斷面與基本水尺斷面重合，最大洪水水面寬不超過150m，上、下浮標斷面間距200m，上、下浮標斷面分別位于流速儀測流斷面上下各100m處，流速儀測流斷面上、下300m處分別為上比降斷面和下比降斷面。

本站為復式河槽，斷面垂直于水流方向，左岸為灘地，蔓延到高山，右岸為公路。河底為大卵石及塊石組成，間有粗砂卵石；左岸表層為黑黏土，其下為粗砂和卵石，右岸靠岸邊表層為黑砂土，其下有零星大塊石，高壩以上為黃砂土。

水位流量關系變化的基本特性：本站歷年水位流量關系比較穩(wěn)定。由于下游有彎道，形成測站控制，且河道比較穩(wěn)定，即使是大洪水也無嚴重的河道變形，因此水位流量關系曲線分布在一個帶狀內，其變幅不到5%。水位91.00m以上形成單一的水位—流量關系曲線。

4應用資料

1953年設荒溝雨量站，1958年6月改為水文站，觀測項目比較齊全，是牡丹江地區(qū)建站比較早的一個國家基本水文站。本次產匯流分析采用的降水、蒸發(fā)、水位、流量等資料系列為1960年—2004年，雨量資料共計7220站年，流量資料共計31站次。本次產、匯流分析共采用了31年流量資料中的48次孤立洪峰及峰型較好的復式洪峰參加計算。

5產流部分

從荒溝站以上流域特性看，除了個別年份(前期比較干旱)有超滲產流現(xiàn)象外，大部分是屬于蓄滿產流，我們進行產流分析時，確定該流域為蓄滿產流。

降雨徑流計算采用下列公式：

Wm′=(1+b)Wm

當P+a

當P+a≥Wm′，R=P+W0-Wm

(1)

式中：Wm′為蓄水量中最大值；b為反映流域內蓄水容量不均勻系數(shù)；W0為流域前期蓄水量；Wm為流域平均最大蓄水量。

5.1產流參數(shù)的選用

荒溝站與長汀子水文站同屬海林市的山區(qū)，河流特性、流域植被覆蓋率、土壤含水能力、蒸散發(fā)等情況基本類似。依據(jù)長汀子水文站Im=P+Pa-R-Es計算的計算結果，選用lm為80 mm。根據(jù)流域特點，選用b為0.30。

5.2流域蓄水量Pa的計算

流域內氣包帶土壤含水量的變化是影響降雨形成徑流過程的一個主要因素，但又沒有全流域土壤含水量的實測資料，只能用間接的方法求得流域前期蓄水量Pa：

有雨時：Pa，t+1=K(Pt+Pa，t-Rt)

無雨時：Pa，t+1=K Pa，t

令5月1日Pa=Im/2，當Pa，t+1≥Im時，則以Pa，t+1=Im為控制。

5.3實測徑流深的計算

5.3.1基流的分割

本方案用綜合退水曲線，將前次洪峰過程和本次洪峰退至基流，然后平割基流。根據(jù)當年汛期最小流量，其范圍在5～15m3/s。前期比較干旱的年份基流采用53/s，一般水年份基流采用10m3/s。

5.3.2綜合退水曲線的采用

本站綜合退水曲線采用的是退水曲線的外包線。參照相應的降水量產生的徑流，進行前峰退水及峰后雨造成的徑流分割，退至基流為止。

5.3.3實測徑流深R的計算

R實=3.6×∑QIΔt/A

(2)

式中：∑QI為時段徑流量之和，m3/s；A為流域面積，km2；Δt為時段，h。

6匯流部分

由于暴雨洪水資料較少，我們從歷年實測資料中選出1978、1981、1997、2000、2002年共5次洪水，進行單位線分析。由于流域集水面積小，下墊面也較為一致，所以不進行分塊計算，又考慮應用時主要是預報洪峰值及主要過程，所以沒有進行水源劃分，把地面徑流和壤中流合并起來進行匯流計算。

6.1流域平均時段凈雨量計算

1960-1975年流域內只有出口斷面荒溝一處雨量站，代表性較差，計算流域平均雨量時，引用了流域外雨量站，采用算術平均法求值。1976-2004年流域內有5處雨量站，用面積加權法求出時段雨量。采用造洪雨量開始時的Pa，用累積時段雨量法計算凈雨過程。

6.2地面徑流量計算

用綜合退水曲線結合時段雨量進行前峰退水及峰后雨造成的徑流分割，退至基流后，把地面及壤中流合并計算出地面徑流量。

6.3單位線計算

單位線時段Δt選取采用6 h，能滿足推算地面徑流過程的要求，根據(jù)時段凈雨及假定單位線，求出地面徑流過程，并與實測徑流過程擬合，著重于峰值及其附近過程的吻合，此假定的單位線即為所求的單位線。

7誤差評定

根據(jù)《水文情報預報規(guī)范》GB/T 22482—2008(6.5.4)中的規(guī)定，徑流深預報以實測值得20%作為許可誤差，當該值>20 mm時，取20 mm；當<3 mm時，取3 mm。本次計算共31個點據(jù)22次合格，方案的合格率為22/31×100%=71.0%，為乙級方案。

8結語

本方案在應用時，因本流域較小，暴雨中心不明顯，主要是根據(jù)雨型及其時空分布與單位線原型進行對照，近似于哪次徑雨過程對應的單位線就利用其單位線進行作業(yè)預報。上述預報模型還存在一定誤差，尚需作進一步研究和完善。

參考文獻：

[1]黑龍江省水文局.黑龍江省實用洪水預報方案[R].哈爾濱：黑龍江省水文局，2009.

[2]莊一鸰，林三益.水文預報[M].北京：水利電力出版社，1986.

Analysis of Runoff Scheme of Huanggou Hydrologic Station in Sandaohezi

XU Min-liang

(Mudanjiang Hydrologic Bureau Heilongjiang Province，Mudanjiang 157000，China)

Abstract：Mudanjiang City is located in the north to the Changbai Mountain，the east to the Zhangguangcai Ridge，its topographic features take the mountainous areas and hilly landform as the main，there are many rivers，large rainfall，frequent rainstorm and rich water resources，and the hydrologic regime is very complex.In order to ensure the flood forecast job well done in the survey area，Mudanjiang Hydrologic Bureau works out a practical scheme of flood forecast within the survey area，taking the runoff yield and concentration scheme of Huanggou Station in Sandaohezi as the example，the paper describes in details the natural geographic characteristics，hydrological and meteorologic conditions and the section forecast of Huanggou Hydrologic Station，by analyzing the scheme of runoff yield and runoff concentration at the upstream of the Huanggou Hydrologic Station in Sandaohezi，to build a forecast scheme model for the Huanggou Station in Sandaohezi to evaluate the scheme，realizing to supply scientific accordance of short-time flood forecast for the station.

Key words：runoff；Saodaohezi；depth of runoff；unit line；evaluation of results

文章編號：1007-7596(2016)04-0096-03

[收稿日期]2016-03-11

[作者簡介]徐敏亮(1978-)，男，黑龍江望奎人，工程師，從事水文情報預報，洪水預報，水資源研究工作。

中圖分類號：TV121+.7

文獻標識碼：B