(57)【要約】【課題】連続波を利用して目標位置から生体内情報を抽出する技術において、生体内における超音波の周波数依存性減衰を補償する。【解決手段】変調処理部２０は、変調信号に基づいて変調処理された連続波の送信信号を出力する。受信ミキサ３０は、生体内の目標位置との間の相関関係が調整された参照信号を用いて、受信信号に対して復調処理を施すことにより、その目標位置に対応した復調信号を得る。この送受信処理において、受信信号補償部１８は、受信信号に対して、目標位置の深さに応じて、周波数依存性減衰を補償するための処理を実行する。また、参照信号補償部２８は、参照信号に対して、目標位置の深さに応じて、周波数依存性減衰を補償するための処理を実行する。これにより、生体内における周波数依存性減衰が補償される。【選択図】図１
2011036599.tif 000002
【特許請求の範囲】
【請求項１】
周期的な信号に基づいて変調された連続波の送信信号を出力する送信信号出力部と、
前記送信信号に対応した超音波の送信波を生体に送波してその送信波に伴う受信波を生体から受波することにより受信信号を得る送受波部と、
生体内の目標位置との間の相関関係が調整された参照信号を用いて、受信信号に対して復調処理を施すことにより、当該目標位置に対応した復調信号を得る受信信号処理部と、
前記受信信号と前記参照信号のうちの少なくとも一方の信号に対して、前記目標位置に応じて、生体内における超音波の周波数依存性減衰の補償処理を施す信号補償部と、
前記目標位置に対応した復調信号から生体内情報を抽出する生体内情報抽出部と、
を有する、
ことを特徴とする超音波診断装置。
【請求項２】
請求項１に記載の超音波診断装置において、
前記信号補償部は、前記受信信号に対して、前記目標位置の深さに応じた周波数特性のフィルタ処理を施す、
ことを特徴とする超音波診断装置。
【請求項３】
請求項１または２に記載の超音波診断装置において、
前記信号補償部は、前記参照信号に対して、前記目標位置の深さに応じて、前記周波数依存性減衰に相当するフィルタ処理を施すことにより、前記目標位置から得られる受信信号と前記参照信号との間の相関を強める、
ことを特徴とする超音波診断装置。
【請求項４】
請求項１に記載の超音波診断装置において、
前記信号補償部は、前記受信信号に対して、前記目標位置の深さに応じて、前記周波数依存性減衰の逆特性に相当する補償処理を施すことにより、前記目標位置から得られる受信信号と前記参照信号との間の相関を強める、
ことを特徴とする超音波診断装置。
【請求項５】
請求項１から４のいずれか１項に記載の超音波診断装置において、
前記送信信号出力部は、周期的な変調信号を用いて搬送波信号に対して周波数変調処理を施すことにより形成された連続波の送信信号を出力する、
ことを特徴とする超音波診断装置。
【請求項６】
請求項１から４のいずれか１項に記載の超音波診断装置において、
前記送信信号出力部は、周期的な変調信号を用いて搬送波信号に対して位相変調処理を施すことにより形成された連続波の送信信号を出力する、
ことを特徴とする超音波診断装置。
【請求項７】
請求項１から４のいずれか１項に記載の超音波診断装置において、
前記送信信号出力部は、周期的な信号列に基づいた位相シフトキーイングにより位相を変化させて形成された連続波の送信信号を出力する、
ことを特徴とする超音波診断装置。
【請求項８】
請求項１から４のいずれか１項に記載の超音波診断装置において、
前記送信信号出力部は、周期的な信号列に基づいた周波数シフトキーイングにより周波数を変化させて形成された連続波の送信信号を出力する、
ことを特徴とする超音波診断装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、超音波診断装置に関し、特に、連続波を利用する超音波診断装置に関する。
【背景技術】
【０００２】
超音波診断装置の連続波を利用した技術として、連続波ドプラが知られている。連続波ドプラでは、例えば、数ＭＨｚの正弦波である送信波が生体内へ連続的に放射され、生体内からの反射波が連続的に受波される。反射波には、生体内における運動体（例えば血流など）によるドプラシフト情報が含まれる。そこで、そのドプラシフト情報を抽出して周波数解析することにより、運動体の速度情報を反映させたドプラ波形などを形成することができる。
【０００３】
連続波を利用した連続波ドプラは、パルス波を利用したパルスドプラに比べて一般に高速の速度計測の面で優れている。こうした事情などから、本願の出願人は、連続波ドプラに関する研究を重ねてきた。その成果の一つとして、特許文献１において、周波数変調処理を施した連続波ドプラ（ＦＭＣＷドプラ）に関する技術を提案している。
【０００４】
一方、連続波ドプラでは、連続波を利用していることにより位置計測が困難である。例えば、従来の一般的な連続波ドプラの装置（ＦＭＣＷドプラを利用しない装置）では、位置計測を行うことができなかった。これに対し、本願の出願人は、特許文献２において、ＦＭＣＷドプラにより選択的に生体内組織の所望の位置からドプラ情報を抽出することができる極めて画期的な技術を提案している。
【先行技術文献】
【特許文献】
【０００５】
【特許文献１】特開２００５−２５３９４９号公報
【特許文献２】特開２００８−２８９８５１号公報
【発明の概要】
【発明が解決しようとする課題】
【０００６】
特許文献１や特許文献２に記載されたＦＭＣＷドプラの技術は、それまでにない超音波診断の可能性を秘めた画期的な技術である。本願発明者は、この画期的な技術の改良についてさらに研究を重ねてきた。特に、生体内における超音波の周波数依存性減衰に注目した。
【０００７】
本発明は、このような背景において成されたものであり、その目的は、連続波を利用して目標位置から生体内情報を抽出する技術において、生体内における超音波の周波数依存性減衰を補償することにある。
【課題を解決するための手段】
【０００８】
上記目的にかなう好適な超音波診断装置は、周期的な信号に基づいて変調された連続波の送信信号を出力する送信信号出力部と、前記送信信号に対応した超音波の送信波を生体に送波してその送信波に伴う受信波を生体から受波することにより受信信号を得る送受波部と、生体内の目標位置との間の相関関係が調整された参照信号を用いて、受信信号に対して復調処理を施すことにより、当該目標位置に対応した復調信号を得る受信信号処理部と、前記受信信号と前記参照信号のうちの少なくとも一方の信号に対して、前記目標位置に応じて、生体内における超音波の周波数依存性減衰の補償処理を施す信号補償部と、前記目標位置に対応した復調信号から生体内情報を抽出する生体内情報抽出部と、を有することを特徴とする。
【０００９】
望ましい具体例として、前記信号補償部は、前記受信信号に対して、前記目標位置の深さに応じた周波数特性のフィルタ処理を施す、ことを特徴とする。
【００１０】
望ましい具体例として、前記信号補償部は、前記参照信号に対して、前記目標位置の深さに応じて、前記周波数依存性減衰に相当するフィルタ処理を施すことにより、前記目標位置から得られる受信信号と前記参照信号との間の相関を強める、ことを特徴とする。
【００１１】
望ましい具体例として、前記信号補償部は、前記受信信号に対して、前記目標位置の深さに応じて、前記周波数依存性減衰の逆特性に相当する補償処理を施すことにより、前記目標位置から得られる受信信号と前記参照信号との間の相関を強めることを特徴とする。
【００１２】
望ましい具体例として、前記送信信号出力部は、周期的な変調信号を用いて搬送波信号に対して周波数変調処理を施すことにより形成された連続波の送信信号を出力する、ことを特徴とする。
【００１３】
望ましい具体例として、前記送信信号出力部は、周期的な変調信号を用いて搬送波信号に対して位相変調処理を施すことにより形成された連続波の送信信号を出力する、ことを特徴とする。
【００１４】
望ましい具体例として、前記送信信号出力部は、周期的な信号列に基づいた位相シフトキーイングにより位相を変化させて形成された連続波の送信信号を出力する、ことを特徴とする。
【００１５】
望ましい具体例として、前記送信信号出力部は、周期的な信号列に基づいた周波数シフトキーイングにより周波数を変化させて形成された連続波の送信信号を出力する、ことを特徴とする。
【発明の効果】
【００１６】
本発明により、連続波を利用して目標位置から生体内情報を抽出する技術において、生体内における超音波の周波数依存性減衰を補償することが可能になる。
【図面の簡単な説明】
【００１７】
【図１】本発明の実施において好適な超音波診断装置の全体構成を示す図である。
【図２】ＰＳＫにより形成される連続波の送信信号を説明するための図である。
【図３】送信信号と受信信号と復調信号の周波数スペクトラムを示す図である。
【図４】ＰＳＫの場合における位置選択性を説明するための図である。
【図５】ＦＳＫにより形成される連続波の送信信号を説明するための図である。
【図６】周波数依存性減衰の影響を受けた受信信号を示す図である。
【図７】受信信号補償部におけるフィルタ処理を説明するための図である。
【図８】受信信号補償部における補償処理を説明するための図である。
【発明を実施するための形態】
【００１８】
図１は、本発明の実施において好適な超音波診断装置の全体構成を示す図である。送信用振動子１０は生体内へ送信波を連続的に送波し、また、受信用振動子１２は生体内からの反射波を連続的に受波する。このように、送信および受信がそれぞれ異なる振動子で行われて、いわゆる連続波ドプラ法による送受信が実行される。
【００１９】
図１の超音波診断装置で利用される連続波は、周期的な信号に基づいた変調処理により形成される。変調処理された連続波は、変調処理部２０において形成される。変調処理部２０は、変調信号発生部２４から供給される周期的な信号に基づいて、ＲＦ波発振器２２から供給されるＲＦ波に対して変調処理を施して連続波を発生する。変調処理部２０における変調処理としては、周波数変調（ＦＭ）および位相変調（ＰＭ）などのアナログ変調処理と、位相シフトキーイング（ＰＳＫ）や周波数シフトキーイング（ＦＳＫ）などのデジタル変調処理が好適である。変調処理部２０において形成される連続波の波形等については後に詳述する。変調処理部２０は、変調処理された連続波を送信ビームフォーマ（送信ＢＦ）１４へ出力する。
【００２０】
送信ビームフォーマ（送信ＢＦ）１４は、送信用振動子１０が備える複数の振動素子に対して送信信号を出力する。送信ビームフォーマ１４には、変調処理された連続波が入力される。送信ビームフォーマ１４は、その連続波に対して、各振動素子に応じた遅延処理を施して各振動素子に対応した送信信号を形成する。なお、送信ビームフォーマ１４において形成された各振動素子に対応した送信信号に対して、必要に応じて電力増幅処理が施されてもよい。こうして変調処理された連続波による送信ビームが形成される。
【００２１】
受信ビームフォーマ（受信ＢＦ）１６は、受信用振動子１２が備える複数の振動素子に対応した複数の受波信号を整相加算処理して受信ビームを形成する。受信ビームフォーマ１６は、各振動素子から得られる受波信号に対してその振動素子に応じた遅延処理を施し複数の振動素子から得られる複数の受波信号を加算処理することにより受信ビームを形成する。なお、複数の振動素子から得られる複数の受波信号の各々に対して低雑音増幅等の処理を施してから受信ビームを形成してもよい。こうして受信ビームに沿った受信ＲＦ信号が得られる。
【００２２】
受信信号補償部１８は、受信ビームフォーマ１６から得られる受信ＲＦ信号に対して、生体内における超音波の周波数依存性減衰に関する補償処理を施す。この補償処理については後に詳述する。
【００２３】
受信ミキサ３０は、補償処理後の受信ＲＦ信号に対して直交検波を施して複素ベースバンド信号を生成する回路であり、２つのミキサ３２，３４で構成される。各ミキサは受信ＲＦ信号を所定の参照信号と混合する回路である。
【００２４】
受信ミキサ３０の各ミキサに供給される参照信号は、変調処理部２０から出力される連続波に基づいて生成される。つまり、変調処理部２０から出力される連続波が遅延回路２５において遅延処理され、さらに、参照信号補償部２８において補償処理されてから、ミキサ３２には補償処理後の連続波が直接供給され、一方、ミキサ３４には補償処理後の連続波がπ／２シフト回路２６を経由して供給される。
【００２５】
なお、遅延回路２５は、目標位置の深さに応じた遅延量だけ連続波に遅延処理を施し、遅延された参照信号を出力する。これにより、後に詳述するように、目標位置に対応した復調信号が得られる。また、参照信号補償部２８は、遅延回路２５において遅延処理された連続波である参照信号に対して、生体内における超音波の周波数依存性減衰に関する補償処理を施す。この補償処理についても後に詳述する。
【００２６】
π／２シフト回路２６は、遅延処理と補償処理を施された連続波の位相をπ／２だけずらす回路である。この結果、２つのミキサ３２，３４の一方から同相信号成分（Ｉ信号成分）が出力され、他方から直交信号成分（Ｑ信号成分）が出力される。なお、受信ミキサ３０の後段に設けられたＬＰＦ（ローパスフィルタ）３６，３８により、同相信号成分および直交信号成分の各々の高周波数成分がカットされ、検波後の必要な帯域のみの復調信号が抽出される。ＬＰＦ３６，３８から出力される復調信号は、ＦＦＴ処理部５０へ出力される。
【００２７】
ＦＦＴ処理部（高速フーリエ変換処理部）５０は、復調信号（同相信号成分および直交信号成分）の各々に対してＦＦＴ演算を実行する。その結果、ＦＦＴ処理部５０において復調信号が周波数スペクトラムに変換される。
【００２８】
ドプラ情報解析部５２は、周波数スペクトラムに変換された復調信号からドプラ信号を抽出する。後に詳述するが、図１の超音波診断装置では、遅延回路２５における遅延処理により目標位置が設定され、ドプラ情報解析部５２において目標位置からのドプラ信号が選択的に抽出される。ドプラ情報解析部５２は、例えば、時間的に変化するドプラ信号の表示波形を形成する。なお、生体内の各深さ（各位置）ごとにドプラ信号を抽出して、例えば、超音波ビーム（音線）上の各深さごとに生体内組織の速度を算出し、リアルタイムで出力してもよい。また、超音波ビームを走査させて二次元的あるいは三次元的に生体内組織の各位置の速度を算出してもよい。
【００２９】
表示部５４は、ドプラ情報解析部５２において形成されたドプラ信号の波形などを表示する。なお、図１に示す超音波診断装置内の各部は、システム制御部６０によって制御される。つまり、システム制御部６０は、送信制御や受信制御や表示制御などを行う。
【００３０】
次に、図１の超音波診断装置における目標位置の選択性について説明する。図１の超音波診断装置では、変調処理された連続波に対応した超音波を送受波して受信信号を得て、生体内の目標位置の深さに応じて参照信号と受信信号との間の遅延関係を調整し、目標位置からの受信信号と参照信号との間の相関を強めて復調処理を施すことにより、目標位置からの生体内情報としてドプラ情報を選択的に抽出する。変調処理部２０における変調方式としては、周波数変調（ＦＭ）や位相変調（ＰＭ）などのアナログ変調処理と、位相シフトキーイング（ＰＳＫ）や周波数シフトキーイング（ＦＳＫ）などのデジタル変調処理が好適である。そこで、各変調方式ごとに位置選択性について説明する。なお、既に図１に示した構成（部分）については、以下の説明において図１の符号を利用する。
【００３１】
＜周波数変調（ＦＭ）および位相変調（ＰＭ）による位置選択性＞
図１の変調処理部２０においてアナログ変調処理を施す場合、例えば周波数変調（ＦＭ変調）や位相変調（ＰＭ変調）を施す場合、変調処理部２０は、ＲＦ波発振器２２から供給されるＲＦ波、および、変調信号発生部２４から供給される正弦波や鋸歯状波などの変調波に基づいて連続波を発生する。
【００３２】
ここでは、変調波として正弦波を用いてＦＭ変調処理する場合を例にとり、位置選択性について説明する。周波数ｆ０のＲＦ波（搬送波）に対して、周波数ｆｍの正弦波によりＦＭ変調を施したＦＭＣＷ送信波は次式のように表現できる。次式において、Δｆは周波数変動幅の０−Ｐ値（ゼロピーク値：最大周波数偏移）であり、最大周波数偏移Δｆと変調周波数ｆｍの比であるβはＦＭの変調指数（変調度）である。
【００３３】
【数１】
2011036599.tif 000003
【００３４】
また、ドプラシフトを伴う場合のＦＭＣＷ受信波は、生体における往復の減衰をαとすると次式で表現できる。なお、次式において、ｆｍに対するドプラシフトは、ｆ０のシフト分ｆｄに比べて小さいので無視している。
【００３５】
【数２】
2011036599.tif 000004
【００３６】
数２式で表される受信波形は、超音波振動子を介して受信される信号波形（受信ＲＦ信号）である。ＦＭＣＷドプラでは、受信ＲＦ信号に対する復調処理において、ＦＭＣＷ送信波を参照信号として受信波と乗算を行う。図１を利用して説明したように、変調処理部２０から出力されるＦＭ連続波が参照信号として利用され、遅延回路２５において遅延処理され、参照信号補償部２８を経て、ミキサ３２には遅延処理されたＦＭ連続波が直接供給され、一方、ミキサ３４には遅延処理されたＦＭ連続波がπ／２シフト回路２６を経由して供給される。したがって、ミキサ３２へ供給される参照信号ｖｒＩ（ｔ）と、ミキサ３４へ供給される参照信号ｖｒＱ（ｔ）は、次式のように表現できる。
【００３７】
【数３】
2011036599.tif 000005
【００３８】
数３式において、φｍｒは、遅延回路２５における遅延処理により任意に設定できる参照信号の位相を示しており、φ０ｒは、任意に設定した参照信号の位相に対応して決まる搬送波の位相変化量を示している。
【００３９】
受信ミキサ３０では、復調処理として直交検波が行われる。つまり、ミキサ３２において、受信ＲＦ信号ｖＲ（ｔ）と参照信号ｖｒＩ（ｔ）の乗算に相当する処理が実行され、また、ミキサ３４において、受信ＲＦ信号ｖＲ（ｔ）と参照信号ｖｒＱ（ｔ）の乗算に相当する処理が実行される。
【００４０】
ミキサ３２における受信ＲＦ信号ｖＲ（ｔ）と参照信号ｖｒＩ（ｔ）の乗算ｖＤＩ（ｔ）は次式のように表現される。なお、次式の計算途中において、周波数２ｆ０の成分が消去されている。これは、ＬＰＦ３６によって除去される。
【００４１】
【数４】
2011036599.tif 000006
【００４２】
ここで、ベッセル関数に関する次の公式を利用する。
【００４３】
【数５】
2011036599.tif 000007
【００４４】
数５式の公式を用いると、数４式はさらに次式のように計算される。
【００４５】
【数６】
2011036599.tif 000008
【００４６】
一方、ミキサ３４における受信ＲＦ信号ｖＲ（ｔ）と参照信号ｖｒＱ（ｔ）の乗算ｖＤＱ（ｔ）は次式のように表現される。なお、次式の計算途中において、周波数２ｆ０の成分が消去されている。これは、ＬＰＦ３８によって除去される。
【００４７】
【数７】
2011036599.tif 000009
【００４８】
ここで、数６式のｖＤＩ（ｔ）と数７式のｖＤＱ（ｔ）とに基づいて、複素ベースバンド信号を定義する。まず、ｖＤＩ（ｔ）とｖＤＱ（ｔ）に含まれている直流（ＤＣ）成分、変調周波数ｆｍの偶数次高調波成分を次式のように表現する。
【００４９】
【数８】
2011036599.tif 000010
【００５０】
次に、ｖＤＩ（ｔ）とｖＤＱ（ｔ）に含まれている変調周波数ｆｍの成分、変調周波数ｆｍの奇数次高調波成分を次式のように表現する。
【００５１】
【数９】
2011036599.tif 000011
【００５２】
数８式と数９式から、直交検波後のベースバンド信号において、ドプラシフトｆｄを含んだドプラ信号は、ＤＣ成分と変調周波数ｆｍの成分と変調周波数ｆｍの高調波成分とからなる複数の成分の各々についての両側帯波として出現することがわかる。
【００５３】
ここで、受信信号と参照信号の位相を互いに揃えた場合、つまり、遅延回路２５における遅延処理によりφｍｒを調整してφｍと一致させた場合（φｍｒ＝φｍ）を考える。φｍｒとφｍを一致させた場合には、数４式におけるｋが０となる。この結果をベッセル関数に適用すると、第１次ベッセル関数の性質により、次式のように、０次のベッセル関数の値のみが１となり、それ以外のベッセル関数の値は０となる。
【００５４】
【数１０】
2011036599.tif 000012
【００５５】
数１０式に示す結果を数８式と数９式に適用すると次式のとおりとなる。
【００５６】
【数１１】
2011036599.tif 000013
【００５７】
数１１式は、参照波（参照信号）の位相φｍｒを送受信間の位相差φｍに設定すると、圧縮変換により、ＤＣ成分（直流信号成分）に対応したドプラ信号のみが抽出できることを示している。
【００５８】
上述した数２式の受信波形は、ある深さからの受信信号の波形である。これに対し、ＦＭＣＷ送信波を利用して、実際に受信用振動子１２において得られる受信信号は、複数の深さからの信号が合成された受信信号である。受信ミキサ３０においては、複数の深さからの信号が合成された受信信号と参照信号との乗算に相当する処理が実行される。
【００５９】
数８式などに現れた直流信号成分に対応したドプラ信号の振幅を支配するＪ０（ｋβ）は、第１次ベッセル関数の性質により、ｋβが０のときに最大値である１となり、ｋβが０からずれると急激に小さくなる。そのため、遅延回路２５においてφｍｒを調整し、目標位置から得られる受信信号のφｍと一致させると、目標位置におけるＪ０（ｋβ）が最大値である１となり、目標位置以外におけるＪ０（ｋβ）が極端に小さな値となる。したがって、遅延回路２５においてφｍｒを調整し、目標位置から得られる受信信号のφｍと一致させることにより、目標位置におけるドプラ信号（直流信号成分）を選択的に抽出することができる。
【００６０】
このように、周波数変調処理を利用した場合において、目標位置からの受信信号に対応するように参照信号の位相を調整することにより、その目標の深さからのドプラ情報を選択的に検出する位置選択性が実現される。
【００６１】
なお、周波数変調処理に換えて、周波数変調処理と同じ角度変調の方式として一般的によく知られている位相変調処理（ＰＭ処理）を利用してもよい。つまり、搬送波信号に対して位相変調処理を施すことにより、変調処理部２０から出力されるＦＭ連続波と同じ波形あるいは同等な波形を形成してもよい。
【００６２】
＜位相シフトキーイング（ＰＳＫ）による位置選択性＞
図２は、位相シフトキーイング（ＰＳＫ）により形成される連続波の送信信号を説明するための図である。図２（Ａ）には、ＲＦ波発振器２２（図１）から出力されるＲＦ信号（ＲＦ波）の波形が示されている。図２（Ｂ）には、変調信号発生部２４（図１）から出力される周期的な信号列の一例が示されている。変調信号発生部２４は、例えば図２（Ｂ）に示すような、ランダムに値を変化させた２値符号（擬似ランダム信号）を発生する。
【００６３】
図２（Ｃ）には、ＰＳＫ変調器として機能する変調処理部２０（図１）において形成される変調された連続波（送信信号）が示されている。変調処理部２０は、図２（Ａ）のＲＦ信号に対して、図２（Ｂ）の２値符号に基づいて、位相シフトキーイング（ＰＳＫ）の変調処理を施す。こうして、例えば図２（Ｃ）の送信信号に対応した連続波の超音波が送信用振動子１０（図１）から出力され、受信用振動子１２（図１）を介して生体内から受信信号が得られる。
【００６４】
図３は、送信信号と受信信号と復調信号の周波数スペクトラムを示す図である。図３（Ａ）には、ＰＳＫ変調器として機能する変調処理部２０（図１）において形成される送信信号、つまりＰＳＫ変調された連続波の周波数スペクトラムが示されている。周波数ｆ０はＲＦ信号の周波数である。ＲＦ信号の周波数ｆ０を中心として広がっている側帯波の周波数間隔は、擬似ランダム信号（図２（Ｂ）の２値符号）の繰り返し周波数ｆｐである。また、周波数ｆ０を中心として広がっている側帯波の電力が０（ゼロ）となる、いわゆるヌル（ｎｕｌｌ）点が存在する。周波数ｆ０からヌル点までの周波数間隔は、擬似ランダム信号（図２（Ｂ）の２値符号）の１ビットの時間間隔Ｔの逆数となる。
【００６５】
図３（Ｂ）には、受信信号の周波数スペクトラムが示されている。受信信号は、生体内における減衰を無視すると、送信信号に対応した波形となる。したがって、図３（Ｂ）に示す受信信号の周波数スペクトラムは、図３（Ａ）に示す送信信号の周波数スペクトラムとほぼ同じである。但し、生体内における超音波の伝搬時間に応じて、送信信号と受信信号との間では位相が異なる。
【００６６】
ＰＳＫの場合においても、変調処理部２０で形成された送信信号に対して遅延処理を施して参照信号を形成し、受信ミキサ３０（図１）においてその参照信号を用いて受信信号に対してミキサ処理（参照信号と受信信号の乗算）が行われる。後に詳述するが、このミキサ処理において、遅延処理された参照信号の位相に対応する深さ（目標位置の深さ）からの受信信号と参照信号との間の相関が強められて最大となり、その他の深さからの受信信号と参照信号との間の相関が極端に小さくなる。
【００６７】
図３（Ｃ）には、ミキサ処理により得られる復調信号の周波数スペクトラムが示されている。図３（Ｃ）の復調信号は、相関が最大の場合における参照信号と受信信号の乗算結果に相当する。つまり、目標位置からの受信信号と、目標位置の深さに位相を合わせた参照信号との間の乗算結果が、図３（Ｃ）の復調信号となる。
【００６８】
図３（Ｃ）に示す復調信号には、直流信号成分と、ＲＦ信号の周波数ｆ０の２倍の高調波成分が含まれている。ドプラ信号は、これらの成分に付着した形で出現する。なお、ＬＰＦ３６，３８（図１）において、高調波成分がカットされて直流信号成分のみが抽出されるため、ＦＦＴ処理部５０（図１）においては、図３（Ｃ）に示す直流信号成分の周波数スペクトラムのみが形成される。そして、ドプラ情報解析部５２（図１）において、図３（Ｃ）に示す直流信号成分の周波数スペクトラムからドプラ信号が抽出され、ドプラシフト量などに基づいて、目標位置に存在する血流の流速などが算出される。受信ミキサ３０において、直交検波を施しているため、流速の極性を判断することもできる。直流信号成分の周波数スペクトラムからクラッタ信号を抽出して、目標位置に存在する血管壁の位置などを算出してもよい。
【００６９】
図４は、ＰＳＫの場合における位置選択性を説明するための図である。受信信号と参照信号の相関の鋭さは、変調信号発生部２４（図１）において形成される周期的な信号列のシーケンスに依存する。相関性を鋭くするためには、周期的な信号列である擬似ランダム信号の符号系列として、ＰＮ系列などのパルス圧縮などで実用化されている符号系列を用いればよい。簡単な例として、ｎ＝３のＰＮ符号を用いた場合の位置選択性について、図４を用いて説明する。
【００７０】
ｎ＝３の場合のＰＮ符号の長さは７（＝２３―１）ビットである。このシーケンスが際限なく繰り返すので、この擬似ランダムパターンは繰り返し周期の逆数の線スペクトラムを持っている。この信号を用いて周波数ｆ０の搬送波に０−πの２相のＰＳＫ変調をかけると、その時間波形は、先に説明した図２（Ｃ）のようになる。
【００７１】
受信信号は、送信信号が目標の深さに応じた遅延時間だけ遅れ、また組織によって減衰した信号である。その減衰を無視すると、例えば図４の受信信号の波形が得られる。送信信号を遅延処理して得られる参照信号の位相をφ１〜φ６まで変化させて受信信号と乗算した結果（乗算器出力）が図４に示されている。
【００７２】
図４から、参照信号と受信信号の位相が一致したφ３の場合に、乗算器出力（ミキサ出力）の直流成分が最大になる。また、交流成分としては、搬送波およびその高調波成分のみとなるのも、参照信号と受信信号の位相が一致した場合の特徴である。この信号の周波数スペクトラムは、図３（Ｃ）に示したとおりである。また、図４から、位相がφ３以外の場合には、乗算器出力として正と負の電圧がランダムに発生するので、これらの平均電圧は非常に小さくなる。
【００７３】
このように、位相シフトキーイングを利用した場合においても、目標位置からの受信信号に対応するように参照信号の位相を調整することにより、その目標の深さからのドプラ情報を選択的に検出する位置選択性が実現される。
【００７４】
＜周波数シフトキーイング（ＦＳＫ）による位置選択性＞
図５は、周波数シフトキーイング（ＦＳＫ）により形成される連続波の送信信号を説明するための図である。図５（Ｉ）には、変調信号発生部２４（図１）から出力される周期的な信号列の一例が示されている。変調信号発生部２４は、例えば図５（Ｉ）に示すような、ランダムに値を変化させた２値符号（擬似ランダム信号）を発生する。周期的な信号列である擬似ランダム信号の符号系列としてはＰＮ系列などの符号系列を用いればよい。
【００７５】
図５（ＩＩ）には、ＦＳＫ変調器として機能する変調処理部２０（図１）において形成される変調された連続波（送信信号）が示されている。変調処理部２０は、図５（Ｉ）の２値符号に基づいた周波数シフトキーイングにより、ＲＦ波（搬送波）に対してデジタル変調処理を施して連続波の送信信号を形成する。変調処理部２０は、例えば２値符号が「１」の期間において周波数ｆ１とし、２値符号が「０」の期間において周波数ｆ２とすることにより、図５（ＩＩ）の送信信号を形成する。
【００７６】
こうして、例えば図５（ＩＩ）の送信信号に対応した連続波の超音波が送信用振動子１０（図１）から出力され、受信用振動子１２（図１）を介して生体内から受信信号が得られる。そして、周波数シフトキーイング（ＦＳＫ）の場合においても変調処理部２０で形成された送信信号に対して遅延処理を施して参照信号を形成し、受信ミキサ３０（図１）においてその参照信号を用いて受信信号に対してミキサ処理（参照信号と受信信号の乗算）が行われる。このミキサ処理において、ＰＳＫの場合と同様に、遅延処理された参照信号の位相に対応する深さ（目標位置の深さ）からの受信信号と参照信号との間の相関が強められて最大となり、その他の深さからの受信信号と参照信号との間の相関が極端に小さくなる。
【００７７】
このように、周波数シフトキーイングを利用した場合においても、目標位置からの受信信号に対応するように参照信号の位相を調整することにより、その目標の深さからのドプラ情報を選択的に検出する位置選択性が実現される。
【００７８】
さらに、図１の超音波診断装置は、受信信号と参照信号のうちの少なくとも一方の信号に対して、目標位置に応じて、生体内における超音波の周波数依存性減衰の補償処理を施している。そこで、その補償処理について説明する。
【００７９】
＜周波数依存性減衰の補償処理＞
超音波が生体内を伝搬する際の減衰は、周波数と伝搬距離にほぼ比例することが知られており、周波数依存性減衰（ＦＤＡ：Frequency dependent attenuation）と呼ばれる。この特性は、送受信される超音波の周波数スペクトラムの形に影響を与える。
【００８０】
図６は、周波数依存性減衰の影響を受けた受信信号を示す図である。図６には、生体内の浅い位置（浅部）から深い位置（深部）までの複数の深さの各々に対応した受信信号の電力スペクトラムが図示されている。
【００８１】
図６に示すように、深くなるにつれて受信信号が全体的に減衰しており、それに加え、深くなるにつれて低域の周波数成分に比べて高域の周波数成分が大きく減衰している。そして、各受信信号の周波数軸方向に関する対称性も崩れている。
【００８２】
こうした、周波数依存性減衰を補償するために、図１の受信信号補償部１８と参照信号補償部２８において以下に説明する処理が行われる。
【００８３】
図７は、受信信号補償部１８におけるフィルタ処理を説明するための図である。図６を利用して説明したように、周波数依存性減衰の影響により、受信信号の周波数特性が深さに応じて異なるため、受信信号補償部１８は、受信信号に対して、図７に示すように深さに応じて異なる周波数特性のフィルタ処理を施す。図７には、生体内の浅い位置（浅部）から深い位置（深部）までの複数の深さの各々に対応したバンドパスフィルタ（ＢＰＦ）の周波数特性が図示されている。
【００８４】
受信信号補償部１８は、目標位置の深さに応じたバンドパスフィルタを選択して、受信ビームフォーマ１６から得られる受信信号に対してフィルタ処理を施す。その目標位置は遅延回路２５に設定される目標位置である。例えば、システム制御部６０により、受信信号補償部１８と遅延回路２５に対して同一の目標位置（同じ深さ）が指示される。
【００８５】
受信信号補償部１８において目標位置の深さに応じたバンドパスフィルタが選択されるため、浅部においては高域を含む比較的広い帯域に亘って受信信号が利用され、一方、深部においては、低域側の比較的狭い帯域のみが利用される。これにより、周波数依存性減衰の影響を受けて減衰した受信信号の有効な帯域のみを利用することができ、信号対雑音比（ＳＮＲ）を向上させることができる。なお、目標位置の深さに応じて、受信信号の周波数スペクトラムの周波数軸方向に関する対称性を補償するようにバンドパスフィルタを構成してもよい。
【００８６】
各変調方式ごとの位置選択性の説明でも述べたように、図１の超音波診断装置では、目標位置から得られる受信信号と参照信号との間の相関関係に基づいて、目標位置に対応した復調信号を抽出している。そのため、目標位置から得られる受信信号と参照信号との間の相関関係が崩れてしまうと、目標位置に対応した復調信号を適切に抽出できなくなる。例えば、周波数依存性減衰の影響により、図６に示したように受信信号が減衰してしまうと、減衰のない参照信号との間の相関関係が崩れてしまい、例えば、復調信号内に、直流信号成分以外の高次不要波が発生してしまう。
【００８７】
そこで、参照信号補償部２８は、参照信号に対して、目標位置の深さに応じて、周波数依存性減衰に相当するフィルタ処理を施す。これにより、目標位置から得られる受信信号と参照信号が共に周波数依存性減衰を受けた信号となり、これら二つの信号間の相関が強められる。なお、この場合においては、例えば、システム制御部６０により、受信信号補償部１８と遅延回路２５と参照信号補償部２８に対して同一の目標位置（同じ深さ）が指示される。
【００８８】
また、参照信号は減衰させずにそのままとし、受信信号の減衰を補償して、参照信号と受信信号の相関を強めるようにしてもよい。この場合、受信信号補償部１８は、受信信号に対して、目標位置の深さに応じて、周波数依存性減衰の逆特性に相当する補償処理を行い、目標位置から得られる受信信号と参照信号との間の相関を強める。
【００８９】
図８は、受信信号補償部１８における補償処理を説明するための図である。図６に示したように、周波数依存性減衰の影響により、深くなるにつれて受信信号が全体的に減衰して、さらに、深くなるにつれて低域の周波数成分に比べて高域の周波数成分が大きく減衰する。
【００９０】
そこで、受信信号補償部１８は、図６に示した周波数依存性減衰の逆特性に相当する補償処理を行う。つまり、図８に示すように、深さに応じて異なる重み特性の補償処理を施す。図８には、生体内の深い位置（深部）から浅い位置（浅部）までの複数の深さの各々に対応した重み特性が図示されており、深部ほど受信信号が全体的に大きくなるように補償され、それに加え、深部ほど低域の周波数成分に比べて高域の周波数成分が大きく補償される。また、各受信信号の周波数軸方向に関する対称性を補償するようにしてもよい。
【００９１】
受信信号補償部１８は、目標位置の深さに応じた重み特性を選択して、受信ビームフォーマ１６から得られる受信信号に対して補償処理を施す。その目標位置は遅延回路２５に設定される目標位置である。例えば、システム制御部６０により、受信信号補償部１８と遅延回路２５に対して同一の目標位置（同じ深さ）が指示される。こうして、目標位置から得られる受信信号が周波数依存性減衰を受ける前のスペクトラム形状に戻され、参照信号との間の相関が強められる。
【００９２】
なお、図１においては、１つの受信ミキサ３０を備えた構成を示したが、複数の受信ミキサ３０を並列的に設けて、複数の目標位置に対応した受信信号を並列的に処理する並列構成としてもよい。この並列構成の場合には、複数の受信ミキサ３０の各々に対して受信信号が供給され、また、複数の受信ミキサ３０の各々に対して参照信号が供給される。そのため、並列構成の場合には、複数の受信ミキサ３０の各々に対応した受信信号補償部１８が設けられ、また、複数の受信ミキサ３０の各々に対応した参照信号補償部２８が設けられる。そして、各受信ミキサ３０ごとに、受信信号と参照信号が適宜補償され、複数の受信ミキサ３０から複数の目標位置に対応した復調信号が並列的に出力される。例えば、１本の超音波ビームに沿った複数の復調信号が並列的に出力される。
【００９３】
以上、本発明の好適な実施形態を説明したが、上述した本発明の好適な実施形態等は、あらゆる点で単なる例示にすぎず、本発明の範囲を限定するものではない。本発明は、その本質を逸脱しない範囲で各種の変形形態を包含する。
【符号の説明】
【００９４】
１８  受信信号補償部、２０  変調処理部、２２  ＲＦ波発振器、２４  変調信号発生部、２５  遅延回路、２８  参照信号補償部、３０  受信ミキサ、５０  ＦＦＴ処理部、５２  ドプラ情報解析部。
【図１】
2011036599.tif 000014
【図２】
2011036599.tif 000015
【図３】
2011036599.tif 000016
【図４】
2011036599.tif 000017
【図５】
2011036599.tif 000018
【図６】
2011036599.tif 000019
【図７】
2011036599.tif 000020
【図８】
2011036599.tif 000021