地理標(biāo)志農(nóng)產(chǎn)品質(zhì)量安全風(fēng)險(xiǎn)評(píng)估及預(yù)警研究

張 彪，胡晨鈺，李 晶

（昆明理工大學(xué)管理與經(jīng)濟(jì)學(xué)院，云南昆明 650500）

0 引言

中國(guó)作為一個(gè)農(nóng)業(yè)大國(guó)，農(nóng)業(yè)是國(guó)民經(jīng)濟(jì)的根本，而“三農(nóng)問題”一直是我國(guó)國(guó)民經(jīng)濟(jì)發(fā)展所面臨的重要挑戰(zhàn)。2021 年是“十四五”開局之年，《2020 年農(nóng)產(chǎn)品質(zhì)量安全工作要點(diǎn)》中明確指出各企業(yè)要提升監(jiān)管能力，完善質(zhì)量監(jiān)管體系，防范質(zhì)量安全隱患，穩(wěn)步提升農(nóng)產(chǎn)品質(zhì)量安全水平，做好農(nóng)產(chǎn)品質(zhì)量安全工作對(duì)于實(shí)施鄉(xiāng)村振興戰(zhàn)略具有特殊且重要的意義［1］。與一般的農(nóng)產(chǎn)品相比，地理標(biāo)志農(nóng)產(chǎn)品在其原產(chǎn)地獨(dú)特的氣候和特殊的土壤條件等自然生態(tài)環(huán)境以及長(zhǎng)期沉淀的歷史人文因素影響下，形成了自己獨(dú)特的品質(zhì)和特征。然而，近些年來(lái)頻頻爆出地理標(biāo)志農(nóng)產(chǎn)品質(zhì)量安全事故，如贛南臍橙染色事件、五常大米香精事件、金華火腿敵敵畏事件等讓人瞠目結(jié)舌［2］。這些事故足以表明當(dāng)前我國(guó)在地理標(biāo)志農(nóng)產(chǎn)品質(zhì)量安全方面存在著明顯問題，嚴(yán)重破壞了地理標(biāo)志農(nóng)產(chǎn)品的品牌效益與市場(chǎng)信譽(yù)，阻礙了農(nóng)業(yè)的規(guī)模化與產(chǎn)業(yè)化發(fā)展。云南省作為地理標(biāo)志農(nóng)產(chǎn)品大省，截至2021 年，擁有地理標(biāo)志農(nóng)產(chǎn)品86 個(gè)［3］。因此，對(duì)地理標(biāo)志農(nóng)產(chǎn)品的質(zhì)量安全風(fēng)險(xiǎn)進(jìn)行評(píng)估與預(yù)警，對(duì)于降低農(nóng)產(chǎn)品質(zhì)量安全風(fēng)險(xiǎn)，保證品牌影響力具有重要意義。

很多學(xué)者在農(nóng)產(chǎn)品質(zhì)量安全相關(guān)理論與實(shí)踐方面進(jìn)行了大量探索，尤其是在地理標(biāo)志農(nóng)產(chǎn)品品牌建設(shè)、農(nóng)產(chǎn)品質(zhì)量安全風(fēng)險(xiǎn)因素研究、農(nóng)產(chǎn)品質(zhì)量安全風(fēng)險(xiǎn)評(píng)估等方面取得了豐碩成果。學(xué)者們大多基于不同視角，采用不同方法展開研究。在農(nóng)產(chǎn)品質(zhì)量安全風(fēng)險(xiǎn)評(píng)估方法方面，可分為單一風(fēng)險(xiǎn)因子評(píng)估和多指標(biāo)綜合評(píng)價(jià)［4］。在單一風(fēng)險(xiǎn)因子評(píng)估中，影響農(nóng)產(chǎn)品質(zhì)量安全的主要來(lái)源是化學(xué)性污染和生物性污染。齊艷麗等［5］對(duì)玉米及其秸稈中殘留的戊唑醇和吡唑醚菌酯進(jìn)行風(fēng)險(xiǎn)評(píng)估，并針對(duì)玉米生產(chǎn)中農(nóng)藥的安全使用提出建議；何祥祥等［6］通過(guò)評(píng)估草莓的革蘭氏陰性細(xì)菌、霉菌等微生物污染情況，為控制草莓微生物數(shù)量和開發(fā)殺菌抑菌技術(shù)提供了參考。但由于單一風(fēng)險(xiǎn)因子的評(píng)估因子較為單一且多用于化學(xué)、生物等自然科學(xué)領(lǐng)域，因此社會(huì)學(xué)科相關(guān)研究大多使用多指標(biāo)綜合評(píng)價(jià)法。多指標(biāo)綜合評(píng)價(jià)常用方法有：模糊評(píng)價(jià)法、Logistic 回歸模型、人工神經(jīng)網(wǎng)絡(luò)法等。祁南南［7］建立了靜態(tài)模糊綜合評(píng)價(jià)模型來(lái)評(píng)估水果質(zhì)量安全風(fēng)險(xiǎn)，雖能夠有效判斷水果的風(fēng)險(xiǎn)等級(jí)，但計(jì)算復(fù)雜，指標(biāo)權(quán)重的確立缺乏依據(jù)；張東玲等［8］提出基于Logistic 回歸方程模型的農(nóng)產(chǎn)品質(zhì)量風(fēng)險(xiǎn)評(píng)估方法，并對(duì)風(fēng)險(xiǎn)進(jìn)行了判別，但其可輸入的指標(biāo)數(shù)量過(guò)少，很難擬合出數(shù)據(jù)的真實(shí)分布；陳婷婷［9］基于BP 神經(jīng)網(wǎng)絡(luò)建立肉制品冷鏈物流質(zhì)量風(fēng)險(xiǎn)評(píng)估模型，并對(duì)風(fēng)險(xiǎn)等級(jí)進(jìn)行分類，但算法的學(xué)習(xí)速度較慢，且訓(xùn)練精度較低。因此，本文擬采用的評(píng)估模型為RF-DBN 模型，主要原因有：①RF 因具有強(qiáng)大的特征抽取和表達(dá)能力，被廣泛用于高維、海量數(shù)據(jù)的降維處理，去除冗余特征，防止過(guò)擬合現(xiàn)象出現(xiàn)；②DBN 模型前期采用無(wú)監(jiān)督學(xué)習(xí)，最后一層采用反向傳播神經(jīng)網(wǎng)絡(luò)自頂向下有監(jiān)督地微調(diào)整個(gè)模型，減少前向無(wú)監(jiān)督學(xué)習(xí)的整體誤差，使網(wǎng)絡(luò)結(jié)構(gòu)達(dá)到最優(yōu)。通過(guò)樣本數(shù)據(jù)的輸入，不斷學(xué)習(xí)并改進(jìn)網(wǎng)絡(luò)模型權(quán)值，能夠減少人為主觀性對(duì)結(jié)果造成的影響。

國(guó)內(nèi)外學(xué)者在地理標(biāo)志農(nóng)產(chǎn)品品牌建設(shè)及發(fā)展相關(guān)理論與實(shí)踐方面取得了豐碩成果，但在地理標(biāo)志農(nóng)產(chǎn)品質(zhì)量安全風(fēng)險(xiǎn)評(píng)估及預(yù)警方面的研究仍存在很大的提升空間，如構(gòu)建系統(tǒng)性的質(zhì)量安全風(fēng)險(xiǎn)指標(biāo)，而非采用單一性評(píng)價(jià)指標(biāo)等。鑒于此，本文以云南普洱茶為例，運(yùn)用全面質(zhì)量管理中“人、機(jī)、料、法、環(huán)”5 要素，建立地理標(biāo)志農(nóng)產(chǎn)品質(zhì)量安全風(fēng)險(xiǎn)指標(biāo)評(píng)估體系，再運(yùn)用隨機(jī)森林模型對(duì)指標(biāo)體系進(jìn)行降維，最終構(gòu)建云南省地理標(biāo)志農(nóng)產(chǎn)品質(zhì)量安全風(fēng)險(xiǎn)評(píng)估及預(yù)警的RF 和DBN 模型，并進(jìn)行實(shí)證研究。

1 地理標(biāo)志農(nóng)產(chǎn)品質(zhì)量安全風(fēng)險(xiǎn)評(píng)價(jià)指標(biāo)體系構(gòu)建

1.1 評(píng)價(jià)依據(jù)及分析

全面質(zhì)量管理是指一種以產(chǎn)品質(zhì)量為核心，全員參與以實(shí)現(xiàn)對(duì)產(chǎn)品進(jìn)行有效質(zhì)量控制和質(zhì)量改進(jìn)的管理體系，其中“人、機(jī)、料、法、環(huán)”是全面質(zhì)量管理理論中影響產(chǎn)品質(zhì)量的主要因素。因此，本文擬從供應(yīng)鏈角度出發(fā)，結(jié)合上述5 個(gè)因素對(duì)地理標(biāo)志農(nóng)產(chǎn)品在種植、加工、物流和銷售4個(gè)階段的質(zhì)量安全風(fēng)險(xiǎn)因素進(jìn)行分析，具體如下：

（1）種植階段質(zhì)量安全風(fēng)險(xiǎn)。與一般農(nóng)產(chǎn)品相比，地理標(biāo)志農(nóng)產(chǎn)品獨(dú)特的品質(zhì)和特征大多來(lái)源于原產(chǎn)地特有的氣候與特殊的土壤條件。在種植階段，“人”的因素可分為種植人員的種植經(jīng)驗(yàn)和種植技能培訓(xùn)水平，若種植人員缺乏種植知識(shí)及相關(guān)技能，則可能導(dǎo)致其在設(shè)備操作、病蟲害防治等對(duì)專業(yè)性要求較高的操作中出現(xiàn)問題，影響地理標(biāo)志農(nóng)產(chǎn)品質(zhì)量安全；“機(jī)”的因素可分為農(nóng)產(chǎn)品種植機(jī)械化水平和農(nóng)產(chǎn)品質(zhì)量檢測(cè)設(shè)備水平，種植機(jī)械化設(shè)備的缺乏將直接導(dǎo)致地理標(biāo)志農(nóng)產(chǎn)品生產(chǎn)效率低下，產(chǎn)量較低，質(zhì)量水平參差不齊；“料”的因素主要包括地理標(biāo)志農(nóng)產(chǎn)品種苗選擇及農(nóng)藥、化肥使用情況，不同品質(zhì)的種苗會(huì)生長(zhǎng)成不同質(zhì)量的農(nóng)產(chǎn)品，若濫用農(nóng)藥和化肥，將會(huì)在農(nóng)產(chǎn)品內(nèi)部殘留有害物質(zhì)，影響其質(zhì)量安全；“法”的因素主要包括病蟲害防治方法和地理標(biāo)志農(nóng)產(chǎn)品采摘方法，不同方法也會(huì)導(dǎo)致不同的質(zhì)量問題；“環(huán)”的因素包括地理標(biāo)志農(nóng)產(chǎn)品生長(zhǎng)環(huán)境和灌溉水質(zhì)，地理標(biāo)志農(nóng)產(chǎn)品的生長(zhǎng)對(duì)種植環(huán)境和水質(zhì)要求較高，這些都會(huì)對(duì)產(chǎn)品造成安全隱患。

（2）加工階段質(zhì)量安全風(fēng)險(xiǎn)。對(duì)地理標(biāo)志農(nóng)產(chǎn)品采摘之后進(jìn)行初加工，可增加其附加值，提升產(chǎn)品在市場(chǎng)上的競(jìng)爭(zhēng)力和價(jià)格水平。在加工階段，“人”的因素主要包括工作人員從業(yè)資質(zhì)和知識(shí)技能培訓(xùn)水平，地理標(biāo)志農(nóng)產(chǎn)品獨(dú)特的加工工藝也是其具有特殊品質(zhì)的原因之一，因此工作人員需要具有專業(yè)的加工知識(shí)和技能；“機(jī)”的因素主要包括加工機(jī)器自動(dòng)化程度和加工成品檢測(cè)設(shè)備情況，在加工過(guò)程中，機(jī)器的自動(dòng)化水平會(huì)直接影響產(chǎn)品的標(biāo)準(zhǔn)化，若使用劣質(zhì)的加工設(shè)備和質(zhì)量檢測(cè)設(shè)備，會(huì)導(dǎo)致加工工藝不完善，質(zhì)量不合格；“料”的因素主要包括原材料新鮮程度和食品添加劑情況，在加工中使用變質(zhì)的原材料或大量防腐劑和食品添加劑會(huì)影響農(nóng)產(chǎn)品質(zhì)量安全；“法”的因素主要包括地理標(biāo)志農(nóng)產(chǎn)品獨(dú)特的加工工藝和溫濕度控制技術(shù)；“環(huán)”的因素主要包括加工場(chǎng)所環(huán)境和倉(cāng)庫(kù)環(huán)境，地理標(biāo)志農(nóng)產(chǎn)品對(duì)溫度、濕度的敏感度較高，過(guò)濕或溫度過(guò)高極易導(dǎo)致其發(fā)霉變質(zhì)。此外，倉(cāng)庫(kù)不衛(wèi)生也會(huì)增加其質(zhì)量安全隱患。

（3）物流階段質(zhì)量安全風(fēng)險(xiǎn)。地理標(biāo)志農(nóng)產(chǎn)品在物流階段產(chǎn)生質(zhì)量安全風(fēng)險(xiǎn)主要來(lái)源于運(yùn)輸過(guò)程中的存儲(chǔ)問題和其包裝產(chǎn)生的安全隱患?！叭恕钡囊蛩刂饕獮楣ぷ魅藛T缺乏存儲(chǔ)經(jīng)驗(yàn)造成的質(zhì)量安全隱患；“機(jī)”的因素主要為質(zhì)量檢測(cè)設(shè)備情況，由于地理標(biāo)志農(nóng)產(chǎn)品具有保質(zhì)期短、易變質(zhì)的特點(diǎn)，若無(wú)法及時(shí)清理破損、變質(zhì)的農(nóng)產(chǎn)品，會(huì)威脅其質(zhì)量安全；“料”的因素主要包括地理標(biāo)志農(nóng)產(chǎn)品最終銷售成品質(zhì)量和農(nóng)產(chǎn)品外部包裝材質(zhì)質(zhì)量；“法”的因素主要包括包裝技術(shù)和溫濕度控制技術(shù)，因包裝不當(dāng)導(dǎo)致微生物含量超標(biāo)等會(huì)造成安全隱患；“環(huán)”的因素主要包括倉(cāng)庫(kù)存儲(chǔ)環(huán)境和運(yùn)輸環(huán)境，嚴(yán)格監(jiān)控物流階段的環(huán)境衛(wèi)生情況可大幅降低農(nóng)產(chǎn)品質(zhì)量安全風(fēng)險(xiǎn)［4］。

（4）銷售階段質(zhì)量安全風(fēng)險(xiǎn)。地理標(biāo)志農(nóng)產(chǎn)品不同于一般農(nóng)產(chǎn)品具有流動(dòng)銷售的特點(diǎn)，因此無(wú)法輕易辨別其質(zhì)量，但與一般農(nóng)產(chǎn)品類似，地理標(biāo)志農(nóng)產(chǎn)品具有易變質(zhì)、難儲(chǔ)存的特點(diǎn)。因此，在銷售階段，“人”的因素主要為工作人員存儲(chǔ)經(jīng)驗(yàn)；“機(jī)”的因素為質(zhì)量檢測(cè)設(shè)備情況，如果商品質(zhì)量檢驗(yàn)水平不高，難以剔除出已變質(zhì)、損耗的產(chǎn)品，則會(huì)造成質(zhì)量安全風(fēng)險(xiǎn)；“料”的因素主要為商品質(zhì)量水平，需要從商品的源頭控制其質(zhì)量安全風(fēng)險(xiǎn)；“法”的因素主要包括產(chǎn)品冷藏保鮮技術(shù)及商品標(biāo)識(shí)與可追溯管理水平；“環(huán)”的因素主要包括倉(cāng)庫(kù)存儲(chǔ)環(huán)境和銷售場(chǎng)所環(huán)境，需要從產(chǎn)品存儲(chǔ)到銷售環(huán)節(jié)都保證其質(zhì)量安全。

為更清晰地描述“人、機(jī)、料、法、環(huán)”5 要素與種植、加工、物流、銷售4 階段之間的關(guān)系，系統(tǒng)化地繪制質(zhì)量安全風(fēng)險(xiǎn)體系魚骨圖如圖1所示。

Fig.1 Quality and safety risk analysis fish bone map圖1 質(zhì)量安全風(fēng)險(xiǎn)分析魚骨圖

1.2 質(zhì)量安全風(fēng)險(xiǎn)評(píng)估指標(biāo)收集

目前，關(guān)于農(nóng)產(chǎn)品質(zhì)量安全風(fēng)險(xiǎn)評(píng)估指標(biāo)構(gòu)建的研究成果豐富，本文的主要研究對(duì)象是地理標(biāo)志農(nóng)產(chǎn)品，其本身就屬于農(nóng)產(chǎn)品范疇，因此在指標(biāo)構(gòu)建上可參考農(nóng)產(chǎn)品質(zhì)量安全相關(guān)文獻(xiàn)來(lái)確定。運(yùn)用文獻(xiàn)收集法對(duì)農(nóng)產(chǎn)品質(zhì)量安全評(píng)估體系指標(biāo)進(jìn)行初步的收集與選取，查找國(guó)內(nèi)外關(guān)于農(nóng)產(chǎn)品質(zhì)量安全評(píng)估體系的相關(guān)文獻(xiàn)，最后選取了10篇文獻(xiàn)資料，具體如表1所示。

1.3 質(zhì)量安全風(fēng)險(xiǎn)評(píng)估指標(biāo)體系構(gòu)建

首先，采用文獻(xiàn)收集法對(duì)質(zhì)量安全風(fēng)險(xiǎn)因素進(jìn)行收集與歸納，建立質(zhì)量安全風(fēng)險(xiǎn)評(píng)估體系；其次，運(yùn)用專家打分法，邀請(qǐng)7 名質(zhì)量風(fēng)險(xiǎn)管理領(lǐng)域的專家，對(duì)質(zhì)量安全風(fēng)險(xiǎn)評(píng)估體系進(jìn)行論證與篩選；最終，確定地理標(biāo)志農(nóng)產(chǎn)品質(zhì)量安全風(fēng)險(xiǎn)評(píng)估體系，其中包含4 個(gè)二級(jí)指標(biāo)，34 個(gè)三級(jí)指標(biāo)，具體如表2所示。

2 地理標(biāo)志農(nóng)產(chǎn)品質(zhì)量安全風(fēng)險(xiǎn)評(píng)估模型構(gòu)建

2.1 隨機(jī)森林模型

2.1.1 模型簡(jiǎn)介

隨機(jī)森林（Random Forests，RF）通常又稱為組合決策樹，其是通過(guò)集成學(xué)習(xí)思想將多個(gè)決策樹集成的一種算法，因其具有強(qiáng)大的特征抽取與表達(dá)能力，被廣泛用于高維、海量數(shù)據(jù)的降維處理。隨機(jī)森林的基本原理是從初始訓(xùn)練集M 中隨機(jī)且有放回地抽取m 個(gè)訓(xùn)練樣本，抽取出的每個(gè)訓(xùn)練樣本便是一棵決策樹（弱分類器），m 個(gè)樣本便組成了隨機(jī)森林。

Table 1 Information on indicators of agricultural product quality and safety evaluation system表1 農(nóng)產(chǎn)品質(zhì)量安全評(píng)價(jià)體系指標(biāo)資料

Table 2 Quality and safety risk assessment system表2 質(zhì)量安全風(fēng)險(xiǎn)評(píng)估體系

2.1.2 模型構(gòu)建

本文采用隨機(jī)森林模型對(duì)指標(biāo)體系進(jìn)行降維，對(duì)各個(gè)指標(biāo)特征的重要性進(jìn)行評(píng)估，計(jì)算隨機(jī)森林中每棵決策樹上的特征作了多少貢獻(xiàn)，依據(jù)特征重要性剔除相應(yīng)比例的特征，選擇最優(yōu)特征并對(duì)指標(biāo)重要性進(jìn)行排序，最終篩選出對(duì)結(jié)果影響較大的指標(biāo)實(shí)現(xiàn)數(shù)據(jù)降維。本文采用隨機(jī)森林模型對(duì)影響地理標(biāo)志農(nóng)產(chǎn)品質(zhì)量的安全風(fēng)險(xiǎn)因素重要性進(jìn)行排序，篩選掉對(duì)地理標(biāo)志農(nóng)產(chǎn)品質(zhì)量安全影響較小的指標(biāo)因素。RF 結(jié)構(gòu)如圖2所示。

Fig.2 RF structure圖2 RF結(jié)構(gòu)

2.2 深度置信網(wǎng)絡(luò)模型

2.2.1 模型簡(jiǎn)介

深度置信網(wǎng)絡(luò)（Deep Belief Networks，DBN）是2006 年由“深度學(xué)習(xí)之父”Hinton 提出的深度學(xué)習(xí)模型，其是一個(gè)將概率統(tǒng)計(jì)學(xué)與機(jī)器學(xué)習(xí)及神經(jīng)網(wǎng)絡(luò)相融合的概率生成模型。深度置信網(wǎng)絡(luò)是由若干個(gè)受限玻爾茲曼機(jī)（Re?stricted Boltzmann Machines，RBM）元件堆疊而成，通過(guò)不斷從底層特征中提取出抽象的高層特征后進(jìn)行分類，具有良好的學(xué)習(xí)能力。單層RBM 是包含一層可視層v 和一層隱藏層h 的無(wú)向圖模型，且層內(nèi)的單元之間不存在連接。RBM 結(jié)構(gòu)如圖3所示。

Fig.3 RBM structure圖3 RBM結(jié)構(gòu)

首先，在RBM 訓(xùn)練過(guò)程中，利用對(duì)比散度（Contrastive Divergence，CD）算法進(jìn)行參數(shù)更新，使得訓(xùn)練樣本的概率達(dá)到最大。每一層RBM 經(jīng)過(guò)多次迭代訓(xùn)練均得到初始化權(quán)值參數(shù)，形成無(wú)監(jiān)督的DBN 網(wǎng)絡(luò)。然后，采用BP 算法將最后一層RBM 的輸出值作為輸入值，自頂向下有監(jiān)督地微調(diào)整個(gè)模型，減少前向無(wú)監(jiān)督學(xué)習(xí)的整體誤差，達(dá)到最優(yōu)網(wǎng)絡(luò)結(jié)構(gòu)。最后，將得到的特征矩陣輸入softmax 分類器中，利用隨機(jī)梯度下降（Stochastic Gradient Descent，SGD）算法對(duì)損失函數(shù)進(jìn)行優(yōu)化，最終達(dá)到分類的目的。

2.2.2 模型構(gòu)建

本文將收集的數(shù)據(jù)集按照8∶2 的比例劃分為訓(xùn)練集和測(cè)試集，利用訓(xùn)練集訓(xùn)練RBM 模型，確定模型初始權(quán)值，將最后一層RBM 的輸出值當(dāng)作BP 神經(jīng)網(wǎng)絡(luò)的輸入層，對(duì)模型進(jìn)行有監(jiān)督的調(diào)優(yōu)。softmax 分類器在微調(diào)后對(duì)地理標(biāo)志農(nóng)產(chǎn)品的質(zhì)量安全風(fēng)險(xiǎn)進(jìn)行判別分類，并輸出風(fēng)險(xiǎn)等級(jí)。之后利用測(cè)試樣本集檢驗(yàn)DBN 模型分類效果，并將此模型用于地理標(biāo)志農(nóng)產(chǎn)品質(zhì)量安全風(fēng)險(xiǎn)分類與識(shí)別。DBN 模型結(jié)構(gòu)如圖4所示。

Fig.4 DBN structure圖4 DBN結(jié)構(gòu)

3 實(shí)證分析

3.1 數(shù)據(jù)收集

普洱茶作為云南省的地理標(biāo)志農(nóng)產(chǎn)品，在省內(nèi)被普遍種植，尤其位于瀾滄江中下游地區(qū)的普洱市茶區(qū)、西雙版納茶區(qū)、臨滄茶區(qū)三大茶區(qū)更是貢獻(xiàn)了90%以上的普洱茶原料。本文為驗(yàn)證模型，將選取300 個(gè)普洱茶企業(yè)，對(duì)其普洱茶質(zhì)量安全風(fēng)險(xiǎn)等級(jí)進(jìn)行評(píng)價(jià)。將其風(fēng)險(xiǎn)評(píng)價(jià)等級(jí)分為5 級(jí)，最低風(fēng)險(xiǎn)為1 級(jí)，最高風(fēng)險(xiǎn)為5 級(jí)。然后邀請(qǐng)7位質(zhì)量安全風(fēng)險(xiǎn)評(píng)估專家，根據(jù)企業(yè)相關(guān)情況采用問卷形式對(duì)指標(biāo)體系進(jìn)行打分，分?jǐn)?shù)范圍為0～10。分?jǐn)?shù)越接近10，代表該指標(biāo)的情況越優(yōu)異，所對(duì)應(yīng)的風(fēng)險(xiǎn)越低，如表3所示。最后根據(jù)得到的數(shù)據(jù)構(gòu)建風(fēng)險(xiǎn)矩陣表。

Table 3 Risk rating table表3 風(fēng)險(xiǎn)等級(jí)表

通過(guò)多輪匿名方式征詢質(zhì)量安全風(fēng)險(xiǎn)評(píng)估專家意見，不斷反饋、調(diào)整并綜合多數(shù)專家的主觀經(jīng)驗(yàn)和判斷，統(tǒng)計(jì)、處理、分析專家意見，對(duì)選取的300 家普洱茶企業(yè)的34 個(gè)普洱茶質(zhì)量安全風(fēng)險(xiǎn)指標(biāo)等級(jí)進(jìn)行客觀、合理估算，最終得到300 份普洱茶質(zhì)量安全風(fēng)險(xiǎn)指標(biāo)體系調(diào)查數(shù)據(jù)。將收集得到的專家打分?jǐn)?shù)據(jù)進(jìn)行歸一化處理，得到最終影響云南普洱茶質(zhì)量安全風(fēng)險(xiǎn)的指標(biāo)體系數(shù)據(jù)。

3.2 基于RF的指標(biāo)降維

將RF 模型導(dǎo)入影響云南普洱茶質(zhì)量安全風(fēng)險(xiǎn)的指標(biāo)體系數(shù)據(jù)，借助Python 軟件得到34個(gè)影響因素的重要性分布，如圖5所示。

Fig.5 Indicator importance distribution map圖5 指標(biāo)重要性分布圖

為更直觀地觀察數(shù)據(jù)，將各個(gè)指標(biāo)的重要性數(shù)據(jù)匯總成表格，如表4 所示。通過(guò)對(duì)表4 中的指標(biāo)重要性進(jìn)行分析，可得出：①重要性指數(shù)最高的指標(biāo)。最重要的因素是倉(cāng)庫(kù)存儲(chǔ)環(huán)境情況，重要性指數(shù)為0.11。茶葉需嚴(yán)格控制存儲(chǔ)的溫度、濕度以保證其口感、色澤等，因此物流期間的倉(cāng)庫(kù)存儲(chǔ)環(huán)境對(duì)于保證茶葉質(zhì)量至關(guān)重要；②重要性指數(shù)達(dá)到0.09～0.1 區(qū)間的指標(biāo)。在此區(qū)間的指標(biāo)包括加工場(chǎng)所環(huán)境、物流階段成品質(zhì)量檢測(cè)設(shè)備和銷售階段商品質(zhì)量檢測(cè)設(shè)備情況。普洱茶的加工工序較為繁瑣，對(duì)加工環(huán)境要求較為嚴(yán)苛，因此對(duì)加工場(chǎng)所環(huán)境的嚴(yán)格把控會(huì)提高普洱茶加工質(zhì)量。物流過(guò)程中的頻繁裝卸會(huì)給茶葉帶來(lái)一定損傷，采用良好的茶葉質(zhì)檢設(shè)備可規(guī)避質(zhì)量風(fēng)險(xiǎn)，而銷售階段作為供應(yīng)鏈的最后一環(huán)，采用良好的茶葉質(zhì)檢設(shè)備更顯得尤為重要；③重要性指數(shù)達(dá)到0.07～0.8 區(qū)間的指標(biāo)。在此區(qū)間的指標(biāo)包括銷售管理人員存儲(chǔ)經(jīng)驗(yàn)和農(nóng)產(chǎn)品加工工藝。若銷售管理人員對(duì)茶葉的存儲(chǔ)經(jīng)驗(yàn)不足，將會(huì)導(dǎo)致普洱茶變質(zhì)，不同的加工工藝也會(huì)大幅影響茶葉的風(fēng)味和口感，從而影響其質(zhì)量；④重要性指數(shù)在0.06～0.07區(qū)間的指標(biāo)。在此區(qū)間的指標(biāo)包括工作人員種植經(jīng)驗(yàn)、運(yùn)輸環(huán)境、生長(zhǎng)環(huán)境及加工人員知識(shí)技能培訓(xùn)情況。種植階段工作人員的種植經(jīng)驗(yàn)和茶葉生長(zhǎng)環(huán)境將直接影響到茶葉原材料質(zhì)量，良好的運(yùn)輸環(huán)境可減少因碰撞或溫濕度變化帶來(lái)的茶葉破碎或變質(zhì)。普洱茶的加工工藝特殊，因此需要對(duì)加工人員進(jìn)行技能培訓(xùn)，才能保證生產(chǎn)出符合質(zhì)量的茶葉；⑤重要性指數(shù)在0.05～0.06 之間的指標(biāo)。在此區(qū)間的指標(biāo)為農(nóng)產(chǎn)品采摘方式、倉(cāng)庫(kù)存儲(chǔ)環(huán)境及食品添加劑超標(biāo)情況。茶葉采摘方式分為機(jī)器和人工兩種，人工采摘的茶葉質(zhì)量較高，而加工階段的倉(cāng)庫(kù)存儲(chǔ)情況則確保從原葉、半成品到成品不會(huì)因存儲(chǔ)不當(dāng)而變質(zhì)。若食品添加劑過(guò)量，將會(huì)對(duì)人體造成嚴(yán)重影響，也必須重視。

Table 4 Indicator importance table表4 指標(biāo)重要性表

通過(guò)對(duì)表4 的進(jìn)一步分析，可得出重要性指數(shù)超過(guò)0.05 的指標(biāo)共有13 個(gè)，分別為：倉(cāng)庫(kù)存儲(chǔ)環(huán)境情況（C6）、加工場(chǎng)所環(huán)境情況（B9）、成品質(zhì)量檢測(cè)設(shè)備（C2）、商品質(zhì)量檢測(cè)設(shè)備（D2）、銷售管理人員存儲(chǔ)經(jīng)驗(yàn)（D1）、農(nóng)產(chǎn)品加工工藝（B7）、工作人員種植經(jīng)驗(yàn)（A1）、運(yùn)輸環(huán)境情況（C7）、生長(zhǎng)環(huán)境情況（A9）、人員知識(shí)技能培訓(xùn)（B2）、農(nóng)產(chǎn)品采摘方式（A8）、倉(cāng)庫(kù)存儲(chǔ)環(huán)境情況（B10）、食品添加劑超標(biāo)情況（B6）。由于其余21 個(gè)指標(biāo)的重要性指數(shù)不足0.01，與上述13 個(gè)指標(biāo)相比重要性較低，因此選取這13 個(gè)指標(biāo)對(duì)本文分類模型進(jìn)行訓(xùn)練，探尋各影響因素與普洱茶質(zhì)量安全風(fēng)險(xiǎn)等級(jí)之間的關(guān)系。降維后得到質(zhì)量安全風(fēng)險(xiǎn)評(píng)估體系，如表5所示。

3.3 基于DBN的風(fēng)險(xiǎn)評(píng)估分析

在以上運(yùn)用RF 降維之后選取的樣本中，將其中80%的樣本數(shù)據(jù)用于訓(xùn)練，剩余20%用于測(cè)試，并匯總出相應(yīng)等級(jí)的訓(xùn)練集和測(cè)試集數(shù)量，如表6所示。

Table 5 Quality and safety risk assessment system after dimension re?duction表5 降維后的質(zhì)量安全風(fēng)險(xiǎn)評(píng)估體系

Table 6 Sample distribution表6 樣本分布情況

根據(jù)3.2 中降維后的指標(biāo)體系，可確定模型的輸入層有13 個(gè)神經(jīng)元，輸出層有5 個(gè)神經(jīng)元，即等級(jí)1、等級(jí)2、等級(jí)3、等級(jí)4和等級(jí)5。本文根據(jù)得到的樣本數(shù)據(jù)對(duì)DBN 模型參數(shù)不斷進(jìn)行調(diào)整優(yōu)化，分別選用兩個(gè)RBM 模型進(jìn)行堆疊，其隱含層中的神經(jīng)元個(gè)數(shù)依次為50、35 個(gè)，并設(shè)RBM 前向?qū)W習(xí)率為0.01。根據(jù)訓(xùn)練樣本測(cè)試不同迭代次數(shù)對(duì)DBN 模型精度的影響，不同迭代次數(shù)的DBN 精度如圖6所示。

Fig.6 DBN accuracy for different iterations圖6 不同迭代次數(shù)的DBN精度

由圖6 可以得出，在迭代次數(shù)小于170 次時(shí)，模型精度隨著迭代次數(shù)的增加而增加。當(dāng)?shù)螖?shù)為170 次時(shí)，精度達(dá)到最高為96.67%；當(dāng)?shù)螖?shù)為180 次時(shí)，模型處于過(guò)擬合狀態(tài)，所以精度略有下降，且增加了迭代時(shí)間。因此，本文選擇迭代次數(shù)為170次。

接下來(lái)將全部數(shù)據(jù)輸入上述調(diào)整好的DBN 模型中，DBN 模型預(yù)測(cè)結(jié)果如圖7 所示。由圖7 可知，在普洱茶的質(zhì)量安全風(fēng)險(xiǎn)分類中，實(shí)際分類結(jié)果與模型分類預(yù)測(cè)結(jié)果相差不大。

Fig.7 DBN model prediction results圖7 DBN模型預(yù)測(cè)結(jié)果

分類預(yù)測(cè)匯總結(jié)果如表7 所示。由表7 可知，從等級(jí)1到等級(jí)3，分類精度高達(dá)100%，而分類錯(cuò)誤主要出現(xiàn)在等級(jí)4 和等級(jí)5 中，有兩組等級(jí)4 的樣本被錯(cuò)分為等級(jí)5。該結(jié)果說(shuō)明DBN 風(fēng)險(xiǎn)評(píng)估模型對(duì)普洱茶的質(zhì)量安全風(fēng)險(xiǎn)有較強(qiáng)的分類預(yù)測(cè)能力。

Table 7 Categorical predictions summary results表7 分類預(yù)測(cè)匯總結(jié)果

3.4 模型對(duì)比分析

反向傳播神經(jīng)網(wǎng)絡(luò)（BP）具有很強(qiáng)的非線性映射能力，是應(yīng)用最為廣泛的神經(jīng)網(wǎng)絡(luò)，而支持向量機(jī)（SVM）是一種被廣泛應(yīng)用于數(shù)據(jù)分類問題的監(jiān)督學(xué)習(xí)算法。因此，為進(jìn)一步驗(yàn)證DBN 模型的可靠性，本文選取傳統(tǒng)機(jī)器學(xué)習(xí)算法中的BP 神經(jīng)網(wǎng)絡(luò)和SVM 分類器分別對(duì)普洱茶質(zhì)量安全風(fēng)險(xiǎn)等級(jí)進(jìn)行分類測(cè)試，并對(duì)測(cè)試結(jié)果進(jìn)行對(duì)比。BP 神經(jīng)網(wǎng)絡(luò)選取3 層神經(jīng)網(wǎng)絡(luò)結(jié)構(gòu)，分別為輸入層、隱含層和輸出層。其中，輸入層和輸出層的神經(jīng)元數(shù)量與DBN 模型一致，分別為13 和5，且迭代次數(shù)為4 000 次，激活函數(shù)為Sig?moid 函數(shù)。在SVM 分類器中，設(shè)置內(nèi)核函數(shù)為RBF，懲罰系數(shù)為1，分類策略選擇“ovr（one versus rest）”，即一對(duì)多分類。3種模型分類結(jié)果比較如表8所示。

Table 8 Comparison of classification results of three models表8 3種模型分類結(jié)果比較

由表8可知，在測(cè)試樣本相同的情況下，可以得出：

（1）DBN 模型在普洱茶質(zhì)量安全風(fēng)險(xiǎn)分類的60 個(gè)測(cè)試樣本中，可正確分類出58 個(gè)樣本，評(píng)估精度較高，可達(dá)96.67%，明顯優(yōu)于BP 神經(jīng)網(wǎng)絡(luò)和SVM 分類器。該結(jié)果表明，DBN 模型可有效實(shí)現(xiàn)對(duì)地理標(biāo)志農(nóng)產(chǎn)品質(zhì)量安全的風(fēng)險(xiǎn)評(píng)估，且相較于BP 神經(jīng)網(wǎng)絡(luò)和SVM 分類器，該模型在質(zhì)量安全風(fēng)險(xiǎn)等級(jí)評(píng)估與分類中具有一定優(yōu)越性。

（2）BP 神經(jīng)網(wǎng)絡(luò)模型在60 個(gè)測(cè)試樣本中，可正確分類出54 個(gè)樣本，評(píng)估精度一般，僅為90%。BP 神經(jīng)網(wǎng)絡(luò)只有輸入層、隱含層和輸出層3 層網(wǎng)絡(luò)結(jié)構(gòu)，對(duì)于復(fù)雜函數(shù)關(guān)系的學(xué)習(xí)與訓(xùn)練能力明顯低于深度神經(jīng)網(wǎng)絡(luò)DBN，所以其對(duì)普洱茶質(zhì)量安全風(fēng)險(xiǎn)等級(jí)的評(píng)估結(jié)果一般。

（3）SVM 分類器在60 個(gè)測(cè)試樣本中，僅正確分類出46個(gè)樣本，評(píng)估精度較低，僅為85%。原因在于其最初設(shè)計(jì)為解決二值分類問題，因此在處理多分類問題時(shí)，只能間接構(gòu)造合適的多分類器，無(wú)法直接解決多分類問題，所以評(píng)估結(jié)果不理想。

4 結(jié)論與建議

近年來(lái)頻頻爆出地理標(biāo)志農(nóng)產(chǎn)品質(zhì)量安全事故，嚴(yán)重破壞了地理標(biāo)志農(nóng)產(chǎn)品的品牌效益與市場(chǎng)信譽(yù)，阻礙了農(nóng)業(yè)規(guī)?；?、產(chǎn)業(yè)化發(fā)展。本文充分利用深度學(xué)習(xí)在海量、高維數(shù)據(jù)中優(yōu)秀的特征學(xué)習(xí)與提取能力，將RF 和DBN 模型引入地理標(biāo)志農(nóng)產(chǎn)品質(zhì)量安全風(fēng)險(xiǎn)評(píng)估。通過(guò)實(shí)證分析，相對(duì)于傳統(tǒng)方式，該模型具有較高的評(píng)估準(zhǔn)確率。該研究結(jié)果為我國(guó)地理標(biāo)志農(nóng)產(chǎn)品質(zhì)量評(píng)估方法的應(yīng)用提供了參考。但本文研究中也存在一些不足：①由于數(shù)據(jù)獲取和標(biāo)量化的限制，本文只選取了部分影響地理標(biāo)志農(nóng)產(chǎn)品質(zhì)量安全的因素，實(shí)際上仍存在大量未被提及的因素有待繼續(xù)研究；②本文采用專家打分法獲取調(diào)查數(shù)據(jù)，具有一定程度上的主觀性。

后續(xù)可從以下兩方面進(jìn)行深入研究：①基于深度學(xué)習(xí)的地理標(biāo)志農(nóng)產(chǎn)品質(zhì)量安全風(fēng)險(xiǎn)評(píng)估方法優(yōu)化。RF 和DBN 只是眾多機(jī)器學(xué)習(xí)算法的一個(gè)，還可進(jìn)一步探索概率神經(jīng)網(wǎng)絡(luò)、無(wú)監(jiān)督學(xué)習(xí)等方法；②多種風(fēng)險(xiǎn)評(píng)估方法的有機(jī)融合。任何單一的風(fēng)險(xiǎn)評(píng)估方法都有自身的優(yōu)勢(shì)和局限，可通過(guò)組合評(píng)價(jià)方法充分發(fā)揮各方法的優(yōu)點(diǎn)，取長(zhǎng)補(bǔ)短，從而提高風(fēng)險(xiǎn)評(píng)估的準(zhǔn)確度。