电阻抗成像电极系统优化设计仿真探究-

1、电阻抗成像电极系统优化设计仿真探究*【摘要】eit电极系统是整个eit系统最为敏感和关键 的部分之一。本文在线电极强制等势点模型的基础上，以复 合电极为例，提出了一种电极结构及参数优化设计方法，并 建立了电极结构参数优化设计仿真研究平台。可针对不同电 极结构和各种参数变化给出其对eit成像质量和检测灵敏度 的影响，进行优化设计，从而为实用化eit系统电极结构设 计提供理论依据。【关键词】eit电极；电极结构；强制等势点模型；仿 真研究；优化设计abstract】 the electrode system is one of the most sensitive and crucial part

2、s in eit system. in this paper a method for eit electrode structure and parameter optimizing design ， taking example for compound electrode, has been presented, which is based on the coercive equipotential node model of the line electrode. an imitation research platform of electrode strueture and pa

3、rameter optimizing design has been developed by the methods， the variety influences of electrode structure and parameters on reconstrueted image and measurement sensitivity can be obtained this will provide theory basis for the electrode construetion optimizing design of real eit system.key words】 e

4、it electrode； electrode structure； coercive equipotential node model； imita/tion research； optimizing design电阻抗断层成像技术(eit , electrical impedance tomography)是当今生物医学工程学重大研究课题之一。它 是继形态、结构成像之后，于近二十年才出现的新一代更为 有效的功能成像技术1, 2,具有无损伤、功能成像和医 学图像监护三大突出优势3,是一种理想的、具有诱人应 用前景的无损伤医学成像技术。eit测量中，电极直接与人体接触，位于系统的最前端。 在电

5、极上发生的事件，包括有用信息、噪声、伪差、接触阻 抗、极化电压等，都会作为信号进入后续电路被放大、传输, 参与信号处理过程，影响图像重建结果。电极系统结构及其 性能对于eit前端信息的有效提取、系统适时性和图像分辨 率的影响，特别是对eit检测灵敏度较差的中心区有用信息 的提取影响非常大，是整个eit系统最为敏感和关键的部分 之一，也是eit技术走向临床应用，向实用化研究发展必须 解决的问题。在电极材料和电极数确定后，如何根据确定的应用目标 和成像要求，合适地选择电极结构参数，是构建真实eit系 统时必须认真解决的问题。由于eit成像区域多为圆形或球 型，一些重要的电极结构参数，如电极宽度和电

6、极间距，还 相互制约，影响因素复杂，给eit系统电极结构参数优化设 计带来了困难。1993年ping hua等采用有限元方法研究了电极接触阻抗 的影响，认为阵列电极总宽度为测量圆域周长的80%90% 时效果最好4。2002年王超等采用强制等势点有限元模型, 考虑激励电极和测量电极宽度的影响，进行了仿真研究，认 为电极覆盖比率为57. 1时效果最佳5。关于eit电极结 构及其参数选择，虽然国内外学者的探索性研究取得了一些 进展，但至今还没有看到较全面、令人满意的结果。在实际 eit电极结构设计中往往仍采用经验数据。eit的正问题和图像重建大都以点电极为基础。但在实际 的eit系统中，电极不可能是

7、一个点，而是具有一定大小的 面。电极面和被测对象接触，电极区域将被强制为等电势， 从而改变场域的电场及其分布规律。激励电极和测量电极面 积的大小会影响eit检测灵敏度和图像重建质量。显然在进 行eit电极结构设计时应采用具有一定宽度的电极模型，以 使研究结果更加符合真实情况。本文将在线电极强制等势点模型的基础上，建立了复合 电极有限元模型，提出了一种电极结构及参数优化设计方 法，并建立电极结构参数优化设计仿真研究平台。采用本文 的方法和研究平台，可针对不同电极结构和各种参数变化给 出其对eit成像质量和检测灵敏度的影响，进行优化设计， 从而为实用化eit系统电极结构设计提供理论依据。1 eit

8、电极有限元仿真模型eit系统中通常使用圆形或矩形片状电极，其在断层成像 平面的投影为与其宽度相等的线段，简称为线电极。医学 eit常采用结构更为复杂的复合电极，其在断层平面的投影 由3条线段组成，中间线段为测量电极，两边的线段为激励 电极。如图1所示。图1线电极和复合电极示意图ke = b(1)式中，k为有限元方程的系数矩阵；e为所有剖分 节点的电势矩阵；b为有限元方程右侧常数项，包含有限 元方程的边界条件。设no为有限元剖分的节点总数；对于j个激励电极，有j 组激励电极强制等势节点，构成j个集合equli (owiwj, iwn)。每个集合的元素为该组激励电极等势节点的节点编 号，每个集合存

9、在ml个元素，其中，每组中最小的元素为 min_equl i。类似的，对于测量电极，也有j组测量电极强制等势节点, 构成j个集合equ2i (owiwj, i£n)o每个集合的元素为 该组测量电极等势节点的节点编号，每个集合存在m2个元 素，其中，每组中最小的元素为min_equ2i。首先，进行列合并列合并完成后，将klj(l=l , 2 , . , no , juequli-min_equli, iej)和 ki j (1=1, 2, . . . , no, jeequ2i-min_equ2i, i£j)删除，未被删除的列前 移，补进删除后的空列。然后，进行行合并,行合并

10、后，将 kjl(l=l ,2 ,. , no ,j£equli-min_equli, iej)和 kjl(l=l, 2, . . . , no, jwequ2i-min_equ2i, iej)删除，对空位进行前移 补充。最后，右侧常数项合并，右侧常数项合并后，将 bj(jeequli-min_equli, iej)和 bj(jgequ2i-min _equ2i, iej)删除, 对空位进行前移补充。经过上述合并过程后，有限元方程变为线性方程组(2)即为考虑了测量电极和激励电极宽度的 eit复合电极的有限元模型。解上述线性方程组，即可获得场域中各节点的电势。本方法还可推广应用于其他复杂

11、结构的电极系统，以建 立相应的有限元模型。2评价函数d和s为了评价电极结构参数变化对eit检测灵敏度和图像重 建质量的影响，本文定义以下检测灵敏度函数s和图像重建 质量函数d。2. 1检测灵敏度函数s式中m为每一次激励所对应的测量次数，n为激励次数， sij表示第i次激励第j次测量的检测灵敏度：表示场域中某部分的电导率由。0变为0 1时，电极测 量电压由v0变为vi。s表示了 n次激励的检测电压灵敏度的平均值，它只与电 场分布相关，与硬件误差和成像算法误差无关。s越大表示 检测灵敏度越高。2.2图像重建质量函数d为定量评价eit图像重建质量本文定义如下函数d：式中，m为成像区域中剖分单元总数，

12、gs (p)为成像模 型中第p单元的灰度值，g i (p)为重建图像中第p单元 的灰度值。d表示各单元重建图像灰度值与成像模型灰度值 的差经归一化后的平均值。采用目标函数d,不需对图像进 行二值化处理，即可直接获得重建图像和成像模型之间的差 别，可灵敏地反映重建图像的质量。d值越小表示成像结果 与模型的差别越小，图像重建质量越高。3 eit仿真和图像重建软件平台本文相关的研究工作是在作者开发的医学eit仿真和图 像重建软件的基础上进行的，如图2所示。整个软件采用 visual c+ +语言编写8, 9,界面友好，功能完善，可 扩展性强，与其他程序接口良好。可以实现有限元的自动剖 分，可对有限元

13、方程解法、激励电极数、测量电极数、电极 类型、电极结构参数、激励测量模式等多个参数进行设置， 并根据剖分和参数设置结果进行电磁场数值计算、显示电场 分布、保存正向计算结果、进行eit图像重建、显示结果并 保存图像、进行图像后处理并计算目标函数等功能。图2 eit仿真和图像重建软件平台界面4结果与讨论采用式(2)表示的eit电极有限元模型，通过eit仿真 和图像重建软件平台，借助评价函数s和d,即可对eit电 极进行电极结构及参数影响和优化设计仿真研究。作为本文 方法的应用，作者对线电极，复合电极等进行了电极结构参 数影响与优化设计研究，并对线电极与复合电极结构进行了 性能对比研究。以下是根据这

14、些研究获得的有关eit电极设 计与优化的共性问题。具体研究内容和相关结果将另文详 述。4.1电极宽度、间距与成像质量和检测灵敏度关系的研 究与窄的激励电极相比，较宽的激励电极可以增大敏感场 较弱区域的电场强度，而且可以通过增大接触面积减少电极 -皮肤接触阻抗，有利于提高系统的灵敏度。但电极过宽也 会产生一些负面效应：(1)宽电极的使用，必然会使电极面 强制为等电势，从而影响场域内部的电场分布，使电场更加 扭曲。激励电极越宽，强制等势面积越大，电场分布与点电 极相比扭曲的程度就越大，这将直接导致成像质量的下降。(2)电极间距增大，可使测量电流更加深入生物组织内部, 即探测深度会增加。这也是激励测

15、量模式的改变会显著影响 中心区灵敏度的主要原因之一，eit问题一般采用圆形边界 区域，场域和电极数确定后，电极的宽度和相邻电极间距之 和为一定值，使用宽电极必然会使相邻电极间距减小，这一 因素会在一定程度上减小测量电流的探测深度，使eit系统, 特别是中心区的灵敏度下降。4. 2电极数与成像质量的关系增大电极数可以增加独立 测量数，采集更多的测量数据，从而提高系统的总体成像质 量。但电极数增大至一定程度后，其作用就非常微小了。另 一方面。电极数增加势必提高对测量系统的要求，在硬件测 量精度一定的情况下，反而会使成像质量下降。4.3复合电极的结构参数优化医学应用eit可采用结构 较为复杂的复合电

16、极。复合电极包括激励电极的宽度，测量 电极的宽度，激励与测量电极间距和相邻电极间距等四个结 构参数。应该综合考虑各参数的影响，不能简单地用电极覆 盖率来处理电极优化问题。【参考文献】1 barber dc, brown bh. applied potential tomography. j.phsy. e. sci. instrum, 1984, (17):723-733.2 brown bh. medical impedance tomography and processimpedance tomography :a brief review. measurementscience and technology, 2001, (12): 991-996.3任超世生物电阻抗测量技术中国医疗器械信息， 2004, 10(1)： 21-254 ping hua, eung jw. using compound electrodes in electrical impedance tomography. ieee trans. biome eng, 1993,40(1)： 29-34.5王超，王化祥