薛敏敏,劉芳芳?,朱 禮,袁列江,吳海智,陳 瑤
(湖南省產(chǎn)商品質(zhì)量檢驗研究院,湖南 長沙 410007)
豆奶粉是一種粉狀植物蛋白飲料,營養(yǎng)價值較高[1-2],既含有人體所必需的賴氨酸、蛋氨酸等多種氨基酸,又含有大量的不飽和脂肪酸與膳食纖維,以及多種礦物質(zhì)和維生素,且膽固醇含量較低[3],因此在我國居民食品構(gòu)成中占據(jù)重要地位。由于豆奶粉的主要原料之一大豆對鎳等重金屬具有富集作用,且豆制品生產(chǎn)加工過程也可能引入鎳元素,因此豆奶粉中存在鎳超標(biāo)的風(fēng)險。攝入過多鎳元素會造成心血管和腎臟疾病等多種危害[4-5],迫切需要通過準(zhǔn)確測定豆奶粉中鎳含量來保障我國食品安全。
豆奶粉中鎳含量測定可以采樣GB 5009.138—2017《食品安全國家標(biāo)準(zhǔn) 食品中鎳的測定》中規(guī)定的石墨爐原子吸收光譜法[6-7],但該方法穩(wěn)定性較差,易受載氣、樣品酸度等多種個因素的影響。GB 5009.268—2016《食品安全國家標(biāo)準(zhǔn) 食品中多元素的測定》規(guī)定了電感耦合等離子體質(zhì)譜法(inductively coupled plasma mass spectrometry,ICP-MS)[8-9]與電感耦合等離子體發(fā)射光譜法[10-11],其中ICP-MS具有靈敏度高、檢出限低、線性范圍寬等優(yōu)勢。然而國標(biāo)中使用的壓力罐消解法與微波消解法等樣品前處理技術(shù)普遍存在趕酸時間長、效率低等不足[12-13],超級微波消解具有空白低、污染少、酸加入量少等優(yōu)勢,且消解完全后無需趕酸,消解能力更強(qiáng)[14]。因此,超級微波消解-電感耦合等離子體質(zhì)譜法是測定豆奶粉中鎳含量極具潛力的方法。
不確定度是根據(jù)測定方法所用到的信息表征賦予被測量值分散性的非負(fù)參數(shù),能夠反映測量結(jié)果受測量誤差影響的不確定程度。豆奶粉中鎳的不確定度評定與優(yōu)化對于評估超級微波消解–電感耦合等離子體質(zhì)譜法的不確定程度、提升豆奶粉中鎳安全風(fēng)險評價的準(zhǔn)確度等具有現(xiàn)實意義。黃姍、韓偉立、楊愛君等研究了ICP-MS法測定食品中的重金屬的不確定度[15-17],馬玉婷、陳璐、康靜艷等具體探討了ICP-MS法測定食品中鎳的不確定度[18-20],侯建榮等研究了 ICP-MS法測定豆奶粉中鎘的不確定度[21],然而上述ICP-MS法的前處理技術(shù)主要采用國標(biāo)方法,且關(guān)于豆奶粉中鎳含量測定不確定度的相關(guān)研究尚未見報道。本文依據(jù)JJF1059.1—2012《測量不確定度評定與表示》與 CNAS—GL006:2019《化學(xué)分析中不確定度的評估指南》探討超級微波消解–電感耦合等離子體質(zhì)譜法測定豆奶粉中鎳含量的不確定度評定與優(yōu)化問題。一方面,建立不確定度評定的數(shù)學(xué)模型,分析不確定度的主要來源并評定不確定度分量,進(jìn)而計算合成相對標(biāo)準(zhǔn)不確定度與擴(kuò)展不確定度;另一方面,分析不同來源對不確定度的影響大小,并提出減小不確定度的可能優(yōu)化方法,為提升豆奶粉中鎳含量測定與安全風(fēng)險評價的準(zhǔn)確度提供參考依據(jù)。
豆奶粉樣品為隨機(jī)采集于當(dāng)?shù)爻械氖惺鄱鼓谭郏ㄆ放疲壕S維豆奶粉,生產(chǎn)廠家:維維食品飲料股份有限公司,生產(chǎn)日期:2022年 2月 21日,生產(chǎn)地點:江蘇徐州);豆奶粉中鎳質(zhì)控樣(ZKQC5210):批號P25295,河南標(biāo)準(zhǔn)物質(zhì)研發(fā)中心。
濃硝酸(痕量金屬級):美國Fisher化學(xué)公司;鎳標(biāo)準(zhǔn)儲備液(1 000 μg/mL):國家鋼鐵材料測試中心;內(nèi)標(biāo)儲備液:銠溶液(100 μg/mL):國家有色金屬及電子材料分析測試中心;質(zhì)譜調(diào)諧液:鈰、鈷、鋰、鈦、釔混合液(均為10 μg/mL):美國Agilent公司;實驗用水均為Milli-Q超純水。
UltraWAVE超級微波消解儀:北京萊伯泰科儀器股份有限公司;7700X電感耦合等離子體質(zhì)譜儀:美國 Agilent公司;BSA224S-CW電子天平:德國Sartorius公司;Milli-Q超純水處理系統(tǒng):美國 Millipore公司;100-1000 μL移液器、0.5-5 mL移液器:德國Eppendorf公司。
1.3.1 樣品前處理
準(zhǔn)確稱取0.5 g樣品(精確至0.001 g)于15 mL消解管中,加入3 mL濃硝酸,加蓋放置1 h后放入超級微波消解儀中,按照優(yōu)化后的升溫程序(見表1)進(jìn)行消解。消解完畢冷卻后取出,并用少量超純水多次洗滌消解管,轉(zhuǎn)移至50 mL容量瓶中,用超純水稀釋至刻度,混勻備用,同時制備試劑空白。
表1 超級微波消解升溫程序Table 1 Heating program of ultra-wave digestion
1.3.2 標(biāo)準(zhǔn)工作溶液的配置
準(zhǔn)確吸取0.1 mL 鎳標(biāo)準(zhǔn)儲備溶液,用5%硝酸稀釋定容至100 mL,得到質(zhì)量濃度為1 000 μg/L的鎳標(biāo)準(zhǔn)使用液。分別移取鎳標(biāo)準(zhǔn)使用溶液0.00、1.00、2.00、3.00、4.00、5.00 mL于100 mL容量瓶中,用 5%硝酸定容至刻度,得到質(zhì)量濃度為0.0、10.0、20.0、30.0、40.0、50.0 μg/L 的標(biāo)準(zhǔn)工作液。
1.3.3 樣品測定
參考 GB 5009.268—2016《食品安全國家標(biāo)準(zhǔn) 食品中多元素的測定》第一法對消解液進(jìn)行測定,其中 ICP-MS工作條件見表2。選取鎳的質(zhì)量數(shù)60,采用八級桿氦模式及四級桿質(zhì)量過濾器消除相應(yīng)的質(zhì)譜型干擾,在ICP-MS的動能歧視(kinetic energy discrimination,KED)模式下利用碰撞/反應(yīng)池技術(shù)除去多原子離子的干擾,同時采用銠為內(nèi)標(biāo)元素消除非質(zhì)譜干擾。
表2 ICP-MS工作條件Table 2 Working conditions of ICP-MS
1.3.4 鎳含量的數(shù)學(xué)模型
豆奶粉中鎳含量的測定結(jié)果按以下公式計算:
式中,X為豆奶粉中鎳的含量,mg/kg;ρ與0ρ分別表示豆奶粉樣品消解液與空白溶液中鎳的質(zhì)量濃度,μg/L;V為樣品消解液定容體積,mL;f為稀釋倍數(shù);m為樣品稱取質(zhì)量,g。
利用超級微波消解–電感耦合等離子質(zhì)譜法測定豆奶粉中的鎳含量,測定方法的標(biāo)準(zhǔn)曲線為Y=0.0830X+0.0114,在質(zhì)量濃度 0~50 μg/L 范圍內(nèi)線性關(guān)系良好,相關(guān)系數(shù)為0.999 9。方法的檢出限為 0.002 mg/kg,加標(biāo)回收率為 91.37%~106.46%,相對標(biāo)準(zhǔn)偏差為1.19%~3.22%。
同時,按照與豆奶粉樣品相同的方法,對豆奶粉中鎳質(zhì)控樣(ZKQC5210)進(jìn)行消解和測定,質(zhì)控樣品的測定值為2.23 mg/kg,測定結(jié)果在標(biāo)定范圍(2.24±0.40) mg/kg之內(nèi),表明方法準(zhǔn)確度較高。
超級微波消解–電感耦合等離子質(zhì)譜法測定豆奶粉中鎳的不確定度來源依照測定過程具體可分為6個方面,具體見表3,且根據(jù)不確定度來源的特點,分別選取不同的評定分類,其中A類評定是指對獨立重復(fù)測定的量值用統(tǒng)計分析方法進(jìn)行不確定度的評定,B類評定主要是指通過判斷量值的可能值區(qū)間與概率分布對不確定度進(jìn)行的評定。
表3 不確定度來源與評定分類Table 3 Sources and evaluation types of uncertainty
2.3.1 樣品稱重
測定使用的天平為特種準(zhǔn)確度級,以其檢定證書標(biāo)示的結(jié)果為依據(jù)值確定不確定度。
天平偏載引入的標(biāo)準(zhǔn)不確定度采用B類評定方法,證書規(guī)定的最大允許誤差為±0.5 mg,概率采用矩形分析,即包含因子取,則標(biāo)準(zhǔn)不確定度為。
天平重復(fù)性引入的標(biāo)準(zhǔn)不確定度也采用B類評定方法,由證書規(guī)定的最大允許誤差±0.1 mg以及矩形分布對應(yīng)的包含因子,計算其標(biāo)準(zhǔn)不確定度為。
天平示值引入的標(biāo)準(zhǔn)不確定度同樣采用B類評定方法最大允許誤差為±0.2 mg,包含因子取,則標(biāo)準(zhǔn)不確定度為。
樣品稱重量為0.5 g,需要2次獨立稱重(去皮與樣重),因此樣品稱重的相對標(biāo)準(zhǔn)不確定度為
2.3.2 消解液定容
JJG196—2006《常用玻璃量器檢定規(guī)程》規(guī)定在20 ℃時50.0 mL容量瓶允差為±0.05 mL,采用三角分布,即,則其體積標(biāo)準(zhǔn)不確定度為。同時,考慮溫度影響,這里采用水的體積膨脹系數(shù)β= 0 .000 208℃-1近似溶液的體積膨脹系數(shù),實驗室溫度變化為±4 ℃,并采用矩形分布,即,則其50 mL容量瓶溫度標(biāo)準(zhǔn)不確定度為。因此,消解液定容引入的相對標(biāo)準(zhǔn)不確定度為
2.3.3 消解不均勻性
利用超級微波消解儀對加入硝酸后的樣品進(jìn)行消解時,儀器的功率、溫度、時間、壓力等因素可能影響樣品消解的均勻性,進(jìn)而引入不確定度。加標(biāo)回收實驗的統(tǒng)計結(jié)果見表4,由于每組的實驗次數(shù)相同,則合并標(biāo)準(zhǔn)偏差為sp=(2.37% + 1 .19% + 3 .22%)/3 = 2 .26%,此時標(biāo)準(zhǔn)不確定度為。由表4還可計算合并均值:
表4 加標(biāo)回收實驗結(jié)果統(tǒng)計(n=6)Table 4 Statistics results of standard additionexperiments(n=6)
因此相對標(biāo)準(zhǔn)不確定度為
2.3.4 標(biāo)準(zhǔn)工作液的配置
2.3.4.1 鎳標(biāo)準(zhǔn)儲備液 鎳標(biāo)準(zhǔn)儲備液證書給出鎳濃度值與擴(kuò)展不確定度分別為1 000 μg/mL與10 μg/mL,其中k=2,則由鎳標(biāo)準(zhǔn)儲備液引入的相對標(biāo)準(zhǔn)不確定為
2.3.4.2 鎳標(biāo)準(zhǔn)使用液 鎳標(biāo)準(zhǔn)使用液的不確定度主要源自鎳標(biāo)準(zhǔn)儲備液、溶液移取與溶液定容。溶液移取使用1 000 μL的移液器,在20 ℃時其最大允許誤差為±0.001 mL,故按照矩形分布時其體積標(biāo)準(zhǔn)不確定度為;同時考慮溫度影響,參考 2.2.2消解液定容可知,1 000 μL的移液器的溫度標(biāo)準(zhǔn)不確定度為,因此相對標(biāo)準(zhǔn)不確定度為
溶液定容參考JJG196—2006《常用玻璃量器檢定規(guī)程》,其中規(guī)定在20 ℃時100.0 mL容量瓶允差為±0.1 mL,同樣采用三角分布,即,則其體積標(biāo)準(zhǔn)不確定度為。同時,100 mL容量瓶的溫度標(biāo)準(zhǔn)不確定度為。因此,溶液定容引入的相對標(biāo)準(zhǔn)不確定度為因此綜合上述三個方面可計算鎳標(biāo)準(zhǔn)使用液的相對標(biāo)準(zhǔn)不確定度:
2.3.4.3 不同濃度的鎳標(biāo)準(zhǔn)工作液 不同濃度的鎳標(biāo)準(zhǔn)工作液由鎳標(biāo)準(zhǔn)使用液經(jīng)稀釋定容配制而成,因此其不確定度主要源于鎳標(biāo)準(zhǔn)使用液、溶液移取的移液管以及溶液定容的容量瓶。鎳標(biāo)準(zhǔn)使用液的相對標(biāo)準(zhǔn)不確定度已經(jīng)計算獲得,溶液定容采用的是100 mL容量瓶,其相對標(biāo)準(zhǔn)不確定度在計算鎳標(biāo)準(zhǔn)使用液的相對標(biāo)準(zhǔn)不確定度時也已經(jīng)計算獲得,即urel(V100)=0.00063,因此下面主要考慮移液管引入的不確定度。
溶液移取使用的5 mL移液器,其在不同檢定點的允許誤差如表5所示,引入的不確定度主要考慮移液器的體積標(biāo)準(zhǔn)不確定度與溫度標(biāo)準(zhǔn)不確定度,其中概率采用矩形分布,即,溶液的體積膨脹系數(shù)近似采用水的體積膨脹系數(shù)β= 0 .000 208℃-1,同時實驗室溫度變化為±4 ℃,具體分析結(jié)果見表5,分別給出了移液器對不同濃度的鎳標(biāo)準(zhǔn)工作溶液引入的相對標(biāo)準(zhǔn)不確定度。
表5 移液器引入的相對標(biāo)準(zhǔn)不確定度Table 5 Relative standard uncertainty introduced by pipette
按照相對標(biāo)準(zhǔn)不確定度的合成公式,不同濃度的鎳標(biāo)準(zhǔn)工作液的相對標(biāo)準(zhǔn)不確定度分別為:
此時,鎳標(biāo)準(zhǔn)工作液引入的相對標(biāo)準(zhǔn)不確定度為
2.3.5 標(biāo)準(zhǔn)曲線的擬合
標(biāo)準(zhǔn)曲線擬合時,針對每個濃度的鎳標(biāo)準(zhǔn)工作液進(jìn)行3次平行測定,共進(jìn)行p=6組測定,此時n=3×6= 18。經(jīng)過最小二乘擬合求得線性方程為Y=0.0830X+0.0114,其中相關(guān)系數(shù)為 0.999 9。測定結(jié)果以及初步數(shù)據(jù)統(tǒng)計結(jié)果見表6,則據(jù)此計算標(biāo)準(zhǔn)曲線的擬合剩余標(biāo)準(zhǔn)差為
表6 鎳標(biāo)準(zhǔn)工作溶液測定結(jié)果及統(tǒng)計Table 6 Statistics and determination results of nickel standard working solution
式中,Yi,j為第i組測定中第j個平行的測定結(jié)果,Xi為第i組測定的鎳標(biāo)準(zhǔn)工作液的濃度。樣品中鎳含量的平均濃度可由表7中測定結(jié)果求得,即= 3 8.25 μg/L,且鎳標(biāo)準(zhǔn)工作液濃度的均值為= 2 5.00 μg/L,因此標(biāo)準(zhǔn)曲線擬合引入的標(biāo)準(zhǔn)不確定度為
表7 豆奶粉中鎳含量測定結(jié)果Table 7 Determination results of nickel content in soy milk powder
則標(biāo)準(zhǔn)曲線擬合引入的相對標(biāo)準(zhǔn)不確定度為
2.3.6 樣品測定重復(fù)性
利用超級微波消解–電感耦合等離子體質(zhì)譜法重復(fù)測定豆奶粉樣品中鎳含量6次,即m=6,測定結(jié)果見表7,則按照A類評定方法,首先計算樣品中鎳含量的標(biāo)準(zhǔn)差:
則鎳含量的標(biāo)準(zhǔn)不確定度為
因此樣品重復(fù)性引入的相對標(biāo)準(zhǔn)不確定度為
依據(jù)不確定度分量的評定,測定方法的相對標(biāo)準(zhǔn)不確定度為
則合成標(biāo)準(zhǔn)不確定度為u=urel×=0.1111 mg/kg,并取置信概率為95%,即包含因子為k=2,擴(kuò)展不確定度為U=k×u= 0.222 2 mg/kg。因此超級微波消解–電感耦合等離子體質(zhì)譜法測定豆奶粉中的鎳含量為:(3.825±0.222) mg/kg(k=2)。
基于超級微波消解–電感耦合等離子體質(zhì)譜法測定豆奶粉中鎳含量的不確定度計算過程,不同不確定度分量的相對標(biāo)準(zhǔn)不確定度如圖1所示,可知影響測定不確定度的主要因素由強(qiáng)到弱依次為:標(biāo)準(zhǔn)工作溶液的配置>標(biāo)準(zhǔn)曲線的擬合>消解不均勻性>樣品測定重復(fù)性>樣品稱重>消解液定容。盡管降低實驗儀器與設(shè)備的允許誤差以及提升試劑與材料的含量精度是降低方法不確定度的根本途徑,但改進(jìn)測定方法也可以實現(xiàn)不確定度的優(yōu)化。下面在不改變實驗儀器設(shè)備與試劑材料的前提下,探討不確定度優(yōu)化的可能途徑。
圖1 不同不確定度分量的相對標(biāo)準(zhǔn)不確定度Fig.1 Relative standard uncertainty of different uncertainty components
2.5.1 優(yōu)化標(biāo)準(zhǔn)工作液的配置
標(biāo)準(zhǔn)工作溶液的配置主要包括由標(biāo)準(zhǔn)儲備液配置標(biāo)準(zhǔn)使用液,以及由標(biāo)準(zhǔn)使用液配置標(biāo)準(zhǔn)工作液兩個過程。通過分析不確定度分量可知,溶液移取是溶液配置中引入不確定度的主要來源,特別是當(dāng)移取的溶液體積較小時,引入的不確定度較大。因為由標(biāo)準(zhǔn)儲備液配置標(biāo)準(zhǔn)使用液過程僅移取了0.1 mL溶液,因此主要通過優(yōu)化由標(biāo)準(zhǔn)儲備液配置標(biāo)準(zhǔn)使用液的過程,降低方法不確定度。
由標(biāo)準(zhǔn)儲備液配置標(biāo)準(zhǔn)使用液的原過程為:準(zhǔn)確吸取0.1 mL鎳標(biāo)準(zhǔn)儲備溶液,用5%硝酸稀釋定容至100 mL,得到質(zhì)量濃度為1 000 μg/L的鎳標(biāo)準(zhǔn)使用液。優(yōu)化后的配置過程為:準(zhǔn)確吸取1 mL 標(biāo)準(zhǔn)儲備溶液,用5%硝酸稀釋定容至100 mL,得到質(zhì)量濃度為10 000 μg/L的標(biāo)準(zhǔn)中間液;再準(zhǔn)確吸取10 mL標(biāo)準(zhǔn)中間液,用5%硝酸稀釋定容至100 mL,得到質(zhì)量濃度為1 000 μg/L的鎳標(biāo)準(zhǔn)使用液。
計算優(yōu)化后引入的不確定度。吸取1 mL標(biāo)準(zhǔn)儲備溶液使用1 000 μL的移液器,在20 ℃時其最大允許誤差為±0.001 mL;吸取10 mL標(biāo)準(zhǔn)中間液使用10 mL移液器,在20 ℃時其最大允許誤差為±0.03 mL。綜合考慮定容等引入的不確定度,計算標(biāo)準(zhǔn)使用液的相對標(biāo)準(zhǔn)不確定度urel(N2)由0.007 7下降至0.005 2,進(jìn)一步計算標(biāo)準(zhǔn)工作液的相對標(biāo)準(zhǔn)不確定度urel(N3)由0.020 5下降至0.015 5。
2.5.2 優(yōu)化標(biāo)準(zhǔn)曲線的擬合
由標(biāo)準(zhǔn)曲線擬合的不確定度計算公式可知,在標(biāo)準(zhǔn)曲線及其擬合剩余標(biāo)準(zhǔn)差不變的前提下,增加每個濃度標(biāo)準(zhǔn)工作液的平行測定次數(shù)可以有效降低標(biāo)準(zhǔn)曲線擬合引入的不確定度。例如將平行測定次數(shù)增加到6次,則標(biāo)準(zhǔn)曲線的擬合引入的不確定度urel(S)由0.016 7下降至0.015 9。
調(diào)整標(biāo)準(zhǔn)工作液的濃度使得其濃度均值與樣品濃度均值相接近,會對降低標(biāo)準(zhǔn)曲線擬合的不確定度產(chǎn)生積極作用,但是會影響標(biāo)準(zhǔn)工作液的配置引入的不確定度,因此需要具體分析以實現(xiàn)合成標(biāo)準(zhǔn)不確定度的下降。
2.5.3 優(yōu)化消解不均勻性
由消解不均勻性引入的不確定度計算公式可知,增加每組加標(biāo)回收實驗的平行實驗次數(shù),在平均回收率與平均實驗偏差基本不變的前提下可以降低其不確定度。當(dāng)平行次數(shù)由6次增加到12次時,消解不均勻性引入的不確定度由 0.009 4下降至0.006 7。
2.5.4 優(yōu)化樣品測定重復(fù)性
由樣品測定重復(fù)性引入的不確定度計算公式可知,在鎳含量均值與標(biāo)準(zhǔn)差基本不變的前提下,增加樣品的測定次數(shù)能夠降低其不確定度。當(dāng)測定次數(shù)從6次增加到12次時,樣品測定重復(fù)性引入的不確定度由0.007 4下降至0.005 2。
2.5.5 不確定度優(yōu)化結(jié)果
采取上述不確定度優(yōu)化方法后,超級微波消解–電感耦合等離子體質(zhì)譜法測定豆奶粉中鎳含量的相對標(biāo)準(zhǔn)不確定度由0.029 0下降至0.023 8。特別需要注意的是上述優(yōu)化措施僅考慮測定方法的改進(jìn),如果同時采用允許誤差更小的儀器設(shè)備以及含量精度更高的材料試劑,相對標(biāo)準(zhǔn)不確定度還會可進(jìn)一步下降,在此不再贅述。
本研究對超級微波消解–電感耦合等離子體質(zhì)譜法測定豆奶粉中鎳的不確定度進(jìn)行了評定與優(yōu)化,分析了不確定度的來源并評定了其不確定度分量,進(jìn)而計算了合成相對標(biāo)準(zhǔn)不確定度與擴(kuò)展不確定度。結(jié)果表明,豆奶粉中鎳含量的擴(kuò)展不確定度為 0.222 mg/kg,即測定結(jié)果為(3.825±0.222) mg/kg(k=2)。同時影響不確定度的主要因素由強(qiáng)到弱依次為:標(biāo)準(zhǔn)工作溶液的配置>標(biāo)準(zhǔn)曲線的擬合>消解不均勻性>樣品測定重復(fù)性>樣品稱重>消解液定容,因此優(yōu)化標(biāo)準(zhǔn)工作溶液的配置與標(biāo)準(zhǔn)曲線的擬合是不確定度的主要來源,通過改進(jìn)溶液移取方式、增加平行測定次數(shù)等方式可以將不確定度由0.029 0下降至0.023 8。通過對測定方法不確定度的評定與優(yōu)化,可為提升豆奶粉中鎳含量測定與安全風(fēng)險評價的準(zhǔn)確度提供參考依據(jù)。