信息多維采集的農藥殘留大數據評估市售果蔬安全水平
|
|
風險評估體系的構建
風險評估模型
果蔬農藥殘留膳食暴露風險評估模型
本文采用我國民眾平均體重（bw）60?kg，每日攝入食品（Fi）0.38 kg/（人·天），以及本研究組?2012—2015?年對全國果蔬的農藥殘留偵測結果（Ri），計算農藥?c?的食品安全指數?IFSc，并以此表征該果蔬的膳食風險程度：
其中：EDIc為農藥?c?的實際日攝入量估算值；EDIc=∑(Ri×Fi×Ei×Pi)，i?為食品種類，Ri為食品?i?中農藥?c?的殘留水平（單位：mg/kg），Fi為食品?i?的日攝入量（單位：g/(人·天)），Ei為食品?i?的可食用部分因子（Ei=1），Pi為食品?i?的加工處理因子（Pi=1）；SIc為安全攝入量，采用《食品安全國家標準食品中農藥最大殘留限量》（GB 2763-2016）中的?ADI?值，此時校正因子?f=1。
值?IFSc≤0.1、∈(0.1, 1]?和?>1?分別代表含農藥?c?的果蔬的膳食暴露風險為：對人體沒有影響、影響可接受和影響不可接受。當果蔬中存在多種農藥時，使用?IFSc的平均值?IF—S?評價果蔬安全情況，即該水果或蔬菜的膳食風險程度。值?IF—S??≤0.1、∈(0.1, 1]?和?>1，分別代表某一果蔬樣品對人體的膳食暴露風險程度為：沒有影響、影響可接受和影響不可接受。
果蔬農藥殘留超標預警風險評估模型
本研究進而計算了農藥殘留的預警風險系數?R：
其中，P?為該種危害物的超標率；F?為對危害物的施檢頻率，本文涉及的偵測數據?F=1；S?為危害物的敏感因子；a、b?分別為相應的權重系數。據相關資料，本次評價中?S=1，a=100，b=0.1[21,28]。R?值?≤1.5、∈(1.5, 2.5]?和?＞2.5?分別表示該種危害物的預警風險程度為低度、中度和高度。
計算?P?時，首先判斷農藥是否為禁用農藥。若為禁用農藥，檢出即為超標，則?P=?能檢出的樣品數/總樣品數。若為非禁用農藥，則?P=?超標的樣品數/總樣品數。超標是指危害物含量高于最大殘留限量標準值（MRLs），本研究中?MRLs?又分別采用中國國家標準和歐盟標準，這一方面因為中國的?MRLs?值尚不健全必須引入歐盟標準才可全面計算，另一方面二者對比以期支撐我國已有?MRLs?標準的進一步優化。
風險評估軟件開發
軟件開發基本要求
膳食暴露風險評估專業軟件。軟件應能分別計算?GC-Q-TOF/MS?和?LC-Q-TOF/MS?方法偵測出的果蔬樣品中每種農藥的?IFSc值，進而按?IFSc值大小排序，重點分析膳食暴露風險高的樣本，并對風險程度進行多維度解析。
超標預警風險評估專業軟件。分別采用?MRLs?的中國國家標準和歐盟標準，逐個計算不同果蔬中不同農藥的風險系數?R，并將結果按照禁用農藥和非禁用農藥分別列表，對風險程度進行多維度解析。
軟件開發方法
Python?語言是兼具解釋性、編譯性、互動性和面向對象的腳本語言，本研究選用?Python?語言進行編程。風險評估專用程序主要功能包括：分別讀入每例樣品?GC-Q-TOF/MS?和?LC-Q-TOF/MS?偵測的農藥殘留狀況數據；根據風險評估工作要求，對不同農藥種類、限量標準、果蔬種類、地域等的?IFSc值和?R?值進行計算；與?MRLs?的中國國家標準和歐盟標準限值進行對比，篩選出超標農藥，并單獨分析；再分別按不同農藥、不同果蔬種類分類處理，編寫自動計算和排序的計算機代碼；最后將生成的膳食暴露風險評估和超標預警風險評估結果列入設計好的各個表格中；定性判斷風險對目標的影響程度。
將本研究開發出的上述定制軟件命名為“風險值自動計算—信息多維采集分析”專業軟件。該軟件通過多維大數據融合技術，自動獲得不同維度的果蔬膳食暴露風險和超標預警風險評估結果。