引言

心力衰竭是一种临床综合症，其产生的根本原因为心室泵血能力衰弱，导致每一次心脏搏动所排出的血液量无法满足人体新陈代谢的正常需要，因此也被称作为充血性心力衰竭。心力衰竭会导致人体其他组织和器官的血液灌注不足，正常功能受到影响，甚至是威胁患者生命安全。根据美国心脏病基金会(ACCF)以及美国心脏病协会(AHA)共同发布的充血性心力衰竭指南[1]，心力衰竭的等级一共分为四级，只有三、四级的患者才会出现明显的症状。在我国，心力衰竭的患病率已经达到了1.61%，大约一般的心力衰竭患者会在确诊后的五年内死亡[2]。因此，心力衰竭的自动诊断具有十分重要的临床价值和社会意义。

有研究表明，心力衰竭可以基于连续心率的量化指标进行检测[3]。该研究发现，连续心率的标准偏差(SDNN)可以被用于预测心力衰竭患者死亡的风险：SDNN 为50ms的受试分组死亡率(51.4%)约为SDNN 为100ms的受试分组死亡率(5.5%)的9.3倍。从此，不断有研究者尝试通过对连续心率的分析来识别心力衰竭患者。

近年来，深度学习算法不断被应用在图像识别、语音分析、自然语言处理等诸多领域[4,5]，并取得了显著的效果[6,7,8]。在众多深度学习模型中，深度残差网络(Residual Network，ResNet)[9]的表现尤其突出，迅速成为多种计算机视觉任务中最流行的网络框架之一。ResNet 允许原始输入信息跳过中间卷积层直接传输到后面的层中，使得神经网络学习的目标发生改变，从之前的完整输出转变成为输出与输入之间的差值，也就是残差。本文基于深度残差网络架构和短时连续心率信号进行心力衰竭的识别，有效提升心力衰竭检测的准确度和实时性，帮助心力衰竭患者进行自我管理，提升其生存率和生活质量。

1 心力衰竭自动诊断研究现状

Yu和Lee[10] 基于双谱分析方法对连续心率进行特征提取，并利用遗传算法进行特征选择，最终选取了约三分之一的特征输入到支持向量机(SVM)中进行心力衰竭的识别；Chen[11]等人从连续心率中提取了180个特征，其中包括54个经典统计特征和126个动态指标，利用非平衡决策树及SVM分类器进行心力衰竭的检测，最终结果达到了96.61%的准确度。

近年来，也有研究者基于短时连续心率数据进行各类心血管疾病的检测，其中也包括心力衰竭。研究表明[12]，短时连续心率可以用作检测心房颤动等心血管疾病，并且取得了良好的效果。在心力衰竭识别方法，短时连续心率分析也具有一定的诊断意义。Liu[13]等人提出一种基于模糊熵的方法对心力衰竭患者进行识别。之后，Liu[14]等人还尝试利用多尺度熵分析的方法对短时连续心率进行特征提取，并同样基于SVM分类器进行心力衰竭的诊断。Chen[15]等人基于稀疏自编码方法(Sparse auto-encoder)对短时心率信号进行特征提取，并利用深度学习网络对心力衰竭进行检测。

2 实验数据

3.FN(False negative)是假阴性样本数量，代表被错误划分为负例的个数(实际为正例，且被分类器划分为负例的样本数)；

BIDMC 心力衰竭数据库包括15位患有严重心力衰竭(NYHA等级为3～4级)患者的连续心电数据，其中包括11位年龄在22至71岁之间的男性以及4位年龄在54至63之间的女性；MIT-BIH 正常窦性心律数据库包括18位健康人群的连续心电数据，其中包括5位年龄在26至45岁之间的男性以及13位年龄在20至50岁之间的女性；而幻想曲数据库则包括20名年轻人(21到34岁)和20名老年人(68到85岁)经过严格筛查的健康受试者经过连续 2个小时的休息后测试的生理特征数据。在本研究中，我们首先根据受试者的档案信息将数据库分成训练集和测试集，这样做的目的在于保证测试集中的数据从未在训练阶段出现，防止过拟合对算法评估的影响。之后，将数据集中的连续心率数据以长度为500、1000和2000个连续心跳划分成等长的数据，表1总结了以上三种数据库在经过数据分割处理之后得到的样本数量。表2是测试集中受试者的信息(训练集受试者较多，在此不列出)。图1 表示的是在分割长度为500心跳时两个种类连续心率信号示意图。

3 深度网络结构

3.1 深度残差网络

深度残差网络是何恺明在2015年提出的一种深度网络模型，它对网络连接结构进行了改进，并在实际应用中取得了不错的效果。这种全新的网络连接方式也被称作为恒等快捷连接(Identify shortcut connection)，如图2所示，它直接跳过一个或多个层，从而将前若干层的数据输出直接引入到后面数据层的输入部分。

 

表1 不同数据库的样本量Table1 The number of signal for different database

  

分割长度500 1000 2000 BIDMC-CHF 3214 1607 803 NSR 3579 1739 869 FD 500 250 125数据库

 

表2 测试集中受试者信息Table2 Information of subjects in testing dataset

  

63, M, #13 38, F, #19140 61, M, #14 34, M, #19093受试者信息 分割长度CHF Normal 500 1000 2000 54, F, #11 50, F, #19830 686 339 164

  

图1 样本长度为500个心跳时两类数据Fig.1 RR intervals of CHF and normal subjects: a, CHF subjects; b, normal subjects

4.TN(True negative)是真阴性样本数量，代表被正确划分为负例的个数(实际为负例，且被分类器划分为负例的样本数)。

3.2 基于深度残差网络的心力衰竭诊断

整整一个月，西双没有回公司上班。他在度他的蜜月。在医院里在殡仪馆里度他的蜜月。西双焦头烂额。西双行色匆匆。日子兵荒马乱。生活水深火热。

在本研究中，我们采用准确性、灵敏性和特异性作为验证指标，具体的定义如下：

相关抵消对于信号估计是最佳线性处理方法，利用线性变换去掉信号x(N维)、y(M维)之间相关的部分.假设对y进行线性变换的矩阵是H，且y与x相关的部分为

统计结果显示，出版专著与教材主要是副高级及以上职称、45周岁以下教师，具有一定的科研能力，多数具有教研室主任或系部领导任职经历。

  

图2 残差构建块示意图Fig.2 A building block of residual learning

  

图3 LSTM 单元结构Fig.3 A LSTM cell with three gates

  

图4 深度网络架构Fig.4 Deep learning network structure

如图4所示，在本研究中，我们利用 LSTM 网络代替残差块之间的卷积网络。我们选择自适应矩估计方法(Adaptive Moment Estimation，adam)作为优化器，其参数设置为学习率为0.001, β1= 0.9以及β2=0.999。

信息评估准备阶段，根据水库影响区各乡镇居民实际信息，以及工程相关信息，制定移民安置规划目标。为了保证信息正确性及信息结构合理性，需要成立工程专家组，对移民安置评估对象进行资料收集整理。

4 实验结果及分析

4.1 评价指标

如图3所示，长短时记忆单元的原型是循环神经网络(Recurrent neural network, RNN)。它通过加入输入门、记忆门和遗忘门，来改进对长期时间序列时间信息的处理过程。LSTM 网络非常适合基于时间序列数据的分类、处理和预测任务，并且已经成功地应用在旅游时间预测和音乐生成等领域。

式中：cF为折减后的总黏聚力；φF为折减后的内摩擦角；kF为折减后的DP5总黏聚力；αF为折减后的DP5内摩擦角；Ftrial为折减系数。

 

其中，TP，FP，FN 和TN 的意义如下：

1.TP(True positive)是真阳性样本数量，代表被正确划分为正例的个数(实际为正例，且被分类器划分为正例的样本数)；

我想起回乡后第一次见到她，在半边街上骑一辆高档踏板车的情景，心里默默算计过去了多少日子，说，是啊，一晃又是好几个月过去了。

如上文所述，连续心率信号其本身是一种时间序列信号，它同时具有空间信息和时间信息。因此，我们在深度残差网络的基础上，引入长短时记忆(Long short-term memory, LSTM)单元对时间序列特征进行更有效的提取。

本研究的实验数据来源于国际公认的开源数据库PhysioNet，该数据库主要分为临床数据库和波形数据库两大类。其中，波形数据库包含多种非重症监护环境的生理信号的连续记录。在本研究中，我们采用了BIDMC 心力衰竭数据库(BIDMC Congestive Heart Failure Database，BIDMC-CHF)、MIT-BIH正常窦性心律数据库(MIT-BIH Normal Sinus Rhythm，NSR)和幻想曲数据库(Fantasia database，FD)。

①设计程序简单化。年度小农水项目投资确定后，省即下发通知编报实施方案，各项目县从接到通知至要求的编制完成上报时间一般为30~40天，由于时间紧迫，加之前期工作费用非常有限，勘察、测量等工作不能完全做到位，基础资料不足，设计依据不可靠，致使设计与实际存在脱节现象。如：部分井因前期未探明水文地质情况，在实施中为保证设计出水量，井深调整较多；桥涵设计跨径与其所在的渠沟断面不配套等。

深度残差网络结构的出现有效地解决了深度学习中的退化问题，这是因为残差网络学习的是残差函数F(x)=H(x)–x 而不是H(x)，尽管这两种形式的目标函数原则上都能够逼近所需的函数，但学习的难易程度并不相同。如果最优函数更趋近于恒等映射而不是零映射，那么对于求解器来说寻找关于恒等映射的干扰比学习一个新函数要容易的多。事实上，经过大量的实践验证，深度残差网络比其他网络，例如VGGNet和GoogLeNet 具有更高的准确率。

需要指出的是，在医学诊断中，灵敏性是指算法能够正确检测出患有某种疾病的患者的能力，而特异性则是代表算法能够正确识别没有疾病的人的能力。

2.FP(False positive)是假阳性样本数量，代表被错误划分为正例的个数(实际为负例，且被分类器划分为正例的样本数)；

4.2 实验结果

在实际的临床应用中，衡量一个模型是否成熟有效主要是看它在面对位置患者时，是否能做出准确可靠的检测。因此，在本研究中，我们采取无偏见测试(Blind test)对算法的有效性进行评估。表3是本方法在无偏见测试集上的表现。

 

表3 无偏见测试集上的表现Table3 Performance with blind test

  

信号长度 准确性 灵敏性 特异性500 99.67% 99.34% 100%1000 98.84% 97.53% 100%2000 96.63% 100% 93.64%

当输入的连续心率长度为500时，模型的识别整体准确度为99.67%，灵敏性为99.34%，特异性100%；当输入的连续心率长度为1000时，模型的识别整体准确度为98.84%，灵敏性为97.53%，特异性100%；而当输入的连续心率长度为2000时，模型的识别整体准确度为96.63%，灵敏性为100%，特异性93.64%。

5 讨论

与传统的心力衰竭识别算法相比，基于深度学习的心力衰竭识别算法能够在高维空间中进行可靠的特征提取，并且不需要人工操作，避免可能存在的人为误差。深度学习系统通过组合低级特征来发现数据的分布式特征表示，从而形成更加抽象的特征表示。

本方法具有两个优点，首先，基于残差神经网络进行心力衰竭的检测，基于该方法的决策系统可以从全部数据中自动获取有用的信息，而不是通过特征提取来进行人工的数据降维。这样做能够最大范围地保有数据的有用信息，避免潜在误差；其次，我们基于LSTM 网络单元对网络架构进行了修改，使得该模型更加适用于时间序列信号的分类。

然而，本研究也存在一些局限性。首先，本研究没有就不同严重程度的心力衰竭进行多分类识别，并且没有深入讨论心力衰竭药物对数据库中受试者连续心率的影响；其次，本方法需要大数据对模型进行训练，以获得最佳性能，并且基于LSTM修改后的网络模型在训练阶段的计算消耗有所增加。

6 总结

在文献 [9]的启发下，本文提出一种基于LSTM改进的残差网络，并将其应用在利用短期连续心率检测心力衰竭的研究中。基于3个开源数据库和4种不同长度的输入数据对所提出的方法进行了评估，结果表明，在长度为500、1000和2000时，本方法的准确度分别达到了99.67%、98.84%以及96.63%。

使用连续心率来进行心力衰竭检测对于医疗保健应用十分重要，特别是对于智能手表、手环等穿戴式设备。在下一步工作中，我们会将本模型部署到保健应用中，作为心力衰竭患者日常监控的辅助手段，同时尝试加入注意力机制进一步提高准确度。

参考文献

[1]Clyde W.Yancy, Mariell Jessup, Biykem Bozkurt, et al.2013 ACCF/AHA Guideline for the Management of Heart Failure: Executive Summary[J].Journal of the American College of Cardiology, 2013, 62(16): 1459-1539.

[2]Jun Huang, Hongjun Yin, Milun Zhang, et al.Understanding the economic burden of heart failure in China: impact on disease management and resource utilization[J].Journal of Medical Economics, 2017,20(5):5.

[3]Nolan J , Batin P D , Andrews R , et al.Prospective Study of Heart Rate Variability and Mortality in Chronic Heart Failure : Results of the United Kingdom Heart Failure Evaluation and Assessment of Risk Trial(UK-Heart)[J].Circulation, 1998, 98(15):1510-1516.

[4]Thakre T P, Smith M L.Loss of lag-response curvilinearity of indices of heart rate variability in congestive heart failure[J].BMC Cardiovascular Disorders, 2006, 6(1):27-0.

[5]Yoon K H, Thap T, Jeong C W, et al.Analysis of Statistical Methods for Automatic Detection of Congestive Heart Failure and Atrial Fibrillation with Short RR Interval Time Series.[C]In: International Conference on Innovative Mobile & Internet Services in Ubiquitous Computing.Blumenau: IEEE, 2015.

[6]Schmidhuber, Jürgen.Deep learning in neural networks:An overview[J].Neural Networks, 2015, 61:85-117.

[7]Potes C, Parvaneh S, Rahman A, et al.Ensemble of Feature-based and Deep learning-based Classifiers for Detection of Abnormal Heart Sounds[C]In: Computing in Cardiology Conference.Vancouver: IEEE, 2017.

[8]Hwang B, You J, Vaessen T, et al.Deep ECGNet: An Optimal Deep Learning Framework for Monitoring Mental Stress Using Ultra Short-Term ECG Signals[J].Telemed J E Health, 2018.

[9]He K, Zhang X, Ren S, et al.Deep Residual Learning for Image Recognition[C].In: IEEE Conference on Computer Vision and Pattern Recognitiaon.Las Vegas: IEEE, 2015.

[10]Yu S N, Lee M Y.Bispectral analysis and genetic algorithm for congestive heart failure recognition based on heart rate variability[J].Computers in Biology and Medicine, 2012, 42(8):816-825.

[11]Wenhui C, Lianrong Z, Kunyang L, et al.A Novel and Effective Method for Congestive Heart Failure Detection and Quantification Using Dynamic Heart Rate Variability Measurement[J].PLOS ONE, 2016, 11(11):e0165304-.

[12]Lake D E, Moorman J R.Accurate estimation of entropy in very short physiological time series: the problem of atrial fibrillation detection in implanted ventricular devices[J].Am J Physiol Heart Circ Physiol, 2011, 300(1):H319.

[13]Liu C, Li K, Zhao L, et al.Analysis of heart rate variability using fuzzy measure entropy[J].Computers in Biology and Medicine, 2013, 43(2):100-108.

[14]Liu C, Gao R.Multiscale Entropy Analysis of the Differential RR Interval Time Series Signal and Its Application in Detecting Congestive Heart Failure[J].Entropy, 2017, 19: 251-270.

[15]Chen W, Liu G Z, Su S W, et al.A CHF Detection Method based on Deep Learning with RR Intervals[C].In: 39th Annual International Conference of the IEEE Engineering in Medicine and Biology Society.Seogwipo: IEEE, 2017.