CN105282956A - 一种强流回旋加速器高频系统智能自启动方法 - Google Patents
一种强流回旋加速器高频系统智能自启动方法 Download PDFInfo
- Publication number
- CN105282956A CN105282956A CN201510649568.0A CN201510649568A CN105282956A CN 105282956 A CN105282956 A CN 105282956A CN 201510649568 A CN201510649568 A CN 201510649568A CN 105282956 A CN105282956 A CN 105282956A
- Authority
- CN
- China
- Prior art keywords
- cavity
- value
- radio frequency
- kalman filter
- frequency system
- 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
Landscapes
- Particle Accelerators (AREA)
Abstract
本发明涉及回旋加速器的高频控制技术,具体涉及一种强流回旋加速器高频系统智能自启动方法。该方法以PWM调制的高频信号驱动高频腔体,通过Sample-Hold技术取得高频腔体取样信号在脉冲内的幅度和失谐角,得到高频腔体在次级电子倍增效应影响下的谐振位置,采用卡尔曼滤波器对调谐环输出的失谐角进行滤波,得到高频腔体的近似实际谐振位置。DSP根据卡尔曼滤波器的输出为目标,设定微调电容到达指定的位置,并在该目标位置开始启动过程。这种高频系统智能自启动方法可以提高强流回旋加速器高频系统的启动效率,缩短启动时间,有效避免次级电子倍增效应对高频系统的影响。
Description
技术领域
本发明涉及回旋加速器的高频控制技术,具体涉及一种强流回旋加速器高频系统智能自启动方法。
背景技术
高频腔体是回旋加速器高频系统的基本组成之一,粒子在回旋加速器中回旋运动,经过高频腔体的加速间隙时获得能量被加速。对于高频腔体,当腔体结构、频率及功率满足一定条件时将产生多电子效应。该效应阻止高频功率馈入高频腔体,阻碍加速电场的建立,并有可能对腔体、耦合窗等结构产生破坏。一般情况下,在高频系统整体调试运行前对腔体等结构做清洁、锻炼等预处理,减小次级电子倍增效应发生的可能性。但此效应并非不能单次非可逆克服,当腔体再次与大气接触或污损将导致预处理过程的逆变。因此,在启动阶段克服腔体内多电子效应,以最快的速度建立高频电场是高频系统关注重点。从传输线的角度看,次级电子倍增效应发生时等效于改变了高频系统的等效负载,高频匹配的恶化将导致高频信号驻波比增大,触发低电平系统或者发射机的反射保护,甚至对高频系统设备造成损害,对高频系统的启动提出较高的要求。
为防止启动期间次级电子倍增效应引起的反射功率过大对设备造成损害,回旋加速器一般采用PWM调制高频信号的方式启动。在PWM脉冲内,幅度检波和失谐角检测器能够根据腔体反馈信号对当前腔体所处状态进行评估,确定腔体当前所处状态,得到腔体当前失谐角。在PWM脉冲间隙内,由于采用Sample-Hold单元,幅度检波和失谐角监测电路输出值与PWM脉冲内一致,使后续电路持续工作。然而由于次级电子倍增效应的影响,幅度检波和失谐角监测均受到影响,该影响可以归结为符合高斯分布的噪声影响,因此腔体精确谐振位置不能够准确测量。卡尔曼滤波器是一个最优化自回归数据处理算法,该算法的一个典型实例是从一组有限的,包含噪声的,对物体位置的观察序列(可能有偏差)预测出物体的位置的坐标及速度。卡尔曼滤波利用目标的动态信息,设法去掉噪声的影响,得到一个关于目标位置的概率上最好的估计。这个估计可以是对当前目标位置的估计(滤波),也可以是对于将来位置的估计(预测),也可以是对过去位置的估计(插值或平滑)。本发明中采用的是对当前目标位置的估计,实现滤除次级电子倍增效应的目的。在得到卡尔曼滤波器的最优结果后,控制微调电容运动到指定位置,在该位置的启动速度将快于其他位置的启动速度,并可以最大概率的避免打火,反射功率过大等异常情况。
强流回旋加速器高频系统智能自启动方法是解决高频系统启动过程中由于次级电子倍增影响而导致的启动失败,反射功率过大等异常情况的一种途径,是实现无人干预下自动处理各种异常情况,快速安全启动高频系统、建立加速电压的方法。
发明内容
本发明的目的在于针对目前强流回旋加速器高频腔体存在次级电子倍增效应影响高频系统正常启动的问题,通过对腔体和多电子效应的研究,提出了一种强流回旋加速器高频系统智能自启动方法,该方法能够减少次级电子倍增效应对调谐环输出结果的影响,加速高频系统的启动过程,有效解决次级电子倍增效应区导致的大功率反射和暗电流导致的失配等加速器高频工程技术困难。
本发明的技术方案如下:一种强流回旋加速器高频系统智能自启动方法,所述方法包括如下步骤:
(1)以占空比为D=1/10的PWM调制高频信号驱动高频腔体。以PWM调制信号的上升沿触发幅度和失谐角Sample-Hold单元采样得到高频腔体在脉冲内的幅度Amp和脉冲内的失谐角φ0。
(2)将步骤(1)中得到的失谐角φ0通过数模转换器(ADC)量化为16位数字量,输入DSP中经过φ-P计算得到腔体在多电子效应影响下的谐振位置P0,该值即为测量值ZK。
(3)根据高频腔体锻炼的数据,统计高频腔体冷状态下极小脉冲功率驱动下无多电子效应的谐振点P,作为卡尔曼滤波器的初始化值Χ(-1|-1)。
(4)根据卡尔曼滤波器的理论,t时刻的预测值由(t-1)时刻的值和t时刻的测量值估计获得,t时刻的测量值为步骤(2)中得到的测量值ZK,由此,t时刻的卡尔曼滤波器输出值为:
Pk|k=(I-KkHk)Pk|k-1
其中Kk为最优卡尔曼增益,为测量余量,由下式定义:
上式中,测量值zk满足zk=HkΧk+Vk,Vk是观测噪声,其均值为零,协方差矩阵为Rk,且服从正态分布。
为测量余量协方差,其中Hk是观测模型,将真实状态空间映射成观测空间。
下一状态预测值为:
其中,Fk是作用在Xk-1上的状态变换模型。Bk是作用在控制器向量uk上的输入-控制模型。
最终,以卡尔曼滤波器的输出作为微调电容的目的谐振位置P1。
(5)以步骤(4)中的滤波器输出为目标地址,DSP控制微调电容运动到指定位置P1。等待Amp在正反馈作用下超过状态转变阈值Ampsetpoint,系统由PWM调制状态进入CW状态;如果等待时间大于timeout上限,则返回步骤(2)循环执行算法。
(1)在CW状态下,逐步提升功率至加速电压达到束流加速要求值V0,系统
幅度闭环控制。高频系统启动过程完成。
本发明针对强流回旋加速器高频系统的需求,已在我国目前能量最高的质子回旋加速器CYCIAE-100工程实践中显示出十分有益的效果。本装置能够应对CYCIAE-100许多个多电子倍增效应区导致的大功率反射和暗电流导致的失配等大型加速器高频工程技术困难,实现高频系统的智能自启动以及运行过程中的异常保护。因为本发明的成功应用,使得CYCIAE-100能在毫秒量级时间范围内从打火等故障中自动恢复到正常加速状态,是我国大型带电粒子回旋加速器的顺利调试取得成功的核心技术之一。
附图说明
图1是CYCIAE-100高频系统启动过程波形图;
图2是采样-保持单元的工作原理图;
图3是卡尔曼滤波器的递推算法流程图;
图4为强流回旋加速器高频系统智能自启动方法示意图。
具体实施方式
下面结合附图和实施例对本发明进行详细的描述。
以100MeV强流质子回旋加速器CYCIAE-100为例,介绍其高频系统智能自启动方法的设计及实际验证。该加速器实际运行所需功率约64kW,包括腔体损耗62kW及2kW的束流功率。腔体取样耦合度为-50.0dB,由于传输线存在损耗,实际测量得额定功率运行时,腔体取样为8.4Vpp。智能自启动方法以高频脉冲方式驱动腔体,占空比1:10,周期20ms。根据锻炼数据,该回旋加速器高频腔体在取样值大于5Vpp之后,次级电子倍增效应影响可以忽略。智能自启动方法设定状态转变阈值Ampsetpoint=5.5Vpp作为PWM模式和CW模式切换的工作点,以保证在没有多电子效应的区间上在PWM模式和CW模式之间相互切换。根据高频腔体锻炼数据,冷状态下高频腔体谐振位置为-8V,设定卡尔曼滤波器的初始化值Χ(-1|-1)为此值。设定等待超时上限timeout为2s,即指定位置2s后高频系统未转到连续状态运行,则重新计算谐振位置,重复以上过程。
图1是高频系统智能自启动方法在100MeV回旋加速器上真机启动过程测试图。其中,黄色信号为反射功率取样波形,绿色信号为腔体电压取样波形。该波形图明确反映出PWM信号占空比逐渐展宽的过程,在此过程中反射信号逐渐减小,直至最终由PWM模式切换到CW模式。在由PWM模式切换到CW模式的瞬间,RF功率突变为原来的2倍。
图2是智能自启动方法采样-保持单元的原理示意图。该单元的作用是在脉冲内快速充电,使输出信号迅速达到腔体取样信号的峰值;在脉冲间隙内,没有RF激励信号,腔体反馈信号为0,采样-保持单元通过高阻回路缓慢放电,使其输出与脉冲内输出维持一致。这一特性使智能自启动装置在工作状态切换时,保持模式切换工作点一致,避免了在CW模式下进入多电子效应区而引起的反射功率过大而对设备造成损害。
图3是智能自启动方法采用的卡尔曼滤波器的递推算法流程图,该滤波算法通过对预测值的校正,得到当前时刻的微调电容位置的最优估计。
图4为强流回旋加速器高频系统智能自启动方法示意图。
显然,本领域的技术人员可以对本发明进行各种改动和变型而不脱离本发明的精神和范围。这样,倘若对本发明的这些修改和变型属于本发明权利要求及其同等技术的范围之内,则本发明也意图包含这些改动和变型在内。
Claims (8)
1.一种强流回旋加速器高频系统智能自启动方法,其特征在于,所述方法包括如下步骤:
(1)以占空比为D=1/10的PWM调制高频信号驱动高频腔体,以PWM调制信号的上升沿触发幅度和失谐角Sample-Hold单元采样得到高频腔体在脉冲内的幅度Amp和脉冲内的失谐角φ0;
(2)将步骤(1)中得到的失谐角φ0通过数模转换器(ADC)量化为16位数字量,输入DSP中经过φ-P计算得到腔体在多电子效应影响下的谐振位置P0,该值即为测量值Zk;
(3)根据高频腔体锻炼的数据,统计高频腔体冷状态下极小脉冲功率驱动下无多电子效应的谐振点P,作为卡尔曼滤波器的初始化值Χ(-1|-1);
(4)根据卡尔曼滤波器的理论,t时刻的预测值由t-1时刻的值和t时刻的测量值估计获得,t时刻的测量值为步骤(2)中得到的测量值Zk,该值作为卡尔曼滤波器的输入值,由此,t时刻的卡尔曼滤波器输出值为:
Pk|k=(I-KkHk)Pk|k-1
其中Kk为最优卡尔曼增益,为测量余量,由下式定义:
上式中,测量值zk满足zk=HkΧk+Vk,Vk是观测噪声,其均值为零,协方差矩阵为Rk,且服从正态分布;
为测量余量协方差,其中Hk是观测模型,将真实状态空间映射成观测空间;
下一状态预测值为:
其中,Fk是作用在Xk-1上的状态变换模型,Bk是作用在控制器向量uk上的输入-控制模型;
最终,以卡尔曼滤波器的输出作为微调电容的目的谐振位置P1;
(5)以步骤(4)中的滤波器输出为目标地址,DSP控制微调电容运动到指定位置P1;等待Amp在正反馈作用下超过状态转变阈值Ampsetpoint,系统由PWM调制状态进入CW状态;如果等待时间大于timeout上限,则返回步骤(2)循环执行算法;
(6)在CW状态下,逐步提升功率至加速电压达到束流加速要求值V0,系统幅度闭环控制,高频系统启动过程完成。
2.根据权利要求1所述的一种强流回旋加速器高频系统智能自启动方法,其特征在于:步骤(1)中所述的PWM调制信号频率为20KHz。
3.根据权利要求1所述的一种强流回旋加速器高频系统智能自启动方法,其特征在于:步骤(1)中所述的Sample-Hold单元,在PWM脉宽内通过低阻回路快速充电,在PWM脉宽间隙内通过高阻回路缓慢放电,放电回路时间常数远大于PWM周期T,以此保证整个PWM周期内的腔体电压信号检波输出值与连续波模式一致。
4.根据权利要求1所述的一种强流回旋加速器高频系统智能自启动方法,其特征在于:步骤(2)中所述的φ-P计算为线性映射计算P=kφ+b,φ∈[31,341],b∈R,其中比例系数k和偏移量b由实际系统测量计算得到。
5.根据权利要求1所述的一种强流回旋加速器高频系统智能自启动方法,其特征在于:步骤(3)所述的卡尔曼滤波器的初始化值Χ(-1|-1),是高频腔体冷状态下,以D=1/10的PWM调制高频信号驱动高频腔体,高频驱动信号幅度极小,此时高频腔体处在次级电子倍增谐振场强区域外,该测量结果不受次级电子倍增效应影响,但会随腔体温度变化而变化。
6.根据权利要求1所述的一种强流回旋加速器高频系统智能自启动方法,其特征在于:步骤(4)所述的卡尔曼滤波器为一维卡尔曼滤波器,计算输入和结果均为一维标量,其初始化值不采用常规的Χ(-1|-1)=μs,其中Χ(-1|-1)~N(μs,Cs),而采用高频腔体锻炼数据的统计结果作为初始化值。
7.根据权利要求1所述的一种强流回旋加速器高频系统智能自启动方法,其特征在于:步骤(5)所述的状态转变阈值Ampsetpoint,与具体的高频腔体多电子效应区域分布有关,该值由高频腔体锻炼数据总结得出。
8.根据权利要求1所述的一种强流回旋加速器高频系统智能自启动方法,其特征在于:步骤(5)所述的等待时间上限timeout,为人为经验设定值,与高频腔体实际锻炼情况有关。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510649568.0A CN105282956B (zh) | 2015-10-09 | 2015-10-09 | 一种强流回旋加速器高频系统智能自启动方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510649568.0A CN105282956B (zh) | 2015-10-09 | 2015-10-09 | 一种强流回旋加速器高频系统智能自启动方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105282956A true CN105282956A (zh) | 2016-01-27 |
CN105282956B CN105282956B (zh) | 2018-08-07 |
Family
ID=55151051
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510649568.0A Active CN105282956B (zh) | 2015-10-09 | 2015-10-09 | 一种强流回旋加速器高频系统智能自启动方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105282956B (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108650771A (zh) * | 2018-07-12 | 2018-10-12 | 北京航天广通科技有限公司 | 加速器超导腔的自激稳定控制方法、系统和存储介质 |
CN108815723A (zh) * | 2018-07-26 | 2018-11-16 | 中国科学院近代物理研究所 | 高频腔失谐检测单元及其检测方法 |
CN109100567A (zh) * | 2018-06-27 | 2018-12-28 | 中国原子能科学研究院 | 同步回旋加速器调制频率测试方法 |
CN109561567A (zh) * | 2018-11-27 | 2019-04-02 | 中国原子能科学研究院 | 一种大功率高稳定度变负载高频加速系统 |
CN110536536A (zh) * | 2019-09-03 | 2019-12-03 | 中国原子能科学研究院 | 回旋加速器低电平系统重启动自动锻炼装置和控制方法 |
CN112947363A (zh) * | 2021-01-29 | 2021-06-11 | 中国工程物理研究院流体物理研究所 | 一种基于cRIO的强流加速器机器快保护系统 |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH05144600A (ja) * | 1991-11-25 | 1993-06-11 | Toshiba Corp | ビーム信号模擬装置 |
US20040155206A1 (en) * | 2001-06-08 | 2004-08-12 | Bruno Marchand | Device and method for regulating intensity of beam extracted from a particle accelerator |
CN101420247A (zh) * | 2008-11-20 | 2009-04-29 | 长沙华诺星空电子技术有限公司 | 开环体制超宽带脉冲源阵列时域波束合成与扫描方法 |
EP2259664A2 (en) * | 2004-07-21 | 2010-12-08 | Still River Systems, Inc. | A programmable radio frequency waveform generator for a synchrocyclotron |
CN102933021A (zh) * | 2012-11-28 | 2013-02-13 | 中国原子能科学研究院 | 医用回旋加速器的射频启动方法 |
CN102984878A (zh) * | 2012-11-28 | 2013-03-20 | 中国原子能科学研究院 | 医用回旋加速器的多态调谐方法 |
CN103025044A (zh) * | 2012-12-27 | 2013-04-03 | 中国原子能科学研究院 | 医用回旋加速器腔体打火保护的方法及所用装置 |
CN103079334A (zh) * | 2013-01-04 | 2013-05-01 | 中国原子能科学研究院 | 回旋加速器射频谐振腔体自动锻炼系统 |
WO2014088958A1 (en) * | 2012-12-03 | 2014-06-12 | Varian Medical Systems, Inc. | Charged particle accelerator systems including beam dose and energy compensation and methods therefor |
-
2015
- 2015-10-09 CN CN201510649568.0A patent/CN105282956B/zh active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH05144600A (ja) * | 1991-11-25 | 1993-06-11 | Toshiba Corp | ビーム信号模擬装置 |
US20040155206A1 (en) * | 2001-06-08 | 2004-08-12 | Bruno Marchand | Device and method for regulating intensity of beam extracted from a particle accelerator |
EP2259664A2 (en) * | 2004-07-21 | 2010-12-08 | Still River Systems, Inc. | A programmable radio frequency waveform generator for a synchrocyclotron |
CN101420247A (zh) * | 2008-11-20 | 2009-04-29 | 长沙华诺星空电子技术有限公司 | 开环体制超宽带脉冲源阵列时域波束合成与扫描方法 |
CN102933021A (zh) * | 2012-11-28 | 2013-02-13 | 中国原子能科学研究院 | 医用回旋加速器的射频启动方法 |
CN102984878A (zh) * | 2012-11-28 | 2013-03-20 | 中国原子能科学研究院 | 医用回旋加速器的多态调谐方法 |
WO2014088958A1 (en) * | 2012-12-03 | 2014-06-12 | Varian Medical Systems, Inc. | Charged particle accelerator systems including beam dose and energy compensation and methods therefor |
CN103025044A (zh) * | 2012-12-27 | 2013-04-03 | 中国原子能科学研究院 | 医用回旋加速器腔体打火保护的方法及所用装置 |
CN103079334A (zh) * | 2013-01-04 | 2013-05-01 | 中国原子能科学研究院 | 回旋加速器射频谐振腔体自动锻炼系统 |
Non-Patent Citations (3)
Title |
---|
刘太联等: "卡尔曼滤波在加速器电源控制中的应用研究", 《核电子学与探测技术》 * |
唐政等: "基于自适应卡尔曼滤波的机动目标跟踪算法", 《电子科技》 * |
纪彬等: "100MeV负氢回旋加速器高频冷模腔的理论与实验研究", 《CHINESE PHYSICS C(HEP & NP)》 * |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109100567A (zh) * | 2018-06-27 | 2018-12-28 | 中国原子能科学研究院 | 同步回旋加速器调制频率测试方法 |
CN108650771A (zh) * | 2018-07-12 | 2018-10-12 | 北京航天广通科技有限公司 | 加速器超导腔的自激稳定控制方法、系统和存储介质 |
CN108815723A (zh) * | 2018-07-26 | 2018-11-16 | 中国科学院近代物理研究所 | 高频腔失谐检测单元及其检测方法 |
CN108815723B (zh) * | 2018-07-26 | 2024-01-30 | 中国科学院近代物理研究所 | 高频腔失谐检测单元及其检测方法 |
CN109561567A (zh) * | 2018-11-27 | 2019-04-02 | 中国原子能科学研究院 | 一种大功率高稳定度变负载高频加速系统 |
CN110536536A (zh) * | 2019-09-03 | 2019-12-03 | 中国原子能科学研究院 | 回旋加速器低电平系统重启动自动锻炼装置和控制方法 |
CN110536536B (zh) * | 2019-09-03 | 2020-06-23 | 中国原子能科学研究院 | 回旋加速器低电平系统重启动自动锻炼装置和控制方法 |
CN112947363A (zh) * | 2021-01-29 | 2021-06-11 | 中国工程物理研究院流体物理研究所 | 一种基于cRIO的强流加速器机器快保护系统 |
CN112947363B (zh) * | 2021-01-29 | 2023-06-20 | 中国工程物理研究院流体物理研究所 | 一种基于cRIO的强流加速器机器快保护系统 |
Also Published As
Publication number | Publication date |
---|---|
CN105282956B (zh) | 2018-08-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105282956A (zh) | 一种强流回旋加速器高频系统智能自启动方法 | |
CN103079334B (zh) | 回旋加速器射频谐振腔体自动锻炼系统 | |
CN104104112B (zh) | 用于两级拓扑结构的光伏并网逆变器的mppt控制方法 | |
CN110380601A (zh) | 一种数字llc谐振变换器的软启动系统及方法 | |
CN103414418B (zh) | 一种直流电机电流pi控制系统的控制方法 | |
CN104765403B (zh) | 基于混沌优化搜索的光伏阵列最大功率点跟踪方法 | |
CN104965411A (zh) | 一种调平撑腿落地检测阀值的自适应调节方法 | |
Cho et al. | The KOMAC accelerator facility | |
CN111081506B (zh) | 一种速调管测试和老炼系统及相应方法 | |
CN114094901A (zh) | 一种动力电池脉冲加热过程中的电机转子位置控制方法 | |
CN211928161U (zh) | 一种速调管测试和老炼系统 | |
CN209994271U (zh) | 一种改善电源瞬间加载过功率的电源 | |
CN111060912A (zh) | 一种天气雷达速调管仿真的方法 | |
CN106851886B (zh) | 确定磁控管的退灯丝电压的方法、装置及磁控管、加速器 | |
Chong et al. | 550-W Ka-band pulsed helix TWT for radar applications | |
CN205405139U (zh) | 一种用于电动汽车单片机程序远程更新的硬件电路 | |
CN110536536B (zh) | 回旋加速器低电平系统重启动自动锻炼装置和控制方法 | |
CN114734845A (zh) | 超级充电枪温度控制方法、系统和基于其的超级充电枪 | |
Shimozuma et al. | Optimization of high power and high efficiency operation of 77GHz gyrotrons for ECRH in the large helical device | |
CN104393477A (zh) | 种子注入式激光器的正弦扫描谐振探测装置及其探测方法 | |
CN103427748B (zh) | 一种单相感应电机励磁电流pi控制方法 | |
CN116292080B (zh) | 风力发电机组的变桨距控制方法及装置 | |
Bustinduy et al. | The ess mebt rf buncher cavities conditioning process | |
CN107067614A (zh) | 一种红外对射装置 | |
Li | An Optimal Procedure for RF Conditionning at the FREIA Laboratory |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |