MRI用超电导磁石设计
理论：如何产生均匀磁场
概述
我以为MRI磁体就是长螺线管。简单加长不就均匀了吗？
加长确实更接近均匀，但MRI要求的均匀度完全不是一个量级。在直径40～50 cm的球内，磁场偏差要控制在几ppm到十几ppm，即百万分之几。1.5 T时1 ppm为1.5 µT。用简单螺线管满足这一点，会长到病人进不去。所以实际磁体在有限长度内排列多个线圈，调整各自的位置和匝数，依次抵消偏离中心时磁场变化的分量。
轴上磁场
半径 $a$、长 $L$、单位长度匝数 $n$、电流 $I$ 的螺线管中心轴上的磁场。
把中心附近的磁场展开为 $z$ 的幂级数，对称布置只保留偶次项。选择线圈位置和电流使 $c_2$、$c_4$、$c_6$ … 依次为零，是均匀磁体设计的基础。轴上系数确定了，球内的磁场也就确定了。
超导走进医院的那一天
MRI的原理核磁共振于1946年被发现，20世纪70年代推进了将其用于人体成像的研究。最早的装置用常导磁体，磁场也弱，但80年代用液氦冷却铌钛超导线的磁体实用化，1.5 T装置普及到医院。超导线圈电阻为零，一旦通电，即使切断电源也能持续流动多年，保持极其稳定的磁场。MRI成为超导技术规模最大的实际应用之一。
算例：用线圈布置实现均匀度
例1：简单螺线管
半径0.5 m的螺线管，求中心沿轴向±0.2 m范围内的磁场偏差：
| 长度 | ±0.2 m内偏差 |
|---|---|
| 1.0 m | 60,800 ppm |
| 1.5 m | 23,582 ppm |
| 2.0 m | 9,939 ppm |
| 3.0 m | 2,454 ppm |
例2：分成线圈对
把半径同为0.5 m的细线圈左右对称排列，第一对电流取1，其余各对的电流按从 $c_2$ 起依次抵消的原则确定：
| 布置 | 抵消的项 | 电流比 | ±0.2 m内偏差 |
|---|---|---|---|
| 3对（±0.15、0.45、0.75 m） | z², z⁴ | 1 : 0.40 : 2.80 | 1,591 ppm |
| 4对（±0.1、0.3、0.55、0.8 m） | z², z⁴, z⁶ | 1 : 0.96 : 2.02 : 0.97 | 77 ppm |
4对就降到77 ppm了。为什么两端电流更大？
线圈端部磁场减弱，所以加强外侧线圈来补偿——与在长螺线管端部多绕匝数是同一思路。每多抵消一项，偏差大约降一个数量级。实际磁体用6对以上线圈设计到几ppm，制造误差造成的偏差在安装后测量磁场，通过贴铁片（无源匀场）或给校正线圈通电（有源匀场）来修正。
例3：安匝与漏磁场
| 项目 | 值 |
|---|---|
| 1.5 T所需安匝（半径0.5 m、长1.5 m） | 2.15 MA匝（无限长近似为1.79） |
| 1.5 T磁场能量密度 | 0.895 MJ/m³ |
| 无屏蔽降到0.5 mT的距离（1.5 T，轴上） | 8.8 m |
| 同（3 T） | 11.1 m |
设计步骤
- 确定中心磁场、均匀度（球直径与ppm）、漏磁场范围、磁体长度和内径。
- 优化线圈位置和匝数，依次抵消轴上展开系数。
- 在外侧放置反向线圈抑制漏磁场（有源屏蔽）。
- 确认线圈所受电磁力，以及相对导线最大磁场和临界电流的裕度。
- 分析失超（超导态破坏）时的温度和电压，设计保护电路。
“设计5 ppm，安装后300 ppm”
在工厂确认过均匀度的MRI磁体安装到医院后偏差大幅增加，这种情况并不少见。原因是周围的磁性体——楼板钢筋、附近钢柱、隔壁房间的设备——扰乱了磁场。因此安装后要在球面上数百个点测量磁场，通过计算确定铁片布置（匀场）。在设计阶段就把房间钢结构的影响纳入分析，可以避免匀场无法校正的大偏差。
常见错误
错误与对策
| 错误 | 后果 | 对策 |
|---|---|---|
| 靠长螺线管均匀化 | 磁体太长 | 用多线圈抵消各项 |
| 忽略周围磁性体 | 安装后均匀度变差 | 把房间钢结构纳入分析 |
| 只按中心磁场选导线 | 线圈上最大磁场超过临界 | 按线圈上最大磁场确认 |
| 忽略电磁力 | 线圈移动引发失超 | 分析支撑结构 |
| 只在轴上评价漏磁场 | 误判径向范围 | 求全周等磁场线 |
我想试试相关计算。
可以用螺线管磁场计算和毕奥-萨伐尔定律模拟器确认线圈磁场，用MRI信噪比与线圈工具确认图像信噪比。超导方面参考失超分析、持续电流模式、临界电流密度和高温超导与低温超导。