（测试计量技术及仪器专业论文）电容层析成像系统三维数据场可视化及计算机仿真.pdf

电容层析成像系统三维数据场可视化及计算机仿真摘要电容层析成像技术是一种多相流参数检测新技术，三维数据场可视化是运用计算机图形学和图像处理技术，将科学计算过程及计算结果转换为图形及图像在屏幕上显示出来并进行交互处理的理论、方法和技术。本文针对石油工业广泛应用的油水液液两相流系统，利用电容传感器检测到的时间离散的一维投影数据，以电容传感器空间内不同点似，y ) 处的灵敏度分布函数为先验数据，采用m a r c h i n gc u b e s 算法，将液液两相流系统的流动、变化过程再现在屏幕上。本文首先分析了电容层析成像系统传感器原理、系统工作原理及以有限元模型为基础的场域剖分、敏感场分析与计算。在此基础上，利用三维数据场可视化原理，采用三维空间数据场构造等值面方法，将m c 算法中的体元用灰度值矩阵中的一个元素替代；各体元之间等值面图像用一系列相互连接的三角面片替代，运用最大体积法，将相邻两横断面的轮廓线连接起来；利用人眼视两膜停留时间原理，实现电容层析成像系统三维数据场的可视化。本文将以上方法用于圆形管道内油水两相流的实例分析。完成了油水两相流系统的灰度图显示、三维空间数据显示以及三维空间动态显示设计。通过仿真试验，达到了满意效果。通过对工业现场的油水两相流的可视化仿真，将测量数据通过可视化变成图像，从而将工业现场的“黑箱”系统变成透明系统。从海量数据中，找到其中隐藏的规律，为工业测量和控制决策提供有效、直观的数据。算法关键词电容层析成像；两相流：仿真；三维数据场可视化；m a r c h i n gc u b e st h e3 一dd a t af i e l dv i s u a l i z a t i o na n dc o m p u t e rs i m u l a t i o no fe l e c t r i c a lc a p a c i t a n c et o m o g r a p h ys y s t e ma b s t r a c te l e c t r i c a lc a p a c i t a n c et o m o g r a p h y ( e c t ) i sak i n do fn e wt e c h n o l o g yo fm e a s u r i n gp a r a m e t e r so fm b l t i c o m p o n e n tf l o w t h r e e d i m e n s i o n a ld a t af i e l dv i s u a l i z a t i o ni st h et h e o r y ，m e t h o da n dt e c h n o l o g ya b o u th o wt oc o n v e r ts c i e n t i f i cc o m p u t a t i o np r o c e s s e sa n dr e s u l t si n t og r a p ha n di m a g eu s i n gc o m p u t e rg r a p h i c sa n di m a g ep r o c e s s i n gt e c h n o l o g y ，a n dt h e nd i s p l a yt h e mo ns c r e e ns oa st os u p p o r tt h ei n t e r a c t i v e l yp r o c e s s i n g i nt h i sp a p e r ，w ea i m e da tt h ew i d e l yu s e do i lw a t e rl i q u i d l i q u i d2 - p h a s ef l o ws y s t e mi nt h eo i li n d u s t r y m a r c h i n gc u b e sa l g o r i t h mi su s e dt od i s p l a yt h ep r o c e s so ff l o w i n ga n dv a r i a t i o no ns c r e e nu s i n g1 - d i m e n s i o n a lp r o j e c t i o nd a t ad e t e c t e db ye l e c t r i c a lc a p a c i t a n c es e n s o r si nd i s c r e t et i m ea n dt h es e n s i t i 啊t yd i s t r i b t i t i o no fd i f f e r e n tp o i n t ( x ，y ) i nt h es e n s o ri su s e da sp r i o rd a t a i nt h ep a d e l - w ef i r s ta n a l y z e dt h ep r i n c i p l eo ft h es e n s o r si ne c ts y s t e m ，t h ew o r k i n gp r i n c i p l eo ft h es y s t e ma n dt h ef i e l dm e s h i n gb a s e do nf i n i t ee l e m e n tm o d e la sw e l la st h ec a l c u l a t i o no f s e n s i t i v i t yf i l e d o nt h eb a s i so f 也ea n a l y s i s w er e p l a c ee v e r ys o l i de l e m e n to ft h em ca l g o f i t h mb yo n ee l e m e n to fg r a yl e v e lm a t r i x ；t h ei s o v a l u es u r f a c ei m a g eo ft h es o l i de l e m e n t sc a l lb ec o n s t r u c t e db yas e r i e so fc o n n e c t e dt r i a n g u l a rf a c e t s ；t h ec o n t o u rl i n e so ft h ea d j a c e n tt w oc r o s ss e c t i o ni m a g ea r ec o n n e c t e db yu s i n gt h em a x i m u mv o l u m em e t h o d ，a n dt h e nw eu s et h ep r i n c i l ) l eo fr e t i n at i m ed e l a yt oi m p l e m e n tt h ev i s u a l i z a t i o no ft h r e e d i m e n s i o n a ld a t af i e l df o re c ts y s t e m w eu s e 也em e t h o d sm e n t i o n e da b o v ei nt h ea n a l y s i so fo i l w a t e rt w o - p h a s ef l o wi nc i r c u l a rp i p l i n e s a t i s l y i n gr e s u l ti so b t a i n e di nc o m p u t e rs i m u l a t i o ne x p e r i m e n t b yt h ev i s u a l i z a t i o no ft h eo i l w a t e r2 - p h a s ef l o wi n 也ei n d u s t r i a ls i t e ，w ec a l lc o n v e r tt h em e a s u r e dd a t ai n t oi m a g e a n dc h a n g et h eb l a c k b o xs y s t e mi n t ot r a n s p a r e n ts y s t e m w ec a nf i n dt h eu n d e r l y i n gr e g u l a r i t yf r o mm a g n i t u d ed a t as oa st op r o v i d ee f f e c t i v ea n dv i s u a ld a t af o ri n d n s t r i a lm e a s u r ea n dd e c i s i o no f c o n t r 0 1 k e yw o r d se l e c t r i c a lc a p a c i t a n c et o m o g r a p h y ，t w o - p h a s ef l o w ，c o m p u t e rs i m u l a t i o n ，3 - dd a t af i e l dv i s u a l i z a t i o n ，m a r c h i n gc u b e sa l g o r i t h mi i哈尔滨理1 大学工学硕士学位论文第1 章绪论在现代工业中，多相流检测系统得到广泛的应用，发展多相流检测技术是现代工业发展的迫切要求，该技术在解决多相流检测问题上有巨大的发展潜力和广阔的工业应用前景；而科学计算可视化( v i s u a l i z a t i o ni ns c i e n t i f i cc o m p u t i n g )是发达国家在2 0 世纪8 0 年代后期提出并发展起来的一个新的研究领域，它利用计算机图形学和图像处理技术，将科学计算过程中及计算结果转换为图形及图像在屏幕上显示出来并进行交互处理。将多相流检测系统与科学计算可视化结合在一起，可大大提高检测与科学计算的质量，使科学研究的面貌发生根本的变化。1 1 引言“相”是指具有相同成分和相同物理、化学性质的均匀物质部分，或者说相是物质的单一状态，如固态、液态和气态等，在两相流动的研究中往往将其称为固相、液相和气相。一般来说，各相需有明显可分的界面。单相物质的流动称之为单相流，两相流指的是两种相同时存在的流动，多相流是指同时存在多种不同的物质流动i ”。在多相流的研究中，常常将物质分为连续介质和离散( 非连续) 介质。气体和液体属于连续介质，因此又被称为连续相或流体相；固相颗粒、液滴、气泡等属于离散介质，常被称为离散相或颗粒相。连续相和离散相组成的流动叫做两相流动或两相流。由于离散相可以是不同物态、不同化学组成、不同形状的颗粒，如气泡、液滴等，因此可以将复杂的流动概括为两相流动，大大简化了数学模型1 3 。工业中常见的多相流动系统大体上可分为以下五类，其中以两相流最为普遍，尤其是气固两相流、气液两相流和液液两相流【4 】。( 1 ) 气液两相流。热力设备( 如锅炉) 与制冷设备的工作过程，石油、天然气的传送过程等，都存在气液两相流。( 2 ) 气固两相流。炼油中的催化裂化、冶金工业的粉粒体的风力输送，固体火箭发动机的燃气流动等，都属于气固两相流。( 3 ) 液固两相流。采矿、冶金中的矿石选矿，水利工程中的泥沙运动和治理，化工中的液固硫化床等，都包含液固两相流。( 4 ) 液液两相流。冶金、化工中的萃取过程，石油工业中油水系统的油水哈尔滨理工大学工学硕士学位论文分离、洗涤和排放等，都与液液两相流密切相关。( 5 ) 气液液和气液固多相流。在原油开采过程中，油品加氢和精制中的滴流床、淤浆反应器以及化学合成和生化反应的悬浮床工作过程中，气液液和气液固多相流的例子大量存在。从这五大类可以看出多相流系统具有广泛的应用，本文所研究的对象是液液两相流，即石油工业中油水系统的液液两相流。但是多相流检测技术发展现状和水平与飞速发展的现代工业应用之间的矛盾也同益突出，因此，发展多相流检测技术是现代工业发展的迫切要求，研究多相流技术对国民经济的发展具有重要意义。1 2 多相流检测技术现状与发展趋势多相流检测是指采用合适的方法和手段检测多相流的各项参数，其主要参数有：流型( 又称流态，即流体流动的形式和结构) 、分相含率、总流量、分相流量、密度、速度等，典型的流型有层流、环流、多滴流、核心流等，如图1 1所示。其中，流型、分相含率( 对两相流而言为离散像浓度) 和流量是三个最重要同时也是最难以测量的参数，因此，这三个参数的测量也是目前多项流检测技术研究的重点”】。图1 - 1四种典型的流型f i g1 - 14t y p i c a lf l o wp a t t e r n s由于多相流系统各相之间存在着界面效应和相对速度，多相流体流动特性复杂多变，因此，其参数的检测存在很大的难度。虽然国内外诸多的研究人员为此做了大量工作，研究工作不仅涉及各种新技术，如放射线技术、激光多普勒技术、核磁共振技术、超声波技术、光纤技术、流动层析成像技术等，而且也有很多是采用传统的单相流仪表和多相流检测模型进行多参数组合辨识来测哈尔滨理工大学工学硕士学位论文t量的，但这些测量方法和手段基本处于实验室探索阶段，实用化的工业型测量仪表和方法尚未见有关报导。目前，关于多相流检测技，杞的研究发展趋势主要有以下几个方面( 6 ：( 1 ) 针对多相流动过程是一多变量的随机过程，随着随机过程理论和信号处理技术的不断完善和发展，应用数理统计、参数估计和系统辩识的理论和技术，进行多相流参数的估计将成为重要的发展趋势。( 2 ) 随着计算机技术和图像处理技术的发展，获取关于多相流体的多个信息，应用流动层析成像技术，对多相流体局部空间区域进行微观和瞬态的测量将是一个重要的发展方向。( 3 ) 借助电子技术、半导体技术的发展，研制商灵敏度、准确度和可靠性的多相流参数检测仪表是一项重要的基础工作。( 4 ) 在速度测量方面，相关法、空间滤波法和激光多普勒法等测量技术和方法将受到广泛的应用。( 5 ) 如何将成熟的传统单相流参数检测技术和测量仪表用于多相流参数检测当中，这将是个重要的酽究方面。在多相流检测技术发展的几个主要趋势中，流动层析成像技术的研究是一个引人注目的发展方向。本论文就是利用电容流动层析成像技术，通过对两相流体局部空间区域进行瞬态的数据测量，经科学计算，利用计算机图形学和图像处理技术，再现流动的两相流体的流型及离散相浓度。1 3 流动层析成像技术1 3 1 简述层析成像( t o m o g r a p h y ) 技术也称计算机断面( 层析) 成像( c o m p u t e r i s e dt o m o g r a p h y ) ，简称c t 技术或计算机辅助层析成像( c o m p u t e ra s s i s t e dt o m o g r a p h y ，简记为c a t ) ，是指在不损伤研究对象内部结构的条件下，利用某种探测源，根据从被探测对象的外部获得的投影数据，运用一定的数学模型和重建技术，使用计算机生成对象内部的二维图像，重现对象内部特征。层析成像在图像处理学上称为图像重建，它不同于图像处理学中由图像到图像的处理方法，而是一种由投影数据( 在医学工程中投影数据由计算机断层扫描装置提供) 获得图像的特殊图像处理方法，即根据物理对象某一截面的一组投影数据，采用图像重建算法经过计算机处理获得该物体截面的图像，从而直观地反映截面区域的有关信息m 。哈尔滨理工大学工学硕士学位论文8 0 年代以来，c t 技术工业应用的发展引人注目。工业c t 不但成为工业设备或部件的无损检测新军，而且它们正在迈入工业生产过程的在线实州质量监控和工业漫备在役安全检查等领域。鉴于工业c t 具有其它检测手段无法比拟的优点，科研人员将c t 技术“移植”到多相流检测领域，称为过程层析成像，或称为流动层析成像，简称流动成像，以区别于医学工程中的层析成像技术。流动层析成像技术同其他多相流参数检测技术相比，比较明显的特点和优势主要表现在以下几方面m ：( 1 ) 能提供多相流混合体流动过程中局部区域瞬态微观的信息，即管道某一截面处的多相流断层图像，基于这些信息可对多相流系统的流型进行辨别；( 2 ) 与其它测量技术相配合可应用于多相流流量及流速等的实时检测：( 3 ) 通过对图像的处理和分析，可以获得多相流混合体各分相局部浓度的分布，进一一步处理可得各分相的总浓度即离散相浓度；( 4 ) 流动层析成像技术具有适用范围广的优点，各种工业生产过程中常见的多相流均能应用该技术。1 3 2 国内外相关领域研究状况流动层析成像技术的研究开始于2 0 世纪7 0 年代后期，在医学c t 技术所获巨大成功的激励下，研究人员开始尝试将c t 技术应用于多相流参数检测领域，以获得更加直观和更加准确的测量结果。从2 0 世纪8 0 年代中期开始，以英国的曼彻斯特大学理工学院( u n i v e r s i t yo fm a n c h e s t e ri n s t i t u t eo fs c i e n c ea n dt e c h n o l o g y ) m s b e c k 教授为首的研究小组开始了基于电容、超声波等适合于工业应用传感机理的层析成像技术研究，并正式提出了流动层析成像这一概念，以区别于医学的c t 。由于先期的研究成果极大地鼓舞了科技工作者，不同学科领域的研究人员开始联合起来，积极探讨基于不同传感机理的层析成像系统，并将这一多学科交叉的高新技术概括地称为“流动层析成像”技术 9 1 。进入2 0世纪9 0 年代，随着流动层析成像技术的研究和发展，令人鼓舞的研究成果不断涌现，该技术在工业应用中的巨大潜力逐步为人们所认识，这方面的研究也得到了有关政府的重视和资助，同时，流动层析成像技术也得到了工业界的关注，一些工业企业也丌始资助和积极参与该技术方面的研究。在我国，清华大学和天津大学于2 0 世纪8 0 年代率先进行了流动层析成像技术的研究，而后浙江大学、东北大学、武汉大学、浙江工业大学、中国科学院和哈尔滨理工大学等高等院校和科研单位也相继开展研究，已经开发出几种基于不同传感机理的流动层析成像样机系统。这些科研成果得到国家自然科学基金委员会及有关工业部门给予的大力的资金资助，以支持流动层析成像技术的研究和开发。由于多相流的复杂性，难以完全用理论解决问题，欲将理论计算结果应用于工程实践，必须依靠试验来验证或修正，而多相流参数检测技术又是实验研究的重点之一。国内外众多的研究人员为此做了大量的研究工作，已经取得了较大的进展。通常，可将多栩流参数检测技术分为以下几类：( 1 ) 将传统的单相流仪：莨，采用实验校正方法应用于多相流参数测量。( 2 ) 以现有检测手段为项目，采用新的信号处理技术进行多相流的参数测试。将模糊数学、神经网络、小波变换、过程辨识、近代谱估计等理论引入多相流参数测试领域中，来解决复杂的、不确定的，且难以用数学模型精确描述的多相流系统参数测量。例如，应用参数估计测量多相流流速，可克服用传统渡越时间相关法测速存在的随机误差大、分辨率低、动态响应慢等不足，并已研制成仪表雏形。采用电容传感器阵列交流抗杂散微小电容测量电路，基于8 0 5 1单片机及嵌入式等相关装置，在测量气固两相流浓度和质量流率已取得初步的研究结果。( 3 ) 采用近代新技术应用于两相流参数测量。将辐射线技术，激光技术、核磁共振技术、微波技术、光谱技术、新型示踪技术、相关技术、过程层析成像技术等应用于两相流参数测量。本世纪8 0 年代中期兴起的过程层析成像技术( p r o c e s st o m o g r a p h y ，简称p t ) ，主要包括核子层析成像、光学层析成像、电学( 电阻、电容、电磁) 层析成像、微波层析成像、超声层析成像等，能在线连续提供多相流体截面状况的二维可视化信息。过程层析成像技术在多相流测试中存在巨大的潜力，将发挥越来越大的作用。尽管上述这些方法在解决多相流参数检测方面取得很大进展，但是面对飞速发展的现代工业应用的需求，多相参数检测技术还远未满足工业应用和发展的要求，在国内外都属于探索研究的领领域，世界各工业发达国家除形成许多多相流研究中心外，美国中西部2 0 所大学与阿贡国家实验室及工矿企业联合成立了多相流研究所，开展多相流方面的研究工作，涉及多相流问题的国际学术会议平均每年t 多次，刊登有关多相流学术论文的国际杂志已超过3 0 0 种。我国从1 9 8 6 年起已连续四次召开了全国多相流测试学术会议，并于1 9 9 2 年正式成立了中国计量测试学会多相流测试专业委员会。国家在八五、九五科研规划中对两相流的研究给予了高度重视，并在西交大建立了动力工程多相流国家重点实验室，我国多相流测试技术研究从分散到集中并逐步缩小与世界各发达国家的差距。哈尔浜理工大学工学硕士学位论文1 4 三维数据场可视化概述1 4 1 引言科学计算可视化( v i s u a l i z a t i o ni ns c i e n t i f i cc o m p u t i n g ) 是运用计算机图形学和图像处理技术，将科学一算过程中及计算结果转换为图形及图像在屏幕上显示出来并进行交互处理的理论、方法和技术。它涉及到计算机图形学、图像处理、讨算机辅助设计、计算机视觉及人机交互技术等多个领域。数据可视化首先来自科学计算可视化( v i s u a l i z a t i o ni ns c i e n t i f i cc o m p u t i n g ) ，科学家们不仅需要通过图形图像来分析由计算机计算出的数据，而且需要了解在计算过程中数据的变化。随着计算机技术的发展，数据可视化的概念大大扩展，它不仅包括科学计算数据的可视化，而且包括工程数据和测量数据的可视化。学术界通常将这种分布在三维空间内的空间数据的可视化称为体视化，也叫做三维数据场可视化m j 。c t 技术f 是数据可视化技术的一个应用，它是物理学和计算机科学发展的产物。电容层析成像虽与c 1 技术在传感器检测方法上不同，但其后在层析成像、运用特定的重建算法原理上是相同的，都是将传感器检测的一组离散的一维投影数据经计算机处理后，再现被测对象真实写真。尽管计算机断层扫描已广泛应用于疾病的诊断，但这些医疗仪器只能提供人体内部的二维图像。医生们只能凭经验由多幅二维图像去估计病灶的大小及形状，“构思”病灶与周围组织的三维几何关系 i l l 。与之相对应，目前，在工业上，电容层析成像系统的研究也局限于离散相分布的二维图像重建上，将电容传感器测量的电容量经计算机处理后，利用计算机图形学和图像处理技术，变成二维图像的灰度值，再现流动的两相流体的流型及离散相浓度。组测量数据对应一幅二维图像，多幅二维图像间没有相关性，无法展现被测对象内在运行的动态过程 12 】。科学计算可视化技术的核心是三维空间数据的可视化m 】，它可以由一系列二维图像重构出三维形体，并在计算机上显示出来。在医学上，可实现矫形手术、放射治疗等的计算机模拟及手术规划；在多相流检测系统，可实现将测量数据变成图形数据、将工业现场的“黑箱”系统变成透明系统，将不稳定流场的迹线直观的表现出来，激发人的形象思维，为科学发现、工程开发、工程诊断和业务决策提供依据；供控制者对被测对象内在运行情况、变化规律进行深入的了解和掌握。哈尔滨理工大学工学硕士学位论文1 4 2 三维数据场可视化技术的特点计算机用于科学计算和数据处理已有近5 0 年的历史。但是，由于计算机技术的限制，数据只能以批处理而不能进行交互处理：不能对计算过程进行干预和引导，只能被动地等待计算结果的输出；而大量的输出数据也只能采用人工方式进行处理。不仅不能及时地得到有关数据的直观、形象的整体概念，而且还有可能丢失大量信息。三维数据可视化技术与传统的计算机数据处理相比，具有三大明显特点：11 4 ，2 1 交互性一组基础数据j 乏源于基层或一次仪表的测量数据，它可供系统完成各种管理、控制及输出、打印等功能，这是将数据转化为信息的过程。知识可分成两大类，一一类是不可啦用语言和文字来描述的，称之为隐知识：另一类是可以用语言和文字来描述能，称之为显知识。显然，信息就是显知识。而数据可视化系统怎样来分析和处理大量复杂的数据呢? 答案是要提供象人眼一样的直觉的、交互的、反应灵敏的可视化环境，将数据和信息用图形和图像表示出来，为用户提供获得隐知u 的创造条件，使用户可以方便地以交互的方式管理和开发数据。1 4 2 2 多维性用户可以看到表示对象或时间的多个属性或变量，而数据可以按其每一维的值，将其分类、排序、组合和显示。1 4 2 3 可视性数据可以用图像、曲线、二维图形、三维体和动画等形式来显示，并可对其模式和相互关系进行可视化分折。人的创造性不仅取决于人的逻辑思维，而且取决于人的形象思维。海量的数据只有通过可视化变成形象，才能激发人的形象思维。从在表面上看视杂乱无章的海量数据中，找出其中的规律，为科学发现、工程开发、医疗诊断和业务决策提供依据。1 4 3 三维数据场可视化现状与发展趋势科学计算可视化( v i s u a l i z a t i o ni as c i e n t i f i cc o m p u t i n g ) 是发达国家在2 0世纪8 0 年代后期提出并发展起来的一个新的研究领域。1 9 8 7 年2 月，美国国家科学基金会在华盛顿召开了有关科学计算可视化的首次会议。会议指出：科学家们不仅需要分析计算机得出的计算数据，而且需要了解在计算过程中数据的变化，并要将这种变化显示出来，激发人们的潜在意识，为科学发现提供感性认识。这些都要借助计算机图形学及图像处理技术，因而“将图形与图像处理技术应用于科学计算是一个全新的领域”。会议将这一领域定名为“v i s u a l i z a t i o ni ns c i e n t i f i cc o m p u t i n g ”：自1 9 9 0 年起，每年举行多次“信息可视化”国际学术会 义；自1 9 9 5 年起，每年举行一次国际学术会议。近年来，科学计算可视化又扩展到工程计算可视化及测量数据可视化；随着网络技术和电子商务的发展，信息可视化又被作为新的课题提出。通过数据可视化技术，发现大量测量数据中的隐知识、发现金融、通信和商业数据中隐含的规律，从而为决策提供依据，已成为数据可视化技术中新的热点。目前，科学计算可视化的应用领域己涉及医学、地质勘探、气象学、分子模型构造、计算流体力学、有限元计算等各个领域m ) 。我国科学计算可视化( v i s u a l i z a t i o ni ns c i e n t i f i cc o m p u t i n g ) 研究自2 0 世纪9 0 年代，最早清华大学、北西中科院计算所、软件所和浙江大学分别开展了科学计算可视化研究，在医学、石油、地质勘探、地震预报等领域取得了可喜的成绩。在这两校及中国科学院的推动下，天津大学、第三军医大学、国家气象中心和哈尔滨理工大学等高等院校和科研单位也相继开展研究。迄今，已有较成熟、针对不同应用领域的三维数据场可视化算法、应用实例及丰硕研究成果。最为代表的有2 0 0 3 年7 月笫三军医大学报发表的“中国数字化可视人体”m 】，文中介绍了计算机三维数据重构方法及心脏各部位可视化结果，图像形象、逼真、立体效果好，这无疑为科学计算可视化技术的进步发展打下了坚实的基石出。信息可视化应用依赖于人的视觉和感知系统天生就能迅速反映、理解、判断的基本特征。信息系统的设计人员利用这些特征来提高显示信息的数据密度，通过人天生的感知系统，发现数据间的内在规律、变化趋势、数据的异常值或错误，准确找到最小值和最：t 值以及识别聚类( c l u s t e r ) 。再通过人的认知系统的分析、判断，更好地理解复杂系统，提高决策能力，挖掘原本可能不知道的信息m i 。随着信息化社会的到来，信息处理量越来越大，未来信息可视化向以下五个方面发展一：( i ) 三维数据的可视亿。科学可视化的许多应用是三维可视化，因为科学可视化的用途主要就是表示实际的三维物体，这种计算机模型为科学家提供了进行操纵及实验以及预测真实物体实际行为的方法。近年来，三维可视化应用于众多方面，特别是在地质、建筑和医学领域、虚拟现实建模语言( v r m l ) 和数字成像等技术用来构建能够逼真显示三维数据的系统。有了这些系统，人们就能分析及研究三维物体和空间，这种方式比使用实际物体来得更实用、更有效。( 2 ) 多维数据的可视化。在信息可视化环境中，多维数据是描述有三种以上属性的物体的数据，每一种属性在视图上大体相同。譬如一个数据库，列有图书馆各种图书信息，每一图书具有多种属性。在诸多属性中，其应用目的是让用户可以用其中任何一种属性对图书信息进行查找、排序、组合，那么这些数据就是多维数据。将这些多维数据按其属性进行可视化处理，称多维数据的可视化。马里兰大学人机交! i 实验室的研究人员开发出了一种名为动态查询的框架，使用者就可以对这种多维数据进行可视化。( 3 ) 时态数据的可视化。以图形方式显示随着时间不断产生的数据是可视化信息最常见、最有用的方法之一。随着时间线，使用户一眼就能够看出事件前后发生时的持续情况，以及哪些事件与其他事件相关，哪些事件与其他事件无关。( 4 ) 层次数据的可视化。层次数据即树形数据就是有这样一种内在结构的数据，每个项日或节点都有个父节点( 最上面的节点即根节点除外、。节点分兄弟节点( 拥有闻一个父节点的节点) 和子节点( 从属某个父节点的：声点) 。层次结构相当常见，商业组织、计算机数据存储系统都是按树形结构排列的层次数据。层次结构数据可视化很常见，但层次数据的三维视图有时会遇到这样一个问题：每个节点的大小和内容都是隐藏的。为了克服层次结构数据可视化所遇到的这问题，可将树形图用嵌套矩形来表示数据层次。所有矩形都放在一个很大的边界矩形罩面，里面的每个矩形表示一个节点。如果是父节点，里面还会有子节点，分别对嵌套矩形表示数据可视化，这种设计使树形图用一一个视图就能够显示数据层次罩面的所有节点。( 5 ) 网络数据的可视化。网络数据指与任意数量的其他项目有着关系的项目( 有时又叫节点) 。因为网络数据集里面的节点不受与它有关系的数量有限的其他节点的约束( 不像层次节点，它们都有惟一的父节点) ，所以网络数据与层次数据不同，没有固有的层次结构。两个节点之间可以有多条路径，项目与项目之间的关系其属性数量部是可变的。因此属性和项目之间的关系可能非常复杂，所以与其它数据可视化相比，如果不用某种号门可视化方法，网络数据很难显示。人的创造性不仅取决于人的逻辑思维，而且取决于人的形象思维。信息可视化能够以多种方式帮助人们更有效地获取信息。可视化系统使人们能够浏览大量的复杂数据，迅速找到所需信息，更加轻松地浏览数据及与之互动、识别的模式及趋势，以及更好地理解信息。信息可视化的未来，定会使人类社会变得更加丰富多彩。哈尔滨理上大学工学硕士学位论文1 5 课题来源及研究的主要内容随着流动层析成像技术的逐步发展，其诸多优点决定了该技术在解决多相流检测问题上有巨大的发展潜力和广阔的工业应用前景，但是电容测量的精度同题、软场问题( 即灵敏度分布受被测多相流相分布的影响) 以及图像重建质量不高的问题，严重阻碍电容层析成像技术实用化进一步发展，为此进行研究。本文研究得到黑龙江省自然科学基金( f 0 1 2 5 ) 、黑龙江省教育厅计划基金( 9 5 5 1 1 0 6 ) 的资助。本课题以8 电极电容层析成像系统为对象，研究了油水两相流的测问题、敏感场的分析与计算问题、8 电极电容层析成像系统三维数据场的可视化问题。本文主要研究内容包括以下几方面：( 1 ) 电容层析成像系统传感器原理及系统工作原理；( 2 ) 以有限元模型为基础的场域剖分、敏感场分析及计算；( 3 ) 以电容测量数据及敏感场分析数据为基础数据进行可视化映射；( 4 ) 对油水两相流离散帽分布进行动态显示。第2 章电容层析成像系统的组成及原理2 1 电容层析成像系统的组成电容层析成像( e l e c t r i c a lc a p a c i t a n c et o m o g r a p h y ，缩写为e c t ) 系统主要由电容传感器、数据采集及信号处理、图像重构三部分组成，如图2 1 所示。f i 9 2 【a r c h i t e c t u r e o f f a s t s y s t e m2 1 1 电容传感器电容层析成像系统的电容传感器与普通测量用的电容传感器有所不同，它的结构通常采用在绝缘管道外侧均匀地安装多对铜极板所构成，为防止外界电磁场的干扰，电极板的外部设置接地的屏蔽罩。同时，为了扩大系统的动态范围，降低相邻极板间的固有电容，在各极板间设置了接地的径向保护极板，结构如图2 2 所示。电容极板径向屏蔽极板绝缘管道外屏蔽景图2 - 2电容层析成像系统传感器示意图f i 9 2 2 t h e i l l u s w a l i o n o f e c ts y s t e ms e o s o r目前常用的e c t 系统电容极板数目是8 、1 2 、1 6 及2 0 等，本文以典型的8电极电容传感器为研究对象。以8 个极板中的任意一个极板为起点，逆时针对8个极板编号1 ，2 ，8 。在一个完整的测量过程中，首先选择电极板1 为源极板，并给源极板1 加一定的正电压e ，分别以电极板2 ，3 ，8 为检测极板，测量极板对1 与2 ，1 与3 ，1 与8 之间的电容值，每次测量时其它闲置电极都接地。接着选择电极2 为源电极，测量极板对2 与3 ，2 与4 ，2 与8的电容值，由于上一步测量中电极1 和2 之问的电容已经测量过了，因此这次测量中电极1 是冗余的。再阱电极3 作为源极板，依此类推，直至测量完电极板7 与8 之间的电容值。这样，在8 电极板系统中总共可获得2 8 个独立的测量值，为便于整体计算，将上述测量值编号为c ，c ：，c ：。1 1 9 。除了用序号表示所测量电极对的电容值外，本文还采用c 来表示极板对叫问的电容值。对具有个电极板的系统，可得到的独立电极对总数n 为：n = c ：= 掣( 2 - 1 )当n = 8 时，一个采样周期的电容测量数据如表2 1 所示：表2 1 电容测量数据( 一个采样周期)t a b l e 2 - 1t h e m e a s u r e d d a t a o f e l e c i c a lc a p a c i t m a c e ( i n a s a m p l i n g c y c l e )c , jc 】2c 1 ，3c l ，i c 1 5c 1 ，6c 1 7c 1 ，8c 2 ，3c 24c 2 ，5c 2 ，6c 2 ，7c 2 ，8c 34c 3 ，5c 36c 3 7c 3 ，8c 4 ，5c 4 ，6c 4 7c a 8c 56c 57c 58c 67c 68c 7 ，8管道内不同位置的介质对同一极板对的电容值影响程度不同，这一性质称为软场特性。某一个电容极板对之间的电容值，实际上可近似地看作是管道截面内所有点对该电容不同程度的迭加，因此，每对电极电容的测量实质上是对管截面内多相流体的扫描，其电容测量值就可看成管内多相流体在某一方向或角度的投影数据m 。本课题的研究，就是对管内多相流体进行2 8 个不同角度的扫描，从而获得2 8 个投影数据后，利用这些投影数据来获取多相流的各相参数，从而进行图像重构。哈尔滨理工大学工学硕士学位论文2 1 2 数据采集及信号处理该系统包括多通道数据的巡回采集控制、电容电压( c v ) 转换、a d 转换及通讯接口等，其中电容l 乜压( c v ) 转换电路的设计是最关键也是最困难的。成像所依据的数据就是由该系统采集并传输给计算机的数据，成像的准确性很大程度上依赖于该系统的准确性。这也是电容层析成像技术发展过程中的主要难题之一，也是研究的一个主要方向。数据采集及信号处理系统中的多通道数据采集控制由多通道扫描电路完成，该电路的作用是根据讨。算机指令的要求选择相应的两个电极作为测量电极，测量其电容值，同时使其它闲置电极可靠按地。该系统的a d d 转换及通讯接1 ：3的作用是完成a d 转换，并且负责对原始数据的模拟滤波和预处理，阻及与计算机进行通信。2 1 3 图像重构这部分工作是由数据采集及信号处理的软、硬件完成的。采样的数据被送到计算机，经处理后，为了展现工业现场( 黑箱系统) 的真实动态过程，必须由这些样本重新产生原始的连续信号，这一过程称为图像重构m t 。图像重构部分采用的是p c 机，它主要负责对接口电路发出指令，控制数据采集系统实时采集数据，并将数据予以存储，然后采用相应的算法进行图像重构和图像显示。p c 机显示器屏幕上的二：维图像是有由存放在帧缓存中的二维离散信号经图形重构而成。因而，可视化映射算法就应将离散分布的三维数据场，按照一定的规则转换为图形显示设备帧缓存中的每个像素点的灰度值，以实现电容层析成像系统三维数据场可视化及计算机仿真。2 2 电容层析成像系统的工作原理在实际测量时，极板间的电容测量值实际是由两部分所组成的，一部分是极板与管道内部的电介质形成的电容，另一部分是极板与屏蔽罩间屏蔽层的电介质所形成的电容。通过合理设计电容传感器，使第二部分的电容值很小，以至于可忽略其对整个电容测量值的影响m 】。以两相流为例，设两相流体的等价介电常数为s ，离散相和连续相的介电常数分别为和s ，两种介质所占的体积分别是k 和。假设离散相均匀地分布在连续相中，则两相流体的等价介电常数s 为：s = ( 1 ，v ) 毛+ ( ) s ，( 2 2 )哈尔滨理工大学工学硕士学位论文式中、, f f ，为己知常数，、，为两相流的总体积。令扩= _ + ( 2 3 )则由上式可推导出：f 2 ( v j 一7 功e 1 + ( 弘v i ) v l f2=e2 ( v t 功( f 1 一e2 )( 2 4 )又令风= k v( 2 - 5 )即口。一表示离散相浓度。这样，可得电容测量值c 为：c = k 馏= k - 厂( 嵋v ) = k - 厂( 风)( 2 6 )式中芷一由电容传感器结构尺寸所决定的特征常数。由此可见，电容测量值( 的大小可作为浓度厌测量的依据。显然，要获得多相流体流经管道时相浓度分布的图像，就必须获得关于多相流体的多个电容测量值。因此，在电容流动层析成像系统中的电容传感器采用的是多电极阵列式电容传感器，传感器各电极之间的相互组合可提供反映多相流体相浓度分布的多个电容测量值，再以这些电容测量值为投影数据，采用某种图像算法，即可重建管道截面的图像。在研究电容层析成像系统中，如果忽略管道径向的多相流体各分相的分布变化，即认为沿管道在极板长度的范围内，每一个截面的介质分布均相同。在极板长度较短时，这种假设具有其合理性。在这一假设前提下，根据电学原理，并忽略管道壁以外部分屏蔽层的影响，对于8 电极阵列式电容传感器中任意两个电极构成的电容，其测量值c ，为【2 3 j ：c ，= 怍( x ，y ) s 廖，y ，( x ，y ) ) a x 咖( 2 7 )分式中7 = 1 ，2 ，2 8 ；卜一管道截面；e ( x ，y ) 一管道截面内电介质分布函数：s ，( x ，y ，占( t y ) ) 一该测量电容c 1 的灵敏度分布函数，即电容c ，对点( x ，y ) 处介质的敏感程度。上式实际上是把电容作为管截面上各点介电常数的加权和，而某点的加权值的大小反映了电容对此点的介质变化的灵敏程度。对某电容来讲，它对传感器空间内不同的点有不同的灵敏度，它们构成了一个“灵敏场”。从灵敏度分布函数的定义中可以看出电容c ，的灵敏度分布函数受到电介质分布e ( x ，y ) 的影响，即当管道内电介质分布不同时，管道内同一点对电容c ，的影响程度也不同，这就是电容层析成像系统的软场特性。基于测量电容c ，的严格表达式( 2 7 ) ，为了便于求解，在管道截面非常小的区域内，假设电介质分布函数e ( x ，y ) 对测量电容值c 的灵敏度分布函数s ，( x ，y ，e ( x ，y ) ) 的影响很小，可以忽略。这样，在8 电极系统中，有2 8 个独立的电容，可用一个2 8 维的向量 c 1 c 2 ，c 2 8 表示，相应地有2 8 个灵敏度函数s 1 ( x ，y ) ，s 2 ( x ，y ) ，s 2 8 ( x ，y ) ，则在8 电极系统中，根据式( 2 7 ) 有【2 4 1 ：阶c ：，( ? ：。】- 肛( ) s 。y ) a x d y= 肛( m 鼢y ) d x d y2 _ 8 = j ( 训) s ：。( x ，y ) d , ：d y图像重构的目的就是根据式( 2 - 8 ) 求解出e ( x ，y ) 。由于e ( x ，y ) = 厂( 卢s ) ，所以根据2 8 个已知电容测量值c 1 c 2 ，q 8 ，利用灵敏度作为先验信息求解电介质分布e ( x ，y ) ，即求出离散相浓度，最终也就实现了两相流系统的图像重构。由于一个采样周期电容测量值的个数较少，只有2 8 个，少于成像像素个数，需通过合理的敏感场变换进行插值，由此计算出管道截面内的介质分布图像。2 3 本章小结本章主要介绍了电容层析成像系统中电容传感器、数据采集及信号处理、图像重构三个部分的组成原理及在系统中所起的作用。在此基础上，介绍了电容层析成像系统的工作原理，分析了电容传感器检测数据、电容传感器灵敏度分布数据以及两相流离散相浓度之间的关系，推导出计算公式。3 1 引言第3 章敏感场的分析与计算敏感场( s e n s i t i v i t yf i e l d ) 又称为电容测量值c ，的灵敏度( s e n s i t i v i t y ) 分布函数，即电容c ，对点( x ，y ) 处介质的敏感程度。它具有软场特性，也就是说敏感场受被测介质分布的影响，对某电容c 来讲，它对传感器空间内不同的点有不同的灵敏度。一旦电容传感器设计完成、两相流流体确定，则软场特性就确定。所以，灵敏度( s e n s i t i v i t 5 ，) 分布数据作为图像重构所需的先验数据被使用，在电容层析成像系统三维数据场可视化及计算机仿真系统软件运行前，应通过传感器结构上各微元分布电窖c ，( 七) 及离散相和连续相的介电常数和，计算出灵敏度( s e n s i t i v i t y ) 分佰函数。灵敏度( s e n