Health monitoring of a Brayton cycle-based power conversion unit Thesis submitted for the degree Doctor of Philosophy at the Potchefstroom Campus of the North-West University Carel P. du Rand Promoter: Prof. G. van Schoor December 2007 Abstract The next generation nuclear power plants like the Pebble Bed Modular Reactor (PBMR) permit for the design of advanced health monitoring (fault diagnosis) systems to improve safety, system reliability and operational performance. Traditionally, fault diagnosis has been performed by applying limit value checking techniques. Although simple, the inability of these techniques to model parameter dependencies and detect incipient fault behaviour renders them unfavourable. More recent approaches to fault diagnosis can be attributed to the advances in computational intelligence. Data driven methods like artificial neural networks are more widely used when modelling complex nonlinear systems, using only historical plant data. These methods are however dependent on the quality and amount of data used for model development. The key to developing an advanced fault diagnosis system is to adopt an integrated approach for monitoring the different aspects of the total process. Within this context, this goal is realized by presenting a new integrated architecture for sensor fault diagnosis in addition to the enthalpy-entropy graph approach for process fault diagnosis. The integrated architecture for sensor fault diagnosis named SENSE, exploits the strengths of several existing techniques whilst reducing their individual shortcomings. A novel approach for process fault diagnosis is proposed based on the characteristics inherent in the design of the PBMR. Power control by means of an inventory control system and no bypass valve operation facilitates a reference model that remains invariant over the power range. Consequently, the devised reference fault signatures remain static during steady state and transient variations of the normal process. In the thesis, both single and multiple fault conditions are considered during steady state and transient variations of the normal process. It is demonstrated that by applying SENSE, the fused variable estimates are consistent and more accurate than the individual sensor readings. Test cases corresponding to 32 single and multiple fault conditions confirmed that it is possible to use the enthalpy-entropy graph approach for process fault diagnosis. In addition, the proposed fault diagnosis approach is validated through an application to real data from the prototype Pebble Bed Micro Model (PBMM) plant. This application demonstrated that the proposed approach is ideally suited for early detection of faults and greatly reduces the amount of plant data required for model development. Acknowledgements To God with Love I would like to extend my sincere gratitude to everybody that helped me to complete the work presented in the thesis: To my supervisor, Prof. George van Schoor for his persistent guidance and valuable inputs. I would like to thank MTech Industrial (Pty.) Ltd. for the use of the simulation software Flownex. I greatly appreciate the assistance from the PBMR (Pty.) Ltd. team for all the resources they made available to me. Furthermore, I am deeply grateful to Chris Nieuwoudt for all his valuable remarks, suggestions and the long discussions. To all my colleagues at the McTronX Research Group who contributed in various degrees. The financial support of THRIP and MTech Industrial is gratefully acknowledged. Lastly, I would like to thank my family and all my friends for their constant support, encouragement and inspiration. Contents List of figures iv List of tables vii Notation and symbols ix Acronyms x 1. Introduction 1 1.1 Motivation 1.2 Problem description 2 1.3 Thesis objectives 1.4 The diagnostic methodology 3 1.5 Thesis layout 4 1.6 Original contributions 5 1.7 Publications 2. Fault detection and isolation 7 2.1 Introduction 2.2 The scope of health monitoring 8 2.2.1 Fault classification 2.2.2 The health monitoring tasks 2.2.3 The health monitoring method 9 2.3 Sensor fault detection and isolation 10 2.3.1 Independent component analysis 1 2.3.2 Instrumentation and calibration monitoring program 1 2.3.3 Nonlinear partial least squares 12 2.3.4 Multivariate state estimation 2.3.5 Auto-associative neural networks 3 2.3.6 Comparisons and limitations 1 2.3.7 Uncertainty analysis 5 2.3.8 Conclusions 1 i 2.4 Process fault detection and isolation 16 2.4.1 Model-based methods 17 2.4.2 Process history based methods 8 2.4.3 Desirable qualities of the fault diagnosis approach 1 2.4.4 The proposed process fault diagnosis approach 9 2.4.5 Conclusions 20 2.5 Component and system degradation in HTGRs 2 2.6 Summary and conclusions 22 3. The plant model and reference system faults 23 3.1 Introduction 2 3.2 Theoretical analysis of the PBMR plant model 4 3.3 Sensitivity analysis of the PBMR MPS 29 3.3.1 The component and system performance parameters 2 3.3.2 Simplified simulation model of the Brayton cycle 30 3.4 Fault classification in the PBMR MPS 33 3.4.1 Fault class 1: Flow bypass 3.4.2 Fault class 2: Main flow resistance (resistive losses) 34 3.4.3 Fault class 3: Effectiveness or efficiency 3 3.5 The PBMR simulation model 36 3.6 Summary and conclusions 9 4. Sensor fault detection and isolation 40 4.1 Introduction 4 4.2 Sensor malfunctions 1 4.3 Sensor validation 4.3.1 Measurement redundancy 42 4.3.2 Non-temporal parity space analysis 4 4.3.3 Statistical shape analysis 4 4.3.4 Maximum process change 6 4.3.5 Principle component analysis 47 4.4 Sensor fusion 53 4.5 Sensor validation and fusion module architecture 55 4.6 Application of sensor FDI in the PBMM and PBMR 60 4.6.1 Case study 1 6 4.6.2 Case study 2 4 4.7 Summary and conclusions 7 5. Application of traditional fault detection techniques 70 5.1 Introduction 7 5.2 Model-free methods - Limit value checking 7 5.3 Model-based methods - Mathematical models 3 5.3.1 Model development and assumptions 5.3.2 The linear turbine model 74 5.4 Summary and conclusions 8 n 6. Implementation of the h-s graph approach for FDI 84 6.1 Introduction 8 6.2 Enthalpy and entropy - An overview 85 6.3 Attributes and construction of the h-s graph 87 6.4 Effect of faults on the h-s graph 9 6.5 Creating reference fault signatures with the h-s graph 95 6.5.1 The error method 9 6.5.2 The area error method 6 6.6 Fault detection and isolation with the h-s graph approach 99 6.6.1 Noise properties 101 6.6.2 Fault detection 6.6.3 Fault isolation 2 6.6.4 Single fault extraction method 104 6.6.5 The h-s graph and reference fault signatures at different power levels 105 6.7 Process variations and the reference h-s graph 108 6.8 Application of the h-s graph approach in the PBMR MPS 113 6.8.1 Fault conditions during steady state operation of the plant 11 6.8.2 Fault conditions during load following of the plant 116 6.9 Summary and conclusions 120 7. Validation of the h-s graph approach for FDI 122 7.1 Introduction 12 7.2 Differences between the PBMR models and the PBMM 123 7.3 Modelling the PBMM plant in Flownex 126 7.4 The PBMM data used for fault emulation 8 7.5 Simulation of the emulated faults in Flownex 132 7.6 Fault detection in the PBMM 135 7.7 The T-P and h-s models applied to the PBMM 137 7.8 Summary and conclusions 8 8. Conclusions and recommendations 140 8.1 Introduction 14 8.2 Conclusions 1 8.2.1 Component degradation and suitable fault classes 14 8.2.2 Sensor fault diagnosis 14 8.2.3 Process fault diagnosis 2 8.2.4 Validation of the process fault diagnosis approach 144 8.3 Original contributions 145 8.4 Recommendations for future research 14 Bibliography 147 Appendices 153 A. Normal distribution and the central limit theorem 15 B. The temperature-pressure versus enthalpy-entropy graphs 154 B.l The h-s graph shape at different power levels 155 B.2 The T-P and h-s models applied to the PBMR 6 in List of figures Fig. 2.2.1 The health monitoring system 10 Fig. 2.3.1 A simple sensor monitoring system Fig. 2.3.2 The structure of a NLPLS 2 Fig. 2.3.3 The structure of an AANN 4 Fig. 3.2.1 PBMR MPS layout 2 Fig. 3.2.2 Solid model of the PBMR MPS 25 Fig. 3.2.3 Modes and states for normal operation of the PBMR 28 Fig. 3.3.1 A network diagram of a simplified Brayton cycle-based MPS 31 Fig. 3.5.1 A network diagram of the PBMR Flownex model 37 Fig. 4.2.1 Sensor malfunctions 4 Fig. 4.3.1 Probability density function for the normal distribution 45 Fig. 4.3.2 Varying the maximum process change 46 Fig. 4.3.3 Measured and PCA reconstructed values for T and P sensors 50 Fig. 4.3.4 The data segments used for training and testing the PCA models 51 Fig. 4.3.5 The SPE index for the pressure PCA model 5 Fig. 4.3.6 The SPE indices for fault detection together with the reconstructed sensors... 52 Fig. 4.4.1 Fusion algorithm applied to random measurements 55 Fig. 4.5.1 A flow chart illustrating the SENSE architecture 6 Fig. 4.5.2 A reasoning map illustrating the SENSE architecture 57 Fig. 4.5.3 Sensor configurations based on the amount of faulty sensors 58 Fig. 4.5.4 Decision-tree illustrating the expert system reasoning 59 Fig. 4.6.1 Sensor notation utilized in the two case studies 60 Fig. 4.6.2 The SPE index for the temperature PCA model 1 Fig. 4.6.3 The SPE index for the healthy sensor configurations 62 Fig. 4.6.4 The sensor residuals for the eight measurement channels 63 Fig. 4.6.5 Fused estimates for the PBMM temperature measurements 4 Fig. 4.6.6 The PBMR sensor data 65 Fig. 4.6.7 The SPE index for the healthy sensor configuration 6 Fig. 4.6.8 The sensor residuals for the eight measurement channels 66 Fig. 4.6.9 Fused estimates for the PBMR pressure measurements 68 IV Fig. 5.2.1 Plant measurements at different nodes during normal fault free conditions.... 72 Fig. 5.3.1 Turbine efficiency and pressure ratio to non dimensional mass flow 75 Fig. 5.3.2 The input-output measurements of the turbine model 7 Fig. 5.3.3 Inlet pressure transient for the turbine 80 Fig. 5.3.4 Flownex and the linear turbine model response for an inlet pressure transient 81 Fig. 5.3.5 Diagram illustrating the interaction between the individual turbine models.... 82 Fig. 6.2.1 Conversion of enthalpy for a steady flow turbine 86 Fig. 6.3.1 The h-s graphs of a closed Brayton cycle 88 Fig. 6.3.2 Theoretical h-s graph of the Brayton cycle for full as well as reduced power. 88 Fig. 6.3.3 The h-s graph for the PBMR, shown at full and reduced power 89 Fig. 6.4.1 The h-s graphs for normal power operation 91 Fig. 6.4.2 The h-s graphs for normal power operation 2 Fig. 6.4.3 The h-s graphs for normal power operation 3 Fig. 6.4.4 The h-s graphs for normal power operation 94 Fig. 6.5.1 Fault signatures for fault 23 for different fault magnitudes 96 Fig. 6.5.2 The normalized error signatures for the PBMR for normal power operation.. 97 Fig. 6.5.3 Defining the areas between the reference and fault h-s graphs 98 Fig. 6.5.4 Deriving h and s fault signatures with the area error method 9 Fig. 6.5.5 The normalized area error signatures for the PBMR 100 Fig. 6.6.1 Flow diagram of the single fault extraction method 4 Fig. 6.6.2 The extracted h and s error signatures for the two single faults 106 Fig. 6.6.3 The normalized signatures for a 1 % and 10 % change at 40/100 % MCR... 107 Fig. 6.6.4 The reference h-s graph 108 Fig. 6.7.1 The h-s graph for GBPC valve operation at different power levels 109 Fig. 6.7.2 Variation surface for valve operation 110 Fig. 6.7.3 The normalized signatures at 100 % MCR for GBPC valve operation 112 Fig. 6.8.1 The isolated faults during steady state operation of the PBMR 114 Fig. 6.8.2 Load following transient during normal power operation of the PBMR 117 Fig. 6.8.3 Single fault detection during load following 11 Fig. 6.8.4 The normalized signatures for fault 18 during load following 119 Fig. 6.8.5 Multiple fault detection during load following 11 Fig. 7.2.1 Schematic showing the three-shaft MPS of the first PBMR configuration.... 123 Fig. 7.2.2 Gas flow path through the three-shaft MPS 124 Fig. 7.2.3 The theoretical h-s graph for the three-shaft PBMR model 12 Fig. 7.2.4 Simplified schematic showing the PBMM three-shaft MPS 125 Fig. 7.2.5 Solid model of the PBMM 126 Fig. 7.3.1 Schematic diagram showing the PBMM Flownex model 127 Fig. 7.3.2 The practical h-s graphs for the PBMM steady state simulation at 95 kPa.... 128 Fig. 7.4.1 Turbo machinery parameters for bearing failure 130 Fig. 7.4.2 Turbo machinery parameters for thrust test 131 Fig. 7.5.1 The h-s graphs for the steady state PBMM datasets 133 Fig. 7.5.2 Normalized fault signatures for the 6 PBMM faults 4 Fig. 7.6.1 Normal probability plots for measurements in the PBMM 135 Fig. 7.6.2 FII for the emulated faults in the PBMM after fault detection 137 v Fig. 7.7.1 Practical graphs of the Brayton cycle for full and reduced power 138 Fig. B.l.l Theoretical and practical graphs of the Brayton cycle 154 Fig. B.2.1 Normalized error signatures for fault 15 157 vi List of tables Table 3.2.1 Typical parameters of the PBMR MPS for normal operation at full power.. 25 Table 3.3.1 Simplified Flownex model reference values 31 Table 3.3.2 Summary of sensitivity analysis results 2 Table 3.4.1 Summary of faults in the PBMR MPS 5 Table 3.5.1 Summary of the Flownex model results for normal power operation 38 Table 4.4.1 MSE results obtained for the sensor fusion algorithm 55 Table 4.6.1 Summary of PBMM temperature sensor faults: Case study 1 61 Table 4.6.2 MSE results obtained for sensor fusion with no faults present 62 Table 4.6.3 MSE results obtained for sensor fusion with the eight faults present 62 Table 4.6.4 Summary of PBMR pressure sensor faults: Case study 2 65 Table 4.6.5 MSE results obtained for sensor fusion with no faults present 67 Table 4.6.6 MSE results obtained for sensor fusion with the eight faults present 67 Table 5.2.1 Fault alarms for the 25 single faults in the PBMR 72 Table 5.3.1 Summary of turbine model variables 76 Table 5.3.2 Summary of turbine model coefficients 9 Table 6.2.1 Relationships for h and s in open systems 86 Table 6.4.1 Flownex results for fault 23 90 Table 6.6.1 Classification results with the FII for multiple faults 105 Table 6.7.1 Results for the h and s calculations Ill Table 6.8.1 Summary of isolated faults during steady state 114 Table 6.8.2 Average (A) and maximum (M) isolation percentage for 3000 samples .... 115 Table 6.8.3 The isolation (I) and rejection (R) averages of the FII for multiple faults.. 116 Table 6.8.4 Summary of isolated faults (fl|f2) during load following 118 Table 7.3.1 PBMM steady state conditions 126 Table 7.3.2 Steady state values for the two operating points 127 Table 7.3.3 Results for the two operating points 129 Table 7.5.1 Parameters for the two steady state simulations 132 Table 7.5.2 Steady state results for the two PBMM datasets 3 Table 7.5.3 PBMM emulated fault conditions modelled in Flownex 13 vii Table 7.5.4 Flownex results for the six emulated fault conditions 134 Table 7.6.1 Statistical properties of the PBMM measurement noise 136 Table 7.6.2 Alarms during fault detection for the PBMM datasets 13 Table 7.6.3 Average FII for the PBMM fault conditions after fault detection 137 Table 7.7.1 Results for the reference models (PBMM operating at 40 % MCR) 139 Table B.2.1 Temperature and pressure results for the Flownex and reference models.. 156 Table B.2.2 Enthalpy and entropy results for the Flownex and reference models 157 viii Notation and symbols General notation PBMR PBMR (Pty) Ltd. Flownex Flownex® Simulink Simulink® Matlab Matlab® List of symbols I Unit matrix v(t) Noise vector O(T) Correlation matrix E expectation ofx(t) t discrete time T discrete time shift QPT work deliverd by turbine QHPC work absorbed by HPC QLPC work absorbed by LPC QRU heat supplied by reactor T]cyck cycle efficiency rjsy switchyard efficiency r]m mechanical efficiency rjgen generator efficiency Phouse house load power Pgnd grid power mRu reactor mass flow cp constant pressure specific heat R gas constant h enthalpy ix s entropy T temperature P pressure M normal distribution mean a standard deviation J normal distribution variance y normal distribution skew Ax maximum process change T normal distribution confidence value Xf fused estimate Sc sensor configuration matrix S fault signature Pref reference parameter value P measured parameter value Snorm normalized fault signature r Euclidean distance r error vector TN noise vector Sarea Area error signature ¥ covariance matrix AtCC hypothesis threshold A^o pipeline pressure drop £ heat exchanger effectiveness n turbo machinery efficiency AU product of the heat transfer area and the overall heat transfer coefficient PR turbo machinery pressure ratio P parity vector V projection matrix X measurement matrix E residual space matrix P principles components vector c2 °SPE SPE threshold h fault direction vector X Acronyms AANN Auto-associative neural network CBCS Core barrel conditioning system CC Correlation coefficient CFD Computational fluid dynamics CWT Cooling water temperature EPRI Electric power research institute FDI Fault detection and isolation FII Fault isolation index GBPC Gas cycle bypass control valve HP High pressure HPC High pressure compressor HTGR High-temperature gas-cooled reactor HV High voltage ICA Independent component analysis ICMP Instrument and calibration monitoring program ICS Inventory control system LP Low pressure LPC Low pressure compressor MCR Maximum continuous rating MCRI Maximum continuous rating inventory MPS Main power system MSET Multivariate state estimation NLPLS Nonlinear partial least squares NN Neural network NPP Nuclear power plant NRC National regulatory commission OLM On-line monitoring PBMM Pebble bed micro model PBMR Pebble bed modular reactor PCA Principle component analysis PCU Power conversion unit PLS Partial least squares PPB Primary pressure boundary XI PR Pressure ratio PT Power turbine RBP Recuperator bypass valve RMSE Root mean squared error ROT Reactor outlet temperature RSQ r-square statistic RU Reactor unit RUCS Reactor unit conditioning system SENSE Sensor validation and fusion module SPE Square prediction error VM Variance and mean index VRE Variance of the reconstruction error VS Variation surface Xll CHAPTER 1 Introduction This chapter lists the primary objectives of the study together with an overview of the chapters presented. 1.1 Motivation Advanced system diagnostics have been extensively researched the past few years to support nuclear power plant (NPP) utilities in plant supervision. The most important tasks of these diagnostic systems are fault detection and isolation. Even though research shows that these diagnostic systems are essential to prolong the lifespan of the plant, only a few real systems are actually installed in operating units [1], [2]. For the next generation type NPPs, it is expected that these diagnostic systems will become a necessity. From a theoretical point of view, fault diagnosis of nonlinear systems is particularly difficult [2]. In addition, obtaining a sufficient accurate analytical model for complex processes like NPPs could take years. Traditionally, limit value checking techniques have been proven to perform well if the plant operates close to steady state. However, implementing a diagnostic system that performs well only during steady state conditions is not a desirable trait. Another traditional approach to fault diagnosis is signal processing. The difficulties with these techniques are distinguishing between changes in the signal properties due to faults or transient variations of the process. More recent approaches to fault diagnosis can be attributed to the advances in computational intelligence [3]. The methods are however data-driven and dependent on the quality and amount of data used for model development. Acquiring such data for the entire operating range in the next generation NPPs will be very difficult due to economical impacts on normal operation. All these factors motivate the development of a new approach to NPP supervision. The goal is to realize a total health monitoring system that is simplistic, reliable and most important, accurate for different variations of the supervised process. 1 Introduction 1.2 Problem description Given the preceding motivation, the goal of the study is to develop an advanced fault diagnosis approach for NPP supervision. The approach should facilitate a novel sensor and process fault diagnosis technique that functions independently. To address these problems, the following solutions are proposed: • The goal of advanced sensor fault diagnosis is realized by integrating existing techniques in a new way to reduce their individual shortcomings. Measurement redundancy is exploited to allow early detection of instrument drift. • A novel approach to process fault diagnosis is accomplished by developing a new method based on a graphical representation of the supervised process. This technique aims to minimize the amount of monitored variables necessary to quantify the overall health of the system without any knowledge of the mathematical structure of the nonlinear dynamic process. 1.3 Thesis objectives The goal of the thesis is to develop a novel approach to fault diagnosis in a nonlinear high-temperature gas-cooled reactor (HTGR) NPP. To address the shortcomings of current fault diagnosis techniques, the main objectives of the study are: 1. Determine the most relevant mechanisms for component degradation in an HTGR main power system (MPS) and formulate suitable fault classes. 2. Develop and implement a comprehensive fault diagnosis approach for health monitoring in an HTGR MPS. The approach must comprise independent sensor and process fault diagnosis methods. Specifically, the following areas are addressed: 2.1 Propose and implement a novel integrated architecture for sensor fault diagnosis to take advantage of the strengths of existing techniques. For this goal, the independent detection of instrument drift is emphasized. 2.2 Propose and implement a novel approach for process fault diagnosis. The goal is to develop a method that adheres to the strengths of existing techniques without incorporating their general deficiencies. The following desirable qualities must be realized: 2.2.1 Robustness with regard to fault propagation, noise and modelling errors. 2.2.2 Model development and re-training should be simplistic. 2.2.3 Early detection of small faults with incipient time behaviour. 2 Introduction 2.2.4 Supervision of the process during transient variations of the normal process. 2.2.5 Isolation of single faults for multiple fault symptoms. 3. Validate the proposed process fault diagnosis approach through application in the Pebble Bed Micro Model (PBMM). Since there are many ambiguities inherent from directly inducing faults in the real system with regard to control and safety concerns, faults are only simulated. The following constraints are imposed on the simulated HTGR NPP: 1. Since the HTGR NPP is mostly operated at full power, only normal power operation of the plant is investigated. This includes steady state operation and transient variations of the normal process. 2. The number of system faults is limited. Also, critical system faults that cause mode and state transitions of the plant are not applicable (discussed in Chapter 3). These faults are accommodated in the automated plant protection systems. From this constraint, it is concluded that the faults will typically be characterized with incipient time behaviour caused by plant degradation. 1.4 The diagnostic methodology The engineering aspects of the study commence in Chapter 3. Firstly, a simplified model of an HTGR is developed in Flownex® comprising the key MPS components. Through a sensitivity analysis of the model, the most relevant system faults are identified (caused by the component degradation mechanisms). The fault parameters are grouped with regard to cause and effect and the final listing of probable single system faults is summarized. The choices for the fault symptoms are motivated and their importance was confirmed with engineers at PBMR (Pty.) Ltd. Following this, an optimized design of the PBMR that includes the inventory control system (ICS) is modelled in Flownex which serves as the reference NPP. The sensor and process fault diagnosis system is developed using Matlab®, Simulink® and Flownex. By means of a Flownex and Simulink interface, data is collectively transferred between the Flownex and Simulink models. Firstly, random noise with different variance is added to the Flownex measurements, which is then passed through a filter model that infers the appropriate sensor malfunctions on the signals. Next, the signals are evaluated by SENSE (Simulink m-file) and the fused estimates are passed to the process fault diagnosis module (Simulink m-file). After signal transformation, the model residuals are checked for consistency. If a discrepancy is detected, fault signatures are extracted from the residuals and matched to the reference fault database with a statistical classifier. Finally, the relevant information regarding process status and sensor health is collectively displayed. 3 Introduction 1.5 Thesis layout An overview of health monitoring techniques is presented in Chapter 2 together with some basic terminology. For the purpose of advanced sensor fault diagnosis and parameter estimation, two redundant and three non-redundant techniques are investigated. To incorporate techniques that will be readily acceptable to regulatory bodies, the Nuclear Regulatory Commission regulations pertaining to on-line monitoring techniques are examined. For the second part of the chapter, process monitoring methods are discussed with their advantages and shortcomings. Lastly, the overall structure of the proposed process fault diagnosis approach is presented. Chapter 3 describes the general topology of the PBMR MPS and stipulates the relevant operating conditions. Through a sensitivity analysis of a simplified PBMR model, the type and origin of the component performance parameters are identified that are synonymous to the probable fault parameters caused by the component degradation mechanisms. Chapter 4 describes the development of a comprehensive sensor fault diagnosis methodology. The relevant monitoring techniques identified in Chapter 2 are integrated to reduce their individual shortcomings and improve measurement integrity. The fault detection and isolation capabilities of the proposed methodology are demonstrated through application to PBMR and PBMM data. With regard to the latter, real plant data obtained from the PBMM prototype plant is used for the validation. Chapter 5 applies two traditional process fault diagnosis techniques to the PBMR. These methods are based on limit value checking of the monitored variables and mathematical modelling of the plant for the purpose of residual generation. The implementation of the methods in the PBMR highlights their general limitations. Chapter 6 derives the h-s graph approach for process fault diagnosis. Firstly, two analytical techniques are utilized to generate reference fault signatures (with Flownex) for the related fault symptoms; whereafter the correlation among the fault signatures is established. It is demonstrated that each of the different reference fault signatures are highly correlated during transient variations of the normal process with negligible variation. Following this, the fault detection and isolation tasks are developed by means of a statistical hypothesis test and classifier that decide whether a given set of process observations contains any faults. In the presence of multiple fault symptoms, a single fault subtraction procedure is formulated to extract and classify the contributing single faults. To incorporate normal process variations like valve changes into the reference system model, the variation surface is proposed. In the final part of the chapter, the proposed methodology is applied to the PBMR. Fault detection and isolation is demonstrated for both steady state and transient conditions. Chapter 7 validates the h-s graph approach for process fault diagnosis through application in the prototype PBMM plant. Firstly, plant measurements captured during test runs are used to validate the integrity of the FLOWNEX simulation model. Following this, a reference system model together with fault signatures is derived for two emulated fault conditions. Lastly, the h-s graph approach is utilized for process fault diagnosis to identify the emulated fault symptoms in the plant data. 4 Introduction Chapter 8 summarizes the most important conclusions reached in the thesis and documents the original scientific contributions of the study. Recommendations and suggestions for future research are also presented. Appendix A lists the central limit theorem. In the thesis, assumptions regarding the measurement noise are based on the theorem. In Appendix B, the prove for the constant shape of the h-s graph at different power levels (bounded by constraints) is derived. This idea forms the basis for the proposed methodology. Following this, the improvement in model prediction is demonstrated through transformation of the measured variables. 1.6 Original contributions The main scientific contributions of the thesis are summarized as follows: • A novel approach is proposed for process fault diagnosis in an HTGR NPP. Plant supervision is realized with a graphical model-based process model (h-s graph) that remains invariant over the power range. The proposed error and area error methods provide static reference h-s fault signatures that remain invariant to operating point changes, transient variations of the normal process and changes in the fault magnitude. There was no reference found to such an approach for process fault diagnosis in an HTGR NPP. • In addition, a new integrated architecture is proposed for sensor fault diagnosis that forms a comprehensive methodology of existing techniques. In a multi-sensor environment, the unique reasoning structure of this approach produces more accurate and reliable estimates of the sensed variables. A literature survey revealed that this unique and integrated reasoning structure has not been developed for application in an HTGR NPP. 1.7 Publications "Enthalpy-entropy graph approach for the classification of faults in the main power system of a closed Brayton cycle HTGR", Annals of Nuclear Energy, Article in press. Article abstract: An enthalpy-entropy (h-s) graph approach for the classification of faults in a new generation type high temperature gas-cooled reactor (HTGR) is presented. The study is performed on a 165 MW model of the main power system (MPS) of the pebble bed modular reactor (PBMR) that is based on a single closed-loop Brayton thermodynamic cycle. In general, the h-s graph is a useful tool in order to understand and characterise a thermodynamic process. It follows that it could be used to classify system malfunctions from fault patterns (signatures) based on a comparison between actual plant graphs and reference graphs. It is demonstrated that by applying the h-s graph approach, different 5 Introduction fault signatures are derived for the examined fault conditions. The fault conditions that are considered for the MPS are categorized in three fault classes and comprise the main flow bypass of the working fluid, an increase in main flow resistance, and a decrease in component effectiveness or efficiency. The proposed approach is specifically illustrated for four single and two multiple fault conditions during normal power operation of the plant. The simulation of the faults suggests that it is possible to classify all of the examined system malfunctions correctly with the h-s graph approach, using only single reference fault signatures. 6 CHAPTER 2 Fault detection and isolation This chapter gives a comprehensive review of advanced health monitoring techniques used in NPPs. In addition, the mechanisms for component degradation in HTGRs are discussed. 2.1 Introduction Health monitoring is an important component in any large scale engineering plant to improve safety, reliability and overall plant performance. With this in mind, the next generation HTGR NPPs offer more complex challenges for advanced system diagnostics. This chapter gives a summary of the different health monitoring techniques that are currently either implemented or proposed for implementation in NPPs. In section 2.2, the fundamental concepts and basic terminology of health monitoring are introduced together with the total health monitoring framework. The framework includes several different tasks and comprises fault detection, fault isolation and fault identification. Section 2.3 presents an overview of different techniques that are applicable to NPP sensor fault diagnosis and process state estimation. Furthermore, the motivations for the techniques used in the thesis are also discussed. Section 2.4 summarizes some of the most relevant process fault diagnosis techniques. These techniques are mainly model-and process history based methods, each with their own unique strengths and weaknesses. In the final part of this section, the desirable qualities of an advanced fault diagnostic approach are presented, together with the general framework of the proposed approach. In order to determine the specific fault classes in the PBMR, section 2.5 presents an outline of the most relevant mechanisms for component degradation in HTGRs. These mechanisms include component corrosion, erosion, fouling and leakage. A summary of this chapter is presented in section 2.6. 7 Fault detection and isolation 2.2 The scope of health monitoring In modern NPPs, information about the current health of the system is essential to improve plant safety and operational levels [4]. Therefore, it is important to detect component faults and irregular system operation promptly. 2.2.1 Fault classification In general, unpermitted deviations from the normal behavior of the components or process are termed faults or failures. Faults are caused by physical defects or imperfections that occur within the component, whilst a failure suggest complete breakdown of the component [5]. The faults that are applicable for this investigation can be divided into the following categories [6]: • Additive process faults: These faults are caused by unknown inputs acting on the plant, which results in a shift in the plant outputs, independent of the measured inputs. A plant leak is a typical example of an additive fault. • Multiplicative process faults: These faults result in system parameter changes, where the outputs are dependent on the magnitude of the inputs. Such faults are mostly associated with component degradation and include fouling and efficiency changes. • Sensor faults: Any discrepancies between the measured and the expected values of the process variables are considered to be sensor malfunctions. • Actuator faults: These faults are described by discrepancies between the intended control actions and the actual realization of these commands by the actuators. 2.2.2 The health monitoring tasks In order to identify and characterize the faults, the health monitoring system should comprise the following tasks [6]: • Fault detection: The identification of an irregularity in the monitored system. • Fault isolation!classification: The origin and the type of fault are determined. • Fault identification: The magnitude of the fault is established. The fault identification phase generally does not justify the additional computation it requires, and for this reason, most monitoring systems only comprise the fault detection and isolation (FDI) phases [6]. For the proposed diagnostic system, real-time computational complexity is reduced by implementing only the fault detection and 8 Fault detection and isolation isolation tasks. The drawback from this restriction is that the fault magnitude will not be directly calculated, but rather estimated from the fault signature magnitude before normalization. With the help of early fault detection and accurate fault isolation, process and component malfunctions can therefore be identified at an early stage to reduce the risk of sudden failure as well as allow enough time for maintenance or repair. Given the complex and safety critical nature of NPPs, the advanced FDI tasks should adhere to the following requirements [7]: • Reduce the occurrence of false alarms during operation due to normal transient variations of the process. • Original fault detection in the event of multiple fault conditions and propagation across subsystems. • Early detection of small faults with abrupt or incipient time behavior. • Reducing misdiagnosed faults due to modelling uncertainties and noise. 2.2.3 The health monitoring method Designing a system for advanced fault diagnosis is a challenging engineering task, particularly in fields related to nuclear processes, owing to the stringent safety and environmental regulations. For these reasons, it is important that the diagnostic method meet the following performance requirements [6]: • Fault detection: - Fault sensitivity: The method must detect incipient faults with a small magnitude. - Detection time: The method must be able to detect faults with the smallest time delay after induction. - Robustness: The method must be able to function in the presence of noise, modelling uncertainties and disturbances, with minimal false alarms. • Fault isolation: The method must be able to distinguish between the different types of faults (single or multiple simultaneous faults) in the presence of noise and modelling errors. It is important to note that some faults, single or multiple, might be non-isolable, since their influence on the system is undistinguishable [6]. The health monitoring system, which constitutes sensor and process FDI, is illustrated in Fig. 2.2.1. 9 Fault detection and isolation System information: Control information: Plant System inputs raw data 3.L.. Sensor FDI Validation and flision II 3i Process FDI features Advanced diagnosis Monitoring and diagnosis Alarms Fault location Fault cause Corrective actions System outputs Fig. 2.2.1 The health monitoring system. 2.3 Sensor fault detection and isolation This section reviews current sensor fault diagnosis techniques ranging from basic, well established methods to the latest reported advanced strategies. Throughout the literature, various methods are proposed by the industry and academia [8] - [15]. Advanced in this context signifies methods that will allow nuclear power utilities to diverge from utilizing a periodic based maintenance approach to condition based strategies. In general, these advanced methods aim at describing the sensors health whilst the plant is operational. A simple block diagram illustrating a sensor monitoring system is depicted in Fig. 2.3.1. Based on research applications, the most relevant techniques used in NPPs include the parity space method, principle component analysis (PCA), independent component analysis (ICA), instrumentation and calibration monitoring program (ICMP), nonlinear partial least squares (NLPLS), multivariate state estimation (MSET) and auto-associative neural networks (AANN) [16]. These techniques can roughly be classified into two categories: techniques that model a redundant group of sensors to obtain the estimate (first four) and models that include non-redundant measurements that are correlated, but not redundant (last three). A study conducted by [16] concluded that the simplicity of redundant techniques and the tractability of their uncertainty calculations could favour them for acceptance by regulatory bodies. For this reason, the non-temporal parity space and PCA techniques are adopted for sensor fault diagnosis based on their relative simplicity and individual strong points (discussed in Chapter 4). ? ? Model Comparison Decision h t Xl —? : Xn Sensed inputs Model estimates Residuals Status Fig. 2.3.1 A simple sensor monitoring system [16]. 10 Fault detection and isolation The following sections discuss and compare the applicability of the remaining techniques for sensor fault diagnosis and highlight their general limitations. In this thesis, it is important to note that only a general overview of each technique is presented together with the basic notation. For the complete mathematical formulations, the reader is referred to the literature referenced. 2.3.1 Independent component analysis In ICA, the sensed variables are described by means of a linear transformation of independent components that is maximally non-Gaussian (normally) distributed [16]. An important characteristic of this technique is its ability to separate the true signal (includes process noise) from the independent measurement noise [17]. This makes ICA a notable candidate for signal pre-processing and filtering. The ICA model is given by (2.3.1) X = AS (2.3.1) with X the observed data of n samples from m sensors, S is the matrix of m independent components and A the mixing matrix. The linear transformation of the observed data into non-Gaussian distributed components Y is Y = WX (2.3.2) where W is the weight matrix. The parameter estimate, which denotes the true signal with process noise, is then given by one of the independent components. 2.3.2 Instrumentation and calibration monitoring program The ICMP was developed by the Electric Power Research Institute (EPRI). The ICMP algorithm is based on a weighted average of each sensor, which is denoted by consistency values c, [16]. The consistency values signify how much a sensor reading contributes to the final estimate. If the value correlates within the defined limits to another, they are consistent. The consistency value c, of the z'-th sensor is calculated with \x, -Xj\<^ + dt, then c, = c, +1 (2.3.3) with x» Xj the values of the z-th and j-th sensors and d,-, d j the consistency check allowance for sensor z and j respectively. The values are checked for consistency in an iterative way against the remaining sensors. The ICMP parameter estimated is calculated with n x = ^-n (2.3.4) 11 Fault detection and isolation where Wi is the weight related to the z'-th sensor. The influence of more reliable sensors within a redundant group can therefore be increased by varying their weights. If there is no preference between the sensors, the weights are set to 1. Following this, the performance of each sensor is determined in relation to x through an acceptance criterion |*-*,|