Development of dynamic power control system in proton therapy facility
-
摘要: 针对质子治疗装置中主环动态电源多平台能量的引出需求，研制了基于开源平台的高速实时动态电源控制系统，该控制系统以开源平台Beaglebone作为顶层硬件接口，以现场可编程逻辑门阵列（FPGA）为核心的控制器作为底层硬件接口，采用分布式的实验物理及工业控制系统（EPICS）进行远程控制。该系统可实时传输任意动态电源的输出参考电流波形数据，并结合定时系统与联锁系统，控制动态电源按预设电流波形进行输出，并实现多平台能量的引出。实验结果显示该控制系统能实现每秒最高十万条指令传输，百万次数据传输零误码率。同时，该系统结构灵活、扩展性强，可作为通用控制平台。
-
关键词:
- 质子治疗 /
- 电源控制 /
- Beaglebone /
- GPMC /
- FPGA
Abstract: A high-speed real-time dynamic power control system based on an open source platform has been developed to satisfy the multi-platform energy extraction requirements of the main loop dynamic power supply in a proton therapy facility. The control system uses a Beaglebone based open-source platform as the top-level hardware interface, the controller with Field Programmable Gate Array (FPGA) as the core as the underlying hardware interface, and uses a distributed Experiment Physics and Industrial Control System (EPICS) for remote control. This system can transmit the output reference current waveform data of any dynamic power source in real time, and controls the dynamic power source to output according to the preset current waveform combined with the timing system and the interlocking system, and realizes the energy extraction of multiple platforms. Experimental results show that the control system can achieve a maximum of 100 000 instructions per second and zero error rates for millions of data transmissions. In addition, the system structure is flexible and extensible, and it can be served as a universal control platform.-
Key words:
- proton therapy /
- power supply control /
- Beaglebone /
- GPMC /
- FPGA
-
新型核燃料的研发过程中，必须对燃料进行辐照考核。在研究堆内建立模拟真实压水堆运行环境的高温高压辐照考验回路，可以为辐照考核提供广泛、全面、真实、可信的试验环境。目前国内外许多研究堆都配套建设有堆内高温高压试验回路，例如中国高通量工程试验堆(HFETR)的500 kW高温高压回路[1]，比利时研究堆(BR2)内的CALLISTO回路[2]，韩国HANARO研究堆内的高温高压回路等[3]。辐照考验回路因为其本身的特殊性，须对其进行安全分析。张毅对中国先进研究堆(CARR)内回路进行了失流事故及失水事故的系统分析[4]；H.Ahmed针对CALLISTO回路，计算了辐照装置内压力管大破口失水事故瞬态；H.Bae等人对不同破口位置的失水事故进行了分析[5]。回路基于研究堆而建，在安全分析中不仅要考虑回路系统，还应考虑辐照装置与研究堆，但目前将系统分析与辐照装置及研究堆相结合的研究较少。本文介绍基于中国绵阳研究堆(CMRR)的堆内高温高压辐照考验回路初步设计方案，并对回路失水事故(冷/热管段失水事故)和失流事故(全部失流/卡轴事故)两类典型设计基准事故进行计算分析，以验证回路及研究堆在这两种设计基准事故下的安全性。
1. CMRR堆内辐照考验回路设计方案
CMRR堆内高温高压辐照考验回路初步设计方案由中国工程物理研究院完成。在设计过程中，主要考虑了以下几个方面：(1)回路的稳态运行环境应与真实压水堆主回路运行环境保持一致；(2)回路在设计基准事故下的安全性应得到验证，燃料棒不发生熔化；(3)回路不应对研究堆产生较大影响，研究堆冷却水不发生沸腾。
回路设计方案如图 1所示。设计方案包括三条回路。中间回路(二回路)作用是将主回路热量传递至三回路，三回路连接冷却塔，最终将热量通过冷却塔导出。主回路包含堆内辐照装置，放置于CMRR竖直孔道内，内部放有16根待辐照考验的燃料棒。主回路其他主设备还有：稳压器，控制主回路压力；热交换器，将主回路热量传递给中间回路；一台主泵及一台备用泵；电加热器，用以在回路运行前加热回路冷却剂，抬升冷却剂温度平台；安注箱及低压安注系统；化容控制系统。主回路主要设计参数见表 1。
表 1 主回路主要设计参数Table 1. Designed value of key parameters of primary looppower/kW primary loop pressure/MPa argon gap pressure/MPa coolant flow rate/(kg·s-1) inlet coolant temperature/K outlet coolant temperature/K 60 15.5 0.165 3.5 553 583 辐照装置俯视图如图 2所示，其主体由两根同轴压力管组成。双层压力管结构内部间隙填充氩气，填充压强为0.165 MPa。相比于CALLISTO回路内辐照装置填充的氦气[2]，氩气具有更低的热导率[6-7]，可以充当绝热层，以减小内压力管热应力和对研究堆冷却剂的影响。填充气隙还可以用来监视内压力管的完整性。内压力管破裂会导致气隙压强迅速升高，通过监测间隙压强变化可以判断内压力管完整性。设计时考虑到上述因素，在间隙上方连接有一安全阀，当压力超过1.0 MPa时，自动开启安全阀向水箱泄压以降低间隙压强，从而保证外压力管的完整性，避免辐照装置破损对研究堆造成影响。
辐照装置与主回路的接口均设在装置顶端。冷却剂自主泵流向辐照装置后，从内压力管与围板之间向下流动，在辐照装置底部转向后从围板内向上流经燃料棒表面，最后从顶部流出辐照装置。
2. RELAP5建模
图 3为高温高压回路的RELAP5节点图。节点图包含主回路主要设备，包括辐照装置、稳压器、热交换器、主回路主泵及管道等，电加热系统仅在回路功率运行前开启以抬升冷却剂温度平台，在功率运行及事故过程中无动作，无需对其建模。二回路和三回路不重点考虑。其中辐照装置的建模包含了内压力管和外压力管、氩气间隙、安全阀和竖直孔道内研究堆冷却水。破口用一个触发阀模拟[8]，在事故发生0 s时开启。触发阀分别放置在冷管段和热管段，模拟两种不同位置的破口事故。
系统稳态参数见表 2。RELAP5计算值与设计值近乎一致，在事故瞬态计算中是可以接受的。
表 2 系统稳态参数Table 2. Key parameters in steady-statepower/ kW pressure/ MPa inlet coolant temperature/℃ outlet coolant temperature/℃ main loop coolant flow rate/(kg·s-1) argon gap pressure/ MPa cladding temperature/ K design value 600 15.500 280.00 310.00 3.500 0.165 630.00 RELAP5 calculation value 600 15.556 278.89 310.41 3.476 0.175 640.63 辐照装置放置在研究堆孔道内，其裂变功率降低依托研究堆负反应性引入。当回路停堆信号触发后，研究堆停堆信号同时被触发，研究堆停堆系统介入，反应堆达到次临界状态，辐照装置裂变功率也随之降低。
3. 小破口失水事故计算
小破口失水事故瞬态计算基于以下假设：(1)回路在0 s前处于100%功率运行状态，0 s时发生小破口失水事故；(2)触发稳压器低压力信号后0.3 s，停堆信号触发；(3)事故发生后，所有设备失去厂外电，主泵与主给水系统停止运行。
0 s前，回路处于600 kW稳态运行状态，小破口失水事故在0 s时发生。回路破口流量如图 4所示，冷管段失水事故破口流量在瞬间达到峰值4.844 kg/s，而热管段破口流量峰值为3.686 kg/s。这是由两方面原因造成的：一是冷管段破口位置在主泵后，压力较高，与外环境压差更大，导致流量更高；二是热管段冷却剂温度更高，破口处闪蒸现象来得更早，导致流动阻力增加，流量较小[9]。
图 5为稳压器压力变化图。事故发生后稳压器压力迅速降低，触发停堆信号，随后触发中压安注系统，下降趋势变缓，最后低压安注系统介入，主回路压力保持稳定。
整个过程中，热管段事故破口流量几乎均小于冷管段事故破口流量。在热管段事故中，冷却剂要先流经辐照装置才会从破口流出，而冷管段事故中冷却剂反而是从辐照装置流出的。因此，热管段事故辐照装置水位要比冷管段事故水位高，如图 6所示。
图 6 失水事故辐照装置水位Figure 6. Coolant level in the irradiation facility during loss-of-coolant accident图 7是包壳热点温度变化图。冷管段事故中包壳温度在事故发生后迅速升高，约2.8 s时达到峰值880.6 ℃，这是因为破口位置处于冷管段，冷却剂从破口大量流出，流动方向与稳态运行时相反，在辐照装置内会发生流量反转现象，导致燃料棒产出余热无法被带走，包壳温度升高。热管段事故则没有类似现象，包壳温度在事故发生后迅速降低。这是因为破口发生在热管段，冷却剂流出方向与稳态相同，会增加通过辐照装置内部流量，使得导出热量增加，包壳温度降低。整个过程中包壳温度低于安全限值1204 ℃[10]，符合验收准则。
图 8为氩气间隙压力变化图。冷管段事故中压力在事故开始时有小幅度上涨，这是由冷却剂温度上升引起的。热管段事故中则没有这种现象。图 9为CMRR孔道冷却剂温度变化图，两种失水事故下CMRR孔道冷却剂温度都不会发生明显上升。CMRR是池式研究堆，要求冷却剂不发生沸腾，冷却剂温度限值为90 ℃。显然在这两类事故中研究堆冷却剂温度均低于安全限值，可以保证研究堆安全。
图 9 失水事故CMRR冷却水温Figure 9. Temperature of China's Mianyang Research Reactor coolant during loss-of-coolant accident4. 失流事故
失流事故是所计算的另一类设计基准事故。本章节计算了两种事故，分别是全部失流事故和卡轴事故。事故发生时，主回路主泵停止运转，回路冷却剂失去驱动力。两种事故区别在于，卡轴事故中主泵没有发生惰转，而失流事故中有主泵惰转。计算结果验证了回路在失流事故中的安全性。
图 10是在事故初期50 s内通过辐照装置的冷却剂流量图。全部失流事故中主泵发生堕转，冷却剂流量下降速度较为缓慢，在21.8 s时降低为2.85 kg/s，触发低流量停堆信号，功率迅速下降。卡轴事故中流量在事故发生后迅速下降，极短时间内触发了低流量停堆信号，回路功率下降，如图 11所示。
图 10 卡轴事故与失流事故主回路冷却剂流量Figure 10. Flow rate through the irradiation facility during loss-of-flow accidents (LOFAs)
图 11 卡轴事故与失流事故回路功率Figure 11. Power during clamp shaft accident and loss-of-flow accidents图 12为稳压器压力变化图。全部失流事故发生后回路冷却能力降低，冷却剂出口温度上升，稳压器压力上涨。稳压器压力在22.6 s时达到峰值15.624 MPa，随后由于功率迅速降低，稳压器压力回落。反观卡轴事故，由于通过辐照装置流量下降迅速，导致回路短时间内几乎丧失冷却能力，因此辐照装置出口冷却剂温度增高明显，稳压器压力上升幅度大。可以看到稳压器压力在事故开始后大幅上升，在6.4 s达到峰值压力16.21 MPa。随后由于功率降低，稳压器压力回落。整个过程中压力峰值均低于110%设计压力。
图 12 卡轴事故与失流事故主回路压力Figure 12. Pressure in the main loop during clamp shaft accident and loss-of-flow accidents全部失流事故被归为II类工况，而卡轴事故则被归为IV类工况，二者验收准则有所不同[10]。图 13为全部失流事故中DNBR变化图，整个过程中MDNBR大于1.5，出现在事故发生后22.0 s，意味着没有发生偏离泡核沸腾，且回路压力低于110%设计压力，符合验收准则[10]。图 14为卡轴事故包壳热点温度变化图。在事故发生后2.4 s，包壳温度达到峰值，约为734.1 ℃，随后由于功率下降及回路自然循环带走余热，包壳温度回落。整个过程中包壳热点温度小于温升极限1482 ℃[10]，且回路压力低于120%设计值，满足验收准则。
5. 结论
本文基于CMRR堆内高温高压辐照考验回路初步设计方案，建立了RELAP5模型，模型稳态运行关键参数与设计值保持一致。在此基础上进行了冷/热管段失水事故及卡轴/失流事故瞬态计算。结果表明，冷/热管段失水事故中包壳热点温度最高分别为880.6 ℃及367.6 ℃，低于限值1204 ℃；CMRR孔道内冷却水温变化幅度很小；双层压力管间隙压力在冷管段事故中涨幅很小，不会对双层压力管结构产生冲击；全部失流事故中MDNBR大于1.5，燃料元件表面不会发生偏离泡核沸腾；卡轴事故中包壳热点温度最高为734.1 ℃，低于1482 ℃，回路压力低于120%设计压力。计算结果表明，回路设计方案能应对所计算事故瞬态，回路与CMRR的安全都能得到保障。
致谢: 感谢中国工程物理研究院核物理与化学研究所202科室郭斯茂老师及王冠博老师给予的指导和帮助。 -
[1] 汤启升, 王志山, 李雪军. 上海质子治疗装置同步环真空布局及真空室设计[J]. 原子能科学技术, 2015, 49(s2):529-532. (Tang Qisheng, Wang Zhishan, Li Xuejun. Vacuum layout and chamber design for synchrotron ring of Shanghai Proton Therapy Facility[J]. Atomic Energy Science and Technology, 2015, 49(s2): 529-532 [2] Qian Xiangping, Yao Zeen, Wang Qiang. Model-predictive control of power supply for particle accelerators[J]. Nuclear Science and Techniques, 2014, 25(5): 25-28. [3] 蒋舸扬, 谭松清, 赵欢, 等. 上海质子治疗装置动态电源控制系统[J]. 核技术, 2018, 41(2):42-46. (Jiang Geyang, Tan Songqing, Zho Huan, et al. Dynamic power controls in Shanghai Proton Therapy Facility[J]. Nuclear Techniques, 2018, 41(2): 42-46 [4] Ke Xinhua, Lu Songlin, Li Rui, et al. A novel testing approach for SSRF digital power supply controllers[J]. Nuclear Science and Techniques, 2008(4): 209-212. [5] 缪亚运, 谷鸣, 陈志豪, 等. 质子治疗装置脉冲电源研制[J]. 核技术, 2016, 39(4):32-36. (Miao Yanyun, Gu Ming, Chen Zhihao, et al. Development of pulsed power supply in proton therapy[J]. Nuclear Technology, 2016, 39(4): 32-36 [6] Tan S, Li R, Guo C, et al. A novel input power control strategy for high-power dynamic dipole power supply for proton therapy[J]. Nuclear Instruments and Methods in Physics Research Section A: Accelerators, Spectrometers, Detectors and Associated Equipment, 2018, 911: 25-29. [7] 丁建国, 赵欢, 朱海君, 等. 上海同步辐射装置电源控制系统[J]. 核技术, 2012, 35(8):573-577. (Ding Jianguo, Zhao Huan, Zhu Haijun, et al. The SSRF power supply control system[J]. Nuclear Technology, 2012, 35(8): 573-577 [8] 赵江. 重离子加速器数字电源实时控制方法的研究与实现[D]. 兰州: 中国科学院研究生院(近代物理研究所), 2014: 14-31.Zhao Jiang. Research and Implementation of real-time control method for heavy-ion accelerator digital power supply[D]. Lanzhou: Graduate school of Chinese Academy of Sciences (Institute of Modern Physics), 2014: 14-31 [9] 焦喜香, 敬岚, 顾可伟, 等. 重离子加速器注入器电源控制器的设计[J]. 强激光与粒子束, 2010, 22(9):2155-2159. (Jiao Xixiang, Jing Lan, Gu Kewei, et al. Design of digital power control system applied in sector focusing cyclotron[J]. High Power Laser and Particle Beams, 2010, 22(9): 2155-2159 doi: 10.3788/HPLPB20102209.2155 [10] 王永鹏, 郭玉辉, 罗冰峰, 等. 基于FPGA技术的加速器切束控制系统设计[J]. 强激光与粒子束, 2016, 28:105103. (Wang Yongpeng, Guo Yuhui, Luo Bingfeng, et al. Design of accelerator beam cut control system based on FPGA[J]. High Power Laser and Particle Beams, 2016, 28: 105103 doi: 10.11884/HPLPB201628.160004 [11] 李瑞, 谭松清, 郭春龙. 数字化超低纹波可编程动态电源[J]. 电力电子技术, 2014, 48(12):22-24. (Li Rui, Tan Songqing, Guo Chunlong. Digital programmable dynamic power supply with low ripple[J]. Power Electronics, 2014, 48(12): 22-24 doi: 10.3969/j.issn.1000-100X.2014.12.006 [12] Molloy D. Exploring Beaglebone: Tools and techniques for building with embedded Linux[M]. 2nd ed. Indiana: John Wiley & Sons Inc, 2019 [13] Li Fenghan, Yong Zhongchen, Jun Cai, et al. The application of EPICS in TMSR radiation protection and access control system[J]. Nuclear Science and Techniques, 2016, 27(02): 60-65. [14] 张根灿, 蒋舸扬, 李德明, 等. 基于EPICS的质子注入器远程控制及监测系统的构建[J]. 核技术, 2018, 41(03):53-57. (Zhang Gencan, Jiang Geyang, Li Deming, et al. Construct of remote control and monitor system for proton injector based on EPICS[J]. Nuclear Technology, 2018, 41(03): 53-57 [15] 蔡袁琦, 唐雷雷, 周泽然. 基于嵌入式EPICS的合肥光源储存环束流损失监测系统[J]. 强激光与粒子束, 2019, 31:085103. (Cai Yuanqi, Tang Leilei, Zhou Zeran. Embedded EPICS based beam loss monitor system of HLS-Ⅱ storage ring[J]. High Power Laser and Particle Beams, 2019, 31: 085103 doi: 10.11884/HPLPB201931.190041 [16] 朱黎黎, 黄连生, 傅鹏, 等. 基于实验物理和工业控制系统的电源操作节点开发[J]. 强激光与粒子束, 2017, 29:046002. (Zhu Lili, Huang Liansheng, Fu Peng, et al. Development of console system for EAST poloidal field power supply[J]. High Power Laser and Particle Beams, 2017, 29: 046002 doi: 10.11884/HPLPB201729.160438 [17] 刁彦华, 贾宝青, 王晓君. GPMC总线在AM3354与FPGA接口应用中的Android驱动实现[J]. 电子器件, 2017, 40(6):1502-1505. (Diao Yanhua, Jia Baoqing, Wang Xiaojun. Implementation of Android driver in AM3354 and FPGA interface application based on GPMC bus[J]. Chinese Journal of Electron Devices, 2017, 40(6): 1502-1505 doi: 10.3969/j.issn.1005-9490.2017.06.033 [18] Cleva S, Pivetta L, Sigalotti P. Beaglebone for embedded control system applications[C]//Proc ICALEPCS2013. 2013. 期刊类型引用(1)
1. 万臻博，丁卫东，孙凤举. 直线变压器驱动源磁芯磁化过程的仿真研究. 高电压技术. 2024(04): 1792-1800 . 
百度学术其他类型引用(2)
-
下载:
下载:
点击查看大图
计量
- 文章访问数: 1233
- HTML全文浏览量: 339
- PDF下载量: 72
- 被引次数: 3
