基于隨機前沿生產(chǎn)函數(shù)法的飼草生產(chǎn)技術(shù)效率測算及影響因素分析——以西藏“一江兩河”地區(qū)為例

曲云鶴，余成群，孫維，武俊喜，李少偉

(中國科學(xué)院地理科學(xué)與資源研究所，生態(tài)系統(tǒng)網(wǎng)絡(luò)觀測與模擬重點實驗室，北京 100101)

基于隨機前沿生產(chǎn)函數(shù)法的飼草生產(chǎn)技術(shù)效率測算及影響因素分析——以西藏“一江兩河”地區(qū)為例

曲云鶴*，余成群，孫維，武俊喜，李少偉

(中國科學(xué)院地理科學(xué)與資源研究所，生態(tài)系統(tǒng)網(wǎng)絡(luò)觀測與模擬重點實驗室，北京 100101)

摘要：以西藏“一江兩河”地區(qū)農(nóng)戶調(diào)研數(shù)據(jù)為基礎(chǔ)，運用隨機前沿生產(chǎn)函數(shù)估算農(nóng)戶飼草生產(chǎn)技術(shù)效率及影響因素，并從飼草類別、農(nóng)業(yè)規(guī)模和地區(qū)分布3個維度進一步分析技術(shù)效率的增產(chǎn)空間。結(jié)果顯示，飼草生產(chǎn)存在顯著效率損失，提高生產(chǎn)效率，單產(chǎn)還有近30%的增長空間；投入要素勞動、種子、燃料、農(nóng)家肥、化肥和土地產(chǎn)出彈性分別為0.1423，0.2197，0.1309，0.2600，0.2265和0.0206；品種差異是飼草技術(shù)效率最主要影響因素，農(nóng)業(yè)規(guī)模和流動資金要素對飼草生產(chǎn)效率提高產(chǎn)生間接負(fù)向影響，人力資本和農(nóng)業(yè)綜合生產(chǎn)能力要素對飼草生產(chǎn)效率影響不顯著。釋放飼草生產(chǎn)技術(shù)效率增產(chǎn)空間，可采取推廣普及豆科飼草生產(chǎn)技術(shù)、擴大飼草種植規(guī)模、各地區(qū)(特別是拉薩地區(qū))進一步明確適宜引種的飼草品種等策略。

關(guān)鍵詞：隨機前沿生產(chǎn)函數(shù)；飼草生產(chǎn)技術(shù)效率；西藏“一江兩河”地區(qū)

Use of the stochastic frontier production function in a technical efficiency analysis of grass production: An example from the “One River and Two Tributaries” region of Tibet

QU Yun-He*, YU Cheng-Qun, SUN Wei, WU Jun-Xi, LI Shao-Wei

KeyLaboratoryofEcosystemNetworkObservationandModeling,InstituteofGeographic,SciencesandNaturalResourcesResearch,ChineseAcademyofSciences,Beijing100101,China

Abstract:Based on the agricultural household on-the-spot investigation data collected in the “One River and Two Tributaries” Region of Tibet, we estimated the technical efficiency of grass production of agricultural households using the stochastic frontier production function. This method estimates the contributions of grass family differences, farm size, and region distribution to the technical efficiency of grass production. The results indicated that there is a loss of technical efficiency in grass production in the “One River and Two Tributaries” Region of Tibet. The yield per mu could increase by approximately 30% with improved technical efficiency. The output elasticity coefficients with respect to input factors of labor, seeds, fuel, organic fertilizer, chemical fertilizer, and land were 0.1423, 0.2197, 0.1309, 0.2600, 0.2265, and 0.0206, respectively. Grass family difference was identified as the key factor in improving technical efficiency. Farm size and liquidity strongly affected the technical efficiency of grass production via their effects on other input factors such as labor and capital. This was because of the disadvantageous position of resource use, compared with animal husbandry and off-farm industries. Human capital and agricultural comprehensive production capability did not affect the technical efficiency of grass production. Increasing the grass yield by improving technical efficiency can include measures such as extending the production technology to legume forage grasses, enlarging grass planting areas, and identifying which grass families are best adapted for growth in different regions, especially Lhasa.

Key words:stochastic frontier production function; technical efficiency of grass production; “One River and Two Tributaries” Region in Tibet

我國從20世紀(jì)90年代初開始發(fā)展糧食作物—經(jīng)濟作物—飼料作物三元種植結(jié)構(gòu)，經(jīng)過20多年發(fā)展，2005年我國青飼料種植面積達最大值337.4萬hm2，占總耕地面積的2.17%，近年來面積有所降低，但仍保持在200萬hm2以上[1]。我國飼料供需總體基本平衡，但結(jié)構(gòu)差異較大，籽食性飼料供求平衡，豐年有余，但營養(yǎng)體飼草料供應(yīng)缺口逐年加大，尤其是優(yōu)質(zhì)飼草料(如紫花苜蓿)對外依存度逐年提高，缺口達66.7%[1]。我國飼草供應(yīng)不足形勢嚴(yán)峻。大力發(fā)展人工種植牧草是解決我國畜牧業(yè)生產(chǎn)瓶頸問題的必然方向[2]，也是解決草畜供應(yīng)不平衡矛盾的主要途徑[3]。作為我國五大牧區(qū)之一的西藏，近年來地方政府加大科技、人才和資金投入力度，花大力氣扶持飼草業(yè)發(fā)展，為我國高原飼草業(yè)發(fā)展奠定良好基礎(chǔ)。西藏“一江兩河”地區(qū)(拉薩河、年楚河和雅魯藏布江河谷地區(qū))是全區(qū)農(nóng)牧業(yè)發(fā)展的核心地帶，也是人口主要聚居區(qū)。近年來，隨著地方政府引導(dǎo)和扶持力度加大，以及退牧還草、草原生態(tài)保護補助獎勵機制等國家政策措施持續(xù)推進，西藏畜牧業(yè)生產(chǎn)方式由依賴天然牧場靠天養(yǎng)畜向舍飼和半舍飼方式轉(zhuǎn)變步伐加快，畜牧業(yè)專業(yè)化、集約化、商品化程度顯著提高，農(nóng)牧民牲畜養(yǎng)殖對天然飼草依賴度降低，對種植飼草需求逐年增加。再加上西藏災(zāi)害性天氣爆發(fā)頻繁，調(diào)節(jié)季節(jié)余缺和防抗災(zāi)儲備飼草需求也不斷增加。從長遠(yuǎn)發(fā)展以及節(jié)省財政資金角度考慮，西藏政府把飼草供應(yīng)的基點立足于區(qū)內(nèi)自給，20世紀(jì)70年代開始從區(qū)外引種適宜飼草[4]，以項目為依托開展人工種草，在“一江兩河”地區(qū)構(gòu)建三元種植結(jié)構(gòu)，把工作重點轉(zhuǎn)到為農(nóng)牧民提供優(yōu)質(zhì)草種、種植技術(shù)以及鼓勵農(nóng)牧民種植飼草上?！段鞑亟y(tǒng)計年鑒》數(shù)據(jù)顯示，1998年以來青飼料年產(chǎn)量以10%以上速度增長，2012年產(chǎn)量達到31.6萬t。但全區(qū)飼草料仍供不應(yīng)求，優(yōu)質(zhì)飼草料缺口達50%[5-6]。目前，在政府引導(dǎo)下，西藏農(nóng)牧民飼草種植意愿顯著提高，但出于口糧需求考慮，農(nóng)牧民不愿意把好地用來種草。為確保地區(qū)糧食安全，避免與糧食作物爭地，在資源硬約束條件下，發(fā)揮農(nóng)牧民生產(chǎn)潛力，提高飼草生產(chǎn)效率就成為挖掘西藏飼草增產(chǎn)潛力的可行途徑。

國內(nèi)外對飼草的研究多集中在飼草理化生長習(xí)性上[3,6-7]，而從生產(chǎn)要素資源配置及生產(chǎn)效率角度的研究還不多見。生產(chǎn)效率提高可通過技術(shù)進步和技術(shù)效率改進實現(xiàn)[8-9]。從投入產(chǎn)出曲線來看，技術(shù)進步指整個生產(chǎn)前沿(一定技術(shù)水平下能夠?qū)崿F(xiàn)的最大產(chǎn)出曲線)的外移，而技術(shù)效率則是指向生產(chǎn)前沿的逼近，越接近生產(chǎn)前沿，說明技術(shù)效率越高[10]。農(nóng)戶飼草生產(chǎn)技術(shù)會隨地區(qū)社會經(jīng)濟發(fā)展而自然提高，技術(shù)進步具有全局性和不可逆性，而技術(shù)效率往往受到農(nóng)戶個體生產(chǎn)能力、地區(qū)自然狀況和政府扶持等因素影響，農(nóng)戶個體技術(shù)效率差異較大；與提高技術(shù)進步相比，通過平衡農(nóng)戶技術(shù)效率差異提高飼草產(chǎn)量更易于實現(xiàn)。因此，在現(xiàn)有投入產(chǎn)出情況下，以飼草生產(chǎn)技術(shù)效率高的農(nóng)戶為樣板優(yōu)化農(nóng)戶生產(chǎn)資源配置，以及找出效率影響因素，控制不利因素的發(fā)生是本文重點關(guān)注和研究的內(nèi)容。

生產(chǎn)技術(shù)效率測算的隨機前沿生產(chǎn)函數(shù)理論最早由Meeusen和Broeck[11]、Lovell和Schmidt[12]于1977年分別提出。初期采用兩步法來測算，由于很難確保技術(shù)無效率變量獨立且同分布假定存在，被1993年Battese和Coelli[13]提出的一步法所替代。一步法就是同時估計隨機前沿生產(chǎn)函數(shù)和技術(shù)效率影響因素的方法，從而避免兩步法可能出現(xiàn)的有偏估計。目前，這一理論已經(jīng)相對成熟，被國內(nèi)外學(xué)者廣泛用于農(nóng)業(yè)領(lǐng)域應(yīng)用研究中。20世紀(jì)90年代初，這一方法被應(yīng)用于我國農(nóng)業(yè)生產(chǎn)領(lǐng)域，特別是糧食作物生產(chǎn)技術(shù)效率測算，隨后推廣應(yīng)用到經(jīng)濟作物(如蔬菜、水果、茶葉、煙草等)以及畜禽養(yǎng)殖等領(lǐng)域，但在飼草料生產(chǎn)領(lǐng)域的研究還不多[14-19]。學(xué)者們從時空、地區(qū)發(fā)展、氣候環(huán)境、品種、規(guī)模等維度，研究我國農(nóng)業(yè)生產(chǎn)技術(shù)效率發(fā)展變化規(guī)律，探尋提高農(nóng)產(chǎn)品生產(chǎn)技術(shù)效率的路徑。從宏觀層面來看，我國農(nóng)業(yè)生產(chǎn)技術(shù)效率在1979至2005年間先增加后降低[20]；東部地區(qū)高、中西部地區(qū)低[21]；適宜作物生長的主產(chǎn)區(qū)比非主產(chǎn)區(qū)技術(shù)效率高[22]；糧食作物、經(jīng)濟作物和畜產(chǎn)品生產(chǎn)都存在技術(shù)效率損失[23-26]。從微觀層面來看，經(jīng)營規(guī)模是農(nóng)戶生產(chǎn)技術(shù)效率的限定因素，園藝作物規(guī)模與技術(shù)效率呈“倒U型”關(guān)系[27]，油料作物規(guī)模與技術(shù)效率呈“U型”關(guān)系[28]。本研究擬在國內(nèi)外現(xiàn)有研究基礎(chǔ)上，進一步探尋西藏農(nóng)戶種植飼草的生產(chǎn)技術(shù)效率變化規(guī)律及影響因素，為促進飼草產(chǎn)量提高提供技術(shù)支持和指導(dǎo)。

1材料與方法

1.1　數(shù)據(jù)來源說明

本文采用2011年中國科學(xué)院拉薩農(nóng)業(yè)生態(tài)試驗站西藏高原草業(yè)工程技術(shù)研究中心對西藏“一江兩河”地區(qū)覆蓋拉薩、日喀則、山南11個村鎮(zhèn)農(nóng)戶農(nóng)牧業(yè)生產(chǎn)和消費狀況實地調(diào)研數(shù)據(jù)，樣本覆蓋了主要引種飼草品種的示范推廣地區(qū)。采用隨機抽樣的方法入戶調(diào)研299戶農(nóng)民家庭，樣本蓋度平均為16.71%(在4.37%到49.03%之間)，并從中抽調(diào)出年度種植飼草的185戶為研究對象(表1)。樣本農(nóng)戶種植的飼草品種涵蓋了燕麥、紫花苜蓿、箭筈豌豆和青飼玉米4種西藏重點引種的人工飼草品種，樣本值分別為85、50、49和1。由于青飼玉米種植戶數(shù)量較少缺乏代表性，本文剔除這一數(shù)據(jù)。

表1　樣本分布地區(qū)和數(shù)量

1.2　理論模型

本研究采用1993年Battese和Coelli[13,29]提出的基于隨機前沿生產(chǎn)函數(shù)的生產(chǎn)技術(shù)效率測算方法，理論模型表示如下：

Yit=Xitβit+Vit-Uiti=1,…,N，t=1,…,T

(1)

式中，Yit為產(chǎn)量，Xit為投入要素，βit為待估參數(shù)，Vit為代表誤差項的隨機變量，與Uit完全獨立，滿足正態(tài)分布N(0，σv2)。Uit為用來解釋生產(chǎn)中技術(shù)無效率的非負(fù)隨機變量，相互獨立且服從期望為zitδit、方差為σu2在0點截取的半正態(tài)分布N+(zitδit，σu2)。zit為無效影響因素，δit為待估參數(shù)。

(1)式殘差項不滿足獨立同分布假定，不適合用最小二乘法估計，Battese和Coelli[13,29]根據(jù)γ=σu2/(σv2+σu2) 中兩個方差變化與技術(shù)無效率之間的關(guān)系建立了似然函數(shù)，采用最大似然法估計各參數(shù)，并根據(jù)技術(shù)效率表達式(2)，估計單位i技術(shù)效率。

TEit=E(Yit*︱Uit,Xit)/E(Yit*︱Ui=0,Xit)

(2)

式中,Yit*為i單位t期的產(chǎn)量。TEit取值范圍為從0到1。(1)式采用對數(shù)形式時，Yit*為EXP(Yit)，TEit為EXP(-Uit)。當(dāng)Uit的期望為0時，i單位技術(shù)效率最大，即TEit等于1。Battese和Coelli[13,29]把除投入要素以外的所有因素放入U中，因此，估算出的結(jié)果是總技術(shù)效率損失(gross technical inefficiency)[29]。Kumbhakar等[30]在生產(chǎn)函數(shù)中直接加入生產(chǎn)技術(shù)效率影響變量，估算出的技術(shù)效率結(jié)果是凈技術(shù)效率損失(net technical inefficiency)。本文采用Battese和Coelli[13,29]的U設(shè)定方法，因此估計出的是總技術(shù)效率損失。

1.3　實證模型

1)隨機前沿生產(chǎn)函數(shù)模型選擇

農(nóng)戶層面隨機前沿生產(chǎn)函數(shù)生產(chǎn)技術(shù)效率測算多采用柯布—道格拉斯生產(chǎn)函數(shù)(簡稱CD函數(shù))和超對數(shù)生產(chǎn)函數(shù)。后者的優(yōu)點是放寬了技術(shù)中性和產(chǎn)出彈性固定的假設(shè)，缺點是易產(chǎn)生多重共線問題。本文采用似然比檢驗確定函數(shù)形式，通過估算獲得CD函數(shù)似然值L(βR)為-212.45，超對數(shù)函數(shù)似然值L(βUR)為-190.23，似然比為44.86(LR=-2[L(βR)-L(βUR)])大于5%顯著水平下的卡方臨界值25，所以拒絕原假設(shè)，接受超對數(shù)函數(shù)形式。

2)超對數(shù)隨機前沿生產(chǎn)函數(shù)模型設(shè)定及變量說明

本研究假定農(nóng)戶飼草生產(chǎn)規(guī)模報酬不變，根據(jù)飼草單位面積投入產(chǎn)出情況，建立超對數(shù)隨機前沿生產(chǎn)函數(shù)估算飼草生產(chǎn)技術(shù)效率，模型具體設(shè)定如下：

(3)

式中，Yi表示i農(nóng)戶家庭生產(chǎn)飼草的單產(chǎn)(kg/hm2)。飼草生產(chǎn)投入要素包括勞動力(X1，d/hm2)、種子(X2，kg/hm2)、燃料(X3，元/hm2)、農(nóng)家肥(X4，kg/hm2)和化肥(X5，kg/hm2)(表2)。把農(nóng)家肥與化肥區(qū)分開，主要為比較有機農(nóng)業(yè)與化肥農(nóng)業(yè)對產(chǎn)量的差異影響。

3)技術(shù)效率損失函數(shù)，即U的設(shè)定如下：

(4)

式中，δm(m=0,1,2,…,9)和θn(n=1,2,3,…,6)為待估參數(shù)，μi為殘差項，生產(chǎn)技術(shù)效率影響因素(Zm和Dn)主要包括5個方面(表2)。

① 人力資本情況。主要考察人力資本的質(zhì)量和數(shù)量狀況，分別用家庭平均受教育年限(Z1)、代表戶主農(nóng)業(yè)生產(chǎn)經(jīng)驗的戶主年齡(Z2)、表示勞動力充裕程度的家庭人口規(guī)模(Z3)和戶主性別二元選擇變量(D1，戶主為男D1=1，否則為0)來表示。

② 品種差異。分別用燕麥、紫花苜蓿和箭筈豌豆3個農(nóng)戶二元種植變量(D4、D5和D6)衡量禾本科和豆科飼草品種差異的影響(本文不涉及農(nóng)戶交叉種植情況)。

③ 農(nóng)業(yè)規(guī)模。設(shè)定耕地面積(Z6)和牲畜飼養(yǎng)量(Z7)兩個變量來分析農(nóng)牧業(yè)生產(chǎn)規(guī)模對飼草生產(chǎn)率的影響。農(nóng)業(yè)規(guī)模對技術(shù)效率的影響是間接的[31-33]，在規(guī)模效應(yīng)不變的情況下，種植業(yè)和養(yǎng)殖業(yè)的規(guī)?；a(chǎn)與飼草生產(chǎn)將形成要素競爭。

④ 農(nóng)業(yè)綜合生產(chǎn)能力。設(shè)定小麥糧單產(chǎn)(Z5)、奶牛日均最大產(chǎn)奶量(Z4)和以前是否種植過飼草二元選擇變量(D2)共3個變量來分析農(nóng)戶飼草生產(chǎn)技術(shù)率高低情況。

⑤ 流動資金限制。設(shè)定農(nóng)戶獲得信貸與否二元變量(D3)、畜牧業(yè)收入和非農(nóng)收入占家庭收入(Z8，Z9)比重等3個變量，從不同收入來源分析它們對飼草生產(chǎn)資金投向的驅(qū)動差異。

此外，農(nóng)戶參加培訓(xùn)、參與合作社經(jīng)營情況和受災(zāi)情況理論上都是生產(chǎn)效率的重要影響因素，但調(diào)研問卷中這些指標(biāo)數(shù)據(jù)缺失較多，所以沒有納入分析中。

4)投入要素產(chǎn)出彈性，如下所示：

(5)

2結(jié)果與分析

本文采用Coelli[34]于1996年更新的Frontier 4.1軟件綜合測算飼草種植戶的生產(chǎn)技術(shù)效率，估計結(jié)果見表2。

2.1　生產(chǎn)函數(shù)模型的估計結(jié)果及分析

首先，對模型進行檢驗。估計結(jié)果顯示，大部分變量都通過了顯著性檢驗。生產(chǎn)函數(shù)參數(shù)的確定結(jié)果為變差率γ=0.5518，且通過顯著性檢驗。γ=0.5518的結(jié)果說明在生產(chǎn)技術(shù)效率的影響中，技術(shù)效率損失項μ是大于隨機誤差ε的，但是與1接近度相對較低，主要原因是研究存在缺失效率損失變量，如缺少農(nóng)戶專項培訓(xùn)和參加合作社情況變量。在變差率γ=0的零假設(shè)條件下，約束條件為1，顯著性概率為1%的χ2臨界值是30.58，模型(3)的單邊似然比檢驗統(tǒng)計量LR=63.52>30.58，變差率的零假設(shè)被拒絕，即生產(chǎn)非效率項μ是存在的，所以農(nóng)戶飼草生產(chǎn)存在效率損失情況。

其次，生產(chǎn)函數(shù)投入要素彈性分析。把估計結(jié)果顯著變量帶入式(5)獲得勞動、種子、燃料、農(nóng)家肥和化肥各要素產(chǎn)出彈性分別為0.1423，0.2197，0.1309，0.2600和0.2265。根據(jù)規(guī)模報酬不變的假定，推出土地產(chǎn)出彈性為0.0206。在各要素年增長率變化相同情況下，農(nóng)家肥和化肥對飼草增產(chǎn)貢獻最大，種子其次，勞動和燃料貢獻第三，而土地最小，表明目前西藏飼草生產(chǎn)處于技術(shù)密集型發(fā)展階段，而非勞動和土地密集型階段，這與西藏以科技帶動飼草業(yè)發(fā)展的實際情況相符合。農(nóng)家肥產(chǎn)出彈性大于化肥，農(nóng)家肥投入每增加1%，飼草產(chǎn)出較等增化肥投入多增0.0335%，表明在西藏大力發(fā)展有機飼草產(chǎn)業(yè)，發(fā)展生態(tài)農(nóng)業(yè)大有潛力可挖。

2.2　技術(shù)效率損失函數(shù)的估計結(jié)果及影響因素分析

技術(shù)效率損失函數(shù)估計結(jié)果表明：第一，人力資本數(shù)量和質(zhì)量狀況不對飼草生產(chǎn)技術(shù)效率變化產(chǎn)生顯著影響。與預(yù)期不同的是，家庭人口受教育水平(Z1)影響不顯著，提高文化教育對飼草效率改善作用不明顯；而戶主年齡變量(Z2)顯著，但與預(yù)期相反，戶主農(nóng)業(yè)生產(chǎn)經(jīng)驗豐富對飼草技術(shù)效率提高是不利的，因為年齡相對較大的戶主受傳統(tǒng)思想影響較深，相對比較保守，對新技術(shù)和新產(chǎn)品接受度較低；戶主性別變量(D1)不顯著，表明男性戶主和女性戶主在飼草生產(chǎn)決策力上不存在顯著差異；家庭規(guī)模(Z3)變量弱顯著，勞動力資源豐富的家庭對技術(shù)效率提高起到一定效果。第二，飼草類別差異是飼草技術(shù)效率重要影響因素。禾本科飼草品種(以燕麥為代表)比豆科飼草品種(以紫花苜蓿和箭筈豌豆為代表)效率高(D4顯著)，前者較后者對單產(chǎn)貢獻高2.36%。豆科以下細(xì)分品種差異，如紫花苜蓿和箭筈豌豆，對效率影響不顯著(D5和D6都不顯著)。第三，農(nóng)業(yè)生產(chǎn)規(guī)模，尤其是種植業(yè)規(guī)模是飼草技術(shù)效率重要影響因素。家庭耕地面積越小(Z6顯著)，飼草技術(shù)效率越高，與理論預(yù)期相一致，進一步驗證這一理論在飼草生產(chǎn)領(lǐng)域也適用。但牲畜飼養(yǎng)規(guī)模對效率提高無作用(Z7不顯著)。第四，農(nóng)戶農(nóng)業(yè)綜合生產(chǎn)能力不是飼草生產(chǎn)技術(shù)效率的影響因素。研究結(jié)果顯示，有種草經(jīng)驗并沒有對飼草生產(chǎn)效率提高產(chǎn)生積極影響(D2不顯著)，糧食生產(chǎn)能力對飼草生產(chǎn)不起作用(Z5不顯著)，牲畜飼養(yǎng)能力對飼草生產(chǎn)起到一定作用(Z4弱顯著)。出現(xiàn)這樣的結(jié)果，一方面表示飼草生產(chǎn)是一項技術(shù)密集型工作，農(nóng)戶需要接受專業(yè)指導(dǎo)；另一方面表示農(nóng)戶對飼草生產(chǎn)重視程度不夠。第五，農(nóng)戶資金狀況是飼草生產(chǎn)技術(shù)效率的影響因素，但影響是負(fù)面的。無論是非農(nóng)收入增加(δ9=2.8774，顯著)，還是牧業(yè)收入增加(δ8=3.1509，顯著)，都會使飼草生產(chǎn)效率降低，非農(nóng)收入和牧業(yè)收入占家庭收入比重每提高1%，將使飼草單產(chǎn)分別降低3.09%和3.39%，表明飼草業(yè)與非農(nóng)業(yè)和牧業(yè)發(fā)展在勞動和資金投入上處于資源競爭劣勢地位。另外，家庭能否獲得貸款則與飼草生產(chǎn)效率無關(guān)(D3不顯著)。

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

注：***，**，*和+分別表示1%，5%，10%和20%的顯著水平。

Notes:***,**,* and + express the significant level of one percent, five percent, ten percent and twenty percent.

2.3　生產(chǎn)技術(shù)效率測算值分析

樣本農(nóng)戶飼草生產(chǎn)技術(shù)效率測算值從7.44%到93.86%不等，平均為72.78%。如果保持現(xiàn)有技術(shù)和投入水平，進一步消除技術(shù)無效率因素影響，飼草產(chǎn)出仍有27.22%的增長空間，因此，提高生產(chǎn)技術(shù)效率可以顯著地擴大飼草供給量，增加飼草自給率。

2.3.1不同科類飼草生產(chǎn)技術(shù)效率分析根據(jù)技術(shù)效率的樣本分布，本文把技術(shù)效率分為5個區(qū)間：小于30%、大于等于30%到60%、大于等于60%到80%、大于等于80%到90%、大于等于90%。各區(qū)間樣本容量分別為13、30、51、64和37。樣本技術(shù)效率分布較分散，表明飼草生產(chǎn)技術(shù)差異較大。絕大多數(shù)禾本科飼草生產(chǎn)技術(shù)效率在80%以上，禾本科飼草農(nóng)戶生產(chǎn)技術(shù)相似度很高，但豆科飼草農(nóng)戶生產(chǎn)技術(shù)則差異很大(效率值方差較大)(圖1)。在這種情況下，要提高豆科飼草生產(chǎn)技術(shù)效率的前提條件是進一步推廣豆科生產(chǎn)技術(shù)，提高農(nóng)戶飼草生產(chǎn)管理技能。

圖1　不同品種飼草生產(chǎn)技術(shù)效率分布Fig.1　Technical efficiency distribution of different families’ grasses

2.3.2不同農(nóng)業(yè)規(guī)模農(nóng)戶技術(shù)效率分析從圖2可以看出，從整體來看，農(nóng)戶農(nóng)業(yè)規(guī)模(指種植業(yè)規(guī)模，由耕地面積表示)與生產(chǎn)技術(shù)效率并不完全呈單調(diào)的相反關(guān)系，準(zhǔn)確地說應(yīng)該呈似“N型”。對于種植面積大戶(2 hm2以上，占11.56%)和小戶(0.67 hm2以下，占25.12%)來說，兩者呈同向變動關(guān)系，而在0.667~2.000 hm2區(qū)間段，兩者呈反向變動關(guān)系。由于多數(shù)農(nóng)戶(占63.3%)處于中間段，所以回歸結(jié)果顯示兩者呈相反關(guān)系。如果放松規(guī)模報酬不變假定，飼草種植規(guī)模與生產(chǎn)技術(shù)效率也呈似“N”型關(guān)系，拐點為0.133和0.266 hm2，可見，飼草規(guī)模與技術(shù)效率的變動關(guān)系與種植業(yè)規(guī)模與技術(shù)效率變動相似。

圖2　農(nóng)戶農(nóng)業(yè)種植規(guī)模與飼草生產(chǎn)技術(shù)效率變化趨勢Fig.2　Trends of technical efficiencies and farm size

調(diào)研農(nóng)戶飼草戶均種植規(guī)模為0.093 hm2，占家庭農(nóng)作物總播種面積的7.43%；但農(nóng)戶間差異較大，種植大戶接近0.667 hm2，而種植小戶不到0.033 hm2。在耕地面積小于0.667 hm2、0.667~2.000 hm2和2.000 hm2以上3個區(qū)間段內(nèi)，農(nóng)戶家庭飼草種植規(guī)模平均為0.067,0.113和0.133 hm2，戶均技術(shù)效率分別為0.69、0.66和0.70，戶均單產(chǎn)分別為883.0,834.0和898.5 kg。如果以“N型”變動趨勢相對照，可以發(fā)現(xiàn)：不同規(guī)模農(nóng)戶飼草種植面積不合理，多數(shù)低于理論值，小規(guī)模農(nóng)戶和中等規(guī)模農(nóng)戶(88.42%的農(nóng)戶)在0.133 hm2極值點內(nèi)繼續(xù)擴大面積，更有助于發(fā)揮土地規(guī)模效應(yīng)，提高飼草生產(chǎn)技術(shù)效率，釋放單產(chǎn)潛力；種植大戶可以繼續(xù)擴大飼草種植規(guī)模，但種植規(guī)模擴大到0.133~0.266 hm2區(qū)間，會面臨技術(shù)效率降低的風(fēng)險。

另外，從大、中和小不同種植規(guī)模戶均生產(chǎn)狀況來看，小規(guī)模農(nóng)戶技術(shù)效率較高是由于施用農(nóng)家肥和化肥在3種規(guī)模農(nóng)戶中最高。大規(guī)模農(nóng)戶技術(shù)效率高是由于燃料費用投入多，即機械化水平最高。而中等規(guī)模農(nóng)戶的相對低效率是由于相對農(nóng)機和肥料更傾向于使用勞動，勞動時間投入較其他規(guī)模農(nóng)戶多6%，但農(nóng)業(yè)現(xiàn)代科技水平較低。

2.3.3不同地區(qū)農(nóng)戶技術(shù)效率分析從各地區(qū)來看，日喀則地區(qū)飼草生產(chǎn)效率在西藏各地區(qū)中最高，山南地區(qū)次之，拉薩最低。日喀則處于首位的原因是選種了適宜優(yōu)良品種，96.6%農(nóng)戶種植燕麥。近年來日喀則地區(qū)結(jié)合項目開展，以培育和推廣高產(chǎn)且便于管理的禾本科燕麥為主，而山南和拉薩地區(qū)則以推廣牲畜適口性好但產(chǎn)量相對低的豆科紫花苜蓿和箭筈豌豆為主(表3)。

山南生產(chǎn)效率高于拉薩，則由勞動和資金投入差異造成。受農(nóng)業(yè)規(guī)模和非農(nóng)收入間接因素影響，拉薩農(nóng)戶平均耕地面積是山南的2倍多，且拉薩非農(nóng)收入占家庭收入比重較山南高出5個百分點。拉薩相對便利的交通環(huán)境，充足的非農(nóng)就業(yè)機會，繁榮的農(nóng)產(chǎn)品交易市場，使農(nóng)戶資金和勞動力資源的機會成本較其他地區(qū)高。對拉薩農(nóng)戶來說，非農(nóng)就業(yè)、糧食種植和牧業(yè)生產(chǎn)的資源投入動力要明顯大于飼草生產(chǎn)，這就導(dǎo)致拉薩地區(qū)飼草生產(chǎn)勞動、資金投入嚴(yán)重不足，生產(chǎn)效率低下。

表3　2011年各地區(qū)農(nóng)戶飼草生產(chǎn)技術(shù)效率差異

來源：根據(jù)估計結(jié)果整理。

Source: Based on the estimation results.

同品種不同地區(qū)飼草生產(chǎn)效率差異也主要由勞動和資金投入差異造成。燕麥生產(chǎn)日喀則農(nóng)戶技術(shù)效率全區(qū)最高；箭筈豌豆和紫花苜蓿生產(chǎn)山南效率最高(表3)。需要指出的是，拉薩紫花苜蓿生產(chǎn)相對山南的低效率是由于非農(nóng)就業(yè)、糧食種植與飼草生產(chǎn)形成資源競爭而導(dǎo)致勞動力和資金投入不足造成的，而拉薩箭筈豌豆生產(chǎn)相對日喀則的低效率則由于牧業(yè)生產(chǎn)、糧食種植與飼草生產(chǎn)形成資源競爭而勞動力和資金投入不足造成。與拉薩相比，日喀則飼草種植戶非農(nóng)就業(yè)機會成本更高。

3結(jié)論

本文采用隨機前沿生產(chǎn)函數(shù)法，以2011年西藏“一江兩河”地區(qū)農(nóng)戶調(diào)研數(shù)據(jù)為基礎(chǔ)，測算了該區(qū)農(nóng)戶飼草生產(chǎn)技術(shù)效率及估算出效率影響因素，并在飼草科類、農(nóng)戶規(guī)模和地區(qū)分布3個維度，探究飼草生產(chǎn)技術(shù)效率變化規(guī)律，得到如下結(jié)論：

第一，農(nóng)戶飼草生產(chǎn)存在顯著的效率損失，提高生產(chǎn)技術(shù)效率，單產(chǎn)還能提高近30%。目前，飼草生產(chǎn)在西藏處于技術(shù)密集型發(fā)展階段，較勞動、資金和土地投入，肥料、種子等技術(shù)要素投入增產(chǎn)更顯著。需要指出的是，農(nóng)家肥投入對飼草增產(chǎn)效果比化肥好，這對西藏開展有機農(nóng)業(yè)和生態(tài)農(nóng)業(yè)提供了良好的理論支持。農(nóng)戶禾本科飼草種植技術(shù)相似度較高且處于較高水平，而豆科飼草技術(shù)效率分布分散，農(nóng)戶飼草生產(chǎn)技術(shù)差異很大。

第二，依靠改善人力資本狀況和提高農(nóng)戶農(nóng)業(yè)生產(chǎn)能力，對提高飼草生產(chǎn)技術(shù)效率作用不明顯。飼草科類差異會對飼草技術(shù)效率產(chǎn)生直接影響，禾本科飼草種植效率高于豆科飼草，前者會使單產(chǎn)較后者增加2.36%；農(nóng)戶資金狀況和農(nóng)業(yè)規(guī)模通過作用于勞動力、資金等投入要素對生產(chǎn)技術(shù)效率產(chǎn)生間接負(fù)向影響。另外，研究發(fā)現(xiàn)農(nóng)業(yè)規(guī)模與技術(shù)效率呈似“N”型關(guān)系，絕大多數(shù)農(nóng)戶(大于88.42%)飼草種植規(guī)模偏低，如果能在0.133 hm2內(nèi)擴大種植規(guī)模，能更好地發(fā)揮規(guī)模效應(yīng)。農(nóng)戶農(nóng)業(yè)綜合生產(chǎn)能力與飼草生產(chǎn)技術(shù)效率高低無關(guān)。

第三，西藏飼草生產(chǎn)效率地區(qū)差異顯著，日喀則地區(qū)最高，山南地區(qū)次之，拉薩地區(qū)最低。地區(qū)效率差異的主要原因是品種差異；而同品種不同地區(qū)的效率差異則是由勞動和資金要素投入差異引起。燕麥種植日喀則地區(qū)技術(shù)效率最高，紫花苜蓿和箭筈豌豆種植山南地區(qū)技術(shù)效率最高。

4政策建議

西藏飼草生產(chǎn)突破技術(shù)效率瓶頸的首要條件是破除資源約束，其次要因地制宜。政府部門要意識到農(nóng)戶飼草生產(chǎn)的積極性不高，這需要政府從源頭抓起，一靠市場引導(dǎo)，培育草產(chǎn)品市場，提高牧民飼草商品意識；二靠政府合理的資金支持和技術(shù)指導(dǎo)。完善肥料、機械等生產(chǎn)資料購置補貼策略，使之更趨合理和實現(xiàn)常態(tài)化；拓寬農(nóng)業(yè)貸款渠道，為飼草生產(chǎn)提供專項貸款；建立飼草生產(chǎn)培訓(xùn)機制，提高農(nóng)民飼草(主要是豆科飼草)生產(chǎn)技術(shù)。所謂的因地制宜就是指政府要根據(jù)各地區(qū)資源狀況，在各地區(qū)推廣宜種、效率高、性能好的飼草品種。在這方面，日喀則和山南地區(qū)提供了較好的經(jīng)驗，但拉薩地區(qū)由于定位不清，飼草生產(chǎn)處于探索階段，這需要政府在拉薩地區(qū)加大科研投入，找出拉薩地區(qū)適宜種植的飼草品種。另外，由于農(nóng)民飼草種植規(guī)模普遍偏低，可以在確保糧食安全情況下，鼓勵農(nóng)民適當(dāng)擴大播種面積，能更有效地發(fā)揮土地規(guī)模效應(yīng)。

References：

[1]The Whole Country Grass Product Quality Supervision and Inspection Center of Agriculture Ministry. It’s time to pay attention to grass industry[N]. Farmer’s Daily, 10th, Aug, 2011(eight).

[2]Ma Q Z. Protect grassland ecology and contribute beautiful pastoral area. Acta Prataculturae Sinica, 2014, 23(1): 1-2.

[3]Zhang M Y, Xue S M, Cui L Y,etal. Effect of harvest time and harvest method on seed yield and quality ofBrachiariadecumbens‘basilisk’. Acta Prataculturae Sinica, 2014, 23(4): 351-356.

[4]Xin S P, Jia Y D P. Development measures of artificial grassland of Tibet river valley region. China Animal Husbandry, 2012, 6: 32-36.

[5]Xie H, Lu J J. Thoughts on structure adjustment of Tibet agricultural and animal husbandry industries. Journal of Tibet University, 2005, 20(3): 29-32.

[6]Zhang Y F. The role of alfalfa in the structure adjustment of three-component cropping system in Chinese agriculture. Pratacultural Science, 1999,4: 10-13.

[7]Yue Y H, Qi X, Wang Y R,etal. Persistence of 35Medicagosativavarieties at the 10th year after establishment. Acta Prataculturae Sinica, 2014, 23(1): 58-64.

[8]Farrell M J. The measurement of productive efficiency. Journal of the Royal Statistical Society, Series A (General), 1957, 120(2): 253-290.

[9]Leibenstein H. Allocative efficiency vs. “x-efficiency”. The American Economic Review, 1966,56: 392-415.

[10]RolfFare S G, Mary N, Zhang Z Y. Productivity growth technical progress and efficiency changes in industrialised countries. The American Economic Review, 1994, 84: 1033-1041.

[11]Meeusen W, Broeck J. Efficiency estimation from Cobb-Douglas production functions with composed error. International Economic Review, 1977, 18: 435-444.

[12]Lovell A, Schmidt. figuretion and estimation of stochastic frontier production function. Journal of Econometrics, 1977, 6: 21-37.

[13]Battese G E, Coelli T J. A stochastic frontier production function incorporating, a model for technical inefficiency effects[A]. Working Papers in Econometrics and Applied Statistics, No. 69, Department of Econometrics, University of New England, 1993.

[14]Wang G H. The estimation of stochastic frontier production function and average technical efficiency. Journal of Agrotechnical Economics, 1988, 1: 21-25.

[15]Wang G H. The theory and application of frontier production function. Journal of Agrotechnical Economics, 1988, 1: 21-25.

[16]Liao Y X. The discussion and survey of planting industry production technical efficiency. Technology Economics, 1994, 9: 38-42.

[17]Su J C. Estimation and measures chosen of grain production technical efficiency of Guangdong Province. China Rural Survey, 1996, 6: 14-17.

[18]Meng L J. The dynamic study on Chinese agricultural production technical efficiency. Journal of Agrotechnical Economics, 2000, 5: 1-4.

[19]Liang X M. Technical efficiency and economy efficiency[N]. China Economic Times, 23th, Feb, 2004, 5.

[20]Du W J. Policy differences of agricultural production technical efficiency, based on panal data stochastic frontier method. Journal of Quantitative & Technical Economics, 2009, 9: 107-119.

[21]Fang H. Study on China agricultural production technical efficiency, based on province level estimation, found and explanation. Journal of Agrotechnical Economics, 2010, 1: 34-42.

[22]Yu J B, Qiao J, Gong C G. The analysis of Chinese soybean production technical efficiency. Journal of Agrotechnical Economics, 2007, 4: 41-48.

[23]Kang X, Li X M. The analysis of our country’ grain production technical efficiency, based on stochastic frontier method. China Rural Survey, 2005, 4: 25-33.

[24]Meng L J, Zhang H M. Regional technical efficiency differences in wheat production in China. Journal of Nanjing Agricultural University (Social Sciences edition), 2004, 4(2): 13-17.

[25]Cao J, Sun D Q, Tan X Y. Farmer’ cow production technical efficiency and effect factors. China Rural Economy, 2005,10: 42-49.

[26]Mang H Y, Tang H C. The total factors productivity analysis of Chinese milk production. Chinese Rural Economy, 2007, 2: 40-48.

[27]Qu X B. The production technical efficiency difference of different farmer sizes and effect factors. Journal of Nanjing Agricultural University (Social Sciences edition), 2009, 9(3): 27-36.

[28]Jin F L, Wang L, Li G C. Winter rapeseed’s technical efficiency and its influence factors: Based on the model of stochastic frontier production function and 1707 micro-datas of farmers. Journal of China Agricultural University, 2013, 18(1): 210-217.

[29]Battese G E, CoelliT J. A model for technical inefficiency effects in a stochastic frontier production function for panel data. Empirical Economics, 1995, 20: 325-332.

[30]Kumbhakar S C, Ghosh S, McGuckin J T. A generalized production frontier approach for estimating determinants of inefficiency in U. S. dairy farms. Journal of Business & Economics Statistics, 1991, 9: 279-286.

[31]Lucas R EJr. On the mechanics of economic development. Journal of Monetary Economics, 1988, 22(3): 3-42.

[32]Schultz T W.Transforming Traditional Agriculture[M]. Liang X M, Translated. Beijing: Beijing Commerical Press, 1987: 133.

[33]Huang C J, Bagi F S. Technical efficiency on individual farms in Northwest India. Southern Economic Journal, 1984, 51: 108-115.

[34]Coelli T J. A guide to frontier version 4.1: a computer program for stochastic frontier production and cost function estimation[A]. Working Paper, Centre for Efficiency and Productivity Analysis, University of New England, Armidale, Australia, 1996.

參考文獻：

[1]農(nóng)業(yè)部全國草業(yè)產(chǎn)品質(zhì)量監(jiān)督檢驗測試中心. 到該重視草業(yè)的時候了[N]. 農(nóng)民日報, 2011-08-10, 008.

[2]馬啟智. 保護草原生態(tài)、建設(shè)美麗牧區(qū). 草業(yè)學(xué)報, 2014, 23(1): 1-2.

[3]張美艷, 薛世明, 崔玲艷, 等. 不同收種時間和收種方式對貝斯莉斯克伏生臂形草種子產(chǎn)量和質(zhì)量的影響. 草業(yè)學(xué)報, 2014, 23(4): 351-356.

[4]辛盛鵬, 加央旦培. 西藏河谷地帶人工種草發(fā)展與對策. 中國畜牧業(yè), 2012, 6: 32-36.

[5]謝紅, 盧金金. 關(guān)于西藏農(nóng)牧業(yè)結(jié)構(gòu)調(diào)整和發(fā)展的思考. 西藏大學(xué)學(xué)報, 2005, 20(3): 29-32.

[6]張玉發(fā). 試論苜蓿生產(chǎn)在我國農(nóng)業(yè)三元種植結(jié)構(gòu)調(diào)整中的地位和作用. 草業(yè)科學(xué), 1999, 4: 10-13.

[7]岳彥紅, 齊曉, 王彥榮, 等. 35個10齡紫花苜蓿品種的持久性比較. 草業(yè)學(xué)報, 2014, 23(1): 58-64.

[14]萬廣華. 隨機邊界生產(chǎn)函數(shù)和平均技術(shù)效率的估算. 農(nóng)業(yè)技術(shù)經(jīng)濟, 1988, 1: 21-25.

[15]萬廣華. 邊界生產(chǎn)函數(shù)理論及其應(yīng)用. 農(nóng)業(yè)技術(shù)經(jīng)濟, 1989, 2: 25-36.

[16]廖耀先. 種植業(yè)生產(chǎn)技術(shù)效率的討論與度量. 技術(shù)經(jīng)濟, 1994, 9: 38-42.

[17]蘇基才. 廣東省糧食生產(chǎn)技術(shù)效率的測定分析與政策取向. 中國農(nóng)村觀察, 1996, 6: 14-17.

[18]孟令杰. 中國農(nóng)業(yè)產(chǎn)出技術(shù)效率動態(tài)研究. 農(nóng)業(yè)技術(shù)經(jīng)濟, 2000, 5: 1-4.

[19]梁小民. 技術(shù)效率與經(jīng)濟效率[N]. 中國經(jīng)濟時報, 2004-02-23,05.

[20]杜文杰. 農(nóng)業(yè)生產(chǎn)技術(shù)效率的政策差異研究——基于時不變閾值面板隨機前沿分析. 數(shù)量經(jīng)濟技術(shù)經(jīng)濟研究,2009,9: 107-119.

[21]方鴻. 中國農(nóng)業(yè)生產(chǎn)技術(shù)效率研究：基于省級層面的測度、發(fā)現(xiàn)與解釋. 農(nóng)業(yè)技術(shù)經(jīng)濟, 2010, 1: 34-42.

[22]余建斌, 喬娟, 龔崇高. 中國大豆生產(chǎn)的技術(shù)效率和技術(shù)效率分析. 農(nóng)業(yè)技術(shù)經(jīng)濟, 2007, 4: 41-48.

[23]亢霞, 李秀梅. 我國糧食生產(chǎn)的技術(shù)效率分析——基于隨機前沿分析方法. 中國農(nóng)村觀察, 2005, 4: 25-33.

[24]孟令杰, 張紅梅. 中國小麥生產(chǎn)的技術(shù)效率地區(qū)差異. 南京農(nóng)業(yè)大學(xué)學(xué)報(社會科學(xué)版), 2004, 4(2): 13-17.

[25]曹暕, 孫頂強, 譚向勇. 農(nóng)戶奶牛生產(chǎn)技術(shù)效率及影響因素分析. 中國農(nóng)村經(jīng)濟, 2005, 10: 42-49.

[26]馬恒運, 唐華倉. 中國牛奶生產(chǎn)的全要素生產(chǎn)率分析. 中國農(nóng)村經(jīng)濟, 2007, 2: 40-48.

[27]屈小博. 不同規(guī)模農(nóng)戶生產(chǎn)技術(shù)效率差異及其影響因素分析. 南京農(nóng)業(yè)大學(xué)學(xué)報(社會科學(xué)版), 2009, 9(3): 27-36.

[28]金福良, 王璐, 李谷成, 等. 不同規(guī)模農(nóng)戶冬油菜生產(chǎn)技術(shù)效率及影響因素分析——基于隨機前沿函數(shù)與1707個農(nóng)戶微觀數(shù)據(jù). 中國農(nóng)業(yè)大學(xué)學(xué)報, 2013, 18(1): 210-217.

[32]西奧多·W·舒爾茨. 改造傳統(tǒng)農(nóng)業(yè)[M]. 梁小民，譯. 北京：北京商務(wù)印書館, 1987: 133.

http://cyxb.lzu.edu.cn

曲云鶴，余成群，孫維，武俊喜，李少偉. 基于隨機前沿生產(chǎn)函數(shù)法的飼草生產(chǎn)技術(shù)效率測算及影響因素分析——以西藏“一江兩河”地區(qū)為例. 草業(yè)學(xué)報, 2015, 24(10): 70-79.

QU Yun-He, YU Cheng-Qun, SUN Wei, WU Jun-Xi, LI Shao-Wei. Use of the stochastic frontier production function in a technical efficiency analysis of grass production: An example from the “One River and Two Tributaries” region of Tibet. Acta Prataculturae Sinica, 2015, 24(10): 70-79.

通訊作者*Corresponding author. E-mail:quyh@igsnrr.ac.cn

作者簡介：曲云鶴(1979-)，女，黑龍江佳木斯人，在站博士后。E-mail: quyh@igsnrr.ac.cn

基金項目：國家科技支撐課題“藏北退化草地綜合整治技術(shù)與示范”(2011BAC09B03)和國家科技支撐課題“西北和青藏地區(qū)優(yōu)質(zhì)牧草豐產(chǎn)栽培及草畜耦合技術(shù)集成與產(chǎn)業(yè)化示范”(2011BAD17B05)資助。

收稿日期：2014-08-25；改回日期：2015-04-24

DOI：10.11686/cyxb2014360