黑龍江省糧食水足跡時(shí)空分布規(guī)律與影響因素分析

劉 巍 李雨晴 付 強(qiáng) 孟 軍 李天霄 戴 爽

(1.東北農(nóng)業(yè)大學(xué)水利與土木工程學(xué)院， 哈爾濱 150030； 2.東北農(nóng)業(yè)大學(xué)經(jīng)濟(jì)管理學(xué)院， 哈爾濱 150030)

0 引言

水足跡概念最先于2002年由HOEKSTRA提出，將其定義為：在一定時(shí)期內(nèi)(一般是1年)生產(chǎn)某一產(chǎn)品或提供某種服務(wù)所消耗的淡水資源量[1-2]。糧食作為一種耗水密集型產(chǎn)品，其水足跡能夠反映某一區(qū)域糧食生產(chǎn)過(guò)程中水資源的利用情況[3-4]。黑龍江省是我國(guó)重要的商品糧基地[5]，自2010年起，全省糧食產(chǎn)量連續(xù)9年位列全國(guó)第一[6]，然而，黑龍江省單位面積灌溉水量比全國(guó)平均值高5～6倍[7-8]，且水資源時(shí)空分布不均，農(nóng)業(yè)水資源利用效率低[9-10]。因此，有必要研究糧食藍(lán)水足跡、綠水足跡及總水足跡的時(shí)空變化規(guī)律及影響因素，使黑龍江省在穩(wěn)產(chǎn)的同時(shí)提高農(nóng)業(yè)水資源利用率。

隨著水足跡理論的發(fā)展，關(guān)于糧食水足跡的研究日漸增多。孫世坤等[11]利用水足跡理論分析了中國(guó)大陸小麥水足跡空間分布特征，并對(duì)造成水足跡區(qū)域差異的因素進(jìn)行歸因分析。CHAPAGAIN等[12]對(duì)世界主要生產(chǎn)國(guó)的水稻藍(lán)水足跡、綠水足跡、灰水足跡構(gòu)成進(jìn)行了量化和分析。徐猛等[13]利用CROPWAT模型評(píng)估了哈爾濱地區(qū)大豆的水足跡，為黑龍江省大豆主產(chǎn)區(qū)水肥統(tǒng)籌管理提供了依據(jù)。劉聰[14]計(jì)算了中國(guó)31個(gè)省的小麥、玉米和水稻3種主要糧食作物的水足跡，并對(duì)其空間分布特征進(jìn)行評(píng)價(jià)。張凡凡等[15]對(duì)全國(guó)31個(gè)省糧食水足跡的影響因素進(jìn)行了分析。通過(guò)對(duì)比可以發(fā)現(xiàn)，不同地區(qū)的糧食水足跡差異顯著。

針對(duì)黑龍江省區(qū)域特點(diǎn)及存在的農(nóng)業(yè)用水問(wèn)題，根據(jù)水足跡理論，本文分別對(duì)2008—2018年黑龍江省11個(gè)市級(jí)地區(qū)的主要糧食作物(水稻、大豆、小麥)的平均藍(lán)水足跡、綠水足跡及總水足跡[16]進(jìn)行評(píng)價(jià)，明晰糧食藍(lán)水足跡、綠水足跡及總水足跡的時(shí)空分布規(guī)律，分析其影響因素，并對(duì)未來(lái)5年糧食藍(lán)水足跡、綠水足跡及總水足跡進(jìn)行預(yù)測(cè)，以期為黑龍江省農(nóng)業(yè)水資源利用效率的提高及保障糧食安全提供科學(xué)參考。

1 材料與方法

1.1 研究區(qū)概況

黑龍江省位于中國(guó)東北部，全省年平均氣溫-5～5℃，年降水量400～650 mm[17-18]，多年灌溉水利用系數(shù)平均值為0.53，低于全國(guó)平均水平[19]。糧食生產(chǎn)用水是黑龍江省主要的農(nóng)業(yè)水資源消耗，占農(nóng)業(yè)用水總量的89.63%[20]。黑龍江省水資源總量并不豐富，且時(shí)空差異顯著[21-22]。黑龍江省共13個(gè)下轄地區(qū)，包括12個(gè)地級(jí)市和1個(gè)地區(qū)[23-24]，本文選取2008—2018年黑龍江省11個(gè)市區(qū)(大興安嶺、黑河地區(qū)由于氣候等因素不適合作物生長(zhǎng)，暫不予以研究)的主要糧食作物(水稻、大豆和小麥)為研究對(duì)象。計(jì)算糧食水足跡所需的地區(qū)年降水量，作物播種面積以及產(chǎn)量數(shù)據(jù)來(lái)自《黑龍江省統(tǒng)計(jì)年鑒》、黑龍江省水資源公報(bào)，作物灌溉定額來(lái)自黑龍江省灌溉信息網(wǎng)(http:∥www.hljggxx.gov.cn)。研究區(qū)域概況見(jiàn)圖1。

圖1 黑龍江省概況Fig.1 Location of Heilongjiang Province

1.2 研究方法

1.2.1斯皮爾曼相關(guān)系數(shù)法

斯皮爾曼相關(guān)系數(shù)法是根據(jù)2個(gè)集合中的元素在各自集合的等級(jí)來(lái)計(jì)算二者之間的相關(guān)性[25-26]。斯皮爾曼相關(guān)系數(shù)法對(duì)數(shù)據(jù)的要求較低，只需2個(gè)變量的觀測(cè)值是成對(duì)的等級(jí)，或是通過(guò)連續(xù)變量觀測(cè)資料轉(zhuǎn)換得到的等級(jí)數(shù)據(jù)，不論2個(gè)變量的總體分布形態(tài)、樣本容量，都可以用斯皮爾曼相關(guān)系數(shù)法進(jìn)行相關(guān)性研究。計(jì)算公式為

(1)

式中rs——斯皮爾曼系數(shù)

di——不同集合對(duì)應(yīng)變量之間的等級(jí)差

n——樣本容量

1.2.2時(shí)間序列模型

時(shí)間序列分析法[27-28]于1927年由英國(guó)統(tǒng)計(jì)學(xué)家YULE提出。時(shí)間序列分析法是通過(guò)曲線擬合和參數(shù)估計(jì)對(duì)觀測(cè)得到的時(shí)間序列數(shù)據(jù)建立數(shù)學(xué)模型的理論和方法[29]。時(shí)間序列模型包含Arima模型、Simple模型等7種模型。因本研究中不涉及趨勢(shì)和季節(jié)成分，故采用時(shí)間序列Simple模型，具體公式為

t+1=αxt+α(1-α)xt-1+α(1-α)2xt-2+…+α(1-α)t-1x1+(1-α)tl0

(2)

xt——t時(shí)刻的觀測(cè)值

α——平滑參數(shù)l0——初始值

1.2.3糧食水足跡

糧食水足跡在水足跡[30-31]概念上得來(lái)，其定義為某個(gè)區(qū)域生產(chǎn)糧食作物單位產(chǎn)量所消耗的水資源量。糧食水足跡主要包括糧食綠水足跡、糧食藍(lán)水足跡和糧食總水足跡[32-35]，為三者的統(tǒng)稱。糧食藍(lán)水足跡是指農(nóng)田中可被作物利用的灌溉地表水和地下水的消耗量，反映了糧食生產(chǎn)過(guò)程中消耗的灌溉水資源；糧食綠水足跡是指雨水落入農(nóng)田中可被作物利用的水資源量，反映了糧食生產(chǎn)過(guò)程中降水資源的利用情況。糧食總水足跡為糧食藍(lán)水足跡和綠水足跡的總和。本文中涉及的糧食藍(lán)水足跡、綠水足跡及總水足跡均為3種作物藍(lán)水足跡、綠水足跡及總水足跡平均值。糧食水足跡的具體公式為

wf=wfb+wfg

(3)

其中

(4)

式中wf——糧食總水足跡，m3/t

wfb——糧食藍(lán)水足跡，m3/t

wfg——糧食綠水足跡，m3/t

Ig——糧食作物生育期內(nèi)的毛灌溉水量，mm

Pe——糧食作物生育期內(nèi)的有效降水量，mm

y——單位面積糧食產(chǎn)量，t/hm2

其中10是常量因子，是將水深轉(zhuǎn)換為單位面積水量的轉(zhuǎn)換系數(shù)。

圖3 糧食綠水足跡、藍(lán)水足跡及總水足跡空間分布圖Fig.3 Spatial distribution maps of crop green, blue and total water footprints

2 結(jié)果與分析

2.1 黑龍江省糧食水足跡時(shí)空分布規(guī)律

由圖2可知，2008—2018年間黑龍江省糧食總水足跡整體呈波動(dòng)下降趨勢(shì)：2008—2010年間糧食總水足跡呈緩慢上升趨勢(shì)，而2011年驟降至1 750 m3/t后，一直處在2 000 m3/t上下，比2008—2010年的糧食總水足跡平均值下降20%，直至2018年恢復(fù)上升趨向。由趨勢(shì)線可以看出，糧食總水足跡的趨勢(shì)線斜率為-27.487，其下降趨勢(shì)最為顯著。造成黑龍江省糧食總水足跡下降的原因包括兩方面：①自然環(huán)境因素。由于全球氣候變暖對(duì)糧食種植產(chǎn)生了較大影響，降水減少而蒸發(fā)加快，降低了可用于糧食生產(chǎn)的有效降水量，使得糧食總水足跡下降[36]。②社會(huì)發(fā)展因素。隨著工業(yè)制造業(yè)用水增多，農(nóng)業(yè)用水逐漸被擠占，促使農(nóng)業(yè)現(xiàn)代化程度提高，水資源利用效率得到提升，使得糧食總水足跡呈下降趨勢(shì)。

圖2 2008—2018年黑龍江省糧食綠水足跡、藍(lán)水 足跡及總水足跡時(shí)間分布規(guī)律Fig.2 Temporal distribution of crop green, blue and total water footprints in Heilongjiang Province in 2008—2018

由圖3可知，黑龍江省糧食總水足跡整體呈現(xiàn)南北高，東西低的變化趨勢(shì)，其中，糧食綠水足跡及藍(lán)水足跡的空間分布趨勢(shì)具有一致性，糧食藍(lán)水足跡、綠水足跡高值區(qū)在黑龍江省的牡丹江、伊春和齊齊哈爾，而哈爾濱、大慶等中東部地區(qū)屬糧食藍(lán)水足跡、綠水足跡低值區(qū)。糧食總水足跡空間分布規(guī)律與降水量有密切聯(lián)系，黑龍江省南、北部和西部降水量較多，屬于濕潤(rùn)地區(qū)，糧食總水足跡高；而中部和東部地區(qū)降水量偏少，屬于半濕潤(rùn)半干旱地區(qū)，糧食總水足跡較低[37]。

2.1.1時(shí)間分布

由圖2可知，黑龍江省糧食藍(lán)水足跡、綠水足跡的變化趨勢(shì)與糧食總水足跡一致：2008—2010年糧食藍(lán)水足跡、綠水足跡呈上升趨勢(shì)，自2011年驟降后，期間雖有波動(dòng)但未能恢復(fù)，直至2018年才有所上升。從糧食藍(lán)水足跡、綠水足跡的對(duì)比情況來(lái)看，糧食綠水足跡趨勢(shì)線斜率為-14.422，糧食藍(lán)水足跡趨勢(shì)線斜率為-13.065，糧食綠水足跡下降趨勢(shì)比藍(lán)水足跡顯著；糧食藍(lán)水足跡與綠水足跡比例始終趨于1∶3，表明黑龍江省糧食用水結(jié)構(gòu)10年間基本保持不變。

2008—2018黑龍江省糧食綠水足跡波動(dòng)在1 750 m3/t，整體有下降趨勢(shì)。由于2011年黑龍江省糧食作物生產(chǎn)經(jīng)歷了春早、伏旱等氣象災(zāi)害，使得當(dāng)年黑龍江省糧食綠水足跡呈現(xiàn)顯著下降，之后的幾年持續(xù)波動(dòng)始終未超過(guò)1 750 m3/t，因此將2008—2018年黑龍江省糧食綠水足跡分成兩個(gè)時(shí)期：2008—2010年(糧食綠水足跡平均值1 838 m3/t)和2011—2018年(糧食綠水足跡平均值1 562 m3/t)。黑龍江省糧食藍(lán)水足跡也有著類似的情況，自然災(zāi)害同樣影響到2011年黑龍江省糧食藍(lán)水足跡，由于旱災(zāi)使得降水量減少，導(dǎo)致可灌溉供水量的驟降，進(jìn)而引起糧食藍(lán)水足跡下降。

2.1.2空間分布

利用ArcGIS軟件繪制2008—2018年黑龍江省11個(gè)地區(qū)糧食綠水足跡、藍(lán)水足跡及總水足跡的空間分布圖，如圖4～6所示。黑龍江省綠水足跡的低值區(qū)主要集中在黑龍江省的中部和東部地區(qū)，包括大慶、哈爾濱、雞西、雙鴨山、佳木斯、綏化，其中大慶最低，占總綠水足跡的5.3%。綠水足跡的高值區(qū)主要分布在南部、北部和西部地區(qū)，包括伊春、牡丹江、七臺(tái)河、鶴崗市、齊齊哈爾，其中伊春市最高，占糧食總綠水足跡比例為15.2%。黑龍江省藍(lán)水足跡的低值區(qū)主要集中在黑龍江省的中部和東部地區(qū)，特別是大慶、哈爾濱、綏化和佳木斯。而高值區(qū)則分布在南部、北部和西部地區(qū)，包括齊齊哈爾、伊春、牡丹江，北部地區(qū)齊齊哈爾的糧食藍(lán)水足跡最高。黑龍江省糧食總水足跡整體呈現(xiàn)南北高、東西低的變化趨勢(shì)。為了更清晰地表征黑龍江省各個(gè)地區(qū)的糧食總水足跡等級(jí)分布情況，將糧食總水足跡分為4個(gè)等級(jí)：伊春、牡丹江糧食總水足跡為Ⅰ級(jí)，糧食總水足跡過(guò)剩；齊齊哈爾、七臺(tái)河和鶴崗為Ⅱ級(jí)，糧食總水足跡較為適宜；雙鴨山、佳木斯、雞西、綏化為Ⅲ級(jí)，糧食總水足跡相對(duì)正常；而大慶地區(qū)、哈爾濱地區(qū)為Ⅳ級(jí)，糧食總水足跡十分稀缺。由此可以看出，各個(gè)地區(qū)水足跡分布不均，差異較大，存在著糧食水足跡過(guò)剩與稀缺并存的情況?？刹扇」こ檀胧?shí)現(xiàn)省內(nèi)的“南北水，東西調(diào)”，將糧食水足跡過(guò)剩地區(qū)如伊春、牡丹江的水資源，調(diào)入水足跡稀缺地區(qū)如哈爾濱、大慶，以改善糧食水足跡的空間分布結(jié)構(gòu)。

2.2 黑龍江省糧食水足跡影響因素分析

利用斯皮爾曼相關(guān)系數(shù)法對(duì)黑龍江省糧食水足跡的影響因素進(jìn)行分析，得到各影響因素相關(guān)系數(shù)(表1)。對(duì)于糧食藍(lán)水足跡，人口密度對(duì)其影響最大，其次是GDP、播種面積、降水量、糧食產(chǎn)量、單位面積產(chǎn)量，其中，GDP、單位面積產(chǎn)量、糧食產(chǎn)量、播種面積的相關(guān)系數(shù)為負(fù)值。GDP的相關(guān)系數(shù)為-0.573 2，負(fù)相關(guān)程度最顯著。糧食藍(lán)水足跡反映了糧食生產(chǎn)過(guò)程中消耗的灌溉水資源，在灌溉用水比例一定的情況下，糧食藍(lán)水足跡越低，灌溉水資源的利用效率越高。GDP是指一定時(shí)期內(nèi)生產(chǎn)活動(dòng)的最終成果，代表著經(jīng)濟(jì)發(fā)展水平，一個(gè)地區(qū)的經(jīng)濟(jì)發(fā)展水平越高，工業(yè)及第三產(chǎn)業(yè)越發(fā)達(dá)，相應(yīng)的用水量越大，導(dǎo)致農(nóng)業(yè)可供水被工業(yè)及第三產(chǎn)業(yè)用水?dāng)D占，因而GDP與糧食藍(lán)水足跡呈現(xiàn)負(fù)相關(guān)；人口密度對(duì)糧食藍(lán)水足跡影響最大，相關(guān)系數(shù)為0.718 2，表明人口越多，糧食產(chǎn)品的需求量越高，生產(chǎn)過(guò)程中消耗的農(nóng)業(yè)灌溉水資源總量就越多；此外，降水量可間接影響灌溉可供水量，降水量越大，灌溉可供水量越大，因而降水量與糧食藍(lán)水足跡呈正相關(guān)。

圖4 糧食綠水足跡逐年空間分布圖Fig.4 Spatial distribution maps of crop green water footprint

圖5 糧食藍(lán)水足跡逐年空間分布圖Fig.5 Spatial distribution maps of crop blue water footprint

圖6 糧食總水足跡逐年空間分布圖Fig.6 Spatial distribution maps of crop total water footprint

對(duì)于糧食綠水足跡，降水量對(duì)其影響最大，其次是播種面積、GDP、糧食產(chǎn)量、單位面積產(chǎn)量、人口密度，其中，GDP、單位面積產(chǎn)量、糧食產(chǎn)量、播種面積與糧食綠水足跡呈負(fù)相關(guān)。糧食綠水足跡反映了糧食生產(chǎn)過(guò)程中降水資源的可利用情況，由于降水量直接影響作物生長(zhǎng)蒸騰作用中可利用雨水資源量，其對(duì)糧食綠水足跡影響最大，相關(guān)系數(shù)為0.672 7。播種面積與糧食綠水足跡的負(fù)相關(guān)性最顯著，相關(guān)系數(shù)為-0.263 6，在作物生育期有效降水量穩(wěn)定的狀態(tài)下，隨著播種面積的增大，單位播種面積上的可利用降水量減少，進(jìn)而導(dǎo)致糧食綠水足跡下降。

對(duì)于糧食總水足跡，GDP、播種面積、糧食產(chǎn)量、降水量與其呈負(fù)相關(guān)，單位面積產(chǎn)量、人口密度與糧食總水足跡呈正相關(guān)。根據(jù)糧食水足跡的計(jì)算公式可知，單位面積產(chǎn)量是計(jì)算糧食藍(lán)水足跡、綠水足跡及總水足跡量的重要變量之一，因而單位面積產(chǎn)量的影響最大，相關(guān)系數(shù)為0.663 6；而GDP越高，農(nóng)業(yè)工業(yè)化程度越大，水資源利用效率越高，使得糧食總水足跡呈下降趨勢(shì)。因而GDP與糧食總水足跡的負(fù)相關(guān)性顯著。

2.3 黑龍江省糧食水足跡預(yù)測(cè)結(jié)果分析

結(jié)合當(dāng)?shù)匾?guī)劃目標(biāo)，根據(jù)式(3)、(4)及2.2節(jié)的結(jié)果選取對(duì)糧食水足跡影響較大的3個(gè)因素進(jìn)行預(yù)測(cè)，即播種面積、降水量、糧食產(chǎn)量(1998—2018年)，得出2021—2025年糧食水足跡3個(gè)因素的預(yù)測(cè)趨勢(shì)圖(圖7)。由圖可以看出，糧食產(chǎn)量和播種面積在未來(lái)5年整體都有上升趨勢(shì)，而降水量波動(dòng)幅度較大。根據(jù)以上預(yù)測(cè)結(jié)果對(duì)糧食藍(lán)水足跡、綠水足跡及總水足跡進(jìn)行計(jì)算，由表2可知，2021—2025年黑龍江省糧食綠水足跡遠(yuǎn)大于糧食藍(lán)水足跡，糧食藍(lán)水足跡平均值為563 m3/t，糧食綠水足跡平均值為1 487 m3/t，糧食總水足跡平均值為2 050 m3/t，糧食藍(lán)水足跡、綠水足跡、總水足跡均有下降趨勢(shì)。

結(jié)合上述結(jié)果分析造成糧食藍(lán)水足跡、綠水足跡及總水足跡下降趨勢(shì)的原因：①氣候環(huán)境的變化使降水量減少，同時(shí)氣候變暖使蒸發(fā)加快，因而糧食藍(lán)水足跡、綠水足跡及總水足跡下降。②根據(jù)當(dāng)前的政策形勢(shì)來(lái)看，農(nóng)業(yè)逐漸工業(yè)化，形成規(guī)模效應(yīng)。雖然播種面積和產(chǎn)量在未來(lái)幾年略有上升，預(yù)計(jì)2021—2025年播種面積從8.15×106hm2增長(zhǎng)到8.74×106hm2，漲幅為7.24%；產(chǎn)量由3.319×107t增長(zhǎng)到3.474×107t，漲幅為4.67%，然而工業(yè)等其他行業(yè)用水?dāng)D占農(nóng)業(yè)用水的現(xiàn)象凸顯[38]，農(nóng)業(yè)工業(yè)化程度越大，水資源利用效率越高，進(jìn)而糧食藍(lán)水足跡、綠水足跡及總水足跡呈現(xiàn)降低趨勢(shì)?？刹扇】茖W(xué)方法合理制定農(nóng)業(yè)發(fā)展規(guī)劃，持續(xù)提高農(nóng)業(yè)水資源利用率，合理促進(jìn)農(nóng)業(yè)發(fā)展。

3 結(jié)論

(1)黑龍江省不同地區(qū)糧食水足跡差異顯著，>存在糧食水足跡過(guò)剩與稀缺并存的現(xiàn)象。糧食水足跡在空間上呈現(xiàn)南北高、東西低的趨勢(shì)，伊春、牡丹江地區(qū)的糧食水足跡較豐富，而哈爾濱、大慶地區(qū)的糧食水足跡相對(duì)匱乏；在時(shí)間分布上，2011、2016年糧食總水足跡最低。

(2)人口密度對(duì)糧食藍(lán)水足跡影響最大，降水量對(duì)糧食綠水足跡影響最大，GDP、播種面積、糧食產(chǎn)量、降水量與糧食總水足跡呈負(fù)相關(guān)，單位面積產(chǎn)量、人口密度與糧食總水足跡呈正相關(guān)。

(3)未來(lái)5年(2021—2025年)，黑龍江省糧食藍(lán)水足跡平均值為563 m3/t，糧食綠水足跡平均值為1 487 m3/t，糧食總水足跡平均值為2 050 m3/t，糧食藍(lán)水足跡、綠水足跡及總水足跡呈下降趨勢(shì)。

(4)黑龍江省糧食生產(chǎn)供水量不足的問(wèn)題不容樂(lè)觀，一方面應(yīng)科學(xué)制定農(nóng)業(yè)發(fā)展規(guī)劃，提高農(nóng)業(yè)水資源利用率，合理促進(jìn)農(nóng)業(yè)發(fā)展；另一方面應(yīng)采取工程措施，實(shí)現(xiàn)省內(nèi)的“南北水，東西調(diào)”，將糧食水足跡過(guò)剩地區(qū)的水資源調(diào)入水足跡稀缺地區(qū)，以改善糧食水足跡的空間分布結(jié)構(gòu)。