近50年清洋河水沙變化影響因素定量分析

孫曉蕾

(遼寧省鞍山水文局，遼寧 鞍山 114039)

我國現有水土流失面積近300萬km2，是世界上水土流失情況最嚴重的國家之一[1]。水土流失具有侵蝕耕地、淤塞河道、影響農業生產等危害，是我國經濟發展過程中的重大生態環境問題，是生態文明建設過程中的重要阻礙[2]。近年來，極端氣候頻現，人類農林牧漁等生產及工程建設等活動導致的河流水沙災害日益嚴重，因此，進行降雨變化和人類活動與水沙變化之間的定量分析，開展水沙災害的影響因素、預測研究，在水土保持理論研究和應用中都具有重要意義。

大沽夾河，位于山東省煙臺市，分為東支流和西支流。西支流亦為清洋河，干流長75km，流域面積為 1224km2。大沽夾河山區植被較少，水土流失嚴重，且河口段大潮進沙量大，容易造成河道泥沙淤積，汛期暴雨時洪水下泄，河道內水位急速上升而產生洪澇災害。清洋河流域內最重要的水利工程是門樓水庫，本研究選取門樓水庫水文站的水文資料作為清洋河流域水文要素情況的代表站[3]。

1 分析方法

1.1 Mann-Kendall檢驗模型

Mann-Kendall檢驗是分析水文序列趨勢的有效工具，對于任意時間序列xi(i=1，2，…，n)，Mann-Kendall[4-6]趨勢檢驗的統計量如下：

(1)

式中

(2)

當數據樣本數n≥8時，統計量S服從正態分布，方差為

(3)

標準化統計量服從正態分布，即

(4)

衡量趨勢大小指標為

(5)

式中，1

1.2 啟發式分割法

設x(t0)為時間序列x(t)的分割點，分割后的兩個子序列的平均值設為μ1(i)和μ2(i)，標準差設為S1(i)和S2(i)，則i點的合并偏差SD(i)為[7-9]

(6)

i點分割的兩個子序列均值差異檢驗標準TD(i)為

(7)

Tmax(i)對應兩個子序列差異最大的情況，統計概率P(Tmax)為

(8)

設定臨界值P0，當P(Tmax)≥P0，則在該點將數列分為兩個具有差異的子序列，并對分割得到的子序列進行循環分割，直至不可分割。設0為最小分割尺度，當序列長度小于0時，將不再進行分割。一般0.5≤P0≤0.95，0≥25。

1.3 有序樣本聚類檢驗法

有序樣本聚類檢驗法是對多變量的研究樣本數據進行多元統計的分析方法，可根據研究樣本數據的離散程度或相似系數，推求最優的分段點，將樣本數據分成若干子序列[10]，使同一子序列的數據之間離差平方和較小，而不同子序列數據之間的離差平方和較大[11]。

假設τ為序列xi(t=1，2，…，n)的分段點，則子序列的離差平方和為

(9)

與

(10)

總離差平方和為

Sn(τ)=Vτ+Vn-τ

(11)

式中，Sn取最小值時對應的τ即為最優分段點。

1.4 雙累積曲線法

Searcy[12]在Merriam[13]研究的基礎上，對雙累積曲線法進行了拓展研究，使其廣泛應用于水文要素的一致性檢驗、變化規律檢驗、缺值插補等水文要素趨勢性變化分析。

(12)

(13)

式中，X′，Y′分別為X，Y的累積量；Xi為第i年的參考變量；Yi為第i年的檢驗變量。

2 氣候變化對清洋河徑流影響硏究分析

2.1 降水、蒸發和徑流趨勢

1963—2016年歷年實測年降水量、徑流量和輸沙量如圖1所示。

圖1 1963—2016年清洋河實測年降水量、徑流量、輸沙量圖

利用Mann-Kendall模型對降水量、徑流量和輸沙量的趨勢進行檢驗，結果顯示降水量、徑流量均呈下降趨勢，但年際變化不明顯，而輸沙量年際間呈顯著下降的趨勢，詳見表1。

表1 清洋河年降水量、徑流量、輸沙量檢驗結果

2.2 年際要素變點分析

為了更準確確定年際輸沙量的突變年份，采用啟發式分割法和有序樣本聚類檢驗法進行計算分析，相關結果如圖2所示。由此可以看出，采用啟發式分割法1976年的T值最大，且P(Tmax)≥P0，采用有序樣本聚類檢驗法S0(τ)最小值發生在1976年。兩種方法均推得年際輸沙量突變年份為1976年。

圖2 清洋河徑流量的啟發式分割法和有序聚類分割法計算結果示意圖

同理，采用啟發式分割法和有序樣本聚類檢驗法計算確定1976年為徑流量明顯變點。則年際輸沙量和徑流量序列分割為1963—1976年和1977—2016年兩個階段，相關水文氣候要素統計結果見表2。

表2 清洋河不同階段水文氣候變量統計結果

2.3 降雨變化和人類活動對水沙影響的定量區分

利用Excel軟件繪制1963—2016年實際測量值的累積降水量-累積徑流量和累積降水量-累積輸沙量的雙累積曲線，如圖3—4所示。

圖3 清洋河1963—2016年累積降水量-累積徑流量雙累積曲線

圖4 清洋河1963—2016年累積降水量-累積輸沙量雙累積曲線

將累積降水量-累積徑流量曲線和累積降水量-累積輸沙量的雙累積曲線在計算突變點1976年分別分成兩段數據散點，并繪制趨勢線，計算回歸方程和方差。根據圖3—4，兩條雙累積曲線均在1976年發生轉折，趨勢線斜率均變小。

根據1963—1976年回歸方程預測1977—2016年的徑流量和輸沙量理論預測均值，則降水量變化導致的定量影響為不同時間跨度的理論均值差值，人類活動導致的定量影響為理論預測均值與實測均值的差值。

由表3計算結果可知，1977—2016 年的徑流量實測均值比1963—1976年減少1.03×108m3，其中降水量變化因素占76.70%，人類活動因素占23.30%；1977—2016 年的輸沙量實測均值比1963—1976年減少7.06×104t/a，其中降水量變化因素占29.89%，人類活動因素占70.11%。

表3 人類活動與降水量變化對清洋河水文要素影響

3 結論

(1)Mann-Kendall檢驗顯示年際間降水量、徑流量均呈微弱下降趨勢，輸沙量呈顯著下降趨勢，并使用啟發式分割法和有序樣本聚類檢驗法確定1976年為徑流量和輸沙量趨勢變化明顯變點。

(2)大沽夾河流域降雨呈逐年遞減趨勢，降雨分布極其不均，特別是1986、1989和1999年的降水量達到歷史極低值，對年際間徑流量降低趨勢起主要影響；另外，生活、工業、農業用水的取用量逐年遞增，是徑流量減少的次要因素。徑流量年際間呈下降趨勢對輸沙量的變化起一定影響，但主要是由取水攔沙、河道整治、退耕還林等人類活動對輸沙量的顯著下降趨勢起決定性作用。

(3)河流水沙變化是復雜的動態變化過程，本研究為具有相似氣候和地理條件的河流體系水沙變化分析做出一次有益的嘗試，提供一種可供借鑒的定量分析方法，能否廣泛推廣還需進一步研究探討。在后續的研究中，可將地下水開采、水利工程興建等因素納入分析范圍，更加全面準確地分析清洋河水沙變化，使該方法具有更廣的借鑒性和更強的實用性。