城鎮(zhèn)化影響化肥施用強度了嗎？

伍國勇 包瀅暉

摘 要：本文測度城鎮(zhèn)化對化肥施用強度的影響，為降低化肥施用強度和實現糧食生產綠色化、現代化提供指引。本文基于2001—2017年全國31省面板數據，采用熵值法對構建的城鎮(zhèn)化綜合指標進行賦權。構建面板VAR模型和面板Tobit模型分別對城鎮(zhèn)化與化肥施用強度進行脈沖響應與回歸分析。結果表明：（1）人口城鎮(zhèn)化對化肥施用強度影響顯著。其中，在全國層面、主產區(qū)、平衡區(qū)、中部和西部地區(qū)有正向影響，在主銷區(qū)和東部地區(qū)有負向影響。（2）土地城鎮(zhèn)化對化肥施用強度影響差異顯著。其中，在全國層面和西部地區(qū)影響在統計學意義上并不成立;在主銷區(qū)、平衡區(qū)和東部地區(qū)有正向影響;在主產區(qū)耕地面積和人均綠地面積對其有正向影響;在中部地區(qū)耕地面積和城市征用土地面積對其有負向影響。（3）經濟城鎮(zhèn)化對化肥施用強度影響差異顯著。其中，在全國層面、主產區(qū)和中部地區(qū)人均地區(qū)生產總值與二三產業(yè)值比重對其有正向影響，人均非農產值對其有負向影響;在主銷區(qū)、平衡區(qū)、東西部地區(qū)有正向影響。基于實證分析結果，提出實現糧食施肥減量投入和糧食生產綠色化與農業(yè)現代化的建議。

關鍵詞：化肥施用強度;城鎮(zhèn)化;糧食生產;脈沖響應

中圖分類號：F326.11;S147.2

文獻標識碼：A

文章編號：1000-5099（2022）02-0049-12

化肥作為糧食生產要素投入的重要組成部分，對糧食增產發(fā)揮了至關重要的作用，但隨著化肥施用強度的增大[1]，使得糧食生產施肥過量[2-4]，對糧食生產和農業(yè)生態(tài)環(huán)境造成負面影響。一方面，會破壞土壤有機元素平衡結構，導致土壤存肥能力下降，使得化肥增產邊際效用下降[5]。另一方面，導致土壤酸化與水質下降，嚴重污染了生態(tài)環(huán)境[6]。隨著城市化的快速推進，大量農業(yè)勞動力非農化成為趨勢，城市用地規(guī)模擴大引起的耕地面積減少，經濟發(fā)展帶來的收入增加與工業(yè)反哺農業(yè)等都影響了糧食生產的化肥投入。那么，城鎮(zhèn)化快速發(fā)展對化肥施用強度的具體影響是提高了還是降低了？在糧食功能區(qū)與地理分區(qū)的影響是否存在差異？因此，深入探討二者關系既回答了上述問題，也有利于優(yōu)化糧食生產要素投入結構，為實現化肥減量投入，達到糧食生產綠色化與現代化相統一有重要意義。

當前，學者們對城鎮(zhèn)化與化肥施用強度的影響關系研究主要體現在如下層面：一是勞動力結構對化肥施用的影響。顏姜慧和朱舜[7]、劉穎等[8]認為，城市群發(fā)展吸引了眾多農村人口向城市聚集，實現了農業(yè)勞動力兼業(yè)化與非農化。進一步，劉華軍和劉傳明[9]、向晶和鐘甫寧[10]研究指出，農村青壯年勞動力的配置非農化使得農村農業(yè)勞動力老齡化速度加快，對農業(yè)有序生產有消極影響。李光泗和朱麗莉[11]、史常亮等[12]認為，大量農村人口外流，導致農業(yè)勞動力對糧食生產的邊際產出效應開始上升，農業(yè)勞動力不足使化肥、機械等要素開始替代勞動。胡浩和楊泳冰[13]、張艦等[14]、欒江[15]進一步研究指出，隨著城市群的發(fā)展，勞動力市場交易機制逐漸完善，引起農業(yè)勞動力用工成本上升，且其增長速度遠遠高于化肥價格的增長速度，這使得農戶愿意提高化肥施用強度以替代必要的勞動。張聰穎和霍學喜[16]則認為，農村勞動力轉移導致的農業(yè)生產勞動力不足會促使農戶采取測土配方施肥技術，從而減少了農業(yè)生產化肥施用量。不僅如此，張軍偉等[17]研究指出，糧食生產勞動力投入數量也會對化肥施用存在抑制效應。二是耕地面積變化對化肥施用的影響。潘丹[18]認為，城鎮(zhèn)化擴大了城市用地規(guī)模，導致耕地面積減少，對化肥消費有重要影響。諸培新等[19]、紀龍等[20]研究指出，農地經營規(guī)模越大、土地集中程度越高，越有利于促使農戶減少化肥投入。而黎孔清和馬豆豆[21]則認為，地塊破碎程度越高，農地經營越不集中，越不利于減少化肥投入。三是經濟發(fā)展對化肥施用的影響。城鎮(zhèn)化發(fā)展促進三產結構調整、經濟發(fā)展、技術進步和收入的增加，這些因素也會對化肥施用產生顯著影響。潘丹[18]研究指出，三產結構的調整將提高化肥施用強度。張軍偉等[17]研究發(fā)現，地區(qū)經濟發(fā)展水平也顯著促進了化肥投入的增加。而楊建輝[23]則認為，農業(yè)經濟發(fā)展水平與農業(yè)化學品投入存在絕對脫鉤關系，即農業(yè)經濟發(fā)展對其沒有影響。這與郭相平等[24]研究結論相反，他認為隨著經濟發(fā)展可以減少農作物化肥施用量，進而減少農作物水足跡。

綜上所述，雖然相關研究較為豐富，但城鎮(zhèn)化對化肥施用強度影響的文獻較少，且多數文獻的研究視角較為單一，也未指出城鎮(zhèn)化對化肥施用強度的影響存在區(qū)域差異。因此，本文構建城鎮(zhèn)化水平（人口、土地、經濟城鎮(zhèn)化）對糧食生產化肥施用強度影響的綜合分析框架，基于2001—2017年全國31省糧食生產的面板數據，采用面板VAR模型和Tobit模型對城鎮(zhèn)化與化肥施用強度進行脈沖響應分析和回歸分析。以期為化肥減量投入，實現糧食生產現代化、綠色化相統一和保障糧食安全提供參考。

一、研究方法與數據來源

（一）研究方法

1.熵值法。其基本原理是基于各觀測指標的差異度，合理確定其權重。賦權步驟如下：

（1）數據標準處理。若i個對象有j個觀測指標Zj，對其標準處理后有（假定標準化后指標值用Xij表示）：

Xij=Zij-minZijmaxZij-minZij（1）

（2）信息熵計算。對于i個對象有j個觀測指標Zj，則其信息熵為：

βij=Xij∑ni=1XijIj=-lnn-1∑ni=1βijlnβij（2）

式（2）中，I為信息熵。當且僅當β=0時，有βijlnβij收斂于0值。

（3）權重確定。假定ωj為各觀測指標權重，由式（2）可得：

ωj=1-Ij∑1-Ij（3）

2.面板VAR模型。相較于VAR模型，該模型引入面板結構，注重個體間的差異，動態(tài)估計參數方程。模型如下：

Hit=Hi，t-1B1+Hi，t-2B2+...+Hi，t-k+1Bk-1+Hi，t-kBk+RitD+ui+εit（4）

Ehitεit=Eritεit=Euiεit=0（5）

式（4）中，i為省份，t為年份，H為n個內生變量組成的行向量，R為m個外生變量組成的行向量，B、D為待估參數，ui為個體效應，εit為誤差項。對于面板VAR模型，采用動態(tài)GMM估計，使結果更為穩(wěn)健。

3.Tobit模型。該模型運用MLE估計的結果更穩(wěn)健，而傳統OLS估計忽視了因變量的約束條件，難以避免估計偏差。故采用如下Tobit模型：

yit=α0+x′itγ+ui+εit，若α0+x′itγ+ui+εit>00，若α0+x′itγ+ui+εit≤0 （6）

（二）數據來源

綜合考慮，截取2001—2017年作為研究時段。文中數據包括：《中國統計年鑒（2002—2018）》《中國人口和就業(yè)統計年鑒（2002—2018）》《中國農村統計年鑒（2002—2018）》《新中國六十年統計資料匯編》，部分數據來源于《中國城市統計年鑒（2002—2017）》和歷年各省統計年鑒及統計公報等。文中缺失數據采用歷年增長率進行補充，為消除各變量統計單位差異，對數據標準化。

二、城鎮(zhèn)化綜合指標體系構建

（一）指標選取

城鎮(zhèn)化正在經歷前所未有的大浪潮，為經濟高質量發(fā)展注入了強勁動力，推動了三次產業(yè)結構調整和收入增加。但城鎮(zhèn)化發(fā)展吸引了大量農村青壯年向城市轉移，實現農業(yè)勞動力就業(yè)非農化，影響了農業(yè)勞動力生產要素投入;城市用地規(guī)模擴大，不可避免地會影響農用可耕地面積。因此，選取人口城鎮(zhèn)化、土地城鎮(zhèn)化、經濟城鎮(zhèn)化綜合衡量城鎮(zhèn)化水平。

1.人口城鎮(zhèn)化。城市經濟的快速發(fā)展，使得農業(yè)比較利益降低，出于理性決策，大量農村青壯年向城市聚集，使得城市人口比重上升，客觀上實現了農業(yè)勞動力非農化配置，但影響了農業(yè)生產要素投入結構平衡。對于農業(yè)現代化程度較低、經濟落后或偏遠地區(qū)，其糧食生產對勞動力和化肥等傳統要素投入依賴程度大，勞動力非農化使得糧食生產要素投入不足。一方面，為彌補勞動力投入不足帶來的不利影響，以保證基本糧食生產，農戶會在一定程度上增施化肥。另一方面，農戶會減少糧食生產或實現耕地用途非農化，也影響了化肥施用。

2.土地城鎮(zhèn)化。城市用地規(guī)模增加，會占用耕地，使得糧食生產耕地面積減少，進而減少化肥施用。同時，城市人口增加對糧食需求增大，在有限的耕地面積內進行糧食生產面臨重要挑戰(zhàn)。因此，農業(yè)現代化程度較高的發(fā)達地區(qū)會調整化肥等要素投入結構以保證糧食產量，而對于農業(yè)現代化程度較低的偏遠落后地區(qū)，由于難以優(yōu)化糧食生產要素投入結構，且受糧食生產要素投入慣性影響，農戶可能增加化肥投入，影響了化肥施用。

3.經濟城鎮(zhèn)化。城鎮(zhèn)化發(fā)展優(yōu)化了產業(yè)結構，提高了人均可支配收入水平。一方面，人均可支配收入增加，有利于優(yōu)化消費結構，提升消費水平，對糧食生產要素消費結構產生了影響。另一方面，產業(yè)結構升級帶來的是二三產業(yè)產值比重增加，有利于實現工業(yè)反哺農業(yè)、城市支持農村的戰(zhàn)略目標，使傳統農業(yè)向現代農業(yè)轉變，客觀上也影響了糧食化肥施用。

（二）指標賦權

基于熵值法，我們對城鎮(zhèn)化各項指標進行賦權（表1）。總體而言，土地城鎮(zhèn)化（51.8%）>經濟城鎮(zhèn)化（33.31%）>人口城鎮(zhèn)化（14.89%）。具體來看，城市征用土地面積占比最大（32.5%），人均非農產值其次（17.66%），耕地面積和人均地區(qū)生產總值占比緊隨其后（15.5%），二三產業(yè)產值比占比最低（0.16%）。

三、城鎮(zhèn)化對化肥施用強度分析

（一）城鎮(zhèn)化對化肥施用強度的脈沖響應分析

基于赤池信息準則（Akaike information? criterion，AIC）、貝葉斯信息準則（Bayesian information criterion，BIC）、HQ準則（hannan-quinn information criterion，HQIC）等信息準則，綜合確定最優(yōu)滯后期數，并將研究時段劃分為17期，對面板VAR模型進行了脈沖響應分析（圖1）。

就全國而言，城鎮(zhèn)化對糧食生產化肥施用強度脈沖響應總體正向顯著。其中，在城鎮(zhèn)化開始沖擊后，化肥施用強度顯正響應;在第1期至第3期城鎮(zhèn)化對化肥施用強度的顯正響應水平下降;第3期開始變得不顯著;第3期之后又轉變?yōu)檎蝻@著響應;在第10期之后又開始緩慢下降，并逐漸趨于不顯著。這表明，在研究時段的前幾年城鎮(zhèn)化綜合水平提高了化肥施用強度，但在2001至2003年城鎮(zhèn)化綜合水平對化肥施用強度的增幅效應開始下降，并在2003年下降為0，之后緩慢增長，最后在2010年開始下降。

在主產區(qū)，當城鎮(zhèn)化綜合水平開始沖擊后，化肥施用強度并沒有響應，而是隨著時間推移開始顯正響應;在第2期顯正響應最大;在第3期化肥施用強度沒有響應城鎮(zhèn)化的沖擊，之后又顯正響應;在第10期開始下降，逐漸趨于不顯著。這表明，在初期城鎮(zhèn)化綜合水平的提高對化肥施用強度沒有影響，隨著時間推移，提高了化肥施用強度;在2002年城鎮(zhèn)化的增幅效應達到最大;在2003年其增幅效應下降為0，之后緩慢增長;最后在2010年又開始下降，并逐漸趨于0。

在主銷區(qū)，當城鎮(zhèn)化綜合水平開始沖擊，化肥施用強度響應水平不斷減弱;在第1期減弱為0，于第2期顯著負向沖擊;于第4期沖擊為0，隨后又開始正向沖擊;在第9期沖擊力度減弱，逐漸趨于0。這表明，在城鎮(zhèn)化沖擊初期，城鎮(zhèn)化綜合水平的提高促進了化肥施用強度的提高;在2001至2004年，化肥施用強度下降，并在2002年其減幅效應達到最大;在2004年之后開始轉為增幅效應，但在2009年其增幅效應又開始減弱，并逐漸趨于0。

在平衡區(qū)，當城鎮(zhèn)化綜合水平開始沖擊，化肥施用強度響應水平不斷減弱;在第2期不響應，于第3期顯著負向沖擊;于第4期沖擊為0，隨后正向沖擊;在第7期顯著正向沖擊，隨后減弱并逐漸趨于0。這表明，在城鎮(zhèn)化沖擊初期，城鎮(zhèn)化綜合水平的提高促進了化肥施用強度的提高;在2002至2004年，化肥施用強度下降，并在2003年其減幅效應達到最大;在2004年之后開始轉為增幅效應，但在2007年其增幅效應開始減弱，并逐漸趨于0。

此外，地理區(qū)位優(yōu)勢造就的經濟發(fā)展水平對化肥施用強度的沖擊也會有差異。因此，基于上述三種信息準則，綜合確定最優(yōu)滯后期數，同樣將研究時段劃分為17期，對面板VAR模型繼續(xù)進行脈沖響應分析（圖2）。

就東部地區(qū)而言，當城鎮(zhèn)化綜合水平開始沖擊后，化肥施用強度響應水平有小幅增加，在第1期達到最大;隨后開始減弱，在第2至第3期化肥施用強度響應水平顯著為負;之后繼續(xù)上升，并在第5期出現第二個正向響應峰值;但隨后又開始減弱，逐漸趨于0。這表明，在沖擊初期，城鎮(zhèn)化水平提高能增大化肥施用強度，但在2002至2003年降低了化肥施用強度;在2005年又顯著增大了化肥施用強度;2005年之后化肥施用強度增幅效應開始減弱，并逐漸變得不顯著。

就中部地區(qū)而言，在第1至第17期內，當城鎮(zhèn)化開始沖擊后，化肥施用強度基本沒有響應。這說明在研究時段內，城鎮(zhèn)化水平的提高對化肥施用強度的影響不顯著或沒有影響。

就西部地區(qū)而言，當城鎮(zhèn)化開始沖擊后化肥施用強度的響應速度減弱，在第1至第2期響應速度基本為0，在第3期響應速度顯著為負，在第5期及之后其響應速度顯著為正。這說明在研究時段的前幾年隨著城鎮(zhèn)化水平不斷提高，推動了化肥施用強度的增加，但增幅效應有所下降;在2003年其對化肥施用強度的影響顯著為負;在2005年之后又轉變?yōu)轱@著正向影響。

（二）城鎮(zhèn)化對化肥施用強度的回歸分析

1.Tobit模型構建

基于城鎮(zhèn)化綜合評價指標體系，構建如下面板Tobit模型：

Pit=0.176 1P1it+0.149 9P2it+0.674 0P3it（7）

Lit=0.627 5L1it+0.298 1L2it+0.074 4L3it（8）

Eit=0.464 9E1it+0.530 3E2it+0.004 9E3it（9）

URit=0.148 9Pit+0.518 0Lit+0.333 1Eit（10）

Fit=β0+βitURit+ui+εit（11）

式中，i為省份;t為年份;P為人口城鎮(zhèn)化綜合指數，其中，P1為城鎮(zhèn)人口比重，P2為非農就業(yè)人口比重，P3為城市人口密度;L為土地城鎮(zhèn)化綜合指數，其中，L1為城市征用土地面積，L2為耕地面積，L3為人均公共綠地面積;E為經濟城鎮(zhèn)化綜合指數，其中，E1為人均地區(qū)生產總值，E2為人均非農產值，E3為二三產業(yè)產值比重;UR為城鎮(zhèn)化綜合水平;F為化肥施用強度（用單位面積糧食生產化肥施用量表示）;β為常系數;ui為個體效應;εit為誤差項。各變量描述性統計如下（表2）。

2.單位根檢驗與協整檢驗

為保證數據回歸結果真實，先采用HT與IPS檢驗方法進行面板單位根檢驗。由于部分變量在零階和一階處不能拒絕單位根存在的假設，故進行二階單位根測試。結果表明，各變量存在二階單整（表3）。

采用Pedroni協整檢驗以驗證城鎮(zhèn)化與化肥施用強度的長期均衡關系，在協整檢驗中必須保證各變量同階單整。同時，由于變量較多，將所有變量分為三組進行檢驗（第一組：F與P1、P2、P3、L1、L2;第二組：F與P2、P3、L1、L2、L3、E1;第三組F與P3、L1、L2、L3、E1、E2、E3），檢驗結果均強烈拒絕沒有協整關系的原假設，表明城鎮(zhèn)化與化肥施用強度存在長期均衡關系（表4）。

3.Tobit回歸分析

首先進行OLS回歸作為Tobit回歸結果的對照（表5、表7）。關于Tobit模型是否采取固定效應或隨機效應模型回歸，通過似然比LR檢驗（表6、表8），結果表明強烈拒絕固定效應假設，認為有個體效應存在，所以采取包含個體的隨機效應面板回歸。不同城鎮(zhèn)化水平在三大糧食產區(qū)對化肥施用強度也會有差異，結果如表6所示。

（1）全國層面。人口城鎮(zhèn)化對化肥施用強度有顯著正向影響，即，隨著城鎮(zhèn)人口比重上升會增加化肥施用強度。這說明城市化的快速發(fā)展吸引了大量農業(yè)勞動力向城市轉移，農村勞動力就業(yè)非農化，導致糧食生產勞動力要素投入不足。對于全國農業(yè)現代化水平總體偏低的局面，且受傳統農業(yè)生產要素投入慣性影響，農戶會通過增加化肥投入以盡可能減少由于勞動力非農轉移對糧食生產帶來的不利影響。土地城鎮(zhèn)化對化肥施用強度沒有影響，說明城市用地規(guī)模增加對糧食生產的耕地投入影響不顯著，對化肥施用沒有產生影響。經濟城鎮(zhèn)化對化肥施用強度影響顯著。其中，人均地區(qū)生產總值與二三產業(yè)產值比重對其有顯正影響，人均非農產值對化肥施用強度有顯著負向影響。前者說明，經濟發(fā)展帶來的人均收入增加在理性決策下可能被用于非糧領域，并沒有改善糧食生產要素投入結構;二三產業(yè)產值比重上升帶來的工業(yè)反哺農業(yè)的激勵效果不明顯，使得農戶仍采用經驗性的施肥習慣，農業(yè)現代化之路任重道遠;城市化快速發(fā)展使得城市人口增加對糧食需求增加，在有限耕地面積約束條件下，糧食生產任務重大。因此，農戶會通過大量增加化肥、農藥等要素投入刺激糧食生產，以保障糧食供需平衡。后者說明，農業(yè)比較利益較低，且二三產業(yè)產值比重增加，農戶出于理性決策，會實現勞動力非農轉移或兼業(yè)化，從而減少或放棄一部分糧食生產，減緩了化肥施用強度。

（2）主產區(qū)層面。人口城鎮(zhèn)化對化肥施用強度有顯著正向影響。這說明在主產區(qū)勞動力要素配置對糧食生產至關重要，隨著大量農業(yè)人口非農轉移，影響了要素投入平衡結構，農戶為確保糧食生產會在一定程度上增施化肥彌補勞動力投入不足的影響。土地城鎮(zhèn)化對化肥施用強度影響顯著。其中，人均公共綠地面積對其有顯著正向影響，而耕地面積對其也有顯著正向影響。前者說明主產區(qū)作為全國糧食安全保障重心，人均公共綠地面積擴大帶來的耕地面積減少使其在短時間內難以優(yōu)化要素投入結構，為了保障糧食安全，客觀上會增施化肥促進生產。后者說明城鎮(zhèn)化水平的提高引起城市人口快速增加，人均公共綠地面積減少。為促進城鎮(zhèn)均衡發(fā)展，會擴大城市用地規(guī)模，導致耕地面積與糧食生產用地減少，在一定程度上減少了化肥施用。經濟城鎮(zhèn)化對化肥施用強度影響顯著。其中，人均地區(qū)生產總值與全國相比其正向影響減弱。這說明隨著收入增加主產區(qū)農民會改善糧食生產條件，增加產出效率更高的要素投入，客觀上替代了化肥、農藥等部分傳統要素，使得化肥施用強度減弱。人均非農產值負向影響減弱，這說明非農領域的豐厚收益吸引了大批農業(yè)勞動力。然而，主產區(qū)較豐厚的種糧補貼在一定程度上減緩勞動力的外流規(guī)模，對糧食生產要素投入減少有一定抑制作用，使得化肥減施速度放緩。二三產業(yè)產值比重正向影響顯著，與全國層面一樣，非農產值比重越低，化肥施用強度越大。

（3）主銷區(qū)層面。人口城鎮(zhèn)化對化肥施用強度有較為顯著的負向影響，這表明農業(yè)勞動力非農轉移，會引起糧食生產勞動力投入不足，使得農民減少糧食生產，降低了化肥施用強度。同時，主銷區(qū)大多處于經濟發(fā)達地區(qū)，人均購買力較高。因此，農民也會減少糧食生產，也降低了化肥施用強度。土地城鎮(zhèn)化對化肥施用強度影響較為顯著。其中，耕地面積對其有較為顯著的負向影響，人均公共綠地面積對化肥施用強度正向影響顯著。由于主銷區(qū)經濟發(fā)展程度高，土地資源寶貴，耕地面積減少影響了基本糧食生產。隨著城市人口不斷增加，對糧食需求不斷擴大，促使農民增施化肥增加糧食產量。經濟城鎮(zhèn)化對化肥施用強度有顯著正向影響。與二三產業(yè)相比農業(yè)比重下降，此時農戶可能會增施化肥以維護本地糧食生產。

（4）平衡區(qū)層面。人口城鎮(zhèn)化對化肥施用強度正向影響顯著，如前所述，農戶在一定程度上增施化肥保障糧食生產。在土地城鎮(zhèn)化中，耕地面積對化肥施用強度有較為顯著的負向影響。由于平衡區(qū)土地城鎮(zhèn)化減少了耕地面積，其糧食生產受阻。且受其農業(yè)現代化程度低和農業(yè)生產慣性的雙重影響，農戶可能會增加化肥施用，以保障糧食正常生產。經濟城鎮(zhèn)化對化肥施用強度正向影響顯著，如前所述，農戶也會增加化肥投入。

地理區(qū)位優(yōu)勢也會影響化肥施用強度，綜合地區(qū)經濟發(fā)展水平與糧食生產特征進行回歸，分析不同地理區(qū)位對化肥施用強度的影響，結果如表8所示。

（1）東部地區(qū)。人口城鎮(zhèn)化對化肥施用強度有顯著的負向影響?？赡茉颍阂皇菛|部地區(qū)靠近沿海，經濟發(fā)展程度高，勞動力以非農就業(yè)為主，使得農業(yè)生產勞動力少，在一定程度上減少了糧食生產。二是人均收入高，購買力較強，也減少了糧食生產。三是東部地區(qū)大多處于主銷區(qū)，糧食以外購為主。因此，綜合原因促使東部地區(qū)人口城鎮(zhèn)化水平越高而化肥施用強度越低。土地城鎮(zhèn)化對化肥施用強度影響顯著。其中，耕地面積對其有較顯著的負向影響，表明城市規(guī)模的擴大，使得耕地減少，糧食生產面臨重要挑戰(zhàn)，此時農戶會增施化肥以求短時間穩(wěn)定糧食生產。同樣，人均公共綠地面積的增加也使得化肥投入增加。經濟城鎮(zhèn)化對化肥施用強度正向影響顯著，如前所述，農戶也會增施化肥，以保證糧食基本生產。

（2）中部地區(qū)。人口城鎮(zhèn)化對化肥施用強度有正影響（該結果得到統計數據支持），原因是中部地區(qū)大多位于糧食主產區(qū)，農業(yè)勞動力資源寶貴，如前所述，其大量非農轉移會促使農戶使用化肥替代，客觀上增加了化肥施用強度。土地城鎮(zhèn)化對化肥施用強度有顯著影響。其中，城市征用土地面積對其有負向影響（該結果得到統計數據支持），耕地面積對其有負向影響（該結果得到統計數據支持）。一方面，城市用地規(guī)模擴大，導致耕地面積減少，使得糧食生產下降，減少了化肥施用。另一方面，城市人口大量增加對糧食需求增大，促使農民增施化肥提高糧食產出效率。經濟城鎮(zhèn)化對化肥施用強度影響顯著，其中人均地區(qū)生產總值和二三產業(yè)產值比重對其有正向影響（該結果得到統計數據支持），人均非農產值對其有負向影響（該結果得到統計數據支持）。如前所述，前者引起糧食供需失衡和農業(yè)產值比重下降，會推動化肥投入的增加。后者引起生產要素的非農配置，從而減少化肥投入。

（3）西部地區(qū)。人口城鎮(zhèn)化對化肥施用強度正向影響顯著，原因是西部地區(qū)經濟發(fā)達程度和農業(yè)現代化程度低，農業(yè)勞動力大量非農化，使得該地區(qū)糧食生產勞動力投入不足，客觀上推動了農戶利用化肥替代勞動力，引起化肥投入增大。土地城鎮(zhèn)化對化肥施用強度影響不顯著，表明城市用地規(guī)模增加沒有影響到糧食生產，對其要素投入影響不顯著。經濟城鎮(zhèn)化對化肥施用強度有正影響（該結果得到統計數據支持），與糧食產區(qū)和東中部地區(qū)一樣，工業(yè)反哺農業(yè)力度有待提高，且資本的非農配置，導致化肥施用增加。

四、結論與啟示

本文基于省際面板數據，構建包含人口城鎮(zhèn)化、土地城鎮(zhèn)化、經濟城鎮(zhèn)化的綜合分析框架，采用熵值法對城鎮(zhèn)化綜合指標賦權，構建面板VAR模型和面板Tobit模型對城鎮(zhèn)化與化肥施用強度進行實證分析。結果表明：1.人口城鎮(zhèn)化對化肥施用強度影響顯著。就全國層面而言，人口城鎮(zhèn)化對化肥施用強度有正向影響（該結果得到統計數據支持）。就糧食功能區(qū)而言，在主產區(qū)和平衡區(qū)人口城鎮(zhèn)化對化肥施用強度有正向影響（該結果得到統計數據支持），在主銷區(qū)人口城鎮(zhèn)化對化肥施用強度有負向影響（該結果得到統計數據支持）。就地理區(qū)位而言，在中部和西部地區(qū)人口城鎮(zhèn)化對化肥施用強度正向影響，在東部地區(qū)人口城鎮(zhèn)化對化肥施用強度有負向影響（該結果得到統計數據支持）。2.土地城鎮(zhèn)化對化肥施用強度影響存在顯著差異。就全國層面而言，土地城鎮(zhèn)化對化肥施用強度影響未取得統計數據支持。就糧食功能區(qū)而言，土地城鎮(zhèn)化在主銷區(qū)和平衡區(qū)均對其有正向影響，但在主產區(qū)存在顯著差異（該結果取得統計數據支持）。具體表現為：耕地面積減少降低了化肥施用強度，人均公共綠地面積增加增大了化肥施用強度。就地理分區(qū)而言，在東部地區(qū)土地城鎮(zhèn)化對化肥施用強度正向影響顯著;在中部地區(qū)土地城鎮(zhèn)化對其影響存在顯著差異。具體表現為：城市征用土地面積引起化肥施用強度降低，耕地面積減少引起化肥施用強度增大;在西部地區(qū)土地城鎮(zhèn)化對化肥施用強度影響不顯著（該結果未取得統計數據的支持）。3.經濟城鎮(zhèn)化對化肥施用強度影響差異顯著。就全國層面而言，經濟城鎮(zhèn)化對化肥施用強度影響顯著，其中人均地區(qū)生產總值與二三產業(yè)產值比重對化肥施用強度有正向影響，人均非農產值對化肥施用強度有負向影響（該結果得到統計數據支持）。就糧食功能區(qū)而言，在主產區(qū)經濟城鎮(zhèn)化對化肥施用強度影響顯著（該結果取得統計數據的支持）。表現為：人均地區(qū)生產總值的正向影響減弱，人均非農產值負向影響減弱;在主銷區(qū)與平衡區(qū)經濟城鎮(zhèn)化均對化肥施用強度有正向影響（該結果得到統計數據支持）。就地理區(qū)位而言，在東西部地區(qū)經濟城鎮(zhèn)化均對化肥施用強度正向影響顯著，但在中部地區(qū)存在差異。表現為：人均地區(qū)生產總值和二三產業(yè)產值比重對其有正向影響，人均非農產值對其有負向影響（該結果得到統計數據支持）。

綜上所述，為了降低化肥施用強度，實現糧食生產綠色化和農業(yè)現代化，可從以下幾個方面入手：一是合理引導農業(yè)勞動力非農轉移，對勞動力資源配置進行科學規(guī)劃，推動實現農業(yè)勞動力就地就業(yè)，防止農業(yè)勞動力過度轉移引起化肥替代勞動力現象的發(fā)生。二是適度推進城市化進程，在規(guī)劃城市用地應優(yōu)先保障糧食生產用地，防止過多占用田地導致糧食播種面積減少，引發(fā)糧食供需矛盾而使得化肥施用增加。三是加大力度實施工業(yè)反哺農業(yè)，城市支持農村的戰(zhàn)略。加大政府支農財政力度，大力實施科技興農、藏糧于地、藏糧于技的重大戰(zhàn)略，研發(fā)優(yōu)質、高產、綠色的糧食作物新品種，推動傳統農業(yè)向現代農業(yè)轉變。

參考文獻：

[1]欒江，仇煥廣，井月，等.我國化肥施用量持續(xù)增長的原因分解及趨勢預測[J].自然資源學報，2013（11）：1869-1878.

[2]仇煥廣，欒昊，李瑾，等.風險規(guī)避對農戶化肥過量施用行為的影響[J].中國農村經濟，2014（3）：85-96.

[3]紀月清，張惠，陸五一，等.差異化、信息不完全與農戶化肥過量施用[J].農業(yè)技術經濟，2016（2）：14-22.

[4]孔凡斌，郭巧苓，潘丹.中國糧食作物的過量施肥程度評價及時空分異[J].經濟地理，2018（10）：201-210.

[5]陸文聰，劉聰.化肥污染對糧食作物生產的環(huán)境懲罰效應[J].中國環(huán)境科學，2017（5）：1988-1994.

[6]周亮，徐建剛，蔡北溟，等.淮河流域糧食生產與化肥消費時空變化及對水環(huán)境影響[J].自然資源學報，2014（6）：1053-1064.

[7]顏姜慧，朱舜.農村人口轉移趨勢及空間指向研究[J].中國人口·資源與環(huán)境，2017（5）：146-152.

[8]劉穎，鄧偉，宋雪茜，等.基于綜合城鎮(zhèn)化視角的省際人口遷移格局空間分析[J].地理科學，2017（8）：1151-1158.

[9]劉華軍，劉傳明.城鎮(zhèn)化與農村人口老齡化的雙向反饋效應：基于中國省際面板數據聯立方程組的經驗估計[J].農業(yè)經濟問題，2016（1）：45-52.

[10]向晶，鐘甫寧.農村人口轉移、工業(yè)化和城鎮(zhèn)化[J].農業(yè)經濟問題，2018（12）：51-56.

[11]李光泗，朱麗莉.農村勞動力流動背景下我國糧食生產技術變動分析[J].中國科技論壇，2014（7）：143-148.

[12]史常亮，朱俊峰，欒江.我國小麥化肥投入效率及其影響因素分析：基于全國15個小麥主產省的實證[J].農業(yè)技術經濟，2015（11）：69-78.

[13]胡浩，楊泳冰.要素替代視角下農戶化肥施用研究：基于全國農村固定觀察點農戶數據[J].農業(yè)技術經濟，2015（3）：84-91.

[14]張艦，亞伯拉罕·艾賓斯坦，瑪格麗特·麥克米倫，等.農村勞動力轉移、化肥過度使用與環(huán)境污染[J].經濟社會體制比較，2017（3）：149-160.

[15]欒江.農業(yè)勞動力轉移與化肥施用存在要素替代關系嗎？：來自我國糧食主要種植省份的經驗證據[J].西部論壇，2017（4）：12-21.

[16]張聰穎，霍學喜.勞動力轉移對農戶測土配方施肥技術選擇的影響[J].華中農業(yè)大學學報（社會科學版），2018（3）：65-72.

[17]張軍偉，張錦華，吳方衛(wèi).糧食生產中化肥投入的影響因素研究：基于Durbin模型的分析[J].經濟地理，2018（11）：174-182.

[18]潘丹.中國化肥消費強度變化驅動效應時空差異與影響因素解析[J].經濟地理，2014（3）：121-126.

[19]諸培新，蘇敏，顏杰.轉入農地經營規(guī)模及穩(wěn)定性對農戶化肥投入的影響：以江蘇四縣（市）水稻生產為例[J].南京農業(yè)大學學報（社會科學版），2017（4）：85-94.

[20]紀龍，徐春春，李鳳博，等.農地經營對水稻化肥減量投入的影響[J].資源科學，2018（12）：2401-2413.

[21]黎孔清，馬豆豆.生態(tài)脆弱區(qū)農戶化肥減量投入行為及決策機制研究：以山西省4縣421戶農戶為例[J].南京農業(yè)大學學報（社會科學版），2018（5）：138-145.

[22]潘丹.中國化肥施用強度變動的因素分解分析[J].華南農業(yè)大學學報（社會科學版），2014（2）：24-31.

[23]楊建輝.農業(yè)化學投入與農業(yè)經濟增長脫鉤關系研究：基于華東6省1市數據[J].自然資源學報，2017（9）：1517-1527.

[24]郭相平，高爽，吳夢洋，等.中國農作物水足跡時空分布與影響因素分析[J].農業(yè)機械學報，2018（5）：295-302.

（責任編輯：楊 洋）

The Influence of Urbanization on the Intensity of Fertilizer Usage：

the Analysis Based on the Fact of Grain Production in 31

Provinces（Autonomous Ragions and Municipalities）

WU Guoyong1，2，BAO Yinghui3

（1.China Institute of Rural Revitalization in Karst Areas，Guizhou University，Guiyang，Guizhou，China，550025;2.Guizhou High-end Think Tank for Grassroots Social Governance，Guizhou University，Guiyang，Guizhou，China，550025;3.Subsidiary Unit of Groupe ESC Clermont-Graduate School of Management，Clermont-Ferrand，Pue-de-Dome，France，63000）

Abstract：

This paper aims to measure the impact of urbanization on fertilizer usage intensity so as to provide guidance for reducing fertilizer usage intensity to achieve green and modern grain production.Based on the panel data of 31 provinces in China from 2001 to 2017，this paper adopts the entropy method to weight the constructed comprehensive urbanization indicators，and constructs panel VAR model and panel Tobit model to analyze respectively the impulse response and regression of the relation between urbanization and fertilizer usage intensity.The results find that： firstly，population urbanization has a significant influence on fertilizer usage intensity. positive effects are across the whole country，main production areas，balance areas，central and western regions，while the negative effects can be seen in the main sales areas and eastern regions.Secondly，land urbanization has a significant impact on the intensity of fertilizer usage.The influence was not statistically tenable in the whole country and the western regions，while the main sales areas，balance areas and eastern regions are influenced in a positive way.Besides，the cultivated land area and per capita green land area in the main producing areas have a positive impact on the fertilizer usage intensity，while the cultivated land area，as well as the urban land requisition area in central China has negative influence on it.Thirdly，the influence of economic urbanization on fertilizer usage intensity differs significantly.The ratio of GDP per capita to the value of secondary and tertiary industries in the whole country，major producing areas and the central regions has a positive impact，while per capita non-agricultural value has a negative impact.It has positive influence in main sales area，balance areas，eastern and western areas.Based on the results of empirical analysis，some suggestions are put forward to realize the reduction of fertilizer input，the green grain production and agricultural modernization.

Key words：

fertilizer usage intensity; urbanization; food production; impulse response

3185500589256