JP5697716B2 - 超音波測定装置 - Google Patents

Description

本発明は、超音波を用いて測定対象物の断層像若しくは三次元像を取得する超音波測定装置に関し、特に、試料の歪特性や弾性特性を取得する超音波測定装置に関する。

従来の一般的な超音波を用いて断層像を得る装置は、超音波を測定対象物（試料）に送信する送信部と、反射波を受信する受信部と、送受信波を走査するための走査手段と、受信した反射信号を輝度信号に変換、可視化するための手段によって構成されている。そして、これらの手段により取得された時系列断層画像を用いて試料の内部を観察することが行われている。また、上記装置の一つの形態においては、上記走査手段によって超音波を上下左右に走査し三次元像を得ることも行われている。

最近では、試料の表面を圧迫することによって人為的に試料内部に歪を生じさせ、超音波信号の処理によって歪を計測する方法や、ヤング率等の弾性係数を画像化する方法も検討されてきている。

特許文献１においては、試料を縦方向に圧迫し各点の相対位置変化を測定することによって試料の弾性特性を測定することが記載されている。また、弾性体である試料が自ら動きを持つ場合には、上記圧迫動作を行うことなく超音波信号の処理のみによって歪を計測することが可能となっている。

また、特許文献２には、心筋の厚み変動を心臓壁の２点の相対位置変動に基づいて計測する方法が開示されている。

米国特許第５，１７８，１４７号明細書 特開平９−３１３４８６号公報

前記特許文献１及び２において開示されている技術の共通する点は、一つの超音波プローブによって得られた超音波信号を演算することにより試料内各点の位置ずれ量（変位）を測定していることである。そして、前述の文献によれば、試料の圧迫及び超音波信号の測定は一方向に関してのみ行われている。

試料の変位や歪が、圧迫方向（例えば、超音波の送受信方向）にのみ生ずる場合には、前述のように一つの超音波プローブによる測定で精度の高い情報が得られるが、実際には、試料は圧迫方向と異なる方向にも変位したり、歪む場合が多い。

従って、上記従来例を用いて圧迫方向と異なる方向にも歪む試料を測定した場合等には、当該試料の変位や歪みに関して十分な情報が得られていなかった。

そこで、本発明は、測定対象物に関して、一方向の変位や歪みに関する情報のみならず、複数方向の歪みに関する情報を出力することのできる超音波測定装置を提供することを目的とする。

上記目的を達成するために本発明は、以下の構成を採用する。

本発明は、
超音波測定装置であって、
第１の超音波プローブと第２の超音波プローブを少なくとも含む複数の超音波プローブと、
測定対象物に加わる力もしくは測定対象物の変形量の時間変化が、増加傾向を示す増加期間であるか、または前記時間変化が減少傾向を示す減少期間であるかに関する情報を取得する情報取得部と、
演算処理部と、
を備え、
前記第１の超音波プローブは、該第１の超音波プローブの超音波送信領域と前記第２の超音波プローブの超音波送信領域との重複箇所に位置する前記測定対象物によって反射された超音波を受信して、第１の受信データを生成し、
前記第２の超音波プローブは、前記重複箇所に位置する前記測定対象物によって反射された超音波を受信して、第２の受信データを生成し、
前記演算処理部は、
前記第１の受信データの内、前記増加期間あるいは前記減少期間のいずれか一方の期間におけるデータを用いて、前記測定対象物の第１の方向に関する第１の変位関連情報を演算し、且つ
前記第２の受信データの内、前記第１の変位関連情報を演算するために用いた期間と同じ期間におけるデータを用いて、前記測定対象物の第２の方向に関する第２の変位関連情報を演算し、
前記第１の変位関連情報と、前記第２の変位関連情報と、
前記第１及び第２の超音波プローブの相対的な位置関係情報と、を用いて多次元の変位関連情報を生成することを特徴とする超音波測定装置を用いる。

本発明によれば、複数方向の変位や歪みに関する情報を出力することのできる超音波測定装置が提供される。

特に、本発明は、第１及び第２の変位関連情報を共に増加期間のデータ（あるいは共に減少期間のデータ）から算出するため、測定対象物の弾性特性をより正確に反映した情報を得ることができる。

本発明の実施形態に係る超音波測定装置の構成を示す図。 測定対象物に加わる力又は測定対象物の変形量の、増加期間と減少期間を説明するための図。 本発明の実施例１、２、３に係る超音波測定装置の構成を示す図。 実施例１，２における、プローブの構成を示す図。 プローブの圧迫動作と取得信号を説明する図。 プローブ面と視野、測定点の関係を表した図。 実施例３における、プローブ及び圧迫手段の構成を示す図。 実施例３における、プローブ及び圧迫手段の配置を示す図。 実施例１における圧迫動作と超音波データの取得タイミングを説明するための図。 複数のプローブの送受信タイミングを説明する図。 プローブの送信タイミングを切り替える制御のフロー図。

まず、本発明の実施形態を説明するに際して、使用単語を定義する。

「変位関連情報」とは、測定対象物の変位（位置の変化）そのもの、あるいは変位から算出される情報であり、「変位から算出される情報」には、たとえば、測定対象物の歪み、弾性係数、硬さ（柔らかさ）などが含まれる。これらの変位関連情報は、測定対象物の硬さ（変形しにくさ）に依存する量ということもできる。

次に、本発明の実施形態に係る超音波測定装置を図面を用いて説明する。

本実施形態では超音波測定装置の一例として、医用超音波診断装置の例を示しているが、本発明は、医用超音波診断装置に限らず、弾性体計測のための探傷用装置などのその他の装置にも適用できることはいうまでもない。

図１は、本実施形態に係る超音波測定装置１０００を説明するための模式図である。

１０１０は第１の超音波プローブであり、１０２０は第２の超音波プローブである。各超音波プローブは複数の振動子を含み構成される。１０１５は第１の超音波プローブ１０１０の超音波送信領域を示し、１０１１は第１の超音波プローブ１０１０の超音波の音線方向を示している。また、１０２５は第２の超音波プローブ１０２０の超音波送信領域を示し、１０２１は第２の超音波プローブ１０２０の超音波の音線方向を示している。図１の例では、音線が交差し、且つ超音波送信領域１０１５と１０２５とが一部で重複するように、第１及び第２の超音波プローブ１０１０、１０２０が配置されている。超音波測定装置１０００は、超音波送信領域１０１５と１０２５の重複箇所に位置する測定対象物１０８９（たとえば、腫瘍、しこり、臓器など）について、変位関連情報を測定することができる。

１０４０は、測定対象物１０８９に加わる力あるいは測定対象物の変形量の時間変化が増加傾向を示す増加期間であるか、あるいは前記時間変化が減少傾向を示す減少期間であるかに関する情報を取得する情報取得部である。

図２は、情報取得部が出力する情報（信号）の一例を示している。図２の縦軸は測定対象物に加わる力もしくは測定対象物の変形量を示し、横軸は時間を示している。なお、測定対象物の変形量は、測定対象物上に設定された点（たとえば重心、端点など）や境界線の変位量として捉えてもよいし、測定対象物の高さや幅の変化量として捉えることもできる。図２におけるｔ_０からｔ_１が増加期間、ｔ_１からｔ_２が減少期間、ｔ_２からｔ_３が増加期間である。変位関連情報を測定する際は、超音波プローブあるいは別途設けられた圧迫手段によって、人体１０９９の外側から測定対象物１０８９に対し間接的に周期的な力を加えることで、測定対象物１０８９を周期的に変形させる。ただし、心臓の心筋運動などのように外部からの加圧力が無くても周期的な変形をしている物を測定する場合には、外部から力を加える必要はない。

１０６０は、第１及び第２の超音波プローブ１０１０、１０２０に、超音波を送信するための信号を入力する送信回路である。１０９０は、第１及び第２の超音波プローブ１０１０、１０２０が超音波を受信することにより生成する信号を受信する受信回路である。
受信回路１０９０は、超音波プローブで受信した反射信号を増幅し、各素子の信号に対して適当な遅延量制御を行った後、整相加算を行う。１０５０は、後述する演算処理を行なうための演算処理部である。１０８０は、演算処理部１０５０により得られた情報を出力するための出力部であり、例えばディスプレイである。１０７０は、送信回路１０６０と受信回路１０９０と演算処理部１０５０と出力部１０８０とを制御するためのシステム制御部である。システム制御部１０７０は、例えば送信回路１０６０と受信回路１０９０の送受信タイミングの制御や送受信周波数の制御を行うものである。

第１の超音波プローブ１０１０は、超音波送信領域１０１５、１０２５の重複箇所に位置する測定対象物１０８９によって反射された超音波を受信して、第１の受信データを生成する。また、第２の超音波プローブ１０２０は、測定対象物１０８９によって反射された超音波を受信して、第２の受信データを生成する。第１及び第２の受信データは演算処理部１０５０にそれぞれ入力される。

演算処理部１０５０は、前記第１の受信データの内、増加期間あるいは減少期間のいずれか一方の期間におけるデータを用いて、測定対象物１０８９の第１の方向に関する第１の変位関連情報を演算する。更に、演算処理部１０５０は、第２の受信データの内、第１の変位関連情報を演算するために用いた期間と同じ期間におけるデータを用いて、測定対象物１０８９の第２の方向に関する第２の変位関連情報を演算する。たとえば、演算処理部１０５０は、図２のｔ_０からｔ_１の増加期間において得られた受信データを用いて第１の変位関連情報を演算したら、第２の変位関連情報もｔ_０からｔ_１の増加期間において得られた受信データから演算するのである。具体的には、時刻ｔ_０からｔ_１の範囲内の互いに異なる時刻に受信した受信データから変位を求めることができる。そして、必要に応じて、変位から歪みを求めたり、当該測定対象物に加わる力が分かれば弾性係数を求めることもできる。ここで、第１の方向は、第１の超音波プローブ１０１０から発信される超音波の進行方向であり、典型的には音線方向１０１１に一致する。また、第２の方向は、第２超音波プローブ１０２０から発信される超音波の進行方向であり、典型的には音線方向１０２１に一致する。なお、前記第１及び第２の変位関連情報は互いに共通する傾向を示す期間内同士で求めることができる。つまり、「第１の変位関連情報を演算するために用いた期間と同じ期間におけるデータ」には、前記第１の変位関連情報を時刻ｔ_０からｔ_１の範囲で求めて、前記第２の変位関連情報を時刻ｔ_０からｔ_１に含まれる期間の一部のデータを用いて取得する場合をも含む。また、増加期間と減少期間との境界に位置する時刻は、両方の期間に含まれるものとする。

そして、出力部１０９０は、前記第１及び第２の変位関連情報を用いて測定対象物１０８９の情報を出力する。具体的には、出力部１０９０は、第１及び第２の変位関連情報を同時に出力することもできるし、選択的に一方の変位関連情報のみを出力することもできる。さらに、出力部１０９０は、第１及び第２の変位関連情報に基づき生成した多次元の変位関連情報を出力することもできる。変位関連情報は、数値データの形式で出力することもできるし、画像情報の形式で出力することも可能である。また、出力部１０９０は、変位関連情報を、受信データから生成された超音波断層像（Ｂモード像など）に合成（オーバーラップ）表示することができる。たとえば、超音波断層像中の測定対象物１０８９の該当位置に数値データを表示したり、変位、歪み、弾性係数などの数値の大きさを擬似色で表した超音波断層像を表示したりすることができる。このような出力結果により、人体１０９９の内部にある測定対象物１０８９（腫瘍など）の歪特性や弾性特性を容易に検査でき、診断に役立てることができる。

第１の超音波プローブ１０１０の送受信と第２の超音波プローブ１０２０の送受信が異なるタイミングで行われる場合、厳密には異なる時刻に受信した第１及び第２の受信データを用いざるを得ない。このような場合に、本実施形態では、増加期間に得られた受信デ
ータと減少期間に得られた受信データとを混在して用いることを禁止し、必ず同じ（もしくは同じ種類の）期間において得られた受信データから変位関連情報を演算する。これにより、測定対象物１０８９の変形状態や挙動が共通する受信データから各方向の変位を算出することができ、極めて有用性の高い、高精度な歪情報及び弾性パラメータ情報が得られる。

なお、不図示の治具で予め第１及び第２の超音波プローブが固定され、両者の位置関係（プローブの相対位置と向き）が分かっている場合には、特定の二方向に関する歪み情報が得られることになる。勿論、超音波プローブが治具により固定されていない場合であっても、例えばそれぞれの超音波プローブに、その位置や向きに関する情報を出力するセンサを備えさせておけば、当該センサからの情報を用いて、所定方向に関する歪み情報に変換することができる。詳細は後述する。

なお、増加期間の受信データのみを用いる場合に、各プローブからの超音波の送信は以下のように行うことが可能である。ｉ）増加期間か減少期間かにかかわらず、常時、超音波信号の送信及び受信を行う。そして、演算処理部１０５０が増加期間の受信データのみを処理の対象として抽出する。ｉｉ）増加期間か減少期間かにかかわらず、常時、超音波信号は送信する。ただし、受信は、前記第１及び第２の超音波プローブ共に、増加期間でしか行わない。ｉｉｉ）増加期間においてしか、第１及び第２の超音波プローブによる超音波信号の送信及び受信を行わない。なお、減少期間の受信データのみを用いる場合についても、同様に考えることができる。勿論、本発明は、これらｉ）、ｉｉ）、ｉｉｉ）に記載した超音波信号の送受制御に限定されるものではない。

（第１及び第２の超音波プローブ）
本発明においては、２以上の複数の超音波プローブを備えていれば良く、使用する周波数帯域としては、例えば１ＭＨｚから１５ＭＨｚである。

超音波プローブ同士が治具などで相対的に固定されている構成では、プローブ同士の相対的な位置関係情報（相対位置とプローブの向きに関する情報）は、予めわかっている。ただし、プローブの配置に関して高い自由度を確保するという点からは、プローブにセンサを設けておき、当該センサからの情報を用いて前述の位置関係情報を取得することが望ましい。ここでいうセンサには、６軸方向の情報を取得できる磁気センサや光学センサやエンコーダなどが含まれる。特に、センサの精度は１ｍｍ以下であることが好ましい。

また、超音波プローブ同士が治具により固定されていない場合は、両者の測定領域（図１の１０１５，１０２５）や音線方向が交差していることを前述のセンサを用いて確認したりすることは望ましい形態である。

超音波プローブ同士を治具により固定する場合には、複数のプローブは各々、ｎ回対称（ｎは２以上の整数）な位置に配置されているのがよい。非対称配置はプローブの位置精度に対する非対称性を招き、誤差論上不利だからである。

また、複数のプローブの相対位置は、例えば１０ｍｍ以上（４００ｍｍ以下）離れているように固定しておくことも望ましい形態である。相対位置が近すぎると精度上の問題が生じるためである。また複数のプローブを用いる効果が薄れるからである。

これらの超音波プローブを用いることにより、同定された方向に関する測定対象の変位情報、歪み情報の取得が可能となる。なお、測定対象に対して加える圧力の情報を用いれば、変位情報と圧力情報とから弾性係数に関する情報の取得も可能である。

複数の超音波プローブを使用する超音波測定装置においては、受信時の混信などを防ぐために、それぞれのプローブから出力される超音波の周波数を異ならせておくこともできる。受信信号から周波数分離により、いずれのプローブから出力された超音波により生じた信号であるかを同定できれば、誤った信号処理が回避できるからである。このように互いに異なる周波数で超音波の送受を行うことは、予め設定しておいてもよい。あるいはシステム制御部１０７０が、第１及び第２の超音波プローブが互いに異なる周波数の超音波を送受信するように、送信回路１０６０及び受信回路１０９０を制御することもできる。

また、受信時の混信は、プローブ同士の超音波送信領域が重複している場合に生じると考えられ、当該混信回避の為には、超音波プローブによる超音波の送受信を図１０に示すように、交互に行うことが考えられる。一方、それぞれのプローブで取得された受信信号を元に測定対象物の断層像を生成し表示している場合には、複数のプローブ間の信号の送受を図１０のように交互に行うと、必然的にフレームレートを落とさざるを得ず、滑らかな画像表示が難しい。そのため、複数の超音波プローブを使って多次元の歪みに関する情報を取得する時、即ち、プローブ同士の測定領域を重複させて信号の送受を行う際に、選択的に図１０のようにプローブ間で交互の信号送受を行う構成にすることも可能である。そして、多次元の歪みに関する情報を取得する必要が無い場合や、プローブ同士の超音波送信領域が重複していない場合には、フレームレートを高めるべく、一つのプローブのみで、フレームレートの高い画像情報を取得することができる。具体的なフローを図１１に示す。まず、第１の及び第２の超音波プローブの相対位置と向きに関する情報を取得する（Ｓ１）。そして、第１及び第２の超音波プローブの送信領域を入力し（Ｓ２）、互いの超音波送信領域が重複しているか否かを判断する（Ｓ３）。重複している場合には（Ｓ３；ＹＥＳ）、２本のプローブを送信タイミングをずらして用いる（Ｓ４）。一方、重複していない場合には（Ｓ３；ＮＯ）、２本のプローブの一方を用いて送受信する（Ｓ５）。各プローブによる超音波の送信タイミングは、システム制御部１０７０によって制御可能である。

演算処理部１０５０が、前記第１の方向と前記第２の方向の相対関係に基づいて、前記第１及び第２の変位関連情報から多次元の変位関連情報を生成し、出力部１０８０に出力することができる。前述のように第１及び第２の超音波プローブ１０１０、１０２０の相対的な位置と向きに関する情報を取得するセンサを備えている場合、演算処理部１０５０は、該センサから取得される情報を用いて、第１の方向と第２の方向との相対関係を取得できる。演算処理部１０５０は、前記第１及び第２の方向から構成される座標系（斜行座標系）から直交座標系へと、第１及び第２の変位関連情報を座標変換することにより、多次元の変位関連情報を生成することができる。なお、第１及び第２の超音波プローブが、後述する実施例のように、プローブの位置と向きに関する情報を出力するためのセンサを備えていることは好ましい形態である。

なお、第１の超音波プローブ１０１０の超音波送信領域１０１５と第２の超音波プローブ１０２０の超音波送信領域１０２５とが重複箇所を有しているか否かの情報を出力する手段（例えば、画面表示や通知音である）を有することも好ましい形態である。また、第１の超音波プローブ１０１０の超音波の音線１０１１と第２の超音波プローブ１０２０の超音波の音線１０２１とが、互いに交差しているか否かの情報を出力する手段（例えば、画面表示や通知音である）を有することも好ましい形態である。このような情報を出力することで、プローブの位置や向きが、多次元の変位関連情報を測定するのに適切か否かを操作者に知らせることができ、操作性の向上を図ることができる。なお、超音波送信領域が重複しているか否か、及び、音線が交差しているか否かは、演算処理部１０５０が、プローブの位置と向きに関する情報に基づいて判定することができる。

（情報取得部）
情報取得部１０４０は、圧力センサ、速度センサ、加速度センサ、あるいは位置センサで構成可能である。たとえばプローブ（あるいは圧迫手段）に圧力センサを設けることで、プローブ等が人体１０９９を押圧する力を検知することができるため、測定対象物１０８９に加わる力の時間変化を間接的に検知可能となる。同様に、プローブ等に速度センサ、加速度センサ、位置センサなどを設けることで、プローブ等の変位を検知でき、測定対象物１０８９の変形量の時間変化を間接的に検知することができる。

なお、情報取得部は、第１の受信データあるいは第２の受信データから、測定対象物の変形量の時間変化を算出することも可能である。すなわち、異なる時刻の受信データを解析することにより、測定対象物の変位量（変形量）を取得することができる。あるいは、受信データから生成された画像情報（エコー像）を解析してもよい。たとえば、時系列画像の差分をとることで、測定対象物の変位量（変形量）を求めることができる。このように受信データに基づき増加期間か減少期間かを判断する場合、情報取得部は、演算処理部１０５０の一機能として実現される。

（圧迫手段）
測定対象物１０８９に対して力を加える為の圧迫手段を、第１及び第２の超音波プローブ１０１０、１０２０とは別に設けることができる。勿論、圧迫手段として、第１及び第２の超音波プローブ１０１０、１０２０の一方あるいは両方を利用することもできる。

圧迫手段を前記超音波プローブとは別に設ける場合には、第１及び第２の超音波プローブ１０１０、１０２０の間に配置することができる。３つ以上の超音波プローブを用いる場合には、それらの超音波プローブの形成する多角形の内側に圧迫手段を配置することが好ましい形態である。斯かる配置は、圧迫による各プローブの測定変位量がほぼ等しくなり、測定精度を向上させる効果がある。

また、複数の超音波プローブが治具により相対的に固定されている場合には、前記圧迫手段は前記第１及び第２の超音波プローブの中間、または上記３つ以上の超音波プローブで構成される多角形の中心若しくは重心に配置することができる。勿論、圧迫手段は、複数具備されていてもよい。

次に、具体的な実施例を詳細に説明する。

（実施例１）
以下、実施例１では試料（測定対象物）を圧迫することによって試料の弾性特性を計測する方法について説明する。

図３は実施例１の超音波測定装置の概略図である。図３を使用し信号の流れを以下説明する。超音波測定装置は、複数の超音波プローブ１ａ、１ｂ、１ｃ、送信回路系６、受信回路系７、プローブ位置処理系８、断層信号処理系９、変位信号処理系１０、画像処理系１１、画像表示装置１２、システム制御部１３、圧迫データ処理部１７を備える。なお、図３における変位信号処理系１０が、上述の実施形態における演算処理部に該当する。

図３において、システム制御部１３によって送信回路系６に超音波送信の命令が転送され、送信回路系６から超音波プローブ１ａ、１ｂ、１ｃに超音波信号が送信される。超音波プローブ１ａ、１ｂ、１ｃは試料２に超音波を照射し、反射した超音波を受信する。この際、超音波プローブ１ａ、１ｂ、１ｃそれぞれの６軸位置は各超音波プローブ１ａ、１ｂ、１ｃに取り付けられた位置センサ５ａ、５ｂ、５ｃによって随時検出されている。受信された超音波信号及び超音波プローブ位置信号はそれぞれ受信回路系７、プローブ位置処理系８によって処理される。また、プローブ位置処理系８によって処理された各プロー
ブ位置は受信回路系７に入力され、超音波信号の補正に使用される。さらに、プローブ１ａには圧力センサ１６が取付けられており、プローブ１ａに加えられる力の状態をモニタできるようになっている。この圧力センサ１６と圧迫データ処理部１７が、上述した実施形態の情報取得部に該当する。

次に受信信号は断層信号処理系９及び変位信号処理系１０に入力され、断層信号生成処理、変位信号処理がなされる。断層信号処理系９は、得られた超音波信号を包絡線検波し、その強度を輝度信号として画像処理系１１に入力する。変位信号処理系１０は試料中の各点に対応する超音波信号の時間変動を変位信号として演算処理し、輝度信号もしくは色信号として画像処理系１１に入力する。これら処理結果は画像処理系１１に入力され画像表示装置１２よって断層信号及び変位信号が表示される。

次に試料の弾性特性を計測するプロセス及び変位信号処理系１０の処理プロセスについて説明する。図４は本実施例にかかる超音波測定装置のプローブ部分を拡大した図である。本実施例の超音波測定装置は、３つのプローブ１ａ，１ｂ，１ｃを有しており、これらプローブより試料２に超音波の照射を行っている。試料の弾性特性の計測は、試料に歪を生じさせることによって実現される。本実施例においては、プローブ１ａにより試料を圧迫し、圧迫前及び圧迫後の試料中各点の移動を測定することによって試料の変位や歪を演算する。さらに圧力センサ１６によって測定された圧力を圧迫データ処理部１７によって解析することによって、加えられた圧力が増加傾向にあるか減少傾向にあるかを判別している。

試料中各点の歪量は、３本の超音波プローブ１ａ，１ｂ，１ｃによって取得された試料からの反射超音波信号を処理することによって求めることができる。信号処理の概念を図５を用いて説明する。図５では圧迫動作前後の試料の動きと、各々の超音波信号を示している。試料中に２ａの如く周囲と音響インピーダンス、弾性パラメータの異なる物質が存在する場合、超音波プローブ１ａによって得られる信号は１５ａの如く物質界面にて強い反射信号を伴う。試料圧迫前と圧迫後に其々の超音波信号１５ａ、１５ｂを取得すると、界面間の距離Ｌ，Ｌ´を測定することができる。従って超音波プローブ１ａによって得られる歪量は（Ｌ−Ｌ´）のように求めることができる。同様にプローブ１ｂ、１ｃについても同様の演算によって各軸方向の歪を求めることができる。

実際の信号処理においては信号の位相変化もしくは得られた画像の処理によって歪量を演算するが、本実施例では信号処理方法としては自己相関法と呼ばれる方法を用いて説明する。空間相関法もしくはその他の演算処理を用いてもほぼ同様の結果を得ることができる。自己相関法は異なる二つの時間における超音波信号の位相差を計測するものであり、特開平５-１６１６４７号公報によれば下記の原理による。

散乱体から反射してくる超音波受信信号Ｓ_１（ｔ）は、次の（１）式で表される。ここに、Ａ_１（ｔ）は振幅、Φ_１（ｔ）は位相、ω_０は送信信号の角周波数を示す。

Ｓ_１（ｔ）＝Ａ_１（ｔ）ｅｘｐ〔−ｊ（ω_０ｔ＋Φ_１（ｔ））〕 …（１）

散乱媒質全体がΔｘだけ変位したときの受信信号をＳ_２（ｔ）とすると、これは次の（２）式で表される。

Ｓ_２（ｔ）＝Ｓ_１（ｔ−Δｔ）＝Ａ_１（ｔ−Δｔ）ｅｘｐ〔−ｊ（ω_０・（ｔ−Δｔ）
＋Φ_１（ｔ−Δｔ））〕 …（２）

ここで、Δｔは散乱媒質がΔｘだけ変位したときの伝搬時間差（変位時間）であり、次の（３）式で表される。ただし、ｃは音速である。

Δｔ＝２・Δｘ／ｃ …（３）

この受信信号Ｓ_１（ｔ）、Ｓ_２（ｔ）を角周波数ω_０で直交検波した出力をそれぞれＩ_１（ｔ）、Ｑ_１（ｔ）、Ｉ_２（ｔ）、Ｑ_２（ｔ）とすると、位相θ_１（ｔ）、θ_２（ｔ）は次の（４）、（５）式で与えられる。

θ_１（ｔ）＝Φ_１（ｔ）
＝tan^−１（Ｑ_１（ｔ）／Ｉ_１（ｔ）） …（４）
θ_２（ｔ）＝Φ_１（ｔ−Δｔ）＋ω_０Δｔ
＝tan^−１（Ｑ_２（ｔ）／Ｉ_２（ｔ）） …（５）

従って、変位前後の両者の位相差Δθは、次の（６）式で与えられる。

Δθ＝θ_２（ｔ）−θ_１（ｔ）
＝Φ_１（ｔ−Δｔ）−Φ_１（ｔ）＋ω_０Δｔ
＝tan^−１（Ｑ_１（ｔ）／Ｉ_１（ｔ））
−tan^−１（Ｑ_２（ｔ）／Ｉ_２（ｔ）） …（６）

ここで、Φ_１（ｔ−Δｔ）はΦ_１（ｔ）にほぼ等しいことから、（６）式は次の（７）式となる。

Δｘ＝（ｃ／２ω_０）Δθ
＝（ｃ／２ω_０）｛tan^−１（Ｑ_１（ｔ）／Ｉ_１（ｔ））−tan^−１（Ｑ_２（ｔ）／Ｉ_２（ｔ））｝ …（７）

従って、変位量Δｘは超音波信号を（７）式を用いた処理を行うことによって容易に求められる。

次に試料中各点の移動より歪量を求める方法について説明する。圧迫による試料中各点の移動は３本の超音波プローブ１ａ，１ｂ，１ｃによって計測される。図６は図４における、超音波プローブ面と各プローブの視野及びある計測点を模式的に示したものである。超音波プローブ１ａ，１ｂ，１ｃの各プローブ面を３ａ，３ｂ，３ｃ、一つの計測点を４ａとする。各超音波プローブ面から試料中の１点に向かう音線によって形成される座標系は斜交座標系であり、各プローブによって計測される点の三次元的な移動は上記（１）〜（７）式の演算を用いて斜交座標系上の射影として記録される。ここで、各プローブの相対位置が６軸センサ５ａ，５ｂ，５ｃによって随時計測されている為、座標変換によって直交座標系に変換することができる。また、変換の容易さ及び使い勝手を考慮すると斜交座標系、直交座標系の一つの軸が圧迫手段の圧迫方向、もしくはいずれかのプローブの中心軸に一致することが望ましい。

従って本実施例によれば、試料中各点の移動を３次元的に計測することができ、それに
伴い歪量Ｓも以下の式にて求めることができる。

ｕ_ｘ，ｕ_ｙ，ｕ_ｚは試料中各点の変位量を示す。つまり式（８）の処理を試料中の必要な領域に対して行うことによって、領域の歪量分布を精確に知ることが可能となる。

試料の弾性係数分布は下式の通り試料内部の歪分布と応力分布から求められる。

Ｔ＝Ｃ：Ｓ …（９）

しかしながら試料内部の応力分布Ｔを直接計測することは現状では困難である為、歪分布と試料圧縮の際の境界条件から逆問題的に弾性係数分布を計算することとなる。また、幾つかの仮定を置くことにより弾性係数分布を計算する方法も試みられている。以下に弾性係数分布を求める手法の一例を示す。

応力分布を求める際に無限均質媒体の無限一様応力場を仮定した場合には、以下の式が成り立つ。

Ｔ_１は超音波プローブに圧力センサを付加することにより測定可能である。従って弾性係数Ｃも求めることが可能である。また応力場が有限長である場合においても均質媒体である場合には応力分布式を理論的に計算することができる。つまり式（１０）において０としたＴ_２〜Ｔ_６はＴ_１の関数となる為、同様に応力場を計算することが可能である。

以上の計算によって弾性係数Ｃを求めることができた為、ヤング率Ｅ、ポアソン比σなどの弾性パラメータは以下の式を用いて導くことが可能となる。

以上の方式に限らず、歪から弾性係数を導出する方法は試料が非圧縮であることを仮定する方法、有限要素法を用いた方法等考案されているが、何れの方法を用いても同様の演算によって弾性パラメータを導出できる。

次に超音波信号の取得タイミングについて詳述する。超音波データの取得は圧迫動作中に取得される。取得可能な超音波データには、プローブ１ａ，１ｂ，１ｃの３種類及び圧迫前後の時系列データが存在する。本実施例では、データの取得は圧力の増加期間、もしくは減少期間のいずれかの期間でなされる。この状態を説明したものが図９である。図９では横軸に時間、縦軸に圧迫圧力を示している。図示するように、プローブ１ａ，１ｂ，１ｃの超音波データが交互に取得されていく。本方式の原理上、圧迫方向が反転する時点（つまり増加期間と減少期間が切り替わる時点）の前後のデータでは歪量の差が非常に小さく、有意なデータを取得することができない。よって、圧迫データ処理部１７によって圧迫方向の判別を行い、圧迫方向が反転する際にはデータ取得を行わない。これにより常に良好な歪量計測が可能となる。

以上述べたように本実施例によれば、超音波信号から試料の変位、歪、弾性パラメータを導出することが可能となる。

上記方法を用いて歪及び弾性分布を求める場合には各超音波プローブの配置によって得られる精度が異なることに注意する必要がある。これは上に述べたように各超音波プローブの音線が形成される座標系は斜交座標系であるためであり、各超音波プローブの間隔が狭い場合大きな誤差を伴う。医用超音波診断装置であることを考慮した場合には、試料は人体であり観察深さは数ｍｍ〜数十ｃｍ程度であるため、最大深さを２０ｃｍ、精度の悪化を直交座標の５％までとすると、プローブ間隔は約１０ｍｍ以上離れていることが好ましい。また、センサ５ａ，５ｂ，５ｃに要求される精度も同様に測定点の測定精度より小さいことが必要である。医用超音波装置の解像度が１ｍｍ以下であることから、センサ精度も１ｍｍ以下であることが望ましい。さらに各超音波プローブの計測誤差を最も小さくするには超音波プローブは３回対称であることが好ましい。

以上の方法による歪及び弾性分布を求める方法では圧迫手段は一つのプローブであるが、二本、もしくは三本のプローブを同時もしくは交互に圧迫し、各々の圧迫によって得られた歪量から演算によって歪量を推定しても良い。平均値、最大、最小値、中央値等、対象に応じた演算結果を出力することが可能である。

また、本実施例においては、増加期間か減少期間かを判定する目的にて圧力センサを用いている。しかし、圧力センサに限らず速度センサ、加速度センサ、位置センサ等のセンサ、プローブ１ａ，１ｂ，１ｃによる超音波信号の解析、ディスプレイ出力の解析等の方法によって圧迫方向を演算することによっても同様の効果を得ることができる。あるいは、プローブ間の位置及び向きを計測する目的にて配置されている位置センサ５ａ，５ｂ，５ｃの出力を解析して、増加期間か減少期間かを判定することもできる。

なお、本実施例においては超音波プローブが３本である場合について説明したが、４本またそれ以上の複数のプローブを用いた場合においても同様の類比演算によって同様の効果を得ることができる。プローブが２本である場合についても精度は落ちるものの等方性などの仮定を入れることによって同様の効果を得ることができる。また、各精度及び対称性に関しては本実施例から外れた場合には精度が落ちることとなるが同様の効果を得ることが出来る。

（実施例２）
実施例２においては実施例１と異なる方法にて試料の歪を計測する方法を説明する。実施例１では試料を圧迫することによって試料の歪計測を行ったが、試料が自ら動く場合には測定点の移動を伴うため圧迫動作を行うことなく歪計測を行うことができる。たとえば、心臓や心臓周辺の部位は、心筋運動によって周期的に変形しているため、外部からの圧迫動作は不要である。このような場合、図３において超音波プローブ１ａ，１ｂ，１ｃを動かすことなく歪計測を行う。また計測方法及び演算も実施例１と全く同じ方式であり、同様の効果を得ることが出来る。

（実施例３）
実施例３においては実施例１と異なる圧迫手段を用いて弾性係数分布計測を行った例について説明する。図７に本実施例を用いた試料周辺の図を示す。実施例１と同様に超音波プローブ１ａ，１ｂ，１ｃ、位置センサ５ａ，５ｂ，５ｃ及び試料２が配置されているが、本実施例においてはプローブ以外に圧迫手段１４が配置されている。超音波装置の全体システムは実施例１（図３）と同様である。測定時は、圧迫手段１４で試料を圧迫し、圧迫前及び圧迫後の超音波信号を超音波プローブ１ａ，１ｂ，１ｃによって取得する。また、各超音波プローブ１ａ，１ｂ，１ｃ及び圧迫手段１４は対称性を持つことが最も測定精度が高くなる。よって、図８の上面図に示す如く、超音波プローブ１ａ，１ｂ，１ｃを圧迫手段１４を対称軸とする３回対称性を持つ配置とするのがよい。本実施例によっても実施例１と同様の歪量及び弾性量を取得することが可能となる。

さらに実施例１と同様に十分に精度の高い計測を行うにはプローブ間隔は約１０ｍｍ以上離れていることが必要であり、センサの精度は１ｍｍ以下が必要である。

以上の方法による歪及び弾性分布を求める方法では圧迫手段は一つの圧迫手段であるが、二本、もしくは三本以上の圧迫手段もしくはプローブを同時もしくは交互に圧迫し、各々の圧迫によって得られた歪量から演算によって歪量を推定しても良い。平均値、最大、最小値、中央値等、対象に応じた演算結果を出力することが可能である。

なお、本実施例においては超音波プローブが３本である場合について説明したが、２本、４本またそれ以上の複数のプローブを用いた場合においても同様の類比演算によって同様の効果を得ることが出来る。また、各精度及び対称性に関しては本実施例から外れた場合には精度が落ちることとなるが同様の効果を得ることが出来る。

１ａ、１ｂ、１ｃ 超音波プローブ
２、２ａ 試料
３ａ、３ｂ、３ｃ 超音波プローブ面
４、４ａ 超音波プローブの作る視野と測定点
５ａ、５ｂ、５ｃ 位置センサ
６ 送信回路系
７ 受信回路系
８ プローブ位置処理系
９ 断層信号処理系
１０ 変位信号処理系
１１ 画像処理系
１２ 画像表示装置
１３ システム制御部
１４ 圧迫手段
１５ａ、１５ｂ 超音波信号
１６ 圧力センサ
１７ 圧迫データ処理部

Claims (12)

1. 超音波測定装置であって、
   第１の超音波プローブと第２の超音波プローブを少なくとも含む複数の超音波プローブと、
   測定対象物に加わる力もしくは測定対象物の変形量の時間変化が、増加傾向を示す増加期間であるか、または前記時間変化が減少傾向を示す減少期間であるかに関する情報を取得する情報取得部と、
   演算処理部と、
   を備え、
   前記第１の超音波プローブは、該第１の超音波プローブの超音波送信領域と前記第２の超音波プローブの超音波送信領域との重複箇所に位置する前記測定対象物によって反射された超音波を受信して、第１の受信データを生成し、
   前記第２の超音波プローブは、前記重複箇所に位置する前記測定対象物によって反射された超音波を受信して、第２の受信データを生成し、
   前記演算処理部は、
   前記第１の受信データの内、前記増加期間あるいは前記減少期間のいずれか一方の期間におけるデータを用いて、前記測定対象物の第１の方向に関する第１の変位関連情報を演算し、且つ
   前記第２の受信データの内、前記第１の変位関連情報を演算するために用いた期間と同じ期間におけるデータを用いて、前記測定対象物の第２の方向に関する第２の変位関連情報を演算し、
   前記第１の変位関連情報と、前記第２の変位関連情報と、
   前記第１及び第２の超音波プローブの相対的な位置関係情報と、を用いて多次元の変位関連情報を生成することを特徴とする超音波測定装置。
2. 前記情報取得部は、圧力センサ、速度センサ、加速度センサ、あるいは位置センサによって、前記測定対象物に加わる力もしくは前記測定対象物の変形量の時間変化を検知することを特徴とする請求項１に記載の超音波測定装置。
3. 前記測定対象物に対して力を加える圧迫手段を有する
   請求項１または２に記載の超音波測定装置。
4. 前記圧迫手段は、前記第１及び第２の超音波プローブの間、若しくは前記複数の超音波プローブの形成する多角形の内側に位置する
   ことを特徴とする請求項３に記載の超音波測定装置。
5. 前記複数の超音波プローブは治具により相対的に固定されていることを特徴とする請求項４に記載の超音波測定装置。
6. 前記変位関連情報とは、前記測定対象物の変位、歪みあるいは弾性係数に関する情報のことである
   請求項１から５のいずれか１項に記載の超音波測定装置。
7. 前記第１の超音波プローブの超音波送信領域と前記第２の超音波プローブの超音波送信領域とが前記重複箇所を有しているか否かの情報を出力する手段を有する
   ことを特徴とする請求項１から６のいずれか１項に記載の超音波測定装置。
8. 前記第１の超音波プローブの超音波の音線と前記第２の超音波プローブの超音波の音線とが、互いに交差しているか否かの情報を出力する手段を有する
   ことを特徴とする請求項１から６のいずれか１項に記載の超音波測定装置。
9. 前記第１及び第２の超音波プローブに、超音波を送信するための信号を入力する送信回路と、前記第１及び第２の超音波プローブが超音波を受信することにより生成する信号を受信する受信回路と、前記送信回路と前記受信回路とを制御するシステム制御部と、をさらに備え、
   当該システム制御部は、前記第１及び第２の超音波プローブが交互に超音波を送受信するように、前記送信回路及び前記受信回路を制御する
   ことを特徴とする請求項１から８のいずれか１項に記載の超音波測定装置。
10. 前記第１及び第２の超音波プローブに、超音波を送信するための信号を入力する送信回路と、前記第１及び第２の超音波プローブが超音波を受信することにより生成する信号を受信する受信回路と、前記送信回路と前記受信回路と、を制御するシステム制御部と、をさらに備え、
    当該システム制御部は、前記第１及び第２の超音波プローブが互いに異なる周波数の超音波を送受信するように、前記送信回路及び前記受信回路を制御する
    ことを特徴とする請求項１から８のいずれか１項に記載の超音波測定装置。
11. 前記演算処理部は、前記第１及び第２の方向から構成される座標系から直交座標系へと、前記第１の受信データを用いて求められる前記第１の方向に関する第１の変位関連情報と、第２の受信データを用いて求められる前記第２の方向に関する第２の変位関連情報と、を座標変換することにより、前記多次元の変位関連情報を生成する
    ことを特徴とする請求項１から１０のいずれか１項に記載の超音波測定装置。
12. 前記第１及び第２の超音波プローブは離れて固定されていることを特徴とする請求項１から１０のいずれか１項に記載の超音波測定装置。

Images