(57)【特許請求の範囲】
【請求項１】
周波数を周期的に変化させる連続波の変調送信信号を出力する送信信号処理部と、
変調送信信号に基づいて超音波を生体に送波して生体からの反射波を受波することにより受信信号を得る送受波部と、
実質的に変調送信信号に等しい波形の参照信号を用いて受信信号に対して復調処理を施すことにより復調信号を得る受信信号処理部と、
復調信号からドプラ情報を抽出するドプラ情報抽出部と、
を有し、
生体内の目標位置の深さに応じた遅延処理を施して参照信号と受信信号との間の遅延関係を調整し、目標位置からの受信信号と参照信号との相関を強めて復調処理を施すことにより、目標位置からのドプラ情報を選択的に抽出するにあたり、
前記ドプラ情報抽出部は、ドプラ情報として、直流成分と変調送信信号の変調周波数に関する基本波成分と２次以上の高調波成分とを含む復調信号から、直流成分に対応したドプラ信号を抽出する、
ことを特徴とする超音波診断装置。
【請求項２】
請求項１に記載の超音波診断装置において、
前記送信信号処理部から出力される変調送信信号と受信信号処理部で用いられる参照信号のうちの少なくとも一方の信号を遅延処理し、目標位置からの受信信号と参照信号との相関を強めることにより、目標位置からのドプラ情報を選択的に抽出する、
ことを特徴とする超音波診断装置。
【請求項３】
請求項２に記載の超音波診断装置において、
目標位置からの受信信号の位相と参照信号の位相が等しくなるように前記少なくとも一方の信号を遅延処理する、
ことを特徴とする超音波診断装置。
【請求項４】
請求項３に記載の超音波診断装置において、
目標位置からの受信信号の位相と参照信号の位相が等しくなるように、目標位置の深さに応じた遅延量だけ、前記変調送信信号または前記参照信号を遅延処理する、
ことを特徴とする超音波診断装置。
【請求項５】
請求項４に記載の超音波診断装置において、
前記遅延量を変化させることにより深さ方向に沿って目標位置を移動させつつ目標位置からのドプラ情報を抽出する、
ことを特徴とする超音波診断装置。
【請求項６】
請求項５に記載の超音波診断装置において、
前記遅延量を周期的に変化させ、深さ方向に沿った対象区間内において目標位置を周期的に移動させることにより、対象区間内に亘って複数の位置からドプラ情報を抽出する、
ことを特徴とする超音波診断装置。
【請求項７】
請求項６に記載の超音波診断装置において、
前記対象区間内に亘って複数の位置から抽出されるドプラ情報に基づいて、対象区間内における流体の速度分布を形成する、
ことを特徴とする超音波診断装置。
【請求項８】
請求項７に記載の超音波診断装置において、
生体内の深さに対応した軸と時間に対応した軸と流体の速度に対応した軸とを含み、前記速度分布の時間変化の様子を示した表示画像を形成する、
ことを特徴とする超音波診断装置。
【請求項９】
請求項１から８のいずれか１項に記載の超音波診断装置において、
前記ドプラ情報抽出部は、ドプラ情報として、復調信号に含まれる直流信号成分に対応したドプラ信号成分を抽出する、
ことを特徴とする超音波診断装置。
【請求項１０】
請求項１から９のいずれか１項に記載の超音波診断装置において、
前記送信信号処理部は、搬送波信号に対して周波数変調処理を施すことにより変調送信信号を生成し、
前記送信信号処理部による周波数変調処理の変調指数が調整されることにより位置分解能が設定される、
ことを特徴とする超音波診断装置。
【請求項１１】
請求項１０に記載の超音波診断装置において、
前記送信信号処理部で生成された変調送信信号を目標位置の深さに応じた遅延量だけ遅延処理することにより前記参照信号が形成される、
ことを特徴とする超音波診断装置。
【請求項１２】
請求項１０に記載の超音波診断装置において、
前記周波数変調処理で用いられる変調信号を目標位置の深さに応じた遅延量だけ遅延処理して得られる遅延変調信号を用いて前記搬送波信号を周波数変調処理することにより前記参照信号が形成される、
ことを特徴とする超音波診断装置。
【請求項１３】
請求項１０から１２のいずれか１項に記載の超音波診断装置において、
前記送信信号処理部による周波数変調処理の変調指数は、最大周波数偏移と変調周波数との比であり、その値が１以上に設定される、
ことを特徴とする超音波診断装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、超音波診断装置に関し、特に、変調された連続波を利用する超音波診断装置に関する。
【背景技術】
【０００２】
超音波診断装置の連続波を利用した技術として、連続波ドプラが知られている。連続波ドプラでは、例えば、数ＭＨｚの正弦波として構成される送信波が生体内へ連続的に放射され、生体内からの反射波が連続的に受波される。反射波には、生体内における運動体（例えば血流）によるドプラシフト情報が含まれる。そこで、そのドプラシフト情報を抽出して周波数解析することにより、運動体の速度情報を反映させたドプラ波形などを形成することができる。
【０００３】
連続波を利用した連続波ドプラは、パルス波を利用したパルスドプラに比べて一般に高速の速度計測の面で優れている。こうした事情などから、本願発明者は、連続波ドプラに関する研究を重ねてきた。その成果の一つとして、特許文献１において、周波数変調処理を施した連続波ドプラ（ＦＭＣＷドプラ）に関する技術を提案している。
【０００４】
一方、連続波ドプラでは、連続波を利用していることにより位置計測が困難である。例えば、従来の一般的な連続波ドプラの装置（ＦＭＣＷドプラを利用しない装置）では、位置計測を行うことができなかった。これに対し、本願発明者は、特許文献２において、ＦＭＣＷドプラにより生体内組織の速度に加えて生体内組織の位置を計測することができる技術を提案している。
【０００５】
【特許文献１】特開２００５−２５３９４９号公報
【特許文献２】特開２００６−１４９１６号公報
【発明の開示】
【発明が解決しようとする課題】
【０００６】
特許文献１や特許文献２に記載されたＦＭＣＷドプラの技術は、それまでにない超音波診断の可能性を秘めた画期的な技術である。本願発明者は、この画期的な技術の改良についてさらに研究を重ねてきた。
【０００７】
本発明は、このような背景において成されたものであり、その目的は、連続波を利用して目的の位置からドプラ情報を取得する技術を提供することにある。
【課題を解決するための手段】
【０００８】
上記目的を達成するために、本発明の好適な態様である超音波診断装置は、周波数を周期的に変化させる連続波の変調送信信号を出力する送信信号処理部と、変調送信信号に基づいて超音波を生体に送波して生体からの反射波を受波することにより受信信号を得る送受波部と、実質的に変調送信信号に等しい波形の参照信号を用いて受信信号に対して復調処理を施すことにより復調信号を得る受信信号処理部と、復調信号からドプラ情報を抽出するドプラ情報抽出部と、を有し、生体内の目標位置の深さに応じた遅延処理を施して参照信号と受信信号との間の遅延関係を調整して復調処理を施すことにより、目標位置からのドプラ情報を選択的に抽出することを特徴とする。
【０００９】
上記態様では、参照信号を用いて受信信号に対して復調処理を施しているため、参照信号との相関が比較的大きい信号成分を含んだ復調信号を得ることができる。そして、その復調処理にあたり、目標位置の深さに応じて参照信号と受信信号との間の遅延関係（時間軸方向の関係）が調整される。例えば参照信号と受信信号との間の位相関係が調整される。そのため、例えば、目標位置からの受信信号の位相と参照信号の位相とを揃えることにより、参照信号との相関が比較的大きい信号成分として、目標位置からの受信信号を抽出することができる。さらに、その受信信号から例えばバンドパスフィルタあるいはローパスフィルタなどを用いてドプラ情報を抽出することにより、目標位置からのドプラ情報を選択的に抽出することが可能になる。なお、上記態様において、参照信号と変調送信信号は、完全に等しい波形であることが望ましい。但し、参照信号と変調送信信号は、実質的に等しい波形とみなせる程度の対応関係であってもよい。
【００１０】
望ましい態様において、前記超音波診断装置は、前記送信信号処理部から出力される変調送信信号と受信信号処理部で用いられる参照信号のうちの少なくとも一方の信号を遅延処理し、目標位置からの受信信号と参照信号との相関を強めることにより、目標位置からのドプラ情報を選択的に抽出することを特徴とする。
【００１１】
望ましい態様において、前記超音波診断装置は、目標位置からの受信信号の位相と参照信号の位相が等しくなるように前記少なくとも一方の信号を遅延処理することを特徴とする。
【００１２】
望ましい態様において、前記超音波診断装置は、目標位置からの受信信号の位相と参照信号の位相が等しくなるように、目標位置の深さに応じた遅延量だけ、前記変調送信信号または前記参照信号を遅延処理することを特徴とする。
【００１３】
望ましい態様において、前記超音波診断装置は、前記遅延量を変化させることにより深さ方向に沿って目標位置を移動させつつ目標位置からのドプラ情報を抽出することを特徴とする。
【００１４】
望ましい態様において、前記超音波診断装置は、前記遅延量を周期的に変化させ、深さ方向に沿った対象区間内において目標位置を周期的に移動させることにより、対象区間内に亘って複数の位置からドプラ情報を抽出することを特徴とする。
【００１５】
望ましい態様において、前記超音波診断装置は、前記対象区間内に亘って複数の位置から抽出されるドプラ情報に基づいて、対象区間内における流体の速度分布を形成することを特徴とする。
【００１６】
望ましい態様において、前記超音波診断装置は、生体内の深さに対応した軸と時間に対応した軸と流体の速度に対応した軸とを含み、前記速度分布の時間変化の様子を示した表示画像を形成することを特徴とする。
【００１７】
望ましい態様において、前記ドプラ情報抽出部は、ドプラ情報として、復調信号に含まれる直流信号成分に対応したドプラ信号成分を抽出することを特徴とする。
【００１８】
望ましい態様において、前記送信信号処理部は、搬送波信号に対して周波数変調処理を施すことにより変調送信信号を生成し、前記送信信号処理部による周波数変調処理の変調指数が調整されることにより位置分解能が設定されることを特徴とする。望ましい態様において、前記送信信号処理部による周波数変調処理の変調指数は、最大周波数偏移と変調周波数との比でありその値が例えば１以上に設定される。なお、変調指数は２０以上に設定されてもよいし、３０以上に設定されてもよい。
【００１９】
望ましい態様において、前記送信信号処理部で生成された変調送信信号を目標位置の深さに応じた遅延量だけ遅延処理することにより前記参照信号が形成されることを特徴とする。望ましい態様において、前記周波数変調処理で用いられる変調信号を目標位置の深さに応じた遅延量だけ遅延処理して得られる遅延変調信号を用いて前記搬送波信号を周波数変調処理することにより前記参照信号が形成されることを特徴とする。
【発明の効果】
【００２０】
本発明により、連続波を利用して目的の位置からドプラ情報を取得することが可能になる。
【発明を実施するための最良の形態】
【００２１】
以下、本発明の好適な実施形態を図面に基づいて説明する。
【００２２】
図１には、本発明に係る超音波診断装置の好適な実施形態が示されており、図１はその全体構成を示す機能ブロック図である。送信用振動子１０は生体内へ送信波を連続的に送波し、また、受信用振動子１２は生体内からの反射波を連続的に受波する。このように、送信および受信がそれぞれ異なる振動子で行われて、いわゆる連続波ドプラ法による送受信が実行される。
【００２３】
電力増幅器１４は、送信用振動子１０に対し、電力増幅されたＦＭ連続波（ＦＭＣＷ）を供給する。電力増幅器１４には、例えば正弦波によるＦＭ変調処理が施されたＦＭ連続波（ＦＭＣＷ波）が入力され、このＦＭ連続波に対応する送信波が送信用振動子１０から送波される。ＦＭ変調器２０は、ＦＭ連続波を電力増幅器１４に出力する。ＦＭ変調器２０は、ＲＦ波発振器２２から供給されるＲＦ波、および、ＦＭ変調波発振器２４から供給される正弦波の変調波に基づいてＦＭ連続波を発生する。このＦＭ連続波の波形については後の原理説明で詳述する。
【００２４】
前置増幅器１６は、受信用振動子１２から供給される受波信号に対して低雑音増幅等の受信処理を施し、受信ＲＦ信号を形成して受信ミキサ３０へ出力する。受信ミキサ３０は受信ＲＦ信号に対して直交検波を施して複素ベースバンド信号を生成する回路であり、２つのミキサ３２，３４で構成される。各ミキサは受信ＲＦ信号を所定の参照信号と混合する回路である。
【００２５】
受信ミキサ３０の各ミキサに供給される参照信号は、ＦＭ変調器２０から出力されＦＭ連続波に基づいて生成される。つまり、ＦＭ変調器２０から出力されるＦＭ連続波が遅延回路２５において遅延処理され、ミキサ３２には遅延処理されたＦＭ連続波が直接供給され、一方、ミキサ３４には遅延処理されたＦＭ連続波がπ／２シフト回路２６を経由して供給される。π／２シフト回路２６は遅延処理されたＦＭ連続波の位相をπ／２だけずらす回路である。この結果、２つのミキサ３２，３４の一方から同相信号成分（Ｉ信号成分）が出力され、他方から直交信号成分（Ｑ信号成分）が出力される。そして、受信ミキサ３０の後段に設けられるＬＰＦ（ローパスフィルタ）３６，３８によって、同相信号成分および直交信号成分の各々の高周波数成分がカットされ、検波後の必要な帯域のみの復調信号が抽出される。
【００２６】
なお、ＦＭ変調波発振器２４から供給される正弦波の変調波に遅延処理を施して遅延変調波を形成し、その遅延変調波によってＲＦ波発振器２２から供給されるＲＦ波を周波数変調することにより、参照信号を生成してもよい。
【００２７】
後の原理説明で詳述するが、各ミキサで実行される受信ＲＦ信号と参照信号との混合処理の結果である受信ミキサ出力信号（復調信号）には、ＦＭ変調波発振器２４から供給される変調波の変調波周波数ｆｍに関する複数の第ｎ次波成分（ｎは０以上の自然数）が含まれている。つまり、第０次波成分である直流成分、第１次波成分である基本波成分、さらに、ｎが２以上の複数の高調波成分が含まれている。つまり、これら複数の第ｎ次波成分を含んだ復調信号が、ＬＰＦ３６，３８の各々から出力される。
【００２８】
ＦＦＴ回路（高速フーリエ変換回路）４０，４２は、復調信号（同相信号成分および直交信号成分）の各々に対してＦＦＴ演算を実行する。その結果、ＦＦＴ回路４０，４２において復調信号が周波数スペクトラムに変換される。なお、ＦＦＴ回路４０，４２から出力される周波数スペクトラムは、回路の設定条件などにより周波数分解能δｆの周波数スペクトラムデータとして出力される。ＦＦＴ回路４０，４２から出力される周波数スペクトラムについては、後に図２などを利用して詳述する。
【００２９】
ドプラ情報解析部４４は、周波数スペクトラムに変換された復調信号からドプラ情報を抽出する。その際、予め遅延回路２５によって生体内の目標位置の深さに応じて参照信号と受信信号との間の位相関係が調整されているため、目標位置からのドプラ情報が選択的に抽出される。位相の調整と目標位置からのドプラ情報の抽出との関連については、後の原理説明において詳述する。ドプラ情報解析部４４は、生体内の各深さ（各位置）ごとにドプラ情報を抽出して、例えば、超音波ビーム（音線）上の各深さごとに生体内組織の速度を算出し、リアルタイムで出力する。なお、超音波ビームを走査させて二次元的あるいは三次元的に生体内組織の各位置の速度を算出してもよい。
【００３０】
表示処理部４６は、生体組織の深さ（位置）ごとの速度に基づいて、例えばドプラ波形や、深さ速度の情報を含むグラフなどを形成し、形成したドプラ波形やグラフなどを表示部４８にリアルタイムで表示させる。なお、図１に示す超音波診断装置内の各部は、システム制御部５０によって制御される。つまり、システム制御部５０は、送信制御や受信制御や表示制御などを行う。
【００３１】
以上、概説したように、本実施形態では、連続波（ＣＷ）を変調波でＦＭ変調した超音波（ＦＭＣＷ波）を送受波して受信信号が得られて、生体内の目標位置の深さに応じて参照信号と受信信号との間の位相関係が調整されて復調処理されることにより、目標位置からのドプラ情報が選択的に抽出される。そこで、目標位置からのドプラ情報が選択的に抽出される原理について詳述する。
【００３２】
周波数ｆ０のＣＷに変調周波数ｆｍの正弦波によりＦＭ変調を施したＦＭＣＷ送信波は次式のように表現できる。
【００３３】
【数１】
0005459963.tif 000002
【００３４】
数１において、Δｆは周波数変動幅の０−Ｐ値（ゼロピーク値：最大周波数偏移）であり、最大周波数偏移Δｆと変調周波数ｆｍの比であるβはＦＭの変調指数である。また、ドプラシフトを伴わない場合のＦＭＣＷ受信波は生体による減衰を無視すると次式で表現できる。
【００３５】
【数２】
0005459963.tif 000003
【００３６】
ＦＭＣＷ送信波の周波数スペクトラムは、数１をベッセル級数を用いて展開することで得られる。数１に示すＦＭＣＷ送信波は次式のように展開できる。
【００３７】
【数３】
0005459963.tif 000004
【００３８】
数３において、J０（β），J２ｎ（β），J２ｎ＋１（β）は、第1種ベッセル関数である。各項の振幅は、変調指数βおよびそれに対応するベッセル関数によって決定される。
【００３９】
また、ドプラシフトを伴わない場合の受信波ｖR（ｔ）の周波数スペクトラムは、数２を展開することで得られる。数２に示すＦＭＣＷ受信波は次式のように展開できる。
【００４０】
【数４】
0005459963.tif 000005
【００４１】
数４に示されるように、受信波の周波数スペクトラムは送信波と同じ周波数成分を持っている。しかし受信波の各周波数成分の振幅は、位相差φ０とφｍに応じて変化している。
【００４２】
さらに、ドプラシフトを伴う場合、数２のｖR（ｔ）は以下のように書き換えられる。  
【００４３】
【数５】
0005459963.tif 000006
【００４４】
なお、数５においてｆｍに対するドプラシフトは、ｆ０のシフト分ｆｄに比較して小さいので無視している。
【００４５】
上述の数２や数５で表される受信波形は、超音波振動子が受信する信号波形（受信ＲＦ信号）である。超音波診断装置は受信ＲＦ信号に対して復調処理を実行する。ＦＭＣＷの受信ＲＦ信号を復調する場合、復調系ではＦＭＣＷ送信波を参照信号として受信波と乗算を行う。復調系における受信ミキサ出力は、ｖT（ｔ）とｖR（ｔ）を乗算した結果として次式のように算出される。
【００４６】
【数６】
0005459963.tif 000007
【００４７】
ここで、数３、数４、数６の算出にはベッセル関数に関する次の公式を利用する。
【００４８】
【数７】
0005459963.tif 000008
【００４９】
数７の公式を用いると、数６はさらに次式のように計算される。なお、数８では数６における係数１／２を省略する。
【００５０】
【数８】
0005459963.tif 000009
【００５１】
数８で表現される受信ミキサ出力の周波数スペクトラム、つまり、図１のＦＦＴ回路４０，４２から出力される周波数スペクトラムを図示すると図２のようになる。
【００５２】
図２は、復調信号の周波数スペクトラムを示す模式図である。数８にも示されるように、復調信号には、変調波周波数に関する複数の第ｎ次高調波成分（ｎは０以上の自然数）が含まれている。つまり、図２において、原点Ｏ付近に存在する第０次波成分である直流成分、第１次波成分である基本波成分（ｆｍ）、第２次波成分である第２高調波成分（２ｆｍ）、第３次波成分である第３高調波成分（３ｆｍ）が含まれている。なお、図示省略しているｎが７以上の高調波成分も存在する。そして、第ｎ次波成分の各々には、固定物エコー６４とドプラエコー６２が含まれている。
【００５３】
固定物エコー６４は、生体内の強反射体である静止体からのエコー（クラッタエコー）であり、ドプラ情報を観測する場合においては妨害要因となる。これに対し、ドプラエコー６２が必要としているドプラ信号である。ドプラシフトを伴う場合、ドプラ信号の第ｎ次波成分の各々は、ＦＭ変調周波数が抑圧されたＤＳＢ−ＳＣ（Double  Sideband−Suppressed  Carrier）の形で周波数スペクトル上に出現する。
【００５４】
本実施形態は、上述したＦＭＣＷの基本原理を応用したものである。本実施形態では、図１に示す遅延回路２５（移相器あるいは遅延線）により、生体内の目標位置の深さに応じて参照信号と受信信号との間の位相関係が調整され、目標位置（目標となる深さ）からのドプラ情報が選択的に抽出される。つまり、遅延回路２５による遅延処理の遅延時間を、超音波が生体内を往復する伝搬時間に設定することにより、目標となる深さからの受信信号と参照信号との相関が最大となり目標からの信号のみが選択的に抽出される。
【００５５】
図３は、参照信号と受信信号との間の位相関係（時間関係）を説明するための図である。（ａ）はＦＭＣＷ送信信号（送信波）の時間波形を示しており、（ｂ）はその送信信号に対応した受信信号（受信波）を示している。受信信号は送信時からτだけ遅れて受信される。音速をｃとすると、深さｄから得られる受信信号はτ＝２ｄ／ｃだけ遅れて受信される。また、図３において、ｆ０は搬送波周波数（超音波周波数に相当）を示しており、ｆｍは搬送波を変調する周波数（変調周波数）を示している。受信信号の周波数スペクトラム電力を（ｅ）に示す。ちなみに、送信信号の周波数スペクトラムも電力として表現すると、（ｅ）と同じ形をしているが、各電力スペクトラムの位相は遅延時間τに対応して、送信スペクトラムとは異なっている。
【００５６】
受信信号と直交検波する参照波（参照信号）（ｃ）は、送信波を送受信時間差（往復の伝搬時間τ）だけシフトさせた信号である。直交検波器入力、つまり受信ミキサ（図１の符号３０）に入力される受信信号と参照信号は、このシフト操作により周波数、位相ともまったく等しい時間波形となる。したがって、受信波（受信信号）を二乗した波形が直交検波器出力（ｄ）となる。以降、直交検波器出力信号をベースバンド信号と呼ぶ。この信号は、振幅が１／２の直流成分と、搬送波周波数が２ｆ０で変調度が２βのＦＭ信号の和として表現できる。受信波を二乗することにより、受信波が、直流付近と搬送波の２倍の周波数帯に移行するという特徴がある。（ｆ）におけるベースバンド信号の周波数スペクトラムは、この様子を表している。
【００５７】
ベースバンド信号は、受信波と参照波間の遅延時間が完全に一致している場合は、ｆｍおよびその高調波成分（図２参照）は出現せず、図３（ｆ）のように直流成分のみとなる。なお、図３（ｆ）における２ｆ０や−２ｆ０の成分は、例えば受信ミキサの後段のＬＰＦ（図１の符号３６，３８）によって除去される。受信波と参照波間の遅延時間が一致していないと、両者の時間差により、変調波の高調波成分、つまりｆｍおよびその高調波成分（図２参照）が発生する。高調波成分は、受信波と参照波間の時間差が“０”からわずかでもずれると発生する。
【００５８】
こうした特徴から、本実施形態のＦＭＣＷ方式では、受信波（受信信号）と参照波（参照信号）との間の位相関係を調整することにより、ベースバンド信号の直流および直流付近の信号成分に基づいて対象とする組織の速度情報を得ることができる。この意味において、本実施形態に係る超音波診断装置を位相シフト型ＦＭＣＷ超音波ドプラシステムと称することができる。
【００５９】
図４は、ベースバンド信号の直流および高調波成分の生体内深さ（体表からの距離）依存性を説明するための図である。図４には、複数の変調指数βについて、各変調指数βごとにベースバンド信号（復調信号）の周波数スペクトラムが示されている。図４に示す各変調指数βごとの周波数スペクトラムは、図２の周波数スペクトラムの表示に距離方向の次元を追加したものに相当する。なお、β＝０からβ≒３０までの各周波数スペクトラムは、固定組織からからの反射電力を表している。図４においては、組織における減衰の効果は省略した。なお、ＦＭ変調度の大きさは、一般的に、変調指数βにより定量化される。βは、ＦＭ変調による搬送波の最大周波数偏移Δｆと変調周波数ｆｍの比として定義され、β＝Δｆ／ｆｍにより定義される。
【００６０】
β＝０の場合は無変調であるため、これは通常のＣＷドプラ速度計測システムと等価である。この場合は、どんな深さからの反射電力にも位置依存性は無い。また、送受信ともに無変調なので、ベースバンド信号に変調波成分が出現する余地は無い。連続超音波（ＣＷ）にＦＭ変調をかけ、周波数偏移Δｆを徐々に増加させると、ＦＭＣＷ送受信波はＦＭ変調された信号となり、その電力は搬送波から側帯波に移行してゆく。送受信間の遅延時間差が無い場合はベースバンド信号に変調波の高調波成分は発生しない（図３参照）。直交検波器の２つの入力（つまり受信信号と参照信号）に時間差を生じさせないためには、当該深さからの反射波の遅延時間に相当する遅延時間を参照波に与えてやればよい。
【００６１】
図４の例は、距離ｄ０＝７．５ｃｍにおいて、参照波に与える遅延時間を受信波の遅延時間と一致するように設定した場合を示している。したがって、距離ｄ０からの反射電力は、直流成分のみで、変調波の高調波成分はまったく発生しない。距離がｄ０以外の場所では、送信波と受信波間に時間差が生じるので、両者の相関性が少なくなってくる。したがって、直流成分が減少し、同時に高調波成分が出現しはじめる。この現象は、図３の（ｂ），（ｃ）間に時間差が生じたときに、瞬時周波数の差が、変調周波数の周期で繰り返す現象が発生するという概念からも説明できる。
【００６２】
変調度（変調指数β）が大きくなるにつれて、距離ｄ０における反射電力の位置ずれに対する変化は敏感になってくる。すなわち、距離ｄ０における反射電力の選択性が増加する。この傾向はβが増加するにしたがってますます顕著になる。図４では、その様子をβが０から３０の場合について示している。そして、βが３０以上になると、この選択性はＰＷ（パルス波）ドプラにおけるレンジゲートの役割と類似の機能に近づいてくる。
【００６３】
次に、距離ｄ０からの反射電力が、ドプラシフトｆｄを伴っていると仮定する。この場合は、固定物と相似形の距離依存性が出現するが、ドプラ周波数ｆｄだけ、直流成分からシフトしてあらわれる。ドプラスペクトラムは、変調波あるいはその高調波の両側帯波にも同時に出現するが、距離ｄ０からのドプラエコーは直流からｆｄだけシフトした周波数成分のみとなる。この様子を図４（ＤＰ）に示す。
【００６４】
図４（ＤＰ）に示すように、距離ｄ０からのドプラエコーが直流からｆｄだけシフトした周波数成分のみとなるのは、搬送波周波数ｆ０からｆｄだけずれた周波数において、参照信号とドプラ信号との相関が最も強くなるからである。このドプラ信号は、直流付近にのみ出現し、変調波の高調波成分（ｆｍ，２ｆｍ，３ｆｍ，・・・）の付近には出現しない。したがって、直流付近のこの成分だけをローパスフィルタ（低域濾波器）によって抽出することにより、通常のＣＷと同様のＳＮＲ（信号体ノイズ比）を保った状態で、位置情報の特定されたドプラ情報を得ることができる。
【００６５】
具体的には、例えば、図１の遅延回路２５によって、目標となる深さ（位置）ｄに対応した遅延時間τ（τ＝２ｄ／ｃ：音速ｃ，深さｄ）で参照信号に対して遅延処理が施されることにより、図１のＦＦＴ回路４０，４２から図４（ＤＰ）に相当するドプラ周波数スペクトラム情報が出力され、図１のドプラ情報解析部４４によって、図４（ＤＰ）の直流付近のドプラ信号が抽出される。
【００６６】
次に、本実施形態に係る超音波診断装置（位相シフト型ＦＭＣＷ超音波ドプラシステム）に関するいくつかの特徴について説明する。
【００６７】
「クラッタ電力」
深さ（位置）情報の特定の程度は変調度（変調指数）βに依存するので、必要な位置分解能を確保するためには、βを大きく設定する必要がある。βを大きく設定すると、位置の選択性がシャープになる。つまり、図４に示したようにβが大きくなると目標の距離ｄ０において急峻な山の波形が得られる。そのため、距離ｄ０付近の反射電力だけが受信信号に大きく反映される。この性質は、固定目標でも移動目標でも共通である。
【００６８】
音線上のすべての固定目標からの反射電力の和がクラッタ電力となる従来のＣＷ法に比べて、本実施形態の方式（位相シフト型ＦＭＣＷ）では、選択された位置からのクラッタしか生じないので、クラッタ電力を極端に減少させることができる。
【００６９】
図５は、ドプラエコーが発生する深さ（体表からの深さ）とクラッタ積算電力の関係を示す計算結果の図である。パラメータはＦＭ変調度（β）である。図５には、各ＦＭ変調度（β）ごとに、クラッタ積算電力８４の波形とドプラ信号電力８２の波形が示されている。また、クラッタ改善度９０は、従来のＣＷと比較した場合のクラッタ積算電力の低減を示している。変調度（変調指数）βが増加するに従って、クラッタ電力（積算電力）が全体的に減少する傾向が示されている。例えばβ＝１００に設定すると、従来のＣＷに比べてクラッタ電力を２５ｄＢ程度低減できる。
【００７０】
「信号対雑音比（ＳＮＲ）」
図６は、ＰＷ（パルス波）と通常のＣＷとＦＭＣＷの各々のドプラ速度計測システムについてのＲＦ帯とベースバンド帯における受信信号のスペクトルを示す図である。ＰＷドプラのＳＮＲを支配している雑音は、直交検波後のサンプリング時における信号帯域により決定される。この雑音は、通信システムと同様に、ＲＦ信号の帯域に分布している白色ガウス雑音として扱うことができる。サンプリング時の信号帯域は、ベースバンドにおけるパルス幅の逆数で与えられる。たとえば、パルス幅を２μｓｅｃ．と仮定すると、帯域は５００ｋＨｚとなる。
【００７１】
一方、通常のＣＷドプラ速度計測システムでは、ＰＲＦが存在せず、信号の帯域は測定するドプラ信号の最高周波数程度である。これをたとえば５ＫＨｚとすると、ＰＷとＣＷのＳＮＲの差は帯域比の差で表現されるので、この例では、５００ｋＨｚ／５ｋＨｚ＝１００となる。すなわち、ＣＷのＳＮＲはＰＷに比べて２０ｄＢ良好となる。
【００７２】
ＦＭＣＷの信号帯域は、ベースバンド信号におけるドプラ信号の最高周波数により決定される。ただし、ＦＭＣＷでは、ＲＦ（超音波）周波数帯では、変調度に依存した周波数帯域を持っており、信号は広い帯域に拡散している。しかし、直交検波により、この広帯域信号を直流付近の狭帯域ベースバンド信号に圧縮変換するので、雑音も狭帯域となり、ＳＮＲが大幅に改善される。改善の程度は、ＣＷと同様で、５００ｋＨｚ／５ｋＨｚ＝１００の場合は、ＰＷに比べて２０ｄＢの改善が予測される。
【００７３】
「位置の選択性（位置分解能）」
ＰＷドプラでは、特定位置からのドプラ情報を得られるが、通常のＣＷドプラ方式では位置情報が得られない。この特徴を相補するのが本実施形態の位相シフト型ＦＭＣＷということができる。本実施形態の復調過程では、直交検波用参照信号の遅延時間を目標までの距離に整合させることにより位置の選択性を持たせ、この選択性は変調度βを大きくすることにより向上するという性質がある。変調度βを大きくすると、搬送波帯の帯域が広がる。必要とする占有帯域幅ＢＷは近似的に下式で与えられる。
【００７４】
【数９】
0005459963.tif 000010
【００７５】
図７は、変調度に応じた位置選択性の計算結果を示す図である。図７に示す半値幅は、目標となる位置のスペクトル波形の半値幅（図４に示す各周波数スペクトラムの距離ｄ０におけるスペクトル波形の山の半値幅）である。
【００７６】
図７に示すように、βを大きくすると、その位置選択性を数ｍｍ以下に設定することができる。この結果は、搬送波の周波数帯を広帯域にするほど、目標の位置分解能が向上することを意味している。これは、ＰＷにおいてパルス幅をせまくするほど位置分解能が向上するという特性とも矛盾しない。
【００７７】
つまり、本実施形態の位相シフト型ＦＭＣＷ方式では、位置分解能を満足する程度までβを大きく設定すれば、すなわち、数９で与えられるＦＭＣＷの帯域をＰＷと同程度に設定すれば、位置情報を持ち且つＳＮＲが良好なドプラ情報を得ることが可能となる。
【００７８】
「最高血流速度」
ＦＭＣＷ方式は、ＰＷと同様に、最高速度に制限がある。ＦＭＣＷ方式では、周波数ｆｍの変調波で搬送波を変調しており、折り返し現象により、最高測定可能周波数は、ｆｍ／２に制限される。この制限を低減するためにはｆｍを大きく設定すればよい。しかし、この場合は、選択した位置のほかに、ｆｍの周期に相当する位置からのドプラ信号が本来の位置からの信号に重畳される。この性質は、ＰＷのＨＰＲＦモードの場合とまったく同様である。本実施形態の位相シフト型ＦＭＣＷドプラ方式によれば、計測可能な最高血流速度に関して、ＰＷ方式の場合と同程度の計測能力が実現可能である。
【００７９】
「送受信部の構成」
ＰＷ、ＣＷ、ＦＭＣＷドプラ速度計測システムで、探触子入力における平均電力を等しいと仮定する。この場合、前述したようなクラッタおよび雑音に対する改善効果を達成するためのピーク電力は、ＰＷに比べて、本実施形態のＦＭＣＷ方式では、１／１００程度、ピーク電圧は１／１０に減少できる。一方、探触子に入力できる送信電力がピーク値により制限されていると仮定する。この場合、本実施形態のＦＭＣＷは、ＰＷに比べて、クラッタおよび雑音に対する改善効果に加え、電力を増加できる比率で、さらにクラッタ抑圧とＳＮＲの改善が期待できる。
【００８０】
以上、本発明の好適な実施形態を説明したが、本実施形態のいくつかの特徴（効果）をまとめると図８のようになる。
【００８１】
また、上述した実施形態では、変調波として、図９に示すような正弦波を利用した。この正弦波に換えて、変調波として図９に示す対称三角波を利用してもよい。対称三角波を用いると、位置情報とドプラ情報を切り分けて計測できる利点がある（特許文献２参照）。なお、対称三角波は、正弦波と似た周期性を有しているため、正弦波に換えて対称三角波を用いても、位相シフト型ＦＭＣＷ超音波ドプラシステムを構築できる。
【００８２】
さらに、図９に示すように、変調周波数変化の変換点における周波数変化を時間的にゆるやかにした対称三角波（角を丸めた対称三角波）を変調波として利用してもよい。この波形によれば、対称三角波が変調周波数変化の変換点で発生する高周波成分を除去できる特長がある。すなわち、この変換点における周波数変化を時間的にゆるやかにすることで、高周波成分の発生を低減し、ＲＦ帯域を必要以上に広げない効果が期待できる。
【００８３】
また、受信信号と参照信号の位相を調整するための変形態様として、図１において、遅延回路２５を電力増幅器１４の直前（図１のＰの位置）に移してもよい。つまり、電力増幅器１４へ供給される送信用のＦＭＣＷに遅延処理を施し、受信ミキサ３０へ供給される参照信号に遅延処理を施さずに、参照信号と受信信号の位相を調整してもよい。なお、図１に示す遅延回路２５をそのまま残し、さらに電力増幅器１４の直前にも遅延回路２５を設けてもよい。つまり、電力増幅器１４へ供給される送信用のＦＭＣＷに遅延処理を施し、さらに受信ミキサ３０へ供給される参照信号にも遅延処理を施して、参照信号と受信信号の位相を調整してもよい。
【００８４】
また、周波数変調処理で用いられる変調信号（変調波）を目標位置の深さに応じた遅延量だけ遅延処理して遅延変調信号を生成し、そして、遅延変調信号を用いて搬送波信号（ＲＦ波）を周波数変調処理することにより参照信号を形成して、受信信号と参照信号の位相を調整してもよい。なお、この場合には、例えば、遅延処理されていない変調信号を用いて搬送波信号を周波数変調処理することより変調送信信号（ＦＭＣＷ送信波）が形成される。
【００８５】
図１に示した超音波診断装置では、遅延回路２５によって生体内の目標位置の深さに応じて参照信号と受信信号との間の位相関係が調整されているため、目標位置からのドプラ情報が選択的に抽出される。そのため、生体内の各深さ（各位置）ごとにドプラ情報を抽出して、例えば、超音波ビーム（音線）上の各深さごとに生体内組織の速度を算出することができる。さらに、遅延回路２５における遅延量を周期的に変化させることにより、以下に説明するように、超音波ビームに沿った対象区間内において目標位置を周期的に移動させ、対象区間内に亘って複数の位置からドプラ情報を抽出することも可能になる。
【００８６】
図１０には、本発明に係る超音波診断装置の別の好適な実施形態が示されており、図１０はその全体構成を示す機能ブロック図である。図１０に示す超音波診断装置は、図１の装置を改良したものであり、第２変調波発振器５２が挿入されている点において図１の装置とは異なる。以下、図１０の超音波診断装置について、第２変調波発振器５２の挿入に伴う改良点を中心に説明し、図１の装置と同じ部分については説明を省略する。
【００８７】
図１０の超音波診断装置においても、図１の場合と同様に、遅延回路２５によって生体内の目標位置の深さに応じて参照信号と受信信号との間の位相関係が調整され、目標位置からのドプラ情報が選択的に抽出される。つまり、遅延回路２５による遅延処理の遅延時間（遅延量）を、超音波が生体内を往復する伝搬時間に設定することにより、目標となる深さからの受信信号と参照信号との相関が最大となり目標からの信号のみが選択的に抽出される。
【００８８】
図１０においては、遅延回路２５における遅延時間を周期的に変化させ、深さ方向に沿った対象区間内において目標位置を周期的に移動させることにより、対象区間内に亘って複数の位置からドプラ情報を抽出する。遅延回路２５は、第２変調波発振器５２が出力する信号に基づいて、遅延時間を周期的に変化させる。
【００８９】
第２変調波発振器５２は、ＦＭ変調波発振器２４から出力される変調波（第１変調波）よりも低い周波数の第２変調波を出力する。第１変調波の周波数が、例えば５ｋＨｚ程度であれば、第２変調波の周波数が、例えば５０Ｈｚ程度に設定される。遅延回路２５は、この比較的低い周波数の第２変調波を利用して遅延時間を周期的に変化させる。
【００９０】
遅延時間が周期的に変化することにより、参照信号との相関が最大となる目標位置が超音波ビーム方向（生体の深さ方向）に沿って周期的に変化する。つまり、深さ方向に沿ったある範囲（対象区間）内において目標位置が周期的に移動し、対象区間内に亘る複数の位置からドプラ情報を抽出することができる。
【００９１】
なお、図１の場合と同様に、図１０においても、遅延回路２５を電力増幅器１４の直前（図１０のＰの位置）に移してもよい。つまり、電力増幅器１４へ供給される送信用のＦＭＣＷに周期的な変化の遅延処理を施し、受信ミキサ３０へ供給される参照信号に遅延処理を施さずに、参照信号と受信信号の位相を調整してもよい。
【００９２】
周期的に移動する目標位置は、遅延回路２５における遅延時間（遅延量）に応じて決定される。システム制御部５０は、第２変調波発振器５２が出力する第２変調波に基づいて、あるいは、遅延回路２５における遅延時間に基づいて、周期的に移動する目標位置を確認する。そして、システム制御部５０は、ドプラ情報解析部４４において得られるドプラ情報（ドプラシフト量やドプラ成分の電力など）と、そのドプラ情報が得られた位置（深さ）とを対応付ける。その対応関係は、表示処理部４６などに伝えられる。
【００９３】
図１１は、第２変調波が正弦波の場合のドプラ情報と深さの対応関係を説明するための図である。図１１には、底面に時間軸と深さ軸を設けて高さ方向のＺ軸にドプラ情報（ドプラシフト量またはドプラ成分の電力）を示したグラフが図示されている。
【００９４】
第２変調波が正弦波の場合、時間の経過と共に遅延回路２５（図１０）における遅延時間が正弦波状に変化し、ドプラ情報が得られる目標位置（深さ）も時間の経過と共に正弦波状に変化する。
【００９５】
図１１には、時間の経過と共に、目標位置が深さ−ｄ／２から＋ｄ／２の間（対象区間）を正弦波に沿って変化する様子が示されている。なお、図１１に示すグラフは、第２変調波の周波数がｆｍ２の場合であり、目標位置が周期１／ｆｍ２の正弦波に沿って変化している。そして、正弦波に沿って変化する各位置において、その位置から得られるドプラ情報がＺ軸方向に示されている。
【００９６】
このように、深さ−ｄ／２から＋ｄ／２までの対象区間内において、各位置（深さ）ごとにドプラ情報が得られるため、例えば、そのドプラ情報から各位置（深さ）ごとに速度を算出することにより、図１１に示すように、深さ−ｄ／２から＋ｄ／２までの対象区間内における速度分布７０を得ることができる。
【００９７】
なお、図１１においては、深さ−ｄ／２と＋ｄ／２の中間位置において、比較的大きなドプラ情報が得られている。血管内においては、通常、血管壁の近傍に比べて血管の中心付近において大きな血流が得られる。そのため、例えば、血管と交差するように超音波ビームを設定し、超音波ビーム上の血管に対応した部分に対象区間（深さ−ｄ／２から＋ｄ／２）を設定することにより、例えば、図１１のような計測結果が得られ、血管内の血流の速度分布７０が得られる。
【００９８】
図１２は、速度分布の表示態様を示す図であり、図１２には、底面に時間軸と深さ軸を設けて、ドプラ情報から得られる流速を高さ方向に示したグラフが図示されている。図１２に示すグラフは、図１０の表示処理部４６において形成されて表示部４８に表示される画像の一例である。表示処理部４６は、例えば、システム制御部５０から得られるドプラ情報とその位置（深さ）との対応関係に基づいて図１２に示す画像を形成する。
【００９９】
つまり、時間の経過と共に移動する目標位置（深さ）と、図１０のドプラ情報解析部４４から得られるドプラ情報が対応付けられ、ドプラ情報から速度（流速）が算出される。そして、時間の経過と共に、目標位置が深さ−ｄ／２から＋ｄ／２までの対象区間内を繰り返し移動することにより、対象区間内における流速の時間変化を示す図１２のグラフが形成される。
【０１００】
図１２のグラフは、深さ−ｄ／２から＋ｄ／２までの対象区間内における流速分布が、時間の経過と共に変化する様子を示している。例えば、血管と交差するように超音波ビームを設定し、超音波ビーム上の血管に対応した部分に対象区間（深さ−ｄ／２から＋ｄ／２）を設定することにより、その血管の部分における血流の流速分布の変化として図１２に示すような計測結果が得られる。
【０１０１】
なお、図１１および図１２においては、第２変調波が正弦波の場合について説明したが、第２変調波は鋸波や三角波などでもよい。
【０１０２】
図１３は、第２変調波が鋸波の場合のドプラ情報と深さの対応関係を説明するための図である。図１３には、図１１と同様に、底面に時間軸と深さ軸を設けて高さ方向にドプラ情報を示したグラフが示されている。第２変調波が鋸波の場合、時間の経過と共に遅延回路２５（図１０）における遅延時間が鋸波状に変化し、ドプラ情報が得られる目標位置（深さ）も時間の経過と共に鋸波状に変化する。
【０１０３】
図１３には、時間の経過と共に、目標位置が深さ−ｄ／２から＋ｄ／２の間（対象区間）を鋸波に沿って変化する様子が示されている。なお、図１３に示すグラフにおいても、第２変調波の周波数がｆｍ２であり、目標位置が周期１／ｆｍ２の鋸波に沿って変化している。そして、鋸波に沿って変化する各位置において、その位置から得られるドプラ情報が高さ方向に示されている。また、図１１の場合と同様に、深さ−ｄ／２から＋ｄ／２までの対象区間内における速度分布７０を得るようにしてもよい。
【０１０４】
図１４は、第２変調波が対称三角波の場合のドプラ情報と深さの対応関係を説明するための図である。図１４には、図１１，図１３と同様に、底面に時間軸と深さ軸を設けて高さ方向にドプラ情報を示したグラフが示されている。第２変調波が対称三角波の場合、ドプラ情報が得られる目標位置（深さ）が時間の経過と共に対称三角波に沿って変化する。図１４には、対称三角波に沿って変化する各位置において、その位置から得られるドプラ情報が高さ方向に示されている。なお、図１１，図１３の場合と同様に、深さ−ｄ／２から＋ｄ／２までの対象区間内における速度分布７０を得るようにしてもよい。
【０１０５】
以上、本発明の好適な実施形態といくつかの変形形態を説明したが、上述した本発明の好適な実施形態等は、あらゆる点で単なる例示にすぎず、本発明の範囲を限定するものではない。本発明は、その本質を逸脱しない範囲で各種の変形形態を包含する。
【０１０６】
例えば、上述した実施形態においては、周波数を周期的に変化させる連続波の変調送信信号を形成するにあたって、搬送波信号（ＲＦ波発振器２２から供給されるＲＦ波）に対して周波数変調処理を施している。この周波数変調処理に換えて、周波数変調処理と同じ角度変調の方式として当業者において明らかな位相変調処理（ＰＭ処理）を利用してもよい。つまり、搬送波信号（ＲＦ波発振器２２から供給されるＲＦ波）に対して位相変調処理を施すことにより、ＦＭ変調器２０から出力されるＦＭ連続波と同じ波形あるいは同等な波形を形成してもよい。なお、周波数を周期的に変化させる連続波のデータをメモリなどに記憶しておき、このメモリから読み出されるデータに基づいて、当該連続波を生成してもよい。
【図面の簡単な説明】
【０１０７】
【図１】本発明に係る超音波診断装置の全体構成を示す機能ブロック図である。
【図２】復調信号（ベースバンド信号）のドプラエコーと固定物エコーの周波数スペクトラムを示す図である。
【図３】参照信号と受信信号との間の位相関係およびベースバンド信号の周波数スペクトラムを示す図である。
【図４】ベースバンド信号電力の深さ依存性を説明するための図である。
【図５】ドプラエコーが発生する深さとクラッタ積算電力の関係を示す計算結果の図である。
【図６】各方式についての受信信号と雑音の周波数スペクトルを示す図である。
【図７】変調度に応じた位置選択性の計算結果を示す図である。
【図８】本実施形態のいくつかの特徴を示す図である。
【図９】変調波の波形を説明するための図である。
【図１０】本発明に係る超音波診断装置の別の好適な実施形態を示す図である。
【図１１】第２変調波が正弦波の場合のドプラ情報と深さと時間の関係を説明するための図である。
【図１２】速度分布の表示態様を示す図である。
【図１３】第２変調波が鋸波の場合のドプラ情報と深さの対応関係を説明するための図である。
【図１４】第２変調波が対称三角波の場合のドプラ情報と深さの対応関係を説明するための図である。
【符号の説明】
【０１０８】
２０  ＦＭ変調器、２２  ＲＦ波発振器、２４  ＦＭ変調波発振器、２５  遅延回路、４０，４２  ＦＦＴ回路、４４  ドプラ情報解析部、５２  第２変調波発振器。
【図１】
0005459963.tif 000011
【図２】
0005459963.tif 000012
【図３】
0005459963.tif 000013
【図６】
0005459963.tif 000014
【図７】
0005459963.tif 000015
【図８】
0005459963.tif 000016
【図９】
0005459963.tif 000017
【図１０】
0005459963.tif 000018
【図４】
0005459963.tif 000019
【図５】
0005459963.tif 000020
【図１１】
0005459963.tif 000021
【図１２】
0005459963.tif 000022
【図１３】
0005459963.tif 000023
【図１４】
0005459963.tif 000024