何桂華，陳晞 孫志成 張彤羽

（濟南出入境檢驗檢疫局，濟南 250014） （山東省實驗中學，濟南 250000） （青島第二中學，山東青島 266071）

應用多元統(tǒng)計技術研究動物產(chǎn)品中重金屬的分布特征

何桂華，陳晞 孫志成 張彤羽

（濟南出入境檢驗檢疫局，濟南 250014） （山東省實驗中學，濟南 250000） （青島第二中學，山東青島 266071）

采用直方圖頻率說明、參數(shù)分布擬合假設檢驗、非參數(shù)密度估計等多元統(tǒng)計技術，研究了鉛、鎘、汞、砷在不同動物源性產(chǎn)品中的分布特征，通過非參數(shù)方法給出了4種元素在不同產(chǎn)品中的密度分布函數(shù)曲線和經(jīng)驗分布函數(shù)曲線，以此解釋了元素的分布特征，掌握了重金屬的累積富集規(guī)律。

多元統(tǒng)計；重金屬；動物產(chǎn)品；分布特征

食品安全是關系到人民健康和國計民生的重大問題，也是社會關注的熱點、焦點問題。重金屬廣泛存在于水、空氣、土壤和各種生活、生產(chǎn)資料中，并通過食物鏈在人體內(nèi)累積富集，嚴重影響人們的身體健康。目前，我國已建立比較完善的動物產(chǎn)品殘留監(jiān)控體系，每年在殘留監(jiān)控檢測獲得大量的檢測數(shù)據(jù)，但所獲得的數(shù)據(jù)都是孤立、零散的，無法科學準確地揭示殘留污染的規(guī)律性。

目前，有關食品中有害元素的報道大多數(shù)仍停留在檢測方法的改進和樣品的分析檢測方面［1–4］，應用多元統(tǒng)計技術研究蔬菜、蜂蜜中的重金屬已有報道［5–7］，但研究不同動物產(chǎn)品中重金屬的分布特征尚未見報道。筆者應用多元統(tǒng)計技術研究了不同動物產(chǎn)品重金屬的分布特征，揭示重金屬的累積富集規(guī)律，有利于發(fā)揮檢測數(shù)據(jù)對檢驗檢疫監(jiān)管的有效指導，提高產(chǎn)品質量。

1 實驗部分

1.1 主要儀器與試劑

石墨爐原子吸收光譜儀：SPECTRAA 220Z型，美國Varian公司；

原子熒光光譜儀：AFS-2202型，北京海光儀器公司；

微波消解儀：Mars5型，美國CEM公司；

鉛、鎘、汞標準溶液：濃度均為1 000 μg／mL，國家標準物質研究中心；

硝酸：優(yōu)級純;

基體改進劑磷酸銨溶液：20 g／L，稱取磷酸銨2.0 g溶于水中，稀釋至100 mL；

氫氧化鉀溶液：2 g／L，稱取氫氧化鉀2 g溶于水中，稀釋至1 000 mL；

硼氫化鉀溶液：7 g／L，稱取3.5 g硼氫化鉀溶于500 mL 2 g／L氫氧化鉀溶液中；

實驗用水均為二次去離子水。

1.2 測定方法

依據(jù)GB／T 5009.11–2003，GB 5009.12–2010，GB／T 5009.15–2003，GB／T 5009.17–2003［8–11］分別測定樣品中的無機砷、鉛、鎘和汞。其中鉛、鎘采用石墨爐原子吸收光譜儀測定，汞和無機砷采用原子熒光光譜儀測定。

1.3 研究對象與數(shù)據(jù)來源

研究對象：兔腎、雞腎、豬腸衣、羊腸衣、蝦、甲殼類、鰻魚。

被測重金屬元素：鉛、鎘、汞、無機砷。

研究數(shù)據(jù)：山東省3年來出口動物產(chǎn)品殘留監(jiān)控檢測結果。

2 結果與討論

2.1 鉛、鎘、汞、砷在腎中的分布特征研究

以雞腎、兔腎為例，通過數(shù)據(jù)的均值、方差、標準差等數(shù)字特征和用參數(shù)或非參數(shù)的方法給出概率密度或分布密度函數(shù)直方圖的形式研究鉛、鎘、汞、砷在腎中的分布特征，由此給出4種重金屬在產(chǎn)品腎中的分布規(guī)律，為動物產(chǎn)品檢驗監(jiān)管提供依據(jù)。

2.1.1 鉛在腎中的分布特征

（1）對雞腎、兔腎中的鉛含量數(shù)據(jù)進行數(shù)理統(tǒng)計分析，考察其樣本均值、標準偏差（s)、樣本方差、極大值、中位數(shù)等統(tǒng)計量，結果見表1。SAS數(shù)據(jù)分析軟件輸出結果見圖1。

表1 鉛在雞腎、兔腎中的分布情況(n=93)

由表1可知，雞腎、兔腎中鉛含量的樣本均值為555.40 μg／kg，標準差為2 121.41 μg／kg，最大值為18 000μg／kg，遠高于限量值500 μg／kg，說明存在嚴重不合格樣品。在所檢測的93個樣品中，絕大多數(shù)樣品鉛含量小于500 μg／kg，合格率為91.4%，不合格樣品8個，不合格率8.6%，并且超標樣品鉛含量均值為5 410 μg／kg，遠大于限量值500 μg／kg，超標嚴重，說明鉛比較容易通過食物鏈在雞、兔腎臟中富集，應重點加強監(jiān)控檢測。

圖1 腎中鉛分布的直方圖及指數(shù)、正態(tài)分布擬合曲線圖

（2）為進一步研究雞腎、兔腎中鉛含量分布，繪制雞腎、兔腎中鉛含量分布直方圖及分布擬合曲線圖，如圖1所示。觀察頻率直方圖及擬合曲線圖，根據(jù)經(jīng)驗假定雞腎、兔腎中鉛含量近似服從正態(tài)或指數(shù)分布，并對其參數(shù)分布函數(shù)（正態(tài)分布、指數(shù)分布）進行假設檢驗，檢驗結果見表2。

從直方圖中可看出，正態(tài)分布密度函數(shù)和指數(shù)分布密度函數(shù)均不能較好地擬合鉛分布。表2中鉛含量正態(tài)分布假設檢驗和指數(shù)分布假設檢驗p值結果均遠小于檢驗水平0.05，因此可以判定鉛含量不服從正態(tài)分布和指數(shù)分布。

表2 腎中鉛含量正態(tài)、指數(shù)分布假設檢驗結果

（3）采用非參數(shù)統(tǒng)計方法，給出腎中鎘、汞、砷含量分布的密度函數(shù)估計。假定隨機變量X具有密度函數(shù)f(x)，則可求其數(shù)學期望（均值）方差也可計算隨機變量X落在任一區(qū)間[a,b]的概率因此給出鉛含量的分布密度函數(shù)估計，就可以得到對其感興趣的分布特征。圖2和圖3分別為腎中鉛的密度函數(shù)估計曲線和經(jīng)驗分布函數(shù)曲線，此二曲線描述了腎中鉛含量的分布規(guī)律。從密度函數(shù)估計曲線可以預測元素鉛在任一濃度下的密度或發(fā)生幾率，而從經(jīng)驗分布函數(shù)曲線則可以得到元素鉛在任一濃度下的累積概率。

圖2 腎中鉛含量密度函數(shù)估計曲線

2.1.2 鎘、汞、砷在腎中的分布特征

同上述方法，研究了腎中鎘、汞、砷含量的分布，通過統(tǒng)計分析可知，兔腎樣本61個，鎘含量超標樣本11個,不合格率為18%，超標樣品鎘的平均含量為7 575 μg/kg，是限量值的15倍，說明鎘極易在兔腎中富集，應重點加強監(jiān)控；腎中汞、砷和雞腎中鎘含量較低，合格率為100%；3種元素的分布特征均不能用參數(shù)分布合理表示，采用非參數(shù)統(tǒng)計方法給出了3種元素的密度函數(shù)估計曲線和經(jīng)驗分布函數(shù)曲線，由曲線可以計算3種元素在任一濃度下的分布密度和累積概率。

圖3 腎中鉛含量經(jīng)驗分布函數(shù)曲線

2.2 鉛、鎘、汞、砷在腸衣中的分布特征研究

腸衣中元素鉛容易富集，在所檢測的57個腸衣樣品中，10個樣品鉛超標，不合格率為17.5%。從鉛含量分布密度估計可以直觀看到鉛的分布特征。鎘、汞、無機砷含量均較低，沒有超標樣品，合格率為100%，從密度函數(shù)估計曲線可以預測每個元素在任一濃度下的密度或發(fā)生幾率，從經(jīng)驗分布函數(shù)曲線可以得到每種元素在任一濃度下的累積概率。

2.3 鉛、鎘、汞、砷在水產(chǎn)品中的分布特征

大多數(shù)水產(chǎn)品中鉛、鎘含量較低，鉛不合格樣品1個（養(yǎng)殖蝦），不合格率為2.3%；鎘不合格樣品2個（蝦、甲殼類各1個），不合格率為4.6%；水產(chǎn)品中汞、無機砷含量低，無超標樣品。利用直方圖可直觀地表征4種元素的分布特征，采用非參數(shù)統(tǒng)計方法得到4種元素的密度函數(shù)曲線和經(jīng)驗分布函數(shù)曲線，從中可以計算水產(chǎn)品中各元素在任一濃度點的密度和累積概率。

3 結語

采用直方圖頻率說明、參數(shù)分布擬合假設檢驗、非參數(shù)密度估計等多元統(tǒng)計技術，研究了鉛、鎘、汞、砷在不同動物源性產(chǎn)品中的分布特征，通過非參數(shù)方法給出了4種元素在不同產(chǎn)品中的分布密度函數(shù)曲線和經(jīng)驗分布函數(shù)曲線，良好地解釋了元素的分布特征，為深入研究動物源性產(chǎn)品中重金屬的防控技術、優(yōu)化動物源性產(chǎn)品的監(jiān)督管理模式、提高動物源性產(chǎn)品質量提供了良好的技術支持和決策依據(jù)。

［1］朱芳，劉翅，楊中藝.微波消解—原子熒光光譜法測定番茄不同部位中微量汞［J］.分析測試學報，2005(1):89–92.

［2］夏偉，謝濤.論用ICP–MS對食品中重金屬元素的測定［J］.大眾科技，2005(6):36–37.

［3］黃志勇，楊妙峰，劉波.茶葉和蔬菜中鉛的同位素稀釋電感耦合等離子體質譜的測定［J］.分析實驗室，2005(6):65–68.

［4］Yvers Roggo，Ludovic Duponchel，Jean-pierre Huvenne.Quality Evaluation of Sugar Beet （Beta Vulgaris）by Near-Infrared Spectroscopy［J］. J Agric Food chem，2004，52(5):1055–1061.

［5］柳勇，何江華，王少毅，等，廣州市蔬菜地重金屬劑量對蔬菜富集重金屬的影響——以菜心為例［J］.生態(tài)環(huán)境，2003，12(3):273–276.

［6］Devillers J，Dore J C，marenco M，et al. Chemometrical aanalysis of 18 metallic and nonmetallic elements Found in honeys sold in France［J］. J Argric Food chem，2002，50(21):5 998–6 007.

［7］程繼滿，周韞韜，劉羿，等.海南島龍灣港濱珊瑚重金屬含量的年際變化［J］.海洋地質與第四紀地質，2005，25(1):11–18.

［8］GB／T 5009.11–2003 食品中總砷和無機砷的測定[S].

［9］GB 5009.12–2010 食品中鉛的測定[S].

［10］GB／T 5009.15–2003 食品中鎘的測定[S].

［11］GB／T 5009.17–2003 食品中總汞及有機汞的測定[S].

新型金屬玻璃場發(fā)射電極問世 可用于電鏡等多種儀器

據(jù)美國物理學家組織網(wǎng)報道，澳大利亞科學家在最新一期的《物理評論快報》雜志上報告稱，他們研制出一種屬性與玻璃類似的新型金屬化合物，并用其替代塑料與碳納米管結合制成新的場發(fā)射電極。該場發(fā)射電極能制造出穩(wěn)定的電子束，有望用在消費電子和電子顯微鏡等領域。

以前，科學家們主要通過將碳納米管和其它納米材料內(nèi)嵌于塑料中來制造場發(fā)射電極。這些場發(fā)射電極盡管種類繁多且容易制造，擁有很大應用潛力，但其瑕疵也很多，比如，塑料的導電能力太弱;塑料的熱穩(wěn)定性很低，無法對抗長時間操作產(chǎn)生的大量熱量。

現(xiàn)在，澳大利亞莫納什大學的科研團隊和澳大利亞聯(lián)邦科學與工業(yè)研究組織下屬的過程科學和工程研究院的科學家攜手，研發(fā)出了一種新的應用潛力很大且容易制造的材料——非晶塊金屬玻璃(ABM)，并用其代替塑料制造出場發(fā)射電極。當這些非晶塊金屬玻璃合金冷卻時會形成非晶材料，讓它們的一舉一動更像玻璃。

這種非晶塊金屬玻璃合金由鎂、銅和稀土族元素釓制造而成，擁有很多塑料特有的特性;可順應很多形狀、大批量地制造并能作為碳納米管的有效基體。除了具有優(yōu)良的導電性之外，這種金屬玻璃也擁有非常穩(wěn)定的熱性能，這意味著即使經(jīng)受高溫，它也能保持其形狀和耐用性?？茖W家們表示，以上諸多優(yōu)勢和其卓越的電子發(fā)射屬性，使得這種非晶塊金屬玻璃成為制造電子發(fā)射設備的最好材料之一。

(中國科技網(wǎng))

Study on Distribution of Heavy Metal in Animal Products by Multivate Statiscal Technique

He Guihua， Chen Xi
（Jinan Entry-Exit Inspection and Quarantine Burea， Jinan 250014，China）
Sun Zhicheng Zhang Tongyu
(Shandong Shiyan High Middle School， Jinan 250010， China） （Qingdao No.2 Middle School, Qingdao 266071， China)

Distribution of lead，cadium，mercury and arsenic in different animal products was studied by multivariate statistical technique of frequency histogram，parameter distribution test of hypotheses and non-parameter density estimate. Function curves of density estimate and empirical distribution of four kinds of elements were given，which explained the distribution of elements. The regular of heavy metal cumulation and enrichment was grasped.

multivariate statistical analysis; heavy metal; animal product; distribution characteristic

O212.6

A

1008-6145(2012)01-0070-03

聯(lián)系人：何桂華；E-mail:heguihuas@126.com

2011- 10-12

10.3969/j.issn.1008-6145.2012.01.022