周家俊　周德

[摘要]基于全國232個地級市數(shù)據(jù)，運用有限混合模型（FMM）分析了農(nóng)業(yè)機械投入對化肥減量化的作用效果。結(jié)果表明：總樣本由農(nóng)業(yè)機械化低水平樣本和高水平樣本組成;機械強度與化肥強度之間存在非線性關(guān)系;在機械化低水平階段，機械強度與化肥強度呈正相關(guān)關(guān)系;在機械化高水平階段，機械強度與化肥強度呈負相關(guān)關(guān)系;要通過增加農(nóng)業(yè)機械投入促進化肥減量化，關(guān)鍵在于轉(zhuǎn)變農(nóng)戶施肥需求。建議加快研發(fā)推廣大型農(nóng)業(yè)機械，推廣秸稈還田、化肥機械深施等技術(shù)，改善土壤理化性質(zhì)，提高化肥利用率，實現(xiàn)化肥減量化。

[關(guān)鍵詞]化肥減量化;農(nóng)業(yè)機械化;施肥方式;作用效果;有限混合模型

[中圖分類號]F323.3[文獻標(biāo)識碼]A

近年來，化肥過度施用造成了資源浪費、環(huán)境污染等一系列問題，制約了農(nóng)業(yè)可持續(xù)發(fā)展。研究農(nóng)業(yè)化肥減量化的有效途徑具有重要現(xiàn)實意義。為了加快實現(xiàn)化肥減量化目標(biāo)，農(nóng)業(yè)農(nóng)村部于2015年2月發(fā)布了《到2020年化肥使用量零增長行動方案》，為化肥減量化工作提供了政策條件。同時，秸稈還田、化肥機械深施、種肥同播等技術(shù)的大型農(nóng)業(yè)機械相繼出現(xiàn)，為農(nóng)業(yè)化肥施用方式的轉(zhuǎn)變提供了技術(shù)條件。由此，增加農(nóng)業(yè)機械投入是否能促進傳統(tǒng)化肥施用方式轉(zhuǎn)向精細化施肥進而實現(xiàn)化肥減量化？針對該問題的思考，有助于為實現(xiàn)化肥減量化提供有益參考。

1 文獻綜述

目前，對于我國化肥減量化實現(xiàn)路徑的研究仍較少，已有研究大多數(shù)為定性研究，從化肥施用對象、施用總量等角度分析總結(jié)我國化肥施用現(xiàn)狀，根據(jù)當(dāng)前不同作物的化肥施用強度提出減量化的作物對象優(yōu)先順序，也有學(xué)者關(guān)注中國化肥減量化的目標(biāo)值，基于農(nóng)產(chǎn)品供給和水資源壓力的雙重約束設(shè)定中國化肥減量化總體目標(biāo)。對于農(nóng)戶化肥施用行為影響因素的研究，主要集中在農(nóng)戶的認知風(fēng)險、個人特征、家庭特征層面，大多直接采用傳統(tǒng)線性回歸分析進行影響因素識別，或從結(jié)構(gòu)驅(qū)動效應(yīng)和效率驅(qū)動效應(yīng)的角度對化肥施用強度的變化進行分解;有學(xué)者研究了水稻種植這一單一作物生產(chǎn)類型中的農(nóng)業(yè)機械化技術(shù)對稻農(nóng)化肥減量化施用行為的影響，認為生產(chǎn)方式的轉(zhuǎn)變是實現(xiàn)化肥減量化的關(guān)鍵途徑。

關(guān)于機械化發(fā)展帶來的系列影響，從已有成果來看，主要集中在機械化與農(nóng)業(yè)總產(chǎn)出、糧食產(chǎn)出、生產(chǎn)效率、糧食播種面積、農(nóng)業(yè)生產(chǎn)中的勞動力要素替代等方面，認為農(nóng)業(yè)機械化發(fā)展能夠增加農(nóng)業(yè)產(chǎn)出，對糧食生產(chǎn)存在正的空間溢出效應(yīng)，對農(nóng)業(yè)經(jīng)濟發(fā)展方式轉(zhuǎn)變具有正向影響，但對糧食播種面積沒有顯著影響。少量研究從定性角度關(guān)注了機械化與生態(tài)環(huán)境的關(guān)系，認為農(nóng)業(yè)機械化發(fā)展與農(nóng)業(yè)資源利用效率提高、農(nóng)村環(huán)境改善和農(nóng)業(yè)可持續(xù)發(fā)展密切相關(guān)，應(yīng)大力研發(fā)節(jié)能、低碳農(nóng)業(yè)機械，發(fā)展精準(zhǔn)耕整、播種、施肥等技術(shù)，但仍未明確涉及對化肥減量化影響的作用效果。

總體而言，關(guān)于化肥減量化實現(xiàn)路徑的研究仍處于起步階段，存在較大空白，鮮有針對農(nóng)業(yè)機械投入對于化肥減量化效果的討論，因此，及時展開這方面的研究，除了為農(nóng)業(yè)生產(chǎn)安排和政策制定提供參考，對于補充新的研究案例也具有一定意義。

為分析農(nóng)業(yè)機械投入能否促進化肥減量化，本文假設(shè)一個地區(qū)的農(nóng)業(yè)機械化水平發(fā)展具有階段性，并且不同的階段存在一定概率混合，在不同階段農(nóng)戶的化肥施用方式和施用需求不同，農(nóng)業(yè)機械投入與化肥施用存在多階段演變關(guān)系。在此基礎(chǔ)上，以2015年232個地級市數(shù)據(jù)為樣本，運用有限混合模型（FMM）識別我國農(nóng)業(yè)機械化發(fā)展水平的異質(zhì)性構(gòu)成，分析農(nóng)業(yè)機械投入與化肥減量化兩者間關(guān)系的階段性變化，為加快實現(xiàn)我國化肥減量化目標(biāo)提供經(jīng)驗證據(jù)和政策建議。

2 模型與方法

化肥施用屬于農(nóng)業(yè)生產(chǎn)中的投入行為，經(jīng)濟學(xué)上可用追求利潤最大化的農(nóng)戶生產(chǎn)模型來描述農(nóng)戶選擇行為：在保證同等產(chǎn)出的情況下投入最少生產(chǎn)要素。然而農(nóng)戶的生產(chǎn)投入受到給定的技術(shù)條件限制，不同的機械水平帶來耕作模式變化，也導(dǎo)致化肥投入方式發(fā)生變化。由于農(nóng)業(yè)機械化發(fā)展水平的階段性，本文認為農(nóng)業(yè)機械投入與化肥施用兩者之間存在多階段演變關(guān)系。

農(nóng)戶層面，農(nóng)戶施肥方式在不同農(nóng)業(yè)機械化水平階段有所不同。在農(nóng)業(yè)機械化低水平階段，農(nóng)戶使用小型簡易機械進行表土疏松、作物收割，化肥施用方式以撒施、表施為主，屬于傳統(tǒng)施肥方式，養(yǎng)分損失較大，化肥利用率較低。在此階段，隨著機械投入的增加，農(nóng)戶為獲得最大利潤，提高化肥強度，直至化肥投入的邊際收益為零。隨著農(nóng)業(yè)機械化發(fā)展，進入農(nóng)業(yè)機械大型化高水平階段，大型農(nóng)業(yè)機械出現(xiàn)，農(nóng)戶可采用秸稈還田、土壤深耕、種肥同播、根施肥料等技術(shù)，能夠改良土壤營養(yǎng)特性，減少肥料流失，提高化肥利用率。因此，在農(nóng)機大型化高水平階段，隨著機械投入的增加，施肥方式轉(zhuǎn)變，為獲得利潤最大化，農(nóng)戶所需要投入的化肥更少，化肥強度隨機械強度提高而降低。

然而農(nóng)業(yè)生產(chǎn)與機械化發(fā)展水平往往情況多樣，較為復(fù)雜。不同農(nóng)戶的農(nóng)業(yè)機械投入不同，其施肥方式不同，其發(fā)展轉(zhuǎn)型的階段并不能清晰劃分。因此，可將一個地市的化肥施用情況理解為不同機械化水平施肥行為的混合：農(nóng)業(yè)機械化低水平的傳統(tǒng)施肥行為和農(nóng)業(yè)機械化高水平的施肥行為按照一定概率比例的混合，且每一類行為都有各自的分布。如果一個地市機械化低水平傳統(tǒng)施肥行為的概率比例大于機械化高水平施肥行為的概率比例，那么該市處于農(nóng)業(yè)機械化低水平階段;反之，處于高水平階段。同理，一個地市層面的樣本，是由處于不同機械化水平階段的地市按照一定概率比例的混合，各階段地市服從各自分布，地市之間具有異質(zhì)性。

上述分析可用包含兩個子分布的有限混合模型（FMM）表示，總樣本的概率密度函數(shù)可以寫成：

其中，fi（y|x'βi）是第i個子分布的概率密度函數(shù);x'是解釋變量向量;βi是第i個子分布的待估計解釋變量參數(shù);πi是第i個子分布的概率比例，0≤πi≤1，且。假定各子分布服從正態(tài)分布，πi服從多項logistic分布。

3 實證分析

本文收集了2015年232個地市數(shù)據(jù)，通過有限混合模型，定量考察在不同農(nóng)業(yè)機械化水平階段，農(nóng)業(yè)機械投入與化肥減量化之間的關(guān)系，以期揭示農(nóng)業(yè)機械投入水平與化肥減量化之間的關(guān)系。

3.1 指標(biāo)選取與數(shù)據(jù)來源

本文被解釋變量選取化肥強度，即單位農(nóng)作物播種面積對應(yīng)的化肥折純使用量。關(guān)鍵解釋變量選取機械強度，即單位農(nóng)作物播種面積對應(yīng)的農(nóng)業(yè)機械總動力投入。根據(jù)農(nóng)戶生產(chǎn)模型，化肥施用和機械投入還受到種植經(jīng)濟結(jié)構(gòu)、城鎮(zhèn)化程度和人力資本的影響，因此本文將以上變量納入控制變量。對于種植經(jīng)濟結(jié)構(gòu)，經(jīng)濟作物單位耕種面積帶來的經(jīng)濟效益普遍高于糧食作物，選取單位耕種面積對應(yīng)的種植業(yè)產(chǎn)值作為種植經(jīng)濟結(jié)構(gòu)的代理變量。對于城鎮(zhèn)化程度，選取常住城鎮(zhèn)人口占常住總?cè)丝诒壤鳛槌擎?zhèn)化率。對于人力資本，采用每萬人在校學(xué)生數(shù)代理人力資本。

根據(jù)2015年2月開始實施《到2020年化肥使用量零增長行動方案》，本文以此作為化肥減量化工作標(biāo)志性時間節(jié)點。本文通過各省統(tǒng)計年鑒收集了2015年全國232個地級市數(shù)據(jù)，涵蓋糧食主產(chǎn)省份中的內(nèi)蒙古、吉林、安徽、山東、江蘇、河北、河南、湖南、黑龍江和非主產(chǎn)省份中的廣東、浙江、海南、廣西、新疆、甘肅、貴州、陜西。

3.2 描述性分析

表1顯示了232個地級市5個指標(biāo)的描述性統(tǒng)計結(jié)果，圖1顯示了總樣本中機械強度的概率密度分布。對于機械強度，總樣本分布并不完全服從正態(tài)分布，考慮存在至少兩個分布的疊加。

3.3 回歸結(jié)果

為進行對照，先對包含232個地市的總樣本進行傳統(tǒng)OLS線性回歸，系數(shù)估計結(jié)果見表2，殘差分布見圖2。傳統(tǒng)OLS線性回歸顯示，在其他條件保持不變的情況下，隨著農(nóng)業(yè)機械投入增加，化肥強度下降，但不具有統(tǒng)計顯著性，且殘差分布不符合正態(tài)分布，總樣本分布存在起碼兩個分布的疊加。為此，有必要采用FMM作進一步分析。

總樣本可能由兩個及以上子分布混合組成，但當(dāng)子分布數(shù)目較多時，會增加有限混合模型估計的難度，且本文數(shù)據(jù)中，當(dāng)子分布數(shù)目大于兩個時，估計結(jié)果并不顯著。結(jié)合機械化水平存在兩階段特征，本文選擇子分布數(shù)目為兩個。

FMM回歸結(jié)果顯示（見表2），總樣本分布由子分布一占0.779，子分布二占0.221的概率比例混合組成，且子分布一的機械強度均值（7.652）小于子分布二（9.618），故子分布一對應(yīng)農(nóng)業(yè)機械化低水平階段，子分布二對應(yīng)農(nóng)業(yè)機械化高水平階段。

在農(nóng)業(yè)機械化低水平階段，農(nóng)戶使用低水平農(nóng)業(yè)機械，增加低水平機械投入未能轉(zhuǎn)變農(nóng)戶的傳統(tǒng)施肥方式，化肥強度與機械強度呈正相關(guān)關(guān)系。從邊際效應(yīng)看，在保持其他因素不變的情況下，農(nóng)業(yè)機械強度每提高1個單位（kW/hm2），化肥強度提高10.070個單位（kg/hm2），并在1%的水平下具有統(tǒng)計顯著性。

在農(nóng)業(yè)機械化高水平階段，大型農(nóng)業(yè)機械出現(xiàn)，帶來秸稈還田、化肥機械深施技術(shù)的改進，有助于改良土壤肥力，轉(zhuǎn)變農(nóng)戶的施肥需求，化肥強度與機械強度呈負相關(guān)關(guān)系。從邊際效應(yīng)看，在保持其他因素不變的情況下，農(nóng)業(yè)機械強度每提高1個單位（kW/hm2），化肥強度降低9.306個單位（kg/hm2）。

實證表明，在農(nóng)業(yè)機械化發(fā)展過程中，機械強度與化肥強度之間存在非線性兩階段演變關(guān)系，增加農(nóng)業(yè)機械投入能否促進化肥減量化，關(guān)鍵在于能否轉(zhuǎn)變農(nóng)戶的施肥需求。在機械化低水平階段，機械投入以簡易型機械為主，可實現(xiàn)的技術(shù)層次較低，如簡單的表土翻耕、收割，并未轉(zhuǎn)變農(nóng)戶施肥方式，機械強度與化肥強度呈正相關(guān)關(guān)系。在機械化高水平階段，大型農(nóng)業(yè)機械出現(xiàn)，秸稈還田、化肥深施、精準(zhǔn)根施等技術(shù)得到應(yīng)用，化肥利用率提高，化肥施用需求下降，機械強度與化肥強度呈負相關(guān)關(guān)系。

3.4 穩(wěn)健性檢驗

為減少農(nóng)業(yè)種植經(jīng)濟結(jié)構(gòu)不同導(dǎo)致化肥強度存在差異而對實證結(jié)果造成影響，本文對糧食主產(chǎn)省份的123個地市進行FMM分析（見表3）。結(jié)果表明，在機械化低水平階段，機械強度與化肥強度之間呈正相關(guān)關(guān)系，從邊際效應(yīng)看，保持其他因素不變，農(nóng)業(yè)機械強度每增加1個單位（kW/hm2），化肥強度增加17.253個單位（kg/hm2）;在機械化高水平階段，機械強度與化肥強度呈負相關(guān)關(guān)系，從邊際效應(yīng)看，保持其他因素不變，農(nóng)業(yè)機械強度增加1個單位（kW/hm2），化肥強度減少9.364個單位（kg/hm2），與理論設(shè)想一致。

4 研究結(jié)論

本文假設(shè)機械化發(fā)展水平具有階段性，由于機械化水平發(fā)展會導(dǎo)致農(nóng)戶化肥施用方式發(fā)生轉(zhuǎn)變，然而由于實際情況的復(fù)雜性，一個地區(qū)的農(nóng)業(yè)生產(chǎn)是農(nóng)業(yè)機械化高水平階段生產(chǎn)行為和低水平階段生產(chǎn)行為的概率混合;在機械化低水平階段，農(nóng)戶沿用傳統(tǒng)施肥方式;在機械化高水平階段，農(nóng)戶逐漸轉(zhuǎn)向精細化施肥。增加農(nóng)業(yè)機械投入對于化肥減量化的作用是非線性的兩階段演變過程。通過有限混合模型對232個地市組成的總樣本進行分類，分別探討在不同農(nóng)業(yè)機械化水平下農(nóng)業(yè)機械投入對化肥減量化的作用效果，結(jié)果表明：總樣本由77.9%的農(nóng)業(yè)機械化低水平階段樣本和22.1%的機械化高水平階段樣本組成;在機械化低水平階段，機械強度與化肥強度呈正相關(guān)關(guān)系，保持其他因素不變，機械強度提高1個單位（kW/hm2），化肥強度提高10.070個單位（kg/hm2）;在機械化高水平階段，機械強度與化肥強度呈負相關(guān)關(guān)系，保持其他因素不變，農(nóng)業(yè)機械強度增加1個單位（kW/hm2），化肥強度減少9.306個單位（kg/hm2）。該實證結(jié)果符合理論設(shè)想。

本文的啟示在于，在不同農(nóng)業(yè)機械化水平階段，機械投入對于化肥減量化的作用效果可能截然相反，要通過增加農(nóng)業(yè)機械投入促進化肥減量化，關(guān)鍵在于依靠高層次機械耕作技術(shù)，改良土壤性質(zhì)，提高化肥利用率，減少化肥施用需求。為此，本文建議加快研發(fā)推廣大型農(nóng)業(yè)機械，提高農(nóng)業(yè)機械化水平，推廣秸稈還田、化肥機械深施等技術(shù)，改善土壤理化性質(zhì)，提高化肥利用率，實現(xiàn)化肥減量化。

[參考文獻]

[1] 金書秦，張惠，吳娜偉.2016年化肥、農(nóng)藥零增長行動實施結(jié)果評估[J].環(huán)境保護，2018，46（01）：45-49.

[2] 金書秦，周芳，沈貴銀.農(nóng)業(yè)發(fā)展與面源污染治理雙重目標(biāo)下的化 肥減量路徑探析[J].環(huán)境保護，2015，43（08）：50-53.

[3] 虞祎，楊泳冰，胡浩，等.中國化肥減量目標(biāo)研究—基于滿足農(nóng)產(chǎn)品供給與水資源的雙重約束[J].農(nóng)業(yè)技術(shù)經(jīng)濟，2017（02）：102-110.

[4] 徐衛(wèi)濤，張俊飚，李樹明，等.循環(huán)農(nóng)業(yè)中的農(nóng)戶減量化投入行為分析—基于晉、魯、鄂三省的化肥投入調(diào)查[J].資源科學(xué)，2010，32（12）：2407-2412.

[5] 鞏前文，穆向麗，田志宏.農(nóng)戶過量施肥風(fēng)險認知及規(guī)避能力的影響因素分析—基于江漢平原284個農(nóng)戶的問卷調(diào)查[J].中國農(nóng)村經(jīng)濟，2010（10）：66-76.

[6] 鄭微微，徐雪高.江蘇省化肥施用強度變化驅(qū)動因子分解及其影響因素分析[J].華中農(nóng)業(yè)大學(xué)學(xué)報（社會科學(xué)版），2017（04）：55-62.

[7] 楊萬江，李琪.稻農(nóng)化肥減量施用行為的影響因素[J].華南農(nóng)業(yè)大學(xué)學(xué)報（社會科學(xué)版），2017，16（03）：58-66.

[8] 王新利，趙琨.黑龍江省農(nóng)業(yè)機械化水平對農(nóng)業(yè)經(jīng)濟增長的影響研究[J].農(nóng)業(yè)技術(shù)經(jīng)濟，2014（06）：31-37.

[9] 伍駿騫，方師樂，李谷成，等.中國農(nóng)業(yè)機械化發(fā)展水平對糧食產(chǎn)量的空間溢出效應(yīng)分析—基于跨區(qū)作業(yè)的視角[J].中國農(nóng)村經(jīng)濟，2017（06）：44-57.

[10] 劉成，馮瀟.中國農(nóng)業(yè)機械投入效率的三階段DEA的分析[J].中國農(nóng)機化學(xué)報，2018，39（03）：79-86.

[11] 楊進，吳比，金松青，等.中國農(nóng)業(yè)機械化發(fā)展對糧食播種面積的影響[J].中國農(nóng)村經(jīng)濟，2018（03）：89-104.

[12] 林善浪，葉煒，張麗華.農(nóng)村勞動力轉(zhuǎn)移有利于農(nóng)業(yè)機械化發(fā)展嗎—基于改進的超越對數(shù)成本函數(shù)的分析[J].農(nóng)業(yè)技術(shù)經(jīng)濟，2017（07）：4-17.

[13] 王旭東.我國農(nóng)業(yè)機械技術(shù)進步與農(nóng)業(yè)可持續(xù)發(fā)展關(guān)聯(lián)性探究[J].南方農(nóng)機，2017，48（17）：27-30.

[14] 羅錫文，廖娟，胡煉，等.提高農(nóng)業(yè)機械化水平促進農(nóng)業(yè)可持續(xù)發(fā)展[J].農(nóng)業(yè)工程學(xué)報，2016，32（01）：1-11.

[15] De Z， Yu X. Calorie Elasticities with Income Dynamics： Evidence from the Literature. Applied Economic Perspectives and Policy， 2015， 37（4）： 575–601.

[16] Marron J， Wand M. Exact Mean Integrated Squared Error. The Annals of Statistics， 1992， 2（2）： 712-736.