(57)【要約】【課題】基準深度より深い位置にサンプルゲートが設定される場合について、ドプラ測定の精度を向上させることを目的とする。【解決手段】サンプルゲートの深さから基準深度（１繰り返し周期で超音波パルスが往復する深さ）を引いた深さの位置にはサブゲートと呼ばれる領域が生じる。パルスドプラ装置は、サンプルゲートで反射して生じる観測対象超音波を受信しながらも、サブゲートで反射して生じる不要超音波を減衰させる超音波ビームを形成する。すなわち、整相重み付け部２４が実行する処理によって、サンプルゲートに受信ビームの受信焦点が設定され、サブゲートに受信ビームの受信減衰点が設定される。整相重み付け部２４は、受信焦点で反射した超音波に基づく検波信号を強め合うように合成し、かつ、受信減衰点で反射した超音波に基づく検波信号を弱め合うように合成した合成信号を出力する。【選択図】図１
2018007952.tif 000002
【特許請求の範囲】
【請求項１】
複数の振動素子を含むプローブから被検体に所定の繰り返し周波数で超音波パルスを送信する送信部と、
前記被検体内で反射した超音波を前記プローブを介して受信する受信部と、
各前記振動素子で受信された超音波に基づく各受信信号に重み付け処理を施す重み付け処理部と、
前記重み付け処理が施された各前記受信信号を合成する合成処理部と、
前記合成処理部で生成された合成信号に基づいて前記被検体に対するドプラ測定を行うドプラ測定部と、を備え、
前記重み付け処理部は、
前記プローブを介して受信される超音波の受信ビームに受信減衰点が生じる重み付け処理を、各前記受信信号に対して施すことを特徴とする超音波診断装置。
【請求項２】
請求項１に記載の超音波診断装置において、
前記重み付け処理部は、
前記受信ビームに受信焦点が生じる重み付け処理を各前記受信信号に対して施し、
前記受信減衰点の位置を、前記受信焦点の位置および前記繰り返し周波数に基づいて定めることを特徴とする超音波診断装置。
【請求項３】
請求項１または請求項２に記載の超音波診断装置において、
前記重み付け処理部は、
各前記受信信号に対して重み付け係数を作用させる係数作用部を備え、
各前記重み付け係数は、
前記受信減衰点で反射した超音波に基づく前記合成信号が所定値未満または最小となるという条件に基づいて定められることを特徴とする超音波診断装置。
【請求項４】
請求項１から請求項３のいずれか１項に記載の超音波診断装置において、
前記受信部から前記合成処理部に至る信号経路に設けられ、各前記受信信号を個別に遅延させることで前記受信ビームを形成する遅延部を備えることを特徴とする超音波診断装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、超音波診断装置に関し、特に、被検体に対してドプラ測定を行う装置に関する。
【背景技術】
【０００２】
ドプラ法を用いて被検体の血流速度等を測定する超音波診断装置として、パルスドプラ装置が広く用いられている。パルスドプラ装置は、パルス変調された超音波を所定のパルス繰り返し周波数（PRF：Pulse Repetition Frequency：以下、単に繰り返し周波数という。）で被検体に向けて送信する。測定対象の領域としてサンプルゲートが設定され、パルスドプラ装置は、サンプルゲートに応じたタイミングで受信された超音波のドプラシフト周波数成分に基づいてドプラ波形（速度対時間波形）を表示する。
【０００３】
パルスドプラ装置では、測定可能なドプラシフト周波数、すなわち、測定可能な速度が超音波パルスの繰り返し周波数に応じて制限される。繰り返し周波数を大きくすると測定可能な速度の上限は大きくなり、繰り返し周波数を小さくすると測定可能な速度の上限は小さくなる。血流速度等の測定対象の速度が測定可能な速度の範囲を超えると、いわゆる折り返し現象が生じ、測定誤差が生じる。したがって、原則として繰り返し周波数で規定される範囲内で血流速度等の測定が行われる。
【０００４】
以下の特許文献１および２にはパルスドプラ装置が記載されている。特許文献３には、本発明に関連する処理として、超音波診断装置で送受信される超音波についてのビームフォーミング処理が記載されている。
【先行技術文献】
【特許文献】
【０００５】
【特許文献１】特開平１１−１６９３７４号公報
【特許文献２】特開平０２−３６８５２号公報
【特許文献３】特開２０１５−１９８６２号公報
【発明の概要】
【発明が解決しようとする課題】
【０００６】
測定可能な速度の上限を大きくするため、繰り返し周波数を高くするＨＰＲＦ測定法が考えられている。ＨＰＲＦ測定法では、繰り返し周期（繰り返し周波数の逆数）の１周期で超音波パルスが往復する基準深度より深い位置にサンプルゲートが設定される。そのため、先に送信した超音波パルスがサンプルゲートで反射し、この反射による観測対象超音波が受信される前に次の超音波パルスが送信される。
【０００７】
ＨＰＲＦ測定法では、後発の超音波パルスの一部が基準深度より浅い位置で反射して生じる不要超音波が十分小さい場合や、不要超音波が観測対象超音波と時間帯を重ねて受信されない場合には、十分な測定精度が得られる。しかし、被検体によっては、この条件が必ずしも成立せず、測定精度が低下することがある。
【０００８】
本発明は、基準深度より深い位置にサンプルゲートが設定される場合について、ドプラ測定の精度を向上させることを目的とする。
【課題を解決するための手段】
【０００９】
本発明は、複数の振動素子を含むプローブから被検体に所定の繰り返し周波数で超音波パルスを送信する送信部と、前記被検体内で反射した超音波を前記プローブを介して受信する受信部と、各前記振動素子で受信された超音波に基づく各受信信号に重み付け処理を施す重み付け処理部と、前記重み付け処理が施された各前記受信信号を合成する合成処理部と、前記合成処理部で生成された合成信号に基づいて前記被検体に対するドプラ測定を行うドプラ測定部と、を備え、前記重み付け処理部は、前記プローブを介して受信される超音波の受信ビームに受信減衰点が生じる重み付け処理を、各前記受信信号に対して施すことを特徴とする。
【００１０】
本発明における重み付け処理によれば、プローブを介して受信される超音波の受信ビームに受信減衰点が生じる。したがって、受信減衰点の位置で反射して生じる不要な超音波に基づく成分が受信信号において低減される。
【００１１】
望ましくは、前記重み付け処理部は、前記受信ビームに受信焦点が生じる重み付け処理を各前記受信信号に対して施し、前記受信減衰点の位置を、前記受信焦点の位置および前記繰り返し周波数に基づいて定める。
【００１２】
本発明によれば、受信焦点の位置および繰り返し周波数に基づいて受信減衰点の位置が定められ、受信焦点の位置および繰り返し周波数に対して適切な位置に受信減衰点を有する受信ビームが形成される。
【００１３】
望ましくは、前記重み付け処理部は、各前記受信信号に対して重み付け係数を作用させる係数作用部を備え、各前記重み付け係数は、前記受信減衰点で反射した超音波に基づく前記合成信号が所定値未満または最小となるという条件に基づいて定められる。
【００１４】
本発明においては、受信減衰点で反射した超音波に基づく合成信号が所定値未満または最小となるという条件に基づいて重み付け係数が定められ、この重み付け係数を各受信信号に作用させる。重み付け係数を定める処理には、例えば、ＭＶ法（Minimum Variance法）、最小ノルム法、線形予測法、ＭＵＳＩＣ法（Multiple Signal Classification法）等の数値解析法を用いることができるため、重み付け係数を定めるに際して数値計算による処理が可能となる。
【００１５】
望ましくは、前記受信部から前記合成処理部に至る信号経路に設けられ、各前記受信信号を個別に遅延させることで前記受信ビームを形成する遅延部を備える。
【００１６】
本発明においては、受信部から合成処理部に至る信号経路に遅延部が設けられていればよく、設計に自由度がある。
【発明の効果】
【００１７】
本発明によれば、基準深度より深い位置にサンプルゲートが設定される場合について、ドプラ測定の精度を向上させることができる。
【図面の簡単な説明】
【００１８】
【図１】パルスドプラ装置の構成を示す図である。
【図２】受信ビーム、サンプルゲート、サブゲート、受信焦点および受信減衰点を断層画像と共に概念的に示す図である。
【図３】従来のパルスドプラ装置において測定精度の低下が生じることを説明するタイミングチャートである。
【図４】整相重み付け部の構成を示す図である。
【図５】シミュレーション結果を示す図である。
【図６】シミュレーション結果を示す図である。
【図７】複数ゲート測定法を説明する図である。
【図８】整相重み付け部の変形例を示す図である。
【発明を実施するための形態】
【００１９】
（１）パルスドプラ装置の構成
図１には、本発明の実施形態に係るパルスドプラ装置の構成が示されている。パルスドプラ装置は、送信部１８、プローブ１０、受信部２２、整相重み付け部２４、ドプラ測定部２６および表示部２８を備える。パルスドプラ装置は、パルス変調された超音波を所定の繰り返し周波数（ＰＲＦ）でプローブ１０から被検体１４に向けて送信し、被検体１４内で反射した超音波を受信する。測定対象の領域としてサンプルゲートがユーザの操作によって設定され、パルスドプラ装置は、サンプルゲートに応じたタイミングで受信された超音波のドプラシフト周波数成分に基づいてドプラ波形画像を生成し、表示部２８に表示する。
【００２０】
パルスドプラ装置の具体的な構成について説明する。被検体１４を診断するに際してプローブ１０は被検体１４の表面に接触した状態とされる。プローブ１０は複数の振動素子１２を備える。送信部１８は、制御部２０による制御に基づいてプローブ１０の各振動素子１２に送信信号を出力する。この送信信号は、例えば、矩形の時間波形を有するパルス信号に応じてオンオフされた正弦波信号、すなわち、パルス変調された電気信号である。各振動素子１２が送信信号に応じて超音波を発生することで、プローブ１０から超音波が送信される。
【００２１】
超音波の送信に際して送信部１８は、各振動素子１２に出力する送信信号の遅延時間またはレベルを制御部２０による制御に従って調整し、プローブ１０から送信される超音波による送信ビームを形成する。すなわち、制御部２０は送信部１８に対する制御によって、プローブ１０において送信ビームを形成する。
【００２２】
なお、送信部１８は、各振動素子１２に出力する送信信号の遅延時間またはレベルを制御部２０による制御に従って調整し、送信ビームの方向または放射位置を調整してもよい。すなわち、制御部２０は送信部１８に対する制御によって、プローブ１０において形成した送信ビームの方向または放射位置を調整してもよい。
【００２３】
被検体１４内で反射した超音波がプローブ１０の各振動素子１２で受信されると、各振動素子１２は、受信された超音波に応じた電気信号を受信部２２に出力する。受信部２２は、制御部２０による制御に従い、各振動素子１２から出力された受信信号を増幅し、複数の振動素子１２に対応する複数チャネルの受信信号を整相重み付け部２４に出力する。
【００２４】
送信部１８は、送信信号と同一の周波数を有するローカル信号Ｌを整相重み付け部２４に出力する。さらに、制御部２０は、各受信信号に対する処理を実行するための各種情報を整相重み付け部２４に出力する。
【００２５】
整相重み付け部２４は、後述するように、複数チャネルの受信信号のそれぞれを適切な遅延時間だけ遅延させ、さらに、重み付け係数を作用させることで合成信号を生成しドプラ測定部２６に出力する。合成信号は、受信ビーム１６上に設定された受信焦点で反射した超音波に基づく受信信号が強め合うように合成され、かつ、受信ビーム１６上に設定された受信減衰点で反射した超音波に基づく受信信号が弱め合うように合成されることで得られる信号である。
【００２６】
ドプラ測定部２６は、合成信号に対して周波数解析処理を施し、ドプラ波形を表す画像データを生成する。ドプラ波形は、横軸に時間を取り、縦軸に速度、ドプラシフト周波数、またはこれらに相当する値を取った時間波形である。ドプラ測定部２６は、ドプラ波形画像データを表示部２８に出力する。表示部２８は、ドプラ波形画像データに基づいてドプラ波形を表示する。なお、ドプラ測定部２６は、ドプラ波形を数値によって表すデータを生成してもよい。この数値データは、例えば、時間に対して速度、ドプラシフト周波数、またはこれらに相当する値を対応づけたデータである。ドプラ測定部２６は、ドプラ波形を数値によって表示する画像データを生成し表示部２８に出力する。表示部２８は、この画像データに基づいてドプラ波形を数値によって表示する。
【００２７】
（２）整相重み付け部が実行する処理
整相重み付け部２４が実行する処理の技術的意義について説明する。図２には、受信ビーム３０、サンプルゲート３２、サブゲート３４、受信焦点３６および受信減衰点３８が断層画像４０と共に概念的に示されている。サンプルゲート３２は、基準深度Ｄ０よりも深い位置にある循環器４２上に設定されている。基準深度Ｄ０は、１繰り返し周期Ｔ（＝１／ＰＲＦ）で超音波が往復する距離である。
【００２８】
このように、基準深度Ｄ０よりも深い位置にサンプルゲート３２が設定され、血液等の運動物がサンプルゲート３２にある場合、先に送信した超音波パルスがサンプルゲート３２で反射し、この反射による超音波がプローブで受信される前に次の超音波パルスが送信される。
【００２９】
プローブから送信された超音波が反射しプローブに戻るまでの時間が、サンプルゲート深度時間（プローブとサンプルゲートとの間を超音波が往復する時間）から繰り返し周期Ｔを減算した時間と一致するような位置には、サブゲートと呼ばれる領域が生じる。別の表現では、プローブからの深さが、基準深度Ｄ０からサンプルゲート３２までの距離に相当する深さである位置には、サブゲートと呼ばれる領域が生じる。
【００３０】
サブゲート３４に超音波を反射する組織があった場合には、従来のパルスドプラ装置では次のような問題が生じていた。すなわち、後発の超音波パルスの一部がサブゲート３４で反射して生じる不要超音波が、観測対象の超音波と時間帯を重ねて受信され、この不要超音波が観測対象超音波に干渉して測定精度が低下するという問題が生じていた。
【００３１】
図３には、従来のパルスドプラ装置において測定精度の低下が生じることを説明するタイミングチャートが示されている。図３（ａ）には、パルスドプラ装置が送信する超音波パルスが示されている。パルスドプラ装置は、繰り返し周期Ｔで超音波パルスを繰り返し送信する。サンプルゲートは基準深度よりも深い位置に設定されている。サンプルゲートの位置に対応するサンプルゲート時間帯Ｓ０、Ｓ１、Ｓ２、・・・・・は、送信パルスが送信されてからサンプルゲート深度時間ＤＰＴだけ遅れた時間に設定されている。すなわち、送信パルスＰ１が送信されてからサンプルゲート深度時間ＤＰＴが経過した時間にサンプルゲート時間帯Ｓ１が設定されている。また、送信パルスＰ２が送信されてからサンプルゲート深度時間ＤＰＴが経過した時間にサンプルゲート時間帯Ｓ２が設定され、送信パルスＰ３が送信されてからサンプルゲート深度時間ＤＰＴが経過した時間にサンプルゲート時間帯Ｓ３が設定されている。
【００３２】
図３（ｂ）には、従来のパルスドプラ装置で受信される超音波が示されている。この図には、後発の超音波パルスの一部がサブゲートで反射して不要超音波が生じ、観測対象の超音波と時間帯を重ねて受信された場合が示されている。図３（ｃ）は、観測対象の超音波と不要超音波とを区別して示したものである。観測対象の超音波が実線で示され、不要超音波が破線で示されている。この図は、説明の便宜上、観測対象の超音波と不要超音波とを区別して描いた概念図であり，実際にはこれらの超音波を区別することは困難である。
【００３３】
図３（ｃ）に示されている観測対象超音波Ａ０は、送信パルスＰ１より１つ前に送信された送信パルスがサンプルゲートで反射して受信されたものであり、不要超音波Ｂ１は、送信パルスＰ１がサブゲートで反射して受信されたものである。観測対象超音波Ａ１は、送信パルスＰ１がサンプルゲートで反射して受信されたものであり、不要超音波Ｂ２は、送信パルスＰ２がサブゲートで反射して受信されたものである。観測対象超音波Ａ２は、送信パルスＰ２がサンプルゲートで反射して受信されたものであり、不要超音波Ｂ３は、送信パルスＰ３がサブゲートで反射して受信されたものである。観測対象超音波Ａ３は、送信パルスＰ３がサンプルゲートで反射して受信されたものであり、不要超音波Ｂ４は、送信パルスＰ４がサブゲートで反射して受信されたものである。
【００３４】
図３（ｂ）に示されている受信超音波Ｒ０は、観測対象超音波Ａ０および不要超音波Ｂ１が重畳されたものである。受信超音波Ｒ１は、観測対象超音波Ａ１および不要超音波Ｂ２が重畳されたものである。受信超音波Ｒ２は、観測対象超音波Ａ２および不要超音波Ｂ３が重畳されたものである。受信超音波Ｒ３は、観測対象超音波Ａ３および不要超音波Ｂ４が重畳されたものである。
【００３５】
このように、従来のパルスドプラ装置では、サンプルゲート時間帯に観測対象超音波に加えて不要超音波が受信され、不要超音波の影響によって測定精度が低下するという問題が生じていた。
【００３６】
そこで、本実施形態に係るパルスドプラ装置では、サンプルゲートで反射して生じる観測対象超音波を受信しながらも、サブゲートで反射して生じる不要超音波を減衰させる超音波ビームを形成する。すなわち、図３（ｃ）に示されている不要超音波Ｂ１、Ｂ２、Ｂ３、・・・・・に基づく受信信号成分を減衰させるような超音波ビームを形成する。
【００３７】
すなわち、図１の整相重み付け部２４が実行する処理によって、サンプルゲートに受信ビームの受信焦点が設定され、サブゲートに受信ビームの受信減衰点が設定される。
【００３８】
整相重み付け部２４は、各受信信号を個別に遅延させた上で直交検波し、直交検波によって得られた検波信号に対して重み付け係数を作用させて加算合計することで、合成信号を出力する。合成信号は、受信ビーム上に設定された受信焦点で反射した超音波に基づく検波信号が強め合うように合成され、かつ、受信ビーム上に設定された受信減衰点で反射した超音波に基づく検波信号が弱め合うように合成された信号である。図３（ｃ）に示されている不要超音波Ｂ１、Ｂ２、Ｂ３、・・・・・に基づく成分は、合成信号によって抑制される。受信ビームおよび受信焦点は各受信信号に対する遅延処理によって形成される。また、受信ビーム上の受信減衰点は、各受信信号に対する重み付け処理（各検波信号に対する重み付け係数の複素共役値の乗算）によって形成される。
【００３９】
（３）整相重み付け部の構成
図４には、整相重み付け部２４の構成が示されている。整相重み付け部２４は、遅延部４４−１〜４４−ｎ、直交検波部４６−１〜４６−ｎ、重み付け係数決定部４８、複素乗算部５０−１〜５０−ｎ、および合成処理部５２を備える。直交検波部４６−１〜４６−ｎのうちディジタル信号処理を行う部分、重み付け係数決定部４８、複素乗算部５０−１〜５０−ｎ、および合成処理部５２は、例えば、プロセッサによって構成される。これらの構成要素は、プロセッサが実行するプログラムによって構成してもよいし、ハードウエアである電子回路によって個別に構成してもよい。
【００４０】
遅延部４４−１〜４４−ｎには、ｎチャネルの受信信号ｒ１〜ｒｎが入力される。遅延部４４−ｐ（ｐ＝１〜ｎ）は、制御部から与えられた遅延時間τｐだけ受信信号ｒｐを遅延させ、直交検波部４６−ｐに出力する。
【００４１】
直交検波部４６−ｐは、送信超音波と同一の周波数を有するローカル信号Ｌによって受信信号ｒｐに対して直交検波を施すと共に、サンプルゲート時間帯外の時間帯でのレベルを強制的に０とした同相成分信号Ｉｐおよび直交成分信号Ｑｐを生成する。そして、複素数の検波信号ｖｐ＝Ｉｐ＋ｊ・Ｑｐを複素乗算部５０−ｐに出力する。同相成分信号Ｉｐは、ローカル信号Ｌと同位相の成分を示し、直交成分信号Ｑｐは、ローカル信号Ｌとの位相差が９０°である成分を示す。ｊは虚数単位であり、同相成分信号Ｉｐおよび直交成分信号Ｑｐは、プローブで受信された超音波のドプラシフト周波数成分である。すなわち、同相成分信号Ｉｐおよび直交成分信号Ｑｐの周波数は、プローブで受信された超音波のドプラシフト周波数に一致する。
【００４２】
重み付け係数決定部４８および複素乗算部５０−１〜５０−ｎは、各受信信号に重み付け処理を施す重み付け処理部としての機能を有する。
【００４３】
重み付け係数決定部４８は後述の係数決定処理によって重み付け係数Ｗ１〜Ｗｎを決定し、それぞれ、複素乗算部５０−１〜５０−ｎに出力する。
【００４４】
複素乗算部５０−１〜５０−ｎは、各検波信号（各受信信号）に対して重み付け係数を作用させる係数作用部として機能する。ここで、重み付け係数を作用させるとは、重み付け係数に基づいて処理対象の信号の値を変化させることをいい、本実施形態では、重み付け係数の共役複素数を検波信号に乗じることをいう。
【００４５】
複素乗算部５０−ｐ（ｐ＝１〜ｎ）は、検波信号ｖｐに重み付け係数Ｗｐの共役複素数である共役重み付け係数Ｗｐ＊を乗じて合成処理部５２に出力する。合成処理部５２は、Ｗ１＊・ｖ１〜Ｗｎ＊・ｖｎを加算合計して合成信号を生成し、図１のドプラ測定部２６に出力する。この合成信号は、サンプルゲートに受信焦点が設定され、サブゲートに受信減衰点が設定された受信ビームに従って受信された超音波に基づく信号である。合成信号では、受信減衰点において反射した不要超音波が抑制されるため、合成信号によって生成されるドプラ波形に含まれる測定誤差が減少する。
【００４６】
（４）係数決定処理
重み付け係数決定部４８が実行する係数決定処理について説明する。本実施形態に係る係数決定処理は、ＭＶ法（Minimum Variance法）に基づく処理である。演算対象の重み付け係数は、（数１）のような列ベクトルＷで表される。上付きのＴはベクトルの転置を表す。
【００４７】
【数１】
2018007952.tif 000003
【００４８】
サンプルゲートに受信焦点を設定し、サブゲートに受信減衰点を設定するようなフィルタ重み付けベクトルＷχは、（数２）および（数３）が成立するという条件の下で、（数４）で表される評価値Ｐ（Ｗ）を所定値未満または最小値とする列ベクトルＷを探索することで求められる。上付きのＨは共役転置行列であることを示す。
【００４９】
【数２】
2018007952.tif 000004
【００５０】
【数３】
2018007952.tif 000005
【００５１】
【数４】
2018007952.tif 000006
【００５２】
（数２）に含まれる列ベクトルａは（数３）によって表される基準ベクトルである。基準ベクトルａは、各遅延部によって受信焦点からの伝搬時間が揃えられた検波信号ｖ１〜ｖｎを成分とする列ベクトルに対応する。（数２）は、重み付け係数ベクトルＷが、受信焦点から到来した検波信号に対して合成信号の値を維持するという意義を有する。
【００５３】
（数４）に含まれる行列Ｒは、（数５）および（数６）によって得られる。
【００５４】
【数５】
2018007952.tif 000007
【００５５】
【数６】
2018007952.tif 000008
【００５６】
ここで、ｕｐはサブゲートに設定された受信減衰点で反射し、ｐ番目の振動素子で受信される超音波に対応する。（数６）に含まれるｔは時間を表しｆは超音波の周波数を表す。ｔｐａｔｔは、受信減衰点から第ｐ番目の振動素子に至るまでの距離を超音波が伝搬する時間である。ｔｐｆは、受信焦点から第ｐ番目の振動素子に至るまでの距離を超音波が伝搬する時間である。Δｔは受信減衰点と受信焦点との間の距離を超音波が伝搬する時間である。
【００５７】
先に送信された超音波パルスのうち受信焦点で反射してプローブに戻る超音波と、次に送信された超音波パルスのうち受信減衰点で反射してプローブに戻る超音波とが、プローブで同時に受信されるという条件では、受信焦点で反射した超音波は（数７）のように表される。
【００５８】
【数７】
2018007952.tif 000009
【００５９】
すなわち、（数６）は、受信焦点で反射した超音波がｐ番目の振動素子で受信されると同時に、受信減衰点で反射した超音波がｐ番目の振動素子で受信される場合において、受信減衰点で反射した超音波が、受信焦点で反射した超音波に対して位相が２πｆ（ｔｐａｔｔ＋Δｔ−ｔｐｆ）だけ遅れていることを示している。
【００６０】
ここでは、図１に示されているように複数の振動素子の配列方向をｘ軸方向とし、被検体１４の深さ方向をｙ軸正方向とする。また、受信ビーム１６の終点のｘ座標を０とする。被検体１４の深さ方向を示すｙ軸正方向と受信ビーム１６とのなす角度θ、ｐ番目の振動素子のｘ座標ｘｐ、受信ビーム方向における受信焦点の位置ｒｆ（サンプルゲートに設定される受信焦点の位置）、および、受信ビーム方向における受信減衰点の位置ｒａｔｔ（サブゲートに設定される受信減衰点の位置）を与えることで、ｔｐａｔｔおよびｔｐｆは、それぞれ、（数８）および（数９）のように求められる。ここで、ｃは超音波の伝搬速度である。
【００６１】
【数８】
2018007952.tif 000010
【００６２】
【数９】
2018007952.tif 000011
【００６３】
受信減衰点から受信焦点に至るまでの距離を超音波が伝搬する時間Δｔは、（数１０）で表される。
【００６４】
【数１０】
2018007952.tif 000012
【００６５】
また、受信減衰点はサブゲートに設定されため、受信減衰点の位置は（数１１）で表される。なお、（数１１）に含まれるｃ／（２ＰＲＦ）は基準深度Ｄ０に等しい。
【００６６】
【数１１】
2018007952.tif 000013
【００６７】
したがって、初めに、既知である受信焦点の位置ｒｆ、超音波の伝搬速度ｃ、および繰り返し周波数ＰＲＦを（数１１）に与えることで、サブゲートに設定される受信減衰点の位置ｒａｔｔが求められる。次に、受信焦点の位置ｒｆ、（数１１）によって既知となった受信減衰点の位置ｒａｔｔおよび超音波の伝搬速度ｃを（数１０）に与えることで受信減衰点から受信焦点に至るまでの距離を超音波が伝搬する時間Δｔが求められる。そして、受信減衰点の位置ｒａｔｔ、ｐ番目の振動素子のｘ座標ｘｐ、深さ方向と受信ビームとのなす角度θ、および超音波の伝搬速度ｃを（数８）に与えることで、受信減衰点から第ｐ番目の振動素子に至るまでの距離を超音波が伝搬する時間ｔｐａｔｔが求められる。さらに、受信焦点の位置ｒｆ、ｐ番目の振動素子のｘ座標ｘｐ、深さ方向と受信ビームとのなす角度θ、および超音波の伝搬速度ｃを（数９）に与えることで、受信焦点から第ｐ番目の振動素子に至るまでの距離を超音波が伝搬する時間ｔｐｆが求められる。
【００６８】
このように、受信焦点の位置ｒｆ、繰り返し周波数ＰＲＦ、１番目〜ｎ番目の各振動素子のｘ座標ｘ１〜ｘｎ、深さ方向と受信ビーム１６とのなす角度θ、および超音波の伝搬速度ｃが既知であれば、これら既知の値に基づいて不要超音波の遅延時間（ｔｐａｔｔ＋Δｔ−ｔｐｆ）が得られ、さらには（数６）で示されるｕｐが求められる。ｐ＝１〜ｎのそれぞれについてｕｐを求めることで、（数５）および（数４）に基づいて、ある重み付け列ベクトルＷについて評価値Ｐ（Ｗ）が求められる。フィルタ重み付けベクトルＷχは、重み付け列ベクトルＷの各成分を変化させて、各評価値Ｐ（Ｗ）が所定値未満または最小値となるような重み付け列ベクトルＷを探索することで求められる。
【００６９】
（５）重み付け係数決定部が実行する係数決定処理
重み付け係数決定部４８が実行する処理という観点から係数決定処理について説明する。重み付け係数決定部４８には、ユーザの操作に基づいて受信焦点の位置ｒｆ、繰り返し周波数ＰＲＦ、１番目〜ｎ番目の各振動素子のｘ座標ｘ１〜ｘｎ、深さ方向と受信ビームとのなす角度θ、および超音波の伝搬速度ｃが入力される。
【００７０】
重み付け係数決定部４８は、（数６）および（数８）〜（数１１）に基づいて、ｐ＝１〜ｎのそれぞれについて超音波ｕｐを求め、さらには、（数５）で表される相関行列Ｒを求める。重み付け係数決定部４８は、（数２）および（数３）が成立するという条件の下で、（数４）で表される評価値Ｐ（Ｗ）を所定値未満または最小値とする列ベクトルＷを探索することで、フィルタ重み付けベクトルＷχを求める。
【００７１】
係数決定処理では、検波信号ｖ１〜ｖｎを用いた適応的な処理を実行する必要がないため、迅速にフィルタ重み付けベクトルＷχが求められる。
【００７２】
（６）シミュレーション結果
図５には本実施形態に係るパルスドプラ装置についてのシミュレーション結果が示されている。このシミュレーション結果は、受信焦点を深さ８０ｍｍの位置に設定し、受信減衰点を深さ４０ｍｍの位置に設定したものである。図５（ａ）には、受信ビームの深さ方向について求められた受信レベルが示されている。横軸は、方位角θ＝０°の方向についての深さを示し、縦軸は、各深さの位置に点波源を配置した場合における合成信号の大きさを示す。図５（ｂ）には、受信ビームの方位角方向について求められた受信レベルが示されている。横軸は、深さｙ＝８０ｍｍについての方位角を示し、縦軸は、深さ８０ｍｍの各方位に点波源を配置した場合における合成信号の大きさを示す。
【００７３】
図６には従来技術におけるパルスドプラ装置についてのシミュレーション結果が示されている。このシミュレーション結果は、本実施形態に係るパルスドプラ装置において、フィルタ重み付けベクトルＷをハミング窓とした場合に相当する。図６（ａ）には、受信ビームの深さ方向について求められた受信レベルが示されている。横軸は、方位角θ＝０°の方向についての深さを示し、縦軸は、各深さの位置に点波源を配置した場合における合成信号の大きさを示す。図６（ｂ）には、受信ビームの方位角方向について求められた受信レベルが示されている。横軸は、深さｙ＝８０ｍｍについての方位角を示し、縦軸は、深さ８０ｍｍの各方位に点波源を配置した場合における合成信号の大きさを示す。
【００７４】
図５および図６から明らかように、本実施形態に係るパルスドプラ装置では、受信焦点で反射した超音波の受信レベルに比べて、受信減衰点で反射した超音波の受信レベルは小さい。したがって、サンプルゲートに受信焦点を設定し、サブゲートに受信減衰点を設定することで、サブゲートで反射する超音波の影響が抑制され、測定精度が向上する。
【００７５】
（７）複数ゲート測定法
本発明の実施形態に係るパルスドプラ装置は、時分割で複数の方向のそれぞれに送信ビームおよび受信ビームを向けて、各送受信ビームに対応するサンプルゲートについて血液等の運動物を測定する複数ゲート測定法を行ってもよい。複数ゲート測定法は、例えば、特許文献２に記載されている。この方法では、時分割で複数のサンプルゲートに対して運動物の速度を測定するため、複数の箇所における運動物の速度が迅速に測定される。
【００７６】
図７には、送信ビームの方向を切り換え時間Ｔａごとに方位角−θ１の方向と、方位角θ２の方向に切り換える場合の送信ビーム５４−１および５４−２が示されている。パルスドプラ装置は、送信ビーム５４−１に対して超音波パルスを送信してから切り換え時間Ｔａが経過した時に送信ビーム５４−２に対して超音波パルスを送信し、さらに、切り換え時間Ｔａが経過した時に送信ビーム５４−１に対して超音波パルスを送信し、切り換え時間Ｔａが経過した時に送信ビーム５４−２に対して超音波パルスを送信するという処理を繰り返す。
【００７７】
この場合、基準深度Ｄａは、切り換え時間Ｔａの間に超音波パルスが往復する距離Ｄａ＝ｃ・Ｔａ／２となる。送信ビーム５４−１に対するサンプルゲート５６−１、および、送信ビーム５４−２に対応するサンプルゲート５６−２が基準深度Ｄａよりも深い位置にある場合、各送受信ビームについてサブゲートが生じる。
【００７８】
一般に、送信ビームおよび受信ビームの指向性には広がりがあるため、送信ビームおよび受信ビームの指向性が最大でない方向についても超音波が送受信される。そのため、各受信ビームのサブゲートで超音波が反射する場合には不要超音波が受信されることがある。
【００７９】
そこで、各受信ビームについてサブゲートに受信減衰点を設け、パルスドプラ装置で受信される不要超音波を抑制することで、測定誤差が低減される。図７には、送信ビーム５４−１に対応するサンプルゲート５６−１およびサブゲート５８−１に、それぞれ、受信焦点および受信減衰点が設けられ、送信ビームに対応するサンプルゲート５６−２およびサブゲート５８−２に、それぞれ、受信焦点および受信減衰点が設けられた例が示されている。受信焦点および受信減衰点の位置は送受信ビームごとに異なるため、送受信ビームごとに異なるフィルタ重み付けベクトルＷχが用いられる。
【００８０】
（８）変形例
（８−１）サンプルゲートよりも深い位置にあるサブゲート
上記では、サンプルゲートよりも浅い位置に生じるサブゲートに受信減衰点を設定する例について説明した。サブゲートは、サンプルゲートよりも深い位置にも生じる。すなわち、サンプルゲートを中心にして、サンプルゲートよりも浅い位置にあるサブゲートと対称である位置にもサブゲートが生じる。サンプルゲートよりも深い位置にあるサブゲートに超音波を反射する組織があると、先発の超音波パルスがサブゲートで反射して生じる不要超音波が、後発の超音波パルスがサンプルゲートで反射して生じる観測対象超音波と時間帯を重ねて受信される。そこで、サンプルゲートよりも深い位置にあるサブゲートに受信減衰点を設けてもよい。サンプルゲートよりも深い位置にあるサブゲートに設定する受信減衰点の位置は、（数１２）によって求められる。
【００８１】
【数１２】
2018007952.tif 000014
【００８２】
サンプルゲートよりも深い位置にあるサブゲートに受信減衰点を設ける場合、サンプルゲートよりも深い位置にあるサブゲートに設定される受信減衰点につき（数６）に従って求められるｕｐを（数５）に代入して、評価値Ｐ（Ｗ）を求めればよい。
【００８３】
また、サンプルゲートよりも浅い位置および深い位置の両方のサブゲートに受信減衰点を設ける場合、サンプルゲートよりも浅い位置にあるサブゲートに設定される受信減衰点につき（数６）に従って求められるｕｐｓと、サンプルゲートよりも深い位置にあるサブゲートに設定される受信減衰点につき（数６）に従って求められるｕｐｄとを加算したｕｐ＝ｕｐｓ＋ｕｐｄを（数５）に代入して、評価値Ｐ（Ｗ）を求めればよい。
【００８４】
（８−２）フィルタ重み付けベクトルＷχを求めるその他の方法
上記では、ＭＶ法に基づく係数決定処理について説明した。受信ビーム上に受信焦点および受信減衰点が設けられるような重み付け係数を求める方法には、最小ノルム法、線形予測法、ＭＵＳＩＣ法（Multiple Signal Classification法）等、様々なものがある。
【００８５】
（８−３）整相重み付け部の変形例
遅延部４４−１〜４４−ｎは、受信部２２から合成処理部５２に至る信号経路のいずれかの位置に設けられることで、受信ビームおよび受信焦点を形成し得る。図８には、図４に示される整相重み付け部２４に対し、各遅延部の位置を変更した変形例が示されている。図４に示されている整相重み付け部２４では、各直交検波部の前段に遅延部が設けられているのに対し、図８に示されている整相重み付け部６０では、各複素乗算部の後段に遅延部が設けられている。
【００８６】
整相重み付け部６０では、重み付け係数決定部４８は、（数２）および（数３）で表される基準ベクトルａに代えて、（数１３）で表される基準ベクトルｂが用いられる。
【００８７】
【数１３】
2018007952.tif 000015
【００８８】
この基準ベクトルｂは、各遅延部で遅延処理が施されていないこと基づき、基準ベクトルａの各成分に対して、各成分の位相が進んでいるベクトルを基準とすることを意味する。基準ベクトルｂを用いることで、受信焦点から到来した検波信号に対して合成信号の値が維持される。また、ｔｐｆは、（数９）によらず、一定値０として係数決定処理が実行される。
【符号の説明】
【００８９】
１０  プローブ、１２  振動素子、１４  被検体，１６，３０  受信ビーム、１８  送信部、２０  制御部、２２  受信部、２４，６０  整相重み付け部、２６  ドプラ測定部、２８  表示部、３２，５６−１，５６−２  サンプルゲート、３４，５８−１，５８−２  サブゲート、３６  受信焦点、３８  受信減衰点、４０  断層画像、４２  循環器、４４−１〜４４−ｎ  遅延部、４６−１〜４６−ｎ  直交検波部、４８  重み付け係数決定部、５０−１〜５０−ｎ  複素乗算部、５２  合成処理部、５４−１，５４−２  送信ビーム。
【図１】
2018007952.tif 000016
【図２】
2018007952.tif 000017
【図３】
2018007952.tif 000018
【図４】
2018007952.tif 000019
【図５】
2018007952.tif 000020
【図６】
2018007952.tif 000021
【図７】
2018007952.tif 000022
【図８】
2018007952.tif 000023