付 強 劉 燁 李天霄 崔 嵩 劉 東 成 琨

(東北農(nóng)業(yè)大學水利與土木工程學院， 哈爾濱 150030)

水足跡視角下黑龍江省糧食生產(chǎn)用水分析

付 強 劉 燁 李天霄 崔 嵩 劉 東 成 琨

(東北農(nóng)業(yè)大學水利與土木工程學院， 哈爾濱 150030)

運用水足跡理論，計算黑龍江省各地區(qū)2007—2012年的糧食水足跡，分析其糧食作物用水情況。結(jié)果表明：黑龍江省糧食生產(chǎn)用水量巨大，糧食作物灌溉用水量存在逐年增加的趨勢，用于糧食生產(chǎn)的灌溉用水量占黑龍江省水資源消耗量的絕大多數(shù)。黑龍江省各地區(qū)的糧食水足跡在空間分布上存在顯著差異，大體上呈東北部高，西南部低的趨勢；糧食藍水足跡明顯低于綠水足跡，糧食生產(chǎn)對綠水的依賴度較高。黑龍江省糧食水足跡整體上小于單產(chǎn)需水量，部分地區(qū)糧食作物生產(chǎn)過程中缺水現(xiàn)象嚴重?？臻g自相關(guān)分析顯示，黑龍江省2007—2012年的糧食綠水足跡和藍水足跡均存在相似值之間的空間聚集性。黑龍江省農(nóng)墾總局的糧食產(chǎn)量較高，但糧食藍水足跡比例較高，糧食生產(chǎn)比較依賴灌溉水資源，水資源消耗量較大。由于2009年糧食作物單位面積產(chǎn)量顯著降低，使得2007—2012年被分成2007—2008年和2009—2012年2個不同的時期，糧食水足跡在2個時期均存在逐年降低的趨勢，糧食生產(chǎn)的用水效率逐年升高。

糧食； 水足跡； 水資源； 空間分布； 灌溉； 黑龍江省

引言

水足跡是一個國家、一個地區(qū)或一個人在一定時期內(nèi)消費的產(chǎn)品和服務(wù)所需要的淡水資源的數(shù)量[1]。它也可以表示生產(chǎn)某一產(chǎn)品或提供某種服務(wù)所消耗的淡水資源量[2]，它是虛擬水理論的繼承和發(fā)展[3]。糧食作為一種水資源密集型產(chǎn)品，其水足跡能夠有效地反映某一區(qū)域內(nèi)糧食生產(chǎn)過程中的水資源利用情況。對于糧食水足跡的研究既有灌區(qū)尺度研究[4]，也有以全國各省區(qū)為研究對象的全國尺度研究[5]，還有全球尺度研究[6]。

黑龍江省作為糧食大省，當前對其糧食生產(chǎn)的研究較多，主要是基于全省尺度的糧食產(chǎn)量研究[7-8]和全省尺度的糧食生產(chǎn)用水研究[9]。與以用水量為主導的水資源評價和水資源利用分析視角不同，王斌[10]利用產(chǎn)量和水分利用效率估算糧食總體耗水量的方法分析了黑龍江省2002—2011年的糧食耗水量，但同樣只能從整體上分析黑龍江省糧食生產(chǎn)用水狀況，并沒有對其在黑龍江省各地區(qū)的空間分布進行分析。操信春等[5]采用水足跡對全國各省的糧食生產(chǎn)用水進行了分析和評價，分析了各省糧食生產(chǎn)用水的空間分布差異。而在黑龍江省省區(qū)尺度，利用水足跡理論分析糧食生產(chǎn)用水的研究并不多，尚不能滿足黑龍江省糧食生產(chǎn)用水管理的需要。

本文分析黑龍江省的糧食生產(chǎn)用水概況，并運用水足跡理論，分別對黑龍江省各地區(qū)2007—2012年的糧食藍水足跡和綠水足跡進行計算，分析黑龍江省糧食水足跡的時間變化規(guī)律和空間分布特征，在水足跡的視角下分析糧食作物生產(chǎn)過程中的水資源利用情況，比較糧食水足跡與單產(chǎn)需水量的差異，為黑龍江省各地區(qū)糧食生產(chǎn)的用水管理提供理論依據(jù)，為黑龍江省可用水資源的分配提供技術(shù)支持。

1 材料與方法

1.1 研究區(qū)域概況

黑龍江省耕地面積占全國耕地面積的1/9，糧食產(chǎn)量排在全國第2位，糧食商品量占全國的1/6，是我國最大的商品糧生產(chǎn)基地[11-12]。糧食作物生產(chǎn)用水作為黑龍江省主要的水資源消耗方式之一，嚴重制約著黑龍江省的可用水資源分配。黑龍江省糧食作物生產(chǎn)用水的高效利用，關(guān)系到國家的糧食安全、水資源安全以及生態(tài)安全[13]。因此，對于黑龍江省糧食生產(chǎn)中水資源利用情況的研究具有非常巨大的現(xiàn)實意義。黑龍江省除了下轄13個地級市外，還存在著農(nóng)墾總局這一特殊的行政管理體制。黑龍江省農(nóng)墾總局下轄9個管理局、113個農(nóng)牧場，分布在全省12個市，遍布小興安嶺山麓、松嫩平原和三江平原地區(qū)。

1.2 數(shù)據(jù)來源

文中的糧食是谷物、豆類和薯類等的統(tǒng)稱[14]。黑龍江省各地區(qū)降水量數(shù)據(jù)來自中國氣象數(shù)據(jù)網(wǎng)(http:∥data.cma.cn)。農(nóng)業(yè)用水數(shù)據(jù)來自《黑龍江省水利建設(shè)統(tǒng)計年鑒》[15]。糧食產(chǎn)量、播種面積等農(nóng)業(yè)數(shù)據(jù)來自《黑龍江省統(tǒng)計年鑒》[16]。

1.3 計算方法

1.3.1 糧食水足跡

研究區(qū)域為黑龍江省的13個地級市行政區(qū)，根據(jù)糧食生產(chǎn)過程中的實際用水量來計算糧食水足跡。為了避免不同年份間的差異對研究結(jié)果的影響，同時分析糧食水足跡的時間變化規(guī)律，分別計算2007—2012年黑龍江省各地區(qū)的糧食水足跡。糧食水足跡包括糧食藍水足跡和糧食綠水足跡，其計算公式[5]分別為

wf=wfb+wfg

(1)

其中

(2)

式中wf——單位產(chǎn)量的糧食水足跡，m3/kgwfb——單位產(chǎn)量的糧食藍水足跡，m3/kgwfg——單位產(chǎn)量的糧食綠水足跡，m3/kgIg——糧食作物生育期內(nèi)單位面積的毛灌溉水量，mm

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

y——單位面積的糧食作物產(chǎn)量，kg/hm2

糧食藍水足跡比例計算式為

(3)

糧食作物毛灌溉總水量計算式為

(4)

式中Wf——糧食作物總用水量，m3Wa——農(nóng)作物總用水量，m3Sfi——第i種糧食作物種植面積，hm2Saj——第j種農(nóng)作物種植面積，hm2mfi——第i種糧食作物灌溉定額，m3/hm2maj——第j種農(nóng)作物灌溉定額，m3/hm2n——糧食作物種類數(shù)l——農(nóng)作物種類數(shù)

黑龍江省2007—2012年6年間，糧食作物種植面積均占農(nóng)作物種植面積的90%以上，雖然各糧食作物和非糧食作物的灌溉定額差異較大，但非糧食作物的種植面積遠小于糧食作物種植面積，因此將糧食作物和非糧食作物的灌溉定額假設(shè)相同，其計算不會產(chǎn)生較大的誤差[17]。

黑龍江省糧食作物的生育期主要分布在5—9月份[18]，在計算糧食綠水足跡時，生育期內(nèi)有效降水量取5—9月份的有效降水量之和。

有效降水量的計算采用美國農(nóng)業(yè)部土壤保持局推薦的方法[13,19]，以旬為步長，計算式為

(5)

其中

(6)

式中Pei——糧食作物生育期內(nèi)第i旬的有效降水量，mm

Pi——糧食作物生育期內(nèi)第i旬降水量，mm

k——糧食作物生育期內(nèi)總旬數(shù)

1.3.2 空間自相關(guān)性分析的全局Moran’sI

全局自相關(guān)是反映全部研究對象間的相關(guān)性，Moran’sI是一個常用的度量全局空間自相關(guān)的指標[20]。Moran’sI反映空間鄰接或空間鄰近的區(qū)域單元屬性值的相似程度，范圍在-1～1之間，1表示空間數(shù)據(jù)之間存在強烈的正自相關(guān)，-1表示空間數(shù)據(jù)之間存在強烈的負自相關(guān)，0則表示空間數(shù)據(jù)是隨機模式，獨立分布，不存在相關(guān)性[21]。Moran’sI的計算公式[22]為

(7)

式中Wij——空間權(quán)重矩陣元素，通過空間鄰接方式和空間距離來確定

xi——區(qū)域單元i的屬性值

t——區(qū)域單元數(shù)

常用標準化統(tǒng)計量Z來檢驗區(qū)域是否存在空間自相關(guān)關(guān)系，Z的計算公式為

(8)

式中E(I)——Moran’sI的期望值 Var(I)——Moran’sI的方差

當Z值為正時，表明存在正的空間自相關(guān)，也就是說相似的觀測值趨于空間集聚；當Z值為負時，表明存在負的空間自相關(guān)，相似的觀測值趨于分散分布；當Z值為0時，觀測值呈現(xiàn)獨立隨機分布。

1.3.3 糧食作物單產(chǎn)需水量

各地區(qū)的糧食作物單產(chǎn)需水量計算式[23]為

(9)

其中

(10)

式中wr——主要糧食作物加權(quán)平均得到的單產(chǎn)需水量，m3/kg

wri——第i種糧食作物單產(chǎn)需水量，m3/kg

Wri——第i種糧食作物單位面積需水量，可看作是糧食作物生育期的累積蒸發(fā)蒸騰總量，可由彭曼公式計算得到[24]，m3/hm2

yi——第i種糧食作物單位面積產(chǎn)量，kg/hm2

2 結(jié)果與分析

2.1 黑龍江省糧食生產(chǎn)及用水分析

圖1 黑龍江省各地區(qū)糧食生產(chǎn)及用水指標分布Fig.1 Distributions of grain production and water utilization in different areas of Heilongjiang Province

2007—2012年黑龍江省各地區(qū)糧食生產(chǎn)及用水指標分布情況如圖1所示?？臻g上，黑龍江省各地區(qū)糧食作物種植面積的分布差異較大(圖1a)，農(nóng)墾總局、齊齊哈爾市、哈爾濱市、綏化市和佳木斯市等地區(qū)的糧食作物種植面積較大；而大興安嶺地區(qū)、七臺河市、鶴崗市和伊春市等地區(qū)的糧食作物種植面積較小。黑龍江省各地區(qū)糧食作物總產(chǎn)量的分布規(guī)律(圖1b)和糧食作物種植面積的分布規(guī)律基本一致。

黑龍江省各地區(qū)糧食作物灌溉用水量分布差異也較大(圖1c)，農(nóng)墾總局最大，占全省糧食作物灌溉用水量的30.53%，哈爾濱市、佳木斯市、綏化市和齊齊哈爾市等地區(qū)也均較大；而大興安嶺地區(qū)、黑河市和七臺河市等地區(qū)均較小。黑龍江省各地區(qū)糧食作物用水比例分布差異不明顯(圖1d)，黑龍江省農(nóng)墾總局糧食作物用水比例最高；大興安嶺地區(qū)最低；由于石油工業(yè)存在著較大的耗水量，大慶市也明顯低于全省平均值。

黑龍江省各地區(qū)的灌溉面積比例的分布差異顯著(圖1e)，其中，大慶市最高，農(nóng)墾總局也較高；大興安嶺地區(qū)、黑河市最低；而鶴崗市、佳木斯市6年間存在明顯的增加趨勢。黑龍江省各地區(qū)的水田面積比例的分布差異也較為顯著(圖1f)，其中，農(nóng)墾總局最高，鶴崗市、雞西市、佳木斯市、哈爾濱市均較高；大興安嶺地區(qū)沒有水田種植，黑河市的水田種植面積極低。

黑龍江省2007—2012年糧食作物生產(chǎn)及用水總體概況如表1所示。時間上，2009年相比2008年糧食播種面積呈現(xiàn)顯著增加的趨勢，而其他年份則呈平穩(wěn)增加的趨勢。黑龍江省灌溉面積比例6年間呈上升趨勢，僅2009年呈下降趨勢，這主要是由于2009年糧食播種面積顯著增加，而灌溉面積沒有明顯增加，導致灌溉面積比例在2009年呈現(xiàn)一個較大的下降趨勢；同時2009年也是豐水年，降水量較大，對灌溉的需求量較小。黑龍江省水田面積比例6年間同樣呈現(xiàn)平穩(wěn)增加趨勢。黑龍江省糧食作物用水比例6年間變化幅度很小，6年平均值為86.27%。通過糧食作物用水比例可以看出，黑龍江省作為一個農(nóng)業(yè)大省，糧食作物消耗的灌溉水資源量占黑龍江省水資源消耗的絕大多數(shù)，較高的糧食作物用水比例嚴重制約著黑龍江省除了農(nóng)業(yè)以外其他產(chǎn)業(yè)的發(fā)展，阻礙了黑龍江省經(jīng)濟的良性發(fā)展。在2007—2012年的6年間，黑龍江省糧食作物的播種面積逐年增長，糧食作物總產(chǎn)量也穩(wěn)步上升，同時灌溉面積比例大體上也呈逐年增加的趨勢，糧食生產(chǎn)消耗的灌溉水量也隨之逐年增加，導致黑龍江省的水資源壓力逐年增加。而糧食作物灌溉水資源投入的增加對于糧食作物生產(chǎn)過程中灌溉水資源的利用效率具有一定的影響，這種影響關(guān)系到黑龍江省可用水資源分配的決策，這就需要從糧食水足跡的角度來分析水資源的利用效率。

表1 黑龍江省糧食作物生產(chǎn)及用水總體概況

圖2 黑龍江省各地區(qū)糧食水足跡分布Fig.2 Distributions of grain water footprint in different areas of Heilongjiang Province

2.2 糧食水足跡分析

2.2.1 空間分析

運用水足跡理論，分別計算黑龍江省各地區(qū)的糧食藍水足跡和綠水足跡，并計算黑龍江省各地區(qū)的糧食水足跡和藍水比例，分析黑龍江省糧食作物灌溉過程中灌溉水資源和降水資源的利用情況。計算得到黑龍江省13個地區(qū)2007—2012年的糧食水足跡空間分布狀況如圖2所示。

糧食藍水足跡反映了糧食生產(chǎn)過程中消耗的灌溉水資源，在灌溉水比例一定的情況下，糧食藍水足跡越低，灌溉水資源的利用效率越高。空間上，黑龍江省各地區(qū)糧食藍水足跡的分布差異顯著(圖2a)，大體上呈北部低、南部高、東部高、西部低的趨勢。其中，鶴崗市最高，最高值為全省平均值的3.57倍，雞西市、佳木斯市、哈爾濱市均較高；大興安嶺地區(qū)最低，最低值僅為全省平均值的2%，其次為黑河市；黑龍江省農(nóng)墾總局6年的糧食藍水足跡均高于全省平均水平。

糧食綠水足跡反映了糧食生產(chǎn)過程中降水資源的利用情況，黑龍江省各地區(qū)糧食綠水足跡的分布具有一定的差異性(圖2b)，大體上呈北部高、南部低、東部高、西部低的趨勢。其中，大興安嶺地區(qū)最高，最高值為全省平均值的2.33倍，其次為黑河市；大慶市最低，最低值僅為全省平均值的44%，哈爾濱市、綏化市、雞西市、雙鴨山市均較低；黑龍江省農(nóng)墾總局6年的糧食綠水足跡含量均低于全省平均水平。

空間上，黑龍江省各地區(qū)的糧食水足跡的分布差異顯著(圖2c)，全省大體上呈北部高、南部低、東部高、西部低的趨勢。其中，鶴崗市的糧食水足跡最高，最高值為全省平均值的1.99倍，大興安嶺地區(qū)、黑河市均較高；綏化市最低，最低值僅為全省平均值的54%，大慶市也較低；黑龍江省農(nóng)墾總局6年的糧食水足跡接近全省平均水平。

藍、綠水足跡不足以表示糧食生產(chǎn)用水的組成狀況，而藍水比例可以反映不同區(qū)域水足跡的構(gòu)成存在差異。黑龍江省各地區(qū)的糧食水足跡藍水比例的分布差異顯著(圖2d)，大體上呈北部低、南部高、東部高、西部低的趨勢。絕大多數(shù)地區(qū)6年平均藍水比例均在50%以下，黑龍江省的灌溉系統(tǒng)與中國南方省區(qū)相比較為不發(fā)達，糧食生產(chǎn)比較依賴降水資源。藍水比例略高于50%的地區(qū)僅有三江平原的雞西市和鶴崗市，這主要是由于水田面積比例較高，而水稻生產(chǎn)的耗水量明顯大于其他糧食作物[25]，因此應(yīng)適當減少水田面積或采用先進的節(jié)水灌溉技術(shù)，降低糧食灌溉水引用量，緩解該地區(qū)的水資源壓力。而同處于三江平原的佳木斯市、七臺河市、雙鴨山市等地區(qū)的藍水比例相對較低。主要處于松嫩平原的哈爾濱市、大慶市、齊齊哈爾市、綏化市等地區(qū)的糧食水足跡均較低，藍水比例處于全省中等水平，灌溉水和降水的利用效率均較高，結(jié)合這些地區(qū)較高的糧食作物單產(chǎn)水平，較為適合進行糧食生產(chǎn)。藍水比例最低的大興安嶺地區(qū)和黑河市均不足5%，糧食生產(chǎn)用水幾乎完全依賴于降水資源。且兩地的糧食單產(chǎn)不高，不適合發(fā)展灌溉農(nóng)業(yè)生產(chǎn)，兩地應(yīng)控制本地區(qū)水田的面積，以達到水資源可持續(xù)發(fā)展的需要；同時通過適當?shù)挠盟呗?，緩解可用水資源短缺的情況，促進區(qū)域經(jīng)濟發(fā)展。牡丹江市和伊春市的糧食綠水足跡均處于中等水平，而糧食藍水足跡較低，藍水比例較低。農(nóng)墾總局的藍水比例高于全省平均水平，在全省屬于較高的水平。由于耕地規(guī)模大、現(xiàn)代化程度高、綜合生產(chǎn)能力強，黑龍江省農(nóng)墾總局僅用全省19.94%的糧食作物種植面積生產(chǎn)出全省35.40%的糧食作物產(chǎn)量，糧食作物生產(chǎn)能力較強。但農(nóng)墾總局的糧食作物生產(chǎn)比較依賴灌溉水資源，灌溉水資源的消耗量較大，占全省糧食灌溉用水量的30.53%。這勢必會增加周邊地區(qū)的水資源壓力，不利于周邊地區(qū)的發(fā)展。

2.2.2 空間自相關(guān)性分析

對2007—2012年黑龍江省各地區(qū)的糧食藍水足跡和綠水足跡進行空間自相關(guān)分析，得到各年糧食水足跡的Moran’sI值如表2所示。黑龍江省2007—2012年6年間的綠水足跡空間分布存在著顯著的正的空間自相關(guān)性，即黑龍江省各地區(qū)糧食綠水足跡空間分布并非表現(xiàn)出完全的隨機性，而是表現(xiàn)出相似值之間的空間集聚，糧食綠水足跡高的地區(qū)相對地趨于和糧食綠水足跡高的地區(qū)相鄰，糧食綠水足跡低的地區(qū)相對地趨于和糧食綠水足跡低的地區(qū)相鄰。這主要是由于糧食綠水足跡主要受降水量的影響，而降水量在地理位置相鄰近的地區(qū)空間上呈相似集聚[26]。2009—2011年黑龍江省的藍水足跡空間分布呈現(xiàn)為顯著的、正的空間自相關(guān)；而2008年和2012年的藍水足跡空間分布呈現(xiàn)為不顯著的正的空間自相關(guān)性，只有2007年的藍水足跡空間分布表現(xiàn)為不顯著的負的空間自相關(guān)性，且Moran’sI值只有-0.020。6年間藍水足跡的空間分布整體上也呈現(xiàn)為正的空間自相關(guān)性，表現(xiàn)為相似值之間的空間集聚，糧食藍水足跡高的地區(qū)相對地趨于和糧食藍水足跡高的地區(qū)相鄰，糧食藍水足跡低的地區(qū)相對地趨于和糧食藍水足跡低的地區(qū)相鄰。這主要是由于黑龍江省各相鄰地區(qū)的區(qū)域灌溉特征、糧食作物生產(chǎn)管理方式和糧食生產(chǎn)投入狀況等社會經(jīng)濟條件在地理空間上呈現(xiàn)出一定的相似性[5]。黑龍江省各地區(qū)的藍水足跡雖然差異明顯，但在相鄰區(qū)域具有一定的相似性，黑龍江省各地區(qū)應(yīng)根據(jù)自身實際情況，借鑒其相鄰地區(qū)的糧食作物種植策略和灌溉管理方式經(jīng)驗，因地制宜地發(fā)展符合地區(qū)特點的糧食生產(chǎn)方式。

2.2.3 時間分析

黑龍江省各地區(qū)2007—2012年糧食水足跡的全省平均值如表3所示。2007—2012年黑龍江省糧食藍水足跡整體上呈現(xiàn)逐年下降的趨勢。糧食藍水足跡在2009年呈現(xiàn)出異常的上升趨勢，這主要是由于雖然2009年相比2008年糧食播種面積激增，但是2009年黑龍江省農(nóng)業(yè)作物生產(chǎn)經(jīng)歷了春旱、夏澇、伏旱、寡照和低溫等多重自然災(zāi)害[27]，使得糧食作物的單位面積產(chǎn)量顯著降低，最終導致2009年的糧食藍水足跡異常增多，同時也使得2007—2012年被分成2007—2008年和2009—2012年2個不同的時期。

表2 黑龍江省糧食水足跡Moran’s I及其統(tǒng)計參數(shù)值

注：p為假設(shè)概率。

表3 黑龍江省各地區(qū)糧食水足跡平均值

黑龍江省的糧食綠水足跡在2007—2008年和2009—2012年2個不同的時期分別呈逐年下降趨勢，而在2009年有顯著的增加。2009年糧食作物單位面積產(chǎn)量的顯著降低同樣也影響到了2009年的糧食綠水足跡，導致2009年糧食綠水足跡增多；另一方面，糧食綠水足跡主要來源于作物生育期的有效降水量，由于2009年是一個典型的豐水年，降水量異常豐富也是導致2009年糧食綠水足跡較2008年增多較多的一個主要原因。2007—2012年黑龍江省各地區(qū)糧食作物生育期(5—9月份)平均降水量分別為271.66、296.16、355.74、324.62、300.33、338.50 mm。2008年的生育期平均降水量高于2007年，2008年的糧食綠水足跡卻較2007年低；同樣2012年的降水量高于2010年和2011年，糧食綠水足跡在2009—2012年卻呈逐年下降的趨勢。

黑龍江省糧食水足跡在2007—2008年呈現(xiàn)減小的趨勢，2009—2012年的糧食水足跡同樣呈現(xiàn)逐年減小的趨勢，而2009年的糧食水足跡由于藍水足跡和綠水足跡均異常增多而較2008年呈增加趨勢。在2007—2008年和2009—2012年2個時期內(nèi)，糧食作物灌溉用水量的逐年增加雖然消耗了更多的水資源，但同時也生產(chǎn)了更多的糧食，且生產(chǎn)單位質(zhì)量的糧食消耗的水資源量減少，獲得的糧食作物增產(chǎn)的增益高于水資源消耗的增加。節(jié)約水資源并不僅是要節(jié)約水資源的用量，更要提高水資源的利用效率，而糧食水足跡能夠有效地反映水資源的利用效率。

黑龍江省2007—2012年糧食藍水足跡比例變化呈平穩(wěn)趨勢，6年的均值為31.31%，說明黑龍江省糧食生產(chǎn)用水結(jié)構(gòu)6年間基本保持不變。灌溉水量雖然逐年增加，并沒有明顯地影響到糧食作物的生產(chǎn)用水結(jié)構(gòu)。僅2009年由于生育期降水量充沛，糧食藍水足跡比例相比其他年份較低。

2.2.4 糧食作物水足跡與單產(chǎn)需水量比較

圖3 黑龍江省各地區(qū)的單產(chǎn)需水量與糧食水足跡Fig.3 Water requirement per unit yield and grain water footprint in different areas of Heilongjiang Province

2012年黑龍江省各地區(qū)糧食作物的單產(chǎn)需水量見圖3。2012年黑龍江省糧食作物水足跡整體上小于單產(chǎn)需水量，糧食作物生產(chǎn)過程中實際投入的水資源量不足，存在著一定的缺水現(xiàn)象[28]。其中，鶴崗市的糧食水足跡大于單產(chǎn)需水量較多，哈爾濱市、雞西市、綏化市、農(nóng)墾總局糧食水足跡略大于單產(chǎn)需水量；其他9個地區(qū)糧食水足跡均不同程度地小于單產(chǎn)需水量，其中伊春市缺水程度最大，糧食水足跡不到單產(chǎn)需水量的1/2。這一結(jié)果顯示，雖然黑龍江省水資源量相對比較充沛，但是在部分地區(qū)的糧食作物生產(chǎn)過程中仍存在較為嚴重的缺水現(xiàn)象，水資源的高效利用應(yīng)該被充分重視。水資源短缺嚴重的地區(qū)可采取以下措施降低糧食作物生產(chǎn)中水資源的消耗量，提高水資源的利用效率：改善糧食作物種植結(jié)構(gòu)減少作物需水量；提高渠道襯砌率減少灌溉水資源的滲漏損失；修建蓄水建筑物存蓄未被作物直接利用的降水資源。

3 結(jié)論

(1)黑龍江省糧食作物的播種面積逐年增長，糧食生產(chǎn)消耗的灌溉水量也逐年增加，糧食作物用水比例巨大，嚴重制約著黑龍江省除農(nóng)業(yè)以外其他產(chǎn)業(yè)的發(fā)展。

(2)黑龍江省各地區(qū)的糧食水足跡在空間分布上存在顯著差異，大體上呈東北部高，西南部低的趨勢。黑龍江省的糧食藍水足跡明顯低于綠水足跡，綠水資源是糧食作物生產(chǎn)過程中消耗的最主要的水資源類型。2007—2012年黑龍江省各地區(qū)綠水足跡和藍水足跡的空間分布整體上都呈現(xiàn)為相似值之間的空間集聚，說明相鄰地區(qū)間糧食作物生產(chǎn)的社會經(jīng)濟條件具有一定的相似性，值得互相交流借鑒糧食作物種植和節(jié)水經(jīng)驗。

(3)由于2009年黑龍江省糧食作物單位面積產(chǎn)量的顯著降低，使得2007—2012年被分成2007—2008年和2009—2012年2個不同的時期。糧食水足跡在2個時期均存在逐年降低的趨勢，雖然黑龍江省糧食作物的灌溉用水量逐年增加，糧食生產(chǎn)的用水效率卻得到了提升。

(4)黑龍江省農(nóng)墾總局由于糧食生產(chǎn)技術(shù)先進，僅用較少的種植面積生產(chǎn)出較多的糧食作物，但其糧食藍水足跡比例較高，糧食生產(chǎn)比較依賴灌溉水資源，水資源消耗量較大。

(5)黑龍江省糧食作物水足跡整體上小于單產(chǎn)需水量，部分地區(qū)糧食作物生產(chǎn)過程中存在較為嚴重的缺水現(xiàn)象。黑龍江省水資源空間分布不均，各地區(qū)水資源情況差異較大。

1 HOEKSTRA A Y, CHAPAGAIN A K, ALDAYA M M, et al. The water footprint assessment manual: setting the global standard[M].London: Earthscan, 2011.

2 孫世坤，王玉寶，吳普特，等.小麥生產(chǎn)水足跡區(qū)域差異及歸因分析[J].農(nóng)業(yè)工程學報，2015，31(13)：142-148. SUN Shikun, WANG Yubao, WU Pute, et al. Spatial variability and attribution analysis of water footprint of wheat in China[J]. Transactions of the CSAE, 2015,31(13):142-148. (in Chinese)

3 諸大建, 田園宏.虛擬水與水足跡對比研究[J].同濟大學學報：社會科學版，2012，23(4)：43-49. ZHU Dajian, TIAN Yuanhong. Comparative research of virtual water and water footprint[J]. Journal of Tongji University: Social Science Section, 2012,23(4):43-49. (in Chinese)

4 劉靜，吳普特，王玉寶，等.河套灌區(qū)糧食水足跡與虛擬水凈輸出時空演變[J].排灌機械工程學報，2014，32(5)：435-440. LIU Jing, WU Pute, WANG Yubao, et al. Temporal and spatial variation of grain production water footprint and net virtual water export in Hetao Irrigation District[J]. Journal of Drainage and Irrigation Machinery Engineering, 2014,32(5):435-440. (in Chinese)

5 操信春，吳普特，王玉寶，等.中國灌區(qū)糧食生產(chǎn)水足跡及用水評價[J].自然資源學報，2014，29(11)：1826-1835. CAO Xinchun, WU Pute, WANG Yubao, et al. Agriculture water utilization assessment based on water footprint of grain production evaluation in irrigated farmland of China[J]. Journal of Natural Resources, 2014,29(11):1826-1835. (in Chinese)

6 MEKONNEN M M, HOEKSTRA A Y. A global and high-resolution assessment of the green, blue and grey water footprint of wheat[J].Hydrology and Earth System Sciences, 2010,14(7):1259-1276.

7 杜國明，劉彥隨，劉閣.黑龍江省近30年來糧食生產(chǎn)變化及增產(chǎn)因素分析[J].農(nóng)業(yè)現(xiàn)代化研究，2014，35(5)：519-524. DU Guoming, LIU Yansui, LIU Ge. The analysis of changes in grain production and driving factors in Heilongjiang Province recent 30 years[J]. Research of Agricultural Modernization, 2014,35(5):519-524. (in Chinese)

8 張德華，婁廈，劉慧萍.黑龍江省糧食生產(chǎn)地位及潛在問題研究[J].農(nóng)業(yè)現(xiàn)代化研究，2014，35(5)：525-530. ZHANG Dehua, LOU Sha, LIU Huiping. Analysis on the status and potential problems of grain production in Heilongjiang Province[J]. Research of Agricultural Modernization, 2014,35(5):525-530. (in Chinese)

9 盧鐵光，付強.黑龍江省水資源態(tài)勢與糧食生產(chǎn)研究[C]∥中國水利學會2002學術(shù)年會論文集, 2002:202-205.

10 王斌.黑龍江省糧食生產(chǎn)與耗水問題探討[J].節(jié)水灌溉，2015(12)：77-80. WANG Bin. A study on grain production and its water consumption in Heilongjiang Province[J]. Water Saving Irrigation, 2015(12):77-80. (in Chinese)

11 付強，劉巍，劉東,等.黑龍江省灌溉水利用率分形特征與影響因素分析[J/OL].農(nóng)業(yè)機械學報，2016，47(9)： 147-153. http:∥www.j-csam.org/jcsam/ch/reader/view_abstract.aspx?flag=1. DOI:10.6041/j.issn.1000-1298.2016.09.022. FU Qiang, LIU Wei, LIU Dong, et al. Fractal characteristics and influence factors analysis of irrigation water use efficiency in Heilongjiang Province[J/OL]. Transactions of the Chinese Society for Agricultural Machinery, 2016,47(9):147-153. (in Chinese)

12 CHENG K, FU Q, LI T X, et al. Regional food security risk assessment under the coordinated development of water resources[J].Natural Hazards, 2015,78(1):603-619.

13 操信春，吳普特，王玉寶，等.水分生產(chǎn)率指標的時空差異及相關(guān)關(guān)系[J].水科學進展，2014，25(2)：268-274. CAO Xinchun, WU Pute, WANG Yubao, et al. Temporal and spatial variation and correlativity of water productivity indexes in irrigated land of China[J]. Advances in Water Science, 2014,25(2):268-274. (in Chinese)

14 王玉寶，吳普特，孫世坤，等.我國糧食虛擬水流動對水資源和區(qū)域經(jīng)濟的影響[J/OL].農(nóng)業(yè)機械學報，2015，46(10)：208-215. http:∥www.j-csam.org/jcsam/ch/reader/view_abstract.aspx?flag=1. DOI:10.6041/j.issn.1000-1298.2015.10.027. WANG Yubao, WU Pute, SUN Shikun, et al. Impact of virtual water flows of grain on water resources and regional economy in China[J/OL]. Transactions of the Chinese Society for Agricultural Machinery, 2015,46(10):208-215. (in Chinese)

15 黑龍江省水利廳. 黑龍江省水利建設(shè)統(tǒng)計年鑒[M].

16 黑龍江省統(tǒng)計局. 國家統(tǒng)計局黑龍江調(diào)查總隊.黑龍江統(tǒng)計年鑒2015[M]. 北京: 中國統(tǒng)計出版社,2015.

17 劉寶勤.我國糧食虛擬水流動空間格局及其調(diào)控政策[J].水利發(fā)展研究，2010，10(2)：16-20.

18 楊飛，姚作芳，宋佳，等.松嫩平原作物生長季氣候和作物生育期的時空變化特征[J].中國農(nóng)業(yè)氣象，2012，33(1)：18-26. YANG Fei, YAO Zuofang, SONG Jia, et al. Temporal and spatial changes characteristics of the agricultural meteorological factors and crops growth stages in Songnen Plain[J]. Chinese Journal of Agrometeorology, 2012,33(1):18-26. (in Chinese)

19 SUN S, WU P, WANG Y, et al. The impacts of interannual climate variability and agricultural inputs on water footprint of crop production in an irrigation district of China[J].Science of the Total Environment, 2013, 444:498-507.

20 王雪青，陳媛，劉炳勝.中國區(qū)域房地產(chǎn)經(jīng)濟發(fā)展水平空間統(tǒng)計分析——全局Moran’sI、Moran散點圖與LISA集聚圖的組合研究[J].數(shù)理統(tǒng)計與管理，2014，33(1)：59-71. WANG Xueqing, CHEN Yuan, LIU Bingsheng. Exploratory spatial data analysis about the development level of the regional real estate economy in China—the research based on global Moran’sI, Moran scatter plots and LISA cluster map[J]. Journal of Applied Statistics and Management, 2014,33(1):59-71. (in Chinese)

21 付強，劉巍，劉東，等.黑龍江省灌溉用水效率指標體系空間格局研究[J/OL].農(nóng)業(yè)機械學報，2015，46(12)：127-132. http:∥www.j-csam.org/jcsam/ch/reader/view_abstract.aspx?flag=1. DOI:10.6041/j.issn.1000-1298.2015.12.018. FU Qiang, LIU Wei, LIU Dong, et al. Spatial distribution of irrigation water use efficiency index system in Heilongjiang Province[J/OL]. Transactions of the Chinese Society for Agricultural Machinery, 2015,46(12):127-132. (in Chinese)

22 冒小棟，劉瓊芳.江西省人口分布空間自相關(guān)分析[J].華東交通大學學報，2014，31(2)：131-138. MAO Xiaodong, LIU Qiongfang. Spatial autocorrelation analysis of population distribution in Jiangxi Province[J]. Journal of East China Jiaotong University, 2014,31(2):131-138. (in Chinese)

23 關(guān)博，付強，李婷.基于虛擬水理論的三江平原農(nóng)畜產(chǎn)品需水研究[J].農(nóng)業(yè)現(xiàn)代化研究，2008，29(4)：457-459. GUAN Bo, FU Qiang, LI Ting. Research on water requirement of agricultural and livestock products in Sanjiang Plain based on virtual water[J]. Research of Agricultural Modernization, 2008,29(4):457-459. (in Chinese)

24 劉小剛，符娜，李闖，等.河南省主糧作物需水量變化趨勢與成因分析[J/OL].農(nóng)業(yè)機械學報，2015，46(9)：188-197. http:∥www.j-csam.org/jcsam/ch/reader/view_abstract.aspx?flag=1. DOI: 10.6041/j.issn.1000-1298.2015.09.027. LIU Xiaogang, FU Na, LI Chuang, et al. Trends and causes analysis of water requirement for main grain crops in Henan Province[J/OL]. Transactions of the Chinese Society for Agricultural Machinery, 2015,46(9):188-197. (in Chinese)

25 曹連海，吳普特，趙西寧，等.近50年河套灌區(qū)種植系統(tǒng)演化分析[J/OL].農(nóng)業(yè)機械學報，2014，45(7)：144-150. http:∥www.j-csam.org/jcsam/ch/reader/view_abstract.aspx?flag=1. DOI: 10.6041/j.issn.1000-1298.2014.07.023. CAO Lianhai, WU Pute, ZHAO Xining, et al. Plantation system evolution analysis of Hetao Irrigation District in recent 50 years[J/OL]. Transactions of the Chinese Society for Agricultural Machinery, 2014,45(7):144-150. (in Chinese)

26 孫才志, 陳麗新, 劉玉玉.中國農(nóng)作物綠水占用指數(shù)估算及時空差異分析[J].水科學進展，2010，21(5)：637-643. SUN Caizhi, CHEN Lixin, LIU Yuyu. Spatial and temporal variation of crops green water occupancy index in China[J]. Advances in Water Science, 2010,21(5):637-643. (in Chinese)

27 馬德全. 強基礎(chǔ) 調(diào)結(jié)構(gòu) 上科技 提標準 實現(xiàn)墾區(qū)農(nóng)業(yè)連續(xù)大豐收[J].農(nóng)場經(jīng)濟管理，2009(8)：13-14.

28 殷志強，秦小光，李長生.東北三省主要農(nóng)作物耗水量與缺水量研究[J].科技導報，2009，27(13)：42-49. YIN Zhiqiang, QIN Xiaoguang, LI Changsheng. Study on water consumption and deficiency of main crops in northeastern China[J]. Science & Technology Review, 2009,27(13):42-49. (in Chinese)

Analysis of Water Utilization in Grain Production from Water Footprint Perspective in Heilongjiang Province

FU Qiang LIU Ye LI Tianxiao CUI Song LIU Dong CHENG Kun

(SchoolofWaterConservancyandCivilEngineering,NortheastAgriculturalUniversity,Harbin150030,China)

In order to provide decision-making references for the allocation of available water resources in Heilongjiang Province, water footprint theory was used to calculate the grain water footprint during 2007—2012 in different areas of Heilongjiang Province, and the water utilization of food crops was analyzed in Heilongjiang Province during 2007—2012. The results showed that water consumption of grain production was huge, which was increased instantly over the years in Heilongjiang Province. Grain water footprint was less than the water requirement for per unit yield in different areas of Heilongjiang Province, and water shortage was serious in the process of grain crop production in some areas. Through spatial autocorrelation analysis, it showed that the grain green water footprint and the grain blue water footprint of Heilongjiang Province in 2007—2012 both had spatial clustering between similar values. The grain yield of Heilongjiang Farms & Land Reclamation Administration was high, but the ratio of grain blue water footprint was also high, and the grain production was more dependent on irrigation water resources, the water consumption was great. There were significant differences in the spatial distribution of grain water footprint in different areas of Heilongjiang Province, which was generally high in the northeast and low in the southwest. In the six years, the average grain water footprints were 1.15 m3/kg, 1.01 m3/kg, 1.32 m3/kg, 1.04 m3/kg, 0.91 m3/kg and 0.83 m3/kg, respectively. The blue water footprint of grain was significantly lower than that of the green water footprint in Heilongjiang Province, and grain production was highly dependent on the green water resources. In the six years, the average grain blue water footprint ratios were 33.91%, 31.68%, 25.76%, 29.81%, 32.97% and 33.73%, respectively. Due to the significant reduction of grain yield per unit area in 2009, the years of 2007—2012 was divided into two distinct periods of 2007—2008 and 2009—2012. The grain water footprint had the trend of decreasing year by year during the two periods, and the water use efficiency of grain production in the two periods were increased year by year. From the perspective of grain water footprint, the situation of water use in grain production in different areas of Heilongjiang Province was discussed. It can effectively reflect the grain crop production and water resources utilization efficiency, and provide technical support for the water management in grain production in Heilongjiang Province.

grain; water footprint; water resources; spatial distribution; irrigation; Heilongjiang Province

10.6041/j.issn.1000-1298.2017.06.024

2016-10-18

2016-12-03

國家自然科學基金項目(51479032、51279031、51609039)、黑龍江省高校長江學者后備支持計劃項目、黑龍江省水利科技項目(201318、201503)和中國博士后科學基金項目(2016M601410)

付強(1973—)，男，教授，博士生導師，主要從事農(nóng)業(yè)水土資源系統(tǒng)分析等研究，E-mail： fuqiang0629@126.com

S273; TV213.9

A

1000-1298(2017)06-0184-09