孤山川洪水特性及遭遇洪水分析

王 鋒，楊彩云

(1.陜西省府谷縣防汛抗旱指揮部辦公室，陜西 府谷 719400；2.陜西省府谷縣皇甫川流域治理指揮部辦公室，陜西 府谷 719400)

孤山川發源于內蒙古準格爾旗川掌，自內蒙古羊市塔入陜西府谷境，經4個鄉、鎮，由府谷鎮高家灣村匯入黃河，全長79 km，流域面積1272 km2。其中府谷縣內長57 km，縣內流域面積1018 km2，河道比降4.1‰。孤山川流域在府谷縣境內較大的支溝有保兔溝、紅石巖溝、安山溝、陽灣川、新城川、木瓜川、后河川、沙川溝等，洪水組成復雜。孤山川河流除縣城設有部分堤防外，總體無設防。為提高縣城防洪能力，需對孤山川干支流匯合處洪水進行分析。

1 暴雨及洪水特性

1.1 暴雨特性

孤山川屬中溫帶半干旱大陸性季風氣候，區內降水量較少，多年平均降水量453 mm，年降水多集中在7、8兩個月內，占全年降水總量的51.6%，降水量年內分配極不均勻，汛期多暴雨。具有春旱多風，秋季涼爽，冬長夏短等氣候特征，陣雨集中于夏季后期，較多且急驟，多為陣雨，常伴有冰雹發生。

1.2 洪水特性

孤山川具有黃河流域洪水特性，洪水跟隨暴雨出現，主要集中在7月、8月，一次洪水過程持續約1天時間，主峰約12小時。具有洪峰流量大、洪水歷時短、陡漲陡落、含沙量大特點。峰高量小、峰型尖瘦，洪峰和枯季流量相差較大。水、沙具有變幅大、比較集中，年內、年際分配不均衡情況。河道大多處于枯期，全年來水、來沙多集中于一次或幾次洪水過程之中。孤山川正常泄洪能力7610 m3/s（50年一遇），水位823.43 m（高石崖水文站），實測最大洪峰流量10300 m3/s（1977年8月2日），最小流量為0，大部分年份均可出現。

從20世紀80年代以來，孤山川來水處于不斷減小而含沙量有所增大的趨勢，1953年～1980年和1981年～2007年兩時段的年最大洪峰流量平均值分別為1760 m3/s和1031 m3/s，1997年之后10年最為明顯，平均徑流量不到多年平均值一半，最大洪峰流量平均值僅為815 m3/s。且常有斷流現象出現，最嚴重的是在1999年，實測年最大洪峰流量僅89 m3/s，為實測系列最小值。

2 水文分析計算

2.1 孤山川洪水資料

孤山川流域設有高石崖水文站，于1953年7月建站，控制流域面積1263 km2，有1954年～2017年64年連續實測徑流、洪水資料。多年平均流量3.49 m3/s，平均徑流量1.101億m3，最大年徑流量2.374億m3（1959年），最小年徑流量0.3024億m3（1965年），相差7.9倍，平均年輸沙量2760萬t。7月～9月徑流、輸沙量分別占全年的70%和95%。實測最大含沙量1300 kg/m3。因此，采用高石崖站水文資料進行洪水分析。

在高石崖水文站64年實測最大洪峰流量系列中，實測最大值為1977年8月2日的10300 m3/s洪峰，水位為823.009 m，這也是孤山川1893年以來調查的最大值。洪水伴隨大暴雨出現，暴雨中心位于內蒙古自治區伊克昭盟的烏審旗、伊金霍洛旗與陜西省神木、府谷兩縣的交界地區，暴雨歷時達8 h～10 h。孤山川流域位于暴雨中心的東南方向，流域平均暴雨量達144 mm。暴雨具有時間長、雨強大、雨區集中、暴雨與河流走向一致等特點，使得洪峰迅速形成。

2.2 洪水頻率分析計算

根據皮爾遜Ⅲ型理論曲線與實測及調查點據之間配合，采用工程中常用的目估適線法確定最佳參數。高石崖水文站年最大洪峰系列的適線成果其均值為1384 m3/s，CV=1.06，重現期為111年，其洪水配線參數CV=1.35，CS=3.5CV，洪水量配線結果見為100年一遇洪峰為9840 m3/s，50年一遇洪峰為7560 m3/s，20年一遇洪峰為4790 m3/s，10年一遇洪峰為2980 m3/s，5年一遇洪峰為1570 m3/s。

2.3 設計防洪標準

根據國家《防洪標準》（GB50201-94）中 2.0.3條規定，“位于山丘區的城市，當城區分布高程相差較大時，應分析不同量級洪水可能淹沒的范圍，并根據淹沒區非農業人口和損失的大小，按規定確定其防洪標準”。府谷縣城共有人口8.2萬人，其中常住人口5.6萬人，是府谷縣社會政治和經濟發展的中心。根據陜西省縣城堤防建設標準，府谷縣城堤防應具有防御50年一遇洪水的能力，相應的孤山川高石崖站洪峰流量應為7560 m3/s。

2.4 洪水遭遇分析

黃河府谷水文站、孤山川高石崖水文站分別至匯合口距離為3182 m和2121 m。府谷縣城段為孤山川與黃河交匯處，在暴雨洪水出現時，也是孤山川洪水與黃河洪水遭遇之時，洪水一般呈正交水流態勢。根據歷史洪水調查數據：1977年，當高石崖水文站出現洪峰10300 m3/s時，黃河府谷站洪水僅為620 m3/s，遭遇后合成流量為10900 m3/s，匯合口壅水高度達1.0 m左右。為合理確定匯合口壅水高度及堤防高度，必須對孤山川與黃河洪水遭遇情況進行分析。根據高石涯水文站和府谷水文站監測資料，現以孤山川為主洪峰和黃河為主洪峰的兩種情況，分別對1964年、1967年、1977年、1992年、2003年5次最大洪水的遭遇情況，進行分析。分析結果見表1。

表1 干支流遭遇分析表

在分析中，采用洪峰倒退辦法，即先推算主洪峰由水文站演進至匯合口的時間，再從匯合口沿另一條河向上，反推至另一條河的最近水文站，求得其相應的流量和出現的時間。將相應流量與主洪峰疊加，即為匯合口兩條河洪水遭遇后的洪峰流量。

分析時，以最大洪峰出現時間前后與洪峰流量值最接近的一次實測斷面最大流速，近似地作為水文站至匯合口的洪峰傳播速度。這種方法只考慮洪峰沿流程的平移，忽略了河槽對洪水速度的阻塞和變形的影響（水文站距入河口距離短，洪水傳播時間短，故其忽略不計）。

從表1可得到以孤山川為主洪峰時的遭遇后洪峰放大系數，結果見表2。由于洪水編號為20030730的洪峰對合成流量貢獻太小，故不予采用，將其它4次洪水放大系數取平均值為1.906。所以取高石崖設計洪峰流量7560 m3/s放大1.9倍作為匯合口的洪峰流量，即14420 m3/s。

根據西安理工大學所作的《陜西省府谷縣城區黃河堤防工程初步設計報告》資料表明，黃河府谷站到孤山川匯合口洪峰放大系數平均值為1.224，取1.2，以黃河來水為主洪峰時，匯合口50年一遇的流量為16470 m3/s。

表2 干支流遭遇后洪峰放大系數

同時假定發生高石崖50年一遇7560 m3/s（頻率為2%）和黃河府谷50年一遇13730 m3/s（頻率為2%）的洪水與發生高石崖10年一遇2900 m3/s和黃河府谷50年一遇13730 m3/s的洪水兩種情況下，干支流匯合口的洪峰流量分別為16470 m3/s和 14420 m3/s。

根據以上分析，最后按較不利的組合，偏于安全的原則，確定干支流匯合口的洪峰流量為16470 m3/s，以此流量計算匯合口壅水高度對孤山川高石崖站的壅水影響，確定匯合處的設計洪水位。

3 匯合處設計洪水位

設計的府谷縣城孤山川河段，高石崖水文站具有1953年～2003年共有51年的實測最大洪水資料。高石崖水文站1977年8月2日流量為10300 m3/s的洪峰時水位為823.009 m，以及實測洪水水位資料，做出水位～流量圖（見圖1），從圖1中可查得設計洪水7560 m3/s的水位為822.60 m。

圖1 高石崖站水位～流量曲線

根據水面比降、糙率，可求出相應于設計流量的水深和水位。所以，設計流量時斷面2+121（高石崖水文站）的水位為822.60 m，即為控制水位。已知設計流量求相應水位是一個試算過程，求出的各斷面正常水深及斷面水力要素見表3，依據2003年《陜西省府谷縣城區黃河堤防工程設計報告》結果，黃河50年一遇和孤山川50年一遇設計洪水匯合后，在河口產生的洪水位為815.52 m，由水力學法推求各斷面的設計壅水洪水位見表3。

表3 河段各斷面水力要素及計算成果

控制斷面設計洪水位是考慮匯合口壅水和兩河洪水遭遇時頂托壅水影響的洪水位，即正常洪水位＋壅水高度。

通過前文干支流遭遇分析得知干支流匯合口最不利組合的洪峰流量為16470 m3/s，以此流量計算由匯合口至設計斷面的壅水水面曲線。采用水力學推求水面曲線，方程如下：

表4 最不利組合下洪峰回水線計算

水面曲線計算的起始斷面為干支流匯合口0#斷面，逐漸向上游推求。從表4可看出，干支流發生都50年一遇的洪水時，匯合口產生的壅水高度為0.228 m，回水長度距匯合口1300 m。干流發生50年一遇的洪水，支流發生10年一遇洪水時，匯合口產生的壅水高度為2.88 m，回水長度距匯合口1180 m，但過流流量較支流50年一遇小，水位低。由此說明，在最不利情況下，匯合口產生的壅水對距匯合口2121多米上的設計段無影響，因此，控制斷面壅水水位即為設計洪水位。

4 結語

孤山川支流眾多，洪水陡漲陡落，泥沙含量大，危害大。孤山川與黃河在府谷縣城交匯，使得此處洪水變化多端，遭遇情況復雜。文章通過探討分析孤山川洪水特性及遭遇的不同情況，計算得出遭遇時50年一遇的流量和水位，為研究黃河陜北段與支流匯合處的洪水遭遇情況提供借鑒。