CN106551691B

CN106551691B - 一种心率变异性分析方法、装置及用途

Info

Publication number: CN106551691B
Application number: CN201611094149.6A
Authority: CN
Inventors: 李路明; 刘洪运; 杨曌
Original assignee: Tsinghua University; Beijing Pins Medical Co Ltd
Current assignee: Tsinghua University; Beijing Pins Medical Co Ltd
Priority date: 2016-12-02
Filing date: 2016-12-02
Publication date: 2020-01-21
Anticipated expiration: 2036-12-02
Also published as: US10973455B2; US11793449B2; US20210219906A1; US20190099121A1; CN106551691A; WO2018099120A1

Abstract

本发明涉及一种心率变异性(Heart Rate Variability，HRV)分析的方法、装置及用途。利用低成本、便携、可穿戴式信号采集设备获取癫痫患者的术前24小时心电(Electrocardiography，ECG)信号，通过程式化的HRV分析方法计算长程、短程的ECG的时域、频域和非线性指标，并基于表征迷走神经对心率调制作用程度即迷走神经活性的特征参数，准确、高效地对迷走神经刺激(Vagus Nerve Stimulation，VNS)手术的药物难治性癫痫患者进行疗效预测，避免不必要的经费支出和避免耽误最佳治疗时机，同时通过ECG的HRV特征参数明确选择VNS手术适应症患者，可以从整体上提高VNS疗法的疗效。

Description

一种心率变异性分析方法、装置及用途

技术领域

本发明涉及一种心率变异性分析方法、装置及用途，尤其涉及心率变异性分析方法及其装置在迷走神经刺激(Vagus Nerve Stimulation，VNS)治疗药物难治性癫痫领域的疗效预测。

背景技术

癫痫作为一种疾病会影响患者的生活，大部分患者可以用一种或多种药物联合治疗对发病加以控制，但仍然会有部分患者对药物治疗并不敏感，这些患者被称为药物难治性癫痫患者。迷走神经刺激(Vagus Nerve Stimulation，VNS)作为一种辅助治疗手段可以有效控制药物难治性癫痫患者的癫痫发作，但其疗效的个体差异非常大，而且具有很高的不确定性。统计分析结果显示，仅有5％-9％接受VNS手术的药物难治性癫痫患者癫痫发作得到完全控制，另有约10％的患者完全无效，其余表现为不同程度的发作频率减少。总体而言，50％-60％接受VNS手术的药物难治性癫痫患者可以达到发作频率减少50％的治疗效果。针对VNS治疗药物难治性癫痫疗效的不确定性和个体差异大的问题，研究人员一直试图通过术前评估筛选出适合VNS手术的患者。

现在尚无明确的、量化的VNS治疗药物难治性癫痫的疗效预测方法应用于临床，而基于脑电(Electroencephalography，EEG)、核磁影像资料(Magnetic Resonance Imaging，MRI)、患者人口统计学特征(性别、年龄等)、临床病史(病程)、发作特征(包括发作类型、发作频率、病灶位置等)进行VNS疗效相关性因素的研究，结论也不相一致，甚至相互矛盾。

现有技术中的思路如附图1所示：药物难治性癫痫患者在接受VNS手术之前首先需进行系统、全面的术前评估(包括人口统计学特征、病史、发作特征、MRI、EEG等)，然后进行VNS植入手术，术后2周左右开机，之后根据每位患者对VNS急性刺激的具体反应情况按照VNS产品的性能、技术特点逐步调整刺激参数并定期统计所有患者发作减少情况。在相应的随访周期(一般为1年)结束后，根据患者发作减少情况即疗效对患者进行分类，最后对不同疗效分类患者的术前评估数据进行统计学分析，寻找不同疗效患者组之间有统计学差异的参数作为VNS治疗药物难治性癫痫疗效预测的敏感因子。

现有的EEG和MRI的方法存在成本高、操作复杂和对分析人员专业知识水平要求高的缺点，总体而言，基于上述方法的研究都没有在临床得到应用，主要原因是对同一问题的研究结论不相一致，甚至相互矛盾。本发明提出的基于24小时动态心电信号的方法，仅需要给被试佩戴一个便携式动态心电记录盒，在被试自由活动的状态下采集ECG，操作简单。一般医院门诊对24小时动态心电检查的收费为240元，与长程视频脑电及核磁共振影像检查1000元左右的费用相比，成本相对较低。最为重要的是24小时动态心电采集不受活动限制，与EEG和MRI相比，采集相对简单，一致性较好。

心率变异性(Heart Rate Variability,HRV)是指心电(Electrocardiography，ECG)信号中的相邻心搏间期随时间的变化，它起源于自主神经系统对窦房结自律性的调制，使心搏间期存在几十毫秒甚至更大的差异或波动。HRV蕴含着神经体液调节的大量信息，是目前评价心血管系统自主神经系统活性及其调节功能的定量、无创、可重复指标，通过HRV分析可间接反映交感神经及副交感神经之间的相互作用。

目前还没有利用心电信号HRV分析技术对VNS治疗药物难治性的疗效进行预测的研究和技术方案。由于癫痫疾病与心脏自主神经功能失调有着密切的联系。伴随癫痫的发病和进展，患者心脏自主神经系统平衡被打破，一般表现为交感神经活性增强和迷走神经活性降低。基于这个结论，本发明考虑利用低成本、便携、可穿戴式信号采集设备获取癫痫患者的术前24小时ECG信号，通过程式化的HRV分析方法计算ECG的时域、频域和非线性特征参数，并选取表征迷走神经活性的特征参数，准确、高效地对VNS治疗药物难治性癫痫的疗效进行预测。

发明内容

申请人发现，癫痫疾病与心脏自主神经功能失调有着密切的联系。伴随癫痫的发病和进展，患者心脏自主神经系统平衡被打破，一般表现为交感神经活性增强和迷走神经活性降低。而且，还发现其中心率变异性(Heart Rate Variability,HRV)这一指标非常重要，HRV是指心电(Electrocardiography，ECG)信号中的相邻心搏间期随时间的变化，它起源于自主神经系统对窦房结自律性的调制，使心搏间期存在几十毫秒甚至更大的差异或波动。HRV蕴含着神经体液调节的大量信息，是目前评价心血管系统自主神经系统活性及其调节功能的定量、无创、可重复指标，通过HRV分析可间接反映交感神经及副交感神经之间的相互作用。

基于这一研究发现，本发明利用低成本、便携、可穿戴式信号采集设备获取癫痫患者的术前24小时ECG信号，通过程式化的HRV分析方法计算长程、短程ECG的时域、频域及非线性的特征参数，并通过表征自主神经系统中迷走神经活性的特征参数，准确、高效地预测VNS治疗药物难治性癫痫的疗效，避免不必要的经费支出和避免耽误最佳治疗时机，同时也可以利用ECG的HRV分析得到的迷走神经活性特征参数明确选择VNS手术适应症患者，可以从整体上提高VNS疗法的疗效。

一种心率变异性分析方法，包括如下步骤：

1)体外采集心电数据；

2)对心电数据进行数字化、去噪处理；

3)将处理后的心电数据形成窦性NN间期序列；

4)分别选取长程、清醒、睡眠状态下短程的窦性NN间期数据；

5)对长程、清醒、睡眠状态下各的窦性NN间期数据进行时域、频域、非线性指标的计算(各项指标的具体计算见下表)；

表1时域指标

注：RR为相邻两次心跳之间的时间间隔或间期；N为RR间期的总数；RR_i为第i个RR间期。

频域指标：功率谱分析得出，典型的HRV功率谱上有三个分离的谱峰,大致集中于0.04Hz、0.1Hz和大于0.15Hz频段，将各频段下的面积可作为功率谱能量的测定值，即为HRV的频域指标。频域指标反映了自主神经在各频率段对心率的调制，而且是将各种生理因素作适当的分离，因而能准确反映HRV的真实情况，临床应用价值较大。

表2频域指标

非线性指标：近年来非线性动力学的研究取得了很大的进展，而且在物理学、化学、信息科学领域中引起了人们很大的兴趣。用非线性动力学方法来描述和分析自主神经系统对心率调控这一复杂过程，可以获得更多的有意义信息，以提高HRV作为预测指标的准确性。本发明采用的HRV非线性指标主要是庞卡莱散点图(Poincare Plot)的SD1和SD2指标(如图4所示)。以NN间期作为研究分析对象，并采用坐标形式显示前后两个NN间期的变化特征散点。分别测定散点图中心点处的长轴及短轴上的标准差，以短轴为轴心测定的标准差称为SD1，反应瞬时的心搏与心搏间差异，以长轴为轴心测得的标准差称为SD2，则反映长时间内的NN间期的差异性。

附图说明：

图1是现有技术筛选VNS适应症患者流程图

图2 12导联ECG采集示意图

图3 ECG信号具体处理流程图

图4 HRV非线性指标SD1和SD2计算示意图

图5判断流程图

图6响应组Responder₅₀和非响应组Non-responder₅₀的参数图

图7受试者工作特征曲线

具体实施方式：

实施例1

如图2所示，于术前的24小时标准12导联ECG采集：要求心电采集装置的采样率要大于或等于500Hz，ECG记录期间避免剧烈运动、服药等可能影响心脏功能的活动，记录时间长度为24小时。受试对象、个体的记录环境和条件应基本类似。并要求确保用于HRV分析的数据是正常窦性NN间期。在进行HRV分析时，提取24小时时域、频域和非线性指标的同时，提取清醒、睡眠状态下各5min短程心电的时域、频域和非线性指标。ECG信号具体处理流程如图3所示。

1)采集ECG信号，并对信号数字化处理；

2)对数字信号进行去噪声、去伪迹处理；

3)对其中的QRS波自动检测

4)对检测后的信号进行QRS波人工检视；

5)再对异位起搏QRS波信号进行剔除；

6)形成窦性NN间期序列；

7)分别选取被试24小时长程、清醒、睡眠状态下各5min短程的窦性NN间期数据；

8)按照表1、表2及图4的方法对被试24小时长程、清醒、睡眠状态下各5min短程的窦性NN间期数据进行时域、频域、非线性指标的计算；

9)选取表征迷走神经活性的特征参数RMSSD、pNN50、SD1和HF。

对药物难治性癫痫患者进行术前24小时心电采集，采集到的24小时心电数据按上述处理方法得到24小时长程和5min清醒、睡眠状态下的短程正常窦性NN间期序列。对上述长程、短程NN间期序列按表1、表2和图4的方法进行HRV时域、频域和非线性分析并选取表征迷走神经活性的特征参数RMSSD、pNN50、SD1和HF。通过相应的阈值判断进行综合判断选择，进而对VNS治疗药物难治性癫痫的疗效进行预测判断(如图5所示)。利用表征迷走神经活性的夜间5min的特征参数RMSSD、pNN50、SD1和HF分别单独进行VNS治疗药物难治性癫痫的疗效预测，其对应阈值的选择及其相应的预测准确率如下：

RMSSD＝37ms时，患者术前ECG分析结果大于该值，其VNS术后发作改善将≥50％，其预测的准确率为88.2％，同时也可以作为筛选患者的依据；

pNN50＝27％时，患者术前ECG分析结果大于该值，其VNS术后发作改善将≥50％，其预测的准确率为82.7％，同时也可以作为筛选患者的依据；

SD1＝35ms时，患者术前ECG分析结果大于该值，其VNS术后发作改善将≥50％，其预测的准确率为70.6％，同时也可以作为筛选患者的依据；

HF＝798ms²时，患者术前ECG分析结果大于该值，其VNS术后发作改善将≥50％，其预测的准确率为76.5％，同时也可以作为筛选患者的依据；

实施例2

实施例1的HRV分析方法中其他可能表征迷走神经活性的标亦可用于上述VNS治疗药物难治性癫痫的疗效预测和患者筛选。

本发明通过药物难治性癫痫患者术前24小时ECG的采集和HRV的时域、频域及非线性分析，可以在VNS术前对药物难治性癫痫患者的疗效进行预测，指导药物难治性癫痫患者正确选择VNS手术以避免不必要的经费支出和避免耽误最佳治疗时机，同时通过ECG的HRV分析得到的表征迷走神经活性的特征参数进行VNS疗效的预测及明确选择VNS手术适应症患者，可以从整体上提高VNS疗法的疗效。

实施例3

根据上述筛选方法，选取于2014年8月13日至2014年12月31日期间在北京天坛医院完成VNS手术的32例药物难治性癫痫患者进行验证。32例药物难治性癫痫患者在VNS术前进行了完整的评估(包括人口统计学特征、临床病史、抗癫痫药物史、24小时视频脑电、MRI以及24小时动态心电等)。

通过术前24小时动态心电数据按上述ECG信号处理方法进行HRV的时域、频域和非线性分析，提取了每位患者表征迷走神经活性的特征参数RMSSD、pNN50、SD1和HF。术后1年随访结束时，32例接受VNS治疗的药物难治性癫痫患者中，有17例患者的发作频率减少50％或以上(其中6例患者发作完全得到了控制)视为响应组Responder₅₀，15例患者的发作频率减少50％以下(其中4例患者完全无效)视为非响应组Non-responder₅₀。响应组Responder₅₀和非响应组Non-responder₅₀的RMSSD、pNN50、SD1和HF如图6所示，两组患者的上述表征迷走神经活性的参数存在显著统计学差异，提示VNS疗效和患者术前迷走神经活性有关，同时也说明通过HRV分析方法可以进行VNS的疗效预测。进一步对24小时长程及清醒、睡眠状态下各5min短程的RMSSD、pNN50、SD1和HF的疗效预测能力进行评估，通过ROC曲线分析，结果表明上述4个参数中睡眠状态下5min的RMSSD预测的最为准确(如图7所示)：将其阈值设定为37ms时，17例响应组的患者中仅有2例患者的RMSSD低于37ms，其余患者的均大于37ms，预测的准确率超过85％，从而证明上述HRV分析的获取的表征迷走神经活性的特征参数可以准确、有效的进行VNS疗效的预测及VNS适应症的筛选。

以上所述仅是本发明的优选实施方式，应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明原理的前提下，还可以做出若干改进和润饰，这些改进和润饰也应视为本发明的保护范围。此外，尽管本说明书中使用了一些特定的术语，但这些术语仅仅是为了方便说明，并不对本发明构成任何限制。

Claims

1.一种心率变异性分析方法，包括如下步骤：

1)体外采集心电数据；

2)对心电数据进行数字化、去噪处理；

3)将处理后的心电数据形成窦性NN间期序列；

4)分别选取被试20小时以上长程、清醒、睡眠状态下各2-10分钟短程的窦性NN间期数据；

对被试者20小时以上长程、清醒、睡眠状态下各2-10分钟短程的窦性NN间期数据进行时域、频域、非线性指标至少之一的计算，所述计算结果用于设定药物难治性癫痫患者术前疗效阈值；

其中时域指标计算公式如下，分析时选取其中一种或多种进行计算：

频域HF指标计算公式如下：

指标名称指标定义单位频率范围 HF 高频段的功率 ms<sup>2</sup> 0.15～0.4Hz

非线性指标为庞卡莱散点图的SD1指；所述RMSSD阈值为37ms、pNN50阈值为27％、SD1阈值为35ms、HF阈值为798ms。

2.如权利要求1所述的心率变异性分析方法，其特征在于：长程时间为24小时，短程时间为5分钟。

3.如权利要求1或2所述的方法，其特征在于：用于VNS治疗药物难治性癫痫的分析。

4.一种心率变异性分析装置，其特征在于包括一计算模块，所述计算模块采用如权利要求1或2所述的方法进行计算。

5.如权利要求4所述心率变异性分析装置，其特征在于还包括数据采集模块、数字化处理模块、去噪模块、判断模块中的一个或多个。

6.如权利要求4或5所述的装置，其特征在于：用于VNS治疗药物难治性癫痫的分析。

7.一种可穿戴心电监测设备，其特征在于包含权利要求4-6之一所述的心率变异性分析装置。