摘 要：以2000-2018年黑龍江省13個(gè)地級(jí)市(區(qū))為樣本，運(yùn)用DEA-Windows模型評(píng)價(jià)黑龍江各地市糧食生產(chǎn)綜合效率，采用ArcGIS分析空間差異的時(shí)空演變，運(yùn)用面板Tobit模型對(duì)影響區(qū)域糧食生產(chǎn)綜合效率的因素進(jìn)行分析研判。結(jié)果表明：黑龍江糧食生產(chǎn)綜合效率呈現(xiàn)波動(dòng)上升的趨勢，各地市之間的差異明顯但呈縮小態(tài)勢;糧食播種面積和人均GDP對(duì)糧食生產(chǎn)綜合效率具有顯著的促進(jìn)作用，化肥折純施用量對(duì)糧食生產(chǎn)綜合效率具有顯著的負(fù)向影響，農(nóng)業(yè)從業(yè)人員和農(nóng)業(yè)機(jī)械總動(dòng)力的作用效果不顯著。合理配置糧食生產(chǎn)資源，提高農(nóng)業(yè)科技的創(chuàng)新能力，完善糧食貿(mào)易和糧食補(bǔ)貼制度，是提升糧食生產(chǎn)綜合效率的主要對(duì)策。

　　關(guān)鍵詞：糧食生產(chǎn)效率 空間差異 影響因素 DEA-Windows 黑龍江省

　　民以食為天，隨著社會(huì)經(jīng)濟(jì)的持續(xù)發(fā)展，糧食的供需結(jié)構(gòu)呈現(xiàn)出不同的發(fā)展態(tài)勢。現(xiàn)有糧食生產(chǎn)方式不能保證糧食產(chǎn)量穩(wěn)定增長，并對(duì)生態(tài)環(huán)境造成破壞。認(rèn)清糧食生產(chǎn)現(xiàn)狀、理清糧食生產(chǎn)效率的影響因素、優(yōu)化資源配置才是促進(jìn)糧食生產(chǎn)能力可持續(xù)發(fā)展的重要舉措。研究黑龍江省13個(gè)地級(jí)市(區(qū))的糧食生產(chǎn)綜合效率，分析影響黑龍江省糧食生產(chǎn)綜合效率的因素對(duì)穩(wěn)定全國糧食儲(chǔ)備和維護(hù)糧食安全至關(guān)重要。

　　目前學(xué)界有關(guān)生產(chǎn)效率測算方法和影響因素方面的研究成果豐富。Farrell[1]最早提出了測量糧食生產(chǎn)效率方法，通過測算農(nóng)業(yè)生產(chǎn)效率得到生產(chǎn)效率前沿面;黃金波[2]等運(yùn)用隨機(jī)前沿生產(chǎn)函數(shù)估算了中國糧食生產(chǎn)的技術(shù)效率，發(fā)現(xiàn)農(nóng)業(yè)基礎(chǔ)設(shè)施建設(shè)和制度因素是影響技術(shù)效率的關(guān)鍵;肖紅波[3]等利用 DEA 模型測算了中國糧食生產(chǎn)的技術(shù)效率，認(rèn)為技術(shù)創(chuàng)新在影響因素中呈減退趨勢，提出應(yīng)加強(qiáng)科技創(chuàng)新，倡導(dǎo)規(guī)?；a(chǎn);曾福生、高鳴[4]通過非徑向SBM模型測算了中國的糧食生產(chǎn)效率，得出糧食播種面積、化肥施用量和有效灌溉面積對(duì)糧食生產(chǎn)效率具有顯著影響;唐建[5]采用隨機(jī)前沿生產(chǎn)函數(shù)估算中國31個(gè)省糧食生產(chǎn)的技術(shù)效率，得出家庭經(jīng)濟(jì)水平和勞動(dòng)力受教育水平與技術(shù)效率呈現(xiàn)顯著正向相關(guān)的結(jié)論;楊勇[6]等人以華北平原縣域?yàn)檠芯繂卧?，建立隨機(jī)前沿生產(chǎn)函數(shù)估算糧食的生產(chǎn)效率，表明提高耕地的復(fù)種指數(shù)可以顯著提升糧食的生產(chǎn)效率;曹慧、趙凱[7]通過超越對(duì)數(shù)生產(chǎn)函數(shù)模型分析全國13個(gè)主產(chǎn)區(qū)糧食生產(chǎn)的技術(shù)效率和時(shí)空演變，發(fā)現(xiàn)糧食的產(chǎn)出與勞動(dòng)力、機(jī)械動(dòng)力、農(nóng)藥等都有顯著的相關(guān)關(guān)系;田紅宇、祝志勇[8]基于DEA-Tobit兩步法研究中國省域糧食生產(chǎn)效率，發(fā)現(xiàn)支農(nóng)政策、化肥、農(nóng)業(yè)機(jī)械等影響因素存在地區(qū)差異;何悅、漆雁斌[9]研究城鎮(zhèn)化如何影響糧食生產(chǎn)的技術(shù)效率，發(fā)現(xiàn)人口城鎮(zhèn)化的影響呈“倒U型”。

　　本文在已有研究的基礎(chǔ)上，基于2000-2018年黑龍江省糧食生產(chǎn)的面板數(shù)據(jù)，通過DEA-Windows模型測算不同地區(qū)糧食生產(chǎn)的綜合效率，分析糧食生產(chǎn)綜合效率的時(shí)空演變，提出提升糧食生產(chǎn)綜合效率的有效途徑，為黑龍江省糧食生產(chǎn)效率的研究提供理論支持。

　　1 材料與方法

　　1.1 研究區(qū)概況

　　黑龍江省位于121°11′E-135°05′E、43°25′N-53°33N′，屬于寒溫帶大陸性季風(fēng)氣候，四季分明，夏季雨熱同期，冬季漫長寒冷。全省年平均氣溫多在-4℃～5℃之間，年降水量400～650 mm。自然資源豐富，土壤有機(jī)質(zhì)含量高，是世界著名的三大黑土帶之一，包括三江平原、松嫩平原，是中國耕地面積最大的省份。作為重要的商品糧基地，盛產(chǎn)大豆、水稻、玉米、小麥、馬鈴薯等糧食作物。糧食的播種面積達(dá)到農(nóng)作物面積的95 %左右，2018年全省糧食的播種面積為21321.8萬畝，糧食產(chǎn)量達(dá)到7507萬噸，位居全國首位。

　　1.2 測算指標(biāo)與模型

　　1.2.1 測算指標(biāo)

　　本文基于糧食產(chǎn)出視角對(duì)黑龍江糧食生產(chǎn)的綜合效率進(jìn)行研究。根據(jù)指標(biāo)選取的可獲取性和代表性等原則，構(gòu)建糧食生產(chǎn)綜合效率的評(píng)價(jià)體系如表1。選取農(nóng)業(yè)從業(yè)人員、糧食播種面積、化肥折純施用量、農(nóng)用機(jī)械總動(dòng)力，分別表征糧食生產(chǎn)過程中勞動(dòng)力、土地、資本方面的投入方面;產(chǎn)出方面則選取地區(qū)糧食作物產(chǎn)量和糧食作物產(chǎn)值來表征。由于統(tǒng)計(jì)年鑒未對(duì)糧食生產(chǎn)中的各項(xiàng)指標(biāo)做專項(xiàng)統(tǒng)計(jì)，需從農(nóng)林牧副漁業(yè)中剝離出來。參考已有研究，將農(nóng)業(yè)從業(yè)人員和糧食作物產(chǎn)值采用權(quán)重A進(jìn)行轉(zhuǎn)化，化肥折純施用量和農(nóng)業(yè)機(jī)械總動(dòng)力使用權(quán)重B進(jìn)行調(diào)整。

　　權(quán)重A=(農(nóng)業(yè)總產(chǎn)值/農(nóng)林牧漁業(yè)總產(chǎn)值)×(糧食播種面積/農(nóng)作物總播種面積) 。

　　權(quán)重B=糧食播種面積/農(nóng)作物總播種面積。

　　1.2.2 DEA-Window模型

　　傳統(tǒng)DEA方法主要基于靜態(tài)數(shù)據(jù)，無法分析黑龍江省各地級(jí)市在不同時(shí)間序列的效率變動(dòng)，處理面板數(shù)據(jù)的時(shí)候存在問題。參考Charnes 等人的DEA-Windows 分析方式[10]，將黑龍江省各地級(jí)市的不同時(shí)段看作不同的研究單元，從而實(shí)現(xiàn)了各單元效率值的動(dòng)態(tài)比較。假設(shè)黑龍江省T時(shí)間序列內(nèi)存在N個(gè)研究單元(DMU)，視窗寬度為D，則第i個(gè)視窗包含i到d+i-1時(shí)間段內(nèi)的N × D個(gè)DMU，其中1≤i≤T-D+1。模型設(shè)定DMUnt(1≤n≤N，1≤t≤T)，視窗P(1≤p≤T-D+1)的投入產(chǎn)出矩陣為：

　　(1)

　　(2)

　　在規(guī)模報(bào)酬不變的情況下，根據(jù)公式(3)計(jì)算得出效率值。

　　(3)

　　1.2.3 面板Tobit模型

　　由于運(yùn)用 DEA-Windows 模型測算的效率值介于 0～1 ，存在明顯的截?cái)鄶?shù)據(jù)特征，用普通最小二乘法進(jìn)行回歸的參數(shù)估計(jì)可能會(huì)造成偏差，因此本文通過 Tobit 模型，研究影響黑龍江省糧食生產(chǎn)綜合效率的因素。

　　模型具體形式如下：

　　(4)

　　其中，n代表黑龍江各地級(jí)市，t代表年份，Y表示糧食生產(chǎn)效率值，X是影響黑龍江糧食生產(chǎn)效率的各解釋變量，μn為表示個(gè)體異質(zhì)性的截距項(xiàng)，εnt是隨機(jī)擾動(dòng)項(xiàng)。

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

　　糧食生產(chǎn)的綜合效率是勞動(dòng)力、耕地資源、原料和技術(shù)投入與所產(chǎn)生的經(jīng)濟(jì)效益之間配置情況的表現(xiàn)，反映糧食生產(chǎn)中各要素整體利用程度。本研究以2000-2018年黑龍江省13個(gè)地級(jí)市(包括大興安嶺地區(qū))為研究對(duì)象，由于統(tǒng)計(jì)口徑不一致未將農(nóng)墾總局包含在內(nèi)，個(gè)別年份缺失數(shù)據(jù)用插值法補(bǔ)齊，研究所需數(shù)據(jù)均來源于《黑龍江省統(tǒng)計(jì)年鑒》(2000-2018年)。

　　2 結(jié)果與分析

　　2.1 黑龍江省糧食生產(chǎn)效率測算與空間差異分析

　　2.1.1 測算結(jié)果

　　基于2000-2018年黑龍江省13個(gè)地級(jí)市的面板數(shù)據(jù)，以產(chǎn)出為導(dǎo)向構(gòu)建規(guī)模報(bào)酬不變的DEA-Windows 模型，設(shè)置窗寬為10[11]，計(jì)算得出黑龍江省的糧食生產(chǎn)綜合效率如表2。

　　從長時(shí)間序列上的變化來看，黑龍江省糧食生產(chǎn)效率分為2個(gè)階段，第1階段是2000-2012年黑龍江省糧食生產(chǎn)效率總體上呈現(xiàn)波動(dòng)上升的趨勢。在這期間，國家實(shí)行糧食生產(chǎn)補(bǔ)貼政策，充分調(diào)動(dòng)了糧食從業(yè)者種糧積極性，促進(jìn)黑龍江省糧食生產(chǎn)效率的提高。受SARS疫情的影響，2003年糧食生產(chǎn)效率大幅度下降。第2個(gè)階段2014年至今，糧食生產(chǎn)效率逐漸穩(wěn)定并趨于平均水平0.83，值得注意的是，為了進(jìn)一步調(diào)整糧食種植結(jié)構(gòu)，2016年國家取消了東北三省的玉米臨時(shí)收儲(chǔ)政策，導(dǎo)致2017年黑龍江省大部分地區(qū)糧食生產(chǎn)效率出現(xiàn)一個(gè)小幅度的銳減。

　　2000-2018年黑龍江省糧食生產(chǎn)綜合效率的平均水平為0.756，沒有達(dá)到DEA有效，表明黑龍江省的糧食綜合生產(chǎn)效率還有繼續(xù)提升的空間。從各年份的效率發(fā)展趨勢分析，2004年鶴崗;2006年佳木斯;2008年七臺(tái)河;2010年大慶、伊春、綏化、大興安嶺;2011年大興安嶺;2012年哈爾濱、雞西、大慶、伊春、佳木斯、綏化、大興安嶺;2014、2015、2016年的大興安嶺以及2018年雞西、雙鴨山、伊春分別在各年份期間的綜合效率值為1，表明糧食生產(chǎn)的投入和產(chǎn)出達(dá)到平衡，綜合效率達(dá)到最佳，糧食生產(chǎn)要素配置合理。

　　2.1.2 空間差異分析

　　為了更加有效直觀表現(xiàn)黑龍江省各地級(jí)市糧食生產(chǎn)綜合效率的空間演變，本文將綜合效率分為5組：0.95以上為高效率組、0.85～0.95為中高效率組、0.75～0.85為中效率組、0.65～0.75為中低效率組、0.65以下為低效率組，分別研究2003、2008、2013和2018年黑龍江省糧食生產(chǎn)綜合效率的空間分布態(tài)勢(如表3所示)?？傮w而言，4個(gè)年份黑龍江省糧食生產(chǎn)效率的平均值分別為0.524、0.832、0.810、0.842，說明糧食生產(chǎn)效率穩(wěn)步提高，空間差距逐步減小，仍具提升空間。

　　推薦閱讀：吉林糧食高等?？茖W(xué)校學(xué)報(bào)農(nóng)業(yè)技術(shù)類期刊征稿