Loboda Igor

Dynamic Model Identification Procedure.

Introduction

Modern gas turbine engines represent powerful and efficient sources of mechanical energy for numerous industry branches and transport objects. For a gas turbine that is considered as a complex and expensive system, condition monitoring technologies have the potential to save millions of dollars per year [1], through lowering fuel consumption, reducing catastrophic failures, decreasing mean time to repair, and optimizing maintenance planning. This will significantly reduce life cycle cost and improve competitive position of companies that maintain engines equipped with condition (health) monitoring systems.

Aircraft engine monitoring systems have become increasingly standard in the last two decades [2]. Above traditional aircraft application and usage in power plants of natural gas pumping and electrical power production, condition monitoring is used, for example, in shipboard [1] and battle tank [3] propulsion systems.

The thermodynamic (temperature, pressure, RPM, etc.) sensors are located at strategic points along the gas flow in the engine to provide more detailed thermodynamic picture of the engine’s state. The algorithms of failure detection that analyze registered thermodynamic parameters and use gas turbine models for diagnostic aims can be considered as principal algorithms of monitoring systems [2]. Gaspath failures can be detected by these algorithms, for example, compressor blade contamination, aerodynamic surfaces distortion, and seal wear as well as measurement system malfunctions.

Wide usage of mathematical models in diagnosing process is explained by high cost of physical failure modelling, an infrequent display of gas turbine failures, and other causes. Steady state regimes and corresponding static models are traditionally chosen for diagnostic needs [4]; however transient regime analysis [5] is also involved in the diagnosis now. For feather diagnosing enhancement, a dynamic model and its identification procedure had been developed later and adapted to a stationary gas turbine power plant [6]. As every new tool perspective for practical implementation in active monitoring system, this procedure must be carefully verified. Therefore, a statistical testing of the dynamic model identification (DMI) procedure was carried out. The testing is described and the results are discussed in this paper as well as perspectives of DMI-procedure incorporating into the monitoring system.

Dynamic Model

Dynamic nonlinear gaspath model describes behavior of thermodynamic parameters on transient regimes and can be presented by the common expression

,                     (1)

where the vector of regime and atmospheric conditions  is given as a function of time, and a separate influence of the time variable t is explained by an inertia nature of gas turbine dynamic processes. Every engine component (compressor, turbine, combustion chamber etc.) in this model is presented by its performance. The vector of state parameters  is used for describing and simulating the engine failures. These parameters are able to displace component performances in different directions and simulate various failures by this mode. The state parameters have a relative form; and the normal value is one.

The thermodynamic parameters of the model (1) are computed numerically as a solution of the system of differential equations in which the right parts are calculated from a system of algebraic equations reflecting the conditions of the components combined work on transient regimes. The measured values  differ from the model generated ones due to the model errors  and the measurement errors , therefore

.                     (2)

Similarly, the expression for regime and atmospheric conditions

                           (3)

is formed.

It is suggested to divide the total measurement error level into three components influencing on identification process in different ways: ε0Y, ε0U - levels of systematic errors; ε1Y, ε1U - levels of long-term random fluctuations that are changed from one transient regime to another; ε2Y, ε2U - levels of short-term random fluctuations that are changed during a transient regime.

Model Identification

The objective of model identification consists in finding such values of model internal parameters which minimize a discrepancy between model external parameters and measured ones. The state parameters  are chosen as such internal parameters that must be estimated. This is explained by the known fact that engine components theoretical performances used into the model are not sufficiently certain and may be specified on real data . So, the state parameter estimations  may be expressed as follows

.       (4)

Besides the better model accuracy resulting from such an adjustment, the simplification of the diagnosing process is reached because the found state parameters contain information of current technical condition of each component.

For any iteration of number n+1 of the identification procedure, the current estimation may be written as follows

,                       (5)

in which the current correction  presents a regularized solution of the linear system

,                   (6)

where C - matrix of influence coefficients of state parameters on thermodynamic parameters calculated in registration points 1, 2,…, NT of the transient regime chosen for the identification;  - vector of discrepancies between the model values  and measured ones  formed in NT corresponding time-points.

To solve the linear system (6), a standard technique is used that selects an optimal value of the regularization coefficient α. The variation boundary αmax may be changed. The iterations are repeated until a moment when current increments of state parameters and thermodynamic parameters will be sufficiently low or on reaching the established cycle number NI.

The software of the developed DMI-procedure includes about 70 program modules; around 90 percents of modules were imported from static model identification procedure and are time-proved. The software was tested on simulated data and real information and had demonstrated a correct functioning and quick convergence (Loboda, 2002).

Identification Procedure Statistical Testing

Although the DMI-procedure is critical to computer operating speed, a stochastic measurement simulation and numerous repetitions of the identification in the course of procedure statistical testing are still possible. In this paper, an influence of various factors on identification accuracy is investigated by such a testing.

The statistical testing was carried out in the following sequence: 1) the failures are simulated by state parameters  that displace component performances; 2) corresponding thermodynamic parameters are generated by the model with changed component performances; 3) a random measurement noise is added to these thermodynamic parameters and the inlet and control parameters; 4) the DMI-procedure is executed on the simulated data  and ; 5) error of the estimations  is determined.

In an external testing cycle this simulation and identification sequence is repeated NS times for reliable determination of obtained average accuracy. Behavior and accuracy of tested procedure are checked with the following cycle - averaged values: δΘj,  - discrepancies between simulated and estimated values of every state parameter ΘJ and on average.

The state parameters present a primary interest for detection algorithms therefore these criteria are chosen to define final identification accuracy.

The factors affecting the behavior and accuracy and analyzed in the paper may be listed as follows:

-  cycle number NS;

-  regularization coefficient boundary αmax;

-  iteration number NI;

-   measured thermodynamic parameters structure and number m;

-  state parameters structure and number r;

-   measurement error levels ε0Y, ε1Y, ε2Y, ε0U, ε1U, ε2U;

-   dynamic process profile  given in NT points;

-  simulated state parameters .

A base calculation was established; factors were varied later relatively the base calculation conditions and analyzed independently; most interesting and important results were verified by additional calculations in different conditions.

The base calculation of identification procedure testing was executed in the following conditions:

NS = 1000; αmax = 5600; NI = 3;

m = 8; r = 6;

ε1Y = ε1U = 0.008;

vector  is given in 5 time-points and includes a linear change of high pressure rotor speed;

Θmod1 = - 0.03 (compressor flow parameter) and Θmod4= - 0.02 (high pressure turbine efficiency parameter).

The first task of factors analysis was to establish proper cycle number NS to provide necessary precision level of the testing calculations and, with the number 1000, the level [±(2-3)% for δΘj and ±1% for ] was reached. This level is sufficient for reliable determination of the factor influence on the criteria and presents a compromise between the achieved precision and required computer time.

To determine the identification procedure parameters, iteration increments and current accuracies were plotted versus iteration number variable. It was noted that more and more short steps follow after the largest first one, and the increments are stabilized to the iterations 7-8, but accuracy stabilization is reached already on the iterations 2-3. So, the value NI = 3 is accepted as basic.

From the investigations related with gas turbine model identification, it is known that measured parameters number increase and/or estimated parameters number decrease lead to better estimation accuracy.  To verify and specify the noted influence five calculations were fulfilled in addition to base calculation.

Consequent comparison of these calculations with the basic one shows that the mentioned rule is kept generally. It may be noted too that the parameter structure also influences on the accuracy.

During one transient regime, the long-term noise induces a constant displacement of all registered values of a thermodynamic parameter and may be considered as a systematic error which is transferred wholly on state parameter estimations. In contrast, the short short-term noise produces the random displacements and the averaged error only is transmitted on the estimations. So, it may be supposed beforehand that the short-term noise will induce lower identification error than the long-term one. To verify this supposition, calculations with short-term noise simulation only were executed and the results are presented in Table 1 as well as the base calculation result. From comparison of the calculations D1 and A follows that the short-term noise causes more than two times lower estimation error. Comparison of the calculations D1, D2, and D3 confirms noted influence of the measured parameters number m.

Table 1. Short-term noise calculations

Calculation designations	Calculation conditions
A	basic (ε1Y=ε1U=0.008)	0.00557
D1	ε2Y=ε2U=0.008	0.00256
D2	ε2Y=ε2U=0.008, m = 6	0.00835
D3	ε2Y=ε2U=0.008, m=10	0.00254

In constant identification conditions, systematic measurement errors will induce identification systematic errors. On the other hand, the condition change may lead to a shift of the estimation and estimation errors. To evaluate this effect, three sets of fixed systematic measurement errors were formed and DMI-procedure was repeated for every set in the seven different condition variants. Related with variant change scatters of six estimated parameters are presented in Table 2. It is seen that the scatters may reach the level of simulated state parameters. So, systematic measurement errors may induce the random estimation errors when random changes of identification conditions take place. Systematic estimation errors are not considered here because it is supposed that in the diagnosis, the time-series of estimation will be analyzed and a relative estimative change only will be considerable.

Table 2. Systematic measurement error influence

Sets                     Estimation scatters

1        0.0363  0.0139  0.0170  0.0076  0.0262  0.0146

2        0.0099  0.0010  0.0059  0.0035  0.0074  0.0032

3        0.0122  0.0064  0.0121  0.0047  0.0143  0.0072

To estimate dynamic process and failure type on the identification accuracy, the process profile (regime parameter change during the process and process total time) and failure development were modified. The results demonstrate the estimation error stability to these factors (the scatter does not exceed the testing error level).

So, an attempt had been made to analyze systematically a maximal number of factors affecting the identification accuracy and to determine numerically this influence. Applied statistical testing included the model-based simulation of failures and generation of random measurement errors with theoretical distribution.

Conclusions

Thus, to verify the developed dynamic model identification procedure for including into gas turbine health monitoring systems, the statistical testing had been carried out. A lot of factors affecting the identification accuracy were analyzed. Sufficient factors had been determined, for instance, the systematic and random long-term measurement errors. An invariance of the accuracy to other group of factors such as dynamic process profile and simulated failure type was also noted.

Acknowledgments

The work has been carried out with the support of the National Polytechnic Institute of Mexico (projects 20021022 and 20030903).