摘要:高频谐振腔在回旋加速器实际运行过程中受束流负载、重力和热损耗引起的腔体变形等因素的影响, 其谐振频率会发生一定的偏移, 导致高频工作频率会随着谐振腔的谐振频率而变化. 为满足等时性加速的要求, 当高频工作频率发生改变时磁场强度也应进行相应的变化, 即励磁电流的大小需要相应改变, 使得粒子回旋频率与高频谐振频率相匹配以克服滑相. 首先通过有限元仿真软件建立静磁场模型模拟不同励磁电流下回旋加速器的平均磁场, 然后理论分析磁场与谐振频率的关系, 最后得到励磁电流在小区间变化时与谐振频率的关系; 根据计算的不同谐振频率对应的最佳励磁电流, 完成励磁电流的自动跟频. 在保证最大碳膜束流的情况下, 实验得到不同谐振频率对应的最佳励磁电流, 使理论得以验证. 根据其关系实现了励磁电流自动调节, 克服了滑相, 保证了法拉第束流的稳定输出. 该方法使得励磁电流能够快速、准确地寻找并跟踪谐振腔频率, 克服了频率偏移导致的滑相, 完成束流的稳定输出.
关键词: 高频谐振腔/
谐振频率/
频率跟踪/
稳定出束

English Abstract

全文HTML

--> --> -->

正电子发射断层显像(positron emission tomography, PET)是一种先进的核医学影像技术. 在美国、欧洲和日本, ¹⁸F-FDG (氟代脱氧葡萄糖) PET已被纳入医疗保险范围^[1,2], 其实现方法是将正电子核素标记到某种参与体内代谢的物质上, 形成放射性药物(也称显像剂、放射性示踪剂), 其具有很高的灵敏度和分辨率^[3,4]. 然而只有PET是完不成成像工作的, 它需要配套回旋加速器. 在我国, 生产用于PET诊断的放射性同位素的回旋加速器几乎都是从国外进口的^[5].
回旋加速器正常工作时, 高频腔应处于谐振状态. 然而, 由于高频功率的大小、温度的升降、热损耗的差异、束流的强弱都是随时间而不停地变化的, 这导致腔体频率偏移, 不能有效地发挥加速器的性能^[6,7]. 因此, 必须采取相应的措施来消除这些不确定因素带来的影响, 以保持高频谐振腔与高频功率源始终处于稳定的匹配状态^[8].
目前, 针对射频腔体会由于腔的热不稳定性产生微小形变, 从而导致谐振腔的固有谐振频率发生变化的问题, 张天爵等^[9]、雷钰^[10]、答嘉曦等^[11,12]、刘毅^[13]、田瑞霞等^[14]、郝焕锋^[15]以及美国通用电气公司(GE) MINItrace医用回旋加速器等都是通过控制器驱动频率微调马达, 带动微调电容片前后动作, 调整腔体的谐振频率, 使其固有谐振频率与高频功率源产生的信号频率一致.
中国工程物理研究院流体物理研究所建造的紧凑型医用回旋加速器商业样机, 其引出能量为11 MeV, 引出束流为50 μA, 工作频率为72.12 MHz, 平均磁感应强度为1.18 T, 占空比可调节^[16]. 其高频结构采用非对称的半波长腔, 高频腔用于产生高频加速电场, 高频腔主要设计参数有加速电压、频率、腔壁损耗功率、品质因数等, 回旋加速器没有为谐振腔设置自动稳频装置, 从结构上简化了腔体设计, 降低了成本. 针对回旋加速器运行过程中谐振腔频率会逐渐变化的问题, 高频系统采用高频功率源的频率实时与高频腔中加速电场的频率匹配的方案. 磁铁系统为带电粒子作回旋运动提供磁场, 同时实现粒子的聚焦^[17]. 为了保证粒子的共振加速, 磁场需根据实际的高频运行频率进行调节, 即励磁电流需根据高频频率实时调节. 磁铁产生的轴向磁场须是等时性磁场, 使得被加速的离子在不同半径处的回旋频率相等, 最终在某一高频频率下得到共振加速.
本文采用数值仿真、理论分析、实验研究等方法, 对回旋加速器谐振频移的问题使用新方法得以解决.

为了实现高频功率源对谐振频率的寻找和跟踪, 在加速器启动之前, 应调节高频谐振腔处于合适的谐振频率点. 对于高频功率源工作频率与变化的腔体频率实时匹配的方法, 采用了锁相型自动频率控制(automatic frequency control, AFC)系统, 该系统高频腔的入射波信号与反射波信号之间的相位差$\Delta \varphi $是功率源输出频率f与高频腔谐振频率${f_0}$之差$\Delta f$的函数:

${\rm{tan}}\varphi \left( \omega \right) = 2{Q_{\rm{L}}}\frac{{\Delta f}}{{{f_0}}}, $

其中, ${Q_{\rm{L}}}$是有载品质因数, $\omega $是谐振角频率. 驻波锁相型AFC系统通过检测$\Delta \varphi $来控制$\Delta f$, 使其频率一致, 即达到加速器正常工作的要求.
在加速器运行过程中, 高频功率源不断改变信号源的输出频率并实时检测前向功率P₁和反射功率P₂的反馈信号, P₁与P₂的比值大小直接反映频率的匹配情况, 当0.5P₁ $ \leqslant$ P₂ $\leqslant $ 0.8P₁时, 视为严重失谐; 当P₂ < 0.5P₁, 视为轻度失谐; 当P₂ > 0.8P₁时, 判断为打火; 从而寻找高频谐振腔的谐振频率. 如果在加速器运行过程中谐振腔发生轻度失谐, 控制系统则会根据相位以寻找具有最小驻波比的频率^[18]; 如果在加速器运行过程中谐振腔发生严重失谐, 则降低占空比至10%, 升高频率, 频率升高步长为(f₁ – f₀) × 0.1 (f₁为当前频率, f₀为初始频率), 若频率调节到f₀时, 还处于严重失谐(相位变正常, 但反射未降低), 则将占空比降低为1%, 并重新启动扫描程序以确定频率; 如果在加速器运行过程中谐振腔发生打火, 那么高频功率源的控制模块将会切断高频功率源的输出, 一段时间之后, 高频功率源会以最终的工作频率再次运行. 图1—图3分别为回旋加速器开机并连续工作2 h时高频相位、高频腔压、谐振频率的实时变化曲线.
图 1 高频相位随时间的变化
Figure1. High frequency phase variation with time.

图 2 高频腔压随时间的变化
Figure2. High frequency cavity voltage variation with time.

图 3 谐振频率随时间的变化
Figure3. Resonant frequency variation with time.

从图1可以看出, 在打靶的初始阶段, 由于占空比的迅速提升, 加速器内部分部件的温度发生显著变化, 导致高频谐振频率发生变化, 控制系统会自动调整高频信号源的频率以保证谐振关系^[19], 此时相位也会发生变化, 而随着高频工作频率趋于稳定, 相位也稳定在208°左右. 从图2可以看出, 在占空比提升的初始阶段, 由于谐振频率发生偏移, 高频功率源工作频率还没稳定跟频, 而为了保证腔压, 需要强制提升激励电压, 待工作频率与谐振频率一致时, 激励电压还没下调, 因此出现短时间的波动, 之后腔压一直稳定在48 kV左右. 从图3可以看出, 当高频功率源开机, 占空比增加时, 谐振频率急剧下降, 约下降50 kHz、10 min后, 频率会基本稳定在72.52 MHz, 之后变化的范围约为10 kHz.

对于电流产生静磁场, 导体磁场强度可直接由积分得到^[20]:

$H = \int_{{\varOmega _J}} {\frac{{J \times R}}{{{{\left| R \right|}^3}}}} {\rm{d}}{\varOmega _J}, $

其中, H为磁场强度, J为$\varOmega _J $域内的体分布电流, R为从源点到场点的矢径, $ \varOmega _J$表示分布区域.
根据设计好的回旋加速器, 采用有限元软件, 建立静磁场的1/8模型, 如图4所示, 得到的后处理模型如图5所示.
图 4 静磁场模型
Figure4. Model of static magnetic field.

计算得到不同半径的平均磁场变化曲线如图6, 励磁电流与磁感应强度的曲线见图7. 由图7可知, $\Delta I$与$\Delta B$成正比, 当$\Delta I =$$ [\sum\limits_{i = 0}^n {({I_i} - {I_{i - 1}})} ]/n $= 1 A时, $\Delta B = [\sum\limits_{i = 0}^n {({B_i} - {B_{i - 1}})}]/n $ = 0.0025 T.
图 5 后处理模型
Figure5. Post-processing model.

图 6 不同半径的平均磁场变化曲线
Figure6. Average magnetic field curve at different radii.

图 7 励磁电流与磁感应强度的关系
Figure7. Relation between the excitation current and the magnetic induction intensity.

磁路中磁场强度的计算公式为

$H = \frac{{N \times I}}{{{L_{\rm{e}}}}}, $

而

$B = \mu H, $

所以由(3)和(4)式有

$B = \mu H = \frac{{N \times I}}{{{L_{\rm{e}}}}}, $

其中, H为磁场强度(A/m), B为磁感应强度(T), μ为介质的绝对磁导率(H/m), L_e为有效磁路长度(m), N为线圈数, I为励磁电流(A).
由(5)式可得, 线圈圈数确定, 有效磁路长度确定, 则

$\frac{{\Delta I}}{I} \propto \frac{{\Delta B}}{B}, $

在励磁电流很小的变化范围内, $\dfrac{{\Delta I}}{I}$与$\dfrac{{\Delta B}}{B}$呈近线性关系, 因此可得到

$\frac{{\Delta I}}{I} = \frac{1}{{199}} = k\frac{{\Delta B}}{B} = k\frac{{0.0025}}{{1.1874}},$

计算得k = 2.3867.
一个电荷量为q (C)、质量为m (kg)的粒子, 在恒定磁场以一定的速度v (m/s)在与磁场相垂直的平面上运动, 将受到磁场洛伦兹力的作用而做圆周运动, 圆周运动的曲率半径为r (m), 由下式可求得回旋加速器的回旋频率f (Hz):

$f = \frac{\omega }{{2{\text{π}}}} = \frac{{qB}}{{2{\text{π}}m}}.$

由(8)式可得, 粒子确定, 在频率小变化时, 则

$\frac{{\Delta f}}{f} \propto \frac{{\Delta B}}{B}.$

磁场发生变化必然引起相位的变化(°), 而磁感应强度B需要励磁电流I来体现. 由于由磁场引起的相移为

$\Delta {\varPhi _B} = N\frac{{\Delta f}}{f}{360^ \circ }, $

那么,

$\Delta {\varPhi _B} \propto \frac{{\Delta f}}{f} \propto \frac{{\Delta B}}{B} \propto \frac{{\Delta I}}{I},\;\frac{{\Delta I}}{I} = k\frac{{\Delta f}}{f}.$

根据励磁电流得到束流的相移曲线如图8所示, 得到的谐振频率变化曲线如图9所示.
图 8 励磁电流变化导致的相移
Figure8. Phase shift caused by change of magnet current.

图 9 谐振频率随励磁电流变化量的变化
Figure9. Relation between resonant frequency and magnet current.

从图8可以看出, 励磁电流增加0.5 A, 相移度数增大了37.8985°, 那么励磁电流增加0.01 A, 相移度数增大0.7580°; 励磁电流减少0.01 A, 相移度数减小0.7580°.
从图9可以看出, 励磁电流增加0.5 A, 谐振频率增加了76.429 kHz, 那么谐振频率每增大1.529 kHz, 励磁电流需要增大0.01 A; 谐振频率每减少1.529 kHz, 励磁电流需要减小0.01 A.

在高频功率源启动的初始阶段, 占空比提升, 谐振频率减小, 通过手动调节励磁电流, 以期得到最大碳膜束流. 通过30余次的可靠性打靶实验, 每次在保证最大束流的条件下, 记录不同谐振频率对应的励磁电流, 比较了每次的拟合曲线, 其斜率都是在1.55 kHz附近稍有浮动, 选取其中的一条实验结果, 得到的关系曲线如图10所示.
图 10 谐振频率与励磁电流的对应关系
Figure10. Corresponding relation between resonant frequency and magnet current.

此前, 为了保证靶上束流的稳定性, 通过手动调节磁场大小. 整个打靶过程中, 虽然励磁电流的变化范围约只有0.3 A, 但是对束流的影响是巨大的, 如果不去调节, 可能只产生极小的束流, 甚至可能无法产生束流; 即使励磁电流的变化范围是0.02 A, 束流也会降低约10%; 而且在束流调节过程中, 经常性的调大或调小, 导致束流不够稳定. 基于此, 为了保证束流的足够稳定, 理论分析了励磁电流与谐振频率的关系, 设计了稳定靶流下关于励磁电流与谐振频率的实验, 并设计了励磁电流跟踪谐振频率变化的稳流算法, 在连续2 h的打靶实验中, 得到了如图11所示的靶流变化曲线.
图 11 靶流随时间的变化
Figure11. Variation of target beam current with time.

从图10可以看出, 当励磁电流变化很小时, 励磁电流与高频谐振频率呈近线性关系: y = 0.1587x + 40.94, 吻合度R² = 0.9686. 在此条件下, 励磁电流每变化0.01 A, 谐振频率变化0.001587 MHz (1.587 kHz), 与理论计算结果符合得很好, 验证了计算的正确性.
图11中曲线在开始上升时, 靶流有一个下降, 这是由于占空比迅速上升, 虽然励磁电流能够实现对谐振频率的跟踪, 但是磁场响应偏慢. 在第1次到达目标束流值时, 靶流有小量的下降过程, 这是由于过快的占空比调整, 控制系统通过相位不断的调整谐振频率, 虽然励磁电流能够实现与谐振频率的匹配, 但是工作频率却并不是最佳频率, 此时最大偏差也不到稳态值的7%. 约7 min后, 腔压、相位、谐振频率、励磁电流、前向功率、反射功率等参数将稳定下来, 靶流能够很好地保持在24.2 μA左右, 误差值小于 ± 0.2 μA.
高频功率源从开机到打靶结束, 谐振频率总变化约为50 kHz, 因此励磁电流最大约需调节0.33 A, 也就是说, 在加速器实际过程中进行磁场调节时, 能量从最大值11 MeV会减小为10.985 MeV, 减小量为0.015 MeV, 变化率为0.1379%. 显然, 这是可以忍受的, 基本保证了输出能量的稳定性, 保证了束流的稳定性, 输出能量的变化曲线如图12所示.
图 12 输出能量的变化情况
Figure12. Change of output energy

本文通过高频功率源工作频率与变化的腔体频率实时匹配的方法, 实现了加速器的稳定运行, 简化了加速器的腔体结构, 节约了医用回旋加速器的建造成本; 通过仿真模拟得到了励磁电流与磁场的匹配关系; 通过实验验证了理论计算的励磁电流变化对谐振频率影响的正确性; 通过改进自动控制算法, 使得励磁电流自动调节匹配高频功率源工作频率的方法实现了回旋加速器的自动出束、自动稳束.

非常感谢回旋加速器研究组与其他参与者在问题讨论方面的帮助.