隨著天然氣計量技術(shù)的發(fā)展和對天然氣貿(mào)易交接計量要求的提高，流量計的正確計量顯得至關(guān)重要。在此背景下，某省組織開展了全省氣體渦輪流量計計量比對工作，以確保全省量值統(tǒng)一、正確、可靠，促進計量技術(shù)機構(gòu)能力提升。
比對選用了一臺DN150mm氣體渦輪流量計作為傳遞標(biāo)準(zhǔn)，流量范圍為(80~1600)m³/h,精度等級為1.0級。其中0.4qmax=640m³/h、.0.7qmax=1120m³/h為比對流量點，每點檢測6次，測得其儀表系數(shù)K,取6次平均值作為比對值。
本次比對，各參比實驗室完成率為100%，比對完成后，主導(dǎo)實驗室對各參比實驗室的數(shù)據(jù)進行了匯總統(tǒng)計，各實驗室測量結(jié)果如圖1所示。由圖1可以看出，實驗室3的測量結(jié)果隨著流量增大而增大，與其他參比實驗室的測量結(jié)果產(chǎn)生了較大偏移，需進行異常值判別。
一、異常值判別常用的統(tǒng)計方法
判別測量結(jié)果異常值常用的統(tǒng)計方法有拉依達.準(zhǔn)則(3σ準(zhǔn)則)、格拉布斯準(zhǔn)則、狄克遜準(zhǔn)則等。-般來說，在重復(fù)觀測次數(shù)n>50的前提下，可以考慮使用3σ準(zhǔn)則，但在GB/T4883-2008《數(shù)據(jù)的統(tǒng)計處理和解釋正態(tài)樣本離群值的判斷和處理》中已不采用此方法;在3<n<50的情況下，格拉布斯準(zhǔn)則效果較好，適用于單個異常值;有多于一個異常值時，選用狄克遜準(zhǔn)則較好。
二、氣體渦輪流量計計量比對數(shù)據(jù)中異常值判別
在此次氣體渦輪流量計計量比對工作過程中，異常值需要慎重處理。為避免異常值的遺漏和測量真值的誤判，可以綜合運用多種判定方法處理異常值。考慮比對工作的重要性，此次異常值判別選用格拉布斯準(zhǔn)則和狄克遜準(zhǔn)則同時進行，選.用格拉布斯準(zhǔn)則判別異常值，并用狄克遜準(zhǔn)則判別以作驗證。
1.格拉布斯準(zhǔn)則
設(shè)在一組重復(fù)觀測值xi中，其殘差vi的絕對值lvil最大者為可疑值xd，在給定的包含概率為p=0.95或p=0.99,也就是顯著性水平為α=1-p=0.01或0.05時，如果滿足式(1)，可以判定xa為異常值。
式中:G(α，n)一與顯著性水平α和重復(fù)觀測次數(shù)n有關(guān)的格拉布斯臨界值，查格拉布斯準(zhǔn)則臨界值G(α,n)表可得;s一測量結(jié)果的實驗標(biāo)準(zhǔn)偏差。
本次計算中，測量結(jié)果平均值`x和實驗標(biāo)準(zhǔn)偏差s以算術(shù)平均值計算，步驟如下:
(1)0.4qmx流量點
計算過程如表1所示。參比實驗室`x值的算術(shù)平均值為x=5261.28。
實驗標(biāo)準(zhǔn)偏差為
計算各個測量結(jié)果的殘差xj-`x，其中絕對值最大的殘差為32.28，相應(yīng)的測量結(jié)果可疑值為x3=5293.6。按置信水平P=95%=0.95,即α=1-0.95=0.05，n=13,查格拉布斯準(zhǔn)則的臨界值表得到G(0.05，13)=2.331，則
可以判定實驗室3在0.4qmax.流量點的測量結(jié)果不屬于異常值，不需要剔除。
(2)0.7qms流量點.
計算過程如表2所示。參比實驗室K值的算術(shù)平均值為`x=5266.87。
實驗標(biāo)準(zhǔn)偏差為
計算各個測量結(jié)果的殘差，其中絕對值最大的殘差為48.24，相應(yīng)的測量結(jié)果可疑值為X3=5315.11。按置信水平P=95%=0.95,即α=1-0.95=0.05，n=13,查格拉布斯準(zhǔn)則的臨界值表得到
可以判定實驗室3在0.7qmax流量點的測量結(jié)果屬于異常值，需要剔除。剔除后，12組數(shù)據(jù)重復(fù)上述過程，計算結(jié)果如表3所示，置信水平P=95%，即α=1-0.95=0.05，n=12，查格拉布斯準(zhǔn)則的臨界值表得到G(0.05，12)=2.285，可以判斷無異常值。
2.狄克遜準(zhǔn)則驗證
設(shè)，所得的重復(fù)觀測值按由小到大的順序排列為X1，X2.....，xn。按以下4種情況計算統(tǒng)計量γij,或γij'
(1)在n=3~7的情況下:
(2)在n=8~10的情況下:
判斷異常值的狄克遜準(zhǔn)則為
若γij>γij'，ry>D(α，),則xn為異常值。
若γij>γij'，ry>D(α，),則x1為異常值。
式中:D(α，n)一狄克遜臨界值，查狄克遜檢驗的臨界值D(α,n)表可得。
本次計算的具體步驟如下:
測量次數(shù)n=13，將各實驗室測量結(jié)果按從小到大的順序排列為x1，x2，....X13。按置信水平P=95%，即α=1-0.95=0.05,n=13，查狄克遜準(zhǔn)則的臨界值表得到G(0.05,13)=0.557。
在n=11~13的情況下:
若γ21>γ21',γ21>D(α，n)，則xn為異常值。
若γ21>γ21'，γ21>D(α，n)，則x1為異常值。
否則沒有異常值。
(1)0.4qmax流量點
0.4qmax.流量點的測量結(jié)果如表4所示。
由于是屬于n=11~13的情況，所以統(tǒng)計量計算如下
γ21>γ21',則
γ21=0.3655<D(0.0513)=0.557
可以判定實驗室3在0.4qmax流量點的測量結(jié)果X13=5293.56不屬于異常值，不需要剔除。
(2)0.7qmax流量點
0.7qmax流量點的測量結(jié)果如表5所示
由于是屬于n=11~13的情況，所以統(tǒng)計量計算如下
γ2>γ21'，則
γ21=0.5549<D(0.05，13)=0.557
可以判定實驗室3在0.7.流量點的測量結(jié)果X13=5315.11不屬于異常值，不需要剔除。
三、異常值判別結(jié)果分析
綜上，考慮比對工作的嚴(yán)謹(jǐn)性，在此次異常值判別時，選用格拉布斯準(zhǔn)則和狄克遜準(zhǔn)則綜合進行。通過異常值判別結(jié)果分析可知:在0.4qmax流量點，實驗室3的測量結(jié)果經(jīng)過兩種判別準(zhǔn)則判斷均不屬于異常值，可正常參與計算;在0.7max.流量點，實驗室3的測量結(jié)果用格拉布斯準(zhǔn)則判別屬于異常值，用狄克遜準(zhǔn)則判別不屬于異常值，兩種判別方法結(jié)論出現(xiàn)矛盾。比對工作對統(tǒng)計數(shù)值要求較高，出現(xiàn)既可能是異常值又可能不是異常值的情況時，考慮將此次測量結(jié)果作為異常值考慮。最終將0.7qmax流量點實驗室3的測量結(jié)果作為異常值,予以剔除。
四、結(jié)束語.
對于異常值判別方法的選擇，很多學(xué)者做了大量研究，但大多集中在研究如何利用工具輔助快速正確地識別異常值，以節(jié)省時間提高工作效率。目前對于判別方法的選擇還沒有定論，達成比較-.致的意見是在涉及數(shù)據(jù)要求較高的場合，可以綜合運用多種判別方法進行異常值判別。在實際工作中，計量檢定人員除了提升工作技能，盡量避免人為因素?fù)饺水惓Ｖ狄酝猓煜げ⒄莆斩喾N異常值判別方法，以保證日常檢定工作中數(shù)據(jù)的有效性.和科學(xué)性。
本文來源于網(wǎng)絡(luò),如有侵權(quán)聯(lián)系即刪除！
