一种基于h细化的电阻层析成像有限元模型的制作方法
【专利摘要】本发明公开了一种基于h细化的电阻层析成像有限元模型，属于电学层析成像技术领域。针对电阻层析成像技术中灵敏度矩阵的病态性，在可有效提高正问题计算精度的有限元模型a的基础上，以h细化区域的起始层数、终止层数以及三角形有限元内部所插入节点的横坐标、纵坐标为变量，以敏感场均匀分布时灵敏度矩阵条件数的倒数为适应度函数，利用改进粒子群算法h细化有限元模型a，提出一种基于h细化的电阻层析成像有限元模型b。本发明在不影响正问题计算精度的前提下，改善了敏感场均匀分布时灵敏度矩阵的病态性，有效提高了图像重建质量，适用于以灵敏度理论为基础的电阻层析成像图像重建算法。
【专利说明】
一种基于h细化的电阻层析成像有限元模型
技术领域
[0001 ] 本发明涉及电阻层析成像(Electrical Resistance Tomography，ERT)有限元模 型，具体是一种基于h细化的电阻层析成像有限元模型，属于电学层析成像技术领域。
【背景技术】
[0002] 作为过程层析成像技术的重要分支，电学层析成像技术可分为基于电磁感应原 理、可重建敏感场域内磁导率分布状态的电磁层析成像(Electromagnetic Tomography， EMT)技术、基于电容敏感原理、可重建敏感场域内介电常数分布状态的电容层析成像 (Electrical Capacitance Tomography，ECT)技术、基于阻抗敏感机理、可重建敏感场域内 复导纳分布状态的电阻抗层析成像(Electrical Impedance Tomography，EIT)技术以及基 于电阻传感机理、可重建敏感场域内电阻率/电导率分布状态的电阻层析成像技术。
[0003] 为了提高电阻层析成像有限元计算精度，通常采用细化有限元的方法。目前常用 的有限元细化方法可简单概括为两类:h细化与p细化。其中h细化的思想是:始终采用同一 种单元，但是不断改变网格尺寸的大小，使其在部分区域变得更粗，而在另外区域变得更 细，从而提供达到所要求结果的最经济方案。根据具体实施方法，h细化又可细分为三类：
[0004] 第一种h细化方法是单元细分(Element Subdivision)。这种细化方法很容易实 施，得到了非常广泛的应用。单元细分可以在保持原有单元边界完整的前提下，简单的将那 些出现较大误差的单元细分为更小的单元。但是在带有中间节点的单元和没有中间节点的 单元相连处，由于细分会出现许多悬空节点，这时，需要在悬空节点处施加局部约束，增加 了计算量，从而使单元细分的实施显得繁琐。另外，在单元合并操作中所采取的复杂数据处 理将降低单元细分方法的效率。
[0005] 第二种h细化方法是一种完全的网格重构(Mesh Regeneration)。这种细化方法根 据已获得的结果，对所有区域的单元尺寸进行新的预测，再进行一次全新的网格划分。在网 格重构方法中，虽然网格细分与单元合并可同时进行，但计算量很大，特别是在三维问题 中，采用某些类型单元进行网格划分会很困难，而且从一种网格划分到另外一种网格划分 传输数据时也存在问题。
[0006] 第三种h细化方法有时也被称为γ细化（γ refinement)。γ细化并非一种真正的 细化方法，其基本思想是:在保证节点总数不变的前提下，通过调整节点位置以获得最佳的 近似结果。
[0007] 电阻层析成像逆问题病态的根源是灵敏度矩阵条件数很大。为了提高算法实时 性，目前基于灵敏度理论的图像重建算法通常采用固定敏感场均匀分布时灵敏度矩阵不变 的策略，而不同拓扑结构的有限元模型对应的灵敏度矩阵并不相同。图1中的有限元模型a 是在传统按等间隔原理剖分的有限元模型基础上，利用改进遗传算法优化得到的、可有效 提高电阻层析成像正问题的计算精度有元限模型。国内研究学者在有限元模型a基础上，在 每个三角形有限元形心的位置插入一个节点，将其对应的灵敏度矩阵应用于修正牛顿-拉 夫逊算法，并在算法重建过程中，遵循"计算正问题时采用细化前有限元模型，修正电阻率 分布时采用细化后有限元模型及其对应的灵敏度矩阵"的原则，改善了 Hessian矩阵的病态 性，在不影响算法实时性的前提下，有效提高了图像重建质量，但由于并未改善敏感场均匀 分布时灵敏度矩阵的病态性，因此，不适用于其他以灵敏度理论为基础的图像重建算法，且 三角形有限元形心并不是所插入节点的最优位置。
【发明内容】
[0008] 针对上述现存的技术问题，本发明提出一种基于h细化的电阻层析成像有限元模 型，在不影响正问题计算精度的前提下，可改善敏感场均匀分布时灵敏度矩阵的病态性，有 效提高图像重建质量，适用于以灵敏度理论为基础的电阻层析成像图像重建算法。
[0009] 为了实现上述目的，本发明采用的一种基于h细化的电阻层析成像有限元模型，有 限元模型分为采取h细化的区域与未采取h细化的区域两部分，具体是在可有效提高正问题 计算精度的有限元模型a的基础上，以h细化区域的起始层数、终止层数以及三角形有限元 内部所插入节点的横坐标、纵坐标为变量，以敏感场均匀分布时灵敏度矩阵条件数的倒数 为适应度函数，利用改进粒子群算法h细化有限元模型a，得到基于h细化的电阻层析成像有 限元模型b。
[0010] 所述基于h细化的电阻层析成像有限元模型建立的具体步骤是：
[0011] 步骤一:建立可有效提高电阻层析成像正问题计算精度的有限元模型a;
[0012]步骤二：以h细化区域的起始层数、终止层数以及三角形有限元内部所插入节点的 横坐标、纵坐标为变量，以敏感场均匀分布时灵敏度矩阵条件数的倒数为适应度函数，利用 改进粒子群算法h细化有限元模型a，得到基于h细化的电阻层析成像有限元模型b;
[0013] 步骤三:完成有限元模型b节点与有限元的编号。
[0014] 其中，有限元模型a是在传统按等间隔原理剖分的有限元模型基础上，以除最外层 之外的每一层半径为变量，以敏感场均匀分布时模型均方根值的倒数为适应度函数，并引 入三角形最长边与最短边的比值作为惩罚函数，利用改进遗传算法离线优化有限元模型拓 扑结构得到的。
[0015] 所述采取h细化的区域与未采取h细化的区域是由改进粒子群算法确定的。
[0016] 所述采取h细化是指在确定采取细化措施区域中的每个三角形有限元内部插入一 个节点。
[0017] 所述插入节点的横坐标与纵坐标是由改进粒子群算法确定的。
[0018] 所述改进粒子群算法是以h细化区域的起始层数、终止层数以及三角形有限元内 部所插入节点的横坐标、纵坐标为变量，以敏感场均匀分布时灵敏度矩阵条件数的倒数为 适应度函数。
[0019] 与现有技术相比，本发明针对电阻层析成像技术中灵敏度矩阵的病态性，在可有 效提高正问题计算精度的有限元模型a的基础上，以h细化区域的起始层数、终止层数以及 三角形有限元内部所插入节点的横坐标、纵坐标为变量，以敏感场均匀分布时灵敏度矩阵 条件数的倒数为适应度函数，利用改进粒子群算法h细化有限元模型a，提出一种基于h细化 的电阻层析成像有限元模型b。本发明在不影响正问题计算精度的前提下，改善了灵敏度矩 阵的病态性，有效提高了图像重建质量，适用于以灵敏度理论为基础的电阻层析成像图像 重建算法。
【附图说明】
[0020] 图1为有限元模型a即可有效提高正问题计算精度的有限元模型示意图，包含537 个节点、880个三角形有限元；
[0021] 图2为本发明提出的基于h细化的电阻层析成像有限元模型b示意图；
[0022]图3为有限元模型c即在有限元模型a基础上，以全部三角形有限元内部所插入节 点的横坐标、纵坐标为变量，以敏感场均匀分布时灵敏度矩阵条件数的倒数为适应度函数， 利用改进粒子群算法优化所得的有限元模型；
[0023]图4为计算敏感场边界电压测量值时所采用的有限元模型d示意图，包含3697个节 点、6944个三角形有限元；
[0024]图5为有限元模型a、b、c在敏感场均匀分布时正问题计算精度的比较示意图；
[0025] 图6为仿真实验中设置的介质分布示意图。
[0026] 图中：1、未采取h细化的区域，2、采取h细化的区域，3、采取h细化区域中三角形有 限元内部所插入的节点。
【具体实施方式】
[0027]下面结合附图对本发明作进一步说明。
[0028] 如图1和2所示，一种基于h细化的电阻层析成像有限元模型，有限元模型主要由采 取h细化的区域与未采取h细化的区域两部分组成，具体是在可有效提高正问题计算精度的 有限元模型a的基础上，以h细化区域的起始层数、终止层数以及三角形有限元内部所插入 节点的横坐标、纵坐标为变量，以敏感场均匀分布时灵敏度矩阵条件数的倒数为适应度函 数，利用改进粒子群算法h细化有限元模型a，得到基于h细化的电阻层析成像有限元模型b。
[0029] 建立基于h细化的电阻层析成像有限元模型b的具体步骤是：
[0030] 首先，建立可有效提高电阻层析成像正问题计算精度的有限元模型a;
[0031] 其次，以h细化区域的起始层数、终止层数以及三角形有限元内部所插入节点的横 坐标、纵坐标为变量，以敏感场均匀分布时灵敏度矩阵条件数的倒数为适应度函数，利用改 进粒子群算法对有限元模型a进行h细化，从而获得基于h细化的电阻层析成像有限元模型 b；
[0032] 最后，完成有限元模型b节点与有限元的编号。
[0033] 上述有限元模型a是在传统按等间隔原理剖分的有限元模型基础上，以除最外层 之外的每一层半径为变量，以敏感场均匀分布时模型均方根值的倒数为适应度函数，并引 入三角形最长边与最短边的比值作为惩罚函数，利用改进遗传算法离线优化有限元模型拓 扑结构得到的。
[0034] 所述采取h细化的区域与未采取h细化的区域是由改进粒子群算法确定的。
[0035] 所述采取h细化是指在确定采取细化措施区域中的每个三角形有限元内部插入一 个节点;所述插入节点的横坐标与纵坐标可由改进粒子群算法确定。
[0036] 所述改进粒子群算法是以h细化区域的起始层数、终止层数以及三角形有限元内 部所插入节点的横坐标、纵坐标为变量，以敏感场均匀分布时灵敏度矩阵条件数的倒数为 适应度函数。
[0037]如图5所示，采用图1中的有限模型a、图2中的有限模型b和图3中的有限模型c计算 电阻层析成像正问题，在敏感场均匀分布时，包含激励电极所在节点，图1中有限元模型a、 图2中有限元模型b和图3中有限元模型c均方根值均为1.4428%，去除激励电极所在节点， 三种有限元模型均方根值均为0.7374%，可见，与图1中有限元模型a和图3中有限元模型c 相比，图2中本发明提出的基于h细化的电阻层析成像有限元模型b不影响正问题计算精度。 [0038] 经验证，在敏感场均匀分布时，图1中有限元模型a、图2中有限元模型b和图3中有 限元模型c所对应的灵敏度矩阵的条件数分别为1.4661 X 106、9.4777 X 105、1.4108 X 106， 与图1中有限元模型a和图3中有限元模型c相比，图2中本发明提出的基于h细化的电阻层析 成像有限元模型b对应的灵敏度矩阵条件数分别降低了35.3543%、32.8204%，有效改善了 灵敏度矩阵的病态程度，从而有利于提高图像重建质量。
[0039] 为了验证本发明提出的基于h细化的电阻层析成像有限元模型b在提高图像重建 质量方面的有效性，设置如图6所示六种不同介质分布，并以改进Landweber预迭代算法为 例，在相同实验条件下(Duo T8100 CPU 3.OOGB内存2. IOGHz MATLAB 7.0)，将三种不同有 限元模型及其对应的敏感场均勾分布时灵敏度矩阵应用于改进Landweber预迭代算法。在 实时性方面，图1中有限元模型a、图2中有限元模型b和图3中有限元模型c对应的改进 Landweber预迭代算法耗时分别为0.15-0.16晕秒、0.31-0.32晕秒、0.46 -0.47晕秒，与图1 中有限元模型a和图3中有限元模型c相比，图2中本发明提出的基于h细化的电阻层析成像 有限元模型b基本不影响算法实时性。
[0040] 在电阻层析成像技术中，目前通常采用相关系数与图像相对误差评价算法图像重 建质量，表达式如式（1)、（2)所示。相关系数越大、图像相对误差越小，表明图像重建质量越 尚。
[0041] 三种不同有限元模型相关系数与图像相对误差比较如表1、2所示：
[0042] (1)
[0043] (2)
[0044] 式中：g为设置介质分布；i为重建结果;L为有限元数目；f与I分别为g与#的平均 值。
[0045] 表1三种不同有限元模型相关系数比较
[0047]表2三种不同有限元模型图像相对误差比较（％ )
LOO49」由表1、2?知，在相问买验条仵卜，在相夫糸数万_，三柙有I很兀模型相夫糸数的 平均值分别为0.6433、0.7460、0.6690，与图1中有限元模型&和图3中有限元模型(：相比，图2 中本发明提出的基于h细化的电阻层析成像有限元模型b相关系数平均提高了 15.9646%、 1 1.5097% ;在图像相对误差方面，三种有限元模型图像相对误差的平均值分别为 48.1012%、40.1110%、45.7794%，与图1中有限元模型a和图3中有限元模型c相比，图2中 本发明提出的基于h细化的电阻层析成像有限元模型b图像相对误差平均降低了 16.6112%、12.3820%，有效提高了图像重建质量。
[0050]由于实际中存在测量误差，因此将图4所示的有限元模型d计算正问题所得的敏感 场边界电压测量值加入5%幅度噪声，三种不同有限元模型相关系数与图像相对误差比较 如表3、4所示。
[0051 ]表3噪声干扰下三种不同有限元模型相关系数比较
[0055]由表3、4可知，在相同实验条件下，在相关系数方面，三种有限元模型相关系数的 平均值分别为0.6419、0.7438、0.6671，与图1中有限元模型&和图3中有限元模型(：相比，图2 中本发明提出的基于h细化的电阻层析成像有限元模型b相关系数平均提高了 15.8747%、 1 1.4975% ;在图像相对误差方面，三种有限元模型图像相对误差的平均值分别为 48.4246%、40.3814%、45.9447%，与图1中有限元模型a和图3中有限元模型c相比，图2中 本发明提出的基于h细化的电阻层析成像有限元模型b图像相对误差平均降低了 16.6097%、12.1087%，有效提高了图像重建质量，与无噪声干扰时结论相同。
[0056]综上所述，本发明提出的基于h细化的电阻层析成像有限元模型，在不影响正问题 计算精度的前提下，改善了灵敏度矩阵的病态性，有效提高了图像重建质量，适用于以灵敏 度理论为基础的电阻层析成像图像重建算法。
【主权项】
1. 一种基于h细化的电阻层析成像有限元模型，其特征在于，有限元模型分为采取h细 化的区域与未采取h细化的区域两部分，具体是在可有效提高正问题计算精度的有限元模 型a的基础上，以h细化区域的起始层数、终止层数以及三角形有限元内部所插入节点的横 坐标、纵坐标为变量，以敏感场均匀分布时灵敏度矩阵条件数的倒数为适应度函数，利用改 进粒子群算法h细化有限元模型a，得到基于h细化的电阻层析成像有限元模型b。2. 根据权利要求1所述的一种基于h细化的电阻层析成像有限元模型，其特征在于，所 述基于h细化的电阻层析成像有限元模型建立的具体步骤是： 步骤一:建立可有效提高电阻层析成像正问题计算精度的有限元模型a; 步骤二：以h细化区域的起始层数、终止层数以及三角形有限元内部所插入节点的横坐 标、纵坐标为变量，以敏感场均匀分布时灵敏度矩阵条件数的倒数为适应度函数，利用改进 粒子群算法h细化有限元模型a，得到基于h细化的电阻层析成像有限元模型b; 步骤三:完成有限元模型b节点与有限元的编号。3. 根据权利要求1所述的一种基于h细化的电阻层析成像有限元模型，其特征在于，所 述采取h细化的区域与未采取h细化的区域是由改进粒子群算法确定的。4. 根据权利要求1所述的一种基于h细化的电阻层析成像有限元模型，其特征在于，所 述采取h细化是指在确定采取细化措施区域中的每个三角形有限元内部插入一个节点。5. 根据权利要求4所述的一种基于h细化的电阻层析成像有限元模型，其特征在于，所 述插入节点的横坐标与纵坐标是由改进粒子群算法确定的。6. 根据权利要求3或5所述的一种基于h细化的电阻层析成像有限元模型，其特征在于， 所述改进粒子群算法是以h细化区域的起始层数、终止层数以及三角形有限元内部所插入 节点的横坐标、纵坐标为变量，以敏感场均匀分布时灵敏度矩阵条件数的倒数为适应度函 数。
【文档编号】G06F17/50GK105843984SQ201610143207
【公开日】2016年8月10日
【申请日】2016年3月14日
【发明人】肖理庆, 唐翔
【申请人】肖理庆