CN112985613B - 一种THz回波高温温度测量方法 - Google Patents
一种THz回波高温温度测量方法 Download PDFInfo
- Publication number
- CN112985613B CN112985613B CN202110176882.7A CN202110176882A CN112985613B CN 112985613 B CN112985613 B CN 112985613B CN 202110176882 A CN202110176882 A CN 202110176882A CN 112985613 B CN112985613 B CN 112985613B
- Authority
- CN
- China
- Prior art keywords
- thz
- echo
- temperature
- calibration
- reference signal
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 33
- 238000009529 body temperature measurement Methods 0.000 claims abstract description 80
- 230000003287 optical effect Effects 0.000 claims abstract description 8
- 238000001228 spectrum Methods 0.000 claims description 109
- 238000005070 sampling Methods 0.000 claims description 14
- 238000004422 calculation algorithm Methods 0.000 claims description 11
- 230000009466 transformation Effects 0.000 claims description 7
- 238000004616 Pyrometry Methods 0.000 claims description 2
- 238000004364 calculation method Methods 0.000 claims description 2
- 230000001131 transforming effect Effects 0.000 claims description 2
- 238000001514 detection method Methods 0.000 abstract description 10
- 238000005516 engineering process Methods 0.000 description 15
- 238000005259 measurement Methods 0.000 description 15
- XEEYBQQBJWHFJM-UHFFFAOYSA-N Iron Chemical compound [Fe] XEEYBQQBJWHFJM-UHFFFAOYSA-N 0.000 description 12
- XOFYZVNMUHMLCC-ZPOLXVRWSA-N prednisone Chemical compound O=C1C=C[C@]2(C)[C@H]3C(=O)C[C@](C)([C@@](CC4)(O)C(=O)CO)[C@@H]4[C@@H]3CCC2=C1 XOFYZVNMUHMLCC-ZPOLXVRWSA-N 0.000 description 12
- 239000002184 metal Substances 0.000 description 10
- 229910052751 metal Inorganic materials 0.000 description 10
- 230000005855 radiation Effects 0.000 description 9
- 238000010586 diagram Methods 0.000 description 7
- 239000000428 dust Substances 0.000 description 7
- 229910052742 iron Inorganic materials 0.000 description 6
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 description 5
- 230000008859 change Effects 0.000 description 5
- 229910002804 graphite Inorganic materials 0.000 description 5
- 239000010439 graphite Substances 0.000 description 5
- 230000007704 transition Effects 0.000 description 4
- 229910001218 Gallium arsenide Inorganic materials 0.000 description 3
- 230000007547 defect Effects 0.000 description 3
- 239000010410 layer Substances 0.000 description 3
- 229910052782 aluminium Inorganic materials 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 230000005670 electromagnetic radiation Effects 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 230000004927 fusion Effects 0.000 description 2
- 229910052733 gallium Inorganic materials 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 230000000149 penetrating effect Effects 0.000 description 2
- 230000035515 penetration Effects 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 230000004044 response Effects 0.000 description 2
- 239000000779 smoke Substances 0.000 description 2
- 238000012549 training Methods 0.000 description 2
- 229910000980 Aluminium gallium arsenide Inorganic materials 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000036760 body temperature Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 230000009970 fire resistant effect Effects 0.000 description 1
- 239000012535 impurity Substances 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 238000012821 model calculation Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000011241 protective layer Substances 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
- 238000004861 thermometry Methods 0.000 description 1
- 230000001052 transient effect Effects 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Chemical compound O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
- G01J5/00—Radiation pyrometry, e.g. infrared or optical thermometry
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S17/00—Systems using the reflection or reradiation of electromagnetic waves other than radio waves, e.g. lidar systems
- G01S17/02—Systems using the reflection of electromagnetic waves other than radio waves
- G01S17/06—Systems determining position data of a target
- G01S17/08—Systems determining position data of a target for measuring distance only
Landscapes
- Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开了一种THz回波高温温度测量方法，通过生成指定THz频率的THz源信号与THz参考信号，其中THz源信号用于通过THz波采集光学镜头传输并聚焦到被测高温对象上，检测THz源信号经待测高温对象反射的THz回波，测量被测高温对象与THz回波高温温度测量装置之间的距离以及根据THz参考信号、THz回波以及距离获得被测高温对象的温度值，解决了现有非接触式测温装置在复杂恶劣环境下对高温物体温度检测精度低的技术问题，巧妙利用太赫兹波穿透性强的特点，在保证复杂环境中还能接受到太赫兹信号，从而对被测高温源进行精确测温，同时具有非接触的特点，这对于高温物体复杂环境内的测温具有重要意义。
Description
技术领域
本发明主要涉及测温技术领域，特指一种THz(TeraHertz太赫兹)回波高温温度测量方法。
背景技术
太赫兹(又称THz波、T射线)通常是指频率在0.1-10THz范围内的电磁辐射，在电磁波谱上位于微波和红外线之间。THz频段是一个非常有科学价值但尚未被完全认识和利用的最后一个电磁辐射区域。许多年来，由于缺乏切实可行的THz波产生方法和检测手段，人们对THz波段的特性知之甚少，以致于该波段被称为电磁波谱中的“THz空隙”，太赫兹波的频率范围处于电子学与光子学的交叉蔓域。在长波方向，它与毫米波有重叠；在短波方向，它与红外线有重叠。在频域上，太赫兹处于宏观经典理论向微观量子理论的过渡区。这就导致太赫兹波即具有类似微波的强穿透性又具有红外线可测温的特性。而THz波这一奇特的特性，使基于THz波的测温技术是实现强干扰环境下、特别是强粉尘的工业环境中，非接触测量高温物体温度的绝佳测温技术。尤其是近年来随着THz量子级联激光器技术(QuantumCascade Lasers,QCL)和基于GaAs/(Al,Ga)As量子阱子带跃迁效应的THz量子阱探测器(Quantum Well Photodetectors，QWP)技术的高速发展和成熟，使得利用THz波的相关特性，来设计具有实用价值的高端前沿工业传感器成为可能。
在高温对象非接触式测温领域，基于近红外、红外的测温方法和技术，由于其传感器件和理论方法均比较成熟，一直是该领域的首选技术方案。但该技术缺陷也十分明显，近红外线和红外线的波长均较短，光子频率高能量高，穿透能力差，极易被环境中杂物阻挡，且由于任意物体均向外辐射红外线，这导致当高温对象处于强粉尘、强干扰等复杂恶劣环境中时，其测量测温精度低、波动大实用性差。尤其是在环境极端恶劣，而对测温技术又极端依赖的工业领域，红外测温技术的缺陷更为明显，因此，研发基于THz波的不受环境影响的非接触式测温技术和方法极为迫切。
专利公布号CN107941357发明专利是一种高温金属熔液温度测量方法及其装置。其通过将热电偶包裹在耐温体中，在测量时将内含热电偶的耐火体浸渍到待测金属熔液中，通过给热电偶套上保护层，隔绝热电偶与避高温金属熔液直接接触，短时起到保护热电偶作用，同时高温金属熔液的热量经由耐火体的传导至热电偶，实现对高温金属熔液进行接触测温。该方法是当前工业界测量高温金属熔液最常采用的方案，但该测温技术路线的缺陷不足也十分明显。首先，对高温金属熔液进行接触测量，必极易损坏耐火体和热电偶，因此这类测温传感器通常是一次性的，俗称“快偶”测温，即每测量一次均需更换，使用成本高；其次，使用快偶测温需要人工操作，而在高温金属熔液旁边进行测温操作是比较危险的工作，极易出现人员安全事故；最后，由于是人工进行测温操作，因此无法实时连续的获取高温金属熔液的温度，难以全面反应高温金属熔液的变化趋势。
专利公布号CN101545808A发明专利是一种高温铁水的红外辐射测温***。该专利采用了一种隶属于红外测温的比色红外辐射测温技术，其主要工作原理是通过测量被测物红外辐射中相邻波段的能量大小来确定被测物的温度，因此其在烟雾、水汽等环境中测量精度仍有保证。但红外测温的致命问题仍然存在，即当空间布满粉尘时，两个对比红外光线被阻挡，而造成接受的红外能量降低，导致测量温度极低，完全偏离了被测对象真实温度。因此，该装置不适用于测量粉尘浓度大的出铁口处铁水温度，在实际使用中，常将其安装在粉尘量小的铁沟槽的位置，用于测量铁水进溜槽时的温度，其实际使用受到了严重制约。
专利公布号CN103033282A发明专利是当前少有的利用THz波测温高温对象的温度的专利，该专利的测量对象是磁约束聚变装置偏滤器石墨瓦瞬态温度，该专利利用太赫兹波射入石墨反射回的太赫兹谱的相位和频谱会随着温度的变化而发生变化和频移这一特性，通过理论计算不同温度与标定温度(如常温300K)特征谱线位置的频移距离，从而确定出温度随频移量的变化函数关系，作为测温曲线,集成为计算机数据库，将石墨瓦工作温度与标定温度太赫兹频域谱的频移距离和数据库相对比，以此可计算出石墨瓦所处环境的温度，从而实现对石墨瓦的温度进行实时测量。该专利的思路较为新颖，其技术路线从表面上看与本发明专利有些许类似，但两者却属于完全不同技术思路和途经，不可相互借鉴。具体而言：①从THz波的理论研究可知，要使THz波的特征频率出现可检测的偏移，需要极高的温度(通常几万℃以上)，在核聚变装置中存在这样的高温，但通常的民用及工业领域几乎不存在如此高温，换句话说在通常情况下，无法采用该原理进行测温；②从目前的THz源产生技术和检测技术而言，难以产生THz波的功率通常在(pw～mw)级别，因此，检测的THz波信号时域谱完全取淹没在环境噪声中，信噪比极低，而要在如此低信噪比的信号中，提取THz波信号的特征频率极为困难；③即使提取了THz波信号的特征频率，但由于大量噪声的存在与THz波信号时域谱中，当对其进行快速傅立叶变化后，其频域谱还存在大量重叠峰的问题，也难于分离和区分THz特征谱的特征峰所在位置，测量精度也较低。因此该专利所提的基于THz波的高温测量技术，短期还难应用于真正工业实际。
发明内容
本发明提供的THz回波高温温度测量方法，解决了现有非接触式测温装置在复杂恶劣环境下对高温物体温度检测精度低的技术问题。
为解决上述技术问题，本发明提出的THz回波高温温度测量方法包括：
生成指定THz频率的THz源信号与THz参考信号，其中THz源信号用于通过THz波采集光学镜头传输并聚焦到被测高温对象上；
检测THz源信号经待测高温对象反射的THz回波；
测量被测高温对象与THz回波高温温度测量装置之间的距离；
根据THz参考信号、THz回波以及距离获得被测高温对象的温度值。
进一步地，根据THz参考信号、THz回波以及距离获得被测高温对象的温度值包括：
根据THz参考信号获得THz参考信号时域谱，根据THz回波获得THz回波信号时域谱；
对THz参考信号时域谱和THz回波信号时域谱进行线性快速傅立叶变换，获得THz参考信号频域谱和THz回波信号频域谱；
根据THz参考信号频域谱、THz回波信号频域谱及距离获得被测高温对象的温度值。
进一步地，对THz回波信号时域谱进行线性快速傅立叶变换，获得THz回波信号频域谱包括：
利用传统快速傅立叶变换算法，对THz回波信号时域谱进行变换；
对采用传统快速傅立叶变换后的THz回波信号时域谱，进行线性快速傅立叶变换，获得THz回波信号频域谱，其中线性快速傅立叶变换的计算公式为：
其中，
代表k次谐波的复振幅，
为传统的FFT算法傅立叶变换后的THz回波信号时域谱结果，x0,xn为信号离散的采样点，Zk为变换辅助复数，uk为辅助复数的实部，vk辅助复数的虚部，为N为采样点数，F为采样步长，k代表第k个采样变换点。
进一步地，根据THz参考信号频域谱、THz回波信号频域谱及距离获得被测高温对象的温度值包括：
根据THz参考信号频域谱，获得THz参考信号频域幅值；
根据THz回波信号频域谱，获得THz回波信号频域幅值；
计算指定工作THz频率下THz参考信号频域幅值和THz回波信号频域幅值之间的偏差，获得THz参考信号与THz回波之间的幅值偏差量；
根据不同的温度对应的幅值偏差量和距离不同，拟合构建THz回波测温模型；
基于THz回波测温模型，获得被测高温对象的温度值。
进一步地，根据不同的温度对应的幅值偏差量和距离不同，拟合构建THz回波测温模型包括：
实时获得THz回波高温温度测量装置与黑体之间的标定距离，黑体用于产生所需的可控高温；
采集黑体在标定工作温度下的THz标定参考信号时域谱以及THz标定回波信号时域谱，同步记录下黑体上显示的标定温度；
对THz标定参考信号时域谱和THz标定回波信号时域谱进行线性快速傅立叶变换，获得THz标定参考信号频域谱和THz标定回波信号频域谱；
根据THz标定参考信号频域谱，获得THz标定参考信号频域幅值；
根据THz标定回波信号频域谱，获得THz标定回波信号频域幅值；
计算THz标定参考信号频域幅值和THz标定回波信号频域幅值之间的偏差，获得THz标定参考信号与THz标定回波之间的标定幅值偏差量；
根据标定温度与标定幅值偏差量和标定距离之间的关系，拟合构建THz回波测温模型。
进一步地，根据THz回波获得THz回波信号时域谱包括：
采用前置源放大器和锁相放大器提取THz回波的THz回波信号时域谱。
进一步地，生成指定THz频率的THz源信号与THz参考信号之前包括：
确定THz波测温频谱段，THz波测温频谱段为1.3THz。
进一步地，生成指定THz频率的THz源信号与THz参考信号具体为通过THz量子级联激光器生成。
进一步地，检测THz源信号经待测高温对象反射的THz回波具体为通过THz量子阱探测器检测。
进一步地，测量被测高温对象与THz回波高温温度测量装置之间的距离具体为通过测距激光雷达测量。
与现有技术相比，本发明的优点在于：本发明提供的THz回波高温温度测量方法，通过生成指定THz频率的THz源信号与THz参考信号，其中THz源信号用于通过THz波采集光学镜头传输并聚焦到被测高温对象上，检测THz源信号经待测高温对象反射的THz回波，测量被测高温对象与THz回波高温温度测量装置之间的距离以及根据THz参考信号、THz回波以及距离获得被测高温对象的温度值，解决了现有非接触式测温装置在复杂恶劣环境下对高温物体温度检测精度低的技术问题，巧妙利用太赫兹波穿透性强的特点，在保证复杂环境中还能接受到太赫兹信号，从而对被测高温源进行精确测温，同时具有非接触的特点，这对于高温物体复杂环境内的测温具有重要意义。
此外，本发明通过先主动发射功率可控可调的待测高温对象测温敏感的THz频段波，照射到高温对象上，由于高温对象本身也会辐射该频段的THz波，且该频段THz波的辐射强度随温度正相关变化，导致从高温对象上反射回来的THz回波能量叠加了高温对象在该频段自身辐射THz波的能量，使THz回波能量达到了足够的能量裕度，能被QWP探测器探测到。同时为了提高THz回波信号反应被测对象温度信息的信噪比，在探测到指定工作THz频率Ts下THz回波信号的幅值(能量强度)后，将其与参考THz波信号的幅值进行做差，提取该频段发射的参考THz波和测量的THz回波幅值偏差量ΔE，作为反应被测对象温度高低的关键量，从而达到提高THz回波信号反应被测对象温度信息的信噪比的目的，进而提高测温精度。
本发明的目的在于提出了一种全新的基于THz量子级联激光器技术(QCL)和THz量子阱探测器(QWP)技术的THz工作频率下的回波信号幅值偏移量提取技术，用于稳定、高效高精度提取回波信号幅值偏移量。
本发明的目的在于构建了基于测温距离及幅值偏差量的THz回波测温模型，实现对高温对象的温度进行实时测量。
本发明的目的在于为高温对象测温提供了一种基于THz回波的一种全新的主动式非接触式测温方法，利用THz波的强穿透性，实现了高粉尘、强干扰等恶劣环境下实时高精度测温，解决了恶劣环境中非接触式测温精度不高、波动大、实用性差的难题。
附图说明
图1为本发明实施例二的THz回波高温温度测量方法的流程图；
图2为本发明实施例二的THz参考信号时域谱和THz回波信号时域谱的示意图；
图3为本发明实施例二的THz参考与回波信号在指定频率下的幅值偏差量示意图；
图4为本发明实施例二的指定测温距离及工作频率下幅值偏差量－温度函数关系模型示意图。
附图标记：
U0：被测高温对象；U11：THz参考信号时域谱功能模块；U12：THz回波信号时域谱功能模块；U13：测距激光雷达测距信息功能模块；U21：线性快速傅立叶变换功能模块；U31：参考信号频域幅值谱功能模块；U32：THz回波信号频域幅值谱功能模块；U33：被测高温对象与测温设备间距离功能模块；U41：指定工作THz频率下幅值偏差量ΔE功能模块；U51：温度可控高温黑体；U52：依照标定步长设定黑体温度作为标定温度功能模块；U53：黑体达到某一标定温度并稳定后作为实验工作温度功能模块；U61：测温距离功能模块；U62：标定温度功能模块；U63：标定温度参考THz信号时域谱功能模块；U64：标定温度THZ回波信号时域谱功能模块；U65：快速傅立叶变换并计算工作频率下频域谱幅值偏差ΔE功能模块；U71：基于测温距离及幅值偏差量构建THZ回波测温模型功能模块；U81：实现对高温对象温度的实时测量功能模块。
具体实施方式
为了便于理解本发明，下文将结合说明书附图和较佳的实施例对本发明作更全面、细致地描述，但本发明的保护范围并不限于以下具体的实施例。
以下结合附图对本发明的实施例进行详细说明，但是本发明可以由权利要求限定和覆盖的多种不同方式实施。
实施例一
本发明实施例一提供的THz回波高温温度测量方法包括：
步骤S101，生成指定THz频率的THz源信号与THz参考信号，其中THz源信号用于通过THz波采集光学镜头传输并聚焦到被测高温对象上；
步骤S102，检测THz源信号经待测高温对象反射的THz回波；
步骤S103，测量被测高温对象与THz回波高温温度测量装置之间的距离；
步骤S104，根据THz参考信号、THz回波以及距离获得被测高温对象的温度值。
本发明提供的THz回波高温温度测量方法，包括生成指定THz频率的THz源信号与THz参考信号，其中THz源信号用于通过THz波采集光学镜头传输并聚焦到被测高温对象上，检测THz源信号经待测高温对象反射的THz回波，测量被测高温对象与THz回波高温温度测量装置之间的距离以及根据THz参考信号、THz回波以及距离获得被测高温对象的温度值，解决了现有非接触式测温装置在复杂恶劣环境下对高温物体温度检测精度低的技术问题，巧妙利用太赫兹波穿透性强的特点，在保证复杂环境中还能接受到太赫兹信号，从而对被测高温源进行精确测温，同时具有非接触的特点，这对于高温物体复杂环境内的测温具有重要意义。
实施例二
如图1所示，本发明实施例二提供的THz回波高温温度测量方法包括：
S0：THz波测温频谱段确定。即依据被测高温对象所处的环境，确定待测高温对象U0测温敏感的THz波频谱段Ts，本实施例设定Ts＝1.3THz。
S1：THz参考信号时域谱、THz回波信号时域谱以及被测高温对象与THz回波高温温度测量装置之间的距离信息采集。该功能需三个功能模块参与，分别为THz量子级联激光器(QCL)产生的THz参考信号时域谱功能模块U11、THz量子阱探测器(QWP)探测到的THz回波信号时域谱功能模块U12和激光测距仪测量到被测高温对象与测温设备间距离的测距激光雷达测距信息功能模块U13。其中量子级联激光器(QCL)是一种采用GaAs/AlGaAs材料体系，基于多量子阱子带间跃迁的单极性半导体激光器，是最有效的电泵浦THz辐射源，具有结构紧凑、易集成、输出功率高和转换效率高等优点。而THz量子阱探测器(QWP)是一种基于GaAs/(Al，Ga)As材料体系工作于太赫兹频段的光子型探测器，其由上电极层、多量子阱层和下电极层构成，利用束缚电子吸收太赫兹辐射发生子带跃迁，来实现对太赫兹辐射的探测和响应的。而为了实现THz波时域谱及测量距离信息的采集的目的，上述三个功能模块的工作流程为：首先，由THz量子级联激光器(QCL)源产生频率段为Ts的THz波，并将其分成全同的两部分，一部分作为THz参考信号时域谱，另一部分通过THz传输光路聚焦到被测高温对象上；其次，高温对象反射回来的THz回波，经过THz量子阱探测器(QWP)的感知，再经过前置源放大器和锁相放大器提取THz回波信号时域谱；最后，再由测距激光雷达测量被测高温对象与测温设备间距离d，完成THz波时域谱及测量距离信息的采集工作。
S2：THz参考信号频域谱和THz回波信号频域谱的获取。对获取的THz参考信号及回波信号的时域谱，如图2中所示，观察图2可知，THz参考及回波信号的时域谱均存在大量的干扰和噪声。因此为了能高精度的提取THz波频域谱，本实施例针对THz时域谱信号高噪声的特点，将传统快速傅立叶变换对变换信号简单离散采样，修改成基于线性插值原理的数据采样方法，从而提高快速傅立叶的精度，修改后的线性快速傅立叶变换(LFFT)公式如下式(1)所示：
对于uk,vk的计算可由下式式(2)表达的形式来进行计算，
其中，
代表k次谐波的复振幅，
为传统的FFT算法傅立叶变换后的THz回波信号时域谱结果，x0,xn为信号离散的采样点，Zk为变换辅助复数，uk为辅助复数的实部，vk辅助复数的虚部，为N为采样点数，F为采样步长，k代表第k个采样变换点。而采用线性快速傅立叶变换(LFFT)算法提取THz波频域谱的流程为：先利用传统快速傅立叶变换(FFT)算法，对THz波时域谱进行变换，再将变换结果代入(LFFT)算法的式(1)、(2)中，对k次谐波的复振幅
进行修正，从而达到提高变换精度的目的。综上，在S2中通过采用线性快速傅立叶变换功能模块U21中的线性快速傅立叶变换(LFFT)算法，实现对THz参考信号及回波信号的时域谱的高速、高精度频率变换，获得其频域谱。
S3：THz波信号频域幅值谱提取。本实施例通过参考信号频域幅值谱功能模块U31和THz回波信号频域幅值谱功能模块U32，对获取到的THz参考信号频域谱及回波信号频域谱，取出高温对象敏感的THz波频谱段Ts的幅值，从而依次获得THz参考信号频域幅值，THz回波信号频域幅值，以及通过被测高温对象与测温设备间距离功能模块U32获得由激光测距雷达测量的被测高温对象与THz回波高温温度测量装置之间的距离。
S4：获得指定工作THz频率Ts下THz参考信号与回波信号间幅值偏差量ΔE。该步骤通过指定工作THz频率下幅值偏差量ΔE功能模块U41实现，是整个方案技术的关键点，经过当前最新关于THz波性质的理论研究可知，高温对象均会自发性对外辐射出全谱电磁波，而且对象温度越高，辐射电磁波能量越大。但辐射的电磁波强度在相同温度下，不同的频段强度是不同的，尤其是在THz频段，能量尤其低。因此，要利用THz波的强穿透性，来克服环境中粉尘、烟雾、噪音等干扰，实现对高温对象进行测温，就必须提高能反应被测对象温度高低的THz波信号能量强度和信噪比，使其达到能被QWP探测的程度。因此，本实施例先主动发射功率可控可调的待测高温对象测温敏感的THz频段波，照射到高温对象上，由于高温对象本身也会辐射该频段的THz波，且该频段THz波的辐射强度随温度正相关变化，导致从高温对象上反射回来的THz回波能量叠加了高温对象在该频段自身辐射THz波的能量，使THz回波能量达到了足够的能量裕度，能被QWP探测器探测到。同时为了提高THz回波信号反应被测对象温度信息的信噪比，在探测到指定工作THz频率Ts下THz回波信号的幅值(能量强度)后，将其与参考THz波信号的幅值进行做差，提取该频段发射的参考THz波和测量的THz回波幅值偏差量ΔE，作为反应被测对象温度高低的关键量，从而达到提高THz回波信号反应被测对象温度信息的信噪比的目的。上述原理的示意图如图2所示，图2中展示的为某一黑体标定工作温度下实测的THz参考信号时域谱和THz回波信号时域谱，比较两图可发现THz回波信号的振幅最大值为1.5(a.u.)大于THz参考信号的振幅最大值0.8(a.u.)。而在图3所示的THz参考信号及回波信号频域谱幅值图，显然在指定工作频率为Ts＝1.3THz处，参考THz波和测量的THz回波幅值偏差量为ΔE＝0.60(a.u.)，因此，从实验的实测数据来看，进一步验证了，采用本实施例的技术路线，确实能提高THz回波信号反应被测对象温的能量强度和信噪比。
S5：实验室选定黑体标定温度，并利用黑体产生所需的实验工作温度。该功能，首先需要一个温度可控高温黑体U51，用于产生所需的可控高温。其次，还需根据被测高温对象实际温度所处范围，确定标定步长，并依据标定步长，设定标定温度，完成依照标定步长设定黑体温度作为标定温度功能模块U52的功能。再次，通过控制让黑体达到某一标定温度并稳定后，即可作为实验工作温度，实现黑体达到某一标定温度并稳定后作为实验工作温度功能模块U53的功能。
S6：获取构建THz回波测温模型所需的实验室数据。具体包括：首先，在实验室通过测温距离功能模块U61利用激光测距雷达实时获得测温设备与黑体之间的距离信息；其次，通过标定温度参考THz信号时域谱功能模块U63和标定温度THZ回波信号时域谱功能模块U64，利用步骤S1中所提的基于THz量子级联激光器(QCL)和THz量子阱探测器(QWP)构建的THz波时域信号采集器，采集黑体在标定工作温度下的标定温度参考THz信号时域谱以及THz回波信号时域谱，同步在标定温度功能模块U62中记录下黑体上显示的标定温度；最后，利用步骤S2中所述的线性快速傅立叶变换(LFFT)算法，实现快速傅立叶变换并计算工作频率下频率谱幅值偏差量ΔE，完成快速傅立叶变换并计算工作频率下频域谱幅值偏差ΔE功能模块U65的功能。经过该步骤为建立基于测温距离及幅值偏差量构建THz回波测温模型，奠定了扎实的实验数据支撑。
S7：将步骤S6中获取的测温距离数据、黑体标定温度数据和工作频率下频率谱幅值偏差量数据，按照时间精准匹配后，构建模型拟合所需的标准训练数据样本集(Tk,ΔEk,dk),k＝1,2,…,l，其中k代表实验中获得的第k个标定温度下的实验数据，l代表标定温度的总个数。则通过基于测温距离及幅值偏差量构建THZ回波测温模型功能模块U71，可获得基于测温距离和幅值偏差量的THz回波测温模型可表示为如式(3)形式，
其中p,q表示拟合多项式的次数，通常人为给定，而对于系数aij则可再利用最小二乘算法，结合标准训练数据样本集(Tk,ΔEk,dk)对式(3)进行拟合得到，从而完成基于测温距离和幅值偏差量的THz回波测温模型T＝F(ΔE,d)的构建工作。
S8：实现对高温对象温度的实时测量。通过实现对高温对象温度的实时测量功能模块U81，将步骤S4获取的被测高温对象与THz回波高温温度测量装置之间的距离信息d和指定工作THz频率下幅值偏差量ΔE，输入步骤S7中构建的基于测温距离和幅值偏差量的THz回波测温模型T＝F(ΔE,d)中，即可通过模型计算获得被测高温对象的温度。在图4中，指定测温距离及工作频率下幅值偏差量－温度函数关系模型示意图展示出了步骤S8的这一个过程，当被测高温对象与测温设备间距离信息d确定后，被测对象的温度与幅值偏差量ΔE间的函数关系图即如图4中曲线所示，只需将幅值偏差量ΔE代入到图示的函数关系曲线中，即可获得当前高温被测对象的工作温度T，实现实时测温。
本发明提供的THz回波高温温度测量方法，通过生成指定THz频率的THz源信号与THz参考信号，其中THz源信号用于通过THz波采集光学镜头传输并聚焦到被测高温对象上，检测THz源信号经待测高温对象反射的THz回波，测量被测高温对象与THz回波高温温度测量装置之间的距离以及根据THz参考信号、THz回波以及距离获得被测高温对象的温度值，解决了现有非接触式测温装置在复杂恶劣环境下对高温物体温度检测精度低的技术问题，巧妙利用太赫兹波穿透性强的特点，在保证复杂环境中还能接受到太赫兹信号，从而对被测高温源进行精确测温，同时具有非接触的特点，这对于高温物体复杂环境内的测温具有重要意义。
此外，本实施例通过采集的被测高温源发射的太赫兹实现对被测高温源的温度检测，能借助太赫兹的穿透特性减少环境中遮挡物对测温结果造成的误差，进一步提高被测高温源的温度检测精度。且本实施例通过先主动发射功率可控可调的待测高温对象测温敏感的THz频段波，照射到高温对象上，由于高温对象本身也会辐射该频段的THz波，且该频段THz波的辐射强度随温度正相关变化，导致从高温对象上反射回来的THz回波能量叠加了高温对象在该频段自身辐射THz波的能量，使THz回波能量达到了足够的能量裕度，能被QWP探测器探测到。同时为了提高THz回波信号反应被测对象温度信息的信噪比，在探测到指定工作THz频率Ts下THz回波信号的幅值(能量强度)后，将其与参考THz波信号的幅值进行做差，提取该频段发射的参考THz波和测量的THz回波幅值偏差量ΔE，作为反应被测对象温度高低的关键量，从而达到提高THz回波信号反应被测对象温度信息的信噪比的目的，进而提高测温精度。
以上实施方式仅用于说明本发明，而非对本发明的限制。尽管参照实施例对本发明进行了详细说明，本领域的普通技术人员应当理解，对本发明的技术方案进行各种组合、修改或者等同替换，都不脱离本发明技术方案的精神和范围，均应涵盖在本发明的权利要求范围当中。
Claims (7)
1.一种THz回波高温温度测量方法，其特征在于，所述方法包括：
生成指定THz频率的THz源信号与THz参考信号，其中所述THz源信号用于通过THz波采集光学镜头传输并聚焦到被测高温对象上；
检测所述THz源信号经待测高温对象反射的THz回波；
测量被测高温对象与THz回波高温温度测量装置之间的距离；
根据所述THz参考信号获得THz参考信号时域谱，根据所述THz回波获得THz回波信号时域谱，对所述THz参考信号时域谱和THz回波信号时域谱进行线性快速傅立叶变换，获得THz参考信号频域谱和THz回波信号频域谱，根据所述THz参考信号频域谱，获得THz参考信号频域幅值，根据所述THz回波信号频域谱，获得THz回波信号频域幅值，计算指定工作THz频率下所述THz参考信号频域幅值和THz回波信号频域幅值之间的偏差，获得THz参考信号与THz回波之间的幅值偏差量；
根据不同的温度对应的幅值偏差量和距离不同，拟合构建THz回波测温模型，其中根据不同的温度对应的幅值偏差量和距离不同，拟合构建THz回波测温模型包括实时获得THz回波高温温度测量装置与黑体之间的标定距离，所述黑体用于产生所需的可控高温，采集黑体在标定工作温度下的THz标定参考信号时域谱以及THz标定回波信号时域谱，同步记录下黑体上显示的标定温度，对所述THz标定参考信号时域谱和THz标定回波信号时域谱进行线性快速傅立叶变换，获得THz标定参考信号频域谱和THz标定回波信号频域谱，根据所述THz标定参考信号频域谱，获得THz标定参考信号频域幅值，根据所述THz标定回波信号频域谱，获得THz标定回波信号频域幅值，计算所述THz标定参考信号频域幅值和THz标定回波信号频域幅值之间的偏差，获得THz标定参考信号与THz标定回波之间的标定幅值偏差量，根据所述标定温度与所述标定幅值偏差量和标定距离之间的关系，拟合构建THz回波测温模型；
将被测高温对象与THz回波高温温度测量装置之间的距离信息和指定工作THz频率下的THz参考信号与THz回波之间的幅值偏差量，输入所述THz回波测温模型，从而获得被测高温对象的温度。
2.根据权利要求1所述的THz回波高温温度测量方法，其特征在于，对所述THz回波信号时域谱进行线性快速傅立叶变换，获得THz回波信号频域谱包括：
利用传统快速傅立叶变换算法，对所述THz回波信号时域谱进行变换；
对采用传统快速傅立叶变换后的THz回波信号时域谱，进行线性快速傅立叶变换，获得THz回波信号频域谱，其中线性快速傅立叶变换的计算公式为：
其中，
代表k次谐波的复振幅，
为传统的FFT算法傅立叶变换后的THz回波信号时域谱结果，x0,xn为信号离散的采样点，Zk为变换辅助复数，uk为辅助复数的实部，vk辅助复数的虚部，为N为采样点数，F为采样步长，k代表第k个采样变换点。
3.根据权利要求2所述的THz回波高温温度测量方法，其特征在于，根据所述THz回波获得THz回波信号时域谱包括：
采用前置源放大器和锁相放大器提取所述THz回波的THz回波信号时域谱。
4.根据权利要求3所述的THz回波高温温度测量方法，其特征在于，生成指定THz频率的THz源信号与THz参考信号之前包括：
确定THz波测温频谱段，所述THz波测温频谱段为1.3THz。
5.根据权利要求4所述的THz回波高温温度测量方法，其特征在于，生成指定THz频率的THz源信号与THz参考信号具体为通过THz量子级联激光器生成。
6.根据权利要求5所述的THz回波高温温度测量方法，其特征在于，检测所述THz源信号经待测高温对象反射的THz回波具体为通过THz量子阱探测器检测。
7.根据权利要求6所述的THz回波高温温度测量方法，其特征在于，测量被测高温对象与THz回波高温温度测量装置之间的距离具体为通过测距激光雷达测量。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110176882.7A CN112985613B (zh) | 2021-02-07 | 2021-02-07 | 一种THz回波高温温度测量方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110176882.7A CN112985613B (zh) | 2021-02-07 | 2021-02-07 | 一种THz回波高温温度测量方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112985613A CN112985613A (zh) | 2021-06-18 |
| CN112985613B true CN112985613B (zh) | 2021-12-28 |
Family
ID=76392543
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110176882.7A Active CN112985613B (zh) | 2021-02-07 | 2021-02-07 | 一种THz回波高温温度测量方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112985613B (zh) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114355463B (zh) * | 2021-12-30 | 2024-04-26 | 中国科学院沈阳自动化研究所 | 一种使用太赫兹波进行声音探测与感知的方法及*** |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103033282A (zh) * | 2012-10-29 | 2013-04-10 | 大连理工大学 | 一种检测磁约束聚变装置偏滤器石墨瓦瞬态温度的方法 |
| CN103048061A (zh) * | 2012-10-29 | 2013-04-17 | 大连理工大学 | 反射太赫兹谱技术检测偏滤器石墨瓦瞬态温度的装置 |
| CN103115893A (zh) * | 2013-01-30 | 2013-05-22 | 大连理工大学 | 一种检测托卡马克钨第一壁灰尘沉积层成分及厚度的装置 |
| JP2016142557A (ja) * | 2015-01-30 | 2016-08-08 | いすゞ自動車株式会社 | 立体装置の温度及び濃度計測装置、燃焼機関の温度及び濃度計測装置、燃焼機関及び立体装置の温度及び濃度計測方法 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP2031374B1 (en) * | 2007-08-31 | 2012-10-10 | Canon Kabushiki Kaisha | Apparatus and method for obtaining information related to terahertz waves |
-
2021
- 2021-02-07 CN CN202110176882.7A patent/CN112985613B/zh active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103033282A (zh) * | 2012-10-29 | 2013-04-10 | 大连理工大学 | 一种检测磁约束聚变装置偏滤器石墨瓦瞬态温度的方法 |
| CN103048061A (zh) * | 2012-10-29 | 2013-04-17 | 大连理工大学 | 反射太赫兹谱技术检测偏滤器石墨瓦瞬态温度的装置 |
| CN103115893A (zh) * | 2013-01-30 | 2013-05-22 | 大连理工大学 | 一种检测托卡马克钨第一壁灰尘沉积层成分及厚度的装置 |
| JP2016142557A (ja) * | 2015-01-30 | 2016-08-08 | いすゞ自動車株式会社 | 立体装置の温度及び濃度計測装置、燃焼機関の温度及び濃度計測装置、燃焼機関及び立体装置の温度及び濃度計測方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112985613A (zh) | 2021-06-18 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Yuan et al. | Accuracy improvement of quantitative analysis for major elements in laser-induced breakdown spectroscopy using single-sample calibration | |
| CN108180999B (zh) | 基于激光扫描的红外探测器绝对响应度定标装置及方法 | |
| CN101566501B (zh) | 一种光纤光谱协同放电电流测量等离子体电子密度的方法 | |
| EP3136065B1 (en) | Method for measuring actual temperature of flame using all information about radiation spectrum and measurement system | |
| Fu et al. | Temperature measurements using multicolor pyrometry in thermal radiation heating environments | |
| CN108709871B (zh) | 一种基于tdlas的气体浓度检测方法 | |
| CN112985613B (zh) | 一种THz回波高温温度测量方法 | |
| CN108760635A (zh) | 高灵敏度的土壤重金属元素探测装置与探测方法 | |
| CN111289496A (zh) | 一种远距离变焦距激光诱导击穿光谱的检测方法及装置 | |
| CN106248653B (zh) | 一种提高激光诱导击穿光谱定量分析长期稳定性的方法 | |
| CN106768471A (zh) | 一种基于压力展宽的非接触式温度测量方法 | |
| Qingdong et al. | Portable fiber-optic laser-induced breakdown spectroscopy system for the quantitative analysis of minor elements in steel | |
| Xiaojian et al. | Atomic emission dual-spectrum thermometry for laser-induced Cu plasma temperature | |
| GB2602887A (en) | Thermal measurement system | |
| CN103792214A (zh) | 一种提高钢铁中碳元素含量测量精度的方法 | |
| CN112985610B (zh) | 一种THz回波高温温度测量装置 | |
| RU2752809C1 (ru) | Способ спектротермометрии | |
| Shu et al. | Research on the normal emissivity of graphite between 150 and 500° C by an infrared camera for nuclear fusion devices | |
| CN114623932A (zh) | 一种基于红外辐射成像的mpd阴极温度监测***及方法 | |
| Xu | Calibration technology and application of laser power meter | |
| CN114061766A (zh) | 一种颗粒燃烧过程多光谱重建测温装置与方法 | |
| Chase et al. | A search for the Sunyaev–Zeldovich effect at millimetre wavelengths | |
| CN204807232U (zh) | 温度测量*** | |
| Baolin et al. | Investigation of the effective wavelength applied on the active dual wavelength infrared temperature measurement | |
| Lu et al. | Quantitative analysis of Cr in soil using laser induced breakdown spectroscopy combined with relevance vector machine method |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |