全 文 ：收稿日期：!""#$"%$!# 接受日期：!""&$"’$"(
基金项目：中央级公益性科研院所基本科研业务费专项资金（!""#$’’）资助。 作者简介：范洪黎（’&(&—），女，浙江海宁人，副研究员，博士，主要从事光谱分析、农业资源与环境研究。 )*+："’"$#!’"#(,-，./012+：3456+2715!,891344: ;40: ;5。! 通讯作者 )*+："’"$#!’"#(%" !"#$%&’法测定水溶肥料中硼含量
的不确定度分析
范洪黎，韩岩松，刘 密，王 旭!
（中国农业科学院农业资源与农业区划研究所，北京 ’"""#’）
摘要：对 <=>/?.@测定水溶肥料中硼的分析方法的测量不确定度进行了分析；建立了测定过程中各分量的数学模
型，并识别了测量过程中不确定度来源；估算了各不确定度分量对总不确定度的影响，确定了测量结果的置信区
间；给出了水溶肥料中硼的含量及其不确定度：!（A）B（"C--, D "C"E!）F（G B !）。
关键词：不确定度分析；<=>/?.@；硼含量
中图分类号：@’E$,, 文献标识码：? 文章编号：’""#$-"-H（!""&）"-$’!E"$"-
!"#$%&’& () *"+,-.#’".% ’" /,.,-0’"#.’(" () 1(-(" +(".,". ’" 2#.,-3&($*1$, ),-.’$’4,-& 1% 5673!89
I?J K456/+2，K?J L15/M456，N（ "#$%&%’%( )* +,-&.’/%’-0/ 1($)’-.($0#2 1(,&)#0/ 3/0##&#,，45&#($( +.02(67 )* +,-&.’/%’-0/ 8.&(#.($，9(&:&#, ;<<<=;，45�） !1&.-#+.：?51+9M2M 47 S5;*TU125U9 25 V*U*T0251U245 47 W4T45 ;45U*5U 25 M4+SW+* 7*TU2+2X*T W9 <=>/?.@ Y1M 25UT4VS;*V: ? 01U3*01U2;1+ 04V*+ U4 *Z1+S1U* U3* S5;*TU125U9 25 0*1MST*0*5U Y1M ;45MUTS;U*V W1M*V 45 1 ;40[T*3*5M2Z* S5V*TMU15V256 47 U3* Y34+* [T4;*VST* 47 U3* V*U*T0251U245 : )3* M4ST;*M 47 S5;*TU125U9 15V U3*2T *77*;UM 45 U3* U4U1+ S5;*TU125U9 Y*T* 1+M4 151+9X*V，U3* 0*1MST*0*5U T*MS+U 47 W4T45 ;45U*5U 25 M4+SW+* 7*TU2+2X*T Y1M 62Z*5 1M !（A）B（"C--, D "C"E!）F（G B !）: :,% 2(-/&：S5;*TU125U9 25 0*1MST*0*5U；<=>/?.@；W4T45 ;45U*5U 化学定量分析的结果可用于判定某些产品是 否符合要求，某些材料是否符合特定规范。当使用 分析结果作为决策依据的时候，必须对这些结果的 质量有所了解，换句话说，这些结果在多大程度上是 可靠的。证明结果质量的一个有效方法就是测量不 确定度。自 RA \ ) ’-E#’$!"""《检测和校准实验室能
力的通用要求》［’］发布以来，在实验室认可和计量
认证的评审中，测量不确定度的评定已成为一项重
要工作。测量不确定度是表征被赋予测量之值的分
散性，是与测量结果相联系的参数［’$! ］，是评定测 量结果质量的一个重要指标［,］。在一个完整的测量 过程中，引起测量不确定度的因素有很多，包括测量 系统的所有要素。测量系统的概念不只局限于测量 仪器，测量设备的范畴，而是指用来对被测量值赋值 的测量操作程序，测量人员、设备，环境及软件等要 素的综合，是获得测量结果的整个过程［E$- ］。因此，
一切测量结果都具有不确定度。一份完整的检测报
告应包括对检测结果的不确定度的分析［’］。本研究
根据不确定度的评估程序，以 <=>/?.@法测定水溶
肥料中硼含量为例，对试验过程中可能产生的不确
定度来源进行了识别和分析，确定了测量结果的不
确定度。
; 材料与方法
水溶肥料中硼含量的测定按照 JL ’E!#$!""% 《微量元素水溶肥料》［(］规定的方法进行。 所用仪器为等离子体发射光谱仪（日本岛津 <=>@$’"""!型真空单通道扫描型）。光学系统为双
植物营养与肥料学报 !""&，’-（-）：’!E"$’!EE """""""""""""""""""""""""""""""""""""""""""""""""""""""""""""" >+15U JSUT2U245 15V I*TU2+2X*T @;2*5;* 光栅 !"#$%&’()$%#$装置，入口狭缝 *+!,，出口狭缝
-+!,。射频发生器频率 *./0* 12"，入射功率 0/*
34，观察高度 05 ,,，提升量 */6 ,7 8 ,9%。工作气
体氩气流量：冷却气 0: 7 8 ,9%，等离子气 0/*
7 8 ,9%，载气 0/+ 7 8 ,9%，净化气 -/5 7 8 ,9%。积分时
间 5 ;。电子分析天平（瑞士梅特勒 <=0++型）。
常用玻璃器皿使用前用稀硝酸溶液浸泡 *: >
以上，自来水反复冲洗，再用去离子水冲洗，晾干后
备用。
试剂用硼的标准储备液 0+++!? 8 ,7。用此储备
液稀释至所需浓度。
! 结果与分析
测量不确定度通常由测量过程的数学模型和不
确定的传播律来评定［-’:］。因此，要评定一种测量
方法的不确定度，必须首先建立测量过程的数学模
型。
!"# 测量过程概述
利用等离子体发射光谱法（@!AB<=C）测定肥料
中的硼含量，其工作原理是：试样溶液中的硼在 @!A
光源中原子化并激发至高能态，处于高能态的原子
跃迁至基态时产生具有特征波长的电磁辐射，辐射
强度与硼原子浓度成正比。在测量过程中，首先是
配制一系列不同硼浓度的标准溶液，用等离子体发
射光谱仪在波长 *:D/E %, 处测定各标准溶液的辐
射强度。以各标准溶液的硼的质量浓度（!? 8 ,7）为
横坐标，相应的辐射强度为纵坐标，绘制校准曲线。
然后按照标准溶液的测定方法进行未知样品的测
定，并通过校准曲线求出未知样品的浓度。在整个
测量过程中，影响不确定度的主要因素有标准溶液
的配制、浓度的测定、样品称量等。下面就基于测量
过程建立的数学模型进行讨论。
!"! 建立数学模型
试样中硼的含量按下列公式进行计算：
! "
（! #!+）$% & 0+6 & 0++ （0） 式中!—由工作曲线查出的试样溶液硼的质量浓 度，单位为微克每毫升（!? 8 ,7）；!+—由工作曲线查 出的空白溶液中硼的质量浓度，单位为微克每毫升 （!? 8 ,7）；$—试样溶液总体积，单位为毫升（,7）；
%—试料的质量，单位为克（?）。
!"$不确定度来源分析 按照数学模型式（0）以及测量过程，分析其不 确定度来源有：浓度测定引入的不确定度；空白溶 液引入的不确定度；定容体积引入的不确定度；样 品称量引入的标准不确定度以及测量的重复性的不 确定度等。 !"% 相对标准不确定度分量的评定 */:/0 !的相对标准不确定度 */:/0/0 校准曲线的标准不确定度 在化学实验 室中，很多分析方法都采用线性回归计算校准曲线。 校准曲线包含了仪器本身的不确定度，求回归方程 的目的之一是解决实际当中的预测问题，即通过自 变量（&）的值求因变量（F）及其取值范围。标准溶液 配制中也会产生误差，由于标准溶液是标准系列的 基础，因此，这些随机误差将作为系统误差传递到校 准过程。而目前普遍采用的简单线性回归将 &作为 常数，&的系统误差又不能反映到 F的随机误差中， 因此，必须将标准溶液配制中的不确定度加入到 & G H I JF中进行合成。对于回归方程 & G H I JF，求 被测量 F的数学模型为；F G & ’ HJ ，因此，由 &、H、J的 -个值的不确定度，从而求出被测量 F 的不确定 度［.’E］： )*K（F）G K*（&）)*（&）I K*（J）)*（J）I K*（H）)*（H） （*） 其中，灵敏度系数： K（&）G!F !& G 0J；K（J）G !F !J G ’ "& ’ H J* ； K（H）G!F !H G ’ 0J 进行测定前，选择仪器最佳工作条件。然后依 据 */0 的测量过程完成校准曲线的绘制。用 @!AB <=C法测得硼的校准曲线 ’ G +/+-.- I +/5*6:!，! 对应 F，’ 对应 &，H G +/+-.-，J G +/5*6:，相关系数$
G +/DDDD。硼校准曲线的具体数值见表 0。
依据式（*），利用所得实验数据计算被测量 F的
不确定度。
0）求 )（&）及 K（&）：
)（&）G C& 8 F G
#
%
9 G 0
#
,
L G 0
（&9L ’ &9）*
(%$’ * G +/0-5: 其中，%和 ,分别为工作曲线测量点数及测量次数。 灵敏度系数：K（&）G!F !& G 0J G 0/D++ *）求 )（J）及 K（J）： )（J）G C& 8 F #（F9 ’ F）$ *
G */D-E M 0+’ -
灵敏度系数：K（J）G!F
!J
G
"& ’ H
J*
G ’ -E/6-
0:*05期 范洪黎，等：@!AB<=C法测定水溶肥料中硼含量的不确定度分析
!）求 "（#）及 $（#）： "（#）% &’ ( ) !)*+ !"!（)+ ,")）# * % -./*01 灵敏度系数：$（#）%
$)$#
% , 12 % , 1.3--
/）"$（)）%$*（’）"*（’）4 $*（2）"*（2）4$*（#）"*（## ）
% 1.3--* 5 -.1!0/* 4（ , !678!）* 5（*73!6 5 1-, !）* 4（ , 173--）* 5 -./*01# * % -.600*
表 ! 硼校准曲线
"#$%& ! ’ (#%)$*#+),- (.*/&
# !
（!9 ( :;）
- * 1- *- /- 0-
$1 -.-18*1* 1.-0<68 0.*!8!0 1-./!61 *-.30<1 *8.186*$* -.-1<<1 1.--<6- 0./61*0 1-.<6*8 *1.*813 *8./<*8
$-.-18381 1.-!*6! 0.!066 1-.81-! *1.1-30 *8.!*-/ 如果在校准曲线测量和未知样品测定中作使用 相同的空白溶液，可不考虑空白的影响［3］，可忽略 !-的不确定度。 *./.1.* 硼标准储备溶液的相对标准不确定度 试验所用的 1---!9 ( :;硼标准溶液购自国家钢铁 材料测试中心，证书给出的不确定度为 -.!=。按 正态分布考虑，置信水平为 33=，> % !。则硼标准 溶液的标准不确定度为 -.!= ( ! % 1 5 1-, !，相对标 准不确定度为 1 5 1-, ! ( 1--- % 1-, 8。 *./.1.! 配制标准溶液所用玻璃量器的相对标准 不确定度 配制标准溶液所用玻璃量器为：1- :; 分度吸管（? 级）0 支，1-- :; 容量瓶（? 级）8 个。 根据 @@A138 ,*--8［1- ］玻璃量器最大允差，在没有特 别指明的情况下，一般可以认为是均匀分布， #> % !。 玻璃量器的相对标准不确定度的具体数值详见表 *。 合成配制标准溶液的所用玻璃量器的相对标准 不确定度为： 0 5（* 763 5 1-, !）* 4 8 5（0 7<< 5 1-, /）# * % 8.810 5 1-, ! 将 *./.1.1，*./.1.*和 *./.1.!分量合成即为! 的标准不确定度： "BCD（!）% -.600* * 4（1-, 8）* 4（8 7810 5 1-, !）# * % -.600* 将试样溶液在与测定标准溶液相同的条件下， 测得硼的辐射强度，在工作曲线上查出相应的硼质 量浓度。1-次试样溶液中硼的浓度!见表 !。 !的相对标准不确定度：-.600* ( !8.-8-8 % -.-*!<*。 表 0 配制校准曲线所用玻璃量器的相对标准不确定度 "#$%& 0 "1& *&%#+)/& 2+#-3#*3 .-(&*+#)-+4 ,5 6,*7)-8 8%#22 (,-+#)-&* .2&3 )- (#%)$*#+),- (.*/& 玻璃量器 EFB>+G9 9D#HH$FGI#+GCB
最大允差（:;）
JFHI KCB:+IICL LCM+#I+FG
标准不确定度（:;）
&I#GL#BL "G$CBI#+GI’ 相对标准不确定度 NCD#I+MC HI#GL#BL "G$CBI#+GI’
1- :;分度吸管 1- :; 9B#L"#ICL K+KCIIC O -.-0 -.-*63 "1 % *.63 5 1-, !
1-- :;容量瓶 1-- :; MFD":C PD#H> O -.1- -.-0<< "* % 0.<< 5 1-, /
注：标准不确定度为最大允差 #( !，相对标准不确定度为标准不确定度 (标称总容量。
QFIC：&I#GL#BL "G$CBI#+GI’ % JFHI KCB:+IICL LCM+#I #+FG ( !；NCD#I+MC HI#GL#BL "G$CBI#+GI’ % &I#GL#BL "G$CBI#+GI’ ( RGL+$#ICL IFI#D $#K#$+I’7
表 9 试样溶液中硼的浓度（!8 : ;<）
"#$%& 9 ’ (,-(&-+*#+),- )- 2#;=%& 2,%.+),- 浓度 测量次数 SCHI+G9 HCT"CG$C 平均值
UFG$7 1 * ! / 0 8 < 6 3 1- ?MCB#9C ! /-.*11* //.-68 !6.*08 *0.68!! *<.13 /*.08!1 /-.*!80 !3.1*!6 *8.!360 !8.1*0< !8.-8-8<br*/*1 植 物 营 养 与 肥 料 学 报 10卷 !"#"! 提取液总体积 ! 的相对标准不确定度 "#$%（!） 样品经过预处理后定容于 !$% &’容量瓶 中，查 (() *+,-!%%,［*% ］，标称总容量为 !$% &’ 的 .
级容量瓶容量允差为 / %"*$&’，仍按均匀分布，0 !1 2。标准不确定度为 !%"*$ 3 2 1 %"%4,,%，相对标
准不确定度为 %"%4,, 3 !$% 1 2"#,# 5 *%- #。 !"#"2 样品质量 & 的相对标准不确定度 "#$%（&）
& 的不确定度主要由天平带来（将称量中的变动性
归并至重复性中考虑），查天平检定证书得 ’++ 1
%",, &6。按均匀分布，并考虑是两次称量（分别为
去皮质量和样品质量）获得，& 的标准不确定度：
"（&）1 %",, 5 *%
- 2 (
!2
!5 ! 1 $)24+ 5 *%- # ( 相对标准不确定度： "#$%（&）1
"（&）
"& 1
$)24+ 5 *%- # ( * ),$!! ( 1 2 )!,! 5 *%
- #
*%次样品质量的见表 #。
表 ! 样品质量（"）
#$%&’ ! ($)*&’ )$++ 样品质量 测量次数 789:;<6 98=>8?8 平均值<br@A&BC8 &A99 * ! 2 #$ , D 4 + *% .E8FA68
& *"+D4+ *"44%4 *"+!#2 *"%!+, *"2*2$*"D4$# !"%##$*"+$,! *"*%2, *"$%$! *",$!! !"#"# 测量重复性 F8B 的相对标准不确定度 "#$%
（ #$*） 该样品 *%次硼含量（表$）是由独立检验过
程完成的，即每一次测定都是经过样品称量、振荡溶
解、定容、上机测定等步骤来完成。因此，样品均匀
性、检验人员操作的随机性、仪器的漂移均能引起测
定结果的差异。根据 *%次样品中硼含量计算，单次
测量结果的实验标准偏差 +（,）1 %"%$!4!G，平均 值的标准偏差 +（#,）1 +（,） !- 1 +（,） !*% （ - 为测量次 数），平均值的相对标准偏差 ./ 1 +（#,） , 1 %"%*,D%0 %"$$20 1 %"%2%!%，即为相对标准不确定度。 表 , 样品中硼含量 #$%&’ , - ./01’01 20 +$)*&’ 含量（G） 测量次数 789:;<6 98=>8?8 平均值<brHI<:8<: * ! 2 #$ , D 4 + *% .E8FA68
, %"$%4 %"$4, %"#+D %",!4 %"$!4 %"$+, %"#+! %"$%% %"$+4 %",%% %" $$2 !"#"$ 合成标准不确定度 由于各项误差分量彼 此相互独立，总的相对标准不确定度由上述各个分 量合成而得，即： "1#$%（,）1 "!#$%（!）J " ! #$%（!）J "!#$%（&）J "!#$%（ #$*! ） 1 （%"%!2D!）! J（2)#,#5 *%-#）! J（2)!,!5 *%-#）! J（%"%2%!%）! ! 1 %"%24 "1（,）2 "1#$%（,）·#, 2 %"%24 3 %"$$ 2（0）2
%"%!*（0）
!"#", 扩展不确定度 ’ 当置信概率 * 1 ++G，取
包含因子 0 1 !［!，4］，则测量结果的扩展不确定度为：
’ 1 !"1（,）1 ! 5 %"%!*（0）1 %"%#!（0）
硼含量的测量结果可表述为：
,（4）1（%" $$2 / %"%#!）G。 3 结论 KHLM.N@法测定水溶肥料中硼含量的测量结果 不确定度由多个分量所组成，这些分量基本考虑了 测量过程中系统误差和随机误差引起的不确定度， 最后给出测量结果及其不确定度为：O（P）1（%"$$ 2
/ %"%#!）G。根据上面不确定度评定过程可以看
出，标准系列点数、样品测量次数越多，不确定度越
小。虽然在不确定度公式中没有直接反映分析仪
器、方法原理等方面带来的影响，但实际上这些不确
定度已经包括在实验的标准偏差之中。因此，加强
仪器校准和维护，并保证浓度范围符合方法原理的
要求是十分必要的。掌握与不确定度有关的因素，
可以更好地做好实验设计，降低不确定度，提高分析
结果的准确度。
参 考 文 献：
［*］ 国家质量技术监督局 Q )P 3 7*$#4*-!%%%" 检测和校准实验室能
力的通用要求［@］Q 北京：国家标准出版社，!%%* Q
HR;F8A> IS T>AC;:U AB8FE;9;I2#!*$期 范洪黎，等：KHLM.N@法测定水溶肥料中硼含量的不确定度分析 !"###$ %&’&()* (&+,-(&.&’/0 12( /3& 42.5&/&’4& 21 /&0/-’6 )’7 4)*-8
9()/-2’ *)92()/2(-&0［:］; <&-=-’6：:/)’7)(70 >(&00 21 ?3-’)，"##@;
［"］ 中国实验室国家认可委员会 ; 化学分析中不确定度的评估指
南［A］; 北京：中国计量出版社，"##"$B!C，B#!BD$
?3-’) E)/-2’)* F44(&7-/)/-2’ 12( G)92()/2(H; %,-7& /2 /3& &I)*,)/-2’ 21
,’4&(/)-’/H -’ 43&.-4)* )’)*H0-0［A］; <&-=-’6：?3-’) A&/(2*26H >,98
*-03-’6 J2,0&，"##"$B!C，B#!BD$
［K］ 国家质量监督检验检疫总局 ; 测量不确定度评定与表示［A］;
北京：中国计量出版社，@DDD$CK!CL$
%&’&()* F7.-’-0/()/-2’ 21 M,)*-/H :,5&(I-0-2’，-’05&4/-2’ M,)()’/-’& 21
/3& >&25*&’0 N&5,9*-4 21 ?3-’); OI)*,)/-2’ )’7 &P5(&00-2’ 21 ,’4&(8
/)-’/H -’ .&)0,(&.&’/［A］; <&-=-’6：?3-’) A&/(2*26H >,9*-03-’6
J2,0&，@DDD$CK!CL$
［B］ 施昌彦 ; 测量不确定度评定与表示指南［A］; 北京：中国计量
出版社，"###$"B!LL$
:3- ? Q; %,-7& /2 &I)*,)/-2’ )’7 &P5(&00-2’ 21 ,’4&(/)-’/H -’ .&)0,(&8
.&’/［A］; <&-=-’6：?3-’) A&/(2*26H >,9*-03-’6 J2,0&，"###$"B!LL$
［C］ R:S T RO? @U#"C："##C$%&’&()* (&+,-(&.&’/0 12( /3& 42.5&/&’4& 21 4)*-9()/-2’ )’7 /&0/-’6 *)92()/2(-&0［:］; %&’&I)：R:S，"##C$ B!@L$［L］ EQ @B"V!"##U$ 微量元素水溶肥料［:］; 北京：中华人民共和国
农业部，"##U$@K!@C$
EQ @B"V!"##U$W)/&(802*,9*& 1&(/-*-X&(0 42’/)-’-’6 .-4(2’,/(-&’/0 ［:］; <&-=-’6：A-’-0/(H 21 F6(-4,*/,(& 21 /3& >&25*&’0 N&5,9*-4 21 ?3-8 ’)，"##U$ @K!@C$［U］ 陈德文，张卫中，潘思轶 ; 原子荧光光谱法测定米发糕中砷的 测量不确定度的评定［Y］; 农产品加工学刊，"##U，@@C（@#）：U" !UB$
?3&’ Z W，[3)’6 W [，>)’ : Q; OI)*,)/-2’ 2’ /3& ,’4&(/)-’/H 21
)(0&’-4 -’ (-4& 4)\& 9H )/2.-4 1*,2(&04&’4& 05&4/(2.&/(H［Y］; F4)7;
>&(-27; ])(. >(27; >(24;，"##U，@@C（@#）：U"!UB$［V］ 徐为霞，郭智广，王亮，罗俊霞 ; 火焰法—原子吸收测定蔬菜中 铜的测量不确定度分析［Y］; 微量元素与健康研究，"##V，"C （@）：B"!BK$
^, W ^，%,2 [ %，W)’6 G，G,2 Y ^; F’)*H0-0 21 ,’4&(/)-’/H -’
.&)0,(&.&’/ 12( 7&/&(.-’)/-2’ 21 4255&( 9H 1*).&8FF: -’ I&6&/)9*&0
［Y］; :/,7; _()4& O*&.; J&)*/3，"##V，"C（@）：B"!BK$［D］ 王彦刚 ; 化学实验室仪器分析法不确定度的计算［Y］; 计量技 术，"##K，（@@）：CL!CV$
W)’6 Q %; ?)*4,*)/-2’ 21 ,’4&(/)-’/H -’ -’0/(,.&’/)* )’)*H0-0 -’
43&.-4)* *)92()/2(H［Y］; A&/(2* ; _&43;，"##K，（@@）：CL!CV$［@#］ YY% @DL!"##L$ 常用玻璃量器检定规程［:］; 北京：国家质量监
督检验检疫总局，"##L$L!V$
YY% @DL!"##L$‘&(-1-4)/-2’ (&6,*)/-2’ 21 a2(\-’6 6*)00 42’/)-’&(［:］; <&-=-’6：%&’&()* F7.-’-0/()/-2’ 21 M,)*-/H :,5&(I-0-2’， R’05&4/-2’ M,)()’/-’& 21 /3& >&25*&’0 N&5,9*-4 21 ?3-’)，"##L$ L!V$
BB"@ 植 物 营 养 与 肥 料 学 报 @C卷