摘要：应用功率键合图方法，建立了一种多分支血液循环系统的计算机仿真模型，即描述血液循环系统内血流动力学变量变化规律的状态方程。该仿真模型较为细致地刻画了血液循环系统的生理特性，形成了较完整的人体血液循环系统的计算机模型，此模型可模拟血液循环系统的生理和病理特性，得出相应的心血管动力学仿真数据和波形，为进行血液循环系统生理和病理的医学研究提供了新的研究手段。

关键词：血液循环系统　计算机仿真　功率键合图法

0  引　言

　　功率键合图法是一种系统动力学建模方法，它以图形方法来表示、描述系统动态结构，是对流体系统进行动态数字仿真时有效的建模工具。通过已有的研究工作表明，功率键合图方法可以较好地应用于生物流体系统仿真，特别是人体循环系统的建模和数字仿真^[10]。

          我们在以前的工作当中，建立了一个简化的血液循环系统模型^[10]，验证了功率键合图法的可行性和有效性。键合图建模方法的优点是直观形象，便于获得状态空间方程，有利于数值化计算，避免了电模拟方法中推导状态方程困难的弱点 。本文对血液循环系统进行了较细致和全面的划分，建立了一个包括动脉系统、静脉系统、心脏（左、右心室和心房）以及冠脉循环、外周循环的多分支血液循环系统仿真模型。

　　应用功率键合图方法对血液循环系统进行建模和仿真的基本规则是，(1)把血液循环系统的结构及各主要动态影响因素以图示模型形式，即功率键合图加以表示，(2)从功率键合图出发，建立系统的动态数学模型——状态空间方程，(3)在数字计算机上对状态方程进行求解。

1  多分支血液循环系统模型的建立

1.1 系统描述

　　血液循环系统模型如图１所示^[4]。在心血管循环系统中，血液在心脏“泵”的作用下所进行的循环流动，可以看作是一种功率流的流动、传输、分配和转换的过程。血液在左右心室有节律地收缩作用下，被泵向人体的各个部分，其中包括：体循环区（血液由左心室经主动脉、大动脉、外周循环区和腔静脉，回到右心房），肺循环区（血液由右心室流经肺动脉和肺静脉到左心房。），腹部内循环，颈部和头部循环，以及冠脉循环等。在心房和心室、心室和主动脉之间存在着防止血液倒流的膜瓣，如二尖瓣、三尖瓣、主动脉瓣等。

 图1　血液循环系统模型

1.2  功率键合图模型

　　应用功率键合图建模方法的第一步是将原系统表达为功率键合图的图示模型。功率键合图由功率键、结点和作用元等主要元素构成，多分支血液循环系统的功率键合图如图2所示。

 R_nv     C_hv     R_hh     C_ha      R_na

图２　多分支血液循环系统功率键合图模型(此图有省略) 　

    参考图2，绘制多分支血液循环系统功率键合图的步骤可简述如下：

    (1)根据对多分支循环系统各个功率流程分支的分析，依次确定各0结点和1结点。

    ０结点表示集总的流容容腔，如心室腔、主动脉弹性腔，在0结点处血液压力为等值，而该结点输入的血流量等于输出的血流量。１结点表示集总的流阻管路或流感管路，如大动脉血管，在1结点处血流量为等值，而该结点的压力降等于上流压力值减去下流压力值。在图2 的循环系统模型中共有15个0结点和21个1结点。

    (2)画上各结点周围的功率键，并标注功率流向。

    功率键是带有箭头和因果线表示功率的线段。本模型中构成功率的两个变量是血压和血流。箭头表示系统作用元中的功率流向，即循环血液的流动方向。

    (3)在功率键的一端标注上相应的C、R、L作用元。

　　为了能够全面、细致地刻画系统特性，本模型中应用了三种作用元：流容、流阻和流感。

    流容反映血管的顺应性，画在０结点上，用C来表示，简称C元。例如，图2 中的C_ta、C_ar、C_vn、C_pa、C_pv是分别表示与图1相对应部分的胸主动脉、大动脉、腔静脉、肺动脉和肺静脉顺应性的流容。

    流感反映血流的惯性特性，画在1结点上，用L来表示，简称L元。如图2中的L_ta、L_ar、L_vn、L_pa、L_pv、L_co是分别表示相对应的胸主动脉、大动脉、腔静脉、肺动脉、肺静脉及冠状动脉血流惯性的流感。

    流阻反映血流粘滞阻力的特性，简称R元，画在１结点上。例如图2中R_ta、R_ar、R_vn、R_pa、R_pv和R_co是分别表示胸主动脉、大动脉、腔静脉、肺动脉、肺静脉及冠状动脉血流粘滞阻力的阻性作用元。

    (4)在各功率键上标注因果线，以便于建立系统的数学模型。

    功率键上的因果线表示各作用元上流量与压力两变量之间的因果关系，确定了自变量和因变量，便于建立系统的状态方程。对于C元，其功率键上两个变量间，自变量是流量，因变量是压力；对于L元和R元，其功率键上两个变量间压力是自变量，流量是因变量。

    经过以上步骤，就完成了循环系统的功率键合图模型。可以看出，键合图模型就是通过结点、功率键和作用元这些元素对心血管循环系统直观而形象的描述和反映。在将循环系统翻译成键合图模型后，就可以方便、有条不紊地推导系统数学模型。

2  系统数学模型

    功率键合图建模方法的第二步是推导系统的数学模型。在推导系统动态过程的数学模型——状态方程时，首先要确定状态变量。应用键合图方法建模的方便之处就在于对状态变量的确定有一定之规，可遵循固定的法则。

    由于系统的状态方程是一阶微分方程组，在其变量间有导数关系，而在键合图中，只有流容C和流感L作用元中的两个变量间才有导数或积分关系，所以应当从C元和L元各自的变量间取一个变量作为状态变量。

    对于C元，自变量为流量，因变量为压力，其关系为：

　　         　　　　　　(1)

    对于L元，自变量为压力，因变量为流量，其关系为：

                                               (2)

    对于R元，流量和压力之间的关系有：

                                                     (3)

    根据规则，取C元功率键上的压力变量p和L元功率键上的流量变量Q为状态变量,状态变量的一阶导数即为状态方程。

    因此，对于0结点，由(1)式两边取导数可得：

                             (4)

其中， 是第i个0结点处的压力， 为输入血流量， 为输出血流量， 是第i个0结点处的流容。

    对于1结点，由(2)式和(3)式可得：

                         (5)

其中， 是第i个１结点处的血流量， 为上流压力， 为下流压力， 和 分别是第i个１结点处的流阻和流感。

　　对每个0节点和1结点都建立类似(4)和(5)的关系式，则可以得到系统的数学模型。本模型的数学模型是36阶的状态空间方程，即模型由36个一阶微分方程组成。下面列出了主动脉循环部分的状态方程：

                                      (6)

                                               (7)

                                                          (8)

                                  (9)

                                (10)

    (11)

    其中，C_ta、C_aa、C_ar分别是胸主动脉、腹主动脉、外周动脉的流容；L_ta、L_aa、L_ar、L_vn分别是胸主动脉、腹主动脉、外周动脉和腔静脉的流感；R_ta、R_aa、R_sa、R_pc和R_sv是分别表示胸主动脉、腹主动脉、外周动脉、外周循环和腔静脉的流阻。p_tao、p_aao、p_sar和Q_tao、Q_aao、Q_sar分别是动脉循环中的胸主动脉、腹主动脉、外周动脉部分的压力和流量。

　　血液循环是由心脏的舒张－收缩动作推动的，本文采用了心室时变流容 来表示这种舒张－收缩动作， 是时间的周期函数。

    对于循环系统中的膜瓣作用，可以作为模型的约束条件加入到系统数学模型当中：当血液正向流动时，膜瓣阻力为一较小的数值；当血液反向流动时，膜瓣阻力为无穷大，即阻止血液倒流。

　　本模型中的流容、流阻和流感参数参照文献[4]。

3  计算机仿真

　　本文采用４阶定步长Runge－Kutta法来求解模型的状态方程，设定仿真步长为0.0001s，在奔腾586 PC机上进行数字仿真。

    当加入边界约束条件，设置各状态变量初始参数之后，状态变量便以状态方程为基础被同步地展开。在每一步，血液循环系统各部分的压力和流量值根据状态方程被分别计算出来。待仿真数据变化稳定后，由系统输出方程可以得到每个心动周期内系统各部分的血压p、血流量Q、血液容量V以及心输出量CO和射血分数EF等各项生理参数数值，从而可以对多项生理特性进行计算机仿真。本文进行了正常生理条件下和高血压、血管刚性的病理条件下的生理特性仿真。

3.1 正常生理状态仿真

  　设定各状态变量的初始参数为正常值^[4,5]，对系统模型进行计算，即可得到正常生理条件下，血液循环系统血流动力学参数的仿真数据。

    图３给出了在正常状态时，三个心动周期（每个心动周期为0.8秒）内的左心室压力和主动脉血的仿真波形压的仿真波形。从压力仿真波形图中可以看出，心室压力和主动脉压力在每个心动周期内的压力脉动是十分显著的。图４是肺动脉血压和肺静脉血压的仿真波形。肺动脉压的压力脉动也较为显著，而在肺静脉中，血液的压力脉动就不很明显。

 
图3 左心室和主动脉的压力变化仿真

140                               
01.6
t/s
(a)左心室血液容量的周期变化 

140                               
01.6
t/s
 (b)右心室血液容量的周期变化  

图4 肺动脉和肺静脉的压力变化仿真

    在表1中给出了血液循环系统主要血流动力学变量在正常状态时条件下的仿真数值。由生理学规律可知 ，左心室收缩压范围一般在17~18 kPa，主动脉压力范围在12~17 kPa，肺动脉压在2 kPa左右。因此，仿真所得波形和数据与实际的生理规律是相符的。

    表1中还给出了评定心脏功能的两个有用的指标：心输出量CO和射血分数EF，仿真所得到的数据为：心输出量5256 ml/min，射血分数61%，都符合实际的生理规律 。

                            表1 血液循环系统主要血流动力学变量计算机仿真数值             

3.2 高血压仿真

    由于动脉管径窄缩，或是动脉壁增厚等原因常常会使动脉血管的阻力增大，使得心脏在收缩期向主动脉喷血时耗费更多的功，从而引起高血压症状。因此在本实验中，增大键合图模型中的主动脉和外周动脉的流阻R_ta、R_aa、R_ar的数值，可以实现高血压的仿真。

    表1中给出了高血压时各血流动力学变量的仿真数据。从仿真数据中可以看到，左心室压和主动脉压分别达到21.28 kPa和18.63 kPa，血压值明显升高，但是心输出量4989 ml/min和射血分数54%的数值却比正常状态显著降低，这表明高血压时心脏的功能在减弱。

3.3 血管刚性仿真 

    血管顺应性的倒数1/C被称为血管刚性，血管刚性越大，血管顺应性则降低，使心室射血阻抗增大，导致心室喷射压力和动脉血压升高，心输出量和射血分数降低。在本实验中，将主动脉与外周血管的流容C_ta、C_aa、C