李岩++杜永斌++宋海丰++满志强++任相花

摘要：针对在处理大规模样本集的ECT系统数据时，SVM算法存在的图像重建的成像精度不高及速度慢的问题，采用了轮换对称分块支持向量机CSPSVM算法.算法利用ECT系统模型的轮换对称性，将大样本矩阵按照成像单元某一层按轮换对称性进行简化，并选择性分块，形成多个小样本矩阵；然后分别采用SVM算法进行训练，用得出的决策函数进行样本预测；最后将各成像单元组合成像.图像重建实验结果表明使用CSPSVM改进算法要比单独使用SVM算法重建图像具有更高的分类精度和更短的成像时间.

关键词：电容层析成像；支持向量机；轮换对称性；选择分块；图像重建

DOI：10.15938/j.jhust.2015.02.008

中图分类号：TP391.4

文献标志码：A

文章编号：1007-2683（2015）03-0040-05

0 引 言

电容层析成像（electrical capacitance tomo-graphy，ECT）技术是近年来发展起来的一种多相流可视化检测技术，具有非侵入、响应速度快、结构简单、成本低、安全性好及适用广泛等优点，而由于结构本身所固有的，独立电容测量值少、系统敏感场的软场特性及求解问题的非线性等因素的存在，使得ECT系统图像重建算法成为了实际应用中极其关键的环节.

ECT系统图像重建已经采用了多种算法，例如基于迭代的图像重建算法、基于线性反向投影的图像重建算法、基于Tikhonov正则化及TV正则化的图像重建法和基于神经网络的图像重建算法等.但是图像重建算法作为ECT技术应用中的重要环节，随着工业生产的迅速发展，对成像的实时性要求日益提高，输入数据的规模也不断扩大，不管在图像重建速度还是在图像重建质量方面都对电容层析成像算法有了新的要求.

支持向量机（support vector machines，SVM）作为机器学习的一个热点领域，具有良好的泛化能力以及较好的分类精确性，为ECT图像重建提供了有力的手段.

本文就SVM算法在ECT图像重建的规舰模样本问题中训练时间长以及成像质量低的问题，提出了轮换对称分块支持向量机CSPSVM（cyclic sym-metrical partition support vector machine）算法，并使其应用到ECT系统，采用CSPSVM算法，将大样本矩阵按ECT系统结构旋转对称分块，针对某一分块单元，从全局样本矩阵中选出适合这一单元的小样本组合，通过训练可得到最佳模型.

1 基于SVM的ECT图像重建原理

1.1

ECT系统电容测量传感器模型

将ECT系统电容传感器模型中电容极板设置为12极，用ANSYS有限元分析软件建立模型，对传感器设置的电极板编号，分别为1至12，然后进行网格剖分成有限元，采用ADPL语言CMA-TRIX宏命令求解电容值，最终得到66个极间电容计算结果.

ECT系统电容传感器中，极板间电容值是随着敏感场区域内的介电常数（即管道内的介质分布情况）变化而变化.对于12电极的ECT系统，因结构的轮换对称性，只有6种不同类型的敏感场，分别对应于极板对1-2、1-3、1-4、1-5、1-6及1-7的相对敏感场.

1.2基于SVM的图像重建原理

SVM用于分类问题就是要寻找一个最优分类超平面，以此平面作为决策超平面进行分类，同时它还通过引入核函数，而巧妙地解决了将低维空间向量映射成高维空间向量时所带来的维灾难问题.

ECT系统求解的样本数据矩阵是由66个电容值和66个敏感度值共132个属性值构成的一个非线性样本集.通过非线性映射，把低维样本空间向量映射到高维特征空间，然后，在高维特征空间中采用线性支持向量机算法，构造出一个最优分类超平面，解决非线性分类问题.在保证训练样本分类误差极小化的前提下，使分隔特征空间两类属性点的两个超平面之间距离最大，保证所求解收敛于全局最优，

将ECT电容传感器模型剖分为n个单元，可获得m个样本，将其分为训练样本集和预测样本集.SVM通过训练样本集进行学习后，便可对预测样本集进行分类.

假设训练样本集 ，其中 为训练样本集 表示第i个剖分单元为高或低介电常数介质，i=l，2，…，N，N为训练集样本总数.

为了构造出最优超平面，引入非负变量 ，则非线性不可分SVM原问题为：

设核函数 ，文中采用高斯径向基核函数，则得到决策函数表达式为

2 轮换对称分块SVM图像重建算法

2.1建立模型

实验中将ECT电容传感器模型管道内截面剖分成192个成像单元，首先将管道内截面以极板中心和隔板中心为剖分线分成24个全等扇形，再以管道内截面中心为圆心作7个同心圆，半径分别为 ，其中R为管道内截面半径，于是将每个扇形都剖分成相等的8份，这样管道内截面就被剖分成192份等面积成像单元（见图1）.

ECT电容传感器模型是以其管道截面圆心为中心对称的，12电极均匀布置在管道外侧，且设置相同的参数，即12电极参数完全相同，那么每一极都可以进行互换，是具有轮换对称性的.

按逆时针顺序给12电极编号，分别为1至12；按逆时针顺序对应极板给各成像单元编号，扇形分块顺序l至24，同心圆分块顺序外向内为1至8，若把同心圆分块称为层，层数设为X，则X=1，2，…，8.那么对应于每个成像单元编号分别为X01至X24.而由于管道按照12电极的结构对称，扇形分块单元中，奇数编号1，3，5，…，23的单元是轮换对称的，偶数编号2，4，6，…，24的单元也是轮换对称的，而且，奇数编号与偶数编号的单元也分别为镜像对称的，应用传感器模型的这些特性，可以使样本数据进行充分的简化.