CN105282956B - 一种强流回旋加速器高频系统智能自启动方法 - Google Patents

一种强流回旋加速器高频系统智能自启动方法 Download PDF

Info

Publication number
CN105282956B
CN105282956B CN201510649568.0A CN201510649568A CN105282956B CN 105282956 B CN105282956 B CN 105282956B CN 201510649568 A CN201510649568 A CN 201510649568A CN 105282956 B CN105282956 B CN 105282956B
Authority
CN
China
Prior art keywords
frequency
value
cavity
pulse width
kalman filter
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.)
Active
Application number
CN201510649568.0A
Other languages
English (en)
Other versions
CN105282956A (zh
Inventor
殷治国
付晓亮
纪彬
赵振鲁
张天爵
魏俊逸
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China Institute of Atomic of Energy
Original Assignee
China Institute of Atomic of Energy
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by China Institute of Atomic of Energy filed Critical China Institute of Atomic of Energy
Priority to CN201510649568.0A priority Critical patent/CN105282956B/zh
Publication of CN105282956A publication Critical patent/CN105282956A/zh
Application granted granted Critical
Publication of CN105282956B publication Critical patent/CN105282956B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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。所述的采样保持单元,在脉冲宽度调制的脉宽内通过低阻回路快速充电,在脉冲宽度调制的脉宽间隙内通过高阻回路缓慢放电,放电回路时间常数远大于脉冲宽度调制周期T,以此保证整个脉冲宽度调制周期内的腔体电压信号检波输出值与连续波模式一致。
(2)将步骤(1)中得到的失谐角φ0通过数模转换器(ADC)量化为16位数字量,输入数字信号处理器(DSP)中经过φ-P计算得到腔体在多电子效应影响下的谐振位置P0,P0即为卡尔曼滤波器的测量值ZK,φ-P计算为线性映射计算P=kφ+b,φ∈[31,341],b∈R,其中比例系数k和偏移量b由实际系统的测量值经过计算得到。
(3)根据高频腔体锻炼的数据,统计高频腔体冷状态下极小脉冲功率驱动下无多电子效应的谐振点的微调电容位置P,作为卡尔曼滤波器的初始化值X(-1|-1),所述的卡尔曼滤波器的初始化值X(-1|-1),是高频腔体冷状态下,以D=1/10的脉冲宽度调制高频信号驱动高频腔体,高频驱动信号幅度极小,此时高频腔体处在次级电子倍增谐振场强区域外,该测量结果不受次级电子倍增效应影响,但会随腔体温度变化而变化。所述的卡尔曼滤波器为一维卡尔曼滤波器,计算输入和结果均为一维标量,其初始化值采用高频腔体锻炼数据的统计结果作为初始化值。
(4)根据卡尔曼滤波器的理论,k时刻的预测值由(k-1)时刻的值和k时刻的测量值共同计算获得,k时刻的测量值为步骤(2)中得到的测量值ZK,由此,k时刻的卡尔曼滤波器输出值为:
Pk|k=(I-KkHk)Pk|k-1
其中,是在k时刻的状态估计值,即k时刻的滤波器输出,
是k时刻的状态预测值,由k-1时刻的输出和k时刻的控制输入Bkuk计算得来:
Kk为最优卡尔曼增益,该值是使卡尔曼滤波器的均方误差最小的值,由下式定义:
为测量余量,由下式定义:
Pk|k是后验估计误差协方差矩阵,用来度量估计值的精确程度,
I为单位矩阵,
Pk|k-1为k时刻的预测估计协方差矩阵,由下式定义:
Qk为过程噪声的协方差矩阵,可由系统参数计算得出,
其中,测量值zk满足zk=HkXk+Vk,Vk是观测噪声,其均值为零,协方差矩阵为Rk,且服从正态分布。
为测量余量协方差,其中Hk是观测模型,将真实状态空间映射成观测空间,是Sk的逆矩阵,是Hk的转置矩阵。
Fk是状态变换矩阵。Bk是作用在控制器向量uk上的输入-控制模型。
最终,以卡尔曼滤波器k时刻的输出作为微调电容的目的位置P1
(5)以步骤(4)中的滤波器输出为目标地址,DSP控制微调电容运动到目的位置P1。等待脉冲幅度Amp在正反馈作用下超过状态转变阈值Ampsetpoint,系统由PWM调制状态进入连续波状态,状态转变阈值Ampsetpoint与具体的高频腔体多电子效应区域分布有关,该值由高频腔体锻炼数据总结得出;如果等待时间大于等待时间上限,则返回步骤(2)循环执行算法,等待时间上限为人为经验设定值,与高频腔体实际锻炼情况有关。
(6)在连续波状态下,逐步提升功率至加速电压达到束流加速要求值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,设定卡尔曼滤波器的初始化值X(-1|-1)为此值。设定等待超时上限timeout为2s,即指定位置2s后高频系统未转到连续状态运行,则重新计算谐振位置,重复以上过程。
图1是高频系统智能自启动方法在100MeV回旋加速器上真机启动过程测试图。其中,黄色信号为反射功率取样波形,绿色信号为腔体电压取样波形。该波形图明确反映出PWM信号占空比逐渐展宽的过程,在此过程中反射信号逐渐减小,直至最终由PWM模式切换到CW模式。在由PWM模式切换到CW模式的瞬间,RF功率突变为原来的2倍。
图2是智能自启动方法采样-保持单元的原理示意图。该单元的作用是在脉冲内快速充电,使输出信号迅速达到腔体取样信号的峰值;在脉冲间隙内,没有RF激励信号,腔体反馈信号为0,采样-保持单元通过高阻回路缓慢放电,使其输出与脉冲内输出维持一致。这一特性使智能自启动装置在工作状态切换时,保持模式切换工作点一致,避免了在CW模式下进入多电子效应区而引起的反射功率过大而对设备造成损害。
图3是智能自启动方法采用的卡尔曼滤波器的递推算法流程图,该滤波算法通过对预测值的校正,得到当前时刻的微调电容位置的最优估计。
图4为强流回旋加速器高频系统智能自启动方法示意图。
显然,本领域的技术人员可以对本发明进行各种改动和变型而不脱离本发明的精神和范围。这样,倘若对本发明的这些修改和变型属于本发明权利要求及其同等技术的范围之内,则本发明也意图包含这些改动和变型在内。

Claims (5)

1.一种强流回旋加速器高频系统智能自启动方法,其特征在于,所述方法包括如下步骤:
(1)以占空比为D=1/10的脉冲宽度调制高频信号驱动高频腔体,以脉冲宽度调制信号的上升沿触发幅度和失谐角采样保持单元采样得到高频腔体在脉冲内的幅度Amp和脉冲内的失谐角φ0
(2)将步骤(1)中得到的失谐角φ0通过数模转换器量化为16位数字量,输入数字信号处理器中经过φ-P计算得到腔体在多电子效应影响下的谐振位置P0,P0即为卡尔曼滤波器的测量值Zk,φ-P计算为线性映射计算P=kφ+b,φ∈[31,341],b∈R,其中比例系数k和偏移量b由实际系统的测量值经过计算得到;
(3)根据高频腔体锻炼的数据,统计高频腔体冷状态下无多电子效应的谐振点的微调电容位置P,作为卡尔曼滤波器的初始化值X(-1|-1);
(4)根据卡尔曼滤波器的理论,k时刻的预测值由k-1时刻的值和k时刻的测量值共同计算获得,k时刻的测量值为步骤(2)中得到的测量值Zk,Zk作为卡尔曼滤波器的输入值,由此,k时刻的卡尔曼滤波器输出值为:
Pk|k=(I-KkHk)Pk|k-1
其中,是在k时刻的状态估计值,即k时刻的滤波器输出,
是k时刻的状态预测值,由k-1时刻的输出和k时刻的控制输入Bkuk计算得来:
Kk为最优卡尔曼增益,由下式定义:
为测量余量,由下式定义:
Pk|k是后验估计误差协方差矩阵,用来度量估计值的精确程度,
I为单位矩阵,
Pk|k-1为k时刻的预测估计协方差矩阵,由下式定义:
Qk为过程噪声的协方差矩阵,可由系统参数计算得出,
其中,测量值Zk满足Zk=HkXk+Vk,Vk是观测噪声,其均值为零,协方差矩阵为Rk,且服从正态分布;
为测量余量协方差,其中Hk是观测模型,将真实状态空间映射成观测空间,是Sk的逆矩阵,是Hk的转置矩阵;
Fk是状态变换矩阵,Bk是作用在控制器向量uk上的输入-控制模型;
最终,以卡尔曼滤波器k时刻的输出作为微调电容的目的位置P1
(5)以步骤(4)中的滤波器输出为目标地址,数字信号处理器控制微调电容运动到目的位置P1;等待脉冲幅度在正反馈作用下超过状态转变阈值Ampsetpoint,系统由脉冲宽度调制状态进入连续波状态,状态转变阈值Ampsetpoint与具体的高频腔体多电子效应区域分布有关,该值由高频腔体锻炼数据总结得出;如果等待时间大于等待时间上限,则返回步骤(2)循环执行算法,等待时间上限为人为经验设定值,与高频腔体实际锻炼情况有关;
(6)在连续波状态下,逐步提升功率至加速电压达到束流加速要求值V0,系统进行幅度闭环控制,高频系统启动过程完成。
2.根据权利要求1所述的一种强流回旋加速器高频系统智能自启动方法,其特征在于:步骤(1)中所述的脉冲宽度调制信号频率为20KHz。
3.根据权利要求1所述的一种强流回旋加速器高频系统智能自启动方法,其特征在于:步骤(1)中所述的采样保持单元,在脉冲宽度调制的脉宽内通过低阻回路快速充电,在脉冲宽度调制的脉宽间隙内通过高阻回路缓慢放电,放电回路时间常数远大于脉冲宽度调制周期T,以此保证整个脉冲宽度调制周期内的腔体电压信号检波输出值与连续波模式一致。
4.根据权利要求1所述的一种强流回旋加速器高频系统智能自启动方法,其特征在于:步骤(3)所述的卡尔曼滤波器的初始化值X(-1|-1),是高频腔体冷状态下,以D=1/10的脉冲宽度调制高频信号驱动高频腔体,高频驱动信号幅度极小,此时高频腔体处在次级电子倍增谐振场强区域外,该测量结果不受次级电子倍增效应影响,但会随腔体温度变化而变化。
5.根据权利要求1所述的一种强流回旋加速器高频系统智能自启动方法,其特征在于:步骤(4)所述的卡尔曼滤波器为一维卡尔曼滤波器,计算输入和结果均为一维标量,其初始化值采用高频腔体锻炼数据的统计结果作为初始化值。
CN201510649568.0A 2015-10-09 2015-10-09 一种强流回旋加速器高频系统智能自启动方法 Active CN105282956B (zh)

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 CN105282956A (zh) 2016-01-27
CN105282956B true 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)

Families Citing this family (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109100567B (zh) * 2018-06-27 2020-06-23 中国原子能科学研究院 同步回旋加速器调制频率测试方法
CN108650771A (zh) * 2018-07-12 2018-10-12 北京航天广通科技有限公司 加速器超导腔的自激稳定控制方法、系统和存储介质
CN108815723B (zh) * 2018-07-26 2024-01-30 中国科学院近代物理研究所 高频腔失谐检测单元及其检测方法
CN109561567B (zh) * 2018-11-27 2020-01-24 中国原子能科学研究院 一种大功率高稳定度变负载高频加速系统
CN110536536B (zh) * 2019-09-03 2020-06-23 中国原子能科学研究院 回旋加速器低电平系统重启动自动锻炼装置和控制方法
CN110913560B (zh) * 2019-12-09 2024-05-31 中国原子能科学研究院 带电粒子回旋加速器的腔体锻炼加速装置与方法、回旋加速器
CN112947363B (zh) * 2021-01-29 2023-06-20 中国工程物理研究院流体物理研究所 一种基于cRIO的强流加速器机器快保护系统
CN115426761A (zh) * 2022-09-14 2022-12-02 中国原子能科学研究院 一种回旋加速器优化启动装置

Family Cites Families (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH05144600A (ja) * 1991-11-25 1993-06-11 Toshiba Corp ビーム信号模擬装置
EP1265462A1 (fr) * 2001-06-08 2002-12-11 Ion Beam Applications S.A. Dispositif et méthode de régulation de l'intensité d'un faisceau extrait d'un accélérateur de particules
JP5046928B2 (ja) * 2004-07-21 2012-10-10 メヴィオン・メディカル・システムズ・インコーポレーテッド シンクロサイクロトロン及び粒子ビームを生成する方法
CN101420247B (zh) * 2008-11-20 2012-06-13 湖南华诺星空电子技术有限公司 开环体制超宽带脉冲源阵列时域波束合成与扫描方法
CN102984878B (zh) * 2012-11-28 2015-04-29 中国原子能科学研究院 医用回旋加速器的多态调谐方法
CN102933021B (zh) * 2012-11-28 2015-04-29 中国原子能科学研究院 医用回旋加速器的射频启动方法
US9119281B2 (en) * 2012-12-03 2015-08-25 Varian Medical Systems, Inc. Charged particle accelerator systems including beam dose and energy compensation and methods therefor
CN103025044B (zh) * 2012-12-27 2015-04-29 中国原子能科学研究院 医用回旋加速器腔体打火保护的方法及所用装置
CN103079334B (zh) * 2013-01-04 2016-06-22 中国原子能科学研究院 回旋加速器射频谐振腔体自动锻炼系统

Also Published As

Publication number Publication date
CN105282956A (zh) 2016-01-27

Similar Documents

Publication Publication Date Title
CN105282956B (zh) 一种强流回旋加速器高频系统智能自启动方法
CN105826154B (zh) 针对脉冲射频电源的阻抗匹配方法及装置
EP2206215B1 (en) Method and apparatus to determine battery resonance
CN103079334B (zh) 回旋加速器射频谐振腔体自动锻炼系统
CN115777167A (zh) 用于电池胞元的阻抗测量的系统和方法
CN105245151B (zh) 表贴式永磁同步电机转子位置的检测方法
US20230386790A1 (en) Impedance-matching method, impedance-matching device, and semiconductor process apparatus
CN103414418B (zh) 一种直流电机电流pi控制系统的控制方法
US20240282553A1 (en) Base bias adjustment apparatus and method, and semiconductor process device
CN110536536B (zh) 回旋加速器低电平系统重启动自动锻炼装置和控制方法
Engelen et al. Gain scheduling control of a bidirectional dc–dc converter with large dead‐time
CN116936328A (zh) 阻抗匹配方法、装置和介质
WO2014173363A1 (zh) 开关整流器启动控制方法及其装置、存储介质
CN104767290A (zh) 一种感应式无线供电系统双闭环恒定输出电流控制方法
JP2021501452A (ja) マグネトロン温度調節方法、装置及びシステム、可変周波数電源及びマイクロ波機器
Malla et al. DTC-SVM of induction motor by applying two fuzzy logic controllers
CN105188178B (zh) 一种磁控管退灯丝电压的控制系统及方法
CN114513871B (zh) 一种微波功率调节方法、装置、电子设备及存储介质
CN106851886B (zh) 确定磁控管的退灯丝电压的方法、装置及磁控管、加速器
CN111060912A (zh) 一种天气雷达速调管仿真的方法
CN108718169A (zh) 一种异步电机分级变频软起动装置机控制方法及系统
CN110456190B (zh) 一种磁控管状态检测方法、控制器、电路及系统
CN111130309B (zh) 一种提升采样精度的控制方法、存储介质及逆变电源
CN114725887A (zh) 过流保护方法、装置、设备及存储介质
CN209994271U (zh) 一种改善电源瞬间加载过功率的电源

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