🧲 Fizyka MRI / NMR — równania Blocha
Interaktywny symulator MRI/NMR oparty na równaniach Blocha. Obserwuj, jak impulsy RF odchylają wektor magnetyzacji, relaksację T1 i T2 oraz sygnał swobodnej precesji. Presety tkanek: woda, tłuszcz, mięsień, istota szara/biała i płyn mózgowo-rdzeniowy.
🧲 Co pokazuje symulacja
Ten symulator modeluje jądrowy rezonans magnetyczny (NMR) — fizyczną podstawę MRI — za pomocą równań Blocha. Wektor magnetyzacji netto M = (Mx, My, Mz) zaczyna w linii ze statycznym polem B₀ wzdłuż osi z. Impuls RF przechyla M do płaszczyzny poprzecznej, gdzie precesuje z częstotliwością Larmora i zanika wskutek dwóch niezależnych procesów relaksacji: T1 (spin-sieć, powrót podłużny) i T2 (spin-spin, zanik poprzeczny). Zanikający sygnał poprzeczny to sygnał swobodnej precesji (FID).
Jak korzystać
- Kliknij Impuls RF 90°, aby przechylić M do płaszczyzny poprzecznej — obserwuj pojawiający się sygnał FID.
- Kliknij Inwersja 180°, aby uzyskać pełne odzyskiwanie inwersji — Mz przechodzi do −1, a następnie wraca poprzez T1.
- Wybierz presety tkanek, aby wczytać realistyczne wartości T1/T2 dla wody, tłuszczu, mięśni, istoty szarej/białej i płynu mózgowo-rdzeniowego przy 1,5 T.
- Dostosuj B₀, aby zobaczyć, jak natężenie pola zmienia częstotliwość Larmora (γ·B₀, gdzie γ = 42,577 MHz/T dla wodoru).
Czy wiesz, że…
MRI zostało niezależnie opracowane przez Paula Lauterbura i Petera Mansfielda w latach 70. XX wieku; obaj otrzymali w 2003 roku Nagrodę Nobla w dziedzinie medycyny. W przeciwieństwie do CT czy rentgena, MRI nie wykorzystuje promieniowania jonizującego — opiera się wyłącznie na impulsach radiowych i polach magnetycznych. Kliniczny magnes MRI 1,5 T jest około 30 000 razy silniejszy niż ziemskie pole magnetyczne. Różne tkanki mają charakterystyczne czasy T1 i T2, co pozwala MRI tworzyć kontrast między istotą szarą, białą, płynem mózgowo-rdzeniowym i patologią — bez ani jednej kropli środka kontrastowego w wielu sekwencjach.
O fizyce MRI / NMR i równaniach Blocha
Ta symulacja modeluje jądrowy rezonans magnetyczny, fizykę leżącą u podstaw MRI, całkując równania Blocha dla pojedynczego wektora magnetyzacji netto M = (Mx, My, Mz). W spoczynku M leży wzdłuż statycznego pola B₀ na osi z. Każdy krok animacji stosuje relaksację wykładniczą — składowe poprzeczne zanikają jako exp(−dt/T2), podczas gdy składowa podłużna wraca ku M₀=1 jako Mz = 1 + (Mz−1)·exp(−dt/T1) — oraz precesję wizualną.
Przyciski impulsu RF obracają M wokół osi x: impuls 90° przechyla je w pełni do płaszczyzny poprzecznej, impuls 180° je odwraca. Suwaki T1, T2 i B₀ ustawiają czasy relaksacji i natężenie pola, a presety tkanek wczytują realistyczne wartości przy 1,5 T dla wody, tłuszczu, mięśni, istoty szarej i białej oraz płynu mózgowo-rdzeniowego. Zanikający sygnał poprzeczny to sygnał swobodnej precesji, surowy sygnał rejestrowany przez każdy kliniczny skaner MRI.
Najczęściej zadawane pytania
Czym są równania Blocha?
Równania Blocha to zestaw równań różniczkowych opisujących, jak wektor magnetyzacji jądrowej netto ewoluuje w polu magnetycznym. Łączą precesję z dwoma członami relaksacji: zanikiem poprzecznym rządzonym przez T2 i powrotem podłużnym rządzonym przez T1. Ten symulator całkuje je numerycznie w każdej klatce, aby przesuwać wektor magnetyzacji M.
Jaka jest różnica między relaksacją T1 i T2?
T1, czyli relaksacja spin-sieć, to powrót składowej podłużnej Mz do stanu równowagi wraz z uwalnianiem energii do otoczenia. T2, czyli relaksacja spin-spin, to zanik składowej poprzecznej Mxy w miarę utraty spójności fazowej poszczególnych spinów. T1 jest zawsze dłuższe lub równe T2, dlatego suwaki obejmują tu zakres 100 do 5000 ms dla T1 i 10 do 2000 ms dla T2.
Co robi impuls RF 90 stopni?
Przycisk Impuls RF 90° obraca magnetyzację o 90 stopni wokół osi x, przechylając ją z osi podłużnej z do płaszczyzny poprzecznej. Maksymalizuje to wykrywalny sygnał poprzeczny Mxy, który następnie precesuje i zanika, tworząc sygnał swobodnej precesji pokazany na zielono po prawej stronie płótna.
Czym jest częstotliwość Larmora pokazana w statystykach?
Częstotliwość Larmora to tempo, z jakim spiny jądrowe precesują wokół pola statycznego, równe iloczynowi współczynnika żyromagnetycznego i natężenia pola. Symulator używa γ = 42,577 MHz/T dla wodoru, więc przy 1,5 T wyświetlana wartość wynosi około 63,9 MHz. Zwiększenie suwaka B₀ proporcjonalnie zwiększa tę częstotliwość.
Czym jest sygnał swobodnej precesji (FID)?
Sygnał swobodnej precesji to oscylujące, zanikające napięcie indukowane w cewce odbiorczej przez precesującą magnetyzację poprzeczną po impulsie RF. W symulatorze jest wykreślany jako moduł Mxy modulowany cosinusem, zanikający wraz z relaksacją T2 zmniejszającą spójność poprzeczną. To podstawowy surowy sygnał, z którego rekonstruowane są obrazy MRI.
Co robi impuls inwersji 180 stopni?
Przycisk Inwersja 180° obraca magnetyzację o pełne 180 stopni wokół osi x. Zaczynając od równowagi, sprowadza to Mz do −1, po czym wraca ono ku +1 wzdłuż krzywej T1. Odzyskiwanie inwersji jest szeroko stosowane w rzeczywistym MRI do wygaszania konkretnych tkanek, na przykład tłuszczu lub płynu mózgowo-rdzeniowego, poprzez taktowanie odczytu wtedy, gdy ich sygnał przechodzi przez zero.
Dlaczego presety tkanek mają różne wartości T1 i T2?
Każda tkanka ma charakterystyczne środowisko molekularne, więc woda i płyn mózgowo-rdzeniowy relaksują się wolno (długie T1 i T2), podczas gdy tłuszcz i mięśnie relaksują się szybciej. Presety wczytują zmierzone wartości przy 1,5 T — na przykład istota szara mózgu ma około T1 1000 ms i T2 100 ms, a płyn mózgowo-rdzeniowy około T1 4300 ms i T2 2200 ms. Te kontrasty są dokładnie tym, co tworzy różnice obrazu między tkankami.
Czy ta symulacja jest fizycznie dokładna?
Fizyka relaksacji jest wierna: wykorzystuje dokładne rozwiązania wykładnicze członów relaksacji Blocha z poprawnym współczynnikiem żyromagnetycznym i realistycznymi czasami tkankowymi. Dla przejrzystości jednak tempo precesji na ekranie jest wizualnie przeskalowane, a nie odpowiada rzeczywistej częstotliwości Larmora rzędu dziesiątek megaherców, a przesunięcie chemiczne, dyfuzja i gradienty przestrzenne nie są modelowane. To koncepcyjne narzędzie edukacyjne, a nie symulator sekwencji klasy skanera.
Dlaczego MRI nie wykorzystuje promieniowania jonizującego?
MRI opiera się wyłącznie na silnym statycznym polu magnetycznym i impulsach radiowych do manipulowania i wykrywania jąder wodoru, głównie z wody i tłuszczu w organizmie. W przeciwieństwie do rentgena czy CT nie są zaangażowane fotony o wysokiej energii, więc nie ma dawki promieniowania jonizującego. Czyni to MRI szczególnie cennym przy powtarzanych badaniach oraz przy obrazowaniu dzieci i tkanek miękkich.
Jak kontrast relaksacyjny tworzy obraz MRI?
Poprzez taktowanie impulsów RF i odczytu sygnału skaner waży obraz w kierunku różnic T1 lub T2 między tkankami. Sekwencja ważona T1 sprawia, że tłuszcz jest jasny, a płyn ciemny, podczas gdy sekwencja ważona T2 sprawia, że płyn jest jasny. Odmienne wartości T1 i T2 wczytywane przez presety są przyczyną tego kontrastu między istotą szarą, białą, płynem mózgowo-rdzeniowym i patologią.
Co w praktyce zmienia natężenie pola B0?
B₀ ustala natężenie statycznego magnesu, tu regulowane od 0,5 do 7 T. Wyższe pole podnosi częstotliwość Larmora i zwiększa stosunek sygnału do szumu, na ogół dając ostrzejsze obrazy, ale też nasila pewne artefakty i ogranicza nagrzewanie tkanek. Skanery kliniczne to zwykle 1,5 T lub 3 T, natomiast systemy 7 T są używane w badaniach naukowych.