(57)【要約】 【課題】  血圧測定を必要としない、手首や指先等の部位で測定可能な小型の血液レオロジー測定装置及び血液レオロジー測定方法を提供する。【解決手段】  生体表面から生体内の動脈血流に対する超音波の送受信を行う超音波送受信素子からなるセンサを用いて、時間変化する血流情報を超音波の送受信で得られるドップラー信号を観測する。さらに、得られたドップラー信号からドップラーシフト量成分毎の強度を求め、ドップラーシフト量成分毎の強度を積分する信号強度を求める。信号強度が血流によって変化することを利用し、血流の最大時の信号強度と最小時の信号強度の変化割合を血液レオロジーの指標として演算する。【選択図】    図１
2006181268.tif 000002
【特許請求の範囲】
【請求項１】
生体表面から生体内の動脈血流に対して超音波を送信する発信素子と生体中で反射された反射波を受信する受信素子とからなるセンサ部と、前記発信素子を駆動する送信回路と前記受信素子で受信したドップラー信号を検波する受信検波回路とからなる回路部と、前記受信検波回路によって検波された前記ドップラー信号を演算する演算部と、前記演算部で得られた結果を出力する出力部とを備える血液レオロジー測定装置において、
前記演算部は、前記受信検波回路から出力する前記ドップラー信号を高速フーリエ変換して演算を行うＦＦＴ演算処理部と、該ＦＦＴ演算処理部から出力する演算結果に基づいて、血流速度に対応する信号波形を演算する血流速度波形演算処理部と、該血流速度波形演算処理部において演算された脈動毎の前記信号波形の最大値と最小値に対応する時刻における前ドップラーシフト量成分毎の信号強度を積算し、前記最大値に対応する時刻における第１の積算値と前記最小値に対応する時刻における第２の積算値との変化割合から血液レオロジーの指標値を算出し、前記出力部に出力する血液レオロジー指標演算処理部と、からなることを特徴とする血液レオロジー測定装置。
【請求項２】
前記センサ部は、前記発信素子と受信素子とを１対とする複数対の超音波センサからなり、前記演算部は、前記受信検波回路を介して得られたそれぞれの対の前記超音波センサの受信素子からの前記ドップラー信号に基づいて前記信号強度を積算し、それぞれの対の対応する所定時刻における前記積算値の平均値を前記第１の積算値及び第２の積算値として前記変化割合を算出し、血液レオロジーの指標値を算出することを特徴とする請求項１に記載の血液レオロジー測定装置。
【請求項３】
前記変化割合を心拍の拍動毎に求め、拍動毎の前記変化割合の平均値から血液レオロジー指標値を算出することを特徴とする請求項１に記載の血液レオロジー測定装置。
【請求項４】
前記変化割合は、前記第２の積算値を前記第１の積算値で除算することにより算出することを特徴とする請求項１〜３のいずれかに記載の血液レオロジー測定装置。
【請求項５】
前記変化割合は、前記第１の積算値と前記第２の積算値との差を前記第１の積算値で除算することにより算出することを特徴とする請求項１〜３のいずれかに記載の血液レオロジー測定装置。
【請求項６】
血液の流動性を表す血液レオロジーを測定する血液レオロジー測定方法であって、
生体表面から生体内の動脈血流に対して超音波を送受信し、受信したドップラー信号を検波する第１の工程と、前記ドップラー信号を高速フーリエ変換して血流速度に対応する信号波形を演算する第２の工程と、脈動毎の前記信号波形の最大値と最小値に対応する時刻におけるドップラーシフト量成分毎の信号強度を積算する第３の工程と、前記最大値に対応する時刻における第１の積算値と前記最小値に対応する時刻における第２の積算値との変化割合から血液レオロジーの指標値を算出する工程と、からなることを特徴とする血液レオロジー測定方法。
【請求項７】
発信素子と受信素子を１対とする複数対のセンサを用いて前記第１から第３の工程までを行い、前記第４の工程において、それぞれの対の対応する所定時刻における前記積算値の平均値を前記第１の積算値及び第２の積算値とすることを特徴とする請求項6に記載の血液レオロジー測定方法。
【請求項８】
前記第４の工程において、前記変化割合を心拍の拍動毎に求め、拍動毎の前記変化割合の平均値から血液レオロジー指標値を算出することを特徴とする請求項１に記載の血液レオロジー測定方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、一般に血液のサラサラ／ドロドロ度と称される流動性を表す血液レオロジー測定装置及び血液レオロジー測定方法にかかわり、特に動脈を流れる血流を測定し、人体組織の活動の基になる微小循環血流を見極め、健康の評価、疾患の診断、薬品の効果の評価等を行う技術に関する。
【背景技術】
【０００２】
人間の健康の評価、疾患の診断、人体への薬品の効果の評価、食品の健全性・機能性の評価等を行うために、血液レオロジーを計測して、その結果から前記の評価や診断を行っていくことが従来から行われている。従来の技術としては、血液流動性測定装置として被検者から血液を採取し、リソグラフィックな手法で製作されたマイクロチャネルアレイを用いて、定圧下の血流の通過時間から血液レオロジーを計測する方法が知られている。（非特許文献１参照）
しかし、従来のようなマイクロチャネルアレイを用いた血液レオロジー測定法では、どうしても被検者から血液を採取するために、注射針を用いて肘部に針を刺し、採血を行わなければならない。従って、食品成分の血液レオロジーへの影響をみるためのin vitro試験を行うとしても、同じ人から１日何回も血液採取を行うことができず、連続試験が困難であるという問題がある。また、医療機関を離れて個人が自宅等で自ら採血をして血液レオロジー測定を行おうとしても、従来例のような血液流動性測定装置を用いた一方法では、自宅に機器を置くこともできず、適切な処理もできないため、医療機関でしか測定ができないという問題もあった。
【０００３】
ところで、血液レオロジーと生体内の血流速度は強い相関があると考えられている。すなわち、血液の粘性が高い場合、血流速度は遅く、一方、粘性が低い場合は血流速度が速いと考えられている。そのため、生体内の血流速度を計測することで、間接的に血液レオロジーを知ることが可能となる。（例えば、特許文献１参照。）
一方、血管内の血流速度から血液レオロジーの指標を算出するためには、前記特許文献１に記載されているように、血流速度の計測以外に、カフを用いて生体の血圧を測定する必要があるが、この血圧値と血流速度をもちいて血液レオロジー、すなわち、血液の運動粘性率の指標を算出する方法として、対象とする動脈内部の血流圧力を血圧値で近似するといった概念に基づく方法がある。
【特許文献１】特開２００３−１５９２５０号公報  「血液レオロジー測定装置」
【非特許文献１】菊池佑二「毛細血管モデルを用いた全血流動性の測定」（食品研究成果情報，NO.11  1999年発行）
【発明の開示】
【発明が解決しようとする課題】
【０００４】
しかしながら、血圧値と血流速度をもちいて血液レオロジー、すなわち、血液の運動粘性率の指標を算出する方法は、血圧測定のための機構及び煩雑さ等の観点から、手首や指先等の部位で血液レオロジーを測定するために必要不可欠な装置の小型化が困難であるという問題があった。
【０００５】
本発明の目的は、手首や指先等の部位で測定可能あり、血圧測定を必要としない小型の血液レオロジー測定装置及び血液レオロジー測定方法を提供する事である。
【課題を解決するための手段】
【０００６】
本発明に係る血液レオロジー測定装置、及び血液レオロジー測定方法は、上記課題を解決するために、生体表面から生体内の動脈血流に対する超音波の送受信を行う超音波送受信素子からなるセンサを用いて、時間変化する血流情報を超音波の送受信で得られるドップラー信号によって観測する。さらに、得られたドップラー信号からドップラーシフト量成分毎の強度を求め、ドップラーシフト量成分毎の強度を積分する信号強度を求める。信号強度が血流によって変化することを利用し、血流の最大時の信号強度と最小時の信号強度の変化割合を血液レオロジーの指標として、演算する手段を備えることによって、上述の課題を解決するものである。
信号強度の変化割合とは、以下に記すものである。
送信信号強度をIoとおき、ある時間ｔでの受信信号強度Iｔとおくと、I０とIｔの関係は式1で表せる。
【０００７】
【数１】
2006181268.tif 000003
ここで、F(x)は主に皮膚・生体内を伝播するときの減衰項、Gt(y)は主に血流からの反射項である。血流の最大時をImax、最小時をIminとすると、それぞれ式2、式3で表せる。
【０００８】
【数２】
2006181268.tif 000004
【０００９】
【数３】
2006181268.tif 000005
式2・式3では減衰項の項があるため、ImaxやIminと他の血液レオロジー指標との相関をとると、減衰項の個人差が影響するために、良い相関は得られない。
そこで、Itの変化割合I ’を求める式を式4のようにすることで、F(x)を除算でき、個人差による減衰項の違いを除去できる。
【００１０】
【数４】
2006181268.tif 000006
Gt(y)の成分としては、赤血球からの反射成分をGt(r)、白血球からの反射成分をGt(w)、血小板からの反射成分をGt(p)、骨・筋肉等の反射成分をGt(o)とすると、Gt(y)は式5で表せる。
【００１１】
【数５】
2006181268.tif 000007
となる。Itの信号の中で、ドップラーシフト成分を抽出する場合、骨等からの反射にドップラーシフトは発生しないとみなせるので、式6が成立つ。
【００１２】
【数６】
2006181268.tif 000008
また、血流中で最も多くしめるものは赤血球であり、白血球・血小板は赤血球に比較して、非常に少ないので、式7が成立つ。
【００１３】
【数７】
2006181268.tif 000009
従って、式4に式5から式7を代入して、整理すると、変化割合I ’は、式8となる。
【００１４】
【数８】
2006181268.tif 000010
Itの変化割合I ’を求めることは、式8から、赤血球の移動量に応じたもの（速度分布）の変化割合を表すこととなる。血液の粘性が高ければ、速度分布の変化が小さくなり、血液の粘性が低ければ速度分布の変化は大きくなる。ドップラー信号をFFT（Fast  Fourier  Transform：高速フーリエ変換）解析し、ドップラーシフト量成分の信号強度を積分して求めることは、ドップラー信号が主に赤血球の移動によって生じていることから、赤血球の速度分布Gt(r)を求める一つの表現方法である。
従って、信号強度の変化割合を求めることは、赤血球速度分布の成分の変化割合を求めることであり、血液の粘性が高ければ、信号強度の変化割合の値は小さくなり、血液の粘性が低ければ、信号強度の変化割合の値は大きくなる。
【発明の効果】
【００１５】
図８は、本発明の効果を示す特性図であって、本発明に係る血液レオロジー装置にて測定した血液レオロジーの指標値と血液流動性測定装置を用いた採血方式による血液レオロジーの指標である全血通過時間の相関を示している。詳細は後述するが、縦軸は本発明での指標値を示しており、指標値が小さいということは血液の粘性が大きいということを意味する。横軸は全血通過時間を示しており、全血通過時間の値が小さいということは、血液の粘性が小さいことを意味する。即ち血液の粘性が大きくなると、本発明の指標値では小さくなり、全血通過時間は大きくなるということである。これらの関係から、図８は、本発明での指標値と全血通過時間とは有意な相関を有しているといえる。
【００１６】
従って、本発明に係る血液レオロジー測定装置は、血圧測定を必要とせずに手首や指先で血液レオロジーを測定する事が可能となり、小型の血液レオロジー測定装置が供給できる。その結果、被験者から採血を行うことなく、専門家以外の誰でも手軽に正確なレオロジーを調べることができ、健康状態の確認に利用することができるようになる。
【発明を実施するための最良の形態】
【００１７】
図１に、本発明に係る血液レオロジー測定装置の構成を示すブロック図を示す。
センサ部１は、２対の超音波センサ１１及び１２と圧力センサ１３と温度センサ１４から構成されている。超音波センサ１１は発信素子１１１と受信素子１１２、超音波センサ１２は発信素子１２１と受信素子１２２から構成されている。これらの送信素子及び受信素子はPZTなどの圧電素子が適している。
【００１８】
センサ部１において、超音波センサを2対用いる理由であるが、２対の超音波センサを血流に対して異なる角度に配置することにより、生体中の見えない血管の流れる方向を特定し、接触位置によらず安定し、高精度に測定が可能になるからである。
【００１９】
回路部２は発信素子１１１及び１２１を駆動する送信回路２１と、受信素子１１２及び１２２からの受信信号を受信し、ドップラー信号を検波する受信検波回路２２と、圧力センサ１３にかかる荷重を検出する検出回路２３と、温度センサ１４が感知する温度を検出する検出回路２４から構成されている。受信検波回路２２と、検波回路２３と、検波回路２４から得られる信号は演算部３に送られる。
【００２０】
演算部３はFFT（Fast  Fourier  Transform：高速フーリエ変換）演算処理部３１と血流速度波形演算処理部３２と血液レオロジー指標演算処置部３３から構成される。FFT演算処理部３１は受信検波回路２２で得られるドップラー信号に対してFFT演算を行い、その結果を血流速度波形演算処理部３２及び血液レオロジー指標演算処置部３３に送る。血流速度波形演算処理部３２ではFFT演算処理部３１の演算結果を更に演算処理することで血流速度及び血流速度波形を求める。血液レオロジー指標演算処置部３３ではFFT演算処理部３１の演算結果を更に演算処理することで血液レオロジーを求める。これらの結果は、出力部４に送られる。
出力部４はモニター或いはスピーカーであり、演算結果をモニター上に表示或いは音により表現する。
【００２１】
図２に本発明での測定方法を模式図で示す。センサ部１は本体の筐体に組み込まれており、図示しないが、フレキシブル配線基板等によりセンサ部１と回路部２は接続されている。図２に示すとおり、超音波センサ１１及び超音波センサ１２の角度及び間隔などの配置位置は、超音波センサ１１によって送受信される超音波と、超音波センサ１2によって送受信される超音波が交差する部分と、測定する動脈血管５２が重なるように設計することで、様々な測定部位での正確な測定が可能となる。
【００２２】
温度センサ１４は測定部位の温度周辺を測定するために設けている。図２では、２つの超音波センサ間に埋め込むことで、小型化している。本方法では小型化を実現するために、温度センサ１４にはチップサーミスタを用いている。また、圧力センサ１３上に2つの超音波センサを構成することで測定部位にかかる圧力を測定できる。圧力センサ１３及び温度センサ１４を備える理由は、測定が適正に行われているかを判断するためのものであり、その判断は血流速度波形演算処理部３２及び血液レオロジー指標演算処置部３３で行われる。
【００２３】
図２に示すように測定部位の皮膚５１とセンサ部１を密着させて測定すると、血流速度波形演算処理部３２での演算結果として図３に示す血流速度波形がえられる。図３は横軸が時間であり、縦軸が血流速度である。心臓は収縮と弛緩を周期的に行っているので、心拍に対応した血流速度の変化が測定される。血流速度は、心臓が収縮するときに上昇し、心臓が弛緩するときに下降するので、心拍毎に血流速度のピークとボトムが得られ、図３において血流速度のボトムが時刻T1・T3・T5に表れ、血流速度のピークが時刻T2・T4・T6に表れる。
【００２４】
このとき、血流波形のボトム時とピーク時即ち図３のT1及びT2の時のドップラー信号に対して、FFT演算処理部３１では図４から図７が得られる。図4から図７において、横軸はドップラーシフト量、縦軸は信号強度であり、図４は時刻T1のときの受信素子１１２から得られた結果であり、図５は時刻T1のときの受信素子１２２から得られた結果であり、図６は時刻T2のときの受信素子１１２から得られた結果であり、図７は時刻T2のときの受信素子１２２から得られた結果である。
【００２５】
図４において図に示すようにドップラーシフト量成分の信号強度を積分すると、積分値S1が得られる。同様に図５から図７においても積分値S2からS4が得られる。ここで、血流速度のボトム時の信号強度をS1とS2の平均値ST1を求め、ピーク時の信号強度をS3とS4の平均値ST2を求める。ST1をSminとし、ST2をSmaxとし、SminとSmaxの変化割合ηを血液レオロジーの指標値とすると、採血方式の血液レオロジー指標である全血通過時間と相関がある。測定は心拍1拍分だけを行うのでなく、数秒間行うので、心拍毎のピーク時とボトム時の信号強度を求め、更にピーク時毎に平均した値をS maxとし、ボトム時毎に平均した値をSminとする。図３を用いて説明しなおすと、血流速度のボトム時としてT1・T3・T5が観測されているので、それぞれの信号強度をST1・ST3・ST5と求める。このとき、Sminは、式9で表される。
【００２６】
【数９】
2006181268.tif 000011
同様にSmaxは、式10で表される。
【００２７】
【数１０】
2006181268.tif 000012
となる。本発明での血液レオロジーの指標値、即ち、変化割合ηを求める式としては、式11とした。
【００２８】
【数１１】
2006181268.tif 000013
変化割合ηを求める式としてはこれに限らなく、例えば式12でもよい。
【００２９】
【数１２】
2006181268.tif 000014
これらの演算を血液レオロジー指標演算処置部３３で行う。
【００３０】
図８に本発明の方式と採血方式との相関を示す。図8において横軸は採血方式の血液レオロジーの指標である全血通過時間であり、縦軸は本発明での血液レオロジーの指標値である。測定は本発明の方式で測定した直後に採血方式での測定を行っている。また、本発明の方式では指尖部位の血流を測定した。ドップラーシフト量成分の信号強度の積分範囲は、0Hz〜5kHzとした。この範囲は測定部位や、発信素子の駆動周波数によって、適宜変える必要がある。また、測定時間は数秒間行うので、心拍毎に血流速度のボトム時の信号強度とピーク時の信号強度を求め、それらを平均してSmin・Smaxを求めた。
信号強度の変化割合ηと全血通過時間の相関については、以下のように、説明できる。
送信信号強度をIoとおき、ある時間ｔでの受信信号強度Itとおくと、IoとItの関係は式1で表せる。
【００３１】
【数１３】
2006181268.tif 000015
ここで、F(x)は主に皮膚・生体内を伝播するときの減衰項、Itは主に血流からの反射項である。血流の最大時をImax、最小時をIminとすると、それぞれ式2、式3で表せる。
【００３２】
【数１４】
2006181268.tif 000016
【００３３】
【数１５】
2006181268.tif 000017
式2・式3では減衰項の項があるため、ImaxやIminと他の指標との相関をとると、減衰項の個人差が影響するために、良い相関は得られない。
そこで、Itの変化割合I ’を求める式を式4もしくは式13のようにすることで、F(x)を除算でき、個人差による減衰項の違いを除去できる。
【００３４】
【数１６】
2006181268.tif 000018
【００３５】
【数１７】
2006181268.tif 000019
Gt(y)の成分としては、赤血球からの反射成分をGt(r)、白血球からの反射成分をGt(w)、血小板からの反射成分をGt(p)、骨・筋肉の反射成分をGt(o)とすると、Gt(y)は式5で表せる。
【００３６】
【数１８】
2006181268.tif 000020
となる。Itの信号の中で、ドップラーシフト成分を抽出する場合、骨等からの反射にドップラーシフトは発生しないとみなせるので、式6が成立つ。
【００３７】
【数１９】
2006181268.tif 000021
また、血流中で最も多くしめるものは赤血球であるので、式7が成立つ。
【００３８】
【数２０】
2006181268.tif 000022
従って、式4もしくは式13に式5から式7を代入して、整理すると、変化割合I ’は、式8もしくは式14となる。
【００３９】
【数２１】
2006181268.tif 000023
【００４０】
【数２２】
2006181268.tif 000024
Itの変化割合I ’を求めることは、式8もしくは式14から、赤血球の移動量に応じたもの（速度分布）の変化割合を表すこととなる。血液の粘性が高ければ、速度分布の変化が小さくなり、血液の粘性が高ければ速度分布の変化は大きくなる。ドップラー信号をFFT解析し、ドップラーシフト量成分の信号強度を積分して求めるStは赤血球の速度分布Gt(r)の一つの表現方法である。従って、信号強度の変化割合ηを求めることは、赤血球速度分布の成分の変化割合を求めることであり、血液の粘性が高ければ、信号強度の変化割合ηは小さくなり、血液の粘性が小さければ信号強度の変化割合ηは大きくなる。
【００４１】
一方、血液流動性測定装置の全血通過時間は血液のサラサラ・ドロドロを示す指標となっており、血液がドロドロしている、即ち、血液の粘性が高いほど全血通過時間が大きくなることが知られている。従って、信号強度の変化割合ηと全血通過時間には、右下がりの相関がある。
【００４２】
ここまで、ドップラー信号をFFT処理することによって、血液レオロジーの指標値を求める方法について説明してきたが、速度分布を算出することでも血液レオロジーの指標値を求められる。図３に示す血流波形を求めるときに、同時に図９に示す速度分布を血流速度波形演算処理部３２で算出する。図９は時刻T1における速度分布を算出しており、横軸が速度、縦軸が、速度成分の信号強度を示す。この時、速度成分に対して積分を実施し、ST1 ’を求める。同様に時刻T2・T3・T4・T5・T6において、速度成分の信号強度を求め、先ほど同様にSmin及びSmaxを求めることが可能である。そして、SminとSmaxの変化割合を血液レオロジー指標演算処置部３３で演算することで、本発明での血液レオロジーの指標値とすることが可能である。
【産業上の利用可能性】
【００４３】
本発明は、医療および健康維持・増進を目的として、体液の流動性を示す指標としての血液レオロジーを採血することなく求めることができる。また、食品が血液に及ぼす影響を容易に把握できるようになるので、食品開発における検査装置として利用可能である。
【図面の簡単な説明】
【００４４】
【図１】血液レオロジー測定装置の構成を示すブロック図
【図２】測定状態の模式図
【図３】血流速度波形を示す図
【図４】時刻T1での受信素子１１２から得られるドップラーシフト成分と信号強度を示す図
【図５】時刻T1での受信素子１２２から得られるドップラーシフト成分と信号強度を示す図
【図６】時刻T2での受信素子１１２から得られるドップラーシフト成分と信号強度を示す図
【図７】時刻T2での受信素子１２２から得られるドップラーシフト成分と信号強度を示す図
【図８】本発明の効果を示す特性図
【図９】速度成分と強度を示す図
【符号の説明】
【００４５】
１  センサ部
１１  超音波センサ
１１１  発信素子
１１２  受信素子
１２  超音波センサ
１２１  発信素子
１２２  受信素子
１３  圧力センサ
１４  温度センサ
２  回路部
２１  超音波送信回路
２２  超音波受信検波回路
２３  圧力検波回路
２４  温度検波回路
３  演算部
３１  FFT演算処理部
３２  血流速度波形演算処理部
３３  血液レオロジー指標演算処置部
４  出力部
５  生体
５１  皮膚
５２  動脈血管
【図１】
2006181268.tif 000025
【図２】
2006181268.tif 000026
【図３】
2006181268.tif 000027
【図４】
2006181268.tif 000028
【図５】
2006181268.tif 000029
【図６】
2006181268.tif 000030
【図７】
2006181268.tif 000031
【図８】
2006181268.tif 000032
【図９】
2006181268.tif 000033