磁共振成像原理(理论)30:基本自旋回波成像 (Basic Spin-Echo Imaging)
本专题主要参考《Principles of Magnetic Resonance Imaging A Signal Processing Perspective 》-Sec 7.4

饱和恢复自旋回波序列
在磁共振成像中,普通的饱和恢复和反转恢复序列产生的是T1加权的自由感应衰减信号。然而,为了实现k空间的对称覆盖,通常采用回波信号进行数据采集。本节将详细讨论将自旋回波技术引入这些成像序列后所产生的图像对比度特性。

上图展示了饱和恢复自旋回波序列的时序图。该序列结合了饱和恢复的磁化准备机制和自旋回波的信号采集方式。
第1次饱和恢复:0≤t<TR0 \leq t< T_R0≤t<TR
在t=0−t = 0_-t=0−时刻,即施加90°射频脉冲 (假设是x′x'x′方向的)前,系统处于热平衡状态,纵向和横向磁化强度分别是:
Mz′(0−)=Mz0Mx′y′(0−)=0()
\begin{aligned}
M_z'(0_-) &= M_z^0 \\
M_{x'y'}(0_-) &=0
\end{aligned}\tag{}
Mz′(0−)Mx′y′(0−)=Mz0=0()
在t=0+t = 0_+t=0+时刻,即施加90°射频脉冲后,纵向和横向磁化强度分别为:
Mz′(0+)=0Mx′y′(0+)=Mz0()
\begin{aligned}
M_z'(0_+) &= 0 \\
M_{x'y'}(0_+) &= M_z^0
\end{aligned}\tag{}
Mz′(0+)Mx′y′(0+)=0=Mz0()
根据弛豫方程(式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)
在0≤t<TI0 \leq t < T_I0≤t<TI 时刻,横向和纵向磁化强度分别为:
Mz′(t)=Mz0(1−e−t/T1)(0<t<TI)Mx′y′(t)=Mz0e−t/T2(0<t<TI)()
\begin{aligned}
M_{z'}(t) &=M_z^0 \left(1 - e^{-t/T_1}\right) &\quad (0 < t < T_I) \\
M_{x'y'}(t) &= M_z^0e^{-t/T_2} &\quad (0 < t < T_I)
\end{aligned}
\tag{}
Mz′(t)Mx′y′(t)=Mz0(1−e−t/T1)=Mz0e−t/T2(0<t<TI)(0<t<TI)()
在t=TI−t = T_{I-}t=TI−时刻,即180°脉冲施加前,横向和纵向磁化强度分别为:
Mz′(TI−)=Mz0(1−e−TI/T1)Mx′y′(TI−)=Mz0e−TI/T2()
\begin{aligned}
M_{z'}(T_{I-}) &= M_z^0 \left(1 - e^{-T_I/T_1} \right) &\quad\\
M_{x'y'}(T_{I-}) &= M_z^0e^{-T_I/T_2} &\quad
\end{aligned}
\tag{}
Mz′(TI−)Mx′y′(TI−)=Mz0(1−e−TI/T1)=Mz0e−TI/T2()
理想180°射频脉冲给横向和纵向磁化带来的影响如下:
Mz′→180∘x′−Mz′()
M_{z'} \xrightarrow{{180^\circ}_{x'}} -M_{z'} \tag{}
Mz′180∘x′−Mz′()
Mx′y′→180∘x′−Mx′y′()
M_{x'y'} \xrightarrow{{180^\circ}_{x'}} -M_{x'y'}
\tag{}
Mx′y′180∘x′−Mx′y′()
在t=TI+t = T_{I+}t=TI+时刻,即180°脉冲施加后,横向和纵向磁化强度分别为:
Mz′(TI+)=−Mz′(TI−)=−Mz0(1−e−TI/T1)Mx′y′(TI+)=−Mx′y′(TI−)=−Mz0e−TI/T2()
\begin{aligned}
M_{z'}({T_I}+) &= -{M_{z'}({T_I}_-)} = -M_z^0 \left(1 - e^{-T_I/T_1} \right)\\
M_{x'y'}({T_I}+) &= -{M_{x'y'}({T_I}_-)} = -M_z^0e^{-T_I/T_2}
\end{aligned}
\tag{}
Mz′(TI+)Mx′y′(TI+)=−Mz′(TI−)=−Mz0(1−e−TI/T1)=−Mx′y′(TI−)=−Mz0e−TI/T2()
在t>TI+t > T_{I+}t>TI+时刻,横向和纵向磁化强度分别为:
Mz′(t)=Mz0(1−e−(t−TI)/T1)−Mz0(1−e−TI/T1)e−(t−TI)/T1=Mz0(1−e−(t−TI)/T1)−Mz0(e−(t−TI)/T1−e−t/T1)=Mz0(1−2e−(t−TI)/T1)+Mz0e−t/T1=Mz0(1−2e−(t−TI)/T1+e−t/T1)Mx′y′(t)=−Mz0e−t/T2()
\begin{aligned}
M_{z'}(t) &= M_z^0 (1 - e^{-(t-T_I)/T_1}) -M_z^0 \left(1 - e^{-T_I/T_1} \right) e^{-(t-T_I)/T_1} &\quad\\
&=M_z^0 (1 - e^{-(t-T_I)/T_1}) -M_z^0 \left(e^{-(t-T_I)/T_1} - e^{-t/T_1} \right) \\
&=M_z^0 (1 - 2e^{-(t-T_I)/T_1}) + M_z^0 e^{-t/T_1} \\
&=M_z^0 (1 - 2e^{-(t-T_I)/T_1}+e^{-t/T_1}) \\
\\
M_{x'y'}(t) &= -M_z^0e^{-t/T_2} &\quad
\end{aligned}
\tag{}
Mz′(t)Mx′y′(t)=Mz0(1−e−(t−TI)/T1)−Mz0(1−e−TI/T1)e−(t−TI)/T1=Mz0(1−e−(t−TI)/T1)−Mz0(e−(t−TI)/T1−e−t/T1)=Mz0(1−2e−(t−TI)/T1)+Mz0e−t/T1=Mz0(1−2e−(t−TI)/T1+e−t/T1)=−Mz0e−t/T2()
第2次饱和恢复:TR≤t<2TRT_R \leq t< 2T_RTR≤t<2TR
在t=TR−t = {T_R}_-t=TR−时刻,横向和纵向磁化强度分别为:
Mz′(TR−)=Mz0(1−2e−(TR−TI)/T1−e−TR/T1)Mx′y′(TR−)=−Mz0e−TR/T2()
\begin{aligned}
M_{z'}({T_R}_-) &=M_z^0 (1 - 2e^{-(T_R-T_I)/T_1}-e^{-T_R/T_1}) \\
\\
M_{x'y'}({T_R}_-) &= -M_z^0e^{-T_R/T_2} &\quad
\end{aligned}
\tag{}
Mz′(TR−)Mx′y′(TR−)=Mz0(1−2e−(TR−TI)/T1−e−TR/T1)=−Mz0e−TR/T2()
这里将TE=2TIT_E = 2 T_ITE=2TI代入上式,可得在时间t=TR−t = {T_R}_-t=TR−时刻,横向和纵向磁化强度分别为:
Mz′(TR−)=Mz0(1−2e−(TR−TE/2)/T1−e−TR/T1)Mx′y′(TR−)=−Mz0e−TR/T2()
\begin{aligned}
M_{z'}({T_R}_-) &=M_z^0 (1 - 2e^{-(T_R-T_E/2)/T_1}-e^{-T_R/T_1}) \\
\\
M_{x'y'}({T_R}_-) &= -M_z^0e^{-T_R/T_2} &\quad
\end{aligned}
\tag{}
Mz′(TR−)Mx′y′(TR−)=Mz0(1−2e−(TR−TE/2)/T1−e−TR/T1)=−Mz0e−TR/T2()
为了符号方便阅读,这里将TR−→0−{T_R}_- \to 0_-TR−→0−,可得:
Mz′(0−)=Mz0(1−2e−(TR−TE/2)/T1−e−TR/T1)Mx′y′(0−)=−Mz0e−TR/T2()
\begin{aligned}
M_{z'}({0}_-) &=M_z^0 (1 - 2e^{-(T_R-T_E/2)/T_1}-e^{-T_R/T_1}) \\
\\
M_{x'y'}({0}_-) &= -M_z^0e^{-T_R/T_2} &\quad
\end{aligned}
\tag{}
Mz′(0−)Mx′y′(0−)=Mz0(1−2e−(TR−TE/2)/T1−e−TR/T1)=−Mz0e−TR/T2()
由于T2<TE<<TRT_2 < T_E<<T_RT2<TE<<TR,所以上面的式子中横向磁化可以忽略不计。因此
在t=0−t = 0_-t=0−时刻,即施加90°射频脉冲前,纵向和横向磁化强度分别为
Mz′(0−)=Mz0(1−2e−(TR−TE/2)/T1−e−TR/T1)Mx′y′(0−)~0()
\begin{aligned}
M_{z'}({0}_-) &=M_z^0 (1 - 2e^{-(T_R-T_E/2)/T_1}-e^{-T_R/T_1}) \\
\\
M_{x'y'}({0}_-) &~ 0 &\quad
\end{aligned}
\tag{}
Mz′(0−)Mx′y′(0−)=Mz0(1−2e−(TR−TE/2)/T1−e−TR/T1)~0()
在t=0+t = 0_+t=0+时刻,即施加90°射频脉冲后,纵向和横向磁化强度分别为:
Mz′(0+)=0Mx′y′(0+)=Mz0(1−2e−(TR−TE/2)/T1−e−TR/T1)()
\begin{aligned}
M_z'(0_+) &= 0\\
M_{x'y'}(0_+) &= M_z^0 (1 - 2e^{-(T_R-T_E/2)/T_1}-e^{-T_R/T_1})
\end{aligned}\tag{}
Mz′(0+)Mx′y′(0+)=0=Mz0(1−2e−(TR−TE/2)/T1−e−TR/T1)()
根据弛豫方程(式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)
在0≤t<TI0 \leq t < T_I0≤t<TI 时刻,横向和纵向磁化强度分别为:
Mz′(t)=Mz0(1−e−t/T1)(0<t<TI)Mx′y′(t)=Mz0(1−2e−(TR−TE/2)/T1−e−TR/T1)e−t/T2(0<t<TI)()
\begin{aligned}
M_{z'}(t) &=M_z^0 \left(1 - e^{-t/T_1}\right) &\quad (0 < t < T_I) \\
M_{x'y'}(t) &=M_z^0 (1 - 2e^{-(T_R-T_E/2)/T_1}-e^{-T_R/T_1}) e^{-t/T_2} &\quad (0 < t < T_I)
\end{aligned}
\tag{}
Mz′(t)Mx′y′(t)=Mz0(1−e−t/T1)=Mz0(1−2e−(TR−TE/2)/T1−e−TR/T1)e−t/T2(0<t<TI)(0<t<TI)()
在t=TI−t = T_{I-}t=TI−时刻,即180°脉冲施加前,横向和纵向磁化强度分别为:
Mz′(TI−)=Mz0(1−e−TI/T1)Mx′y′(TI−)=Mz0(1−2e−(TR−TE/2)/T1−e−TR/T1)e−TI/T2()
\begin{aligned}
M_{z'}({T_I}_-) &=M_z^0 \left(1 - e^{-{T_I}/T_1}\right) &\quad \\
M_{x'y'}({T_I}_-) &=M_z^0 (1 - 2e^{-(T_R-T_E/2)/T_1}-e^{-T_R/T_1}) e^{-{T_I}/T_2} &\quad
\end{aligned}
\tag{}
Mz′(TI−)Mx′y′(TI−)=Mz0(1−e−TI/T1)=Mz0(1−2e−(TR−TE/2)/T1−e−TR/T1)e−TI/T2()
理想180°射频脉冲给横向和纵向磁化带来的影响如下:
Mz′→180∘x′−Mz′()
M_{z'} \xrightarrow{{180^\circ}_{x'}} -M_{z'} \tag{}
Mz′180∘x′−Mz′()
Mx′y′→180∘x′−Mx′y′()
M_{x'y'} \xrightarrow{{180^\circ}_{x'}} -M_{x'y'}
\tag{}
Mx′y′180∘x′−Mx′y′()
在t=TI+t = T_{I+}t=TI+时刻,即180°脉冲施加后,横向和纵向磁化强度分别为:
Mz′(TI+)=−Mz′(TI−)=−Mz0(1−e−TI/T1)Mx′y′(TI+)=−Mx′y′(TI−)=−Mz0(1−2e−(TR−TE/2)/T1−e−TR/T1)e−TI/T2()
\begin{aligned}
M_{z'}({T_I}+) &= -{M_{z'}({T_I}_-)} = -M_z^0 \left(1 - e^{-{T_I}/T_1}\right)\\
M_{x'y'}({T_I}+) &= -{M_{x'y'}({T_I}_-)} = -M_z^0 (1 - 2e^{-(T_R-T_E/2)/T_1}-e^{-T_R/T_1}) e^{-{T_I}/T_2}
\end{aligned}
\tag{}
Mz′(TI+)Mx′y′(TI+)=−Mz′(TI−)=−Mz0(1−e−TI/T1)=−Mx′y′(TI−)=−Mz0(1−2e−(TR−TE/2)/T1−e−TR/T1)e−TI/T2()
在t>TI+t > T_{I+}t>TI+时刻,横向和纵向磁化强度分别为:
Mz′(t)=Mz0(1−e−(t−TI)/T1)−Mz0(1−e−TI/T1)e−(t−TI)/T1=Mz0(1−e−(t−TI)/T1)−Mz0(e−(t−TI)/T1−e−t/T1)=Mz0(1−2e−(t−TI)/T1)+Mz0e−t/T1=Mz0(1−2e−(t−TI)/T1+e−t/T1)Mx′y′(t)=−Mz0(1−2e−(TR−TE/2)/T1−e−TR/T1)e−t/T2()
\begin{aligned}
M_{z'}(t) &= M_z^0 (1 - e^{-(t-T_I)/T_1}) -M_z^0 \left(1 - e^{-T_I/T_1} \right) e^{-(t-T_I)/T_1} &\quad\\
&=M_z^0 (1 - e^{-(t-T_I)/T_1}) -M_z^0 \left(e^{-(t-T_I)/T_1} - e^{-t/T_1} \right) \\
&=M_z^0 (1 - 2e^{-(t-T_I)/T_1}) + M_z^0 e^{-t/T_1} \\
&=M_z^0 (1 - 2e^{-(t-T_I)/T_1}+e^{-t/T_1}) \\
\\
M_{x'y'}(t) &= -M_z^0 (1 - 2e^{-(T_R-T_E/2)/T_1}-e^{-T_R/T_1})e^{-t/T_2} &\quad
\end{aligned}
\tag{}
Mz′(t)Mx′y′(t)=Mz0(1−e−(t−TI)/T1)−Mz0(1−e−TI/T1)e−(t−TI)/T1=Mz0(1−e−(t−TI)/T1)−Mz0(e−(t−TI)/T1−e−t/T1)=Mz0(1−2e−(t−TI)/T1)+Mz0e−t/T1=Mz0(1−2e−(t−TI)/T1+e−t/T1)=−Mz0(1−2e−(TR−TE/2)/T1−e−TR/T1)e−t/T2()
由上式可得:
-
第一次循环后系统就已经进入了稳态,即后续的循环的起始时刻t=0−t=0_-t=0−时的纵向磁化强度都是:
Mz′(0−)=Mz0(1−2e−(TR−TE/2)/T1−e−TR/T1)(7.24) \begin{aligned} M_{z'}({0}_-) &=M_z^0 (1 - 2e^{-(T_R-T_E/2)/T_1}-e^{-T_R/T_1}) \\ \end{aligned} \tag{7.24} Mz′(0−)=Mz0(1−2e−(TR−TE/2)/T1−e−TR/T1)(7.24) -
自旋回波信号的幅度可表示为:
AE=Mz0(1−2e−(TR−TE/2)/T1+e−TR/T1)e−TE/T2(7.25) A_E = M_z^0 \left(1 - 2e^{-(T_R - T_E/2)/T_1} + e^{-T_R/T_1}\right)e^{-T_E/T_2} \tag{7.25} AE=Mz0(1−2e−(TR−TE/2)/T1+e−TR/T1)e−TE/T2(7.25)
在实际成像中,通常满足TE≪TRT_E \ll T_RTE≪TR的条件(即回波时间远小于重复时间)。在此条件下,可对公式(7.25)进行合理简化:当TE≪TRT_E \ll T_RTE≪TR时,有TR−TE/2≈TRT_R - T_E/2 \approx T_RTR−TE/2≈TR,所以:
e−(TR−TE/2)/T1≈e−TR/T1
e^{-(T_R - T_E/2)/T_1} \approx e^{-T_R/T_1}
e−(TR−TE/2)/T1≈e−TR/T1
代入公式(7.25)得:
AE=Mz0(1−2e−TR/T1+e−TR/T1)e−TE/T2=Mz0(1−e−TR/T1)e−TE/T2(7.26)
A_E = M_z^0 \left(1 - 2e^{-T_R/T_1} + e^{-T_R/T_1}\right)e^{-T_E/T_2} = M_z^0 \left(1 - e^{-T_R/T_1}\right)e^{-T_E/T_2} \tag{7.26}
AE=Mz0(1−2e−TR/T1+e−TR/T1)e−TE/T2=Mz0(1−e−TR/T1)e−TE/T2(7.26)
选择性对比度强调
公式(7.26)表明,饱和恢复自旋回波序列的图像强度同时具有三种加权特性:
T1T_1T1加权对比度:采用短TET_ETE和适当的TRT_RTR
- 短TET_ETE使e−TE/T2→1e^{-T_E/T_2} \rightarrow 1e−TE/T2→1,消除T2T_2T2加权效应
- 适当TRT_RTR确保(1−e−TR/T1)(1 - e^{-T_R/T_1})(1−e−TR/T1)产生明显的T1T_1T1对比
T2T_2T2加权对比度:采用长TET_ETE和长TRT_RTR
- 长TRT_RTR使(1−e−TR/T1)→1(1 - e^{-T_R/T_1}) \rightarrow 1(1−e−TR/T1)→1,消除T1T_1T1加权效应
- 长TET_ETE使e−TE/T2e^{-T_E/T_2}e−TE/T2产生明显的T2T_2T2对比
自旋密度加权对比度:采用短TET_ETE和长TRT_RTR
- 短TET_ETE消除T2T_2T2加权
- 长TRT_RTR消除T1T_1T1加权
- 剩余信号正比于自旋密度
下表总结了不同序列参数值对图像对比度的影响。

反转恢复自旋回波序列

上展示了反转恢复自旋回波序列的时序图。该序列首先施加180°反转脉冲,随后在反转时间TIT_ITI后施加90°激发脉冲,最后通过180°重聚脉冲产生自旋回波。
对于反转恢复自旋回波序列,自旋回波信号的幅度可推导为:
AE=Mz0(1−2e−TI/T1+2e−(TR−TE/2)/T1−e−TR/T1)e−TE/T2(7.27)
A_E = M_z^0 \left(1 - 2e^{-T_I/T_1} + 2e^{-(T_R - T_E/2)/T_1} - e^{-T_R/T_1}\right)e^{-T_E/T_2} \tag{7.27}
AE=Mz0(1−2e−TI/T1+2e−(TR−TE/2)/T1−e−TR/T1)e−TE/T2(7.27)
同样基于TE≪TRT_E \ll T_RTE≪TR的实用假设,公式(7.27)可简化为:
AE=Mz0(1−2e−TI/T1+e−TR/T1)e−TE/T2(7.28)
A_E = M_z^0 \left(1 - 2e^{-T_I/T_1} + e^{-T_R/T_1}\right)e^{-T_E/T_2} \tag{7.28}
AE=Mz0(1−2e−TI/T1+e−TR/T1)e−TE/T2(7.28)
该表达式表明,反转恢复自旋回波序列能够产生与饱和恢复自旋回波序列相似的T1T_1T1、T2T_2T2和自旋密度加权对比度,但具有更强的T1T_1T1对比度调控能力。
结论
基本自旋回波成像序列通过将自旋回波技术引入饱和恢复和反转恢复序列,实现了对图像对比度的精细调控。公式(7.24)至(7.28)完整描述了这些序列的信号特性,为临床MRI扫描参数优化提供了理论基础。
通过合理选择TRT_RTR、TET_ETE和TIT_ITI等参数,可以获得满足特定诊断需求的图像对比度,这是自旋回波序列在临床磁共振成像中得以广泛应用的重要原因。
