複素アドミッタンス値による、呼吸気道の閉塞性および開放性無呼吸の識別
【課題】閉塞性または開放性の無呼吸の発生を確定するための方法および器具が、開示される。
【解決手段】空気流量信号を得るために、患者からの呼吸気流量が測定される。患者の気道に既知の周波数の振動圧力波形を印加し、患者のアドミッタンスを表す複素量を計算し（１２）、その値を開放性または閉塞性の無呼吸を示す範囲（１４、１６）と比較することにより、無呼吸の確定が行われる。該方法は、アドミッタンスを計算するために使用されるモデルが、呼吸器の詳細に基づいていない場合でも、開放性の無呼吸を閉塞性の無呼吸と識別する。さらに、気道を特徴付ける必要性を回避するために、患者のアドミッタンスは、正常呼吸時のアドミッタンスと比較され得る。
【解決手段】空気流量信号を得るために、患者からの呼吸気流量が測定される。患者の気道に既知の周波数の振動圧力波形を印加し、患者のアドミッタンスを表す複素量を計算し（１２）、その値を開放性または閉塞性の無呼吸を示す範囲（１４、１６）と比較することにより、無呼吸の確定が行われる。該方法は、アドミッタンスを計算するために使用されるモデルが、呼吸器の詳細に基づいていない場合でも、開放性の無呼吸を閉塞性の無呼吸と識別する。さらに、気道を特徴付ける必要性を回避するために、患者のアドミッタンスは、正常呼吸時のアドミッタンスと比較され得る。
【発明の詳細な説明】
【関連出願】
【０００１】
本願は、２００６年８月３０日出願の米国仮出願第６０／８２３，９７３号、および２００７年５月４日出願の米国仮出願第６０／９１６，１４７号の優先権を主張し、それらの明細書および図面は、参照することにより本明細書に組み込まれる。
【技術分野】
【０００２】
本発明は、開放性および閉塞性の無呼吸の識別（すなわち、呼吸の完全停止）に関する。開放性の無呼吸において、気道は開存している一方で、閉塞性の無呼吸においては、気道は完全に閉塞している。そのような無呼吸の識別は、ヒトの健康に悪影響を及ぼす呼吸状態の診断および治療に有利である。
【背景技術】
【０００３】
本明細書に使用される場合、「気道」という表現は、気管を含む鼻孔と気管支との間の呼吸器系の解剖学的部位として理解されるものとする。「呼吸」という表現は、吸気（吸い込み）、続いて呼気（吐き出し）の連続的繰り返し事象として理解されるものとする。
【０００４】
睡眠時無呼吸症候群において、ヒトは睡眠中に呼吸を停止する。１０秒を超える空気流の停止は、「無呼吸」と呼ばれる。無呼吸は、血液酸素化の減少を引き起こし、したがって睡眠障害を引き起こす。無呼吸は、従来、呼吸努力がない中枢性、または呼吸努力がある閉塞性のいずれかに分類される。ある種の中枢性無呼吸では気道は開存しており、対象は、単に呼吸をしようとしないだけである。逆に、他の中枢性無呼吸およびすべての閉塞性無呼吸では、気道は開存していない（すなわち、閉塞している）。閉塞は通常、舌または軟口蓋のレベルで存在する。気道もまた、部分的に閉塞している（すなわち、狭窄または部分的に開存している）。これもまた、換気の減少（呼吸低下）、血液酸素化の低下、および睡眠障害を引き起こす。
【０００５】
睡眠中の閉塞性呼吸の危険性は、閉塞性睡眠時無呼吸（ＯＳＡ）症候群との関連でよく知られている。無呼吸、呼吸低下、および重度のいびきは、睡眠障害の原因、および、ある種類の心臓病の危険因子として認識される。いびきまたは睡眠時無呼吸を伴わない、睡眠中の上気道抵抗の増加（上気道抵抗症候群）もまた、睡眠断片化および日中の眠気を引き起こす可能性がある。
【０００６】
これらの症候群の一般的治療形態は、持続的気道陽圧法（ＣＰＡＰ）の実行である。簡単に述べると、ＣＰＡＰ治療は、通常４〜２０ｃｍＨ２Ｏの範囲の陽圧の提供により、気道の空気の副木としての機能を果たす。出口が、空気送達ホースを介して、患者の顔に密閉して係合される、鼻（または鼻および／もしくは口腔）マスクまで通る、電動ブロワまたは他のフロージェネレータ（ＦＧ）によって、空気が気道に供給される。排気ポートは、マスク近接の送達チューブ内に設けられる。バイレベルＣＰＡＰおよびオートＣＰＡＰ等のより高度な形態のＣＰＡＰは、米国特許第５，１４８，８０２号、および第５，２４５，９９５号にそれぞれ記載されている。
【０００７】
前述のように、中枢性無呼吸は、必ずしも気道の閉塞を伴うものではなく、非常に浅い睡眠中に、かつ上気道の状態とは無関係の様々な心臓、脳血管、および内分泌疾患を有する患者において発生する場合が多い。無呼吸が気道の閉塞を伴わずに発生するそれらの場合において、ＣＰＡＰ等の技術による疾患の治療における利点はほとんどない。自動ＣＰＡＰシステムにおいて、ＣＰＡＰ空気圧副木を不適切に増加させることを回避するために、開放性気道を有する無呼吸を、閉塞性気道を有する無呼吸と正確に識別することが重要である。そのような不必要の圧力増加は、呼吸を反射的に阻害し、さらに、呼吸障害を悪化させる。
【０００８】
米国特許第５，２４５，９９５号は、睡眠中の吸気および呼気の圧力測定により、いびきおよび異常呼吸のパターンの検出が可能であり、それにより、閉塞の前症状の出現または他の型の呼吸障害の早期指摘をもたらす方法について記載している。特に、理想的には、閉塞症状の出現または他の型の呼吸障害の発生を防止するために、呼吸パラメータのパターンが監視され、所定パターンの検出によって、ＣＰＡＰ圧力が上昇し、気道内圧の増加を提供する。
【０００９】
強制振動法の従来における使用
米国特許第５，７０４，３４５号、発明の名称「Ｄｅｔｅｃｔｉｏｎ Ｏｆ Ａｐｎｅａ Ａｎｄ Ｏｂｓｔｒｕｃｔｉｏｎ Ｏｆ Ｔｈｅ Ａｉｒｗａｙ Ｉｎ Ｔｈｅ Ｒｅｓｐｉｒａｔｏｒｙ Ｓｙｓｔｅｍ」は、既知の周波数の振動圧力波形が、患者の気道に印加され、該振動圧力波形により誘導される既知の周波数における空気流量信号成分の振幅が計算され、閾値と比較される強制振動法（ＦＯＴ）による、気道開存性の測定を含む、閉塞性呼吸を示す異常呼吸パターンを検知および検出するための様々な技術について記載している。本発明は、第５，７０４，３４５号の特許に開示される方法および器具の改良形態である。
【００１０】
気道開存性を確定するための約４Ｈｚの周波数での圧力振動の使用は、ＲｅｓＭｅｄ ＡｕｔｏＳｅｔ Ｃｌｉｎｉｃａｌ自動ＣＰＡＰ装置およびＰＩＩ Ｐｌｕｓで用いられた。ＦＯＴを使用したこれらの機械において、圧力は、マスクにおいて測定され、流量は、マスクの通気孔の患者側で、マスクに非常に近接して測定された。本発明は、少なくとも強制振動法の効果を分析するために、フロージェネレータで、またはフロージェネレータの近傍での圧力および流量の測定に利点を見出す。
【００１１】
フロージェネレータと患者との間の空気通路に、中程度の漏出および中程度の「抵抗」が存在する場合、ＦＯＴの従来技術の実施は、閉塞性および開放性の無呼吸の識別においてそれほど正確ではない。例えば、３〜４ｃｍＨ２Ｏ／（ｌ／ｓ）抵抗から成る受動的患者シミュレーションは、調整可能な漏出により、１５ｌ／ｍｉｎの漏出における開放性気道、および２０ｌ／ｍｉｎの漏出における閉塞性気道を示すであろう。求められているものは、開放性および閉塞性の無呼吸をより正確に識別するシステムである。具体的には、求められているものは、単なる非線形抵抗としての、空気通路の成分を扱うだけではなく、空気通路のインピーダンスの容量分および誘導分を考慮するアルゴリズムを利用するシステムである。
【発明の概要】
【００１２】
２００６年８月３０日出願の、先出願された米国仮特許出願第６０／８２３，９７３号において、呼吸気道の閉塞性および開放性無呼吸を識別するための効果的な方法は、患者の気道の複素アドミッタンスを確定し、複素アドミッタンスの絶対値を閾値と比較することによる、と開示された。その開示は、患者の圧力および流量を確定するために、患者回路を考慮し、各成分（患者回路、換気流、および漏出）をモデリングするアルゴリズムを用いた。具体的には、小信号ホース圧力降下が、２ポートネットワークとして、モデリングされ、そこでは、パラメータが平均流量の関数であった。本発明は、開放性および閉塞性の無呼吸の識別において、その複素量からのさらなる情報、すなわち、アドミッタンスの位相角、または複素量の実成分および虚数成分の両方を含むことによる、その技術の改良形態である実施形態を含む。具体的には、意外にも、アドミッタンスの位相角が考慮されると、データの分離が著しく優れることが、開放性および閉塞性の無呼吸で、すでに識別されたデータを分析することによって見出された。さらに、流動範囲にわたって患者回路のＡＣインピーダンスを測定する必要なく、圧力降下をモデリングすることにより、これまでのシステムが単純化される。
【００１３】
アドミッタンスの位相に対するアドミッタンスの大きさのグラフは、開放性または閉塞性気道を示すグラフのはっきりと区別できる領域を提供する。あるいは、アドミッタンスの実部は、アルガン図内の虚部に対して示され得て、この場合もやはり、開放性および閉塞性の無呼吸に関連している、図のはっきりと区別できる領域をもたらす。
【００１４】
本発明は、無呼吸になりやすい患者を治療するための、改善された方法および器具を提供する。具体的には、本発明は、それぞれ開放性または閉塞性の無呼吸の特徴を示す、患者の複素アドミッタンスの値の組を識別することにより、無呼吸が閉塞性であるか、または開放性であるかを確定するための改善された方法に依存する。アドミッタンスは、例えば、流量インピーダンスの２ポートモデルとしてパラメータ化されるか、または実験的に導かれた抵抗定数および理論的に決定された流量のイナータンス定数を用いる理論モデルから、抵抗および慣性成分を確定する、器具の複素インピーダンスの逆数である。
【００１５】
正弦波（例えば、４Ｈｚ）圧力振動が、空気通路への入口で印加される一方で、流量および圧力は、空気通路の入口および出口の両方で測定される。２ポートモデルまたはより理論的なホース圧力降下モデルに基づき、アドミッタンスは、患者の空気流のＡＣ成分（換気流および漏出のＡＣ成分を引くことにより得られる）、およびマスク圧力のＡＣ成分から確定される。また、換気流は、次いで、２ポートモデルの場合、マスク圧力への改善された二次関係から、またはより理論的なホース圧力降下モデルの場合、線形化計算から確定される。漏出は、推定される漏れ係数から確定される。これらの量のＡＣ成分の計算は、方形波で正規直交関数の組を近似するのではなく、入力振動数でフーリエ正弦および余弦成分を使用することにより、従来の推定よりも改善される。
【００１６】
本発明は、患者の気道の開存性を確定するための方法を開示し、方法は、
患者の気道に既知の周波数の振動圧力波形を印加するステップと、
フロージェネレータにおいて、呼吸気流量および圧力を測定するステップと、
２ポートインピーダンスモデルから、マスクで流量および圧力のＡＣ値を計算するステップと、
複素アドミッタンスが、開存性の特徴を示す領域内にあるかどうか確定することにより、気道が開存していることを確定するステップと、
を含む。
【００１７】
アドミッタンスは、患者の流量とマスク圧力とのＡＣ値の比から有利に確定され、また、複素アドミッタンスの値を、気道が開存していることが明らかな値の範囲と比較するステップがある。
【００１８】
本発明は、患者の気道の開存性の程度を確定するための方法をさらに開示し、方法は、
患者の気道への入口で、既知の周波数、大きさ、および位相の振動圧力波形を印加するステップと、
患者からの呼吸気流量を測定するステップと、
該振動圧力波形により誘導される、該既知の周波数における該空気流の成分の大きさおよび位相を確定するステップと、
開存性を示す領域内の複素アドミッタンスの位置として、開存性の程度を確定するステップと、
を含む。
【００１９】
上記の方法において、所定の波形を供給しない送達システムの影響を最小限に抑えるために、実際に生成される圧力波形が測定される。この技術において、最初に、フロージェネレータにおいて特定の振幅の波形を印加し、マスクで圧力波形の振幅を計算または観察し、次に、マスクで望ましい振幅を生成するために、フロージェネレータで振幅を調整する（典型的は増加させる）。システムが、これらの小信号に対してほぼ直線であると仮定すると、駆動波形をどれほど増加させるかに関する計算は、マスクでの実際の圧力の大きさに対する望ましい圧力の大きさの比の計算である。この方法の利点は、信号対ノイズ比の改善である。
【００２０】
本発明は、大気圧を上回るように上昇させた選択可能な圧力で、継続的に患者の気道に呼吸可能な空気を制御可能に供給することにより、患者の気道へのＣＰＡＰ治療の実行を制御するための方法をさらに開示し、方法は、
時間の関数として患者からの呼吸気流量を測定するステップと、
開存性の領域の中心から、該アドミッタンスの偏差を確定するステップと、
により、無呼吸が発生していると確定される場合、ＣＰＡＰ治療圧を開始するか、または増加させるステップと、
を含む。
【００２１】
本発明は、患者の気道の開存性を決定するための器具をさらに開示し、器具は、
患者の気道に既知の周波数の振動圧力波形を印加するための手段と、
患者からの呼吸気流量を測定するための手段と、
アドミッタンスが、該振動圧力波形により誘導される該既知の周波数で、開存性領域内にある場合、気道が開存していることを確定するための手段と、
を備える。
【００２２】
本発明は、患者の気道の開存性の程度を決定するための器具をさらに開示し、器具は、
患者の気道に、既知の周波数および振幅の振動圧力波形を印加するための手段と、
患者からの呼吸気流量の複素アドミッタンスを確定するための手段と、
該誘導されたアドミッタンスの、開存性領域の中心からの偏差として、開存性の程度を確定するための手段と、
を備える。
【００２３】
本発明は、患者の気道の開放性および閉塞性無呼吸を識別する方法をさらに開示し、
（ｉ）空気送達チューブと患者との接触面を介して、患者に呼吸装置を接続するステップと、
（ｉｉ）空気送達チューブを介して呼吸装置から患者に、患者への交流圧力波形を送達するステップと、
（ｉｉｉ）呼吸装置における空気の流量および圧力を測定するステップと、
（ｉｖ）空気送達チューブのインピーダンスの容量分を確定するステップと、
（ｖ）該測定された空気の流量と圧力ならびに該容量分から、患者のアドミッタンスを確定するステップと、
（ｖｉ）該患者のアドミッタンスに基づき、気道の開放性および閉塞性無呼吸を識別するステップと、
を含む。
【００２４】
本発明は、患者の気道の開放性および閉塞性無呼吸を識別する方法をさらに開示し、
（ｉ）空気送達チューブと患者との接触面を介して、患者に呼吸装置を接続するステップと、
（ｉｉ）空気送達チューブを介して呼吸装置から患者に、患者への交流圧力波形を送達するステップと、
（ｉｉｉ）呼吸装置における空気の流量および圧力を測定するステップと、
（ｉｖ）空気送達チューブのインピーダンスの誘導分を確定するステップと、
（ｖ）該測定された空気の流量と圧力ならびに該容量分から、患者のアドミッタンスを確定するステップと、
（ｖｉ）該患者のアドミッタンスに基づき、気道の開放性および閉塞性無呼吸を識別するステップと、
を含む。
【００２５】
本発明の別の形態において、気道の開存性は、未定義、または「未決定」領域内にあることが確定される。本発明の別の態様は、開存性は、開放性から閉塞性への段階的な尺度を有すると定義されることである。別の態様は、無呼吸の識別領域は、複素平面内に画定され、領域は必ずしも隣接していないことである。さらなる改良点は、識別領域は、平均漏出レベルの関数であり得て、具体的には、「未決定」領域は、恐らく漏出の増加とともにその大きさが増大することである。
【００２６】
本発明のさらなる形態は、特許請求の範囲に示される通りである。
【００２７】
本発明の実施形態が、ここで添付図面を参照して説明される。
【図面の簡単な説明】
【００２８】
【図１】実施形態の基本的技法のフロー図を示す。
【図２】本発明を具体化する器具を図式で示す。
【図３】本発明の２ポート解析の電子的な類似形態を示す。
【図４】詳細なアルゴリズムに基づいた、複素アドミッタンスの値を示す。
【図５】単純化アルゴリズムに基づいた、複素アドミッタンスの値を示す。
【発明を実施するための形態】
【００２９】
本発明を実施するシステム
本発明は、全体として参照することによって組み込まれる、米国特許第５，７０４，３４５号に開示される実施形態の改良形態である。図１は、本発明の一実施形態の基本的技法のフロー図である。本発明の第１のステップ１０は、無呼吸が発生している間の、フロージェネレータ近傍の点での、経時的な呼吸流量および圧力の測定である。この情報は、その後の処理のための定性的測定値として使用されるアドミッタンス値を生成するために、ステップ１２において処理される。ステップ１４〜１６は、時間窓内の複素アドミッタンス値を開存性領域と比較することにより、閉塞性、開放性、または混合性無呼吸が発生しているかどうかを検出する。
【００３０】
無呼吸が進行中である場合、次に、無呼吸が開放性であるか、または閉塞性であるかの確定が続く。開放性気道を有する無呼吸が発生している場合、かつ必要に応じて、事象は、ステップ１８で記録される。ステップ１６の結果が、閉塞性気道を有する無呼吸が発生しているとする場合、ＣＰＡＰ治療圧の増加が、ステップ２０で行われる。必要に応じて、ステップ２０は、検出された異常を任意に記録することを含み得る。
【００３１】
開放性気道を有する無呼吸の場合、無呼吸の発生を未然に防ぐか、または少なくとも軽減するために必要とされる、最小圧力を設定しようとする一般的方法論に従い、ＣＰＡＰ治療圧は減少される。ステップ２２における減少量は、必要に応じて、ゼロであり得る。
【００３２】
図１に表される技法は、患者のＣＰＡＰ圧力が、適宜に経時的に制御される、臨床的実施形態の技法である。純診断的実施形態は、それぞれ、ステップ２０およびステップ２２における、ＣＰＡＰ圧力増加および圧力減少の動作を省略する以外は、同一の方法で機能する。
【００３３】
図２は、図１の技法を実施するための一実施形態に従った、臨床ＣＰＡＰ器具を図式で示す。マスク３０は、鼻マスク、および／または顔マスクのいずれであっても、患者の顔に密閉して係合される。新鮮な空気または酸素富化空気は、可撓性管３２を介してマスク３０に入り、可撓性管３２は、次いで、電動タービン（フロージェネレータ）３４に接続され、そこに、空気入口３６が提供される。ＣＰＡＰ治療として、マスク３０に供給される空気の圧力を増加あるいは減少させるために、タービン用のモータ３８は、モータサーボユニット４０により制御される。マスク３０はまた、管３４とマスク３０との接合部に近接する排気ポート４２を含む。
【００３４】
流動抵抗要素４４は、フロージェネレータ３４に隣接する。これは、差圧が存在するアイリス絞り形状を取ることが可能である。流動抵抗要素４４のマスク側は、小口径管４６を介して、圧力トランデューサ４８と、差圧トランデューサ５０の入力に接続される。流動抵抗要素４４の反対側の圧力は、もう１つの小口径管５２を介して、差圧トランデューサ５０の他の入力に移動する。
【００３５】
圧力トランデューサ４８は、流動圧力に比例した電気信号を生成し、それは、増幅器５６により増幅され、マルチプレクサ／ＡＤＣユニット５８とモータサーボユニット４０との両方に移動する。モータサーボユニット４０に提供される信号は、静圧が、定値圧力に近似するように制御されることを確実にするための、フィードバックとしての機能を果たす。
【００３６】
流動抵抗要素４４にわたって検知される差圧は、差圧トランデューサ５０から電気信号として出力され、もう１つの増幅器６０により増幅される。したがって、増幅器５６からの出力信号は、マスクまたは呼吸の空気流量の測定値を表す。制御装置６２は、多くの処理機能を実行するようにプログラムされる。
【００３７】
圧力および流量は、定常値およびＡＣ値から成ると考えられ得て、後者は、４Ｈｚの周波数を有する圧力への強制振動信号の影響を示す。以下において、「定常」量は、（ａ）例えば、フィルタリング操作により、振動成分が意図的に除去されている量、または（ｂ）振動成分が除去されている量に部分的または完全に基づいた、計算結果、または（ｃ）フロージェネレータで測定される圧力および流量等の瞬間量（したがって、振動成分を含む）に基づき計算される量のいずれか、およびシステムの無効分を部分的あるいは完全に無視し、それを、単に（場合により非線形の）抵抗システムとして扱う、空気通路および患者の漏出のモデルである。管３２に沿った定常圧力損失は、チューブを通る流量、および例えばテーブル索引による、管の静圧流量特性の知識から計算される。次に、マスクにおける定常圧力は、チューブの圧力損失を引くことにより計算される。次に、チューブ３２に沿った圧力損失は、圧力ジェネレータ３４で望ましい瞬間圧力を得るために、マスクにおける望ましい設定圧力に加えられる。排気口４２を通る流量は、例えばテーブル索引による、排気口の圧力流量特性からの、マスクにおける圧力から計算される。定常マスク流量は、管３２を通る流量から、排気口４２を通る流量を引くことにより計算される。次に、定常患者流量は、例えば、入力が瞬間マスク流量である、１０秒の時定数を有する一次ローパスフィルタにより、ＣＰＡＰ装置において確定され得る推定定常漏出を、定常マスク流量から引くことにより計算される。
【００３８】
制御装置６２により実行される技法を、ここで図２の器具を参照して説明する。患者呼吸流量が、非常に低い、またはゼロである場合（漏出があれば、患者が無呼吸である場合、マスクの流動は停止しないことに留意する）、外部から誘導された振動技術を使用することにより、気道開存性の確定（ステップ１４〜１６）が行われる。気道が開放性であるが、呼吸筋が弛緩している場合（すなわち、開放性気道を有する中枢性無呼吸）、外部に起因するマスク圧力の小さい変動が、肺を膨張および収縮させることにより、かつ肺内のガスを圧縮および減圧することにより、小さい呼吸気流を誘導する。反対に、気道が閉塞性である場合、空気流は誘導されない。これは、以下のように定量化される。
【００３９】
アドミッタンスの閾値による無呼吸の識別
アドミッタンスＹは、
Ｙ＝Ｇ＋ｉＢ
により与えられ、Ｇは、コンダクタンスであり、Ｂは、サセプタンスである。
【００４０】
気道が開放性であるかどうかを判断するために、患者のアドミッタンスＹの複素数値は、値の領域と比較される。この閾値の値は、以下の考察に基づき選択され得る。
【００４１】
マグニチュードよりもよりよい分類指標である角度に関する説明は、回路内の複素インピーダンスの相互作用に関連付けられ得る。漏出および換気流は、本質的に複素（ｊ）成分がない非線形抵抗である。患者は、基本的に接地された、直列の抵抗およびコンプライアンス（キャパシタンス）である。ホース圧力降下もまた、流動のイナータンスの形で複素成分を有する。したがって、インピーダンス全体の大きさは、すべての成分（ホース、漏出、換気流、患者の気道）が実部を有するため、それらに影響される。しかしながら、虚数成分は、慣性ホース圧力降下と容量肺との間の相互作用である。それは、慣性と容量性リアクタンスとの相対的な「強度」を比較している。
【００４２】
また、角度および振幅ではなく、インピーダンスの実部および虚部を使用することも可能である。これは、ＣＰＵ負荷の減少という利点を有する。振幅の計算は、平方根を必要とし、角度の計算は、逆タンジェントを伴う。（これらの両方は、計算コストの低い、ルックアップテーブルを伴う技術により、必要とされる精度で容易に行われ得る。）実および虚の両方、または振幅および角度の両方を使用することにより、よりロバストな分類が得られる。開放性および閉塞性の値のグループ化は、はるかに明らかとなる。
【００４３】
図４の値から確定される通り、可能な閾値の組は、
角度＝ＡＲＧ（患者＿アドミッタンス）
ｉｆ（角度＜−２ ＯＲ 角度＞１．２）
気道状態＝閉塞性；
ｅｌｓｅ ｉｆ（角度＜０．９ ＡＮＤ 角度＞−１．５）
軌道状態＝開放性；
ｅｌｓｅ
気道状態＝未定義；
ｅｎｄ ｉｆ
となり、角度は、ラジアン単位である。閾値に対する許容値は、図４から確定される開放性および閉塞性の無呼吸を表す値の区別を可能にする任意の値であり得る。
【００４４】
４Ｈｚにおける患者抵抗も、気道の状態を示す。この抵抗は、上記の方程式におけるコンダクタンスの逆数ではない。４Ｈｚにおける患者のインピーダンスは、
Ｚ＝Ｒ＋ｉＸ
であり、Ｙ＝１／Ｚである。患者抵抗Ｒは、
Ｒ＝Ｒｅ（Ｚ）＝Ｒｅ（１／Ｙ）
により与えられる。
【００４５】
代替実施形態では、Ｒが、｜Ｙ｜よりも、気道の状態を特徴付けるために使用され得る。
【００４６】
アドミッタンスの計算
患者のアドミッタンスは、方程式
【数１】
により計算され、分子は、ＡＣ患者流量のマグニチュードであり、実質的に患者流量の（時間に関する）微分の測定値である。分母は、マスク圧力の同一の微分である。
【００４７】
換気流および漏出の補正
患者の空気流量は、マスクへの流入量から、換気口を通る流量と漏れ流量を引くことにより確定される。
【数２】
【００４８】
平均換気流量は、測定中の平均マスク圧力、またはより好ましくは、測定中のマスク圧力の平方根の平均の平方（漏出の計算を行うためにも使用され得る）のいずれかから、確定される。マスク圧力自体は、空気通路のＡＣ動作を無視して、従来の方法で計算される。
【００４９】
ＡＣ換気流は、動作点の線形近似を使用して計算される。
【数３】
【００５０】
マスクにおけるＡＣ流量を確定するための、マスク圧力のモデリング
マスク圧力は、換気流の関数として換気口を通じた圧力降下を得る方程式において、係数ｋ１およびｋ２を求めることにより確定される。
【数４】
したがって、特定のマスク圧力における換気流を計算するために、そのマスク圧力における換気流に対する二次方程式を解く。
【００５１】
【数５】
は、先の方程式の微分により得られ、
【数６】
を与える。
【００５２】
【数７】
であるため、小信号の近似は、
【数８】
を得て、そこから、
【数９】
が、妥当な近似として得られる。
【００５３】
漏れ流量の計算
患者の空気流の計算を完全にするために、マスクへの流入からの漏れ流量を計算することも必要である。実際の漏出は、患者の４Ｈｚのアドミッタンスが計算される間、計算される。典型的な漏れ係数推定は、少し前に起こったことに関連し、従来技術のＡｕｔｏＳｅｔ ＣＳおよびＡｕｔｏＶＰＡＰ装置におけるように、
【数１０】
から漏れ係数Ｋｌｅａｋを推定し得て、上線は、測定期間にわたる平均値を示す。無呼吸中、漏出は、全非換気流量に等しく、
【数１１】
呼吸中、漏出は、平均非換気流に（よい近似で）等しい。上記のＫｌｅａｋの推定において、すべての量は、瞬間値に基づき、空気通路のＡＣ特性を参照することなく計算され、空気通路を非線形抵抗器として扱う。
【００５４】
ＤＣ（Ｋｌｅａｋに対する上記の式に基づく）とＡＣとの両方の、瞬間漏出に対するモデルは、
【数１２】
であり、換気流と同様に、漏れ流量の（４Ｈｚにおける）ＡＣ成分は、小信号近似
【数１３】
を使用して得られ、上記の方程式を直接微分することにより、
【数１４】
を計算し得て、
【数１５】
を得る。
【００５５】
【数１６】
の平均は、
【数１７】
または
【数１８】
のいずれかであり得る。後者が好ましく（Ｐｍａｓｋが定数ではない場合、それが、
【数１９】
の平均値の不偏推定値を得るため、かつＫｌｅａｋを推定するために、それをすでに計算したため）、
【数２０】
を得る。
【００５６】
ＡＣ値の確定
この説明を通して、一貫して参照されるＡＣ値は、測定から、励振周波数における三角フーリエ成分として確定される。これは、正規直交関数の組として方形波を使用していた強制振動法の、従来技術の使用の改良形態である。
【００５７】
すべてのＡＣ量は、複素数として計算される。具体的には、マスクにおけるＡＣ圧力および流量は、正弦および余弦関数のそれぞれを有する積和を使用して、測定中、瞬間圧力および流量の正弦および余弦成分を得ることにより計算され、それにより、成分係数ｃｓおよびｃｃを得て、次に、ＡＣ値としてｃｓ＋ｉｃｃまたはｃｃ＋ｉｃｓのいずれかを記述する。（これらの２つの式のいずれも使用され、２つの式は、単に位相シフトの点で異なるが、１つの式は、初めから終わりまで一貫して使用される）。あるいは、正弦および余弦の最小二乗法等の、励振周波数における圧力および流量の正弦波成分の振幅および位相を推定する、他の標準的方法が使用され得る。
【００５８】
明確には、
【数２１】
を得、励振周波数は、Ｔ／２Ｂであり、Ｔは、正規化係数であり、積分は、正弦および余弦の離散的な標本値を使用した総和により置き換えられ得る。ｆ（ｔ）は、フーリエ係数が必要とされるいかなる関数でもある。
【００５９】
すべての量は、ｔ＝０として象徴的に示される、６秒スライディングウィンドウ（「アドミッタンスウィンドウ」）上で確定される。アドミッタンス計算は、原則として、アルゴリズムのサンプリング周波数（例えば、５０Ｈｚ）で行われ得るが、臨床目的では、これは必須ではなく、計算は、毎４Ｈｚサイクルの終わりに、効果的に実行され得る。あるいは、２Ｈｚ、あるいは１Ｈｚ周波数においてでも、アドミッタンスの計算は、臨床目的では妥当である。
【００６０】
モータ制御装置の遅延により、圧力波形の第１のサイクル（２５０ｍｓ）は、正弦波ではない場合があり、空気通路内の定常の設定に多少の遅延が出るため、データの最初の２５０ｍｓは、無視され得る。
【００６１】
正弦および余弦値は、起動時に生成されるテーブル内に記憶されるか、またはソースコード内の定数として特定され得る。もちろん、合計のオーバーフローが発生しないと仮定して、乗算は、３２ビットの結果を有する、１６ビット＊１６ビットの固定小数点であり、かつこの方法の数値安定性において問題がない場合、十分な精度が提供される。したがって、計算コストは、比較的小さい。
【００６２】
２ポートネットワーク上のインピーダンスのモデリング
マスクにおけるＡＣ圧力およびマスクに入るＡＣ流量の計算は、一実施形態では、特定の周波数における、システム内の流動の２ポート電気ネットワークへの類似から確定される。
【００６３】
図３から、ルーチン回路解析は、
【数２２】
を得て、もちろん、これらは、マスクにおける圧力、ならびにマスクに入る流量のそれぞれに対応する。明確には、
【数２３】
である。
【００６４】
インピーダンスＺｉｎｏｕｔおよびＺｏｕｔｇｎｄは、平均流量の関数とみなされる。これらのインピーダンスは、多くの平均流量レベルで、特定の空気通路に対して測定される（詳細は以下を参照）。アドミッタンスウィンドウ中の、またＺｉｎｏｕｔおよびＺｏｕｔｇｎｄのそれぞれに対する、フロージェネレータの総流量の平均が、測定されたインピーダンス間を線形補間してその平均流量レベルにおけるインピーダンスを確定するために使用される。
【００６５】
入口および出口のＡＣ流量およびＡＣ圧力は、フーリエ係数の計算に関して上記に記載される標準的な積和方法を使用して、（ノイズを減少させるために）３０秒間にわたって確定され、複素数値を得る。これらの値は、以下により２ポートパラメータを計算するために使用される（図１を参照。その図中の電圧および電流は、以下における圧力および流量のそれぞれに対応する）。
【数２４】
分母が非常に小さいか、またはゼロである場合、インピーダンスは、典型的な空気通路のインピーダンスとの関連で、特定の数値的に非常に大きい値を有するとみなされる。
【００６６】
ホース圧力降下の理論モデル
ホース圧力降下は、マスク内圧とＦＧ圧力との間の圧力差である。それは、マスクの前に配置されるすべての要素にわたる圧力損失である。これは、マフラ、加湿器、ＡＢフィルタ、およびマスク接続ホースを含む。
【００６７】
完全なホース圧力降下は、抵抗ホース圧力降下および慣性圧力降下で構成される。
【数２５】
Ｋ１、Ｋ２＝実験的に導かれた定数
ＫＬ＝理論的に決定された流量のイナータンス定数
ＱＦＧ＝ＦＧ流量
ＫＬ＝ρｌ／Ａ
ρ＝空気密度（１．１９ｋｇ／ｍ３）
ｌ＝チューブ長（２または３ｍ）
Ａ＝断面積（πｄ２／４、ｄ＝０．０１９ｍ）
【００６８】
ＲｅｓＭｅｄ ＡｃｔｉｖａおよびＳｗｉｆｔマスクなどのマスクの場合、主要ホースとマスクとの間に別のチューブがある。これもまた、慣性定数において考慮されなければならない。
【００６９】
平均ホース圧力降下は、抵抗部のみで構成される（上線は平均を示す）。
【数２６】
Ｋ１、Ｋ２＝実験的に導かれた定数
【数２７】
【００７０】
ＡＣホース圧力降下を計算するために、ある種の線形化が必要である。抵抗および慣性成分は分離される。抵抗成分は、動作点に関して、上記の平均ホース圧力降下の式を線形化することにより得られる。
【数２８】
【００７１】
したがって、ＡＣホース圧力降下は、ＡＣ流量に、周囲に微小振動が発生する位置における、ホースの非線形抵抗の勾配を乗じたものである。導出のために、振動は、ホース抵抗の勾配の変化と比較して、小さい振幅の振動であると仮定する。
【００７２】
慣性成分は、ＦＧ流量を導関数で乗じる定数である。標準的な線形回路理論により、インダクタンスＬは、インピーダンスｓＬを有し、周波数ωにおける正弦波信号に対して、ｊωＬである。次に、慣性ホース圧力降下の成分は、
【数２９】
となり、２つの組み合わせは、
【数３０】
を得る。
【００７３】
マスク圧力
マスク圧力は、ＦＧ圧力からホース圧力降下を減じたものである。
Ｐｍａｓｋ＝ＰＦＧ−Ｐｈｏｓｅｄｒｏｐ
Ｐｍａｓｋ，ＡＣ＝ＰＦＧ，ＡＣ−Ｐｈｏｓｅｄｒｏｐ，ＡＣ
【実施例】
【００７４】
図４では、太線が、開放性無呼吸を表し、細線が、閉塞性無呼吸を表すグラフが表される。横軸が、アドミッタンスの絶対マグニチュードを表す一方で、縦軸は、複素アドミッタンスの位相を表す。図面から理解され得るように、１つは起点近傍、他方はｙ軸の最上値にある、２つのより細い曲線領域により表される閉塞性無呼吸と、グラフの中心の太い曲線により表される開放性無呼吸が完全に分離されている。実質的に、データのすべてをｘ軸に投影することにより、アドミッタンスの絶対値のみを見ると、細い（閉塞性無呼吸）曲線が、太い（開放性無呼吸）データと接触し、データの分離は、失われる。
【００７５】
図４は、上記の式に従って計算された患者のアドミッタンスデータの、プロッティング結果を示す。ｘ軸の値は、患者のアドミッタンスの平均絶対値である。ｙ軸の値は、患者のアドミッタンスの位相角である。図４から理解され得るように、閉塞性無呼吸および開放性無呼吸に対する値の完全な分離が見られる。
【００７６】
図５は、先のアドミッタンスの計算の代わりに、複素アドミッタンスに対する値が、モデリングにおけるマスクまたは換気回路の任意の特定構成の詳細から、本質的に独立したモデルから確定される、患者のアドミッタンスのプロットの結果を示す。具体的には、図５のデータは、単純漏出オリフィスモデル、Ｐ＝（ＫＱ）２を使用することにより確定された。開放性および閉塞性の無呼吸は、複素アドミッタンスのグラフにおいて分離した領域をもたらすことが、図５から理解され得る。
【００７７】
適度な（近似された患者のアドミッタンス）アルゴリズム
第３のアルゴリズムは、患者回路を総称的にモデリングし、漏出および換気流を組み合わせることにより、患者のアドミッタンスを近似しようと試みる。
【００７８】
患者回路は、上記の複素アルゴリズムと同様の方法でモデリングされるが、定数は、ＡＢフィルタを有する２ｍホースに対する定数で固定される。
【００７９】
漏出および換気流のモデリングは、非患者流量である１つのパラメータに組み合わされる。これは、
漏出＋換気流
としてモデリングされる。
【００８０】
漏出および換気流は、組み合わされ、以下の方程式を使用する際にモデリングされる。
【００８１】
漏出オリフィスは、
【数３１】
を使用して、ランタイムで近似され得る。
【００８２】
次に、瞬間漏出は、
【数３２】
を使用して計算される。
【００８３】
小信号近似を使用すると、
【数３３】
であると言うことができ、式中
【数３４】
であり、Ｋｌｅａｋに置換して、
【数３５】
を得て、小さい平均Ｐｍａｓｋの変化を扱っているので、
【数３６】
で近似することができ、したがって、ＡＣ漏出は、
【数３７】
となる。
【００８４】
近似された患者流量
次に、患者流量は、
Ｑｐａｔｉｅｎｔ，ＡＣ＝ＱＦＤ，ＡＣ−Ｑｌｅａｋ，ＡＣ
である。
【００８５】
Ｑｌｅａｋは、換気流と組み合わされる漏れ流量である（すべての非患者流量）。
【００８６】
分類
近似アドミッタンスを使用して、気道状態を分類する多くの方法がある。アドミッタンスの角度対マグニチュードをプロットすると、角度のみにおけるよい分離を見出し得る。したがって、単純分類は、２つの閾値を有するものとなるが、純粋に角度の上のみである。しかしながら、角度の計算は、除算および逆タンジェントを伴い、さらなる計算を必要とする。代わりに、使用されるアドミッタンスの実部および虚部を有することが好ましいことが示唆される。
【００８７】
閾値
よい分類精度を提供する最も単純な閾値は、単純な対角直線である。これらは、ｙ＞ｘ＋ｂおよびｙ＜ｘ＋ｃの式である。この場合、アドミッタンスの虚部（ｙ）対実部（ｘ）をプロットすると、閉塞性は、ｙ＞ｘ＋ｂであり、開放性は、ｙ＜ｘ＋ｃである。浮動小数点アルゴリズムに基づき、定数に対する値は、ｂ＝０、ｃ＝−０．０３であった。これは、何の誤分類も与えない。データ領域を分離させる他の任意の曲線が使用され得るが、線形分離は、閾値計算の複雑さを減少させ、その理由から好ましい。
【関連出願】
【０００１】
本願は、２００６年８月３０日出願の米国仮出願第６０／８２３，９７３号、および２００７年５月４日出願の米国仮出願第６０／９１６，１４７号の優先権を主張し、それらの明細書および図面は、参照することにより本明細書に組み込まれる。
【技術分野】
【０００２】
本発明は、開放性および閉塞性の無呼吸の識別（すなわち、呼吸の完全停止）に関する。開放性の無呼吸において、気道は開存している一方で、閉塞性の無呼吸においては、気道は完全に閉塞している。そのような無呼吸の識別は、ヒトの健康に悪影響を及ぼす呼吸状態の診断および治療に有利である。
【背景技術】
【０００３】
本明細書に使用される場合、「気道」という表現は、気管を含む鼻孔と気管支との間の呼吸器系の解剖学的部位として理解されるものとする。「呼吸」という表現は、吸気（吸い込み）、続いて呼気（吐き出し）の連続的繰り返し事象として理解されるものとする。
【０００４】
睡眠時無呼吸症候群において、ヒトは睡眠中に呼吸を停止する。１０秒を超える空気流の停止は、「無呼吸」と呼ばれる。無呼吸は、血液酸素化の減少を引き起こし、したがって睡眠障害を引き起こす。無呼吸は、従来、呼吸努力がない中枢性、または呼吸努力がある閉塞性のいずれかに分類される。ある種の中枢性無呼吸では気道は開存しており、対象は、単に呼吸をしようとしないだけである。逆に、他の中枢性無呼吸およびすべての閉塞性無呼吸では、気道は開存していない（すなわち、閉塞している）。閉塞は通常、舌または軟口蓋のレベルで存在する。気道もまた、部分的に閉塞している（すなわち、狭窄または部分的に開存している）。これもまた、換気の減少（呼吸低下）、血液酸素化の低下、および睡眠障害を引き起こす。
【０００５】
睡眠中の閉塞性呼吸の危険性は、閉塞性睡眠時無呼吸（ＯＳＡ）症候群との関連でよく知られている。無呼吸、呼吸低下、および重度のいびきは、睡眠障害の原因、および、ある種類の心臓病の危険因子として認識される。いびきまたは睡眠時無呼吸を伴わない、睡眠中の上気道抵抗の増加（上気道抵抗症候群）もまた、睡眠断片化および日中の眠気を引き起こす可能性がある。
【０００６】
これらの症候群の一般的治療形態は、持続的気道陽圧法（ＣＰＡＰ）の実行である。簡単に述べると、ＣＰＡＰ治療は、通常４〜２０ｃｍＨ２Ｏの範囲の陽圧の提供により、気道の空気の副木としての機能を果たす。出口が、空気送達ホースを介して、患者の顔に密閉して係合される、鼻（または鼻および／もしくは口腔）マスクまで通る、電動ブロワまたは他のフロージェネレータ（ＦＧ）によって、空気が気道に供給される。排気ポートは、マスク近接の送達チューブ内に設けられる。バイレベルＣＰＡＰおよびオートＣＰＡＰ等のより高度な形態のＣＰＡＰは、米国特許第５，１４８，８０２号、および第５，２４５，９９５号にそれぞれ記載されている。
【０００７】
前述のように、中枢性無呼吸は、必ずしも気道の閉塞を伴うものではなく、非常に浅い睡眠中に、かつ上気道の状態とは無関係の様々な心臓、脳血管、および内分泌疾患を有する患者において発生する場合が多い。無呼吸が気道の閉塞を伴わずに発生するそれらの場合において、ＣＰＡＰ等の技術による疾患の治療における利点はほとんどない。自動ＣＰＡＰシステムにおいて、ＣＰＡＰ空気圧副木を不適切に増加させることを回避するために、開放性気道を有する無呼吸を、閉塞性気道を有する無呼吸と正確に識別することが重要である。そのような不必要の圧力増加は、呼吸を反射的に阻害し、さらに、呼吸障害を悪化させる。
【０００８】
米国特許第５，２４５，９９５号は、睡眠中の吸気および呼気の圧力測定により、いびきおよび異常呼吸のパターンの検出が可能であり、それにより、閉塞の前症状の出現または他の型の呼吸障害の早期指摘をもたらす方法について記載している。特に、理想的には、閉塞症状の出現または他の型の呼吸障害の発生を防止するために、呼吸パラメータのパターンが監視され、所定パターンの検出によって、ＣＰＡＰ圧力が上昇し、気道内圧の増加を提供する。
【０００９】
強制振動法の従来における使用
米国特許第５，７０４，３４５号、発明の名称「Ｄｅｔｅｃｔｉｏｎ Ｏｆ Ａｐｎｅａ Ａｎｄ Ｏｂｓｔｒｕｃｔｉｏｎ Ｏｆ Ｔｈｅ Ａｉｒｗａｙ Ｉｎ Ｔｈｅ Ｒｅｓｐｉｒａｔｏｒｙ Ｓｙｓｔｅｍ」は、既知の周波数の振動圧力波形が、患者の気道に印加され、該振動圧力波形により誘導される既知の周波数における空気流量信号成分の振幅が計算され、閾値と比較される強制振動法（ＦＯＴ）による、気道開存性の測定を含む、閉塞性呼吸を示す異常呼吸パターンを検知および検出するための様々な技術について記載している。本発明は、第５，７０４，３４５号の特許に開示される方法および器具の改良形態である。
【００１０】
気道開存性を確定するための約４Ｈｚの周波数での圧力振動の使用は、ＲｅｓＭｅｄ ＡｕｔｏＳｅｔ Ｃｌｉｎｉｃａｌ自動ＣＰＡＰ装置およびＰＩＩ Ｐｌｕｓで用いられた。ＦＯＴを使用したこれらの機械において、圧力は、マスクにおいて測定され、流量は、マスクの通気孔の患者側で、マスクに非常に近接して測定された。本発明は、少なくとも強制振動法の効果を分析するために、フロージェネレータで、またはフロージェネレータの近傍での圧力および流量の測定に利点を見出す。
【００１１】
フロージェネレータと患者との間の空気通路に、中程度の漏出および中程度の「抵抗」が存在する場合、ＦＯＴの従来技術の実施は、閉塞性および開放性の無呼吸の識別においてそれほど正確ではない。例えば、３〜４ｃｍＨ２Ｏ／（ｌ／ｓ）抵抗から成る受動的患者シミュレーションは、調整可能な漏出により、１５ｌ／ｍｉｎの漏出における開放性気道、および２０ｌ／ｍｉｎの漏出における閉塞性気道を示すであろう。求められているものは、開放性および閉塞性の無呼吸をより正確に識別するシステムである。具体的には、求められているものは、単なる非線形抵抗としての、空気通路の成分を扱うだけではなく、空気通路のインピーダンスの容量分および誘導分を考慮するアルゴリズムを利用するシステムである。
【発明の概要】
【００１２】
２００６年８月３０日出願の、先出願された米国仮特許出願第６０／８２３，９７３号において、呼吸気道の閉塞性および開放性無呼吸を識別するための効果的な方法は、患者の気道の複素アドミッタンスを確定し、複素アドミッタンスの絶対値を閾値と比較することによる、と開示された。その開示は、患者の圧力および流量を確定するために、患者回路を考慮し、各成分（患者回路、換気流、および漏出）をモデリングするアルゴリズムを用いた。具体的には、小信号ホース圧力降下が、２ポートネットワークとして、モデリングされ、そこでは、パラメータが平均流量の関数であった。本発明は、開放性および閉塞性の無呼吸の識別において、その複素量からのさらなる情報、すなわち、アドミッタンスの位相角、または複素量の実成分および虚数成分の両方を含むことによる、その技術の改良形態である実施形態を含む。具体的には、意外にも、アドミッタンスの位相角が考慮されると、データの分離が著しく優れることが、開放性および閉塞性の無呼吸で、すでに識別されたデータを分析することによって見出された。さらに、流動範囲にわたって患者回路のＡＣインピーダンスを測定する必要なく、圧力降下をモデリングすることにより、これまでのシステムが単純化される。
【００１３】
アドミッタンスの位相に対するアドミッタンスの大きさのグラフは、開放性または閉塞性気道を示すグラフのはっきりと区別できる領域を提供する。あるいは、アドミッタンスの実部は、アルガン図内の虚部に対して示され得て、この場合もやはり、開放性および閉塞性の無呼吸に関連している、図のはっきりと区別できる領域をもたらす。
【００１４】
本発明は、無呼吸になりやすい患者を治療するための、改善された方法および器具を提供する。具体的には、本発明は、それぞれ開放性または閉塞性の無呼吸の特徴を示す、患者の複素アドミッタンスの値の組を識別することにより、無呼吸が閉塞性であるか、または開放性であるかを確定するための改善された方法に依存する。アドミッタンスは、例えば、流量インピーダンスの２ポートモデルとしてパラメータ化されるか、または実験的に導かれた抵抗定数および理論的に決定された流量のイナータンス定数を用いる理論モデルから、抵抗および慣性成分を確定する、器具の複素インピーダンスの逆数である。
【００１５】
正弦波（例えば、４Ｈｚ）圧力振動が、空気通路への入口で印加される一方で、流量および圧力は、空気通路の入口および出口の両方で測定される。２ポートモデルまたはより理論的なホース圧力降下モデルに基づき、アドミッタンスは、患者の空気流のＡＣ成分（換気流および漏出のＡＣ成分を引くことにより得られる）、およびマスク圧力のＡＣ成分から確定される。また、換気流は、次いで、２ポートモデルの場合、マスク圧力への改善された二次関係から、またはより理論的なホース圧力降下モデルの場合、線形化計算から確定される。漏出は、推定される漏れ係数から確定される。これらの量のＡＣ成分の計算は、方形波で正規直交関数の組を近似するのではなく、入力振動数でフーリエ正弦および余弦成分を使用することにより、従来の推定よりも改善される。
【００１６】
本発明は、患者の気道の開存性を確定するための方法を開示し、方法は、
患者の気道に既知の周波数の振動圧力波形を印加するステップと、
フロージェネレータにおいて、呼吸気流量および圧力を測定するステップと、
２ポートインピーダンスモデルから、マスクで流量および圧力のＡＣ値を計算するステップと、
複素アドミッタンスが、開存性の特徴を示す領域内にあるかどうか確定することにより、気道が開存していることを確定するステップと、
を含む。
【００１７】
アドミッタンスは、患者の流量とマスク圧力とのＡＣ値の比から有利に確定され、また、複素アドミッタンスの値を、気道が開存していることが明らかな値の範囲と比較するステップがある。
【００１８】
本発明は、患者の気道の開存性の程度を確定するための方法をさらに開示し、方法は、
患者の気道への入口で、既知の周波数、大きさ、および位相の振動圧力波形を印加するステップと、
患者からの呼吸気流量を測定するステップと、
該振動圧力波形により誘導される、該既知の周波数における該空気流の成分の大きさおよび位相を確定するステップと、
開存性を示す領域内の複素アドミッタンスの位置として、開存性の程度を確定するステップと、
を含む。
【００１９】
上記の方法において、所定の波形を供給しない送達システムの影響を最小限に抑えるために、実際に生成される圧力波形が測定される。この技術において、最初に、フロージェネレータにおいて特定の振幅の波形を印加し、マスクで圧力波形の振幅を計算または観察し、次に、マスクで望ましい振幅を生成するために、フロージェネレータで振幅を調整する（典型的は増加させる）。システムが、これらの小信号に対してほぼ直線であると仮定すると、駆動波形をどれほど増加させるかに関する計算は、マスクでの実際の圧力の大きさに対する望ましい圧力の大きさの比の計算である。この方法の利点は、信号対ノイズ比の改善である。
【００２０】
本発明は、大気圧を上回るように上昇させた選択可能な圧力で、継続的に患者の気道に呼吸可能な空気を制御可能に供給することにより、患者の気道へのＣＰＡＰ治療の実行を制御するための方法をさらに開示し、方法は、
時間の関数として患者からの呼吸気流量を測定するステップと、
開存性の領域の中心から、該アドミッタンスの偏差を確定するステップと、
により、無呼吸が発生していると確定される場合、ＣＰＡＰ治療圧を開始するか、または増加させるステップと、
を含む。
【００２１】
本発明は、患者の気道の開存性を決定するための器具をさらに開示し、器具は、
患者の気道に既知の周波数の振動圧力波形を印加するための手段と、
患者からの呼吸気流量を測定するための手段と、
アドミッタンスが、該振動圧力波形により誘導される該既知の周波数で、開存性領域内にある場合、気道が開存していることを確定するための手段と、
を備える。
【００２２】
本発明は、患者の気道の開存性の程度を決定するための器具をさらに開示し、器具は、
患者の気道に、既知の周波数および振幅の振動圧力波形を印加するための手段と、
患者からの呼吸気流量の複素アドミッタンスを確定するための手段と、
該誘導されたアドミッタンスの、開存性領域の中心からの偏差として、開存性の程度を確定するための手段と、
を備える。
【００２３】
本発明は、患者の気道の開放性および閉塞性無呼吸を識別する方法をさらに開示し、
（ｉ）空気送達チューブと患者との接触面を介して、患者に呼吸装置を接続するステップと、
（ｉｉ）空気送達チューブを介して呼吸装置から患者に、患者への交流圧力波形を送達するステップと、
（ｉｉｉ）呼吸装置における空気の流量および圧力を測定するステップと、
（ｉｖ）空気送達チューブのインピーダンスの容量分を確定するステップと、
（ｖ）該測定された空気の流量と圧力ならびに該容量分から、患者のアドミッタンスを確定するステップと、
（ｖｉ）該患者のアドミッタンスに基づき、気道の開放性および閉塞性無呼吸を識別するステップと、
を含む。
【００２４】
本発明は、患者の気道の開放性および閉塞性無呼吸を識別する方法をさらに開示し、
（ｉ）空気送達チューブと患者との接触面を介して、患者に呼吸装置を接続するステップと、
（ｉｉ）空気送達チューブを介して呼吸装置から患者に、患者への交流圧力波形を送達するステップと、
（ｉｉｉ）呼吸装置における空気の流量および圧力を測定するステップと、
（ｉｖ）空気送達チューブのインピーダンスの誘導分を確定するステップと、
（ｖ）該測定された空気の流量と圧力ならびに該容量分から、患者のアドミッタンスを確定するステップと、
（ｖｉ）該患者のアドミッタンスに基づき、気道の開放性および閉塞性無呼吸を識別するステップと、
を含む。
【００２５】
本発明の別の形態において、気道の開存性は、未定義、または「未決定」領域内にあることが確定される。本発明の別の態様は、開存性は、開放性から閉塞性への段階的な尺度を有すると定義されることである。別の態様は、無呼吸の識別領域は、複素平面内に画定され、領域は必ずしも隣接していないことである。さらなる改良点は、識別領域は、平均漏出レベルの関数であり得て、具体的には、「未決定」領域は、恐らく漏出の増加とともにその大きさが増大することである。
【００２６】
本発明のさらなる形態は、特許請求の範囲に示される通りである。
【００２７】
本発明の実施形態が、ここで添付図面を参照して説明される。
【図面の簡単な説明】
【００２８】
【図１】実施形態の基本的技法のフロー図を示す。
【図２】本発明を具体化する器具を図式で示す。
【図３】本発明の２ポート解析の電子的な類似形態を示す。
【図４】詳細なアルゴリズムに基づいた、複素アドミッタンスの値を示す。
【図５】単純化アルゴリズムに基づいた、複素アドミッタンスの値を示す。
【発明を実施するための形態】
【００２９】
本発明を実施するシステム
本発明は、全体として参照することによって組み込まれる、米国特許第５，７０４，３４５号に開示される実施形態の改良形態である。図１は、本発明の一実施形態の基本的技法のフロー図である。本発明の第１のステップ１０は、無呼吸が発生している間の、フロージェネレータ近傍の点での、経時的な呼吸流量および圧力の測定である。この情報は、その後の処理のための定性的測定値として使用されるアドミッタンス値を生成するために、ステップ１２において処理される。ステップ１４〜１６は、時間窓内の複素アドミッタンス値を開存性領域と比較することにより、閉塞性、開放性、または混合性無呼吸が発生しているかどうかを検出する。
【００３０】
無呼吸が進行中である場合、次に、無呼吸が開放性であるか、または閉塞性であるかの確定が続く。開放性気道を有する無呼吸が発生している場合、かつ必要に応じて、事象は、ステップ１８で記録される。ステップ１６の結果が、閉塞性気道を有する無呼吸が発生しているとする場合、ＣＰＡＰ治療圧の増加が、ステップ２０で行われる。必要に応じて、ステップ２０は、検出された異常を任意に記録することを含み得る。
【００３１】
開放性気道を有する無呼吸の場合、無呼吸の発生を未然に防ぐか、または少なくとも軽減するために必要とされる、最小圧力を設定しようとする一般的方法論に従い、ＣＰＡＰ治療圧は減少される。ステップ２２における減少量は、必要に応じて、ゼロであり得る。
【００３２】
図１に表される技法は、患者のＣＰＡＰ圧力が、適宜に経時的に制御される、臨床的実施形態の技法である。純診断的実施形態は、それぞれ、ステップ２０およびステップ２２における、ＣＰＡＰ圧力増加および圧力減少の動作を省略する以外は、同一の方法で機能する。
【００３３】
図２は、図１の技法を実施するための一実施形態に従った、臨床ＣＰＡＰ器具を図式で示す。マスク３０は、鼻マスク、および／または顔マスクのいずれであっても、患者の顔に密閉して係合される。新鮮な空気または酸素富化空気は、可撓性管３２を介してマスク３０に入り、可撓性管３２は、次いで、電動タービン（フロージェネレータ）３４に接続され、そこに、空気入口３６が提供される。ＣＰＡＰ治療として、マスク３０に供給される空気の圧力を増加あるいは減少させるために、タービン用のモータ３８は、モータサーボユニット４０により制御される。マスク３０はまた、管３４とマスク３０との接合部に近接する排気ポート４２を含む。
【００３４】
流動抵抗要素４４は、フロージェネレータ３４に隣接する。これは、差圧が存在するアイリス絞り形状を取ることが可能である。流動抵抗要素４４のマスク側は、小口径管４６を介して、圧力トランデューサ４８と、差圧トランデューサ５０の入力に接続される。流動抵抗要素４４の反対側の圧力は、もう１つの小口径管５２を介して、差圧トランデューサ５０の他の入力に移動する。
【００３５】
圧力トランデューサ４８は、流動圧力に比例した電気信号を生成し、それは、増幅器５６により増幅され、マルチプレクサ／ＡＤＣユニット５８とモータサーボユニット４０との両方に移動する。モータサーボユニット４０に提供される信号は、静圧が、定値圧力に近似するように制御されることを確実にするための、フィードバックとしての機能を果たす。
【００３６】
流動抵抗要素４４にわたって検知される差圧は、差圧トランデューサ５０から電気信号として出力され、もう１つの増幅器６０により増幅される。したがって、増幅器５６からの出力信号は、マスクまたは呼吸の空気流量の測定値を表す。制御装置６２は、多くの処理機能を実行するようにプログラムされる。
【００３７】
圧力および流量は、定常値およびＡＣ値から成ると考えられ得て、後者は、４Ｈｚの周波数を有する圧力への強制振動信号の影響を示す。以下において、「定常」量は、（ａ）例えば、フィルタリング操作により、振動成分が意図的に除去されている量、または（ｂ）振動成分が除去されている量に部分的または完全に基づいた、計算結果、または（ｃ）フロージェネレータで測定される圧力および流量等の瞬間量（したがって、振動成分を含む）に基づき計算される量のいずれか、およびシステムの無効分を部分的あるいは完全に無視し、それを、単に（場合により非線形の）抵抗システムとして扱う、空気通路および患者の漏出のモデルである。管３２に沿った定常圧力損失は、チューブを通る流量、および例えばテーブル索引による、管の静圧流量特性の知識から計算される。次に、マスクにおける定常圧力は、チューブの圧力損失を引くことにより計算される。次に、チューブ３２に沿った圧力損失は、圧力ジェネレータ３４で望ましい瞬間圧力を得るために、マスクにおける望ましい設定圧力に加えられる。排気口４２を通る流量は、例えばテーブル索引による、排気口の圧力流量特性からの、マスクにおける圧力から計算される。定常マスク流量は、管３２を通る流量から、排気口４２を通る流量を引くことにより計算される。次に、定常患者流量は、例えば、入力が瞬間マスク流量である、１０秒の時定数を有する一次ローパスフィルタにより、ＣＰＡＰ装置において確定され得る推定定常漏出を、定常マスク流量から引くことにより計算される。
【００３８】
制御装置６２により実行される技法を、ここで図２の器具を参照して説明する。患者呼吸流量が、非常に低い、またはゼロである場合（漏出があれば、患者が無呼吸である場合、マスクの流動は停止しないことに留意する）、外部から誘導された振動技術を使用することにより、気道開存性の確定（ステップ１４〜１６）が行われる。気道が開放性であるが、呼吸筋が弛緩している場合（すなわち、開放性気道を有する中枢性無呼吸）、外部に起因するマスク圧力の小さい変動が、肺を膨張および収縮させることにより、かつ肺内のガスを圧縮および減圧することにより、小さい呼吸気流を誘導する。反対に、気道が閉塞性である場合、空気流は誘導されない。これは、以下のように定量化される。
【００３９】
アドミッタンスの閾値による無呼吸の識別
アドミッタンスＹは、
Ｙ＝Ｇ＋ｉＢ
により与えられ、Ｇは、コンダクタンスであり、Ｂは、サセプタンスである。
【００４０】
気道が開放性であるかどうかを判断するために、患者のアドミッタンスＹの複素数値は、値の領域と比較される。この閾値の値は、以下の考察に基づき選択され得る。
【００４１】
マグニチュードよりもよりよい分類指標である角度に関する説明は、回路内の複素インピーダンスの相互作用に関連付けられ得る。漏出および換気流は、本質的に複素（ｊ）成分がない非線形抵抗である。患者は、基本的に接地された、直列の抵抗およびコンプライアンス（キャパシタンス）である。ホース圧力降下もまた、流動のイナータンスの形で複素成分を有する。したがって、インピーダンス全体の大きさは、すべての成分（ホース、漏出、換気流、患者の気道）が実部を有するため、それらに影響される。しかしながら、虚数成分は、慣性ホース圧力降下と容量肺との間の相互作用である。それは、慣性と容量性リアクタンスとの相対的な「強度」を比較している。
【００４２】
また、角度および振幅ではなく、インピーダンスの実部および虚部を使用することも可能である。これは、ＣＰＵ負荷の減少という利点を有する。振幅の計算は、平方根を必要とし、角度の計算は、逆タンジェントを伴う。（これらの両方は、計算コストの低い、ルックアップテーブルを伴う技術により、必要とされる精度で容易に行われ得る。）実および虚の両方、または振幅および角度の両方を使用することにより、よりロバストな分類が得られる。開放性および閉塞性の値のグループ化は、はるかに明らかとなる。
【００４３】
図４の値から確定される通り、可能な閾値の組は、
角度＝ＡＲＧ（患者＿アドミッタンス）
ｉｆ（角度＜−２ ＯＲ 角度＞１．２）
気道状態＝閉塞性；
ｅｌｓｅ ｉｆ（角度＜０．９ ＡＮＤ 角度＞−１．５）
軌道状態＝開放性；
ｅｌｓｅ
気道状態＝未定義；
ｅｎｄ ｉｆ
となり、角度は、ラジアン単位である。閾値に対する許容値は、図４から確定される開放性および閉塞性の無呼吸を表す値の区別を可能にする任意の値であり得る。
【００４４】
４Ｈｚにおける患者抵抗も、気道の状態を示す。この抵抗は、上記の方程式におけるコンダクタンスの逆数ではない。４Ｈｚにおける患者のインピーダンスは、
Ｚ＝Ｒ＋ｉＸ
であり、Ｙ＝１／Ｚである。患者抵抗Ｒは、
Ｒ＝Ｒｅ（Ｚ）＝Ｒｅ（１／Ｙ）
により与えられる。
【００４５】
代替実施形態では、Ｒが、｜Ｙ｜よりも、気道の状態を特徴付けるために使用され得る。
【００４６】
アドミッタンスの計算
患者のアドミッタンスは、方程式
【数１】
により計算され、分子は、ＡＣ患者流量のマグニチュードであり、実質的に患者流量の（時間に関する）微分の測定値である。分母は、マスク圧力の同一の微分である。
【００４７】
換気流および漏出の補正
患者の空気流量は、マスクへの流入量から、換気口を通る流量と漏れ流量を引くことにより確定される。
【数２】
【００４８】
平均換気流量は、測定中の平均マスク圧力、またはより好ましくは、測定中のマスク圧力の平方根の平均の平方（漏出の計算を行うためにも使用され得る）のいずれかから、確定される。マスク圧力自体は、空気通路のＡＣ動作を無視して、従来の方法で計算される。
【００４９】
ＡＣ換気流は、動作点の線形近似を使用して計算される。
【数３】
【００５０】
マスクにおけるＡＣ流量を確定するための、マスク圧力のモデリング
マスク圧力は、換気流の関数として換気口を通じた圧力降下を得る方程式において、係数ｋ１およびｋ２を求めることにより確定される。
【数４】
したがって、特定のマスク圧力における換気流を計算するために、そのマスク圧力における換気流に対する二次方程式を解く。
【００５１】
【数５】
は、先の方程式の微分により得られ、
【数６】
を与える。
【００５２】
【数７】
であるため、小信号の近似は、
【数８】
を得て、そこから、
【数９】
が、妥当な近似として得られる。
【００５３】
漏れ流量の計算
患者の空気流の計算を完全にするために、マスクへの流入からの漏れ流量を計算することも必要である。実際の漏出は、患者の４Ｈｚのアドミッタンスが計算される間、計算される。典型的な漏れ係数推定は、少し前に起こったことに関連し、従来技術のＡｕｔｏＳｅｔ ＣＳおよびＡｕｔｏＶＰＡＰ装置におけるように、
【数１０】
から漏れ係数Ｋｌｅａｋを推定し得て、上線は、測定期間にわたる平均値を示す。無呼吸中、漏出は、全非換気流量に等しく、
【数１１】
呼吸中、漏出は、平均非換気流に（よい近似で）等しい。上記のＫｌｅａｋの推定において、すべての量は、瞬間値に基づき、空気通路のＡＣ特性を参照することなく計算され、空気通路を非線形抵抗器として扱う。
【００５４】
ＤＣ（Ｋｌｅａｋに対する上記の式に基づく）とＡＣとの両方の、瞬間漏出に対するモデルは、
【数１２】
であり、換気流と同様に、漏れ流量の（４Ｈｚにおける）ＡＣ成分は、小信号近似
【数１３】
を使用して得られ、上記の方程式を直接微分することにより、
【数１４】
を計算し得て、
【数１５】
を得る。
【００５５】
【数１６】
の平均は、
【数１７】
または
【数１８】
のいずれかであり得る。後者が好ましく（Ｐｍａｓｋが定数ではない場合、それが、
【数１９】
の平均値の不偏推定値を得るため、かつＫｌｅａｋを推定するために、それをすでに計算したため）、
【数２０】
を得る。
【００５６】
ＡＣ値の確定
この説明を通して、一貫して参照されるＡＣ値は、測定から、励振周波数における三角フーリエ成分として確定される。これは、正規直交関数の組として方形波を使用していた強制振動法の、従来技術の使用の改良形態である。
【００５７】
すべてのＡＣ量は、複素数として計算される。具体的には、マスクにおけるＡＣ圧力および流量は、正弦および余弦関数のそれぞれを有する積和を使用して、測定中、瞬間圧力および流量の正弦および余弦成分を得ることにより計算され、それにより、成分係数ｃｓおよびｃｃを得て、次に、ＡＣ値としてｃｓ＋ｉｃｃまたはｃｃ＋ｉｃｓのいずれかを記述する。（これらの２つの式のいずれも使用され、２つの式は、単に位相シフトの点で異なるが、１つの式は、初めから終わりまで一貫して使用される）。あるいは、正弦および余弦の最小二乗法等の、励振周波数における圧力および流量の正弦波成分の振幅および位相を推定する、他の標準的方法が使用され得る。
【００５８】
明確には、
【数２１】
を得、励振周波数は、Ｔ／２Ｂであり、Ｔは、正規化係数であり、積分は、正弦および余弦の離散的な標本値を使用した総和により置き換えられ得る。ｆ（ｔ）は、フーリエ係数が必要とされるいかなる関数でもある。
【００５９】
すべての量は、ｔ＝０として象徴的に示される、６秒スライディングウィンドウ（「アドミッタンスウィンドウ」）上で確定される。アドミッタンス計算は、原則として、アルゴリズムのサンプリング周波数（例えば、５０Ｈｚ）で行われ得るが、臨床目的では、これは必須ではなく、計算は、毎４Ｈｚサイクルの終わりに、効果的に実行され得る。あるいは、２Ｈｚ、あるいは１Ｈｚ周波数においてでも、アドミッタンスの計算は、臨床目的では妥当である。
【００６０】
モータ制御装置の遅延により、圧力波形の第１のサイクル（２５０ｍｓ）は、正弦波ではない場合があり、空気通路内の定常の設定に多少の遅延が出るため、データの最初の２５０ｍｓは、無視され得る。
【００６１】
正弦および余弦値は、起動時に生成されるテーブル内に記憶されるか、またはソースコード内の定数として特定され得る。もちろん、合計のオーバーフローが発生しないと仮定して、乗算は、３２ビットの結果を有する、１６ビット＊１６ビットの固定小数点であり、かつこの方法の数値安定性において問題がない場合、十分な精度が提供される。したがって、計算コストは、比較的小さい。
【００６２】
２ポートネットワーク上のインピーダンスのモデリング
マスクにおけるＡＣ圧力およびマスクに入るＡＣ流量の計算は、一実施形態では、特定の周波数における、システム内の流動の２ポート電気ネットワークへの類似から確定される。
【００６３】
図３から、ルーチン回路解析は、
【数２２】
を得て、もちろん、これらは、マスクにおける圧力、ならびにマスクに入る流量のそれぞれに対応する。明確には、
【数２３】
である。
【００６４】
インピーダンスＺｉｎｏｕｔおよびＺｏｕｔｇｎｄは、平均流量の関数とみなされる。これらのインピーダンスは、多くの平均流量レベルで、特定の空気通路に対して測定される（詳細は以下を参照）。アドミッタンスウィンドウ中の、またＺｉｎｏｕｔおよびＺｏｕｔｇｎｄのそれぞれに対する、フロージェネレータの総流量の平均が、測定されたインピーダンス間を線形補間してその平均流量レベルにおけるインピーダンスを確定するために使用される。
【００６５】
入口および出口のＡＣ流量およびＡＣ圧力は、フーリエ係数の計算に関して上記に記載される標準的な積和方法を使用して、（ノイズを減少させるために）３０秒間にわたって確定され、複素数値を得る。これらの値は、以下により２ポートパラメータを計算するために使用される（図１を参照。その図中の電圧および電流は、以下における圧力および流量のそれぞれに対応する）。
【数２４】
分母が非常に小さいか、またはゼロである場合、インピーダンスは、典型的な空気通路のインピーダンスとの関連で、特定の数値的に非常に大きい値を有するとみなされる。
【００６６】
ホース圧力降下の理論モデル
ホース圧力降下は、マスク内圧とＦＧ圧力との間の圧力差である。それは、マスクの前に配置されるすべての要素にわたる圧力損失である。これは、マフラ、加湿器、ＡＢフィルタ、およびマスク接続ホースを含む。
【００６７】
完全なホース圧力降下は、抵抗ホース圧力降下および慣性圧力降下で構成される。
【数２５】
Ｋ１、Ｋ２＝実験的に導かれた定数
ＫＬ＝理論的に決定された流量のイナータンス定数
ＱＦＧ＝ＦＧ流量
ＫＬ＝ρｌ／Ａ
ρ＝空気密度（１．１９ｋｇ／ｍ３）
ｌ＝チューブ長（２または３ｍ）
Ａ＝断面積（πｄ２／４、ｄ＝０．０１９ｍ）
【００６８】
ＲｅｓＭｅｄ ＡｃｔｉｖａおよびＳｗｉｆｔマスクなどのマスクの場合、主要ホースとマスクとの間に別のチューブがある。これもまた、慣性定数において考慮されなければならない。
【００６９】
平均ホース圧力降下は、抵抗部のみで構成される（上線は平均を示す）。
【数２６】
Ｋ１、Ｋ２＝実験的に導かれた定数
【数２７】
【００７０】
ＡＣホース圧力降下を計算するために、ある種の線形化が必要である。抵抗および慣性成分は分離される。抵抗成分は、動作点に関して、上記の平均ホース圧力降下の式を線形化することにより得られる。
【数２８】
【００７１】
したがって、ＡＣホース圧力降下は、ＡＣ流量に、周囲に微小振動が発生する位置における、ホースの非線形抵抗の勾配を乗じたものである。導出のために、振動は、ホース抵抗の勾配の変化と比較して、小さい振幅の振動であると仮定する。
【００７２】
慣性成分は、ＦＧ流量を導関数で乗じる定数である。標準的な線形回路理論により、インダクタンスＬは、インピーダンスｓＬを有し、周波数ωにおける正弦波信号に対して、ｊωＬである。次に、慣性ホース圧力降下の成分は、
【数２９】
となり、２つの組み合わせは、
【数３０】
を得る。
【００７３】
マスク圧力
マスク圧力は、ＦＧ圧力からホース圧力降下を減じたものである。
Ｐｍａｓｋ＝ＰＦＧ−Ｐｈｏｓｅｄｒｏｐ
Ｐｍａｓｋ，ＡＣ＝ＰＦＧ，ＡＣ−Ｐｈｏｓｅｄｒｏｐ，ＡＣ
【実施例】
【００７４】
図４では、太線が、開放性無呼吸を表し、細線が、閉塞性無呼吸を表すグラフが表される。横軸が、アドミッタンスの絶対マグニチュードを表す一方で、縦軸は、複素アドミッタンスの位相を表す。図面から理解され得るように、１つは起点近傍、他方はｙ軸の最上値にある、２つのより細い曲線領域により表される閉塞性無呼吸と、グラフの中心の太い曲線により表される開放性無呼吸が完全に分離されている。実質的に、データのすべてをｘ軸に投影することにより、アドミッタンスの絶対値のみを見ると、細い（閉塞性無呼吸）曲線が、太い（開放性無呼吸）データと接触し、データの分離は、失われる。
【００７５】
図４は、上記の式に従って計算された患者のアドミッタンスデータの、プロッティング結果を示す。ｘ軸の値は、患者のアドミッタンスの平均絶対値である。ｙ軸の値は、患者のアドミッタンスの位相角である。図４から理解され得るように、閉塞性無呼吸および開放性無呼吸に対する値の完全な分離が見られる。
【００７６】
図５は、先のアドミッタンスの計算の代わりに、複素アドミッタンスに対する値が、モデリングにおけるマスクまたは換気回路の任意の特定構成の詳細から、本質的に独立したモデルから確定される、患者のアドミッタンスのプロットの結果を示す。具体的には、図５のデータは、単純漏出オリフィスモデル、Ｐ＝（ＫＱ）２を使用することにより確定された。開放性および閉塞性の無呼吸は、複素アドミッタンスのグラフにおいて分離した領域をもたらすことが、図５から理解され得る。
【００７７】
適度な（近似された患者のアドミッタンス）アルゴリズム
第３のアルゴリズムは、患者回路を総称的にモデリングし、漏出および換気流を組み合わせることにより、患者のアドミッタンスを近似しようと試みる。
【００７８】
患者回路は、上記の複素アルゴリズムと同様の方法でモデリングされるが、定数は、ＡＢフィルタを有する２ｍホースに対する定数で固定される。
【００７９】
漏出および換気流のモデリングは、非患者流量である１つのパラメータに組み合わされる。これは、
漏出＋換気流
としてモデリングされる。
【００８０】
漏出および換気流は、組み合わされ、以下の方程式を使用する際にモデリングされる。
【００８１】
漏出オリフィスは、
【数３１】
を使用して、ランタイムで近似され得る。
【００８２】
次に、瞬間漏出は、
【数３２】
を使用して計算される。
【００８３】
小信号近似を使用すると、
【数３３】
であると言うことができ、式中
【数３４】
であり、Ｋｌｅａｋに置換して、
【数３５】
を得て、小さい平均Ｐｍａｓｋの変化を扱っているので、
【数３６】
で近似することができ、したがって、ＡＣ漏出は、
【数３７】
となる。
【００８４】
近似された患者流量
次に、患者流量は、
Ｑｐａｔｉｅｎｔ，ＡＣ＝ＱＦＤ，ＡＣ−Ｑｌｅａｋ，ＡＣ
である。
【００８５】
Ｑｌｅａｋは、換気流と組み合わされる漏れ流量である（すべての非患者流量）。
【００８６】
分類
近似アドミッタンスを使用して、気道状態を分類する多くの方法がある。アドミッタンスの角度対マグニチュードをプロットすると、角度のみにおけるよい分離を見出し得る。したがって、単純分類は、２つの閾値を有するものとなるが、純粋に角度の上のみである。しかしながら、角度の計算は、除算および逆タンジェントを伴い、さらなる計算を必要とする。代わりに、使用されるアドミッタンスの実部および虚部を有することが好ましいことが示唆される。
【００８７】
閾値
よい分類精度を提供する最も単純な閾値は、単純な対角直線である。これらは、ｙ＞ｘ＋ｂおよびｙ＜ｘ＋ｃの式である。この場合、アドミッタンスの虚部（ｙ）対実部（ｘ）をプロットすると、閉塞性は、ｙ＞ｘ＋ｂであり、開放性は、ｙ＜ｘ＋ｃである。浮動小数点アルゴリズムに基づき、定数に対する値は、ｂ＝０、ｃ＝−０．０３であった。これは、何の誤分類も与えない。データ領域を分離させる他の任意の曲線が使用され得るが、線形分離は、閾値計算の複雑さを減少させ、その理由から好ましい。
【特許請求の範囲】
【請求項１】
無呼吸の発生を確定するための方法であって、
空気流量信号を提供するために、患者の気道から呼吸気流量を測定するステップと、
前記患者の気道に既知の周波数の振動圧力波形を印加するステップと、
患者のアドミッタンスを表す複素量を計算するステップと、
前記複素量の値を、開存性または閉塞性気道を示す領域と比較するステップと、を含む、方法。
【請求項２】
患者のアドミッタンスを、正常呼吸時のアドミッタンスと比較するステップをさらに含む、請求項１に記載の無呼吸の発生を確定するための方法。
【請求項３】
マスクとフロージェネレータとを有するＣＰＡＰ器具を使用して、患者における呼吸器官の閉塞性および開放性無呼吸を識別するための方法であって、
前記フロージェネレータで、特定の振幅の波形を有する圧力波を印加するステップと、
前記マスクで、前記圧力波形の振幅を確定するステップと、
前記マスクで望ましい振幅を生成するために、前記フロージェネレータで、前記圧力振幅を調整するステップと、
前記患者のアドミッタンスを確定し、それを閉塞性または開放性の無呼吸を示す値と比較するステップと、を含む、方法。
【請求項４】
前記圧力振幅を調整する前記ステップは、前記マスクにおける前記実際の圧力のマグニチュードに対する前記望ましい圧力のマグニチュードの比で、前記フロージェネレータにおける前記圧力を増加させることにより達成される、請求項３に記載の呼吸器官の閉塞性および開放性無呼吸を識別するための方法。
【請求項５】
患者の気道の開放性および閉塞性無呼吸を識別する方法であって、
（ｉ）空気送達チューブと患者との接触面を介して、患者に呼吸装置を接続するステップと、
（ｉｉ）前記空気送達チューブを介して前記呼吸装置から前記患者に、前記患者への交流圧力波形を送達するステップと、
（ｉｉｉ）前記呼吸装置において空気の流量と圧力とを測定するステップと、
（ｉｖ）空気送達チューブのインピーダンスの容量分を確定するステップと、
（ｖ）前記測定された空気の流量と圧力ならびに前記容量分から、患者のアドミッタンスを確定するステップと、
（ｖｉ）前記患者のアドミッタンスに基づいて、気道の開放性および閉塞性無呼吸を識別するステップと、
を含む、方法。
【請求項６】
漏出に対して前記流量を修正するステップをさらに含む、請求項５に記載の方法。
【請求項７】
換気流に対して前記流量を修正するステップをさらに含む、請求項５に記載の方法。
【請求項８】
患者の気道の開放性および閉塞性無呼吸を識別する方法であって、
（ｉ）空気送達チューブと患者との接触面を介して、患者に呼吸装置を接続するステップと、
（ｉｉ）前記空気送達チューブを介して前記呼吸装置から前記患者に、前記患者への交流圧力波形を送達するステップと、
（ｉｉｉ）前記呼吸装置における空気の流量と圧力とを測定するステップと、
（ｉｖ）空気送達チューブのインピーダンスの誘導分を確定するステップと、
（ｖ）空気の前記測定された流量と圧力ならびに前記容量分から、患者のアドミッタンスを確定するステップと、
（ｖｉ）前記患者のアドミッタンスに基づいて、気道の開放性および閉塞性無呼吸を識別するステップと、
を含む、方法。
【請求項９】
漏出に対して前記流量を修正するステップをさらに含む、請求項８に記載の方法。
【請求項１０】
換気流に対して前記流量を修正するステップをさらに含む、請求項８に記載の方法。
【請求項１１】
前記アドミッタンスの位相角を確定するステップをさらに含む、請求項５または請求項８に記載の方法。
【請求項１２】
前記アドミッタンスの大きさを確定するステップをさらに含む、請求項５または請求項８に記載の方法。
【請求項１３】
前記アドミッタンスの実成分を確定するステップをさらに含む、請求項５または請求項８に記載の方法。
【請求項１４】
前記アドミッタンスの虚数成分を確定するステップをさらに含む、請求項５または請求項８に記載の方法。
【請求項１５】
前記交流圧力波形は、約４Ｈｚの周波数を有する、請求項５または請求項８に記載の方法。
【請求項１６】
前記交流圧力波形は、約２〜８Ｈｚの範囲の周波数を有する、請求項５または請求項８に記載の方法。
【請求項１７】
前記交流圧力波形は、約１〜１６Ｈｚの範囲の周波数を有する、請求項５または請求項８に記載の方法。
【請求項１８】
前記圧力波形の印加に続く、最初の約２５０ｍｓの流量または圧力データは、測定されない、請求項５または請求項８に記載の方法。
【請求項１】
無呼吸の発生を確定するための方法であって、
空気流量信号を提供するために、患者の気道から呼吸気流量を測定するステップと、
前記患者の気道に既知の周波数の振動圧力波形を印加するステップと、
患者のアドミッタンスを表す複素量を計算するステップと、
前記複素量の値を、開存性または閉塞性気道を示す領域と比較するステップと、を含む、方法。
【請求項２】
患者のアドミッタンスを、正常呼吸時のアドミッタンスと比較するステップをさらに含む、請求項１に記載の無呼吸の発生を確定するための方法。
【請求項３】
マスクとフロージェネレータとを有するＣＰＡＰ器具を使用して、患者における呼吸器官の閉塞性および開放性無呼吸を識別するための方法であって、
前記フロージェネレータで、特定の振幅の波形を有する圧力波を印加するステップと、
前記マスクで、前記圧力波形の振幅を確定するステップと、
前記マスクで望ましい振幅を生成するために、前記フロージェネレータで、前記圧力振幅を調整するステップと、
前記患者のアドミッタンスを確定し、それを閉塞性または開放性の無呼吸を示す値と比較するステップと、を含む、方法。
【請求項４】
前記圧力振幅を調整する前記ステップは、前記マスクにおける前記実際の圧力のマグニチュードに対する前記望ましい圧力のマグニチュードの比で、前記フロージェネレータにおける前記圧力を増加させることにより達成される、請求項３に記載の呼吸器官の閉塞性および開放性無呼吸を識別するための方法。
【請求項５】
患者の気道の開放性および閉塞性無呼吸を識別する方法であって、
（ｉ）空気送達チューブと患者との接触面を介して、患者に呼吸装置を接続するステップと、
（ｉｉ）前記空気送達チューブを介して前記呼吸装置から前記患者に、前記患者への交流圧力波形を送達するステップと、
（ｉｉｉ）前記呼吸装置において空気の流量と圧力とを測定するステップと、
（ｉｖ）空気送達チューブのインピーダンスの容量分を確定するステップと、
（ｖ）前記測定された空気の流量と圧力ならびに前記容量分から、患者のアドミッタンスを確定するステップと、
（ｖｉ）前記患者のアドミッタンスに基づいて、気道の開放性および閉塞性無呼吸を識別するステップと、
を含む、方法。
【請求項６】
漏出に対して前記流量を修正するステップをさらに含む、請求項５に記載の方法。
【請求項７】
換気流に対して前記流量を修正するステップをさらに含む、請求項５に記載の方法。
【請求項８】
患者の気道の開放性および閉塞性無呼吸を識別する方法であって、
（ｉ）空気送達チューブと患者との接触面を介して、患者に呼吸装置を接続するステップと、
（ｉｉ）前記空気送達チューブを介して前記呼吸装置から前記患者に、前記患者への交流圧力波形を送達するステップと、
（ｉｉｉ）前記呼吸装置における空気の流量と圧力とを測定するステップと、
（ｉｖ）空気送達チューブのインピーダンスの誘導分を確定するステップと、
（ｖ）空気の前記測定された流量と圧力ならびに前記容量分から、患者のアドミッタンスを確定するステップと、
（ｖｉ）前記患者のアドミッタンスに基づいて、気道の開放性および閉塞性無呼吸を識別するステップと、
を含む、方法。
【請求項９】
漏出に対して前記流量を修正するステップをさらに含む、請求項８に記載の方法。
【請求項１０】
換気流に対して前記流量を修正するステップをさらに含む、請求項８に記載の方法。
【請求項１１】
前記アドミッタンスの位相角を確定するステップをさらに含む、請求項５または請求項８に記載の方法。
【請求項１２】
前記アドミッタンスの大きさを確定するステップをさらに含む、請求項５または請求項８に記載の方法。
【請求項１３】
前記アドミッタンスの実成分を確定するステップをさらに含む、請求項５または請求項８に記載の方法。
【請求項１４】
前記アドミッタンスの虚数成分を確定するステップをさらに含む、請求項５または請求項８に記載の方法。
【請求項１５】
前記交流圧力波形は、約４Ｈｚの周波数を有する、請求項５または請求項８に記載の方法。
【請求項１６】
前記交流圧力波形は、約２〜８Ｈｚの範囲の周波数を有する、請求項５または請求項８に記載の方法。
【請求項１７】
前記交流圧力波形は、約１〜１６Ｈｚの範囲の周波数を有する、請求項５または請求項８に記載の方法。
【請求項１８】
前記圧力波形の印加に続く、最初の約２５０ｍｓの流量または圧力データは、測定されない、請求項５または請求項８に記載の方法。
【図１】
【図２】
【図３】
【図４】
【図５】
【図２】
【図３】
【図４】
【図５】
【公開番号】特開２０１２−８６０２５（Ｐ２０１２−８６０２５Ａ）
【公開日】平成２４年５月１０日（２０１２．５．１０）
【国際特許分類】
【外国語出願】
【出願番号】特願２０１１−２４３７８７（Ｐ２０１１−２４３７８７）
【出願日】平成２３年１１月７日（２０１１．１１．７）
【分割の表示】特願２００９−５２５８５８（Ｐ２００９−５２５８５８）の分割
【原出願日】平成１９年８月３０日（２００７．８．３０）
【出願人】（５０００４６４５０）レスメド・リミテッド (192)
【氏名又は名称原語表記】ＲＥＳＭＥＤ ＬＴＤ
【Ｆターム（参考）】
【公開日】平成２４年５月１０日（２０１２．５．１０）
【国際特許分類】
【出願番号】特願２０１１−２４３７８７（Ｐ２０１１−２４３７８７）
【出願日】平成２３年１１月７日（２０１１．１１．７）
【分割の表示】特願２００９−５２５８５８（Ｐ２００９−５２５８５８）の分割
【原出願日】平成１９年８月３０日（２００７．８．３０）
【出願人】（５０００４６４５０）レスメド・リミテッド (192)
【氏名又は名称原語表記】ＲＥＳＭＥＤ ＬＴＤ
【Ｆターム（参考）】
[ Back to top ]