城鎮(zhèn)化與糧食生態(tài)效率
——基于異質(zhì)性城鎮(zhèn)化的理論視角與實證檢驗

王帥，王亞靜

（中國農(nóng)業(yè)科學(xué)院農(nóng)業(yè)資源與農(nóng)業(yè)區(qū)劃研究所，北京 100081）

“十三五”期間，我國糧食產(chǎn)量持續(xù)穩(wěn)定在6.65億t以上，糧食數(shù)量安全得到基本保障。然而，糧食生產(chǎn)對化肥和農(nóng)藥的依賴程度較高，生態(tài)環(huán)境壓力不斷加大[1]，推進(jìn)綠色發(fā)展已成為普遍共識。與此同時，我國正在經(jīng)歷快速的城鎮(zhèn)化過程，城鎮(zhèn)常住人口占總?cè)丝诘谋戎匾阅昃?.04個百分點的速度快速增長，從1978年的17.92%提高到2021年的64.72%，人口城鎮(zhèn)化主導(dǎo)了中國城鎮(zhèn)化發(fā)展方向。尤其是自改革開放以來，農(nóng)村人口不斷流向城市，農(nóng)村勞動力結(jié)構(gòu)性短缺等問題日趨嚴(yán)重，威脅著糧食安全。在此背景下，我國開始探索以縣域為核心的中國特色城鎮(zhèn)化道路。就地城鎮(zhèn)化以中心村或小城鎮(zhèn)為依托，通過發(fā)展農(nóng)村二三產(chǎn)業(yè)實現(xiàn)農(nóng)村居民的就地非農(nóng)就業(yè)和市民化。相較于人口城鎮(zhèn)化下人口的單向流動，就地城鎮(zhèn)化一方面吸引進(jìn)城務(wù)工人員返鄉(xiāng)創(chuàng)業(yè)就業(yè)，另一方面帶動當(dāng)?shù)厝丝诘木偷鼐蜆I(yè)，淡化了城鄉(xiāng)二元結(jié)構(gòu)的同時加快了鄉(xiāng)村產(chǎn)業(yè)振興。人口城鎮(zhèn)化和就地城鎮(zhèn)化通過人口的不同流動影響著其他生產(chǎn)要素的配置，對糧食生產(chǎn)要素的綠色化利用水平產(chǎn)生異質(zhì)性影響。因此，本文以城鎮(zhèn)化對糧食生態(tài)效率的異質(zhì)性影響為核心問題，通過對這一問題的嘗試性解釋為城鎮(zhèn)化和糧食安全的協(xié)調(diào)發(fā)展提供參考。

糧食生態(tài)效率從生態(tài)效率延伸而來，通過將糧食生產(chǎn)的環(huán)境代價納入分析框架，反映了更真實的效率水平[2]。關(guān)于糧食生態(tài)效率測算的研究主要集中在兩方面。一是測算邊界。已有研究關(guān)于糧食生產(chǎn)投入要素的界定基本相差無幾，但對產(chǎn)出的界定則有所不同。期望產(chǎn)出有糧食產(chǎn)量和糧食產(chǎn)值兩種界定方法；非期望產(chǎn)出主要包含碳排放和面源污染兩種[3]，亦或兼而有之[4]。二是測算方法。早期測算方法主要基于“索洛余值”法，隨著生產(chǎn)前沿面分析方法的發(fā)展，逐漸過渡到隨機前沿生產(chǎn)函數(shù)[5]和數(shù)據(jù)包絡(luò)分析法[6]，但這兩種方法因需要提前假定生產(chǎn)函數(shù)形式和技術(shù)非效率項分布而受到質(zhì)疑。直到方向性距離函數(shù)演化到非徑向、非角度的松弛形式才從技術(shù)上解決了松弛量對效率評價的影響，此時納入了環(huán)境成本的生產(chǎn)效率被稱為綠色全要素生產(chǎn)率或生態(tài)效率[7]。糧食生產(chǎn)過程也會帶來一定的碳排放和面源污染，將二者進(jìn)行計量并納入到非期望產(chǎn)出中以反映糧食生產(chǎn)的環(huán)境負(fù)外部性，所測得的生產(chǎn)效率被稱為糧食生態(tài)效率。

已有研究對城鎮(zhèn)化與糧食生產(chǎn)之間的關(guān)系進(jìn)行了深入探討，但觀點尚不統(tǒng)一。部分學(xué)者認(rèn)為“城鎮(zhèn)化阻礙糧食生產(chǎn)”。即城鎮(zhèn)化的要素配置效應(yīng)使糧食生產(chǎn)投入要素減少，影響糧食生產(chǎn)。從勞動力角度看，城鎮(zhèn)化造成了農(nóng)業(yè)生產(chǎn)部門勞動力數(shù)量性和結(jié)構(gòu)性短缺[8]；從耕地角度看，工商資本下鄉(xiāng)投資農(nóng)村二三產(chǎn)業(yè)，使得農(nóng)地利用趨于“非糧化”[9]。而城鎮(zhèn)化失衡發(fā)展導(dǎo)致建設(shè)用地?zé)o序擴張，更是直接造成了耕地數(shù)量減少和質(zhì)量損害[10]。還有學(xué)者通過論證城鎮(zhèn)化所導(dǎo)致的水資源短缺[11]、種植結(jié)構(gòu)變動[12]、政策空間外部性[13]來證明城鎮(zhèn)化對糧食生產(chǎn)的阻礙作用。與此相反，部分學(xué)者認(rèn)為“城鎮(zhèn)化促進(jìn)糧食生產(chǎn)”，即城鎮(zhèn)化與糧食安全保障的目標(biāo)是可以相容的[14]，城鎮(zhèn)化帶來的農(nóng)業(yè)技術(shù)進(jìn)步、農(nóng)地規(guī)?；?jīng)營、經(jīng)營方式轉(zhuǎn)變等效應(yīng)也不容忽視。城鎮(zhèn)化對糧食生產(chǎn)的資本投入與技術(shù)進(jìn)步效應(yīng)能有效替代要素利用效應(yīng)，從長期看對糧食生產(chǎn)能力提高具有積極作用[15]。甚至單就城鎮(zhèn)化導(dǎo)致農(nóng)村勞動力短缺而影響糧食生產(chǎn)這一觀點也有學(xué)者提出批駁，如王躍梅等[16]認(rèn)為城鎮(zhèn)化引導(dǎo)了農(nóng)村剩余勞動力向城鎮(zhèn)流動，緩解了糧食生產(chǎn)的內(nèi)卷化現(xiàn)象，提高了勞動效率，促進(jìn)了糧食生產(chǎn)。

綜上，已有研究為本文提供了理論參考，但仍存在以下兩點不足。第一，對關(guān)鍵概念理解不夠深入，導(dǎo)致模型估計可能存在偏誤。事實上，當(dāng)前的糧食生產(chǎn)效率會受到前期影響，僅構(gòu)建靜態(tài)模型可能會造成估計的非一致性。第二，缺乏對城鎮(zhèn)化類型的考察，造成研究結(jié)論存在較大分歧。既有研究多從人口城鎮(zhèn)化視角審視城鎮(zhèn)化與糧食生產(chǎn)的關(guān)系，鮮有文獻(xiàn)關(guān)注城鎮(zhèn)化發(fā)展類型的差異。本文基于“城鎮(zhèn)化—要素配置—糧食綠色生產(chǎn)”之間的關(guān)系構(gòu)建理論分析框架，在利用超效率SBM模型測算糧食生態(tài)效率的基礎(chǔ)上，運用雙向固定效應(yīng)模型剖析城鎮(zhèn)化對糧食生態(tài)效率的影響機制。

1 理論分析與研究假設(shè)

1.1 人口城鎮(zhèn)化與糧食生態(tài)效率

傳統(tǒng)的人口城鎮(zhèn)化以農(nóng)村人口向城鎮(zhèn)的單向流動為主要路徑，通過對勞動力資源的重新配置對其他要素配置產(chǎn)生影響。人口城鎮(zhèn)化下大量農(nóng)村青壯年進(jìn)城轉(zhuǎn)為非農(nóng)就業(yè)，直接導(dǎo)致種糧勞動力結(jié)構(gòu)性短缺。這種情況對糧食生態(tài)效率的負(fù)向影響主要表現(xiàn)在以下幾個方面。第一，為了緩解勞動力短缺約束，農(nóng)戶會加大化肥、農(nóng)藥等替代性要素的投入強度，造成投入要素的不合理配置和浪費，提高糧食生產(chǎn)碳排放和面源污染強度的同時降低生產(chǎn)效率，進(jìn)而降低糧食生態(tài)效率。第二，從農(nóng)業(yè)技術(shù)角度來看，農(nóng)村青壯年勞動力不足，“老人農(nóng)業(yè)”成為糧食生產(chǎn)的常態(tài)，其直接后果是降低生產(chǎn)效率，同時還制約了綠色生產(chǎn)技術(shù)的采納和應(yīng)用，阻礙糧食生態(tài)效率提高。第三，從農(nóng)地角度來看，人口城鎮(zhèn)化的無序發(fā)展侵占耕地、耕地占優(yōu)補劣等投機行為對農(nóng)地數(shù)量和質(zhì)量的雙重影響，不僅會直接減少期望產(chǎn)出，同時還可能導(dǎo)致農(nóng)地面源污染，增加農(nóng)地非期望產(chǎn)出，降低糧食生態(tài)效率。基于此提出假說1：

H1：人口城鎮(zhèn)化對糧食生態(tài)效率具有顯著的負(fù)向影響。

1.2 就地城鎮(zhèn)化與糧食生態(tài)效率

就地城鎮(zhèn)化以中心鎮(zhèn)或村為依托，以相關(guān)現(xiàn)代產(chǎn)業(yè)體系吸納農(nóng)民的就地就近就業(yè)，通過人口回流帶來的系列效應(yīng)對糧食生態(tài)效率產(chǎn)生積極影響。在此情境下農(nóng)戶的農(nóng)業(yè)經(jīng)營存在兩種選擇，分別是兼業(yè)經(jīng)營和流轉(zhuǎn)土地。首先，兩種情況都能夠?qū)崿F(xiàn)農(nóng)村勞動力數(shù)量上的補充，一定程度上減少了勞動替代性要素的投入，促進(jìn)農(nóng)機、化肥、農(nóng)藥等要素合理使用，提升要素配置效率同時減少糧食生產(chǎn)碳排放和面源污染，提高了糧食生態(tài)效率；其次，兩種方式也都能夠促進(jìn)農(nóng)村勞動力質(zhì)量上的提升，返鄉(xiāng)就業(yè)創(chuàng)業(yè)人員具有更高的人力資本水平、更強的環(huán)保意識和生態(tài)自覺性，能夠加快糧食綠色生產(chǎn)的新知識、新技術(shù)和新模式的普及應(yīng)用，提高生產(chǎn)效率同時減少非期望產(chǎn)出，提升糧食生態(tài)效率。而從收入效應(yīng)上看，就地城鎮(zhèn)化對兼業(yè)農(nóng)戶與流轉(zhuǎn)農(nóng)戶的影響則不同。對于兼業(yè)農(nóng)戶而言，就地城鎮(zhèn)化增加了其工資性收入，緩解了糧食生產(chǎn)資金約束，促進(jìn)其綠色生產(chǎn)技術(shù)投資，提高糧食生態(tài)效率；選擇土地流轉(zhuǎn)則促進(jìn)了糧食規(guī)?；?jīng)營，發(fā)揮規(guī)模優(yōu)勢降低綠色技術(shù)投資成本，促進(jìn)新技術(shù)、新模式的應(yīng)用，減少分散農(nóng)戶不合理施肥用藥等行為帶來的污染，提高糧食生態(tài)效率。基于此提出假說2：

H2：就地城鎮(zhèn)化對糧食生態(tài)效率具有顯著的正向影響。

1.3 異質(zhì)性城鎮(zhèn)化對糧食生態(tài)效率的影響

城鎮(zhèn)化的本質(zhì)是人口在城鄉(xiāng)之間流動，無論以何種經(jīng)濟現(xiàn)象或指標(biāo)所表征的城鎮(zhèn)化，都不可避免地通過人口在城鄉(xiāng)之間的流動所引發(fā)的生產(chǎn)要素重新配置而對糧食生態(tài)效率產(chǎn)生影響。人口城鎮(zhèn)化下農(nóng)村人口向城鎮(zhèn)單向轉(zhuǎn)移，影響了糧食生產(chǎn)勞動力、技術(shù)、農(nóng)地等要素的合理配置，降低效率的同時增加了碳排放和面源污染等非期望產(chǎn)出，降低了糧食生態(tài)效率。就地城鎮(zhèn)化通過開發(fā)當(dāng)?shù)刭Y源發(fā)展鄉(xiāng)村二三產(chǎn)業(yè)，引導(dǎo)轉(zhuǎn)移人口返鄉(xiāng)就業(yè)創(chuàng)業(yè)和未轉(zhuǎn)移人口就地就近就業(yè)，激發(fā)了鄉(xiāng)村生產(chǎn)要素內(nèi)生活力。通過改善勞動力、技術(shù)、農(nóng)地、資金等要素的配置方式，降低了碳排放和面源污染水平，提高了糧食生態(tài)效率（圖1）?；诖颂岢黾僬f3：

圖1 城鎮(zhèn)化影響糧食生態(tài)效率的理論框架Fig. 1 Theoretical framework of the impacts of urbanization on grain eco-efficiency

H3：城鎮(zhèn)化對糧食生態(tài)效率具有顯著影響，且因城鎮(zhèn)化發(fā)展路徑的異質(zhì)性而有所差異。

2 研究方法

2.1 模型設(shè)計

2.1.1 超效率SBM-Undesirable模型 超效率SBMUndesirable模型既克服了傳統(tǒng)DEA模型預(yù)設(shè)生產(chǎn)函數(shù)和關(guān)鍵技術(shù)變量分布形式的不足又提高了普通SBM模型的解釋力，還能夠?qū)⒎瞧谕a(chǎn)出引入到模型當(dāng)中對綠色生產(chǎn)效率進(jìn)行測度[17]，因此本文選擇包含非期望產(chǎn)出的超效率SBM模型測算糧食生態(tài)效率。糧食生產(chǎn)投入采用廣義權(quán)重系數(shù)法計算。超效率SBM- Undesirable模型可設(shè)置為：

式（1）表示效率測算模型，式（2）是效率測算模型的系列約束條件的集合。其中，π表示糧食生態(tài)效率；決策單元DUM由投入要素、期望產(chǎn)出和非期望產(chǎn)出組成，分別由m、c1、c2表示；n表示決策單元數(shù)，n=30；xik、yαk、bβk分別表示投入矩陣、期望產(chǎn)出矩陣、非期望產(chǎn)出矩陣的元素。

2.1.2 雙向固定效應(yīng)模型 為進(jìn)一步研究不同城鎮(zhèn)化類型對糧食生態(tài)效率的影響，避免遺漏變量問題，采用雙向固定效應(yīng)模型進(jìn)行檢驗，形式如下：

式（3）、（4）中i表示省份，t表示年份，Y表示糧食生態(tài)效率，X表示傳統(tǒng)城鎮(zhèn)化，X’表示就地城鎮(zhèn)化，C表示系列控制變量，pi表示地區(qū)固定效應(yīng)項，yt表示時間固定效應(yīng)，ρ表示待估參數(shù)，ε表示隨機擾動項。

2.2 變量選擇

糧食生態(tài)效率測算變量的選取及解釋如表1所示；城鎮(zhèn)化對糧食生態(tài)效率影響機制的變量選取及解釋如表2所示。中國保障糧食安全的基本目標(biāo)是“口糧絕對安全”，其基礎(chǔ)仍是提高糧食自給率。因此糧食生態(tài)效率測算中所用到的產(chǎn)量數(shù)據(jù)指的是三大主糧（稻谷、小麥和玉米）的產(chǎn)量。

表1 糧食生態(tài)效率測算相關(guān)指標(biāo)選取與解釋Table 1 Selection and interpretation of relevant indicators for grain eco-efficiency measurement

表2 回歸模型變量選擇與解釋Table 2 Regression model variable selection and interpretation

2.3 數(shù)據(jù)來源與說明

2006—2019年中國31個?。▍^(qū)、市）（不包括港澳臺）的面板數(shù)據(jù)來源于《中國統(tǒng)計年鑒》《中國農(nóng)村統(tǒng)計年鑒》《中國水利統(tǒng)計年鑒》《中國農(nóng)業(yè)機械工業(yè)年鑒》《中國人口和就業(yè)統(tǒng)計年鑒》以及各省市統(tǒng)計年鑒。樣本時間間隔的選取基于三個原因：第一，2006年以前部分省份的相關(guān)指標(biāo)數(shù)據(jù)缺失較為嚴(yán)重，且統(tǒng)計口徑差異大，而2019年后的數(shù)據(jù)受國際局勢等的干擾，用于模型分析可能會產(chǎn)生難以消除的數(shù)據(jù)噪音干擾，影響估計結(jié)果的可靠性；第二，2006年以后就地城鎮(zhèn)化快速發(fā)展并逐漸成為獨立于人口城鎮(zhèn)化的新型城鎮(zhèn)化發(fā)展路徑，并進(jìn)入政策視野和學(xué)術(shù)研究；第三，2006年以后國家糧食供求開始出現(xiàn)緊平衡狀態(tài)，使用該節(jié)點后的數(shù)據(jù)可更為深入地考察不同城鎮(zhèn)化路徑對糧食生態(tài)效率的影響。為便于比較分析和穩(wěn)健性檢驗，本文對各糧食生產(chǎn)功能區(qū)進(jìn)行細(xì)分考察。根據(jù)《國家糧食安全中長期規(guī)劃綱要（2008—2020）》的標(biāo)準(zhǔn)劃分為糧食主產(chǎn)區(qū)、主銷區(qū)和產(chǎn)銷平衡區(qū)，糧食主產(chǎn)區(qū)包括遼寧、吉林、黑龍江、內(nèi)蒙古、河北、山東、安徽、江蘇、江西、河南、湖南、四川和湖北13個省（區(qū)）；主銷區(qū)包括：北京、上海、天津、浙江、海南、廣東和福建7個?。ㄊ校划a(chǎn)銷平衡區(qū)包括：山西、廣西、重慶、云南、貴州、西藏、陜西、甘肅、青海、寧夏和新疆11個?。▍^(qū)、市）。

3 結(jié)果與分析

3.1 糧食生態(tài)效率的時空差異分析

3.1.1 全國層面 2006—2019年間，我國糧食生態(tài)效率總體呈“U”型波動變化。2006—2012年先下降后提高，波動范圍較大，總體水平較低；2013年以后，經(jīng)歷了緩慢的爬坡期，至2019年達(dá)到近14年來的最高水平（圖2）?？傮w上，我國糧食生態(tài)效率的變化趨勢為“波動中略有提升”，效率值從2006年的0.803提升到2019年的0.908。隨著糧食生產(chǎn)的資源環(huán)境約束趨緊，國家至各級地方政府都制定了農(nóng)業(yè)綠色發(fā)展規(guī)劃及相關(guān)政策措施，加大力度推廣水肥一體化、保護(hù)性耕作、秸稈綜合利用等綠色技術(shù)，推動了糧食生態(tài)效率的提升。

圖2 2006—2019年我國糧食生態(tài)效率變化趨勢Fig. 2 Trend of grain eco-efficiency in China from 2006 to 2019

3.1.2 分區(qū)層面 三大糧食生產(chǎn)功能區(qū)糧食生態(tài)效率的變動同全國層面的波動基本一致（圖2）。從“十一五”到“十三五”期間，三大糧食生產(chǎn)功能區(qū)的糧食生態(tài)效率均有所提升，但呈現(xiàn)出明顯的地區(qū)差異性。從均值來看，主產(chǎn)區(qū)的糧食生態(tài)效率最高，產(chǎn)銷平衡區(qū)、主銷區(qū)次之。三大區(qū)域的糧食生態(tài)效率自2010年降至最低值之后開始了緩慢發(fā)展態(tài)勢，經(jīng)歷了近7年的爬坡期后才有了顯著提高。到2017年以后才開始有顯著的提高?？傮w來看，各區(qū)域的糧食生態(tài)效率水平都得到了穩(wěn)步推進(jìn)。

3.2 城鎮(zhèn)化對糧食生態(tài)效率的影響機制

模型1～4基于雙向固定效應(yīng)模型（表3）。從模型1和模型2的結(jié)果來看，在未引入其他控制變量的情況下，人口城鎮(zhèn)化對糧食生態(tài)效率具有顯著的負(fù)向影響，就地城鎮(zhèn)化對糧食生態(tài)效率具有顯著的正向影響。模型3和模型4的估計結(jié)果則表明，引入關(guān)鍵控制變量以后，這種影響關(guān)系依然顯著，且影響程度有所提升。模型3～4的擬合優(yōu)度與模型1～2相比也有明顯提升。模型5～8首先對人口城鎮(zhèn)化和就地城鎮(zhèn)化進(jìn)行一階差分處理，將一階滯后項作為工具變量進(jìn)行兩階段最小二乘回歸（表4）。IV-2SLS的回歸結(jié)果與雙向固定效應(yīng)模型回歸結(jié)果基本一致。尤其是引入控制變量后，運用工具變量法得到的城鎮(zhèn)化水平估計系數(shù)顯著上升，表明忽視內(nèi)生性問題將出現(xiàn)明顯的估計偏差。因此，本文以工具變量-兩階段最小二乘法的回歸結(jié)果作為基準(zhǔn)回歸結(jié)果對其進(jìn)行解釋。

表3 雙向固定效應(yīng)模型回歸結(jié)果Table 3 Two-way fixed effect model regression results

表4 工具變量-兩階段最小二乘法回歸結(jié)果Table 4 Instrumental variables-two-stage least squares regression results

3.2.1 基準(zhǔn)回歸結(jié)果分析 1）人口城鎮(zhèn)化對糧食生態(tài)效率的影響。表4中模型7結(jié)果表明，人口城鎮(zhèn)化對糧食生態(tài)效率的影響彈性為1.071。我們可從人口城鎮(zhèn)化對農(nóng)戶生產(chǎn)決策的影響對此做出解釋。人口城鎮(zhèn)化導(dǎo)致農(nóng)村勞動力減少，勞動力市場供求關(guān)系的改變使得農(nóng)業(yè)雇工價格居高不下，糧食生產(chǎn)面臨勞動力和雇工成本的雙重壓力[20]，理性農(nóng)戶將從兩個環(huán)節(jié)對此壓力做出響應(yīng)：首先調(diào)整種植結(jié)構(gòu)，進(jìn)而通過替代性生產(chǎn)資料投入以最大限度地減少勞動投入。當(dāng)雇工成本過高時，農(nóng)戶會更傾向于種植耗費人力勞動更少的糧食作物，減少經(jīng)濟作物種植[21]，同時提高各生產(chǎn)環(huán)節(jié)的機械化水平以減少人力勞動。例如，在整地、播種、收割等環(huán)節(jié)提高機械作業(yè)比例，在田間管理中大量施用農(nóng)藥以減少人工除草成本[22]、提高化肥施用強度以降低施肥頻率，其結(jié)果一方面導(dǎo)致了糧食生產(chǎn)碳排放和面源污染程度增加[23]，另一方面還可能因不合理的施肥施藥降低糧食產(chǎn)量，降低糧食生產(chǎn)生態(tài)效率。此外，城鎮(zhèn)化的就業(yè)篩選機制只給人力資本水平較高者流入城市的機會[24]，年齡較大、受教育水平較低者只能繼續(xù)留守農(nóng)村，這也將制約綠色生產(chǎn)技術(shù)在糧食生產(chǎn)中的應(yīng)用，阻礙糧食生態(tài)效率的提高。

2）就地城鎮(zhèn)化對糧食生態(tài)效率的影響。表4中模型8的結(jié)果表明，就地城鎮(zhèn)化對糧食生態(tài)效率的影響彈性為0.195?？梢詮某擎?zhèn)化對農(nóng)民投資決策和農(nóng)村人力資本水平的影響做出解釋。兩種城鎮(zhèn)化的發(fā)展能夠增加農(nóng)民收入，但在人口城鎮(zhèn)化路徑作用下，進(jìn)城農(nóng)民已經(jīng)脫離或半脫離農(nóng)業(yè)生產(chǎn)，收入增加對農(nóng)戶的農(nóng)業(yè)投資決策影響較小，對教育、醫(yī)療等其他生活投資影響更為明顯[25]。就地城鎮(zhèn)化通過發(fā)展鄉(xiāng)村工業(yè)、旅游業(yè)、康養(yǎng)業(yè)等二、三產(chǎn)業(yè)實現(xiàn)農(nóng)民就地就業(yè)，提高了農(nóng)民非農(nóng)收入[26]，對農(nóng)戶農(nóng)業(yè)生產(chǎn)投資的促進(jìn)作用更為明顯，能有效促進(jìn)農(nóng)戶綠色技術(shù)投資，加快節(jié)水、節(jié)肥、節(jié)藥等節(jié)約型生產(chǎn)技術(shù)和測土配方施肥、水肥一體化等增效型生產(chǎn)技術(shù)的應(yīng)用[27]，提高了要素使用效率同時降低了污染水平，提高了糧食生態(tài)效率。此外，人力資本水平較高、富有創(chuàng)新精神的農(nóng)戶返鄉(xiāng)將改善農(nóng)村人力資本的整體水平，能夠直接提高要素配置效率，促進(jìn)農(nóng)戶綠色技術(shù)采納，推動糧食生產(chǎn)生態(tài)效率前沿面外移。

無論是否引入控制變量，人口城鎮(zhèn)化對糧食生態(tài)效率的影響都要遠(yuǎn)高于就地城鎮(zhèn)化。究其原因，我國自改革開放以來，城鎮(zhèn)化的推進(jìn)路徑長期以傳統(tǒng)的人口城鎮(zhèn)化為主導(dǎo)，存在明顯的路徑依賴現(xiàn)象，就地城鎮(zhèn)化雖然在黨的十八大以后逐漸為政策所倡導(dǎo)，但僅在近幾年推進(jìn)新型城鎮(zhèn)化過程中才逐漸發(fā)展起來，其發(fā)展程度遠(yuǎn)低于人口城鎮(zhèn)化。描述性統(tǒng)計結(jié)果可以作為驗證，2019年人口城鎮(zhèn)化發(fā)展水平的均值為53.67%，而就地城鎮(zhèn)化僅有18.91%。因此，人口城鎮(zhèn)化與就地城鎮(zhèn)化對糧食生態(tài)效率影響程度存在明顯差異，這與其發(fā)展程度是吻合的。

3）控制變量對糧食生態(tài)效率的影響?？刂谱兞康幕貧w結(jié)果符合一般性的經(jīng)濟預(yù)期，與已有研究成果基本一致。農(nóng)地經(jīng)營規(guī)模對糧食生態(tài)效率有顯著的正向影響，即農(nóng)地經(jīng)營規(guī)模的擴大能顯著提高糧食生態(tài)效率。財政支農(nóng)水平對糧食生態(tài)效率也具有顯著正向影響，隨著農(nóng)業(yè)綠色發(fā)展戰(zhàn)略的提出以及國家對糧食生態(tài)安全重視程度的提高，財政支持農(nóng)業(yè)綠色發(fā)展的資金投入逐年提高，有利于糧食生態(tài)效率的提升。農(nóng)業(yè)機械化水平對糧食生態(tài)效率的提升有顯著的負(fù)向影響，農(nóng)業(yè)機械化水平的提升導(dǎo)致農(nóng)業(yè)生產(chǎn)中化石能源消耗增加，對糧食生態(tài)效率產(chǎn)生負(fù)向影響。受災(zāi)情況帶來了糧食生態(tài)效率的顯著降低，自然災(zāi)害會給糧食的生產(chǎn)造成破壞性影響，阻礙了糧食生態(tài)效率的提高。農(nóng)村人力資本水平對糧食生態(tài)效率的影響雖然不顯著，但其符號方向仍具有參考意義。相對較高的農(nóng)村人力資本水平會增強農(nóng)戶對綠色生產(chǎn)的認(rèn)知，提高綠色生產(chǎn)技術(shù)的采納和應(yīng)用水平，對糧食生態(tài)效率產(chǎn)生正向影響。

3.2.2 分糧食生產(chǎn)功能區(qū)的考察 現(xiàn)實中，不同糧食生產(chǎn)功能區(qū)在城鎮(zhèn)化水平、種植結(jié)構(gòu)、產(chǎn)業(yè)結(jié)構(gòu)、自然地理條件以及資源稟賦等諸多方面都存在明顯的系統(tǒng)性差異，其對糧食生態(tài)效率的影響路徑也可能存在區(qū)域性差異，因此有必要對城鎮(zhèn)化影響糧食生態(tài)效率的區(qū)域異質(zhì)性進(jìn)行重新估計（表5）。

表5 區(qū)域異質(zhì)性分析Table 5 Regional heterogeneity analysis

1）主產(chǎn)區(qū)。通過對糧食主產(chǎn)區(qū)進(jìn)行分樣本回歸可以發(fā)現(xiàn)，人口城鎮(zhèn)化在主產(chǎn)區(qū)對糧食生態(tài)效率表現(xiàn)出顯著的正向影響，這一結(jié)果與全國樣本中的影響效應(yīng)不一致。由于糧食作物的種植比例本身就較高，“趨糧化”結(jié)構(gòu)調(diào)整策略在應(yīng)對勞動力短缺時彈性不足[28]。與此同時，雇工成本居高不下時農(nóng)戶分散經(jīng)營的比較收益下降，刺激了農(nóng)戶土地流轉(zhuǎn)，加速了耕地規(guī)?；?、集約化經(jīng)營[29]，提高了要素配置效率和綠色生產(chǎn)水平，促進(jìn)了糧食生態(tài)效率的提高。另外，糧食主產(chǎn)區(qū)是中國糧食生產(chǎn)的核心區(qū)域，承擔(dān)著70%以上的糧食產(chǎn)量和80%左右的商品糧供給，國家實行了最嚴(yán)格的耕地保護(hù)政策，并出臺了糧食主產(chǎn)區(qū)利益補償、產(chǎn)糧大縣獎勵、新增農(nóng)業(yè)補貼向糧食主產(chǎn)區(qū)和優(yōu)勢產(chǎn)區(qū)集中等多項政策，一定程度上也紓解了城鎮(zhèn)化發(fā)展對主產(chǎn)區(qū)農(nóng)地、技術(shù)、資金的負(fù)向影響。

2）非主產(chǎn)區(qū)。糧食非主產(chǎn)區(qū)的回歸結(jié)果與全國層面基本一致。在主銷區(qū)，人口城鎮(zhèn)化對糧食生態(tài)效率的負(fù)向作用更強。從現(xiàn)實來看，我國糧食主銷區(qū)多集中在東部經(jīng)濟較為發(fā)達(dá)的省份，其糧食播種面積和糧食產(chǎn)量均占到全國水平的10%左右。但2006年主銷區(qū)內(nèi)各省人口城鎮(zhèn)化平均水平為66.39%，到2019年達(dá)到75.07%，年均增速0.67%，城鎮(zhèn)化發(fā)展水平和速度都明顯高于另外兩區(qū)?？焖俚某擎?zhèn)化進(jìn)程使得耕地面積不斷下降，農(nóng)業(yè)人口短缺也更為嚴(yán)重，勞動力短缺導(dǎo)致相關(guān)生產(chǎn)要素配置效率低、面源污染和碳排放嚴(yán)重，對糧食生態(tài)效率的負(fù)向影響更為顯著。

在產(chǎn)銷平衡區(qū)，人口城鎮(zhèn)化、就地城鎮(zhèn)化對糧食生態(tài)效率的作用方向與全國層面均保持一致，但程度都略低。從現(xiàn)實來看，產(chǎn)銷平衡區(qū)主要集中在我國中西部省份，糧食產(chǎn)量約占全國總產(chǎn)量的20%左右，主要目標(biāo)是“保持應(yīng)有的自給率”，其糧食安全形勢頗具特殊性。產(chǎn)銷平衡區(qū)多為我國邊疆地區(qū)以及丘陵地區(qū)等生態(tài)環(huán)境脆弱區(qū)，糧食生產(chǎn)的自然條件惡劣與經(jīng)濟發(fā)展薄弱相疊加，鄉(xiāng)村二、三產(chǎn)業(yè)發(fā)展落后，兩種類型的城鎮(zhèn)化發(fā)展水平都較低，城鎮(zhèn)化發(fā)展對糧食生態(tài)效率的影響程度較小。

3.3 穩(wěn)健性檢驗

本文進(jìn)行穩(wěn)健性檢驗主要有兩個思路。一是從樣本角度，通過縮尾處理剔除掉離群值，避免樣本異常值導(dǎo)致的回歸曲線與實際的偏離。另外，對樣本進(jìn)行分組并重新回歸也是一種從樣本角度檢驗穩(wěn)健性的有效方法，這在上一節(jié)對各糧食生產(chǎn)功能區(qū)的考察中已經(jīng)完成，在此不再贅述。二是從模型角度，通過更換計量經(jīng)濟模型對樣本進(jìn)行重新回歸。由于受自然因素的影響以及農(nóng)業(yè)技術(shù)進(jìn)步和生產(chǎn)方式變遷的“路徑依賴”特征，糧食生態(tài)效率的變化是一個持續(xù)調(diào)整的動態(tài)過程。因此，通過將糧食生態(tài)效率的滯后項納入模型中進(jìn)行控制，同時為了避免內(nèi)生性問題的干擾，采用系統(tǒng)廣義矩估計（GMM）的方法以被解釋變量滯后兩期及以上作為工具變量進(jìn)行回歸。

3.3.1 剔除奇異樣本：非刪失縮尾處理 由于糧食生態(tài)效率的值域較小，回歸結(jié)果容易受到異常值的影響，有必要對離群值進(jìn)行檢驗和處理。為避免人為剔除樣本可能導(dǎo)致的選擇偏誤問題，通過繪制箱線圖來檢驗樣本異常值，并將超出箱體邊界線的觀測值視為異常值。根據(jù)箱線圖的指示，對樣本進(jìn)行非刪失雙邊1%縮尾處理后重新進(jìn)行回歸。從表6模型1～2的結(jié)果可以看出，在縮尾處理后影響關(guān)系與基準(zhǔn)回歸基本一致，表明基準(zhǔn)模型的估計結(jié)果是穩(wěn)健可靠的。

表6 穩(wěn)健性檢驗Table 6 Robustness test

3.3.2 更換計量模型：系統(tǒng)廣義矩估計 考慮糧食生態(tài)效率的動態(tài)特征，采用系統(tǒng)廣義矩估計（GMM）作為新的計量模型，以被解釋變量滯后兩期及以上作為工具變量重新回歸。模型3～4的AR（1）和AR（2）檢驗均表明模型殘差項存在一階序列相關(guān)，二階序列不相關(guān)。Sargan-Hansen檢驗表明工具變量是有效的，證明GMM方法較好地處理了模型內(nèi)生性問題。糧食生態(tài)效率的時間滯后項系數(shù)顯著為正，表明糧食生態(tài)效率的變化存在明顯的路徑依賴性，表明其增長受到歷史因素影響，進(jìn)一步驗證了前文的理論分析結(jié)果，檢驗了基準(zhǔn)模型和研究結(jié)論的穩(wěn)健性。

4 結(jié)論與對策建議

4.1 研究結(jié)論

本文通過超效率SBM模型測算我國31個?。▍^(qū)、市）（不包括港澳臺）的糧食生態(tài)效率，在此基礎(chǔ)上運用雙向固定效應(yīng)模型等方法對不同城鎮(zhèn)化類型影響糧食生態(tài)效率的作用機制進(jìn)行分析和檢驗，得出如下結(jié)論：

1）從糧食生態(tài)效率測算結(jié)果來看，中國糧食生態(tài)效率水平較低，存在明顯的時間波動和區(qū)域差異。從時間角度看，糧食生態(tài)效率總體呈“W”型波動，各糧食生產(chǎn)功能區(qū)總體呈“V”型波動；從空間角度看，三大糧食生產(chǎn)功能區(qū)糧食綠色生產(chǎn)水平存在一定差距，由高到低依次是產(chǎn)銷平衡區(qū)、主產(chǎn)區(qū)、主銷區(qū)。

2）基于異質(zhì)性城鎮(zhèn)化路徑的視角，人口城鎮(zhèn)化阻礙了糧食生態(tài)效率的提高，就地城鎮(zhèn)化促進(jìn)了糧食生態(tài)效率的提高。人口城鎮(zhèn)化的影響程度高于就地城鎮(zhèn)化，這可能與中國長期以來城鎮(zhèn)化路徑以人口城鎮(zhèn)化為主導(dǎo)有關(guān)，城鎮(zhèn)化發(fā)展存在一定的路徑依賴。

3）基于糧食生產(chǎn)功能分區(qū)的視角，在主銷區(qū)和產(chǎn)銷平衡區(qū)，人口城鎮(zhèn)化對糧食生態(tài)效率均表現(xiàn)出顯著的負(fù)向影響，而在主產(chǎn)區(qū)則表現(xiàn)為顯著的正向影響，這可能與糧食主產(chǎn)區(qū)糧食種植面積較大、糧食種植結(jié)構(gòu)調(diào)整彈性不足有關(guān)；在主銷區(qū)，人口城鎮(zhèn)化對糧食生態(tài)效率的負(fù)向作用更強，這可能與主銷區(qū)耕地面積少，糧食生產(chǎn)吸納就業(yè)能力有限，人口城鎮(zhèn)化進(jìn)程早、程度高，對糧食綠色生產(chǎn)要素吸納作用更強有關(guān)。

4.2 對策建議

本文的研究結(jié)論對理解城鎮(zhèn)化影響糧食生態(tài)效率的作用機制具有參考意義，同時也為提升糧食綠色生產(chǎn)水平、保障國家糧食生態(tài)安全提供了新思路。因地制宜推進(jìn)城鎮(zhèn)化，引導(dǎo)城鎮(zhèn)化與糧食綠色生產(chǎn)能力協(xié)調(diào)發(fā)展是城鎮(zhèn)化快速發(fā)展背景下保障我國糧食安全的基本思路。具體對策建議如下：

1）主產(chǎn)區(qū)承擔(dān)著糧食保供壓力，而大多省份又是經(jīng)濟發(fā)展滯后區(qū)，人口城鎮(zhèn)化發(fā)展水平較低，城鎮(zhèn)吸納農(nóng)村剩余勞動力轉(zhuǎn)移的能力不足。因此，主產(chǎn)區(qū)可以通過引導(dǎo)農(nóng)民就地創(chuàng)業(yè)就業(yè)、扶持鄉(xiāng)村產(chǎn)業(yè)，推動就地城鎮(zhèn)化發(fā)展，以提高鄉(xiāng)村承載力和內(nèi)在發(fā)展動能。同時，還必須要警惕就地城鎮(zhèn)化無序發(fā)展占用土地引起的風(fēng)險，嚴(yán)格落實土地用途管制、基本農(nóng)田保護(hù)以及土地利用總體規(guī)劃等措施，嚴(yán)守耕地紅線。

2）主銷區(qū)經(jīng)濟發(fā)展水平較高，人口城鎮(zhèn)化發(fā)展已不可避免，需要通過發(fā)展社會化服務(wù)補齊農(nóng)村勞動力短板。同時增加社會化服務(wù)中綠色技術(shù)的供給，推廣節(jié)約增效型農(nóng)業(yè)技術(shù)和裝備，提高生產(chǎn)要素的合理利用水平，減少糧食生產(chǎn)過程中的污染物排放，在保證糧食產(chǎn)出增長的同時降低污染水平。

3）產(chǎn)銷平衡區(qū)多為生態(tài)脆弱的邊遠(yuǎn)省份，城鎮(zhèn)化發(fā)展水平低，糧食生產(chǎn)能力有限。產(chǎn)銷平衡區(qū)應(yīng)通過加強高標(biāo)準(zhǔn)農(nóng)田建設(shè)、提升技術(shù)服務(wù)水平、提高糧食品種適宜性等途徑確保糧食自給能力。同時，還要關(guān)注農(nóng)村留守勞動力的人力資本水平，通過新型職業(yè)農(nóng)民教育培訓(xùn)提升勞動力質(zhì)量以提高農(nóng)戶綠色生產(chǎn)參與程度和生態(tài)自覺性，保障城鎮(zhèn)化發(fā)展與糧食生態(tài)效率的協(xié)調(diào)平衡。