于春平

（新疆兵團(tuán)勘測(cè)設(shè)計(jì)院（集團(tuán)）有限責(zé)任公司）

0 引言

阿克蘇柯克亞河出山口至末端之間洪災(zāi)泛濫面積為298.50 km2，區(qū)域夏季暴雨匯入柯克亞河和依來克溝，兩洪溝在地區(qū)實(shí)驗(yàn)林場(chǎng)匯合為卡爾斯亞溝，由西北向東南，流經(jīng)阿克蘇市東城區(qū)，跨過多浪渠，退入新大河。卡爾斯亞溝多次給阿克蘇市東城區(qū)帶來洪水災(zāi)害。

1 流域概況

柯克亞河發(fā)源于溫宿縣境內(nèi)天山托木爾峰的南坡，是一條山區(qū)冰川融水和降水補(bǔ)給為主的河流；依來克溝發(fā)源于天山南坡，位于溫宿縣柯克亞河和臺(tái)蘭河之間，由三條較大洪水溝組成，其中上游玻璃廠洪水溝的洪水在臺(tái)蘭河老龍口經(jīng)泄洪渡槽全部退入臺(tái)蘭河，其余兩條洪水溝出山口后由北向南穿越戈壁礫石帶，在實(shí)驗(yàn)林場(chǎng)與柯克亞河匯合成為卡爾斯亞溝；卡爾斯亞溝流經(jīng)阿克蘇市東城區(qū)，最后匯入新大河。

2 洪水類型及特點(diǎn)

柯克亞河洪水的成因可分為暴雨型洪水、冰雪消融型洪水、暴雨與冰雪消融混合型洪水，其中以暴雨型洪水最為突出。據(jù)分析，柯克亞河洪峰流量在30 m3/s 以上的洪水均為暴雨型洪水，峰值在50 m3/s 以上的洪水，其主峰段均由暴雨形成；依來克溝洪水為暴雨型洪水。

3 歷史洪水

3.1 歷史洪水訪問

2002 年對(duì)柯克亞河歷史洪水進(jìn)行訪問，訪問對(duì)象為紅旗坡水管站2 名退休職工，年齡為62 歲和68 歲（1958 年到2002年一直在柯克亞龍口從事水管工作），兩人均表示1971、1984、1989、1998、2002年柯克亞河河水較大。

3.2 洪水調(diào)查情況

自1984年起，水文工作者多次對(duì)柯克亞河洪水進(jìn)行調(diào)查，調(diào)查情況見表1。

3.3 重現(xiàn)期

根據(jù)柯克亞河歷史洪水訪問和調(diào)查情況，確定2017 年洪水是1958年來最大的一場(chǎng)洪水，洪水重現(xiàn)期為60 a一遇。

表1 洪水調(diào)查表

3.4 設(shè)計(jì)洪峰流量

柯克亞河柯克亞水文站只有1 a 不完整實(shí)測(cè)資料，且洪水期缺測(cè)，因此水文分析計(jì)算選取鄰近流域臺(tái)蘭河臺(tái)蘭水文站為參證站。

3.4.1 參證站設(shè)計(jì)洪峰流量

對(duì)臺(tái)蘭水文站1957－2013 年實(shí)測(cè)洪水系列進(jìn)行頻率分析，得到臺(tái)蘭水文站設(shè)計(jì)洪峰流量，成果見表2。

表2 臺(tái)蘭水文站洪水分析成果表 單位：m3/s

3.4.2 柯克亞河與依來克溝設(shè)計(jì)洪峰流量

柯克亞河與依來克溝無實(shí)測(cè)洪水資料，采用無資料地區(qū)洪水計(jì)算方法得到各河溝設(shè)計(jì)洪水。

3.4.2.1 洪峰模數(shù)法

選用鄰近流域臺(tái)蘭水文站為參證站，利用洪峰模數(shù)法計(jì)算各河溝設(shè)計(jì)洪峰流量，成果見表3。

式中：、Q設(shè)計(jì)、F設(shè)計(jì)為設(shè)計(jì)站洪峰流量與流域面積；Q參正、F參正為參證站洪峰流量與流域面積（km2）；np為參數(shù)，新疆無資料地區(qū)小流域一般取2/3。

表3 洪峰模數(shù)法計(jì)算各河溝設(shè)計(jì)洪水成果表 （單位：m3/s）

3.4.2.2 地區(qū)綜合經(jīng)驗(yàn)公式法

利用南疆巴州至克州有實(shí)測(cè)洪水資料的8 個(gè)中小河流水文站建立暴雨洪峰流量經(jīng)驗(yàn)公式（1）計(jì)算各河溝多年平均年最大洪峰流量，移用臺(tái)蘭水文站設(shè)計(jì)參數(shù)，得到各河溝設(shè)計(jì)洪峰流量，成果見表4。

式中：Qm為多年平均年最大洪峰流量（m3/s）；F 為控制斷面流域面積（km2）；L為河流長(zhǎng)度（km）；J為河道縱坡（‰）；F為流域形狀系數(shù)，f=F/L2；Pm為流域多年平均最大24 h 降雨量（mm）。

3.4.2.3 地區(qū)洪峰流量綜合頻率曲線法

選取與柯克亞河流域臨近的臺(tái)蘭、卡木魯克、卡拉蘇、黑孜水文站作為參證站，用地區(qū)洪峰流量綜合頻率曲線法，推求調(diào)查斷面不同頻率洪峰流量，公式如下：

表4 地區(qū)綜合經(jīng)驗(yàn)公式法計(jì)算各河溝設(shè)計(jì)洪峰流量成果表

式中：QP為設(shè)計(jì)洪峰流量（m3/s）；KP為設(shè)計(jì)洪峰流量模比系數(shù)；Qd為調(diào)查歷史洪水洪峰流量（m3/s）；Kd為調(diào)查洪水洪峰流量模比系數(shù)。

根據(jù)2017年9月洪水調(diào)查成果，實(shí)驗(yàn)林場(chǎng)大橋斷面洪峰流量241 m3/s，為60 a一遇洪水，模比系數(shù)Kd=3.95。

表5 地區(qū)洪峰流量綜合頻率曲線法計(jì)算實(shí)驗(yàn)林場(chǎng)大橋斷面設(shè)計(jì)洪峰流量成果表

3.4.2.4 推理公式法

3.4.2.4.1 最大24 h降水量

臺(tái)蘭水文站是距柯克亞河最近的水文氣象站，有實(shí)測(cè)的最大1 d降水量資料，由臺(tái)蘭水文站最大1 d降水量推算柯克亞河和依來克溝最大24 h降水量，計(jì)算見公式4～6。

式中：p1r為臺(tái)蘭水文站最大1 d 降雨量（mm）；K1為項(xiàng)目區(qū)內(nèi)河溝與臺(tái)蘭水文站降雨量折算系數(shù)，臺(tái)蘭水文站高程為1 648 m，柯克亞河的山前暴雨區(qū)域平均高程在3 000 左右，通過分析柯克亞河K1取1.10，依來克溝取1；h24為最大24 h點(diǎn)降雨量（mm）；P1r為最大1 d降雨量（mm）；K2為最大1 d降雨量與最大24 h 降雨量折算系數(shù)，取K2=1.13；H24為最大24 h面降雨量（mm）；K3為最大24 h點(diǎn)、面降雨量折算系數(shù)，取K3=0.99。

對(duì)臺(tái)蘭水文站1965-2007年最大1 d降水資料進(jìn)行頻率分析，得到臺(tái)蘭水文站設(shè)計(jì)最大1 d降水成果，見表6。由表6和公式4～6計(jì)算計(jì)算得到各流域不同頻率設(shè)計(jì)降水成果，見表7。

表6 臺(tái)蘭水文站最大1 d降水量計(jì)算成果表 單位：mm

表7 各河溝最大24 h降水計(jì)算成果表 單位：mm

續(xù)表7

3.4.2.4.2 設(shè)計(jì)洪峰流量

設(shè)計(jì)洪峰流量計(jì)算公式如下：

式中：Qm為設(shè)計(jì)洪峰流量（m3/s）；ψ為洪峰徑流系；SP為設(shè)計(jì)頻率的雨力（mm/h）；τ為匯流時(shí)間（h）；tc為產(chǎn)流時(shí)間（h）；m為匯流參數(shù)；F為流域面積（km2）；L為河長(zhǎng)（km）；J為河道比降（用小數(shù)表示）；θ為流域特征參數(shù)；μ為產(chǎn)流參數(shù)；R 為損失系數(shù)，取1.08；r1為損失指數(shù)，取0.75；n為暴雨衰減指數(shù)，根據(jù)新疆《可能最大暴雨圖集》，確定n取0.75。

用上述公式，采用迭代試算法求設(shè)計(jì)洪峰流量，見表8。

表8 各河溝設(shè)計(jì)洪峰流量成果表

卡爾斯亞溝由柯克亞河與依來克溝在實(shí)驗(yàn)林場(chǎng)匯合而成，實(shí)驗(yàn)林場(chǎng)大橋斷面洪峰流量由柯克亞河與依來克溝出山口洪峰流量衰減疊加而來。柯克亞河、依來克1溝、依來克21溝、依來克22 溝出山口距離實(shí)驗(yàn)林場(chǎng)大橋斷面分別為36.70 km、34.90，37.10，39.40 km，移用臺(tái)蘭河洪峰衰減系數(shù)1.42%/km，得到實(shí)驗(yàn)林場(chǎng)大橋斷面設(shè)計(jì)洪峰流量見表9。

表9 實(shí)驗(yàn)林場(chǎng)大橋斷面設(shè)計(jì)洪峰流量成果表

4 結(jié)語(yǔ)

4.1 成果選定及合理性分析

文章選取無資料地區(qū)中小河流洪水計(jì)算常用的4 種方法分析柯克亞河與依來克溝設(shè)計(jì)洪水，經(jīng)對(duì)比分析推理公式法與地區(qū)洪峰流量綜合頻率曲線法成果接近。2017年8月25日實(shí)驗(yàn)林場(chǎng)大橋斷面調(diào)查洪峰流量241 m3/s，洪水重現(xiàn)期為60 a 一遇，與推理公式法推算得到的60 a一遇洪峰流量為251 m3/s 接近，說明2017年9月調(diào)查洪峰流量和重現(xiàn)期的確定較為合理。選取與柯克亞河流域臨近的臺(tái)蘭、卡木魯克、卡拉蘇、黑孜水文站作為參證站，用地區(qū)洪峰流量綜合頻率曲線法推求的設(shè)計(jì)洪峰流量成果可靠。因此推薦選用地區(qū)洪峰流量綜合頻率曲線法推求的設(shè)計(jì)洪峰流量成果。

4.2 卡爾斯亞溝代表斷面設(shè)計(jì)洪水

實(shí)驗(yàn)林場(chǎng)大橋至喬格塔電站之間工程斷面設(shè)計(jì)洪峰流量采用2017 年8 月25 日調(diào)查洪水衰減率推算，實(shí)驗(yàn)林場(chǎng)大橋斷面洪峰流量241 m3/s，喬格塔電站洪峰流量193 m3/s，兩斷面間距離16.99 km，洪水衰減系數(shù)1.18%/km。代表斷面設(shè)計(jì)洪峰流量見表10。

表10 卡爾斯亞溝代表斷面設(shè)計(jì)洪峰流量計(jì)算表