结合有限元与边界元方法求解心电图逆问题的研究.pdfVIP

电工数学进展 =aaa＆殛!堡签迪呈iQ睦堑M趔楚g碰壁 !塑 结合有限元和边界元方法求解心电图逆问题理论研究 何为 姚德贵 (重庆大学电气工程学院．重庆400044) 摘要 本文摄出了一种结台有限元和边再元方法求解心电田逆闫曩的曩学模型瓤方法；文中推导了有曩元翔边界 元站台的数学公式，并应用一十三维球型攮塑进}亍了模拟计算． 美t嗣z逆向囊有限元法边界元法 1 引言 随着描述电磁场基本性质的方程组的出现，电磁场问题的求解很大程度归结为求解一组给 定边界条件下的偏微分方程组，即所谓边值问题．按电磁场经典理论划分，有三类边值问置： Di触、№珊姗和混台边值问题。如果场域确定，边界条件按E述三类边值闻题限定，则可 获得这组偏微分方程的唯一解．偏离电磁场经典理论划分的边值问题，则被认为是无解或无意 义【u。但是。随着科学技术的发展．在地质天文热物理以及生物医学工程诸如心电圈逆问曩等， 如当场域的部分边界条件是可测量的，另一部分边界是不可触及的，而恰恰这部分无法■及的 边界上的电位分布或电流密度是待求量(心电图逆问题)，这是第一类特殊类型的边值问l妒1． 第二种为已知部分边界上的边界条件，求解界定场域内部的不同物体的导电率分布(电磁法物 探和电阻抗成象)州．这些类型的边值问题均属于电磁场的逆问愿．本文重点就第一类型涉及的 心电图逆闻题的数值计算方法和模型进行探讨．根据心电图逆问题中人体模型包含各向异性导 电性肌肉层的特殊结构，文中提出了一种结合有限元和边界元求解心电图逆问蜃的方法，并应 用一个三维球域的例子进行了模拟计算，其结果令人满意． 2胸腔电流场数学模型 ．人体胸腔削面可以帮助我们理解场域的划分和界定．图l是一个经过适当简化的人体胸腔 削面．在图l所示场域内，＆代表人体表面，s．代表肌肉层内某一表面，s。代表肌肉层表面， 而＆代表心外膜表面． ． 本文提出的结合有限元和边界元的方法是根据组成场域媒质的导电特性不同，把人体胸腔 划分为两个区间，分别应用不同的方法进行场域分割，再进行组装求解，这样可以充分发挥有 限元和边界元方法不同的优势． 2．1有限元膏教模型 有限元方法处理包含各向异性导电媒质的场域有巨大的优点，而人体胸腔被一层具有备向 异性导电性的骨骼肌肉层包围，所以，我们选取有限元法对骨骼肌肉层进行处理．骨骼肌的各 ．91． 图1．人体胸腔音尊面 向异性导电性体现在肌肉纤维捧列方向，沿纤维方向具有较高的导电奉，其切线方向期导电奉 较低，两者差别为15倍．我们可以应用一种特殊的方法处理各向异性导电性。其方法可以参考 文献阎．在每一个剖分单元里，其纤维的排列可以近似处理为一致．只要单元捌分足够搪细和纤 维方向分块足够准确。这种处理方法完全可以满足模拟计算要求． 描述骨骼肌肉层里的电流场方程为； V·pⅣ矿(r，只力=o G咖方程) (1) 以=O C人体表面釉 (2) 这里电导率【o】是—个张量．甲是电位函数，Jb是体表法向电流密度．方程(2]lf已裹人体表面和空气 是绝缘的。即传导电流不能从人体表面流入空中．应用有限元法对肌肉层区问的离散形式为： B·． 矗一 苫L叫．【盯】vNj吼dQ=善L虬厶凼 (3) fi=1．2．3．．．…．n。：j2j．2，3．．．…，n．}。《p=1．2．3，…．。ns：q2l。2，3，…．．H1) N为插值函数．民为空间三维有限单元总数，～为单元内的节点数，吼代表单元节点电位． So为包围臂骼肌肉层的全体表面(和S1)边界单元数，～为表面单元的节点数，而JI代表