設(shè)施蔬菜農(nóng)戶病蟲害綜合防控行為風(fēng)險(xiǎn)評估*
——基于貝葉斯分類統(tǒng)計(jì)方法的實(shí)證分析

田家榛，孫煒琳

(中國農(nóng)業(yè)科學(xué)院農(nóng)業(yè)經(jīng)濟(jì)與發(fā)展研究所，北京 100081)

0 引言

病蟲害是蔬菜生產(chǎn)過程中最主要的生物災(zāi)害，對蔬菜的品質(zhì)和數(shù)量造成直接影響[1]。作為農(nóng)業(yè)生產(chǎn)發(fā)展中一種先進(jìn)的生產(chǎn)方式，設(shè)施農(nóng)業(yè)存在著違背自然規(guī)律的作業(yè)形式，制約設(shè)施農(nóng)業(yè)發(fā)揮優(yōu)勢的最大因素就是病蟲害問題[2]。由于溫室的生產(chǎn)環(huán)境和反季節(jié)的生產(chǎn)方式增加了防治的難度，致使設(shè)施農(nóng)業(yè)本身就為病蟲害的生存創(chuàng)造了適宜的條件[3]。目前依靠單一化學(xué)藥劑防治蔬菜病蟲害效率低下，防治效果欠佳[4]。農(nóng)作物病蟲害專業(yè)化綜合防治與產(chǎn)業(yè)化服務(wù)是現(xiàn)代農(nóng)業(yè)的重要標(biāo)志，也是蔬菜病蟲害綠色防控技術(shù)試驗(yàn)研究基礎(chǔ)上，現(xiàn)代化植保積極發(fā)展的方向。因此，設(shè)施蔬菜病蟲害的防控應(yīng)從保護(hù)菜田生態(tài)和考慮經(jīng)濟(jì)效益的角度出發(fā)，推行以農(nóng)業(yè)防治為基礎(chǔ)，生物、物理、化學(xué)防治密切配合的綜合措施，把化學(xué)農(nóng)藥使用量控制在國家最低標(biāo)準(zhǔn)以下。近年來，農(nóng)業(yè)農(nóng)村部積極引導(dǎo)各級農(nóng)業(yè)部門加大農(nóng)作物病蟲害綜合防控技術(shù)，尤其是綠色防控技術(shù)研發(fā)、集成、示范和推廣力度。像殺蟲燈、性誘劑等病蟲害綠色防控技術(shù)得到了大面積的推廣應(yīng)用。

設(shè)施蔬菜種植區(qū)域非化學(xué)藥劑防控行為的缺失和農(nóng)戶不規(guī)范的農(nóng)藥施用行為引發(fā)了農(nóng)產(chǎn)品質(zhì)量安全的風(fēng)險(xiǎn)，這種狀況較為普遍且短時(shí)間很難改變。在我國，設(shè)施農(nóng)業(yè)生產(chǎn)者缺乏農(nóng)藥施用的正確認(rèn)知、過量施用農(nóng)藥、施用高毒高殘留農(nóng)藥、不按照使用說明用藥、藥械不能及時(shí)維護(hù)這些現(xiàn)象十分普遍[5]。高效防治手段的實(shí)施可以替代化學(xué)藥劑的病蟲害防治效果，有效降低病蟲害抗藥性，保障蔬菜品質(zhì)和安全[6]。根據(jù)設(shè)施蔬菜生長特性和栽培要求，剔除化學(xué)藥劑防治這種傳統(tǒng)低利用率的防治方式[4]，文章選取代表生物防治和物理防治的生物農(nóng)藥和綠色防控措施，以及代表栽培防治的土壤消毒和污染修復(fù)?？紤]到蔬菜品種的多樣性以及農(nóng)戶層面的防治問題，該文并未將品種防治和行政上的植物檢疫納入研究的綜合防治體系中[7]。目前來看，推廣綜合防控體系僅靠政府自上而下的監(jiān)管很難讓農(nóng)戶做出實(shí)質(zhì)性改變，應(yīng)該落實(shí)在微觀農(nóng)戶層面進(jìn)行推廣[8]。

該文以山東、河北、浙江3省設(shè)施蔬菜種植區(qū)域的調(diào)研數(shù)據(jù)為基礎(chǔ)，以設(shè)施蔬菜種植戶的生產(chǎn)性特征作為風(fēng)險(xiǎn)因素進(jìn)行評估，構(gòu)建樸素的貝葉斯分類模型并深入分析影響農(nóng)戶病蟲害綜合防控行為的重要因素。為設(shè)施蔬菜產(chǎn)業(yè)良好有序的發(fā)展，推動蔬菜病蟲害專業(yè)化防治服務(wù)進(jìn)程提供政策啟示。

1 相關(guān)文獻(xiàn)梳理及其評述

1.1 農(nóng)戶施用生物農(nóng)藥行為的缺失風(fēng)險(xiǎn)評估

生物農(nóng)藥與化學(xué)藥劑相比在推廣過程中本身固有一定的缺陷，另外，農(nóng)戶家庭收入、人均耕地、種植規(guī)模與年限、政府補(bǔ)貼、生產(chǎn)合作社成員資格等都是形成生產(chǎn)者這類防控行為缺失風(fēng)險(xiǎn)的關(guān)鍵因素。朱淀[9]等研究結(jié)果顯示，具有家庭蔬菜種植規(guī)模在0.05hm2以下、非農(nóng)業(yè)性收入在3萬及以上、參加農(nóng)業(yè)生產(chǎn)合作組織等生產(chǎn)性特征的蔬菜種植戶更傾向于施用生物農(nóng)藥，而政府補(bǔ)貼對蔬菜種植戶施用生物農(nóng)藥行為意愿的驅(qū)動作用不明顯。傅新紅[10]等總結(jié)農(nóng)戶施用生物農(nóng)藥的優(yōu)點(diǎn)為低毒、良好的環(huán)境兼容性以及消除病蟲害的抗藥性。在農(nóng)戶收入無法達(dá)到預(yù)期的條件下，生物農(nóng)藥的價(jià)格過高會致使施用生物農(nóng)藥行為的缺失。張?jiān)迫A[11]等指出人均耕地面積越大說明農(nóng)戶勞動力強(qiáng)度越高，更傾向于施用無公害農(nóng)藥和綠色農(nóng)藥。農(nóng)戶與涉農(nóng)企業(yè)、農(nóng)業(yè)技術(shù)協(xié)會聯(lián)系的緊密程度也是評估這一防控行為缺失風(fēng)險(xiǎn)的重要因素。OleSkovmand[12]指出生物農(nóng)藥市場競爭力弱的原因在于其購買成本、質(zhì)量和性能。馬永生[13]等分析指出政府層面的宣傳力度不夠以及無公害農(nóng)產(chǎn)品認(rèn)證制度的缺乏是影響生物農(nóng)藥推廣的重要因素。

1.2 農(nóng)戶土壤消毒和污染修復(fù)行為的缺失風(fēng)險(xiǎn)評估

定期的土壤消毒可以鈍化、除去土壤中的微生物和毒素，除了常見的藥劑處理之外，利用日光、生石灰做消毒處理是值得提倡的栽培防治方式[14]。與農(nóng)戶土壤消毒處理的影響因素有關(guān)的文獻(xiàn)基本是從自然學(xué)科的角度去研究，涉及到農(nóng)戶行為理論和特征因素的文獻(xiàn)極少。楊志海[15]等分析得出農(nóng)戶耕地塊數(shù)、人均耕地面積、是否有租入土地、兼業(yè)與家庭非農(nóng)就業(yè)比例都是造成土壤保護(hù)行為缺失風(fēng)險(xiǎn)的顯著性因素。鄭紀(jì)芳[16]等研究表明農(nóng)戶改良土壤的意愿比較強(qiáng)烈，但具體土壤保護(hù)行為的落實(shí)不到位，土壤處理行為缺失帶來有毒物質(zhì)殘留和酸堿度不平衡等風(fēng)險(xiǎn)。劉洪彬[17]等研究表明：土地利用方式、程度和投入強(qiáng)度的不同會獲得不同質(zhì)量土壤，土壤有機(jī)質(zhì)的含量差異性明顯； 另外，農(nóng)戶兼業(yè)化程度不同會對農(nóng)戶的土壤處理和土地利用方式產(chǎn)生重要影響。盧彥琦[18]等認(rèn)為農(nóng)戶純收入是土壤污染防治支付意愿的主要因素。

1.3 綠色防控行為的缺失風(fēng)險(xiǎn)評估

收入結(jié)構(gòu)、耕地類型、種植規(guī)模、土地細(xì)碎化、政府財(cái)政補(bǔ)貼等都是形成生產(chǎn)者這類防控行為缺失風(fēng)險(xiǎn)的關(guān)鍵因素。徐小俊[19]等認(rèn)為良好防控的效果和低成本是海南農(nóng)戶采納綠色防控措施的主要條件。在綠色防控技術(shù)推行方面，多數(shù)農(nóng)戶希望加大技術(shù)和資金的支持力度。Bonabana-WabbiJ[20]提出影響生產(chǎn)者綠色防控行為缺失的最主要的生產(chǎn)性特征因素為種植規(guī)模。耿宇寧[21]等分析了經(jīng)濟(jì)激勵對果農(nóng)綠色防控技術(shù)采納的影響，市場激勵機(jī)制的作用明顯大于政府激勵機(jī)制的作用。另外研究還得出當(dāng)前土地細(xì)碎化制約了需要大面積連片采用的物理防治技術(shù)的推廣。Pietola[22]同樣提出財(cái)政補(bǔ)貼在農(nóng)戶采納有機(jī)技術(shù)方面有著重要的促進(jìn)作用。此外，SherlundSM[23]研究指出直接影響到農(nóng)戶病蟲害綠色防控技術(shù)采納的成本收益、風(fēng)險(xiǎn)感知、資源獲取和勞力需求，也是病蟲害綠色防控行為缺失的重要風(fēng)險(xiǎn)因素。

由于數(shù)據(jù)獲取的局限性、特質(zhì)性等因素的影響，農(nóng)業(yè)政策和農(nóng)戶偏好對綜合防控行為影響的相關(guān)研究較多。而從農(nóng)戶種植環(huán)節(jié)入手，將農(nóng)戶的生產(chǎn)性特征作為風(fēng)險(xiǎn)評估的因素涉及到的文獻(xiàn)不多。相關(guān)文獻(xiàn)研究中變量的選取大多集中體現(xiàn)在單一生產(chǎn)性特征因素與單一農(nóng)戶防控行為的關(guān)系。該文研究對象選取設(shè)施蔬菜種植戶這類采用先進(jìn)生產(chǎn)方式的生產(chǎn)者群體，根據(jù)設(shè)施蔬菜生長特性和栽培要求將病蟲害高效防控行為進(jìn)行總結(jié)并納入綜合防治體系中，深入剖析農(nóng)戶生產(chǎn)性特征因素對農(nóng)戶病蟲害綜合防控體系的影響機(jī)理并進(jìn)行系統(tǒng)完整的比對分析。探究制約高效病蟲害防控技術(shù)在生產(chǎn)者群體中推廣的因素，以期為政府實(shí)行病蟲害專業(yè)化防治、保障蔬菜質(zhì)量提供一定的理論借鑒。

2 問卷設(shè)計(jì)和數(shù)據(jù)來源

2.1 樣本選擇

為了保證樣本的全面性和準(zhǔn)確性，選取的調(diào)研對象包括散戶和種植大戶且都有1年以上大棚蔬菜種植經(jīng)驗(yàn)。在選擇調(diào)研區(qū)域時(shí)，主要采取的是判斷性抽樣調(diào)查的方法，調(diào)查區(qū)域都是大棚蔬菜種植戶的集聚區(qū)域，即選取的是河北省張家口市、浙江省嘉興市、山東省濰坊市設(shè)施蔬菜種植戶集聚區(qū)域。為保證樣本的合理性，根據(jù)當(dāng)?shù)亟?jīng)濟(jì)發(fā)展?fàn)顩r和地形地貌特征對3個(gè)省采取分層抽樣與隨機(jī)抽樣相結(jié)合的抽樣方式。每個(gè)地級市或縣下面選取5個(gè)鎮(zhèn)，每個(gè)鎮(zhèn)下面選取全部有代表性的行政村，每個(gè)行政村選取具有設(shè)施蔬菜種植經(jīng)驗(yàn)的農(nóng)戶，盡量保證各個(gè)行政村抽取的樣本數(shù)量相等。

2.2 問卷內(nèi)容和調(diào)研方式

調(diào)查問卷都是在咨詢自然學(xué)科技術(shù)專家和充分閱讀文獻(xiàn)的基礎(chǔ)上設(shè)計(jì)的，考慮設(shè)計(jì)指標(biāo)的可操作性和量化性，內(nèi)容涉及到設(shè)施蔬菜種植戶的病蟲害防控行為和行為的風(fēng)險(xiǎn)因素。調(diào)研對象均全程參與蔬菜種植環(huán)節(jié)，調(diào)研工作都是在一對一情況下進(jìn)行農(nóng)戶訪談，問卷由調(diào)研人員填寫，這樣可避免農(nóng)戶文化層次和理解能力的偏差對問卷的真實(shí)性和有效性的影響。調(diào)研工作在2018年3月結(jié)束，累計(jì)發(fā)放問卷的總數(shù)為734份，其中回收有效問卷為721份，有效回收率為98%。問卷的主要內(nèi)容包括農(nóng)戶的基本特征、具體防控行為和意愿等問題。

3 方法選擇和模型設(shè)定

3.1 方法選擇

貝葉斯分類器是各種分類器中分類錯(cuò)誤概率最小或者在預(yù)先給定代價(jià)的情況下平均風(fēng)險(xiǎn)最小的分類器，設(shè)計(jì)方法是一種最基本的統(tǒng)計(jì)分類方法。其分類原理是通過某對象的先驗(yàn)概率，利用貝葉斯公式計(jì)算出其后驗(yàn)概率，即該對象屬于某一類的概率，選擇具有最大后驗(yàn)概率的類作為該對象所屬的類。

該文選取樸素的貝葉斯分類器對設(shè)施蔬菜種植戶病蟲害綜合防控行為缺失風(fēng)險(xiǎn)進(jìn)行評估，該分類器使用廣泛、結(jié)構(gòu)簡單、處理數(shù)據(jù)高效。樸素貝葉斯分類器進(jìn)行分析時(shí)具備較高的分類準(zhǔn)確度，保持最大概率值的特征。貝葉斯分類器主要描述了設(shè)施蔬菜種植戶病蟲害綜合防控行為的先驗(yàn)概率、條件概率和后驗(yàn)概率之間的關(guān)系。

3.2 模型設(shè)定

貝葉斯定理的概述是已知條件概率如何得到兩個(gè)事件交換后的概率，即已知P(A/B)的前提下求解P(B/A)。事件A代表設(shè)施蔬菜種植戶防控風(fēng)險(xiǎn)行為，B代表影響設(shè)施蔬菜種植戶防控行為缺失風(fēng)險(xiǎn)的因素，主要選取設(shè)施蔬菜種植戶生產(chǎn)環(huán)節(jié)的特征。該文的關(guān)注點(diǎn)在P(A/B)這一條件概率的數(shù)值，即在B發(fā)生條件下A的概率，條件概率求解公式即貝葉斯定理公式如(1)式所示：

(1)

在設(shè)施蔬菜種植戶防控行為數(shù)據(jù)分類的應(yīng)用中，P(A)是作為設(shè)施蔬菜種植戶防控行為的先驗(yàn)概率；P(B)是設(shè)施蔬菜種植戶防控行為調(diào)研數(shù)據(jù)的基礎(chǔ)概率，即歸一化因子；P(B/A)是作為沒有分類設(shè)施蔬菜種植戶防控行為調(diào)研數(shù)據(jù)B的似然函數(shù)；P(A/B)稱作后驗(yàn)概率，即B發(fā)生條件下A發(fā)生的概率?？偟脕碚f貝葉斯公式是通過對未分類設(shè)施蔬菜種植戶防控行為數(shù)據(jù)的條件概率分布進(jìn)行計(jì)算，實(shí)現(xiàn)對設(shè)施蔬菜種植戶病蟲害防控行為缺失風(fēng)險(xiǎn)評估數(shù)據(jù)的貝葉斯統(tǒng)計(jì)推斷。設(shè)施蔬菜種植戶防控行為缺失風(fēng)險(xiǎn)評估的貝葉斯分類定理公式如(2)式所示：

(2)

i=1， 2，…， 8分別表示設(shè)施蔬菜種植戶防控行為的各類風(fēng)險(xiǎn)因素，即農(nóng)戶生產(chǎn)性特征。具體事件主要有種植規(guī)模和種植年限，勞動力強(qiáng)度，生產(chǎn)經(jīng)營類型，生產(chǎn)過程中補(bǔ)貼、投保、蔬菜認(rèn)證情況，生產(chǎn)合作社成員資格。j=1時(shí)，表示設(shè)施蔬菜種植戶綠色防控措施缺失情況的風(fēng)險(xiǎn)因素評估；j=2時(shí)，表示設(shè)施蔬菜種植戶土壤消毒或污染修復(fù)缺失的風(fēng)險(xiǎn)因素評估；j=3時(shí)，表示設(shè)施蔬菜種植戶施用生物農(nóng)藥缺失情況的風(fēng)險(xiǎn)因素評估。

4 實(shí)證分析

4.1 防控行為及其風(fēng)險(xiǎn)因素的統(tǒng)計(jì)分析

表1 設(shè)施蔬菜種植戶防控行為風(fēng)險(xiǎn)因素及評估風(fēng)險(xiǎn)行為的參數(shù)

事件標(biāo)號變量名稱參數(shù)標(biāo)號評估參數(shù)頻數(shù)頻率(%)B1種植年限B115年以下12417.20B125～10年26436.62B1310～15年7610.54B1415年以上25735.64B2種植規(guī)模B210.01hm2以下7410.26B220.011～0.03hm231043.00B230.031～0.05hm215321.22B240.05hm2以上18425.52B3勞動力強(qiáng)度B310～20131.80B3221%～40%12717.62B3341%～60%21129.26B3461%～80%19226.63B3581%～100%17824.69B4農(nóng)業(yè)生產(chǎn)經(jīng)營類型B41純設(shè)施蔬菜種植戶56878.78B42有其他收入來源15321.22B5農(nóng)業(yè)生產(chǎn)的政府補(bǔ)貼B51有補(bǔ)貼24834.40B52無補(bǔ)貼47365.60B6農(nóng)業(yè)生產(chǎn)的風(fēng)險(xiǎn)投保B61有投保9613.31B62無投保62586.69B7農(nóng)產(chǎn)品認(rèn)證B71有認(rèn)證22030.51B72無認(rèn)證50169.49B8生產(chǎn)專業(yè)組織的成員資格B81加入20328.16B82沒有加入51871.84

表1表示從農(nóng)戶生產(chǎn)性特征的角度出發(fā)，對設(shè)施蔬菜種植戶防控行為缺失風(fēng)險(xiǎn)的影響因素進(jìn)行統(tǒng)計(jì)的結(jié)果。從農(nóng)戶的生產(chǎn)性特征來看，設(shè)施蔬菜種植戶的種植年限集中在5～10年和15年以上這兩個(gè)時(shí)間段，生產(chǎn)規(guī)模集中在0.011～0.03hm2這一分段。該文衡量農(nóng)業(yè)生產(chǎn)勞動力強(qiáng)度的指標(biāo)是指從事農(nóng)業(yè)生產(chǎn)的人口占家庭總?cè)丝诘谋戎?，比重越小，代表農(nóng)業(yè)勞動力強(qiáng)度越大。除去20%以下的分段，樣本的農(nóng)業(yè)勞動力強(qiáng)度在劃分的4個(gè)分段中大致平均分布。農(nóng)戶生產(chǎn)經(jīng)營兩種類型所占比例相差較大。

在設(shè)施生產(chǎn)過程中，不同地區(qū)設(shè)施蔬菜種植戶面臨的政策環(huán)境不同，農(nóng)業(yè)生產(chǎn)上的補(bǔ)貼類型也是不同的。統(tǒng)計(jì)結(jié)果顯示，約65%的農(nóng)戶在設(shè)施生產(chǎn)過程中沒有政府補(bǔ)貼。農(nóng)業(yè)生產(chǎn)上的風(fēng)險(xiǎn)投保主要涵蓋蔬菜價(jià)格和災(zāi)害方面的投保，有85%以上的農(nóng)戶缺乏生產(chǎn)上的投保行為，總體上風(fēng)險(xiǎn)投保的意識不強(qiáng)。農(nóng)產(chǎn)品認(rèn)證指農(nóng)戶生產(chǎn)的設(shè)施蔬菜是否經(jīng)過有關(guān)權(quán)威機(jī)構(gòu)的有機(jī)認(rèn)證和綠色認(rèn)證，經(jīng)過認(rèn)證的農(nóng)戶僅占比30%。統(tǒng)計(jì)結(jié)果顯示，有70%以上的農(nóng)戶沒有加入當(dāng)?shù)氐纳a(chǎn)專業(yè)合作組織，說明設(shè)施蔬菜種植區(qū)域相關(guān)合作經(jīng)營組織的建設(shè)機(jī)制不健全，對種植戶生產(chǎn)行為缺乏有效的規(guī)整。

該文運(yùn)用MATLAB軟件對721個(gè)種植戶樣本數(shù)據(jù)進(jìn)行分析，計(jì)算得到先驗(yàn)概率和條件概率，運(yùn)用樣本數(shù)據(jù)構(gòu)建出貝葉斯分類評估模型并進(jìn)行訓(xùn)練，對訓(xùn)練后的評估模型進(jìn)行驗(yàn)證，計(jì)算后驗(yàn)概率并對結(jié)果進(jìn)行對比分析。

表2 設(shè)施蔬菜種植戶防控風(fēng)險(xiǎn)行為的先驗(yàn)概率

事件標(biāo)號變量名稱參數(shù)標(biāo)號評估參數(shù)頻率(%)A1是否采用綠色防控措施A11經(jīng)常使用5.13A12有時(shí)使用23.44A13很少使用35.78A14不使用35.65A2是否對土壤做消毒處理或污染修復(fù)A21是71.98A22否28.02A3設(shè)施蔬菜上施用生物農(nóng)藥的比重A31一半以上29.13A32一半以下52.43A33幾乎沒有18.44

該文根據(jù)設(shè)施蔬菜的栽培要求和生長特性，選取的3個(gè)指標(biāo)較充分地概括了設(shè)施蔬菜種植戶的病蟲害綜合防控風(fēng)險(xiǎn)行為。問卷涉及防蟲網(wǎng)、黃板、殺蟲燈等高效綠色的防控措施和阿維菌素、多抗霉素、農(nóng)用鏈霉菌等常見生物農(nóng)藥供農(nóng)戶選擇，土壤的處理包括消毒和污染修復(fù)處理。從農(nóng)戶病蟲害防控風(fēng)險(xiǎn)行為來看，80%以上的農(nóng)戶在設(shè)施蔬菜種植過程中施用生物農(nóng)藥， 但生物農(nóng)藥的施用在全部農(nóng)藥中的占比并不高。在生產(chǎn)過程中采用綠色防控措施和土壤處理防治病蟲害的農(nóng)戶在70%左右(表2)。

4.2 條件概率分析

表3計(jì)算的是生產(chǎn)者在采用綠色防控措施頻率不同的條件下，設(shè)施蔬菜種植戶防控行為缺失的風(fēng)險(xiǎn)因素概率分布，根據(jù)種植戶采取綠色防控措施的行為頻率，劃分為4種行為表現(xiàn)：即經(jīng)常使用(A11)、有時(shí)使用(A12)、很少使用(A13)、不使用(A14)。分析結(jié)果顯示，綠色防控措施采用頻率較高的農(nóng)戶多為種植經(jīng)驗(yàn)豐富、勞動力強(qiáng)度在20% ～60%的農(nóng)業(yè)生產(chǎn)者； 在采用綠色防控措施的頻率較高的條件下，農(nóng)戶在生產(chǎn)過程中傾向于風(fēng)險(xiǎn)投保，生產(chǎn)的農(nóng)產(chǎn)品大多經(jīng)過權(quán)威認(rèn)證； 在采用綠色防控措施的頻率較高的條件下，農(nóng)業(yè)生產(chǎn)者一般具有生產(chǎn)專業(yè)組織成員資格； 農(nóng)戶這類防控行為在不同的生產(chǎn)經(jīng)營類型和是否有生產(chǎn)性補(bǔ)貼的情況下差異不明顯。

表3 生產(chǎn)者在采用綠色防控措施條件下風(fēng)險(xiǎn)因素的條件概率

變量評估參數(shù)參數(shù)編號經(jīng)常使用P(Bi/A11)有時(shí)使用P(Bi/A12)很少使用P(Bi/A13)不使用P(Bi/A14)種植年限5年以下B110.1620.0890.1550.2455～10年B120.1890.3310.3720.40910～15年B120.1080.1240.0970.10115年以上B140.5410.4560.3760.245種植規(guī)模0.01hm2以下B210.1350.1010.1050.0970.011～0.03hm2B220.4850.5680.4460.3150.031～0.05hm2B230.1350.1660.2130.2530.05hm2以上B240.2430.1660.2360.335勞動力強(qiáng)度0～20%B3100.0120.0120.03121%～40%B320.2160.1950.1740.16041%～60%B330.3510.3080.2910.27661%～80%B340.2970.3020.2980.20681%～100%B350.1350.1830.2250.327農(nóng)業(yè)生產(chǎn)經(jīng)營類型純設(shè)施蔬菜種植戶B410.7030.7810.7950.798有其他收入來源B420.2970.2190.2050.202農(nóng)業(yè)生產(chǎn)的政府補(bǔ)貼有補(bǔ)貼B510.2430.2370.3330.440無補(bǔ)貼B520.7570.7630.6670.560農(nóng)業(yè)生產(chǎn)的風(fēng)險(xiǎn)投保有投保B610.3240.1780.1320.078無投保B620.6760.8220.8680.922農(nóng)產(chǎn)品認(rèn)證有認(rèn)證B710.5950.3200.3060.253無認(rèn)證B720.4050.6800.6940.747生產(chǎn)專業(yè)組織的成員資格加入B810.4050.3430.2790.226沒有加入B820.5950.6570.7210.774

表4 生產(chǎn)者在土壤消毒和污染修復(fù)的條件下風(fēng)險(xiǎn)因素的條件概率

變量評估參數(shù)參數(shù)編號進(jìn)行處理P(Bi/A21)不作處理P(Bi/A22)種植年限5年以下B110.1980.1045～10年B120.4030.27210～15年B120.0890.14915年以上B140.310.475種植規(guī)模0.01hm2以下B210.0910.1340.011～0.03hm2B220.3740.5740.031～0.05hm2B230.2240.1830.05hm2以上B240.3120.109勞動力強(qiáng)度0～20%B310.0150.02521%～40%B320.1540.23341%～60%B330.2930.29261%～80%B340.2580.28781%～100%B350.2790.263農(nóng)業(yè)生產(chǎn)經(jīng)營類型純設(shè)施蔬菜種植戶B410.8090.733有其他收入來源B420.1910.267農(nóng)業(yè)生產(chǎn)的政府補(bǔ)貼有補(bǔ)貼B510.430.124無補(bǔ)貼B520.570.876農(nóng)業(yè)生產(chǎn)的風(fēng)險(xiǎn)投保有投保B610.1370.124無投保B620.8630.876農(nóng)產(chǎn)品認(rèn)證有認(rèn)證B710.3060.302無認(rèn)證B720.6940.698生產(chǎn)專業(yè)組織的成員資格加入B810.3160.193沒有加入B820.6840.807

表4計(jì)算的是生產(chǎn)者在土壤處理與否的條件下，設(shè)施蔬菜種植戶防控行為缺失的風(fēng)險(xiǎn)因素概率分布，表中的概率數(shù)據(jù)反映的是進(jìn)行土壤消毒和污染修復(fù)的設(shè)施生產(chǎn)者群體與不進(jìn)行土壤處理的設(shè)施生產(chǎn)者群體之間的生產(chǎn)性特征差異。分析結(jié)果顯示，農(nóng)業(yè)生產(chǎn)者在進(jìn)行土壤消毒和污染修復(fù)與否的條件下，種植規(guī)模、生產(chǎn)性政府補(bǔ)貼和生產(chǎn)的專業(yè)組織成員資格等生產(chǎn)性風(fēng)險(xiǎn)因素的條件概率差異較大。設(shè)施蔬菜種植規(guī)模大于0.05hm2的種植大戶傾向于采取這種防控行為，而蔬菜種植規(guī)模不足0.03hm2的散戶則極少通過土壤處理來預(yù)防病蟲害。

表5計(jì)算的是生產(chǎn)者在生物農(nóng)藥施用比例不同的條件下，設(shè)施蔬菜種植戶防控行為缺失的風(fēng)險(xiǎn)因素概率分布。根據(jù)種植戶在作物一個(gè)完整生長周期內(nèi)施用的生物農(nóng)藥占全部農(nóng)藥的比重，劃分為3種行為表現(xiàn)：即比重在50%以上(A31)、比重在50%以下(A32)、比重接近于0(A33)。表中的概率數(shù)據(jù)反映的是生物農(nóng)藥施用比重不同的農(nóng)業(yè)設(shè)施生產(chǎn)者群體之間的生產(chǎn)性特征差異。分析結(jié)果顯示，在農(nóng)業(yè)生產(chǎn)者施用生物農(nóng)藥占全部農(nóng)藥比重不同的條件下，生產(chǎn)者種植年限、種植規(guī)模、農(nóng)業(yè)生產(chǎn)經(jīng)營類型、農(nóng)業(yè)生產(chǎn)的政府補(bǔ)貼等生產(chǎn)性風(fēng)險(xiǎn)因素的條件概率波動都相對較大。生物農(nóng)藥施用占比較高的農(nóng)業(yè)生產(chǎn)者家中勞動力人口占總?cè)丝诘谋戎丶性?1% ～60%，可見生物農(nóng)藥施用對勞動力的替代作用是可觀的。農(nóng)業(yè)生產(chǎn)者施用生物農(nóng)藥占全部農(nóng)藥比重在50%以下的條件下，生產(chǎn)出的農(nóng)產(chǎn)品經(jīng)過認(rèn)證的概率不足30%； 農(nóng)業(yè)生產(chǎn)者施用生物農(nóng)藥占全部農(nóng)藥比重在50%以上的條件下，生產(chǎn)出的農(nóng)產(chǎn)品經(jīng)過認(rèn)證的概率在30%以上?？梢?，生物農(nóng)藥施用是保障農(nóng)產(chǎn)品質(zhì)量安全的重要環(huán)節(jié)。

表5 生產(chǎn)者在生物農(nóng)藥施用的條件下風(fēng)險(xiǎn)因素的條件概率

變量評估參數(shù)參數(shù)編號一半以上P(Bi/A31)一半以下P(Bi/A32)幾乎沒有P(Bi/A33)種植年限5年以下B110.1330.1930.1735～10年B120.3140.4580.18810～15年B120.1050.0950.13515年以上B140.4480.2540.504種植規(guī)模0.01hm2以下B210.1190.0630.1880.011～0.03hm2B220.4480.3940.5040.031～0.05hm2B230.2050.2410.1430.05hm2以上B240.2290.3020.165勞動力強(qiáng)度0～20%B310.0050.0210.03021%～40%B320.1620.1610.24141%～60%B330.3000.2990.26361%～80%B340.2900.2410.30181%～100%B350.2430.2780.165農(nóng)業(yè)生產(chǎn)經(jīng)營類型純設(shè)施蔬菜種植戶B410.8100.8310.632有其他收入來源B420.1900.1690.368農(nóng)業(yè)生產(chǎn)的政府補(bǔ)貼有補(bǔ)貼B510.3140.4370.128無補(bǔ)貼B520.6860.5630.872農(nóng)業(yè)生產(chǎn)的風(fēng)險(xiǎn)投保有投保B610.1520.1160.150無投保B620.8480.8840.850農(nóng)產(chǎn)品認(rèn)證有認(rèn)證B710.3760.2780.271無認(rèn)證B720.6240.7220.729生產(chǎn)專業(yè)組織的成員資格加入B810.2760.2940.256沒有加入B820.7240.7060.744

4.3 計(jì)算并比較后驗(yàn)概率

圖1表示在設(shè)施蔬菜種植戶病蟲害防控行為在生產(chǎn)性特征風(fēng)險(xiǎn)因素條件下的后驗(yàn)概率。橫坐標(biāo)具體含義為既定的10個(gè)生產(chǎn)性特征風(fēng)險(xiǎn)因素，縱坐標(biāo)的具體含義為在具備確定生產(chǎn)性特征的條件下，種植戶防控風(fēng)險(xiǎn)行為的后驗(yàn)概率值。確定的農(nóng)戶生產(chǎn)性特征條件分別是種植經(jīng)驗(yàn)豐富、種植大戶、高勞動力強(qiáng)度、純設(shè)施蔬菜種植戶、政府補(bǔ)貼、風(fēng)險(xiǎn)投保、經(jīng)過農(nóng)產(chǎn)品認(rèn)證和加入生產(chǎn)專業(yè)組織。

圖1 設(shè)施蔬菜種植戶病蟲害防控行為風(fēng)險(xiǎn)的后驗(yàn)概率分布

對比分析3類設(shè)施蔬菜種植戶病蟲害防控風(fēng)險(xiǎn)行為的后驗(yàn)概率圖，可以得出如下結(jié)論。

(1)在10個(gè)生產(chǎn)性特征風(fēng)險(xiǎn)因素已知條件下，農(nóng)業(yè)生產(chǎn)者的施用生物農(nóng)藥和進(jìn)行土壤處理的后驗(yàn)概率變化趨勢大致相同，說明兩者的生產(chǎn)性特征風(fēng)險(xiǎn)因素的影響機(jī)制相近。B1(種植年限)對于綠色防控風(fēng)險(xiǎn)行為的影響較大，對3種農(nóng)業(yè)生產(chǎn)者病蟲害防控風(fēng)險(xiǎn)行為的影響具有較小的差異性，因?yàn)榻?jīng)驗(yàn)豐富的農(nóng)業(yè)生產(chǎn)者會有意識地統(tǒng)籌安排各種防控手段并納入綜合管理體系。

(2)在10個(gè)生產(chǎn)性特征風(fēng)險(xiǎn)因素已知條件下，生產(chǎn)者施用生物農(nóng)藥的概率最高，在0.85的概率值上下波動。表示農(nóng)戶施用生物農(nóng)藥的防控行為比較普遍且易受到風(fēng)險(xiǎn)因素的影響，因此評估農(nóng)業(yè)生產(chǎn)者施用生物農(nóng)藥行為缺失的風(fēng)險(xiǎn)具有較好的現(xiàn)實(shí)意義?？傮w來說B6(風(fēng)險(xiǎn)投保)、B7(農(nóng)產(chǎn)品認(rèn)證)、B8(生產(chǎn)專業(yè)組織成員資格)等風(fēng)險(xiǎn)因素對生產(chǎn)者防控行為缺失影響的概率水平較穩(wěn)定，其中B6(風(fēng)險(xiǎn)投保)對3種農(nóng)業(yè)生產(chǎn)者病蟲害防控行為風(fēng)險(xiǎn)的影響具有較小的差異性。主要原因在于進(jìn)行風(fēng)險(xiǎn)投保的農(nóng)戶風(fēng)險(xiǎn)防范意識更好，對有利于減少病蟲害風(fēng)險(xiǎn)的各類防控措施接受意愿更高。農(nóng)產(chǎn)品的認(rèn)證和生產(chǎn)合作社的成員資格對農(nóng)戶種植環(huán)節(jié)提出更高標(biāo)準(zhǔn)的要求，規(guī)范農(nóng)戶的防控行為并促使其更高效綜合地防控病蟲害。

(3)B5(農(nóng)業(yè)生產(chǎn)的政府補(bǔ)貼)對3種農(nóng)業(yè)生產(chǎn)者病蟲害防控行為風(fēng)險(xiǎn)的影響具有較大的差異性。在諸多風(fēng)險(xiǎn)因素中，農(nóng)戶施用生物農(nóng)藥和土壤處理的防控行為受到政府生產(chǎn)性補(bǔ)貼影響的概率較高，而對于農(nóng)戶綠色防控風(fēng)險(xiǎn)行為的影響概率較低。政府生產(chǎn)性補(bǔ)貼對生產(chǎn)者病蟲害防控行為風(fēng)險(xiǎn)概率影響一致性較弱的原因比較復(fù)雜： ①一方面，政府生產(chǎn)性補(bǔ)貼可能會促使農(nóng)戶擴(kuò)大生產(chǎn)規(guī)模，采取更高效的防控行為。另一方面，生產(chǎn)規(guī)模的擴(kuò)大和生產(chǎn)效率的提高會促使政府加大補(bǔ)貼和獎勵力度； ②當(dāng)期補(bǔ)貼政策的實(shí)施具有一定的滯后效應(yīng)，對農(nóng)戶防控行為的改變不是立竿見影的。

(4)B2(種植規(guī)模)、B3(勞動力強(qiáng)度)、B4(生產(chǎn)經(jīng)營類型)對生產(chǎn)者3類防控行為風(fēng)險(xiǎn)的影響具有一定的差異性。種植大戶在設(shè)施蔬菜栽培環(huán)節(jié)會實(shí)現(xiàn)病蟲害通防統(tǒng)治，自負(fù)盈虧和承擔(dān)風(fēng)險(xiǎn)。他們按照基本種植技術(shù)要求規(guī)范耕種，注重土壤防治和生物農(nóng)藥防治。但是受限于種植面積，綠色防控措施的普及度不高。家庭從事農(nóng)業(yè)生產(chǎn)的勞動力人口占總?cè)丝诒戎刈畹蜁r(shí)，勞動力強(qiáng)度最大，此時(shí)農(nóng)戶生產(chǎn)蔬菜很少進(jìn)入市場，使用綠色防控措施減少產(chǎn)量損失風(fēng)險(xiǎn)的概率較低。純設(shè)施蔬菜經(jīng)營的農(nóng)戶收入來源單一，有更多的時(shí)間和精力從事農(nóng)業(yè)生產(chǎn)并更看重農(nóng)業(yè)生產(chǎn)的長期效益，統(tǒng)籌推行高效綜合防控的力度加強(qiáng)。

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

5.1 研究結(jié)論

該文對設(shè)施蔬菜典型種植區(qū)域721個(gè)生產(chǎn)者的樣本數(shù)據(jù)進(jìn)行分析，建立樸素的貝葉斯分類模型并進(jìn)行樣本的訓(xùn)練。選取有利于減少病蟲害風(fēng)險(xiǎn)的高效綜合防控行為，進(jìn)一步對設(shè)施蔬菜種植戶病蟲害綜合防控行為缺失進(jìn)行了風(fēng)險(xiǎn)評估。

結(jié)果顯示：農(nóng)戶施用生物農(nóng)藥的行為比較普遍且易受到風(fēng)險(xiǎn)因素的影響。近些年，設(shè)施蔬菜種植戶對生物農(nóng)藥的無害性、強(qiáng)選擇性等優(yōu)勢更加明晰，接受度也越來越高，從生產(chǎn)環(huán)節(jié)入手增加農(nóng)戶家全部農(nóng)藥中的生物農(nóng)藥比重是切實(shí)可行的藥劑防治手段。在現(xiàn)代農(nóng)業(yè)生產(chǎn)中，農(nóng)戶對藥劑防治有著高依賴性，一定程度上忽略了物理防治和生物防治等防治手段的重要性。設(shè)施蔬菜在栽培過程中容易受到土傳病害的影響，故調(diào)研區(qū)域農(nóng)戶土壤處理的意識普遍較強(qiáng)。綠色防控是通過生態(tài)控制、生物防治等環(huán)境友好型的措施來防治病蟲的有效行為，目標(biāo)是減少化學(xué)農(nóng)藥的用量，促進(jìn)農(nóng)產(chǎn)品安全生產(chǎn)。設(shè)施蔬菜種植戶采用綠色防控措施防控病蟲害的概率不高，下一步應(yīng)加強(qiáng)引導(dǎo)設(shè)施蔬菜種植戶推廣和應(yīng)用綠色防控措施。

從農(nóng)戶生產(chǎn)性特征出發(fā)可知影響農(nóng)戶防控行為的風(fēng)險(xiǎn)因素是多方面的，其中種植規(guī)模、勞動力強(qiáng)度等風(fēng)險(xiǎn)因素的影響機(jī)制相近，主要由于種植規(guī)模的擴(kuò)大需要有更多的勞動力投入。農(nóng)產(chǎn)品認(rèn)證、風(fēng)險(xiǎn)投保、生產(chǎn)專業(yè)組織資格這類風(fēng)險(xiǎn)因素的影響具有一定的平穩(wěn)性，因?yàn)榧尤肷a(chǎn)專業(yè)組織的農(nóng)戶風(fēng)險(xiǎn)規(guī)避動機(jī)更強(qiáng)，對產(chǎn)品質(zhì)量標(biāo)準(zhǔn)要求更高。綜合平穩(wěn)性特征和后驗(yàn)概率值比對，政府的生產(chǎn)性補(bǔ)貼對農(nóng)戶土壤處理和生物農(nóng)藥的施用影響較大，考慮到土壤消毒的成本和生物農(nóng)藥的價(jià)格，外部的政策環(huán)境對此類防控行為的約束力更強(qiáng)。生產(chǎn)者種植年限和風(fēng)險(xiǎn)投保意愿對農(nóng)戶綠色防控行為的影響較大，主要由于綠色防控的田間防治費(fèi)用低廉，種植經(jīng)驗(yàn)豐富、風(fēng)險(xiǎn)承受能力強(qiáng)的農(nóng)戶更重視綠色防控的長期效益而不會受到推廣成本的影響。另外，農(nóng)戶積極申請農(nóng)產(chǎn)品質(zhì)量認(rèn)證和加入專業(yè)生產(chǎn)合作組織也會增強(qiáng)其采用綠色防控措施的意愿?？傮w來說，農(nóng)業(yè)生產(chǎn)者的種植年限、是否享有政府補(bǔ)貼、風(fēng)險(xiǎn)投保、農(nóng)產(chǎn)品認(rèn)證、生產(chǎn)專業(yè)組織成員資格等因素都可能是致使防控行為的缺失的風(fēng)險(xiǎn)因素。

5.2 政策啟示

(1)病蟲害的防治應(yīng)當(dāng)貫徹“預(yù)防為主，綜合防治”的植保工作方針。合作社應(yīng)引導(dǎo)農(nóng)戶在藥劑防治過程中適時(shí)適量用藥，強(qiáng)化藥劑輪換施用的意識，減弱病蟲害的抗藥性。農(nóng)資店應(yīng)加強(qiáng)農(nóng)藥銷售市場的規(guī)范性，加大生物農(nóng)藥的存貨比例，嚴(yán)格按照國家要求禁止高毒高殘留的農(nóng)藥品種在銷售市場流通。同時(shí)，公共部門應(yīng)大力推廣高效防控技術(shù)的使用，尤其注重綜合防控中物理防治和生物防治的推廣和使用，逐步改變藥劑防治在綜合防治體系中的主導(dǎo)作用。政府在整個(gè)統(tǒng)防統(tǒng)治工作中，應(yīng)明確任務(wù)目標(biāo)、制定詳盡的工作方案、強(qiáng)化工作的督查指導(dǎo)，同時(shí)提供強(qiáng)有力的技術(shù)支持，確保統(tǒng)防統(tǒng)治工作順利有序的開展。

(2)推行和實(shí)施病蟲害綜合防控有關(guān)的獎懲制度，約束農(nóng)戶低效高風(fēng)險(xiǎn)的防控行為。擴(kuò)大農(nóng)藥補(bǔ)貼范圍和親環(huán)境技術(shù)的財(cái)政補(bǔ)貼力度，例如政府應(yīng)當(dāng)給予使用一定量環(huán)保型生物農(nóng)藥的農(nóng)戶相應(yīng)的補(bǔ)貼。加大懲罰力度，對于農(nóng)殘過高和土壤污染未能及時(shí)修復(fù)的農(nóng)戶應(yīng)當(dāng)上交一定的罰款。擴(kuò)大農(nóng)業(yè)保險(xiǎn)的險(xiǎn)種范圍，積極引導(dǎo)農(nóng)戶在災(zāi)害和價(jià)格上進(jìn)行投保，為防范農(nóng)業(yè)風(fēng)險(xiǎn)提供強(qiáng)有力的保障。建設(shè)規(guī)范化的信息平臺給農(nóng)戶適時(shí)傳遞正確的病蟲害防治信息。政府應(yīng)首先積極搭建農(nóng)村鄉(xiāng)鎮(zhèn)之間信息交流和服務(wù)平臺、電商平臺，形成當(dāng)?shù)厣舷侣?lián)動的農(nóng)戶防控行為網(wǎng)絡(luò)體系。

(3)保障蔬菜質(zhì)量，實(shí)現(xiàn)蔬菜從田間到餐桌的全過程管理，需要政府建立嚴(yán)格的農(nóng)產(chǎn)品質(zhì)量安全的監(jiān)管和控制體系[24]。定期進(jìn)行種植監(jiān)管和農(nóng)產(chǎn)品質(zhì)量檢測，積極引導(dǎo)農(nóng)戶申請農(nóng)產(chǎn)品更高級別的認(rèn)證。同時(shí)通過建立一套相對完整的農(nóng)產(chǎn)品質(zhì)量安全數(shù)據(jù)庫和信息服務(wù)平臺，逐步公開和透明農(nóng)產(chǎn)品質(zhì)量安全的相關(guān)信息。在農(nóng)產(chǎn)品質(zhì)量安全管理中促使農(nóng)業(yè)生產(chǎn)者綠色防控和更高效地防控病蟲害。發(fā)揮基層生產(chǎn)專業(yè)組織和龍頭企業(yè)對農(nóng)戶生產(chǎn)經(jīng)營活動的規(guī)整作用，圍繞優(yōu)良品種、高效防控、知識更新三大過程加強(qiáng)優(yōu)質(zhì)蔬菜基地建設(shè)。鼓勵龍頭企業(yè)面向基地、蔬菜種植戶推廣先進(jìn)防控技術(shù)，派遣技術(shù)人員指導(dǎo)田間工作，幫助農(nóng)戶調(diào)控棚室環(huán)境、改進(jìn)防控技術(shù)、建立合理種植制度[25]。充分發(fā)揮龍頭企業(yè)和示范基地標(biāo)準(zhǔn)化病蟲害防控工作的帶動作用和輻射效應(yīng)，積極引導(dǎo)農(nóng)戶樹立設(shè)施蔬菜質(zhì)量意識和品牌意識。