[发明专利]基于血流储备分数和血流阻力模型的微循环阻力快速计算方法

说明书

本发明提出一种结合血流储备分数(FFR)和血流阻力模型的快速计算冠状动脉微循环阻力的方法，包括：基于冠状动脉CT造影图像计算感兴趣血管段血流阻力模型参数；将冠状动脉微循环阻力模拟为纯阻力单元并串联于计算得到的血流阻力模型之后，构成集总参数模型；以FFR测量的近端压力为已知条件并预设冠状动脉微循环阻力数值，求解方程得到远端压力值；对比仿真的远端压力值与FFR测量的远端压力值，迭代的调整冠状动脉微循环阻力数值，直至计算与仿真的差值小于预设阈值；提取调整后的模型微循环阻力数值作为冠状动脉微循环阻力计算值；以侵入式FFR测量结果为已知条件，基于血流阻力模型，快速仿真获得与FFR测量结果相匹配的血流量及微循环阻力。本发明结合血流阻力模型，在避免了同步测量血流量的同时，提高了准确性；同时相较于其他迭代求解方法，本发明使用血流阻力模型将迭代求解流体控制方程简化为求解一元二次方程，速度得到巨大提升。

技术领域

本发明涉及微循环计算领域，尤其涉及应用在结合血流储备分数(FractionalFlow Reserve, FFR)和血流阻力模型计算微循环阻力的方法。

背景技术

本发明是基于本人公开号为CN107411767A，名为“一种基于冠状动脉CT血管造影评估狭窄病灶血流阻力的非侵入式方法”的专利发明的发明。

冠状动脉微循环阻力(microcirculatory resistance,MR)是反映冠状动脉微循环功能的有效指标，其定义为冠脉远端动脉压力(P_d)除以最大充血状态下的冠脉血流量(Q_max)。现有技术主要分为侵入式和非侵入式两种：

第一类方法通过集成了压力传感器和血流量传感器的导丝置入冠状动脉远端，在药物诱导最大充血状态的条件下同时测量压力和血流量，进而计算MR。依据采用的血流量传感器的不同，又可分为热稀释法和超声多普勒测速法(Williams,R.P.,et al.(2018).Doppler Versus Thermodilution-Derived Coronary Microvascular Resistance toPredict Coronary Microvascular Dysfunction in Patients With Acute MyocardialInfarction or Stable Angina Pectoris.Am J Cardiol 121(1):1-8.)。

第二类方法主要通过序列造影图像估算Q_max，基于血流动力学模型仿真计算最大充血状态下远、近端压力差(ΔP)，基于测量的静息态冠脉入口压力和深度学习模型估算最大充血状态下冠脉入口压力(P_a)(“快速计算微循环阻力的方法与系统”，中国201711258493.9[P]) 或直接测量最大充血状态下P_a(“基于造影图像和流体力学模型的微循环阻力指数计算方法”，中国201810413391.8[P])。

上述技术尽管从不同角度、不同计算方法中给出了确定MR的方法，但其都至少具有以下这个技术缺陷，即无法精确的测量(或计算)最大充血状态下的Q_max：第一类方法由于测量原理以及导管置入对血流量的影响等原因，而造成无法精确测量Q_max；第二类方法采用估算Q_max的方式，精度难以保证。

发明内容

有鉴于此，本发明提供一种结合FFR和血流阻力模型的快速计算微循环阻力的方法。该方法将侵入式FFR测量远、近端压力值作为已知条件，基于血流阻力模型快速求解冠状动脉微循环阻力(MR)。所采取的技术方案如下：

1.基于冠状动脉CT血管造影图像计算感兴趣血管段血流阻力模型参数；

2.将MR模拟为纯阻力单元并串联于计算得到的血流阻力(压降-血流曲线)模型之后，构成集总参数模型；