JP7019962B2 - 流体解析装置、血流解析装置および流体解析方法 - Google Patents

Description

本発明は、血液等の流体に関する情報を生成する技術に関する。

生体の血流量を測定する技術が従来から提案されている。例えば特許文献１には、生体の血管を通過した光を受光素子により受光し、その受光強度を表す検出信号のパワースペクトルと周波数の各数値との積を２００Ｈｚ以上かつ１５ｋＨｚ以下の範囲内で積算することで、生体の血流量を測定する構成が開示されている。

特開２０１２－２１０３２１号公報

ところで、周波数軸上の広い範囲にわたり均等に分布するショットノイズが、検出信号には不可避的に発生し得る。特許文献１の技術のもとでは、検出信号のパワースペクトルと周波数の各数値との積が積算されるから、周波数が高い範囲ほどショットノイズが強調される。したがって、積算範囲を厳格に選定しないと、血流量を高精度に測定できないという問題がある。なお、以上の説明では血流量の測定に着目したが、血液に代表される各種の流体を解析する多様な場面で同様の問題が想定される。以上の事情を考慮して、本発明の好適な態様は、検出信号におけるショットノイズの影響を低減して、血液等の流体を高精度に解析することを目的とする。

以上の課題を解決するために、本発明の好適な態様に係る液体解析装置は、流体を通過した光の強度を表す検出信号に対し、所定の処理帯域内で周波数が高い成分ほど抑制されるようにフィルター処理を実行する信号処理部と、前記フィルター処理後の信号の強度スペクトルから、前記流体の速度に関する情報を生成する演算処理部とを具備する。以上の態様では、所定の処理帯域内で周波数が高い成分ほど抑制されるように検出信号に対するフィルター処理が実行されるから、例えば高周波側において特に支配的となるショットノイズの影響が低減される。したがって、流体の速度に関する情報を高精度に生成することが可能である。

本発明の好適な態様において、前記演算処理部は、前記強度スペクトルにおける各周波数の強度と当該周波数との積を所定の演算範囲について積算することで、前記流体の速度に関する情報を生成し、前記処理帯域と前記演算範囲とは部分的に相互に重複する。強度スペクトルにおける各周波数の強度と当該周波数との積（加重強度）を積算する構成では、強度スペクトルにおける高周波側のショットノイズが強調される。処理帯域と演算範囲とが部分的に相互に重複する本発明の好適な態様によれば、検出信号のうちショットノイズが支配的となる高周波側の帯域の一部を含むように演算範囲を充分に広く確保した場合でも、高周波側のショットノイズの影響が低減される。したがって、流体の速度に関する情報を高精度に生成することが可能である。

本発明の好適な態様において、前記演算範囲は、第１周波数と前記第１周波数を上回る第２周波数との間の範囲であり、前記処理帯域の下端の周波数は、前記第２周波数を下回る。以上の態様では、処理帯域の下端の周波数が、演算範囲の上端である第２周波数を下回るから、演算範囲が充分に確保されるように第２周波数を高目に設定した場合でも、高周波側のショットノイズの影響が低減される。したがって、流体の速度に関する情報を高精度に生成することが可能である。

本発明の好適な態様において、前記処理帯域は、前記信号処理部による前記検出信号の抑制が６ｄＢ／Ｏｃｔ以上となる範囲である。以上の態様では、処理帯域において６ｄＢ／Ｏｃｔの周波数特性で検出信号が抑制される。したがって、強度スペクトルの周波数毎の強度に対して周波数を乗算することに起因して強調されるショットノイズの影響を効果的に低減することが可能である。

本発明の好適な態様に係る血流解析装置は、血管を通過した光の強度を表す検出信号に対し、所定の処理帯域内で周波数が高い成分ほど抑制されるようにフィルター処理を実行する信号処理部と、前記フィルター処理後の信号の強度スペクトルから前記血管内の血液の速度に関する情報を生成する演算処理部とを具備する。以上の態様では、所定の処理帯域内で周波数が高い成分ほど抑制されるように検出信号に対するフィルター処理が実行されるから、例えば高周波側において特に支配的となるショットノイズの影響が低減される。したがって、血液の速度に関する情報を高精度に生成することが可能である。

本発明の好適な態様に係る流体解析方法は、血管を通過した光の強度を表す検出信号に対し、所定の処理帯域内で周波数が高い成分ほど抑制されるようにフィルター処理を実行し、前記フィルター処理後の信号の強度スペクトルから前記血管内の血液の速度に関する情報を生成する。以上の態様では、所定の処理帯域内で周波数が高い成分ほど抑制されるように検出信号に対するフィルター処理が実行されるから、例えば高周波側において特に支配的となるショットノイズの影響が低減される。したがって、血液の速度に関する情報を高精度に生成することが可能である。

本発明の第１実施形態に係る血流解析装置の側面図である。 血流解析装置の機能に着目した構成図である。 受光部および出力回路に着目した構成図である。 演算処理部の動作を例示するフローチャートである。 信号処理部によるフィルター処理の周波数特性である。 検出信号の強度スペクトルである。 対比例におけるショットノイズの問題の説明図である。 第１実施形態による効果の説明図である。 処理帯域の下端周波数を特定する方法を説明するための図である。 処理帯域の下端周波数を求める例を示す図である。 第２実施形態における受光部および出力回路に着目した構成図である。 変形例における受光部および出力回路の着目した構成図である。

＜第１実施形態＞
図１は、本発明の第１実施形態に係る血流解析装置１００の側面図である。第１実施形態の血流解析装置１００は、被験者（生体の例示）の血流に関する生体情報（以下「血流情報」という）を非侵襲的に生成する生体測定機器であり、被験者の身体のうち測定対象となる部位（以下「測定部位」という）Ｍに装着される。第１実施形態の血流解析装置１００は、図１に例示される通り、筐体部１２とベルト１４とを具備する腕時計型の携帯機器である。すなわち、測定部位Ｍの例示である手首にベルト１４を巻回することで、血流解析装置１００は被験者の手首に装着される。第１実施形態の血流情報は、被験者の血流速度（例えば単位時間内に赤血球が動脈内を移動する距離）を血流情報として生成する。

図２は、血流解析装置１００の機能に着目した構成図である。図２に例示される通り、第１実施形態の血流解析装置１００は、制御装置２０と記憶装置２２と表示装置２４と検出装置３０とを具備する。制御装置２０および記憶装置２２は、筐体部１２の内部に設置される。図１に例示される通り、表示装置２４（例えば液晶表示パネル）は、例えば筐体部１２における測定部位Ｍとは反対側の表面に設置され、測定結果を含む各種の画像を制御装置２０による制御のもとで表示する。

図２の検出装置３０は、測定部位Ｍの状態に応じた検出信号Ｓdを生成する光学センサーモジュールである。図２に例示される通り、第１実施形態の検出装置３０は、発光部３１と受光部３２と駆動回路３３と出力回路３４とを具備する。発光部３１および受光部３２は、例えば筐体部１２において測定部位Ｍに対向する位置（典型的には測定部位Ｍに接触する表面）に設置される。なお、駆動回路３３および出力回路３４の一方または双方を検出装置３０とは別体の外部回路として設置することも可能である。

発光部３１は、測定部位Ｍに光を照射する光源である。第１実施形態の発光部３１は、狭帯域でコヒーレントなレーザー光を測定部位Ｍに照射する。例えば共振器内の共振によりレーザー光を出射するＶＣＳＥＬ（Vertical Cavity Surface Emitting LASER）等の発光素子が発光部３１として好適に利用される。第１実施形態の発光部３１は、例えば近赤外域内の所定の波長λ（λ＝８００ｎｍ～１３００ｎｍ）の光を測定部位Ｍに照射する。図２の駆動回路３３は、制御装置２０による制御のもとで発光部３１を発光させる。なお、相異なる波長の光を出射する複数の発光素子を発光部３１として利用することも可能である。また、発光部３１が出射する光の波長λは近赤外域に限定されない。

発光部３１から出射して測定部位Ｍに入射した光は、測定部位Ｍの内部で反射および散乱を繰返したうえで筐体部１２側に出射して受光部３２に到達する。具体的には、測定部位Ｍの内部に存在する動脈（例えば橈骨動脈または尺骨動脈）等の血管と血管内の血液とを通過した光が受光部３２に到達する。受光部３２は、測定部位Ｍから到来する光を受光する。第１実施形態の受光部３２は、測定部位Ｍから到達する光の強度を表す検出信号Ｓaを生成する。例えば、図３に例示される通り、受光強度に応じた電荷を発生するフォトダイオード（ＰＤ：Photo Diode）等の受光素子３２１が受光部３２として好適に利用される。検出信号Ｓaは、測定部位Ｍからの受光強度に応じたアナログの電流信号である。以上の説明から理解される通り、第１実施形態の検出装置３０は、発光部３１と受光部３２とが測定部位Ｍに対して片側に位置する反射型の光学センサーである。

受光部３２に到達する光は、測定部位Ｍの内部において静止する組織（静止組織）で反射した成分と、測定部位Ｍの内部の動脈内で移動する物体（典型的には赤血球）で反射した成分とを含有する。静止組織での反射の前後において光の周波数は変化しないが、赤血球での反射の前後においては、赤血球の移動速度（すなわち血流速度）に比例した変化量（以下「周波数シフト量」という）Δｆだけ光の周波数が変化する。すなわち、測定部位Ｍを通過して受光部３２に到達する光は、発光部３１が出射する光の周波数に対して周波数シフト量Δｆだけ変動（周波数シフト）した成分を含有する。以上の説明から理解される通り、第１実施形態の検出信号Ｓaは、測定部位Ｍの内部の血流による周波数シフトが反映された光ビート信号である。

図２の出力回路３４は、受光部３２が生成した検出信号Ｓaから検出信号Ｓdを生成する。検出信号Ｓdは、受光部３２による受光の強度に応じたデジタルの電圧信号である。前述の通り、測定部位Ｍに照射された光は、測定部位Ｍの内部の血管および血液を通過してから受光部３２に到達する。したがって、検出信号Ｓdは、被験者の血液を通過した光の強度を表す信号とも換言され得る。

図３は、第１実施形態における受光部３２および出力回路３４の構成図である。図３に例示される通り、第１実施形態の出力回路３４は、信号増幅部５１と信号処理部５２とＡ/Ｄ変換部５３とを含んで構成される。信号増幅部５１は、受光部３２から供給される検出信号Ｓaを電圧信号に変換するとともに増幅することで検出信号Ｓbを生成する。例えば、信号増幅部５１は、検出信号Ｓaを電圧信号に変換する電流/電圧変換回路と、当該電圧信号を増幅する電圧増幅回路とを含んで構成される。

信号処理部５２は、信号増幅部５１から供給される検出信号Ｓb（すなわち、流体を通過した光の強度を表す信号）に対して所定のフィルター処理を実行することで検出信号Ｓcを生成する。信号処理部５２が実行するフィルター処理の具体例については後述する。Ａ/Ｄ変換部５３は、信号処理部５２が生成したアナログの検出信号Ｓcを所定のサンプリング周波数Ｆsでデジタルの検出信号Ｓdに変換する。以上の説明から理解される通り、検出信号Ｓ（Ｓa，Ｓb，Ｓc，Ｓd）は、測定部位Ｍの内部の血流による周波数シフトが反映された光ビート信号である。

図２の制御装置２０は、ＣＰＵ（Central Processing Unit）またはＦＰＧＡ（Field-Programmable Gate Array）等の演算処理装置であり、血流解析装置１００の全体を制御する。記憶装置２２は、例えば不揮発性の半導体メモリーで構成され、制御装置２０が実行するプログラムと制御装置２０が使用する各種のデータとを記憶する。なお、制御装置２０の機能を複数の集積回路に分散した構成、または、制御装置２０の一部または全部の機能を専用の電子回路で実現した構成も採用され得る。また、図２では制御装置２０と記憶装置２２とを別個の要素として図示したが、記憶装置２２を内包する制御装置２０を例えばＡＳＩＣ（Application Specific Integrated Circuit）等により実現することも可能である。

第１実施形態の制御装置２０は、記憶装置２２に記憶されたプログラム（アプリケーションプログラム）を実行することで演算処理部６１として機能する。演算処理部６１は、検出装置３０の出力回路３４が生成した検出信号Ｓd（すなわち信号処理部５２による処理後の信号）から被験者の血流情報を生成する。第１実施形態の演算処理部６１は、前述の通り、測定部位Ｍの内部の動脈における血流速度を血流情報として算定する。

図４は、演算処理部６１が血流速度を算定するための処理のフローチャートである。例えば所定の時間毎に図４の処理が実行される。図４の処理を開始すると、演算処理部６１は、検出信号Ｓdから強度スペクトルＸを算定する（Ｓ1）。強度スペクトルＸは、検出信号Ｓdにおいて周波数軸上の各周波数ｆに対応する信号成分の強度（パワーまたは振幅）Ｐ(f)の分布である。強度スペクトルＸの算定には、離散フーリエ変換等の公知の周波数解析が任意に採用され得る。

演算処理部６１は、検出信号Ｓdの強度スペクトルＸから血流速度Ｖ（血流情報）を算定する（Ｓ2）。具体的には、第１実施形態の演算処理部６１は、以下の数式(1)の演算により血流量Ｑ（すなわち単位時間内に動脈内を移動する血液の体積）を算定し、測定部位Ｍの血管について別途に推定された断面積Ａにより血流量Ｑを除算することで血流速度Ｖ（Ｖ＝Ｑ／Ａ）を算定する。なお、相異なる時点について算定された複数の血流速度Ｖの平均値を血流情報として演算処理部６１が生成することも可能である。

数式(1)は、検出信号Ｓdの各周波数ｆと当該周波数ｆにおける強度Ｐ(f)とから血流量Ｑを算定するための演算式である。数式(1)の記号Ｋ1は所定の比例定数であり、記号＜Ｉ^２＞は、検出信号Ｓdの全帯域にわたるパワーである。

数式(1)から理解される通り、第１実施形態の演算処理部６１は、強度スペクトルＸにおける各周波数ｆの強度Ｐ(f)と当該周波数ｆとの積（ｆ×Ｐ(f)）を所定の範囲（以下「演算範囲」という）について積算することで血流量Ｑを算定する。強度スペクトルＸの強度Ｐ(f)と周波数ｆとの積（ｆ×Ｐ(f)）は、周波数ｆにより重み付けされた強度（以下「加重強度」という）を意味する。演算範囲は、加重強度の積分範囲に相当し、周波数軸上の所定の周波数（以下「下端周波数」という）ｆ1と、下端周波数ｆ1を上回る所定の周波数（以下「上端周波数」という）ｆ2との間の範囲である。下端周波数ｆ1は第１周波数の例示であり、上端周波数ｆ2は第２周波数の例示である。以上の説明から理解される通り、第１実施形態の演算処理部６１は、検出装置３０が生成した検出信号Ｓdの強度スペクトルＸから血流情報（血流速度Ｖ）を生成する。図２の表示装置２４は、演算処理部６１が生成した血流情報（血流速度Ｖ）を表示する。

図５は、第１実施形態の信号処理部５２が実行するフィルター処理の説明図である。具体的には、信号処理部５２が実行するフィルター処理の周波数特性（すなわち、周波数領域におけるゲインの分布）が図５には例示されている。図５から理解される通り、第１実施形態の信号処理部５２は、所定の周波数帯域（以下「処理帯域」という）Ｂ内で周波数が高い成分ほど抑制されるようにフィルター処理を実行する。例えばローパスフィルターまたはバンドパスフィルターが信号処理部５２として好適に利用される。具体的には、１次のアナログフィルター回路の所定個（単数または複数）を組合せることで信号処理部５２が構成される。

処理帯域Ｂは、信号処理部５２による検出信号Ｓbの抑制の度合（すなわち、フィルター処理の周波数特性）が６ｄＢ／Ｏｃｔ（オクターブ）以上となる範囲であり、所定の周波数Ｆxに対して高周波側の範囲である。すなわち、処理帯域Ｂ内では、検出信号Ｓを構成する各周波数ｆの成分が、周波数ｆに比例した度合（６ｄＢ／Ｏｃｔ）以上の度合で抑制される。例えば、処理帯域Ｂ内でフィルター処理の周波数特性が６ｄＢ／Ｏｃｔに設定された場合を想定すると、処理帯域Ｂ内の特定の周波数ｍ×ｆ（ｍは自然数）でのゲインは、処理帯域Ｂ内の周波数ｆでのゲインの１／ｍに設定される。周波数Ｆxは、検出信号Ｓbの抑制の度合が６ｄＢ／Ｏｃｔ以上となる範囲（すなわち処理帯域Ｂ）の下端に相当する。

図５に例示される通り、血流量Ｑの算定に適用される演算範囲Ｒと信号処理部５２によるフィルター処理の処理帯域Ｂとは、部分的に相互に重複する。具体的には、演算範囲Ｒにおける高周波側の部分と処理帯域Ｂにおける低周波側の部分とが相互に重複する。すなわち、処理帯域Ｂの下端の周波数Ｆxは、演算範囲Ｒの下端周波数ｆ1を上回り、かつ、演算範囲Ｒの上端周波数ｆ2以下である。図５においては、処理帯域Ｂの下端の周波数Ｆxが演算範囲Ｒの上端周波数ｆ2を下回る場合が例示されている。

ここで、上端周波数ｆ2は、Ａ/Ｄ変換部５３のナイキスト周波数（Ｆs／２）以下である必要がある。すなわち、演算範囲Ｒの上端周波数ｆ2は、処理帯域Ｂの下端の周波数ＦxとＡ/Ｄ変換部５３のナイキスト周波数（Ｆs／２）との間の数値に設定される（Ｆx≦ｆ2≦Ｆs／２）。例えば、処理帯域Ｂの周波数Ｆxが４５ｋＨｚであり、Ａ/Ｄ変換部５３が１００ｋＨｚのサンプリング周波数Ｆsで動作する構成を想定すると、上端周波数ｆ2は、４５ｋＨｚ（＝Ｆx）以上かつ５０ｋＨｚ（＝Ｆs／２）以下の範囲内の適切な数値（例えば５０ｋＨｚ）に設定される。他方、下端周波数ｆ1は、上端周波数ｆ2と比較して充分に小さい数値（例えば２００Ｈｚ程度）に設定される。以上の説明から理解される通り、第１実施形態の演算処理部６１は、信号処理部５２によるフィルター処理で抑制される処理帯域Ｂの一部を含む演算範囲Ｒについて加重強度（ｆ×Ｐ(f)）を積算することで血流量Ｑを算定する。

図６には、検出信号Ｓdの強度スペクトルＸが図示されている。受光部３２や出力回路３４等の回路要素に起因して不可避的に発生するショットノイズが検出信号Ｓdには含有される。ショットノイズは、周波数ｆの広い範囲にわたり均等に分布するホワイトノイズである。他方、測定部位Ｍを通過した光（すなわち本来の解析対象）に由来する信号成分の強度は、高周波側ほど低いという傾向がある。したがって、図６から理解される通り、強度スペクトルＸのうち高周波側の周波数帯域Ｎでは、測定部位Ｍを通過した光に対してショットノイズの影響が支配的となる。ショットノイズの影響が支配的となる周波数帯域Ｎの下端の周波数（以下「高ノイズ周波数」という）Ｆnは、測定部位Ｍにおける血管または血液の状態に応じて変動し得る。

図７および図８は、周波数軸上における加重強度（ｆ×Ｐ(f)）の分布である。図７は、信号処理部５２によるフィルター処理を省略した構成（以下「対比例」という）における加重強度の分布である。図８は、処理帯域Ｂ内の信号成分を抑制するフィルター処理を信号処理部５２が実行する第１実施形態における加重強度の分布である。

いま、数式(1)の演算により血流量Ｑを高精度に測定するという観点からは、演算範囲Ｒを広く確保する必要がある。すなわち、上端周波数ｆ2を高目に設定することが要求される。他方、強度スペクトルＸの強度Ｐ(f)に対する周波数ｆの乗算により、対比例では、図７から理解される通り、強度スペクトルＸにおける高周波側のショットノイズが強調される。したがって、上端周波数ｆ2を高目に設定した構成では、演算範囲Ｒの高周波側においてショットノイズが支配的となり、結果的に血流量Ｑの高精度な測定が阻害される。以上の問題を解決するために、上端周波数ｆ2を低目に設定した場合には、演算範囲Ｒが過度に狭くなる結果、血流量Ｑの測定の精度が却って低下する可能性がある。以上の説明から理解される通り、血流量Ｑを常に高精度に測定するためには、測定部位Ｍにおける血管または血液の状態に応じて上端周波数ｆ2を変化させる必要がある。例えば、高ノイズ周波数Ｆnが高い状況では上端周波数ｆ2を高目に設定し、高ノイズ周波数Ｆnが低い状況では上端周波数ｆ2を低目に設定するといった具合である。

以上に説明した対比例とは対照的に、第１実施形態では、検出信号Ｓbのうち演算範囲Ｒに部分的に重複する処理帯域Ｂ内の信号成分がフィルター処理により抑制される。すなわち、図８から理解される通り、上端周波数ｆ2の高低に関わらず、演算範囲Ｒ内については、強度Ｐ(f)に対する周波数ｆの乗算により強調されるべきショットノイズの影響がフィルター処理で低減される。したがって、例えば演算範囲Ｒが充分に確保されるように上端周波数ｆ2を高目に設定した構成でも、検出信号Ｓbのショットノイズの影響を低減して血流量Ｑを高精度に測定することが可能である。また、測定部位Ｍにおける血管または血液の状態に応じて上端周波数ｆ2を変化させる前述の処理が不要であるから、血流情報の生成のための負荷が低減されるという利点もある。

次に、演算範囲Ｒと処理帯域Ｂとが部分的に相互に重複する前述の状況で成立する条件について説明する。いま、血流解析装置１００が測定可能な仕様上の血流速度Ｖの範囲が最小値Ｖ1から最大値Ｖ2までの範囲である場合を想定する。前述の通り、血液中の赤血球での光反射に起因した周波数シフト量Δｆは血流速度Ｖに比例する。具体的には、周波数シフト量Δｆは以下の数式(2)で表現される。

数式(2)の記号λは発光部３１が測定部位Ｍに照射する光の波長であり、記号θは、発光部３１から測定部位Ｍに入射する光の入射角である。実際の血流解析装置１００を想定すると、波長λは発光部３１が出射する光の波長として既知であり、入射角θは測定部位Ｍの表面に対する発光部３１の光軸の角度から確定される。また、記号ｎは測定部位Ｍ（特に動脈および血液）の屈折率であり、概略的には１．３３～１．３４の範囲内の既知の数値である。血流解析装置１００による測定が想定される血流速度Ｖの最大値Ｖ2を以上の各定数（λ，θ，ｎ）とともに数式(2)に代入すると、血流解析装置１００が測定可能な範囲内で最大の周波数シフト量Δｆ2を得ることができる。

なお、血流速度Ｖの最大値Ｖ2は、超音波を利用した流速計により測定することができる。また、動脈の血流を測定する場合には、例えば、血流速度Ｖの最大値Ｖ2は０．８ｍ／ｓｅｃ以上かつ１．２ｍ／ｓｅｃ以下であり、毛細血管の血流を測定する場合には、例えば、血流速度Ｖの最大値Ｖ2は２ｍｍ／ｓｅｃ以上かつ１２ｍｍ／ｓｅｃ以下であることが知られている。

次に、処理帯域Ｂの下端周波数Ｆxを決定する方法を説明する。図９は、処理帯域Ｂの下端周波数Ｆｘを特定する方法を説明するための図である。図９に例示される通り、スペクトラムアナライザーＳＡの端子Ｉから例えば正弦波の入力信号を信号処理部５２に入力し、信号処理部５２からの電力または電圧を出力信号として端子Ａに入力する。ここで、入力信号の周波数を変化させ、各周波数に対する出力信号を取得することで、図５から図８に例示したような周波数特性を得ることができる。処理帯域Ｂの下端周波数Ｆｘは、入力信号の周波数が高くなるにしたがって出力信号が低下する周波数である。具体的には、図１０に示すように、入力信号の周波数が高くなるにしたがって出力信号が３ｄＢだけ低下する周波数を処理帯域Ｂの下端周波数Ｆｘとする。なお、図９では、信号処理部５２の入力信号と出力信号との関係から周波数特性を得る場合を例示したが、信号増幅部５１に入力信号を入力してもよいし、信号処理部５２よりも後段の要素（例えばＡ／Ｄ変換部５３）から出力信号を取得してもよい。

また、より多くの血流情報を得るには、血流速度Ｖの最大値Ｖ2までの血流情報を得ることが望ましい。換言すると、上端周波数ｆ2が周波数シフト量Δｆ2以上である（ｆ2≧Δｆ2）ことが望ましい。この場合、周波数シフト量Δｆ2が周波数Ｆxを上回るという関係（Ｆx＜Δｆ2）が確認できれば、上端周波数ｆ2が周波数シフト量Δｆ2以上である（Δｆ2≦ｆ2）ので、上端周波数ｆ2が周波数Ｆxを上回る（Ｆx＜ｆ2）。上端周波数ｆ2が周波数Ｆxを上回る（Ｆx＜ｆ2）ということは、演算範囲Ｒと処理帯域Ｂとが部分的に重複するという条件が成立していることになる。

＜第２実施形態＞
本発明の第２実施形態を説明する。なお、以下に例示する各形態において作用または機能が第１実施形態と同様である要素については、第１実施形態の説明で使用した符号を流用して各々の詳細な説明を適宜に省略する。

図１１は、第２実施形態における受光部３２および出力回路３４の構成図である。図１１に例示される通り、第２実施形態の受光部３２は、相異なる位置に配置された受光素子３２１と受光素子３２２とを含んで構成される。受光素子３２１は、測定部位Ｍからの受光強度に応じた検出信号Ｓa1を生成し、受光素子３２２は、測定部位Ｍからの受光強度に応じた検出信号Ｓa2を生成する。

第２実施形態の信号増幅部５１は、受光素子３２１が生成した検出信号Ｓa1と受光素子３２２が生成した検出信号Ｓa2との差分に相当する検出信号Ｓbを生成する。したがって、検出信号Ｓa1と検出信号Ｓa2とに共通に含まれる定常的なノイズを低減した検出信号Ｓbが生成される。例えば差動増幅回路が信号増幅部５１として好適に利用される。信号処理部５２は、信号増幅部５１から供給される検出信号Ｓbに対して第１実施形態と同様のフィルター処理を実行する。他の要素の機能および作用は第１実施形態と同様である。

第２実施形態においても第１実施形態と同様の効果が実現される。また、第２実施形態では、受光素子３２１が生成した検出信号Ｓa1と受光素子３２２が生成した検出信号Ｓa2との差分に相当する検出信号Ｓbが生成される。すなわち、検出信号Ｓa1と検出信号Ｓa2とに共通に含まれるノイズを低減したＳＮ比が高い検出信号Ｓbが生成される。したがって、血流情報を高精度に生成できるという効果は格別に顕著である。

＜変形例＞
以上に例示した各形態は多様に変形され得る。具体的な変形の態様を以下に例示する。以下の例示から任意に選択された２以上の態様を適宜に併合することも可能である。

（１）相互に別体で構成された複数の機器により血流解析装置１００を実現することも可能である。例えば、前述の各形態で例示した演算処理部６１を、携帯電話機またはスマートフォン等の汎用の情報端末により実現することも可能である。また、演算処理部６１が生成した血流情報を、情報端末が具備する表示装置２４に表示させる構成も採用され得る。

（２）出力回路３４を構成する複数の要素の順序は前述の各形態の例示に限定されない。例えば、前述の各形態では、信号処理部５２が生成した検出信号ＳcをＡ/Ｄ変換部５３がＡ/Ｄ変換したが、例えば図１２に例示される通り、信号処理部５２とＡ/Ｄ変換部５３との順序を逆転することも可能である。図１２の構成では、信号増幅部５１が増幅した検出信号ＳbをＡ/Ｄ変換部５３がアナログからデジタルに変換し、変換後の検出信号Ｓcに対するフィルター処理で信号処理部５２が検出信号Ｓdを生成する。したがって、検出信号Ｓcのうち処理帯域Ｂ内の成分を抑制するデジタルフィルターが信号処理部５２として利用される。また、制御装置２０がプログラムを実行することで信号処理部５２を実現することも可能である。

（３）前述の各形態では、血流速度Ｖを血流情報として例示したが、血液の速度に関する情報（血流情報）の種類は以上の例示に限定されない。例えば、前述の数式(1)で算定される血流量Ｑを血流情報として被験者に提示することも可能である。また、血流量Ｑおよび血流速度Ｖ等の血流情報から他の生体情報を生成することも可能である。例えば、血圧、酸素飽和度（ＳｐＯ２）または血管年齢（血管の硬さ）等の各種の生体情報が、血流量Ｑおよび血流速度Ｖ等の血管情報から推定され得る。

（４）前述の形態では、筐体部１２とベルト１４とから構成される腕時計型の血流解析装置１００を例示したが、血流解析装置の具体的な形態は任意である。例えば、被験者の身体に貼付可能なパッチ型、被験者の耳介に装着可能なイヤリング型、被験者の指先に装着可能な指装着型（例えば着爪型）、または、被験者の頭部に装着可能なヘッドマウント型など、任意の形態の血流解析装置１００が採用され得る。

（５）前述の各形態では、被験者の血流を解析する血流解析装置１００を例示したが、本発明が適用される範囲は血流の解析に限定されない。例えば、血液以外の各種の液体（例えば管内を流動する薬液）の流動を解析する装置にも本発明を適用することが可能である。以上の説明から理解される通り、本発明の好適な態様は、流体を解析する装置（流体解析装置）であり、前述の各形態で説明した血流解析装置１００は、本発明の好適な態様に係る流体解析装置の例示である。

１００…血流解析装置、１２…筐体部、１４…ベルト、２０…制御装置、２２…記憶装置、２４…表示装置、３０…検出装置、３１…発光部、３２…受光部、３２１，３２２…受光素子、３３…駆動回路、３４…出力回路、５１…信号増幅部、５２…信号処理部、５３…Ａ/Ｄ変換部、６１…演算処理部。

Claims (4)

1. 流体を通過した光の強度を表す検出信号に対し、所定の処理帯域内で周波数が高い成分ほど抑制されるようにフィルター処理を実行する信号処理部と、
   前記フィルター処理後のアナログ信号を所定のサンプリング周波数でデジタル信号に変換するＡ/Ｄ変換部と、
   前記デジタル信号の強度スペクトルにおける各周波数の強度と当該周波数との積を、前記処理帯域に部分的に重複する所定の演算範囲について積算することで、前記流体の速度に関する情報を生成する演算処理部とを具備し、
   前記演算範囲は、第１周波数と前記第１周波数を上回る第２周波数との間の範囲であり、
   前記第２周波数は、前記処理帯域の下端の周波数を上回り、かつ、前記Ａ/Ｄ変換部の前記サンプリング周波数に対応するナイキスト周波数以下である
   流体解析装置。
2. 前記処理帯域は、前記信号処理部による前記検出信号の抑制が６ｄＢ／Ｏｃｔ以上となる範囲である
   請求項１の流体解析装置。
3. 血管を通過した光の強度を表す検出信号に対し、所定の処理帯域内で周波数が高い成分ほど抑制されるようにフィルター処理を実行する信号処理部と、
   前記フィルター処理後のアナログ信号を所定のサンプリング周波数でデジタル信号に変換するＡ/Ｄ変換部と、
   前記デジタル信号の強度スペクトルにおける各周波数の強度と当該周波数との積を、前記処理帯域に部分的に重複する所定の演算範囲について積算することで、前記血管内の血液の速度に関する情報を生成する演算処理部とを具備し、
   前記演算範囲は、第１周波数と前記第１周波数を上回る第２周波数との間の範囲であり、
   前記第２周波数は、前記処理帯域の下端の周波数を上回り、かつ、前記Ａ/Ｄ変換部の前記サンプリング周波数に対応するナイキスト周波数以下である
   血流解析装置。
4. 血管を通過した光の強度を表す検出信号に対し、所定の処理帯域内で周波数が高い成分ほど抑制されるようにフィルター処理を実行し、
   前記フィルター処理後のアナログ信号を所定のサンプリング周波数でデジタル信号にＡ/Ｄ変換し、
   前記デジタル信号の強度スペクトルにおける各周波数の強度と当該周波数との積を、前記処理帯域に部分的に重複する所定の演算範囲について積算することで、前記血管内の血液の速度に関する情報を生成し、
   前記演算範囲は、第１周波数と前記第１周波数を上回る第２周波数との間の範囲であり、
   前記第２周波数は、前記処理帯域の下端の周波数を上回り、かつ、前記サンプリング周波数に対応するナイキスト周波数以下である
   流体解析方法。
   

Images