基于SWAT的灃河流域水文特征研究

楊靜媛，郗 林

灃河屬黃河流域渭河水系，是渭河右岸的一級支流，長安八水之一，也是西安市第三大河流。灃河發源于西安市長安區喂子坪鄉雞窩子以南的秦嶺北側，流經西安市長安區、戶縣，咸陽市秦都區等十余個鄉鎮，于咸陽北部魚王村匯入渭河，流域總面積1386 km2[1]。灃河多年平均年徑流量為2.381億m3，徑流年際變化大，年內分配也極不均勻[2]。灃河流域在灃河干流設有秦渡水文站，在大峪河峪口設有大峪水文站。秦渡鎮水文站控制面積566 km2。大峪水文站控制流域面積58.6 km2。

目前，灃河處于西安城市建設西邊緣區域，流域防洪形勢嚴峻，水資源開發利用不夠合理，生態環境有待進一步改善。隨著西安建造國際化大都市的步伐及旅游事業、文化產業、經濟發展和生態環境建設，對灃河流域水污染分析及水文特征分析研究就顯得十分必要。

1 不同代表年選擇

圖1 秦渡鎮水文站年最大徑流量頻率曲線由頻率曲線可以得到不同頻率的徑流量，其中P=20%

以秦渡鎮水文站1956年～2009年共54年實測多年年徑流量作為基礎數據，采用采用P-Ⅲ型水文頻率分析軟件對徑流量進行頻率分析，選取不同代表年，結果見圖1。（豐水年）時徑流量為3.30億 m3，P=50%（平水年）時徑流量為2.11億 m3，P=75%（枯水年）時徑流量為 1.71億 m3。最終確定灃河流域在2003年為豐水年，2005年為平水年，2004年為枯水年，然后進行不同代表年降雨量和泥沙量計算。

2 灃河全流域不同代表年計算結果

運用SWAT模型對灃河全流域2003年（豐水年）、2005年（平水年）、2004年（枯水年）三個不同代表年降雨量和泥沙量進行計算，得到的降雨量和泥沙量模擬結果見表1。

表1 不同代表年的計算結果

由表1可以看出在不同代表年，隨降水量和徑流量變化，泥沙量也在變化。其中豐水年、平水年、枯水年的降水量、泥沙量依次減小；泥沙量的變化趨勢和降水量有密切關系。

3 時間和空間分布

3.1 時間分布

灃河流域三個代表年不同月份降雨量和泥沙量模擬結果分別見表2～4。

由表2看出2003年的降雨量和泥沙量都集中在汛期（6～10月），各項模擬值在2003年的年內分布非常不均。汛期降水量和泥沙量在全年的比例分別為0.79和0.95，汛期降雨量和泥沙量的年內變化具有明顯對應關系，明顯大于非汛期。

表2 豐水年（2003）逐月計算結果

表3 平水年（2005）逐月計算結果

表4 枯水年（2004）逐月計算結果

由表3看出2005年的降雨量和泥沙量都集中在汛期（6～10月），各項模擬值在2005年的年內分布非常不均。汛期降水量和泥沙量在全年的比例分別為0.82和0.92，降雨量和泥沙量的年內變化具有明顯對應關系，明顯大于非汛期。

由表4看出2004年的降雨量、泥沙量和非點源污染物負荷都集中在汛期（6～10月），各項模擬值在2004年的年內分布非常不均。汛期降水量和泥沙量在全年的比例分別為0.69和0.78，汛期降雨量和泥沙量的年內變化具有明顯對應關系，明顯大于非汛期。

綜合以上分析：泥沙量是伴隨著降雨徑流過程特別是暴雨過程產生的，所以泥沙年內主要在汛期產生，集中在降雨量較大的月份。一般在這段時間內，會有多場大型暴雨發生，隨暴雨過程會產生很大的輸沙過程。降雨量、徑流量和泥沙量是明顯正相關的，降雨量越大，地表徑流越大，泥沙量也越大。以上不同代表年分析可以看出，泥沙在降雨較多的月份都比降雨較少的月份多，這符合泥沙的產生規律，即伴隨降雨徑流過程產生的。

3.2 空間分布

本文選擇的模擬年份為2003年（豐水年）、2004年（枯水年）和2005年（平水年），這三年期間灃河流域的土地利用、植被覆蓋情況不會有明顯大的變化，所以非點源污染負荷的空間變化將由流域內降雨量大小及其在空間上的分布情況決定。以下對照灃河流域圖對三個代表年的流域內降雨、徑流和泥沙的空間變化規律進行分析。

圖5 灃河流域圖

3.2.1 流域內降雨量的空間分布

圖6～圖8分別是2003年（豐水年）、2005年（平水年）和2004年（枯水年）各年降雨量分布圖。由圖可以看出，在該研究區內降雨量的分布是很不均勻的，降水量隨著地形由南向北逐漸減少。降雨量較大的地區集中在靠近秦嶺周邊植被覆蓋良好的山地地區，中下游居住和耕地地區降水量偏小。從整個流域來看，這三個代表年降雨量分布情況大致相同，各年降雨量最少的地區主要集中在長安區的杜曲、韋曲、郭杜三個村鎮，其次是長安區的斗門、細柳、高橋和戶縣的三個村鎮；降水量最大的區域集中在上游秦嶺周邊的山地地區以及長安區的五星、灤鎮、子午地區。

SWAT模型要求所使用的降雨量站點盡可能的均勻分布在研究區域內，否則系統就會使用相近站點的降雨數據代替模擬，這樣就會使模擬出來的降雨量數據產生較大偏差。

圖6 2003年降雨量分布圖 （單位：mm）

圖7 2005年降雨量分布圖 （單位：mm）

圖8 2004年降雨量分布圖 （單位：mm）

3.2.2 流域內徑流的空間分布

圖9～圖11分別是2003年（豐水年）、2005年（平水年）和2004年（枯水年）各年徑流分布圖。由圖可以看出，在該研究區內徑流分布隨著地形由南向北逐漸減少，與降雨量呈正相關關系。徑流最大的區域集中在上游秦嶺周邊的山地地區以及長安區的五星、灤鎮、子午地區，因為這一部分子流域的降雨量很大，且地形坡度大，所以相應的徑流很大，是主要的集水區域。另外從河段流量來看，在灃河干流上潏河入灃地區，因為有子午河、抱龍峪河、澗子河、滈河和潏河等支流的匯入，使得干流水量迅速增加，所以在潏河入灃這一段流域徑流很大。

圖9 2003年徑流深分布圖 （單位：mm）

圖10 2005年徑流深分布圖 （單位：mm）

圖11 2004年徑流深分布圖 （單位：mm）

3.2.3 流域內泥沙負荷的空間分布

圖12～圖14分別是2003年（豐水年）、2005年（平水年）和2004年（枯水年）各年泥沙負荷分布圖。不同的土地利用具有不同的泥沙負荷特性，這不僅與其自身性質有關，而且與土地利用在流域內的空間分布特征有關，每種土地利用都與地形或土壤類型相關聯。由圖看出，灃河流域的產沙區主要分布在流域中上游長安區的五臺、太乙、子午和灤鎮。這些地區靠近秦嶺的周邊山地，山地灌木草叢分布廣泛，土壤類型主要為棕壤性土、棕壤和褐土，地形坡度大于下游地區，土壤侵蝕量大，所以該區域泥沙產出是最大的。下游耕地和區民區地形平坦，土壤侵蝕量小，所以泥沙產量很小。

圖12 2003年泥沙負荷分布圖 （單位：t/ha）

圖13 2005年泥沙負荷分布圖 （單位：t/ha）

圖14 2004年泥沙負荷分布圖 （單位：t/ha）

4 結論

通過對灃河流域三個代表年不同月份降雨量和泥沙量模擬結果進行分析，發現汛期降雨量和泥沙量的年內變化具有明顯對應關系。從時間分布來看，隨降雨徑流特別是暴雨徑流過程產生的泥沙負荷在豐水年較大，而枯水年相對較小，年內泥沙負荷主要集中在汛期；通過空間分析指出了灃河流域內主要水文因子產出關鍵源區。根據模型模擬計算結果，結合國內外經驗提出灃河流域整體規劃以及防治措施與建議。

[1]高榕.西安市灃河流域水質變化特征及污染控制對策研究[D].西安：西安建筑科技大學，2003.

[2]朱記偉，雒望余，劉建林.灃河干流生態治理的必要性研究[J].新西部，2009，（22）：63.