江蘇省稻米鎘含量調(diào)查及其膳食暴露評估

宋　雯，李優(yōu)琴，呂　康，劉賢金*（1.浙江省農(nóng)業(yè)科學(xué)院農(nóng)產(chǎn)品質(zhì)量標(biāo)準(zhǔn)研究所/浙江省食品安全重點實驗室，杭州31001；.江蘇省農(nóng)業(yè)科學(xué)院食品質(zhì)量安全與檢測研究所/農(nóng)業(yè)部農(nóng)產(chǎn)品質(zhì)量安全風(fēng)險評估實驗室（南京），南京10014）

?

江蘇省稻米鎘含量調(diào)查及其膳食暴露評估

宋雯1，2，李優(yōu)琴2，呂康2，劉賢金2*
（1.浙江省農(nóng)業(yè)科學(xué)院農(nóng)產(chǎn)品質(zhì)量標(biāo)準(zhǔn)研究所/浙江省食品安全重點實驗室，杭州310021；2.江蘇省農(nóng)業(yè)科學(xué)院食品質(zhì)量安全與檢測研究所/農(nóng)業(yè)部農(nóng)產(chǎn)品質(zhì)量安全風(fēng)險評估實驗室（南京），南京210014）

摘要：為評價江蘇省稻米重金屬鎘污染現(xiàn)狀，于2013年在全省13市調(diào)研抽檢了1101份稻谷，采用石墨爐原子吸收光譜法測定稻米中的鎘含量。檢測結(jié)果顯示，全省稻米鎘的檢出率為97.6%，平均含量為0.036 4 mg·kg-1，檢出范圍是0.000 3～0.431 0 mg·kg-1，僅有0.27%的樣本超出0.2 mg·kg-1的限量。整合我國居民20個性別年齡組人群的稻米消費和體重信息，采用非參數(shù)概率模型對稻米鎘的膳食暴露量進(jìn)行評估，并與JECFA推薦的鎘的暫定每月耐受攝入量（PTMI）25 μg·kg-1體重相比，評估結(jié)果表明，我國居民食用江蘇地區(qū)稻米產(chǎn)生的平均鎘暴露風(fēng)險尚可接受，但在P95的風(fēng)險水平下，14歲以下人群的攝入量相對較高，占PTMI的86.3%～130.7%，其中11歲以下人群中有6.8%～9.0%的個體攝入量超過PTMI，潛在風(fēng)險較大。建議對稻米中的重金屬鎘進(jìn)行持續(xù)監(jiān)測。

關(guān)鍵詞：鎘；稻米；膳食暴露；概率評估；江蘇

宋雯，李優(yōu)琴，呂康，等.江蘇省稻米鎘含量調(diào)查及其膳食暴露評估［J］.農(nóng)業(yè)環(huán)境科學(xué)學(xué)報, 2016, 35（5）：886-893.

SONG Wen, LI You-qin, Lü Kang, et al. Survey and dietary exposure assessment of cadmium in milled rice in Jiangsu Province［J］. Journal of Agro-Environment Science, 2016, 35（5）:886-893.

鎘是一種天然穩(wěn)定地存在于地殼中，含量甚微的重金屬元素，但其在電鍍、化工、電子和核工業(yè)等領(lǐng)域得到廣泛應(yīng)用，近年來在環(huán)境中急劇釋放［1］。鎘對動物和植物都是非必需的，且蓄積在人體內(nèi)的半衰期很長［2］，慢性鎘中毒會損傷腎臟和骨骼［3-9］?？紤]到鎘長期暴露對人體健康的危害而更需要引起重視，JECFA （Joint FAO/WHO Expert Committee on Food Additives，聯(lián)合國糧農(nóng)組織和世界衛(wèi)生組織下的食品添加劑聯(lián)合專家委員會）在2010年將鎘的健康指導(dǎo)值由暫定每周耐受攝入量（Provisional Tolerable Weekly Intake，PTWI）7 μg·kg-1體重（相當(dāng)于每天1 μg·kg-1體重）改為暫定每月耐受攝入量（Provisional Tolerable Monthly Intake，PTMI），并調(diào)低數(shù)值至25 μg·kg-1體重（相當(dāng)于每天0.8 μg·kg-1體重）［2］。

膳食是人體攝入鎘的最主要途徑［10］。大米是我國居民膳食鎘的主要來源［11-13］，稻米鎘污染現(xiàn)已成為研究熱點問題［9，14-15］。以往我國對鎘攝入量的研究大都基于總膳食調(diào)查開展［12，16-18］。高俊全等［12］以成年男子為例，比對了我國3次總膳食研究的結(jié)果（1990、1992、 2000年），得出了鎘的膳食攝入量（全國平均值）呈上升趨勢的結(jié)論。但基于總膳食研究的評估往往對單種農(nóng)產(chǎn)品的針對性不足，關(guān)注的人群也較為單一，且立足于點評估方法，忽視了評估結(jié)果的變異性和不確定性。而另一方面，針對稻米鎘的人群食用風(fēng)險評價往往以市售大米為研究對象［19-20］，無法追溯產(chǎn)地，缺乏對源頭安全的關(guān)注。

江蘇地處長三角，是我國水稻播種面積和產(chǎn)量較大的省份之一?；诋?dāng)前經(jīng)濟發(fā)展與環(huán)境保護(hù)的現(xiàn)實矛盾，在保證糧食數(shù)量安全的同時還必須兼顧人群健康，開展針對原產(chǎn)地的風(fēng)險監(jiān)測與評價是必要手段。為此，本文特以江蘇稻米為研究對象，結(jié)合我國居民稻米消費和體重信息，采用非參數(shù)概率模型對20個性別年齡組人群的稻米鎘膳食暴露量進(jìn)行評估，并將評估結(jié)果與以往相關(guān)研究得到的結(jié)果進(jìn)行比較，以期為相關(guān)質(zhì)量標(biāo)準(zhǔn)和公共衛(wèi)生政策的制定提供有益的參考。

1　材料與方法

1.1樣品采集

于2013年的收獲季節(jié)，在江蘇蘇州、無錫、常州、南京、鎮(zhèn)江、揚州、泰州、南通、淮安、宿遷、鹽城、徐州和連云港，通過田間現(xiàn)場抽樣的方式，分別抽樣86、90、95、83、70、60、97、80、80、90、90、90、90份，共計1101份稻谷樣品。調(diào)查所選大部分縣市的水稻產(chǎn)量均在10萬t以上，各抽檢田塊的種植面積都在3.33 hm2以上。采樣時，在每個田塊的中心位置采用5點取樣法，每點取3株，共150株混合成一份稻谷樣品。在核對確認(rèn)樣品信息無誤、包裝完好后，分為待檢樣和復(fù)檢樣入庫待化學(xué)分析。

1.2樣品的前處理與分析測定

稻谷經(jīng)出糙、磨精之后，得到精米。再以小型粉碎機將精米粉碎，過40目篩后貯存于干凈塑料袋或塑料瓶中備用。

參考NY/T1100—2006和GB/T5009.15—2003中的濕法消解，稱取0.5 g粉碎均勻的樣品，與8 mL硝酸溶液混合，放置過夜，次日于消解爐上緩慢加熱消化，劇烈反應(yīng)過后，稍冷加入2 mL過氧化氫溶液，繼續(xù)消化至消解液體積為1～2 mL（呈無色透明或略帶黃色），取下冷卻，用純水轉(zhuǎn)移并定容至25 mL，上機待測。檢測儀器使用德國耶拿分析儀器股份公司生產(chǎn)的石墨爐原子吸收光譜儀novAA400。

本研究所建分析方法的回收率為83%～112%，檢出限（Limit of detection，LOD）為0.000 1 mg·kg-1，定量限（Limit of quantification，LOQ）為0.000 3 mg·kg-1。

1.3質(zhì)量控制

檢測過程中使用的酸和過氧化氫為優(yōu)級純，其他試劑為分析級。測試用水符合GB/T 6682—2008二級用水標(biāo)準(zhǔn)。所有器皿在硝酸溶液中浸泡24 h后用去離子水沖洗。每批樣品分析都做試劑空白和大米粉標(biāo)準(zhǔn)樣品（GBW080684）測定，以保證分析結(jié)果的可靠性。

1.4未檢出值的處理

由于檢測方法和儀器的局限，以及痕量數(shù)據(jù)的客觀存在，檢測結(jié)果中常含有一定比例的未檢出值。這些未檢出值可能是真實的“0”，也可能是低于LOD的痕量值。為盡可能有效地解讀利用這部分信息，本研究采用WHO推薦的替代法對未檢出值進(jìn)行處理，即選用較為通行的以1/2LOD替換未檢出值［21］。

1.5各城市樣本的一致性分析

因各城市樣本的總體都呈右偏分布（非正態(tài)），故采用非參數(shù)的Kruskal-Wallis檢驗，對13個獨立樣本來自的13個總體的分布作一致性分析，以檢驗多個總體的分布是否存在顯著差異。上述分析使用SPSS 22.0完成。

1.6稻米中鎘的膳食暴露評估

稻米中鎘的膳食暴露評估，首先需要整合稻米中鎘含量數(shù)據(jù)（mg·kg-1）和目標(biāo)人群每日稻米消費量數(shù)據(jù)（g·kg-1體重），得到膳食暴露量的估計值，再將結(jié)果與相應(yīng)的健康指導(dǎo)值（PTMI）進(jìn)行比較，權(quán)衡風(fēng)險［22］。

每日膳食暴露估計值可通過下式計算：

式中：y代表個體每日由稻米攝入鎘的量，μg·kg-1體重；x代表個體每日的稻米消費量，g·d-1；c代表稻米中鎘的含量，mg·kg-1；w代表個體的體重，kg。

評估所需的居民稻米消費量調(diào)查數(shù)據(jù)和體重數(shù)據(jù)分別引自《中國居民營養(yǎng)與健康狀況調(diào)查報告之二——2002膳食與營養(yǎng)素攝入狀況》［23］和《中國居民營養(yǎng)與健康狀況調(diào)查報告之三——2002居民體質(zhì)與營養(yǎng)狀況》［24］公布的第4次總膳食調(diào)查數(shù)據(jù)。本研究將這兩套數(shù)據(jù)進(jìn)行整合，把目標(biāo)人群分為20個性別年齡組［25］，見表1。

由于JECFA推薦使用PTMI來衡量人體鎘的攝入安全［2］，需要對公式（1）稍作調(diào)整，將其右側(cè)部分乘以系數(shù)30，以此計算個體的月暴露量。

表1　目標(biāo)人群體重及其稻米消費量Table 1 Body weights and daily rice consumption of different subpopulations

1.7評估模型構(gòu)建

本研究采用概率模型作為評估方法［26-28］，將稻米鎘含量和人群稻米消費量視為兩個獨立分布的總體。在獲得兩總體的分布特征和總體參數(shù)后，基于描述暴露過程的公式（2），利用計算機模擬對兩個總體進(jìn)行bootstrap抽樣，再對每組bootstrap樣本作100 000次的Monte Carlo隨機抽樣，用每一輸入分布的隨機數(shù)值來模擬計算膳食暴露。通過2000次的bootstrap迭代重復(fù)上述過程，即可得到人群暴露量的概率分布。根據(jù)該分布可以算得一系列的統(tǒng)計量，如平均值與P50、P90和P95等百分位數(shù)可作為目標(biāo)人群攝入量的估計值，以及相應(yīng)攝入量估計值的置信區(qū)間，從而實現(xiàn)對變異性和不確定性的定量描述。上述模擬過程在SAS/IML環(huán)境下完成。

模型總體分布的構(gòu)建有非參數(shù)和參數(shù)兩種方法［22］。非參數(shù)方法是將目標(biāo)參數(shù)觀察值的數(shù)據(jù)集直接定義為離散的均勻分布作為該參數(shù)的分布［22，29］，以保證每一數(shù)據(jù)在有放回抽樣時有相同的概率被抽到。參數(shù)方法則是通過對目標(biāo)參數(shù)觀察值的分布擬合，將其定義為某一特殊分布［22］。本研究中，稻米鎘含量的數(shù)據(jù)源于多個獨立樣本，宜在對多樣本一致性分析的基礎(chǔ)上，選擇適宜的分布構(gòu)建方法。

2　結(jié)果與分析

2.1檢測結(jié)果與分析

1101份稻米樣品鎘含量分布見圖1，檢測結(jié)果的描述性統(tǒng)計見表2。結(jié)果顯示：稻米中鎘的檢出率較高，累計達(dá)97.6%，檢出范圍為0.0003～0.4310mg·kg-1。比照GB 2762—2012對于稻米中鎘的限量（Maximum limit，ML）為0.2 mg·kg-1，僅有3份樣品超出限量，累計超標(biāo)率為0.27%。未檢出值經(jīng)1/2LOD替代后，所抽稻米樣品的鎘含量平均值為0.036 4 mg·kg-1，中位數(shù)為0.022 3 mg·kg-1。

圖1　稻米樣品中鎘含量分布Figure 1 Frequency distribution of cadmium concentrations in rice

Kruskal-Wallis檢驗顯示，檢驗統(tǒng)計量的漸進(jìn)顯著性取值遠(yuǎn)小于0.000 1，故可以認(rèn)為13個城市稻米的鎘含量并不都服從具有相同分布的總體，分布差異具有統(tǒng)計學(xué)意義。經(jīng)多重比較得到相應(yīng)齊性子集的劃分，見表3。需要注意的是，該檢驗及相應(yīng)多重比較結(jié)果針對的是各城市稻米鎘含量的數(shù)學(xué)分布，而非基于數(shù)值高低對質(zhì)量優(yōu)劣的判定。

表2　1101份稻米樣品中鎘含量的分析結(jié)果（2013年）Table 2 Statistical summary of cadmium concentrations in 1101 milled rice samples（in 2013）

表3　基于Kruskal-Wallis檢驗的齊性子集劃分Table 3 Homogeneous subsets based on Kruskal-Wallis test

2.2江蘇稻米中鎘的膳食暴露評估

基于Kruskal-Wallis檢驗結(jié)果，本研究選用非參數(shù)方法構(gòu)建膳食暴露評估模型。評估結(jié)果采用變異性和不確定性描述，即以20個性別年齡組人群基于各自攝入分布下的均數(shù)、P50、P75、P90和P95的估計值描述變異性，以其90%的置信區(qū)間描述不確定性，詳見表4。比照J(rèn)ECFA推薦的鎘的PTMI（25μg·kg-1體重），食用江蘇稻米產(chǎn)生的平均鎘暴露風(fēng)險尚可接受，但在低年齡組以及高端暴露人群中均存在健康隱患。

表4顯示，20個性別年齡組人群的平均稻米鎘攝入量介于4.225～9.339 μg·kg-1體重，未超過PTMI（占比16.9%～37.4%），總體安全。但在暴露量分布的P95下，14歲以下人群（2～14歲）呈現(xiàn)高暴露，攝入量達(dá)21.579～32.684 μg·kg-1體重，占PTMI的86.3%～130.7%，不容樂觀。14歲以上人群面臨的風(fēng)險則相對較小，P95風(fēng)險水平的暴露量介于14.786～20.103 μg· kg-1體重，占PTMI的59.1%～80.4%。

表4　目標(biāo)人群的稻米鎘月暴露量的估計（μg·kg-1體重）Table 4 Estimated monthly cadmium exposures of different subpopulations via milled rice（μg·kg-1體重）

20個性別年齡組人群的稻米鎘攝入量分布見表5。在各組人群中，相應(yīng)有0.2%～9.0%的個體的稻米鎘攝入量大于PTMI，其中11歲以下的人群中，有6.8%～9.0%的個體攝入量超過PTMI，應(yīng)該引起重視。

2.3與總膳食研究的比較

2007年全國總膳食研究的結(jié)果表明，成年男子膳食鎘的平均攝入量為每人每周5.27 μg·kg-1體重［13］，即為每人每月22.586μg·kg-1體重。因居民膳食鎘的來源主要是谷類，尤其是大米，貢獻(xiàn)率近50%［11，13］，由此可估算出2007年我國成年男子因消費米制品而產(chǎn)生的鎘平均暴露量約為11.293 μg·kg-1體重。

在本研究中，18歲以上的成年男子因食用江蘇稻米而產(chǎn)生的每月鎘暴露量的平均值為4.503 μg·kg-1體重。在所使用的人群稻米消費和體重信息大體一致的前提下，與總膳食研究結(jié)果相對照可知，本研究結(jié)果一方面佐證了大米是膳食鎘的主要來源之一，另一方面也從側(cè)面反映出江蘇省稻米的鎘含量要低于全國平均水平。

表5　20個性別年齡組人群的稻米鎘攝入量分布Table 5 Distributions of estimated cadmium exposure for populations of interest via milled rice consumption

3　討論

水稻是我國重要的糧食作物，全國約有65%以上的人口以稻米為主糧［30］，稻米的安全性直接關(guān)系到我國居民的健康。因此，立足產(chǎn)地開展稻米鎘的膳食暴露評估，提出稻米的安全生產(chǎn)及消費對策，有其必要性和緊迫性。本研究對江蘇全省各市2013年當(dāng)季收獲的稻米進(jìn)行了調(diào)查，整合全國居民的稻米消費和體重信息，采用非參數(shù)概率模型對20個性別年齡組人群的稻米鎘的膳食暴露量及分布情況做了全面的評估。研究結(jié)果顯示，所檢稻米樣品中僅有部分受到鎘污染，但仍需關(guān)注低年齡組以及高暴露人群的稻米鎘攝入風(fēng)險。

本研究1101份稻米樣品中僅有0.27%的樣品鎘含量超出我國制定的限量標(biāo)準(zhǔn)，說明江蘇稻米的鎘含量總體處于安全水平，但仍存在安全隱患。本研究的調(diào)查結(jié)果顯示，超標(biāo)樣品在空間上呈點狀散落分布于省內(nèi)（常州、揚州和宿遷），與礦區(qū)污染斑塊狀或條帶狀的分布格局不同，可能是由于省內(nèi)部分稻區(qū)的環(huán)境已受工業(yè)點源污染影響。據(jù)第四次全國總膳食研究指出，我國南方主要稻區(qū)均存在不同程度的鎘環(huán)境污染［13］。針對江蘇水稻產(chǎn)區(qū)環(huán)境，金亮等［31］采集測定了蘇北6縣市的稻田土壤及稻米樣品，分析認(rèn)為供試地區(qū)土壤未發(fā)現(xiàn)嚴(yán)重的重金屬污染，滿足綠色食品和無公害農(nóng)產(chǎn)品生產(chǎn)基地要求；而劉洪蓮等［32］采樣分析了太湖地區(qū)某冶煉廠周圍的稻田土壤和稻米，發(fā)現(xiàn)由于工業(yè)的點源污染，供試農(nóng)田土壤和水稻已經(jīng)受到多種重金屬復(fù)合污染，不適于農(nóng)產(chǎn)品的安全生產(chǎn)。由此可見，控制大米中的鎘含量應(yīng)從源頭做起，聚焦土壤-水稻系統(tǒng)中鎘的遷移與籽粒積累規(guī)律，同時制定稻田鎘的限量標(biāo)準(zhǔn)，合理劃分產(chǎn)區(qū)。

降低稻米鎘攝入風(fēng)險，除了從生產(chǎn)上控制稻米鎘含量，還需從消費層面引導(dǎo)，特別應(yīng)加強對高消費人群和敏感人群的關(guān)注。涂杰峰等［33］對福建省的稻米鎘進(jìn)行了健康風(fēng)險評價，結(jié)果表明，食用福建省稻米，有8.67%的居民月鎘暴露量大于PTMI，健康面臨風(fēng)險。Masayuki Ikeda等［34］回顧了日本居民膳食鎘攝入量近40年的變化，經(jīng)水稻產(chǎn)區(qū)調(diào)整和消費引導(dǎo)，污染區(qū)居民（日本居民的體重按50 kg計算［35］）的鎘攝入量已由過去最高360 μg·kg-1體重下降到60 μg·kg-1體重。相較而言，食用江蘇省稻米，鎘的攝入風(fēng)險較低。依本研究樣本數(shù)據(jù)分析，18歲以上成年男子食用江蘇地區(qū)稻米，在高暴露位點P95下的暴露量估計值為15.759 μg·kg-1體重，占PTMI的63.04%，風(fēng)險可接受。而11歲以下人群中有6.8%～9.0%的個體攝入量超過PTMI，面臨風(fēng)險。同時，在高暴露位點P95下，11～14歲人群和14～18歲男性的暴露量已非常接近PTMI，有潛在風(fēng)險。因此，基于攝入風(fēng)險的考量，有必要維持現(xiàn)行大米中鎘限量衛(wèi)生標(biāo)準(zhǔn)不放松，同時關(guān)注并合理調(diào)整低年齡組及高暴露人群的膳食結(jié)構(gòu)，才能在消費層面盡可能地保證人民健康。

在前人對稻米鎘的暴露風(fēng)險研究中，對風(fēng)險的性別指向存有一定爭議。如Mirei等［36］研究了日本神通川流域不同程度的鎘污染區(qū)內(nèi)，稻米鎘含量與70歲以上人口的性別比率的關(guān)系，認(rèn)為鎘暴露對女性的致死風(fēng)險高于男性。韓雪梅等［37］采用年齡系數(shù)法和消費量擬合法估算人群累積稻米鎘暴露水平，其結(jié)果顯示同年齡組男性的鎘暴露水平要高于女性。而本研究針對稻米鎘的暴露分析和全國總膳食鎘的調(diào)查結(jié)果［13］則都未發(fā)現(xiàn)鎘的暴露風(fēng)險與性別差異有關(guān)聯(lián)。由表4可見，2～4、30～45、60～70歲這3組人群的男性鎘暴露量較同組女性人群略低，其余年齡組的暴露量則呈男高女低趨勢，且各年齡組的暴露量在性別間的差異并不顯著。這一分歧可能是由于各研究針對的樣本和采用的方法均不一樣，結(jié)論的普適性還有待進(jìn)一步針對性的研究論證。

膳食暴露評估是一項對風(fēng)險定性定量的系統(tǒng)工程，不可避免會有不確定性。本研究中不確定性主要源于建模過程和相關(guān)參數(shù)援引［38］。為降低模型中隨機抽樣帶來的不確定性，設(shè)置了100 000次Monte Carlo模擬［39］和2000次的bootstrap迭代［40］。而參數(shù)的不確定性則主要源于最新膳食數(shù)據(jù)的匱乏。本研究中目標(biāo)人群體重及其稻米消費量信息都來源于衛(wèi)生部門2002年公布的數(shù)據(jù)。隨著生活水平的提高，若某些性別年齡組人群的稻米消費量減少或體重增加，則本研究的評估結(jié)果有高估風(fēng)險的可能；反之則有可能低估。因此，更為符合實際情況的評估結(jié)果還有賴于多平臺數(shù)據(jù)的共享和整合利用。

4　結(jié)論

江蘇地區(qū)稻米的鎘污染情況總體較輕，僅有0.27%的樣本超出限量。消費這一區(qū)域的稻米，目標(biāo)人群面臨的平均鎘暴露風(fēng)險處在可接受范圍內(nèi)。但低年齡組和高暴露人群面臨的潛在風(fēng)險隱患應(yīng)引起重視，本研究中14歲以下人群（2～14歲）為敏感性群體，其中11歲以下人群中，有6.8%～9.0%的個體攝入量超過PTMI，面臨的潛在風(fēng)險較大。

參考文獻(xiàn)：

［1］J?rup L. Hazards of heavy metal contamination［J］. British Medical Bulletin, 2003, 68（1）：167-182.

［2］World Health Organization（WHO）. Safety evaluation of certain food additives and contaminants（WHO food additives series; 64）［R］. Geneva：World Health Organization, 2011：372-374.

［3］J?rup L, Berglund M, Elinder C G, et al. Health effects of cadmium exposure：A review of the literature and a risk estimate［J］. Scandinavian Journal of Work, Environment & Health, 1998, 24（Suppl 1）：1-51.

［4］Buchet J P, Lauwerys R, Roels H, et al. Renal effects of cadmium body burden of the general population［J］. The Lancet, 1990, 336（8717）：699-702.

［5］Friberg L. Health hazards in the manufacture of alkaline accumulators withspecialreferencetochroniccadmiumpoisoning［J］.ActamedicaScandinavica, 1950, 240（Suppl）：1-124.

［6］Bernard A, Roels H, Buchet J P, et al. Cadmium and health：The Belgian experience［J］. IARC Scientific Publications, 1992, 118：15-33.

［7］Hellstr?m L, Elinder C G, Dahlberg B, et al. Cadmium exposure and end-stage renal disease［J］. American Journal of Kidney Diseases, 2001, 38（5）：1001-1008.

［8］Staessen J A, Roels H A, Emelianov D, et al. Environmental exposure to cadmium, forearm bone density, and risk of fractures：Prospective population study［J］. The Lancet, 1999, 353（9159）：1140-1144.

［9］Nordberg G, Jin T, Bernard A, et al. Low bone density and renal dysfunction following environmental cadmium exposure in China［J］. Ambio, 2002, 31（6）：478-481.

［10］World Health Organization（WHO）. Cadmium. Environmental health criteria 134［EB/OL］.［2015-10-08］. http：//www. inchem. org/documents/ehc/ehc/ehc134. htm.

［11］中華人民共和國國家衛(wèi)生和計劃生育委員會.《食品中污染物限量》（GB 2762—2012）問答［EB/OL］.［2015-10-08］. http：//www. moh. gov.cn/wsb/pzcjd/201301/58113f0b463e4e1d89769e47b9555c03.shtml. National Health and Family Planning Commission of the People’s Republic of China. Maximum limit of contaminants in foods（GB 2762—2012）Q&A［EB/OL］.［2015-10-08］. http：//www. moh. gov. cn/wsb/ pzcjd/201301/58113f0b463e4e1d89769e47b9555c03. shtml.

［12］高俊全,李筱薇,趙京玲. 2000年中國總膳食研究：膳食鉛、鎘攝入量［J］.衛(wèi)生研究, 2006, 35（6）：750-754. GAO Jun-quan, LI Xiao-wei, ZHAO Jing-ling. 2000 Chinese total diet study：The dietary lead and cadmium intakes［J］. Journal of Hygiene Research, 2006, 35（6）：750-754.

［13］李筱薇.中國總膳食研究應(yīng)用于膳食元素暴露評估［D］.北京：中國疾病預(yù)防控制中心, 2012：99-105. LI Xiao-wei. The application of total diet study to evaluate the dietary intakes of trace elements［D］. Beijing：Chinese Center for Disease Control and Prevention, 2012：99-105.

［14］Shi J, Li L Q, Pan G X. Variation of grain Cd and Zn concentrations of 110 hybrid rice cultivars grown in a low-Cd paddy soil［J］. Journal of Environmental Sciences, 2009, 21（2）：168-172.

［15］Li Z W, Li L Q, Pan G X, et al. Bioavailability of Cd in a soil-rice system in China：Soil type versus genotype effects［J］. Plant and Soil, 2005, 271（1）：165-173.

［16］梁高道,革麗亞,譚慧,等.武漢市城市成人居民總膳食中鎘的健康風(fēng)險評價［J］.中國衛(wèi)生檢驗雜志, 2008, 18（12）：2715-2716. LIANG Gao-dao, GE Li-ya, TAN Hui, et al. Health risk assessment of cadmium via diet for adults in Wuhan［J］. Chinses Journal of Health Laboratory Technology, 2008, 18（12）：2715-2716.

［17］洪華榮,王娟,陳劍鋒,等.廈門市居民重金屬膳食攝入水平評估［J］.現(xiàn)代預(yù)防醫(yī)學(xué), 2015, 42（9）：1580-1583. HONG Hua-rong, WANG Juan, CHEN Jian-feng, et al. Assessment on dietary intake level of heavy metals among residents, Xiamen［J］. Modern Preventive Medicine, 2015, 42（9）：1580-1583.

［18］楊巍,蔣冬梅,王定勇.重慶市居民重金屬膳食攝入水平研究［J］.生態(tài)毒理學(xué)報, 2009, 4（1）：50-56. YANG Wei, JIANG Dong-mei, WANG Ding-yong. Investigation on dietary intake of heavy metals intaken by inhabitants in Chongqing［J］. Asian Journal of Ecotoxicology, 2009, 4（1）：50-56.

［19］Qian Y Z, Chen C, Zhang Q, et al. Concentrations of cadmium, lead, mercury and arsenic in Chinese market milled rice and associated population health risk［J］. Food Control, 2010, 21：1757-1763.

［20］李優(yōu)琴,李榮林,石志琦.市售大米重金屬污染狀況及健康風(fēng)險評價［J］.江蘇農(nóng)業(yè)學(xué)報, 2008, 24（6）：977-978. LI You-qin, LI Rong-lin, SHI Zhi-qi. Pollution and risk evaluation of heavy metals in rice［J］. Jiangsu Journal of Agricultural Sciences, 2008, 24（6）：977-978.

［21］World Health Organization Regional Office for Europe（GEMS/Food-EURO）. Second workshop on reliable evaluation of low-level contamination of food：Report on a workshop in the frame of GEMS/Food-EURO［R/OL］.［2015-10-08］. http：//www. who. int/foodsafety/publications/chem/en/lowlevel_may1995. pdf, 1995.

［22］World Health Organization（WHO）. Principles and methods for the risk assessment of chemicals in food［M］. Geneva：WHO, 2009.

［23］翟鳳英,楊曉光.中國居民營養(yǎng)與健康狀況調(diào)查報告之二：2002膳食與營養(yǎng)素攝入狀況［M］.北京：人民衛(wèi)生出版社, 2006：21-23. ZHAI Feng-ying, YANG Xiao-guang. A survey on the Chinese national health and nutritionⅡ：The national diet and nutrition in 2002［M］. Beijing：People’s Medical Publishing House, 2006：21-23.

［24］楊曉光,翟鳳英.中國居民營養(yǎng)與健康狀況調(diào)查報告之三：2002居民體質(zhì)與營養(yǎng)狀況［M］.北京：人民衛(wèi)生出版社, 2006：51-52. YANG Xiao-guang, ZHAI Feng-ying. A survey on the Chinese national health and nutritionⅢ：The national health and nutrition in 2002 ［M］. Beijing：People’s Medical Publishing House, 2006：51-52.

［25］宋雯,陳志軍,朱智偉,等.南方6省稻米總汞含量調(diào)查及其膳食暴露評估［J］.農(nóng)業(yè)環(huán)境科學(xué)學(xué)報, 2011, 30（5）：817-823. SONG Wen, CHEN Zhi-jun, ZHU Zhi-wei, et al. Survey and dietary exposure assessment of total mercury in milled rice farmed in 6 provinces of Southern China［J］. Journal of Agro-Environment Science, 2011, 30（5）：817-823.

［26］European Food Safety Authority Scientific Committee. Guidance of the scientific committee on a request from EFSA related to uncertainties in dietary exposure assessment［J］. The EFSA Journal, 2006, 438：1-2.

［27］van Asselt E D, Sterrenburg P, Noordam M Y, et al. Overview of available methods for risk based control within the european union［J］. Trends in Food Science & Technology, 2012, 23（1）：51-58.

［28］van Der Voet H, Slob W. Integration of probabilistic exposure assessment and probabilistic hazard characterization［J］. Risk Analysis, 2007, 27（2）：351-371.

［29］Vose D. Risk analysis：A quantitative guide［M］. 3rd Ed. Chichester：John Wiley & Sons Ltd, 2008：263-297.

［30］龔偉群,潘根興.中國水稻生產(chǎn)中Cd吸收及其健康風(fēng)險的有關(guān)問題［J］.科技導(dǎo)報, 2006, 24（5）：43-48. GONG Wei-qun, PAN Gen-xing. Issues of grain Cd uptake and the potential health risk of rice production sector of China［J］. Science and Technology Review, 2006, 24（5）：43-48.

［31］金亮,李戀卿,潘根興,等.蘇北地區(qū)土壤-水稻系統(tǒng)重金屬分布及其食物安全風(fēng)險評價［J］.生態(tài)與農(nóng)村環(huán)境學(xué)報, 2007, 23（1）：33-39. JIN Liang, LI Lian-qing, PAN Gen-xing, et al. Distribution of heavy metals in the soil-rice system and food exposure risk assessment of North Jiangsu, China［J］. Journal of Ecology and Rural Environment, 2007, 23（1）：33-39.

［32］劉洪蓮,李艷慧,李戀卿,等.太湖地區(qū)某地農(nóng)田土壤及農(nóng)產(chǎn)品中重金屬污染及風(fēng)險評價［J］.安全與環(huán)境學(xué)報, 2006, 6（5）：60-63. LIU Hong-lian, LI Yan-hui, LI Lian-qing, et al. Pollution and risk evaluation of heavy metals in soil and agro-products from an area in the Taihu Lake region［J］. Journal of Safety and Environment, 2006, 6（5）：60-63.

［33］涂杰峰,劉蘭英,羅欽,等.福建省稻米鎘含量及其健康風(fēng)險［J］.農(nóng)業(yè)環(huán)境科學(xué)學(xué)報, 2015, 34（4）：695-701. TU Jie-feng, LIU Lan-ying, LUO Qin, et al. Contents and health risk assessment of cadmium in milled rice in Fujian Province［J］. Journal of Agro-Environment Science, 2015, 34（4）：695-701.

［34］Masayuki Ikeda, Takafumi Ezaki, Teruomi Tsukahara, et al. Dietary cadmium intake in polluted and non-polluted areas in Japan in the past and in the present［J］. Int Arch Occup Environ Health, 2004, 77 （4）：227-234.

［35］Takeya Inaba, Etsuko Kobayashi, Yasushi Suwazono, et al. Estimation of cumulative cadmium intake causing Itai-Itai disease［J］. Toxicology Letters, 2005, 159（2）：192-201.

［36］Mirei U, Etsuko K, Yasushi S, et al. Cadmium exposure aggravates mortality more in women than in men［J］. International Journal of Environmental Health Research, 2006, 16（4）：273-279.

［37］韓雪梅,張文麗,尚琪.鎘污染區(qū)人群累積米鎘暴露水平兩種估算方法的比較［J］.衛(wèi)生研究, 2015, 44（6）：939-969. HAN Xue-mei, ZHANG Wen-li, SHANG Qi. Comparison of two methods in calculating the cumulative cadmium intake（CCI）from rice in the cadmium polluted area［J］. Journal of Hygiene Research, 2015, 44（6）：939-969.

［38］US Environmental Protection Agency（US EPA）. Risk assessment guidance for superfund Volume I：Human health evaluation manual（Part A）［R］. Washington DC：Office of Emergency and Remedial Response, 1989：6-47.

［39］Waldo J de B, Voet van der H. MCRA, Release 6 a web-based program for Monte Carlo risk assessment［R/OL］.［2015-10-08］. https：//mcra. rivm. nl/Documentation/Manual6. pdf.

［40］Efron B, Tibshirani R J. An introduction to the bootstrap［M］. NewYork：Chapman & Hall, 1993：14-15, 275.

中圖分類號：X820.3

文獻(xiàn)標(biāo)志碼：A

文章編號：1672-2043（2016）05-0886-08

doi:10.11654/jaes.2016.05.010

收稿日期：2015-11-18

基金項目：國家自然科學(xué)基金項目（31401780）；江蘇省自然科學(xué)基金項目（BK20140746）；國家農(nóng)產(chǎn)品質(zhì)量安全風(fēng)險評估項目（GJFP2015003001）

作者簡介：宋雯（1984—），女，江蘇常州人，博士，助理研究員，主要從事農(nóng)產(chǎn)品質(zhì)量安全風(fēng)險評估研究。E-mail：wengweng101@gmail.com

*通信作者：劉賢金E-mail：jaasliu@jaas.ac.cn

Survey and dietary exposure assessment of cadmium in milled rice in Jiangsu Province

SONG Wen1,2, LI You-qin2, Lü Kang2, LIU Xian-jin2*
（1.Institute of Quality and Standards for Agricultural Products, Zhejiang Academy of Agricultural Sciences/Zhejiang Province Key Lab for Food Safety，Hangzhou 310021, China; 2.Institute of Food Safety and Monitoring Technology, Jiangsu Academy of Agricultural Sciences/ Laboratory of Quality & Safety Risk Assessment for Agro-Products（Nanjing）, Ministry of Agriculture, Nanjing 210014, China）

Abstract：A survey of cadmium levels in milled rice in Jiangsu Province was carried out in 2013. A total of 1101 milled rice samples were collected from 13 cities. Cadmium was measured by graphite furnace atomic absorption spectroscopy. Results showed that 97.6%of the samples contained detectable concentrations of cadmium, which ranged from 0.000 3 to 0.431 0 mg·kg-1, with the mean value of 0.036 4 mg·kg-1, but only 0.27%of the samples had cadmium concentrations exceeding the maximum leve（lML）, 0.2 mg·kg-1. Based on recent rice consumption and body weight, an estimated monthly human intake of cadmium was calculated for 20 age-sex categories of Chinese population using a non-parameter probabilistic model. The provisional tolerable monthly intake（PTMI）of cadmium, as recommended by the Joint FAO/WHO Expert Committee on Food Additives（JECFA）was 25 μg·kg-1bw. The estimated exposure values for populations of interest were compared to the PTMI. This study confirmed a low average health risk of cadmium via the milled rice in Jiangsu Province. For the population below 14 years old, however, exposure to cadmium represented 86.3%～130.7%of the PTMI at P95 risk level. For population below 11 years old, 6.8%～9.0%individuals had cadmium exposure exceeding the PTMI. The present study suggests that it is necessary to establish a program to monitor cadmium dynamics in milled rice continuously.

Keywords：cadmium; milled rice; dietary exposure; probabilistic assessment; Jiangsu Province