The present invention relates to an abnormality detection device and an abnormality detection method for determining whether or not equipment is operating abnormally.
A conventional abnormality detection method for detecting abnormality of equipment compares abnormality detection time-series data indicating states of equipment at a plurality of times in time series with time-series data at normal time collected when the equipment is operating normally.
The conventional abnormality detection method detects abnormality of equipment by detecting time-series data of a part whose behavior is different from that of time-series data at normal time (hereinafter, referred to as “subsequence data”) from among pieces of abnormality detection time-series data.
However, the subsequence data is time-series data in a time period in which abnormality may have occurred in the equipment, but abnormality has not necessarily occurred in the equipment, and the equipment may be operating normally.
The following Patent Literature 1 discloses an abnormality detection system for detecting abnormality of equipment by combining a conventional abnormality detection method and a method for analyzing event information in order to avoid occurrence of an erroneous determination indicating that abnormality has occurred in the equipment when the equipment is operating normally.
Examples of the event information include information indicating an event related to operation of equipment by a worker and information indicating an event related to replacement of parts of the equipment.
The abnormality detection system disclosed in Patent Literature 1 determines that no abnormality has occurred in equipment even when detecting subsequence data as long as the detected subsequence data is synchronized with an event indicated by event information.
Patent Literature 1: JP 2013-218725 A
The abnormality detection system disclosed in Patent Literature 1 needs to hold event information in advance.
In a case where the abnormality detection system disclosed in Patent Literature 1 cannot prepare event information in advance, when the abnormality detection system detects subsequence data while equipment is operating normally, the abnormality detection system erroneously determines that abnormality has occurred in the equipment disadvantageously.
The present invention has been achieved in order to solve the above-described problem, and an object of the present invention is to obtain an abnormality detection device and an abnormality detection method capable of avoiding occurrence of an erroneous determination indicating that abnormality has occurred in equipment without preparing event information in advance.
An abnormality detection device according to the present invention includes: processing circuitry to calculate, from abnormality detection time-series data indicating states of equipment which is an abnormality detection target at a plurality of times in time series, the degree of abnormality of the equipment at each of the plurality of times as an abnormality detection outlier score; to extract, from among pieces of the abnormality detection time-series data, a piece of abnormality detection time-series data in a time period in which an abnormality may have occurred in the equipment as abnormality detection outlier data on the basis of the abnormality detection outlier score at each of the plurality of times; to collate a waveform of the abnormality detection outlier data with a waveform condition for determining that a waveform indicating a change in the abnormality detection outlier data is a waveform obtained when the equipment is operating normally, to determine whether or not the equipment is operating abnormally on the basis of a collation result between the waveform condition and the waveform of the abnormality detection outlier data; to calculate a feature amount of the abnormality detection outlier data, and to determine a waveform type of the abnormality detection outlier data from the feature amount; to select a waveform condition corresponding to the type from among one or more waveform conditions; to collate the waveform condition with the waveform of the abnormality detection outlier data, and to determine whether or not the equipment is operating abnormally on a basis of a collation result between the selected waveform condition and the waveform of the abnormality detection outlier data; to calculate, from each of one or more pieces of learning time-series data indicating states of the equipment at a plurality of times when the equipment is operating normally in time series, a degree of abnormality of the equipment at each of the plurality of times as a learning outlier score; to extract, from among the pieces of learning time-series data, learning time-series data in a time period in which an abnormality may have occurred in the equipment as learning outlier data on a basis of the learning outlier score at each of the plurality of times; to calculate a feature amount of each of the pieces of learning outlier data, and to determine a waveform type of each of the pieces of learning outlier data from the feature amount of each of the pieces of learning outlier data, and to generate, from among waveforms of one or more pieces of learning outlier data whose waveforms have been determined to be of the same type out of the pieces of learning outlier data, a waveform condition corresponding to the type.
According to the present invention, the abnormality detection device is configured in such a manner that the processing circuitry collates a waveform of the abnormality detection outlier data with a waveform condition for determining that a waveform indicating a change in the abnormality detection outlier data is a waveform obtained when the equipment is operating normally, and determines whether or not the equipment is operating abnormally on the basis of a collation result between the waveform condition and the waveform of the abnormality detection outlier data. Therefore, the abnormality detection device according to the present invention can avoid occurrence of erroneous determination indicating that an abnormality has occurred in the equipment without preparing event information in advance.
Hereinafter, in order to describe the present invention in more detail, embodiments for performing the present invention will be described by referring to the attached drawings.
In
The learning data inputting unit 1 receives input of N pieces of learning time-series data DG,n,t (n=1, 2, . . . , N) indicating states of equipment which is an abnormality detection target at a plurality of times tin time series when the equipment is operating normally. N is an integer equal to or more than 1.
The learning time-series data DG,n,t includes an observed value of a sensor at each time t, and the observed value of the sensor indicates a state of the equipment.
The learning data inputting unit 1 outputs the received learning time-series data DG,n,t to each of an outlier score calculating unit 3 and an outlier data extraction processing unit 7.
As the equipment which is an abnormality detection target, equipment such as a power plant, a chemical plant, or a water and sewage plant is conceivable. In addition, as the equipment which is an abnormality detection target, air conditioning equipment, electrical equipment, lighting equipment, water supply and drainage equipment, or the like in an office building or a factory is conceivable. In addition, equipment such as a conveyor constituting a production line of a factory, equipment installed in an automobile, or equipment installed in a railway vehicle is conceivable. Furthermore, as the equipment which is an abnormality detection target, equipment of an information system related to economy or equipment of an information system related to management is also conceivable.
An abnormality detection data inputting unit 2 is achieved by, for example, an input interface circuit 22 illustrated in
The abnormality detection data inputting unit 2 receives input of abnormality detection time-series data DU,t indicating states of equipment which is an abnormality detection target at a plurality of times tin time series.
The abnormality detection time-series data DU,t includes an observed value of a sensor at each time t, and the observed value of the sensor indicates a state of the equipment.
The abnormality detection data inputting unit 2 outputs the received abnormality detection time-series data DU,t to each of the outlier score calculating unit 3 and the outlier data extraction processing unit 7.
The outlier score calculating unit 3 is achieved by, for example, an outlier score calculating circuit 23 illustrated in
The outlier score calculating unit 3 calculates the degree of abnormality of the equipment at each time t as a learning outlier score SG,n,t from each of the N pieces of learning time-series data DG,n,t output from the learning data inputting unit 1. The outlier score calculating unit 3 outputs the calculated learning outlier score SG,n,t at each time t to an outlier data extracting unit 4.
The outlier score calculating unit 3 calculates the degree of abnormality of the equipment at each time t as an abnormality detection outlier score SU,t from the abnormality detection time-series data DU,t output from the abnormality detection data inputting unit 2. The outlier score calculating unit 3 outputs the calculated abnormality detection outlier score SU,t at each time t to the outlier data extracting unit 4.
The outlier data extracting unit 4 includes a threshold calculating unit 5, a threshold storing unit 6, and the outlier data extraction processing unit 7.
The outlier data extracting unit 4 extracts time-series data in a time period in which an abnormality may have occurred in the equipment as learning outlier data ODG,n from among pieces of the learning time-series data DG,n,t on the basis of the learning outlier score SG,n,t calculated by the outlier score calculating unit 3. The outlier data extracting unit 4 outputs the extracted learning outlier data ODG,n to each of an abnormality determining unit 8 and a waveform condition generating unit 12.
The outlier data extracting unit 4 extracts abnormality detection time-series data in a time period in which an abnormality may have occurred in the equipment as abnormality detection outlier data ODU,ts′-te′ from among pieces of the abnormality detection time-series data DU,t on the basis of the abnormality detection outlier score SU,t calculated by the outlier score calculating unit 3. The outlier data extracting unit 4 outputs the extracted abnormality detection outlier data ODU,ts′-te′ to the abnormality determining unit 8.
The threshold calculating unit 5 is achieved by, for example, a threshold calculating circuit 24 illustrated in
The threshold calculating unit 5 calculates a threshold Sth from the learning outlier score SG,n,t calculated by the outlier score calculating unit 3, and outputs the threshold Sth to the threshold storing unit 6.
The threshold storing unit 6 is achieved by, for example, a threshold storing circuit 25 illustrated in
The threshold storing unit 6 stores the threshold Sth output from the threshold calculating unit 5.
The outlier data extraction processing unit 7 is achieved by, for example, an outlier data extraction processing circuit 26 illustrated in
The outlier data extraction processing unit 7 compares the learning outlier score SG,n,t calculated by the outlier score calculating unit 3 at each time t with the threshold Sth stored by the threshold storing unit 6.
The outlier data extraction processing unit 7 extracts learning outlier data ODG,n,ts-te from among pieces of the learning time-series data DG,n,t on the basis of a comparison result between the learning outlier score SG,n,t at each time t and the threshold Sth. The outlier data extraction processing unit 7 outputs the extracted learning outlier data ODG,n,ts-te to each of a type determining unit 9, a waveform condition selecting unit 10, a waveform classifying unit 13, and a waveform condition generation processing unit 14.
The outlier data extraction processing unit 7 compares the abnormality detection outlier score SU,t calculated by the outlier score calculating unit 3 at each time t with the threshold Sth stored by the threshold storing unit 6.
The outlier data extraction processing unit 7 extracts abnormality detection outlier data ODU,ts′-te′ from among pieces of the abnormality detection time-series data DU,t on the basis of a comparison result between the abnormality detection outlier score SU,t at each time t and the threshold Sth. The outlier data extraction processing unit 7 outputs the extracted abnormality detection outlier data ODU,ts′-te′ to each of the type determining unit 9, the waveform condition selecting unit 10, and an abnormality determination processing unit 11.
The abnormality determining unit 8 includes the type determining unit 9, the waveform condition selecting unit 10, and the abnormality determination processing unit 11.
The abnormality determining unit 8 collates a waveform condition Wp with a waveform of the abnormality detection outlier data ODU,ts′-te′ extracted by the outlier data extracting unit 4. The waveform condition Wp is a condition for determining that a waveform indicating a change in the abnormality detection outlier data ODU,ts′-te′ extracted by the outlier data extracting unit 4 is a waveform obtained when the equipment is operating normally.
The abnormality determining unit 8 determines whether or not the equipment is operating abnormally on the basis of a collation result between the waveform condition Wp and the waveform of the abnormality detection outlier data ODU,ts′-te′, and outputs a determination result indicating whether or not the equipment is operating abnormally to a detection result outputting unit 16.
The type determining unit 9 is achieved by, for example, a type determining circuit 27 illustrated in
The type determining unit 9 calculates a feature amount CG,n of the learning outlier data ODG,n,ts-te extracted by the outlier data extraction processing unit 7, and determines the waveform type of the learning outlier data ODG,n,ts-te from the feature amount CG,n. The type determining unit 9 outputs the determined waveform type of the learning outlier data ODG,n,ts-te to the waveform classifying unit 13.
The type determining unit 9 calculates a feature amounts Cu of the abnormality detection outlier data ODU,ts′-te′ extracted by the outlier data extraction processing unit 7, and determines the waveform type of the abnormality detection outlier data ODU,ts′-te′ from the feature amount Cu. The type determining unit 9 outputs the determined waveform type of the abnormality detection outlier data ODU,ts′-te′ to the waveform condition selecting unit 10.
The waveform condition selecting unit 10 is achieved by, for example, a waveform condition selecting circuit 28 illustrated in
The waveform condition selecting unit 10 selects a waveform condition Wp corresponding to the type determined by the type determining unit 9 from among one or more waveform conditions Wp stored in a waveform condition storing unit 15, and outputs the selected waveform condition Wp to the abnormality determination processing unit 11.
The abnormality determination processing unit 11 is achieved by, for example, an abnormality determination processing circuit 29 illustrated in
The abnormality determination processing unit 11 collates the waveform condition Wp selected by the waveform condition selecting unit 10 with the waveform of the abnormality detection outlier data ODU,ts′-te′ extracted by the outlier data extraction processing unit 7.
The abnormality determination processing unit 11 determines whether or not the equipment is operating abnormally on the basis of a collation result between the waveform condition Wp and the waveform of abnormality detection outlier data ODU,ts′-te′, and outputs a determination result indicating whether or not the equipment is operating abnormally to the detection result outputting unit 16.
The waveform condition generating unit 12 includes the waveform classifying unit 13, the waveform condition generation processing unit 14, and the waveform condition storing unit 15.
The waveform condition generating unit 12 generates, from waveforms of one or more pieces of learning outlier data ODG,n,ts-te whose waveforms have been determined to be of the same type by the type determining unit 9 out of the pieces of learning outlier data ODG,n,ts-te extracted by the outlier data extracting unit 4, a waveform condition corresponding to the type. The waveform condition generating unit 12 stores the generated waveform condition.
The waveform classifying unit 13 is achieved by, for example, a waveform classifying circuit 30 illustrated in
The waveform classifying unit 13 calculates the degree of similarity between one or more pieces of learning outlier data ODG,n,ts-te whose waveforms have been determined to be of the same type by the type determining unit 9 out of the pieces of learning outlier data ODG,n,ts-te extracted by the outlier data extracting unit 4.
The waveform classifying unit 13 classifies one or more pieces of learning outlier data ODG,n,ts-te whose waveforms have been determined to be of the same type by the type determining unit 9 into groups on the basis of the calculated degree of similarity.
The waveform classifying unit 13 outputs a classification result of one or more pieces of learning outlier data ODG,n,ts-te to the waveform condition generation processing unit 14.
The waveform condition generation processing unit 14 is achieved by, for example, a waveform condition generation processing circuit 31 illustrated in
The waveform condition generation processing unit 14 generates, for each of the groups provided by the waveform classifying unit 13, a waveform condition Wp corresponding the group from the waveforms of the one or more pieces of learning outlier data ODG,n,ts-te classified into the same group by the waveform classifying unit 13. The waveform condition generation processing unit 14 outputs the generated waveform condition Wp to the waveform condition storing unit 15.
The waveform condition storing unit 15 is achieved by, for example, a waveform condition storing circuit 32 illustrated in
The waveform condition storing unit 15 stores the waveform condition Wp generated by the waveform condition generation processing unit 14.
The detection result outputting unit 16 is achieved by, for example, a detection result outputting circuit 33 illustrated in
The detection result outputting unit 16 displays the determination result output from the abnormality determination processing unit 11 on, for example, a display (not illustrated).
In
Here, for example, to each of the threshold storing circuit 25 and the waveform condition storing circuit 32, a nonvolatile or volatile semiconductor memory such as random access memory (RAM), read only memory (ROM), flash memory, erasable programmable read only memory (EPROM), or electrically erasable programmable read only memory (EEPROM), a magnetic disk, a flexible disk, an optical disc, a compact disc, a mini disc, or a digital versatile disc (DVD) is applicable.
For example, to each of the input interface circuit 21, the input interface circuit 22, the outlier score calculating circuit 23, the threshold calculating circuit 24, the outlier data extraction processing circuit 26, the type determining circuit 27, the waveform condition selecting circuit 28, the abnormality determination processing circuit 29, the waveform classifying circuit 30, the waveform condition generation processing circuit 31, and the detection result outputting circuit 33, a single circuit, a composite circuit, a programmed processor, a parallel programmed processor, an application specific integrated circuit (ASIC), a field-programmable gate array (FPGA), or a combination thereof is applicable.
The constituent elements of the abnormality detection device are not limited to those achieved by dedicated hardware, and the abnormality detection device may be achieved by software, firmware, or a combination of software and firmware.
The software or the firmware is stored as a program in a memory of a computer. The computer means hardware for executing a program. For example, to the computer, a central processing unit (CPU), a central processing device, a processing device, an arithmetic device, a microprocessor, a microcomputer, a processor, or a digital signal processor (DSP) is applicable.
When the abnormality detection device is achieved by software, firmware, or the like, the threshold storing unit 6 and the waveform condition storing unit 15 are configured on a memory 41 of a computer. A program for causing the computer to execute a processing procedure performed in the learning data inputting unit 1, the abnormality detection data inputting unit 2, the outlier score calculating unit 3, the threshold calculating unit 5, the outlier data extraction processing unit 7, the type determining unit 9, the waveform condition selecting unit 10, the abnormality determination processing unit 11, the waveform classifying unit 13, the waveform condition generation processing unit 14, and the detection result outputting unit 16 is stored in the memory 41. A processor 42 of the computer executes the program stored in the memory 41.
Next, an operation of the abnormality detection device illustrated in
First, an operation during learning in the abnormality detection device will be described.
First, the learning data inputting unit 1 receives input of N pieces of learning time-series data DG,n,t (n=1, 2, . . . , N) indicating states of equipment which is an abnormality detection target at a plurality of times tin time series when the equipment is operating normally (step ST1 in
The learning data inputting unit 1 outputs the received learning time-series data DG,n,t to each of the outlier score calculating unit 3 and the outlier data extracting unit 4.
In
When receiving N pieces of learning time-series data DG,n,t from the learning data inputting unit 1, the outlier score calculating unit 3 calculates the degree of abnormality of the equipment at each time t as a learning outlier score SG,n,t from each of the N pieces of learning time-series data DG,n,t (step ST2 in
A known technique is applied to a process for calculating the learning outlier score SG,n,t. For example, the following Non-Patent Literature 1 discloses a process for calculating an outlier score. The “Matrix Profile” disclosed in Non-Patent Literature 1 corresponds to an outlier score.
Non-Patent Literature 1:
Chin-Chia Michael Yeh, Yan Zhu, Liudmila Ulanova, Nurjahan Begum, Yifei Ding, Hoang Anh Dau, Diego Furtado Silva, Abdullah Mueen, Eamonn Keogh (2016). Matrix Profile I: All Pairs Similarity Joins for Time Series: A Unifying View that Includes Motifs, Discords and Shapelets.
In the abnormality detection device illustrated in
However, this is only an example, and for example, the outlier score calculating unit 3 may calculate a residual between an observed value of a sensor at each time t included in the learning time-series data DG,n,t and a predicted value at time t as the learning outlier score SG,n,t.
The outlier score calculating unit 3 outputs the calculated learning outlier score SG,n,t at each time t to each of the threshold calculating unit 5 and the outlier data extraction processing unit 7.
The threshold calculating unit 5 calculates the threshold Sth as illustrated in
The threshold calculating unit 5 outputs the calculated threshold Sth to the threshold storing unit 6.
The threshold storing unit 6 stores the threshold Sth output from the threshold calculating unit 5.
Hereinafter, an example of a process for calculating the threshold Sth by the threshold calculating unit 5 will be described.
First, the threshold calculating unit 5 calculates a mean value SG,ave of all the learning outlier scores SG,n,t calculated from the respective N pieces of learning time-series data DG,n,t by the outlier score calculating unit 3.
In addition, the threshold calculating unit 5 calculates a standard deviation a of all the learning outlier scores SG,n,t calculated from the respective N pieces of learning time-series data DG,n,t by the outlier score calculating unit 3.
Next, the threshold calculating unit 5 calculates the threshold Sth from the mean value SG,ave and the standard deviation σ as illustrated in the following formula (1).
S
th
=S
G,ave+3σ (1)
In the abnormality detection device illustrated in
However, this is only an example, and the threshold calculating unit 5 may separately calculate the threshold Sth used during learning and the threshold Sth used during abnormality detection.
As the threshold Sth used during learning, for example, a threshold in a range of (SG,ave+σ) to (SG,ave+2σ) is calculated as a threshold less than the threshold Sth illustrated in formula (1) in such a manner that the outlier data extraction processing unit 7 can extract many pieces of learning outlier data ODG,n,ts-te.
As the threshold Sth used during abnormality detection, for example, the threshold Sth illustrated in formula (1) is calculated.
The outlier data extraction processing unit 7 acquires the learning outlier score SG,n,t calculated by the outlier score calculating unit 3 at each time t and acquires the threshold Sth stored by the threshold storing unit 6.
The outlier data extraction processing unit 7 compares the learning outlier score SG,n,t at each time t with the threshold Sth.
The outlier data extraction processing unit 7 detects a period ts-te in which a learning outlier score SG,n,t is equal to or more than the threshold Sth by specifying a learning outlier score SG,n,t equal to or more than the threshold Sth among the learning outlier scores SG,n,t at respective times t on the basis of a comparison result between the learning outlier score SG,n,t and the threshold Sth.
The outlier data extraction processing unit 7 extracts learning time-series data DG,n,ts to DG,n,te in the detection period ts-te as learning outlier data ODG,n,ts-te from among pieces of learning time-series data DG,n,t (step ST4 in
The outlier data extraction processing unit 7 outputs the extracted learning outlier data ODG,n,ts-te to each of the type determining unit 9, the waveform condition selecting unit 10, the waveform classifying unit 13, and the waveform condition generation processing unit 14.
When receiving the learning outlier data ODG,n,ts-te from the outlier data extraction processing unit 7, the type determining unit 9 calculates a feature amount CG,n of the learning outlier data ODG,n,ts-te, and determines the waveform type of the learning outlier data ODG,n,ts-te from the feature amount CG,n (step ST5 in
The type determining unit 9 outputs the determined waveform type of the learning outlier data ODG,n,ts-te to the waveform classifying unit 13.
Hereinafter, a process for determining a waveform type by the type determining unit 9 will be specifically described.
Here, an example in which the type determining unit 9 classifies the waveforms of pieces of learning outlier data ODG,n,ts-te into six groups of an upper peak type waveform, a lower peak type waveform, an upper and lower peak type waveform, a transient ascending type waveform, a transient descending type waveform, and a vibration type waveform will be described.
In
[Upper Peak Type]
In the upper peak type waveform, as illustrated in
[Lower Peak Type]
In the lower peak type waveform, as illustrated in
[Upper and Lower Peak Type]
In the upper and lower peak type waveform, as illustrated in
In addition, in the upper and lower peak type waveforms, a value of the learning outlier data ODG,n,ts-te rises sharply to a maximum value, then falls sharply to a minimum value, and then returns to the vicinity of the value observed before the value of the learning outlier data ODG,n,ts-te rises sharply.
[Transient Ascending Type]
In the transient ascending type waveform, as illustrated in
[Transient Descending Type Waveform]
In the transient descending type waveform, as illustrated in
[Vibration Type Waveform]
In the vibration type waveform, as illustrated in
First, the type determining unit 9 calculates a mean value DG,n,ave of the pieces of learning outlier data ODG,n,ts-te output from the outlier data extraction processing unit 7.
The type determining unit 9 counts the number of intersections CN, which is the number of times the learning outlier data ODG,n,ts-te intersects with the mean value DG,n,ave, as one of the feature amounts CG, n.
The learning outlier data ODG,n,ts-te illustrated in
The type determining unit 9 focuses on the first intersection counting from the start point of the learning outlier data ODG,n,ts-te among one or more intersections where the learning outlier data ODG,n,ts-te intersects with the mean value DG,n,ave.
When the learning outlier data ODG,n,ts-te changes from a value lower than the mean value DG,n,ave to a value higher than the mean value DG,n,ave after the first intersection, the type determining unit 9 takes “first intersection=positive” as one of the feature amounts CG,n.
When the learning outlier data ODG,n,ts-te changes from a value higher than the mean value DG,n,ave to a value lower than the mean value DG,n,ave after the first intersection, the type determining unit 9 takes “first intersection=negative” as one of the feature amounts CG,n.
In the learning outlier data ODG,n,ts-te illustrated in
In addition, the type determining unit 9 calculates, as one of the feature amounts CG,n, an absolute value Δs-e of a difference between the start point of the learning outlier data ODG,n,ts-te and the end point of the learning outlier data ODG,n,ts-te.
Furthermore, the type determining unit 9 calculates, as one of the feature amounts CG,n, an absolute value Δmax-min of a difference between a maximum value out of pieces of learning outlier data ODG,n,ts-te and a minimum value out of pieces of learning outlier data ODG,n,ts-te.
When the number of intersections CN is 2 and “first intersection=positive”, the type determining unit 9 determines that the waveform type is “upper peak type”.
When the number of intersections CN is 1, “first intersection=positive”, and Δs-e≤Δmax-min×α, the type determining unit 9 determines that the waveform type is “upper peak type”. Provided that α is an arbitrary constant, and 0≤α≤1. The constant α may be stored in an internal memory of the type determining unit 9 or may be given from the outside.
When the number of intersections CN is 2 and “first intersection=negative”, the type determining unit 9 determines that the waveform type is “lower peak type”.
When the number of intersections CN is 1, “first intersection=negative”, and Δs-e≤Δmax-min×α, the type determining unit 9 determines that the waveform type is “lower peak type”.
When the number of intersections CN is 3 and Δs-e≤Δmax-min×β, the type determining unit 9 determines that the waveform type is “upper and lower peak type”. Provided that β is an arbitrary constant, and 0≤β≤1. The constant β may be stored in an internal memory of the type determining unit 9 or may be given from the outside.
When the number of intersections CN is 1, “first intersection=positive”, and Δs-e>Δmax-min×α, the type determining unit 9 determines that the waveform type is “transient ascending type”.
When the number of intersections CN is 1, “first intersection=negative”, and Δs-e>Δmax-min×α, the type determining unit 9 determines that the waveform type is “transient descending type”.
When the number of intersections CN is 4 or more, the type determining unit 9 determines that the waveform type is “vibration type”.
When the number of intersections CN is 3 and Δs-e>Δmax-min×α, the type determining unit 9 determines that the waveform type is “vibration type”.
The waveform classifying unit 13 classifies one or more pieces of learning outlier data ODG,n,ts-te whose waveforms have been determined to be of the same type by the type determining unit 9 out of the pieces of learning outlier data ODG,n,ts-te extracted by the outlier data extracting unit 4 into groups.
Next, the waveform classifying unit 13 calculates, for each of the provided groups, the degree of similarity between one or more pieces of learning outlier data ODG,n,ts-te included in the group.
As the degree of similarity between one or more pieces of learning outlier data ODG,n,ts-te, a distance between the waveforms of one or more pieces of learning outlier data ODG,n,ts-te may be calculated. As the distance to be calculated, a Euclidean distance, a 1-correlation coefficient, a Manhattan distance, a dynamic time warping (DTW) distance, and the like are conceivable. The shorter the distance, the higher the degree of similarity.
Since a process itself for calculating a distance between the waveforms of one or more pieces of learning outlier data ODG,n,ts-te is a known technique, detailed description thereof is omitted.
The waveform classifying unit 13 further classifies one or more pieces of learning outlier data ODG,n,ts-te classified into the same group into groups on the basis of the calculated degree of similarity (step ST6 in
Specifically, the waveform classifying unit 13 performs clustering of learning outlier data ODG,n,ts-te in such a manner that pieces of learning outlier data ODG,n,ts-te having the calculated high degree of similarity to each other are included in the same group among one or more pieces of learning outlier data ODG,n,ts-te classified into the same group. The waveform classifying unit 13 determines, for example, that pieces of learning outlier data ODG,n,ts-te having the calculated degree of similarity higher than or equal to a threshold are pieces of learning outlier data ODG,n,ts-te having a high degree of similarity to each other.
As a clustering method, a k-means method can be used. However, the clustering method is not limited to the k-means method, and spectral clustering, hierarchical clustering, or the like may be used.
The threshold to be compared with the calculated degree of similarity may be stored in an internal memory of the type determining unit 9 or may be given from the outside.
The waveform classifying unit 13 outputs a classification result of one or more pieces of learning outlier data ODG,n,ts-te to the waveform condition generation processing unit 14.
The waveform condition generation processing unit 14 generates, for each of the groups provided by the waveform classifying unit 13, a waveform condition Wp corresponding to the group from the waveforms of the one or more pieces of learning outlier data ODG,n,ts-te included in the group (step ST7 in
The waveform condition generation processing unit 14 generates, for example, a band model indicating a normal range of a waveform as the waveform condition Wp.
The waveform condition generation processing unit 14 outputs the generated waveform condition Wp to the waveform condition storing unit 15.
The waveform condition storing unit 15 stores the waveform condition Wp output from the waveform condition generation processing unit 14.
Hereinafter, a process for generating a band model by the waveform condition generation processing unit 14 will be specifically described.
Here, for convenience of explanation, it is assumed that one or more pieces of learning outlier data ODG,n,ts-te included in one group are represented by P1, P2, . . . , Pm. It is assumed that a value of Pi at time t is represented by Pi[t]. i=1, 2, . . . , m. The time t is any time in the period ts-te, and specifically, the time t is any time when the time ts is replaced with 1 (t=1, 2, . . . , (te-ts)).
The waveform condition generation processing unit 14 calculates a mean value Pmean[t] of m pieces of Pi[t] at time t as illustrated in the following formula (2), and calculates a standard deviation Pstd[t] of m pieces of Pi[t] at time t as illustrated in the following formula (3).
The waveform condition generation processing unit 14 calculates an upper limit value Bupper[t] of a normal range indicated by a band model by using the mean value Pmean[t], the standard deviation Pstd[t], and a constant λ (1≤λ) as illustrated in the following formula (4). The constant λ may be stored in an internal memory of the waveform condition generation processing unit 14 or may be given from the outside.
B
upper[t]=Pmean[t]+Pstd[t]×λ (4)
The waveform condition generation processing unit 14 calculates a lower limit value Blower[t] of a normal range indicated by a band model by using the mean value Pmean[t], the standard deviation Pstd[t], and a constant λ (1≤λ) as illustrated in the following formula (5).
B
lower[t]=Pmean[t]−Pstd[t]×λ (5)
Here, the waveform condition generation processing unit 14 calculates the upper limit value Bupper[t] and the lower limit value Blower[t] of a normal range indicated by a band model by using the mean value Pmean[t] and the standard deviation Pstd[t]. However, this is only an example, and the waveform condition generation processing unit 14 may calculate the upper limit value Bupper[t] and the lower limit value Blower[t] of a normal range indicated by a band model by using a maximum value Pmax[t] and a minimum value Pmin[t] out of m pieces of Pi[t] at time t.
The waveform condition generation processing unit 14 determines the maximum value Pmax[t] out of m pieces of Pi[t] at time t as illustrated in the following formula (6), and determines the minimum value Pmin[t] out of m pieces of m Pi[t] at time t as illustrated in the following formula (7).
P
max[t]=max(P1[t],P2[t], . . . ,Pm[t]) (6)
P
min[t]=min(P1[t],P2[t], . . . ,Pm[t]) (7)
The waveform condition generation processing unit 14 calculates the upper limit value Bupper[t] of a normal range indicated by a band model by using the maximum value Pmax[t], the minimum value Pmin[t], and a constant δ (1≤δ≤m) as illustrated in the following formula (8).
In formula (8), Pmax[t−δ/2: t+δ/2] is a maximum value Pmax[t] at each time t included in time (t−δ/2) to time (t+δ/2).
The waveform condition generation processing unit 14 calculates the lower limit value Blower[t] of a normal range indicated by a band model by using the maximum value Pmax[t], the minimum value Pmin[t], and a constant δ (1≤δ≤m) as illustrated in the following formula (9).
In formula (9), Pmin[t−δ/2: t+δ/2] is a minimum value Pmin[t] at each time t included in time (t−δ/2) to time (t+δ/2).
In
The solid line part indicates learning outlier data ODG,n,ts-te, and the broken line part indicates learning time-series data DG,n,t before and after the learning outlier data ODG,n,ts-te.
In
In the example of
Next, an operation during abnormality detection in the abnormality detection device will be described.
First, the abnormality detection data inputting unit 2 receives input of abnormality detection time-series data DU,t indicating states of equipment which is an abnormality detection target at a plurality of times tin time series (step ST11 in
The abnormality detection data inputting unit 2 outputs the received abnormality detection time-series data DU,t to each of the outlier score calculating unit 3 and the outlier data extraction processing unit 7.
When receiving the abnormality detection time-series data DU,t output from the abnormality detection data inputting unit 2, the outlier score calculating unit 3 calculates an abnormality detection outlier score SU,t at each time t from the abnormality detection time-series data DU,t (step ST12 in
A process for calculating the abnormality detection outlier score SU,t is similar to the process for calculating a learning outlier score SG,n,t.
The outlier score calculating unit 3 outputs the calculated abnormality detection outlier score SU,t at each time t to the outlier data extraction processing unit 7.
The outlier data extraction processing unit 7 acquires the abnormality detection outlier score SU,t calculated by the outlier score calculating unit 3 at each time t and acquires the threshold Sth stored by the threshold storing unit 6.
The outlier data extraction processing unit 7 compares the abnormality detection outlier score SU,t at each time t with the threshold Sth.
The outlier data extraction processing unit 7 detects a period ts′-te′ in which an abnormality detection outlier score SU,t is equal to or more than the threshold Sth by specifying an abnormality detection outlier score SU,t equal to or more than the threshold Sth among the abnormality detection outlier scores SU,t at respective times t on the basis of a comparison result between the abnormality detection outlier score SU,t and the threshold Sth.
The outlier data extraction processing unit 7 extracts abnormality detection time-series data DU,ts′ to DU,te′ in the detection period ts′-te′ as abnormality detection outlier data ODU,ts′-te′ from among pieces of abnormality detection time-series data DU,t (step ST13 in
The outlier data extraction processing unit 7 outputs the extracted abnormality detection outlier data ODU,ts′-te′ to each of the type determining unit 9, the waveform condition selecting unit 10, and the abnormality determination processing unit 11.
In the abnormality detection device illustrated in
When receiving the abnormality detection outlier data ODU,ts′-te′ from the outlier data extraction processing unit 7, the type determining unit 9 calculates a feature amount Cu of the abnormality detection outlier data ODU,ts′-te′.
A process for calculating the feature amount Cu in the abnormality detection outlier data ODU,ts′-te′ is similar to the process for calculating a feature amount CG,n in the learning outlier data ODG,n,ts-te.
The type determining unit 9 determines the waveform type of the abnormality detection outlier data ODU,ts′-te′ from the feature amount Cu of the abnormality detection outlier data ODU,ts′-te′ (step ST14 in
A process for determining the waveform type of the abnormality detection outlier data ODU,ts′-te′ is similar to the process for determining the waveform type of the learning outlier data ODG,n,ts-te.
The type determining unit 9 outputs the determined waveform type to the waveform condition selecting unit 10.
The waveform condition selecting unit 10 calculates the degree of similarity between the abnormality detection outlier data ODU,ts′-te′ output from the outlier data extraction processing unit 7 and each of N pieces of learning outlier data ODG,n,ts-te output from the outlier data extraction processing unit 7.
As the degree of similarity between the abnormality detection outlier data ODU,ts′-te′ and the learning outlier data ODG,n,ts-te, a distance between the waveform of abnormality detection outlier data ODU,ts′-te′ and the waveform of the learning outlier data ODG,n,ts-te may be calculated. As the distance to be calculated, a Euclidean distance, a 1-correlation coefficient, a Manhattan distance, a DTW distance, and the like are conceivable. Since a process itself for calculating the distance is a known technique, detailed description thereof is omitted.
The waveform condition selecting unit 10 searches for a piece of learning outlier data ODG,n,ts-te having the highest degree of similarity to the abnormality detection outlier data ODU,ts′-te′ among N pieces of learning outlier data ODG,n,ts-te. The waveform type of the piece of learning outlier data ODG,n,ts-te having the highest degree of similarity to the abnormality detection outlier data ODU,ts′-te′ is the same as the waveform type of the abnormality detection outlier data ODU,ts′-te′.
The waveform condition selecting unit 10 selects a waveform condition Wp corresponding to a group including the piece of learning outlier data ODG,n,ts-te that has been searched for from among waveform conditions Wp corresponding to the one or more groups stored by the waveform condition storing unit 15 (step ST15 in
The waveform condition selecting unit 10 outputs the selected waveform condition Wp to the abnormality determination processing unit 11.
The abnormality determination processing unit 11 collates the waveform condition Wp selected by the waveform condition selecting unit 10 with the waveform of the abnormality detection outlier data ODU,ts′-te′ extracted by the outlier data extraction processing unit 7.
The abnormality determination processing unit 11 determines whether or not the equipment is operating abnormally on the basis of a collation result between the waveform condition Wp and the waveform of abnormality detection outlier data ODU,ts′-te′ (step ST16 in
The abnormality determination processing unit 11 outputs a determination result indicating whether or not the equipment is operating abnormally to the detection result outputting unit 16.
The detection result outputting unit 16 displays the determination result output from the abnormality determination processing unit 11 on, for example, a display (not illustrated) (step ST17 in
Hereinafter, a process for determining abnormality of equipment by the abnormality determination processing unit 11 will be specifically described.
In
When the waveform of abnormality detection outlier data ODU,ts′-te′ is equal to or more than the lower limit value Blower[t] of the bandpass and equal to or less than the upper limit value Bupper[t] of the bandpass over the entire period ts′-te′, the abnormality determination processing unit 11 determines that the equipment is operating normally because the waveform is included in the normal range.
The waveform of abnormality detection outlier data ODU,ts′-te′ illustrated in
When the waveform of abnormality detection outlier data ODU,ts′-te′ is less than the lower limit value Blower[t] at any time tin the period ts′-te′, or more than the upper limit value Bupper[t] at any time t, the abnormality determination processing unit 11 determines that the equipment is operating abnormally because the waveform deviates from the normal range.
The waveform of abnormality detection outlier data ODU,ts′-te′ illustrated in
Here, when the waveform of abnormality detection outlier data ODU,ts′-te′ is equal to or more than the lower limit value Blower[t] and equal to or less than the upper limit value Bupper[t] over the entire period ts′-te′, the abnormality determination processing unit 11 determines that the equipment is operating normally. However, this is only an example. Even when the waveform of abnormality detection outlier data ODU,ts′-te′ deviates from the normal range indicated by the band model, the abnormality determination processing unit 11 may determine that the equipment is operating normally as long as the outlier is within an allowable range.
This will be specifically described as follows.
The abnormality determination processing unit 11 prepares a variable K having an initial value of 0.
When a value of abnormality detection outlier data ODU,ts′-te′ is more than the upper limit value Bupper[t] at each time tin the period ts′-te′, the abnormality determination processing unit 11 adds “1” to the variable K. Therefore, for example, when there are three times as time t at which a value of abnormality detection outlier data ODU,ts′-te′ is more than the upper limit value Bupper[t], the abnormality determination processing unit 11 adds “3” to the variable K.
When a value of abnormality detection outlier data ODU,ts′-te′ is less than the lower limit value Blower[t] at each time tin the period ts′-te′, the abnormality determination processing unit 11 adds “1” to the variable K. Therefore, for example, when there are two times as time t at which a value of abnormality detection outlier data ODU,ts′-te′ is less than the lower limit value Blower[t], the abnormality determination processing unit 11 adds “2” to the variable K.
As illustrated in the following formula (10), when a value obtained by multiplying the period ts′-te′ by a coefficient ζ (0≤ζ<1) is equal to or more than the variable K, the abnormality determination processing unit 11 determines that the equipment is operating normally.
K≤|ts′-te′|×ζ (10)
When the value obtained by multiplying the period ts′-te′ by the coefficient ζ is less than the variable K, the abnormality determination processing unit 11 determines that the equipment is operating abnormally.
Note that the constant ζ may be stored in an internal memory of the abnormality determination processing unit 11 or may be given from the outside. When ζ=0, the allowable range is zero.
Here, as an example in which even when the waveform of abnormality detection outlier data ODU,ts′-te′ deviates from the normal range indicated by the band model, the abnormality determination processing unit 11 determines that the equipment is operating normally as long as the outlier is within an allowable range, an example is described in which the abnormality determination processing unit 11 determines that the equipment is operating normally when formula (10) is satisfied.
However, this is only an example, and the following specific examples are also conceivable.
Even when the number of outliers of the waveform of abnormality detection outlier data ODU,ts′-te′ from the normal range indicated by the band model is small, the width of each outlier may be large.
Meanwhile, even when the number of outliers of the waveform of abnormality detection outlier data ODU,ts′-te′ from the normal range indicated by the band model is large, the width of each outlier may be small.
For example, it is conceivable that possibility that the equipment is operating normally is higher in a case where the number of outliers when the width of outlier is about 1% of the width of the band model is 2 to 3 times than in a case where the number of outliers when the width of outlier is about the same as the width of the band model is one time.
The abnormality determination processing unit 11 prepares a variable G having an initial value of 0.
The abnormality determination processing unit 11 subtracts the upper limit value Bupper[t] from a value of abnormality detection outlier data ODU,ts′-te′ at each time tin the period ts′-te′, and adds the value obtained by the subtraction to the variable G when the value obtained by the subtraction is positive.
The abnormality determination processing unit 11 subtracts a value of abnormality detection outlier data ODU,ts′-te′ from the lower limit value Blower[t] at each time tin the period ts′-te′, and adds the value obtained by the subtraction to the variable G when the value obtained by the subtraction is positive.
The abnormality determination processing unit 11 determines that the equipment is operating normally when the variable G is equal to or less than the threshold Gth, and determines that the equipment is operating abnormally when the variable G is more than the threshold Gth.
The threshold Gth may be stored in an internal memory of the abnormality determination processing unit 11 or may be given from the outside.
As the threshold Gth, such a threshold Gth as illustrated in the following formula (11) or (12) can be used.
Gth=(max(Bupper[t])−min(Blower[t]))×θ (11)
In formula (11), max(Bupper[t]) represents a maximum value out of the upper limit values Bupper[t] in the period ts′-te′, min(Blower[t]) represents a minimum value out of the lower limit values Blower[t] in the period ts′-te′, and θ represents a coefficient equal to or more than 0. The coefficient θ may be stored in an internal memory of the abnormality determination processing unit 11 or may be given from the outside.
In formula (12), h represents the number of times tin the period ts′-te′.
In the abnormality detection device illustrated in
However, this is only an example, and the outlier data extraction processing unit 7 may extract two or more pieces of abnormality detection outlier data ODU,ts′-te′ having different detection periods ts′-te′ from each other from among pieces of abnormality detection time-series data DU,t.
When the outlier data extraction processing unit 7 extracts two or more pieces of abnormality detection outlier data ODU,ts′-te′, the type determining unit 9, the waveform condition selecting unit 10, and the abnormality determination processing unit 11 perform the process described above for each of the pieces of abnormality detection outlier data ODU,ts′-te′.
In the first embodiment described above, the abnormality detection device is configured in such a manner that the abnormality determining unit 8 collates a waveform of the abnormality detection outlier data extracted by the outlier data extracting unit 4 with a waveform condition for determining that a waveform indicating a change in the abnormality detection outlier data is a waveform obtained when equipment is operating normally, and determines whether or not the equipment is operating abnormally on the basis of a collation result between the waveform condition and the waveform of the abnormality detection outlier data. Therefore, the abnormality detection device can avoid occurrence of erroneous determination indicating that an abnormality has occurred in the equipment without preparing event information in advance.
In addition to an event that can be predicted in advance, there is an event that is difficult to predict. Therefore, event information cannot be prepared in advance in some cases.
Meanwhile, in the abnormality detection device illustrated in
In the abnormality detection device illustrated in
However, the lengths of the waveforms of one or more pieces of learning outlier data ODG,n,ts-te are not necessarily the same, but may be different.
For example, when the lengths of the waveforms of two pieces of learning outlier data ODG,n,ts-te are different, the waveform classifying unit 13 first aligns the beginning of a waveform having a shorter length with the beginning of a waveform having a longer length, and calculate a distance between the waveform having a shorter length and the waveform having a longer length.
The waveform classifying unit 13 repeatedly calculates a distance between the waveform having a shorter length and the waveform having a longer length while sliding the waveform having a shorter length in parallel to the waveform having a longer length until the end of the waveform having a shorter length coincides with the end of the waveform having a longer length.
The waveform classifying unit 13 selects a minimum distance out of all the calculated distances, and determines the degree of similarity corresponding to the selected distance as the degree of similarity between a piece of learning outlier data ODG,n,ts-te having a longer waveform length and a piece of learning outlier data ODG,n,ts-te having a shorter waveform length. As the degree of similarity corresponding to the distance, for example, an integral multiple of a reciprocal of the distance is conceivable.
When classifying pieces of learning outlier data ODG,n,ts-te having the degree of similarity equal to or higher than the threshold into the same group, the waveform classifying unit 13 specifies a slide position at which the degree of similarity of a piece of learning outlier data ODG,n,ts-te having a shorter waveform length with respect to a piece of learning outlier data ODG,n,ts-te having the longest waveform length is maximum.
The waveform classifying unit 13 disposes the piece of learning outlier data ODG,n,ts-te having a shorter waveform length at the slide position specified with respect to the piece of learning outlier data ODG,n,ts-te having the longest waveform length.
By disposing the piece of learning outlier data ODG,n,ts-te having a shorter waveform length at the specified slide position, the beginning of the piece of learning outlier data ODG,n,ts-te having a shorter waveform length may be located closer to the end than the beginning of the piece of learning outlier data ODG,n,ts-te having the longest waveform length.
By adding a piece of learning time-series data DG,n,t at a time earlier than the piece of learning outlier data ODG,n,ts-te having a shorter waveform length to a beginning side of the piece of learning outlier data ODG,n,ts-te having a shorter waveform length, the waveform classifying unit 13 aligns the beginning of the piece of learning outlier data ODG,n,ts-te having a shorter waveform length with the beginning of the piece of learning outlier data ODG,n,ts-te having the longest waveform length.
In addition, by disposing the piece of learning outlier data ODG,n,ts-te having a shorter waveform length at the specified slide position, the end of the piece of learning outlier data ODG,n,ts-te having a shorter waveform length may be located closer to a beginning side than the end of the piece of learning outlier data ODG,n,ts-te having the longest waveform length.
By adding a piece of learning time-series data DG,n,t at a time later than the piece of learning outlier data ODG,n,ts-te having a shorter waveform length to an end side of the piece of learning outlier data ODG,n,ts-te having a shorter waveform length, the waveform classifying unit 13 aligns the end of the piece of learning outlier data ODG,n,ts-te having a shorter waveform length with the end of the piece of learning outlier data ODG,n,ts-te having the longest waveform length.
The waveform classifying unit 13 classifies the same pieces of learning outlier data ODG,n,ts-te having the same waveform length into the same group.
In the abnormality detection device illustrated in
An observed value of a sensor may be the outside air temperature or the seawater temperature, or the observed value of the sensor may be affected by external factors from other equipment. In these cases, since a waveform related to an event appears in a long-term trend of the learning outlier data ODG,n,ts-te, even when pieces of the learning outlier data ODG,n,ts-te have similar waveforms or change widths to each other, ranges of observed values may be different from each other.
When the ranges of observed values included in the pieces of learning outlier data ODG,n,ts-te are different from each other, the waveform classifying unit 13 may classify the pieces of learning outlier data ODG,n,ts-te into different groups because the pieces of learning outlier data ODG,n,ts-te are not similar to each other.
Therefore, the waveform classifying unit 13 calculates a mean value M of waveforms of each of the one or more pieces of learning outlier data ODG,n,ts-te whose waveforms have been determined to be of the same type by the type determining unit 9.
The waveform classifying unit 13 subtracts the mean value M of waveforms of each of the one or more pieces of learning outlier data ODG,n,ts-te from a value at each time t.
When the waveform classifying unit 13 subtracts the mean value M of waveforms of each of the one or more pieces of learning outlier data ODG,n,ts-te from a value at each time t, the ranges of observed values included in the one or more pieces of learning outlier data ODG,n,ts-te can be the same.
In addition, when a change width of the one or more pieces of learning outlier data ODG,n,ts-te is also affected by external factors, the waveform classifying unit 13 may divide a value of each of the one or more pieces of learning outlier data ODG,n,ts-te at each time t by a standard deviation of the pieces of learning outlier data ODG,n,ts-te.
By dividing a value of each of the one or more pieces of learning outlier data ODG,n,ts-te at each time t by the standard deviation, the influence of external factors can be reduced.
In addition, the one or more pieces of learning outlier data ODG,n,ts-te may fluctuate in a time direction. For example, in an event waveform that appears in temperature data, the speed of temperature rise is high and the speed of temperature fall is slow in summer. On the contrary, the speed of temperature rise is low, and the speed of temperature fall is high in winter.
When the one or more pieces of learning outlier data ODG,n,ts-te fluctuate in the time direction, the waveform classifying unit 13 calculates a DTW distance between the one or more pieces of learning outlier data ODG,n,ts-te by using a dynamic time warping method.
By expanding and contracting each of waveforms of the one or more pieces of learning outlier data ODG,n,ts-te according to an expansion and contraction path obtained by calculating the DTW distance, the waveform classifying unit 13 can eliminate the fluctuation of the learning outlier data ODG,n,ts-te in the time direction. The expansion and contraction path indicates time corresponding to one or more pieces of learning outlier data ODG,n,ts-te obtained when a distance between the one or more pieces of learning outlier data ODG,n,ts-te is a minimum. Since a process itself for expanding and contracting the waveform of learning outlier data ODG,n,ts-te according to an expansion and contraction path is a known technique, detailed description thereof is omitted.
In the abnormality detection device illustrated in
However, this is only an example, and instead of using a mean value Pmean[t] at time t, the waveform condition generation processing unit 14 may use an observed value at time t included in a representative piece of learning outlier data ODG,n,ts-te out of one or more pieces of learning outlier data ODG,n,ts-te included in a group.
As the representative piece of learning outlier data ODG,n,ts-te, a piece of learning outlier data ODG,n,ts-te having the highest degree of similarity to mean outlier data of one or more pieces of learning outlier data ODG,n,ts-te included in a group can be used.
In the abnormality detection device illustrated in
The waveform condition generation processing unit 14 may extend the normal range indicated by the band model by calculating a margin of the normal range from a width of the normal range, and adding the margin to the normal range.
This will be specifically described as follows.
As illustrated in the following formula (13), the waveform condition generation processing unit 14 calculates a margin r of the normal range from the width of the normal range indicated by the band model.
r=(max(Bupper[t])−min(Blower[t]))×η (13)
In formula (13), max(Bupper[t]) represents a maximum value out of upper limit values Bupper[t] in the period ts-te, min(Blower[t]) represents a minimum value out of lower limit values Blower[t] in the period ts-te, and η represents a coefficient equal to or more than 0. The coefficient η may be stored in an internal memory of the waveform condition generation processing unit 14 or may be given from the outside.
The waveform condition generation processing unit 14 extends the normal range by adding the margin r to the upper limit value Bupper[t] as illustrated in the following formula (14) and subtracting the margin r from the lower limit value Blower[t] as illustrated in the following formula (15).
B
upper[t]←Bupper[t]+r (14)
B
lower[t]←Blower[t]−r (15)
Here, the waveform condition generation processing unit 14 calculates the margin r of the normal range according to formula (13). However, this is only an example, and the waveform condition generation processing unit 14 may calculate the margin r of the normal range according to the following formula (16).
In formula (16), p represents the number of times tin the period ts-te.
In the abnormality detection device illustrated in
In the second embodiment, an abnormality detection device will be described in which the waveform condition generation processing unit 14 generates a histogram indicating a time period in which learning outlier data ODG,n,ts-te is generated when equipment is operating normally, as a waveform condition Wp.
The configuration of the abnormality detection device of the second embodiment is similar to the configuration of the abnormality detection device of the first embodiment, and the configuration diagram of the abnormality detection device of the second embodiment is illustrated in
The waveform condition generation processing unit 14 generates, for each of groups provided by the waveform classifying unit 13, a histogram indicating a time period in which one or more pieces of learning outlier data ODG,n,ts-te included in the group are generated, as a waveform condition Wp.
The learning outlier data ODG,n,ts-te includes period information indicating a period ts-te in which a learning outlier score SG,n,t is equal to or more than a threshold Sth. The period information includes information indicating a start time when the learning outlier score SG,n,t becomes equal to or more than the threshold Sth, and information indicating an end time when the learning outlier score SG,n,t becomes equal to or less than the threshold Sth.
The information indicating the start time and the information indicating the end time each include not only information indicating a so-called time but also information indicating a date and information indicating a day of the week.
Since a process itself for generating a histogram is a known technique, detailed description thereof is omitted, but a histogram can be generated on the basis of the period ts-te indicated by the period information included in the learning outlier data ODG,n,ts-te.
In
As in the first embodiment, the waveform condition selecting unit 10 calculates the degree of similarity between abnormality detection outlier data ODU,ts′-te′ output from the type determining unit 9 and each of N pieces of learning outlier data ODG,n,ts-te.
As in the first embodiment, the waveform condition selecting unit 10 searches for a piece of learning outlier data ODG,n,ts-te having the highest degree of similarity to the abnormality detection outlier data ODU,ts′-te′ among N pieces of learning outlier data ODG,n,ts-te.
As in the first embodiment, the waveform condition selecting unit 10 selects a waveform condition Wp corresponding to a group including the learning outlier data ODG,n,ts-te that has been searched for from among waveform conditions Wp corresponding to one or more groups stored by the waveform condition storing unit 15.
The waveform condition Wp selected by the waveform condition selecting unit 10 is a histogram generated by the waveform condition generation processing unit 14.
The waveform condition selecting unit 10 outputs the selected waveform condition Wp to the abnormality determination processing unit 11.
The abnormality determination processing unit 11 refers to period information included in the abnormality detection outlier data ODU,ts′-te′ output from the outlier data extraction processing unit 7, and recognizes a period ts′-te′ which is a time period in which abnormality detection outlier data ODU,ts′-te′ occurs.
The abnormality determination processing unit 11 collates the period ts′-te′ in which abnormality detection outlier data ODU,ts′-te′ is generated with a generation time period indicated by a histogram which is a waveform condition Wp output from the waveform condition selecting unit 10.
The abnormality determination processing unit 11 determines that equipment is operating normally when the period ts′-te′ in which the abnormality detection outlier data ODU,ts′-te′ is generated is included in the generation time period indicated by the histogram.
In the example of
The abnormality determination processing unit 11 determines that equipment is operating abnormally when the period ts′-te′ in which the abnormality detection outlier data ODU,ts′-te′ is generated is not included in the generation time period indicated by the histogram.
In the example of
In the second embodiment described above, the abnormality detection device is configured in such a manner that the abnormality determining unit 8 determines that equipment is operating normally when a time period in which abnormality detection outlier data extracted by the outlier data extracting unit 4 is generated is included in a generation time period indicated by a histogram, and determines that the equipment is operating abnormally when the time period in which the abnormality detection outlier data is generated is not included in the generation time period indicated by the histogram. Therefore, the abnormality detection device can avoid occurrence of erroneous determination indicating that an abnormality has occurred in the equipment without preparing event information in advance.
In the abnormality detection device of the second embodiment, the abnormality determination processing unit 11 determines that equipment is operating normally when a time period in which the abnormality detection outlier data ODU,ts′-te′ is generated is included in the generation time period indicated by the histogram.
In the abnormality detection device of the second embodiment, as in the abnormality detection device of the first embodiment, the abnormality determination processing unit 11 determines whether or not the waveform of the abnormality detection outlier data ODU,ts′-te′ is within a normal range of a bandpass over the entire period ts′-te′.
Then, the abnormality determination processing unit 11 may determine that equipment is operating normally when the time period in which the abnormality detection outlier data ODU,ts′-te′ is generated is included in the generation time period indicated by the histogram, and the waveform of the abnormality detection outlier data ODU,ts′-te′ is included in the normal range of the bandpass over the entire period ts′-te′.
In a third embodiment, an abnormality detection device including a selection accepting unit 17 for presenting waveform conditions Wp generated by the waveform condition generation processing unit 14 and accepting user's selection of an effective waveform condition Wp from among the presented waveform conditions Wp will be described.
In
The selection accepting unit 17 is achieved by, for example, a selection accepting circuit 34 illustrated in
The selection accepting unit 17 presents waveform conditions Wp generated by the waveform condition generation processing unit 14 and accepts user's selection of an effective waveform condition Wp from among the presented waveform conditions Wp.
The selection accepting unit 17 leaves only an effective waveform condition Wp whose selection has been accepted as a waveform condition Wp generated by the waveform condition generation processing unit 14, and discards a waveform condition Wp whose selection has not been accepted.
In
Here, for example, to each of the input interface circuit 21, the input interface circuit 22, the outlier score calculating circuit 23, the threshold calculating circuit 24, the outlier data extraction processing circuit 26, the type determining circuit 27, the waveform condition selecting circuit 28, the abnormality determination processing circuit 29, the waveform classifying circuit 30, the waveform condition generation processing circuit 31, the detection result outputting circuit 33, and the selection accepting circuit 34, a single circuit, a composite circuit, a programmed processor, a parallel programmed processor, ASIC, FPGA, or a combination thereof is applicable.
The constituent elements of the abnormality detection device are not limited to those achieved by dedicated hardware, and the abnormality detection device may be achieved by software, firmware, or a combination of software and firmware.
When the abnormality detection device is achieved by software, firmware, or the like, the threshold storing unit 6 and the waveform condition storing unit 15 are configured on a memory 41 of a computer. A program for causing the computer to execute a processing procedure performed in the learning data inputting unit 1, the abnormality detection data inputting unit 2, the outlier score calculating unit 3, the threshold calculating unit 5, the outlier data extraction processing unit 7, the type determining unit 9, the waveform condition selecting unit 10, the abnormality determination processing unit 11, the waveform classifying unit 13, the waveform condition generation processing unit 14, the detection result outputting unit 16, and the selection accepting unit 17 is stored in the memory 41 illustrated in
Next, an operation of the abnormality detection device illustrated in
Provided that constituent elements other than the selection accepting unit 17 among the constituent elements of the abnormality detection device illustrated in
As illustrated in
A user can evaluate appropriateness of each of the waveform conditions Wp by confirming the list confirmation screen.
The list confirmation screen illustrated in
The list confirmation screen illustrated in
The selection accepting unit 17 accepts user's selection of a waveform condition Wp whose check box has been checked by a user among the one or more waveform conditions Wp generated by the waveform condition generation processing unit 14, as an effective waveform condition Wp.
The selection accepting unit 17 causes the waveform condition storing unit 15 to store only an effective waveform condition Wp whose selection has been accepted as a waveform condition Wp generated by the waveform condition generation processing unit 14.
The selection accepting unit 17 discards a waveform condition Wp whose selection has not been accepted, and does not causes the waveform condition storing unit 15 to store the waveform condition Wp whose selection has not been accepted.
The selection accepting unit 17 has a function of displaying learning outlier data ODG,n,ts-te from which the waveform conditions Wp displayed on the list confirmation screen have been generated on a display (not illustrated).
When a user clicks on any waveform condition Wp among the one or more waveform conditions Wp displayed on the list confirmation screen, the selection accepting unit 17 displays one or more pieces of learning outlier data ODG,n,ts-te from which the waveform condition Wp has been generated on a display (not illustrated).
The list confirmation screen illustrated in
By confirming the list confirmation screen, a user can determine a piece of learning outlier data ODG,n,ts-te which is considered to be unnecessary for generating a waveform condition Wp out of the 12 pieces of learning outlier data ODG,n,ts-te.
The list confirmation screen illustrated in
In the example of
The selection accepting unit 17 accepts user's selection of a piece of learning outlier data ODG,n,ts-te whose check box is not unchecked out of the 12 pieces of learning outlier data ODG,n,ts-te.
The waveform condition generation processing unit 14 regenerates a waveform condition Wp from a piece of learning outlier data ODG,n,ts-te whose selection has been accepted by the selection accepting unit 17.
The list confirmation screen illustrated in
In the generation method selecting box, a generation method for calculating an upper limit value Bupper[t] and a lower limit value Blower[t] of a normal range indicated by a band model which is a waveform condition Wp can be selected by using a mean value Pmean[t] and a standard deviation Pstd[t].
In addition, in the generation method selecting box, the generation method for calculating an upper limit value Bupper[t] and a lower limit value Blower[t] of a normal range indicated by a band model can be selected by using a maximum value Pmax[t] and a minimum value Pmin[t].
Therefore, a user can select a method for generating a waveform condition Wp by operating the generation method selecting box.
The selection accepting unit 17 accepts user's selection of a method for generating a waveform condition Wp, the selection being caused by an operation of the generation method selecting box by a user.
The waveform condition generation processing unit 14 generates a waveform condition Wp from a piece of learning outlier data ODG,n,ts-te whose selection has been accepted by the selection accepting unit 17 on the basis of a generation method whose selection has been accepted by the selection accepting unit 17.
The list confirmation screen illustrated in
Therefore, a user can select a margin by operating the margin selecting box.
The selection accepting unit 17 accepts user's selection of a margin, the selection being caused by an operation of the margin selecting box by a user.
The waveform condition generation processing unit 14 extends the normal range by adding a margin whose selection has been accepted by the selection accepting unit 17 to the normal range.
The list confirmation screen illustrated in
When a user clicks the “reflect” button, the waveform condition generation processing unit 14 regenerates a waveform condition Wp from a piece of learning outlier data ODG,n,ts-te whose selection has been accepted by the selection accepting unit 17, and operates so as to display the regenerated waveform condition Wp on the list confirmation screen.
When the user clicks the “save” button, it is operated in such a manner that the waveform condition Wp regenerated by the selection accepting unit 17 is stored in the waveform condition storing unit 15.
When the user clicks the “add” button, it is operated in such a manner that a piece of learning outlier data ODG,n,ts-te included in a group different from the group of the pieces of learning outlier data ODG,n,ts-te displayed on the list confirmation screen illustrated in
In the list confirmation screen illustrated in
In the third embodiment described above, the abnormality detection device is configured in such a manner that the selection accepting unit 17 presents a waveform condition Wp generated by the waveform condition generation processing unit 14, accepts user's selection of an effective waveform condition Wp from among the presented waveform conditions Wp, leaves only the effective waveform condition Wp whose selection has been accepted as a waveform condition Wp generated by the waveform condition generation processing unit 14, and discards a waveform condition Wp whose selection has not been accepted. Therefore, the abnormality detection device can generate a waveform condition Wp reflecting determination of a user.
In a fourth embodiment, an abnormality detection device will be described in which the waveform condition generating unit 12 uses, as learning outlier data ODG,n,ts-te, a piece of abnormality detection outlier data ODU,ts′-te′ collated with a waveform condition Wp when the abnormality determining unit 8 determines that equipment is operating abnormally.
In
A type determining unit 18 is achieved by, for example, a type determining circuit 35 illustrated in
As in the type determining unit 9 illustrated in
As in the type determining unit 9 illustrated in
The type determining unit 18 acquires, as learning outlier data ODG,n,ts-te, a piece of abnormality detection outlier data ODU,ts′-te′ collated with a waveform condition Wp from a detection result outputting unit 19 when the abnormality determination processing unit 11 determines that equipment is operating abnormally.
The type determining unit 18 calculates a feature amount of the acquired abnormality detection outlier data ODU,ts′-te′, and determines the waveform type of the abnormality detection outlier data ODU,ts′-te′ from the calculated feature amount. The type determining unit 18 outputs the determined waveform type of the abnormality detection outlier data ODU,ts′-te′ to the waveform classifying unit 13.
The detection result outputting unit 19 is achieved by, for example, a detection result outputting circuit 36 illustrated in
As in the detection result outputting unit 16 illustrated in
The detection result outputting unit 19 displays a piece of abnormality detection outlier data ODU,ts′-te′ collated with a waveform condition Wp and abnormality detection time-series data DU,t on, for example, a display when the abnormality determination processing unit 11 determines that equipment is operating abnormally.
The detection result outputting unit 19 accepts user's selection of a piece of abnormality detection outlier data ODU,ts′-te′ used as learning outlier data ODG,n,ts-te among pieces of abnormality detection outlier data ODU,ts′-te′ collated with a waveform condition Wp when the abnormality determination processing unit 11 determines that equipment is operating abnormally.
The detection result outputting unit 19 outputs, as learning outlier data ODG,n,ts-te, the piece of abnormality detection outlier data ODU,ts′-te′ whose selection has been accepted to each of the type determining unit 18, the waveform classifying unit 13, and the waveform condition generation processing unit 14.
In
Here, for example, to each of the input interface circuit 21, the input interface circuit 22, the outlier score calculating circuit 23, the threshold calculating circuit 24, the outlier data extraction processing circuit 26, the type determining circuit 35, the waveform condition selecting circuit 28, the abnormality determination processing circuit 29, the waveform classifying circuit 30, the waveform condition generation processing circuit 31, and the detection result outputting circuit 36, a single circuit, a composite circuit, a programmed processor, a parallel programmed processor, ASIC, FPGA, or a combination thereof is applicable.
The constituent elements of the abnormality detection device are not limited to those achieved by dedicated hardware, and the abnormality detection device may be achieved by software, firmware, or a combination of software and firmware.
When the abnormality detection device is achieved by software, firmware, or the like, the threshold storing unit 6 and the waveform condition storing unit 15 are configured on a memory 41 of a computer. A program for causing the computer to execute a processing procedure performed in the learning data inputting unit 1, the abnormality detection data inputting unit 2, the outlier score calculating unit 3, the threshold calculating unit 5, the outlier data extraction processing unit 7, the type determining unit 18, the waveform condition selecting unit 10, the abnormality determination processing unit 11, the waveform classifying unit 13, the waveform condition generation processing unit 14, and the detection result outputting unit 19 is stored in the memory 41 illustrated in
Next, an operation of the abnormality detection device illustrated in
As in the first embodiment, the abnormality determination processing unit 11 collates a waveform condition Wp selected by the waveform condition selecting unit 10 with a waveform of the abnormality detection outlier data ODU,ts′-te′ extracted by the outlier data extraction processing unit 7.
As in the first embodiment, the abnormality determination processing unit 11 determines whether or not equipment is operating abnormally on the basis of a collation result between the waveform condition Wp and the waveform of abnormality detection outlier data ODU,ts′-te′.
As in the first embodiment, the abnormality determination processing unit 11 outputs a determination result indicating whether or not the equipment is operating abnormally to the detection result outputting unit 19.
When determining that the equipment is operating abnormally, the abnormality determination processing unit 11 outputs a piece of abnormality detection outlier data ODU,ts′-te′ collated with the waveform condition Wp to the detection result outputting unit 19.
The detection result outputting unit 19 displays the determination result output from the abnormality determination processing unit 11 on, for example, a display (not illustrated).
As illustrated in
In
In the enlarged diagrams, the solid line part indicates abnormality detection outlier data ODU,ts′-te′ and the broken line part indicates abnormality detection time-series data DU,t before and after the abnormality detection outlier data ODU,ts′-te′.
In
The data display screen illustrated in
In the example of
The detection result outputting unit 19 accepts, as learning outlier data ODG,n,ts-te, user's selection of the piece of abnormality detection outlier data ODU,ts′-te′ whose check box has been checked by a user.
The detection result outputting unit 19 outputs, as learning outlier data ODG,n,ts-te, the piece of abnormality detection outlier data ODU,ts′-te′ whose selection has been accepted to each of the type determining unit 18, the waveform classifying unit 13, and the waveform condition generation processing unit 14.
As in the type determining unit 9 illustrated in
As in the type determining unit 9 illustrated in
The type determining unit 18 acquires, as learning outlier data ODG,n,ts-te, the abnormality detection outlier data ODU,ts′-te′ output from the detection result outputting unit 19.
The type determining unit 18 calculates a feature amount of the acquired abnormality detection outlier data ODU,ts′-te′, and determines the waveform type of the abnormality detection outlier data ODU,ts′-te′ from the calculated feature amount.
A process for determining the waveform type of the abnormality detection outlier data ODU,ts′-te′ is similar to the process for determining the waveform type of the learning outlier data ODG,n,ts-te.
The type determining unit 18 outputs the determined waveform type of the abnormality detection outlier data ODU,ts′-te′ to the waveform classifying unit 13.
Operations of the waveform classifying unit 13 and the waveform condition generation processing unit 14 are similar to those of the first embodiment except that abnormality detection outlier data ODU,ts′-te′ output from the detection result outputting unit 19 is used as learning outlier data ODG,n,ts-te.
In the fourth embodiment described above, the abnormality detection device is configured in such a manner that when the abnormality determining unit 8 determines that equipment is operating abnormally, the type determining unit 18 calculates a feature amount of abnormality detection outlier data collated with a waveform condition, and determines the waveform type of the abnormality detection outlier data collated with the waveform condition from the feature amount, and then, the waveform condition generating unit 12 generates, from waveforms of one or more pieces of outlier data whose waveforms have been determined to be of the same type by the type determining unit 18 out of the pieces of learning outlier data extracted by the outlier data extracting unit 4 and the pieces of abnormality detection outlier data collated with waveform conditions, a waveform condition corresponding to the type. Therefore, the abnormality detection device can increase the number of pieces of learning outlier data and improve the accuracy of waveform conditions corresponding to the types thereof as compared with the abnormality detection device of the first embodiment.
Note that in the present invention, it is possible to freely combine the embodiments to each other, modify any constituent element in each of the embodiments, or omit any constituent element in each of the embodiments within the scope of the invention.
The present invention is suitable for an abnormality detection device and an abnormality detection method for determining whether or not equipment is operating abnormally.
This application is a Continuation of PCT International Application No. PCT/JP2018/044643 filed on Dec. 5, 2018, which is hereby expressly incorporated by reference into the present application.
Number | Date | Country | |
---|---|---|---|
Parent | PCT/JP2018/044643 | Dec 2018 | US |
Child | 17232565 | US |