種植規(guī)模與細(xì)碎化對(duì)小農(nóng)戶耕地質(zhì)量保護(hù)行為的影響
——以蔬菜種植中農(nóng)藥、化肥施用為例

李 昊，銀敏華，馬彥麟，康燕霞，賈 瓊，齊廣平，汪精海

（1.蘭州大學(xué)經(jīng)濟(jì)學(xué)院，甘肅 蘭州 730000；2.甘肅農(nóng)業(yè)大學(xué)水利水電工程學(xué)院，甘肅 蘭州 730070）

1 引言

家庭聯(lián)產(chǎn)承包責(zé)任制的推廣，使農(nóng)民生產(chǎn)積極性明顯提高，農(nóng)業(yè)生產(chǎn)力水平不斷增強(qiáng)，但限于當(dāng)時(shí)正處于計(jì)劃經(jīng)濟(jì)向市場(chǎng)經(jīng)濟(jì)過渡的起始階段，農(nóng)業(yè)產(chǎn)出的增加在很大程度上依賴于勞動(dòng)的投入。隨著改革開放的深入推進(jìn)以及中國(guó)工業(yè)體系的建立，農(nóng)藥、化肥等現(xiàn)代農(nóng)業(yè)生產(chǎn)要素的產(chǎn)量和投入量不斷增加[1]，一方面替代了傳統(tǒng)勞動(dòng)投入，使現(xiàn)代農(nóng)業(yè)生產(chǎn)要素成為農(nóng)業(yè)生產(chǎn)的重要投入，增加了農(nóng)業(yè)產(chǎn)出，另一方面也形成了中國(guó)耕地污染不斷加重的基本格局。雖然近年來從中央到地方相繼出臺(tái)了一系列的耕地污染治理政策，在一定程度上抑制了農(nóng)藥、化肥的大量投入，但耕地污染問題仍未從根本上緩解[2-3]，目前仍對(duì)食品安全和人體健康構(gòu)成較大威脅。因此，如何緩解耕地污染已成為亟待解決的核心問題[4-6]。

從當(dāng)前學(xué)術(shù)界對(duì)小農(nóng)戶、規(guī)模農(nóng)業(yè)經(jīng)營(yíng)戶與耕地污染關(guān)系的討論來看，大體可分為小農(nóng)低效論和規(guī)模效率論兩種觀點(diǎn)。小農(nóng)低效論觀點(diǎn)認(rèn)為，人多地少的基本國(guó)情以及傳統(tǒng)小農(nóng)生產(chǎn)經(jīng)營(yíng)方式，使農(nóng)業(yè)生產(chǎn)很難實(shí)現(xiàn)規(guī)模經(jīng)濟(jì)，導(dǎo)致現(xiàn)代農(nóng)業(yè)生產(chǎn)要素的不當(dāng)施用[7-8]。規(guī)模效率論觀點(diǎn)認(rèn)為，以家庭農(nóng)場(chǎng)為代表的現(xiàn)代農(nóng)業(yè)經(jīng)營(yíng)主體在提高農(nóng)業(yè)生產(chǎn)規(guī)模效益的同時(shí)，由于規(guī)范化、規(guī)?；慕?jīng)營(yíng)模式，使其也兼具耕地保護(hù)的環(huán)境效益[9-10]。在實(shí)踐中，規(guī)模效率論的觀點(diǎn)得到了部分經(jīng)驗(yàn)證據(jù)的支撐，在當(dāng)前中國(guó)農(nóng)藥用量遠(yuǎn)高于世界平均水平以及化肥用量遠(yuǎn)超國(guó)際公認(rèn)安全施用上限的情況下，農(nóng)場(chǎng)制的經(jīng)營(yíng)方式能顯著降低農(nóng)藥施用強(qiáng)度，農(nóng)場(chǎng)規(guī)模每增加1%，農(nóng)藥用量可降低0.2%[11]；對(duì)大規(guī)模種植戶而言，種植規(guī)模越大，施用有機(jī)肥的可能性越高[12]；家庭農(nóng)場(chǎng)特別是規(guī)范化和適度規(guī)模經(jīng)營(yíng)的家庭農(nóng)場(chǎng)往往具有較高的技術(shù)效率[13]，能明顯降低農(nóng)藥施用概率，緩解現(xiàn)代農(nóng)業(yè)生產(chǎn)要素的密集施用[14]。

規(guī)?；?jīng)營(yíng)在一定程度上提高了農(nóng)業(yè)生產(chǎn)的規(guī)范化水平，降低了生產(chǎn)成本，這也是當(dāng)前較多學(xué)者嘗試推動(dòng)農(nóng)業(yè)規(guī)?；?jīng)營(yíng)以求化解食品安全風(fēng)險(xiǎn)和保障耕地生產(chǎn)潛力的重要原因之一。但也不可否認(rèn)，包干到戶后，以公平為實(shí)質(zhì)的土地分配造成了耕地細(xì)碎化問題并延續(xù)至今，無疑增加了農(nóng)戶農(nóng)業(yè)經(jīng)營(yíng)成本。此外，雖然向規(guī)模化農(nóng)業(yè)經(jīng)營(yíng)轉(zhuǎn)變是中國(guó)農(nóng)業(yè)現(xiàn)代化發(fā)展的必經(jīng)之路，但在可預(yù)見的時(shí)間范圍內(nèi)，受城鎮(zhèn)就業(yè)容納能力以及農(nóng)村勞動(dòng)力供給與市場(chǎng)需求結(jié)構(gòu)性失衡的影響，小農(nóng)戶仍將占據(jù)中國(guó)農(nóng)業(yè)生產(chǎn)的主體地位。在小農(nóng)戶農(nóng)業(yè)生產(chǎn)仍有強(qiáng)大生命力，強(qiáng)制、快速土地流轉(zhuǎn)實(shí)現(xiàn)規(guī)模經(jīng)營(yíng)已被證明違反中國(guó)農(nóng)業(yè)農(nóng)村發(fā)展規(guī)律，且以自愿為基礎(chǔ)的農(nóng)戶土地流轉(zhuǎn)仍未有效解決耕地細(xì)碎化的前提下，兩個(gè)邏輯相關(guān)的現(xiàn)實(shí)問題便亟待回答：對(duì)于小農(nóng)戶①根據(jù)第三次全國(guó)農(nóng)業(yè)普查公報(bào)對(duì)種植業(yè)規(guī)模農(nóng)業(yè)經(jīng)營(yíng)戶標(biāo)準(zhǔn)的界定，一年一熟制地區(qū)露地種植農(nóng)作物的土地達(dá)到100畝及以上，一年二熟及以上地區(qū)露地種植農(nóng)作物的土地達(dá)到50畝及以上。遠(yuǎn)高于本文戶均6.6畝的種植規(guī)模（最大值25畝），本文小農(nóng)戶指?jìng)鹘y(tǒng)意義上以家庭為生產(chǎn)單位，農(nóng)產(chǎn)品自用并將多余部分出售的小規(guī)模種植戶。而言，種植規(guī)模的擴(kuò)大是否有利于促進(jìn)其耕地質(zhì)量保護(hù)行為？若上述問題成立，則在當(dāng)前耕地細(xì)碎化的現(xiàn)實(shí)背景下，種植規(guī)模擴(kuò)大對(duì)小農(nóng)戶耕地質(zhì)量保護(hù)行為的促進(jìn)作用是否會(huì)趨于弱化？對(duì)上述問題的回答無疑是更為契合中國(guó)農(nóng)村農(nóng)業(yè)生產(chǎn)實(shí)際的核心議題，也是轉(zhuǎn)型期中國(guó)農(nóng)業(yè)現(xiàn)代化建設(shè)的應(yīng)有之義。

遺憾的是，目前學(xué)界和管理界探討的農(nóng)業(yè)規(guī)模經(jīng)營(yíng)所形成的環(huán)境效益主要聚焦于家庭農(nóng)場(chǎng)等現(xiàn)代農(nóng)業(yè)經(jīng)營(yíng)主體，對(duì)傳統(tǒng)“人均一畝三分、戶均不足十畝”的小農(nóng)戶關(guān)注不足，且以耕地細(xì)碎化問題切入農(nóng)戶生產(chǎn)領(lǐng)域的研究主要集中于農(nóng)業(yè)生產(chǎn)效率[15-17]，農(nóng)戶耕地質(zhì)量保護(hù)行為的研究亟待探討。鑒于此，本文以小農(nóng)戶為研究對(duì)象，在構(gòu)建種植規(guī)模、耕地細(xì)碎化與小農(nóng)戶耕地質(zhì)量保護(hù)行為理論關(guān)系的基礎(chǔ)上，以小農(nóng)戶農(nóng)藥、化肥施用行為為例，基于山東、河南和甘肅3省微觀調(diào)查數(shù)據(jù)，采用貝葉斯似不相關(guān)回歸模型、基于遺傳算法的傾向得分匹配模型和斜率分析等對(duì)所構(gòu)建的理論分析框架進(jìn)行驗(yàn)證。本文的邊際貢獻(xiàn)在于：從研究對(duì)象來看，著眼于仍占據(jù)中國(guó)農(nóng)業(yè)生產(chǎn)主體地位的小農(nóng)戶耕地質(zhì)量保護(hù)行為，可與現(xiàn)有關(guān)于規(guī)模種植戶帶來環(huán)境效益的相關(guān)研究互補(bǔ)，構(gòu)成中國(guó)農(nóng)戶耕地質(zhì)量保護(hù)行為研究的閉環(huán)；就研究方法而言，考慮了在中國(guó)現(xiàn)實(shí)農(nóng)情下小農(nóng)戶農(nóng)業(yè)生產(chǎn)不同行為之間存在的關(guān)聯(lián)性，以及農(nóng)戶調(diào)查數(shù)據(jù)不能完全符合傳統(tǒng)頻率統(tǒng)計(jì)關(guān)于變量服從正態(tài)分布的假定，采用更為靈活的貝葉斯似不相關(guān)回歸模型進(jìn)行分析，使分析結(jié)果有更高精度。以期為中國(guó)耕地污染治理特別是傳統(tǒng)小農(nóng)生產(chǎn)經(jīng)營(yíng)模式下的耕地污染治理提供一個(gè)較為契合農(nóng)情和農(nóng)業(yè)生產(chǎn)實(shí)際的理論依據(jù)和經(jīng)驗(yàn)參考。

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

探討種植規(guī)模、耕地細(xì)碎化與小農(nóng)戶耕地質(zhì)量保護(hù)行為的關(guān)系是本文的核心。因此，在分別論述種植規(guī)模和耕地細(xì)碎化對(duì)小農(nóng)戶耕地質(zhì)量保護(hù)行為影響的基礎(chǔ)上，進(jìn)一步探索在耕地細(xì)碎化條件下，種植規(guī)模對(duì)小農(nóng)戶耕地質(zhì)量保護(hù)行為的影響，構(gòu)成本文的理論邏輯脈絡(luò)。

2.1 種植規(guī)模對(duì)小農(nóng)戶耕地質(zhì)量保護(hù)行為的影響

在傳統(tǒng)的小農(nóng)生產(chǎn)方式下，種植規(guī)模至少以三種方式影響小農(nóng)戶耕地質(zhì)量保護(hù)行為。第一，在新古典經(jīng)濟(jì)學(xué)體系下，小農(nóng)生產(chǎn)的邊際產(chǎn)出呈倒“U”型結(jié)構(gòu)，在邊際技術(shù)替代率遞減規(guī)律作用下，種植規(guī)模擴(kuò)大導(dǎo)致生產(chǎn)效率提高[18]，特別是在小農(nóng)戶種植規(guī)模遠(yuǎn)低于家庭農(nóng)場(chǎng)種植規(guī)模的情況下，形成以農(nóng)業(yè)邊際產(chǎn)出遞增為實(shí)質(zhì)的生產(chǎn)過程會(huì)降低現(xiàn)代農(nóng)業(yè)生產(chǎn)要素的投入，進(jìn)而促進(jìn)農(nóng)戶耕地質(zhì)量保護(hù)行為；第二，由提高家庭收入動(dòng)機(jī)所驅(qū)動(dòng)，在當(dāng)前小農(nóng)戶留鄉(xiāng)務(wù)農(nóng)和進(jìn)城務(wù)工存在較大收入差距的條件下，受比較收益的影響，種植規(guī)模較小的農(nóng)戶往往有動(dòng)機(jī)脫離農(nóng)業(yè)生產(chǎn)束縛[19]，選擇進(jìn)城務(wù)工以增加家庭收入，在農(nóng)業(yè)生產(chǎn)中更傾向于采用節(jié)省勞動(dòng)時(shí)間的化肥、農(nóng)藥投入來替代勞動(dòng)力投入[20]，甚至形成化肥對(duì)有機(jī)肥的反向替代，但隨著種植規(guī)模的擴(kuò)大與農(nóng)戶農(nóng)業(yè)收入的提高，非農(nóng)收入比較優(yōu)勢(shì)的式微會(huì)逐步緩解這一趨勢(shì)，促進(jìn)農(nóng)戶耕地質(zhì)量保護(hù)行為；第三，從比較靜態(tài)分析的方法論角度來看，農(nóng)戶種植規(guī)模的提高強(qiáng)化了家庭對(duì)農(nóng)業(yè)收入的依賴，使農(nóng)戶更注重耕地地力的培育以實(shí)現(xiàn)長(zhǎng)期增產(chǎn)的目標(biāo)，從而提高農(nóng)戶有機(jī)肥和綠色農(nóng)藥的施用[21-22]?；诖耍疚奶岢觯?/p>

H1：種植規(guī)模的提高能使小農(nóng)戶通過更合理地施用農(nóng)業(yè)生產(chǎn)要素而促進(jìn)其耕地質(zhì)量保護(hù)行為。

2.2 耕地細(xì)碎化對(duì)小農(nóng)戶耕地質(zhì)量保護(hù)行為的影響

改革開放后，以公平為導(dǎo)向的中國(guó)農(nóng)村土地分配制度客觀上形成了當(dāng)前耕地細(xì)碎化的現(xiàn)實(shí)，增加了農(nóng)業(yè)生產(chǎn)成本，但鑒于包干到戶之初小農(nóng)戶生產(chǎn)積極性的釋放，以及在城鄉(xiāng)二元結(jié)構(gòu)體制下，農(nóng)民進(jìn)城務(wù)工的通道尚未完全打開，耕地細(xì)碎化所增加的農(nóng)業(yè)生產(chǎn)成本被農(nóng)村充足的勞動(dòng)力隱性消解。但隨著20世紀(jì)90年代初中國(guó)逐步承接國(guó)際產(chǎn)業(yè)資本轉(zhuǎn)移，以及隨后糧票制度的取消和市場(chǎng)經(jīng)濟(jì)體制改革的深化，限制農(nóng)村勞動(dòng)力流動(dòng)的制度藩籬被沖破，越來越多的農(nóng)民離土離鄉(xiāng)、進(jìn)廠進(jìn)城，使原本農(nóng)村充裕的勞動(dòng)力資源具有了稀缺性，耕地細(xì)碎化所形成的生產(chǎn)成本不斷顯化[23-24]。一方面，隨著近年來全國(guó)勞動(dòng)力市場(chǎng)的放開和農(nóng)民工隊(duì)伍的形成，在農(nóng)村勞動(dòng)力農(nóng)業(yè)生產(chǎn)成本不斷顯化的背景下，耕地細(xì)碎化無疑增加了農(nóng)業(yè)勞動(dòng)成本[25]，導(dǎo)致越來越多的農(nóng)業(yè)勞動(dòng)力舍棄了傳統(tǒng)農(nóng)業(yè)的精耕細(xì)作，轉(zhuǎn)而采用化肥、農(nóng)藥等現(xiàn)代生產(chǎn)要素代替勞動(dòng)投入；另一方面，耕地細(xì)碎化條件下的農(nóng)業(yè)生產(chǎn)方式與現(xiàn)代農(nóng)業(yè)生產(chǎn)方式難以契合，導(dǎo)致更具規(guī)范化、統(tǒng)一化的農(nóng)業(yè)社會(huì)化服務(wù)體系難以介入[26]，在降低農(nóng)業(yè)生產(chǎn)效率的同時(shí)，也在一定程度上導(dǎo)致農(nóng)戶化肥、農(nóng)藥投入的隨意性，增加了化肥、農(nóng)藥不當(dāng)施用的風(fēng)險(xiǎn)。基于此，本文提出：

H2：耕地細(xì)碎化的現(xiàn)實(shí)背景制約了小農(nóng)戶耕地質(zhì)量保護(hù)行為。

2.3 種植規(guī)模、耕地細(xì)碎化與小農(nóng)戶耕地質(zhì)量保護(hù)行為

中國(guó)農(nóng)村土地制度改革的歷史路徑與農(nóng)村不同發(fā)展階段的現(xiàn)實(shí)情境，共同決定了小農(nóng)生產(chǎn)方式下種植規(guī)模的擴(kuò)大與耕地細(xì)碎化交織并存，形成當(dāng)前農(nóng)業(yè)生產(chǎn)中以規(guī)?；癁榇淼默F(xiàn)代性與以細(xì)碎化為代表的歷史性的對(duì)抗。在此背景下，種植規(guī)模對(duì)小農(nóng)戶耕地質(zhì)量保護(hù)行為的影響必然受耕地細(xì)碎化的制約而與預(yù)期存在偏差。當(dāng)耕地細(xì)碎化程度較低且種植規(guī)模較大的情況下，連片的生產(chǎn)經(jīng)營(yíng)模式使小農(nóng)戶農(nóng)業(yè)生產(chǎn)更接近于規(guī)?；纳a(chǎn)方式，能有效與現(xiàn)代農(nóng)業(yè)生產(chǎn)方式對(duì)接[27]，降低了農(nóng)業(yè)社會(huì)化服務(wù)體系準(zhǔn)入的自然條件制約，在提高小農(nóng)戶標(biāo)準(zhǔn)化、規(guī)范化生產(chǎn)的同時(shí)，降低了農(nóng)藥、化肥的不合理施用行為[28]；反之，即便有較大的種植規(guī)模，也往往由于耕地細(xì)碎化程度的提高而增加農(nóng)業(yè)勞動(dòng)成本，削減甚至瓦解農(nóng)業(yè)生產(chǎn)的規(guī)模效率，促使小農(nóng)戶通過現(xiàn)代農(nóng)業(yè)生產(chǎn)要素的投入以緩解勞動(dòng)的高強(qiáng)度投入，增加農(nóng)藥、化肥不當(dāng)施用的概率[29]。故本文提出：

H3：源于種植規(guī)模擴(kuò)大所形成的小農(nóng)戶耕地質(zhì)量保護(hù)行為的現(xiàn)代性優(yōu)勢(shì)會(huì)由于歷史因素形成的耕地細(xì)碎化而趨于弱化。

根據(jù)上述分析，本文構(gòu)建的理論分析框架見圖1。

圖1 理論分析框架Fig.1 Theoretical analysis framework

3 數(shù)據(jù)來源與研究方法

3.1 數(shù)據(jù)來源

本文數(shù)據(jù)來源于2021年對(duì)山東、河南和甘肅3省微觀農(nóng)戶調(diào)查，數(shù)據(jù)內(nèi)容由經(jīng)過培訓(xùn)的調(diào)查員負(fù)責(zé)填寫。鑒于經(jīng)濟(jì)作物種植戶特別是蔬菜種植戶現(xiàn)代農(nóng)業(yè)生產(chǎn)要素的不當(dāng)施用更為常見，本文調(diào)查對(duì)象聚焦于菜農(nóng)。省份的選擇綜合考慮了2020年《中國(guó)統(tǒng)計(jì)年鑒》中各省蔬菜播種面積（不低于2×105hm2）和區(qū)域經(jīng)濟(jì)發(fā)展水平，在中國(guó)東中西部3類地區(qū)分別隨機(jī)抽取1個(gè)省，而后按照類似的邏輯在每個(gè)省選取3個(gè)市（縣）。最后采用隨機(jī)抽樣的方式，在每個(gè)市（縣）選取2～3個(gè)鄉(xiāng)（鎮(zhèn)）、每個(gè)鄉(xiāng)（鎮(zhèn)）選取3～4個(gè)村，每村根據(jù)小農(nóng)戶數(shù)量隨機(jī)發(fā)放問卷15～20份，共發(fā)放問卷1 500份，收回有效問卷1 448份，有效率為96.5%。

3.2 變量選取與樣本描述性統(tǒng)計(jì)

因變量。現(xiàn)代農(nóng)業(yè)生產(chǎn)要素特別是農(nóng)藥、化肥的不當(dāng)施用是導(dǎo)致耕地污染的主要原因，也是近年來包括中央1號(hào)文件在內(nèi)的涉農(nóng)文件多次提出有機(jī)肥替代化肥、合理施用農(nóng)藥的原因。故本文小農(nóng)戶耕地質(zhì)量保護(hù)行為重點(diǎn)關(guān)注農(nóng)藥和肥料施用。參考已有研究的處理方式[30-31]，本文將農(nóng)藥施用行為設(shè)定為施用劇毒高毒化學(xué)農(nóng)藥、低毒化學(xué)農(nóng)藥、綠色農(nóng)藥、不施藥，分別賦值1～4；肥料施用行為設(shè)定為完全施用化肥、主要施用化肥、化肥有機(jī)肥混用、主要施用有機(jī)肥、完全施用有機(jī)肥，分別賦值1～5；將受訪者農(nóng)藥和肥料施用行為的賦值加總構(gòu)成小農(nóng)戶耕地質(zhì)量保護(hù)行為（表1）。

表1 變量定義Tab.1 Variable definition

核心自變量。種植規(guī)模和耕地細(xì)碎化是本文的核心自變量，參照相關(guān)研究的測(cè)量方式，分別用小農(nóng)戶蔬菜種植面積和蔬菜種植地塊數(shù)來衡量[32-33]。調(diào)查情況表明，受訪戶平均蔬菜種植面積6.6畝，較為符合當(dāng)前小農(nóng)戶“人均一畝三分、戶均不足十畝”的經(jīng)營(yíng)現(xiàn)狀；平均蔬菜種植地塊數(shù)接近4塊，農(nóng)戶耕種較為不便，對(duì)于本文抽樣地點(diǎn)為地勢(shì)相對(duì)較為平坦的非山區(qū)而言，耕地細(xì)碎化現(xiàn)象相對(duì)較為明顯。

控制變量。依據(jù)農(nóng)戶行為分析的基本范式[31]，本文選取受訪者性別、年齡、受教育程度等人口統(tǒng)計(jì)學(xué)特征[34]和家庭年收入、農(nóng)業(yè)收入占比、務(wù)農(nóng)時(shí)間等社會(huì)經(jīng)濟(jì)特征[35]作為控制變量。從調(diào)查情況來看，男性受訪者占比68.3%，平均年齡49.9歲，較為符合中國(guó)“老人農(nóng)業(yè)”的基本生產(chǎn)格局；受教育程度普遍偏低，近60.0%的受訪者為初中及初中以下學(xué)歷；家庭勞動(dòng)力人數(shù)普遍集中在3人以下，占比82.7%；平均家庭年收入4.4萬元；平均農(nóng)業(yè)收入占比超過50.0%。從人口統(tǒng)計(jì)學(xué)特征和社會(huì)經(jīng)濟(jì)特征來看，受訪對(duì)象較為符合當(dāng)前中國(guó)小農(nóng)戶生產(chǎn)的基本狀況。

3.3 研究方法

3.3.1 貝葉斯似不相關(guān)回歸模型

小農(nóng)戶農(nóng)藥、肥料施用行為均為有序分類變量，通常的研究方法是采用有序Logit模型或OLS回歸分別進(jìn)行分析。但小農(nóng)戶自身往往由于個(gè)體或家庭稟賦因素而呈現(xiàn)出行為間的相似性，導(dǎo)致農(nóng)藥、肥料施用行為具有內(nèi)在關(guān)聯(lián)，忽略這一因素將造成結(jié)果偏誤。此外，傳統(tǒng)的頻率統(tǒng)計(jì)依據(jù)大數(shù)定律和中心極限定理，認(rèn)為隨樣本量的增加，變量更傾向于服從正態(tài)分布，但現(xiàn)實(shí)中往往由于問卷測(cè)量方式的不同，使部分變量無法滿足這一要求。因此，本文根據(jù)李昊等[31]的建議，選擇更為靈活的貝葉斯似不相關(guān)回歸模型（Bayesian Seemingly Unrelated Regression，BSUR）對(duì)數(shù)據(jù)進(jìn)行分析，并將柯西分布作為參數(shù)估計(jì)的先驗(yàn)分布，公式如下：

式（1）—式（4）中：yni為第i個(gè)受訪者的第n種行為，即農(nóng)藥或肥料施用行為；xji為第i個(gè)受訪者的第j個(gè)控制變量；x規(guī)模和x細(xì)碎化分別表示種植規(guī)模和耕地細(xì)碎化程度；εn為殘差項(xiàng)；β為回歸系數(shù)；P（θ|yni，x）為后驗(yàn)概率密度函數(shù)；P（yni|θ，x）為似然函數(shù)；P（θ，x）為參數(shù)θ與自變量x的聯(lián)合分布；C（yni|μ，σ）為柯西分布的概率密度函數(shù)；μ和σ分別為均值和標(biāo)準(zhǔn)差。當(dāng)殘差協(xié)方差矩陣非對(duì)角元素不顯著時(shí)，退化為傳統(tǒng)單模型回歸結(jié)果。

3.3.2 基于遺傳算法的傾向得分匹配模型

在不同耕地細(xì)碎化程度下，小農(nóng)戶由于人口統(tǒng)計(jì)學(xué)特征和社會(huì)經(jīng)濟(jì)特征的差異可能出現(xiàn)對(duì)土地流轉(zhuǎn)的不同選擇，造成結(jié)果的內(nèi)生性。故本文采用基于遺傳算法的傾向得分匹配模型平衡耕地細(xì)碎化程度高和低（按均值劃分）兩種情況下受訪者人口統(tǒng)計(jì)學(xué)特征和社會(huì)經(jīng)濟(jì)特征，構(gòu)建反事實(shí)分析框架：

式（5）中：ATT表示實(shí)驗(yàn)組平均處理效應(yīng)；E（Y1|Ii= 1）和E（Y0|Ii= 1）分別表示同一受訪者在耕地細(xì)碎化程度高和低時(shí)的耕地質(zhì)量保護(hù)行為。對(duì)同一個(gè)體而言，只能觀測(cè)到其耕地細(xì)碎化程度高或低，二者不會(huì)同時(shí)出現(xiàn)。學(xué)術(shù)界普遍采用廣義線性回歸模型結(jié)果作為傾向得分[36]，而后基于核匹配或局部線性回歸等方式平衡實(shí)驗(yàn)組與對(duì)照組間的差異，雖然為解決可能存在的內(nèi)生性提供了條件，但也可能導(dǎo)致不存在顯著差異的變量在匹配后反而有顯著差異，引入新的估計(jì)誤差。因此，本文采用精度更高的基于遺傳算法的傾向得分匹配模型規(guī)避上述問題：

式（6）中：D為距離函數(shù)；與分別表示耕地細(xì)碎化程度高和低時(shí)受訪者的人口統(tǒng)計(jì)學(xué)特征或社會(huì)經(jīng)濟(jì)特征；W為權(quán)重；s為組間協(xié)方差矩陣；T為轉(zhuǎn)置矩陣。

4 結(jié)果與分析

4.1 種植規(guī)模、細(xì)碎化對(duì)小農(nóng)戶耕地質(zhì)量保護(hù)行為影響的初步證據(jù)

為驗(yàn)證種植規(guī)模、耕地細(xì)碎化程度對(duì)小農(nóng)戶耕地質(zhì)量保護(hù)行為的影響，以及驗(yàn)證種植規(guī)模對(duì)小農(nóng)戶耕地質(zhì)量保護(hù)行為的影響過程是否隨耕地細(xì)碎化程度的提高而弱化，本文采用貝葉斯似不相關(guān)回歸模型進(jìn)行分析。同時(shí)，為保持結(jié)果穩(wěn)健性，將小農(nóng)戶農(nóng)藥施用行為和肥料施用行為的賦值加總構(gòu)成整體行為，采用OLS進(jìn)行回歸。其中，貝葉斯似不相關(guān)回歸模型采用HMC（Hamiltonian Monte Carlo）抽樣算法抽樣20 000次，“熱身期”10 000次后馬爾科夫雙鏈迭代估計(jì)，如圖2所示（僅報(bào)告核心變量）；在引入交互項(xiàng)前，本文采用方差膨脹因子（VIF）檢驗(yàn)了多重共線性，結(jié)果表明，VIF均小于1.5，說明本文回歸結(jié)果不存在多重共線性。為避免引入交互項(xiàng)后出現(xiàn)多重共線性，本文對(duì)種植規(guī)模和細(xì)碎化進(jìn)行中心化處理。結(jié)果如表2所示。

圖2 馬爾科夫雙鏈迭代路徑圖Fig.2 Path diagram of Markov double-chain iteration

表2 貝葉斯似不相關(guān)回歸和OLS回歸結(jié)果Tab.2 Results of Bayesian seemingly unrelated regression and ordinary least squares regression

由圖2分析可知，馬爾科夫雙鏈迭代路徑重合，模型已收斂。由表2分析可知，小農(nóng)戶農(nóng)藥施用行為和肥料施用行為回歸結(jié)果的殘差協(xié)方差矩陣非對(duì)角元素顯著為正，說明對(duì)小農(nóng)戶而言，作為農(nóng)業(yè)生產(chǎn)的要素投入，農(nóng)藥和肥料的施用行為有相關(guān)性，采用貝葉斯似不相關(guān)回歸模型比傳統(tǒng)的單模型方法更有優(yōu)勢(shì)。

從結(jié)果來看，種植規(guī)模的擴(kuò)大使小農(nóng)戶更傾向于合理施用農(nóng)藥和肥料，促進(jìn)小農(nóng)戶耕地質(zhì)量保護(hù)行為，從而緩解耕地污染，假設(shè)H1初步得證。耕地細(xì)碎化程度的提高，增加了小農(nóng)戶農(nóng)業(yè)生產(chǎn)成本，明顯降低了小農(nóng)戶合理施用農(nóng)藥、肥料的可能性，阻礙了小農(nóng)戶耕地質(zhì)量保護(hù)行為的出現(xiàn)，假設(shè)H2初步得證。實(shí)質(zhì)上，耕地細(xì)碎化程度的提高，使現(xiàn)代農(nóng)業(yè)服務(wù)體系難以充分發(fā)揮助農(nóng)的作用，導(dǎo)致規(guī)范化、統(tǒng)一化的農(nóng)業(yè)生產(chǎn)要素投入的可能性下降，增加了小農(nóng)戶農(nóng)藥、化肥施用的隨意性，這也可能是小農(nóng)戶不合理施用農(nóng)藥、肥料的原因。種植規(guī)模與耕地細(xì)碎化程度的交互項(xiàng)顯著為負(fù)，說明雖然隨著種植規(guī)模的擴(kuò)大，小農(nóng)戶合理施用農(nóng)藥、肥料的耕地質(zhì)量保護(hù)行為更為明顯，但較高的耕地細(xì)碎化程度顯著降低了種植規(guī)模擴(kuò)大對(duì)小農(nóng)戶耕地質(zhì)量保護(hù)行為的正向影響，假設(shè)H3初步得證。

就控制變量而言，與女性相比，男性更傾向于采用綠色農(nóng)藥甚至不施藥等合理施藥行為，但卻并不影響其施肥行為，可能與男性更傾向以化肥替代勞動(dòng)投入，將節(jié)省的勞動(dòng)力投入到比較收益更高的非農(nóng)產(chǎn)業(yè)有關(guān)。隨著年齡的增大，小農(nóng)戶合理施藥與合理施肥的可能性增大。隨著勞動(dòng)力人數(shù)和家庭人數(shù)的增加，小農(nóng)戶合理施用農(nóng)藥、肥料進(jìn)而保護(hù)耕地質(zhì)量的可能性增大，為前文性別和年齡的解釋提供了佐證。隨著農(nóng)業(yè)收入占家庭總收入比例的增加，小農(nóng)戶對(duì)農(nóng)業(yè)生產(chǎn)的依賴性增大，促使其更可能通過合理施藥和施肥來保護(hù)耕地質(zhì)量。種植年限的增加促進(jìn)了小農(nóng)戶合理施藥行為，但對(duì)施肥行為無影響，可能是隨著種植年限的增加，小農(nóng)戶更了解不合理的農(nóng)藥施用對(duì)食品安全造成較大威脅，但肥料施用對(duì)食品安全的影響則沒有農(nóng)藥明顯所致。受教育程度和耕地污染認(rèn)知程度的提高均促進(jìn)了小農(nóng)戶耕地質(zhì)量保護(hù)行為。

4.2 基于遺傳算法的傾向得分匹配模型再分析

隨著中國(guó)農(nóng)村“三提五統(tǒng)”的取消，以村集體為單位的農(nóng)村耕地“五年一小調(diào)、十年一大調(diào)”這一傳統(tǒng)分地形式由于村集體資源配置功能弱化而無法持續(xù)，在一定程度上固化了當(dāng)前耕地細(xì)碎化的基本格局。隨著近年來農(nóng)業(yè)現(xiàn)代化的推進(jìn)，耕地細(xì)碎化與農(nóng)業(yè)現(xiàn)代化生產(chǎn)方式難以契合，催生了具有一定稟賦特征的小農(nóng)戶自發(fā)形成以耕地互換甚至“送人耕種”為實(shí)質(zhì)的土地流轉(zhuǎn)，即在耕地細(xì)碎化程度存在差異的情況下，由于小農(nóng)戶人口統(tǒng)計(jì)學(xué)特征和社會(huì)經(jīng)濟(jì)特征的不同，形成了種植規(guī)模的“自選擇”過程，從而影響本文種植規(guī)模這一核心變量對(duì)小農(nóng)戶耕地質(zhì)量保護(hù)行為的因果推斷。鑒于此，本文以耕地細(xì)碎化程度為分類變量，按其均值劃分為高和低兩組，采用基于遺傳算法的傾向得分匹配模型平衡小農(nóng)戶人口統(tǒng)計(jì)學(xué)特征和社會(huì)經(jīng)濟(jì)特征的差異，同時(shí)，為避免數(shù)據(jù)排序的不同影響抽樣結(jié)果，本文采用有放回抽樣的方式進(jìn)行匹配，共形成新樣本1 392個(gè)。傾向得分匹配前后的結(jié)果對(duì)比如表3所示，以匹配后新樣本為基礎(chǔ)的貝葉斯似不相關(guān)（BSUR）回歸結(jié)果和OLS回歸結(jié)果如表4所示。

表3 傾向得分匹配前后結(jié)果對(duì)比Tab.3 Comparison of results before and after propensity score matching

表4 種植規(guī)模、細(xì)碎化對(duì)小農(nóng)戶耕地質(zhì)量保護(hù)行為影響的再檢驗(yàn)（n=1 392）Tab.4 Re-examination of the effect of planting scale and fragmentation on the behavior of smallholders’ farmland quality protection（n=1 392）

由表3分析可知，與匹配前相比，匹配后在耕地細(xì)碎化程度高和低兩組間，小農(nóng)戶人口統(tǒng)計(jì)學(xué)特征和社會(huì)經(jīng)濟(jì)特征變量的標(biāo)準(zhǔn)均值差均有所下降，且匹配后組間標(biāo)準(zhǔn)均值差均不顯著，表明匹配效果較好。

由表4分析可知，在傾向得分匹配后采用新樣本進(jìn)行回歸，小農(nóng)戶農(nóng)藥施用行為和肥料施用行為的模型回歸殘差協(xié)方差矩陣非對(duì)角元素仍顯著為正，再次說明了采用貝葉斯似不相關(guān)（BSUR）回歸模型的合理性。從回歸結(jié)果來看，即便在考慮可能存在的內(nèi)生性后，種植規(guī)模對(duì)小農(nóng)戶農(nóng)藥施用行為、肥料施用行為和整體行為仍存在顯著的正向影響，再次驗(yàn)證了隨著種植規(guī)模的擴(kuò)大，小農(nóng)戶更傾向于保護(hù)耕地質(zhì)量的假設(shè)（H1）；隨著耕地細(xì)碎化程度的增加，地塊數(shù)量的增多使小農(nóng)戶農(nóng)業(yè)生產(chǎn)更為不便，降低了合理施藥與合理施肥的可能性，抑制了小農(nóng)戶耕地質(zhì)量保護(hù)行為，假設(shè)H2得證；此外，在種植規(guī)模對(duì)小農(nóng)戶耕地質(zhì)量保護(hù)行為的正向影響中，由于耕地細(xì)碎化的存在，顯著阻礙了種植規(guī)模擴(kuò)大對(duì)小農(nóng)戶耕地質(zhì)量保護(hù)行為的正向影響，假設(shè)H3得證。至此，本文所提出的三個(gè)假設(shè)全部得證。

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

為保證分析結(jié)果的穩(wěn)健性，本文將農(nóng)戶農(nóng)藥和化肥減量作為因變量，再次采用貝葉斯似不相關(guān)回歸模型進(jìn)行分析，結(jié)果如表5所示。

表5 穩(wěn)健性檢驗(yàn)結(jié)果Tab.5 Result of robustness test

由表5分析可知，種植規(guī)模的擴(kuò)大顯著促進(jìn)了小農(nóng)戶農(nóng)藥和化肥減量施用，但耕地細(xì)碎化程度的提高卻對(duì)小農(nóng)戶農(nóng)藥和化肥減量施用有顯著的負(fù)向影響，同時(shí)，在耕地細(xì)碎化的條件下，種植規(guī)模擴(kuò)大對(duì)小農(nóng)戶農(nóng)藥和化肥減量的影響明顯受到抑制。在一定程度上說明了前文分析結(jié)果的穩(wěn)健性。

4.4 斜率分析

為進(jìn)一步探討在種植規(guī)模對(duì)小農(nóng)戶耕地質(zhì)量保護(hù)行為的影響過程中，耕地細(xì)碎化程度的調(diào)節(jié)效應(yīng)，本文采用斜率分析法將耕地細(xì)碎化程度低于均值的1個(gè)標(biāo)準(zhǔn)差（-1 SD）和高于均值的1個(gè)標(biāo)準(zhǔn)差（+1 SD）分別定義為耕地細(xì)碎化程度低和耕地細(xì)碎化程度高，分別對(duì)小農(nóng)戶的農(nóng)藥施用行為、肥料施用行為和整體行為進(jìn)行分析。同時(shí)，為保證結(jié)果的穩(wěn)健性，本文分別以傾向得分匹配前和匹配后兩組數(shù)據(jù)為基礎(chǔ)進(jìn)行檢驗(yàn)，結(jié)果如表6、圖3和圖4所示。

表6 傾向得分匹配前后斜率分析結(jié)果Tab.6 Slope analysis results before and after propensity score matching

圖3 傾向得分匹配前斜率分析結(jié)果Fig.3 Results of slope analysis before propensity score matching

圖4 傾向得分匹配后斜率分析結(jié)果Fig.4 Results of slope analysis after propensity score matching

由表6、圖3和圖4分析可知，從農(nóng)藥施用行為來看，無論耕地細(xì)碎化程度低或高，種植規(guī)模對(duì)小農(nóng)戶農(nóng)藥施用行為均有顯著正向影響，但相比之下，在耕地細(xì)碎化程度低時(shí)，種植規(guī)模對(duì)小農(nóng)戶農(nóng)藥施用行為的影響更大；就肥料施用行為而言，耕地細(xì)碎化程度低時(shí)，種植規(guī)模的擴(kuò)大促進(jìn)了小農(nóng)戶合理的肥料施用行為，但當(dāng)耕地細(xì)碎化程度高時(shí)，種植規(guī)模的擴(kuò)大對(duì)小農(nóng)戶肥料施用行為無影響；從整體行為來看，無論耕地細(xì)碎化程度低或高，種植規(guī)模的擴(kuò)大均促進(jìn)了小農(nóng)戶耕地質(zhì)量保護(hù)行為的提升，但相比之下，在耕地細(xì)碎化程度低時(shí)，種植規(guī)模擴(kuò)大對(duì)小農(nóng)戶耕地質(zhì)量保護(hù)行為的促進(jìn)作用更為明顯。上述結(jié)果無論以匹配前的原始樣本或匹配后的新樣本為基礎(chǔ)，均成立，進(jìn)一步表明隨著耕地細(xì)碎化程度的提高，種植規(guī)模擴(kuò)大對(duì)小農(nóng)戶合理施用農(nóng)藥、肥料等耕地質(zhì)量保護(hù)行為的促進(jìn)作用降低，同時(shí)再次說明前文分析結(jié)果的穩(wěn)健性。

5 結(jié)論與政策啟示

本文以小農(nóng)戶為研究對(duì)象，在構(gòu)建種植規(guī)模、耕地細(xì)碎化和小農(nóng)戶耕地質(zhì)量保護(hù)行為理論關(guān)系的基礎(chǔ)上，以小農(nóng)戶農(nóng)藥、化肥施用行為為例，基于山東、河南和甘肅3省1 448份微觀農(nóng)戶調(diào)查數(shù)據(jù)，采用貝葉斯似不相關(guān)回歸模型、基于遺傳算法的傾向得分匹配模型和斜率分析模型等驗(yàn)證本文所提出的研究假設(shè)，主要得出如下結(jié)論：（1）種植規(guī)模的擴(kuò)大促進(jìn)了小農(nóng)戶耕地質(zhì)量保護(hù)行為，種植規(guī)模擴(kuò)大對(duì)農(nóng)戶耕地質(zhì)量保護(hù)行為的促進(jìn)作用在小農(nóng)戶層面存在。（2）在當(dāng)前傳統(tǒng)農(nóng)業(yè)向現(xiàn)代農(nóng)業(yè)過渡過程中，耕地細(xì)碎化對(duì)小農(nóng)戶農(nóng)業(yè)生產(chǎn)的影響明顯，在增加小農(nóng)戶農(nóng)業(yè)生產(chǎn)成本的同時(shí)，降低了小農(nóng)戶合理施藥與施肥等耕地質(zhì)量保護(hù)行為的傾向。（3）種植規(guī)模擴(kuò)大對(duì)小農(nóng)戶耕地質(zhì)量保護(hù)行為有促進(jìn)作用，特別是耕地細(xì)碎化程度低時(shí)，該促進(jìn)作用更為明顯，但也因當(dāng)前耕地細(xì)碎化較為普遍的現(xiàn)實(shí)而趨于弱化。

基于上述分析，本文提出：（1）在不違反法律法規(guī)的前提下，合理推進(jìn)小農(nóng)戶通過土地互換等方式擴(kuò)大單塊耕地種植規(guī)模，降低耕地細(xì)碎化程度，從而提升農(nóng)業(yè)生產(chǎn)便利性，降低生產(chǎn)成本，促進(jìn)小農(nóng)戶耕地質(zhì)量保護(hù)行為。（2）在當(dāng)前城市經(jīng)濟(jì)和就業(yè)無法迅速容納中國(guó)小農(nóng)的現(xiàn)實(shí)背景下，穩(wěn)定農(nóng)村土地承包關(guān)系長(zhǎng)久不變，有序推動(dòng)土地承包經(jīng)營(yíng)權(quán)流轉(zhuǎn)，完善適度規(guī)模經(jīng)營(yíng)，在使小農(nóng)戶進(jìn)城與返鄉(xiāng)實(shí)現(xiàn)“進(jìn)可攻、退可守”的同時(shí)，促進(jìn)小農(nóng)戶保護(hù)耕地質(zhì)量。（3）在小農(nóng)經(jīng)濟(jì)仍將承載中國(guó)農(nóng)業(yè)生產(chǎn)主體地位的形勢(shì)下，耕地質(zhì)量保護(hù)的政策作用對(duì)象在考慮規(guī)模種植戶的基礎(chǔ)上，也應(yīng)著眼于小農(nóng)戶，不可忽視小農(nóng)戶耕地質(zhì)量保護(hù)的主體性，從而改善耕地污染。