基于磁共振弛豫时间实现加权像和压脂技术

尹朝阳,孙骏力,姚红英

(复旦大学 物理学系,上海 200433)

磁共振成像是根据生物磁性核(如氢核)在磁场中表现的共振特性进行成像的新技术. 1973年,美国科学家Paul Lauterbur首次得到核磁共振图像,随后英国科学家Peter Mansfield又进一步验证和改进了这种方法,并将其定量化,从此磁共振成像(MRI)成为核磁共振技术的一个重要分支,得到了空前的发展[1]. 随着磁体技术、超导技术、低温技术、电子技术和计算机技术等相关技术的不断进步,MRI技术得到了飞速发展,如今不仅在现代医学影像领域中大放异彩,在食品、石油测井等工业领域以及材料科学等领域都有广泛的应用. 弛豫时间是MRI中的重要物理量[2-4],利用样品的不同弛豫特性,可以通过调节序列参数实现加权成像,也可以得到不同组织成分的相关信息.

MRI基于物理现象核磁共振和它的理论基础,是物理理论和技术在医学等领域的重要应用,为此在近代物理实验加入了磁共振成像的内容,可让本科生对新技术、新仪器、新设备有所了解,开阔眼界. 此教学内容也在医学物理实验中开展,根据不同的课时设计不同的内容,本工作是学生在医学物理实验设计性、研究性实验教学模块所做的,可以让学生深入理解弛豫时间及其应用和反转恢复序列成像及加权像.

本文通过使用磁共振成像仪,进行样品的纵、横向弛豫时间的测量,并实现样品不同组织在磁共振成像过程中的区分,以及实现磁共振抑制脂肪技术.

1 原理

1.1 核磁共振

原子核系统在外磁场中发生能级分裂,在一定射频场作用下吸收能量发生能级跃迁的现象,叫做核磁共振. 其物理基础是原子核的自旋.

当核自旋系统处在z方向的外磁场Bz中时,角动量J和磁矩μ以角频率ω0=γBz绕Bz方向进动,这种进动称为拉莫尔进动,同时引起核能级的劈裂. 此时若在垂直于外磁场Bz的方向上施加一个角频率为ω的射频电磁场B1(B1<

ω=ω0=γBz

(1)

对于由大量I=1/2的微观磁矩组成的宏观物质,布洛赫提出用体磁化强度M来描述原子核系统被磁化的程度.M定义为单位体积内N个核磁矩的矢量和,即

分析M的运动,可以比较容易理解核磁共振、弛豫时间及其成像的物理内容. 如图1所示,氢原子核在外磁场B0作用下,核磁矩在上下两个圆锥面上做拉莫尔进动,对应于氢核能级的塞曼分裂高、低两个能级,形成平衡态的纵向磁矩M0. 如果在横向(比如x轴)施加射频脉冲B1,磁矩也绕x轴做拉莫尔进动,结果会产生能级跃迁,使M0向xy平面翻转. 这种沿两个轴的进动,叫做章动,也叫拉比振荡. 章动引起的M的运动轨迹变化如图2所示.

图1 氢核(I=1/2)系统在外磁场B0中核磁矩在两个塞曼能级的分布及矢量和M,平衡态时只有纵向分量M0

图2 体磁化强度M在共振跃迁过程中的运动轨迹迹,从z轴向xy平面翻转

射频脉冲的能量不同可以使M翻转到不同的位置,刚好使M翻转到xy平面的脉冲叫90°脉冲,使M翻转到负z轴方向的脉冲叫180°脉冲. 当射频脉冲停止施加,由于Bz的作用,M又会恢复到平衡态M0处,这个过程为弛豫过程,是复杂的能量交换过程. 其中自旋-自旋的相互作用可以对应到磁化强度横向分量Mxy的变化,称为横向弛豫,对应的时间为横向弛豫时间T2. 自旋-晶格的相互作用可以对应到磁化强度纵向分量Mz的变化,称为纵向弛豫,对应的时间为纵向弛豫时间T1.T2和T1这两个量及核自旋密度都对成像信号的大小有直接的影响,临床上为了更好的诊断、治疗,需要清晰的不同条件下的磁共振图像,因而需要突出某一参数的影响而弱化其他参数的影响,即加权像,比如质子密度加权像、T2加权像和T1加权像.

本文采用的核磁共振成像技术实验仪所施加的软脉冲对应的频带较窄[5],只能激发较小进动频率范围的质子,选择性较好,通常用于成像.

1.2 弛豫时间

弛豫是粒子受到激发后,以非辐射的方式回到基态而达到玻尔兹曼平衡的过程. 弛豫过程中,同时而又独立地发生两方面磁矢量变化,它的定义为90°脉冲后,z轴方向的纵向磁化矢量Mz由小恢复到M0的63%时所需的时间,称为纵向弛豫时间,反映自旋系统粒子数差从非平衡态恢复到平衡态的特征时间常数.xy平面的横向磁化矢量Mxy由大衰减到M0的37%时所需的时间,称为横向弛豫时间,即相位一致的氢质子磁矩发生相位离散,进而导致失相位,表征由于非平衡态进动位相相关产生的不为零的磁化强度横向分量恢复到平衡态时相位无关(相位随机分布)所需要的特征时间.

对于多成分样品,不同成分的弛豫特性不同,往往会有不同的纵向、横向弛豫时间,基于此可以实现磁共振成像并分析样品的组成和分布信息,以及对于特定成分进行信号抑制或加强.

1.3 SE(Spin Echo)自旋回波成像序列

自旋回波成像序列是磁共振成像的基本序列,其他序列都是在其之上发展变化而来的. 其原理如图3[5]所示.

图3 自旋回波序列时序原理图(图中横轴为时间轴,RF表示射频脉冲,Gs为选层梯度场,Gp为相位梯度场,Gr为读出梯度,也就是频率梯度场,出现回波信号时采样.)

在每个周期内,自旋回波的产生和采样过程为“90°脉冲—τ时间—180°脉冲—τ时间—采样”. 在样品线圈里感应出“自旋回波”信号,回波的幅度通常小于原始的FID信号. 原因是由于热弛豫及扰动核磁矩进动的局域场随机波动的影响,使磁化强度的幅度稍有损失,回波信号正是我们所关注和需要采集的. 信号的强度表达式为[5]

SSE(TE,TR)=AN(H)(1-e-TR/T1)e-TE/T2

(2)

其中TR是脉冲序列重复时间,TE是回波时间,A为增益,N(H)是自旋氢核密度,T1和T2分别是样品的纵向和横向弛豫时间.

1.4 IRSE(Inversion Recover Spin Echo)反转恢复自旋回波成像序列

反转恢复成像序列通常由180°X-90°X-180°YRF脉冲序列和三个正交梯度脉冲(选层、相位编码、频率编码)组成,序列时序原理如图4[5]所示.

IRSE序列先施加180°脉冲,使纵向磁化矢量Mz=+M0反转到Z轴的负方向变为Mz=-M0,然后Mz以T1时间常数进行自由弛豫向+M0恢复,经过一段时间t=TI后,在水平方向上施加90°脉冲,施加90°脉冲及之后的脉冲施加与自旋回波序列相同,都是用180°脉冲得到回波,采集回波得到成像数据. 其中第一个180°脉冲与90°脉冲之间的时间间隔记为反转时间TI.

1.5 自旋回波序列加权像

由公式2可知,通过改变脉冲序列重复时间TR和回波时间TE可以改变信号幅值,进而改变样品组织在图像上的灰度. 因此,参数加权图像即为通过选择合适的TR和TE,来实现该参数对样品组织最终信号的影响权重,以突出或者强调该参数对图像的影响.

当回波时间很短,即TE<>T1时,1-e-TR/T1≈1. 此时,T1和T2对信号的影响都近似为1,最终信号强度SSE就正比于组织所含的质子密度N(H). 这种条件下获取的图像亮度差别主要体现了质子密度的差别,故称为质子密度加权像.

当回波时间很短,即TE<

SSE=AN(H)(1-e-TR/T1)

(3)

这种条件下获取的图像亮度差别除了质子密度的影响,还受到组织T1的影响,随着TR时间的缩短,T1的影响程度增大,故称为T1加权像T1W1.

当脉冲序列重复时间TR很长,即TR>>T1时,1-e-TR/T1≈1;当回波时间不短时,T1对信号的影响近似为1,最终信号强度:

SSE=AN(H)e-TE/T2

(4)

这种条件下获取的图像亮度差除了质子密度的影响,还受到组织T2的影响,随着TE时间的延长,T2的影响程度增大,故称为T2加权像T2W1.

1.6 磁共振脂肪抑制技术

脂肪抑制技术是MRI 检查中非常重要的技术,合理利用脂肪抑制技术不仅可以明显改善图像的质量(如防止水脂肪信号重叠引起的伪影),发现高脂肪信号掩盖下的病变,从而提高病变的检出率,还可以为鉴别诊断提供重要信息[6].

根据磁共振成像原理,在一般情况下,脂肪组织在磁共振T1WI上呈很高信号,在T2WI上也呈现较高信号,即无论选择突出哪个参数进行加权成像,脂肪信号都会呈现高亮度,从而影响对其它组织的分析,此时需要对脂肪信号进行抑制.

本文通过反转恢复成像序列实现对脂肪的抑制. 在反转恢复成像序列中,设置TR>>T1,使得每采集完一个回波,Mz尽量恢复到+M0,从而不产生残余的横向磁化强度. 在施加90°脉冲之前,磁共振信号的纵向磁化分量以T1为时间常数从-M0向+M0进行弛豫,弛豫过程模拟如图5所示.

图5 反转恢复序列施加第一个180°脉冲后的弛豫曲线

其表达式可写为

SIR=AN(H)(1-2e-TI /T1)

(5)

如果信号的读取在脂肪组织的弛豫曲线过零点时进行,则脂肪对纵向磁化矢量没有贡献,无法在数据采集时产生信号,即脂肪信号在成像过程中得到抑制.TI是实现影响脂肪抑制效果的关键参数,根据公式(5),只要TI=T1ln 2,则脂肪信号在采集时为0,得到抑制和去除.

2 方法及结果

本文所用仪器为纽迈科技磁共振成像谱仪,如图6所示. 图中左侧为磁体系统(样品试管放在磁体中央);中间为电脑显示器;右侧从上往下依次为谱仪系统、主机、射频单元和梯度单元. 磁场由永磁铁提供,大小约为0.5 T,磁场均匀性在15 PPM以下,匀场范围在1 cm直径的球内. 为保证磁场的稳定性,永磁铁恒温在32 ℃,温度由程序控温表头在射频单元中显示. 实验样品为水、大豆油、肥肉、瘦肉等,放在直径15 mm的试管中,样品在试管中的高度不超过1.5 cm,免得超出匀场空间范围,引起伪影. 另外实验仪器配有相应的核磁共振分析应用软件、核磁共振成像软件.

图6 磁共振成像实验仪器

2.1 弛豫时间的测量

本文用CPMG序列测量横向弛豫时间T2;用反转恢复IR序列测量纵向弛豫时间T1.

CPMG序列是在施加一个90°射频脉冲后又施加了很多个180°射频脉冲,这样可以产生多个自旋回波,它的峰值包络线体现了纯T2的衰减规律,读出回波信号强度采用一定算法按照指数衰减规律进行反演,即可得到横向弛豫时间T2. 图7为CPMG序列示意图. 图中NECH为180°脉冲个数,P1为90°脉冲宽度,P2为180°脉冲宽度,DL1为180°脉冲和90°脉冲时间间隔,ACQ为采样时间. 图8为反演曲线示意图,由峰顶位置得到横向弛豫时间T2.

图7 CPMG序列示意图

图8 反演曲线示意图

反转恢复IR序列由180°脉冲和90°脉冲组成,IR序列脉冲时序如图9所示. 首先施加180°脉冲,使M翻转到负z轴,脉冲关闭之后M进行弛豫. 在弛豫过程中产生了纵向分量Mz,当要测量纵向分量Mz时,用90°脉冲把它翻转到xy平面,马上进行测量,该横向分量大小与Mz相等,这样就测得了Mz的值。改变IR序列中DL1的值,就得到一系列不同时刻Mz的值. 根据样品不同的纵向弛豫时间,设定不同的时间间隔DL1和采集数据个数,比如采集20个或者30个数据点,按照指数增加的规律对数据进行反演,就可以得到T1的值. 反演曲线与图8类似.采用以上方法测量得到相关成分的纵向、横向弛豫时间,如表1所示.

表1 相关成分的纵向T1、横向T2弛豫时间

图9 反转恢复序列示意图

其中,由于脂肪和瘦肉均为一整块猪肉样品中切下,为防止在脂肪/瘦肉组织中混有其他成分,采用多次取样,测量弛豫时间取平均值得到脂肪和瘦肉的纵向和横向弛豫时间,以减小实验误差.

可以发现,成分的纵向弛豫时间T1均大于横向弛豫时间T2. 因为T2是总磁化矢量M在XY平面的分量的衰减过程,当总磁化矢量弛豫到接近热力学平衡态的时候,XY平面的磁化量已经弛豫到0,质子群内部能量交换达到平衡,T2衰减已结束,而此时纵向磁化分量仍在与外界晶格进行热交换,未弛豫到最大值,即仍然在进行T1弛豫[1].

2.2 实现加权像并区分水-油样品

将水和油依次装入试管,由于水的密度比油大,选择水在下、油在上的放置方式,使两者不会混合到一起.

选择SE成像序列,按照成像的三个标准步骤1-PRESCAN、2-SOUT、3-SCAN进行成像. 选取xz层面进行成像,调节选层位置位于样品中央,层厚不要太大,防止部分容积效应而导致的MRI的层间污染,空间分辨率的下降[8]. 这是因为虽然层厚增加时,相对于较小的层厚,体素内质子数量增加,信号强度增加,图像的信噪比将会增加,图像的表观改善,但是在4 mm以下时,空间分辨率还是主要取决于部分容积效应的影响,即空间分辨率与层厚成负相关.

先做T1加权像,此时要求回波时间TE<

将成像软件中的共振频率调节为水样品的共振频率,固定TE为20 ms,改变TR分别为100、150、300、600、1 200、2 500 ms. 成像如图10所示,其中左半部分为水,右半部分为大豆油.

图10 TE=20 ms时改变TR得到的T1加权像(从左往右、从上至下依次为TR为100、150、300、600、1 200、2 500 ms)

可以发现,T1加权像中,水的亮度小于油的亮度,且随着TR的增大,水的亮度逐渐增大,而油的亮度基本不变.

对水和油的纵向弛豫过程进行模拟,得到弛豫曲线示意图如图11所示.

图11 水和油的T1弛豫曲线示意图

可以看到,在T1弛豫曲线中,大豆油的信号强度始终大于水的信号强度. 此外,由于TR的初始值100 ms已经接近大豆油的T1(152 ms),大豆油已经基本完成弛豫,故大豆油的亮度基本不发生变化;而水在该参数下还没有完成弛豫,且随着TR的增大,水的信号强度增大,则在成像上反映为亮度增强,与实验结果一致.

再做T2加权像,此时要求脉冲序列重复时间很长,即TR>>T1,弱化纵向弛豫的影响,而回波时间TE不短,体现横向弛豫的影响.

由于水的T1较长,这里设置TR远大于大豆油T1=152 ms,主要来分析大豆油的T2加权像. 先用核磁共振分析软件得到大豆油样品的共振频率,并在成像软件中设置为激发脉冲的中心频率. 固定TR为2 500 ms,改变TE分别为20、100、200、300 ms. 成像如图12所示,其中左半部分为水,右半部分为大豆油. 可以看到,随着TE的增加,大豆油部分的亮度逐渐降低直至基本为0.

图12 TR=2 500 ms时改变TE得到的T2加权像(从左至右、从上至下TE依次为20、100、200、300 ms)

对水和油的横向弛豫过程进行模拟,得到弛豫曲线示意图如图13所示.

图13 水和油的T2弛豫曲线示意图

在T2弛豫曲线中,横向磁化矢量逐渐弛豫为0,则T2加权像中,随着TE的增大,大豆油的信号强度逐渐减小直至减为0,这与实验结果一致.

然而,从横向弛豫曲线中可以看到,在T2加权像中,水的信号强度应始终大于大豆油的信号强度,且在选择的扫描参数下(TE最大为300 ms),认为水的信号强度基本不变,所以期望成像时,左半部分(水的部分)应当亮度始终为最大,且基本不变,但是从图12中可以看到,水的亮度有一定变化,甚至会大于大豆油部分的亮度.

这是因为之前设置TR时,主要关注大豆油的T1加权像,设置了TR为2 500 ms,远大于大豆油的T1,却不满足远大于水的T1,此时与T2加权像的理论公式(4)不符合. 重新回到原始公式(2),将水的T1=2 477 ms代入调节T1参数权重的项,计算得1-e-TR/T1≈0.636(TR=2 500 ms);将水的T2=1 874 ms代入调节T2参数权重的项,计算得e-TE/T2的变化范围为0.852～0.989,此时对于水,T1参数权重较大,磁共振成像不满足T2加权像的条件,T1导致的图像亮度改变需要考虑,纵向分量在弛豫中恢复的程度和横向分量在弛豫中衰减的程度二者相结合对信号强度的产生的影响,水的亮度从左起1、2、3依次变量到4又变暗.

2.3 基于IR成像序列实现脂肪抑制

对于脂肪-瘦肉样品不仅实现不同成分组织图像的区分,而进一步要求抑制脂肪信号,保留其他组织的信号. 由表1中脂肪和瘦肉的纵向、横向弛豫时间绘制T1、T2弛豫曲线如图14所示.脂肪的T1=132 ms小于瘦肉的T1=351 ms,而脂肪T2=76 ms大于瘦肉的T2=38 ms,在SE序列下,无论是T1加权像还是T2加权像,脂肪的信号都显著强于瘦肉. 所以我们采用反转恢复自旋回波序列来实现磁共振脂肪信号的抑制.

图14 脂肪-瘦肉样品纵向、横向弛豫曲线(上为纵向,下为横向)

首先将肥、廋肉样品装入试管,如图15所示. 其中上半部分为脂肪,下半部分为瘦肉.

图15 脂肪-瘦肉样品

为实现脂肪信号的抑制,使用IRSE反转恢复自旋回波成像序列成像. 由公式5,为消除脂肪信号,设置反转时间TI=T1ln 2=92 ms(这里的T1为脂肪的纵向弛豫时间132 ms),由于脂肪的T1比瘦肉的小,其纵向分量先衰减到0,此时瘦肉的纵向分量还未衰减到0,仍参与后续成像. 设置TR=2 400 ms保证其远大于T1,使得每采完一个回波,Mz尽量恢复到+M0,从而不产生残余的横向磁化强度分量. 设置TE=20 ms,使T2参数在成像中的权重接近于0. 设置在xz层面上进行扫描,按照成像的三个标准步骤1-PRESCAN、2-SOUT、3-SCAN进行成像.

测量结果如图16所示. 其中左边为IRSE序列下的成像,右边为SE序列下的成像. 在每张像中,左半部分为瘦肉组织,右半部分为脂肪组织.

图16 脂肪-瘦肉样品的磁共振成像(左为IRSE序列成像,右为SE序列成像)

可以看到,SE序列的成像中,脂肪明显亮于瘦肉,而在IRSE序列的成像中,脂肪的亮度小于瘦肉的亮度,即脂肪信号得到抑制.

但是,由图16可以发现,在IR序列成像中,虽然脂肪信号亮度较低,但是仍有一定的信号强度,此时脂肪信号没有完全被消除. 由于TI对脂肪信号的抑制起决定性作用,下面调整参数TI进行成像质量的改善.

保持TR、TE不变,而只改变反转时间TI为84 ms、88 ms、92 ms、96 ms、100 ms,分别用IR序列进行成像,结果如图17所示.

图17 IRSE序列只改变TI得到的脂肪-瘦肉样品成像(从左至右、从上至下依次为TI=84、88、92、96、100 ms)

由图17可以看到,设置TI≈84～88 ms时,脂肪信号得到最佳的抑制效果,而非理论计算出的TI=92 ms. 即理论值相较于实际的TI偏大.

对于前面推导使得TI=T1ln 2成立的公式(5),其适用条件为TR远大于T1而TE趋近于0. 而在真正的实验过程中,TR无法做到无穷大,且受实验仪器的限制TE可设置的最小值为20 ms,此时考虑这两个参数的影响,则修正后的理论公式为[5]

SIR(TI,TR,TE)=

(6)

代入实际实验时TR和TE的值,计算得修正后的TI=91.629 ms,略小于之前的理论值TI=92 ms,但是仍然大于实际值88 ms.

再考虑自旋-晶格弛豫时间T1受温度的影响. 在磁共振弛豫中,一般温度上升,T1值延长. 实验室的温度约为25 ℃,肉类样品从冷藏中取出,低于10 ℃,实验仪器样品室的温度为32 ℃,实验时对猪肉切割,分几次测量,实验也进行了几个小时之久,所以对样品的温度控制得不是很好,因而也会对猪肉成像时的实际T1值有影响,导致理论预测值与实验值有出入.

最后,仔细观察图像,可以看到,在脂肪信号的抑制过程中,瘦肉组织的成像亮度也受到些许的影响. 这是采用的核磁共振谱仪中激发的软脉冲序列频带宽度不够窄,导致除脂肪信号外,瘦肉组织中一些与脂肪信号频率相近的成分被激发参与抑制过程;以及由于IR序列的特性,某些被增强扫描的成分的T1值缩短到与脂肪组织相近,从而导致信号被抑制. 此时为实现磁共振抑制脂肪,可以采用频率选择饱和法(连续施加单一频率的预脉冲使脂肪组织中的质子被连续激发达到饱和,从而在施加真正的激发射频脉冲时不再接收能量而不产生信号)[7,8]等等其他压脂方法.

3 结论

不同成分具有不同的纵向、横向弛豫时间,一般情况下同一成分的纵向弛豫时间T1大于其横向弛豫时间T2. 基于不同弛豫时间实现水-油样品在磁共振成像中的区分,发现T1加权像要求回波时间TE<>T1而回波时间TE不很短. 在T1加权像中,随着TR的增大,样品的亮度增大直至完全弛豫;在T2加权像中,随着TE的增加,样品的亮度逐渐降低直至基本为0. 采用IRSE成像序列实现磁共振抑制脂肪技术,分析温度等参数对抑制效果的影响. 在磁共振压脂技术中,为确保成像质量,要求脉冲序列重复时间TR大于5～6个T1.