本发明属于温度测量技术领域，尤其涉及一种基于声学技术的空冷岛翅片管束温度场测量装置及其方法。
背景技术：
空冷是指利用空气冷却电站汽轮机乏汽的一种冷却方式。为适应我国部分地区水资源不充裕的特点，火力发电中的空冷机组得到了较快的发展。但在其运行过程中，特别是在冬季低温环境下，由于冷凝蒸汽流量太小或者冷却空气流动不均匀等情况，容易造成翅片管束局部结冰的问题，影响凝气设备的正常运行，进而降低机组发电效率和经济效益。
现有技术一般通过测量凝结水联箱内的凝结水温度和抽真空温度来对空冷岛进行监测，或将测温电缆按一定规律布置于翅片管束的内测侧与外侧，根据测量点的温度变化完成对空冷岛的监测，但这些手段均不能完整地反映翅片管束内的温度场情况，也无法有效预防翅片管束局部冻结的问题。
技术实现要素：
针对上述提到的目前空冷岛翅片管束温度场测量的不足，本发明提出一种基于声学技术的空冷岛翅片管束温度场测量装置及其方法。
本发明一方面提出一种基于声学技术的空冷岛翅片管束温度场测量装置，包括：空心管道、声波导管、声波发生器、声波传感器、功率放大器、信号调理器、接线盒、输入/输出设备、工控机，其特征在于：所述声波导管一端安装在所述空心管道上，另一端和所述声波发生器连接；沿所述空心管道间隔一定距离依次布置第1—n个声波传感器m1-n；其中n＝1,2,3,…,n，从而形成多段声波路径；所述声波发生器先后发出固定频率的声波信号，所述声波传感器将接收到的声波信号传递给所述信号调理器，通过所述输入/输出设备传递给所述工控机；所述工控机将接收到的声波信号进行计算和分析，并经过整理得到空冷岛翅片管束温度场信息。
进一步，声波传感器m1-n根据翅片管束布置结构和测量精度要求确定安装的数量，其中n为正整数。
进一步，空心管道表面开有小孔，内介质通过经过小孔的空气流动与翅片管束周围温度相平衡。
本发明另一方面提供一种基于声学技术的空冷岛翅片管束温度场测量方法，包括：
1)以声波传感器不失真地接收到声波信号为前提在空冷岛翅片管束上布置本发明提出的一种基于声学技术的空冷岛翅片管束温度场测量装置对空冷岛翅片管束温度场进行重建；
2)声波发生器产生伪随机信号中的m序列声波信号；
3)计算声波信号在两个声波传感器之间的飞渡时间τ：对于两个声波传感器接收到的声波信号的时延值，采用时延估计算法得到不同时间声波信号形成的互相关函数达到最大峰值所对应的时间延迟，即为声波信号在两个声波传感器之间的飞渡时间τ；
4)计算两个传感器之间气体介质的平均温度t：
相邻两个传感器之间的距离为l，翅片管束内气体介质的气体常数为z。
进一步，采用时延估计算法得到不同时间声波信号形成的互相关函数达到最大峰值所对应的时间延迟为：
假设某两个声波传感器接收到的声波信号的数字模型如下：
式中：x1(n)、x2(n)分别为声波传感器1和2接收到的声波信号函数；s(n)为声源信号函数；w1(n)和w2(n)分别为2个声波传感器接收到的高斯白噪声函数；d为2个声波传感器之间的相对时间延迟；α为声波相对衰减系数；
假设s(n)、w1(n)和w2(n)为互不相关的平稳随机过程，则x1(n)和x2(n)的二次相关函数为：
上述式子可转换为：
rrr(τ)＝αrrs(τ-d)
令α＝1，则rrr(τ)＝rrs(τ-d)为信号s(n)的二次相关函数，由相关函数性质可知，r(τ-d)≤r(0)，当τ＝d时，rrr(τ)取得最大值，计算两声波传感器收到的信号的互相关函数求出峰值出现的时刻，则该时刻就是对应的时间延迟d。
更进一步，通过二次相关phat-β算法来得到时间延迟，具体为通过求相邻两个声波传感器得到的信号的自相关函数和互相关函数之间的互功率谱，并在频域内给予一定的加权，来对声波信号和背景噪声进行白化处理，增强信号中信噪比比较高的频率成分，从而抑制噪声的影响，再将结果通过傅里叶逆变换到时域，得到两信号之间的广义互相关函数，其表达式为：
式中：ψ12(ω)为广义互相关加权函数，
本发明的有益效果为：非接触式测量，可应用于高温、复杂的环境；安装方便简单，不影响原有系统的正常运行，设备成本较低；可以根据翅片管束尺寸和测量精度要求调整传感器数量，测量精度高，灵敏性好；能够及时整体反馈翅片管束温度场变化情况，便于调节和控制。
附图说明
图1为本发明的空冷岛翅片管束温度场测量装置的声波发生器和声波传感器布置图；
图2为本发明的信号处理流程图；
图3为本发明的信号m序列制取原理；
图4为本发明的测量路径平面示意图；
图5为本发明的phat-β算法的计算流程。
具体实施方式
本发明提出一种基于声学技术的空冷岛翅片管束温度场测量装置及其方法，下面根据附图对本发明作进一步详细说明。
如图1所示为本发明的空冷岛翅片管束温度场测量装置的声波发生器和传感器布置图。该空冷岛翅片管束温度场测量装置，包括空心管道4、声波导管2、声波发生器1、声波传感器3、功率放大器、信号调理器、接线盒、输入/输出设备、工控机；其中空心管道4表面开有小孔5，空心管道4内介质通过经过小孔的空气流动与翅片管束周围温度相平衡；声波导管2一端安装在空心管道4上，另一端和声波发生器1连接；沿空心管道4间隔一定距离依次布置第1—n个声波传感器m1-n；其中n＝1,2,3,…,n，n为正整数，从而形成多段声波路径；根据翅片管束布置结构和测量精度要求确定声波传感器m1-n安装的数量。测量精度主要受信号调理器的采样频率以及声波传感器布置间隔距离影响。当采样频率为102400、两声波传感器距离为0.2m时，工控机上显示的温度将以10℃左右阶跃变化，即例如实际温度为20℃时，工控机显示温度为20℃；实际温度为25℃时，显示温度还是20℃；但实际温度升为30℃时，显示温度则才变为30℃。当两声波传感器距离为1m时，工控机显示的温度将以2℃左右阶跃变化，精度上升。同型号的空冷岛其翅片管束长度也不尽相同。例如需要测量长度为5m的管束，我们采用声波传感器间隔1m的布置方式，则安装6个声波传感器。声波发生装置先后发出固定频率的声波信号，声波传感器将接收到的信号传递给信号调理器，通过输入/输出装备传递给工控机；工控机将接收到的信号进行计算和分析，并经过整理得到温度场信息。该装置在测量翅片管束温度场信息时，位于翅片管束上方，靠近翅片管束但不碰触。
如图2所示，声波传感器m1-n的输出与信号调理器连接，信号调理器通过接线盒、输入/输出设备与工控机连接，声波发生器与功率放大器连接，功率放大器通过接线盒、输入/输出设备与工控机连接。工控机发出的声波信号转换成模拟信号后，通过功率放大器放大输出至声波发生器，声波发生器发射出声波信号，声波传感器接收到声波信号后通过信号调理器、接线盒和输入/输出设备将采集到的数据传输到工控机内，工控机利用labview对采集的数据进行分析和处理得到声波飞渡时间，从而计算得到温度场信息。声波发生装置和声波传感器的布置应根据翅片管束结构和以声波传感器能够不失真地接收到声波信号为前提，从而准确地对温度场进行重建。当声波传感器距离声源过远时，由于混响现象，声波传感器采集到的信号会失真而导致测量结果错位；运用a/d转换，数值滤波技术，快速傅里叶变换和阿达玛变换等数值信号处理技术，可以大幅提升整个测量过程的速度。
声波发生器产生的声源信号为伪随机信号中的m序列。如果一个序列，一方面它是可以预先确定的，并且可以重复生产和复制的；另一方面它又具有某种随机序列的随机特性(即统计特性)，这种序列便为伪随机序列。m序列制取原理如图3，当反馈移位寄存器的联级多项式f(x)＝1+c1x+c2x2+c3x3+……+cnxn为本原多项式时，输出即为m序列。m序列是一种周期序列，具有很好的均衡性、自相关性和抗噪能力。由于输入的m序列信号与待测系统的背景噪声是不相关的，便可通过多次测量平均等手段来减少背景噪声的影响。一般来讲，在信噪比低于-20db的环境中仍可以用这种方法来进行有效的测量；而常规的方法要求信噪比不小于10db。利用信号处理的方法将m序列信号的频带变窄。
如图4所示，应用声波飞渡时间时延估计算法计算温度，设相邻两个声波传感器之间的距离为l(m)，空心管道内气体介质的气体常数为z，当声波发生器发出声波信号后，声波经过上述两个声波传感器的飞渡时间为τ(ms)，则两个声波传感器之间气体介质的平均温度t(℃)为：
如图5所示为二次相关phat-β算法流程。先求得原始信号x1(n)的自相关函数
对于两个传感器接收到的信号的时延值，采用时延估计算法得到不同时间信号形成的互相关函数达到最大峰值所对应的时间延迟，即为声波信号在两个传感器之间的飞渡时间τ；
假设某两个声波传感器接收到的声波信号的数字模型如下：
式中：x1(n)、x2(n)分别为声波传感器1和2接收到的声波信号函数；s(n)为声源信号函数；w1(n)和w2(n)分别为2个声波传感器接收到的高斯白噪声函数；d为2个接收器之间的相对时间延迟；α为声波相对衰减系数；
假设s(n)、w1(n)和w2(n)为互不相关的平稳随机过程，则x1(n)和x2(n)的二次相关函数为：
上述式子可转换为：
rrr(τ)＝αrrs(τ-d)
令α＝1，则rrr(τ)＝rrs(τ-d)为信号s(n)的二次相关函数，由相关函数性质可知，r(τ-d)≤r(0)，当τ＝d时，rrr(τ)取得最大值，计算两声波传感器收到的信号的互相关函数求出峰值出现的时刻，则该时刻就是对应的时间延迟d。
对声源信号和背景噪声进行白化处理，通过求相邻两个传感器得到的信号的自相关函数和互相关函数之间的互功率谱，并在频域内给予一定的加权，来对声源信号和背景噪声进行白化处理，增强信号中信噪比比较高的频率成分，从而抑制噪声的影响，再将结果通过傅里叶逆变换到时域，得到两信号之间的广义互相关函数，其表达式为：
式中：ψ12(ω)为广义互相关加权函数，
若取权函数
则称之为二次相关phat-β算法，二次相关phat-β算法时延估计峰值尖锐，具有较强的抗混响与抗噪能力。
通过上述方法能够及时整体反馈翅片管束温度场变化情况，便于调节和控制。