(57)【要約】 【課題】  手首や指先等の部位で測定可能であり、かつ血圧測定を必要としない簡便、高精度及び小型の血液レオロジー測定装置及び血液レオロジー計測方法を提供する。【解決手段】  生体表面から生体内の動脈血流に対する超音波の送受信を行う超音波送受信素子からなるセンサを用いて、時間変化する動脈の血流速度、拍動変位量、動脈径、動脈壁厚、心拍数（心拍周波数）、及び該血流速度と該拍動変位量の位相差または振幅比を検出し、前記位相差及び前記振幅比のうちいずれか一方と前記血管径と前記心拍周波数とを用いて演算を行い血液の運動粘性率を算出して血液レオロジーの指標値を得る【選択図】    図１
2006218169.tif 000002
【特許請求の範囲】
【請求項１】
生体表面から該生体内を測定し、前記生体内の動脈を流れる血液の血流速度に対応した血流速度波形を出力する第１の測定部と、同じく前記生体表面から該生体内を測定し、前記動脈の血管径の計測値、及び心臓の拍動運動にともなって伸縮する前記動脈血管径の変位量に対応した血管拍動変位波形を出力する第２の測定部と、前記第１の測定部と前記第２の測定部とから出力された前記血流速度波形と前記血管拍動変位波形と前記血管径の計測値から前記動脈内の血液レオロジーを算出するレオロジー演算処理装置とからなることを特徴とする血液レオロジー測定装置。
【請求項２】
前記レオロジー演算処理装置は、前記第１の測定部から出力された前記血流速度波形、または前記第２の測定部から出力された前記血管拍動変位波形から心拍周波数を算出する心拍数演算処理装置と、前記血流速度波形と前記血管拍動変位波形とからそれぞれの波形の位相差または振幅比を算出するパラメータ演算処理装置と、算出された前記位相差及び前記振幅比のうちいずれか一方と前記血管径と前記心拍周波数とを用いて演算を行い血液の運動粘性率を算出するレオロジー指標演算処理装置とを有することを特徴とする請求項１に記載の血液レオロジー測定装置。
【請求項３】
前記レオロジー指標演算処理装置は、前記運動粘性率と前記心拍周波数と前記血管径との対応関係を記した特性表または関係式を記憶する記憶部を有し、該記憶部に記憶した前記特性表または関係式を用いて、前記運動粘性率を算出することを特徴とする請求項２に記載の血液レオロジー測定装置。
【請求項４】
前記第１の測定部は、前記生体表面から前記生体内の動脈血流に対して超音波連続信号を送受信する第１の送受信部と、受信した超音波連続信号からドップラー偏移を受けたドップラーシフト信号を演算処理して前記前記血流速度波形を出力する血流速度演算処理装置とからなることを特徴とする請求項１〜３のいずれかに記載の血液レオロジー測定装置。
【請求項５】
前記第２の測定部は、前記生体表面から前記生体内の動脈対して超音波バースト信号を送受信する第２の送受信部と、受信した超音波バースト信号を演算処理して前記動脈の血管径の計測値及び前記血管拍動変位波形を出力する血管情報演算処理装置とからなることを特徴とする請求項１〜４のいずれかに記載の血液レオロジー測定装置。
【請求項６】
前記第１の送受信部は、前記超音波の出射方向および受信感度の指向方向が互いに平行でない角度に配置された２対の超音波センサ素子からなる第１のセンサ部と、前記第１のセンサ部を駆動し、受信した超音波を検波する第１の送受信回路部とからなることを特徴とする請求項１〜５のいずれかに記載の血液レオロジー測定装置。
【請求項７】
前記第２の送受信部は、前記生体内に対して前記超音波バースト信号を送受信する１対の超音波センサ素子からなる第２のセンサ部と、前記第２のセンサ部を駆動し、受信した超音波を検波する第２の送受信回路部とからなることを特徴とする請求項１〜６のいずれかに記載の血液レオロジー測定装置。
【請求項８】
前記第１のセンサ部を構成する前記２対の超音波センサ素子のうちの１対の超音波センサ素子が前記第２のセンサ部を形成し、該第２のセンサ部に接続される前記第１の送受信回路部及び前記第２の送受信回路部との切り替えを行う切り替え回路部を有することを特徴とする請求項７に記載の血液レオロジー測定装置。
【請求項９】
前記第１のセンサ部と前記第２のセンサ部が、同一基板上に配置されていることを特徴とする請求項６〜８のいずれかに記載の血液レオロジー測定装置。
【請求項１０】
超音波を用いて生体表面から生体内の動脈を流れる血液の血液レオロジーを計測する血液レオロジー計測方法であって、
前記生体内の動脈を流れる血液の血流速度に対応した血流速度波形を計測する工程と、前記動脈の血管径を計測する工程と、心臓の拍動運動にともなって伸縮する前記動脈の血管径の変位量に対応した血管拍動変位波形を計測する工程と、計測された前記血流速度波形と前記血管拍動変位波形と前記血管径の計測値から前記動脈内の血液レオロジーを算出する工程とからなることを特徴とする血液レオロジー計測方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、生体中を循環する体液の測定装置及び計測方法にかかわり、特に血液の状態を把握し健康の評価、疾患の診断、薬品の効果の評価等のために使用する血液レオロジー測定装置および血液レオロジー計測方法に関する。
【背景技術】
【０００２】
人体の健康状態を判断する検査項目のひとつとして、血液の流動性に着目した血液レオロジー測定が注目されている。血液レオロジーを測定する手段として、被験者より採血した一定量の血液が微小流路（マイクロチャネルアレイ）を通過する時間を計測するマイクロチャネルアレイ型血液流動性測定装置が開発されている（例えば、非特許文献１参照。）。現在においては、マイクロチャネルアレイ型血液流動性測定装置は、血液レオロジー測定における標準機とされている。
【０００３】
しかし、マイクロチャネルアレイ型血液流動性測定装置による測定においては上記のように必ず採血を行う必要があり、測定が行えるのは医療機関に限られ、いつでもだれでもが手軽に健康状態を検査するというわけにはいかない。また、採血は被験者に対する肉体的および心理的な負担も大きく、１日あたりに計測できる回数もせいぜい数回までが限界であるため、時系列的に連続したデータが容易に得られないという問題がある。
ところで、血液レオロジーと生体内の血流速度は強い相関があると考えられている。すなわち、血液の粘性が高い場合、血流速度は遅く、一方、粘性が低い場合は血流速度が速いと考えられている。そのため、生体内の血流速度を計測することで、間接的に血液レオロジーを知ることが可能となる（例えば、特許文献１参照。）。
一方、血管内の血流速度から血液レオロジーの指標を算出するためには、前記特許文献１に記載されているように、血流速度の計測以外に、カフを用いて生体の血圧を測定する必要があるが、この血圧値と血流速度をもちいて血液レオロジー、すなわち、血液の運動粘性率の指標を算出する方法として、対象とする動脈内部の血流圧力を血圧値で近似するといった概念に基づく方法がある。
【特許文献１】特開２００３−１５９２５０号公報
【非特許文献１】菊池佑二「毛細血管モデルを用いた全血流動性の測定」（食品研究成果情報，NO.11  1999年発行）
【発明の開示】
【発明が解決しようとする課題】
【０００４】
しかしながら、血圧値と血流速度をもちいて血液レオロジー、すなわち、血液の運動粘性率の指標を算出する方法は、動脈内部の血流圧力を血圧値で近似している事が原因で、計測誤差が大きいという問題があった。さらには、血圧測定のための機構及び煩雑さ等の観点から、手首や指先等の部位で血液レオロジーを測定するために必要不可欠な装置の小型化が困難であるという問題もあった。
【０００５】
そこで本願発明の目的は、手首や指先等の部位で測定可能かつ血圧測定を必要としない簡便、高精度及び小型の血液レオロジー測定装置及び血液レオロジー計測方法を提供する事である。
【課題を解決するための手段】
【０００６】
本発明は、上記課題を解決するために、生体表面から生体内の動脈血流に対する超音波の送受信を行う超音波送受信素子からなるセンサを用いて、時間変化する動脈の血流速度、拍動変位量、動脈径、動脈壁厚、心拍数（心拍周波数）、及び該血流速度と該拍動変位量の位相差または振幅比を検出し、前記位相差及び前記振幅比のうちいずれか一方と前記血管径と前記心拍周波数とを用いて演算を行い血液の運動粘性率を算出して血液レオロジーの指標値を得ることを特徴とするものである。
【発明の効果】
【０００７】
図１１は、本発明の効果を示す特性図であって、本発明に係る血液レオロジー装置にて測定した血流速度と拍動変位量の位相差から算出した血液レオロジーの指標値である運動粘性率νとマイクロチャネルアレイを用いた採血方式による血液レオロジーの指標である全血通過時間Tの相関を示している。縦軸はνを示しており、図１１を見て縦軸の原点に近いほうがνは小さく、上方に行くほどνは大きな値となる。詳細については後述するが、νが小さいということは血液の粘性率が大きいということを意味している。
【０００８】
一方、横軸は全血通過時間Tを示しており、図１１を見て縦軸の原点に近いほうがＴの値は小さく、図中、右方に行くほどＴは大きな値となる。すなわち、全血通過時間Tの値が小さいということは、測定している血液がさらさらの血液であることを意味し、すなわち、νが大きな値となる。一方、全血通過時間Tの値が大きいということは、測定している血液が粘性率の高いどろどろの血液であることを意味している。すなわち、粘性率が高いということはνが小さいということであり、これらの関係を考慮して図１１を見た場合、νと全血通過時間Tとは有意な相関を有しているとみなすことが可能である。
【０００９】
従って、この図１１からわかる通り、本発明に係る血液レオロジー測定装置は、血圧測定を必要せずに精度よく手首や指先で血液レオロジーを測定する事が可能となるので、簡便、高精度及び小型の血液レオロジー測定装置が供給でき、その結果、被験者から採血を行うことなく、専門家以外の誰でも手軽に正確なレオロジーを調べることができ、健康状態の確認に利用することができるようになる。
【発明を実施するための最良の形態】
【００１０】
図１に、本発明に係る血液レオロジー測定装置の構成を示すブロック図を示す。センサ部は、２対の超音波センサすなわち、超音波センサ１と超音波センサ４と前記2対の超音波センサとは独立な超音波センサ７より構成されている。超音波センサ１は発信素子２と受信素子３、超音波センサ４は発信素子５と受信素子６から構成されており、超音波センサ７は発信素子８と受信素子９から構成されている。発信素子２及び発信素子５は連続超音波発生回路１０と接続しており、連続超音波発生回路１０にて発生する電気信号を機械的超音波に変換して、生体内に超音波を発信させる。
【００１１】
動脈内の血流に反射してドップラー信号を伴った超音波信号は、受信素子３及び受信素子６で電気信号に変換され、連続超音波検波回路１１に入力され、ドップラー電気信号が検波される。この連続超音波発生回路１０と連続超音波検波回路１１の２種類の回路で超音波回路１２が構成されている。バースト発生回路１３は接続された発信素子８に対して電気的バースト信号を出力して、発信素子８を駆動し、該発信素子８は前記電気的バースト信号を超音波バースト信号に変換し、その超音波バーストを生体内に照射する。この照射された超音波バースト信号は、動脈にて反射して受信素子９にて電気信号に変換された後にバースト検波回路１４に入力され、バースト検波回路１４において動脈による反射バースト超音波が電気信号として検波される。このバースト発生回路１３とバースト検波回路１４の２種類の回路で超音波バースト回路１５が構成されている。
【００１２】
連続超音波検波回路１１にて検波されたドップラー電気信号は、生体心拍と同期した周期的時間変化を伴う血流速度の信号成分を含んでいる。このドップラー電気信号から血流速度に対応する電気信号、即ち、血流速度信号を分離抽出するための装置が血流速度演算処理装置１６である。また、バースト検波回路１４にて検波された反射バースト超音波に対応する電気信号は、生体心拍と同期した周期的時間変化を伴う動脈の拍動変位と形状情報を含んでいる。電気信号から動脈の拍動変位量及び形状情報を分離抽出するための装置が血管情報演算処理装置１７である。
【００１３】
血流速度演算処理装置１６から出力された血流速度信号と、血管情報演算処理装置１７から出力された血管情報信号は、本発明に係るレオロジー演算処理装置１８に入力される。このレオロジー演算処理装置１８は、波形情報演算処理装置１９、心拍周波数演算処理装置２０、パラメータ演算処理装置２１及びレオロジー指標演算処理装置２２より構成されている。このレオロジー演算処理装置１８にて検出されたレオロジー指標が出力装置２３にて出力される。以上、図１は血流速度を測定する超音波センサ部と血管情報を測定する超音波センサ部が独立に配置されている実施形態である。
【００１４】
図２に、本発明に係る第二の血液レオロジー測定装置の構成を示すブロック図を示す。センサ部は、２対の超音波センサすなわち、超音波センサ１と超音波センサ２４で構成されている。超音波センサ２４を構成している発信素子２５と受信素子２６は図１記載の超音波回路１２と超音波バースト回路１５の双方に切り替え回路２７を介して接続されている。この切り替え回路２７は、超音波センサ２４の接続回路を周期的に超音波回路１２と超音波バースト回路１３に切り替える。
すなわち、超音波センサ２４は、血流速度検出と血管情報検出の双方の機能を持つ事になる。血流速度検出は超音波周波数のドップラー偏移量測定、血管情報検出は反射バースト波の遅延時間等の測定に基づいて行われる。それゆえ、この図２記載の切り替え回路２７における切り替えの周期の設定は、これらの測定物理量が測定可能な時間幅であれば良いので、単なる設計的事項にすぎない。この図２は、血流速度を測定する超音波センサ部が、血管情報を測定するセンシング機能も兼ね備えている実施形態である。
【００１５】
以上の実施形態においては、血流速度を測定するための超音波センサを複数個使用した場合について説明したが、特に複数個の超音波センサを使用することに拘るものではなく、例えば1個の超音波センサ、即ち、発信素子２と受信素子３とからなる超音波センサ１のみであってもよい。しかし、好ましくは、超音波センサを本実施形態のように２個使用することが望ましい。それは、超音波センサを本実施形態（図８及び図９）のように、超音波の出射方向および受信感度の指向方向が互いに平行でない角度に配置された２個の超音波センサとして使用すると、見えない血管の流れる方向を特定し、例えば、指の接触位置によらず安定して、しかも高精度に測定が可能になるからである。
【００１６】
まず、本発明に係るレオロジー指標演算処理装置１８でおこなわれる演算処理の理論的背景を以下で説明する。図３は、心拍に同期して拍動する動脈の模式図であって、動脈２８における圧力分布２９によって、血流は動脈の軸方向であるZ軸方向と半径方向に血流速度分布３０を生じる。この圧力分布２９が血圧値と相関する事はいうまでもない。さらに、動脈壁３１が弾性を持っているために、動脈壁３１はZ軸方向と半径方向に振動変位を起こす。この振動変位が図中の拍動変位３２である。さらに、この拍動変位３２は拍動にともなってZ軸方向に波動として、動脈壁を伝搬する。この波動が脈波３３である。  また、図３において、圧力分布２９、血流速度分布３０及び拍動変位３２は流体力学におけるナビエ・ストークスの方程式および動脈壁の力学的運動方程式を用いて解析的に求める事が可能である。すなわち、心拍の角振動数をω、脈波の波数をｋ、動脈２８の内径をＲ、動脈壁３１の厚み、ヤング率、密度及びポアソン比をそれぞれｈ、Ｅ、ρ０、σ、動脈２８の圧力分布２９をP、血流速度分布３０の軸方向速度成分をＶ、動脈半径方向の拍動変位をξとすると、Ｐ、Ｖ、ξはベッセル関数J０を用いて次式にて決定される。
【００１７】
【数１】
2006218169.tif 000003
【００１８】
【数２】
2006218169.tif 000004
【００１９】
【数３】
2006218169.tif 000005
ここで、φ、Fは次式で定義された無次元関数である。
【００２０】
【数４】
2006218169.tif 000006
【００２１】
【数５】
2006218169.tif 000007
また、γは、動脈壁厚みｈ、動脈壁密度ρ０及び動脈半径Rで以下のように定義された無次元パラメーターであり、次式で与えられている。
【００２２】
【数６】
2006218169.tif 000008
さらにαは、血液の運動粘性率ν、心拍周波数ω及び半径Ｒを用いて
【００２３】
【数７】
2006218169.tif 000009
と定義できる。
【００２４】
図１、図２記載の実施形態で検出される血流速度は、（２）式で与えられる血流速度分布Ｖにおける最大速度成分を振幅とした波形として検出される。この最大速度成分Ｖｍは、近似的に次式で与えられる。
【００２５】
【数８】
2006218169.tif 000010
となる。ここで、Ｋに関する詳細はのべないが、αの値に依存する値であって、０．６５〜1の値をとる。
以上の解析結果から、動脈内の心拍角振動数ωに同期して周期変動する軸方向の最大速度分布Vｍと拍動変位ξの位相差δは、（３）式と（８）式より、（９）式より決定される。
【００２６】
【数９】
2006218169.tif 000011
この位相差δは動脈内の圧力振幅Pｍとは無関係に無次元定数α、σ、γのみで決定される事が判明する。なぜならば、（９）式中に現れる、無次元関数F、Φはα、σ、γの関数であるからである。
【００２７】
また、同様に動脈内の心拍角振動数ωに同期して周期変動する軸方向の最大速度分布Vｍと拍動変位ξの振幅比μは、（３）式と（８）式より、（１０）及び（１１）式で決定される。
【００２８】
【数１０】
2006218169.tif 000012
【００２９】
【数１１】
2006218169.tif 000013
ここで、Γは規格化振幅比である。この振幅比μは、（９）式で与えられる位相差δと同様に動脈内の圧力振幅Pｍとは無関係に、動脈内径Ｒ、動脈壁厚ｈ、血液密度ρ、動脈ヤング率Ｅ及び無次元定数α、σ、γのみで決定される事が判明する。
【００３０】
図４は、（９）式で計算される最大血流速度Vmと拍動変位ξの位相差δのα依存性を理論計算した特性図であって、本発明に係る血液レオロジーの指標の演算処理を説明するための特性図である。縦軸はα、横軸は位相差δである。この特性図の特性曲線３４―Ａは、生体組織での代表的なσ、γの値
σ＝０．４〜０．６    γ＝０．０〜０．３
では大きな変化がなく、ほぼαの値の変化で決定される事が判明した。このαは（７）式で定義されているように、心拍振動数ω、血液の運動粘性率ν及び動脈内半径Rで決定される。よってこの位相差δ、心拍振動数数ω及び動脈内径Rを測定すれば、動脈内の圧力Pに無関係に、血液の運動粘性率νが検出できることになる。すなわち、血圧を測定せずに血液の運動粘性率νが検出できる事に他ならない。
【００３１】
図５は、（１１）式で計算される最大血流速度Vmと拍動変位ξの規格化振幅比Γと（７）式で定義されたαの関係を理論計算した特性図であって、本発明に係る血液レオロジーの指標の演算処理を説明するための特性図である。縦軸はα、横軸は規格化振幅比Γである。この特性図の特性曲線３４−Ｂは、生体組織での代表的な値σ、γの値
σ＝０．４〜０．６    γ＝０．０〜０．３
では大きな変化がなく、ほぼαの値の変化で決定される事が判明した。このαは（７）式で定義されているように、心拍振動数ω、血液の運動粘性率ν及び動脈内半径Rで決定される。よってこの規格化振幅比Γ、心拍振動数数ω動脈内径R、動脈壁厚ｈを測定すれば、通常、生体組織でのヤング率Ｅと血液密度ρは個体差もほとんどなく、一定値とみなせるので、動脈内の圧力Pに無関係に、血液の運動粘性率νが検出できることになる。すなわち、位相差と同様に血圧を測定せずに血液の運動粘性率νが検出できる事に他ならない。
【００３２】
ちなみに、従来の血液レオロジーの検出方法は、（２）式より計算される最大血流速度Vｍを圧力振幅項Pｍで除算するかわりに、血圧値で除算するものであった。以上が本発明に係る演算処理の理論的背景である。また、（１）式から（９）式の式中に現れるｊは虚数（−１の平方根）であり、計算値は複素数となるが、実際に意味のある物理量は、電子工学における交流理論と同様で、その実数部である事は言うまでもない事である。
【００３３】
図６、図７は、本発明に係る血液レオロジ測定装置の血管情報演算処理装置１７における演算処理を説明するための原理図である。図６は、図１記載の本発明に係る血液レオロジー測定装置の血管情報計測用の超音波センサ７を用いた、原理図である。
皮膚３５に装着された超音波センサ７の超音波の発信素子８より、発信バースト波３６が動脈３７に向けて放射される。この放射された超音波バースト波は、動脈３７の動脈外壁３８と動脈内壁３９に反射されるので、四種類の反射バースト波が受信素子９にて検出される。この四種類の反射バースト波がそれぞれ、第一反射バースト波４０、第二反射バースト波４１、第三反射バースト波４２、第四反射バースト波４３である。これら、四種類の反射バースト波において、第一反射バースト波４０と第四反射バースト波４３は、動動脈３７の動脈外壁３８にて反射する波であり、第二反射バースト波４１と第三反射バースト波４２が動脈内壁３９にて反射する波である。また、図６において、動脈３７と皮膚３５の距離は動脈離間距離４４、動脈３７と皮膚３５の角度は、動脈角度４５である。
【００３４】
これら四種類の反射バースト波は、それぞれ異なった時刻において受信素子９で検出される。図７は、受信素子９にて観測される四種類のバースト波の出力強度を図示した特性図である。図７において、縦軸は反射バースト信号波形４６と発信バースト信号波形５１の信号強度、横軸は時間である。図６記載の発信バースト波３６は発信素子８より生体中に対して周期的に放射される。この放射周期εは、心拍周期に比較して十分小さい。ｎ回目の発信バースト波３６の放射時刻が図７記載の発信バースト波形５１の出現時刻Tｎである。また、発信バースト波形５１の信号強度はVｎである。ここで、第一回目のバースト波放射時刻を０とすると、ｎ回目のバースト波放射時刻Ｔｎは、先の放射周期εを用いて、
【００３５】
【数１２】
2006218169.tif 000014
となる。
動脈３７から反射する反射バースト信号波形４６は、四個のピーク出力をもっている。すなわち、第一反射バースト信号４７、第二反射バースト信号４８、第三反射バースト信号４９及び第四反射バースト信号５０である。この時、それぞれ、第一反射バースト波４０は第一反射バースト信号４７に、第二反射バースト波４１は第二反射バースト信号４８に、第三反射バースト波４２は第三反射バースト信号４９に、第四反射バースト波４３は第三反射バースト信号５０に対応している。また、第一から第四反射バースト信号は、前述の発信バースト波形５１の出現時刻Tｎより遅れて出現し、それぞれ、Tｎ（１）、Tｎ（２）、Tｎ（３）、Tｎ（４）の時刻で観測される。ただし、この時間の遅れ量は、発信バースト波３６の放射周期より十分小さい。また、信号強度はそれぞれ、Vｎ（１）、Vｎ（２）、Vｎ（３）、Vｎ（４）である。
【００３６】
図６記載の動脈３７は、心臓の周期的運動にともない拍動運動を行っている。それゆえ、前記信号波形の出現時刻Tｎ（１）〜Tｎ（４）及び信号強度Vｎ（１）〜Vｎ（４）も時間に対して周期的に変化する。この時の時間周期が生体の心拍周期に他ならない。
本発明に係る図１及び図２記載の血管情報演算処理装置１７においては、図７で説明した反射バースト信号波形４６と発信バースト信号波形５１を任意の時間にて計測し、その波形情報として、それぞれの信号波形の出現時刻Tｎ（１）〜Tｎ（４）及びVｎ（１）〜Vｎ（４）が計測され、これらの計測値から、表１記載のパラメータすなわち、遅延時間差Δτと振幅強度比が検出される。
【００３７】
【表１】
2006218169.tif 000015
表１記載の数値は、発信バースト波５１の放射回数毎に計測される。この放射回数がサンプリング数である。この表１記載のパラメータを用いて、動脈３７の血管形状情報と拍動変位量が検出される。以下で、血管形状情報を検出する演算処理と拍動変位量を検出する演算処理を説明する。
本発明に係る血管情報演算処理装置１７にて検出される血管形状情報は、図６記載の動脈内径５２、動脈外径５３、動脈壁厚５４、動脈壁厚５５及び動脈内径と動脈壁厚の比（動脈壁厚／動脈内径）である。
前述の理論式（１）〜（１１）における動脈半径Ｒは、本発明に係る血管情報演算処理装置１７にて検出される動脈内径５２の時間平均値に等しい。心拍周波数に同期して動脈内径５２は周期的変化をする。動脈内径５２は図７記載の第二反射バースト信号４８と第三反射バースト信号４９の遅延時間差Δτ３（ｎ）に比例する。すなわち、発信バースト波３６のｎ回目放射時における動脈内径５２をＤ１（n）とすれば、
【００３８】
【数１３】
2006218169.tif 000016
となる。このＤ１（n）の時間平均値より、動脈内径Ｒは、
【００３９】
【数１４】
2006218169.tif 000017
となる。ここで、ｍはサンプリング数である。
同様に動脈外径５３は図７記載の第一反射バースト信号４７と第四反射バースト信号５０の遅延時間差に比例する。すなわち表１記載のΔτ２（ｎ）、Δτ３（ｎ）及びΔτ４（ｎ）の総和に比例する。よって、発信バースト波３６のｎ回目放射時における動脈外径５３をＤ２（n）とすれば、
【００４０】
【数１５】
2006218169.tif 000018
となる。このＤ２（n）の時間平均値より、動脈外径５３の平均値Ｄ２は、
【００４１】
【数１６】
2006218169.tif 000019
となる。
さらに動脈肉厚５４は、第一反射バースト信号４７と第二反射バースト信号４８の遅延時間差の時間平均に比例し、動脈肉厚５５は、第三反射バースト信号４９と第四反射バースト信号５０の遅延時間差の時間平均に比例する。すなわち、動脈肉厚５４の平均値をｈ１、動脈肉厚５５の平均値をｈ２とすれば、
【００４２】
【数１７】
2006218169.tif 000020
【００４３】
【数１８】
2006218169.tif 000021
として検出される。通常、ｈ１とｈ２は等しいので、動脈壁厚３１の値ｈとしてどちらの値を採用しても問題ない。もし大きく異なる場合は、ｈ１とｈ２の平均値を採用しても問題ない。すなわち、
【００４４】
【数１９】
2006218169.tif 000022
としてもよい。さらに、動脈内径Ｒと動脈壁厚ｈの比、ｈ／Ｒは、（１４）式、（１７）〜（１９）式を用いて、
【００４５】
【数２０】
2006218169.tif 000023
または、
【００４６】
【数２１】
2006218169.tif 000024
または、
【００４７】
【数２２】
2006218169.tif 000025
次に、拍動変位量を検出する演算処理を以下で説明する。拍動変位量は、動脈内径５２の時間変化量または動脈外径の時間変化量から検出される。すなわち、（１３）〜（１６）式を用いて、
【００４８】
【数２３】
2006218169.tif 000026
または
【００４９】
【数２４】
2006218169.tif 000027
として検出される。通常、ξ１とξ２は等しいので、拍動変位量の値ξとしてどちらの値を採用しても問題ない。もし大きく異なる場合は、ξ１とξ２の平均値を採用しても問題ない。すなわち
【００５０】
【数２５】
2006218169.tif 000028
としてもよい。
また他の方法として、図６記載の動脈離間距離４４が、表１記載の第一バースト信号４７の出現時刻Tｎ（１）と発信バースト信号波形５１の出現時刻Tｎの遅延時間差Δτ１（ｎ）の時間変化に比例する事を利用して拍動変位量を決定しても問題ない。
【００５１】
実際の生体中の反射バースト波は、生体中の各組織からの反射波が存在するので、図７にて示す動脈からの反射波以外に他の多くの反射波が存在する。それゆえ、数多くの反射波から、動脈からの反射波を検出する事が必要となる。まず、動脈からの反射波の大きな特徴は、心拍周波数に同期して周期的な変動を伴う事である。この周期的な変動は、前述の遅延時間差Δτにて観測されるだけでなく、表１記載の振幅強度比にても観測される。さらに、この振幅強度比の内、第一反射バースト波の振幅強度比は、動脈外形Ｄ２の大きさに比例する。
【００５２】
本発明に係る血液レオロジー測定装置にて計測される血流速度は、計測部位中の最大の外径（内径）を持つ動脈の中を流れる血流の速度である。それゆえ、心拍周波数に同期する遅延時間差及び振幅強度を持つ信号を検出し、それらの信号から、最大の振幅強度比を持つ反射バースト波を選択することで、動脈からの反射波が選択される。
【００５３】
次に、（１０）〜（１８）式にある動脈角度４５の大きさθは、後述するように図１記載の血流速度演算処理装置１６にて検出される。図８は、本発明に係るセンサ部の構造と生体内の動脈の位置関係を説明するための模式図であって、本発明に係る血液レオロジー測定装置に関する図１記載の実施例に対応する模式図である。２対の超音波センサ１と超音波センサ４及び超音波センサ７が同一のセンサ基板５６に配置されている。該２対の超音波センサは、連続超音波の射出および受信の指向性の方向が互いに平行にならないように、センサ角度５７をなすように傾けて配置されている。該センサ角度５７の大きさはψである。またバースト波を送受信する超音波センサ７は、超音波センサ１及び超音波センサ４の中間位置に配置されている。
【００５４】
また、超音波センサ１及び超音波センサ４は図１記載の超音波回路１２に接続されている。超音波センサ１、超音波センサ２は動脈６０に対する連続超音波の送信と該動脈６０中を流れる血流に対する反射連続超音波の受信を行い、その受信された超音波の周波数成分から、血流速度成分が血流速度演算処理装置１６にて検出される。また超音センサ７はバースト波の送受信を行い、バースト検波回路１４を経て血管情報演算処理装置１７にて拍動変位と血管形状情報の検出を行う。超音波センサ１、超音波センサ２及び超音波センサ７を構成している発信素子２と受信素子３、発信素子５と受信素子６及び発信素子８と受信素子９の材質は圧電セラミックスである。
【００５５】
センサ基板５６は音響整合層５８を介在して生体表面６１に装着される。生体組織５９の内部に存在する動脈６０とセンサ基板５６とのなす角度が動脈角度４５であり、その大きさはθである。動脈６０に対して超音波センサ１と超音波センサ４で発信された連続超音波は動脈６０の内部を流れる血流６２に反射し、ドップラー効果によるドップラーシフト（周波数シフト）をともなった反射波超音波として超音波センサ１及び超音波センサ４の受信素子に受信される。
【００５６】
受信された連続超音波は連続超音波検波回路１１と血流速度演算処理装置１６にてドップラーシフト量が検出される。さらに、この血流速度演算処理装置１６においては、動脈角度θと血流速度Ｖも決定される。すわわち、超音波センサ１にて観測されるドップラーシフト量をΔｆ１、超音波センサ４にて観測されるドップラーシフト量をΔｆ２、血流６２の血流速度をＶとすると、Δｆ１、Δｆ２、は、
【００５７】
【数２６】
2006218169.tif 000029
【００５８】
【数２７】
2006218169.tif 000030
となるので、この二式を血流速度Ｖとθの連立方程式として、Ｖとθが決定できる。その結果、（１０）式〜（１８）式中のｓｉｎθは、
【００５９】
【数２８】
2006218169.tif 000031
と決定できる。また、この（２６）式よりθの値を求め、（２４）式または（２５）式に代入する事で、血流速度Ｖを検出する事ができる。
【００６０】
図９は、本発明に係るセンサ部の構造と生体内の動脈の位置関係を説明するための第二の模式図であって、本発明に係る血液レオロジー測定装置に関する図２記載の実施例に対応する模式図である。２対の超音波センサ１と超音波センサ２４が同一のセンサ基板６３に配置されている。該２対の超音波センサは、超音波の射出および受信の指向性の方向が互いに平行にならないように、図８と同様にセンサ角度５７をなすように傾けて配置されている。該センサ角度５７の大きさはψである。さらに、超音波センサ１と超音波センサ２４を構成している発信素子２と受信素子３、発信素子２５と受信素子２６の材質は圧電セラミックスである。
【００６１】
図２の説明で述べたように、図９記載のセンサ構造にては、超音波センサ２４を構成している発信素子２５と受信素子２６は図１記載の超音波回路１２と超音波バースト回路１５の双方に切り替え回路２７を介して接続されている。この切り替え回路２７は、超音波センサ２４の接続回路を周期的に超音波回路１２と超音波バースト回路１３に切り替える。すなわち、超音波センサ２４は、血流速度検出と血管情報検出の双方の機能を持つ事になる。この場合、求めるべき動脈角度６４の大きさθは、超音波センサ２４にて観測されるドップラーシフト量Δｆ２、超音波センサ１にて観測されるドップラーシフト量Δｆ１を用いて、
【００６２】
【数２９】
2006218169.tif 000032
から決定できる。さらに、血流速度Ｖは、
【００６３】
【数３０】
2006218169.tif 000033
と決定される。
【００６４】
図１０は、血流速度演算処理装置１６において演算処理された心拍信号に同期して周期的変化する血流速度信号の波形（血流速度波形）６５と血管情報演算処理装置１７において演算処理された同じく心拍信号に同期して周期的変化する拍動変位波形６６を示す説明図である。縦軸は出力強度であって、血流速度波形６５においては、速度強度を表し、拍動変位波形６６においては、（２１）または（２２）または（２３）式にて、検出された拍動変位量ξを表している。横軸は時間である。図中の血流速度波形６５及び拍動変位波形６６は双方ともに時間間隔ΔTの範囲でN個のピーク値を持っている。
【００６５】
すなわち、血流速度波形６５においては、ＶＰ(1)からVＰ(N)のN個のピーク値、拍動変位波形６６においては、ξＰ(1)からξＰ(Ｎ)のN個のピーク値である。  これらのピーク値は、血流速度波形６５においては、そのベースライン６７からの計測値であり、拍動変位波形６６においては、そのベースライン６８からの計測値である。それゆえ、これらのピーク値は、両波形の振幅強度である。
表２に図１０記載の血流速度波形６５のピーク値VＰ(1)とその出現時刻τＶ(ｎ)及び拍動変位波形６６のピーク値ξＰ(n)とその出現時刻τｈ(n)をまとめて記す。
【００６６】
【表２】
2006218169.tif 000034
表２及び血流速度波形６５のN個のピーク値は（２）式で与えられた血流速度分布の最大血流速度Vｍ対応している。
【００６７】
次に、本発明に係るレオロジー演算処理装置１８に内蔵されている波形情報演算処理装置１９、心拍周波数演算処理装置２０、パラメータ演算処理装置２１及びレオロジー指標演算処理装置２２について以下説明する。第一の演算処理装置である波形情報演算処理装置１９は、表２に記載の波形情報（ピーク値）を検出するものであり、この波形情報演算処理装置１９における波形に関する演算処理方法は、コンパレーターを用いたピーク検出方法等を用いた演算処理方法を採用している。
次に、第二の演算処理装置である心拍周波数演算処理装置２０は、表２に記載の波形情報のうち、血流速度波形６５または拍動変位波形６６のピーク値出現の時刻間隔を以下の演算処理式（１）によって
【００６８】
【数３１】
2006218169.tif 000035
または
【００６９】
【数３２】
2006218169.tif 000036
をもとめ、さらに、以下の第二の演算処理式（２）
【００７０】
【数３３】
2006218169.tif 000037
によって心拍周波数Fを検出する。この演算処理方法はピーク値を基準にした演算処理方法であるが、図１０に記載の血流速度波形６５と拍動変位波形６６の最小値を基準にしてもなんら本質的に問題ない。
【００７１】
本発明に係るレオロジー演算処理装置１８に内蔵される第三の演算処理装置であるパラメータ演算処理装置２１においては、血流速度波形６５と拍動変位波形６６の両波形から、演算処理によって、位相差または振幅比の少なくともどちらか一方が検出される。
【００７２】
まず、位相差の検出について説明する。位相差の検出においては、表２記載のＮ個の血流速度波形６５と拍動変位波形６６のピーク出現時刻の時間差を、以下の演算処理式（３）、
【００７３】
【数３４】
2006218169.tif 000038
によって求めると共に、この時間差ΔτＶｈ（n）の平均値と心拍周波数演算処理装置１９にて検出された心拍周波数Fを用いて、以下の演算処理式（４）
【００７４】
【数３５】
2006218169.tif 000039
によって血流速度波形６５と拍動変位波形６６の位相差δを検出している。
【００７５】
また、振幅比μにおいては、表２記載のＮ個の血流速度波形６５と拍動変位波形６６のピーク値の比ＶＰ（ｎ）／ξＰ（ｎ）の平均値として検出される。すなわち、
【００７６】
【数３６】
2006218169.tif 000040
として検出される。
以上の演算処理のうち、位相差に関する演算処理はピーク値を基準にした演算処理方法であるが、心拍周波数演算処理装置２０と同様に図１０に記載の血流速度波形６５と拍動変位波形６６の最小値を基準にしてもなんら本質的に問題ない。
【００７７】
さらに本発明に係る波形情報演算処理装置１９において、波形情報の演算処理手段としてフーリエ解析法、フェーズロックトループ（PLL）法などの一般的演算処理方法を採用することで、血流速度波形６５及び拍動変波形６６から直接、心拍周波数Fと位相差δを検出する構成であってもなんら問題はなく、適宜、変更可能である。
【００７８】
本発明に係るレオロジー演算処理装置１８に内蔵される第四の演算処理装置であるレオロジー指標演算処理装置２２においては、血管情報演算処理装置１７にて検出された動脈半径Ｒ等の血管情報、心拍周波数演算処理装置２０にて検出された心拍周波数Ｆ及びパラメータ演算処理装置２１にて検出された、血流速度波形６５と拍動変波形６６の位相差δまたは振幅比μより、運動粘性率νが検出される。
【００７９】
以下で、血流速度波形６５と拍動変波形６６の位相差δより血液の運動粘性率νを決定する場合について説明する。この場合は、レオロジー指標演算処理装置２２に記憶されている図４記載の特性曲線３４−Ａ、パラメータ演算処理装置２１（演算処理式（４））にて検出された位相差δ、心拍周波数演算処理装置２０（演算処理式（２））にて検出された心拍周波数Ｆ、及び血管情報演算処理装置１７（（１４）式）にて検出された動脈半径Ｒより、血液運動粘性率νが決定される。まず、記憶されている特性曲線３４−Ａとパラメータ演算処理装置２１にて検出された位相差δから、αの値が検出される。この検出されたαの値をα１とすると、心拍周波数Ｆ及び動脈半径Ｒより、以下の演算処理式（６）より、求めるべき血液の運動粘性率νが決定される。
【００８０】
【数３７】
2006218169.tif 000041
この演算処理式（６）を用いて、検出された運動粘性率νとマイクロチャネルアレイを用いた採血方式による血液レオロジーの指標である全血通過時間Tの相関を示した特性図が、図１１である。先に説明したように、本発明に係る血液レオロジー測定装置にて検出した血液の運動粘性率νと全血通過時間Ｔは、高度の相関を持っている事が判明している。
【００８１】
次に、血流速度波形６５と拍動変波形６６の振幅比μより血液の運動粘性率νを決定する場合について説明する。この場合は、図５記載の特性曲線３４−Ｂ、血液密度ρ及び血管ヤング率Ｅがレオロジー指標演算処理装置２２に記憶されている。さらにパラメータ演算処理装置２１（演算処理式（５））にて検出された振幅比μ、心拍周波数演算処理装置２０（演算処理式（２））にて検出された心拍周波数Ｆ、及び血管情報演算処理装置１７（（１４）式、（２０）式）にて検出された動脈半径Ｒ及びｈ／Ｒより、血液運動粘性率νが決定される。
まず、検出された振幅比μ、動脈半径Ｒ及びｈ／Ｒ、さらに記録されている血液密度ρ及び血管ヤング率Ｅより、（１０）、（１１）式で定義された規格化振幅比Γが検出される。
この検出された規格化振幅比をΓ０とすると、Γ０は
【００８２】
【数３８】
2006218169.tif 000042
となる。
さらにこのΓ０と記録されている特性曲線３４−Ｂをもちいて、αの値が検出される。この検出されたαの値をα２とすると、心拍周波数Ｆ及び動脈半径Ｒより、以下の演算処理式（８）を用いて、求めるべき血液の運動粘性率νが決定される。
【００８３】
【数３９】
2006218169.tif 000043
ちなみに、振幅比より運動粘性率νを算出するために、記憶されている血液密度ρ及び血管ヤング率Ｅの値は、それぞれ、生体の個体差はほとんどないので、図１１記載の特性図と同等な検出精度が得られることが判明している。
【産業上の利用可能性】
【００８４】
本発明は、医療および健康維持・増進を目的として、体液の流動性を示す指標としての血液レオロジーと強い相関がある生体内の血流速度を計測することが可能であるだけでなく、生体（人体）の活動状況と生体各部における血流状態の相関を知るための計測においても利用可能である。
【図面の簡単な説明】
【００８５】
【図１】本発明に係る血液レオロジー測定装置の構成を示すブロック図
【図２】本発明に係る血液レオロジー測定装置の構成を示す第二のブロック図
【図３】本発明に係る心拍に同期して拍動する動脈の模式図
【図４】本発明に係る血液レオロジーの指標算出のための演算処理を説明する特性図
【図５】本発明に係る血液レオロジーの指標算出のための演算処理を説明する特性図
【図６】本発明に係る血液レオロジーの指標算出のための演算処理を説明する模式図
【図７】本発明に係る血液レオロジーの指標算出のための演算処理を説明する模式図
【図８】本発明に係るセンサ部の構造と血液レオロジー測定を説明するための模式図
【図９】本発明に係るセンサ部の構造と血液レオロジー測定を説明するための第二の模式図
【図１０】本発明に係る血流速度波形と拍動変位波形を示す説明図
【図１１】本発明の効果を示す特性図
【符号の説明】
【００８６】
１  超音波センサ
２  発信素子
３  受信素子
４  超音波センサ
５  発信素子
６  受信素子
７  超音波センサ
８  発信素子
９  受信素子
１０  連続超音波発信回路
１１  連続超音波検波回路
１２  超音波回路
１３  バースト発生回路
１４  バースト検波回路
１５  超音波バースト回路
１６  血流速度演算処理装置
１７  血管情報演算処理装置
１８  レオロジー演算処理装置
１９  波形情報演算処理装置
２０  心拍数演算処理装置
２１  パラメータ演算処理装置
２２  レオロジー指標演算処理装置
２３  出力装置
【図１】
2006218169.tif 000044
【図２】
2006218169.tif 000045
【図３】
2006218169.tif 000046
【図４】
2006218169.tif 000047
【図５】
2006218169.tif 000048
【図６】
2006218169.tif 000049
【図７】
2006218169.tif 000050
【図８】
2006218169.tif 000051
【図９】
2006218169.tif 000052
【図１０】
2006218169.tif 000053
【図１１】
2006218169.tif 000054