特許 5954841 超音波医用装置

(19)【発行国】日本国特許庁(JP)

(12)【公報種別】特許公報(B2)

(11)【特許番号】5954841

(24)【登録日】2016年6月24日

(45)【発行日】2016年7月20日

(54)【発明の名称】超音波医用装置

(51)【国際特許分類】

   A61B 17/00 20060101AFI20160707BHJP

   A61B 8/08 20060101ALI20160707BHJP

【ＦＩ】

   A61B17/00 700

   A61B8/08

【請求項の数】6

【全頁数】20

(21)【出願番号】特願2014-168139(P2014-168139)

(22)【出願日】2014年8月21日

(65)【公開番号】特開2016-42944(P2016-42944A)

(43)【公開日】2016年4月4日

【審査請求日】2015年7月23日

(73)【特許権者】

【識別番号】504137912

【氏名又は名称】国立大学法人 東京大学

(73)【特許権者】

【識別番号】000005108

【氏名又は名称】株式会社日立製作所

(74)【代理人】

【識別番号】110001210

【氏名又は名称】特許業務法人ＹＫＩ国際特許事務所

(72)【発明者】

【氏名】東 隆

(72)【発明者】

【氏名】杉山 隆介

(72)【発明者】

【氏名】清水 悠輝

(72)【発明者】

【氏名】藤原 圭祐

(72)【発明者】

【氏名】射谷 和徳

【審査官】 村上 聡

(56)【参考文献】

【文献】 特開２０１３−２４８１４１（ＪＰ，Ａ）

【文献】 特開２０１３−０５５９８４（ＪＰ，Ａ）

【文献】 特開２０１３−０４３０８２（ＪＰ，Ａ）

(58)【調査した分野】（Int.Cl.，ＤＢ名）

Ａ６１Ｂ １７／００

Ａ６１Ｂ ８／０８

(57)【特許請求の範囲】

【請求項1】

組織を周期的に変位させる変位用超音波の送信信号を出力する送信部と、
周期的に変位する組織から測定用超音波の受信信号を得る受信部と、
測定用超音波の受信信号に基づいて組織の変位を測定する変位測定部と、
測定された変位に基づいて、当該変位の周期的な変化を示す変位データを得る変位処理部と、
変位データに基づいて組織の凝固を診断する凝固診断部と、
を有し、
前記送信部は、変調処理された変位用超音波の送信信号を出力し、
前記変位処理部は、前記変位データとして、前記変調処理の周期に対応した時間軸に沿って変位の時間変化を示した変位マップを形成し、
前記凝固診断部は、前記変位マップに基づいて組織の凝固を診断する、
ことを特徴とする超音波医用装置。

【請求項2】

請求項１に記載の超音波医用装置において、
前記変位処理部は、サンプリング方向に並ぶ複数のサンプリング点について、各サンプリング点における変位の時間変化を示した変位マップを形成し、
前記凝固診断部は、前記変位マップに基づいて、組織のサンプリング方向の凝固サイズを導出する、
ことを特徴とする超音波医用装置。

【請求項3】

請求項２に記載の超音波医用装置において、
前記凝固診断部は、前記変位マップを時間軸方向に周波数解析することにより前記変調処理の周期に対応した周波数成分を抽出し、抽出した周波数成分に基づいて組織の凝固サイズＶｓを導出する、
ことを特徴とする超音波医用装置。

【請求項4】

請求項３に記載の超音波医用装置において、
前記凝固診断部は、複数のサンプリング点について抽出した前記周波数成分のサンプリング方向における位相の変化に基づいて、組織のサンプリング方向の凝固サイズＶｓを導出する、
ことを特徴とする超音波医用装置。

【請求項5】

請求項２から４のいずれか１項に記載の超音波医用装置において、
前記凝固診断部は、複数の診断時刻に亘って各診断時刻ごとに得られる前記変位マップに基づいて、各サンプリング点ごとに変位の減少率を算出し、複数のサンプリング点についての変位の減少率に基づいて、組織のサンプリング方向の凝固サイズＡｒを導出する、
ことを特徴とする超音波医用装置。

【請求項6】

請求項２から５のいずれか１項に記載の超音波医用装置において、
前記凝固診断部は、前記変位マップから得られる前記変調処理の周期に対応した周波数成分に基づいて組織の凝固サイズＶｓを導出し、前記変位マップから得られる変位の減少率に基づいて組織の凝固サイズＡｒを導出し、凝固サイズＶｓと凝固サイズＡｒに基づく総合的な判定により組織の凝固サイズを決定する、
ことを特徴とする超音波医用装置。

【発明の詳細な説明】

【技術分野】

【0001】

本発明は、組織の凝固を診断する超音波医用装置に関する。

【背景技術】

【0002】

強力集束超音波（ＨＩＦＵ：High Intensity Focused Ultrasound）を例えば生体に照射し、その音響エネルギーを利用して腫瘍などの治療部位を加熱して凝固させる治療方法が知られている。

【0003】

組織が加熱されて凝固する際に、その凝固の前後において、組織の弾性率（ヤング率）が増大することが知られている。また、ＨＩＦＵなどの比較的強力な超音波は、その進行方向に向かって放射力を発生させるため、例えば、ＨＩＦＵの超音波ビームの焦点部位の組織に１０〜１００μｍ（マイクロメートル）程度の変位を与えることができる。

【0004】

そのため、ＨＩＦＵなどの比較的強力な超音波により組織に変位を与えて、弾性率の増大による変位の減少を測定して、組織の凝固を観察することができる。例えば、ＨＩＦＵの超音波を変調周波数ｆ_Ｍで振幅変調し、放射力の強度を変動させることにより、焦点部位の組織に振動を励起させ、その振動の変位または速度を超音波診断装置によって測定する。

【0005】

この原理を用いて組織の凝固を検出し、その検出結果を画像上にマッピングする手法がＨＭＩ（Harmonic Motion Imaging）である（特許文献１，２参照）。放射力や組織の振動の周波数が上記変調周波数ｆ_Ｍの２倍になることからＨＭＩと呼ばれている。

【先行技術文献】

【特許文献】

【0006】

【特許文献1】米国特許出願公開第２００５／０００４４６６号明細書

【特許文献2】米国特許出願公開第２００７／０２７６２４２号明細書

【発明の概要】

【発明が解決しようとする課題】

【0007】

ところで、ＨＩＦＵを利用した治療では、治療部位における凝固の状況に応じてＨＩＦＵを適切に制御することが望ましい。これは生体へのＨＩＦＵの照射においては、その伝搬経路の音響特性が患者個々人によって異なることに起因して減衰や位相歪が焦点音圧に与える影響が異なることや、焦点近傍での血流の違いによる冷却効果の違いに起因する凝固に必要なエネルギーの変化などのために、治療のためのパラメータを事前に最適化することが難しいためである。例えば、治療部位において凝固が開始されるタイミングや、目的とする大きさまで凝固が終了するタイミングなどを知ることができれば、それらのタイミングなどに応じてＨＩＦＵの照射を制御することが可能になる。

【0008】

このような状況において、本願の発明者らは、超音波を利用して組織の凝固を診断する技術について研究開発を重ねてきた。

【0009】

本発明は、その研究開発の過程において成されたものであり、その目的は、超音波を利用して組織の凝固を診断する改良技術を提供することにある。

【課題を解決するための手段】

【0010】

上記目的にかなう好適な超音波医用装置は、組織を周期的に変位させる変位用超音波の送信信号を出力する送信部と、周期的に変位する組織から測定用超音波の受信信号を得る受信部と、測定用超音波の受信信号に基づいて組織の変位を測定する変位測定部と、測定された変位に基づいて、当該変位の周期的な変化を示す変位データを得る変位処理部と、変位データに基づいて組織の凝固を診断する凝固診断部と、を有することを特徴とする。

【0011】

上記装置において、測定用超音波は、一般的な超音波診断装置における診断用の超音波と同程度の特性（周波数，波形，強度）であり、診断用の超音波振動子を利用して送受することができる。これに対し、変位用超音波は、放射力により組織を変位させる程度の比較的強度の大きい超音波である。変位用超音波は、測定用超音波に比べて強度が大きく、例えば強力集束超音波（ＨＩＦＵ：High Intensity Focused Ultrasound）が変位用超音波の好適な具体例である。なお、変位用超音波として強力集束超音波（ＨＩＦＵ）を利用する場合には、その強力集束超音波（ＨＩＦＵ）により組織を変位させつつ当該組織を加熱して凝固させてもよい。

【0012】

上記装置は、変位用超音波により組織を周期的に変位させ、その変位の周期的な変化を示す変位データを得ている。これにより、組織に係る変位の周期的な変化を示した変位データに基づいて組織の凝固を診断することが可能になる。

【0013】

望ましい具体例において、前記送信部は、変調処理された変位用超音波の送信信号を出力し、前記変位処理部は、前記変位データとして、前記変調処理の周期に対応した時間軸に沿って変位の時間変化を示した変位マップを形成し、前記凝固診断部は、前記変位マップに基づいて組織の凝固を診断する、ことを特徴とする。

【0014】

望ましい具体例において、前記変位処理部は、サンプリング方向に並ぶ複数のサンプリング点について、各サンプリング点における変位の時間変化を示した変位マップを形成し、前記凝固診断部は、前記変位マップに基づいて、組織のサンプリング方向の凝固サイズを導出する、ことを特徴とする。

【0015】

望ましい具体例において、前記凝固診断部は、前記変位マップを時間軸方向に周波数解析することにより前記変調処理の周期に対応した周波数成分を抽出し、抽出した周波数成分に基づいて組織の凝固サイズＶｓを導出する、ことを特徴とする。

【0016】

望ましい具体例において、前記凝固診断部は、複数のサンプリング点について抽出した前記周波数成分のサンプリング方向における位相の変化に基づいて、組織のサンプリング方向の凝固サイズＶｓを導出する、ことを特徴とする。

【0017】

望ましい具体例において、前記凝固診断部は、複数の診断時刻に亘って各診断時刻ごとに得られる前記変位マップに基づいて、各サンプリング点ごとに変位の減少率を算出し、複数のサンプリング点についての変位の減少率に基づいて、組織のサンプリング方向の凝固サイズＡｒを導出する、ことを特徴とする。

【0018】

望ましい具体例において、前記凝固診断部は、前記変位マップから得られる前記変調処理の周期に対応した周波数成分に基づいて組織の凝固サイズＶｓを導出し、前記変位マップから得られる変位の減少率に基づいて組織の凝固サイズＡｒを導出し、凝固サイズＶｓと凝固サイズＡｒに基づく総合的な判定により組織の凝固サイズを決定する、ことを特徴とする。

【発明の効果】

【0019】

本発明により、超音波を利用して組織の凝固を診断する改良技術が提供される。例えば本発明の好適な態様によれば、組織に係る変位の周期的な変化を示した変位データに基づいて組織の凝固を診断することが可能になる。

【図面の簡単な説明】

【0020】

【図1】本発明の実施において好適な超音波医用装置の全体構成を示す図である。

【図2】図１の超音波医用装置の動作を示すタイミングチャートである。

【図3】変調周波数の設定の具体例を説明するための図である。

【図4】位相１から位相１３までのデータが収集される様子を示す図である。

【図5】図１の超音波医用装置の動作の具体例を示すフローチャートである。

【図6】変位マップの具体例を示す図である。

【図7】凝固の進行に伴う変位マップの変化を説明するための図である。

【図8】変位マップに基づいた凝固サイズＶｓの測定例を示す図である。

【図9】変位マップに基づいた凝固サイズＡｒの測定例を示すフローチャートである。

【図10】変位減少率マップの具体例を示す図である。

【図11】治療領域内における治療部位Ｐの具体例を示す図である。

【図12】凝固サイズＡｒの信頼性を説明するための図である。

【発明を実施するための形態】

【0021】

図１は、本発明の実施において好適な超音波医用装置（本超音波医用装置）の全体構成図である。本超音波医用装置は、複合型の超音波振動子１０を有しており、超音波振動子１０は、ＨＩＦＵ用振動子１０Ｈと診断用振動子１０Ｄを備えている。

【0022】

ＨＩＦＵ用振動子１０Ｈは、強力集束超音波（ＨＩＦＵ）を送波する振動子であり、例えば二次元的に配列された複数の振動素子を備えている。ＨＩＦＵ用振動子１０Ｈは、例えば癌や腫瘍などの治療部位Ｐに向けて治療用超音波ビームＴＢを形成して強力集束超音波を送波し、その治療部位Ｐを加熱して治療するために利用される。

【0023】

また、ＨＩＦＵ用振動子１０Ｈは、治療部位Ｐに向けて変位用超音波ビームＥＢを形成して変位発生用の超音波を送波し、治療部位Ｐにおいて放射力を発生させ組織を変位させる。変位用超音波ビームＥＢは、治療部位Ｐにおいて有効な放射力を発生させる程度の強度で形成されるビームであり、例えば、治療用超音波ビームＴＢを変位用超音波ビームＥＢとして利用してもよい。もちろん、治療用超音波ビームＴＢとは異なる変位用超音波ビームＥＢが利用されてもよい。

【0024】

一方、診断用振動子１０Ｄは、例えば二次元的に配列された複数の振動素子を備えており、例えば治療部位Ｐを有する被検体（患者）に対して、超音波画像を形成するための比較的弱い超音波を送受する。つまり、公知の一般的な超音波診断装置と同じ程度の強度（エネルギー）の超音波を送受する。

【0025】

また、診断用振動子１０Ｄは、治療部位Ｐに向けて測定用超音波ビームＭＢを形成して測定用の超音波を送受し、測定用超音波ビームＭＢに沿って受信信号を得る。測定用超音波ビームＭＢに沿って得られる受信信号は、変位用超音波ビームＥＢの放射力による治療部位Ｐにおける変位の測定に利用される。

【0026】

なお、超音波振動子１０は、例えば、お椀（どんぶり）状に凹ませた内部の表面を振動子面とする。そして、例えば、お椀状に凹んだ内部の中央に位置する底の部分に診断用振動子１０Ｄが設けられ、診断用振動子１０Ｄを取り囲むようにＨＩＦＵ用振動子１０Ｈが設けられる。なお、超音波振動子１０の振動子面の形状は、お椀状に限定されず、例えば治療の用途等に応じた形状とされることが望ましい。また、全ての振動素子またはいくつかの振動素子が、ＨＩＦＵ用と診断用の両用途に併用されてもよい。

【0027】

測定診断ブロック２０は、診断用振動子１０Ｄの送受信を制御する送受信部２２を備えている。送受信部２２は、診断用振動子１０Ｄを構成する複数の振動素子の各々に対応した送信信号を出力することにより、診断用振動子１０Ｄを制御して送信ビームを形成し、さらに、それら複数の振動素子の各々から得られる受信信号に対して整相加算処理などを施すことにより、受信ビームに沿って受信信号を得る。つまり、送受信部２２は、送信部（送信ビームフォーマ）の機能と、受信部（受信ビームフォーマ）の機能を備えている。

【0028】

送受信部２２は、治療部位Ｐを含んだ三次元空間内または断面内で診断用の超音波ビームを走査させて画像用の受信信号を収集する。そして、収集された受信信号に基づいて、超音波画像形成部２８が三次元の超音波画像または二次元の断層画像の画像データを形成し、その画像データに対応した超音波画像が表示部５０に表示される。

【0029】

ユーザ（検査者）は、表示部５０に表示される超音波画像から、治療部位Ｐの位置等を確認し、図示しない操作デバイス等を利用して治療部位Ｐの位置情報を本超音波医用装置に入力する。もちろん、本超音波医用装置が、超音波画像に対する画像解析処理等により治療部位Ｐの位置を確認して位置情報を得るようにしてもよい。

【0030】

また、送受信部２２は、診断用振動子１０Ｄを制御して測定用超音波ビームＭＢを形成し、測定用超音波ビームＭＢに沿って受信信号を得る。そして、変位測定部２４は、測定用超音波ビームＭＢに沿って得られる受信信号に基づいて、治療部位Ｐにおける変位を測定する。また、変位マップ形成部２５は、変位測定部２４において測定された変位に基づいて、その変位の周期的な変化を示す変位データを得る変位処理部として機能する。変位マップ形成部２５は、変位データとして、後に詳述する変位マップを形成する。そして、凝固測定部２６は、変位マップに基づいて治療部位Ｐにおける組織の凝固を測定する。変位測定部２４と変位マップ形成部２５と凝固測定部２６における具体的な処理については後に詳述する。

【0031】

一方、治療放射ブロック３０は、治療用送信部３２を備えており、治療用送信部３２はＨＩＦＵ用振動子１０Ｈを構成する複数の振動素子の各々に対応した送信信号を出力することにより、ＨＩＦＵ用振動子１０Ｈを制御して治療用超音波ビームＴＢを形成する。治療用送信部３２は、制御部４０によって制御され、例えば治療部位Ｐ内に焦点を設定した治療用超音波ビームＴＢが形成される。

【0032】

また、治療放射ブロック３０は、変位用送信部３４を備えており、変位用送信部３４はＨＩＦＵ用振動子１０Ｈを構成する複数の振動素子の各々に対応した送信信号を出力することにより、ＨＩＦＵ用振動子１０Ｈを制御して変位用超音波ビームＥＢを形成する。変位用超音波ビームＥＢは変調処理を施され、その変調処理における変調周波数は、変調周波数制御部３６によって制御される。なお、変調周波数制御部３６は、制御部４０により制御される。

【0033】

治療用超音波ビームＴＢに沿って強力集束超音波（ＨＩＦＵ）が送波されて治療部位Ｐが加熱されると治療部位Ｐの組織が凝固する。その凝固の前後において、組織の弾性率（ヤング率）が増大することが知られている。そして、組織の弾性率の変化を知るために、本超音波医用装置は、変位用超音波ビームＥＢに沿って超音波を送波して放射力を発生させ、その放射力による治療部位Ｐにおける組織の変位を測定する。その変位の測定は、測定用超音波ビームＭＢに沿って得られる受信信号に基づいて行われる。

【0034】

なお、測定診断ブロック２０内の各部と治療放射ブロック３０内の各部は、それぞれ、例えばプロセッサや電子回路等のハードウェアを利用して実現することができる。制御部４０は、例えば、演算機能を備えたハードウェアとその動作を規定するソフトウェア（プログラム）によって構成される。表示部５０は、例えば液晶ディスプレイなどである。

【0035】

また、測定診断ブロック２０は、一般的な超音波診断装置により実現されてもよい。そして、治療放射ブロック３０に対応した超音波治療装置と、測定診断ブロック２０に対応した超音波診断装置とを組み合わせたシステムにより、本超音波医用装置が具現化されてもよい。

【0036】

本超音波医用装置の全体構成は以上のとおりである。次に、本超音波医用装置により実現される機能等について説明する。なお、図１に示した構成（符号を付した各部）については、以下の説明において図１の符号を利用する。

【0037】

図２は、図１の超音波医用装置（本超音波医用装置）の動作を示すタイミングチャートである。

【0038】

メイントリガは、強力集束超音波（ＨＩＦＵ）による治療の開始タイミングを示す信号であり、例えば、ユーザ（検査者）による治療開始の操作に応じて、制御部４０から本超音波医用装置内の各部へ出力される。

【0039】

フレームトリガは、測定用超音波ビームＭＢのフレーム開始を示す信号である。送受信部２２は、例えばフレームトリガの立ち上がりのタイミングから、治療部位Ｐに向けて複数本の測定用超音波ビームＭＢを順に形成する。例えば、２つのフレームトリガの間において、治療部位Ｐに向けて１０本の送信ビームが形成され、１つの送信ビームにつき２本ずつ、２０本の受信ビームが形成される。もちろん、送信ビームと受信ビームの本数は上記の具体例に限定されない。

【0040】

加熱期間信号は、治療用超音波ビームＴＢによる治療部位Ｐの加熱処理期間を示す信号であり、加熱期間信号の立ち上がりから立下りまでの期間において、例えば治療部位Ｐを焦点として治療用超音波ビームＴＢが形成される。

【0041】

測定期間トリガは、変位測定の期間を示す信号であり、測定期間トリガの立ち上がりから立下りまでの期間が変位の測定期間となる。

【0042】

変調信号は、治療用超音波ビームＴＢの変調処理に利用される変調信号であり、例えば変調周波数制御部３６から変位用送信部３４へ出力される。

【0043】

ＨＩＦＵ信号は、治療用超音波ビームＴＢの送信信号であり、変位用送信部３４が、例えば周波数２ＭＨｚ程度の連続波を変調信号に従って振幅変調して得られる。

【0044】

変調信号は、測定期間トリガの立ち上がりから立下りまでの測定期間において、振幅が０（ゼロ）とされ、その結果、測定期間においてＨＩＦＵ信号の振幅も０となり、測定期間において治療用超音波ビームＴＢの送信が停止される。

【0045】

なお、測定期間トリガは、フレームトリガに対して遅延量（Ｄｅｌａｙ）を伴って出力される。この遅延量は、例えば、ユーザによって適宜に調整される。これにより、２つのフレームトリガの間において形成される複数本の測定用超音波ビームＭＢ（例えば２０本の受信ビーム）のうち、測定期間内に収まる１本又はいくつかの測定用超音波ビームＭＢが選択的に測定に利用される。

【0046】

本超音波医用装置は、変位を発生させる変位用超音波ビームＥＢを変調周波数ｆ_Ｍで振幅変調し、放射力の強度を変動させることにより、焦点部位の組織に振動を励起させ、その振動の変位を測定用超音波ビームＭＢを利用して測定する。その際に、変調周波数制御部３６により、変位用超音波ビームＥＢの変調周波数ｆ_Ｍが制御される。なお、振幅変調に代えて、変調周波数ｆ_Ｍの周波数変調が利用されてもよい。

【0047】

図３は、変調周波数の設定の具体例を説明するための図である。本超音波医用装置においては、例えば次式により変調周波数を決定する。

【0048】

［数１］変調周波数（Ｈｚ）＝｛フレームレート（Ｈｚ）／素数｝×自然数Ｎ

【0049】

数１式におけるフレームレート（Ｈｚ）と素数は、例えば装置の仕様や治療対象等に応じて適宜に設定される。以下においては、フレームレートを５００Ｈｚ、素数を１３とした具体例について説明する。

【0050】

フレームレートが５００Ｈｚであり素数を１３とすると、数１式により得られる変調周波数（Ｈｚ）は、自然数Ｎをその変調周波数の識別番号（変調周波数Ｎ）とすると、変調周波数１（３８．４６Ｈｚ），変調周波数２（７６．９２Ｈｚ），・・・，変調周波数５（１９２．３０Ｈｚ），・・・となる。

【0051】

フレームレートは、フレームトリガ（図２）の周期であり、これは測定用超音波ビームＭＢにより同じ位置において変位の測定が繰り返しされる周期である。つまり、フレームレートが変位の測定におけるサンプリングレートとなる。

【0052】

図３において、＜Ａ＞は、変調信号の１周期の位相を示しており、位相番号（位相１〜位相１３）は、変調信号の１周期を素数１３により等分した場合の位相位置（位相角度）に対応している。＜Ｂ＞は、各変調周波数ごとに、位相番号（位相１〜位相１３）とデータの取得順であるサンプリング番号（ＳＰ１〜ＳＰ１３）の対応関係を示している。

【0053】

素数が１３であれば、数１式により、変調周波数１（３８．４６Ｈｚ）は、フレームレート（５００Ｈｚ）の１／１３となる。つまり、変位の測定におけるサンプリングレートは変調周波数１の１３倍となる。したがって、＜Ａ＞に示す位相番号の配置で、位相１においてサンプリング番号１（ＳＰ１）のデータが得られると、位相２においてサンプリング番号２（ＳＰ２）のデータが得られ、その後も＜Ｂ＞に示すように、位相３，位相４，位相５，・・・の順にデータが得られ、位相１３においてサンプリング番号１３（ＳＰ１３）のデータが得られて、１周期分の１３個のデータが収集される。この様子を図４（Ｉ）に示す。

【0054】

図３に戻り、変調周波数２（７６．９２Ｈｚ）は、フレームレート（５００Ｈｚ）の２／１３となる。つまり、変位の測定におけるサンプリングレートは、変調周波数２の１３／２倍となる。したがって、＜Ａ＞に示す位相番号の配置で、位相１においてサンプリング番号１（ＳＰ１）のデータが得られると、位相３においてサンプリング番号２（ＳＰ２）のデータが得られ、その後も＜Ｂ＞に示すように、位相５，位相７，位相９，・・・の順にデータが得られ、位相１３においてサンプリング番号７（ＳＰ７）のデータが得られると、次の周期の位相２においてサンプリング番号８（ＳＰ８）のデータが得られる。さらに、その後も＜Ｂ＞に示すように、位相４，位相６，位相８，・・・の順にデータが得られ、位相１２においてサンプリング番号１３（ＳＰ１３）のデータが得られる。つまり＜Ｂ＞に示すように、サンプリング番号１〜１３（ＳＰ１〜ＳＰ１３）までのデータを得ることにより、位相１から位相１３までの１周期分のデータ（波形上は２周期分）が収集される。この様子を図４（ＩＩ）に示す。

【0055】

図３に戻り、変調周波数５（１９２．３０Ｈｚ）は、フレームレート（５００Ｈｚ）の５／１３となる。つまり、変位の測定におけるサンプリングレートは、変調周波数５の１３／５倍となる。したがって、＜Ａ＞に示す位相番号の配置で、位相１においてサンプリング番号１（ＳＰ１）のデータが得られると、その後も＜Ｂ＞に示すように順にデータが得られて、サンプリング番号１〜１３（ＳＰ１〜ＳＰ１３）までのデータを得ることにより、位相１から位相１３までの１周期分のデータが収集される。

【0056】

同様に、図３に例示していない他の変調周波数においても、サンプリング番号１〜１３（ＳＰ１〜ＳＰ１３）までのデータを得ることにより、位相１から位相１３までの１周期分のデータを収集することができる。

【0057】

図３に示す変調周波数の設定の具体例によれば、比較的少ないサンプリング数（例えば１３）で位相の偏りの無いデータの収集が可能になり、また、エリアシングの問題も回避できる。なお、位相偏りの無いデータの収集により、ＲＭＳ（Root Mean Squre）の値が十分に細かくサンプリングした時のＲＭＳの値と大きく乖離しない値となる。

【0058】

本超音波医用装置においては、例えば上述した数１式により変調周波数が決定され、決定された変調周波数を用いて変位用超音波ビームＥＢが変調処理される。

【0059】

図５は、図１の超音波医用装置（本超音波医用装置）の動作の具体例を示すフローチャートである。まず、フレームレートが設定される（Ｓ５０１）。例えば診断範囲の深さを１５ｃｍとすると、超音波の往復に要する時間が１５（ｃｍ）×２／１５００（ｍ／ｓ）＝１００μｓ（マイクロ秒）となる。測定用超音波ビームＭＢが例えば２０本の場合にはフレームレートは１／（２０×１００μｓ）＝５００Ｈｚとなる。なお、測定用超音波ビームＭＢの２０本のうち、例えば４本の測定用超音波ビームＭＢが利用される。例えば、図２に示したように、測定用超音波ビームＭＢが２０本あり、そのうち測定用に４本の測定用超音波ビームＭＢが用いられる場合、温度が上昇する時間と下降する時間の比率が１６：４となり、上昇する時間が下降する時間の４倍あるので、効率的に温度上昇させることが出来る。（実際には変調波の位相によっては必ずしも上記の上昇と下降の比率とはならない時もあるが、治療時間全体をみれば、大よそ上記の関係が成立している。）

【0060】

次に、変調周波数が設定される（Ｓ５０２）。変調の周期を測定期間（図２参照）よりも十分に低くすることにより、測定期間が振動に与える影響を軽減または回避できる。例えば、測定期間が４００μｓであると、変調周波数を２．５ｋＨｚよりも十分に低くする必要がある。また、秒オーダーの時間分解能を確保するためには、１秒間に複数回以上の振動を繰り返すことが望ましいため変調周波数は数Ｈｚ以上となる。以上の条件を満たしつつ、図３を利用して説明した具体例に従って、例えば、変調周波数が３８．４６Ｈｚ，７６．９２Ｈｚ，１９２．３０Ｈｚのいずれかに設定される。

【0061】

そして、変調周波数で変調処理された変位用超音波の送信信号に基づく変位用超音波ビームＥＢにより治療部位Ｐの組織に振動を与え、測定用超音波ビームＭＢを介して変位測定用の受信信号データが得られる（Ｓ５０３）。例えば図３を利用して説明したように、時間方向（周期方向）においてサンプリング番号１〜１３（ＳＰ１〜ＳＰ１３）のデータが収集される。なお、サンプリング番号１〜１３の１セットのデータのみでも変調周波数の１周期分のデータを得ることができるものの、例えば、ノイズ等の影響を低減または除去するために、２セット分のデータが収集される。もちろん、２セット以上のデータが収集されてもよい。

【0062】

変位測定用の受信信号データが得られると、変位測定部２４により、受信信号データに基づいて変位が測定される（Ｓ５０４）。変位は、例えば測定用超音波ビームＭＢの深さ方向に並ぶ複数のサンプリング点について、各サンプリング点（各深さ）ごとに測定される。また、各サンプリング点ごとに、例えば、１つのデータセット２６点（２６時相）の中の隣り合うデータ同士（ｎとｎ＋１，ｎ＝１，２，・・・，２５）が相互相関演算等により比較され、各サンプリング点（各深さ）ごとに変位が算出される。例えば、深さ方向に並ぶ１０２４個のサンプリング点について、相関演算の相関窓を６４サンプリング点として、各サンプリング点ごとに変位が算出される。なお、例えば加熱治療前の基準となる時相と最新時相との比較により変位が算出されてもよい。また、変位の算出に先立って、必要に応じて、ベースバンド除去処理やノイズ除去処理などが行われてもよい。

【0063】

変位が測定されると、変位マップ形成部２５は、測定された変位に基づいて、その変位の周期的な変化を示す変位マップを形成する（Ｓ５０５）。さらに、凝固測定部２６は、変位マップに基づいて、治療部位Ｐにおける組織の凝固を測定する（Ｓ５０６）。

【0064】

そして、治療部位Ｐに対して強力集束超音波（ＨＩＦＵ）が照射される（Ｓ５０７）。例えば、治療用超音波ビームＴＢにより０．５〜１．０秒間程度ＨＩＦＵが照射される。さらに、加熱期間（図２参照）が終了していなければ（Ｓ５０８）、次の測定期間（図２参照）において、Ｓ５０３〜Ｓ５０６までの測定処理が実行され、測定期間後にＳ５０７においてＨＩＦＵが照射される。なお、治療用超音波ビームＴＢが変位用超音波ビームＥＢとして利用される場合には、Ｓ５０７において、変調処理されたＨＩＦＵが照射され、ＨＩＦＵ照射後に残る振動成分がＳ５０３〜Ｓ５０６において測定される。

【0065】

一方、加熱期間（図２参照）が終了したならば（Ｓ５０８）、治療部位Ｐにおける治療が終了する。また、Ｓ５０６における凝固の測定において、目標とする凝固の大きさが確認された時点で、治療部位Ｐにおける治療を終了してもよい。治療部位Ｐにおける治療が終了したならば、別の位置にある治療部位Ｐに対して治療を行うようにしてもよい。

【0066】

図６は、変位マップの具体例を示す図である。変位マップ形成部２５は、変位測定部２４において測定された変位に基づいて、例えば、図６に示す変位マップを形成する。図６には、横軸を時間として縦軸を深さとした変位マップが図示されている。

【0067】

図６の縦軸方向には、測定用超音波ビームＭＢの深さ方向に並ぶ複数のサンプリング点（例えば１０２４個のサンプリング点）が並べられ、各サンプリング点（各深さ）ごとに２周期分（図３のサンプリング番号ＳＰ１〜ＳＰ１３のセットが２セット分）の変位が図６の横軸方向に並べられる。

【0068】

なお、変位マップを形成する際には、時間軸方向において変位の並べ替えが行われる。例えば、変調周波数として図３において説明した１９２．３０Ｈｚ（変調周波数５）を利用した場合には、サンプリング番号ＳＰｎの取得順、つまりＳＰ１，ＳＰ２，ＳＰ３，・・・，ＳＰ１３の順と、位相番号の順番がずれているため、図３に示すように、位相番号の順となるように、サンプリング番号ＳＰｎが時間方向において並べ替えられて、図６の横軸に対応付けられる。

【0069】

そして、図６においては、深さ方向の各サンプリング点ごとに、つまり各深さごとに、時間軸方向の２周期に亘る変位の大きさ（変位量）が輝度の大きさに対応付けてマッピングされている。つまり、図６において、正方向に大きな変位ほど高い輝度で明るく表現され、負方向に大きな変位ほど低い輝度で暗く表現されている。なお、輝度に代えて、例えば色等により変位の大きさが表現されてもよい。また、変位マップが表示部５０に表示されてもよい。

【0070】

図６の具体例において、破線で示す深さｒにおける変位の変化を時間軸方向に見ると、変位が２周期に亘って正弦波状に変化している。例えば、深さｒの位置を焦点として、変調処理されたＨＩＦＵ（変調処理された変位用超音波）が照射されると、そのＨＩＦＵの振動に追従するように深さｒの位置における組織が変位し、図６に示す具体例のような変位マップが得られる。

【0071】

なお、変位マップに対して、有効振動領域が設定されてもよい。例えばメディアンフィルタ等の処理を施した変位マップ内において、変位が２周期分存在し且つ変位の大きさが閾値以上となる領域が有効振動領域とされる。そして、変位マップに基づく後述する凝固の測定において、例えば、有効振動領域内におけるデータのみを利用するようにしてもよい。これにより、有効振動領域外のデータに伴う測定精度の悪化を抑えることができる。

【0072】

ＨＩＦＵの照射により組織の凝固が進行すると変位マップに変化が現れる。凝固測定部２６は、変位マップに基づいて組織の凝固を測定する。

【0073】

図７は、凝固の進行に伴う変位マップの変化を説明するための図である。図７（Ａ）には、ＨＩＦＵの照射開始直後における変位マップの具体例が図示されている。例えば、深さ方向のビーム長がＨＬ（例えば焦点領域がＨＬ）であるＨＩＦＵを照射した直後において、図７（Ａ）の変位マップが得られる。ビーム長ＨＬの領域（例えば焦点領域）においては、ＨＩＦＵの振動に追従するように組織が変位する。

【0074】

また、ＨＩＦＵの照射開始直後の図７（Ａ）では、組織が凝固していない又は凝固が小さいため、ビーム長ＨＬの領域外において、ビーム長ＨＬから離れるに従って、ＨＩＦＵの振動から徐々に遅れて組織が変位する。つまり、ＨＩＦＵの振動により、ビーム長ＨＬの領域においてせん断波（ずり波）が発生し、そのせん断波が深さ方向に（浅い方と深い方のそれぞれに向かって）進み、変位マップ内においてせん断波の伝播Ｓａに伴う変位の傾きが現れる。

【0075】

これに対し、図７（Ｂ）には、凝固後における変位マップの具体例が図示されている。例えば、図７（Ａ）に示すビーム長ＨＬのＨＩＦＵを照射し続けると、ビーム長ＨＬの領域から凝固が進行し、例えば図７（Ｂ）に示す凝固部分Ｃの範囲まで組織が凝固する。凝固部分Ｃにおいては、凝固により組織が硬く、せん断波が即伝播するため、例えば図７（Ｂ）に示すように、凝固部分Ｃが全体的にＨＩＦＵの振動に追従するように変位する。また、図７（Ｂ）において、凝固していない部分では、せん断波の伝播Ｓｂに伴う変位が現れる。

【0076】

図８は、変位マップに基づいた凝固サイズＶｓの測定例を示す図である。図８（ａ）には、変位マップの具体例が図示されている。凝固測定部２６は、変位マップを時間軸方向に周波数解析することにより、変位用超音波（ＨＩＦＵ）に対する変調処理の周期に対応した周波数成分を抽出し、抽出した周波数成分に基づいて組織の凝固サイズＶｓを導出する。

【0077】

図８（ｂ）は、図８（ａ）の変位マップに対する時間軸方向の周波数解析により得られる処理結果の具体例である。凝固測定部２６は、図８（ａ）の変位マップにおいて、深さ方向に並ぶ複数のサンプリング点について、各サンプリング点ごとに時間軸方向にＦＦＴ処理を行って、図８（ｂ）の処理結果を得る。

【0078】

図８（ｂ）の縦軸は、図８（ａ）と同じ深さに対応しており、図８（ｂ）の横軸は周波数を示している。そして、図８（ｂ）内における輝度がＦＦＴ処理後の絶対値に対応している。

【0079】

凝固測定部２６は、図８（ｂ）の処理結果において、変位用超音波（ＨＩＦＵ）の変調処理における変調周波数に対応した周波数成分（変調周波数成分）を抽出する。そして、凝固測定部２６は、深さ方向に並ぶ複数のサンプリング点において変調周波数成分を抽出し、変調周波数成分のサンプリング方向における位相の変化に基づいて、組織のサンプリング方向の凝固サイズＶｓを導出する。

【0080】

図８（ｃ）には、ＦＦＴ処理により得られる変調周波数成分の位相が図示されている。図８（ｃ）の横軸は、深さ方向つまり図８（ａ）（ｂ）の縦軸に対応しており、図８（ｃ）の縦軸は、位相をラジアン単位で示している。

【0081】

せん断波の伝播速度が早い部分、つまり組織が凝固している部分においては、その部分内でほぼ同位相で組織が振動するため、図８（ｃ）に示す位相波形の傾きが小さくなる。一方、せん断波の伝播速度が遅い部分、つまり組織が凝固していない部分においては、せん断波の伝播に時間がかかるため、場所ごとに（深さごとに）位相にずれが生じて、図８（ｃ）に示す位相波形の傾きが大きくなる。

【0082】

そこで、凝固測定部２６は、図８（ｃ）に示す位相波形に基づいて、ＨＩＦＵの焦点位置を含む波形部分において、波形の傾きが小さい部分（例えば閾値よりも傾きが小さい部分）を凝固領域と判定する。これにより、例えば、図８（ｃ）に示すように、ＨＩＦＵの焦点位置を含む凝固領域が判定され、深さ方向の凝固サイズＶｓが導出される。

【0083】

なお、図８に示す具体例においては、ＦＦＴ処理を行っているため信号対雑音比の向上が期待される。また、図８に示す具体例においては、凝固部分の凹凸を、この区間を構成するデータ点の標準偏差または分散等の統計的な指標等を用いて、推定精度を評価すること等も可能であり、その推定精度を後に詳述する総合判定に利用してもよい。

【0084】

また、凝固測定部２６は、変位マップに基づいて、深さ方向の各サンプリング点ごとに変位の減少率を算出し、複数のサンプリング点についての変位の減少率に基づいて、組織の凝固サイズＡｒを測定してもよい。

【0085】

図９は、変位マップに基づいた凝固サイズＡｒの測定例を示す図である。図９のフローチャートは、例えば、治療部位Ｐに対する強力集束超音波（ＨＩＦＵ）の照射と並行して実行される。

【0086】

まず、初期の診断時刻において変位マップが形成される（Ｓ９０１）。例えば、強力集束超音波（ＨＩＦＵ）の照射直前または照射直後が初期の診断時刻とされ、変位マップ（図６参照）が形成される。

【0087】

さらに、次の診断時刻において変位マップが形成される（Ｓ９０２）。例えば、初期（１時刻前）の診断時刻から一定時間Ｔが経過した次の診断時刻において、変位マップ（図６参照）が形成される。

【0088】

そして、Ｓ９０２の診断時刻における変位マップ内の変位量最大値が、それよりも１時刻前の診断時刻における変位マップの変位量最大値よりも小さいか否かが確認される（Ｓ９０３）。例えば、変位マップの有効振動領域（図６参照）内において最も大きな変位が変位量最大値とされる。

【0089】

１時刻前の診断時刻における変位マップの変位量最大値よりも小さくない場合には、Ｓ９０２に戻り、さらに次の診断時刻において変位マップが形成され、Ｓ９０３において変位量最大値が確認される。Ｓ９０２とＳ９０３における処理は、１時刻前の診断時刻における変位マップの変位量最大値よりも小さくなるまで繰り返し実行される。

【0090】

Ｓ９０３において１時刻前の診断時刻における変位マップの変位量最大値よりも小さいことが確認されると、リファレンスとなる変位マップが決定される（Ｓ９０４）。例えばＳ９０３において確認された１時刻前の診断時刻における変位マップ、つまり変位がピークに到達した時刻における変位マップがリファレンスとされる。なお、例えば、変位がピークに到達した時刻とその前後の時刻における変位マップを平均化して、リファレンスとなる変位マップが決定されてもよい。

【0091】

治療組織の初期温度や組織性状等によっては、ＨＩＦＵが照射されることにより、変位量最大値が上昇してから減少する場合もある（下記文献１参照）。つまり、ＨＩＦＵの照射直前または直後に変位量最大値のピークがあるとは限らない。Ｓ９０２からＳ９０４までの処理によれば、ＨＩＦＵが照射されることにより変位量最大値が上昇してから減少する場合においても、変位量最大値がピークに到達した時刻における変位マップをリファレンスとすることができる。
文献１「E Sapin-de Brosses, J Gennisson, M Pernot, M Fink and M Tanter, “Temperature dependence of the shear modulus of soft tissues assessed by ultrasound.” Phys Med Biol. 2010 Mar 21; 55(6):1701-18」（文献１には、生体組織の温度を変化させたときの横弾性係数を測定した結果が記載されている。生体組織によっては、温度を上昇させていくと、ある温度までは横弾性係数が徐々に小さくなり（組織が柔らかくなる＝組織の変位量が大きくなる）、ある温度から高温になるにつれて横弾性係数が大きくなる（組織が硬くなる＝組織の変位量が小さくなる）という結果が文献１に示されている。）

【0092】

そして、複数の診断時刻に亘って各診断時刻において変位マップが形成され（Ｓ９０５）、各診断時刻における変位マップとリファレンスとなる変位マップとの比較から得られる変位の減少率に基づいて、治療部位Ｐにおける組織の凝固サイズＡｒが測定される（Ｓ９０６）。凝固サイズＡｒの測定において、凝固測定部２６は、変位減少率マップを形成することが望ましい。

【0093】

図１０は、変位減少率マップの具体例を示す図である。凝固測定部２６は、例えば、図１０に示す変位減少率マップを形成する。図１０には、ＨＩＦＵが照射される時間つまり診断時刻を横軸とし、縦軸を深さとした変位減少率マップが図示されている。縦軸の深さは、変位マップ（図６）の縦軸に対応している。

【0094】

凝固測定部２６は、各診断時刻ごとに得られる変位マップについて、その変位マップの各深さごとにその深さに対応した変位量を導出する。例えば、各診断時刻ごとに得られる変位マップ内において、各深さごとに２周期分（図６の横軸）に亘る変位の変化から、その変位の二乗平均平方根（ＲＭＳ）つまり実効値を算出し、その実効値をその深さの変位量とする。もちろん、二乗平均平方根以外の演算により各深さの変位量が決定されてもよい。

【0095】

また、凝固測定部２６は、リファレンスとなる変位マップに基づいて、各深さごとにその深さのリファレンスとなる変位量を導出する。例えば、リファレンスとなる変位マップの各深さごとに２周期分（図６の横軸）に亘る変位の変化から、その変位の二乗平均平方根（ＲＭＳ）つまり実効値を算出し、その実効値をその深さにおけるリファレンスの変位量とする。もちろん、二乗平均平方根以外の演算により各深さにおけるリファレンスの変位量が決定されてもよい。

【0096】

そして、凝固測定部２６は、各診断時刻において、各深さごとに得られる変位量（Ｄｍ）と、その深さにおけるリファレンスの変位量（Ｄｒ）との比率（Ｄｍ／Ｄｒ）を算出して、その深さの変位減少率とする。そして、複数の深さについて複数の診断時刻に亘って変位減少率が算出され、変位減少率に応じた大きさの輝度をマッピングすることにより、図１０の変位減少率マップが形成される。なお、輝度に代えて、例えば色等により変位減少率が表現されてもよい。また、変位減少率マップが表示部５０に表示されてもよい。

【0097】

図１０の変位減少率マップは、リファレンスの変位量からの減少が小さいほど、つまりリファレンスの変位量に近いほど、高い輝度で明るく表現され、リファレンスの変位量からの減少が大きいほど、つまりリファレンスの変位量よりも小さいほど、低い輝度で暗く表現されている。図１０において暗く表現される部分は、凝固の進行により組織が硬くなり、変位量が減少していることが予想される。つまり、図１０において暗く表現される部分において組織が凝固していると考えられる。

【0098】

なお、変位減少率マップに対して、有効振動領域が設定されてもよい。例えば変位マップの有効振動領域（図６参照）を変位減少率マップに適用してもよいし、変位減少率マップ内において減少率の有効性を判定するなどして有効な減少率が得られる領域を有効振動領域としてもよい。

【0099】

凝固測定部２６は、変位減少率マップに基づいて、各診断時刻ごとに、組織の深さ方向の凝固サイズＡｒを導出する。例えば、図１０に示す変位減少率マップの有効振動領域内において減少率が小さい部分（例えば閾値よりも減少率が小さい部分）を凝固領域と判定し、凝固領域の深さ方向の長さを凝固サイズＡｒとする。

【0100】

以上に説明したように、凝固測定部２６は、変位マップに基づいて、凝固サイズＶｓ（図８参照）と凝固サイズＡｒ（図１０参照）を導出する。さらに、凝固測定部２６は、複数の診断時刻において、各診断時刻ごとに凝固サイズＶｓと凝固サイズＡｒの両方を算出し、各診断時刻において凝固サイズＶｓと凝固サイズＡｒに基づく総合的な判定により組織の凝固サイズＣｏａｇを決定してもよい。例えば、凝固サイズＶｓと凝固サイズＡｒの各々に対する重み付け係数Ｗ１とＷ２（Ｗ１＋Ｗ２＝１）を用いて、総合的な組織の凝固サイズＣｏａｇが「Ｃｏａｇ＝Ｗ１・Ｖｓ＋Ｗ２・Ａｒ」とされる。

【0101】

重み付け係数Ｗ１とＷ２は、例えば、デフォルト値をＷ１＝Ｗ２＝０．５とし、以下に説明する条件に応じて適宜変更されてもよい。

【0102】

＜条件１＞ＨＩＦＵの深さ方向のビーム長よりも小さな凝固サイズを計測する場合
凝固サイズＶｓは、ＨＩＦＵの深さ方向のビーム長ＨＬ（図７）よりも凝固サイズが大きくならないと信頼性が低いため、条件１の場合には、凝固サイズＡｒの方を尊重する。例えば「Ｃｏａｇ＝Ａｒ（Ｗ１＝０，Ｗ２＝１）」とする。

【0103】

＜条件２＞ＨＩＦＵが照射される目標部位の隣接部位が既に凝固している場合
ＨＩＦＵを利用した治療においては、治療の対象となる領域において、複数の治療部位Ｐに対して次々にＨＩＦＵが照射される場合があり、治療部位Ｐの近傍において既に凝固した部位があると、治療部位Ｐにおける変位の減少率の導出において信号対雑音比が低くなる可能性がある。

【0104】

図１１は、治療領域内における治療部位Ｐの具体例を示す図である。図１１（ａ）においては、これからＨＩＦＵの照射を予定している治療部位Ｐに隣接する部位が既に凝固（凝固部位）している。そのため、治療部位Ｐにおいて減少率に基づいて得られる凝固サイズＡｒよりも凝固サイズＶｓの方を尊重する。例えば「Ｃｏａｇ＝Ｖｓ（Ｗ１＝１，Ｗ２＝０）」とする。

【0105】

一方、図１１（ｂ）においては、これからＨＩＦＵの照射を予定している治療部位Ｐが凝固部位から比較的離れているため、治療部位Ｐにおいて、凝固サイズＶｓと凝固サイズＡｒの両方を利用し、例えば「Ｃｏａｇ＝Ｗ１・Ｖｓ＋Ｗ２・Ａｒ（Ｗ１＋Ｗ２＝１）」とする。

【0106】

なお、治療部位Ｐと凝固部位との間の距離Ｌに応じて、重み付け係数Ｗ１とＷ２が変更されてもよい。例えば、距離Ｌが小さいほどＷ２を０．５よりも小さくしてＷ１を０．５よりも大きくする。また、距離ＬがＨＩＦＵビームの方位方向におけるビーム幅の何倍かに応じて重み付け係数Ｗ１とＷ２が変更されてもよい。

【0107】

また、複数の治療部位Ｐに対して次々にＨＩＦＵを照射するにあたり、常に凝固部位に隣接する治療部位Ｐを対象とする場合には、最初の治療部位Ｐに対する照射において、例えば「Ｃｏａｇ＝Ｗ１・Ｖｓ＋Ｗ２・Ａｒ（Ｗ１＋Ｗ２＝１）」とし、２つ目以降の治療部位Ｐに対する照射において常に「Ｃｏａｇ＝Ｖｓ（Ｗ１＝１，Ｗ２＝０）」としてもよい。

【0108】

＜条件３＞凝固サイズＶｓと凝固サイズＡｒについて信頼性が評価できる場合
凝固サイズＶｓと凝固サイズＡｒの少なくとも一方について、計測結果の信頼性が評価できる場合には、その信頼性に応じて、重み付け係数Ｗ１とＷ２が変更されてもよい。

【0109】

例えば、図８を利用して説明した具体例においては、凝固部分の凹凸を、この区間を構成するデータ点の標準偏差または分散等の統計的な指標等を用いて、凝固サイズＶｓの推定精度を評価すること等も可能であるため、その推定精度に応じて、重み付け係数Ｗ１とＷ２が設定されてもよい。

【0110】

また、図１０を利用して説明した具体例においては、減少率が小さい部分（例えば閾値よりも減少率が小さい部分）を凝固領域と判定し、凝固領域の深さ方向の長さを凝固サイズＡｒとして計測しているが、減少率の変化の度合いに応じて、凝固サイズＡｒの信頼性が評価されてもよい。

【0111】

図１２は、凝固サイズＡｒの信頼性を説明するための図である。図１２には、変位減少率マップの具体例が図示されている。つまり、ＨＩＦＵが照射される時間（診断時刻）を横軸とし、縦軸を深さとした変位減少率マップが図示されている。この変位減少率マップの時刻Ｔにおける変位減少率の具体例が（ａ）（ｂ）に図示されている。（ａ）（ｂ）には、横軸を深さとして縦軸を変位減少率とした波形Ａ，波形Ｂが示されている。

【0112】

波形Ａは、凝固部分とそれ以外の部分の境界における変位減少率の変化が比較的急峻であり、例えば閾値により比較的高い精度で境界を検出することができる。そのため、波形Ａに基づいて測定される凝固サイズＡｒの信頼性が高い。

【0113】

これに対し、波形Ｂは、凝固部分とそれ以外の部分の境界における変位減少率の変化が緩やかであり、例えば閾値の大きさに応じて境界の位置が比較的大きく変動する。そのため、波形Ｂに基づいて測定される凝固サイズＡｒは、閾値の大きさに応じて比較的大きく変動してしまい信頼性が低い。

【0114】

そこで、凝固測定部２６は、変位減少率マップの診断時刻における深さ方向の変位減少率の変化の度合い、例えば、凝固部分とそれ以外の部分の境界付近における変化の急峻さに応じて、凝固サイズＡｒの信頼性を評価し、その信頼性に応じて、重み付け係数Ｗ１とＷ２を変更するようにしてもよい。例えば、境界付近における変化が急峻であるほど、凝固サイズＡｒの信頼性が高いと判定され、重み付け係数Ｗ２の値が大きくされる。

【0115】

また、凝固測定部２６は、組織の凝固サイズＣｏａｇを決定するにあたり、治療対象に応じて凝固サイズＣｏａｇを決定するようにしてもよい。

【0116】

例えば、治療対象が癌の場合には、癌細胞（癌組織）を確実に凝固壊死させることが望ましいため、「Ｃｏａｇ＝Ｖｓ∩Ａｒ」とする。つまり、凝固サイズＶｓの範囲と凝固サイズＡｒの範囲の積集合（両方に属する領域）を凝固サイズＣｏａｇとする。なお、凝固サイズＶｓと凝固サイズＡｒのうちのいずれか小さい方の値を凝固サイズＣｏａｇとしてもよい。これにより、癌細胞の残存を低減して望ましくは癌細胞の残存を無くして、治療の安全性と信頼性を高めることが可能になる。

【0117】

これに対し、治療対象が子宮筋腫のような良性腫瘍においては、癌の場合とは異なり、治療領域を縮小することが望ましい場合がある。つまり、癌細胞の場合のように確実に凝固壊死に導く必要性が低く、また、まばらに照射点を設定してＨＩＦＵを照射することも考えられる。そこで、治療対象が例えば良性腫瘍の場合には「Ｃｏａｇ＝Ｖｓ∪Ａｒ」とする。つまり、凝固サイズＶｓの範囲と凝固サイズＡｒの範囲の和集合（少なくとも一方に属する領域）を凝固サイズＣｏａｇとする。なお、凝固サイズＶｓと凝固サイズＡｒのうちのいずれか大きい方の値を凝固サイズＣｏａｇとしてもよい。

【0118】

以上、本発明の好適な実施形態を説明したが、上述した実施形態は、あらゆる点で単なる例示にすぎず、本発明の範囲を限定するものではない。本発明は、その本質を逸脱しない範囲で各種の変形形態を包含する。なお、本発明に係る超音波医用装置を利用した治療等は、医師等の専門家の指導の下で十分に慎重に行われるべきことは言うまでもない。

【符号の説明】

【0119】

１０ 超音波振動子、２０ 測定診断ブロック、２２ 送受信部、２４ 変位測定部、２５ 変位マップ形成部、２６ 凝固測定部、２８ 超音波画像形成部、３０ 治療放射ブロック、３２ 治療用送信部、３４ 変位用送信部、３６ 変調周波数制御部、４０ 制御部、５０ 表示部。

【図1】

【図2】

【図3】

【図4】

【図5】

【図9】

【図11】

【図12】

【図6】

【図7】

【図8】

【図10】