(57)【特許請求の範囲】
【請求項１】
変位用超音波ビームを形成して関心部位において組織を変位させる変位波処理部と、
測定用超音波ビームを形成して前記関心部位から受信信号を得る測定波処理部と、
前記変位用超音波ビームに対する変調処理を制御する変調制御部と、
前記測定用超音波ビームを介して得られる受信信号に基づいて、前記関心部位における組織の変位を測定する変位測定部と、
前記変位の測定結果に基づいて、前記関心部位における組織の凝固を測定する凝固測定部と、
を有し、
前記変調制御部は、比較的高い変調周波数と比較的低い変調周波数を用いて変位用超音波ビームを変調処理するように前記変位波処理部を制御し、
前記変位測定部は、前記各変調周波数ごとに関心部位における組織の変位を測定し、
前記凝固測定部は、前記比較的高い変調周波数による変位の測定結果に基づいて、関心部位における局所的な凝固を測定し、前記比較的低い変調周波数による変位の測定結果に基づいて、関心部位における広域的な凝固を測定する、
ことを特徴とする超音波医用装置。
【請求項２】
請求項１に記載の超音波医用装置において、
前記凝固測定部は、前記各変調周波数ごとに得られる変位の測定結果に基づいて、関心部位における凝固のサイズを測定する、
ことを特徴とする超音波医用装置。
【請求項３】
請求項２に記載の超音波医用装置において、
治療用超音波ビームを形成して前記関心部位の組織を加熱して治療する治療波処理部をさらに有し、
前記凝固測定部は、前記加熱の期間内において複数の時相に亘って各時相ごとに関心部位における凝固のサイズを測定するにあたり、前記比較的高い変調周波数による変位の測定結果に基づいて、凝固発生初期の時相における局所的な凝固のサイズを測定し、前記比較的低い変調周波数による変位の測定結果に基づいて、凝固進行後の時相における広域的な凝固のサイズを測定する、
ことを特徴とする超音波医用装置。
【請求項４】
請求項３に記載の超音波医用装置において、
複数の時相を一方の軸に示し、各時相ごとに測定された凝固のサイズを他方の軸に示した凝固状態画像を形成する画像形成部をさらに有する、
ことを特徴とする超音波医用装置。
【請求項５】
請求項１から４のいずれか１項に記載の超音波医用装置において、
前記変位波処理部は、前記比較的高い変調周波数による変調処理と前記比較的低い変調周波数による変調処理とを合成した変位用超音波ビームを形成し、
前記変位測定部は、前記測定用超音波ビームを介して得られる受信信号から、前記各変調周波数に対応した周波数成分を抽出することにより、前記各変調周波数ごとに関心部位における組織の変位を測定する、
ことを特徴とする超音波医用装置。
【請求項６】
請求項１から４のいずれか１項に記載の超音波医用装置において、
前記変位波処理部は、前記比較的高い変調周波数により変調処理された変位用超音波ビームと、前記比較的低い変調周波数により変調処理された変位用超音波ビームと、を互いに異なる時相において形成し、
前記測定波処理部は、前記各変調周波数ごとにその変調周波数に対応した時相において測定用超音波ビームを形成し、
前記変位測定部は、前記各変調周波数ごとに形成される測定用超音波ビームを介して得られる受信信号に基づいて、前記各変調周波数ごとに関心部位における組織の変位を測定する、
ことを特徴とする超音波医用装置。
【請求項７】
請求項１から４のいずれか１項に記載の超音波医用装置において、
前記変調制御部は、前記比較的高い変調周波数による変位の測定結果に基づいて測定される凝固のサイズが閾値に達した場合に、前記比較的高い変調周波数から前記比較的低い変調周波数に切り替えるように、前記変位波処理部を制御する、
ことを特徴とする超音波医用装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、組織の凝固を測定する超音波医用装置に関する。
【背景技術】
【０００２】
強力集束超音波（ＨＩＦＵ：High Intensity Focused Ultrasound）を例えば生体に照射し、その音響エネルギーを利用して腫瘍などの治療部位を加熱して凝固させる治療方法が知られている。
【０００３】
組織が加熱されて凝固する際に、その凝固の前後において、組織の弾性率（ヤング率）が増大することが知られている。また、ＨＩＦＵなどの比較的強力な超音波は、その進行方向に向かって放射力を発生させるため、例えば、ＨＩＦＵの超音波ビームの焦点部位の組織に１０〜１００μｍ（マイクロメートル）程度の変位を与えることができる。
【０００４】
そのため、ＨＩＦＵなどの比較的強力な超音波により組織に変位を与えて、弾性率の増大による変位の減少を測定して、組織の凝固を観察することができる。例えば、ＨＩＦＵの超音波を変調周波数ｆＭで振幅変調し、放射力の強度を変動させることにより、焦点部位の組織に振動を励起させ、その振動の変位または速度を超音波診断装置によって測定する。
【０００５】
この原理を用いて組織の凝固を検出し、その検出結果を画像上にマッピングする手法がＨＭＩ（Harmonic Motion Imaging）である（特許文献１，２参照）。放射力や組織の振動の周波数が上記変調周波数ｆＭの２倍になることからＨＭＩと呼ばれている。
【先行技術文献】
【特許文献】
【０００６】
【特許文献１】米国特許出願公開第２００５／０００４４６６号明細書
【特許文献２】米国特許出願公開第２００７／０２７６２４２号明細書
【発明の概要】
【発明が解決しようとする課題】
【０００７】
ところで、ＨＩＦＵを利用した治療では、治療部位における凝固の状況に応じてＨＩＦＵを適切に制御することが望ましい。これは生体へのＨＩＦＵの照射においては、その伝搬経路の音響特性が患者個々人によって異なることに起因して減衰や位相歪が焦点音圧に与える影響が異なることや、焦点近傍での血流の違いによる冷却効果の違いに起因する凝固に必要なエネルギーの変化などのために、治療のためのパラメータを事前に最適化することが難しいためである。例えば、治療部位において凝固が開始されるタイミングや、目的とする大きさまで凝固が終了するタイミングなどを知ることができれば、それらのタイミングなどに応じてＨＩＦＵの照射を制御することが可能になる。
【０００８】
このような状況において、本願の発明者らは、超音波を利用して組織の凝固を測定する技術について研究開発を重ねてきた。
【０００９】
本発明は、その研究開発の過程において成されたものであり、その目的は、超音波を利用した組織の凝固の測定において測定精度を高めることにある。
【課題を解決するための手段】
【００１０】
上記目的にかなう好適な超音波医用装置は、変位用超音波ビームを形成して関心部位において組織を変位させる変位波処理部と、測定用超音波ビームを形成して前記関心部位から受信信号を得る測定波処理部と、前記変位用超音波ビームに対する変調処理を制御する変調制御部と、前記測定用超音波ビームを介して得られる受信信号に基づいて、前記関心部位における組織の変位を測定する変位測定部と、前記変位の測定結果に基づいて、前記関心部位における組織の凝固を測定する凝固測定部と、を有し、前記変調制御部は、比較的高い変調周波数と比較的低い変調周波数を用いて変位用超音波ビームを変調処理するように前記変位波処理部を制御し、前記変位測定部は、前記各変調周波数ごとに関心部位における組織の変位を測定し、前記凝固測定部は、前記比較的高い変調周波数による変位の測定結果に基づいて、関心部位における局所的な凝固を測定し、前記比較的低い変調周波数による変位の測定結果に基づいて、関心部位における広域的な凝固を測定する、ことを特徴とする。
【００１１】
上記構成において、測定用超音波ビームは、例えば一般的な超音波診断装置における診断用の超音波ビームであり、診断用の超音波振動子を利用して形成することができる。これに対し、変位用超音波ビームは、放射力により組織を変位させる程度の比較的強度の大きい超音波により形成される。変位用超音波ビームは、診断用の超音波ビームに比べて強度が大きく、例えば強力集束超音波（ＨＩＦＵ：High Intensity Focused Ultrasound）により形成されてもよい。さらに、その強力集束超音波（ＨＩＦＵ）により組織を加熱して凝固させてもよい。この場合、例えば、加熱による治療の対象となる治療部位が関心部位となる。
【００１２】
そして、上記構成によれば、比較的高い変調周波数により、比較的狭い領域において限定的に組織を変位させることができ、変位の領域が狭いほど、より微小な（局所的な）凝固を検出することができる。そのため、比較的高い変調周波数により、例えば発生直後の局所的な凝固の有無や凝固のタイミングなどを高精度に測定できる。また、比較的低い変調周波数により、比較的広い領域において組織を変位させることができ、変位の領域が広いほど、より大きな（広域的な）凝固を検出することができる。そのため、比較的低い変調周波数により、例えば進行後の広域的な凝固の大きさや加熱による治療の終了タイミングなどを高精度に測定できる。
【００１３】
なお、比較的高い変調周波数と比較的低い変調周波数は、２つの変調周波数に限定されない。例えば、３つ以上の互いに異なる変調周波数が利用され、最も高い変調周波数により最も局所的な凝固が測定され、変調周波数が低くなるにつれて徐々に広い（広域的な）凝固が測定されるようにしてもよい。また、凝固の測定においては、凝固の有無や凝固の大きさ（サイズ）の他に凝固の程度（組織の歪や硬さ）などが測定されてもよい。
【００１４】
望ましい具体例において、前記凝固測定部は、前記各変調周波数ごとに得られる変位の測定結果に基づいて、関心部位における凝固のサイズを測定する、ことを特徴とする。
【００１５】
望ましい具体例において、治療用超音波ビームを形成して前記関心部位の組織を加熱して治療する治療波処理部をさらに有し、前記凝固測定部は、前記加熱の期間内において複数の時相に亘って各時相ごとに関心部位における凝固のサイズを測定するにあたり、前記比較的高い変調周波数による変位の測定結果に基づいて、凝固発生初期の時相における局所的な凝固のサイズを測定し、前記比較的低い変調周波数による変位の測定結果に基づいて、凝固進行後の時相における広域的な凝固のサイズを測定する、ことを特徴とする。
【００１６】
望ましい具体例において、複数の時相を一方の軸に示し、各時相ごとに測定された凝固のサイズを他方の軸に示した凝固状態画像を形成する画像形成部をさらに有する、ことを特徴とする。
【００１７】
望ましい具体例において、前記変位波処理部は、前記比較的高い変調周波数による変調処理と前記比較的低い変調周波数による変調処理とを合成した変位用超音波ビームを形成し、前記変位測定部は、前記測定用超音波ビームを介して得られる受信信号から、前記各変調周波数に対応した周波数成分を抽出することにより、前記各変調周波数ごとに関心部位における組織の変位を測定する、ことを特徴とする。
【００１８】
望ましい具体例において、前記変位波処理部は、前記比較的高い変調周波数により変調処理された変位用超音波ビームと、前記比較的低い変調周波数により変調処理された変位用超音波ビームと、を互いに異なる時相において形成し、前記測定波処理部は、前記各変調周波数ごとにその変調周波数に対応した時相において測定用超音波ビームを形成し、前記変位測定部は、前記各変調周波数ごとに形成される測定用超音波ビームを介して得られる受信信号に基づいて、前記各変調周波数ごとに関心部位における組織の変位を測定することを特徴とする。
【００１９】
望ましい具体例において、前記変調制御部は、前記比較的高い変調周波数による変位の測定結果に基づいて測定される凝固のサイズが閾値に達した場合に、前記比較的高い変調周波数から前記比較的低い変調周波数に切り替えるように、前記変位波処理部を制御することを特徴とする。
【発明の効果】
【００２０】
本発明により、超音波を利用した組織の凝固の測定において測定精度が高められる。例えば、本発明の好適な態様によれば、比較的高い変調周波数により、例えば発生直後の局所的な凝固の有無などを高精度に測定でき、比較的低い変調周波数により、例えば進行後の広域的な凝固の大きさなどを高精度に測定できる。
【００２１】
また、例えば、本発明の好適な態様によれば、凝固の開始時間が分かるようになることで、体表から焦点までの減衰量や、伝搬経路上の音響特性の不均一の効果など、外部から入力したエネルギーが焦点での音圧ピークを形成するときに患者毎に異なる非理想的な効果を補正することや、焦点での超音波吸収量や熱的な特性や血流量の個体差など焦点で一定の音圧ピークが形成されたとしても温度上昇の値が異なる効果を補正するためのデータとして使うことが可能となる。
【図面の簡単な説明】
【００２２】
【図１】本発明の実施において好適な超音波医用装置の全体構成を示す図である。
【図２】図１の超音波医用装置の動作を示すタイミングチャートである。
【図３】変調処理された変位発生用の超音波による組織の振動を説明するための図である。
【図４】変調周波数と凝固の関係を説明するための図である。
【図５】各変調周波数ごとに検出された凝固サイズの実験結果を示す図である。
【図６】変調周波数の設定の具体例を説明するための図である。
【図７】位相１から位相１３までのデータが収集される様子を示す図である。
【図８】複数の変調周波数の合成波を利用した具体例１のフローチャートである。
【図９】ＤＣ成分を持たない変調法とＤＣ成分を持つ変調法を示す図である。
【図１０】ＤＣ成分を持たない変調におけるＮＨとＮＬの対応関係を示す図である。
【図１１】ＮＬ＝１，ＮＨ＝５の場合の二変調周波数を示す図である。
【図１２】ＮＬ＝１，ＮＨ＝４の場合の二変調周波数を示す図である。
【図１３】ＤＣ成分を持つ変調におけるＮＨとＮＬの対応関係を示す図である。
【図１４】複数の変調周波数を切り替える具体例２を示すフローチャートである。
【図１５】複数の変調周波数を段階的に変更する具体例３のフローチャートである。
【図１６】凝固の判定結果に基づいて変調周波数を切り替える具体例４を示すフローチャートである。
【図１７】凝固サイズと変調周波数の対応関係を示す図である。
【図１８】図１の超音波医用装置により形成される凝固状態画像の具体例を示す図である。
【発明を実施するための形態】
【００２３】
図１は、本発明の実施において好適な超音波医用装置（本超音波医用装置）の全体構成図である。本超音波医用装置は、超音波プローブ１０を有しており、その超音波プローブ１０は、ＨＩＦＵ用振動子１０Ｈと診断用振動子１０Ｄを備えている。
【００２４】
ＨＩＦＵ用振動子１０Ｈは、強力集束超音波（ＨＩＦＵ）を送波する振動子であり、例えば二次元的に配列された複数の振動素子を備えている。ＨＩＦＵ用振動子１０Ｈは、例えば癌や腫瘍などの治療部位Ｐに向けて治療用超音波ビームＴＢを形成して強力集束超音波を送波し、その治療部位Ｐを加熱して治療するために利用される。
【００２５】
また、ＨＩＦＵ用振動子１０Ｈは、治療部位Ｐに向けて変位用超音波ビームＥＢを形成して変位発生用の超音波を送波し、治療部位Ｐにおいて放射力を発生させ組織を変位させる。変位用超音波ビームＥＢは、治療部位Ｐにおいて有効な放射力を発生させる程度の強度で形成されるビームであり、例えば、治療用超音波ビームＴＢを変位用超音波ビームＥＢとして利用してもよい。
【００２６】
一方、診断用振動子１０Ｄは、例えば二次元的に配列された複数の振動素子を備えており、例えば治療部位Ｐを有する被検体（患者）に対して、超音波画像を形成するための比較的弱い超音波を送受する。つまり、公知の一般的な超音波診断装置と同じ程度の強度（エネルギー）の超音波を送受する。
【００２７】
また、診断用振動子１０Ｄは、治療部位Ｐに向けて測定用超音波ビームＭＢを形成して測定用の超音波を送受し、測定用超音波ビームＭＢに沿って受信信号を得る。測定用超音波ビームＭＢに沿って得られる受信信号は、変位用超音波ビームＥＢの放射力による治療部位Ｐにおける変位の測定に利用される。
【００２８】
なお、超音波プローブ１０は、例えば、お椀（どんぶり）状に凹ませた内部の表面を振動子面とする。そして、例えば、お椀状に凹んだ内部の中央に位置する底の部分に診断用振動子１０Ｄが設けられ、診断用振動子１０Ｄを取り囲むようにＨＩＦＵ用振動子１０Ｈが設けられる。なお、超音波プローブ１０の振動子面の形状は、お椀状に限定されず、例えば治療の用途等に応じた形状とされることが望ましい。また、全ての振動素子またはいくつかの振動素子が、ＨＩＦＵ用と診断用の両用途に併用されてもよい。
【００２９】
測定診断ブロック２０は、診断用振動子１０Ｄの送受信を制御する送受信部２２を備えている。送受信部２２は、診断用振動子１０Ｄを構成する複数の振動素子の各々に対応した送信信号を出力することにより、診断用振動子１０Ｄを制御して送信ビームを形成し、さらに、それら複数の振動素子の各々から得られる受信信号に対して整相加算処理などを施すことにより、受信ビームに沿って受信信号を得る。
【００３０】
送受信部２２は、治療部位Ｐを含んだ三次元空間内または断面内で診断用の超音波ビームを走査させて画像用の受信信号を収集する。そして、収集された受信信号に基づいて、超音波画像形成部２８が三次元の超音波画像または二次元の断層画像の画像データを形成し、その画像データに対応した超音波画像が表示部５０に表示される。
【００３１】
ユーザ（検査者）は、表示部５０に表示される超音波画像から、治療部位Ｐの位置等を確認し、図示しない操作デバイス等を利用して治療部位Ｐの位置情報を本超音波医用装置に入力する。もちろん、本超音波医用装置が、超音波画像に対する画像解析処理等により治療部位Ｐの位置を確認して位置情報を得るようにしてもよい。
【００３２】
また、送受信部２２は、診断用振動子１０Ｄを制御して測定用超音波ビームＭＢを形成し、測定用超音波ビームＭＢに沿って受信信号を得る。そして、変位測定部２４は、測定用超音波ビームＭＢに沿って得られる受信信号に基づいて、治療部位Ｐにおける変位を測定する。また、凝固測定部２５は、治療部位Ｐにおける変位の測定結果に基づいて、治療部位Ｐにおける組織の凝固を測定する。さらに、凝固画像形成部２６は、治療部位Ｐにおける凝固の測定結果に基づいて凝固状態画像を形成し、その凝固状態画像が表示部５０に表示される。変位測定部２４と凝固測定部２５と凝固画像形成部２６における処理については後に詳述する。
【００３３】
一方、治療放射ブロック３０は、治療用送信部３２を備えており、治療用送信部３２はＨＩＦＵ用振動子１０Ｈを構成する複数の振動素子の各々に対応した送信信号を出力することにより、ＨＩＦＵ用振動子１０Ｈを制御して治療用超音波ビームＴＢを形成する。治療用送信部３２は、制御部４０によって制御され、例えば治療部位Ｐ内に焦点を設定した治療用超音波ビームＴＢが形成される。
【００３４】
また、治療放射ブロック３０は、変位用送信部３４を備えており、変位用送信部３４はＨＩＦＵ用振動子１０Ｈを構成する複数の振動素子の各々に対応した送信信号を出力することにより、ＨＩＦＵ用振動子１０Ｈを制御して変位用超音波ビームＥＢを形成する。変位用超音波ビームＥＢは変調処理を施され、その変調処理における変調周波数は、変調周波数制御部３６によって制御される。なお、変調周波数制御部３６は、制御部４０により制御される。
【００３５】
治療用超音波ビームＴＢに沿って強力集束超音波（ＨＩＦＵ）が送波されて治療部位Ｐが加熱されると治療部位Ｐの組織が凝固する。その凝固の前後において、組織の弾性率（ヤング率）が増大することが知られている。そして、組織の弾性率の変化を知るために、本超音波医用装置は、変位用超音波ビームＥＢに沿って超音波を送波して放射力を発生させ、その放射力による治療部位Ｐにおける組織の変位を測定する。その変位の測定は、測定用超音波ビームＭＢに沿って得られる受信信号に基づいて行われる。
【００３６】
なお、測定診断ブロック２０内の各部と治療放射ブロック３０内の各部は、それぞれ、例えばプロセッサや電子回路等のハードウェアを利用して実現することができる。制御部４０は、例えば、演算機能を備えたハードウェアとその動作を規定するソフトウェア（プログラム）によって構成される。表示部５０は、例えば液晶ディスプレイなどである。
【００３７】
また、測定診断ブロック２０は、一般的な超音波診断装置により実現されてもよい。そして、治療放射ブロック３０に対応した超音波治療装置と、測定診断ブロック２０に対応した超音波診断装置とを組み合わせたシステムにより、本超音波医用装置が具現化されてもよい。
【００３８】
図２は、図１の超音波医用装置（本超音波医用装置）の動作を示すタイミングチャートである。図１に示した部分（構成）については、以下の説明において図１の符号を参照する。
【００３９】
メイントリガは、強力集束超音波（ＨＩＦＵ）による治療の開始タイミングを示す信号であり、例えば、ユーザ（検査者）による治療開始の操作に応じて、制御部４０から本超音波医用装置内の各部へ出力される。
【００４０】
フレームトリガは、測定用超音波ビームＭＢのフレーム開始を示す信号である。送受信部２２は、例えばフレームトリガの立ち上がりのタイミングから、治療部位Ｐに向けて複数本の測定用超音波ビームＭＢを順に形成する。例えば、２つのフレームトリガの間において、治療部位Ｐに向けて１０本の送信ビームが形成され、１つの送信ビームにつき２本ずつ、２０本の受信ビームが形成される。もちろん、送信ビームと受信ビームの本数は上記の具体例に限定されない。
【００４１】
加熱期間信号は、治療用超音波ビームＴＢによる治療部位Ｐの加熱処理期間を示す信号であり、加熱期間信号の立ち上がりから立下りまでの期間において、例えば治療部位Ｐを焦点として治療用超音波ビームＴＢが形成される。
【００４２】
測定期間トリガは、変位測定の期間を示す信号であり、測定期間トリガの立ち上がりから立下りまでの期間が変位の測定期間となる。
【００４３】
変調信号は、治療用超音波ビームＴＢの変調処理に利用される変調信号であり、例えば変調周波数制御部３６から変位用送信部３４へ出力される。
【００４４】
ＨＩＦＵ信号は、治療用超音波ビームＴＢの送信信号であり、変位用送信部３４が、例えば周波数２ＭＨｚ程度の連続波を変調信号に従って振幅変調して得られる。
【００４５】
変調信号は、測定期間トリガの立ち上がりから立下りまでの測定期間において、振幅が０（ゼロ）とされ、その結果、測定期間においてＨＩＦＵ信号の振幅も０となり、測定期間において治療用超音波ビームＴＢの送信が停止される。
【００４６】
なお、測定期間トリガは、フレームトリガに対して遅延量（Ｄｅｌａｙ）を伴って出力される。この遅延量は、例えば、ユーザによって適宜に調整される。これにより、２つのフレームトリガの間において形成される複数本の測定用超音波ビームＭＢ（例えば２０本の受信ビーム）のうち、測定期間内に収まる１本又はいくつかの測定用超音波ビームＭＢが選択的に測定に利用される。
【００４７】
本超音波医用装置では、変位を発生させる変位用超音波ビームＥＢを変調周波数ｆＭで振幅変調し、放射力の強度を変動させることにより、焦点部位の組織に振動を励起させ、その振動の変位を測定用超音波ビームＭＢを利用して測定する。その際に、変調周波数制御部３６により、変位用超音波ビームＥＢの変調周波数ｆＭが制御される。なお、振幅変調に代えて、変調周波数ｆＭの周波数変調が利用されてもよい。
【００４８】
図３は、変調処理された変位発生用の超音波による組織の振動を説明するための図である。図３に矢印で示す進行方向（上下方向）に沿って変位用超音波ビームＥＢが形成されて超音波が照射されると、組織内において、変位用超音波ビームＥＢの中心から、左右両方向に向かって進行する、ずり波と呼ばれる横波が発生する。ずり波の振動の周波数は変位発生用超音波の変調周波数ｆＭの２倍となる。
【００４９】
図３には、比較的高い変調周波数の超音波によって発生するずり波と、比較的低い変調周波数の超音波によって発生するずり波の２つの波形が示されている。変調周波数が高いほど、組織の振動部位が局在するため、変位測定における位置分解能が高められ、微小な凝固域の検出に適する。しかし、変調周波数が高いと、変位の領域が狭いため、その変位の領域を超える広域的な凝固の検出には不向きであり、広域的な凝固の検出には低い変調周波数が望ましい。
【００５０】
図４は、変調周波数と凝固の関係を説明するための図である。図４には、超音波プローブ１０から放射される変位用超音波ビームＥＢと、組織内の凝固領域と振動領域が模式的に示されている。
【００５１】
＜Ａ＞は、変位用超音波ビームＥＢに対して比較的高い変調周波数（例えば２００Ｈｚ程度）の変調処理を施した場合を示している。比較的高い変調周波数においては、振動領域が比較的小さく局所的になり、その局所的な振動領域内における組織の平均弾性率の変化が大きく、小さな凝固領域の検出に適している。
【００５２】
一方、＜Ｂ＞は、変位用超音波ビームＥＢに対して比較的低い変調周波数（例えば３０Ｈｚ程度）の変調処理を施した場合を示している。比較的低い変調周波数においては、振動領域が比較的大きく広域的になり、広域的な振動領域内における組織の平均弾性率の変化を計測でき、大きな凝固領域の検出に適している。
【００５３】
図５は、各変調周波数ごとに検出された凝固サイズの実験結果を示す図である。図５には、変調周波数が３４Ｈｚ，６７Ｈｚ，１０２Ｈｚの３つの実験結果が示されており、各実験結果の横軸は測定部位を加熱した時間であり、縦軸は測定された凝固サイズを示している。また、各実験結果内には、測定用超音波ビームＭＢを利用して複数回得られた測定結果Ｕと、光学的な測定結果Ｐが示されている。光学的な測定結果Ｐは、実際の凝固サイズの参考値である。
【００５４】
測定用超音波ビームＭＢを利用して複数回得られた測定結果Ｕが示すように、変調周波数３４Ｈｚでは１０〜１８秒において、変調周波数６７Ｈｚでは１０〜１５秒において、変調周波数１０２Ｈｚでは４〜１０秒において、凝固の開始が検出されている。つまり、変調周波数が高いほど、凝固の開始が適切に検出されている。特に変調周波数１０２Ｈｚでは、実際の凝固サイズの参考値である光学的な測定結果Ｐと、ほぼ同じ時刻に凝固の開始が検出されている。
【００５５】
また、各実験結果内に示される直線Ｒは、振動（変位）を与えている範囲（振動範囲）を示しており、振動範囲は、変調周波数３４Ｈｚで１３ｍｍ程度、変調周波数６７Ｈｚで８ｍｍ程度、変調周波数１０２Ｈｚで７ｍｍ程度である。振動範囲内、つまり直線Ｒの下方における凝固サイズの測定結果が信頼できる領域であり、変調周波数が低いほど、広域的な凝固の検出に適している。
【００５６】
本超音波医用装置においては、比較的高い変調周波数と比較的低い変調周波数を用いて変位用超音波ビームＥＢを変調処理し、比較的高い変調周波数により発生直後の局所的な凝固の有無などを高精度に測定し、比較的低い変調周波数により進行後の広域的な凝固の大きさを高精度に測定する。したがって、本超音波医用装置においては、少なくとも２つの変調周波数を利用するが、本超音波医用装置においては、例えば、フレームレートを一定としつつ変調周波数を以下の手法で選択する。
【００５７】
図６は、変調周波数の設定の具体例を説明するための図である。本超音波医用装置においては、次式により変調周波数を決定する。
【００５８】
［数１］変調周波数（Ｈｚ）＝｛フレームレート（Ｈｚ）／素数｝×自然数Ｎ
【００５９】
数１式におけるフレームレート（Ｈｚ）と素数は、例えば装置の仕様や治療対象等に応じて適宜に設定される。以下においては、フレームレートを５００Ｈｚ、素数を１３とした具体例について説明する。
【００６０】
フレームレートが５００Ｈｚであり素数を１３とすると、数１式により得られる変調周波数（Ｈｚ）は、自然数Ｎをその変調周波数の識別番号（変調周波数Ｎ）とすると、変調周波数１（３８．４６Ｈｚ），変調周波数２（７６．９２Ｈｚ），・・・，変調周波数５（１９２．３０Ｈｚ），・・・となる。
【００６１】
フレームレートは、フレームトリガ（図２）の周期であり、これは測定用超音波ビームＭＢにより同じ位置において変位の測定が繰り返しされる周期である。つまり、フレームレートが変位の測定におけるサンプリングレートとなる。
【００６２】
図６において、＜Ａ＞は、変調信号の１周期の位相を示しており、位相番号（位相１〜位相１３）は、変調信号の１周期を素数１３により等分した場合の位相位置（位相角度）に対応している。＜Ｂ＞は、各変調周波数ごとに、位相番号（位相１〜位相１３）とデータの取得順であるサンプリング番号（ＳＰ１〜ＳＰ１３）の対応関係を示している。
【００６３】
素数が１３であれば、数１式により、変調周波数１（３８．４６Ｈｚ）は、フレームレート（５００Ｈｚ）の１／１３となる。つまり、変位の測定におけるサンプリングレートは変調周波数１の１３倍となる。したがって、＜Ａ＞に示す位相番号の配置で、位相１においてサンプリング番号１（ＳＰ１）のデータが得られると、位相２においてサンプリング番号２（ＳＰ２）のデータが得られ、その後も＜Ｂ＞に示すように、位相３，位相４，位相５，・・・の順にデータが得られ、位相１３においてサンプリング番号１３（ＳＰ１３）のデータが得られて、１周期分の１３個のデータが収集される。この様子を図７（Ｉ）に示す。
【００６４】
図６に戻り、変調周波数２（７６．９２Ｈｚ）は、フレームレート（５００Ｈｚ）の２／１３となる。つまり、変位の測定におけるサンプリングレートは、変調周波数２の１３／２倍となる。したがって、＜Ａ＞に示す位相番号の配置で、位相１においてサンプリング番号１（ＳＰ１）のデータが得られると、位相３においてサンプリング番号２（ＳＰ２）のデータが得られ、その後も＜Ｂ＞に示すように、位相５，位相７，位相９，・・・の順にデータが得られ、位相１３においてサンプリング番号７（ＳＰ７）のデータが得られると、次の周期の位相２においてサンプリング番号８（ＳＰ８）のデータが得られる。さらに、その後も＜Ｂ＞に示すように、位相４，位相６，位相８，・・・の順にデータが得られ、位相１２においてサンプリング番号１３（ＳＰ１３）のデータが得られる。つまり、＜Ｂ＞に示すように、サンプリング番号１〜１３（ＳＰ１〜ＳＰ１３）までのデータを得ることにより、位相１から位相１３までの１周期分のデータ（波形上は２周期分）が収集される。この様子を図７（ＩＩ）に示す。
【００６５】
図６に戻り、変調周波数５（１９２．３０Ｈｚ）は、フレームレート（５００Ｈｚ）の５／１３となる。つまり、変位の測定におけるサンプリングレートは、変調周波数５の１３／５倍となる。したがって、＜Ａ＞に示す位相番号の配置で、位相１においてサンプリング番号１（ＳＰ１）のデータが得られると、その後も＜Ｂ＞に示すように順にデータが得られて、サンプリング番号１〜１３（ＳＰ１〜ＳＰ１３）までのデータを得ることにより、位相１から位相１３までの１周期分のデータが収集される。
【００６６】
同様に、図６に例示していない他の変調周波数においても、サンプリング番号１〜１３（ＳＰ１〜ＳＰ１３）までのデータを得ることにより、位相１から位相１３までの１周期分のデータを収集することができる。
【００６７】
図６に示す変調周波数の設定の具体例によれば、比較的少ないサンプリング数（例えば１３）で位相の偏りの無いデータの収集が可能になり、また、エリアシングの問題も回避できる。なお、位相偏りの無いデータの収集により、ＲＭＳ（Root Mean Squre）の値が、十分に細かくサンプリングした時のＲＭＳの値と大きく乖離しない値となる。
【００６８】
本超音波医用装置においては、上述した数１式により変調周波数が決定され、比較的高い変調周波数と比較的低い変調周波数を用いて変位用超音波ビームＥＢを変調処理し、比較的高い変調周波数により発生直後の局所的な凝固の有無などを高精度に測定し、比較的低い変調周波数により進行後の広域的な凝固の大きさを高精度に測定する。つまり、本超音波医用装置は、比較的高い変調周波数と比較的低い変調周波数を含む複数の変調周波数を利用する。そこで、以下に、複数の変調周波数を利用した具体例について説明する。
【００６９】
図８は、複数の変調周波数の合成波を利用した具体例１を示すフローチャートである。まず、フレームレートが設定される（Ｓ７０１）。例えば診断範囲の深さを１５ｃｍとすると、超音波の往復に要する時間が１５（ｃｍ）×２／１５００（ｍ／ｓ）＝１００μｓ（マイクロ秒）となる。測定用超音波ビームＭＢが例えば２０本の場合には、フレームレートは１／（２０×１００μｓ）＝５００Ｈｚとなる。なお、測定用超音波ビームＭＢの２０本のうち、例えば４本の測定用超音波ビームＭＢが利用される。例えば、図２に示したように、測定用超音波ビームＭＢが２０本あり、そのうち測定用に４本の測定用超音波ビームＭＢが用いられる場合、温度が上昇する時間と下降する時間の比率が１６：４となり、上昇する時間が下降する時間の４倍あるので、効率的に温度上昇させることが出来る。（実際には変調波の位相によっては必ずしも上記の上昇と下降の比率とはならない時もあるが、治療時間全体をみれば、大よそ上記の関係が成立している。）
【００７０】
次に、変調周波数が設定される（Ｓ７０２）。変調の周期を測定期間（図２参照）よりも十分に低くすることにより、測定期間が振動に与える影響を軽減または回避できる。例えば、測定期間が４００μｓであると、変調周波数を２．５ｋＨｚよりも十分に低くする必要がある。また、数ｍｍ〜数ｃｍの大きさの凝固域に対して十分な感度を得るために、ずり波の波長が１０ｃｍ程度より短いことが望ましい。ずり波（横波）の音速を１ｍ／ｓとすると、ずり波の波長が１０ｃｍ以下となる変調周波数は１０Ｈｚ以上となる。また、秒オーダーの時間分解能を確保するためには、１秒間に複数回以上の振動を繰り返すことが望ましいため変調周波数は数Ｈｚ以上となる。以上の条件を満たしつつ、図６を利用して説明した具体例に従って、例えば、比較的低い変調周波数（低変調周波数）が３８．４６Ｈｚ、比較的高い変調周波数（高変調周波数）が１９２．３０Ｈｚに設定される。
【００７１】
そして、高変調周波数（１９２．３０Ｈｚ）と低変調周波数（３８．４６Ｈｚ）を合成した変調処理により得られる変位用超音波ビームＥＢにより治療部位Ｐの組織に振動を与え、測定用超音波ビームＭＢを介して変位測定用のデータが収集される（Ｓ７０３）。図６を利用して説明したように、サンプリング番号１〜１３（ＳＰ１〜ＳＰ１３）のデータが収集される。なお、サンプリング番号１〜１３の１セットのデータのみでも変調周波数の１周期分のデータを得ることができるものの、例えば、ノイズ等の影響を低減または除去するために、２セット分のデータが収集される。もちろん、２セット以上のデータが収集されてもよい。
【００７２】
以下、素数Ｍの例に関してフレームレートをＦＲ、二つの変調周波数をＦＬ、ＦＨ（但しＦＨ＞ＦＬ）、ＦＲ／Ｍ＝Ｆ１として説明する。更にＦＬ／Ｆ１＝ＮＬ、ＦＨ／Ｆ１＝ＮＨと表記する。放射力自体は音圧の二乗に比例するので、二つの変調周波数の合成波が生成されることや、これらの波のエリアシングが存在することを考慮して、選択出来る変調周波数の条件を整理しておく。変調の仕方は図９に示すＤＣ成分を持たない変調法と、ＤＣ成分を持つ変調法があり、放射力の周波数成分の現れ方が異なるので、わけて議論する。
【００７３】
まずＤＣ成分を持たない変調を与えた場合に関して説明する。駆動波形の変調周波数はＮＬ，ＮＨ。これに対して計測される振動成分は、２ＮＬ，２ＮＨ，ＤＣ，ＮＨ−ＮＬ，ＮＨ＋ＮＬ，ＮＬ−ＮＨの６種類の放射力の振動成分と、Ｍ−２ＮＨ，Ｍ−２ＮＬ，Ｍ−ＮＭ−ＮＬ，Ｍ−ＮＨ＋ＮＬと４種類のＭでサンプリングしたことに起因するエリアシング成分が存在する。（ちなみに上記ＤＣはＮＬ同士もしくはＮＨ同士の差周波成分である。また、ＮＬ−ＮＨは通常は負となるがエリアシングとの一致を考慮する場合にはこれも考える必要があるので、ここに加えてある。）これらの中でＮＨとＮＬが区別して独立に計測できる条件のみが適用可能なので、これを調べてみる。
【００７４】
まずパターン（１）として、放射力の振動成分同士が一致することがある条件を調べると、ＮＨ＞ＮＬなどを考慮して整理すると
・ＮＨ−ＮＬ＝２ＮＬとなるＮＨ＝３ＮＬとなる条件、例えばＭ＝１３の場合は（ＮＬ，ＮＨ）＝（１，３）、（２，６）、（３，９）、（４，１２）
がこの条件に該当する。
【００７５】
次にパターン（２）として放射力の振動成分とエリアシング成分が一致する条件を調べる。（エリアシング成分同士が一致する条件はパターン（１）と同じになる。）ここではＭが素数（つまり他の数の倍数にならないこと）などを使って整理すると、一致する条件が存在するのは下記の５つの条件である。
・Ｍ−２ＮＬ＝ＮＨ＋ＮＬとなるＭ＝３ＮＬ＋ＮＨとなるときで、Ｍ＝１３では（ＮＬ，ＮＨ）＝（１，１０）、（２，７）、（３，４）、（５，１１）、（６，８）。但し、あとは二つはＭの２倍。
・Ｍ−２ＮＬ＝ＮＨ−ＮＬとなるＭ＝ＮＬ＋ＮＨとなるとき、Ｍ＝１３では（ＮＬ，ＮＨ）＝（１，１２）、（２，１１）、（３，１０）、（４，９）、（５，８）、（６，７）。
・Ｍ−２ＮＨ＝ＮＨ＋ＮＬとなるＭ＝３ＮＨ＋ＮＬとなるとき、Ｍ＝１３では（ＮＬ，ＮＨ）＝（１，４）、（２，８）、（３，１２）。但し、あとは二つはＭの２倍、３倍と一致。
・Ｍ−２ＮＨ＝ＮＨ−ＮＬとなるＭ＝３ＮＨ−ＮＬとなるとき、Ｍ＝１３では（ＮＬ，ＮＨ）＝（２，５）、（５，６）。
・Ｍ−２ＮＬ＝ＮＬ−ＮＨとなるＭ＝３ＮＬ−ＮＨとなるとき、Ｍ＝１３では（ＮＬ，ＮＨ）＝（７，８）、（８，１１）。
【００７６】
これらを表に纏めたものが図１０である。図１０に示す表の中でＮＨ／ＮＬが１に近いものは二変調周波数を使う効果が少ないので、二変調周波数としては、１と５〜９、１１や２と９，１０，１２、３と８や１１などが適していることがわかる。典型例を図１１と図１２に示す。
【００７７】
次にＤＣ成分を持つ変調を与えた場合に関して説明する。駆動波形の変調周波数は先ほどと同じくＮＬ，ＮＨ。これに対して計測される振動成分は、ＮＬ，ＮＨ，２ＮＬ，２ＮＨ，ＤＣ，ＮＨ−ＮＬ，ＮＨ＋ＮＬ，ＮＬ−ＮＨの８種類の放射力の振動成分と、Ｍ−ＮＨ，Ｍ−ＮＬ，Ｍ−２ＮＨ，Ｍ−２ＮＬ，Ｍ−ＮＭ−ＮＬ，Ｍ−ＮＨ＋ＮＬ，Ｍ＋ＮＨ−ＮＬと７種類のＭでサンプリングしたことに起因するエリアシング成分が存在する。これらの中でＮＨとＮＬが区別して独立に計測できる条件のみが適用可能なので、これを調べてみる。
【００７８】
まずパターン（１）として、放射力の振動成分同士が一致することがある条件を調べると、ＮＨ＞ＮＬなどを考慮して整理すると
・ＮＨ＝２ＮＬとなる条件、例えばＭ＝１３の場合は（ＮＬ，ＮＨ）＝（１，２）、（２，４）、（３，６）、（４，８）、（５，１０）、（６，１２）。
・ＮＨ−ＮＬ＝２ＮＬとなるＮＨ＝３ＮＬとなる条件、例えばＭ＝１３の場合は（ＮＬ，ＮＨ）＝（１，３）、（２，６）、（３，９）、（４，１２）
がこの条件に該当する。
【００７９】
次にパターン（２）として放射力の振動成分とエリアシング成分が一致する条件を調べる。計算の過程は上記のＤＣ成分無の場合と類似するので省略して、結果のみ纏めると
・Ｍ＝ＮＬ＋ＮＨとなるときで、Ｍ＝１３では（ＮＬ，ＮＨ）＝（１，１２）、（２，１１）、（３，１０）、（４，９）、（５，８）、（６，７）。
・Ｍ＝２ＮＬ＋ＮＨとなるとき、Ｍ＝１３では（ＮＬ，ＮＨ）＝（１，１１）、（２，９）、（３，７）、（４，５）、（７，１２）、（８，１０）。但し、あとは二つはＭの２倍と一致。
・Ｍ＝２ＮＨ＋ＮＬとなるとき、Ｍ＝１３では（ＮＬ，ＮＨ）＝（３，５）、（１，６）、（８，９）（６，１０）、（４，１１）、（２，１２）。
・Ｍ＝２ＮＨ−ＮＬとなるとき、Ｍ＝１３では（ＮＬ，ＮＨ）＝（１，７）、（３，８）、（５，９）、（７，１０）、（９，１１）、（１１，１２）。
・Ｍ＝３ＮＬ＋ＮＨとなるとき、Ｍ＝１３では（ＮＬ，ＮＨ）＝（１，１０）、（２，７）、（３，４）、（５，１１）、（６，８）、（９，１２）。
・Ｍ＝３ＮＬ−ＮＨとなるとき、Ｍ＝１３では（ＮＬ，ＮＨ）＝（７，８）、（８，１１）。
・Ｍ＝３ＮＨ＋ＮＬとなるとき、Ｍ＝１３では（ＮＬ，ＮＨ）＝（１，４）、（５，７）、（２，８）、（９，１０）（６，１１）、（３，１２）。
・Ｍ＝３ＮＨ−ＮＬとなるとき、Ｍ＝１３では（ＮＬ，ＮＨ）＝（２，５）、（５，６）、（１０，１２）。
【００８０】
これらを表に纏めたものが図１３である。図１３に示す表の中でＮＨ／ＮＬが１に近いものは二変調周波数を使う効果が少ないので、二変調周波数としては、１と５，８，９や２と１０、３と１１などが適していることがわかる。ＤＣ成分が有る場合と無い場合を比べるとＤＣ成分が無い場合の方が変調周波数の選択肢が多いため使いやすいが、送波アンプの非線形性などにより、ＤＣ成分が入ってしまう場合があるので、アンプの性能の制約を厳しく設定するのが難しい場合は、上記のＤＣ成分が有る場合を選択した方が良い。
【００８１】
図８に戻り、Ｓ７０３で変位データが収集されると、次に、収集されたデータに基づいて変位が測定される（Ｓ７０４）。変位は、例えば測定用超音波ビームＭＢの深さ方向について、各深さごとに測定される。また、各深さごとに、例えば、１つのデータセット２６点（２６時相）の中の隣り合うデータ同士（ｎとｎ＋１，ｎ＝１，２，・・・，２５）が相互相関演算等により比較され、各深さごとに変位が算出される。なお、例えば加熱治療前の基準となる時相と最新時相との比較により変位が算出されてもよい。また、変位の算出に先立って、必要に応じて、ベースバンド除去処理やノイズ除去処理などが行われてもよい。
【００８２】
次に、高変調周波数成分と低変調周波数成分が分離される（Ｓ７０５）。例えば、バンドパスフィルタ等を利用して、高変調周波数（１９２．３０Ｈｚ）に対応した振動成分（変位成分）と低変調周波数（３８．４６Ｈｚ）に対応した振動成分（変位成分）が分離して抽出される。なお、Ｓ７０４におけるの変位の測定に先立って高変調周波数に対応したデータ成分と低変調周波数に対応したデータ成分を分離し、Ｓ７０４において各変調周波数ごとに変位を測定してもよい。ここで二周波成分の分離とＲＭＳの順序は重要である。ＲＭＳでは一度二乗されるので、二つの周波数の和周波や差周波が発生してしまい、これらがエリアシングも起こすので、図１３の利用可能な条件が更に狭まってしまう可能性があるためである。
【００８３】
次に、各変調周波数成分ごとに得られた振動成分に基づいて凝固が測定される（Ｓ７０６）。例えば、各深さごと且つ各変調周波数成分ごとに、２セット分（２フレーム分）に亘って得られた振動成分（変位成分）から、変位の二乗平均平方根（ＲＭＳ）、つまり実効値を算出し、例えばその実効値を現時相（最新時相）の振幅値とする。さらに、各深さごと且つ各変調周波数成分ごとに、例えば、現時相の振幅値が、加熱治療前の基準となる時相における振幅値の７０パーセントとなった場合に、当該現時相の当該深さにおいて凝固が始まったと判定する。また、現時相において凝固が始まったと判定された複数の深さについて、それら複数の深さの範囲から、凝固の大きさ（凝固サイズ）が算出される。なお、各深さ（各位置）ごとに凝固を検出して、その検出結果を例えば画像等にマッピングする手法をＬＭＩ（Localized Motion Imaging）と称する。
【００８４】
また、Ｓ７０６において、実効値に代えて、フィッティング（fitting）やロックイン検波を利用して現時相の振幅値を得るようにしてもよい。フィッティングやロックイン検波を利用する場合には、これらの利用の前に、データの並べ替えが行われる。つまり、図６を利用して説明したように、例えば高変調周波数（変調周波数５：１９２．３０Ｈｚ）においてはサンプリング番号ＳＰｎの取得順、つまりＳＰ１，ＳＰ２，ＳＰ３，・・・，ＳＰ１３の順と、位相番号の順番がずれているため、図６に示すように、位相番号の順となるように、サンプリング番号ＳＰｎが並べ替えられる。
【００８５】
次に、治療部位Ｐに対して強力集束超音波（ＨＩＦＵ）が照射される（Ｓ７０７）。例えば、治療用超音波ビームＴＢにより、０．５〜１．０秒間程度ＨＩＦＵが照射される。
【００８６】
そして、加熱期間（図２参照）が終了していなければ（Ｓ７０８）、次の測定期間（図２参照）において、Ｓ７０３〜Ｓ７０６までの測定処理が実行され、測定期間後にＳ７０７においてＨＩＦＵが照射される。なお、治療用超音波ビームＴＢが測定用超音波ビームＭＢとして利用される場合には、Ｓ７０７において、高変調周波数と低変調周波数を合成したＨＩＦＵが照射され、ＨＩＦＵ照射後に残る振動成分がＳ７０３〜Ｓ７０６において測定される。
【００８７】
一方、加熱期間（図２参照）が終了したならば（Ｓ７０８）、治療部位Ｐにおける治療が終了する。また、Ｓ７０６における凝固の測定において、目標とする凝固の大きさが確認された時点で、治療部位Ｐにおける治療を終了してもよい。治療部位Ｐにおける治療が終了したならば、別の位置にある治療部位Ｐに対して治療を行うようにしてもよい。
【００８８】
図１４は、複数の変調周波数を切り替える具体例２を示すフローチャートである。図１４の具体例２においても、まず、フレームレートが設定されて（Ｓ８０１）、変調周波数が設定される（Ｓ８０２）。Ｓ８０１とＳ８０２における処理は、図８のＳ７０１とＳ７０２における処理と同じである。つまり、Ｓ８０１においてフレームレートが５００Ｈｚに設定され、Ｓ８０２において低変調周波数が３８．４６Ｈｚ、高変調周波数が１９２．３０Ｈｚに設定される。
【００８９】
図１４の具体例２においては、低変調周波数と高変調周波数について、別々にデータが収集されて変位が測定される。つまり、まず、低変調周波数を利用した変調処理により得られる変位用超音波ビームＥＢにより治療部位Ｐの組織に振動を与え、測定用超音波ビームＭＢを介して変位測定用のデータが収集される（Ｓ８０３）。例えば、図８のＳ７０３と同様に、２セット分のデータが収集される。次に、収集されたデータに基づいて変位（振動成分）が測定される（Ｓ８０４）。例えば、図８のＳ７０４と同様な処理により、各深さごとに変位が算出される。
【００９０】
続いて、高変調周波数を利用した変調処理により得られる変位用超音波ビームＥＢにより治療部位Ｐの組織に振動を与え、測定用超音波ビームＭＢを介して変位測定用のデータが収集される（Ｓ８０５）。例えば、図８のＳ７０３と同様に、２セット分のデータが収集される。次に、収集されたデータに基づいて変位（振動成分）が測定される（Ｓ８０６）。例えば、図８のＳ７０４と同様な処理により、各深さごとに変位が算出される。
【００９１】
そして、各変調周波数成分ごとに得られた振動成分に基づいて凝固が測定される（Ｓ８０７）。例えば、図８のＳ７０６と同様な処理により、現時相において凝固が始まった深さの判定や、現時相における凝固の大きさ（凝固サイズ）が算出される。
【００９２】
次に、治療部位Ｐに対して強力集束超音波（ＨＩＦＵ）が照射される（Ｓ８０８）。例えば、治療用超音波ビームＴＢにより、０．５〜１．０秒間程度ＨＩＦＵが照射される。
【００９３】
そして、加熱期間（図２参照）が終了していなければ（Ｓ８０９）、次の測定期間（図２参照）において、Ｓ８０３〜Ｓ８０７までの測定処理が実行され、測定期間後にＳ８０８においてＨＩＦＵが照射される。なお、治療用超音波ビームＴＢが測定用超音波ビームＭＢとして利用される場合には、低変調周波数により変調処理されたＨＩＦＵ照射後にＳ８０３とＳ８０４の測定処理を実行し、高変調周波数により変調処理されたＨＩＦＵ照射後にＳ８０５とＳ８０６の測定処理を実行すればよい。
【００９４】
一方、加熱期間（図２参照）が終了したならば（Ｓ８０９）、治療部位Ｐにおける治療が終了する。また、Ｓ８０７における凝固の測定において、目標とする凝固の大きさが確認された時点で、治療部位Ｐにおける治療を終了してもよい。治療部位Ｐにおける治療が終了したならば、別の位置にある治療部位Ｐに対して治療を行うようにしてもよい。
【００９５】
図１５は、複数の変調周波数を段階的に変更する具体例３を示すフローチャートである。図１５の具体例３においても、まず、フレームレートが設定されて（Ｓ９０１）、変調周波数が設定される（Ｓ９０２）。Ｓ９０１における処理は、図８のＳ７０１における処理と同じである。つまり、Ｓ９０１においてフレームレートが５００Ｈｚに設定される。またＳ９０２においては、図８のＳ７０２における処理と同様な処理により、複数の変調周波数が設定される。例えば、数１式に基づいて、変調周波数６（２３０．７６Ｈｚ），変調周波数５（１９２．３０Ｈｚ），変調周波数４（１５３．８４Ｈｚ），変調周波数３（１１５．２８Ｈｚ），変調周波数２（７６．９２Ｈｚ），変調周波数１（３８．４６Ｈｚ）の６つの変調周波数が設定される。なお、初期状態においては、変調周波数６に設定される。
【００９６】
そして、現在設定されている変調周波数を利用した変調処理により得られる変位用超音波ビームＥＢにより治療部位Ｐの組織に振動を与え、測定用超音波ビームＭＢを介して変位測定用のデータが収集される（Ｓ９０３）。例えば、図８のＳ７０３と同様に、２セット分のデータが収集される。次に、収集されたデータに基づいて変位（振動成分）が測定される（Ｓ９０４）。例えば、図８のＳ７０４と同様な処理により、各深さごとに変位が算出される。
【００９７】
さらに、現在設定されている変調周波数成分について得られた振動成分に基づいて凝固が測定される（Ｓ９０５）。例えば、図８のＳ７０６と同様な処理により、現時相において凝固が始まった深さの判定や、現時相における凝固の大きさ（凝固サイズ）が算出される。
【００９８】
次に、治療部位Ｐに対して強力集束超音波（ＨＩＦＵ）が照射される（Ｓ９０６）。例えば、治療用超音波ビームＴＢにより、０．５〜１．０秒間程度ＨＩＦＵが照射される。
【００９９】
続いて、全ての変調周波数に関する処理が終了したか否かが確認される（Ｓ９０７）。つまり、６つの変調周波数に関する処理が全て終了したか否かが確認され、終了していなければ、Ｓ９０８において１つだけ低い変調周波数に変更され、その変更後の変調周波数によりＳ９０３〜Ｓ９０６までの処理が行われる。
【０１００】
Ｓ９０３〜Ｓ９０８の処理が繰り返し実行され、Ｓ９０７において、６つ全ての変調周波数に関する処理の終了が確認されると、Ｓ９０９に処理が進められる。
【０１０１】
そして、加熱期間（図２参照）が終了していなければ（Ｓ９０９）、Ｓ９０２において変調周波数が初期状態の変調周波数６に設定され、Ｓ９０３〜Ｓ９０８までの処理が実行される。なお、治療用超音波ビームＴＢが測定用超音波ビームＭＢとして利用される場合には、Ｓ９０６において、１つだけ低い変調周波数に変更され、その変更後の変調周波数により変調処理されたＨＩＦＵが照射される。
【０１０２】
一方、加熱期間（図２参照）が終了したならば（Ｓ９０９）、治療部位Ｐにおける治療が終了する。また、Ｓ９０５における凝固の測定において、目標とする凝固の大きさが確認された時点で、治療部位Ｐにおける治療を終了してもよい。治療部位Ｐにおける治療が終了したならば、別の位置にある治療部位Ｐに対して治療を行うようにしてもよい。
【０１０３】
図１６は、凝固の判定結果に基づいて変調周波数を切り替える具体例４を示すフローチャートである。図１６の具体例４においても、まず、フレームレートが設定されて（Ｓ１００１）、変調周波数が設定される（Ｓ１００２）。Ｓ１００１における処理は、図８のＳ７０１における処理と同じである。つまり、Ｓ１００１においてフレームレートが５００Ｈｚに設定される。またＳ１００２においては、図８のＳ７０２における処理と同様な処理により、複数の変調周波数が設定される。例えば、数１式に基づいて、変調周波数６（２３０．７６Ｈｚ），変調周波数５（１９２．３０Ｈｚ），変調周波数４（１５３．８４Ｈｚ），変調周波数３（１１５．２８Ｈｚ），変調周波数２（７６．９２Ｈｚ），変調周波数１（３８．４６Ｈｚ）の６つの変調周波数が設定される。なお、初期状態においては、変調周波数６に設定される。
【０１０４】
そして、現在設定されている変調周波数を利用した変調処理により得られる変位用超音波ビームＥＢにより治療部位Ｐの組織に振動を与え、測定用超音波ビームＭＢを介して変位測定用のデータが収集される（Ｓ１００３）。例えば、図８のＳ７０３と同様に、２セット分のデータが収集される。次に、収集されたデータに基づいて変位（振動成分）が測定される（Ｓ１００４）。例えば、図８のＳ７０４と同様な処理により、各深さごとに変位が算出される。
【０１０５】
さらに、現在設定されている変調周波数により得られた振動成分に基づいて凝固が測定される（Ｓ１００５）。例えば、図８のＳ７０６と同様な処理により、現時相において凝固が始まった深さの判定や、現時相における凝固の大きさ（凝固サイズ）が算出される。
【０１０６】
次に、治療部位Ｐに対して強力集束超音波（ＨＩＦＵ）が照射される（Ｓ１００６）。例えば、治療用超音波ビームＴＢにより０．５〜１．０秒間程度ＨＩＦＵが照射される。
【０１０７】
そして、現在設定されている変調周波数について、凝固サイズがその変調周波数に対応した閾値に達したか否かが確認される（Ｓ１００７）。具体例４においては、各変調周波数ごとに、凝固サイズに関する閾値が設定される。
【０１０８】
図１７は、凝固サイズと変調周波数の対応関係を示す図である。例えば、凝固サイズが０（未検出）から２ｍｍまでの範囲においては、変調周波数６（２３０．７６Ｈｚ）が利用され、凝固サイズが２ｍｍから５ｍｍまでの範囲においては、変調周波数５（１９２．３０Ｈｚ）が利用される。他の凝固サイズについても、図１７に示すとおりである。
【０１０９】
図１６に戻り、変調周波数６（２３０．７６Ｈｚ）を利用した凝固の測定においては、凝固サイズの２ｍｍが閾値とされ、Ｓ１００７において、凝固サイズが２ｍｍ以上または凝固サイズが２ｍｍより大きくなると、凝固が閾値に達したと判定される。
【０１１０】
Ｓ１００７において凝固が閾値に達すると、全ての変調周波数に関する処理が終了したか否かが確認される（Ｓ１００８）。つまり、６つの変調周波数に関する処理が全て終了したか否かが確認され、終了していなければ、Ｓ１００９において１つだけ低い変調周波数に変更され、その変更後の変調周波数によりＳ１００３〜Ｓ１００７までの処理が行われる。
【０１１１】
変調周波数５（１９２．３０Ｈｚ）を利用した凝固の測定においては、凝固サイズの５ｍｍが閾値とされ（図１７参照）、Ｓ１００７において、凝固サイズが５ｍｍ以上または凝固サイズが５ｍｍより大きくなると、凝固が閾値に達したと判定される。
【０１１２】
また、変調周波数４（１５３．８４Ｈｚ）を利用した凝固の測定においては、凝固サイズの８ｍｍが閾値とされ、変調周波数３（１１５．２８Ｈｚ）を利用した凝固の測定においては、凝固サイズの１５ｍｍが閾値とされ、変調周波数１（３８．４６Ｈｚ）を利用した凝固の測定においては、治療目標サイズが閾値とされる（図１７参照）。
【０１１３】
こうして、Ｓ１００３〜Ｓ１００９までの処理が繰り返され、変調周波数１（３８．４６Ｈｚ）を利用した凝固の測定が行われ、Ｓ１００７において治療目標サイズに達したことが確認されると、Ｓ１００８において全変調周波数の終了が確認され、治療部位Ｐにおける治療が終了する。なお、加熱期間（図２参照）が終了した場合にも、治療部位Ｐにおける治療を終了してもよい。治療部位Ｐにおける治療が終了したならば、別の位置にある治療部位Ｐに対して治療を行うようにしてもよい。
【０１１４】
図１８は、図１の超音波医用装置（本超音波医用装置）により形成される凝固状態画像の具体例を示す図である。凝固画像形成部２６は、凝固測定部２５から得られる凝固の測定結果に基づいて、＜Ｍ＞に示す凝固状態画像を形成する。なお、図１８において、横軸はＨＩＦＵによる加熱時間を示しており、縦軸は深さを示している。焦点はＨＩＦＵを照射する治療用超音波ビームＴＢの焦点である。また、斜線で示される部分は凝固が確認された領域を示している。
【０１１５】
＜Ｌ＞は、低変調周波数（３８．４６Ｈｚ）による凝固の測定結果を示している。低変調周波数では、振動領域が比較的広いため、測定できる凝固サイズ、つまり、焦点を中心とした縦軸方向の広がりが比較的大きい。但し、＜Ｌ＞内において破線丸印で示す領域において凝固が検出できていない。
【０１１６】
一方、＜Ｈ＞は、高変調周波数（１９２．３０Ｈｚ）による凝固の測定結果を示している。高変調周波数では、焦点を中心とした振動領域が比較的小さく局所的になり、その局所的な振動領域内における組織の平均弾性率の変化が大きく、早期の小さな凝固の検出に適している。そのため、＜Ｈ＞の高変調周波数では、＜Ｌ＞内の破線丸印で示す領域に対応する時間において、凝固が検出されている。但し、＜Ｈ＞内において破線丸印で示す振動領域を超える部分においては、振動が与えられていないため、凝固の検出に適していない。
【０１１７】
凝固画像形成部２６は、低変調周波数の測定結果＜Ｌ＞と高変調周波数の測定結果＜Ｈ＞に基づいて、例えば測定結果＜Ｌ＞と測定結果＜Ｈ＞を合成して、＜Ｍ＞の凝固状態画像を形成する。
【０１１８】
＜Ｍ＞の凝固状態画像では、高変調周波数（１９２．３０Ｈｚ）の成分に基づいて得られる測定結果から早期の小さな凝固が検出されており、さらに、低変調周波数（３８．４６Ｈｚ）の成分に基づいて得られる測定結果から比較的広い領域に亘って凝固が検出されている。
【０１１９】
したがって、＜Ｍ＞の凝固状態画像から、ユーザ（検査者）は、例えば発生直後の局所的な凝固の有無や凝固のタイミングなどを確認でき、さらに、例えば進行後の広域的な凝固の大きさなどを確認することができる。
【０１２０】
なお、３つ以上の変調周波数を利用した場合には、３以上の変調周波数の各々から得られる測定結果を合成して、凝固状態画像が形成される。
【０１２１】
これまで、本発明の好適な実施形態として、振幅に変調を与えることで放射力に変調を与える手法に関して説明を行ってきた。しかし放射力を変調する手法としては他にも次の方法がある。各素子に与える電圧の包絡線振幅は時間軸上で一定として、ビームの焦点位置を振動させる方法である。つまり、音の伝搬方向をｘ方向として、焦点位置のｙ座標をｆｙとするとｆｙ＝ｆＡｓｉｎ（ωｔ）となるように駆動する。ｆＡがビーム幅と同程度もしくはそれより大きいと各点における放射力はωの角周波数で変調されたことになる。駆動電圧を振幅変調する手法では送波回路にトランスを用いる手法があるが、この段落で説明した方法では駆動電圧は一定で良いので、トランスが不要となる。この結果回路規模が小さく出来るので、装置実用化上はメリットが大きい。ゆえに、これまでに説明した放射力の変調方法として上記の手法を用いることも実用上は有効である。
【０１２２】
以上、本発明の好適な実施形態を説明したが、上述した実施形態は、あらゆる点で単なる例示にすぎず、本発明の範囲を限定するものではない。本発明は、その本質を逸脱しない範囲で各種の変形形態を包含する。なお、本発明に係る超音波医用装置を利用した治療等は、医師等の専門家の指導の下で十分に慎重に行われるべきことは言うまでもない。
【符号の説明】
【０１２３】
１０  超音波プローブ、２０  測定診断ブロック、２２  送受信部、２４  変位測定部、２５  凝固測定部、２６  凝固画像形成部、２８  超音波画像形成部、３０  治療放射ブロック、３２  治療用送信部、３４  変位用送信部、３６  変調周波数制御部、４０  制御部、５０  表示部。
【要約】
【課題】超音波を利用した組織の凝固の測定において測定精度を高める。
【解決手段】変調周波数制御部３６は、比較的高い変調周波数と比較的低い変調周波数を用いて変位用超音波ビームＥＢを変調処理するように変位用送信部３４を制御する。変位測定部２４は、各変調周波数ごとに治療部位Ｐにおける組織の変位を測定し、凝固測定部２５は、比較的高い変調周波数による変位の測定結果に基づいて、治療部位Ｐにおける局所的な凝固を測定し、比較的低い変調周波数による変位の測定結果に基づいて、治療部位Ｐにおける広域的な凝固を測定する。これにより、例えば発生直後の局所的な凝固の有無などを高精度に測定でき、さらに、例えば進行後の広域的な凝固の大きさなどを高精度に測定できる。
【選択図】図１【図１】
0005590493.tif 000002
【図２】
0005590493.tif 000003
【図３】
0005590493.tif 000004
【図４】
0005590493.tif 000005
【図５】
0005590493.tif 000006
【図６】
0005590493.tif 000007
【図７】
0005590493.tif 000008
【図８】
0005590493.tif 000009
【図９】
0005590493.tif 000010
【図１０】
0005590493.tif 000011
【図１１】
0005590493.tif 000012
【図１２】
0005590493.tif 000013
【図１３】
0005590493.tif 000014
【図１４】
0005590493.tif 000015
【図１５】
0005590493.tif 000016
【図１６】
0005590493.tif 000017
【図１７】
0005590493.tif 000018
【図１８】
0005590493.tif 000019