蔣燕舒，劉玲君，李道霞，王 健，李仁偉

（四川省食品檢驗(yàn)研究院，四川 成都 611731）

0 引 言

食品安全關(guān)系人民群眾的身體健康和生命安全，關(guān)系中華民族的未來(lái).在《中共中央國(guó)務(wù)院關(guān)于深化改革加強(qiáng)食品安全工作的意見(jiàn)》［1］中指出，要“完善問(wèn)題導(dǎo)向的抽檢監(jiān)測(cè)機(jī)制”，其中，食品安全監(jiān)督抽檢是排查風(fēng)險(xiǎn)，倒逼食品生產(chǎn)經(jīng)營(yíng)者落實(shí)主體責(zé)任，撬動(dòng)稽查執(zhí)法，保障食品安全的一項(xiàng)重要監(jiān)管手段［2］.監(jiān)管部門在開(kāi)展食品安全監(jiān)督抽檢工作的過(guò)程中，收集了大量的抽檢數(shù)據(jù)，需要采取有效方法對(duì)數(shù)據(jù)加以分析利用，有針對(duì)性地制定監(jiān)管計(jì)劃［3］.現(xiàn)有的研究情況顯示：食品安全監(jiān)督抽檢數(shù)據(jù)分析較多局限于描述性統(tǒng)計(jì)分析［4-5］，或是針對(duì)于單一食品品種的空間分析［6-7］，缺乏對(duì)某一地區(qū)整體食品安全監(jiān)督抽檢數(shù)據(jù)進(jìn)行統(tǒng)計(jì)學(xué)空間分布和時(shí)空聚集分布的分析研究.隨著大數(shù)據(jù)時(shí)代的到來(lái)，分析食品抽檢數(shù)據(jù)，反映其空間規(guī)律、地區(qū)聯(lián)系和時(shí)空特征，將為食品安全風(fēng)險(xiǎn)防控提供針對(duì)性的決策［8］.本研究以 2017—2020年四川省省級(jí)食品監(jiān)管部門公布的監(jiān)督抽檢數(shù)據(jù)為例進(jìn)行分析研究，擬采用非參數(shù)檢驗(yàn)、單樣本t檢驗(yàn)、空間自相關(guān)性和時(shí)空掃描統(tǒng)計(jì)量等多種統(tǒng)計(jì)學(xué)方法粗略評(píng)估當(dāng)?shù)氐氖称钒踩珷顩r，從食品品種、不合格指標(biāo)、年度和地域4個(gè)維度分別探討當(dāng)?shù)刂饕媾R的食品安全風(fēng)險(xiǎn)，為靶向制訂食品安全監(jiān)管方案，提供科學(xué)依據(jù).

1 材料與方法

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

食品安全監(jiān)督抽檢結(jié)果來(lái)源于四川省人民政府持續(xù)在網(wǎng)站發(fā)布的數(shù)據(jù).下載2017—2020年發(fā)布的食品安全監(jiān)督抽檢信息，作為研究對(duì)象，分析當(dāng)前該省的食品安全狀況及存在的主要問(wèn)題.

1.2 方法

本研究應(yīng)用Microsoft Excel 2016和SPSS 26.0對(duì)下載的數(shù)據(jù)進(jìn)行錄入處理和統(tǒng)計(jì)分析，使用Geoda和SatScan軟件進(jìn)行空間統(tǒng)計(jì)分析和時(shí)空掃描統(tǒng)計(jì)量分析，尋找熱點(diǎn)聚集區(qū)域及相對(duì)風(fēng)險(xiǎn)度，并采用Q-gis軟件對(duì)時(shí)空分析結(jié)果進(jìn)行可視化呈現(xiàn).P＜0.05為差異有統(tǒng)計(jì)學(xué)意義.

1.2.1 線性擬合

假設(shè)某類食品4年間的合格率存在函數(shù)關(guān)系，作線性擬合曲線y=ax+b，式中：y為合格率；a為斜率，表示某類食品4年間合格率的趨勢(shì)變化，a＞0表明合格率整體呈上升趨勢(shì)，a＜0表明合格率整體呈下降趨勢(shì)；x為年度；b為截距；R2為線性相關(guān)系數(shù).

1.2.2 非參數(shù)檢驗(yàn)

采用Kruskal-Wallis檢驗(yàn)進(jìn)行多個(gè)獨(dú)立樣品之間相同分布研究［9-11］.零假設(shè)分別為同一類別食品不同年份之間的合格率無(wú)顯著性差異，不同類別食品的合格率無(wú)顯著性差異.

1.2.3 柯?tīng)柲曷宸?斯米諾夫（Kolmogorov-Smirnov，K-S）檢驗(yàn)

K-S檢驗(yàn)是研究樣本觀察值的分布和設(shè)定的理論分布是否相一致.零假設(shè)是各年度不合格樣品中各類不合格指標(biāo)項(xiàng)次服從正態(tài)分布.

1.2.4 單樣本t檢驗(yàn)

單樣本t檢驗(yàn)是用t分布理論來(lái)推論差異發(fā)生的概率［12-13］.零假設(shè)是各類不合格指標(biāo)項(xiàng)次平均值與某類不合格指標(biāo)項(xiàng)次數(shù)的差異無(wú)顯著性.

1.2.5 空間自相關(guān)分析

空間自相關(guān)是指同一變量在相應(yīng)空間位置上的相關(guān)性與相互依賴性［14-15］，選用莫蘭指數(shù)（Moran’s I）作為空間自相關(guān)指數(shù)［16］，利用 Moran’s I對(duì)食品抽檢合格率的空間聚集性進(jìn)行識(shí)別，應(yīng)用檢驗(yàn)統(tǒng)計(jì)量Z來(lái)檢驗(yàn) Moran’s I.當(dāng)||Z ＞1.96，則 P＜0.05，全局 Moran’s I＞0，表示空間正相關(guān)，其值越大，空間相關(guān)性越明顯；全局 Moran’s I＜0，表示空間負(fù)相關(guān)，其值越小，空間差異越大；全局Moran’s I=0，表示無(wú)空間相關(guān)性.

采用 Anselin［17］提出的局部空間自相關(guān)（local indicators of spatial association，LISA）分析方法尋找可能被掩蓋的局部空間自相關(guān)的位置，進(jìn)一步探索每個(gè)空間單元上食品抽檢合格率與其鄰近區(qū)域食品抽檢合格率的空間聚集性.一個(gè)空間單元i的局部Moran統(tǒng)計(jì)量可定義為

1.2.6 時(shí)空掃描分析

時(shí)空掃描分析同時(shí)考慮了時(shí)間和空間2個(gè)因素［18-19］，以四川省21個(gè)被抽檢市（州）為對(duì)應(yīng)研究區(qū)域，建立一個(gè)可活動(dòng)的掃描窗口，利用掃描窗口內(nèi)和外的實(shí)際抽檢不合格樣品數(shù)和理論不合格樣品數(shù)構(gòu)造檢驗(yàn)統(tǒng)計(jì)量對(duì)數(shù)似然比（log likelihood radio，LLR）［20］，用 LLR來(lái)評(píng)價(jià)掃描窗口內(nèi)抽檢不合格數(shù)的異常程度，然后采用蒙特卡羅模擬法計(jì)算P值對(duì)該窗口的統(tǒng)計(jì)學(xué)意義進(jìn)行評(píng)價(jià)［21］.本研究在分析過(guò)程中，掃描最大半徑為50%，時(shí)間間隔為年，假設(shè)該省食品安全監(jiān)督抽檢結(jié)果在時(shí)間和空間均呈柏松分布［22］，不合格批次數(shù)作為事件數(shù)，抽檢總批次作為總案例數(shù).當(dāng)LLR值具有統(tǒng)計(jì)學(xué)意義時(shí)，用相對(duì)危險(xiǎn)度（risk ratio，RR）來(lái)表示掃描窗口內(nèi)外的相對(duì)危險(xiǎn)度，首先掃描探測(cè)到LLR值為最大者，稱之為主要聚集區(qū)，其他掃描探測(cè)出的具有統(tǒng)計(jì)學(xué)意義的地區(qū)為次要聚集區(qū).

2 結(jié) 果

2.1 整體情況

2017—2020年四川省省級(jí)監(jiān)管部門共完成77 054批次食品監(jiān)督抽檢，其中合格樣品74 918批次，不合格樣品2 136批次，樣品總體合格率為97.23%.近4年抽檢合格率、合格和不合格批次變化趨勢(shì)如圖1所示.2017—2020年食品安全監(jiān)督抽檢完成數(shù)量整體呈上升趨勢(shì)，2020年完成的監(jiān)督抽檢數(shù)量比2017年增長(zhǎng)45.64%.4年間的樣品合格率呈現(xiàn)先降低后增高的趨勢(shì)，合格率最高與最低相差僅為0.41%，4年間合格率標(biāo)準(zhǔn)偏差為0.18%，各年度合格率均保持在97.00%以上.

圖1 2017—2020年抽檢合格率及批次變化趨勢(shì)

2.2 食品類別

2017—2020年監(jiān)督抽檢共覆蓋34類食品，其抽檢合格率相關(guān)結(jié)果如圖2所示.其中，可可及焙烤咖啡產(chǎn)品進(jìn)行了2年抽檢，特殊醫(yī)學(xué)用途配方食品、嬰幼兒配方食品和食鹽進(jìn)行了3年抽檢，其余食品類別均進(jìn)行了4年抽檢.可可及焙烤咖啡產(chǎn)品、特殊醫(yī)學(xué)用途配方食品和嬰幼兒配方食品平均合格率為100.00%，有24類食品平均合格率＞97.23%，酒類、糕點(diǎn)、食用農(nóng)產(chǎn)品、餐飲食品、炒貨食品及堅(jiān)果制品、冷凍飲品、蔬菜制品、飲料、特殊膳食食品和蜂產(chǎn)品10個(gè)類別的食品平均合格率＜97.23%，蜂產(chǎn)品平均合格率最低，為90.77%.

圖2 2017—2020年不同食品類別的抽檢合格率及變化趨勢(shì)

2017—2020年34類食品的合格率線性擬合結(jié)果表明：有3類食品4年間合格率持平，為100.00%；有13類食品4年來(lái)合格率整體呈上升趨勢(shì)，其中肉質(zhì)品呈逐年上升趨勢(shì)（R2=0.942）；有18類食品4年來(lái)合格率整體呈下降趨勢(shì)，其中蔬菜制品呈現(xiàn)逐年下降趨勢(shì)（R2=0.923）.采用非參數(shù)檢驗(yàn)法對(duì)2017—2020年34類食品的合格率差異進(jìn)行檢驗(yàn)，表明樣本中不同年份之間的合格率不存在顯著性差異（χ2=3.00，P＞0.05）；對(duì)不同食品類別之間的合格率差異進(jìn)行檢驗(yàn)，樣本中不同食品類別之間不存在顯著差異（χ2=3.65，P＞0.05）.

2.3 不合格項(xiàng)目

2017—2020年共計(jì)檢出2 339項(xiàng)次不合格指標(biāo)，將不合格指標(biāo)分為食品添加劑，微生物，農(nóng)藥獸藥，質(zhì)量指標(biāo)，重金屬、無(wú)機(jī)和有機(jī)污染物，非食用物質(zhì)和生物毒素污染7類，這些指標(biāo)項(xiàng)不合格次數(shù)列于表1.由K-S檢驗(yàn)可知各年度7類不合格指標(biāo)項(xiàng)次的分布為正態(tài)分布（P=0.20＞0.05），故可用單樣本t檢驗(yàn)來(lái)分析不合格指標(biāo)平均值與各類不合格指標(biāo)之間的差異.以各年度平均值為對(duì)照來(lái)研究不合格指標(biāo)的顯著性，食品添加劑不合格指標(biāo)項(xiàng)次連續(xù)4年顯著高于年度平均值（t=-3.75、-4.47、-3.32和-3.78，P＜0.05），生物毒素污染指標(biāo)項(xiàng)次在 2017、2018和2020年顯著低于年度平均值（t=2.83、2.65和2.64，P＜0.05）；微生物不合格指標(biāo)項(xiàng)次在 2017和2018年顯著高于年度平均值（t=-3.01和-2.59，P＜0.05）；農(nóng)藥獸藥不合格指標(biāo)項(xiàng)次在2019年顯著高于該年度平均值（t=-3.56，P＜0.05）；重金屬、無(wú)機(jī)和有機(jī)污染物（t=2.46，P＜0.05）和非食用物質(zhì)（t=2.76，P＜0.05）不合格指標(biāo)項(xiàng)次分別在2019和2020年顯著低于年度平均水平.

表1 2017—2020年7類指標(biāo)不合格次數(shù)統(tǒng)計(jì)

2.4 空間聚集性

2.4.1 全局空間自相關(guān)

2017—2020年抽檢合格率全局空間自相關(guān)分析結(jié)果列于表2.2017—2020年四川省食品抽檢合格率Moran’s I均不為0，經(jīng)統(tǒng)計(jì)學(xué)檢驗(yàn)，4年間總體食品抽檢合格率呈正向空間相關(guān)分布特征（Moran’s I＞0，Z＞1.96，P＜0.05），即高合格率市（州）在空間上彼此聚集，且低合格率市（州）在空間上也彼此聚集.其余各年度呈隨機(jī)分布.

表2 2017—2020年抽檢合格率全局空間自相關(guān)結(jié)果

2.4.2 局部空間自相關(guān)

2017—2020年食品抽檢合格率的局部空間自相關(guān)分析結(jié)果如圖3所示.在P=0.05檢驗(yàn)水平上，2017年，涼山州呈高-高聚集分布（Ii=0.184，P＜0.05），表明該州的合格率高，同時(shí)周邊的城市合格率也高；廣安市呈低-低聚集分布（Ii=1.465，P＜0.05），表明該市的合格率低，同時(shí)周邊的城市合格率也低.2018年，自貢市和雅安市呈高-低聚集分布（Ii=-0.304和-0.383，P＜0.05），其合格率高，但周邊的城市合格率低；德陽(yáng)市和甘孜州呈低-高分聚集布（Ii=-0.026和-1.378，P=0.05），該2市（州）的合格率低，但周邊的城市合格率高.2019年，瀘州市和遂寧市呈高-低聚集分布（Ii=-0.220和-1.589，P＜0.05），雅安市呈低-高聚集分布（Ii=-0.262，P＜0.05）.2020年，眉山市和資陽(yáng)市呈低-低聚集分布（Ii=0.783和0.935，P＜0.05）.2017—2020年總體來(lái)看，自貢市、內(nèi)江市、眉山市和資陽(yáng)市呈低-低聚集分布（Ii=0.550、1.617、1.164和 0.474，P＜0.05），表明該 4座城市的合格率低，同時(shí)周邊的城市合格率也低.

圖3 近4年四川省各地區(qū)抽檢合格率局部空間自相關(guān)結(jié)果

2.4.3 時(shí)空掃描統(tǒng)計(jì)

對(duì)2017—2020年食品安全監(jiān)督抽檢數(shù)據(jù)進(jìn)行時(shí)空掃描分析，相應(yīng)掃描結(jié)果如圖4和表3所示.2017—2020年四川省食品安全抽檢情況具有統(tǒng)計(jì)學(xué)意義的時(shí)空聚類區(qū)域有3個(gè).主要聚集區(qū)位于四川中偏南部，涉及樂(lè)山市和眉山市，聚集時(shí)間為2018年，該聚集區(qū)內(nèi)食品抽檢不合格樣品風(fēng)險(xiǎn)是其他區(qū)的2.20倍；次要聚集區(qū)Ⅰ位于四川中偏東部，涉及資陽(yáng)市和內(nèi)江市，聚集時(shí)間為2019—2020年，該聚集區(qū)內(nèi)食品抽檢不合格樣品風(fēng)險(xiǎn)是其他區(qū)的1.60倍；次要聚集區(qū)Ⅱ位于四川東部，涉及廣安市、南充市、遂寧市和達(dá)州市，聚集時(shí)間為2017年，該聚集區(qū)內(nèi)食品抽檢不合格樣品風(fēng)險(xiǎn)是其他區(qū)的1.58倍.

表3 2017—2020年抽檢不合格率時(shí)空掃描分析結(jié)果

圖4 2017—2020年抽檢不合格率時(shí)空掃描分析風(fēng)險(xiǎn)聚集區(qū)

3 結(jié)論與討論

3.1 總體形勢(shì)分析

2017—2020年四川省省級(jí)的食品監(jiān)督抽檢數(shù)量顯著上升，表明監(jiān)管部門在制定抽檢計(jì)劃時(shí)能夠按照《“十三五”國(guó)家食品安全規(guī)劃》［23］要求，不斷加大食品檢驗(yàn)量.2017—2020年間整體合格率顯示，四川省的食品安全形勢(shì)總體保持穩(wěn)定態(tài)勢(shì)，可見(jiàn)在堅(jiān)持以問(wèn)題為導(dǎo)向，以排查風(fēng)險(xiǎn)為目的開(kāi)展食品監(jiān)督抽檢的要求下，監(jiān)管部門能夠保障該省的整體食品安全合格率高位運(yùn)行，食品安全監(jiān)管工作效果顯著.食品安全抽樣檢驗(yàn)?zāi)甓扔?jì)劃要重點(diǎn)選擇風(fēng)險(xiǎn)程度高以及污染水平呈上升趨勢(shì)的食品，由此可知上一年度食品安全抽檢結(jié)果中發(fā)現(xiàn)問(wèn)題較多、風(fēng)險(xiǎn)較高的食品品種、檢測(cè)指標(biāo)和抽樣地區(qū)等將會(huì)納入下一年度的食品抽檢計(jì)劃當(dāng)中，由于計(jì)劃的逐年調(diào)整和動(dòng)態(tài)變化，由該4年的抽檢結(jié)果推測(cè)之后該省的年度食品抽檢合格率也將在97.00%以上小幅波動(dòng).

3.2 食品類別分析

區(qū)別于劉歡［24］僅用合格率判定當(dāng)?shù)厥称钒踩珷顩r的研究，本文采用非參數(shù)檢驗(yàn)得出，2017—2020年，食品安全監(jiān)督抽檢年度間的合格率與各食品類別間的合格率均沒(méi)有顯著性差異，且平均合格率均保持在97.00%以上，進(jìn)一步研究表明該省4年間的食品抽檢合格率一直保持較高水平，說(shuō)明在該省流通銷售的這34類食品，絕大部分生產(chǎn)加工者對(duì)食品品質(zhì)的把控基本都保持在較高水平.從歷年抽檢合格率的趨勢(shì)來(lái)看，可可及焙烤咖啡產(chǎn)品、特殊醫(yī)學(xué)用途配方食品和嬰幼兒配方食品4年間合格率總體保持在100.00%，其中，可可及焙烤咖啡產(chǎn)品、特殊醫(yī)學(xué)用途配方食品抽檢批次較少，不具備代表性，而嬰幼兒配方乳粉由于國(guó)家采取全國(guó)實(shí)施月月抽檢，嚴(yán)格嬰幼兒配方食品生產(chǎn)許可審查發(fā)證等最強(qiáng)監(jiān)管措施，確保嬰幼兒配方乳粉的食品安全風(fēng)險(xiǎn)最低.食鹽、食用油、油脂及其制品、乳制品、豆制品、食糖、蛋制品、肉制品、保健食品、酒類、食用農(nóng)產(chǎn)品、餐飲食品、炒貨食品及堅(jiān)果制品和冷凍飲品共13類食品抽檢合格率呈整體上升趨勢(shì)，表明該省監(jiān)管部門對(duì)油、肉、蛋、乳和菜等居民大宗消費(fèi)食品的品質(zhì)把控更為嚴(yán)格.其余18類食品抽檢合格率整體呈下降趨勢(shì)，合格率較低且合格率下降趨勢(shì)明顯的食品類別有蜂產(chǎn)品、蔬菜制品、特殊膳食食品、飲料和糕點(diǎn)等，其中特殊膳食食品抽檢批次較少，不具有代表性.蜂產(chǎn)品平均合格率最低且下降趨勢(shì)最為明顯，以此展開(kāi)分析，蜂產(chǎn)品不合格的原因主要有：（1）蜂蜜中檢出禁用獸藥諾氟沙星、氧氟沙星和氯霉素等；（2）蜂王漿中質(zhì)量指標(biāo)10-羥基-2-癸烯酸不達(dá)標(biāo).蜂產(chǎn)品是一種動(dòng)物源性食品，雖然中華人民共和國(guó)農(nóng)業(yè)部公告第2292號(hào)已明確規(guī)定“自2016年12月31日起，停止經(jīng)營(yíng)、使用用于食品動(dòng)物的洛美沙星、培氟沙星、氧氟沙星和諾氟沙星4種原料藥的各種鹽、酯及其各種制劑”［25］，但在 2019—2020年的抽檢蜂產(chǎn)品中仍有禁用獸藥的檢出情況，主要由于目前還無(wú)有效的低毒、低殘留可替代抗生素的藥物來(lái)治療蜂病，且蜂農(nóng)使用純化工藝不達(dá)標(biāo)的氟苯尼考藥物也會(huì)造成蜂產(chǎn)品中氯霉素檢出［26］.10-羥基-2-癸烯酸是又稱王漿酸，是蜂王漿特有的成分，是衡量蜂王漿品質(zhì)和辨別其真?zhèn)蔚闹匾笜?biāo).該省蜂王漿中不合格的情況主要為10-羥基-2-癸烯酸含量不滿足《蜂王漿：GB 9697—2008》［27］中產(chǎn)品等級(jí)要求，10-羥基-2-癸烯酸易受空氣、光線、水分和酸堿的影響［28］，由于生產(chǎn)企業(yè)對(duì)原料品質(zhì)把控和投料控制不嚴(yán)、生產(chǎn)工藝設(shè)計(jì)不合格，生產(chǎn)經(jīng)營(yíng)企業(yè)貯存條件不達(dá)標(biāo)等易導(dǎo)致其含量不達(dá)標(biāo).

3.3 不合格項(xiàng)目分析

區(qū)別于郭添榮等［29］僅采用不合格項(xiàng)目占比情況進(jìn)行描述性分析，本文采用以年度平均值和單個(gè)不合格項(xiàng)目次數(shù)進(jìn)行單樣本t檢驗(yàn)的方法，精準(zhǔn)地確定超范圍、超限量使用食品添加劑連續(xù)4年成為該省食品安全問(wèn)題中最為突出的問(wèn)題，其中較為多發(fā)的問(wèn)題有蔬菜制品干筍中SO2殘留量超標(biāo)，餐飲食品油條中Al的殘留量超標(biāo)，酒類散裝白酒中檢出甜蜜素（以C6H13NNaO3S計(jì)）等.根據(jù)我國(guó)食品添加劑的使用規(guī)定，SO2可作為漂白劑、防腐劑和抗氧化劑加入干制蔬菜中，KAl（SO4）2和NH4Al（SO4）2可作為膨松劑、穩(wěn)定劑加入油炸面制品中.由于食品添加劑具有防腐、護(hù)色、漂白、增香、乳化和增稠等多種功能，食品添加劑的使用一直是食品生產(chǎn)加工中必不可少的重要一環(huán)，且食品添加劑容易獲取，食品生產(chǎn)加工者的違法成本低且需要檢測(cè)才能確定其含量，是導(dǎo)致該指標(biāo)一直突顯的重要原因.此外，食品添加劑指標(biāo)不合格分為超出《食品安全國(guó)家標(biāo)準(zhǔn) 食品添加劑使用標(biāo)準(zhǔn)：GB 2760—2014》［30］規(guī)定的使用范圍和超過(guò)使用限量2種情形，該指標(biāo)不合格的原因一般還包括生產(chǎn)加工企業(yè)或個(gè)人：（1）專業(yè)技術(shù)缺乏，對(duì)標(biāo)準(zhǔn)規(guī)定不熟悉、理解不到位［31］；（2）為了使產(chǎn)品達(dá)到更好的效果而追求更多的經(jīng)濟(jì)利益，人為故意過(guò)量添加；（3）原材料進(jìn)貨檢驗(yàn)把關(guān)不嚴(yán)，原料中的食品添加劑超標(biāo)從而導(dǎo)致終產(chǎn)品中食品添加超標(biāo)；（4）所使用的食品添加劑本身不合格.

次要突出的問(wèn)題為微生物指標(biāo)不合格，連續(xù)在2017和2018年顯著較高.微生物指標(biāo)是反映食品衛(wèi)生品質(zhì)的指示性指標(biāo).展開(kāi)來(lái)看，最突出的問(wèn)題是飲料中桶裝飲用水銅綠假單胞菌不合格，其次是餐飲食品中復(fù)用餐飲具大腸菌群不合格以及肉制品、糕點(diǎn)食品中菌落總數(shù)不合格等.銅綠假單胞菌是一種條件致病菌，對(duì)于免疫力較弱的人群健康風(fēng)險(xiǎn)大，包裝飲用水銅綠假單胞菌不合格的主要原因是：（1）生產(chǎn)環(huán)境控制不嚴(yán)；（2）回收水桶消毒不徹底；（3）桶裝水密封不嚴(yán)等［32］.大腸菌群作為檢測(cè)食品是否被糞便污染，即是否被致病菌污染，食用后對(duì)人體健康造成危害程度的重要衛(wèi)生指標(biāo)，復(fù)用餐飲具大腸菌群不合格充分表明復(fù)用餐飲具消毒不達(dá)標(biāo)［33］.菌落總數(shù)是食品被微生物污染的程度，肉制品、糕點(diǎn)等食品中菌落總數(shù)不合格的原因有：（1）沒(méi)有嚴(yán)格控制生產(chǎn)加工過(guò)程的衛(wèi)生條件；（2）沒(méi)有合理使用防腐劑；（3）包裝材質(zhì)及厚度不合適；（4）貯存溫度不合適［34］.針對(duì)最為突出的桶裝飲用水銅綠假單胞菌問(wèn)題，該省監(jiān)管部門采取對(duì)桶裝水生產(chǎn)企業(yè)食品安全管理員開(kāi)展專題培訓(xùn)等一系列行之有效的措施［35］，因此從2019年起，微生物指標(biāo)不合格問(wèn)題已經(jīng)不具有顯著性.

農(nóng)藥獸藥指標(biāo)不合格的情況僅在2019年顯著高于平均水平，在其余年度沒(méi)有顯著性，可采取普通關(guān)注.農(nóng)藥獸藥指標(biāo)不合格的食品類別中，89.80%為食用農(nóng)產(chǎn)品、4.00%為茶葉及相關(guān)制品、3.60%為蜂產(chǎn)品、2.40%為肉制品和0.20%為餐飲食品；農(nóng)藥獸藥指標(biāo)不合格問(wèn)題中，55.80%為檢出禁用獸藥、25.00%為獸藥殘留超標(biāo)、15.60%為農(nóng)藥殘留超標(biāo)和3.60%為檢出禁用農(nóng)藥；禁用獸藥中，有69.18%為抗生素類藥物，究其原因主要是抗生素類獸藥價(jià)格便宜、效果好，養(yǎng)殖戶追求高的經(jīng)濟(jì)利益而使用.中華人民共和國(guó)農(nóng)業(yè)農(nóng)村部公告第194號(hào)已規(guī)定“自2020年7月1日起，飼料生產(chǎn)企業(yè)停止生產(chǎn)含有促生長(zhǎng)類藥物飼料添加劑（中藥類除外）的商品飼料.”［36］可見(jiàn)我國(guó)已采取源頭治理的方式避免抗生素類藥物在動(dòng)物性食品中的檢出.

食品中生物毒素，重金屬、無(wú)機(jī)和有機(jī)污染物與非食用物質(zhì)均為顯著低于年度平均水平的不合格指標(biāo)，表明這3類不合格情況在該省各年度抽檢結(jié)果中均較少發(fā)生，可予以一般關(guān)注.

3.4 時(shí)空分析

區(qū)別于王博遠(yuǎn)等［7］僅采用空間自相關(guān)分析研究食用油抽檢結(jié)果，王舟等［21］采用時(shí)空掃描分析麻痹性貝類毒素的監(jiān)測(cè)數(shù)據(jù)，本研究建立的時(shí)空掃描分析結(jié)果與各年度局部空間自相關(guān)分析結(jié)果可相互驗(yàn)證，得出時(shí)空掃描風(fēng)險(xiǎn)聚集區(qū)與局部空間自相關(guān)的合格率低聚集區(qū)基本一致的結(jié)論.2017年，廣安市呈合格率低-低聚集分布，且廣安市、南充市、遂寧市和達(dá)州市為次要聚集區(qū)Ⅱ；2018年，與自貢市和雅安市均接壤并位于2座城市之間的眉山市和樂(lè)山市為合格率低聚集區(qū)，且為主要聚集區(qū)；2019年，位于瀘州市和遂寧市之間的資陽(yáng)市和內(nèi)江市為合格率低聚集區(qū)，2020年，眉山市和資陽(yáng)市為合格率低-低聚集區(qū)，且2019和2020年，資陽(yáng)市和內(nèi)江市為次要聚集區(qū)Ⅰ.以上結(jié)果說(shuō)明，2017—2020年四川省食品安全風(fēng)險(xiǎn)聚集區(qū)從四川東部逐步往中部遷移，最終聚集在四川中偏東部.

3.5 監(jiān)管建議

食品抽檢的類別共有34大類，其細(xì)類超過(guò)450種，涉及的具體食品品種數(shù)以萬(wàn)計(jì)，每種食品檢測(cè)的指標(biāo)不同，且每年抽檢的具體食品品種和檢測(cè)指標(biāo)會(huì)根據(jù)年度計(jì)劃發(fā)生動(dòng)態(tài)調(diào)整，故本文以統(tǒng)計(jì)分析方法探索出問(wèn)題較為突出的食品類別和問(wèn)題項(xiàng)目后，以例舉形式繼續(xù)展開(kāi)深入分析，綜合得出該省存在的主要食品安全問(wèn)題有：（1）突出的問(wèn)題食品類別和突出的不合格指標(biāo)相互交織；（2）食品生產(chǎn)環(huán)節(jié)把控不嚴(yán)；（3）食品安全風(fēng)險(xiǎn)區(qū)域集中.根據(jù)本研究結(jié)果，提出加強(qiáng)該省食品安全的監(jiān)管建議：（1）根據(jù)食品類別和不合格指標(biāo)的綜合分析，建議加強(qiáng)對(duì)蔬菜制品中的食品添加劑SO2，蜂產(chǎn)品中的禁用獸藥諾氟沙星等，飲料中微生物指標(biāo)銅綠假單胞菌的分類監(jiān)管和專項(xiàng)整治，加大對(duì)問(wèn)題品種和指標(biāo)的食品抽檢力度；（2）建議對(duì)歷年抽檢合格率低于平均水平和總體合格率呈下降趨勢(shì)的23類食品的監(jiān)管進(jìn)行清單式管理，逐年解決部分問(wèn)題；（3）建議加強(qiáng)對(duì)食品生產(chǎn)環(huán)節(jié)的監(jiān)管力度，加強(qiáng)對(duì)生產(chǎn)加工者的宣傳引導(dǎo)和培訓(xùn)教育，從而提升從業(yè)人員的素質(zhì)，利用專家學(xué)者幫助生產(chǎn)企業(yè)改進(jìn)生產(chǎn)工藝，加強(qiáng)食品生產(chǎn)許可審查、體系檢查，從生產(chǎn)源頭解決問(wèn)題；（4）禁用藥物屢禁不止，多由于禁用藥物對(duì)食品性動(dòng)物的治病、防病效果好，但由于禁用藥物可能對(duì)人體產(chǎn)生致癌致畸的風(fēng)險(xiǎn)，建議監(jiān)管部門以多種形式鼓勵(lì)科研機(jī)構(gòu)開(kāi)展低毒獸藥的研究和應(yīng)用，做好養(yǎng)殖源頭治理；（5）針對(duì)食品生產(chǎn)者的人為故意導(dǎo)致食品抽檢不合格的，建議加大對(duì)違法處罰力度，起到震懾作用；（6）根據(jù)時(shí)空掃描的結(jié)果，建議重點(diǎn)加強(qiáng)對(duì)該省中偏東部地區(qū)的食品安全治理.

4 結(jié)束語(yǔ)

以四川省食品安全監(jiān)督抽檢結(jié)果為基礎(chǔ)，本文采用5種統(tǒng)計(jì)分析方法從食品類別、不合格項(xiàng)目、時(shí)間和地域等方面探索該省面臨的突出食品安全風(fēng)險(xiǎn)，創(chuàng)新性地建立了一種食品抽檢不合格情況的時(shí)空掃描分析模型，時(shí)空掃描分析結(jié)果可由局部空間自相關(guān)分析結(jié)果相驗(yàn)證，并發(fā)現(xiàn)該省的食品安全風(fēng)險(xiǎn)隨時(shí)間推移從東部向中部遷移的特征，為早期的食品安全風(fēng)險(xiǎn)預(yù)警提供基礎(chǔ)研究依據(jù).通過(guò)分析研究，科學(xué)識(shí)別該省食品安全問(wèn)題的高發(fā)品種和高發(fā)項(xiàng)目為蔬菜制品中的SO2、桶裝飲用水中的銅綠假單胞菌和蜂產(chǎn)品中檢出禁用獸藥，高發(fā)環(huán)節(jié)為生產(chǎn)環(huán)節(jié)，對(duì)科學(xué)制定食品抽檢計(jì)劃，進(jìn)一步解決食品安全重點(diǎn)難點(diǎn)問(wèn)題具有借鑒意義.期望在本研究的基礎(chǔ)上，進(jìn)一步收集環(huán)境污染、人口學(xué)、社會(huì)經(jīng)濟(jì)學(xué)和營(yíng)養(yǎng)膳食等方面的資料，采用多元的空間分析尋找其相互關(guān)系，將時(shí)間單位從年細(xì)化為月并建立預(yù)警模型，進(jìn)一步深入探討食品安全問(wèn)題的影響因素和解決辦法.