电阻抗成像的数值模拟和分析计算数学专业论文.docxVIP

Shandong Shandong University Doctoral Dissertation 电阻抗成像的数值模拟和分析 李久平 (山东大学数学与系统科学学院，济南，2sol00) 摘要 电阻抗成像(Electrical Impedance Tomography,简称EIT)是一种新兴的生物医学 成像技术，它的具体做法是在人体表面设置大量电极，注入安全电流，测量人体表皮上 的电势分布，把这些数据传输给计算机，通过计算，得到人体内部的电阻率分布，并将其 转化为灰度或彩色图像．由于人体各部分电阻率不同，因此这实际上是一个人体内部的 图像．EIT在地球物理探矿，工业无伤检测，水下目标探测等领域也有着重要的应用前 景(【5，10，11，12，13】)，从而引起了人们广泛的兴趣，其研究涉及理论，数值计算和实验 等几个方面，横跨数学，物理学，电子信息，生物医学等多个学科． 电阻抗成像实际上是一类椭圆型偏微分方程反问题．反问题有别于传统的数学物 理方程的定解问题(通常称为正问题)，反问题研究由解的部分已知信息来求定解问题中 的某些未知量，正问题通常都是适定的，而反问题往往是不适定的，该特点也是反问题的 难点所在．基于实际应用问题的推动，20世纪60年代中期苏联科学院院士Tikhonov([2]) 提出了处理不适定问题的正则化方法，为不适定问题和反问题的研究开辟了道路．在中 国，反问题的研究始于20世纪80年代冯康院士的大力倡导． 正是由于非线性和严重不适定性，EIT问题的数值求解是非常复杂和困难的．目前 求解EIT问题的算法大体可分为三类： 1．基于全局线性化的非迭代算法， 2．直接法 3．解完全非线性问题的迭代法． 其中基于全局线性化的非迭代法可看作迭代法的一步迭代即终止而来，如Calderon算 法，等势线反投影法和一步牛顿法(NOSER)，其局限性在于只能求解和均匀分布较为接 近的电阻率分布，而且求解精度较差；直接法是指刻画在均匀分布介质上嵌入的分解法 或实现唯一性证明的直接方法，由于其需要在目标的整个边界上施加电流，测量电压， 用数学的语言来讲，就是要观测Neumann．to-Dirichlet映射，不适用于有限电极系统，因 而在实际应用上受到限制．基于有限元的迭代法目前是求解EIT问题的一类比较有效 的方法，但其收敛性分析一直没有得到解决．除了该问题的严重不适定性所带来的数值 解的稳定性差这一难题外，由于计算量大，难以实时成像是其另一难点．为解决这一问 题，Cheney([6])等人提出了著名的一步牛顿法(NOSER)．在该方法中，把均匀分布作为 Shandong Shandong University Doctoral Dissertation 初始猜测解，这样，均匀分布下的Jacobi矩阵可预先计算好并存储起来，实际计算时只 进行一步迭代，从而有助于快速求解．但这样的方法也限制了其适用范围和求解精度． 本文作者在袁益让教授的精心指导下，在大量的数值试验的基础上，选择Levenberg- Marquardt方法与偏微分方程数值解法相结合，对电阻抗成像进行数值模拟研究．L—M 方法是一种信赖域方法，也是一种正则化的迭代算法，因而可以用来稳定地求解非线性 反问题(【18，19】)．本文主要考虑两种电阻抗成像模型：连续模型(continuum model)和 电极模型(electrode model)．分别应用有限元、体积元和块中心有限体积方法在规则和 不规则区域，对二维和三维EIT问题进行求解．其中体积元和块中心有限体积方法是 首次用于求解电阻抗成像问题．对于连续模型，我们用精确解来进行数值模拟，验证模 型的正确性和算法的可靠性及可行性，然后应用到电极模型上去，并在电极模型下给出 了适用于以上三种数值方法(对称格式)的一类电流模式，将计算量最大的Jacobi矩阵 的计算彻底简化． 全文共分三章． 第一章介绍基于有限元方法的电阻抗成像的数值模拟和分析．§1．2—1．5研究了基 于有限元的二维EIT问题的数值模拟的方法和技术．采用了两种不同的网格：三角形 网格和四边形网，给出了伪单元刚度矩阵的计算公式，进行分析比较．在矩形区域上，采 用连续模型，对边界条件施以不同幅度的扰动，分别用不同空间步长的数值解和精确 解进行对比，验证了连续模型的正确性和算法的可靠性及可行性．数值实验表明，当步 长缩小时，反问题的数值解能逼近精确解；数值解关于Neumann边界条件的稳定性优 于Diriehlet边界条件，四边形网格优于三角形网格．§1．5在电极模型上的应用，研究两 种网格上EIT问题的各种数值表现．特别比较研究了固定发射电极模式和轮换发射电 极模式下的数值表现，结果表明轮换发射电极模式优于固定电极模式．在§1．6-1．9研究 了三维E