(57)【要約】【課題】CVHRの特徴を拍動間隔から簡単な演算で検出して無呼吸状態の有無を把握できる無呼吸状態検出システムの提供。【解決手段】対象者の心臓の拍動間隔を測定する拍動間隔測定手段と、短時間平均拍動間隔を算出する短時間平均拍動間隔演算手段と、長時間平均拍動間隔を算出する長時間平均拍動間隔演算手段と、短時間平均拍動間隔と長時間平均拍動間隔との差を演算する拍動間隔谷演算手段と、短時間平均拍動間隔と長時間平均拍動間隔との差の時系列のうち、所定の谷深さ閾値Ｖよりも小さい値が存在している箇所を１個の拍動間隔谷と判定する拍動間隔谷判定手段と、所定の判定時間内の拍動間隔谷の数を積算する所定判定時間内谷数積算手段と、拍動間隔谷の数を所定の異常谷数閾値Ｐと比較し、拍動間隔谷の数が異常谷数閾値Ｐ以上ならば無呼吸状態であると判定する谷数比較手段とを備える。【選択図】図１
2020092907.tif 000002
【特許請求の範囲】
【請求項１】
対象者の無呼吸状態を検出する無呼吸状態検出システムであって、対象者の心臓の拍動間隔を測定する拍動間隔測定手段と、前記拍動間隔測定手段により測定した前記拍動間隔から、呼吸の周期程度の所定の時間Ｘで平均した短時間平均拍動間隔を算出する短時間平均拍動間隔演算手段と、前記時間Ｘよりも長い所定の時間Ｙで平均した長時間平均拍動間隔を算出する長時間平均拍動間隔演算手段と、前記短時間平均拍動間隔と前記長時間平均拍動間隔との差を演算する拍動間隔谷演算手段と、前記拍動間隔谷演算手段で得られた前記短時間平均拍動間隔と前記長時間平均拍動間隔との差の時系列のうち、所定の谷深さ閾値Ｖよりも小さい値が存在している箇所を１個の拍動間隔谷と判定する拍動間隔谷判定手段と、所定の判定時間内に前記拍動間隔谷判定手段で判定された拍動間隔谷の数を積算する所定判定時間内谷数積算手段と、前記所定判定時間内谷数積算手段で積算した拍動間隔谷の数を所定の異常谷数閾値Ｐと比較し、前記拍動間隔谷の数が前記異常谷数閾値Ｐ以上ならば無呼吸状態であると判定する谷数比較手段とを備えたことを特徴とする無呼吸状態検出システム。
【請求項２】
請求項１記載の無呼吸状態検出システムにおいて、前記拍動間隔谷判定手段は、前記短時間平均拍動間隔と前記長時間平均拍動間隔との差の時系列のうち、前記谷深さ閾値Ｖよりも小さい値が、所定の谷判定閾値Ｑ以上連続して存在している箇所を１個の拍動間隔谷と判定するように構成されていることを特徴とする無呼吸状態検出システム。
【請求項３】
請求項１，２いずれか１項に記載の無呼吸状態検出システムにおいて、対象者の姿勢を測定する姿勢測定手段と、前記姿勢測定手段から得られた姿勢から立位か伏臥位か側臥位かを判定する体位判定手段と、前記拍動間隔測定手段により測定した前記拍動間隔の時系列から、前記体位判定手段によって得られた体位情報に基づき、立位，伏臥位若しくは側臥位のいずれかの体位のときの拍動間隔を除外し、当該拍動間隔が除外された前記拍動間隔の時系列を前記短時間平均拍動間隔演算手段及び前記長時間平均拍動間隔演算手段に転送する所定体位拍動間隔除外手段とを備えたことを特徴とする無呼吸状態検出システム。
【請求項４】
請求項１，２いずれか１項に記載の無呼吸状態検出システムにおいて、対象者の姿勢を測定する姿勢測定手段と、前記姿勢測定手段から得られた姿勢の時系列から体動を検出する体動検出手段と、前記拍動間隔測定手段により測定した前記拍動間隔の時系列から、前記体動検出手段によって検出された体動が発生した時刻から所定の前後の時間の拍動間隔を除外し、当該拍動間隔が除外された前記拍動間隔の時系列を前記短時間平均拍動間隔演算手段及び前記長時間平均拍動間隔演算手段に転送する体動時拍動間隔除外手段とを備えたことを特徴とする無呼吸状態検出システム。
【請求項５】
請求項１，２いずれか１項に記載の無呼吸状態検出システムにおいて、対象者の姿勢を測定する姿勢測定手段と、
前記姿勢測定手段から得られた姿勢から立位か伏臥位か側臥位かを判定する体位判定手段と、前記拍動間隔測定手段により測定した前記拍動間隔の時系列から、前記体位判定手段によって得られた体位情報に基づき、立位，伏臥位若しくは側臥位のいずれかの体位のときの拍動間隔を除外し、当該拍動間隔が除外された前記拍動間隔の時系列を前記短時間平均拍動間隔演算手段及び前記長時間平均拍動間隔演算手段に転送する所定体位拍動間隔除外手段と、
前記姿勢測定手段から得られた姿勢の時系列から体動を検出する体動検出手段と、前記拍動間隔測定手段により測定した前記拍動間隔の時系列から、前記体動検出手段によって検出された体動が発生した時刻から所定の前後の時間の拍動間隔を除外し、当該拍動間隔が除外された前記拍動間隔の時系列を前記短時間平均拍動間隔演算手段及び前記長時間平均拍動間隔演算手段に転送する体動時拍動間隔除外手段とを備えたことを特徴とする無呼吸状態検出システム。
【請求項６】
請求項１〜５いずれか１項に記載の無呼吸状態検出システムにおいて、前記谷深さ閾値Ｖは、前記長時間平均拍動間隔演算手段から算出された長時間平均拍動間隔の関数であることを特徴とする無呼吸状態検出システム。
【請求項７】
対象者の無呼吸状態を解消して安眠を提供する安眠提供システムであって、対象者の心臓の拍動間隔を測定する拍動間隔測定手段と、前記拍動間隔測定手段により測定した前記拍動間隔から、無呼吸状態であることを判定する解析器と、音または振動を発生し対象者をわずかに覚醒させることによって無呼吸状態を解消し眠りを改善する刺激手段とを備えたことを特徴とする安眠提供システム。
【請求項８】
対象者の無呼吸状態を解消して安眠を提供する安眠提供システムであって、対象者の心臓の拍動間隔を測定する拍動間隔測定手段と、前記拍動間隔測定手段により測定した前記拍動間隔から、呼吸の周期程度の所定の時間Ｘで平均した短時間平均拍動間隔を算出する短時間平均拍動間隔演算手段と、前記時間Ｘよりも長い所定の時間Ｙで平均した長時間平均拍動間隔を算出する長時間平均拍動間隔演算手段と、前記短時間平均拍動間隔と前記長時間平均拍動間隔との差を演算する拍動間隔谷演算手段と、前記拍動間隔谷演算手段で得られた前記短時間平均拍動間隔と前記長時間平均拍動間隔との差の時系列のうち、所定の谷深さ閾値Ｖよりも小さい値が存在している箇所を１個の拍動間隔谷と判定する拍動間隔谷判定手段と、所定の判定時間内に前記拍動間隔谷判定手段で判定された拍動間隔谷の数を積算する所定判定時間内谷数積算手段と、前記所定判定時間内谷数積算手段で積算した拍動間隔谷の数を所定の異常谷数閾値Ｐと比較し、前記拍動間隔谷の数が前記異常谷数閾値Ｐ以上ならば無呼吸状態であると判定する谷数比較手段と、前記谷数比較手段が無呼吸状態であると判定したときに、音または振動を発生し対象者をわずかに覚醒させることによって無呼吸状態を解消し眠りを改善する刺激手段とを備えたことを特徴とする安眠提供システム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、無呼吸状態検出システム及び安眠提供システムに関するものである。
【背景技術】
【０００２】
睡眠時無呼吸症候群(SAS)の潜在患者数は、治療が必要なものに限定しても256万人いると言われている。SAS患者の交通事故発生率は、健常者の約７倍と有意に高く、また薬剤抵抗性高血圧、心不全、心房細動、高血圧、冠動脈疾患、糖尿病などの様々な生活習慣病が合併する。無呼吸状態を改善する装置で逐次無呼吸状態を通知し、眠りを改善すれば交通事故を避けることができるし、職場の健康診断などによる検診でSASを早期発見できれば、成人病などによる経済損失を減少できる。
【０００３】
ところで、SASには、閉塞型睡眠時無呼吸症候群(OSAS)および中枢性睡眠時無呼吸症候群があるが、いずれの場合も拍動間隔に周期的な変動があり、CVHR(Cyclic Variation of Heart Rate)と呼ばれている。このCVHRは25から120秒の周期性があり、心拍変動のパワースペクトラム密度を0.008から0.04Hzまで積分して得られるVLFの帯域と一致している。実際にこの特徴に基づいてSASを検出する方法が提案されている（特許文献１参照）。
【０００４】
このようなCVHRの特徴を利用してSASを検出するには長時間心電図を用いる方法があるが、長時間心電図を取得するには高価なホルター心電計が必要であり、心電図から拍動間隔を導出しさらにVLFを計算するには扱うデータ量が膨大で、多くの人が同時に受診する健康診断では用いることができない。
【０００５】
また、人は睡眠中寝返りが生じ、このようなときに心電図にアーチファクトが混入するため、VLFを正しく計算できないことがある。さらに、OSASではうつぶせ（伏臥位）や横向き（側臥位）で寝ているとき無呼吸状態にならない。
【０００６】
またSASは、前述の通り心房細動も合併することがあるので、健康診断においてこれを同時に検出することも期待される。
【先行技術文献】
【特許文献】
【０００７】
【特許文献１】特開２００５−１６０６５０号公報
【特許文献２】特許第６１５０８２５号公報
【発明の概要】
【発明が解決しようとする課題】
【０００８】
本発明は、上述のような現状に鑑みなされたもので、CVHRの特徴を拍動間隔から簡単な演算で検出して無呼吸状態の有無を把握でき、例えば、人の姿勢も同時に解析してアーチファクトによる誤検出なく眠りを改善したり、健康診断で用いることのできる（特許文献２等に開示される）心房細動検出システムと併用したりすることも可能な実用的な無呼吸状態検出システム及び安眠提供システムを提供するものである。
【課題を解決するための手段】
【０００９】
添付図面を参照して本発明の要旨を説明する。
【００１０】
対象者の無呼吸状態を検出する無呼吸状態検出システム１であって、対象者の心臓の拍動間隔を測定する拍動間隔測定手段４と、前記拍動間隔測定手段４により測定した前記拍動間隔から、呼吸の周期程度の所定の時間Ｘで平均した短時間平均拍動間隔を算出する短時間平均拍動間隔演算手段８と、前記時間Ｘよりも長い所定の時間Ｙで平均した長時間平均拍動間隔を算出する長時間平均拍動間隔演算手段９と、前記短時間平均拍動間隔と前記長時間平均拍動間隔との差を演算する拍動間隔谷演算手段10と、前記拍動間隔谷演算手段10で得られた前記短時間平均拍動間隔と前記長時間平均拍動間隔との差の時系列のうち、所定の谷深さ閾値Ｖよりも小さい値が存在している箇所を１個の拍動間隔谷と判定する拍動間隔谷判定手段11と、所定の判定時間内に前記拍動間隔谷判定手段11で判定された拍動間隔谷の数を積算する所定判定時間内谷数積算手段14と、前記所定判定時間内谷数積算手段14で積算した拍動間隔谷の数を所定の異常谷数閾値Ｐと比較し、前記拍動間隔谷の数が前記異常谷数閾値Ｐ以上ならば無呼吸状態であると判定する谷数比較手段16とを備えたことを特徴とする無呼吸状態検出システムに係るものである。
【００１１】
また、請求項１記載の無呼吸状態検出システムにおいて、前記拍動間隔谷判定手段11は、前記短時間平均拍動間隔と前記長時間平均拍動間隔との差の時系列のうち、前記谷深さ閾値Ｖよりも小さい値が、所定の谷判定閾値Ｑ以上連続して存在している箇所を１個の拍動間隔谷と判定するように構成されていることを特徴とする無呼吸状態検出システムに係るものである。
【００１２】
また、請求項１，２いずれか１項に記載の無呼吸状態検出システムにおいて、対象者の姿勢を測定する姿勢測定手段19と、前記姿勢測定手段19から得られた姿勢から立位か伏臥位か側臥位かを判定する体位判定手段23と、前記拍動間隔測定手段４により測定した前記拍動間隔の時系列から、前記体位判定手段23によって得られた体位情報に基づき、立位，伏臥位若しくは側臥位のいずれかの体位のときの拍動間隔を除外し、当該拍動間隔が除外された前記拍動間隔の時系列を前記短時間平均拍動間隔演算手段８及び前記長時間平均拍動間隔演算手段９に転送する所定体位拍動間隔除外手段24とを備えたことを特徴とする無呼吸状態検出システムに係るものである。
【００１３】
また、請求項１，２いずれか１項に記載の無呼吸状態検出システムにおいて、対象者の姿勢を測定する姿勢測定手段19と、前記姿勢測定手段19から得られた姿勢の時系列から体動を検出する体動検出手段25と、前記拍動間隔測定手段４により測定した前記拍動間隔の時系列から、前記体動検出手段25によって検出された体動が発生した時刻から所定の前後の時間の拍動間隔を除外し、当該拍動間隔が除外された前記拍動間隔の時系列を前記短時間平均拍動間隔演算手段８及び前記長時間平均拍動間隔演算手段９に転送する体動時拍動間隔除外手段26とを備えたことを特徴とする無呼吸状態検出システムに係るものである。
【００１４】
また、請求項１，２いずれか１項に記載の無呼吸状態検出システムにおいて、対象者の姿勢を測定する姿勢測定手段19と、
前記姿勢測定手段19から得られた姿勢から立位か伏臥位か側臥位かを判定する体位判定手段23と、前記拍動間隔測定手段４により測定した前記拍動間隔の時系列から、前記体位判定手段23によって得られた体位情報に基づき、立位，伏臥位若しくは側臥位のいずれかの体位のときの拍動間隔を除外し、当該拍動間隔が除外された前記拍動間隔の時系列を前記短時間平均拍動間隔演算手段８及び前記長時間平均拍動間隔演算手段９に転送する所定体位拍動間隔除外手段24と、
前記姿勢測定手段19から得られた姿勢の時系列から体動を検出する体動検出手段25と、前記拍動間隔測定手段４により測定した前記拍動間隔の時系列から、前記体動検出手段25によって検出された体動が発生した時刻から所定の前後の時間の拍動間隔を除外し、当該拍動間隔が除外された前記拍動間隔の時系列を前記短時間平均拍動間隔演算手段８及び前記長時間平均拍動間隔演算手段９に転送する体動時拍動間隔除外手段26とを備えたことを特徴とする無呼吸状態検出システムに係るものである。
【００１５】
また、請求項１〜５いずれか１項に記載の無呼吸状態検出システムにおいて、前記谷深さ閾値Ｖは、前記長時間平均拍動間隔演算手段９から算出された長時間平均拍動間隔の関数であることを特徴とする無呼吸状態検出システムに係るものである。
【００１６】
また、対象者の無呼吸状態を解消して安眠を提供する安眠提供システムであって、対象者の心臓の拍動間隔を測定する拍動間隔測定手段４と、前記拍動間隔測定手段４により測定した前記拍動間隔から、無呼吸状態であることを判定する解析器３と、音または振動を発生し対象者をわずかに覚醒させることによって無呼吸状態を解消し眠りを改善する刺激手段とを備えたことを特徴とする安眠提供システムに係るものである。
【００１７】
また、対象者の無呼吸状態を解消して安眠を提供する安眠提供システムであって、対象者の心臓の拍動間隔を測定する拍動間隔測定手段４と、前記拍動間隔測定手段４により測定した前記拍動間隔から、呼吸の周期程度の所定の時間Ｘで平均した短時間平均拍動間隔を算出する短時間平均拍動間隔演算手段８と、前記時間Ｘよりも長い所定の時間Ｙで平均した長時間平均拍動間隔を算出する長時間平均拍動間隔演算手段９と、前記短時間平均拍動間隔と前記長時間平均拍動間隔との差を演算する拍動間隔谷演算手段10と、前記拍動間隔谷演算手段10で得られた前記短時間平均拍動間隔と前記長時間平均拍動間隔との差の時系列のうち、所定の谷深さ閾値Ｖよりも小さい値が存在している箇所を１個の拍動間隔谷と判定する拍動間隔谷判定手段11と、所定の判定時間内に前記拍動間隔谷判定手段11で判定された拍動間隔谷の数を積算する所定判定時間内谷数積算手段14と、前記所定判定時間内谷数積算手段14で積算した拍動間隔谷の数を所定の異常谷数閾値Ｐと比較し、前記拍動間隔谷の数が前記異常谷数閾値Ｐ以上ならば無呼吸状態であると判定する谷数比較手段16と、前記谷数比較手段16が無呼吸状態であると判定したときに、音または振動を発生し対象者をわずかに覚醒させることによって無呼吸状態を解消し眠りを改善する刺激手段とを備えたことを特徴とする安眠提供システムに係るものである。
【発明の効果】
【００１８】
本発明は上述のように構成したから、CVHRの特徴を拍動間隔から簡単な演算で検出して無呼吸状態の有無を把握できる実用的な無呼吸状態検出システム及び安眠提供システムとなる。
【図面の簡単な説明】
【００１９】
【図１】本実施例の構成概略説明図である。
【図２】別例１の構成概略説明図である。
【図３】別例２の構成概略説明図である。
【図４】別例３の構成概略説明図である。
【図５】CVHRが生じている拍動間隔の時系列を示すグラフである。
【図６】無呼吸状態の検出方法を説明するためのグラフである。
【図７】体動によって生じる拍動間隔の変化を説明するためのグラフである。
【発明を実施するための形態】
【００２０】
好適と考える本発明の実施形態を、図面に基づいて本発明の作用を示して簡単に説明する。
【００２１】
睡眠時無呼吸症候群(SAS)のために無呼吸状態になると、CVHRと呼ばれる拍動間隔が周期的に大きく増大減少する変動が生じる（変化に山と谷が生じる。）。図５はCVHRが生じている拍動間隔の時系列を示すもので、横軸は時刻、縦軸は拍動間隔であり、CVHRの発生箇所を括弧で示している。従って、拍動間隔の周期的な変動を探索し（CVHRが生じていることが確認できれば）、SASによる無呼吸状態を検出できる。具体的には、拍動間隔の時系列中に所定深さの谷（深い谷）が周期的に出現するかどうか調べれば良い。
【００２２】
本発明者等は、短時間平均拍動間隔と長時間平均拍動間隔との差の時系列中で、所定の谷深さ閾値Ｖよりも小さい値が存在している箇所（好ましくは谷深さ閾値Ｖよりも小さい値が所定の谷判定閾値Ｑ以上連続している箇所）を１個の拍動間隔谷とし、所定の判定時間内の拍動間隔谷の数が所定の異常谷数閾値Ｐ以上ある場合にCVHRが生じていると判断することが可能であることを確認した。
【００２３】
従って、本発明によれば、CVHRが生じていること、即ち、無呼吸状態であることを拍動間隔を用いた簡単な演算で検出可能となる。
【実施例】
【００２４】
本発明の具体的な実施例について図面に基づいて説明する。
【００２５】
本実施例は、対象者の無呼吸状態を検出する無呼吸状態検出システム１であって、対象者の心臓の拍動間隔を測定する拍動間隔測定手段４と、前記拍動間隔測定手段４により測定した前記拍動間隔から、呼吸の周期程度の所定の時間Ｘで平均した短時間平均拍動間隔を算出する短時間平均拍動間隔演算手段８と、前記時間Ｘよりも長い所定の時間Ｙで平均した長時間平均拍動間隔を算出する長時間平均拍動間隔演算手段９と、前記短時間平均拍動間隔と前記長時間平均拍動間隔との差を演算する拍動間隔谷演算手段10と、前記拍動間隔谷演算手段10で得られた前記短時間平均拍動間隔と前記長時間平均拍動間隔との差の時系列のうち、所定の谷深さ閾値Ｖよりも小さい値が存在している箇所を１個の拍動間隔谷と判定する拍動間隔谷判定手段11と、所定の判定時間内に前記拍動間隔谷判定手段11で判定された拍動間隔谷の数を積算する所定判定時間内谷数積算手段14と、前記所定判定時間内谷数積算手段14で積算した拍動間隔谷の数を所定の異常谷数閾値Ｐと比較し、前記拍動間隔谷の数が前記異常谷数閾値Ｐ以上ならば無呼吸状態であると判定する谷数比較手段16とを備えたものである。
【００２６】
具体的には、本実施例は、図１に図示したように、拍動間隔測定手段４が設けられた拍動間隔測定用のセンサ２と、短時間平均拍動間隔演算手段８と長時間平均拍動間隔演算手段９と拍動間隔谷演算手段10と拍動間隔谷判定手段11と所定判定時間内谷数積算手段14と谷数比較手段16とが設けられた解析器３とで構成されている。
【００２７】
各部を具体的に説明する。
【００２８】
センサ２には、拍動間隔測定手段４と、この拍動間隔測定手段４によって測定された拍動間隔データを一時的に保存する拍動間隔保存手段５と、この一時的に保存した拍動間隔データを解析器３の拍動間隔受信手段７に送信する拍動間隔送信手段６とが設けられている。
【００２９】
拍動間隔測定手段４は、例えばマイコン等を用いて電極から得られた電圧の変化をもとにした心電図から一つのＲ波とこれに隣り合う他のＲ波との間隔、または、一つのＳ波とこれに隣り合う他のＳ波との間隔から拍動間隔を測定するように構成されている。従って、センサを小型化して電極を介して皮膚に貼り付ければ、センサを衣服の下に隠すことができ、睡眠に支障を来すことなく測定できる。なお、拍動間隔測定手段４は、例えば赤外線の反射光から脈波を測定し、そのピーク間隔などから拍動間隔を測定するように構成してもよい。この場合、耳たぶや手首や腕などにクリップやバンドでセンサを固定するだけで良く、装着しやすいものとなる。
【００３０】
拍動間隔測定手段４によって測定された拍動間隔データは、拍動間隔保存手段５へ逐次転送され一時保存される。検査のためにセンサ２のみを対象者に適用する場合、一時保存する拍動間隔データは半日間程度分となる。この場合、拍動間隔保存手段５は、半導体メモリやテープなどを採用できる。ホルター心電計とは異なり、心電図波形を保存する必要がなく、拍動間隔保存手段５として半導体メモリを用いた場合、非常に小型で低消費電力のメモリを用いることができるので、小型かつ軽量で対象者の負担にならないセンサ２を構成することが可能となる。検査のため、若しくは、無呼吸を改善するための装置を動作させるため、対象者の拍動間隔データからリアルタイムで無呼吸状態か否か判定したい場合は、拍動間隔保存手段５は外部へ拍動間隔データを転送するための一時バッファであるので、例えば、拍動間隔測定手段４のために用いられるマイコン内のランダムアクセスメモリでも良い。この場合は、センサ２をさらに小型かつ軽量化でき、対象者の負担を一層低減できる。
【００３１】
拍動間隔保存手段５で一時保存された拍動間隔データは、拍動間隔送信手段６へ転送される。
【００３２】
拍動間隔送信手段６は、拍動間隔保存手段５から転送された拍動間隔データを受け取り、解析器３に設けた拍動間隔受信手段７へ転送する。転送方法としては、電波や光を用いた無線や、USBやRS-232Cなどの有線接続を用いることができる。センサ２と解析器３との接続に無線接続を用いたとき、センサ２と解析器３とを物理的に分離でき、センサ２を取り付けられた対象者の負担が小さくなる。従って、検査のため、若しくは、無呼吸を改善するための装置を動作させるため対象者の拍動間隔データからリアルタイムで無呼吸状態か否か判定したい場合に好適である。また、電話回線やインターネットなどの公衆回線を用いることも可能であり、この場合、対象者に対して遠隔に無呼吸状態を検出することが可能となる。また、センサ２と解析器３との接続に有線接続を用いた場合は、拍動間隔データを確実かつ高速に転送できるので、拍動間隔保存手段５に保存された半日間程度分の拍動間隔データを、拍動間隔送信手段６を介して拍動間隔受信手段７へ転送するのに好適である。
【００３３】
解析器３には、拍動間隔受信手段７と、短時間平均拍動間隔演算手段８と長時間平均拍動間隔演算手段９と拍動間隔谷演算手段10と拍動間隔谷判定手段11と所定拍動間隔谷深さ保存手段12と所定谷内拍動間隔数保存手段13と所定判定時間内谷数積算手段14と谷数比較手段16と所定谷数保存手段17とが設けられている。解析器３は一連の計算、比較、表示を行う電子計算機や計測器であり、上記各手段を備えた専用の機器、パーソナルコンピュータ、タブレット型のコンピュータ、スマートフォン、携帯電話もしくはインターネット上のサーバなどを採用することができる。
【００３４】
拍動間隔受信手段７は、センサ２の拍動間隔送信手段６から拍動間隔データを受信し、短時間平均拍動間隔演算手段８と長時間平均拍動間隔演算手段９とへ転送する。
【００３５】
短時間平均拍動間隔演算手段８は、呼吸の周期程度の時間である所定の短時間（時間Ｘ）で平均した短時間平均拍動間隔SRを算出する。
【００３６】
長時間平均拍動間隔演算手段９は、呼吸の周期程度の時間よりも十分長い所定の長時間（時間Ｙ）で平均した長時間平均拍動間隔LRを算出する。
【００３７】
短時間平均拍動間隔演算手段８及び長時間平均拍動間隔演算手段９のいずれの手段も算出したデータを拍動間隔谷演算手段10へ転送する。
【００３８】
また、詳細は後述するが、本実施例では、時間Ｘは例えば４秒から８秒程度に設定する。時間Ｙを例えばCVHRの周期である25秒から120秒程度に設定する。従って、例えば拍動間隔測定手段４で拍動が検出された毎に測定され蓄積された拍動間隔データを用い、演算時点から過去４秒分の平均拍動間隔を短時間平均拍動間隔SRとし、過去25秒分の平均拍動間隔を長時間平均拍動間隔LRとすることができる。また、短時間平均拍動間隔演算手段８及び長時間平均拍動間隔演算手段９の演算間隔は、例えば拍動が検出された毎に測定された拍動間隔と同間隔とすることができる。
【００３９】
拍動間隔谷演算手段10は、短時間平均拍動間隔演算手段８で算出された短時間平均拍動間隔SRと長時間平均拍動間隔演算手段９で算出された長時間平均拍動間隔LRの平均拍動間隔差SR-LRを演算し、その結果を拍動間隔谷判定手段11へ転送する。
【００４０】
拍動間隔谷判定手段11は、拍動間隔谷演算手段10で演算された平均拍動間隔差SR-LR（の時系列）のうち、所定拍動間隔谷深さ保存手段12に保存された所定の谷深さ閾値Ｖよりも小さい値（拍動間隔）が、所定谷内拍動間隔数保存手段13に保存された所定の谷判定閾値Ｑ以上連続しているとき、これを１個の拍動間隔谷と判定し、拍動間隔谷判定結果データを所定判定時間内谷数積算手段14へ転送する。所定拍動間隔谷深さ保存手段12に保存された所定の谷深さ閾値Ｖは解析器３内で保存されているが、この谷深さ閾値Ｖはユーザが設定しても良い。そうすれば、無呼吸状態の検出感度を調整できる。
【００４１】
なお、拍動間隔谷判定手段11は、谷深さ閾値Ｖよりも小さい値が所定数連続する場合に限らず、１つだけ存在した場合にこれを１個の拍動間隔谷と判定するように構成しても良い。
【００４２】
谷深さ閾値Ｖは、拍動間隔の長時間平均拍動間隔LRの関数Ｖ＝Ｖ(LR)としてもよい。関数Ｖ(LR)としては、例えば、αを定数としてＶ(LR)＝αLRやＶ(LR)＝αLR２などを採用できる。所定谷内拍動間隔数保存手段13に保存された所定の拍動間隔数は、後述するように、突発的な拍動間隔谷の検出を防止し、誤検出を防止するために用いられる。従って、対象者の状態によって適宜設定すれば、感度と特異度のバランスのとれた無呼吸状態検出システムを構成できる。
【００４３】
所定判定時間内谷数積算手段14は、所定判定時間内に拍動間隔谷判定結果データに幾つの拍動間隔谷が含まれるかカウントし、その結果を谷数比較手段16へ転送する。
【００４４】
谷数比較手段16は、所定谷数保存手段17に保存された所定の異常谷数閾値Ｐと比較し、拍動間隔谷の数が異常谷数閾値Ｐ以上ならば無呼吸状態であると判定する。異常谷数閾値Ｐを小さくすれば、感度は上がるが誤検出の恐れが生じる。また、異常谷数閾値Ｐを大きくすれば、無呼吸状態における拍動間隔変動周期の短い対象者の無呼吸状態を、より正確に検出することができる。従って、対象者の状態に応じて異常谷数閾値Ｐを適宜設定すれば、無呼吸状態検出システムの誤検出を避けることができる。
【００４５】
谷数比較手段16で判定した無呼吸状態の有無は、無呼吸状態の有無を通知する通知手段18へ転送される。谷数比較手段16と通知手段18は同じ筐体（解析器３）内に配置されていてもよい。この場合、転送方法はプリント配線基板上の配線による。あるいは、谷数比較手段16が含まれる解析器３と、通知手段18とは別体であってもよい。この場合、転送方法としては、電波や光を用いた無線や、USBやRS-232Cなどの有線接続を用いることができる。無線を用いたとき、通知手段18と解析器３を物理的に分離でき、通知手段18のみを枕元に置くことで対象者にとって利便性の高い無呼吸状態検出システムとなる。また、電話回線やインターネットなどの公衆回線を用いることも可能であり、この場合、対象者に対して遠隔に無呼吸状態を検出しつつ、結果のみを対象者へ通知することができるので、対象者は解析器３の配置などの手間がない。有線接続を用いたとき、解析器３から通知手段18へ電力供給もできるので、無呼吸状態を検出した場合にベッドを揺らすなどが可能になり、通知方法の多様性が増す。
【００４６】
通知手段18としては、例えば文字や画像などを表示するディスプレイを用いることができる。この場合、拍動間隔データと無呼吸状態の有無をグラフ表示し、いつ無呼吸状態になったのかをわかりやすく示すことができる。また、通知手段18として、光や音や振動を用いることもできる。この場合、無呼吸状態を直ちに通知できる。
【００４７】
また、通知手段18と共に、若しくは、通知手段18を設けずに、無呼吸状態を検出した際に対象者をわずかに覚醒する程度（睡眠状態を継続できる程度）に、音や振動を発生させたり、枕やベッドをわずかに傾けたりする刺激手段を設けた場合には、無呼吸状態を検出した際に対象者の無呼吸状態を解消し眠りを改善する効果が得られることになり、対象者の無呼吸状態を解消して安眠を提供する安眠提供システムを実現可能となる。
【００４８】
本実施例による無呼吸状態の検出方法について詳述する。
【００４９】
上述したように、SASによる無呼吸状態を検出するには、拍動間隔の時系列中に所定深さの谷（深い谷）が周期的に出現するかどうか調べれば良い。
【００５０】
この方法を図６に従って説明する。「深い谷」を検出するには、谷の深さの測り方が必要である。谷の深さを測るための深さの基準として、拍動間隔の日内変動を考慮し、所定の長時間における拍動間隔の平均を用いる。所定の長時間（時間Ｙ）は、CVHRの周期である25秒から120秒程度が好適である。拍動間隔と、谷の深さを測るための基準とした所定の長時間における拍動間隔の平均との差が、所定の値より小さいひとまとまりを１つの谷とすることができる。しかし、拍動間隔には常に呼吸性不整脈が含まれるため、これに伴う拍動間隔の増減がいちいち谷として識別される恐れがある。そこで、呼吸性不整脈の影響を小さくするため、所定の短時間（時間Ｘ）における拍動間隔の平均を用いる。所定の短時間は、呼吸の周期が４秒程度なので、４秒から８秒が好適である。前述したいずれの拍動間隔の平均も、次のように求めることができる。
【００５１】
つまり、時刻ｔの拍動間隔の平均Ｒバー(ｔ)を求める方法は単純平均でも良く、拍動間隔の時系列をＲ(ｔｉ)、平均時間をＴとして、
【００５２】
【数１】
2020092907.tif 000003
【００５３】
を採用することができる。ここで、インデックスｉは何番目に測定されたかを表し、ｔｉはｉ番目に測定された時刻である。また、ｎ(ｔ,Ｔ)は、時刻ｔ−Ｔから時刻ｔまでに測定された拍動間隔の個数を表す。また、片側ガウス分布を用いた加重平均でも良い。この場合時刻ｔの平均拍動間隔Ｒバー(ｔ)は、
【００５４】
【数２】
2020092907.tif 000004
【００５５】
で表され、時間的により一層滑らかな平均拍動間隔Ｒバー(ｔ)が得られる。図６では、所定の短時間における拍動間隔の平均、短時間平均拍動間隔SRを、Ｔ＝５ｓとした式(1)を使って求め、所定の長時間における拍動間隔の平均、長時間平均拍動間隔LRを、Ｔ＝６０ｓとした式(1)を使って求めたものを示している。
【００５６】
無呼吸によって生じる拍動間隔の「深い谷」は、短時間平均拍動間隔SRと長時間平均拍動間隔LRとの差が、所定の値Ｖ（谷深さ閾値Ｖ）よりも小さい、SR−LR＜Ｖの部分である。ところで、拍動間隔の変動は、心拍数が大きければ（平均拍動間隔が小さければ）小さく、心拍数が小さければ（平均拍動間隔が大きければ）大きい。従って、「深い谷」を判定するための所定の値Ｖは、長時間における拍動間隔の平均LRに応じて変化する、関数Ｖ＝Ｖ(LR)でもよい。関数Ｖ(LR)としては、例えば、αを定数としてＶ(LR)＝αLRやＶ(LR)＝αLR２などを採用できる。
【００５７】
図６では、Ｖ(LR)＝−0.00009LR２として、ＡとＢの２カ所の谷を検出している。さらに、SR−LR＜Ｖによって探索した谷であっても、レム睡眠時によって生じた突発的な谷である可能性がある。これを避けるために、SR−LR＜Ｖに含まれる拍動間隔のデータが、所定の数Ｑ（谷判定閾値Ｑ）以上連続しているときに、無呼吸状態によって生じた「深い谷」の候補と判定する。所定の数Ｑは、突発的な谷を検出しないためであるから、２以上が好適である。図６のＡとＢの箇所は、拍動間隔データ（SR−LRの時系列）が連続で２以上含まれている。最後に、SASの特徴であるCVHRは、拍動間隔が周期的に大きく増大減少することであるので、所定の判定時間内に、所定の数Ｐ（異常谷数閾値Ｐ）以上の「深い谷」を検出したときに無呼吸状態と判定すればよい。所定の判定時間は、CVHRの周期が25秒から120秒程度であることから、120秒以上が好適である。所定の数Ｐは、「深い谷」が繰り返しているか調べたいので、２以上が好適である。図６では、120秒の間にＡとＢの２カ所の「深い谷」があり、無呼吸と判定できる。
【００５８】
図２に本実施例の別例１を示す。
【００５９】
閉塞型睡眠時無呼吸症候群(OSAS)である対象者が、側臥位や伏臥位で睡眠しているときは、無呼吸状態にならない。従って、対象者の姿勢が仰臥位かどうか検知できれば、より一層誤検出の少ない無呼吸状態検出システムを構成できる。
【００６０】
図２に示す別例１は、図１の無呼吸状態検出システム１に、拍動間隔を測定するのと同時に対象者の姿勢を測定する姿勢測定手段19と、姿勢測定手段19によって得られた姿勢から、立位か伏臥位か側臥位かを判定する体位判定手段23と、拍動間隔測定手段４により測定した拍動間隔の時系列から、体位判定手段23によって得られた姿勢状態（体位情報）に基づき、立位，伏臥位若しくは側臥位のいずれか、または、立位，伏臥位及び側臥位すべての状態のとき（仰臥位でないとき）の拍動間隔を除外し、当該拍動間隔が除外された除外後の前記拍動間隔の時系列を短時間平均拍動間隔演算手段８及び長時間平均拍動間隔演算手段９に転送する所定体位拍動間隔除外手段24とを備えたものである。
【００６１】
具体的に、図２に示すセンサ２には、本実施例（図１）のセンサ２に、姿勢測定手段19と姿勢保存手段20と姿勢送信手段21とが付加されている。また、図１の解析器３に、姿勢受信手段22と体位判定手段23と所定体位拍動間隔除外手段24とが付加されている。
【００６２】
センサ２の姿勢測定手段19は、対象者の姿勢を測定できるデバイスであり、例えばMEMSの３軸加速度センサを使用することができる。このようなセンサは軽量、小型、低消費電力なので、被験者の負担は小さいものになる。
【００６３】
姿勢測定手段19によって測定された姿勢データは、姿勢保存手段20へ逐次転送され一時保存される。センサ２を軽量、小型かつ低消費電力とするため、姿勢保存手段20は拍動間隔保存手段５と共通でもよい。
【００６４】
姿勢保存手段20で一時保存された姿勢データは、姿勢送信手段21へ転送される。この転送経路は、センサ２を軽量、小型かつ低消費電力とするため、拍動間隔保存手段５と拍動間隔送信手段６の転送経路と共通でもよい。
【００６５】
姿勢送信手段21は、姿勢保存手段20から転送された姿勢データを受け取り、解析器３に設けた姿勢受信手段22へ転送する。センサ２を軽量、小型かつ低消費電力とするため、姿勢送信手段21と拍動間隔送信手段６とを一体とし、更に、姿勢受信手段22と拍動間隔受信手段７とを一体としても良い。このとき、姿勢送信手段21から姿勢受信手段22への転送経路は、拍動間隔送信手段６から拍動間隔受信手段７への転送経路と共通となる。
【００６６】
姿勢受信手段22で受信した姿勢データは体位判定手段23へ転送される。体位判定手段23は、対象者が仰臥位で寝ているかどうかを判定する。例えば、姿勢測定手段19が重力に対しても感度を持つ３軸加速度センサで構成され、ｘ，ｙ，ｚ軸の出力それぞれを、対象者の右から左、対象者の足から頭、対象者の胸から背中の方向へ一致させ、出力値を(Ａｘ，Ａｙ，Ａｚ)とする。もし、対象者が臥位であれば、対象者の足から頭の方向には重力加速度は０に近いので、
【００６７】
【数３】
2020092907.tif 000005
【００６８】
を満たす。対象者が仰臥位かそれ以外かは、ＡｘとＡｚの合成ベクトルの向く角度から判定できる。合成ベクトルの角度θを
【００６９】
【数４】
2020092907.tif 000006
【００７０】
として、
【００７１】
【数５】
2020092907.tif 000007
【００７２】
ならば仰臥位と判定できる。つまり式(3)と式(5)を満たせば仰臥位である。
【００７３】
体位判定手段23は前述のように、姿勢データが仰臥位かそれ以外かを判定し、所定体位拍動間隔除外手段24へ判定結果を転送する。
【００７４】
所定体位拍動間隔除外手段24は、体位判定手段23から姿勢データが仰臥位かどうかを受けとると共に、拍動間隔受信手段７から拍動間隔データを受け取り、拍動間隔データから対象者の姿勢が仰臥位でない時間の拍動間隔を除外し、除外後の拍動間隔データを短時間平均拍動間隔演算手段８及び長時間平均拍動間隔演算手段９へ転送する。
【００７５】
従って、別例１においては、一層誤検出の少ない無呼吸状態検出システムを実現可能となる。
【００７６】
図３に本実施例の別例２を示す。
【００７７】
対象者が寝返りなどの大きな体動を生じたとき、拍動間隔は大きく増減する。図７は対象者の寝返りによって生じた拍動間隔の変化と、体動の大きさを示していて、横軸は時刻、左縦軸は体動、右縦軸は拍動間隔である。ここで、体動の大きさは重力加速度を１Ｇとした単位で表している。寝返りが生じたとき（図７中、Ｐで図示）、体動は0.3Ｇを超えていて、同時におよそ1000ms程度だった拍動間隔が700ms程度まで低下している。このような拍動間隔の変化はCVHRで生じる拍動間隔の増減に酷似している。大きな体動時の拍動間隔を除外すれば、誤検出の少ない無呼吸状態検出システムを構成することができる。
【００７８】
図３に示す別例２は、図１に示す無呼吸状態検出システム１に、拍動間隔を測定するのと同時に対象者の姿勢を測定する姿勢測定手段19と、姿勢測定手段から得られた姿勢の時系列から体動を検出する体動検出手段25と、拍動間隔測定手段4により測定した心臓の拍動間隔の時系列から、体動検出手段25によって検出された体動が発生した時刻から所定の前後の時間の拍動間隔を除外し、除外後の拍動間隔の時系列を短時間平均拍動間隔演算手段８と長時間平均拍動間隔演算手段９とに転送する体動時拍動間隔除外手段26とを備えたものである。
【００７９】
具体的には、センサ２は別例１（図２）と同様であり、解析器３には、図１の解析器3に対して、姿勢受信手段22と体動検出手段25と体動時拍動間隔除外手段26とが付加されている。
【００８０】
センサ２の姿勢測定手段19は、対象者の体動を測定できるデバイスであり、例えばMEMSの３軸加速度センサを使用するのが好適である。このようなセンサは軽量、小型、低消費電力なので、被験者の負担は小さいものになる。
【００８１】
センサ２の姿勢測定手段19で測定された姿勢データは、姿勢保存手段20と姿勢送信手段21とを通して、解析器３の姿勢受信手段22へ転送される。
【００８２】
解析器３の姿勢受信手段22で受信した姿勢データは、体動検出手段25へ転送される。
【００８３】
体動検出手段25は、姿勢データから体動の有無を探索する。対象者が、例えば寝返りをうつとき、対象者は加速度運動する。センサ２が対象者に常に密着していれば姿勢測定手段19が、寝返りに伴う加速度を検出できる。姿勢測定手段19が３軸加速度センサで構成され、ｘ，ｙ，ｚ軸の出力それぞれを(Ａｘ，Ａｙ，Ａｚ)とし、静止状態の加速度センサが重力のみを検出するため√(Ａｘ２＋Ａｙ２＋Ａｚ２)＝１であることを使って、
【００８４】
【数６】
2020092907.tif 000008
【００８５】
と定義する。図７からもわかるように、Ｍが所定の値よりも大きければ体動があると判定できる。ここで言う所定の値は0.1から0.3が好適である。体動検出手段25は、このように体動の有無を検出し、その結果を体動時拍動間隔除外手段26へ転送する。
【００８６】
体動時拍動間隔除外手段26は、体動検出手段25から体動の有無を受け取ると共に、拍動間隔受信手段７から拍動間隔データを受け取り、体動のあった時刻から所定の前後の時間にある拍動間隔データを除外する。体動のあった前後は一定時間拍動間隔が変動するためである。ここで言う所定の時間は60秒以上が好適である。体動時拍動間隔除外手段26は体動に関する拍動間隔を除外した後、この拍動間隔データを短時間平均拍動間隔演算手順８及び長時間平均拍動間隔演算手段９へ転送する。
【００８７】
従って、別例２においては、一層誤検出の少ない無呼吸状態検出システムを実現可能となる。
【００８８】
図４は、別例１及び別例２の構成を共に備えた別例３である。
【００８９】
別例１（または別例２）のセンサ２と、図１における解析器３に姿勢受信手段22と体位判定手段23と所定体位拍動間隔除外手段24と体動検出手段25と体動時拍動間隔除外手段26とを備えたものとを用いれば、より一層誤検出の少ない無呼吸状態検出システム１を構成できる。
【００９０】
なお、解析器３の所定体位拍動間隔除外手段24と体動時拍動間隔除外手段26の順序（拍動間隔を除外する順序）は当然入れ替えても良い。
【００９１】
本実施例は上述のように構成したから、CVHRの特徴を拍動間隔から簡単な演算で検出して無呼吸状態の有無を把握できる実用的なものとなる。
【符号の説明】
【００９２】
１  無呼吸状態検出システム
３  解析器
４  拍動間隔測定手段
８  短時間平均拍動間隔演算手段
９  長時間平均拍動間隔演算手段
10  拍動間隔谷演算手段
11  拍動間隔谷判定手段
14  所定判定時間内谷数積算手段
16  谷数比較手段
19  姿勢測定手段
23  体位判定手段
24  所定体位拍動間隔除外手段
25  体動検出手段
26  体動時拍動間隔除外手段
【図１】
2020092907.tif 000009
【図２】
2020092907.tif 000010
【図３】
2020092907.tif 000011
【図４】
2020092907.tif 000012
【図５】
2020092907.tif 000013
【図６】
2020092907.tif 000014
【図７】
2020092907.tif 000015