郭衛(wèi)華,周永章

基于PLSR的中山市水資源壓力演變特征與趨勢

郭衛(wèi)華1,2,3,周永章2

(1.中山大學(xué)地理科學(xué)與規(guī)劃學(xué)院,廣東廣州 510275;2.中山大學(xué)地球環(huán)境與地球資源研究中心,廣東廣州 510275; 3.中山市環(huán)境保護局,廣東中山 528400)

以珠江三角洲腹地城市中山市為研究目標,分析城市社會經(jīng)濟發(fā)展多個指標對水資源壓力的影響,并利用偏最小二乘回歸(PLSR)方法對1991年以來中山市水資源壓力的演變特征和趨勢進行分析和預(yù)測。結(jié)果表明,節(jié)水水平和第三產(chǎn)業(yè)發(fā)展是中山市水資源壓力的主要影響因素。1991—2011年間,中山市的水資源壓力整體上呈先降后升的變化趨勢,并將可能在今后20年的時間內(nèi)保持上升趨勢。中山市水資源壓力與人均GDP之間則呈現(xiàn)出N形環(huán)境庫爾茲涅曲線特征,表明隨著經(jīng)濟發(fā)展到一定程度,中山市水資源壓力將有所下降,但同時也面臨著重新上升的壓力。

水資源壓力;偏最小二乘回歸方法;珠江三角洲;環(huán)境庫爾茲涅曲線

水資源壓力是指在一定的時空尺度和自然地理背景下,人類社會經(jīng)濟活動對其區(qū)域內(nèi)水資源和水生態(tài)環(huán)境所產(chǎn)生的影響和沖擊。水資源壓力受區(qū)域人口、社會、經(jīng)濟發(fā)展所驅(qū)動,但又在一定程度上,隨著科學(xué)技術(shù)和城市管理水平的提高而得到緩解。因此,水資源壓力是一個由水資源-社會經(jīng)濟-生態(tài)環(huán)境等方面因素綜合作用的結(jié)果[1-2]。

近年來,國內(nèi)外學(xué)者基于不同的研究區(qū)域,使用不同的評價方法從水資源總量、區(qū)域水資源安全、水資源承載力等不同角度對水資源壓力開展研究,建立了相應(yīng)的指標體系和評價方法[3-5]。但目前關(guān)于區(qū)域水資源壓力的諸多研究中,側(cè)重點在于評價指標體系的構(gòu)建[4];而筆者認為,為深入研究區(qū)域經(jīng)濟發(fā)展過程各影響因素與水資源壓力的內(nèi)在聯(lián)系,需解決以下兩個問題：①如何實現(xiàn)區(qū)域水資源壓力的定量化評價;②如何確定不同影響因素對水資源壓力的貢獻。

筆者以珠江三角洲中山市1991—2011年的相關(guān)數(shù)據(jù)為基礎(chǔ),提出基于城市水資源開發(fā)強度-水環(huán)境污染物相結(jié)合的水資源壓力指標,并利用偏最小二乘回歸(PLSR)方法,研究人口、社會、經(jīng)濟、水資源及環(huán)境管理水平等因素對城市水資源壓力的影響,以期為區(qū)域水資源管理提供參考。

1 研究對象及方法

1.1 研究區(qū)域概況

中山市地處廣東省中南部,珠江口西岸,全市行政區(qū)面積為1 800.14 km2,多年平均水資源總量為17.38億m3。盡管位于雨量充沛的珠三角河網(wǎng)地區(qū),但由于社會經(jīng)濟的迅猛發(fā)展和人口的快速增加,中山市的人均水資源量僅為1100m3/a,低于全國平均水平。此外,中山市水資源年內(nèi)分配不均,汛期(4—10月)降雨量約占全年的85%,且由于地處珠江出?？?每年枯水期都要遭受咸潮上溯侵襲,頻繁出現(xiàn)局部季節(jié)性淡水資源短缺的現(xiàn)象。

1.2 評價指標

1.2.1 水資源壓力指標

在珠三角地區(qū),水資源供需壓力和水質(zhì)性缺水是水資源利用中的兩個主要問題,水資源壓力集中體現(xiàn)在由水量—水質(zhì)所導(dǎo)致的各種問題中。筆者以年度水資源開發(fā)強度和入河污染物等標污染負荷的乘積作為區(qū)域水環(huán)境壓力的表征指標,以實現(xiàn)水環(huán)境壓力的定量化評價。水資源壓力指標y的計算公式為

式中：E為水資源開發(fā)強度(年度總?cè)∷颗c年度水資源總量之比);Pn為入河污染物等標污染負荷,其計算方法見文獻[6],參評項目及標準分別為：COD(10mg/L)、Cr6+(0.05 mg/L)、氰化物(0.1 mg/L)、石油類(0.5 mg/L)、NH3-N(0.5 mg/L)。

1.2.2 水資源壓力的驅(qū)動指標

水資源壓力驅(qū)動指標應(yīng)全面、綜合反映區(qū)域內(nèi)社會經(jīng)濟發(fā)展狀況對水資源系統(tǒng)的影響,因此按照區(qū)域性(能夠體現(xiàn)珠三角地區(qū)城市水資源情勢特征和利用特點)、規(guī)范性(指標科學(xué)、規(guī)范,在統(tǒng)計、水利部門中得到廣泛使用)、可量化(便于定量化,為進一步的數(shù)據(jù)分析提供條件)原則,將以下指標作為水資源壓力的驅(qū)動指標：①社會發(fā)展壓力指標包括：總?cè)丝趚1,農(nóng)林牧漁業(yè)總產(chǎn)值x2,工業(yè)總產(chǎn)值x3,第三產(chǎn)業(yè)生產(chǎn)總值x4;②用水特征指標,即人均用水量x5;③環(huán)境管理指標,即單位產(chǎn)值COD排放量x6。中山市1991—2011年各相關(guān)指標年度統(tǒng)計值見表1。

表1 中山市1991—2011年水資源壓力相關(guān)指標數(shù)據(jù)

1.3 PLSR模型算法

由于城市水資源壓力各指標之間通常存在相關(guān)性,利用傳統(tǒng)的最小二乘法建立多元線性回歸方程常會破壞參數(shù)估計,擴大模型誤差,難以保證預(yù)測模型的精確性和可靠性[7-8],因此本研究中選擇PLSR方法以解決回歸方程中變量的多重共線性問題。PLSR方法是一種基于高維投影思想的非參數(shù)回歸方法,通過提取主成分能很好地概括自變量系統(tǒng),可較好地解釋因變量,并排除系統(tǒng)中的噪聲干擾[9]。近年來,PLSR方法在城市水資源承載能力研究、城市生活用水量預(yù)報等方面得到良好的應(yīng)用[10-11]。PLSR方法的原理是：設(shè)有單因變量y和p個自變量x1、x2、…、xp,觀測n個樣本點,構(gòu)成自變量與因變量的數(shù)據(jù)矩陣X=(xij),其中i=1,2,…,n,j=1,2,…, p和Y=y1,y2,…,yn。PLSR分別在X和Y中提取成分t1和u1,這兩個主成分分別為X和Y的線性組合,它們應(yīng)滿足兩個要求：①t1和u1應(yīng)盡可能大地提取它們各自原變量系統(tǒng)中的變異信息;②t1和u1的相關(guān)程度能夠達到最大。PLSR方法的具體步驟如下[9,11]：

步驟1：標準化處理X和Y,得到的標準化矩陣,分別記為E0和F0。

步驟2：從F0中提取一個成分u1=F0c1,‖c1‖=1。由于F0只是一個變量,所以c1為常數(shù)。同時從E0中提取一個成分t1,t1=E0w1,w1是第一主軸,w1=ET0F0/‖ET0F0‖,并且‖w1‖=1。分別求E0和F0對t1的兩個回歸方程：

其中,回歸系數(shù)向量為

步驟3：用殘差矩陣E1和F1取代E0和F0,用同樣方法求第2個w2軸和c2以及第2個成分u2和t2,如此計算下去,如果進行了m次運算,則會有

由于t1,t2,…,tm均可表示成E01,E02,…,E0p的線性組合,因此,上式可還原成F0關(guān)于E0p的回歸方程形式,即

在PLSR方法中確定主成分個數(shù)的常用方法有舍一交叉驗證方法、分批交叉驗證方法、分裂樣本交叉驗證方法、隨機樣本交叉驗證方法等[9,12]。在本文中采用舍一交叉驗證方法確定主成分的選取個數(shù)。

2 結(jié)果與討論

2.1 水資源壓力驅(qū)動指標多重相關(guān)性診斷

本研究中利用方差膨脹因子對x1～x6的相關(guān)數(shù)據(jù)進行多重相關(guān)性診斷。自變量xj的方差膨脹因子記為SJVIF,其表達式為

式中,R2J為以xj為自變量時對其他自變量回歸的復(fù)測定系數(shù)。

若最大的SJVIF＞10,則認為自變量之間存在高度的相關(guān)現(xiàn)象[9]。本文中水資源壓力驅(qū)動指標的多重相關(guān)性診斷結(jié)果見表2。由表2可見,最大的SJVIF=401＞10,自變量間存在高度的多重相關(guān)性。

表2 水資源壓力驅(qū)動指標的多重相關(guān)性診斷結(jié)果

2.2 PLSR分析

根據(jù)上述原理及步驟,通過統(tǒng)計軟件SAS 9.1 for Windows(SAS Institute Inc.)實現(xiàn)數(shù)據(jù)的PLSR分析。經(jīng)過數(shù)據(jù)標準化、提取成分并作交叉有效性判別,計算結(jié)果提取了2個主成分t1和t2,這兩個主成分解釋了原自變量系統(tǒng)中92.47%的信息,同時解釋了因變量系統(tǒng)中80.89%的變異信息,顯示模型精度較好,其模型效應(yīng)負荷量及標準偏回歸系數(shù)見表3。

表3 模型效應(yīng)負荷量及標準偏回歸系數(shù)

結(jié)果顯示,在提取的兩個主成分中,x1～x4對t1有較大負載,反映了社會經(jīng)濟發(fā)展對水資源造成的壓力;而x5和x6對t2有較大負載,反映了城市節(jié)水水平和環(huán)境管理水平對壓力的影響。而各指標的標準偏最小二乘系數(shù)表明,人均用水量x5對城市水資源壓力的影響最大,顯示節(jié)水水平對于水資源壓力有直接的作用。代表不同產(chǎn)業(yè)發(fā)展特征的x2～x4項指標對水資源壓力的影響基本相同,其中以第三產(chǎn)業(yè)總產(chǎn)值x4的影響最大。而在各項指標中,人口因素和環(huán)境管理因素對水資源壓力的貢獻較小。

通過建立y對提取的主成分的回歸,以及數(shù)據(jù)的逆標準化操作,最終得到的回歸方程為

經(jīng)檢驗,所建方程的決定系數(shù)R2=0.81,p＜0.01,方程具有顯著性。通過多重共線性檢驗得到最大的SJVIF=1.01,多重共線性問題已基本解決。

2.2 水資源壓力趨勢及演變特征

根據(jù)中山市水資源綜合規(guī)劃等資料,獲得2020年及2030年中山市水資源壓力各影響指標的數(shù)據(jù),代入式(10)后得到水資源壓力預(yù)測,結(jié)果見表4?？梢婋S著人口及城市工農(nóng)業(yè)的發(fā)展,今后一段時期內(nèi)中山市水資源壓力仍將保持一定的上升態(tài)勢。2020年中山市水資源壓力為2011年的1.54倍,而到了2030年將達到2.32倍。

表4 中山市水資源壓力預(yù)測

1991—2030年中山市水資源壓力數(shù)據(jù)的演變趨勢表明,在不同的時間段內(nèi),中山市的水資源壓力增長呈現(xiàn)不同的特點。1991—2000年間,中山市水資源壓力在往復(fù)變化中整體呈下降趨勢,年均增長率為-6.89%,主要原因在于該時期內(nèi)環(huán)境保護工作的大力開展,廢水污染物排放得到一定程度的緩解。2001—2010年間,中山市水資源壓力呈迅速上升趨勢,年均增長率為17.3%,增長的主要貢獻來自于工業(yè)總產(chǎn)值以及第三產(chǎn)業(yè)總產(chǎn)值,顯示該時期內(nèi)城鎮(zhèn)化、工業(yè)化過程中所帶來的用水量增長及工業(yè)廢水污染物排放的增加對水資源帶來極大的壓力。而2011—2020年間以及2020—2030年間,由于人口及各項經(jīng)濟指標的上升,水資源壓力不斷增加,但年均增長速度有所下降,分別為4.9% (2011—2020年)和4.2%(2020—2030年)。

2.3 水資源壓力與城市經(jīng)濟發(fā)展的關(guān)系

經(jīng)驗數(shù)據(jù)顯示,發(fā)達國家城市環(huán)境質(zhì)量指標隨經(jīng)濟發(fā)展的變化存在著一個由上升轉(zhuǎn)而下降的轉(zhuǎn)折點。因此,城市環(huán)境質(zhì)量隨收入增長的演變模式可以用環(huán)境庫茲涅茨曲線(environment Kuznets curve, EKC)形式來表示[13-14]。筆者以環(huán)境庫茲涅茨理論為基礎(chǔ),研究中山市人均GDP與水資源壓力之間的關(guān)系。為了消除異方差的問題,并使模型估計更準確,將水資源壓力指標值作自然對數(shù)變換,并采用三次多項式形式建立中山市經(jīng)濟發(fā)展與水資源壓力關(guān)系的EKC分析模型。擬合結(jié)果顯示,方程的一、二、三次項系數(shù)分別為7.548、-1.344、0.0648,R2為0.91,p＜0.01,方程有顯著意義。

圖1顯示,1991—2030年,中山市水資源壓力與人均GDP之間呈現(xiàn)出典型的N形曲線特征。這與國外環(huán)境庫茲涅茨曲線理論一般假說中認為的倒U形曲線形式有所不同[15-16],但與張捷等[17]發(fā)現(xiàn)廣東省人均GDP增長與工業(yè)“三廢”排放量之間的關(guān)系一致。該曲線特征表明,在中山市城市發(fā)展過程中,經(jīng)濟增長和水資源壓力存在復(fù)雜的相關(guān)關(guān)系。當(dāng)社會經(jīng)濟發(fā)展到一定程度后,隨著節(jié)水、環(huán)保意識的提高以及相關(guān)公共管理措施的加強,城市水資源壓力將得到一定的緩解,但并非就從此必然一路下降。隨著人口的增長和社會經(jīng)濟的發(fā)展,水資源壓力的各種內(nèi)在驅(qū)動因素仍然存在,可能導(dǎo)致城市水資源壓力的改善出現(xiàn)曲折反復(fù)。根據(jù)中山市統(tǒng)計資料,2011年中山市人均GDP約為10.5萬元,按照上述EKC曲線,這段時期中山市水資源壓力正處于該曲線二次上升的轉(zhuǎn)折點(人均GDP 9.9萬元)。

圖1 中山市人均GDP與水資源壓力關(guān)系的EKC曲線

3 結(jié) 論

a.利用偏最小二乘回歸方法在自變量存在高度相關(guān)性的條件下進行回歸建模,適用于水資源壓力與城市社會經(jīng)濟發(fā)展各指標間的相互關(guān)系研究。研究結(jié)果顯示節(jié)水水平和第三產(chǎn)業(yè)發(fā)展是中山市水資源壓力的主要影響因素。

b.1991—2011年間,中山市的水資源壓力整體上呈先降后升的演變特征。對中山市水資源壓力的預(yù)測結(jié)果顯示,未來20年中山市的水資源壓力仍然將處于上升的過程。中山市水資源壓力與人均GDP之間呈現(xiàn)出N形EKC曲線特征,顯示隨著城市經(jīng)濟的發(fā)展,水資源壓力有經(jīng)過改善后再次加劇的傾向。在今后的城市發(fā)展過程中,應(yīng)更加重視水資源利用和水環(huán)境保護,通過轉(zhuǎn)變產(chǎn)業(yè)發(fā)展模式、提高水資源利用水平、加大環(huán)境保護工作力度等措施以緩解未來面臨的水資源壓力。

參考文獻：

[1]王穎,邵磊,周孝德,等.山西省水資源系統(tǒng)壓力綜合評價[J].水力發(fā)電學(xué)報,2011,30(6)：189-198.(WANG Ying,SHAO Lei,ZHOU Xiaode,et al.Comprehensive evaluation of water resources system stress in Shanxi Province[J].Journal of Hydroelectric Engineering,2011, 30(6)：189-198.(in Chinese))

[2]郭梅,許振成,彭曉春.水資源安全問題研究綜述[J].水資源保護,2007,23(3)：40-43.(GUO Mei,XU Zhengcheng,PENG Xiaochun.Progress in water security [J].Water Resources Protection,2007,23(3)：40-43.(in Chinese))

[3]賈紹鳳,張軍巖,張士鋒.區(qū)域水資源壓力指數(shù)與水資源安全評價指標體系[J].地理科學(xué)進展,2002,21 (6)：538-545.(JIA Shaofeng,ZHANG Junyan,ZHANG Shifeng.Regionalwaterresourcesstressandwater resources security appraisement indicators[J].Progress in Geography,2002,21(6)：538-545.(in Chinese))

[4]吳佩林.我國區(qū)域發(fā)展的水資源壓力分析[J].西北農(nóng)林科技大學(xué)學(xué)報：自然科學(xué)版,2005,33(10)：143-149. (WU Peilin.A quantificational and comparative analysis on the regional water resource pressure and water scarcity in China[J].Journal of Northwest A&F University：Natural Science Edition,2005,33(10)：143-149.(in Chinese))

[5]文琦,何彤慧.近10年來我國水資源承載力研究綜述[J].水資源保護,2005,21(6)：19-22.(WEN Qi,HE Tonghui.Review on the study of water resources carrying capacity in China in the last ten years[J].Water Resources Protection,2005,21(6)：19-22.(in Chinese))

[6]鐘定勝,張宏偉.等標污染負荷法評價污染源對水環(huán)境的影響[J].中國給水排水,2005,21(5)：101-103. (ZHONG Dingsheng,ZHANG Hongwei.Equal standard pollution load method for evaluating the effect of pollution source on aquatic environment[J].China Water& Wastewater,2005,21(5)：101-103.(in Chinese))

[7]高惠璇.兩個多重相關(guān)變量組的統(tǒng)計分析(2)[J].數(shù)理統(tǒng)計與管理,2002,21(2)：60-64.(GAO Huixuan. Statistical analyses for multiple correlation variables of two sets(2)[J].Application of Statistics and Management, 2002,21(2)：60-64.(in Chinese))

[8]WOLD S,RUHE A,WOLD H,et al.The collinearity problem in linear regression：the partial least squares (PLS)approach to generalized inverses[J].SIAM Journal on Scientific and Statistical Computing,1984,5(3)：735-743.

[9]王慧文,吳載斌,孟潔.偏最小二乘回歸的線性與非線性方法[M].北京：國防工業(yè)出版社,2006：111-116.

[10]王帥,孫月峰.基于耦合逐步回歸的PLS模型城市用水量預(yù)測[J].安全與環(huán)境學(xué)報,2012,12(4)：170-173. (WANG Shuai,SUN Yuefeng.Urban water consumption prediction based on apartial least square model coupled with stepwise regression[J].Journalof Safety and Environment,2012,12(4)：170-173.(in Chinese))

[11]李林,付強.偏最小二乘回歸模型的城市水資源承載能力研究[J].水科學(xué)進展,2005,16(6)：822-825.(LI Lin,FU Qiang.Study on the city water resources capacity model with partial least squares regression[J].Advances in Water Science,2005,16(6)：822-825.(in Chinese))

[12]高惠璇.兩個多重相關(guān)變量組的統(tǒng)計分析(3)：偏最小二乘回歸與PLS過程[J].數(shù)理統(tǒng)計與管理,2002,21 (3)：58-64.(GAO Huixuan.Statistical analyses for multiple correlation variables of two sets(3)：partial leastsquares regression and PLS procedure[J].Application of Statistics and Management,2002,21(3)：58-64.(in Chinese))

[13]陳華文,劉康兵.經(jīng)濟增長與環(huán)境質(zhì)量：關(guān)于環(huán)境庫茲涅茨曲線的經(jīng)驗分析[J].復(fù)旦學(xué)報：社會科學(xué)版,2004 (2)：87-94.(CHEN Huawen,LIU Kangbing.Economic growth and environmental quality：an empirical analysis of the environmentel Kuznets curve[J].Fudang Journal：Social Sciences,2004(2)：87-94.(in Chinese))

[14]DINDA S.Environmental Kuznets curve hypothesis：a survey[J].Ecological Economics,2004,49(4)：431-455.

[15]DUARTE R,PINILLA V,SERRANO A.Is there an environmental Kuznets curve for water use?a panel smooth transition regression approach[J].Economic Modelling,2013,31：518-527.

[16]CHENG H,HU Y.Economic transformation,technological innovation,and policy and institutional reforms hold keys to relieving China's water shortages[J].Environmental Science&Technology,2010,45(2)：360-361.

[17]張捷,張玉媚.廣東省的庫茲涅茨環(huán)境曲線及其決定因素：廣東省工業(yè)化進程中經(jīng)濟增長與環(huán)境變遷關(guān)系的實證研究[J].廣東社會科學(xué),2006(3)：17-23. (Zhang Jie,Zhang Yumei.Environmental Kuznets curve in Guangdong and its determinants：an empirical study of relationship between economic growth and environmental changes in the progress of industrialization in Guangdong Province[J].Social Sciences in Guangdong,2006(3)：17-23.(in Chinese)

Evolution characteristics and trend of water resources stress of Zhongshan City based on PLSR method

GUO Weihua1,2,3,ZHOU Yongzhang2
(1.Geography and Planning School,Sun Yat-sen University,Guangzhou 510275,China; 2.Center for Earth Environment and Resources,Sun Yat-sen University,Guangzhou 510275,China; 3.Zhongshan Environmental Protection Bureau,Zhongshan 528400,China)

Taking Zhongshan City,a hinterland city in the Pearl River Delta,as a case study,the influence of urban socioeconomic development indicators on the water resources stress was examined.The partial least squares regression method was used to analyze the evolution characteristics and predict the trend of the water resources stress of Zhongshan City since 1991.The results show that the water conservation level and the development of the tertiary industry were the principal factors influencing the water resources stress of the city.During the period from 1991 to 2011,the water resources stress of the city first decreased,then increased on the whole,and may maintain an increasing trend over the coming 20 years.The relationship between the water resources stress and per-capita GDP of Zhongshan City was characterized by an N-shaped environmental Kuznets curve,indicating that when the economy develops to a certain extent,the water resources stress of the city may decease,but,meanwhile,rising pressure on water resources in Zhongshan City may still exist.

water resources stress;partial least squares regression method;Pearl River Delta;environmentel Kuznets curve

TV213.4

A

1004-6933(2014)01-0023-05

201307-15 編輯：徐 娟)

10.3969/j.issn.1004-6933.2014.01.005

廣東省科技計劃項目(2009B030801003)

郭衛(wèi)華(1971—),男,博士研究生,研究方向為資源環(huán)境一體化與區(qū)域可持續(xù)發(fā)展。E-mail：guoweihua01@sina.com