JP5695598B2 - 閉塞睡眠時無呼吸の改良された流量制限検出の方法及び装置 - Google Patents

Description

本発明は、睡眠時無呼吸又は減呼吸を含む呼吸状態を診断及び治療する改良された方法及び装置に関する。具体的には、本発明は、呼吸気流曲線の吸息部分の時間展開の解析に基づき、患者の気道の部分的又は全面的閉塞を検出することに関する。本発明は、更に、ＣＰＡＰ人工呼吸装置によって適用される治療圧力を調節する装置に関する。

閉塞睡眠時無呼吸（ＯＳＡ）、即ち、無呼吸、減呼吸、及び猛烈ないびきを含む症状は、可能性として心臓疾患型を含む重大な健康問題をもたらす睡眠中断を引き起す。ＯＳＡは、人の気道通路部分の虚脱によって引き起こされる。ＯＳＡに選択される治療は、患者の気道を開いたままに保つ持続的気道陽圧法（ＣＰＡＰ）の施術である。通常、圧力範囲４〜２０ｃｍ Ｈ_２Ｏにある空気がモータ駆動送風機によって供給される。送風機の出力は空気配送機器を通過し、患者の気道へ密閉的に係合する。マスク、気管切開チューブ、気管内チューブ、鼻枕、又は他の適切な機器が使用されてもよい。空気配送機器に隣接した配送管の中に排気口が設けられる。二層ＣＰＡＰと呼ばれるＣＰＡＰ機器は、吸息及び呼息の息サイクルを感知し、吸い込み及び吐き出しの間に異なる正の圧力レベルを提供するものがある。幾つかの自己滴定ＣＰＡＰ機器は、ＯＳＡを解消する最小圧力を選択することによって、個々の患者について適切な圧力レベルを決定するように設計される。そのような機器において、いつＯＳＡが存在するかを決定するために呼吸パラメータのパターンが監視され、適切なパターンが検出されると、ＣＰＡＰ圧力が上げられて増加した気道圧力が提供され、理想的には、閉塞エピソード及び他の形態の息障害の生起が覆される。そのような機器は、米国特許第５，１４８，８０２号及び５，２４５，９９５号で説明されている。

典型的には、ＯＳＡを患っている人は、睡眠検査室で診断及び治療される。睡眠検査室では、最初の夜間睡眠期間中に不快の存在が確認され、２番目の夜間睡眠期間中に適切な治療圧力が決定される。生じる１つの問題は、人が異なる段階の睡眠を経るにつれて、適切な圧力が夜間中に変動することである。したがって、ＯＳＡの存在を決定し、圧力レベルを最適圧力へ修正することのできる簡単な装置を、患者自身の家庭で使用する必要があると長い間思われてきた。時間の関数としての気流曲線の形状解析からＯＳＡの存在を決定するため、幾つかの試行が行われたが成功はまちまちであった。

閉塞睡眠時無呼吸における上方気道圧力・流量関係の監視は、Ｓｍｉｔｈ ｅｔ ａｌ．，１９８８，Ｊ．Ａｐｐｌ Ｐｈｙｓｉｏｌ．６４：７８９−７９５で説明されている。この論文の図１は、増加する鼻圧力の変動レベルにおけるポリグラフ睡眠記録を示す。吸息容量流量は或る息で安定水準に達し、これは気流制限の存在を示唆することが特筆された。マスク圧力又は食道内の圧力に対する中間吸息気流を描くことによって、圧力・流量曲線が構築された。次に、平均中間吸息流量に対する鼻圧力の圧力・流量図表が最小二乗線形回帰によって当て嵌められ、虚脱可能場所の上流抵抗が計算された。

上方気道圧力・流量関係への正の鼻圧力の効果は、Ｓｃｈｗａｒｔｚ ｅｔ ａｌ．，１９８９，Ｊ．Ａｐｐｌ Ｐｈｙｓｉｏｌ．６６：１６２６−１６３４で説明されている。この論文の図４は、圧力・流量の記録が低圧力レベルで安定水準に達すること示す。更に、圧力が増加されたとき、流量は安定水準に達しなかったことが示された。

米国特許第５，３３５，６５４号（Ｒａｐｏｐｏｒｔ）は、ＯＳＡを患っている患者について気流対時間曲線へのＣＰＡＰの効果を示す。Ｒａｐｏｐｏｒｔの図１〜図５は、圧力が２ｃｍ Ｈ_２Ｏの階段で１０ｃｍ Ｈ_２Ｏから２ｃｍ Ｈ_２Ｏへ低減されるにつれて、気流対時間の曲線は、ほとんど滑らかな正弦パターンから、初期及び終期に流量スパイクを有して曲線の吸息部分が平坦化されたパターンへ変化することを示す。２ｃｍ Ｈ_２Ｏにおいて、曲線は、いわゆるＭ形（即ち、中間に脈動波を有する）を展開し、更に平坦化された中間領域の各々の端に行き過ぎ（即ち、ピーク）を展開した。閉塞を示す流量形状を特性化する試みの中で、Ｒａｐｏｐｏｒｔは流量制限及び／又は部分閉塞パターンを示すと称される幾つかの指数を列挙している。これらの指数の中には、（１）流量信号の導関数がゼロであること、（２）延長された区間について流量信号のピーク間の第２の導関数がゼロであること、（３）中間吸息流量に対する早期吸息流量の比が１より小さいか等しいこと、が含まれる。この特許は、更に、閉塞の指標と称される事象を列挙している。即ち、（１）ピーク吸息流量をピーク呼息流量へ接続する線の傾斜が低減すること、（２）流量信号の上向き又は下向き行程（ｄＶ／ｄｔ）が急峻であること、及び（３）呼息流量に対する吸息流量の比が０．５を超えることである。

ＣＰＡＰ治療の管理に関しては、閉塞を示す異常息パターンを感知及び検出する様々な手法が公知である。例えば、米国特許第５，２４５，９９５号（Ｓｕｌｌｉｖａｎ ｅｔ ａｌ．）は、睡眠中の吸息及び呼息圧力を測定することによって、いびき及び異常息パターンを検出し、閉塞前エピソード又は他の形態の息障害の早期表示を導く方法を説明している。具体的には、呼吸パラメータのパターンが監視され、既定のパターンが検出されると、ＣＰＡＰ圧力が上げられて増加気道圧力が提供され、理想的には、閉塞エピソード及び他の形態の息障害の生起が防止される。

米国特許第５，６４５，０５３号（Ｒｅｍｍｅｒｓ）は、平坦度指数の計算を説明する。この場合、平坦度とは、平均気流からの観察気流の相対的偏差であると定義される。Ｒｅｍｍｅｒｓにおいて、気流の個々の値は、吸息期間の４０％と８０％との間で取得される。吸息流量の個々の値から、平均値が計算及び減じられる。個々の差は二乗され、観察総数マイナス１によって割られる。この結果の平方根は、相対変分を決定するために使用される。相対変分は平均吸息気流によって割られ、当該息について相対変分又は変分係数が与えられる。

米国特許第５，７０４，３４５号（Ｂｅｒｔｈｏｎ−Ｊｏｎｅｓ）は、患者の監視された呼吸気流の吸息部分の平坦化をパラメータ化する２つの閉塞指数値を計算することによって、患者の気道の部分閉塞を検出する方法を開示する。次に、いずれかの閉塞指数が閾値と比較される。第１の形状因子は、息サイクルの吸息気流の中間部分の平均と吸息気流の平均との比を包含する。

ここで、ｆｓ（ｔ）は患者吸息気流の標本であり、Ｍは

によって与えられる吸息気流の平均である。

第２の形状因子は、次式に従って、吸息気流の中間部分のルート平均平方偏差と平均吸息気流との比を包含する。

米国特許第６，８１４，０７３号（Ｗｉｃｋｈａｍ）は、吸息気流に基づいて或る形態の閉塞を検出する方法及び装置を開示する。この方法において、中間吸息に対応する吸息気流標本が同定される。１つの実施形態において、吸息気流標本が、閾値レベル、例えば、平均気流よりも小さいか大きいかに基づいて、加重因子が適用される。他の実施形態では、息における標本の時間位置に基づいて、異なる加重因子が標本へ適用される。吸息中の或る事象に先行して取られた標本、例えば、吸息の中間点に先立つ標本は、事象に後続する標本よりも低い加重因子を割り当てられる。次に、対応する加重因子を有するこれらの標本を使用して、閉塞指数が計算される。

これらの開示された全ての手法は、或る型の流量パターン、特に、Ｍ波パターンにおける流量制限を検出することはできない。更に一般的には、先導又は遅延する行き過ぎを呈示す吸気パターンにおける流量制限が適切に検出されない。検出可能度は、２つの因子、即ち、（１）行き過ぎが吸息気流の中間部分にわたる範囲、及び（２）行き過ぎの大きさ、に依存する。もし行き過ぎの大きさが比較的に小さく、中間吸息の大きな部分に広がらず、吸息気流の後の半分に存在しなければ、Ｗｉｃｋｈａｍによって呈示された解決法は満足に働く。しかしながら、Ｗｉｃｋｈａｍ方法からの平坦化指数は、時には、Ｂｅｒｔｈｏｎ−Ｊｏｎｅｓによって開示された方法からの指数ほど正確ではない。本発明の１つの態様は、米国特許第５，７０４，３４５号（Ｂｅｒｔｈｏｎ−Ｊｏｎｅｓ）で言及されたアルゴリズムを単純化することである。この米国特許の開示は、参照して組み入れられる。こうして、本発明の目的は、言及された制限、特に、行き過ぎを有するＭ波パターンが存在するときに流量制限を検出できる方法を呈示することである。

従来技術によって説明された流量制限検出／推定手法は、更に、ディジタル処理能力及び要求される流量検出正確度の点で費用が高い。したがって、低コストの電子／ソフトウェア・プラットフォームにおける手法の利用性が制限される。したがって、本発明の他の目的は、アルゴリズムを単純化して、低価格帯の電子及びソフトウェア・プラットフォームがアルゴリズムを受け入れやすくすることである。

本発明は、吸息気流の平坦化に基づいて、呼吸閉塞の幾つかの形態を検出する改良方法及び装置に関連する。

本発明の１つの態様は、患者の気道の部分閉塞を検出する方法である。この方法は、患者からの呼吸気流を測定するステップと、気流の吸息部分を検出するステップと、吸息部分を整形して呼息休止を除去するステップと、派生物（ｏｆｆｓｈｏｏｔ）の効果を正規化するステップと、吸息流量を参照値へ尺度変更するステップと、吸息気流の中心部分の平均偏差を計算するステップと、平均偏差から平坦化指数を計算するステップと、移動平均フィルタを用いて平坦化指数をフィルタリングするステップとを備える。

部分閉塞を検出する方法に関しては、行き過ぎ（ｏｖｅｒｓｈｏｏｔ）の正規化補正が行われる。正規化補正が行われない行き過ぎは、吸息気流の円形推定を誇張する。

本発明の実施形態において、行き過ぎの正規化は、吸息気流の中心部分の外側にあって吸息気流の中心部分の平均値を超過する気流値を打ち切ることによって達成される。

本発明の他の実施形態では、吸息気流波を複数の部分へ分割し、各々の部分について平均振幅を計算し、平均振幅を解析して行き過ぎの存在を決定し、もし行き過ぎが示されたならば、平坦化を低として取り扱うことによって、行き過ぎが評価される。

本発明の更なる実施形態において、最大平坦度の探求は、吸息波形の中間部分の各々の可能な標本について平坦化指数の集合を計算し、吸息波形の平坦化指数を、計算された指数の最小値とすることによって実行される。

こうして、本発明の態様は、ピーク流量の位置が吸息流量曲線の中心から外れる範囲によって、患者の息における流量制限の範囲を推定することである。

本発明の更なる態様は、息の異なる部分の相対平均振幅が、正常な息によって示されるパターンから外れる範囲によって、流量制限を推定することである。

本発明の更なる態様は、最大平坦度を示す息部分を探求することによって、流量制限の範囲を推定する方法を提供することである。

本発明の目的は、患者の部分的又は全体的気道閉塞が、患者の気道の部分閉塞を検出する方法により呼吸気流を解析することによって検出される装置を提供することである。前記方法が、患者からの呼吸気流を測定するステップと、この気流の吸息部分を検出するステップと、吸息部分を整形して呼息休止を除去するステップと、行き過ぎの効果を正規化するステップと、吸息流量を参照値へ尺度変更するステップと、吸息気流の中心部分の平均偏差を計算するステップと、平均偏差から平坦化指数を計算するステップと、移動平均フィルタを用いて平坦化指数をフィルタリングするステップとを備える。

更なる目的は、追加のコンポーネントを使用することなく、又は、既存の呼吸装置の構造へ実質的な変更を施すことなく、気道閉塞を検出する改良アルゴリズムが実現される装置を提供することである。

したがって、患者の呼吸気流が継続的に監視され、吸息に関連づけられた呼吸気流部分が同定及び標本化される呼吸装置が提供される。

１つの態様において、本発明は呼吸装置に関連する。この呼吸装置は、加圧されて呼吸可能な空気を患者へ選択的に提供するようになっている空気源と、患者からの呼吸気流を感知し、気流を示す気流信号を生成する流量検出機器と、患者の気道の部分閉塞を検出する方法であって、患者からの呼吸気流を測定するステップと、この気流の吸息部分を検出するステップと、吸息部分を整形して呼息休止を除去するステップと、派生物の効果を正規化するステップと、吸息流量を参照値へ尺度変更するステップと、吸息気流の中心部分の平均偏差を計算するステップと、平均偏差から平坦化指数を計算するステップと、移動平均フィルタを用いて平坦化指数をフィルタリングするステップを含む方法を実現し、前記流量センサへ結合された閉塞検出器と、及び圧力センサ等価物へ接続された圧力コントローラとを含む。この圧力コントローラは、空気源の動作を制御し、閉塞信号を受け取り、閉塞信号に応答して空気源の動作を変更するように、物理的又は論理的に配列される。流量検出機器の用語は、流量センサ又は代替の流量検出機器又はアルゴリズム、例えば、モータ電流から流量を決定するものを含む一般的な用語であることを意図される。

本発明の他の態様は、睡眠障害を有する患者を監視及び／又は治療する装置に関係する。この装置は、患者の呼吸を感知し、対応する流量信号を生成する流量検出機器と、及び流量検出機器へ結合され、患者の気道の部分閉塞を検出する方法を実現するようになっている閉塞検出器とを含む。前記方法が、患者からの呼吸気流を測定するステップと、この気流の吸息部分を検出するステップと、吸息部分を整形して呼息休止を除去するステップと、派生物の効果を正規化するステップと、吸息流量を参照値へ尺度変更するステップと、吸息気流の中心部分の平均偏差を計算するステップと、平均偏差から平坦化指数を計算するステップと、移動平均フィルタを用いて平坦化指数をフィルタリングするステップを備え、閉塞検出器は、平坦化指数に基づいて気道閉塞を示す信号を生成する信号生成器を含む。

本発明の更なる態様は、睡眠障害を有する患者を治療する装置に関係する。この装置は、マスク、このマスクを介して、加圧されて呼吸可能な空気を患者へ選択的に供給する気体源と、気流を感知し、呼吸を示す流量信号を生成する流量検出機器と、流量検出機器へ結合され、患者の気道の部分閉塞を検出する方法であって、患者からの呼吸気流を測定するステップと、この気流の吸息部分を検出するステップと、吸息部分を整形して呼息休止を除去するステップと、派生物の効果を正規化するステップと、吸息流量を参照値へ尺度変更するステップと、吸息気流の中心部分の平均偏差を計算するステップと、平均偏差から平坦化指数を計算するステップと、移動平均フィルタを用いて平坦化指数をフィルタリングするステップを含む方法を実現するようになっている閉塞検出器と、閉塞信号を受け取り、応答して気体源を活性化する指令を生成するコントローラとを備える。

本発明の他の態様は、患者の気道の中の閉塞を検出する方法に関係する。この方法は、患者の気流を測定するステップと、前記気流の所定の区分を検出するステップと、患者の気道の部分閉塞を検出する方法であって、患者からの呼吸気流を測定する工程と、この気流の吸息部分を検出する工程と、吸息部分を整形して呼息休止を除去する工程と、派生物の効果を正規化する工程と、吸息流量を参照値へ尺度変更する工程と、吸息気流の中心部分の平均偏差を計算する工程と、平均偏差から平坦化指数を計算する工程と、移動平均フィルタを用いて平坦化指数をフィルタリングする工程を備える方法を実現するステップを含む。

本発明に従って構築された呼吸装置のブロック図を示す。 図１の装置の動作を図解するフローチャートを示す。 健康人及び部分気道閉塞を有する人について、典型的な呼吸信号の吸息段階を示す。 正常な息サイクルを示す。 正常及び流量制限の流量対時間曲線を示し、左方のパネルは正常曲線を示し、右方のパネルは流量制限曲線を示す。 ２つの可能な息について、時間に対する吸息流量を示す。 Ｍ形指数の計算に使用される基底ベクトルのグラフである。 典型的なＭ形息のグラフである。

装置及び方法論
図１は、本発明に従って構築された例示的呼吸装置１０を示す。呼吸装置１０は、可撓管１６によって送風機１４へ接続されたマスク１２を含む。マスク１２は、患者へ嵌められ、鼻マスク又はフェースマスクであってもよい。送風機１４は空気出口２２を有し、空気出口２２は、サーボコントローラ２０からの制御信号に従って、モータ１８によって駆動される。この配列によって、呼吸装置１０は加圧された空気（又は、図示されない源から酸素を強化された空気）を配送することができる。加圧された空気は、管１６によってマスク１２へ配送される。管１６は狭い排気口２６を設けられ、患者によって吐き出された空気は排気口２６を介して追い出される。

所定の規準を使用してサーボコントローラ２０及びモータ１８の動作を制御し、それによって装置１０の動作モードを規定するため、制御回路２４が使用される。好ましくは、本発明によれば、制御回路２４は、装置１０を動作させてＣＰＡＰを患者へ提供するように適応される。

制御回路２４は流量限定要素２８を含む。管３０及び３１は、それぞれ限定要素２８の異なる側からマスク圧力トランスデューサ３２及び微分圧力トランスデューサ３４へ導かれる。マスク圧力トランスデューサ３２は、更に、他の管３３を介して微分圧力トランスデューサ３４へ接続される。

マスク圧力トランスデューサ３２は第１の電気信号を生成する。第１の電気信号は増幅器３６によって増幅され、マスク１２の中の空気圧力に比例した出力Ｐ（ｔ）が生成される。この出力はサーボコントローラ２０へ直接送られる。

微分圧力トランスデューサ３４は、流量限定要素２８を横切る微分圧力を感知する。この微分圧力は、流量限定要素２８及び管１６を通る気流速度に関連する。微分圧力トランスデューサ３４は第２の電気信号を生成する。第２の電気信号は増幅器３８によって増幅される。この増幅された信号Ｆ（ｔ）は気流信号と呼ばれる。なぜなら、Ｆ（ｔ）は管１６を通る気流を表すからである。代替として、流量はモータ電流から引き出されてもよい。

気流信号Ｆ（ｔ）はフィルタ４０へ送られる。フィルタ４０はプリセット範囲内で信号をフィルタリングする。フィルタ４０及び増幅器３６の出力はＡＤＣ（アナログ／ディジタル）変換器４２へ送られる。変換器４２は、対応する信号をマイクロプロセッサ４４へ生成する。マイクロプロセッサ４４はアナログ制御信号を生成し、このアナログ制御信号は、ＤＡＣ４６によって、対応するディジタル制御信号へ変換され、サーボ２０の参照信号Ｐｓｅｔ（ｔ）として使用される。

呼吸装置１０の１つの動作方法は、図２のフローチャートで示される。当業者は、本明細書で開示される改良された流量平坦化指数を利用する他の方法論を認識する。図２の方法論の実施形態は、更に、米国特許第５，７０４，３４５号（‘３４５特許）で詳細に説明される。最初のステップ１００は、時間にわたる呼吸流量を測定する。この情報はステップ１０２で処理され、後続処理で質的尺度として使用される指数値が生成される。こうして、ステップ１０２は、本明細書で開示される平均方法に基づいて、閉塞指数値を生成することを含む。ステップ１０４は、息指数と閾値とを比較することによって、無呼吸が起こっているかどうかを検出する。

もしステップ１０４の答えが「イエス」であれば、無呼吸が進行中であり、Ｗｉｃｋｈａｍ又はＢｅｒｔｈｏｎ−Ｊｏｎｅｓ特許で開示される方法によって、ステップ１１０の開存性の決定が続く。もし気道の開存性が存在すれば、開放気道を有する中心無呼吸が起こっており、もし所望されるならば、この事象はステップ１１２で記録される。もしステップ１１０の結果が、気道が開存していないことであれば、全閉塞無呼吸又は閉鎖気道を有する中心無呼吸が起こっており、これはステップ１０８におけるＣＰＡＰ治療圧力の開始又は増加を結果する。もし所望されるならば、ステップ１０８は、検出された異常の任意的記録を含んでもよい。代替として、平坦化は、無呼吸が起こるとき単純に拒絶される。なぜなら、無呼吸が存在するとき、息パターンは、標本数、又は要求される一回呼吸量 − これは平坦化のための有効な息として要求される −のいずれかの規準を満足させないからである。

もしステップ１０４の答えが「ノー」であれば、１つ又は複数の閉塞指数、例えば、改良された流量平坦化指数が、ステップ１０６で閾値と比較される。それによって、気道閉塞の決定が取得される。ステップ１０６で、もし答えが「イエス」であれば、部分閉塞が存在し、もし「ノー」であれば、閉塞は存在しない（正常）。

気道の完全又は部分閉塞の場合、ステップ１０８が適用され、ＣＰＡＰ治療圧力の増加が続いて起こる。閉塞を有しない正常息の場合、無呼吸の生起を回避又は少なくとも低減するために要求される最小圧力を設定しようとする通常の方法論に従って、ＣＰＡＰ治療圧力が低減される。ステップ１０７における低減量は、もし所望されるならば、ゼロであってもよい。同様に、開存性気道を有する中心無呼吸の事象（ステップ１１０、１１２）では、治療圧力は増加されない。そのような圧力増加は、反射的に息を妨げ、更には息障害を悪化させる。

改良された流量平坦化指数
図３は、典型的な息サイクルの吸息部分に関する気流信号を描写する。健康人の息サイクルの吸息部分では（実線）、気流は吸息と共に滑らかに上昇し、ピークに達し、滑らかにゼロへ降下する。しかしながら、部分閉塞気道を有する患者は（ダッシュ線）、吸息の間の顕著な平坦ゾーンによって特性化された息パターンを示す。理論的には、閉塞された流量の場合、部分閉塞の程度が増加するにつれて、吸息気流信号は方形波になる傾向がある。

息における吸息及び呼息半サイクルの検出
図４で示されるように、正常な呼吸サイクルは吸息、呼息、及び短い呼息終期休止から成り立つ。呼息終期休止の間に、心拍に起因して気流の小さな振動が起こる。吸息のスタートを正しく検出するため、心拍に起因する空気の小さなパルスを無視することが必要である。同様に、嚥下に起因する短い停止を無視することが必要である。

吸息のスタートは、暫定的に、吸息気流が約０．０５Ｌ／秒の閾値を超過した時点とされてもよい。全容量が約０．１リットルに到達する前に、もし流量が再びゼロよりも下に降下するならば、呼息は継続していると仮定される。

呼息のスタートは、暫定的に、呼息気流がゼロよりも下に降下する時点とされる。呼息された全容量が約０．１リットルへ到達する前に、もし気流が吸息方向で約０．０５Ｌ／秒よりも上に進むならば、吸息が継続していると仮定される。

上記のアルゴリズムの目的は、各々の呼吸半サイクルのスタートと終わりを検出することである。この情報は形状検出器に必要である。この情報は、嚥下及び心原性気流（これらは０．１リットル容量の要件を達成しない）のような事象から影響を受けない。正常な息は約０．３Ｌ／秒のピーク流量及び０．５リットルの一回呼吸量を有する。こうして、正常な息を検出するための３〜５倍の予備能力が存在する。タイミング情報は、非常に小さい息は関係しない形状アルゴリズムに使用されるから、極端に小さい息に応答する必要はない。閾値交差点は、更なる低域フィルタリングなしに検出されるから、アルゴリズムは位相の遅延を導入しない。

曲率又は平坦化指数の計算
正常な吸息流量・時間曲線は円形又は疑似正弦の形状である。気流を生み出すために利用可能な横隔膜の筋肉運動は、息の始まりと終わりで低く、息の中間で高い。十分なＣＰＡＰ圧力を用いると、気道は本質的に剛性管として働き、流量は運動の関数となり、増加する運動は、増加する流量を生み出す。

もしＣＰＡＰ圧力が、無呼吸及びいびきを防止するには十分であるが、完全開存性を達成するには不十分であれば、気道は柔らかい伸縮管のように振る舞い、増加する運動は大気常態値以下の圧力を増加的に引き起こし、気道の累進狭窄を導く。或る閾値運動よりも上では、気道は狭くなる。即ち、増加狭窄は増加運動を正確に補償し、流量は運動から独立して一定になる。したがって、流量対時間曲線は、図５の中間パネルで示されるように、方形波を近似する。

実際には、もし運動が、運動が最大になる息の中間部分の間でのみ閾値平坦化振る舞いよりも上にあれば、平坦化は吸息の中間半分にわたって最も目立つ。逆に、もし運動が非常に高いならば、慣性運動が息の早期に高い流量の短い期間を導き、図５の右方枠で示される初期行き過ぎを引き起こす。こうして、本発明のアルゴリズムは、息全体ではなく、特に吸息の中間半分に注目する。

図５において、吸息流量対時間曲線は、単位持続時間及び面積へ尺度変更される。したがって、単位流量からのＲＭＳ偏差は吸息時間の中間５０％にわたって計算される。左方のパネルは、ＲＭＳ偏差が大きい正常曲線を示し、右方のパネルは、小さいＲＭＳ偏差を有する流量制限曲線を示す。

流量・時間曲線の中間吸息平坦化の程度を定量化するため、次のステップを使用して曲率又は平坦化指数が引き出される。

１．図６を参照すると、吸息半サイクルに含まれた呼息終期休止は、次のようにして削除される。即ち、流量がピーク吸息流量の７５％へ最初に到達する場所「ａ」を発見し、流量がピークの２５％へ最後に到達した点「ｂ」を逆方向に探索し、逆方向に外挿して、半直線「ａｂ」がｘ軸に到達した時間を突き止める。即ち、流量がゼロになるべき時間「ｃ」を発見する。これは、吸息がスタートした最も早い時点とされる。更に、もし「ｃ」が原点「ｄ」の左方にあれば、原点「ｄ」を使用する。

２．流量・時間曲線に対する０．１Ｈｚ低域フィルタの近似効果は、次のように計算及び反転されてもよい。（ａ）同一の０．１Ｈｚ低域フィルタへ波を２回通す。（ｂ）１回目にフィルタリングされた波から２回目にフィルタリングされた波を減じて、フィルタ効果の一次近似値を与える。（ｃ）この差を元の波へ加える。代替として、この漏れフィルタに対する補償は、結果を顕著には変更しないと予期してもよく、したがって計算効率を目的として除去されてもよい。

対象の息の場合、吸息持続時間は約２秒である。０．１Ｈｚフィルタは波形状の或る歪みを生み出す。例えば、２秒方形波の立ち下がり縁は、振幅の２０％だけ落ち込み、曲線の下の面積は１０％だけ低減される。補正の後、立ち下がり縁は振幅の２％だけ落ち込み、曲線の下の面積は０．７％へ補正される。フィルタは正弦波の下の面積を１０％だけ低減し、補正の後、面積は０．６％へ補正される。これは発明者らの現在の目的に十分である。

３．行き過ぎの効果は、下記で説明されるように正規化される。

４．吸息流量・時間曲線は、単位長さ及び単位平均高さへ尺度変更される。尺度変更は、発明者らが流量・時間曲線の振幅又は持続時間ではなく流量・時間曲線の形状（円形対平坦化）に関心を有するからである。正しいＣＰＡＰ圧力を用いる場合でも、振幅及び持続時間は非常に変動する（具体的には、ＲＥＭ睡眠中だけでなく、段階１の睡眠及び吐息を伴う目覚め、など）。したがって、長さ及び振幅はＣＰＡＰ圧力の妥当性に関して重要な情報を伝達しない。

５．甚だしく異常な形状を有する息は、テンプレートと比較することによって検出される。中間半分の上の任意の点で、閾値よりも大きくテンプレートから異なる息は拒絶される。デフォルトの閾値は１．０単位である。

６．波の中間部分について、参照振幅（Ｍ）からの吸息波の平均偏差が計算される。例えば、もし中間部分が吸い込み波の１／４であって中心にあれば、平均偏差は次式のようになる。

７．平坦化指数（ＦＩ）は、次式のように計算される。

実際の患者の場合、平坦化指数の典型的な値は正常曲線について約０．２１を与える。０．１５は閾値であり、それよりも下で機器はＣＰＡＰ圧力を増加すべきである。

８．もし息が拒絶されなければ、複数息点別移動平均流量・時間曲線が、例えば、５息平均で更新される。

ステップ８の目的（５つの息にわたる波形の平均化）は、心原性気流の効果を低減することである。心原性気流は、そうでないときの平坦化された流量・時間曲線を偽装する。嚥下、咳、会話などは、非常に異常な形状の息を生み出す。ステップ５の目的は、そのような息が５息平均に含まれるのを防止することである。そのような息は、真性の流量制限の認識を遅延させる。

ステップ８は必要でないかも知れない。なぜなら、流量制限の間、心原性圧力振動は、運動依存性に起因して、流速に重要な変化を引き起こさないからである。この場合、ステップ５も必要でない。なぜなら、咳、嚥下などは吸息の中間半分で平坦でなく、したがって非流量制限として正しく分類されるからである。

ステップ３、即ち、行き過ぎの正規化は、幾つかの実施形態で達成されてもよい。本発明のこの態様は、（１）ピーク流量は波の中間部分の近くで起こること、（２）波の最初の１／３は一般的に正の勾配を有すること、（３）最後の１／３は一般的に負の勾配を有すること、（４）これらの条件からの偏差は行き過ぎを示すこと、と仮定する。

１つの実施形態において、方法は、（１）波の中間部分、例えば、波の中間の１／３でピーク流量を推定するステップと、（２）その部分の外側、例えば、前記波の最初及び／又は最後の１／３において、前記ピーク流量を超過する任意の点が、ピーク流量に再設定されるステップとを含む。したがって、この方法は行き過ぎの影響を低減する。

第２の実施形態において、方法は、（１）波の中間部分（例えば、１／３）でピーク流量を推定するステップと、（２）その部分の外側（例えば、前記波の最初及び／又は最後の１／３）でピーク流量を超過する点が、前後の標本の平均に再設定されるステップと、（３）ピーク流量よりも大きい外側部分の標本がなくなるまでプロセスを継続するステップとを含む。この方法は再び行き過ぎの影響を低減し、同時に波形の形状を維持する。

第３の（好ましい）実施形態において、方法は、（１）波をｎ個の部分に分割し、各々の部分について平均振幅を計算するステップと、（２）パターンが、どのような種類の波を示すかを決定する論理を適用するステップと、（３）もし行き過ぎが示されたならば、平坦化を低にマークするステップとから構成される。

第４の実施形態において、方法は、（１）波の中間半分においてｎ個の標本ごとに平坦化の推定値を計算する（例えば、標本１６．．２４、１７．．２５、などについて平坦化を計算する）ステップと、（２）全ての標本順列について取得された最低値に対応する値として、平坦化指数をマークするステップとから構成される。

Ｍ形検出の代替の統計計算
本発明の代替の実施形態は、異なる方策を使用して、吸息波パターンがＭ形を有するかどうかを決定する。手法は、吸息波を特徴ベクトル（浮動小数点数）及び直交する基底ベクトル行列対へ分解することである。時間区間を正規化して、吸息がゼロで始まってＢラジアンで終わるようにし、基底ベクトルはｓｉｎ（Ｂｔ）及びｓｉｎ（３Ｂｔ）であるように選択されてもよい。こうして、吸息波形は曲線ａ ｓｉｎ（Ｂｔ）＋ｂ ｓｉｎ（３Ｂｔ）へ当て嵌められる。ここで、ａ及びｂは、特異値分解を使用し、最小二乗の意味で基底ベクトルを吸息波形へ当て嵌めることによって決定される。これは、基底行列の疑似逆を発見し、この疑似逆数を吸息波形へ乗じて因子を発見することと等価である。疑似逆は、オフラインで、即ち、人工呼吸装置に格納されたＣＰＵによらずに、計算可能である。

この統計計算に含まれる計算の量を低減するため、多くの患者について［ａ，ｂ］空間をマップし、個々の患者を表す点が落ちる集団を記録することが可能である。Ｍ形息は、非Ｍ形息から離れて集まることが発見される。こうして、息がＭ形であるかどうかは、Ｍ集団の中心からのユークリッド距離測定によって決定可能である。

Ｍ形指数計算の例
各々の吸息は、Ｎ個、好ましくはＮ＝６５の点の格子の上で補間される。図７で示される２つの基礎関数が、ｔ＝ｉ／（Ｎ−１）について計算される。ここで、ｉは０からＮ−１へ進む。
Ｂ１＝ｓｉｎ（πｔ）
Ｂ２＝ｓｉｎ（３πｔ）

次に、これらの基礎関数は、Ｍ形特徴の全ての後続計算で使用されるように記憶可能である。

次に、各々の吸息は、Ｎ個の点の格子の上で抽出及び補間される。次に、２つの因子が次のように計算される。
Ｆ_１＝ｓｕｍ（Ｂ１・ｆｓ）
Ｆ_２＝ｓｕｍ（Ｂ２・ｆｓ）
ここで、ｆｓは補間された吸息点を表し、・は点乗積演算子である。

最終の形状値は、次のように正規化することによって取得される。

次に、この形状因子は、０（純粋に正弦）と１（非常にＭ形）との間で変動するように制限される。

図８で示されるように、典型的なＭ形息の場合、上記の計算は次の結果を与える。
Ｆ１＝４．６０８２
Ｆ２＝２．６５３８
形状指数＝０．５０である。これは、単に約０．２のＭ形指数を有する典型的な非流量制限息と対比されてもよい。

本発明は特定の実施形態を参照して説明されたが、この実施形態は、本発明の原理の適用を例証するだけであることを理解すべきである。したがって、本発明の例証的実施形態で多くの修正が行われてもよく、本発明の趣旨及び範囲から逸脱することなく他の配列が考え出されてもよいことを理解すべきである。例えば、本発明の好ましい実施形態は、波形を特定の解析切片へ分割し、中心値に専念したが、他の分割及び中心外領域が、同じように適用されてもよい。

１０…呼吸装置、１２…マスク、１４…送風機、１６…可撓管、１８…モータ、２０…サーボコントローラ、２２…空気出口、２４…制御回路、２６…排気口、２８…流量限定要素、３０…管、３２…マスク圧力トランスデューサ、３３…管、３４…微分圧力トランスデューサ、３６…増幅器、３８…増幅器、４０…フィルタ、４２…変換器、４４…マイクロプロセッサ。

Claims (28)

1. 吸息波パターンにおける流量制限された形の存在をプロセッサ内で決定する方法であって、
   前記吸息波パターンを、互いに直交する２つの基底ベクトルへ前記プロセッサが分解するステップと、
   前記吸息波パターンに近似的に当て嵌まるように、前記基底ベクトルの関数としてのデータ曲線を前記プロセッサが構築するステップと、
   前記関数から、前記基底ベクトルの１つにそれぞれ関連付けられた２つの因子を前記プロセッサが決定するステップと、
   前記因子の関数に基づいて前記流量制限された形の存在を前記プロセッサが特定するステップと、
   を含み、
   それにより、Ｍ形を示す前記流量制限された形が検出され得る、方法。
2. 各基底ベクトルが、前記吸息波パターンの時間区間を正規化して、吸息がゼロで始まって所定ラジアンで終わるようにする、請求項１に記載の方法。
3. 前記基底ベクトルが正弦関数である、請求項１に記載の方法。
4. 前記基底ベクトルが、前記所定ラジアンの正弦関数である、請求項２に記載の方法。
5. 前記基底ベクトルに関連付けられた因子を空間の第１の位置に前記プロセッサがマップするステップを更に含む、請求項１に記載の方法。
6. 前記空間が、Ｍ形息を表す第２の位置と、非Ｍ形息を表す第３の位置とを含む、請求項５に記載の方法。
7. 前記第１の位置と前記Ｍ形息を表す第２の位置との間の距離の測定値を前記プロセッサが計算するステップを更に含む、請求項６に記載の方法。
8. 流量制限された息の存在をプロセッサ内で決定する方法であって、
   所定数の点の上で吸息を前記プロセッサが補間するステップと、
   第１の基礎関数と前記補間された吸息との関数である第１の因子を前記プロセッサが計算するステップと、
   第２の基礎関数と前記補間された吸息との関数である第２の因子を前記プロセッサが計算するステップと、
   前記第１及び第２の因子に基づいて、前記流量制限された息の存在を前記プロセッサが決定するステップと、
   を含み、
   Ｍ形を示す前記流量制限された息が検出され得る、方法。
9. 前記第１及び第２の基礎関数が正弦関数である、請求項８に記載の方法。
10. 前記第１の基礎関数がｓｉｎ（πｔ）であり、前記第２の基礎関数がｓｉｎ（３πｔ）である、請求項８に記載の方法。
11. 後続計算のために前記基礎関数を前記プロセッサが記憶するステップを更に含む、請求項８に記載の方法。
12. 前記第１及び第２の因子に基づいて形状因子を前記プロセッサが決定するステップを更に含む、請求項８に記載の方法。
13. 前記形状因子が次式
    によって算出され、
    Ｆ_１及びＦ_２がそれぞれ前記第１の因子及び前記第２の因子を表す、請求項１２に記載の方法。
14. 前記決定するステップが、前記形状因子を、非流量制限息を示す閾値と比較するステップを含む、請求項１２に記載の方法。
15. 呼吸を感知し、吸息波パターンを生成する流量検出器と、
    前記流量検出器に接続されたコントローラと、
    を備え、
    前記コントローラが、
    前記吸息波パターンを、互いに直交する２つの基底ベクトルへ分解し、
    前記吸息波パターンに近似的に当て嵌まるように、前記基底ベクトルの関数としての曲線を構築し、
    前記関数から、前記基底ベクトルの１つにそれぞれ関連付けられた２つの因子を決定し、
    前記因子の関数に基づいて、前記吸息波パターンにおいて流量制限された形の存在を特定し、
    それにより、Ｍ形を示す前記流量制限された形が検出され得るように構成される、呼吸装置。
16. 各基底ベクトルが、前記吸息波パターンの時間区間を正規化して、吸息がゼロで始まって所定ラジアンで終わるようにする、請求項１５に記載の装置。
17. 前記基底ベクトルが正弦関数である、請求項１５に記載の装置。
18. 前記基底ベクトルが、前記所定ラジアンの正弦関数である、請求項１６に記載の装置。
19. 前記コントローラが、前記基底ベクトルに関連付けられた因子を空間の第１の位置にマップするように構成される、請求項１５に記載の装置。
20. 前記空間が、Ｍ形息を表す第２の位置と、非Ｍ形息を表す第３の位置とを含む、請求項１９に記載の装置。
21. 前記コントローラが、前記第１の位置と前記Ｍ形息を表す第２の位置との間の距離の測定値を計算するように構成される、請求項２０に記載の装置。
22. 息を感知する流量検出器と、
    前記流量検出器に接続されたコントローラと、
    を備え、
    前記コントローラが、
    所定数の点の上で吸息を補間し、
    第１の基礎関数と前記補間された吸息との関数である第１の因子を計算し、
    第２の基礎関数と前記補間された吸息との関数である第２の因子を計算し、
    前記第１及び第２の因子に基づいて、前記息が流量制限されているかどうかを特定し、
    Ｍ形を示す前記流量制限された息が検出され得るように構成される、呼吸装置。
23. 前記基礎関数が正弦関数である、請求項２２に記載の装置。
24. 前記第１の基礎関数がｓｉｎ（πｔ）であり、前記第２の基礎関数がｓｉｎ（３πｔ）である、請求項２２に記載の装置。
25. 前記コントローラが、後続計算で使用されるように前記基礎関数を記憶するように構成される、請求項２２に記載の装置。
26. 前記コントローラが、前記第１及び第２の因子に基づいて形状因子を導出するように構成される、請求項２２に記載の装置。
27. 前記形状因子が次式
    によって算出され、
    Ｆ_１及びＦ_２がそれぞれ前記第１の因子及び前記第２の因子を表す、請求項２６に記載の装置。
28. 前記コントローラが、前記形状因子を、非流量制限息を示す閾値と比較することによって、前記息が流量制限されているかどうかを決定するように構成される、請求項２６に記載の装置。