李曉菲

摘? ?要：基于統(tǒng)計(jì)模型PCA-AHP構(gòu)建評(píng)估水資源稀缺的綜合指數(shù)WSCEI（Water Scarcity Comprehensive Evaluation Indicator），并嵌入虛擬水流動(dòng)核算，從時(shí)間和空間視角綜合評(píng)價(jià)水資源的稀缺性。結(jié)果表明：一是商品中內(nèi)嵌的虛擬水理應(yīng)起到調(diào)節(jié)水資源稀缺分布不均衡的作用，然而與之相反，虛擬水是從缺水地區(qū)輸送至富水地區(qū)，這使中國(guó)水資源稀缺分布的不平衡問(wèn)題更加突出。二是我國(guó)的水資源稀缺壓力逐漸由“北高南低、西高東低”逐漸轉(zhuǎn)變?yōu)椤澳媳苯愿?、東西皆高”。三是經(jīng)濟(jì)較發(fā)達(dá)地區(qū)的水資源情況并未因其水資源利用相對(duì)效率和水資源經(jīng)濟(jì)效率的高效而有所改善；經(jīng)濟(jì)欠發(fā)達(dá)地區(qū)的缺水問(wèn)題則更為嚴(yán)重。

關(guān)鍵詞：水資源；稀缺性；虛擬水

中圖分類號(hào)：C82? ? ? ? 文獻(xiàn)標(biāo)志碼：A? ? ? 文章編號(hào)：1673-291X（2023）13-0049-03

一、研究背景

關(guān)于水資源評(píng)估許多學(xué)者提出了一些綜合性評(píng)價(jià)指標(biāo)。例如，WTA比率、水應(yīng)力指標(biāo)、供水壓力指數(shù)。另一方面，由于水資源的污染和浪費(fèi)，中國(guó)的水質(zhì)問(wèn)題也十分嚴(yán)峻。關(guān)于水質(zhì)問(wèn)題，雖然有一些研究，但僅僅評(píng)估了水質(zhì)狀況，并未對(duì)水資源短缺情況進(jìn)行評(píng)估。如李連香等（2015）基于分層構(gòu)權(quán)主成分分析，對(duì)皖北地下水水質(zhì)進(jìn)行了評(píng)估，但并未進(jìn)一步深入研究該地區(qū)的水資源短缺問(wèn)題。此外，水資源環(huán)境與社會(huì)經(jīng)濟(jì)系統(tǒng)也有著密不可分的關(guān)聯(lián)性。劉寧等（2020）基于虛擬水核算方法對(duì)京津冀地區(qū)農(nóng)產(chǎn)品虛擬水含量進(jìn)行計(jì)算，并通過(guò)建立回歸模型分析農(nóng)產(chǎn)品虛擬水與經(jīng)濟(jì)增長(zhǎng)的關(guān)系。

近年來(lái)，隨著水足跡理論的成熟和發(fā)展，許多學(xué)者將其用于水資源評(píng)估中，例如評(píng)價(jià)指標(biāo)IW、QQE等。但這些研究并未考慮到地區(qū)間的經(jīng)濟(jì)需求與交流，而且只是局限于一個(gè)區(qū)域的水資源評(píng)估。此外，大多數(shù)文章僅從時(shí)間角度討論了一個(gè)地區(qū)的水資源變化，而忽略了水資源可以通過(guò)產(chǎn)品運(yùn)輸從一個(gè)地區(qū)流向另一個(gè)地區(qū)。虛擬水理論的發(fā)展則從空間角度為水資源的評(píng)估提供了新思路。孫思奧、王晶等（2020）測(cè)算了青藏高原與中國(guó)其他區(qū)域之間的虛擬水貿(mào)易關(guān)系，建立了中國(guó)區(qū)域城鎮(zhèn)與農(nóng)村地區(qū)的虛擬水貿(mào)易網(wǎng)絡(luò)。魏怡然、邵玲等（2019）測(cè)度北京市與國(guó)內(nèi)其他省市和其他國(guó)家的虛擬水貿(mào)易往來(lái)。朱智洺、李秀琴等（2020）基于彭曼公式對(duì)“一帶一路”沿線國(guó)家糧食的單位面積虛擬水含量進(jìn)行測(cè)算。然而有關(guān)虛擬水的文章大多關(guān)注于它的基本理論、它在貿(mào)易中的重要性以及使用不同方法來(lái)計(jì)算其值，例如王勇（2016）基于水資源擴(kuò)展型MRIO模型和中國(guó)區(qū)域間投入產(chǎn)出表，對(duì)中國(guó)區(qū)域間貿(mào)易隱含虛擬水轉(zhuǎn)移進(jìn)行了測(cè)算，并發(fā)現(xiàn)虛擬水的流動(dòng)并未緩解中國(guó)的缺水問(wèn)題，但卻并未進(jìn)一步將虛擬水引入水資源評(píng)價(jià)體系中。事實(shí)上，僅有少數(shù)文章將虛擬水視為評(píng)估水資源的一個(gè)重要指標(biāo)。因此，可以基于統(tǒng)計(jì)模型PCA-AHP構(gòu)建評(píng)估水資源稀缺的綜合指數(shù)WSCEI（Water Scarcity Comprehensive Evaluation Indicator），并嵌入?yún)^(qū)域中的虛擬水流動(dòng)核算，從時(shí)間和空間視角綜合評(píng)價(jià)水資源稀缺性。

二、數(shù)據(jù)與方法

（一）數(shù)據(jù)來(lái)源

研究數(shù)據(jù)主要包括：2010年和2019年各區(qū)域的水質(zhì)、水量以及社會(huì)經(jīng)濟(jì)水平數(shù)據(jù)，來(lái)源于《中國(guó)統(tǒng)計(jì)年鑒》和《中國(guó)環(huán)境統(tǒng)計(jì)年鑒》，用于計(jì)算WSCEI。為了便于計(jì)算和分析，根據(jù)產(chǎn)業(yè)結(jié)構(gòu)、經(jīng)濟(jì)水平和地理位置的相似性，將我國(guó)劃分為8個(gè)區(qū)域：東北地區(qū)、京津地區(qū)、北部沿海地區(qū)、東部沿海地區(qū)、南部沿海地區(qū)、中部地區(qū)，西北地區(qū)和西南地區(qū)。

（二）研究方法

本文基于PCA-AHP統(tǒng)計(jì)模型構(gòu)建了一個(gè)三層次的指標(biāo)體系WSCEI，它綜合考慮了水質(zhì)（QL）、水量（QT）和社會(huì)經(jīng)濟(jì)水平（SC）的影響，從水質(zhì)、水量和社會(huì)經(jīng)濟(jì)水平三方面選擇了15個(gè)基本指標(biāo)。水質(zhì)水平選擇了3個(gè)指標(biāo)，分別為廢水排放比例、城市污水處理率和水土流失治理面積比例；水量水平選擇了9個(gè)指標(biāo)，分別為水資源總量、總用水量、水庫(kù)總庫(kù)容量、降水量、節(jié)水量、重復(fù)利用率、森林覆蓋率、濕地面積比例和虛擬水凈流出；社會(huì)經(jīng)濟(jì)水平選擇了3個(gè)指標(biāo)，分別為單位用水量產(chǎn)生的GDP、人均用水量和城市化比率。

為了分析中國(guó)水資源情況的動(dòng)態(tài)變化，選擇2010年和2019年各區(qū)域的水資源數(shù)據(jù)進(jìn)行橫向?qū)Ρ?。同時(shí)，通過(guò)比較是否引入NOVW的WSCEI值的變化，確認(rèn)虛擬水對(duì)水資源評(píng)估的影響以及區(qū)域間水資源分布的空間相關(guān)性。

三、結(jié)果分析與討論

（一）WSCEI等級(jí)劃分

在確保數(shù)據(jù)信息完整性的前提下，對(duì)于評(píng)估水量的基本指標(biāo)，使用PCA來(lái)降低數(shù)據(jù)維度。在未引入NOVW和引入NOVW的兩種情況下，特征值的累積方差貢獻(xiàn)率分別為85.406％和82.606%。根據(jù)累積方差貢獻(xiàn)率大于80％的原則，本文選擇了前三個(gè)主成分進(jìn)行之后的分析?；谒x的三個(gè)主成分，得到了兩種情況下的成份矩陣（因子載荷陣）以及相應(yīng)的得分系數(shù)。使用所獲得的得分系數(shù)分別構(gòu)建引入變量NOVW前后的三個(gè)主成分，從而實(shí)現(xiàn)水量數(shù)據(jù)降維的目的。

之后基于AHP計(jì)算各項(xiàng)指標(biāo)的權(quán)重，構(gòu)建出綜合評(píng)價(jià)指標(biāo)WSCEI。由于虛擬水僅對(duì)水量產(chǎn)生影響，因此QL，SC和WSCEI的方程式與未引入NOVW時(shí)相同。根據(jù)WSCEI值和我國(guó)水資源短缺的實(shí)際情況，將WSCEI劃分為四個(gè)層次，以顯示8個(gè)地區(qū)的缺水相對(duì)程度。為保證輕度缺水地區(qū)占比控制在50%，選取WECEI指標(biāo)的下四分之一分位數(shù)和上四分之一分位數(shù)，經(jīng)計(jì)算下四分之一分位數(shù)約為0.4 和上四分之一分位數(shù)為0.7，同時(shí)將WSCEI<0界定為嚴(yán)重缺水，WSCEI>0.7（上四分之一分位數(shù)）界定為不缺水，進(jìn)而制定WSCEI等級(jí)劃分標(biāo)準(zhǔn)。

總體來(lái)看，我國(guó)的缺水情況并未得到有效改善，甚至部分地區(qū)缺水問(wèn)題愈加嚴(yán)重。從空間角度看，由于虛擬水的引入，部分地區(qū)的水資源狀況確實(shí)發(fā)生了變化。然而，虛擬水的流動(dòng)并沒(méi)有改善我國(guó)水資源分布不平衡的問(wèn)題。相反，它使缺水地區(qū)的水資源短缺問(wèn)題更加嚴(yán)重，富水地區(qū)的水資源狀況卻得到明顯改善。

（二）WSCEI的時(shí)空特征及其因素分解

盡管京津地區(qū)的WSCEI值有所提高，但其仍處于嚴(yán)重缺水狀態(tài)，是八個(gè)地區(qū)中缺水最為嚴(yán)重的區(qū)域。西北地區(qū)似乎脫離了嚴(yán)重缺水狀態(tài)，其WSCEI值從-0.238增加到0.067。但0.067幾乎落在嚴(yán)重缺水和中度缺水的邊界，說(shuō)明這種緩解是非常脆弱的，存在回到嚴(yán)重缺水狀態(tài)的可能性。西南地區(qū)一直在嚴(yán)重缺水和中度缺水之間的邊界線上移動(dòng)，表明西南地區(qū)的缺水問(wèn)題十分嚴(yán)峻，存在著進(jìn)入嚴(yán)重缺水狀態(tài)的可能性。南部沿海地區(qū)和中部地區(qū)下降幅度最大，其中，南部沿海地區(qū)的WSCEI從0.518下降到0.208；中部地區(qū)的WSCEI則從1.160降至0.957，說(shuō)明這兩個(gè)地區(qū)水資源情況惡化較為嚴(yán)重。

根據(jù)WSCEI的準(zhǔn)則層指標(biāo)，可以深入分析各地區(qū)水資源短缺的具體原因，這是其他水資源評(píng)估指標(biāo)很難實(shí)現(xiàn)的。盡管整體水質(zhì)略有改善，但多數(shù)地區(qū)的水量明顯降低，這是造成我國(guó)缺水問(wèn)題愈發(fā)嚴(yán)重的主要原因。此外，東北地區(qū)和中部地區(qū)的SC值變化最為明顯。部分地區(qū)的QL值有所下降（北部沿海地區(qū)、東部沿海地區(qū)、西北地區(qū)及西南地區(qū)）。相對(duì)而言，其他地區(qū)的水質(zhì)有所改善，且這些地區(qū)QL值的變化更為明顯。例如，在這十年中，中部地區(qū)的水質(zhì)指標(biāo)增長(zhǎng)最快，QL值從0.72增加到0.86；而在北部和東部沿海地區(qū)，QL下降最快，但其減少量也僅為0.04。此外，京津地區(qū)的水質(zhì)最差，其QL值從0.35增加到0.44；沿海地區(qū)和中部地區(qū)水質(zhì)較好，其QL值穩(wěn)定在區(qū)間（0.6，0.8）中。由此可見(jiàn)，我國(guó)的水質(zhì)問(wèn)題雖然有所改善，但變化程度很弱，因此對(duì)各地區(qū)的缺水情況影響不明顯。

我國(guó)某些地區(qū)的水量有所增加，但總體而言，水量仍在下降。這主要體現(xiàn)在我國(guó)水量增長(zhǎng)地區(qū)的變化程度較弱，而水量下降地區(qū)的變化則十分明顯。因此，相較于2010年，我國(guó)2019年的水量大大減少。具體而言，西北地區(qū)的QT值增長(zhǎng)最快，增幅接近0.6；而南部沿海地區(qū)則下降最快，其QT值從0.47減至-0.16。此外，京津地區(qū)的水量最少，2010年和2019年的QT值分別為-1.76和-1.47。相比之下，中部地區(qū)QT值最大，但也從2010年的1.60降至2019年的1.12?？梢?jiàn)，QT值的降幅明顯高于其增幅，這也是我國(guó)水資源短缺情況一直沒(méi)有明顯改善的原因之一。

多數(shù)地區(qū)社會(huì)經(jīng)濟(jì)水平的變化并不明顯，京津地區(qū)的SC值最大（0.83），西北地區(qū)的SC值最?。?.15）；東北地區(qū)和中部地區(qū)SC值變化最為明顯：東北地區(qū)的SC值從0.36降至0.30，中部地區(qū)的SC值則從0.49降至0.45，雖然這兩個(gè)地區(qū)SC的變化很弱，但確實(shí)影響到它們的水資源狀況。

（三）NOVW-WSCEI分析及其影響因素

相較引入NOVW前后的WSCEI，發(fā)現(xiàn)西北地區(qū)、東部沿海地區(qū)和南部沿海地區(qū)的缺水情況變化程度很大，這是由地區(qū)間虛擬水流動(dòng)引起的。而虛擬水的流出，使得西北地區(qū)和東北地區(qū)的缺水問(wèn)題更為嚴(yán)重。其中，西北地區(qū)WSCEI值減少最多，從-0.238降至-0.587；東北地區(qū)的WSCEI則從0.490下降到0.270。相比之下，由于虛擬水的流入，東部和南部沿海地區(qū)的水資源狀況則得到明顯改善。其中，南部沿海地區(qū)的WSCEI增長(zhǎng)最快，其值從0.518增加到0.903；京津地區(qū)的WSCEI則從-0.712增至-0.453。以上事實(shí)表明，虛擬水流動(dòng)大大影響了我國(guó)的缺水分布。

總體來(lái)看，NOVW的引入使得水量指標(biāo)QT變化很大，大多數(shù)地區(qū)QT值的變化與虛擬水的流向一致。其中，由于虛擬水流入，京津地區(qū)的QT值從-1.76增加到-1.26；東部沿海地區(qū)從0.61增加到0.89；南部沿海地區(qū)則從0.47增加到1.19。同時(shí)，由于虛擬水流出，東北地區(qū)的QT值從0.48降至0.06，而西北地區(qū)的QT值從-0.89降至-1.55。

根據(jù)上述分析，發(fā)現(xiàn)我國(guó)大部分地區(qū)的缺水問(wèn)題并未得到有效解決，部分地區(qū)的情況甚至惡化，這與多數(shù)的研究結(jié)果基本一致。表明WSCEI可以用于評(píng)估一個(gè)地區(qū)的缺水情況，其結(jié)果是準(zhǔn)確可靠的?；赪SCEI的變化和本文構(gòu)建的指標(biāo)體系，可以分析影響我國(guó)水資源短缺問(wèn)題日益嚴(yán)重的主要原因。

從指標(biāo)體系的準(zhǔn)則層來(lái)看，我國(guó)水質(zhì)總體上有所改善，但因其變化太小，無(wú)法有效提高WSCEI的值。即使某些地區(qū)的水質(zhì)明顯改善，由于其水量銳減，該地區(qū)的缺水問(wèn)題也并未得到緩解，而是越來(lái)越嚴(yán)重。社會(huì)經(jīng)濟(jì)水平的變化也并不明顯。然而，它確實(shí)影響到地區(qū)的水資源狀況（例如東北地區(qū)）。另外，根據(jù)NOVW引入前后的WSCEI變化，發(fā)現(xiàn)虛擬水流動(dòng)對(duì)水量分布的影響很大，這可能是緩解我國(guó)水資源短缺問(wèn)題的一個(gè)有效途徑。對(duì)于東部沿海和西北地區(qū)水資源情況的改善，我們進(jìn)行深入探討并分析了原因。作為中國(guó)最發(fā)達(dá)的地區(qū)，東部沿海地區(qū)擁有最先進(jìn)的水資源再利用技術(shù)，再加上虛擬水的大量流入，使得其缺水狀況得到有效改善。相比之下，我國(guó)政府雖然發(fā)布了一系列政策來(lái)幫助西北地區(qū)，特別是南水北調(diào)工程，使得西北地區(qū)的缺水問(wèn)題有所改善，但由于地理位置的限制和虛擬水的流出，西北地區(qū)的缺水問(wèn)題依然相當(dāng)嚴(yán)峻。

水量的大幅下降導(dǎo)致中國(guó)水資源問(wèn)題日益嚴(yán)重。近年來(lái)，由于水資源的過(guò)度消耗和全球氣候變暖的影響，我國(guó)水資源量明顯減少。雖然我國(guó)的節(jié)水技術(shù)和水資源再利用技術(shù)有所進(jìn)步，但隨著生活水平的提高，人們對(duì)水資源的需求越來(lái)越大，導(dǎo)致用水量明顯增加。此外，由于虛擬水的引入，各地區(qū)的水量也發(fā)生了很大變化，這使得每個(gè)地區(qū)的水資源短缺程度有了明顯變化。它表明虛擬水流動(dòng)確實(shí)會(huì)影響到水資源分布。然而，虛擬水的流動(dòng)并未改善我國(guó)水資源分布的不均衡，反而使缺水地區(qū)的缺水問(wèn)題越發(fā)嚴(yán)重，“富水”地區(qū)的缺水問(wèn)題卻得到改善。這是因?yàn)樘摂M水從缺水地區(qū)（例如西北地區(qū)）流向“富水”地區(qū)（例如東部沿海地區(qū)），使得我國(guó)水資源分布更加不均勻，但有虛擬水流入的地區(qū)缺水情況的確明顯改善。這些事實(shí)表明，水資源的狀況與虛擬水的流動(dòng)密切相關(guān)。

參考文獻(xiàn)：

[1]? ?李連香，許迪，程先軍，等.基于分層構(gòu)權(quán)主成分分析的皖北地下水水質(zhì)評(píng)價(jià)研究[J].資源科學(xué)，2015，37（1）：61-67.

[2]? ?劉寧，沙景華，張宏亮，等.京津冀地區(qū)農(nóng)產(chǎn)品虛擬水與區(qū)域可持續(xù)發(fā)展關(guān)系研究[J].資源與產(chǎn)業(yè)，2016，18（3）：80-85.

[3]? ?孫思奧，王晶，戚偉.青藏高原地區(qū)城鄉(xiāng)虛擬水貿(mào)易格局與影響因素[J].地理學(xué)報(bào)，2020，75（7）：46-58.

[4]? ?魏怡然，邵玲，張寶剛，等.北京市虛擬水消費(fèi)與貿(mào)易[J].自然資源學(xué)報(bào)，2019，34（9）：62-73.

[5]? ?朱智洺，李秀琴，黃永春.“一帶一路”沿線國(guó)家糧食虛擬水貿(mào)易網(wǎng)絡(luò)研究[J].資源與產(chǎn)業(yè)，2020，22（3）：31-42.

[6]? ?王勇.全行業(yè)口徑下中國(guó)區(qū)域間貿(mào)易隱含虛擬水的轉(zhuǎn)移測(cè)算[J].中國(guó)人口·資源與環(huán)境，2016，26（4）：107-114.

Research on Evaluation of Water Resources Scarcity and Accounting of Embedded Virtual Water Based on Industry Wide Statistics

Li Xiaofei

（Dalian Neusoft University of Information， Dalian 116023， China）

Abstract： Based on the statistical model PCA-AHP， the comprehensive index WSCEI for evaluating water scarcity is constructed， and virtual water flow accounting is embedded to comprehensively evaluate water scarcity from the perspective of time and space.The research reveals the following points： first， virtual water embedded in the commodity ought to play an important role in regulating the imbalance of water scarcity distribution， but in contrast， Chinas virtual water flows from water-scarce areas to water-abundant areas， which makes the imbalance of water scarcity distribution in China more severe; second， the scarcity pressure of water resources in China has gradually shifted from “high in the north and low in the south， high in the west and low in the east” to “high in both the north and south， and high in both the east and west”; third， the situation of water resources in economically developed regions has not been improved due to the relative efficiency of water resources utilization and the high efficiency of economic efficiency of water resources; the problem of water shortage in economic underdevelopment region is even more serious.

Key words： water resource; water scarcity; virtual water

[責(zé)任編輯? ?白? ?雪]