陳雯，黃長生，王寧濤

CHEN Wen，HUANG Chang-Sheng，Wang Ning-Tao

（中國地質(zhì)調(diào)查局武漢地質(zhì)調(diào)查中心，武漢 430205）（Wuhan Center of China Geological Survey，Wuhan 430205）

廣州市水資源承載力的主成分分析

陳雯，黃長生，王寧濤

CHEN Wen，HUANG Chang-Sheng，Wang Ning-Tao

（中國地質(zhì)調(diào)查局武漢地質(zhì)調(diào)查中心，武漢 430205）（Wuhan Center of China Geological Survey，Wuhan 430205）

水資源承載力是衡量經(jīng)濟社會可持續(xù)發(fā)展的重要指標(biāo)。本文以廣州市水資源承載力為研究對象，在現(xiàn)有資料的基礎(chǔ)上，運用主成分分析法從14個指標(biāo)中提取出影響水資源承載力動態(tài)變化的兩個主成分，并對廣州市水資源承載力年際變化進(jìn)行了綜合評價。結(jié)果表明：2005-2010年，廣州市水資源承載力呈現(xiàn)逐年下降的趨勢。但就目前情況來看，通過擴展供水工程、加強節(jié)水措施等能緩解供需矛盾，水資源承載力還存在一定的提升空間。

主成分分析；水資源承載力；廣州市

水資源是人類社會賴以生存和發(fā)展的物質(zhì)基礎(chǔ)。近年來，由于水資源不合理開發(fā)利用以及工農(nóng)業(yè)污染造成的水環(huán)境惡化等問題，不僅進(jìn)一步加劇了水資源供需矛盾，也給水資源可持續(xù)利用帶來威脅。因此，水資源承載力作為衡量區(qū)域可持續(xù)發(fā)展的一項重要指標(biāo)，成為當(dāng)前水資源科學(xué)領(lǐng)域研究的重點和熱點。

廣州市地處廣東省中南部，珠江三角洲北緣，全市面積7434.4 km2，下轄越秀區(qū)、荔灣區(qū)、海珠區(qū)、天河區(qū)、白云區(qū)、黃埔區(qū)、番禺區(qū)、花都區(qū)、南沙區(qū)、蘿崗區(qū)10個市轄區(qū)及增城市、從化市2個縣級市。廣州市本地水資源較少，全市多年本地平均水資源總量為79.79億m3，人均水資源占有量僅為1375 m3，不到全國人均水平的1/2。入境水量相對豐富，多年平均徑流量為1141億m3，為本地水資源的14倍[1]。在城市化進(jìn)程不斷推進(jìn)的背景下，經(jīng)濟快速增長和人口劇增使得廣州市用水需求不斷增加，水資源供需矛盾日益突出，已成為制約經(jīng)濟社會可持續(xù)發(fā)展的“瓶頸”因素。因此，研究廣州市水資源承載力，全面分析水資源承載力現(xiàn)狀，對促進(jìn)當(dāng)?shù)厮Y源可持續(xù)利用以及經(jīng)濟社會可持續(xù)發(fā)展具有十分重要的現(xiàn)實意義。

自20世紀(jì)80年代末以來，我國許多學(xué)者在不同流域、區(qū)域、城市開展了水資源承載力研究，已有的研究成果主要集中在以下三個方面[2-14]：（1）對水資源承載力評價指標(biāo)建立的探索；（2）對水資源承載力著重概念、內(nèi)涵等方面的理論研究；（3）對水資源承載力量化方法、評價模型上的探索。存在的主要問題在于，目前水資源承載力研究理論基礎(chǔ)不完善，尚未形成公認(rèn)的基礎(chǔ)理論和研究方法。綜合評價水資源承載力的量化方法有很多，常見的包括常規(guī)趨勢法、系統(tǒng)動力學(xué)法、模糊綜合評價法等，各種方法的應(yīng)用中均存在不足之處[13]。主成分分析法基于降維思想，通過線性變換，將多個指標(biāo)進(jìn)行最佳綜合與簡化，通過少數(shù)幾個主成分來揭示多個指標(biāo)間的內(nèi)部結(jié)構(gòu)關(guān)系，并盡可能多地保留原始指標(biāo)的信息。本文采用主成分分析法，把影響廣州市水資源承載力的多個因子轉(zhuǎn)化為少數(shù)幾個綜合指標(biāo)，并綜合分析廣州市水資源承載力的年際變化趨勢。

1 主成分分析法

本次研究運用SPSS16.0統(tǒng)計分析軟件進(jìn)行數(shù)據(jù)處理。設(shè)水資源承載力評價指標(biāo)體系原來有p個評價指標(biāo)，用X1、X2、…、Xp表示，經(jīng)主成分分析后得到m（m＜p）個主成分Z1、Z2、…、Zp，則主成分分析的主要步驟為：（1）原始數(shù)據(jù)進(jìn)行標(biāo)準(zhǔn)化處理；（2）計算相關(guān)矩陣R；（3）計算R的J個特征根λ1、λ2、…、λJ，以及每個特征根對應(yīng)的特征向量U1、U2、…、UJ；（4）通過計算累計貢獻(xiàn)率，貢獻(xiàn)率E＞85%的原則確定主成分個數(shù)m；（5）求主成分權(quán)重ek；（6）求綜合評價指數(shù)Z。綜合評價指數(shù)Z越大，反映了水資源承載力狀況越強，反之，評價值越小，水資源承載力越弱。

本文采用的數(shù)據(jù)資料來自《廣州市統(tǒng)計年鑒》、《廣東省水資源公報》、《廣州市水資源公報》，選取2005-2010年統(tǒng)計數(shù)據(jù)作為評價的原始數(shù)據(jù)。

2 廣州市水資源承載力分析

2.1 指標(biāo)選取

水資源承載力既受到自然地理條件的約束，又受到區(qū)域人口、社會、經(jīng)濟、生產(chǎn)力水平等多方面的影響。在充分考慮廣州市水資源承載力現(xiàn)狀及各類影響因素的基礎(chǔ)上，指標(biāo)體系的建立主要包含以下四個方面：（1）反映水資源條件；（2）反映水資源供需關(guān)系；（3）反映社會經(jīng)濟發(fā)展?fàn)顩r；（4）反映污水排放處理狀況。同時，也參考了前人的研究工作，基于科學(xué)性、可比性等原則，在2005-2010年基礎(chǔ)數(shù)據(jù)資料中，選取了14個指標(biāo)（表1）對廣州市水資源承載力進(jìn)行動態(tài)綜合評價（表1）。具體包括：

表1 廣州市水資源承載力評價指標(biāo)Table 1 Water resources carrying capacity evaluation of Guangzhou City

GDP X1，億元；固定資產(chǎn)投資額X2，億元；總?cè)丝赬3，萬人；水資源總量X4，億m3；全市入境水量X5，億m3；總供水量X6，億m3；工業(yè)用水量X7，億m3；農(nóng)業(yè)用水量X8，億m3；生活用水量X9，億m3；城鎮(zhèn)居民生活人均日用水量X10，L/d；萬元GDP用水量X11，m3；萬元工業(yè)增加值用水量X12，m3；農(nóng)田灌溉畝均用水量X13，m3；污水處理率X14，%。

2.2 主成分分析

運用SPSS16.0軟件進(jìn)行數(shù)據(jù)分析計算，得出各評價指標(biāo)的標(biāo)準(zhǔn)化矩陣、相關(guān)系數(shù)矩陣、特征值、主成分貢獻(xiàn)率和累計貢獻(xiàn)率（表2）。

表2 評價指標(biāo)體系標(biāo)準(zhǔn)化矩陣Table 2 Standardized evaluation system matrix

從表3可以看出，所選取的評價指標(biāo)如X6和X7、X8存在一定的相關(guān)性，相關(guān)系數(shù)分別為0.953、0.925。指標(biāo)間的相關(guān)性是進(jìn)行主成分分析的基礎(chǔ)，也進(jìn)一步說明了進(jìn)行主成分分析的必要性和科學(xué)性。

表3 評價指標(biāo)體系相關(guān)系數(shù)矩陣Table 3 correlation matrix evaluation system

表4列出了計算之后所有主成分的特征根、主成分貢獻(xiàn)率及累計貢獻(xiàn)率。前兩個主成分的累積貢獻(xiàn)率已達(dá)到88.328%，根據(jù)累積貢獻(xiàn)率＞85%的原則，選擇前兩個主成分進(jìn)行數(shù)據(jù)分析。

表4 主成分特征值及貢獻(xiàn)率Table 4 Principal component eigenvalues and contribution rates

運用SPSS16.0軟件可直接得到主成分荷載矩陣（表5），它是主成分與對應(yīng)變量的相關(guān)系數(shù)。

表5 主成分荷載矩陣Table 5 Principal component loading matrix

從表5分析可知，第一主成分與X1、X2、X3、X14之間存在著強負(fù)相關(guān)關(guān)系，與X6、X7、X8、X10、X11、X12之間存在較強的正相關(guān)關(guān)系；第二主成分與X4、X5之間存在較強的正相關(guān)關(guān)系。這兩個主成分概括了人口、經(jīng)濟社會發(fā)展、水資源利用情況和區(qū)域水資源條件方面的內(nèi)容，因此可以用這些因素的變化來反映廣州市水資源承載力情況：

（1）社會經(jīng)濟發(fā)展與水資源承載力

2005-2010年，廣州市GDP總量翻了一番，經(jīng)濟快速增長伴隨著城市化進(jìn)程不斷加快，導(dǎo)致用水需求不斷增長，但總供水量卻總體上逐年下降，使得現(xiàn)有供水能力與經(jīng)濟增長速度不匹配，對水資源承載力構(gòu)成很大壓力。

（2）人口與水資源承載力

廣州市人口總量從2005年750.53萬增長到2010年1270.96萬，在常住人口持續(xù)增長，外地務(wù)工人員大量涌入的情況下，生活用水量劇增，但實際生活供水能力并未隨著人口的增長幅度相應(yīng)增長。2010年相比2005年，供水能力僅從8.69億m3增長到9.81億m3。同時，人類活動造成的水資源污染不斷加劇，從而引起水資源承載能力降低。

（3）水資源利用情況與水資源承載力

廣州市用水結(jié)構(gòu)以工業(yè)、農(nóng)業(yè)、生活用水為主。近年來，通過調(diào)整產(chǎn)業(yè)結(jié)構(gòu)，逐漸由耗水型向節(jié)水型轉(zhuǎn)變，取得了一定效果，萬元GDP用水量、萬元工業(yè)增加值用水量、農(nóng)田灌溉畝均用水量等指標(biāo)都逐年有所下降（表1），但與全國平均水平相比（2005-2010年中國水資源公報），仍存在一定的節(jié)水空間。

（4）水資源條件與水資源承載力

廣州市入境水量的補給在一定程度上緩解了水資源承載壓力。目前，全市過境河流主要有東江、西江、北江。東江水質(zhì)有局部惡化現(xiàn)象，枯水期不能滿足用水需求。西江、北江水質(zhì)較好。西江開發(fā)利用率較高，且入境水量大，整個流域尚存在較大開發(fā)利用空間；北江全流域開發(fā)利用率小于20%，尚有一定的開發(fā)潛力[15]。

2.3 主成分綜合評價分析

根據(jù)公式（4）、（5）、（6）利用SPSS16.0計算出廣州市水資源承載力主成分單因子得分，再根據(jù)單因子與其特征根貢獻(xiàn)率來確定兩個主成分的綜合得分，由此分析廣州市水資源承載力的變化情況（表6，圖1）。從表6和圖1可以看出，主成分得分有正有負(fù)，正負(fù)值并不代表水資源承載力的真實水平，而是反映其所處的相對位置，負(fù)值表示該年份處于平均水平之下，正值說明該年份處于平均水平之上。綜合指數(shù)Z值越大，說明水資源承載力越大，反之越小?？傮w而言，2005-2010年廣州市水資源承載力基本呈現(xiàn)逐年下降的趨勢，但在2007年出現(xiàn)波動，主要是由于當(dāng)年入境水量相對偏少。

表6 2005-2010年水資源承載力綜合得分Table 6 Water Resources Carrying Capacity composite score in 2005-2010

從圖1可以看出，第一主成分的變化趨勢與水資源承載力變化總體趨勢基本一致。第二主成分則波動較大，經(jīng)歷了先上升再下降再上升的過程，說明廣州市水資源條件具有不穩(wěn)定性。

圖1 主成分趨勢圖Fig.1 Main components of Trend

3 結(jié)論與討論

（1）運用主成分分析法，在保證原始數(shù)據(jù)信息損失最小的情況下，通過線性變換，用少數(shù)幾個具有代表意義的新變量代替原有變量，避免了主觀隨意性，能夠比較客觀全面地反映當(dāng)?shù)厮Y源承載力現(xiàn)狀。

（2）通過主成分分析，確定了人口、社會經(jīng)濟發(fā)展?fàn)顩r、水資源利用情況以及水資源條件為影響廣州市水資源承載力的主要因素。

（3）從單個因子分析，GDP和固定資產(chǎn)投資額的荷載數(shù)值分別為0.984、0.985，可見經(jīng)濟高速發(fā)展對水資源承載力造成了巨大壓力。廣州市作為廣東省的省會城市和經(jīng)濟重鎮(zhèn)，經(jīng)濟一直處于高速發(fā)展?fàn)顟B(tài)，那么在擴展供水工程（地表水、地下水）的同時加強節(jié)水措施，走“集約”型用水道路能很好地延緩供需矛盾加劇問題，促進(jìn)水資源可持續(xù)利用。

（4）本次研究由于資料有限，選取的評價指標(biāo)不夠全面，為了更加全面、客觀地掌握廣州市水資源承載力狀況，還需要進(jìn)一步收集數(shù)據(jù)資料分析研究，尤其是水資源污染相關(guān)數(shù)據(jù)資料，并開展不同評價方法所得結(jié)果的對比分析。

（5）本次研究僅僅對廣州市2005-2010年水資源承載力的年際變化情況進(jìn)行了初步評價，未能深入分析各影響因素對水資源承載力的具體影響，這有待進(jìn)一步研究。

[1]朱小娟,普智曉.廣州市水資源問題及發(fā)展對策研究[J].水資源與水工程學(xué)報,2004,15(4)：67-70.

[2]施雅風(fēng),曲耀光.烏魯木齊河流域水資源承載力及其合理利用[M].北京：科學(xué)出版社,1992：210-220.

[3]許有鵬.干旱地區(qū)水資源承載能力綜合評價[J].自然資源學(xué)報,1993,8(3)：229-237.

[4]阮本青.區(qū)域水資源適度承載能力計算模型研究[J].土壤侵蝕與水土保持學(xué)報,1998,4(3)：57-61.

[5]曲耀光,樊勝岳.黑河流域水資源承載力分析計算與對策[J].中國沙漠,2000,20(1)：1-8.

[6]王在高,梁虹.巖溶地區(qū)水資源承載力指標(biāo)體系及其理論模型初探[J].中國巖溶,2001,20(2)：143-148.

[7]秦莉云.淮河流域水資源承載能力的評價分析[J].河南農(nóng)業(yè)大學(xué)學(xué)報,2001(35)：98-100.

[8]魏斌.城市水資源合理利用與水資源承載力研究[J].干旱區(qū)地理,2002,25(4)：377-381.

[9]吳九紅,曾開華.城市水資源承載力的系統(tǒng)動力學(xué)研究[J].水利經(jīng)濟,2003,21(3)：36-39.

[10]惠泱河,蔣曉輝,黃強,薛小杰.水資源承載力評價指標(biāo)體系研究[J].水土保持通報,2001,20(6)：30-34.

[11]陳洋波,李長興,馮智瑤,陳俊合.深圳市水資源承載能力模糊綜合評價[J].水力發(fā)電,2004,30(3)：10-14.

[12]杜敏,徐征和,彭利民,朱運海,徐秀鳳.濟寧市水資源承載力評價與分析[J].濟南大學(xué)學(xué)報（自然科學(xué)版）,2011,25 (4)：340-343.

[13]寧立波,劉新宇,仝曉霞.基于物元分析的武漢市水資源承載力分析[J].安全與環(huán)境學(xué)報,2008,8(4)：76-80.

[14]李濱勇,史正濤,董銘.水資源承載力研究現(xiàn)狀與發(fā)展趨勢[J].節(jié)水灌溉,2007(2)：40-42.

[15]張曄.淺析廣州市水資源承載力[J].中國水運,2010,10 (6)：143-144.

Evaluation on water resources carrying capacity in Guangzhou City using principal components analysis.

Chen W,Huang C S and Wang N T.

Water resources carrying capacity is an important indicator to measure sustainable economic and social development.In this paper,by taking the water resources carrying capacity of Guangzhou City as the research object,based on the available data,using principal components analysis to extract two principal components that affect the dynamic changes of water resources carrying capacity,and making a comprehensive evaluation of annual variation of Guangzhou City.The results show that Guangzhou water resources carrying capacity from 2005 to 2010 was declining,but so far as the status,by extending the water supply project and strengthening water saving measures to relieve the imbalance between supply and demand,the water resources carrying capacity of Guangzhou still exists some rising space.

principal components analysis;water resources carrying capacity;Guangzhou City

P641.1

A

1007-3701(2013)04-322-05

10.3969/j.issn.1007-3701.2013.04.008

2013-07-10；

2013-11-15.

中國地質(zhì)調(diào)查局地質(zhì)調(diào)查項目“珠三角-環(huán)北部灣經(jīng)濟區(qū)地質(zhì)環(huán)境調(diào)查評價綜合研究”（1212011220012）資助。

陳雯（1985—），女，工程師，研究方向環(huán)境地質(zhì)、水文地質(zhì)，E-mail：diandian259@tom.com.