[发明专利]一种基于心肌灌注的无创冠状动脉微循环检测方法

权利要求书

1.一种基于心肌灌注的无创冠状动脉微循环检测方法，其特征在于，所述方法包括以下步骤：

S01：冠脉造影：使用血管造影仪，以改良Seldinger法穿刺右桡动脉或右股动脉，按标准方法进行选择性冠状动脉造影；

S02：测定FFR、CFR、IMR：造影结束后，采用压力导丝结合温度稀释法的方式测量冠脉，当氯化钠液流入冠状动脉口时，压力导丝轴记录到第一条触发的温度曲线；当氯化钠液流至导丝头端感受器时，记录到第二条触发的温度曲线；两条曲线触发的时间差即为平均传导时间Tmn；

S03：重复上述操作2-4次，得到平均传导时间bTmn，记录静息及充血状态时的Pa和冠状动脉狭窄远端的压力Pd、FFR及CFR值，计算IMR值：

S04：进行冠状动脉CTA扫描及CT-MPI心肌动态灌注扫描，心肌动态灌注扫描在冠脉扫描后5min进行，得到冠脉血管三维扫描数据；

S05：图像后处理，在CTA扫描后：将获取的心脏收缩期或舒张期血管数据信息调入软件，进行VR彩图、血管树和MIP血管分析，对右冠状动脉RCA、左冠状动脉前降支LAD及左冠状动脉回旋支CX行管腔狭窄状况分析评估；

S06：计算左室壁厚度、曲率和体积，采用心内与心外膜表面垂直线法测定心肌室壁厚度T，多次测量取平均值，通过拟合心内膜表面的圆弧计算出该点的心室壁曲率R，并计算左室体积，室壁厚度T与曲率半径R的比值为左室壁压力指数σ；

S07：快速重建主动脉及相关分支血管，建立基于深度神经网络的冠脉三维重建，生成三维几何模型，选取流体方程，根据流体的雷诺数判断流体的流态，血流模拟使用非线性弹性理论；

S08：计算心肌血流量Q和血容量，，采用非结构网格加密与放粗技术，同时还采用多水平型区域分解技术来保证子区域的平衡性采用拉格朗日欧拉方法Arbitrary Lagrangian-Eulerian、间断Discontinuous Galerkin有限元方法和基于Lagrange乘子的新型方法，处理流体、固体以及流固耦合方程，采用NKS算法对大规模的稀疏线性或非线性系统进行计算心肌血流量Q；

S09:计算无创的CT-IMR，计算公式如下：

式中：Q为上述步骤计算得出的的心肌血流量，FFR_myo＝P_d/P_a，P_d和P_a可以通过计算流体方法无创获得，Tmn表示热稀释曲线记录液体两温敏感受器之间的时间，它与血流速度呈反比，即V＝L/Tmn，V表示冠脉血流平均速度，L表示两温敏感受器之间的距离，r为血管的半径。

2.根据权利要求1所述的一种基于心肌灌注的无创冠状动脉微循环检测方法，其特征在于，所述IMR计算方式为：当冠状动脉轻中度狭窄、血流储备分数FFR＞0.80时，使用简化公式：IMR＝Pd*Tmn；当冠状动脉重度狭窄、FFR≤0.80，IMR＝Pa*Tmn[(Pd-Pw)/(Pa-Pw),其中Pw表示冠状动脉楔压。

3.根据权利要求1所述的一种基于心肌灌注的无创冠状动脉微循环检测方法，其特征在于，所述血管造影仪采用GE公司IGS520或Phillips血管造影仪。

4.根据权利要求1所述的一种基于心肌灌注的无创冠状动脉微循环检测方法，其特征在于，所述冠脉CTA扫描：若病人心率次≤70beat/min，全剂量窗设在65％～75％的R-R间期；若病人心率次70best/min，则全剂量窗设在35％～45％的R-R间期。

5.根据权利要求1所述的一种基于心肌灌注的无创冠状动脉微循环检测方法，其特征在于，所述CT-MPI扫描：启动扫描程序，心电门控触发点设置在收缩末期，在R峰后250ms。动态穿梭扫描过程历时30s，灌注图像层厚3mm，间隔2mm，卷积取B25。