封闭管路内散装固体物料流速的相关测量.docx

1、0钟国徘"七流疽，拿河童I。w封闭管路内散装固体物料'.流速的相关测量*中国矿业大学机械系黄民肖世德/魏任之本澡题为国家赦委博士点基金课题、(91年6月通过鉴定.1. 引言长期以来，两相流(气/液、气/固、液/固)的流动参数的测量和控制是许多生产过程中一个急需解决而又未能很好解决的问题.对连续单相气体和液体介质的流动用经典的孔板差压等方法来监测，可以达到足够的精度。而对于两相流和多相流尤其固体介质的流动用各种传统方法来达到实时监测是难以奏效的。而在很多工业生产过程中，为了实现生产过程的连续化、自动化和高速化，迫切需要能够实时在线连续监测和控制两相流的流动情况.本文利用相关测量

2、技术研制了一种用于管道内气固两相流流速在线监测的实时相关测量系统。2. 相关测速的基本原理相关测速系统原理如图Lup%-传哆!一ir.I*IMJUIh.:b卜.:ffjfi物鞘图1互相关测速系统原理图当被测流体在管道内作稳定流动时，随机流动噪声信号x(t)y(t)可以认为是符合各态历经平稳随机过程x(t)和y(t)的两个样本函数.只要两个传感器的间距L合理，两个传感器和变送器的静动态性能完全一样，可以认为两个随机流动噪声信号x(t)和y(t)是相似的，其互相关函数峰值对应的滞后时间t,就是被测流体从截面aa'和bb'的流动时间，在理想流动条件下，被测流体的体积平均流速Vcp可以

3、用相关速度Vc来表示，即：Vcp=Vc=L/l从而，两相流的流速问题转化为随机流动噪声信号的检测和互相关峰值位置的确定问题，因此一般相关测量系统可分为以下二个部分：(1) 传感器和变送器系统，用于获取随机流动噪声信号，并转化为可供后续处理的电信号，其信号检取的正确与否是相关测量系统的基础。(2相关处理系统.用于确定相关函数及其最大峰值位置.相关处理过程的在线实时性和正确可靠性是相关测量系统实用化的前提c同时应建立一个实验标定系统，用于测定由(1)、(2)两部分构成的相关测量系统的性能和精度以及修正常数.实验标定系统是整个相关测量系统的精度保证。3. 传感器的设计与制作3.1 传感器类型的选择迄

4、今为止，国内外已研制出不同类型的流动噪声传感器，有超声式，光学式，电动力学式，核放射式，热线探针式，电容式等等.对于气固两相流而言，由于超声式存在衰减这个难以克服的障碍，核辐射式存在安全性和响应过低问题，光学式存在光污染和透光性问题等等，使电容式传感器成为唯一的最佳选择，它具有以下优点：(1) 电容式传感器原理简单、可靠；(2) 可实现非接抵式测量，便于“卡钳式”安装；(3) 适用性广，响应快；(4) 价格低廉，便于应用；3.2 电容传翌器的结构设计和制造传感器的制作：(试验管道为。100mm内径的拿叙乙烯绝缘塑料管道).用两片宽b=10mm的薄铜皮贴于管道的外壁，形成一对电容传感器电极板.相

5、距L=64mm安装两对同样的极板，作为上、下游传感器.如图2所示.92传结构图用细密铁丝网对上下游传感器作超宽屏蔽并接于良性实地，传感器极扳到变送器采用优质屏蔽信号线连接传输.4. 变送器(电容-电压转换器)的设计变送器(电容-电压转换器)的功能是电容传感器的电容信号转换成电压信号,由于电容变化量较小.因而要求变送电路具有较高的灵敏度，同时能对变换得到的技弱电压信号加以放大，达到二次仪表所要求的输入辐值.鉴于集成芯片具有灵敏、稳定、高度集成等优点，本系统的变送电路主要由集成芯片枸成，采用高灵敏度电容测量电路(即调制脉宽电路)作为第一级转换，后经真有效值转换电路，放大电路，滤波电路等四级信号处理

6、，最后得到所需要的理想的电压信号,如图3.声味冲AT号必丽丽有效值电塑明溢号I土四丽矿宅压僖号-压信号S3变送原理椎图整个变送器的电路原理图如图4.图4变送电路图所研制的传感器及其变送系统的抗干扰措施主要体现在以下几个方面：(1) 采用脉宽调制式电容检测电路和真有效值转换电路-均化电容具有抗干扰作用.(2) 传感器采用铁丝网超宽屏蔽，消除环境电容的影响.(3) 变送器与传感器本体尽量接近，减少介入的导线电容。(4) 传感器和变送器及屏蔽系统接于良好实地。(5) 采用单一稳压电源供电，并消除电源的级间交流藕合效应.(6) 切比雪夫保通混波器的使用。(7) 采用集成元件、减少分立元件。(8) 级间

7、采用隔直电容。5. 相关处理系统5.1相关算法的原理从理论上讲，信号x(t)和y(t)的互相关函数表述了信号在时延域上相似性，其计算公式如下：Rxy(t)=l;Tjx(t)xy(t+t)dt对应离散算式如下：Rxy(m)=l/N£x(k)xy(k+m)m=0,1,、，N-1式中x(kxy(k)为信号x(ny(t)采样值，采样频率为f«可以看出，计算一次相关函数要做N次乘法和N次加法，我们的目标是确定相关函数最大峰值所对应的时间T,计算量正比于N。一方面为了保证相关函数的计算精度，要求采样点数N足够多，以满足Toc条件；另一方面为了保证时延测量的精度，又要求采样频率f足够高，

8、可以想象计算量是巨大的。要提高相关处理的速度，可从两个方面入手：(1)提高互相关函数Rxy(m)的计算速度.(2) 提高寻求峰值时滞t的速度.把信号x(ty(t)加以Ibit极化、化乘法运算为异或非逻辑运算,是对Rxy(m)计算的重大改进，理论上已经证明，随机流动噪声是高斯型的，而高斯型随机过程及其极性互相关函数Rxyp(m)与互相关函数Rxy(m)有如下关系：Rxy(m)2、Rxyp(m)=ancsinn<Rx(0)xRy(0)可见，Rxyp(m)与Rxy(m)的峰值位置是一致的，而极化带来的误差可以通过增加计算点数N来弥补。以上极性互相关函数是各种硬件相关器和软件相关器的理论基础。5

9、、2相对过竿快速寻峰极性相关算法由于电容传感器空间滤波效应的影响，所测得的随机流动噪声信号的频带一般在800Hz左右，可实际工作中又常常要求时延精度在ms或俸级，对这种低信号频带高分辨率要求的互相关处理而言，采用Henry过零时刻算法是最有效的，其计算量和存储量仅与信号过年频率有关.Henry算法中采用绝对过年时刻效据，存在时钟溢出中断计数困难，实际上直接采用相对过竿时刻数据，计算互相关函数是可以的。(所谓相对过竿指下一个过年时刻的度量，不是以时钟起点为基准，而是以上一个过零时刻为基准)此时计算极性互相关函数Rxyp(z)转化为统计过竿时刻信号极性一致和不一致情况下过年时刻累计和之差。同Hen

10、ry绝对过年时刻算法相比：(1) 采用相对过早时刻算法，时延修正T时，不必每项X(t)过竿数据都减去T,节省了计算时间。(2) 本算法只累加极性一致时的过零时刻，节省了计算时间。(3) 以过竿时刻个数做为互相关函数计算时循环控制参数，不以累加总时钟长度为循环终止条件，有自适应功能。(4) 采用相对过零时刻易于实现环形存储和动态跟踪.时钟一经起动，不用再打扰.寻找互相关函数最大峰值对应的时间滞后T实际上是一个一元函数求极值问题.Henry采用第一询扫描求峰，第二遍再用对分法细化求峰的方法，对256个延时只须计算37次(扫描步长为8)就可求出峰值，速度提高了7倍.假设互相关函数Rxy(t)为单峰函

11、数，根据最优化理论，采用0.618法每次收缩区间只须计算一个新函数值而对分法则要计算两个，因此可以把Henry算法寻峰速度再提高一步，这个优点在高分辨率时更加明显.当然真实的极性互相关函数不可能是单峰函数，经常的情况是存在一个最大峰和很多辅助峰(非最大峰).对软件相关而言，这个问题可以这样解决：(1) 根据经验给定寻峰范围Cas,ae，减少不必要的寻峰搜索.(2) 采用扫描法先求出一个粗峰，然后在细化求峰时再用0.618法，扫描步长可以调整,尽量取得长一些.(3) 考虑到稳定流动时，流速不可能突然发生很大变化，因此还可以进行自适应寻峰，即第一轮时用扫描寻峰，此后在前一轮寻峰位置前后自适应搜峰，

12、从而大大减少寻峰工作量.依照以上观点，提出了扫描自适应0.618寻峰算法.根据以上算法编制的相关软件，仿真结果和实际验证表明：软件相关在实时性方面达到了Is,与硬件相关器不相上下，并且在灵活性、适用性方面大大优于硬件相关器.5.3 实时软件相关器的硬件与软件设计研制软件相关器采用8031单片机作为软件相关器硬件基础，整个实时相关软件采用汇编语言开发，软件设计遵循模块化结构程序设计原理，可以分五个层次,.软件模块调用如图5.其中：中模块和眼务程序库模块,用执行.可供软件随时随炮调图5程序堵构框图5.4 极化过零电路的设计与中断处理从电容传感器和变送器输出的随机流动噪声信号是模拟信号(其直流分量已

13、滤除)，为了配合8031单片机完成在线实时极性互相关，需要用两路极化过零电路、对上下游传感器输出信号x(以y(t)进行极化.成为符合逻辑电平的矩形脉冲串，再在矩形脉冲串上下沿处(过零时刻)产生过零中断脉冲.极化过零电路包括：低通滤波、限幅、极化、整形、单稳、组合等部分，其电路原理图如图6.图6极化过零电路原理图8031单片机有两个外部中断INTO,INTI.把两路过零中断脉冲接到其输入端，中断触发方式采用后沿触发.由于过零诊断的发生是随机的，因此设定其优先级为高.中断服务总是花费一定的时间.在一路过零中断正处理时又来了新的过零中断，若中断响应不及时，可能造成过零中断算法失败。Henry在其设计

14、的过零中断电路中采用双稳来保证正确性.本研究设计中用软件来保证：把极化后信号也接到8031的P1接口两个位上。每次每路过零中断时极性相比，若极性相反，本次过零信号有效.否则若极性相同，证明有中断丢失和重迭发生,不执行过零中断处理.若CTC采用单片机本身机器时钟，计时周期为2口,溢出周期为0.131s,有可能相对过零时刻超过0.131s.为了防止出现这一情况，决定采用外部时钟.因为分辨率到2俸根本没有必要，实用时达到1ms和0.1ms足够了.为此设计了晶振分频电路作为外部时钟.为了便于相关处理系统调试和检查，还设计了专门的伪随机信号发生器。6. 动态试验及系统标定6.1动态试验装置的设计与制作为

15、了对所设计的整套传感器、变送器及软件测试系统进行模拟实际的动态试验，特制作一模拟试验架，管道是内径为100mm聚氯乙烯塑料管，长度为L2m,两传感器安装位置以a、b、0表示，a=587mm,b=64mm,0=0.试验物料不失一般性，采用0.5立方毫米左右的碎白石，如图7所示.6.2动态测试系统动态测试系统原理图如图8所示。直流稳压电源土15V作为变送器的直流电源,磁带机记录来自变送器的两路信号，同时用双踪示波器进行观察.调整两路信号的幅值大小等，然后将磁带机记录的两路信号转放给HP3562A进行信号分析，并做软记录和绘图输出.fTEACMR3磴蒂记承仅HP3562A动态信号5仪软盘记承HP74

16、75A一!HP7475A一!图8动恋测试系统原理图6.3系统的标定及实验结果'为了考核电容传感器的性能(信噪比和频带)，论证相关测速系统的可行性和精度，必须对相关测速系统进行标定实验.(1)自由落体实验T=E+L)Vg/2HVg图7动态试验装置用HP3562A跟踪显示自由落体的信号曲线如图9.S9自由落体信号曲线把物料下落处理成自由落体，即用一个徐码或一团物料作自由落体，估算物料经历两个传感器中线需要的时间差匚设两个传感器间距L,上游传感器距物料下落点为H,忽略空气阻力的影响.则有：理想延时估算：H=506mmL=64mmt19.71(ms)HP3562A相关延时：r=19.53(ms

17、)相对误差：0.96%结果表明用HP3562A进行测试有足够的精度.(2)用HP3562A作为参照标定用磁带机录下真实物料流动噪声信号，分别用HP3562A和梢关处理系统进行寻峰，由于HP3562A相关用真实的模拟信号采样数据作的，其缩果可以认为是真实的，而相关处理系统用极化信号过零相关的，两个结果的一致性可以证实相关测速的可行性和极化相关的误差.图10为随机流动噪声的相关曲线。(50次平均结果)ffi10随机流动值号龄相关曲线HP3562A相关时滞：Tavg=18.359ms衰1实时相关软件寻峥跟踪的结果序导12345678Lt_10H12n】4T*竺竺2S227J2S9塑289292230

18、72873001792M平均t*-M4J8S9对比：时钟频率f=15.619kHzHP3562A50次平均：t=18.359ms实时相关结果：t=18.233ms相对误差：0.7%可见：实时相关的精度在1%左右.7. 总结相关测量技术作为一种有广阔应用前景的技术，在本研究工作中为实现相关测量的在线实时性做了一定工作，理论和实验均证明：(1) 采用电容传感器和脉宽调制真有效值变送电路检测微弱的随机流动噪声信号是可行的.对离散物料(气固两相流)的流动噪声检测是较好的，研制成功的电容传感器和变送器是正确实用的，信噪比达40dB左右.(2) 本研究工作中提出的相对过零极性互相关算法和扫描自适应0.618快速寻峰算法是可行的，依此研制的在线实时全软件相关器是成功的，实时性达Is左右，(3) 采用HP3562A的相关结果作为实时相关系统标定依据的办法是实用的.实验结果表明相关测速系统精度在1%左右，(4) 本研究工作证明了采用相关法对离散物料(气/固两相流)进行流速检测是可行的.8. 希考文献1. 相关淹量测景技术，徐苓安，天津大学出版社.2. 日本电子电咤精选，常玉燕译，电子工业出版社，3. 555定时器原理及实用电路集锦.杨兆远，天津大学出版社.4. AnImprovodAlgorithmAllowingFastO