蘇　晉，　袁小燕，　劉小利，　石春花

(1.長治醫(yī)學(xué)院 基礎(chǔ)醫(yī)學(xué)院， 山西 長治　046000;2.長治醫(yī)學(xué)院 生物醫(yī)學(xué)工程系， 山西 長治　046000)

?

不穩(wěn)定場對磁共振信號影響分析

蘇晉1，袁小燕1，劉小利1，石春花2

(1.長治醫(yī)學(xué)院 基礎(chǔ)醫(yī)學(xué)院， 山西 長治046000;2.長治醫(yī)學(xué)院 生物醫(yī)學(xué)工程系， 山西 長治046000)

摘要：利用Bloch方程，選擇適當(dāng)?shù)膮?shù)，通過數(shù)值模擬的方法詳細(xì)討論了主磁場隨時(shí)間的變化以及渦流對射頻場的擾動環(huán)境下，磁共振信號隨時(shí)間的變化情況，為臨床磁共振成像技術(shù)的應(yīng)用提供了理論參考。

關(guān)鍵詞：磁共振成像； Bloch方程； 不穩(wěn)定磁場； 磁共振信號

0引言

磁共振成像中，主磁場的穩(wěn)定性和均勻性、射頻場的穩(wěn)定性以及渦流等因素會影響到圖像處理過程，并最終影響到圖像質(zhì)量。造成磁場擾動的因素很多，除了主磁場的穩(wěn)定性之外，射頻場與植入人體內(nèi)的金屬物有相互作用，根據(jù)法拉第電磁感應(yīng)原理，變化的磁場會使得閉合的導(dǎo)體產(chǎn)生電流，從而產(chǎn)生渦流，它一方面會產(chǎn)生熱效應(yīng)[1]，但對人體影響很小，另一方面閉合線圈內(nèi)的渦流所產(chǎn)生的信號，也會對磁共振圖像帶來偽影。其中植入人體內(nèi)的金屬物品，包括外科和介入治療過程中用到人工血管、假牙、血管夾，以及一些電子植入器件如心臟起搏器、骨增長刺激器等[2-4]，之前的研究中，射頻場導(dǎo)致的渦流大多是定性的和實(shí)驗(yàn)性的討論，較少提到它對磁共振信號的具體影響。為了定量分析，以下討論了主磁場隨時(shí)擾動下，以及渦流對射頻場在時(shí)間上的擾動情況下對磁共振信號的影響。

1Bloch方程的形式和應(yīng)用

磁共振成像具有空間分辨率高、無電離輻射、多參數(shù)成像等優(yōu)點(diǎn)[5]。Bloch方程描述了在主磁場和射頻場作用下，磁化強(qiáng)度矢量隨時(shí)間演化的動力學(xué)行為，是經(jīng)典力學(xué)描述磁共振現(xiàn)象的重要理論基礎(chǔ)之一。布洛赫假設(shè)，磁化強(qiáng)度矢量M的運(yùn)動受到兩種力的支配，一種是磁力矩的作用，另一種是弛豫力，在實(shí)驗(yàn)室坐標(biāo)系下，磁化強(qiáng)度矢量M隨時(shí)變化方程如下[6]：

式中:γ----旋磁比;

M0----主磁場B0中熱平衡狀態(tài)時(shí)磁化強(qiáng)度矢量的大小。

Bloch方程的求解較為復(fù)雜，直接求解很難得到精確結(jié)果，一般采用數(shù)值求解的方法來研究具體模型。下面具體討論不穩(wěn)定磁場環(huán)境下磁共振信號的變化規(guī)律。

2穩(wěn)定均勻的主磁場環(huán)境

理想的情況為主磁場在空間和時(shí)間上都是均勻分布的。磁共振的信號是在射頻脈沖的激勵(lì)下生成，考慮在90°脈沖作用下磁共振信號隨時(shí)變化情況，具體參數(shù)選擇為：B0=1.0T，T1=1.0 s，T2=0.1 s，B1=0.001T，M0=0.02 A/m，其中B0為主磁場強(qiáng)度大小，B1為射頻脈沖磁場分量的幅值大小，T1，T2分別為縱向和橫向馳豫時(shí)間，對生物組織來說，縱向弛豫時(shí)間在100 ms到1 s，橫向馳豫時(shí)間在50～100 ms的范圍，如一個(gè)特斯拉的磁場中肌肉的橫向和縱向馳豫時(shí)間分別為750 ms和50 ms左右[7]。這里選擇T1=1.0 s，T2=0.1 s。

根據(jù)Bloch方程，通過直接求解一階微分方程組，得到磁共振信號在上述條件下隨時(shí)演化情況，如圖1所示。

圖1　90°脈沖作用下磁化強(qiáng)度矢量M在y方向投影值的變化

在90°射頻脈沖作用下，圖1顯示了直角坐標(biāo)系下磁化強(qiáng)度矢量M沿y方向的分量My(t)，它反映了磁共振信號隨時(shí)間的變化情況。顯然My(t)隨著時(shí)間的變化不斷振蕩，同時(shí)逐漸增大，由于頻率較高，所以圖像比較密集，圖中左下為t=0～1.2×10-7s時(shí)間段的局部放大圖，同時(shí)這與M由主磁場B0方向偏轉(zhuǎn)到xy平面的過程相對應(yīng),它反映了信號隨時(shí)變化趨勢。

3主磁場隨時(shí)擾動情況對磁共振信號的影響

磁共振成像中，受磁體所處的環(huán)境或者勻強(qiáng)線圈漂移的影響，磁場的穩(wěn)定會發(fā)生變化，所以，磁場的均勻性是磁共振成像的一個(gè)重要指標(biāo)。主磁場B0(設(shè)為直角坐標(biāo)系下z方向)不均勻意味著不同空間位置的自旋質(zhì)子旋進(jìn)頻率不一樣，造成了質(zhì)子的自旋磁矩之間相位不一致，從而減弱磁共振信號。在磁共振成像設(shè)備中，主磁場的不均勻性通常隨空間位置的變化被討論的較多[8-9]。磁場隨時(shí)間的變化同樣會加劇自旋磁矩的分布，影響磁共振信號的大小。為研究方便，設(shè)在射頻脈沖作用期間，隨時(shí)間周期性變化的幅度很小磁場疊加在主磁場中，并考慮其對磁共振信號的影響。

設(shè)定主磁場為:

Bz=B0+kB0cos(ωt)

式中：kB0cos(ωt)----隨時(shí)間擾動的磁場分量；

ω----擾動頻率；

k----疊加在擾動場的幅度值，通常很小，這里選擇k=0.001。

頻率選擇

式中:B1----射頻脈沖磁場分量的幅值。

磁化強(qiáng)度矢量M在xy平面投影的大小如圖2所示。

從圖中可以看到，含時(shí)微擾下的主磁場環(huán)境中，磁化強(qiáng)度矢量在xy平面的分量相對較小，這也反映出磁共振信號會在不穩(wěn)定的主磁場中出現(xiàn)信號衰減的情況。從幅度上看，這種衰減很小。在變化趨勢上，90°脈沖作用期間的中末期階段衰減較多。另外，從數(shù)值計(jì)算的結(jié)果來看，含時(shí)微擾場的頻率較大時(shí)，其對磁共振信號的影響較小，尤其是接近磁共振頻率時(shí)，幾乎沒有影響。

4RF脈沖場不均勻?qū)Υ殴舱裥盘柕挠绊?/p>

在一般的情況下，射頻場的不均勻和渦流的產(chǎn)生也會影響到磁共振信號的穩(wěn)定性，如在高場和超高場的環(huán)境下[10]。射頻場的不均勻被討論的較多[11-12]，而渦流的產(chǎn)生對磁共振信號的影響較少被提及，為了定量討論其影響，設(shè)定磁場環(huán)境為[13]：

主磁場:

脈沖射頻場:

它表示在磁共振成像區(qū)域中，射頻脈沖作用時(shí)，由于金屬物體的存在，導(dǎo)致渦流的出現(xiàn)，同時(shí)也激勵(lì)出與射頻脈沖頻率相同的電磁波分量。為了討論方便，不考慮主磁場和射頻場的均勻性問題，僅對射頻擾動場的影響進(jìn)行討論，其它參數(shù)不變。設(shè)u1=9.1×10-4,u2=1.1×10-3,v1=5×10-4,v2=8.3×10-4。

磁化強(qiáng)度矢量M在xy平面分量的模的平方，在90°射頻脈沖作用隨時(shí)間的變化趨勢如圖3所示。

圖3　渦流對射頻場擾動下|Mxy|2隨時(shí)間的變化

從數(shù)值模擬結(jié)果來看，90°脈沖作用期間，|Mxy|2會逐漸增加，渦流對射頻場的擾動(虛線所示)導(dǎo)致|Mxy|2變化率也增加，比同一時(shí)間下的穩(wěn)定射頻場環(huán)境也有了更大的|Mxy|2值。而在射頻脈沖結(jié)束前的一小段時(shí)間開始逐漸減小，90°脈沖結(jié)束時(shí)，|Mxy|2變小。這表明在選擇脈沖序列產(chǎn)生磁共振信號時(shí)，如果考慮到渦流的影響，相同的脈寬作用下，磁共振信號會有所不同，這個(gè)結(jié)論為初始磁化強(qiáng)度矢量的正確選擇提供了理論參考。

5結(jié)語

磁共振信號在不穩(wěn)定的磁場環(huán)境中，主要考慮主磁場不穩(wěn)定以及渦流對射頻場不均勻的影響，分析了其在90°脈沖作用下，磁化強(qiáng)度矢量的變化情況。結(jié)果顯示，一方面，主磁場的不穩(wěn)定將會減弱磁化強(qiáng)度矢量在xy平面的分量，從而減弱磁共振信號，但這種影響較小，因?yàn)閷?shí)際情況中擾動場幅度很小，同時(shí)從數(shù)值模擬的結(jié)果來看，擾動場的頻率較高時(shí)，比如接近磁共振頻率范圍時(shí)，幾乎沒有影響。另一方面，渦流對射頻場的影響導(dǎo)致相同的脈寬作用下，會產(chǎn)生不同的Mxy分量，從而得到不同的磁共振信號的初始值，如在脈沖作用期間Mxy會增大，而在脈沖作用末期會減小。這些變化規(guī)律為臨床磁共振成像終脈沖序列的設(shè)計(jì)和應(yīng)用提供了理論參考。

參考文獻(xiàn)：

[1]馬超,褚旭,李燁,等.磁共振成像中渦流效應(yīng)的仿真[J].清華大學(xué)學(xué)報(bào):自然科學(xué)版,2007,47(4):457-461.

[2]賀李,余深平,孟悛非.金屬物對磁共振檢查的影響及相關(guān)研究進(jìn)展[J].國際醫(yī)學(xué)放射學(xué)雜志,2009,32(1):53-57.

[3]楊長偉,李明,魏顯招.矯形外科金屬植入物磁共振成像偽影分析專用水模:CN202505350U[P].2012.

[4]王傳兵,鄒月芬,儲斌.磁共振場強(qiáng)及序列選擇對金屬植入物偽影大小影響探討[J].生物醫(yī)學(xué)工程與臨床,2014(1):5-9.

[5]潘洪凱.新型磁共振成像造影劑的合成及性質(zhì)研究[D].長春:長春工業(yè)大學(xué),2012.

[6]吉強(qiáng).醫(yī)學(xué)影像物理學(xué)[M].北京:人民衛(wèi)生出版社,2010:82.

[7]Hendee W R, Ritenour E R. Medical imaging physics[J]. Fourth Edition Medical Physics,2003,30(4):730-730.

[8]駱國麗,肖亮.主磁場不均勻情況下的磁共振虛擬成像初步研究[J].磁共振成像,2015(6):471-475.

[9]彭璽.參考制導(dǎo)的快速磁共振成像和場不均勻性修正[D].武漢:武漢大學(xué),2012:74-80.

[10]Yulin, Robert, Grossman,等.超高場MRI:腦部成像結(jié)果探討[J].磁共振成像,2010(3):169-176.

[11]蔣瑞瑞.磁共振脈沖優(yōu)化設(shè)計(jì)[D].長春:中國科學(xué)院研究生院(長春光學(xué)精密機(jī)械與物理研究所),2013:20-21.

[12]邵汀汀.高場磁共振下并行激發(fā)技術(shù)的相關(guān)研究[D].杭州:浙江大學(xué),2012:10.

[13]Balac S, Chupin L. Fast approximate solution of Bloch equation for simulation of RF artifacts in magnetic resonance imaging[J]. Mathematical & Computer Modelling,2008,48(11/12):1901-1913.

Influence of unstable field on MRI signal

SU Jin1,YUAN Xiaoyan1,LIU Xiaoli1,SHI Chunhua2

(1.College of Basic Medical, Changzhi Medical College, Changzhi 046000, China;2.Department of Biomedicine Engineering, Changzhi Medical College, Changzhi 046000, China)

Abstract：Bloch equation is applied to study the change of MRI signal with time by means of numerical simulation, considering main magnetic field variation and eddy induced RF field. The results can offer some theoretic reference for MRI application.

Key words：MRI; Bloch equation; unstability field; MRI signal.

收稿日期：2016-01-25

基金項(xiàng)目：山西省高?？萍奸_發(fā)項(xiàng)目(20091025)； 2015年長治醫(yī)學(xué)院博士啟動基金項(xiàng)目(BS15015); 長治醫(yī)學(xué)院科研啟動基金項(xiàng)目(QDZ201530)

作者簡介：蘇晉(1980-)，男，漢族，山西長治人，長治醫(yī)學(xué)院講師,碩士，主要從事理論物理方向研究,E-mail:littlesujin@163.com.

DOI:10.15923/j.cnki.cn22-1382/t.2016.3.13

中圖分類號：O 571.22

文獻(xiàn)標(biāo)志碼：A

文章編號：1674-1374(2016)03-0273-04