METHOD FOR DETECTING AN ANOMALY OF A ROLLING EQUIPMENT EXPLOITING A DEFORMATION SIGNAL FROM A RAIL SUPPORT

Abstract

The invention relates to a computer-implemented method for detecting an anomaly of a rolling equipment rolling on rails of a railway resting on a rail support. This method comprises a decomposition (DECOMP) by discrete wavelet transform of a measurement signal (S) transmitted by a strain sensor detecting the deformation of the rail support into an approximation signal (A.sub.J) and a residual signal (R.sub.J) and a search (RECH-PA) for outliers (PA) in the residual signal (R.sub.J) in order to detect an anomaly of the rolling equipment.

Claims

1. A computer-implemented method for detecting an anomaly of a rolling equipment rolling on railway rails resting on a rail support, comprising the steps of: applying a wavelet transform to a measurement signal transmitted by a strain sensor detecting a deformation of the rail support thereby decomposing said measurement signal into an approximation signal and a series of detail signals, summing all or part of the detail signals to form a residual signal; searching for outliers in the residual signal in order to detect an anomaly of the rolling equipment.

2. The computer-implemented method according to claim 1, wherein searching for outliers in the residual signal consists of searching for points of the residual signal which have an absolute value of the amplitude |r.sub.i| satisfying |r.sub.i|>.sub.,R+.sub.,R, where .sub.,R is the average noise contained in the residual signal, .sub.,R is the standard deviation of the noise contained in the residual signal and is a parameter for adjusting a detection sensitivity.

3. The computer-implemented method according to claim 1, further comprising a prior step of determining a level of decomposition of the wavelet transform, said level of decomposition minimising a square error given by w(.sub.,R.sub.,S).sup.2+(.sub.R.sub.,S).sup.2, where w is a weighting parameter, .sub.,R is the standard deviation of the noise contained in the residual signal, .sub.,S is the standard deviation of the noise contained in the measurement signal and .sub.R is the standard deviation of the residual signal.

4. The computer-implemented method according to claim 1, further comprising in the event that an anomaly of the rolling equipment is detected, classifying the detected anomaly as an anomaly of a first type or of a second type.

5. The computer-implemented method according to claim 4, wherein the detected anomaly is classified as an anomaly of the first type when it is associated with one single peak of the residual signal and is classified as an anomaly of the second type when it is associated with at least two single peaks of the residual signal of opposite signs.

6. The computer-implemented method according to claim 5, wherein the detected anomaly is classified as an anomaly of the second type when it is associated with outliers, one whereof has an amplitude that is less than a first negative threshold and another whereof has an amplitude that is greater than a second positive threshold.

7. The computer-implemented method according to claim 1, further comprising a step of determining a severity of a detected anomaly.

8. The computer-implemented method according to claim 1, further comprising a step of detecting peaks in the approximation signal.

9. A data processing system configured to implement the method according to claim 1.

10. A computer program product comprising instructions which, when the program is executed by a computer, cause same to implement the method according to claim 1.

Description

BRIEF DESCRIPTION OF THE FIGURES

[0015] Other aspects, purposes, advantages and characteristics of the invention will be better understood upon reading the following detailed description given of the non-limiting preferred embodiments of the invention, provided for illustration purposes, with reference to the accompanying figures, in which:

[0016] FIG. 1 shows one example of a measurement signal transmitted by a strain sensor detecting deformation of a rail support, this signal carrying a transient induced by an anomaly of a rolling equipment;

[0017] FIG. 2 is a close-up view of the signal in FIG. 1 focused on the transient that is characteristic of the anomaly;

[0018] FIG. 3 shows the measurement signal and the approximation signal at the top and the residual signal at the bottom, in which outliers were detected;

[0019] FIG. 4 is a close-up view of FIG. 3 focused on the transient that is characteristic of the anomaly;

[0020] FIG. 5 shows a residual signal carrying one single peak that is characteristic of an anomaly of the axle overload type;

[0021] FIG. 6 shows a residual signal carrying two single peaks of opposite signs that are characteristic of an anomaly of the wheel-flat type;

[0022] FIG. 7 shows one possible embodiment of the detection of an anomaly of the wheel-flat type;

[0023] FIG. 8 shows a measurement signal and the approximation signal thereof from which axle peaks are detected;

[0024] FIG. 9 is a flow chart of a method according to the invention.

DETAILED DESCRIPTION OF SPECIFIC EMBODIMENTS

[0025] The invention relates to a computer-implemented method for detecting an anomaly of a rolling equipment rolling on rails of a railway resting on a rail support. The invention exploits a measurement signal transmitted by a sensor capable of measuring the deformation of a rail support of a railway. The sensor can be joined to the surface of the support or integrated into the support.

[0026] The description below considers the example of a strain sensor of the optical fibre Bragg grating type, the measurement signal whereof is, for example, sampled at 500 or 1,000 Hz. FIG. 1 shows a recording of some thirty seconds of the measurement signal transmitted by such a sensor during which a train passes over the support. This recording carries a transient T induced by an anomaly of the train. FIG. 2 is a close-up view of the signal in FIG. 1 focused on this transient T. In particular, the characteristic M shape of the peaks linked to the passage of the axles of the locomotives and wagons over the support equipped with the sensor can be recognised. This recording can undergo preliminary processing including drift compensation and potentially normalisation. In order to detect a transient characteristic of a rolling equipment anomaly, the method according to the invention includes, with reference to FIG. 9, a step DECOMP consisting, where appropriate after the aforementioned preliminary processing, of decomposing, by discrete wavelet transform, the measurement signal S transmitted by the strain sensor into an approximation signal A.sub.J and a residual signal R.sub.J. This decomposition is carried out up to a level of decomposition J according to the following conventional procedure. The measurement signal S is decomposed into an approximation signal A.sub.1 which provides a smooth view of the original signal and into a detail signal D.sub.1 amplifying the high-frequency components of the signal. The approximation signal A.sub.1 is in turn decomposed into an approximation signal A.sub.1 and into a detail signal D.sub.2 (which thus appears as a corrective term between the two successive approximations A.sub.1 and A.sub.2). This procedure is repeated until the desired level of decomposition J is obtained. Ultimately, an approximation signal A.sub.J and a series of detail signals D.sub.1, D.sub.2, . . . , D.sub.J are obtained and the signal S is thus decomposed according to S=A.sub.1+.sub.jJD.sub.j. A residual signal R.sub.j is formed with the sum of all or part of the detail signals. Thus, the residual signal can consist of the complete sum from 1 to J of the detail signals, but can also constitute a choice of specific detail signals (for example for J=3, R.sub.3=D1+D3, D2 not included).

[0027] The decomposition of the signal S more specifically produces a set of coefficients denoted cA.sub.J and cD.sub.j (1jJ) respectively for approximation and details. Using these coefficients, the signal S can be reconstructed at the desired level on the one hand in order to obtain the approximation A.sub.J (low-frequency content) and the sum of the detail signals .sub.iJD.sub.j (high-frequency content) such that the signal S is decomposed according to S=A.sub.J+.sub.jJD.sub.j, with the relation A.sub.J-1=A.sub.J+D.sub.J between two successive approximations and

D.sub.j=cD(j,k).sub.j,k, where the variable k represents the temporal phase shift and .sub.j,k is the wavelet at the level j by phase-shifted k samples.

[0028] Unlike the sine and cosine functions used in the Fourier transform, the functions .sub.j,k are localised in time: only part of the samples is non-zero. By appropriately selecting the level of decomposition, the micro-deformation signal (i.e. the axle peaks) can be separated from the measurement noise and anomalies of the rolling equipment which are characterised by sudden transient phenomena localised in time.

[0029] The selection of the wavelet is guided by the form of the signal S to be decomposed, typically by selecting a wavelet resembling this signal. In the examples described below, the Symlet 5 wavelet is chosen.

[0030] The choice J of the level of decomposition is made so as to ensure the best possible separation compromise. In one possible embodiment, the method according to the invention comprises a prior step of determining a level of decomposition of the discrete wavelet transform decomposition of the measurement signal, said level of decomposition minimising the square error w(.sub.,R.sub.,S).sup.2+(.sub.R.sub.,S).sup.2, where w is a weighting parameter, .sub.,R is the standard deviation of the noise contained in the residual signal, .sub.,S is the standard deviation of the noise contained in the measurement signal S and .sub.R is the standard deviation of the residual signal.

[0031] More specifically, a good separation requires:

[0032] on the one hand, that the standard deviation of the noise contained in the residual signal .sub.,R and the standard deviation of the noise in the measurement signal .sub.,S are as close as possible, i.e. .sub.,R.sub.,S; and

[0033] on the other hand, that the standard deviation of the residual signal .sub.R remains close to the standard deviation of the noise in the measurement signal .sub.,S, generally slightly greater since a contribution from the axle peaks can subsist and because of the presence of an anomaly, i.e. .sub.R.sub.,S.

[0034] In doing so, the detail signals are guaranteed to primarily only contain the noise and the transients resulting from the presence of an anomaly.

[0035] The invention is not exclusive to this example of the chosen level of decomposition, and can be carried out according to other methods such as methods based on the energy contained in the detail signals for example.

[0036] The standard deviation of the noise contained in the measurement signal .sub.,S can be easily estimated over the part of the signal recorded before the passage of the train, for example over the first n seconds in the example in FIG. 1 (where n=3 seconds for example).

[0037] In one example embodiment exploiting the Symlet 5 wavelet, the selected level of decomposition is J=3. The top portion of FIG. 3 shows the measurement signal S and the approximation signal A.sub.3 and the bottom portion of the same FIG. 3 shows the residual signal R.sub.3. FIG. 4 shows a close-up view of FIG. 3 focused on the transient T that is characteristic of the anomaly. The standard deviation of the noise affecting the measurement signal .sub.,S is estimated to equal 1.08592 over the first 3 seconds of the signal. The level of decomposition J=3 gives a standard deviation of the noise contained in the residual signal .sub.,R equal to 1.03472 and a standard deviation of the residues .sub.R equal to 1.20381.

[0038] In one possible embodiment, a thresholding of the coefficients of the decomposition (for example a so-called soft coefficient thresholding) can be implemented in order to reduce the noise level in the reconstruction and obtain a residual signal that ideally only contains the transients characteristic of potential anomalies.

[0039] After the decomposition DECOMP of the measurement signal into an approximation signal and a residual signal, the method according to the invention comprises a step RECH-PA of searching for outliers in the residual signal in order to detect an anomaly of the rolling equipment. These outliers PA appear in the form of full circles in FIGS. 3 and 4. They correspond to points that get out from noise. The detection of an anomaly occurs for a succession of outliers, for example in a time window following an outlier, the size whereof is variable as a function of the speed of the rolling equipment.

[0040] The search for outliers in the residual signal can in particular consist of searching for the points of the residual signal, the absolute value of the amplitude |r.sub.i| thereof satisfies |r.sub.i|>+.sub.,R+.sub.,R, where .sub.,R is the average noise contained in the residual signal during the first n seconds (i.e. before the passage of the train), .sub.,R is the standard deviation of the noise contained in the residual signal and a is a parameter for adjusting a detection sensitivity. The average .sub.,R is generally substantially equal to zero as a result of the shift compensation carried out during preliminary processing. For example, =8 is chosen.

[0041] In one possible embodiment, an anomaly detected is classified CLAS, for example as a first anomaly type or as a second anomaly type, by means of an analysis of the residual signal. The first anomaly type is, for example, an axle overload which, as shown in FIG. 5, causes a transient Tps of the residual signal in the form of one single peak. The second anomaly type is, for example, a wheel-flat which, as shown in FIG. 6, causes a transient Tpd of the residual signal of the type of a sudden alternation about the M-shaped axle peak curve in the form of at least two single peaks of opposite signs. Thus, the anomaly detected can be classified as an anomaly of the first type when it is associated with one single peak of the residual signal and can be classified as an anomaly of the second type when it is associated with at least two single peaks of opposite signs.

[0042] The analysis of the residual signal used to carry out this classification can exploit the outliers previously detected in order to differentiate between the different types of anomaly. For example, the anomaly detected is classified as an anomaly of the second type when it is associated with outliers, one whereof has an amplitude that is less than a first negative threshold and another whereof has an amplitude that is greater than a second positive threshold. As shown in FIG. 7, symmetric thresholds s.sub.d and s.sub.d can be selected. The anomaly is thus classified as characteristic of a wheel-flat if the minimum of the outliers of the anomaly is less than s.sub.d and if the maximum of the outliers of the anomaly is greater than sd. Conversely, the anomaly is classified as being characteristic of an axle overload.

[0043] The method can further comprise a dating of an anomaly detected, for example as a function of the time of the maximum, in absolute value form, of the outliers of the anomaly, as a function of the time of the median, in absolute value form, of the outliers of the anomaly, or even as a function of the time of the first outlier of the anomaly.

[0044] The method can further comprise the determination of a severity of an anomaly detected. For example, for an anomaly of the axle overload type, this severity can correspond to the maximum of the outliers of the anomaly. For an anomaly of the wheel-flat type, this severity can, for example, correspond to the largest deviation in amplitude y (see FIG. 7) between outliers of the anomaly.

[0045] It should be noted that by having an annotated database of cases of anomalies, a supervised classification algorithm can be trained and used to effectively recognise different types of anomalies.

[0046] The detection of the transit of the axles on the support and the strain sensor thereof is generally based on peak detection algorithms which search for fast variations in the deformation measurement signal S. The adjustment parameters can be chosen so as to make these algorithms less sensitive to noise, for example by setting a minimum distance between two axle peaks or by setting a minimum variation in deformation. However, this detection is not robust against anomalies contained in the signal since these are transient phenomena that also have fast variation.

[0047] Within the scope of the invention, the reconstruction of the approximation signal A.sub.J provides a signal from which measurement noise and the anomalies detected have been removed, on which the detection of axle peaks can be carried out. Thus, in one possible embodiment of the invention, the method further comprises a step RECH-Ep of detecting peaks in the approximation signal. In this respect, FIG. 8 shows the result of the detection of axle peaks EP in the area of the anomaly used in the preceding example. This robust detection of the axle peaks in particular allows for automatic recognition of the rolling equipment. More specifically, the locomotives and wagons have known and documented features, in particular the car length and location of the trucks which allow the inter-axle distances to be calculated. If two strain sensors are used, separated by a known distance, the speed of the train can be easily determined by the time difference between the deformation measurements. Exploiting the approximation signal with a low noise or anomaly level prevents distortion in the assessment of the inter-axle distances and deteriorating the automatic recognition of the rolling equipment.

[0048] The invention is not limited to the method described hereinabove, but also extends to a data processing system configured to implement same, as well as to a computer program product comprising instructions which, when the program is executed by a computer, result in the former implementing this method.