睡眠時無呼吸症候群の測定装置
【課題】無拘束・無意識状態で、被験者の覚醒／睡眠を簡単に判定すると共に睡眠時における異常呼吸を高精度に判定する。
【解決手段】睡眠時無呼吸症候群測定装置は、エアマット１０からマイクロフォン３０ａ、圧力センサ３０ｂを介して、被験者の生体信号４０ｃと圧力信号４０ｉを検出し、ＣＰＵ４０ｇで解析する。ＲＯＭ４０ｊに格納されている覚醒／睡眠判定部４０ｋは圧力信号の体動波形から前回値との差分（動き量）波形を抽出し、差分波形の標準偏差を算出し、その一晩分のトレンドから波形の振幅と閾値を比較し、覚醒／睡眠を判定する。通気量推定部４０ｍは生体信号４０ｃから呼吸波形４０ｅを読み出し、その１サイクル分を検出し呼吸の通気量を推定し、通気量が閾値より小さい場合に異常呼吸と判定する。なお、包絡形状解析部４０ｌで１サイクル毎の呼吸波形の包絡線を測定し、包絡線が紡錘形の場合に、異常呼吸の可能性が大きいので通気量を推定する。
【解決手段】睡眠時無呼吸症候群測定装置は、エアマット１０からマイクロフォン３０ａ、圧力センサ３０ｂを介して、被験者の生体信号４０ｃと圧力信号４０ｉを検出し、ＣＰＵ４０ｇで解析する。ＲＯＭ４０ｊに格納されている覚醒／睡眠判定部４０ｋは圧力信号の体動波形から前回値との差分（動き量）波形を抽出し、差分波形の標準偏差を算出し、その一晩分のトレンドから波形の振幅と閾値を比較し、覚醒／睡眠を判定する。通気量推定部４０ｍは生体信号４０ｃから呼吸波形４０ｅを読み出し、その１サイクル分を検出し呼吸の通気量を推定し、通気量が閾値より小さい場合に異常呼吸と判定する。なお、包絡形状解析部４０ｌで１サイクル毎の呼吸波形の包絡線を測定し、包絡線が紡錘形の場合に、異常呼吸の可能性が大きいので通気量を推定する。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、人体にセンサを取り付けない無拘束・無意識状態で、睡眠時無呼吸症候群を測定する装置に関する。
【背景技術】
【０００２】
無拘束・無意識状態で睡眠時無呼吸症候群を測定する装置として特許文献１、２がある。特許文献１の検査装置は、就寝者の呼吸による体動に応じた荷重変化を呼吸信号として生成し、この呼吸信号の周波数の変化に基づいて、無呼吸状態もしくは低呼吸状態を判定する。特に、無呼吸もしくは低呼吸状態からの回復呼吸動作は通常呼吸に比較して速いため、呼吸信号の周波数が変化することから判定が可能になる。特許文献２の測定装置は、被験者の生体情報である呼吸数と心拍数といびきから特徴的要因を抽出し、その要因数に応じて無呼吸症候群を判定する。
【０００３】
【特許文献１】特開２００４−２４６８４号公報
【特許文献２】特開２００６−１６７２１１号公報
【発明の開示】
【発明が解決しようとする課題】
【０００４】
特許文献１では呼吸動曲線の周波数変化から無呼吸状態もしくは低呼吸状態を判定する際に、寝姿が仰臥位、腹臥位、または横臥位かを推定する必要があり、その寝姿によって検出精度に影響が出る。また特許文献２では呼吸数と心拍数といびきの特徴的要因から無呼吸を判定するが、構造が複雑でありかつ検出精度に問題がある。さらに、特許文献１では被験者の荷重を検出して在離床を判定し、離床のときは判定を行わない。しかし、被験者が在床であっても覚醒状態の場合が有り、このような場合の検出精度は大幅に低下してしまう恐れがある。
【０００５】
本発明の目的は、上記従来技術の問題点に鑑み、呼吸被験者の覚醒／睡眠を体動から判定するとともに、呼吸動波形および／またはいびき波形の通気度から異常呼吸（無呼吸状態もしくは低呼吸状態）を検出し、簡単かつ高精度に異常呼吸状態を判定する測定装置を提供することにある。
【課題を解決するための手段】
【０００６】
上記目的を達成するための本発明は、無拘束・無意識で被験者の無呼吸乃至低呼吸を判定する睡眠時無呼吸症候群測定装置において、被験者が横たわるエアマットの挙動から生体信号を検出し、前記生体信号である体動波形について、次々と前回検出分との差分を取って差分波形（動き量）を抽出し、該差分波形の標準偏差を所定時間にわたって求め、該標準偏差のトレンドにおける波形の振幅と閾値を比較して覚醒状態または睡眠状態を判定する覚醒／睡眠判定部を設け、前記睡眠状態の場合に、前記生体信号から抽出した体動波形および／またはいびき波形に基づいて、無呼吸乃至低呼吸となる異常呼吸を判定することを特徴とする。
【０００７】
あるいは、無拘束・無意識で被験者の無呼吸乃至低呼吸を判定する睡眠時無呼吸症候群測定装置において、被験者が横たわるエアマットの挙動から生体信号を検出し、前記生体信号である体動波形の差分を取って睡眠状態を判定し、前記睡眠状態の場合に、前記生体信号から呼吸波形を読み出し、その１サイクル分の通気量を推定し、前記通気量が閾値より小さい場合に異常呼吸と判定する通気量推定部を有することを特徴とする。
【０００８】
さらに、無拘束・無意識で被験者の無呼吸乃至低呼吸を判定する睡眠時無呼吸症候群測定装置において、被験者が横たわるエアマットの挙動から生体信号を検出し、前記生体信号である体動波形について、次々と前回検出分との差分を取って差分波形（動き量）を抽出し、該差分波形の標準偏差を所定時間にわたって求め、該標準偏差のトレンドにおける波形の振幅と閾値を比較して覚醒状態または睡眠状態を判定する覚醒／睡眠判定部と、前記生体信号から呼吸波形を読み出し、その１サイクル分の通気量を推定し、通気量が閾値より小さい場合に異常呼吸と判定する通気量推定部と、を有することを特徴とする。
【０００９】
また、前記生体信号からいびき波形を読み出し、その１サイクル分の通気量を推定し、呼吸波形といびき波形の通気量を足し合わせた通気量が閾値より小さい場合に異常呼吸と判定する通気量推定部と、を有することを特徴とする睡眠時無呼吸症候群測定装置。
【００１０】
また、前記１サイクル毎の呼吸波形の包絡線を測定し、該包絡線が紡錘形をなす場合に、前記通気量の推定を行うことを特徴とする。
【発明の効果】
【００１１】
本発明によれば、被験者の体動の変化量（動き量）の標準偏差から睡眠状態を判定するので被験者の覚醒／睡眠状態を簡単に判定できる。特に、標準偏差を取ることで動き量の大きさのばらつき具合を数値として求めることができるので、判定が正確になる。また、呼吸波形のサイクル毎に通気量を推定することで異常呼吸の有無が判定できる。この場合、呼吸動波形の包絡線の形状が紡錘型の場合は異常呼吸である可能性が高いので、この判定を踏まえることで高精度の判定が可能になる。
【発明を実施するための最良の形態】
【００１２】
本発明の実施形態について、以下、図１から図８を参照しながら説明する。図１は本発明の一実施例である睡眠時無呼吸症候群測定装置の全体構成を示している。睡眠時無呼吸症候群測定装置はエアマット１０を有しており、このエアマット１０には生体信号を精度よく検出するために、等間隔に配置された穴１０ａが設けられている。エアマット１０は中空チューブ２０を介してセンサ検出部３０に接続されている。
【００１３】
センサ検出部３０は中空チューブ２０を介して検出する。すなわち、呼吸時に生じる呼吸音「スーハー」やいびきによる「グゴー」という所謂「音」ではなく、呼吸時の胸や腹の動きやいびきによる微振動などを生体情報に関わる音圧信号として検出する。この生体情報から、呼吸信号といびき信号をマイクロフォン３０ａで検出し、体位の動きは圧力センサ３０ｂで検出する。マイクロフォン３０ａの信号は、増幅器（ＯＰＡＭＰ１）４０ａで増幅され音圧信号４０ｂになる。圧力センサ３０ｂの信号は増幅器（ＯＰＡＭＰ２）４０ｈで増幅され圧力信号４０ｉとして出力される。
【００１４】
音圧信号４０ｂはＡＤ変換器４０ｃにより呼吸デジタル信号４０ｅに変換される。他方、音圧信号４０ｂはバンドパスフィルタ４０ｄでフィルタリングされ、いびき信号４０ｆとして出力される。
【００１５】
呼吸信号４０ｅ、いびき信号４０ｆ及び体動信号４０ｉは、コンピュータで形成される異常呼吸判定装置（ＣＰＵ）４０ｇに入力され、ＲＯＭ４０ｊに格納されるプログラムによりデジタル処理される。ＲＯＭ４０ｊのプログラムには、覚醒／睡眠判定部４０ｋ、呼吸波形の包絡形状解析部４０ｌ、通気量推定部４０ｍがあり、バス４０ｓを介してＣＰＵ４０ｇに呼び出される。
【００１６】
ＣＰＵ４０ｇのバス４０ｓには、計算結果を一時記憶するＲＡＭ４０ｎ、アラーム用通信デバイスＳＩＯ４０ｏ、上位コンピュータとデータ交換する通信デバイスＬＡＮ４０ｑが接続されている。上位コンピュータは（図示省略）健康管理センターなどに置かれ、患者データを集中的に管理し、異常の監視が行われる。従って、患者はエアマット１０に横たわっているのみで生体情報に関わる音圧信号が観測され、睡眠中の異常呼吸を自動で検知できるので、患者負担が少なく、かつ高性能な睡眠時無呼吸症候群の測定装置を提供できる。
【００１７】
次に、異常呼吸判定装置４０ｇの各部の処理を説明する。図２は覚醒／睡眠判定処理のフローチャートを、図３は覚醒／睡眠判定処理に応じた波形の変化を示している。図２と図３を用いて、本発明による覚醒／睡眠判定処理を説明する。
【００１８】
まず、圧力センサ３０ｂを介してエアマット１０にかかる圧力を体動波形として読み出し（ｓ１０１）、次々と前回検出分との差分を取って差分波形を抽出する（ｓ１０２）。図３の波形データ(１)に示すように、所定サンプリング時間ごとに体動波形の測定データを観測し、前回分との差分である動き量波形(２)を抽出する。図の例では、０．００５秒ごとに体動等の圧力信号をサンプリングし、そのたびに前回（０．００５秒前）の体動の大きさとの差分を求めている。体動（マット圧）の動き量とは圧力差、すなわち速度の関数である。体動波形の差分を求めることで、バイアス成分を除去でき、閾値による判定が容易になる。
【００１９】
次に差分波形の１分間の標準偏差を（１）式により算出し（ｓ１０３）、これを所定期間、たとえば一晩分のデータについて求める（ｓ１０４）。
標準偏差＝√（（動き量ｋ−ｎ点の動き量平均）２／（ｎ点の動き量平均）） …（１）
ここで、ｋ＝０、１、２、・・・、ｎ−１である。標準偏差をとることで、動き量の大きさのばらつき具合を数値として求めることができる。従って、トレンドを求めるときに、動き量の平均値をプロットするよりも、動き量の標準偏差をプロットする方がより小さな動き量まで検知できる。
【００２０】
波形データ(３)は、１晩分の体動波形の動き量に対し標準偏差トレンドを求めた例で、覚醒と睡眠の違いが明瞭に表れている。
【００２１】
さらに、この標準偏差トレンドに対して移動平均を求めることで（ｓ１０５）、覚醒と睡眠の振幅の違いがより顕著になる。その移動平均された波形データ(４)に対し、閾値を設けることで、振幅より大きい場合は覚醒、振幅より小さい場合は睡眠と判定することができる（ｓ１０６）。
【００２２】
本実施例によれば、睡眠と判定された期間について後述する通気量推定を行い、異常呼吸の有無を判定できる。また、被験者がマットに横たわった状態で覚醒している場合には「覚醒」と判定し、誤って「睡眠」とすることはないので、誤判定を確実に防止できる顕著な効果がある。
【００２３】
なお、移動平均には単純移動平均、修正移動平均、荷重移動平均、それらの組み合わせ等、いずれの方法であっても良い。さらに、１分間の動き量波形(２)から体動の回数を求め、その一晩のトレンドを上記判定方法に加味しても良い。
【００２４】
また、判定の精度は落ちるが、体動動き量トレンド（体動動き量に対し、１分間の振幅平均値を一晩分求めたもの）、あるいは脈動波形や体動波形の振幅平均トレンドに対して閾値を設けても判定は可能である。図１の測定装置では脈動波形は検出していないが、ＡＤ変換器４０ｃを２チャンネル設けることにより、生体情報に関わる音圧信号４０ｂから脈動デジタル信号を検出し、脈動波形の１分間の振幅平均値を一晩分求めればよい。
【００２５】
図４は呼吸波形の正常・異常を判定する通気量推定部のフローチャートを示し、図５は通気量推定部の各処理における波形を示すグラフである。図４と図５を用いて、呼吸波形の正常・異常判定法を説明する。
【００２６】
被験者が息を吸う（吸気）と、胸や腹が膨らむためエアマット１０を押す力が大きくなり、また息を吐（呼気）くと、胸や腹が元の状態に戻るのでエアマット１０を押す力が弱まる。その結果、振幅が周期的に変化する呼吸波形を検知できる。一方、エアマット１０上に横たわった状態で、息をしないで胸や腹を膨張・収縮すると、息をしているときと同様の波形を得ることが知られている。このことから、従来のように波形の振幅の大きさを計測するだけでは、睡眠時無呼吸症候群のように胸や腹が動いても息をしていない症状を検知することはできない。
【００２７】
本発明では「通気」に着目して呼吸波形を観測する。通常の呼吸では十分に通気ができているため、安定した胸や腹の膨張・収縮を行い、呼吸波形の包絡線は平坦形となる。一方、数秒間に渡って低呼吸や無呼吸が続く異常呼吸では、被験者は息苦しいため、息をしようと胸や腹の膨張・収縮を繰り返し、十分な通気ができたところで正常の呼吸に戻る。また、一刻も早く息をしようとするため、周期が変動することもある。このような結果、異常呼吸の時間帯では、呼吸波形の包絡線は紡錘形となる。
【００２８】
通気量推定部４０ｍは、マイクロフォン３０ａを介して呼吸波形を読み出し（ｓ２０１）、呼吸波形の１サイクルを検出する（ｓ２０２）。図５の波形(１)に、呼吸波形の１サイクルを示す。１周期は呼気（山側）から吸気（谷側）の１サイクルである。
【００２９】
次に呼吸波形の各サイクルの包絡線を求め（ｓ２０３）、包絡線が平坦形ならば異常呼吸はないとみて処理を終了する（ｓ２０４）。この詳細については後述する。
【００３０】
一方、包絡線が紡錘形であれば、通気量を推定する（ｓ２０５）。通気量は呼吸波形１サイクルの波形率（横幅／高さ）と面積から求まり、これを全てのサイクルにおいて求める。図５の波形（２）に、呼吸１サイクルの通気量推定用パラメータ算出のしかたを示す。波形率は（２）式、通気量は（３）式から求まる。
【００３１】
波形率＝波幅／波高値 …（２）
通気量＝波形率×積算値（面積） …（３）
波形率は呼吸波形が丸いか尖っているかを判定する。波形率が小であれば波形が尖っていて、異常呼吸の可能性が高い。波形率が大であれば波形は丸みを帯び、異常呼吸の可能性は低い。図５の波形(３)、(４)に、１分間に含まれる呼吸波形全サイクルに対して求めた通気量の遷移を示す。波形(３)が正常呼吸時、波形(４)が異常呼吸時の通気量を示し、枠線の中が異常時呼吸である。
【００３２】
ここで、呼吸波形の包絡線解析部４０ｌの処理を説明する。図６は包絡線解析部のフローチャートを、図７はその解析波形を示している。以下、図６と図７に基づいて包絡線解析処理を説明する。
【００３３】
まず呼吸波形を読み出し（ｓ３０１）、山側のピークを線で連結する（ｓ３０２）。同様に、谷側のピークを線で連結する（ｓ３０３）。上記ｓ３０２、ｓ３０３の処理を任意の期間について実施し、包絡線の作成を行う（ｓ３０４、ｓ３０５）。
【００３４】
図７の呼吸波形データ(１)は１分間の呼吸波形で、山側および谷側の包絡線を求めている。これより山側および谷側の包絡線に対して極大値および極小値の有無を判定し、極大値および極小値を検知する毎に両数値の差を所定の閾値と比較し、包絡線の極大値および極小値の差が閾値より大きい場合は紡錘形、差が閾値より小さい場合は平坦形と判定する（ｓ３０６、ｓ３０７）。図７の波形(２)は平坦形であり、波形(３)、(４)は紡錘形の例で、右上がり矢印と右下がり矢印が連続する場合、紡錘形と判定される。このように、山側の包絡線または谷側の包絡線が紡錘形と判定される場合は、異常呼吸の可能性が高いので、包絡線が紡錘形であるか否かを判定することで、異常呼吸判定の精度を向上できる。
【００３５】
再び図４に戻る。次に、推定した通気量と閾値の比較から、閾値を上回れば正常呼吸、閾値を下回れば異常呼吸と判定する（ｓ２１０）。図５の波形(４)は異常呼吸が任意時間（ここでは１０秒）続いたとき、四角の枠で囲んで表示し、アラーム用通信デバイス４０ｏやＬＡＮ４０ｑに出力する。
【００３６】
上記の説明では、呼吸の通気量のみを推定している。しかし、バンドパスフィルタ４０ｄを介していびき波形を読み出し（ｓ２０６）、いびき波形の１波分を検出し（ｓ２０７）、そのサイクル毎に通気量を推定し（ｓ２０８）、その通気量を呼吸の通気量と足し合わせ（ｓ２０９）、この加算した通気量を閾値と比較するようにしても良い。
【００３７】
図８はいびき波形を示すグラフである。波形(１)はいびき波形の１サイクルを示す。いびき波形には周期性が無いため、便宜上、呼吸１サイクル中のいびきを１波として考える。波形(２)はいびき１波の通気量推定用パラメータ算出のしかたを示す。推定式は呼吸波形のものと同様で、上記（２）、（３）式に示した通りである。ただし、いびきの波形は１波の中に複数のピークが存在するため、波高には平均波高値を用いる。なお、波形(１)のように、呼吸の開始時刻といびきの開始時刻は必ずしも一致しない。そこで、位相成分を（４）式により求め、通気量を（５）式で求めるようにしても良い。
位相＝（（いびき開始時刻−呼吸開始時刻）／呼吸周期）×π …（４）
通気量＝波形率×積算値×位相 …（５）
【図面の簡単な説明】
【００３８】
【図１】本発明の一実施例による睡眠時無呼吸症候群測定装置の構成図。
【図２】一実施例による覚醒／睡眠を判定するフローチャート。
【図３】覚醒／睡眠判定処理における各処理の波形データを示すグラフ。
【図４】一実施例による通気量を推定するフローチャート。
【図５】通気量推定処理における各処理の波形データを示すグラフ。
【図６】一実施例による呼吸波形の包絡形状を解析するフローチャート。
【図７】呼吸包絡線の形状解析処理における各処理の波形データを示すグラフ。
【図８】いびき波形の１波と通気量推定用のパラメータを示すグラフ。
【符号の説明】
【００３９】
１０…エアマット、２０…チューブ、３０ａ…マイクロフォン、３０ｂ…圧力センサ、４０ａ，４０ｈ…オペアンプ、４０ｃ…ＡＤ変換器、４０ｄ…バンドパスフィルタ、４０ｇ…異常呼吸判定装置（ＣＰＵ）、４０ｊ…ＲＯＭ、４０ｋ…覚醒／睡眠判定部、４０ｌ…呼吸波形の包絡形状解析部、４０ｍ…通器量推定部、４０ｓ…バス、４０ｎ…ＲＡＭ、４０ｏ…アラーム用通信デバイス（ＳＩＯ）、４０ｑ…通信デバイスＬＡＮ。
【技術分野】
【０００１】
本発明は、人体にセンサを取り付けない無拘束・無意識状態で、睡眠時無呼吸症候群を測定する装置に関する。
【背景技術】
【０００２】
無拘束・無意識状態で睡眠時無呼吸症候群を測定する装置として特許文献１、２がある。特許文献１の検査装置は、就寝者の呼吸による体動に応じた荷重変化を呼吸信号として生成し、この呼吸信号の周波数の変化に基づいて、無呼吸状態もしくは低呼吸状態を判定する。特に、無呼吸もしくは低呼吸状態からの回復呼吸動作は通常呼吸に比較して速いため、呼吸信号の周波数が変化することから判定が可能になる。特許文献２の測定装置は、被験者の生体情報である呼吸数と心拍数といびきから特徴的要因を抽出し、その要因数に応じて無呼吸症候群を判定する。
【０００３】
【特許文献１】特開２００４−２４６８４号公報
【特許文献２】特開２００６−１６７２１１号公報
【発明の開示】
【発明が解決しようとする課題】
【０００４】
特許文献１では呼吸動曲線の周波数変化から無呼吸状態もしくは低呼吸状態を判定する際に、寝姿が仰臥位、腹臥位、または横臥位かを推定する必要があり、その寝姿によって検出精度に影響が出る。また特許文献２では呼吸数と心拍数といびきの特徴的要因から無呼吸を判定するが、構造が複雑でありかつ検出精度に問題がある。さらに、特許文献１では被験者の荷重を検出して在離床を判定し、離床のときは判定を行わない。しかし、被験者が在床であっても覚醒状態の場合が有り、このような場合の検出精度は大幅に低下してしまう恐れがある。
【０００５】
本発明の目的は、上記従来技術の問題点に鑑み、呼吸被験者の覚醒／睡眠を体動から判定するとともに、呼吸動波形および／またはいびき波形の通気度から異常呼吸（無呼吸状態もしくは低呼吸状態）を検出し、簡単かつ高精度に異常呼吸状態を判定する測定装置を提供することにある。
【課題を解決するための手段】
【０００６】
上記目的を達成するための本発明は、無拘束・無意識で被験者の無呼吸乃至低呼吸を判定する睡眠時無呼吸症候群測定装置において、被験者が横たわるエアマットの挙動から生体信号を検出し、前記生体信号である体動波形について、次々と前回検出分との差分を取って差分波形（動き量）を抽出し、該差分波形の標準偏差を所定時間にわたって求め、該標準偏差のトレンドにおける波形の振幅と閾値を比較して覚醒状態または睡眠状態を判定する覚醒／睡眠判定部を設け、前記睡眠状態の場合に、前記生体信号から抽出した体動波形および／またはいびき波形に基づいて、無呼吸乃至低呼吸となる異常呼吸を判定することを特徴とする。
【０００７】
あるいは、無拘束・無意識で被験者の無呼吸乃至低呼吸を判定する睡眠時無呼吸症候群測定装置において、被験者が横たわるエアマットの挙動から生体信号を検出し、前記生体信号である体動波形の差分を取って睡眠状態を判定し、前記睡眠状態の場合に、前記生体信号から呼吸波形を読み出し、その１サイクル分の通気量を推定し、前記通気量が閾値より小さい場合に異常呼吸と判定する通気量推定部を有することを特徴とする。
【０００８】
さらに、無拘束・無意識で被験者の無呼吸乃至低呼吸を判定する睡眠時無呼吸症候群測定装置において、被験者が横たわるエアマットの挙動から生体信号を検出し、前記生体信号である体動波形について、次々と前回検出分との差分を取って差分波形（動き量）を抽出し、該差分波形の標準偏差を所定時間にわたって求め、該標準偏差のトレンドにおける波形の振幅と閾値を比較して覚醒状態または睡眠状態を判定する覚醒／睡眠判定部と、前記生体信号から呼吸波形を読み出し、その１サイクル分の通気量を推定し、通気量が閾値より小さい場合に異常呼吸と判定する通気量推定部と、を有することを特徴とする。
【０００９】
また、前記生体信号からいびき波形を読み出し、その１サイクル分の通気量を推定し、呼吸波形といびき波形の通気量を足し合わせた通気量が閾値より小さい場合に異常呼吸と判定する通気量推定部と、を有することを特徴とする睡眠時無呼吸症候群測定装置。
【００１０】
また、前記１サイクル毎の呼吸波形の包絡線を測定し、該包絡線が紡錘形をなす場合に、前記通気量の推定を行うことを特徴とする。
【発明の効果】
【００１１】
本発明によれば、被験者の体動の変化量（動き量）の標準偏差から睡眠状態を判定するので被験者の覚醒／睡眠状態を簡単に判定できる。特に、標準偏差を取ることで動き量の大きさのばらつき具合を数値として求めることができるので、判定が正確になる。また、呼吸波形のサイクル毎に通気量を推定することで異常呼吸の有無が判定できる。この場合、呼吸動波形の包絡線の形状が紡錘型の場合は異常呼吸である可能性が高いので、この判定を踏まえることで高精度の判定が可能になる。
【発明を実施するための最良の形態】
【００１２】
本発明の実施形態について、以下、図１から図８を参照しながら説明する。図１は本発明の一実施例である睡眠時無呼吸症候群測定装置の全体構成を示している。睡眠時無呼吸症候群測定装置はエアマット１０を有しており、このエアマット１０には生体信号を精度よく検出するために、等間隔に配置された穴１０ａが設けられている。エアマット１０は中空チューブ２０を介してセンサ検出部３０に接続されている。
【００１３】
センサ検出部３０は中空チューブ２０を介して検出する。すなわち、呼吸時に生じる呼吸音「スーハー」やいびきによる「グゴー」という所謂「音」ではなく、呼吸時の胸や腹の動きやいびきによる微振動などを生体情報に関わる音圧信号として検出する。この生体情報から、呼吸信号といびき信号をマイクロフォン３０ａで検出し、体位の動きは圧力センサ３０ｂで検出する。マイクロフォン３０ａの信号は、増幅器（ＯＰＡＭＰ１）４０ａで増幅され音圧信号４０ｂになる。圧力センサ３０ｂの信号は増幅器（ＯＰＡＭＰ２）４０ｈで増幅され圧力信号４０ｉとして出力される。
【００１４】
音圧信号４０ｂはＡＤ変換器４０ｃにより呼吸デジタル信号４０ｅに変換される。他方、音圧信号４０ｂはバンドパスフィルタ４０ｄでフィルタリングされ、いびき信号４０ｆとして出力される。
【００１５】
呼吸信号４０ｅ、いびき信号４０ｆ及び体動信号４０ｉは、コンピュータで形成される異常呼吸判定装置（ＣＰＵ）４０ｇに入力され、ＲＯＭ４０ｊに格納されるプログラムによりデジタル処理される。ＲＯＭ４０ｊのプログラムには、覚醒／睡眠判定部４０ｋ、呼吸波形の包絡形状解析部４０ｌ、通気量推定部４０ｍがあり、バス４０ｓを介してＣＰＵ４０ｇに呼び出される。
【００１６】
ＣＰＵ４０ｇのバス４０ｓには、計算結果を一時記憶するＲＡＭ４０ｎ、アラーム用通信デバイスＳＩＯ４０ｏ、上位コンピュータとデータ交換する通信デバイスＬＡＮ４０ｑが接続されている。上位コンピュータは（図示省略）健康管理センターなどに置かれ、患者データを集中的に管理し、異常の監視が行われる。従って、患者はエアマット１０に横たわっているのみで生体情報に関わる音圧信号が観測され、睡眠中の異常呼吸を自動で検知できるので、患者負担が少なく、かつ高性能な睡眠時無呼吸症候群の測定装置を提供できる。
【００１７】
次に、異常呼吸判定装置４０ｇの各部の処理を説明する。図２は覚醒／睡眠判定処理のフローチャートを、図３は覚醒／睡眠判定処理に応じた波形の変化を示している。図２と図３を用いて、本発明による覚醒／睡眠判定処理を説明する。
【００１８】
まず、圧力センサ３０ｂを介してエアマット１０にかかる圧力を体動波形として読み出し（ｓ１０１）、次々と前回検出分との差分を取って差分波形を抽出する（ｓ１０２）。図３の波形データ(１)に示すように、所定サンプリング時間ごとに体動波形の測定データを観測し、前回分との差分である動き量波形(２)を抽出する。図の例では、０．００５秒ごとに体動等の圧力信号をサンプリングし、そのたびに前回（０．００５秒前）の体動の大きさとの差分を求めている。体動（マット圧）の動き量とは圧力差、すなわち速度の関数である。体動波形の差分を求めることで、バイアス成分を除去でき、閾値による判定が容易になる。
【００１９】
次に差分波形の１分間の標準偏差を（１）式により算出し（ｓ１０３）、これを所定期間、たとえば一晩分のデータについて求める（ｓ１０４）。
標準偏差＝√（（動き量ｋ−ｎ点の動き量平均）２／（ｎ点の動き量平均）） …（１）
ここで、ｋ＝０、１、２、・・・、ｎ−１である。標準偏差をとることで、動き量の大きさのばらつき具合を数値として求めることができる。従って、トレンドを求めるときに、動き量の平均値をプロットするよりも、動き量の標準偏差をプロットする方がより小さな動き量まで検知できる。
【００２０】
波形データ(３)は、１晩分の体動波形の動き量に対し標準偏差トレンドを求めた例で、覚醒と睡眠の違いが明瞭に表れている。
【００２１】
さらに、この標準偏差トレンドに対して移動平均を求めることで（ｓ１０５）、覚醒と睡眠の振幅の違いがより顕著になる。その移動平均された波形データ(４)に対し、閾値を設けることで、振幅より大きい場合は覚醒、振幅より小さい場合は睡眠と判定することができる（ｓ１０６）。
【００２２】
本実施例によれば、睡眠と判定された期間について後述する通気量推定を行い、異常呼吸の有無を判定できる。また、被験者がマットに横たわった状態で覚醒している場合には「覚醒」と判定し、誤って「睡眠」とすることはないので、誤判定を確実に防止できる顕著な効果がある。
【００２３】
なお、移動平均には単純移動平均、修正移動平均、荷重移動平均、それらの組み合わせ等、いずれの方法であっても良い。さらに、１分間の動き量波形(２)から体動の回数を求め、その一晩のトレンドを上記判定方法に加味しても良い。
【００２４】
また、判定の精度は落ちるが、体動動き量トレンド（体動動き量に対し、１分間の振幅平均値を一晩分求めたもの）、あるいは脈動波形や体動波形の振幅平均トレンドに対して閾値を設けても判定は可能である。図１の測定装置では脈動波形は検出していないが、ＡＤ変換器４０ｃを２チャンネル設けることにより、生体情報に関わる音圧信号４０ｂから脈動デジタル信号を検出し、脈動波形の１分間の振幅平均値を一晩分求めればよい。
【００２５】
図４は呼吸波形の正常・異常を判定する通気量推定部のフローチャートを示し、図５は通気量推定部の各処理における波形を示すグラフである。図４と図５を用いて、呼吸波形の正常・異常判定法を説明する。
【００２６】
被験者が息を吸う（吸気）と、胸や腹が膨らむためエアマット１０を押す力が大きくなり、また息を吐（呼気）くと、胸や腹が元の状態に戻るのでエアマット１０を押す力が弱まる。その結果、振幅が周期的に変化する呼吸波形を検知できる。一方、エアマット１０上に横たわった状態で、息をしないで胸や腹を膨張・収縮すると、息をしているときと同様の波形を得ることが知られている。このことから、従来のように波形の振幅の大きさを計測するだけでは、睡眠時無呼吸症候群のように胸や腹が動いても息をしていない症状を検知することはできない。
【００２７】
本発明では「通気」に着目して呼吸波形を観測する。通常の呼吸では十分に通気ができているため、安定した胸や腹の膨張・収縮を行い、呼吸波形の包絡線は平坦形となる。一方、数秒間に渡って低呼吸や無呼吸が続く異常呼吸では、被験者は息苦しいため、息をしようと胸や腹の膨張・収縮を繰り返し、十分な通気ができたところで正常の呼吸に戻る。また、一刻も早く息をしようとするため、周期が変動することもある。このような結果、異常呼吸の時間帯では、呼吸波形の包絡線は紡錘形となる。
【００２８】
通気量推定部４０ｍは、マイクロフォン３０ａを介して呼吸波形を読み出し（ｓ２０１）、呼吸波形の１サイクルを検出する（ｓ２０２）。図５の波形(１)に、呼吸波形の１サイクルを示す。１周期は呼気（山側）から吸気（谷側）の１サイクルである。
【００２９】
次に呼吸波形の各サイクルの包絡線を求め（ｓ２０３）、包絡線が平坦形ならば異常呼吸はないとみて処理を終了する（ｓ２０４）。この詳細については後述する。
【００３０】
一方、包絡線が紡錘形であれば、通気量を推定する（ｓ２０５）。通気量は呼吸波形１サイクルの波形率（横幅／高さ）と面積から求まり、これを全てのサイクルにおいて求める。図５の波形（２）に、呼吸１サイクルの通気量推定用パラメータ算出のしかたを示す。波形率は（２）式、通気量は（３）式から求まる。
【００３１】
波形率＝波幅／波高値 …（２）
通気量＝波形率×積算値（面積） …（３）
波形率は呼吸波形が丸いか尖っているかを判定する。波形率が小であれば波形が尖っていて、異常呼吸の可能性が高い。波形率が大であれば波形は丸みを帯び、異常呼吸の可能性は低い。図５の波形(３)、(４)に、１分間に含まれる呼吸波形全サイクルに対して求めた通気量の遷移を示す。波形(３)が正常呼吸時、波形(４)が異常呼吸時の通気量を示し、枠線の中が異常時呼吸である。
【００３２】
ここで、呼吸波形の包絡線解析部４０ｌの処理を説明する。図６は包絡線解析部のフローチャートを、図７はその解析波形を示している。以下、図６と図７に基づいて包絡線解析処理を説明する。
【００３３】
まず呼吸波形を読み出し（ｓ３０１）、山側のピークを線で連結する（ｓ３０２）。同様に、谷側のピークを線で連結する（ｓ３０３）。上記ｓ３０２、ｓ３０３の処理を任意の期間について実施し、包絡線の作成を行う（ｓ３０４、ｓ３０５）。
【００３４】
図７の呼吸波形データ(１)は１分間の呼吸波形で、山側および谷側の包絡線を求めている。これより山側および谷側の包絡線に対して極大値および極小値の有無を判定し、極大値および極小値を検知する毎に両数値の差を所定の閾値と比較し、包絡線の極大値および極小値の差が閾値より大きい場合は紡錘形、差が閾値より小さい場合は平坦形と判定する（ｓ３０６、ｓ３０７）。図７の波形(２)は平坦形であり、波形(３)、(４)は紡錘形の例で、右上がり矢印と右下がり矢印が連続する場合、紡錘形と判定される。このように、山側の包絡線または谷側の包絡線が紡錘形と判定される場合は、異常呼吸の可能性が高いので、包絡線が紡錘形であるか否かを判定することで、異常呼吸判定の精度を向上できる。
【００３５】
再び図４に戻る。次に、推定した通気量と閾値の比較から、閾値を上回れば正常呼吸、閾値を下回れば異常呼吸と判定する（ｓ２１０）。図５の波形(４)は異常呼吸が任意時間（ここでは１０秒）続いたとき、四角の枠で囲んで表示し、アラーム用通信デバイス４０ｏやＬＡＮ４０ｑに出力する。
【００３６】
上記の説明では、呼吸の通気量のみを推定している。しかし、バンドパスフィルタ４０ｄを介していびき波形を読み出し（ｓ２０６）、いびき波形の１波分を検出し（ｓ２０７）、そのサイクル毎に通気量を推定し（ｓ２０８）、その通気量を呼吸の通気量と足し合わせ（ｓ２０９）、この加算した通気量を閾値と比較するようにしても良い。
【００３７】
図８はいびき波形を示すグラフである。波形(１)はいびき波形の１サイクルを示す。いびき波形には周期性が無いため、便宜上、呼吸１サイクル中のいびきを１波として考える。波形(２)はいびき１波の通気量推定用パラメータ算出のしかたを示す。推定式は呼吸波形のものと同様で、上記（２）、（３）式に示した通りである。ただし、いびきの波形は１波の中に複数のピークが存在するため、波高には平均波高値を用いる。なお、波形(１)のように、呼吸の開始時刻といびきの開始時刻は必ずしも一致しない。そこで、位相成分を（４）式により求め、通気量を（５）式で求めるようにしても良い。
位相＝（（いびき開始時刻−呼吸開始時刻）／呼吸周期）×π …（４）
通気量＝波形率×積算値×位相 …（５）
【図面の簡単な説明】
【００３８】
【図１】本発明の一実施例による睡眠時無呼吸症候群測定装置の構成図。
【図２】一実施例による覚醒／睡眠を判定するフローチャート。
【図３】覚醒／睡眠判定処理における各処理の波形データを示すグラフ。
【図４】一実施例による通気量を推定するフローチャート。
【図５】通気量推定処理における各処理の波形データを示すグラフ。
【図６】一実施例による呼吸波形の包絡形状を解析するフローチャート。
【図７】呼吸包絡線の形状解析処理における各処理の波形データを示すグラフ。
【図８】いびき波形の１波と通気量推定用のパラメータを示すグラフ。
【符号の説明】
【００３９】
１０…エアマット、２０…チューブ、３０ａ…マイクロフォン、３０ｂ…圧力センサ、４０ａ，４０ｈ…オペアンプ、４０ｃ…ＡＤ変換器、４０ｄ…バンドパスフィルタ、４０ｇ…異常呼吸判定装置（ＣＰＵ）、４０ｊ…ＲＯＭ、４０ｋ…覚醒／睡眠判定部、４０ｌ…呼吸波形の包絡形状解析部、４０ｍ…通器量推定部、４０ｓ…バス、４０ｎ…ＲＡＭ、４０ｏ…アラーム用通信デバイス（ＳＩＯ）、４０ｑ…通信デバイスＬＡＮ。
【特許請求の範囲】
【請求項１】
無拘束・無意識で被験者の無呼吸乃至低呼吸を判定する睡眠時無呼吸症候群測定装置において、
被験者が横たわるエアマットの挙動から生体信号を検出し、
前記生体信号である体動波形について、次々と前回検出分との差分を取って差分波形（動き量）を抽出し、該差分波形の標準偏差を所定時間にわたって求め、該標準偏差のトレンドにおける波形の振幅と閾値を比較して覚醒状態または睡眠状態を判定する覚醒／睡眠判定部を設け、
前記睡眠状態の場合に、前記生体信号から抽出した呼吸波形および／またはいびき波形に基づいて、無呼吸乃至低呼吸となる異常呼吸を判定することを特徴とする睡眠時無呼吸症候群測定装置。
【請求項２】
無拘束・無意識で被験者の無呼吸乃至低呼吸を判定する睡眠時無呼吸症候群測定装置において、
被験者が横たわるエアマットの挙動から生体信号を検出し、
前記生体信号である体動波形の差分を取って睡眠状態を判定し、
前記睡眠状態の場合に、前記生体信号から呼吸波形を読み出し、その１サイクル分の通気量を推定し、前記通気量が閾値より小さい場合に異常呼吸と判定する通気量推定部を有することを特徴とする睡眠時無呼吸症候群測定装置。
【請求項３】
無拘束・無意識で被験者の無呼吸乃至低呼吸を判定する睡眠時無呼吸症候群測定装置において、
被験者が横たわるエアマットの挙動から生体信号を検出し、
前記生体信号である体動波形について、次々と前回検出分との差分を取って差分波形（動き量）を抽出し、該差分波形の標準偏差を所定時間にわたって求め、該標準偏差のトレンドにおける波形の振幅と閾値を比較して覚醒状態または睡眠状態を判定する覚醒／睡眠判定部と、
前記生体信号から呼吸波形を読み出し、その１サイクル分の通気量を推定し、通気量が閾値より小さい場合に異常呼吸と判定する通気量推定部と、を有することを特徴とする睡眠時無呼吸症候群測定装置。
【請求項４】
無拘束・無意識で被験者の無呼吸乃至低呼吸を判定する睡眠時無呼吸症候群測定装置において、
被験者が横たわるエアマットの挙動から生体信号を検出し、
前記生体信号である体動波形について、次々と前回検出分との差分を取って差分波形（動き量）を抽出し、該差分波形の標準偏差を所定時間にわたって求め、該標準偏差のトレンドにおける波形の振幅と閾値を比較して覚醒状態または睡眠状態を判定する覚醒／睡眠判定部と、
前記生体信号から呼吸波形を読み出し、その１サイクル分の通気量を推定し、前記生体信号からいびき波形を読み出し、その１サイクル分の通気量を推定し、両者の通気量を足し合わせた通気量が閾値より小さい場合に異常呼吸と判定する通気量推定部と、を有することを特徴とする睡眠時無呼吸症候群測定装置。
【請求項５】
請求項３または４において、前記１サイクル毎の呼吸波形の包絡線を測定し、該包絡線が紡錘形をなす場合に、前記通気量の推定を行うことを特徴とする睡眠時無呼吸症候群測定装置。
【請求項６】
請求項１、３乃至５のいずれかにおいて、前記標準偏差のトレンドについて移動平均を求めた後、該移動平均の振幅と閾値を比較して覚醒状態または睡眠状態を判定することを特徴とする睡眠時無呼吸症候群測定装置。
【請求項７】
請求項３乃至６のいずれかにおいて、前記通気量は、前記１サイクルの波形の波形率（波幅／波高値）×積算値から算出することを特徴とする睡眠時無呼吸症候群測定装置。
【請求項１】
無拘束・無意識で被験者の無呼吸乃至低呼吸を判定する睡眠時無呼吸症候群測定装置において、
被験者が横たわるエアマットの挙動から生体信号を検出し、
前記生体信号である体動波形について、次々と前回検出分との差分を取って差分波形（動き量）を抽出し、該差分波形の標準偏差を所定時間にわたって求め、該標準偏差のトレンドにおける波形の振幅と閾値を比較して覚醒状態または睡眠状態を判定する覚醒／睡眠判定部を設け、
前記睡眠状態の場合に、前記生体信号から抽出した呼吸波形および／またはいびき波形に基づいて、無呼吸乃至低呼吸となる異常呼吸を判定することを特徴とする睡眠時無呼吸症候群測定装置。
【請求項２】
無拘束・無意識で被験者の無呼吸乃至低呼吸を判定する睡眠時無呼吸症候群測定装置において、
被験者が横たわるエアマットの挙動から生体信号を検出し、
前記生体信号である体動波形の差分を取って睡眠状態を判定し、
前記睡眠状態の場合に、前記生体信号から呼吸波形を読み出し、その１サイクル分の通気量を推定し、前記通気量が閾値より小さい場合に異常呼吸と判定する通気量推定部を有することを特徴とする睡眠時無呼吸症候群測定装置。
【請求項３】
無拘束・無意識で被験者の無呼吸乃至低呼吸を判定する睡眠時無呼吸症候群測定装置において、
被験者が横たわるエアマットの挙動から生体信号を検出し、
前記生体信号である体動波形について、次々と前回検出分との差分を取って差分波形（動き量）を抽出し、該差分波形の標準偏差を所定時間にわたって求め、該標準偏差のトレンドにおける波形の振幅と閾値を比較して覚醒状態または睡眠状態を判定する覚醒／睡眠判定部と、
前記生体信号から呼吸波形を読み出し、その１サイクル分の通気量を推定し、通気量が閾値より小さい場合に異常呼吸と判定する通気量推定部と、を有することを特徴とする睡眠時無呼吸症候群測定装置。
【請求項４】
無拘束・無意識で被験者の無呼吸乃至低呼吸を判定する睡眠時無呼吸症候群測定装置において、
被験者が横たわるエアマットの挙動から生体信号を検出し、
前記生体信号である体動波形について、次々と前回検出分との差分を取って差分波形（動き量）を抽出し、該差分波形の標準偏差を所定時間にわたって求め、該標準偏差のトレンドにおける波形の振幅と閾値を比較して覚醒状態または睡眠状態を判定する覚醒／睡眠判定部と、
前記生体信号から呼吸波形を読み出し、その１サイクル分の通気量を推定し、前記生体信号からいびき波形を読み出し、その１サイクル分の通気量を推定し、両者の通気量を足し合わせた通気量が閾値より小さい場合に異常呼吸と判定する通気量推定部と、を有することを特徴とする睡眠時無呼吸症候群測定装置。
【請求項５】
請求項３または４において、前記１サイクル毎の呼吸波形の包絡線を測定し、該包絡線が紡錘形をなす場合に、前記通気量の推定を行うことを特徴とする睡眠時無呼吸症候群測定装置。
【請求項６】
請求項１、３乃至５のいずれかにおいて、前記標準偏差のトレンドについて移動平均を求めた後、該移動平均の振幅と閾値を比較して覚醒状態または睡眠状態を判定することを特徴とする睡眠時無呼吸症候群測定装置。
【請求項７】
請求項３乃至６のいずれかにおいて、前記通気量は、前記１サイクルの波形の波形率（波幅／波高値）×積算値から算出することを特徴とする睡眠時無呼吸症候群測定装置。
【図１】
【図２】
【図３】
【図４】
【図５】
【図６】
【図７】
【図８】
【図２】
【図３】
【図４】
【図５】
【図６】
【図７】
【図８】
【公開番号】特開２００９−２８４２３（Ｐ２００９−２８４２３Ａ）
【公開日】平成２１年２月１２日（２００９．２．１２）
【国際特許分類】
【出願番号】特願２００７−１９７５２２（Ｐ２００７−１９７５２２）
【出願日】平成１９年７月３０日（２００７．７．３０）
【出願人】（０００２３３０４４）株式会社日立エンジニアリング・アンド・サービス (276)
【Ｆターム（参考）】
【公開日】平成２１年２月１２日（２００９．２．１２）
【国際特許分類】
【出願日】平成１９年７月３０日（２００７．７．３０）
【出願人】（０００２３３０４４）株式会社日立エンジニアリング・アンド・サービス (276)
【Ｆターム（参考）】
[ Back to top ]