(57)【要約】【課題】冷凍機の稼動音が大きい環境下でもＭＲＩ装置に組み込まれたコイルから発生する音の到着時刻を正確に計測して、音源位置を推定する。【解決手段】ＭＲＩ装置を稼動する前の試験で、冷凍機７が稼動中、かつ、コイル５、６が非励磁状態のときに、冷凍機音用のセンサ１１で計測した音響信号とコイル音用センサ１２(１)〜(ｎ)で計測した音響信号から、冷凍機音の伝達関数を予め算出しておく。ＭＲＩ装置の稼動中には、冷凍機音用センサ１１で計測した冷凍機音に、予め算出しておいた伝達関数を乗じて、コイル音計測用センサ位置１２(１)〜(ｎ)における冷凍機ノイズを推定し、コイル音計測用センサ１２(１)〜(ｎ)で受信した音響信号から推定ノイズを減ずることによってコイルで発生した音を明瞭化してから、音響信号の到着時刻を正確に計測し、音源位置を推定する。【選択図】図１
2007044348.tif 000002
【特許請求の範囲】
【請求項１】
冷凍機を有する磁気共鳴イメージング装置内に配置されたコイルから発生する音の音源位置を推定する音源位置推定方法において、
冷凍機から発生する音を受信する冷凍機音検出手段及びコイルから発生する音を受信する複数のコイル音検出手段を磁気共鳴イメージング装置の所定位置に配置し、
冷凍機稼動時かつコイル非励磁時に前記冷凍機音検出手段及びコイル音検出手段でそれぞれ音響信号を収録し、
収録された音響信号に基づいて前記冷凍機から発生する音の伝播音の伝達関数を算出し、
冷凍機稼動時かつコイル励磁時に前記コイル音検出手段によって受信した実測音から冷凍機音と前記伝達関数に基づいて推定した冷凍機ノイズを除去したコイル音の受信時刻を算出し、
コイルで発生した音の音源位置を推定すること
を特徴とする音源発生位置推定方法。
【請求項２】
冷凍機を有する磁気共鳴イメージング装置内に配置されたコイルから発生する音の音源位置を推定する音源位置推定方法において、
前記冷凍機が稼動中かつ前記コイル非励磁状態のときに、冷凍機音検出手段で受信した音響信号と複数のコイル音用検出手段で受信した音響信号から前記冷凍機から発生する音の伝播音の伝達関数を予め算出しておき、
前記冷凍機が稼動中かつ前記コイル励磁状態のときに、前記冷凍機音検出手段で受信した冷凍機音と前記予め算出しておいた伝達関数とから前記コイル音用検出手段の設置位置における冷凍機ノイズを推定し、
前記コイル音用検出手段で受信した音響信号と前記推定された冷凍機ノイズに基づいて前記コイルで発生したコイル音を明瞭化し、
前記明瞭化したコイル音の到着時刻を計測し、
前記到着時刻に基づいてコイル音の音源位置を推定すること
を特徴とする音源位置推定方法。
【請求項３】
冷凍機を有する磁気共鳴イメージング装置内に配置されたコイルから発生する音の音源位置を推定する音源位置推定方法において、
前記冷凍機が稼動中かつ前記コイル非励磁状態のときに、冷凍機音検出手段で受信した音響信号と複数のコイル音用検出手段で受信した音響信号から前記冷凍機音の伝達関数を予め算出しておき、
前記冷凍機が稼動中かつ前記コイル励磁状態のときに、前記冷凍機音検出手段で受信された冷凍機音に前記予め算出しておいた伝達関数を乗じて前記コイル音用検出手段の設置位置における冷凍機ノイズを推定し、
前記コイル音用検出手段で受信した音響信号から前記コイル音用検出手段の設置位置における冷凍機ノイズを減じてコイル音を推定し、
前記推定されたコイル音の到着時刻を計測し、
前記到着時刻に基づいてコイル音の音源位置を推定すること
を特徴とする音源位置推定方法。
【請求項４】
前記伝達関数は各コイル音検出手段毎に算出されることを特徴とする請求項１ないし３のいずれか１項に記載の音源位置推定方法。
【請求項５】
前記音源位置の推定は、前記コイル音検出手段で受信した各コイル音の到着時刻の時刻差に基づいて行われることを特徴とする請求項１ないし４のいずれか１項に記載の音源位置推定方法。
【請求項６】
前記冷凍機は液体へリムを冷却し、前記各音の受信は液体ヘリウムの充填状態で実行されることを特徴とする請求項１ないし５のいずれか１項に記載の音源位置推定方法。
【請求項７】
冷凍機を有する磁気共鳴イメージング装置内に配置されたコイルから発生する音の音源位置を推定する音源位置推定装置において、
前記冷凍機から発生した音を受信する冷凍機音検出手段と、
前記コイルから発生した音を受信する複数のコイル音検出手段と、
前記各検出手段で受信した音響波をデジタル信号として記録する波形収録手段と、
前記冷凍機が稼動中かつ前記コイル非励磁状態のときに前記波形収録部で記録した冷凍機音と前記コイルからの伝播音の伝達関数を算出する伝達関数算出手段と、
前記冷凍機音と前記算出された前記伝達関数とから冷凍機ノイズを推定する冷凍機ノイズ推定手段と、
前記コイル音検出手段で検出した実測音から冷凍機ノイズを差し引いてコイル音を求めるコイル音算出手段と、
前記コイル音算出手段によって算出され、冷凍機ノイズが除去された各コイル音の受信時刻を測定する受信時刻測定手段と、
前記受信時刻測定手段によって測定された各コイル音の受信時刻差から音源位置を推定する音源位置推定手段と、
を備えていることを特徴とする音源位置推定装置。
【請求項８】
冷凍機を有する磁気共鳴イメージング装置内に配置されたコイルから発生する音の音源位置を推定する音源位置推定装置において、
前記冷凍機から発生した音を受信する冷凍機音検出手段と、
前記コイルから発生した音を受信する複数のコイル音検出手段と、
前記各検出手段で受信した音響波をデジタル信号として記録する波形収録手段と、
前記冷凍機が稼動中かつ前記コイル非励磁状態のときに前記波形収録部で記録した冷凍機音と前記コイルからの伝播音の伝達関数を算出する伝達関数算出手段と、
前記伝達関数算出手段によって算出された伝達関数を記憶する記憶手段と、
前記冷凍機音と前記記憶手段に記憶された前記伝達関数とから冷凍機ノイズを推定する冷凍機ノイズ推定手段と、
前記コイル音検出手段で検出した実測音から前記冷凍機ノイズ推定手段によって推定された冷凍機ノイズを差し引いてコイル音を求めるコイル音算出手段と、
前記コイル音算出手段によって算出され、冷凍機ノイズが除去された各コイル音の受信時刻を測定する受信時刻測定手段と、
前記受信時刻測定手段によって測定された各コイル音の受信時刻差から音源位置を推定する音源位置推定手段と、
を備えていることを特徴とする音源位置推定装置。
【請求項９】
コイルが非励磁状態で冷凍機が稼動状態にある試験モード及びコイルが励磁状態で冷凍機が稼働状態にある稼動モードのいずれかを選択するモード選択手段を備えていることを特徴とする請求項７又は８記載の音源位置推定装置。
【請求項１０】
前記モード選択手段が前記モードをいずれかに切り替える切替スイッチからなることを特徴とする請求項９記載の音源位置推定装置。
【請求項１１】
前記切替スイッチが前記波形収録手段に設けられていることを特徴とする請求項１０記載の音源位置推定装置。
【請求項１２】
前記冷凍機音検出手段が冷凍機本体に取り付けられていることを特徴とする請求項７又は８記載の音源位置推定装置。
【請求項１３】
前記コイル音検出手段が前記コイルを胴部に対して支持する支持部材の前記胴部への取り付け位置近傍に設けられていることを特徴とする請求項７又は８記載の音源位置推定装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、複数の音検出手段で受信した信号の到着時刻から、その音の発生源の位置を推定する音源位置推定装置に関するものであり、特に磁気共鳴イメージング装置（以下、ＭＲＩ装置と称す）に組み込まれている超電導コイルから発生する音を対象とした磁気共鳴イメージング装置の音源位置推定方法及び音源位置推定装置に関する。
【背景技術】
【０００２】
複数の音響信号検出センサ（以下、単にセンサと称す）で音響信号を受信し、受信信号の到着時刻差から音源位置を推定する技術としては、例えば、特許文献１に開示された発明が知られている。この発明では、音の伝播経路全体の音速を測定できない場合に、既知の位置で故意に発生させた音を、前記複数のセンサで受信し、その受信信号の到着時刻から音源位置を演算し、その演算された音源位置と前記既知の位置との誤差を格納し、前記複数のセンサで受信した受信信号の到着時刻に基づいて演算された音源位置を、前記格納された誤差に基づき補正する補正手段とを設けたもので、既知の音源位置で発生させた音の測定結果を用いて経路全体の平均音速を補正することにより、未知の音源位置の推定誤差を小さくすることを可能としている。
【特許文献１】特開平６−３２４１３９号公報
【発明の開示】
【発明が解決しようとする課題】
【０００３】
前記特許文献１記載の発明は、音の伝播経路内で音速が変化する場合があることを想定しているが、音の到着時刻は明瞭に計測できることを前提としている。したがって、音源位置を推定したい着目音に雑音が重畳し、着目音のセンサへの到着時刻が明瞭に計測できない場合には、当然、未知の音源位置を推定することはできない。すなわち、前記特許文献１記載の発明では、着目音に雑音が重畳し、着目音のセンサへの到着時刻が明瞭に計測できない場合のことは想定していない。このような状況は、例えば、冷凍機を備えているＭＲＩ装置において発生する。ＭＲＩ装置において、当該ＭＲＩ装置に組み込まれたコイルから発生する音を受信しようとする場合、冷凍機稼動音のレベルがコイルから発生する音より大きくなってしまう場合があり、冷凍機稼働音のレベルがノイズとして着目音に重畳されることになる。このようなノイズの大きい環境で着目音の到着時刻を正確に計測する方法については、特許文献１は考慮していない。
【０００４】
本発明の目的は、上記従来技術の現状に鑑み、冷凍機の稼動音が大きい環境の下でもＭＲＩ装置に組み込まれたコイルから発生する音の到着時刻を正確に計測し、この正確な到着時刻に基づいてコイルの音源位置を推定する方法及び音源位置推定装置を提供することにある。
【課題を解決するための手段】
【０００５】
本発明では、ＭＲＩ装置のコイルが非励磁状態で冷凍機が稼動状態にあるときを試験モード、コイルが励磁状態にあるときを稼動モードと定義する。試験モードは、コイルからは音が発生せず、冷凍機音は発生している状況となる。稼動モードでは、コイルから音が発生する可能性がある。本発明は、試験モードにおいて、冷凍機音用検出手段で計測した冷凍機音とコイル音検出手段で計測した冷凍機ノイズをＭＲＩ装置の稼動前に収録し、稼動モードでＭＲＩ装置を操作する前に、伝達関数を算出・記録し、この記録された伝達関数を使用して冷凍機ノイズを推定し、更に、推定された冷凍機ノイズを使用してコイル音検出手段によって検出される音響信号から冷凍機ノイズを除去し、除去したコイル音の受信時刻から音源位置を推定するようにしたものである。
【０００６】
具体的には、本発明は、冷凍機を有する磁気共鳴イメージング装置内に配置されたコイルから発生する音の音源位置を推定する音源位置推定方法において、前記冷凍機が稼動中かつ前記コイル非励磁状態のときに、冷凍機音検出手段で受信した音響信号と複数のコイル音用検出手段で受信した音響信号から前記冷凍機音の伝達関数を予め算出しておき、前記冷凍機が稼動中かつ前記コイル励磁状態のときに、前記冷凍機音検出手段で受信した冷凍機音と前記予め算出しておいた伝達関数とから前記コイル音用検出手段の設置位置における冷凍機ノイズを推定し、前記コイル音用検出手段で受信した音響信号と前記推定された冷凍機ノイズに基づいて前記コイルで発生したコイル音を明瞭化した後、音響信号の到着時刻を計測し、前記到着時刻に基づいて音源位置を推定することを特徴とする。
【０００７】
また、本発明は、前記方法を実施する装置として、冷凍機を有する磁気共鳴イメージング装置内に配置されたコイルから発生する音の音源位置を推定する音源位置推定装置において、前記冷凍機から発生した音を受信する冷凍機音検出手段と、前記コイルから発生した音を受信する複数のコイル音検出手段と、前記各検出手段で受信した音響波をデジタル信号として記録する波形収録手段と、前記波形収録部で記録した冷凍機音と前記コイルからの伝播音の伝達関数を算出する伝達関数算出手段と、前記冷凍機音と前記算出された前記伝達関数とから冷凍機ノイズを推定する冷凍機ノイズ推定手段と、前記コイル音検出手段で検出した実測音から冷凍機ノイズを差し引いてコイル音を求めるコイル音算出手段と、前記コイル音算出手段によって算出され、冷凍機ノイズが除去された各コイル音の受信時刻を測定する受信時刻測定手段と、前記受信時刻測定手段によって測定された各コイル音の受信時刻差から音源位置を推定する音源位置推定手段とを備えていることを特徴とする。
【発明の効果】
【０００８】
本発明によれば、冷凍機の稼動音が大きい環境下でも磁気共鳴イメージング装置に組み込まれたコイルから発生する音の到着時刻を正確に計測できるので、コイルから発生した音の音源位置を正確に推定することができる。その際、冷凍機ノイズを除去し、小さなコイル音の発生まで検知できるので、コイルの状態を細かく把握することが可能になる。
【発明を実施するための最良の形態】
【０００９】
以下、本発明の実施形態について図面を参照して説明する。
【００１０】
本発明に係る音源位置推定方法は、冷凍機を有するＭＲＩ装置に組み込まれたコイルから発生する音を対象としている。図１は本発明の一実施形態に係る音源位置推定装置の構成を示すブロック図である。図１において、ＭＲＩ装置は上部胴１、下部胴２、両者を連結する連結管３，４、及び上部胴の上部に設けられた冷凍機７を備えている。上部胴１には上コイル５が、下部胴２には下コイル６がそれぞれ組み込まれている。冷凍機７には冷凍機７から発する音を検出する冷凍機音検出手段としての冷凍機用出センサ１１が、上部胴１及び下部胴２にはそれぞれに収容されたコイルから発する音を検出する複数のコイル音検出手段としてのコイル音用センサ１２(１) 〜１２(ｎ)がそれぞれ設けられている（ｎは２以上の整数、以下、同様）。これらのセンサとして、アコースティク・エミッション（ＡＥ）センサが使用される。また、アコースティク・エミッション（ＡＥ）センサに代えて超音波探触子を使用することもできる。なお、本発明では「音」は固体中を伝搬するものを対象としている。この実施形態では、センサはｎ（４≦ｎ）個設けられている。
【００１１】
また、アンプ１３、１４（１）〜１４（ｎ）には、使用センサに適したＡＥセンサ用アンプが使用され、これらのアンプ１３、１４（１）〜１４（ｎ）は波形収録部１５の前段に設けられる。そして、冷凍機用センサ１１からの検出出力がアンプ１３に、各コイル音用センサ１２（１）〜１２（ｎ）からの検出出力がアンプ１４（１）〜１４（ｎ）に入力され、アンプ１３，１４（１）〜１４（ｎ）で増幅された出力が波形収録部１５に入力される。波形収録部１５の後段には伝達関数算出部１７及びコイル音算出部２０が設けられている。コイル音算出部２０の前段には冷凍機ノイズ推定部１９が設けられ、冷凍機ノイズ推定部１９には前記伝達関数算出部１７の後段に設けられた伝達関数記憶部１８に格納された伝達関数が入力される。更に、前記コイル音算出部２０の後段には、受信時刻推定部２１及び音源位置推定部２２が設けられ、音源位置推定部２２には伝播時間記憶部２５から伝播時間が入力され、表示部２３に前記音源位置推定部２２から推定された音源位置が出力される。
【００１２】
波形収録部１５は、Ａ／Ｄボードを搭載したコンピュータなどの情報処理装置、ここではパーソナルコンピュータ（以下、パソコンと称す）、または、デジタルオシロスコープなどを使用する。波形収録部１５にはモード切替スイッチ１６が備えられている。このモード切替スイッチ１６は、多接点の機械的な切換器を使用したハード的なものと、パソコンにインストールしたプログラムで信号処理の流れを切り替えるソフト的なもので実現できる。ソフト的な方法で切り替える場合は、パソコンの表示画面を見ながらキーボードかマウスを操作してモードを選択し、意図とおりのモードを選択していることを画面表示で確認することになる。伝達関数算出部１７、冷凍機ノイズ推定部１９、コイル音算出部２０、受信時刻測定部２１、音源位置推定部２２はパソコンとパソコンにインストールしたプログラムによって実現でき、伝達関数記憶部１８、伝播時間差記憶部２５はハードディスクなどの記憶媒体によって実現できる。表示部２３には、ＣＲＴ、液晶モニタなどのパソコン画面を使用する。波形収録部１５をＡ／Ｄボードを搭載したパソコンとし、モード切替スイッチ１６をソフト的なものにすると、音響センサ１１，１２とアンプ１３，１４以外は、パソコンとＡ／Ｄボード、ハードディスクなどの周辺機器とパソコンにインストールしたプログラムなどで本発明の装置を実現できる。
【００１３】
被検査体は磁気共鳴イメージング装置（ＭＲＩ）であり、前述のようにそれぞれの内部に上コイル５と下コイル６を組み込んだ上部胴１と下部胴２が連結管３，４で接続されている。上部胴１に取り付けられた冷凍機７は、稼動中は機械音を発生する。冷凍機７には、冷凍機音を計測するためにセンサ１１が設置され、計測した音響信号をアンプ１３で増幅する。上部胴１および下部胴２の周囲には、コイル５，６からの音を計測するための複数のセンサ１２(１)〜１２(ｎ)が設置してあり、計測した音響信号をアンプ１４(１)〜１４(ｎ)で増幅する。上部胴１、下部胴２の内部には、図示しない容器が設置されており、上コイル５、下コイル６はこの容器内に取り付けられている。コイル５，６は図１には図示していない数本のコイル支持部材（図４符号８（１）〜８（ｐ））で保持されており、コイル５，６からの音を計測するためのセンサ１２(１)〜１２(ｎ)は、この支持部材と上部胴１、下部胴２の接続部近傍に取り付けてある。また、上部胴１と下部胴２には前記コイル５，６の容器を内包するように密封容器が配置され、液体ヘリウムが充填される。なお、液体ヘリウムの密封容器や密封構造は公知のものなので、説明は省略する。また、前述の支持構造材と上部胴１、下部胴２の接続部近傍とは当該接続部が音の伝達部であることから、コイル５，６側から胴部１，２への音の伝達部近傍である。従って、接続形状が異なれば、センサの取り付け位置も異なってくる。
【００１４】
冷凍機７が停止しているときには、大きな機械音ノイズがないため、コイル５，６で発生する音響信号の受信時刻を正確に評価できる。このため、正確な時刻差に基づく音源位置推定が可能である。しかしながら、冷凍機７の稼動中では、冷凍機音が発生し、上部胴１、連結管３，４、下部胴２等の構造材を伝播してコイル音用センサ１２(１)〜１２(ｎ)の設置位置まで達してノイズとなる。コイルで発生した音と冷凍機ノイズが重畳してしまうと、コイル音の受信時刻を正確に判断することができない。
【００１５】
そこで、本実施形態では、ＭＲＩ装置のコイルが非励磁状態で冷凍機が稼動状態にあるときと、コイルが励磁状態で冷凍機が稼働状態のときのそれぞれについて前記センサ１３、１４（１）〜１４（ｎ）から音響信号を検出する。この実施形態では、前者を試験モード、後者を稼動モードと定義する。試験モードでは、コイル５，６は非励磁状態であるのでコイル５，６から音は発生せず、冷凍機７から音が発生している状況となる。稼動モードでは、コイル５，６から音が発生する可能性がある。
【００１６】
そのため、冷凍機７が稼動中で上コイル５、下コイル６が非励磁中の試験モードで、冷凍機音の伝播特徴を把握するためのデータを収録して、伝達関数を算出・記憶しておく。稼動モードでは、コイル音用センサ１２(１)〜１２(ｎ)の受信信号から冷凍機ノイズを除去し、コイル音の受信時刻を測定する。これにより、正確な時刻差に基づく音源位置推定が可能となる。
【００１７】
まず、試験モードにおける信号処理方法について説明する。  
試験モードでは、モード切替スイッチ１６を試験モード（図１の矢印方向）にしておく。この状態で、センサ１１、１２(１)〜(ｎ)で冷凍機７の振動によって生じる音響信号を受信し、アンプ１３、１４(１)〜(ｎ)で増幅した信号を波形収録部１５でデジタル信号に変換して収録する。モード切替スイッチ１６が試験モードにしてあるので、収録波形は伝達関数算出部１７に送られる。
【００１８】
伝達関数算出部１７における信号処理方法を図２に模式的に示す。冷凍機音は連続音ではなく間歇的な繰り返し音の場合があるので、オペレータはセンサ１１で受信した冷凍機音波形を観察し、冷凍機音の特徴を含む信号振幅の大きい部分を時間ゲートＧＴで切り取る。また、オペレータはセンサ１２(１)で受信した伝播音Ｓ１の波形（横軸に時間、縦軸に電圧を取った波形、以下、同様）を観察し、前述の時間ゲートＧＴで切り取った部分の冷凍機波形がセンサ１２(１)位置まで伝播した部分を推定して、この部分の波形を冷凍機音と同じ時間幅のゲートＧＴで切り取る。冷凍機音Ｓ０と伝播音Ｓ１の時間ゲートＧＴの開始時刻が異なる可能性があるので、時間ゲートＧＴ開始の時刻差をむだ時間ｄ1として記録しておく。ゲートＧＴで切り取った冷凍機音波形にウェーブレット変換、または、高速フーリエ変換（ＦＥＴ）を適用し、周波数毎の振幅Ａ0と位相φ0を算出する。伝播音Ｓ１の波形も同様に変換し、周波数毎の振幅Ａ1と位相φ1を算出する。冷凍機音Ｓ０を入力信号、伝播音Ｓ１を出力信号と考え、次式に示す周波数毎の振幅比Ｇ（ｆ）と位相差θ（ｆ）で表現できる伝達関数
Ｇ１（ｆ）＝Ａ１（ｆ）／Ａ０（ｆ）                    ……  （１）
θ１（ｆ）＝φ１（ｆ）−φ０（ｆ）                    ……  （２）
を算出する。位相差θ(f)は周波数毎の信号応答の遅れを表すものなので、次式
τ１（ｆ）＝ｄ１＋φ１（ｆ）／２πｆ                  ……  （３）
を適用して、周波数毎の時間遅れτとして表現することもできる。
【００１９】
算出した伝達関数は、周波数ｆ、振幅比Ｇ１（ｆ）、位相差φ１（ｆ）、むだ時間ｄ１からなる行列（記録形式１）、または、周波数ｆ、振幅比Ｇ１（ｆ）、遅れ時間τ１（ｆ）からなる行列（記録形式２）などの形式で、伝達関数記憶部１８に記録する。
【００２０】
冷凍機７からセンサ１２（１）までの音響波伝播経路とセンサ１２（２）までの伝播経路は異なるので、伝達関数も異なる可能性がある。センサ１２（２）〜センサ１２（ｎ）で受信した伝播音Ｓ２〜Ｓｎについても、上記センサ１２（１）の場合と同様の手順で伝達関数を算出し。伝達関数記憶部１８に記録する。なお、切り取るゲートＧＴの数は図２に示す１つのゲートＧＴの波形が繰り返されるので、観察に基づいて立ち上がりから減衰まで１つのゲートＧＴで切り取れば充分である。
【００２１】
次に、稼動モードにおける信号処理について説明する。  
試験モードでは、モード切替スイッチ１６を試験モードにしておく。センサ１１、１２（１）〜（ｎ）で受信し、アンプ１３、１４（１）〜（ｎ）で増幅した信号を、波形収録部１５でデジタル信号に変換して収録する。モード切替スイッチ１６が稼動モードにしてあるので、センサ１１で計測した冷凍機音Ｓ０は冷凍機ノイズ推定部１９に、センサ１２（１）〜（ｎ）で受信した実測音Ｓ１’〜Ｓｎ’はコイル音算出部２０に送られる。冷凍機ノイズ推定部１９とコイル音算出部２０における信号処理法を図３に模式的に示す。センサ１１で計測した冷凍機音波形ＷＲ０にウェーブレット変換を適用し、ｆ１〜ｆｋの周波数成分毎の信号変化Ｗｆ１〜Ｗｆｋに分離する。伝達関数記憶部１８に記録した伝達関数の振幅比Ｇ１（ｆ）と遅れ時間τ１（ｆ）を用いて、
Ｗ’ｆ１（ｔ）＝Ｇ１（ｆ１）・Ｗｆ１（ｔ−τ１（ｆ１））      ……（４．１）
Ｗ’ｆ２（ｔ）＝Ｇ２（ｆ２）・Ｗｆ２（ｔ−τ２（ｆ２））      ……（４．２）
：                                        ：
Ｗ’ｆｋ（ｔ）＝Ｇｋ（ｆｋ）・Ｗｆｋ（ｔ−τｋ（ｆｋ））      ……（４．ｋ）
によりセンサ１２（１）位置での各周波数成分の応答Ｗ’ｆ１〜Ｗ’ｆｋを予測する。
【００２２】
各周波数成分の応答Ｗ’ｆ１〜Ｗ’ｆｋに逆ウェーブレット変換を適用し、冷凍機ノイズＮ１の予測波形Ｗ’Ｒ０を算出する。コイル音算出部２０では、波形収録部１５から送られてきた実測音Ｓ１’の波形ＷＲ１から冷凍機ノイズＮ１の予測波形Ｗ’Ｒ０を
ＷＣ１（ｔ）＝ＷＲ１（ｔ）−Ｗ’Ｒ０（ｔ）                  ……（５）
のように差し引き、コイル音Ｃ１の波形ＷＣ１を算出する。式（５）で算出されたコイル音Ｃ１の波形ＷＣ１は、実測音Ｓ１’の波形から冷凍機ノイズＮ１の波形を減じることにより冷凍機ノイズＮ１が除去されたものとなっており、コイル５で発生する音のＳＮ比が著しく向上している。センサ１２（２）〜センサ１２（ｎ）で受信した伝播音についても、上記センサ１２（１）の場合と同様の手順で冷凍機ノイズＮ２〜Ｎｎを除去し、ＳＮ比の良いコイル音Ｃ２ないしＣｎの波形ＷＣ２〜  ＷＣｎを算出する。
【００２３】
受信時刻測定部２１では、コイル音Ｃ１〜Ｃｎの波形ＷＣ１〜  ＷＣｎの受信時刻を測定し、センサ１２（１）〜（ｎ）で受信したコイル音の時間差を算出する。コイル音ＷＣ１〜  ＷＣｎは冷凍機ノイズＮ１〜Ｎｎを除去した信号なので、受信時刻を正確に測定することができる。音源位置推定部２２では、受信時刻測定部２１で測定した時間差と、予め伝播時間差記憶部２５に記録しておいた時間差データベースを比較し、コイル音の発生位置を推定する。
【００２４】
ここで、仮想音源走査法に基づいて伝播時間差記憶部２５に記録する時間差データベースの算出する方法について、図４を用いて説明する。上コイル５は、複数のコイル支持材８(１)〜８(ｐ)によって上部胴１に固定されている。上コイル５の周方向、径方向、高さ方向をそれぞれ適切な間隔で分割し、コイル５をメッシュ状に細分して考える。このメッシュの交点に仮想音源２４があると仮定し、仮想音源２４からコイル音用センサ１２(１)〜１２(ｐ)までの伝播経路を考える。上コイル５がコイル支持材８(１)〜８(ｐ)で支持されているので、仮想音源２４から発した音は、まず上コイル５とコイル支持材８(１)〜８(ｐ)の接合点５（１）〜５（ｐ）まで上コイル５の内部を伝播する。上コイル５の中心孔を経由しないルートで仮想音源２４から上コイル５とコイル支持材８(１)〜８(ｐ)の接合点５（１）〜５（ｐ）までの最短経路を考え、コイル音速の情報を加えると、接合点５（１）〜５（ｐ）までの伝播時間を算出できる。
【００２５】
次に、コイル支持材８(１)〜８(ｐ)、および、上部胴１を伝播する経路を考え、コイル支持材８(１)〜８(ｐ)中の音速、上部胴１構造材中の音速の情報を加えると、接合点５（１）〜５（ｐ）からコイル音用センサ１２までの伝播時間を算出できる。その際、上部胴１の上コイル５を収納した空間には上部胴１との間に真空層を挟んで密閉構造で液体ヘリウムが充填されているので、コイル５，６の中心部の空間部にも液体ヘリウムが満たされている。音は固体（ここでは金属）の方が伝達特性に優れており、固相から液相に伝播することは少なく、よしんば液相に伝播しても液相から固相に伝播する間に減衰するので、液相部分（液体ヘリウム）に伝達される音は実質的に無視することができる。
【００２６】
そこで、仮想音源２４から接合点５(１)〜５(ｐ)までの伝播時間と、接合点５(１)〜５(ｐ)からコイル音用センサ１２までの伝播時間を加算すると、全体の伝播時間を算出できる。もし、仮想音源２４からコイル音用センサ１２までの経路が複数考えられる場合には、考えられる経路の伝播時間を全て算出し、最短の伝播時間を選択する。各コイル音用センサ１２(１)〜(ｐ)の伝播時間から最短伝播時間を差し引くと、仮想音源２４の位置で発生した音の時間差を算出できる。上コイル５のメッシュ上で仮想音源２４を走査しながら各走査点での上記のように伝播時間差を算出する操作を繰り返し、時間差データベースを算出する。音源位置推定部２２では、受信時刻算出部２１で測定した時間差と伝播時間差記憶部２５に記録しておいた時間差データベースを比較し、時間差の差が最小となるメッシュ交点位置を音源位置と推定する。推定結果は、表示部２３に表示する。
【００２７】
このように処理すると、冷凍機ノイズを除去したＳＮ比の良いコイル音の受信時刻差に基づいて音源位置を推定しているので、確からしい音源位置推定が可能となる。
【００２８】
また、冷凍機ノイズＮを除去しない一般的な方法では、発生したコイル音が冷凍機音より小さい場合には検出は困難であり、コイル音の発生を検知することができないため音源位置の評定を行うこともできない。これに対し、本実施形態では、前記式（５）のようにしてコイル検出音から冷凍機ノイズを除去できるので、小さなコイル音の発生まで検知することが可能となり、概略位置を評定できる。その結果、コイルの状態を細かく把握することが可能になる。
【図面の簡単な説明】
【００２９】
【図１】本発明の一実施形態に係る音源位置推定装置の構成を示すブロック図である。
【図２】試験モードで収録した音響信号の伝達関数算出手順を模式的に示す説明図である。
【図３】稼動モードでコイル音信号から冷凍機ノイズを除去する信号処理手順を模式的に示す説明図である。
【図４】時間差データベースの算出方法を示す説明図である。
【符号の説明】
【００３０】
１  上部胴
２  下部胴
５  上コイル
６  下コイル
７  冷凍機
８(１)〜８(ｐ)  コイル支持材
１１  冷凍機音用センサ
１２(１)〜１２(ｎ)  コイル音用センサ
１５  波形収録部
１６  モード切替スイッチ
１７  伝達関数算出部
１８  伝達関数記憶部
１９  冷凍機ノイズ推定部
２０  コイル音算出部
２１  受信時刻算出部
２２  音源位置推定部
２３  表示部
２４  仮想音源
【図１】
2007044348.tif 000003
【図２】
2007044348.tif 000004
【図３】
2007044348.tif 000005
【図４】
2007044348.tif 000006