吳奕 宋瑞鵬 張紅衛(wèi) 陳莉 李洋 肖航

摘 要：以河南省1956—2016年降水量、地表水資源量數(shù)據(jù)為基礎(chǔ)，采用距平累積法、Mann-Kendall檢驗(yàn)法、雙累積曲線法、累積量斜率變化率比較法判定河南省及其各流域降水量、地表水資源量變化趨勢(shì)，分析降水量、地表水資源量演變關(guān)系。結(jié)果表明：河南省1956—2016年降水量、地表水資源量均呈減小趨勢(shì)，地表水資源量減小趨勢(shì)更顯著，北部海河、黃河流域地表水資源量減小趨勢(shì)比南部淮河、長(zhǎng)江流域顯著;20世紀(jì)90年代，該省海河、黃河流域一些站降水—徑流關(guān)系發(fā)生改變。以海河流域新村站為例，量化降水和人類(lèi)活動(dòng)對(duì)徑流量的影響程度，發(fā)現(xiàn)降水對(duì)徑流量變化的影響較小、人類(lèi)活動(dòng)對(duì)徑流量變化的影響較大。

關(guān)鍵詞：降水量;地表水資源量;Mann-Kendall檢驗(yàn)法;累積量斜率變化率比較法;河南省

中圖分類(lèi)號(hào)：TV213.4

文獻(xiàn)標(biāo)志碼：A

doi：10.3969/j.issn.1000-1379.2021.11.017

引用格式：吳奕，宋瑞鵬，張紅衛(wèi)，等.河南省降水量、地表水資源量變化趨勢(shì)及演變關(guān)系[J].人民黃河，2021，43（11）：92-96.

Variation Trend and Relationship of Annual Precipitation and Surface Water Resources in Henan Province

WU Yi1， SONG Ruipeng2， ZHANG Hongwei1， CHEN Li1， LI Yang1， XIAO Hang1

（1.Henan Hydrology and Water Resources Bureau， Zhengzhou 45000 ?China;

2.Hydrological Bureau， Yellow River Conservancy Commission， Zhengzhou 450004， China）

Abstract： According to the data of precipitation and surface water resources in Henan Province from 1956 to 2016， the accumulated deviation， Mann-Kendall test， rainfall-runoff double mass curve， comparative method of cumulative slope change ratio were used to indicate the variation trend of precipitation and surface water resources and the evolution relations between them. Both precipitation and surface water resources had a decreasing trend， while the decreasing trend was more significant in surface water resources. Surface water resources in the Haihe River and the Yellow River Basin decreased more significant than that in the Huaihe River and the Yangtze River Basin. Since the 1990s， rainfall-runoff relationship had changed in some hydrological stations of the Haihe River and the Yellow River Basin . Taking Xincun hydrological station as an example， quantified the influence of precipitation and human activities on runoff. The results show that precipitation has little influence on runoff while human activities is the main factor of runoff variation.

Key words： precipitation; surface water resources; Mann-Kendall test; comparative method of cumulative slope change ratio; Henan Province

河南省地處中原，總面積16.7萬(wàn)km2。地勢(shì)西高東低，北、西、南三面環(huán)山，間有斷陷盆地，東部為寬闊平原。河南省地跨海河、黃河、淮河、長(zhǎng)江四大流域，其中：海河流域1.53萬(wàn)km2，占全省總面積的9.2%;黃河流域3.62萬(wàn)km2，占全省總面積的21.7%;淮河流域8.83萬(wàn)km2，占全省總面積的52.8%;長(zhǎng)江流域2.72萬(wàn)km2，占全省總面積的16.3%。豫南、豫西南水資源相對(duì)豐沛，占全省水資源總量的一半以上;豫東、豫北水資源匱乏，與耕地及人口分布不相匹配。

近年來(lái)，河南省許多河流的徑流衰減現(xiàn)象明顯。河川徑流源于降水，氣候變化、人類(lèi)活動(dòng)對(duì)區(qū)域水資源狀況有直接影響[1-3]。一些學(xué)者對(duì)徑流變化趨勢(shì)、徑流變化成因進(jìn)行了研究[4-8]。目前，尚未有關(guān)于河南省境內(nèi)各流域降水量、地表水資源量演變趨勢(shì)及其關(guān)系的分析。本文基于1956—2016年河南省各流域年降水量、地表水資源量，分析降水量、地表水資源量變化趨勢(shì)，厘清兩者之間的演變關(guān)系，以期為全省水資源調(diào)度、防汛抗旱提供參考。

1 研究方法與數(shù)據(jù)來(lái)源

1.1 Mann-Kendall檢驗(yàn)

Mann-Kendall（M-K）檢驗(yàn)法是由Mann和Kendall提出的，用于檢測(cè)序列的變化趨勢(shì)，被廣泛應(yīng)用于水文和氣象序列的趨勢(shì)檢驗(yàn)中[9-11]。對(duì)于時(shí)間序列X，M-K趨勢(shì)檢驗(yàn)統(tǒng)計(jì)量的表達(dá)式[12-13]為

S=∑n-1i=1∑nj=i+1sgn（xj-xi）（1）

式中：S為M-K檢驗(yàn)法統(tǒng)計(jì)量;xi、xj分別為時(shí)間序列的第i個(gè)、第j個(gè)數(shù)據(jù)值;n為數(shù)據(jù)樣本數(shù)量。

當(dāng)n>10時(shí)，統(tǒng)計(jì)量S大致服從正態(tài)分布，方差為

Var（S）=n（n-1）（2n-5）/18（2）

其正態(tài)分布的檢驗(yàn)統(tǒng)計(jì)量為Z：

Z=S-1Var（S） （S>0）

0（S=0）

S+1Var（S）（S<0）（3）

雙邊趨勢(shì)檢驗(yàn)中，Z為正表示序列具有增大趨勢(shì)，Z為負(fù)表示序列具有減小趨勢(shì)。在給定的α置信水平上，若Z的絕對(duì)值大于或等于1.28、1.64和2.32，則分別表示通過(guò)了置信度90%、95%和99%的顯著性趨勢(shì)檢驗(yàn)，分別達(dá)到顯著、十分顯著、極顯著水平。

1.2 距平累積法

距平累積法是一種根據(jù)曲線直觀判斷變化趨勢(shì)的方法[14]。對(duì)于序列X，某一時(shí)刻t的距平累積值為

x∧t=∑ti=1（xi-x-） （t= …，n）（4）

式中：x-=1n∑ni=1xi。

將序列的距平累積值全部算出，即可繪出距平累積曲線[15]。距平值隨年份累加，曲線上升階段表示距平值為正，序列呈上升趨勢(shì);曲線下降階段表示距平值為負(fù)，序列呈下降趨勢(shì)[16]。

1.3 降水量—徑流量雙累積曲線

假定降水量與徑流量有同步成因關(guān)系，在直角坐標(biāo)系中繪制同期降水量累積值與徑流量累積值關(guān)系曲線，曲線斜率發(fā)生明顯變化處對(duì)應(yīng)的年份就是降水量、徑流量關(guān)系變化的年份[17]。

1.4 累積量斜率變化率比較法

累積量斜率變化率比較法是近年提出的定量分析降水及人類(lèi)活動(dòng)對(duì)徑流變化貢獻(xiàn)率的一種方法[18]。在直角坐標(biāo)系中分別擬合累積降水量—年份、累積徑流量—年份的線性關(guān)系，假設(shè)累積徑流量—年份線性關(guān)系的斜率突變前后兩個(gè)時(shí)期分別為SRb和SRa，累積降水量—年份線性關(guān)系的斜率突變前后兩個(gè)時(shí)期分別為SPb和SPa，則累積徑流量斜率變化率為

RSR=（SRa-SRb）/SRb（5）

累積降水量斜率變化率為

RSP=（SPa-SPb）/SPb（6）

降水量變化對(duì)徑流量變化的貢獻(xiàn)率為

CP=RSP/RSR （7）

人類(lèi)活動(dòng)對(duì)徑流量變化的貢獻(xiàn)率為

Cm=1-CP （8）

自變量為年份，因變量為累積量，累積量的引入在一定程度上消除了實(shí)測(cè)數(shù)據(jù)年際波動(dòng)的影響;自變量的客觀性提高了累積量—年份之間的相關(guān)性，克服了降水—徑流回歸曲線相關(guān)系數(shù)不高的缺陷[19]。

1.5 數(shù)據(jù)來(lái)源

1956—2016年序列降水量、地表水資源量數(shù)據(jù)來(lái)源于河南省水文部門(mén)相關(guān)成果。本文以1956—2016年序列為研究對(duì)象，并與1956—2000年序列對(duì)比，分析河南省降水量和地表水資源量變化趨勢(shì)及兩者之間的演變關(guān)系。

2 降水量、地表水資源量變化趨勢(shì)

2.1 降水量變化趨勢(shì)

1956—2016年河南省多年平均降水量768.5 mm，與1956—2000年的多年平均降水量774.7 mm相比減小0.80%。河南省各流域都有減小，減小幅度最大的為長(zhǎng)江流域（1.44%），其次為海河流域（1.04%），淮河流域減小0.75%，黃河流域減小0.24%。降水量多年均值變化情況見(jiàn)表1。

繪制河南省及其各流域降水量、地表水資源量1956—2016年序列距平累積曲線（見(jiàn)圖1～圖5），并結(jié)合全省及各流域降水量數(shù)據(jù)可知，河南省及其各流域1956—1964年降水量偏豐，全省降水量較多年（1956—2016年）均值768.5 mm增大10.30%，各流域較多年均值增大6.33%～25.32%（海河流域增幅最大，長(zhǎng)江流域增幅最小）;海河流域1978—2002年以及黃河流域1986—2002年是河南省各流域1956—2016年序列中降水量減小最明顯且持續(xù)時(shí)間最長(zhǎng)的時(shí)段，分別比流域多年均值偏小9.27%、9.07%;20世紀(jì)80年代中期至90年代中期是淮河、長(zhǎng)江流域降水量波動(dòng)減小相對(duì)較長(zhǎng)的時(shí)期，分別比多年均值偏小7.54%、6.42%;2011—2016年淮河、長(zhǎng)江流域降水量持續(xù)低于或接近多年均值，分別比多年均值偏小11.73%、13.33%。

利用M-K趨勢(shì)檢驗(yàn)法對(duì)河南省及其各流域1956—2000年、1956—2016年序列降水量變化趨勢(shì)進(jìn)行顯著性檢驗(yàn)，結(jié)果見(jiàn)表2。兩序列統(tǒng)計(jì)量Z均為負(fù)值，表明兩序列總體呈減小趨勢(shì)。1956—2000年序列河南省及其各流域檢驗(yàn)統(tǒng)計(jì)量Z均沒(méi)有達(dá)到90%的置信度，表明減小趨勢(shì)不顯著;1956—2016年序列海河流域檢驗(yàn)統(tǒng)計(jì)量Z達(dá)到90%的置信度，表明減小趨勢(shì)顯著，但全省及其他流域降水量減小趨勢(shì)未達(dá)到顯著性水平。1956—2016年序列比1956—2000年序列降水量減小趨勢(shì)更明顯，但也未達(dá)到顯著性減小的水平。

2.2 地表水資源量變化趨勢(shì)

1956—2016年河南省多年平均地表水資源量289.30億m 與1956—2000年多年平均地表水資源量299.65億m3相比減小3.45%（見(jiàn)表1）。河南省各流域中，海河流域地表水資源量減小幅度最大（7.24%），黃河流域減小2.57%，淮河流域減小3.19%，長(zhǎng)江流域減小3.92%。

由圖1～圖5可以看出，河南省及其各流域地表水資源量變化趨勢(shì)與降水量相似，但變幅大于降水量。結(jié)合河南省及其各流域地表水資源量數(shù)據(jù)可知，1956—1964年為持續(xù)偏豐較長(zhǎng)的時(shí)段，較多年（1956—2016年）均值289.30億m3偏大31.49%，各流域較多年均值偏大24.05%～99.94%（海河流域增幅最大，淮河流域增幅最?。?。海河流域1978—2016年以及黃河流域1986—2002年是地表水資源量減小最明顯、持續(xù)時(shí)間最長(zhǎng)的時(shí)段，分別比流域多年均值偏小24.57%、23.25%。20世紀(jì)80年代中期至90年代中期是淮河、長(zhǎng)江流域地表水資源量波動(dòng)減小相對(duì)較長(zhǎng)的時(shí)期，分別比多年均值偏小18.90%、21.29%。2011—2016年是全省地表水資源量最枯的時(shí)段，較多年均值偏小37.25%，各流域較多年均值偏小17.48%～47.58%（長(zhǎng)江流域減幅最大，黃河流域減幅最小）。

利用M-K趨勢(shì)檢驗(yàn)法對(duì)河南省及其各流域1956—2016年、1956—2000年序列地表水資源量變化趨勢(shì)進(jìn)行顯著性檢驗(yàn)，結(jié)果見(jiàn)表2。兩序列統(tǒng)計(jì)量Z均為負(fù)值，表明兩序列總體呈減小趨勢(shì)。全省1956—2000年序列檢驗(yàn)統(tǒng)計(jì)量Z沒(méi)有達(dá)到90%的置信度，表明地表水資源量減小趨勢(shì)不顯著;海河、黃河流域檢驗(yàn)統(tǒng)計(jì)量Z達(dá)到95%的置信度，表明減小趨勢(shì)十分顯著。1956—2016年序列河南省及其海河、黃河、長(zhǎng)江流域檢驗(yàn)統(tǒng)計(jì)量Z均達(dá)到95%的置信度，表明減小趨勢(shì)十分顯著，其中海河流域減小趨勢(shì)最為明顯;淮河流域檢驗(yàn)統(tǒng)計(jì)量達(dá)到90%的置信度，較其他流域減小趨勢(shì)平緩。1956—2016年序列比1956—2000年序列地表水資源量減小趨勢(shì)總體更明顯，減小趨勢(shì)由不顯著演變?yōu)槭诛@著。

3 降水量和地表水資源量演變關(guān)系

從1956—2016年降水量、地表水資源量變化趨勢(shì)可看出，降水量較大的年份地表水資源量也較大。降水量和地表水資源量變化趨勢(shì)總體上同步，但地表水資源量的年際變幅大于降水量的，減小趨勢(shì)也較降水量更顯著，且具有明顯的區(qū)域特征。河南省海河、黃河流域比其他流域地表水資源量減小趨勢(shì)更顯著、持續(xù)時(shí)間更長(zhǎng)，這兩個(gè)流域1956—2016年與1956—2000年相比，地表水資源量減小幅度分別為降水量減小幅度的7倍和11倍，而全省及其他兩個(gè)流域在3～4倍之間。20世紀(jì)90年代海河、黃河流域一些水文站降水—徑流雙累積曲線發(fā)生明顯拐點(diǎn)，即降水—徑流關(guān)系發(fā)生改變，主要原因是經(jīng)濟(jì)快速發(fā)展、城市化進(jìn)程加快引起區(qū)域水循環(huán)特征改變，以及人類(lèi)活動(dòng)對(duì)水資源的開(kāi)發(fā)利用引起產(chǎn)匯流條件改變。

以海河流域新村水文站為例[20]，根據(jù)其多年天然河川徑流量和降水量，繪制降水—徑流雙累積曲線（見(jiàn)圖6），雙累積曲線突變年份為1998年，表明1998年以后降水—徑流關(guān)系發(fā)生了改變。分別繪制累積徑流量曲線、累積降水量曲線，對(duì)1998年前后曲線進(jìn)行分段線性擬合（見(jiàn)圖7、圖8），可以看出，擬合關(guān)系式的判定系數(shù)較大，擬合效果較好。利用累積量斜率變化率比較法，量化降水和人類(lèi)活動(dòng)對(duì)徑流量變化的影響程度，降水對(duì)徑流量變化的貢獻(xiàn)率為13%，人類(lèi)活動(dòng)對(duì)徑流量變化的貢獻(xiàn)率為87%，表明降水對(duì)徑流量變化影響較小，徑流量的變化主要受人類(lèi)活動(dòng)影響，流域下墊面條件發(fā)生了明顯變化。

4 結(jié) 論

基于河南省1956—2016年序列降水量與地表水資源量相關(guān)資料，采用距平累積法、M-K趨勢(shì)檢驗(yàn)法，分析了河南省及其各流域降水量、地表水資源量變化趨勢(shì)及演變關(guān)系。以新村站為例，利用雙累積曲線法、累積量斜率變化率比較法對(duì)地表水資源量減小幅度大于降水量減小幅度進(jìn)行了成因量化分析，結(jié)論如下。

（1）河南省多年降水量、地表水資源量均呈減小趨勢(shì)，1956—2016年與1956—2000年序列相比，降水量、地表水資源量的減小趨勢(shì)進(jìn)一步加劇，但降水量減小趨勢(shì)沒(méi)有達(dá)到顯著性水平，地表水資源量減小趨勢(shì)顯著。

（2）1956—2016年河南省海河、黃河流域地表水資源量減小趨勢(shì)較淮河、長(zhǎng)江流域更顯著。

（3）河南省海河、黃河流域一些水文站20世紀(jì)90年代降水—徑流關(guān)系發(fā)生改變，降水對(duì)徑流量變化的影響較小，人類(lèi)活動(dòng)對(duì)徑流量和地表水資源量減小的影響較大。

參考文獻(xiàn)：

[1] 王國(guó)慶，張建云，管曉祥，等.中國(guó)主要江河徑流變化成因定量分析[J].水科學(xué)進(jìn)展，2020，31（3）：313-323.

[2] 任國(guó)玉.氣候變化與中國(guó)水資源[M].北京：氣象出版社，2007：47-323.

[3] 張國(guó)棟，張照璽，余韻，等.汾河上游土地利用變化對(duì)徑流的影響研究[J].人民黃河，2020，42（10）：29-33.

[4] 王國(guó)慶，李迷，金君良，等. 涪江流域徑流變化趨勢(shì)及其對(duì)氣候變化的響應(yīng)[J].水文，2012，32（1）：22-28.

[5] 劉洪波，菅浩然.氣候和下墊面要素對(duì)涇河流域徑流變化的影響[J].人民黃河，2021，43（1）：22-28.

[6] 師潤(rùn)，田鵬，趙廣舉，等.南北過(guò)渡帶典型流域徑流變化歸因?qū)Ρ妊芯縖J].人民黃河，2020，42（12）：29-35.

[7] 譚云娟，邱新法，曾燕，等. 近50 a來(lái)中國(guó)不同流域降水的變化趨勢(shì)分析[J].氣象科學(xué)，2016，36（4）：494-501.

[8] 劉曉麗，陳明哲，汪子雄，等.淮河流域中上游徑流變化歸因分析[J].人民黃河，2020，42（10）：16-22.

[9] 徐宗學(xué)，張楠.黃河流域近50年降水變化趨勢(shì)分析[J].地理研究，2006，25（1）：27-32.

[10] 王樂(lè)，劉德地，李天元，等.基于多變量M-K檢驗(yàn)的北江流域降水趨勢(shì)分析[J].水文，2015，35（4）：85-90.

[11] 高彥春，王金鳳，封志明.白洋淀流域氣溫、降水和徑流變化特征及其相互響應(yīng)關(guān)系[J].中國(guó)生態(tài)農(nóng)業(yè)學(xué)報(bào)，2017，25（4）：467-476.

[12] 丁晶，鄧育仁.隨機(jī)水文學(xué)[M]成都：成都科技大學(xué)出版社，1988：25-35.

[13] 郝璐，李王成，李晨，等.近51 a銀川市降雨變化特征分析[J].中國(guó)農(nóng)村水利水電，2020（8）：153-158.

[14] 姜瑤，徐宗學(xué)，王靜.基于年徑流序列的五種趨勢(shì)檢驗(yàn)方法性能對(duì)比[J].水利學(xué)報(bào)，2020，51（7）：845-857.

[15] 傅小城，王芳，王浩，等.柴達(dá)木盆地氣溫降水的長(zhǎng)序列變化及與水資源關(guān)系[J].資源科學(xué)，2011，33（3）：408-415.

[16] 張應(yīng)華，宋獻(xiàn)方.水文氣象序列趨勢(shì)分析與變異診斷的方法及其對(duì)比[J].干旱區(qū)地理，2015，38（4）：652-665.

[17] 穆興民，張秀琴，高鵬，等.雙累積曲線方法理論在水文氣象領(lǐng)域應(yīng)用中應(yīng)注意的問(wèn)題[J].水文，2010，30（4）：47-50.

[18] 王彥君，王隨繼，蘇騰.1955—2010年松花江流域不同區(qū)段徑流量變化影響因素定量評(píng)估[J].地理科學(xué)進(jìn)展，2012，33（1）：65-75.

[19] 王隨繼，閻云霞，顏明.皇甫川流域降水和人類(lèi)活動(dòng)對(duì)徑流量變化的貢獻(xiàn)分析：累積量斜率變化率比較方法的提出及應(yīng)用[J].地理學(xué)報(bào)，2012，67（3）：388-397.

[20] 吳奕，張紅衛(wèi)，趙清虎.徑流量突變點(diǎn)的檢驗(yàn)及其成因分析[G]//水利部信息中心.水文水資源監(jiān)測(cè)與評(píng)價(jià)應(yīng)用技術(shù)論文集.南京：河海大學(xué)出版社，2019：136-141.

【責(zé)任編輯 張華興】