基于贝叶斯优化的高能质子空间传输精确控制.pdf

粒子束及加速器技术基于贝叶斯优化的高能质子空间传输精确控制*申诗雨1，杨晓虎1，张国博1，赵子琦1，马燕云2（1.国防科技大学理学院核科学与技术系，长沙410073；2.国防科技大学前沿交叉学科学院，长沙410073）摘要：综合考虑地磁场、高能质子的相对论效应以及同步辐射的影响，建立了质子在空间传输的单粒子运动模型。基于该模型，提出利用贝叶斯优化方法，实现给定质子能量下，质子从空间初始点传输到目标点的精确控制方法，获得了出射方向随出射位置的变化规律，当位置径向角取 0和 180时，位置轴向角的取值不会改变粒子的最优出射方向。研究结果为质子束在空间环境中的长程传输提供理论支撑。关 键 词：高能质子；传输控制；空间传输；贝叶斯优化中图分类号：O53文献标志码：Adoi:10.11884/HPLPB202335.230231Precisecontrolofhigh-energyprotonstransportinspaceenvironmentbyusingbayesianoptimizationShenShiyu1，YangXiaohu1，ZhangGuobo1，ZhaoZiqi1，MaYanyun2（1.Department of Nuclear Science and Technology,National University of Defense Technology,Changsha 410073,China;2.College of Advanced Interdisciplinary Studies,National University of Defense Technology,Changsha 410073,China）Abstract：Consideringthegeomagneticfield,therelativisticeffectandbremsstrahlungradiationofhigh-energyprotons,asingleparticlemotionmodelofprotontransportinthespaceenvironmentisestablished.Basedonthismodel,theBayesianoptimizationmethodisproposedtorealizetheprecisecontrolofprotonstransportfromtheinitialpositiontothetargetunderagivenprotonenergy.Thedependenceoftheprotonlaunchangleonthelaunchheightisobtained,thatis,whenthecoordinateradialangleis0and180,thevalueofthecoordinateaxialanglewillnotchangetheoptimalemissiondirectionoftheparticles.Theresultscanprovidetheoreticalreferencesforlong-distancetransportofprotonbeamsinthespaceenvironment.Keywords：highenergyproton,transportcontrol,spacetransmission,Bayesianoptimization在目前的空间主动实验中，已经完成的形式包括发射带电粒子束、释放等离子体、释放化学物质和核爆炸等。其中发射带电粒子对于辐射防护以及空间碎片清除等方面具有工程价值。由于探寻带电粒子在空间传输中的精确控制实验成本昂贵，因此有必要利用理论方法、数值模拟等简易、经济的研究手段，对质子束在空间环境中的长程传输控制开展预研和探索。由于高能电子束加速相比起高能质子束加速要相对容易，在国内外开展的研究中主要集中在电子束传输研究，对高能质子的研究相对较少。Neubert 等1-2使用蒙特卡罗模拟方法，分析了同向地磁场对束包络的影响；Mironychev 等3则结合蒙特卡罗方法计算了电子束水平射入大气环境的情况，指出垂直方向的地磁场会改变电子束的传输方向；Porazik 等4-5研究了电子束损失锥的计算以及电子束注入位置的选择，之后 Powis6利用单粒子运动模型模拟了电子束从轨道卫星发射后，沿磁力线传输并最终抵达地球附近的全过程，发现磁层环境的变化会显著影响电子束在大气层内束斑的落点。Hao 等7-8利用二维 PIC 方法，建立了适用于真空环境内相对论电子束长程传输的计算模型，总结了地磁场及自生电磁场对电子束长程传输的影响。姚海波等人的研究表明通过调制相对论电子束的束/背景等离子体密度比值、相对论因子、束包络和密度分布等，可以抑制电子束传输中的束流不稳定*收稿日期：2023-07-26；修订日期：2023-09-19基金项目：国家自然科学基金项目(12175309、11975308、12005297)；中国科学院战略先导专项(XDA25050200)；湖南省自然科学基金项目(CX20190001)；国防科技大学青年创新奖资助课题(20180104)联系方式：申诗雨，。通信作者：杨晓虎，。第35卷第10期强激光与粒子束Vol.35，No.102023年10月HIGHPOWERLASERANDPARTICLEBEAMSOct.，2023104005-1性，如成丝不稳定性和水龙管不稳定性，保证电子束稳定的长距离传输9。然而，人们对高能质子的研究较少，而高能质子的传输对空间武器和空间垃圾清理等的研究有重要作用，故有必要对高能质子的空间传输精确控制进行研究。高能质子在空间中的传输与地磁场的强度及方向、本身能量等因素密切相关，已有的研究并未充分考虑地磁场对于高能质子的偏转，磁场中运动粒子的同步辐射，高能质子的相对论效应对传输过程的影响，以及高能质子传输的准确控制等。本文细致研究了不同能量的高能质子在空间传输的精确控制问题，综合考虑地磁场、相对论以及同步辐射等影响，建立了高能质子的单粒子运动模型。给定初始能量，利用贝叶斯优化方法，实现了高能质子在空间中从起始点到目标点的精确传输。1高能质子在地磁场中的传输首先，我们建立地磁场的笛卡尔坐标系，认为地磁北极与地理北极重合，以地球球心为圆心，以从球心向地磁北极的连线为 z 轴方向，向赤道与零度经线交点的连线为 x 轴方向建立右手笛卡尔坐标系。高能质子在空间中传输时受到地磁场的作用，由于受到洛伦兹力而产生速度的偏转，在某种情况下地球可认为是一个巨大的磁偶极子。所以，可以用偶极场来近似地磁场，其模型为10Bdip(r)=04r33(M r)r M（1）r=x x+y y+z z r=r/rMM=M zx=ReRey=z=0B0=3.07105T0M/4=B0R3e式中：为矢径，为单位矢量，为地球的磁矩。对于地球，反平行于 z 轴。在地球赤道处（为地球半径），地磁场强度约为，得。于是在笛卡尔坐标中，磁场的表达式为Bdip=B0R3er53xz x+3yz y+(2z2x2y2)z（2）R2=0.973 9选取零度经线与赤道交点处高度为 0km 至 1000km每隔 10km 取一个值的数据与 IGRF 模型数据比较，如图 1所示，可以看出理论模型计算所得数据与 IGRF 数据库11数据拟合效果较好，拟合优度计算得，验证较好，可以代替实际地磁场使用。m0近地空间中，带电量为 q、质量为、速度为 v 的质子，在电场强度为 E、磁感应强度为 B 的场中运动时，其所受的磁场作用力由洛伦兹力公式给出，即dm0vdt=qE(r)+qvB(r)（3）1/10当粒子能量大于本身静止能量的时，应当考虑相对论对于粒子速度、质量等的影响，相对论对于粒子的影响根据洛伦兹变换公式给出，即dmvdt=qE(r)+qvB(r)（4）=(1v2/c2)m=m0/(1v2/c2)式中：为相对论因子，。给定质子能量为 100MeV，出射位置为高度 500km、经度 0、北纬 30处，出射角度为(60,0)，即出射径向角为 60，出射轴向角为 0，运动时间为 0.01s。是否考虑相对论的高能质子在地磁场中运动轨迹图如图 2 所示。可以发现，质子在地磁场中呈回旋运动，总体轨迹呈沿磁场方向螺旋前进的螺旋线，而由于相对论的作用，质子质量发生变化，变化约为 11%，最终作用于洛伦兹力偏转半径时相差量级约为 100m，虽然运动时间很短，但由于质子速度接近光速，导致实际运动距离较长，最终导致差距较大。如图 2 所示，在运动时间为 0.01s 时考虑相对论效应和未考虑相对论效应对于质子运动有较大影响，经过计算得同一时刻下两种情况间最大距离可达 20km 以上。2同步辐射对高能质子传输的影响磁场中的高能粒子运动时受到洛伦兹力的作用产生垂直于运动方向的加速度，从而发生同步辐射，粒子在发02004006008001 000h/km2.02.53.03.5B/(105 T)theoretical model dataIGRF database dataFig.1MagneticinductionintensitydiagramoftheoreticalmodelandIGRFdatabase图1理论模型与 IGRF 数据库磁感应强度图强激光与粒子束104005-2射同步辐射后损失能量，影响后期洛伦兹偏转半径和运动轨迹，故需予以考虑。粒子在磁场中发射的同步辐射功率由下式确定12P(t)=e2 v260c31(1v2/c2)2（5）给定质子能量为 100MeV，出射位置为高度 500km、经度 0，北纬 30处，出射方向为(60,0)，即出射径向角为 60，出射轴向角为 0，运动时间为 0.01s。是否考虑同步辐射的高能质子在地磁场中运动轨迹图如图 3 所示。可以看出由于地磁场的磁镜效应，质子出现弹跳运动；同步辐射对质子运动轨迹的影响与相对论的影响相似，是否考虑同步辐射对于质子的总体运动轨迹影响不大，其对于质子的影响在运动初期并不明显，但由于质子的速度接近光速，仅 0.01s 的时间两者之间的距离已经相差了 80km 以上，相差了将近一个螺旋，可以说明同步辐射的影响会随着运动距离的增加而逐渐增加，予以考虑。3空间中的单粒子运动模型建立及应用3.1空间中的单粒子运动模型的建立由于粒子在运动过程中同时受到地磁场偏转、相对论效应以及同步辐射的影响，因此在综合计算中将时间离散化后需要在一个时间步长里同时进行三个效果的计算。由于时间步长取得很小，在某一个时间步长内的计算可认为粒子的状态不发生变化。首先使用上一时刻的最终速度计算由相对论效应引发的质量变化；之后由于地磁场的作用发生洛伦兹偏转及同步辐射：受到洛伦兹偏转会影响粒子的加速度，受到同步辐射会影响粒子能量进而影响速度，且偏转和辐射应是同步发生，故选择上一时刻的速度及本时刻计算所得的质量同时计算偏转和辐射，在计算的过程中需将粒子的运动在三维进行展开来分别计算每个方向的速度及加速度。再使用通过偏转计算得到的加速度以及通过辐射计算得到的速度来计算并保存本时刻的位移、计算并保存本时刻的最终速度并作为下一时刻的初始速度。3.2空间中的单粒子运动模型的应用X0Xtdt=108s建立笛卡尔坐标系，并通过坐标变换将给定坐标中的经纬度换算为笛卡尔坐标系下的坐标，给定初始点的坐标、目标点的坐标、质子的初始能量及出射角度，将质子的运动在时间上进行离散化，以为一个时间步长，共走 N 步，应用上述空间中的单粒子运动模型，计算每个时刻粒子的坐标及速度并保存，最后得到粒子运动坐标（Xp)及运动速度（vp)的 3*N 维矩阵。XpXtN通过计算中每个坐标点与目标点的距离，找到其中的最小值，得到此出射速度方向下与目标点的最小距离；通过使用每时刻粒子的速度即可计算得到每时刻粒子的能量，与初始能量相比即可得到每时刻的能量损耗；再找到距离最小值的坐标后，其对应的步数即对应粒子运动的时间。4基于贝叶斯优化方法对高能质子传输的优化控制高能质子在空间传输中由于受到地磁场影响，其运动轨迹比较复杂，难以实现准确控制。因此，我们提出利用500y/km3 4203 4403 460503 480 x/km3 5005 4003 520z/km5 5003 5405 6003 5605 7003 5805 8003 6003 6205 9006 000without relativistic effectwith relativistic effectstarting pointfinal point without relativistic effectfinal point with relativistic effectFig.2Protontrajectorywithandwithoutrelativisticeffect图2考虑相对论效应和未考虑相对论效应的运动轨迹图0200y/km5 0006 0007 0002 0008 000 x/km4009 00010 0000z/km2 0004 000without synchrotron radiationwith synchrotron radiationstarting pointFig.3Protontrajectorywithandwithoutsynchrotronradiation图3考虑同步辐射和未考虑同步辐射的质子运动轨迹图申诗雨，等：基于贝叶斯优化的高能质子空间传输精确控制104005-3P(A|B)=P(B|A)*P(A)/P(B)P(A|B)P(B|A)P(A)P(B)贝叶斯优化方法对高能质子空间传输进行精确控制。贝叶斯优化是一种非常有效的全局优化算法，可以快速高效的找到相对最优的输入参数。目标函数通常是一个黑盒函数，它的精确表达式或导数是未知。高斯过程是贝叶斯优化的基础13。一个高斯过程由其均值函数和协方差函数唯一定义14。高斯过程假设给定参数为正态分布，并将具有不确定性的未知函数建模为替代函数。通过应用贝叶斯公式，可以从先验分布得到后验分布。表示为后验分布，代表需要优化的分布。同时，代表目标函数，是未知函数的先验分布，代表证据或边际可能数据15。经过多次迭代，基于高斯过程的贝叶斯优化器会以一定的置信水平逼近函数，并逐渐增强置信水平，收敛到最优解。y=f(X)以到目标点距离和质子出射角度的优化为例。将目标点坐标，出射角度及其对应的路径中每个点的坐标之间的关系建模为，其中 X 代表相关参数X=X1，Xn（6）y=(XtXp)2设为目标点与路径点之间的距离，可以将问题转化为寻找出射角度对应到目标点最小距离的优化问题 A，即X=argminXAf(X)（7）X0Xf(X)贝叶斯优化需要初始采样来获得数据集。给定均值函数和核函数，高斯过程回归然后可以通过使用均值向量和协方差矩阵来提供数据集的简洁表示。X0(X)使用均值和核函数获得初始函数 f1，但由于有限的数据和高度不确定性，它可能不是最佳的。为了改进估计，除了先前的数据点，还引入了额外的观察数据点。这降低了不确定性，并提供了更好的设计变量。可以使用各种方法来选择下一个数据点，这可以使用高斯过程回归中的采集函数来确定16。在这项工作中，采用置信界限标准方法，以避免局部最优，并根据当前估计找到新的最优解。贝叶斯优化流程如图 4 所示。100 MeV100 km现给定质子能量为，初始位置为距离地球表面高度 500km、经度 0、北纬 30处，目标位置为距离地球表面高度 500km、经度 0、北纬 30.9处，两者距离为，将初始条件代入上述的综合模型中，使用贝叶斯优化后得到结果如图 5 所示。图 5 给出了此情况下运动过程中与目标点的最小距离与出射方向的关系图，可以发现最小距离与出射方向之间不存在突变，其变化较平滑，最小值的区间比较明显，据此可以缩小出射方向的范围，进行下一步优化，优化得到三种最优出射情况。(v,v)图 6（a）给出了三个最小范围的最优出射方向，其分别为(85.7,120.8)、(99.5,120.7)、(63.3,98.6)，能initial point,targetcoordinates and emergentdirectionselect initialsamplecalculate theminimum distancefrom the targetIf the minimum distanceless than 10 moutput minimumdistance and exitdirectionGaussionregressionfitting the valuecorresponding to theminimum acquisition functionselect newsampleFig.4Bayesianoptimizationflowchart图4贝叶斯优化流程图20406360min distance/km8010024v/radv/rad210observed pointmodel meanbest pointFig.5Relationshipbetweentheminimumdistanceandtheexitvelocitydirection图5最小距离与出射速度方向关系图强激光与粒子束104005-4MeVMeVMeV1.078105s1.252105s2.895105s1.078105s(v,v)量损耗约为 4.99、4.99，4.98；运动时间分别为、；距离目标点最小距离分别约为 5.55m、1.15m、6.56m，均可满足与目标点的最小距离小于 10m 的限制条件。结合同步辐射的能量损耗，粒子运动时间以及与目标点的最小距离三种评价标准，我们可以分别给出能量损耗最小的出射方向、运动时间最短的出射方向以及与目标点有最小距离的出射方向：结果表明第三个方向为能量损耗最小的方向，能量仅损耗 4.98MeV，即图中黄色的曲线；第一个方向为运动时间最短的出射方向，耗时约，即图中蓝色的曲线；第二个方向为与目标点距离最近的出射方向，最小距离为 1.15m，即图中红色的曲线。对于三个方向的综合评价，首先，我们在研究本方法的应用后认为能否击中目标对效果影响最大，而与目标点的距离越小代表着击中目标的概率越大，故认为最小距离应优先考虑；其次，在击中目标后粒子的能量代表了毁伤能力，故粒子损耗能量应当第二考虑；最后，由于粒子的运动速度基本接近光速，且通过实际结果发现运动绝对时间非常小，故认为运动时间应最后考虑。对三种评价标准进行分权，其中对最小距离赋权 70%，对能量损耗赋权 20%，对运动时间赋权 10%。以最优值为标准值计算相对偏差来综合评价三种出射方向的优劣，最终得第二个出射方向为最优出射方向，其加权平均相对偏差仅为 3.24%。图 6（b）给出了三个无法满足限制条件的出射方向的粒子运动轨迹图，即在这三种出射方向下，粒子的运动路径与目标点的最小距离大于 10m，其分别为(160.5,120.4)、(91.7,114.6)、(63.0,103.2)，距离目标点最小值分别约为 81711m、86844m、25786m。图 6（c）给出了初始点与目标点相同高度下最小距离与出射方向的关系图，可以发现小于 10m 最小距离的范围呈类椭圆形堆叠分布，且满足条件的角度范围只有 0.01量级，说明对于出射方向的掌控需要十分精细。图 6（d）给出了初始点较目标点高 10km 时最小距离与出射方向的关系图，在出射高度上升后仍呈类椭圆形堆叠分布，出射的角度相较相同高度时增加，但类椭圆形数目减小，范围减小，说明在改变目标点与出射点的相对位置时最优出射方向的改变是连续变化的。以第二个最优出射方向为例，可以发现在迭代 5 次之后即可得到预期结果，效率较高，其优化过程中最小距离y/km5 950(a)three situations that can hit the target point3 440505 90003 46050 x/km3 480z/km3 5003 520first casesecond casethird casestarting pointtarget3 4005 950 x/km3 450(b)three situations that can not hit the target point50z/km3 50005 900y/km3 550502.060 462.060 502.060 54v/radv/rad1.737 81.738 01.738 21.738 41.738 61.738 8v/radv/rad(c)the diagram of minimum distanceand exit angle at the same height123456789102.131 252.131 302.131 351.930 31.930 41.930 51.930 61.930 71.930 8(d)the diagram of minimum distanceand exit angle at a higher height2345678910first casesecond casethird casestarting pointtargetmin distance/mmin distance/mFig.6Schematicdiagramofoptimizationresults图6优化结果示意图申诗雨，等：基于贝叶斯优化的高能质子空间传输精确控制104005-5与迭代次数的变化关系如图 7 所示。现以目标点为坐标原点建立与原坐标系方向相同的球坐标系，在以目标点为球心的 50km 半径圆上设置初始点，限制出射方向量级为 0.01得出的最佳出射方向，可以发现当位置径向角取 0和 180时，位置轴向角的取值不会改变粒子的最优出射方向，出射方向与出射位置关系如表 1所示。5结论本文综合考虑地磁场、相对论效应以及同步辐射之后建立了高能质子的单粒子空间运动模型，并通过贝叶斯优化的方式来实现给定能量下高能质子传输的精确控制，获得了能量损耗最少的最优解，即传输的最优路径，给出了满足限制条件的出射角度范围及变化趋势；并通过在目标点附近设置与目标点距离相同的一系列出射点，得到了出射位置改变对出射方向的影响规律。研究结果将为空间粒子束长程传输的相关引用提供参考和借鉴。参考文献：NeubertT,GilchristB,WildermanS,etal.RelativisticelectronbeampropagationintheEarthsatmosphere:modelingresultsJ.GeophysicalResearchLetters,1996,23（9）:1009-1012.1KrauseLH.Theinteractionofrelativisticelectronbeamswiththenear-earthspaceenvironmentD.AnnArbor:UniversityofMichigan,1998.2MironychevPV,BabichLP.Propagationofanelectronbeaminatmosphereataltitudesfrom15to100km:numericalexperimentJ.HighTemperature,2000,38（6）:834-842.3PorazikP,JohnsonJR,KaganovichI,etal.ModificationofthelossconeforenergeticparticlesJ.GeophysicalResearchLetters,2014,41（22）:8107-8113.4WillardJM,JohnsonJR,SnellingJM,etal.Effectoffield-linecurvatureontheionosphericaccessibilityofrelativisticelectronbeamexperimentsJ.5表1不同出射位置的最佳出射方向及最小距离Table1Optimalemissiondirectionandminimumdistanceatdifferentemissionpositionscoordinateradialangle/()coordinateaxialangle/()radialangleofexit/()axialangleofexit/()minimumdistancefromtarget/m00146.34271.762.8450126.45214.084.890014.54259.391.9135050.89155.381.8180033.9892.103.2090146.33271.772.84590101.28274.373.6909056.7986.264.71359011.32284.692.91809033.9992.113.20180146.33271.772.845180128.89334.472.79018091.025.003.213518054.1034.934.518018033.9992.113.20270146.33271.762.845270168.63104.922.49027079.0294.642.6135270123.9296.313.918027033.9992.113.21 266.8935.8211.597.481.1502004006008001 0001 2001 4000123456min distance/miteration timesFig.7Variationofmindistanceanditerationtimes图7最小距离与迭代次数的关系图强激光与粒子束104005-6FrontiersinAstronomyandSpaceSciences,2019,6:56.PowisAT,PorazikP,Greklek-MckeonM,etal.EvolutionofarelativisticelectronbeamfortracingmagnetosphericfieldlinesJ.FrontiersinAstronomyandSpaceSciences,2019,6:69.6郝建红,王希,张芳,等.随移动窗推进的带电粒子束团长程传输模拟分析J.国防科技大学学报,2021,43（5）：168-174.（HaoJianhong,WangXi,ZhangFang,etal.Simulationanalysisoflong-rangepropagationofchargedparticlebeamspropelledbymovingwindowJ.JournalofNationalUniversityofDefenseTechnology,2021,43（5）:168-174）7HaoJianhong,WangXi,ZhangFang,etal.Theinfluenceofmagneticfieldonthebeamqualityofrelativisticelectronbeamlong-rangepropagationinnear-EarthenvironmentJ.PlasmaScienceandTechnology,2021,23:115301.8YaoHaibo,YangXiaohu,ZhangGB,etal.StabletransportofrelativisticelectronbeamsinplasmasJ.JournalofPlasmaPhysics,2022,88:905880105.9钟海坚,陈宗华,赵炳炎.基于 MATLAB 的地磁场中带电粒子运动模拟分析J.大学物理实验,2021,34（1）：83-86.（ZhongHaijian,ChenZonghua,ZhaoBingyan.ThemotionsimulationofchargedparticlesingeomagneticfieldbasedonMATLABJ.PhysicalExperimentofCollege,2021,34（1）:83-86）10MausS.IGRFEB/OL.2023-07-15.https:/ccmc.gsfc.nasa.gov/models/IGRF13/.11李承祖,银燕,赵晶,等.电动力学M.长沙:国防科技大学出版社,2022:305-307.（LiChengzu,YinYan,ZhaoJing,etal.ElectrodynamicsM.Changsha:NationalUniversityofDefenseTechnologyPress,2022:305-307）12SeegerM.GaussianprocessesformachinelearningJ.InternationalJournalofNeuralSystems,2004,14（2）:69-106.13SnoekJ,LarochelleH,AdamsRP.PracticalBayesianoptimizationofmachinelearningalgorithmsC/Proceedingsofthe25thInternationalConferenceonNeuralInformationProcessingSystems.2012:2951-2959.14SrinivasN,KrauseA,KakadeS,etal.Gaussianprocessoptimizationinthebanditsetting:noregretandexperimentaldesignC/Proceedingsofthe27thInternationalConferenceonMachineLearning(ICML).2010:1015-1022.15BrochuE,CoraVM,deFreitasN.AtutorialonBayesianoptimizationofexpensivecostfunctions,withapplicationtoactiveusermodelingandhierarchicalreinforcementlearningDB/OL.arXivpreprintarXiv:1012.2599,2010.16申诗雨，等：基于贝叶斯优化的高能质子空间传输精确控制104005-7