技术领域：
本发明属于血流动力学数值模拟领域。尤其涉及一种基于零维血流动力学模型的ffr快速计算系统。
背景技术：
：
冠脉的狭窄是导致心肌缺血从而引发心血管疾病的主要原因。目前临床诊断功能性心肌缺血的“金标准”是血流储备分数(ffr)。它被定义为冠脉狭窄时最大血流量与同一支冠脉无狭窄时最大血流量的比值。基于在心肌微循环阻力恒定且小到忽略不计时的心肌血流量与灌注压成正比的概念，可被化简为最大充血状态下狭窄远端的平均压力与主动脉根部的平均压力的比值。临床上通过侵入式导管插入术来直接测量压力。但是，这些侵入式测量技术存在一定的风险，且费用高。因此无创计算ffr对临床应用具有重大意义。
目前常用的无创计算ffr的方法是基于cta图像，结合个性化冠脉零维血流动力学模型，构建0d-3d耦合的几何多尺度血流动力学模型，采用计算流体动力学方法计算静息和最大充血状态下冠脉内的流场和压力，实现ffr的无创计算。但是这种计算方法由于需要零维模型和三维模型进行耦合迭代计算，模型复杂，计算量大，仿真时间长，因此提出一种快速计算ffr的方法很有必要。
零维血流动力学模型由于参数少、生理意义明确，计算简单而常被用来模拟人体心血管循环系统。最早的零维血流动力学模型是frank(1899)基于弹性腔理论提出的二元件模型。弹性腔理论把大动脉血管类比为一个弹性腔，将小动脉血管及毛细血管类比为弹性腔的外周阻力，心脏作为动力源，首先把血液压入弹性腔，之后血液再通过弹性腔的弹性压缩以及外周阻力的作用下，流入到人体各个组织。为了能够更真实地模拟心血管循环系统，学者们先后提出了三元件模型、四元件模型、五元件模型、改进的五元件模型等各种多元件的零维血流动力学模型。通过建立血流参数和电路参数的等效关系，将复杂的心血管系统简化为零维血流动力学模型(由各种电学元件构成的电路回路，即零维血流动力学模型)，进而通过零维血流动力学模型提供所代表区域量化的压力和流量信息实现对血液循环系统中血液流动特性定量的分析。因此可以通过建立冠脉零维血流动力学模型来模拟静息和充血状态下冠脉的血流动力学，提供冠脉各分支的压力流量波形，实现ffr的快速计算。
技术实现要素：
：
本发明提出一种基于零维血流动力学模型的ffr快速计算系统，该系统的优势在于基于零维血流动力学模型模拟冠脉的血流动力学，大大的简化了模型，减少了模型计算量，提高了仿真速度，对于快速计算ffr具有重要意义。
技术方案如下：
一种基于零维血流动力学模型的ffr快速计算系统模型，其包括如下模型：
(1)构建个性化零维血流动力学模型；
(2)构建狭窄阻力理论计算模型；
(3)确定狭窄阻力模型；
(4)计算ffr的模型。
其中所述(1)包括如下：
1.1基于冠脉cta图像提取出每一支冠脉血管长度l、血管直径d、横截面积a等冠脉的个性化解剖参数以及冠脉狭窄处直径ds、狭窄远端直径dd和狭窄长度ls等狭窄病变的个性化解剖参数；
1.2通过建立电学参数与血流动力学参数的等效关系，将复杂的心血管循环系统简化为电路结构(即零维血流动力学模型)；心血管系统的零维血流动力学模型主要由心脏模块、大动脉模块和冠脉模块组成。每个模块中电阻r模拟血流阻力，电容c模拟血管顺应性，电感l模拟血流惯性；心脏模块中，单向导电的二极管模拟心脏瓣膜，可变电容模拟心室的收缩舒张。冠脉模块中，每个冠脉分支由冠脉血管部分和冠脉远端负载(即冠脉后负荷)部分组成，在冠脉模块远端增加左心室压力模拟心肌收缩对冠脉血流的影响(实现舒张期供血)；心脏模块中用压力-容积关系描述一个心动周期内可变电容的变化：
e(t)为时变弹性(mmhg/ml)，是可变电容的倒数，p(t)为时变心室压力(mmhg)，v(t)为时变心室容积(ml)，v0是参考心室容积(ml)。数学上描述e(t)的公式为：
e(t)＝(emax-emin)·en(tn) emin式(1-2)
其中，emax是收缩末期心室压力-容积比，emin是舒张末期心室压力-容积比，en(tn)是归一化的时变弹性，表达式为：
其中，
1.3根据冠脉cta图像，确定冠脉分支结构；基于步骤1.2，结合零维血流动力学模型中冠脉标准16段模型，针对冠脉模块，根据个性化人自身情况删去多余的分支，添加个性化人自身情况特有的分支，实现对零维血流动力学模型结构(电路结构)的个性化；
1.4假设冠脉总流量是心输出量的4％，左右冠状动脉的流量分别占冠脉总流量的60％和40％；将平均动脉压(由公式1-4计算)作为主动脉入口压力；由公式(1-5)确定每支冠脉的阻力；然后根据q∝d3(d是血管直径)对冠脉各分支进行流量分配；由公式(1-6)计算冠脉各分支节点的压力；最后由流量和压力共同确定静息状态下(假设冠脉无狭窄)各分支的冠脉后负荷；
平均动脉压＝(收缩压 2*舒张压)/3式(1-4)
其中，μ是动态粘度(取0.0035pa·s)，l是血管长度，a是血管横截面积。
p＝q·r式(1-6)
其中，p是压力，q是流量，r是阻力；
上述可应用于不同的血管；
1.5使用遗传算法优化包括心脏在内的8个敏感参数(r0，c0，l0，ca，rdoap，rdoad，emax，emin)，以正常人的压力流量波形作为标准波形，根据标准波形结合个性化的收缩压、舒张压、心率和心输出量得到一个目标模型，当仿真波形与目标波形最大程度匹配(仿真波形与目标波形的均方根误差小于目标波形均值的5％)时，得到优化的模型参数；
1.6基于步骤1.4和1.5实现对零维血流动力学模型参数的个性化，结合步骤1.3中对零维血流动力学模型结构的个性化，实现对静息状态下(冠脉无狭窄)零维血流动力学模型的个性化。冠脉存在狭窄时，充血状态下的冠脉后负荷用静息状态下冠脉后负荷的0.24倍来模拟，并计算狭窄阻力，由此模拟充血状态下(冠脉存在狭窄)的零维血流动力学模型；
1.7通过fortran程序语言对零维血流动力学模型(电路结构)进行编译，在程序中定义所有电路元器件的参数值，并给定初值，结合kirchhoff定律，用常微分方程组表示电路结构中每个节点对应的电压或电流，用显示欧拉法对零维血流动力学模型进行求解；设置心动周期、步长和周期数，在编译软件中对零维血流动力学模型(fortran程序)进行仿真，输出压力和流量波形；
其中所述的(2)具体包括如下：
2.1霍云龙等人基于能量守恒提出计算狭窄分支压降的理论模型，理论模型中压降主要由扩散和膨胀导致的能量损失形成，计算步骤见公式(2-1)-(2-8)；在此基础上，增加公式(2-9)用来计算狭窄阻力；
2.2根据lstenosis和狭窄分支流量q(由零维血流动力学模型输出)，结合公式(2-1)和(2-2)确定α；
其中，lentrance是入口长度，lstenosis是狭窄长度，μ是动态粘度(取0.0035pa·s)，ρ是血液密度(取1050kg/cm3)，q是流量，α是定义的非粘核半径；
2.3当α≥0.05时,由公式(2-3)结合公式(2-4)、(2-5)和(2-6)确定δpstenosis；
其中，
当α<0.05时,先由公式(2-7)重新确定lentrance，再由公式(2-8)结合公式(2-6)确定δpstenosis。
其中，δpstenosis是通过狭窄分支的压降，csastenosis是狭窄处的横截面积，csadistal是狭窄远端的横截面积，
2.4由公式(2-9)确定rstenosis。
其中，rstenosis为狭窄阻力；
2.5由上述理论公式可知，狭窄阻力rstenosis主要与狭窄处横截面积csastenosis、狭窄远端横截面积csadistal和狭窄长度lstenosis以及狭窄分支流量q相关(即狭窄阻力理论计算模型是关于狭窄处直径、狭窄远端直径、狭窄长度和狭窄分支流量的函数)。用matlab对上述理论公式进行求解，每给定一个狭窄分支流量q可以计算出一个相应的狭窄阻力rstenosis；
其中所述的(3)具体包括如下：
3.1将提取的相关解剖参数赋给狭窄阻力理论计算模型(matlab程序)，将静息状态下的零维血流动力学模型(fortran程序)仿真计算达到收敛的最后一步(例如一共计算500步，取第500步计算结果)的狭窄处流量的4倍作为狭窄阻力理论计算模型(matlab程序)中狭窄分支流量的初值，由matlab程序运行可以得到一个相应的狭窄阻力；
3.2由充血状态下的零维血流动力学模型(fortran程序)和狭窄阻力理论计算模型(matlab程序)耦合求解个性化狭窄阻力：充血状态下的零维血流动力学模型每输出一个狭窄处流量(取仿真结果中最后一步的狭窄处流量)，狭窄阻力理论计算模型将重新计算狭窄阻力；狭窄阻力理论计算模型每计算出一个新的狭窄阻力，充血状态下的零维血流动力学模型将重新分配冠脉各分支的流量(即充血状态下的零维血流动力学模型为狭窄阻力理论计算模型提供流量的边界条件，狭窄阻力理论计算模型为充血状态下的零维血流动力学模型提供阻力的边界条件)。两个模型进行多次迭代直到流量与阻力匹配(即系统达到平衡，流量与阻力不再改变)，最终确定个性化狭窄阻力；
其中所述的(4)包括：
4.1根据步骤3.2确定的个性化狭窄阻力，结合步骤1.6，模拟充血状态下的零维血流动力学模型。对充血状态下的零维血流动力学模型(fortran程序)进行仿真计算，输出冠脉各分支压力流量波形；
4.2根据充血状态下的零维血流动力学模型输出的压力波形计算狭窄远端的压力平均值和主动脉根部的压力平均值，由公式(4-1)计算ffr；
其中，pd为狭窄远端的压力平均值，pa为主动脉根部的压力平均值。
本发明系统模型可以快速地计算ffr，根据本发明系统模型建立的快速计算ffr的零维血流动力学模型具有很高的准确性。
附图说明：
图1：基于零维血流动力学模型的ffr快速计算系统框图；
图2：心血管循环系统的零维血流动力学模型示意图；
图3：静息和充血状态下的压力波形图(临床实测ffr为0.91).
具体实施方式：
下面结合实施例对本发明做进一步说明，但本发明并不限于以下实施例。
实施例1
临床基于cta图像提取个性化解剖参数(每一支冠脉血管长度l、血管直径d、横截面积a及冠脉狭窄处直径ds、狭窄远端直径dd和狭窄长度ls等)，并采集个性化生理参数(收缩压、舒张压、心率和心输出量)。计算冠脉后负荷，优化模型参数，构建静息状态下(假设冠脉无狭窄)的个性化零维血流动力学模型。冠脉存在狭窄时，计算个性化狭窄阻力，用静息状态下冠脉后负荷的0.24倍模拟充血状态下(冠脉存在狭窄)的冠脉后负荷，构建充血状态下(冠脉存在狭窄)的个性化零维血流动力学模型，实现对血液循环系统的模拟。构建的心血管循环系统的零维血流动力学模型示意图(电路结构)如图2所示。
通过充血状态下的零维血流动力学模型(fortran程序)和狭窄阻力理论计算模型(matlab程序)耦合求解个性化狭窄阻力。将静息状态下的个性化零维血流动力学模型(fortran程序)仿真计算最后一步的狭窄处流量的4倍作为狭窄阻力理论计算模型狭窄处流量的初值，然后狭窄阻力理论计算模型(matlab程序)计算得到一个相应的狭窄阻力。之后，由充血状态下的零维血流动力学模型输出的狭窄处流量作为狭窄阻力理论计算模型的边界条件，可计算得到相应的狭窄阻力；然后将狭窄阻力理论计算模型中计算的狭窄阻力作为充血状态下的零维血流动力学模型的边界条件，重新分配冠脉各分支流量。两个模型进行多次迭代(这里选取输出结果中最后一步的狭窄处流量作为参考值)，当迭代结果与迭代前一次结果一致时，流量与阻力匹配，最终确定个性化狭窄阻力。
将个性化狭窄阻力赋给充血状态下的零维血流动力学模型(fortran程序)，仿真运行得到冠脉各分支的压力和流量波形。个性化零维血流动力学模型(fortran程序)仿真输出的静息和充血状态下的压力波形图如图3所示，提取计算收敛后3个周期的波形。取充血状态下狭窄远端的压力平均值与主动脉根部的压力平均值，两者的比值即为ffr。
选取10例进行试验，将临床测得的实际的ffr与本发明系统模型计算的ffr进行比较分析，准确率为86.4％，本发明系统具有一定的准确性。
1.一种基于零维血流动力学模型的ffr快速计算系统模型，其特征在于，包括如下模型：
(1)构建个性化零维血流动力学模型；
(2)构建狭窄阻力理论计算模型；
(3)确定狭窄阻力模型；
(4)计算ffr的模型。
2.按照权利要求1所述的一种基于零维血流动力学模型的ffr快速计算系统模型，其特征在于，
其中所述(1)包括如下：
1.1基于冠脉cta图像提取出每一支冠脉血管长度1、血管直径d、横截面积a等冠脉的个性化解剖参数以及冠脉狭窄处直径ds、狭窄远端直径dd和狭窄长度ls等狭窄病变的个性化解剖参数；
1.2通过建立电学参数与血流动力学参数的等效关系，将复杂的心血管循环系统简化为电路结构即零维血流动力学模型；心血管系统的零维血流动力学模型主要由心脏模块、大动脉模块和冠脉模块组成；每个模块中电阻r模拟血流阻力，电容c模拟血管顺应性，电感l模拟血流惯性；心脏模块中，单向导电的二极管模拟心脏瓣膜，可变电容模拟心室的收缩舒张；冠脉模块中，每个冠脉分支由冠脉血管部分和冠脉远端负载(即冠脉后负荷)部分组成，在冠脉模块远端增加左心室压力模拟心肌收缩对冠脉血流的影响(实现舒张期供血)；心脏模块中用压力-容积关系描述一个心动周期内可变电容的变化：
e(t)为时变弹性(mmhg/ml)，是可变电容的倒数，p(t)为时变心室压力(mmhg)，v(t)为时变心室容积(ml)，v0是参考心室容积(ml)。数学上描述e(t)的公式为：
e(t)＝(emax-emin)·en(tn) emin式(1-2)
其中，emax是收缩末期心室压力-容积比，emin是舒张末期心室压力-容积比，en(tn)是归一化的时变弹性，表达式为：
其中，
1.3根据冠脉cta图像，确定冠脉分支结构；基于步骤1.2，结合零维血流动力学模型中冠脉标准16段模型，针对冠脉模块，根据个性化人自身情况删去多余的分支，添加个性化人自身情况特有的分支，实现对零维血流动力学模型结构(电路结构)的个性化；
1.4假设冠脉总流量是心输出量的4％，左右冠状动脉的流量分别占冠脉总流量的60％和40％；将平均动脉压(由公式1-4计算)作为主动脉入口压力；由公式(1-5)确定每支冠脉的阻力；然后根据q∝d3(d是血管直径)对冠脉各分支进行流量分配；由公式(1-6)计算冠脉各分支节点的压力；最后由流量和压力共同确定静息状态下(假设冠脉无狭窄)各分支的冠脉后负荷；
平均动脉压＝(收缩压 2*舒张压)/3式(1-4)
其中，μ是动态粘度(取0.0035pa·s)，l是血管长度，a是血管横截面积；
p＝q·r式(1-6)
其中，p是压力，q是流量，r是阻力；
1.5使用遗传算法优化包括心脏在内的8个敏感参数(r0，c0，l0，ca，rdoap，rdoad，emax，emin)，以正常人的压力流量波形作为标准波形，根据标准波形结合个性化的收缩压、舒张压、心率和心输出量得到一个目标模型，当仿真波形与目标波形最大程度匹配(仿真波形与目标波形的均方根误差小于目标波形均值的5％)时，得到优化的模型参数；
1.6基于步骤1.4和1.5实现对零维血流动力学模型参数的个性化，结合步骤1.3中对零维血流动力学模型结构的个性化，实现对静息状态下(冠脉无狭窄)零维血流动力学模型的个性化。冠脉存在狭窄时，充血状态下的冠脉后负荷用静息状态下冠脉后负荷的0.24倍来模拟，并计算狭窄阻力，由此模拟充血状态下(冠脉存在狭窄)的零维血流动力学模型；
1.7通过fortran程序语言对零维血流动力学模型(电路结构)进行编译，在程序中定义所有电路元器件的参数值，并给定初值，结合kirchhoff定律，用常微分方程组表示电路结构中每个节点对应的电压或电流，用显示欧拉法对零维血流动力学模型进行求解；设置心动周期、步长和周期数，在编译软件中对零维血流动力学模型(fortran程序)进行仿真，输出压力和流量波形。
3.按照权利要求1所述的一种基于零维血流动力学模型的ffr快速计算系统模型，其特征在于，
其中所述的(2)具体包括如下：
2.1霍云龙等人基于能量守恒提出计算狭窄分支压降的理论模型，理论模型中压降主要由扩散和膨胀导致的能量损失形成，计算步骤见公式(2-1)-(2-8)；在此基础上，增加公式(2-9)用来计算狭窄阻力；
2.2根据lstenosis和狭窄分支流量q(由零维血流动力学模型输出)，结合公式(2-1)和(2-2)确定α；
其中，lentrance是入口长度，lstenosis是狭窄长度，μ是动态粘度(取0.0035pa·s)，ρ是血液密度(取1050kg/cm3)，q是流量，α是定义的非粘核半径；
2.3当α≥0.05时，由公式(2-3)结合公式(2-4)、(2-5)和(2-6)确定δpstenosis；
其中，
当α＜0.05时，先由公式(2-7)重新确定lentrance，再由公式(2-8)结合公式(2-6)确定δpstenosis；
其中，δpstenosis是通过狭窄分支的压降，csastenosis是狭窄处的横截面积，csadistal是狭窄远端的横截面积，
2.4由公式(2-9)确定rstenosis；
其中，rstenosis为狭窄阻力；
2.5由上述理论公式可知，狭窄阻力rstenosis主要与狭窄处横截面积csastenosis、狭窄远端横截面积csadistal和狭窄长度lstenosis以及狭窄分支流量q相关即狭窄阻力理论计算模型是关于狭窄处直径、狭窄远端直径、狭窄长度和狭窄分支流量的函数；用matlab对上述理论公式进行求解，每给定一个狭窄分支流量q可以计算出一个相应的狭窄阻力rstenosis。
4.按照权利要求1所述的一种基于零维血流动力学模型的ffr快速计算系统模型，其特征在于，
其中所述的(3)具体包括如下：
3.1将提取的相关解剖参数赋给狭窄阻力理论计算模型(matlab程序)，将静息状态下的零维血流动力学模型(fortran程序)仿真计算达到收敛的最后一步的狭窄处流量的4倍作为狭窄阻力理论计算模型(matlab程序)中狭窄分支流量的初值，由matlab程序运行可以得到一个相应的狭窄阻力；
3.2由充血状态下的零维血流动力学模型(fortran程序)和狭窄阻力理论计算模型(matlab程序)耦合求解个性化狭窄阻力：充血状态下的零维血流动力学模型每输出一个狭窄处流量即取仿真结果中最后一步的狭窄处流量，狭窄阻力理论计算模型将重新计算狭窄阻力；狭窄阻力理论计算模型每计算出一个新的狭窄阻力，充血状态下的零维血流动力学模型将重新分配冠脉各分支的流量，即充血状态下的零维血流动力学模型为狭窄阻力理论计算模型提供流量的边界条件，狭窄阻力理论计算模型为充血状态下的零维血流动力学模型提供阻力的边界条件；两个模型进行多次迭代直到流量与阻力匹配即系统达到平衡，流量与阻力不再改变，最终确定个性化狭窄阻力。
5.按照权利要求1所述的一种基于零维血流动力学模型的ffr快速计算系统模型，其特征在于，
其中所述的(4)包括：
4.1根据步骤3.2确定的个性化狭窄阻力，结合步骤1.6，模拟充血状态下的零维血流动力学模型。对充血状态下的零维血流动力学模型(fortran程序)进行仿真计算，输出冠脉各分支压力流量波形；
4.2根据充血状态下的零维血流动力学模型输出的压力波形计算狭窄远端的压力平均值和主动脉根部的压力平均值，由公式(4-1)计算ffr；
其中，pd为狭窄远端的压力平均值，pa为主动脉根部的压力平均值。
技术总结