一种基于阵列式静电传感器的速度场层析成像方法
【专利摘要】本发明公开了一种基于阵列式静电传感器的速度场层析成像方法，具体按照以下步骤实施：步骤1：建立阵列式静电传感器的三维仿真模型，获取动态灵敏场S和比例因子矩阵K；步骤2：阵列式静电传感器从管道的不同角度检测管道内带电颗粒流动所产生的静电信号；步骤3：调理电路将静电信号转换为?5V?+5V的交流电压信号，然后经数据采集卡传输到计算机进行处理；步骤4：获取速度频率场P；步骤5：根据速度频率场P和比例因子K得到速度场V。本发明一种基于阵列式静电传感器的速度场层析成像方法，同现有的静电传感器测速方法相比，计算效率更高，速度更快，精度更高。
【专利说明】
-种基于阵列式静电传感器的速度场层析成像方法
技术领域
[0001] 本发明属于气固两相流测量技术领域，具体设及一种基于阵列式静电传感器的速 度场层析成像方法。
【背景技术】
[0002] 气固两相流动过程中固体颗粒速度的大小是气固两相流的重要特征参数，其无论 是实现流动过程的在线监测，还是建立气固两相流动的数学模型都具有重要意义。目前，基 于不同的测量原理，人们已经研究和开发了多种实现气固两相流中固相颗粒速度测量的方 法，如多普勒法，示踪法，相关法和空间滤波法等。多普勒测速WDoppler频移定理作为固体 速度测量的基础，可W使用激光、微波或者超声波作为能源，具有传感器形式多样，测量精 度较高的优点，但是对于恶劣工业环境下的复杂气固两相流动，存在安装复杂，受测量工况 影响较大，且成本较高的缺点;示踪法的测量原理是在两相流中的任一相加入示踪物，根据 从下游示踪物的采样或检测的时间来确定混合物的流速，此方法的原理简单，其一般要求 加入放射性物质，具有不需要校准的优点而且装置安装方便，但即使放射源采用具有较短 半衰期的放射性同位素，也会导致严重的污染问题，且成本较高;相关法利用上下游随机信 号的相似性实现速度测量，具有测量范围宽，适应性强的优点，只需选择合适的传感器便可 W实现对不同对象的测量，但测量时要求流动稳定，固相弥散度尽可能均匀，并且计算量较 大，不适于速度场的获取;空间滤波法作为一种光学测速方法是上个世纪六十年代由Ator 教授提出的，目前已延伸到其他传感器空间滤波效应上，其利用传感器的空间频率响应特 性实现速度参数的测量，具有适用范围广，使用灵活的突出优点，但是其测量精度容易受到 流型，粒径等因素的影响。
[0003] 在气固两相流中，固体颗粒在气体携带下流动，流动过程中由于颗粒与颗粒之间， 颗粒与管道之间的接触摩擦，颗粒与气体之间的相对滑移W都可能使颗粒产生自然荷电现 象，从而使运动的颗粒形成了静电噪声信号。静电噪声信号包含了大量的颗粒流动参数信 息，对其加W适当的信息处理方法，即可获得两相流流动参数(速度、浓度、流量等）。静电传 感器能够有效地捕捉到运些静电噪声信号，其具有结构简单，灵敏度高，适应各种工业环境 等多种优点。
【发明内容】
[0004] 本发明的目的是提供一种基于阵列式静电传感器的速度场层析成像方法，解决了 现有的静电传感器测速精度低的问题。
[0005] 本发明所采用的技术方案是，一种基于阵列式静电传感器的速度场层析成像方 法，具体按照W下步骤实施：
[0006] 步骤1:建立阵列式静电传感器的=维仿真模型，获取动态灵敏场S和比例因子矩 阵K;
[0007] 步骤2:阵列式静电传感器从管道的不同角度检测管道内带电颗粒流动所产生的 静电信号；
[0008] 步骤3:调理电路将静电信号转换为-5V-巧V的交流电压信号，然后经数据采集卡 传输到计算机进行处理；
[0009] 步骤4:获取速度频率场P;
[0010] 步骤5:根据速度频率场巧P比例因子K得至幡度场V。
[0011] 本发明的特点还在于：
[001 ^ 步骤巧体为：
[0013] 动态灵敏场S的获取
[0014] 建立阵列式静电传感器=维仿真模型，将阵列式静电传感器检测区域的截面划分 为M个像素单元，并依次计算阵列式静电传感器各电极对每个像素单元Im/s时的动态灵敏 度Dij，然后对化j进行FFT变化，电极的序号i = 1，2，3-'H，H为电极的个数，像素单元的序号j =1，2，3-M，依次获取第i个电极对第j个像素单元的峰值频率f U，最终建立阵列式静电传 感器的动态灵敏场S:
[001 引
（1)
[0016] 其中，S是HXM矩阵；
[0017] 比例因子K的获取
[0018] =维仿真模型中，设定所有像素单元存在单位电荷，计算当所有像素单元上的电 荷W单位速度Im/s运动时每个电极上的动态灵敏度Di，i = 1，2，3…H，然后对Di进行FFT变 换，从而获取不同电极上的峰值频率化i，用矩阵RJ表示各个电极峰值频率化i的集合，RJ为H Xl矩阵，比例因子K为：
[0019] K = ST 冲 U (2)
[0020] K为M Xl的矩阵麻为矩阵S的转置。
[0021 ]步骤4速度频率场P的获取方法为：
[0022] 计算机对采集到的阵列式静电传感器每个电极上的静电信号进行功率谱分析，根 据每个电极的功率谱分析结果计算每个电极的等效峰值频率Fi:
[0023]
(3) 1234 其中，P表示功率谱分析中的幅值，f表示功率谱分析中的频率，i = l，2,3-'H，H为 电极的个数，n = l，2,3-'N，n表示功率谱分析中的频率的序号，N表示总的频率数； 2 采用LBP算法计算截面的速度频率场P: 3 P = ST 冲 （4) 4 P是MXl的矩阵，M为截面划分的像素单元数，ST为动态灵敏场S的转置，F是Fi的集 合，是HXl矩阵。
[0028] 步骤5中速度场V为：
[0029] V = P/K (5)
[0030] V是M X 1矩阵，每个元素代表一个像素单元的速度。
[0031] 本发明的有益效果是：本发明一种基于阵列式静电传感器的速度场层析成像方 法，将静电传感器的空间滤波特性与层析成像技术相结合，充分利用静电传感器空间滤波 测速的简单，可靠，灵活，低成本，高适用性等优点，结合已有的层析成像技术将测量对象进 行像素单元细分，通过建立不同像素单元的空间滤波灵敏度(动态灵敏场），从而克服现有 静电传感器空间滤波测速受管内颗粒分布影响较大的缺点，同时发挥层析成像技术快速直 观的优点实现整个流动截面速度场参数的在线测量。本发明同现有的静电传感器测速方法 相比，计算效率更高，速度更快，精度更高。
【附图说明】
[0032] 图1是本发明中采用的阵列式静电传感器的俯视图；
[0033] 图2是本发明中基于阵列式静电传感器的气固两相流测量系统图；
[0034] 图3是本发明中信号调理电路的电路图。
[0035] 图中，1.计算机，2.法兰，3.绝缘管道，4.电极，5.接地电极，6.外屏蔽层，7.同轴屏 蔽电缆，8.信号调理电路，9.数据采集卡。
【具体实施方式】
[0036] 下面结合附图和【具体实施方式】对本发明进行详细说明。
[0037] 实施本发明的基于阵列式静电传感器的速度场层析成像方法时的阵列式静电传 感器的形式多样，可W将电极4布置在绝缘管道3外表面实现非接触式测量，也可W将电极4 镶嵌在管道内壁，即通过在金属管道内壁开槽铺垫绝缘层将电极4嵌入的方式，使得电极4 与管道内壁平齐，实现非接触式测量。实现速度场层析成像时均为阵列式静电传感器获取 静电信号后，经信号调理电路8和数据采集卡9进入计算机1，在计算机1内计算得到速度场 V，并通过IW图形的形式直观的显示。
[0038] 本发明中具体实施时采用的阵列式静电传感器的结构如图1、基于阵列式静电传 感器的气固两相流测量系统图如图2所示，包括在圆形绝缘管道3中通过法兰2安装阵列式 矩阵铜制电极4,电极4的两侧均设置有环状接地电极5,电极4的外部设置有外屏蔽层6,将 阵列电极4和接地电极5通过同轴屏蔽电缆7与信号调理电路8连接，信号调理电路8如图3所 示，其中，图3(a)为静电传感器的电压源等效电路图，其中Rt和C分别代表电压放大电路的 等效输入电阻和电容，其中等效电阻Rt由静电传感器的绝缘电阻和放大器的输入电阻两部 分组成；等效电容由静电传感器的对地电容、电缆的等效电容和运算发大器的输入电容组 成；图3(b)是电压跟随放大电路，信号调理电路8依次与数据采集卡9、计算机1连接，阵列式 静电传感器将绝缘管道3内带电颗粒的流动信息捕获后，经同轴屏蔽电缆7传送到信号调理 电路8,经电压跟随和放大后转换为-5-+5V的电压信号送数据采集卡9,经过A/D转换后送入 计算机1，计算机1利用采集到的静电信号完成速度场信息的反演计算。
[0039] 本发明基于阵列式静电传感器的速度场层析成像方法，具体步骤为：
[0040] 步骤I:建立阵列式静电传感器的=维仿真模型，获取动态灵敏场S和比例因子矩 阵K，具体为：
[0041] 动态灵敏场S的获取
[0042] 将阵列式静电传感器检测区域的截面划分为M个像素单元，并依次计算阵列式静 电传感器各电极4对每个像素单元Im/s时的动态灵敏度Du,然后对Du进行FFT变化，电极4 的序号i = l，2,3-H，H为电极4的个数，像素单元的序号j = l，2,3-M，依次获取第i个电极4 对第j个像素单元的峰值频率fu，最终建立阵列式静电传感器的动态灵敏场S:
[0043]
(1)：
[0044] 其中，S是H XM矩阵；
[0045] 比例因子K的获取
[0046] =维仿真模型中，设定所有像素单元存在单位电荷，计算当所有像素单元上的电 荷W单位速度Im/s运动时每个电极上的动态灵敏度Di，i = 1，2，3…H，然后对Di进行FFT变 换，从而获取不同电极上的峰值频率化i，用矩阵RJ表示各个电极峰值频率化i的集合，RJ为H Xl矩阵，比例因子K为：
[0047] K = ST 冲 U (2)
[004引K为M Xl的矩阵麻为矩阵S的转置。
[0049] 步骤2:阵列式静电传感器从管道3的不同角度检测管道3内带电颗粒流动所产生 的静电信号；
[0050] 步骤3:调理电路則尋静电信号转换为-5V-巧V的交流电压信号，然后经数据采集卡 9传输到计算机1进行处理；
[0051 ]步骤4:获取速度频率场P，具体为：
[0052] 计算机1对采集到的阵列式静电传感器每个电极4上的静电信号进行功率谱分析， 根据每个电极4的功率谱分析结果计算每个电极4的等效峰值频率Fi:
[0053]
(3)
[0054] 其中，P表示功率谱分析中的幅值，f表示功率谱分析中的频率，i = l，2,3-'H，H为 电极4的个数，n = l，2,3-'N，n表示功率谱分析中的频率的序号，N表示总的频率数；
[0055] 采用LBP算法计算截面的速度频率场P:
[0056] P = S^F (4)
[0057] P是MXl的矩阵，M为截面划分的像素单元数，ST为动态灵敏场S的转置，F是Fi的集 合，是HXl矩阵。
[005引步骤5:根据速度频率场P和比例因子K得到速度场V:
[0059] V = P/K (5)
[0060] V是M X I矩阵，每个元素代表一个像素单元的速度。
[0061] 本发明一种基于阵列式静电传感器的速度场层析成像方法，将静电传感器的空间 滤波特性与层析成像技术相结合，充分利用静电传感器空间滤波测速的简单，可靠，灵活， 低成本，高适用性等优点，结合已有的层析成像技术将测量对象进行像素单元细分，通过建 立不同像素单元的空间滤波灵敏度(动态灵敏场），从而克服现有静电传感器空间滤波测速 受管内颗粒分布影响较大的缺点，同时发挥层析成像技术快速直观的优点实现整个流动截 面速度场参数的在线测量。本发明同现有的静电传感器测速方法相比，具有计算效率更高， 速度更快，精度更高的优点，为实现实时在线测量速度场参数提供了全新的手段。
【主权项】
1. 一种基于阵列式静电传感器的速度场层析成像方法，其特征在于，具体按照以下步 骤实施： 步骤1:建立阵列式静电传感器的三维仿真模型，获取动态灵敏场S和比例因子矩阵K; 步骤2:阵列式静电传感器从管道(3)的不同角度检测管道(3)内带电颗粒流动所产生 的静电信号； 步骤3:调理电路(8)将静电信号转换为-5V-+5V的交流电压信号，然后经数据采集卡 (9)传输到计算机(1)进行处理； 步骤4:获取速度频率场P; 步骤5:根据速度频率场P和比例因子K得到速度场V。2. 根据权利要求1所述的一种基于阵列式静电传感器的速度场层析成像方法，其特征 在于，所述步骤1具体为： 动态灵敏场S的获取 建立阵列式静电传感器的三维仿真模型，将阵列式静电传感器检测区域的截面划分为 M个像素单元，并依次计算阵列式静电传感器各电极(4)对每个像素单元lm/s时的动态灵敏 度Dij，然后对D ij进行FFT变化，电极⑷的序号i = 1，2，3…H，H为电极⑷的个数，像素单元 的序号」= 1，2,3···Μ，依次获取第i个电极(4)对第j个像素单元的峰值频率fij，最终建立阵 列式静电传感器的动态灵敏场S:(1) 其中，S是HXM矩阵； 比例因子K的获取 三维仿真模型中，设定所有像素单元存在单位电荷，计算当所有像素单元上的电荷以 单位速度lm/s运动时每个电极上的动态灵敏度Di，i = 1，2，3···Η，然后对Di进行FFT变换，从 而获取不同电极上的峰值频率f m，用矩阵FU表示各个电极峰值频率f m的集合，FU为H X 1矩 阵，比例因子K为： K = ST*FU (2) K为M Xl的矩阵;St为矩阵S的转置。3. 根据权利要求1所述的一种基于阵列式静电传感器的速度场层析成像方法，其特征 在于，所述步骤4速度频率场P的获取方法为： 计算机(1)对采集到的阵列式静电传感器每个电极(4)上的静电信号进行功率谱分析， 根据每个电极(4)的功率谱分析结果计算每个电极(4)的等效峰值频率F1: (3) 其中，P表示功率谱分析中的幅值，f表示功率谱分析中的频率，i = I，2，3···Η，H为电极 (4)的个数，η = 1，2，3···Ν，η表示功率谱分析中的频率的序号，N表示总的频率数； 采用LBP算法计算截面的速度频率场P: P = ST*F (4) P是M X 1的矩阵，M为截面划分的像素单元数，St为动态灵敏场S的转置，？是巧的集合，是 HXl矩阵。4.根据权利要求1所述的一种基于阵列式静电传感器的速度场层析成像方法，其特征 在于，所述步骤5中速度场V为： V = P/K (5) V是M X 1矩阵，每个元素代表一个像素单元的速度。
【文档编号】G01P5/08GK106018872SQ201610311039
【公开日】2016年10月12日
【申请日】2016年5月11日
【发明人】高鹤明
【申请人】西安理工大学