王丹丹楊陳玉何鑫王振昌操信春

摘要：基于廣義水資源利用過(guò)程，同時(shí)考慮區(qū)域灌溉效率、水分生產(chǎn)率、水資源來(lái)源及消耗途徑構(gòu)建廣義水資源利用效率綜合評(píng)價(jià)指數(shù)（C指數(shù)）?；谥袊?guó)主要灌區(qū)數(shù)據(jù)，在量化1998年-2010年各省區(qū)C指數(shù)的基礎(chǔ)上探索了其時(shí)空分異特征。結(jié)果顯示：省區(qū)尺度C指數(shù)總體呈增長(zhǎng)態(tài)勢(shì)，全國(guó)值由0550增加到0704，中國(guó)糧食生產(chǎn)用水結(jié)構(gòu)趨于合理；空間自相關(guān)分析結(jié)果顯示，[JP2]代表年省區(qū)尺度C指數(shù)在空間上呈現(xiàn)顯著的聚集現(xiàn)象，全局莫蘭I的檢驗(yàn)值均超過(guò)了001的置信水平；C指數(shù)較大省區(qū)集中于黃淮海平原和長(zhǎng)江流域北部，低值省區(qū)則廣泛分布于長(zhǎng)江以南、東北及西北地區(qū)；各省區(qū)局部自相關(guān)屬性穩(wěn)定，HH和LL的省區(qū)在20個(gè)左右。從自然氣候特點(diǎn)、農(nóng)業(yè)生產(chǎn)條件等方面分析了C指數(shù)的時(shí)空分異特征。研究成果為農(nóng)業(yè)水資源利用評(píng)價(jià)方法和水資源管理宏觀政策制定提供參考。

關(guān)鍵詞：藍(lán)綠水；農(nóng)業(yè)用水綜合指數(shù)；空間格局；灌區(qū)；糧食安全；中國(guó)

中圖分類(lèi)號(hào)：TV213文獻(xiàn)標(biāo)志碼：A文章編號(hào)：

16721683（2016）06005606

Spatial and temporal differences of agricultural water utilization composite index in irrigated farmland of China

WANG Dandan1，YANG Chenyu2，HE Xin2，WANG Zhenchang2，CAO Xinchun2

（1.School of Resource and Environmental Sciences，Chifeng University，Chifeng 024000，China；2.Key Laboratory of HighEffective Irrigation and Drainage and Agricultural Water and Soil Environment in Southern China，Ministry of Education，Hohai University，Nanjing 210098，China）

Abstract：

Based on the use of generalized water resources and the existing evaluation index system of agricultural water，regional irrigation efficiency，water productive efficiency，source of water resources and consumptive ways to build agricultural water utilization composite index（Cindex）were considered.Based on data collected from major irrigation districts of China，this paper quantified irrigation Cindex from 1998 to 2010，explored the spatial differentiation characteristics.The results showed that：the Cindex at provincial level presented an increasing trend，and the national Cindex increased from 0.550 to 0.704，the structure of water for food production trended to be rationalized；Spatial autocorrelation analysis results showed that，C index on provincial scale gathered significantly in space，Moran′s I inspection values were higher than the confidence level（0.01）.High value provinces were concentrated in HuangHuaiHai Plain and north of Changjiang River，low value provinces were widely distributed in the south of Changjiang River，Northeast and Southeast China；Attribute of local autocorrelation was stable，and about 20 provinces showed highhigh positive autocorrelation（HH）and lowlow positive autocorrelation（LL）.Spatial differentiation characteristics of Cindex from natural features，agricultural production conditions were also analyzed.This paper provides references for development of evaluation methods for agricultural water use and macro policies of water resources management.

Key words：bluegreen water resources；agricultural water utilization composite index（Cindex）；spatial distribution pattern；irrigation district；food security；China

1研究背景

農(nóng)業(yè)水資源利用效率區(qū)域差異解析是進(jìn)行水資源管理的重要依據(jù)。傳統(tǒng)研究主要針對(duì)農(nóng)業(yè)生產(chǎn)過(guò)程，通過(guò)建立灌溉效率、水分生產(chǎn)率及綜合生產(chǎn)效率指標(biāo)來(lái)分析農(nóng)業(yè)水資源利用問(wèn)題[12]?，F(xiàn)代水資源管理研究通過(guò)水足跡來(lái)衡量人類(lèi)生產(chǎn)或消費(fèi)活動(dòng)對(duì)水資源的影響[36]。兩者的共同基礎(chǔ)為廣義水資源利用和水資源量平衡，傳統(tǒng)的農(nóng)業(yè)用水效用評(píng)價(jià)中注重區(qū)域水文循環(huán)下的水量平衡，而水足跡理論更強(qiáng)調(diào)在消費(fèi)社會(huì)中區(qū)域間的水資源量產(chǎn)銷(xiāo)平衡。水足跡理論強(qiáng)調(diào)區(qū)分藍(lán)、綠水資源和虛擬水的重要作用，在更大尺度上評(píng)價(jià)水資源的利用效益，使得水資源管理成為兼顧區(qū)域內(nèi)外部相互影響的系統(tǒng)工程。建立水資源利用與農(nóng)業(yè)產(chǎn)出之間關(guān)系的評(píng)價(jià)指標(biāo)是兩類(lèi)方法契合點(diǎn)。學(xué)者在基于傳統(tǒng)研究指標(biāo)的特定區(qū)域農(nóng)業(yè)用水評(píng)價(jià)上做了大量工作，主要包括作物田間需水、耗水規(guī)律[79]，不同角度灌溉效率及水分生產(chǎn)率評(píng)價(jià)指標(biāo)構(gòu)建、量化及其尺度效應(yīng)[1012]，農(nóng)業(yè)水資源高效優(yōu)化配置及農(nóng)業(yè)節(jié)水潛力計(jì)算方法的建立等[1315]。在水足跡與農(nóng)業(yè)用水評(píng)價(jià)結(jié)合方面，分析區(qū)域間農(nóng)業(yè)用水狀況差異，并進(jìn)一步評(píng)價(jià)消費(fèi)水足跡和虛擬水流動(dòng)現(xiàn)象帶來(lái)的水資源節(jié)約效益的相關(guān)研究也取得一定進(jìn)展[1618]。然而，當(dāng)前研究?jī)H以田間尺度水分生產(chǎn)率作為比較區(qū)域間農(nóng)業(yè)用水效率的指標(biāo)，這使得研究結(jié)論具有局限性。因?yàn)橛绊戅r(nóng)業(yè)用水狀況的因素眾多，包括氣候條件、作物品種、灌溉設(shè)施、技術(shù)推廣等，某一類(lèi)指標(biāo)很難全面地體現(xiàn)[1920]。同時(shí)，從水資源投入糧食產(chǎn)出關(guān)系衡量的角度講，田間尺度水分生產(chǎn)率也并非最具代表性的指標(biāo)[21]。本文參考傳統(tǒng)農(nóng)業(yè)用水效用評(píng)價(jià)和水足跡理論中的相關(guān)范式，選擇從不同角度衡量農(nóng)業(yè)用水狀況指標(biāo)來(lái)建立廣義水資源利用效率綜合評(píng)價(jià)指數(shù)（C指數(shù)），量化中國(guó)灌區(qū)C指數(shù)并利用空間統(tǒng)計(jì)學(xué)方法分析其時(shí)空分異特征。為相關(guān)研究和宏觀政策制定提供參考。[HJ]

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

[JP+1]空間自相關(guān)分析方法反映一個(gè)地理屬性與鄰近區(qū)域同一屬性相關(guān)程度，可分為全局和局部空間自相關(guān)。全局空間自相關(guān)對(duì)地理屬性在整個(gè)區(qū)域的空間特征進(jìn)行描述，判斷該屬性值在總體空間上是否存在聚集現(xiàn)象[23]。在此，通過(guò)對(duì)全局Moran′s I指標(biāo)的計(jì)算來(lái)分析C指數(shù)的總體空間關(guān)聯(lián)程度。通過(guò)局部Moran′s I測(cè)算局部地理單元C指數(shù)的空間異質(zhì)性，并推算出聚集地的地理位置、范圍及其隨時(shí)間變化情況，并以LISA聚集圖的形式表征。與總體的平均水平相比，可以判斷出一個(gè)省區(qū)C指數(shù)在總體中所處的位置。LISA聚集圖中屬性分為四種類(lèi)型[7]：HH，表示某該省區(qū)和周?chē)^(qū)的地理屬性值都較高；LH，表示該省區(qū)較低而周?chē)^(qū)較高；HL，表示該省區(qū)較高而周?chē)^低；LL，表示該省區(qū)和周?chē)^(qū)都較低。在給定顯著性水平時(shí)，可以判斷各區(qū)域局部空間自相關(guān)屬性的顯著性。

2.3數(shù)據(jù)來(lái)源

降水量數(shù)據(jù)由中國(guó)氣象科學(xué)數(shù)據(jù)共享服務(wù)網(wǎng)下載；各灌區(qū)代表年的灌溉效率、糧食產(chǎn)量、灌溉面積、復(fù)種指數(shù)、灌溉用水量等數(shù)據(jù)來(lái)源于中國(guó)大型灌區(qū)網(wǎng)站（http：//www.dxgq.org.cn）及部分的灌區(qū)管理局；還參考了研究時(shí)段內(nèi)的中國(guó)統(tǒng)計(jì)年鑒、中國(guó)水利年鑒、中國(guó)水資源公報(bào)、各省區(qū)統(tǒng)計(jì)年鑒及水資源公報(bào)等資料；因數(shù)據(jù)所限，選取1998年-2010年的起始年、中間年（2005年）及末尾年為代表年，來(lái)分析C指數(shù)的時(shí)空分布特征。

3結(jié)果與分析

[BT（3]3.1[ZK（][JP+1]不同年份廣義水資源利用效率綜合評(píng)價(jià)指數(shù)（C指數(shù)）[ZK）][BT）]

將各省區(qū)C指數(shù)進(jìn)行對(duì)有效灌溉面積的加權(quán)得到全國(guó)灌區(qū)不同年份C指數(shù)值，見(jiàn)圖1。為便于分析和對(duì)比，圖1還給出了全國(guó)灌區(qū)灌溉效率（η）、廣義水資源水分生產(chǎn)率（Pwu）、綠水比例（Rg）及有效耗水比例（RET）四個(gè)指標(biāo)在不同年份的表現(xiàn)。

由圖1可知，中國(guó)灌區(qū)η、Pwu、Rg及RET均在0300～1000之間，且不同年份四個(gè)指標(biāo)在數(shù)值上大小關(guān)系為：Pwu>RET>η>Rg。η、Pwu、Rg、RET和C指數(shù)代表年均值分別為0452、0860、0365、0652以及0619，農(nóng)業(yè)生產(chǎn)用水綜合指數(shù)與有效耗水比例在數(shù)值最為接近。C指數(shù)在三個(gè)代表年的值分別為0550、0603及0704，呈穩(wěn)步增長(zhǎng)態(tài)勢(shì)，說(shuō)明綜合灌溉效率、水分生產(chǎn)率、水資源來(lái)源構(gòu)成以及消耗途徑看，中國(guó)灌區(qū)糧食生產(chǎn)用水向更加合理的趨勢(shì)發(fā)展。η、Pwu、Rg及RET均呈隨時(shí)間增大趨勢(shì)，其中Pwu最為明顯，這也直接導(dǎo)致了C指數(shù)的隨時(shí)間增加。四個(gè)指標(biāo)及C指數(shù)的增長(zhǎng)主要由糧食單產(chǎn)水平的提高和灌溉設(shè)施改進(jìn)使得灌溉引水量減少引起。如1998年中國(guó)灌區(qū)的平均灌水量約為8 900 m3/hm2，對(duì)應(yīng)耕地上的糧食單產(chǎn)約為9 900 kg/hm2，2010年這兩個(gè)參數(shù)分別變化為7 300 m3/hm2和12 200 kg/hm2，在降水量變化不明顯的情況下，η、Pwu、Rg及RET分別由0420、0725、0350和0623增長(zhǎng)到0492、1020、0390以及0690。為了C指數(shù)隨時(shí)間變化的區(qū)域差異，計(jì)算出各省區(qū)C指數(shù)的年均變化幅度，列于表2。

結(jié)果顯示，各省區(qū)C指數(shù)均隨時(shí)間增加，增加速度存在一定差異。北京和內(nèi)蒙古的年均增幅在300%以上，為中國(guó)糧食生產(chǎn)用水合理性變強(qiáng)最快的區(qū)域；天津、吉林、湖北和重慶的年均增加幅度均超過(guò)250%，屬于灌區(qū)農(nóng)業(yè)生產(chǎn)用水綜合指數(shù)增長(zhǎng)較快的省區(qū)；還有13個(gè)省區(qū)在全國(guó)平均水平（200%）之上。12個(gè)省區(qū)的增長(zhǎng)速率在全國(guó)值以下，其中安徽最小，為唯一不足100%的省份。總體來(lái)看，北方地區(qū)的C指數(shù)增長(zhǎng)較快而南方省區(qū)相對(duì)較慢，如增幅最大前4個(gè)省區(qū)均在北方而最小的4省區(qū)均位于長(zhǎng)江流域。C指數(shù)的全面增加說(shuō)明中國(guó)農(nóng)業(yè)水資源利用效率的顯著提升，對(duì)于緩解區(qū)域水資源壓力、提升水資源利用效益以及促進(jìn)區(qū)域經(jīng)濟(jì)發(fā)展均具有重要意義。

3.2農(nóng)業(yè)生產(chǎn)用水綜合指數(shù)空間差異

為分析C指數(shù)的空間分散程度及其隨時(shí)間變化情況，計(jì)算出不同年份省級(jí)行政區(qū)尺度C指數(shù)的最大值、最小值、極差及標(biāo)準(zhǔn)差等統(tǒng)計(jì)量，列于表3。同時(shí)，由于C指數(shù)值在各代表年省區(qū)間大小關(guān)系的一致性較好，用各省區(qū)C指數(shù)三個(gè)年份的平均值來(lái)分析空間格局，見(jiàn)圖2。

由表3可知，省級(jí)行政區(qū)尺度C指數(shù)的最大值、最小值及平均值（算術(shù)平均值）均隨時(shí)間穩(wěn)步增長(zhǎng)，這與圖2、表2中糧食生產(chǎn)用水合理性全面升高的研究結(jié)果相吻合。與C指數(shù)全國(guó)加權(quán)值值相比，不同年份的省區(qū)算術(shù)平均值均偏小，說(shuō)明總體上農(nóng)業(yè)生產(chǎn)規(guī)模較大省區(qū)的糧食生產(chǎn)用水合理性更占優(yōu)勢(shì)。省區(qū)C指數(shù)的極差和標(biāo)準(zhǔn)差也均在1998年-2005年-2010年區(qū)間內(nèi)逐漸增大，說(shuō)明中國(guó)糧食生產(chǎn)用水合理性的區(qū)域差距正在不斷擴(kuò)大。

中國(guó)灌區(qū)農(nóng)業(yè)生產(chǎn)用水綜合指數(shù)的省區(qū)差異較大，總體空間分布格局為黃河中下游及華北地區(qū)的省區(qū)較大而東北、華南、東南及西北部分省區(qū)較小。圖2顯示，陜西的C指數(shù)最高，達(dá)0924，陜西也是唯一超過(guò)0900的省份；山東、河南、河北及北京分別為0893、0867、0859及0859，分列2位至5位；其他省區(qū)均在0800以下，其中7個(gè)省區(qū)大于全國(guó)值的0619；19個(gè)小于全國(guó)值的省區(qū)中的10個(gè)不足0500，包括東南沿海的廣東、廣西、福建、海南，西北的新疆、青海和東北的遼寧、吉林、黑龍江；黑龍江的C指數(shù)全國(guó)最低，為0389，不足陜西、山東、河南、河北及北京五省市的一半；四川、湖北和江西的C指數(shù)值分別為0761、0695及0647均為南方為數(shù)不多的高值省區(qū)之一；除四川外的西南諸省區(qū)C指數(shù)均接近全國(guó)值，處于中等水平。

氣候狀況、作物播種類(lèi)型、農(nóng)業(yè)生產(chǎn)條件、經(jīng)濟(jì)發(fā)展水平等影響區(qū)域灌溉效率、水分生產(chǎn)率、水資源的因素共同決定了省區(qū)度C指數(shù)的大?。宏P(guān)中平原和黃淮海地區(qū)地理?xiàng)l件優(yōu)越、土地肥沃，糧食生產(chǎn)效率高的同時(shí)面臨比較嚴(yán)重的缺水問(wèn)題，單位面積上產(chǎn)出的糧食可達(dá)130 ton/hm2，灌溉用水量在4 000 m3/hm2上下，使得水分生產(chǎn)率（Pwu）在1200 m3/kg左右，雖然灌溉效率（η）與其他省份差異不大，但綠水比例（Rg）和有效耗水比例（RET）均分別在40%和70%左右，明顯高于其他省區(qū)和全國(guó)水平，故C指數(shù)在省區(qū)中占明顯優(yōu)勢(shì)；經(jīng)濟(jì)發(fā)展較快的華南和東南地區(qū)以種植水田作物為主，降水量和灌溉用水同時(shí)較多且糧食單產(chǎn)也不占優(yōu)勢(shì)，這些地區(qū)的Pwu僅在0500 m3/kg左右，Rg和RET分別為低于全國(guó)值的30%和60%，所以農(nóng)業(yè)生產(chǎn)綜合用水C指數(shù)較低。東北地區(qū)也種植了大量水稻，灌溉引用水量較大，雖然糧食單產(chǎn)水平較高，但因復(fù)種程度、單位耕地糧食產(chǎn)出受到限制，加上η較低，糧食生產(chǎn)用水綜合效率還有很大的提升空間。西北的青海、新疆及寧夏等地區(qū)的C指數(shù)較小主要因降水資源有限、灌水量大從降低了Rg和RET導(dǎo)致；西南地區(qū)將水資源豐富、灌溉水資源投入較少，使得Rg與RET較大，雖然糧食單產(chǎn)和Pwu較低，C指數(shù)仍能處于全國(guó)中上游位置。

[BT（3]3.3[ZK（][JP+1]農(nóng)業(yè)生產(chǎn)用水綜合指[JP+2]數(shù)的空間自相關(guān)分析[ZK）][BT）]

表4列出了不同時(shí)間C指數(shù)全局空間自相關(guān)的計(jì)算結(jié)果。不同年份全局Moran′s I估計(jì)值均明顯大于零，同時(shí)對(duì)應(yīng)的ZScore均大于標(biāo)準(zhǔn)正態(tài)分布函數(shù)概率為001的臨界值，即通過(guò)了顯著性檢驗(yàn)。代表年C指數(shù)均值的全局Moran′s I、ZScore及Pvalue分別為0264 8、3818 3以及0000 1。說(shuō)明灌區(qū)C指數(shù)在省級(jí)行政區(qū)間的存在正自相關(guān)關(guān)系，總體上表現(xiàn)為相似大小的值在地理上顯著地集中。灌區(qū)C指數(shù)與區(qū)域的氣候條件、作物生產(chǎn)類(lèi)型、水土資源等自然因素關(guān)系密切，同時(shí)也受灌溉水平、糧食單產(chǎn)、生產(chǎn)資料投入及管理方式等經(jīng)濟(jì)社會(huì)條件的影響，而以上各因素大體上均能表現(xiàn)出地理上的類(lèi)似性，這是大小相似C指數(shù)省區(qū)在地理空間上呈現(xiàn)聚集現(xiàn)象的重要原因。

大部分省區(qū)的局部空間自相關(guān)屬性在不同時(shí)間無(wú)明顯變化；三個(gè)年份表現(xiàn)為HL的省區(qū)均不多于4個(gè)，其余均為HH、LL及LH省區(qū)；各年份表現(xiàn)為HH、LL的省區(qū)個(gè)數(shù)之和分別為21、19及20，占省區(qū)的絕大多數(shù)；HH省區(qū)以華北平原為中心連片分布并涉及到部分的西南省區(qū)，LL省區(qū)則主要分布于長(zhǎng)江以南和東北地區(qū)，均表現(xiàn)為地理上的正相關(guān)；表現(xiàn)為L(zhǎng)H及HL省區(qū)在空間上分布相對(duì)散亂；表5顯示，HH（11個(gè)）、LL（9個(gè)）省區(qū)個(gè)數(shù)在各年份維持不變，只是顯著區(qū)減少了一個(gè)；LH省區(qū)的數(shù)量由7增至8，HL省區(qū)相應(yīng)減少1個(gè)?？傮w說(shuō)來(lái)，13年來(lái)C指數(shù)的局部分布特征沒(méi)有發(fā)生明顯的變化，這與省區(qū)間自然條件和社會(huì)經(jīng)濟(jì)狀況差異的穩(wěn)定有關(guān)。具體來(lái)說(shuō)，天津一直處于LH區(qū)，其C指數(shù)明顯低于與之相鄰的河北與北京，不僅因?yàn)榫o鄰HH省區(qū)，而且因?yàn)樵谧匀粭l件、灌溉設(shè)施、種植結(jié)構(gòu)和糧食產(chǎn)量相當(dāng)相似的情況下，天津灌水量較大幅度大于與之相鄰的北京、河北，使得Pwu、Rg及RET均小于兩省市；內(nèi)蒙古與江蘇的C指數(shù)處全國(guó)中下游水平，同時(shí)因受HH區(qū)輻射作用明顯，故均長(zhǎng)期位于LH區(qū)；江西表現(xiàn)為HL，其代表年均C指數(shù)為0647，明顯大于鄰省浙江、福建、廣東的0450左右，江西與鄰省相比經(jīng)濟(jì)發(fā)展滯后而有較好的農(nóng)業(yè)基礎(chǔ)，局部自相關(guān)屬性受與其社會(huì)經(jīng)濟(jì)條件類(lèi)似的湖南、湖北影響較大；此外，黑龍江與遼寧分別由HL和LL變?yōu)長(zhǎng)L與LH，而在云南HL與HH之間搖擺。不同時(shí)期C指數(shù)的全局及局部空間自相關(guān)特性差異較小，說(shuō)明中國(guó)糧食生產(chǎn)用水狀況的空間格局將穩(wěn)定存在。

由于C指數(shù)能夠更加全面衡量區(qū)域間農(nóng)業(yè)用水狀況差異，各地區(qū)應(yīng)結(jié)合自身糧食生產(chǎn)能力與水資源狀況等條件，針對(duì)各農(nóng)業(yè)用水評(píng)價(jià)指標(biāo)（η、Pwu、Rg及RET）的具體狀況制定相應(yīng)對(duì)策以提升C指數(shù)，從而緩解區(qū)域水資源壓力、提高農(nóng)業(yè)水資源利用效率和促進(jìn)經(jīng)濟(jì)社會(huì)持續(xù)發(fā)展。例如南方地區(qū)降水資源豐富、糧食生產(chǎn)任務(wù)小、經(jīng)濟(jì)較為發(fā)達(dá)的同時(shí)其他行業(yè)對(duì)水資源需求較大，應(yīng)在充分利用降水資源的基礎(chǔ)上提高η與Rg，將節(jié)約的灌溉水用于工業(yè)生產(chǎn)與服務(wù)業(yè)，提升水資源經(jīng)濟(jì)效益；東北糧食主產(chǎn)區(qū)的Pwu較高，而面臨灌溉水資源管理粗放的同時(shí)降水有效利用率不高，所以提高η與RET是提高C指數(shù)的主要途徑；而華北和西北省區(qū)面臨嚴(yán)峻的水資源壓力，農(nóng)業(yè)水資源供應(yīng)有限，不僅要重視農(nóng)業(yè)水分生產(chǎn)率Pwu，還應(yīng)該重視有效耗水比例RET的提高以減少水資源的浪費(fèi)。

4結(jié)論

由于灌溉效率、廣義水資源水分生產(chǎn)率、綠水比例及有效耗水比例均不同程度地增加，中國(guó)及各省區(qū)灌區(qū)糧食生產(chǎn)用水C指數(shù)呈全面上升態(tài)勢(shì)，全國(guó)值由1998年的0550增長(zhǎng)到2010年的0704，糧食生產(chǎn)用水合理性在增強(qiáng)。省級(jí)行政區(qū)尺度C指數(shù)的差距較大，多年值在黑龍江的0398至陜西的0924之間變化；C指數(shù)的省區(qū)間有隨時(shí)間增大趨勢(shì)，中國(guó)區(qū)域間農(nóng)業(yè)用水狀況的改善步調(diào)一致性不強(qiáng)?？臻g自相關(guān)方法分析發(fā)現(xiàn)，其在空間上呈顯著聚集現(xiàn)象，高值省分緊密分布在黃淮海平原及其周邊區(qū)域，低值省區(qū)則主要集中于長(zhǎng)江以南和東北地區(qū)。局部自相關(guān)空間屬性為HH、LL、LH和HL的省區(qū)個(gè)數(shù)均相對(duì)穩(wěn)定，顯著省區(qū)個(gè)數(shù)在減少。全局和局部自相關(guān)分析結(jié)果隨時(shí)間變化不大，中國(guó)糧食生產(chǎn)用水的時(shí)空分異格局穩(wěn)定。

本文從區(qū)域用水狀況角度進(jìn)行了評(píng)估并分析其時(shí)空分異格局，可為區(qū)域間農(nóng)業(yè)水資源利用評(píng)價(jià)提供新思路，也為全國(guó)及區(qū)域農(nóng)業(yè)生產(chǎn)及水資源宏觀管理提供理論參考。由于水資源問(wèn)題涉及社會(huì)經(jīng)濟(jì)系統(tǒng)的各個(gè)方面，提升農(nóng)業(yè)水資源綜合利用效率是任何區(qū)域均需重視的問(wèn)題。然而，不同地區(qū)的種植結(jié)構(gòu)、氣候條件、水資源條件存在差異，導(dǎo)致對(duì)不同農(nóng)業(yè)用水評(píng)價(jià)指標(biāo)的重視程度不盡相同，針對(duì)農(nóng)業(yè)用水過(guò)程本身的水資源利用效率評(píng)價(jià)指數(shù)提升的急迫性和策略制定還需與區(qū)域自然及農(nóng)業(yè)生產(chǎn)條件相結(jié)合，以為區(qū)域農(nóng)業(yè)用水效率提升明確方向。而基于現(xiàn)代水資源管理理念的農(nóng)業(yè)水資源具體管理政策的建立與實(shí)施，需要在水資源利用相對(duì)狀況評(píng)價(jià)的基礎(chǔ)上，進(jìn)行以提升全國(guó)水資源利用效益和保障全國(guó)范圍內(nèi)糧食安全生產(chǎn)的水土資源優(yōu)化配置研究是需要考慮的問(wèn)題。

參考文獻(xiàn)（References）：

[1]Molden D，Sakthivadivel R，Christopher JP，et al.Indicators for Comparing Performance of Irrigated Agricultural Systems.IWMI Research Report no.20[R].International Water Management Institute（IWMI）：Colombo，Sri Lanka.DOI：10.1002/ird.1816

[2][JP2]Heydari N.Water productivity in agriculture：challenges in concept，terms and values[J].Irrig.and Drain.，2014，63：2228.

[3]Allan JA.Virtual water eliminates water wars？ A case study from the Middle East.In Hoekstra AY edited.Virtual Water Trade：Proceedings of the International Expert Meeting on Virtual Water Trade.Value of Water Research Report Series No12[R].IHE Delft，the Netherlands，2003：137145.

[4]Cao X C，Wu P T，Wang Y B，et al.Water footprint of grain product in irrigated farmland of China[J].Water Resources Management，2014，28（8）：22132227.DOI：10.1007/s1126901406071

[5]孫克，徐中民.基于地理加權(quán)回歸的中國(guó)灰水足跡人文驅(qū)動(dòng)因素分析[J].地理研究，2016，35（1）：3748.（SUN Ke，XU Zhongmin.The impacts of human driving factors on grey water footprint in China using a GWR model[J].Geographical Research，2016，35（1）：3748.（in Chinese）） DOI：10.11821/dlyj201601004

[6]Aldaya M M，Allan J A，Hoekstra A Y.Strategic importance of green water in international crop trade[J].Ecological Economics，2010，69（4）：887894.DOI：10.1016/j.ecolecon.2009.11.001

[7]操信春，吳普特，王玉寶，等.中國(guó)灌區(qū)水分生產(chǎn)率及其時(shí)空差異分析[J].農(nóng)業(yè)工程學(xué)報(bào)，2012，28（13）：17.（CAO Xinchun，WU Pute，WANG Yubao，et al.Analysis on temporal and spatial differences of water productivity in irrigation districts in China[J].Transactions of the Chinese Society of Agricultural Engineering，2012，28（13）：17.（in Chinese））DOI：103969/j.issn.10026819.2012.13.001

[8]張芊，任理.應(yīng)用根系層水質(zhì)模型分析冬小麥夏玉米輪作體系的農(nóng)田水氮利用效率Ⅱ：模型的驗(yàn)證與情景分析[J].水利學(xué)報(bào)，2012，43（3）：354362.（ZHANG Qian，REN Li.Applying application of root zone water quality model to simulate water and nitrogen use efficiency of winter wheatsummer maize double cropping system II.Model validation and scenario analysis[J].Journalof Hydraulic Engineering，2012，43（3）：354362.（in Chinese）） DOI：05599350（2012）03035409

[9]Cai X，Yang Y，Ringler C，et al.Agricultural water productivity assessment for the Yellow River Basin[J].Agricultural Water Management，2011，98（8）：12971306.[ZK）]

[10][ZK（#]陳皓銳，黃介生，伍靖?jìng)?，?井渠結(jié)合灌區(qū)用水效率指標(biāo)尺度效應(yīng)研究框架[J].農(nóng)業(yè)工程學(xué)報(bào)，2009，25（8）：17.（CHEN Haorui，HUANG Jiesheng，WU Jingwei，et al.Scale effect research framework of irrigation water use efficiency indices in wellcanal combined irrigation area[J].Transactions of the Chinese Society of Agricultural Engineering，2009，25（8）：17.（in Chinese））DOI：103969/J.issn.10026819.2009.08.001

[11]崔遠(yuǎn)來(lái)，董斌，李遠(yuǎn)華.水分生產(chǎn)率指標(biāo)隨空間尺度變化規(guī)律[J].[JP+2]水利學(xué)報(bào)，2006，37（1）：4551.（CUI Yuanlai，DONG Bin，LI Yuanhua.Variation of water productivity in different spatial scales[J].Journal of Hydraulic Engineering，2006，37（1）：4551.（in Chinese））DOI：0559-9350（2006） 01004507

[12]Zwart S，Bastiaanssen W，F(xiàn)raiture C，et al.A global benchmark map of water productivity for rainfed and irrigated wheat[J].Agricultural Water Management，2010，97（10）：16171627.DOI：10.1016/j.agwat.2010.05.018

[13]Juliet C S，Heather C，Peter H G.Potential water savings associated with agricultural water efficiency improvements：a case study of California，USA[J].Water Policy，2012，14（2）：194213.DOI：10.2166/wp.2011.017

[14][JP2]Keramatzadeh A，Chizari A H，Moore R.Economic optimal allocation of agriculture water：mathematical programming approach[J].Journal of Agricultural Science and Technology，2011，13（4）：477490.

[15]Chapagain A K，Hoekstra A Y，savenije H H G.Water saving through international trade of agricultural products[J].Hydrol.Earth Syst.Sci.，2006，10（3）：455468.

[16][ZK（#][JP3]Fader M，Gerten D，Thammer M，et al.Internal and external greenblue agricultural water footprints of nations，and related water and land savings through trade[J].Hydrol.Earth Syst.Sci.，2011，15（5）：16411660.DOI：10.5194/hess1516412011.

[17]Sun S K，Wu P T，Wang Y B，et al.The virtual water content of major grain crops and virtual water flows between regions in China[J].Journal of the Science of Food and Agriculture，2013，93（6）：14271437.DOI：10.1002/jsfa.5911.

[18][JP3]馬濤，陳家寬.虛擬水貿(mào)易在解決中國(guó)和全球水危機(jī)中的作用[J].生態(tài)經(jīng)濟(jì)，2006（11）：2226.（MA Tao，CHEN Jiakuan.The role of virtual water trade in efforts to resolve china and global water crisis[J].Ecological Economy，2006（11）：2226.（in Chinese））

[19]許迪，龔時(shí)宏，李益農(nóng)，等.作物水分生產(chǎn)率改善途徑與方法研究綜述[J].水利學(xué)報(bào)，2010，41（6）：631639.（XU Di，GONG Shihong，LI Yinong，et al.Overview of recent study on improvement approaches and methods for crop water productivity[J].Journal of Hydraulic Engineering，2010，41（6）：631639.（in Chinese））DOI：05599350（2010）06063109

[20]Jensen ME.Beyond irrigation efficiency[J].Irrig.Sci.2007，25（3），233245.DOI：10.1007/s0027100700605

[21]操信春，吳普特，王玉寶，等.水分生產(chǎn)率指標(biāo)的時(shí)空差異及相關(guān)關(guān)系[J].水科學(xué)進(jìn)展，2014，25（2）：268274.（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）：268274.（in Chinese）） DOI：10.14042/j.cnki.32.1309.2014.02.011.

[22]Sun S K，Wu P T，Wang Y B，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：498507.DOI：10.1016/j.scitotenv.2012.12.016.

[23]孟斌，王勁峰，張文忠，等.基于空間分析方法的中國(guó)區(qū)域差異研究[J].地理科學(xué)，2005，25（4）：393400.（MENG Bin，WANG JinFeng，ZHANG WenZhong，et al.Evaluation of Regional Disparity in China Based on Spatial Analysis[J].Scientia Geographical Sinica，2005，25（4）：393400.（in Chinese））