CN115758735A - A scanning magnet dynamic identification real-time segmented slope feedback method and system - Google Patents
A scanning magnet dynamic identification real-time segmented slope feedback method and system Download PDFInfo
- Publication number
- CN115758735A CN115758735A CN202211454473.XA CN202211454473A CN115758735A CN 115758735 A CN115758735 A CN 115758735A CN 202211454473 A CN202211454473 A CN 202211454473A CN 115758735 A CN115758735 A CN 115758735A
- Authority
- CN
- China
- Prior art keywords
- slope
- power supply
- triangular wave
- zero
- positive
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 41
- 238000009826 distribution Methods 0.000 claims abstract description 40
- 230000005284 excitation Effects 0.000 claims abstract description 21
- 230000011218 segmentation Effects 0.000 claims abstract description 13
- XEEYBQQBJWHFJM-UHFFFAOYSA-N Iron Chemical compound [Fe] XEEYBQQBJWHFJM-UHFFFAOYSA-N 0.000 claims description 43
- 238000005457 optimization Methods 0.000 claims description 34
- 238000012937 correction Methods 0.000 claims description 22
- 229910052742 iron Inorganic materials 0.000 claims description 21
- 238000004590 computer program Methods 0.000 claims description 11
- 238000003860 storage Methods 0.000 claims description 6
- 238000010586 diagram Methods 0.000 description 28
- 230000000694 effects Effects 0.000 description 6
- 230000005540 biological transmission Effects 0.000 description 4
- 238000004422 calculation algorithm Methods 0.000 description 4
- 238000002474 experimental method Methods 0.000 description 4
- 238000012545 processing Methods 0.000 description 4
- 238000011160 research Methods 0.000 description 4
- 206010028980 Neoplasm Diseases 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 3
- 201000011510 cancer Diseases 0.000 description 3
- 230000006870 function Effects 0.000 description 3
- 230000001133 acceleration Effects 0.000 description 2
- 230000003247 decreasing effect Effects 0.000 description 2
- 239000002245 particle Substances 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 238000001959 radiotherapy Methods 0.000 description 2
- 230000002457 bidirectional effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 230000005684 electric field Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000010223 real-time analysis Methods 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
Images
Classifications
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E30/00—Energy generation of nuclear origin
- Y02E30/10—Nuclear fusion reactors
Landscapes
- Radiation-Therapy Devices (AREA)
Abstract
Description
技术领域technical field
本发明涉及一种扫描磁铁动态识辨实时分段斜率反馈方法和系统,属于同步加速器技术领域。The invention relates to a scanning magnet dynamic identification real-time segmented slope feedback method and system, belonging to the technical field of synchrotrons.
背景技术Background technique
如今,加速器在基础研究和拓展应用方面都有了十分迅猛的发展。加速器有非常多的种类:直线加速器、回旋加速器、同步加速器等。各种加速器将粒子加速到需要的能量之后,通常会选择通过输运线传输到束流终端进行使用。为了保持较大的束流传输效率,同时减小加速器元件尺寸以节约成本,传输过程中束流尺寸一般保持在较小的水平。如果最终需要的束流尺寸远大于传输过程中的尺寸,通常使用扫描磁铁增大束斑以达到终端所需的束流尺寸指标。扫描磁铁是一种应用广泛的磁铁元件,结构和原理类似普通的二极磁铁,通过铁芯和一对励磁线圈产生垂直于线圈平面的二极磁场,对束流产生正比于励磁电流的偏转力。Today, accelerators have developed rapidly in both basic research and expanded applications. There are many types of accelerators: linear accelerators, cyclotrons, synchrotrons, etc. After various accelerators accelerate the particles to the required energy, they are usually selected to be transported to the end of the beam through the transport line for use. In order to maintain a large beam transmission efficiency and reduce the size of accelerator components to save costs, the beam size is generally kept at a small level during transmission. If the final required beam size is much larger than the size in the transmission process, usually a scanning magnet is used to increase the beam spot to achieve the desired beam size index at the end. The scanning magnet is a widely used magnetic component. Its structure and principle are similar to ordinary two-pole magnets. A two-pole magnetic field perpendicular to the coil plane is generated through an iron core and a pair of excitation coils, and a deflection force proportional to the excitation current is generated on the beam current. .
由于扫描磁铁电源电流需要持续双向且线性可调,输出的电流波形往往存在过零畸变的问题,也就是每次电流在正负之间转换时都会产生一定的跟踪偏差。此外,由于励磁电流工作在高频率的三角波模式下,扫描磁铁的震动、涡流等效应也会造成实际产生的磁场偏离理想值,最终得到一个均匀性较差的大尺寸束斑,降低终端最终获得的束流品质、影响放射治疗、辐照实验等加速器束流应用的效果。Since the current of the scanning magnet power supply needs to be continuously bidirectional and linearly adjustable, the output current waveform often has the problem of zero-crossing distortion, that is, a certain tracking deviation will be generated every time the current switches between positive and negative. In addition, since the excitation current works in the high-frequency triangular wave mode, the vibration and eddy current effects of the scanning magnet will also cause the actual generated magnetic field to deviate from the ideal value, and finally a large-sized beam spot with poor uniformity will be obtained, which will reduce the terminal's final gain. The quality of the beam, affecting the effect of accelerator beam applications such as radiation therapy and irradiation experiments.
发明内容Contents of the invention
针对上述问题,本发明的目的是提供一种扫描磁铁动态识辨实时分段斜率反馈方法和系统,其将束流均匀性提高到90%以上,可以优化高终端均匀性。In view of the above problems, the object of the present invention is to provide a scanning magnet dynamic identification real-time segmented slope feedback method and system, which can improve the beam uniformity to more than 90%, and can optimize high terminal uniformity.
为实现上述目的,本发明提出了以下技术方案:一种扫描磁铁动态识辨实时分段斜率反馈方法,包括以下步骤:获取束流在水平和垂直方向的束流位置分布信息;根据束流位置分布信息获得不同种类的误差信息;根据误差信息进行分段斜率反馈,得到优化后的电源三角波励磁曲线;将电源三角波励磁曲线输入磁铁电源,束流再次经过后获得新的分段斜率反馈,直至终端的扫描束斑满足预设的均匀性。In order to achieve the above object, the present invention proposes the following technical solutions: a real-time segmental slope feedback method for dynamic identification of scanning magnets, comprising the following steps: obtaining beam position distribution information of the beam in the horizontal and vertical directions; Different types of error information are obtained from the distribution information; the segmented slope feedback is performed according to the error information, and the optimized triangular wave excitation curve of the power supply is obtained; the triangular wave excitation curve of the power supply is input into the magnet power supply, and a new segmented slope feedback is obtained after the beam passes through again, until The scanning beam spot of the terminal satisfies the preset uniformity.
进一步,误差信息的种类包括:分段斜率偏差、正负半轴偏差和过零区域凸起。Further, the types of error information include: segmental slope deviation, positive and negative semi-axis deviation, and zero-crossing area bulge.
进一步,分段斜率反馈的方法为:通过电源波形的正负零漂解决正负半轴偏差;通过电源波形的过零截距解决过零区域凸起;将经过正负零漂和过零截距校正后的电源三角波按照斜率进行分段,并根据分段结果对电源三角波进行优化;将经过优化的电源三角波发送至扫描磁铁电源,束流再次经过后,对分段斜率优化的结果进行反馈。Further, the method of segmented slope feedback is: solve the positive and negative semi-axis deviation through the positive and negative zero drift of the power waveform; solve the bulge of the zero crossing area through the zero crossing intercept of the power waveform; The corrected power triangular wave is segmented according to the slope, and the power triangular wave is optimized according to the segmentation results; the optimized power triangular wave is sent to the scanning magnet power supply, and after the beam passes through again, the result of the segmented slope optimization is fed back .
进一步,正负零漂是指扫描铁电流波形的正负幅度不再相等,可以分别增大或减小。Further, positive and negative zero drift means that the positive and negative amplitudes of the scanning ferrocurrent waveform are no longer equal, and can be increased or decreased respectively.
进一步,过零截距是指在扫描铁波形的正负半边的绝对值上,分别整体增加一个固定值。Further, the zero-crossing intercept refers to adding a fixed value as a whole to the absolute value of the positive and negative halves of the scanning iron waveform.
进一步,电源三角波按照斜率进行分段的方法为:将束流位置分布信息中的计数曲线按位置分段,并将电源波形斜率曲线按电流值进行对应范围同样数量的分段,将两个步骤获得的分段曲线的每一段进行对应,获得最终的分段曲线。Further, the method of segmenting the triangular wave of the power supply according to the slope is as follows: the counting curve in the beam position distribution information is segmented according to the position, and the slope curve of the power supply waveform is divided into the same number of segments corresponding to the current value, and the two steps Each segment of the obtained segmented curve is corresponding to obtain the final segmented curve.
进一步,终端束流的计算公式为:Further, the calculation formula of the terminal beam current is:
其中,Amax和Amin分别是分条电离室束流分布数据中最大和最小的计数,U是终端束流均匀性。Among them, A max and A min are the maximum and minimum counts in the beam distribution data of the striped ionization chamber, respectively, and U is the terminal beam uniformity.
本发明还公开了一种扫描磁铁动态识辨实时分段斜率反馈系统,包括:束流位置分布信息获取模块,用于获取束流在水平和垂直方向的束流位置分布信息;误差信息分类模块,用于根据束流位置分布信息获得不同种类的误差信息;误差优化模块,用于根据误差信息进行分段斜率反馈,得到优化后的电源三角波励磁曲线;输出反馈模块,用于将电源三角波励磁曲线输入磁铁电源,束流再次经过后获得新的分段斜率反馈,直至终端的扫描束斑满足预设的均匀性。The invention also discloses a scanning magnet dynamic identification real-time segmental slope feedback system, which includes: a beam position distribution information acquisition module for acquiring beam position distribution information in the horizontal and vertical directions of the beam; an error information classification module , used to obtain different types of error information according to the beam current position distribution information; the error optimization module is used to perform segmental slope feedback according to the error information, and obtain the optimized triangular wave excitation curve of the power supply; the output feedback module is used to excite the triangular wave of the power supply The curve is input to the magnet power supply, and the new segmental slope feedback is obtained after the beam passes through again, until the scanning beam spot at the terminal meets the preset uniformity.
进一步,误差优化模块包括:正负零漂模块,用于通过电源波形的正负零漂解决正负半轴偏差;过零截距模块,用于通过电源波形的过零截距解决过零区域凸起;分段斜率优化模块,用于将经过正负零漂和过零截距校正后的电源三角波按照斜率进行分段,并根据分段结果对电源三角波进行优化;反馈模块,用于将经过优化的电源三角波发送至电源和扫描磁铁,束流再次经过后,对分段斜率优化的结果进行反馈。Further, the error optimization module includes: a positive and negative zero drift module, which is used to solve the positive and negative semi-axis deviation through the positive and negative zero drift of the power waveform; a zero-crossing intercept module, which is used to solve the zero-crossing area through the zero-crossing intercept of the power waveform Convex; the segmented slope optimization module is used to segment the triangular wave of the power supply after positive and negative zero drift and zero-crossing intercept correction according to the slope, and optimize the triangular wave of the power supply according to the segmentation results; the feedback module is used to The optimized triangular wave of the power supply is sent to the power supply and the scanning magnet, and after the beam current passes through again, the result of segmented slope optimization is fed back.
本发明还公开了一种计算机可读存储介质,计算机可读存储介质上存储有计算机程序,计算机程序被处理器执行以实现如上述任一项扫描磁铁动态识辨实时分段斜率反馈方法。The present invention also discloses a computer-readable storage medium, on which a computer program is stored, and the computer program is executed by a processor to realize the real-time segmental slope feedback method for dynamic identification of any scanning magnet described above.
本发明还一种计算设备,包括:一个或多个处理器、存储器以及一个或多个程序,其中一个或多个程序存储在存储器中并被配置为由一个或多个处理器执行,一个或多个程序包括用于执行如上述任一项扫描磁铁动态识辨实时分段斜率反馈方法。The present invention is also a computing device, comprising: one or more processors, memory and one or more programs, wherein one or more programs are stored in the memory and configured to be executed by the one or more processors, one or more A plurality of programs include a real-time segmented slope feedback method for performing any scanning magnet dynamic identification as described above.
本发明由于采取以上技术方案,其具有以下优点:The present invention has the following advantages due to the adoption of the above technical scheme:
1、本发明不管是小型化的治癌加速器还是庞大的同步加速器,只要在终端具备扫描铁和分条电离室类似的束流分布探测器,就可以应用此发明。本发明可以直接检测并优化束流均匀性,具有高度兼容性和可靠性,无需引入额外的硬件设备,成本低廉使用便捷。1. No matter whether the present invention is a miniaturized cancer treatment accelerator or a huge synchrotron, as long as the terminal has a beam distribution detector similar to a scanning iron and a striped ionization chamber, this invention can be applied. The invention can directly detect and optimize the uniformity of the beam current, has high compatibility and reliability, does not need to introduce additional hardware equipment, and is cheap and convenient to use.
2、本发明直接使用分条电离室采集的束流分布数据,实时计算优化波形并发送到扫描磁铁电源,电源更新波形后立即可以再次用分条电离室采集优化后的束流分布数据并再次优化,重复数次即可反复迭代来进行反馈控制。计算优化波形在数秒内即可完成,主要花费时间来源于分条电离室数据采集和传输、扫描铁电源波形更新以及束流本身的脉冲间隔。2. The present invention directly uses the beam distribution data collected by the striped ionization chamber, calculates the optimized waveform in real time and sends it to the scanning magnet power supply, and immediately after the power supply updates the waveform, the optimized beam distribution data can be collected by the striped ionization chamber again and again Optimization, repeating several times can be iterated repeatedly for feedback control. The calculation of the optimized waveform can be completed within a few seconds, and the time spent mainly comes from the data acquisition and transmission of the striped ionization chamber, the update of the scanning iron power supply waveform, and the pulse interval of the beam itself.
3、本发明使用的算法经过几个月的在线运行,能够高效分析终端均匀性的各类误差来源,并分别在扫描磁铁电源波形上予以最有效的校正手段和校正量。扫描磁铁波形与终端均匀性之间有一定的对应关系,但是最终产生误差的来源大不相同,如果直接在波形上采用同样的校正方法往往无法获得良好的校正效果。通过总结大量的束流测试数据,此算法能够简洁高效地优化束流均匀性。3. After several months of online operation, the algorithm used in the present invention can efficiently analyze various error sources of terminal uniformity, and provide the most effective correction means and correction amount on the scanning magnet power waveform respectively. There is a certain correspondence between the waveform of the scanning magnet and the uniformity of the terminal, but the source of the final error is very different. If the same correction method is directly used on the waveform, it is often impossible to obtain a good correction effect. By summarizing a large amount of beam test data, this algorithm can optimize beam uniformity concisely and efficiently.
4、本发明可以更少地考虑磁铁复杂的涡流、震动效应,以及电源的过零点校正等优化,直接从结果上获取最终的束流均匀性,并反馈到源头的电源波形,可以大幅降低中间过程优化的资金和时间成本,廉价高效地优化终端束流均匀性,成为癌症治疗、材料辐照、航空航天研究等领域应用和实验加速器的最佳选择。4. The present invention can take less consideration of the complex eddy current and vibration effects of the magnet, and the zero-crossing correction of the power supply, etc., and directly obtain the final beam current uniformity from the result, and feed back to the source power waveform, which can greatly reduce the intermediate The capital and time cost of process optimization, and the low-cost and efficient optimization of the uniformity of the terminal beam have become the best choice for applications and experimental accelerators in the fields of cancer treatment, material irradiation, and aerospace research.
附图说明Description of drawings
图1是本发明一实施例中分条电离室的束流位置与励磁电流之间的关系图;Fig. 1 is a relation diagram between the beam current position and the excitation current of a striped ionization chamber in an embodiment of the present invention;
图2是本发明一实施例中校正前的分条电离室的束流位置分布信息图;Fig. 2 is an information map of the beam position distribution of the striped ionization chamber before correction in an embodiment of the present invention;
图3是本发明一实施例中电流波形的分段斜率校正图,图3(a)为扫描铁电源输出的三角波形图,图3(b)为扫描铁电源输出的三角波的斜率图;Fig. 3 is the segmented slope correction diagram of current waveform in an embodiment of the present invention, Fig. 3 (a) is the triangular waveform diagram that scanning iron power supply outputs, Fig. 3 (b) is the slope diagram of the triangular wave output of scanning iron power supply;
图4是本发明一实施例中分段斜率反馈的方法的流程图;Fig. 4 is a flowchart of a method for segmented slope feedback in an embodiment of the present invention;
图5是本发明一实施例中分段斜率反馈的方法的示意图,图5(a)为分条电离室采集的束流位置分布信息,图5(b)为扫描铁电源输出的三角波形图,图5(c)为过零区域凸起的示意图,图5(d)为过零点截距优化的示意图,图5(e)为正负半轴偏差的示意图,图5(f)为正负零漂优化的示意图,图5(g)为分段计数曲线的示意图,图5(h)为分段斜率优化的示意图,图5(i)为校正完成后的曲线图;Fig. 5 is a schematic diagram of the method of segmented slope feedback in an embodiment of the present invention, Fig. 5 (a) is the beam position distribution information collected by the striped ionization chamber, and Fig. 5 (b) is a triangular waveform diagram output by the scanning iron power supply , Figure 5(c) is a schematic diagram of the bulge in the zero-crossing area, Figure 5(d) is a schematic diagram of the optimization of the zero-crossing intercept, Figure 5(e) is a schematic diagram of the positive and negative semi-axis deviations, Figure 5(f) is the positive A schematic diagram of negative zero drift optimization, Figure 5(g) is a schematic diagram of a segmented counting curve, Figure 5(h) is a schematic diagram of a segmented slope optimization, and Figure 5(i) is a graph after calibration;
图6是本发明一实施例中校正后的分条电离室的束流位置分布信息图。Fig. 6 is an information map of the beam position distribution of the striped ionization chamber after correction in an embodiment of the present invention.
具体实施方式Detailed ways
为了使本领域技术人员更好的理解本发明的技术方案,通过具体实施例对本发明进行详细的描绘。然而应当理解,具体实施方式的提供仅为了更好地理解本发明,它们不应该理解成对本发明的限制。在本发明的描述中,需要理解的是,所用到的术语仅仅是用于描述的目的,而不能理解为指示或暗示相对重要性。In order to enable those skilled in the art to better understand the technical solution of the present invention, the present invention is described in detail through specific examples. However, it should be understood that specific embodiments are provided only for better understanding of the present invention, and they should not be construed as limiting the present invention. In describing the present invention, it should be understood that the terms used are for the purpose of description only, and should not be understood as indicating or implying relative importance.
为了解决现有技术中存在的扫描磁铁电源电流输出的电流波形往往存在过零畸变的问题,也就是每次电流在正负之间转换时都会产生一定的跟踪偏差;以及由于励磁电流工作在高频率的三角波模式下,扫描磁铁的震动、涡流等效应造成实际产生的磁场偏离理想值,最终使得束流终端得到一个均匀性较差的大尺寸束斑,降低终端最终获得的束流品质,影响放射治疗、辐照实验等加速器束流应用的效果等问题。本发明提出了一种扫描磁铁动态识辨实时分段斜率反馈方法、系统、可读介质和计算设备,实时分析分条电离室获得的束流分布,可以得到上述不同种类的误差信息,称之为“动态识辨”。分条电离室获得的计数曲线通过以分段斜率优化为主的优化算法后,加载优化波形到扫描铁电源,并通过下一次分条电离室的测量结果得到新的优化波形,数次迭代反馈优化后可以得到足够高均匀性的扫描均匀性,称之为“实时分段斜率反馈”。该动态识辨实时分段斜率反馈系统,通过分析终端分条电离室的束流分布,可以实时反馈优化扫描铁励磁电流,将终端束流均匀性提升到90%以上,满足研究、实验和应用的需求。下面结合附图通过实施例对本发明的方案进行详细阐述。In order to solve the problem that the current waveform of the scanning magnet power supply current output in the prior art often has zero-crossing distortion, that is, a certain tracking deviation will be generated every time the current is switched between positive and negative; and because the excitation current works at high In the triangular wave mode of the frequency, the vibration and eddy current effects of the scanning magnet cause the actual magnetic field to deviate from the ideal value, and finally make the beam terminal obtain a large-sized beam spot with poor uniformity, which reduces the final beam quality obtained by the terminal and affects Effects of accelerator beam applications such as radiation therapy and irradiation experiments. The present invention proposes a scanning magnet dynamic identification real-time segmental slope feedback method, system, readable medium and computing equipment, real-time analysis of the beam distribution obtained by the striped ionization chamber, and the above-mentioned different types of error information can be obtained, which is called For "dynamic identification". After the counting curve obtained by the striped ionization chamber passes through the optimization algorithm based on segmented slope optimization, the optimized waveform is loaded to the scanning iron power supply, and a new optimized waveform is obtained through the measurement results of the next striped ionization chamber, and several iterations are fed back After optimization, scan uniformity with high enough uniformity can be obtained, which is called "real-time segmented slope feedback". This dynamic identification real-time segmented slope feedback system, by analyzing the beam current distribution of the terminal strip ionization chamber, can feed back and optimize the excitation current of the scanning iron in real time, and improve the uniformity of the terminal beam current to more than 90%, which meets the needs of research, experiment and application demand. The solutions of the present invention will be described in detail below through embodiments in conjunction with the accompanying drawings.
实施例一Embodiment one
本实施例公开了一种扫描磁铁动态识辨实时分段斜率反馈方法,包括以下步骤:This embodiment discloses a scanning magnet dynamic identification real-time segmented slope feedback method, including the following steps:
S1获取束流在水平和垂直方向的束流位置分布信息;S1 obtains the beam position distribution information of the beam in the horizontal and vertical directions;
本实施例采用分条电离室获取束流在水平和垂直方向的束流位置分布信息,如果其他束诊设备能起到相同作用也可以进行替代。分条电离室是一种常用的束诊设备,其中有一对电极板并充有稀薄的气体,束流经过气体时电离出正负离子对,电荷沿着垂直于束流运动方向的电场运动,被分割成一定数量条状的极板接收,在对应的极板上产生电信号。不同的极板对应不同的空间位置,水平和垂直方向的条状极板组合在一起,通过分析电信号可以得到束流的位置分布信息,电流信号的强弱正比于该处粒子的数量。对于束流终端或治癌加速器,分条电离室通常在验收时放在靶室或病床附近,得到束流实际应用时的束流位置分布信息,分条电离室的束流位置与励磁电流之间的关系如图1所示。In this embodiment, the striped ionization chamber is used to obtain the beam position distribution information of the beam in the horizontal and vertical directions, and other beam diagnosis equipment can be substituted if it can play the same role. The striped ionization chamber is a commonly used beam diagnostic equipment, which has a pair of electrode plates and is filled with thin gas. When the beam passes through the gas, positive and negative ion pairs are ionized, and the charges move along the electric field perpendicular to the direction of the beam. The polar plates divided into a certain number of strips receive and generate electrical signals on the corresponding polar plates. Different plates correspond to different spatial positions. The horizontal and vertical strip plates are combined, and the position distribution information of the beam can be obtained by analyzing the electrical signal. The strength of the current signal is proportional to the number of particles at the location. For beam terminals or cancer treatment accelerators, the striped ionization chamber is usually placed near the target room or the hospital bed during acceptance, and the beam position distribution information in the actual application of the beam is obtained. The relationship between the beam position of the striped ionization chamber and the excitation current The relationship between them is shown in Figure 1.
S2根据束流位置分布信息获得不同种类的误差信息;S2 obtains different types of error information according to the beam position distribution information;
由于图1中的扫描磁铁加载三角波,理想情况下束流在分条电离室上的位置也随时间均匀变化,分条电离室会显示平整的束流分布。但由于电源、磁铁和束流本身各方面的影响,分条电离室实际显示的束流分布如图2所示,束流并非均匀分布,而本发明的目的就是校正这种不均匀性。本发明首先提出,通过调节扫描铁波形斜率来校正束流被扫描铁展开时的不均匀性。Since the scanning magnet in Figure 1 is loaded with a triangular wave, ideally the position of the beam on the striped ionization chamber also changes uniformly with time, and the striped ionization chamber will display a flat beam distribution. However, due to the influence of power supply, magnet and beam current itself, the beam current distribution actually displayed by the striped ionization chamber is shown in Figure 2, and the beam current is not evenly distributed, and the purpose of the present invention is to correct this inhomogeneity. The present invention first proposes to correct the inhomogeneity when the beam is spread by the scanning iron by adjusting the slope of the scanning iron waveform.
图3是本发明一实施例中电流波形的分段斜率校正图,图3(a)为扫描铁电源输出的三角波形图,图3(b)为扫描铁电源输出的三角波的斜率图;可以看出在未做校正的情况下,图3(b)中虚线代表原波形斜率为正负变化的恒定值。此时将波形分段以便于调节,每一小段波形的斜率都可以单独调节,图3(b)中的实线就是根据束流分布反馈调节的分段斜率曲线,原本波形与调节后的校正量在图3(a)中分别用实线和虚线表示。Fig. 3 is the subsection slope correction figure of current waveform in an embodiment of the present invention, and Fig. 3 (a) is the triangular waveform figure that scanning iron power supply outputs, and Fig. 3 (b) is the slope figure of the triangular wave output of scanning iron power supply; It can be seen that in the case of no correction, the dotted line in Figure 3(b) represents the constant value of the slope of the original waveform as a positive or negative change. At this time, the waveform is segmented for easy adjustment, and the slope of each small segment of the waveform can be adjusted individually. The solid line in Figure 3(b) is the segmented slope curve adjusted according to the beam current distribution feedback. The original waveform and the adjusted correction Quantities are represented by solid and dashed lines in Fig. 3(a), respectively.
S3分条电流室接收束流信号后,根据误差信息进行分段斜率反馈,得到优化后的电源三角波励磁曲线;After receiving the beam current signal, the S3 sub-strip current chamber performs segmented slope feedback according to the error information, and obtains the optimized triangular wave excitation curve of the power supply;
误差信息的种类包括:分段斜率偏差、正负半轴偏差和过零区域凸起。如图4所示,分段斜率反馈的方法为:The types of error information include: segmental slope deviation, positive and negative semi-axis deviation, and zero-crossing area bulge. As shown in Figure 4, the method of segmented slope feedback is:
S3.1通过电源波形的正负零漂解决正负半轴偏差;S3.1 Solve the positive and negative semi-axis deviation through the positive and negative zero drift of the power waveform;
S3.2通过电源波形的过零截距解决过零区域凸起;S3.2 Solve the zero-crossing region bulge through the zero-crossing intercept of the power waveform;
S3.3将经过正负零漂和过零截距校正后的电源三角波按照斜率进行分段,并根据分段结果对电源三角波进行优化;S3.3 Segment the triangular wave of the power supply after positive and negative zero drift and zero-crossing intercept correction according to the slope, and optimize the triangular wave of the power supply according to the segmentation results;
S3.4将经过优化的电源三角波发送至电源和扫描磁铁,束流再次经过后,对分段斜率优化的结果进行反馈。S3.4 sends the optimized triangular wave of the power supply to the power supply and the scanning magnet, and after the beam current passes through again, the result of segmented slope optimization is fed back.
图5是通过算法分析分条电离室计数曲线并在扫描铁波形上进行校正的示意图。图5(a)为分条电离室采集的束流位置分布信息,图5(b)为扫描铁电源输出的三角波形图,图5(c)为过零区域凸起的示意图,图5(d)为过零点截距优化的示意图,图5(e)为正负半轴偏差的示意图,图5(f)为正负零漂优化的示意图,图5(g)为分段计数曲线的示意图,图5(h)为分段斜率优化的示意图,图5(i)为校正完成后的曲线图;Fig. 5 is a schematic diagram of analyzing the counting curve of the striped ionization chamber by an algorithm and performing correction on the scanning iron waveform. Figure 5(a) is the beam position distribution information collected by the striped ionization chamber, Figure 5(b) is the triangular waveform diagram output by the scanning iron power supply, Figure 5(c) is a schematic diagram of the bulge in the zero-crossing area, and Figure 5( d) is a schematic diagram of zero-crossing intercept optimization, Figure 5(e) is a schematic diagram of positive and negative semi-axis deviations, Figure 5(f) is a schematic diagram of positive and negative zero drift optimization, and Figure 5(g) is a schematic diagram of a segmented counting curve Schematic diagram, Fig. 5 (h) is the schematic diagram of subsection slope optimization, and Fig. 5 (i) is the curve diagram after correction;
扫描铁波形分段斜率控制是主要的优化方法,对应于图5(g)和图5(h)中校正优化过程。电源三角波按照斜率进行分段的方法为:将束流位置分布信息中的计数曲线,即图2中的计数曲线,按位置(横轴)分段后(如图5(g)所示),并将电源波形斜率曲线按三角波电流值(纵轴)对应范围同样进行同样数量的分段(如图5(h)所示),将两个步骤获得的分段曲线的每一段进行对应,获得最终的分段曲线。最终的校正结果如图5(i)所示,整体的经过校正优化的扫描铁波形图如图3所示。但仅通过分段斜率控制不能校正所有的不均匀性,还需要另外两个优化方法,这两种优化方法没有先后之别,但两种优化方法需要在分段斜率控制之前进行。The segmental slope control of the scanning iron waveform is the main optimization method, corresponding to the correction optimization process in Fig. 5(g) and Fig. 5(h). The method of segmenting the triangular wave of the power supply according to the slope is as follows: after the counting curve in the beam position distribution information, that is, the counting curve in Figure 2, is segmented according to the position (horizontal axis) (as shown in Figure 5(g)), And the slope curve of the power supply waveform is also divided into the same number of segments according to the corresponding range of the triangular wave current value (vertical axis) (as shown in Figure 5(h)), and each segment of the segmented curve obtained in the two steps is correspondingly obtained. The final segmented curve. The final correction result is shown in Fig. 5(i), and the overall corrected and optimized scanning iron waveform is shown in Fig. 3 . However, all inhomogeneities cannot be corrected by segmental slope control alone, and two other optimization methods are required. There is no sequence difference between these two optimization methods, but the two optimization methods need to be performed before segmental slope control.
对于图2中过零区域的凸起,在扫描铁波形中用过零点截距来解决。过零截距是指在扫描铁波形的正负半边的绝对值上,可以分别整体增加一个固定值,对应于图5(d)中校正优化。For the bulge in the zero-crossing area in Figure 2, the zero-crossing intercept is used in the scanning iron waveform to solve it. The zero-crossing intercept means that a fixed value can be added to the absolute value of the positive and negative halves of the scanning iron waveform, corresponding to the correction optimization in Figure 5(d).
对于图2中正负半轴一高一低的情况,在扫描铁波形中通过加正负零漂来解决。正负零漂是指扫描铁电流波形的正负幅度不再相等,可以分别增大或减小,对应于图5(f)中校正优化。For the situation in Figure 2 that the positive and negative half axes are one high and one low, it can be solved by adding positive and negative zero drift in the scanning iron waveform. Positive and negative zero drift means that the positive and negative amplitudes of the scanning ferrocurrent waveform are no longer equal, and can be increased or decreased respectively, corresponding to the correction optimization in Figure 5(f).
S4将电源三角波励磁曲线输入磁铁电源,获得新的分段斜率反馈,直至终端束流的均匀性满足预设条件的终端扫描束斑。S4 inputs the triangular wave excitation curve of the power supply to the magnet power supply to obtain a new segmental slope feedback until the uniformity of the terminal beam current meets the preset condition of the terminal scanning beam spot.
该动态识辨实时分段斜率反馈系统,通过分析终端分条电离室的束流分布,可以实时反馈优化扫描铁励磁电流,将终端束流均匀性提升到90%以上,满足研究、实验和应用的需求。This dynamic identification real-time segmented slope feedback system, by analyzing the beam current distribution of the terminal strip ionization chamber, can feed back and optimize the excitation current of the scanning iron in real time, and improve the uniformity of the terminal beam current to more than 90%, which meets the needs of research, experiment and application demand.
如图6所示,校正前后的分条电离室的束流分布数据,束流均匀性从76.6%提升到93.6%。图中可以看出过零点畸变、正负半轴的平均偏差以及整体的偏差水平已经得到很好的校正。其中,终端束流的计算公式为:As shown in Figure 6, the beam distribution data of the striped ionization chamber before and after correction, the beam uniformity increased from 76.6% to 93.6%. It can be seen from the figure that the zero-crossing distortion, the average deviation of the positive and negative semi-axis, and the overall deviation level have been well corrected. Among them, the calculation formula of terminal beam current is:
其中,Amax和Amin分别是分条电离室束流分布数据中最大和最小的计数,U是终端束流均匀性。Among them, A max and A min are the maximum and minimum counts in the beam distribution data of the striped ionization chamber, respectively, and U is the terminal beam uniformity.
实施例二Embodiment two
基于相同的发明构思,本实施例公开了一种扫描磁铁动态识辨实时分段斜率反馈系统,包括:Based on the same inventive concept, this embodiment discloses a scanning magnet dynamic identification real-time segmented slope feedback system, including:
束流位置分布信息获取模块,用于获取束流在水平和垂直方向的束流位置分布信息;The beam position distribution information acquisition module is used to obtain the beam position distribution information of the beam in the horizontal and vertical directions;
误差信息分类模块,用于根据束流位置分布信息获得不同种类的误差信息;The error information classification module is used to obtain different types of error information according to the beam position distribution information;
误差优化模块,用于根据误差信息进行分段斜率反馈,得到优化后的电源三角波励磁曲线;The error optimization module is used to perform segmented slope feedback according to the error information to obtain the optimized power triangular wave excitation curve;
输出反馈模块,用于将电源三角波励磁曲线输入磁铁电源,束流再次经过后,获得新的分段斜率反馈,直至终端的扫描束斑满足预设的均匀性。The output feedback module is used to input the triangular wave excitation curve of the power supply into the magnet power supply. After the beam current passes through again, a new segmented slope feedback is obtained until the scanning beam spot at the terminal meets the preset uniformity.
误差优化模块包括:正负零漂模块,用于通过正负零漂解决正负半轴偏差;过零截距模块,用于通过过零截距解决过零区域凸起;分段斜率优化模块,用于将经过正负零漂和过零截距后的电源三角波按照斜率进行分段,并根据分段结果对电源三角波进行优化;反馈模块,用于将经过优化的电源三角波发送至电源和扫描磁铁,束流再次经过后,对分段斜率优化的结果进行反馈。The error optimization module includes: positive and negative zero drift module, used to solve positive and negative semi-axis deviation through positive and negative zero drift; zero crossing intercept module, used to solve zero crossing area bulge through zero crossing intercept; segmental slope optimization module , used to segment the triangular wave of the power supply after positive and negative zero drift and zero-crossing intercept according to the slope, and optimize the triangular wave of the power supply according to the segmentation results; the feedback module is used to send the optimized triangular wave of the power supply to the power supply and After the magnet is scanned and the beam passes through again, the result of segmented slope optimization is fed back.
实施例三Embodiment Three
基于相同的发明构思,本实施例公开了一种存储一个或多个程序的计算机可读存储介质,一个或多个程序包括指令,指令当由计算设备执行时,使得计算设备执行根据上述任一项非谐振快加速全波形动态补偿方法。Based on the same inventive concept, this embodiment discloses a computer-readable storage medium storing one or more programs, one or more programs include instructions, and the instructions, when executed by a computing device, cause the computing device to execute any one of the above-mentioned Term non-resonant fast acceleration full waveform dynamic compensation method.
实施例四Embodiment four
基于相同的发明构思,本实施例公开了一种计算设备,包括:一个或多个处理器、存储器以及一个或多个程序,其中一个或多个程序存储在存储器中并被配置为由一个或多个处理器执行,一个或多个程序包括用于执行根据上述任一项非谐振快加速全波形动态补偿方法。Based on the same inventive concept, this embodiment discloses a computing device, including: one or more processors, memory, and one or more programs, wherein one or more programs are stored in the memory and configured to be executed by one or more Executed by multiple processors, one or more programs are used to execute the full waveform dynamic compensation method for non-resonant fast acceleration according to any one of the above.
本领域内的技术人员应明白,本申请的实施例可提供为方法、系统、或计算机程序产品。因此,本申请可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且,本申请可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。Those skilled in the art should understand that the embodiments of the present application may be provided as methods, systems, or computer program products. Accordingly, the present application may take the form of an entirely hardware embodiment, an entirely software embodiment, or an embodiment combining software and hardware aspects. Furthermore, the present application may take the form of a computer program product embodied on one or more computer-usable storage media (including but not limited to disk storage, CD-ROM, optical storage, etc.) having computer-usable program code embodied therein.
本申请是参照根据本申请实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器,使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。The present application is described with reference to flowcharts and/or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of the present application. It should be understood that each procedure and/or block in the flowchart and/or block diagram, and a combination of procedures and/or blocks in the flowchart and/or block diagram can be realized by computer program instructions. These computer program instructions may be provided to a general purpose computer, special purpose computer, embedded processor, or processor of other programmable data processing equipment to produce a machine such that the instructions executed by the processor of the computer or other programmable data processing equipment produce a An apparatus for realizing the functions specified in one or more procedures of the flowchart and/or one or more blocks of the block diagram.
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品,该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。These computer program instructions may also be stored in a computer-readable memory capable of directing a computer or other programmable data processing apparatus to operate in a specific manner, such that the instructions stored in the computer-readable memory produce an article of manufacture comprising instruction means, the instructions The device realizes the function specified in one or more procedures of the flowchart and/or one or more blocks of the block diagram.
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上,使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理,从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。These computer program instructions can also be loaded onto a computer or other programmable data processing device, causing a series of operational steps to be performed on the computer or other programmable device to produce a computer-implemented process, thereby The instructions provide steps for implementing the functions specified in the flow chart or blocks of the flowchart and/or the block or blocks of the block diagrams.
最后应当说明的是:以上实施例仅用以说明本发明的技术方案而非对其限制,尽管参照上述实施例对本发明进行了详细的说明,所属领域的普通技术人员应当理解:依然可以对本发明的具体实施方式进行修改或者等同替换,而未脱离本发明精神和范围的任何修改或者等同替换,其均应涵盖在本发明的权利要求保护范围之内。上述内容仅为本申请的具体实施方式,但本申请的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本申请揭露的技术范围内,可轻易想到变化或替换,都应涵盖在本申请的保护范围之内。因此,本申请的保护范围应以权利要求的保护范围。Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention and not to limit them. Although the present invention has been described in detail with reference to the above embodiments, those of ordinary skill in the art should understand that: the present invention can still be Any modifications or equivalent replacements that do not depart from the spirit and scope of the present invention shall fall within the protection scope of the claims of the present invention. The above content is only the specific implementation of the application, but the scope of protection of the application is not limited thereto. Any person familiar with the technical field can easily think of changes or replacements within the technical scope disclosed in the application, and should cover Within the protection scope of this application. Therefore, the protection scope of the present application should be defined as the protection scope of the claims.
Claims (10)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211454473.XA CN115758735B (en) | 2022-11-21 | 2022-11-21 | Scanning magnet dynamic identification real-time segmentation slope feedback method and system |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211454473.XA CN115758735B (en) | 2022-11-21 | 2022-11-21 | Scanning magnet dynamic identification real-time segmentation slope feedback method and system |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115758735A true CN115758735A (en) | 2023-03-07 |
CN115758735B CN115758735B (en) | 2023-07-07 |
Family
ID=85333360
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202211454473.XA Active CN115758735B (en) | 2022-11-21 | 2022-11-21 | Scanning magnet dynamic identification real-time segmentation slope feedback method and system |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115758735B (en) |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20020081788A1 (en) * | 2000-09-15 | 2002-06-27 | Varian Semiconductor Equipment Associates, Inc. | Ion implanter optimized scan waveform retention and recovery |
CN101764029A (en) * | 2008-12-04 | 2010-06-30 | 北京中科信电子装备有限公司 | Method and device for accurately detecting and correcting parallelism of ion beam |
CN114377307A (en) * | 2021-12-28 | 2022-04-22 | 中国原子能科学研究院 | Device and method for quickly correcting beam position in pencil beam scanning |
-
2022
- 2022-11-21 CN CN202211454473.XA patent/CN115758735B/en active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20020081788A1 (en) * | 2000-09-15 | 2002-06-27 | Varian Semiconductor Equipment Associates, Inc. | Ion implanter optimized scan waveform retention and recovery |
CN101764029A (en) * | 2008-12-04 | 2010-06-30 | 北京中科信电子装备有限公司 | Method and device for accurately detecting and correcting parallelism of ion beam |
CN114377307A (en) * | 2021-12-28 | 2022-04-22 | 中国原子能科学研究院 | Device and method for quickly correcting beam position in pencil beam scanning |
Non-Patent Citations (1)
Title |
---|
张文铎;邱瑞昌;: "扫描磁铁电源的电流跟踪控制策略研究", 机械与电子, no. 01 * |
Also Published As
Publication number | Publication date |
---|---|
CN115758735B (en) | 2023-07-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN100552865C (en) | ion implanter | |
US6822405B2 (en) | Deceleration of hadron beams in synchrotrons designed for acceleration | |
JP5336496B2 (en) | Particle beam transport apparatus and particle beam transport method having a small beam spot size | |
US10361066B2 (en) | Ion implantation apparatus | |
JP5490310B2 (en) | Sensitivity correction method for dose monitoring apparatus and particle beam therapy apparatus | |
Zhang et al. | Design and commissioning of the beam test facility at the spallation neutron source | |
Pommerenke | Compact radio-frequency quadrupoles for industrial and medical applications | |
JP2002058750A (en) | Charged beam irradiation method and apparatus, and computer readable storage medium | |
CN115758735B (en) | Scanning magnet dynamic identification real-time segmentation slope feedback method and system | |
JP2015041440A (en) | Ion implantation device, magnetic field measuring device, and ion implantation method | |
Gu et al. | The electron lens test bench for the relativistic heavy ion collider at Brookhaven National Laboratory | |
JP2022122112A (en) | Ion implanting device and ion implanting method | |
Kazantseva et al. | Accurate Taylor transfer maps for large aperture iron dominated magnets used in charged particle separators and spectrometers | |
JP2019195408A (en) | Scanning irradiation device and particle beam therapy system, and adjustment method of scanning irradiation device | |
JP2018160695A (en) | Electron beam processing machine | |
WO2017221360A1 (en) | Particle therapy apparatus | |
JP2020069302A (en) | METHOD FOR DETERMINING OPERATING CONDITIONS OF PARTICLE BEAM TREATMENT DEVICE, AND PARTICLE BEAM TREATMENT DEVICE | |
Wang et al. | Commissioning of the 123 MeV injector for 12 GeV CEBAF | |
Witte et al. | Halbach magnets for CBETA and eRHIC | |
Rosenthal | Investigation of beam and spin dynamics for EDM measurements at COSY | |
US20240359035A1 (en) | Particle beam control system and particle beam control method | |
Amirkhanov et al. | Mathematical simulation of the main operation mode of the AIC-144 multipurpose isochronous cyclotron | |
Zhang et al. | Analysis and mitigation of bump magnetic field tracking errors in a multiturn injection system | |
EP4180089A1 (en) | Particle therapy system and irradiation control apparatus | |
Gabor et al. | Design Report of a Non–Destructive Emittance Instrument for Rutherford Appleton Laboratory’s Front End Test Stand FETS |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |