趙自陽(yáng)，李王成,2,3

(1．寧夏大學(xué)土木與水利工程學(xué)院，銀川 750021; 2．寧夏節(jié)水灌溉與水資源調(diào)控工程技術(shù)研究中心，銀川 750021; 3. 旱區(qū)現(xiàn)代農(nóng)業(yè)水資源高效利用教育部工程研究中心，銀川 750021)

隨著近些年社會(huì)的快速發(fā)展，水資源緊缺形勢(shì)愈加嚴(yán)峻。為了對(duì)區(qū)域水資源狀況進(jìn)行合理的評(píng)價(jià)分析，國(guó)內(nèi)外學(xué)者在許多領(lǐng)域進(jìn)行了相關(guān)研究。根據(jù)原理的不同大致可分為數(shù)理統(tǒng)計(jì)方法[1,2]、運(yùn)籌決策方法[3]、灰色系統(tǒng)方法[4,5]、模糊數(shù)學(xué)方法[6,7]和機(jī)器學(xué)習(xí)方法[8]等。但是這些方法的適用條件各不相同，計(jì)算結(jié)果也有一定誤差。1992年，William和Wackernagel提出了生態(tài)足跡模型，可以對(duì)自然資源利用狀況進(jìn)行相關(guān)評(píng)價(jià)。由于其綜合性強(qiáng)和易于操作等特點(diǎn)，在水資源領(lǐng)域得到迅速應(yīng)用，并逐漸發(fā)展成為水資源生態(tài)足跡模型[9]和基于虛擬水理論的水足跡模型[10]。兩者的區(qū)別主要有以下三點(diǎn)[11]：①核算視角：水資源生態(tài)足跡是從水資源供需角度出發(fā)，而虛擬水足跡則是從生產(chǎn)和消費(fèi)的角度出發(fā)；②核算終點(diǎn)：水資源生態(tài)足跡是把水資源的使用換算為水資源用地面積；而虛擬水足跡則是把人類對(duì)各種產(chǎn)品和服務(wù)的消費(fèi)換算為相應(yīng)的水資源體積。③核算目的：水資源生態(tài)足跡與原始生態(tài)足跡模型最為接近，主要探討區(qū)域水資源承載力情況，而虛擬水足跡主要是為了探討生產(chǎn)、消費(fèi)和產(chǎn)業(yè)結(jié)構(gòu)之間的關(guān)聯(lián)。

目前，國(guó)內(nèi)學(xué)者對(duì)水資源生態(tài)足跡的研究主要集中在水資源生態(tài)足跡賬戶的劃分[12,13]、水資源生態(tài)足跡研究與預(yù)測(cè)[14,15]和水資源生態(tài)承載力分析[16,17]等幾個(gè)方面。河南省是我國(guó)的農(nóng)業(yè)大省，水資源緊缺問(wèn)題日益突出，目前對(duì)于河南省水資源生態(tài)足跡的研究一般僅集中在淡水資源生態(tài)足跡方面。而2011年中央水利會(huì)議在1號(hào)文件中提出用水總量控制紅線、用水效率控制紅線和限制納污紅線的“三條紅線”規(guī)定[18]?；诖耍疚膶⑺廴炯{入水資源生態(tài)足跡模型，使得淡水生態(tài)足跡和水污染生態(tài)足跡分別對(duì)應(yīng)“三條紅線”中的用水控制紅線和限制納污紅線。從而對(duì)河南省2005-2014年間水資源生態(tài)足跡特征進(jìn)行分析，衡量河南省水資源利用現(xiàn)狀，為政府部門評(píng)價(jià)水資源狀況，科學(xué)管理水資源提供依據(jù)。

1 研究區(qū)概況

河南省位于我國(guó)中部偏東、黃河中下游地區(qū)，全省面積約16.7 萬(wàn)km2。由南向北平均氣溫為15.7～12.1 ℃，年均降雨量1 380.6～532.5 mm，降雨以6-8月份最多，年均日照1 848.0～2 488.7 h，河南地跨長(zhǎng)江、淮河、黃河、海河四大流域。流域面積100 km2以上的河流有493條。全省多年平均水資源總量405 億m3，居全國(guó)第19位，人均水資源占有量不足420 m3，相當(dāng)于全國(guó)平均水平的1/5。

2 研究方法

2.1 水資源生態(tài)足跡模型

根據(jù)水資源生態(tài)足跡的內(nèi)涵，本研究將水資源生態(tài)足跡賬戶[19]劃分為淡水生態(tài)足跡和水污染生態(tài)足跡。

(1)淡水生態(tài)足跡。根據(jù)河南省水資源利用現(xiàn)狀，將淡水生態(tài)足跡分為農(nóng)業(yè)用水、工業(yè)用水、生活用水和生態(tài)環(huán)境補(bǔ)水4類。水資源生態(tài)足跡的計(jì)算模型[20]如下表示：

(2)

式中：i為水資源利用類型；EFa為淡水生態(tài)足跡，hm2；EFi為第i類水資源生態(tài)足跡，hm2；N為人口數(shù)；efi為第i類水資源人均生態(tài)足跡，hm2/cap；γ為水資源全球均衡因子；W為消耗的水資源量，m3；pw為水資源全球平均生產(chǎn)能力，m3/hm2。

(2)水污染生態(tài)足跡[19]。水污染生態(tài)足跡主要是基于將污染物稀釋至標(biāo)準(zhǔn)所需水資源量的“灰色水”理論。本文以化學(xué)需氧量(COD)排放指標(biāo)對(duì)水污染生態(tài)足跡進(jìn)行計(jì)算，在考慮水資源全球均衡因子情況下的計(jì)算模型如下：

(3)

所以最終總的水資源生態(tài)足跡就是淡水生態(tài)足跡和水污染生態(tài)足跡之和，如下所示：

EF=EFa+EFb

(4)

2.2 水資源生態(tài)承載力計(jì)算模型

水資源承載力是指一個(gè)地區(qū)在某個(gè)時(shí)期的水資源量可供當(dāng)?shù)厣鐣?huì)、經(jīng)濟(jì)和環(huán)境可持續(xù)發(fā)展的支撐能力。計(jì)算模型[21]如下：

(5)

式中：ECw為水資源承載力，hm2；N為人口數(shù)；ecw為人均水資源承載力，hm2/cap；γ為水資源全球均衡因子；ψ為區(qū)域水資源的產(chǎn)量因子；Q為水資源總量，m3；pw為水資源全球平均生產(chǎn)能力，m3/hm2；而一個(gè)地區(qū)的水資源中用于人類活動(dòng)的只占40%，所以在計(jì)算水資源承載力的時(shí)候要乘以系數(shù)0.4。

2.3 水資源生態(tài)盈余和生態(tài)赤字

為了評(píng)價(jià)一個(gè)地區(qū)的水資源狀況，有關(guān)專家學(xué)者提出了生態(tài)赤字和生態(tài)盈余的概念[21]，主要是將區(qū)域水資源生態(tài)足跡和生態(tài)承載力進(jìn)行比較，看其結(jié)果大小。具體如下所示：

水資源生態(tài)盈余(赤字)=ECw-EF

(6)

當(dāng)結(jié)果為正時(shí)，水資源生態(tài)盈余,表明區(qū)域水資源良性發(fā)展；當(dāng)為零時(shí)，水資源生態(tài)平衡；當(dāng)小于零時(shí)，水資源生態(tài)赤字,表明不利于區(qū)域可持續(xù)發(fā)展。

2.4 萬(wàn)元GDP水資源生態(tài)足跡

萬(wàn)元GDP水資源生態(tài)足跡[21]主要用于衡量水資源的利用效率，它是指區(qū)域水資源生態(tài)足跡和當(dāng)?shù)谿DP的比值。并且與水資源利用效率成反比，如下所示：

(7)

2.5 水資源負(fù)載指數(shù)

水資源負(fù)載指數(shù)[21]主要是將一地區(qū)的水資源量、人口、經(jīng)濟(jì)和降雨量聯(lián)系起來(lái)，反映區(qū)域水資源利用概況，并對(duì)今后的水資源開發(fā)難易程度做出評(píng)價(jià)。計(jì)算公式如下：

(8)

式中：c為水資源負(fù)載指數(shù)；p為人口數(shù)，萬(wàn)人；G為地區(qū)生產(chǎn)總值，億元；Q為水資源總量，億m3；k為與降雨量有關(guān)的系數(shù)，它的取值如下：

(9)

2.6 模型中相關(guān)參數(shù)的確定

本研究所需的參數(shù)[20]主要是水資源全球均衡因子γ、水資源全球平均生產(chǎn)能力pw和區(qū)域水資源的產(chǎn)量因子ψ三個(gè)。文章中水資源全球均衡因子取5.19；水資源全球平均生產(chǎn)能力取3 140 m3/hm2；河南省水資源產(chǎn)量因子根據(jù)黃林楠[20]的計(jì)算取0.78。

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

本研究中計(jì)算所需的數(shù)據(jù)來(lái)源于2006-2015年的《河南省統(tǒng)計(jì)年鑒》和2005-2014年的《河南省水資源公報(bào)》，如表2所示。

表2 2005-2014年河南省各項(xiàng)數(shù)據(jù)

3 結(jié)果與分析

3.1 河南省水資源生態(tài)足跡分析

對(duì)河南省2005-2014年淡水生態(tài)足跡和水污染生態(tài)足跡進(jìn)行計(jì)算，計(jì)算結(jié)果見表3和圖1。可以看出，相比其他幾種水資源生態(tài)足跡，農(nóng)業(yè)用水生態(tài)足跡占據(jù)絕對(duì)的主導(dǎo)地位，2005-2014年間農(nóng)業(yè)用水生態(tài)足跡平均為21.28 Mhm2，占總水資源生態(tài)足跡的57.61%。其中最大值出現(xiàn)在2012年和2013年，為22.40 Mhm2，最小值出現(xiàn)在2014年為18.63 Mhm2。從總體看，農(nóng)業(yè)用水生態(tài)足跡雖有所波動(dòng)，但穩(wěn)中有升，這主要是由于河南省是一個(gè)農(nóng)業(yè)大省，農(nóng)業(yè)相比其他產(chǎn)業(yè)處于領(lǐng)導(dǎo)地位。其次為工業(yè)用水生態(tài)足跡，平均為8.89 Mhm2，占總水資源生態(tài)足跡的24.08%。其中最大值出現(xiàn)在2013年，為10.05 Mhm2，最小值出現(xiàn)在2005年，為7.56 Mhm2。工業(yè)用水除2008和2014年有所下降之外，總體上升趨勢(shì)比較明顯，這主要得益于近些年河南省經(jīng)濟(jì)的快速發(fā)展。生活用水生態(tài)足跡平均為5.68 Mhm2，占總水資源生態(tài)足跡的15.48%，其中最大值出現(xiàn)在2008年，為6.97 Mhm2，最小值出現(xiàn)在2006年，為4.33 Mhm2。生活用水生態(tài)足跡整體上穩(wěn)中有降，在人口整體增長(zhǎng)的情況下，說(shuō)明人們的節(jié)水意識(shí)有所提高。生態(tài)環(huán)境補(bǔ)水生態(tài)足跡平均為1.05 Mhm2，占總水資源生態(tài)足跡的2.81%，其中最大值出現(xiàn)在2012年，為1.75 Mhm2。伴隨著河南省越來(lái)越嚴(yán)重的環(huán)境污染問(wèn)題，一定要加大生態(tài)環(huán)境補(bǔ)水的力度。水污染生態(tài)足跡平均為0.63 萬(wàn)hm2，2005-2014年間，水污染生態(tài)足跡整體增速明顯，由2005年的0.48 萬(wàn)hm2增加到2014年的0.87 萬(wàn)hm2，增幅達(dá)到1.81倍。

表3 河南省歷年水資源生態(tài)足跡

圖1 河南省2015-2014年生態(tài)足跡趨勢(shì)圖

從圖1可以看出，河南省總水資源生態(tài)足跡歷年都會(huì)有波動(dòng)，但總體呈緩慢上升趨勢(shì)。2005-2014年間總水資源生態(tài)足跡平均為36.92 Mhm2。最小年份出現(xiàn)在2005年為32.70 Mhm2；最大年份出現(xiàn)在2012年為39.45 Mhm2，增加了20.64%。而淡水生態(tài)足跡和總水資源生態(tài)足跡的變化趨勢(shì)一致，這主要是由于淡水生態(tài)足跡在總水資源生態(tài)足跡中占絕大部分。

3.2 河南省水資源各類生態(tài)指數(shù)分析

為了進(jìn)一步研究河南省水資源利用概況，計(jì)算了河南省人均水資源生態(tài)足跡、水資源生態(tài)承載力、水資源生態(tài)盈余(赤字)、萬(wàn)元GDP水資源生態(tài)足跡和水資源負(fù)載指數(shù)，具體見表4。

表4 河南省歷年水資源指數(shù)

在不考慮黃河客水情況下，2005—2014年間河南省人均水資源生態(tài)足跡相對(duì)較低，在0.31～0.36 hm2波動(dòng)，最大的0.36 hm2出現(xiàn)在2006年、2008年和2009年，最小的0.31 hm2出現(xiàn)在2005年和2014年；水資源生態(tài)承載力用來(lái)描述水資源對(duì)生態(tài)系統(tǒng)和經(jīng)濟(jì)系統(tǒng)健康可持續(xù)發(fā)展的支撐能力，河南省水資源生態(tài)承載力呈現(xiàn)波浪形，分別在2007年、2010年和2014年有所上升，其余年份相對(duì)下降，2005年最大為27.68 Mhm2，2013年最小為9.11 Mhm2；10年間河南省水資源生態(tài)赤字波動(dòng)較大，最小年份出現(xiàn)在2005年為-1.42 Mhm2，最大年份出現(xiàn)在2013年為-23.58 Mhm2。

在考慮黃河客水情況下，2005-2014年間，河南省人均水資源生態(tài)足跡保持在0.35～0.42 hm2之間。各年份之間變化不是特別大，和總水資源生態(tài)足跡變化趨勢(shì)相似，人均水資源生態(tài)足跡在2005年最小為0.35 hm2，在2012年達(dá)到最大為0.42 hm2。

10年間河南省水資源生態(tài)承載力變化趨勢(shì)和不考慮黃河客水情況相似，2005年最大為28.80 Mhm2，2013年最小為11.10 Mhm2，下降幅度達(dá)到61.5%。河南省水資源承載力與降雨量之間的相關(guān)系數(shù)為0.925 5，表明水資源承載力的變化與降雨量之間密切相關(guān)。河南省水資源生態(tài)盈余(赤字)在2005-2014年間均為負(fù)值，表明生態(tài)赤字，不利于經(jīng)濟(jì)社會(huì)的良性發(fā)展。和水資源承載力變化趨勢(shì)相似，水資源生態(tài)盈余(赤字)在2005最大為-3.90 Mhm2，2013年最小為-27.96 Mhm2，下降了7.2倍。隨著河南省近幾年經(jīng)濟(jì)的迅速發(fā)展，水資源利用概況持續(xù)惡化，但工業(yè)經(jīng)濟(jì)的發(fā)展不應(yīng)該以生態(tài)為代價(jià)，政府部門應(yīng)該調(diào)整產(chǎn)業(yè)結(jié)構(gòu)，實(shí)現(xiàn)生態(tài)環(huán)境的修復(fù)和保護(hù)。

從萬(wàn)元GDP水資源生態(tài)足跡來(lái)看，2005-2014年間持續(xù)穩(wěn)定下降，由2005年的3.82×10-5hm2下降到2014年的0.99×10-5hm2，下降幅度達(dá)到74.08%。而萬(wàn)元GDP水資源生態(tài)足跡與水資源利用效率是反比關(guān)系，表明河南省水資源利用效率逐年提高。這主要是由于在社會(huì)和生活中，人們普遍意識(shí)到節(jié)約水資源的重要性。河南省水資源生態(tài)足跡和GDP之間的相關(guān)系數(shù)為0.579 2，表明水資源生態(tài)足跡和GDP之間有著一定的密切聯(lián)系。

水資源負(fù)載指數(shù)將區(qū)域的降雨量、人口、GDP和水資源總量聯(lián)系起來(lái)，反應(yīng)了一個(gè)地區(qū)的水資源利用程度。由表4知，2005-2014年間，河南省水資源負(fù)載指數(shù)增加明顯，2013年是2005年的4.97倍。即便是水資源負(fù)載指數(shù)最小的2005年也屬于Ⅰ類。表明河南省水資源利用程度很高，水資源前景不容樂(lè)觀，有條件時(shí)急需外流域的調(diào)水，魯豐先[22]的研究中也得到相似的結(jié)論。

4 結(jié)論與建議

生態(tài)足跡主要用來(lái)評(píng)價(jià)一個(gè)地區(qū)的可持續(xù)利用狀態(tài)和生態(tài)壓力。河南省2005-2014年間，水資源生態(tài)足跡整體上呈不斷上升趨勢(shì)，對(duì)水資源的需求量也不斷加大。其中對(duì)總水資源生態(tài)足跡的貢獻(xiàn)大小依次為農(nóng)業(yè)用水、工業(yè)用水、生活用水、生態(tài)環(huán)境補(bǔ)水和“灰色水”，這也與河南省的發(fā)展現(xiàn)狀相符合，為了實(shí)現(xiàn)社會(huì)經(jīng)濟(jì)的可持續(xù)發(fā)展，應(yīng)該提高水資源特別是農(nóng)業(yè)用水的利用效率和生態(tài)環(huán)境補(bǔ)水的力度。

從總體來(lái)看，河南省水資源比較緊缺，2005-2014年間都出現(xiàn)不同程度的生態(tài)赤字，水資源負(fù)載指數(shù)在10年間也全為Ⅰ ，并逐年增大，水資源的利用前景不容樂(lè)觀。人均水資源生態(tài)足跡相對(duì)較低，在考慮黃河客水情況下一直在0.35～0.42 hm2之間波動(dòng)，但萬(wàn)元GDP水資源生態(tài)足跡逐年下降，表明河南省水資源的利用效率不斷提高。河南省水資源主要來(lái)自于降雨，水資源承載力和降雨量之間關(guān)聯(lián)性比較強(qiáng)。為了社會(huì)經(jīng)濟(jì)的良性發(fā)展，對(duì)于水利工程應(yīng)該加大建設(shè)，不斷提高對(duì)干旱年的應(yīng)對(duì)能力。

[1] 凌子燕,劉 銳.基于主成分分析法的廣東省區(qū)域水資源緊缺風(fēng)險(xiǎn)評(píng)價(jià)[J].資源科學(xué),2010,32(12):2 324-2 328.

[2] 高媛媛,王紅瑞,許新宜,等.水資源安全評(píng)價(jià)模型構(gòu)建與應(yīng)用----以福建省泉州市為例[J].自然資源學(xué)報(bào),2012,27(2):204-214.

[3] 萬(wàn)坤揚(yáng),胡其昌. 基于層次分析法的杭州市水資源安全現(xiàn)狀評(píng)價(jià)及趨勢(shì)[J]. 水電能源科學(xué),2013,(1):21-25,222.

[4] 韓 美,杜 煥,張 翠,等. 黃河三角洲水資源可持續(xù)利用評(píng)價(jià)與預(yù)測(cè)[J]. 中國(guó)人口·資源與環(huán)境,2015,(7):154-160.

[5] 王 壬,陳 瑩,陳興偉. 區(qū)域水資源可持續(xù)利用評(píng)價(jià)指標(biāo)體系構(gòu)建[J]. 自然資源學(xué)報(bào),2014,(8):1 441-1 452.

[6] 向文英,楊 靜,張 雪. 模糊綜合評(píng)價(jià)法的改進(jìn)及其在水庫(kù)水質(zhì)評(píng)價(jià)中的應(yīng)用[J]. 安全與環(huán)境學(xué)報(bào),2015,(6):344-348.

[7] 欽麗娟,曹劍峰,平建華,等. 模糊數(shù)學(xué)在鄭州市水資源價(jià)值評(píng)價(jià)中的應(yīng)用[J]. 吉林大學(xué)學(xué)報(bào)(地球科學(xué)版),2005,(4):487-490,495.

[8] 陳守煜,李慶國(guó). 一種新的模糊聚類神經(jīng)網(wǎng)絡(luò)及其在水資源評(píng)價(jià)中的應(yīng)用[J]. 水利學(xué)報(bào),2005,(6):662-666.

[9] 王文國(guó),何明雄,潘 科,等. 四川省水資源生態(tài)足跡與生態(tài)承載力的時(shí)空分析[J]. 自然資源學(xué)報(bào),2011,(9):1 555-1 565.

[10] 劉 梅,許新宜,王紅瑞,等. 基于虛擬水理論的河北省水足跡時(shí)空差異分析[J]. 自然資源學(xué)報(bào),2012,(6):1 022-1 034.

[11] 馬 晶,彭 建. 水足跡研究進(jìn)展[J]. 生態(tài)學(xué)報(bào),2013,3(18):5 458-5 466.

[12] 張義,張合平,李豐生,等. 基于改進(jìn)模型的廣西水資源生態(tài)足跡動(dòng)態(tài)分析[J]. 資源科學(xué),2013,(8):1 601-1 610.

[13] 楊建軍,洪 輝,付 娜,等. 水資源生態(tài)足跡消費(fèi)賬戶及其計(jì)算模型----以西安市為例[J]. 安全與環(huán)境學(xué)報(bào),2010,(1):122-126.

[14] 譚秀娟,鄭欽玉. 我國(guó)水資源生態(tài)足跡分析與預(yù)測(cè)[J]. 生態(tài)學(xué)報(bào),2009,(7):3 559-3 568.

[15] 李玉平,王曉妍,朱 琛,等.邢臺(tái)市水資源生態(tài)足跡核算與預(yù)測(cè)研究[J].水土保持研究,2014,(3):227-230.

[16] 張建龍,馮慧敏. 基于生態(tài)足跡法的山西省水資源承載能力研究[J]. 黑龍江大學(xué)工程學(xué)報(bào),2013,(1):61-65.

[17] 陳棟為,陳曉宏,孔 蘭. 基于生態(tài)足跡法的區(qū)域水資源生態(tài)承載力計(jì)算與評(píng)價(jià)----以珠海市為例[J]. 生態(tài)環(huán)境學(xué)報(bào),2009,(6):2 224-2 229.

[18] 孫宇飛,王建平,王曉娟. 關(guān)于“三條紅線”指標(biāo)體系的幾點(diǎn)思考[J]. 水利發(fā)展研究,2010,(10):62-65,72.

[19] 于 強(qiáng),王金龍,王亞南. 基于水資源足跡分析的河北省城鎮(zhèn)化發(fā)展路徑[J]. 經(jīng)濟(jì)地理,2014,(11):69-73,106.

[20] 黃林楠,張偉新,姜翠玲,等. 水資源生態(tài)足跡計(jì)算方法[J]. 生態(tài)學(xué)報(bào),2008,(3):1 279-1 286.

[21] 孟麗紅,葉志平,袁素芬,等. 江西省2007-2011年水資源生態(tài)足跡和生態(tài)承載力動(dòng)態(tài)特征[J]. 水土保持通報(bào),2015,(1):256-261.

[22] 魯豐先. 河南省綜合生態(tài)承載力研究[D]. 河南開封：河南大學(xué),2009.