[發明專利]基于拋物線插值法的燃料燃耗確定方法及裝置

在審

• 申請號: 202011209168.5
• 申請日: 2020-11-03
• 公開（公告）號: CN112380680A
• 公開（公告）日: 2021-02-19
• 發明（設計）人: 郭鳳晨;蘆韡;柴曉明;李慶;劉東;安萍;尹強;盧宗健;曾輝;涂曉蘭;馬永強;鄭勇;賀濤
• 申請（專利權）人: 中國核動力研究設計院
• 主分類號: G06F30/20
• 分類號: G06F30/20;G21C17/104;G06F119/06
• 代理公司: 核工業專利中心 11007
• 代理人: 王朋

1.一種基于拋物線插值法的燃料燃耗確定方法，其特征在于，所述方法包括：

針對多個核素當前燃耗步以及緊鄰當前燃耗步的前兩個燃耗步中的每個燃耗步進行以下操作：

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

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

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

根據確定的所述每個燃耗步末時刻的各核素的微觀反應率的預測值和核子密度的第二預測值，采用拋物線插值法，確定各核素在當前燃耗步末時刻的微觀反應率的修正值；

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

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

2.根據權利要求1所述的方法，其特征在于，根據確定的所述每個燃耗步末時刻的各核素的微觀反應率的預測值和核子密度的第二預測值，采用拋物線插值法，確定各核素在當前燃耗步末時刻的微觀反應率的修正值，包括：

針對每個核素，進行以下操作：

將式一作為該核素的拋物線插值法公式；

通過式二確定a_i、b_i和c_i的值；

根據確定的a_i、b_i、c_i，N_i,n，以及式一，確定各核素在當前燃耗步末時刻的微觀反應率的修正值；

其中，n為所述當前燃耗步的序號，n為自然數，N_i,n為核素i在當前燃耗步的前一個燃耗步末時刻的核子密度，N_i,n+1為核素i在當前燃耗步末時刻的核子密度，N_i,n-1為核素i在當前燃耗步的前二個燃耗步末時刻的核子密度，R_i,n為核素i在當前燃耗步的前一個燃耗步末時刻的微觀反應率，R_i,n+1為核素i在當前燃耗步末時刻的微觀反應率，R_i,n-1為核素i在當前燃耗步的前二個燃耗步末時刻的微觀反應率，R_i為核素i在當前燃耗步末時刻的微觀反應率的修正值。

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

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

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