[发明专利]一种用于压水堆系统动态模型参数评价的灵敏度分析方法

权利要求书

1.一种用于压水堆系统动态模型参数评价的灵敏度分析方法，其特征在于：包括如下步骤，

步骤S1：针对压水堆系统动态模型，按照其内部物理边界分解为多个子模块模型；该些子模块包括：堆芯中子动态模块、堆芯燃料及冷却剂温度模块、热线温度模块、冷线温度模块、一回路平均温度模块、蒸汽发生器模块、反应堆功率控制系统模块和冷却剂主泵模块；

步骤S2：考虑模型参数之间的相关性，确定各子模块变量集、参数集；

步骤S3：进行子模块变量初始化，计算轨迹灵敏度指标，从而分析确定影响输出变量稳态值的全部参数；

步骤S4：评估得到参数的可辨识难易程度，以利于模型参数的辨识及其准确性评价；

所述步骤S2中，各子模块变量集、参数集的确定方法如下：

对于堆芯中子动态模块，其变量集为

{N_r,C_r,T_F,T_av,ρ_ext}

其中，N_r为中子通量密度；C_r为等效单组缓发中子先驱核密度，为中间变量；T_F为堆芯燃料温度；T_av为反应堆冷却剂平均温度；ρ_ext为控制棒引入的反应性；

对于堆芯中子动态模块，其参数集为

A＝{l,β,λ,α_F,α_C}

其中，l为平均中子寿命；β为缓发中子组的总份额；λ为等效缓发中子组的延时常数；α_F为燃料温度反应性系数；α_C为冷却剂温度反应性系数；

对于堆芯燃料及冷却剂温度模块，其变量集为

{T_F,T_av,T_θ1,T_θ2}

其中，T_θ1为反应堆冷却剂入口温度；T_θ2为反应堆冷却剂出口温度；

对于堆芯燃料及冷却剂温度模块，其参数集为

B＝{b₁,b₂,b₃,b₄,b₅}

b₁＝P₀×F_f/μ_f

b₂＝P₀×(1-F_f)/μ_c

b₃＝Ω/μ_f

b₄＝Ω/μ_c

b₅＝M/μ_c

其中，P₀为堆芯热功率，F_f为燃料发热份额，Ω为堆芯中燃料与冷却剂传热系数，μ_f为堆芯中燃料比热容，μ_c为堆芯冷却剂比热容，M＝D_sp×C_pc×m_Cn，D_sp为冷却剂主泵流量，C_pc为冷却剂比热容，m_Cn为额定工况下冷却剂质量流量；

对于蒸汽发生器模块，其变量集为

{T_p,P_s,T_m,T_HL,T_CL,Q_s}

其中，T_p为一回路冷却剂平均温度；P_s为蒸汽发生器出口蒸汽压力；T_m为U形传热管温度；T_HL为蒸汽发生器一回路冷却剂入口温度；T_CL为蒸汽发生器一回路冷却剂出口温度；Q_s为二回路蒸汽流量；

对于蒸汽发生器模块，其参数集为

C＝{c₁,c₂,c₃,c₄,c₅,c₆}

c₁＝M/μ_p

c₂＝Ω_p/μ_p

c₃＝Ω_p/μ_m

c₄＝Ω_S/μ_m

c₅＝Ω_S

c₆＝(h_s-h_fw)×G_sn

其中，Ω_p为蒸汽发生器中一回路冷却剂与U形传热管传热系数，Ω_S为U形传热管与二回路蒸汽传热系数，μ_p为蒸汽发生器冷却剂比热容，μ_m为U形传热管比热容，G_sn为额定工况下二回路主蒸汽流量，h_fw和h_s分别为二回路给水入口温度比焓和出口蒸汽比焓。