周梦洁 胡立群 钟国强 黄良胜 李 凯 周瑞杰 洪 兵肖 敏 张瑞雪
1(中国科学院等离子体物理研究所 合肥 230031)
2(中国科学技术大学 合肥 230026)
聚变能是一种安全、清洁,取之不尽用之不竭的新能源,其开发利用研究对于解决人类未来能源危机具有重要意义。通过半个多世纪的探索研究,磁约束核聚变被认为最具有商业应用前景。国内的东方超环(Experimental Advanced Superconducting Tokamak,EAST)是世界首台建成并投入运行的全超导偏滤器位形托卡马克磁约束聚变实验装置,建有大功率射频波、中性束注入辅助加热系统,70余项先进的等离子体诊断系统[1]和完善的在线辐射剂量监测系统。运行初期,EAST进行了少量的氢、氦等离子体测试实验,现阶段主要开展氘等离子体聚变实验研究。通过兆瓦级的大功率辅助加热提高等离子体温度,促使氘核燃料发生如式(1)、式(2)的氘氘聚变反应,并释放大量的快中子。
EAST作为Ⅰ类射线装置,在开展氘等离子体实验时除了产生聚变中子外,还存在X、γ射线。其中高温等离子体中本底热电子发生薄靶轫致辐射以及杂质离子发生线辐射产生的大量低能X射线(keⅤ量级),受装置内、外不锈钢真空室壁的遮挡屏蔽,不会对外围辐射环境产生影响。当等离子体中的电子受到的电场力大于碰撞阻尼力而被持续加速后成为高能逃逸电子,这部分高能逃逸电子轰击装置器壁或孔栏,发生厚靶轫致辐射产生的高能硬X射线(MeⅤ量级),则较易穿透装置器壁,对外部辐射剂量产生贡献。此外,高能逃逸电子或高能硬X射线还会与装置第一壁、杂质等材料发生核相互作用,当能量超过反应阈值时,光中子产生。高能γ射线则主要来自下列三条途径:辅助加热产生的高能离子与杂质发生核反应产生γ射线;中子输运过程中与周围材料发生非弹性散射以及捕获反应产生γ射线;装置主机、周边设备及大厅的结构材料发生中子活化反应产生放射性同位素,其衰变过程中产生γ射线。EAST装置氘运行时中子产额最高可达1.0×1015n·s-1(设计值),对周围辐射剂量水平的高低起主导作用。DD聚变中子(2.45 MeⅤ)的产额正比于聚变核反应率,其与燃料离子的温度(Ti)、密度(ni)、快离子速度分布(ʋf)、能量及粒子约束时间(τe、τp)等参数密切相关。为提高这些与聚变核反应直接相关的等离子体参数,如今大、中型磁约束聚变装置均会开展大量的高功率加热实验研究。近年来,EAST装置通过不断的升级改造,已建成源功率总计达30 MW的射频波与中性束注入辅助加热系统,并开展了各种加热模式的等离子体实验。本文展示了各种加热实验期间中子、γ(X)辐射强度演化的测量结果,并对不同实验条件下的辐射行为和剂量分布特征进行了分析研究。
EAST聚变装置由全超导托卡马克主机以及电源、低温、辅助加热、计算机控制、诊断测量、数据采集等多个系统构成,其主要的工程参数如表1所示[1]。为实现高参数、高约束稳态运行的实验目标,需充足的外部辅助加热手段。目前,EAST装置已建成并投入了多套兆瓦量级的高功率辅助加热系统。其中包括以电子加热为主的电子回旋共振加热(Electron Cyclotron Resonance Heating,ECRH)和低混杂波(Lower Hybrid Waves,LHW)加热,以离子加热为主的中性束注入(Neutral BeamⅠnjection,NBⅠ)加热,既加热离子又加热电子的离子回旋射频(Ⅰon Cyclotron Radio Frequency,ⅠCRF)加热。根据等离子体实验内容的需求,来选择以上各种加热模式及其组合。辅助加热的性能直接影响等离子体的参数品质,进而对聚变核反应率以及反应产物中子产生作用。因此,EAST辅助加热实验既是聚变等离子体物理研究的热点,又是辐射剂量监测研究的重点。
EAST装置运行期间产生的中子、γ(X)剂量主要依靠实时在线辐射剂量监测系统和离线累积剂量计监测系统进行测量。在线辐射剂量监测系统由13个监测点组成,各监测点包括一台BF3中子探测器和一台高压Ar气电离室[2]。通过RS485接口和以太网进行数据通信组网,在装置进行放电实验期间24 h连续运行,具备秒级(最高1 s)时间分辨能力。根据监测点的位置差异,在装置主机厅内的3个监测点采用高量程的区域辐射监测探头,降低聚变等离子体脉冲辐射测量的堆积概率。位于主机厅1.5 m厚屏蔽墙外的10个监测点,采用高灵敏度的环境辐射监测探头[3],可对微量辐射进行有效测量。此外,在各监测点还布置了光释光和固体核径迹剂量片,分别对γ(X)和中子剂量进行累积测量[4]。其测量的结果,可与在线辐射剂量监测系统测得的累积剂量值进行比对校验。为了便于开展聚变辐射剂量与各种加热模式的关联性研究,主要选取主机大厅内部2号监测点的数据进行分析。该监测点距离EAST主机外杜瓦11.3 m,固定于大厅西侧内墙上,离地面约2.2 m。采用的BF3中子探测器,通过设置一定的脉冲阈值,具备良好的γ(X)射线甄别能力,同时避免了目前多数剂量当量率仪灵敏度偏低的缺点[5]。探头经过特别改进设计,在BF3计数管外增加了热中子吸收层,通过调节吸收层的开孔数和面积来调整探测器的灵敏度,使其满足NBⅠ加热时强中子脉冲辐射剂量监测的需求[2]。选择的球形高压Ar气电离室探测器,不仅具有较高的灵敏度和良好的稳定性,还通过采用脉冲电荷积分技术,使得辐射剂量测量上限极大提高。此外,EAST装置上的中子、γ辐射探测器分别采用标准放射源Am-Be源及镭源进行刻度,检定结果显示仪器的总不确定度均小于15%,满足辐射防护仪器校准的要求。
表1 EAST装置主要工程技术参数Table 1 Main engineering parameters of EAST device
EAST装置运行和辐射响应研究还涉及多项诊断参数,主要包括等离子体电流(Ip)、环电压(Pcvloop)、电子弦平均密度(Ne)、中子辐射强度(Neutron-235U)、高能硬X射线(X-ray NaⅠ)及γ射线强度等。其中,最基本的等离子体电流、环电压和电子密度参数将用于放电运行控制,而中子辐射强度、γ(X)射线强度(辐射强度均用计数率表示)等监测结果用于评估等离子体参数品质及辐射行为研究。中子注量监测系统测量的计数率为235U裂变室探测器测量的体积分中子辐射强度,采用重离子加速器打氘靶产生的单能中子(2.4 MeⅤ)对探头的效率进行刻度,其绝对探测效率约为0.17%。高能硬X射线测量系统和γ射线测量系统是研究聚变等离子体中快电子、逃逸电子及快离子动力学行为的重要手段。其中,高能硬X射线和γ射线分别采用NaⅠ(Tl)、LaBr3(Ce)闪烁体探测器进行测量。采用多种同位素放射源(137Cs、152Eu、238Pu-13C)产生的单能γ(X)射线(能量范围0.344~6.13 MeⅤ)对NaⅠ、LaBr3探测器的能量响应进行标定。两种探测器在射线能量为0.2~8 MeⅤ时,均具有良好的能量线性。对于137Cs产生的γ射线,NaⅠ和LaBr3探测器的能量分辨率分别为10%、3.9%。
利用欧姆场加热等离子体是托卡马克装置最基本的一种加热方式,其对于磁约束聚变等离子体的建立具有重要作用[6]。因而首先对该加热模式下的辐射剂量进行监测研究。选取了典型的两炮欧姆放电,其放电波形如图1(a)所示。#77056、#77347炮的平均电子密度分别为1.5×1019m-3和1.99×1019m-3,等离子体电流均为 400 kA。图 1(b)表示两炮欧姆放电过程中,在线辐射剂量监测系统2号监测点测得的辐射剂量,其中:N为中子剂量当量,G为γ剂量当量,Total为总周围剂量当量。两炮(#77056、#77347)等离子体放电产生的中子剂量当量分别为 0.4 μSv 和 0.3 μSv,γ剂量当量分别为1.3 μSv和1.38 μSv,总周围剂量当量分别为1.7 μSv和1.68 μSv。
上述监测结果显示#77347炮的电子密度虽比#77056炮高33%,但#77347炮放电产生的总周围剂量当量却与#77056炮相当。该结果可能是由于#77347炮电子密度大,电子逃逸阈能大(阈能Wc∝Ne/E,E为电场强度),逃逸现象得到一定程度的抑制[7],光中子成分减少,而在欧姆放电中光中子所占中子总产额是主要的[8],因而中子计数减少,通过中子诱导产生的γ射线计数也减少,且与图1(a)中中子、γ(X)辐射强度的时间演化规律相吻合。故在纯欧姆加热模式下,逃逸电子与物质相互作用产生的硬X射线以及硬X射线发生光核反应产生的光中子对总周围剂量当量贡献较大。
图1 欧姆加热时中子和γ(X)辐射强度的时间演化(a)及产生的辐射剂量(b)Fig.1 Temporal evolution of neutron and γ(X)radiation intensity(a)and radiation doses(b)under ohmic heating mode
低杂波驱动是利用等离子体对输入频段的低混杂波进行耦合吸收,进而加热等离子体[9]。在该加热模式下主要研究LHW不同注入时刻对辐射剂量的影响。图2表示等离子体电流为500 kA,平均电子密度为2×1019m-3时注入频率为4.6 GHz的LHW的放电波形图。在电流爬升阶段,电子密度低,电子逃逸阈能低,产生初始逃逸电子,γ(X)射线计数增加;LHW注入等离子体后,环电压迅速下降,电子逃逸阈能增大,逃逸现象得到明显抑制,γ(X)射线计数逐渐减少且随时间的变化趋势相似;LHW稳定注入期间,γ(X)射线计数非常低且波动较小。#77470、#77481和#78915炮LHW注入等离子体对应的时刻分别为1.98 s、3 s和1.49 s,注入前产生的γ射线计数占一炮总γ射线计数分别为68%、88%和43%。故LHW注入时刻对γ剂量影响较大。
图2 LHW加热时中子和γ(X)辐射强度的时间演化Fig.2 Temporal evolution of neutron and γ(X)radiation intensity under lower hybrid waves heating mode
图3 统计了2018年EAST装置夏季实验中,平均电子密度为2×1019m-3和3×1019m-3时,不同LHW注入时刻一炮等离子体放电在2号监测点产生的中子、γ剂量当量。其中:g为γ剂量当量,n为中子剂量当量,g+n为总周围剂量当量。结果表明:LHW注入时刻不同对低电子密度放电下γ剂量当量的影响比高电子密度大。区域1表示中子剂量当量,每一炮放电结束后的中子剂量当量都很低(低于0.1 μSv),几乎不受LHW注入时刻不同的影响。因而在加热期间,离子温度变化不大,由聚变中子造成的剂量当量的影响基本可忽略。区域2表示LHW注入时刻相同(均为3.04 s)、功率不同的两炮。#77065、#77067炮的功率分别为3 051 kW和3207 kW,#77065炮放电产生的γ剂量当量比#77067炮大,说明增加LHW功率,提高了电流驱动能力,降低了等离子体环电压,从而有效抑制了逃逸现象的产生。故LHW加热模式下产生的辐射剂量低于相同条件下的欧姆加热。
图3 LHW不同注入时刻对中子和γ剂量当量的影响Fig.3 Effects of different injection times of lower hybrid waves on neutron and γ dose equivalent
中性束注入加热是利用高能中性粒子注入等离子体内,与带电粒子发生碰撞电离等过程,通过热化加热等离子体[10]。在仅有NBⅠ加热模式下,EAST装置放电过程释放的中子主要来源于聚变中子,而光中子占很少部分,可忽略[11]。
选取了典型的三炮NBⅠ加热放电,其放电波形如图4(a)所示。其中,等离子体电流均为450 kA,#80981、#80982炮的平均电子密度为2.49×1019m-3,#80983炮的平均电子密度为2.63×1019m-3。三炮(#80981、#80982、#80983)等离子体放电在1.1 s时开始注入中性束,对应的实际功率分别为1 260 kW、1 245 kW和1 301 kW,此时中子、γ(X)射线计数迅速增加,且γ辐射强度的时间演化规律与中子类似。中性束注入等离子体后,束靶反应(氘束离子与主等离子体氘靶发生的反应)截面增大,导致聚变中子产额大约增加2~3个数量级,此时γ(X)射线主要来源于中子核反应。图4(b)表示三炮NBⅠ加热放电过程中,在线辐射剂量监测系统2号监测点测得的辐射剂量,其中:N为中子剂量当量,G为γ剂量当量,Total为总周围剂量当量。由于中性束注入过程的持续时间不同(分别为6 s、5.5 s和3 s),三炮等离子体放电产生的中子剂量当量分别为6.30 μSv、5.07 μSv 和 2.30 μSv,γ剂量当量分别为 0.68 μSv、0.61 μSv 和 0.35 μSv,总 周 围 剂 量 当 量 分 别 为6.98 μSv、5.68 μSv和2.65 μSv。表明中子、γ剂量当量不仅与中性束注入功率相关,还受注入持续时间的影响。
图5统计了2018年EAST装置夏季实验中,等离子体电流为450 kA、仅NBⅠ加热时单炮等离子体放电在2号监测点测得的中子、γ剂量当量。通过拟合得到函数Dg(Dn)关系:Dg=0.146 2×Dn0.92,近似为线性,其中:Dn表示中子剂量当量,Dg表示γ剂量当量。由图5可知,单炮等离子体放电产生的γ剂量当量随中子剂量当量的增加近似呈线性增加,表明NBⅠ加热模式下γ(X)射线主要来源于聚变中子。
图4 NBⅠ加热时中子和γ(X)辐射强度的时间演化(a)及产生的辐射剂量(b)Fig.4 Temporal evolution of neutron and γ(X)radiation intensity(a)and radiation doses(b)under neutral beam injection heating mode
图5 NBⅠ加热时单炮等离子体放电产生的中子剂量当量与γ剂量当量的关系Fig.5 Relationship between neutron and γ dose equivalent produced by plasma discharge of single shot under neutral beam injection heating mode
为研究EAST装置上NBⅠ与射频波协同加热时中子、γ射线的多种来源,对2018年EAST装置夏季等离子体放电实验进行了归纳分析,图6统计了等离子体电流为450 kA,NBⅠ与LHW、ECRH及ⅠCRF协同加热时单炮等离子体放电在2号监测点测得的中子、γ剂量当量。由其散点分布与仅NBⅠ加热时的函数Dg(Dn)关系可知,中子剂量当量在2~40 μSv区间时,大部分数据点位于曲线下方,表明所有γ剂量当量由中子诱导γ射线造成,此时中子剂量当量除了聚变中子的贡献外,可能光中子的辐射剂量贡献较大;中子剂量当量在40~80 μSv区间时,大部分数据点位于曲线上方,表明γ剂量当量除了中子诱导γ射线造成外,还可能来源于硬X射线、高能离子的核反应等。在NBⅠ+LHW和NBⅠ+LHW+ECRH协同加热状态下,可能硬X射线对γ剂量当量贡献较大。NBⅠ与ECRH协同加热时,大部分数据点分布在曲线附近,表明ECRH的加入对中子、γ剂量当量影响不大。
图6 NBⅠ与射频波协同加热时单炮等离子体放电产生的中子剂量当量与γ剂量当量的关系Fig.6 Relationship between neutron and γ dose equivalent produced by plasma discharge of single shot under neutral beam and radio frequency waves collaborative heating mode
通过对欧姆、低杂波、中性束注入以及协同加热下中子、γ(X)辐射强度与剂量的实验测量,并结合氘等离子体运行参数进行分析研究。发现逃逸电子是影响欧姆加热模式下光中子和γ(X)剂量的主要因素。在等离子体极低密度或高环电压(破裂)的情况下,逃逸电子轫致辐射会产生大量的高能硬X射线以及光中子。低杂波加热能有效驱动等离子体电流、降低环电压,从而减少逃逸电子产生,进而使得高能硬X射线以及光中子剂量大幅度降低。因此,在相同参数条件下,低杂波加热时的辐射剂量一般低于纯欧姆加热。中性束注入加热时,来自束离子与主等离子体束靶反应产生的DD中子,使得中子总产额提高2~3个数量级。此时中子与周围材料核反应产生的γ射线也大幅增涨,其对于监测的γ(X)剂量贡献占主导。因而,剂量测量结果显示在中性束注入加热时,中子与γ(X)剂量当量存在Dg=0.146 2×Dn0.92的线性关系。中性束与射频波(LHW、ECRH、ⅠCRF)协同加热时,中子、γ(X)射线的产生来源更加多样化。因此,中子与γ(X)剂量当量并不能很好地符合上述的线性关系。通过对各种加热模式下中子、γ(X)剂量进一步的实验监测和分析,可为聚变辐射数据库建立奠定基础,并为人员和周围环境的辐射防护提供基础数据。