□方 蘭 池營(yíng)營(yíng)

一、引言

(一)研究背景

自上世紀(jì)90年代起，山東省大力推廣蔬菜種植，經(jīng)過(guò)十幾年發(fā)展探索，蔬菜產(chǎn)業(yè)已成為山東省經(jīng)濟(jì)發(fā)展的支柱產(chǎn)業(yè)，成為中國(guó)第一蔬菜大省。蔬菜產(chǎn)業(yè)的高收益受到農(nóng)戶青睞，并已成為農(nóng)戶增收的主要源泉。但是蔬菜種植過(guò)程中的資本、勞動(dòng)密集型投入是實(shí)現(xiàn)高產(chǎn)出的主要原因，而且蔬菜種植過(guò)程中需要反復(fù)灌溉，耗水量較大。雖然山東省位于中國(guó)東部沿海地區(qū)，地理位置優(yōu)越、資源豐富，但同時(shí)又是中國(guó)嚴(yán)重缺水的省份之一。山東省人均水資源占有量?jī)H334m3，不到全國(guó)平均水平的六分之一，居全國(guó)倒數(shù)第4位。同時(shí)，黃河上中游各省農(nóng)業(yè)、工業(yè)耗水量迅速增加，導(dǎo)致下游水量急劇減少，黃河斷流次數(shù)、持續(xù)時(shí)間越來(lái)越長(zhǎng)，這也是山東省水資源短缺的主要原因。先進(jìn)的節(jié)水灌溉技術(shù)在山東省已有了很大程度的應(yīng)用，但并未得到廣泛普及，農(nóng)村地區(qū)土渠占很大比例，灌溉過(guò)程中滲透嚴(yán)重，浪費(fèi)了大量水資源。不僅如此，山東省蔬菜大部分用于外銷，無(wú)形中又導(dǎo)致了虛擬水的流失，這加劇了山東省水資源危機(jī)。面對(duì)蔬菜種植過(guò)程中水資源的高消耗狀況，如何在有限的水資源制約下，提高蔬菜產(chǎn)業(yè)水資源利用效率，實(shí)現(xiàn)水資源短缺地區(qū)農(nóng)業(yè)的高效用水是本文研究的主要目標(biāo)。

(二)文獻(xiàn)綜述

國(guó)外學(xué)者對(duì)灌溉水利用效率的研究時(shí)間較早，理論及技術(shù)應(yīng)用都比較成熟。Karagiannis，G.，Tzouvelekas，V［2］等在 1998 － 1999 年希臘克里特地區(qū)50個(gè)農(nóng)戶數(shù)據(jù)的基礎(chǔ)上，運(yùn)用SFA超對(duì)數(shù)模型衡量灌溉用水效率，并分析了影響灌溉效率的影響因素。Dhehibi，Lachaal［3］等運(yùn)用BC—1995隨機(jī)前沿分析模型，以尼泊爾和突尼斯144個(gè)柑橘農(nóng)戶數(shù)據(jù)為基礎(chǔ)，測(cè)算農(nóng)戶的技術(shù)和灌溉用水效率。Speelman et al［4］基于2005年南非西北省Zeerust市60個(gè)調(diào)查農(nóng)戶的截面數(shù)據(jù)，采用數(shù)據(jù)包絡(luò)分析方法分別對(duì)規(guī)模收益不變和可變情況下的灌溉用水效率進(jìn)行估計(jì)，結(jié)果表明耕地面積、土地所有權(quán)、灌溉方式、灌溉項(xiàng)目類型及種植結(jié)構(gòu)將會(huì)顯著影響灌溉用水效率。Abeer MunlaHasan［5］詳細(xì)研究了敘利亞的灌溉用水效率，敘利亞灌溉用水效率仍舊很低，政府部門需要增加資金投入、推廣先進(jìn)灌溉技術(shù)，改善灌溉用水效率。

在農(nóng)作物和灌溉用水效率方面，王會(huì)肖、劉昌明［6］詳細(xì)分析了水分利用效率的內(nèi)涵，分別就葉片水平上的水分利用效率、群體水平上的水分利用效率、產(chǎn)量水平上水分利用效率給出了明確計(jì)算公式。在政府政策制定及農(nóng)戶行為方面，王金霞等［7］從農(nóng)戶行為角度研究中國(guó)北方地區(qū)農(nóng)戶對(duì)日益短缺的水資源的反映。張仁田等［8］從水價(jià)制定角度建立了不同水價(jià)制定方法與水資源效率及公平性的數(shù)學(xué)關(guān)系模型，但是定價(jià)方法僅由土地面積影響水資源公平性。還有學(xué)者采用指標(biāo)體系評(píng)價(jià)法、數(shù)理模型、計(jì)量模型研究農(nóng)戶灌溉效率，崔遠(yuǎn)來(lái)、熊佳［9］通過(guò)回顧灌溉水利用效率評(píng)價(jià)指標(biāo)的發(fā)展歷史，詳細(xì)分析了效率評(píng)價(jià)指標(biāo)的不足，研究發(fā)現(xiàn)，目前情況下還不存在一套任何條件下都適用的灌溉水利用效率評(píng)價(jià)指標(biāo)。王曉娟、李周［10］以河北省石津灌區(qū)為例，采用超越對(duì)數(shù)隨機(jī)前沿生產(chǎn)函數(shù)和農(nóng)戶調(diào)查資料，對(duì)灌區(qū)農(nóng)業(yè)用水效率及影響因素進(jìn)行了實(shí)證研究。趙連閣、王學(xué)淵［11］根據(jù)2007年甘肅省玉內(nèi)蒙古自治區(qū)包頭市典型灌區(qū)的農(nóng)戶調(diào)查數(shù)據(jù)，運(yùn)用數(shù)據(jù)包絡(luò)分析方法計(jì)算并評(píng)價(jià)了農(nóng)戶灌溉用水效率。

綜上所述，國(guó)內(nèi)外學(xué)者普遍關(guān)注灌溉效率問(wèn)題，對(duì)水資源利用效率的研究已比較成熟，但是以精細(xì)作物蔬菜為研究對(duì)象來(lái)分析農(nóng)戶行為及灌溉用水效率的文章尚少。本研究在國(guó)內(nèi)外學(xué)者研究的基礎(chǔ)上，運(yùn)用效率評(píng)價(jià)工具隨機(jī)前沿分析模型實(shí)證研究山東省肥城縣蔬菜產(chǎn)業(yè)農(nóng)戶灌溉用水效率。肥城縣作為山東省蔬菜產(chǎn)業(yè)的一個(gè)代表性生產(chǎn)基地，其研究結(jié)果亦適用于整個(gè)黃河下游水資源短缺地區(qū)。

二、數(shù)據(jù)及模型方法

(一)實(shí)地調(diào)研情況

本研究應(yīng)用的數(shù)據(jù)來(lái)自對(duì)山東省肥城縣王莊鎮(zhèn)的實(shí)地調(diào)查。肥城縣地屬溫帶季風(fēng)氣候，氣候溫暖、四季分明，光照條件較好，農(nóng)業(yè)基礎(chǔ)比較好，是全國(guó)的糧食大縣，也是山東省頗具規(guī)模的蔬菜產(chǎn)業(yè)基地。調(diào)研地區(qū)王莊鎮(zhèn)位于肥城縣西南部，總面積91.44×106m2，經(jīng)過(guò)十幾年蔬菜種植的探索，目前80%的農(nóng)戶種植馬鈴薯、白菜，同時(shí)套種玉米，這種種植方式科學(xué)合理，收益高，已被廣大農(nóng)戶所接受。王莊鎮(zhèn)馬鈴薯種植面積超過(guò)66.7×106m2，已經(jīng)形成相當(dāng)?shù)囊?guī)模經(jīng)濟(jì)效益。王莊鎮(zhèn)的灌溉水源以地下水為主，采用明渠輸水，其中以土渠為主。由于先進(jìn)的灌溉技術(shù)在該地區(qū)還沒有得到廣泛普及，因此王莊鎮(zhèn)渠道的輸水狀況不容樂觀，渠道滲漏嚴(yán)重，水資源在運(yùn)輸過(guò)程中產(chǎn)生大量浪費(fèi)。

本研究采用隨機(jī)抽樣調(diào)查方式，在王莊鎮(zhèn)下屬村莊前于村、白屯村隨機(jī)抽取了112個(gè)農(nóng)戶進(jìn)行調(diào)查，其中前于村55戶，白屯村57戶。表1為此次調(diào)研的基本情況。

表1 馬鈴薯投入產(chǎn)出變量的統(tǒng)計(jì)描述

調(diào)研結(jié)果顯示，王莊鎮(zhèn)馬鈴薯平均產(chǎn)量為2790kg/667m2，而我國(guó)馬鈴薯平均產(chǎn)量為2000kg/667m2左右，王莊鎮(zhèn)馬鈴薯產(chǎn)量明顯高于全國(guó)平均水平。這主要是因?yàn)楫?dāng)?shù)剞r(nóng)戶在馬鈴薯種植過(guò)程中，不僅注重優(yōu)良品種的選擇，而且采用先進(jìn)的地膜覆蓋種植技術(shù)及科學(xué)的灌溉方式。

除此之外，化肥、農(nóng)藥的大量投入也是實(shí)現(xiàn)馬鈴薯高產(chǎn)的主要原因。值得注意的是灌溉成本很高，由于王莊鎮(zhèn)主要通過(guò)抽取地下水的方式進(jìn)行灌溉，農(nóng)民用于灌溉的支出幾乎已經(jīng)占到了總成本投入的20%，水資源的稀缺性已經(jīng)凸顯。表1還顯示，由于青年勞動(dòng)力轉(zhuǎn)向城市，王莊鎮(zhèn)主要的勞動(dòng)力是中老年人，每戶一般有2個(gè)勞動(dòng)力，3個(gè)勞動(dòng)力的農(nóng)戶僅占少數(shù)?？梢娹r(nóng)業(yè)發(fā)達(dá)的山東蔬菜基地同樣存在農(nóng)村青壯年勞力流失的問(wèn)題。由于本研究的重點(diǎn)在于水資源的利用，我們以每667m2土地灌溉用水的資本投入為變量衡量水資源利用效率。

(二)模型方法

1.隨機(jī)前沿分析模型

研究農(nóng)戶灌溉用水效率，可以通過(guò)設(shè)立相應(yīng)評(píng)價(jià)指標(biāo)分析灌溉水利用效率，但是指標(biāo)評(píng)價(jià)法在實(shí)際操作中可能面臨不易操作或者指標(biāo)之間存在相互交叉的現(xiàn)象，而且數(shù)據(jù)不易獲取，這必然會(huì)影響研究結(jié)果。通過(guò)隨機(jī)前沿分析方法(Stochastic Frontier Analysis，簡(jiǎn)稱 SFA)建立隨機(jī)前沿模型能夠把類似于氣候、地理等隨機(jī)游走因素與技術(shù)效率對(duì)產(chǎn)出影響的部分分離出來(lái)，可以避免誤差影響，測(cè)定效率值比較直觀。自SFA理論被提出后，眾多學(xué)者對(duì)其進(jìn)行了深入研究，產(chǎn)生了很多模型和估計(jì)處理方法，但是最常用的模型是巴特斯與科里(Battese、Coeili，1992)［12］于 1992年發(fā)展的模型(簡(jiǎn)稱 BC－1992模型)和兩人于1995年發(fā)展的模型，簡(jiǎn)稱BC － 1995 模型(Battese、Coeili，1995)［13］。簡(jiǎn)要而言，BC－1992模型函數(shù)設(shè)定如下所示:

公式(1)中，N為樣本個(gè)數(shù)，yi為樣本中第i個(gè)個(gè)體的產(chǎn)出，xi為該個(gè)體的投入量，β為待估參數(shù)向量。Vi代表隨機(jī)誤差項(xiàng)，通常被假定服從IID(獨(dú)立一致分布)，且服從正態(tài)分布N(0ui為管理誤差項(xiàng)，被用來(lái)解釋技術(shù)非效率影響，ui也屬于 IID，且服從正態(tài)分布 N(0)。Vi與 ui共同構(gòu)成了誤差項(xiàng)，并且兩者之間相互獨(dú)立。除此之外，為了方便實(shí)踐應(yīng)用，設(shè)定了一個(gè)參數(shù)BC－1992模型規(guī)定，當(dāng)γ=0時(shí)，同時(shí)意味著=0，此時(shí)，樣本中便不存在非效率狀態(tài)，也就是說(shuō)，在既定的投入下實(shí)現(xiàn)了產(chǎn)出最大化或者達(dá)到目前的產(chǎn)出水平時(shí)實(shí)現(xiàn)了投入最小化，應(yīng)用OLS方法即可進(jìn)行相關(guān)分析。反之，當(dāng)樣本中存在技術(shù)非有效率，意味著既定投入沒有實(shí)現(xiàn)產(chǎn)出最大化或者達(dá)到目前的產(chǎn)出水平?jīng)]有實(shí)現(xiàn)投入最小化，此時(shí)使用SFA就非常必要了。SFA方法的優(yōu)勢(shì)在于不僅能夠告訴樣本中每個(gè)個(gè)體的效率值，而且還能就個(gè)體間效率差異進(jìn)行定量分析［14］。

2.方法及變量設(shè)計(jì)

本研究選用SFA方法，采用BC－1992模型對(duì)農(nóng)戶灌溉用水效率進(jìn)行分析。選取單位面積灌溉用水成本、化肥、農(nóng)藥的資金投入、勞動(dòng)力人數(shù)作為影響產(chǎn)出的投入變量，馬鈴薯每667m2土地的產(chǎn)量作為產(chǎn)出變量，利用 Frontier4.1軟件分別對(duì)C－D生產(chǎn)函數(shù)和超對(duì)數(shù)生產(chǎn)函數(shù)進(jìn)行估計(jì)，SFA模型結(jié)果表明，C－D生產(chǎn)函數(shù)擬合度更好。因此，本研究采用C－D隨機(jī)前沿生產(chǎn)函數(shù)測(cè)定農(nóng)業(yè)生產(chǎn)技術(shù)效率。C－D生產(chǎn)函數(shù)，可以表示為如下方式:

(2)式中i和t分別表示農(nóng)戶的序號(hào)和年份，這里t=2010，表示年份是2010年;Y表示單位面積馬鈴薯的產(chǎn)量;X表示單位面積馬鈴薯的資本投入，包括化肥、農(nóng)藥的投入成本;L表示單位面積的勞動(dòng)力投入;W表示單位面積馬鈴薯灌溉用水的資本投入;Vit表示農(nóng)戶不可控的影響因素，服從正態(tài)分布;Uit表示技術(shù)效率損失，服從半正態(tài)分布，Vit與Uit是相互獨(dú)立的;β0、β1、β2、β3為估計(jì)參數(shù)。通過(guò)公式(2)，農(nóng)業(yè)生產(chǎn)技術(shù)效率的估計(jì)公式為:

農(nóng)業(yè)生產(chǎn)技術(shù)效率衡量的是實(shí)際產(chǎn)出與生產(chǎn)可能性邊界產(chǎn)出的比值，所謂農(nóng)戶灌溉用水效率就是指產(chǎn)出和其他投入已確定情況下，可能達(dá)到的最小灌溉水量與實(shí)際灌溉用水量之比。

參照其他學(xué)者的科研成果及公式(3)，進(jìn)一步建立如下灌溉用水效率模型:

三、研究結(jié)果分析

基于所調(diào)查的112戶農(nóng)戶馬鈴薯投入產(chǎn)出面板數(shù)據(jù)，應(yīng)用軟件Frontier4.1對(duì)模型(2)進(jìn)行估計(jì)，模型結(jié)果見表2。

表2 模型(2)的估計(jì)結(jié)果

表2中給出了生產(chǎn)函數(shù)的估計(jì)參數(shù)值，資本投入、勞動(dòng)力投入及灌溉用水投入的參數(shù)都是正值，與預(yù)期結(jié)果一致，資本、勞動(dòng)力投入對(duì)農(nóng)業(yè)生產(chǎn)有正影響。γ的值為0.7457，在1%的顯著水平下不為0，說(shuō)明農(nóng)業(yè)生產(chǎn)過(guò)程中的低效率主要來(lái)自于技術(shù)無(wú)效，而且技術(shù)無(wú)效對(duì)農(nóng)業(yè)產(chǎn)出有顯著影響，其余部分為農(nóng)戶不可控因素造成的，且影響較小?？梢姴捎肧FA方法計(jì)算農(nóng)業(yè)生產(chǎn)技術(shù)效率是非常有必要的。

根據(jù)公式(3)與(4)，基于112個(gè)農(nóng)戶數(shù)據(jù)，運(yùn)用Frontier4.1軟件對(duì)其生產(chǎn)技術(shù)效率與灌溉用水效率進(jìn)行測(cè)算，結(jié)果如表3、表4所示。

表3 農(nóng)戶生產(chǎn)技術(shù)效率與灌溉用水效率頻度分布表

表3顯示，與農(nóng)戶生產(chǎn)技術(shù)效率相比，灌溉用水效率值遠(yuǎn)低于生產(chǎn)技術(shù)效率值。所有的農(nóng)戶灌溉用水效率值都小于1，表明該地區(qū)農(nóng)業(yè)用水處于低效率狀態(tài)，都存在一定的提高空間。而且，農(nóng)戶的灌溉水利用效率最大值為91.54%，最小值僅為31.80%，而農(nóng)戶之間生產(chǎn)技術(shù)效率相差不大，灌溉用水效率與生產(chǎn)技術(shù)效率相比有更大的提高空間。通過(guò)表3我們也注意到，王莊鎮(zhèn)馬鈴薯灌溉過(guò)程中，8.04%的農(nóng)戶灌溉用水效率低于50%，也就是說(shuō)這些農(nóng)戶在灌溉過(guò)程中浪費(fèi)了50%以上的水資源。更有3.57%的農(nóng)戶灌溉用水效率低于40%，這與王莊鎮(zhèn)主要采用土渠灌溉有很大關(guān)系。

由表4可以看出，調(diào)研地區(qū)農(nóng)業(yè)生產(chǎn)技術(shù)效率平均值為94.61%，這與化肥農(nóng)藥的大量投入及山東大力推進(jìn)農(nóng)業(yè)機(jī)械化程度有很大關(guān)系?？茖W(xué)種植同樣起到了很大作用，馬鈴薯種植過(guò)程中普遍與其他糧食作物進(jìn)行田間套種，并且山東省95%的馬鈴薯采用地膜覆蓋種植，地膜的使用在抗寒、抗病等方面有積極作用，不僅提高了馬鈴薯的產(chǎn)量，還使馬鈴薯生長(zhǎng)期提前，搶先占領(lǐng)市場(chǎng)，增加了農(nóng)戶經(jīng)濟(jì)收益。所有這些技術(shù)因素導(dǎo)致了山東省農(nóng)戶生產(chǎn)技術(shù)效率普遍較高。本研究中的王莊鎮(zhèn)農(nóng)業(yè)生產(chǎn)平均技術(shù)效率已經(jīng)達(dá)到了94.61%，成果顯著。反觀農(nóng)戶灌溉用水效率，我們發(fā)現(xiàn)，所調(diào)研的112戶農(nóng)戶的灌溉用水效率平均值為70.08%，說(shuō)明就目前技術(shù)條件下，在其他投入保持不變的情況下，達(dá)到當(dāng)前馬鈴薯產(chǎn)量可以減少29.92%灌溉用水投入，即該地區(qū)灌溉用水具有顯著節(jié)約潛力。

表4 農(nóng)戶生產(chǎn)技術(shù)效率與灌溉用水效率分布表

通過(guò)對(duì)王莊鎮(zhèn)馬鈴薯灌溉水利用現(xiàn)狀的模型分析，我們認(rèn)為山東省蔬菜種植過(guò)程中存在的主要問(wèn)題是:(1)與生產(chǎn)技術(shù)效率相比，王莊鎮(zhèn)灌溉用水效率還有相當(dāng)?shù)奶嵘臻g;(2)土渠的普遍使用導(dǎo)致輸水過(guò)程中的大量水資源浪費(fèi)，這是灌溉水利用效率低下的原因之一;農(nóng)戶田間灌溉主要使用畦灌方式，先進(jìn)的節(jié)水灌溉技術(shù)使用較少，諸多原因可能導(dǎo)致灌溉水利用效率低下;(3)農(nóng)戶的節(jié)水意識(shí)薄弱、資金投入不足也影響了灌溉水利用效率。

四、政策建議

基于以上模型結(jié)果分析，結(jié)合山東省馬鈴薯灌溉水利用現(xiàn)狀，本研究認(rèn)為可以從以下幾方面提高蔬菜灌溉用水效率，緩解水資源短缺壓力，在有限的水資源總量情況下，實(shí)現(xiàn)水資源的最優(yōu)配置和有效利用。

1.發(fā)展精準(zhǔn)農(nóng)業(yè)技術(shù)，建設(shè)科技農(nóng)業(yè)示范田

山東省蔬菜產(chǎn)業(yè)規(guī)模巨大且日趨成熟，推廣以3S(全球定位系統(tǒng)(GPS)、地理信息系統(tǒng)(GIS)、遙感系統(tǒng)(RS))技術(shù)為基礎(chǔ)的精準(zhǔn)農(nóng)業(yè)成為提升蔬菜產(chǎn)業(yè)競(jìng)爭(zhēng)力的有效途徑。根據(jù)目前蔬菜產(chǎn)業(yè)發(fā)展現(xiàn)狀，建立以精準(zhǔn)技術(shù)為基礎(chǔ)的農(nóng)業(yè)示范田，運(yùn)用信息技術(shù)及先進(jìn)的管理模式充分挖掘農(nóng)田的生產(chǎn)潛力、提升農(nóng)業(yè)的科技水平;通過(guò)推廣精準(zhǔn)農(nóng)業(yè)技術(shù)，大幅提高蔬菜產(chǎn)量和質(zhì)量。同時(shí)組織動(dòng)員農(nóng)戶到示范田參觀學(xué)習(xí)，引導(dǎo)農(nóng)戶學(xué)習(xí)掌握先進(jìn)的3S技術(shù)及高效的種田方法，加快新技術(shù)的應(yīng)用速度，深化新品種的普及程度。通過(guò)推廣科學(xué)生產(chǎn)，帶動(dòng)普通農(nóng)戶科技興農(nóng)的意識(shí)和積極性，實(shí)現(xiàn)蔬菜產(chǎn)業(yè)的高效生產(chǎn)。

2.明確水權(quán)，通過(guò)適度提高水價(jià)達(dá)到促使農(nóng)戶節(jié)水的目的

通過(guò)明確的法律文件規(guī)定水權(quán)使用，防止水資源濫用及過(guò)度開采，加強(qiáng)水權(quán)、水法的宣傳力度，培養(yǎng)農(nóng)民合理、正確利用水資源的意識(shí)。在宣傳法律法規(guī)同時(shí)，適當(dāng)提高水價(jià)，把水資源作為經(jīng)濟(jì)物品，按灌溉用水量進(jìn)行收費(fèi)，以價(jià)格機(jī)制促使農(nóng)民意識(shí)到水資源的稀缺性，達(dá)到節(jié)水目的。同時(shí)，鼓勵(lì)農(nóng)戶成立水資源利用協(xié)會(huì)，加強(qiáng)水資源利用的監(jiān)督，防止過(guò)度開采、浪費(fèi)水資源及水資源的不公平分配，提高灌溉用水效率。

3.實(shí)現(xiàn)灌溉效率提高的同時(shí)注重公平分配

在水資源短缺及利用效率低下的現(xiàn)狀下，可把農(nóng)戶灌溉效率評(píng)價(jià)指標(biāo)與經(jīng)濟(jì)發(fā)展指標(biāo)(GDP)同時(shí)納入村干部考核機(jī)制，可以更好地促使村干部加強(qiáng)水資源的合理配置，促使政策制定者在追求經(jīng)濟(jì)發(fā)展的同時(shí)也重視提高灌溉用水效率，實(shí)現(xiàn)水資源有效利用。結(jié)合各地區(qū)實(shí)際情況制定灌溉效率目標(biāo)，將此目標(biāo)作為獎(jiǎng)勵(lì)村干部政績(jī)的一部分，可有效促進(jìn)農(nóng)村水資源利用效率。同時(shí)，在農(nóng)田灌溉過(guò)程中，應(yīng)該加強(qiáng)監(jiān)督、制定透明機(jī)制，公布灌溉明細(xì)，通過(guò)公眾監(jiān)督實(shí)現(xiàn)水資源的公平分配。建立農(nóng)村用水協(xié)會(huì)，對(duì)灌溉水量進(jìn)行分配、監(jiān)督，確保公平、合理、有效利用水資源。

［1］Lester R.Brown，Brian Halweil.China’s Water Shortage Could Shake World Grain Markets［R］ .World watch Institute，Washington.C.，1998.

［2］Karagiannis，G.;Tzouvelekas，V.Xepapadeas，A.:Measuring Irrigation Water Efficiency with a Stochastic Production Frontier［J］.Environmental and Resource Economics，2003，26:57－72.

［3］Dhehibi B，Lachaal L，Elloumi M，etal.Measuring irrigation water use efficiency using stochastic production frontier:An application on citrus producing farms in Tunisia［J］.African Journal of Agricultural and Resource E-conomics，2007，1(2):1 －15.

［4］Stijn Speelman，Marijke D＇Haese，Jeroen Buysse，etal.A measure for the efficiency of water use and its determinants，a case study of small－ scale irrigation schemes in North－ West Province，South Africa［J］.Agricultural Systems，2008:31 －39.

［5］Abeer MunlaHasan .Water Use Efficiency in Syrian Agriculture［J］.NAPC，2001.1:1 －11.

［6］王會(huì)肖，劉昌明.作物水分利用效率內(nèi)涵及研究進(jìn)展［J］.水科學(xué)進(jìn)展，2000，11(1):99 －104.

［7］王金霞，黃季焜，張麗娟等.北方地區(qū)農(nóng)民對(duì)水資源短缺的反應(yīng)［J］.水利經(jīng)濟(jì)，2008，26(5):1 －5.

［8］張仁田，童立忠，陸小偉等.水定價(jià)方法對(duì)農(nóng)業(yè)用水效率及公平性的影響分析［J］.水利科技與經(jīng)濟(jì)，2005，11(5).257 －261.

［9］崔遠(yuǎn)來(lái)，熊佳.灌溉水利用效率指標(biāo)研究進(jìn)展［J］.水科學(xué)進(jìn)展，2009，20(4):590 －598.

［10］王曉娟，李周.灌溉用水效率及影響因素分析［J］.中國(guó)農(nóng)村經(jīng)濟(jì)，2005，7:11 －18.

［11］趙連閣，王學(xué)淵.農(nóng)戶灌溉用水的效率差異——基于甘肅、內(nèi)蒙古兩個(gè)典型灌區(qū)實(shí)地調(diào)查的比較分析［J］.農(nóng)業(yè)經(jīng)濟(jì)問(wèn)題，2010，3:71 －78.

［12］Battese and Coelli battese，G.E.and T.J.Coelli.Frontier Production Functions，Technical Efficiency and Panel Data.:With Application to Paddy Farmers in India［J］.Journal of Productivity Analysis，June 1992，3，153－169.

［13］Battese and Coelli battese，G.E.and T.J.Coelli.A model for Technical Inefficiency Effects in a Stochastic Frontier for Panel Data.［J］ Empirical Economics，1995，Vol.20，325 －332.

［14］水資源經(jīng)濟(jì)學(xué)［M］.北京:中國(guó)環(huán)境科學(xué)出版社，2008:110－111.