[發(fā)明專利]基于非線性最小二乘法的燃料燃耗確定方法及裝置

在審

• 申請?zhí)枺?/td>: 202011210059.5
• 申請日: 2020-11-03
• 公開（公告）號: CN112347632A
• 公開（公告）日: 2021-02-09
• 發(fā)明（設(shè)計）人: 郭鳳晨;蘆韡;柴曉明;李慶;劉東;安萍;涂曉蘭;盧宗健;曾輝;馬永強(qiáng);潘俊杰;湯琪芬;尹強(qiáng);陳定勇;馬黨偉
• 申請（專利權(quán)）人: 中國核動力研究設(shè)計院
• 主分類號: G06F30/20
• 分類號: G06F30/20
• 代理公司: 核工業(yè)專利中心 11007
• 代理人: 王朋

1.一種基于非線性最小二乘法的燃料燃耗確定方法，其特征在于，所述方法包括：

針對多個核素從起始時刻至當(dāng)前燃耗步的每個燃耗步進(jìn)行以下操作：

根據(jù)各核素在該燃耗步初時刻的微觀反應(yīng)率的預(yù)測值，以及各核素在該燃耗步起始時刻的核子密度，采用燃耗方程，確定各核素在該燃耗步末時刻的核子密度的第一預(yù)測值；

根據(jù)各核素在該燃耗步末時刻的核子密度的第一預(yù)測值，采用中子輸運方程，確定各核素在該燃耗步末時刻的微觀反應(yīng)率的預(yù)測值；

根據(jù)各核素在該燃耗步末時刻的微觀反應(yīng)率的預(yù)測值，以及各核素在該燃耗步起始時刻的核子密度，采用燃耗方程，確定各核素在該燃耗步末時刻的核子密度的第二預(yù)測值；

根據(jù)確定的所述每個燃耗步末時刻的各核素的微觀反應(yīng)率的預(yù)測值和核子密度的第二預(yù)測值，采用非線性最小二乘法，確定各核素在當(dāng)前燃耗步末時刻的微觀反應(yīng)率的修正值；

根據(jù)各核素在當(dāng)前燃耗步末時刻的微觀反應(yīng)率的修正值，各核素在起始時刻的核子密度，采用燃耗方程，確定各核素在當(dāng)前燃耗步末時刻的核子密度的修正值；

根據(jù)各核素在當(dāng)前燃耗步末時刻的核子密度的修正值，以及各核素在當(dāng)前燃耗步末時刻的核子密度的第一預(yù)測值，確定各核素在當(dāng)前燃耗步末時刻的核子密度的最終值。

2.根據(jù)權(quán)利要求1所述的方法，其特征在于，根據(jù)確定的所述每個燃耗步末時刻的各核素的微觀反應(yīng)率的預(yù)測值和核子密度的第二預(yù)測值，采用非線性最小二乘法，確定各核素在當(dāng)前燃耗步末時刻的微觀反應(yīng)率的修正值，包括：

針對每個核素，進(jìn)行以下操作：

將式一作為該核素的非線性最小二乘法公式；

R(x_i)＝a₀+a₁x_i+a₂x_i²+…+a_kx_i^k 式一

通過式二確定a₀,a₁,…a_k的值；

根據(jù)確定的a₀,a₁,…a_k，x_i，以及式一，確定各核素在當(dāng)前燃耗步末時刻的微觀反應(yīng)率的修正值；

其中，i為燃耗步的序號，n為所述多個燃耗步的數(shù)量，1≤i≤n，i為自然數(shù)，k為式一中系數(shù)的序號，k為大于2的正整數(shù)，x_i為該核素在第i燃耗步末時刻的核子密度的第二預(yù)測值，y_i為該核素在第i燃耗步末時刻的微觀反應(yīng)率的預(yù)測值，R(x_i)為該核素在當(dāng)前燃耗步末時刻的微觀反應(yīng)率的修正值。

3.根據(jù)權(quán)利要求1所述的方法，其特征在于，根據(jù)各核素在當(dāng)前燃耗步末時刻的核子密度的修正值，以及各核素在當(dāng)前燃耗步末時刻的核子密度的第一預(yù)測值，確定各核素在當(dāng)前燃耗步末時刻的核子密度的最終值，包括：

將每個核素在當(dāng)前燃耗步末時刻的核子密度的修正值，以及該核素在當(dāng)前燃耗步末時刻的核子密度的第一預(yù)測值的平均值，作為該核素在當(dāng)前燃耗步末時刻的核子密度的最終值。

4.根據(jù)權(quán)利要求1所述的方法，其特征在于，所述多個核素包括：¹⁵²Gd、¹⁵⁴Gd、¹⁵⁵Gd、¹⁵⁶Gd、¹⁵⁷Gd、¹⁵⁸Gd、¹⁶⁰Gd。