莫崇勛, 趙梳坍, 阮俞理, 莫細(xì)喜, 孫桂凱, 黃 亞

(1.廣西大學(xué) 土木建筑工程學(xué)院, 廣西 南寧530004; 2.廣西大學(xué) 工程防災(zāi)與結(jié)構(gòu)安全教育部重點(diǎn)實(shí)驗(yàn)室, 廣西 南寧530004; 3.廣西大學(xué) 廣西防災(zāi)減災(zāi)與工程安全重點(diǎn)實(shí)驗(yàn)室,廣西 南寧530004; 4.中國(guó)水利水電科學(xué)研究院 流域水循環(huán)模擬與調(diào)控國(guó)家重點(diǎn)實(shí)驗(yàn)室, 北京 100038)

水資源生態(tài)特征時(shí)空變化規(guī)律及其驅(qū)動(dòng)因素是當(dāng)前水文水資源學(xué)科領(lǐng)域關(guān)注的熱點(diǎn)之一[1]。生態(tài)足跡模型因具有概念直觀、操縱性強(qiáng)、區(qū)域可比性強(qiáng)等優(yōu)點(diǎn)[2-3]而被廣泛運(yùn)用于水資源生態(tài)特征時(shí)空變化規(guī)律研究，可用于客觀地揭示自然資本與經(jīng)濟(jì)發(fā)展間的關(guān)系，如張洺也和趙春子[4]采用該模型將水資源劃分為生活、生產(chǎn)、生態(tài)子賬戶，并結(jié)合壓力指數(shù)、泰爾指數(shù)分析了延邊州的水資源生態(tài)環(huán)境；WANG等[5]采用結(jié)合了均勻均衡和收益因子的改進(jìn)生態(tài)足跡模型，分析了湖北省的水資源可持續(xù)利用狀況；Li H.等[6]采用結(jié)合土壤水分的改進(jìn)生態(tài)足跡法，分析了黃河下游的水資源綜合利用情況。已有研究側(cè)重于從水資源生態(tài)足跡的賬戶劃分及預(yù)測(cè)出發(fā)，對(duì)其驅(qū)動(dòng)因素的研究相對(duì)較少。同時(shí)，目前用于定量分析資源利用驅(qū)動(dòng)機(jī)制的主要方法有結(jié)構(gòu)分解法(SDA)和指數(shù)分解法(IDA)[7]，其中IDA中的迪氏指數(shù)分解模型(LMDI)具有更好消除分解因素殘差[8]，達(dá)到降低結(jié)果誤差的優(yōu)點(diǎn)，已被廣泛應(yīng)用于研究能源和水資源的驅(qū)動(dòng)因素[9-11]。上述研究主要集中在干旱[1,12]、半干旱[4,13]及經(jīng)濟(jì)發(fā)達(dá)地區(qū)[14-15]，而對(duì)于水資源豐富、經(jīng)濟(jì)落后地區(qū)的研究相對(duì)較少。廣西壯族自治區(qū)地處我國(guó)西南邊疆，降水豐沛的同時(shí)面臨著水資源時(shí)空分布不均且局部地區(qū)干旱和洪澇日益加劇、區(qū)內(nèi)生態(tài)—經(jīng)濟(jì)用水不協(xié)調(diào)的問(wèn)題[16]。以往針對(duì)廣西水資源生態(tài)特征的研究，側(cè)重于全區(qū)[17]或某一城市[18-19]在時(shí)間序列上的變化，而鮮少涉及到地區(qū)差異性及驅(qū)動(dòng)因素的分析。因此，論文基于生態(tài)足跡模型分析廣西地區(qū)水資源生態(tài)足跡、水資源生態(tài)承載力、水資源生態(tài)容量的時(shí)空變化規(guī)律，在此基礎(chǔ)上采用LMDI模型分析引起生態(tài)足跡變化的驅(qū)動(dòng)因素及各地級(jí)市在關(guān)鍵驅(qū)動(dòng)因素上的貢獻(xiàn)比例，以期為廣西或同類區(qū)域水資源合理開(kāi)發(fā)利用與管理保護(hù)提供參考依據(jù)。

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

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

(1) 水資源生態(tài)足跡(EFW)為水資源生態(tài)足跡模型中用于衡量水資源消耗狀況的指標(biāo)，是指某一區(qū)域在某一時(shí)間段內(nèi)，人類為維系正常生產(chǎn)生活所需占用的水資源生態(tài)生產(chǎn)性用地面積，可用于反映現(xiàn)有水資源的可持續(xù)性及其對(duì)生態(tài)系統(tǒng)的施壓狀況。

EFW=N·efW=N·γW·〔W/PW〕

(1)

式中：EFW水資源總生態(tài)足跡(108hm2);N為人口數(shù)(104人); efw為人均水資源生態(tài)足跡(hm2/人);γw為水資源的全球均衡因子;W為人均水資源消耗量(m3/人);PW為水資源全球平均生產(chǎn)能力(m3/hm2)。根據(jù)黃林楠和張偉新[20]基于WWF 2000對(duì)全球范圍水資源生態(tài)足跡模型內(nèi)關(guān)鍵參數(shù)的核算結(jié)果,γw取5.19,PW取3 140 m3/hm2。

(2) 水資源生態(tài)承載力(ECW)為水資源生態(tài)足跡模型中用于衡量水資源供給狀況的指標(biāo)，是指某一區(qū)域在某一時(shí)間段內(nèi)水資源最大供給量所能維系該區(qū)域資源、環(huán)境及社會(huì)可持續(xù)發(fā)展的能力。引入了產(chǎn)量因子的概念，使水資源生態(tài)生產(chǎn)性面積具備區(qū)域可比性，可用于反映不同區(qū)域間現(xiàn)有水資源對(duì)生態(tài)系統(tǒng)的支撐狀況。計(jì)算公式為：

ECW=N·ecW=0.4×Ψ·〔Q/PW〕

(2)

式中：ECW為水資源生態(tài)承載力(108hm2);N為人口數(shù)(104人); ecW為人均水資源生態(tài)承載力(hm2/人);Ψ為區(qū)域水資源產(chǎn)量因子;Q為區(qū)域水資源總量(108m3);PW為水資源全球平均生產(chǎn)能力(m3/hm2)。

(3) 水資源生態(tài)容量(Budget)為水資源生態(tài)足跡模型中用于衡量水資源供需狀況的指標(biāo)，由水資源生態(tài)承載力與水資源生態(tài)足跡兩者間的差值直接表示，可用于反映現(xiàn)有水資源供需的客觀狀況及區(qū)域的生產(chǎn)生活是否屬于該區(qū)域生態(tài)系統(tǒng)可承載的范圍。計(jì)算如下：

Budget=ECW-EFW=N×(ecW-efW)

(3)

式中：Budget為水資源生態(tài)容量(108hm2),根據(jù)與0的大小關(guān)系劃分水資源生態(tài)處境為水資源生態(tài)盈余、水資源生態(tài)平衡或水資源生態(tài)赤字[18]。

1.2 LMDI模型

借鑒Kaya恒等式在能源因素分解中的應(yīng)用，運(yùn)用迪氏指數(shù)分解法(LMDI)[21]構(gòu)建廣西地區(qū)水資源生態(tài)足跡LMDI分解模型，具體模型為[13]：

(4)

式中：EFW為廣西水資源生態(tài)總足跡(108hm2);i為水資源生態(tài)足跡賬戶類型; EFWij為j地區(qū)i賬戶的水資源生態(tài)足跡(108hm2);Tj為j地區(qū)國(guó)民生產(chǎn)總值(108元);Pj為j地區(qū)年人口總數(shù)(104人),代表人口發(fā)展水平;Sij為j地區(qū)第i類用水生態(tài)足跡對(duì)總用水生態(tài)足跡的占比,代表j地區(qū)用水結(jié)構(gòu);Ij為j地區(qū)單位國(guó)民生產(chǎn)總值的水資源生態(tài)足跡(hm2/元),代表j地區(qū)的技術(shù)發(fā)展水平;Gj為j地區(qū)人均國(guó)民生產(chǎn)總值(104元/人),代表j地區(qū)的經(jīng)濟(jì)發(fā)展水平。

(5)

(6)

(7)

(8)

(9)

式中：ΔEFW為水資源生態(tài)足跡在從基準(zhǔn)年到目標(biāo)年 間的變化量(108hm2/a); ΔS,ΔI,ΔG,ΔP分別為結(jié)構(gòu)效應(yīng)、技術(shù)效應(yīng)、經(jīng)濟(jì)效應(yīng)和人口效應(yīng)(108hm2/a)，是指水資源生態(tài)足跡中由于結(jié)構(gòu)因素、技術(shù)因素、經(jīng)濟(jì)因素和人口因素引起的水資源生態(tài)足跡變化量,通過(guò)LMDI加和分解即公式(6)—(7)進(jìn)行計(jì)算。

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

研究數(shù)據(jù)包括廣西各地級(jí)市人口數(shù)、國(guó)民生產(chǎn)總值、水資源總量及用水賬戶明細(xì)，人口數(shù)和國(guó)民生產(chǎn)總值來(lái)源于《廣西統(tǒng)計(jì)年鑒》(2008—2017年)；各地級(jí)市的水資源總量及用水?dāng)?shù)據(jù)來(lái)源于《廣西水資源公報(bào)》(2008—2017年)。其中，國(guó)民生產(chǎn)總值以2008年為基期進(jìn)行實(shí)際國(guó)民生產(chǎn)總值的計(jì)算，以消除價(jià)格波動(dòng)對(duì)計(jì)算的誤差影響。

2 結(jié)果與分析

2.1 水資源生態(tài)特征時(shí)間變化規(guī)律

2008—2017年廣西地區(qū)的水資源生態(tài)特征表現(xiàn)為水資源生態(tài)足跡整體下降，水資源生態(tài)承載力和水資源生態(tài)容量歷年波動(dòng)整體上升(圖1)，對(duì)歷年結(jié)果做線性趨勢(shì)擬合可知：人均水資源生態(tài)足跡的平均下降速率為-0.017 4 hm2/人，其中，2013—2017年呈現(xiàn)較大的下降幅度;人均水資源生態(tài)承載力的總上升量達(dá)0.116 hm2/人，平均上升速率為0.325 hm2/(人· a)，最高值出現(xiàn)在2008和2015年，分別達(dá)到9.139，9.513 hm2/人，最小值出現(xiàn)在2011年，為3.109 hm2/人，最值間相差6.404；人均水資源生態(tài)容量的平均上升速率為0.342 1 hm2/(人· a)，最高值出現(xiàn)在2008年和2015年，達(dá)到8.124和8.617 hm2/人，最小值出現(xiàn)在2011年，為2.146 hm2/人，最值間相差6.41。

圖1 廣西地區(qū)2008-2017年人均水資源生態(tài)特征

通過(guò)水資源生態(tài)容量直觀地反映各項(xiàng)水資源供需狀況及其與生態(tài)系統(tǒng)的關(guān)系(表1)，全區(qū)人均水資源生態(tài)容量值在2008，2012，2017年分別為8.124，6.409和8.182，處于生態(tài)盈余狀態(tài)且呈微弱增加趨勢(shì)；人均生態(tài)容量在時(shí)間上變化不顯著，中間年份的變異系數(shù)值較始末年份大，中間年份的最大值與極差的排序?yàn)橹?；地?jí)市間人均水資源生態(tài)容量的相對(duì)差異度隨時(shí)間無(wú)顯著變化，但在空間上存在較大差異。

表1 地級(jí)市尺度人均水資源生態(tài)容量的主要統(tǒng)計(jì)量

2.2 水資源生態(tài)特征空間變化規(guī)律

由圖2可以看出，2008—2017年廣西地區(qū)的水資源生態(tài)特征在空間上呈現(xiàn)一定的差異性和不均衡性。核算期內(nèi)廣西地區(qū)的年均水資源生態(tài)足跡為0.938 hm2/人，有7個(gè)地級(jí)市的人均水資源生態(tài)足跡高于全區(qū)均值；其中，來(lái)賓市以1.566 hm2/人居全區(qū)之首，除此外，桂林市、賀州市、防城港市、北海市的年均水資源生態(tài)足跡達(dá)1.384，1.123，1.130，1.114 hm2/人，超過(guò)1.100，屬于人均用水消耗較大的地級(jí)市；年均水資源生態(tài)足跡最高的來(lái)賓市和最低的玉林市(0.613 hm2/人)在數(shù)值上相差2.55倍。廣西的年均水資源生態(tài)承載力為6.875 hm2/人，有7個(gè)地級(jí)市的人均水資源生態(tài)承載力高于全區(qū)均值；其中，防城港市以27.105 hm2/人居全區(qū)之首，除此外，桂林市、柳州市、賀州市的年均水資源生態(tài)承載力達(dá)18.127，14.061，11.451 hm2/人，超過(guò)10，屬于水資源支撐社會(huì)可持續(xù)發(fā)展能力較大的地級(jí)市；年均水資源生態(tài)承載力最高的防城港市和最低的貴港市(2.793 hm2/人)在數(shù)值上相差9.70倍。廣西各地級(jí)市的年均水資源生態(tài)容量均大于0，生態(tài)處境上屬于水資源生態(tài)盈余。其中，柳州市、桂林市、防城港市、賀州市、河池市的年均生態(tài)容量高于全區(qū)均值(5.937 hm2/人)；年均水資源生態(tài)容量最高值出現(xiàn)在防城港市，達(dá)25.975 hm2/人，最小值出現(xiàn)在南寧市，為1.310 hm2/人，二者在水資源供給量反向相差1.50倍，年均GDP值反向相差5.54倍，年均人口密度反向相差2.23倍的情況下，水資源生態(tài)容量數(shù)值上相差19.08倍；反映了人均水資源生態(tài)容量在人口密度大且經(jīng)濟(jì)發(fā)展快的水資源相對(duì)匱乏城市，較易出現(xiàn)水資源供給矛盾，超過(guò)區(qū)域生態(tài)系統(tǒng)承載的范圍。

2.3 水資源生態(tài)足跡驅(qū)動(dòng)因素分析

由圖3可以看出，廣西地區(qū)水資源生態(tài)足跡驅(qū)動(dòng)因素的分解結(jié)果表現(xiàn)為總效應(yīng)、結(jié)構(gòu)效應(yīng)歷年波動(dòng)整體下降，技術(shù)效應(yīng)整體下降，經(jīng)濟(jì)效應(yīng)、人口效應(yīng)整體上升。其中，經(jīng)濟(jì)效應(yīng)和人口效應(yīng)的分解值呈正數(shù)，正向拉動(dòng)全區(qū)水資源生態(tài)足跡的增長(zhǎng)，而結(jié)構(gòu)效應(yīng)和技術(shù)效應(yīng)的分解值呈負(fù)數(shù)，反向抑制全區(qū)水資源生態(tài)足跡的增長(zhǎng)。對(duì)歷年結(jié)果做線性趨勢(shì)擬合可知：總效應(yīng)的下降速率為4.40×106hm2/a，時(shí)間變化上在2008—2011年呈正值緩慢上升，在2011—2014年由正轉(zhuǎn)負(fù)大幅下降，而后呈負(fù)值又較快上升；結(jié)構(gòu)效應(yīng)的下降速率為4.10×106hm2/a，核算期內(nèi)的變化趨勢(shì)與總效應(yīng)一致；技術(shù)效應(yīng)的下降速率為5.04×106hm2/a，在核算期內(nèi)始終呈負(fù)值；經(jīng)濟(jì)效應(yīng)的上升速率為4.13×106hm2/a，在核算期內(nèi)始終呈正值；人口效應(yīng)的上升速率相對(duì)較小，僅為6.10×105hm2/a，可忽略不計(jì)。在正向拉動(dòng)水資源生態(tài)足跡的效應(yīng)中，經(jīng)濟(jì)效應(yīng)的上升速率最大，在反向抑制水資源生態(tài)足跡的效應(yīng)中，各效應(yīng)下降速率的排序?yàn)椋杭夹g(shù)效應(yīng)>結(jié)構(gòu)效應(yīng)>人口效應(yīng)。

圖2 廣西地區(qū)2008-2017年各地級(jí)市人均水資源生態(tài)特征

圖3 廣西水資源生態(tài)足跡驅(qū)動(dòng)因素分解效應(yīng)

通過(guò)累計(jì)變化量和累計(jì)貢獻(xiàn)率進(jìn)一步確定水資源生態(tài)足跡驅(qū)動(dòng)因素中的關(guān)鍵驅(qū)動(dòng)因素，為結(jié)構(gòu)效應(yīng)、技術(shù)效應(yīng)和經(jīng)濟(jì)效應(yīng)(表2)。核算期內(nèi)，經(jīng)濟(jì)效應(yīng)引起的水資源生態(tài)足跡累計(jì)增加3.78×107hm2，累計(jì)貢獻(xiàn)率為191.8%，是拉動(dòng)廣西水資源生態(tài)足跡的第一因素；人口因素引起的水資源生態(tài)足跡累計(jì)增加5.00×106hm2，是拉動(dòng)廣西水資源生態(tài)足跡的第二因素，對(duì)于水資源生態(tài)足跡變化的作用不顯著；技術(shù)效應(yīng)引起的水資源生態(tài)足跡累計(jì)減少4.69×107hm2，累計(jì)貢獻(xiàn)率為-238.0%，是抑制廣西水資源生態(tài)足跡增長(zhǎng)的第一因素；結(jié)構(gòu)效應(yīng)引起的水資源生態(tài)足跡累計(jì)減少-1.50×107hm2，累計(jì)貢獻(xiàn)率為-79.2%，是作為抑制水資源消耗的第二因素。

表2 廣西水資源生態(tài)足跡驅(qū)動(dòng)因素值

2.4 各地級(jí)市在關(guān)鍵驅(qū)動(dòng)因素上的貢獻(xiàn)率

圖4為統(tǒng)計(jì)對(duì)廣西水資源生態(tài)足跡關(guān)鍵驅(qū)動(dòng)因素，即結(jié)構(gòu)效應(yīng)、技術(shù)效應(yīng)和經(jīng)濟(jì)效應(yīng)貢獻(xiàn)最大的5個(gè)地級(jí)市。由圖4可知，結(jié)構(gòu)效應(yīng)的增加效應(yīng)主要分布在百色市、來(lái)賓市(2014—2017年)、欽州市，年均貢獻(xiàn)比例為125.3%，12.4%，12.4%。其中，百色市和欽州市的增加效應(yīng)整體呈減少趨勢(shì)，來(lái)賓市的增加效應(yīng)呈穩(wěn)定趨勢(shì)；減少效應(yīng)主要分布在南寧市、桂林市(2012—2017年)，平均貢獻(xiàn)比例為-35%，-25%；技術(shù)效應(yīng)的增加效應(yīng)主要分布在梧州市，年均貢獻(xiàn)比例為19.3%，減少效應(yīng)主要分布在南寧市、柳州市、桂林市、玉林市，年均貢獻(xiàn)比例為-17.1%，-10.6%，-27.2%，-12.1%，整體呈現(xiàn)減少效應(yīng)；經(jīng)濟(jì)效應(yīng)的增加效應(yīng)主要分布在南寧市、柳州市、桂林市、貴港市，年均貢獻(xiàn)比例為18.5%，10.9%，20.5%，11.7%。全區(qū)各地級(jí)市的經(jīng)濟(jì)效應(yīng)整體呈增長(zhǎng)趨勢(shì)，減少效應(yīng)集中表現(xiàn)在梧州市，年均貢獻(xiàn)率為-22.2%，且在2009—2011年后梧州市的較少效應(yīng)呈穩(wěn)定變小。

3 結(jié) 論

(1) 廣西2008—2017年人均水資源生態(tài)足跡整體呈下降趨勢(shì)，人均水資源生態(tài)承載力整體呈波動(dòng)上升趨勢(shì)，人均水資源生態(tài)容量均處于生態(tài)盈余狀態(tài)且整體呈微弱上升趨勢(shì)，平均上升速率0.342 hm2/(人·a)，全區(qū)年均水資源生態(tài)盈余為5.94×108hm2/人；區(qū)內(nèi)水資源生態(tài)容量分布不均，年均水資源生態(tài)盈余最大的防城港市與最小的南寧市在數(shù)值上相差19.08倍。

圖4 廣西各地級(jí)市在關(guān)鍵驅(qū)動(dòng)因素上的貢獻(xiàn)比例

(2) 經(jīng)濟(jì)效應(yīng)、人口效應(yīng)是拉動(dòng)廣西地區(qū)水資源生態(tài)足跡增長(zhǎng)的第一、第二因素，核算期內(nèi)對(duì)水資源生態(tài)足跡的貢獻(xiàn)效率分別達(dá)191.80%和25.38%，人口效應(yīng)的作用不顯著；技術(shù)效應(yīng)、結(jié)構(gòu)效應(yīng)是抑制水資源生態(tài)足跡增長(zhǎng)的第一、第二因素，核算期內(nèi)對(duì)水資源生態(tài)足跡的貢獻(xiàn)效率達(dá)-238.0%，-79.18%；反映了廣西地區(qū)水資源利用效率的不斷提高，該區(qū)正處于經(jīng)濟(jì)優(yōu)化產(chǎn)業(yè)結(jié)構(gòu)、發(fā)展綠色經(jīng)濟(jì)的關(guān)鍵節(jié)點(diǎn)。

(3) 水資源生態(tài)足跡驅(qū)動(dòng)因素中的關(guān)鍵驅(qū)動(dòng)因素為結(jié)構(gòu)效應(yīng)、技術(shù)效應(yīng)和經(jīng)濟(jì)效應(yīng)。南寧市、桂林市在結(jié)構(gòu)效應(yīng)、技術(shù)效應(yīng)中表現(xiàn)出較大的減少效應(yīng)，在結(jié)構(gòu)效應(yīng)中對(duì)水資源生態(tài)足跡增長(zhǎng)的平均貢獻(xiàn)比例分別為-35%和-25%，在技術(shù)效應(yīng)中對(duì)水資源生態(tài)足跡增長(zhǎng)的平均貢獻(xiàn)比例為-17.1%，-10.6%。除梧州市外，全區(qū)各地級(jí)市在技術(shù)效應(yīng)上呈減少效應(yīng)。全區(qū)城市在發(fā)展與規(guī)劃上存在較大差異，建議針對(duì)地區(qū)間不同的生態(tài)—經(jīng)濟(jì)協(xié)調(diào)度因地制宜地采取對(duì)策。