关富荣 李成乾 邓敏艺†

1) (广西师范大学物理科学与技术学院,桂林 541004)

2) (广西核物理与核技术重点实验室,桂林 541004)

1 引言

螺旋波是非线性时空系统远离热力学平衡态下自组织形成的波斑图[1],它广泛存在于物理[2]、化学[3]、生物系统[4]中,研究表明:生物系统中的螺旋波对生物的功能有着重要影响,如大脑皮层中出现的螺旋波在介观尺度上能组织和调节神经元的集体活动,既有助于皮层中神经元的正常信息处理,也能诱发癫痫病理性活动斑图[5].在老化的心脏中,螺旋波可自发地出现,但由于螺旋波是自维持波,其频率高于窦房结的起搏频率,因此当心脏组织中出现螺旋波时会导致心动过速,而螺旋波破碎成时空混沌将导致心颤,危及生命[6].大量研究表明,螺旋波是引起人类心房颤动[7,8]、心室颤动[9,10]等复杂心律失常的原因之一,由于存在潜在应用价值,心脏等可激发介质中螺旋波的动力学[11]、螺旋波的产生及其控制[12-14]已成为斑图动力学研究的热点问题.

心脏是典型的可激发系统,在进行体力活动时,心脏的跳动频率高于正常频率,此时心脏收缩和舒张的相对持续时间随之调整,以确保心室血液的充盈和正常循环.如果心跳频率增加而收缩期没有改变,会导致舒张期不呈比例减少,意味着下一次收缩前心室血液不充盈,结果使动作电位持续时间减少,原因是所有的离子过程还未完全恢复到静息态,减少的跨膜电流会导致第二个动作电位的持续时间比第一个动作电位短[15].心肌细胞的动作电位时长(action potential duration,APD)随上一个周期的静息时间长度(也称舒张间隔,diastolic duration,DI)变化关系(曲线)通常称为APD 恢复关系(曲线)[16-19],当舒张间隔很长时对应的动作电位时长本文称为稳态动作电位时长(即稳态APD).由于动作电位时长恢复曲线的变化在生物体的正常节律调节中起着有重要作用,它与心律失常、螺旋波破碎机制的关系受到了研究者的极大关注[20-29].实验和理论研究发现,APD 恢复曲线倾斜度的不均匀分布对心律失常、螺旋波失稳起着重要作用[20-23],陡峭的APD 恢复曲线容易引起室颤[24,25],通过药物平缓APD 恢复曲线的斜率能达到缓解心律失常的治疗效果[25,26];倾斜度大于1 的APD 恢复曲线使得细胞动作电位在高频刺激下容易诱发出持续的APD 长短交替,这时很小的舒张间隔变化都会导致APD 的剧烈变化,进而引起细胞不应期的变化[27,28],还会导致螺旋波破碎[29].

2005 年,Clayton 和Taggart[30]发现斜率大于1 的APD 恢复曲线不是引发螺旋波破碎的必要条件,空间上具有两个恢复曲线斜率都小于1 的区域且这两个区内元胞连续两次激发产生的APD 差别很大时也会导致螺旋波破碎;2006 年,Nash 等[31]测量了14 名患者的整个心室外膜组织上256 个点的动作电位恢复曲线斜率,观察到动作电位恢复曲线斜率分布极不均匀,其中54%的记录点处的斜率分布在[0,1)范围内,20%的记录点的斜率分布在[1,1.5)范围内,约11%的记录点的斜率分布在[1.5,2)范围内,其余记录点斜率分布在[2,4.1]范围内,他们还通过数值模拟观察到螺旋波在恢复曲线斜率大于1 情况下既可能破碎、也可能不破碎,但是未能解释其原因.从上述研究可以看出,APD恢复曲线的倾斜度对螺旋波演化行为的影响尚未完全清楚,需进一步探讨,而APD 恢复关系中的稳态APD 对螺旋波的影响仍缺乏研究.

元胞自动机是研究复杂系统的重要辅助手段,具有物理图像清晰、无截断误差、完全的并行性等优点,在心肌组织螺旋波动力学和电生理特性建模的研究中得到广泛应用[32-35],先后用于研究非均匀心肌组织中螺旋波自发产生和终止以及靶波到时空混沌态的转变[32];正常和心肌缺血情况下的心电图走势[33,34];心房颤动的随机产生和终止的机制[35],这些研究丰富了心肌组织中螺旋波动力学的认知并拓宽元胞自动机的应用.本文将心肌细胞APD恢复关系引入到元胞自动机模型中,研究不同倾斜程度的APD 恢复曲线和稳态APD 对螺旋波演化行为的影响,发现APD 恢复曲线斜率大于1 并不一定导致螺旋波破碎,也可使螺旋波漫游,甚至还可以使螺旋波保持稳定,对产生这些现象的物理机制做了解释,并根据数值模拟结果提出消除螺旋波的方法.

2 考虑APD 恢复关系的元胞自动机模型

取具有500×500 个格点的二维正方形点阵,可激发心肌细胞均匀分布在这个点阵的格点上,称为元胞,相邻元胞间距为1.元胞在正方形点阵上的位置用坐标 (i,j) 表示,其t时步的状态用ui,j(t) 表示,对应心肌细胞的膜电位,采用扩展型Moore 邻居,元胞邻域半径为r(即在邻域半径内有(2r+1)2-1个邻居),元胞状态总数为n+1,即元胞 (i,j) 于t时刻的状态ui,j(t) 在集合{0,1,2,···,n(n≥2)}中取值,其中ui,j(t)=0 表示静息态,ui,j(t)=1 表示激发态,ui,j(t)=2,3,···,n表示不应态.不同状态之间的演化规则为

式中Mi,j代表t时步以静息态元胞 (i,j) 为中心、半径为r的邻域范围内处于激发态邻居元胞的个数,K为元胞激发阈值,n称为元胞的APD,初始元胞的APD 记作n0.在元胞自动机模型中引入文献[17,31,36]给出的APD 恢复关系:

式中α就是稳态APD,α,β,γ为常数.不同DI 对应的恢复曲线斜率定义为

为了对(2)式的APD 恢复曲线有直观印象,图1给出了α=40 时不同β,γ下的APD 恢复曲线和斜率λ随舒张间隔DI 的变化.从图1 可以看出,不同参数的曲线在DI 较小情况下的倾斜程度不同(即最大斜率不同),在γ不变的情况下增大β,曲线的倾斜度增大,最大斜率相应增大;对于同一条APD 恢复曲线,不同DI 值对应不同的斜率.所有曲线有相同的特点:曲线斜率随着DI 增加而单调减少,当 D I ≥8 个时步时,曲线的斜率小于1;当DI ≥30个时步时,曲线斜率趋于0,APD 达到稳态APD,γ越小越容易达到稳态APD.在数值模拟中,α,β和γ都是可调参数,APD 对应元胞自动机模型中的参数n;DI 为元胞在相继两次激发之间处于静息态的时长,简称静息时长,其影响着元胞在下一激发周期的APD.数值模拟时我们始终保持r=3,K=6 不变,采用映射型边界条件,即分别以系统的4 条边界为镜面向外拓展宽为r=3的虚拟区域,将系统内距边界小于3 的元胞的状态镜面反射到边界外,并以顶点为对称点将系统4 个边角区域的元胞状态映射出去.在执行映射型边界条件后,边界处的元胞就与系统内其余元胞具有相同的邻居个数.

图1 在 α=40 和不同 β,γ 下APD 恢复曲线(a)和斜率 λ 随舒张间隔DI 的变化(b)Fig.1.APD restitution curve (a) and the change of its slope λ with diastolic interval (DI) (b) for α=40 and different values of β and γ .

3 在不同倾斜度的APD 恢复曲线下的数值模拟结果

在不考虑APD 恢复关系(A PD=n0=40,r=3,K=6)情况下产生螺旋波初态,其是漫游螺旋波,波头轨迹接近长条型矩形闭合回路,长度约40 个格点,绝大部分元胞的静息时长DI=1 时步.然后在系统中考虑APD 恢复关系并以此时为计时起点,研究不同倾斜程度的APD 恢复曲线和稳态APD 对螺旋波动力学的影响.

为研究在不同倾斜度的APD 恢复曲线下的螺旋波演化行为,保持α=40,γ=10 不变,通过调节β来改变APD 恢复曲线的倾斜度,观察不同β下初态螺旋波的演化行为.从图1 可以看出,引入APD 恢复关系后,元胞的APD 最初从n0=40 将减小到DI=1 对应的最小值,之后元胞的APD动态变化.数值模拟发现,当β≤10 时,初态漫游螺旋波演化成稳定螺旋波,说明APD 恢复关系可以使不稳定螺旋波变得稳定;当 11 ≤β≤14 时,螺旋波发生漫游,且β越大,螺旋波漫游的幅度越大,漫游情况如图2 所示,图2 中白色对应元胞处于静息态,黑色线对应元胞处于激发态,灰色区对应元胞处于不应态,黑色线与白色区的边界线之间的灰色区宽度除以每时步波传播的距离对应元胞的APD,白色条纹越宽代表DI 值越大.从图2 可以看出,螺旋波漫游导致DI,APD 分布略有不均匀,但没有引起传导障碍,所以螺旋波波臂没有折断;当β≥15 时,螺旋波发生波臂断裂和破碎,其破碎方式有三种.

图2 在参数 α=40 ,γ=10 ,β=14 下不同时刻的膜电位斑图 (a) t=0 ;(b) t=20000 ;(c) t=60000 ;(d)t=100000Fig.2.Patterns of the membrane potential at different time moments for α=40 ,γ=10 and β=14 :(a) t=0 ;(b) t=20000 ;(c) t=60000 ;(d) t=100000 .

第一种破碎方式是:螺旋波先发生漫游,然后波臂在波头附近断裂,最终导致螺旋波完全破碎,这种由于漫游导致的螺旋波破碎称为多普勒失稳[37],如图3 所示.在图3 中,为了同时显示处于不应态的元胞APD,在作膜电位斑图时,将元胞APD 在[31,40]范围的不应态元胞统一用红色显示(记为A 态),元胞APD＜31 的不应态元胞用绿色显示(记为B 态),激发态和静息态元胞分别用黑色和白色显示,之后的斑图也以相同方式显示.从图3 可以看出,引入APD 恢复关系后,白色线逐渐变宽,因为元胞的APD 瞬间减小到n ≈26 (参见图1),致使元胞在演化过程中提前返回到静息态.一方面由于元胞APD 减小,使波头加快旋转而漫游,白色线逐渐演化成粗细极度不均匀的线(DI 分布极度不均匀);另一方面,波头加快旋转使螺旋波的中心区黑色线曲率增大,致使大部分元胞都处于B 态,而处于A 态的元胞则呈零散分布,即在螺旋波波头附近,A 态元胞与激发态相邻且呈小的辐射状分布,导致传导障碍使波臂折断.而在螺旋波外围,位于波后的A 态元胞呈小块状分布,不与激发元胞相遇,因此不引起传导障碍(参见图3(c),(d)).随着螺旋波波臂多次折断,所有A 态元胞都呈辐射状分布,最终使螺旋波完全破碎(参见图3(e)—(h)),可见多普勒失稳是螺旋波波头漫游使其附近大量A 态元胞分布呈辐射状导致传导障碍造成的.

图3 在参数 α=40 ,γ=10,β=15 下不同时刻的膜电位斑图,白色代表静息态;黑色代表激发态;绿色(红色)代表元胞APD在[26,30] ([31,40])范围的不应态 (a) t=0 ;(b) t=400 ;(c) t=490 ;(d) t=510 ;(e) t=560 ;(f) t=1800 ;(g) t=2000 ;(h)t=3070Fig.3.Patterns of the membrane potential at different time moments for α=40 ,γ=10 and β=15 .The white and black represent the rest-state and excited state,respectively,the green (red) represents the refractory states in which the APD of the cell is in the range of [26,30] ([31,40]):(a) t=0 ;(b) t=400 ;(c) t=490 ;(d) t=510 ;(e) t=560 ;(f) t=1800 ;(g) t=2000 ;(h) t=3070 .

第二种破碎方式是:螺旋波波臂先在外围折断,最终导致螺旋波完全破碎,这种螺旋波破碎方式称为爱克豪斯失稳[37],如图4 所示.从图4 可以看出,引入APD 恢复关系后,白色线逐渐变宽(参见图4(b)),是因为元胞的APD 瞬间减小到n ≈23(参见图1),元胞的APD 变化相比于图3 的情况进一步增大.从图4(b)—(e)可以看出,在元胞APD变小后,螺旋波波臂保持均匀,说明螺旋波是稳定螺旋波,正因为是稳定螺旋波,使得由处于B 态元胞组成的螺旋波首先在其外围出现呈小三角形分布的A 态元胞,导致螺旋波波臂折断,随后在折断处形成新的螺旋波(见图4(f)).随着螺旋波继续演化,在t=310 时步,螺旋波波头周围才出现A 态元胞大范围呈小块状分布或小块辐射状分布,这种分布不会导致波臂折断,因螺旋波不漫游,螺旋波中心区保持完好.螺旋波最终破碎是在外围形成的螺旋波与原螺旋波相互碰撞的结果,系统最终演化为多螺旋波态(参见图4(h)).

图4 在参数 α=40 ,γ=10,β=18 下不同时刻的膜电位斑图,作图方式与图3 一致 (a) t=0 ;(b) t=40 ;(c) t=220 ;(d) t=250 ;(e) t=255 ;(f) t=310 ;(g) t=2850 ;(h)t=34300Fig.4.Patterns of the membrane potential at different time moments for α=40 ,γ=10 and β=18,the drawing method is consistent with Fig.3:(a) t=0 ;(b) t=40 ;(c) t=220 ;(d) t=250 ;(e) t=255 ;(f) t=310 ;(g) t=2850 ;(h) t=34300 .

第三种破碎方式是:APD 和DI 大幅度振荡变化导致螺旋波破碎成时空混沌,如图5 所示.从图5可以看出,引入APD 恢复关系后,元胞APD 从n0=40瞬间减小到n ≈10 (根据方程(2)可以得出),白色线逐渐变得很宽(螺旋波波头附近除外),直到所有元胞处于B 态.在下一次激发时,除螺旋波中心区外,其他位置的元胞又恢复至A 态,形成螺旋状分布,在螺旋波波头附近与B 态元胞分布区域相遇(参见图5(b)),导致传导障碍使螺旋波波头消失,出现新的断点,此时螺旋波不再旋转,而是断点回缩(参见图5(c)),直到传导障碍消失后,螺旋波又开始旋转(参见图5(d)),从效果上看,螺旋波类似作大幅度漫游.从图5(d),(e)可以看出,由于A 态元胞和B 态元胞分布区在波头附近相遇,导致螺旋波波臂在波头附近折断,脱离主体这段小波演化成由处于B 态的元胞组成的小螺旋波.随着系统演化,A 态和B 态元胞分布区交替出现,且总在螺旋波波头附近相遇,产生传导障碍,最终使螺旋波破碎成时空混沌.我们把元胞APD 大幅度振荡变化称为APD 交替,可认为螺旋波破碎是由于APD 按A→B→A→B 方式交替(记为AB 交替)变化导致的,APD 交替变化在心律失常中是常见的现象[27,28].

图5 在参数 α=40 ,γ=10 ,β=33 下不同时刻的膜电位斑图,作图方式与图3 一致 (a) t=0 ;(b) t=50 ;(c) t=60 ;(d)t=100;(e) t=130 ;(f) t=150 ;(g) t=200 ;(h)t=500Fig.5.Patterns of the membrane potential at different time moments for α=40 ,γ=10 and β=33,the drawing method is consistent with Fig.3:(a) t=0 ;(b) t=50 ;(c) t=60 ;(d) t=100 ;(e) t=130 ;(f) t=150 ;(g) t=200 ;(h) t=500 .

以上结果表明,在不同β值下,螺旋波破碎方式存在明显的不同:第一、第二种螺旋波破碎不存在大范围APD 交替变化,只是局部存在APD 非均匀分布导致传导障碍使螺旋波破碎,由于形成的APD 非均匀分布具有偶然性,因此螺旋波完全破碎需要经过很长时间.第三种螺旋波破碎由于APD 按AB 方式交替导致,螺旋波完全破碎需要的时间很短,这种破碎方式可以使心脏由心动过速快速转变为心室纤维性颤动,危害很大.

为了解参数β,γ同时变化对螺旋波的影响,图6给出了在β-γ参数平面上的相图.为了解释图6的结果,在给定参数β,γ下计算总平均斜率,该值用如下方式得到:首先将t时步各激发元胞的DI 值求和后求平均,得到t时步的平均DI,此时可通过(3)式算出平均斜率,接着将 1-t0(t0为螺旋波完全破碎前的某个时刻)这段时间内各时刻平均DI 求和再平均,得到总平均DI,最后利用(3)式算出作为该组参数的总平均斜率.对于不破碎螺旋波,取t0=2000 时步,如果螺旋波很快漫游出系统,则取螺旋波消失前某个时刻为t0值,图7 给出了随β,γ变化的直方图.

图6 在 β-γ 参数平面上的相图Fig.6.Phase diagram on β-γ parameter plane.

从图6 和图7 可以看出:760 组参数中有216 组参数出现稳定螺旋波且集中分布在一个大的区域,稳定螺旋波对应参数的总平均斜率一般在[0.0476,0.8671]范围内变化,有2 组参数的总平均斜率超过1,这两组参数是 (β,γ)=(3,1),(4,1),表明总平均斜率大于1 也能使螺旋波稳定.有85 组参数出现螺旋波漫游且几乎分布在一个大的区域,漫游螺旋波的总平均斜率在[0.5297,2.2073]范围内变化,其中有4 组参数的总平均斜率高于1,这四组参数是 (β,γ)=(4,2),(5,1),(5,2),(6,1),可见总平均斜率大于1 不一定使螺旋波破碎.当β≥18时,出现第一种、第三种螺旋波破碎和螺旋波漫游出边界消失,不存在稳定和漫游螺旋波,其中有167 组参数出现第一种螺旋波破碎方式,这种破碎方式的总平均斜率在[1.1198,1.2754]范围内变化,有98 组参数出现第三种螺旋波破碎方式,这种破碎方式的总平均斜率一般在[1.0014,1.0233]范围内变化,有4 组参数的总平均斜率小于1,这四组参数是 (β,γ)=(19,19),(19,20),(22,18),(22,19),说明总平均斜率小于1 也可以使螺旋波破碎;从图6 可以看出,对于第一、第三种螺旋波破碎方式,其参数混合地分布在γ≥10 的区域,表明这两种破碎方式具有偶然性,对参数变化比较敏感.有168 组参数出现螺旋波漫游出边界消失且几乎分布在一个大的区域(对应γ ＜10),螺旋波消失对应的总平均斜率约在[1.2738,2.9874]范围内变化,这表明要让螺旋波漫游出系统消失,需要的恢复曲线斜率非常大.有26 组参数出现第二种螺旋波破碎方式,而且主要以孤立方式出现,这种破碎方式的总平均斜率在[1.0301,1.1196]范围内变化,可见第二种螺旋波破碎少见.

上述结果表明:用恢复曲线总平均斜率是否大于1 来判断螺旋波的演化行为是否稳定和漫游存在不可靠因素,为了解释其原因,图8 给出了几组螺旋波稳定、漫游、破碎和消失对应的参数下λt随时间的变化.从图8(a),(b)可以看出,即使总平均斜率大于1,只要平均斜率保持不变(即元胞的APD和DI 保持不变),螺旋波也可以稳定,同理对漫游螺旋波也成立(参见图8(d)),这个结果说明要消除心律失常应尽可能使心肌细胞的APD 和DI 保持不变.图8(a),(b)为总平均斜率大于1 的特殊情况,一般情况下稳定螺旋波对应的平均斜率变化如图8(c)所示.

从图8(e)可以看出,即使总平均斜率大于1,只要平均斜率随时间大幅度快速振荡,也可以保证螺旋波不破碎而处于漫游态.对于一般漫游螺旋波,其总平均斜率小于1,平均斜率随时间在某个固定值附近小幅振荡(见图8(f)).

图8 不同参数(螺旋波态)下恢复曲线的平均斜率随时间变化 (a) β=3 ,γ=1 (稳定);(b) β=4 ,γ=1 (稳定);(c) β=7,γ=17(稳 定);(d) β=4 ,γ=2 (漫 游);(e) β=5 ,γ=1 (漫 游);(f) β=10 ,γ=9 (漫 游);(g) β=13,γ=5 (第一种破碎);(h) β=18,γ=11 (第二种破碎);(i) β=38,γ=15 (第三种破碎);(j) β=19 ,γ=19 (总平均斜率小于1 的破碎);(k) β=31,γ=5(消失);(l) β=32 ,γ=8 (消失)Fig.8.Evolution of the average slope of the restitution curve for different parameters (spiral wave states):(a) β=3 ,γ=1(stable);(b) β=4 ,γ=1 (stable);(c) β=7 ,γ=17 (stable);(d) β=4 ,γ=2 (meandering);(e) β=5,γ=1 (meandering);(f) β=10,γ=9 (meandering);(g) β=13,γ=5 (first type of breakup);(h) β=18 ,γ=11 (second type of breakup);(i) β=38,γ=15(third type of breakup);(j) β=19 ,γ=19 (breakup with total average slope less than 1);(k) β=31 ,γ=5 (disappear);(l) β=32,γ=8 (disappear).

对于破碎的螺旋波,其对应平均斜率随时间变化一般在某个大于或等于1 的值附近小幅振荡(见图8(g)—(i)),但存在平均斜率在小于1 的范围内振荡的情况(参见图8(j)),其总平均斜率也小于1.为了对总平均斜率小于1 下螺旋波的破碎有直观印象,图9 给出了在图8(j)参数下不同时刻的膜电位斑图.从图9 可以看出,引入APD 恢复关系后,白色线没有加宽,而是出现元胞APD 相差很大的两个螺旋状区域(A,B 状态区),这两块区通过细的白色线分离开(参见图9(b),(c)),但在波头附近相遇,导致传导障碍使螺旋波波臂在波头附近折断,脱离主体的这段小波演化成螺旋波(参见图9(c),(d)),而原螺旋波停止旋转,直到波头附近恢复为绿色区.由于两块元胞APD 相差很大的螺旋状区域大小交替变化,且都在螺旋波波头附近相遇,致使螺旋波波臂在波头附近陆续出现折断现象,脱离主体的这段小波同样演化成螺旋波,这些螺旋波相互作用,最终导致螺旋波破碎.与图3 所示的第一种破碎不同的是,图9 显示的这种破碎在初始阶段螺旋波不发生漫游.与图5 所示的动作电位按AB 交替变化导致的第三种破碎不同的是,图9 显示动作电位按A→AB→B→BA→A 交替变化.存在两条白色线,第一条白色线位于激发态波前,第二条白色线位于A,B 状态分布区之间,如图9(b)所示,其实际上是t=0 时刻的静息态线所在位置,它是这样产生的:在波传播的过程中,激发波线前方的静息态元胞相继变成B 态,而最先变成B 态的元胞返回到静息态时,有与之相邻元胞仍处于A 态,导致第二条白色线出现,它是由于AB 两种状态的元胞APD 相差一倍多引起的.当然,如果分别计算螺旋波波头周围和螺旋波外围的APD 恢复曲线总平均斜率,可以发现:在螺旋波中心区的总平均斜率大于1,因为波头周围元胞处于B 态,APD 较小,而静息时长是1 个时步;在螺旋波其他区的总平均斜率接近0,使得整体总平均斜率比1 小很多,这就是为什么A,B 两种元胞态分布区总在波头附近相遇导致传导障碍.

图9 在参数 α=40 ,γ=19 ,β=19 下不同时刻的膜电位斑图,作图方式与图3 相同 (a) t=0 ;(b) t=27 ;(c) t=37 ;(d) t=40 ;(e) t=57 ;(f) t=80 ;(g) t=90 ;(h) t=113 ;(i)t=1900Fig.9.Patterns of the membrane potential at different time moments for α=40 ,γ=19 and β=19 .The drawing method is consistent with Fig.3:(a) t=0 ;(b) t=27 ;(c) t=37 ;(d) t=40 ;(e) t=57 ;(f) t=80 ;(g) t=90 ;(h) t=113 ;(i) t=1900 .

对于消失的螺旋波,其参数的平均斜率随时间大幅度慢振荡,最大平均斜率可达到5,最小接近0,如图8(k),(l)所示,正是这种大幅度慢振荡,螺旋波出现呼吸效应(即APD 大幅度交替变化),从而引起螺旋波波头间歇性遇到传导障碍而作漂移运动,其效果表现为螺旋波在短时间内大幅度漫游,迫使螺旋波运动出边界消失,因此可以通过使心肌细胞动作电位大幅度交替变化来消除螺旋波.

上述研究结果表明,螺旋波的演化与恢复曲线倾斜度有关,螺旋波破碎与恢复曲线总平均斜率大于或等于1 有密切联系,但是恢复曲线总平均斜率大于1 并不一定导致螺旋波破碎,只要恢复曲线平均斜率随时间大幅度高频振荡或者保持不变,也可以只导致螺旋波漫游,甚至稳定;此外,恢复曲线总平均斜率小于1,螺旋波也不一定保持不破碎态,只要出现元胞的APD 大幅度交替变化,A 态和B 态元胞分布区在螺旋波波头附近相遇,波在传播过程中就会因遇到不应态元胞而导致传导障碍,诱发螺旋波波臂折断,这种现象频繁出现就会导致螺旋波破碎.

4 在不同稳态APD 恢复曲线下的数值模拟结果

前面只探究了在不同倾斜度的APD 恢复曲线下的螺旋波演化行为,为了探究不同稳态APD 对螺旋波的影响,首先取不同n0来产生初态螺旋波,其可以是稳定或漫游螺旋波,然后在固定值γ=9及不同参数α和β情况下观察螺旋波的演化行为,图10 给出了在β-α参量平面上的相图.从图10可以看出,在给定β的情况下,逐渐增大稳态APD(即α),螺旋波出现从破碎→漫游→消失的转变过程,而且都有这样的特点:1) 当α比n0小足够多时,螺旋波会按第一种方式破碎,当α比n0大足够多时,螺旋波会消失,这是由于引入APD 恢复关系后,元胞APD 剧增,元胞返回到静息态的时间延长,导致整体出现传导障碍,因此可以通过增大心肌细胞稳态动作电位来消除心脏中的螺旋波;2)当n0=25,40,55 时,126 组参数中分别有19,44,69组参数处于漫游螺旋波态,比例依次为0.1508,0.3492,0.5476,说明增大初态螺旋波各元胞的APD,有利于避免螺旋波破碎.

图10 在不同 n0 情况下 β-α 参数平面上的相图 (a) n0=25 ;(b) n0=40 ;(c) n0=55 .Fig.10.Phase diagram on β-α parameter plane for different values of n0 :(a) n0=25 ;(b) n0=40 ;(c) n0=55 .

5 结论

本文采用二维可激发介质元胞自动机模型,研究了元胞APD 按不同的APD 恢复曲线动态变化后对螺旋波的影响,探讨螺旋波失稳的机制,研究发现:不同倾斜度的APD 恢复曲线可以使不稳定螺旋波变得稳定,也可以使螺旋波漫游、破碎和消失.对于稳定和漫游螺旋波,对应的APD 恢复曲线总平均斜率一般小于1,但也存在总平均斜率大于1 的情况下仍不出现传导障碍,是因为APD 恢复曲线平均斜率保持不变(即元胞APD 和DI 保持不变),或者平均斜率大幅度高频振荡.螺旋波破碎通常与APD 恢复曲线总平均斜率大于1 有关,但是少数情况下出现APD 恢复曲线总平均斜率比1 小很多时也会导致螺旋波的破碎,原因是复杂的A 态、B 态交替引发传导障碍而使螺旋波破碎,如果A 态、B 态交替使螺旋波波头消失,新出现的断点只作漂移运动,直到传导障碍消失才继续旋转,这种情况会导致螺旋波漫游出系统消失,其与APD 恢复曲线的平均斜率大幅度慢振荡有关.如果初态元胞APD 比APD 恢复曲线的稳态APD小很多时,还会导致螺旋波直接消失,这时恢复曲线平均斜率不再振荡,因此增大稳态APD,有利于螺旋波稳定不破碎和消失,但是初态APD 比APD恢复曲线的稳态APD 大很多时,可导致螺旋波破碎.观察到三种螺旋波破碎方式,即多普勒失稳、爱克豪斯失稳、APD 交替变化导致螺旋波破碎,其中多普勒失稳和APD 交替变化导致的螺旋波破碎所占比率最大,爱克豪斯失稳占比很小,这些现象在心脏系统中都观察到了,我们的结果有助于理解心脏系统中螺旋波的演化机制,为心律失常的治疗和螺旋波的控制提供有用信息.