Кейс: магнитно-резонансная томография — объясните физику ядерного магнитного резонанса, роль магнитного поля и радиочастотных импульсов, и почему различные ткани дают различный контраст в МРТ-изображениях
Кратко о физике ядерного магнитного резонанса (ЯМР) и о том, почему в МРТ разные ткани дают разный контраст. 1) Основные явления - Протоны (в основном ядра водорода) имеют спин и связанный магнитный момент. В однородном статическом магнитном поле B0B_0B0 их суммарная магнитизация ориентируется вдоль поля и равна M0M_0M0. - Индивидуальные спины прецессируют с угловой частотой Лармора ω0=γB0\displaystyle \omega_0=\gamma B_0ω0=γB0, где γ\gammaγ — гиромагнитное отношение протона. 2) Возбуждение радиочастотой - Радиочастотное (РЧ) поле B1B_1B1, перпендикулярное B0B_0B0 и с частотой ω0\omega_0ω0, резонансно взаимодействует со спинами и «переводит» магнитизацию в поперечную плоскость (флип-угол) α=γB1τ\displaystyle \alpha=\gamma B_1\tauα=γB1τ, где τ\tauτ — длительность импульса. - После импульса возникает сигнал в катушке при детектировании поперечной составляющей магнитизации. 3) Релаксация (возвращение к равновесию) - Два ключевых времени релаксации: - продольное (спин-решётка) T1T_1T1: восстановление продольной компоненты Mz(t)=M0(1−e−t/T1)\displaystyle M_z(t)=M_0\big(1-e^{-t/T_1}\big)Mz(t)=M0(1−e−t/T1); - поперечное (спин-спин) T2T_2T2: затухание поперечной компоненте Mxy(t)=Mxy(0)e−t/T2\displaystyle M_{xy}(t)=M_{xy}(0)e^{-t/T_2}Mxy(t)=Mxy(0)e−t/T2. - Также существует T2∗T_2^*T2∗ — более короткое время, включающее неоднородности поля: 1T2∗=1T2+1Tinhom\displaystyle \frac{1}{T_2^*}=\frac{1}{T_2}+\frac{1}{T_{\text{inhom}}}T2∗1=T21+Tinhom1. 4) Уравнения движения (суммарно) - Динамика магнитизации описывается уравнениями Блоха: dMdt=γM×B−(Mx/T2My/T2(Mz−M0)/T1).\displaystyle \frac{d\mathbf{M}}{dt}=\gamma\mathbf{M}\times\mathbf{B}-\begin{pmatrix}M_x/T_2\\[4pt]M_y/T_2\\[4pt](M_z-M_0)/T_1\end{pmatrix}.dtdM=γM×B−Mx/T2My/T2(Mz−M0)/T1. 5) Пространственное кодирование (получение изображения) - Наложением градиентов поля GGG делается частота Лармора зависимой от координаты, например ω(x)=γ(B0+Gxx)\displaystyle \omega(x)=\gamma\big(B_0+G_x x\big)ω(x)=γ(B0+Gxx). - Используют селекцию среза, частотную и фазовую кодировку; затем сигнал преобразуют по Фурье в пространственное изображение. 6) Причины контрастов между тканями - Контраст определяется главным образом: - плотностью протонов (PD): количество возбужденных ядер в объёме; - различиями T1T_1T1 и T2T_2T2 у тканей (зависит от молекулярной подвижности, вязкости, взаимодействий с макромолекулами и т.д.); - химическим сдвигом и магнитной восприимчивостью (особенно в градиент-эхо последовательностях). - Практически: короткие TRT_RTR и короткие TET_ETE дают T1-взвешивание (жировая ткань обычно ярче), длинные TET_ETE дают T2-взвешивание (жидкости, например ликвор, ярче), а короткий TET_ETE и длинный TRT_RTR дают PD-взвешивание. Обобщённо параметры последовательности TRT_RTR (repetition time) и TET_ETE (echo time) управляют вкладом T1T_1T1, T2T_2T2 и PD в итоговый контраст. 7) Технические приёмы повышения/изменения контраста - Спин-эхо (SE) с 180°-пульсом восстанавливает дефазирование от неоднородностей -> измеряется близко к T2T_2T2. - Градиент-эхо (GRE) более чувствителен к T2∗T_2^*T2∗ и восприимчивости. - Дополнительные методы (инверсия-восстановление, контрастные препараты) позволяют выделять специфические ткани или процессы. Коротко: статическое поле B0B_0B0 упорядочивает спины и задаёт частоту резонанса; РЧ-импульсы резонансно переворачивают магнитизацию; возвращение к состоянию равновесия описывается временами T1T_1T1 и T2T_2T2; различия в плотности протонов и в этих временах у разных тканей и при разных последовательностях дают различный МР-контраст.
1) Основные явления
- Протоны (в основном ядра водорода) имеют спин и связанный магнитный момент. В однородном статическом магнитном поле B0B_0B0 их суммарная магнитизация ориентируется вдоль поля и равна M0M_0M0 .
- Индивидуальные спины прецессируют с угловой частотой Лармора
ω0=γB0\displaystyle \omega_0=\gamma B_0ω0 =γB0 ,
где γ\gammaγ — гиромагнитное отношение протона.
2) Возбуждение радиочастотой
- Радиочастотное (РЧ) поле B1B_1B1 , перпендикулярное B0B_0B0 и с частотой ω0\omega_0ω0 , резонансно взаимодействует со спинами и «переводит» магнитизацию в поперечную плоскость (флип-угол)
α=γB1τ\displaystyle \alpha=\gamma B_1\tauα=γB1 τ,
где τ\tauτ — длительность импульса.
- После импульса возникает сигнал в катушке при детектировании поперечной составляющей магнитизации.
3) Релаксация (возвращение к равновесию)
- Два ключевых времени релаксации:
- продольное (спин-решётка) T1T_1T1 : восстановление продольной компоненты
Mz(t)=M0(1−e−t/T1)\displaystyle M_z(t)=M_0\big(1-e^{-t/T_1}\big)Mz (t)=M0 (1−e−t/T1 );
- поперечное (спин-спин) T2T_2T2 : затухание поперечной компоненте
Mxy(t)=Mxy(0)e−t/T2\displaystyle M_{xy}(t)=M_{xy}(0)e^{-t/T_2}Mxy (t)=Mxy (0)e−t/T2 .
- Также существует T2∗T_2^*T2∗ — более короткое время, включающее неоднородности поля:
1T2∗=1T2+1Tinhom\displaystyle \frac{1}{T_2^*}=\frac{1}{T_2}+\frac{1}{T_{\text{inhom}}}T2∗ 1 =T2 1 +Tinhom 1 .
4) Уравнения движения (суммарно)
- Динамика магнитизации описывается уравнениями Блоха:
dMdt=γM×B−(Mx/T2My/T2(Mz−M0)/T1).\displaystyle \frac{d\mathbf{M}}{dt}=\gamma\mathbf{M}\times\mathbf{B}-\begin{pmatrix}M_x/T_2\\[4pt]M_y/T_2\\[4pt](M_z-M_0)/T_1\end{pmatrix}.dtdM =γM×B− Mx /T2 My /T2 (Mz −M0 )/T1 .
5) Пространственное кодирование (получение изображения)
- Наложением градиентов поля GGG делается частота Лармора зависимой от координаты, например
ω(x)=γ(B0+Gxx)\displaystyle \omega(x)=\gamma\big(B_0+G_x x\big)ω(x)=γ(B0 +Gx x).
- Используют селекцию среза, частотную и фазовую кодировку; затем сигнал преобразуют по Фурье в пространственное изображение.
6) Причины контрастов между тканями
- Контраст определяется главным образом:
- плотностью протонов (PD): количество возбужденных ядер в объёме;
- различиями T1T_1T1 и T2T_2T2 у тканей (зависит от молекулярной подвижности, вязкости, взаимодействий с макромолекулами и т.д.);
- химическим сдвигом и магнитной восприимчивостью (особенно в градиент-эхо последовательностях).
- Практически: короткие TRT_RTR и короткие TET_ETE дают T1-взвешивание (жировая ткань обычно ярче), длинные TET_ETE дают T2-взвешивание (жидкости, например ликвор, ярче), а короткий TET_ETE и длинный TRT_RTR дают PD-взвешивание.
Обобщённо параметры последовательности TRT_RTR (repetition time) и TET_ETE (echo time) управляют вкладом T1T_1T1 , T2T_2T2 и PD в итоговый контраст.
7) Технические приёмы повышения/изменения контраста
- Спин-эхо (SE) с 180°-пульсом восстанавливает дефазирование от неоднородностей -> измеряется близко к T2T_2T2 .
- Градиент-эхо (GRE) более чувствителен к T2∗T_2^*T2∗ и восприимчивости.
- Дополнительные методы (инверсия-восстановление, контрастные препараты) позволяют выделять специфические ткани или процессы.
Коротко: статическое поле B0B_0B0 упорядочивает спины и задаёт частоту резонанса; РЧ-импульсы резонансно переворачивают магнитизацию; возвращение к состоянию равновесия описывается временами T1T_1T1 и T2T_2T2 ; различия в плотности протонов и в этих временах у разных тканей и при разных последовательностях дают различный МР-контраст.