4. SQUID磁强计与梯度计设计

SQUID 本身只是一个对磁通极其敏感的“探测器”,真正决定系统能测到什么、测得多准的,是它与拾取线圈构成的整体结构。同一只 SQUID,配上不同的线圈几何,就得到磁强计、一阶梯度计、二阶梯度计、平面梯度计或轴向梯度计。它们的差别不在电子学,而在“线圈如何对空间磁场加权求和”。理解这一点,才能在设计阶段做出正确的取舍。

从磁强计到梯度计:结构差异的本质

磁强计使用单个拾取线圈,直接测量该位置的磁通密度。它的结构最简单、灵敏度最高,但同时也把环境磁场噪声(地磁扰动、电梯、地铁、电源线)一并接收进来。梯度计则用两个或多个反向串联的线圈,让均匀磁场在总输出中相互抵消,只保留随空间变化的磁场梯度。

设线圈匝数为 N、面积为 A,沿某一方向 z 放置。磁强计输出正比于 B(z)。一阶梯度计由两个反向串联、间距为基线 b 的线圈组成,输出正比于 B(z) − B(z+b),在小基线近似下正比于 ∂B/∂z·b。二阶梯度计则再对两个一阶梯度做差分,输出正比于 ∂²B/∂z²·b²。每提高一阶,对均匀场的抑制就增强一次,但对真实源信号的衰减也更快。

平面梯度计与轴向梯度计的区别在于线圈法线方向与基线方向的关系。轴向梯度计的两个线圈共轴、沿同一法线方向排列,基线沿轴向,对轴向梯度最敏感;平面梯度计的两个线圈位于同一平面内、法线平行,基线在平面内,测量的是横向梯度。平面结构便于用薄膜工艺在基片上集成,适合多通道阵列;轴向结构通常用绕线或叠层实现,基线可以做得更长,但通道密度受限。

磁强计与梯度计线圈结构对比 磁强计 单线圈 · 输出 ∝ B 灵敏度最高 环境噪声全接收 一阶梯度计 b 反向串联 · 输出 ∝ ∂B/∂z·b 抑制均匀场 二阶梯度计 输出 ∝ ∂²B/∂z²·b² 抑制更强 · 信号衰减大 平面梯度计 基线在平面内 测横向梯度 薄膜集成 · 多通道 轴向梯度计 b 共轴 · 基线沿轴向 共模抑制与基线权衡 基线 b 小 基线 b 大 共模抑制比 CMRR ↑ 空间分辨率 ↓ 基线增大 → 均匀场抑制增强,但源信号空间平均效应加剧 需在 CMRR 与分辨率之间取平衡点

共模抑制比:梯度计的核心指标

共模抑制比(CMRR)衡量梯度计对均匀磁场的抑制能力,定义为均匀场下的输出与梯度场下输出的比值,通常用 dB 表示。理想一阶梯度计对均匀场的响应为零,但实际线圈的匝数、面积、位置不可能完全对称,残余不平衡会直接限制 CMRR。

设两个线圈的等效面积分别为 A₁、A₂,理想情况下 A₁ = A₂。若存在面积失配 ΔA = A₁ − A₂,则均匀场 B₀ 产生的残余输出正比于 ΔA·B₀,而梯度信号正比于 A·b·∂B/∂z。于是 CMRR 近似为:

CMRR ≈ 20·log10( A·b / ΔA )   [dB]

示例:A = 100 mm²,b = 50 mm,ΔA = 0.01 mm²
CMRR ≈ 20·log10( 100×50 / 0.01 ) = 20·log10(5×10⁵) ≈ 114 dB

这个示例说明,面积失配必须控制在万分之一量级,才能获得 100 dB 以上的 CMRR。实际系统中,线圈绕制误差、基片厚度不均、引线耦合都会引入失配。因此高 CMRR 梯度计往往需要可调平衡结构,例如在其中一个线圈上并联微调电容或增加补偿匝,通过实测均匀场下的残余输出进行校准。

校准方法:将梯度计置于均匀场线圈(Helmholtz 线圈)中,施加已知均匀场,调节平衡元件使输出最小。记录残余输出与施加场之比,即为实际 CMRR。建议在系统安装后和每次更换杜瓦位置后重新校准。

基线选择与空间分辨率权衡

基线 b 是梯度计两个线圈之间的间距,它同时影响两个关键指标。增大 b 会提高对均匀场的抑制能力,因为差分信号对空间变化的敏感度随 b 增大;但 b 增大也意味着线圈对源信号的空间平均范围变大,空间分辨率下降。对于靠近探测器的浅层源,大基线会把相邻特征“糊”在一起;对于深层源,适当增大基线反而能提高信噪比。

选择基线时,通常先估计目标源的深度 d 和空间尺度 L。经验上,一阶梯度计基线取 b ≈ 0.5d ~ d 较为合理;二阶梯度计基线可略小,因为二阶差分本身已提供更强抑制。若目标是心磁图(源深度约 100 mm),一阶轴向梯度计基线常取 50~70 mm;若目标是脑磁图(源深度约 30~50 mm),基线常取 20~40 mm,并优先采用平面梯度计以贴近头皮。

类型典型基线CMRR 量级空间分辨率适用场景
磁强计——最高屏蔽室、低噪声环境
一阶轴向梯度计50~70 mm80~120 dB中等心磁、地质勘探
一阶平面梯度计10~30 mm60~100 dB较高脑磁、多通道阵列
二阶梯度计30~50 mm100~140 dB较低强干扰环境、移动测量

需要强调的是,梯度计并不能“消除”噪声,它只是改变了对不同空间频率磁场的响应权重。均匀场被抑制,但梯度噪声仍然存在。因此在高梯度噪声环境中,二阶梯度计未必优于一阶;此时更有效的做法是改善屏蔽或采用参考通道做自适应抵消。

常见误区:认为梯度计阶数越高越好。实际上每增加一阶,对源信号的衰减约增加一个 b 因子,若源信号本身较弱,高阶梯度计可能把有用信号压到本底噪声以下。设计时应先估算源信号强度与噪声谱,再决定阶数和基线。

设计流程与参数确定

实际设计可按以下顺序推进:首先明确目标源深度、信号频带和预期噪声水平;然后根据噪声的空间相关性选择梯度计阶数——若环境噪声以均匀场为主,一阶即可;若存在明显梯度噪声,考虑二阶或参考通道。接着确定基线:用源深度估算最佳 b,并校核空间分辨率是否满足要求。最后进行线圈几何设计,计算匝数、面积和失配容差,并预留平衡调节结构。

以心磁测量为例,假设源深度 d ≈ 100 mm,目标信号为 10~100 pT 量级的 MCG 信号,环境噪声在屏蔽室内约 1~10 pT/√Hz。选择一阶轴向梯度计,基线 b = 60 mm,线圈直径 20 mm,匝数 20 匝。估算面积 A ≈ 314 mm²,若要求 CMRR ≥ 100 dB,则面积失配需满足 ΔA/A ≤ b/(10⁵) ≈ 6×10⁻⁴,即失配面积小于约 0.19 mm²。这一容差在手工绕制中较难保证,因此需要加入平衡调节。

完成设计后,应通过均匀场测试验证 CMRR,通过已知梯度场测试验证灵敏度,并在实际环境中对比磁强计与梯度计的噪声谱。若梯度计噪声谱在低频段显著低于磁强计,说明共模抑制有效;若两者接近,则说明残余失配或梯度噪声占主导,需要重新调整平衡或改变基线。