1. 引言
呼吸机是重症监护室中至关重要的生命支持医疗设备之一,但呼吸机的不当使用易导致管路积水等问题。不仅不利于治疗,而且还会对患者肺部造成伤害,使患者出现一系列并发症 [1] 。同时,呼吸机管路积水也会造成流量传感器和呼出阀等部件损坏的问题,导致呼吸机故障或失效,影响病患治疗效果。因此,对呼吸机管路积水故障进行监测具有重要意义 [2] 。
刘香君 [3] 等根据不同故障类别下端口时序性电信号及对应征兆现象,将预处理后的特征编码进行归一化处理并融合,筛选特征值并建立数据特征集,构建长短时记忆网络故障监测模型,将训练后多模态特征输入模型中,实现医疗设备故障智能诊断。陈艳 [4] 等分别收集造成呼吸机故障的多方面因素数据,采用收集到的数据训练搭建的8输入3输出粗糙神经网络,将待分析数据输入至训练后模型中,实现呼吸机故障监测。虽然上述方法均能够有效实现医疗设备和呼吸机故障的智能诊断和监测,但上述方法未对输入信号进行去噪处理,监测效果还需不断提升。
为了解决上述方法中存在的问题,本文提出了呼吸机管路积水故障自动监测算法研究。
2. 呼吸机管路信号去噪
监测呼吸机管路积水故障的前提是采集呼吸机相关信号数据,在采集信号过程中,外界环境和设备自身因素会导致采集到的呼吸机管路信号中包含一定噪声,因此为保障后续故障信号监测效果,需要去除信号中的噪声,提升呼吸机管路信号的真实性和可靠性。
所提方法结合实际去噪需求,采用小波分析中的小波变换去噪法去除呼吸机管路信号中的噪声 [5] [6] 。假设
表示母小波,
表示
的傅里叶变换 [7] [8] [9] ,
表示平方可积实数空间,t表示
时间,
表示频率,R表示空间半径,a和b表示平移和伸缩因子,当
满足
时,经平移和伸缩可得到小波序列
如下所示:
(1)
假设存在待处理信号
,将该信号与
做内积,可得到连续小波变换表达式
如下所示:
(2)
时频窗能够依据
和
进行自适应调整,通过小波变换实现信号高频和低频成分的自适应描述。
在小波变换去噪中,首先选取合适的小波函数将采集到的呼吸机管路信号分解,然后依据设定的阈值处理分解后信号参数,最后对参数执行逆小波变换,重构分解后信号,实现呼吸机管路信号去噪。小波分析去噪可以看作是信号的滤波问题,而小波降噪则是将低通滤波与特征提取相结合 [10] 。小波去噪和小波分解流程如图1所示。
![](//html.hanspub.org/file/136-2571684x23_hanspub.png?20240604170829550)
Figure 1. Flow chart of wavelet denoising and wavelet decomposition
图1. 小波去噪和小波分解流程图
小波去噪主要流程如图1(a)。用
代表信号,
代表含噪信号,
代表噪声系数标准偏差,
代表高斯白噪声 [11] ,则
,噪声通常表现为高频信号形式,而有价值信号往往为低频或平稳信号形式,因此通过图1(b)中小波分解处理信号能够达到信号去噪的效果,其中
、
、
表示低频信号,
、
、
表示高频信号,对
、
、
作相应的小波系数处理并重构小波系数,即可实现呼吸机管路积水信号去噪。
3. 多输出支持向量机故障自动监测
将呼吸机管路积水信号进行去噪后,对其进行故障监测。冷凝水处理不当、呼吸机整理和放置不合理、接水杯未及时倾倒、湿化罐水位超线、痰液堵塞、气道痉挛等因素是造成呼吸机管路积水的常见原因,为了有效区分不同原因造成的呼吸机管路积水问题并及时解决,所提方法采用多输出支持向量机构建呼吸机管路积水故障自动监测模型 [12] ,引入健康度指标量化呼吸机管路健康程度,通过监测健康度的劣化过程评估呼吸机管路系统性能下降情况。
3.1. 支持向量机回归算法
支持向量机建立在VC维理论和SRM准则上 [13] ,是一种智能机器学习算法,在提升泛化能力、避免陷入局部最优、解决维数灾难等方面均具有不俗效果。
用
表示训练集,其中
表示输入的去噪后信号指标向量,
,
表示输出指标,
,
,
,
表示法向量,
表示位移项,
表示预测值,由此可将支持向量机回归问题视为采用
搜索
中实值函数
问题,进而实现全部输入指标向量对应输出指标的计算。
在支持向量机中引入
不敏感损失函数,设置不敏感损失系数
取值,若在
处存在
,则认为
为无损预测值。
依据SRM准则,将松弛变量
、惩罚因子
、
添加到支持向量机中,转化支持向量机优化问题为下式优化问题:
(3)
通过对偶理论将公式(3)所得最优值设定为原始问题最优值下界 [14] 。用
表示拉格朗日乘子 [15] [16] ,
,依据拉格朗日函数转化公式(3)问题为对偶问题,如下所示:
(4)
对上式的凸二次规划问题求解可得到最优解
,
,结合
求解法向量和位移项的最优解
和
,如下所示:
(5)
支持向量机多用于解决非线性回归问题,引入核函数
变换输入指标向量
,生成非线性回归函数,如下所示:
(6)
选取合适的核函数和相关参数后构造支持向量机,即可实现输入数据的准确识别。
3.2. 多输出支持向量机回归算法
准确识别输入数据后,采用多输出支持向量机回归算法改进原有的支持向量机回归算法,提高监测精度。支持向量机多为多输入单输出的支持向量机,但呼吸机管路健康度是持续改变的过程值,多输出的支持向量机更利于呼吸机管路健康度的多步监测,因此所提方法改进支持向量机回归算法,将多输出支持向量机用于呼吸机管路健康度监测之中,从而实时掌控呼吸机管路的健康度下降情况。
与单输出支持向量机不同,多输出支持向量机采用矩阵替代单输出支持向量机参数,即将法向量
转变为矩阵
,位移项
转变为矩阵
,且输出指标
为向量,多输出支持向量机回归函数
,其中,
,
,
,在此基础上,利用
值不敏感的损失函数来求解凸规划问题,如下所示:
(7)
引入松弛变量
和惩罚因子
,
,构造优化目标如下:
(8)
同样依据拉格朗日函数转化上式问题为对偶问题,用
表示拉格朗日乘子,
,得到对偶问题如下所示:
(9)
求解上式可得到最优解
,
且
,
,
,求解
和
的最优解
和
,如下所示:
(10)
用
表示核函数,所提方法引入高斯径向基核作为核函数,用
表示宽度参数,得到多输出支持向量机非线性回归函数如下所示:
(11)
由以上过程可以看出,多输出支持向量机的监测精度主要由
、
和
三个参数决定。
3.3. 多输出支持向量机参数优化
为了进一步提高监测算法的精准度,所提方法引入遗传算法 [17] [18] 改进人工蜂群算法 [19] [20] ,进而提升种群适应性和多样性,使算法快速收敛于全局最优解,将改进后人工蜂群算法用于多输出支持向量机的
、
和
三个参数优化之中。
蜜源和参数最优解一一对应,用
表示蜜源
的第
维位置,
表示蜜源
的第
维位置,
表示蜜源总数,
,
表示蜜源维数,即目标函数参数总数,
,其中
,
表示在[−1, 1]中均匀分布的随机数,用于引领蜂搜索范围的调控,引领蜂基于旧位置
生成临近竞争新蜜源位置
的规则如下所示:
(12)
跟随蜂依据引领蜂在蜜源
处释放花蜜信息和花蜜量概率
选取蜜源,用
表示蜜源
花蜜量,即适应度值,
表示当前循环次数,
表示对应适应度值,则
计算方式如下所示:
(13)
用
表示在上次迭代中侦查蜂获取到新蜜源位置经历的迭代次数,依据
确定侦查蜂是否需要搜索蜜源。
是否成功以适应度值更新成功次数在总循环次数中的占比决定。
![](//html.hanspub.org/file/136-2571684x135_hanspub.png?20240604170829550)
Figure 2. Improved artificial bee colony algorithm flowchart
图2. 改进的人工蜂群算法流程图
改进的人工蜂群算法与人工蜂群算法不同,改进后侦查蜂首先在搜索新蜜源前依据遗传算法选择因子再次筛选跟随蜂选取的优质蜜源,选择最优蜜源直接遗传至下代,然后借鉴遗传算法中交叉因子将重组后父代某段结构互相替换,使后代蜜源同时继承上代双亲优点,提升搜索性能,最后在交叉操作趋近于最优解区域时,该方法利用遗传算法的变异运算,实现了局部优化,提高了算法的全局优化能力和收敛速度,并在一定程度上,保证了种群的多样性,避免了过早的收敛。
完成遗传算法相关操作后设定循环数
,若达到
次循环后
仍未得到优化,则舍弃
,由侦查蜂随机生成新蜜源替代原始蜜源,用
和
表示第
维变量上界和下界,新蜜源位置生成规则如下所示:
(14)
最终,依据适应度的变化情况确定新种群。改进的人工蜂群算法主要流程如图2所示。
图2中,
表示取整函数,
和
表示
的更新步长。
通过改进的人工蜂群算法得到多输出支持向量机最优参数后对支持向量机加以训练,建立多输出支持向量机拟合监测信号与健康度变化之间的非线性映射关系,将呼吸机管路数据输入训练后的多输出支持向量机中,对呼吸机管路数据多步监测,输出健康度监测结果,实现呼吸机管路积水故障自动监测。
4. 实验与结果
为了验证呼吸机管路积水故障自动监测算法研究的有效性,设定多输出支持向量机参数取值范围分别为:
、
和
。在MATLAB环境下对改进的人工蜂群算法优化多输出支持向量机参数模型进行仿真,通过对训练样本进行训练,得到多输出支持向量机参数优化结果。
![](//html.hanspub.org/file/136-2571684x151_hanspub.png?20240604170829550)
Figure 3. Waveform diagram of monitoring water accumulation faults in the ventilator pipeline using three methods
图3. 三种方法的呼吸机管路积水故障监测波形图
为了验证所提方法对呼吸机管路积水故障监测效果,分别采用文献 [3] 方法和文献 [4] 方法与所提方法进行对比。若呼吸机管路内存在积水,则在呼吸机管路积水故障监测中,易出现高频波动噪声,压力幅值变化越大,说明检测到的噪声频率越大,表明方法的监测效果越好。三种方法的呼吸机管路积水故障监测波形如图3所示。
根据图3可知,文献 [3] 方法和文献 [4] 方法的故障监测压力幅值变化较小,而所提方法的故障监测压力幅值变化较大,说明所提方法能够有效、全面地监测到呼吸机管路中的积水。由此可知,所提方法的呼吸机管路积水故障监测效果较好。
进一步验证三种方法在不同信噪比下的故障监测精度,每次实验训练批次为20次,迭代次数为50次,每组实验在完成10次后,求取故障监测准确率均值,故障监测准确率越高,表明方法的故障监测精度越高。实验训练系统及三种方法在不同信噪比下故障监测准确率如图4、表1所示。
![](//html.hanspub.org/file/136-2571684x152_hanspub.png?20240604170829550)
Figure 4. Screenshot of experimental training system
图4. 实验训练系统截图
由图4可以看出实验训练系统设定了声量控制器的参数。其中,基本参数包括:潮汐量(Tidal Volume) 500 ml,呼吸率(Resp.Rate)为15 b/min,呼气暂停时间(PEEP)为5 cm H2O,氧气浓度(O2 conc.)为100%。吸入率的参数包括:吸气时间(Ti) 0.90 s,吸气暂停时间(T pause)为0.00 s,吸气上升时间(T insp.rise)为0.15 s。触发流量(Trigg.Flow)为5。
![](Images/Table_Tmp.jpg)
Table 1. Accuracy of fault monitoring for three methods under different signal to noise ratios
表1. 三种方法在不同信噪比下故障监测准确率
由表1可以看出,随着信噪比的降低,所提方法、文献 [3] 方法和文献 [4] 方法的呼吸机管路积水故障监测准确率均有所下降,但所提方法明显高于文献 [3] 方法和文献 [4] 方法。在信噪比为−10 dB时,三种方法的故障监测准确率均达到最低值,但所提方法的故障监测准确率仍保持在97%以上,而文献 [3] 方法和文献 [4] 方法的故障监测准确率已低于90%,说明所提方法在噪声条件下具有更高的故障监测精度。因为所提方法在将采集到的呼吸机管路信号输入监测模型前通过小波变换去噪信号,滤除信号中噪声部分,将剩余有意义信号输入监测模型,从而有效提升故障监测准确率。
为了进一步验证所提方法、文献 [3] 方法和文献 [4] 方法的故障监测能力,将包含六种不同故障类型的实验样本信号特征降维至三维空间后作为三种方法的输入值,对比三种方法的故障监测能力,实验系统及实验结果如图5~7和图8所示。
![](//html.hanspub.org/file/136-2571684x153_hanspub.png?20240604170829550)
Figure 5. Screenshot of the experimental monitoring system
图5. 实验监测系统截图
![](//html.hanspub.org/file/136-2571684x154_hanspub.png?20240604170829550)
Figure 6. Reference [3] method fault monitoring capability
图6. 文献 [3] 方法故障监测能力
![](//html.hanspub.org/file/136-2571684x155_hanspub.png?20240604170829550)
Figure 7. Reference [4] method fault monitoring capability
图7. 文献 [4] 方法故障监测能力
![](//html.hanspub.org/file/136-2571684x156_hanspub.png?20240604170829550)
Figure 8. Fault monitoring capability of the proposed method
图8. 所提方法故障监测能力
其中图5实验监测系统为监护仪界面,它显示了患者信息以及时间、病床号等基本信息。同时波形图区域展示的是患者的生命体征,包括心电图波形(绿色)、血氧饱和度波形(蓝色)、血流波形(红色)以及呼吸波形(黄色)等。这些波形能够直观地反映患者心脏活动、血氧供应、血流情况以及呼吸节律等生理状态。其中,参数区是监护仪界面上最为关键的部分,它实时显示患者的各项生命体征数据,如心率(HR)、血压(BP)、呼吸频率(RR)、血氧饱和度(SpO2)、体温(TEMP)等。
由图4、图5和图6可以看出,文献 [3] 方法无法清晰区分类别1、类别3和类别5,使三种类别之间存在混淆,类别4和类别6之间也存在类样本间间距较小情况,且类别4和类别6的类样本内间距较大。文献 [4] 方法对类别1、类别2、类别3以及类别2、类别3和类别5的区分界限不够明显,类别之间存在明显的连接,类样本间间距过小,类别4和类别6的类内间距过大,且两者间也存在连接部分。而所提方法对六种故障类别的区分十分明显,类样本间距较大且类样本内间距较小。由此可知,所提方法未出现误判问题,具有较好的故障监测能力。
5. 讨论
本文研究了呼吸机管路积水故障自动监测算法,利用小波分析中的小波变换去噪方法,对呼吸机管路信号中的噪声进行去除,从而提升呼吸机管路信号可靠性和故障监测准确率。通过构建多输出支持向量机,结合改进的人工蜂群算法,对其相关参数进行优化,生成最优多输出支持向量机并加以训练。将去噪后信号输入至训练后多输出支持向量机中,输出故障监测结果,确保故障监测效果。
与文献 [3] 方法和文献 [4] 方法相比,本次研究利用小波变换去噪方法,去除了呼吸机管路中的噪声干扰,采用改进的人工蜂群算法,优化多输出支持向量机相关参数。而文献 [3] 方法和文献 [4] 方法仅对数据进行预处理,将收集到的信息放入训练好的模型中,容易存在监测误差。在不同信噪比下故障监测准确率测试中,文献 [3] 方法和文献 [4] 方法的故障监测准确率最高值分别为97.68%和96.81%,而所提方法的故障监测准确率最高值为99.81%,高于文献 [3] 方法和文献 [4] 方法,证明所提方法的故障监测效果更好。
本次研究仅对恺得VADI呼吸机进行监测,在今后的工作当中,需对不同型号的呼吸机进行监测,以便取得更加准确的监测结果。
6. 结论
为了解决目前呼吸机存在的监测精度较低,监测效果较差的问题,提出呼吸机管路积水故障自动监测算法研究。通过采用小波变换去噪技术,能够显著提高呼吸机管路积水故障的监测精度。这一改进有助于更准确地识别和预警潜在的故障,从而提升呼吸机使用的安全性。同时采用多输出支持向量机模型,并引用改进的人工蜂群算法对模型参数进行优化,使得模型能够更好地适应和识别各种故障信号,增强了其在实际应用中的可靠性。另外,该算法为呼吸机的更安全、更稳定使用提供了有力保障。有助于减少因故障导致的停机时间,提高设备的整体运行效率。
基金项目
上海市2021年度“科技创新行动计划”科学仪器和化学试剂项目申报指南/仪器共享配套操作与应用技术研究“大型多导睡眠监测仪模块拓展及共享配套应用技术研究”(21142203600)。
NOTES
*第一作者。
#通讯作者。