随着天然氣(qì)計量技術的發展(zhan)和對天然氣貿易(yì)交接計量要求的(de)提高,流量計的正(zheng)确計量顯得至關(guān)重要。在此背景下(xia),某省組織開展了(le)全省氣體渦輪流(liú)量計
計量比對工(gong)作,以确保全省量(liàng)值統一、正确、可靠(kao),促進計量技術機(jī)構能力提升。
比對(duì)選用了一台DN150mm氣體(ti)渦輪流量計作爲(wei)傳遞标準,流量範(fan)圍爲(80~1600)m³/h,精度等級爲(wèi)1.0級。其中0.4qmax=640m³/h、.0.7qmax=1120m³/h爲比對流(liu)量點,每點檢測6次(ci),測得其儀表系數(shu)K,取6次平均值作爲(wei)比對值。
本次比對(dui),各參比實驗室完(wan)成率爲100%,比對完成(cheng)後,主導實驗室對(duì)各參比實驗室的(de)數據進行了彙總(zong)統計,各實驗室測(cè)量結果如圖1所示(shi)。由圖1可以看出,實(shi)驗室3的測量結果(guo)随着流量增大而(er)增大,與其他參比(bi)實驗室的測量結(jie)果産生了較大偏(pian)移,需進行異常值(zhí)判别。
一、異常值判(pàn)别常用的統計方(fang)法
判别測量結果(guǒ)異常值常用的統(tǒng)計方法有拉依達(dá).準則(3σ準則)、格拉布(bù)斯準則、狄克遜準(zhǔn)則等。-般來說,在重(zhong)複觀測次數n>50的前(qian)提下,可以考慮使(shi)用3σ準則,但在GB/T4883-2008《數據(ju)的統計處理和解(jiě)釋正态樣本離群(qún)值的判斷和處理(li)》中已不采用此方(fang)法;在3<n<50的情況下,格(ge)拉布斯準則效果(guo)較好,适用于單個(ge)異常值;有多于一(yi)個異常值時,選用(yòng)狄克遜準則較好(hǎo)。
二、氣體渦輪流量(liàng)計計量比對數據(jù)中異常值判别
在(zài)此次氣體渦輪流(liú)量計計量比對工(gong)作過程中,異常值(zhí)需要慎重處理。爲(wèi)避免異常值的遺(yí)漏和測量真值的(de)誤判,可以綜合運(yun)用多種判定方法(fǎ)處理異常值。考慮(lǜ)比對工作的重要(yào)性,此次異常值判(pàn)别選用格拉布斯(sī)準則和狄克遜準(zhǔn)則同時進行,選.用(yong)格拉布斯準則判(pan)别異常值,并用狄(dí)克遜準則判别以(yi)作驗證。
1.格拉布斯(si)準則
設在一組重(zhòng)複觀測值xi中,其殘(cán)差vi的絕對值lvil最大(da)者爲可疑值xd,在給(gěi)定的包含概率爲(wei)p=0.95或p=0.99,也就是顯著性(xìng)水平爲α=1-p=0.01或0.05時,如果(guǒ)滿足式(1),可以判定(ding)xa爲異常值。
式中:G(α,n)一(yī)與顯著性水平α和(hé)重複觀測次數n有(yǒu)關的格拉布斯臨(lín)界值,查格拉布斯(si)準則臨界值G(α,n)表可(kě)得;s一測量結果的(de)實驗标準偏差。
本(běn)次計算中,測量結(jié)果平均值`x和實驗(yan)标準偏差s以算術(shu)平均值計算,步驟(zhou)如下:
(1)0.4qmx流量點
計算(suan)過程如表1所示。參(can)比實驗室`x值的算(suan)術平均值爲x=5261.28。
實驗(yàn)标準偏差爲
計算(suàn)各個測量結果的(de)殘差xj-`x,其中絕對值(zhi)最大的殘差爲32.28,相(xiàng)應的測量結果可(kě)疑值爲x3=5293.6。按置信水(shui)平P=95%=0.95,即α=1-0.95=0.05,n=13,查格拉布斯(sī)準則的臨界值表(biǎo)得到G(0.05,13)=2.331,則
可以判定(ding)實驗室3在0.4qmax.流量點(dian)的測量結果不屬(shǔ)于異常值,不需要(yao)剔除。
(2)0.7qms流量點.
計算(suàn)過程如表2所示。參(cān)比實驗室K值的算(suàn)術平均值爲`x=5266.87。
實驗(yan)标準偏差爲
計算(suan)各個測量結果的(de)殘差,其中絕對值(zhí)最大的殘差爲48.24,相(xiang)應的測量結果可(ke)疑值爲X3=5315.11。按置信水(shuǐ)平P=95%=0.95,即α=1-0.95=0.05,n=13,查格拉布斯(si)準則的臨界值表(biao)得到
可以判定實(shí)驗室3在0.7qmax流量點的(de)測量結果屬于異(yì)常值,需要剔除。剔(tī)除後,12組數據重複(fú)上述過程,計算結(jie)果如表3所示,置信(xìn)水平P=95%,即α=1-0.95=0.05,n=12,查格拉布(bu)斯準則的臨界值(zhi)表得到G(0.05,12)=2.285,可以判斷(duan)無異常值。
2.狄克遜(xun)準則驗證
設,所得(dé)的重複觀測值按(an)由小到大的順序(xù)排列爲X1,X2.....,xn。按以下4種(zhǒng)情況計算統計量(liàng)γij,或γij'
(1)在n=3~7的情況下:
(2)在(zai)n=8~10的情況下:
判斷異(yì)常值的狄克遜準(zhun)則爲
若γij>γij',ry>D(α,),則xn爲異常(cháng)值。
若γij>γij',ry>D(α,),則x1爲異常值(zhi)。
式中:D(α,n)一狄克遜臨(lin)界值,查狄克遜檢(jiǎn)驗的臨界值D(α,n)表可(kě)得。
本次計算的具(jù)體步驟如下:
測量(liàng)次數n=13,将各實驗室(shi)測量結果按從小(xiao)到大的順序排列(liè)爲x1,x2,....X13。按置信水平P=95%,即(ji)α=1-0.95=0.05,n=13,查狄克遜準則的(de)臨界值表得到G(0.05,13)=0.557。
在(zai)n=11~13的情況下:
若γ21>γ21',γ21>D(α,n),則xn爲(wei)異常值。
若γ21>γ21',γ21>D(α,n),則x1爲異(yì)常值。
否則沒有異(yi)常值。
(1)0.4qmax流量點
0.4qmax.流量(liang)點的測量結果如(ru)表4所示。
由于是屬(shu)于n=11~13的情況,所以統(tǒng)計量計算如下
γ21>γ21',則(ze)
γ21=0.3655<D(0.0513)=0.557
可以判定實驗室(shì)3在0.4qmax流量點的測量(liang)結果X13=5293.56不屬于異常(cháng)值,不需要剔除。
(2)0.7qmax流(liú)量點
0.7qmax流量點的測(ce)量結果如表5所示(shì)
由于是屬于n=11~13的情(qing)況,所以統計量計(ji)算如下
γ2>γ21',則
γ21=0.5549<D(0.05,13)=0.557
可以判(pàn)定實驗室3在0.7.流量(liàng)點的測量結果X13=5315.11不(bú)屬于異常值,不需(xu)要剔除。
三、異常值(zhí)判别結果分析
綜(zong)上,考慮比對工作(zuò)的嚴謹性,在此次(cì)異常值判别時,選(xuan)用格拉布斯準則(zé)和狄克遜準則綜(zōng)合進行。通過異常(chang)值判别結果分析(xī)可知:在0.4qmax流量點,實(shi)驗室3的測量結果(guo)經過兩種判别準(zhǔn)則判斷均不屬于(yu)異常值,可正常參(can)與計算;在0.7max.流量點(dian),實驗室3的測量結(jie)果用格拉布斯準(zhun)則判别屬于異常(cháng)值,用狄克遜準則(ze)判别不屬于異常(cháng)值,兩種判别方法(fǎ)結論出現矛盾。比(bi)對工作對統計數(shù)值要求較高,出現(xiàn)既可能是異常值(zhi)又可能不是異常(chang)值的情況時,考慮(lü)将此次測量結果(guǒ)作爲異常值考慮(lǜ)。最終将0.7qmax流量點實(shi)驗室3的測量結果(guǒ)作爲異常值,予以(yǐ)剔除。
四、結束語.
對(duì)于異常值判别方(fāng)法的選擇,很多學(xue)者做了大量研究(jiū),但大多集中在研(yán)究如何利用工具(jù)輔助快速正确地(di)識别異常值,以節(jie)省時間提高工作(zuò)效率。目前對于判(pan)别方法的選擇還(hái)沒有定論,達成比(bi)較-.緻的意見是在(zai)涉及數據要求較(jiào)高的場合,可以綜(zōng)合運用多種判别(bie)方法進行異常值(zhi)判别。在實際工作(zuo)中,計量檢定人員(yuán)除了提升工作技(ji)能,盡量避免人爲(wèi)因素摻人異常值(zhi)以外,更要熟悉并(bing)掌握多種異常值(zhi)判别方法,以保證(zhèng)日常檢定工作中(zhōng)數據的有效性.和(he)科學性。
本文來源(yuán)于網絡,如有侵權(quan)聯系即删除!