确定cpap治疗期间的泄漏的制作方法

专利名称：确定cpap治疗期间的泄漏的制作方法
技术领域：
本发明涉及呼吸暂停和其它呼吸紊乱的治疗。具体而言，本发明涉及 确定机械施加气道正压过程中的泄漏气流量和真实呼吸气流量的方法和装 置。
背景技术：
为了治疗呼吸暂停和其它呼吸紊乱，在可以在吸气期间较高而在呼气 期间较低的压力下，例如通过面罩从机械呼吸器或呼吸机提供可呼吸的气 体。(在本说明书中对"面罩"的任何提及应被理解为包括用于将可呼吸的气 体传送到人的气道的任何形式的设备，包括鼻罩、口鼻面罩、鼻塞/鼻枕和 气管内插管或气管切开插管。术语"呼吸机"用于描述完成呼吸的部分工作的 任何设备。)一般要在机械通气期间测量对象的呼吸气流量，以评定治疗的 适合性，或控制呼吸机的操作。
通常使用放置在面罩和压力源之间的气体输送路径中的呼吸速度描记 器来测量呼吸气流量。面罩和对象之间的泄漏是不可避免的。呼吸速度描 记器测量呼吸气流量加上泄漏流量加上通过通气口的流量(也称为"故意泄 漏(deliberateleak)")的总和。如果泄漏的瞬时流量已知，则可通过从呼吸速 度描记器处的流量减去泄漏流量和通过通气口的流量来计算呼吸气流量。 一般，通过通气口的流量是通气口处的压力的已知函数，且假定能以合理 的精确性估计通气口处的压力，于是可直接计算通过通气口的流量。另一 方面，如果通气口特征对所使用的泄漏模型适当，则可将通气口流量和非 故意泄漏聚集在一起，并估计为单一量。在下文中将假定使用通气口处的 压力来直接估计通气口流量，并且在没有明确提到时，假定已发生了将该 通气口流量从总的呼吸机流出量减去。
校正泄漏流量的一些方法假设(i)泄漏实质上是不变的，以及(ii)在足 够长的吋间内，吸入和呼出的呼吸气流量将抵消。如果满足这些假定，则 在足够长的一段时期内，通过呼吸速度描记器的平均流量将等于泄漏的数量，且可接着如所述地计算真实的呼吸气流量。
只有当面罩处的压力不变时，已知的方法才是正确的。另一方面，如
果面罩压力随着时间变化(例如，在呼吸机的情况下)，则上面的假设(i)将是 无效的，且所计算的呼吸气流量因此将不正确。这在图la-lf中显著示出。
图la示出在双水平CPAP(持续气道正压通气)治疗中测得的面罩压力的 迹线，压力在呼气时的约4 cm H20和吸气时的约12 cm H:jO之间。图lb 示出与面罩压力同步的真实的呼吸气流量的迹线。在时间=21秒处，面罩泄 漏出现，导致由于泄漏而产生的泄漏流量，其为治疗压力的函数，如图lc 所示。在图ld中示出的所测得的面罩流量现在包括由于泄漏流量而产生的 偏移。现有技术方法接着在多次呼吸后确定所计算的泄漏流量，如图le所 示。如图lf所示，得到的所计算的呼吸流量为所测得的流量减去计算的泄 漏流量，其已经返回到正确的平均值，然而在量值上并没有被正确地估算， 给出峰值正气流量和负气流量的虚假指示。
在欧洲公布号0 714 670 A2中公开了另一现有技术装置，包括计算压 力依赖性泄漏分量。此方法依赖于准确地知道吸气事件开始和下一吸气事 件开始的出现。换句话说，泄漏计算形成为在已知呼吸中的平均值，并应 用于随后的呼吸。
如果前面呼吸开始和结束的时刻是未知的，则不能使用该方法。通常， 可能很难精确地计算呼吸开始的时间。在紧接着泄漏突然变化之后的情况 中尤其如此。
而且，该方法在对象没有作出呼吸努力并且例如在呼吸暂停期间根本 没有立刻通气的情况下将不起作用，因为对于呼吸暂停的持续时间来说， 没有可在其间进行计算的呼吸的开始或结束。
在美国专利6，152，129(Berthon-Jones)中，通过首先从长期孔流量估计泄 漏路径的传导率(conductance)来确定泄漏
<formula>formula see original document page 5</formula>
其中Gfl/RL是传导率(L表示泄漏)，Q是瞬时流量，p是瞬时压力， 而o表示例如通过使用IIF或具有长时间常数的其它滤波器来低通滤波而 计算的长期平均值。注意，这里使用的词"平均值"包括包含低通滤波步骤的结果的一般意义，且不限于算术平均值或其它标准平均值，例如RMS平均 值。
基于通过孔的流量模型，瞬时泄漏流量于是为&=丄^。注意，瞬时
呼吸气流量因而为Qr=Q-Ql。 l &
Berthon-Jones试图通过以下步骤来处理瞬时泄漏流量的突然变化使 用模糊逻辑动态地调节滤波器的时间常数，如果确定泄漏是稳定的就延长 时间常数，如果确定泄漏突然变化就减小时间常数，并且如果在中间确定 泄漏是稳定的就在中间使用较长或较短的时间常数。
Berthon-Jones还通过模糊逻辑开发了干扰指数(jamming index),以处理 泄露传导率的较大和突然增加的情况，在这种情况下，被计算的呼吸气流 量将是不正确的。特别是在明显的吸气期间，被计算的呼吸气流量在比正 常吸气的预期持续时间长的时间内是大的正值。相反，如果有泄露传导率 的突然减小，则在明显的呼气期间，被计算的呼吸气流量在比正常呼气的 持续时间长的时间内是大的负值。
因此，导出了干扰指数，gP，确定泄漏突然变化的程度的指数，使得 气流量离零越远且量越大，指数就越大。在"29专利中描述了通过模糊逻 辑清楚地计算干扰指数，该专利在这里通过弓I用并入本申请。
低通滤波器的时间常数接着被调节以与干扰指数相反地变化。在操作 中，如果在泄漏中有突然和较大的变化，则该指数将较大，且用于泄露传 导率的计算的时间常数将较小，从而允许泄露传导率快速收敛到新值。相 反，如果泄漏在较长时间内是稳定的，则该指数将较小，且用于泄露传导 率的计算的时间常数将较大，实现了瞬时呼吸气流量的精确计算。在被计 算的瞬吋呼吸气流量较大并持续较长的时期的中间情况的谱中，指数将逐 渐变大，且用于泄漏的计算的时间常数将逐渐减小。例如，在不确定是否 泄漏实际上不变，以及对象仅仅开始大叹气，或是否实际上有泄漏的突然 增加的时刻，指数将具有中间值，且用于泄漏阻力的计算的时间常数也具 有中间值。

发明内容
本发明快速地确定CPAP系统中的瞬时泄漏，而不用详细地模拟泄漏 源，且不必确定呼吸周期中出现泄漏的精确相位。代替地，本发明依赖于限定呼吸周期的计时器的使用，以及确保瞬时流量与在长到足以包括整个 呼吸的一段时间内的流量比较的计算。在本发明，这通过向后查找以包括 整个相位周期来完成。这避免了必须采用在多次呼吸中的长期平均值，或 必须具有识别呼吸开始和结束的模型。
识别泄漏中的突然变化，且泄漏正在快速变化的程度表示为干扰指数
值，其接着用作参数来调节组成泄漏估计值的分量的份额(contribution),且 在伺服呼吸机的情况下，临时减慢或暂停对伺服呼吸机控制器输出参数(一 般是压力支持水平)的调节。
具体实施例方式
本发明以产生泄漏模型参数的持续升级估计值的愿望为动机，该泄漏 模型参数在实际泄漏参数稳定时非常稳定，但当实际泄漏参数快速变化时， 逐渐和适度地降级到较不稳定、较不精确但较快响应的估计值。所述泄漏 模型参数一般是孔常数(orifice constant)(相当于泄露传导率)，但不必需是孔 常数。
泄露传导率的持续升级估计(特别是在每次呼吸期间的持续升级)可通 过对非通气口流量(non-vent flow)(等于呼吸流量和瞬时泄漏流量的总和)和 对面罩压力的平方根执行某种低通滤波操作来计算，例如一般具有固定窗 口长度的一阶低通滤波器或移动的平均滤波器，产生泄露传导率估计值Gp 如在Berthon-Jones中所述。该方法具有优于一些其它方法的优点它不依 赖于呼吸相位(在当前呼吸中的位置)的确定。因此，可能出现泄漏的突然变 化，这使呼吸流量估计值以及因此呼吸相位估计值明显出错，然而泄漏参 数估计值的升级继续。这样的不依赖于呼吸的估计值的缺点是这些估计值 在呼吸中是不稳定的，除非出现特别偶然的事件，例如碰巧窗口平均滤波 器的持续时间精确地包括N次呼吸，其中N为整数。通过考虑当使用面罩 流量和面罩压力的一阶低通滤波器估计时的情况来说明这种不稳定性的特 别情况。为了简单起见，假定面罩压力是不变的，且真实的泄露传导率为 零。于是泄漏流量估计值仅仅为呼吸流量的一阶低通滤波形式。该估计值 在呼吸流量高于泄漏流量估计值时上升，并在呼吸流量低于泄漏流量估计 值时下降。特别是，在使用合理的滤波器时间常数的情况下，泄漏流量估 计值在大部分吸气期间上升，并在大部分吸气期间下降，在呼气暂停期间缓慢上升，以及在正常情况下在呼气暂停期间决定性地低于零。因为在呼 气暂停期间真实的呼吸流量为零，所以等于面罩流量(假定等于呼吸流量) 和被估计的泄漏流量之间的差的被估计的呼吸流量在呼气暂停期间为正，
比如等于a^ 。如果呼吸机设计成当被估计的呼吸流量超过某个真实的呼吸 流量阈值2,，,A^时触发吸气，则为了以期望的真实的呼吸流量触发，使用 这种泄漏估计的呼吸机必须将其触发阈值设定为较高的值仏,p一^+^ 。不 过，込p,.是呼吸流量形状的函数，如果不是不可能，也很难实时地确定呼吸 速率和低通滤波器时间常数。因此，触发实际出现在可变阈值处，且在最 坏的情况下，可能出现自动触发(在为零的真实的呼吸流量处触发)。应注意， 在产生失真的呼吸流量时，泄漏参数的非恒定估计的效应存在于整个呼吸 期间，且不管是否有可识别的呼气暂停，存在对周期(呼气触发)以及对(吸 气)触发的潜在的不利效应，以及对被估计的呼吸流量操作的其它算法。
如Berthon-Jones所述的，"干扰"是泄漏还没有得到补偿的程度，并通 常由泄漏的快速变化产生。它在这里被认为是模糊逻辑量。
如上所述计算泄露传导率估计值G,。注意，滤波器的时间常数优选地 随着干扰的增加而减小，如Berthon-Jones描述的。
在与呼吸频率比较仍然高的算术取样频率或较低的频率(例如10Hz)下 连续计算第二泄露传导率估计值G2。以下面描述的方式，该算法识别在当 前呼吸中的位置，然后试图找到与前面呼吸中的相同位置相关的时间。如 果它未能找到这样的位置，则它改为使用过去的时间10秒。在过去的这个 时间和当前之间，建立窗口。低通滤波的面罩流量和面罩压力的低通滤波 平方根(通过以不依赖与呼吸的方法滤波，例如一阶LPF卜一般为用于确定 Q的低通滤波值接着通过在该窗口内取平均值而被进一步低通滤波。这些 窗口平均值的比率是泄露传导率估计值G2，其在稳定泄漏的条件下是非常 稳定的。
因为G2对泄露传导率的变化响应相当慢，所以它不适合在泄漏变化快 时使用。因此若到了存在干扰的程度，则使用G,而不是G2。在优选实施方 案中，如果/是干扰依
中的量)，则使用由下式给出的传导率估计G^
接着通过aw=^7^:直接计算瞬时泄漏。确定G2的平均窗口开始的位置
目的是确定与患者在当前呼吸中位置相同的前面呼吸中的位置。为此， 我们需要呼吸相位的概念，该呼吸相位不仅仅是较小的一组类别中的一个， 例如吸气和呼气，而是概念上的实值(实际上为有理数)变量，其从吸气开始
到呼气结束逐渐增加，可能有小的有限数量的跳跃。在Berthon-Jones Cheyne-Stokes专利WO98/012965中提供了这样的概念，该专利在这里通过 引用并入本申请。在WO98/012965中，呼吸相位在吸气开始时被定义为0， 在呼气开始时为0.5，且在呼气结束时接近于l。给定这样的呼吸相位，我 们找到在当前时刻的呼吸相位，然后在时间上向后搜索以找到在前面呼吸 中相同的呼吸相位。虽然一般呼吸相位在呼吸期间随着时间增加，但因为 被如Berthon-Jones所述的系统估计的呼吸相位在呼吸期间不一定随着时间 增加(忽略呼吸相位必须降低的呼气到吸气转变)，所以必须有这样一种算 法，其在时间上向后搜索，使得具有与当前值相同的呼吸相位的相同呼吸 中的一点不被识别为在前面的呼吸中。下面描述这样的搜索算法；该算法 可能在例外情况下失败，但在大部分时间相当耐用。由于相位中的跳跃， 前面的呼吸中可能不存在这样的一个点，该点的相位等同于与当前时刻相 关的相位，因此取而代之地使用前面的呼吸中的其相位小于或等于当前时 刻的相位的最近的时间。
另一方面，为了在前面呼吸中找到与当前呼吸中相同的位置，使用传 统流量阈值来触发和周期化的系统不需要使用用于确定呼吸相位的模糊逻 辑系统。假定在吸气期间，在当前直到吸气结束之间的最大时间是已知的(一 般在吸气开始时确定，但不一定)，那么可给每个取样间隔处的呼吸相位增 加一个量，以便通过使该量相等地递增而使得相位在吸气结束时将达到0.5。 例如，在在吸气开始时确定最大吸气时间为1.6秒的简单情况下，相位将以 0.5/1.6相位单位/%!>的稳定速率增加。如果周期(到呼气的转变)在早些时候 出现，例如因为呼吸流量落到周期阈值之下，相位将在该点跳到0.5。类似 的考虑因素在呼气期间适用，相位增加的速率为1减去当前相位与到最大 呼气时间的剩余时间的商所得的差值。如果因为以这种方式确定的呼吸相 位一般仅用于确定与当前呼吸中相同的前面呼吸中的位置的目的，则它被 称为"簿记(book-keeping)"相位。
不管所使用的相位确定方法是Berthon-Jones的方法，如上所述的"簿记"相位，还是一些其它方法，都如下执行在时间上向后搜索，以在具有小 于或等于当前时刻处的相位的相位的前面呼吸中找到最迟的时间(但应认识
至IJ，对于"簿记"相位，可利用较简单的方法)。
以当前相位例如cpo开始，本发明在时间上向后査找在间隔[cpo-0.75， cpo-0.25]内最近的相位。目的是在时间上寻找在当前之前的呼吸的至少0.25 的点。当找到这样的相位时，本发明计算cp产(f)(r0.25，并向后查找在间隔 [cp广0.75，(p广0.25〗内的相位。这种操作继续下去， 一次0.25，艮p，化+1=9广0.25。 当找到在[93-0.75， cp3-0.25]内的相位时，重复停止，因为这恰恰是[(po-0.5, (p0]。 如果相位连续变化，则将准确地找到(Pq。实际上，最可能找到(P()-s，其中 希望s小。通过以这种方式继续，我们有些把握，相位向后而不是向前移 动，因为我们在当前之前找到在4个相位象限中的相位。该算法将0.5的两 个相位转变连续地视为向后移动，但实际方向当然实际上是不确定的。如 果该算法未能找到当前时刻和当前之前的时间之间内满足此标准的点，则 我们采用平均窗口的开始作为当前之前的某个合理的最大时间，例如10秒。 作为具体实施方式
，为了减少计算要求，泄漏流量值可在最后0.1秒内取平 均值并储存在附有相关的呼吸相位的缓冲器中，以便对最后呼吸的搜索在 IOO个点的缓冲器中执行，并在每0.1秒完成。接着可通过在最近的平均泄 露传导率估计值和正好在它之前的平均传导率泄漏估计值之间的线性插值 来计算在瞬时泄漏计算频率例如100Hz处的平均泄漏估计值。
在伺服呼吸机或使用某种通气测量来调节通气支持(例如呼吸流量的绝 对值的一半、或总肺泡通气量、或峰值流量、或在吸气或呼气期间确定的 流量的某个加权平均值)的其它系统中，当观察到干扰时，系统慢下来或中 止压力支持的变化。这是因为呼吸流量估计在干扰存在时不可靠，且基于 呼吸流量的各种通气测量可能过高估计通气，例如导致伺服呼吸机在通气 支持中没有保证的撤回，因为通气出现为高于目标通气。呼吸支持的调节 慢下来的程度优选地是干扰的某个递增函数。例如，如果从在以前时间计 算变化的呼吸支持中被计算的变化是某个值M，则在支持中被调节的变化 将为AAS，其中例如对于化O.l， t为l，对于^0.3， A:为O，且采用线性插 值作为/的中间值。
权利要求
1.一种确定患者呼吸设备中的面罩泄露传导率的方法，包括步骤确定相对快地响应于泄露传导率的变化的第一泄露传导率估计值，确定较慢地响应于泄露传导率的变化的第二泄露传导率估计值，以及确定作为所述第一泄露传导率估计值和第二泄露传导率估计值的函数的所述面罩泄露传导率。
2. 如权利要求l所述的方法，进一步包括确定所述面罩中的压力的步骤。
3. 如权利要求2所述的方法，进一步包括确定作为所述面罩泄露传导 率和面罩压力的函数的非故意面罩泄漏的估计值。
4. 如权利要求3所述的方法，其中所述非故意面罩泄漏被计算为-其中Q,eak为非故意面罩泄漏， Gj为面罩泄露传导率，以及Pmask为面罩压力。
5. 如权利要求1所述的方法，进一步包括确定无补偿泄漏指数的步骤， 所述无补偿泄漏指数表示面罩泄漏未得到补偿的程度。
6. 如权利要求5所述的方法，其中计算作为所述无补偿泄漏指数以及 所述第一泄露传导率估计值和第二泄露传导率估计值的函数的所述面罩泄 露传导率。
7. 如权利要求6所述的方法，其中所述面罩泄露传导率被计算为 其中Gj是面罩泄露传导率估计值， G,是第一泄露传导率估计值， G2是第二泄露传导率估计值，以及J是无补偿面罩泄漏指数。
8. —种向患者提供通气支持的方法，包括步骤 确定表示非故意泄漏未得到补偿的程度的指数；根据所述指数调节通气支持的水平。
9. 如权利要求8所述的方法，其中以作为所述指数的函数的速率调节 所述通气支持的水平。
10. 如权利要求9所述的方法，其中调节的所述速率随着无补偿非故 意泄漏的增加而降低。
全文摘要
在CPAP或正压通气系统中，在当单次呼吸建立平均窗口的开始时的时间内通过回顾在每个时刻确定泄漏，在该时间内对在确定泄漏模型的参数时使用的低通滤波值取平均值。干扰指数指示泄漏是否快速变化。到了干扰高的程度，所使用的泄漏估计从使用滑动呼吸窗口平均逐渐变化到更坚实和更快的响应低通滤波方法，且基于使用估计的呼吸流量进行测量来调节通气支持慢下来或停止。
文档编号A61M16/00GK101541366SQ200780031961
公开日2009年9月23日 申请日期2007年8月30日 优先权日2006年8月30日
发明者D·J·巴辛 申请人:雷斯梅德有限公司