鄭永輝
(永春縣環境監測站,福建 泉州 362600)
在不確定度評估方法中,bottom-up和top-down是各種標準或指南中主要應用的方法,但大多以bottom-up方法對測量不確定度進行評定。該方法對于所有的物力測量都可以,其在測量中更加注重細節,但存在計算和操作步驟繁瑣的問題,降低了評估測量的精準度,而且,在對較為復雜的化學測量時,存在測量局限性;同時,對一些不確定度分量是很難確保量化的精準性。2018年發布實施的《基于質控數據環境檢測測量不確定度評定指南》(CNAS-GL022)[1],是GUM的簡化和延伸。相比于GUM,其適用性更加廣泛、技術操作性更強,它為環境監測實驗室開展不確定度評定提供了重要的參考。 該指南采用國際上較為流行的四種方法,即Top-Down技術,主要有精密度法、線性擬合法、控制圖法,以及經驗模型法。長期累計的質控數據,是該方法的技術手段,可以實現對評估過程中不確定度的潛在來源的全面反應,加大了精準度率。應用top-down方法開展化學測量結果不確定度的評估,比bottom-up方法更加具有實用性。本文采用《基于質控數據環境檢測測量不確定度評定指南》(CNAS-GL022)中線性擬合法,以實驗室日常質控數據作為信息支持,開展了水中總磷測定結果的不確定度的評定。
以《水質 總磷的測定 鉬酸銨分光光度法》(GB 11893-89)[2]作為實驗室主要技術依據,使用環境保護標準樣品研究所的質量濃度為500 mg/L的有證標準物質,經稀釋得到2.00 mg/L的磷標準使用液;再經稀釋,獲得6個質量濃度水平的標準溶液。具體為:0.04 mg/L,0.08 mg/L,0.24 mg/L,0.40 mg/L,0.80 mg/L,1.20 mg/L。需要進行4次重復測量工作,以保障6個質量濃度水平的精準度。以一個月的時間段作為間隔時間,保障測量的獨立性。表1給出了擬合實驗數據的詳細匯總。

表1 擬合實驗數據匯總

表1(續)
1)假設標準偏差為常數模型。表1中數據組按公式(1)進行擬合工作曲線。
yn=β0+β1RQV
(1)
式中:yn為ynk的估計值,ynk為第n個水平(n=1,……,N)的第k次測量值(k=1,……,K);β0為截距估計值;β1為斜率估計值。
2)由公式(1)計算可得常數擬合工作曲線為:
yn=β0+β1RQV=0.0004+0.9991RQV
3)使用表1所列RQV值代入常數擬合工作曲線中,求得擬合值yn,列于表2。
4)殘差值為每次測量值ynk與擬合值yn之差,即:
enk=ynk-yn
(2)
殘差值計算結果列于表2。

表2 常數模型殘差值 mg/L
5)將殘差對應于擬合值,繪制出相應的分析圖進行數據分析(見圖1)。觀察圖1中信息發現,各點的分布情況主要是以0點為中心隨機分布。由此可知,沒有偏離線性的假定,常數模型的擬合假定成立。

圖1 常數模型假定下的殘差作圖
模型誤差計算如式(3)、(4)。
(3)
(4)

根據公式(3)和公式(4),yn=β0+β1RQV=0.0004+0.9991RQV的方差分析結果見表3。
根據表3中的統計:F=0.976 表3 失擬誤差與試驗誤差的比較 控制上限(UCL)和控制下限(LCL)的計算如(5)、(6): (5) (6) 根據公式(5)和公式(6),計算得: 變換值的計算如式(7)、(8)。 (7) (8) 實驗室在《水質 總磷的測定 鉬酸銨分光光度法》(GB 11893-89)實施期間,對于測量的范圍區間的覆蓋情況,可以在所用的6個質量濃度的標準溶液中,選擇最低和最高的兩個水平(0.04 mg/L和 1.20 mg/L)開展,每天需要進行重復雙試實驗。實驗的前提條件為,確保在期間精密度條件下進行,連續7天。表4給出了一周的測量數據。 表4 后續測量值的變換與監控 根據表4計算結果,按照時間順序建立相應的控制圖(見圖2)。圖2中顯示,di值均落在控制限內,表明yn=β0+β1RQV=0.0004+0.9991RQV處于統計受控狀態。 圖2 yn=β0+β1RQV=0.0004+0.9991RQV的控制圖 6.1 標準不確定度 根據控制值di,按下列公式計算出標準不確定度。 (9) 式中:dlj和dmj分別為最低和最高標準物質/標準樣品的di值,其中,j為測量時間點;J為測量次數。 根據表4數據,計算得:SR=0.006 mg/L 在0.04~1.20 mg/L范圍內,在包含概率95%的概率下,方法擴展不確定度U=0.006×2=0.012 mg/L。 基于實驗室長期積累的質控數據來評定測量不確定度的方法,具有簡便與操作性強等特點。相對GUM評定法結果只能對應一個測量值,該方法能夠實現不確定度在一定濃度范圍內的實時評定,具有更廣泛的應用場景。同時,實驗室在開展質控數據檢測的時候,將top-down方法與測量不確定評估有效結合起來,能夠滿足“頂層設計”的實際需求。通過線性擬合法評定水中總磷測定的不確定度時,測定結果在校準曲線0.04~1.20 mg/L 范圍內,擴展不確定度為 0.012 mg/L。
4 控制限的確定
5 后續測量值的變換



6 不確定度的評定
6.2 擴展不確定度
7 結語