西南巖溶地區徑流計算研究
——以云南某水庫計算為例

姜姍姍，韓冬梅

（中水北方勘測設計研究有限責任公司，天津 300222）

擬建的某水庫位于珠江流域西江水系南盤江右岸一級支流清水江上游南丘河，云南省文山州境內。清水江發源于硯山縣的回龍水庫，依次流經硯山、丘北、廣南3 個縣，于丘北縣的壩達村匯入南盤江，全長229 km，集水面積5 117 km2，總落差790 m，河道平均坡降3.4‰。清水江南支為南丘河，北支為清水河，南丘河、清水河匯合后稱清水江。南丘河發源于硯山縣回龍水庫，河流呈南北流向，流經硯山城區、者臘、干河、維摩等鄉鎮，進入丘北縣天星鄉，在法白堂上村與清水河交匯后進入清水江，全長110 km，集水面積1 746 km2，落差140 m，河道平均比降4.1‰。境內多為巖溶地貌，地表水展布格局受巖性和地質構造的控制，徑流存在異向排泄的問題，分水嶺難以確切劃分，如長嶺街、岔路口等地，大面積溶斗網格排列十分典型。某水庫壩址以上流域巖溶發育，存在地表水分水嶺與地下水分水嶺不一致的情況，因此壩址流域面積的計算要結合水文地質圖扣除地表水分水嶺與地下水分水嶺不一致的巖溶區面積。經量算，壩址地表集水面積為1 539 km2，扣除巖溶區面積后，壩址徑流面積取1 171 km2。

壩址斷面附近設有清水河水文專用站，有2017—2019 年3 a 實測資料，壩址下游清水江干流上革雷站有1970—2005 年長系列水文觀測資料，革雷站下游清水江水文站有2001—2019 年觀測資料。因此，以革雷站為設計參證站，同時參考清水河水文專用站實測成果。

1 壩址年徑流分析計算

某水庫壩址處清水河水文專用站僅有3 a 的實測流量資料，本次分別采用徑流深等值線圖法、水文比擬法、降水徑流關系法3種方法推算壩址徑流。

1.1 徑流深等值線圖法

根據文山州1956—2013 年徑流深等值線圖，通過讀圖確定壩址多年平均徑流深。先分析1956—2013、1970—2019 年時段多年平均降水量情況，再分析是否需要對讀圖成果進行修正，利用壩址以上集水面積乘以徑流深計算壩址以上流域1970—2019 年多年平均徑流量。經分析計算，壩址多年平均徑流深為330 mm，由此計算壩址以上流域1970—2019年多年平均徑流量為3.86億m3。

1.2 水文比擬法

某水庫所在的清水江流域有北橋、革雷、清水江3個國家基本水文站，但只有革雷站資料系列較長，以革雷站為設計依據站，采用下式計算各壩址徑流量：

式中：W革雷,i，W革雷,i分別為壩址，革雷站第i年徑流量（億m3）；F壩址，F革雷分別為壩址，革雷站徑流計算面積（km2）；P壩址，i，P革雷，i分別為壩址以上流域，革雷站以上流域第i年降水量（mm），采用泰森多邊形法計算。

其中，壩址以上流域面平均雨量計算采用天星、樹皮、清水河、五家寨、大新寨和阿野共6 個雨量站資料，革雷站以上流域面平均雨量計算采用膩腳、大戛勒、樹皮等共11個雨量站資料。

根據式（1）計算，壩址1970—2019 年多年平均徑流量為3.77億m3、徑流深為322 mm。

1.3 降水徑流關系法

根據革雷站降水徑流關系，按下式計算壩址逐年徑流量：

式中：W壩址采用，i，W壩址查圖，i分別為壩址第i年采用的徑流量（億m3），壩址第i年根據壩址降水量由革雷站以上流域P-R關系查得的徑流量（億m3）；W革雷站實際，i，W革雷站查圖，i分別為革雷站第i年實際天然徑流量（億m3），革雷站第i年根據革雷站降水量由革雷站以上流域P-R關系查得的天然徑流量（億m3）。

按此方法，先根據革雷站以上流域P-R 關系確定壩址逐年徑流查算值，然后乘以革雷站同一年徑流實際值與徑流查圖值的比例，以對壩址徑流進行離散化處理。經計算，壩址1970—2019 年多年平均徑流量為3.63億m3、徑流深為311 mm。

從計算成果來看，徑流深等值線圖法、水文比擬法、降水徑流關系法3 種方法推算壩址年徑流量呈依次減少趨勢，但總體較為接近，相對誤差在7%以內，詳見表1。

表1 不同方法壩址1970—2019年多年平均徑流量對比 億m3

2 壩址徑流年內分配

通過對比清水河水文專用站徑流月分配與革雷站徑流月分配過程發現，兩站點徑流年內分配差異較大，因此直接采用革雷站徑流月分配成果不合理。

造成壩址和革雷站徑流月分配過程不一致的原因可能是地下暗河的調節作用。革雷站以上的6條暗河中，出露點在壩址以上的有4條，流量合計值為0.87 m3/s，占壩址多年平均流量的7.7%；出露點在壩址—革雷站之間的有2 條，革雷站以上暗河流量合計值為8.1 m3/s，占革雷站多年平均天然流量的24.2%。與革雷站相比而言，壩址徑流中由暗河形成的地下徑流占比明顯偏低，地表徑流占比明顯偏高。由暗河形成的地下徑流由于具有多年調節作用，相對穩定，年內變化不大，地表徑流由于主要受當年降水影響，年內變化較大。由于大泉、匯流富集帶無明顯地帶性差異，總體來看，壩址徑流年內分配應比革雷站更加不均勻，汛期6—10 月占比更高，非汛期11—5月占比更低。因此，采用2種方法分析壩址徑流月分配情況。

2.1 根據革雷站月分配過程進行修正

在壩址1970—2019 年逐年年徑流量的基礎上，按革雷站逐年各月徑流占比計算壩址同一年的逐月徑流，作為初步成果。在此基礎上，對清水河水文專用站徑流進行還原計算，得到壩址天然徑流直接計算成果。基于革雷站逐月天然徑流，采用水文比擬法計算壩址月徑流，得到壩址徑流間接計算成果。枯期（11—5 月）直接計算值占間接計算值的比例為0.43。因此，在壩址徑流年內分配初步計算成果的基礎上，以枯期徑流乘以0.43 的系數得到枯期修正后的徑流，然后將枯期原徑流與修正后的徑流差值分配到汛期，得到壩址徑流年內分配的采用成果。

2.2 根據水文專用站插補延長計算

雖然清水河水文專用站觀測年限較短，但作為壩址專用水文站其觀測成果對確定壩址徑流有非常重要作用。鑒于數據較少的情況，本次對清水河水文專用站和下游清水江水文站2017—2019 年逐日徑流進行相關分析，根據相關關系插補延長清水河水文專用站徑流，兩者相關關系如圖1所示。從圖1可以看出，兩者相關性較好，以此插補延長清水河水文專用站徑流作為壩址徑流，插補后計算壩址1970—2019年多年平均徑流量為3.33億m3。

圖1 清水河和清水江水文站徑流相關關系

通過上述2 種方法計算的壩址徑流年內分配，如圖2所示。

圖2 兩站相關法和修正法計算的壩址徑流年內分配過程

從圖2可以看出，2種方法計算的年內分配成果較為接近，說明對參證站月分配過程進行修正是合理的。

3 結論

本次計算清水江流域某水庫壩址徑流時，充分考慮了巖溶地區水文地質特性，同時確定壩址集水面積時充分考慮了灰巖區漏水情況及流域地表水分水嶺與地下水分水嶺不一致情況。在此基礎上，采用徑流深等值線圖法、水文比擬法以及降雨徑流關系法3 種方法計算了壩址年徑流，同時對壩址徑流年內分配計算方法進行了探索。

徑流深等值線圖法綜合考慮了降水量、高程、區域平衡等因素，是規劃階段常用的方法，工程設計階段還需在此基礎上結合擬建工程的具體情況做進一步分析復核；水文比擬法綜合考慮了壩址降水量、革雷站降水量及徑流量等因素，比較充分地反映了壩址與設計依據站的降水及下墊面條件；降水徑流關系法在水文比擬法基礎上又考慮了徑流系數等因素，充分地反映了壩址與設計依據站降水及下墊面條件；徑流相關法雖然出現年份較短，但是水文專用站是直接反映壩址徑流的站點，對其進行延長分析，成果也較為可靠。實際進行工程設計時，根據具體情況合理選擇計算方法和成果。