微波消解-电感耦合等离子体质谱法测定土壤样品中重金属含量的不确定度评定

刘 卫，段路路，翟 帅

(上海化工院检测有限公司/上海化学品公共安全工程技术研究中心 上海 200062)

随着工业和经济的发展，大量的污染物进入土壤环境，土壤中的重金属已对生态环境造成了严重的污染[1-2]。重金属不易被微生物分解，易与土壤中的有机质发生螯合作用并稳定存在于土壤中，通过农作物被人体吸收，危害极大[3-4]。土壤中重金属元素的测定主要采用石墨炉原子吸收光谱法、电感耦合等离子体原子发射光谱法、原子荧光光谱法、火焰原子吸收光谱法、电感耦合等离子体质谱法(ICP-MS)等[5]。ICP-MS法具有检出限低、线性范围宽、干扰少、灵敏度高、可同时测定几十种元素等优点[6]，痕量分析检测下限可达10-12，近年来随着该技术的不断推广，ICP-MS法已被许多发达国家应用于环境监测中。

依据《测量不确定度评定与表示》(JJF 1059.1—2012)[7]，测量不确定度定义是根据所获信息，表征赋予被测量值分散性的非负参数。在实际工作中，测量不确定度对评定结果具有重要意义[8-12]。本文采用微波消解-ICP-MS法测定土壤中重金属元素铬、铜、砷、铅的含量，根据JJF1059.1—2012建立针对ICP-MS法测定土壤中重金属含量的较为完整的不确定度评定方案，以期找到影响不确定度的主要因素，保证测量结果的准确性，为土壤安全监测和控制提供科学的指导。

1 材料与方法

1.1 主要仪器与试剂

主要仪器：ICAP RQ型电感耦合等离子体质谱仪，赛默飞世尔科技公司；CEM MARS6型微波消解仪，美国培安公司；Simplicty型超纯水器，德国默克密理博公司；ME204-2型电子天平，瑞士梅特勒公司。

主要试剂：超纯水，电阻率18.2 MΩ·cm；硝酸，优级纯，国药集团化学试剂有限公司；氩气，体积分数为99.996%；26种元素质量控制ICP-MS混合标准储备溶液，100 mg/L，CFGG-160008-01-01。

1.2 仪器工作参数

等离子体功率：1 550 W；冷却气(氩气)流量：15 L/min；辅助气(氩气)流量：0.8 L/min；雾化气(氩气)流量：1 L/min；总测试时间：180 s；停留时间：0.05 s。

1.3 试验方法

1.3.1 标准溶液和内标溶液的制备

移取100 mg/L混合标准储备溶液1.00 mL，用2%(体积分数，下同)硝酸溶液稀释，配制成1.00 mg/L标准溶液。分别移取0.050、0.100、0.500、1.00、2.00、5.00、10.00 mL标准溶液于100 mL容量瓶中，用2%硝酸溶液稀释配制成质量浓度为0.50、1.00、5.00、10.0、20.0、50.0、100.0 μg/L的标准溶液。内标溶液为73Ge 20 μg/L、193Ir 20 μg/L。

1.3.2 样品制备

用万分之一天平称取0.25 g(精确至0.000 1 g)土壤样品置于微波消解罐中，加入3 mL硝酸、1 mL盐酸，进行微波消解(消解程序为室温升温至185 ℃，升温时间15 min，保持40 min)。待消解完全后，将消解后的溶液转移至50 mL容量瓶中，用2%硝酸溶液洗涤消解罐3～5次，洗液合并至容量瓶(洗涤过程中防止溶液损失)，用2%硝酸溶液稀释至刻度、混匀。用0.22 μm水相滤膜过滤，即为待测溶液。同时做空白试验。

2 不确定度评定

2.1 数学模型的建立

数学模型如式(1)所示：

(1)

式中：w——试样中被测重金属元素的质量分数，mg/kg；

ρ——由标准曲线计算得到的试样溶液中被测重金属元素的质量浓度，mg/L；

V——试样消解液的定容体积，mL；

f——试样溶液的稀释因子，无量纲；

m——试样称样量，g。

2.2 不确定度来源分析

由数学模型式(1)可知，影响ICP-MS法测定土壤中重金属含量的不确定度主要来源如下。

(1)试样称量：由称量质量的重复性、电子天平分辨力和示值误差引入的不确定度。

(2)被测元素的质量浓度：由最小二乘法拟合引起的不确定度。

(3)样品消解回收率：在土壤样品微波消解过程中，因受样品性质、消解条件、元素属性等影响，会引起被测金属元素在容器内附着、污染、挥发以及消解不完全等情况，使得金属元素损失。

(4)试样消解液的定容体积：由容量瓶的误差以及因温度变化导致溶液体积变化引入的误差。

(5)标准溶液配制：包括标准储备溶液和标准溶液配制引入的不确定度，以及在配制过程中容量瓶和移液器引入的不确定度。

(6)测量重复性：由仪器稳定性引起的测定平行样品的不确定度。

综上所述，采用ICP-MS法测定土壤中重金属含量的不确定度来源因果关系如图1所示。

图1 ICP-MS法测定土壤中重金属含量的不确定度来源因果关系

2.3 标准不确定度分量评定

2.3.1 试样称样量m引入的标准不确定度分量u(m)

影响称样量m的因素主要有称量质量的重复性mr和电子天平的分辨力δm。修正后的质量见式(2)：

m=mr+δm

(2)

(1)称量质量的重复性mr引入的标准不确定度u(mr)

重复称量质量为0.250 0 g的样品10次，所得结果见表1。

表1 试样重复称量10次的结果

根据表1结果计算可得，称样量的标准偏差为0.12 mg，则u(mr)=s(m)=0.12 mg。

(2)电子天平的分辨力引入的标准不确定度u(δm)

以上两项合成得到u(m)：

试样称样量m的平均值为0.250 0 g，则其相对标准不确定度为：

2.3.2 消解液的定容体积引入的标准不确定度分量u(V1)

影响试样消解液定容体积的因素主要有容量瓶示值校准VE和容量瓶与溶液温度差引入的误差δVt，修正后的体积见式(3)：

V1=VE+δVt

(3)

(1)容量瓶示值校准引入的标准不确定度u(VE)

(2)温度影响引入的标准不确定度u(δVt)

以上两项合成，得到u(V1)：

2.3.3 试样溶液中被测元素质量浓度的标准不确定度分量u(ρ)

2.3.3.1 标准曲线拟合引入的标准不确定度分量u(ρ1)

配制的标准工作溶液质量浓度依次为0.50、1.00、5.00、10.0、20.0、50.0、100.0 μg/L，按照ICP-MS分析条件分别进样3次，共计21次，测定结果(见表2)用最小二乘法拟合得到线性方程为y=bρ+a(a为截距，b为斜率，y为ICP-MS测得的信号强度)。

表2 各标准溶液质量浓度与对应的信号强度

拟合线性回归方程的试验标准偏差u(y)按式(4)计算，回归系数a、b的试验标准偏差s(a)、s(b)分别按式(5)和式(6)计算：

(4)

(5)

(6)

式中：yi——标准曲线各个浓度点对应的信号强度；

yfi——各浓度点通过拟合线性方程计算出的数值；

n——配制7个标准溶液浓度点，每个浓度点测试3次，共计21次，即n=21；

实测质量浓度ρ1由最小二乘法拟合引入的标准不确定度u(ρ1)按式(7)计算：

(7)

式中：p——试样平行测定次数，p=10；

ρp——试液中待测元素的测定值，μg/L。

标准不确定度u(ρ1)计算较为复杂，利用Excel软件中的Linest函数，依次选择信号强度y和质量浓度ρ对应的单元格，以公式单元格开始的区域，按F2，再按Ctrl+Shift+Enter。铬、铜、砷和铅的Linest函数参数见表3。

表3 各元素的Linest函数参数

2.3.3.2 样品消解引入的不确定度分量urel(R)

由于土壤样品消解不完全或消解过程中导致待测金属元素损失等，依据极差法进行消解过程引入的不确定度评定，按式(8)和式(9)计算：

s(R)=(R大-R小)/C

(8)

(9)

式中：R大——回收率最大值，%；

R小——回收率最小值，%；

n′——平行测定次数，n′=6;

C——极差系数，平行测定6次时C=2.53；

s(R)——回收率标准偏差，%。

称取土壤平行样共计6份，按1.3.2的步骤进行微波消解，用ICP-MS检测。试验样品中铬、铜、砷和铅的回收率及不确定度见表4。

表4 铬、铜、砷和铅的回收率及相对标准不确定度

2.3.3.3 标准溶液配制引入的标准不确定度分量u(ρ2)

标准溶液配制中包含标准储备溶液配制的不确定度u(ρ储)和标准溶液配制的不确定度u(ρ曲)，建立的经验模型见式(10)：

(10)

式中：ρ曲——标准溶液的质量浓度，mg/L；

ρ储——标准储备溶液的质量浓度，mg/L；

Vi——移液器取液量，mL；

V曲——标准溶液定容体积，V曲=50 mL。

(1)标准储备溶液配制引入的标准不确定度u(ρ储)

(2)标准溶液配制引入的不确定度u(ρ曲)

①定容引入的标准不确定度

②移液器引入的标准不确定度

表5 移液器不同取液量引入的不确定度

移液器7次分取标准溶液的相对标准不确定度urel(V移)为：

urel(V移)=[(0.008 66)2+(0.001 15)2+

(0.002 89)2+(0.005 77)2+(0.001 15)2+

(0.001 73)2+(0.003 46)2]1/2=0.011 6

标准溶液配制引入的相对标准不确定度为：

urel(ρ曲)=[(0.002 5)2+(0.000 7)2+(0.011 6)2]1/2=0.011 9

2.3.3.4 样品测量重复性的标准不确定度分量u(x)

根据ICP-MS的分析条件，对土壤样品进行10次测定，以最小二乘法拟合数据，结果见表6。

根据表6计算得铬、铜、砷、铅质量浓度的平均值依次为69.37、21.08、12.08、30.05 μg/L。

相对标准不确定度的计算见式(11)：

(11)

表6 铬、铜、砷、铅的重复性测定结果mg/kg

式中：xi——第i个样品的测定结果，mg/kg;

q——测定次数，q=10;

s(x)——测定结果的标准偏差，mg/kg。

铬的相对标准不确定度urel(xCr)=0.008 79；

铜的相对标准不确定度urel(xCu)=0.012 50；

砷的相对标准不确定度urel(xAs)=0.009 15；

铅的相对标准不确定度urel(xPb)=0.010 95。

综上所述，被测重金属元素质量浓度ρ的相对标准不确定度计算见式(12)：

(12)

各元素质量浓度的相对标准不确定度：

2.4 合成标准不确定度的评定

合成相对标准不确定度的计算见式(13)：

(13)

则可计算合成相对标准不确定度：

ucrel(wCr)=[(0.019 0)2+(0.000 7)2+(0.000 5)2]1/2=0.019 0

ucrel(wCu)=[(0.030 3)2+(0.000 7)2+(0.000 5)2]1/2=0.030 3

ucrel(wAs)=[(0.042 1)2+(0.000 7)2+(0.000 5)2]1/2=0.042 1

ucrel(wPb)=[(0.021 2)2+(0.000 7)2+(0.000 5)2]1/2=0.021 2

各元素合成标准不确定度为：

uc(wCr)=ucrel(wCr)·wCr=0.019 0×13.874=

0.264(mg/kg)

uc(wCu)=ucrel(wCu)·wCu=0.030 3×4.217=

0.128(mg/kg)

uc(wAs)=ucrel(wAs)·wAs=0.042 1×2.416=

0.102(mg/kg)

uc(wPb)=ucrel(wPb)·wPb=0.021 2×6.010=

0.127(mg/kg)

以铬为例，不确定度概算见表7。重金属元素铬、铜、砷、铅的不确定度分量对合成标准不确定度的贡献如图2所示。

表7 铬的不确定度概算

图2 铬、铜、砷、铅不确定度分量对合成标准不确定度的贡献

2.5 扩展不确定度评定

取包含因子k=2，则扩展不确定度如下：

U(wCr)=k·uc(wCr)=2×0.264=

0.528(mg/kg)

U(wCu)=k·uc(wCu)=2×0.128=

0.256(mg/kg)

U(wAs)=k·uc(wAs)=2×0.102=

0.204(mg/kg)

U(wPb)=k·uc(wPb)=2×0.127=

0.254(mg/kg)

2.6 测定结果报告

wCr=(13.874±0.528) mg/kg，k=2；

wCu=(4.217±0.256) mg/kg，k=2；

wAs=(2.416±0.204) mg/kg，k=2；

wPb=(6.010±0.254) mg/kg，k=2。

3 结语

通过建立ICP-MS法测定土壤中重金属元素含量的不确定度评价模型，对土壤中铬、铜、砷、铅的测定结果进行了不确定度评定。结果表明，土壤中重金属铬的含量为(13.874±0.528) mg/kg，k=2；铜的含量为(4.217±0.256) mg/kg，k=2；砷的含量为(2.416±0.204) mg/kg，k=2；铅的含量为(6.010±0.254) mg/kg，k=2。铬和铅的测量不确定度来源主要是样品消解过程；铜的测量不确定度来源主要是样品测定过程中标准曲线的拟合；砷的测量不确定度来源主要是样品测定过程中标准曲线的拟合和样品消解过程。因此在实际测定过程中，可通过增加标准曲线的浓度水平数以减少标准曲线拟合带来的影响。为了减少样品消解过程对试验结果的影响，消解前应将所使用的消解罐和玻璃容器用体积分数约为10%的硝酸溶液浸泡，然后用自来水和超纯水冲洗3～5次，置于干净的环境中晾干；检验所采购硝酸的质量，避免带入金属杂质；溶液转移过程中应将消解罐冲洗3～5次，同时避免消解液逸出造成损失。