CN113536679B - Point source dose rate correction method based on artificial neural network - Google Patents

Point source dose rate correction method based on artificial neural network Download PDF

Info

Publication number
CN113536679B
CN113536679B CN202110816233.9A CN202110816233A CN113536679B CN 113536679 B CN113536679 B CN 113536679B CN 202110816233 A CN202110816233 A CN 202110816233A CN 113536679 B CN113536679 B CN 113536679B
Authority
CN
China
Prior art keywords
data
neural network
artificial neural
dose rate
energy spectrum
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
CN202110816233.9A
Other languages
Chinese (zh)
Other versions
CN113536679A (en
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.)
Nanjing University of Aeronautics and Astronautics
Original Assignee
Nanjing University of Aeronautics and Astronautics
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 Nanjing University of Aeronautics and Astronautics filed Critical Nanjing University of Aeronautics and Astronautics
Priority to CN202110816233.9A priority Critical patent/CN113536679B/en
Publication of CN113536679A publication Critical patent/CN113536679A/en
Application granted granted Critical
Publication of CN113536679B publication Critical patent/CN113536679B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/27Design optimisation, verification or simulation using machine learning, e.g. artificial intelligence, neural networks, support vector machines [SVM] or training a model
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/16Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/18Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/02Neural networks
    • G06N3/04Architecture, e.g. interconnection topology
    • G06N3/044Recurrent networks, e.g. Hopfield networks
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/02Neural networks
    • G06N3/08Learning methods
    • G06N3/084Backpropagation, e.g. using gradient descent

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • Mathematical Physics (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Mathematical Optimization (AREA)
  • Evolutionary Computation (AREA)
  • Mathematical Analysis (AREA)
  • Pure & Applied Mathematics (AREA)
  • Computational Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Computing Systems (AREA)
  • Artificial Intelligence (AREA)
  • Molecular Biology (AREA)
  • Health & Medical Sciences (AREA)
  • Databases & Information Systems (AREA)
  • Computational Linguistics (AREA)
  • Biomedical Technology (AREA)
  • Algebra (AREA)
  • General Health & Medical Sciences (AREA)
  • Biophysics (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Probability & Statistics with Applications (AREA)
  • Operations Research (AREA)
  • Evolutionary Biology (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Medical Informatics (AREA)
  • Computer Hardware Design (AREA)
  • Geometry (AREA)
  • Measurement Of Radiation (AREA)

Abstract

The application discloses a point source dose rate correction method based on an artificial neural network, which comprises the following steps: s1, energy spectrum data and corresponding dose rate values of the unmanned aerial vehicle radioactivity monitoring system at different heights are obtained and are respectively used as input parameters and output parameters; s2, dividing part of energy spectrum data at different heights into training data, dividing the other part of energy spectrum data into test data, constructing an artificial neural network model by utilizing input parameters and output parameters in the training data, and importing the training data into the artificial neural network model for training to obtain a trained artificial neural network model; s3, respectively importing the test data into the trained artificial neural network model to obtain an ideal output result, and comparing errors between the ideal output result and the corresponding test data; and if the error is greater than or equal to the set precision expected value, repeating the steps S2 and S3, and if the error is less than the set precision expected value, training the debugged artificial neural network model to be a point source dose rate correction algorithm.

Description

一种基于人工神经网络的点源剂量率修正方法A point source dose rate correction method based on artificial neural network

技术领域Technical Field

本申请涉及无人机放射性监测技术领域,特别地,涉及一种基于人工神经网络的点源剂量率修正方法。The present application relates to the technical field of radioactivity monitoring by unmanned aerial vehicles, and in particular, to a point source dose rate correction method based on an artificial neural network.

背景技术Background Art

无人机放射性监测是指以固定翼或者旋翼无人机为平台,依托辐射探测器等机载放射性监测载荷设备,进行空中放射性测量,适用于核事故、环境放射性监测、矿产地质勘查等领域。无人机放射性监测系统搭载的常规探测器获取的是核素的能谱信息,在进行辐射评价时,仅能分辨出不同的放射性核素的种类,不能同时测量出核素的空气吸收剂量率,不具备多参数评价手段。若要直观测量剂量率则需要搭载剂量率仪,这大大提高了无人机作业的难度,降低了工作效率,不利于使用无人机平台进行辐射监测工作。同时,在无人机机载系统在作业时,为了直观的观测各环境的辐射水平以及核事故发生时的各点源的位置信息和周围的辐射水平,需要对探测获取的剂量率进行高度修正,将放射源的空气吸收剂量率转换到同一水平高度下,绘制剂量率热力图,即可直观的监测出各个位置的辐射水平。图1示出了实际应用中使用无人机辐射监测系统的示意图。如图1所示,无人机上搭载碘化钠探测器(铊激发),形成无人机辐射监测系统。无人机飞行在不同高度时,探测地面上放射性点源在空气中的剂量率。但此时由于无人机飞行高度会变化,探测到的在同一高度下的剂量率往往不准确,因此需要进行修正。Radioactivity monitoring by drones refers to the use of fixed-wing or rotary-wing drones as platforms, relying on radiation detectors and other airborne radioactivity monitoring payload equipment to conduct aerial radioactivity measurements. It is applicable to the fields of nuclear accidents, environmental radioactivity monitoring, mineral geological exploration, etc. The conventional detectors carried by drone radioactivity monitoring systems obtain the energy spectrum information of nuclides. When conducting radiation evaluation, they can only distinguish the types of different radioactive nuclides, and cannot simultaneously measure the air absorption dose rate of the nuclides, and do not have a multi-parameter evaluation method. If the dose rate is to be measured intuitively, a dose rate meter is required, which greatly increases the difficulty of drone operations, reduces work efficiency, and is not conducive to the use of drone platforms for radiation monitoring. At the same time, when the drone airborne system is operating, in order to intuitively observe the radiation levels of various environments and the location information of each point source and the surrounding radiation levels when a nuclear accident occurs, it is necessary to perform a height correction on the dose rate obtained by the detection, convert the air absorption dose rate of the radioactive source to the same horizontal height, draw a dose rate heat map, and intuitively monitor the radiation level of each location. Figure 1 shows a schematic diagram of the use of drone radiation monitoring system in actual applications. As shown in Figure 1, a sodium iodide detector (thallium excitation) is carried on the drone to form a drone radiation monitoring system. When the drone flies at different altitudes, it detects the dose rate of radioactive point sources on the ground in the air. However, since the flight altitude of the drone changes at this time, the dose rate detected at the same altitude is often inaccurate, so it needs to be corrected.

为了解决上述问题,国内外研究人员做了相关研究工作。针对点源模式的航空辐射监测,当前已经有研究人员使用大型NaI(Tl)探测器(探测系统由18条4.2L的NaI(Tl)晶体组成,总探测体积为75.6L)装备在大型飞机上对地面点源产生的剂量率进行高度修正,目的是寻找点源模式下的高度衰减系数,修正使用的方法是根据IAEA提供的经验公式进行校准。此方法的缺点首先是装载的NaI(Tl)探测器的体积太大,无法适用于小型无人机,导致作业环境受限;其次使用的经验公式的条件较为理想化,未考虑实际应用场景,适用性较差;In order to solve the above problems, researchers at home and abroad have done relevant research work. For aviation radiation monitoring in point source mode, some researchers have used large NaI (Tl) detectors (the detection system consists of 18 4.2L NaI (Tl) crystals, with a total detection volume of 75.6L) to equip large aircraft to perform altitude corrections on the dose rate generated by ground point sources. The purpose is to find the altitude attenuation coefficient in the point source mode. The correction method used is calibrated according to the empirical formula provided by the IAEA. The disadvantages of this method are that the volume of the loaded NaI (Tl) detector is too large to be suitable for small drones, resulting in a limited operating environment; secondly, the conditions of the empirical formula used are relatively idealized, and the actual application scenarios are not considered, so the applicability is poor;

其次也有研究人员采用能谱剂量率转换函数(GE函数),实现探测获得的能谱计数与相对应的剂量率的转换。此方法工作较为繁琐,实现起来步骤繁多,对无人机辐射监测系统时进行实时在线的数据处理造成较大困难;所以,根据无人机放射性监测原理以及实际应用需求,设计一种可以实时在线监测的针对点源的放射性核素剂量率高度修正方案是十分必要的。Secondly, some researchers have used the energy spectrum dose rate conversion function (GE function) to realize the conversion between the energy spectrum counts obtained by detection and the corresponding dose rate. This method is relatively cumbersome and requires many steps to implement, which makes it difficult to perform real-time online data processing in the UAV radiation monitoring system; therefore, according to the UAV radioactivity monitoring principle and actual application requirements, it is very necessary to design a radionuclide dose rate height correction scheme for point sources that can be monitored online in real time.

发明内容Summary of the invention

针对现有技术的不足,本发明提供一种基于人工神经网络的的剂量率修正方法,通过采用基于人工神经网络的预测模型,对探测器获取的能谱中的关键数据进行训练,进而实现预测。In view of the deficiencies in the prior art, the present invention provides a dose rate correction method based on an artificial neural network, which trains key data in the energy spectrum acquired by the detector by adopting a prediction model based on an artificial neural network, thereby achieving prediction.

本发明提供了一种基于人工神经网络的点源剂量率修正方法,包括:S1,获取无人机放射性监测系统在不同高度下的能谱数据对应剂量率值,分别作为输入参数和输出参数;S2,将部分上述不同高度下能谱数据划分为训练数据,另一部分划分为测试数据,利用上述训练数据中的上述输入参数和输出参数构建人工神经网络模型,将上述训练数据导入上述人工神经网络模型中进行训练,获得训练后的人工神经网络模型;S3,将上述测试数据分别导入上述训练后的人工神经网络模型,得到理想输出结果,比较上述理想输出结果与对应的上述测试数据之间的误差;若上述误差大于或等于设定的精度期望值,重复S2、S3,若上述误差小于设定的精度期望值,则训练调试后的人工神经网络模型为点源剂量率修正算法。The present invention provides a point source dose rate correction method based on an artificial neural network, comprising: S1, obtaining dose rate values corresponding to energy spectrum data of an unmanned aerial vehicle radioactivity monitoring system at different altitudes, respectively serving as input parameters and output parameters; S2, dividing part of the energy spectrum data at different altitudes into training data, and the other part into test data, using the input parameters and output parameters in the training data to construct an artificial neural network model, importing the training data into the artificial neural network model for training, and obtaining the trained artificial neural network model; S3, importing the test data into the trained artificial neural network model respectively, obtaining an ideal output result, and comparing the error between the ideal output result and the corresponding test data; if the error is greater than or equal to a set expected accuracy value, repeating S2 and S3; if the error is less than the set expected accuracy value, the trained and debugged artificial neural network model is a point source dose rate correction algorithm.

进一步地,获取无人机放射性监测系统在不同高度下的能谱的方法为:使用蒙特卡洛软件模拟无人机放射性监测系统在不同高度下的能谱。Furthermore, the method for obtaining the energy spectrum of the UAV radioactivity monitoring system at different altitudes is: using Monte Carlo software to simulate the energy spectrum of the UAV radioactivity monitoring system at different altitudes.

进一步地,获取上述不同高度能谱数据中的输入参数、输出参数,并构建人工神经网络模型的方法为:将上述不同高度下的能谱数据进行归一化处理;提取归一化处理后对剂量率沉积起主要贡献的主成分个数作为输入参数;放射源在空气中沉积的空气吸收剂量率为输出参数;利用上述输入参数和上述输出参数计算隐含层;上述人工神经网络模型的结构为,输入层-隐含层-输出层,利用上述人工神经网络模型的结构构建人工神经网络的模型。Furthermore, the method for obtaining the input parameters and output parameters in the above-mentioned energy spectrum data at different heights and constructing an artificial neural network model is as follows: normalizing the energy spectrum data at different heights; extracting the number of principal components that mainly contribute to the dose rate deposition after normalization as input parameters; the air absorption dose rate deposited by the radioactive source in the air is the output parameter; the hidden layer is calculated using the above-mentioned input parameters and the above-mentioned output parameters; the structure of the above-mentioned artificial neural network model is input layer-hidden layer-output layer, and the artificial neural network model is constructed using the structure of the above-mentioned artificial neural network model.

进一步地,提取归一化处理后对剂量率沉积起主要贡献的主成分个数的方法为:提取上述不同高度下的能谱数据中与空气沉积剂量率相关的数据;将上述相关数据进行标准化处理后获得标准化数据;计算上述标准化数据之间的相关系数,并组成相关系数矩阵;计算上述相关系数矩阵的各个特征值;利用上述各个特征值计算每个上述相关数据的贡献率和累计贡献率,当上述累计贡献率大于阈值时,对应的最少的特征值数量,作为起主要贡献的主成分个数,即为输入参数。Furthermore, the method for extracting the number of principal components that make the main contribution to the dose rate deposition after normalization is: extracting data related to the air deposition dose rate in the energy spectrum data at different heights; obtaining standardized data after standardization of the above-mentioned related data; calculating the correlation coefficients between the above-mentioned standardized data and forming a correlation coefficient matrix; calculating each eigenvalue of the above-mentioned correlation coefficient matrix; using the above-mentioned eigenvalues to calculate the contribution rate and cumulative contribution rate of each of the above-mentioned related data, when the above-mentioned cumulative contribution rate is greater than the threshold, the corresponding minimum number of eigenvalues is used as the number of principal components that make the main contribution, which is the input parameter.

进一步地,与空气沉积剂量率相关数据包括:各核素对应的全能峰计数、能量、单逃逸峰、双逃逸峰、湮灭峰和对应的测量时间和高度。Furthermore, the data related to the air deposition dose rate include: total energy peak count, energy, single escape peak, double escape peak, annihilation peak and corresponding measurement time and height corresponding to each nuclide.

进一步地,将上述相关数据进行标准化处理获得上述标准化数据的方法为,利用如下公式:Furthermore, the method for normalizing the above-mentioned related data to obtain the above-mentioned standardized data is to use the following formula:

其中,其中xij表示的是第i个能谱的第j个指标;表示的是xij经过标准化操作之后的数据;其中是表示的是所有上述能谱数据的平均值;为数据的方差,代表数据差异性,sj是代表的两个数据之间的标准差。Among them, x ij represents the jth index of the i-th energy spectrum; It represents the data after xij is standardized; It represents the average value of all the above energy spectrum data; is the variance of the data, representing the data difference, and sj is the standard deviation between the two data represented.

进一步地,计算上述标准化数据之间的相关系数的方法为,利用如下公式:Furthermore, the method for calculating the correlation coefficient between the above standardized data is to use the following formula:

其中,rij为上述能谱数据xi与xj的相关系数;xi为第i个能谱数据的平均值;xj为第j个能谱数据的平均值;n表示的是用来进行计算的能谱的个数,xkj表示的是第k个能谱的第j个指标。Among them, rij is the correlation coefficient of the above energy spectrum data xi and xj ; xi is the average value of the i-th energy spectrum data; xj is the average value of the j-th energy spectrum data; n represents the number of energy spectra used for calculation, and xkj represents the j-th index of the k-th energy spectrum.

进一步地,上述精度期望值为10-3Furthermore, the expected value of the above accuracy is 10 -3 .

进一步地,上述方法还包括:若上述误差小于上述精度期望值,在实际场景中使用无人机放射性监测系统测量相同上述放射源在不同高度下的剂量率,并与上述实际输出结果进行对比,以验证上述人工神经网络模型的质量。Furthermore, the above method also includes: if the above error is less than the above expected value of accuracy, using an unmanned aerial vehicle radioactivity monitoring system in an actual scenario to measure the dose rate of the same above radiation source at different heights, and comparing it with the above actual output result to verify the quality of the above artificial neural network model.

本申请中公开地方法解决现有点源模式下剂量率高度修正过程中剂量验证工作存在花费时间长、人力物力成本高的问题,能提高剂量修正的效率和质量,其结果有助于对验证结果进行分析,减少验证成本。使用BPNN人工神经网络,可以较为有效的提高整个剂量率高度修正的效率,实现核素剂量率高度修正的一体式修正。The method disclosed in this application solves the problem that the dose verification work in the process of dose rate height correction in the existing point source mode takes a long time and has high manpower and material costs, and can improve the efficiency and quality of dose correction. The result is helpful to analyze the verification results and reduce the verification cost. Using the BPNN artificial neural network can effectively improve the efficiency of the entire dose rate height correction and realize the integrated correction of the nuclide dose rate height correction.

本发明的附加方面和优点将在下面的描述中部分给出,部分将从下面的描述中变得明显,或通过本发明的实践了解到。Additional aspects and advantages of the present invention will be given in part in the following description and in part will be obvious from the following description, or will be learned through practice of the present invention.

附图说明BRIEF DESCRIPTION OF THE DRAWINGS

图1为实际应用中使用无人机辐射监测系统的示意图。Figure 1 is a schematic diagram of a UAV radiation monitoring system used in actual applications.

图2为基于神经人工神经网络的点源剂量率修正方法流程图。FIG2 is a flow chart of a point source dose rate correction method based on a neural artificial neural network.

图3为Nal(TI)探测器模型示意图。Figure 3 is a schematic diagram of the Nal (TI) detector model.

具体实施方式DETAILED DESCRIPTION

以下结合附图和实施例,对本发明的具体实施方式进行更加详细的说明,以便能够更好地理解本发明的方案以及其各个方面的优点。然而,以下描述的具体实施方式和实施例仅是说明的目的,而不是对本发明的限制。The following is a more detailed description of the specific embodiments of the present invention in conjunction with the accompanying drawings and examples, so that the scheme of the present invention and its advantages in various aspects can be better understood. However, the specific embodiments and examples described below are only for the purpose of illustration, rather than for limiting the present invention.

下文的公开提供了许多不同的实施方式或例子用来实现本发明的不同结构。为了简化本发明的公开,下文中对特定例子的部件和设置进行描述。当然,它们仅仅为示例,并且目的不在于限制本发明。此外,本发明可以在不同例子中重复参考数字和/或参考字母,这种重复是为了简化和清楚的目的,其本身不指示所讨论各种实施方式和/或设置之间的关系。此外,本发明提供了的各种特定的工艺和材料的例子,但是本领域普通技术人员可以意识到其他工艺的应用和/或其他材料的使用。The disclosure below provides many different embodiments or examples to realize different structures of the present invention. In order to simplify the disclosure of the present invention, the parts and settings of specific examples are described below. Of course, they are only examples, and the purpose is not to limit the present invention. In addition, the present invention can repeat reference numbers and/or reference letters in different examples, and this repetition is for the purpose of simplicity and clarity, which itself does not indicate the relationship between the various embodiments and/or settings discussed. In addition, the present invention provides various specific examples of processes and materials, but those of ordinary skill in the art can be aware of the application of other processes and/or the use of other materials.

首先,需解释下BPNN(反向传播神经网络算法)的由来。人工神经网络(ArtificialNeural Network,即ANN),是20世纪80年代以来人工智能领域兴起的研究热点。最近十多年来,人工神经网络的研究工作不断深入,已经取得了很大的进展,而反向传播算法(Back-Propagation)作为最早且最普遍使用的多层感知机网络训练方法,早在70、80年代就被提出来,并一直沿用至今。以误差反向传播算法为核心的反向传播人工神经网络(BPNN)也同样应用广泛。First of all, we need to explain the origin of BPNN (back propagation neural network algorithm). Artificial Neural Network (ANN) is a research hotspot in the field of artificial intelligence since the 1980s. In the past decade, the research on artificial neural networks has continued to deepen and has made great progress. The back-propagation algorithm (Back-Propagation) as the earliest and most commonly used multi-layer perceptron network training method was proposed as early as the 1970s and 1980s and has been used ever since. The back-propagation artificial neural network (BPNN) with the error back-propagation algorithm as the core is also widely used.

人工神经网络(BP)模型拓扑结构包括输入层、隐含层以及输出层,输入层的神经元由样本属性的维度决定,输出层的神经元个数由样本分类个数决定,隐含层则可以依据实际情况自己决定。从本质上讲,BP算法就是以网络误差平方目标函数、采用梯度下降法来计算目标函数的最小值。基本BP算法包括信号的前向传播和误差的反向传播两个过程。BPNN(反向传播神经网络)具有诸多优点,比如非线性映射能力强、高度的自学习和自适应能力、容错能力强等。使用BPNN,可以较为有效的提高整个剂量率高度修正的效率,实现核素剂量率高度修正的一体式修正。The topological structure of the artificial neural network (BP) model includes an input layer, a hidden layer, and an output layer. The neurons in the input layer are determined by the dimension of the sample attributes, the number of neurons in the output layer is determined by the number of sample classifications, and the number of hidden layers can be determined by the actual situation. In essence, the BP algorithm uses the squared network error objective function and the gradient descent method to calculate the minimum value of the objective function. The basic BP algorithm includes two processes: forward propagation of the signal and back propagation of the error. BPNN (back propagation neural network) has many advantages, such as strong nonlinear mapping capabilities, high self-learning and adaptive capabilities, and strong fault tolerance. Using BPNN can effectively improve the efficiency of the entire dose rate height correction and realize the integrated correction of the nuclide dose rate height correction.

此外,还需简单介绍下本申请中利用的MCNP模拟软件。在MCNP软件是一种利用蒙特卡洛模拟程方法的软件。核技术的学科特点导致很多情况下无法进行实际的测量或实验工作,此时计算机模拟计算就显示出了独特的优势。同时对于结构和反应机制复杂的核反应问题,一般的数值方法难以求解,而MC(Monte Carlo)方法能够准确地模拟实际中的物理过程,解决以往数值方法难以解决的问题,因此该方法广泛地应用于涉核领域的研究中。上世纪40年代中期随着原子能事业的发展,MC方法逐渐成长起来,其基本思想是基于随机数选择的统计抽样方法,由于传统的经验方法不能逼近真实的物理过程,难以得到满意的结果,在解决粒子输运问题方面,蒙特卡罗方法有着显著的优势,随后人们逐渐开发了众多的MC模拟程序,目前主流的有EGS、MCNP、GEANT4等,其中MCNP是由美国洛斯阿拉莫斯国家实验室开发,主要用于三维几何结构中光子、电子、中子等在物质中输运问题的蒙特卡罗模拟程序,其适用的光子能量范围在1E-3MeV—1E5MeV,电子能量范围在1E-3MeV—1E3MeV,中子能量范围1E-11MeV—20MeV。MCNP的程序功能较全,各种物质反应截面数据丰富,减小方差的方法多样,通用性强且使用简单。In addition, the MCNP simulation software used in this application needs to be briefly introduced. MCNP software is a software that utilizes the Monte Carlo simulation method. The subject characteristics of nuclear technology lead to the inability to carry out actual measurement or experimental work in many cases, and computer simulation calculations have just shown unique advantages. At the same time, for nuclear reaction problems with complex structures and reaction mechanisms, general numerical methods are difficult to solve, and the MC (Monte Carlo) method can accurately simulate the physical process in practice and solve the problems that were difficult to solve with previous numerical methods. Therefore, this method is widely used in research related to nuclear fields. In the mid-1940s, with the development of atomic energy, the MC method gradually grew up. Its basic idea is a statistical sampling method based on random number selection. Since the traditional empirical method cannot approximate the real physical process and it is difficult to obtain satisfactory results, the Monte Carlo method has significant advantages in solving particle transport problems. Subsequently, people gradually developed a large number of MC simulation programs. The mainstream ones are EGS, MCNP, GEANT4, etc. Among them, MCNP was developed by the Los Alamos National Laboratory in the United States. It is mainly used for the Monte Carlo simulation program of photon, electron, neutron and other transport problems in materials in three-dimensional geometric structures. Its applicable photon energy range is 1E-3MeV-1E5MeV, electron energy range is 1E-3MeV-1E3MeV, and neutron energy range is 1E-11MeV-20MeV. MCNP has a complete program function, rich reaction cross-section data of various substances, various methods to reduce variance, strong versatility and simple use.

图2示出了本发明基于人工神经网络的点源剂量率修正方法流程图。FIG2 shows a flow chart of a point source dose rate correction method based on an artificial neural network according to the present invention.

如图2所示,首先获取无人机放射性监测系统在不同高度下的能谱,即模拟数据,本申请采取的是利用MCNP软件进行模拟。提取不同高度下能谱数据和剂量率值作为输入参数和输出参数。将部分所述不同高度下能谱数据划分为训练数据,另一部分划分为测试数据。需注意的是测试数据不得少于50组,训练数据需远大于测试数据的数量,如本申请实施例中,测试数据为50组,训练数据为450组。As shown in Figure 2, first obtain the energy spectrum of the unmanned aerial vehicle radioactivity monitoring system at different altitudes, that is, simulated data. The present application adopts the simulation using MCNP software. Extract the energy spectrum data and dose rate values at different altitudes as input parameters and output parameters. The energy spectrum data at different altitudes described in part are divided into training data, and the other part is divided into test data. It should be noted that the test data must not be less than 50 groups, and the training data must be much larger than the number of test data. For example, in the present application embodiment, the test data is 50 groups and the training data is 450 groups.

用训练数据中的输入参数和输出参数构建人工神经网络模型,即BPNN算法模型,将训练数据分别导入人工神经网络模型中进行训练,并获得训练后的人工神经网络模型。An artificial neural network model, namely a BPNN algorithm model, is constructed using the input parameters and output parameters in the training data, and the training data are respectively imported into the artificial neural network model for training, and the trained artificial neural network model is obtained.

再将上述的测试数据分别导入所述训练好的人工神经网络模型中,将得到理想输出结果,将同等高度下理想输出与对应的测试数据进行比较,若误差小于精度期望值,则该人工神经网络训练模型,即为需要的点源空气剂量率修正算法。若误差大于或等于精度期望值,则该人工神经网络训练模型不符合要求,需要重复以上步骤对人工神经网络模型继续训练。此外若误差小于精度期望值,则还需无人机机载辐射监测系统在现场实测到的能谱数据对本修正算法进行验证。Then import the above test data into the trained artificial neural network model respectively, and get the ideal output result. Compare the ideal output at the same height with the corresponding test data. If the error is less than the expected value of accuracy, the artificial neural network training model is the required point source air dose rate correction algorithm. If the error is greater than or equal to the expected value of accuracy, the artificial neural network training model does not meet the requirements, and it is necessary to repeat the above steps to continue training the artificial neural network model. In addition, if the error is less than the expected value of accuracy, the energy spectrum data measured on site by the drone airborne radiation monitoring system is required to verify the correction algorithm.

下面将对每一步骤进行详细的描述。本申请采用MCNP模拟地面点源在不同高度的空气中的放射率情况。需首先第一步选用选择合适的无人机飞行高度区间以及核素的能量区间,选择合适的机载探测系统。Each step will be described in detail below. This application uses MCNP to simulate the radioactivity of ground point sources in the air at different altitudes. The first step is to select the appropriate UAV flight altitude range and the energy range of the nuclide, and select the appropriate airborne detection system.

第二步,使用MCNP模拟获取无人机放射性监测系统在不同高度下的能谱数据。并将该能谱数据划分为测试数据和训练数据。提取训练数据中与空气沉积剂量率相关的数据,对所有与空气沉积剂量率相关的信息数据进行归一化处理,挑选其中对剂量率沉积起主要贡献的主成分个数作为输入参数。放射源在空气中沉积的空气吸收剂量率作为输出参数。In the second step, MCNP simulation is used to obtain the energy spectrum data of the UAV radioactivity monitoring system at different altitudes. The energy spectrum data is divided into test data and training data. The data related to the air deposition dose rate in the training data is extracted, and all the information data related to the air deposition dose rate are normalized, and the number of principal components that contribute to the dose rate deposition is selected as the input parameter. The air absorption dose rate deposited by the radioactive source in the air is used as the output parameter.

第三步,根据上述的输入参数和输出参数,计算隐含层节点数。,构建人工神经网络模型,。将训练数据分别导入上述构建的人工神经网络模型中进行模型训练。得到训练后的人工神经网络模型。The third step is to calculate the number of hidden layer nodes according to the above input parameters and output parameters. Then, an artificial neural network model is constructed. The training data is respectively imported into the artificial neural network model constructed above to perform model training. Then, the trained artificial neural network model is obtained.

构建人工神经网络模型的方法如下:将不同高度下的能谱数据进行归一化处理。提取归一化处理后对剂量率沉积起主要贡献的主成分个数作为输入参数。放射源在空气中沉积的空气吸收剂量率为输出参数;利用输入参数和输出参数计算隐含层节点数。人工神经网络模型的结构为,输入参数-隐含层节点数-输出参数,利用人工神经网络模型的结构构建人工神经网络的模型。The method of constructing an artificial neural network model is as follows: the energy spectrum data at different heights are normalized. The number of principal components that contribute to the dose rate deposition after normalization is extracted as an input parameter. The air absorption dose rate deposited by the radioactive source in the air is the output parameter; the number of hidden layer nodes is calculated using the input parameters and the output parameters. The structure of the artificial neural network model is input parameters-number of hidden layer nodes-output parameters, and the artificial neural network model is constructed using the structure of the artificial neural network model.

其中提取归一化处理后对剂量率沉积起主要贡献的主成分个数的方法为,提取不同高度下的能谱数据中与空气沉积剂量率相关的数据。将相关数据进行标准化处理后获得标准化数据。计算标准化数据之间的相关系数,并组成相关系数矩阵。计算相关系数矩阵的各个特征值;利用各个特征值计算每个所述相关数据的贡献率和累计贡献率,确定当累计贡献率超过阈值时,对应的最少特征值个数,作为起主要贡献的主成分个数。The method for extracting the number of principal components that contribute mainly to the dose rate deposition after normalization is to extract data related to the air deposition dose rate from the energy spectrum data at different heights. Standardized data is obtained after the related data is standardized. The correlation coefficients between the standardized data are calculated and a correlation coefficient matrix is formed. Each eigenvalue of the correlation coefficient matrix is calculated; the contribution rate and cumulative contribution rate of each of the related data are calculated using each eigenvalue, and the minimum number of eigenvalues corresponding to the cumulative contribution rate exceeds the threshold is determined as the number of principal components that contribute mainly.

与空气沉积剂量率相关的数据包括各核素对应的全能峰计数、能量、单逃逸峰、双逃逸峰、湮灭峰和对应的测量时间和高度。The data related to the air deposition dose rate include the total energy peak count, energy, single escape peak, double escape peak, annihilation peak and the corresponding measurement time and height of each nuclide.

将所述相关数据进行标准化处理获得所述标准化数据的方法为,利用如下公式: The method of performing standardization processing on the relevant data to obtain the standardized data is to use the following formula:

其中,其中xij表示的是第i个能谱的第j个指标;表示的是xij经过标准化操作之后的数据;其中是表示的是所有所述能谱数据的平均值;为数据的方差,代表数据差异性,sj是代表的两个数据之间的标准差。n表示的是用来计算的能谱个数。Among them, x ij represents the jth index of the i-th energy spectrum; It represents the data after xij is standardized; It represents the average value of all the energy spectrum data; is the variance of the data, representing the data variability, sj is the standard deviation between the two data represented, and n represents the number of energy spectra used for calculation.

其中平均值和标准差的计算可依照如下公式:The mean and standard deviation can be calculated according to the following formula:

公式中各项的含义同上。The meaning of each term in the formula is the same as above.

计算所述标准化数据之间的相关系数的方法为,利用如下公式:The method for calculating the correlation coefficient between the standardized data is to use the following formula:

其中,rij为所述能谱数据xi与xj的相关系数;xi为第i个能谱数据的平均值;xj为第j个能谱数据的平均值;n表示的是用来进行计算的能谱的个数,xkj表示的是第k个能谱的第j个指标。Among them, rij is the correlation coefficient of the energy spectrum data xi and xj ; xi is the average value of the i-th energy spectrum data; xj is the average value of the j-th energy spectrum data; n represents the number of energy spectra used for calculation, and xkj represents the j-th index of the k-th energy spectrum.

将相关系数组成相关系数矩阵RThe correlation coefficients are combined into a correlation coefficient matrix R

计算相关系数矩阵,并对获得的特征值和特征向量进行排序,即解特征方程:Calculate the correlation coefficient matrix and sort the obtained eigenvalues and eigenvectors, that is, solve the characteristic equation:

|λI-R|=0|λI-R|=0

其中I为对应维度的单位矩阵,计算上式获取特征值λ1,λ2,λ3,...,λ17与其对应的特征向量,并按照大小进行排序。Where I is the unit matrix of the corresponding dimension. The above formula is used to obtain the eigenvalues λ 1 , λ 2 , λ 3 , ..., λ 17 and their corresponding eigenvectors, and they are sorted by size.

接下来计算每一个成分的贡献率以及累计贡献值进行主成分的筛选,第l个主成分的贡献率αl的计算公式为:Next, the contribution rate and cumulative contribution value of each component are calculated to screen the principal components. The calculation formula for the contribution rate α l of the lth principal component is:

k个主成分的累计贡献率Gk的计算公式为:The calculation formula for the cumulative contribution rate G k of k principal components is:

当累计贡献率大于80%时,对应的最少的k值,即为主成分数目,也是输入信息。When the cumulative contribution rate is greater than 80%, the corresponding minimum k value, that is, the number of principal components, is also the input information.

第四步,对训练后的人工神经网络模型进行测试,将测试数据分别导入训练后的人工神经网络模型中,得到理想的输出结果,比较其与对应的测试数据之间的误差,若误差小于设定的精度期望值,则该训练后的网络模型即为所需的修正方法。The fourth step is to test the trained artificial neural network model. The test data are imported into the trained artificial neural network model to obtain the ideal output result, and the error between it and the corresponding test data is compared. If the error is less than the set expected accuracy value, the trained network model is the required correction method.

比较理想的输出结果与对应的测试数据之间的误差的方法依照如下公式:The method for comparing the error between the ideal output result and the corresponding test data is based on the following formula:

其中,mi表示经所述训练好的人工神经网络模型获得的理想输出结果,oi表示经训练好的人工神经网络模型获得的实际输出结果,n表示样本数量。Wherein, mi represents the ideal output result obtained by the trained artificial neural network model, o i represents the actual output result obtained by the trained artificial neural network model, and n represents the number of samples.

按照经验,一般设定模型精度期望值为10-3According to experience, the expected value of model accuracy is generally set at 10 -3 .

此外,可将MCNP软件中提取到的输入层、输出层、隐含层信息一部分用于训练人工神经网络模型,另一部分可在训练网络模型前,用于测试该BN网络模型的识别能力。In addition, part of the input layer, output layer, and hidden layer information extracted from the MCNP software can be used to train the artificial neural network model, and the other part can be used to test the recognition ability of the BN network model before training the network model.

另外如图2所示,利用训练好的人工神经网络,对使用无人机放射性监测系统测量质量的放射源在一定距离处的数据样本进行验证,验证人工神经网络模型的训练质量。In addition, as shown in Figure 2, the trained artificial neural network is used to verify the data samples of the radioactive source at a certain distance whose mass is measured using the UAV radioactivity monitoring system to verify the training quality of the artificial neural network model.

实施例Example

为了便于理解,本申请详细公开了应用本方法进行修正的实施例一。For ease of understanding, the present application discloses in detail Example 1 in which the present method is applied for correction.

步骤一:选择合适的机载探测系统。本发明选用NaI(T1)探测器系统,由该探测器系统为2×2英寸NaI(T1)探测器,多道幅度分析仪及配套MAESTRO软件组成,其中探测器能量范围为30keV至3MeV,能量分辨率为7.7%(Cs-137,662KeV),多道幅度分析仪为1024道。由NaI(T1)探测器探测原理可知,γ粒子与探测器作用主要在晶体部分,图3该NaI(T1)探测器模型图,其晶体几何尺寸及周围包覆材料结构如图,其中铝壳厚度为2.5mm,玻璃材料为SiO2,厚度为2mm,反射层材料为MgO,厚度为0.5mm。在MCNP模拟过程中使用的模型也依据上述数据进行建模。Step 1: Select a suitable airborne detection system. The present invention uses a NaI (T1) detector system, which is composed of a 2×2 inch NaI (T1) detector, a multi-channel amplitude analyzer and a matching MAESTRO software. The detector energy range is 30keV to 3MeV, the energy resolution is 7.7% (Cs-137, 662KeV), and the multi-channel amplitude analyzer is 1024 channels. It can be seen from the detection principle of the NaI (T1) detector that the gamma particles and the detector mainly interact in the crystal part. Figure 3 is a model diagram of the NaI (T1) detector, and its crystal geometric dimensions and surrounding coating material structure are shown in the figure, wherein the aluminum shell thickness is 2.5mm, the glass material is SiO2 , the thickness is 2mm, and the reflective layer material is MgO, the thickness is 0.5mm. The model used in the MCNP simulation process is also modeled based on the above data.

步骤二:根据本发明所提供的方法,模拟选择的无人机飞行高度区间是0米至100米,取值在0-100米之间间隔0.2米取一个高度值,一共选取500组高度数据,相对应可以获取500组的模拟数据;对应的能量信息表示为入射光子的能量,设置为代表不同的核素的能量值,内容包括:Am-241(59.5KeV)、Cs-137(662KeV)、Co-60(1173KeV和1332KeV)、K-40(1460KeV)。Step 2: According to the method provided by the present invention, the simulated selected UAV flight altitude range is 0 meters to 100 meters, and a height value is taken at an interval of 0.2 meters between 0-100 meters. A total of 500 sets of altitude data are selected, and 500 sets of simulation data can be obtained accordingly; the corresponding energy information is expressed as the energy of the incident photon, which is set to represent the energy values of different nuclides, including: Am-241 (59.5KeV), Cs-137 (662KeV), Co-60 (1173KeV and 1332KeV), K-40 (1460KeV).

步骤三:根据本发明所提供的方法,MCNP模拟选用的探测器模型:点源在地面0米处,无人机放射性监测系统的探测器在点源的正上方,距离采取上述的距离区间(如图1中所示)。模拟的源项是各向同性,发射权重为1的点源。用相应的曲面卡和栅元卡建立模型,使用F8卡获取碘化钠晶体探测器的能谱图。根据本发明提供的方法,在MCNP模拟过程中,模拟的能谱的真实性与发射粒子的数量设置相关,可以采取发射粒子数为1×1011进行模拟。Step 3: According to the method provided by the present invention, the detector model selected for MCNP simulation: the point source is at 0 meters on the ground, and the detector of the UAV radioactivity monitoring system is directly above the point source, and the distance takes the above distance interval (as shown in Figure 1). The simulated source term is isotropic, and the emission weight is 1. Use the corresponding surface card and gate element card to establish the model, and use the F8 card to obtain the energy spectrum of the sodium iodide crystal detector. According to the method provided by the present invention, in the MCNP simulation process, the authenticity of the simulated energy spectrum is related to the number of emitted particles, and the number of emitted particles can be set to 1×10 11 for simulation.

步骤四:将模拟获得的500组数据中的450组数据作为训练数据进行BP人工神经网络的训练,将剩下的50组数据作为测试数据测试BP人工神经网络的识别能力。Step 4: Use 450 sets of data from the 500 sets of data obtained by simulation as training data to train the BP artificial neural network, and use the remaining 50 sets of data as test data to test the recognition ability of the BP artificial neural network.

步骤五:根据本发明所提供的方法,首先对MCNP获得的能谱进行归一化处理,再进行能谱特征提取,获得能谱中的Am-241核素对应的全能峰计数和能量;Cs-137核素对应的全能峰计数和能量;Co-60的1173KeV峰对应的全能峰计数和能量;Co-60的1332KeV峰对应的全能峰计数、能量以及单逃逸峰和双逃逸峰;K-40的1460KeV峰对应的全能峰计数、能量以及单逃逸峰和双逃逸峰;再加上能谱中的湮灭峰,对应的测量时间和高度这一共17个与核素的空气沉积剂量率相关的数据。由于能谱信息繁多,提取获得的能谱数据是多维数据,需要对提取到的数据进行降维处理,选出对应剂量率贡献最大的数据。Step 5: According to the method provided by the present invention, the energy spectrum obtained by MCNP is first normalized, and then the energy spectrum features are extracted to obtain the total energy peak count and energy corresponding to the Am-241 nuclide in the energy spectrum; the total energy peak count and energy corresponding to the Cs-137 nuclide; the total energy peak count and energy corresponding to the 1173KeV peak of Co-60; the total energy peak count, energy, single escape peak and double escape peak corresponding to the 1332KeV peak of Co-60; the total energy peak count, energy, single escape peak and double escape peak corresponding to the 1460KeV peak of K-40; plus the annihilation peak in the energy spectrum, the corresponding measurement time and height, a total of 17 data related to the air deposition dose rate of the nuclide. Due to the large amount of energy spectrum information, the extracted energy spectrum data is multidimensional data, and the extracted data needs to be reduced in dimension to select the data with the largest contribution to the corresponding dose rate.

步骤六:根据本发明所提供的方法,选取训练集的450组数据,再根据步骤五所挑选的17个数据作为变量,构建450×17的数据矩阵X。Step 6: According to the method provided by the present invention, 450 sets of data are selected as training sets, and then the 17 data selected in step 5 are used as variables to construct a 450×17 data matrix X.

式中i=1,2,3,...,450,j=1,2,3,...,17。i表示训练集的450个数据,j表示的是提取获得的17个能谱特征值数据。Wherein, i=1, 2, 3, ..., 450, j=1, 2, 3, ..., 17. i represents 450 data of the training set, and j represents 17 energy spectrum eigenvalue data extracted.

步骤七:对步骤六的数据进行处理,获得标准化矩阵。目的是为了消除各个数据之间的量纲化和数量级上的差异,获得标准化的数据其中矩阵内的数据可由下式计算:Step 7: Process the data from step 6 to obtain a standardized matrix. The purpose is to eliminate the differences in dimensionality and magnitude between the data and obtain standardized data. The data in the matrix It can be calculated by the following formula:

其中xij表示的是第i个能谱的第j个指标;表示的是xij经过标准化操作之后的数据;其中是表示的是所有数据的平均值;为数据的方差,代表数据差异性,sj是代表的两个数据之间的标准差。Where x ij represents the jth index of the i-th energy spectrum; It represents the data after xij is standardized; It represents the average value of all data; is the variance of the data, representing the data difference, and sj is the standard deviation between the two data represented.

上式中n表示的是用来进行计算的能谱的个数,n=450。In the above formula, n represents the number of energy spectra used for calculation, n=450.

为了获得每个成分的贡献值,计算矩阵R(相关系数矩阵):In order to obtain the contribution value of each component, the matrix R (correlation coefficient matrix) is calculated:

其中rij的计算公式为:The calculation formula of r ij is:

上式中,rij为原变量xi与xj的相关系数;xi为第i个能谱;xj为第j个能谱特征值;n表示的是用来进行计算的能谱的个数。In the above formula, rij is the correlation coefficient between the original variables xi and xj ; xi is the i-th energy spectrum; xj is the j-th energy spectrum eigenvalue; n represents the number of energy spectra used for calculation.

计算相关系数矩阵获得的特征值和特征向量进行排序,即解特征方程:Calculate the eigenvalues and eigenvectors obtained from the correlation coefficient matrix and sort them, that is, solve the characteristic equation:

|λI-R|=0|λI-R|=0

其中I为对应维度的单位矩阵,计算上式获取特征值λ1,λ2,λ3,...,λ17与其对应的特征向量,并按照大小进行排序。Where I is the unit matrix of the corresponding dimension. The above formula is used to obtain the eigenvalues λ 1 , λ 2 , λ 3 , ..., λ 17 and their corresponding eigenvectors, and they are sorted by size.

接下来计算每一个成分的贡献率以及累计贡献值进行主成分的筛选,第l个主成分的贡献率αl的计算公式为:Next, the contribution rate and cumulative contribution value of each component are calculated to screen the principal components. The calculation formula for the contribution rate α l of the lth principal component is:

k个主成分的累计贡献率Gk的计算公式为:The calculation formula for the cumulative contribution rate G k of k principal components is:

由上式以及本发明提供的实施例计算可得,当主成成分的个数等于12时,累计贡献率大于80%,可以认为,这十二个数据为本次实施例的主成分。It can be calculated from the above formula and the embodiment provided by the present invention that when the number of the main components is equal to 12, the cumulative contribution rate is greater than 80%. It can be considered that these twelve data are the main components of this embodiment.

其中12个主成分分别是模拟中Am-241(59.5KeV)、Cs-137(662KeV)、Co-60(1173KeV和1332KeV)、K-40(1460KeV)四个放射源发射射线的个能量、能谱中5个能量对应的全能峰的计数以及测量时的时间和高度。The 12 main components are the energies of the rays emitted by the four radiation sources Am-241 (59.5KeV), Cs-137 (662KeV), Co-60 (1173KeV and 1332KeV), and K-40 (1460KeV) in the simulation, the counts of the full-energy peaks corresponding to the five energies in the energy spectrum, and the time and height during measurement.

步骤八:将步骤七获得的12个特征参数应用于人工神经网络的输入参数,同时输出参数为各个核素在空气中的沉积空气吸收剂量率,一共是4个输出参数,因此对用的输入节点数N=12,输出节点数M=4。确定隐含层节点数Q可以通过经验公式获得:Step 8: Apply the 12 characteristic parameters obtained in step 7 to the input parameters of the artificial neural network. At the same time, the output parameters are the deposition air absorption dose rate of each nuclide in the air, which is a total of 4 output parameters. Therefore, the number of input nodes used is N = 12, and the number of output nodes is M = 4. The number of hidden layer nodes Q can be determined by the empirical formula:

本方法选择的最优的节点数为10。所以,本发明中的BPNN的结构为12-10-4。The optimal number of nodes selected by this method is 10. Therefore, the structure of the BPNN in the present invention is 12-10-4.

步骤九:将随机分组完毕后的450组训练数据导入人工神经网络模型中,对模型进行训练;Step 9: Import the 450 sets of training data after random grouping into the artificial neural network model to train the model;

步骤十:将测试数据导入已经训练好的人工神经网络模型中,比较模型获取的输出结果与对应的测试数据之间的误差MSE:按照公式1计算,其中n=450;Step 10: Import the test data into the trained artificial neural network model and compare the error MSE between the output results obtained by the model and the corresponding test data: calculated according to formula 1, where n = 450;

步骤十一:将训练获得的测试误差MSE和模型精度期望值进行比较,如果MSE小于10-3,则认为模型训练完成;若MSE大于精度期望值,则需要继续修改网络参数,直到MSE小于精度期望值,训练结束;Step 11: Compare the test error MSE obtained from training with the expected value of model accuracy. If the MSE is less than 10 -3 , the model training is considered complete. If the MSE is greater than the expected value of accuracy, the network parameters need to be modified until the MSE is less than the expected value of accuracy and the training is completed.

步骤十二:根据本发明提供的方法,最后利用已经训练好的人工神经网络模型,对使用无人机放射性监测系统测量Am-241、Cs-137、Co-60这三种人工放射源在1米处的数据样本进行验证,验证网络训练质量。Step 12: According to the method provided by the present invention, the trained artificial neural network model is finally used to verify the data samples of three artificial radiation sources, Am-241, Cs-137, and Co-60, measured at 1 meter using the drone radioactivity monitoring system to verify the quality of network training.

显然,上述实施例仅仅是为清楚地说明本发明所作的举例,而并非对实施方式的限定。对于所属领域的普通技术人员来说,在上述说明的基础上还可以做出其它不同形式的变化或变动。这里无需也无法对所有的实施方式予以穷举。而由此所引申出的显而易见的变化或变动仍处于本发明的保护范围之中。Obviously, the above embodiments are merely examples for clearly explaining the present invention, and are not intended to limit the implementation methods. For those skilled in the art, other different forms of changes or modifications can be made based on the above description. It is not necessary and impossible to list all the implementation methods here. The obvious changes or modifications derived therefrom are still within the protection scope of the present invention.

Claims (7)

1.一种基于人工神经网络的点源剂量率修正方法,其特征在于,包括:1. A point source dose rate correction method based on artificial neural network, characterized by comprising: S1,获取无人机放射性监测系统在不同高度下的能谱数据和对应剂量率值,分别作为输入参数和输出参数;S1, obtaining the energy spectrum data and corresponding dose rate values of the UAV radioactivity monitoring system at different altitudes as input parameters and output parameters respectively; S2,将部分所述不同高度下能谱数据划分为训练数据,另一部分划分为测试数据,利用所述训练数据中的所述输入参数和输出参数构建人工神经网络模型,将所述训练数据导入所述人工神经网络模型中进行训练,获得训练后的人工神经网络模型;S2, dividing part of the energy spectrum data at different heights into training data, and dividing the other part into test data, using the input parameters and output parameters in the training data to build an artificial neural network model, importing the training data into the artificial neural network model for training, and obtaining a trained artificial neural network model; 其中,获取不同高度能谱数据中的输入参数、输出参数,并构建人工神经网络模型的方法为,The method of obtaining input parameters and output parameters in energy spectrum data at different heights and constructing an artificial neural network model is as follows: a,将不同高度下的能谱数据进行归一化处理;a, normalize the energy spectrum data at different heights; b,提取归一化处理后对剂量率沉积起主要贡献的主成分个数作为输入参数,方法为,提取不同高度下的能谱数据中与空气沉积剂量率相关的数据;将所述相关数据进行标准化处理后获得标准化数据;计算所述标准化数据之间的相关系数,并组成相关系数矩阵;计算所述相关系数矩阵的各个特征值;利用所述各个特征值计算每个所述相关数据的贡献率和累计贡献率,当所述累计贡献率大于阈值时,对应的最少的特征值数量,作为起主要贡献的主成分个数,即为输入参数;b. Extracting the number of principal components that contribute mainly to the dose rate deposition after normalization as an input parameter, the method is as follows: extracting data related to the air deposition dose rate from the energy spectrum data at different heights; obtaining standardized data after normalization of the related data; calculating the correlation coefficients between the standardized data and forming a correlation coefficient matrix; calculating each eigenvalue of the correlation coefficient matrix; using the eigenvalues to calculate the contribution rate and cumulative contribution rate of each of the related data, when the cumulative contribution rate is greater than a threshold, the corresponding minimum number of eigenvalues is used as the number of principal components that contribute mainly, which is the input parameter; c,放射源在空气中沉积的空气吸收剂量率为输出参数;c. The air absorbed dose rate deposited by the radiation source in the air is the output parameter; d,利用所述输入参数和所述输出参数计算隐含层;d. calculating a hidden layer using the input parameters and the output parameters; e,所述人工神经网络模型的结构为,输入参数-隐含层-输出参数,利用所述人工神经网络模型的结构构建人工神经网络的模型;e. The structure of the artificial neural network model is input parameter-hidden layer-output parameter, and the artificial neural network model is constructed using the structure of the artificial neural network model; S3,将所述测试数据分别导入所述训练后的人工神经网络模型,得到理想输出结果,比较所述理想输出结果与对应的所述测试数据之间的误差;S3, importing the test data into the trained artificial neural network model respectively to obtain an ideal output result, and comparing the error between the ideal output result and the corresponding test data; 若所述误差大于或等于设定的精度期望值,重复S2、S3,若所述误差小于设定的精度期望值,则训练调试后的人工神经网络模型为点源剂量率修正算法。If the error is greater than or equal to the set expected value of accuracy, S2 and S3 are repeated. If the error is less than the set expected value of accuracy, the trained and debugged artificial neural network model is a point source dose rate correction algorithm. 2.如权利要求1所述的基于人工神经网络的点源剂量率修正方法,其特征在于,获取无人机放射性监测系统在不同高度下的能谱的方法为:使用蒙特卡洛软件模拟无人机放射性监测系统在不同高度下的能谱。2. The point source dose rate correction method based on artificial neural network as described in claim 1 is characterized in that the method for obtaining the energy spectrum of the unmanned aerial vehicle radioactivity monitoring system at different altitudes is: using Monte Carlo software to simulate the energy spectrum of the unmanned aerial vehicle radioactivity monitoring system at different altitudes. 3.如权利要求2所述的基于人工神经网络的点源剂量率修正方法,其特征在于,与空气沉积剂量率相关数据包括:3. The point source dose rate correction method based on artificial neural network according to claim 2, characterized in that the data related to the air deposition dose rate includes: 各核素对应的全能峰计数、能量、单逃逸峰、双逃逸峰、湮灭峰和对应的测量时间和高度。The total energy peak count, energy, single escape peak, double escape peak, annihilation peak and corresponding measurement time and height of each nuclide. 4.如权利要求3所述的基于人工神经网络的点源剂量率修正方法,其特征在于,将所述相关数据进行标准化处理获得所述标准化数据的方法为,利用如下公式:4. The point source dose rate correction method based on artificial neural network according to claim 3 is characterized in that the method of normalizing the relevant data to obtain the standardized data is to use the following formula: 其中,xij表示的是第i个能谱的第j个指标;表示的是xij经过标准化操作之后的数据;其中是表示的是所有所述能谱数据的平均值;sj是代表的两个数据之间的标准差。Among them, x ij represents the jth index of the i-th energy spectrum; It represents the data after xij is standardized; represents the average value of all the energy spectrum data; sj represents the standard deviation between the two data. 5.如权利要求4所述的基于人工神经网络的点源剂量率修正方法,其特征在于,计算所述标准化数据之间的相关系数的方法为,利用如下公式:5. The point source dose rate correction method based on artificial neural network according to claim 4, characterized in that the method for calculating the correlation coefficient between the standardized data is to use the following formula: 其中,rij为所述能谱数据xi与xj的相关系数;为第i个能谱数据的平均值;为第j个能谱数据的平均值;n表示的是用来进行计算的能谱的个数,xkj表示的是第k个能谱的第j个指标,xki表示的是第k个能谱的第i个指标。Wherein, r ij is the correlation coefficient between the energy spectrum data x i and x j ; is the average value of the i-th energy spectrum data; is the average value of the j-th energy spectrum data; n represents the number of energy spectra used for calculation, x kj represents the j-th index of the k-th energy spectrum, and x ki represents the ith index of the k-th energy spectrum. 6.如权利要求5所述的基于人工神经网络的点源剂量率修正方法,其特征在于,所述精度期望值为10-36 . The point source dose rate correction method based on artificial neural network according to claim 5 , wherein the expected accuracy value is 10 −3 . 7.如权利要求1所述的基于人工神经网络的点源剂量率修正方法,其特征在于,还包括:7. The point source dose rate correction method based on artificial neural network according to claim 1, characterized in that it also includes: 若误差小于精度期望值,在实际场景中使用无人机放射性监测系统测量相同放射源在不同高度下的剂量率,并与实际输出结果进行对比,以验证所述人工神经网络模型的质量。If the error is less than the expected value of accuracy, the dose rate of the same radiation source at different altitudes is measured using a UAV radioactivity monitoring system in an actual scenario and compared with the actual output results to verify the quality of the artificial neural network model.
CN202110816233.9A 2021-07-20 2021-07-20 Point source dose rate correction method based on artificial neural network Active CN113536679B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110816233.9A CN113536679B (en) 2021-07-20 2021-07-20 Point source dose rate correction method based on artificial neural network

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110816233.9A CN113536679B (en) 2021-07-20 2021-07-20 Point source dose rate correction method based on artificial neural network

Publications (2)

Publication Number Publication Date
CN113536679A CN113536679A (en) 2021-10-22
CN113536679B true CN113536679B (en) 2024-09-24

Family

ID=78100324

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110816233.9A Active CN113536679B (en) 2021-07-20 2021-07-20 Point source dose rate correction method based on artificial neural network

Country Status (1)

Country Link
CN (1) CN113536679B (en)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113835111B (en) * 2021-11-26 2022-02-11 常州微亿智造科技有限公司 Unmanned aerial vehicle-based radioactive source positioning method and device
CN114966800A (en) * 2022-04-22 2022-08-30 南京航空航天大学 Movable internal irradiation detection method and equipment
CN116010785A (en) * 2022-12-08 2023-04-25 中国辐射防护研究院 ResNet-based gamma energy spectrum data analysis method, storage medium and system

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110133703A (en) * 2019-05-15 2019-08-16 中国辐射防护研究院 A kind of gamma-rays Dosimetry based on BP neural network and Monte carlo algorithm
CN111539324A (en) * 2020-04-23 2020-08-14 重庆建安仪器有限责任公司 Novel nuclide identification method

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR102051576B1 (en) * 2018-03-07 2019-12-05 한국과학기술원 Apparatus and method for identifying multi-radioisotope based on plastic scintillator using Artificial Neural Network

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110133703A (en) * 2019-05-15 2019-08-16 中国辐射防护研究院 A kind of gamma-rays Dosimetry based on BP neural network and Monte carlo algorithm
CN111539324A (en) * 2020-04-23 2020-08-14 重庆建安仪器有限责任公司 Novel nuclide identification method

Also Published As

Publication number Publication date
CN113536679A (en) 2021-10-22

Similar Documents

Publication Publication Date Title
CN113536679B (en) Point source dose rate correction method based on artificial neural network
Miller et al. MCNPX-PoliMi post-processing algorithm for detector response simulations
He et al. Rapid radionuclide identification algorithm based on the discrete cosine transform and BP neural network
CN103135125A (en) Energy spectrum analytical method of X-ray filtering reference radiation
CN110133703A (en) A kind of gamma-rays Dosimetry based on BP neural network and Monte carlo algorithm
CN107657119B (en) Measurement method for improving gamma energy spectrum data quality
CN104599302B (en) The method for obtaining PET crystal energies peak value and setting energy frequency discriminator
CN105182422A (en) Stratum element logging element standard spectrum making and application method
CN108287357A (en) A kind of source peak detection efficient acquisition methods of cylinder bulk detector to point source
CN102713677B (en) Method for identifying a nuclear explosion based on krypton and xenon isotopes
Qin et al. Calibration of an airborne γ-ray spectrometer based on an unmanned aerial vehicle using a point source
CN110687584B (en) LSTM-based rapid nuclide identification method
CN107092028A (en) A kind of computational methods of closed container nuclear material quality
CN116598030A (en) A Neutron Flux Distribution Inversion Method Based on PGNAA Technology
Ghawaly Jr A Datacentric Algorithm for Gamma-ray radiation anomaly detection in unknown Background Environments
Li et al. Quantitative algorithm for airborne gamma spectrum of large sample based on improved shuffled frog leaping–particle swarm optimization convolutional neural network
Wu et al. Design of radionuclides identification algorithm based on sequence Bayesian method
Cai et al. Rapid determination of radon progeny concentration based on artificial neural networks
CN108008437A (en) Method and device for X-ray detection X
Wang et al. Explainable neural network algorithm for rapid radionuclide identification under low count gamma-ray spectrum data
Walt et al. Small-Angle Elastic Scattering of Fast Neutrons and the Electric Polarizability of the Neutron
CN110887853B (en) SGS tomographic efficiency scaling method based on weighted average of voxel attenuation efficiency
Rezaei et al. Numerical evaluation of gamma radiation monitoring
Xiao-Feng et al. Pulse pile-up correction by particle swarm optimization with double-layer parameter identification model in x-ray spectroscopy
CN109669207A (en) Utilize the method for lanthanum bromide detector measurement Neutron Radiation Field neutron energy spectrum

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
CB03 Change of inventor or designer information

Inventor after: Wang Zeyu

Inventor after: Mo Hong

Inventor after: Tang Xiaobin

Inventor after: Gong Pin

Inventor after: Wang Peng

Inventor after: Liang Dajian

Inventor after: Zhou Cheng

Inventor after: Zhu Xiaoxiang

Inventor after: Jiang Ruocheng

Inventor before: Tang Xiaobin

Inventor before: Gong Pin

Inventor before: Mo Hong

Inventor before: Wang Peng

Inventor before: Liang Dajian

Inventor before: Wang Zeyu

Inventor before: Zhou Cheng

Inventor before: Zhu Xiaoxiang

Inventor before: Jiang Ruocheng

CB03 Change of inventor or designer information
GR01 Patent grant
GR01 Patent grant