1.2 布洛赫方程（Bloch Equations）

1.2 布洛赫方程：旋钮的物理法典

为什么同一台机器，换个序列就能扫出完全不同的图像？答案藏在一条方程里。1945 年布洛赫与珀塞尔两组人几乎同时独立观察到核磁共振现象，而布洛赫顺手做了一件更影响深远的事：他把观察到的磁化行为写成了一条微分方程。七十多年过去，几乎所有序列设计软件的内核里，跑的都是这条方程或它的变体。

法典的全文与逐项翻译

布洛赫方程写作 dM/dt = γM×B − (Mx·x̂ + My·ŷ)/T2 − (Mz − M₀)·ẑ/T1。看着吓人，拆开只有三件事。第一项 γM×B 是运动项：只要磁化感受到的磁场是 B，它就绕着 B 进动，转多快由 γ 决定。第二项是横向弛豫：Mx、My 以时间常数 T2 指数衰减，代表自旋之间互相"传染"失步。第三项是纵向弛豫：Mz 以时间常数 T1 向平衡值 M₀ 恢复，代表磁化把能量交还给环境（晶格）。

这三项各自对应工作台上一类旋钮。B 是你施加的一切磁场之和：静磁场 B₀、梯度磁场 G·r、射频磁场 B₁——序列里写下的每一行波形，本质都是在按时程指定 B(t)。T1、T2 不是旋钮，是每种组织自带的"材质参数"，序列设计师改不了它们，但可以通过等待时间（TR、TE、TI）决定方程解出来时各分量剩多少。换句话说：布洛赫方程的右端，一半是你写的，一半是组织定的，图像就是这两半在时间轴上博弈的结果。

转动旋钮的物理后果

先看运动项。当 B = B₀·ẑ 且无射频时，方程退化为纯进动：横向分量以 ω = γB₀ 匀速旋转，幅值不变。打开梯度 Gx 后，B 变成 B₀ + Gx·x，不同 x 位置进动频率出现差异——这一项就是第二章空间编码的全部物理基础，此处只需确认它确实来自方程第一项。再看射频项：B₁ 是旋转磁场，与进动同步（共振）时，磁化绕 B₁ 缓慢翻转，翻转角 α = γB₁τ，其中 τ 是脉冲持续时间。这就是"翻转角旋钮"的来源：要 90° 还是 180°，靠的是脉冲强度与时间的乘积。

弛豫两项则在所有时间里持续生效，与你在做什么无关。这带来一个重要推论：任何序列参数的改变，同时改变"进动路径"与"弛豫预算"两条线。把 TE 从 10 ms 拉到 100 ms，不只是"读数晚了一点"，而是横向磁化在 e^(-TE/T2) 上多衰减了一个数量级；把 TR 从 4000 ms 压到 400 ms，纵向磁化来不及回满，下一个激发拿到的本钱按 1−e^(-TR/T1) 打了折扣。旋钮从来不是单变量开关，它们都在同一条方程里记账。

对图像的影响

方程解出来的形状直接决定图像气质。短 TR 序列里，长 T1 组织（如脑脊液）恢复慢、每次被激发时本钱薄，信号被压低——这就是 T1 加权的来源之一。长 TE 序列里，短 T2 组织（如骨皮质附近的韧带）横向磁化先衰减殆尽，图像上变黑——T2 加权由此而来。而如果介质内磁场不均匀，实际失步比 T2 更快，有效时间常数是 T2*：梯度回波家族（第三章）对此格外敏感，因为它们不像自旋回波那样有"纠错脉冲"。

一个完整的参数推演案例

背景：设计一个简单实验，测量某体模的 T1。操作：选一组反转恢复式的采样点，用同一个序列只改恢复时间，分别取 TR = 200/400/800/1600/3200 ms，记录每次激发读到的信号幅度 S。结果：S 随 TR 增长，呈 S = S₀(1 − e^(−TR/T1)) 的饱和曲线，拟合得 T1 ≈ 600 ms。解读：这验证了纵向弛豫项——磁化恢复速度确实由组织自身决定，序列只是在不同时刻"拍照留证"。变式：若换成 T2 测量，则固定 TR 足够长（保证每次都满本钱），只改 TE，拟合 S = S₀e^(−TE/T2)。两个实验合起来说明：弛豫参数是材质属性，序列参数是观测安排，两者的乘积关系就是对比度的配方。

边界与坑

第一个坑：把布洛赫方程当成精确的量子描述。它是宏观经验方程，在射频功率极高或组织磁化率剧烈变化时会出现偏离，工程上够用但要知边界。第二个坑：忽略 T1 与 T2 的独立性。两者同源于分子运动，但数值可以差一个数量级，混为一谈会导致对比度推演全错。第三个坑：忘了方程是逐体素成立的。每个体素都有自己的 M、T1、T2 和局部 B，图像是方程在几十万个体素上并行求解后的快照——任何"整体上"的直觉，都要经得起"逐体素"的拷问。

从方程到仿真器

布洛赫方程的数值求解并不神秘：把时间切成微秒级小步，每步内先按旋转算子处理进动项（一个三维旋转矩阵），再按指数规律处理两种弛豫，循环往复。这正是所有布洛赫仿真器的内核，也是第八章验证链第一站用到的工具。自己写一个百行级别的仿真器是检验理解程度的最好练习：给它输入任意的射频波形与梯度时程，输出磁化三维轨迹与信号曲线。做一遍之后你会发现，时序图上每个"陡沿"在仿真里都是一段真实的数值演化，直觉与方程从此对得上账。工程上更快的路径是用旋转矩阵的解析解处理硬脉冲近似，只有软脉冲与复杂波形才回退到数值积分——这是精度与速度的经典折中。

本节要点回顾：布洛赫方程 = 进动 + 两种弛豫；旋钮改变的是进动路径与弛豫预算的交点；对比度配方 = 组织参数 × 观测安排。下一节看这些磁化如何变成可测量的电压，以及噪声从哪里来。