(57)【要約】【課題】試料の対流が生じている状況下において、試料の対流により生じた位相成分を除外しつつ、静磁場における不均一磁場成分の分布を表す磁場マップを取得する。【解決手段】磁場マップの取得のために二重スピンエコーパルスシーケンスが利用される。二重スピンエコーパルスシーケンスは、９０°パルス７２、位相展開期間７３、勾配磁場印加期間７４、第１の１８０°パルス７８、勾配磁場印加期間８０、第２の１８０°パルス８４、受信期間８８及び勾配磁場印加期間８６を有する。第１の１８０°パルス７８の照射後のスピンエコー発生時点で残留する位相成分（試料の対流により生じる位相成分）が第２の１８０°パルス８４の照射後のスピンエコー検出時点でゼロに収束する。【選択図】図５
2019035716.tif 000002
【特許請求の範囲】
【請求項１】
静磁場における不均一磁場成分の分布を表す磁場マップを取得するための少なくとも１つのパルスシーケンスに従って、送受信を行う送受信手段と、
前記送受信の実行結果に基づいて前記磁場マップを生成する磁場マップ生成手段と、
前記磁場マップに基づいて前記静磁場に対して補正磁場を加えることにより前記静磁場を補正するシミング手段と、
を含み、
前記パルスシーケンスは、試料の対流により生じる位相成分がゼロに収束するように構成された二重スピンエコーパルスシーケンスである、
ことを特徴とするＮＭＲ測定装置。
【請求項２】
請求項１記載の装置において、
前記二重スピンエコーパルスシーケンスは、第１スピンエコー用の第１の１８０度パルス、及び、第２スピンエコー用の第２の１８０度パルス、を含む、
ことを特徴とするＮＭＲ測定装置。
【請求項３】
請求項２記載の装置において、
前記二重スピンエコーパルスシーケンスは、
横磁化を生じさせる９０度パルス、
前記９０度パルス後に設けられ、時間長δを有する第１勾配磁場印加時間、
前記第１勾配磁場印加時間後に設けられた前記第１の１８０度パルス、
前記第１の１８０度パルス後に設けられ、時間長２δを有する第２勾配磁場印加期間、
前記第２勾配磁場印加期間後に設けられた前記第２の１８０度パルス、
前記第２の１８０度パルスの後に設けられた受信期間、及び、
前記受信期間と並列に設けられた第３勾配磁場印加期間、
を含むことを特徴とするＮＭＲ測定装置。
【請求項４】
請求項２記載の装置において、
前記二重スピンエコーパルスシーケンスとして、短い位相展開時間を有する第１パルスシーケンス及び長い位相展開時間を有する第２パルスシーケンスが用意され、
前記磁場マップ生成手段は、前記第１パルスシーケンスの実行により得られた第１位相マップ及び前記第２パルスシーケンスの実行により得られた第２位相マップに基づいて前記磁場マップを生成する、
ことを特徴とするＮＭＲ測定装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明はＮＭＲ測定装置に関し、特に、不均一磁場成分の補正に関する。
【背景技術】
【０００２】
ＮＭＲ（Nuclear Magnetic Resonance）測定装置は、一般に、シム（Shim）システムを備えている。シムシステムは、静磁場に対して補正磁場を加えることにより静磁場の空間的な均一性を向上させるシミング（Shimming）のためのシステムである。勾配磁場を利用してシミングを行う手法はグラジエントシミングと呼ばれている。
【０００３】
通常、シムシステムは、シムコイル群を備えるシムコイルユニット、シムコイル群に電流群を供給するドライバ、シムコイル群に流す電流群を制御するシムシステムコントローラ等を有する。シムシステムを利用して、不均一磁場成分の空間的な分布（磁場マップ）がフラットになるように、補正磁場を規定する複数のシム値の最適な組み合わせが探索される。ここで、複数のシム値は、複数のシム項（Z1，Z2，Z3，・・・）の実体をなす複数の基本マップ（キャリブレーションマップ、シム関数）に与える複数の係数値である。より詳しくは、複数の基本マップからなる合成マップを磁場マップに加算した後において、両者の差分に相当する残差マップがそれ全体としてゼロに近付くように、複数の基本マップに与える複数のシム値が調整される。探索された複数のシム値の組み合わせから複数の電流値が定められる。なお、静磁場と平行な方向（ｚ方向）についてだけシミングを行う一次元グラジエントシミングの他、二次元グラジエントシミング及び三次元グラジエントシミングも知られている。
【０００４】
以下に、一次元グラジエントシミングを例にとり、スピンエコー法を用いた磁場マップ生成方法について説明する。
【０００５】
スピンエコー法を用いた磁場マップ生成方法においては、通常、第１パルスシーケンスと第２パルスシーケンスが利用される。例えば、第１パルスシーケンスは、短い位相展開時間（以下「shortτ」と表現する。）を有するパルスシーケンスであり、第２のパルスシーケンスは、長い位相展開時間（以下「longτ」と表現する。）を有するパルスシーケンスである。いずれのパルスシーケンスにおいても、ｚ方向の位置を特定するためにｚ方向に強さが変化する勾配磁場が利用される。
【０００６】
第１パルスシーケンスの実行結果として得られた受信信号を周波数解析して得られる複素信号から、shortτ位相マップが生成される。また、第２のパルスシーケンスの実行結果として得られた受信信号を周波数解析して得られる複素信号から、longτ位相マップが生成される。各位相マップの内訳は以下のとおりである。
shortτ位相マップ
＝shortτ＊ΔB0inhomo磁場マップ＋位相オフセットマップ  …(1)
longτ位相マップ
＝longτ＊ΔB0inhomo磁場マップ＋位相オフセットマップ   …(2)
【０００７】
上記(1)式及び上記(2)式は、それぞれ、ｚ方向の各位置における位相の内訳を示すものである。ここで、ΔB0inhomo磁場マップは、ｚ方向の各位置における不均一磁場成分（歪み成分）を示している。そのような静磁場不均一性が生じている状況下において、位相展開時間τにわたって、不均一磁場成分の大きさに依存して位相が広がる（横磁化が回転する回転座標系上でのオンレゾナンスからのずれ（位相）が増大する）。なお、位相オフセットマップは、ｚ方向の各位置における位相オフセットからなるものであり、それはｚ方向の各位置における位相の初期値とも言えるものである。そのような位相オフセットは、ＮＭＲ測定装置内での信号処理上、不可避的に生じてしまうものであり、それは、不均一磁場成分には依存せず且つτにも依存しない未知成分である。
【０００８】
longτ位相マップからshortτ位相マップを引くことにより、以下に示すΔτ位相マップが求められる。
Δτ位相マップ
＝longτ位相マップ−shortτ位相マップ
＝（longτ＊ΔB0inhomo磁場マップ＋位相オフセットマップ）
−（shortτ＊ΔB0inhomo磁場マップ＋位相オフセットマップ）
＝（longτ−shortτ）＊ΔB0inhomo磁場マップ            …(3)
【０００９】
上記のΔτ位相マップは、ｚ方向の各位置において、位相展開時間差Δτにおいて広がる位相を示すものである。ここで、Δτは（longτ−shortτ）である。上記(3)式により演算されたΔτ位相マップをΔτで割ることにより、ΔB0inhomo磁場マップ（ｚ方向の各位置における不均一磁場成分からなる磁場マップ）が求められる。もっとも、実際の装置上、一定条件下で、Δτによる除算を省略し得る。Δτ位相マップをΔB0inhomo磁場マップに相当するものとして取り扱えるからである。上記の差分演算の過程で、上記(1)式及び(2)式に含まれる位相オフセットマップがキャンセルされている。なお、上記のように求められた磁場マップに基づいて複数のシム値の最適な組み合わせが探索されることについては上述したとおりである。
【００１０】
磁場マップ生成方法としては、上記スピンエコー法の他、グラジエントエコー法も知られている。いずれの方法も、ｚ方向に傾斜した勾配磁場を利用するものである。また、いずれの方法も、通常、２種類の位相展開時間（τ）を利用するものである。
【００１１】
非特許文献１の図１には、スピンエコー法が示されている。同文献の図２には、複数のシム項（Ｚ１，Ｚ２，Ｚ３，・・・）の実体をなす複数の基本マップが示されている。非特許文献２には三次元グラジエントシミングについて記述されている。非特許文献３及び特許文献１には、二重刺激パルスシーケンス（Double-Stimulated-echo pulse sequence）が開示されている。しかし、それはシミングのためのものではない。
【先行技術文献】
【特許文献】
【００１２】
【特許文献１】特開２０１４−６６６５０号公報
【非特許文献】
【００１３】
【非特許文献１】Herve Barjat, Paul B Chilvers, Bayard K Fetler, Timothy J Horne, Gareth A Morris, A Practical Method for Automated Shimming with Normal Spectrometer Hardware, Journal of Magnetic Resonance, Volume 125, Issue 1, 1997, pp.197-201.
【非特許文献２】Guangcao Liu, Xiaobo Qu, Shuhui Cai, Zhiyong Zhang, Zhiwei Chen,  Congbo Cai, Zhong Chen, Fast 3D gradient shimming by only 2 × 2 pixels in XY plane for NMR-solution samples, Journal of Magnetic Resonance, Volume 248, November 2014, pp.13-18.
【非特許文献３】Alexej Jerschow and Norbert Muller, Suppression of Convection Artifacts in Stimulated-Echo Diffusion Experiments. Double-Stimulated-Echo Experiments, Journal of Magnetic Resonance, Volume 125, 1997, pp.372-375.
【発明の概要】
【発明が解決しようとする課題】
【００１４】
位相マップの取得のためのパルスシーケンス（第１パルスシーケンス、第２パルスシーケンス）には、位相分散（Dephase）過程と再収束（Rephase）過程とが含まれる。パルスエコー法は、位相分散過程と再収束過程との間での対称性を利用するものである。
【００１５】
しかし、試料の対流が生じていると、位相分散期間での試料状態と再収束期間での試料状態とが異なってしまうことから、スピンエコー法が正しく機能しなくなる。具体的には、測定中において測定対象である原子核がｚ方向に移動すると、位相分散期間と再収束期間とで原子核のｚ方向の位置が異なることになり、勾配磁場下において位相分散期間と再収束期間とで原子核が感じる磁場の強さが異なってしまう。このためスピンエコー法が正しく機能しなくなる。対流の影響を受けた位相マップに基づいて磁場マップを求めた場合、磁場マップにも対流の影響が現れてしまう。そのような不適正な磁場マップに基づいてシミングを行うならば、不均一磁場成分を十分にキャンセルすることができなくなる。
【００１６】
本発明の目的は、シミングの信頼性を高めることにある。あるいは、本発明の目的は、試料の対流が生じている状況下において、試料の対流により生じた位相成分が除外又は軽減された正確な磁場マップを取得することにある。
【課題を解決するための手段】
【００１７】
実施形態に係るＮＭＲ測定装置は、静磁場における不均一磁場成分の分布を表す磁場マップを取得するための少なくとも１つのパルスシーケンスに従って、送受信を行う送受信手段と、前記送受信の実行結果に基づいて前記磁場マップを生成する磁場マップ生成手段と、前記磁場マップに基づいて前記静磁場に対して補正磁場を加えることにより前記静磁場を補正するシミング手段と、を含み、前記パルスシーケンスは、試料の対流により生じる位相成分がゼロに収束するように構成された二重スピンエコーパルスシーケンスである。この構成によれば、二重スピンエコーパルスシーケンスにより、試料の対流により生じる位相成分をゼロに収束させることができ、これによって適正な位相マップひいてはシミング用の適正な磁場マップを得られる。例えば、位相展開時間が異なる２つの二重スピンエコーパルスシーケンスが順次実行される。
【００１８】
実施形態において、前記二重スピンエコーパルスシーケンスは、第１スピンエコー用の第１の１８０度パルス、及び、第２スピンエコー用の第２の１８０度パルス、を含む。第１の１８０度パルス後のスピンエコー発生時において試料の対流による位相成分が生じても、その位相成分が第２の１８０度パルス後のスピンエコー検出時においてゼロに収束する。
【００１９】
実施形態において、前記二重スピンエコーパルスシーケンスは、横磁化を生じさせる９０度パルス、前記９０度パルス後に設けられ、時間長δを有する第１勾配磁場印加時間、前記第１勾配磁場印加時間後に設けられた前記第１の１８０度パルス、前記第１の１８０度パルス後に設けられ、時間長２δを有する第２勾配磁場印加期間、前記第２勾配磁場印加期間後に設けられた前記第２の１８０度パルス、前記第２の１８０度パルスの後に設けられた受信期間、及び、前記受信期間と並列に設けられた第３勾配磁場印加期間、を含む。各１８０度パルスの前後においては、対称性をもって勾配磁場が印加される。
【００２０】
実施形態において、前記二重スピンエコーパルスシーケンスとして、短い位相展開時間を有する第１パルスシーケンス及び長い位相展開時間を有する第２パルスシーケンスが用意され、前記磁場マップ生成手段は、前記第１パルスシーケンスの実行により得られた第１位相マップ及び前記第２パルスシーケンスの実行により得られた第２位相マップに基づいて前記磁場マップを生成する。すなわち、この構成は、第１の二重スピンエコーパルスシーケンスと第２の二重スピンエコーパルスシーケンスとを実行するものである。
【発明の効果】
【００２１】
本発明によれば、シミングの信頼性を高められる。あるいは、本発明によれば、試料の対流が生じている状況下において、試料の対流に起因して生じる位相成分を除外又は軽減できる。
【図面の簡単な説明】
【００２２】
【図１】実施形態に係るＮＭＲ測定装置の構成例を示すブロック図である。
【図２】一般的なシミング方法を示すフローチャートである。
【図３】試料の対流を示す図である。
【図４】従来のスピンエコーパルスシーケンスを示す図である。
【図５】実施形態に係る二重スピンエコーパルスシーケンスを示す図である。
【図６】従来のスピンエコーパルスシーケンスを実行して得られた磁場マップを示す図である。
【図７】実施形態に係る二重スピンエコーパルスシーケンスを実行して得られた磁場マップを示す図である。
【発明を実施するための形態】
【００２３】
以下、本発明の実施形態を図面に基づいて説明する。
【００２４】
図１には、実施形態に係るＮＭＲ測定装置の構成例がブロック図として示されている。以下においては、説明簡略化のため、ｚ方向の一次元グラジエントシミングについて説明する。本発明は、二次元グラジエントシミング及び三次元グラジエントシミングに対しても適用され得る。
【００２５】
ＮＭＲ測定装置は、本実施形態において、演算制御部１０、送受信部１２及び測定部１４を有する。それらに跨ってシムシステムが構成されている。シムシステムは、シムシステムコントローラ４０、シムドライバ４３、シムコイルユニット２０の他、磁場マップを演算するための構成を含むものである。演算制御部１０は、演算部及び制御部として機能する情報処理装置等によって構成される。測定部１４は、超電導コイルを含む静磁場発生器１６、そのボア１６Ａ内に設置されるシムコイルユニット２０、及び、プローブ１８を有している。プローブ１８は、試料管２２を備えており、その内部には測定対象となる試料が入れられている。試料は溶液である。プローブ１８は、ボア１６Ａ内に挿入される挿入部と、ボア１６Ａの外部に設けられる基部と、からなる。挿入部には、試料管の他、検出コイルを含む電気回路が設けられている。
【００２６】
シムコイルユニット２０は、シミングのための複数のシムコイル（シムコイル群）を含むものである。それには、Ｚ１項に対応するＺ１シムコイル、Ｚ２項に対応するＺ２シムコイル、Ｚ３項に対応するＺ３シムコイル、等が含まれる。通常、これらのシムコイルは、それぞれ論理的なコイルであり、実際にはより多くの物理的なコイルにより構成される。シムコイルユニット２０内の特定のシムコイルを利用して勾配磁場が生成されてもよい（その場合、そのコイルはHomospiolコイルと呼ばれる）。プローブ１８内に勾配磁場生成用のコイル（ＦＧコイル）を設けてもよい。
【００２７】
測定分解能を高めるために、測定対象となる試料が存在する領域の全体にわたって、静磁場の均一度を高めることが必要であり、そのための構成がシムシステムである。プローブごとに且つ試料ごとに、必要に応じて、シミングが実行される。
【００２８】
温度可変機構４４は、試料の温度を所望温度にするための機構である。温度可変機構４４は、ボア１６Ａの上側開口又は下側開口を通じて、設定された温度を有するガスを送り込む配管４６を有する。そのガスにより試料管２２内の試料が温められ、又は、冷やされる。そのような場合、試料の対流が生じ得る。検出コイルを含む電気回路を冷却する構造及び装置を利用してもよい。そのような場合にも、試料の対流が生じ得る。試料の対流による問題とその対策については後に詳述する。
【００２９】
演算制御部１０は、送受信制御部２４、信号処理部３０、磁場マップ演算部３８、シムシステムコントローラ４０、記憶部４２等を有する。送受信制御部２４は、パルスシーケンスに従って、送信信号の生成及び受信信号の処理を制御するものである。送受信制御部２４はシーケンサとしての機能を備えている。図示されていない情報処理装置においてパルスシーケンスが生成され、それを示すデータが送受信制御部２４に与えられている（符号２５を参照）。送受信制御部２４がパルスシーケンス生成機能を備えていてもよい。
【００３０】
信号処理部３０は、受信信号を処理するモジュールであり、それは、周波数解析部３２を有している。周波数解析部３２は、複素ＦＦＴ演算により、時間軸上の複素信号（受信信号）から周波数軸上の複素信号を生成する。
【００３１】
なお、信号処理部３０は位相マップ生成手段として機能し、後述する磁場マップ演算部３８が磁場マップ生成手段として機能し、シムシステムコントローラ４０がシミング手段として機能する。
【００３２】
記憶部４２には、複数のシム項の実体をなす複数の基本マップ（Ｚ１マップ等）が格納される。各基本マップはキャリブレーション工程の実施により取得される。なお、三次元グラジエントシミングを行う場合、Ｚ１マップ〜Ｚ６マップの他、Ｘ１マップ、ＸＺ１マップ、Ｘ２マップ、ＸＺ２マップ、・・・等が記憶部４２に格納される。
【００３３】
磁場マップ演算部３８は、短い位相展開時間shortτを有する第１パルスシーケンスの実行によって得られた位相マップと、長い位相展開時間longτを有する第２パルスシーケンスの実行によって得られた位相マップと、に基づいて、磁場マップを演算するものである。第１パルスシーケンス及び第２パルスシーケンスの間では、位相展開時間τだけが相違している。本実施形態において、第１パルスシーケンス及び第２パルスシーケンスは、いずれも、対流による影響を消去することが可能な二重スピンエコーパルスシーケンスである。これに関しては後に詳述する。
【００３４】
シムシステムコントローラ４０は、シミングを制御するものであり、具体的には、磁場マップから合成マップを減算した後の残差マップが最小となるように、合成マップを構成する複数の基本マップに与える複数のシム値の最適な組合せを求めるものである。その際には、記憶部４２に記憶された複数の基本マップが参照される。通常、複数のシム値の組み合わせを変えながら、残差マップを演算及び評価する過程（つまり補正磁場調整）が繰り返される。
【００３５】
送受信部１２は、シムドライバ４３、送信部２６及び受信部２８からなる。送信部２６は、送受信制御部２４の制御の下、パルスシーケンスに従うＲＦ送信パルスを生成し、それをプローブ１８へ供給する回路である。受信部２８は、送受信制御部２４の制御の下、プローブ１８から出力されたＲＦ受信信号（ＦＩＤ（Free Induction Decay）信号）を処理する回路である。処理後の受信信号が信号処理部３０へ送られる。なお、勾配磁場を形成する際に、送受信制御部２４からシムシステムコントローラ４０に対して制御信号が与えられてもよい。
【００３６】
図２には、一般的なシミング動作が示されている。Ｓ１０Ａでは、shortτパルスシーケンス（第１パルスシーケンス）に基づく送受信が実行される。これにより、Ｓ１２Ａにおいて、shortτ位相マップが生成される。位相マップの横軸はｚ軸であり、位相マップの縦軸は位相軸である。shortτ位相マップは上記(1)式で表される。Ｓ１０Ａとは別のタイミングにおいて、Ｓ１０Ｂが実行される。Ｓ１０Ｂでは、longτパルスシーケンス（第２パルスシーケンス）に基づく送受信が実行される。これにより、Ｓ１２Ｂにおいて、longτ位相マップが生成される。longτ位相マップは上記(2)式で表される。
【００３７】
Ｓ１６では、上記(3)式に従って、longτ位相マップからshortτ位相マップを減算することにより、Δτ位相マップが演算される。Δτ位相マップをΔτで割ることにより、磁場マップが演算される。その磁場マップは、上記(3)式中のΔB0inhomo磁場マップである。このように求められた磁場マップに対して合成マップを加算した結果が、できるだけ平坦になるように、Ｓ１８において、合成マップを構成する複数の基本マップに与える複数のシム値（シム値列）が探索される。Ｓ２０では、複数のシム値に基づいて複数の電流値が演算され、それらの電流値を有する複数の電流がシムコイル群へ供給される。これにより、補正磁場が生成され、あるいは、既に加えられている補正磁場が変更される。補正磁場を可変しながら上記過程を繰り返すことによって、複数のシム値の最適な組合せが探索される。なお、Ｓ１６において、Δτ位相マップに対するΔτによる除算を省略し、Δτ位相マップをΔB0inhomo磁場マップに相当するものとして取り扱うようにしてもよい。なお、shortτを事実上、０[s]にしてもよい。
【００３８】
図３には、試料管５０が模式的に示されている。試料管５０の内部には溶液としての試料が入れられている。ｚ方向においてＮＭＲ測定対象が範囲５４である。少なくとも範囲５４の全体にわたって静磁場（磁場の向きはｚ方向に並行）を均一にすることが求められる。
【００３９】
試料管５０に対して、例えば、符号５１で模式的に示すように、温度可変用のガスが送り込まれる。すると、図３の左側に示すように温度勾配５２が生じる（横軸は温度Ｔを示しており、縦軸はｚ方向の高さを示している）。この温度勾配５２又は他の事象を原因として、試料管５０内において、試料の対流５６が生じ得る。ｚ方向における特定の高さに着目した場合、あるタイミングにおける試料状態とその後の他のタイミングにおける試料状態とが異なることになる。位相マップを取得する場合において、位相分散過程と再収束過程との間で試料状態が異なるならば、正しい位相マップを取得できなくなる。換言すれば、対流に起因する歪成分を含んだ位相マップが得られてしまう。それに基づいて磁場マップを演算するならば、その磁場マップにも対流の影響が現れてしまう。結果として、シミングにより不均一磁場成分を適正にキャンセルすることが困難となる。
【００４０】
なお、図３の右側には勾配磁場５８が示されている（横軸はｚ方向に向く磁場の強さを示しており、縦軸はｚ方向の高さを示している）。勾配磁場は、ｚ方向の各高さにおいて不均一磁場成分を識別するために、具体的には、ｚ方向の各高さにおいて共鳴周波数を異ならせるために、加えられる磁場である。符号５９は範囲５４の中央高さを示している。例えば、中央高さ５９の一方側において正の磁場が加えられ、中央高さ５９の他方側において負の磁場が加えられる。勾配磁場５８はそれ全体として線形である。
【００４１】
図４には、従来のスピンエコーパルスシーケンスが示されている。第１位相マップを取得する際に位相展開時間τとしてshortτが設定された第１スピンエコーパルスシーケンスが実行され、第２位相マップを取得する際に位相展開時間τとしてlongτが設定された第２スピンエコーパルスシーケンスが実行される。第１スピンエコーパルスシーケンスと第２スピンエコーパルスシーケンスの間では位相展開時間τの大きさだけが相違する。以下においては、それら２つのスピンエコーパルスシーケンスを区別することなく、それらの内容について具体的に説明する。なお、図４の上段には、２つの送信パルス６０，６４と、受信期間６８が示されている。図４の下段には、勾配磁場印加期間６２，６６が示されている。
【００４２】
最初に、９０°パルス６０により、横磁化が生成される。その後に位相展開期間６１が設けられている。その時間長はτ（shortτ又はlongτ）である。位相展開期間６１において、ｚ方向の各位置において、不均一磁場成分の大きさに応じて、横磁化の位相が広がる（換言すれば、回転座標系において、オンレゾナンスからのずれが拡大していく）。次の期間は、受信期間６８の半分の時間長δを有する勾配磁場印加期間６２である。勾配磁場印加期間６２においては、不均一磁場成分及び勾配磁場の２つにより位相が広がる。その後、１８０°パルス６４により、ｚ方向の各位置において位相が反転する。受信期間６８においては、その中間時点７０で、期間６１以外で広がった位相成分がゼロに収束し、スピンエコーが観測される。すなわち、その時に観測される位相は、不均一磁場成分の影響を受けて広がった位相だけとなる（以下に具体的に説明するように試料の対流がある場合にはその条件は成立しない）。受信期間６８は、勾配磁場印加期間６２の開始時点から時間長Δをおいて開始される期間である。ここではΔ＝δ＋２ｑの関係がある。受信期間６８は時間長２δを有する。受信期間６８に合わせて勾配磁場印加期間６６が設定されている。なお、周波数解析後の複素信号における実数部と虚数部とを用いて偏角を演算することにより、位相が特定される。その場合には、例えばφ＝atan2（imag, real）の計算が実行される。
【００４３】
対流が存在する場合におけるスピンエコーについて説明する。上記のようにDephase期間の時間長をδとし、Dephase期間とRephase期間の間の時間長をΔとする。流速ｖの対流が存在する場合、対流項（対流に起因する位相成分）は以下の(4)式で表される（この(4)式については上記の非特許文献３にも記載されている）。
【数１】
2019035716.tif 000003
【００４４】
上記(4)式にはｇ＊が含まれる。それはｇ＊＝ｐ（ｔ）×ｇ（ｔ）で定義される。ここでｐ（ｔ）はコヒーレンス係数であり、図４に示した例ではｐ（ｔ）は１か−１をとる。ｇ（ｔ）は勾配磁場強度である。便宜上、図４において、勾配磁場印加期間６２の開始時点をｔ＝０とする。以下においては、上記(4)式中における以下の(5)式の部分に着目する。
【数２】
2019035716.tif 000004
【００４５】
ｔ＝０からｔ＝δの期間において、上記(5)式は以下の(6)式のようになる。
【数３】
2019035716.tif 000005
【００４６】
ｔ＝δからｔ＝Δの期間において、上記(5)式は以下の(7)式のようになる。
【数４】
2019035716.tif 000006
【００４７】
ｔ＝Δからｔ＝Δ＋δの期間において、上記(5)式は以下の(8)式のようになる。
【数５】
2019035716.tif 000007
【００４８】
結局、ｔ＝０からｔ＝Δ＋δまでの全期間において、上記(5)式は以下の(9)式のようになる。
【数６】
2019035716.tif 000008
【００４９】
上記(９)式に示されるように、対流が生じている場合、対流項は０にはならない。その状態で、位相マップを得た場合、その位相マップには対流による位相成分（不均一磁場成分による位相成分とは違う位相成分）が含まれてしまうことになる。
【００５０】
図５には、上記であげた問題を解消する二重スピンエコーパルスシーケンスが示されている。そのシーケンスに含まれる位相展開時間τを異ならせることにより、第１位相マップ取得用の第１二重スピンエコーパルスシーケンスと、第２位相マップ取得用の第２二重スピンエコーパルスシーケンスと、が作成される。以下においては、それらを区別することなく、二重スピンエコーパルスシーケンスについて具体的に説明する。なお、図５の上段には、３つの送信パルス７２，７８，８４と、受信期間８８が示されている。図５の下段には、勾配磁場印加期間７４，８０，８６が示されている。
【００５１】
最初に、９０°パルス７２により、横磁化が生成される。その後に位相展開期間７３が設けられている。その時間長はτ（shortτ又はlongτ）である。位相展開期間７３において、ｚ方向の各位置において、不均一磁場成分の大きさに応じて、横磁化の位相が広がる（換言すれば、回転座標系において、オンレゾナンスからのずれが拡大していく）。次の期間は、時間長δを有する勾配磁場印加期間７４である。勾配磁場印加期間７４においては、不均一磁場成分及び勾配磁場の２つにより位相が広がる。その後、期間７６内で生じる１８０°パルス７８により、ｚ方向の各位置において位相が反転する。
【００５２】
期間７６の後、時間長２δを有する勾配磁場印加期間８０が設けられている。勾配磁場印加期間８０は前半期間８０Ａと後半期間８０Ｂとからなる。いずれも時間長δを有する。前半期間８０Ａと後半期間８０Ｂとの間で（つまり勾配磁場印加期間８０の中間時点で）、期間７３以外で広がった位相成分がゼロに収束し、スピンエコーが生じる。しかし、図４を用いて説明したように、その時点では対流項は０にならず、対流による位相成分が残留する。
【００５３】
続いて、期間８２内において、２つ目の１８０°パルス８４が照射される。これにより、ｚ方向の各位置において位相が再反転する。その後の受信期間８８と並行して勾配磁場印加期間８６が設けられている。それは、時間長δを有する前半期間８６Ａと、時間長δを有する後半期間８６Ｂと、からなる。受信期間８８の中間時点９０において、２回目のスピンエコーが生じる。その時に観測される位相は、期間７３において不均一磁場成分の影響を受けて広がった位相成分だけとなる。すなわち、１回目のスピンエコーにおいて残留していた対流による位相成分はゼロに収束する。なお、上記のように、Δ＝δ＋２ｑの関係がある。
【００５４】
二重スピンエコーパルスシーケンスにおいて、上記(5)式は以下のように表される。具体的には、ｔ＝０からｔ＝δの期間において、上記(5)式は以下の(10)式のようになる。
【数７】
2019035716.tif 000009
【００５５】
ｔ＝δからｔ＝Δの期間において、上記(5)式は以下の(11)式のようになる。
【数８】
2019035716.tif 000010
【００５６】
ｔ＝Δからｔ＝Δ＋δの期間において、上記(5)式は以下の(12)式のようになる。
【数９】
2019035716.tif 000011
【００５７】
ｔ＝Δ＋δからｔ＝Δ＋２δの期間において、上記(5)式は以下の(13)式のようになる。
【数１０】
2019035716.tif 000012
【００５８】
ｔ＝Δ＋２δからｔ＝２Δ＋δの期間において、上記(5)式は以下の(14)式のようになる。
【数１１】
2019035716.tif 000013
【００５９】
ｔ＝２Δ＋δからｔ＝２Δ＋２δの期間において、上記(5)式は以下の(15)式のようになる。
【数１２】
2019035716.tif 000014
【００６０】
結局、ｔ＝０からｔ＝２Δ＋２δまでの全期間において、上記(5)式は以下の(16)式のようになる。
【数１３】
2019035716.tif 000015
【００６１】
上記(16)式に示されるように、二重スピンエコー法によれば、対流項をゼロにすることができ、つまり、対流により生じる位相成分をゼロに収束させることが可能である。これにより、不均一磁場成分による位相成分だけを検出することが可能である。
【００６２】
２つの１８０°パルス７８，８４の一方又は両方を２つの９０°パルスで構成してもよい。スピンエコー法において１８０°パルスの代わりに２つの９０°パルスを用いる方法は刺激エコー（Stimulated Echo）法と称されている。もっとも、その方法を採用する場合、シーケンス全体の時間長が増大してしまう。位相マップ取得時間、ひいてはシミングに要する時間を短縮化するためには、２つの１８０°パルスを利用するのが望ましい。
【００６３】
磁場マップ取得に際しては通常、shortτを有する第１パルスシーケンス及びlongτを有する第２パルスシーケンスが実行されるが、上記で説明した手法は、いずれかのパルスシーケンスのみを実行する場合においても適用され得る。例えば、(1)式及び(2)式に示したshortτ位相マップ及びlongτ位相マップの相互関係を事前に特定できる場合、あるいは、それらの式に含まれる位相オフセットマップを事前に特定できる場合、１つのパルスシーケンスの実行により磁場マップを求めることが可能である。
【００６４】
図６には、図４に示したスピンエコーパルスシーケンスを利用して生成された磁場マップ（ｚ方向における不均一磁場成分の分布）９２が示されている。図７には、図５に示した二重スピンエコーパルスシーケンスを利用して生成された磁場マップ９４が示されている。各マップ９２，９４における横軸はｚ軸であり、縦軸は不均一磁場成分の強さΔB0inhomoを示している。なお、磁場マップの取得に際して、測定対象核は２Ｈであった。静磁場の強さは１４．１Ｔ（１Ｈを基準とした場合に６００ＭＨｚ）であった。
【００６５】
磁場マップ９２においては、対流項の位相成分を原因とするｚ方向に沿った大きな変動が認められるのに対し、磁場マップ９４においては、ｚ方向に沿った変動がかなり抑えられている。磁場マップ９４においては、対流の影響が除外又は軽減されている。
【００６６】
上記の手法を二次元グラジエントシミング及び三次元グラジエントシミングに適用してもよい。いずれの場合においても、ｚ方向に二重スピンエコーパルスシーケンスを適用すれば、対流による影響を除外又は軽減できる。検出コイルを冷却する場合、その影響を受けて試料に対流が生じることもある。そのような場合にも、二重スピンエコーパルスシーケンスを利用することが可能である。なお、一重スピンエコーパルスシーケンスによる測定結果と二重スピンエコーパルスシーケンスによる測定結果とを突き合わせて、対流項の実態を解明する変形例も考えられる。
【００６７】
従来のＤＯＳＹ測定においても二重スピンエコーパルスシーケンスが利用されているが、それは、シミングのためのものではない点で、また、そこには受信期間において勾配磁場の印加が存在しない点で、更に、shortτ及びlongτに対応した２種類のシーケンスを利用しない点で、実施形態に係る二重スピンエコーパルスシーケンスとは相違する。
【符号の説明】
【００６８】
１０  演算制御部、１２  送受信部、１４  測定部、１６  静磁場発生器、１８  プローブ、２０  シムコイルユニット、２２  試料管、３８  磁場マップ演算部、４０  シムシステムコントローラ、４３  シムドライバ、４４  温度可変機構。
【図１】
2019035716.tif 000016
【図２】
2019035716.tif 000017
【図３】
2019035716.tif 000018
【図４】
2019035716.tif 000019
【図５】
2019035716.tif 000020
【図６】
2019035716.tif 000021
【図７】
2019035716.tif 000022