The present invention relates to the field of monitoring a degradation of an on-board device of an aircraft, such as a system for lubricating a turbojet engine, a feeder for supplying fuel to a turbojet engine, a redundant acquisition chain of a turbojet engine, a particle filter, etc.
In order to monitor an on-board device, it is known to form an indicator which is characteristic of a degradation of the on-board device. This indicator is known to a person skilled in the art as an abnormality score. Conventionally, an abnormality score is formed from measurements of physical parameters of the on-board device such as a geometric position, a control current, an opening angle, a temperature, etc. The abnormality score is characteristic of the degree of damage of the degradation. Preferably, an abnormality score is formed during each flight of the aircraft.
By way of example, the patent application filed in France on 5 Dec. 2011 under the number 1161175 by the company SNECMA discloses a method for monitoring a fuel metering device of a turbojet engine, wherein the degradation relative to the drift of the standing current of a servo-valve controlling the fuel metering device is monitored. In order to form the abnormality score which is characteristic of the drift of the standing current of the servo-valve, a plurality of measurements of the average control current are carried out, and the norm is taken.
In order to determine whether the on-board device is actually degraded, the monitoring method includes a step of comparing an abnormality score obtained for a given flight of the aircraft to a decision threshold and a step of emitting an alarm when the decision threshold is exceeded. Therefore, by following the evolution of the abnormality score, it can be detected whether the degree of degradation increases, the risk of malfunction of the on-board device can be anticipated, and the management of the maintenance operations can be improved.
Setting the decision threshold for a given degradation is crucial given that, firstly, an overly low decision threshold leads to the emission of frequent alarms when the degradation is not genuine (false alarm) and, secondly, an overly high decision threshold inhibits the emission of alarms when the degradation is genuine (non-detection).
Conventionally, for each degradation that is to be detected, the value of the decision threshold is defined empirically. In order to ensure maximum safety, the value of the decision thresholds is generally undervalued in order to minimise the risk of non-detection. As a result, the number of false alarms remains high, thereby presenting a disadvantage for airlines, who are obligated to frequently carry out maintenance operations when the on-board aeronautical device is not degraded. In order to eliminate this disadvantage, airlines impose specifications on the manufacturers of on-board devices imposing limits on the risk of errors. Taking into account the precision required, any empirical definition of a decision threshold for a method for monitoring malfunctions is thus ruled out.
There is therefore a need to determine in a reliable and precise manner the value of a decision threshold for monitoring a degradation of an on-board device of an aircraft.
The invention relates to a method for monitoring a degradation of an on-board device of an aircraft, implemented by a computer, the degree of the degradation of the on-board device being defined by an abnormality score formed from measurements of physical parameters of the on-board device, the monitoring method including a step of comparing an abnormality score obtained for a given flight of the aircraft to a decision threshold and a step of emitting an alarm when the decision threshold is exceeded, the decision threshold being determined automatically for a given alarm probability Pa, corresponding to the probability of an alarm being emitted during the monitoring method when the on-board device is undamaged, by means of the following steps:
By means of the invention, the decision threshold is automatically determined with a great deal of precision according to the alarm probability Pa corresponding to the requirements of the airlines. Such a method is reliable by comparison with the methods according to the prior art, thereby making it possible to detect in a precise manner any degradation of an on-board device and to anticipate any malfunction of said device during monitoring.
Preferably, the non-parametric estimator of the probability density is a Parzen window estimator. Such an estimator is very advantageous, as it gives the same importance to the central portion of the distribution as to the end portions which are the most critical for monitoring degradation in the field of aeronautics, for which the requirements in terms of reliability are significant.
Preferably, a Parzen window estimator is defined by its window width h. Preferably, the distribution being formed from a number N of abnormality scores for a plurality of flights of the aircraft without degradation, the distribution having an estimated standard deviation {circumflex over (σ)}, the Parzen window has a width h defined according to the following formula:
Such a Parzen window width h makes it possible to guarantee precision of the adjustment for monitoring on-board devices, thereby ensuring precision of the determined decision threshold.
According to a preferred aspect of the invention, the alarm probability Pa is obtained by the following formula:
in which formula:
Pe corresponds to the probability of the on-board device being undamaged when an alarm is emitted;
Pd corresponds to the probability of a degradation for a given flight; and
(1−β)a priori corresponds to the probability of detecting a degradation a priori.
Advantageously, the alarm probability Pa is determined on the basis of an error probability Pe which corresponds to the practical requirements of airlines who wish to reduce the number of unnecessary maintenance steps for the on-board devices.
More preferably, the method includes:
These steps advantageously make it possible to refine the value of the decision threshold quickly and reliably when the alarm probability Pa is determined on the basis of an error probability Pe.
The invention will be better understood upon reading the following description, given purely by way of example and with reference to the accompanying drawings, in which:
It should be noted that the drawings show the invention in a detailed manner in order to implement the invention, said drawings of course being able to be used to better define the invention where necessary.
The automatic determination of a decision threshold will be set out for a method for monitoring a degradation of an on-board device of an aircraft such as is known from the prior art, for example, from the patent application under the number 1161175 by the company SNECMA.
In a known manner, the degree of degradation of the on-board device is defined by an abnormality score obtained by calculating using measurements of physical parameters of the on-board device as set out above. Preferably, the abnormality score is formed by an on-board computer of the aircraft which is connected to the on-board device by a communication link. In a known manner, as set out above, the monitoring method includes a step of comparing an abnormality score obtained for a given flight of the aircraft to a decision threshold and a step of emitting an alarm when the decision threshold is exceeded.
When an alarm is emitted, a step of maintenance of the on-board device must be carried out by the airline in order to avoid a malfunction of said on-board device.
Alarm Probability
Automatically setting a decision threshold is defined for a given alarm probability Pa. Alarm probability is understood to mean the probability of an alarm being emitted during the monitoring method when the on-board device is undamaged, that is to say has no genuine degradation. In other words, an alarm probability Pa of 1% corresponds to one alarm being emitted per 100 flights without degradation.
In practice, the alarm probability Pa is not directly known. Indeed, the specifications of the airlines impose as a requirement an error probability Pe which corresponds to the probability of the on-board device being undamaged when an alarm is emitted. In other words, an error probability Pe of 1% corresponds to one undamaged on-board device per 100 alarms emitted.
In practice, when an alarm is emitted by the monitoring method, a maintenance operation is carried out by the airline. Limiting the error probability Pe thus allows airlines to limit the number of maintenance operations carried out on an on-board device which is undamaged, such operations leading to unnecessary operating costs.
Preferably, the alarm probability Pa is determined on the basis of the error probability Pe by means of the conditional probability formula (F1) below.
The formula (F1) includes the following parameters:
The formula (F1) is derived from the equations of conditional probability in which the probability Pa corresponds to an alarm when the device is undamaged (Pa=P(Alarm|Undamaged)) and the probability Pe corresponds to a undamaged device when an alarm is emitted (Pe=P(Undamaged|Alarm)).
In other words
In addition, as
P(Alarm)=P(Alarm∩Undamaged)+P(Alarm∩Degraded)
P(Alarm)=P(Alarm|Undamaged)·P(Undamaged)+P(Alarm|Degraded)·P(Degraded)
Where
P(Degraded)=Pd
P(Undamaged)=1−Pd
P (Alarm|Degraded)=1−β
P(Alarm|Undamaged)=Pa
P(Undamaged|Alarm)=Pe
P(Alarm)=P(Alarm|Undamaged)·P(Undamaged)+P(Alarm|Degraded)·P(Degraded)
P(Alarm)=Pa·(1−Pd)+(1−β)·Pd
P(Alarm)=Pa·(1−Pd)+(1−β)·Pd
By means of the two equalities relating to P (Alarm), (F1) is thus derived.
In the formula (F1), the probability of a degradation for a given flight Pd is known by experimentation or estimation and can, for example, be approximately 10−7. The probability of detecting a degradation (1−β)a priori is fixed a priori at 1 and refined by iteration as will be explained in detail below in order to improve the precision of the decision threshold.
By way of example, conventionally an alarm probability Pa of approximately 10E-9 for a required error probability Pe of approximately 1%, a probability of detecting a degradation (1−β)a priori of approximately 1 and a probability of a degradation for a flight Pd of approximately 10−7 are obtained.
Abnormality Score, Distribution and Allocation Function
The on-board device on the aircraft conventionally comprises a plurality of sensors in order to measure physical parameters of the on-board device (control current, geometric positions, temperatures, etc.). Conventionally, some of the physical parameters measured are used to form an abnormality score which is based on the degree of degradation of the on-board device (positioning error, drift of the control current, overheating, etc.). An abnormality score is obtained by transformation of one or more measurements of physical parameters of the on-board device. For example, the measurements obtained are standardised by an absolute value or by the Mahalanobis distance.
As shown in
The distribution of the probability density of the abnormality score is shown in
The alarm probability Pa, which as a reminder corresponds to the probability of an alarm being emitted during the monitoring method when the on-board device is undamaged, is conventionally approximately 10E-9. In other words, the decision threshold S according to the alarm probability Pa is an abnormality score value above which are 10E-9 of the abnormality scores obtained for flights without degradation.
Even if a very large number of abnormality scores are available, the distribution in
Adjustment of the Distribution Function
In order to improve the interpretation of the distribution, with reference to
This non-parametric approach makes it possible to obtain an adjustment close to the distribution because a Parzen window adjustment advantageously gives the same weight to the ends of the distribution as to the central portion, thereby making it possible to interpret the end quantiles in a precise manner.
As a reminder, the Parzen window adjustment is a non-parametric adjustment for estimating the probability density of a random variable. It is commonly referred to as kernel density estimation because functions referred to as “kernels” are used to estimate the probability density of a random variable. The analytical expression of the Parzen probability density is:
Where K and h are, respectively, the kernel function and the width or depth of the Parzen window.
The concept of the Parzen window adjustment is that of estimating the probability density of a random variable for each value of x by means of a kernel function K( ) which most of the time is based on probability density. The closer the observation x is to the training samples xi, the more significant the contribution of x is in the kernel function {circumflex over (ƒ)}h(x) centred on xi. Inversely, the less close the observation x is to the training samples xi, the less significant the contribution of x is in the kernel function {circumflex over (ƒ)}h(x) centred on xi. Then, the estimator of the probability density function is formed by an average of the values of the kernel functions. This estimator is governed by a smoothing parameter “h”, which is referred to as a width or depth of the window. The estimation of probability density, which is dependent on the smoothing parameter “h”, has good statistical properties. Under certain non-restrictive conditions on “h”, the Parzen window estimator is consistent.
There are several kernel functions (Gaussian, box, triangle, etc.) but the performance of the Parzen window is mainly based on the selection of the width or depth of the window “h”. A compromise between the bias and the variance of the estimator has to be made in order to select the width “h”.
The width of the window can be selected by cross validation, by maximising the probability of the kernel function. In this example, a Gaussian kernel (F2) has been used. As for the width or depth of the window, it is based on an estimation of the standard deviation {circumflex over (σ)} and of the size N of all of the training samples according to formula (F3).
A Parzen adjustment is very close to the elementary distribution observed for several flights, as will be set out in the following, which makes an optimal adjustment by comparison with the other known adjustments.
Adjusted Allocation Function
In order to determine the value of the decision threshold S in a simple and quick manner, the adjusted distribution function Daj is conventionally integrated in order to obtain an adjusted allocation function Raj which is also continuous.
The adjusted reallocation function Raj of the adjusted distribution function Daj from
In addition, in order to respect the requirement of an alarm probability Pa of approximately 10E-9, it is necessary to precisely read the precedent of the adjusted allocation function Raj, for which said function is equal to 1-Pa, the precedent corresponding to the desired decision threshold S. Indeed, this signifies that 1-10E-9 of the abnormality scores for an undamaged device are less than the decision threshold S. In other words, 10E-9 of the abnormality scores for an undamaged device which can generate an alarm are greater than the decision threshold S. The precedent corresponds to the desired decision threshold S for the determined alarm probability Pa.
Determination of the Decision Threshold
With reference to
By way of information, a comparison of several types of adjustments of the allocation function (Gaussian adjustment, gamma adjustment and Parzen adjustment) has been shown in
This precision of the Parzen window adjustment translates into a precise determination of the decision threshold. Indeed, as shown in
Refining the Decision Threshold Value
A determination of the decision threshold has previously been set out in which the alarm probability Pa is known or estimated from the error probability Pe. When the alarm probability Pa is estimated, a step of refining the decision threshold S can optionally be carried out by refining the probability of detecting a degradation (1−β) as shown in the schematic diagrams in
As shown in
As shown in
In an implementation of the invention, the method for monitoring a degradation of an on-board device of an aircraft, the method for automatically determining a decision threshold and the method for generating degradations are implemented using a computer, preferably using a plurality of processors of the computer.
By way of example, the monitoring method is implemented in order to monitor a fuel metering device of an aircraft turbojet engine.
The processor for determining the decision threshold receives a plurality of abnormality scores without degradation from the monitoring processor and determines, for a given error probability Pe and a probability for detecting a priori a degradation using the monitoring method, the decision threshold S. Once the decision threshold S is determined, this is transmitted to the monitoring processor, which can then compare the calculated abnormality scores to said decision threshold S in order to monitor the evolution of the degradation of the fuel metering device during the flights of the aircraft.
In order to refine the value of the decision threshold S, the processor for generating degradations simulates degraded flight data which are submitted to the monitoring processor which emits a certain number of alarms on the basis of the data received, thereby making it possible to deduce therefrom a posteriori the probability of detecting a degradation using the monitoring method. This value is then communicated to the processor for determining the decision threshold S which supplies a new value of the decision threshold S for the new detection probability obtained.
The method is iterated until a convergence of the value of the decision threshold S is obtained. In practice, a satisfactory convergence is obtained from two iterations.
The invention has been set out for a fuel metering device but it can be applied to any on-board device of an aircraft.
Number | Date | Country | Kind |
---|---|---|---|
12 54506 | May 2012 | FR | national |
Filing Document | Filing Date | Country | Kind |
---|---|---|---|
PCT/FR2013/051059 | 5/15/2013 | WO | 00 |
Publishing Document | Publishing Date | Country | Kind |
---|---|---|---|
WO2013/171427 | 11/21/2013 | WO | A |
Number | Name | Date | Kind |
---|---|---|---|
5455777 | Fujiyama | Oct 1995 | A |
6049741 | Kawamura | Apr 2000 | A |
6409661 | Murphy | Jun 2002 | B1 |
6532421 | Miwa | Mar 2003 | B2 |
6599241 | Murphy | Jul 2003 | B1 |
6880967 | Isozumi | Apr 2005 | B2 |
6892163 | Herzog | May 2005 | B1 |
7243042 | Plotts | Jul 2007 | B2 |
7440811 | Giebels | Oct 2008 | B2 |
7933754 | Goebel | Apr 2011 | B2 |
8321157 | Omori | Nov 2012 | B2 |
9014918 | Hagen | Apr 2015 | B2 |
9251502 | Schoonveld | Feb 2016 | B2 |
20010034581 | Miwa | Oct 2001 | A1 |
20050143873 | Wilson | Jun 2005 | A1 |
20060142976 | Bonanni et al. | Jun 2006 | A1 |
20080255776 | Beard | Oct 2008 | A1 |
20080312783 | Mansouri | Dec 2008 | A1 |
20100030519 | Collier | Feb 2010 | A1 |
20100100259 | Geiter | Apr 2010 | A1 |
20110265026 | Tarassenko | Oct 2011 | A1 |
20130274991 | Cheriere | Oct 2013 | A1 |
20130274992 | Cheriere | Oct 2013 | A1 |
20130325219 | Pitard | Dec 2013 | A1 |
20140149325 | Clifton | May 2014 | A1 |
20140257716 | Lajnef | Sep 2014 | A1 |
20150227837 | Clifton | Aug 2015 | A1 |
Number | Date | Country |
---|---|---|
1 677 172 | Jul 2006 | EP |
Entry |
---|
International Search Report Issued Jul. 9, 2013 in PCT/FR13/051059 Filed May 15, 2013. |
Number | Date | Country | |
---|---|---|---|
20150134198 A1 | May 2015 | US |