-
砷(arsenic)是一种广泛存在于自然界中的有毒有害元素[1]. 人体砷暴露主要来源于受污染的食物和水的摄入. 长期的砷暴露会带来癌症、皮肤病、血管疾病、糖尿病等健康风险[2-4]. 砷的毒性具有形态依赖性,其形态主要包括无机态的砷酸盐(arsenate, As(Ⅴ))、亚砷酸盐(arsenite, As(Ⅲ)),以及有机态的一甲基砷酸盐(monomethylated arsenic, MMA)、二甲基砷酸盐(dimethylated arsenic, DMA)、砷甜菜碱(arsenobetaine, AsB)、砷胆碱(arsenocholine, AsC)、砷糖和砷脂等[5]. 世界卫生组织国际癌症研究机构公布的致癌物清单中,砷和无机砷化合物为一类致癌物, MMA和DMA为潜在的致癌物质;无机砷可以抑制酶的活性,因此比有机砷具有更大的毒性[6];As(Ⅲ)的易迁移性、活性和毒性都远远高于As(Ⅴ),而AsB、AsC、砷糖等有机砷几乎没有毒性[7]. 因此,需结合砷形态来科学地评估食品中砷的健康风险.
寿司(sushi)是一种起源于日本,以海苔卷和米饭为主要原料,结合生鱼片、肉松、芥末、甜虾等配料而制成的吃食[8]. 食品添加剂联合专家委员会(JECFA)第72次会议相关数据显示,海藻、鱼类、贝类、香菇和食用真菌、稻米和米制品及一些肉制品中总砷含量较高[9],是人类膳食砷摄入的主要来源. 海苔作为寿司的主要原料之一,是传统的海产加工品,由新鲜条斑紫菜(Pyropia yezoensis)经清洗、切割、干燥和调味而成. 海藻中砷的含量和形态因种类不同而存在差异,如褐藻(Hijiki)是一种广泛用于亚洲烹饪的可食用海藻,总砷含量高且大多数是无机形式[10]. 此外,水稻由于容易吸收和积累砷,对人类健康构成重大威胁[11]. 鉴于寿司主要食材即是海苔和米饭,可能存在一定的砷污染和人体健康危害问题,然而目前对寿司中砷人体暴露风险的研究仍缺乏.
本实验在南京市场随机选购20份寿司产品,分别测定海苔、米饭以及整个寿司卷中总砷的含量,并采用高效液相色谱-电感耦合等离子体质谱联用技术测定代表性样品中砷的形态,旨在揭示食用寿司的人体健康危害.
-
于2021年8月份在南京市场购买不同品牌、不同实体店的新鲜寿司样品共20份. 所购买的寿司样品主要来自南京市栖霞区、玄武区、秦淮区,分属于15家不同品牌的寿司店,在一定程度上能代表南京市场所售卖寿司的整体状况. 所购买的寿司包含的主要食材有:海苔、米饭、肉松、黄瓜、胡萝卜、大根条、香酥丝等. 考虑到寿司在制作过程中会在米饭中添加少许寿司醋等以调味,本研究中将添加调味剂的寿司米整体作为米饭样品. 根据实验需要,采购时部分样品单独从店家购买了原材料海苔和所用米饭.
-
本实验从市场购买的寿司为新鲜即食食品,在4 h内将所采集的样品带回实验室,编号装入自封袋,从每份寿司中随机选取3小段,分别剥离海苔、寿司米、寿司馅,确保各成分剥离完全,米饭中不沾有海苔和馅料. 各个成分利用万分之一天平称量鲜重后,立即放入–20 ℃冰箱冷冻24 h,再放入真空冷冻干燥机中干燥72 h,利用食物磨碎机磨成粉末.
-
利用USEPA 3050B方法消解样品[12],具体流程为:取寿司样品、米饭样品干重约0.5 g、取海苔样品约0.1 g于50 mL消解管中,首先加入10 mL 50%硝酸溶液(65%优级纯浓硝酸与超纯水体积比1:1 配制),利用石墨炉消解仪在105 ℃下消解,待剩余溶液体积为2—3 mL时,补加5 mL 50%硝酸溶液,继续在105 ℃下消解至近干后,取出消解管,在通风橱内冷却至室温,加入2 mL 30%优级纯过氧化氢溶液,待反应平稳后,继续在105 ℃下消解,直至溶液剩余2 mL左右,取出冷却至室温并用超纯水定容至30 mL. 消解液过0.45 μm滤膜至10 mL塑料离心管中,使用电感耦合等离子体质谱仪(ICP-MS, NexION300X, PerkinElmer, USA)测定消解样中砷的含量. 用1000 mg·L−1的多元素储备液配置标准曲线溶液,以In(0.05 mg·L−1)为内标,确保 ICP-MS 信号的稳定性.
为确保数据的可靠性,当标准曲线的R2>0.9995时,才接受该标准曲线用于实验样品的测定,每个样品,平行消解3 份,3次平行消解测定的砷含量的相对标准偏差小于5%. 对于QA/QC,实验过程中采用大米标准物质GBW10010a(中国地质科学院地球物理地球化学勘查研究所)进行消解和提取,砷含量为0.08±0.01 mg·kg−1,As回收率为86.61%±5.55%. 此外,在每批次消解时,另做3个消解空白,结果证明消解空白中砷含量接近于零,排除人为污染对实验结果造成的影响.
-
部分寿司样品从店家采购到单独的食材(海苔卷和米饭),选取其中的S02、S06、S12、S14样品进行砷形态提取分析(n=3). 将研磨均匀的样品,分别称取1.000 g的米饭和寿司样品粉末,称取0.100 g的海苔样品粉末于50 mL离心管中,准确记录质量,加入20 mL 甲醇/水(3:1, V/V)溶液,55 ℃条件下超声提取1 h;然后在4000 r·min−1下,离心10 min,取上清液到50 mL刻度管中;再次向样品中加入10 mL甲醇/水(1:1, V/V)的溶液,按照第一次的提取条件,重复操作1次,离心后取出上清液,再重复一次后续操作. 合并3次提取的上清液,用超纯水定容到50 mL,过0.22 µm滤膜到进样瓶中,在−80 ℃冰箱保存. 本实验选用的砷形态分离柱为阴离子交换柱(PRP-X100, 250 mm × 4.1 mm, 10 µm, Hamilton, UK),并配有保护柱(Hamilton, UK),以分析样品不同形态砷的质量分数. 流动相由1%的碳酸氢铵和5%的甲醇溶液组成,pH采用浓氨水调至8.72,流速为1.5 mL·min−1. 根据砷形态混合标液分离色谱图出峰情况,该流动相能很好地实现不同砷形态的分离.
-
本研究采用健康风险评价模型评估食用寿司的砷暴露风险,无机砷的人体平均日摄入量(DI)的计算公式为:DI=C×I/W,其中C表示样品中无机砷含量(mg·kg−1),I表示寿司的每日平均摄入量(g),W为人体体重(kg),本研究中成人体重以60 kg计算.
-
本研究中,使用Excel软件进行砷含量、平均值、标准偏差等数值的计算,数据以平均值±标准差表示(
¯x±s ),所有数据采用SigmaPlot 12.5、Origin 2021处理并作图,使用SPSS软件进行相关性分析和t检验分析. -
图1显示了从南京市场采集得到的20种寿司样品中海苔、米饭和其他馅料所占的鲜重相对比例. 可见,从寿司卷上面剥离下来的海苔鲜重仅占整个寿司鲜重的2.20%—5.38%,从寿司卷上面剥离下来的米饭鲜重占整个寿司鲜重的45.64%—87.31%,剩余其他馅料占整个寿司鲜重的8.67%—52.16%,海苔、米饭和其他馅料占寿司鲜重的平均值分别为3.31%±0.88%、69.48%±9.57%、27.20%±10.10%. 显然,在所采集的不同品牌的寿司样品中,米饭均是最主要的成分,而海苔所占比例最小.
-
图2显示了从南京市场采集得到的20种寿司样品、各寿司中的海苔和米饭样品中砷的含量. 由图2可见,不同样品中砷的含量存在着显著差异. 由图2(A)可得,海苔样品中砷含量为3.14—27.18 mg·kg−1,平均值为(9.67±6.65)mg·kg−1. 虽然海水中砷浓度较低且均匀(1—5 µg·L−1)[13],但由于近海海域受到不同程度的重金属污染,加上藻类对包括砷元素在内的微量元素具有特异性吸收作用,导致海藻类制品重金属污染问题受到广泛关注[14]. 已有研究发现紫菜中总砷含量很高,为14.0—42.1 mg·kg−1[15],与本研究中测得的海苔样品总As含量相似. 然而,除褐藻(Hijiki)之外的海藻,JECFA报告显示无机砷占总砷的比例不到15%,无机砷含量通常低于2 mg·kg−1[16]. 欧盟对藻类饲料中总砷限量为40 mg·kg−1,但对食用藻类中总砷含量限值未作明确规定.
由图2(B)可得米饭样品中砷含量为0.10—0.85 mg·kg−1,平均值为(0.31±0.22)mg·kg−1. 我国食品中污染物限量国家标准GB 2762—2017《食品安全国家标准 食品中污染物限量》规定,米饭中无机砷限量为0.20 mg·kg−1[17]. 由于稻米中砷(尤其是无机砷)的过度积累,米饭成为人体砷摄入的主要来源,大约占我国人群无机砷平均摄入量的60%[18-19]. 相关文献报道称米饭中的总砷含量为0.09—0.33 mg·kg−1,在砷污染地区种植的水稻,其籽粒中砷可达到1.50 mg·kg−1,甚至更高[11]. 本研究测得经过烹饪的熟米中总砷含量有11组样品低于 0.2 mg·kg−1,属于正常范围内. 可能受米饭中添加不同调味剂的影响,部分样品中总砷含量接近甚至超过0.3 mg·kg−1.
对寿司整体进行总砷含量测定,由图2(C)可得,寿司样品中总砷含量为0.34—1.57 mg·kg−1,平均值为(0.63±0.27)mg·kg−1. 显然如图2(D)所示,海苔样品中的总砷含量远高于米饭以及寿司样品中的总砷含量. 本研究所采集的寿司样品其馅料占一定比重,但馅料的鲜重受水分含量较多的黄瓜、胡萝卜、大根条等果蔬影响,对各个寿司样品其馅料进行砷含量测定,馅料样品中总砷含量为0.04—0.12 mg·kg−1,均值为(0.06±0.03)mg·kg−1,相较于海苔和米饭总砷含量要低的多. 由于市售寿司的馅料种类复杂多样,其砷含量受馅料成分影响,故寿司馅料在本研究中不做重点分析.
结合不同组分在寿司样品中的相对重量比以及不同组分砷含量,舍弃掉S06、S07、S14海苔和米饭砷贡献率超过100%的3组样品,剩余17组样品不同组分对寿司总砷含量贡献值如图3所示,海苔、米饭和其他成分对寿司总砷浓度贡献均值分别为41.68%±22.43%、36.96%±21.15%、21.36%±15.58%,即海苔在寿司中相对重量占比虽然最小,但对寿司总砷含量起主要贡献作用.
-
将4种形态As标准储备液稀释成浓度为1、2、5、10、20 µg·L−1的砷化合物混合标准溶液,在确立的实验条件下进行分析,以各种砷形态的峰面积对浓度进行线性回归,得出4种砷化物的质量浓度(x)与峰面积(y)的线性关系. 如表1所示,4种砷形态的线性相关系数均在0.9990以上. 20 µg·L−1砷化合物混合标准溶液色谱图如图4A所示,显然As(Ⅲ)、DMA、MMA、As(Ⅴ)能很好地实现分离,各形态峰到达最大值点的时刻分别为2.82、4.51、9.69、27.23 min,并且各峰尖锐,峰形较为对称.
采用建立的上述方法对S02、S06、S12、S14样品(海苔、米饭、寿司)的砷形态进行测定,结果如表2、3、4所示. 4个海苔样品提取液色谱图均出现3个峰,出峰时间分别为2.16、4.51、24.2 min,其中仅第2个峰出峰时间与砷形态标准物质DMA吻合(图4B),另外2个出峰位置与4种砷标准物质出峰时间均不吻合. 研究表明,对于包括条斑紫菜在内的藻类,海水中的砷酸盐很容易通过磷通道跨膜运输被吸收到植物体内,吸收之后,藻类会通过还原和甲基化作用快速解毒砷酸盐,从而转化为砷糖、DMA、甲基砷化合物等[20]. 以往对紫菜砷形态的分析鉴定出在As(Ⅴ)前出峰的砷形态为砷糖-OH(arsenosugar-OH)[21],基于此,判定本研究在24.2 min左右出峰的砷为砷糖-OH. 此外,以往研究发现,利用HPLC-ICP-MS分离和鉴定砷形态时,2.00 min左右会出现AsC、AsB、As(Ⅲ)[21],鉴于AsC是砷甜菜碱的代谢前体,能被快速吸收并转化为AsB[22],故判定2.16 min出峰的砷为AsB. 基于以上形态判定,发现海苔样品中含有AsB、DMA、砷糖-OH这3种形态砷,未检测出无机砷. 本研究方法As形态提取的回收率在67.0%—86.4%(表2),无机砷总量明显低于藻类调味品无机砷限量国际标准0.5 mg·kg−1[17]. 各形态砷在海苔样品可提取总砷含量中占比如图5(A)所示,砷甜菜碱(AsB)、二甲基砷酸盐(DMA)、砷糖-OH占比分别为8.02%±6.57%、37.40%±12.70%和54.59%±17.31%,这些都是毒性极低甚至无毒的有机砷[22],且已有实验证明,小鼠食用紫菜消化的总砷只有0.12%—0.78%在肌肉/器官中积累,约65%—77%都通过粪便排除,不会存在砷健康风险[23].
米饭样品提取液色谱图均出现4个峰,其中第1、第2和第4个峰分别与标准物质As(Ⅲ)、DMA和As(Ⅴ)的出峰时间吻合,而在14.5 min出现的峰与砷标准物质均不吻合(图4C). 研究发现稻田土壤中除存在无机砷(如无机三价砷和五价砷)与甲基砷(如一甲基砷和二甲基砷)外,还广泛存在着巯基砷化合物,包含无机巯基砷和甲基巯基砷,其中二甲基一巯基砷(DMMTA)具有高毒性,是土壤孔隙水中主要的甲基巯基砷形态,且会在水稻籽粒中积累,对食品安全和人体健康构成威胁[24]. 最新研究发现,在大米中除了As(Ⅲ)、DMA和As(Ⅴ),还广泛存在DMMTA,其浓度为DMA的46.6%—74.5%,代表着米饭总砷含量的8.74%—13.6%. [25]. 因此,判定14.5 min出现的峰为DMMTA. 基于此,大米中砷存在形态有As(Ⅲ)、DMA、DMMTA、As(Ⅴ),其中As(Ⅲ)、DMA含量较高,As形态提取的回收率均大于90%(表3). 各形态砷在米饭样品可提取总砷含量占比如图5(B)所示,低毒性的DMA以及高毒性的As(Ⅲ)、DMMTA、As(V)占比分别为22.44%±5.79%、49.15%±9.52%、12.48%±2.36%、11.59±9.65%. 米饭样品中无机砷含量在43.31—88.87 µg·kg−1,均低于米饭无机砷限量国际标准0.20 mg·kg−1.
整个寿司样品包含了海苔和米饭的特征砷形态,具体来说,寿司中含有AsB、DMA、砷糖-OH、As(Ⅲ)、As(Ⅴ)等6种形态砷(图4D),As形态提取的回收率在47.0%—88.0%(表4),各形态砷在寿司样品可提取总砷含量占比如图5(C)所示,AsB、DMA、砷糖-OH、As(Ⅲ)、As(Ⅴ)这5种砷形态占比分别为10.32%±7.91%、34.40%±22.04%、32.87%±22.87%、15.94%±9.19%、2.44%±4.87%,DMA和砷糖是寿司样品最主要的砷形态. 由于继承了海苔和米饭中的砷形态特征,虽然海苔中砷总量特别高,但是寿司中可能产生毒害作用的无机砷,均主要继承自米饭,占寿司砷含量较大部分的有机砷均继承自海苔样品,且不同样品间各种形态有机砷含量相对高低与海苔样品保持一致.
-
目前,已有大量研究针对不同砷富集食品构成的潜在健康风险进行评估,如分析食用药用真菌(39科82属164个物种)中总砷和不同形态砷含量[26];分析不同种类蔬菜的As污染特征[27];测定中国市场生熟鸡肉中砷的浓度,并进行砷形态分布研究等[28];旨在评估人体从该类食品中摄入无机砷相关的癌症风险. 对本研究中S02、S06、S12、S14样品进行食用寿司的相关癌症风险评估,结合砷形态化合物检测回收率数据,4组样品中无机砷含量分别取各自最高值0.10、0.17、0.08、0.07 mg·kg−1,以每日成人摄入一份寿司为准,南京市寿司一份大约为300 g左右鲜重,根据实验室称量的寿司含水率数据,寿司每日平均摄入干重取120 g,估算得到成人食用该4组寿司导致的每天无机砷摄入量为0.20、0.34、0.16、0.14 µg·kg−1·d−1 bw, 平均为(0.21±0.09) µg·kg−1·d−1. 世界卫生组织规定的肺癌发病率增加0.5%的基准剂量(BMDL0.5)为3 µg·kg−1·d−1 bw [29]. 基于此安全限制,食用寿司导致的无机砷摄入量远低于BMDL0.5,表明食用寿司带来砷摄入的健康风险较低. 市场销售寿司虽然由于原材料海苔导致砷总量较高,但无机砷含量较低,因此健康风险较低.
-
(1)南京市场所采集的不同品牌的寿司样品中,米饭均是寿司的最主要成分,达到整个寿司的近70%,而海苔所占重量比极低,不超过5%,这符合市场大部分普通寿司的相对重量比情况.
(2)南京市场销售的寿司不同样品中砷的含量存在着差异,总体趋势是海苔总砷含量近乎是米饭总砷含量的30倍,导致海苔和米饭对寿司总砷含量的贡献率不同,海苔在寿司中相对重量占比虽然非常小,但对寿司总砷含量起主要贡献作用.
(3)市售寿司海苔中含有AsB、DMA、砷糖-OH共3种形态砷,未检测出无机砷;米饭中含有As(Ⅲ)、DMA、DMMTA、As(Ⅴ)等5种形态砷,其中As(Ⅲ)、DMA含量较高;继承自海苔和米饭的砷形态,寿司中相应地主要含AsB、DMA、砷糖-OH、As(Ⅲ)、As(Ⅴ)等形态砷,且有毒的砷形态主要来自于米饭,近乎无毒的有机砷主要来自海苔.
(4)采用砷摄入量评估模型评估人体食用寿司导致的无机砷摄入量显示,食用寿司带来的健康风险较低,市面上销售的寿司可以安全食用.
市售寿司中砷的人体健康风险
Human health risk of arsenic in commercial sushi
-
摘要: 砷是一种在自然界中广泛存在的有毒有害污染物,摄入受污染的食物是人体砷暴露的重要途径. 食品添加剂联合专家委员会(Joint FAO/WHO Expert Committee on Food Additives,JECFA)第72次会议指出,海苔、米饭和米制品是总砷含量较高的食物. 寿司是一种日本传统美食,近年来颇受中国大众喜爱,其主要食材即是海苔和米饭,可能存在一定的砷污染和人体健康危害问题. 本研究从我国南京市20家店面采集了共20份寿司样品,测定样品中总砷含量,并对部分样品进行砷形态分析,进而评估市售寿司中砷的人体健康风险. 对寿司进行食物成分分类,发现海苔、米饭和其他馅料的鲜重在整个寿司中的重量占比平均值分别为3.31%±0.88%、69.48%±9.57%、27.20%±10.10%;海苔、米饭和整个寿司样品(n = 20)中总砷的含量分别为3.14—27.18、0.10—0.85、0.34—1.57 mg·kg−1,平均值分别为(9.67±6.65)、(0.31±0.22)、(0.63±0.27)mg·kg−1. 对4份寿司及对应成分海苔和米饭进行砷形态分析,发现海苔中砷主要为低毒的砷甜菜碱(arsenobetaine,AsB)、二甲基砷酸盐(dimethylated arsenic,DMA)、砷糖(arsenosugar),占比分别为8.02%±6.57%、37.40%±12.70%和54.59%±17.31%,未检测出高毒性的无机砷;米饭中除了含有低毒性的DMA(22.44%±5.79%)外,含有高毒性的亚砷酸盐(arsenite,As(Ⅲ),49.15%±9.52%)、二甲基一巯基砷(dimethylated monothioarsenate,DMMTA,12.48%±2.36%)和砷酸盐(arsenate,As(Ⅴ),11.59%±9.65%);整个寿司样品包含了海苔和米饭的特征砷形态,包括AsB(10.32%±7.91%)、DMA(34.40%±22.04%)、砷糖(32.87%±22.87%)、As(Ⅲ)(15.94%±9.19%)、As(Ⅴ)(2.44%±4.87%)等. 基于寿司中砷含量和形态,成人每天通过食用寿司导致的无机砷摄入量平均为( 0.21±0.09) μg·kg−1·d−1 bw,远低于世界卫生组织规定的BMDL0.5的基准剂量3 μg·kg−1·d−1 bw,表明食用寿司的癌症风险较低. 本研究对南京市场销售的寿司样品进行砷含量、砷形态测定以及开展健康风险评价,对理解食用寿司导致的砷暴露和健康风险具有重要意义.Abstract: Arsenic (As) is a toxic element that exists widely in the environment, with food consumption being an important route of human As exposure. According to the Joint FAO/WHO Expert Committee on Food Additives (JECFA) No. 72, nori, rice, and rice products tend to contain high As concentrations. Sushi is a traditional Japanese food, which has gained popularity in the Chinese population in recent years. With nori and rice being the main ingredients, sushi may contain a considerable amount of As and pose health risk to human health. In this study, a total of 20 sushi samples were collected from 20 stores in Nanjing, China. Total As concentration and species in samples (n = 20 and 4) were measured to assess the health risk associated with sushi consumption. Dividing sushi into different ingredients, we observed that nori, rice, and other fillings contributed to 3.31%±0.88%, 69.48%±9.57%, and 27.20%±10.10% of fresh weight of sushi, respectively. Arsenic concentrations in the nori, rice, and the whole sushi samples were 3.14—27.18, 0.10—0.85, 0.34—1.57 mg·kg−1, with average of (9.67±6.65), (0.31±0.22), (0.63±0.27) mg·kg−1, respectively. Speciation analysis showed that As in the nori mainly occurred as low toxicity arsenobetaine (AsB), dimethylated arsenic (DMA), and arsenosugar, which contributed to 8.02%±6.57%, 37.40%±12.70%, and 54.59%±17.31% of total As, with inorganic As being not detected. In addition to DMA (22.44%±5.79%), the rice contained high toxicity arsenite (As(Ⅲ), 49.15%±9.52%), dimethylated monothioarsenate (DMMTA, 12.48%±2.36%), and arsenate (As(Ⅴ), 11.59%±9.65%). Consisting of nori and rice, sushi contained all As species observed in nori and rice, including AsB (10.32%±7.91%), DMA (34.40%±22.04%), arsenosugar (32.87%±22.87%), As(Ⅲ) (15.94% ± 9.19%), and As(Ⅴ) (2.44%±4.87%). Based on As concentration and species, health risk associated with sushi consumption was assessed. The estimated daily intake of inorganic As for adults was averagely (0.21±0.09) µg·kg−1·d−1 bw, which was far lower than the BMDL0.5 baseline dose of 3 µg·kg−1·d−1 bw proposed by the World Health Organization, indicating that eating sushi had a low cancer risk. In this study, As concentration, As species, and health risk of sushi samples collected from Nanjing market were studied, advancing our understanding of As exposure and health risk associated with sushi consumption.
-
Key words:
- sushi /
- nori /
- rice /
- arsenic /
- health risk assessment.
-
挥发性有机物(volatile organic compounds,VOCs)是城市大气光化学氧化剂和有机气溶胶的重要前体物,对臭氧和细颗粒物的生成起到重要作用[1-3]。工业源VOCs排放占全国人为源排放总量的56.8%,主要来自石化、化工、工业涂装、包装印刷等行业[4]。根据北京、上海和天津3个城市的VOCs源解析结果,汽车尾气和汽油挥发对VOCs的贡献率分别为(40.2%±16.9%)和(19.7%±7.5%)[5-6]。为了加强油气污染的控制,生态环境部先后颁布了《排污许可证申请与核发技术规范 储油库、加油站》(HJ 1118-2020)[7]和《加油站大气污染物排放标准》(GB 20952-2020)[8],对加油和卸油过程的油气回收、油气处理装置和在线监测系统等提出明确要求。
加油站VOCs排放如图1所示包括5个环节:卸油排放、加油排放、呼吸排放、加油枪滴油和胶管渗透排放[9]。其中呼吸排放是加油站埋地油罐内油气压力在达到排放管压力/真空阀(pressure/vacuum valve,P/V阀)的开启压力后排放的VOCs。卸油排放和加油排放分别通过一次和二次油气回收技术控制,埋地油罐呼吸排放通过安装油气处理装置控制,油枪滴油和胶管渗透分别采取防滴油加油枪和低渗透胶管等措施进行控制。目前,国内加油站已实施一次、二次油气回收技术。然而,针对呼吸排放是否需要安装油气处理装置仍存在较大争议,因此开展加油站呼吸排放因子检测可为油气排放控制提供依据。
针对加油站不同工艺过程的排放因子,美国环保署(Environmental Protection Agency,EPA)、加州空气资源委员会(California Air Resources Board,CARB)和欧洲环境署(European Environment Agency,EEA)等机构开展了相关研究[10-12]。我国研究者也在北京、上海等地开展了加油站排放因子的研究工作[13-16]。黄玉虎[17]对比了国内外未控制排放因子(Uncontrolled Emission Factor,UEF),即未采取油气回收措施时,每加注1 L汽油所排放的VOCs质量,中国、美国和欧洲加油站的总UEF分别为3 743 mg·L−1、2 902 mg·L−1、1 787 mg·L−1。
国内加油站排放因子的测算是基于典型加油过程现场测试数据,并类比国外的加油排放因子,通过统计学的方法获得。呼吸排放因子与油品特性、气液比、P/V阀设定值、日均加油量、环境温度、埋地油罐内油品温度等因素相关。如我国颁布的国家标准GB 17930-2016《车用汽油国家标准》[18]中规定,每年11月1日—4月30日,国V油品雷氏蒸气压为45~85 kPa;每年5月1日—10月31日,油品雷氏蒸气压为40~65 kPa;美国EPA在AP-42中计算机动车加油过程排放因子时雷氏蒸气压设定为78.6 kPa。我国普遍采用真空辅助式二次油气回收系统,气液比在1.0~1.2,高于国外接近于1.0的设定值[17,19]。GB 20952-2020要求P/V阀正压开启压力为2.2~3.0 kPa[8],高于国外的设定值0.75~1.5 kPa [9]。本研究采用实验测试和数据统计拟合相结合的方法,测定加油站呼吸排放因子随加油量变化的规律,可为不同规模加油站分级管理提供参考。
1. 研究方法
为了准确检测加油站呼吸排放因子,本研究借鉴美国CARB相关的测试规程[20],开发了如图2所示专用的排放采集系统。通过在每台加油机内部的油气回收管路上设置涡街流量计,采集加油过程中的加油机脉冲信号和流量信号,获得每次加油作业的加油量、回气量和气液比;在油罐区的排放管上安装罗茨流量计和压力变送器,连续监测油罐压力和呼吸排放累计流量[21-23]。流量计采用Dresser公司专门用于测量油气挥发物的罗茨CRM8C175型,误差值为±0.75%;压力变送器为北京昆仑海岸公司定制JYB-315MG,测量结果的最大允许误差为0.2%。
为确保实验测试结果具有代表性,本研究选取北京某年销汽油量约9 000 t的二级加油站为测试站点。测试站点配有4个30 m³地下汽油埋地油罐,使用电子式调节系统使加油气液比约维持在1.10,且卸油排放和加油排放油气回收设施性能良好;油罐通气管联通后共用1根通气立管和P/V,测试站点埋地油罐呼吸排放未安装油气回收系统。
在实验前,按GB 20952-2020附录A-C要求,检测站点加油油气回收系统的气液比和密闭性、油气回收管线液阻液阻等指标;同时,校准所有监测仪表,并在随后的实验过程中定期开展校准工作。排放监测流量计和压力表的安装位置如图3(a)和(b)所示,控制界面见图3(c)。从2021年7月起,对测试站点的加油量、回气量、油罐压力、呼吸排气流量等开展了近1年的连续监测。
根据气液两相平衡原理,埋地油罐和通气管形成的气相空间内油气浓度与温度、压力等物理量有关,呼吸阀排放的油气为罐内饱和汽油蒸气。因此,为降低采样难度和提高准确性,在通气管距离地面1 m处设置浓度检测采样口 (图3(b)) 。通过真空负压采样法,用TEDLAR气袋采集立管内的油气;气样经预处理后,采用安捷伦7890气相色谱检测气体中非甲烷总烃(non-methane hydrocarbon,NMHC)的浓度。
2. 结果与讨论
2.1 日加油量和呼吸排放量的相关性
《加油站大气污染物排放标准》(GB 20952-2020)中规定,各种加油油气回收系统的气液比均应为1.00~1.20。每次加油过程通过油气回收型加油枪回收油气多于加油的体积,导致埋地油罐油气压力增大,发生呼吸排放。在气液比保持稳定的情况下,站点超量回收的油气体积与日加油量相关。因此,本研究根据站内监控系统的加油数据,以及所采集的罗茨流量计数据,对2021年10月15日—11月10日期间的单日累计加油量与单日累计P/V阀排放油气量进行了统计分析,结果见图4。
使用SPSS软件进行相关性分析,由于样本数据小于50个,选择Shapro-Wilk检验样本数据的正态分布[24],计算可得单日加油量Po=0.003<0.05,单日排放量Pg=0.276>0.05,这表明单日加油量不具备正态分布特征。对经过现场观测发现,在加油站的日常作业中,油气回收检测、计量检定、设备维修维护等偶发性作业均可能影响油罐压力。虽然此类作业根据操作规程会不定期开展,但每次均会关停加油设施,对单日加油量产生较大影响,为53 000~92 000 L。
选用Spearman相关系数来度量单日加油量(x)和单日排放量(y)2个变量之间的相关性。计算得到r值为0.319,P值为0.105>0.05。这说明单日加油量和单日排放量之间没有显著的相关关系。因此,对于实际加油站而言,由于单日作业过程存在较多干扰因素,不具备总结单日加油量和单日排放量规律的条件。
2.2 累计加油量和呼吸排放量的相关性
影响加油站VOCs排放清单的不确定性因素较多,包括活动水平的可靠性及VOCs 排放因子的适用性和准确性。
为了减少偶发性作业对单日排放量的波动影响,本研究进一步对多日累计加油量和多日累计排放量进行分析。图5为2021年10月15日—11月10日期间多日累计加油量和多日呼吸排放量数据的变化图,累计排放量与累计加油量之间呈现明显线性关系。
通过SPSS计算Pearson相关系数r=0.998 8,这说明2个变量显著线性相关。使用最小二乘估计法对图5的数据序列做线性回归,即使每个样本点
与回归线上的对应点(xi,yi) 在垂直方向上的距离偏差总和最小。使用SPSS内置工具对上述数据做线性回归,可得公式(1)。(xi,f(xi)) y=0.0366x+1052.7 (1) 式(1)线性相关系数R2=0.997 5,回归方程对样本数据点的拟合优度很高。回归方程的截距为1 052.7 L,对比样本总累计呼吸排放量65 419 L ,仅为1.6 %;对比样本总累计加油量1.73×106 L,仅为0.061%。加油站不进行加油作业时,埋地油罐没有大呼吸,且小呼吸明显小于大呼吸,不会引起超压排放,需要对公式(1)进行修正。根据《大气挥发性有机物源排放清单编制技术指南(试行)》[25],油品储运源主要包括油品储运和加油站,VOCs排放总量可按公式Ei=Pi×EFi计算,仅与排放因子EFi和活动水平Pi有关,其中活动水平Pi为油品的周转量,公式中没有设置截距。因此,对回归方程进一步优化,重新拟合所得回归方程见公式(2)。
y=0.0379x (2) 式(2)的线性相关系数R2=0.995 3,即在特定时间段内累计排放量与累计加油量具有固定斜率的显著线性关系。
《石化行业VOCs污染源排查工作指南》[26]、US EPA参考AP-42[8]文件中(石油液体运输和销售)公式计算加油过程VOCs排放量均与加油站汽油销售量或者排放源的活动水平呈线性相关,这与本研究得到的结论一致。
2.3 排放斜率的影响因素
为了验证上述拟合方程在累计加油量变量上的通用性,本研究选取检测周期内的若干个时间段,使用最小二乘估计法对累计加油和累计排放数据进行线性回归,回归结果如表1所示。
表 1 不同时间段累计加油量与累计排放量线性回归参数Table 1. Linear regression parameters of cumulative refueling volume and cumulative emissions in different test periods时间段 排放斜率 拟合优度 气液比 日均加油量Q /L 2021-10-15至2021-10-23 0.039 7 0.956 9 1.062 70 718 2021-10-24至2021-11-10 0.039 2 0.978 2 1.059 62 746 2022-01-05至2022-01-12 0.035 4 0.982 8 1.066 68 365 2022-02-09至2022-03-22 0.041 2 0.989 8 1.133 29 896 2022-04-15至2022-04-28 0.045 0 0.922 2 1.082 33 791 2022-04-29至2022-05-11 0.015 5 0.851 5 1.080 21 075 2022-05-12至2022-05-31 0.007 6 0.644 2 1.102 15 918 2022-07-01至2022-07-31 0.038 3 0.953 2 1.061 33 474 该加油站所有时间段的加油气液比均为 (1.08±0.05) ,气液比处于非常稳定的状态。气液比与拟合所得斜率的相关系数r=−0.060 5,说明这2个变量的相关性不显著,故暂不考虑气液比对排放量的影响。
对于指定加油站,加油站内外部环境、油非互动等是影响油品销售量的重要因素。在2021年10月至2022年4月期间,测试站点日均加油量较大,排放斜率稳定在0.035 4~0.045;2022年4月底到6月期间,北京因新型冠状病毒肺炎疫情发布系列社会面防控要求社会面流动和活动的显著减少直接导致了加油站日均加油量的减少,排放斜率相应明显降低。测试站点5月12日至31日,日均加油量降至15 918 L,排放斜率下降至0.007 6。根据现场测试流量和压力数据,2022年5月份仅有9天发生了P/V阀超压开启情况,且油气总排放量仅为6 034 L。
根据EPA在AP-42中油罐车卸油过程排放因子LL和机动车加油过程排放因子LD(1995年)计算公式[9],排放因子与装卸方式、环境温度、油品的雷氏蒸气压等因素有关。2022年4月15日到28日周期内统计所得排放斜率为0.045,高于前4个测试周期的数据。这与环境温度升高,油品的蒸发速率加快,雷氏蒸气压升高有关,因此在以后的研究中,还需进一步细化分析这些参数对排放斜率的影响。同时,由于2022年4月29日以后的3个测试周期内,出现29个样本P/V阀单日排放量为0,这3个周期的排放斜率变化也需从温度、加油量等方面开展进一步的分析。
测试站点不同时间段日均加油量与排放斜率的散点图如图6所示,排放斜率k先随着日均加油量的升高而增大;当日均加油量升至约28 800 L时,k的变化逐渐变小。因此,k是与某个年均销售规模Q加油站排放因子相关的参数。使用最小二乘估计法对日均加油量和排放斜率进行拟合,见式(3)。
k=6.66×10−16×Q3−1.17×10−11×Q2+6.44×10−6×Q−0.069 (3) 式(3)线性相关系数R2=0.935 2。对式(3)求解,日均加油量Q为15 920 L时,排放斜率k=0。上述现象是由于埋地油罐P/V阀的开启压力设定为+(2.2-3.0) kPa,具有一定的保压效果,当日均加油量低于15 920 L时,埋地油罐P/V阀基本没有呼吸排放。这与实际加油站作业中检测的数据基本一致,说明可通过日平均加油量Q(L)分段计算不同规模加油站的排放斜率k,并以此计算排放量,如式(4)所示。
k={6.66×10−16Q3−1.17×10−11Q2+6.44×10−6Q−0.069Q≥159200Q<15920 (4) 2.4 排放斜率与排放因子换算
在实验测试过程中,在加油站通气管采样检测油气的体积浓度,并根据《石化行业VOCs污染源排查工作指南》[26]附表二-19,15.6 ℃时油气分子量为68 g
,得到呼吸排放的非甲烷总烃NHMC平均质量浓度为777 mg·L−1,即呼吸排放因子·mol−1 。UEF=k·777mg·L−1 根据《HJ 1118-2020》中汽油密度为760 g·L−1[18],计算典型规模加油站的排放因子如表2。日加油量44 717 L,即年销售量12 400 t时,呼吸排放因子最高为34.65 mg·L−1;当日加油量增加到72 098 L时,呼吸排放因子与28 839 L的计算结果基本接近。
表 2 不同年汽油销售规模加油站的排放因子Table 2. Emission factors of the gasoline filling station with different annual sales scale年汽油销售规模Qa/t 日均汽油销售量Q/L 呼吸排放因子 / (mg·L-1)UEF 4 416 15 918 0 5 000 18 024 9.50 8 000 28 839 29.45 10 000 36 049 32.83 12 400 44 717 34.65 20 000 72 098 29.07 表1说明,当日加油量Q超过29 000L时,采集足够多样本且在气象因素基本一致情况下,排放斜率k取值为0.037 9~0.041 2;取其平均值0.038 3,也可估算排放因子,
=0.038 3×0.777 g·L−1=29.76 mg·L−1。UEF 对比表2,当日加油量Q为28 839 L时,即年销售量为8 000 t时,排放因子为29.45 mg·L−1。这与上述所得排放因子29.76 mg·L−1基本一致,说明实验测试与拟合公式计算值可以相互验证。
参考黄玉虎等[13,17,19]的研究结果,本研究得出的排放因子为给非ORVR (Onboard Refueling Vapor Recovery,车载加油油气回收系统)车辆加油,且加油站内具备一次和二次油气回收、在线监测系统时,加油站正常运行条件下,埋地油罐呼吸的排放因子约为30 mg·L−1。文献报道的国内外研究机构得出的呼吸排放因子如表3所示。
表 3 国内外不同研究机构得出的排放因子Table 3. Emission factors obtained by different research institutions at home and abroad各研究机构得出的数据为8~122 mg·L−1,与本研究通过实验测试所得的呼吸排放因子相比均不超过同一个数量级,检测结果较为接近。美国EPA、欧洲EEA和上海市环科院是未采安装油气处理装置时的呼吸排放因子,国外的P/V阀正压开启压力设定值普遍低于国内。这说明加装油气处理装置可显著降低呼吸排放因子。美国CARB和北京市环科院呼吸排放因子较本研究得出的结果要低。
2.5 油气处理装置安装的可选择性
《加油站大气污染物排放标准》(GB 20952-2020)实施后,加油站是否需要安装油气处理装置由各省级环保主管部门自行决定,在《排污许可证申请与核发技术规范 储油库、加油站》(HJ 1118-2020)执行过程中,多个省市环保部门均强制要求所有加油站必须安装油气油气处理装置。
与现有文献中基于统计采样分析得出的排放因子相比,本研究进一步明确了加油规模对呼吸排放因子的影响,年汽油销售规模超过4 500 t的加油站才有必要安装油气处理装置对埋地油罐的呼吸排放进行控制。以我国东部某省为例,全省共计有10 000多座加油站,年销汽油量超过4 500 t的油站数量占比仅约10%,有近9 000余座站点并不需要安装油气处理装置,油气处理设施安装后减排效益不明显,而且存在着自身能耗较高、容易产生危废等问题,造成了较大的环保投入浪费,并没有体现“协同增效、减污降碳”的发展思路[27]。按油气处理装置单价10万元计算,全省可减少油气处理装置投资约8亿元;将这些费用投入到其它油气回收环节的技术升级,能起到更好的油气排放控制效果。
3. 结论
1) 加油站的埋地油罐呼吸排放因子与日均加油量密切相关。对于常规配置4台汽油埋地油罐的标准化加油站,加油气液比设定为1.10、P/V阀正压开启压力为2.2~3.0 kPa条件下,日均加油量小于15 920 L时,呼吸排放因子为0;随着日均加油量的升高,呼吸排放因子将逐步升高,超过28 800 L左右后,呼吸排放因子变化趋于稳定,呼吸排放因子约为30 mg·L−1。
2) 综合考虑油气处理装置安装成本和产生的环境效率,年汽油销售规模超过4 500 t的加油站才有必要安装油气处理装置,以降低呼吸排放的油气量。
3) 目前实验加油站的气液比、P/V阀开启压力均按照国家标准要求设定,需要将进一步研究不同气液比设定范围、不同P/V阀开启压力对呼吸排放因子的影响,降低加油站呼吸排放因子,减少呼吸排放的油气量。
4) 我国在2020年开始实施GB 18352.6-2016标准,具备车载油气回收ORVR功能的车辆将越来越多。有必要对ORVR与二次油气回收系统的不兼容排放因子开展研究,可为ORVR兼容型油气回收控制技术的研发提供数据支撑。
-
表 1 4种砷化物的质量浓度与峰面积的线性关系
Table 1. Linear relation between mass concentration and peak area of four arsenic species
砷形态Arsenic species 线性方程Linear regression equation 相关系数R2Correlation coefficient 亚砷酸盐As(Ⅲ) y=610.65x-139.35 0.9991 二甲基砷酸DMA y=1137.51x-222.29 0.9995 一甲基砷酸MMA y=941.02x+39.89 0.9997 砷酸盐As(Ⅴ) y=579.71x-105.93 0.9996 表 2 海苔样品中砷形态化合物检测结果(µg·kg−1)
Table 2. Test results of arsenic speciation in nori samples (µg·kg−1)
样品Sample As(Ⅲ) As(Ⅴ) AsB DMA MMA 砷糖-OH Arsenosugar-OH 总砷 Total arsenic 回收率/%Recovery S02海苔 ND ND 457.9 8498.6 ND 14196.3 27175.5 85.2 S06海苔 ND ND 1117.7 4540.5 ND 2609.9 12347.1 67.0 S12海苔 ND ND 408.8 3756.8 ND 10980.0 20598.9 73.5 S14海苔 ND ND 953.7 2279.1 ND 3640.6 7957.9 86.4 表 3 米饭样品中砷形态化合物检测结果(µg·kg−1)
Table 3. Test results of arsenic speciation in rice samples (µg·kg−1)
样品Sample As(Ⅲ) DMA DMMTA As(Ⅴ) iAs Unknow 总砷Total arsenic 回收率/%Recovery S02米饭 60.5 38.4 17.9 16.9 77.4 11.2 158.3 91.5 S06米饭 50.4 20.8 15.5 26.8 77.2 ND 116.4 97.5 S12米饭 75.6 20.0 11.1 13.3 88.9 ND 127.0 94.5 S14米饭 43.3 25.8 13.4 ND 43.3 8.8 98.4 92.8 表 4 寿司样品中砷形态化合物检测结果(µg·kg−1)
Table 4. Test results of arsenic speciation in sushi samples (µg·kg−1)
样品 Sample AsB As(Ⅲ) DMA 砷糖-OH Arsenosugar-OH As(Ⅴ) iAs Unknow 总砷Total arsenic 回收率/%Recovery S02寿司 35.2 45.0 488.2 143.6 ND 45.0 24.6 1567.7 47.0% S06寿司 57.4 75.2 49.0 30.6 26.7 101.9 35.0 416.4 61.1% S12寿司 19.3 58.9 109.0 315.3 ND 58.9 ND 654.6 76.8% S14寿司 36.2 57.2 98.2 117.7 ND 57.2 ND 351.3 88.0% 注:ND, 未检出. ND, not detected. -
[1] 陈保卫, LE X C. 中国关于砷的研究进展 [J]. 环境化学, 2011, 30(11): 1936-1943. CHEN B W, LE X C. Recent progress in arsenic research in China [J]. Environmental Chemistry, 2011, 30(11): 1936-1943(in Chinese).
[2] NAUJOKAS M F, ANDERSON B, AHSAN H, et al. The broad scope of health effects from chronic arsenic exposure: Update on a worldwide public health problem [J]. Environmental Health Perspectives, 2013, 121(3): 295-302. doi: 10.1289/ehp.1205875 [3] ANETOR J I, WANIBUCHI H, FUKUSHIMA S. Arsenic exposure and its health effects and risk of cancer in developing countries: Micronutrients as host defence [J]. Asian Pacific Journal of Cancer Prevention:APJCP, 2007, 8(1): 13-23. [4] RAI P K, LEE S S, ZHANG M, et al. Heavy metals in food crops: Health risks, fate, mechanisms, and management [J]. Environment International, 2019, 125: 365-385. doi: 10.1016/j.envint.2019.01.067 [5] 胥佳佳, 冯鑫, 汤静, 等. 超声辅助提取-高相液相色谱-电感耦合等离子体质谱法测定香菇中6种形态砷化合物 [J]. 食品科学, 2016, 37(24): 216-221. doi: 10.7506/spkx1002-6630-201624034 XU J J, FENG X, TANG J, et al. Arsenic speciation analysis of Lentinus edodes by ultrasonic-assisted extraction-high performance liquid chromatography-inductively coupled plasma mass spectrometry [J]. Food Science, 2016, 37(24): 216-221(in Chinese). doi: 10.7506/spkx1002-6630-201624034
[6] ALEXANDER J, BENFORD D, BOOBIS A, et al. Scientific opinion on arsenic in food. EFSA panel on contaminants in the food chain (CONTAM)[J]. EFSA Journal, 2009, 7(10: 1351): 1-198. [7] 张维, 齐丽娟, 宁钧宇, 等. 砷的健康危害评估 [J]. 毒理学杂志, 2021, 35(5): 367-372,378. doi: 10.16421/j.cnki.1002-3127.2021.05.002 ZHANG W, QI L J, NING J Y, et al. Health hazard assessment of arsenic [J]. Journal of Toxicology, 2021, 35(5): 367-372,378(in Chinese). doi: 10.16421/j.cnki.1002-3127.2021.05.002
[8] 刘蕊芬. 中日两国饮食疗法(药膳)的源流与异同的研究[D]. 广州: 广州中医药大学, 2007. LIU R F. Study on the origin of and similarities and differences between Chinese and Japanese dietotherapy[D]. Guangzhou: Guangzhou University of Chinese Medicine, 2007(in Chinese).
[9] JECFA (Joint FAO/WHO Expert Committee on Food Additives). Safety evaluation of certain food additives and contaminants. WHO Food Additives Series No. 63, Prepared by the Seventy-second Meeting of JECFA[R]. World Health Organization, Geneva, 2011. [10] TAYLOR V, GOODALE B, RAAB A, et al. Human exposure to organic arsenic species from seafood [J]. Science of the Total Environment, 2017, 580: 266-282. doi: 10.1016/j.scitotenv.2016.12.113 [11] FRANSISCA Y, SMALL D M, MORRISON P D, et al. Assessment of arsenic in Australian grown and imported rice varieties on sale in Australia and potential links with irrigation practises and soil geochemistry [J]. Chemosphere, 2015, 138: 1008-1013. doi: 10.1016/j.chemosphere.2014.12.048 [12] U. S. EPA. Method 3050B: Acid Digestion of Sediments, Sludges, and Soils[S]. 1996. [13] CAUMETTE G, KOCH I, REIMER K J. Arsenobetaine formation in plankton: A review of studies at the base of the aquatic food chain [J]. Journal of Environmental Monitoring, 2012, 14(11): 2841-2853. doi: 10.1039/c2em30572k [14] YANG G D, ZHENG J P, CHEN L, et al. Speciation analysis and characterisation of arsenic in lavers collected from coastal waters of Fujian, south-eastern China [J]. Food Chemistry, 2012, 132(3): 1480-1485. doi: 10.1016/j.foodchem.2011.12.006 [15] 王亚, 张春华, 葛滢. 高效液相色谱-氢化物发生-原子荧光光谱法检测紫菜中的砷形态 [J]. 分析试验室, 2013, 32(5): 34-38. doi: 10.13595/j.cnki.issn1000-0720.2013.0125 WANG Y, ZHANG C H, GE Y. Determination of arsenic speciation in Porphyra using HPLC-(UV)-HG-AFS [J]. Chinese Journal of Analysis Laboratory, 2013, 32(5): 34-38(in Chinese). doi: 10.13595/j.cnki.issn1000-0720.2013.0125
[16] MOLIN M, ULVEN S M, MELTZER H M, et al. Arsenic in the human food chain, biotransformation and toxicology - Review focusing on seafood arsenic [J]. Journal of Trace Elements in Medicine and Biology, 2015, 31: 249-259. doi: 10.1016/j.jtemb.2015.01.010 [17] 国家卫生和计划生育委员会, 国家食品药品监督管理总局. 食品安全国家标准 食品中污染物限量: GB 2762—2017[S]. 北京: 中国标准出版社, 2017. National Health and Family Planning Commission of the People's Republic of China, China Food and Drug Administration. National Food Safety Standard. Maximum Levels of Contaminants in Food: GB 2762—2017[S]. Beijing: Standards Press of China, 2017(in Chinese).
[18] SUN S K, XU X J, TANG Z, et al. A molecular switch in sulfur metabolism to reduce arsenic and enrich selenium in rice grain [J]. Nature Communications, 2021, 12: 1392. doi: 10.1038/s41467-021-21282-5 [19] MAWIA A M, HUI S Z, ZHOU L, et al. Inorganic arsenic toxicity and alleviation strategies in rice [J]. Journal of Hazardous Materials, 2021, 408: 124751. doi: 10.1016/j.jhazmat.2020.124751 [20] BORAK J, HOSGOOD H D. Seafood arsenic: Implications for human risk assessment [J]. Regulatory Toxicology and Pharmacology, 2007, 47(2): 204-212. doi: 10.1016/j.yrtph.2006.09.005 [21] 黄东仁. 福建省紫菜中砷的形态及含量[J]. 中国渔业质量与标准, 2015, 5(1): 52-60. HUANG D R. Arsenic speciations and their contents in Porphyra in Fujian Province[J]. Chinese Fishery Quality and Standards, 2015, 5(1): 52-60(in Chinese). HUANG D R. Arsenic speciations and their contents in Porphyra in Fujian Province[J]. Chinese Fishery Quality and Standards, 2015, 5(1): 52-60(in Chinese).
[22] YANG G D, XIE W Y, ZHU X, et al. Effect of arsenite-oxidizing bacterium B. laterosporus on arsenite toxicity and arsenic translocation in rice seedlings [J]. Ecotoxicology and Environmental Safety, 2015, 120: 7-12. doi: 10.1016/j.ecoenv.2015.05.023 [23] WU Y C, ZHANG H, WANG K T, et al. Metabolic and residual characteristic of different arsenic species contained in laver during mouse digestion [J]. Science of the Total Environment, 2021, 793: 148434. doi: 10.1016/j.scitotenv.2021.148434 [24] DAI J, CHEN C, GAO A X, et al. Dynamics of dimethylated monothioarsenate (DMMTA) in paddy soils and its accumulation in rice grains [J]. Environmental Science & Technology, 2021, 55(13): 8665-8674. [25] DAI J, TANG Z, GAO A X, et al. Widespread occurrence of the highly toxic dimethylated monothioarsenate (DMMTA) in rice globally [J]. Environmental Science & Technology, 2022, 56(6): 3575-3586. [26] ZHANG J, BARAŁKIEWICZ D, WANG Y Z, et al. Arsenic and arsenic speciation in mushrooms from China: A review [J]. Chemosphere, 2020, 246: 125685. doi: 10.1016/j.chemosphere.2019.125685 [27] 程家丽, 任硕, 刘婷婷, 等. 2001—2017年我国部分地区蔬菜中砷和重金属累积特征及膳食暴露风险 [J]. 中国食品卫生杂志, 2018, 30(2): 187-193. doi: 10.13590/j.cjfh.2018.02.013 CHENG J L, REN S, LIU T T, et al. Accumulation and dietary exposure risk of arsenic and heavy metals in the vegetables from some areas of China, 2001-2017 [J]. Chinese Journal of Food Hygiene, 2018, 30(2): 187-193(in Chinese). doi: 10.13590/j.cjfh.2018.02.013
[28] ZHAO D, WANG J Y, YIN D X, et al. Arsanilic acid contributes more to total arsenic than roxarsone in chicken meat from Chinese markets [J]. Journal of Hazardous Materials, 2020, 383: 121178. doi: 10.1016/j.jhazmat.2019.121178 [29] JECFA (Joint FAO/WHO Expert Committee on Food Additives). Evaluation of certain contaminants in food[Z]. 2017. -