ISSN 
下载:
-
重离子束治癌是生物医学与核技术深度交叉和融合的产物。由于重离子具有倒转的深度剂量分布曲线(Bragg曲线)和较高的相对生物效应，重离子束治癌技术可以克服常规放射疗法的局限性[1-2]。在重离子治癌装置[3-5]中，12C束流轰击靶组织(肿瘤)会引发弹核碎裂产生β+发射体(11C、15O、10C)，因此在径迹末端遇负电子会产生湮灭反应，同时释放两个在同一直线上沿相反方向发射的511 keV γ光子。在患者受辐射后立即对该反应空间进行测量，可以测得有效的照射剂量信息[6]。重离子治癌装置要求采用符合探测技术进行γ-γ符合测量，能够从不同探测器的输出脉冲中挑选相关的时间脉冲，进而确认关联事件，甄别排除无效事件，提高探测器的探测效率。符合测量系统通过一对相对放置的探测器测量两个511 keV γ光子击中探测器的时间信息和能量信息，并基于时间窗和能量窗的判选，进而获取符合事件和精确推断出湮灭位置[7]。能量符合的目的是为了消除正电子发射断层扫描成像(positron emission tomography, PET)系统中由康普顿散射造成的散射符合事件，时间符合的目的为了消除系统中的随机符合事件。因此，设计用于重离子治癌装置的符合测量系统必须满足高分辨率和高精度时间测量。
时间量在核物理与粒子物理实验中具有重要意义，时间−数字转换(time digital conversion, TDC)是时间间隔测量的基本方法之一[8-9]，广泛应用于现代物理测量系统，可采用模拟方法、数字方法以及专用的集成电路实现[10]。采用专用集成电路(ASIC)获取高分辨率时间间隔测量，如GP2、HPTDC等专用TDC ASIC芯片[11]，虽然简化了设计复杂度，但是其开发周期长、通用性差。采用现场可编程逻辑阵列(FPGA)技术，由于其编程灵活性强、精度高、可靠性高、通道多特性，同时拥有丰富的延迟链单元、时钟、触发器等资源，使其成为实现数字TDC的主要手段[12]。
传统的符合测量电路通常采用NIM(nuclear instrument module)机箱、线性放大插件、符合插件等，存在着设计复杂，抗干扰性、测量再现、易操作性和系统升级等方面的不足，探索开发新型的符合测量系统迫在眉睫[13]。本文设计基于Altera公司Cyclone Ⅴ系列FPGA实现重离子治癌装置中γ-γ符合时间测量系统，实现高分辨率及高精度的时间间隔测量，本征时间测量精度达百皮秒量级，符合时间谱分辨达1.5 ns (FWHM)。
-
研制的γ-γ符合时间测量系统结构如图1所示，包括晶体阵列探测器、前端读出电路和FPGA电路和PC。前端读出电路包括时间链、能量链、数模转换器(DAC)和时钟模块组成，其中时间链负责处理晶体阵列探测器打拿极输出的时间信号，产生定时脉冲信号；能量链完成电荷信号的滤波成形和模数转换；时钟模块提供FPGA所需的系统时钟。FPGA设计实现了控制模块、TDC模块、积分面积算法模块、数据打包模块及通信模块(USB)，其中控制模块用于配置修正表(Cali_RAM)及控制DAC的阈值电平。TDC模块标记信号的时间戳信息，即通过记录每个击中信号发生的时刻，利用修正表对记录的时间信息进行修正，得到精确的“细”时间，由组帧模块将“细”时间与“粗”时间测量结果合成时间戳。组帧模块根据帧格式将时间戳(T)与能量(E)数据打包暂存于先进先出缓冲器(FIFO)中，而后通过USB[14]进行数据传输和配置。在PC上利用MATLAB编写分析代码读取上传的数据包，顺序提取各个通道的时间和能量信息，基于设置的时间窗筛选符合的事件同时做差处理，得到时间谱。
该系统电路结构简单，具有较强的灵活性和可扩展性，易于系统升级。
图 1 γ-γ符合时间测量系统结构框图
-
符合时间电路中时间链和FPGA-TDC逻辑框图如图2所示。光电倍增管打拿极输出信号为固定的脉冲波形，其宽度、幅度、周期及边沿特性一致性好，适合选用前沿定时甄别。其中时间链采用LEMO接头，阻抗50 Ω，交流耦合至由超高速比较器构成的定时甄别电路，同向端是DAC产生的阈值电平输入，经甄别给出触发信号(start/stop)。由于比较器是PECL电平标准输出，选用的FPGA的IO端口为LVDS电平标准，故电路采用Micrel公司的高速电平转换芯片SY55855将PECL转为LVDS电平送给FPGA，提高信号的抗干扰性。
图 2 符合时间电路时间链和FPGA-TDC模块
FPGA-TDC设计采用“粗细”时间测量相结合的方法实现。“粗”时间测量采用计数器型TDC构建，分辨率为系统时钟f2的周期，实际工作中为了消除亚稳态现象，计数器系统时钟f2与f1有180°的相位差。“细”时间测量单元分为3个部分：缓冲器和寄存器阵列构成的级联延迟链、编码逻辑和Cali_RAM修正模块。在触发电平信号通道中插入一条延迟链，信号在延迟链上传播，在其到来后的第一个时钟上升沿进行锁存寄存器阵列上的电平信号，将输出的多bit电平称为温度码。在编码逻辑中对温度码进行二进制编码，输出“细”时间结果。考虑到由于器件工艺、布局布线、温度以及硬件电路等因素影响，延迟链单元存在非线性，每个单元延迟时间不一致，从而降低“细”时间测量精度。故在本设计中，通过将修正查找表载入Cali_RAM来修正非线性以提高测量精度，得到最后的“细”时间，最终输出“细”和“粗”的时间测量结果，并将时间测量数据按数据帧格式打包暂存于先进先出缓冲器FIFO中；之后通过USB数据通信模块将数据传输至PC机完成后续的处理。
-
时间间隔测量时序如图3所示，时间间隔
$\Delta t$ 是stop与start信号的时间戳差值。时间戳信息${t_{{\rm{start}}}}$ 和${t_{{\rm{stop}}}}$ 由“粗细”时间测量单元实现，如图4所示，外部使能信号连接至二进制计数器使能控制端En，输入固定频率f2(周期T=1/f2)的系统时钟信号，在系统复位后，使能信号En控制计数器启动计数。击中信号(start/stop)分别经过缓冲器和寄存器阵列构成的级联延迟链完成“细”时间的测量，“细”时间测量的分辨为一个延迟单元的延时。Start信号进入延迟链并在其中传输，在start到来之后的第一个系统时钟f1信号的上升沿，寄存器阵列对延迟单元的输出电平进行锁存。随后寄存器锁存的输出作为“细”时间码进入编码模块，采用折半查找算法[12]进行自然二进制数编码，至此“细”时间测量获得$\Delta {t_1}$ ，并输出Selector信号(图2中标记)记录二进制计数器的输出值${N_{{\rm{start}}}}$ ；同理，当stop到来后，获得$\Delta {t_2}$ 和二进制计数器的输出值${N_{{\rm{stop}}}}$ 。根据文献[12]得到“粗”时间测量公式为：
图 3 时间间隔测量时序图
图 4 粗时间测量和细时间测量逻辑框图
$${t_{{\rm{start}}}} = {N_{{\rm{start}}}}T\quad\quad {t_{{\rm{stop}}}} = {N_{{\rm{stop}}}}T$$ (1) 同时测量的时间间隔为：
$$\Delta t = ({t_{{\rm{stop}}}} + \Delta {t_2}) - ({t_{{\rm{start}}}} + \Delta {t_1})$$ (2) -
为减少非线性的影响，本文设计Cali_RAM存储码密度查找表进行“细”时间码的修正。实验中start或stop的频率与系统时钟频率呈弱相关，累计大量数据，将有输出的延迟单元个数作统计分布图，即码密度图，如图5所示。
图 5 码密度示意图
在呈弱相关且大量数据前提下，认为信号的前沿到来时刻均匀分布在一个时钟周期内部，故有输出的延迟单元的个数的统计差异性取决于自身的延迟时间大小，单元的延时计算为：
$${\rm{dela}}{{\rm{y}}_i} = T{S_i}\Big/\sum {{S_i}} $$ (3) 式中，
${S_i}$ 为延迟单元的输出；$\displaystyle\sum {{S_i}} $ 为延迟单元的输出的统计数；T为时钟周期。求得每个延迟单元的延迟时间后，当延迟链有n个单元，即有n个输出时，对应的延迟时间为：$$t = {\rm{dela}}{{\rm{y}}_1} + {\rm{dela}}{{\rm{y}}_2} + \cdots + {\rm{dela}}{{\rm{y}}_n}$$ (4) 把有效输出的延迟单元对应的总延迟时间计算出来后，将其转化为二进制并写入到Cali_RAM中。正常工作状态下得到未修正的“细”计数编码值作为Cali_RAM的地址查找到对应的真实延时值，即可完成时间测量的修正。
-
符合时间测量系统采用Altera高性能Cyclone V 系列。底层架构由4输入查找表(LUTs)构成的自适应逻辑模块(ALM)、存储器模块以及乘法器组成。图6a为实现延迟链底层布线图，每一级均由ALM左端的输入引脚到右端的寄存器输出。图6b为10个ALM组成的1个LAB单元。ALM是Cyclone V中最小的逻辑单元，每一个ALM中包括6个LUT，4个可编程寄存器。
图 6 延迟链底层布局布线
-
通过基于电缆延迟方法的延迟时间测量来研究定时测量的精度[15]。电缆延迟方法旨在消除信号发生器本身的时间抖动：输入大量延迟时间固定的成对快脉冲信号作为触发信号，测量其时间间隔，将测试得到的数据进行统计分析，得到固有时间分辨。利用Tektronix AFG3252双通道信号发生器模拟探测器输出信号，在同样的输出脉冲频率和恒定的室温下，采用不同长度的同轴电缆，选择固定延迟时间为10 ns进行时间间隔的测量，经过统计分析使用高斯拟合[16]方法得到FPGA-TDC的时间分辨。固定的线缆延迟时间下测得时间分布，分辨为255 ps (FWHM)，如图7所示。输入不同延迟时间测得FPGA-TDC的线性及分辨好于276 ps (FWHM)，如图8所示。
图 7 固定延迟时间的分布图
图 8 不同延迟间的线性拟合及均方根值
-
γ-γ符合时间测量系统联合测试实验装置如图9所示。每个探头采用H8500耦合LYSO晶体阵列组成，LYSO阵列用白色的特氟龙胶带和黑色胶带避光包裹起来，并连接到H8500C光电倍增管上。在每个光电倍增管之后，采用DPC桥式电路通过LEMO电缆连接到符合时间测量系统板上[17-18]。 USB电缆将数据传送到PC进行后期数据分析处理。其中的探测器晶体阵列采用由22×22个2 mm×2 mm×15 mm LYSO 晶体条组成，在950 V的正高压偏置下，H8500的增益高达106。光电倍增管打拿极输出信号作为时间信号，用于符合测量。测得有效符合事件时间分布，系统时间测量精度优于1.12 ns (FWHM)，如图10所示。统计分析有效符合事件得到的LYSO晶体阵列位置映射散点图清晰，如图11所示。
图 9 γ-γ符合时间测量系统装置图
图 10 符合时间间隔的时间谱
图 11 LYSO阵列晶体位置映射图
-
本文设计并实现了一种用于重离子治癌装置的γ-γ符合时间测量系统，该系统具有结构简单、易实现、低成本的特点，可方便实现通道数的扩展，且各个通道间独立。经过试验测试，本征时间分辨好于276 ps (FWHM)，具有良好线性。LYSO晶体阵列构成的符合测量系统对511 keV γ射线的测试，时间测量精度优于1.12 ns (FWHM)，略好于文献[19]提出的多相时钟测得符合时间精度。而且晶体阵列位置映射散点图清晰，能够满足γ-γ符合时间测量系统的研究和测试需求。
Design of γ-γ Coincident Time Measurement System in Heavy-Ion Cancer Therapy Device
-
摘要: 为实现重离子治癌装置中γ-γ符合时间的高精度和高分辨率测量，设计实现了符合时间的测量系统，该系统主要由高速比较器构成的定时甄别电路和基于现场可编程逻辑阵列(FPGA)的时间−数字转换(TDC)电路构成。FPGA-TDC通过“粗”时间和“细”时间结合的方法实现精确的时间测量：“粗”时间测量采用二进制计数器实现，“细”时间测量基于片内缓冲器和寄存器阵列构成的级联延迟链实现。实测结果表明，该符合时间测量系统的本征时间间隔测量分辨好于276 ps(FWHM)。构建了H8500耦合LYSO晶体阵列组成的测量系统，并对22Na 511 keV γ射线进行测试，结果表明时间测量精度优于1.12 ns(FWHM)，且通过统计分析有效符合时间得到的晶体阵列位置映射散点图清晰。
-
关键词:
- 符合时间测量系统 /
- 现场可编程逻辑阵列FPGA /
- 时间间隔测量 /
- 时间-数字转换器TDC
Abstract: In order to achieve high-precision and high-resolution measurement of γ-γ coincidence time in the heavy-ion cancer therapy device, a newly implemented coincident time measurement system is designed. The system is mainly composed of a timing discrimination circuit constructed by the high-speed comparator and FPGA-based time-to-digital converter (TDC). A precision time measurement is performed by TDC through the combination of coarse time and fine time measurement methods. The coarse time measurement is realized depending on a binary counter, and the fine time measurement is realized by using the cascade delay chain developed with the in-chip buffer unit and register array. The measured results show that the intrinsic time interval measurement resolution is better than 276 ps (FWHM). The system were tested combining with a coincidence detecting system built by H8500 coupled to LYSO scintillators arrays and the pairs of γ-rays with 511 keV from 22Na source. The results show that the time accuracy of coincidence detecting system is better than 1.12 ns (FWHM) and the black scattered points on LYSO crystal flood image are clear and easy to be distinguished. -
[1] SCHARDT D, ELSASSER T, SCHULZ-ERTNER D. Heavy-ion tumor therapy: Physical and radiobiological benefits[J]. Reviews of Modern Physics, 2010, 82(1): 383-425. doi: 10.1103/RevModPhys.82.383 [2] 王晓林, 高天欣, 韩潇, 等. 重离子放射治疗技术及临床应用[J]. 北京生物医学工程, 2019(3): 312-318. doi: 10.3969/j.issn.1002-3208.2019.03.015. WANG Xiao-ling, GAO Tian-xin, HAN Xiao, et al. Technology and clinical application of heavy ion radiation therapy[J]. Beijing Biomedical Engineering, 2019(3): 312-318. doi: 10.3969/j.issn.1002-3208.2019.03.015. [3] 夏佳文, 詹文龙, 魏宝文, 等. 兰州重离子加速器研究装置HIRFL[J]. 科学通报, 2016, 61(4-5): 467-477. XIA Jia-wen, ZHAN Wen-long, WEI Bao-wen, et al. Heavy ions research facility in lanzhou[J]. Chinese Science Bulletin, 2016, 61(4-5): 467-477. [4] 李斌, 李思琪. 兰州重离子加速器经济社会效益调研[J]. 工程研究−跨学科视野中的工程, 2015, 7(1): 3-15. LI Bin, LI Si-qi. The investigation report on the economic and social benefits of heavy ion research facility in lanzhou (HIRFL)[J]. Journal of Engineering Studies, 2015, 7(1): 3-15. [5] 肖国青, 张红, 李强, 等. 中国科学院近代物理研究所重离子束治癌进展[J]. 原子核物理评论, 2007, 24(2): 85-88. doi: 10.3969/j.issn.1007-4627.2007.02.001 XIAO Guo-qing, ZHANG Hong, LI Qiang, et al. Progresses of heavy-ion cancer therapy in IMP[J]. Nuclear Physics Review, 2007, 24(2): 85-88. doi: 10.3969/j.issn.1007-4627.2007.02.001 [6] ATTANASI F, BELCARI N, CAMARDA M, et al. Preliminary results of an in-beam PET prototype for proton therapy[J]. Nuclear Instruments and Methods in Physics Research Section A: Accelerators, Spectrometers, Detectors and Associated Equipment, 2008, 591(1): 296-299. [7] 杨龙, 梁国栋, 张如美. 基于FPGA的PET快速符合系统设计[J]. 电子设计工程, 2019, 27(8): 180-184. YANG Long, LIANG Guo-dong, ZHANG Ru-mei. FPGA-based rapid coincidence design for PET system[J]. Electronic Design Engineering, 2019, 27(8): 180-184. [8] 叶春逢. 飞行时间质谱仪数据获取系统的研究与设计[D]. 合肥: 中国科学技术大学, 2014. YE Chun-feng. Research and design of time-of-flight mass spectrometer data acquisition system[D]. Hefei: University of Science and Technology of China, 2014. [9] 郭建华. 北京谱仪(BESIII)飞行时间读出电子学系统设计与实现[D]. 合肥: 中国科学技术大学, 2007. GUO Jian-hua. Design and implementation of time-of-flight readout electronics for Beijing spectrometer (BESIII)[D]. Hefei: University of Science and Technology of China, 2007. [10] 李广军, 孟宪元. 可编程ASIC设计及应用[M]. 成都: 电子科技大学出版社, 2000. LI Guang-jun, MENG Xian-yuan. Design and application of programmable ASIC[M]. Chengdu: University of Electronic Science and Technology of China Press, 2000. [11] 李良辉, 千奕, 孔洁, 等. 基于TDC的时间数字转换模块研制[J]. 核电子学与探测技术, 2015, 35(12): 1263-1266. doi: 10.3969/j.issn.0258-0934.2015.12.025 LI Liang-hui, QIAN Yi, KONG Jie, et al. Development of time to digital convert module based on TDC-Chip[J]. Nuclear Electronics and Detection Technology, 2015, 35(12): 1263-1266. doi: 10.3969/j.issn.0258-0934.2015.12.025 [12] 宋健. 基于FPGA的精密时间−数字转换电路研究[D]. 合肥: 中国科学技术大学, 2006. SONG Jian. Research on precision time−Digital conversion circuit based on FPGA[D]. Hefei: University of Science and Technology of China, 2006. [13] 颜拥军, 赖伟, 周剑良, 等. 核信号数字符合测量的实验研究[J]. 南华大学学报(自然科学版), 2012, 26(1): 5-9. YAN Yong-jun, LAI Wei, ZHOU Jian-liang, et al. Experimental research of digital coincidence measurement[J]. Journal of University of South China(Science and Technology), 2012, 26(1): 5-9. [14] 陈松. 基于USB3.0的高速数据传输系统的研究与设计[D]. 成都: 电子科技大学, 2014. CHEN Song. Research and design of the high speed data transmission system based on USB 3.0[D]. Chengdu: University of Electronic Science and Technology of China, 2014. [15] ZHAO Lei, KANG Long-fei, ZHOU Jia-wen, et al. A 16-Channel high-resolution time and charge measurement module for the external target experiment in the CSR of HIRFL[J]. Nuclear Techniques, 2014, 25(1): 27-32. [16] JEFFERYS W H, FITZPATRICK M J, MCARTHUR B E. GaussFit—A system for least squares and robust estimation[J]. Celestial Mechanics, 1988, 41(1-4): 39-49. [17] 陈金达. 重离子治癌在线影像实时监测技术的研究—双头平面型在束TOF-PET系统[D]. 兰州: 中国科学院近代物理研究所, 2011. CHEN Jin-da. Research on online image real-time monitoring technology for heavy ion treatment of Cancer—Two-head planar Beam TOF-PET system[D]. Lanzhou: Institute of Modern Physics, Chinese Academy of Sciences, 2011. [18] 杜成名. 数字化波形采样技术在粒子探测中的应用与研究[D]兰州: 中国科学院近代物理研究所, 2016. DU Cheng-ming. Application and research of digital waveform sampling technology in particle detection[D]. Lanzhou: Institute of Modern Physics, Chinese Academy of Sciences, 2016. [19] DENG P, ZHAO L, LU J, et al. Prototype design of singles processing unit for the small animal PET[J]. Journal of Instrumentation, 2018, 13(5): T05007. doi: 10.1088/1748-0221/13/05/T05007 -
点击查看大图
图(11)
计量
- 文章访问数: 4673
- HTML全文浏览量: 1459
- PDF下载量: 25
- 被引次数: 0
