需要JoVE订阅才能观看此内容。 请登录或开始免费试用

方法文章

近红外测温技术用于测量感应加热小型磁性球体周围水的温度

8.1K 次观看

⸱

DOI：

10.3791/57407

⸱

2018年4月30日

 ,  ,  ,  ,  , 
1Department of Mechanical Engineering, Tokyo Metropolitan University, 2Department of Electrical and Electronic Engineering, Tottori University, 3Brain Science Inspired Life Support Research Center, The University of Electro-Communications

本文内容

勘误通知

Important: There has been an erratum issued for this article. View Erratum Notice

摘要

介绍了一种利用1150 nm和1412 nm波长测量感应加热下小磁性球体周围水温的技术。

摘要

本文介绍了一种用于测量感应加热的小型磁性球体周围水及非浑浊水相介质温度的技术。该技术利用了水在1150 nm和1412 nm波长处的吸收系数随温度变化的特性。通过窄带通滤光片选择1150 nm或1412 nm的入射光，照射含有直径为2.0 mm或0.5 mm磁性球体的水或非浑浊水凝胶；同时，利用近红外相机获取吸收系数的横向投影——二维吸光度图像。当三维温度分布可假设为球对称时，可通过将逆阿贝尔变换应用于吸光度分布来估算温度分布。实验观察到温度随时间及感应加热功率的变化呈现一致的规律性。

引言

在许多科学研究领域和应用中，都需要一种用于测量介质内靠近小型热源处温度的技术。例如，在磁热疗研究中——这是一种利用磁性颗粒或微小磁性元件的电磁感应进行癌症治疗的方法——准确预测由磁性颗粒产生的温度分布至关重要1,2。然而，尽管微波3,4、超声5,6,7,8、光声9、拉曼10以及基于磁共振11,12的温度测量技术已被研究和开发，目前仍无法准确测量此类内部温度分布。迄今为止，人们主要通过温度传感器测量单点位置或少数几个位置的温度，在感应加热的情况下，通常使用非磁性的光纤温度传感器13,14。另一种方法是利用红外辐射测温仪远程测量介质表面温度，以估算其内部温度14。然而，当含有小型热源的介质为水层或非浑浊的水相介质时，我们已证明近红外（NIR）吸收技术可用于有效测量温度15,16,17,18,19。本文将详细介绍该技术的操作流程及典型实验结果。

近红外（NIR）吸收技术基于水在近红外区域吸收带随温度变化的原理。如图1a所示，水的ν1 + ν2 + ν3吸收带出现在波长（λ）1100 nm至1250 nm范围内，并随着温度升高向更短波长方向移动19。其中，ν1 + ν2 + ν3表示该吸收带对应于三种基本O-H振动模式的组合：对称伸缩振动（ν1）、弯曲振动（ν2）和反对称伸缩振动（ν3）20,21。光谱的这一变化表明，该吸收带中对温度最敏感的波长为λ ≈ 1150 nm。水的其他吸收带也表现出类似的温度依赖性15,16,17,18,20,21。在图1b中展示了波长范围λ = 1350–1500 nm内观测到的水的ν1 + ν3吸收带及其温度依赖性。在水的ν1 + ν3吸收带中，1412 nm是对温度最敏感的波长。因此，可通过使用近红外相机在λ = 1150 nm或1412 nm处采集二维（2D）吸光度图像，从而获得二维温度分布图像。由于水在λ = 1150 nm处的吸收系数小于在λ = 1412 nm处的吸收系数，前者适用于约10 mm厚的水相介质，而后者适用于约1 mm厚的水相介质。最近，我们利用λ = 1150 nm获得了含有一个感应加热的直径1 mm钢球的10 mm厚水层中的温度分布19。此外，已有研究利用λ = 1412 nm测量了0.5 mm厚水层中的温度分布15,17。

基于近红外（NIR）的温度成像技术的一个优势是，其设置和实施较为简便，因为该技术属于透射-吸收测量方法，无需使用荧光染料、磷光体或其他热探针。此外，其温度分辨率达到0.2 K以下15,17,19。其他基于干涉测量法的透射技术通常用于传热与传质研究22,23,24，但难以实现如此高的温度分辨率。然而需要指出的是，当存在显著的局部温度变化时，基于近红外的温度成像技术并不适用，因为较大的温度梯度会引起明显的光偏折效应19。本文将就该问题在实际应用中的影响进行讨论。

本文描述了基于近红外（NIR）的温度成像技术用于感应加热的小型磁性球体的实验装置与操作流程；此外，还展示了两张具有代表性的二维（2D）吸光度图像结果。第一张图像为直径2.0 mm的钢球置于厚度为10.0 mm的水层中，在波长λ = 1150 nm处采集得到。第二张图像为直径0.5 mm的钢球置于厚度为2.0 mm的麦芽糖浆层中，在波长λ = 1412 nm处采集得到。本文还介绍了通过将逆阿贝尔变换（IAT）应用于二维吸光度图像，计算三维（3D）径向温度分布的方法及结果。当假设三维温度分布具有球对称性时（如加热球体的情况，参见图2），逆阿贝尔变换是适用的19。在本研究中，IAT计算采用了多高斯函数拟合方法，因为高斯函数的逆阿贝尔变换可解析求解25,26,27,28,29，且能很好地拟合单调递减的数据，包括单个热源热传导实验中的数据。

访问受限。请登录或开始试用以查看此内容。

方案

1. 实验设置与操作步骤

按照以下步骤准备用于近红外成像的光学导轨，以安装样品和光学元件。

  1. 样品制备。
    注意：使用水或水溶液时，请执行步骤 1.1.1；使用高黏度水凝胶时，请执行步骤 1.1.2。
    1. 水中钢球的设置
      1. 使用少量胶水将一个直径为 2.0 mm 的钢球固定在细塑料绳的末端。
      2. 将钢球悬挂在矩形玻璃池的中心位置，该池的光程长度为 10.0 mm，宽度为 10 mm，高度为 45 mm（图 3）。
      3. 小心地将过滤后的水倒入池中，避免产生气泡。
        ​注意：也可使用少量胶水将钢球固定在细塑料棒的尖端19。
    2. 水凝胶中钢球的设置
      1. 加热水凝胶以降低其黏度，使其足够低以便顺利倾倒。
      2. 使用注射器将水凝胶注入光程长度为 2.0 mm、宽度为 10 mm、高度为 45 mm 的矩形玻璃池至半满，并静置冷却。
      3. 将一个直径为 0.5 mm 的钢球置于凝胶表面的中心位置。
      4. 用水凝胶填满玻璃池。
        注意：不应在凝胶中使用较大的钢球（直径 >~1 mm），因为它们在感应加热过程中会受到重力和/或磁力作用而移动。
    3. 将玻璃池放入塑料支架中，并将其安装在光学导轨上（图 3）。
  2. 近红外成像系统的准备。
    1. 准备一个带有光纤导光管的卤素灯，并使用支架将光纤导光管的末端固定在光学导轨上。
    2. 在光纤导光管与样品池之间放置一个窄带通滤光片（NBPF），其透射峰值位于 λ = 1150 nm 或 λ = 1412 nm（图 3）。
    3. 在卤素灯与窄带通滤光片（NBPF）之间插入另一个带通滤光片（BPF），其透射波长范围应比 NBPF 更宽。
      注意：由于 BPF 直接接收强光，因此需要该滤光片以防止 NBPF 因热损伤而损坏。
    4. 在 NBPF 与样品池支架之间的光路中插入一个或多个光圈，以减少杂散光（图 3）。
    5. 设置一台近红外相机，用于检测透过样品池的光（图 3）。通过数据传输线将相机连接至安装在个人计算机（PC）中的图像采集软件所对应的图像采集卡。
    6. 在样品池与相机之间安装一个远心镜头（图 3）。
      注意：也可使用普通相机镜头，但远心镜头在选择性检测与主光线平行的光（用于 IAT）以及减小衍射影响方面性能更优。
      注意：不应将 NBPF 和 BPF 放置在样品池与相机之间，否则水会因直接吸收卤素灯的高强度光而导致温度升高。
    7. 打开近红外相机，并启动图像采集软件。
    8. 开启卤素灯，并通过观察显示器上的图像调节其输出功率（图 4）。
    9. 调节远心镜头的轴线、位置和焦距，以获得清晰的钢球图像。
      注意：若调节不充分，将出现不规则的强度分布，导致吸光度计算错误。
  3. 感应加热系统的准备。
    1. 准备一套感应加热系统，包括高频发生器（最大输出功率：5.6 kW；频率：780 kHz）、水冷线圈和水冷机。
      注意：适用于小型金属部件钎焊、焊接和锡焊的感应加热系统即可满足本实验需求；参见 材料表。
    2. 如有可能，将线圈安装在可三维移动的平台上，以便调节其位置。
    3. 将线圈放置在样品池附近，使线圈中心与钢球之间的距离约为 15 mm（图 3）。确保线圈周围无其他金属部件。
      注意：距离应根据感应加热功率和钢球尺寸进行适当调整。
    4. 启动冷却水循环。
  4. 图像采集与感应加热。
    1. 在图像采集软件中点击“开始”，以连续存储图像。
    2. 在感应加热控制软件中点击“开始”，启动感应加热。
    3. 数秒后（根据实验条件和目的而定），在图像采集软件中点击“停止”。
    4. 在感应加热控制软件中点击“停止”。
    5. 在图像采集软件中将临时存储的图像保存为 TIFF 序列（或其他无压缩格式）。
      注意：若温度足够高，图像上将出现光偏转效应7。必须通过实验适当降低感应加热功率，使钢球附近温度升高幅度小于约 10 K，该结果可在后续温度估算步骤中验证。

2. 图像处理与温度估算

注意：保存的序列图像表示为 Ii(x, z)，其中 i 为序列帧编号。x、y、z、r 和 r' 等坐标定义如图2所示；z 轴正方向与重力方向相反。以下实验步骤的概要也在补充材料1中进行了图示说明。

  1. 吸光度图像构建。
    1. 使用图像处理软件打开 Ii(x, z)。
    2. 通过实施 3 × 3 像素平均法降低 Ii(x, z) 中的噪声。
    3. 在加热前对 i = 1 至 5（或更多）的 Ii(x, z) 创建平均图像，并将其定义为参考图像 Ir(x, z)。
      注：该平均处理可降低噪声，从而获得比单帧图像更可靠的图像。
    4. 通过以下公式构建吸光度差值的序列图像 ΔAi(x, z)：
      瞬态吸收公式，ΔA(x, τ)=-log10[I₁(x, τ)/I₀(x, τ)]；用于光谱分析。   (1)
      注：ΔAi(x, z) 表示加热前相对于参考吸光度 Ar(x, z) 的吸光度变化量 Ai(x, z)，其推导如下15,16,17,18,19：
      静态平衡方程 ΣFx=0，用于光谱分析的吸收差 ΔA 计算示意图。   (2)
      其中 I0 为入射到样品池的光强。
    5. 使用适当的色彩图（如蓝-红渐变）对 ΔAi 图像进行着色。
      注：用于 ImageJ 软件执行步骤 2.1.2 至 2.1.5 的命令脚本文件见补充材料2。
  2. 温度估算。
    1. 通过目视观察图像，选择 ΔAi(x, z) 关于球体中心呈圆形对称的时间段。
      注：自由对流是破坏圆形对称性的主要原因。先前研究19中已提出基于图像的自由对流发生判断方法；但在实际操作中，目视判断已足够有效。
    2. 在 ΔAi(x, z) 图像上沿 360 条径向线（Δθ = 1°）提取 ΔAi(rʹ, θ) 数据。
    3. 排除球体内部及其附近区域（Δrʹ ≈ 0.2 mm）的 ΔAi(rʹ, θ) 数据。注：由于球体轻微移动，该区域数据异常偏小或偏大。
    4. 对 θ 方向上的 ΔAi(rʹ, θ) 数据进行平均，以确定线型分布 ΔAi(rʹ)。
      注：用于 ImageJ 软件执行步骤 2.2.2 至 2.2.4 的命令脚本文件见补充材料3。
    5. 使用以下多高斯函数对 ΔAi(rʹ) 数据进行拟合：
      静态平衡公式，Σaj[exp(-r²/σ²)-exp(-R²/σ²)]，用于物理方程分析。   (3)
      其中 aj 为权重因子，σj 为弥散参数，R 为可假设 ΔAi(R) = 0 成立时 rʹ 的最大值。
    6. 将拟合得到的 N、aj 和 σj 代入以下反阿贝尔变换（IAT）公式（即公式 (3)）中，计算吸收系数差 Δµi(r)：
      静态分布公式 Δμl(r)，包含求和、指数函数和误差函数。   (4)
      其中 erf 为误差函数。
    7. 通过以下公式将 Δµi(r) 转换为温度：
      热力学方程 ΔTi(r)=Δμi(r)/αf，用于平衡分析的符号公式。   (5) 
      其中水的温度系数 αf 在 λ = 1150 nm 时为 2.8 × 10-4 K-1 mm-119，在 λ = 1412 nm 时为 4.1 × 10-3 K-1 mm-117。
      注：用于执行步骤 2.2.5 至 2.2.7 的命令脚本文件见补充材料4，其中步骤 2.2.5 采用了 Levenberg-Marquardt 非线性最小二乘算法17,19。

访问受限。请登录或开始试用以查看此内容。

结果

水中直径为 2.0 mm 的钢球在 λ = 1150 nm 处以及麦芽糖浆中直径为 0.5 mm 的钢球在 λ = 1412 nm 处的 ΔAi(x, z) 图像分别如图 5a 和图 6a 所示。在这两种情况下，球体均位于线圈底部沿其中心轴线下方 12 mm 处。图 5b 和图 6b 分别展示了 ΔA(rʹ) 数据及其拟合的多高斯函数，其中 R = 3.0 mm 和 R = 1.5 mm。通常仅需两到三个高斯函数（...

访问受限。请登录或开始试用以查看此内容。

讨论

本文介绍的技术是一种利用水的近红外吸收对温度的依赖性的新方法，在所需设备的搭建和实施方面不存在显著困难。入射光可方便地通过卤素灯和窄带通滤光片（NBPF）产生。然而，不能使用激光器，因为相干干涉图样会出现在图像中。常见的用于可见光的光学透镜和玻璃池均可使用，因为它们在 λ = 1150 nm 和 1412 nm 处能透射足够强度的光。此外，目前可相对低成本地购得铟镓砷（InGaAs）相机。

在 λ = 1150 nm 和 1412 nm 处的窄带通滤光片（NBPF）可通过半定制方式订购，但价格并不昂贵。如果已有其他波长的现成窄带通滤光片，只要其波长位于温度依赖性波长范围内（图1），也可替代使用，尽管其温度灵敏度 αf 可能会降低。例如，在 λ = 1175 nm 处的 αf 值仅为 λ = 1150 nm 处的一半。此外，窄带通滤光片的带宽或锐度也会影响 α

访问受限。请登录或开始试用以查看此内容。

披露

作者无任何利益冲突需要披露。

致谢

作者感谢山田健太先生、藤冈亮先生和京田瑞树先生在实验和数据分析方面提供的支持。本研究得到了日本JSPS科研费（编号25630069）、铃木基金会以及精密测量技术促进基金会的资助。

访问受限。请登录或开始试用以查看此内容。

材料

本文使用的材料清单
姓名公司目录编号评论
感应加热系统CEIA，意大利SPW900/56780 kHz，5.6 kW（最大）
线圈SA-日本定制水冷铜管；双匝；外径 28 mm
冷水机松本机械，日本MP-401CT
卤素灯林时计制作所，日本LA-150UE-A
窄带通滤光片，λ = 1150 nmAndover115FS10-25半高全宽（FWHM）：10 nm
窄带通滤光片，λ = 1412 nmAndover半定制半高全宽（FWHM）：10 nm
带通滤光片，λ = 850−1300 nmSpectrogonSP-1300
带通滤光片，λ = 1100−2000 nmSpectrogonSP-2000
近红外相机FLIR SystemsAlpha NIRInGaAs
图像采集软件FLIR SystemsIRvista
图像处理软件NIHImageJ版本 1.51r
图像处理软件MathWorksMatlab版本 2016a
远心镜头Edmond Optics55350-LX1
钢球（0.5 mm 直径）神户制钢，日本Fe-1.5Cr-1.0C-0.4Mn（质量百分比）
钢球（2.0 mm 直径）神户制钢，日本Fe-1.5Cr-1.0C-0.4Mn（质量百分比）
麦芽糖浆（作为水性凝胶）双同，日本Mizuame食品产品

参考文献

  1. Physics of Thermal Therapy. Moros, E. G. , CRC Press. (2012).
  2. Périgo, E. G., Hemery, G., Sandre, O., Ortega, D., Garaio, E., Plazaola, F., Teran, F. J. Fundamentals and advances in magnetic hyperthermia. Appl Phys Rev. 2, 041302(2015).
  3. Bardati, F., Marrocco, G., Tognolatti, P. Time-dependent microwave radiometry for the measurement of temperature in medical applications. IEEE Trans Microwave Theo Tech. 52, 1917-1924 (2004).
  4. Levick, A., Land, D., Hand, J. Validation of microwave radiometry for measuring the internal temperature profile of human tissue. Meas Sci Technol. 22, 065801(2011).
  5. Daniels, M. J., Varghese, T., Madsen, E. L., Zagzebski, J. A. Non-invasive ultrasound-based temperature imaging for monitoring radiofrequency heating-phantom results. Phys Med Biol. 52, 4827(2007).
  6. Daniels, M. J., Varghese, T. Dynamic frame selection for in vivo ultrasound temperature estimation during radiofrequency ablation. Phys Med Biol. 55, 4735(2010).
  7. Seo, C. H., Shi, Y., Huang, S. -W., Kim, K., O'Donnell, M. Thermal strain imaging: A review. Interface Focus. 1, 649-664 (2011).
  8. Bayat, M., Ballard, J. R., Ebbini, E. S. Ultrasound thermography: A new temperature reconstruction model and in vivo results. AIP Conf Proc. 1821, 060004(2017).
  9. Petrova, E., Liopo, A., Nadvoretskiy, V., Ermilov, S. Imaging technique for real-time temperature monitoring during cryotherapy of lesions. J Biomed Opt. 21, 116007(2016).
  10. Gardner, B., Matousek, P., Stone, N. Temperature spatially offset Raman spectroscopy (T-SORS): Subsurface chemically specific measurement of temperature in turbid media using anti-Stokes spatially offset Raman spectroscopy. Anal Chem. 88, 832-837 (2016).
  11. Yoshioka, Y., Oikawa, H., Ehara, S., Inoue, T., Ogawa, A., Kanbara, Y., Kubokawa, M. Noninvasive measurement of temperature and fractional dissociation of imidazole in human lower leg muscles using 1H-nuclear magnetic resonance spectroscopy. J Appl Physiol. 98, 282-287 (2004).
  12. Galiana, G., Branca, R. T., Jenista, E. R., Warren, W. S. Accurate temperature imaging based on intermolecular coherences in magnetic resonance. Science. 322, 421-424 (2008).
  13. Rapoport, E., Pleshivtseva, Y. Optimal Control of Induction Heating Processes. , CRC Press. Boca Raton, FL. (2006).
  14. Lucía, O., Maussion, P., Dede, E. J., Burdío, J. M. Induction heating technology and its applications: Past developments, current technology, and future challenges. IEEE Trans Ind Electron. 61, 2509-2520 (2014).
  15. Kakuta, N., Kondo, K., Ozaki, A., Arimoto, H., Yamada, Y. Temperature imaging of sub-millimeter-thick water using a near-infrared camera. Int J Heat Mass Trans. 52, 4221-4228 (2009).
  16. Kakuta, N., Fukuhara, Y., Kondo, K., Arimoto, H., Yamada, Y. Temperature imaging of water in a microchannel using thermal sensitivity of near-infrared absorption. Lab Chip. 11, 3479-3486 (2011).
  17. Kakuta, N., Kondo, K., Arimoto, H., Yamada, Y. Reconstruction of cross-sectional temperature distributions of water around a thin heating wire by inverse Abel transform of near-infrared absorption images. Int J Heat Mass Trans. 77, 852-859 (2014).
  18. Kakuta, N., Yamashita, H., Kawashima, D., Kondo, K., Arimoto, H., Yamada, Y. Simultaneous imaging of temperature and concentration of ethanol-water mixtures in microchannel using near-infrared dual-wavelength absorption technique. Meas Sci Technol. 27, 115401(2016).
  19. Kakuta, N., Nishijima, K., Kondo, K., Yamada, Y. Near-infrared measurement of water temperature near a 1-mm-diameter magnetic sphere and its heat generation rate under induction heating. J Appl Phys. 122, 044901(2017).
  20. Libnau, F. O., Kvalheim, O. M., Christy, A. A., Toft, J. Spectra of water in the near- and mid-infrared region. Vib Spectrosc. 7, 243-254 (1994).
  21. Siesler, H. W., Ozaki, Y., Kawata, S., Heise, H. M. Near-Infared Spectroscopy. , Wiley-VCH. (2002).
  22. Shakher, C., Nirala, A. K. A review on refractive index and temperature profile measurements using laser-based interferometric techniques. Opt Laser Eng. 31, 455-491 (1999).
  23. Assebana, A., Lallemanda, M., Saulniera, J. -B., Fominb, N., Lavinskaja, E., Merzkirchc, W., Vitkinc, D. Digital speckle photography and speckle tomography in heat transfer studies. Opt Laser Technol. 32, 583-592 (2000).
  24. Ambrosini, D., Paoletti, D., Spagnolo, S. G. Study of free-convective onset on a horizontal wire using speckle pattern interferometry. Int J Heat Mass Trans. 46, 4145-4155 (2003).
  25. Bracewell, R. N. The Fourier Transform and Its Applications. , McGraw-Hill. (2000).
  26. Yoder, L. M., Barker, J. R., Lorenz, K. T., Chandler, D. W. Ion imaging the recoil energy distribution following vibrational predissociation of triplet state pyrazine-Ar van der Waals clusters. Chem Phys Lett. 302, 602-608 (1999).
  27. De Colle, F., de Burgo, C., Raga, A. C. Diagnostics of inhomogeneous stellar jets: convolution effects and data reconstruction. Astron Astrophys. 485, 765-772 (2008).
  28. Green, K. M., Borrás, M. C., Woskov, P. P., Flores, G. J., Hadidi, K., Thomas, P. Electronic excitation temperature profiles in an air microwave plasma torch. IEEE Trans Plasma Sci. 29, 399-406 (2001).
  29. Bendinelli, O. Abel integral equation inversion and deconvolution by multi-Gaussian approximation. Astrophys J. 366, 599-604 (1991).
  30. Dorband, B., Muller, H., Gross, H. Vol. 5 Metrology of Optical Components and Systems. Handbook of Optical System. Gross, H. , Wiley-VCH. (2012).
  31. Sheikholeslami, M., Rokni, H. B. Simulation of nanofluid heat transfer in presence of magnetic field: A review. Int J Heat Mass Trans. 115, 1203-1233 (2017).
  32. Häfeli, U., Schütt, W., Teller, J., Zborowski, M. Scientific and Clinical Applications of Magnetic Carriers. , Springer Science and Business Media. (2013).

访问受限。请登录或开始试用以查看此内容。

重印与许可

勘误


Formal Correction: Erratum: Near-Infrared Temperature Measurement Technique for Water Surrounding an Induction-heated Small Magnetic Sphere
Posted by JoVE Editors on 12/06/2018. Citeable Link.

An erratum was issued for: Near-Infrared Temperature Measurement Technique for Water Surrounding an Induction-heated Small Magnetic Sphere. The Protocol section was updated.

In 2.2.7, the temperature coefficient of water, αf, for λ = 1150 nm has been corrected from:

4.0 x 10-3 K-1 mm-1

to:

2.8 x 10-4 K-1 mm-1

探索更多文章

感应加热水吸收率窄带通滤光片近红外相机阿贝尔变换磁性微球加热热传导自然对流光程长度
近红外测温技术用于测量感应加热小型磁性球体周围水的温度
联系我们
向图书馆推荐

隐私政策

使用条款

政策

版权所有 © 2026 MyJoVE Corporation。保留所有权利。