CN104933245B - 基于神经网络和遗传算法的船用反应堆屏蔽设计优化方法 - Google Patents
基于神经网络和遗传算法的船用反应堆屏蔽设计优化方法 Download PDFInfo
- Publication number
- CN104933245B CN104933245B CN201510329363.4A CN201510329363A CN104933245B CN 104933245 B CN104933245 B CN 104933245B CN 201510329363 A CN201510329363 A CN 201510329363A CN 104933245 B CN104933245 B CN 104933245B
- Authority
- CN
- China
- Prior art keywords
- neural network
- neutral net
- reactor
- shielding
- dose
- 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
Links
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了基于神经网络和遗传算法的船用反应堆屏蔽设计优化方法,利用蒙特卡洛方法模拟计算反应堆中子和光子输运过程,采用MCNP程序对船用反应堆分层屏蔽模型进行模拟计算,得到反应堆屏蔽后的中子和光子通量,并计算出总的当量剂量;根据样本的输入参数个数n1和输出参数个数n2来确定BP神经网络的拓扑结构;根据反应堆屏蔽问题目标确定遗传算法的适应度函数,选择推荐的适应度函数,利用遗传算法强寻优能力,与神经网络相互耦合找到最佳的屏蔽参数。本发明的优点是利用神经网络极强的拟合能力减少了MCNP程序在粒子输运计算过程的时间消耗和利用遗传算法很好的寻优能力,能够在很少的迭代步骤中找到最优屏蔽参数。
Description
技术领域
本发明属于核工程技术领域,涉及一种基于神经网络和遗传算法的船用反应堆屏蔽设计优化方法。
背景技术
船用反应堆屏蔽设计优化主要任务是在保证剂量安全的条件下找到体积小重量轻且屏蔽效果佳的屏蔽方案。在过去的实际工程中,船用反应堆屏蔽设计都是计算许多方案后再从中进行筛选,将花费大量的机器时间和工作准备时间,但是这却往往找不到既满足辐射安全、经济和材料特性等因素最优的屏蔽方案。因此,针对船用反应堆屏蔽优化设计一直是核工程界的热点和挑战。
国际上对于船用反应堆屏蔽设计分析主要有数值计算法、半经验修正的计算模型和智能优化方法三大类。在上个世纪80年代由于计算机方面的限制,反应堆屏蔽设计优化主要集中在数值计算方法上,文献主要有:线性规划、变分法、综合法、梯度法、拉格朗日乘子法和最速下降法等数值方法。近些年来,随着计算机科学的快速发展,为规避数值计算方法的一些问题,国际上相继开展了基于蒙特卡洛程序在反应堆屏蔽计算和利用一些智能算法来对反应堆屏蔽参数进行估算和优化。
对于船用反应堆屏蔽设计优化,数值计算方法主要存在计算模型简单、精度低和计算复杂等缺点;蒙特卡洛程序则存在计算机模拟耗时长和效率低下的缺点。而本发明基于神经网络和遗传算法的船用反应堆屏蔽设计优化方法可以有效解决以上的问题。
发明内容
本发明提供了一种基于神经网络和遗传算法的船用反应堆屏蔽设计优化方法,计算简单、耗时少、精度高。
本发明所采用的技术方案是,基于神经网络和遗传算法的船用反应堆屏蔽设计优化方法,按照以下步骤进行:
步骤1,利用蒙特卡洛方法模拟计算反应堆中子和光子输运过程,采用MCNP程序对船用反应堆分层屏蔽模型进行模拟计算,得到反应堆屏蔽后的中子和光子通量,并计算出总的当量剂量;
步骤2,根据样本的输入参数个数n1和输出参数个数n2来确定BP神经网络的拓扑结构;
步骤3,根据反应堆屏蔽问题目标确定遗传算法的适应度函数,适应度函数有以下推荐类型:
(a)约束型
Fittness=Weight,Dose≤D0
Fittness为适应度值;Weight为反应堆屏蔽层质量;D0为剂量限值。
(b)权重型
Fittness=W1*Weight′+W2*Dose′
Fittness为适应度值;Weight′和Dose′分别为归一化的质量和剂量;W1和W2分别为两者之间的权重,默认均为1;
选择上述推荐的适应度函数,利用遗传算法强寻优能力,与神经网络相互耦合找到最佳的屏蔽参数。
进一步的,步骤2中BP神经网络结构有以下原则:
(a)对于复杂工程问题,神经网络隐含层采用双层神经元
(b)在单层的隐含层神经网络中,整个神经网络神经元个数结构为:
n1→2n1±1→n2
(c)在双层的隐含层神经网络中,整个神经网络神经元个数结构为:
n1→1.5n1→2n1±1→n2
对于BP神经网络训练参数,有以下推荐:
(a)节点传递函数
正切S型传递函数tansig
(b)训练函数
BP算法训练函数trainlm
(c)网络学习函数
梯度下降动量函数learngdm
利用神经网络极强的非线性拟合能力对上述的剂量数据进行训练学习,直到达到误差小于2%的预测精度来替换MCNP程序计算的过程。
本发明的有益效果是:
1.利用遗传算法很好的寻优能力,能够在很少的迭代步骤中找到最优屏蔽参数。
2.利用神经网络极强的拟合能力减少了MCNP程序在粒子输运计算过程的时间消耗。
3.利用整个算法流程搭建了一个GUI可视化界面,降低了设计人员的编程难度,使整个反应堆屏蔽设计工作更加直观、高效。
附图说明
图1是方法流程图;
图2是神经网络预测数据和原始数据的回归曲线图,a、b、c、d分别训练集、验证集、测试集和全部数据的回归图;
图3是神经网络训练过程的误差变化图;
图4是神经网络预测效果对比图;
图5是遗传算法的寻优过程。
具体实施方式
基于神经网络和遗传算法的船用反应堆屏蔽设计优化方法,方法流程如图1所示,其中具体步骤包括:
(1)利用蒙特卡洛(MC)方法模拟计算反应堆中子和光子输运过程,采用反应堆常用的MCNP程序对船用反应堆分层屏蔽模型进行模拟计算,得到反应堆屏蔽后的中子和光子通量,并计算出总的当量剂量。
(2)根据样本的输入参数个数n1和输出参数个数n2来确定BP神经网络的拓扑结构。考虑到实际屏蔽问题的复杂性,BP神经网络结构有以下指导原则:
(a)对于复杂工程问题,神经网络隐含层尽量采用双层神经元
(b)在单层的隐含层神经网络中,整个神经网络神经元个数结构推荐为:
n1→2n1±1→n2
(c)在双层的隐含层神经网络中,整个神经网络神经元个数结构推荐为:
n1→1.5n1→2n1±1→n2
对于BP神经网络训练参数,有以下推荐:
(a)节点传递函数
正切S型传递函数tansig
(b)训练函数
BP算法训练函数trainlm
(c)网络学习函数
梯度下降动量函数learngdm
利用神经网络极强的非线性拟合能力对上述的剂量数据进行训练学习,直到达到误差小于2%的预测精度来替换MCNP程序计算的过程。
(3)根据反应堆屏蔽问题目标确定遗传算法的适应度函数,适应度函数有以下推荐类型:
(a)约束型
Fittness=Weight,Dose≤D0
Fittness为适应度值;Weight为反应堆屏蔽层质量;D0为剂量限值;
(b)权重型
Fittness=W1*Weight′+W2*Dose′
Fittness为适应度值;Weight′和Dose′分别为归一化的质量和剂量;
W1和W2分别为两者之间的权重,默认均为1;
选择上述推荐的适应度函数,利用遗传算法很强的寻优能力,与神经网络相互耦合找到最佳的屏蔽参数(本例中指各个屏蔽层厚度)。
(4)利用Matlab程序设计上述整个流程的GUI界面,编译成一个操作简单、效率高的人机交互软件。
本发明的发明点是:在遗传算法中利用神经网络替换MCNP来计算反应堆屏蔽后的剂量并快速找到反应堆屏蔽最佳参数。
涉及此方法的理论及软件介绍:
1.MCNP程序:MCNP(Monte Carlo N Particle Transport Code)是由美国洛斯阿拉莫斯国家实验室开发的基于蒙特卡罗方法的用于计算三维复杂几何结构中的中子、光子、电子或者耦合中子/光子/电子输运问题的通用软件包。目前此软件已经成为核工程领域反应堆模拟计算不可缺少的软件。
2.BP神经网络:BP(Back Propagation)神经网络是1986年提出的一种按误差逆传播算法训练的多层前馈网络,是目前应用最广泛的神经网络之一。BP网络能学习和存贮大量的输入-输出模式映射关系,而无需事前揭示描述这种映射关系的数学方程。它的学习规则是使用最速下降法,通过反向传播来不断调整网络的权值和阈值,使网络的误差平方和最小。神经网络的非线性映照理论表明神经网络能精度逼近任意非线性连续函数。
3.遗传算法:GA(Genetic algorithm)遗传算法是模仿自然界生物进化机制发展起来的随机全局搜索和优化方法,它借鉴了达尔文的进化论和孟德尔的遗传学说。其本质是一种高效、并行、全局搜索的方法,它能在搜索过程中自动获取和积累有关搜索空间的知识,并自适应的控制搜索过程以求得最 优解。
实施例
本例对四种不同厚度的反应堆屏蔽层进行了最优化计算。假设四种屏蔽材料依次为水、铁、铅和聚乙烯且每种材料的厚度范围是0-15cm。首先通过建模将船用反应堆屏蔽设计优化问题转换为各种不同屏蔽层厚度的优化问题;然后利用MCNP程序计算一系列不同屏蔽层厚度屏蔽后的剂量;进一步使用神经网络进行训练上述样本并达到97%以上的预测精度;最后利用遗传算法进行优化来搜寻最佳屏蔽层厚度方案。
本实例目标方案限制反应堆屏蔽剂量在1e-7条件下,寻找到重量最轻的屏蔽方案。本次迭代找到最佳屏蔽方案厚度是[15 12 1 15]cm,反应堆总重量为3514.2kg,神经网络预测剂量是7.6843e-8,MCNP计算剂量是8.2224E-08。
本发明的效果图
神经网络训练时,把所能利用的数据按照0.7、0.15、0.15的比例分为3个子集。第一个子集是训练样本集,用以计算梯度和修正网络的权值及阈值;第二个子集是验证样本集,在训练的过程中,监控样本集的误差,确认误差在正常地减少。当网络样本集的误差在明显地增加且连续增加的次数达到指定的迭代次数时,训练就会终止,此时,网络返回具有最小样本集误差的权值和阈值;第三个子集是测试样本集,但在网络的训练过程中没有用到测试样本集的误差,而是用于比较各种不同的模型误差。
图2是神经网络预测数据和原始数据的回归曲线图,a,b,c,d四个子图依次为训练集、验证集、测试集和全部的回归图。其中R代表数据的相关系数,越接近1相关性越好,相应的精度越高,从图中看出本次神经网络的训练效果极好。
图3是神经网络训练过程的误差变化图,3个子集的误差随着神经网络迭代次数的增加而逐渐减少,最后趋于平稳的状态。验证集在第128次迭代误差最小达到训练的最佳水平,随后的迭代样本误差在连续增加并达到预设 6次验证则停止神经网络的训练。
图4是神经网络预测效果对比图,神经网络训练完成后预测一组新的数据,然后和mcnp计算数据进行对比,从图中可以看出神经网络预测的准确度比较高,满足实施例的要求。
图5是遗传算法寻优的过程,本实施例适应度函数采取了约束型,在满足剂量条件下寻找到质量最轻的反应堆屏蔽层。随着种群代数的增加,适应度函数逐渐减少并趋于最小值,即寻找到最佳屏蔽厚度。
实施例相关参数
A.MCNP误差:均在2%以内
B.神经网络参数:
(c)网络结构参数
输入层、隐含层、输出层:4:[6 9]:1
(d)节点传递函数
正切S型传递函数tansig
(e)训练函数
BP算法训练函数trainlm
(f)网络学习函数
梯度下降动量函数learngdm
C.遗传算法的参数:
·种群规模:100(20-200)
·字串长度:16(二进制编码)
·交叉概率:0.4(0.3-0.6)
·变异概率:0.16(0.05-0.2)
·中止条件:100代
·适应度函数:约束型
·剂量限值:1e-7。
Claims (2)
1.基于神经网络和遗传算法的船用反应堆屏蔽设计优化方法,其特征在于,按照以下步骤进行:
步骤1,利用蒙特卡洛方法模拟计算反应堆中子和光子输运过程,采用MCNP程序对船用反应堆分层屏蔽模型进行模拟计算,得到反应堆屏蔽后的中子和光子通量,并计算出总的当量剂量;
步骤2,根据样本的输入参数个数n1和输出参数个数n2来确定BP神经网络的拓扑结构;
步骤3,根据反应堆屏蔽问题目标确定遗传算法的适应度函数,适应度函数有以下推荐类型;
(a)约束型
Fittness=Weight,Dose≤D0
Fittness为适应度值;Weight为反应堆屏蔽层质量;D0为剂量限值;Dose为中子和光子穿过反应堆屏蔽层之后的剂量;
(b)权重型
Fittness=W1*Weight′+W2*Dose′
Fittness为适应度值;Weight′和Dose′分别为归一化的质量和剂量;
W1和W2分别为两者之间的权重,默认均为1,
选择上述推荐的适应度函数,利用遗传算法强寻优能力,与神经网络相互耦合找到最佳的屏蔽参数。
2.根据权利要求1所述的基于神经网络和遗传算法的船用反应堆屏蔽设计优化方法,其特征在于,所述步骤2中BP神经网络结构有以下原则:
(a)对于复杂工程问题,神经网络隐含层采用双层神经元;
(b)在单层的隐含层神经网络中,整个神经网络神经元个数结构为:
n1→2n1±1→n2
(c)在双层的隐含层神经网络中,整个神经网络神经元个数结构为:
n1→1.5n1→2n1±1→n2
对于BP神经网络训练参数,有以下推荐:
(d)节点传递函数
正切S型传递函数tansig
(e)训练函数
BP算法训练函数trainlm
(f)网络学习函数
梯度下降动量函数learngdm
利用神经网络极强的非线性拟合能力对上述的剂量数据进行训练学习,直到达到误差小于2%的预测精度来替换MCNP程序计算的过程。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510329363.4A CN104933245B (zh) | 2015-06-15 | 2015-06-15 | 基于神经网络和遗传算法的船用反应堆屏蔽设计优化方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510329363.4A CN104933245B (zh) | 2015-06-15 | 2015-06-15 | 基于神经网络和遗传算法的船用反应堆屏蔽设计优化方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104933245A CN104933245A (zh) | 2015-09-23 |
CN104933245B true CN104933245B (zh) | 2017-11-24 |
Family
ID=54120411
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510329363.4A Active CN104933245B (zh) | 2015-06-15 | 2015-06-15 | 基于神经网络和遗传算法的船用反应堆屏蔽设计优化方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104933245B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US11850445B2 (en) | 2016-09-07 | 2023-12-26 | Elekta, Inc. | System and method for learning models of radiotherapy treatment plans to predict radiotherapy dose distributions |
Families Citing this family (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106442470A (zh) * | 2016-08-31 | 2017-02-22 | 广州博谱能源科技有限公司 | 一种基于激光诱导击穿光谱和遗传神经网络的煤质特性定量分析方法 |
EP3520037B1 (en) | 2016-11-04 | 2024-01-03 | Google LLC | Training neural networks using a variational information bottleneck |
CN106951683B (zh) * | 2017-02-28 | 2020-04-10 | 西安交通大学 | 一种用于核电厂厂房屏蔽计算的高效并行扫描方法 |
CN107230505A (zh) * | 2017-06-21 | 2017-10-03 | 中国核动力研究设计院 | 一种反应堆核功率监测方法及系统 |
FR3085787B1 (fr) * | 2018-09-12 | 2020-10-23 | Framatome Sa | Procede d'exploitation d'un reacteur nucleaire avec calcul du rftc en ligne, reacteur nucleaire correspondant |
CN112925006A (zh) * | 2019-12-05 | 2021-06-08 | 核工业西南物理研究院 | 一种磁约束装置聚变中子实时能谱解谱方法 |
US11077320B1 (en) | 2020-02-07 | 2021-08-03 | Elekta, Inc. | Adversarial prediction of radiotherapy treatment plans |
CN111666718B (zh) * | 2020-06-08 | 2021-10-22 | 南华大学 | 核设施源项活度智能反演方法、装置、设备及存储介质 |
CN111666719B (zh) * | 2020-06-08 | 2022-02-08 | 南华大学 | γ辐射多层屏蔽累积因子计算方法、装置、设备及介质 |
CN112687351B (zh) * | 2021-01-07 | 2023-04-18 | 哈尔滨工业大学 | 一种基于遗传算法-bp神经网络快速预测复合介质微波电磁性能的方法 |
CN113221269B (zh) * | 2021-05-06 | 2022-09-13 | 大连海事大学 | 基于尺寸链的圆柱凸轮接口优化设计方法及其结构 |
CN115577471B (zh) * | 2022-10-14 | 2024-03-05 | 上海交通大学 | 半经验-半定量的小型反应堆轻量化屏蔽方法 |
CN116595885B (zh) * | 2023-05-29 | 2024-03-12 | 上海交通大学 | 小型反应堆轻量化屏蔽深度耦合优化方法 |
CN116595814B (zh) * | 2023-07-19 | 2023-10-31 | 清华大学 | 移动反应堆的屏蔽体设计方法、装置、屏蔽体及介质 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6528316B1 (en) * | 2000-08-02 | 2003-03-04 | Uop Llc | Control of solid catalyst alkylation process using raman spectroscopy |
WO2010054089A1 (en) * | 2008-11-06 | 2010-05-14 | Saudi Arabian Oil Company | Automation and control of energy efficient fluid catalytic cracking processes for maximizing value added products |
CN102426862A (zh) * | 2011-11-17 | 2012-04-25 | 中广核工程有限公司 | Dcs系统中反应堆停堆保护信号可靠性建模方法及系统 |
CN103065056A (zh) * | 2013-01-22 | 2013-04-24 | 中国科学院合肥物质科学研究院 | 一种基于数据场分割的移动人体剂量蒙特卡罗模拟方法 |
CN104298836A (zh) * | 2014-11-06 | 2015-01-21 | 中国科学院合肥物质科学研究院 | 一种基于蒙特卡罗计算的反应堆堆芯迭代设计系统 |
-
2015
- 2015-06-15 CN CN201510329363.4A patent/CN104933245B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6528316B1 (en) * | 2000-08-02 | 2003-03-04 | Uop Llc | Control of solid catalyst alkylation process using raman spectroscopy |
WO2010054089A1 (en) * | 2008-11-06 | 2010-05-14 | Saudi Arabian Oil Company | Automation and control of energy efficient fluid catalytic cracking processes for maximizing value added products |
CN102426862A (zh) * | 2011-11-17 | 2012-04-25 | 中广核工程有限公司 | Dcs系统中反应堆停堆保护信号可靠性建模方法及系统 |
CN103065056A (zh) * | 2013-01-22 | 2013-04-24 | 中国科学院合肥物质科学研究院 | 一种基于数据场分割的移动人体剂量蒙特卡罗模拟方法 |
CN104298836A (zh) * | 2014-11-06 | 2015-01-21 | 中国科学院合肥物质科学研究院 | 一种基于蒙特卡罗计算的反应堆堆芯迭代设计系统 |
Non-Patent Citations (2)
Title |
---|
反应堆主屏蔽的设计与计算;钟文发等;《清华大学学报(自然科学版)》;20010228;第41卷(第2期);第126-128页 * |
屏蔽设计组合优化研究;贾小波等;《核科学与工程》;20101231;第30卷(第4期);第338-343页 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US11850445B2 (en) | 2016-09-07 | 2023-12-26 | Elekta, Inc. | System and method for learning models of radiotherapy treatment plans to predict radiotherapy dose distributions |
Also Published As
Publication number | Publication date |
---|---|
CN104933245A (zh) | 2015-09-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104933245B (zh) | 基于神经网络和遗传算法的船用反应堆屏蔽设计优化方法 | |
Kamali et al. | Automatic calibration of HEC‐HMS using single‐objective and multi‐objective PSO algorithms | |
Long | A constraint handling technique for constrained multi-objective genetic algorithm | |
Martelli et al. | PGS-COM: a hybrid method for constrained non-smooth black-box optimization problems: brief review, novel algorithm and comparative evaluation | |
Lin et al. | Machine learning templates for QCD factorization in the search for physics beyond the standard model | |
Safarzadeh et al. | Loading pattern optimization of PWR reactors using Artificial Bee Colony | |
Poursalehi et al. | Differential harmony search algorithm to optimize PWRs loading pattern | |
CN111666719B (zh) | γ辐射多层屏蔽累积因子计算方法、装置、设备及介质 | |
Chen et al. | Multi-objective optimization strategies for radiation shielding design with genetic algorithm | |
Kattan et al. | Evolving radial basis function networks via gp for estimating fitness values using surrogate models | |
Schlünz et al. | A comparative study on multiobjective metaheuristics for solving constrained in-core fuel management optimisation problems | |
Boddula et al. | Groundwater management using a new coupled model of meshless local Petrov-Galerkin method and modified artificial bee colony algorithm | |
Amigoni et al. | Modeling the management of water resources systems using multi-objective DCOPs | |
Wang et al. | A research on the optimal design of BP neural network based on improved GEP | |
CN111126560A (zh) | 一种基于云遗传算法优化bp神经网络的方法 | |
Liu et al. | A numerical method for interval multi-objective mixed-integer optimal control problems based on quantum heuristic algorithm | |
Yang et al. | Development and application of a master-slave parallel hybrid multi-objective evolutionary algorithm for groundwater remediation design | |
Kim et al. | A preliminary study on applicability of artificial neural network for optimized reflector designs | |
Ziver et al. | Genetic algorithms and artificial neural networks for loading pattern optimisation of advanced gas-cooled reactors | |
CN108133286B (zh) | 一种基于地面沉降替代模型的地下水多目标计算方法 | |
Khoshahval et al. | A new method for multi-objective in core fuel management optimization using biogeography based algorithm | |
Gozalvez et al. | Sensitivity study on determining an efficient set of fuel assembly parameters in training data for designing of neural networks in hybrid genetic algorithms | |
Schirru et al. | Intelligent soft computing in nuclear engineering in Brazil | |
Cross et al. | CRRS: Custom regression and regularisation solver for large-scale linear systems | |
Pevey | Multi-objective Optimization of the Fast Neutron Source by Machine Learning |
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 |