CN106844208A - 用于反应堆物理计算程序适用性验证的方法 - Google Patents
用于反应堆物理计算程序适用性验证的方法 Download PDFInfo
- Publication number
- CN106844208A CN106844208A CN201710041019.4A CN201710041019A CN106844208A CN 106844208 A CN106844208 A CN 106844208A CN 201710041019 A CN201710041019 A CN 201710041019A CN 106844208 A CN106844208 A CN 106844208A
- Authority
- CN
- China
- Prior art keywords
- nuclear reactor
- eff
- multiplication factor
- relative
- effective multiplication
- 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.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F11/00—Error detection; Error correction; Monitoring
- G06F11/36—Preventing errors by testing or debugging software
- G06F11/3604—Software analysis for verifying properties of programs
- G06F11/3608—Software analysis for verifying properties of programs using formal methods, e.g. model checking, abstract interpretation
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Software Systems (AREA)
- Computer Hardware Design (AREA)
- Quality & Reliability (AREA)
- Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Monitoring And Testing Of Nuclear Reactors (AREA)
Abstract
本发明公开了用于反应堆物理计算程序适用性验证的方法,1、对现有的临界实验系统和新型核反应堆系统进行敏感性分析,获得有效增殖系数关于多群截面的相对灵敏度系数向量;2、计算各个临界实验系统和新型核反应堆系统之间的相关性系数,并挑选出满足相似性限值的临界实验系统;3、计算多群截面的相对调整量,使得Monte‑Carlo程序对挑选出来的临界实验系统有效增殖系数的计算结果与实测结果之间的偏差整体达到最小;4、调整Monto‑Carlo程序对新型核反应堆系统有效增殖系数的计算结果,并将其作为“实测结果”的最优估计值,用于验证待验证核反应堆物理计算程序对新型核反应堆系统设计分析的适用性。
Description
技术领域
本发明涉及核反应堆安全技术领域,是一种用于核反应堆物理计算程序适用性验证的方法。
背景技术
核反应堆是一个多尺度、多物理耦合的复杂系统,计算程序广泛地应用于核反应堆的设计分析和研究工作。其中,反应堆物理计算是核反应堆系统设计、安全分析和性能评估的基础,其计算结果将直接影响核反应堆的安全性和经济性。因此,反应堆物理计算程序有着严格的验证和确认方法:首先通过与解析结果对比的方法检验程序中计算模型的求解精度;然后通过与大量临界实验系统实测结果对比的方法检验程序模拟结果的正确性和真实性。
随着核能的不断发展和应用需求的不断扩大,新型核反应堆设计方案被不断地提出,设计方案的可行性和安全性需要利用成熟的反应堆物理计算程序对其进行研究和分析。但是,现有的成熟的反应堆物理计算程序的验证和确认工作是基于传统的核反应堆的临界实验系统建立的,是否适用于新型核反应堆系统的设计分析需要增加对程序的验证和确认工作。因此,将成熟的反应堆物理计算程序用于新型核反应堆系统的设计分析,需要进一步对其进行程序的适用性验证。
根据传统的反应堆物理计算程序的验证和确认方法,需要建造新型核反应堆对应的临界实验系统,并测量得到这些临界实验系统的实测结果用于程序的验证和确认工作。新型核反应堆临界实验系统的建造不仅工程代价高、建设周期长,而且存在一定的风险性,远无法满足核能的快速、安全的发展需求。因此,亟需研究一种反应堆物理计算程序适用性验证的方法,科学快速地评估成熟的核反应堆物理计算程序对新型核反应堆系统设计分析的适用性。
发明内容
为了在不建造新型核反应堆临界实验系统的前提下实现对成熟的反应堆物理计算程序适用性的验证研究,本发明的目的在于提供一种用于核反应堆物理计算程序适用性验证的方法,基于敏感性和不确定性分析技术和已有的核反应堆系统临界实验的实测结果,能够科学快速地验证反应堆物理计算程序对新型核反应堆系统的适用性。
为了实现对反应堆物理计算程序的适用性验证,本发明的技术方案概括如下:
一种用于反应堆物理计算程序适用性验证的方法,包括如下步骤:
步骤1:对现有的临界实验系统和新型核反应堆系统进行敏感性分析,获得临界实验系统和新型核反应堆系统的有效增殖系数keff关于多群截面的相对灵敏度系数向量;将有效增殖系数keff关于多群截面的相对灵敏度系数写成向量的形式,表示为
式中:
σ——所有核素、所有类型反应道和所有能群的多群截面组成的向量;
σi——所有核素、所有类型反应道和所有能群的多群截面向量σ中第i个截面;
——keff关于多群截面σi的相对灵敏度系数;
M——所有核素、所有类型反应道和所有能群数目的加和值,表示为:
式中:
Ng——多群截面的能群数目;
Niso——不同核素的数目;
——核素i所有反应道类型数目;
用和分别表示临界实验系统和新型核反应堆系统的有效增殖系数keff的相对灵敏度系数向量;
步骤2:基于步骤1中的相对灵敏度系数向量和进行不确定性分析计算,获得新型核反应堆系统和各个临界实验系统之间的相关性系数;所述的不确定性分析旨在计算新型核反应堆系统和各个临界实验系统的有效增殖系数keff之间的相对协方差矩阵,由公式(4)计算:
式中:
Cσσ——所有核素、所有类型反应道和所有能群的多群截面相对协方差矩阵;
Ca,e——临界实验系统和新型核反应堆系统的有效增殖系数keff之间的相对协方差矩阵;
var(a,a)——新型核反应堆系统的有效增殖系数keff的相对方差;
var(e,e)——临界实验系统的有效增殖系数keff的相对方差;
cov(a,e)——临界实验系统和新型核反应堆系统的有效增殖系数keff的相对协方差,数值上与cov(e,a)相同;
根据相对协方差矩阵Ca,e,新型核反应堆系统和临界实验系统的相关性系数ca,e按照公式(5)计算得到:
相关性系数ca,e表示了临界实验系统和新型核反应堆系统在中子学层面上的相似程度,该值越接近1.0表示系统之间的相似程度越高,该值等于1.0表示两个系统完全相同;根据相似性挑选限值ξρ,挑选出满足条件ca,e≥ξρ的临界实验作为反应堆物理计算程序适用性验证的实验基础;
步骤3:采用Monte‐Carlo程序对步骤2中挑选出来的临界实验系统进行精确建模,获得临界实验系统的有效增殖系数keff的计算结果k=[k1,k2,…,kI],其中I表示挑选后的临界实验系统的数目;对应的临界实验系统的有效增殖系数keff的实测结果表示为m=[m1,m2,…,mI];采用广义线性最小二乘方法,对多群截面进行调整,使得Monte‐Carlo程序对临界实验系统的有效增殖系数keff的计算结果k’与实测结果m之间的整体偏差达到最小;由此计算得到多群截面的相对调整量,表示为:
式中:
δσ——多群截面的相对调整量;
Cσσ——所有核素、所有类型反应道和所有能群的多群截面相对协方差矩阵;
Sk,σ——所有临界实验系统的有效增殖系数keff关于多群截面的相对灵敏度系数向量组成的矩阵:
d——所有临界实验系统的有效增殖系数keff计算结果和实测结果之间的相对偏差,表示为:
Cdd——所有临界实验系统的有效增殖系数keff偏差的相对协方差矩阵;
——所有临界实验系统的有效增殖系数keff关于多群截面的相对灵敏度系数向量组成的矩阵Sk,σ的转置矩阵;
步骤4:使用步骤3中得到的多群截面的相对调整量δσ,对Monte‐Carlo程序计算的新型核反应堆系统的计算结果keff,a进行调整,作为新型核反应堆系统的有效增殖系数keff“实测结果”的最优估计值ma;所述的新型核反应堆系统的有效增殖系数keff“实测结果”的最优估计值表示为:
式中:
ma——新型核反应堆系统的有效增殖系数keff的“实测结果”的最优估计值;
keff,a——新型核反应堆系统的有效增殖系数keff的计算结果;
δσi——多群截面σi的相对调整量,向量δσ中第i个值;
使用待验证的核反应堆物理计算程序对新型核反应堆系统进行精确建模和模拟,计算结果表示为ka;采用新型核反应堆系统的有效增殖系数keff“实测结果”的最优估计值ma对ka进行验证,定义计算偏差为|ka-ma|,检验计算偏差是否在计算接受误差范围,即检验偏差限值ε之内:若|ka-ma|>ε,则待验证的反应堆物理计算程序不适应于新型核反应堆的设计分析;若|ka-ma|≤ε,则待验证的反应堆物理计算程序适用于新型核反应堆系统的设计分析。
步骤2所述相似性挑选限值ξρ取值为ξρ=0.9。
步骤4所述检验偏差限值ε取值为0.5%。
与传统的反应堆物理计算程序验证和确认方法相比,本发明有如下突出优点:
1、本发明无需建造新型核反应堆相关的临界实验系统,节约工程建设成本,并大大减少程序适用性验证周期。
2、基于敏感性和不确定性分析技术,从中子学角度量化临界实验系统和新型核反应堆系统之间的相关性,从而保障用于反应堆物理计算程序适用性验证的临界实验的可靠性。
3、采用核数据库调整技术,获得新型核反应堆系统有效增殖系数keff“实测结果”的最优估计值,作为反应堆物理计算程序适用性验证的基础。
附图说明
图1为反应堆物理计算程序适用性验证流程图。
具体实施方式
本发明基于敏感性和不确定性分析技术,采用广义线性最小二乘方法,在不建造新型核反应堆临界实验系统的前提下,只使用现有的核反应堆临界实验系统的实测结果,实现对核数据库的调整并用于评估新型核反应堆系统keff“实测结果”的最优估计值,用于对核反应堆物理计算程序适用性的验证工作,本发明包括以下方面:
1)采用敏感性分析技术,计算分析临界实验系统和新型核反应堆系统的有效增殖系数keff关于多群截面的相对灵敏度系数向量;
2)基于敏感性和不确定性分析技术,量化临界实验系统和新型核反应堆系统在中子学层面上的相似性,挑选用于程序适用性验证的临界实验系统,相比传统的专家评估的方法更加严谨科学;
3)采用广义线性最小二乘方法,对核数据库进行调整,实现Monte‐Carlo程序对所有临界实验系统的有效增殖系数keff的计算结果与实测结果整体偏差达到最小;
4)基于调整后的多群截面和新型核反应堆系统的有效增殖系数keff关于多群截面的相对灵敏度系数向量,计算评估新型核反应堆系统的有效增殖系数keff“实测结果”的最优估计值,并用于检验待验证的程序对新型核反应堆系统的有效增殖系数keff的计算精度,验证程序对新型核反应堆系统的适用性。
如图1所示,本发明用于反应堆物理计算程序适用性验证的方法,包括如下步骤:
步骤1:使用核反应堆物理计算敏感性和不确定性分析程序对现有的临界实验系统和新型核反应堆系统进行敏感性分析,计算获得临界实验系统和新型核反应堆系统的有效增殖系数keff关于多群截面的相对灵敏度系数向量;
步骤1中keff关于多群截面的相对灵敏度系数如公式(1)所示:
式中:
keff——有效增殖系数;
i——核素标识;
x——反应道类型;
g——能群第g群;
——核素i的类型为x的反应道第g群的截面大小;
——keff关于核素i的类型为x的反应道第g群截面的相对灵敏度系数;
将上述有效增殖系数keff关于所有核素、所有反应道类型和所有能群的相对灵敏度系数写成向量的形式可以表示如公式(2)所示:
式中:
σ——所有核素、所有类型反应道和所有能群的多群截面组成的向量;
σi——所有核素、所有类型反应道和所有能群的多群截面向量中第i个截面;
——keff关于截面σi的相对灵敏度系数;
M——所有核素、所有类型反应道和所有能群数目加和值;
M可表示如公式(3)所示:
式中:
Ng——能群数目;
Niso——核素数目;
——核素i的反应道类型数目;
分别用和表示临界实验系统和新型核反应堆系统的有效增殖系数keff关于多群截面的相对灵敏度系数向量;
步骤2:基于步骤1中临界实验系统和新型核反应堆系统的有效增殖系数keff关于多群截面的相对灵敏度系数向量和结合多群截面的相对协方差矩阵,对新型核反应堆系统和各个临界实验系统进行不确定性分析,计算得到对应的临界实验和新型核反应堆系统的有效增殖系数keff的相对协方差矩阵;所述的临界实验系统和新型核反应堆系统的有效增殖系数keff的相对协方差矩阵如公式(4)所示:
式中:
Cσσ——所有核素、所有类型反应道和所有能群的多群截面相对协方差矩阵;
Ca,e——临界实验系统和新型核反应堆系统的有效增殖系数keff的相对协方差矩阵;
var(a,a)——新型核反应堆系统的有效增殖系数keff的相对方差;
var(e,e)——临界实验系统的有效增殖系数keff的相对方差;
cov(a,e)——临界实验系统和新型核反应堆系统的有效增殖系数keff的相对协方差,数值上与cov(e,a)相同;
根据相对协方差矩阵Ca,e,新型核反应堆系统和临界实验系统的相关性系数ca,e按照公式(5)计算得到:
相关性系数ca,e表示了临界实验系统和新型核反应堆系统在中子学层面上的相似程度,该值越接近1.0表示系统之间的相似程度越高,该值等于1.0表示两个系统完全相同;根据相似性挑选限值ξρ,挑选出满足条件ca,e≥ξρ的临界实验作为反应堆物理计算程序适用性验证的实验基础;一般取ξρ=0.9作为相似性挑选的限值。
步骤3:使用Monte‐Carlo程序对步骤2中经过挑选的临界实验系统进行精确建模,得到临界实验系统的有效增殖系数keff的计算值k=[k1,k2,…,kI],其中I表示挑选后的临界实验系统的数目;对应的临界实验系统的有效增殖系数keff的实测结果表示为m=[m1,m2,…,mI];采用广义线性最小二乘方法对多群截面进行调整,使得采用调整后的临界实验系统的有效增殖系数keff的计算结果k′=[k′1,k′2,...,k′I]与keff的实测结果之间的偏差整体达到最小;
步骤3中,Monte‐Carlo程序被视为核反应堆物理计算程序的对标程序,其理论模型和建模过程误差可忽略不计,核数据库的误差是其计算结果误差的主要来源,因此Monte‐Carlo程序计算得到的keff可表示为关于核数据库的函数ki(σ)(i=1,2,…,I);当核数据库σ存在一定的相对调整量δσ,在一阶线性条件下,调整后的临界实验系统计算结果k′如公式(6)所示:
上式:
ki——第i个临界实验系统初始keff的计算结果;
δσm——多群截面σm的相对调整量;
——ki关于截面σm的相对灵敏度系数;
临界实验系统keff的初始计算结果k与实测结果之间的偏差如公式(7)所示:
调整后的临界实验系统keff的计算结果k′与实测结果之间的偏差如公式(8)所示:
按照公式(6)~(8),获得如公式(9)所示的关系成立:
y=d+Sk,σδσ 公式(9)
广义线性最小二乘方法用于核数据库调整的目标在于:基于广义线性最小二乘算法,调整核数据库使得公式(9)中y整体达到最小;此时,广义线性最小二乘算法定义二次损失函数如公式(10)所示:
式中:
Q(δσ,y)——二次损失函数;
Cmm——实测结果的相对协方差矩阵;
临界实验实测结果的相对协方差矩阵如公式(11)所示:
式中:
cov(mi,mj)——临界实验系统实测结果mi和mj之间的协方差;
根据广义线性最小二乘算法的核数据库调整目标,可以求解得到核数据库的相对调整量如公式(12)所示:
式中:
Cdd——偏差d的相对协方差矩阵;
临界实验系统的有效增殖系数keff的初始计算结果k与实测结果之间的偏差的相对协方差矩阵可表示如公式(13)所示:
Cdd=Cmm+Ckk 公式(13)
式中:
Ckk——临界实验系统的有效增殖系数keff计算结果相对协方差矩阵,其计算如公式(14)所示:
式中:
Sk,σ——所有临界实验系统的有效增殖系数keff相对灵敏度系数向量组成的矩阵:
根据上述公式计算得到核数据截面的相对调整量如公式(12)所示。
步骤4:基于步骤3中得到的多群截面相对调整量δσ,使用新型核反应堆系统的有效增殖系数keff关于多群截面的相对灵敏度系数向量采用一阶线性关系,可计算得到调整后的新型核反应堆系统的计算结果k′eff,作为新型核反应堆系统的有效增殖系数keff“实测结果”的最佳估计值ma;采用新型核反应堆系统的有效增殖系数keff“实测结果”的最佳估计值ma对待验证的程序进行验证,检验其对新型核反应堆系统的适用性;
步骤4中,新型核反应堆系统的有效增殖系数keff“实测结果”的最优估计值如公式(16)所示:
式中:
keff,a——Monte‐Carlo程序对新型核反应堆系统初始keff计算结果;
ma——新型核反应堆系统调整后的keff,可作为实测结果的最佳估计值;
使用ma对待验证的程序进行检验,设待验证程序对新型核反应堆系统keff的计算结果表示为ka,检验偏差限值为ε:若|ka-ma|>ε,则待验证的反应堆物理计算程序不适应于新型核反应堆的设计分析;若|ka-ma|≤ε,则待验证的反应堆物理计算程序适用于新型核反应堆系统的设计分析。
基于以上的理论模型,本发明采用标准FORTRAN90语言编制了反应堆物理计算程序适用性验证程序,该程序的计算流程如图1所示。
在该程序的计算中,敏感性分析用于计算得到初始(未经挑选的)临界实验系统和新型核反应堆系统keff计算结果关于多群截面的相对灵敏度系数向量;基于keff计算结果关于多群截面的相对灵敏度系数向量和多群截面的相对协方差数据,不确定性分析用于计算各个临界实验系统与新型核反应堆系统之间的相关性系数,对临界实验系统进行筛选,挑选出满足条件的临界实验系统;基于挑选出的临界实验系统,核数据库调整采用广义线性最小二乘算法,获得多群截面的相对调整量,实现所有临界实验系统keff计算结果和实测结果之间的偏差整体达到最小;基于新型核反应堆系统keff关于多群截面的相对灵敏度系数向量和多群截面的相对调整量,计算得到新型核反应堆系统keff“实测结果”的最优估计值,用于对待验证程序进行验证,最终实现核反应堆物理计算程序适用性验证研究。
Claims (3)
1.一种用于反应堆物理计算程序适用性验证的方法,其特征在于:包括如下步骤:
步骤1:对现有的临界实验系统和新型核反应堆系统进行敏感性分析,获得临界实验系统和新型核反应堆系统的有效增殖系数keff关于多群截面的相对灵敏度系数向量;将有效增殖系数keff关于多群截面的相对灵敏度系数写成向量的形式,表示为
式中:
σ——所有核素、所有类型反应道和所有能群的多群截面组成的向量;
σi——所有核素、所有类型反应道和所有能群的多群截面向量σ中第i个截面;
——keff关于多群截面σi的相对灵敏度系数;
M——所有核素、所有类型反应道和所有能群数目的加和值,表示为:
式中:
Ng——多群截面的能群数目;
Niso——不同核素的数目;
——核素i所有反应道类型数目;
用和分别表示临界实验系统和新型核反应堆系统的有效增殖系数keff的相对灵敏度系数向量;
步骤2:基于步骤1中的相对灵敏度系数向量和进行不确定性分析计算,获得新型核反应堆系统和各个临界实验系统之间的相关性系数;所述的不确定性分析旨在计算新型核反应堆系统和各个临界实验系统的有效增殖系数keff之间的相对协方差矩阵,由公式(4)计算:
式中:
Cσσ——所有核素、所有类型反应道和所有能群的多群截面相对协方差矩阵;
Ca,e——临界实验系统和新型核反应堆系统的有效增殖系数keff之间的相对协方差矩阵;
var(a,a)——新型核反应堆系统的有效增殖系数keff的相对方差;
var(e,e)——临界实验系统的有效增殖系数keff的相对方差;
cov(a,e)——临界实验系统和新型核反应堆系统的有效增殖系数keff的相对协方差,数值上与cov(e,a)相同;
根据相对协方差矩阵Ca,e,新型核反应堆系统和临界实验系统的相关性系数ca,e按照公式(5)计算得到:
相关性系数ca,e表示了临界实验系统和新型核反应堆系统在中子学层面上的相似程度,该值越接近1.0表示系统之间的相似程度越高,该值等于1.0表示两个系统完全相同;根据相似性挑选限值ξρ,挑选出满足条件ca,e≥ξρ的临界实验系统作为反应堆物理计算程序适用性验证的基础;
步骤3:采用Monte-Carlo程序对步骤2中挑选出来的临界实验系统进行精确建模,获得临界实验系统的有效增殖系数keff的计算结果k=[k1,k2,…,kI],其中I表示挑选后的临界实验系统的数目;对应的临界实验系统的有效增殖系数keff的实测结果表示为m=[m1,m2,…,mI];采用广义线性最小二乘方法,对多群截面进行调整,使得Monte-Carlo程序对临界实验系统的有效增殖系数keff的计算结果k’与实测结果m之间的整体偏差达到最小;由此计算得到多群截面的相对调整量,表示为:
式中:
δσ——多群截面的相对调整量;
Cσσ——所有核素、所有类型反应道和所有能群的多群截面相对协方差矩阵;
Sk,σ——所有临界实验系统的有效增殖系数keff关于多群截面的相对灵敏度系数向量组成的矩阵:
d——所有临界实验系统的有效增殖系数keff计算结果和实测结果之间的相对偏差,表示为:
Cdd——所有临界实验系统的有效增殖系数keff偏差的相对协方差矩阵;
——所有临界实验系统的有效增殖系数keff关于多群截面的相对灵敏度系数向量组成的矩阵Sk,σ的转置矩阵;
步骤4:使用步骤3中得到的多群截面的相对调整量δσ,对Monte-Carlo程序计算的新型核反应堆系统的计算结果keff,a进行调整,作为新型核反应堆系统的有效增殖系数keff“实测结果”的最优估计值ma;所述的新型核反应堆系统的有效增殖系数keff“实测结果”的最优估计值表示为:
式中:
ma——新型核反应堆系统的有效增殖系数keff的“实测结果”的最优估计值;
keff,a——新型核反应堆系统的有效增殖系数keff的计算结果;
δσi——多群截面σi的相对调整量,向量δσ中第i个值;
使用待验证的核反应堆物理计算程序对新型核反应堆系统进行精确建模和模拟,计算结果表示为ka;采用新型核反应堆系统的有效增殖系数keff“实测结果”的最优估计值ma对ka进行验证,定义计算偏差为|ka-ma|,检验计算偏差是否在计算接受误差范围即检验偏差限值ε之内:若|ka-ma|>ε,则待验证的反应堆物理计算程序不适应于新型核反应堆的设计分析;若|ka-ma|≤ε,则待验证的反应堆物理计算程序适用于新型核反应堆系统的设计分析。
2.根据权利要求1所述的一种用于反应堆物理计算程序适用性验证的方法,其特征在于:步骤2所述相似性挑选限值ξρ取值为ξρ=0.9。
3.根据权利要求1所述的一种用于反应堆物理计算程序适用性验证的方法,其特征在于:步骤4所述检验偏差限值ε取值为0.5%。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710041019.4A CN106844208B (zh) | 2017-01-17 | 2017-01-17 | 用于反应堆物理计算程序适用性验证的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710041019.4A CN106844208B (zh) | 2017-01-17 | 2017-01-17 | 用于反应堆物理计算程序适用性验证的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106844208A true CN106844208A (zh) | 2017-06-13 |
CN106844208B CN106844208B (zh) | 2020-01-14 |
Family
ID=59119218
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710041019.4A Active CN106844208B (zh) | 2017-01-17 | 2017-01-17 | 用于反应堆物理计算程序适用性验证的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106844208B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109063233A (zh) * | 2018-06-15 | 2018-12-21 | 中国核电工程有限公司 | 一种用蒙特卡洛法评定核素含量对keff不确定度影响的方法 |
CN109063231A (zh) * | 2018-06-15 | 2018-12-21 | 中国核电工程有限公司 | 基于GUM导则的核素不确定度对临界系统keff影响的评定方法 |
CN110377639A (zh) * | 2019-07-18 | 2019-10-25 | 西安交通大学 | 一种最小化核反应堆物理不确定性分析抽样样本量的方法 |
CN114065126A (zh) * | 2021-10-26 | 2022-02-18 | 中国核动力研究设计院 | 降低核数据相关计算不确定度方法、装置、设备及介质 |
CN114491905A (zh) * | 2020-10-27 | 2022-05-13 | 中国核动力研究设计院 | 一种基于蒙卡抽样的相似性评价方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20140133620A1 (en) * | 2012-11-13 | 2014-05-15 | Westinghouse Electric Company Llc | Method of validating nuclear reactor in-vessel detector output signals |
CN104850750A (zh) * | 2015-05-27 | 2015-08-19 | 东北大学 | 一种核电站反应堆保护系统可靠性分析方法 |
CN106128518A (zh) * | 2016-06-24 | 2016-11-16 | 西安交通大学 | 一种获取高精度的快中子反应堆组件少群截面的方法 |
CN106126480A (zh) * | 2016-06-24 | 2016-11-16 | 西安交通大学 | 一种获取反应堆多群数据库中的多群Pn散射矩阵的方法 |
-
2017
- 2017-01-17 CN CN201710041019.4A patent/CN106844208B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20140133620A1 (en) * | 2012-11-13 | 2014-05-15 | Westinghouse Electric Company Llc | Method of validating nuclear reactor in-vessel detector output signals |
CN104850750A (zh) * | 2015-05-27 | 2015-08-19 | 东北大学 | 一种核电站反应堆保护系统可靠性分析方法 |
CN106128518A (zh) * | 2016-06-24 | 2016-11-16 | 西安交通大学 | 一种获取高精度的快中子反应堆组件少群截面的方法 |
CN106126480A (zh) * | 2016-06-24 | 2016-11-16 | 西安交通大学 | 一种获取反应堆多群数据库中的多群Pn散射矩阵的方法 |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109063233A (zh) * | 2018-06-15 | 2018-12-21 | 中国核电工程有限公司 | 一种用蒙特卡洛法评定核素含量对keff不确定度影响的方法 |
CN109063231A (zh) * | 2018-06-15 | 2018-12-21 | 中国核电工程有限公司 | 基于GUM导则的核素不确定度对临界系统keff影响的评定方法 |
CN109063233B (zh) * | 2018-06-15 | 2023-05-16 | 中国核电工程有限公司 | 一种用蒙特卡洛法评定核素含量对keff不确定度影响的方法 |
CN109063231B (zh) * | 2018-06-15 | 2023-06-27 | 中国核电工程有限公司 | 基于GUM导则的核素不确定度对临界系统keff影响的评定方法 |
CN110377639A (zh) * | 2019-07-18 | 2019-10-25 | 西安交通大学 | 一种最小化核反应堆物理不确定性分析抽样样本量的方法 |
CN110377639B (zh) * | 2019-07-18 | 2021-06-22 | 西安交通大学 | 一种最小化核反应堆物理不确定性分析抽样样本量的方法 |
CN114491905A (zh) * | 2020-10-27 | 2022-05-13 | 中国核动力研究设计院 | 一种基于蒙卡抽样的相似性评价方法 |
CN114491905B (zh) * | 2020-10-27 | 2023-11-28 | 中国核动力研究设计院 | 一种基于蒙卡抽样的相似性评价方法 |
CN114065126A (zh) * | 2021-10-26 | 2022-02-18 | 中国核动力研究设计院 | 降低核数据相关计算不确定度方法、装置、设备及介质 |
CN114065126B (zh) * | 2021-10-26 | 2023-09-26 | 中国核动力研究设计院 | 降低核数据相关计算不确定度方法、装置、设备及介质 |
Also Published As
Publication number | Publication date |
---|---|
CN106844208B (zh) | 2020-01-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106844208A (zh) | 用于反应堆物理计算程序适用性验证的方法 | |
CN110472846A (zh) | 核电厂热工水力安全分析最佳估算加不确定性方法 | |
CN105224737A (zh) | 一种空间目标轨道改进初值修正方法 | |
CN108090614B (zh) | 一种基于相关系数的空间风场预测模型建立方法 | |
CN106780104A (zh) | 一种基于概率统计的平均风向计算方法 | |
MXPA06015255A (es) | Metodo para determinar imagenes para limites de operación para una operacion de reactor nuclear. | |
CN104102836A (zh) | 一种电力系统快速抗差状态估计方法 | |
CN104820204A (zh) | 一种减小偏差的加权最小二乘定位方法 | |
CN106886649A (zh) | 一种多元热流体吞吐注入参数优化方法 | |
CN104406715A (zh) | 一种遥感估算地表感热/潜热通量的精度评价方法及系统 | |
CN104766250A (zh) | 一种管廊管道的风险因素权重值计算方法 | |
CN103106332B (zh) | 一种测量不确定度的分析方法 | |
Ke et al. | Non‐Gaussian characteristics and extreme distribution of fluctuating wind pressures on large cylindrical–conical steel cooling towers | |
Madadi et al. | Probabilistic real-time dynamic line rating forecasting based on dynamic stochastic general equilibrium with stochastic volatility | |
Huang et al. | Peak factors of non‐Gaussian wind forces on a complex‐shaped tall building | |
CN104021288B (zh) | 用于导管架平台频谱疲劳分析的基本波确定方法 | |
CN104008258A (zh) | 基于位移形态的钢结构火灾温度场逆向推定方法 | |
Shonder et al. | Bayesian analysis of savings from retrofit projects | |
CN111310121A (zh) | 一种新能源出力概率预测方法和系统 | |
CN101929684A (zh) | 混烧低热值煤气机组煤耗计算方法 | |
CN115456325A (zh) | 一种用于非煤矿山的灾害设防能力的分析方法 | |
CN104021311A (zh) | 一种基于Hermite函数约束的数据融合计算方法 | |
CN106407160A (zh) | 含多个风电场的电力系统概率潮流联合分布的计算方法 | |
CN108268694B (zh) | 核事故辐射数据的模拟方法、装置、终端设备和存储介质 | |
Fynan | Uncertainty Quantification for Reactor Safety Analysis. |
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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |