JPS59226841A - 超音波による温度測定装置 - Google Patents

Abstract

(57)【要約】本公報は電子出願前の出願データであるた
め要約のデータは記録されません。

Description

【発明の詳細な説明】 （Ａ）発明の技術分野 本発明は、生体組織等の超音波媒体の温度または温度変
化量の空間的分布を測定する装置に係り２、温度Ｃ２お
よび、密度ρの温度依存性を利用した、超音波による温
度測定装置に関する。

（Ｂ）従来技術とその問題点 従来、生体組織の如き超音波媒体内の温度分布を得る方
法としては、音速Ｃの温度依存性を利用して、音速の空
間分布を測定する事により温度の空間分布を得ようとい
う方法が知られているが、この方法には以下の様な問題
がある。

第１に、音速の測定は、いわゆる透過法で行ない、音速
の空間分布は、Ｘ線ＣＴの場合と同様の像の再構成手法
によって得るので、計算量が膨大なものとなり、時間が
かかる。

第２に、像の再構成を行なうためには、被測定体の全周
（３６０度）にわたって透過法による音速の測定ができ
る必要があるが、実際の生体組織では、骨や体内のガス
の影響で超音波が透過できない部分があり、全周にわた
っての測定が不可能な場合が大部分である。従って、欠
落のあるデータを用いた像の再構成を行なう事になり、
得られた音速分布、従って音速分布から得られる温度分
布、は誤差の多いものとなる。

（Ｃ）発明の目的と構成 本発明の目的は、Ｘ線ＣＴの場合に用いられる如き像の
再構成手法を用いる事なく、高速かつ容易に超音波媒体
内の温度分布または温度変化分布を測定できる装置を提
供するにある。

本発明は、超音波媒体の特性量である等価非線Ｘ線ＣＴ
の如き像の再構成手法を用いる事なく容によって異った
値として観測される事を利用して、温度分布または温度
変化分布を推定するようにしたものである。

以下では、先ず原理的な説明を行ない、続いて実施例に
ついて説明する。

超音波媒体内の音圧がゼロの時の音速をＣｏ　。

密度をρ０とすると測定波と直交する方向からＰなる音
圧が加えられたときの音速Ｃは 形パラメータ」と呼ばれるものであり、不均一媒質にお
いては場所により異った値を取る。

ｉとｊは直交する方向であり、Ｓば等エントロピーであ
ることを示す。

一方、一般に音響の分野で用いられている非線形パラメ
ーターは次のように定義され′ζいる。

一との間には という関係があり、液体、生体組織ではν嬌０．５従っ
て、ボンピング波の音圧Ｐにより、音速Ｃは、 だけ変化する事になる。

今、第１図に示す様に、送信用超音波振動子ＸＡから測
定用の連続波超音波ビームＷｍを図のＺ軸方向に被測定
媒体中に送り込み、送信用超音波振動子ＸＢでこのビー
ムを受信する。この時、測定用ビームＷｍと交叉する方
向から、別の超音波振動子Ｘｐにより平面パルス波Ｗｐ
をボンピング波として図示の様にＸ軸方向に送り込む。

以下簡単のため交叉角は直角であるとする。また、Ｗｍ
の音圧はＷｐの音圧に比べて充分に小さく、式（１）〜
（５）のＰとしてはＷｐによるものだけを考慮すれば良
いものとする。

ボンピング波Ｗｐが測定用ビームＷｍと交叉しいる間に
は、測定用ビームＷｍの各Ｚにおける波形は式（団によ
り、場所により異なる音速変化・・　（６） を受ける事になり、従って場所（Ｚ）により異った△Ｃ
（ｚ）に比例する位相変化を受ける事になる。従って、
受信振動子Ｘｇからの受信信号の位相情報及び予め測定
しておいたＰを用いて逆に第２図は、測定用ビームＷｍ
が進行するにつれてポンピング波Ｗｐからどの様な音圧
を受けるかという事を模式的に示したものである。

ポンピング波の音圧をＰ　（ｆ）とする。時刻を一〇で
ｚ＝ｚＯにあって音圧Ｐ　（ｏ）を受けていた測定波（
第２図（Ａ）参照）は、時刻を一△ｔにはｚ＝ｚ□＋△
ｚ＝ｚｏ　＋Ｃｏ△ｔの位置にあり、この時にはｐ（ｔ
＋△ｔ）の音圧を受ける（第２図（Ｂ）参照）。以下同
様にして時刻ｔ＝、２△ｔにおいて第２図（Ｃ）図示、
時刻ｔ＝３△ｔにおいて第２図（Ｄ）の図示の如き音圧
を受ける事になる。これから判る様に、一般に、ｔ＝Ｑ
においてｚ＝ｚＯにあった測定波が、ｚ＝ｚで受ける音
速変化は ・・　（７） 但し、ここで△Ｃ（ｚ）の値はｃｏに比べて小さく、△
ｔは測定波が２＝ＺＱからｚ＝ｚまで移△ｔ＝□　・　
・　・　（８） Ｏ で与えられると近似した。

従って、１＝０でＺ＝２６にあった測定波が２＝ｚで受
ける位相の変化はｋを比例定数として、となり、この測
定波が送信振動子Ｘθを出てから受信振動子Ｘｔ３に受
信されるまでの間に受ける位相変化の総和は となる。ここで と置くと式（１０ンは となる。

ここで、式（１２）の「＊」はいわゆるたたみ込みｇ　
（ｚ）なる関数が重み関数となってφ（ｚｏ）を与える
事を示している。

式（０）フーリエ変換すれば、 （ｚ）を求める際には必ずしも式（１５）の如きフｅ −リエ逆変換を行う必要はなく、周波数特性が−に 出力として得られる。またＧ　（ｗ）は式１１を用いて
Ｐ　（ｔ）から求めておく事ができる。この様にして得
られた−（２ンを映像として表示すれば。

、Ａ −，（ｚ）の空間分布像を得る事もできる。

以上の説明ではＣｏ、　　ρＯは定数として扱ったさて
、超音波媒体の温度が変化すると、ρ０゜（Ｚ）の変化
を知る事によって、各部位での温度変化を知る事ができ
る。

特別な場合として、密度ρ。は温度によって余り変化せ
ず、音速Ｃｏ及び等価非線形パラメータすると媒体の温
度が未知の値Ｔｏから未知の値Ｔ１に変化した場合 ρ、（Ｚ）ＴＩ−ρｏ（Ｚ）ｔｏ　−ρｏ　（ｚ）Ｃ１
＋　　（Ｚ）ＴＩ　　＝ＣＯ（ｚ）Ｔｈ４α△′ｒ　　
但し△Ｔ＝Ｔ＋　　　　Ｔ。

α、βは定数　・・・（１６ン （２ンの差△は １ ρｏ　　（ｚ）　　Ｃｏ　　（Ｚ）ＴＯ＋α△Ｔ・　・
　・　（１７） となる。但し、式（１７）の計算の途中で、て、 β△Ｔ　　　　　　　　α△Ｔ 前述の様にそれぞれ実測する事ができるから、α。

れば、式（１７ンから八Ｔを求める事ができる。実際、
α、β、　　Ｃｏ　　（ｚ）ｔｏ　＋　　ρ。　（ｚ）
は例えば肝臓や脂肪等の種々の超音波媒体について予め
実測しき、且つρｏ　　（ｚ）、Ｃｏ　（ｚ）−ｒｂを
予め知る事ておく事ができる。従って式（１７）から八
Ｔを求める事が可能となる。

尚、式（１７）で与えられる量と温度変化量との対応が
ついておれば、式（１７）−ｒ与えられる量を直接に出
力しても良い事は言うまでもない。

以上、本発明の詳細な説明を行なったが、以下実施例に
ついて第３図を用いて説明する。

第３図において、１はボンピング波の送信のタイミング
を発生ずるタイミング制御部、２は測定用連続波のだめ
の発振器、３は発振器出力を受ｋＪで振動子を駆動する
ためのドライバ、４は測定波送信用の振動子ＸＡ、５は
被測定超音波媒体Ｍ。

６は測定波受信用の振動子Ｘ８，７は受信増幅器。

義されるフィルタ、ＩＯはボンピング波相のドライバ、
１１はボンピング波発生用の振動子Ｘｐである。１２ば
、図示しない温度操作手段によって１２の出力とフィル
タ９の出力との差を得る減算部である。またｘｌはＸｐ
から測定用ビームまでの距離、２工は被測定区間の距離
、Ｚ＋　は被測定区間とＸＩＩＩとの距離である。第３
図において、発振器２ばタイミング制御部１と無関係に
動作するように措かれているが、同期動作しても良いこ
とは言うまでもない。フィルタ９は、例えばＣＣＤを用
いたトランスバーサルフィルタで実現できる。

記憶部１２は、例えば第４図に示す如く、Ａ／Ｄ変換器
１５．ＩＣメモリ１６．Ｄ／Ａ変換器１７゜制御部１８
により実現でき、ポンプ波の送信に同期して動作するよ
うにタイミング制御部１により制御される。第５図には
第３図の主要部の時間波形を示しており、第５図には第
３図と同じ信号名称を記しである。ボンピング波を送信
してからボンピング波が測定波の位置に到達するまでの
時間た測定波が初めて振動子Ｘ８に到達するのに要する
時間（ル）との和をｔｌとすると、受信増幅器７０出力
ＶＱは、第５図に示す様にボンピング波の送信から時間
Ｌ１の後に位相変調を受は始め、ボンピング波による位
相変調を受けた測定波が被の間位相変調出力を出し続り
る。位相検出器８ば、発振器２の出力と■ρとの位相を
比較し、Ｖ９の位相を時間の関数φ（１）として、従っ
て座標２の関数として出力とする。この出力が、フィル
タして、従って座標２の関数として得られる。図示しな
い温度操作手段によって超音波媒体５の温度が変化する
前のフィルタ出力がタイミング制御部１の制御の下に記
憶部１２に記↑、ａされる。温度操作手段によって媒体
の温度が変化した後、記憶部１２の出力とフィルタ９の
出力との差、つまり式（１７）で示される量が減算１３
で最終出力として得られる。

間分布を、Ｘ線ＣＴの場合の如く複雑で時間のがかる像
の再構成手法を用いる事なく、簡単な構成で高速に得る
事ができる。

特許請求の範囲第２項は、第３図におけるフィルタ９の
挿入位置が、減算部１３の後でも良い事を述べたもので
ある。第３図におけるフィルタ９はボンピング波の周波
数特性の逆特性にほぼ等しい特性を持つ線形フィルタで
あるがら、２つの入力Ｘｏ　　（ｔ）とｘｌ　（ｔ）と
の差Ｘ。（ｔ）　−ｘｌ（１）をフィルタに通した出力
Ｆ　（ｘｏ　　（ｔ）−Ｘ７（ｔ））は、ｘｏ（ｔ）の
フィルタ出力Ｆ〔ｘａ（ｔ））とＸＩ　　（ｔ）のフィ
ルタ出力Ｆ（ｘ）（ｆ）〕との差Ｆ　（ｘｏ　　（ｆ）
　）　−Ｆ　（ｘｔ　　（ｔ）〕に等しくなる。従って
、特許請求の範囲第２項の構成の装置でも、第１項の構
成の装置と同じ効果が得られる事になる。第２項で示さ
れる装置の構成例を第６図に示す。第６図の各構成要素
の番号は、第３図の対応する構成要素と同一の番号を付
してあり、説明は省略する。

特許請求の範囲第３項は、ボンピング波として正しい平
面波が得られない場合についての対策を述べたものであ
る。一般に、完全な平面波を得るには無限に広い平板振
動子を用いる必要があり、これは実際上不可能である。

近似的に平面波をｆ４る方法としは、振動子面から充分
に離れた所で充分に曲率半径の大きくなった球面波の一
部を使う事もあるが、ボンピング波として用いるには、
単位面積当りのエネルギー密度が小さくなり過ぎるので
必らずしも適当でない。従って、充分な音圧を持ったボ
ンピング波は、実際には平面波とはならず、従ってその
音圧は、時間ｔだけでなく座標２によっても変化する事
になり、平面波からのズレが無視できない位大きくなっ
た場合にはその補正が必要となる。この補正は以下の様
にして行う事ができる。

先ず、ボンピング波の音圧分布を、時間も及び座標の関
数Ｐ（ｚ、ｔ）として測定波ビームの線上で予め測定し
ておく。このＰ　（ｚ、ｔ）から、ｔ＝Ｑでｚ＝ｚｏに
あった測定波が、測定波ビーム線上で受ける音圧変化を
２　をパラメータとして２の関数ｇｉｏ（ｚ）として求
めておく事ができするものとなる。これを用いて、ｔ＝
ＱでＺ−Ｚ。

にあった測定波に対する重み関数（式（１１）のｇ（Ｚ
ｏ−ｚ）に相当する）を平面波の場合と同様にして求め
ておく事ができる。従って、この重み関数の逆特性を持
つフィルタを平面波の場合と同様。（′）を得る事がで
きる・但し・平面波の場合と異なり、測定波の注目して
いる部分がｔ＝Ｑで２軸上のどの部分にあるかによって
、つまりＺ。

Ｏ値によって、その注目している部分が２軸上で受ける
音圧分布が異ってくる。この様子を第７図に模式的に示
す。なお、第７図（Ａ）はｔ＝Ｑの時点に球面波が到来
している場合の音圧分布を表わし、第７図（Ｂ）はｔ＝
ｏにおいてｚ−２１にあった測定波に対して本来の平面
波のボンピング波を受けた場合の２軸上での音圧分布、
第７図（Ｃ）はｔ＝Ｑにおいてｚ＝ｚ１にあった測定波
が球面波のボンピング波を受けた場合の２軸−ヒで０音
圧分布、第７図（Ｄ）はｔ＝Ｑにおいてｚ＝Ｑにあった
測定波が球面波のボンピング波を受けた場合の２軸上で
の音圧分布をそれぞれ表わしている。従って、式（１１
ンのｇ’（ｚ）に相当する重み関夕も、２６をパラメー
タとして異なったものとなる。つまり、このフィルタは
、どの２座標の変化するものとなるが、第５図に示す如
く、Ｚ。

の各値に対して注目している測定波が２軸上で受ける音
圧を求める事ができるので、これがら各Ｚ。

の値に対するフィルタ特性を予め計算しておく事ができ
る。第３図における位相検出器８の出力を時間的に見る
と、これは、逐次異なるＺｏに対応する位相変化が出力
されているから、この出力に対し、各Ｚｏに対する上述
のフィルタを順次作用−（Ｚ）の分布を得る事ができる
。この時、順次異なる特性のフィルタを作用させる方法
としては、例えば、位相検出器出力を一旦Ａ／Ｄ変換し
た上でディジタルフィルタリングを施しても良い。

特別な場合として、測定波ビームがボンピング波の遠距
離音場にある場合を考えると、ボンピング波が球面波で
あるので、測定波ビーム上での音圧分布は Ｐ（ｗｔ−に人工＋ｙ”　＋ｚ’　） と書くことができる。ここでｙ＞ｘ、ｚであり、ポンプ
波の十分遠距離音場であるとすると、ビーム上（ｙ＝ｙ
ｏ）で Ｐ（ｗｔ−ｋｒ石１フー） ・　・　・　（１６２ 第１項は（ｘ、　ｙ、　　ｚ）　＝　（０，ｙｏ　、　
　Ｏ）なる映像面の中心の位相であり、第２項は位置に
よる位相のずれを表わしている。従って今、という変数
変換を行えば、式（９）は ・・・（１８） （１７）の関係から変数をＺから２に戻してやると、う
にボンピング波と測定波ビームとをある程度距離をもた
せることにより信号処理が非常に容易になる。

の変化量等の分布のＳ／Ｎ比を向上させる手段を提供し
ている。等価非線形パラメータと音圧との積による測定
波の位相変化量が小さい時には、本来の位相変化量φ（
ｚ）に対して回路内その他で発生ずる雑音を無視できな
くなる。このため、φ（ｚ）から求めた る事になる。この様子を第８図に示す。第８図（Ａ）図
示の無雑音の位相検出器出力φ（Ｚ）に第図（Ｂ）図示
の雑音Ｎが加わる事により、第８図（Ｃ）図示の如き信
号が実際には出力される。この対策として、同一測定部
位をに回測定して第８図（Ｃ）図示のＳｌ　−３ｋの如
き雑音の加わった信号を得て、これらを同一２座標の点
毎に加算すると、各点において信号成分は振幅でに倍さ
れるが、雑音成分は電力でに倍されるに過ぎず、もし雑
音が不規則雑音であれば各点における雑音振幅は５倍さ
れるにとでよる。従ってＳ／Ｎ比は尼倍改善されて第８
図（Ｄ）図示の如き出力が得られる事になる。

この方法を用いたシステム構成例を第９図に示す。第９
図において、第３図と同一の構成要素には同一の番号を
付しており説明は省略する。第９　　　　　　′図にお
いて、２１はボンピング波の送信繰り返し周期Ｔだけ信
号を遅らせる遅延回路であり、例え　　　　　　□ばＣ
ＣＤやＢＢＤ等のアナログ的手段で実現しても良いし、
またＡ／Ｄ変換器とシフトレジスタ又はメモリとを用い
たディジタル的手段で実現して良い事は言うまでもない
。２０は加算器であり、遅延回路の種類に応じてアナロ
グ又はディジタルのいずれのタイプでも良い。第９図図
示の構成の場合、位相検出器８の出力Ｓｉが同一の２軸
座標の各点毎に加算される事は明らかであり、いわゆる
同期加算によりＳ／Ｎ比の改善を行なう事ができる。

の変化量等の２次元又は３次元分布を得る事が可能とな
る。これまでの説明で明らかな様に、特定の分布を得る
事ができるから、送信・受信振動子の対を相対位置を一
定に保ったままでＸまたは１およびｙ方向に移動させ各
Ｘ、または１およびｙ座ρ０ＣＯ ３次元分布を得る事ができる。

第１０図はその構成例で、第９図と同じ構成要素には同
一番号を付しており、説明は省略する。

第１０図において、駆動部２２例えばステンビングモー
タによって振動子４．６の対を移動させる部分であり、
タイミング制御部１からの制御パルスに応じて振動子の
対の移動を行なう。位置検出部２３は例えばステンビン
グモータの軸に取り付けられたロータリ・エンコーダ等
により振動子４゜６の対の位置を検出する部分である。

メモリ制御部２４は、タイミング制御部１からのボンピ
ング波送信同期信号やクロック等を受けて位置検出部２
３からの出力に対応したメモリアドレスを発生する。メ
モリ２５は、メモリ制御部２４からの書込／読出等の制
御信号及びアドレスに従って、減液化量をＡ／Ｄ変換し
た後、記憶する。表示部２の分布をメモリから読め出し
て表示する。

化量の分布の２次元又は３次元分布を得る事ができる事
は明らかでる。

尚、第１０図において、振動子４，６の対の移動は機械
的に行なうものとしたが、振動子の対として、いわゆる
電子スキャン・プローブの対を用いる事により、機械的
走査を行なわなくても得られる様になる事は言うまでも
ない。又、第１０図においてメモリ２５を用いず、表示
管の残光２次元分布像を得ても良い事は言うまでもない
。

特許請求の範囲第６項は、上述の如くして得ら前後の変
化量（式（１７）参照） に然るべき係数、例えば式（１７〕の最終行のΔ′Ｆの
を乗する事によって、温度変化等へＴ　（ｚ）そのもの
を得られる事を述べたもので、前述の如く、Ｋ１は予め
求めておく事ができるので、特許請求の範囲第６項を実
現するには、例えば第１θ図の減算部１３の後にに１倍
の増幅器を設ければ良い。

Ｋ１が２の関数、従って時間の関数である場合には、こ
の増幅器を、時間と共に利得の変化するいわゆるＴＧＣ
（Ｔｉｍｅ　　Ｇａ　ｉｎ　　Ｃａｎｔ−ｒｏｌ）増幅
器とすれば良い事は言うまでもない。

特許請求の範囲第７項は、第１１図に示す如く、ρＯｃ
。

又は被測定超音波媒体の温度変化の前後のは温度変化量
へ１゛を推定しよとするものである。

これを実現するためには、第１１図の曲線に対応する信
号変換回路２７を、第１１図（Ａ）に対応する場合は例
えば第１２図に示す如く、第１１図（Ｂ）に対応する場
合は例えば第１３図に示す如く、挿入すれば良い事は明
らかである。尚、信号変換回路としては、例えば良く知
られたダイオード折線近似回路を作っても良いし、第１
４図に示ず如く、Ａ　／　Ｉ）変換器、符合変換用ＲＯ
Ｍ、Ｄ／Ａ変換器から成る回線を用いても良い事は言う
までもない。

（Ｅ）発明の効果 以上述べた如く、本発明によれば、超音波媒体内の温度
又は温度変化又はそれと密接な関係にある量の空間的分
布を、比較開部Ｑ’ｊな手段で容易に実現する事ができ
る。

【図面の簡単な説明】

第１図および第２図は本発明の詳細な説明する説明図、
第３図は本発明の一実施例構成、第４図は第３図におけ
る記憶部１２の一実施例構成、第５図は第３図図示構成
の動作を説明する説明図。 第６図は本発明の別の実施例構成、第７図は球面波のボ
ンピング波が到来する場合を説明する説明図、第８図ば
Ｓ／Ｎ比を向上する態様を説明する施例構成、第１１図
は温度または温度変化量との例を示す図、第１２図は温
度分布を求めるための一実施例構成、第１３図は温度変
化分布を求めるための一実施例構成、第１４図は第１２
図および第１３図における信号変換回路２７の一実施例
構成である。 図中、１はタイミング制御部、２発振器、３はドライバ
、４は送信用振動子、５は被測定超音波媒体、６ば受信
用振動子、７は受信増幅器、８は位相検出器、９はフィ
ルタ、１０はボンピング波相のドライバ、１１はボンピ
ング波発生用振動子。 １２は記憶部、ｉ３は減算回路、１５はＡ　／　Ｉ）変
換器、１６は■・Ｃメそり、１７ばＤ／Ａ変換器。 １８は制御部、２０は加算部、２１は遅延回路。 ２２は振動子の位置を変えるための駆動部、２３は振動
子の位置検出部、２４ばメモリ制御部、２５はメモリ、
２６は表示部を表わす。 面１・ｌ°″−；’　ｌ　Ｊ・・”　（Ｉ　Ｉｉ　ｌ’
；’　ｉ　：　′’：　：文なし）第１図 第１４図 手続補正書（６式う 昭和　５８年特＋Ｉ’ｌ’腐１第）Ｏ２Ｏ３，脣シン３
　補正をする右 事件との関係　　　　　ｔ１許出願人 ｆ王所　神奈川県用崎市中原区」−小１．１１中１０１
５番地（５２２）名称富士通株式会社

Claims (1)

1. 【特許請求の範囲】 （１１（イ）超音波媒体を透過する測定用連続超音波ビ
   ームを送信・受信する一対の超音波振動子と（ロ）該超
   音波ビームと交差し、被観察部位をカバーするだけのビ
   ーム幅を待つポンピング用超音波パルスを送信する超音
   波振動子と （ハ）受信波の位相の、ポンピング波による遷移量の、
   該測定用ビーム走査線上の空間的分布を検出する手段と
   、 （ニ）上記位相遷移量の、上記測定用ビーム走査線上の
   空間的分布出力を、ポンピング波の周波数特性の逆特性
   にほぼ等しい特性のフィルタに通す事により、該超音波
   媒体の等価非線形バラメーの空間的分布を得てこれを記
   憶する手段と、（ホ）該超音波媒体が温度変化を受けた
   後で上の空間的分布を得て、これと（ニ）で記憶したｆ
   段と、 を持つ事を特徴とする、超音波による温度測定装置。 （２）（イ）超音波媒体を透過する測定用連続超音波ビ
   ームを送信・受信する一対の超音波振動子と、（ロ）該
   超音波ビームと交差し、被観察部位をカバーするだけの
   ビーム幅を持つボンピング用超音波パルスを送信する超
   音波振動子と、（ハ）受信波の位相の、ポンピング波に
   よる遷移量の、該測定用ビーム走査線上の空間的分布を
   検出してこれを記憶する手段と、 （ニ）該超音波媒体が温度変化を受けた後で受信波の位
   相遷移量の上記測定用ビーム走査線上の空間的分布を得
   て、これと（）１）で記憶した位Ｉｌｌ遷移量の空間的
   分布との差を得る手段と、（ホ）上記位相遷移量の差の
   空間的分布出力を、ポンピング波の周波数特性の逆特性
   にほぼ等しし）特性のフィルタに通ず事により、該超音
   波媒体の温度変化に伴う変化量の空間的分布を得る手段
   と、を持つ事を特徴とする、超音波による温度測定装置
   。 （３）　　ポンピング波の音圧分布を場所および時間の
   関数として予め測定しておき、特許請求の範囲第１項（
   ニ）または第２項（ホ）の部分の処理を行う際に、該予
   め測定された音圧分布による補正を行う手段を持つ事を
   特徴とする特許請求の範囲ｍ　１　項ナイし第２項のい
   ずれかに記載の、超音波による温度測定装置。 （４）同一の測定用ビーム走査線上の同一位置に対応す
   る位相遷移量を、複数回のポンピング波の送信毎に加算
   する。いわゆる同期加算機能を持つ事を特徴とする特許
   請求の範囲第１項ないし第３項のいずれかに記載の超音
   波による温度測定装置。 （５）送信・受信用振動子対を一次元的または二次走査
   部・表示部を持つ事を特徴とする特許請求の範囲第１項
   ないし第４項のいずれかに記載の、超音波による温度測
   定装置。 て、適当な係数（必要に応じて場所により異ったもので
   もよい）を乗する事により、温度変化量の分布を得る様
   にした事を特徴とする特許請求の範囲第１項なしル第５
   項のいずれ力・に記載の、超音波による温度測定装置。 （７）　　被測定超音波媒体内の各部の温度または温度
   量を測定する事により、該被測定超音波媒体内の第５項
   のいずれかに記載の、超音波による温度測定装置。