【技术实现步骤摘要】
本专利技术涉及核反应堆物理数值计算，具体涉及一种融合实测数据的压水堆燃料组件弯曲模型反演方法。
技术介绍
1、数值模拟是实现压水堆堆芯燃料管理和运行安全技术支持的重要手段，通常需要根据堆芯设计构建理论燃料组件模型。然而在反应堆运行时，燃料组件由于轴向辐照生长、横向水力学作用等会不可避免地偏离正常的竖直状态而发生弯曲，导致数值模拟建立的理论燃料组件模型和实际堆芯存在一定差异。当燃料组件弯曲量过大时，将会影响反应堆的安全运行，出现如堆内装卸料困难、局部功率过高、控制棒下插困难、象限功率倾斜等问题。因此，针对燃料组件弯曲现象进行高精度数值模拟是十分必要的。但是，在堆芯运行过程中，无法通过测量手段获得真实的燃料组件弯曲状态，从而无法获得真实的燃料组件弯曲模型。值得注意的是，在堆芯运行中会通过测量手段进行关键物理量的测量，如功率分布、临界硼浓度等。因此，如何利用现有实测数据获得真实的燃料组件弯曲模型对于提高压水堆的数值模拟精度及安全运行至关重要。
技术实现思路
1、针对上述真实燃料组件弯曲模型无法获得的问题，本专利技术的目的在于提供一种融合实测数据的压水堆燃料组件弯曲模型反演方法，通过随机扰动获得不同的燃料组件弯曲模型，基于压水堆堆芯计算程序获得对应的实测数据物理量计算值；将随机扰动获得的燃料组件弯曲模型作为输入，对应的实测数据物理量计算值作为输出，构建数据集；基于全连接神经网络训练得到燃料组件弯曲模型与实测数据物理量计算值之间的显式函数关系；基于三维变分数据同化算法构建燃料组件弯曲模型反演
2、为了实现以上目的，本专利技术采取如下的技术方案予以实施：
3、一种融合实测数据的压水堆燃料组件弯曲模型反演方法，其特征在于：包括如下步骤：
4、步骤1：基于堆芯设计建立理论燃料组件模型，对理论燃料组件模型进行随机扰动，获得不同燃料组件弯曲模型；基于压水堆堆芯计算程序对不同的燃料组件弯曲模型进行数值模拟，获得对应的实测数据物理量计算值；
5、步骤2：将步骤1中随机扰动获得的不同燃料组件弯曲模型作为输入，对应的实测数据物理量计算值作为输出，构建数据集；
6、步骤3：基于全连接神经网络对步骤2获得的数据集进行训练，训练完成后基于每一层神经网络的权重矩阵和偏置向量得到燃料组件弯曲模型与实测数据物理量计算值之间的显式函数关系；该步骤具体实施分为以下3个步骤；
7、1)全连接神经网络的输入层和输出层运算分别如公式(1)和公式(2)所示；
8、a(1)＝w(1)x+b(1) 公式(1)
9、y＝w(out)z(p)+b(out) 公式(2)
10、式中：
11、a(1)——第1层隐藏层收到的输入向量
12、w(1)——输入层与第1层隐藏层之间的权重矩阵
13、x——全连接神经网络的输入即燃料组件弯曲模型
14、b(1)——输入层与第1层隐藏层之间的偏置向量
15、y——全连接神经网络的输出即实测数据物理量计算值
16、w(out)——第p层隐藏层与输出层之间的权重矩阵
17、z(p)——第p层隐藏层的输出向量
18、b(out)——第p层隐藏层与输出层之间的偏置向量
19、2)全连接神经网络的第g层隐藏层运算如公式(3)和公式(4)所示；
20、a(g)＝w(g)z(g-1)+b(g) 公式(3)
21、z(g)＝f(g)(a(g)) 公式(4)
22、式中：
23、a(g)——第g层隐藏层收到的输入向量
24、w(g)——第g-1层隐藏层与第g层隐藏层之间的权重矩阵
25、z(g-1)——第g-1层隐藏层的输出向量
26、b(g)——第g-1层隐藏层与第g层隐藏层之间的偏置向量
27、f(g)——第g层隐藏层的激活函数
28、z(g)——第g层隐藏层的输出向量
29、3)训练完成后，得到每一层神经网络的权重矩阵和偏置向量，基于公式(1)至公式(4)获得全连接神经网络输入燃料组件弯曲模型x和输出实测数据物理量计算值y之间的显式函数的关系；
30、步骤4：结合步骤3获得的燃料组件弯曲模型与实测数据物理量计算值之间的显式函数关系，基于三维变分数据同化算法构建燃料组件弯曲模型反演的代价函数；该步骤具体实施分为以下三个步骤；
31、1)采用二阶自回归模型构建燃料组件弯曲模型误差协方差矩阵b，如公式(5)和公式(6)所示；
32、b＝b0c 公式(5)
33、
34、式中：
35、b——燃料组件弯曲模型误差协方差矩阵
36、b0——标量系数
37、c——相关性矩阵
38、ci,j——相关性矩阵第i行、第j列的元素
39、rij——第i个燃料组件和第j个燃料组件的径向距离
40、l——1个燃料组件的径向长度
41、2)假设实测数据之间无相关性，简化实测数据误差协方差矩阵r为对角矩阵，非对角元素均等于0；采用实测数据的误差限值要求作为标准差，实测数据误差协方差矩阵r的对角元素如公式(7)所示；
42、rii＝(αvo,i)2 公式(7)
43、式中：
44、rii——实测数据误差协方差矩阵r第i行、第i列的元素
45、α——实测数据的误差限值要求
46、vo,i——第i个实测数据
47、3)基于三维变分数据同化算法构建燃料组件弯曲模型反演的代价函数，如公式(8)所示；
48、
49、式中：
50、j(x)——燃料组件弯曲模型x反演的代价函数
51、x——燃料组件弯曲模型
52、xb——理论燃料组件模型
53、h——燃料组件弯曲模型与实测数据物理量计算值之间的显式函数关系vo——实测数据
54、r——实测数据误差协方差矩阵
55、步骤5：结合实测数据求解代价函数最小值，得到接近真实状态的燃料组件弯曲模型，完成燃料组件弯曲模型反演。
56、与现有技术相比较，本专利技术具备如下优点：
57、1.本专利技术依赖现有的实测数据即可完成计算，无需额外增加其他数据测量的工作量；
58、2.本专利技术通过全连接神经网络构建燃料组件弯曲模型与实测数据物理量计算值之间的显式函数关系，并基于三维变分数据同化算法将实测数据融合到燃料组件弯曲模型反演中，获得接近真实状态的燃料组件弯曲模型，克服了压水堆运行中真实燃料组件弯曲模型无法获得的困难，进一步提高数值模拟本文档来自技高网...
【技术保护点】
1.一种融合实测数据的压水堆燃料组件弯曲模型反演方法，其特征在于：包括如下步骤：
【技术特征摘要】
1.一种融合实测数据的压水堆燃料组件弯曲...
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。