X射线荧光光谱分析方法的检出限测量方法研究

袁良经， 贾云海, 2， 程大伟, 2*

1. 钢研纳克检测技术股份有限公司， 北京 100094 2. 钢铁研究总院， 北京 100081

引 言

X荧光光谱分析技术由于其无损、 快捷、 准确的特点， 常用于地质[1]、 矿石、 土壤、 合金[2]等样品的快速检测。 XRF仪器能够用于地质勘查的现场分析， 以及岩心样品的无损分析， 进行地球化学研究[3]， 大大缩短工作周期。 像ICP， AAS， ICP-MS等分析仪器一般都需要先消解样品， 不适于现场分析， 但是XRF技术也有局限性， XRF的使用取决于分析物的含量水平和其他因素， 主要检测mg·kg-1级别以上含量的元素。

检出限(LOD)是仪器性能评价或方法建立的重要指标， 在实际样品分析时， 小于检出限的含量， 不可检出； 大于检出限、 小于定量限的含量， 可定性分析； 大于定量限的含量， 可定量分析。 在X荧光光谱测量中， 由于X射线光子的总计数服从泊松分布， 是一种非连续分布， 仅当总计数足够大的时候， 近似于正态分布。 在实际分析中， 往往不会长时间测量去积累足够大的数。 本文以X荧光光谱分析为例， 比较几种检出限计算方法的异同。

1 多种检出限测量及计算方法

检出限的定义和典型计算方法有以下几种： 泊松分布计算法、 K倍标准偏差法、 线性校准法、 RSD法、 SD直线外推法、 环境监测推荐计算法、 海洋监测规范计算法等。

1.1 泊松分布计算法

XRF测量值只有在计数足够大的时候才能接近正态分布， 是离散分布。

在X射线荧光光谱测量中， X射线光子的计数N服从泊松分布[4]， 如式(1)

(1)

式(1)中，P(N)为计数为N的概率；μ为多次测量的平均值； 总计数为N的单次测量标准偏差σN见式(2)。 当总计数足够大时， 泊松分布趋近于正态分布。

(2)

(3)

根据含量和峰强度绘制的工作曲线： c=I/m+b， c为含量； m为灵敏度， 表示单位浓度的计数率， 1/m为斜率； b为截距。 根据误差传递

(4)

(5)

1.2 K倍标准偏差法

图1 检出限LOD和临界含量XCFig.1 Low limit of detection and critical value

1.3 线性校准法

根据GB/T 33260.2—2018[8]， 可采用线性校准情形确定特定物质的存在性， 估计净状态变量的临界值、 响应变量的临界值和净状态变量的最小可检出值(检出限)。 绘制校准直线的标样数I至少是3， 一般取5； 每个标样制备数J相同， 一般是2； 实际被测样品制备数K应和绘制校准直线的标样的制备数J相同； 每个样品重复测量次数L相同， 一般测2次。

临界含量和检出限的计算步骤参照GB/T 22554—2010的假定： 校准函数是线性的； 响应变量的测量值相互独立且服从正态分布， 响应变量的标准差即为“误差标准差”； 误差标准差可分为常数标准差和变量标准差[9]。 本文仅考虑标准差为常数的情形。 线性校准模型为

Yij=a+bxi+εij

(6)

式(6)中：xi为状态i的变量值；εij为误差的随机变量， 相互独立，εij～N(0，σ2)， 响应变量Yij的均值E(Yij)=a+bxi， 方差V(Yij)=σ2不依赖于xi。

a，b和σ2的估计值如式(7)—式(9)

(7)

(8)

(9)

净状态变量的临界值

(10)

式(10)中，t0.95(ν)是自由度为(IJ-2)的t分布的95%分位数。

最小可检出值(检出限)

(11)

式(11)中，δ=δ(ν;α;β) 为非中心参数， 当α=β且ν>3时，δ(ν;α;β)≈2t1-α(ν)； 当α=β=0.05且ν>3时， 最小可检测值xd近似为

(12)

1.4 按相对标准偏差RSD法计算检出限

根据GB/T33260.5—2018《检出能力第5部分非线性校准情形检出限的确定方法》[10]， 在校准曲线是非直线并且检出限的标准偏差是常数时， 可用测量含量的相对标准偏差RSD=30%对应的含量为检出限。 其实当校准曲线为直线时也可以将RSD=30%对应的含量为检出限， 因为直线是曲线的一种特殊形式。

测量值RSD为30%时对应仪器或方法的检出限。 这些定义和推导都是针对多次分析总体呈正态分布的分析， 用于小样本分析时， RSD的数值会发生变化。 贾云海等[5]推导出10次空白样品测量的RSD≤43%时， 被测物有95%的置信度检出， RSD=43%时的含量对应着检出限LOD的上限。 当测量值的RSD>43%时， 不可定性检出， 只有当测量值的RSD≤43%时， 测量值可以被定性检出。 该条件可作为检出限判断的必要条件， 但不是充分条件。

1.5 测量标准偏差SD直线外推法计算检出限

由于固体样品不易形成不含被测元素的样品空白， 可用模拟空白评估法， 做关于标准偏差SD与含量的关系线。 对Pb， As和Cd元素拟合标准偏差SD-含量曲线， 求出模拟空白对应的S0值， 外推求出含量为0时的标准偏差S0， 以3S0为检出限。

1.6 按环境监测分析技术导则计算方法检出限

HJ 168-2020《环境监测分析方法标准制订技术导则》中规定方法检出限是“用特定分析方法在给定的置信度内可从样品中定性检出待测物质的最低浓度或最小量”[6]。 导则强调“定性检出”， 即判定样品中有高于空白的浓度或含量的待测物质。 按照样品分析的步骤， 重复测量n≥7次空白样品， 计算n次平行测定的标准偏差， 按式(13)

MDL=t(n-1, 0.99)×S

(13)

式(13)中， MDL是方法的检出限，n是平行测量次数，t是单侧分布自由度为(n-1)对应的值， 置信度为99%。S是n次测量的标准偏差。

1.7 按海洋监测规范计算方法检出限

GB 17378.2—2007规定检出限(XN)是“通过一次测量， 就能以95%的置信概率定性判定待测物质存在所需要的最小浓度或量”[7]。 检出限按式(14)估算

(14)

2 实验部分

2.1 仪器和样品

实验仪器： 食品重金属检测仪NX-100FA， 程序1测试Pb和As元素， 测量时间10 min； 程序2测试Cd元素， 测量时间10 min。

标准样品： 包括大米粉标样GBW(E)100348—362和AP3参考样品(ICP-MS定值)含量如表1所示。

表1 标准样品和参考样品含量(μg·g-1)Table 1 Content of standard and reference sample (μg·g-1)

2.2 方法

工作曲线绘制： 在预设程序下分别测定表2中的6个参考样品各两次， 采集Pb， As和Cd元素的强度， 考虑元素干扰和基体效应绘制强度-含量工作曲线， 如图2， 决定系数R2>0.99。

图2 大米粉中Pb(a)， As(b)和Cd(c)元素校准曲线Fig.2 Calibration curves of Pb (a), As (b) and Cd (c) in rice

上述七种检出限计算方法需要测试的内容如表2。

表2 七种检出限计算方法及需要测试方法Table 2 Seven LOD methods of calculation and detection

3 结果与讨论

根据检出限计算方法和表2中的测试内容， 分别计算相应方法的检出限。

3.1 泊松分布

本文用强度比绘制校准曲线， 在计算检出限时， 需要考虑参比峰的强度IN， 将式(5)变换为式(15)

(15)

式(15)中，m′为单位浓度的计数率比；IN为参考峰的计数率， s-1；IB为元素的背景强度， 表示每秒的光子计数；T为测量时长， s。

在实际计算X射线荧光方法或仪器检出限时， 可只测试高低两个梯度样品， 计算相应灵敏度系数； 由于先绘制了工作曲线， 本方法直接采用图2工作曲线的灵敏度系数。 计算结果如表3。

表3 泊松分布法计算检出限结果(μg·g-1)Table 3 LOD result in Poisson distribution method (μg·g-1)

3.2 K倍标准偏差法

实际样品中很难找到绝对空白样品， 采用GBW(E)100361作为Pb元素近似空白样品， 采用GBW(E)100362作为As和Cd元素近似空白样品。 该法计算检出限需要先绘制校准曲线， 再测试10次空白样品， 计算结果如表4。

表4 K倍标准偏差法检出限计算结果(μg·g-1)Table 4 LOD result in K times standard deviation method (μg·g-1)

3.3 线性校准法

取5个有梯度的标准或参考样品， 利用线性校准法计算检出限。I=5， 如表5所示， “/”表示该样品不参与绘线。

表5 5个参考样品Pb， As， Cd元素含量(μg·g-1)Table 5 Pb, As, Cd content of five reference sample (μg·g-1)

每种标准或参考样品制备2个，J=2； 每个标准或参考样品测试两次，L=2； 被测试样测试次数相同，K=2。 线性校准法检出限Xd计算结果如表6。

表6 线性校准法检出限计算结果(μg·g-1)Table 6 LOD result in the linear calibration method (μg·g-1)

3.4 相对标准偏差RSD法

选取5个梯度样品， 每个样品分别测试10次， 计算RSD， 如表7。

表7 五个梯度样品10次测试的RSD/%Table 7 RSD of ten measurements for five samples

表8 RSD法检出限计算结果(μg·g-1)Table 8 LOD result in RSD method(μg·g-1)

Pb， As和Cd元素的RSD-C曲线如图3所示。

图3 大米粉中Pb(a), As(b)和Cd(c)元素RSD-C曲线Fig.3 RSD-C curves for Pb (a), As (b), Cd (c) in rice

3.5 标准偏差SD法

与3.4节使用相同测量数据， 计算10次测量数据的标准偏差SD， 如表9。

表9 标准偏差SD外推法检出限计算结果(μg·g-1)Table 9 LOD result in standard deviation extrapolation (μg·g-1)

由表9可知， SD-含量线性较差， 斜率小于截距的一半； 当含量小于1 μg·g-1， 尤其是样品含量段在0.007～0.84 μg·g-1时， 截距比重越大， 决定系数R2越小， 检出限数值的可靠性越差。 SD-含量线性与RSD-含量线性的决定系数R2对比可以看出， RSD-含量拟合的更好， 说明RSD-含量的关系规律更为稳健。

3.6 环境监测分析技术导则计算方法

按照样品分析的步骤， 重复测量10次近似空白样品， 按式(13)计算平行测定的标准偏差， 测量次数n=10时，t=2.821， 计算结果如表10。

表10 环境监测分析检出限计算(μg·g-1)Table 10 LOD result in environmental monitoring and analysis (μg·g-1)

3.7 海洋监测分析

与3.5节使用同组数据， 按照1批10次测试， 计算标准偏差Swb； 在自由度γ=9、α=0.05时， 单边检验t=1.833。 按海洋监测分析计算检出限， 结果如表11。

表11计算的检出限含量相对较高， 主要是因为使用该方法测试时， 每个样品都需单独扣空白。

表11 海洋监测分析检出限计算(μg·g-1)Table 11 LOD result in marine monitoring and analysis(μg·g-1)

3.8 七种检出限测量及计算

选用七种检出限计算方法独立计算， 检出限符号和计算公式区别较大， 但计算结果接近， 结果见表12。 表中除RSD=43%对应含量为95%置信度下检出限的上限外， 其余符号均为检出限。

表12 七种检出限计算结果比较(μg·g-1)Table 12 Comparison of calculation results in seven detection limits cases (μg·g-1)

从表12可以看出， 除海洋监测规范法外， 其他方法测量的检出限数值接近；K倍标准偏差法与泊松分布法和线性校准法相比， 计算的Pb和As略低， Cd检出限相差无几， 可能是由于空白样品Pb和As有一定含量， Cd近似空白； 用这些方法都可以反映X射线光谱分析大米中Pb， As和Cd的检出能力。

4 结 论

(1)X射线荧光分析时， 由于其属于泊松分布， 应优先使用泊松分布法计算其检出限， 而且X荧光方法测试样品， 2个样品测量即可快速计算方法检出限。

(2)线性校准法， 考虑因素最多， 被认为是测量值正态分布时检出限计算较准确方法， 可将该结果作为参考值。

(3)RSD法和SD直线外推法测试次数较多， 可用于没有空白样品或不能直接获得强度时。 RSD法也可作为检出限判断的必要条件， RSD>43%不能定性检出， RSD=30%对应的含量值为检出限。