嚴(yán) 冬，桂東偉，薛 杰，孫懷衛(wèi)，廖衛(wèi)紅

(1.華中科技大學(xué)水電與數(shù)字化工程學(xué)院，武漢 430074；2.中國科學(xué)院新疆生態(tài)與地理研究所 新疆策勒荒漠草地生態(tài)系統(tǒng)國家野外科學(xué)觀測研究站，烏魯木齊 830011；3.中國水利水電科學(xué)研究院 水資源研究所，北京 100038)

由于水資源稟賦以及水資源利用效率不同，國家或地區(qū)間生產(chǎn)相同產(chǎn)品的水資源成本具有較大差異。中國水資源分布和社會需求間存在巨大的矛盾，通過地區(qū)間虛擬水交易，利用各類商品貿(mào)易中蘊含的虛擬水量，可在一定程度上緩解這一矛盾。但由于自然水資源分布和社會生產(chǎn)耗水狀況之間的格局正在被氣候變化所改變，虛擬水流動這樣一種自然水資源間接配置方式所發(fā)揮的作用也在發(fā)生變化。評價這種變化，將有助于更全面地了解未來可能出現(xiàn)的水資源供需狀況。因此，針對我國的地區(qū)間虛擬水流動特征進行分析，并對氣候變化情景下虛擬水流動的響應(yīng)進行評價，具有十分重要的理論與現(xiàn)實意義。

當(dāng)前，已有眾多關(guān)于虛擬水含量[1-4]、虛擬水流動[5-9]、氣候變化對用水[10-13]及水資源影響[14-17]的研究，它們?yōu)橛嬎闾摂M水流向、流量、評估氣候變化下水資源供需狀況提供了堅實的理論和方法。但尚未見到分析氣候變化條件下虛擬水流動響應(yīng)特征的研究。為此，本文以中國地區(qū)間農(nóng)業(yè)灌溉虛擬水流量計算為基礎(chǔ)，運用虛擬水流動強度、輸入份額和不均衡程度等指標(biāo)來評價虛擬水流動狀況，并通過分析各地區(qū)水資源和農(nóng)業(yè)灌溉凈需水量在氣候變化情景下的波動，來評價虛擬水流動的響應(yīng)特征。

1 評價方法

1.1 地區(qū)間虛擬水流動

結(jié)合可獲取的資料，定義和計算地區(qū)間虛擬水流動特征。地區(qū)間投入產(chǎn)出表給出了不同地區(qū)、不同部門(如農(nóng)業(yè)、工業(yè)、建筑業(yè)等)間商品的流動情況，結(jié)合各商品的虛擬水含量，可計算各地區(qū)間的虛擬水流動量。區(qū)域間投入產(chǎn)出平衡關(guān)系可表示為：

(1)

以該式為基礎(chǔ)，計算各部門地區(qū)間虛擬水流量之和為：

(2)

(3)

(4)

1.2 虛擬水流動評價指標(biāo)

為合理評估虛擬水流動特征，提出采用虛擬水流動強度、輸入份額和不均衡程度等3個評價指標(biāo)對其進行評價，定義如下。

(5)

(3) 虛擬水流動不均衡程度。借鑒基尼系數(shù)的概念和計算公式[18]，整體評價地區(qū)間虛擬水流動格局相對水資源和經(jīng)濟總量分布的不均衡程度。首先，分別計算各地區(qū)水資源量占全國水資源總量比例、各地區(qū)國內(nèi)生產(chǎn)總值占全國國內(nèi)生產(chǎn)總值比例、各地區(qū)虛擬水輸出量占全國虛擬水輸出總量比例以及各地區(qū)虛擬水輸入量占全國虛擬水輸入總量比例。然后，分別按各地區(qū)水資源量和國內(nèi)生產(chǎn)總值的大小，從低到高對各地區(qū)進行排序。根據(jù)排序結(jié)果，從水資源數(shù)量最低地區(qū)開始，依次計算各地區(qū)水資源量和虛擬水輸出量的累積比例；從國內(nèi)生產(chǎn)總值最低地區(qū)開始，依次計算各地區(qū)國內(nèi)生產(chǎn)總值和虛擬水輸入量的累積比例。據(jù)此，虛擬水輸出格局相對水資源分布、虛擬水輸入相對經(jīng)濟總量分布的不均衡程度GE、GF分別為：

(6)

(7)

式中：ELk、FLk分別為k地區(qū)的虛擬水輸出、輸入量；RLk、GLk分別為k地區(qū)水資源量和國內(nèi)生產(chǎn)總值的累積比例。

與基尼系數(shù)一樣，不均衡程度也是一個介于0和1之間的小數(shù)。

1.3 氣候變化的影響

虛擬水流量和評價指標(biāo)與地區(qū)耗水和水資源相關(guān)，而這兩者與氣候變化聯(lián)系緊密。因缺乏工業(yè)、建筑業(yè)、第三產(chǎn)業(yè)等部門耗水受氣候變化影響的資料，本文僅針對氣候變化通過農(nóng)業(yè)耗水變化而對地區(qū)間虛擬水流動量及相關(guān)評價指標(biāo)產(chǎn)生的影響。需要指出的是，氣候變化除了影響農(nóng)產(chǎn)品耗水之外，還會影響農(nóng)產(chǎn)品產(chǎn)量[19,20]和需求[21]，改變農(nóng)產(chǎn)品供需格局，進而影響農(nóng)產(chǎn)品流通狀況。因此，氣候變化對虛擬水流動的影響包括了農(nóng)產(chǎn)品耗水和流通兩方面。由于農(nóng)產(chǎn)品流通受到農(nóng)產(chǎn)品價格、運輸費用、糧食購銷政策等多方面因素的影響，其變動機制十分復(fù)雜，因此本文僅從農(nóng)產(chǎn)品耗水方面來進行分析。

本文中農(nóng)業(yè)部門年耗水量僅指各類作物年灌溉需水量之和(相應(yīng)的，后文中的虛擬水及相關(guān)計算結(jié)果都是關(guān)于農(nóng)業(yè)灌溉虛擬水的)。對某一類作物，其年灌溉需水量是作物面積和單位面積年灌溉需水量的乘積，而單位面積年灌溉需水量是各月單位面積需水量的和。對某一個月，其單位面積灌溉需水量C按下式近似估算：

C=ET-P

(8)

ET=kET0

(9)

式中：ET為該月作物騰蒸發(fā)量；E為該月有效降水；k為作物系數(shù)[22]；ET0為月參考作物騰發(fā)量。

本文采用彭曼公式的替代公式[23]計算求得ET0，采用USDA-SCS方法來計算求得P。年徑流深R計算公式[24]為：

R=Pyexp(-Ey/Py)

(10)

式中：Py為年降水量；Ey為當(dāng)年各月潛在蒸發(fā)量的和。

(11)

(12)

2 數(shù) 據(jù)

采用了中國科學(xué)院地理科學(xué)與資源研究所區(qū)域可持續(xù)發(fā)展分析與模擬重點實驗室編制的中國2007年30省(區(qū)、市)(除了西藏、香港、澳門和臺灣)區(qū)域間6部門(農(nóng)林牧漁業(yè)、工業(yè)、建筑業(yè)、交通運輸及倉儲業(yè)、批發(fā)零售業(yè)、其他服務(wù)業(yè))投入產(chǎn)出表[25]。該表單列了各地區(qū)對進口品和出口品的使用情況，其中間投入反映了國內(nèi)生產(chǎn)商品在各地區(qū)各部門的使用情況。各省份水資源量和工農(nóng)業(yè)耗水?dāng)?shù)據(jù)來自2007年中國水資源公報；國內(nèi)生產(chǎn)總值、農(nóng)作物類型和相應(yīng)種植面積數(shù)據(jù)來自2008年中國統(tǒng)計年鑒。

溫室氣體排放情景，是對未來氣候變化預(yù)估的基礎(chǔ)。最新的情景設(shè)計中，采用了“典型濃度目標(biāo)”(Representative Concentration Pathways, RCPs)，其中：RCP8.5情景為最高的溫室氣體排放情景，主要指輻射強迫上升至8.5 W/m2；RCP4.5情景是2100年輻射強迫穩(wěn)定在4.5 W/m2，考慮了與全球經(jīng)濟框架相適應(yīng)的、長期存在的全球溫室氣體和生存期短的物質(zhì)排放；RCP2.6是把全球平均溫度上升限制在2 ℃之內(nèi)的情景。采用北京氣候中心氣候系統(tǒng)模式(BCC-CSM1.1)[26]在CIMIP5試驗情景RCP2.6、RCP4.5和RCP8.5[27]下模擬的2006-2099年降水量和氣溫數(shù)據(jù)；基準(zhǔn)期設(shè)為2006-2010年，并假設(shè)基準(zhǔn)期各省份水資源量和農(nóng)業(yè)單位產(chǎn)值耗水量等于2007年數(shù)值。

3 結(jié)果分析

3.1 農(nóng)業(yè)萬元產(chǎn)值耗水量

農(nóng)業(yè)萬元產(chǎn)值耗水量基準(zhǔn)期數(shù)值、非基準(zhǔn)期相對基準(zhǔn)期的變幅見表1?？梢?，各省市農(nóng)業(yè)萬元產(chǎn)值耗水量彼此差異巨大，在氣候變化條件下的波動也很大，變幅最大值與最小值的差值，最大達(dá)225%(四川)，最小為28%(內(nèi)蒙古)，平均為72%。此外，變幅最大值都大于最小值絕對值，并且僅有兩個省市(吉林、山東)的變幅中位數(shù)小于零，這表明在未來多數(shù)時期農(nóng)業(yè)灌溉耗水量將增加，這一結(jié)果也與兩省都為農(nóng)業(yè)大省的現(xiàn)狀比較相符。

表1 農(nóng)業(yè)萬元產(chǎn)值灌溉耗水量Tab.1 Agricultural water consumption of ten thousand yuan output value

3.2 農(nóng)業(yè)灌溉虛擬水流動量

如圖1所示，基準(zhǔn)期虛擬水輸出量最大的10個省市(自治區(qū))由高到低依次是新疆、黑龍江、河北、內(nèi)蒙古、安徽、廣西、湖南、江蘇、江西和吉林，占全部省市虛擬水輸出總量的75.3%；基準(zhǔn)期虛擬水輸入量最大的10個省市(由高到低)是山東、上海、廣東、江蘇、河北、天津、河南、浙江、北京和吉林，占全部省市虛擬水輸入總量的71.6%。這歸因于兩方面因素。一是這10個省市多為發(fā)達(dá)地區(qū)，對外省市產(chǎn)品需求量大；二是外省市生產(chǎn)耗水量大，產(chǎn)品虛擬水含量高。

圖1 基準(zhǔn)期各省市虛擬水輸入、輸出量Fig.1 Virtual water imports and exports of all provinces in base period

3.3 農(nóng)業(yè)灌溉虛擬水流動強度

匯總RCP26、RCP45和RCP85情景下未來各時期計算成果，統(tǒng)計得到各省份虛擬水流動強度的最大值、中位數(shù)、最小值和基準(zhǔn)值，見圖2和圖3。虛擬水輸出強度較高的省市是天津、河北、寧夏和新疆，其數(shù)值遠(yuǎn)遠(yuǎn)高于其他省份，組成了虛擬水輸出第一梯隊；虛擬水凈輸入強度較高的省市是天津、北京和上海。其中，天津的虛擬水輸出和輸入強度都是全國最大，但因為虛擬水輸入強度太高，以致虛擬水凈輸入強度的各個統(tǒng)計值遠(yuǎn)高于其他省市。

在虛擬水凈輸入強度為負(fù)值的省份中，河北、寧夏和新疆由于虛擬水輸出強度大，其凈輸出強度顯著大于其他省市。結(jié)果表明，擁有高虛擬水流動強度的省市具有更大的變幅，同時，高強度省份存在明顯差異，表明這些省份的虛擬水輸出強度在未來多數(shù)時期內(nèi)存在明顯的增加或下降等變化。這些情況說明，高虛擬水流動強度省市受氣候變化影響的程度更大。

圖2 虛擬水輸出強度Fig.2. The intensity of virtual water export

圖3 虛擬水凈輸入強度Fig.3 The intensity of virtual water net import

圖4 虛擬水輸入份額Fig.4. The scale of virtual water import

3.4 農(nóng)業(yè)灌溉虛擬水輸入份額

各省市虛擬水輸入份額統(tǒng)計值見圖4。由圖4可見，具有極高輸入份額的北京、天津和上海這3個直轄市組成了第一梯隊，河北、吉林和山東組成了第二梯隊，浙江、廣東、河南、重慶和陜西則組成了第三梯隊；這些梯隊輸入份額都在40%以上。第一梯隊的中位數(shù)都高于基準(zhǔn)值，說明未來大部分時期虛擬水輸入對這些地區(qū)的重要性還將進一步增加。同時，高虛擬水輸入強度屬于經(jīng)濟發(fā)達(dá)、人均水資源量較少的地區(qū)，高虛擬水輸出強度省市多屬于經(jīng)濟中等發(fā)達(dá)地區(qū)，輸入份額較大的省市多屬于經(jīng)濟發(fā)達(dá)或中等發(fā)達(dá)地區(qū)。在氣候變化條件下這些省市受到的影響比其他地區(qū)要大得多。

3.5 農(nóng)業(yè)灌溉虛擬水流動不均衡程度

綜合各個情景下的計算結(jié)果，在基準(zhǔn)期年和氣候變化條件下，虛擬水流動不均衡程度見表2?？梢?，虛擬水輸出格局相對水資源分布的不均衡程度數(shù)值較高，但虛擬水輸入格局相對經(jīng)濟總量分布的不均衡程度數(shù)值較低。

表2還列出了在基準(zhǔn)期年和氣候變化條件下水資源分布相對經(jīng)濟總量分布的不均衡程度。對比數(shù)值，發(fā)現(xiàn)農(nóng)業(yè)灌溉虛擬水輸出格局不均衡程度(基準(zhǔn)值為0.511)高于水資源分布(基準(zhǔn)值為0.402)，表明農(nóng)業(yè)灌溉虛擬水輸出加強了水資源分布的不均衡性。但虛擬水輸入格局不均衡程度(基準(zhǔn)值為0.142)低于水資源分布，說明農(nóng)業(yè)灌溉虛擬水輸入緩解了水資源分布不均的影響。以上虛擬水流動對水資源分布均衡性產(chǎn)生的正面和負(fù)面影響，說明我國經(jīng)濟總量占全國比例較高的地區(qū)輸入的虛擬水量多來自水資源量占全國比例較低的地區(qū)；比較氣候變化條件下不均衡程度的變化，氣候變化對該情勢沒有顯著影響。

表2 虛擬水流動不均衡程度Tab.2 The extent of virtual water inflow

4 結(jié) 語

本文建立了農(nóng)業(yè)灌溉虛擬水流動強度、輸入份額和不均衡程度等評價指標(biāo)，基于中國地區(qū)間投入產(chǎn)出表給出的省市間貿(mào)易關(guān)系，結(jié)合氣候模式輸出結(jié)果計算得到的水資源量和農(nóng)業(yè)單位產(chǎn)值耗水量變化情況，分析了中國30個省市間農(nóng)業(yè)灌溉虛擬水的流動特征及其對氣候變化的響應(yīng)。結(jié)果發(fā)現(xiàn)：

(1)農(nóng)業(yè)灌溉虛擬水流動強度、輸入份額和不均衡程度等評價指標(biāo)較好評價了我國地區(qū)間農(nóng)業(yè)灌溉虛擬水流動特征及其對氣候變化的影響。本文所建立的假設(shè)及研究方法，能夠結(jié)合資料的可獲取性，所獲得的結(jié)果與實際情況相符。

(2)從農(nóng)業(yè)灌溉虛擬水流動強度和輸入份額的結(jié)果可見，高虛擬水輸入強度屬于經(jīng)濟發(fā)達(dá)、人均水資源量較少的地區(qū)，高虛擬水輸出強度省市多屬于經(jīng)濟中等發(fā)達(dá)地區(qū)，輸入份額較大的省市多屬于經(jīng)濟發(fā)達(dá)或中等發(fā)達(dá)地區(qū)。

(3)對30個省市整體農(nóng)業(yè)灌溉虛擬水流動不均衡程度的分析表明，農(nóng)業(yè)灌溉虛擬水輸出格局的不均衡程度較水資源分布為高，而農(nóng)業(yè)灌溉虛擬水輸入格局的不均衡程度比水資源分布要低。氣候變化對農(nóng)業(yè)灌溉虛擬水輸出格局有顯著影響，使得未來大部分時期農(nóng)業(yè)灌溉虛擬水輸出格局變得更加不均衡。

需要指出的是，本文在不同氣候變化情景中都使用了同樣的地區(qū)間投入產(chǎn)出關(guān)系，這與實際情況有偏差。地區(qū)間農(nóng)產(chǎn)品供需狀況是影響地區(qū)間投入產(chǎn)出關(guān)系的一個主要因素，而農(nóng)產(chǎn)品的生產(chǎn)和消費都受到氣候變化的顯著影響。此外，市場和政策的波動也會影響地區(qū)間的產(chǎn)品流動。若能量化這些因素的變化對地區(qū)間投入產(chǎn)出關(guān)系的影響，則有望獲得對農(nóng)業(yè)灌溉虛擬水流動狀況更精確的評價結(jié)果。

□

參考文獻(xiàn)：

[1] 龍愛華,徐中民,張志強.虛擬水理論方法與西北4省(區(qū))虛擬水實證研究[J].地球科學(xué)進展,2004,(4):577-584.

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

[3] 郭衛(wèi)寧,孟成才,韓宇平.唐山市 2007 年農(nóng)產(chǎn)品虛擬水測評與分析[J].中國水利水電科學(xué)研究院學(xué)報,2010,(3):177-181.

[4] 黃茄莉,徐中民,王 康.甘州區(qū)社會經(jīng)濟系統(tǒng)水循環(huán)研究[J].水利學(xué)報,2010,(9):1 114-1 120.

[5] Manfred Lenzen.Understanding virtual water flows: a multiregion input-output case study of victoria[J]. Water Resources Research, 2009,45 (9):416-426.

[6] Dabo Guan, Klaus Hubacek.Assessment of regional trade and virtual water flows in China[J].Ecological Economics,2007,61(1):159-170.

[7] Hongrui Wang,Yan Wang.An input-output analysis of virtual water uses of the three economic sectors in Beijing[J].Water International,2009,34 (4):451-467.

[8] Y Wang,H L Xiao,M F Lua.Analysis of water consumption using a regional input-output model: model development and application to Zhangye City,Northwestern China[J].Journal of Arid Environments,2009,73(10):894-900.

[9] Kuishuang Feng,Yim Ling Siu,Dabo Guan,et al.Assessing regional virtual water flows and water footprints in the Yellow River Basin,China:a consumption based approach[J].Applied Geography,2012,32(2):691-701.

[10] 劉曉英,林而達(dá).氣候變化對華北地區(qū)主要作物需水量的影響[J].水利學(xué)報,2004,(2):77-82.

[11] 王小軍,賀瑞敏,尚熳廷.氣候變化對區(qū)域農(nóng)業(yè)灌溉用水影響分析[J].中國農(nóng)村水利水電,2011,(1):29-32.

[12] 吳普特,趙西寧.氣候變化對中國農(nóng)業(yè)用水和糧食生產(chǎn)的影響[J].農(nóng)業(yè)工程學(xué)報,2010,(2):1-6.

[13] 張山清,普宗朝,宋良孌,等.吐魯番地區(qū)氣候變化對參考作物蒸散量的影響[J].中國農(nóng)業(yè)氣象,2009,(4):532-537.

[14] 李峰平,章光新,董李勤.氣候變化對水循環(huán)與水資源的影響研究綜述[J].地理科學(xué),2013,(4):457-464.

[15] 劉春蓁.氣候變化對我國水文水資源的可能影響[J].水科學(xué)進展,1997,(3):220-225.

[16] 陳亞寧,徐宗學(xué).全球氣候變化對新疆塔里木河流域水資源的可能性影響[J].中國科學(xué)(D輯),2004,(11):1 047-1 053.

[17] 丁相毅,賈仰文,王 浩,等.氣候變化對海河流域水資源的影響及其對策[J].自然資源學(xué)報,2010,(4):604-613.

[18] Paul A Samuelson,William D Nordhaus.Economics,19th Edition[M]. McGraw Hill Higher Education,2009.

[19] 李小燕.黃土高原植被對水熱狀況的響應(yīng)研究[J].地理科學(xué), 2013,(7):865-872.

[20] 史本林,朱新玉,胡云川,等.夏玉米氣候適宜度對全球氣候變暖的響應(yīng)----以商丘地區(qū)為例[J].地理科學(xué),2013,(10):1 277-1 284.

[21] 方修琦,鄭景云,葛全勝.糧食安全視角下中國歷史氣候變化影響與響應(yīng)的過程與機理[J].地理科學(xué), 2014,(11):1 291-1 299.

[22] Allen R G,Pereira L S,Raes D,et al.Crop evapotranspiration-Guidelines for computing crop water requirements-FAO irrigation and drainage paper 56[EB/OL]. http:∥www.fao.org/docrep/X0490E/X0490E00.htm.

[23] Willmott C J,Rowe C M,Mintz Y.Climatology of the terrestrial seasonal water cycle[J].Journal of Climatology,1985,5(6):589-606.

[24] Gardner L R.Assessing the effect of climate change on mean annual runoff[J].Journal of Hydrology,2009,379(3-4):351-359.

[25] 劉衛(wèi)東,陳 杰,唐志鵬,等.中國2007年30省市區(qū)域間投入產(chǎn)出表編制理論與實踐[M].北京:中國統(tǒng)計出版社,2012.

[26] 吳統(tǒng)文,宋連春,李偉平,等.北京氣候中心氣候系統(tǒng)模式研發(fā)進展----在氣候變化研究中的應(yīng)用[J].氣象學(xué)報,2014,(1):12-29.

[27] 王紹武,羅 勇,趙宗慈,等.新一代溫室氣體排放情景[J].氣候變化研究進展,2012,(7):305-307.