核磁共振成像MRI原理
关键词
核磁共振; 自旋; 拉莫尔进动; 射频脉冲; 弛豫时间T1和T2; 梯度磁场; k空间; 傅里叶变换; 超导磁铁; 功能MRI
"核"字不等于核辐射——MRI 扫描用的是无线电波
许多患者在听到"核磁共振"中的"核"字时感到恐惧,以为与核辐射有关。这是一个严重的误解:MRI(Magnetic Resonance Imaging)中的"核"指的是原子核(特别是氢原子核,即质子),而非核辐射。MRI既不使用X射线,也不产生任何电离辐射——它利用的是无线电波(射频电磁波,频率约10-300 MHz,属于非电离辐射)和强磁场(通常1.5-3特斯拉,强度约为地球磁场的3万至6万倍)来激发和探测人体中氢原子核的自旋信号。
MRI的信号来源是人体中大量的氢原子(主要以水分子和脂肪分子中存在)。成年人体约60%由水组成,意味着氢原子核(质子)无处不在,是MRI天然的"造影剂"。不同组织(肌肉、脂肪、脑脊液、骨髓)中氢原子的浓度和化学环境不同,在MRI图像中呈现不同的对比度,使MRI成为软组织成像的最佳手段,远优于X射线CT(CT对软组织对比度低)。
射频光子的能量远低于可见光,更到不了辐射防护里用来划分电离与非电离的约 10 电子伏特工作门槛。MRI 的安全议题是强磁场对铁磁性物体的吸力、梯度场引起的周围神经刺激,以及射频的组织加热(比吸收率 SAR),不是电离径迹。把「核」读成核辐射,会把完全不同的相互作用叠在同一个恐惧词上。
氢核几乎遍布软组织,所以对比度首先来自水和脂肪所处的化学环境,而不是像 CT 那样来自电子密度。同一套自旋物理,化学家不加梯度时看的是化学位移,见核磁共振波谱;医生叠上线性梯度后,频率差被翻译成位置。临床成像的物理骨架见核磁共振,本文只把梯度场这一步接到人体。
从化学家的波谱仪到全身成像:布洛赫、劳特伯与那场诺奖风波
核磁共振(NMR,Nuclear Magnetic Resonance)由费利克斯·布洛赫(Felix Bloch)和爱德华·珀塞尔(Edward Purcell)在1946年独立发现,两人因此共享1952年诺贝尔物理学奖。NMR最初是化学家的工具:不同化学环境中的氢原子核(和其他核)有略微不同的共振频率(化学位移),用于分析分子结构——今天,NMR波谱仍是有机化学和生物化学中不可缺少的分析工具(每台高场NMR波谱仪价值数百万美元)。
布洛赫在斯坦福用感应线圈听进动核的磁场,珀塞尔在哈佛测射频被共振吸收的功率,两条实验线在 1946 年几乎同时听见同一件事。他们的对象还是均匀样品里的核磁矩,没有空间编码;劳特伯的梯度场是后来才把「谱」变成「图」。化学位移至今仍是有机化学的主力,成像只是在同一套拉莫尔进动上多加了三个梯度线圈。
将NMR用于人体成像的关键洞察来自保罗·劳特伯(Paul Lauterbur)和彼得·曼斯菲尔德(Peter Mansfield)。劳特伯在1973年证明,通过施加一个随空间线性变化的磁场(梯度磁场),不同位置的质子将有不同的拉莫尔进动频率,从而可以对NMR信号进行空间编码,重建出二维图像。曼斯菲尔德在1970年代开发了快速成像的数学方法(平面回波成像,EPI),使临床实用的MRI成为可能。两人共享2003年诺贝尔生理学或医学奖。
劳特伯 1973 年那篇 Nature 短文的关键不是又一台更灵敏的波谱仪,而是承认:外加随空间线性变化的磁场后,不同位置的质子会唱不同的拉莫尔频率。他把这种方法叫做 zeugmatography(轭合描绘),把化学信息与空间坐标轭在一起。曼斯菲尔德几乎同时讨论固体里的「NMR 衍射」,后来的平面回波成像把一次激励后的 k 空间快采变成临床可承受的扫描时间。
雷蒙德·达马迪安(Raymond Damadian)在1977年建造了第一台全身MRI机器,用于人体成像。达马迪安最初的动机是发现肿瘤组织的弛豫时间与正常组织不同,提出用NMR检测癌症。争议的是,诺贝尔委员会将2003年的奖项授予劳特伯和曼斯菲尔德而未包括达马迪安,后者在各大报纸上发表整版广告表示抗议——这是诺贝尔奖历史上少有的公开争议案例。
达马迪安 1971 年在 Science 上报告肿瘤组织与正常组织的弛豫时间不同,1977 年 7 月用自制全身磁体 Indomitable 得到第一幅人体胸腔图像。诺贝尔委员会把奖项授予空间编码与快速成像的物理,而不是「发现肿瘤弛豫差」这一条动机,争议由此公开化。奖状漏掉谁,并不取消梯度场把频率变成位置这一步在物理上的地位。
拉莫尔进动、T1/T2 弛豫与傅里叶变换:质子如何被编成一张图
氢原子核(质子)是自旋 $I = 1/2$ 的粒子,具有固有磁矩 ,其中 是质子的回旋磁比(gyromagnetic ratio), rad/(s·T)。在外部主磁场 中,质子磁矩沿 方向排列,并以拉莫尔频率 进动(precession)。对于1.5T的MRI机器,质子拉莫尔频率约63.9 MHz;对于3T机器,约127.7 MHz——这正是射频激励脉冲的频率。
1.5 T 对应大约 64 MHz,3 T 对应大约 128 MHz,都落在调频广播与电视之间的射频窗口,所以机房要屏蔽的是无线电干扰,不是 X 射线。回旋磁比是核的属性,梯度场只是让不同体素唱不同的调;选错中心频率,整层都不会共振。临床把场强往上推,信噪比上升,射频加热与磁化率伪影也一起上升,这是同一条拉莫尔账的两面。
梯度线圈快速切换时的洛伦兹力会让梯度罩振动,机房里那种锤击般的噪声由此而来;它是机械耦合,不是射频在“叫”。听力保护与静音序列改的是振动路径和斜率,不是拉莫尔公式本身。
射频激励脉冲(RF pulse)将质子的净磁化矢量从平衡方向(沿 )翻转(例如翻转90°)。停止射频后,质子磁化矢量恢复平衡的过程(弛豫)产生感应信号。T1(纵向弛豫时间)描述磁化矢量纵向分量的恢复速率;T2(横向弛豫时间)描述横向磁化矢量的衰减速率(来自不同质子进动频率的相位扩散)。不同组织的T1和T2值不同,这是MRI软组织对比度的根源。
T1 加权突出脂肪与脑白质的纵向恢复差,T2 加权让脑脊液等横向磁化衰减较慢的组织发亮。对比度是弛豫参数的函数,不是「水越多就越白」的简单照片。射频把净磁化矢量扳离 之后,自由感应衰减和自旋回波把横向磁化变成线圈里的电压;脉冲序列用重复时间 TR 和回波时间 TE 去放大组织之间通常只有数百毫秒量级的差别。没有这两套时间,梯度场只能给出位置,给不出医生要看的软组织对比。
梯度磁场的三个方向(x、y、z)分别负责层面选择(确定成像层面)、相位编码和频率编码,将三维空间中的NMR信号编码为k空间(k-space)数据。k空间是MRI信号的傅里叶空间表示,通过傅里叶逆变换重建出实空间图像。这就是为什么MRI重建需要大量的计算——一次全脑扫描的原始数据重建涉及数千次傅里叶变换。
k 空间里每一个点都是整幅图像的一条空间频率,不是图像上的一个像素。欠采样会按傅里叶变换的规则长出混叠伪影,加速扫描的每一条捷径都是在与采样定理做交易,见傅里叶分析。选层、频率编码和相位编码都是在改局部磁场,从而改进动频率或相位,拉莫尔公式 在三个方向上被用了三次。
液氦告急、fMRI 的死鲑鱼实验:超导磁铁背后的现实代价
临床MRI机器的主磁场(1.5T-7T)几乎全部由超导磁铁产生。低温超导线圈(铌钛合金,Nb-Ti,临界温度约9.8K;或铌三锡,Nb3Sn,约18K)浸泡在液氦(沸点4.2K)中,以极低温维持超导状态,允许大电流无损耗流动,产生强而极稳定的磁场(磁场漂移约每百万分之一每小时)。液氦一旦泄漏("quench",失超),大量液氦瞬间汽化,产生巨量气体,是严重安全隐患。
液氦是氦气(He)在地球上无可再生的资源——主要来自天然气开采的副产品,全球储量有限。氦价格在2010年代多次暴涨,引发了MRI设备供应危机的担忧。解决路径是开发"无液氦"或"液氦零蒸发"MRI系统:使用高温超导材料(如钇钡铜氧YBCO,临界温度约90K,可用液氮冷却)或高效制冷机直接冷却到4K以下(GM制冷机)。西门子和GE已推出采用密闭液氦循环系统(无需补液氦)的MRI机器。
超导磁体把大电流锁在液氦温区的线圈里,才有临床 1.5 T 与 3 T 那种既强又稳的 。失超时磁能转成热,液氦暴沸,需要紧急排气;这是超导体的工程代价,不是核辐射泄漏。零蒸发与密闭循环是在对付氦资源,不是改写拉莫尔公式。
功能磁共振成像(fMRI,functional MRI)利用血氧水平依赖(BOLD,Blood Oxygen Level Dependent)信号间接探测脑神经活动:神经元激活时,局部血流增加,氧合血红蛋白/脱氧血红蛋白比例变化,两者对T2\*弛豫时间的影响不同,产生可探测的信号差异。fMRI是过去30年认知神经科学最重要的实验工具,已发表数万篇论文。但fMRI的时间分辨率(约秒级,相比神经脉冲的毫秒级)和空间分辨率(约毫米级)仍有限制,且许多早期fMRI研究的统计方法存在争议(如"死鲑鱼"实验揭示的多重比较问题)。
BOLD 对比来自脱氧血红蛋白的顺磁性对 T2\* 的扰动,Ogawa 1990 年把这件事写成可成像的信号。它测的是血流代谢的下游,时间分辨率停在秒级,不能当成神经元放电的直接录音。2009–2010 年 Bennett 等人把一条已死亡的大西洋鲑鱼推进扫描仪,在未校正的多重比较下也能标出「激活」,提醒的是统计门槛,不是自旋物理失效。
静息态与任务态设计都仍受血流动力学延迟约束:神经事件在毫秒,BOLD 峰常在数秒后才到。解读相关图时,先问延迟与对照是否匹配,再谈脑区“点亮”。
全球 5 万台机器、每年 1 亿次检查:MRI 如何成为医学的基础设施
MRI今天已成为现代医学不可缺少的诊断工具。全球运行的MRI机器约5万台(2022年),每年进行超过1亿次MRI检查(美国约4000万次)。MRI在神经科(脑肿瘤、脑卒中、多发性硬化症)、骨肌肉系统(关节、韧带、椎间盘)、心脏(心肌病、心脏结构)、腹部(肝脏、胰腺、肾脏)等领域是不可替代的诊断手段。
临床普及靠的是软组织对比,不是把磁场当辐射源来比谁更强。脑、关节、脊髓这些 CT 看不清的地方,T1/T2 把水与脂肪的化学环境写成灰度;选型是物理问题,装机与年检查量仍沿用上段约 5 万台、每年超过 1 亿次,不另造市场数字。
定量 MRI 试图把 T1、T2、扩散系数从“看起来像”推进到可复现的物理量;同一台机器换线圈或换序列,标定不当就会把组织参数测歪。技术进步改的是信息密度与可比性,不是把年检查量再数一遍。
超高场MRI(7T及以上)正在从科研走向临床。7T MRI可提供远高于3T的图像分辨率和信噪比,使皮层柱(cortical column,约500微米宽)的成像成为可能,开启了脑科学研究的新精度。FDA 于 2017 年批准西门子 MAGNETOM Terra 作为首台临床 7 T 系统。法国 CEA NeuroSpin 的 Iseult 装置把人体成像场强做到 11.7 T,2024 年公布活体脑图像;荷兰奈梅亨另有 14 T 全身系统的规划,二者不是同一个项目。
同时,便携式低场MRI(0.05-0.5T,使用永磁体或低场超导,无需固定基础设施)正在走向重症监护室、手术室和发展中国家,有望扩大MRI的普及范围。2020 年 Hyperfine Swoop 获 FDA 许可,永磁约 0.064 T,插墙电、推到床边,不需要液氦。低场换来的是床旁可达,不是 7 T 那种皮层柱分辨率;全球装机与年检查量仍以本节前述约 5 万台、每年超过 1 亿次为口径,此处不另造数字。
连接节点:量子力学诠释 → 核磁共振成像MRI原理(MRI依赖量子自旋概念);核磁共振成像MRI原理 → 超导体(超导磁铁是临床MRI的核心器件)
事实卡
- 卡1
MRI中"核"指原子核(质子),与核辐射无关;它使用无害射频电磁波和强磁场(1.5-3T),无电离辐射。
- 卡2
质子拉莫尔频率与磁场成正比:3T场中约127.7 MHz(射频波段),T1和T2弛豫时间决定软组织对比度。
- 卡3
劳特伯和曼斯菲尔德因MRI的空间编码原理(梯度磁场+傅里叶变换)获2003年诺贝尔生理学或医学奖。
- 卡4
全球运行约5万台MRI机器,每年进行超过1亿次检查,是现代医学最重要的软组织成像手段。
引用
"但愿所有的弯路都能如此富有成果!"(All detours should be so productive!)
跨域连接
- 超导体:临床 MRI 的强磁场几乎全部来自液氦冷却的超导线圈:零电阻允许大电流无损流动,磁场稳定度达到每小时百万分之一量级。可检验推论:一旦失超,储能瞬间转为热量、液氦暴沸——磁场质量与安全风险是同一枚超导硬币的两面。
- CT 与 MRI:两种成像的物理探针完全不同:CT 测 X 射线衰减,MRI 测质子自旋的射频信号,后者无电离辐射。推论:软组织对比来自弛豫时间差异而非密度差异,所以 MRI 在脑与关节上远优于 CT,而 CT 在骨骼与急诊速度上反超——选型是物理问题,不是档次问题。
- 核磁共振波谱:同一种物理,化学家用它看分子结构(化学位移),医生用它看人体解剖(空间编码)。分野只在梯度场:不加梯度,频率差异来自化学环境;叠加线性梯度,频率差异来自空间位置。推论:图像分辨率上限由梯度强度决定,正如波谱分辨由磁场均匀度决定。
- 傅里叶分析:MRI 的空间编码是傅里叶变换的直接工程化:梯度场把位置翻译成频率与相位,采满 k 空间后一次逆变换重建图像。推论:欠采样会在图像上产生可预期的混叠伪影——加速扫描的每一条捷径,都是在与采样定理做交易。
- 统计学:功能磁共振用血氧信号间接追踪脑活动,全脑数万个体素同时检验,多重比较问题随之而来——"死鲑鱼"实验正是用一条死鱼扫描出"显著激活"来示警。推论:不做严格的多重比较校正,噪声必然伪装成功能激活,这是早期相关文献可重复性危机的统计学根源。
参考文献
- Lauterbur, Paul C. "Image Formation by Induced Local Interactions: Examples Employing Nuclear Magnetic Resonance." Nature 242 (1973): 190-191.
- Mansfield, Peter, and P. K. Grannell. "NMR 'Diffraction' in Solids?" Journal of Physics C 6 (1973): L422-L426.
- Haacke, E. Mark, et al. Magnetic Resonance Imaging: Physical Principles and Sequence Design. Wiley-Liss, 1999.
- Ogawa, Seiji, et al. "Brain Magnetic Resonance Imaging with Contrast Dependent on Blood Oxygenation." PNAS 87 (1990): 9868-9872.
- 侯中军, 郑海荣. 《磁共振成像原理与技术》. 科学出版社, 2015.
- Bloch, Felix. "Nuclear Induction." Physical Review 70 (1946): 460-474.
- Purcell, Edward M., H. C. Torrey, and Robert V. Pound. "Resonance Absorption by Nuclear Magnetic Moments in a Solid." Physical Review 69 (1946): 37-38.
- Damadian, Raymond. "Tumor Detection by Nuclear Magnetic Resonance." Science 171 (1971): 1151-1153.
- Bennett, Craig M., Abigail A. Baird, Michael B. Miller, and George L. Wolford. "Neural correlates of interspecies perspective-taking in the post-mortem Atlantic Salmon: An argument for multiple comparisons correction." Journal of Serendipitous and Unexpected Results 1 (2010): 1-5.
- Boulant, Nicolas, et al. "In vivo imaging of the human brain with the Iseult 11.7-T MRI scanner." Nature Methods 21 (2024): 2013-2016. DOI: 10.1038/s41592-024-02472-7