JP7391742B2 - Abnormality analyzer - Google Patents

Description

The present invention includes a detection means for detecting an original time domain signal emitted by a periodically operating device;
The present invention relates to an abnormality analysis device comprising an abnormality analysis means for analyzing an abnormality of a device through learning based on the original time domain signal detected by the detection means.

This type of abnormality analysis device is disclosed in Patent Document 1 and Patent Document 2.
In the technology disclosed in Patent Document 1, in order to detect an abnormality in a rotating machine such as a motor such as a pump or a blower (corresponding to the device of the present invention), a spectrum is extracted from a sound/vibration signal (original time domain signal of the present invention). , skewness, kurtosis, crest factor, minimum value rate, local maximum rate, maximum value rate, stability index, frequency crest factor, transmission band, frequency height ratio, etc. are extracted, and the deviations from normal conditions are detected. A threshold value is set to determine an abnormality in the equipment.

In the technique disclosed in Patent Document 2, cluster analysis is performed based on frequency characteristics acquired according to the rotation speed, and abnormalities are monitored.

These techniques can basically be said to be anomaly analysis techniques for devices equipped with rotating equipment that use machine learning techniques.

Patent No. 4369321 Japanese Patent Application Publication No. 2018-138909

However, in a rotating device that operates with combustion, such as an engine, unlike a device that rotates steadily such as a motor, partial waveform disturbances may occur due to abnormal combustion. When targeting such signals, if anomaly analysis is performed by extracting features in the frequency domain using FFT analysis, time information will be lost in the feature extraction process, and partial or instantaneous events will be lost. changes may be overlooked and abnormality analysis may not be performed properly.

Let me briefly explain this issue here.
The results of determining the degree of abnormality using an autoencoder as an abnormality analyzer are shown in FIGS. 14 and 4.

In abnormality analysis when using an autoencoder, n waveform data (quantized area signals) observed when the device is in normal condition, and n waveform data observed at the point of time to be diagnosed. Diagnose whether or not there is an abnormality based on the difference between the two.

Specifically, an autoencoder consists of an input layer, an intermediate layer, and an output layer, and reproduces a signal input to the input layer at the output layer. The number of nodes in the input layer and the output layer are assumed to be the same. On the other hand, the number of nodes in the middle layer is smaller than that in the input layer. Dimensional compression is performed on the intermediate layer, and feature points that best represent the features of the input signal are extracted. The output layer performs dimension restoration based on the extracted feature points. The dimension compression and dimension restoration methods can be constructed by repeating learning using n waveform data when the device is operating in a normal state.

When waveform data in a normal state is input to the autoencoder that has been trained in this manner, the restoration is performed properly because the feature points and the restoration method match. On the other hand, if abnormal waveform data different from normal waveform data is input, the restoration cannot be performed successfully because the feature points that should originally be extracted and the restoration method are different.

The formula of the intermediate layer of the autoencoder can be, for example, Formula 1. In Equation 1, W represents a weighting coefficient and b represents a bias. f(x) is an activation function, and a Rele function can be used. f(x)=max(0,x), and when the input value is less than or equal to 0, it becomes 0, and when it is greater than 1, the input is output as is.

a _i = f(ΣW _ij x _j + b _i ) Σ is the sum of j=1 to N Formula 1

The autoencoder learns by varying the weighting coefficient W and bias b so that the difference between the input and output is minimized using training data with n waveform data as input and output. can be completed.

In this way, by using an autoencoder constructed using normal waveform data and monitoring the amount of deviation in restoration that involves dimensional compression and dimensional restoration, it is possible to diagnose an abnormality in the device. For example, an abnormality can be diagnosed from the amount of deviation in restoration by determining the degree of abnormality E in Equation 2 below. An abnormality is determined when the degree of abnormality E exceeds a predetermined threshold. When the difference between the input value and the output value is large, the abnormality degree E becomes a large value. In Equation 2, x _i represents an input value, and x _i ' represents an output value.

E=SQRT(Σ(x _i '−x _i ) ² ) Σ is the sum of i=1 to n Formula 2

In this example shown below, the analysis tool as an autoencoder was implemented using python. This autoencoder library uses the neural network library Keras, which is included in the machine learning library TensorFlow developed by Google.

The results shown in Figure 14 are obtained by frequency-converting the original time-domain signal to be analyzed, performing abnormality analysis based on the obtained frequency-domain waveform data, and determining the degree of abnormality (absolute average error). be. In the figure, the horizontal axis is time.

As can be seen from this figure, the degree of abnormality is distributed between 0.01 and 0.02 at almost all times, and it decreases from 0.10 to 0.10 after 15.00, which is written as "startup evaluation point". It has a high value. In this figure, "Learning validation" indicates the point in time when the autoencoder was constructed when the equipment was in a normal state. .015) was secured.

On the other hand, FIG. 4 shows the results of analysis using an autoencoder according to the method of the present invention shown later. The notation in the drawings is the same. As a result, the base level of the abnormality level increases, and in Figure 14, changes in the abnormality level are observed even at times when no increase in the abnormality level was observed (for example, from 8:00 on the first day to the second day). (before 11:00).
Given this situation, there is room for improvement in anomaly analysis in the frequency domain.

As described above, the main object of the present invention is to obtain an abnormality analysis device that can perform abnormality analysis with high sensitivity.

The first characteristic configuration of the present invention is:
detection means for detecting raw time domain signals emitted by periodically operating equipment;
an abnormality analysis means for analyzing an abnormality of the device based on the original time domain signal detected by the detection means,
a storage unit for storing a time domain signal at a specific time interval and a specific phase detected by the detection means and extracted from the original time domain signal emitted by the normal device as a normal time domain signal;
a reference time domain signal generating means for generating a reference time domain signal as a diagnostic standard based on the normal time domain signals sequentially stored in the storage unit;
the reference time domain signal is classified based on the frequency at which the equipment operates;
Regarding the original time domain signal of the diagnosis target detected by the detection means,
The abnormality analysis means analyzes the abnormality of the device based on the reference time domain signal having the same specific time interval and the same specific phase.

The abnormality analysis device includes a detection means, an abnormality analysis means, a storage section, and a reference time domain signal generation means.

The detection means is means for detecting and capturing the original time domain signal. The original time domain signal captured in this manner is used for abnormality analysis by the abnormality analysis means.
Further, the time domain signal obtained when the device is in a normal state is the normal time domain signal, and the signal obtained from the original time domain signal taken in at the time of diagnosis is the time domain signal to be diagnosed.

Now, the above-mentioned normal time domain signal has a specific time interval and a specific phase with respect to the original time domain signal. Here, the specific time interval is assumed to represent the frequency at which the device is operating. Furthermore, by setting the phase to match a specific phase between time-domain signals, it is possible to perform a rational abnormality analysis of a device whose operation is periodic in extracting feature amounts between time-domain signals.

The reference time domain signal generating means generates a reference time domain signal to be used as a diagnostic standard based on the sequentially stored normal time domain signals. In this generation, the time interval between the signals is set to a specific time interval, Since the phase is set as a specific phase, it is possible to obtain a signal generated at the frequency at which the device operates and a signal with less noise as a preferable signal to be used as a diagnostic standard in abnormality analysis.

Then, by obtaining the time domain signal to be diagnosed from the original time domain signal captured at the time of diagnosis, and performing abnormality analysis using the normal time domain signal that is the diagnostic standard, The accuracy of analysis is high by making use of time information. For example, when using an autoencoder as an abnormality analysis means, diagnosis can be performed by obtaining abnormal values with high sensitivity.

As a result, abnormality analysis can be performed with high sensitivity, and it becomes possible to satisfactorily analyze the degree of abnormality in, for example, an operation that involves combustion, such as an engine.

The second characteristic configuration of the present invention is
autocorrelation deriving means for determining the autocorrelation of the original time domain signal;
extracting means for extracting the time domain signal of the specific phase based on the autocorrelation derived by the autocorrelation deriving means,
The abnormality analysis means analyzes the abnormality of the device based on the reference time domain signal having the same specific time interval and the same specific phase.

In this configuration, the autocorrelation deriving means calculates the autocorrelation of the time domain signal. From the autocorrelation obtained in the time domain in this manner, the extraction means extracts time domain signals having the same specific phase between the signals, for example, as information between timings where the autocorrelation is strong. As a result, the abnormality analysis means can properly perform abnormality analysis on the signals having the same specific phase extracted in this way.

The third characteristic configuration of the present invention is
The extraction means extracts the time domain signal between the first extraction timing at which the autocorrelation is maximum and the second extraction timing at which the autocorrelation is the maximum from the original time domain signal at the specific time interval and The point is that it is extracted as a time domain signal of a specific phase.

If the original time domain signal has periodicity, taking the autocorrelation of the signal will calculate its period (for example, if the time domain signal is the sound or vibration emitted by an engine, the number of revolutions (rotational speed) of the engine) corresponds to this period)), the strength of autocorrelation appears depending on the period. Therefore, by sequentially extracting the timing at which the autocorrelation is maximum and the time-domain signals from that timing to the subsequent timing, it is possible to easily extract a signal with a specific phase from a periodic time-domain signal. I can do it. Here, the next timing does not need to be strictly continuous in one cycle, and may be, for example, at a timing interval of several cycles.

Furthermore, although we are talking about the specific phase between signals here, the time interval of the extracted time domain signals (the time interval from the first extraction timing to the second extraction timing mentioned above) is the information corresponding to the period. By making this time interval a specific time interval, the time intervals are unified.

The fourth characteristic configuration of the present invention is
The present invention further comprises an out-of-band signal eliminating means for eliminating, from the original time-domain signal, signal components that are outside a frequency band targeted for anomaly analysis.

For a specific device, it is usually known from experience in which frequency bands abnormalities are likely to occur.
Therefore, it is reasonable and wasteful to focus on this frequency band and perform abnormality analysis. As a result, the speed and reliability of diagnosis improves.
That is, by providing an out-of-band signal eliminating means and eliminating signal components outside the frequency band, a high-performance abnormality analysis device can be obtained.
Noise can also be removed from the reference time domain signal used as the basis for diagnosis, allowing for rational diagnosis.

The fifth characteristic configuration of the present invention is
The classification item of the reference time domain signal includes the operating frequency of the device,
The abnormality analysis means analyzes the abnormality of the equipment using different criteria according to the operating frequency.

The operating frequency of the equipment targeted by the present invention differs depending on the operating frequency.
For example, the frequency of operation at the rated operating frequency is high, and the frequency of operation at frequencies that temporarily pass during startup and stop is low. Furthermore, even at the time of starting and stopping, for example, at the time of starting, the frequency is operated in a relatively low frequency range for a while, and then the frequency is suddenly increased to the rated operating frequency (the opposite occurs when stopping). Therefore, the frequency of operation of the intermediate frequency region through which it passes during this rising process is extremely low.

In this way, between signals with different operating frequencies, the normal time domain signal obtained at the frequently operating frequency is highly reliable due to its extremely high frequency, while the normal time domain signal obtained at the infrequently operating frequency is highly reliable. Time domain oscillations are less reliable.

Therefore, in anomaly analysis, when performing an analysis based on a time domain signal during normal conditions on the high reliability side, the threshold value should be set to the side where the diagnostic evaluation is strict, and the threshold value on the low side to the side where it is lenient. By using different standards, rational anomaly analysis can be performed.

The sixth characteristic configuration of the present invention is
The abnormality analysis means calculates the degree of abnormality of the time domain signal to be diagnosed using an autoencoder constructed based on the reference time domain signal, and analyzes the abnormality of the device.

If this configuration is adopted, appropriate diagnosis can be easily and conveniently performed using the degree of abnormality obtained through processing in the autoencoder.

The seventh characteristic configuration of the present invention is:
The abnormality analysis means calculates the degree of abnormality of the time domain signal to be diagnosed using an autoencoder constructed based on the reference time domain signal classified by frequency, and analyzes the abnormality of the device. .

If this configuration is adopted, highly reliable abnormality analysis can be performed at each frequency by constructing an autoencoder for each frequency and performing abnormality analysis.

The eighth characteristic configuration of the present invention is:
The original time domain signal is sequentially detected by the detection means, and the autoencoder sequentially calculates the degree of abnormality based on the sequentially detected original time domain signal,
The problem lies in analyzing the abnormality of the equipment based on the change trend of the abnormality degree calculated sequentially.

By adopting this configuration, it is possible to perform a qualified abnormality analysis according to the change trend of the degree of abnormality.

The ninth characteristic configuration of the present invention is:
Analyzing equipment abnormalities based on the change trend of the abnormality degree,
The problem lies in analyzing the abnormality of the device by referring to the operation history, which is log data of the device.

If this configuration is adopted, an appropriate abnormality analysis can be performed based on the relationship between the operating history of the equipment and the degree of abnormality.

The tenth characteristic configuration of the present invention is:
The classification item of the reference time domain signal includes the operating frequency of the device,
The point is that the reference time domain signal is classified more finely in a region where the device is operated more frequently.

When this configuration is adopted, the above-mentioned divisions are made finely in areas with high operating frequency, and coarse divisions are made in areas with low operating frequency, so that the abnormality analysis between each division can be made uniform in accordance with the operating frequency. It can be made into something.

The eleventh characteristic configuration of the present invention is:
A frequency detection means for detecting a frequency at which the device operates, and a storage section for storing the frequency detected by the frequency detection means,
There is a point at which the abnormality analysis means analyzes an abnormality in the equipment between the time domain signals at specific time intervals that are calculated based on the frequency detected by the frequency detection means and have the same specific phase. .

When this configuration is adopted, a specific time interval is specified based on the frequency detected by the frequency detection means, and an abnormality is analyzed using a signal whose phase matches a specific phase among the time domain signals of the specific time interval. can.

Functional block diagram of an abnormality analysis device according to the present invention Diagram showing an example of an original time domain signal A diagram showing an example of the autocorrelation of the original time-domain signal shown in FIG. 2 and the form of capturing the time-domain signal. A diagram showing abnormality analysis results by the autoencoder according to the present invention Diagram showing an example of collecting operation frequency Diagram showing the encoding-decoding relationship at times with low abnormalities Diagram showing the encoding-decoding relationship at times with low abnormalities Diagram showing the encoding-decoding relationship at times with high abnormalities Diagram showing the state when the engine starts Diagram converted to 2D using principal component analysis Kernel density estimation diagram at engine startup shown in Figure 10 Functional block diagram of an abnormality analysis device according to another embodiment Explanatory diagram showing an example of the construction of an autoencoder in the abnormality analysis section Diagram showing the results of anomaly analysis by autoencoder performed in the frequency domain

Hereinafter, an abnormality analysis device 100 according to the present invention will be explained based on the drawings.
FIG. 1 is a functional block diagram of an abnormality analysis device 100 according to the present invention, and main functional means are described in the functional blocks.

The abnormality analysis device 100 is configured as an abnormality analysis device that diagnoses the device 1 while learning based on the original time domain signal S emitted by the device 1 that operates periodically. As an example of the device 1, an example in which the device is an engine is shown. As explained above, the engine operates periodically, and the frequency of the signals (sound, vibration, etc.) it emits in response to its start, stop, operating conditions, load, etc. changes. Furthermore, since the engine includes an internal combustion system, it is accompanied by events such as abnormal combustion that only occur for a very short period of time. This type of signal is a time domain signal and is detected by the sensor 2 as a detection means.
Therefore, the abnormality analysis device 100 detects sounds, vibrations, etc. emitted from this device 1, and performs abnormality analysis based on the original time domain signal S thereof.

The main components of the abnormality analysis device 100 include a preprocessing unit 10 that preprocesses the original time domain signal S, and stores the time domain signal s processed by the preprocessing unit 10 as a function of performing diagnosis with learning. an abnormality analysis section 30 that receives a time domain signal s (a first time domain signal s1 to be diagnosed) sent through the preprocessing section 10 when diagnosing the device 1 and performs abnormality analysis. It is configured.

As shown in the figure, the preprocessing unit 10 is a unit that filters the original time domain signal S and processes the filtered original time domain signal S according to a method unique to the present invention. This is a part that extracts, from the domain signal S, a time domain signal s having at least the same phase (specific phase in the present invention).

The time domain signal s extracted in this manner is sent to the storage section 20 and the abnormality analysis section 30.

Here, at the point of time when it is targeted for diagnosis, the time domain signal s sent from the abnormality analysis section 30 becomes the first time domain signal s1 that is targeted for diagnosis. On the other hand, the time-domain signal s sent to the storage unit 20 is associated with the log data and stored in the log-corresponding data storage unit 21, and data by operation frequency, which is one classification criterion, is attached to the log data. It is also stored in the operation frequency classified data storage unit 22 in association with the data storage unit 22.

The time domain signal s stored in the storage unit 20 serves as a reference in diagnosis, and the time domain signal s at a time when it is known to be normal in advance is used as the second time domain signal s2. This second time domain signal s2 is a normal time domain signal in the present invention.

In the example described below, the abnormality analysis unit 30 calculates the degree of abnormality for the first time domain signal s1(s) that is inputted one after another at regular time intervals.

Hereinafter, the preprocessing unit 10 (first preprocessing unit 11, second preprocessing unit 12), storage unit 20 (log corresponding data storage unit 21, operation frequency data storage unit 22), abnormality analysis unit 30 (abnormality degree calculation The section 30a and the abnormality determining section 30b) will be explained in this order.

As shown in FIG. 1, the preprocessing section 10 includes a first preprocessing section 11 and a second preprocessing section 12.
[First pretreatment section]
The first preprocessing unit 11 is a part that performs filter processing on the original time domain signal S, and the filter includes a low pass filter (LPF), a high pass filter (HPF), a band pass filter (BPF), or a band elimination filter (BEF). etc., in combination as appropriate. In FIG. 1, these filters are collectively expressed as out-of-band signal rejection means.

The reason why such filter processing is reasonable is that in the equipment 1 that is the target of anomaly analysis, the frequency at which an anomaly occurs (the frequency band where an anomaly is likely to occur) is often known in advance. This is because by appropriately eliminating signal components that fall outside the frequency band targeted for abnormality analysis, analysis can be performed that excludes unnecessary noise.

[Second pre-processing section]
The second preprocessing unit 12 includes an autocorrelation deriving means 12a constructed as software for signal processing, and a signal having the same phase from the autocorrelation of the time domain signal derived by the autocorrelation deriving means 12a. Extracting means 12b for extracting is provided.

FIG. 2 shows three examples of original time domain signals S (S1, S2, S3) whose acquisition times at specific dates and times differ at intervals of 0 and 01 seconds.
FIG. 3 shows the autocorrelation of the original time-domain signals S (S1, S2, S3) shown in FIG. 2 with different acquisition times using different line types. The autocorrelation R(Δt) is determined by the following equation 3.

R(Δt)=Σ{S(t)・S(t+Δt)}/{S(t)・S(t)}
Σ is the sum of t=1 to N Formula 3
Here, Δt is the time width, and N is the number of samples.

The extraction means 12b extracts the original time domain signal S between the first extraction timings t11, t21, t31 at which the autocorrelation is maximum and the second extraction timings t12, t22, t32 at which the autocorrelation is the maximum. It is extracted as a time domain signal s. The time domain signal s of the timing interval T1 (t11-t12), the time domain signal s of T2 (t21-t22), and the time domain signal s of T3 (t31-t32) shown in the upper left of FIG. This becomes a region signal s. In the present invention, this timing interval is called a specific time interval, and this specific time interval is information corresponding to the frequency when the device 1 is operated.

When adopting such an extraction form, since the time interval between the first extraction timing t11, t21, t31 and the second extraction timing t12, t22, t32 is the same as, for example, 0.067 sec, the time domain signal s Its frequency does not change between. Further, since the timing at which the autocorrelation is maximum is set as the start point and end point of extraction, the phases of the signals (for example, T1 and T2) match. Therefore, in this case, abnormality determination can also be made between these time domain signals s and s. That is, time domain signals s having the same specific time interval and specific phase are extracted.

However, basically, the time domain signal to be diagnosed (the first time domain signal s1) is taken in, and the reference time domain signal (the second time domain signal s2, which is the normal time in the present invention). Since the acquisition time of the domain signal) is far apart, abnormality determination is performed between the time domain signals s1 and s2 having the same specific time interval corresponding to the frequency and the same specific phase.

In the case of the embodiment shown in FIG. 1, frequency data is stored on the log data side, and abnormality determination is performed on both sides using the time domain signal s that satisfies the same frequency and phase conditions as the time domain signal s1 at the time of diagnosis. be able to.

[Storage section]
The configuration of the storage unit 20 that mainly stores such a time domain signal (second time domain signal s2) during normal operation will be described below.

As shown in FIG. 1, the storage unit 20 includes a log-corresponding data storage unit 21 and an operation frequency-based data storage unit 22.

The information referred to as data here is obtained by preprocessing the original time domain signal S that is collected sequentially, and classifying and averaging the time by associating it with predetermined information (log data and operation frequency corresponding to this log data). This is the area signal s. This signal is a reference time domain signal used as a diagnostic standard, and this signal is generated by a reference time domain signal generating means.

As is clear from FIG. 1, the storage unit 20 sequentially stores the log data LD detected according to the operation of the device 1, the time domain signal s processed by the second preprocessing unit 12, and the abnormality determination unit 30b. The judgment result is input.

Specifically, the log data LD includes the target date and time, the engine rotation speed at the date and time (corresponding to the frequency of the time domain signal s), intake air temperature, engine load factor, fuel gas pressure, and the like.
The time domain signal s is a signal obtained through preprocessing from the original time domain signal S previously explained with reference to FIG.

Since the degree of abnormality obtained sequentially and the abnormality determination result based on the degree of abnormality are also input to the storage unit 20, the attribute of whether the past data to be stored is normal or abnormal. It can be classified according to As a result, it is possible to accumulate and average past data, at least regarding the time domain signal s during normal times (second time domain signal s2).

[Log compatible data storage unit]
In the log corresponding data storage unit 21, the time domain signal s accumulated as described above is stored in association with various log data LD (engine speed, intake air temperature, engine load factor, fuel gas pressure, etc.). By specifying the various log data LD, the time domain signal s (averaged signal) associated with the log data LD can be extracted as the second time domain signal s2. Naturally, once a frequency is identified, a second time domain signal s2 of that frequency can be derived based on the corresponding engine speed.
As a result, the information stored in this log-corresponding data storage unit 21 (a second time-domain signal during normal operation having the same specific time interval corresponding to frequency and specific phase as the first time-domain signal s1 to be diagnosed) s2) can be used for abnormality determination by the abnormality determination section 30b.

[Data storage section by operation frequency]
The abnormality analysis device 100 separately collects the operating frequencies of the various classification items mentioned above. Therefore, the data classified and accumulated in the log corresponding data storage section 21 is associated with the operating frequency of that classification. FIG. 5 shows an example in which the horizontal axis is the output, the vertical axis is the rotational speed (rotational speed; corresponds to the frequency explained so far), and the operating frequency is plotted in the height direction. In this example, it can be seen that the operating frequency is high within the range where the output and rotation speed satisfy a certain linear relationship. Therefore, this storage unit 20 includes time-domain signals s during normal times, classified by frequencies at which the device 1 is operated at different operating frequencies.

[Anomaly Analysis Department]
The abnormality analysis section 30 includes an abnormality degree calculation section 30a and an abnormality determination section 30b. This abnormality analysis section 30 serves as abnormality analysis means in the present invention.

The example shown in FIG. 1 is an example in which anomaly analysis is performed using a so-called autoencoder while sequentially calculating the degree of anomaly.

[Construction of autoencoder]
The construction of the autoencoder is executed when the device 1 is in a normal state, such as when shipped from the factory, for example. That is, the output signal from the sensor 2 is taken in in a normal state, and the time domain signal s having the same phase is obtained through the first preprocessing section 11 and the second preprocessing section 12. At this time, by not changing the frequency, time domain signals s having the same phase can be automatically captured sequentially. Then, using such time domain signal s, an autoencoder is initially constructed according to the method described above. At the time of diagnosis, by inputting the time domain signal s to be diagnosed as an input to this autoencoder, the time domain signal s to be diagnosed is calculated from the output obtained through dimensional compression (encoding) and restoration (decoding) processing. You can get the degree of abnormality.

Since the autoencoder is constructed based on the normal time domain signal s that is classified and stored at least for each frequency (for each specific time interval), a separate autoencoder is constructed for each frequency (for example, every 1 Hz). Anomalies can be analyzed by calculating the degree of anomaly shown below using a separate autoencoder for each frequency. This configuration is shown in FIG.

As for the divisions for each frequency, divisions according to the operating frequency of the equipment may be adopted. As explained above, the storage format in the storage unit 20 is to store the frequency and the second time domain signal s2 that corresponds to the frequency and is associated with the operating frequency. By making small divisions, making coarse divisions in areas where operation frequency is low, and building an autoencoder as an abnormality analysis means for each division, the accuracy of abnormality analysis can be improved and the analysis results can be made uniform.

In addition, the construction of this autoencoder can be performed at the time of factory shipment as described above, and since the abnormality analysis device 100 according to the present invention sequentially performs abnormality determination, the determination results are used to determine the time range determined to be normal. A configuration that uses the signal s to sequentially update the autoencoder based on an averaged signal obtained by averaging a plurality of time-domain signals s having the same frequency and phase, and performs highly reliable abnormality determination. shall be.

[Abnormality calculation unit]
The abnormality degree calculation unit 30a is a part that calculates the abnormality degree from the autoencoder constructed as described above and the relationship between its input and output.

FIG. 4 shows the results of abnormality degree calculations calculated according to the above method over a specific 5-day period.
As can be seen from this result, the level of abnormality has increased and the sensitivity has increased with respect to the abnormality degree calculation results shown in FIG. 14.

[Abnormality determination section]
The abnormality determination unit 30b determines the degree of abnormality one by one, and determines that the abnormality is abnormal when the degree of abnormality becomes excessively large.

Examples of deriving the degree of abnormality and determining abnormality will be described below using specific examples shown in FIGS. 6 to 8.
These drawings show the relationship between the input (denoted as input in the diagram) and the output (described as output in the diagram) when the time domain signal s is input/output processed using an autoencoder obtained in advance. ing. The obtained acquisition date and time and degree of abnormality are shown in the upper part of the figure. Further, the upper right of the figure is a reduced view of FIG. 4, and the thick circles indicate the times targeted for abnormality determination.

6 and 7 are diagrams showing the input/output relationship of time domain signals at times when the degree of abnormality is low, and FIG. 8 is a diagram showing the input/output relationships at times when the degree of abnormality is high.

Simply organized and written,
Figure 6 ; Abnormality degree 0.0432
Figure 7 ; Abnormality degree 0.0504
Figure 8 ; Abnormality degree 0.104
It becomes.
It can be seen from FIG. 4 and these figures that, for example, an abnormality determination can be made using an abnormality degree of 0.1 as a threshold value.

Although the above determination is based on a certain threshold value, it is also possible to sequentially calculate the degree of abnormality at intervals of 0.01 seconds and determine that the degree of abnormality is abnormal if the degree of abnormality shows an excessively increasing tendency.

For example, if the abnormality degree calculation interval is 0.01 seconds, a longer time interval, for example, about 5 seconds, is set, and during this time interval, the abnormality degree shows an increasing trend, does not change significantly, or decreases. It is also possible to determine that there is an abnormality only when the abnormality degree shows an increasing tendency, based on the change tendency of the abnormality degree such as shown in FIG. In this case, it is preferable to analyze the abnormality by comparing the operating history of the device 1 and the change trend of the degree of abnormality.

[Weighting of abnormality degree judgment based on operation frequency]
The abnormality degree calculation unit 30a and the abnormality determination unit 30b perform abnormality analysis using different criteria by changing the abnormality degree calculation standard according to the operation frequency based on the time domain signal s for each normal time and operation frequency. It is composed of

In the examples shown in FIGS. 4 to 8, the time domain signal s is detected at the time when the device 1 continues to operate and is operated at a high frequency. Therefore, since the accuracy of diagnosis in this state is high, the threshold value of the degree of abnormality can be determined as 1.00 as shown above.

On the other hand, when operating in a state where the operating frequency is low, the accuracy of diagnosis decreases due to reasons such as a small number of samplings, and the probability of erroneously determining abnormality increases. Therefore, in such a case, it is preferable to set the abnormality degree threshold to a low value (for example, 0.75).

Regarding such a weighting method for determinations related to operating frequency, an example in which weighting is performed using the concept of dimensional compression described in connection with the autoencoder will be specifically described below. In the explanation, the situation when the engine is started will be explained as an example.
FIG. 9 shows the state of such an engine at startup.
In the same figure, the solid line indicates "output," the dashed line indicates "fuel valve opening," and the broken line indicates "rotational speed." The horizontal axis is the elapsed time from the start of activation. This figure shows the state from the start of startup (0 sec) to the time until the output reaches a predetermined value (340 sec). In the operation of the engine, exhaust gas flow rate, exhaust gas temperature, and other factors are also involved (influencing factors), but these factors will be omitted in the explanation.

In dimensional compression, the eigenvectors of principal component analysis are obtained and compression is performed.
The first principal component PCA1 and the second principal component PCA2, which are principal components, are described as follows using the degrees of contribution Wi, Wj of each influencing element x _i, x _j . Here, bi and cj are appropriately determined according to a multivariate analysis method.

First principal component PCA1; x=ΣW _i x _i +b _i Equation 4
Second principal component PCA2; y=ΣW _j x _j +c _j Equation 5

FIG. 10 shows the relationship between the eigenvectors (PCA1, PCA2) obtained in this way, and in this figure, the event at engine startup follows a trajectory of Z1→Z2→Z3 as time passes.

Figure 11 shows the estimated kernel density value calculated by principal component analysis after reducing the dimensionality of the elements during the operating state until 5 hours after engine startup using principal component analysis, and the probability of occurrence thereof. It is. In the figure, the darker side is the side with higher density. For reference, the density distribution is mapped to height and displayed on the right and upper sides of the figure. Therefore, as can be seen from this figure, the weighting described above can be made to be the estimated kernel density value x α, and the separation interval can be calculated as β/estimated kernel density value. Here, α and β are set for each target.

As explained earlier, for areas with a high probability of appearance, the degree of abnormality calculated for this area is used almost as is, and for areas with a low probability of appearance, the degree of abnormality calculated for this area is used as is. can be weighted.

As explained above, in the abnormality analysis device 100 according to the present invention, the second time domain signal s2, which is the time domain signal at the normal time determined in advance (time domain signal at the normal time), and the time domain signal at the time of diagnosis. Abnormality analysis is to be performed in comparison with the first time domain signal s1.

When this configuration is adopted, analysis is performed based on a normal time domain signal that has the same frequency and phase as the time domain signal at the time of diagnosis at a specific time interval and a specific phase. When using a normal time domain signal, obtain an averaged signal of multiple time domain signals s2 and use this as a reference time domain signal. Can perform anomaly analysis.

[Another embodiment]
(1) In the above embodiment, the frequency at which the device is operating is taken in from the log data LD side, and the phase of the time domain signal s is adjusted using the autocorrelation deriving means 12a and the extraction means 12b. However, as shown in FIG. 3, since the signal length (time width) of the time domain signal s according to the present invention is specified, this time width can be used to determine the length of the signal. The frequency at the time of capture can be determined. When the signal time length is Tsec, the frequency is 60/T. Therefore, in another embodiment shown in FIG. 12, the frequency calculation means 12c as a frequency detection means provided in the second preprocessing section performs the calculation.
By doing this, both information regarding frequency and phase is obtained from the original time domain signal S, and the identity of the specific time interval and specific phase corresponding to the frequency and phase in the time domain signal s is ensured, thereby increasing reliability. It is possible to perform highly accurate anomaly analysis.

(2) In the above embodiment, in the anomaly analysis accompanied by learning, the anomaly analysis method executed by the anomaly analysis means is a method of constructing an autoencoder and calculating the degree of anomaly to perform the anomaly analysis. The main focus of the invention is to perform anomaly analysis between signals that correspond to frequency and phase and have a specific time interval and a specific phase in time domain signals. Abnormality analysis methods include methods using autoencoders and other methods. , one-class SWN, etc. may also be adopted.

(3) In the above embodiment, the timings at which the autocorrelation becomes maximum are sequentially used to match the phases between time domain signals. It is possible to adopt various other methods in which the zero cut point on the increasing side of the differential value of is used as the starting point of extraction.

(4) In the above embodiment, an engine is used as an example of the device 1, but abnormality analysis of a gas turbine or the like can also be performed.

1 Engine (equipment)
10 Pre-processing section (pre-processing means)
11 First preprocessing section (out-of-band signal elimination means)
12 Second preprocessing section 12a Autocorrelation deriving means 12b Extracting means 12c Frequency calculation means (frequency detection means)
20 Storage unit 21 Log corresponding data storage unit 22 Data storage unit by operation frequency 30 Abnormality analysis unit (abnormality analysis means)
30a Abnormality degree calculation unit (abnormality degree calculation means)
30b Abnormality determination section (abnormality determination means)
S Original time domain signal s Time domain signal s1 First time domain signal (time domain signal)
s2 Second time domain signal (time domain signal/normal time domain signal)
LD log data

Claims (10)

detection means for detecting raw time domain signals emitted by periodically operating equipment;
an abnormality analysis means for analyzing an abnormality of the device based on the original time domain signal detected by the detection means,
a storage unit for storing a time domain signal at a specific time interval and a specific phase detected by the detection means and extracted from the original time domain signal emitted by the normal device as a normal time domain signal;
A reference time domain signal generating means for generating a reference time domain signal as a diagnostic standard based on the frequency of the device and corresponding to the frequency based on the normal time domain signals sequentially stored in the storage unit,
the reference time domain signal is classified based on the frequency at which the equipment operates;
Regarding the original time domain signal to be diagnosed which the detection means detects based on the frequency of the equipment and corresponding to the frequency, the abnormality analysis means detects the reference time domain signal in which the specific time interval and the specific phase corresponding to the frequency are the same. A configuration that analyzes equipment abnormalities by determining abnormalities based on signals.
An abnormality analysis device, wherein the abnormality analysis means analyzes the abnormality of the equipment based on different criteria according to the operating frequency with respect to frequencies of the equipment having different operating frequencies . detection means for detecting raw time domain signals emitted by periodically operating equipment;
an abnormality analysis means for analyzing an abnormality of the device based on the original time domain signal detected by the detection means,
a storage unit for storing a time domain signal at a specific time interval and a specific phase detected by the detection means and extracted from the original time domain signal emitted by the normal device as a normal time domain signal;
A reference time domain signal generating means for generating a reference time domain signal as a diagnostic standard based on the frequency of the device and corresponding to the frequency based on the normal time domain signals sequentially stored in the storage unit,
the reference time domain signal is classified based on the frequency at which the equipment operates;
Regarding the original time domain signal to be diagnosed which the detection means detects based on the frequency of the equipment and corresponding to the frequency, the abnormality analysis means detects the reference time domain signal in which the specific time interval and the specific phase corresponding to the frequency are the same. A configuration that analyzes equipment abnormalities by determining abnormalities based on signals.
An abnormality analysis device, wherein the abnormality analysis means determines an abnormality based on the reference time domain signal that is classified more finely in a region where the operation frequency of the equipment is higher, with respect to frequencies of the equipment whose operating frequencies are different . autocorrelation deriving means for determining the autocorrelation of the original time domain signal;
extracting means for extracting the time domain signal of the specific phase based on the autocorrelation derived by the autocorrelation deriving means,
3. The abnormality analysis device according to claim 1 , wherein the abnormality analysis means analyzes the abnormality of the device based on the reference time domain signal having the same specific time interval and the same specific phase. The extraction means extracts the time domain signal between the first extraction timing at which the autocorrelation is maximum and the second extraction timing at which the autocorrelation is the maximum from the original time domain signal at the specific time interval and The abnormality analysis device according to claim 3 , wherein the abnormality analysis device extracts the signal as a time domain signal of a specific phase. The abnormality analysis apparatus according to any one of claims 1 to 4 , further comprising an out-of-band signal eliminating means for eliminating signal components outside a frequency band targeted for abnormality analysis from the original time domain signal. Any one of claims 1 to 5, wherein the abnormality analysis means analyzes the abnormality of the device by calculating the degree of abnormality of the time domain signal to be diagnosed using an autoencoder constructed based on the reference time domain signal. Anomaly analysis device as described in section. 2. The abnormality analysis means calculates the degree of abnormality of the time domain signal to be diagnosed using an autoencoder constructed based on the reference time domain signal classified by frequency, and analyzes the abnormality of the device. The abnormality analysis device according to any one of items 6 to 6. The original time domain signal is sequentially detected by the detection means, and the autoencoder sequentially calculates the degree of abnormality based on the sequentially detected original time domain signal,
The abnormality analysis device according to claim 6 or 7, wherein the abnormality analysis device analyzes the abnormality of the device based on the change tendency of the abnormality degree calculated sequentially. Analyzing equipment abnormalities based on the change trend of the abnormality degree,
In a time interval longer than the above-mentioned abnormality degree calculation interval, the change trend of the abnormality degree is determined: whether the abnormality degree shows an increasing tendency, does not change significantly, or shows a decreasing tendency, and if the abnormality degree shows an increasing tendency, an abnormality is detected. 9. The abnormality analysis device according to claim 8, which analyzes if there is an abnormality. A frequency detection means for detecting a frequency at which the device operates, and a storage section for storing the frequency detected by the frequency detection means,
2. The abnormality analysis means analyzes the abnormality of the equipment between time domain signals having a specific phase, which are calculated based on the frequency detected by the frequency detection means at specific time intervals, and which have the same specific phase. The abnormality analysis device according to any one of items 1 to 9 .

Images