汾河中游閘壩區河川徑流演變特性及驅動因素分析*

(山西省水利水電科學研究院，山西 太原 030002)

水是地球萬物生命之源，是生物圈的血液，是維持人類和生態系統的基礎，隨著經濟的發展，水資源已成為重要的戰略資源。河川徑流是水資源存在的重要形式，是人類社會生存和發展最主要的水源。20世紀以來，人類活動直接或間接地影響了河川徑流的質量、數量和徑流過程，而徑流的變化對整個水文水資源系統的演化起著主導作用，并對資源環境和區域經濟發展有著深遠的影響[1]。汾河流域是山西省生產力布局最密集、人口最集中的區域[2]，近年來，在氣候和人類活動影響下，汾河流域水資源情勢與水循環格局已顯著改變，導致中下游部分河段一度出現了斷流、泉水干涸、水質惡化等問題[3]，進一步加劇了該地區水資源供需矛盾。

為蓄滯雨洪，抬升地下水位，改善汾河生態環境，山西省在汾河中游二壩至義棠段堤內建攔河節制蓄水閘15座，閘壩工程運行后，閘壩區將具有工程調節作用，汾河中游面臨河道人工干擾強度大、河流系統完整性受到威脅等問題，該河段上密集閘壩群的合理調度方案以及下游河道水量如何保障等問題亟待研究，而閘壩間水資源的合理調度與優化配置均受上游來水規律的影響，因此研究汾河中游徑流變化規律及成因對解決河道片段化問題及15座閘壩的科學運行十分必要。國內學者針對汾河流域已開展了大量的研究工作，這些研究多集中在汾河上游或下游，且研究多以單一影響因素對徑流的影響為主，研究成果對汾河流域徑流演變特性及影響因素分析缺乏系統性和全面性。因此，本文以汾河二壩站、義棠站實測徑流系列為基礎，從年際與年內兩個時間尺度研究年徑流變化規律，基于Mann-Kendall法突變檢驗結果，建立降水-徑流雙累計曲線模型，對影響徑流的各項驅動因素進行定量分析，探求河川徑流量減少的原因，為緩解汾河流域水資源危機、保護生態環境提供依據。

1 資料與方法

本研究以汾河二壩水文站1965—2015年和義棠水文站1958—2015年實測年徑流量作為基礎數據，采用累計距平法、Mann-Kendall趨勢檢驗和突變檢驗法對徑流的年內分配與年際變化規律進行了分析，采用雙累計曲線模型分析了徑流變化的驅動力。

1.1 Mann-Kendall趨勢檢驗法

Mann-kendal趨勢檢驗法是一種非參數檢驗方法，其具有不需要樣本規律分布、不受異常值干擾等優點，可對降水、徑流、氣溫等非正態分布序列的變化趨勢進行分析和檢驗。其原理為：設(x1,x2,…xn)是n個獨立、隨機樣本，統計量S為：

其中sign()為符號函數。當(xi-xj)小于、等于或大于零時，sign(xi-xj)分別等于-1、0和1。M-K統計量公式S分別為：

Z為正表示趨勢增加，反之表示減小。當|Z|≥1.28、|Z|≥1.64、|Z|≥2.32時，分別表示序列趨勢通過了信度為90%、95%、99%的顯著性檢驗，說明序列趨勢顯著上升或下降。

1.2 Mann-Kendall突變檢驗法

設序列x1,x2,…xn，構造秩序列ri，ri為xi≥xj(1≤j≤i)的累計數，其統計量Sk及其均值和方差分別為：

將Sk標準化：

UFk組成一條UF曲線，其中UF1=0，將此方法用于反序列，計算UB曲線，其中UB1=0，給定顯著水平α，如果UF>0，表示序列趨勢呈現上升，反之呈下降，當UF值超過信度線上下限時，表示序列上升或下降趨勢達到顯著水平。繪制UF和UB曲線圖，則兩條曲線在置信區間內的交點為突變點。

1.3 雙累計曲線

人類活動影響時期實測水文變量與天然時期之間的差值主要包括兩部分，即人類活動影響和氣候變化影響，并提出了二者對徑流影響的分離分析方法[4-5]。計算過程如下：

ΔWT=WHR-WN

ΔWH=WHR-WHN

ΔWC=WHN-WN

式中ΔWT——徑流總量變化值；

ΔWH——徑流受人類活動影響變化量；

ΔWC——徑流受氣候影響變化量；

WN——天然徑流量；

WHR——人類活動影響時段實測的徑流量；

WHN——人類活動影響時段天然徑流量；

ηT、ηC——人類活動與氣候變化分別對徑流影響比例。

2 結果分析

2.1 河川徑流年內分配

河川徑流量的年內分配一般表現為明顯的洪水、枯水交替現象，這種變化因河流所處的地域不同而有所差異。因此，對同一流域，即使年徑流總量相差無幾，但如果該流域的河川徑流量年內分配差異顯著，對從該流域取水的城市生產和生活用水、農業灌溉用水等也會造成較大影響。因此，研究河川徑流的年內分配也是確保社會經濟可持續發展的基礎。

圖1 汾河二壩、義棠站徑流量年內變化曲線

山西省降雨的主要水汽來源是夏季風帶來的暖濕氣流，6—9月降水量占全年降水量的70%以上，7、8兩月降水尤為集中，從汾河二壩站和義棠站的徑流量年內變化曲線(圖1)可知，汾河中游徑流年內分配主要表現為典型的夏雨型，即月徑流量的豐枯程度和降水量的大小對應。汾河中游區從11月到次年5月徑流量均較小，2—3月的徑流量為一年中的最小值，從6月開始，隨著降雨量的增加徑流量迅速上升，8月前后徑流量達到一年中的最大值，之后開始逐漸縮減。總體上，汾河中游閘壩區河川徑流量主要集中在汛期，且除3月外義棠站月徑流量均大于汾河二壩站月徑流量，3月是冬小麥等農作物灌溉的高峰期，上游灌區灌溉用水量的增加是導致下游義棠站流量減少的原因之一。

2.2 年際變化規律分析

汾河二壩站和義棠站徑流年際特征值統計見表1，歷年徑流累計距平曲線如圖2所示。

表1 Mann-Kendall趨勢檢驗結果

注*表示通過了0.01的顯著性檢驗。

圖2 汾河二壩、義棠站年徑流量距平累計曲線

由圖2可以看出，汾河二壩站和義棠站年徑流量在20世紀70年代之前呈上升趨勢，而整個70年代為波動起伏狀態，1980—2010年呈逐漸減小趨勢，特別是1998—2008年的10年間，地表徑流大幅度減少，究其原因，主要是降水量減少所致，加之地下水大量超采，導致地下水位下降、泉水出流量劇減，流域產流也隨之急劇衰退，兩者共同作用造成河川徑流量大幅減少。2008年之后隨著外調水的增加和流域關井壓采、水源逐步置換，汾河中游河川徑流量才有所回升。根據Kendall趨勢檢驗結果(表1)，汾河二壩站1965—2015年、義棠站1958—2015年的河川徑流量整體處于下降趨勢，且這種趨勢通過了0.01的顯著性檢驗，表明汾河中游段河川徑流量下降非常顯著。

通過圖3、圖4顯示的汾河二壩站及義棠站年徑流量Mann-Kendall突變檢驗可以看出，汾河二壩站年徑流量1969—1975年表現為減少的趨勢，這種減少趨勢在1971年通過了0.05顯著性檢驗(U0.05=-1.96)，說明汾河二壩站年徑流量下降趨勢十分顯著；義棠站徑流量1969—1974年表現為減少的趨勢，但趨勢不顯著，1976—1979年徑流量略有回升，但從1980—1987年徑流量又表現為下降趨勢，且這種趨勢在1981年通過了0.05顯著性檢驗(U0.05=-1.96)，表明義棠站年徑流量下降趨勢十分顯著；從2008年開始兩站徑流量均開始回升，但未達到顯著性水平。

圖3 汾河二壩站年徑流量突變檢驗曲線

圖4 義棠站年徑流量突變檢驗曲線

圖3、圖4中UF(k)和UB(k)曲線在1970年和1971年在臨界線之間有一個交點，并且汾河二壩站UF(k)曲線在1971年、義棠站UF(k)曲線在1981年后超出了臨界線，因此汾河二壩站河川徑流量在1971年、義棠站河川徑流量在1981年后顯著下降，汾河二壩站在1970年、義棠站在1971年發生突變。

2.3 驅動因素分析

氣候變化和人類活動是導致徑流發生變化的雙重驅動因素[6]，降水-徑流雙累積曲線可反映出氣候變化與人類活動對徑流的影響。通過趨勢分析和突變檢驗，汾河二壩站、義棠站徑流序列均呈顯著遞減趨勢，突變年份分別為1970年和1971年，因此，汾河二壩站可將1970年以前的徑流量視為基準期，義棠站可將1971年以前的徑流量視為基準期，基準期內徑流變化只受氣候變化的影響，之后的時段視為變化期，變化期內徑流不僅受氣候變化的影響，同時還受人類活動的干擾。根據雙累計曲線建模原理，計算出降水、人類活動對徑流的影響，結果如圖5所示。

圖5 汾河二壩站、義棠站降水-徑流雙累計曲線

汾河二壩站和義棠站，降水因素對徑流的影響量分別為0.38億m3和0.74億m3，影響率分別為10.27%和12.71%；人類活動對徑流的影響量分別為3.32億m3和5.08億m3，影響率分別為89.73%和87.29%，見表2。

表2 降水、人類活動對徑流影響統計結果

流域徑流分布特性主要受氣候因素和人類活動兩個因素影響。氣候因素對徑流產生著持久且基本穩定的影響。人類活動對流域地表徑流量的影響可分為兩方面：?以地表水為水源的工農業用水、生活用水以及水利工程的調蓄，直接改變了下游水文站控制斷面的實測徑流量；?生產實踐和社會活動改變流域下墊面條件，導致入滲、徑流、蒸散發等平衡要素變化，使得流域產匯流條件發生變化，同時改變徑流的時程分配。

對比分析各驅動因素的影響結果得出如下結論：

a.義棠站的人類活動對徑流總影響量遠大于汾河二壩站，但影響程度略小于汾河二壩水文站，主要由于義棠水文站控制流域面積更大，范圍更廣。

b.從汾河二壩站和義棠站降水和人類活動對徑流的影響比例來看，雖然結果有微小差異，但總體上保持一致性，表明氣候和人類活動對汾河流域河川徑流的影響具有一致性。

c.降水量偏少以及時空分布不均是造成汾河中游徑流量減少的原因，但影響程度有限。

d.人類活動干擾是汾河中游徑流減少的主要原因。隨著人口的增加和經濟的發展，人為大量開采地下水，尤其是開采巖溶水，導致巖溶大泉自然出流量嚴重衰減，加之煤礦的無序開采對地下水造成了永久性和不可逆轉性的破壞，導致地下水位持續下降，使淺層地下水疏干，加速了土壤退化和植被衰退，植被破壞又進一步導致涵養水源能力減弱，地表濕地面積萎縮、地下水補給減少，加劇了地下水位快速下降，形成生態系統的惡性循環。

3 結 語

本文采用Mann-Kendall法對汾河中游二壩站1965—2015年、義棠站1958—2015年實測徑流量進行了趨勢檢驗和突變檢驗，系統地分析了汾河中游徑流年內分配和年際變化規律，并根據突變檢驗結果建立了降水-徑流雙累積曲線模型，研究表明:汾河中游段徑流量年內分配不均，年際變化大，年徑流量總體呈現減小趨勢，降水量偏少以及時空分布不均是造成汾河中游徑流量減少的部分原因，而人類活動干擾是汾河中游徑流減少的主要原因，特別是人為大量開采地下水，導致巖溶大泉自然出流量嚴重衰減，加之煤礦的無序開采造成地下水位持續下降，使淺層地下水疏干，加速了土壤退化和植被衰退，植被破壞又進一步導致涵養水源能力減弱，進一步加劇了地表徑流的衰減。