论文部分内容阅读
准确阐明正常与各种异常心电图如何产生于心肌细胞的电活动和考察各种心律失常的形成与持续机制对于心脏疾病的诊断具有重要意义。由于受临床与实验手段的限制,对于心电学领域中存在的许多不清楚和有争议的问题,建模与仿真一直是一个受到高度重视的研究手段。迄今为止,大多数的心电模型以经典的Miller-Geselowitz模型为基础,使用偶极子表达抽象的心肌细胞的电活动,并由偶极矩计算场点电势得到仿真心电图。其它重要的模型种类包括电生理模型和使用Fitzhugh-Nagumo方程考察心律失常动力学性质的模型。已报道的电生理模型都是一维和二维的多细胞模型,旨在详细考察细胞间连接通讯与兴奋传播的关系,未与体表心电图联系起来。大规模并行计算是发展全心脏电生理模型的主要障碍,所涉及的问题有两个方面,即如何高效地执行和如何方便地描述几千至几百万组心肌细胞动作电位方程的并行求解。 细胞自动机是一种完全并行的计算模式,具有固有的并行性,作为离散动力学仿真系统近年来受到广泛注意。尽管按标准定义的细胞自动机是基于规则的,但可以有基于语言的实现,使之具备更灵活的描述能力。作为建模的前提,我们对一种基于语言实现的细胞自动机系统Cellular3.0在语言和可视化设施方面进行了扩展,使之能够描述大规模并行数值计算。 为了在细胞与亚细胞层次精细地和可视化地考察心电图的产生与心律失常的形成,我们使用经扩展的Cellular3.0为工具设计与实现了一个全心脏电生理模型,该模型的两个构成要素分别是细胞自动机式的大规模并行计算和各种心肌细胞的动作电位模型。所包含的心肌组织为窦房结、房室结、心房肌、心室肌、心房传导束以及心室传导束。限于我们现有计算机的性能,模型是二维的,包含五千余个细胞。根据心室的分层特性以及不同层次细胞具有不同电学性质的事实,我们提出了一个与模型解析度和几何特性无关的运行时心壁分层算法,将室壁分层以及使心肌细胞的电活动与其所在层次联系起来。兴奋在同层细胞间是沿轴向端到端传导,在不同层细胞间是沿径向边到边传导,两种传导具有不同的速度。细胞间共有24种跨缝隙连结结构,用于准确描述各种细胞间不同的兴奋传导特性。靠并行求解数千组心肌细胞动作电位方程和计算各个细胞对之间的跨结兴奋传播,心电活动的时空过程,包括许多心律失常的形成与演绎过程,能够得到可视化地展现。单个心肌细胞的电活动可使用跟踪窗口进行跟踪与显示。对于建模中所遇到的问题我们也进行了讨论。 基于全心脏电生理模型,我们还提出了一个新的心电图仿真算法,在细胞层次计算每一个处于相互连接与通讯中的心肌细胞对各种导联位置心电图的贡献,仿真的正常与多种异常心电图波形与临床记录十分一致。根据仿真结果,我们讨论了若干心电图波形的形成机制及诊断意义,包括右胸R波的意义、病理性Q波的形成机制和一 侧 心室电活动对对侧胸表心电场的影响,处于不同位置心肌细胞动作电位波形变井’3 心电图波形变异之间的关系得到形式和定量的说明。已完成的仿真表明该模型对十考 察心律失常的形成具有极大的助益,论文给出了若干心律失常的仿真例,包括缺血诱 导的心动过速及M纠1胞独特的电学性质在心律失常形成中。】能起到的作川。 最后,我们简要讨论了大规模多层次并行i;算在解决广泛存在的复杂性问题方面 的作川,心电问题只是生物医学领域中复杂性问题的一个典型例了。