(57)【特許請求の範囲】
【請求項１】
少なくとも２つの所定波長の電磁波を受信した強さから脈拍数を決定する脈拍数決定方法であって、
受信された前記電磁波の強さを少なくとも１つの第１および第２の電気信号（：Ｌ１［ｉ］、Ｌ２［ｉ］）に変換するステップと、
前記第１の信号の連続した第１の一時平均値（：Ｌｌａ，ｖｇ［ｉ］）を形成するステップと、
前記第２の信号の連続した第２の一時平均値（Ｌ２ａｖｇ［ｉ］）を形成するステップと、
前記第１の信号および前記第１の一時平均値から第１の交番成分（ＡＣ１［ｉ］）を連続的に決定するステップと、
前記第２の信号および前記第２の一時平均値から第２の交番成分（ＡＣ２［ｉ］）を連続的に決定するステップと、
前記２つの交番成分の連続差値（△ＡＣ［ｉ］）を決定するステップと、
前記連続差値（△ＡＣ［ｉ］）のゼロ通過に基づき血液の流れの脈拍数を決定するステップと、
を備えたことを特徴とする脈拍数決定方法。
【請求項２】
人間の組織により反射されたまたは人間の組織を通して伝えられた少なくとも２つの所定波長の電磁波を受信した強さから人間の組織を通る血液の流れの脈拍数を決定する脈拍数決定装置であって、
受信された前記電磁波の強さを少なくとも１つの第１および第２の電気信号（：Ｌ１［ｉ］、Ｌ２［ｉ］）に変換する手段と、
前記第１の信号の連続した第１一時平均値（Ｌ１ａｖｇ［ｉ］）を形成する手段と、
前記第２の信号の連続した第２一時平均値（Ｌ２ａｖｇ［ｉ］）を形成する手段と、
前記第１の信号および前記第１の一時平均値から第１の交番成分（ＡＣ１［ｉ１］）を連続的に決定する手段と、
前記第２の信号および前記第２の一時平均値から第２の交番成分（ＡＣ２［ｉ］）を連続的に決定する手段と、
前記２つの交番成分の連続差値（△ＡＣ［ｉ］）を決定する手段と、
前記連続差値（△ＡＣ［ｉ］）のゼロ通過に基づき血液の流れの脈拍数を決定する手段と、を備えたことを特徴とする脈拍数決定装置。
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】
本発明は、人間の組織により反射された、または人間の組織を通して伝えられた少なくとも二つの所定波長の電磁波の強さに基づいて脈拍数を決定するための脈拍数決定方法および脈拍数決定装置に関する。
【０００２】
【従来の技術】
本発明による方法および装置に関して、例えば酸素飽和度の測定および計算に有利に使用できる従来の技術を説明する。
【０００３】
酸素飽和度は、患者の状態を評価するのに臨床的に非常に適切なパラメータである。特に手術室では、血液の酸素飽和度は、患者状態、その酸素の供給、および他の生理学的因子を示す。
【０００４】
患者の酸素飽和度の、非常に精密な値を得る一つの可能性は、血液試料を採取し、これを血液ガス分析器で分析することである。この方法は非常に精密であるが、侵入的方法であり、これを頻繁に行うことはできない、すなわち連続監視ができないということを意味する。したがって、酸素飽和値の重要な変化を見落とすことがある。よって、この侵入的方法は患者を監視するのに好ましい方法ではないことが理解される。
【０００５】
故に、酸素飽和度を非侵入的に測定することが非常に望まれる。これはパルス酸素計測法といわれる方法により行うことができる。
【０００６】
パルス酸素計は通常、波長の異なる二つ以上の光源を備えている。光は人間の身体に照射され、身体から伝えられた光の強さ、または反射された光の強さが測定される。さらに一般的に言えば、「光」は可視スペクトル内の電磁波を意味するだけではない。例えば、最も普通の酸素計は、可視スペクトル内の一つの波長、および赤外スペクトル内の他の波長を使用している。このようなパルス酸素計については、例えばＳ．Ｋａｅｓｔｌｅ，Ｆ．Ｎｏｌｌｅｒ等がヒューレット・パッカード・ジャーナルの１９９７年２月に、「パルス酸素計測法用センサの新しいファミリー」という表題で説明している。
【０００７】
酸素飽和度測定の理論に関するさらなる詳細については、この主題に関する以前の出版物、例えば米国特許第４，１６７，３３１号または欧州公開特許公報（Ａ）第２６２７７８号に記されている。（後者の資料は理論を完全に分析している。）
【０００８】
パルス酸素計測法で受信した強さは、時間によって変わる重ね合わせ（重畳）
を示している。さらに、既知の方法は、酸素飽和度を受信した強さのピーク値に関して計算する。加えて、前記既知の方法は、受信した強さに、電磁波を送信および受信するのに用いられているＬＥＤおよびホトダイオードの、例えば試験試料としての指に対する変位により生ずるいわゆる「モーション・アーチファクト（ｍｏｔｉｏｎ  ａｒｔｉｆａｃｔｓ）」が重畳しているとき、精密な飽和計算を行うことができない。したがって、動脈酸素飽和度を決定する既知の方法は、受信した強さの時間によって変わる重畳のため不正確であり、同じまたは高い周波数の交番信号時間スペクトルに影響されやすく、既知の決定法はこの場合有効な結果を得ない。
【０００９】
【発明が解決しようとする課題】
本発明の目的は、人間の組織により反射された、または人間の組織を通して伝えられた少なくとも二つの所定波長の電磁波の強さに基づき、人間の組織を通る血液の流れの脈拍数を決定する方法および装置を提供することである。
【００１０】
【課題を解決するための手段】
本発明の目的は、請求項１による方法および請求項２による装置によって達成される。
【００１１】
請求項１に記載の発明は、少なくとも２つの所定波長の電磁波を受信した強さから脈拍数を決定する脈拍数決定方法であって、
受信された前記電磁波の強さを少なくとも１つの第１および第２の電気信号（：Ｌ１［ｉ］、Ｌ２［ｉ］）に変換するステップと、
前記第１の信号の連続した第１の一時平均値（：Ｌｌａ，ｖｇ［ｉ］）を形成するステップと、
前記第２の信号の連続した第２の一時平均値（Ｌ２ａｖｇ［ｉ］）を形成するステップと、
前記第１の信号および前記第１の一時平均値から第１の交番成分（ＡＣ１［ｉ］）を連続的に決定するステップと、
前記第２の信号および前記第２の一時平均値から第２の交番成分（ＡＣ２［ｉ］）を連続的に決定するステップと、
前記２つの交番成分の連続差値（△ＡＣ［ｉ］）を決定するステップと、
前記連続差値（△ＡＣ［ｉ］）のゼロ通過に基づき血液の流れの脈拍数を決定するステップと、
を備えたことを特徴とする。
【００１３】
本発明による方法は、決定された交番成分の連続差値のゼロ通過に基づき、潅流の脈拍数を決定することができる。潅流が低いときおよび別の妨害が発生したとき、脈拍数の決定に隙間がしばしば生ずる既知の方法とは対照的に、本発明による方法は、散乱幅の小さい、したがって標準偏差の小さい連続値を、または１パルス長にわたり平均化された、散乱幅の小さい値、を連続的に得ることができる。
【００１４】
請求項２に記載の発明は、人間の組織により反射されたまたは人間の組織を通して伝えられた少なくとも２つの所定波長の電磁波を受信した強さから人間の組織を通る血液の流れの脈拍数を決定する脈拍数決定装置であって、
受信された前記電磁波の強さを少なくとも１つの第１および第２の電気信号（：Ｌ１［ｉ］、Ｌ２［ｉ］）に変換する手段と、
前記第１の信号の連続した第１一時平均値（Ｌ１ａｖｇ［ｉ］）を形成する手段と、
前記第２の信号の連続した第２一時平均値（Ｌ２ａｖｇ［ｉ］）を形成する手段と、
前記第１の信号および前記第１の一時平均値から第１の交番成分（ＡＣ１［ｉ１］）を連続的に決定する手段と、
前記第２の信号および前記第２の一時平均値から第２の交番成分（ＡＣ２［ｉ］）を連続的に決定する手段と、
前記２つの交番成分の連続差値（△ＡＣ［ｉ］）を決定する手段と、
前記連続差値（△ＡＣ［ｉ］）のゼロ通過に基づき血液の流れの脈拍数を決定する手段と、を備えたことを特徴とするものである。
【００１５】
【発明の実施の形態】
以下に、本発明の好適な実施の形態を付属の図面を参照して詳細に説明する。
【００１６】
本発明による方法を実行するために用いられる装置は、少なくとも二つの送信手段、好適には、発光ダイオード（ＬＥＤ）と、少なくとも一つの光電受信器、例えばホトトランジスタまたはホトダイオードとを備えている。ＬＥＤは異なる波長の光を人間の組織、例えば人間の指または人間の耳に放射する。一方のＬＥＤは、例えば６５０ｍｍ（赤）の波長の光を放射することができるが、他方のＬＥＤは、９４０ｍｍ（赤外）の波長の光を放射することができる。これらＬＥＤにより放射された光は、人間の組織内に伝達、または人間の組織によって反射され、伝達または反射された光は、両波長に応答する受信器とすることもできる光電受信器により受信される。光電受信器により発生された信号は、例えば電流／電圧変換器に送られ、続いて前記信号に基づいてデジタル値を発生するようにアナログ／デジタル変換器に供給される。デジタル値は、脈拍数を検出する方法を、本発明にしたがって実行する処理ユニットに供給される。この処理ユニットを、記憶装置、入力装置および出力装置（例えば、キーボード、出力画面およびプリンタ）の他、さらに他の処理システムに伝えるインターフェースを備えた当該技術において通常使用されている形式のマイクロプロセッサとすることができる。
【００１７】
図１〜図３を参照して、本発明の第１の形態による方法の好適な実施形態を最初に以下に説明する。該方法は、原波（raw-wave）データから最大Ｓｐ０２（酸素飽和度）情報を得ることができるので、脈拍数および潅流などの確実なＳｐ０２値を低潅流状態でも保証することができる。
【００１８】
可能なかぎり「純粋な」動脈交番成分値（ＡＣ値）を得るためには、他の全ての、時間により変わる重畳を、可能な限り前記交番成分値から分離しなければならない。この目的で、本発明による方法は、原波Ｌ１［ｉ］およびＬ２［ｉ］に基づき連続平均値Ｌ１ａｖｇ［ｉ］、Ｌ２ａｖｇ［ｉ］を形成するステップを備えている。Ｌ１［ｉ］は例えば赤信号を、Ｌ２［ｉ］は例えば赤外信号を言う。また、これら連続平均値を基線ということができる。赤信号に対する基線（baseline(red)）は、例えば下記の式（１）によって決定される。
【００１９】
【数１】
0004056697.tif 000002
【００２０】
ここでＴはパルス周期あたりのサンプル数である。脈拍数（ＰＲ）が毎分７５鼓動（ｂｐｍ）であるとき、得られるパルス周期は８００ｍｓである。したがって、Ｔを１００にプリセットすると、サンプルは８ｍｓおきに検出される。赤外信号に対する基線は同様にして決定される。
【００２１】
図１において、赤信号および赤外信号に対して、原波Ｌ１［ｉ］を曲線１０として、原波Ｌ２［ｉ］を曲線１２として、また決定された基線Ｌ１ａｖｇ［ｉ］を曲線１４として、基線Ｌ２ａｖｇ［ｉ］を曲線１６として図示してある。一定成分動力学も収縮期（systole）と拡張期（diastole）との間で有効である。
【００２２】
この点で、図４を除く各図において、横座標はサンプル数を規定するが、縦座標は正規化振幅を規定するということに注意すべきである。
【００２３】
本発明による方法において、交番成分を決定するためのＡＣ分離を、今度は差形成（difference formation）を用いて基線の助けを借りて行う。赤信号について、交番成分ＡＣ１［ｉ］を下記の式（２）にしたがって求める。
【００２４】
【数２】
0004056697.tif 000003
【００２５】
赤外信号に対する交番成分ＡＣ２［ｉ］は、類似した仕方により決定される。
その結果、ＡＣ情報の他に、基線ドリフトが上の差形成により連続的に補償される二つのパルス波が存在する。
【００２６】
原波および基線から得られた交番成分に基づき、下記式（３）にしたがって連続比Ｒ［ｉ］を直接得ることができる。
【００２７】
【数３】
0004056697.tif 000004
【００２８】
しかしながら、この点について、比を計算するために、０／０の状態がある場合もその付近（neighborhood）を使用することもないということは決定的に重要である。すなわち、交番成分ＡＣ１、ＡＣ２およびその付近の全部０の通過は、計算から除外される。例えば、ゼロ通過前およびゼロ通過後のサンプルを、計算から除外することができる。交番成分ＡＣ１である曲線２０、交番成分ＡＣ２である曲線２２、および０点を除く連続比Ｒ［ｉ］である曲線２４を図２に示す。標準偏差が＋０．０３である正規化平均値１．０７が、ここでは５１２サンプルの時間にわたる比に対して得られる。
【００２９】
ゼロ通過およびそれらの付近を除くことにより、値はいくぶん予備フィルタリングされるか、または純粋の動脈情報が入っている値にある程度まで減少する。得られた比の値を、それぞれの連続Ｓｐ０２［ｉ］値を決定するのに使用できる。個別Ｓｐ０２［ｉ］値を、例えば変換器コードによる既知の方法で得ることができる。
【００３０】
一例で得られたＳｐ０２値を、図３に曲線３０として示してある。これらＳｐ０２値は耳センサを用いて、０点とその付近とを除いて、低潅流状態で得られ、連続して得られたものである。図３には、赤信号および赤外信号の交番成分も曲線３２、３４として示されている。図３の縦座標は１００％＝１０００に対して尺度を調節してある。横座標はサンプル数ｉを表している。未フィルタリング拍動間Ｓｐ０２値は、既知のアルゴリズム、すなわち「ＣＭＳ−Ｓｐ０２アルゴリズム」を用いて得られたものであるが、図３に本発明の方法により得られ、かつ図３のＳｐ０２値の２倍の大きさの尺度で、マーカ３６、３８として別に示されている。
【００３１】
図示した例では、５１２サンプルにわたり決定されたＳｐ０２値の平均値は９６．７％であり、標準偏差は±３．０％である。比較のため５，１２０サンプル（約４０秒の事態継続（episode）期間）にわたる平均値を計算した。この平均値は９７．３％であり、標準偏差は±０．６％である。
【００３２】
既知のアルゴリズムと比較して、本発明による方法は、標準偏差に関して、特に平均値を５，１２０サンプルにわたり計算したとき、明確な改善を表わしている。さらに、本発明により得られた連続比およびＳｐ０２値は、全パルス周期より長い期間安定であり、その散乱は心拡張期に到達しても低い。
【００３３】
本発明による上述の方法において比を計算するとき、最大値（心拡張期）の代わりに平均値を考慮するという事実に鑑み、潅流によって変わる小さい補正係数が古典的計算である式（４）、すなわち式（５）との比較を生ずる。
【００３４】
【数４】
0004056697.tif 000005
【００３５】
【数５】
0004056697.tif 000006
【００３６】
ここでＬ１ｍｉｎ、Ｌ１ｍａｘ、Ｌ２ｍｉｎ、Ｌ２ｍａｘは、それぞれ赤色波の最小値、赤外波の最小値、赤色波の最大値、赤外波の最大値を表す。しかしながら、補正係数は、潅流が１０％以上（＞１０％）の値に達するまでＳｐ０２で有意（significant）にならず、これらの場合にはこれを補正することができる。
【００３７】
低潅流および別の妨害の場合には、上述のＣＭＳアルゴリズムを用いるとき、それぞれの拍動およびＳｐ０２に隙間が頻繁に生ずるが、ここで前に説明した方法を用いることにより、値を連続値の形で、または１パルス長にわたって平均した値の形で、連続して得ることができ、既知のアルゴリズムと比較して、前記値の散乱幅は半分である。確認された連続交番成分を別に使用してそこから脈拍数を得ることができ、これは、それぞれのゼロ通過、すなわち符号反転を検出することにより好適に行われる。これにより、中間最大値（「ダイクロチー（Ｄｉｃｈｒｏｔｉｅ）」）で、および古典的「ピーク探索」法を使用するときアーチファクトの影響により、しばしば生ずる可能な不良トリガが減少する。
【００３８】
連続ＤＣ比（DC#Ratio）が下の式（６）にしたがって形成されるとき、興味あるパラメータも得られる。
【００３９】
【数６】
0004056697.tif 000007
【００４０】
この比は、モーション・アーチファクトが生じても、センサがその位置を変えないかぎり非常に安定したままである。短い低潅流またはアーチファクト状態では、監視Ｓｐ０２または保持機能をこのように取り入れて前記状態をつなぐことができる。上の例について、このＤＣ比を５，１２０サンプルにわたり平均して図３に曲線４０として示してある。動脈Ｓｐ０２の突然の降下または変化もこの値を変化させる。
【００４１】
下記では、本発明の第２の形態による方法の実施形態を図４〜図１０に基づき説明する。
【００４２】
上述の連続基線法を用いて、低周波妨害成分（ｆｓ＿ｌｏｗ＜パルスおよび平均割合がそれぞれ１／Ｔｐおよび１／Ｔ）を全て排除した。したがって、ＡＣスペクトル内で妨害（ｄｉｓｔｕｒｂａｎｃｅ）としてなお有効な唯一の成分は、同じまたは高い周波数を有する成分である。主として、ＬＥＤのまたはホトダイオードの指または動脈に対する変位により生ずるモーション・アーチファクトは、このような障害（Ｉｎｔｅｒｆｅｒｅｎｃｅ）である。このような運動を図解する概略図を図４に示す。
【００４３】
このようなアーチファクトが生ずると、同じ変位ａ［ｉ］が両ＡＣ信号、ＡＣ１［ｉ］およびＡＣ２［ｉ］、について生ずる。したがって、第１近似で、それぞれの妨害変調信号ＡＣＳ［ｉ］は、以下の式（７）および式（８）により与えられる。
【００４４】
【数７】
0004056697.tif 000008
【００４５】
【数８】
0004056697.tif 000009
【００４６】
図５は、二つのＡＣ信号、ＡＣ１［ｉ］である曲線５０、およびＡＣ２［ｉ］である曲線５２を５００個のサンプルｉの窓上に示した図である。その周波数が脈拍数より高いか等しいモーション・アーチファクトから生じ、かつホトダイオードの位置で強さの変化を生ずる妨害ａ［ｉ］を図６で曲線６０として図示してある。また、図７に両妨害信号ＡＣ１ｓ［ｉ］である曲線７０、およびＡＣ２ｓ［ｉ］である曲線７２を図示する。図７は、妨害が生じている期間中、純粋ピークも正しいゼロ通過も、特に妨害の振幅が動脈信号振幅に等しいか大きくなる場合に、パルス継続時間の決定に対して存在しないことを示している。これらの場合には、「ピーク探索」法による、飽和の既知の「古典的」計算を適用できない。
【００４７】
妨害信号成分はそれだけでは知られない。しかしながら、信号ＡＣ１ｓ［ｉ］およびＡＣ２ｓ［ｉ］だけにこの妨害信号成分が重畳しているため、妨害が両波長について同じ位相および振れを生ずると仮定すれば、前記妨害を以下の式（９）に示すような差の形成により排除することができる。
【００４８】
【数９】
0004056697.tif 000010
【００４９】
差曲線△ＡＣ［ｉ］を図８に曲線８０として示してある。図８は、差曲線のそれぞれのゼロ通過から得ることができる正しいパルス周期ＴｐＵＬＳＥをも示している。このようなモーション・アーチファクトが生ずるとき、脈拍数を決定できない既知の方法とは対照的に、ここで前に説明した方法は、このようなモーション・アーチファクトが存在しても、人間の組織の潅流の正確な脈拍数を決定できる。
【００５０】
正しいアーチファクト排除の必要条件は、赤および赤外に対するＤＣ成分がほぼ等しいということである。これは、例えば不等ＬＥＤ強度制御を行うことにより達成することができる。比≒１（式（１３））の場合には、差△ＡＣは０であるが、ＳＰ０２は、アーチファクトが重畳していても、例えば８５％に与えられる。この場合について正しい脈拍数を計算するには、例えば異なるＬＥＤ強度を使用することにより、有限差を人工的に計算することができる。
【００５１】
既に述べたとおり、本発明の第１  の形態を参照して前に説明した方法は、低周波妨害信号に対して良い結果を与える。下記では、上に説明した差値△ＡＣ［ｉ］を用いて比を形成することができ、例えばモーション・アーチファクトが存在する場合に成分の濃度、好適な実施形態の場合には動脈酸素飽和度をも決定できる方法を説明する。
【００５２】
比Ｒ［ｉ］を式（１０）のように式（１）と同等の形で書くことができる。
【００５３】
【数１０】
0004056697.tif 000011
【００５４】
比ＡＣ／ＤＣ＜５％（潅流＜１０％）の場合には、目立った誤差を生ぜずに下記の式（１１）のように単純化したものを用いることができる。
【００５５】
【数１１】
0004056697.tif 000012
式（１１）に以下の式（１２）を代入すると、式（１３）のようになる。
【００５６】
【数１２】
0004056697.tif 000013
【００５７】
【数１３】
0004056697.tif 000014
【００５８】
また、ＡＣ１［ｉ］を置き換える代わりに、上記の式（８）および式（９）でＡＣ２［ｉ］を同等に置き換えることもできる。したがって、下に説明する方法を交番成分ＡＣ２［ｉ］の代わりに、交番成分ＡＣ１［ｉ］を用いて行うことができる。
【００５９】
モーション・アーチファクトにより生ずる妨害の場合には、無妨害信号ＡＣ２［ｉ］の代わりにＡＣ２ｓ［ｉ］が有効である。しかし、△ＡＣ［ｉ］は、妨害ばかりでなく無妨害に対しても同じままであり、また係数Ｌ２ａｖｇ［ｉ］／Ｌ１ａｖｇ［ｉ］はアーチファクトの影響があっても不変である。
【００６０】
したがって、アーチファクトが存在するとき、主要な仕事はこのアーチファクトを認識することであり、一方で最も正しい可能な（アーチファクトが全く無い）交番成分ＡＣ２［ｉ］（または代わりにＡＣ１［ｉ］）を見付けることである。
【００６１】
妨害の場合には、ゼロ通過またはピークのようなスポット曲線決定を行う代わりに、計算領域またはサブ周期にわたり積分した値の使用が有利である。この目的で、それぞれの値を、例えば連続ＡＣ差の半周期長にわたり積分することができる。
【００６２】
連続ＡＣ差について、下記の式（１４）が得られる。
【００６３】
【数１４】
0004056697.tif 000015
【００６４】
他の数量ＡＣ２［ｉ］、Ｌ１ａｖｇ［ｉ］、Ｌ２ａｖｇ［ｉ］の積分値を、対応する仕方で確認できる。これらの数量について、値ＡＣ２［Ｔ／２］、Ｌ１ａｖｇ［Ｔ／２］、Ｌ２ａｖｇ［Ｔ／２］が得られる。したがって、以下の式（１５）が各半パルス周期での比Ｒについて得られる。
【００６５】
【数１５】
0004056697.tif 000016
【００６６】
比Ｒ［Ｔ／２］の基礎を成す数量の中で、妨害により変化するＡＣ２［Ｔ／２］だけが不可欠である。妨害の存在が検出されると、本発明は、先行するサブ周期中に決定されているアーチファクト・フリー（ａｒｔｉｆａｃｔｓ  ｆｒｅｅ）のＡＣ２′［Ｔ／２］の値を使用して、式（１５）に従い比Ｒ′［Ｔ／２］を決定する。この方法を、妨害の認識後、少なくとも特定の期間、および差値△ＡＣ［Ｔ／２］が実質的に変わらないかぎり行うことができる。妨害が検出されなければ、瞬時的に検出されたＡＣ２値を使用し続けて積分比Ｒ［Ｔ／２］を決定する。
【００６７】
ＡＣ２に対する無妨害値の存在を、０．６・ＡＣ２＜ＡＣ２ｓ＜１．４・ＡＣ２であるかぎり、仮定することができる。値ＡＣ２ｓがこの範囲外にあれば、これは妨害の存在を示すと判定することができる。妨害はΔＡＣ［ｉ］とＡＣ２ｓ［ｉ］とのゼロ通過の間の相違に基づいて判定することもできる。なぜならば、妨害がなければこれらゼロ通過は一致するからである。値ΔＡＣ（Ｔ／２）の他、ＤＣ比Ｌ２ａｖｇ／Ｌ１ａｖｇも飽和の変化に応答するが、前に述べた人工生成物には応答しない。これにより上述の計算は、積分比値Ｒ’［Ｔ／２］に基づき成分の濃度を決定することにより、成分の濃度、例えば人間の組織の中の動脈酸素飽和度を検出するのに優れた近似を与える。
【００６８】
図９に、上述の近似法にしたがって計算した比を、無妨害信号である曲線２４に基づく第１の形態に従う計算と比較して曲線９０で示してある。
【００６９】
前に記した方法の効率は、図１０からも見ることができる。図１０は、耳センサを用いて検出された弱い信号、および軽い揺れにより生ずる強いアーチファクトの影響のある状態の詳細を示す。ステップ周波数は脈拍数とわずかだけ異なる。
【００７０】
通常のピーク探索アルゴリズムはここでは、「メルリンＳｐ０２アルゴリズム」第４改訂版として知られているアルゴリズムから得られたステップ周波数（ＣＭＳ−Ｓｐ０２値を参照）、ドット１００に応答してトリガし、５１２サンプルにわたる平均飽和値８４．６％を与え、標準偏差は±１４．８％である。図４〜図８を参照して前に説明した方法によれば、５１２サンプルにわたる平均飽和値９４．９％、標準偏差±３．９％が得られ（曲線１０２を参照）、これは正常値に対応する。
【００７１】
図１０のドット１０４は、上述の通常の方法を用いて得られた１２７．４ｂｐｍの平均脈拍数を表している。本発明による方法によれば、さらに精密な脈拍数１１３ｂｐｍ（ドット１０６を参照）が得られる。
【００７２】
図１０は別に、二つの妨害交番信号ＡＣ１ｓ［ｉ］およびＡＣ２ｓ［ｉ］ならびに値ΔＡＣ・５（曲線１０８を参照）を示すが、曲線１１０は、図１０の連続ＤＣ比を表している。本発明の第１の形態に従う方法の場合に得られる脈拍数は、図１０にドット１１２により表されている。
【００７３】
下記に、アーチファクトの場合の比決定の他の方法を説明する。第１のステップは、図１−図３に関して、前に説明した仕方で連続比Ｒ［ｉ］を決定することである。さらに、ＡＣ差分である△ＡＣ［ｉ］を形成し、パルス周期をそのゼロ通過から得ることができる。酸素飽和度を決定するための比として、ＡＣゼロ点を通過せず、振幅の振れが△ＡＣ［ｉ］値と同じ方向（ＳＩＧＮ）であるＡＣ２ｓ［ｉ］値によってだけ形成される１サンプルだけを使用する。これにより無妨害位相を有するさらに多数の信号成分が考慮される可能性が増大するが、反対方向に向かうものは全て考慮外とされる。これは、弱い妨害振幅から中間の大きさの妨害振幅までに対して特に役立つことが証明されている、比較的簡単な方法である。しかしながら、図１０に示してある、５１２サンプルにわたって平均された結果の、飽和の信号曲線である曲線１１４は強い妨害振幅を有する場合を基準としており、この場合、ＡＣ［ｉ］はほぼａ［ｉ］に等しい。この場合には、平均飽和値９０．８％、標準偏差士９．６％が得られた。
【００７４】
以下に本発明の実施の形態を要約する。
１．  人間の組織により反射されたまたは人間の組織を通して伝えられた少なくとも二つの所定波長の電磁波の強さから成分の濃度を決定する成分濃度決定方法であって、
１．１  受信した電磁波の強さを少なくとも一つの第１のおよび一つの第２の電気信号（Ｌ１［ｉ］、Ｌ２［ｉ］）に変換するステップと、
１．２  第１の信号の連続第１一時平均値（Ｌ１ａｖｇ［ｉ］）を形成するステップと、
１．３  第２の信号の連続第２一時平均値（Ｌ２ａｖｇ［ｉ］）を形成するステップと、
１．４  前記第１の信号および前記第１の平均値から第１の交番成分（ＡＣ１［ｉ］）を連続的に決定するステップと、
１．５  前記第２の信号および前記第２の平均値から第２の交番成分（ＡＣ２［ｉ］）を連続的に決定するステップと、
１．６  交番成分のゼロ通過に近い領域を除き、前記第１および第２の信号からの他、前記第１および第２の連続一時的平均値からも連続比（Ｒ［ｉ］）を決定するステップと、
１．７  前記連続比（Ｒ［ｉ］）から成分の濃度を得るステップと、を備えている成分濃度決定方法。
【００７５】
２．  人間の組織を通る血液の流れの酸素飽和度を決定する上記１に記載の成分濃度決定方法。
【００７６】
３．  別に、前記第１の交番成分（ＡＣ１［ｉ］）または前記第２の交番成分（ＡＣ２［ｉ］）のゼロ通過に基づき血液の前記流れの脈拍数を決定するステップを備えている上記２に記載の成分濃度決定方法。
【００７７】
４．  前記第１および第２の連続一時平均値（Ｌ１ａｖｇ［ｉ］、Ｌ２ａｖｇ［ｉ］）はそれぞれ、前記第１の信号または前記第２の信号のパルス周期の継続期間にわたり決定される上記１〜３の一つに記載の成分濃度決定方法。
【００７８】
５．  成分の濃度を決定する基となる連続比（Ｒ［ｉ］）は、前記第１の信号（Ｌ１［ｉ］）、前記第２の信号（Ｌ２［ｉ］）、前記第１の連続一時平均値（Ｌ１ａｖｇ［ｉ］）および前記第２の連続一時平均値（Ｌ２ａｖｇ［ｉ］）を用いて、等式Ｒ［ｉ］＝ｌｎ（Ｌ１［ｉ］／Ｌ１ａｖｇ［ｉ］）／ｌｎ（Ｌ２［ｉ］／Ｌ２ａｖｇ［ｉ］）により決定される上記１〜４の一つに記載の成分濃度決定方法。
【００７９】
６．  前記第１の交番成分（ＡＣ１［ｉ］）は前記第１の信号Ｌ１［ｉ］から前記連続第１一時平均値Ｌ１ａｖｇ［ｉ］を減ずることにより形成され、前記第２の交番成分（ＡＣ２［ｉ］）は前記第２の信号Ｌ２［ｉ］から前記連続第２一時平均値Ｌ２ａｖｇ［ｉ］を減ずることにより形成される上記１〜５の一つに記載の成分濃度決定方法。
【００８０】
７．  人間の組織により反射されたまたは人間の組織を通して伝えられた少なくとも二つの所定波長の電磁波の強さから成分の濃度を決定する成分濃度決定方法であって、
７．１  前記受信された電磁波の強さを少なくとも一つの第１のおよび一つの第２の電気信号（Ｌ１［ｉ］、Ｌ２［ｉ］）に変換するステップと、
７．２  前記第１の信号の連続第１一時平均値（Ｌ１ａｖｇ［ｉ］）を形成するステップと、
７．３  前記第２の信号の連続第２一時平均値（Ｌ２ａｖｇ［ｉ］）を形成するステップと、
７．４  前記第１の信号および前記第１の平均値から第１の交番成分（ＡＣ１［ｉ］）を連続的に決定するステップと、
７．５  前記第２の信号および前記第２の平均値から第２の交番成分（ＡＣ２［ｉ］）を連続的に決定するステップと、
７．６  前記二つの交番成分の連続差値（ΔＡＣ［ｉ］）を決定するステップと、
７．７  前記第１の平均値を前記連続差値のサブ周期にわたり積分することにより第１の積分平均値（Ｌ１ａｖｇ［Ｔ／２］）を決定するステップと、
７．８  前記第２の平均値を前記サブ周期にわたり積分することにより第２の積分平均値（Ｌ２ａｖｇ［Ｔ／２］）を決定するステップと、
７．９  前記第１の交番成分または前記第２の交番成分を前記サブ周期にわたり積分することにより積分交番成分（ＡＣ１［Ｔ／２］、ＡＣ２［Ｔ／２］）を決定するステップと、
７．１０  前記差値を前記サブ周期にわたり積分することにより第１の積分差値（ΔＡＣ［Ｔ／２］）を決定するステップと、
７．１１  瞬時的第１または第２の交番成分を、妨害皆無または妨害を受けていると判定するステップと、
７．１２  前記瞬時的第１または第２の交番成分が妨害皆無と判定されれば、
・前記積分差値（ΔＡＣ［Ｔ／２］）、前記積分交番成分（ＡＣ２［Ｔ／２］）、および前記二つの積分平均値（Ｌ１ａｖｇ［Ｔ／２］、Ｌ２ａｖｇ［Ｔ／２］）に基づき積分比値（Ｒ［Ｔ／２］）を決定し、
・前記積分比値（Ｒ［Ｔ／２］）から成分の濃度を得るステップと、
７．１３  その他の場合には、
・前記積分差値（ΔＡＣ［Ｔ／２］）、決定された前記積分交番成分が妨害皆無と判定された前記先行サブ周期中に決定された積分交番成分（ＡＣ２’［Ｔ／２］）、および比が人工生成物のもとで一定であるが、突然の濃度変化時には一定でない前記二つの積分平均値（Ｌ１ａｖｇ［Ｔ／２］、Ｌ２ａｖｇ［Ｔ／２］）に基づき、積分比値（Ｒ’［Ｔ／２］）を決定し、
・前記積分比値（Ｒ’［Ｔ／２］）から成分の濃度を得るステップと、を備えている成分濃度決定方法。
【００８１】
８．  ステップ７．７〜ステップ７．１０のそれぞれに使用した前記サブ周期は、前記連続差値（ΔＡＣ［ｉ］）の半波である上記７に記載の成分濃度決定方法。
【００８２】
９．  ＡＣを妨害皆無交番成分、ＡＣｓを妨害を受ける交番成分としたとき、ステップ７．１１において、前記第１の交番成分または前記第２の交番成分は、０．６・ＡＣ＜ＡＣｓ＜１．４・ＡＣのとき妨害皆無と判定される上記７または８に記載の成分濃度決定方法。
【００８３】
１０．  ステップ１１で、前記第１の交番成分または前記第２の交番成分は、前記差値（ΔＡＣ［ｉ］）のゼロ通過、および前記第１の交番成分または前記第２の交番成分（ＡＣ１［ｉ］、ＡＣ２［ｉ］）に基づき判定される上記７または８に記載の成分濃度決定方法。
【００８４】
１１．  決定される前記成分の前記濃度は、ガス飽和、特に動脈酸素飽和度である上記１〜１０の一つに記載の成分濃度決定方法。
【００８５】
１２．  人間の組織により反射されたまたは人間の組織を通して伝えられた少なくとも二つの所定波長の電磁波の強さから人間の組織を通る血液の流れの脈拍数を決定する脈拍数決定方法であって、
１２．１  前記受信された電磁波の強さを少なくとも一つの第１のおよび一つの第２の電気信号（Ｌ１［ｉ］、Ｌ２［ｉ］）に変換するステップと、
１２．２  前記第１の信号の連続第１一時平均値（Ｌ１ａｖｇ［ｉ］）を形成するステップと、
１２．３  前記第２の信号の連続第２一時平均値（Ｌ２ａｖｇ［ｉ］）を形成するステップと、
１２．４  前記第１の信号および前記第１の平均値から第１の交番成分（ＡＣ１［ｉ］）を連続的に決定するステップと、
１２．５  前記第２の信号および前記第２の平均値から第２の交番成分（ＡＣ２［ｉ］）を連続的に決定するステップと、
１２．６  前記二つの交番成分の連続差値（ΔＡＣ［ｉ］）を決定するステップと、
１２．７  前記差値（ΔＡＣ［ｉ］）のゼロ通過に基づき血液の流れの脈拍数を決定するステップと、を備えている脈拍数決定方法。
【００８６】
１３．  人間の組織により反射されたまたは人間の組織を通して伝えられた少なくとも二つの所定波長の電磁波の強さから成分の濃度を決定する成分濃度決定装置であって、
前記受信された電磁波の強さを少なくとも一つの第１のおよび一つの第２の電気信号（Ｌ１［ｉ］、Ｌ２［ｉ］）に変換する手段と、
前記第１の信号の連続第１一時平均値（Ｌ１ａｖｇ［ｉ］）を形成する手段と、
前記第２の信号の連続第２一時平均値（Ｌ２ａｖｇ［ｉ］）を形成する手段と、
前記第１の信号および前記第１の平均値から第１の交番成分（ＡＣ１［ｉ］）を連続的に決定する手段と、
前記第２の信号および前記第２の平均値から第２の交番成分（ＡＣ２［ｉ］）を連続的に決定する手段と、
前記交番成分のゼロ通過に近い領域を除き、前記第１および第２の信号からの他、前記第１および第２の連続一時平均値からも連続比（Ｒ［ｉ］）を決定する手段と、
前記連続比（Ｒ［ｉ］）から成分の濃度を得る手段と、を備えている成分濃度決定装置。
【００８７】
１４．  人間の組織により反射されたまたは人間の組織を通して伝えられた少なくとも二つの所定波長の電磁波の強さから成分の濃度を決定する成分濃度決定装置であって、
前記受信された電磁波の強さを少なくとも一つの第１のおよび一つの第２の電気信号（Ｌ１［ｉ］、Ｌ２［ｉ］）に変換する手段と、
前記第１の信号の連続第１一時平均値（Ｌ１ａｖｇ［ｉ］）を形成する手段と、
前記第２の信号の連続第２一時平均値（Ｌ２ａｖｇ［ｉ］）を形成する手段と、
前記第１の信号および前記第１の平均値から第１の交番成分（ＡＣ１［ｉ］）を連続的に決定する手段と、
前記第２の信号および前記第２の平均値から第２の交番成分（ＡＣ２［ｉ］）を連続的に決定する手段と、
前記二つの交番成分の連続差値（ΔＡＣ［ｉ］）を決定する手段と、
前記第１の平均値を前記連続差値のサブ周期にわたり積分することにより第１の積分平均値（Ｌ１ａｖｇ［Ｔ／２］）を決定する手段と、
前記第２の平均値を前記サブ周期にわたり積分することにより第２の積分平均値（Ｌ２ａｖｇ［Ｔ／２］）を決定する手段と、
前記第１の交番成分または前記第２の交番成分を前記サブ周期にわたり積分することにより積分交番成分（ＡＣ１［Ｔ／２］、ＡＣ２［Ｔ／２］）を決定する手段と、
前記差値を前記サブ周期にわたり積分することにより第１の積分差値（ΔＡＣ［Ｔ／２］）を決定する手段と、
瞬時的第１または第２の交番成分を、妨害皆無または妨害を受けていると判定する手段と、
前記積分差値（ΔＡＣ［Ｔ／２］）、前記積分交番値（ＡＣ２［Ｔ／２］）、および前記二つの積分平均値（Ｌ１ａｖｇ［Ｔ／２］、Ｌ２ａｖｇ［Ｔ／２］）に基づき積分比値（Ｒ［Ｔ／２］）を決定し、さらに、前記瞬時的第１または第２の交番成分が、判定する前記手段により妨害皆無と判定されれば、前記積分比値（Ｒ［Ｔ／２］）から成分の濃度を得る手段を備えている手段と、
前記積分差値（ΔＡＣ［Ｔ／２］）に基づき、決定された前記積分交番成分が妨害皆無と判定された前記先行サブ周期中に決定された前記積分交番成分（ＡＣ２’［Ｔ／２］）に基づき、および二つの積分平均値（Ｌ１ａｖｇ［Ｔ／２］、Ｌ２ａｖｇ［Ｔ／２］）に基づき、積分比値（Ｒ’［Ｔ／２］）を決定し、さらに、前記瞬時的第１または第２の交番成分が、判定する前記手段により妨害皆無と判定されれば、前記積分比値（Ｒ’［Ｔ／２］）から成分の濃度を得る手段を備えている手段と、を備えている成分濃度決定装置。
【００８８】
１５．  決定される前記成分の前記濃度は、ガス飽和、特に動脈酸素飽和度である上記１３または１４に記載の成分濃度決定方法。
【００８９】
１６．  人間の組織により反射されたまたは人間の組織を通して伝えられた少なくとも二つの所定波長の電磁波の強さから人間の組織を通る血液の流れの脈拍数を決定する脈拍数決定方法であって、
前記受信された電磁波の強さを少なくとも一つの第１のおよび一つの第２の電気信号（Ｌ１［ｉ］、Ｌ２［ｉ］）に変換する手段と、
前記第１の信号の連続第１一時平均値（Ｌ１ａｖｇ［ｉ］）を形成する手段と、
前記第２の信号の連続第２一時平均値（Ｌ２ａｖｇ［ｉ］）を形成する手段と、
前記第１の信号および前記第１の平均値から第１の交番成分（ＡＣ１［ｉ］）を連続的に決定する手段と、
前記第２の信号および前記第２の平均値から第２の交番成分（ＡＣ２［ｉ］）を連続的に決定する手段と、
前記二つの交番成分の連続差値（ΔＡＣ［ｉ］）を決定する手段と、
前記差値（ΔＡＣ［ｉ］）のゼロ通過に基づき血液の流れの脈拍数を決定する手段と、を備えている脈拍数決定装置。
【００９０】
【発明の効果】
上述のように本発明によれば、低潅流および別の妨害が発生したとき、脈拍数およびＳｐ０２の決定に隙間がしばしば生ずる既知の方法とは対照的に、標準偏差の小さい連続値を得ることができる。また、人間の組織により反射された、または人間の組織を通して伝えられた少なくとも二つの所定波長の電磁波の強さから、人間の組織を通る血液の流れの脈拍数を決定する方法および装置が提供される。
【図面の簡単な説明】
【図１】  記録された信号および関連する連続平均値を示す図である。
【図２】  ゼロ点を除く、決定された連続比を示す図である。
【図３】  本発明による方法の結果と既知の方法との比較を示す図である。
【図４】  送信ＬＥＤおよび受信ホトダイオードを用いる試験を示す概略図である。
【図５】  二つの代表的無妨害交番成分を示す図である。
【図６】  高い周波数を有し、モーション・アーチファクトにより生じた強さ変化を示す概略図である
【図７】  妨害交番成分を示す図である。
【図８】  図７に示した妨害交番成分から得られた差値を示す図である。
【図９】  本発明の第２の形態による方法を用いて得られた結果を示す図である。
【図１０】  本発明の第２の形態による方法を用いて得られた結果を示す図である。
【符号の説明】
１０  Ｌ１［ｉ］
１２  Ｌ２［ｉ］
１４  Ｌ１ａｖｇ［ｉ］
１６  Ｌ２ａｖｇ［ｉ］
２０  ＡＣ１［ｉ］
２４  ＡＣ２［ｉ］
３０  ＳＰ０２値
７０  ＡＣ１ｓ［ｉ］
７２  ＡＣ２ｓ［ｉ］
８０  △ＡＣ［ｉ］
【図１】
0004056697.tif 000017
【図２】
0004056697.tif 000018
【図３】
0004056697.tif 000019
【図４】
0004056697.tif 000020
【図５】
0004056697.tif 000021
【図６】
0004056697.tif 000022
【図７】
0004056697.tif 000023
【図８】
0004056697.tif 000024
【図９】
0004056697.tif 000025
【図１０】
0004056697.tif 000026