用質(zhì)控圖和穩(wěn)健統(tǒng)計-迭代法評估環(huán)境檢測實驗室測量不確定度

狄一安，孫海容，孫培琴，任立軍，劉 巖，周 昊，王婧瑞，李斯明，李玉武*

(1.國家環(huán)境分析測試中心，北京 100029； 2.中國合格評定國家認(rèn)可委員會，北京 100062)

近年來，在中國合格評定國家認(rèn)可委員會(CNAS)的大力推動下，測量結(jié)果不確定度評估越來越受到環(huán)境領(lǐng)域各級檢測實驗室的重視。CNAS要求已獲認(rèn)可的檢測實驗室應(yīng)有能力對每一項有數(shù)值要求的測量結(jié)果進(jìn)行測量不確定度評估。如果客戶有要求，檢測報告必須提供測量結(jié)果的不確定度。

目前各種標(biāo)準(zhǔn)或指南中測量不確定度的評估方法大致可分為兩類：“bottom-up”和“top-down”。“bottom-up”法注重細(xì)節(jié)，分析、計算分析步驟中每一操作環(huán)節(jié)所涉及的不確定度分量，計算合成標(biāo)準(zhǔn)不確定度，然后乘以包含因子得到擴展不確定度。此法可完全覆蓋物理測量，但對與化學(xué)和生物測量相關(guān)的特殊問題(例如關(guān)于樣品前處理步驟)沒有特定的規(guī)定。雖然對分析程序中的單一步驟(稱量、溶液移取等)，標(biāo)準(zhǔn)溶液濃度的計算和校準(zhǔn)曲線擬合不確定度的評估，《化學(xué)分析中不確定度的評估指南》(GUM)[1]都給予一般原則和指導(dǎo)，但應(yīng)用于較為復(fù)雜的化學(xué)測量時仍有一定局限性[2]?！皌op-down”法利用從方法確認(rèn)、實驗室內(nèi)質(zhì)控和實驗室間協(xié)作定值、能力驗證等數(shù)據(jù)，注重從整體上、通過數(shù)月、數(shù)年等一段期間反映樣品檢測全過程的精密度數(shù)據(jù)直接評估測量不確定度，如精密度法[3]、質(zhì)控圖法[4]、線性校準(zhǔn)法[5]和經(jīng)驗?zāi)Ｐ头╗6-7]等。由于使用了長期積累的質(zhì)控數(shù)據(jù),評估過程全面反映不確定度的潛在來源的機率會更大。在評估化學(xué)測量結(jié)果的不確定度時，“top-down”法比“bottom-up”法更為實用。“top-down”法充分利用了實驗室質(zhì)控數(shù)據(jù)并將檢測方法測量不確定度評估與實驗室質(zhì)控工作緊密聯(lián)系在一起，促使實驗室在進(jìn)行質(zhì)控方案“頂層設(shè)計”時就能結(jié)合測量不確定度評估的需要，將測量不確定度評估基礎(chǔ)工作抓實，其意義不可低估。近年來我國已有不少應(yīng)用方法精密度、實驗室能力驗證結(jié)果來評估測量不確定度的報道[8-11]，這些研究工作將測量不確定度與分析方法的質(zhì)量參數(shù)聯(lián)系起來，是“top-down”方法的成功探索。

2011年6月6～7日歐洲化學(xué)聯(lián)合會在葡萄牙里斯本召開了測量不確定度最新進(jìn)展研討會[12]，兩篇特邀報告[13-14]集中反映了測量不確定度評估方法指南第3版結(jié)構(gòu)和內(nèi)容修訂動向，其中最大的變化是增加了“top-down”方法。不少國際組織也發(fā)表相關(guān)技術(shù)文件[2,15-17]。由此可見，在檢測實驗室推廣普及“top-down”理念，將“top-down”法與經(jīng)典的“bottom-up”法相結(jié)合，在國內(nèi)外均是大勢所趨。

基于“top-down”技術(shù)理念的計算方法中，質(zhì)控圖法和線性校準(zhǔn)法主要用于實驗室內(nèi)實驗數(shù)據(jù)處理。質(zhì)控圖法操作簡便，但只涉及一個濃度，并且要求所有數(shù)據(jù)是非離群測量結(jié)果。線性校準(zhǔn)法可用于不同濃度樣品測量結(jié)果不確定度評估，但在建立模型時，要求不同濃度的標(biāo)準(zhǔn)樣品測量次數(shù)一致，需要完整、成套不同濃度水平標(biāo)準(zhǔn)樣品測試數(shù)據(jù)。實驗室以前積累的質(zhì)控數(shù)據(jù)很難滿足這些要求，這使其推廣應(yīng)用受到一定限制。本文嘗試用回收率將同一檢測項目不同濃度水平質(zhì)控數(shù)據(jù)歸一化，然后用質(zhì)控圖方法處理，如果存在離群數(shù)據(jù)時，用穩(wěn)健統(tǒng)計法-迭代法計算期間精密度sR′。探討了下列具體問題：①單一濃度數(shù)據(jù)歸一化前后，同一種計算方法結(jié)果比較；②同一套單一濃度質(zhì)控數(shù)據(jù)未經(jīng)任何前處理，穩(wěn)健統(tǒng)計-迭代法與質(zhì)控圖法結(jié)果比較；③多濃度水平數(shù)據(jù)歸一化條件下，穩(wěn)健統(tǒng)計-迭代法與質(zhì)控圖法結(jié)果比較；④不同濃度水平數(shù)據(jù)歸一化后，質(zhì)控圖法、穩(wěn)健統(tǒng)計-迭代法與線性校準(zhǔn)法計算結(jié)果是否相符。采用作者所在實驗室積累的5套和其他實驗室提供的19套環(huán)境檢測實驗室常規(guī)項目質(zhì)控數(shù)據(jù)驗證了該方法的可行性。

1 實驗數(shù)據(jù)來源

氨氮、總磷、COD，土壤中鉛、土壤中銅5套數(shù)據(jù)來源于本文作者所在實驗室2010～2012年期間積累的質(zhì)控標(biāo)準(zhǔn)樣品實測數(shù)據(jù)(參見表1)，標(biāo)準(zhǔn)樣品均由環(huán)保部標(biāo)準(zhǔn)樣品研究所提供。19套其他環(huán)境檢測實驗室實測數(shù)據(jù)，來源于CNAS組織的用“top-down”技術(shù)評估不確定度培訓(xùn)班教材和學(xué)員提交的報告，參見表2。

表1 作者實驗室不同檢測項目質(zhì)控樣品標(biāo)稱值及測定結(jié)果

表2 其他實驗室提供的質(zhì)控數(shù)據(jù)

2 不確定度評估方法

2.1 質(zhì)控圖法[4]

根據(jù)相關(guān)表格提供的數(shù)據(jù)，可將wi值換算成正態(tài)概率值pi。更方便的方法是利用微軟Excel的函數(shù)公式NORMDIST(x,mean,std dev,cummulative)語句求出。例如，wi=-0.64，NORMDIST(-0.64，0，1，TRUE)=0.261086，與查表值0.2611完全一致。A值和A*值按下列公式計算：

式中：A*—正態(tài)統(tǒng)計量，A的修正值。按s式計算時表示A*(s)，按MR式計算時表示A*(MR)；n—測量次數(shù)。

根據(jù)A*(s)和A*(MR)數(shù)值，可作如下判定。

(a)A*(s)<1.0和A*(MR)<1.0，接受數(shù)據(jù)的正態(tài)性和獨立性的假定。

(b)A*(s)>1.0和A*(MR)>1.0，表明測量系統(tǒng)失控。

(c)A*(s)<1.0和A*(MR)>1.0，表明系列結(jié)果呈非獨立性。

在偏倚受控的期間精密度測量條件下，2×sR′可作為不確定度的評估值。實驗室給出的sR′應(yīng)大于sr，小于sR。

表3以表1中COD歸一化數(shù)據(jù)為例，顯示了整個計算過程及結(jié)果。表3結(jié)果表明，COD歸一化數(shù)據(jù)質(zhì)控圖的計算結(jié)果：平均值為1.000，期間精密度為0.023。正態(tài)性和獨立性均通過AD檢驗。COD回收率及其擴展不確定度為1.000±0.046。

表3 COD的質(zhì)控圖法計算示例

2.2 線性校準(zhǔn)法及數(shù)據(jù)歸一化

線性校準(zhǔn)法[5]原理和主要計算步驟如下：通過帶基體不同濃度水平的標(biāo)準(zhǔn)物質(zhì)測定值與標(biāo)準(zhǔn)偏差建立回歸模型，根據(jù)標(biāo)準(zhǔn)樣品濃度與重復(fù)測量值標(biāo)準(zhǔn)偏差是否顯著相關(guān)確定其屬于比例模型還是常數(shù)模型。若屬于前者應(yīng)進(jìn)行加權(quán)變換。利用單因素的方差分析，檢查所擬合模型的偏倚是否處于統(tǒng)計受控。在期間精密度測量條件下，利用較高和較低兩個標(biāo)準(zhǔn)物質(zhì)，通過相應(yīng)的數(shù)據(jù)變換給出殘差數(shù)據(jù)，并建立質(zhì)控圖。如果圖中的系列數(shù)據(jù)分布呈隨機狀態(tài)，則測量系統(tǒng)給出均方根誤差，即為該實驗室獲得的不確定度估計值。比例模型計算公式如下：

(1)

為處理不同濃度質(zhì)控樣品測量數(shù)據(jù)，可將測試結(jié)果歸一化。期間精密度(測試結(jié)果與標(biāo)準(zhǔn)值相對標(biāo)準(zhǔn)偏差的均方根)計算公式(2)與公式(1)完全相同，與計算標(biāo)準(zhǔn)偏差的公式也基本相同：

(2)

式中，m為不同濃度水平樣品個數(shù)，當(dāng)回收率平均值約等于1時，式(1)和式(2)基本相同。歸一化法計算結(jié)果是所有不同濃度標(biāo)準(zhǔn)樣品測量值與標(biāo)準(zhǔn)值相對標(biāo)準(zhǔn)偏差的均方根。它與線性校準(zhǔn)法計算結(jié)果的區(qū)別是，后者是模型建立后一段期間同時測定高低兩個濃度標(biāo)準(zhǔn)樣品測量值與標(biāo)準(zhǔn)值相對標(biāo)準(zhǔn)偏差的均方根。

公式(2)是計算歸一化后質(zhì)控數(shù)據(jù)sR′的基礎(chǔ)。但實際計算時可采用下列兩條途徑：①如果不存在“離群值”，可用質(zhì)控圖法處理，同時得到sR′和數(shù)據(jù)正態(tài)性和獨立性判據(jù)；②用穩(wěn)健統(tǒng)計-迭代法處理，得到歸一化數(shù)據(jù)的穩(wěn)健性標(biāo)準(zhǔn)偏差sR′。其優(yōu)點是它不需要對離群值進(jìn)行識別和剔除，可將離群值對統(tǒng)計結(jié)果平均值和標(biāo)準(zhǔn)偏差的影響降低到最小。

2.3 穩(wěn)健統(tǒng)計-迭代法

穩(wěn)健統(tǒng)計是指不用事先識別、剔除離群值，直接應(yīng)用全部測量數(shù)據(jù)，將離群值對統(tǒng)計分析結(jié)果影響降低到最小的統(tǒng)計分析方法。迭代法是近年來ISO標(biāo)準(zhǔn)推薦的穩(wěn)健統(tǒng)計方法[18-19]，其應(yīng)用早在20世紀(jì)80年代國外分析化學(xué)文獻(xiàn)中已有報道[20]。它不需要對“離群值”進(jìn)行人為干預(yù)。其原理是對按順序排列，位于數(shù)據(jù)排列兩端遠(yuǎn)離中位值的“可疑值”或“離群值”均以較小權(quán)重予以保留，與中位值接近的值則以較大權(quán)重參與計算，充分利用了全部測量數(shù)據(jù)的信息。分析化學(xué)領(lǐng)域中，有些離群值很明顯，容易刪去，但有時判斷“離群值”的分界線并不明顯，不同的判別方法結(jié)論可能不同。在這種情況下，穩(wěn)健統(tǒng)計方法盡量減小了“離群值”對“平均值”和“標(biāo)準(zhǔn)偏差”的影響，給出這兩個統(tǒng)計量的估計值。穩(wěn)健統(tǒng)計-迭代法的特點是它不需要假設(shè)數(shù)據(jù)是標(biāo)準(zhǔn)正態(tài)分布。在正態(tài)分布基礎(chǔ)上，單峰分布，基本對稱，存在離群值，即使有嚴(yán)重拖尾(heavy tails)現(xiàn)象時，它也能給出較合理的結(jié)果[20]。國內(nèi)能力驗證領(lǐng)域有相關(guān)文獻(xiàn)報道[21-22]，但在實驗室測量數(shù)據(jù)處理方面應(yīng)用報道較少。穩(wěn)健統(tǒng)計-迭代法具體計算步驟如下。

sR′=1.134×s

s*=1.5×sR′

表4 穩(wěn)健統(tǒng)計-迭代法(方法1)計算示例

初始x1和x2也可以按下列方式得到(方法2，參見表5)：測量結(jié)果xi按升序排列后，取數(shù)列的中位值。用排序后的測量數(shù)據(jù)與中位值相減，求出此差值的絕對值數(shù)列的中位值作為標(biāo)準(zhǔn)偏差的估計值s0。以此數(shù)據(jù)為基礎(chǔ)計算下列參數(shù)：

sR′=1.483×s0

s*=1.5×sR′

表5 穩(wěn)健統(tǒng)計-迭代法(方法2)計算示例

表4和表5的計算示例結(jié)果表明，基于兩種不同方法得到的初始參數(shù)進(jìn)行計算的穩(wěn)健統(tǒng)計-迭代法均得到相同計算結(jié)果：平均值為1.000，期間精密度sR′為0.022。COD回收率及其擴展不確定度為1.000±0.044。本示例中，COD數(shù)據(jù)迭代法與質(zhì)控圖法結(jié)果(1.000±0.046)基本相符。質(zhì)控圖法和穩(wěn)健統(tǒng)計-迭代法計算均在Excel軟件上完成，計算示例文件可向本文通訊作者索取。

3 結(jié)果與討論

3.1 “top-down”方法特點、現(xiàn)存問題及解決方案

經(jīng)典的評估測量不確定度方法僅適用于測量過程模型已知的情況，通過分析不確定度來源并對系列觀測值的統(tǒng)計分析(A類評估)或其他方式(B類評估)來完成。方法步驟繁瑣，容易造成不確定度分量遺漏或重復(fù)計算，某些不確定度分量難以準(zhǔn)確量化。大量研究結(jié)果表明，化學(xué)分析中不確定度的主要分量大多為校準(zhǔn)曲線和樣品重復(fù)測定。但這種重復(fù)性測定一般均是當(dāng)天或短時間內(nèi)完成，未能反映實驗室長期質(zhì)控狀態(tài)?！癰ottom-up”法適合新方法建立或驗證期間，實驗室結(jié)合工作曲線、準(zhǔn)確度(回收率)和精密度等實驗數(shù)據(jù)，在短期內(nèi)完成對分析方法測量不確定度的評估。

用于評估測量不確定度的“top-down”技術(shù)核心理念是兩個統(tǒng)計量“期間精密度sR′”和“再現(xiàn)性標(biāo)準(zhǔn)偏差sR”。兩者分別以實驗室內(nèi)部測量數(shù)據(jù)和不同實驗室協(xié)作實驗測試數(shù)據(jù)為基礎(chǔ)。利用這兩個統(tǒng)計量可以分別評估實驗室內(nèi)測量不確定度和分析方法不確定度。實驗室內(nèi)測量不確定度的大小反映了實驗室技術(shù)水平，分析方法或樣品測量不確定度可用于對同一樣品不同實驗室測定結(jié)果的比較。

基于以上分類原理，質(zhì)控圖法、線性校準(zhǔn)法以實驗室內(nèi)部質(zhì)控數(shù)據(jù)為基礎(chǔ)，可計算實驗室測量不確定度。精密度法、經(jīng)驗?zāi)Ｐ头ㄒ阅芰︱炞C數(shù)據(jù)等不同實驗室統(tǒng)計數(shù)據(jù)為基礎(chǔ)，則可以評估分析方法(或樣品不確定度)不確定度。能力驗證數(shù)據(jù)統(tǒng)計結(jié)果“穩(wěn)健標(biāo)準(zhǔn)偏差”是判斷實驗室上報數(shù)據(jù)是否合格的尺度，同時也是方法不確定度的有效估計值。常見的能力驗證數(shù)據(jù)穩(wěn)健統(tǒng)計方法有四分位法和迭代法。四分位數(shù)法假設(shè)數(shù)據(jù)分布是標(biāo)準(zhǔn)正態(tài)分布，且只考慮了50%的上報數(shù)據(jù)。如果實際上報數(shù)據(jù)明顯違背分布假設(shè)，此時直接套用計算公式，其計算結(jié)果穩(wěn)健標(biāo)準(zhǔn)偏差有時會過于“嚴(yán)格”[18]。迭代法也是一種穩(wěn)健統(tǒng)計方法，已得到國際標(biāo)準(zhǔn)認(rèn)可和推薦[19-20]。迭代法對于離群值采取權(quán)重小、位于中位附近數(shù)據(jù)權(quán)重大的方式進(jìn)行計算，穩(wěn)健性強，已受到國內(nèi)專業(yè)機構(gòu)關(guān)注[21-23]。但國內(nèi)文獻(xiàn)中迭代法在處理實驗室測量數(shù)據(jù)應(yīng)用很少有報道。

線性校準(zhǔn)法可以用于評估不同濃度范圍樣品的測量不確定度，但在建立模型時對不同濃度樣品測量次數(shù)嚴(yán)格要求一致，需要成套、完整的實驗數(shù)據(jù)，模型檢驗理論性強，初學(xué)者不易掌握。質(zhì)控圖法雖然只能處理一個濃度，但它可提供數(shù)據(jù)正態(tài)性和獨立性判據(jù)，用于判斷計算結(jié)果可靠性和質(zhì)控數(shù)據(jù)的有效性，但嚴(yán)格要求測量數(shù)據(jù)中不能含離群值。如果將不同濃度數(shù)據(jù)利用回收率進(jìn)行歸一化，這樣就可以同時利用質(zhì)控圖和穩(wěn)健統(tǒng)計-迭代法來處理實驗室內(nèi)長期積累的不同濃度質(zhì)控數(shù)據(jù)，得到sR′數(shù)據(jù)，然后乘以2，可得到實驗室擴展測量不確定度。

3.2 單一濃度數(shù)據(jù)歸一化前后質(zhì)控圖法與穩(wěn)健統(tǒng)計-迭代法結(jié)果比較

為檢驗數(shù)據(jù)歸一化對計算結(jié)果的影響，作者首先考察了單一濃度歸一化前后同一種計算方法計算結(jié)果有無變化。從理論上分析，歸一化對結(jié)果應(yīng)該無影響。用表2的實測數(shù)據(jù)進(jìn)行了驗證。表6顯示了單一濃度測試數(shù)據(jù)歸一化前后質(zhì)控圖法和穩(wěn)健統(tǒng)計-迭代法計算結(jié)果，結(jié)果分別用絕對值和回收率表示。表6(第1～4套數(shù)據(jù))計算結(jié)果表明：①與傳統(tǒng)的質(zhì)控圖法相比，測量數(shù)據(jù)歸一化后，質(zhì)控圖法與歸一化前原結(jié)果完全相符。數(shù)據(jù)歸一化前后穩(wěn)健統(tǒng)計-迭代法計算結(jié)果也完全相符；同套數(shù)據(jù)歸一化后，質(zhì)控圖法和穩(wěn)健統(tǒng)計-迭代法計算結(jié)果sR′(相對值)平均偏差為0.46%，兩種方法計算結(jié)果基本相符；②第3套數(shù)據(jù)苯并[a]芘第8～10行數(shù)據(jù)是不同時間段數(shù)據(jù)合并后計算結(jié)果[24]。第8行是sR，合并計算結(jié)果，第9～10行是測量數(shù)據(jù)合并后計算結(jié)果。不同方法的計算結(jié)果完全一致，但新方法避免了合并計算的繁瑣公式。表7是單一濃度質(zhì)控數(shù)據(jù)，不經(jīng)任何處理，質(zhì)控圖法和計算結(jié)果比較。從表7(第5～11套數(shù)據(jù))計算結(jié)果可以看出，穩(wěn)健統(tǒng)計-迭代法與質(zhì)控圖法結(jié)果基本相符，sR′(相對值)平均偏差為0.15%。表6和表7用于方法比較的11套質(zhì)控數(shù)據(jù)均通過正態(tài)性和獨立性檢驗。這為計算結(jié)果的可靠性提供了基本保障。這是質(zhì)控圖區(qū)別于其他方法的顯著特點。

表6 單一濃度數(shù)據(jù)歸一化前后質(zhì)控圖法及穩(wěn)健統(tǒng)計-迭代法結(jié)果比較

表7 單一濃度數(shù)據(jù)質(zhì)控圖法及穩(wěn)健統(tǒng)計-迭代法結(jié)果比較

3.3 不同濃度水平數(shù)據(jù)歸一化后質(zhì)控圖與穩(wěn)健統(tǒng)計-迭代法結(jié)果比較

實驗室在進(jìn)行質(zhì)控時，質(zhì)控標(biāo)準(zhǔn)樣品濃度需要根據(jù)樣品的大致濃度進(jìn)行選擇。有時一種濃度標(biāo)準(zhǔn)樣品用完了，再次采購時，標(biāo)準(zhǔn)樣品濃度也不會完全與以前購買的相同。因此，實驗室積累了不同濃度標(biāo)準(zhǔn)樣品測量數(shù)據(jù)是普遍現(xiàn)象。不同濃度數(shù)據(jù)數(shù)量不一樣，濃度水平分布不均勻。缺乏完整、成套質(zhì)控數(shù)據(jù)很難采用線性校準(zhǔn)法。充分利用實驗室長期積累的質(zhì)控數(shù)據(jù)評估實驗室測量不確定度是一個非常現(xiàn)實的問題。

本文未涉及偏倚檢驗。日常工作積累的質(zhì)控數(shù)據(jù)已經(jīng)和標(biāo)準(zhǔn)樣品值進(jìn)行了比較，符合質(zhì)控要求。

3.4 不同濃度水平歸一化數(shù)據(jù)質(zhì)控圖、穩(wěn)健統(tǒng)計-迭代法與線性校準(zhǔn)法結(jié)果比較

為了考察新方法與傳統(tǒng)方法線性校準(zhǔn)法(比例模型)計算結(jié)果的區(qū)別，利用不確定度培訓(xùn)班學(xué)員為應(yīng)用線性校準(zhǔn)法使用的實驗室質(zhì)控數(shù)據(jù)進(jìn)行了驗證。表9列出的第17～24套數(shù)據(jù)的特點是這些數(shù)據(jù)已有線性校準(zhǔn)法(比例模型)計算結(jié)果。將測量數(shù)據(jù)經(jīng)歸一化處理后，這些數(shù)據(jù)分別采用質(zhì)控圖法和穩(wěn)健統(tǒng)計-迭代法進(jìn)行了計算，結(jié)果用回收率和相對不確定度表示。表9計算結(jié)果表明，質(zhì)控圖法、穩(wěn)健統(tǒng)計-迭代法與線性校準(zhǔn)法結(jié)果sR′平均偏差分別為0.43%和0.20%，質(zhì)控圖法與穩(wěn)健統(tǒng)計-迭代法結(jié)果sR′平均偏差為0.26%。三種計算方法結(jié)果基本相符，穩(wěn)健統(tǒng)計-迭代法與線性校準(zhǔn)法計算結(jié)果更接近。需要指出的是，參與質(zhì)控圖法和穩(wěn)健統(tǒng)計-迭代法計算的是所有質(zhì)控數(shù)據(jù)測量值與標(biāo)準(zhǔn)值相對標(biāo)準(zhǔn)偏差的均方根(見公式2)。它與線性校準(zhǔn)法的結(jié)果的區(qū)別是，后者是模型建立后的質(zhì)控數(shù)據(jù)測量值與標(biāo)準(zhǔn)值相對標(biāo)準(zhǔn)偏差的均方根。從數(shù)據(jù)的數(shù)量上來講，質(zhì)控圖和穩(wěn)健統(tǒng)計-迭代法利用了全部質(zhì)控數(shù)據(jù)。從表9還可以看到，這8套數(shù)據(jù)中有兩套數(shù)據(jù)用質(zhì)控圖法計算時，沒有通過正態(tài)性和獨立性檢驗，這些數(shù)據(jù)可能在短期內(nèi)完成，不是長期積累的質(zhì)控數(shù)據(jù)。注意到質(zhì)控圖計算表中顯示可疑數(shù)據(jù)的個數(shù)(可疑數(shù)據(jù)判別標(biāo)準(zhǔn)參見表3第5列,2

表8 不同濃度水平質(zhì)控數(shù)據(jù)歸一化后質(zhì)控圖與穩(wěn)健統(tǒng)計-迭代法結(jié)果比較

表9 不同濃度水平質(zhì)控數(shù)據(jù)歸一化后三種計算方法結(jié)果比較

本文建立的方法在用于實際樣品測量不確定度評估時，需要注意樣品濃度是否在質(zhì)控數(shù)據(jù)涵蓋的范圍內(nèi)。對于濃度較低，其標(biāo)準(zhǔn)偏差與濃度無關(guān)，在線性校準(zhǔn)法中屬于常數(shù)模型的質(zhì)控數(shù)據(jù)，不能用本文提出的方法處理，其計算結(jié)果可用測定值與標(biāo)準(zhǔn)值的絕對殘差的均方根來表示，相關(guān)研究工作仍在進(jìn)行中。

4 結(jié)語

利用回收率概念，將實驗室長期積累的不同濃度(或含量)水平的水質(zhì)標(biāo)準(zhǔn)樣品和土壤標(biāo)準(zhǔn)樣品測量結(jié)果等質(zhì)控數(shù)據(jù)進(jìn)行歸一化處理，然后用質(zhì)控圖法、穩(wěn)健統(tǒng)計-迭代法計算期間精密度sR′。用本實驗室積累的5套質(zhì)控數(shù)據(jù)和其他實驗室提供的19套數(shù)據(jù)驗證了新方法的可行性。驗證結(jié)果表明，對單一濃度質(zhì)控數(shù)據(jù)，不經(jīng)任何處理，穩(wěn)健統(tǒng)計-迭代法可得到與質(zhì)控圖法基本相符結(jié)果，sR′(相對值)平均偏差為0.15%。對于多濃度水平質(zhì)控數(shù)據(jù)，經(jīng)歸一化后，質(zhì)控圖法、穩(wěn)健統(tǒng)計-迭代法與線性校準(zhǔn)法結(jié)果平均偏差分別為0.43%和0.20%，質(zhì)控圖法與穩(wěn)健統(tǒng)計-迭代法結(jié)果平均偏差為0.26%，三種方法計算結(jié)果基本相符，穩(wěn)健統(tǒng)計-迭代法更接近于線性校準(zhǔn)法計算結(jié)果。質(zhì)控圖法中正態(tài)性和獨立性判據(jù)可用于判斷質(zhì)控數(shù)據(jù)的有效性。新方法為充分利用實驗室長期積累的質(zhì)控數(shù)據(jù)評估實驗室測量不確定度提供了有用工具。

本文提出的方法原理簡單，計算步驟簡便，適用于線性校準(zhǔn)法比例模型數(shù)據(jù)的處理，值得環(huán)境檢測實驗室試用、深入研究及不斷完善。

致謝：本文為驗證方法的可行性，選用了CNAS組織的“top-down”不確定度評估培訓(xùn)班學(xué)員實習(xí)報告中部分實驗室質(zhì)控數(shù)據(jù)，在此向提供數(shù)據(jù)的學(xué)員及實驗室表示感謝。

5 參考文獻(xiàn)

[1] 中國實驗室國家認(rèn)可委員會.化學(xué)分析中不確定度的評估指南[M].北京:中國計量出版社,2002.

[2] National Association of Testing Authorities,Australia (NATA). Technical Note 33,Guidelines for Estimating and Reporting Measurement Uncertainty of Chemical Test Results [R].

[3] 遼寧出入境檢驗檢疫局，中國標(biāo)準(zhǔn)化研究院.GB/Z 22553—2010，利用重復(fù)性、再現(xiàn)性和正確度的估計值評估測量不確定度的指南[S].北京:中國標(biāo)準(zhǔn)出版社,2010.

[4] 遼寧出入境檢驗檢疫局,中國合格評定國家認(rèn)可委員會.GB/T 27407—2010，實驗室質(zhì)量控制-利用統(tǒng)計質(zhì)量保證和控制圖技術(shù)評價分析測量系統(tǒng)的性能[S].北京:中國標(biāo)準(zhǔn)出版社,2010.

[5] 遼寧出入境檢驗檢疫局,中國合格評定國家認(rèn)可中心.GB/T 22554—2010，基于標(biāo)準(zhǔn)樣品的線性校準(zhǔn)[S].北京:中國標(biāo)準(zhǔn)出版社,2010.

[6] 遼寧出入境檢驗檢疫局，中國合格評定國家認(rèn)可中心，中國計量科學(xué)研究院.GB/T 27411—2012，檢測實驗室中常用不確定度評定方法與表示[S].北京:中國標(biāo)準(zhǔn)出版社,2013.

[7] 李玉武，狄一安，孫海容，曹實，任立軍，周昊，郭婧，楊勇杰，王斗文.用經(jīng)驗?zāi)Ｐ驮u估環(huán)境樣品測量不確定度[J].中國無機分析化學(xué)，2012,2(1)：1-8.

[8] 熊英，郭巨權(quán).利用重復(fù)性和再現(xiàn)性的估計值評估測量不確定度[J].巖礦測試，2012，31(2)：350-354.

[9] 李宣，慕俊澤，張斌，江麗媛.塑料中鎘的測量不確定度評定[J].分析試驗室，2007，26(2)：71-74.

[10] 羅惠明，任春華，梁希楊，蔡穎.應(yīng)用實驗室能力驗證結(jié)果評定蒸餾酒中甲醇濃度的不確定度[J].檢驗檢疫科學(xué)，2007，17(增刊)：8-10．

[11] 尹蘭風(fēng).復(fù)現(xiàn)性限與不確定度評定[J].中國計量，2002(3)：44-45.

[12] Eurachem Workshop.Recent Developments in Measure-ment Uncertainty [R].2011.

[13] The New Draft Eurachem/CITAC Guide: A Summary of Structure and Changes. Quantifying Uncertainty in Analytical Measurement(MU指南第3版) [R].2011.

[14] Van der Veen A M H.Revision of the Guide to the Expression of Uncertainty in Measurement [R].2011.

[15] Magnusson B,Naykki T,Hovind H. NORDTEST Report TR537：Handbook for Calculation of Measurement Uncertainty in Environmental Laboratories [R]. 2003.

[16] European Federation of National Associations of Meas-urement,Testing and Analytical Laboratories(EUROLAB).Measurement Uncertainty Revisited: Alternative Approaches to Uncertainty Evaluation[R].2007.

[17] Ramsey M H,Ellison S L R. Eurachem/EUROLAB/CITAC/Nordtest/AMC Guide: Measurement Uncertainty Arising from Sampling: A Guide to Methods and Approaches[S].Eurachem ISBN 978 0 948926 26 6(http:∥www.eurachem.org/guides/UfS_2007.pdf).2007.

[18] 孫海容，李玉武. 用經(jīng)驗?zāi)Ｐ团袛嗄芰︱炞C上報數(shù)據(jù)統(tǒng)計分析結(jié)果合理性[J].巖礦測試，2012，31(6)：992-996.

[19] ISO 13528—2005,Statistical Methods for Use in Profi-ciency Testing by Inter-laboratory Comparisons [S].

[20] Analytical Methods Committee. Robust Statistics: A Method of Coping with Outliers [R].AMC Technical Brief.2001.

[21] 佟艷春. 拉伸試驗?zāi)芰︱炞C結(jié)果的穩(wěn)健統(tǒng)計指定值及其不確定度[J].理化檢驗(物理分冊)，2010，46(2)，70-83.

[22] 邢小茹，馬小爽，田文，吳忠祥.實驗室間比對能力驗證中的兩種穩(wěn)健統(tǒng)計技術(shù)探討[J].中國環(huán)境監(jiān)測，2011,27(4)：4-8.

[23] 中國合格評定國家認(rèn)可委員會.CNACNAS T0476低合金鋼中化學(xué)成分分析能力驗證計劃結(jié)果報告(附錄D) [R].2010.

[24] 中國實驗室國家認(rèn)可委員會.環(huán)境檢測領(lǐng)域基于質(zhì)控數(shù)據(jù)評定測量不確定度指南[S].2013.