(57)【要約】  非侵襲・高精度・高い時間分解能で心内圧を計測する。  超音波造影剤から反射された信号の減衰曲線、共鳴曲線から圧力を計算する。また、多成分液滴型造影剤を体内に注入し、検出したい圧力範囲を気液平衡になるように組成を設定し、気泡径を圧力に依存させて、減衰曲線、共鳴曲線へ圧力を強く反映させ、高精度な圧力計測を実現する。
2010147055.tif 000002
【特許請求の範囲】
【請求項１】
超音波造影剤を注入した撮像対象に対して超音波を送受信する超音波探触子と、
前記超音波探触子が受信した撮像対象からの超音波信号を信号処理する信号処理部と、
前記信号処理部の処理結果を表示する表示手段とを備え、
前記信号処理部は、撮像対象内の任意の空間位置に焦点を設定し、前記焦点に位置する造影剤が反射した信号から、圧力を検出することを特徴とする超音波撮像装置。
【請求項２】
超音波造影剤を注入した撮像対象に対して超音波を送受信する超音波探触子と、
前記超音波探触子が受信した撮像対象からの超音波信号を信号処理する信号処理部と、
前記信号処理部の処理結果を表示する表示手段とを備え、
前記信号処理部は、前記受信した超音波信号の周波数分布を算出する周波数解析部と、
前記周波数解析部によって算出された周波数分布から、前記受信した超音波信号の減衰率の周波数依存性の、１つ以上の規定の周波数での値、最大値を示す周波数、最大値、半値幅のうち１つ以上である、周波数分布特徴量を検出する周波数分布特徴量検出部と、
前記周波数分布特徴量を圧力に換算する周波数分布特徴量−圧力換算部とを有し、
前記表示手段は、前記周波数分布特徴量−圧力換算部で算出された圧力の値を表示することを特徴とする超音波撮像装置。
【請求項３】
請求項２に記載の超音波撮像装置において、前記周波数分布特徴量−圧力換算部は、前記周波数解析部の算出した周波数分布から、基本波、高調波、分調波、差周波、和周波の強度の周波数分布、値の分岐の大きさのうち１つ以上である周波数分布特徴量を検出することを特徴とする超音波撮像装置。
【請求項４】
請求項２に記載の超音波撮像装置において、前記信号処理部は、体内の圧力を検出した部位の動き量を検出する動き量検出部と、前記周波数分布特徴量−圧力換算部が換算した圧力と前記動き量検出部が検出した動き量から、体内の圧力と動きを検出した部位の硬さである硬さ量を算出する圧力・動き量−硬さ量換算部を備えることを特徴とする超音波撮像装置。
【請求項５】
撮像対象に超音波造影剤を注入する造影剤注入ステップと、
撮像対象に超音波信号を送波する超音波送波ステップと、
撮像対象から反射された超音波信号を受信し、受信した超音波信号の周波数分布を算出する周波数解析ステップと、
前記算出した周波数分布から、前記受信した超音波信号の減衰率の周波数依存性の、１つ以上の規定の周波数での値、最大値を示す周波数、最大値、半値幅のうち１つ以上である、周波数分布特徴量を検出する周波数分布特徴量ステップと、
前記検出した周波数分布特徴量を圧力に換算する周波数分布特徴量−圧力換算ステップと
を有することを特徴とする超音波撮像方法。
【請求項６】
請求項５に記載の超音波撮像方法において、
前記造影剤注入ステップは、複数の液体成分を内包する多成分液滴型造影剤を体内に注入し、
前記超音波送波ステップは、前記多成分液滴型造影剤を気化させるための気化用超音波送波を行った後、圧力検出用の超音波を送波することを特徴とする超音波撮像方法。
【請求項７】
請求項６に記載の超音波撮像方法において、圧力を検出する部位の温度をＴc、検出しようとする圧力の上限をＰu、下限をＰlとするとき、前記複数の液体成分のうち少なくとも１つは温度Ｔcでの蒸気圧がＰu以上であり、他の少なくとも１つは温度Ｔcでの蒸気圧がＰl以下であり、前記複数の液体成分はＰl以上、Ｐu以下の圧力で気液平衡が成り立つモル分率で混合されていることを特徴とする超音波撮像方法。
【請求項８】
請求項５に記載の超音波撮像方法において、さらに、
体内の圧力を検出した部位の動き量を検出するステップと、
前記換算した圧力と前検出した動き量から前記圧力と動きを検出した部位の硬さである硬さ量を算出するステップを有することを特徴とする超音波撮像方法。
【請求項９】
請求項６に記載の超音波撮像方法において、さらに、
体内の圧力を検出した部位の動き量を検出するステップと、
前記換算した圧力と前検出した動き量から前記圧力と動きを検出した部位の硬さである硬さ量を算出するステップを有することを特徴とする超音波撮像方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、医療用の超音波撮像装置及び方法に関し、特に非侵襲・高精度・高い時間分解能で心内圧を計測する超音波撮像装置及び方法に関する。
【背景技術】
【０００２】
心臓疾患は老化にともなう機能不全という意味ではすべての人が患う疾患であり、多くの先進国で３大死因の１つとなっている。このように経過の長い、患者の多い疾患の医療には、予防医療から高度医療まで広いフェーズがあるが、すべてのフェーズに共通する技術的なニーズとして、患者への負担の小ささがある。患者への負担の小さい検査技術があれば、早期診断をより初期へ、経過観察をよりきめ細かくでき、医療の質を向上することができる。【０００３】
血圧計測は心臓疾患の検査の基本である。現在主に行われている血圧計測方法は、予防医療では、上腕などをカフで段階的に圧迫して拍動音が消失する圧迫圧を計測値する圧迫検査、高度医療では、血管から心臓へ管を挿入して体外に残した管の端の圧力を計測値とするカテーテル検査である。【０００４】
また、研究段階であるが、超音波撮像装置によって血圧を非侵襲に計測する方法が報告されている（特許文献１）。この方法は、液体中で超音波を照射された造影剤は、液体の圧力を線形に反映した信号強度の分調波を有するという現象を発見し、利用した方法である。手順としては、超音波造影剤を静脈注射したあと超音波診断装置で心臓を撮像し、分調波の信号強度（ｄＢ）を検出して、あらかじめ調べておいた分調波の信号強度と圧力の関係を用いて受信信号を血圧に換算している。
【先行技術文献】
【特許文献】
【０００５】
【特許文献１】米国特許第６，３０２，８４５号明細書
【発明の概要】
【発明が解決しようとする課題】
【０００６】
しかし、上記従来技術では、非侵襲性と、計測精度、時間分解能が両立できていない。具体的には、圧迫検査は計測値が心内圧ではないため、心内圧とみなすには精度が不足し、確定的な心臓検査として用いることはできない。カテーテル検査は侵襲的で、患者への負担が大きい。【０００７】
特許文献１の方法は、非侵襲に心内圧を計測する方法であるが、精度と時間分解能に問題がある。詳しくは、まず、時間分解能と血圧精度が両立しない点である。この方法では分調波から圧力を算出しているが、一般的に分調波は高いＳＮ比では得られない。そこで、複数心拍分の計測を行って時相の同じ受信信号の平均値をとる、あるいは、時間区間を設定して時間区間内の受信信号の平均値をとることで、ＳＮ比を高めている。すなわち、精度を向上させるかわりに時間分解能を犠牲にしている。しかし、心内圧を医療行為として測定する場合には０．１秒以上の時間分解能が必要であり、かつ、複数心拍の時相平均処理は不整脈患者を検査対象から外すという問題があり、好ましくない。診断に使う値は、左室拡張末期圧など、部位と拍動周期のタイミングを指定した血圧であるからである。また、不整脈患者を検査対象外にしないためにも時間分解能は必要である。非健常者は医療上、重要かつ主な検査対象で、彼らは必然的に不整脈のあることが多い。次に、低い圧力範囲で精度が得られない問題がある。低い圧力範囲（２０〜６０ｍｍＨｇ）は心疾患の診断に重要な役割を果たす左室拡張末期圧の値域で、±５ｍｍＨｇ以上の精度が必要である。このように、特許文献１の方法は、診断に重要な値域で必要な精度が得られないため、医用装置として不十分である。【０００８】
本発明は、非侵襲・高精度・高い時間分解能で心内圧を計測する医用撮像装置を提供するものである。
【課題を解決するための手段】
【０００９】
本発明の超音波撮像装置は、超音波造影剤を注入した撮像対象に対して超音波を送受信する超音波探触子と、超音波探触子が受信した撮像対象からの超音波信号を信号処理する信号処理部と、信号処理部の処理結果を表示する表示手段とを備え、信号処理部は、受信した超音波信号の周波数分布を算出する周波数解析部と、周波数解析部によって算出された周波数分布から、受信した超音波信号の減衰率の周波数依存性の、１つ以上の規定の周波数での値、最大値を示す周波数、最大値、半値幅のうち１つ以上である、周波数分布特徴量を検出する周波数分布特徴量検出部と、周波数分布特徴量を圧力に換算する周波数分布特徴量−圧力換算部とを有し、表示手段は、周波数分布特徴量−圧力換算部で算出された圧力の値を表示する。【００１０】
また、本発明の超音波撮像方法は、撮像対象に超音波造影剤を注入する造影剤注入ステップと、撮像対象に超音波信号を送波する超音波送波ステップと、撮像対象から反射された超音波信号を受信し、受信した超音波信号の周波数分布を算出する周波数解析ステップと、算出した周波数分布から、受信した超音波信号の減衰率の周波数依存性の、１つ以上の規定の周波数での値、最大値を示す周波数、最大値、半値幅のうち１つ以上である、周波数分布特徴量を検出する周波数分布特徴量ステップと、検出した周波数分布特徴量を圧力に換算する周波数分布特徴量−圧力換算ステップとを有する。【００１１】
超音波造影剤は、複数の液体成分を内包する多成分液滴型造影剤であり、圧力を検出する部位の温度をＴc、検出しようとする圧力の上限をＰu、下限をＰlとするとき、液体成分のうち少なくとも１つは温度Ｔcでの蒸気圧がＰu以上であり、他の少なくとも１つは温度Ｔcでの蒸気圧がＰl以下であり、複数の液体成分はＰl以上、Ｐu以下の圧力で気液平衡が成り立つモル分率で混合されていることが好ましい。また、このとき、超音波送波は、多成分液滴型造影剤を気化させるための気化用超音波送波を行った後、圧力検出用の超音波を送波する２段階送波を行う。
【発明の効果】
【００１２】
本発明によれば、非侵襲・高精度・高い時間分解能で体内の圧力を計測することができる。より具体的には、従来技術は規定の１周波数における信号強度であるのに対して、本発明は複数の周波数における物理量を用い、サンプル点を多数用いるため、高い精度が得られる。また、従来例は非線形成分である分調波を抽出するために、１回の受信に２回、送信する必要があるが、本発明はその必要がなく、単純計算でも従来例より２倍高い時間分解能が得られる。【００１３】
さらに、本発明の一態様では、２以上の液体成分を内包する多成分液滴型造影剤を体内に注入し、着目した部位を超音波照射して温度を上昇させ、多成分液滴型造影剤が内包した液体を気化させ、気液平衡状態を成立させる。気液平衡状態では、圧力によって気相の体積が変化するが、この体積変化は気泡径の変化として現れ、その結果、受信信号に含まれる気泡の挙動は圧力の影響が増強され、高い精度で圧力を検出できる。
【図面の簡単な説明】
【００１４】
【図１】本発明による超音波撮像装置の構成例を示すブロック図。
【図２】本発明による超音波撮像方法の処理を表すフローチャート。
【図３】受信信号から圧力を算出する処理の詳細を表すフローチャート。
【図４】減衰量の周波数分布を算出する処理の詳細を表すフローチャート。
【図５】生体に超音波を送信・受信する処理の説明図。
【図６】減衰量の周波数分布を算出する処理の説明図。
【図７】超音波造影剤の気泡径の分布を表す図。
【図８】周波数分布特徴量の例、及び、圧力を変えた場合の減衰曲線の例を表す図。
【図９】生体に超音波を送信・受信する処理を説明する図。
【図１０】減衰量の周波数分布を算出する処理のフローチャート。
【図１１】代表的な超音波造影剤の構造を示す概念図。
【図１２】２成分系の一定温度における気液平衡曲線の概念図。
【図１３】２成分系の一定温度における気液平衡曲の例を示す図。
【図１４】本発明による超音波撮像装置の構成例を示すブロック図。
【図１５】硬さを表示する実施例の処理を示すフローチャート。
【図１６】血管壁の硬さを算出する場合を例に、動き量を検出する方法を説明した図。
【図１７】心臓に超音波を送信・受信する処理の説明図。
【図１８】測定部位の形状と圧力を表示する処理の全体の流れを示すフローチャート。
【図１９】表示の一例を示す図。
【発明を実施するための形態】
【００１５】
図１は、本発明による超音波撮像装置の構成例を示すブロック図である。超音波探触子１は複数の素子から構成される。装置本体２は、送信ビームフォーマ３、増幅手段４、受信ビームフォーマ５、信号処理部６、メモリ７、表示手段８、入力手段９、制御部１０を有する。信号処理部６は、受信信号の周波数分布を算出する周波数解析部６１、周波数解析部６１の算出した周波数分布から、周波数分布特徴量を検出する周波数分布特徴量検出部６２、周波数分布特徴量を圧力に換算する周波数分布特徴量−圧力換算部６３を備える。周波数分布特徴量検出部６２で検出する周波数分布特徴量は、減衰率の周波数依存性の、１つ以上の規定の周波数での値、最大値を示す周波数、最大値、半値幅のうち１つ以上である。周波数分布特徴量−圧力換算部６３は、周波数分布特徴量から圧力へ換算するための知識である周波数分布特徴量−圧力換算知識を格納する周波数分布特徴量−圧力換算知識格納部６３１を備える。なお、通常のＢモード表示やドップラ表示などの構成は簡単のために記していない。
【００１６】
ユーザによって超音波造影剤が撮像対象である生体部位に注入され、送信ビームフォーマ３で生成された超音波パルスが超音波探触子１から生体に送信され、超音波造影剤を含む生体部位から反射した超音波を超音波探触子１が受信する。受信信号は増幅手段４に入力されて増幅され、受信ビームフォーマ５が整相加算する。信号処理部６には受信ビームフォーマ５が出力した受信信号が入力され、画像化と圧力の算出を行う。作成された画像及び算出された圧力はメモリ７に格納された後、読み出し・補間されて表示手段８に表示される。なお、これらの処理は制御手段１０によって制御される。
【００１７】
信号処理部６は、周波数解析部６１が受信信号の周波数分布を算出したあと、この周波数分布から周波数分布特徴量検出部６２が周波数分布特徴量を検出し、周波数分布特徴量−圧力換算知識を用いて、周波数分布特徴量−圧力換算部６３が周波数分布特徴量を圧力に換算する。
【００１８】
図２から図４は、本発明の超音波撮像方法の処理の一例を表す処理フロー図である。本発明では、図２に示したように、最初にユーザが撮像対象に超音波造影剤を注入する（Ｓ１１）。その後、超音波探触子１から撮像対象に超音波信号を送波し（Ｓ１２）、増幅手段が撮像対象から超音波信号を受信して受信ビームフォーマで整相加算し、この受信信号を信号処理部６が処理して圧力を算出する（Ｓ１３）。算出された圧力は前記表示手段に表示される（Ｓ１４）。
【００１９】
ここで、本発明で使用する超音波造影剤はマイクロバブルあるいはナノバブルと呼ばれる直径数１０μｍ以下の微小気泡で、生分解性高分子から構成される殻を有してもよい。
【００２０】
図３に、図２のステップ１３における信号処理部６の処理の詳細を示す。信号処理部６では、受信深度の異なる受信信号をもとに減衰量の周波数分布を算出し（Ｓ２１）、周波数分布から周波数分布特徴量を検出する（Ｓ２２）。その後、圧力をパラメータとした測定ノイズのない理想的な減衰量の周波数分布である、次式（１）に示す周波数分布特徴量−圧力換算知識を用いて周波数分布特徴量を圧力ｐ0に変換する（Ｓ２３）。変換の方法としては、たとえば、周波数分布特徴量を入力し、周波数分布特徴量−圧力換算知識を回帰曲線として最小二乗法でパラメータである圧力を求める方法による。
【数１】
2010147055.tif 000003
【００２１】
ここで、σeは散乱断面積、ωは式（１）の変数である周波数で、共鳴周波数ω0で規格化してΩで表す。ａeは気泡径の平衡値、ρLは周りの流体の密度、κはポリトロープ指数、ＧSは気泡の殻の剛性率、ｄSeは式（３）で定義される殻の厚み、ηSは気泡殻のせん断粘度を表し、造影剤及び血液などの流体の物性値によって決まる。これらは式（１）においては定数であり、あらかじめ調べられて装置内に記憶されているものとする。ｐ0は静水圧で、（１）における唯一のパラメータであり、本発明の装置及び方法で検出すべき圧力である。
【００２２】
図４は、図３のステップ２１における減衰量の周波数分布を算出する処理の詳細を示した図である。まず、受信深度の異なる受信信号Sigi(t)と受信深度ｚi｛ｉ＝１，２，…，ｎ｝を入力する（Ｓ３１）。次に、受信信号Sigi(t)｛ｉ＝１，２，…，ｎ｝の周波数分布|Sigi(f)|＝|Ｆ［Sigi(t)］|を計算する（Ｓ３２）。最後に、受信信号の周波数分布|Sigi(f)|から減衰量の周波数分布
σ(f)＝|Sigi(f)|／|Sigj(f)|
を計算する（Ｓ３３）。ここでＦ［・］はフーリエ変換を表す。
【００２３】
次に、図５から図８を用いて、本発明の超音波撮像装置による撮像方法の詳細について説明する。図５は、超音波探触子が生体に超音波を送信・受信する様子を示す図である。超音波探触子１から生体７２に超音波が送信されると、臓器７３に含まれる超音波造影剤７４から反射された信号７５，７６を受信する。信号７５の受信深度はｚ1であり、信号７６の受信深度はｚ2である。臓器７３は血管や心臓の心腔のように、超音波造影剤より反射輝度がいちじるしく低い流体を含む臓器でもよいし、肝臓のように視認できる程度の反射輝度を持つ柔組織から成る臓器でもよい。後者の場合は、信号処理部は、周波数解析部が、超音波造影剤による減衰のみを処理できるように、超音波造影剤がない場合に、受信信号７５，７６と同様の信号を受信し、超音波造影剤以外に起因する減衰量を差し引く処理を行えばよい。なお、超音波造影剤以外に起因する減衰量は、一般的に、周波数に比例する量であり、たとえば超音波造影剤の注入前に超音波を送受信・解析して求めることができる。
【００２４】
図６は、減衰量の周波数分布を算出する図３のステップ２１の処理を説明した図である。分布８１は受信深度ｚiからの受信信号Sigi(t)をフーリエ変換した信号Sigi(f)＝Ｆ［Sigi(t)］を示し、分布８２は受信深度ｚjでの信号強度の周波数分布である。図４のステップ３２では、複数の受信深度でこのような周波数分布８１，８２を計算する。ただし、受信深度ｚiからの受信信号Sigi(t)は、ｃを音速を表す定数とし、ｄｔを周波数の計算に必要なサンプル数を含む時間とし、時刻は変換撮像対象に超音波を送波した時刻をｔ＝０とするとして、受信深度ｚiに相当する時間（ｔ＝ｚi／ｃ−ｄｔからｔ＝ｚi／ｃ＋ｄｔ）に受信した信号とする。
【００２５】
図６（ｂ）に示した分布は、信号強度の周波数分布８１と８２の比をとって得た、次式で表される減衰量の周波数分布σ(f)である。図４のステップ３３では、この減衰量の周波数分布σ(f)を計算する。ただし、図６ではｄＢ表示で示した。
【００２６】
σ(f)＝|Sigi(f)|／|Sigj(f)|
ここで、式（１）で表される周波数分布特徴量−圧力換算知識について説明する。超音波造影剤は、体内に注入される数μｍ前後の微小気泡で、超音波を照射されると、泡の外部が液体、内部が気体であり、密度差が大きいことから、大きな反射率で超音波信号を反射する。通常、超音波造影剤から反射された信号は、基本波成分及び高周波成分の信号強度が画像化される。このように通常の使い方では、気泡が存在していることを使って、血管壁や肝腫瘍の境界の可視化を行っている。しかし本発明では、気泡の力学的挙動から気泡の存在以上の情報を抽出することを目標とし、気泡の運動を記述する方程式の性質を概観した。気泡の運動は、気泡半径ａ(t)の運動方程式である次のChurch方程式によって記述される。
【数２】
2010147055.tif 000004
【００２７】
ここで、ａ1は気泡の内径、ａ2は気泡の外径、ａ1eは内径の平衡値、ａ2eは外径の平衡値、ρLは周りの流体の密度、ρSは気泡殻の密度、ηLはまわりの流体のせん断粘度、ηSは気泡殻のせん断粘度、ＶSは気泡核の体積、ｐ∞(t)は気泡から十分離れた場所での圧力、ｐ0は静水圧で、検出すべき圧力である。
【００２８】
泡の挙動を記述する方程式は仮定の置き方によって各種あるが、生体内の残存時間を長くするために、２０００年以降は殻のついた気泡が主流であるため、殻のついた気泡の方程式であるChurch方程式を用いた。ここで、式（３）の近似と式（４）に示す変数変換により、気泡径の変位量ｘ(t)を変数とする。ａは気泡径、ａeは気泡径の平衡値である。
【数３】
2010147055.tif 000005
【００２９】
さらに線形化を行って式（５）を得ると、その解として式（６）が得られる。
【数４】
2010147055.tif 000006
【００３０】
ここで、ｐi(t)は超音波パルスの照射による圧変化である。
【００３１】
減衰率の周波数分布は、散乱断面積の式（７）から式（８）と求められる。
【数５】
2010147055.tif 000007
【００３２】
式（８）中のｎ(a)は気泡径の分布である。気泡径は例えば図７に示すような分布をしているが、図４のステップ３３に関する上記の説明では分布の幅が０、すなわち撮像対象内部の気泡は全て同じ値を持つとして説明を簡略化した。実際にはｎ(a)を反映する必要があり、たとえば体内注入前の気泡径の分布を測定しておき、これを用いてもよい。
【００３３】
上記の結果のうち、定性的な把握に必要なものを以下にまとめた。
【数６】
2010147055.tif 000008
【００３４】
すなわち、共鳴曲線（６）の形状は共振周波数ω0と分散の大きさδでほぼ決まる。減衰曲線(７)（８）の形もω0、δ、超音波パルスの周波数分布ｐi(ω)でほぼ決まる。ω0、δを表す式（９）の右辺のうち、ｐ0以外は超音波造影剤の構造（物性定数）によって固定される、測定においては既知とみなされる値だが、これらの物性定数のうち実質的に可変なのは殻の剛性Ｇs及び気泡の平衡径ａeのみである。共鳴曲線（６）及び減衰曲線（７）（８）の形状、すなわち圧ｐ0の感知力はこれら２つの値によって調節でき、気泡殻の剛性率Ｇsが小さく、気泡の平衡径ａeが小さい超音波造影剤が圧検出には適していることがわかる。
【００３５】
参考のために、商品化されている超音波造影剤の物性定数を表１に示す。
【表１】
2010147055.tif 000009
【００３６】
図８に、周波数分布特徴量の例、及び、気泡の平衡径ａeを一定にして、圧力を３通りの値に変えた場合の減衰曲線の一例を示した。減衰曲線は圧力によって有意に変形している。この変形を検出するための周波数分布特徴量として、図８では共鳴周波数ω0、分散δ-1、規定の周波数ω1における減衰強度、規定の２個の周波数ω1，ω2における減衰強度の傾きを示した。
【００３７】
上記の周波数分布特徴量は全て式（１）の周波数分布特徴量−圧力換算知識に入力することで出力として圧力ｐ0を与える。これらの周波数分布特徴量は、共鳴周波数と分散など、異なる物理量であっても同時に式（１）に入力可能であり、入力する値の数が多いほどｐ0の検出精度を向上させることができる。
【００３８】
より詳しくは、規定の周波数ωiにおける減衰強度σi （ｉ＝１，２，…，ｎ）が得られた場合、（７）を回帰式とする回帰問題として圧ｐ0を求めることができる。すなわち、最小二乗法により
J＝Σi＝1〜n|σe（ωi;ｐ0）-σi|2を最小にするｐ0、別の表現では
δJ＝∂/∂ｐ0   [Σi＝1〜n|σe（ωi;ｐ0）-σi|2]
＝０
を満たすｐ0を数値的に求めて圧力の検出値とする。気方径の分布n(ae)の影響を考慮する場合は（８）を回帰式とし、上記のσe（ωi;ｐ0）をΣe（ωi;ｐ0）で置き換えればよい。
【００３９】
さらに、規定の周波数ωiにおける減衰強度σi（ｉ＝１，２，…，ｎ）の他にも共鳴周波数ω0、分散δ-1、規定の２個の周波数ω1，ω2における減衰強度の傾きが得られた場合には、これらを規定の周波数ωiにおける減衰強度σi（ｉ＝ｎ＋１，ｎ＋２，…，ｎ＋ｍ）に換算して、回帰問題の入力値の数を増やし、ｐ0を求めることができる。又は、これらのうち１つのみが得られた場合、回帰問題は１つの式に帰着される。たとえば共鳴周波数ω0のみが得られた場合、回帰問題は式（7’-2）に帰着され、式（7’-2）によりω0をｐ0に変換すればよい。同様に、分散δのみが得られた場合、回帰問題は式（7’-2）及び式（7’-3）に帰着され、式（7’-2）及び式（7’-3）によりδをｐ0に変換できる。
【００４０】
以上のような構成によれば、非侵襲・高精度・高い時間分解能で体内の圧力を計測する医用撮像装置を提供することができる。特に、複数の周波数における信号強度を用い、サンプル点を多数用いるため、ＳＮ比が高くなり、高い精度が得られる。また、分調波の抽出と異なり、１つの空間位置での圧力を検出するために２回以上、送信する必要がないため、高い時間分解能が得られる。さらに、本発明では１つ以上の物理量を検出して血圧に換算することで、精度を向上させることができる。
【００４１】
次に、図９及び図１０を用いて本発明の第２の実施例を説明する。第２の実施例の方法は、図２のステップ１３の処理において、周波数分布特徴量−圧力換算知識として式（１０）に示す共鳴曲線を用いる例である。
【００４２】
図９は、本実施例において、超音波探触子が生体に超音波を送信・受信する処理を説明した図である。超音波探触子１から生体７２に超音波が送信されると、生体に含まれる超音波造影剤７４から反射された信号７５を受信する。臓器７３は血管や心臓の心腔のように、超音波造影剤より反射輝度がいちじるしく低い流体を含む臓器でもよいし、肝臓のように視認できる程度の反射輝度を持つ柔組織から成る臓器でもよい。後者の場合は、図５で説明したような処理を行って超音波造影剤に起因する減衰量のみを抽出すればよい。
【００４３】
図１０は、本実施例の場合、図２のステップ１３において実行する信号処理部６が受信信号から圧力を算出する処理を説明するフローチャートである。ある受信深度での受信信号を入力してその周波数分布を算出し（Ｓ４１）。次に、その周波数分布から、基本波、高調波、分調波、差周波、和周波の強度の周波数分布、値の分岐の大きさのうち１つ以上である周波数分布特徴量を検出し（Ｓ４２）、その後、圧力をパラメータとした測定ノイズのない理想的な共鳴曲線である次式の周波数分布特徴量−圧力換算知識を用いて周波数分布特徴量を圧力に変換する（Ｓ４３）。
【数７】
2010147055.tif 000010
【００４４】
たとえば、基本波の強度の周波数分布ｘB(ω)を測定した場合、ｘB(ω)は共鳴曲線ｘ(ω)になっているはずである。したがって測定値ｘB(ω)を入力値、（１０）を回帰式として最小二乗法を行えばｐ0が求まる。同様にして、ｎ次高調波の強度の周波数分布ｘn(ω)、１／ｍ分調波の強度の周波数分布ｘ1/ｍ(ω)、差周波の強度の周波数分布ｐｍ-n(ω)、和周波の強度の周波数分布ｐｍ+n(ω)を測定した場合は、それぞれ回帰式（１０）のｐi(ω)を、造影剤がない場合に想定されるｎ次高調波、１／ｍ分調波、差周波、和周波の強度の周波数分布として回帰式（１０）によってｐnを計算する。たとえば基本波と高調波など、１つ以上の周波数特徴量が得られた場合は、第１の実施例と同様に、回帰問題の入力数が増えたことになり、より正確にｐ0を求められる。
【００４５】
以上のような構成によれば、構造が一様である領域が小さく、減衰をＳ／Ｎよく検出できない毛細血管や臓器中の超音波造影剤の信号からも、圧力を求めることができる。
【００４６】
次に、図１１から図１３を用いて本発明の第３の実施例を説明する。
【００４７】
第１及び第２の実施例では、式（９）を用いて述べたように、超音波造影剤の信号は、圧力ｐ0と気泡の平衡径ａeを反映しているが、特に西暦２０００年以降、体内の残存時間を長くするために用いられるようになってきた殻のついた気泡、式（９）で言えば、気泡殻の剛性Ｇsが大きい場合、気泡の挙動、ひいては共鳴曲線及び減衰曲線は、外圧ｐ0よりも気泡の平衡径ａeをより強く反映する。そこで、複数の成分を液体状態で含む超音波造影剤を用いることで、気泡の平衡径ａeを外圧ｐ0の俊敏な関数にさせて、圧力の検出精度を向上させた例が以下で述べる第３の実施例である。
【００４８】
第３の実施例では、超音波造影剤の組成が異なり、また、気泡の平衡径ａeと外圧ｐ0の関係
ａe＝ｆ（ｐ0）      (11)
を、式（１）あるいは式（１０）とあわせて用いる。関数ｆ(ｐ0)は実験から求め、あらかじめ装置内に格納しているとする。ただし、特に気液２相共存状態以外では、ファンデルワールスの状態方程式と構成成分のモル分率から理論的に求めて装置内に格納していてもよい。
【００４９】
図１１は、代表的な超音波造影剤の構造を示した概念図である。図１１（ａ）は、内部に気体を含み、殻を有さない、気泡径ｄが１μｍ以上の超音波造影剤である。このような超音波造影剤はすでに商品化されており、第１世代造影剤と呼ばれ、生体内での残存時間が短い。図１１（ｂ）は、内部に気体を含み、殻１４１を有し、気泡径ｄが１μｍ以上の超音波造影剤である。このような超音波造影剤はすでに商品化されており、第２世代造影剤と呼ばれ、生体内での残存時間が長い。殻１４１の厚みｔは径ｄの数％以下であり、親水性の部位１４２と疎水性の部位１４３を持つ高分子をその構成単位とする。図１１（ｃ）は、内部に液体を含み、殻１４４を有し、気泡径ｄが１μｍ以下の超音波造影剤である。このような超音波造影剤はまだ商品化されておらず、研究段階にある。相変化液滴型造影剤と呼ばれる。毛細血管より先など、浸透可能な領域が広く、超音波照射によって意図した部位で気化し、撮像に用いられる。第３の実施例は、複数種類の液体を含む相変化液滴型造影剤を用いる例である。
【００５０】
現在研究されている相変化液滴型造影剤としては日立製作所の川畑らのものがあり、物性値は以下のようである。
【００５１】
殻の材質    ：Lipid / surfactant emulsion
液滴の直径  ：200nm
内包物      ：PEOB
相変化型造影剤を用いる第３の実施例では、図２のステップ１１において、気液相変化型の複数成分の造影剤を撮像対象に注入する。また、ステップ１２の撮像対象に超音波信号を送波する処理において、まず気化用の超音波を撮像対象に送波して造影剤内部の液体を気化し、その後に圧力検出用の超音波を撮像対象に送波する。ただし、造影剤が内包する成分の蒸気圧が測定したい圧力の値域より低い場合、気化用の超音波を撮像対象に送波するこのステップはなくてもよい。送波する場合、気化用の超音波を送波する空間位置は、圧検出用の超音波を送受信する空間位置と同じ、あるいはそれより血流が上流の位置とする。たとえば、心臓の左心室内の圧力を検出する場合、左心室内あるいは左心房内に気化用の超音波を送波すればよい。また、肝臓内の圧力を検出する場合には、肝動脈に気化用の超音波を送波すればよい。
【００５２】
次に、図１２及び図１３を用いて、複数の成分を液体状態で含む超音波造影剤を用いることで、気泡の平衡径ａeを外圧ｐ0の俊敏な関数にできることを説明する。図１２（ａ）は、１成分系の等温曲線、すなわちファンデルワールスの状態方程式で、図１２（ｂ）は、２成分系の一定温度における気液平衡曲線の概念図である。
【００５３】
図１２（ａ）の縦軸は圧力、横軸は体積である。図の下に記したファンデルワースの状態方程式に従う変化を点線で示した。図中の実線は実際の流体が辿る体積−圧力変化である。ただし温度は一定値を仮定している。実際の流体の体積−圧力はファンデルワースの状態方程式と異なり、単調に起こる。状態方程式と異なる部分は、図中の面積Ｓ１とＳ２が等しくなるように引いた直線部で、体積は変化するが圧力が一定であるこの直線部は、気液相変化における気液２相共存状態にあたり、このときの圧力を蒸気圧Ｐｃと呼ぶ。気液２相共存状態より圧力が小さい線上は気体状態、大きい線上は液体状態である。このように１成分系では、温度を一定にして圧力を変化させると、蒸気圧Ｐｃにおいて体積が不連続な飛びをもって大きく変化するが、多成分系では、温度を一定にして圧力を変化させると、体積は常に連続して変化する。このため、圧力変化を体積変化によって定量できる。
【００５４】
多成分系の気液相変化における体積変化の様子を図１２（ｂ）に示す。図１２（ｂ）では２成分系を仮定した。縦軸は圧力、横軸は第１成分のモル分率である。Ｐｃ１が成分１の蒸気圧、Ｐｃ２が成分２の蒸気圧で、太実線が液相線、細実線が気相線を意味する。液相線より圧力の大きい領域では液体状態、細実線より圧力が小さい領域では気体状態、その中間の、斜線部にあたる圧力では気液２相共存状態にある。図１２（ａ）では気液２相共存状態の圧力は１つの値Ｐｃに決まったのに対し、２成分系では、斜線部が縦軸Ｐに対して幅を持つことからわかるように、気液２相共存状態での圧力は成分１の蒸気圧Ｐｃ１から成分２の蒸気圧Ｐｃ２まで幅を持つ。気体状態から圧力を増加していくと、準静過程とみなせるほど変化速度が遅い場合、鎖線のような圧力‐モル分率変化を経て蒸気圧の低い成分である成分１のモル濃度が増加し、凝縮される。本発明の超音波撮像装置及び方法における変化は実践矢印のような変化とみなすことができる。圧力Ｐと気泡径ａeの関係ｐ0＝ｆ(ａe)の式（１１）は、気液２相共存状態（斜線部）以外ではファンデルワールスの状態方程式（ただし関係Ｖ＝πａe3/４により体積Ｖを気泡径ａeに変換する）でよく表される。気液２相共存状態（斜線部）での関係は実験により求めた値を使うとする。
【００５５】
圧力変化の定量精度は、２成分気液混合相の圧力方向の厚み１８が大きいほど高精度になる。２成分気液混合相の圧力方向の厚みは、成分の物性値、具体的には、活性係数と蒸気圧で決まる。モル分率が０．５に近いほど厚みが厚いのはもちろんである。すなわち、多成分系の気液平衡における圧−体積変化を活用することは、測定を意図した圧力範囲において、圧依存性を増強させる効果を持つ。
【００５６】
厚み１８の大きさに影響する活性係数は、２成分の混合溶液において、各成分の相互作用の度合いを表す物性値で、値が１以上の２成分を造影剤の要素とすることが望ましい。１以上の場合の図１２（ｂ）の具体例（メタノール−水系）を図１３（ａ）に、１以下の例（tertブタノール−secブタノール系）を図１３（ｂ）に示した。図１３（ａ）の方が、２成分気液混合相の圧力方向の厚みが大きくなっている。図１３（ａ）、図１３（ｂ）は共に温度が２５度の例であるが、図１３（ａ）では４０ｍｍＨｇから８０ｍｍＨｇの範囲で高精度な圧計測が可能である。
【００５７】
内包物の活量係数をなるべく大きく設定することを考えるとき、相変化型造影剤の内包物としては、Ｃ3Ｆ3，Ｃ3Ｆ8，Ｃ5Ｆ12，ＳＦ6で比べるなら、Ｃ3Ｆ3がよい。そのほか、モル分率の０．５近傍への調節も有効である。
【００５８】
以上のように内包成分を選択すれば、気泡の平衡径ａeを外圧ｐ0の俊敏な関数（１１）にできる。本発明の第3の実施例では、第1あるいは第2の実施例と同様の処理において、それぞれ（１）あるいは（１０）に（１１）を代入した式を回帰式として、ｐ0を求める。  
以上のような構成によれば、体内の圧力を高精度で検出することができる。
【００５９】
次に、図１４から図１６を用いて本発明の第４の実施例を説明する。本実施例は、第１から第３の実施例の処理を含み、超音波を送受する空間位置の圧力のほかに動き量も検出し、圧力と動きから超音波を送受する空間位置の硬さを算出する例である。
【００６０】
図１４は、本実施例の装置構成を示すブロック図である。本実施例の信号処理部６は、動き量検出部６４及び動き量・圧力−硬さ変換部６５を備える。
【００６１】
図１５は、本実施例の処理の流れを示すフローチャートである。第１から第３の実施例で述べた方法で、撮像対象の圧力を算出する（Ｓ５１）。それと同時に、撮像対象の動き量も検出する（Ｓ５２）。動き量の検出方法は既存のどのような方法でもよい。例えば、Ｂ像のスペックルパターンを追跡して動き量とするティッシュトラッキング法などでもよい。その後、圧力と動きから硬さを算出し（Ｓ５３）、硬さを表示する（Ｓ５４）。硬さは圧力に比例し、動き量に逆比例する量で定義される。
【００６２】
図１６は、血管において血管壁の硬さを算出する場合を例に、動き量を検出する方法を説明した図である。図１６（ａ）は血管の全体図である。本発明の方法は、最初に血管１０１を撮像する。撮像断面は径方向断面１０２でも軸方向断面１０３のどちらでもよい。図１６（ｂ）（ｃ）は径方向断面１０２をとった場合、図１６（ｄ）（ｅ）は軸方向断面１０３をとった場合の撮像結果を表す。図１６（ｂ）と図１６（ｃ）、図１６（ｄ）と図１６（ｅ）は隣接する時刻における撮像結果を示している。
【００６３】
図１６（ｂ）（ｃ）について説明すると、図１６（ｂ）は時刻ｔ、図１６（ｃ）は時刻ｔ＋ｄｔにおける撮像結果である。時刻ｔでの位置１０４及び時刻ｔ＋ｄｔでの位置１０４’での圧力が算出されているとする。このとき、ｔ＝ｔでの撮像結果の中に微小領域１０５を設定し、ｔ＝ｔ＋ｄｔで同じ形状を持つと判断される微小領域１０５’を検出し、微小領域１０５と１０５’の重心の移動量１０６を動きとする。また、動き量の精度が低くてよい場合には、たとえば血管径ｄ１とｄ２の差の１／２を大きさ、血管壁の法線方向を方向とするベクトルを動き量としてもよい。図１６（ｄ）及び図１６（ｅ）に示した軸方向断面１０３においても同様に動き量１０６が算出できる。
【００６４】
以上のような構成によれば、非侵襲に、高精度で、生体の硬さを検出することができる。
【００６５】
次に、第１から第４の実施例の処理を心臓に適用した場合の例を、図１７及び図１８を用いて説明する。
【００６６】
図１７は、超音波探触子が心臓に超音波を送信・受信する処理を説明した図である。超音波探触子１から生体７２に超音波を送信し、心臓７３に含まれる超音波造影剤７４から反射された信号７５１，７５２を受信する。
【００６７】
処理の全体の流れは図１８に示した。最初に超音波像により心臓の形状を撮像し（Ｓ６１）、圧力を測定する部位を検出する（Ｓ６２）。次に、測定部位の圧力を実施例１から実施例４の方法で検出し（Ｓ６３）、測定部位の形状と圧力を表示する（Ｓ６４）。
【００６８】
表示の例を図１９に示す。図１９（ａ）は心臓内部の血圧を場所ごとに検出し、形態画像に重ねて表示した例、図１９（ｂ）は心臓の形態画像及び血圧の検出値を解析した結果を表示した例である。図１９（ａ）（ｂ）とも９１は表示画面、９２は撮像領域、９３は撮像対象の形態画像で、図１９（ａ）において、９４は圧力の検出結果を、異なる検出値を異なる色で表示した例である。図１９（ｂ）において、９５は解析結果の表示エリアで、横軸を形態画像を画像処理して抽出した、心臓の左心室の体積ＶLV、縦軸を本発明の撮像方法によって検出した左心室の圧力の空間的な代表値ＰLVとした例を示した。グラフ内に示した９６は (ＶLV, ＰLV)の時間変化である。測定条件の違いにより、複数のループが描かれている。９７は複数のループの接線の傾きで心機能の指標の１つであるＥMAX、９８は前記接線の切片で心機能の指標の１つであるＶ0である。図１９では撮像対象が心臓である例を示したが、撮像対象は心臓に限らない。なお、図１９は静止画でも動画でもよい。
【００６９】
以上のような構成によれば、心臓など撮像対象の内部の血圧を非侵襲に高精度に検出することができ、形態情報とあわせて臨床的に意味のある量をユーザに提示することができる。
【符号の説明】
【００７０】
１：超音波探触子
２：装置本体
３：送信ビームフォーマ
４：増幅手段
５：受信ビームフォーマ
６：信号処理部
６１：周波数解析部
６２：周波数分布特徴量検出部
６３：周波数分布特徴量−圧力換算部
７：メモリ
８：表示手段。
【図１】
2010147055.tif 000011
【図２】
2010147055.tif 000012
【図３】
2010147055.tif 000013
【図４】
2010147055.tif 000014
【図５】
2010147055.tif 000015
【図６】
2010147055.tif 000016
【図７】
2010147055.tif 000017
【図８】
2010147055.tif 000018
【図９】
2010147055.tif 000019
【図１０】
2010147055.tif 000020
【図１１】
2010147055.tif 000021
【図１２】
2010147055.tif 000022
【図１３】
2010147055.tif 000023
【図１４】
2010147055.tif 000024
【図１５】
2010147055.tif 000025
【図１６】
2010147055.tif 000026
【図１７】
2010147055.tif 000027
【図１８】
2010147055.tif 000028
【図１９】
2010147055.tif 000029
【国際調査報告】
2010147055.tif 000030
2010147055.tif 000031
2010147055.tif 000032
2010147055.tif 000033