具体实施方式
如图1所示,本发明技术实施例中,包括睡眠状态检测装置102或者具有控制器104的设备,该控制器具有一个或多个处理器可以用来执行睡眠状态和/或睡眠觉醒的检测方法,在此,给予这样的计算法进行详细的描述。在一个实施例中,该控制器系统含有多个元件或模块,可以控制控制器的多个方面。所述组件或模块可连接有整合的芯片,存储卡和/或者其他控制指令的数据或信息储存方式。例如,包含这种检测方法的程序指令可编码在集成芯片上,作为装置或设备记忆以形成一个特殊用途的集成芯片(ASIC)。这些指令可以或者可选地,以适当的数据存储介质以软件或硬件的形式装载以控制一个或多个程序控制器。
图1的实施例中,所述控制器104可使用呼吸流量信号或其他包括任选的流量测量模块105的数据。通过这些模块,所述控制器可以直接测量呼吸流量。因此, 处理器可以控制睡眠状况(此处详细描述的基于来自先前的睡眠区间的上述呼吸流量数据)的检测或评估。可选择的,所述检测可在同时期的睡眠区间执行,通过检测到呼吸流量信号,使用流量传感器的目前检测的流量数据。因此,在一些实施例中,使用执行方法,根据流量传感器检测的呼吸流量信号,该装置本身随机地被执行。例如,到达或经过鼻管或面罩的流量可被测量,通过使用呼吸速度描记器和不同的压力变换器或者类似装置如软管或输气管的应用包来得到流量信号。选择性地,流量信号可通过其他传感器而被推断,例如2005年11月2号提交的PCT/AU2005/001688 专利申请中描述的马达传感器(motor
current sensor),通过交叉参考,其完全公开的技术在此与本发明结合。
在此,进行更详细的讨论,流量数据可以通过模块处理例如流量限制检测模块106,窒息和/或睡眠不足(AHI)检测模块108,呼吸流量特性或特征提取模块110和泄漏检测组件。所述AHI模块108可以利用睡眠状态检测模块114输入的数据来判定窒息和/或睡眠不足的情况。从流量限制检测模块106,AHI检测模块108和呼吸特征提取模块110输出的信号或数据,可通过觉醒检测模块114加工。同样地,从泄漏检测组件和呼吸特征提取模块110输出的数据可以通过睡眠状态检测模块118处理。进一步的,觉醒检测模块114和病人特征数据存储库112输出的数据可被睡眠稳定性检测模块116处理。
睡眠质量检测/ 评估模块120可处理从睡眠稳定性检测模块116,病人特征数据存储模,块122和睡眠状态检测模块118输出的数据和信号。睡眠质量评估模块120输出的信号和数据可通过用户接口模块124显示和/或者输入于反馈控制模块126,用于修正或者调整治疗控制处理模块128的装置置,例如压力治疗方案装置。
在本技术的一些实施例中,所述检测仪和控制器的模块具备如下特征和功能。
A.病人信息存储模块
该模块可以控制使用睡眠状态检测装置102的病人的数据和信息的存储。这些数据可以通过用户接口模块124界面输入,任意的,该模块包括键盘,面板开关,触摸屏等。例如,在所述模块的控制下,该装置可促进病人或者医生选择或者登陆一个或多个病人信息例如:年龄,性别,体重,身高,BMI和/或者早前出现的情况信息(如,影响心率改变的潜在健康情况,可能导致空气流动模式变化的潜在健康状况,或者其他目前的生理状况等)。这样的健康状况可能包括例如哮喘、肺气肿、慢性心血管疾病。储存信息反过来可以用于所述装置的其他模块。
这些病人的信息可以用来对所述方法学的输出值的加权。例如,可能不存在绝对的空气流量的特征,可以用来作为适应于每一个病人检测不同睡眠状态的标准。因此,病人的信息可以用来调整特征提取模块的敏感度(或者呼吸特征的临界值被比较。这样,对于不同的病人,空气流量特征被更加精确地评估。例如,一个年龄超过50岁的病人的上气道的稳定性就很差。因此,对于这样的病人,在评估气流稳定性时就应该选择本技术中不太剧烈的系统或者方法。
在此,本系统所得病人信息的其他用途也被详细介绍。
B.流量限制检测模块
所述呼吸流量限制检测模块106可通过流量信号数据以检测呼吸流量限制。例如,所述装置可设置用于检测流量限制,例如在2008年5月9日提交的专利申请PCT/AU08/000647中公开的模糊的 流量限制测量(国际专利申请公开号 No.
WO/2008/138040)和2007年4月17日提交的,美国临时专利申请号No. 60/965,172。选择性地,其可用来通过分析流量数据的平缓(flattening)来检测流量限制。例如,所述平缓可以通过美国专利号. 5,704,345的专利中公开的方法判定,其全部的公开通过交叉引用与本发明结合。
C.AHI检测模块
在一些实施例中,流量数据可被用于为睡眠不足和/或者呼吸暂停的检测或者评分(如,阻塞性呼吸暂停和/或者中枢性呼吸暂停)。例如,通过美国专利号Nos . 6,138,675和6,029,665的专利中所述方法(其全部公开的技术通过交叉引用在此被加入),所述装置可以用来检测阻塞性呼吸暂停,部分阻塞性呼吸暂停和/或者中枢性呼吸暂停。
在一些实施例中,所述AHI检测模块108可以被改良用于拒绝基于泄漏检测组件112的评分事件.例如,如果检测到泄漏,该模块可以禁止任何同时期的评分事件,直到检测不到泄漏。这样,相对于在整个带面罩的时间内,AHI可被计算。也即是,这个时间是病人实际上带面罩的时间而不是病人总的睡眠时间。这可以避免在流量发生器(FG)执行过程中的过分评价。
因此,在一些实施例中,选择性地,AHI检测模块可以连接平缓显示器和泄漏显示器以提供简单的睡眠质量打分,这样,可以基本地洞察睡眠质量。
D.泄漏检测模块
相应地,在一些实施例中,流量数据可用来检测泄漏的存在,这样的泄漏与脱去面罩泄漏或高泄漏(high
leak)相关。例如,高泄漏可认为是超过某个临界值的量的泄漏(如,超过0.4米/秒(l/s) 的泄漏)。进一步举例,例如,所述装置可用来检测正在发生的泄漏或已经发生的泄漏,通过美国的专利号Nos .
6,152,129 和/或
6,532,957的专利公开方法,其全部公开的技术通过交叉引用在此被加入。选择性地,泄漏检测可用来帮助因面罩泄漏导致的从睡眠中觉醒的分类,正如前面提到的因泄漏的存在而阻止AHI的打分。
E.特征提取模块
一个或多个流量特性须通过特征提取模块110确定。在典型实施例中,特征可通过处理呼吸流量数据来测定。选择性地,特征可在呼吸-呼吸的基础上或者通过滑动窗口比较几个呼吸周期来计算。所述特征可用来作为该系统的其他模块的指示或者输入。例如,根据一个或者几个检测的特征的数据,在其他模块根据这些流量数据可得出关于睡眠状况的各种结论。以下是可通过流量数据被检测或测定的特征的举例。
1.吸气峰流量的变异度
所述呼吸流量特征可通过计算吸气峰流量的变异数(variance)来测定。所述特征在图2中有详细介绍。所述变异数可用包括多次呼吸(如5次呼吸)的滑动窗口测定。如图2所示,所述计算是采用一组连续呼吸的吸气峰流量值202来计算。当采用具有合适的临界值的状态检测仪来评定,该值可显示不同的睡眠状态,如下所示:
(a) 觉醒状态:显示高水平的吸气峰流量的变异度。
(b) 快速动眼睡眠状态:显示中等水平的吸气峰流量的变异度。
(c)非快速动眼睡眠状态:显示低水平的吸气峰流量的变异度。
此处吸气峰流量的变异度的特征可标记为"IPFV"。
2.呼气峰流量的变异度
该呼吸流量特征可通过计算呼气峰流量的变异度来确定。所述特征在图2中有详细介绍。所述变异度可用包括多次呼吸(如5次呼吸)的滑动窗口测定。如图2所示,所述计算是采用一组连续呼吸的呼气峰流量值204来计算。当采用具有合适的临界值的状态检测仪来评定,该测定值可显示不同的睡眠状态,如下所示:
(a) 觉醒状态:显示高水平的呼气峰流量的变异度。
(b) 快速动眼睡眠状态:显示中等水平的呼气峰流量的变异度。
(c)非快速动眼睡眠状态:显示低水平的呼气峰流量的变异度。
在此,呼气峰流量的变异度特征可标记为"EPFV"。
3.呼气峰流量位置/呼气的时间之比
如图3所示,该呼吸流量特征可被计算,按照从呼气开始到达呼气峰流量(TPE)所需时间与总的呼气时间(TE)之比。该比值如图3所示。当采用具有合适的临界值的状态检测仪来评定,该比值显示不同的睡眠状态,如下所示:
(a) 觉醒状态:相比于睡眠状态显示一个较大的呼气峰流量位置/呼气的时间之比。
(b) 快速动眼睡眠状态:相比于觉醒状态显示小的呼气峰流量位置/呼气的时间之比。
(c)非快速动眼睡眠状态:相比于觉醒状态显示低的呼气峰流量位置/呼气的时间之比。
4. 呼气峰流量位置/呼气的时间之比的变异度
该呼吸特征可按照呼气峰流量位置/呼气的时间之比的变异度(例如变异性)来计算。所述变异度可用包括多次呼吸(如5次呼吸)的滑动窗口测定。当采用具有合适的临界值的状态检测仪来评定,该值显示不同的睡眠状态,如下所示:
(a) 觉醒状态:从一次呼吸到另一次呼吸的显示高的呼气峰流量位置/呼气的时间之比的变异度。
(b) 快速动眼睡眠状态:从一次呼吸到另一次呼吸的显示低的呼气峰流量出现的位置/呼气时间之比的变异度。
(c)非快速动眼睡眠状态:从一次呼吸到另一次呼吸的显示低的呼气峰流量出现的位置/呼气时间之比的变异度。
5.低于75%的呼气峰流量的面积
该呼吸特征可按照呼气峰流量的呼气流量曲线的面积的一定比例值(如75%-100%之间)来计算。图3列举了呼气峰流量面积306例子。当采用具有合适的临界值的状态检测仪来测,该值可显示不同的睡眠状态,如下所示:
(a) 觉醒状态:在低于75%的呼气峰流量时显示大的面积值。
(b) 快速动眼睡眠状态:在低于75%的呼气峰流量时显示小的面积值。
(c)非快速动眼睡眠状态:在低于75%的呼气峰流量时显示小的面积值。
6. 低于75%的呼气峰流量的面积的变异度
该呼吸流量特性可通过对前述特征的面积的变异度(例如方差)计算。所述变异度可用包括多次呼吸(例如5次呼吸)的滑动窗口(sliding window)测定。当采用具有合适的临界值的状态检测仪来评定,该值显示不同的睡眠状态,如下所示:
(a) 觉醒状态:从一次呼吸到另一次呼吸的低于75%的呼气峰流量的面积显示大的变异度。
(b) 快速动眼睡眠状态(REM Sleep State):从一次呼吸到另一次呼吸的低于75%的呼气峰流量的面积显示小的变异度。
(c)非快速动眼睡眠状态(NREM Sleep State):从一次呼吸与另一次呼吸的低于75%的呼气峰流量的面积显示小的变异度。
7. 从75%呼气峰流量到开始吸气的时间
该呼吸流量特征可以按照呼气峰流量的上升部分的流量洪峰的一定比例(例如75%)和下一次开始吸气之间所用的时间来计算。这一特征在图4中得以说明。当采用具有合适的临界值的状态检测仪评估,该值显示不同的睡眠状态,如下所示:
(a)觉醒状态:从75%呼气峰流量到开始吸气会显示一个较缓慢的上升时间。
(b)快速动眼睡眠状态:从75%呼气峰流量到开始吸气会显示一个较快的上升时间。
(c)非快速动眼睡眠状态:从75%呼气峰流量到开始吸气会显示一个较快的上升时间。
8.从上次确认呼吸以来的时间
该呼吸流量特征,可以按照自上一次确认吸气以后的时间来计算。该特征参见图4。所述特征可被具有合适的临界值觉醒检测仪进行评定。
9.从上次确认呼吸以来的时间的变异度
该这种呼吸流量的特征,可以按照上述特征的自上次确认吸气之后的时间的变异度(如方差)进行计算。所述变异度可用包括多次呼吸(如5次呼吸)的滑动窗口测定。当采用具有合适的临界值的状态检测仪评定,该值显示不同的睡眠状态,如下所示:
(a)觉醒状态:从上次确认的呼吸以来,从一个呼吸到另一个呼吸所花的时间显示巨大的变异度。
(b)快速动眼睡眠状态:从上次确认的呼吸以来,一个呼吸到另一个呼吸所花的时间显示中等的的变异度。
(c)非快速动眼睡眠状态:从上次确认的呼吸以来,一个呼吸到另一个呼吸所花的时间显示很小的变异度。
从上次确认呼吸以来时间变异度的特征在此表示为TSLBV。
10.高频率呼吸周期
这种呼吸流量特征可以通过在滑动呼吸时期(如15次)之上扫描任何小部分高频率呼吸来计算。该特征可参见图5。当采用具有合适的临界值的状态检测仪来评定,这个值可以指示快速动眼睡眠的呼吸。
11.吸气时间变异度
这种呼吸流量特征可按照吸气时间的变异度(如方差)被计算。所述变异度可用包括多次呼吸(如5次呼吸)的滑动窗口测定。当采用具有合适的临界值的状态检测仪评定,该值显示不同的睡眠状态,如下所示:
(a)觉醒状态:吸气时间时显示巨大变化
(b)快速动眼睡眠状态:吸气时间时显示中度变化
(c)睡眠状态:吸气时间时显示微小变化。
12.呼吸一致性检测
该特征可以用来测定多次呼吸(如5次)期间的从呼吸到呼吸的呼吸流量特征是否一致。所述一致性考虑到吸气时间,吸气峰流量位置,呼气时间,呼气峰流量位置。所述这些特征的细微差别可以认为是一致的。
13.特征一致性检测
该特征对在多个呼吸时期(如30次呼吸)的上述确定的呼吸流量特征(如特征1—12)的变异度进行检验。例如,通过多个呼吸(如30次)期间的任何特征的自动分析,可以识别任何特征是否一致的或者是否为有意义的变异度。如,在清醒期间,流量峰的变异度特征平均起来,比在睡眠期间要高。每次计算得到的这样一个值,如通过30次呼吸的期间得到,可能不会比觉醒状态的临界值高,但是在评估的临界值以上的部分是有意义的。
14.达到95%吸气峰流量的时间
该呼吸流量特征可以按照从开始吸气到呼吸达到吸气峰流量的一定比例(如95%)所用的时间来计算。该95%的时间的特征,如图6所示,用 T_95表示。所述这种呼吸流量特征通过唤醒检测模块144在评估中被用到。
15. 高于吸气峰流量的7 5%峰面积
该呼吸流量特征按照高于吸气峰流量的一定比例(如75%)的吸气峰流量曲线的面积而计算得到。该呼吸流量特征在评估中可用于觉醒检测模块114。75%吸气峰时间特征(T_75)是指一次一口气达到一个吸气峰流量的比例(如95%),这可以通过计算图6中阴影部分的面积而确定。所述高于7 5%吸气峰流量的峰面积的特征在本文中简写为“IP75A”。
16. 高于7 5%吸气峰流量的峰面积的变异度
这种吸气流量特性可以计算上述吸气流量特征的变化。所述这种变化用包括多次呼吸(如5次)的一个滑动窗口来确定。所述吸气流量特性的特点在评估中可以用于觉醒检测模块114。所述高于7 5%吸气峰流量的峰面积的变化的特征在本文中简写为“I75AV”。
17.
95%吸气流量最峰和呼气峰流量之间的时间
该呼吸流量特性可以按照吸气峰流量曲线上升侧达到一定比例(如95%)到呼气峰流量的位置之间的时间来计算。该呼吸流量的特征在评估中用于醒检测模块114。
18. 当前的呼吸3mvTtot与先前呼吸3mvTtot之比
3mvTtot比率可按照这次呼吸的3mvTtot除以上次呼吸的3mvTtot来计算。3mvTtot由平均每分钟通气量(如三分钟通气量)和总呼吸时间(例如,当前呼吸持续的时间段)的比率决定。每分钟通气量由前几分钟(例如,大约两到五分钟的范围,但最好是三分钟)的平均通气量决定。
19.当前呼吸IPkFlow到当前IPkFlow 的一定比例
所述之比为最大吸气流量(IPkFlow)与当前呼吸的吸气流量中最大吸气流量(IPkFlow)的一定比例 (如75%)的比率。这一特征在本文中标记为“PF75PF”。
F.睡眠/ 觉醒监测模块
图7显示的是一个关于睡眠状态检测模块118检测方法的例子。从本质上讲,所述检测器实现了通过从反应测量透气气体流量的数据监测一种睡眠状态的方法。在770步骤时,所述方法涉及由多个呼吸流量测定多个呼吸特征。尽管呼吸特征的获取可能与睡眠状态检测仪整合,但所述获取与由特征提取模块110获得的数据相吻合。
步骤772时,所述的睡眠监测器检测方法涉及从潜在的睡眠状态来检测当时所处状态,如一个或多个非快速眼动眼睡眠状态,快速动眼睡眠状态,阶段性的快速眼动睡眠状态,补充性快速动眼睡眠状态,一种深度快速动眼睡眠状态和/或浅快速动眼睡眠状态。步骤774时,所述状态的检测基于呼吸特征的测定。
在一些实施例中,所述状态检测器主要基于呼吸特征的测定,因为所有这些状态可通过呼吸流量传感器来测定,而不依靠传统睡眠阶段传感测定技术【如脑电图数据分析(E.E.G.),肌电图(E. M. G),电图(E. 0. G.)传感器】。然后,处理器会指示检测到的状态。例如,所述处理器能够存储记忆检测到的状态,然后可以将它应用于输入另一种模块(如睡眠质量评估模块120)和/或者用于输出,例如显示在用户接口模块124。
在另一个实施例中,所述从流量信号中获得的呼吸特征可以经处理进行睡眠状态的分类,即马尔可夫分类系统。在这样的实施例中,一种非平稳马尔可夫分类系统可以用于鉴定病人的睡眠状态。
所述马尔可夫模型(Markov model)是源于允许该系统进行分析的假设之一,命名为马尔可夫属性。所述马尔可夫属性状态:假定当前状态系统,该系统的未来升级不依赖于历史状态。在每一阶段,所述系统会按照一定的概率,从当前状态变为另一个状态,或者保持相同状态。所述这种状态的改变称为转换,以及相关的各种状态改变的概率称为转换概率。
所述马尔可夫属性的数学表达式:
其中:
X:一个随机变量
N:表示一种特定状态的数字
Pr(Xn+1 = x | Xn = xn)更多地表示为Pi,j( i = Xn, j = x),这些计算可以确定状态从Xn到X的转换概率。状态转换矩阵可以包含一种状态转变为其他任何一种状态的概率(或仅仅保持当前状态)。在目前技术的一些实施例中,所述睡眠状态是一个随机变量,并且可以表现出很多可能的状态(如,以下三种情况:非快速动眼睡眠状态,觉醒状态,快速动眼睡眠状态)。因此,在接下来的这个例子,将计算三维状态的转移矩阵。由于额外的状态也会使矩阵和计算适当增加。
相应地,在一些实施例中,通过从所述特性的获取模型中输入的呼吸特点或特征,以下方法可以应用于确定睡眠状态。
(1)计算状态转换表
所述状态转换表是非平稳的(如随时间的变换而变换),因为所述转换概率的得到有直接依赖于流量信号(也随时间变化而变换)的特点。正如前面所述,一些特征是特定睡眠阶段的特点。任何数量的所述特征(如,来自特征获取模块)的组合可以用于得到这些状态的转换概率。所述这种概率可以根据每种呼吸特点设定的临界值而确定。所述临界值的设定可以通过在设备中实例分析和记录的检测到概率来确定。因此,每次睡眠状态的更新都会重新计算转换表。这就要在连续呼吸的基础上进行。在三种状态的情况下,会出现从一种状态到另一种状态以及本身状态转变的九种转换。
(2)在之前状态的基础上计算当前状态的概率(如,非快速动眼睡眠状态,觉醒状态,快速动眼睡眠状态)
[00108]由泄漏检测组件112的输出,可以为矩阵概率提供一个调节器。例如,如果存在高水平的泄漏检测,可以预计调节器的流量信号不稳定或者不指示任何特定的睡眠状态。甚至可以指示为觉醒状态。因此,在一些实施例中,所述方法可以通过适当改变睡眠状态概率的比重来减少泄漏的影响。
(3)最有可能睡眠状态的输出。
在本技术更进一步的实施例中,一种睡眠监测模块119设置了处理,如图20流程图板块中所示。所述处理代表了病人睡眠状态的一种计算方法。为此,传统的R&K PSG是基于把睡眠分为六步计算:
i . ) 清醒
ii. ) 非快速动眼步骤1 (NREM Stage 1)
iii. ) 非快速动眼步骤2 (NREM Stage 2)
iv. ) 非快速动眼步骤3( NREM Stage 3)
v. ) 非快速动眼步骤4 ( NREM Stage4)
vi . ) 快速动眼步 (REM)
所述这种睡眠阶段的计分是基于不同的生理信号包括:脑电图(EEG),眼电图(EOG)和呼吸努力频带(respiratory
effort bands)。在这个实施例中,睡眠阶段的评估主要基于呼吸流量信号的分析。由于这个原因,R&K
PSG所有阶段的睡眠结构详细的特征在所述这种模型并没有得到实现。相反,本实施例睡眠状态检测压缩为以下四种状态:
i . ) 清醒
ii. ) 非快速动眼睡眠转换(实际上是第一阶段)
iii. ) 非快速动眼睡眠(睡眠2,3,4阶段)
iv. )快速动眼睡眠
最后,所述觉醒状态和快速动眼睡眠状态被认为类似于R&K的睡眠分值。所述这些状态与呼吸流量是有关的,如图22(清醒)和图23(非快速动眼)所示。所述非快速动眼睡眠状态被认为是睡眠第3阶段和第四阶段(慢波睡眠)和部分第二阶段。图23阐明了这一阶段中呼吸流量是平稳的有节律的。图23中非快速眼动睡眠区域(以黑体NRS表示)阐明了非快速眼动睡眠状态2114与呼吸流量有关。这个区域被检测仪描述为NREM
SLEEP状态2114。
所述非快速动眼睡眠T(转换)在睡眠第一阶段是有效的,这一阶段有一个从清醒到睡眠的转变。睡眠开始时,有很典型的睡眠呼吸不稳状态(这是由于从清醒到睡眠转变时CO2设定点改变了)。这经常是与睡眠的第一阶段联系在一起的,有时候也可以延伸至第二阶段(通常可以在脑电图中看到)。在这一时期呼吸流量倾向于不规则的。图22阐示了觉醒状态的不规则的时期(arrhythmic
periods AP)的例子。这种呼吸不稳定的情况不仅限于晚上开始时的睡眠-清醒的过渡阶段。它也有可能发生在后半夜,或者如果病人醒了,准备再次进入睡眠的时候。
在以上所提到的状态的检测过程中,在本实施例中,所述检测仪将要处理一个或多个先前描述的特征提取模块110的特征,其是以来源于流量传感器的表示呼吸气流的数据为根据的。状态的决定取决于病人的信息(例如,年龄,体重指数,性别,当前生理状况)。这些信息来自于病人特征数据库122,或来自于觉醒模块144的数据和/或泄漏检测组件的数据。
在这个实施例中,一个有限状态仪器用来分类睡眠阶段。如图21所示,定义了3个基础状态:觉醒状态2110;快速动眼睡眠状态2112;非快速动眼睡眠状态2114。 在每个呼吸周期结束时,通过计算转换概率来确定一种新的睡眠状态。所述这种从一种状态到另一种状态或到同一种状态的转变,是建立在最有可能的转变机率的基础上的,如图21所示 (标记为转移概率TI1 ,T12等)。
理论上,一个病人在某种程度上及时地从任何一种状态转变到另一种状态是有概率的。然而,睡眠的生理是一种动态的数量过程,整个过程有一种随机性。因此,在每个呼吸周期结束时,应该控制和更新转换概率。在这些转换概率的基础上计算睡眠状态。进一步的,环境条件可以解释在睡眠过程中出现的不稳定(有时是周期性呼吸)。所述有限状态仪器的的另一状态,非快速眼动睡眠转变状态2116,也被定义了,如果检测到相同数量的睡眠不稳定性,状态仪将转变到非快速动眼睡眠的变换状态。所述非快速眼动睡眠转变状态2116与呼吸流量有关,如图23-A。图中,举例说明发生在睡眠开始时的睡眠不稳定区域(显示为方框的SOI)被观察仪描述为非快速眼动睡眠转变状态2116。
在这个例子中,睡眠状态检测模块119的处理方法如图20所示。所述睡眠状态检测模块是基于一个或多个特征提取模块110对呼吸特征的分析,如上所述。在测量呼吸流量的基础上,特征提取模块110可测定一个或多个特征。在这个实施例中,IPFV(吸气峰流量的改变),EPFV(呼气峰流量的改变),TSLBV(从上次确认呼吸改变以来的时间),175AV(高于75%以上的吸气峰流量的峰面积)特征可能会被分析。
所述每个特征提供到临界值单元 2010——拥有一个或多个临界值函数,如图24和25所阐述的例子,它适用于特征子集的每一种特征。然而,根据这些状态数,每一种特征被代表性地应用于每个可能的基础睡眠状态的临界值,如图25中表格所示。
例如,设定IPFV特征临界值,为每一个可能的睡眠状态,图24中,图表2450的函数可能被实施,如输入IA、IB、IC、ID和在表IPFV(图25中所示特征仪2540A)中输出OA, OB, OC 和OD,代表每一种可能的睡眠状态。正如前面这些表中所提到的,这些输入值IA、IB、IC等代表连续输出值范围。因此,对于表IPFV(图25中所示特征2540A)来说,IA是从小于0.001一直到0.001,包括0.001;IB是从大于IA的一直到0.003和包括0.003的序列值;IC是从大于IB一直到0.009,包括0.009的序列值;ID是从大于IC的范围一直到10,包括10的序列值,这是最大可能的输入值。
例如,假设当特征IPFV在确定值0.003的情况下,其输出值,可被认为是睡眠状态权重,当为0.1的特征IPFV值时输出值则被认为是醒着状态。进一步的,输出值为0.8时则被认为是由于非快速动眼睡眠状态的IPFV特征的输出值,当输出值为0.1则被认为是在快速眼动状态下的IPFV特征的输出值。同样的,输出值被选择作为EPFV特征。然而, 在 EPFV特征的情况下,每个状态下图表2450的函数的输入值和输出值都基于EPFV表中的值。同样的,归因于基本的睡眠状态的权重也被测定为表中相应的TSLBV特征和I75AV特征(显示为图25中表TSLBV的特征2540C和图25中表I75AV的特征2540D)。
虽然具体的值如图25的表中所示,但是这些数值会被理解为仅仅只是例子。其它值可采用,或凭经验从无论是一群病人或一个特定的病人中测定。因此,这些值可以通过机器学习(machine
learning)和模式识别技术(pattern recognition techniques)被学会处理。可以采用的技术实例包括聚类分析(cluster analysis),支持向量机(support vector machines)和贝叶斯分类(Bayesian
classification)。
所述输出的权重被提供给一个模式应用元件2014。在此处理过程中,各个不同的独立的特征数据相结合。为此,与当个特征相比,几个特征子集的权重有可能包含更多的睡眠阶段的指示的信息。因此,需考虑一种模式,收集一种或多种特征所显示的信息,以确定睡眠的阶段。进一步的,一系列模式的信息收集起来用来确定睡眠阶段。在这方面,这一系列模式需包含符合下列条件的信息:
i.)所述模式包含睡眠阶段的具体信息。
ii.)任何两个模式不包含有关睡眠阶段的相同信息。
满足第一个条件相对来说比较容易,第二个条件则具有更多的挑战性。假设自然条件下的呼吸流量信号和随后计算特征,创造很少或根本没有重叠的模式是很难的工作。因此,需采用一个明智的组合以实现这些模式,以获取必要的信息。
例如,在一个实施例中,模式应用仪可以适用以下模式((Pl, P2, P3, P4 和P5)。
ID模式的特征子集:
Pl IPFV, EPFV, TSLBV
P2 IPFV, EPFV
P3 TSLBV
P4 IPFV
P5 IPFV, I75AV
在所述模式的基础上,根据为每个睡眠状态设置的临界值得来的特征权重相结合,例如,通过权重总数。因此,对于每个模式(pi…Pn)和每个状态(Si…Si)来说,会产生I*N个模式得分。以上五个模式的举例,以及3个状态的例子,将有20个模式得分。( 即 , Pl_觉醒,Pl_非快速动眼睡眠, Pl_快速动眼睡眠, . . ., P5_觉醒 ,P5_非快速动眼睡眠, P5_快速动眼睡眠)。
总体模式权重计算的例子如下:
使用的特征:IPFV,EPFV,TSLBV
临界值函数的特征权重(wn):
清醒分值:
IPFV: wll
EPFV: wl2
TSLBV: wl3
非快速动眼睡眠分值:
IPFV: w21
EPFV: w22
TSLBV: w23
快速动眼睡眠分值
IPFV: w31
EPFV: w32
TSLBV: w33
所用的模式:
ID模式特征子集:
Pl IPFV, EPFV, TSLBV
P2 IPFV, EPFV
P3 TSLBV
模式权重的计算:
Pl_清醒 = α*(wll
+ wl2 + wl3)
P1_NREM sleep = β* (w21 + w22 +
w23)
Pl_快速动眼睡眠 = μ*
(w31 + w32 + w33)
P2_清醒 = αl*
(wll + w 12)
P2_非快速动眼睡眠= βl*
(w21 + w22)
P2_快速动眼睡眠= μl*
(w31 + w32)
P3_清醒 = α2*
wl3
P3_非快速动眼睡眠= β2*w23
P3快速动眼睡眠= μ2*
w33
在所述计算中,偏向因素(biasing factors)的α,β , μ , αL , βL , μL , α2 , β2和μ2如上所示可被实现。这些因素根据病人信息数据122中一个或多个病人信息允许一些模式得分的调整。例如,如果没有特定病人的信息会采取更多或更少代表一个特别的状态的特别的模式,一个或多个可设置为1。然而,其他值可能会采用,也可能根据病人的特点凭经验确定。
所述原始睡眠状态计算单元2016通过结合各自状态分值,对每一个状态生成一个特别的原始睡眠状态得分。在以上的例子中,原始睡眠状态计算单元2016将生成如下的原始清醒分值:
清醒 = Pl_清醒 +
P2_清醒+ P3_清醒;
非快速动眼_睡眠 =
Pl_非快速动眼_睡眠
+ P2_非快速动眼_睡眠
+非快速动眼_睡眠;
快速动眼_睡眠
= (Pl_快速动眼_睡眠
+ P2_快速动眼_睡眠
+ P3_快速动眼_睡眠) * (快速动眼_睡眠 _增强_指数);
在该实施例中,引入快速动眼增强因子(定义为快速动眼睡眠增强指数)。该快速动眼增强因子是由如图20中所示的快速动眼状态增强子2030生成的。快速动眼状态增强子2030在此被讨论了更多的细节。这实质上为调整(例如.通过乘法)原始快速动眼睡眠得分提供了一个额外的权重。
通常,在这个过程中,所述权重只为基于模式的得分进行计算,例如非快速动眼睡眠状态,快速动眼睡眠状态或觉醒状态。图21中的非快速动眼睡眠T状态,被认为是依赖于得分的“状态”,而不是依赖于得分的模式。在这点上,这个转换状态中的任何值都是基于周期的状态和觉醒状态时概率的计算,图27公开了关于转换概率更详细的论述。
原始状态的分值随后被状态转换概率计算单元2018处理。这种处理一般包括图表中状态转化概率确定,如图21所阐示根据此处所讨论的概率方法的转换概率识别。这些概率可考虑在睡眠阶段的计算过程中的进一步条件步骤(conditioning
step)。图26和27是概率值举例。通常,转换概率与给定状态的原始数据联系在一起,例如,用乘法得到一个修饰的睡眠状态分数。通过这种方式,采用这种转换概率的原始睡眠得分可能有偏差。然而,非快速动眼睡眠T状态,其原始状态的得分可以假定为1,这样转换概率将成为所述状态的修改后的状态得分。
在这方面,由于睡眠的生理性质,在测定睡眠状态时须考虑随机性元素。例如,第一次清醒阶段(假设睡眠前阶段)比随后整个晚上的觉醒状态要长的多(如一个更长时间段)。夜晚第一个睡眠阶段不可能是一个快速动眼睡眠状态阶段,除非病人具有特别的快速动眼状态睡眠紊乱。因此,这些规则可应用于在确定状态转换概率以解释相应的情况。图26和27举例详细说明了这些转换的条件和其相关的概率。
而且,一些转换的概率还可以取决于入睡眠起始指数,所述指数在此有更详细的讨论。通常,睡眠起始指数可能是二进制指数。例如,如果读数值为1 ,这表示病人在睡眠阶段(非快速动眼睡眠或快速动眼睡眠T)。同样的,如果读数是0,这表示病人处于觉醒状态。然而,如果病人处于晚上的第一个觉醒状态,且读数值为0, 所有管理三个睡眠状态中的任何一个状态的变化的转换概率都为零(如,如果睡眠起始指数为0,T12, T13和T14都为0)。
基于修改的睡眠状态得分,状态选择器2020将决定下一个睡眠状态。通过对修改的睡眠状态得分的比较,下一个状态将被选择。例如,修改后的睡眠状态的最大分值可被视为图21中状态机的状态检测所检测到的状态。
所讨论的睡眠状态检测的检测方法可按如下例子考虑:
例子:
有限状态机(Finite Machine) 的当前睡眠状态= 觉醒状态 2110
原始睡眠状态得分计算:
清醒 = 0.1
非快速动眼睡眠 = 0.5
快速动眼睡眠 = 0.2
非快速动眼睡眠转换= 1
基于当前的睡眠状态的转移概率是TIL, T12,T13, T14。TIL ,T13和T14 ,参考表2600A ,假设规则6适用于这个例子。T12可参考表2700A T,假定规则2适用。
TIl =
0.45
T12 = 0
T13 = 0.45
T14 = 0.1
修饰的睡眠状态得分:
清醒= 清醒*Tll
= 0.045
非快速动眼睡眠_T = T12 = 0
非快速动眼睡眠=非快速动眼睡眠*T13
= 0.225
快速动眼睡眠=非快速动眼睡眠*T14
= 0.02
由于非快速动眼睡眠状态是得分最高的,下一步的睡眠状态,可以设置为非快速动眼睡眠状态2114。
G.快速动眼睡眠状态增强仪。
当只分析呼吸气流信号,快速动眼睡眠状态是最难检测的睡眠状态。因此,一些特别的特征可以被特异地应用到捕捉“快速眼动睡眠型”气流的特点,以加强快速动眼睡眠的检测。这些特征的结合输出产生了一个增强的快速动眼状态,可部分地由如图20所示的快速动眼状态增强单元2030实施。
例如,以下附加特征可以被确定,选择性地,至少部分由处理特征提取模块110实施。为此,呼吸气流被捕捉、处理以计算随后的快速动眼所特有的特征。
1)VtRat -当前呼吸吸气潮气量(VTcurr)与之前呼吸吸气潮气量(VTprev)之间的简单比。在一些实施例中,这个比例可根据下列公式计算。
2)正弦曲线的偏差 - 这个特性决定了取决于当前的吸气的幅度和频率的正弦波谱。估计的正弦曲线和实际的吸气量变曲线之间相减,结果被存储在一个缓冲区。然后计算缓冲区中数值的方差。在病人的呼吸文档的基础上,正弦曲线可被修改。例如,这些值可以从仪器学习和模式识别技术中掌握。可以采用的技术实例包括聚类分析,支持向量机和贝叶斯分类。
在一个应用实例中,一个具有可变幅度和频率的简单的平方根的正弦曲线也可使用。因此,基于此配置文件,本实施例计算可精确地表示如下:
其中,α是振幅,ω是频率,t是时间,Φ是相位。
3)特征一致性指数 - 该指数基于一套临界值测定特征的相似度。该指数的计算涉及以下步骤:
i. 收集一些呼吸信息,对每一个呼吸进行以上所有提到的特征的计算,在临界值缓冲区中存储。例如,收集5次呼吸并确定相关的呼吸特征,在临界值缓冲区中存储。
ii. 所述特征应用于临界值函数。临界值函数类似于前面描述的。每次一个特征值超过设定的特征临界值,一致性计算器递增一次(如递增1)。例如,在5次呼吸中,如果一个特征高于呼吸临界值4,一致性计数器将显示4。
iii. 然后根据特征一致性计数器数值和收集的总的呼吸样本数之间的比率推导出特征一致性指数。例如,如果计数器为4,收集呼吸总数为5 ,所述特征一致性指数为80%。
特征一致性指数可由快速动眼睡眠状态增强仪2030输出作为前面讨论的快速动眼睡眠状态增强指数。
H.转换概率规则:
如前所述,为产生修正的睡眠得分,转换概率控制原始睡眠得分。然而,定义如何计算这些概率之前,须考虑一些自然睡眠周期通常的特点:
(a)夜晚(即睡眠开始)第一次清醒阶段可能超过最长的晚上清醒时间。
(b)随后清醒时间不会持续超过120秒,除非有外部原因导致觉醒(如,面罩调整,进行一次如厕中断,被人唤醒或者任何医学措施,例如夜间哮喘发作)。
(c)夜晚第一个睡眠阶段不可能是一个快速动眼睡眠状态阶段,除非病人具有特别的快速动眼状态睡眠紊乱。
(d)作为自然睡眠周期的一部分,病人在深度睡眠中超过20至30分钟,越来越接近一个小的觉醒。
综上所述,以下所述规则已经定义。不过,也有许多方法来解释上述信息( (a)-(d) ),并定义相应的规则。 因此,下面的规则和与其相关的转换概率值仅仅是举例。
其他的规则和价值观念可从仪器学习和模式识别技术中掌握。可以采用的技术实例包括聚类分析,支持向量机和贝叶斯分类。
规则:
规则1:病人处于晚上的第一个清醒时期,并且已经清醒少于或等于10分钟。
规则2:病人处于晚上的第一个清醒时期,并且已经清醒了超过10分钟,但是不到20分钟。
规则3:病人是在晚上的第一个清醒时期,并且已经清醒超过20分钟。
规则4:病人处于晚上的第二个清醒时期,并且已经清醒了少于或等于10分钟。
规则5:病人处于晚上的第二个清醒时期,并且已经清醒了超过10分钟,但是还不到20分钟。
规则6:病人处于晚上的第二个清醒时期,并且已经清醒超过20分钟。
规则7:病人处于晚上的第三个或以后的清醒时期,并且已经处于这个状态少于或等于60秒。
规则8:病人处于晚上的第三个或以后的清醒时期,并且已经处于这个状态超过60秒。
规则9:从治疗开始时的时间少于或等于20分钟。
规则10:从治疗开始时的时间已经超过20分钟,且病人已经处于非快速动眼睡眠状态少于或等于20分钟。
规则11:从治疗开始时的时间已经远远超过20分钟,且病人已经处于非快速动眼睡眠状态超过20分钟但少于或等于30分钟。
规则12:从治疗开始时的时间已经远远多于20分钟,且病人已经处于非快速眼动睡眠状态超过30分钟。
规则13:病人处于快速眼动睡眠状态已经少于或等于10分钟。
规则14:病人处于快速眼动睡眠状态已经多于10分钟,且少于或等于20min。
规则15:病人处于快速眼动睡眠状态已经远远多于20分钟。
图表26总结了所述规则与概率值的相关性,所述概率值可用来管理状态转换概率。在此实施例中,在任何一个时间仅只有一个特定的规则可适用。在规则满足当前条件的基础上(如,当前睡眠状态,目前睡眠状态所花的时间,在晚上睡眠时间内的全部时刻),如前所述,表格的输出可作为转换概率来使用。例如,如果处理器确定第五条规则是正确的,且机器显示当前状态是清醒的,那么转移概率TIl则为0.9 , T13为0.05,T14为0.05。表中其余的转移概率同样可以确定。
I.周期性呼吸﹠基于转变概率的觉醒
如图20所示,另一组利用状态转换概率计算器2018测定状态转变概率的规则是以周期性呼吸探测器和觉醒探测器的数据输出为基础。例如,周期状态检测仪控制从和到非快速动眼睡眠T状态的转换。为此,在图表27详述的以下规则和相关的数值同样可掌控状态转换概率。
规则1:病人当前处于觉醒状态
周期状态计数器数值是0(病人处于周期状态中达零次呼吸)。
规则2:病人当前处于觉醒状态
周期状态计数器数值比0大,但是小于或等于50.(病人处于周期状态中有超过0次,但不多于50次呼吸)。
规则3:病人当前处于觉醒状态
周期状态计数器值大于50(病人处于周期状态中有超过50次呼吸)。
规则4:病人当前处于非快速动眼睡眠状态
周期状态计数器(Periodic State Counter)为0(病人处于周期状态中达零次呼吸)。
规则5:病人处于非快速动眼睡眠状态
周期状态计数器数值比0大,但是小于或等于50.(病人处于大于0而小于等于50次呼吸的周期性状态)。
觉醒状态缓冲总数( Arousal State Buffer Sum) =0(所述缓冲包括源于先前50次呼吸的、先前50次觉醒状态值的数值,所以如果sum=0,则前50次呼吸中没有检测到清醒)。
规则6:病人处于非快速眼动睡眠状态
周期状态计数器是比0大,但是小于或等于50 (病人处于大于0而小于等于50次呼吸的周期性状态)
觉醒状态缓冲总数>0(所述缓冲包括源于先前50次呼吸的、先前50次觉醒状态值的数值,所以如果觉醒状态总数大于0,则前50次呼吸中,一次或更多次清醒 状态被检测到)。
规则7:患者当前状态是非快速动眼睡眠
周期状态计数器数值大于50(患者处于周期状态中有超过50次呼吸)。
觉醒状态缓冲总数 =0(所述缓冲包括源于先前50次呼吸的、先前50次觉醒状态值的数值,所以如果sum=0,则前50次呼吸中没有检测到清醒)。
规则8:患者当前状态是非快速动眼睡眠
周期状态计数器数值大于50(患者处于周期状态中有超过50次呼吸)。
觉醒状态缓冲总数 >0(所述缓冲总数包括源于先前50次呼吸的、先前50次觉醒状态值的数值,所以如果sum> 0,则先前50次呼吸中检测到一次或更多次清醒)。
规则9:患者当前状态是非快速动眼睡眠T
周期状态计数值为0(患者处于周期状态中达零次呼吸)。
规则10:患者当前状态是非快速动眼睡眠转换
周期状态计数器数值比0大,但是小于或等于50(患者处于周期状态中有超过0次,但不多于50次呼吸)。
觉醒状态缓冲总数=0(所述缓冲总数包括源于先前50次呼吸的、先前50次觉醒状态值的数值,所以如果sum=0,则先前50次呼吸中没有检测到清醒)。
规则11:患者当前状态是非快速动眼睡眠转换
周期状态计数器数值比0大,但是小于或等于50(患者处于周期状态中有超过0次,但不多于50次呼吸)。
觉醒状态缓冲总数>0(所述缓冲总数包括源于先前50次呼吸的、先前50次觉醒状态值的数值,所以如果sum> 0,则先前50次呼吸中检测到一次或更多次清醒)。
规则12:患者当前状态是非快速动眼睡眠转换
周期状态计数器数值大于50(患者处于周期状态中有超过50次呼吸)。
觉醒状态缓冲总数=0(所述缓冲总数包括源于先前50次呼吸的、先前50次觉醒状态值的数值,所以如果sum=0,则先前50次呼吸中没有检测到清醒)。
规则13:患者当前状态是非快速动眼睡眠转换
周期状态计数器数值远远大于50(患者处于周期状态中有超过50次呼吸)。
觉醒状态缓冲总数>0(所述缓冲总数包括源于先前50次呼吸的、先前50次觉醒状态值的数值,所以如果sum> 0,则前50次呼吸中检测到一次或更多次清醒)。
规则14:患者当前状态是快速动眼睡眠
周期状态计数器数值是0(患者处于周期状态中达零次呼吸)。
规则15:患者当前状态是快速眼动睡眠
周期状态计数器数值比0大,但是小于或等于100(患者处于周期状态中有超过0次,但不多于100次呼吸)。
规则16:患者当前状态是快速眼动睡眠
周期状态计数是远远大于100(患者处于超过100次呼吸的周期状态中)。
图表27(表2700A,表2700B,表2700C和表 2700D)总结了所述具可能性的这些规则的联系,其可应用于掌控状态转换概率上。基于对一个或多个规则的评估,以及在上一次呼吸中机器所检测到的特殊状态,输出的表格可应用于如前所述的转换概率。例如,如果处理确定规则11是正确的,机器的当前状态是非快速眼动睡眠转变状态2116,那么转换概率T22则为0。同样可确定表中的其它的转换概率。
J.觉醒检测模块
觉醒检测模块114决定呼吸流量数据是否代表从睡眠中醒来。在这点上,从睡眠中醒来可以导致流量信号的干扰。因此,本模块的基本功能是扫描流量数据中的干扰信号,以表明他们是否已或者还没有发生。例如,如果特征提取模块110从呼吸特征中检测到干扰,模块可能会产生一个载有干扰流量的信号或数据传给睡眠稳定性模块。选择性地,包含流量信号干扰的数据包可与干扰之前的周期流量信号产生的数据同时产生。这些数据可能会进行进一步的干扰分类,在此已有详细的描述。
除了如前面所讨论的那样选择性地从特征提取模块110中利用呼吸特征,其他呼吸特征也可以计算。例如,可以测定以下的呼吸特征:
(1)
吸气峰流量到达一定比例(如95%)的时间,
(2)
吸气峰流量高于一定比例(如75%)的面积,
(3)
吸气峰流量高于一定比例(如75%)的面积的变化,
(4)
到达一定比例(如95%)的吸气峰流量和呼气流量之间的时间。
一个检测干扰的实例,使用一个或更多的呼吸特征检测。例如,通过评估上述列举的1-4项呼吸特征检测流量的干扰。这些特征可被认为是流量干扰的指标。可选择的,这些特征,和或者不和之前提到的特征一起,可以比较作为一套经验的临界值来检测干扰。例如,如果比较结果显示发生在流量信号中的干扰信号,模块则会显示,例如,输出一个积极的信号或数据显示(如显示a 1表示有,a 0表示没有)。
因此,一个模块在循环呼吸的基础上的处理方法通过举例如下:
(1)从特征提取仪中输入特征。
(2)计算以上列举的附加特征(1-4)。
(3)从睡眠状态检测模块中输入当前睡眠状态。
(4)确定流量干扰是否存在。
(5)如果存在一个流量干扰,需检查睡眠状态。如果患者在睡眠状态,从所述干扰和先前的呼吸(如15次呼吸)搜集流量数据并输出数据(如输入到睡眠质量模块)。如果存在一个流量干扰,且患者不是处于睡眠状态,不需要采取任何措施。如果呼吸中没有流量干扰,不需要采取任何措施。
在一些实施例中,对代表每个流量干扰的数据进行分析,以确定觉醒存在已否。如果检测到觉醒,须进一步分析,以确定觉醒为以下几种类型之一:
(i)窒息相关的觉醒(Apnea
related arousal)
(ii)呼吸不全相关的觉醒(Hypopnea related arousal)
(iii)呼吸意外相关的觉醒(Respiratory
effort related arousal)
(iv)面罩脱离导致的非呼吸(Non-Respiratory)
(v)周期性腿动(PLM)引起的非呼吸
(vi)高泄漏(High
Leak)引起的非呼吸
(vii)非呼吸-自发性觉醒
例如,有了来自于特征提取模块110中的以上所讨论的呼吸特征的数据,所述这些特征可以应用于临界值函数,这就关系到每一个特征的权重。流量干扰特征源自临界值函数输出的权重。进一步分析流量信号,以确定上面提到的流量干扰之前所发生的事件。随后,所述流量干扰被确定为一种觉醒,模块的输出可以选择性的包括
(a)觉醒模型和(b)觉醒持续模型。
觉醒检测器实施例
进一步的,在本项技术的实施例中,觉醒检测模块114可配套有如图13中的流量图表的方框内所示的处理。呼吸气流或过滤的气流与流量传感器或流量发生器相连,传送给特征提取模块并进行处理。在这个处理过程中,前面所述的呼吸特征的特殊子集被特征提取模块测定并利用。
选择性地,特征子集可包括IP75A(例如,高于75%吸气峰流量的吸气流量曲线的面积)、3mvTtot(当前 3 分钟通气量与Ttot (总呼吸时间)之间的比率)以及PF75PF (例如,吸气峰流量和 75%吸气峰流之间的比率)。每一种特征都提供给临界值单元1310,其含有一个多个临界值函数,如图14所示,含有所有特征的特征子集。虽然这个具体实例中应用了其中一些特定特征,但在可以认为,在其他实例中可应用通过特征获取仪通过测定呼吸流量测定的前面所述特征的任何组合的一套(如全部)或者子集。
在这个实施例中,当用IP75A特征的临界值时,图. 14中图表1450 的函数,输入数据 如IA 和 IB经,输出数据如OA和 OB,显示在IP75A表格中(表中的字符表示如:1440A),在表中的IA,
IB, IC的输入值代表了输出值的连续性范围。因此,表IP75A中,IA的值范围从小于2.5并包括2.5,IB的值比IA大或者比2.5大。比如,IP75A特征值为
2.5,考虑到干扰的权重,IP75A的特征输出值选择在0.05,也就是说输出权重选择为3mvTtot 特征。
然而,在3mvTtot 特征的情况下 ,图 1450函数的输入输出值将基于 3mvTtot 表的值(表中的字符表示如:1440B)。同样的,根据1450函数及表中的值将输出权重作为PF75PF特征(如1440C所示)。需要强调的是, 图14表中特定的值,在这里只是举例说明。其它值可采用,或凭经验从无论是一群患者或一个特定的患者中测定。因此,这些值可以从仪器学习和模式识别技术中掌握。可以采用的技术实例包括聚类分析,支持向量机和贝叶斯分类。
结合三个临界值的权重,比如三个权重的总和,在这个处理过程中,干扰计算器1312会产生一个流量干扰特征。干扰计算器1312比较干扰特征与觉醒临界值的权重,该权重的选择可能是经验性的,以指示干扰特征是否表示觉醒。如果是,设作觉醒标志。
基于觉醒标志的阳性指示,该觉醒事件的数据流经事件分类仪1314分析,以识别可能会导致或引起觉醒的特定事件。例如,通过识别的事件可能是呼吸相关的(如,障碍性或中枢性窒息,障碍性或中枢性呼吸不足或呼吸受限等),或不相关的(如,周期性腿动(PLM)和/或泄漏事件(如口泄漏或高泄漏)。
如果不是这些事件,那么觉醒被认为自发性觉醒。在一些实施例中,识别觉醒事件与事件分类仪1314的分类过程结合。然而,如 图 13所示,它们也可不连续的测定然后将结果传递给事件分类仪1314。最后,觉醒探测器114输出觉醒特征,可以使一个信号或者数据显示觉醒的类别和/或者持续事件。例如,可以是识别为口泄漏相关的觉醒事件,窒息引起的觉醒,流量限制引起的觉醒,周期性腿动事件引起的觉醒,或自发性觉醒等。觉醒持续的时间段是从设定觉醒标志到确定没有觉醒发生。交替地,时间段可以是一些呼吸。
K. 睡眠稳定性检测模块
在本项技术的一些实施例中,睡眠稳定性检测模块116执行多种功能,比如,与觉醒检测模块结合,该模块可以(1)检测干扰流量,(2)分析干扰强度(如觉醒程度等),和/或(3)推断自主神经激活水平。
一种结合觉醒检测模块和睡眠稳定性检测模块的方法实例,如图8所示。本质上,所述模块配置有可以根据呼吸流量的数据检测或分类睡眠觉醒状况的方法。步骤880,所述方法包括通过测定呼吸气体的流量以确定一组呼吸特性。这可能与检测模块结合,或由特征获取模块提供。在882步骤,前述的呼吸特征检测到了干扰。所述干扰可表明觉醒状态。
选择性地在884进行详细说明,所述方法可以评估干扰是否为非呼吸 (non-respiratory)相关的觉醒事件。步骤886,所述模块接着指示检测到的干扰代表觉醒事件而不是因呼吸功能障碍的症状。
例如,它可能储存一些记忆中的测定数据,将其输入另一模块(如睡眠质量评估模块120)和/或作为操作者界面模块124的输出数据。例如,这种非呼吸功能障碍相关事件可能与周期性腿动相关。所以,这种非呼吸相关觉醒可能是周期性腿动或检测到的泄漏引起的。换句话说,周期性腿动出现次数或持续时间可以从检测觉醒或干扰,以及检测呼吸功能障碍(respiratory
dysfunction)的同步消失(例如,没有或无关紧要的流量抑制,没有或无关紧要的流量平缓,没有或无关紧要的呼吸梗阻等)。
如图 9 所示,进一步的实施例,所述方法与觉醒检测模块114和睡眠稳定性检测模块116相结合。本质上,所述模块通过一种方法,可以从可测量呼吸流分析睡眠稳定性。在990步骤,所述方法,包括通过测定呼吸气体的流量以确定一组呼吸特性。这可能与检测模块结合,或由特征提取模块110提供.在992步骤,前述的呼吸特征检测到了干扰。所述干扰可表明觉醒状态。在994步骤,检测干扰程度。干扰程度可指示觉醒状态对睡眠的干扰程度。接着处理器就会显示干扰程度。例如,它可能储存一些记忆中的测定数据,将其输入其他模块(如睡眠质量模块120)和/或作为操作者界面模块124的输出数据。
在下面的实例中,干扰流量已经过分类、分级,并将其作为患者自主神经激活的指标。在该实施例中,以下步骤可被执行:
(1)从觉醒检测器接收干扰流数据包
选择性地,接收干扰流数据的接收模块本身不处理流量信号。
(2)分析在干扰流数据包之前的数据,将干扰进行分类
例如,如果在检测到干扰流之前检测到如窒息、呼吸不足等,该干扰数据就被当做是呼吸相关的觉醒。如果在检测到干扰流之前没有检测到这些事件,那么就被分类为周期性腿动,而不是呼吸功能障碍。
选择性地,如果泄漏检测组件同时检测到泄漏,该仪器就会把所述觉醒或干扰归因于面罩的瓣式泄漏,而不是非呼吸相关事件。
(3)
计算流量干扰的长短(如,连续时间)和强度。
例如,干扰时间的长短可以是测定或计算干扰持续或所需的时间(如,秒或者分)。在一些实施例中,除了确定干扰时间的长短,还要确定强度(如,干扰程度)。例如,强度值可以通过计算干扰前呼吸数据和干扰数据的比例得到。例如,强度值也可以计算(a)数据包(如,干扰期间信号流量样本)中干扰数据的变化和(b)干扰前(干扰前信号流样品)一定数量的呼吸(如15次)数据变化的比率。任何呼吸特征都可用来计算这个比率。
随意地,所述测定的强度值可用于指示患者在干扰期间的自主神经激活。例如所述模块包含或含有一个内存查询表。所述查询表的特征是,包含干扰流强度映射自主神经激活程度的数据(如,心率变动值,脉搏传导值(PTT),等)。因此,强度值可用来指示自主神经激发值。选择性地,所述查询表可能被患者信息存储模块122处理的数据索引。
这种情况下,所述查询表就包括不同类群患者的数据。因此,根据患者的特征,如:年龄、身体质量指数、身高,和/ 或其他一些病症及现在的生理状况等,自主神经激活的测定可能进一步被索引,提供给使用该设备的不同的患者。这样的查询表可通过挖掘大规模数据来扩展,用于比较和关联干扰流强度及之后的自主变化(如,不同类患者的心率变异度(HRV)和脉搏传导时间(PTT))。
(4)显示睡眠稳定性得分给睡眠质量索引模块
干扰强度值、长短、自主激活值及干扰类型等,被检测后被选择性的记录在仪器的存储器上,提供给其他评估模块(如,睡眠质量评估模块120)和/或显示在用户界面模块124。
总之,呼吸治疗(如气管正压)的目的是减少患者呼吸相关的干扰,以此增强睡眠稳定性。因此,睡眠检测模块实例可将患者睡觉期间与睡眠相关的干扰量化。例如,具有这个组件的检测器在每个呼吸周期结束时产生0到1的数值,以指示睡眠干扰水平。此时,1分表明睡眠最稳定,当分数降向0时睡眠稳定性逐渐降低。
最后,在一些实施例中,将特征获取仪确定的特征应用于一个或多个临界值函数,该权重应用于每个特征。临界值函数可能根据患者信息存储模块122提供的患者特征进行适当调整。通过临界值函数输出的数据可计算呼吸干扰值。所述干扰值还要减去泄漏检测仪112和觉醒检测模块114的输出值。通过多次呼吸测定的平均值就可以计算出睡眠稳定分数(如,5-30次范围内的15次呼吸的平均值)。
实施例:睡眠稳定性检测装置
在这个实例中,睡眠检测模块116配置了如图15所示的处理模块。呼吸气流或过滤的气流与流量传感器或流量发生器相连,传送给特征提取装置并进行处理。在这个处理过程中,前面所述的呼吸特征的特殊子集被特征提取装置测定并利用。
选择性地,特征子集可包括如前所述的IPFV(吸气峰流量的变化)功能,EPFV(呼吸峰流量变化)功能,TSLBV(最后的呼吸时间变化,) I75AV(75%以上吸气峰值流量变化)。每一种特征都提供给临界值单元1410,其含有一个或多个临界值函数,如图16所示,含有所有特征的特征子集。虽然本实施例中应用了其中一些特定特征,但在可以认为,在其他实例中可应用通过特征获取仪通过测定呼吸流量测定的前面所述特征的任何组合的一套(如全部特征)或者子集。
在这个实例中,当用IPFV临界值时,图. 16中图表1650 的函数,输入数据 如IA, IB IC 和ID,显示在IP75A表格中则为输出数据如OA,OB,OC 和 OD(可用1640A表示)。在表中的IA, IB, IC的输入值代表了输出值的连续性范围。因此,表IPFV中,IA表示小于等于0.001的范围,IB表示大于IA且大于0.003的范围,IC表示大于IB且大于0.009的范围,ID表示大于IC且大于10的范围,这也是最大输入值。例如,IPFV特征的确定值为0.003,考虑到干扰的权重,IPFV特征输出值为0.8。
同样地,输出权重作为EPFV特征。但是,考虑到IPFV特征中,图 1650函数的输入值和输出值将基于 EPFV表的值(表中的字符表示如:1640B),通过1650图的函数并基于TSLBV表中的值得到TSLBV特征的输出值(如1640C所示)。类似的,根据1650图的函数并基于 I75AV表中的值得到I75AV特征的输出值(如1640D所示)。需要强调的是, 图14表中特定的值,在这里只是举例说明。其它值也可采用,或凭经验从无论是一群患者或一个特定的患者中测定。因此,这些值可以从仪器学习和模式识别技术中掌握。可以采用的技术实例包括聚类分析,支持向量机和贝叶斯分类
通过结合权重确定原始干扰值,如在计分组件1512得到初始数据经处理后临界值仪器1520进行权重相加。计分组件1512原始数据再根据系统情况进行调整或de-weighted。例如,根据泄漏检测,或如图15所示泄漏检测模块和/或觉醒模块的觉醒情况调整计分组件1512的数据。
例如, 通过减重函数对减重因子测定后的减重初始干扰分值 ,来自于当前呼吸的修正后的干扰分值1522可被确定,如图17函数图例所示。因此,图17所示函数可用于考虑目前的泄漏水平。例如,如果泄漏检测装置检测到系统中出现高泄漏(如,泄漏超过0.4 1/s),实际呼吸流信号的信噪比(SNR)就会明显降低。在这种情况下,初始干扰分值将用于分析这种低信噪比。当泄漏达到0.3
1/s,图17临界值利用输出减重因子(de-weighting
factor)逐渐降低睡眠稳定性。一旦泄漏达到 0.5 1/s,睡眠稳定值将降为0。
同样地,这种减重元件(de-weighting component)可用于考虑觉醒检测装置的觉醒状态。例如,呼吸、泄漏及PLM相关的觉醒造成患者睡眠中断,这种类型的觉醒涉及到睡眠稳定性分值。然而,自发性觉醒应该属于患者本身睡眠中的一部分。在这种情况下,对于一些特定的呼吸如果把自发性觉醒考虑进去,此时的睡眠稳定分值是错误的。因此,减重装置(de-weighter) 1514可以分析输入的觉醒特征,如果检测到自发性觉醒,可以阻止修正后的干扰分值被算入当前呼吸的睡眠稳定分值。
然后,过滤器1516,如旋转平均过滤器,处理来自减重 装置1514的干扰数据。例如,对减量器输出的一系列分值求平均值,例如10到20个分值,最好是15个分值。最后,滤器输出值就被当做睡眠稳定性分值1520,该值可代表睡眠稳定性指标。
L.睡眠质量评估模块
一种睡眠质量评估模块120可以用来实施评估睡眠质量。在一些实施例中,所述的评估或检测是基于其他模块的输出。图10详细解释了一种与睡眠质量评估模块120相关实例。从本质上讲,所述方法从可吸入气体流量的测量来指示睡眠质量指标。
在1000位置时,类似于前述模块的输入,通过测定呼吸气体的流量以确定多个呼吸特征。这可能是集成部分模块或其他部分模块(如特征提取模块110)。在1002和1004位置时,所述方法涉及根据所述多个呼吸特征分别确定睡眠的状态和睡眠的稳定性。在1006和1008位置时,通过睡眠状态和稳定性检测,所述方法设计分别确定和指示一种或多种睡眠质量指标/指数。因此,所述这些指标或指数可以储存在设备中,被其他评估模块(如反馈模型126)使用,和/或显示在用户界面模块124。
例如,一个实施例可通过以下步骤实现:
(1)计算睡眠时长(如,快速动眼,非快速动眼,深度快速动眼,轻度快速动眼,阶段性快速动眼,动力快速动眼,所有睡眠状态,等),作为一个索引。
(2)计算睡眠稳定性水平和总睡眠时间的比例作为一个索引,例如,所述一个睡眠稳定性指数除以一个总的睡眠时间。这样有助于稳定指数标准化。
(3)所述比率与夜间清醒时期的长度加权作为一个索引。这个步骤可以考虑到觉醒的频率。例如,全部清醒时间为数字X。但是,如果所述这个数字是由小部分清醒时期Y组成,那么数字Y就很大,即可通过加权因子而计算。换句话说,所述睡眠稳定值通过计算夜间惊醒的频率会更加趋于标准化。
(4)把因单个患者特点而产生的偏差的比率加权作为一个索引;在一些情况下,一种觉醒状态上升时期会有带有一定水平的偏差。基于所述偏差属于特定患者的特点,加权因子有利于比率更进一步趋于标准化。
(5)输出最后睡眠质量得分作为一系列指标。
选择性地,出于反馈的目的,所述模块作为一个索引也可以监视或记录从这一夜到下一夜睡眠时期觉醒的多个故障。所述这样的一个数字或多个数字可以是个简单的数,它可以计算清醒期的数量,总的清醒时间,睡眠期的数量和/或总的睡眠时间。选择性地,这可以是一个时间段的平均数,如一个月。日复一日的数字也可以。所述这样的信息为患者的睡眠结构提供简单的反馈。
所述模型的输出可被记录,用于反馈和/或显示给患者和/或临床医生。例如,它可以通过用户界面模块124得到显示,和/或直接返回或传入治疗算法处理模块128,以便在恰当时间适当控制流量发生器。所述这种反馈有助于实现治疗短期的变化,和/或长期的改变(如几周或几个月)。
睡眠质量检测检测
进一步的,睡眠质量评估模块2010的实施例可参考图28-31。通常,睡眠质量检测装置的处理过程涉及系统的其它的处理数据如睡眠稳定性指数的评估,睡眠状态标识和/或觉醒指数。输入数据可以作为一个动态的力量以在循环呼吸的基础上改变睡眠质量指数(SQI)。评估模块的处理过程如图28所示。在这个实例中,输入数据可应用于临界值单元2810A,
2810B 和2810C
在评估模块处理睡眠状态标识符时,临界值仪2810A可运行一个如图29所示的函数。在这个例子中,临界值仪2810A的输出值可能是在睡眠状态标识符的基础上从-0.05到0.5之间变化的一个值。通常,图29的函数显示,非快速眼动睡眠具有最高的积极影响,而觉醒却有最高的消极影响。如此,如果输入状态标识符代表觉醒状态2010,那么输出值就是-0.5.如果状态标识符代表非快速眼动睡眠转变状态2116,那么输出值是-0.3. 如果状态标识符代表非快速眼动睡眠T状态2114,那么输出值是0.5。这样,状态检测仪通过临界值单元可具有额外的检测状态,额外的加权值。
进一步实施例中,快速动眼状态的检测可以导致对睡眠质量指数的调节。例如,如果检测到快速动眼状态,就可以得出睡眠质量指数。因此,在一些实施例中,睡眠质量指数将会增加作为一个总时间长度的函数,包括一个或更多的快速动眼状态。选择性地,它在检测到的快速动眼时期的数量增加时可以作为一个增加的函数。更进一步,作为一个在晚上快速眼动睡眠状态发生时,睡眠的特定时间点的函数可能增加。在进一步的实施例中,快速眼动状态的检测可能不会导致对睡眠质量指数的调节。
评估模块处理睡眠稳定指数或睡眠稳定分值时,临界值单元2018B可应用如图30所示的函数。在这个实施例中,临界值单元2018B的输出值也可能是在睡眠稳定指数的基础上从-0.05到0.5之间的一个值,输入睡眠稳定指数与输出睡眠稳定指数之间的联系显示在函数线3010.例如,大约0.3的睡眠稳定指数的输入值导致大约0.1的输出值。
所述评估模块处理觉醒数据,如觉醒类型时,临界值单元2018C可应用如图30所示的函数。在这个实施例中,临界值单元2018C的输出值可基于睡眠稳定指数,是从-0.05到0.5之间一个值。这样,如果输入觉醒类型代表呼吸暂停,那么输出值是-0.5。如果输入觉醒类型代表呼吸不全,那么输出值是-0.4.如果输入觉醒类型代表RERA(有关觉醒的呼吸作用),那么输出值是-0.1.如果输入觉醒类型代表泄漏,例如口泄漏或高泄漏,那么输出值是-0.5.最后,如果输入觉醒类型代表自觉的觉醒,那么输出值是0。
基于临界值单元的输出,所述睡眠质量调节装置中设定一个睡眠质量指数,这是介于0和1之间的实数值。例如,通过应用睡眠质量指数作为临界值权重,现有的睡眠质量可能会提高或降低。例如,睡眠质量指数通过大量重复的一个消极权重因子会降低或因一个积极的权重因子而增加。可选的,根据输入的数据,假如临界值函数是正的或负的,则所述临界值仅仅输出对应的正确或错误的值。所述正确或错误的值分别控制对应的递增或递减睡眠质量指数,例如一个复合值。所述这种评估大体上基于以下前提:(a)某些觉醒通常会降低睡眠质量而没有觉醒则会提高睡眠质量。(b)睡眠稳定性差会降低睡眠质量,(c)一些睡眠状态会显示更好的睡眠质量。
如图29-31展示了上述函数,已选择的权重许可睡眠质量指数的输出,输出值在0和1之间,这是睡眠质量的一个恰当指标,其中,1表示最好的睡眠质量,0表示最差的睡眠质量。然而,其他权重值和指标值的选择则要依赖于想得到的睡眠质量指数的性质。
M.反馈控制模块
因此,在一些实施例中,系统所述这些模块的输出可用来控制治疗控制运算的变化,例如,一种由处理器控制流量发生器提供的压力治疗法【如,持续正压通气治疗(“CPAP”)】。例如,睡眠质量评估的输出或睡眠状态的输出会被导入反馈控制模块126.。然后,反馈控制模块就会实施控制,如,给予患者瞬间和长期的大量压力疗法。所述这种反馈控制模块,通过发送指令或控制流量信号给流量发生器或治疗控制处理模块128,而实施动态的或短期的和/或静态的或者长期的变化以完成治疗。
在一个实施例中,一种动态反馈控制是基于不同睡眠状态的输出。例如,如果患者起初检测到觉醒状态,提供的疗法就保持在最低水平。类似的,呼吸检测历程(如AHI模块)可设置禁用,否者对任何流量干扰设置更低敏感度的权重,或者避免或减少压力调整以应对AHI事件。当患者过渡到检测到睡眠状态,可将所述反馈控制模型的控制器设置为对干扰或AHI事件更为敏感。在深度睡眠检测时期(如,深度快速动眼睡眠状态或快速动眼睡眠状态),所述压力疗法算法灵敏度水平可以设置为最高,以充分应对或响应检测到的事件。仍可选的,假如所述设备检测到一位患者醒了,处于一种检测到的觉醒和/或清醒状态,随后,这种疗法就会以逐步稳定的速率下滑或减少,以允许患者更舒适地返回到睡眠。
在一些实施例中,所述这种反馈模块对于治疗控制会表现出静态改变。例如,一个睡眠质量指数可以用于长期反馈。由于患者睡眠模式的改变是在时间上通过改变睡眠质量指数所测得,因此治疗运算的敏感水平也会做相应的调整。例如,改善睡眠质量指数可以控制降低治水平。同样,睡眠质量指数的降低,还能控制治疗水平的提高。这种长期的改变会选择性地表现控制时段,如,仅仅在一定数目的治疗疗程之后,在每个月的最初,或者治疗达到一个月,等。
在更进一步的实施例中,额外的觉醒导致睡眠的中断的打分也可控制治疗的改变。例如,假如在一段较长的时期内(如,一个月)清醒时期很长或比例不断增加,这样,压力控制运算根据分值,通过加权因子提供更大程度的治疗。
基于上述模型的输出,以下额外反馈实施例可实现。
A .清醒呼吸压释放(EPR)
在一些呼吸压力治疗装置中,所述压力控制包括了一种呼气压力释放(EPR)的特征。这样的装置在国际专利申请编号为PCT/US2004019598 (Publ. No. WO 2004/112680)和相应的美国专利编号为7,128,069中已经有了描述,通过参考这些专利所公布的信息与本文相结合。所述EPR特征可为患者提供舒适,以便让患者更容易获得睡眠。
一般来说,根据检测到的每个患者的呼气,所述EPR可自动控制降低输送的治疗压力通过压力疗法的某个数量(如cm H2O),当吸气期间输送治疗压力时。所述压力治疗控制可根据目前的技术进行修改。
通常情况下,一旦患者睡着了,EPR对于患者的治疗被认为是没有必要的,因为EPR不能提供有利于患者最理想的治疗。因此,根据检测到的睡眠状态反过来可以控制EPR,实施一种时间渐变的方案。
例如,无论何时所述装置检测一位患者处于觉醒状态,EPR控制启动,否则的话就会像前面描述那样开始运作,以降低呼气期间的压力。但是,一旦所述睡眠状态检测模块检测到患者进入或者NREM睡眠转变,NREM,或者REM睡眠时,所述EPR功能在呼气时会解除,以阻止降低呼气时的压力。
可选的,在许多呼吸循环或检测到的从一种觉醒状态到一种相关的睡眠状态的转变的一段时期,所述EPR控制功能根据压力治疗中的数量的下降逐渐降低而解除。然而,所述EPR功能的解除在睡眠起始时得到体现,它甚至还会选择性的在一次疗程中进行动态实施。例如,睡眠起始后,所述的睡眠检测模型检测到患者从睡眠状态返回到觉醒状态(如,通过进入觉醒状态并维持一个确定的时期),所述EPR会反复作用(或,通过在多次呼吸循环或一段时期内逐渐增加压力治疗的数量中逐渐反复。)所述EPR可能随后解除(如,逐渐地)正如前面讨论的更进一步检测患者睡眠状态。
B 呼吸事件报告
在一些实施例中,所述睡眠状态检测可以控制如何对检测的事件进行评分。例如,一些呼吸治疗仪可检测的事件,可能不会作为事件得分的一部分进行报告,即使这些得分是由仪器监测到的。例如,检测到呼吸事件,如窒息,呼吸不足,觉醒,面罩脱落,和/或口腔泄漏,同时也检测到了觉醒状态,则这些事件不作为一个特别事件的部分进行报告。
在这样的实施例中,处理装置的控制运算在评分呼吸事件之前,就可确认患者是在NREM睡眠转换,NREM睡眠或REM睡眠状态。可选的,所有的事件都进行评分,但是每个事件的评分需根据事件发生在检测到的不同的睡眠状态进行分类,或者更简单,根据事件发生是在清醒或者睡眠状态下进行分类。例如,一定数量的窒息可以为一个疗程的清醒时间评分,一定数量的窒息可以为一个疗程的睡眠时间(例如睡眠相关状态)分别评分。因此,所述得分事件在检测状态的基础上作为检测设备输出值进行报告,这样与他们发生时所处睡眠状态相关联。
治疗更改控制
在一些呼吸治疗器装置中,治疗压力可被仪器自动设置。例如,一些压力治疗装置可能对病人呼吸异常时的呼吸流量信号进行分析,并调整治疗压力。对于OSA患者,这些异常情况可能被认为是流速限制,打鼾和阻塞性窒息。一旦压力治疗装置发现了一次或多次这样的情况,便会根据情况的严重性,调整压力并上升到预定的最大值。
如果上升调整后异常情况被移除(例如长时间未检测到),压力治疗装置便会调整压力并降低到预定的最小值(有病人的医生所设定)。然而,在某些情况下,这并不是最理想的情况,因为上述的情况会连续并快速的发生,这将导致压力不必要的增加或降低。
因此,根据现有技术,这样的压力控制运算需被修改。例如,压力控制运算可被用于治疗检测到的异常情况而不是简单的降低压力以至于检测不到这些异常情况,在根据临界值比较睡眠稳定性指数的基础上,控制并把握增加的治疗程度直到前面所述的治疗稳定指数被储存在“高的/可接受的”水平。
在治疗压力增加后立即延伸临界值,压力可能被降低或者逐渐倾于向下直至预定的最小值。代表性临界值,例如“高的/可接受的”水平,可根据经验预测,例如,根据现存的临床资料得到数据。因此,用睡眠稳定性指数代替或呼吸异常检测缺失之外,降低治疗压力。
评价心肺结合(
CPC
)
CPC是将心率和呼吸率联合检测的一种有趣的方法。CPC有很多使用方法,其中最相关的是两样是能区别中枢和阻塞性窒息以及检测睡眠质量。呼吸可简单采用呼吸流量测量,而心率测定较难。一种根据呼吸流量测量心律的方法已经发明即根据全部的呼吸流量信号分析心原性流量成分。
例如,美国专利号为5,704,345的专利公开了检测心原性空气流的几种方法,并公开了相关的参考资料。心原性流量是由于心肺间的相互敲击会在呼吸流量信号中产生小小的振动。在疗程中他们不能轻易的被识别。特别地,在觉醒状态,NREM睡眠T状态和REM状态,很难被识别。而在稳定睡眠和NREM睡眠状态可被精确地检测到心源性流动信号。因此,一种自动检测心原性流量的方法相应地产生,即通过等待检测稳定睡眠产生和在睡眠的几个过程中,仪器可完成心原性空气流量的测定和/或估算心肺联合指标。
换言之,特殊睡眠状态(例如NREM睡眠状态)心原性流量的检测可推断或允许检测其他心原性流量,然而,其他状态(如NREM
T, REM或者觉醒)则会阻碍心原性流量的检测程序。因此,心原性空气流量可根据睡眠状态来检测。例如,在一个实施例中,根据睡眠状态,如果检测到区别于阻碍性窒息或中枢性窒息的无呼吸状态,心原性空气流量的测定过程可以被使用。例如,如果治疗仪器的控制器检测到呼吸暂停并且探测到该状态是稳定的(如NREM状态),则CPC探测器或过程控制器可根据探测到心原性空气流量(有或无心原性空气流量)以判断呼吸暂停的状态(中枢性或阻塞性)并采取措施。
周期性呼吸状态检测模型
在一些实施例中,睡眠状态检测装置102可运行周期性呼吸状态模块1880(如图18所绘模型)。就这点而言,在睡眠过程中伴随着呼吸不稳定可能存在睡眠片段。但事情往往不是这样,这是睡眠的周期特性。尤其是病人由觉醒向睡眠转变的时候,人体呼吸控制器失去规律而依次导致不稳定的呼吸形式。周期性呼吸状态检测模块就是为了辨别睡眠中的这一状态。在一些实施例中,这可通过检测处理从呼吸流量传感器传来的呼吸流量信号并输出周期性呼吸状态表明患者是否正在进行周期性呼吸(周期性呼吸状态变量1为是,0为否)。
例如,所述模块可以从特征提取装置110中分析一个或多个上述特征。所述实例的处理,需与图18中详述的模块相结合。通过检测呼吸流量,特征提取模块110可测定一个或多个上述讨论的特征。例如,在一些实施例中,将分析特征3mvTtot和IP75A。在前面的分析中得知,P75A为高于吸气峰流量的7 5%峰面积,而3mvTtot三分钟通气量和总呼吸时间的比率。
所述每个特征提供到临界值单元1810——拥有一个或多个临界值函数,如图19所阐述的例子,它适用于特征子集的每一种特征。例如,设定I75AV特征临界值,图表1950的函数可被运行,如输入IA、IB和IC则在表I75AV(图中所示特征仪1940A)中输出OA, OB, 和OC 。
正如前面这些表中所提到的,这些输入值IA、IB、IC等代表连续输出值范围。因此,对于表I75AV(图19中所示特征仪1940A)来说,IA表示从小于0.01上升到0.01并包括0.01的序列值,IB表示大于IA但小于等于0.035的序列值,IC表示大于IB但小于等于10的序列值,并被认为是输入值中的最大值。例如,例如,假设当I75AV特征具有确定值0.035的情况下,其输出值,可被认为是一个周期性状态权重, “1”作为I75AV特征被选择输出。
同样的,输出值被选为3mvTtot特征。然而,在这种情况下,输入值和输出值都是基于图表1950的函数的3mvTtot特性而来的。值得注意的是,图19表中特定的值,在这里只是举例说明。其它值可采用,或凭经验从无论是一群患者或一个特定的患者中测定。因此,这些值可以从仪器学习和模式识别技术中掌握。可以采用的技术实例包括聚类分析,支持向量机和贝叶斯分类。
然后,这些周期状态权重需经过原始周期状态计算器1882处理。而原始周期状态得分可通过如下方法获得:
IF (O_psTl == 1) 则 原始周期状态= 1
ELSE IF ( (O_psTl == 0) 和 (O_psT2==l) ) 则原始周期状态 = 1
其他时候 原始周期状态 = 0
其中:O_psTl 是周期状态权重根据I75AV特征获得。
O_psT2
是周期状态权重根据3mvTtot特征获得。
周期状态的输出信号(检测是否为周期性呼吸状态)可根据以下方法获得:
如果原始周期状态为1:
a. 设置周期状态输出信号为1
b. 提高周期状态计算器值为1,
c设置稳定睡眠计算器值为 0
如果原始周期状态为0, 检查稳定睡眠计算器是否低于计算器临界值(例如范围为5-15,为10更好)。
如果稳定睡眠计算器值小于等于计算器临界值(e.g., <= 10) 则:
a. 设置周期状态输出值为1
b. 提高周期状态计算器值为1,
c. 提高稳定睡眠计算器值为1
如果稳定睡眠计算器值大于计算器临界值(e.g., > 10) 则:
a. 设置周期状态输出信号为 0
b. 提高稳定睡眠计算器值为1
在这个实施例中,尽管合适的计算器临界值可被设置为10,其它值也可采用,例如,在一些实施例中,这些值可以从仪器学习和模式识别技术中掌握。可以采用的技术实例包括聚类分析,支持向量机和贝叶斯分类。
O.睡眠起始模块
精确地探测睡眠起始对病人睡眠形式和睡眠质量的检测非常有利。对精确的探测睡眠状态也非常重要。因此,在睡眠情况检测仪仪器102中,一个指标可以被应用于检测和计算在特殊睡眠过程中第一次从清醒到睡眠状态的转变,转变的睡眠状态可以是NREM睡眠T状态或者NREM状态。下面将讨论哪一种状态有更高的转变概率。
在实例应用过程中,睡眠起始检测模块可以利用输入数据代替以来自特征特征提取模块110的特征为基础的一个或多个呼吸流量。睡眠起始检测仪可以被用于分析设置以下特征:
1.) EPFL比率——这个特征是呼气最大峰值出现的位置和全部呼吸周期的比值,它可被定义为样品数(例如样品数量)而不是时间单位(例如秒)。
2.)
EPFL变异度特征——这个特征是根据现有的EPFL比率除去一段时期的EPFL比率平均值而得到的值。例如,这段时期的呼吸次数可以在10-20之间但最好为15。
3.) 初期睡眠状态得分——睡眠检测仪根据睡眠状态检测得到的输出值。
通过前面所述输入数据,睡眠起始监测仪在测定睡眠起始时可运行以下方法。这需要在患者进入循环呼吸的基础上进行,直至患者进入睡眠状态:
1.) 获得睡眠-清醒转变的指标(可调整为0):
如果(EPFL变异度特征>= Tl) 或者(原始睡眠状态得分>= T2) 则 (提高睡眠-清醒转变的指标。
2.) 如果(EPFL变异度特征>=
T3) 则睡眠起始 = 1 (i.e., 睡眠起始指标可用于显示患者使用睡眠情况检测模块转变进入睡眠第一时期)。
在这样的实施例中,Tl, T2 and T3的临界值可根据一群患者的经验值确定。也可根据一群特殊患者的来确定和设置仪器。
相应地,睡眠起始指标被认为是预先设置的二元开关或标志(如为0时)例如在启动时,显示睡眠的第一次起始未发生。在实例中,如果紧接其后的指标被设置为1,如有人使用检测仪则被认为处于睡眠状态(睡眠期间第一个睡眠起始)。否则,病人被认为在觉醒状态(在特殊期间使用检测仪后还未进入初始睡眠状态)。在这样情况下,仪器可以被用于辨别睡眠开始(最初转变至睡眠)以及之后在普通或者特殊期间从觉醒状态至睡眠状态的多重转变的睡眠转变。
P.呼吸治疗装置
尽管在检测和数据分析时,控制器可以无反馈运行,但是一些前面讨论的实例可以通过具有流量发生器的控制器实现。选择性地,如图11,睡眠情况检测仪可以应用于流量发生器1110,例如随动系统控制增压器可用合适的传感器控制(例如,压力传感器和/或流量传感器1106)。
呼吸治疗或压迫疗法方案,例如和CPAP治疗方案相关的治疗压力水平,可以通过仪器的控制器传送。治疗压力程度可自动调节以适应先前探测到的睡眠条件,使睡眠呼吸紊乱(SDB)的病人更舒适。其他压力调解方案也可以被应用。压力通过病人接口1108例如套管、插管和导管传递给病人。
在具体实施例中,显示器1120装在睡眠探测仪102的外表面。如图解所述,控制器104可与显示器安装在相同的位置。通过显示器1120,使用者或临床医师可随意看到仪器或者其他仪器、系统或电脑的仪器储存器转换而来的输出信号。因此,控制器可随意配置以联系其他的外部设备,例如,通过网络联系其他的电脑,通过网线接口交换设施(如,以太网卡或连接)(未显示)或无线交换设施(如蓝牙收发器)(未显示)。
控制器结构示例
图12的方块表举例描述了控制器104系统体系结构。在说明中,睡眠条件探测仪102与一个或多个程序控制的处理器1208可被应用于通用计算机。仪器可包括显示器接口1210以显示前面描述的模块的数据(睡眠状态、持续时间、睡眠质量指标、泄漏检测数据或/和稳定性指标等)或者显示器例如监视器或液晶显示器中所描述的结果或者图表。
使用者控制/输入接口1212,例如键盘、操作按钮、控制按钮、鼠标等也包括在内,正如前面所述,以调整流量发生器,输入数据,或者其他激活或者其他所描述的操作方法。所述装置还包括传感器或数据接口1214,如公共汽车一样,接受/传递数据, 例如程序指令、流动数据、压力数据、睡眠质量数据、睡眠状态数据、睡眠稳定性数据、觉醒数据和上述描述的其他模块输出或输入的信号。
所述装置也包括包含控制指令和上述方法和模型的数据存储器/数据存储组件1220。例如,1222包含流动信号处理,例如测量和/或特征提取的存储处理器控制指令。1224包含流速限制检测、AHI检测、泄漏检测、觉醒检测、睡眠状态检测、睡眠稳定性检测和/或睡眠质量检测的存储处理器控制指令。1226步骤,包含呼吸治疗控制,例如压力控制调整和反馈处理等的存储处理器控制指令等。最后,1228包含流动数据、检测呼吸特性、干扰数据、检测呼吸不全和窒息事件(AHI)、睡眠状态、稳定性和觉醒状态的事件和指数、睡眠期时间、睡眠质量指标、病人信息、报告和附图等等存储数据。
在一些实施例中,处理控制指令和控制上述方法的数据包含在电脑可读记录介质中,如通用计算机的应用软件,通过下载软件中的任意描述方法到通用计算机中,可使通用计算机担任专用计算机。
在上述的描述中的附图、专业术语、方程式和图表符号在现有技术中已经有彻底的公开。在某些情况下,术语和符号并非现有技术应用所必需的。例如,检测方法中的处理过程,根据相应的各种模块按字母排序已有详述,这样的排序和模块化是非必须的。
本领域的技术人员将指出这些调整可被修改和/或可被平行执行的特征。这些本领域的技术人员还认可 某些模块中的一些特征可以结合其它一些模块的特征应用于本技术中的非连续的特征。此外,尽管整个系统在图1中被详细描述,但是清楚的特征可单独,或不同组合被应用到其它呼吸治疗和/或检测系统。而且,尽管特殊值和阈值在表格中被描述,但是经验数据和机器学习的决定值也可应用。
因此,尽管本技术已通过具体的实施例进行了详细的描述,但是,可理解的是,这些实施例仅仅用于描述本技术的原理和应用。故需要大量修正以说明具体实例,并发明其他与科技紧密相连的技术。
应当指出的是,对于本领域的普通技术人员来说,在不脱离本发明构思的前提下,还可以做出若干变形和改进,这些都属于本发明的保护范围。