Research on Fault Detection and Identification Method of Small PWR Based on Principal Component Analysis
-
摘要: 故障检测和辨识对于小型压水堆的安全经济运行具有重要意义。反应堆中通常采用基于信号和专家知识经验的故障检测和辨识方法，操纵员往往不能从海量的故障数据信息中及时准确甄别故障类型，追溯故障原因。本文提出了采用主元分析进行小型压水堆故障检测和辨识的方法。首先利用RELAP5程序对小型压水堆建模，获得典型故障的样本数据。其次，基于主元分析理论对样本降维，并计算T2和Q两个统计量，通过判断是否超出阈值来检测反应堆运行状态。然后，利用贡献率图方法分析了过程变量对于统计量的贡献率，从而确定了对故障特征变化起主要作用的变量，实现对不同故障的辨识。最终和实际物理过程分析结果进行对比，验证了该方法的有效性。Abstract: Fault detection and identification are important for the safety and economy of small PWRs. The fault detection and identification method based on signal and expert knowledge and experience is usually applied in nuclear reactors. However, operators are often unable to identify the fault type and trace the fault cause in time and accurately from the massive fault data information. A method of fault detection and identification of small PWR based on principal component analysis is presented in this paper. First, the model of a small PWR is established by RELAP5 code, and the sample data of typical faults is obtained. Second, the dimension of sample data is reduced by using principle component analysis method. T2 and Q statistics are calculated to detect the reactor operation condition by judging whether the thresholds are exceeded. Then, the contribution rate of process variables to statistics is analyzed by using the contribution rate graph method, so as to determine the variables that play a major role in the change of fault characteristics and realize the identification of different faults. Finally, the effectiveness of the method is verified by comparing with the actual physical process analysis results.
-
图 2 PCA方法故障检测和辨识流程图
Figure 2. Fault Detection and Identification Flow Chart Based on PCA Method
图 3 不同事故下T2和Q统计量计算结果
Figure 3. Calculation Results of T2 and Q under Different Accidents
图 4 不同事故工况主成分对T2统计量贡献率
Figure 4. Contribution Rate of Principal Component to T2 under Different Accidents
图 5 LOCA初期不同变量对T2和Q统计量贡献率分析
Figure 5. Analysis of Contribution Rate of Different Variables to T2 and Q under Initial Stage of LOCA
图 6 落棒事故初期不同变量对T2和Q统计量贡献率
Figure 6. Contribution Rate of Different Variables to T2 and Q under Initial Stage of Rod Drop Accident
图 7 LOCA后期不同变量对T2和Q统计量贡献率分析
Figure 7. Analysis of Contribution Rate of Different Variables to T2 and Q under Final Stage of LOCA
图 8 落棒事故后期不同变量对T2和Q统计量贡献率
Figure 8. Contribution Rate of Different Variables to T2 and Q under Final Stage of Rod Drop Accident
图 10 落棒事故下标准化变量随时间变化
Figure 10. Standardized Variables Vary with Time Under Rod Drop Accident
表 1 反应堆中选取的变量
Table 1. Variables Selected for Reactor
编号 变量名称 单位 1 反应堆功率 MW 2/3 冷却剂流量（环路1和环路2） kg/s 4/5 热段和冷段冷却剂温度（环路1） K 6/7 热段和冷段冷却剂温度（环路2） K 8 稳压器压力 MPa 9 稳压器水位 m 10/11 蒸汽发生器（SG1和SG2）蒸汽流量 kg/s 12/13 SG1和SG2给水流量 kg/s 14/15 SG1和SG2出口蒸汽温度 K 16/17 SG1和SG2水位 m 18/19 SG1和SG2出口蒸汽压力 MPa 
下载: 导出CSV
表 2 特征值及累计方差贡献率
Table 2. Contribution Rate of Eigenvalues and Cumulative Variance
主成分 LOCA 落棒事故 λ σ/% λ σ/% PC1 6.86 36.09 6.12 32.20 PC2 6.01 67.74 4.82 57.58 PC3 2.82 82.58 4.02 78.76 PC4 1.53 90.64 1.33 85.78 PC5 0.56 93.61 1.07 91.42 PC6 0.52 96.34 0.76 95.41 PC7 0.39 98.38 0.46 97.86 PC8 0.16 99.24 0.25 99.15 PC9 0.08 99.65 0.09 99.62 PC10 0.05 99.90 0.04 99.82 PC11~PC19 特征较小，此处不继续列出
下载: 导出CSV
表 3 不同事故的不同时期起主要贡献作用的变量
Table 3. Variables with Main Contribution under Different Periods of Different Accidents
事故类型 事故阶段 起主要贡献作用的变量及其编号 LOCA 前期 冷却剂流量（2和3） 稳压器水位（9） 后期 稳压器压力（8） 稳压器水位（9） 落棒事故 前期 功率（1） 冷却剂流量（2和3） 热段冷却剂温度（4和6） 稳压器压力和水位（8和9） 蒸汽发生器蒸汽流量（10） 后期 功率（1） 冷却剂温度（4，5和6） 稳压器压力和水位（8和9） 蒸汽流量（10） 给水流量（12，13） 蒸汽发生器水位（16）
下载: 导出CSV
-
[1] MINATO A, SEKIMOTO H. Design concept and application of small nuclear power reactor[J]. AIP Conference Proceedings, 2009, 1109(1): 171-176. [2] 张新荣. 基于PCA的连续过程性能监控与故障诊断研究[D]. 无锡: 江南大学, 2008. [3] BENAICHA A, MOUROT G, BENOTHMAN K, et al. Determination of principal component analysis models for sensor fault detection and isolation[J]. International Journal of Control, Automation and Systems, 2013, 11(2): 296-305. doi: 10.1007/s12555-012-0142-x [4] AIT-IZEM T, HARKAT M F, DJEGHABA M, et al. Sensor fault detection based on principal component analysis for interval-valued data[J]. Quality Engineering, 2018, 30(4): 635-647. doi: 10.1080/08982112.2017.1391288 [5] TAO E P, SHEN W H, LIU T L, et al. Fault diagnosis based on PCA for sensors of laboratorial wastewater treatment process[J]. Chemometrics and Intelligent Laboratory Systems, 2013, 128: 49-55. doi: 10.1016/j.chemolab.2013.07.012 [6] KARI T, GAO W S. Power transformer fault diagnosis using FCM and improved PCA[J]. The Journal of Engineering, 2017, 2017(14): 2605-2608. doi: 10.1049/joe.2017.0851 [7] GARCIA-ALVAREZ D, FUENTE M J. A comparative study of fault detection techniques based on principal components analysis (PCA)[J]. Revista Iberoamericana de Automática e Informática Industrial (RIAI), 2011, 8(3): 182-195. [8] 宣暨洋. 基于主元分析和贡献图的微小故障诊断研究[D]. 杭州: 浙江大学, 2015. [9] 陈玉昇,杨燕华,林萌,等. 基于主元分析法的核反应堆关键参数提取研究[J]. 核动力工程,2019, 40(S2): 35-38. [10] 艾鑫,刘永阔,蒋利平,等. 基于iForest-Adaboost的核电厂一回路故障诊断技术研究[J]. 核动力工程,2020, 41(3): 208-213. [11] REIS P A L, COSTA A L, PEREIRA C, et al. Assessment of a RELAP5 model for the IPR-R1 TRIGA research reactor[J]. Annals of Nuclear Energy, 2010, 37(10): 1341-1350. doi: 10.1016/j.anucene.2010.05.013 [12] 温冰清. 基于主元分析的故障检测与诊断研究[D]. 南京: 南京师范大学, 2011. [13] 赵春晖, 王福利. 工业过程运行状态智能监控: 数据驱动方法[M]. 北京: 化学工业出版社, 2019: 32-35. -
计量
- 文章访问数: 210
- HTML全文浏览量: 86
- PDF下载量: 31
- 被引次数: 0
