CN103123281A - 非接触式温度测量方法 - Google Patents

Abstract

一种非接触式温度测量方法，包括：拍摄目标点，以获取第一波长的图像数据以及第二波长的图像数据。其中在拍摄该目标点时，调整获取第一波长的图像数据时的快门时间以及获取第二波长的图像数据时的快门时间，以获得对应于第一波长的第一亮度以及对应于第二波长的第二亮度。以及根据下列方程式(1)来计算该目标点的温度：

其中λ₁为第一波长，λ₂为第二波长，L_λ1(T)为第一亮度，L_λ2(T)为第二亮度，A为校正系数，S₁为获取第一波长的图像数据时的快门时间，S₂为获取第二波长的图像数据时的快门时间，常数C＝hc/k，其中k为普朗克(Planck)常数，c为光速以及h为波尔兹曼(Boltzmann)常数。

Description

非接触式温度测量方法

技术领域

本发明是有关于一种温度测量方法，且特别是有关于一种适于测量燃烧炉温度场的非接触式温度测量方法。

背景技术

工业燃烧系统，例如是工业炼钢炉，火力电厂以及燃烧炉等。燃烧过程中常需监控燃烧场的炉炉内温度的分布，例如炉壁温度及钢胚温度等。其中测量炉壁温度是为避免因炉壁温度过高而造成炉体结构损坏，产生危险；而测量钢胚温度可预先判断产品品质，进而改变操作设定，提升产品良率。

燃烧场温度测量装置，分为接触式以及非接触式两种测量方式。其中，传统接触式的温度测量方式以高温热电偶为主。然而，接触式的高温热电偶其温度响应速度慢且仅能单点测量，对燃烧工艺的调整及监控的帮助有限。

此外，非接触式的测量方式常以可见光摄影机吸收目标物的辐射能，进而搭配双色法(two color method)推算出燃烧场的温度。图1为现有的黑体炉中，炉内温度与亮度的关系。图1中可得知当可见光摄影机快门固定在一速度100ms拍摄时，第一波长R在温度超过1273K时，亮度已经饱和。而第二波长G在温度低于1253K时，亮度趋近于0无法使用，造成可测量的温度范围非常小，实用性不佳。另外，也由于参考信息的不足，在燃烧系统的运转上不免造成能源浪费及污染性的增加，若能将燃烧炉内即时的温度变化显示予燃烧系统操作人员，将可提供直接且明确的参考。

发明内容

本发明提供一种具高测量范围的温度测量方法，以非接触式的感测技术，即时性地提供燃烧炉的温度场信息给燃烧系统操作人员。

本发明提供一种温度测量方法，包括：拍摄目标点，以获取第一波长的图像数据以及第二波长的图像数据。其中在拍摄该目标点时，调整获取第一波长的图像数据时的快门时间以及获取第二波长的图像数据时的快门时间，以获得对应于第一波长的第一亮度以及对应于第二波长的第二亮度。以及根据下列方程式(1)来计算该目标点的温度：

$T=\frac{C({{\mathrm{\λ}}_2}^{-1}-{{\mathrm{\λ}}_1}^{-1})}{\ln \frac{L_{{\mathrm{\λ}}_1}\left(T\right)/S_1}{L_{{\mathrm{\λ}}_2}\left(T\right)/S_2}+\ln A-5\ln \frac{{\mathrm{\λ}}_2}{{\mathrm{\λ}}_1}}---\left(1\right)$

其中λ₁为第一波长，λ₂为第二波长，L_λ1(T)为第一亮度，L_λ2(T)为第二亮度，A为校正系数，S₁为获取第一波长的图像数据时的快门时间，S₂为获取第二波长的图像数据时的快门时间，常数C＝hc/k，其中k为普朗克(Planck)常数，c为光速，h为波尔兹曼(Boltzmann)常数。

基于上述，本发明的温度测量方法，先利用快门调节的动作，使图像获取装置可测量的亮度区间增加，并对燃烧场进行图像获取。修正校正系数A之后，以改进式双色法计算得出燃烧场的炉内温度分布信息。

为让本发明的上述特征和优点能更明显易懂，下文特举实施例，并配合附图作详细说明如下。

附图说明

图1是现有的黑体炉中，燃烧温度与亮度的关系图。

图2为本发明的一实施例的温度测量方法的步骤图。

图3是于黑体炉中，快门时间与亮度的结果示意图。

图4是本发明的一实施例的快门调整方法的步骤图。

图5为本发明的一实施例的快门调整方法示意图。

图6为使用带通滤片之后，第一波长与第二波长的亮度与温度关系示意图。

主要元件符号说明

210～230：步骤

410～430：步骤

S_g：快门时间

S_t：快门时间

L_g：参考亮度

L_t：目标亮度

R：第一波长

G：第二波长

具体实施方式

图2为本发明的一实施例的温度测量方法的流程图。在本实施例中，其主要流程包含拍摄一目标点，以获取一第一波长R的图像数据以及一第二波长G的图像数据的步骤210。其中在拍摄该目标点时，还执行步骤220调整获取第一波长R的图像数据时的快门时间以及获取第二波长G的图像数据时的快门时间，以获得对应于第一波长R的一第一亮度以及对应于第二波长G的一第二亮度，最后在步骤230中，经由计算取得目标点的温度。

在本发明的一实施例中，调整获取第一波长R的图像数据时的快门时间以及获取第二波长G的图像数据时的快门时间的步骤，通过图3的测试结果进行调整。图3为利用一黑体炉作为测试目标，并设置一可见光摄影机，用以验证快门时间与单色亮度(以第一波长R为例)关系的结果示意图。由图3可知可见光摄影机在固定温度下的快门与单色亮度呈线性关系，因此可利用快门与亮度呈线性的关系，进而调整获取该第一波长R的图像数据时的快门时间以及获取该第二波长G的图像数据时的快门时间。

图4为本发明的一实施例的快门调整方法的步骤图。图5为基于快门与亮度呈线性关系的快门调整方法示意图。请参考图4与图5，在步骤410时，任意猜测的一测试快门时间S_g拍摄目标点，以取得一参考图像。参考图像对应于第一波长或第二波长的一参考亮度。在本实施例中，以参考图像对应于第二波长G取得参考亮度L_g。在步骤420，比对参考亮度L_g与一目标亮度L_t，以获得一目标快门时间S_t。以及在步骤430时，以目标快门时间S_t来获取第一波长R的图像数据或第二波长G的图像数据。此外，在进行步骤430之前，将会以步骤420所获得的目标快门时间作为另一测试快门时间，并重复步骤410至420至少一次。因此，通过调整快门速度，可使图像获取装置的曝光时间延长或缩短来取得适当亮度，故图1中的可测量的温度范围就可以增加。此外，图像获取装置更具有一光圈。本实施例在固定光圈下进行快门的调整。此外，也可通过调整光圈大小，使可测量的温度范围更大。

完成调整获取第一波长R以及第二波长G的图像数据时的快门时间后，即可采用调整后的快门速度来获取第一波长R的图像数据以及第二波长G的图像数据。

本实施例中可以采用多种可能的方案来获取第一波长R的图像数据以及第二波长G的图像数据。例如，可使用不同的两个图像获取装置分别获取第一波长R的图像数据以及第二波长G的图像数据。取得第一波长R及第二波长G的适当亮度，再将二张图像合成一张图像。

此外，也可只以同一个图像获取装置在不同的时序下获取第一波长R的图像数据以及第二波长G的图像数据。例如先获取第一波长R的图像数据之后，再获取第二波长G的图像数据，故为不同时间点的图像，再将二张图像合成一张图像。

另外，更可用同一个图像获取装置同时获取第一波长R的图像数据以及第二波长G的图像数据。其中，在图像获取装置与目标点之间更设置对应于第二波长G的一带通滤片，以减少图像获取装置获得的第一波长R的光通量。图6为使用第二波长G带通滤片之后，第一波长与第二波长的亮度与温度关系示意图。在图6中，因图像通过一带通滤片，在固定快门下第一波长R及第二波长G的亮度会相近，如图6中温度点1073K～1133K所示，在1153K以后加入快门调整动作，使第二波长G被固定在目标亮度的区域内，而第一波长R的亮度因温度增加，将随着快门时间的缩短而递减。使第一波长的图像数据以及第二波长的图像数据皆在可见光摄影机可拍摄的范围内。

通过获取第一波长R的图像数据以及第二波长G的图像数据，可取得第一波长R以及第二波长G所对应的第一亮度以及第二亮度。再利用第一亮度以及第二亮度计算目标点的温度。

本发明是使用双色法改进而来的计算式来计算目标点的温度。发展过程是由普朗克(Planck)在绝对黑体物体得到的单色辐射公式(普朗克法则，Planck’s Law)推导而得，如方程式(3)表示：

$I_{\mathrm{Planck}}(\mathrm{\λ},T)=\frac{2h{c}^2}{{\mathrm{\λ}}^5}{(e^{\frac{\mathrm{hc}}{\mathrm{\λkT}}}-1)}^{-1}---\left(3\right)$

其中，I_planck(λ，T)是指每单位时间内、每单位表面积、每单位立体角(solid angle)以及每单位波长下的释放能量(Js^-1m^-2sr^-1m^-1)。其中，λ为波长(m)，T为黑体的温度(K)，h为普朗克(Planck)常数(约等于6.62606896×10^-34)，c是光速(约等于2.99792458×10⁸ms^-1)，k为波尔兹曼(Boltzmann)常数(约等于1.3806504×10^-23JK^-1)，e是自然对数的基底。

对于燃煤的火焰辐射，其波长范围在300-1000nm及温度范围800-2000K，或是燃瓦斯的火焰辐射，其波长范围在400-700nm及温度范围在3000K以下，单色辐射公式(Planck’s Law)可用维恩定律(Wien’slaw)来表示，且由于炉膛环境非完全地的黑体，因此需将材料的放射率因子(emissivity)加入维恩定律(Wien’s law)，可以由方程式(4)表示：

$I_{\mathrm{Wien}}(\mathrm{\λ},T)=\mathrm{\ϵ}(\mathrm{\λ},T)\frac{2h{c}^2}{{\mathrm{\λ}}^5}{e}^{-\frac{\mathrm{hc}}{\mathrm{\λkT}}}=\mathrm{\ϵ}(\mathrm{\λ},T)\frac{2C_1}{{\mathrm{\λ}}^5}{e}^{-\frac{C}{\mathrm{\λT}}}---\left(4\right)$

在方程式(4)中，放射率因子(emissivity)的ε(λ，T)跟物体的温度以及辐射波长有关。C₁＝hc²＝0.59552138×10^-16(Wm²)，C＝hc/k＝1.43877516×10^-2(mK)。其中，h为普朗克(Planck)常数(约等于6.62606896×10^-34)，c是光速(约等于2.99792458×10⁸ms^-1)，k为波尔兹曼(Boltzmann)常数(约等于1.3806504×10^-23JK^-1)。若放射率因子(emissivity)不随波长改变则可称的为灰体(gray body)。光谱仪及摄影机的颜色强度通常会与来源的辐射强度线性相依，因此可做出线性响应的假设，同时考虑测量系统的光学和侦测器的效应影响，所以仪器记录的强度与外在来源辐射量的关系为方程式(5)。

$I({\mathrm{\λ}}_i,T)=k_i\·L_{{\mathrm{\λ}}_i}\left(T\right),i=r,g,b---\left(5\right)$

其中，k_i为灰阶尺度转换因子。L_λi(T)是可见光RGB颜色通道(channel)的强度值。

使用双色法的目的，主要是克服放射率因子(emissivity)的影响。因此假设在同一温度下，选择二相近波长的灰体假设成立，利用二个不同且接近波长的辐射能来计算温度，此二波长下的发射系数相比后可相消，故此方法称为双色法。基于同一目标在二个波长上的辐射发射强度，分别代入维恩定律(Wien’s law)后相比，即可由方程式(4)及方程式(6)推得双色法公式，如方程式(7)所示。

$\frac{I({\mathrm{\λ}}_1,T)}{I({\mathrm{\λ}}_2,T)}=\frac{k_1\·L_{{\mathrm{\λ}}_1}\left(T\right)}{k_2\·L_{{\mathrm{\λ}}_2}\left(T\right)}=A\frac{L_{{\mathrm{\λ}}_1}\left(T\right)}{L_{{\mathrm{\λ}}_2}\left(T\right)}---\left(6\right)$

$T=\frac{C({{\mathrm{\λ}}_2}^{-1}-{{\mathrm{\λ}}_1}^{-1})}{\ln \frac{L_{{\mathrm{\λ}}_1}\left(T\right)}{L_{{\mathrm{\λ}}_2}\left(T\right)}+\ln A-5\ln \frac{{\mathrm{\λ}}_2}{{\mathrm{\λ}}_1}}---\left(7\right)$

根据方程式(7)，系数A是唯一待决的参数，若系数A校正后，则方程式(7)就可作为温度场计算的理论依据。因为系数A具有方程式(7)的物理模型。因此若有已知的参考点，也就是T_ref、L_λ1(T_ref)、L_λ2(T_ref)已知，则系数A可利用通过下列方程式(8)求得。

$A=\frac{L_{{\mathrm{\λ}}_2}}{L_{{\mathrm{\λ}}_1}\left(T_{\mathrm{ref}}\right)}{\left(\frac{{\mathrm{\λ}}_2}{{\mathrm{\λ}}_1}\right)}^5\exp (C({{\mathrm{\λ}}_2}^{-1}-{{\mathrm{\λ}}_1}^{-1})/T_{\mathrm{ref}})---\left(8\right)$

此外，由前者黑体时测试可知，快门与亮度呈线性关系。因此在实作上分别调整第一波长R以及第二波长G的快门，以取得第一波长R以及第二波长G适当的亮度是可行的。取得单一波长在某温度下测得的亮度及快门，将亮度除以快门得到单位快门亮度值，此值即为快门与亮度线性关系的斜率值。斜率值的意义为单一波长在同温度下的任意快门与其对应亮度的比例是固定的，代表第一波长R与第二波长G在不同快门下得到的单位快门亮度，只要乘上同一快门，即可得到在同快门下的适当亮度值。因此以单位快门亮度替代原双色法公式，如方程式(7)中的亮度值，将可消除快门因子，解决第一波长R与第二波长G亮度是在不同快门条件下取得的问题，因此以单位快门下的亮度作为计算参考，方程式(7)可改写为方程式(1)所示。

$T=\frac{C({{\mathrm{\λ}}_2}^{-1}-{{\mathrm{\λ}}_1}^{-1})}{\ln \frac{L_{{\mathrm{\λ}}_1}\left(T\right)/S_1}{L_{{\mathrm{\λ}}_2}\left(T\right)/S_2}+\ln A-5\ln \frac{{\mathrm{\λ}}_2}{{\mathrm{\λ}}_1}}---\left(1\right)$

其中，S₁为L_λ1(T)在适当亮度时的快门时间，S₂为L_λ2(T)在适当亮度时的快门时间。因此校正系数A亦可改写为方程式(2)

$A=\frac{L_{{\mathrm{\λ}}_2}\left(T_{\mathrm{ref}}\right)/S_2}{L_{{\mathrm{\λ}}_1}\left(T_{\mathrm{ref}}\right)/S_1}{\left(\frac{{\mathrm{\λ}}_2}{{\mathrm{\λ}}_1}\right)}^5\exp (C({{\mathrm{\λ}}_2}^{-1}-{{\mathrm{\λ}}_1}^{-1})/T_{\mathrm{ref}})---\left(2\right)$

利用图像获取装置取得具有两种波长的图像，例如是以一带通滤片同时获取一具有第一波长R及第二波长G的图像。然而，经过一带通滤片滤波后，图像获取装置积分的波长范围仍有一定间距，无法取得单一波长下所对应的亮度。在双色法的应用上，波长的选择仍是一项问题。因此，本实施例还提出一种波长校正的步骤，在实际测量过程中，若有多点温度可供校正，利用已知的校正点信息，使用数值方法计算多点温度下最适当的二个波长值，将会取得更准确的温度值。

以下步骤为选取已知温度的多个校正点。然后分别拍摄这些校正点的图像，并选定图像中的一校正波长，以取得对应于校正波长的一第一校正亮度以及对应于一目标波长的一第二校正亮度。而本实施例中，可选定两个校正温度，代入方程式(2)。其中，λ₁为校正温度的校正波长，A为校正系数，λ₂为测量的目标温度的目标波长，T_ref为校正温度，C为一常数，L_λ1(T_ref)为校正波长的第一校正亮度，L_λ2(T_ref)为目标波长的第二校正亮度。因方程式(2)中，两个波长的快门时间S₁以及S₂、第一校正亮度L_λ1(T_ref)、第二校正亮度L_λ2(T_ref)、以及校正波长λ₁为已知。而目标波长λ₂以及校正系数A为未知数，所以需要两个校正点温度即可计算取得目标波长λ₂以及校正系数A。

取得校正系数A以及目标波长之后，最后再将校正系数A，以及先前拍摄第一波长以及第二波长的亮度代回方程式(1)，即可计算出温度T。

$T=\frac{C({{\mathrm{\λ}}_2}^{-1}-{{\mathrm{\λ}}_1}^{-1})}{\ln \frac{L_{{\mathrm{\λ}}_1}\left(T\right)/S_1}{L_{{\mathrm{\λ}}_2}\left(T\right)/S_2}+\ln A-5\ln \frac{{\mathrm{\λ}}_2}{{\mathrm{\λ}}_1}}---\left(1\right)$

其中λ₁为该第一波长R，λ₂为该第二波长G，L_λ1(T)为该第一亮度，L_λ2(T)为该第二亮度，A为校正系数，S₁为获取该第一波长的图像数据时的该快门时间，S₂为获取该第二波长的图像数据时的该快门时间，常数C＝hc/k，其中k为普朗克(Planck)常数，c为光速，h为波尔兹曼(Boltzmann)常数。

综上所述，本发明的温度测量方法，先利用快门调节的动作，使图像获取装置可测量的亮度区间增加，并利用图像获取设备以非接触的方式，对燃烧场进行图像获取。本发明还提出一波长选择的步骤，选定多个校正点，拍摄并计算校正系数A，再以改进式双色法计算得出燃烧场的炉内温度分布信息。此温度分布信息可提供给系统操作人员用以判断燃烧效率。

虽然本发明已以实施例揭露如上，然其并非用以限定本发明，任何本领域技术人员，在不脱离本发明的精神和范围内，当可作些许的更动与润饰，故本发明的保护范围当视权利要求书所界定者为准。

Claims (9)

1.一种温度测量方法，包括：

拍摄目标点，以获取一第一波长的图像数据以及一第二波长的图像数据，其中在拍摄该目标点时，调整获取该第一波长的图像数据时的快门时间以及获取该第二波长的图像数据时的快门时间，以获得对应于该第一波长的一第一亮度以及对应于该第二波长的一第二亮度；以及

根据下列方程式(1)来计算该目标点的温度：

$T=\frac{C({{\mathrm{\λ}}_2}^{-1}-{{\mathrm{\λ}}_1}^{-1})}{\ln \frac{L_{{\mathrm{\λ}}_1}\left(T\right)/S_1}{L_{{\mathrm{\λ}}_2}\left(T\right)/S_2}+\ln A-5\ln \frac{{\mathrm{\λ}}_2}{{\mathrm{\λ}}_1}}---\left(1\right)$

其中λ₁为该第一波长，λ₂为该第二波长，L_λ1(T)为该第一亮度，L_λ2(T)为该第二亮度，A为一校正系数，S₁为获取该第一波长的图像数据时的该快门时间，S₂为获取该第二波长的图像数据时的该快门时间，C为一常数，该常数C＝hc/k，其中k为普朗克常数，c为光速，h为波尔兹曼常数。

2.根据权利要求1所述的温度测量方法，其中调整获取该第一波长的图像数据时的快门时间或者调整获取该第二波长的图像数据时的快门时间的方法包括：

1)以一测试快门时间拍摄该目标点，以取得一参考图像，该参考图像对应于该第一波长或该第二波长的一参考亮度；

2)比对该参考亮度与一目标亮度，以获得一目标快门时间；以及

3)以该目标快门时间来获取该第一波长的图像数据或该第二波长的图像数据。

3.根据权利要求2所述的温度测量方法，其中调整获取该第一波长的图像数据时的快门时间或者调整获取该第二波长的图像数据时的快门时间的方法更包括在该步骤3)之前，进行下列步骤：

以该步骤2)所获得的该目标快门时间作为另一测试快门时间；以及

重复步骤1)-2)至少一次。

4.根据权利要求1所述的温度测量方法，其中分别以不同的两个图像获取装置来获取该第一波长的图像数据以及该第二波长的图像数据。

5.根据权利要求1所述的温度测量方法，其中以同一个图像获取装置在不同的时序下获取该第一波长的图像数据以及该第二波长的图像数据。

6.根据权利要求1所述的温度测量方法，其中以同一个图像获取装置同时获取该第一波长的图像数据以及该第二波长的图像数据。

7.根据权利要求6所述的温度测量方法，还包括：减少该图像获取装置所获得的该第一波长的光通量。

8.根据权利要求7所述的温度测量方法，其中减少该图像获取装置所获得的该第一波长的光通量的方法包括：

在该图像获取装置与该目标点之间设置对应于该第二波长的一带通滤片。

9.根据权利要求1所述的温度场测量方法，其中取得该校正系数A的方法包括：

提供已知温度的多个校正点；

分别拍摄所述校正点的图像，并选定该图像中的一校正波长，以取得对应于该校正波长的一第一校正亮度以及对应于一未知波长的一第二校正亮度；以及

依据下列方程式(2)来计算该校正系数A与该未知波长的关系：

$A=\frac{L_{{\mathrm{\λ}}_2}\left(T_{\mathrm{ref}}\right)/S_2}{L_{{\mathrm{\λ}}_1}\left(T_{\mathrm{ref}}\right)/S_1}{\left(\frac{{\mathrm{\λ}}_2}{{\mathrm{\λ}}_1}\right)}^5\exp (C({{\mathrm{\λ}}_2}^{-1}-{{\mathrm{\λ}}_1}^{-1})/T_{\mathrm{ref}})---\left(2\right)$

其中，T_ref为校正点温度，λ₁为该校正波长，λ₂为该未知波长，S₁为获取该第一波长的图像数据时的该快门时间，S₂为获取该第二波长的图像数据时的该快门时间，L_λ1(T_ref)为该第一校正亮度，L_λ2(T_ref)为该第二校正亮度，C为一常数，该常数C＝hc/k，其中k为普朗克常数，c为光速，h为波尔兹曼常数。

Images