CN1213429A - Method of detecting abnormality and abnormality detection system - Google Patents
Method of detecting abnormality and abnormality detection system Download PDFInfo
- Publication number
- CN1213429A CN1213429A CN97192893A CN97192893A CN1213429A CN 1213429 A CN1213429 A CN 1213429A CN 97192893 A CN97192893 A CN 97192893A CN 97192893 A CN97192893 A CN 97192893A CN 1213429 A CN1213429 A CN 1213429A
- Authority
- CN
- China
- Prior art keywords
- vectors
- mrow
- time
- vector
- series data
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 230000005856 abnormality Effects 0.000 title claims abstract description 45
- 238000001514 detection method Methods 0.000 title claims description 22
- 238000000034 method Methods 0.000 title claims description 22
- 239000013598 vector Substances 0.000 claims abstract description 112
- 230000002159 abnormal effect Effects 0.000 claims abstract description 16
- 238000004364 calculation method Methods 0.000 claims description 14
- 238000013519 translation Methods 0.000 claims description 13
- 238000005259 measurement Methods 0.000 claims description 3
- 230000003287 optical effect Effects 0.000 claims description 3
- 238000001228 spectrum Methods 0.000 description 14
- 238000007796 conventional method Methods 0.000 description 7
- 238000004519 manufacturing process Methods 0.000 description 4
- 238000005070 sampling Methods 0.000 description 4
- 230000000052 comparative effect Effects 0.000 description 3
- 238000003745 diagnosis Methods 0.000 description 3
- 238000012731 temporal analysis Methods 0.000 description 3
- 238000012360 testing method Methods 0.000 description 3
- 238000000700 time series analysis Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- 239000000203 mixture Substances 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 238000010183 spectrum analysis Methods 0.000 description 2
- 241000196324 Embryophyta Species 0.000 description 1
- 241001074710 Eucalyptus populnea Species 0.000 description 1
- 238000013528 artificial neural network Methods 0.000 description 1
- 230000001174 ascending effect Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 239000000178 monomer Substances 0.000 description 1
- 238000005312 nonlinear dynamic Methods 0.000 description 1
- 230000000737 periodic effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000012552 review Methods 0.000 description 1
- 230000005236 sound signal Effects 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 230000001131 transforming effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01D—MEASURING NOT SPECIALLY ADAPTED FOR A SPECIFIC VARIABLE; ARRANGEMENTS FOR MEASURING TWO OR MORE VARIABLES NOT COVERED IN A SINGLE OTHER SUBCLASS; TARIFF METERING APPARATUS; MEASURING OR TESTING NOT OTHERWISE PROVIDED FOR
- G01D1/00—Measuring arrangements giving results other than momentary value of variable, of general application
- G01D1/18—Measuring arrangements giving results other than momentary value of variable, of general application with arrangements for signalling that a predetermined value of an unspecified parameter has been exceeded
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01D—MEASURING NOT SPECIALLY ADAPTED FOR A SPECIFIC VARIABLE; ARRANGEMENTS FOR MEASURING TWO OR MORE VARIABLES NOT COVERED IN A SINGLE OTHER SUBCLASS; TARIFF METERING APPARATUS; MEASURING OR TESTING NOT OTHERWISE PROVIDED FOR
- G01D1/00—Measuring arrangements giving results other than momentary value of variable, of general application
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M99/00—Subject matter not provided for in other groups of this subclass
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Testing Or Calibration Of Command Recording Devices (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
- Testing And Monitoring For Control Systems (AREA)
Abstract
A physical quantity which is emitted from a sample to be detected and which changes with the passage of time, is measured at a predetermined time interval, and abnormality of the sample is detected based on the obtained time series data. A plurality of first vectors of a suitable dimension are formed from the time series data, the first vectors are advanced in parallel for a suitable period of time to form second vectors, a deviation is calculated between the second vectors and the first vectors that are formed, the calculated deviation is compared with a preset threshold value, and whether the sample is abnormal or not is determined based upon the result of comparison.
Description
Technical Field
The present invention relates to a method of detecting an abnormality of a manufacturing apparatus, a manufacturing plant, a manufactured product, or the like, and to a system for implementing the method.
Description of the Related Art
Sound or vibration generated by a mechanical device, concentration generated by a reaction in a chemical plant, current or voltage applied to an electronic device are measured at intervals, thereby detecting an abnormality of the mechanical device, an abnormality of reaction control, or an abnormality of product quality.
As an example of such an abnormality detection method, there is a method of detecting an abnormality by spectral analysis of a time-series signal measured at a predetermined period, which is proposed in chapters 12 to 14 of "Numerical method (Numerical receiver in C)" translated by danqing, okun, zotocten, and xiaolin (1993).
In the method, a time-series signal measured at a predetermined cycle is subjected to fast Fourier transform at appropriate frequency intervals in an appropriate frequency range to obtain a coefficient a satisfying the following formula (1)kAnd lkAnd find ak 2+bk 2. At this time, at the start frequency and the end frequency of the fourier transform, in order to prevent the abrupt rise and fall of f (x) in the formula (1), an appropriate data window, such as a Hanning (Hanning) window, is set. Multiplying the Hainin window by f (x) to reduce ak 2+bk 2The calculation error of (2). <math> <mrow> <mi>f</mi> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>=</mo> <munder> <mi>Σ</mi> <mi>k</mi> </munder> <mrow> <mo>(</mo> <msub> <mi>a</mi> <mi>k</mi> </msub> <mi>sin</mi> <msub> <mi>ω</mi> <mi>k</mi> </msub> <mi>t</mi> <mo>+</mo> <msub> <mi>b</mi> <mi>k</mi> </msub> <mi>cos</mi> <msub> <mi>ω</mi> <mi>k</mi> </msub> <mi>t</mi> <mo>)</mo> </mrow> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow> </math>
A to be calculatedk 2+bk 2Plotting creates a power spectrum in the coordinate region with the horizontal axis log ω and the vertical axis log (spectral density). The periodic fluctuation component is subtracted therefrom so that an appropriate straight line is obtained by the least square method. Then, the presence or absence of an abnormality is determined based on the result of comparing the level or slope (power spectrum index) of the straight line thus obtained with a predetermined threshold value.
However, in the above-described method of detecting an abnormality by spectrum analysis, it is difficult to accurately obtain a straight line suitable for a power spectrum, and erroneous judgment occurs due to measurement noise contained in a time-series signal. In addition, it is difficult to set an appropriate data window, so the error of the power spectrum is large.
The main object of the present invention is to provide a method for detecting abnormalities and a system for implementing the same, which can accurately detect various abnormalities by using time-series data and a prediction model thereof.
Summary of the invention
The present inventors constructed a prediction model based on a plurality of past time-series signals, given the prediction model an obtained time-series signal, predicted the trend of the time-series signal, and compared the prediction result with actual time-series data to make the assumption that any abnormality can be detected.
Thus, the inventors have attempted to build many predictive models, such as: (1) linear autoregressive (e.g. "Time Series Analysis: prediction and Control" (Time Series Analysis: forms and Control) ", authors G.E.P.Box and G.M.Jenkins, publisher Holden-Day Inc., 1976), (2) neural networks (e.g." Los Alamos National laboratory report (Los Alamos National laboratory) "LA-UR 87-2662 (1987), LA-UR 88-418 (1988), authors A.S.Lapedes and R.Farber), (3) radius basis function networks (e.g." physical abstracts (physical D) "Vol.35, Vol.356 (1989), M.dagli), and (4) monomer projection methods (e.g." journal of Nature 344, Vol.734 (1990), Ma.1990), all these models are imperfect.
After trial and error, attention was paid to the time series analysis algorithm, which was proposed by Wayland et al, and published on "Physical Review Letters" volume 70, et al, 580-582 (1993), authors R.Wayland, D.Bromley, D.Pickett and A.Passeante.
As an abnormality detection method of the present invention, that is, measuring physical quantities of an object to be detected that change with time at predetermined time intervals and detecting an abnormality of the object to be detected from obtained time-series data, the method of the present invention creates a plurality of first vectors having an appropriate dimension from the time-series data, creates second vectors by shifting (transforming) the first vectors by an appropriate time interval, calculates deviations of the created first and second vectors, compares the calculated values with a predetermined threshold value, and judges whether the object to be detected is abnormal or not based on the comparison result.
In this method, the method of calculating the deviation includes the steps of: selecting an appropriate vector from the first vectors, extracting a predetermined number of vectors from the first vectors in order of distance from small to large near the selected vector, creating translation vectors by translating the extracted vectors and the selected vector by an appropriate time interval, obtaining difference vectors of differences between the created translation vectors and the selected vector and the extracted vectors, and calculating variances of the obtained difference vectors.
The abnormality detection method of the present invention is summarized below. First, a physical quantity of an object to be detected, such as a motor rotating sound or a current or voltage applied to an electronic device, is measured at an appropriate cycle, thereby obtaining time series data { x (t) }t1 N(N is the number of data points). A combined embedding dimension (embedding dimension) D and a sampling time interval Δ t are set, and a vector x (t) represented by the following equation (2) is created from the time series data. It has been demonstrated that the creation of the vectors x (t) separately for the different embedding dimensions D means that the time series data are obtained separately for a plurality of physical quantities. Therefore, the abnormality can be accurately detected. The vector x (T) is translated by an appropriate time interval T Δ T to create a vector y (T) represented by the following equation (3).
X(t)={x(t),x(t-Δt),…,
x(t-(D-1)Δt)} (2)
Y(t)={x(t+TΔt),x(t+TΔt-Δt),…,
x(t+TΔt-(D-1)Δt} (3)
Selecting a suitable vector X (t (0)) from the N vectors X (t),the distances h between the selected vector X (t (0)) and all (N-1) vectors X (t) except X (t (0)) are calculated according to the following formula (4)xAnd extracting the K nearest vectors in Euclidean (Euclidean) distance sense (K < N): <math> <mrow> <mi>hx</mi> <mo>=</mo> <msqrt> <munderover> <mi>Σ</mi> <mrow> <mi>d</mi> <mo>=</mo> <mn>0</mn> </mrow> <mrow> <mi>D</mi> <mo>-</mo> <mn>1</mn> </mrow> </munderover> <mo>[</mo> <mi>X</mi> <mrow> <mo>(</mo> <mi>t</mi> <mrow> <mo>(</mo> <mn>0</mn> <mo>)</mo> </mrow> <mo>-</mo> <mi>dΔt</mi> <mo>)</mo> </mrow> <mo>-</mo> <mi>X</mi> <mrow> <mo>(</mo> <mi>t</mi> <mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo> </mrow> <mo>-</mo> <mi>dΔt</mi> <mo>)</mo> </mrow> <msup> <mo>]</mo> <mn>2</mn> </msup> <mo></mo> </msqrt> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>4</mn> <mo>)</mo> </mrow> </mrow> </math>
the vectors Y (t (K)) corresponding to the thus obtained (K +1) vectors X (t (K)) (K =0,1, …, K) are calculated, and then the difference vectors V (t (K)) are calculated according to equation (5). Next, when the translation error e (trans) which is the variance of the obtained (K +1) vectors V (t (K)) is obtained according to the following equation (6), the degree of clutter of the time series data is quantitatively determined as e (trans).
V(t(k))=X(t(k))-Y(t(k)) (5) <math> <mrow> <mi>E</mi> <mrow> <mo>(</mo> <mi>trans</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mn>1</mn> <mrow> <mi>K</mi> <mo>+</mo> <mn>1</mn> </mrow> </mfrac> <munderover> <mi>Σ</mi> <mrow> <mi>k</mi> <mo>=</mo> <mn>0</mn> </mrow> <mi>K</mi> </munderover> <mfrac> <msup> <mrow> <mo>|</mo> <mi>V</mi> <mrow> <mo>(</mo> <mi>t</mi> <mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo> </mrow> <mo>)</mo> </mrow> <mo>-</mo> <mo><</mo> <mi>V</mi> <mo>></mo> <mo>|</mo> </mrow> <mn>2</mn> </msup> <msup> <mrow> <mo>|</mo> <mo><</mo> <mi>V</mi> <mo>></mo> <mo>|</mo> </mrow> <mn>2</mn> </msup> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>6</mn> <mo>)</mo> </mrow> </mrow> </math> Wherein, <math> <mrow> <mo><</mo> <mi>V</mi> <mo>></mo> <mo>=</mo> <mfrac> <mn>1</mn> <mrow> <mi>K</mi> <mo>+</mo> <mn>1</mn> </mrow> </mfrac> <munderover> <mi>Σ</mi> <mrow> <mi>k</mi> <mo>=</mo> <mn>0</mn> </mrow> <mi>K</mi> </munderover> <mi>V</mi> <mrow> <mo>(</mo> <mi>t</mi> <mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo> </mrow> <mo>)</mo> </mrow> </mrow> </math>
through diligent research by the inventors, it was found that abnormalities can be easily identified. This is because the numerical amount of e (trans) increases when an abnormality occurs in which the noise ratio value increases as contained in the time system data; when an abnormality such as an increase in the power spectrum index of noise included in the time-series data occurs, the value of e (trans) becomes small, and the change in e (trans) is more significant than the power spectrum index. The present method and system have thus been invented which sets a threshold for determining the occurrence of the former anomaly and another threshold for determining the occurrence of the latter anomaly, and detects an anomaly if e (trans) exceeds the area between these two thresholds.
Brief description of the drawings
FIG. 1 is a block diagram showing the composition of an anomaly detection system according to the present invention;
FIG. 2 is a flow chart showing the detection of an abnormality by the computer of FIG. 1;
FIG. 3 is a flow chart showing the detection of an anomaly by the computer of FIG. 1;
FIG. 4 is a flow chart showing the detection of an anomaly by the computer of FIG. 1;
FIG. 5 is a flow chart showing the detection of an abnormality by the computer of FIG. 1;
FIGS. 6(A) and 6(B) are graphs showing time series data used as comparative tests;
FIG. 7 is a graph showing the result of calculating E (trans) from the two time series data of FIGS. 6(A) and 6(B) by the method of the present invention;
FIG. 8 is a graph showing the result of creating a power spectrum from the two time-series data of FIGS. 6(A) and 6(B) by a conventional method;
FIGS. 9(A) and 9(B) are graphs showing time series data used as comparative tests;
FIG. 10 is a graph showing the result of calculating E (trans) from the two time series data of FIGS. 9(A) and 9(B) by the method of the present invention;
fig. 11 is a graph showing the result of creating a power spectrum from the two time-series data of fig. 9(a) and 9(B) by a conventional method.
Description of the preferred embodiments
Embodiments of the present invention are described below with reference to the drawings.
Fig. 1 is a block diagram showing the composition of an abnormality detection system of the present invention. In fig. 1, 1 is a measuring device, such as a sonic sensor, an ultrasonic sensor, a vibration sensor, an optical sensor, a voltmeter, or an ammeter. The type of measuring device 1 is selected according to the physical quantity to be determined. The measuring apparatus 1 measures a physical quantity of an object to be measured, which is a target of abnormality detection, at a predetermined cycle. The output of the measuring device 1 is fed to an analog/digital (a/D) converter 2, the output being converted into a digital signal which is fed to a memory 32 provided to the computer 3 for storage therein as time series data.
The time-series data stored in the memory 32 is sent to the vector deviation calculating section 33 by the central processing unit 31 for a predetermined time. The vector deviation calculating unit 33 calculates e (trans) in the manner described later, and sends e (trans) to the failure diagnosing unit 34, which determines whether e (trans) is abnormal. The first and second thresholds are set at the failure condition input section 35, and the failure diagnosis section 34 compares e (trans) from the vector deviation calculation section 33 with the first and second thresholds set at the failure condition input section 35, respectively. If e (trans) is not less than the first threshold value or not more than the second threshold value, the failure diagnosing section 34 judges that e (trans) is abnormal, and causes the output section 36 to output an alarm signal to the external apparatus 4 such as a system controller or an alarm device.
Fig. 2 to 5 are flowcharts showing abnormality detection by the computer 3 in fig. 1. The vector deviation calculation unit 33 of the computer 3 sets T among the embedding dimension D, the sampling interval Δ T, and the translation time T Δ T, the vector selection number M, and the repetition numbers K and Q, and the failure condition input unit 35 sets the first and second threshold values (step S1). The memory 32 of the computer 3 stores therein the digital signal from the a/D converter 2 as time-series data (step S2), and the time-series data stored in the memory 32 is sent to the vector deviation calculating section 33 via the central processing unit 31 for a predetermined time.
According to known time series data { X (t) }t-1 N(N is the number of data points), the vector deviation calculation unit 33 creates a vector represented by the following expression (2) (step S3).
X(t)={x(t),x(t-Δt),…,
x(t-(D-1)Δt} (2)
Further, the vector deviation calculation section 33 translates the vector by an appropriate time interval T Δ T and creates a vector y (T) represented by the following expression (3) (step S4).
Y(t)={x(t+TΔt),x(t+TΔt-Δt),…,
x(t+TΔt-(D-1)Δt)} (3)
The vector deviation calculation section 33 generates an arbitrary random number (step S5), and selects M (M < N) vectors X (t (0)) from the N vectors X (t) created at S3 by using the generated random number (step S6), and then appropriately specifies one vector X (t (0)) (step S7). The vector deviation calculation unit 33 calculates a predetermined vector X (t (0)) and a division vector according to the following expression (4)Respective distances h between the remaining (N-1) vectors X (t) other than X (t (0)x(step S8), and rearranges these calculated distances in ascending order, and extracts K (K < M) nearest vectors from the head, thereby obtaining a vector X (t (K)) (K =1, …, K) (step S9). <math> <mrow> <mi>hx</mi> <mo>=</mo> <msqrt> <munderover> <mi>Σ</mi> <mrow> <mi>d</mi> <mo>=</mo> <mn>0</mn> </mrow> <mrow> <mi>D</mi> <mo>-</mo> <mn>1</mn> </mrow> </munderover> <mo>[</mo> <mi>X</mi> <mrow> <mo>(</mo> <mi>t</mi> <mrow> <mo>(</mo> <mn>0</mn> <mo>)</mo> </mrow> <mo>-</mo> <mi>dΔt</mi> <mo>)</mo> </mrow> <mo>-</mo> <mi>X</mi> <mrow> <mo>(</mo> <mi>t</mi> <mrow> <mo>(</mo> <mn>0</mn> <mo>)</mo> </mrow> <mo>-</mo> <mi>dΔt</mi> <mo>)</mo> </mrow> <msup> <mo>]</mo> <mn>2</mn> </msup> </msqrt> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>4</mn> <mo>)</mo> </mrow> </mrow> </math>
By shifting the vector X (T (K)) (K =0,1, …, K) obtained at S7 and S9 by the time interval T Δ T, the vector deviation calculation section 33 creates the vector Y (T (K)) thus obtained (step S10).
The vector deviation calculation section 33 calculates (k +1) difference vectors V (t (k)) according to the following expression (5) (step S11), and then obtains an original (primary) translational error e (trans)) from each difference vector V (t (k)) and the following expression (7) (step S12).
V(t(k))=X(t(k))-Y(t(k)) (5) <math> <mrow> <mi>E</mi> <msub> <mrow> <mo>(</mo> <mi>trans</mi> <mo>)</mo> </mrow> <mn>1</mn> </msub> <mo>=</mo> <mfrac> <mn>1</mn> <mrow> <mi>K</mi> <mo>+</mo> <mn>1</mn> </mrow> </mfrac> <munderover> <mi>Σ</mi> <mrow> <mi>k</mi> <mo>-</mo> <mn>0</mn> </mrow> <mi>K</mi> </munderover> <mfrac> <msup> <mrow> <mo>|</mo> <mi>V</mi> <mrow> <mo>(</mo> <mi>t</mi> <mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo> </mrow> <mo>)</mo> </mrow> <mo>-</mo> <mo><</mo> <mi>V</mi> <mo>></mo> <mo>|</mo> </mrow> <mn>2</mn> </msup> <msup> <mrow> <mo>|</mo> <mo><</mo> <mi>V</mi> <mo>></mo> <mo>|</mo> </mrow> <mn>2</mn> </msup> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>6</mn> <mo>)</mo> </mrow> </mrow> </math>
Wherein, <math> <mrow> <mo><</mo> <mi>V</mi> <mo>></mo> <mo>=</mo> <mfrac> <mn>1</mn> <mrow> <mi>K</mi> <mo>+</mo> <mn>1</mn> </mrow> </mfrac> <munderover> <mi>Σ</mi> <mrow> <mi>k</mi> <mo>=</mo> <mn>0</mn> </mrow> <mi>K</mi> </munderover> <mi>V</mi> <mrow> <mo>(</mo> <mi>t</mi> <mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo> </mrow> <mo>)</mo> </mrow> </mrow> </math>
the vector deviation calculation section 33 repeats steps S7 to S12 until it judges at step S13 that M original translational errors E (trans) have been obtained1And is the original translation error E (trans) for the M vectors X (t (0)) selected at step S61. When the vector deviation calculating part 33 judges that M original translation errors E (trans) are obtained1It calculates its median value to obtain the original mean translational error E (trans)11(step S14).
Next, the vector deviation calculation section 33 repeats steps S5 to S14 until it judges at step S15 that Q original average translational errors E (trans) are obtained11Then, when it is judged that Q original average translation errors E (trans) have been obtained11It calculates its average value to obtain e (trans) (step 16). As a result, e (trans) of the time series data can be accurately obtained.
When e (trans) is obtained, the vector deviation calculation section 33 sends it to the failure diagnosis section 34. The failure diagnosis unit 34 compares e (trans) from the vector deviation calculation unit 33 with the first threshold value and/or the second threshold value set in the failure condition input unit 35 (step 17). When e (trans) is not less than the first threshold value or not more than the second threshold value, the failure diagnosing section 34 judges e (trans) as abnormal (step 18), and outputs an alarm signal from the output section 36 to the external apparatus 4 (step 4).
Next, the results of comparison for detecting abnormality by two methods, one using the method of the present invention and the other using a conventional method, are described below.
Fig. 6 and 9 are graphs showing time series data used for comparative tests. The vertical axis represents signal strength and the horizontal axis represents time from the start of measurement. Fig. 6 is time-series data representing vibration of the motor rotating portion of the simulated fan: FIG. 6(A) shows time-series data when no abnormality occurs; fig. 6(B) shows time-series data at the time of occurrence of an abnormality. The signal intensity represented by the vertical axis is obtained by converting the vibration into an electric signal with a vibration sensor, but a reflected wave of the vibrating portion may also be detected with an optical sensor. Further, the vibration component can also be obtained by a voltage value obtained by measuring a voltage characteristic of the motor with a voltmeter or a current value obtained by measuring a current characteristic of the motor with a ammeter.
Fig. 9 is a graph showing time-series data of an analog audio signal: FIG. 9(A) shows time-series data when no abnormality occurs; fig. 9(B) shows time-series data when an abnormality occurs. The signal intensity represented by the vertical axis is obtained by converting sound into an electric signal with an acoustic wave sensor.
FIG. 7 is a graph showing the result of calculating E (trans) from the two time-series data of FIGS. 6(A) and 6(B) by the method of the present invention. ● denotes non-exception. And o indicates an abnormality. Here, the sampling interval Δ T is set to 1, T of the translation time T Δ T is set to 5, the vector selection number M is set to 300, the vector selection number K is set to 4, the repetition number Q is selected to 20, and the embedding dimension D is set to 5, 6, 7, 8, 9, and 10. Further, fig. 8 is a graph showing the result of creating a power spectrum from the two kinds of time-series data shown in fig. 6(a) and 6(B) by a conventional method. The solid line represents non-anomalies and the dashed line represents anomalies.
As is clear from fig. 7, in the method of the present invention, e (trans) in the abnormal case is clearly separated from e (trans) in the non-abnormal case, and by setting 0.05 as a threshold, it is possible to accurately determine whether an abnormality or non-abnormality is present. Meanwhile, as is clear from fig. 8, in the conventional method, many portions of the power spectrum in the abnormal case overlap with the power spectrum in the non-abnormal case, and thus a judgment error occurs.
Fig. 10 is a graph showing the result of calculating e (trans) from the two types of time-series data of fig. 9(a) and 9(B) by the present invention. ● indicates no anomaly and O indicates an anomaly. Here, the sampling interval Δ T is set to 1, T of the translation time T Δ T is set to 5, the vector selection number M is set to 300, the vector selection number K is set to 4, the repetition number Q is set to 20, and the embedding dimension D is set to 11, 12, 13, 14, and 15. Further, fig. 11 is a graph showing the result of creating a power spectrum from the two kinds of time-series data shown in fig. 9(a) and 9(B) by a conventional method. The solid line represents non-anomalies and the dashed line represents anomalies.
As is clear from fig. 10, in the method of the present invention, e (trans) in the abnormal case is clearly separated from e (trans) in the non-abnormal case, and by setting 0.90 as a threshold, it is possible to accurately determine whether an abnormality or non-abnormality is present. Meanwhile, as is clear from fig. 11, in the conventional method, many portions of the power spectrum in the abnormal case overlap with the power spectrum in the non-abnormal case, and thus a judgment error occurs. Industrial applicability
As described in detail above, in the present invention, even if there is a change in the physical quantity to be measured due to the nonlinear dynamics, various abnormalities can be accurately detected. Therefore, an excellent effect is produced such that the operation of the manufacturing plant is accurately managed, the parts of the manufacturing apparatus are accurately replaced, and the quality of the product is accurately managed.
Claims (8)
1. An abnormality detection method that measures a physical quantity of an object to be detected that changes with time at predetermined time intervals and detects an abnormality of the object to be detected based on measured time-series data, the method comprising the steps of:
creating a plurality of first vectors of appropriate dimensions from the time series data;
creating second vectors by translating said first vectors by an appropriate time interval;
calculating a deviation of said created first and second vectors;
comparing said calculated values to a predetermined threshold; and
and judging whether the object to be detected is abnormal or not according to the comparison result.
2. The anomaly detection method in accordance with claim 1, wherein said step of calculating said deviations comprises:
a first step of selecting an appropriate vector from said first vectors;
a second step of extracting a predetermined number of vectors from the first vectors in order of a distance from small to large approaching the selected vector;
a third step of creating translation vectors by translating said extracted vectors and said selected vector by an appropriate time interval;
a fourth step of obtaining difference vectors between said created translation vectors and said selected vector and said extracted vectors; and
and a fifth step of calculating the variance of the obtained difference vectors.
3. The abnormality detection method according to claim 2, further comprising the steps of:
repeating said first to fifth steps while changing said selected vector, obtaining a plurality of intermediate values of the means by repeatedly calculating a plurality of variances; and
obtaining a median or average of the calculated plurality of variances, and setting the average of the obtained plurality of median or average as a calculated value to be compared with the threshold.
4. The abnormality detection method according to claim 2, further comprising the steps of:
calculating a plurality of variances by repeating the first to fifth steps while permuting the selected vector;
obtaining a median or mean of the calculated plurality of variances; and
the obtained median or average value is set as a calculated value to be compared with the threshold value.
5. An abnormality detection system that measures a physical quantity of an object to be detected that changes with time at predetermined time intervals and detects an abnormality of the object to be detected based on measured time-series data, the system comprising:
means for creating a plurality of first vectors of appropriate dimensions from the time series data;
means for creating second vectors by translating said first vector by an appropriate time interval;
calculating means for calculating a deviation between the first and second vectors that have been created;
means for comparing said calculated values with predetermined thresholds; and
and judging whether the object to be detected is abnormal according to the comparison result.
6. The anomaly detection system according to claim 5, wherein said calculation means comprises:
means for selecting an appropriate vector from said first vectors;
means for extracting a predetermined number of vectors from said first vectors in order of distance from small to large near the selected vector;
means for creating translation vectors by translating said extracted vectors and said selected vector for an appropriate time interval;
means for obtaining difference vectors between said created translation vectors and said selected vector and said extracted vectors; and
means for calculating the variance of said obtained difference vectors.
7. An anomaly detection system for detecting anomalies in an object to be detected, the system comprising:
a measuring device that measures a physical quantity of an object to be detected that changes with time at predetermined time intervals;
means for storing time series data of measurement values obtained by the measuring means;
means for creating a plurality of first vectors of appropriate dimensions from said stored time series data;
means for creating second vectors by translating said first vectors by an appropriate time interval;
means for calculating a deviation between said created first and second vectors;
means for comparing said calculated value with a predetermined threshold value; and
and detecting the object to be detected as abnormal according to the comparison result.
8. The anomaly detection system according to claim 7, wherein said measuring means is a device selected from the group consisting of: acoustic wave sensors, ultrasonic sensors, vibration sensors, optical sensors, voltmeters, and ammeters.
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP29733/96 | 1996-02-16 | ||
JP2973396 | 1996-02-16 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN1213429A true CN1213429A (en) | 1999-04-07 |
Family
ID=12284315
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN97192893A Pending CN1213429A (en) | 1996-02-16 | 1997-02-13 | Method of detecting abnormality and abnormality detection system |
Country Status (6)
Country | Link |
---|---|
EP (1) | EP0881467A4 (en) |
KR (1) | KR19990082532A (en) |
CN (1) | CN1213429A (en) |
AU (1) | AU1671697A (en) |
TW (1) | TW327221B (en) |
WO (1) | WO1997030329A1 (en) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101135614B (en) * | 2006-08-29 | 2010-09-08 | 松下电工株式会社 | Anomaly monitoring device and method |
CN101995263A (en) * | 2009-08-07 | 2011-03-30 | 罗伯特.博世有限公司 | Method and apparatus for detecting oscillations |
CN106872168A (en) * | 2017-03-09 | 2017-06-20 | 武汉科技大学 | A kind of method for diagnosing faults for rolling bearing |
CN107703916A (en) * | 2017-09-30 | 2018-02-16 | 英特尔产品(成都)有限公司 | Electrical pressure regulator diagnostic system and method |
CN109861857A (en) * | 2019-01-28 | 2019-06-07 | 网联清算有限公司 | Fault detection method and device |
CN109919495A (en) * | 2019-03-12 | 2019-06-21 | 中国计量大学 | Photovoltaic plant monitoring data exception and degree judgment method based on threshold calculations |
CN111984827A (en) * | 2019-05-24 | 2020-11-24 | 上海东方富联科技有限公司 | Door sensor data anomaly detection method and system, storage medium and terminal |
CN112710446A (en) * | 2020-12-21 | 2021-04-27 | 北京和中普方新能源科技有限公司 | Judgment method and system for vibration test of electric vehicle battery system and storage medium |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6526405B1 (en) * | 1999-12-17 | 2003-02-25 | Microsoft Corporation | Determining similarity between event types in sequences |
US7152068B2 (en) * | 2001-12-21 | 2006-12-19 | Honeywell International Inc. | Method and apparatus for retrieving time series data related to an activity |
KR101029394B1 (en) | 2006-09-28 | 2011-04-15 | 미쓰비시덴키 가부시키가이샤 | Fault detector, fault detection method, and fault detection program |
TWI474023B (en) * | 2008-12-10 | 2015-02-21 | Ind Tech Res Inst | Diagnosis method for diagnosing fault in operation of a motor fault and diagnosis device using the same |
KR101659989B1 (en) * | 2014-12-10 | 2016-09-27 | 부산대학교 산학협력단 | Apparatus and method for analyzing abnormal data using combination of multi-dimensional features |
CN105938011B (en) * | 2016-07-15 | 2017-08-11 | 武汉理工大学 | A kind of hydraulic decoking process drilling rod working condition real-time judgment method |
KR102230463B1 (en) * | 2019-12-12 | 2021-03-19 | 울산대학교 산학협력단 | Diagnosis system and method of defect of equipment component |
CN114184274B (en) * | 2021-11-19 | 2023-06-16 | 哈尔滨工程大学 | Nonlinear characteristic diagnosis method for vibration signal interfered by high-frequency noise |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4453407A (en) * | 1980-04-17 | 1984-06-12 | Hitachi, Ltd. | Vibration diagnosis method and apparatus for rotary machines |
JPH02226024A (en) * | 1989-02-27 | 1990-09-07 | Toshiba Corp | Rubbing detector of rotary machine |
JP2909134B2 (en) * | 1990-04-02 | 1999-06-23 | 株式会社東芝 | Abnormality detector for rotating parts of rotating machinery |
-
1997
- 1997-02-13 EP EP97902662A patent/EP0881467A4/en not_active Withdrawn
- 1997-02-13 AU AU16716/97A patent/AU1671697A/en not_active Abandoned
- 1997-02-13 CN CN97192893A patent/CN1213429A/en active Pending
- 1997-02-13 KR KR1019980706268A patent/KR19990082532A/en active IP Right Grant
- 1997-02-13 WO PCT/JP1997/000367 patent/WO1997030329A1/en not_active Application Discontinuation
- 1997-02-15 TW TW086101800A patent/TW327221B/en active
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101135614B (en) * | 2006-08-29 | 2010-09-08 | 松下电工株式会社 | Anomaly monitoring device and method |
CN101995263A (en) * | 2009-08-07 | 2011-03-30 | 罗伯特.博世有限公司 | Method and apparatus for detecting oscillations |
CN101995263B (en) * | 2009-08-07 | 2015-01-28 | 罗伯特.博世有限公司 | Method and apparatus for detecting oscillations |
CN106872168A (en) * | 2017-03-09 | 2017-06-20 | 武汉科技大学 | A kind of method for diagnosing faults for rolling bearing |
CN107703916A (en) * | 2017-09-30 | 2018-02-16 | 英特尔产品(成都)有限公司 | Electrical pressure regulator diagnostic system and method |
CN109861857A (en) * | 2019-01-28 | 2019-06-07 | 网联清算有限公司 | Fault detection method and device |
CN109919495A (en) * | 2019-03-12 | 2019-06-21 | 中国计量大学 | Photovoltaic plant monitoring data exception and degree judgment method based on threshold calculations |
CN111984827A (en) * | 2019-05-24 | 2020-11-24 | 上海东方富联科技有限公司 | Door sensor data anomaly detection method and system, storage medium and terminal |
CN112710446A (en) * | 2020-12-21 | 2021-04-27 | 北京和中普方新能源科技有限公司 | Judgment method and system for vibration test of electric vehicle battery system and storage medium |
Also Published As
Publication number | Publication date |
---|---|
EP0881467A1 (en) | 1998-12-02 |
KR19990082532A (en) | 1999-11-25 |
AU1671697A (en) | 1997-09-02 |
EP0881467A4 (en) | 1999-05-12 |
TW327221B (en) | 1998-02-21 |
WO1997030329A1 (en) | 1997-08-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN1213429A (en) | Method of detecting abnormality and abnormality detection system | |
Sohn et al. | Structural health monitoring using statistical pattern recognition techniques | |
US11934184B2 (en) | Automated analysis of non-stationary machine performance | |
RU2488815C2 (en) | Method and apparatus for classifying sound-generating processes | |
CN1271594C (en) | Pitch determination method and apparatus on spectral analysis | |
CN1658492A (en) | Encoder and control apparatus for motor | |
Flores-Fuentes et al. | Energy center detection in light scanning sensors for structural health monitoring accuracy enhancement | |
CN111680665B (en) | Motor mechanical fault diagnosis method adopting current signals based on data driving | |
JP2721799B2 (en) | Machine abnormality judgment method | |
Ramasso et al. | Clustering acoustic emission data streams with sequentially appearing clusters using mixture models | |
CN115959549A (en) | Escalator fault diagnosis method based on digital twinning | |
WO2020071066A1 (en) | Abnormality determination device, signal feature value predictor, abnormality determination method, learning model generation method, and learning model | |
Chen et al. | A motor fault diagnosis system based on cerebellar model articulation controller | |
CN117538710B (en) | Intelligent early warning method and system for local dynamic discharge monitoring | |
CN112131781B (en) | Steel structure damage detection method based on fully-connected neural network and transfer rate function | |
CN114862283B (en) | Quality inspection method and quality inspection device for machine tool machining | |
CN112860658A (en) | High-speed train electromagnetic emission prediction method based on data mining | |
CN112880726A (en) | Sensor fault diagnosis method and device based on variational modal decomposition sample entropy | |
CN112881839B (en) | Transformer diagnosis method based on mutual information of frequency concentration and vibration stability | |
CN114034375B (en) | Ultra-high voltage transmission line noise measurement system and method | |
CN113919388B (en) | Electromechanical equipment fault diagnosis method and device integrating signal spectrum amplitude modulation and deep learning | |
JP2667327B2 (en) | OF cable insulation deterioration diagnosis method and apparatus | |
RU2774697C1 (en) | Method for vibration diagnostics for determining the degree and speed of development of nearing defects in industrial equipment | |
Yang et al. | A CUSUM-based condition monitoring algorithm for power electronics based on SOM-MQE feature extraction method | |
CN112001326B (en) | Vibration fingerprint-based high-voltage switch equipment fault diagnosis and prediction method and system |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C02 | Deemed withdrawal of patent application after publication (patent law 2001) | ||
WD01 | Invention patent application deemed withdrawn after publication |