CN107389220A - 一种锅炉内温度场分布的检测方法 - Google Patents
一种锅炉内温度场分布的检测方法 Download PDFInfo
- Publication number
- CN107389220A CN107389220A CN201710547693.XA CN201710547693A CN107389220A CN 107389220 A CN107389220 A CN 107389220A CN 201710547693 A CN201710547693 A CN 201710547693A CN 107389220 A CN107389220 A CN 107389220A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- ultrasonic transducer
- gas
- sound channel
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 51
- 238000001514 detection method Methods 0.000 title claims abstract description 34
- 239000007788 liquid Substances 0.000 claims abstract description 119
- 239000012530 fluid Substances 0.000 claims abstract description 50
- 238000002604 ultrasonography Methods 0.000 claims description 31
- 239000011159 matrix material Substances 0.000 claims description 15
- 238000004364 calculation method Methods 0.000 claims description 10
- 230000017105 transposition Effects 0.000 claims description 6
- 230000008569 process Effects 0.000 claims description 4
- 238000012545 processing Methods 0.000 claims description 4
- 230000002463 transducing effect Effects 0.000 claims description 4
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 3
- 239000007789 gas Substances 0.000 claims 27
- 238000003384 imaging method Methods 0.000 abstract description 6
- 238000009826 distribution Methods 0.000 description 18
- 230000000644 propagated effect Effects 0.000 description 8
- 238000010586 diagram Methods 0.000 description 5
- 230000005540 biological transmission Effects 0.000 description 4
- 238000007639 printing Methods 0.000 description 3
- 238000006243 chemical reaction Methods 0.000 description 2
- 238000005314 correlation function Methods 0.000 description 2
- 238000007689 inspection Methods 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 241000208340 Araliaceae Species 0.000 description 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 1
- 235000003140 Panax quinquefolius Nutrition 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000004888 barrier function Effects 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 239000000919 ceramic Substances 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000000052 comparative effect Effects 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000005284 excitation Effects 0.000 description 1
- 239000004744 fabric Substances 0.000 description 1
- OVBPIULPVIDEAO-LBPRGKRZSA-N folic acid Chemical compound C=1N=C2NC(N)=NC(=O)C2=NC=1CNC1=CC=C(C(=O)N[C@@H](CCC(O)=O)C(O)=O)C=C1 OVBPIULPVIDEAO-LBPRGKRZSA-N 0.000 description 1
- 235000008434 ginseng Nutrition 0.000 description 1
- 238000012806 monitoring device Methods 0.000 description 1
- 230000011514 reflex Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 230000008054 signal transmission Effects 0.000 description 1
- 238000001931 thermography Methods 0.000 description 1
- 230000026683 transduction Effects 0.000 description 1
- 238000010361 transduction Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01K—MEASURING TEMPERATURE; MEASURING QUANTITY OF HEAT; THERMALLY-SENSITIVE ELEMENTS NOT OTHERWISE PROVIDED FOR
- G01K11/00—Measuring temperature based upon physical or chemical changes not covered by groups G01K3/00, G01K5/00, G01K7/00 or G01K9/00
- G01K11/22—Measuring temperature based upon physical or chemical changes not covered by groups G01K3/00, G01K5/00, G01K7/00 or G01K9/00 using measurement of acoustic effects
Landscapes
- Physics & Mathematics (AREA)
- Acoustics & Sound (AREA)
- General Physics & Mathematics (AREA)
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
- Measuring Temperature Or Quantity Of Heat (AREA)
Abstract
本发明公开了一种锅炉内温度场分布的检测方法，涉及超声成像领域，包括步骤：提供锅炉内温度场分布的检测装置；激励液体超声换能器和气体超声换能器发射和接收超声波信号；计算锅炉内气体空间和液体空间的交界面；建立气体空间的温度场到直线声道平均飞行时间的投影方程：求解投影方程中气体空间各像素点的温度；建立液体空间的温度场到直线声道平均飞行时间的投影方程；求解投影方程中液体空间各像素点的温度。本发明提供的锅炉内温度场分布的检测方法，可通过超声成像的方法检测出锅炉内温度场分布，适用范围广。
Description
技术领域
本发明涉及超声成像领域，具体涉及一种锅炉内温度场分布的检测方法。
背景技术
锅炉内温度场直接影响锅炉的安全性，检测锅炉内温度场尤为重要。现有技术，通常采用红外热成像技术进行检测，利用红外探测器接收被测目标的红外辐射能量，实现物体表面温度场的成像，得到锅炉内温度场分布。
然而，上述方法仅适用于管壁透光的锅炉，无法检测出管壁不透光锅炉的内部温度场。
发明内容
针对现有技术中存在的缺陷，本发明的目的在于提供一种锅炉内温度场分布的检测方法，可通过超声成像的方法检测出锅炉内温度场分布，适用范围广。
为达到以上目的，本发明采取的技术方案是：
一种锅炉内温度场分布的检测方法，包括如下步骤：
提供一种锅炉内温度场分布的检测装置，所述检测装置包括：设于锅炉内管壁上的多个液体超声换能器和气体超声换能器，至少部分所述液体超声换能器和气体超声换能器交替布置，且所述液体超声换能器和气体超声换能器位于同一平面上；
激励液体超声换能器和气体超声换能器发射和接收超声波信号，处理所有液体超声换能器和气体超声换能器接收到的超声波信号，得到液体超声换能器和气体超声换能器超声波信号从发射到接收的飞行时间；
计算锅炉内气体空间和液体空间的交界面；
建立气体空间的温度场到直线声道平均飞行时间的投影方程：根据锅炉内气体空间和液体空间的交界面位置得到气体空间内直线声道总数为Mg，将气体空间划分为Ng个像素点，第r个像素点的温度记为Tr，气体空间的温度场到直线声道平均飞行时间的投影方程为：
AXg＝tg
式中，系统矩阵s＝1,2,...,Mg，r＝1,2,...,Ng，lrs为第s个直线声道与第r个像素点的相交长度，cg0为0℃时空气声速，T0＝273.15K。为第Mg个直线声道的平均飞行时间；
求解投影方程中气体空间各像素点的温度：根据已知的系统矩阵A和向量tg求解得到Xg，进而得到气体空间中各像素点的温度；
建立液体空间的温度场到直线声道平均飞行时间的投影方程：根据锅炉内气体空间和液体空间的交界面位置得到液体空间内直线声道总数为Mw，将液体空间划分为Nw个像素点，第r'个像素点的温度记为Tr'，液体空间的温度场到直线声道平均飞行时间的投影方程为：
BXw＝tw
式中
xr'＝1/[cw0+p(Tr'-T0)+q(Tr'-T0)2]，
br's'＝lr's'，s'＝1,2,...,Mw，r'＝1,2,...,Nw，lr's'为第s'个直线声道与第r'个像素点的相交长度，cg0为0℃时空气声速，T0＝273.15K，为第Mw个直线声道的平均飞行时间；T0＝273.15K，cw0为0℃时水中的声速，p，q为已知常数；
求解投影方程中液体空间各像素点的温度：根据已知的系统矩阵B和向量tw求解得到Xw，进而得到液体空间中各像素点的温度。
在上述技术方案的基础上，在计算锅炉内气体空间和液体空间的交界面之后，建立气体空间的温度场到直线声道平均飞行时间的投影方程之前，还包括步骤细分内插直线声道：在相邻的两个液体超声换能器、相邻的两个气体超声换能器之间细分内插多条直线声道，并计算内插的直线声道上超声波信号的平均飞行时间。
在上述技术方案的基础上，细分内插直线声道的具体过程包括：测量所有液体超声换能器和气体超声换能器对应的直线声道与管壁切线的夹角，在相邻的两个液体超声换能器、相邻的两个气体超声换能器之间对应细分内插多条直线声道，并计算内插的直线声道上超声波信号的平均飞行时间t，公式如下：
式中，hv＝αv-αv-1，u＝(α-αv-1)/hv，αv-1＜α＜αv，v＝1,2,...,V0，α为第v个声道和第v-1个声道之间的内插声道与管壁切线的夹角，αv为第v声道与管壁切线的夹角，αv-1为第v-1个声道与管壁切线的夹角，S”(αj)为拟合曲线S(α)的二阶导数，V0为所有液体超声换能器和气体超声换能器对应的直线声道的个数。
在上述技术方案的基础上，求解投影方程中气体空间各像素点的温度的具体步骤包括：采用迭代法求解Xg，计算公式如下：
Xg (b+1)＝Xg (b)+λAT(tg-AXg (b))
式中，AT为系统矩阵A的转置，λ为迭代步长，Xg (b+1)与Xg (b)分别为第b+1步和第b步的迭代结果；
当||Xg (b+1)-Xg (b)||2小于预设值时，迭代结束，得到Xg。
在上述技术方案的基础上，求解投影方程中液体空间各像素点的温度的具体步骤包括：采用迭代法求解Xw，计算公式如下：
Xw (b'+1)＝Xw (b')+λ'BT(tw-BXw (b'))
式中，BT为系统矩阵B的转置，λ'为迭代步长，Xw (b'+1)与Xw (b')分别为第b'+1步和第b'步的迭代结果；
当||Xw (b'+1)-Xw (b')||2小于预设值时，迭代结束，得到Xw。
在上述技术方案的基础上，所述液体超声换能器和气体超声换能器结构相同，所述液体超声换能器包括一半圆筒状的收发体和U型外壳，所述收发体包括位于外侧的匹配层和位于内侧的压电陶瓷层，所述收发体和外壳形成一空腔体，所述空腔体内填满背衬，所述液体超声换能器的压电陶瓷层的中心谐振频率和气体超声换能器的压电陶瓷层的中心谐振频率不同。
在上述技术方案的基础上，所述液体超声换能器的压电陶瓷层的中心谐振频率为1MHz～5MHz。
在上述技术方案的基础上，所述气体超声换能器的压电陶瓷层的中心谐振频率为100KHz～300KHz。
在上述技术方案的基础上，所有所述液体超声换能器和气体超声换能器交替布置，相邻的液体超声换能器和气体超声换能器间距相等。
在上述技术方案的基础上，所述液体超声换能器和气体超声换能器所在的平面经过锅炉的中心点。
与现有技术相比，本发明的优点在于：
(1)本发明的锅炉内温度场分布的检测方法可通过超声成像的方法检测出锅炉内温度场分布，既可适用于管壁透光的情况，也可以管壁不透光的情况，适用范围广；超声波信号在介质中可以传播较远的距离，监测装置的测量尺度大，可同时适用于小型锅炉和大型锅炉。
(2)本发明的气体超声换能器和液体超声换能器均包括一半圆筒状的收发体，使得换能器在横截面上具有180度的发射角和接收角，发射和接收信号广度好，而且换能器在纵截面内的发射角和接收角小于30°，从而使得换能器仅能在换能器的径向上发射和接收信号，有效避免其他平面多重反射信号的干扰。
附图说明
图1为本发明实施例中锅炉内气液二相分布场及温度场分布的检测装置的结构示意图；
图2为本发明实施例中液体超声换能器的结构示意图；
图3为图2的剖视图；
图4本发明实施例中锅炉内气液二相分布场的检测方法流程图；
图5为本发明实施例中计算锅炉内气体空间和液体空间的交界面位置的详细流程图；
图6为本发明实施例中水平声道示意图；
图7为本发明实施例中竖直声道示意图；
图8为本发明实施例中计算液体空间二相分布场的详细流程图；
图9为本发明实施例中一次反射示意图；
图10为本发明实施例中锅炉内温度场分布的检测方法流程图；
图11为本发明实施例中内插前直线声道的分布示意图；
图12为本发明实施例中内插后直线声道的分布示意图。
图中：1-液体超声换能器，2-气体超声换能器，11-收发体，111-匹配层，112-压电陶瓷层，12-外壳，13-空腔体，3-气体空间，4-液体空间，5-水平声道，6-竖直声道，7-气泡。
具体实施方式
以下结合附图及实施例对本发明作进一步详细说明。
参见图1所示，本发明实施例提供一种锅炉内气液二相分布场及温度场分布的检测装置，包括：设于锅炉内管壁上的多个液体超声换能器1和气体超声换能器2，至少部分液体超声换能器1和气体超声换能器2交替布置，且液体超声换能器1和气体超声换能器2位于同一平面上，该平面经过锅炉的中心点。当液体超声换能器1和气体超声换能器2所在的平面经过锅炉的中心点时，探测平面最大，最能反映锅炉内气液二相分布场及温度场分布。
本发明实施例中，所有液体超声换能器1和气体超声换能器2交替布置，相邻的液体超声换能器1和气体超声换能器2间距相等。
参见图2和图3所示，液体超声换能器1和气体超声换能器2结构相同，液体超声换能器1包括一半圆筒状的收发体11和U型外壳12，收发体11包括位于外侧的匹配层111和位于内侧的压电陶瓷层112，收发体11和外壳12形成一空腔体13，空腔体13内填满背衬，液体超声换能器1的压电陶瓷层112的中心谐振频率和气体超声换能器2的压电陶瓷层112的中心谐振频率不同。
本发明实施例中，液体超声换能器1的压电陶瓷层112的中心谐振频率为1MHz～5MHz，气体超声换能器2的压电陶瓷层112的中心谐振频率为100KHz～300KHz。液体超声换能器1和气体超声换能器2分别用于液体空间4和气体空间3内发射或接收超声波信号，液体超声换能器1和气体超声换能器2均为半圆柱形且收发一体式。
本发明实施例中匹配层111和压电陶瓷层112均为半圆筒状，压电陶瓷层112横截面为半圆弧状，使得换能器在横截面上具有180度的发射角和接收角，发射和接收信号广度好；压电陶瓷层112纵截面为矩形结构，在收发体11的轴线方向具有一定长度，使得换能器在纵截面内的发射角和接收角小于30°，从而使得换能器仅能在换能器的径向上发射和接收信号，有效避免其他平面多重反射信号的干扰。
在安装上述半圆筒状的液体超声换能器1和气体超声换能器2时，使液体超声换能器1和气体超声换能器2的径向和锅炉的纵剖面平行，使得液体超声换能器1和气体超声换能器2仅能在同一平面上发射和接收信号，信号传输质量更好。
参见图4所示，本发明还提供了一种锅炉内气液二相分布场的检测方法，包括以下步骤：
S1：激励液体超声换能器1和气体超声换能器2发射和接收超声波信号；
S2：处理所有液体超声换能器1和气体超声换能器2接收到的超声波信号，得到液体超声换能器1和气体超声换能器2超声波信号的幅值、以及超声波信号从发射到接收的飞行时间；
S3：计算超声波在气体空间3的平均声速，并根据气体超声换能器2的位置、以及气体超声换能器2的超声波沿反射路径传播的飞行时间，计算锅炉内气体空间3和液体空间4的交界面位置；
S4：测量并比较液体超声换能器1在无气泡7阻隔的液体环境下和待测液体空间4内接收到的超声波信号的幅值，得到待测液体空间4中气泡7的大小和位置，计算锅炉内液体空间4二相分布场。
步骤S1的具体步骤包括：每次激励一液体超声换能器1和一气体超声换能器2发射超声波信号，所有液体超声换能器1和气体超声换能器2接收超声波信号，直至激励完所有液体超声换能器1和气体超声换能器2。
由于液体超声换能器1和气体超声换能器2的中心谐振频率差异很大，互不干扰，液体超声换能器1和气体超声换能器2超声波信号的发射或接收可以同时进行，可节约探测时间，提高探测效率。然而，为了避免多个同类型超声波信号相互干扰，控制电路每次只能激励一个液体超声换能器1和一个气体超声换能器2发射超声波信号。
步骤S2中得到超声波信号的幅值的具体步骤包括：采用高速AD对液体超声换能器1和气体超声换能器2接收到的超声波信号进行采样，根据超声波的气体声速和液体声速获得超声波信号的到达时间窗，并根据到达时间窗内采集到超声波信号的最大峰峰值确定接收到的超声波信号的幅值。
气体超声换能器2发射的超声波信号只能在气体中飞行，液体超声换能器1发射的超声波信号只能在液体中飞行，超声波信号在气体和液体中飞行时间的测量方法有所不同。
步骤S2中得到气体超声换能器2超声波信号从发射到接收的飞行时间的具体步骤包括：
将气体超声换能器2在已知温度、介质分布的条件下接收到的超声波信号作为参考信号；
对气体超声换能器2接收的超声波信号和参考信号进行互相关处理，得到互相关函数，然后对互相关函数进行Hilbert变换并插值求取过零点获得亚采样精度的飞行时间。
超声波在液体环境下从发射到接收的过程中，锅炉管壁会产生耦合回波，且存在气泡7反射回波等多种传播途径产生的超声波信号干扰，导致采用高速AD在采集液体超声换能器1超声波信号时会采集到很多干扰信号。
步骤S2中得到液体超声换能器1超声波信号从发射到接收的飞行时间的具体步骤包括：
根据已得到的超声波信号的幅值和到达时间窗，从采集到的超声波信号中识别出直线传输的超声波信号，作为目标超声波信号，并利用数字相关法获得目标超声波信号的粗略飞行时间；同时，
利用模拟电路识别目标超声波信号的多个过零点，并通过时间数字转换法高精度测量每个过零点对应的时间间隔；然后，根据已获得的回波信号的粗略飞行时间从所有时间间隔中识别出目标超声波信号的精确飞行时间。
参见图5所示，步骤S3中计算锅炉内气体空间3和液体空间4的交界面位置的具体步骤包括：
S301：参见图6所示，根据所有能接收到超声波信号的气体超声换能器2位置分布，将气体超声换能器2水平传输超声波信号的传输路径作为水平声道5，并获得水平声道5的个数为N1、以及每个水平声道5的长度，每个水平声道5对应两个气体超声换能器2，将这两个气体超声换能器2对向发射和接收的两个超声波信号飞行时间分为超声波正向传播飞行时间和超声波反向传播飞行时间，并从所有液体超声换能器1和气体超声换能器2接收到的超声波信号中筛选每个水平声道5的超声波正向传播飞行时间和超声波反向传播飞行时间；
计算每个水平声道5上的声速，公式如下：
式中，k＝1,2,...,N1，Lk为第k个水平声道5的长度，为超声波正向传播飞行时间，为超声波反向传播飞行时间，ck为第k个水平声道5上的声速；
利用加权平均法计算气体空间3的平均声速计算公式如下：
其中，Lk为第k个水平声道5的长度，ck为第k个水平声道5上的声速，为气体空间3的平均声速；
S302：参见图7所示，根据所有能接收到超声波信号的气体超声换能器2位置分布，将气体超声换能器2竖直传输超声波信号的传输路径作为竖直声道6，并获得竖直声道6的个数为n，每个竖直声道6对应一个气体超声换能器2，从所有接收到的超声波信号中筛选出每个竖直声道6对应的气体超声换能器2从发射经气液界面反射后接收的飞行时间，即为超声波沿反射路径传播的飞行时间，并以锅炉底部中心为坐标原点建立直角坐标系，得到每个气体超声换能器2的横坐标和纵坐标；
计算气体空间3和液体空间4的交界面上若干坐标点Xn,y(Xn)，当坐标点的横坐标为Xn时，纵坐标y(Xn)的计算公式如下：
式中，Xn为第n个竖直声道6对应的气体超声换能器2的横坐标，Yn为第n个竖直声道6对应的气体超声换能器2的纵坐标，为第n个竖直声道6对应的超声波沿反射路径传播的飞行时间，为气体空间3的平均声速；
S303：连接气体空间3和液体空间4的交界面上所有的坐标点，即为锅炉内气体空间3和液体空间4的交界面位置。
参见图8所示，步骤S4的计算液体空间4二相分布场的具体步骤包括：
S401：在液体超声换能器i和液体超声换能器j均处于液体状态，且该液体内没有气泡7的情况下，由液体超声换能器i发射超声波信号，并由液体超声换能器j接收超声波信号，得到液体超声换能器j在无气泡7阻隔的液体环境下接收到的沿直线路径传播的超声波信号的幅值；
S402：在待测的液体空间4内，液体超声换能器i发射超声波，液体超声换能器j接收该超声波，得到待测的液体空间4中液体超声换能器j接收到的沿直线路径传播的超声波信号的幅值，该直线传播路径为直线声道Lij，并计算液体超声换能器i和液体超声换能器j之间超声波直线路径传播的阻塞系数Uij，计算公式为：
式中，A0为无气泡7阻隔时液体超声换能器j接收到的沿直线路径传播的超声波信号的幅值，A为待测的气液二相分布场中液体超声换能器j接收到的沿直线路径传播的超声波信号的幅值，a为阈值常数，根据经验选取a为0.5；Uij表征了液体超声换能器i和液体超声换能器j之间的直线声道Lij上是否存在气泡7阻隔，当A和A0的比值大于设定的阈值a时，则表明不存在气泡7阻隔，当A和A0的比值小于设定的阈值a时，则表明存在气泡7阻隔；
S403：根据经验确定液体空间4的分辨率，并结合坐标系得到液体空间4的所有像素点坐标(x,y)，计算像素点(x,y)与直线声道Lij的相交系数Sij(x,y)，当像素点(x,y)被声道Lij穿过时Sij(x,y)取1，否则取为0；
S404：参见图9所示，在液体超声换能器i和液体超声换能器j均处于液体状态，且该液体内没有气泡7的情况下，由液体超声换能器i发射超声波信号，并由液体超声换能器j接收超声波信号，得到无气泡7阻隔时液体超声换能器j接收到的沿一次反射路径传播的超声波信号的幅值；
S405：在待测的液体空间4内，液体超声换能器i发射超声波，液体超声换能器j接收该超声波，得到待测的液体空间4内液体超声换能器j接收到的沿一次反射路径传播的超声波信号的幅值，该一次反射传播路径为一次反射声道L'ij，并计算液体超声换能器i和液体超声换能器j之间沿一次反射路径传播的超声脉冲的阻塞系数U′ij，计算公式为：
式中，A′0为无气泡7阻隔时液体超声换能器j接收到的沿一次反射路径传播的超声波信号的幅值，A'为待测的气液二相分布场中液体超声换能器j接收到的沿一次反射路径传播的超声波信号的幅值，a'为阈值常数，根据经验选取a'为0.5；U′ij表征了换能器i和换能器j之间的一次反射声道L'ij上是否存在气泡7阻隔，当A'和A'0的比值大于设定的阈值a'时，则表明不存在气泡7阻隔，当A'和A'0的比值小于设定的阈值a'时，则表明存在气泡7阻隔；
S406：根据已得到的气体空间3和液体空间4的交界面位置并利用镜像法确定一次反射声道L′ij的几何位置，计算像素点(x,y)与一次反射声道L′ij的相交系数S′ij(x,y)，当像素点(x,y)被声道L′ij穿过时S′ij(x,y)取1，否则取为0；
S407：采用的反投影重建公式计算像素点(x,y)处为气泡7的概率为：
式中，N2为液体空间4内的液体超声换能器的个数；
S408：基于气液二相分布场的二值特性，对g(x,y)进行二值化处理，公式如下：
当g'(x,y)＝0时，表明像素点(x,y)处为液体，当g'(x,y)＝1时表明像素点(x,y)处为气体，进而得到锅炉内气液二相分布场。
参见图10所示，本发明还提供了一种锅炉内温度场分布的检测方法，包括如下步骤：
A1：提供一种锅炉内温度场分布的检测装置，检测装置包括：设于锅炉内管壁上的多个液体超声换能器1和气体超声换能器2，至少部分液体超声换能器1和气体超声换能器2交替布置，且液体超声换能器1和气体超声换能器2位于同一平面上；
A2：激励液体超声换能器1和气体超声换能器2发射和接收超声波信号，处理所有液体超声换能器1和气体超声换能器2接收到的超声波信号，得到液体超声换能器1和气体超声换能器2超声波信号从发射到接收的飞行时间；
A3：计算锅炉内气体空间3和液体空间4的交界面；
A4：细分内插直线声道：参见图11和图12所示，测量所有液体超声换能器1和气体超声换能器2对应的直线声道与管壁切线的夹角，在相邻的两个液体超声换能器1、相邻的两个气体超声换能器2之间对应细分内插多条直线声道，并计算内插的直线声道上超声波信号的平均飞行时间t，公式如下：
式中，hv＝αv-αv-1，u＝(α-αv-1)/hv，αv-1＜α＜αv，v＝1,2,...,V0，α为第v个声道和第v-1个声道之间的内插声道与管壁切线的夹角，αv为第v声道与管壁切线的夹角，αv-1为第v-1个声道与管壁切线的夹角，S”(αj)为拟合曲线S(α)的二阶导数，V0为所有液体超声换能器1和气体超声换能器2对应的直线声道的个数；
A5：建立气体空间3的温度场到直线声道平均飞行时间的投影方程：根据锅炉内气体空间3和液体空间4的交界面位置得到内插后气体空间3内直线声道总数为Mg，将气体空间3划分为Ng个像素点，第r个像素点的温度记为Tr，气体空间3的温度场到直线声道平均飞行时间的投影方程为：
AXg＝tg
式中，系统矩阵s＝1,2,...,Mg，r＝1,2,...,Ng，lrs为第s个直线声道与第r个像素点的相交长度，cg0为0℃时空气声速，T0＝273.15K。为第Mg个直线声道的平均飞行时间；
A6：求解投影方程中气体空间3各像素点的温度：根据已知的系统矩阵A和向量tg，采用迭代法求解Xg，计算公式如下：
Xg (b+1)＝Xg (b)+λAT(tg-AXg (b))
式中，AT为系统矩阵A的转置，λ为迭代步长，Xg (b+1)与Xg (b)分别为第b+1步和第b步的迭代结果；
当||Xg (b+1)-Xg (b)||2小于预设值时，迭代结束，得到Xg，进而得到气体空间3中每个像素点处的温度，公式如下：
式中，r＝1,2,...,Ng，Tr为第r个像素点的温度；
A7：建立液体空间4的温度场到直线声道平均飞行时间的投影方程：根据锅炉内气体空间3和液体空间4的交界面位置得到内插后液体空间4内直线声道总数为Mw，将液体空间4划分为Nw个像素点，第r'个像素点的温度记为Tr'，得到液体空间4的温度场到直线声道平均飞行时间的投影方程为：
BXw＝tw
式中
xr'＝1/[cw0+p(Tr'-T0)+q(Tr'-T0)2]，
br's'＝lr's'，s'＝1,2,...,Mw，r'＝1,2,...,Nw，lr's'为第s'个直线声道与第r'个像素点的相交长度，cg0为0℃时空气声速，T0＝273.15K，为第Mw个直线声道的平均飞行时间；T0＝273.15K，cw0为0℃时水中的声速，p，q为已知常数；
A8：求解投影方程中液体空间4各像素点的温度：根据已知的系统矩阵B和向量tw，采用迭代法求解Xw，计算公式如下：
Xw (b'+1)＝Xw (b')+λ'BT(tw-BXw (b'))
式中，BT为系统矩阵B的转置，λ'为迭代步长，Xw (b'+1)与Xw (b')分别为第b'+1步和第b'步的迭代结果；
当||Xw (b'+1)-Xw (b')||2小于预设值时，迭代结束，得到Xw，进而得到液体空间4中每个像素点处的温度，公式如下：
式中，r'＝1,2,...,Nw，Tr'为第r'个像素点处的温度。
本发明实施例中的锅炉内气液二相分布场及温度场分布的检测装置，一方面，既可适用于管壁不导电的情况，也可以管壁导电的情况，另一方面，既可适用于管壁透光的情况，也可以管壁不透光的情况，适用范围广。
本发明不局限于上述实施方式，对于本技术领域的普通技术人员来说，在不脱离本发明原理的前提下，还可以做出若干改进和润饰，这些改进和润饰也视为本发明的保护范围之内。本说明书中未作详细描述的内容属于本领域专业技术人员公知的现有技术。
Claims (10)
1.一种锅炉内温度场分布的检测方法，其特征在于，包括如下步骤：
提供一种锅炉内温度场分布的检测装置，所述检测装置包括：设于锅炉内管壁上的多个液体超声换能器和气体超声换能器，至少部分所述液体超声换能器和气体超声换能器交替布置，且所述液体超声换能器和气体超声换能器位于同一平面上；
激励液体超声换能器和气体超声换能器发射和接收超声波信号，处理所有液体超声换能器和气体超声换能器接收到的超声波信号，得到液体超声换能器和气体超声换能器超声波信号从发射到接收的飞行时间；
计算锅炉内气体空间和液体空间的交界面；
建立气体空间的温度场到直线声道平均飞行时间的投影方程：根据锅炉内气体空间和液体空间的交界面位置得到气体空间内直线声道总数为Mg，将气体空间划分为Ng个像素点，第r个像素点的温度记为Tr，气体空间的温度场到直线声道平均飞行时间的投影方程为：
AXg＝tg
式中，系统矩阵lrs为第s个直线声道与第r个像素点的相交长度，cg0为0℃时空气声速，T0＝273.15K。为第Mg个直线声道的平均飞行时间；
求解投影方程中气体空间各像素点的温度：根据已知的系统矩阵A和向量tg求解得到Xg，进而得到气体空间中各像素点的温度；
建立液体空间的温度场到直线声道平均飞行时间的投影方程：根据锅炉内气体空间和液体空间的交界面位置得到液体空间内直线声道总数为Mw，将液体空间划分为Nw个像素点，第r'个像素点的温度记为Tr'，液体空间的温度场到直线声道平均飞行时间的投影方程为：
BXw＝tw
式中
xr'＝1/[cw0+p(Tr'-T0)+q(Tr'-T0)2]，
br's'＝lr's'，s'＝1,2,...,Mw，r'＝1,2,...,Nw，lr's'为第s'个直线声道与第r'个像素点的相交长度，cg0为0℃时空气声速，T0＝273.15K，为第Mw个直线声道的平均飞行时间；T0＝273.15K，cw0为0℃时水中的声速，p，q为已知常数；
求解投影方程中液体空间各像素点的温度：根据已知的系统矩阵B和向量tw求解得到Xw，进而得到液体空间中各像素点的温度。
2.如权利要求1所述的锅炉内温度场分布的检测方法，其特征在于，在计算锅炉内气体空间和液体空间的交界面之后，建立气体空间的温度场到直线声道平均飞行时间的投影方程之前，还包括步骤细分内插直线声道：在相邻的两个液体超声换能器、相邻的两个气体超声换能器之间细分内插多条直线声道，并计算内插的直线声道上超声波信号的平均飞行时间。
3.如权利要求2所述的锅炉内温度场分布的检测方法，其特征在于，细分内插直线声道的具体过程包括：测量所有液体超声换能器和气体超声换能器对应的直线声道与管壁切线的夹角，在相邻的两个液体超声换能器、相邻的两个气体超声换能器之间对应细分内插多条直线声道，并计算内插的直线声道上超声波信号的平均飞行时间t，公式如下：
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>t</mi>
<mo>=</mo>
<mi>S</mi>
<mrow>
<mo>(</mo>
<mi>&alpha;</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>t</mi>
<mrow>
<mi>v</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>&CenterDot;</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<mi>u</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>t</mi>
<mi>v</mi>
</msub>
<mo>&CenterDot;</mo>
<mi>u</mi>
<mo>-</mo>
<mfrac>
<mn>1</mn>
<mn>6</mn>
</mfrac>
<msup>
<mi>S</mi>
<mrow>
<mo>&prime;</mo>
<mo>&prime;</mo>
</mrow>
</msup>
<mrow>
<mo>(</mo>
<msub>
<mi>&alpha;</mi>
<mrow>
<mi>v</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mi>u</mi>
<mo>-</mo>
<mn>3</mn>
<msup>
<mi>u</mi>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mi>u</mi>
<mn>3</mn>
</msup>
<mo>)</mo>
</mrow>
<msubsup>
<mi>h</mi>
<mi>v</mi>
<mn>2</mn>
</msubsup>
<mo>-</mo>
<mfrac>
<mn>1</mn>
<mn>6</mn>
</mfrac>
<msup>
<mi>S</mi>
<mrow>
<mo>&prime;</mo>
<mo>&prime;</mo>
</mrow>
</msup>
<mrow>
<mo>(</mo>
<msub>
<mi>&alpha;</mi>
<mi>v</mi>
</msub>
<mo>)</mo>
</mrow>
<mrow>
<mo>(</mo>
<mi>u</mi>
<mo>-</mo>
<msup>
<mi>u</mi>
<mn>3</mn>
</msup>
<mo>)</mo>
</mrow>
<msubsup>
<mi>h</mi>
<mi>v</mi>
<mn>2</mn>
</msubsup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mfrac>
<mn>1</mn>
<mn>6</mn>
</mfrac>
<msup>
<mi>S</mi>
<mrow>
<mo>&prime;</mo>
<mo>&prime;</mo>
</mrow>
</msup>
<mrow>
<mo>(</mo>
<msub>
<mi>&alpha;</mi>
<mrow>
<mi>v</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<msub>
<mi>h</mi>
<mi>v</mi>
</msub>
<mo>-</mo>
<mfrac>
<mn>1</mn>
<mn>3</mn>
</mfrac>
<msup>
<mi>S</mi>
<mrow>
<mo>&prime;</mo>
<mo>&prime;</mo>
</mrow>
</msup>
<mrow>
<mo>(</mo>
<msub>
<mi>&alpha;</mi>
<mi>v</mi>
</msub>
<mo>)</mo>
</mrow>
<mrow>
<mo>(</mo>
<msub>
<mi>h</mi>
<mi>v</mi>
</msub>
<mo>+</mo>
<msub>
<mi>h</mi>
<mrow>
<mi>v</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mfrac>
<mn>1</mn>
<mn>6</mn>
</mfrac>
<msup>
<mi>S</mi>
<mrow>
<mo>&prime;</mo>
<mo>&prime;</mo>
</mrow>
</msup>
<mrow>
<mo>(</mo>
<msub>
<mi>&alpha;</mi>
<mrow>
<mi>v</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<msub>
<mi>h</mi>
<mrow>
<mi>v</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>t</mi>
<mrow>
<mi>v</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>t</mi>
<mi>v</mi>
</msub>
</mrow>
<msub>
<mi>h</mi>
<mrow>
<mi>v</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
</mfrac>
<mo>-</mo>
<mfrac>
<mrow>
<msub>
<mi>t</mi>
<mi>v</mi>
</msub>
<mo>-</mo>
<msub>
<mi>t</mi>
<mrow>
<mi>v</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
</mrow>
<msub>
<mi>h</mi>
<mi>v</mi>
</msub>
</mfrac>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msup>
<mi>S</mi>
<mrow>
<mo>&prime;</mo>
<mo>&prime;</mo>
</mrow>
</msup>
<mrow>
<mo>(</mo>
<msub>
<mi>&alpha;</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mn>0</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msup>
<mi>S</mi>
<mrow>
<mo>&prime;</mo>
<mo>&prime;</mo>
</mrow>
</msup>
<mrow>
<mo>(</mo>
<msub>
<mi>&alpha;</mi>
<msub>
<mi>V</mi>
<mn>0</mn>
</msub>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mn>0</mn>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
式中，hv＝αv-αv-1，u＝(α-αv-1)/hv，αv-1＜α＜αv，v＝1,2,...,V0，α为第v个声道和第v-1个声道之间的内插声道与管壁切线的夹角，αv为第v声道与管壁切线的夹角，αv-1为第v-1个声道与管壁切线的夹角，S”(αj)为拟合曲线S(α)的二阶导数，V0为所有液体超声换能器和气体超声换能器对应的直线声道的个数。
4.如权利要求1所述的锅炉内温度场分布的检测方法，其特征在于，求解投影方程中气体空间各像素点的温度的具体步骤包括：采用迭代法求解Xg，计算公式如下：
Xg (b+1)＝Xg (b)+λAT(tg-AXg (b))
式中，AT为系统矩阵A的转置，λ为迭代步长，Xg (b+1)与Xg (b)分别为第b+1步和第b步的迭代结果；
当||Xg (b+1)-Xg (b)||2小于预设值时，迭代结束，得到Xg。
5.如权利要求1所述的锅炉内温度场分布的检测方法，其特征在于，求解投影方程中液体空间各像素点的温度的具体步骤包括：采用迭代法求解Xw，计算公式如下：
Xw (b'+1)＝Xw (b')+λ'BT(tw-BXw (b'))
式中，BT为系统矩阵B的转置，λ'为迭代步长，Xw (b'+1)与Xw (b')分别为第b'+1步和第b'步的迭代结果；
当||Xw (b'+1)-Xw (b')||2小于预设值时，迭代结束，得到Xw。
6.如权利要求1所述的锅炉内温度场分布的检测方法，其特征在于：所述液体超声换能器和气体超声换能器结构相同，所述液体超声换能器包括一半圆筒状的收发体和U型外壳，所述收发体包括位于外侧的匹配层和位于内侧的压电陶瓷层，所述收发体和外壳形成一空腔体，所述空腔体内填满背衬，所述液体超声换能器的压电陶瓷层的中心谐振频率和气体超声换能器的压电陶瓷层的中心谐振频率不同。
7.如权利要求6所述的锅炉内温度场分布的检测装置，其特征在于：所述液体超声换能器的压电陶瓷层的中心谐振频率为1MHz～5MHz。
8.如权利要求6所述的锅炉内温度场分布的检测装置，其特征在于：所述气体超声换能器的压电陶瓷层的中心谐振频率为100KHz～300KHz。
9.如权利要求1所述的锅炉内温度场分布的检测方法，其特征在于：所有所述液体超声换能器和气体超声换能器交替布置，相邻的液体超声换能器和气体超声换能器间距相等。
10.如权利要求1所述的锅炉内温度场分布的检测方法，其特征在于：所述液体超声换能器和气体超声换能器所在的平面经过锅炉的中心点。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710547693.XA CN107389220B (zh) | 2017-07-06 | 2017-07-06 | 一种锅炉内温度场分布的检测方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710547693.XA CN107389220B (zh) | 2017-07-06 | 2017-07-06 | 一种锅炉内温度场分布的检测方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107389220A true CN107389220A (zh) | 2017-11-24 |
| CN107389220B CN107389220B (zh) | 2019-07-16 |
Family
ID=60335596
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710547693.XA Expired - Fee Related CN107389220B (zh) | 2017-07-06 | 2017-07-06 | 一种锅炉内温度场分布的检测方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107389220B (zh) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110873612A (zh) * | 2018-08-30 | 2020-03-10 | 恩智浦有限公司 | 用于确定温度的设备 |
| WO2021217407A1 (zh) * | 2020-04-28 | 2021-11-04 | 华为技术有限公司 | 温度测量系统和方法 |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4761650A (en) * | 1985-12-06 | 1988-08-02 | Communications Research Laboratory, Ministry Of Posts And Telecommunications | System for measuring height distributions of atmospheric temperature, wind direction and wind speed |
| CN101699236A (zh) * | 2009-11-19 | 2010-04-28 | 杭州电子科技大学 | 深海热液口温度场原位在线声学检测方法 |
| CN202075062U (zh) * | 2011-04-13 | 2011-12-14 | 沈阳航空航天大学 | 基于声波传感器的炉膛温度场与炉管泄漏一体化检测装置 |
| CN103336145A (zh) * | 2013-06-18 | 2013-10-02 | 清华大学 | 一种针对管道内流体轴向流场的超声成像方法及设备 |
| CN104132749A (zh) * | 2014-08-25 | 2014-11-05 | 哈尔滨理大晟源科技开发有限公司 | 声波测定炉膛温度装置及测定方法 |
| CN105486424A (zh) * | 2014-09-17 | 2016-04-13 | 南京理工大学 | 燃烧室内壁瞬态温度场的超声波非侵入式测量方法 |
| CN106446437A (zh) * | 2016-10-10 | 2017-02-22 | 广东电网有限责任公司电力科学研究院 | 电站锅炉煤粉在变燃烧器出力下燃烧特性模拟方法及装置 |
-
2017
- 2017-07-06 CN CN201710547693.XA patent/CN107389220B/zh not_active Expired - Fee Related
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4761650A (en) * | 1985-12-06 | 1988-08-02 | Communications Research Laboratory, Ministry Of Posts And Telecommunications | System for measuring height distributions of atmospheric temperature, wind direction and wind speed |
| CN101699236A (zh) * | 2009-11-19 | 2010-04-28 | 杭州电子科技大学 | 深海热液口温度场原位在线声学检测方法 |
| CN202075062U (zh) * | 2011-04-13 | 2011-12-14 | 沈阳航空航天大学 | 基于声波传感器的炉膛温度场与炉管泄漏一体化检测装置 |
| CN103336145A (zh) * | 2013-06-18 | 2013-10-02 | 清华大学 | 一种针对管道内流体轴向流场的超声成像方法及设备 |
| CN104132749A (zh) * | 2014-08-25 | 2014-11-05 | 哈尔滨理大晟源科技开发有限公司 | 声波测定炉膛温度装置及测定方法 |
| CN105486424A (zh) * | 2014-09-17 | 2016-04-13 | 南京理工大学 | 燃烧室内壁瞬态温度场的超声波非侵入式测量方法 |
| CN106446437A (zh) * | 2016-10-10 | 2017-02-22 | 广东电网有限责任公司电力科学研究院 | 电站锅炉煤粉在变燃烧器出力下燃烧特性模拟方法及装置 |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110873612A (zh) * | 2018-08-30 | 2020-03-10 | 恩智浦有限公司 | 用于确定温度的设备 |
| WO2021217407A1 (zh) * | 2020-04-28 | 2021-11-04 | 华为技术有限公司 | 温度测量系统和方法 |
| CN114127526A (zh) * | 2020-04-28 | 2022-03-01 | 华为技术有限公司 | 温度测量系统和方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107389220B (zh) | 2019-07-16 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| EP3268701B1 (en) | Hybrid sensing ultrasonic flowmeter | |
| CN100455999C (zh) | 一种超声波测量液位的装置及方法 | |
| US6575043B1 (en) | Method and apparatus for characterizing flows based on attenuation of in-wall propagating wave modes | |
| Zhang et al. | Efficient immersion imaging of components with nonplanar surfaces | |
| CN103842779B (zh) | 超声流量计 | |
| CN113316705B (zh) | 用于确定细长或延伸结构的厚度的方法和系统 | |
| Fan et al. | Non-contact ultrasonic gas flow metering using air-coupled leaky Lamb waves | |
| JP5235112B2 (ja) | 形状検出装置、形状検出方法、及び形状検出装置を用いた骨強度診断装置 | |
| EP2816327B1 (en) | Ultrasonic flowmeter | |
| US10890471B2 (en) | Method and assembly for ultrasonic clamp-on flow measurement, and bodies for implementing off-center flow measurement | |
| JP2001356034A (ja) | 超音波流量測定方法及び超音波流量測定装置 | |
| CN105181793A (zh) | 基于超声波双频信号测量两相流含气率的方法 | |
| CN107389220B (zh) | 一种锅炉内温度场分布的检测方法 | |
| CN114088151A (zh) | 外夹式多声道超声波流量检测装置及检测方法 | |
| CN100380101C (zh) | 多普勒型超声波流量计 | |
| JP7006354B2 (ja) | 計測装置 | |
| JP3814347B2 (ja) | 流体量測定システム | |
| EP2419722B1 (en) | Method and apparatus for investigating a flexible hollow element using acoustic sensors | |
| KR101513697B1 (ko) | 파이프 두께 측정이 가능한 초음파 변환 장치 및 이를 이용한 유속 측정 장치 | |
| CN107478716A (zh) | 一种锅炉内气液二相分布场的检测装置及检测方法 | |
| JP3136002B2 (ja) | 超音波流量計 | |
| JP2001337077A (ja) | コンクリート構造物の剥離の非破壊検査方法 | |
| RU2687086C1 (ru) | Способ ультразвукового контроля толщины стенки трубопровода | |
| CN111220709B (zh) | 一种用于管道内超声相控阵成像的声束偏转时延控制方法 | |
| JP2001201487A (ja) | 超音波を用いたコンクリート構造物の非破壊検査方法及び装置 |
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 | ||
| CF01 | Termination of patent right due to non-payment of annual fee | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20190716 |