家用燃?xì)庠顭峁ば阅軠y量及不確定度評定*

胡偉健，李志宏

0 引言

家用燃?xì)庠钭鳛槲覈彝ブ衅占奥首罡叩囊环N廚房加熱器具，其熱工性能（熱負(fù)荷和熱效率）越來越受到消費者的重視。2008年5月1日開始實施的GB16410-2007要求嵌入式灶熱效率大于等于50%，臺式灶大于等于55%，同時規(guī)定，熱負(fù)荷偏差應(yīng)該在±10%以內(nèi)。而最新版本的GB16410—201X《家用燃?xì)庠罹摺罚ㄕ髑笠庖姼澹﹦t要求嵌入式灶熱效率大于等于55%，臺式灶大于等于58%，同時規(guī)定，熱負(fù)荷偏差應(yīng)該在±10%以內(nèi)。

2012年12月3日，國家質(zhì)量監(jiān)督檢驗檢疫總局發(fā)布《測量不確定度評定與表示》（GUM法）（JJF1059.1-2012）代替JJF1059-1999。2012年12月21日發(fā)布《用蒙特卡洛法評定測量不確定度》（JJF1059.2-2012），將蒙特卡洛法作為對GUM法的重要補充。根據(jù)我國《測量不確定度政策》的規(guī)定，檢測實驗室應(yīng)有能力對每一項有數(shù)值要求的測量結(jié)果進(jìn)行測量不確定度評定。當(dāng)不確定度與檢測結(jié)果的有效性或應(yīng)用有關(guān)、當(dāng)測試方法中有規(guī)定時和認(rèn)可委員會有要求時，檢測報告必須提供測量結(jié)果的不確定度[1]。

熱負(fù)荷和熱效率作為家用燃?xì)庠顭嵝阅艿闹匾笜?biāo)，有必要對其測量結(jié)果進(jìn)行不確定度評定。由于其測量模型比較復(fù)雜，且為非線性，采用GUM法評定不確定度，所求偏導(dǎo)比較困難，于是基于蒙特卡洛法的不確定度評定是有效的代替方法。

1 用蒙特卡洛法評定測量不確定度

1.1 蒙特卡洛法簡介

蒙特卡洛法主要是用相關(guān)函數(shù)產(chǎn)生服從相應(yīng)概率分布的隨機數(shù)對各個輸入量進(jìn)行多次隨機抽樣，然后再進(jìn)行概率分布傳播，最后得到輸出量的模擬結(jié)果，所以只要試驗次數(shù)足夠大，最后的模擬結(jié)果便近似于事件發(fā)生的概率，便可估計輸出量的有關(guān)特征量。近年來，該方法被廣泛的應(yīng)用于測量不確定度的評定，其核心是求取模擬樣本的標(biāo)準(zhǔn)差。其依據(jù)是分布傳播，即能夠通過輸入量的概率分布確定輸出量的概率分布。圖1是概率分布傳遞的示意圖。求解步驟按照用蒙特卡洛法評定測量不確定度[2-4]。

圖1 3個獨立的輸入量的概率分布傳遞示意圖Fig.1 Illustration of propagation of probability distributions for 3 independent input quantities

1.2 家用燃?xì)庠顭嶝?fù)荷測量不確定度評定示例

（1）測量原理和測量模型

熱負(fù)荷是指低熱值與燃?xì)饬髁康某朔e。本次試驗選用某燃?xì)庠钅苄Ъ熬C合性能測試臺對三種不同品牌的三個家用燃?xì)庠睿?2T）的左爐（分別標(biāo)記為1、2、3號灶）進(jìn)行熱負(fù)荷對比測試。如圖2所示，將燃?xì)庠钆c測試平臺連接。燃?xì)馔ㄟ^閥門流經(jīng)濕式氣體流量計，最后到燃?xì)庠睢Ｆ渲袧袷綒怏w流量計的燃?xì)饬髁俊⑷細(xì)鉁囟取⑷細(xì)忪o壓力以及大氣壓力等參數(shù)是通過測試平臺的傳感器采集，并進(jìn)行自動計算，最后將結(jié)果呈現(xiàn)在測試平臺的屏幕上的。測試按照GB30720—2014《家用燃?xì)庠罹吣苄薅ㄖ导澳苄У燃墶穂5]中的步驟進(jìn)行。

圖2 熱負(fù)荷測量原理示意圖Fig.2 Heat load measuring principlediagram

實測熱負(fù)荷用公式（1）計算：

φ實—實測熱負(fù)荷，單位為千瓦(kW)；

Q1—15℃、101.3 kPa狀態(tài)下試驗燃?xì)獾牡蜔嶂担瑔挝粸檎捉苟苛⒎矫?MJ/m3)；

V—測試耗氣量，單位為L；

tg—燃?xì)饬髁坑媰?nèi)的燃?xì)鉁囟龋瑔挝粸閿z氏度(℃)；

Pamb—試驗時的大氣壓力，單位為千帕(kPa)；

Pm——實測燃?xì)饬髁坑媰?nèi)的燃?xì)庀鄬o壓力，單位為千帕(kPa)；

S——溫度為tg時的飽和水蒸氣壓力，單位為千帕(kPa)(當(dāng)使用干式流量計測量時，S值應(yīng)乘以試驗燃?xì)獾南鄬穸冗M(jìn)行修正)；

t——試驗時間，單位秒（s）。

（2）測量數(shù)據(jù)

對三個家用燃?xì)庠畹母鱾€輸入量分別進(jìn)行10次獨立測量，計算得到輸入量的試驗標(biāo)準(zhǔn)差如表1。

表1 1、2、3號灶熱負(fù)荷各輸入量結(jié)果數(shù)據(jù)表Table1 Datasheet of heat input for 1，2 and 3 stoves

其中耗氣量V為10，是一個固定值，這樣只要記錄對應(yīng)時間就可以求得燃?xì)饬髁俊Ｌ烊粴獾牡蜔嶂礠1根據(jù)GBT11062-2014中低熱值體積發(fā)熱量的計算方法計算，本次試驗所用天然氣成分為甲烷含量≥99%，乙烷含量≤0.6%，空氣含量≤0.4%，具體計算過程如下：

已知組成混合物，在燃燒溫度t1，計量溫度t2和壓力p2時的理想氣體體積發(fā)熱量用公式（2）

式 Hˉ°(t1)為混合物的理想摩爾發(fā)熱量可查GB/T 11062-2014表3獲得； R為摩爾氣體常數(shù)（8.31451 J.mol-1.K-1）；T2為絕對溫度（T2=t1+273.15）。

當(dāng)乙烷含量為0.6%，空氣含量為0，甲烷含量為99.4%時低熱值取上限

當(dāng)乙烷含量為0，空氣含量為0.4%，甲烷含量為99.6%時低熱值取下限，

因此天然氣低熱值的取值范圍為[33.80，34.10]，根據(jù)最大熵原理[6-7]，天然氣低熱值為均勻分布，因此可以取值其最大允許誤差為±0.15，則其標(biāo)準(zhǔn)不確定度為

（3）A類測量不確定度評定

將前兩次測量值的均值作為測量結(jié)果時，則A類測量不確定度為u(x)=。結(jié)果如表2所示。

表2 1、2、3號灶熱負(fù)荷各輸入量的標(biāo)準(zhǔn)不確定度Table2 Standard uncertainty of heat input for 1，2 and 3 stoves

（4）分布傳遞的結(jié)果

由最大熵原理可以求得各個輸入量均為正態(tài)分布。當(dāng)把前兩次測量值的均值作為測量結(jié)果時，采用蒙特卡洛法傳遞概率分布[8-11]，利用MATLAB軟件進(jìn)行計算。1、2、3號灶的運行結(jié)果如表3所示，輸出量Y也即實測熱負(fù)荷的概率分布如圖3所示。取三個灶擴(kuò)展不確定度的平均值作為本套設(shè)備熱負(fù)荷測量不確定度的最終值，即為0.76%。

表3 1、2、3號灶熱負(fù)荷的不確定度Table3 Uncertainty of heat load for 1，2 and 3 stoves

1.3 家用燃?xì)庠顭嵝蕼y量不確定度評定示例

（1）測量原理和測量模型

熱效率是指系統(tǒng)有效輸出的能量與輸入的能量之比。同理對三個燃?xì)庠钸M(jìn)行熱效率試驗，測試步驟按照GB30720-2014《家用燃?xì)庠罹吣苄薅ㄖ导澳苄У燃墶分袩嵝试囼灧椒ㄟM(jìn)行。由式（3）計算實測熱效率。

M——實際加水量與鋁鍋換算為當(dāng)量加水量之和，單位為千克（kg）；

M1——加入鍋內(nèi)的水質(zhì)量，單位為千克（kg）；

M2——鋁鍋的質(zhì)量（含蓋子和攪拌器），單位為千克（kg）；

試驗完上限鍋和下限鍋的實測熱效率后用式（5）計算總熱效率。

η——灶頭總熱效率，%；

η實，下——使用下限鍋時的實測熱效率，%；

η實，上——使用上限鍋時的實測熱效率，%；

q下——使用下限鍋試驗時的鍋底熱強度，單位為瓦每平方厘米（W/cm2）；

q上——使用上限鍋試驗時的鍋底熱強度，單位為瓦每平方厘米（W/cm2）；

注：鍋底熱強度等于實測熱負(fù)荷（W）/試驗用鍋在正投影面的面積（cm2）。

根據(jù)式（5），熱效率η的輸入量有實測熱效率η實和鍋底熱強度q，而實測熱負(fù)荷又是鍋底熱強度q的輸入量。因此，可以分別評定實測熱效率和實測熱負(fù)荷的不確定度，以其作為輸入量，再評定燃燒器熱效率的測量不確定度，則式（5）就可以作為燃燒器熱效率的數(shù)學(xué)模型。

圖3 1、2、3號灶熱負(fù)荷概率分布圖Fig.3 Probability distribution of heat load for 1，2 and 3 stoves

（2）測量數(shù)據(jù)

同理對三臺灶的各個輸入量進(jìn)行10次獨立測量，計算得到輸入量的試驗標(biāo)準(zhǔn)差如表4和表5所示。

表4 三個灶的下限鍋熱效率各輸入量數(shù)據(jù)Table4 Theinput dataof thermal efficiency of lower limit

表5 三個灶的上限鍋熱效率各輸入量數(shù)據(jù)Table5 Theinput dataof thermal efficiency of upper limit

（3）A類不確定度評定

同理將前兩次測量值的平均值作為測量結(jié)果，可以得出三個灶的下限鍋和上限鍋熱效率各輸入量的標(biāo)準(zhǔn)不確定度。結(jié)果如表6和表7所示。

表6 下限鍋熱效率各輸入量的標(biāo)準(zhǔn)不確定度Table 6 The standard uncertainty of each input quantity of lower pot thermal efficiency

表7 上限鍋熱效率各輸入量的標(biāo)準(zhǔn)不確定度Table 7 The standard uncertainty of each input quantity of theupper pot heat efficiency

（4）分布傳遞結(jié)果

將兩次測量值的平均值作為測量結(jié)果，采用蒙特卡洛法傳遞概率分布，利用MATLAB軟件進(jìn)行計算，則三個灶的總的實測熱效率的概率分布如圖4所示，總熱效率的計算結(jié)果如表8所示。

表8 三個灶兩次測量的熱效率的不確定度Table8 Uncertainty in thermal efficiency of threefocifor twomeasurements

2 用GUM法評定測量不確定度

2.1 1號灶實測熱負(fù)荷的GUM法評定

根據(jù)實測熱負(fù)荷的模型和不確定度傳遞率，取兩次測量的平均值作為測量結(jié)果時，合成標(biāo)準(zhǔn)不確定度的公式如式（6）所示。

圖4 三個灶總實測熱效率概率分布圖Fig.4 Probability distribution of total measured thermal efficiency of threestoves

其中，

將數(shù)據(jù)代入式（6）得：

2.2 1號灶實測熱效率的GUM法評定

根據(jù)實測熱效率的模型和不確定度傳遞率，取兩次測量的平均值作為測量結(jié)果時，上限鍋合成標(biāo)準(zhǔn)不確定度的公式如式（7）所示：

其中，

代入式（7）得 uc(η實,上)=0.2653。

同理，下限鍋實測熱效率的合成標(biāo)準(zhǔn)不確定度公式為：

其中，

代入式 （8） 得 uc(η實,下)=0.261 0 。

于是，總的熱效率的合成標(biāo)準(zhǔn)不確定度公式為：

代入式（9）得uc(η)=0.1996，則其相對擴(kuò)展不確定度為×100%×2=0.66%，同理可以求得其他兩個灶的相對擴(kuò)展不確定度分別為0.83%、0.73%。

3 MCM和GUM的比較

綜上所述，熱負(fù)荷MCM和GUM計算的結(jié)果如表9所示；熱效率MCM和GUM計算的結(jié)果如表10所示。可以看出兩者結(jié)果是一致的。

表9 三個灶熱負(fù)荷MCM和GUM結(jié)果統(tǒng)計表Table9 Threeheat load MCMand GUMstatistical tables

表10 三個灶的總熱效率MCM和GUM結(jié)果統(tǒng)計表Table 10 The total thermal efficiency of three foci MCM and GUM resultstable

4 結(jié)論

用蒙特卡洛法進(jìn)行概率分布傳遞，最后實測熱負(fù)荷和總熱效率的模擬結(jié)果為正態(tài)分布。

用GUM法和MCM法兩種方法評定的熱負(fù)荷和熱效率的不確定度結(jié)果一致。

在GUM法的評定過程中可知試驗時間的測量不確定度對熱負(fù)荷測量結(jié)果不確定度的影響最大；當(dāng)量加水量的測量不確定度對熱效率測量結(jié)果不確定度的影響最大。

將三個灶用MCM法評定熱負(fù)荷的不確定度的平均值作為本套設(shè)備的最終不確定度，結(jié)果為0.76%，熱效率的最終不確定度為0.74%。