王帥，王亞靜

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

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

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

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

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

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

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

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

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

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

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

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

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

城鎮(zhèn)化的本質(zhì)是人口在城鄉(xiāng)之間流動(dòng)，無(wú)論以何種經(jīng)濟(jì)現(xiàn)象或指標(biāo)所表征的城鎮(zhèn)化，都不可避免地通過(guò)人口在城鄉(xiāng)之間的流動(dòng)所引發(fā)的生產(chǎn)要素重新配置而對(duì)糧食生態(tài)效率產(chǎn)生影響。人口城鎮(zhèn)化下農(nóng)村人口向城鎮(zhèn)單向轉(zhuǎn)移，影響了糧食生產(chǎn)勞動(dòng)力、技術(shù)、農(nóng)地等要素的合理配置，降低效率的同時(shí)增加了碳排放和面源污染等非期望產(chǎn)出，降低了糧食生態(tài)效率。就地城鎮(zhèn)化通過(guò)開(kāi)發(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)生活力。通過(guò)改善勞動(dòng)力、技術(shù)、農(nóng)地、資金等要素的配置方式，降低了碳排放和面源污染水平，提高了糧食生態(tài)效率（圖1）。基于此提出假說(shuō)3：

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

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

2 研究方法

2.1 模型設(shè)計(jì)

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

式（1）表示效率測(cè)算模型，式（2）是效率測(cè)算模型的系列約束條件的集合。其中，π表示糧食生態(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)化類(lèi)型對(duì)糧食生態(tài)效率的影響，避免遺漏變量問(wèn)題，采用雙向固定效應(yīng)模型進(jìn)行檢驗(yàn)，形式如下：

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

2.2 變量選擇

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

表1 糧食生態(tài)效率測(cè)算相關(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ù)來(lái)源與說(shuō)明

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

3 結(jié)果與分析

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

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

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

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

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

模型1～4基于雙向固定效應(yīng)模型（表3）。從模型1和模型2的結(jié)果來(lái)看，在未引入其他控制變量的情況下，人口城鎮(zhèn)化對(duì)糧食生態(tài)效率具有顯著的負(fù)向影響，就地城鎮(zhèn)化對(duì)糧食生態(tài)效率具有顯著的正向影響。模型3和模型4的估計(jì)結(jié)果則表明，引入關(guān)鍵控制變量以后，這種影響關(guān)系依然顯著，且影響程度有所提升。模型3～4的擬合優(yōu)度與模型1～2相比也有明顯提升。模型5～8首先對(duì)人口城鎮(zhèn)化和就地城鎮(zhèn)化進(jìn)行一階差分處理，將一階滯后項(xiàng)作為工具變量進(jìn)行兩階段最小二乘回歸（表4）。IV-2SLS的回歸結(jié)果與雙向固定效應(yīng)模型回歸結(jié)果基本一致。尤其是引入控制變量后，運(yùn)用工具變量法得到的城鎮(zhèn)化水平估計(jì)系數(shù)顯著上升，表明忽視內(nèi)生性問(wèn)題將出現(xiàn)明顯的估計(jì)偏差。因此，本文以工具變量-兩階段最小二乘法的回歸結(jié)果作為基準(zhǔn)回歸結(jié)果對(duì)其進(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)化對(duì)糧食生態(tài)效率的影響。表4中模型7結(jié)果表明，人口城鎮(zhèn)化對(duì)糧食生態(tài)效率的影響彈性為1.071。我們可從人口城鎮(zhèn)化對(duì)農(nóng)戶(hù)生產(chǎn)決策的影響對(duì)此做出解釋。人口城鎮(zhèn)化導(dǎo)致農(nóng)村勞動(dòng)力減少，勞動(dòng)力市場(chǎng)供求關(guān)系的改變使得農(nóng)業(yè)雇工價(jià)格居高不下，糧食生產(chǎn)面臨勞動(dòng)力和雇工成本的雙重壓力[20]，理性農(nóng)戶(hù)將從兩個(gè)環(huán)節(jié)對(duì)此壓力做出響應(yīng)：首先調(diào)整種植結(jié)構(gòu)，進(jìn)而通過(guò)替代性生產(chǎn)資料投入以最大限度地減少勞動(dòng)投入。當(dāng)雇工成本過(guò)高時(shí)，農(nóng)戶(hù)會(huì)更傾向于種植耗費(fèi)人力勞動(dòng)更少的糧食作物，減少經(jīng)濟(jì)作物種植[21]，同時(shí)提高各生產(chǎn)環(huán)節(jié)的機(jī)械化水平以減少人力勞動(dòng)。例如，在整地、播種、收割等環(huán)節(jié)提高機(jī)械作業(yè)比例，在田間管理中大量施用農(nóng)藥以減少人工除草成本[22]、提高化肥施用強(qiáng)度以降低施肥頻率，其結(jié)果一方面導(dǎo)致了糧食生產(chǎn)碳排放和面源污染程度增加[23]，另一方面還可能因不合理的施肥施藥降低糧食產(chǎn)量，降低糧食生產(chǎn)生態(tài)效率。此外，城鎮(zhèn)化的就業(yè)篩選機(jī)制只給人力資本水平較高者流入城市的機(jī)會(huì)[24]，年齡較大、受教育水平較低者只能繼續(xù)留守農(nóng)村，這也將制約綠色生產(chǎn)技術(shù)在糧食生產(chǎn)中的應(yīng)用，阻礙糧食生態(tài)效率的提高。

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

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

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

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

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

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

2）非主產(chǎn)區(qū)。糧食非主產(chǎn)區(qū)的回歸結(jié)果與全國(guó)層面基本一致。在主銷(xiāo)區(qū)，人口城鎮(zhèn)化對(duì)糧食生態(tài)效率的負(fù)向作用更強(qiáng)。從現(xiàn)實(shí)來(lái)看，我國(guó)糧食主銷(xiāo)區(qū)多集中在東部經(jīng)濟(jì)較為發(fā)達(dá)的省份，其糧食播種面積和糧食產(chǎn)量均占到全國(guó)水平的10%左右。但2006年主銷(xiāo)區(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òng)力短缺導(dǎo)致相關(guān)生產(chǎn)要素配置效率低、面源污染和碳排放嚴(yán)重，對(duì)糧食生態(tài)效率的負(fù)向影響更為顯著。

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

3.3 穩(wěn)健性檢驗(yàn)

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

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

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

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

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

4.1 研究結(jié)論

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

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

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

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

4.2 對(duì)策建議

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

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

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

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