专利名称：微型轴流式血泵的优化非恒速控制方法
技术领域：
微型轴流式血泵的优化非恒速控制方法属于人工器官控制技术领域。
目前，微型轴流式血泵在临床应用时主要是工作在恒定转速。据临床实验和动物实验观察，这种工作方式目前尚存在以下几个问题1)在临床应用中，轴流式血泵的转速越高，产生的心输出量越大。但是，泵的转速过高时会导致左心室“塌陷现象”，从而引起血泵的引流导管前端与心室壁接触，妨碍泵的入流状况，甚至可能伴随心律不齐的产生，对左心室产生很大的破坏作用。同时，高速旋转的叶片会产生较大的剪切力，从而对红细胞产生较大的破坏作用。这就提出了这样一个问题对于某一病人特定的生理状况来说，如何确定轴流式血泵的最佳转速？2)尽管Yokihiko在动物实验中曾经证明长期非脉动的双心辅助不会产生有害的生理效果，但是近期的研究结果显示长期的非脉动左心血流可能会对主动脉形态和功能产生一系列的影响，包括血管壁的厚度、血管壁平滑肌的结构和血管收缩功能的改变等。因此，考虑到微型轴流式血泵用于长期的循环辅助的可能性，可以尝试让血泵工作在非恒速工作状态，从而得到与自然生理状况相似的脉动血流。
本发明的特征在于它是一种建立在由目标函数，控制参数和限制条件组成的调节血泵转速的优化模型之上的微型轴流式血泵的优化非恒速控制方法；它依次含有以下几个步骤(1)设定微型轴流血泵Hemopump的工作方式，以一个心拍为单位计，其控制参数如下Hemopump共有7种速度可选，第一级为17000转/分，邻级间隔1500转/分，VB为高速，VA为低速；tAB为从VA到VB的转换时刻；tBA为从VB到VA的转换时刻；Tpb为Hemopump以VA持续运转的时间，Tpe为Hemopump以VB持续运转的时间；设转换过程时间可不计。每一控制方案对应四个控制参数VA，VB，Tpb，Tpe；(2)根据设定的Hemopump控制方案，通过实体实验，利用医疗仪器测出或控制参数计算一个心拍内五个表征心脏辅助状况的独立生理参量，为利用医疗仪器测出每搏心输出量(一个心拍内的心脏搏血量)、心肌供氧量与耗氧量的比值、主动脉最大收缩压和最小舒张压，利用控制参数计算Hemopump在一个心拍内的平均转速；各独立参量的测量和计算方法如下每搏心输出量(SV)三维B超可以测出一个心拍内心室容积变化，每搏心输出量为心室舒张末期容积与收缩末期容积之差；心肌供氧量与耗氧量(Vop为两者之比)可用新型心血管功能测试仪测出；主动脉最大收缩压(Paps)和最小舒张压(Pamd)有创血压计跟踪测出；Hemopump在一个心拍内的平均转速(Speed)speed‾=VA*Tpb+VB*TpeT;Tpb=T-tBA+tAB;Tpe=tBA-tAB.]]>T为一个心拍的持续时间，近似为Tpb+Tpe；(3)把上述SV，Vop，Paps，Pamd四个生理参数和Speed投射到0，1的区域，用它们的隶属函数去反映参数在不同范围内取值时的生理状况；(4)计算各参数的隶属函数并按权重相加得到一组Tpb、Tpe下的目标函数值TfTf＝w1*u(SV)+w2*u(Vop)+w3*u(Paps)+w4*u(Pamd)+w5*u(Speed)；其中u(SV)、u(Vop)、u(Paps)、u(Pamd)和u(Speed)分别代表每博心输出量SV、心肌供氧量与耗氧量的比值Vop、主动脉最高收缩压Paps、主动脉最低舒张压Pamd和血泵平均转速Speed的隶属函数；W1-W5为相应各参数的权重；
(5)改变Tpb和Tpe，重复以上步骤，得到一组VA，VB下目标函数最大值及其相应的Tpb、Tpe；(6)改变VA，VB，Tpb，Tpe，重复以上步骤，得到各种速度组合下目标函数的最大值及相应的Tpb、Tpe；(7)从各种速度搭配的最大值中找一个最大值，其VA，VB，Tpb，Tpe即为最优方案。
所述各参数SV，Vop，Paps，Pamd和Speed的隶属函数如下
所述目标函数的权重因子是通过构造各参数的成对比较矩阵及它的最大特征值得到的。
试验证明这种优化非恒速控制方法在舒张期后期才提高泵的转速，而在左心室体积较小时则减小泵的转速，确保了回心血量能满足血泵的流量供给，从而避免心塌陷现象的发生，同时也减小对红细胞的破坏作用。
图2.优化控制方案设计程序流程框图。
图3.微型轴流式血泵速度转换时间的不同组合。
图4.左心室30％质量心肌缺血，右心正常时，各种控制方案下目标函数的比较。
2)优化控制利用有关模糊理论，构造一个目标函数，使得这个函数可以正确反应心脏辅助效果，函数值的大小反应效果好坏。选择效果作用的各种参量要求作用结果是单一的，经分析，最后选中每博心输出量SV、心肌供氧量与耗氧量的比值Vop、主动脉最高收缩压Paps和主动脉最低舒张压Pamd这四个生理参数以及血泵的平均转速Speed，由于各个参数代表不同的生理量，以及不同的单位和取值范围，构造目标函数时必须对参数实行归一化，采用模糊数学中隶属函数的概念，将参数投射到0，1的区域，各参数的隶属函数反映了参数在不同范围内取值所反映的生理状况。
3)将各项隶属函数值代入式(6)定义的目标函数，计算得到某一控制方案所对应的目标函数的取值；(6Tf＝w1*u(SV)+w2*u(Vop)+w3*u(Paps)+w4*u(Pamd)+w5*u(Speed))其中u(SV)、u(Vop)、u(Paps)、u(Pamd)和u(Speed)分别代表每博心输出量SV、心肌供氧量与耗氧量的比值Vop、主动脉最高收缩压Paps、主动脉最低舒张压Pamd和血泵平均转速Speed的隶属函数；W1-W5为权重因子，表示各个参量在整个目标函数中所起作用的大小，各隶属函数以及权重因子的推导和确定见后式(1)-式(5)及注释。
4)比较各种控制方案下得到的目标函数的取值，目标函数的最大值即对应着最优的控制方案。
本优化非恒速控制方法能够很好的区分处于不同健康水平下的使用者对辅助剂量的不同需求，实现最佳辅助效果。
本发明实施例中针对使用者(犬)所处于的不同的生理状况，设计出各种情况分别所对应的最优控制方案。具体操作过程为1.确定各参量每博心输出量SV、心肌供氧量与耗氧量的比值Vop、主动脉最高收缩压Paps、主动脉最低舒张压Pamd和血泵平均转速Speed的隶属函数，分别如下式(1)-式(5)所示； 其中(1)(2)(5)由自己确定，(3)(4)(参考已发表高血压患者健康状况评估方法的研究/杨德超，白净，王苏中，崔树起——《北京生物医学工程》2002.9第21卷第3期(187))的工作，隶属函数的确定简单介绍如下，以(1)为例，(2)和(5)也同样得出。
隶属度的含义是代表自然心脏生理状况的正常水平，隶属度的确定依据隶属度与代表心脏生理状况的变量相关，隶属度的变化必须符合实际的生理系统的变化情况，在理想状态下为1，最不理想状态为0。
心输出量的隶属函数为例，心室缺血心肌与正常心肌的最大弹性系数比值越大(理想为1)，心脏越正常，左右心室对比，左心室起作用更大，加一指数为心肌缺血率(缺血心肌质量占总心肌质量比率，可以用体积比率近似)，心肌缺血率越高，心脏越不正常，这样设立如下函数U0(sv)=(EesEes0)*(Ees2Ees1)Rm]]>Rm为心肌缺血率，可通过冠状动脉造影技术观测冠脉的堵塞情况，Rm＝堵塞冠脉体积/冠脉总体积；Ees，Ees0为衰竭和正常的右心室心肌的最大弹性系数，Ees2，Ees1为左心室缺血区域和正常区域心肌的最大弹性系数，均可通过三维B超测量近似可取对应区域收缩末期的压强/收缩末期区域体积即为最大弹性系数。
利用医疗仪器测量计算不同条件下的SV，利用医疗仪器测量SV对应的Rm、Ees0、Ees、Ees1、Ees2以计算U0(sv)，然后根据最小误差平方和的准则对心输出量的隶属函数以心输出量为唯一自变量进行拟合，得出最后结论(1)。
心肌供氧量与耗氧量比值Vop的隶属函数，右心缺血时，心肌供氧量和耗氧量的比值不能正确的反映心脏的生理状况(心肌供氧量减少可能比心肌耗氧量减少慢，比值可能增加)，心肌供氧量与耗氧量比值越大，心脏越正常，所以，只考虑右心正常的数据，左心室缺血心肌的最大弹性系数与正常心肌比值越大，缺血区域越少，心脏越正常，设立隶属函数为U0(Vop)=(Ees2Ees1)Rm,]]>利用医疗仪器测量心肌供氧量与耗氧量计算不同条件下的Vop，利用医疗仪器测量Ees2、Ees1和Rm计算U0(Vop)，然后根据最小误差平方和的准则对隶属函数以心肌供氧量与耗氧量比值Vop为唯一自变量进行拟合，得出最后结论(2)。
Hemopump平均转速的隶属函数确定，设旋转叶片对血细胞的破坏作用与泵的转速大小成线性关系，转速0时无破坏，隶属度为1，泵的转速为最大转速的两倍(52000转/分，实验观察得出)，叶片对血细胞的破坏已被认为不可接受，隶属度为0，这样，得到隶属函数为 即结论(5)五个隶属函数均可运用于实体实验。
2.确定在右心正常和右心衰竭两种不同的生理状况下，目标函数权重因子的取值分别如式(6)-式(7)[w1w2w3w4w5]＝
(6)[w1w2w3w4w5]＝
(7)确定方法如下右心健康时，可构造各参量的成对比较矩阵A如下u(SV)u(Vop)u(Paps)u(Pamd)u(Speed)13 5 5 3 u(SV)1/3 1 3 3 1 u(Vop)1/5 1/3 1 1 1/3 u(Paps)1/5 1/3 1 1 1/3 u(Pamd)1/3 1 3 3 1 u(Speed)其中aij代表元素i对元素j的相对重要性。各种重要性对应的数值大小如下
相对重要程度 定义1重要程度相等3略微更重要5非常重要7明显重要9绝对的重要2，4，6，8 介于相邻的两个判断之间上述数值的倒数 j对i相对重要度为i对j相对重要度的倒数可求得A的最大特征值为5.05564，一致性因数(CI)和一致性比值(CR)如下CI＝(5.05564-n)/(n-1)；RI＝1.12(n＝5)；CR＝CI/RI＝0.01242＜10％CR表示关系矩阵对于一致性的偏离程度，若是小于10％，认为矩阵A对一致性的偏离程度是可以接受的，计算其主特征向量(即矩阵最大的特征值所对应的特征向量)并归一，即得(6)。
右心衰竭时，心肌供氧量对耗氧量的比值不能正确反应心脏的生理状况，所以目标函数中去掉Vop，并适当提高心输出量的权重，得到其成对比较矩阵B为u(SV)u(Paps)u(Pamd)u(Speed)15 2 9 u(SV)1/5 1 1/33 u(Paps)1/2 3 1 5 u(Pamd)1/9 1/31/51 u(Speed)按照右心正常时的算法同样计算可得(7)。
3.改变4个控制参数的取值，组合得到17199种不同的控制方案。分别计算在各种控制方案下对应的目标函数的取值。其中，改变各控制参数得到17199种不同控制方式的过程如下(1)速度A在Hemopump的7档速度速度1-速度7之间选择；
(2)当速度A确定时，速度B在速度1-速度7之间选择；(3)速度A和速度B都确定后，改变一个心动周期中速度A到速度B的转换时间tAB。tAB的取值可以从心脏收缩期开始时刻t0变化至心脏舒张期结束时刻t0+T(其中T为一个心动周期的持续时间，当心率设为120搏/分时，T对应0.5秒)，变化间隔为0.02秒；(4)当速度A、速度B和tAB都确定时，改变一个心动周期中速度B到速度A的变化时间tBA。tBA的取值可以从tAB变化到t0+T，变化间隔同样为0.02秒。这样在速度A和速度B的任一种组合下，tAB和tBA共有351种组合，如说明书附2中所示。
图2是本发明提出的优化非恒速控制方法的设计程序流程框图。
通过对图4的分析，可以得到该生理状态下目标函数最大时对应的控制方案。tAB是一个心拍中，轴流泵转速由VA转到VB的时刻；tBA是一个心拍中，轴流泵转速由VB转到VA的时刻，见
图1。由图4知tAB＝220毫秒，tBA＝380毫秒时，目标函数值最大，为0.9098。
当左心轻度缺血时(30％心肌缺血)，生理系统需要的辅助剂量也较小(速度1和速度3的组合)。同时，由于此时左心自身的收缩功能还会对Hemopump的泵流量产生一定的影响，因此，在速度1和速度3的组合下，最佳的转换时刻是在收缩末期由低转速(速度1)转换为高转速(速度3)，在舒张后期由速度3转回为速度1，以减小左心室收缩时的后负荷。这一结果与主动脉内气囊反博的最佳充、排气时间相近。
类似的，可以得到不同的49种速度搭配下的各种控制方案的目标函数的比较，从而选出各自的最优方案，最后，将所得49个最优方案进行比较，得到最终解，此心脏生理状态下的最优控制方案；图3是一部分速度匹配方法。
权利要求
1.微型轴流式血泵的优化非恒速控制方法，其特征在于它是建立在由目标函数，控制参数和限制条件组成的调节血泵转速的优化模型之上的微型轴流式血泵的优化非恒速控制方法；它依次含有以下几个步骤(1)设定微型轴流血泵Hemopump的工作方式，以一个心拍为单位计，其控制参数如下Hemopump共有7种速度可选，第一级为17000转/分，邻级间隔1500转/分，VB为高速，VA为低速；tAB为从VA到VB的转换时刻；tBA为从VB到VA的转换时刻；Tpb为Hemopump以VA持续运转的时间，Tpe为Hemopump以VB持续运转的时间；设转换过程时间可不计，每一控制方案对应四个控制参数VA，VB，Tpb，Tpe；(2)根据设定的Hemopump控制方案，通过实体实验，利用医疗仪器测出或控制参数计算一个心拍内五个表征心脏辅助状况的独立生理参量，为利用医疗仪器测出每搏心输出量(一个心拍内的心脏搏血量)、心肌供氧量与耗氧量的比值、主动脉最大收缩压和最小舒张压，利用控制参数计算Hemopump在一个心拍内的平均转速；各独立参量的测量和计算方法如下每搏心输出量(SV)三维B超可以测出一个心拍内心室容积变化，每搏心输出量为心室舒张末期容积与收缩末期容积之差；心肌供氧量与耗氧量(Vop为两者之比)可用新型心血管功能测试仪测出；主动脉最大收缩压(Paps)和最小舒张压(Pamd)有创血压计跟踪测出；Hemopump在一个心拍内的平均转速(Speed)speed‾=VA*Tpb+VB*TpeT;Tpb=T-tBA+tAB;Tpe=tBA-tAB.]]>T为一个心拍的持续时间，近似为Tpb+Tpe；(3)把上述SV，Vop，Paps，Pamd四个生理参数和Speed投射到0，1的区域，用它们的隶属函数去反映参数在不同范围内取值时的生理状况；(4)计算各参数的隶属函数并按权重相加得到一组Tpb、Tpe下的目标函数值TfTf＝w1*u(SV)+w2*u(Vop)+w3*u(Paps)+w4*u(Pamd)+w5*u(Speed)；其中u(SV)、u(Vop)、u(Paps)、u(Pamd)和u(Speed)分别代表每博心输出量SV、心肌供氧量与耗氧量的比值Vop、主动脉最高收缩压Paps、主动脉最低舒张压Pamd和血泵平均转速Speed的隶属函数；W1-W5为相应各参数的权重；(5)改变Tpb和Tpe，重复以上步骤，得到一组VA，VB下目标函数最大值及其相应的Tpb、Tpe；(6)改变VA，VB，Tpb，Tpe，重复以上步骤，得到各种速度组合下目标函数的最大值及相应的Tpb、Tpe；(7)从各种速度搭配的最大值中找一个最大值，其VA，VB，Tpb，Tpe即为最优方案。
2.根据权利要求1所述的微型轴流式血泵的优化非恒速控制方法，其特征在于，所述各参数SV，Vop，Paps，Pamd和Speed的隶属函数加下
3.根据权利要求1所述的微型轴流式血泵的优化非恒速控制方法，其特征在于所述目标函数的权重因子是通过构造各参数的成对比较矩阵及它的最大特征值得到的。
全文摘要
微型轴流式血泵的优化非恒速控制方法属于人工器官控制技术领域，其特征为它是建立在由目标函数、控制参数、限制条件组成的调节血泵转速的优化模型之上的方法。此非恒速控制方法的控制参数为可选择7种速度的血泵高低转速以及高低速持续运转时间；优化模型的目标函数代表心脏生理健康状况，由五个独立心脏辅助生理参量的隶属函数加权相加构成；优化控制方法是通过实体实验，利用医疗仪器测出计算得到心输出量、心肌供氧量与耗氧量比值、主动脉最高收缩和最低舒张压，利用控制参数计算得到血泵平均转速这样共五个独立参量计算隶属函数经加权后组成目标函数，选取各种速度匹配方案的优化方案比较得到的。它适应性强，使用安全，血泵能耗也低。
文档编号G06F17/00GK1446592SQ03121948
公开日2003年10月8日 申请日期2003年4月18日 优先权日2003年4月18日
发明者白净, 张永红, 李熊, 夏东栋 申请人:清华大学