Медицина · Фізика · Обробка сигналів
📅 Березень 2026 ⏱ ≈ 14 хв читання 🎯 Середній рівень · Останнє оновлення: 28 травня 2026 р.

Як працює МРТ — ядерний магнітний резонанс і медична візуалізація

Щороку у світі виконують понад 100 мільйонів МРТ-сканувань — отримуючи детальні 3-D зображення м'яких тканин без жодного рентгенівського фотона. Апарат використовує цікаву квантово-механічну властивість атомних ядер: їхній спін. Ця стаття простежує фізику від магнітного моменту протона аж до зображення на екрані рентгенолога.

Автор: Команда MySimulator · Редакційна перевірка: Редакція MySimulator

1. Спін протона та магнітні моменти

Кожен протон (ядро водню) має внутрішній кутовий момент, який називають спіном. Хоча спін — це квантова властивість без класичного аналога, макроскопічно він поводиться як крихітний стрижневий магніт: кожен протон несе магнітний момент μ, вирівняний уздовж осі свого спіну.

μ = γ · ℏ · I γ (гіромагнітне відношення для ¹H) = 2.675 × 10⁸ rad s⁻¹ T⁻¹ ℏ = зведена стала Планка = 1.055 × 10⁻³⁴ J·s I = спінове квантове число = ½ для ¹H

Тіло людини на ~60% складається з води (H₂O), що робить водень найпоширенішим ЯМР-активним ядром. Тому клінічна МРТ передусім візуалізує щільність протонів і релаксаційні властивості тканин.

2. Вирівнювання в сильному статичному полі (B₀)

Поза магнітним полем спіни протонів спрямовані у випадкових напрямках і взаємно компенсуються. Усередині тунелю МРТ-сканера надпровідний магніт створює потужне статичне поле B₀ (зазвичай 1.5 T або 3 T — до 50 000 разів сильніше за земне поле).

Квантова механіка дозволяє лише два енергетичні стани для ядра зі спіном ½ у полі B₀:

Паралельний (низька енергія, спін угору): E₊ = −½ γ ℏ B₀ Антипаралельний (висока енергія, спін униз): E₋ = +½ γ ℏ B₀ Різниця енергій: ΔE = γ ℏ B₀ За температури тіла існує невеликий больцманівський надлишок ~7 спінів на мільйон у нижчому за енергією (паралельному) стані. Цей крихітний дисбаланс створює сумарну макроскопічну намагніченість M₀, спрямовану вздовж B₀. M₀ ∝ N · γ² ℏ² B₀ / (4 k_B T) (закон Кюрі) N = щільність числа протонів, k_B = стала Больцмана, T = температура

Сумарна намагніченість M₀ — це те, що зрештою вимірює сканер. Вона пропорційна напруженості поля B₀ — одна з причин, чому сканер 3 T дає чіткіші зображення, ніж сканер 1.5 T.

3. РЧ-імпульс і резонанс Лармора

Щоб отримати реєстрований сигнал, сканер відхиляє сумарну намагніченість M₀ від осі B₀ за допомогою короткого радіочастотного (РЧ) імпульсу. Ключ — це резонанс: частота імпульсу має точно збігатися з частотою Лармора ω₀.

ω₀ = γ · B₀ (рівняння Лармора) Для ¹H за 1.5 T: ω₀ = 2.675×10⁸ × 1.5 ≈ 64 MHz Для ¹H за 3.0 T: ω₀ ≈ 128 MHz Для ¹H за 7.0 T: ω₀ ≈ 298 MHz (FM радіодіапазон)

Імпульс 90° повністю відхиляє M₀ у поперечну (x-y) площину. Потім намагніченість прецесує навколо B₀ з частотою ω₀, проходячи повз приймальну котушку та індукуючи синусоїдальну напругу — сигнал вільного спадання індукції (FID).

Рівняння Блоха описують вектор намагніченості M(t) за комбінованої дії статичного поля B₀, РЧ-поля B₁ та релаксації:

dMx/dt = γ(My·Bz − Mz·By) − Mx/T2
dMy/dt = γ(Mz·Bx − Mx·Bz) − My/T2
dMz/dt = γ(Mx·By − My·Bx) − (Mz − M₀)/T1

4. T1- і T2-релаксація — звідки береться контраст

Після вимкнення РЧ-імпульсу намагніченість повертається до рівноваги через два незалежні процеси. Ці характерні сталі часу є основним джерелом контрасту тканин у МРТ.

Поздовжня (T1) релаксація — спін-ґраткова

z-компонента намагніченості (Mz) експоненційно відновлюється до M₀, оскільки спіни віддають енергію навколишній молекулярній ґратці:

Mz(t) = M₀ · [1 − e^(−t/T1)] Типові значення T1 за 1.5 T: Біла речовина: ~780 ms Сіра речовина: ~920 ms Жир: ~260 ms Вода (ліквор): ~2400 ms

T1 залежить від того, наскільки ефективно молекули обертаються на частоті Лармора. Ліпіди обертаються саме з потрібною швидкістю, що дає жиру короткий T1 (швидке відновлення, виглядає яскравим на T1-зважених зображеннях).

Поперечна (T2) релаксація — спін-спінова

Поперечна намагніченість (Mxy) спадає, оскільки окремі спіни розфазовуються через крихітні локальні варіації поля від сусідніх ядер:

Mxy(t) = Mxy(0) · e^(−t/T2) Типові значення T2 за 1.5 T: Біла речовина: ~90 ms Сіра речовина: ~100 ms Жир: ~80 ms Вода (ліквор): ~1800 ms Пухлини: часто підвищений T2 (більше вільної води)

T2 завжди ≤ T1. Вільна вода має довгі T1 і T2 (повільне обертання, поганий перенос енергії); зв'язана вода в тканині має короткий T2. Пухлини та набряк часто виглядають гіперінтенсивними на T2-зважених зображеннях, оскільки патологічні процеси збільшують вміст води в тканині.

T2* (T2-зірка) — це швидший спад, що враховує неоднорідності поля (ΔB₀). Спін-ехо послідовність перефокусовує ефекти T2* за допомогою імпульсу 180°, відновлюючи справжній сигнал T2. Градієнт-ехо послідовності зберігають чутливість до T2* — це використовується у BOLD-візуалізації fMRI.

Імпульсні послідовності керують контрастом

Обираючи TR (час повторення) та TE (час ехо), рентгенолог придушує або підкреслює внески T1 і T2:

5. Градієнтні котушки та просторове кодування

Після РЧ-імпульсу всі протони дають сигнал на одній і тій самій частоті Лармора — сканер ще не може визначити, звідки надійшов сигнал. Три набори градієнтних котушок (Gx, Gy, Gz) додають невеликі, лінійні варіації до B₀, що кодують просторове положення у частоту та фазу.

Вибір зрізу

Градієнт Gz застосовується під час РЧ-імпульсу, змушуючи B₀ змінюватися вздовж z. Збуджується лише той зріз, де локальне поле збігається з частотою імпульсу:

ω(z) = γ · [B₀ + Gz · z] Положення зрізу: z₀ = (ω_RF − γ·B₀) / (γ·Gz) Товщина зрізу: Δz = BW_pulse / (γ·Gz) BW_pulse = РЧ-смуга пропускання; тонший зріз → вужча смуга або сильніший градієнт

Частотне кодування

Під час зчитування сигналу застосовується градієнт Gx. Протони в різних положеннях x прецесують з різними частотами; отриманий сигнал є суперпозицією синусоїд:

ω(x) = γ · [B₀ + Gx · x] Signal(t) = ∫ ρ(x) · e^(iγGx·x·t) dx → перетворення Фур'є щільності протонів ρ(x)!

Фазове кодування

Градієнт Gy застосовується на короткий період перед зчитуванням. Це надає фазовий зсув, пропорційний положенню y — кодуючи другий просторовий вимір. Уся послідовність повторюється з різними амплітудами Gy, щоб заповнити 2-D матрицю k-простору.

6. k-простір і 2-D перетворення Фур'є

Необроблені дані, зібрані сканером, зберігаються у 2-D матриці, яку називають k-простором. Кожна точка [kx, ky] представляє певну просторову частоту зображення:

kx(t) = γ/(2π) · Gx · t (напрямок частотного кодування) ky(n) = γ/(2π) · Gy(n) · τ_pe (напрямок фазового кодування, n = 1…N) S(kx, ky) = ∫∫ ρ(x,y) · e^(−i2π(kx·x + ky·y)) dx dy → S(kx, ky) є точно 2-D перетворенням Фур'є щільності протонів ρ(x,y)! Реконструкція зображення: ρ(x,y) = 2D-IFFT{ S(kx, ky) }

Центр k-простору (низькі kx, ky) містить низькі просторові частоти — загальну яскравість і контраст. Периферія містить високі просторові частоти — краї та дрібні деталі. Саме тому часткові збирання k-простору (вибірка лише центру) швидші, але розмитіші. 

// Спрощена реконструкція k-простору → зображення
function reconstructMRI(kSpaceData) {
  // kSpaceData: 2-D комплексний масив [Nky][Nkx]
  const N = kSpaceData.length;

  // Крок 1: 2-D обернене FFT (рядок за рядком, потім стовпець за стовпцем)
  const temp = kSpaceData.map(row => ifft1D(row));       // IFFT кожного рядка
  const image = transposeAndIFFT(temp);                   // IFFT кожного стовпця

  // Крок 2: зображення модуля (відкидаємо фазу)
  return image.map(row =>
    row.map(c => Math.sqrt(c.re**2 + c.im**2))
  );
}

// Властивості k-простору, що використовуються у швидких послідовностях:
// Часткове Фур'є: збираємо 60% k-простору, решту заповнюємо нулями (homodyne)  
// Паралельна візуалізація (GRAPPA/SENSE): недовибірка з багатоканальною котушкою
// Стиснене вимірювання: некогерентна недовибірка + розріджена реконструкція
EPI (ехопланарна візуалізація) заповнює весь k-простір після одного РЧ-імпульсу, швидко осцилюючи градієнти Gx і Gy зигзагоподібно. Час збирання: ~50 ms на зріз — досить швидко для fMRI (активність мозку змінюється протягом секунд) і кардіовізуалізації.

7. Механізми контрасту та поширені послідовності

Контрастні агенти на основі гадолінію

Гадоліній (Gd³⁺) є сильно парамагнітним — він різко скорочує T1 сусідніх протонів. Введений внутрішньовенно, Gd накопичується там, де порушений гематоенцефалічний бар'єр (активні пухлини, запалення), створюючи яскраві ділянки на T1-зважених зображеннях.

fMRI — залежна від рівня кисню в крові (BOLD)

Дезоксигемоглобін є парамагнітним і скорочує T2*; оксигемоглобін діамагнітний і не скорочує. Активні ділянки мозку мають вищу оксигенацію крові → довший T2* → яскравіший сигнал на градієнт-ехо послідовностях. Цей BOLD-контраст є основою функціональної МРТ — картування активності мозку без радіоактивних індикаторів.

Дифузійно-тензорна візуалізація (DTI)

Застосовуючи сильні градієнтні імпульси в кількох напрямках, DTI вимірює, як молекули води дифундують крізь тканину. У пучках аксонів білої речовини дифузія анізотропна — швидша вздовж аксона, ніж поперек нього. DTI виявляє нейронні волоконні тракти й виявляє раннє аксональне пошкодження при інсульті та розсіяному склерозі.

Видимий коефіцієнт дифузії: ADC = −(1/b) · ln(S/S₀) b-значення: b = γ² G² δ² (Δ − δ/3) [s/mm²] G = напруженість градієнта, δ = тривалість, Δ = інтервал Високе b (1000 s/mm²) → чутливе до обмеженої дифузії (гострий інсульт: яскравий DWI)

МР-спектроскопія (MRS)

Без просторових градієнтів ЯМР-спектр розкриває хімічну сигнатуру метаболітів. NAA (N-ацетиласпартат) є нейрональним маркером; його зниження сигналізує про втрату нейронів. Підвищення холіну вказує на оновлення мембран (пухлини). Лактат з'являється при анаеробному метаболізмі (ішемія).

8. Безпека, напруженості поля та межі МРТ

Чому МРТ безпечна (для більшості пацієнтів)

На відміну від КТ чи ПЕТ, МРТ не використовує іонізуюче випромінювання. Внесена РЧ-енергія (SAR — питома швидкість поглинання) обмежена 4 W/kg настановами з безпеки. Основні небезпеки:

Тенденції напруженості поля

Стимуляція периферичних нервів, спричинена градієнтами (PNS) обмежує, наскільки швидко можна перемикати градієнтні котушки. Швидке перемикання (високе dB/dt) індукує струми в тілі — відчуття поколювання, яке обмежує продуктивність EPI та дифузійних послідовностей у системах із сильним полем.
🔬 Дослідіть симуляції з медицини та біофізики →