磁共振成像原理(理论)29:反转恢复序列 (Inversion-Recovery Sequence)
反转恢复序列
反转恢复是另一种常用的T1T_1T1加权成像序列,如图7.3所示。该序列以180°预备脉冲(反转脉冲)开始,简写为:
(180∘−TI−90∘−TD)N(7.20)
(180^\circ - T_I - 90^\circ - T_D)_N \tag{7.20}
(180∘−TI−90∘−TD)N(7.20)
其中TIT_ITI为反转时间,TDT_DTD为恢复时间。

第1次反转恢复:0≤t<TR0 \leq t< T_R0≤t<TR
在时间t=0−t = 0_-t=0−时刻,系统处于热平衡状态,纵向磁化强度达到最大值:
Mz(0−)=Mz0()
M_z(0_-) = M_z^0 \tag{}
Mz(0−)=Mz0()
在t=0t = 0t=0时刻施加180°射频脉冲,将纵向磁化强度完全反转:
Mz(0+)=−Mz0(A)
M_z(0_+) = -M_z^0 \tag{A}
Mz(0+)=−Mz0(A)
在反转时间TIT_ITI内,纵向磁化强度按照T1T_1T1弛豫规律从−Mz0-M_z^0−Mz0向Mz0M_z^0Mz0恢复。根据弛豫方程(式3.122)
{Mx′y′(t)=Mx′y′(0+)e−t/T2Mz′(t)=Mz0(1−e−t/T1)+Mz′(0+)e−t/T1(3.122)
\left\{
\begin{array}{l}
M_{x'y'}(t) = M_{x'y'}(0_+) e^{-t/T_2} \\
M_{z'}(t) = M_z^0 (1 - e^{-t/T_1}) + M_{z'}(0_+) e^{-t/T_1}
\end{array}
\right. \tag {3.122}
{Mx′y′(t)=Mx′y′(0+)e−t/T2Mz′(t)=Mz0(1−e−t/T1)+Mz′(0+)e−t/T1(3.122)
可得:
Mz(t)=Mz0(1−2e−t/T1)(0<t<TI)()
M_z(t) = M_z^0 \left(1 - 2e^{-t/T_1}\right) \quad (0 < t < T_I) \tag{}
Mz(t)=Mz0(1−2e−t/T1)(0<t<TI)()
在t=TI−t = T_{I-}t=TI−时刻,即90°脉冲施加前,纵向磁化强度为:
Mz(TI−)=Mz0(1−2e−TI/T1)()
M_z(T_I-) = M_z^0 \left(1 - 2e^{-T_I/T_1}\right) \tag{}
Mz(TI−)=Mz0(1−2e−TI/T1)()
在t=TIt = T_It=TI时刻施加90°射频脉冲,将此时的纵向磁化强度翻转到横向平面。对于理想的90°脉冲,在t=TI+t = T_{I+}t=TI+横向磁化强度为:
Mxy(TI+)=Mz(TI−)=Mz0(1−2e−TI/T1)(B)
M_{xy}(T_I+) = M_z(T_I-) = M_z^0 \left(1 - 2e^{-T_I/T_1}\right) \tag{B}
Mxy(TI+)=Mz(TI−)=Mz0(1−2e−TI/T1)(B)
第2次反转恢复:TR≤t<2TRT_R \leq t< 2T_RTR≤t<2TR
在实际成像序列中,我们需要考虑完整的重复周期TR=TI+TDT_R = T_I + T_DTR=TI+TD。在延迟时间TDT_DTD内,纵向磁化强度继续恢复。
在t=TR−t={T_R}_-t=TR−时刻,即第二个180°脉冲施加前,纵向磁化强度为:
Mz(TR−)=Mz0(1−e−TD/T1)()
M_z(T_R-) = M_z^0 \left(1 - e^{-T_D/T_1}\right) \tag{}
Mz(TR−)=Mz0(1−e−TD/T1)()
第二个180°脉冲将其反转为:
Mz(TR+)=−Mz0(1−e−TD/T1)(C)
M_z(T_R+) = -M_z^0 \left(1 - e^{-T_D/T_1}\right) \tag{C}
Mz(TR+)=−Mz0(1−e−TD/T1)(C)
*注意:式子A和式子C分别是第一次180°脉冲和第二次180°脉冲施加后的纵向磁化强度,两者是不一样的。
在t=(TR+TI)−t=(T_R+T_I)_-t=(TR+TI)−时刻,即第二个90°脉冲施加前,纵向磁化强度为:
Mz((TR+TI)−)=Mz0(1−e−TI/T1)+Mz′(TR+)e−TI/T1=Mz0(1−e−TI/T1)−Mz0(1−e−TD/T1)e−TI/T1=Mz0(1−2e−TI/T1)−Mz0+e−(TD+TI)/T1=Mz0(1−2e−TI/T1)+Mz0e−TR/T1()
\begin{aligned}
M_z((T_R+T_I)_-) &=M_z^0 (1 - e^{-T_I/T_1}) + M_{z'}(T_{R+}) e^{-T_I/T_1} \\
&=M_z^0 (1 - e^{-T_I/T_1}) -M_z^0 \left(1 - e^{-T_D/T_1}\right)e^{-T_I/T_1}\\
&=M_z^0 (1 - 2e^{-T_I/T_1}) -M_z^0 + e^{-(T_D+T_I)/T_1}\\
&=M_z^0 (1 - 2e^{-T_I/T_1}) +M_z^0e^{-T_R/T_1}\\
\end{aligned} \tag{}
Mz((TR+TI)−)=Mz0(1−e−TI/T1)+Mz′(TR+)e−TI/T1=Mz0(1−e−TI/T1)−Mz0(1−e−TD/T1)e−TI/T1=Mz0(1−2e−TI/T1)−Mz0+e−(TD+TI)/T1=Mz0(1−2e−TI/T1)+Mz0e−TR/T1()
在t=(TR+TI)+t=(T_R+T_I)_+t=(TR+TI)+时刻,即第二个90°脉冲施加后,横向磁化强度为:
Mxy((TR+TI)+)=Mz((TR+TI)−)=Mz0(1−2e−TI/T1)+Mz0e−TR/T1(D)
M_{xy}((T_R+T_I)_+) = M_z((T_R+T_I)_-) = M_z^0 (1 - 2e^{-T_I/T_1}) +M_z^0e^{-T_R/T_1} \tag{D}
Mxy((TR+TI)+)=Mz((TR+TI)−)=Mz0(1−2e−TI/T1)+Mz0e−TR/T1(D)
*注意:式子B和式子D分别是第一次90°脉冲和第二次90°脉冲施加后的横向磁化强度,两者是不一样的。
第3次反转恢复:2TR≤t<3TR2T_R \leq t< 3T_R2TR≤t<3TR
在t=2TR−t={2T_R}_-t=2TR−时刻,即第三个180°脉冲施加前,纵向磁化强度为:
Mz(2TR−)=Mz0(1−e−TD/T1)()
M_z(2T_R-) = M_z^0 \left(1 - e^{-T_D/T_1}\right) \tag{}
Mz(2TR−)=Mz0(1−e−TD/T1)()
第三个180°脉冲将其反转为:
Mz(2TR+)=−Mz0(1−e−TD/T1)(E)
M_z(2T_R+) = -M_z^0 \left(1 - e^{-T_D/T_1}\right) \tag{E}
Mz(2TR+)=−Mz0(1−e−TD/T1)(E)
*注意:式子C和式子E分别是第二次180°脉冲和第三次180°脉冲施加后的纵向磁化强度,两者是一样的,重复开始出现。
在t=(2TR+TI)−t=(2T_R+T_I)_-t=(2TR+TI)−时刻,即第三个90°脉冲施加前,纵向磁化强度为:
Mz((2TR+TI)−)=Mz0(1−e−TI/T1)+Mz′(2TR+)e−TI/T1=Mz0(1−e−TI/T1)−Mz0(1−e−TD/T1)e−TI/T1=Mz0(1−2e−TI/T1)−Mz0+e−(TD+TI)/T1=Mz0(1−2e−TI/T1)+Mz0e−TR/T1()
\begin{aligned}
M_z((2T_R+T_I)_-) &=M_z^0 (1 - e^{-T_I/T_1}) + M_{z'}(2T_{R+}) e^{-T_I/T_1} \\
&=M_z^0 (1 - e^{-T_I/T_1}) -M_z^0 \left(1 - e^{-T_D/T_1}\right)e^{-T_I/T_1}\\
&=M_z^0 (1 - 2e^{-T_I/T_1}) -M_z^0 + e^{-(T_D+T_I)/T_1}\\
&=M_z^0 (1 - 2e^{-T_I/T_1}) +M_z^0e^{-T_R/T_1}\\
\end{aligned} \tag{}
Mz((2TR+TI)−)=Mz0(1−e−TI/T1)+Mz′(2TR+)e−TI/T1=Mz0(1−e−TI/T1)−Mz0(1−e−TD/T1)e−TI/T1=Mz0(1−2e−TI/T1)−Mz0+e−(TD+TI)/T1=Mz0(1−2e−TI/T1)+Mz0e−TR/T1()
在t=(2TR+TI)+t=(2T_R+T_I)_+t=(2TR+TI)+时刻,即第三个90°脉冲施加后,横向磁化强度为:
Mxy((2TR+TI)+)=Mz((2TR+TI)−)=Mz0(1−2e−TI/T1)+Mz0e−TR/T1(F)
M_{xy}((2T_R+T_I)_+) = M_z((2T_R+T_I)_-) = M_z^0 (1 - 2e^{-T_I/T_1}) +M_z^0e^{-T_R/T_1} \tag{F}
Mxy((2TR+TI)+)=Mz((2TR+TI)−)=Mz0(1−2e−TI/T1)+Mz0e−TR/T1(F)
*注意:式子D和式子F分别是第二次90°脉冲和第三次90°脉冲施加后的横向磁化强度,两者是一样的,重复开始出现。
因此,我们抛弃第一次循环的数据,从第二次开始采集。
图像强度作为T1T_1T1弛豫时间和脉冲定时参数的函数可表示为:
I(r)∝ρ(r)[1−2e−TI/T1(r)+e−TR/T1(r)](7.21)
I(\mathbf{r}) \propto \rho(\mathbf{r}) \left[1 - 2e^{-T_I / T_1(\mathbf{r})} + e^{-T_R / T_1(\mathbf{r})}\right] \tag{7.21}
I(r)∝ρ(r)[1−2e−TI/T1(r)+e−TR/T1(r)](7.21)
应用举例
为了加强不同组织之间的对比度,必须使得其中一个组织信号为零,即 I(r)=0I(\mathbf{r}) = 0I(r)=0,由于自旋密度 ρ(r)≠0\rho(\mathbf{r}) \neq 0ρ(r)=0,必须有:
1−2e−TI/T1(r)+e−TR/T1(r)=0(1)
1 - 2e^{-T_I/T_1(\mathbf{r})} + e^{-T_R/T_1(\mathbf{r})} = 0 \tag{1}
1−2e−TI/T1(r)+e−TR/T1(r)=0(1)
特别地,对于 T1(r)=T10T_1(\mathbf{r}) = T_1^0T1(r)=T10 的组织,条件变为:
1−2e−TI/T10+e−TR/T10=0(2)
1 - 2e^{-T_I/T_1^0} + e^{-T_R/T_1^0} = 0 \tag{2}
1−2e−TI/T10+e−TR/T10=0(2)
与饱和恢复序列相比,反转恢复序列可通过调整TIT_ITI和TRT_RTR两个参数优化T1T_1T1对比度。特别地,通过合理选择TIT_ITI可使某些组织成分呈负强度甚至零强度。
例如,设TIT_ITI为:
TI=[ln2−ln(1+e−TR/T10)]T10(7.22)
T_I = \left[\ln 2 - \ln(1 + e^{-T_R / T_1^0})\right] T_1^0 \tag{7.22}
TI=[ln2−ln(1+e−TR/T10)]T10(7.22)
将式(7.22)中的 TIT_ITI 表达式代入指数项 e−TI/T10e^{-T_I/T_1^0}e−TI/T10:
e−TI/T10=e−1T10⋅[ln2−ln(1+e−TR/T10)]T10(3)
e^{-T_I/T_1^0} = e^{-\frac{1}{T_1^0} \cdot \left[\ln 2 - \ln\left(1 + e^{-T_R/T_1^0}\right)\right] T_1^0} \tag{3}
e−TI/T10=e−T101⋅[ln2−ln(1+e−TR/T10)]T10(3)
由于 T10T_1^0T10 在分子和分母中抵消,上式简化为:
e−TI/T10=e−[ln2−ln(1+e−TR/T10)](4)
e^{-T_I/T_1^0} = e^{-\left[\ln 2 - \ln\left(1 + e^{-T_R/T_1^0}\right)\right]} \tag{4}
e−TI/T10=e−[ln2−ln(1+e−TR/T10)](4)
应用指数函数的性质 e−(a−b)=e−a⋅ebe^{-(a-b)} = e^{-a} \cdot e^be−(a−b)=e−a⋅eb,得:
e−TI/T10=e−ln2⋅eln(1+e−TR/T10)(5)
e^{-T_I/T_1^0} = e^{-\ln 2} \cdot e^{\ln\left(1 + e^{-T_R/T_1^0}\right)} \tag{5}
e−TI/T10=e−ln2⋅eln(1+e−TR/T10)(5)
利用对数恒等式 elnx=xe^{\ln x} = xelnx=x 和 e−ln2=12e^{-\ln 2} = \frac{1}{2}e−ln2=21,可得:
e−TI/T10=12⋅(1+e−TR/T10)(6)
e^{-T_I/T_1^0} = \frac{1}{2} \cdot \left(1 + e^{-T_R/T_1^0}\right) \tag{6}
e−TI/T10=21⋅(1+e−TR/T10)(6)
将结果(6)代入信号表达式(2):
1−2⋅[12⋅(1+e−TR/T10)]+e−TR/T10=0(7)
1 - 2 \cdot \left[\frac{1}{2} \cdot \left(1 + e^{-T_R/T_1^0}\right)\right] + e^{-T_R/T_1^0} = 0 \tag{7}
1−2⋅[21⋅(1+e−TR/T10)]+e−TR/T10=0(7)
展开计算:
1−(1+e−TR/T10)+e−TR/T10=1−1−e−TR/T10+e−TR/T10=0(8)
1 - \left(1 + e^{-T_R/T_1^0}\right) + e^{-T_R/T_1^0} = 1 - 1 - e^{-T_R/T_1^0} + e^{-T_R/T_1^0} = 0 \tag{8}
1−(1+e−TR/T10)+e−TR/T10=1−1−e−TR/T10+e−TR/T10=0(8)
由此证明,当 TIT_ITI 取公式(7.22)的值时,确实有:
1−2e−TI/T10+e−TR/T10=0(7.23)
1 - 2e^{-T_I/T_1^0} + e^{-T_R/T_1^0} = 0 \tag{7.23}
1−2e−TI/T10+e−TR/T10=0(7.23)
此时,任何T1=T10T_1 = T_1^0T1=T10的组织成分在最终图像中无信号贡献,称为信号抵消效应。该特性使反转恢复序列能产生比饱和恢复序列更强的T1T_1T1对比度,特别适用于区分自旋密度和T2T_2T2值相似但T1T_1T1值略有差异的组织。
