彭雯婧，羅 斌，曹喜春，李越凡，熊 錚，盧素珍

（江西省檢驗(yàn)檢測(cè)認(rèn)證總院食品檢驗(yàn)檢測(cè)研究院，江西南昌 330001)

食品安全與人們的身體健康息息相關(guān)，保障食品安全是社會(huì)穩(wěn)定的重要舉措。農(nóng)產(chǎn)品作為人們?nèi)粘Ｉ钪斜夭豢缮俚氖巢?，每年抽檢量在40%以上，不合格率較高。本研究以2020——2022 年食用農(nóng)產(chǎn)品監(jiān)督抽檢數(shù)據(jù)為例，引入二項(xiàng)分布檢驗(yàn)、卡方檢驗(yàn)等多種統(tǒng)計(jì)學(xué)方法對(duì)食品農(nóng)產(chǎn)品監(jiān)督抽檢數(shù)據(jù)進(jìn)行分析，量化食品品種風(fēng)險(xiǎn)值，結(jié)合場(chǎng)所、購(gòu)進(jìn)/生產(chǎn)月份、不合格項(xiàng)目分析，為食品品種進(jìn)行精準(zhǔn)畫(huà)像，以期為抽檢計(jì)劃執(zhí)行提供科學(xué)支持。

1 數(shù)據(jù)來(lái)源

通過(guò)國(guó)家市場(chǎng)監(jiān)督管理總局食品安全抽檢公布結(jié)果查詢系統(tǒng)收集2020——2022 年食用農(nóng)產(chǎn)品監(jiān)督抽檢數(shù)據(jù)（數(shù)據(jù)來(lái)源：https://spcjsac.gsxt.gov.cn/）作為樣本分析。

2 方法與步驟

本研究用Microsoft Excel、python、SPSS 等軟件對(duì)下載的數(shù)據(jù)進(jìn)行錄入處理和統(tǒng)計(jì)分析。

2.1 食品品種風(fēng)險(xiǎn)值的多維度分析與量化

食品品種的不合格率被廣泛用于食品安全風(fēng)險(xiǎn)評(píng)估，該指標(biāo)是一個(gè)動(dòng)態(tài)變量，在特定時(shí)間段內(nèi)，不合格率存在較大的隨機(jī)性，這給分析帶來(lái)了不確定性。采用二項(xiàng)分布檢驗(yàn)和基于食品品種不合格率的多維度分析方法，全面評(píng)估食品品種的安全風(fēng)險(xiǎn)。

2.1.1 不合格率風(fēng)險(xiǎn)分級(jí)

通常使用不合格率或不合格度等指標(biāo)來(lái)表征食品安全風(fēng)險(xiǎn)。以全年總體不合格率作為衡量的標(biāo)尺，將食品不合格率進(jìn)行分級(jí)賦分，從而體現(xiàn)該食品品種在全年中不同的風(fēng)險(xiǎn)程度。當(dāng)食品品種的不合格率所屬級(jí)別越高，就可以認(rèn)為該食品類(lèi)別的風(fēng)險(xiǎn)較高，并給予相應(yīng)的高分值。

2.1.2 二項(xiàng)分布檢驗(yàn)風(fēng)險(xiǎn)分級(jí)

二項(xiàng)分布檢驗(yàn)是一種適用于二分類(lèi)問(wèn)題的統(tǒng)計(jì)方法，用于檢驗(yàn)二項(xiàng)分布中比例是否等于某個(gè)特定值。本文將合格和不合格作為二分類(lèi)變量，然后使用二項(xiàng)分布檢驗(yàn)來(lái)計(jì)算每個(gè)食品品種合格率的置信區(qū)間。通過(guò)分析置信區(qū)間的寬度和位置，可以將食品品種劃分為低風(fēng)險(xiǎn)、中風(fēng)險(xiǎn)、高風(fēng)險(xiǎn)3 個(gè)等級(jí)。如果置信區(qū)間較寬，表示合格率的不確定性較大，食品品種的風(fēng)險(xiǎn)也較高。

2.1.3 不合格率的趨勢(shì)分級(jí)

將2020——2022 年各類(lèi)食品的不合格率通過(guò)最小二乘法進(jìn)行線性擬合，找到一條與實(shí)際觀測(cè)數(shù)據(jù)的殘差平方和最小的直線，進(jìn)而確定最佳擬合直線的參數(shù)，y＝kx＋b，y是因變量，x是自變量，k是斜率，b是截距。k在一定程度上表達(dá)了該食品品種不合格率在近幾年的整體發(fā)展趨勢(shì)，一般來(lái)說(shuō)k＞0 表示食品不合格率呈上升趨勢(shì)，有風(fēng)險(xiǎn)增高的趨勢(shì)，且k值越高，表示增長(zhǎng)的趨勢(shì)越高；k＜0 表示該類(lèi)食品不合格率整體呈下降趨勢(shì)。

2.1.4 食品品種風(fēng)險(xiǎn)量化

將3 個(gè)因素結(jié)果進(jìn)行分級(jí)量化，綜合評(píng)價(jià)各食品風(fēng)險(xiǎn)性。以比較范圍數(shù)據(jù)的總體不合格率作為不合格率風(fēng)險(xiǎn)分級(jí)的標(biāo)準(zhǔn)，將品種不合格率與總體不合格率的比值以及不合格率2%作為標(biāo)準(zhǔn)進(jìn)行劃分。二項(xiàng)分布檢驗(yàn)可以將食品品種風(fēng)險(xiǎn)分為高、中、低風(fēng)險(xiǎn)，增長(zhǎng)趨勢(shì)按k值進(jìn)行等級(jí)劃分。等級(jí)賦分具體分值情況見(jiàn)表1。

表1 風(fēng)險(xiǎn)指標(biāo)與分值對(duì)應(yīng)表

綜合上述風(fēng)險(xiǎn)因素，根據(jù)實(shí)際賦予不同權(quán)重對(duì)各品種風(fēng)險(xiǎn)值進(jìn)行量化計(jì)算，為抽檢計(jì)劃的制定和執(zhí)行提供指導(dǎo)。風(fēng)險(xiǎn)值計(jì)算公式為

式中：T為食品品種總風(fēng)險(xiǎn)值；vi為各項(xiàng)分值；αi為各分項(xiàng)權(quán)重（根據(jù)工作重點(diǎn)確定，如對(duì)該品種分析的關(guān)注度等），i=1，2，…。

2.2 食品品種抽樣月份、抽樣區(qū)域、檢驗(yàn)項(xiàng)目的風(fēng)險(xiǎn)點(diǎn)分析

卡方檢驗(yàn)是一種基于分布的假設(shè)檢驗(yàn)方法，用于綜合分析食品品種、生產(chǎn)月份、抽樣區(qū)域的不合格率的顯著差異性，并評(píng)估它們?cè)谑称菲贩N中的不同風(fēng)險(xiǎn)表達(dá)方式，可結(jié)合不合格檢驗(yàn)項(xiàng)目為食品品種進(jìn)行預(yù)警。

3 結(jié)果與分析

3.1 食品品種風(fēng)險(xiǎn)值

根據(jù)表1 對(duì)2020——2022 年食用農(nóng)產(chǎn)品的整體數(shù)據(jù)進(jìn)行賦分，假定所有食品抽檢的各分項(xiàng)權(quán)重α均為1，計(jì)算得到各食品品種的風(fēng)險(xiǎn)值見(jiàn)表2。表2 中食莢豌豆和蔥根據(jù)二項(xiàng)分布檢驗(yàn)被歸為低風(fēng)險(xiǎn)品種，不合格批次分別占總批次的1/8、2/39，不合格率得分較高，可能是因?yàn)槌闄z量和不合格批次量較少導(dǎo)致二項(xiàng)分布檢驗(yàn)與不合格率分級(jí)差別較大，后期監(jiān)管中可以適當(dāng)提高抽檢批次進(jìn)行風(fēng)險(xiǎn)驗(yàn)證。說(shuō)明風(fēng)險(xiǎn)分析應(yīng)避免單一指標(biāo)的片面解讀，需要考慮相關(guān)因素的綜合影響。13 個(gè)高風(fēng)險(xiǎn)品種中，結(jié)合不合格率、斜率得分、二項(xiàng)分布檢驗(yàn)分值分析，其他水產(chǎn)品（牛蛙）、海水蟹、豇豆和姜等幾類(lèi)產(chǎn)品在3 個(gè)維度的結(jié)果都屬于高風(fēng)險(xiǎn)分值品種，應(yīng)重點(diǎn)加強(qiáng)監(jiān)管；豆芽的不合格率得分和二項(xiàng)分布檢驗(yàn)分值都屬于高風(fēng)險(xiǎn)組別，斜率得分屬于降低組別，可在重點(diǎn)監(jiān)管范圍內(nèi)降低頻次。

表2 各食品品種風(fēng)險(xiǎn)值

3.2 食品品種時(shí)間、抽樣場(chǎng)所、檢驗(yàn)項(xiàng)目的風(fēng)險(xiǎn)點(diǎn)分析

選取高風(fēng)險(xiǎn)品種進(jìn)行風(fēng)險(xiǎn)點(diǎn)分析，包括豇豆、姜、海水蟹、其他水產(chǎn)品、香蕉、海水蝦、韭菜、淡水魚(yú)、辣椒、獼猴桃、生干籽類(lèi)、食莢豌豆及豆芽?？ǚ綑z驗(yàn)結(jié)果及不合格項(xiàng)目統(tǒng)計(jì)見(jiàn)表3。

表3 食品品種時(shí)間與抽樣場(chǎng)所卡方結(jié)果以及不合格項(xiàng)目匯總表

3.2.1 食品品種時(shí)間、抽樣場(chǎng)所風(fēng)險(xiǎn)點(diǎn)分析

食用農(nóng)產(chǎn)品是非預(yù)包裝食品，一般以購(gòu)進(jìn)/抽樣日期統(tǒng)計(jì)。食用農(nóng)產(chǎn)品抽樣場(chǎng)所主要有餐飲（如餐館和食堂等）和流通（如超市、菜市場(chǎng)等）場(chǎng)所。綜合分析，韭菜、姜、豇豆、香蕉、其他水產(chǎn)品（牛蛙）和淡水魚(yú)等食品品種絕大部分月份甚至全年的不合格率均具有顯著差異性，應(yīng)在全年加強(qiáng)抽檢和監(jiān)管頻率，并在抽樣場(chǎng)所對(duì)有差異性的品種提高抽檢數(shù)量與監(jiān)管力度；食莢豌豆為時(shí)間與場(chǎng)所沒(méi)有顯著性差異的高風(fēng)險(xiǎn)品種，可能是抽檢量和不合格批次量偏少，應(yīng)在全年對(duì)各場(chǎng)所適當(dāng)提高抽檢量與監(jiān)管力度進(jìn)行驗(yàn)證；對(duì)于豆芽、獼猴桃、海水蝦、海水蟹和生干籽類(lèi)等部分月份不合格率具有顯著差異的食品品種，針對(duì)重點(diǎn)時(shí)間和場(chǎng)所適當(dāng)提高抽檢量與監(jiān)管力度。

3.2.2 食品品種檢驗(yàn)項(xiàng)目的風(fēng)險(xiǎn)點(diǎn)分析

農(nóng)藥殘留不合格常檢出于蔬菜水果中。由于農(nóng)藥以及作物特性不同，部分作物與農(nóng)藥間有緊密關(guān)聯(lián)。4-氯苯氧乙酸鈉、6-芐基腺嘌呤對(duì)種子發(fā)芽、細(xì)胞分裂有促進(jìn)作用，且具有抑制胚根生長(zhǎng)的生理作用，不合格樣品只檢出于豆芽[1]。氯吡脲對(duì)瓜果類(lèi)植物有促進(jìn)花芽分化、保花保果、促進(jìn)果實(shí)膨大的作用，能改善獼猴桃的外觀，并延長(zhǎng)獼猴桃的貨架期，只檢出于獼猴桃[2]。腐霉利是一種微毒級(jí)的殺菌劑，被大量用于易被真菌和細(xì)菌感染的韭菜。豇豆與辣椒因成熟期較短且病蟲(chóng)害較多，不合格農(nóng)藥種類(lèi)也較多，需要重點(diǎn)關(guān)注[3]。禁用農(nóng)藥甲胺磷和限用農(nóng)藥（毒死蜱、克百威、氧樂(lè)果等）也在辣椒和豇豆等產(chǎn)品中被檢測(cè)出。以上情況主要是由農(nóng)戶在種植環(huán)節(jié)追求經(jīng)濟(jì)效益，過(guò)度依賴農(nóng)藥造成的。

獸藥殘留主要檢出于水產(chǎn)品中。恩諾沙星是一種用于治療動(dòng)物皮膚感染、呼吸道感染的抗生素，在其他水產(chǎn)品（牛蛙）、淡水魚(yú)、海水蝦中均有檢出，可能是養(yǎng)殖場(chǎng)或者商戶為了減少水產(chǎn)品在養(yǎng)殖過(guò)程或運(yùn)輸過(guò)程中的病死率，使用了含有大量恩諾沙星的飼料或暫養(yǎng)水，導(dǎo)致遷移性污染。此外，禁用的獸藥硝基呋喃類(lèi)藥物和孔雀石綠也被檢測(cè)出于其他水產(chǎn)品（牛蛙）、淡水魚(yú)、海水蝦。

重金屬鉛、鎘被檢出于海水蝦、海水蟹、辣椒、韭菜和姜。重金屬污染主要是由于農(nóng)產(chǎn)品在種植、養(yǎng)殖過(guò)程中受到環(huán)境污染造成的。如海水中的鎘超標(biāo)會(huì)造成甲殼類(lèi)生物重金屬的富集[4]，蔬菜受種植環(huán)境的污染造成重金屬殘留不合格。

質(zhì)量指標(biāo)（酸價(jià)）和生物毒素（黃曲霉毒素B1）不合格只檢出于生干籽類(lèi)（花生）?；ㄉ鳛橹参镉椭闹饕现?，在種植、儲(chǔ)存、運(yùn)輸或加工環(huán)節(jié)的不當(dāng)操作都有可能使花生發(fā)生霉變或者腐敗，造成花生內(nèi)部油脂的酸敗和黃曲霉毒素的產(chǎn)生與富集[5]。

綜上所述，13 個(gè)品種中農(nóng)獸藥殘留和重金屬污染問(wèn)題較為普遍，應(yīng)結(jié)合農(nóng)獸藥殘留與品種的關(guān)聯(lián)性以及農(nóng)業(yè)農(nóng)村部禁限用公告，對(duì)重點(diǎn)項(xiàng)目與品種進(jìn)行重點(diǎn)監(jiān)控。重金屬污染的主要原因是環(huán)境污染和食品品種本身特性，應(yīng)重點(diǎn)關(guān)注產(chǎn)生問(wèn)題品種的種植、養(yǎng)殖環(huán)節(jié)，對(duì)整體過(guò)程加強(qiáng)監(jiān)管。

4 結(jié)語(yǔ)

通過(guò)對(duì)2020——2022 年食用農(nóng)產(chǎn)品監(jiān)督抽檢的分析，將食用農(nóng)產(chǎn)品分為4 個(gè)風(fēng)險(xiǎn)等級(jí)，監(jiān)管部門(mén)可以根據(jù)不同風(fēng)險(xiǎn)等級(jí)的具體情況分配任務(wù)。對(duì)于其他水產(chǎn)品（牛蛙）、豇豆等多維度分析都屬于高風(fēng)險(xiǎn)的品種，應(yīng)加大抽檢與監(jiān)管力度，如開(kāi)展專(zhuān)項(xiàng)行動(dòng)等。對(duì)于部分如蔥、茭白等不同維度分析大部分屬于中高風(fēng)險(xiǎn)情況但不同分析方法風(fēng)險(xiǎn)情況有差異，應(yīng)適當(dāng)提高抽檢批次進(jìn)行風(fēng)險(xiǎn)驗(yàn)證，并結(jié)合各品種在不同月份與環(huán)節(jié)場(chǎng)所不合格率顯著性差異情況，合理執(zhí)行抽檢計(jì)劃。

根據(jù)不同的食品品種對(duì)農(nóng)獸藥殘留、重金屬污染等問(wèn)題的不同敏感程度以及種養(yǎng)殖環(huán)境影響，應(yīng)采取不同措施進(jìn)行干預(yù)和解決，結(jié)合農(nóng)獸藥的禁限用要求，針對(duì)不合格率較高的品種、時(shí)間、區(qū)域分布，合理安排檢查和抽檢工作，避免潛在的食品安全風(fēng)險(xiǎn)，提高監(jiān)管的效率和精確性。對(duì)重金屬污染嚴(yán)重的種植、養(yǎng)殖地進(jìn)行溯源分析。

由于食用農(nóng)產(chǎn)品的特殊性，需要各部門(mén)間加強(qiáng)合作與信息共享，農(nóng)業(yè)、環(huán)境等部門(mén)進(jìn)行聯(lián)動(dòng)監(jiān)管，加強(qiáng)對(duì)農(nóng)業(yè)生產(chǎn)者的管理和監(jiān)督，確保農(nóng)藥合理使用和安全使用，包括正確選擇以及使用農(nóng)藥、控制施藥頻次、遵守安全間隔期和對(duì)重金屬污染源的分析治理等，同時(shí)加大對(duì)農(nóng)產(chǎn)品的質(zhì)量檢測(cè)和監(jiān)督力度，建立完善的農(nóng)產(chǎn)品溯源制度，保障公眾的食品安全。

本文主要以2020——2022 年食用農(nóng)產(chǎn)品安全監(jiān)督抽檢數(shù)據(jù)為例進(jìn)行了分析與預(yù)警，剖析品種、項(xiàng)目、時(shí)間等方面的風(fēng)險(xiǎn)因素，為高風(fēng)險(xiǎn)種類(lèi)的風(fēng)險(xiǎn)畫(huà)像，期望通過(guò)信息化手段，為后期的監(jiān)管與抽檢工作提供方向，為政府監(jiān)管部門(mén)提供決策支持，有針對(duì)性地制定監(jiān)管措施，最大限度地提高監(jiān)管效率和精確性，確保食品安全。