基于PGNAA的化学武器识别
【摘要】 未知化学武器弹药的定性识别在犀护社会安全方面是十分重要的，可指导化学武器的分类处理。瞬发伽马射线中子活化分析（PGNAA）技术利用分析活化产生的伽马射线能谱可以实现对物质中元素的无损，快速检测，在化学武器识别中具有独特的优势。因此，本研究基于PGNAA技术进行了化学武器弹药类型识别装置的设计，同时使用逻辑树判别方法对化学武器样品进行定性分析。首先，基于高纯锗（HPGe）探测器与Cf-252中子源，使用蒙特卡罗MCNP程序对装置结构进行设计优化，主要包括中子源容器尺寸、伽马屏蔽体厚度以及探测器相对位置等。为了最大化样品活化产生的特征伽马射线，需要提高样品位置处的热中子通量，采用聚乙烯作为慢化体，模拟结果显示聚乙烯厚度达到6cm，宽度达到12cm时，样品中热中子通量达到较高水平。为了降低周围材料活化噪声的干扰，选择铅作为屏蔽结构，模拟显示铅屏蔽厚度达到5cm时，可满足屏蔽要求。同时，探测器与样品之间的距离也会影响对伽马射线的探测，最终模拟确定探测器与样品之间的距离为28cm时，特征信号计数最高。根据优化结果搭建测量装置，使用分析纯试剂根据真实化学武器元素含量配制化学武器模拟样品，通过对5种化学武器模拟样品的测量获得伽马能谱。对能谱中的特征峰处理过程中，基于特征峰对元素进行分析，针对计数统计性较好的元素（如H，Cl，S）的特征峰，使用高斯及多项式拟合的方式对特征峰处的高能量康普顿平台进行扣除，获得特征伽马射线的全能峰信息。而对统计性较差的元素特征峰（如N元素的10.829MeV），采用能量区间加和法，对该能量下的全能峰至单逃逸峰之间的计数求和，进而可确定该元素在样品中的存在情况，最后利用建立的逻辑树判别方法根据元素存在信息对样品类型进行判别。实验结果表明，利用该优化的装置可以获得5种模拟样品的能谱，结合能谱分析方法可以得到化学武器模拟样品中的H，Cl，S和N等元素的存在信息，最后使用逻辑树判别方法可以对化学武器样品种类进行判别。关键词化学武器;PGNAA;核素识别;装置优化设计;特征峰拟合分析中图分类号:O582文献标识码:ADOI:10.3964/j.issn.1000-0593（2019）12-3653-06收稿日期:2018-10-13，修订日期:2019-02-06基金项目:国家重点研究发展计划项目（2017YFF0104200），国家自然科学基金项目（11775113），南京航空航天大学研究生创新基地（实验室）开放基金项目（kfjj20180607）和江苏高校优势学科建设工程项目平台资助作者简介:汤亚军，1993年生，南京航空航天大学核科学与技术系硕士研究生e-mail:tangyajun@nuaa.edu.cn*通讯联系人e-mail:heidq@nuaa.edu.cn引言在处置遗弃武器弹药的过程中，常常会因为标识不清而无法分辨弹药的类型，尤其是内部含有危险化学试剂的时候，处理不当会对人员、环境带来严重的后果，所以有必要在处理之前对化学武器类型做出准确判别。常规检测化学武器未知弹药的方法有表面增强拉曼散射技术、声表面波传感技术、离子迁移谱技术[1-5]，这些技术都是通过检测弹药挥发到空气中的微量物质来获取弹药类型信息，对密封性较好的罐体检测效果不理想。利用中子技术检测武器弹药中比较成熟的产品有美国爱达荷实验室开发的PINS（portableiso-topicneutronspectroscopy）系统[6]，自1991年应用以来，已成功检测出数百枚埋藏地下的未知弹药。何铁[7]等设计了一套基于快中子的化学战剂识别装置，采用伴随α粒子技术对能谱本底扣除，大幅度降低噪声的干扰，但探测系统及后端处理系统较为复杂，对探测系统的时间分辨能力要求较高。本研究采用的瞬发伽马射线中子活化分析（promptgam-ma-rayneutronactivationanalysis，PGNAA）技术可以对多元素实时在线检测分析，并且具有非破坏性、高灵敏度等特点，在未知弹药检测方面具有独特的优势[8]。文章设计了一套利用中子作为探针，通过分析中子活化产生的瞬发伽马光谱来对样品核素进行判别的装置，为废置弹药的分类和处理提供依据。1系统设计PGNAA技术是通过检测中子与物质相互作用放出的特征伽马射线来区分元素种类及含量的一种分析方法，为了获得更多的具有元素特征的伽马射线信息，需要对慢化体、样品及探测器位置、必要的屏蔽材料等参数进行优化设计，以便有效降低散射中子、散射光子对测量系统的影响。本研究采用蒙特卡罗MCNP程序，通过中子、伽马射线输运的模拟研究，给出了PGNAA系统的优化设计方案，见图1中所示。下面对模拟设计的相关问题进行简单描述。图1PGNAA系统优化设计方案示意图Fig.1TheschematicdiagramofPGNAAoptimizeddesign1.1慢化体及屏蔽体优化装置使用的是Cf-252裂变中子源，中子源平均能量是2.158MeV，大部分中子能量处于快中子区，为了获得更多的辐射俘获所需要的热中子，需要对慢化体的几何结构进行优化设计。为了增加样品处的热中子通量，采用氢含量较高，结构性能较好的聚乙烯作为慢化体，改变慢化体厚度与宽度，MCNP程序计算结果如图2所示，随着慢化体厚度与宽度的增加，样品中的热中子通量逐渐增大，达到一定厚度与宽度后，样品中热中子通量增加速率变小，考虑慢化效果及装置成本，选择中子慢化体厚度为6cm，宽为12cm。中子与探测器相互作用产生的伽马射线会直接沉积在探测器中，对噪声产生较大的贡献。因此，为降低探测器活化噪声的贡献以及中子损伤，在中子源与探测器之间添加对中子吸收能力较强的含硼聚乙烯作为中子吸收体。然而，添加的中子防护体在与中子反应的过程中，其产生的伽马射线同样为干扰噪声的主要来源之一，为减少此部分对于噪声的贡献，在探测器周围设置伽马屏蔽体。一般而言，在PGNAA装置中，通常使用高原子序数、密度较高的材料对伽马射线进行吸收屏蔽，基于此，选择铅作为伽马屏蔽层材料。同时，根据伽马射线衰减规律，伽马射线衰减强度只与屏蔽层厚度和伽马能量有关，在同样屏蔽层条件下，衰减强度与伽马射线强度呈反比。含硼聚乙烯材料与中子发生反应主要产生H（2.22MeV）以及B（0.478MeV）伽马射线，因此，为有效屏蔽含硼聚乙烯材料中产生的伽马射线，需要以H（2.22MeV）特征伽马射线能量为基准进行屏蔽体厚度的优化选择。随着铅屏蔽厚度变化，H特征伽马射线衰减强度结果如图3所示。图2样品热中子计数与慢化体尺寸关系（a）:随厚度变化;（b）:随宽度变化Fig.2Therelationshipbetweenthenormalizedcountandthesizeofmoderator（a）:Thickness;（b）:Width图3H特征信号与屏蔽体厚度关系Fig.3TherelationshipbetweenthecharacteristicpeakofHandthethicknessofshielding根据图中的模拟数据，拟合伽马射线衰减曲线式（1）y=9.28887×10-7+1.84757×10-5e-0.71389x（R2=0.99992）（1）4563光谱学与光谱分析第39卷探测器中的氢特征峰计数随屏蔽体厚度的增加而减少，当铅屏蔽体厚度达到5cm时，计数减小率达到92.73%，可以对周围材料产生的大多数干扰信号进行屏蔽，最终选择5cm的铅作为屏蔽体。1.2样品位置优化装置设计过程中使用的模拟样品为芥子气，其中含有C，H和Cl等元素，样品初始放置在与中子源平齐的位置处，与慢化体呈45°角放置，如图4（a）所示。沿y轴方向移动待测样品，改变样品与探测器之间的距离，探测器中Cl元素特征伽马计数变化如图4（b）所示。图4样品位置优化（a）:样品移动位置示意图;（b）:Cl特征峰与样品位置关系Fig.4Theoptimizationofsampleposition（a）:Movingthesample;（b）:Therelationshipbetweenthecharac-teristicpeakofClandthesampleposition样品位置发生改变时，样品活化产生的特征伽马射线被探测器探测到的强度也会发生改变。计算结果如图4（b）所示，在距探测器28cm处计数最高，此时探测器可获得更多的特征信号计数，所以样品的最终位置为与探测器表面距离28cm处。2实验部分按照优化设计方案搭建了测量装置，装置由Cf-252中子源、钨屏蔽体、慢化体、中子防护体及铅屏蔽体组成。测量装置平台如图5所示。模拟的样品放置在不锈钢容器内，容器高30cm，半径4.5cm，不锈钢外壳厚度为0.3cm，样品重3kg。使用AR级化学药品（包括氯化钠、聚乙烯、升华硫、三聚氰胺、二氧化硅等）按照特定元素配比制作了化学武器仿真样品，各组样品元素含量如表1所示。图5实验测量装置图Fig.5Experimentalsetup表1化学武器仿真样品元素配比Table1Elementalcompositionofsimulantchemicalweapon编号名称模拟样品元素配比/%CHONSClSiNa样品1HD4.420.74//2045.16/29.68样品2L13.052.18///51.15/33.62样品3CK10.771.8/25.13037.59/24.71样品4VX15.632.6/36.4745.3///样品5沙子//53.33///46.67/实际测量过程中，仅以样品编号代表各化学武器模拟样品，测量样品之前，使用钴-60和铯-137标准源对探测器进行能量响应的标定[9]。探测系统使用的是ORTEC公司生产的效率为30%的N型高纯锗（highpuritygermanium，HPGe）探测器，配合电制冷系统，以及DSPEC-50型号的数字化多道谱仪，中子源使用活度为6×106n·s-1的Cf-252中子源，每个样品测量活时间为36000s。3结果与讨论3.1元素分析对样品的种类识别作要依靠对测量能谱中特征峰计数的分析，中子与样品相互作用产生特征伽马射线能量及截面如表2所示。表2主要元素特征峰及截面Table2Peaksandcrosssectionsofkeyelements元素特征峰能量/keV截面/barnsCl19516.3319594.161116.59H22230.3326S8410.34717200.17522300.56N108290.0113测量结果的部分能谱图如图6所示。5563第12期光谱学与光谱分析图6样品测量能谱图（a）:5种样品能谱图;（b）:部分元素特征峰;（c）:N元素特征峰Fig.6Gamma-rayenergyspectra（a）:Spectraof5samples;（b）:Characteristicpeaksofpartialelements;（c）:CharacteristicpeaksofN图6中可以比较清晰地看到氯元素的辐射俘获峰（1951和1959keV）和硫的特征峰（2230keV），由此可以判断在样品中含有Cl和S元素。在图6（c）中可以看到，N元素的辐射俘获峰（10.83MeV）和它的单逃逸峰（10.31MeV），与Cl和H元素相比，N元素的特征峰较弱，主要原因是作用截面较小，并且探测器在较高能量区域内的探测效率较低，但该特征峰仍然可见，由此可以判断该样品中含有N元素。能谱中的特征峰范围可以认为由两部分组成:全能峰与康普顿平台，当平台计数较高时，会影响对全能峰峰面积的计算，这里采用多项式和高斯拟合的方法对特征峰进行处理，如式（2）f（x）=g（x）+G1（x）+G2（x）+…+Gn（x）（2）式中f（x）是谱函数，为能谱中特征峰范围区间;g（x）是平台函数，可以使用一次多项式、二次多项式或三次多项式拟合，通常采用二次多项式;G（x）是峰形函数，在满足统计涨落的情况下呈正态分布。对能谱中Cl特征峰拟合如图7所示，原始能谱扣除康普顿平台之后呈高斯正态分布，对所有能谱中的Cl和S元素特征峰作拟合扣除康普顿平台后全能峰计数如表3所示。图7Cl特征峰原始能谱与拟合全能峰对比Fig.7ThecomparisonbetweenoriginalspectrumandfittedspectrumofCl表3Cl和S元素拟合扣除本底后全能峰计数Table3ThefullpeakcountsforClandSelementsafterdeductingthebackground样品编号ClS1.951MeV1.959MeV2.230MeV样品165124349569样品2728148480样品3733645720样品400541样品5000对于N元素来说，10.829MeV特征峰能量附近中子相互作用截面小、伽马射线探测效率低导致该能量范围内的计数较较低，统计性较差，无法使用拟合方法扣除康普顿平台。在伽马能谱的高能段部分，只存在少量的电子学噪声，自然本底以及结构噪声的干扰较小，因此可以在全能峰（10.829MeV）至其单逃逸峰（10.31MeV）之间的计数进行加和，利用该能量范围内的计数波动对N元素定性分析（如表4所示），计算结果如表4所示。可以明显看出，样品3和样品4在该能量区域具有较高的计数变化，通过该结果判断样品3和样品4中含有一定量的N元素，与实际相符。表410.31~10.829MeV净计数加和Table4Netcountadditionfrom10.31to10.829MeV样品编号样品1样品2样品3样品4样品5加和计数-23-145122-23.2逻辑树建立及样品判别根据样品中元素的实际存在情况建立样品判别逻辑树，结合实验分析得到的元素存在信息如图8所示（其中0代表无该元素，1代表有该元素），样品1，2和3中存在Cl元素，6563光谱学与光谱分析第39卷样品1和4中存在S元素，样品3和4中含有N元素。已知样品中的元素存在信息，可以利用图8所示的逻辑树判别样品种类，实现判别功能。图8逻辑树判别流程图Fig.8Thelogicaltreeforidentifyingsamples4结论设计了一种基于PGNAA技术对核素进行判别的装置，使用该装置可以对部分化学武器样品种类进行判别。设计过程中主要使用MCNP程序对装置各部分进行了模拟优化，使装置中的样品获得更高的热中子通量，对探测器进行屏蔽防护，尽可能降低周围材料的噪声信号。根据优化结果最终确定慢化体厚度为6cm，宽度为12cm，铅屏蔽体厚度5cm，样品与探测器表面距离为28cm。利用该装置对5种不同的样品进行测量，通过高斯拟合扣除康普顿平台的方法，使用逻辑树及特征峰信号有无对样品种类进行判断，结果表明可以使用该方法识别样品。研究中初步选择的测量时间为36000s，而在实际需求中，对测量时间的要求较高，往往需要在几分钟内判断测量结果，在满足计数统计性的前提下减少测量时间，这将是后续的研究。此外，为了扩大样品的测量范围，可以配置更多种类的样品建立样品库，同时优化谱分析方法，对待测样品元素的检测阈值做研究分析。References[1]LIUPeng，WANGXiao-li，JIANGJin-jun（刘鹏，王小力，姜进军）.ChineseJournalofSpectroscopyLaboratory（光谱实验室），2011，28（2）:955.[2]XuJY，WangJ，KongLT，etal.JournalofRamanSpectroscopy，2011，42（9）:1728.[3]ZHANGPing，MAJin-yi，CHENChuan-zhi，etal（张萍，马晋毅，陈传治，等）.Piezoelectrics&Acoustooptics（压电与声光），2015，37（5）:734.[4]ZHAOJian-jun，LIUWei-wei，PANYong，etal（赵建军，刘卫卫，潘勇，等）.ChemicalSensors（化学传感器），2011，31（1）:47.[5]TrnesJA.InternationalJournalforIonMobilitySpectrometry，2016，19（2-3）:1.[6]CaffreyAJ，ColeJD，GehrkeRI，etal.IEEETransactionsonNuclearScience，1992，39（5）:1422.[7]HETie，ZHENGPu，ANLi（何铁，郑普，安力）.HighPowerLaserandParticleBeams（强激光与粒子束），2016，28（5）:90.[8]JIAWen-bao，ZHANGYan，HEIDa-qian（贾文宝，张焱，黑大千）.SpectroscopyandSpectralAnalysis（光谱学与光谱分析）.2014，34（11）:3123.[9]CanCheng.AppliedRadiationandIsotopes，2018，145:1.ResearchontheIdentificationofChemicalWeaponBasedonPGNAATechnologyTANGYa-jun1，JIAWen-bao1，2，HEIDa-qian1，2*，LIJia-tong1，CHENGCan1，CAIPing-kun1，SUNAi-yun1，ZHAODong1，HUQiang11.InstituteofNuclearAnalyticalTechnology，NanjingUniversityofAeronauticsandAstronautics，Nanjing211106，China2.CollaborativeInnovationCenterofRadiationMedicineofJiangsuHigherEducationInstitutions，Suzhou215000，ChinaAbstractIdentifyingtheunknownchemicalweaponsisanimportantworkformaintainingsocialsecurity，anditcanguidethechemicalweaponsdestruction.Promptgamma-rayneutronactivationanalysis（PGNAA）technologyhastheadvantagesofbeingnon-destructiveandrapid.Inthisresearch，thedevicewasdesignedforidentifyingthechemicalweaponsbasedonPGNAAtechnology，andlogic-tree-baseddiscriminationmethodwasemployedtoconductqualitativeanalysisonsamples.Firstly，withthehighpuritygermanium（HPGe）detectorandCf-252neutronsourceasthecoreinstruments，thestructuresofdevicewereoptimizedusingtheMonteCarloMCNPcode，includingneutronsourcemoderator，thethicknessoftheshieldingbody，andthe7563第12期光谱学与光谱分析relativepositionofthedetector.Tomaximizethecharacteristicgammaraygeneratedbysampleactivation，itisnecessarytoin-creasethethermalneutronfluxinthesample.Inthisresearch，polyethyleneisusedasmoderatortoincreasethescatteringofneutronsbeforethesample，sothatmoreneutronsarethermalized.Thesimulationresultsshowedthatthethermalneutronfluxinthesamplereachesahighlevelwhenthepolyethylenehasathicknessof6cmandawidthof12cm.Inordertoreducethein-terferenceofthesurroundingmaterialactivationnoise，leadisselectedastheshieldingstructure.Andthesimulationshowedthatitcanmeettheshieldingrequirementwhentheleadshieldingthicknessreaches5cm.Atthesametime，thedistancebe-tweenthedetectorandthesamplealsoaffectsthedetectionofgammarays.Thefinalsimulationdeterminesthatthedistancebe-tweenthedetectorandthesampleis28cm，andthecharacteristicsignalcountisthehighest.Accordingtotheoptimizationre-sults，theexperimentaldevicewassetup.Andthesimulatedsamplesofchemicalweaponwerepreparedaccordingtotheactualelementsbyusingtheanalyticalpurereagent，andthegammaspectrumwasobtainedbymeasuringthefivesamples.Intheprocessofanalyzingthecharacteristicpeaksinthespectrum，theelementswereanalyzedbasedonthecharacteristicpeaksoftheelements.Theelementswithwellstatisticwereanalyzedbyusinggaussandpolynomialfitting（suchasH，Cl，S）.Thehigh-energyComptonplatformatthecharacteristicpeakwasdeductedtoobtainthefullpeakinformation.Fortheelementcharacteris-ticpeakwithpoorstatistic（suchas10.829MeVoftheNelement），theenergyintervalsummationmethodwasusedtosumthecountsbetweenthefullpeaksandthesingleescapeenergypeaks，sowecangettheinformationoftheelementsinthesample.Finally，thelogic-tree-baseddiscriminationmethodwasusedforsampleidentification.Theanalysisresultsshowedthatthein-formationofH，Cl，S，Nandotherelementsinthesimulatedsamplesofchemicalweaponscanbeobtainedbyusingthespec-trumfittingmethod.Andthetypeofsimulatedsamplesofchemicalweaponcanbeidentifiedbycombininglogic-tree-baseddis-criminationmethod.KeywordsChemicalweapon;PGNAA;Nuclideidentification;Deviceoptimization;Characteristicpeakfittinganalysis（ReceivedOct.13，2018;acceptedFeb.6，2019）*Correspondingauthor8563光谱学与光谱分析第39卷
