Research on Gas-Path Fault-Diagnosis Method of Marine Gas Turbine Based on Exergy Loss and Probabilistic Neural Network

: In order to improve the accuracy of gas-path fault detection and isolation for a marine three-shaft gas turbine, a gas-path fault diagnosis method based on exergy loss and a probabilistic neural network (PNN) is proposed. On the basis of the second law of thermodynamics, the exergy ﬂow among the subsystems and the external environment is analyzed, and the exergy model of a marine gas turbine is established. The exergy loss of a marine gas turbine under the healthy condition and typical gas-path faulty condition is analyzed, and the relative change of exergy loss is used as the input of the PNN to detect the gas-path malfunction and locate the faulty component. The simulation case study was conducted based on a three-shaft marine gas turbine with typical gas-path faults. Several results show that the proposed diagnosis method can accurately detect the fault and locate the malfunction component.


Introduction
A marine gas turbine operates under hostile ocean environments.The air contains salt-ingested particles, which will have an impact on the gas-path components such as the compressor, combustion chamber, and turbine, and can lead to fouling, erosion, and corrosion [1,2].These faults will change the structure of the components and cause performance degradation, reducing the safety and stability of the gas turbine [3][4][5].The diagnosis of gas-path faults is becoming a major issue [6,7].
Physical failure can be reflected by changes in efficiency and flow of the components, and, in turn, causes changes in gas-path thermal parameters of a gas turbine, such as pressure, temperature, rotational speed, and fuel flow rate.Such relationships were described by Urban who proposes the linear gas-path analysis method [8].In recent years, many gas-path fault diagnosis methods have been proposed to assess gas turbine health status, such as non-linear gas-path analysis [9][10][11], the rule-based fuzzy expert system [2,12,13], Bayesian hierarchical models [14,15], neural networks [16][17][18][19], the genetic algorithm [4,20,21], the multiple-model method [3,22] and exergy analysis [23].Gas-path analysis is an inversely mathematical problem to obtain the deviation of component performance parameters over gas-path measured variables.Most case studies show that artificial intelligence methods, such as artificial neural networks, Bayesian hierarchical models, and fuzzy expert systems, may effectively isolate the faulty components but may not assess the severity of the fault easily [24].
The traditional gas-path analysis method is based on the first law of thermodynamics which may lead to the strong coupling effect when the number of components involved in fault diagnosis is large, and the malfunction may not be effectively identified [10].Exergy analysis is a new branch

Exergy Flow Analysis of Gas Turbine
The general layout of the three-shaft marine gas turbine is shown in Figure 1, and gas-path components include a low-pressure compressor, a high-pressure compressor, a combustion chamber, a high-pressure turbine, a low-pressure turbine, and a power turbine.The marine gas turbine is driven by a low-pressure turbine, and compressed air is then further compressed by a high-pressure compressor.The compressed air enters the combustion chamber, expands, and burns with the fuel.The mixture then expands into the high-and low-pressure turbines, and finally enters the power turbine, where power and torque are output, and the exergy flow of the gas turbine is shown in Figure 2.
There are 11 exergy flows among the components and the external environment, regardless of the generator connected to the power turbine.Flow 0 (atmosphere) is compressed by the low-pressure compressor, generating flow 1 (compressed air).Flow 1 is compressed by a high-pressure compressor, generating flow 2 (compressed air).Flow 2 and flow 10 (fuel) are mixed and burned in the combustion chamber to form flow 3 (gas) and expand into a high-pressure turbine to do work.The energy flowing out of the high-pressure turbine subsystem is split into flow 8 and flow 4. Flow 8 drives the high-pressure compressor to rotate, and flow 4 enters the low-pressure turbine to continue to expand and do work.The energy outflow from the high-pressure turbine subsystem is divided into flow 9 and flow 5. Flow 9 drives the low-pressure compressor to rotate, and flow 5 expands into the power turbine to do work.Flow 6 enters the generator to generate energy flow 7 (electrical energy).Flow 11 is exhaust gas.The operation of a gas turbine is accompanied by exergy flow among components and the environment, and the process of exergy flow will produce irreversible loss, which we call exergy loss.The exergy flow and exergy loss of each component are shown in Table 1.The operation of a gas turbine is accompanied by exergy flow among components and the environment, and the process of exergy flow will produce irreversible loss, which we call exergy loss.The exergy flow and exergy loss of each component are shown in Table 1.
Combustion chamber High-pressure turbine The operation of a gas turbine is accompanied by exergy flow among components and the environment, and the process of exergy flow will produce irreversible loss, which we call exergy loss.The exergy flow and exergy loss of each component are shown in Table 1.

Exergy Models for Gas Turbine
In order to derive the loss of gas turbine systems, it is necessary to establish exergy models of the high-pressure compressor, low-pressure compressor, combustion chamber, high-pressure turbine, low-pressure turbine, and power turbine.
Enthalpy and entropy are the key parameters for calculating exergy and are related to temperature.The enthalpy and entropy of air can be obtained from a table of thermodynamic properties of air.The specific enthalpy and entropy of air can be obtained by the air temperature.T 0 is standard temperature, s 0 is standard-specific entropy, and h 0 is standard-specific enthalpy.
(1) Exergy model of the low-pressure compressor The formula for calculating the specific exergy of energy flow at the outlet of the low-pressure compressor of a marine gas turbine B LCout is as follows: where h LCout is the specific enthalpy of the outlet of the low-pressure compressor, kJ/kg; s LCout is the specific entropy of isentropic variation at the outlet of the low-pressure compressor, kJ/(kg•K); and B 1 is the specific exergy of the energy of compressed air flowing out of the low-pressure compressor and entering the high-pressure compressor.
The formula for calculating the specific exergy of energy flow at the inlet of the low-pressure compressor of a marine gas turbine B LCin is as follows: where B 0 is the specific exergy of air flowing into the low-pressure compressor, B 0 = 0; B 9 is the specific exergy of energy flow (mechanical energy) flowing out of the low-pressure turbine subsystem and driving the rotation of the low-pressure compressor: where J LCLT is the rotary inertia of the shaft connecting the low-pressure compressor and the low-pressure turbine; and ω 1 is the rotation rate of the shaft connecting the low-pressure compressor and the low-pressure turbine.
The exergy loss of the low-pressure compressor I LC is expressed as Equation (4): (2) Exergy model of the high-pressure compressor The specific exergy of energy flow at the outlet of the high-pressure compressor of a marine gas turbine B HCout is calculated by Equation ( 5): where h HCout is the specific enthalpy of the outlet of the high-pressure compressor, kJ/kg; s HCout is the specific entropy of isentropic variation at the outlet of the high-pressure compressor, kJ/(kg•K); and B 2 is the specific exergy of the energy of compressed air flowing out of the high-pressure compressor and entering the combustion chamber.Equation ( 6) is used to calculate the specific exergy of energy flow at the inlet of the high-pressure compressor for marine gas turbine B HCin : Energies 2019, 12, 4701 5 of 17 where B 8 is the specific exergy of the energy flow (mechanical energy) flowing out of the high-pressure turbine subsystem and driving the high-pressure compressor to rotate: where J HCHT is the rotary inertia of the shaft connecting the high-pressure compressor and the high-pressure turbine; and ω 2 is the rotation rate of the shaft connecting the high-pressure compressor and the high-pressure turbine.
Exergy loss of the high-pressure compressor is expressed as Equation ( 8): (3) Exergy model of the combustion chamber Equation ( 9) is used to calculate the specific exergy of energy flow at the outlet of the combustion chamber for a marine gas turbine B CCout : where h CCout is the specific enthalpy at the outlet of the combustion chamber, kJ/kg; s CCout is the specific entropy at the outlet of the combustion chamber, kJ/(kg•K); and B 3 is the specific exergy flow 9, as shown in Figure 2. Equation ( 10) is used to calculate the specific exergy of energy flow at the inlet of the combustion chamber for a marine gas turbine B CCin : where B 10 is the specific exergy of the energy of the fuel entering the combustion chamber: where ∆H f is the specific enthalpy of fuel in the combustion chamber, kJ/kg; f is the fuel-air ratio; G f is the fuel flow, and G Bin is the air flow.The exergy loss of the combustion chamber I CC is expressed as Equation ( 13): (4) Exergy model of the high-pressure turbine Equation ( 14) is used to calculate the specific exergy of energy flow at the outlet of the high-pressure turbine for a marine gas turbine B HTout : where h HTout is the outlet specific enthalpy of the high-pressure turbine, kJ/kg; s HTout is the specific entropy of isentropic variation at the outlet of the high-pressure turbine, kJ/(kg•K); and B 4 is the energy flowing out of the high-pressure turbine entering the low-pressure turbine.
Equation ( 16) is used to calculate the specific exergy of energy flow at the inlet of the high-pressure turbine for a marine gas turbine B HTin : Exergy loss of high-pressure turbine I HT is expressed as Equation ( 17): (5) Exergy model of the low-pressure turbine Equation ( 18) is used to calculate the specific exergy of energy flow at the outlet of the low-pressure turbine of a marine gas turbine B LTout : where h LTout is the outlet specific enthalpy of the low-pressure turbine, kJ/kg; s LTout is the specific entropy of isentropic variation at the outlet of the low-pressure turbine, and B 5 is the specific exergy of the energy flowing out of the high-pressure turbine expanding into the power turbine to do work.Equation ( 20) is used to calculate the specific exergy of energy flow at the inlet of the low-pressure turbine for marine gas turbine B LTin : Exergy loss of the low-pressure turbine I LT is expressed as Equation ( 21): (6) Exergy model of the power turbine A power turbine is a kind of impeller used to drive an external load (propeller, motor, etc.) and convert energy into mechanical work.The gas at the outlet of the low-pressure turbine expands for the last time in the power turbine, outputs external mechanical power, and emits exhaust gas.In this paper, no waste heat boiler is considered, and the exhaust gas is discharged to the atmosphere.
Equation ( 22) is used to calculate the specific exergy of the outlet energy of the power turbine: where B 6 is the specific exergy of energy entering the generator; and W net is the work done by the power turbine.Equation ( 23) is used to calculate the specific exergy of the inlet energy of the power turbine: The exergy loss of the power turbine I PT is expressed as Equation ( 24): The marine gas turbine exergy model can generate large amounts of fault simulation data, which will extract the exergy loss feature of the component under the healthy condition and typical faulty condition.

Gas-Path Fault-Diagnosis Approach Based on Exergy Loss and Probabilistic Neural Network (PNN)
When the gas turbine is in fault status, the exergy loss of each component will change, and its range of change can be used to assess the severity of the fault.However, the outlet temperature of the combustor does not install the temperature sensor for high temperature and costly, B 3 (flow 3) can not be calculated.Therefore, the proposed exergy model in Section 2 can only get the exergy of the low-pressure compressor, the high-pressure compressor, the low-pressure turbine, and the power turbine by Equations ( 4), ( 8), ( 21) and ( 24), To detect and isolate the fault location accurately under different operating conditions, this paper introduces the probabilistic neural network (PNN) classifier in the traditional exergy analysis method.

Probabilistic Neural Network (PNN)
A PNN is a kind of neural network based on statistical principles, which is commonly employed to solve problems of pattern classification [30].The PNN is an efficient and robust classifier, obtained when the Bayes strategy for decision making is combined with a non-parametric estimator for probability density functions [31].Unlike a traditional multilayer feed-forward network, which requires a back-propagation algorithm for back-propagation computation, it is a fully forward computational process.It can be trained very easily and has high classification accuracy [12,32].The fault diagnosis method based on PNN is a widely accepted decision method in probability statistics, which can be described as follows [33]: assuming there are two known fault modes θ A , and θ B , for the test sample X = [X 1 , X 2 , ..., X s ] to be judged, then: where h A , h B is the prior probability of the fault mode θ A , θ B , respectively (h A = N A /N, h B = N B /N); N A , N B is the number of training samples for fault mode θ A , θ B , respectively; and N is the total number of training samples; l A and l B are the loss functions associated respectively with the decision X ∈ θ B and X ∈ θ A ; f A (X) and f B (X) are respectively the probability density functions of the category A and B.
The PNN is composed of an input layer, a pattern layer, a summation layer, and an output layer, as shown in Figure 3 [34].
The input layer does not perform any computation and simply distributes the input to the neurons in the pattern layer.The number of neurons is equal to the dimension of the sample vector.
The task of the pattern layer is to calculate the pattern, matching the relationship between the testing samples and the training samples, and to centralize the categories with high similarity.The number of neurons in the pattern layer is equal to the sum of the number of trained samples in all categories, and the neuron X ij computes the pattern layer output as follows [30]: where X = [X 1 , X 2 , ..., X S ] is the test sample with s dimension; s is the dimension of the input vector X; X ij is the jth training sample of the ith class; and σ = [0, 1] is a smoothing factor.The summation layer adds the cumulative probability of the class to obtain the estimated density function of the fault patterm.Each summation node receives the outputs from pattern nodes associated with a given class [35]: where N i is the number of samples in the ith class (i = 1, 2, ..., n); n is the number of classes; weights w i = (w i0 , w i1 , ..., w iNi ) are created during the training process; and P i (X) ∈ (−1, 1).The output layer receives all kinds of density functions.The class of an unknown pattern X is obtained in the output layer using Equation (27): The output value of the category with the largest function value is 1, and the output value of the remaining categories is 0.
requires a back-propagation algorithm for back-propagation computation, it is a fully forward computational process.It can be trained very easily and has high classification accuracy [12,32].The fault diagnosis method based on PNN is a widely accepted decision method in probability statistics, which can be described as follows [33]: assuming there are two known fault modes θA, and θB, for the test sample X = [X1, X2, ..., Xs] to be judged, then: where hA, hB is the prior probability of the fault mode θA, θB, respectively (hA = NA/N, hB = NB/N); NA, NB is the number of training samples for fault mode θA, θB, respectively; and N is the total number of training samples; lA and lB are the loss functions associated respectively with the decision The PNN is composed of an input layer, a pattern layer, a summation layer, and an output layer, as shown in Figure 3 [34].The input layer does not perform any computation and simply distributes the input to the neurons in the pattern layer.The number of neurons is equal to the dimension of the sample vector.
The task of the pattern layer is to calculate the pattern, matching the relationship between the testing samples and the training samples, and to centralize the categories with high similarity.The

Fault Diagnosis Process
The fault-diagnosis process is shown in Figure 4.
Energies 2019, 12, x FOR PEER REVIEW 8 of 17 number of neurons in the pattern layer is equal to the sum of the number of trained samples in all categories, and the neuron Xij computes the pattern layer output as follows [30]: where X = [X1, X2, ..., XS] is the test sample with s dimension; s is the dimension of the input vector X; Xij is the jth training sample of the ith class; and σ = [0, 1] is a smoothing factor.
The summation layer adds the cumulative probability of the class to obtain the estimated density function of the fault patterm.Each summation node receives the outputs from pattern nodes associated with a given class [35]: where Ni is the number of samples in the ith class (i = 1, 2, ..., n); n is the number of classes; weights wi = (wi0, wi1, ..., wiNi) are created during the training process; and The output layer receives all kinds of density functions.The class of an unknown pattern X is obtained in the output layer using Equation ( 27): The output value of the category with the largest function value is 1, and the output value of the remaining categories is 0.

Fault Diagnosis Process
The fault-diagnosis process is shown in Figure 4.  Step 1: Data acquisition The measured variables that acquire from the gas turbine monitoring and control system include the outlet temperatures of four components, boundary condition (T0, P0), and fule flow wf.The outlet temperatures of the four components make up the set T = {T1, T2, T3, T4}.T1 is the outlet temperature of the low-pressure compressor, T2 is the outlet temperature of the high-pressure compressor, T3 is

Step 1: Data acquisition
The measured variables that acquire from the gas turbine monitoring and control system include the outlet temperatures of four components, boundary condition (T 0 , P 0 ), and fule flow w f .The outlet temperatures of the four components make up the set T = {T 1 , T 2 , T 3 , T 4 }.T 1 is the outlet temperature of the low-pressure compressor, T 2 is the outlet temperature of the high-pressure compressor, T 3 is the outlet temperature of the high-pressure turbine, and T 4 is the outlet temperature of the low-pressure turbine.
Step 2: Parameter estimation under the healthy condition Parameters estimation under healthy condition is on the healthy condition estimation model of a gas turbine that is established by back-propagation neural networks [36].For the gas turbine health estimation model, the inputs are the environment parameters (T 0 , P 0 ) and fuel flow control variable w f , and the outputs are temperature set T h on the healthy condition.
Step 3: Exergy-loss calculation T and T h are input into the exergy model, respectively.The output of the exergy model is the exergy loss of gas turbine I = {I 1 , I 2 , ..., I s }, and exergy loss of its healthy condition estimation model

Step 4: Data normalization
To eliminate the magnitude differences among different input parameters, the relative change of exergy loss The data of one component of the gas turbine in the training samples are gathered together to form a one-dimensional array ∆I m (m = 1, 2, ..., s), where s is the number of components.The maximum ∆I m max and minimum ∆I m min of the array ∆I m are selected, and normalization of the one-dimensional array is carried out according to Equation (28): Step 5: Fault detection and isolation The input X of the PNN that has been trained by the training samples which select from the fault database, and the fault mode C i (i = 0, 1, 2, ..., n) will be output from the output layer of the PNN.Finally, the fault diagnosis results are obtained by consulting the fault discrimination matrix is established for fault detection and isolation and can be expressed as Equation (29): where C i (i = 0, 1, 2, ..., n) is one kind of fault pattern; and n is the total number of fault patterns, which is equal to the dimension of the fault discrimination matrix.E (n+1)×(n+1) is an n + 1 dimensional identity matrix.

Simulation and Result Analysis of Exergy Loss of Typical Gas-Path Faults
In this paper, it is considered that the faulty condition of a gas path will affect some characteristics of components, such as compressor fouling, erosion, turbine wear, and damage, and the main performance parameters are its efficiency and flow rate.This section conducts fault implantation studies for a marine gas turbine low-pressure compressor, high-pressure compressor, high-pressure turbine, low-pressure turbine, and power turbine [21].To simulate the gas-path fault of the marine gas turbine with different severity, the severity of the fault is reflected in the decrease of the efficiency and flow of each component, as shown in Table 2.
In this paper, the exergy loss of the main components of the marine gas turbine is simulated, the performance parameter decline factor is determined, and the variation of exergy loss of each component with the severity of the fault is studied.
The efficiency and flow rate are implanted in the simulation model with severity from 20% to 100% (efficiency and flow rate reduce by 1-5%), respectively.The relative changes of the exergy loss of the main components under the condition of 0.5 to 0.9 are Taking the efficiency decline fault of the low-pressure compressor as an example, the exergy loss of each component of the marine gas turbine is analyzed.
The relative variation of exergy loss of marine gas turbine components with the severity of low-pressure compressor efficiency fault is shown in Figure 5a-d.As can be seen from Figure 5, as the efficiency of the low-pressure compressor is reduced, the relative exergy changes of each component of the marine gas turbine are increased.The relative change in the exergy loss of the low-pressure compressor is the highest in comparison with other components, which can indicate that the low-pressure compressor is the source of failure.It is found that the exergy loss of the component increases with the severity of the gas-path fault, and the change of the exergy loss of the component is clear when the efficiency fault occurs.
In addition, as can be seen from Figure 5a, when the low-pressure compressor malfunctions, the change degree of exergy loss is much greater than that of efficiency.Therefore, exergy loss is more sensitive to fault diagnosis compared with efficiency.

Fault-Detection Index
We select the error rate and the missing rate as the evaluation performance indexes of fault diagnosis.A reliable fault diagnosis system should minimize the error rate and the missing rate.If the total number of test event is T, then the missing rate is defined as the proportion of missed events in the specified events.Equation ( 30) is used to calculate the missing rate Mr.
where M is the number of missed events for detecting a specific fault; and L is the number of specific fault events.The error rate is defined as the proportion of misreported events.Equation ( 31) is used to calculate the error rate Fr.
where F is the number of misreported events for detecting a specific fault.

Example of Marine Gas Turbine Fault Diagnosis under Certain Operating Conditions
According to the established exergy model of the marine gas turbine, the simulation models of health conditions and faulty conditions was established in MATLAB/SIMULINK, and the simulation time set to 500 s, time step set to 0.2 s.Therefore, the total number of test event T is 2500.
The moment of gas path fault (C 1 , C 2 , ..., C 10 ) implantation for the marine gas turbine is shown in Figure 6 with 50 s as an interval.A fault is set at 20 s of each interval, and the system returned to normal at 50 s.
The proposed fault-diagnosis method was verified under the 0.8 perating condition which means the marine gas turbine outputs 80% of its rated power.The ambient temperature was set to a variable varying with the simulation time, as shown in Figure 7.The operating condition setting is shown in Figure 8.
The relative change in the exergy loss of each component of a marine gas turbine over the simulation time is shown in Figure 9.
Figure 9 shows the change of the relative exergy loss of each system when the gas-path fault of the marine gas turbine (C 1 , C 2 , ..., C 10 ) varies with the simulation time, which provides the database for subsequent fault diagnosis.The error rate is defined as the proportion of misreported events.Equation ( 31) is used to calculate the error rate Fr.
where F is the number of misreported events for detecting a specific fault.

Example of Marine Gas Turbine Fault Diagnosis under Certain Operating Conditions
According to the established exergy model of the marine gas turbine, the simulation models of health conditions and faulty conditions was established in MATLAB/SIMULINK, and the simulation time set to 500 s, time step set to 0.2 s.Therefore, the total number of test event T is 2500.
The moment of gas path fault (C1, C2, ..., C10) implantation for the marine gas turbine is shown in Figure 6 with 50 s as an interval.A fault is set at 20 s of each interval, and the system returned to normal at 50 s.
The proposed fault-diagnosis method was verified under the 0.8 perating condition which means the marine gas turbine outputs 80% of its rated power.The ambient temperature was set to a variable varying with the simulation time, as shown in Figure 7.The operating condition setting is shown in Figure 8.
The relative change in the exergy loss of each component of a marine gas turbine over the simulation time is shown in Figure 9.
Figure 9 shows the change of the relative exergy loss of each system when the gas-path fault of the marine gas turbine (C1, C2, ..., C10) varies with the simulation time, which provides the database for subsequent fault diagnosis.The error rate is defined as the proportion of misreported events.Equation ( 31) is used to calculate the error rate Fr.
where F is the number of misreported events for detecting a specific fault.

Example of Marine Gas Turbine Fault Diagnosis under Certain Operating Conditions
According to the established exergy model of the marine gas turbine, the simulation models of health conditions and faulty conditions was established in MATLAB/SIMULINK, and the simulation time set to 500 s, time step set to 0.2 s.Therefore, the total number of test event T is 2500.
The moment of gas path fault (C1, C2, ..., C10) implantation for the marine gas turbine is shown in Figure 6 with 50 s as an interval.A fault is set at 20 s of each interval, and the system returned to normal at 50 s.
The proposed fault-diagnosis method was verified under the 0.8 perating condition which means the marine gas turbine outputs 80% of its rated power.The ambient temperature was set to a variable varying with the simulation time, as shown in Figure 7.The operating condition setting is shown in Figure 8.
The relative change in the exergy loss of each component of a marine gas turbine over the simulation time is shown in Figure 9.
Figure 9 shows the change of the relative exergy loss of each system when the gas-path fault of the marine gas turbine (C1, C2, ..., C10) varies with the simulation time, which provides the database for subsequent fault diagnosis.The first step of diagnosis is to judge the overall state of the system; the diagnosis model is μ0.If μ0 = 1, the system is healthy; if μ0 = 0, the system is under fault.The second step is to judge the specific faults of each system.μ1-μ10 represent typical faults of marine gas turbines (C1, C2, ..., C10), respectively.μi = 1 means fault detected; and μi = 0 means fault undetected.(i = 1, 2, 3, ..., 10).
Using MATLAB Neural Network Toolbox to carry out the PNN fault diagnosis of the marine gas turbine, it can be seen in Figure 10 that the diagnosis effect of each fault is good when the fault occurs, and the error rate is relatively low, but the missing rate is relatively high.Omissions usually occur at the switching time between normal mode and fault mode and decrease with detection time.Error rate, missing rate, and detection time are shown in Table 3.As can be seen in Table 3, the probability neural network diagnosis method based on the exergy loss can diagnose the gas path  The first step of diagnosis is to judge the overall state of the system; the diagnosis model is μ0.If μ0 = 1, the system is healthy; if μ0 = 0, the system is under fault.The second step is to judge the specific faults of each system.μ1-μ10 represent typical faults of marine gas turbines (C1, C2, ..., C10), respectively.μi = 1 means fault detected; and μi = 0 means fault undetected.(i = 1, 2, 3, ..., 10).
Using MATLAB Neural Network Toolbox to carry out the PNN fault diagnosis of the marine gas turbine, it can be seen in Figure 10 that the diagnosis effect of each fault is good when the fault occurs, and the error rate is relatively low, but the missing rate is relatively high.Omissions usually occur at the switching time between normal mode and fault mode and decrease with detection time.Error rate, missing rate, and detection time are shown in Table 3.As can be seen in Table 3, the probability neural network diagnosis method based on the exergy loss can diagnose the gas path The first step of diagnosis is to judge the overall state of the system; the diagnosis model is µ 0 .If µ 0 = 1, the system is healthy; if µ 0 = 0, the system is under fault.The second step is to judge the specific faults of each system.µ 1-µ 10 represent typical faults of marine gas turbines (C 1 , C 2 , ..., C 10 ), respectively.µ i = 1 means fault detected; and µ i = 0 means fault undetected.(i = 1, 2, 3, ..., 10).
Using MATLAB Neural Network Toolbox to carry out the PNN fault diagnosis of the marine gas turbine, it can be seen in Figure 10 that the diagnosis effect of each fault is good when the fault occurs, and the error rate is relatively low, but the missing rate is relatively high.Omissions usually occur at the switching time between normal mode and fault mode and decrease with detection time.Error rate, missing rate, and detection time are shown in Table 3.As can be seen in Table 3, the probability neural network diagnosis method based on the exergy loss can diagnose the gas path failure of the marine gas turbine well and can diagnose the specific fault at a maximum of 3.3 s, which means its sensitivity is good.
Energies 2019, 12, x FOR PEER REVIEW 14 of 17 failure of the marine gas turbine well and can diagnose the specific fault at a maximum of 3.3 s, which means its sensitivity is good.
Figure 10.Fault-diagnosis result of a marine gas turbine system.

Conclusions
In this paper, a gas-path fault-diagnosis method of a marine gas turbine based on exergy loss and a probabilistic neural network was proposed.The exergy model of the three-shaft marine gas turbine was established based the exergy flow analysis to extract the exergy loss feature of the component under the healthy condition and typical faulty condition.A PNN classifier was introduced in the exergy analysis method to achieve good performance in detecting and isolating the gas-path fault.A simulation case study was conducted based on a three-shaft marine gas turbine with typical gas-path faults and environmental change.The simulation results show that the proposed diagnosis method can accurately detect the gas-path fault and locate the malfunction component.

Conclusions
In this paper, a gas-path fault-diagnosis method of a marine gas turbine based on exergy loss and a probabilistic neural network was proposed.The exergy model of the three-shaft marine gas turbine was established based the exergy flow analysis to extract the exergy loss feature of the component under the healthy condition and typical faulty condition.A PNN classifier was introduced in the exergy analysis method to achieve good performance in detecting and isolating the gas-path fault.A simulation case study was conducted based on a three-shaft marine gas turbine with typical gas-path faults and environmental change.The simulation results show that the proposed diagnosis method can accurately detect the gas-path fault and locate the malfunction component.The exergy losses of 10 typical gas-path faults with different severity were analyzed, and we found that the exergy loss relative change of components is the monotone increasing function of the operating condition and fault severity, respectively.Moreover, they are also highly sensitive to the performance parameters of components, which is conducive to locating the fault cause.
In this paper, the proposal was verified by a single gas-path fault.In future work, we will continue to improve and evaluate the performance of the proposed method on a multiple gas-path fault condition.

Figure 2 .
Figure 2. Exergy flow diagram of marine three-shaft gas turbine.

Figure 2 .
Figure 2. Exergy flow diagram of marine three-shaft gas turbine.

Figure 2 .
Figure 2. Exergy flow diagram of marine three-shaft gas turbine.


; fA(X) and fB(X) are respectively the probability density functions of the category A and B.

Figure 5 .
Figure 5.The relative variation of components exergy loss under a given faulty condition.

Energies 2019 , 17 Figure 5 .
Figure 5.The relative variation of components exergy loss under a given faulty condition.

Figure 7 .
Figure 7. Setting of ambient temperature for the marine gas turbine.

Figure 7 .
Figure 7. Setting of ambient temperature for the marine gas turbine.

Figure 7 . 17 Figure 8 .
Figure 7. Setting of ambient temperature for the marine gas turbine.

Figure 9 .
Figure 9. Relative change of exergy loss for marine gas turbine systems under simulated faults.

Figure 9 .
Figure 9. Relative change of exergy loss for marine gas turbine systems under simulated faults.

Figure 9 .
Figure 9. Relative change of exergy loss for marine gas turbine systems under simulated faults.

10 Figure 10 .
Figure 10.Fault-diagnosis result of a marine gas turbine system.

.
Exergy flow and exergy loss of marine gas turbine components.

Table 1 .
Exergy flow and exergy loss of marine gas turbine components.

Table 1 .
Exergy flow and exergy loss of marine gas turbine components.

Table 2 .
Summary of marine gas turbine gas-path fault simulation.