文章信息
- 王玉刚, 吴军, 叶方平, 潘江, 张洪军
- WANG Yugang, WU Jun, YE Fangping, PAN Jiang, ZHANG Hongjun
- Rossini型气体热量计电校正实验与分析
- The experiment and analysis of electrical calibration in Rossini type gas calorimeter
- 中国测试, 2019, 45(2): 83-88
- CHINA MEASUREMENT & TEST, 2019, 45(2): 83-88
- http://dx.doi.org/10.11857/j.issn.1674-5124.2018020004
-
文章历史
- 收稿日期: 2018-02-20
- 收到修改稿日期: 2018-05-08
2. 计量技术实验教学中心,浙江 杭州 310018;
3. 浙江新劲空调设备有限公司,浙江 龙泉 323700
2. Experimental teaching center for metrology & measurement, Hangzhou 310018, China;
3. Zhejiang XinJing Air Conditioning Equipment Co.,Ltd., Longquan 323700, China
气体热量计有氧弹热量计、气体吸收式热量计、水流式热量计、Rossini型热量计,其中Rossini型气体热量计是目前测量准确度最高的气体热值测定装置,在近几年来具有相当高的研究价值。其基本工作原理是:将待测气体定压燃烧,燃烧释放的热量被吸热介质吸收,最后通过测量吸热介质的温升和参与燃烧气体的质量计算得到气体燃烧热值。Rossini型气体热量计[1]电校正实验的目的是测定量热容器的当量热容。20世纪初由美国的Rossini首次发明了纯甲烷气体的热值测定方法以及装置,几十年来相关学者在Rossini的研究基础上不断完善和优化实验装置,大大降低了该装置的测量不确定度。Pittam等[2]基于Rossini热量计的原理,设计了一套纯气体热值测量装置并做了相关实验,实验中所测量的纯气体热值与Rossini的数据相比偏差小于0.017%。Andrew Dale等[3]优化改进了气体样品质量的测定方法,使气体质量测量不确定度更加精确。M. Jaeschke等[4]提出温度和质量是实验总不确定度的主要影响因素,并在此理论基础上研究了一款不确定度可达到0.01%的自动测量与校准装置。F. Haloua等[5]研究出一套等环境气体热量计,并通过实验测得燃烧中大约有90%的热量被吸热介质吸收,只有约10%的热量被周围量热容器和其他辅助设备所吸收。P. Schley等[6]设计了一款基准气体热量计,该装置通过德国物理技术研究院完成测试,测试结果显示该热量计测量纯甲烷的热值不确定度可达0.05%。我国对气体热量计的研究相比国外起步较晚,胡日恒等[7]研究者在20世纪80年代研制了我国第一套精密氧弹热量计。此后,我国还研制了水流式热量计,测定纯甲烷的不确定度为1.0%(k=2)[8]。近年来中国计量科学院改进优化了氧弹热量计,使该热量计的不确定度可以达到0.6%[9]。最近,中国石油西南油田分公司天然气研究院依据Rossini燃烧测量原理,正在研制一套不确定度可达0.25%的气体热值测定装置[10]。
以上研究表明,在测量气体热值方面,提高Rossini型气体热值测量装置的测量不确定度具有重要意义。与德国、法国等欧洲国家相比,我国在气体热值测量方面的研究还存在较大的不足。基于以上背景,本文按照自行设计的一套Rossini型气体热量计展开实验。在改进和优化该装置后做了电校正实验,根据实验中所测得的数据分析该装置的测量不确定度,为以后的实验研究做好铺垫。本实验在一定程度上填补了我国在气体热值计量方面的研究空白。
1 实验装置搭建Rossini型气体热量计的电校正实验示意图如图1所示。整套实验装置放置在恒温间内,燃烧器外壁面缠有电加热带,整个燃烧器位于量热容器的吸热介质(水浴)中,而量热容器又被悬挂浸没在恒温槽中。系统中制冷模块主要功能是调节量热容器内吸热介质的初始温度,与制冷模块相连的是U型管换热器,其作用是给吸热介质降温。其中制冷模块通过硅胶管与U型管换热器连接,硅胶管中流有循环水。实验中为维持恒温间内温度的稳定而设计了风道,通过风道可将恒温槽工作时产生的热风排到恒温间外环境中。图2为电加热电路模块实物连接图。实验中使用Fluke 8508A型参考级81/2位数字多用表和计算机监测、采集实验数据,用Agilent 6634B可编程直流电源为电加热带提供恒定电压,用半导体预冷单元来调节量热容器内吸热介质的温度。表1为上述所选仪器在特定量程内的使用参数。
所需仪器 | 仪器类型 | 仪器参数 |
温度计 | 标准铂电阻温度计 | U=1.5 mK,k=2 |
数字万用表 | Fluke 8508A | 量程2 V,U=2.9×10–6 V,k=2; 量程200 V,U=4.2×10–6 V,k=2 |
直流电源 | Agilent 6634B | 量程100 V,U=50 mV,k=2 |
2 实 验 2.1 实验原理
当量热容采用电校正方法[11]测定,其基本测定原理为:在电加热带两端给予恒定的电压,根据电加热带消耗的电能和吸热介质温升计算得到装置的当量热容,具体计算过程为
${{{C}}_{\rm{eq}}} = \frac{{{{{E}}_{\rm{eq}}}}}{{\Delta {{{T}}_{\rm{eq}}}}}$ | (1) |
式中:Ceq——当量热容,J/K;
Eeq——电加热能量,J;
ΔTeq——吸热介质的温升,K。
${{{E}}_{\rm{eq}}} = \sum\limits_i {\frac{{{U_{{\rm{heat,}}i}}{U_{{\rm{ref,}}i}}}}{{{R_{{\rm{ref,}}i}}}}} \Delta {{{t}}_i}$ | (2) |
式中:Uheat,i——电加热带两端的电压,V;
Uref,i——参考电阻两端的电压,V;
Rref,i——参考电阻的阻值,Ω;
Δti——加热时间,s。
2.2 实验环境温度控制 2.2.1 恒温间内环境温度控制本实验对于温度环境的要求非常高,要求环境温度维持在(25±0.05)℃,因此专门搭建了一个恒温间。恒温间由铝合金型材和铝塑板搭建而成,长3 m、宽2.5 m、高2.5 m。如图3所示,通过恒温间的墙体来削弱和延迟外环境温度的变化,增加恒温间内温度的稳定性[12]。
与此同时还在恒温槽出风口设计了一个排风风道,可将恒温槽工作中冷凝器产生的热风排到外环境中,尽量减少恒温室中的热源。而且在恒温间内还安装了两台风扇,加速恒温间内空气的流动,保证恒温间内温度均匀。恒温间内部与外界主要通过顶部的铝塑板与外界进行换热。为了时刻监测恒温间内的温度变化,还在恒温间内放置了6个温度传感器,确保将环境的影响因素做到最小。
2.2.2 恒温间外环境温度控制为了更精确地控制恒温间内环境温度的稳定,专门设计了恒温间外环境控制模块。恒温间内仪器工作会产生热负荷,这部分热负荷需要通过热传递排到外环境中。故实验中在外环境安装了佳力图机房专用空调,空调采用顶棚孔板式送风。恒温间外环境的温度在20~30℃,控制精度可达±0.5℃,相对湿度控制可达±2%。
实验时为了达到仪器工作产生的热负荷和空调制冷量相对平衡。需要先打开实验仪器和空调并稳定工作12 h,直到恒温间内温度达到稳定值后开始实验。
温度控制结果:外环境温度控制可达(23±0.5)℃,内环境温度控制可达(25±0.05)℃。根据传热设计,内外环境的2℃温差刚好将恒温室内的实验仪器产生的热负荷带走。
2.3 实验过程实验开始前,先连接好所有仪器并检查连接是否正常。接着打开恒温槽,温度设定25℃;打开空调,温度设定23℃;打开量热容器上的搅拌电机,根据前期数值模拟的结果,将转速设定为420 r/min[13]。经过12 h后,当恒温间内环境温度稳定在25℃时,开始电校正实验。
具体实验过程为:1)打开数据采集系统和制冷系统,通过U型管换热器将热量计内吸热介质的温度降至23℃;2)关闭制冷系统,此时热量计内吸热介质通过量热容器壳体与恒温槽水浴缓慢换热,换热时间设定为20 min;3)打开恒流源,在电加热带两端加恒定电压工作20 min;4)关闭恒流源,此时热量计内吸热介质通过量热容器壳体与恒温槽水浴缓慢换热,换热时间20 min;5)保存数据重复以上步骤,记录下9组数据。
3 实验结果与分析 3.1 绝热温升ΔTad的测定如图4所示为电校正实验吸热介质时间-温度曲线。当量热容器内吸热介质温度冷却到23 ℃时开始实验。由于点火前量热容器内吸热介质与恒温槽水浴会缓慢地换热,此时吸热介质温度上升很慢,20 min温升仅为0.2 ℃,这个阶段被称为初始阶段。在初始阶段结束后,电加热带开始工作,吸热介质由于电加热带释放热量,此时温度上升很快,20 min温升达3.8 ℃,这个阶段被称为主阶段。在电加热带停止加热后,吸热介质的温度由于传热还略有上升。在热量完全释放和分配之后,结束阶段开始,此时吸热介质和恒温槽水浴进行缓慢地换热,温度曲线平缓下降。
在初始阶段吸热介质从恒温槽中吸收热量,在结束阶段吸热介质向恒温槽中释放热量,这两部分热量刚好抵消,对实验结果不产生影响。
吸热介质的绝热温升是对环境温度修正后的值,该温升可由下式计算得到:
$\Delta {{{T}}_{\rm{ad}}} = {{{T}}_{\rm{e}}} - {{{T}}_{\rm{b}}} - \delta {{T}}$ | (3) |
式中:ΔTad——吸热介质绝热温升,℃;
Tb——开始加热时的温度,℃;
Te——结束加热时的温度,℃;
热量计在初始阶段和结束阶段的温度变化率可根据牛顿冷却定律推导得到:
$\frac{{{\rm{d}}T}}{{{\rm{d}}t}} = u(t) + {{k}}({T_{\rm{0}}}{\rm{ - }}T(t))$ | (4) |
式中:T(t)——吸热介质的平均温度,℃;
T0——恒温槽水浴温度,℃;
u(t)——温度随时间变化速率,℃/s;
t——时间,s;
k——冷却常数。
当热量计稳定一段时间后,吸热介质将达到平衡温度T∞,该温度高于恒温槽水浴的温度。在这种情况下可近似看作dT/dt=0,根据式(4)可得:
$u(t) = {{k}}({T_\infty }{\rm{ - }}{T_{\rm{0}}})$ | (5) |
${\rm{d}}T/{\rm{d}}t = {{k}}({T_\infty }{\rm{ - }}T(t))$ | (6) |
故初始阶段和结束阶段都可以用指数函数表示为
${T_i}(t) = {T_\infty } - ({T_\infty } - {T_{\rm{b}}}) \cdot \exp ( - {{k}}(t - {t_{\rm{b}}}))$ | (7) |
式中:tb——电加热开始时间,min;
Tb——该时刻吸热介质温度值,℃;
k——冷却常数;
T∞——平衡温度,℃。
其中冷却常数和平衡温度均可通过吸热介质温升曲线的指数回归分析得到。主阶段的吸热介质绝热温度修正项
$\delta T = {{k}}\int_{{t_{\rm{b}}}}^{{t_{\rm{e}}}} {({T_\infty }{\rm{ - }}T(t))} {\rm{d}}t$ | (8) |
本实验为了避免偶然性和粗大误差,做了9次重复性实验,对保存的实验数据处理后,得到实验结果如表2所示。从9组实验中可以发现:吸热介质的平均温升3.863℃,测得的平均当量热容为19 023 J/K。
序号 | 电加热能量Eeq/J | 温升ΔTeq/K | 当量热容Ceq/(J·K–1) |
1 | 73 485 | 3.850 | 19 087 |
2 | 73 498 | 3.854 | 19 069 |
3 | 73 459 | 3.830 | 19 179 |
4 | 73 488 | 3.841 | 19 137 |
5 | 73 472 | 3.861 | 19 034 |
6 | 73 495 | 3.903 | 18 844 |
7 | 73 479 | 3.870 | 18 987 |
8 | 73 465 | 3.909 | 18 789 |
9 | 73 469 | 3.849 | 19 082 |
如图5所示为电加热实验和燃烧实验中热量计主体内吸热介质温度比较图。其中横坐标为时间,纵坐标为电加热与燃烧实验吸热介质温度的差值。由图可知:两种实验吸热介质的温度差小于6 mK,说明电加热实验与燃烧实验吸热介质温升曲线吻合较好。曲线存在明显的下降趋势,主要原因是阶段一中并未电加热或燃烧,吸热介质温升数值较小,此时环境温度影响很大;阶段二中开始电加热或燃烧,此时吸热介质温升数值增大,环境温度对其影响越来越小;阶段三中曲线又有略微上升的趋势,是因为此时电加热或燃烧结束,吸热介质向周围环境放热后温度会略微下降,此时环境温度对其影响又开始慢慢变大。
4 不确定度分析
气体热量计电校正实验测得的当量热容,其不确定度主要与电加热不确定度和绝热温升不确定度相关。根据文献[14]中铂电阻的标准不确定度评定方法,只要计算出这两部分的不确定度,就可以计算出当量热容的不确定度。
4.1 电加热引入的不确定度分析根据电加热能量计算公式和合成标准不确定度的方法,推导出不确定度分量u(Bi)和电加热能量不确定度u(Eeq),具体计算过程为
$u({B_i}) = \sqrt {{{({{\bar U}_{{\rm{ref}}}} \cdot \overline {\Delta t} )}^2}{u^2}({U_{{\rm{heat}},i}}) + {{({{\bar U}_{{\rm{heat}}}} \cdot \overline {\Delta t} )}^2}{u^2}({U_{{\rm{ref}},i}}) + {{({{\bar U}_{{\rm{heat}}}} \cdot {{\bar U}_{{\rm{ref}}}})}^2}{u^2}(\Delta {t_i})} $ | (9) |
$u({E_{{\rm{eq}}}}) = \sqrt {n{{\left(\frac{1}{{{R_{{\rm{ref}}}}}}u({B_i})\right)}^2} + {{\left(\frac{{\displaystyle\sum\limits_1^n {{B_i}u({R_{{\rm{ref}}}})} }}{{R_{{\rm{ref}}}^2}}\right)}^2} + \frac{{n(n - 1)}}{{R_{{\rm{ref}}}^2}}u({B_i},{B_j})} $ | (10) |
电加热能量含多个不确定影响因素。实验中涉及到的重要不确定度分量如表3所示。
不确定度分量 | 不确定度评定方法 | 标准不确定度 |
u1 | A类不确定度评定 | 1.072×10–3 V |
u2 | B类不确定度评定 | 3.941×10–4 V |
u3 | A类不确定度评定 | 1.436×10–5 V |
u4 | B类不确定度评定 | 2.734×10–6 V |
表3中u1指电加热带两端电压重复测量引入的不确定度;u2指使Fluke 8508A测量电加热带两端电压,其示值误差引入的不确定度;u3指参考电阻两端电压重复测量引入的不确定度;u4指Fluke 8508A测量参考电阻两端电压,其示值误差引入的不确定度。
根据表3中的不确定度数据计算得到电加热带的不确定度u(Uheat,i)=1.100×10–3 J;参考电阻的不确定度u(Uref,i)=1.462×10–6 Ω;时间测量不确定度u(Δti)=3.5×10–4 s。最后将测量值
从绝热温升的计算公式中可发现绝热温升的不确定度主要影响因素有:温度修正项
${{u}}(\Delta {{{T}}_{\rm{eq}}}) = \sqrt {{{{u}}^2}({{{T}}_{\rm{e}}}) + {{{u}}^2}({{{T}}_{\rm{b}}}) + {{{u}}^2}(\delta {{T}}) - 2{{u}}({{{T}}_{\rm{e}}},{{{T}}_{\rm{b}}}) - 2{{u}}({{{T}}_{\rm{e}}},\delta {{T}}) + 2{{u}}({{{T}}_{\rm{b}}},\delta {{T}})} $ | (11) |
绝热温升包含的不确定影响因素如表4所示。表中u5指温度Tb重复测量引入的不确定度;u6指Fluke 8508A测量温度Tb时,其示值误差引入的不确定度;u7指吸热介质温度均匀性引入的不确定度;u8指温度Te重复测量引入的不确定度;u9指Fluke 8508A测量温度Te时,其示值误差引入的不确定度。
不确定度分量 | 不确定度评定方法 | 标准不确定度 |
u5 | A类不确定度评定 | 3.142×10–3 ℃ |
u6 | B类不确定度评定 | 2.551×10–4 ℃ |
u7 | B类不确定度评定 | 2.887×10–3 ℃ |
u8 | A类不确定度评定 | 2.873×10–3℃ |
u9 | B类不确定度评定 | 2.574×10–4 ℃ |
根据不确定度计算方法计算得到u(Tb)=4.275×10–3 ℃,u(Te)=4.081×10–3 ℃。式(11)中第3项与前两项相比可忽略不计且3个参数之间相关系数几乎为零。将u(Tb)和u(Te)代入式(11)可以计算得到u(ΔTeq)=5.91×10–3 ℃。
4.3 当量热容不确定度分量分析当量热容合成不确定度为
${u^2}({{{C}}_{\rm{eq}}}) = {\left(\frac{1}{{\Delta {{{T}}_{\rm{eq}}}}}\right)^2}{{{u}}^2}({{{E}}_{\rm{eq}}}) + {\left(\frac{{{{{E}}_{\rm{eq}}}}}{{\Delta {{{T}}_{\rm{eq}}}^2}}\right)^2}{{{u}}^2}({{{T}}_{\rm{eq}}})$ | (12) |
其中:
$\begin{split} {{{u}}^2}({{{E}}_{\rm{eq}}}) = {\left(\frac{{{U_{{\rm{ref,}}i}}\Delta {{{t}}_i}}}{{{R_{{\rm{ref,}}i}}}}\right)^2}{{{u}}^2}({U_{{\rm{heat,}}i}}) +{\left(\frac{{{U_{{\rm{heat,}}i}}\Delta {{{t}}_i}}}{{{R_{{\rm{ref,}}i}}}}\right)^2}{{{u}}^2}({U_{{\rm{ref,}}i}}) + {\left(\frac{{{U_{{\rm{heat,}}i}}{U_{{\rm{ref,}}i}}\Delta {{{t}}_i}}}{{{R_{{\rm{ref,}}i}}^2}}\right)^2}{{{u}}^2}({{{R}}_{{\rm{ref}},i}}) \end{split} $ |
${{u}}(\Delta {{{T}}_{\rm{eq}}}) = \sqrt {{{{u}}^2}({{{T}}_{\rm{e}}}) + {{{u}}^2}({{{T}}_{\rm{b}}}) + {{{u}}^2}(\delta {{T}}) - 2{{u}}({{{T}}_{\rm{e}}},{{{T}}_{\rm{b}}}) - 2{{u}}({{{T}}_{\rm{e}}},\delta {{T}}) + 2{{u}}({{{T}}_{\rm{b}}},\delta {{T}})} $ |
综上所述,可以计算得到u(Ceq)=28 J/K,这里包含因子k取2,则当量热容的扩展不确定度为56 J/K,所得的当量热容则可表示为(19 023±56) J/K(k=2)。
为了更明确分析和对比各不确定度分量,引进了标准不确定度和相对不确定度,结果如表5所示。
不确定度分量 | 实验所测平均值 | 标准不确定度 | 相对不确定度 |
u(Eeq) | 73 479 J | 104 J | 0.14% |
u(ΔTeq) | 3.864 K | 5.91 mK | 0.15% |
u(Ceq) | 19 023 J/K | 28 J/K | 0.15% |
从表5中可发现:电加热能量和绝热温升的相对不确定度差不多,但绝热温升的相对不确定度相比之下较容易降低。绝热温升不确定度较大的原因是:热量计内吸热介质在空间上温度分布不够均匀,造成温度多次测量中数据偏差较大,引入的重复性测量不确定度较大。在今后的实验中会重点考虑该影响因素,提高装置的测量准确度。
5 结束语本文应用Rossini型气体热量计完成了电校正实验。实验中包括实验装置的搭建、实验环境温度的控制、实验方案设计、实验数据处理与分析。电校正实验当量热容测量不确定度为0.15%,不确定度较大的主要原因是测量仪器准确度还有欠缺,在今后的研究中,将考虑引进准确度更高的仪器。本次实验填补了我国在气体热值计量研究领域的部分空白,为之后气体热值测量的研究提供了实验和数据参考。
[1] |
ROSSIINI F D. The heat of formation of water[J].
Proceedings of the National Academy of Sciences, 1930, 16(11): 694-699.
DOI:10.1073/pnas.16.11.694 |
[2] |
PITTAM D A, PELCHER G. Measurements of heats of combustion by flame calorimetry. part 8.-methane, ethane, propane, n-butane and 2-methylpropane[J].
Journal of the Chemical Society, Faraday Transaction1: Physical Chemis-try in Condensed Phases Transactions, 1972(68): 2224-2229.
|
[3] |
DALE A, LYTHALLl C, AUCOTT J, et al. High precision calorimetry to determine the enthalpy of combustion of methane[J].
Thermochimica Acta, 2002, 382(1): 47-54.
|
[4] |
JAESCHKE M, BENITO S A, CREMONESI P L, et al. GERG Project: Development and set-up of a new reference calorimeter[J].
Gartenbauwissens-chaft, 2004, 54(4): 179-184.
|
[5] |
HALOUA F, HAY B, FILTZ J R. New French reference calorimeter for gas calorific value measurements[J].
Journal of Thermal Analysis and Calorimetry, 2009, 97(2): 673-678.
DOI:10.1007/s10973-008-9701-z |
[6] |
SCHLEY P, BECK M, UHRIG M, et al. Measurements of the calorific value of methane with the new GERG reference calorimeter[J].
International Journal of Thermophysics, 2010, 31(4): 665-679.
|
[7] |
胡日恒, 安绪武, 谈夫. 精密氧弹量热计及苯甲酸燃烧热的测定[J].
化学学报, 1981(S1): 18-26.
|
[8] |
苏毅, 骆科东. 我国天然气流量计量技术现状及发展趋势分析[J].
石油仪器, 2015(5): 8-12.
DOI:10.3969/j.issn.1004-9134.2015.05.002 |
[9] |
王海峰, 李佳, 孙国华, 等. 基准气体热量计研究进展[J].
石油与天然气化工, 2014(2): 196-199.
DOI:10.3969/j.issn.1007-3426.2014.02.019 |
[10] |
李克, 潘春锋, 张宇, 等. 天然气发热量直接测量及赋值技术[J].
石油与天然气化工, 2013, 42(3): 297-301.
DOI:10.3969/j.issn.1007-3426.2013.03.021 |
[11] |
HALOUA F, HAY B, FOULON E. Uncertainty analysis of theoretical methods for adiabatic temperature rise determination of methane[J].
Therm Anal Calorim, 2013, 111: 985-994.
DOI:10.1007/s10973-012-2342-2 |
[12] |
廖绍凯, 牛志荣, 庄伟晨. 保温层对墙体温度场的延迟和衰减[J].
建筑节能, 2013(5): 41-44.
DOI:10.3969/j.issn.1673-7237.2013.05.010 |
[13] |
杜军燕, 王玉刚, 潘江, 等. Rossini型气体热量计非稳态温度场的优化分析[J].
中国测试, 2017, 43(8): 119-124.
|
[14] |
陈桂生, 付志勇, 朱育红, 等. 工业铂热电阻不同检定方法检定结果可信度研究[J].
中国测试, 2014, 40(6): 17-21.
|