李書明

(廣州大學(xué),廣州510006)

1 引言

隨著能源形勢日益緊張,人民日益關(guān)注建筑和暖通空調(diào)節(jié)能,而其中的工作如政府決策的數(shù)據(jù)數(shù)據(jù)確定、建筑和設(shè)備效能測評標(biāo)識、相關(guān)技術(shù)研發(fā)領(lǐng)域、室內(nèi)環(huán)境測評、用戶冷熱計量、能源測評審計與計量等都離不開實驗分析來選擇測量方法以及分析測量值的可靠性。而不確定度分析作為實驗分析和實驗結(jié)果處理的重要手段,也受到重視。

國際上早在1980[1]就已開始研究用不確定度(uncertainty)來表示實驗結(jié)果,更自1993年[2,3]大力推廣采用測量不確定度的概念。我國的測量行業(yè)標(biāo)準(zhǔn)JJF1059-1999[4,5]即是對GUM1995[3]的等同采用,對科學(xué)研究、工程技術(shù)及商貿(mào)中大量存在的數(shù)據(jù)的處理和表示均具有適用性。測量不確定度的定義[3～5]是:與測量結(jié)果相關(guān)聯(lián)的一個參數(shù),用以表征合理地賦予被測量之值的分散性。不確定度指導(dǎo)致測量結(jié)果偏離的程度,用分散性尺度來評價測量結(jié)果質(zhì)量,綜合了全部誤差因素對實驗結(jié)果的影響。

而我國當(dāng)前雖然也逐步的重視不確定度分析,如建筑節(jié)能檢測標(biāo)準(zhǔn)JGJ/T132-2009和JGJ/T177-2009均有涉及,但卻一直用誤差理論而非不確定度分析來指導(dǎo)暖通空調(diào)實驗[6]。

而誤差為測量值偏離真值的程度,以標(biāo)準(zhǔn)差來評價測量結(jié)果的質(zhì)量,反應(yīng)的是測量誤差對實驗結(jié)果的影響。若用誤差理論來表示被測量的值,即使對完全已知或猜測的誤差因素進(jìn)行補償、修正后,所得結(jié)果依然只能是被測量的一個估計值。另外,誤差理論通常只考慮實驗數(shù)據(jù)處理階段的誤差或不確定度,而對實驗規(guī)劃及設(shè)計階段的誤差及不確定度分析沒有考慮或不夠重視。然而正是實驗規(guī)劃及設(shè)計階段的不確定分析才能從根本上控制實驗的方案精確度和減少系統(tǒng)不確定度,因此需要推廣不確定度分析。

不確定度分析是優(yōu)化實驗的強大工具,可應(yīng)用于實驗分析的所有階段 (見表1),尤其是實驗的規(guī)劃和設(shè)計階段。根據(jù)待測對象確定實驗方案和測量方法,選擇適當(dāng)?shù)臏y量儀器,搭建試驗臺,實施實驗,經(jīng)過 “實驗—修改—再實驗”,正確得出并報告實驗結(jié)果。

表1 實驗不確定度分析

2 不確定度分析方法

根據(jù)要求不同,不確定度按有標(biāo)準(zhǔn)不確定度、合成不確定度、擴展不確定度三種表述方式。根據(jù)估計方法不同,其分為按統(tǒng)計分布估計的A類不確定度和按非統(tǒng)計分布估計的B類不確定度兩類。

2.1 A類標(biāo)準(zhǔn)不確定度

1)在相同測量條件下,對同一被測量進(jìn)行連續(xù)多次測量的復(fù)制性條件下,采用貝塞爾公式計算得到:樣本均值即測量結(jié)果,樣本標(biāo)準(zhǔn)差即測量結(jié)果的標(biāo)準(zhǔn)不確定度。

2)如果被測量隨時間隨機變化,即為隨機過程,則應(yīng)采用較差法 (阿倫方差法)求標(biāo)準(zhǔn)不確定度。

例如測定f實驗,且f=y(τ)。對被測量進(jìn)行m+1次測量,每次測量取樣時間為τ,共有m個較差 [yi+1(τ)-yi(τ)];阿侖方差:則f的A 類標(biāo)準(zhǔn)不確定度:

3)當(dāng)被測量X的估計值是由實驗數(shù)據(jù)用最小二乘法擬合的一條曲線上得到,則標(biāo)準(zhǔn)不確定度可由擬合函數(shù)和不確定度傳遞率確定。

如對Y=mX+c來說:

式中,s(m),s(c),s(X)分別為m,c,X的估計標(biāo)準(zhǔn)差;r(m,c)為m和c的相關(guān)系數(shù)。

2.2 B類標(biāo)準(zhǔn)不確定度

標(biāo)準(zhǔn)不確定度u(f)可用X的有關(guān)信息或資料來評定。信息來源可以是:以前測量的數(shù)據(jù);經(jīng)驗和對有關(guān)儀器性能或材料的資料;廠家技術(shù)說明書;鑒定證書、校準(zhǔn)證書、測試報告或其他提供數(shù)據(jù)的文件等。

評定方法為:一般根據(jù)經(jīng)驗或有關(guān)信息和資料,分析判斷被測量的可能值的區(qū)間(-α,α),并假設(shè)測量誤差的概率分布,則不確定度為:u(x)=α/k,式中 α為區(qū)間的半寬度,k為包含因子,由要求的置信水平估計。

2.3 合成不確定度

當(dāng)測量結(jié)果為其他被測量求得時,按各被測量的方差和協(xié)方差來計算合成不確定度uc(y)。在實際工作中,當(dāng)被測量接近高斯正態(tài)分布時,可采用測量結(jié)果y及其合成不確定度uc(y)表示,即Y=y±uc(y)

此式是基于y=f(x1,x2,…,xn)的泰勒級數(shù)的一階近似,稱為 “不確定度傳播率”。偏導(dǎo)數(shù)Ci=?f/?xi是在Xi=xi時評定的,稱為靈敏系數(shù)。式中:x,x是X,X的估計值,且u(x,x)=u(x,xi)是xi,xj的協(xié)方差的估計值。u(xi)可以由A類或B類得到;

協(xié)方差與相關(guān)系數(shù)r(xi,xj)的關(guān)系是:r的絕對值大于0.7時,一般稱為強相關(guān),否則為弱相關(guān)。當(dāng)在估計輸入量之值xi時,如果使用了同一個測量標(biāo)準(zhǔn)、測量儀器、參考數(shù)據(jù)或甚至具有相當(dāng)大不確定度的測量方法即會產(chǎn)生明顯相關(guān),就必須考慮其相關(guān)性。

2.4 擴展不確定度

擴展不確定度用U表示,指確定結(jié)果區(qū)間的量,被測量可望大部含于此區(qū)間。U由合成不確定度uc(y)乘包含因子k得到。

由于用合成不確定度表示測量結(jié)果僅有P=68%的置信概率,所以一般采用擴展不確定度來表示測量結(jié)果,即Y=y±U。

求擴展不確定度方法:

1)取擴展不確定度U=kuc,一般取置信概率為95%的包含因子k95=2,有時取置信概率為99%的k99=2.6(大樣本近似計算值,JJF1059取3);在多數(shù)的檢定/校準(zhǔn)結(jié)果中,其合成標(biāo)準(zhǔn)不確定度均可用第一種簡單的方法進(jìn)行計算。

2)如果uc(y)的自由度較小,并要求區(qū)間具有規(guī)定的置信水平時,取U=kPuc。式中k與P由各影響分量自由度Vi及合成標(biāo)準(zhǔn)不確定度的有效自由度Veff來確定。A類評定時,通常對每個Xi測量ni次,自由度為vi=ni-1;如果有n個數(shù)據(jù)點用最小二乘法擬合成曲線來估計m個參數(shù),則每個參數(shù)的標(biāo)準(zhǔn)不確定度自由度為n-m。B類評定時,由相對標(biāo)準(zhǔn)不確定度 (由經(jīng)驗定)帶入式確定。而當(dāng)Veff足夠大時,兩法相同。

3 實驗用不確定度分析

3.1 總體不確定度分析

通過不確定度分析,實驗人員可以從規(guī)劃及設(shè)計階段綜合考慮系統(tǒng)精度、經(jīng)濟(jì)性或簡便性等實驗的各個方面,即在達(dá)到實驗結(jié)果要求的精度下減少資源消耗;通常有不同的實驗方法或多種實驗方案,對于每一種方法,有不同的不確定度傳遞特性,需確定哪一種方法的測量結(jié)果最接近真值,并分析所用方法中哪個量的測量是最重要的。如對同一物理量的測量,應(yīng)根據(jù)研究對象,列出各種可能的實驗方法,分析各種方法的適用條件,比較各種方法的局限性與可能達(dá)到的實驗精確度 (或不確定度)等因素,并考慮方案實施的可能性,最后選出最佳的實驗方法。

若不確定度傳遞方程如:

則實驗結(jié)果的不確定度表述為:

式中UMFi為不確定度放大系數(shù) (正態(tài)敏感系數(shù)),表征變量對結(jié)果不確定度的影響;UPCi為不確定度貢獻(xiàn)百分比,包含有UMFi和變量i的不確定度量級的效果,因此在實驗設(shè)備和方法選擇完成后,便可估計出測量不確定度。

3.2 詳細(xì)不確定度分析

在進(jìn)行詳細(xì)不確定度分析時,根據(jù)誤差源分別評定每個測量量的系統(tǒng)和隨機不確定度,然后再合成為整體不確定度。若測量方程為式(5),則系統(tǒng)不確定度和隨機不確定可以按下式(TSM方法)計算:

式中bX,sX為X變量的系統(tǒng)不確定度和隨機不確定度,分別指每兩個有相同系統(tǒng)誤差源變量的協(xié)方差和每兩個有相同隨機誤差源變量的協(xié)方差。而隨機不確定度相關(guān)項的協(xié)方差一般為0。

在穩(wěn)態(tài)測試中,如果式 (5)中為單變量測量且有足夠多的樣本數(shù),可采用貝塞爾公式計算隨機不確定度:

顯然,直接法考慮了各被測量之間的相關(guān)性,有著可能更精確的結(jié)果。但直接法需要通過多次測量獲得足夠的數(shù)據(jù)信息,測量的時間長度也是很重要的,所以在單次測量中一般不應(yīng)使用直接法。

在實驗的調(diào)試階段,可用不確定度分析來發(fā)現(xiàn)實驗設(shè)備、操作運行、測試方法等方面的問題并加以改進(jìn)。用1階復(fù)制水平進(jìn)行運行逐步檢查,使用N階擴展不確定度評估對接受的結(jié)果和理論值的結(jié)果進(jìn)行比較,如果在結(jié)果中隨機不確定度比預(yù)期值大得多,則表明某些隨機誤差源被低估了,需繼續(xù)調(diào)試。

在試驗的實施階段,可用不確定度分析來進(jìn)行測試點的校正。1)校正測點范圍:因為高階回歸方程需要較多的測點,盡量選用適當(dāng)?shù)淖鴺?biāo)系 (如對數(shù)坐標(biāo)系)和盡量低階的關(guān)系式,可減少必要的測點數(shù)量。2)設(shè)定測試序列:一般情況下,隨機地分布測試點得到的測試結(jié)果更接近實際情況,因此我們設(shè)定測試點是隨機的而不是連續(xù)的。

3.3 實驗結(jié)果處理

通過回歸分析獲得數(shù)學(xué)表達(dá)式的方法大量應(yīng)用于實驗數(shù)據(jù)處理上。實驗數(shù)據(jù)回歸方程Y=mX+c中Y、m、c的回歸標(biāo)準(zhǔn)差:

當(dāng)不考慮X的隨機不確定度時,sY可近似為回歸隨機不確定度,可用來衡量所有隨機因素對Y影響的大小,sY愈小,回歸精度愈高。但sY僅僅表達(dá)了擬合曲線的隨機不確定度,并不包含參考標(biāo)準(zhǔn)、數(shù)據(jù)獲取、試驗臺搭建、概念模型引起的不確定度中的任何一項。

4 應(yīng)用舉例

4.1 整體不確定度分析例子

為了測定某圓管外壁對流換熱系數(shù),分析兩種實驗方案:

1)穩(wěn)態(tài)方法:

根據(jù)能量守恒W=h(πDL)(T-Tα)得到:

式中電功率W、圓管外徑D、長度L和圓柱外壁溫T、來流溫度Tα均通過測量獲得;

2)非穩(wěn)態(tài)方法:

式中圓柱體的質(zhì)量M、直徑D和長度L,測試時間τ,圓柱比熱c通過測量或查資料獲得。

(1)試驗中,忽略輻射和熱傳導(dǎo)的影響,h和c為常數(shù),忽略圓柱內(nèi)部熱阻,即溫度均勻分布。若溫度是單獨測量的,應(yīng)用不確定度傳遞方程(15)有:

表2 兩種方法各式中測量值的不確定度表

(2)用A類或B類方法評定各測量量的不確定度 (表2),并將數(shù)值帶入 (15～16)式得到 (17～18)式:

圖1 穩(wěn)態(tài)測試圖 (左);非穩(wěn)態(tài)測試圖 (右上);對比圖 (右下)

(3)結(jié)果分析

穩(wěn)態(tài)測試總體不確定度分析:由相對不確定度傳遞公式看出,對結(jié)果不確定度影響,△T最大,而D相對于△T的來說可以忽略,W的影響在W值大時較小,即對總體不確定度影響較小,相反影響較大;從測試線圖可以看出,在h較大時,不同溫差條件下的測試均可滿足總體不確定度的要求,這是因為h較大時,W=0.2426亦較大,W相對不確定度較小。同樣原因,在h較小時,W為影響總體不確定度的關(guān)鍵因素,只有在低于0.5W的情況下才能得到h總體不確定度5%的結(jié)果。

非穩(wěn)態(tài)測試總體不確定度分析:由相對不確定度傳遞公式看出,對結(jié)果不確定度影響,c較大,可以通過改善c值來改進(jìn)實驗;由線圖可以看出,隨著h值的增大,τ逐漸成為總體不確定度的決定性因素。并在大h值的領(lǐng)域,要獲得5%的h總體不確定度,需要τ為7ms,8ms或更少時間。

對比分析:從不確定度傳遞方程線圖可以看出,穩(wěn)態(tài)測試h的總體不確定度隨著h的增大而減小,非穩(wěn)態(tài)測試h的總體不確定度隨著h的增大而增大;在h＜150時,宜采用非穩(wěn)態(tài)測試;在h＞200時,宜采用穩(wěn)態(tài)測試;而在150＜h＜200時,兩者基本相當(dāng) [h單位為W/(m2℃)]。

4.2 詳細(xì)不確定度分析例子

對某煤樣品的熱值測量數(shù)據(jù)進(jìn)行處理。其測量數(shù)學(xué)方程為:

H=式中H為熱值 (W/g);(T2-T1)為燃后溫升 (℃);C熱量計熱當(dāng)量(W/℃);e為燃燒放熱量 (W);M為樣本質(zhì)量。此處是一個樣本對樣本實驗。

4.2.1 0階復(fù)制水平分析

表3 0階水平變量不確定度估計表

系統(tǒng)不確定度:

隨機不確定度 (由上公式計算H=18981W/g):

即sH=199 W/g;

合成不確定度:

4.2.2 1階復(fù)制水平分析

對同一樣本進(jìn)行了26次測量,并求出樣本均值=18578W/g和用貝塞爾公式計算出標(biāo)準(zhǔn)差SH=364.6W/g,但可以看出該隨機不確定度較0階水平算出的199W/g大,這是應(yīng)為1階計算結(jié)果考慮了樣本對樣本實驗中測量系統(tǒng)和材料組分的效果。

因樣本變化的隨機不確定度為:

上式結(jié)果說明,因為煤組分的變化而引起的隨機不確定度,實際上比測量系統(tǒng)的隨機不確定度還要大。在樣本對樣本試驗中,可以利用0階復(fù)制水平隨機不確定度和多次測量1階隨機不確定度后,來確定樣本變化對結(jié)果不確定度的影響。

4.2.3 N階復(fù)制水平分析

在N階復(fù)制水平下,使用0階水平計算的系統(tǒng)不確定度和1階直接測量的隨機不確定度來計算合成不確定度和擴展不確定度 (95%置信概率)。

得到置信概率為95%的測量結(jié)果:18758±180W/g;或18800W/g,1%相對不確定度。

5 結(jié)論

在暖通空調(diào)實驗測量領(lǐng)域,包括建筑能耗檢測、建筑設(shè)備測評、暖通實驗室測量等眾多領(lǐng)域的實驗規(guī)劃、實驗設(shè)計、實驗報告等領(lǐng)域應(yīng)推廣應(yīng)用不確定度分析,以在滿足試驗要求的情況下節(jié)省資源和加大試驗的可操作性。

不確定度分析理論是誤差理論的進(jìn)一步發(fā)展,主要按不確定度來源分為系統(tǒng)不確定度和隨機不確定度,系統(tǒng)不確定度可以根據(jù)經(jīng)驗、資料、儀器說明書評定,而隨機不確定度可以根據(jù)直接測量樣本標(biāo)準(zhǔn)差來獲得。

在實驗規(guī)劃階段,可以通過初算不確定度,校核初選實驗方案是否滿足實驗要求的精度,以選定實驗方案;另外可以根據(jù)TSM法不確定度傳遞方程進(jìn)行分析,計算敏感系數(shù)等可觀值來找到對實驗結(jié)果起關(guān)鍵影響的因素,改進(jìn)實驗,選擇實驗設(shè)備,優(yōu)化實驗方案。在試驗實施階段,可以通過監(jiān)測實驗不確定度,監(jiān)測實驗設(shè)備的運行情況,以達(dá)到優(yōu)良的運行效果。在試驗報告階段,以方便報告不確定度的形式組織線圖和回歸方程形式,并按標(biāo)準(zhǔn)要求給出測試線圖和回顧方程各部分的不確定度。

[1] BIP M.Recommendation INC-1.CIPM 81-11,1981

[2] The 1993 ISO Guide to the Expression of Uncertainty?in Measurement(GUM)

[3] The 1995 ISO Guide to the Expression of Uncertainty?in Measurement(GUM)

[4] 中國國家計量技術(shù)規(guī)范——通用計量術(shù)語和定義JJF1001-1998[S]

[5] 中國國家計量技術(shù)規(guī)范——測量不確定度評價與表示JJF1059-1999[S]

[6] 李峰,姬長發(fā).建筑環(huán)境與設(shè)備工程實驗及測試技術(shù)[M].北京:機械工業(yè)出版社,2008,9

[7] The 2005 standard ASME PTC19.1 Test Uncertainty

[8] The 2008 JCGM Supplement 1 to the GUM:Evaluation of Measurement Data—Propagation of Distributions Using a Monte Carlo Method

[9] Coleman H W,Steele W G,Experimentation,Validation,and Uncertainty Analysisfor Engineers.New York:John Wiley&Sons,2009

[10] 劉智敏.不確定度原理[M].北京:中國計量出版社,1993,6

[11] 宋明順.測量不確定度評定與數(shù)據(jù)處理[M].北京:中國計量出版社,2000,1

[12] 葉培德.測量不確定度[M].北京:國防工業(yè)出版社,1996,5

[13] 李慎安.測量誤差與測量不確定度之間的聯(lián)系與異同[J].城市技術(shù)監(jiān)督,1994,2:12-13

[14] 李慶忠.不確定度評估方法的六個關(guān)鍵點[J].工程與試驗,2008,48(B12):16-18

[15] 劉智敏.分布傳播及其在不確定度中的應(yīng)用[J].中國計量學(xué)院學(xué)報,2008,19(3):1-9