JP2010042284A - 血液力学特性測定装置 - Google Patents
血液力学特性測定装置 Download PDFInfo
- Publication number
- JP2010042284A JP2010042284A JP2009244039A JP2009244039A JP2010042284A JP 2010042284 A JP2010042284 A JP 2010042284A JP 2009244039 A JP2009244039 A JP 2009244039A JP 2009244039 A JP2009244039 A JP 2009244039A JP 2010042284 A JP2010042284 A JP 2010042284A
- Authority
- JP
- Japan
- Prior art keywords
- distribution
- blood
- viscosity
- shear stress
- dimensional
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Landscapes
- Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
【解決手段】経皮的に血管内の血流に向けて超音波等の波動を入射させ、その血管内の血液中に含まれる血球成分からの散乱エコーを受信した後、超音波ドプラ計測によって血流速度分布を計測する手段4と、血流速度分布に基づいて、2次元または3次元のナビエ-ストークス方程式より、血液の粘性率分布及び/または血液の粘性率分布平均値を求める手段5,6と、計測された血流速度分布からずり速度分布を算出するずり速度分布算手段7と、算出されたずり速度分布と血液の粘性率分布または血液の粘性率分布平均値とに基づいて血液内のずり応力分布を求める手段9とからなることを特徴とする血液力学特性測定装置である。
【選択図】図1
Description
このように、式(4)を基本式とし、速度成分分布(速度プロファイル)から、粘性率と圧勾配を求める方法については、以下のようなものが提案されている。
この方法は、速度プロファイルを時間軸方向に周波数分解し、その結果得られた拍動成分の速度プロファイルを基に、管内の2点の速度比をとって粘性率を推定するものである。しかし、この方法は時間軸方向の周波数解析を必要としているので、超音波診断技術に求められる実時間性を満たさない。また、速度比をとる2点の選び方によって、結果が変わってしまう問題もある。
この方法は、式(4)が、半径r方向に2回積分できるので、粘性率及び密度を空間的に一様な定数として、2つの異なる位置で積分すると2式からなる連立方程式が得られことから、これを解くことにより粘性率を求めるものである。しかし、この方法は2つの積分路の選び方が非常に微妙であり、積分路の選び方によって結果が変わってくる。この点で、線形化法と同じ問題がある。
この方法は、この非特許文献2に提案されている方法であり、上記(1)(2)のように、2点の計測だけでなく、多数点をとって、最尤推定法(最小二乗法)により粘性率と圧勾配を同時に推定する方法である。空間的に平均化して算出するので推定値の安定化が見込まれる。
これらの方法については、シミュレーションと流体実験により検証され、その結果については手法の妥当性が示されている。
第1の手段は、経皮的に血管内の血流に向けて超音波等の波動を入射させ、その血管内の血液中に含まれる血球成分からの散乱エコーを受信した後、超音波ドプラ計測によって血流速度分布を計測する手段と、前記計測された血流速度分布に基づいて、2次元または3次元のナビエ-ストークス方程式を利用して、前記血液の粘性率分布及び/または前記血液の粘性率分布平均値を求める手段とからなることを特徴とする血液力学特性測定装置である。
また、粘性率分布及びずり応力分布の画像により、血球成分(高粘性)と血漿成分(低粘性)の含有割合など、血液の組成を把握することが可能となる。これにより、例えば、動脈硬化プラークに接する血液の粘性率、あるいは血流のずり応力が通常よりも高い、などの測定結果から、プラーク破綻の危険性を正確かつ容易に予知することが可能となり、血管系疾患の防止に寄与ですることがきる。
同図において、1は生体断層面描出及び血管内の血流速度分布計測に用いる超音波プローブであり、プローブ表面上に、1次元配列または2次元配列された素子から超音波を発生させて生体内に入射し、測定対象血管2内の血球からの散乱波及び血管周辺組織3からの散乱波を、同一の素子を経由して受信するものであり、受信した散乱波は、断層面の構成、及び血流速度分布計測に用いられる。
まず、上記で説明した生体の体表より非侵襲で血液の粘性率分布、ずり速度分布、ずり応力分布、及び粘性率分布平均値、ずり速度分布平均値、ずり応力分布平均値を求めるために、生体断層面描出及び血管内の血流速度分布計測に用いる超音波プローブ1を体表面に配置し接触させる。プローブ表面上に、1次元配列または2次元配列された素子から超音波を発生させて生体内に入射し、測定対象血管2内の血球からの散乱波及び血管周辺組織3からの散乱波を、同一の素子を経由して受信する。受信した散乱波は、断層面の構成及び血流速度分布計測に用いられる。
同図において、12は速度分布サブ領域分割部であり、速度分布測定部4で得られた、2次元速度ベクトル分布、または3次元速度ベクトル分布を、細かなサブ領域13に分割する機能を有する。本発明において2次元または3次元のナビエ-ストークス方程式に基づき粘性率を算出するのに、血液が非圧縮性であると仮定しているため、粘性率は空間的に一意であるとして算出される。しかしながら、血液は赤血球や血小板などの混在流体であり、粘性率は一意とならず、分布を持つ。この分布を得るために、サブ領域13に分割し、該サブ領域13内では粘性率が一様であるとして、粘性率を算出する。これを分割した全てのサブ領域13に対して反復することにより、粘性率分布を得るものである。
同図において、16は速度分布測定部4で得られた、2次元速度ベクトル分布または3次元速度ベクトル分布の断面図、17は粘性率算出部14で得られた全てのサブ領域13における粘性率に対して、粘性率分布構成部(サブ領域統合)15を経て求められた粘性率分布の断面図、18は粘性率分布平均値算出部6で算出された粘性率分布平均値である。
同図に示すように、ずり速度分布算出部9において、粘性率分布算出部5において算出された血液の粘性率分布または粘性率分布平均値算出部6において算出された粘性率分布平均値と、ずり速度分布算出部7において算出されたずり速度分布とを入力することにより、ニュートン粘性法則式に従って、(ずり応力分布)=(粘性率分布)×(ずり速度分布)、または(ずり応力分布)=(粘性率分布平均値)×(ずり速度分布)の演算を行い、ずり応力分布を算出する。
このようにして得られたずり応力分布は、ずり応力分布平均値算出部10にて、ずり応力分布平均値を求めることも可能である。また、粘性率分布時系列の算出と同様な考え方により、ずり応力分布の時系列を求めることも可能である。
同図において、19はずり速度分布算出部7で得られたずり速度分布の断面図、20はずり応力分布の断面図、21はずり応力分布平均値算出部10で得られたずり応力分布平均値、22はニュートン粘性法則式に従い、サブ領域13ごとに得られた粘性率とずり速度の乗算により、サブ領域13ごとにずり応力を算出するずり応力算出部である。
典型的な血液動態では血管内皮側で大きなずり応力を生ずるが、ずり応力分布算出部9において、ずり応力分布を適切に算出することにより、定量的なずり応力分布の評価が可能となる。
同図において、23はずり速度分布算出に必要なひずみ速度分布を得るために、速度成分分布を微分するための速度成分分布微分演算部であり、速度分布測定部4で得られた2次元速度ベクトル分布、または3次元速度ベクトル分布に対して、空間微分演算を行う。
式15a〜式15fのように求められたひずみ速度テンソルから、ずり速度分布抽出部25にてせん断成分を抽出し、3次元の場合のずり速度分布については式15d〜式15f、2次元の場合のずり速度分布については式15dが求められる。すなわち、3次元の場合のずり速度成分として求められるのは、exy、eyz、ezxであり、2次元の場合のずり速度成分として求められるのは、exyである。また、このようにして求められたずり速度分布に対し、ずり速度分布平均値算出部8を適用することにより、ずり速度分布平均値を算出することも可能である。また、粘性率分布時系列の算出と同様な考え方により、ずり速度分布の時系列を求めることも可能である。
同図において、26はサブ領域13ごとにずり速度を算出するずり速度算出部、27はずり速度分布平均値算出部8で得られたずり速度分布平均値である。
速度分布測定部4において求められた2次元速度ベクトル分布、または3次元速度ベクトル分布16に対して、サブ領域13に分割し、各々のサブ領域13に対して、ずり速度算出部26にて、ずり速度を算出する。ここで求められるずり速度は、3次元速度ベクトル分布を対象とした場合は式15d〜式15f、2次元速度ベクトル分布を対象とした場合は式15dに基づいて算出される。このように全てのサブ領域で算出されたずり速度値を、元のサブ領域に再配置し統合することによって、ずり速度分布19を得ることができる。さらに、全てのサブ領域で得られたずり速度を、ずり速度分布平均値算出部8にて平均化することで、空間的に平均化されたずり速度分布平均値27を算出することも可能である。
28は速度分布測定部4にて得られた3次元速度ベクトル分布の時系列、または2次元速度ベクトル分布の時系列の内、1成分のみを表示した血流速度波形である。
2 計測対象血管
3 血管周辺組織
4 速度分布測定部
5 粘性率分布算出部
6 粘性率分布平均値算出部
7 ずり速度分布算出部
8 ずり速度分布平均値算出部
9 ずり応力分布算出部
10 ずり応力分布平均値算出部
11 表示部
12 速度分布サブ領域分割部
13 サブ領域
14 粘性率算出部
15 粘性率分布構成部(サブ領域統合)
16 2次元速度ベクトル分布または3次元速度ベクトル分布の断面図
17 粘性率分布の断面図
18 粘性率分布平均値
19 ずり速度分布の断面図
20 ずり応力分布の断面図
21 ずり応力分布平均値
22 ずり応力算出部
23 速度成分分布微分演算部
24 ひずみ速度テンソル分布算出部
25 ずり速度分布抽出部
26 ずり速度算出部
27 ずり速度分布平均値
28 血流速度波形
29 ずり速度分布、粘性率分布、ずり応力分布の時系列
30 心電同期によるずり速度・粘性率・ずり応力の抽出部
31 抽出されて時系列に並べられたずり速度・粘性率・ずり応力
Claims (2)
- 経皮的に血管内の血流に向けて超音波等の波動を入射させ、その血管内の血液中に含まれる血球成分からの散乱エコーを受信した後、超音波ドプラ計測によって血流速度分布を計測する手段と、前記計測された血流速度分布からずり速度分布を算出するずり速度分布算出手段と、血液の粘性率分布が一様であるとの条件下において、前記算出されたずり速度分布から血液内のずり応力分布を求める手段とからなることを特徴とする血液力学特性測定装置。
- 前記血液の粘性率分布及び前記血液内のずり応力分布を、実時間で計測して画像表示すると同時に、前記血液の粘性率分布及び血液内のずり応力分布の平均値を数値として表示する手段とを備えたことを特徴とする請求項1記載の血液力学特性測定装置。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2009244039A JP4863325B2 (ja) | 2009-10-23 | 2009-10-23 | 血液力学特性測定装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2009244039A JP4863325B2 (ja) | 2009-10-23 | 2009-10-23 | 血液力学特性測定装置 |
Related Parent Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2004359792A Division JP4474530B2 (ja) | 2004-12-13 | 2004-12-13 | 血液力学特性測定装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2010042284A true JP2010042284A (ja) | 2010-02-25 |
JP4863325B2 JP4863325B2 (ja) | 2012-01-25 |
Family
ID=42014118
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2009244039A Expired - Fee Related JP4863325B2 (ja) | 2009-10-23 | 2009-10-23 | 血液力学特性測定装置 |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP4863325B2 (ja) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2011235009A (ja) * | 2010-05-13 | 2011-11-24 | Hitachi Aloka Medical Ltd | 超音波診断装置及び該装置における数値シミュレーション方法 |
CN103917165A (zh) * | 2011-08-26 | 2014-07-09 | Ebm株式会社 | 用于血流性状诊断的系统、其方法及计算机软件程序 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH11316180A (ja) * | 1998-01-23 | 1999-11-16 | Koninkl Philips Electronics Nv | 血管中の粘度及び圧勾配を決定するエコ―検査方法及び装置 |
JP2003144395A (ja) * | 2001-11-09 | 2003-05-20 | Fukuda Denshi Co Ltd | 血管内皮機能測定装置 |
JP2003290226A (ja) * | 2002-04-03 | 2003-10-14 | Seiko Instruments Inc | 循環動態測定装置、循環動態測定方法、血圧測定方法および循環動態センサ |
-
2009
- 2009-10-23 JP JP2009244039A patent/JP4863325B2/ja not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH11316180A (ja) * | 1998-01-23 | 1999-11-16 | Koninkl Philips Electronics Nv | 血管中の粘度及び圧勾配を決定するエコ―検査方法及び装置 |
JP2003144395A (ja) * | 2001-11-09 | 2003-05-20 | Fukuda Denshi Co Ltd | 血管内皮機能測定装置 |
JP2003290226A (ja) * | 2002-04-03 | 2003-10-14 | Seiko Instruments Inc | 循環動態測定装置、循環動態測定方法、血圧測定方法および循環動態センサ |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2011235009A (ja) * | 2010-05-13 | 2011-11-24 | Hitachi Aloka Medical Ltd | 超音波診断装置及び該装置における数値シミュレーション方法 |
CN103917165A (zh) * | 2011-08-26 | 2014-07-09 | Ebm株式会社 | 用于血流性状诊断的系统、其方法及计算机软件程序 |
Also Published As
Publication number | Publication date |
---|---|
JP4863325B2 (ja) | 2012-01-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106028948B (zh) | 超声波摄像装置以及方法 | |
Kanai et al. | Elasticity imaging of atheroma with transcutaneous ultrasound: preliminary study | |
Ferraioli et al. | Ultrasound point shear wave elastography assessment of liver and spleen stiffness: effect of training on repeatability of measurements | |
JP4474530B2 (ja) | 血液力学特性測定装置 | |
KR101398948B1 (ko) | 진폭―위상 변조된 초음파를 이용한 점탄성 측정 | |
Beulen et al. | Toward noninvasive blood pressure assessment in arteries by using ultrasound | |
Balocco et al. | Estimation of the viscoelastic properties of vessel walls using a computational model and Doppler ultrasound | |
JP2006218169A (ja) | 血液レオロジー測定装置および血液レオロジー計測方法 | |
JP5474986B2 (ja) | 血管機能検査装置 | |
JP5224221B2 (ja) | 血管機能検査装置 | |
JP4863324B2 (ja) | 血液力学特性測定装置 | |
JP4863325B2 (ja) | 血液力学特性測定装置 | |
Thomas | Tissue Doppler echocardiography–A case of right tool, wrong use | |
JP6996035B2 (ja) | 超音波診断装置、および、生体組織の物性評価方法 | |
Altamirano-Diaz et al. | Does obesity affect the non-invasive measurement of cardiac output performed by electrical cardiometry in children and adolescents? | |
JP6523918B2 (ja) | 超音波撮像装置 | |
Olesen et al. | Noninvasive estimation of 2-D pressure gradients in steady flow using ultrasound | |
JP6783007B2 (ja) | 組織弾性の測定装置及び測定方法 | |
JP4880971B2 (ja) | 血液循環状態測定装置 | |
Thorne et al. | In vivo Doppler ultrasound quantification of turbulence intensity using a high-pass frequency filter method | |
Padilla et al. | Scattering by trabecular bone | |
Iijima | Approaches to the diagnosis of liver fibrosis | |
CN106963424A (zh) | 检测动脉血管粘弹性的系统及方法 | |
JP5958806B2 (ja) | 超音波診断装置および血流量推定プログラム | |
JP4475457B2 (ja) | 膠原線維割合測定装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
TRDD | Decision of grant or rejection written | ||
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20111101 |
|
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20111102 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20141118 Year of fee payment: 3 |
|
R150 | Certificate of patent or registration of utility model |
Ref document number: 4863325 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20141118 Year of fee payment: 3 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
S533 | Written request for registration of change of name |
Free format text: JAPANESE INTERMEDIATE CODE: R313533 |
|
R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
LAPS | Cancellation because of no payment of annual fees |