徐 秋，雷國(guó)平，楊厚翔

(東北大學(xué)上地管理研究所，遼寧沈陽(yáng) 110189)

0 引言

糧食安全關(guān)系國(guó)家長(zhǎng)治久安，耕地是糧食生產(chǎn)的第一資源。但是隨著城鎮(zhèn)化、工業(yè)化進(jìn)程，耕地非農(nóng)化、粗放利用等現(xiàn)象日益加劇[1-2]，耕地資源正在被過(guò)度消耗和占用。為此我國(guó)先后出臺(tái)了一系列有關(guān)耕地?cái)?shù)量保護(hù)、耕地質(zhì)量改造與建設(shè)等方面的政策文件，形成政策高壓態(tài)勢(shì)，以保護(hù)耕地資源。在糧食安全必要物質(zhì)基礎(chǔ)得以保障的前提下，如何高效利用耕地資源，成為現(xiàn)階段提高糧食安全系數(shù)、維護(hù)國(guó)家經(jīng)濟(jì)社會(huì)健康穩(wěn)定運(yùn)行的重要途徑。

目前國(guó)內(nèi)學(xué)者多圍繞效率測(cè)算[3-5]、區(qū)域差異[6-7]、影響因素[8-10]對(duì)耕地利用效率展開有益探索。譚榮等[11]、陳江龍等[12]的早期研究多基于單要素評(píng)估指標(biāo)對(duì)耕地利用效率進(jìn)行測(cè)算，缺乏對(duì)耕地產(chǎn)出的考量，實(shí)則為耕地生產(chǎn)率的研究。隨后，為修正早期研究評(píng)價(jià)方法的不足，葉浩等[13]、劉玉海等[14]引入隨機(jī)前沿函數(shù)方法、DEA方法，基于全要素生產(chǎn)框架下對(duì)全國(guó)耕地利用效率進(jìn)行研究，耕地利用效率研究漸趨成熟，但仍存在著不足，主要體現(xiàn)在缺乏對(duì)環(huán)境負(fù)效應(yīng)的思量。鑒于化肥污染對(duì)糧食生產(chǎn)具有顯著的“環(huán)境懲罰”效應(yīng)[15]，文章選取化肥污染排放量作為非期望產(chǎn)出，運(yùn)用非期望產(chǎn)出SBM(SBM-Undesirable)模型對(duì)1994～2014年黑龍江省及所轄13個(gè)地級(jí)市的耕地利用效率時(shí)空差異及影響因素開展研究，為提高區(qū)域耕地資源利用效率提供科學(xué)依據(jù)。

1 研究區(qū)概況

黑龍江省位于中國(guó)東北部，介于北緯43°26′～53°33′，東經(jīng)121°11′～135°05′之間，是中國(guó)位置最北、緯度最高的省份。新中國(guó)成立以來(lái)累計(jì)為國(guó)家提供商品糧超過(guò)4 500億kg，是國(guó)家重要的商品糧食基地，被譽(yù)為我國(guó)的“戰(zhàn)略大糧倉(cāng)”。全省耕地面積1 585.9萬(wàn)hm2，占土地總面積33.52%，人均耕地面積為全國(guó)水平4.20倍， 78.3%的耕地分布平原地區(qū)，耕地資源豐富。建國(guó)以后先后對(duì)“五荒”開展了數(shù)次大規(guī)模墾殖，導(dǎo)致后備耕地資源數(shù)量逐漸減少，而且有不少被劃為禁止開發(fā)的沼澤地、灘涂。加之建設(shè)占用、災(zāi)毀、生態(tài)退耕、農(nóng)業(yè)結(jié)構(gòu)調(diào)整等原因?qū)е赂財(cái)?shù)量減少、質(zhì)量下降，不利于糧食產(chǎn)量的進(jìn)一步提升。因此，對(duì)黑龍江省耕地利用效率時(shí)空差異及其影響因素進(jìn)行研究，對(duì)提高耕地利用效率、緩解耕地供需矛盾、提高糧食安全系數(shù)具有重要意義。

2 研究方法

2.1 非期望產(chǎn)出SBM模型

(1)

(2)

2.2 變異系數(shù)

運(yùn)用變異系數(shù)反映黑龍江省耕地利用效率的變動(dòng)差異程度。

(3)

式(3)中,v表示耕地利用效率變異系數(shù)；di表示黑龍江省各市的耕地利用效率樣本值；d為各市耕地利用效率均值；k為樣本的城市數(shù)，v值越大，各城市間耕地利用效率差異越大，均衡性越差,反之則越好。

2.3 Tobit回歸分析

Tobit模型是解決因變量受到限制的一種截?cái)嗍交貧w模型。由于采用SBM模型測(cè)度的效率值介于0～1之間，對(duì)這種切割或截?cái)嗟囊蜃兞?，若采用普通最小二乘?OLS)進(jìn)行回歸系數(shù)分析，則估計(jì)結(jié)果會(huì)存在有偏且不一致的情況[18]。因此，該文選擇基于最大似然法概念的Tobit回歸模型。

(4)

式(4)中,yi、xi分別代表因變量、自變量;β為回歸系數(shù)向量;ei服從正太分布N(0，δ2)。

3 參評(píng)指標(biāo)選取

3.1 耕地利用效率內(nèi)涵

效率包括技術(shù)效率和配置效率[19]。技術(shù)效率是指在產(chǎn)出既定的情況下，實(shí)現(xiàn)最小投入的能力，或者是在投入既定的情況下，擴(kuò)大產(chǎn)出的能力； 配置效率是指在一定的要素價(jià)格條件下以最優(yōu)化比例對(duì)投入組合加以利用的能力。在完全競(jìng)爭(zhēng)市場(chǎng)下，配制有效率意為各要素的產(chǎn)出彈性等于投入要素在總成本中的占比。但在實(shí)際經(jīng)濟(jì)中，多數(shù)產(chǎn)業(yè)處于非完全競(jìng)爭(zhēng)市場(chǎng)，因此，國(guó)內(nèi)外學(xué)者們對(duì)效率的測(cè)評(píng)多集中于技術(shù)效率。依此研究思路，該文所研究的耕地利用效率具體是指在既定農(nóng)業(yè)產(chǎn)出下所能實(shí)現(xiàn)最少耕地資源的投入程度。

3.2 評(píng)價(jià)指標(biāo)選取與處理

通過(guò)耕地利用效率內(nèi)涵可以看出，耕地利用效率測(cè)算包括投入和產(chǎn)出兩方面，因此研究從投入、產(chǎn)出兩方面構(gòu)建耕地利用效率指標(biāo)體系，產(chǎn)出指標(biāo)分為期望產(chǎn)出與非期望產(chǎn)出指標(biāo)(表1)。參考已有文獻(xiàn)，投入指標(biāo)從土地、資本、勞動(dòng)力生產(chǎn)三要素方面選取，土地要素選取農(nóng)作物播種面積，反映農(nóng)戶耕地實(shí)際利用情況； 資本要素選取化肥施用量，反映資本投入水平； 勞動(dòng)力要素選取農(nóng)業(yè)機(jī)械總動(dòng)力、農(nóng)業(yè)勞動(dòng)力，反映耕地勞作過(guò)程中投入的勞動(dòng)力、機(jī)械替代勞動(dòng)力水平； 由于研究為狹義范疇上的種植業(yè)，因此期望產(chǎn)出選取種植業(yè)總產(chǎn)值，未選取農(nóng)業(yè)總產(chǎn)值； 非期望產(chǎn)出選取化肥(TN、TP)排放量。

表1 黑龍江省耕地利用效率指標(biāo)體系

項(xiàng)目投入產(chǎn)出期望產(chǎn)出非期望產(chǎn)出變量I1I2I2I4O1O2土地勞動(dòng)力農(nóng)機(jī)化肥農(nóng)業(yè)總產(chǎn)值化肥污染指標(biāo)解釋農(nóng)作物播種面積種植業(yè)勞動(dòng)力農(nóng)用機(jī)械總動(dòng)力化肥折純量種植業(yè)總產(chǎn)值化肥(TN、TP)排放量

該文選擇單元調(diào)查評(píng)估方法測(cè)算耕地利用過(guò)程中的化肥(TN、TP)排放量，計(jì)算公式為：化肥(TN、TP)排放量=氮肥、磷肥×TN/TP產(chǎn)污系數(shù)×流失率，其中氮肥、磷肥各產(chǎn)污系數(shù)、流失系數(shù)參考封永剛[20]、潘丹[21]等研究并結(jié)合《污染源普查農(nóng)業(yè)源系數(shù)手冊(cè)》獲取。

該文運(yùn)用基礎(chǔ)數(shù)據(jù)為1994～2014年黑龍江省13個(gè)地級(jí)市的耕地利用投入產(chǎn)出及影響因素?cái)?shù)據(jù)，包括農(nóng)作物播種面積、種植業(yè)勞動(dòng)力、農(nóng)機(jī)總動(dòng)力、化肥施用量、種植業(yè)總產(chǎn)值，數(shù)據(jù)來(lái)源于《黑龍江省統(tǒng)計(jì)年鑒(1995～2015)》，為規(guī)避價(jià)格因素影響，該文將種植業(yè)總產(chǎn)值統(tǒng)一折算為以1994年為基期的實(shí)際產(chǎn)值。

圖1 1994～2014黑龍江耕地利用效率變化趨勢(shì)

圖2 耕地綜合技術(shù)效率與純技術(shù)效率、 規(guī)模效率的關(guān)系散點(diǎn)分布

圖3 1994～2014年省耕地利用效率變異系數(shù)

4 結(jié)果與分析

4.1 全省耕地利用效率分析

利用MaxDea Pro 1.0對(duì)數(shù)據(jù)進(jìn)行處理得到1994～2014年全省綜合技術(shù)效率、純技術(shù)效率、規(guī)模效率。如圖1所示。

1994～2014年期間，綜合技術(shù)效率平均值為0.679，從1994年的0.835下降到2014年的0.442，降幅為47.04%，可以看出耕地利用效率不高，總體呈下降態(tài)勢(shì)。

為探析耕地利用效率變化原因，將綜合技術(shù)效率分解為純技術(shù)效率和規(guī)模效率即綜合技術(shù)效率(TE)=純技術(shù)效率(PTE)×規(guī)模效率(SE)。其中純技術(shù)效率反映各決策單元在耕地利用過(guò)程中有效利用資源的能力及管理水平的高低，規(guī)模效率反映各決策單元在耕地利用過(guò)程中的生產(chǎn)規(guī)模是否達(dá)到最優(yōu)狀態(tài)。1994～2014年純技術(shù)效率平均值為0.760，規(guī)模效率均值為0.871，純技術(shù)效率、規(guī)模效率與綜合技術(shù)效率變化趨勢(shì)一致，純技術(shù)效率降幅為24.50%，略高于規(guī)模效率20.30%，可以看出純技術(shù)效率始終是制約耕地利用效率的主要方面。為進(jìn)一步說(shuō)明，研究借助SPSS軟件數(shù)據(jù)分析功能，分析綜合技術(shù)效率與分解效率的相關(guān)性并繪制散點(diǎn)分布(圖2)，可以看出純技術(shù)效率與綜合技術(shù)效率的相關(guān)系數(shù)(0.970)明顯高于規(guī)模效率與綜合技術(shù)效率的相關(guān)系數(shù)(0.927)，表明黑龍江省耕地利用效率與純技術(shù)效率關(guān)系更密切。

圖4 區(qū)域耕地利用效率類型區(qū)及投入要素冗余情況

4.2 耕地利用效率區(qū)域差異分析

利用MaxDea Pro 1.0對(duì)數(shù)據(jù)進(jìn)行處理得到1994～2014年各行政區(qū)綜合技術(shù)效率、純技術(shù)效率和規(guī)模效率。在此基礎(chǔ)上利用式(3)計(jì)算黑龍江省耕地利用時(shí)間序列綜合技術(shù)效率、純技術(shù)效率、規(guī)模效率變異系數(shù)(圖3)。

可以看出雖然部分年份區(qū)域間耕地利用效率分異程度有所波動(dòng)，但總體上呈增大趨勢(shì)，表明區(qū)域間耕地利用效率區(qū)域差異越來(lái)越大，其中1994～2005年各區(qū)域耕地利用效率差異較小，但變幅較大， 2006～2014年耕地利用效率區(qū)域差異較大，但變幅較小。進(jìn)一步分析可知1994～2005年間純技術(shù)效率區(qū)域分異變幅較大，綜合技術(shù)效率區(qū)域分異主要受純技術(shù)效率的影響； 2006～2014年間規(guī)模效率區(qū)域分異變幅較大，該時(shí)期綜合技術(shù)效率區(qū)域分異主要受規(guī)模效率影響。

4.3 耕地利用效率分區(qū)分析

1994～2014年間耕地利用綜合技術(shù)效率變異系數(shù)均值為0.370，純技術(shù)效率為0.292，規(guī)模效率為0.125，表明耕地利用效率區(qū)域差異較大。因此為了解全省耕地利用效率區(qū)域分異規(guī)律，研究匯總求出1997～2004年各行政區(qū)平均綜合技術(shù)效率、純技術(shù)效率和規(guī)模效率(表2)，并以此為基礎(chǔ)將全省所轄13個(gè)地市的綜合效率進(jìn)行聚類。

表2 1994～2014年黑龍江省各地市耕地利用效率及變化情況

DMUTEPTESEDMUTEPTESE哈爾濱0.3610.4170.838佳木斯0.5860.6820.833齊齊哈爾0.9000.9690.925七臺(tái)河0.5660.7670.739雞西0.6060.6630.884牡丹江0.9090.9570.944鶴崗0.4540.5390.810黑河0.5620.7060.792雙鴨山0.6290.7260.855綏化0.7930.8470.919大慶0.5330.6230.839大興安嶺1.0001.0001.000伊春0.9260.9840.939平均值0.6790.7600.871

經(jīng)聚類全省耕地利用效率可分為3類(圖4)，其中I類地區(qū)包括齊齊哈爾、伊春、牡丹江、大興安嶺； Ⅱ類地區(qū)，包括雞西、雙鴨山、大慶、佳木斯、七臺(tái)河、黑河、綏化； Ⅲ類地區(qū)包括哈爾濱、鶴崗。

4.3.1 耕地利用效率分區(qū)合理性分析

基于變異系數(shù)法，利用式(3)計(jì)算3類地區(qū)耕地利用綜合技術(shù)效率、純技術(shù)效率和規(guī)模效率的變異系數(shù)，以此表征3類地區(qū)耕地利用效率分異情況。經(jīng)測(cè)算3類地區(qū)耕地利用綜合技術(shù)效率變異分別為I類地區(qū)(0.101)、Ⅱ類地區(qū)(0.208)、Ⅲ類地區(qū)(0.198)，小于全省的0.370； 純技術(shù)效率變異系數(shù)分別為I類地區(qū)(0.063)、Ⅱ類地區(qū)(0.197)、Ⅲ類地區(qū)(0.245)，小于全省的0.292； 規(guī)模效率變異系數(shù)分別為I類地區(qū)(0.104)、Ⅱ類地區(qū)(0.103)、Ⅲ類地區(qū)(0.065)，小于全省的0.125。可以看出3類地區(qū)之間變異程度較大，3類地區(qū)內(nèi)部耕地利用效率相對(duì)一致，表明聚類分區(qū)合理，有效地將耕地利用效率相似的區(qū)域劃分為同一類型區(qū)，與其他類型區(qū)相異。

4.3.2 不同類型區(qū)耕地利用效率差異分析

I類地區(qū)，綜合技術(shù)效率值高于0.90，高于全省平均水平，為耕地利用效率高值區(qū)，該類型區(qū)純技術(shù)效率、規(guī)模效率值均在0.90之上，純技術(shù)效率略高于規(guī)模效率且接近生產(chǎn)前沿面。Ⅱ類地區(qū)整體耕地綜合技術(shù)效率水平在0.50～0.80之間，低于I類地區(qū)，與全省平均水平相當(dāng)，處于相對(duì)中等效率區(qū)，該類型區(qū)純技術(shù)效率在0.62～0.85之間，規(guī)模效率在0.73～0.92之間，規(guī)模效率高于純技術(shù)效率。Ⅲ類地區(qū)的耕地綜合技術(shù)效率值低于0.50，低于全省平均水平，為耕地利用效率低值區(qū)，該類型區(qū)純技術(shù)效率值為0.41～0.54，規(guī)模效率值為0.80～0.84，規(guī)模效率高于純技術(shù)效率。

3類地區(qū)相較，I類地區(qū)農(nóng)業(yè)生產(chǎn)要素資源投入規(guī)模達(dá)到較高水平并實(shí)行了合理利用，但仍有進(jìn)一步提升空間，兩者相較應(yīng)側(cè)重規(guī)模效率的提升。Ⅱ地區(qū)農(nóng)業(yè)生產(chǎn)要素資源投入規(guī)模、合理利用水平仍有不足，純技術(shù)效率與規(guī)模效率應(yīng)同步優(yōu)化。Ⅲ類地區(qū)純技術(shù)效率遠(yuǎn)低于規(guī)模效率，因此應(yīng)以純技術(shù)效率提升為主。

4.3.3 不同類型區(qū)投入要素分析

如圖4所示，相較于最優(yōu)投入，3類地區(qū)各投入要素均存在不同程度的冗余，I類地區(qū)各投入要素冗余相對(duì)較低，冗余率在1.80%～18.29%范圍內(nèi)波動(dòng)； Ⅱ類地區(qū)次之，各要素冗余率在10.00%～60.90%之間變化，相較I類地區(qū)，勞動(dòng)力、化肥冗余呈大幅增長(zhǎng)，均高于32.75%； Ⅲ類地區(qū)各投入要素冗余均高于其他兩類地區(qū)，冗余率在35.80%～74.42%之間波動(dòng)，表明該區(qū)實(shí)際投入與最優(yōu)投入差距較大，存在大量資源浪費(fèi)。其中勞動(dòng)力、化肥冗余相對(duì)明顯，均高于62.01%。因此在日后耕地利用過(guò)程中要注重農(nóng)村剩余勞動(dòng)力轉(zhuǎn)移及農(nóng)業(yè)生態(tài)可持續(xù)發(fā)展問(wèn)題。

4.4 耕地利用效率的影響因素與對(duì)策建議

表3 黑龍江省耕地利用效率Tobit回歸分析

變量系數(shù)標(biāo)準(zhǔn)誤差Z值P值常數(shù)項(xiàng)0.9924 0.15556.38040.0000農(nóng)村勞動(dòng)力人均播種面積0.0857***0.02173.94190.0001單位面積農(nóng)業(yè)機(jī)械總動(dòng)力-0.1593***0.0496-3.20990.0013單位面積化肥施用用量-0.0322***0.0085-3.78450.0002勞動(dòng)力比例-0.2834**0.1103-2.56870.0102二、三產(chǎn)業(yè)比重0.12210.18450.66160.5083人均種植業(yè)總收入0.19120.26820.71290.4759 注:***、**表示變量值分別在1%、5%水平下顯著

該文參照已有研究，耕地利用效率受耕地資源稟賦、資本投入條件、經(jīng)濟(jì)發(fā)展水平等方面的影響[18，22-23]。其中耕地資源稟賦方面選取農(nóng)村勞動(dòng)力人均播種面積； 資本投入條件方面選取單位播種面積農(nóng)業(yè)機(jī)械總動(dòng)力、單位播種面積化肥施用量、農(nóng)業(yè)勞動(dòng)力比例； 經(jīng)濟(jì)發(fā)展水平方面選取二、三產(chǎn)業(yè)比重與人均種植業(yè)總收入。研究以1994～2014年各市的耕地利用效率作為因變量，運(yùn)用Eview9.0對(duì)21年的面板數(shù)據(jù)進(jìn)行Tobit回歸分析，P<0.05表明該要素對(duì)耕地利用效率影響顯著； 回歸系數(shù)表明影響程度及影響方向，系數(shù)為正表明該要素與耕地利用效率成正向影響，進(jìn)一步加大該投入可促進(jìn)耕地利用效率提升，系數(shù)為負(fù)表明該要素與耕地利用效率成負(fù)向影響，即表明該要素投入冗余，應(yīng)減少該要素的投入，結(jié)果詳見表3。

由表3可知，除二、三產(chǎn)業(yè)比重、人均種植業(yè)總收入外，其他變量均通過(guò)顯著性水平檢驗(yàn)，表明以上變量均對(duì)耕地利用效率存在顯著影響，按影響程度劃分依次為勞動(dòng)力比例、單位播種面積農(nóng)業(yè)機(jī)械總動(dòng)力、農(nóng)村勞動(dòng)力人均播種面積、單位播種面積化肥施用用量。

勞動(dòng)力比例是指種植業(yè)勞動(dòng)力占農(nóng)業(yè)勞動(dòng)力的比例，反映勞動(dòng)力要素投入水平，該指標(biāo)與耕地利用效率回歸系數(shù)為負(fù)，表明現(xiàn)階段種植業(yè)勞動(dòng)力比例越高，耕地利用效率越低，說(shuō)明目前黑龍江省種植業(yè)勞動(dòng)力投入存在冗余，因此應(yīng)推進(jìn)農(nóng)村剩余勞動(dòng)力轉(zhuǎn)移。單位面積農(nóng)業(yè)機(jī)械動(dòng)力反映區(qū)域農(nóng)業(yè)機(jī)械投入水平，該指標(biāo)與耕地利用效率回歸系數(shù)為負(fù)，表明農(nóng)業(yè)機(jī)械投入冗余。一方面，黑龍江省由于農(nóng)機(jī)購(gòu)置補(bǔ)貼政策使得農(nóng)業(yè)機(jī)械裝備大量增加，據(jù)統(tǒng)計(jì)年鑒顯示，農(nóng)業(yè)機(jī)械總動(dòng)力從1994年的1 190萬(wàn)kW增長(zhǎng)到2014年5 155.5萬(wàn)kW，但是由于地形地貌等自然原因、家庭聯(lián)產(chǎn)承包責(zé)任制以及高中低產(chǎn)田等級(jí)平均分配機(jī)制等社會(huì)經(jīng)濟(jì)原因，區(qū)域耕地平整度、集中連片程度低，人均經(jīng)營(yíng)面積少且分散，并不利于大機(jī)械化作業(yè)，導(dǎo)致機(jī)械裝備利用效率低； 另一方面，機(jī)械化作業(yè)和規(guī)模經(jīng)營(yíng)地區(qū)改變了精耕細(xì)作的生產(chǎn)方式，使得耕地生產(chǎn)潛力無(wú)法充分發(fā)揮，因此要在統(tǒng)籌考慮區(qū)域耕作田塊條件、土地流轉(zhuǎn)、城市化水平等因素的基礎(chǔ)上，兼顧機(jī)械作業(yè)與精耕細(xì)作，合理適度引入農(nóng)業(yè)機(jī)械裝備，提高耕地利用效率。農(nóng)村勞動(dòng)力人均播種面積反映經(jīng)營(yíng)規(guī)模，該指標(biāo)與耕地利用效率相關(guān)系數(shù)為正，即經(jīng)營(yíng)規(guī)模越大，耕地利用效率越高，勞均經(jīng)營(yíng)規(guī)模小，既表明人均經(jīng)營(yíng)面積不足，又側(cè)面反映出農(nóng)村勞動(dòng)力存在冗余，因此應(yīng)通過(guò)土地流轉(zhuǎn)，農(nóng)村剩余勞動(dòng)力轉(zhuǎn)移來(lái)增加勞均經(jīng)營(yíng)規(guī)模，提高區(qū)域耕地利用效率。單位面積化肥施用量反映區(qū)域化肥投入水平，該指標(biāo)與耕地利用效率回歸系數(shù)為負(fù)，即單位面積化肥施用量增加，耕地利用效率不升反降，表明化肥投入冗余，冗余部分還會(huì)造成土壤理化性質(zhì)惡化、非點(diǎn)源污染等生態(tài)環(huán)境問(wèn)題，因此耕地利用過(guò)程中要科學(xué)合理施肥，提高耕地利用效率，避免因過(guò)量施肥造成生態(tài)外部性問(wèn)題，促進(jìn)農(nóng)業(yè)生態(tài)可持續(xù)發(fā)展。

5 結(jié)論與討論

運(yùn)用非期望產(chǎn)出SBM模型測(cè)度了1994～2014年間黑龍江省以及13個(gè)地級(jí)市耕地利用效率，利用變異系數(shù)法分析了1994～2014年間耕地利用效率區(qū)域變異情況，經(jīng)聚類對(duì)13個(gè)地級(jí)市的耕地利用效率進(jìn)行分區(qū)，在分區(qū)合理性分析基礎(chǔ)上對(duì)不同類型區(qū)耕地利用效率及投入要素是否存在冗余進(jìn)行分析，通過(guò)Tobit模型診斷耕地利用效率影響因素，并依此提出對(duì)應(yīng)建議。

(1)1994～2004年耕地利用綜合技術(shù)效率總體呈現(xiàn)下降趨勢(shì)，效率值從1994年的0.835下降到2014的0.442，純技術(shù)效率不高是制約綜合技術(shù)效率提升的主要原因。因此，在今后的土地管理中更應(yīng)注重提升管理水平、農(nóng)業(yè)科技改進(jìn)，以此提高耕地利用效率。

(2)1994～2014年間黑龍江省耕地利用效率區(qū)域差異總體呈現(xiàn)增大趨勢(shì)； 1994～2005年各區(qū)域耕地利用效率差異較小，但變幅較大，該時(shí)期綜合技術(shù)效率區(qū)域分異主要受純技術(shù)效率的影響； 2006～2014年耕地利用效率區(qū)域差異較大，但變幅較小，該時(shí)期綜合技術(shù)效率分異主要受規(guī)模效益影響，該趨勢(shì)與2004年全國(guó)范圍內(nèi)實(shí)施農(nóng)機(jī)購(gòu)置補(bǔ)貼政策，以及全省實(shí)行免除農(nóng)業(yè)稅相符，可以看出惠農(nóng)政策的出臺(tái)激勵(lì)了農(nóng)業(yè)資本投入，使得農(nóng)業(yè)資本投入規(guī)模逐年增長(zhǎng)。

(3)黑龍江省耕地利用效率區(qū)域差異較大，通過(guò)聚類將其分為3類地區(qū)，并對(duì)聚類結(jié)果進(jìn)行合理性檢驗(yàn)，經(jīng)檢驗(yàn)，聚類分區(qū)科學(xué)有效，有效地將耕地利用效率相似的區(qū)域劃分為同一類型區(qū)，與其他類型區(qū)相異。其中I類地區(qū)包括齊齊哈爾、伊春、牡丹江、大興安嶺，該類型區(qū)為耕地利用綜合技術(shù)效率高效率區(qū)，純技術(shù)效率略高于規(guī)模效率且接近生產(chǎn)前沿面，表明地區(qū)農(nóng)業(yè)生產(chǎn)要素資源投入規(guī)模達(dá)到較高水平并實(shí)行了合理利用，但仍有進(jìn)一步提升空間，兩者相較應(yīng)側(cè)重規(guī)模效率； Ⅱ類地區(qū)，包括雞西、雙鴨山、大慶、佳木斯、七臺(tái)河、黑河、綏化，該類型區(qū)耕地利用綜合技術(shù)效率處于中等效率區(qū)，純技術(shù)效率和規(guī)模效率與全省平均水平相當(dāng)，表明農(nóng)業(yè)生產(chǎn)要素資源投入規(guī)模、合理利用水平仍有不足，純技術(shù)效率與規(guī)模效率應(yīng)同步優(yōu)化； Ⅲ類地區(qū)包括哈爾濱、鶴崗，為耕地利用綜合技術(shù)效率低效率區(qū)，該類型區(qū)規(guī)模效率與全省平均水平相當(dāng)，但純技術(shù)效率遠(yuǎn)低于規(guī)模效率，因此應(yīng)以純技術(shù)效率提升為主。3類地區(qū)各投入要素均存在不同程度的冗余，投入要素利用率I類地區(qū)>Ⅱ類地區(qū)>Ⅲ類地區(qū)。

(4)耕地利用效率影響因素按影響程度依次為勞動(dòng)力比例、單位播種面積農(nóng)業(yè)機(jī)械總動(dòng)力、農(nóng)村勞動(dòng)力人均播種面積、單位播種面積化肥施用量。其中農(nóng)村勞動(dòng)力人均播種面積對(duì)耕地利用效率為正向影響，其余均為負(fù)向影響。因此今后應(yīng)通過(guò)以下措施提高耕地利用效率： 通過(guò)土地流轉(zhuǎn)，增加勞均經(jīng)營(yíng)規(guī)模； 推進(jìn)農(nóng)村剩余勞動(dòng)力轉(zhuǎn)移； 統(tǒng)籌考慮區(qū)域耕作田塊條件、土地流轉(zhuǎn)、城市化水平等因素的基礎(chǔ)上，兼顧機(jī)械作業(yè)與精耕細(xì)作，合理適度引入農(nóng)業(yè)機(jī)械裝備； 科學(xué)合理施肥，提高耕地利用效率，避免因過(guò)量施肥造成生態(tài)外部性問(wèn)題。

(5)耕地利用效率影響因素方面，僅考慮了農(nóng)村勞動(dòng)力人均播種面積、單位面積農(nóng)業(yè)機(jī)械總動(dòng)力、單位面積化肥施用用量、勞動(dòng)力比例、二、三產(chǎn)業(yè)比重、人均種植業(yè)總收入等方面因素，但是耕地利用效率還受到區(qū)域氣候、地形地貌、水文地質(zhì)和基礎(chǔ)設(shè)施水平等方面的影響，將會(huì)在后續(xù)研究中做進(jìn)一步深入探討。

[1] 楊志海，李鵬，王雅鵬.農(nóng)村勞動(dòng)力老齡化對(duì)農(nóng)戶耕地利用效率的影響.地域研究與開發(fā), 2015，34(5)： 167～171

[2] 樂(lè)容潮， 雷國(guó)平，丁雪，等.哈爾濱市耕地集約利用水平預(yù)測(cè)及障礙因子診斷.中國(guó)農(nóng)業(yè)資源與區(qū)劃, 2017，38(1)： 59～66

[3] 張霞，劉秀華，劉勇.基于DEA的高家鎮(zhèn)耕地利用效率及其影響因素研究.西南師范大學(xué)學(xué)報(bào)(自然科學(xué)版), 2012，37(3)： 72～77

[4] 梁流濤， 曲福田，王春華.基于DEA方法的耕地利用效率分析.長(zhǎng)江流域資源與環(huán)境, 2008，(2)： 242～246

[5] 經(jīng)陽(yáng)， 葉長(zhǎng)盛.基于DEA的江西省耕地利用效率及影響因素分析.水土保持研究, 2015，22(1)： 257～261

[6] 單玉紅， 朱楓，柯新利.湖北省耕地利用效率的區(qū)域差異及提升對(duì)策.水土保持通報(bào)， 2016，36(2)： 288～292,299

[7] 宮攀， 韓振鈴.基于DEA模型的山東省耕地投入產(chǎn)出效率研究.中國(guó)農(nóng)業(yè)資源與區(qū)劃, 2015，36(5)： 123～131

[8] 許恒周， 郭玉燕，吳冠岑.農(nóng)民分化對(duì)耕地利用效率的影響——基于農(nóng)戶調(diào)查數(shù)據(jù)的實(shí)證分析.中國(guó)農(nóng)村經(jīng)濟(jì)， 2012，(6)： 31～39,47

[9] 李鑫， 歐名豪，馬賢磊.基于景觀指數(shù)的細(xì)碎化對(duì)耕地利用效率影響研究——以揚(yáng)州市里下河區(qū)域?yàn)槔?自然資源學(xué)報(bào), 2011，26(10)： 1758～1767

[10]楊俊， 楊鋼橋，胡賢輝.農(nóng)業(yè)勞動(dòng)力年齡對(duì)農(nóng)戶耕地利用效率的影響——來(lái)自不同經(jīng)濟(jì)發(fā)展水平地區(qū)的實(shí)證.資源科學(xué), 2011，33(9)： 1691～1698

[11]譚榮， 曲福田.中國(guó)農(nóng)地非農(nóng)化與農(nóng)地資源保護(hù)：從兩難到雙贏.管理世界， 2006，(12)： 50～59,66

[12]陳江龍， 曲福田，陳雯.農(nóng)地非農(nóng)化效率的空間差異及其對(duì)土地利用政策調(diào)整的啟示.管理世界， 2004，(8)： 37～42,155

[13]葉浩， 濮勵(lì)杰.我國(guó)耕地利用效率的區(qū)域差異及其收斂性研究.自然資源學(xué)報(bào), 2011，(9)： 1467～14

[14]劉玉海， 張麗.耕地生產(chǎn)率與全要素耕地利用效率——基于SBM-DEA方法的省際數(shù)據(jù)比較.農(nóng)業(yè)技術(shù)經(jīng)濟(jì), 2012，(6)： 47～56

[15]陸文聰， 劉聰.化肥污染對(duì)糧食作物生產(chǎn)的環(huán)境懲罰效應(yīng).中國(guó)環(huán)境科學(xué), 2017，37(5)： 1988～1994

[16]Tone K.A Slacksbased measure of efficiency in Data Envelopment Analysis.European Journal of Operational Research， 2001, 130(3): 498～509

[17]Cooper W W，Seiford L M，Tone K.Data Envelopment Analysis(Second Edition).Boston：Kluwer Academic Publishers， 2007

[18]楊朔， 李世平，羅列.陜西省耕地利用效率及其影響因素研究.中國(guó)土地科學(xué), 2011，25(2)： 47～54

[19]Farrell M J.The measurement of productive efficiency.Journal of Royal Statistic Society， 1957，Series A， 120(Ⅲ)： 253～281

[20]封永剛， 彭玨，鄧宗兵，等.面源污染、碳排放雙重視角下中國(guó)耕地利用效率的時(shí)空分異.中國(guó)人口·資源與環(huán)境, 2015，25(8)： 18～25

[21]潘丹. 考慮資源環(huán)境因素的中國(guó)農(nóng)業(yè)生產(chǎn)率研究.南京：南京農(nóng)業(yè)大學(xué)， 2012

[22]李在軍， 管衛(wèi)華，臧磊.山東省耕地生產(chǎn)效率及影響因素分析.世界地理研究, 2013，22(2)： 167～175

[23]張浩， 李世平，靳亞亞，等.河北省耕地利用效率空間差異及影響因素.水土保持研究, 2016，23(3)： 185～191