化肥施用強度對中國糧食單產的影響分析
——基于省級面板數據的分位數回歸

朱滿德，李辛一，徐雪高

（1.貴州大學，a. 經濟學院，b. 中國西部發展能力研究中心，貴州 貴陽 550025；2. 江蘇省農業科學院農業經濟與發展研究所，江蘇 南京 210014）

化肥施用強度對中國糧食單產的影響分析
——基于省級面板數據的分位數回歸

朱滿德1，李辛一1，徐雪高2※

（1.貴州大學，a. 經濟學院，b. 中國西部發展能力研究中心，貴州 貴陽 550025；2. 江蘇省農業科學院農業經濟與發展研究所，江蘇 南京 210014）

在分析中國糧食作物化肥施用特征的基礎上，運用分位數回歸模型探究了1994年以來化肥施用強度對早秈稻、中秈稻、晚秈稻、粳稻、小麥、玉米6種糧食單產水平的影響。研究表明：1）中國糧食作物的化肥施用強度增長具有明顯的階段性特征；6種糧食作物的化肥施用強度持續增加，且均超出國際安全施用標準；小麥與玉米的化肥施用強度顯著高于水稻；部分品種的化肥施用強度已出現優勢產區和非優勢產區、東中西部的區域差異。2）糧食作物的化肥施用強度與其單產水平具有顯著的關聯效應。3）晚秈稻、粳稻和玉米的化肥施用強度對其單產仍存在正向作用，但已呈現邊際效用遞減特征；小麥、早秈稻和中秈稻的這種影響不再顯著，已有過度施用特征。今后在堅持化肥施用零增長前提下，注重科學施用化肥、有效提高化肥利用率，在保障糧食安全與促進農業可持續發展之間達成平衡，是堅持中國特色農業現代化道路的必然要求。

化肥施用強度；糧食單產；分位數回歸；區域分異；邊際效用；化肥施用零增長

Abstract：Based on provincial panel data, and applying the quantile regression model, this paper analyzed the effects of fertilizer application on grain yield in China. Varieties studied in this research include early indica-rice, middle indicarice, late indica-rice, japonica-rice, wheat and corn. Results show that China’s agricultural fertilizer application has been growing continuously and has exceeded the international safety standard. The intensities of fertilizer on wheat and corn were higher than that of rice. There exists intensity difference between advantage producing areas and disadvantage producing areas and it also exists between eastern China and western China. Descriptive statistical analysis indicates that the strength of fertilizer and grain yield have the same trend. Chemical fertilizers still have positive effects on grain yields while showing signifcant declining marginal return. In the future, the Government should encourage zero growth of fertilizer application, pay more attention to improve fertilizer effcacy, and keep the balance between grain security and sustainable agricultural development.

Key words：application intensity of fertilizer; grain yield; quantile regression model; regional difference; marginal utility; zero growth of fertilizer

化肥是重要的生產資料，是糧食的“糧食”。改革開放30余年來，中國化肥施用總量與結構都有顯著變化，施用強度不斷增加。化肥施用為中國持續增加糧食產量、保障食物有效供給發揮了重要作用；但近年來已經出現盲目、過量施用等問題[1]。部分地區因化肥施用不盡合理，造成土壤有機質降低、營養失衡、酸化和板結[2-3]，也導致農業面源污染加劇、作物產量下滑[4-6]，并引起對中國農業可持續發展能力的擔心和憂慮。2015年中國政府出臺《到2020年化肥使用量零增長行動方案》[7]，針對當前化肥施用存在的問題提出了針對性的應對舉措，如推進精準施肥、調整化肥使用結構、改進施肥方式、有機肥替代化肥等，到2020年實現主要農作物化肥施用量“零增長”，明確了未來中國化肥施用的調整方向。進一步分析化肥施用強度對中國糧食產量的影響，對實現化肥“零增長”行動、推進農業及糧食生產方式轉型、確保谷物基本自給意義重大。

化肥施用主要通過作用于糧食單產的方式提高糧食總產量，屬“增產型”技術。有關文獻研究表明，中國化肥施用量快速增長，東部平原和東南沿海地區施用偏多，西部山地和丘陵地區施用相對不足；化肥施用強度不斷提高，單位面積化肥施用量已超越世界大多數發達國家，并高于國際公認的安全標準上限；部分地區、部分作物的化肥投入過量、氮磷鉀搭配不合理等問題日益突出[8-10]。進一步針對糧食的研究表明，化肥施用量對中國糧食產量有顯著正向影響，是播種面積、用工數量、其他物質費用、受災面積等因素中貢獻最大的一項[11]。但近年中國化肥施用量增長對糧食生產的增產效應開始不顯著，表現出邊際效用遞減特征，糧食作物的化肥施用量超過了經濟最優施用水平；其中水稻、小麥、玉米存在嚴重的氮肥過量施用問題，水稻化肥施用效率高于小麥和玉米；小麥化肥施用效率不高；玉米化肥施用效率上升較快[12-13]。

中國化肥施用總量快速增長是不爭的事實，但不同作物品種、不同區域的化肥施用強度則存在明顯差異。是否所有糧食品種都存在化肥施用邊際效用遞減或施用過量特征？本文將運用分位數回歸方法探討化肥施用對不同糧食作物單產的具體影響，探討究竟是增產正效應的邊際遞減，抑或是已呈現負效應，以期形成具體的判斷，從而就實現化肥零增長行動，推動糧食生產向環境友好型方式轉變，促進糧食生產的可持續發展提出政策建議。

1 研究方法與數據來源

1.1 分位數回歸模型構建

分位數回歸模型（Quantile Regression Models，QRM）由Koenker和Bassett在1978年提出的，是線性回歸模型的擴展。分位數回歸模型是基于因變量的條件分位數對自變量進行回歸，進一步得到“所有”分位數條件下的回歸模型，其最大優點是可以關注到全體因變量的分布，對研究自變量對因變量在特定分位點的邊際彈性會有所裨益[14]。為此，本文評估化肥施用強度對中國不同品種糧食單產的影響，使用分位數回歸模型可以探究特定分位（以及所有不同分位條件）的糧食單產水平下化肥施用強度是否存在不同程度的影響。

農業自然再生產和經濟再生產相互交織，它所具有的復雜性和系統性使得影響糧食單產的因素有很多，例如：不同作物乃至同一作物不同品種的生物性狀，耕地、化肥、勞動力等生產要素質量、投入強度及其組合，農業實用技術的應用推廣（如農作物良種、農機深耕作業、測土配方施肥），氣候條件、自然災害等。為探究化肥施用強度對中國早秈稻、中秈稻、晚秈稻、粳稻、小麥、玉米6個主糧品種單產水平的影響，根據農作物生產函數特點，在Cobb-Douglas（C-D）生產函數的基礎上進行拓展。模型所使用的被解釋變量為單位面積糧食平均產量（Y，kg/hm2）；本文所關注的解釋變量為化肥施用強度，用單位面積化肥施用量（FERT，kg/hm2）表示；同時使用單位面積機械作業費（MACH，元/hm2）、單位面積用工數量（WORK，個/hm2）、農戶銷售價格（PRIC，元/50 kg，模型采用幼稚性預期假設，對其進行滯后一期處理）、災害成災率（DISA，%），年降水量（RAIN，mm）、年平均氣溫（TEMP，℃）作為控制性變量[15-16]，構建分位數回歸模型。本文也嘗試了將單位面積種子用量或種子費引入模型進行估計，但模擬效果不佳且模型穩健性受到影響，故未納入最后計量模型。

對擴展的柯布-道格拉斯生產函數取對數處理，形成如下回歸方程：

式中：α為截距項，βn（n=1，2，3，…，7）表示各指標的參數估計值，上標(q)表示特定分位數，下標i代表不同的省區，t代表不同的時間，ε表示隨機誤差項。

一般情況下，Stata軟件可使用“bsqreg”命令對上式進行分位數回歸。但該命令只能適用于時間序列數據或截面數據，并不能直接對面板數據實施分位數估計。解決的方法是，先把省區作截面虛擬變量處理，然后再利用該命令進行分位數回歸運算。在此過程中，使用了自舉法（Bootstrap Method）進行估計，以減少偏誤，得到更加準確有效的估計。

1.2 數據來源

本文所使用的數據中，早秈稻、小麥、玉米等6種糧食的單位面積平均產量，單位面積的化肥施用量、機械作業費、用工數量，以及農戶銷售價格數據來自《建國以來全國主要農產品成本收益資料匯編（1953-1997年）》和《全國農產品成本收益資料匯編》（1999-2016年）。其中，對單位面積機械作業費、農戶銷售價格分別采用農業生產資料價格指數、消費者價格指數（相關數據來自《中國統計年鑒》）進行平減，以剔除價格因素影響。年降水量和年平均氣溫使用各省區的省會城市替代，數據來自《中國農村統計年鑒》和《中國氣象年鑒》；自然災害成災率由成災面積和受災面積估算得到，相關數據來自《中國統計年鑒》。

有關統計資料提供了中國1953年以來的糧食投入產出數據，但經比較分析發現，可能因統計口徑等問題，1978-1993年間，各糧食品種的化肥施用強度多次出現遠高于1994年乃至當前化肥施用強度的情形。例如，1987年全國早秈稻化肥施用量達885 kg/hm2，是2015年全國平均施用強度的2.6倍；1988年全國晚秈稻化肥施用885 kg/hm2，最高省區達1 260 kg/hm2，是2015年全國平均施用強度的3.8倍；1990年全國粳稻化肥施用1 155 kg/hm2，最高省區達1 650 kg/hm2，是2015年全國平均施用強度的4.6倍，這表明化肥施用數據存在“失真”的可能，因此本文選用的時間起點為1994年。

綜合考慮糧食種植區域特征及數據可得性，本文研究的早秈稻包括浙江、安徽、福建、江西、湖北、湖南、廣東、廣西、海南9省區，中秈稻包括江蘇、安徽、福建、河南、湖北、四川、貴州、陜西8省，晚秈稻包括浙江、安徽、福建、江西、湖北、湖南、廣東、廣西、海南9省區，粳稻包括河北、遼寧、吉林、黑龍江、江蘇、浙江、安徽、山東、河南、湖北、云南、寧夏12省區，小麥包括河北、山西、內蒙古、黑龍江、江蘇、安徽、山東、河南、湖北、四川、云南、陜西、甘肅、寧夏、新疆15省區，玉米包括河北、山西、內蒙古、遼寧、吉林、黑龍江、江蘇、安徽、山東、河南、湖北、廣西、四川、貴州、云南、陜西、甘肅、寧夏、新疆19省區，上述省區各糧食品種的產量之和都占到其總產量的95%以上，因而具有代表性。

2 化肥施用強度與糧食單產的描述性統計

2.1 糧食作物的化肥施用強度：區域及品種特征

1994-2015年，中國農用化肥施用量增長迅速，由3 318萬t增至6 023萬t，增長2 705萬t，增幅82%；其中糧食作物的化肥施用強度由253 kg/hm2增至362 kg/hm2（系化肥折純用量，下同），增長43%。分析顯示，糧食作物的化肥施用強度增長具有顯著的階段性波動特征，大致可分為：1）1994-2002年糧食作物化肥施用強度增長階段，期間累計增長61 kg/hm2，增幅24%；2）2002-2004年糧食作物化肥施用強度下降階段，此間下降了27 kg/hm2，降幅約9%；3）自2004年開始，糧食作物化肥施用強度持續穩步增加階段，累計增長26%，其中2010年以來增速開始放緩。

表1列出了1994年和2015年部分省區早秈稻、中秈稻、晚秈稻、粳稻、小麥、玉米6種糧食作物的化肥施用強度。不同糧食品種、不同地區的化肥施用強度差異明顯。

2.1.1 水稻 1994-2015年，粳稻的化肥施用強度高于秈稻。以2015年為例，粳稻化肥施用強度最高，達361 kg/hm2；早秈稻、晚秈稻次之，分別為334 kg/hm2、336 kg/hm2，中秈稻最低，為301 kg/hm2。早秈稻、晚秈稻化肥施用強度未表現出東部、中部、西部等之間的省區差距，而中秈稻化肥施用強度的省區差異則較為明顯，西部省份低于東中部省份。特別是安徽，2015年早秈稻、晚秈稻化肥施用強度分別達到388 kg/hm2和406 kg/hm2，顯著高于其他省份；中秈稻332 kg/hm2，僅次于湖北。粳稻化肥施用強度在優勢產區與非優勢產區，東部、中部、西部與東北地區都存在明顯差距。作為粳稻的優勢產區，東北地區粳稻化肥施用強度較低，如2015年黑龍江251 kg/hm2、吉林330 kg/hm2；而非優勢產區化肥施用強度高，如河南為521 kg/hm2、寧夏為429 kg/hm2。同時，東部省份高于西部和中部省份，西部和中部省份又高于東北省份，如河北542 kg/hm2、江蘇549 kg/hm2，是黑龍江的2倍以上。

2.1.2 小麥 1994年和2015年東部省份的化肥施用強度全面高于中部省份，中部省份則高于西部省份，其中內蒙古和新疆例外。2015年內蒙古小麥化肥施用強度為531 kg/hm2，為全國最高；新疆為461 kg/hm2，也高于大部分省份。作為小麥主產區或優勢產區的河北、江蘇、山東、安徽、河南，其化肥施用強度明顯高于其他省份以及全國平均水平，只有湖北例外，這一點與“水稻優勢產區化肥施用強度低”正好相反。但應注意的是，1994-2015年間，中部和西部省份的小麥化肥施用強度正快速增加，如山西增長了149%，內蒙古增長了255%。

2.1.3 玉米 1994年和2015年玉米的化肥施用強度表現為西北、東北省份高于中部、西南省份，這與中秈稻、粳稻、小麥“東部省份施肥強度高、西部和東北省份施肥強度低”明顯不同。其中，東北和華北兩大玉米優勢產區化肥施用強度低于西南和西北省份，但略高于中部省份。與1994年相比，2015年山西增長161%、安徽增長146%、河南增長102%、新疆增長139%、黑龍江增長111%，5省區玉米化肥施用強度翻了一番有余。

表1 中國各地區不同糧食作物的化肥施用強度（kg/hm2）Table 1 Chemical fertilizers intensities for different grain corps in China (kg/hm2)

從全國看，如果以國際公認的化肥安全施用標準225 kg/hm2為參考，追溯到1994年，中秈稻化肥施用強度尚未超過安全施用標準，但也已經接近上限水平，早秈稻、晚秈稻、粳稻、小麥、玉米5種主糧的化肥施用強度均已超過了安全施用標準；到2015年，6種糧食已經大幅超過安全施用標準，其中小麥最高，是安全施用標準的1.8倍，中秈稻最低，為安全施用標準的1.3倍。如果以省區為基本單位，1994年相當一部分省區的化肥施用強度還沒有超過這一標準，到2015年只有極個別非主產區還沒有超過安全施用標準，大部分省份則遠超這一標準。如果考慮到中國糧食主產區多以兩季農作物為主，單位面積耕地承受的化肥施用強度將會更大。

2.2 化肥施用強度與糧食單產的關聯分析

為判斷各糧食作物化肥施用強度對其單產水平的影響，通過整體的化肥施用強度與其平均單產、不同品種的化肥施用強度與其平均單產、不同品種的單位化肥施用所對應的邊際產量三個方面進行描述性統計分析，以期為后續實證分析提供邏輯支撐。

2.2.1 糧食作物化肥施用強度與糧食平均單產的關聯 1994-2015年中國糧食作物的化肥施用強度（系水稻、小麥和玉米的平均）與其平均單產水平呈同向變動，都具有階段性特征，且按變化特征可劃分的階段也大體相似。圖1顯示，2004年之前化肥施用強度呈先升后降趨勢，糧食單產水平基本穩定，略高于5 000 kg/hm2；2004年之后糧食作物的化肥施用強度與糧食單產水平都平穩增加，二者變化趨勢相近，保持同步變動。數據測算表明，二者相關系數為0.85，屬高度相關，表明糧食作物的化肥施用強度增強與糧食單產水平提高在統計上具有正向相關關系。

圖1 糧食平均單產與化肥施用強度的變化Fig. 1 Change of the grain yield and the intensity of fertilizer

2.2.2 不同品種化肥施用強度與其單產的關聯 數據顯示，早秈稻等6種糧食作物的化肥施用強度和其單產水平都有不同程度增長，且化肥施用強度增長越顯著的品種其單產水平提高也越顯著。具體而言，小麥與玉米的化肥施用強度都有顯著增長，1994-2015年分別由236 kg/hm2、233 kg/hm2增加到406 kg/hm2、365 kg/hm2，增幅為72%與57%；同期小麥、玉米單產的增幅也分別達到73%和33%。與其形成對比，4種水稻的化肥施用強度增長趨勢不甚明顯，早秈稻、中秈稻、晚秈稻、粳稻分別增長了30%、41%、18%和7%，同期4種水稻單產水平的增幅只有18%-26%，較之小麥、玉米單產提高幅度明顯偏低（表2）。數據所反映出的這種“化肥施用強度-糧食單產”同向聯動變化，意味著糧食作物的化肥施用強度對其單產水平或有明顯正向影響，需實證進一步檢驗。

表2 不同糧食品種化肥施用強度與單產變化（kg/hm2）Table 2 Changes of grain yields and fertilizer application intensities for different grain crops（kg/hm2）

2.2.3 不同品種的單位化肥施用量對應的糧食平均產量 圖2顯示，中秈稻單位化肥施用量對應的平均產量較高，每千克化肥所產出的中秈稻產量均保持在25 kg以上；接下來依次是早秈稻、晚秈稻、粳稻、玉米；小麥偏低，每千克化肥所產出的小麥產量只在15 kg左右。2004-2015年特別是2010年以來，6種糧食的單位化肥施用量對應的糧食平均產量變動較小，表明化肥施用增加所產生的邊際產量極小，這意味著當前條件下化肥施用強度對糧食單產增長或有邊際遞減特征，部分品種或已有過度施用特征。

3 化肥施用強度對中國糧食單產影響的實證：模擬結果及其解釋

運用Stata12.0軟件對1994-2013年早秈稻、中秈稻、晚秈稻、粳稻、小麥、玉米6種糧食作物單產水平的影響因素進行分位數回歸，結果如表3、圖3所示。囿于篇幅有限，表3只選擇了0.2、0.5和0.8分位點回歸系數進行匯報；圖3同時繪出了不同糧食作物化肥施用強度系數的分位數圖形（即反映不同分位點的彈性系數）。

圖2 單位化肥對應的平均糧食產量Fig. 2 Grain production per unit of fertilizer

表3 不同分位數下化肥施用強度等因素對糧食單產水平的影響Table 3 Chemical fertilizer intensity effect on grain yields at different quantiles

3.1 化肥施用強度對糧食單產的影響

本文關注的化肥施用強度，對晚秈稻、粳稻、小麥、玉米單產水平都有顯著的正向影響，且隨著分位點由低分位向高分位移動，這種正向影響出現一定程度的減弱趨勢；但它對早秈稻、中秈稻單產水平的影響則不顯著。

3.1.1 晚秈稻 晚秈稻的化肥施用強度對其單產水平影響在0.2分位點、0.5分位點和0.8分位點處分別通過了1%、5%和10%水平的顯著性檢驗，系數均為正值說明了過去20 a間化肥施用強度的提高確實有效提高了晚秈稻的單產水平；伴隨著分位點由低分位到高分位的移動，化肥施用強度的回歸系數（即彈性）由0.2分位點的0.198，下降到0.5分位點的0.142，進一步降低到0.8分位點的0.126，這也意味著化肥施用已經出現明顯的邊際效用遞減特征。

3.1.2 粳稻 與晚秈稻類似，粳稻化肥施用強度的回歸系數由0.2分位點的0.189，下降到0.5分位點的0.130，進一步降低到0.8分位點的0.098，化肥施用也是出現明顯的邊際效用遞減特征。

3.1.3 小麥 小麥化肥施用強度在0.2分位點、0.5分位點分別通過了10%、5%水平的顯著性檢驗，回歸系數亦由0.214下降到0.204；在0.8分位點上，小麥化肥施用強度對單產增長的作用不顯著。這表明小麥的化肥施用不僅存在邊際效用遞減特征，目前其增產效應已經不再顯著。

圖3 不同分位數下糧食單產的化肥施用彈性變化Fig. 3 Changes of chemical fertilizer application elasticities at different quantiles of grain yields

3.1.4 玉米 玉米化肥施用強度在上述3個分位點分別通過了1%、5%、1%水平的顯著性檢驗，回歸系數由0.2分位點的0.180下降到0.5分位點的0.152，在0.8分位點上仍為0.161，明顯高于同分位點上的稻谷和小麥，表明化肥施用強度的增加對玉米單產的提高仍有一定促進作用。

3.1.5 中秈稻 中秈稻化肥施用強度在0.2分位點通過了5%水平的顯著性檢驗，回歸系數為0.200，即，中秈稻化肥施用強度每提高10%，其單產約提高2%；但這在0.5分位點和0.8分位點均未通過顯著性檢驗，表明化肥的增產效應已不再顯著；盡管不顯著，其彈性系數也表現出下降趨勢。

3.1.6 早秈稻 早秈稻化肥施用強度在上述3個分位點均未通過顯著性檢驗，需進一步深入討論。

3.2 其他因素對糧食單產的影響

3.2.1 單位面積機械作業費 單位面積機械作業費對糧食單產的影響整體不顯著，只有少數品種在特定分位點處通過了顯著性檢驗。一種解釋是，農業機械作業偏重的是勞動替代性技術進步，可以提高勞動生產率，促進糧食規模經營，但對于提高土地生產率作用不明顯；另一種可能則是，部分農業機械使用可以提高糧食單產水平，但目前農業機械結構還不夠合理，譬如土壤深耕、秸稈還田、病蟲害防控防治等有助于提高糧食單產的農機配置不足，推廣作業率不高。模型結果顯示，如果單位面積機械作業費用通過了顯著性檢驗，其彈性系數均為正值，這或意味著農業機械作業在一定情形下是可以提高糧食單產的，表明后一種解釋可能更加合理。

3.2.2 農戶銷售價格 農戶銷售價格對糧食單產水平存在正向影響，其中小麥和玉米在上述3個分位點均通過了顯著性檢驗，早秈稻、中秈稻和粳稻則在2個分位點上通過了顯著性檢驗。這意味著，伴隨各種糧食農戶銷售價格的上漲，能夠有效調動農民種糧積極性，激勵其提高糧食單產水平。從回歸系數看，農戶銷售價格對小麥單產提高的影響較明顯，接下來依次是中秈稻、玉米、粳稻、早秈稻，對晚秈稻的影響較弱。

3.2.3 單位面積用工數量 單位面積用工數量對糧食單產的影響方向為負向，其中晚秈稻、粳稻、小麥、玉米4個品種通過了顯著性檢驗。可以理解為，盡管各種糧食的單位用工數量有所下降，但并沒有造成糧食單產的下降。換言之，農業機械化作業推廣和社會化服務發展，可以大幅度減少糧食生產用工數量、降低勞動強度，使得老齡勞動力和兼業勞動力種糧成為可能，并且也不會引起糧食單產和總產的滑坡，這一點在河南、河北等平原地區表現尤為突出；輔以新型農業經營主體的不斷涌現，商品糧生產專業化趨勢日益明顯，大大降低了種糧勞動力老齡化和兼業化對糧食生產和產量的影響。

3.2.4 自然因素 自然災害成災率對糧食單產的影響為負向，在大多數情形下都通過了顯著性檢驗，這符合一般性的經驗認知。此外，年降水量和年平均氣溫對糧食單產水平的影響整體上并不顯著，一方面年降水量和年平均氣溫對糧食單產的影響具有時段性特點，與作物的生物特性緊密有關，如水稻揚花期降水較多會影響產量，其他時期則影響不大；另一方面農田水利基礎設施的改善，也會降低氣候因素對單產水平的影響。

4 結論與對策建議

4.1 結論

本文探討了中國不同糧食品種、不同地區的化肥施用強度差異，進一步運用分位數回歸模型探究了1994-2013年化肥施用強度等因素對早秈稻、中秈稻、晚秈稻、粳稻、小麥、玉米6種糧食單產的影響。研究顯示：1）中國糧食作物的化肥施用強度持續增加，其中小麥化肥施用強度最高，玉米、粳稻次之，秈稻略低，但都已經大幅超過國際公認的安全施用標準。6種糧食作物的化肥施用強度逐步顯示出主產區與非主產區、優勢產區與非優勢產區的差異，中秈稻、粳稻、小麥和玉米同時表現出東部、中部和西部的差異。2）無論是6種糧食整體，還是具體品種，化肥施用強度的增加與其糧食單產水平的提高具有正向相關關系，但目前增加化肥施用所產生的邊際產量已經較小。3）化肥施用強度對晚秈稻、粳稻、玉米單產水平有正向影響，且隨著分位點由低分位向高分位移動，這種正向影響逐步趨于弱化，表明化肥施用已出現明顯的邊際效用遞減特征，但整體上仍有助于提高單產；小麥、早秈稻、中秈稻的正向影響已不再顯著，可能有過度施用的傾向；這可能說明不同糧食品種間化肥施用強度對其單產的邊際效應開始分化。4）若將糧食視為一個整體，農戶銷售價格對糧食單產有正向影響，單位面積用工數量、自然災害成災率為負向影響，單位面積機械作業費、年降水量和年平均氣溫影響不顯著，盡管不同品種略有差異。

4.2 政策建議

過去20多年，增加化肥施用對中國糧食持續增產、保障糧食安全發揮了重要支撐作用。農戶也形成了依仗化肥投入促進糧食增產的經驗傳統和路徑依賴，并導致了現階段化肥施用強度超出國際安全施用標準，不僅增加了糧食生產成本、推高了糧食價格，而且引發了耕地有機質減少、土壤營養失衡，農業面源污染加重、農業發展不可持續等系列問題。結合現階段的實際，提出以下對策建議。

第一，加強糧食綠色增產科技的推廣應用，要破除依仗化肥投入促進糧食增產的“經驗”依賴，轉變為依仗農業科技促進糧食增產的“科技”依賴，逐步實現糧食生產向資源節約、環境友好集約式增長轉變。為此，應堅持科學施肥、經濟施肥、環保施肥的理念，充分依靠科技進步和新型農業經營主體，推進精準施肥、減少盲目施肥，調整肥料使用結構、推廣高效新型肥料，改進施肥方式、推廣測土配方施肥，用有機肥替代化肥等科學施肥技術路徑，實現“減量、節本、提效”，到2020年實現化肥使用量“零增長”[17]。

第二，加快農業支持政策的調整轉型，由補貼農業投入品生產或使用，轉向補貼支持農業綠色生產和環境友好型生產。對農業投入品生產或使用進行補貼（如化肥運輸補貼、化肥淡季儲備補貼、農資綜合補貼），相當于降低投入品價格，進而鼓勵增加投入品使用，化肥亦是如此，此類政策應予以取消或調整支持方向。積極整合現有測土配方施肥補貼、秸稈還田補貼、有機肥補貼等系列政策，并力爭與農業經營主體的綠色生產、環境友好型生產等實際行為掛鉤（如減少化肥、農藥等使用），增強政策的針對性和有效性。對新型農業經營主體的支持和補貼，也應逐步與資源節約型和環境友好型生產掛鉤。

第三，加快革新糧食安全理念。中國糧食安全不僅要保障糧食供給的數量安全，也要保障糧食供給的質量安全；不僅要保障當代人的即期糧食安全，還要保障子孫后代的中長期糧食安全。因此，減少對農業資源的掠奪性索取，降低石油農業對資源和生態的影響，促進農業資源、生態、環境的協調發展，實現農業綠色發展、可持續發展，即成為糧食安全新理念所要堅守的底線。

[1] 史常亮, 郭焱, 朱俊峰. 中國糧食生產中化肥過量施用評價及影響因素研究[J]. 農業現代化研究, 2016, 37(4): 671-679.

Shi C L, Guo Y, Zhu J F. Evaluation of over fertilization in China and its influencing factors[J]. Research of Agricultural Modernization, 2016, 37(4): 671-679.

[2] 張北贏, 陳天林, 王兵. 長期施用化肥對土壤質量的影響[J].中國農學通報, 2010, 26(11): 182-187.

Zhang B Y, Chen T L, Wang B. Effects of long-term uses of chemical fertilizers on soil quality[J]. Chinese Agricultural Science Bulletin, 2010, 26(11): 182-187.

[3] 蔡榮. 農業化學品投入狀況及其對環境的影響[J]. 中國人口?資源與環境, 2010, 20(3): 107-110.

Cai R. Agri-chemicals inputs and its impact on environment[J]. China Population, Resources and Environment, 2010, 20(3): 107-110.

[4] 尚杰, 尹曉宇. 中國化肥面源污染現狀及其減量化研究[J]. 生態經濟, 2016, 32(5): 196-199.

Shang J, Yin X Y. Study on the present situation and reduction of fertilizer non-point source pollution of China[J]. Ecological Economy, 2016, 32(5): 196-199.

[5] 葛繼紅, 周曙東. 要素市場扭曲是否激發了農業面源污染——以化肥為例[J]. 農業經濟問題, 2012, 33(3): 92-96, 112.

Ge J H, Zhou S D. Does factor market distortions stimulate the agricultural non-point source pollution? A case study of fertilizer[J]. Issues in Agricultural Economy, 2012, 33(3): 92-96, 112.

[6] 張利庠, 彭輝, 靳興初. 不同階段化肥施用量對我國糧食產量的影響分析——基于1952-2006年30個省份的面板數據[J].農業技術經濟, 2008(4): 85-94.

Zhang L X, Peng H, Jin X C. The analysis about the effect of fertilizer on China’s grain production at different stages : Based on the panel data of 30 provinces in 1952-2006[J]. Journal of Agrotechnical Economics, 2008(4): 85-94.

[7] 農業部. 農業部關于印發《到2020年化肥使用量零增長行動方案》和《到2020年農藥使用量零增長行動方案》的通知[R].中華人民共和國農業部公報, 2015(3): 19-27.

The Ministry of Agriculture. The notification issued by the Ministry of Agriculture on the plans for zero growth in fertilizer and pesticides utilization[R]. The Communique of Ministry of Agriculture of the People’s Republic of China, 2015(3): 19-27.

[8] 張智峰, 張衛峰. 我國化肥施用現狀及趨勢[J]. 磷肥與復肥, 2008, 23(6): 9-12.

Zhang Z F, Zhang W F. The situation and trend of fertilizer application in China[J]. Phosphate and Compound Fertilizer, 2008, 23(6): 9-12.

[9] 劉忠, 隋曉晨. 中國區域化肥利用特征分析[J]. 資源科學, 2008, 30(6): 822-828.

Liu Z, Sui X C. Regional characteristics of fertilizer use in China[J]. Resources Science, 2008, 30(6): 822-828.

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

Ji Y Q, Zhang H, Lu W Y, et al. Differentiation, incomplete information and farmers excessive application of fertilizer[J]. Journal of Agrotechnical Economics, 2016(2): 14-22.

[11] 王祖力, 肖海峰. 化肥施用對糧食產量增長的作用分析[J]. 農業經濟問題, 2008, 29(8): 65-68.

Wang Z L, Xiao H F. The analysis about the effect of fertilizer on grain production growth[J]. Issues in Agricultural Economy, 2008, 29(8): 65-68.

[12] 曾希柏, 李菊梅. 中國不同地區化肥施用及其對糧食生產的影響[J]. 中國農業科學, 2004, 37(3): 387-392.

Zeng X B, Li J M. Fertilizer application and its effect on grain production in different counties of China[J]. Scientia Agricultura Sinica, 2004, 37(3): 387-392.

[13] 楊增旭, 韓洪云. 化肥施用技術效率及影響因素——基于小麥和玉米的實證分析[J]. 中國農業大學學報, 2011, 16(1): 140-147.

Yang Z X, Han H Y. Technical efficiency of fertilizer and its infuencing factors: Based on wheat and corn empirical study[J]. Journal of China Agricultural University, 2011, 16(1): 140-147.

[14] 郝令昕, Daniel Q N. 分位數回歸模型[M]. 上海: 格致出版社/上海人民出版社, 2012.

Hao L X, Daniel Q N. Quantile Regression Model[M]. Shanghai: Truth and Wisdom Press/The Shanghai People’s Publishing House, 2012.

[15] 朱滿德, 李辛一, 程國強. 綜合性收入補貼對中國玉米全要素生產率的影響分析——基于省級面板數據的DEA-Tobit兩階段法[J]. 中國農村經濟, 2015(11): 4-14, 53.

Zhu M D, Li X Y, Cheng G Q. An analysis of the impacts of comprehensive income subsidies to total factor productivity in maize production of China: A DEA-Tobit two-stage approach with province-level panel data during 2004-2008[J]. Chinese Rural Economy, 2015(11): 4-14, 53.

[16] 房麗萍, 孟軍. 化肥施用對中國糧食產量的貢獻率分析——基于主成分回歸C-D生產函數模型的實證研究[J]. 中國農學通報, 2013, 29(17): 156-160.

Fang L P, Meng J. Application of chemical fertilizer on grain yield in China analysis of contribution rate: Based on principal component regression C-D production function model and its empirical study[J]. Chinese Agricultural Science Bulletin, 2013, 29(17): 156-160.

[17] 李新華, 鞏前文. 從“增量增產”到“減量增效”: 農戶施肥調控政策演變及走向[J]. 農業現代化研究, 2016, 37(5): 877-884.

Li X H, Gong Q W. Trend and direction of China’s development of regulatory policies preventing over-fertilization in farming: From “increasing yield by increasing fertilizer quantity” to“increasing effcacy by reducing fertilizer quantity”[J]. Research of Agricultural Modernization, 2016, 37(5): 877-884.

（責任編輯：王育花）

The effects of fertilizer application intensity on grain yields in China: A quantile regression model analysis based on provincial panel data

ZHU Man-de1, LI Xin-yi1, XU Xue-gao2
（1. Guizhou University, a. School of Economics, b. China Center for Western Development Capacity Research, Guiyang, Guizhou 550025, China; 2. Institute of Agricultural Economics and Development, Jiangsu Province Academy of Agricultural Sciences, Nanjing, Jiangsu 210014, China）

F326.1 文獻標識碼：A 文章編號：1000-0275（2017）04-0649-09

10.13872/j.1000-0275.2017.0074

朱滿德, 李辛一, 徐雪高. 化肥施用強度對中國糧食單產的影響分析——基于省級面板數據的分位數回歸[J]. 農業現代化研究, 2017, 38(4): 649-657.

Zhu M D, Li X Y, Xu X G. The effects of fertilizer application intensity on grain yields in China: A quantile regression model analysis based on provincial panel data[J]. Research of Agricultural Modernization, 2017, 38(4): 649-657.

國家自然科學基金項目（71473052，71273069）；國務院發展研究中心2016年委托項目；貴州省教育廳青年學術創新人才項目（JYT2016004）。

朱滿德（1983-），男，安徽廬江人，博士，教授，主要從事農業經濟與政策、糧食經濟研究，E-mail：mdzhu@gzu.edu.cn；通訊作者：徐雪高（1981-），男，江蘇宜興人，博士，研究員，主要從事農業經濟與政策研究，E-mail：xuxuegaoxxg@163.com。

2017-04-06，接受日期：2017-06-28

Foundation item: National Natural Science Foundation of China (71473052, 71273069); Program of Development Research Center of the State Council; Program of the Academic Innovation Youth Talent of Department of Education at Guizhou Province (JYT2016004).

Corresponding author: XU Xue-gao, E-mail: xuxuegaoxxg@163.com.

Received 6 April, 2017;Accepted 28 June, 2017