CN101738998A - System and method for monitoring industrial process based on local discriminatory analysis - Google Patents

System and method for monitoring industrial process based on local discriminatory analysis Download PDF

Info

Publication number
CN101738998A
CN101738998A CN200910155009A CN200910155009A CN101738998A CN 101738998 A CN101738998 A CN 101738998A CN 200910155009 A CN200910155009 A CN 200910155009A CN 200910155009 A CN200910155009 A CN 200910155009A CN 101738998 A CN101738998 A CN 101738998A
Authority
CN
China
Prior art keywords
data
training data
matrix
module
local
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.)
Granted
Application number
CN200910155009A
Other languages
Chinese (zh)
Other versions
CN101738998B (en
Inventor
荣冈
邵纪东
冯毅萍
吴玉成
许华
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Zhejiang University ZJU
Original Assignee
Zhejiang University ZJU
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Zhejiang University ZJU filed Critical Zhejiang University ZJU
Priority to CN2009101550099A priority Critical patent/CN101738998B/en
Publication of CN101738998A publication Critical patent/CN101738998A/en
Application granted granted Critical
Publication of CN101738998B publication Critical patent/CN101738998B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02PCLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
    • Y02P90/00Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
    • Y02P90/02Total factory control, e.g. smart factories, flexible manufacturing systems [FMS] or integrated manufacturing systems [IMS]

Abstract

The invention discloses a system and a method for monitoring an industrial process based on local discriminatory analysis. The monitoring system comprises a measuring instrument, a distributed control system, a server and an upper computer, wherein the server comprises a real-time database and a relation database; the upper computer comprises an off-line modeling module and an on-line modeling module. The off-line modeling module comprises the following modules: a) a training data extracting module, b) a first standardization module, and c) a model training module. The on-line monitoring module comprises the following modules: a) a signal acquisition module, b) a second standardization module, c) a fault detection and diagnosis module, and d) a result display module. The invention realizes the detection and the diagnosis of a fault at the same time, and has an obviously better effect than that of the system and the method for monitoring the industrial process based on the local discriminatory analysis of the conventional system.

Description

A kind of Industrial Process Monitoring system and method based on local discriminatory analysis
Technical field
The present invention relates to the Industrial Process Monitoring field, relate in particular to a kind of Industrial Process Monitoring system and method based on local discriminatory analysis.
Technical background
The normal even running of keeping complex industrial process is the common objective that each industrial trade such as oil, chemical industry, pharmacy, food is pursued.Its meaning not only is to guarantee the safe and reliable operation produced, also is to realize strict production assurance and controls environment to pollute etc.Effectively the real-time process monitoring system is the key that guarantees large complicated industrial process even running.
The widespread use in industrial process along with Distributed Control System (DCS) and various smart instrumentation, computer resource with low cost and reliable popularizing of Storage Techniques, process and qualitative datas a large amount of in the modern industry process are measured in real time and are noted.These data have accurately been described status of processes, for process monitoring provides reliable foundation.Use various data analysing methods that the representative process historical data of gathering under normal operating condition and under the fault condition is carried out modeling respectively based on the method for data-driven, analyze the process data realization monitoring of on-line measurement then according to institute's established model.These class methods only depend on measurement data, are specially adapted to be difficult to obtain accurately and the monitoring of the large-scale complex process of complete mechanism model.Process monitoring comprises the detection and the diagnosis of fault.The former judges whether current system exists fault, and the latter judges which classification is the fault taken place belong to, for the reparation of fault provides foundation.Existing about based on the patent of the process monitoring system of these class methods as:
Chinese invention patent 200610154826.9 discloses a kind of industrial process nonlinear fault diagnosis system based on FISHER, comprise the field intelligent instrument, DCS system and the host computer that are connected with industrial process object, described DCS system is made of data-interface, control station, database; Intelligence instrument, DCS system, host computer link to each other successively, and described host computer comprises standardization module, FISHER discriminatory analysis module and fault diagnosis module.Can access good fault diagnosis effect.
Chinese invention patent 200610154825.4 discloses a kind of industrial processes fault detection method based on wavelet analysis, comprise the field intelligent instrument, DCS system and the host computer that are connected with industrial process object, described DCS system is made of data-interface, control station, database; Intelligence instrument, DCS system, host computer link to each other successively, and described host computer comprises standardization module, wavelet decomposition module, pivot analysis functional module, wavelet reconstruction functional module, support vector machine classifier functional module and fault judgement module.Can access good diagnosis effect.
But there are two subject matters in current process monitoring method based on data-driven:
A) fault and diagnosis as two tasks independently.Promptly carry out fault detect earlier, after decision process is in malfunction, the classification of failure judgement again.When fault detect, only utilize gathered data under the process normal condition like this, do not utilized fault data.Therefore the effect of fault detect also has the space of further improving.
B) fault diagnosis uses the Fisher discriminatory analysis that fault data is carried out modeling.The Fisher discriminatory analysis only has only two classes in the data of analyzing, and every class data all meet under the situation of Gaussian distribution of identical covariance matrix for optimum.And the classification of procedure fault and corresponding fault data is often more than two classes, and also Gaussian distributed not necessarily of every class data.
Therefore existing method for diagnosing faults can not provide optimum accuracy rate of diagnosis under many circumstances.
Summary of the invention
In order to improve the deficiency that limitation is big, effect is general that existing process monitoring system uses, the invention provides a kind of applied widely, realize the detection and the diagnosis of fault simultaneously, and effect obviously is better than the Industrial Process Monitoring system based on local discriminatory analysis of existing system.
A kind of Industrial Process Monitoring system based on local discriminatory analysis, the host computer that comprises the measuring instrument, the Distributed Control System (DCS) that are connected with industrial process object, is used for storing the server of Distributed Control System (DCS) institute image data and is used for the image data that processing server stores.
Wherein measuring instrument is used to gather the real time data of industrial process object, and Distributed Control System (DCS) (DCS) is controlled industrial process object according to the real time data of measuring instrument collection.
Described server comprises the real-time data base of the real time data that is used for storage industry process object and is used to store the relational database of the data of described industrial process object under normal condition and all kinds of malfunction; Data in the described relational database, can be described as historical data, industrial process object is in normal condition or is in malfunction in the historical data, and the type of malfunction all is clear and definite, can think that each data point all has corresponding class sign.
Described host computer comprises off-line modeling module and on-line monitoring module;
Wherein said off-line modeling module comprises:
A) training data extraction module is used for extracting the class sign of data under normal condition and all kinds of malfunction and each data point correspondence as training data from relational database;
Before extracting, can preestablish the data variable that needs extraction, initial sum termination time, the sample number of sampling interval and each data class.Extraction obtains data set R DThe real number space of expression D dimension and corresponding class sign l 1... l n∈ 0,1,2...c}.The data type of available 0 mark normal condition, 1, the data type under the different fault of 2...c mark.
Each data point that expression is extracted from relational database;
B) first standardized module is used for training data is done standardization, obtains the training data x after the standardization i(i=1 ... n, n are that training data is counted); The step of standardization is as follows:
1) computation of mean values: Wherein Be the data point of extracting from relational database, n is the data point number, and i is the data point call number;
2) calculate variance vectors:
3) translation is flexible: Wherein ./be that vectorial corresponding element is divided by σ xBe the standard deviation vector, The average that obtains for step 1).
The training data that extracts is done standardization, and can make each variable (data point is made up of the value of a plurality of variablees) average is zero, and variance is 1.
C) model training module is used for the training data x after the standardization iCarry out the partial structurtes modeling, find the solution again and obtain optimum projection matrix A, with this optimum projection matrix training data is projected to the subspace and obtain the picture y of training data in the subspace i(i=1,2...n are natural number);
Described on-line monitoring module comprises:
A) signal acquisition module; Be used for gathering the real time data of industrial process object from real-time data base;
B) second standardized module; Be used to utilize the average of training data of first standardized module and variance to collecting to such an extent that real time data does that translation is flexible to be handled, obtain the real time data x of translation after flexible; Still use formula when translation is flexible But wherein Replace with real time data, and σ xWith The still average and the variance of the training data that obtains with first standardized module.
C) fault detection and diagnosis module; Real time data x after flexible projects in the subspace the optimum projection matrix A that obtains with the model training module with translation, obtains the picture y of real time data in the subspace, (y=A TX), use the nearest neighbor search that defines based on Euclidean distance to seek the picture y of training data in the subspace i(i=1 ... n) from picture y nearest some y p, according to y pStatus categories judge the status categories of real time data;
D) display module as a result; Be used to show the status categories of fault detection and diagnosis module judgement.
Wherein said model training module comprises:
The partial structurtes MBM; With the training data x of dispersion matrix between local within class scatter matrix drawn game category after to standardization iLocal geometry and local differentiate structure modeling respectively;
Module is found the solution in the optimum decision projection; Be used for local within class scatter matrix R according to training data w, dispersion matrix R between local class bWith local population variance degree matrix R tTry to achieve between the local class of data after the projection maximum and optimum projection A dispersion minimum in the local class of dispersion;
The training data projection module; With optimum projection matrix A with training data project to obtain in the subspace training data in the subspace the picture y i
Described partial structurtes MBM is with the training data x of dispersion matrix between local within class scatter matrix drawn game category after to standardization iLocal geometry and local differentiate structure modeling respectively, step is as follows:
1) each training data point x under the definition of calculating Euclidean distance iK-arest neighbors set of data points With M (x i) separated into two parts M w(x i) and M b(x i), M wherein w(x i) be and x iBelong to of a sort point, M b(x i) be and x iBelong to inhomogeneous point.
The class here is meant x iBe normal condition or certain malfunction, can identify by class and discern.K is the size of regional area, promptly with the number of each training data point arest neighbors data point.
2) calculate adjacency matrix W in the local class wAdjacency matrix W between the drawn game category b
W in the formula wAnd W bSubscript i, the j respectively line number of representing matrix and the index of columns.
3) calculate local within class scatter matrix R w, dispersion matrix R between local class bWith local population variance degree matrix R t
R w = 1 2 Σ i = 1 n Σ j = 1 n W w ; i , j ( x i - x j ) ( x i - x j ) T ,
R b = 1 2 Σ i = 1 n Σ j = 1 n W b ; i , j ( x i - x j ) ( x i - x j ) T ,
R t=R w+R b
Wherein subscript T represents transposed matrix; N represents training data point number.
Module is found the solution in described optimum decision projection, according to the local within class scatter matrix R of training data w, dispersion matrix R between local class bWith local population variance degree matrix R tAsk between the local class of the data that make after the projection maximum and optimum projection matrix A dispersion minimum in the local class of dispersion, step is as follows:
1) to R tImplement characteristic value decomposition R t=Q Λ Q T, Q=[q wherein 1... q m] form by nonzero eigenvalue characteristic of correspondence vector; Wherein subscript T represents that transposed matrix, Λ represent by R tThe diagonal matrix formed of nonzero eigenvalue;
2) compute matrix Λ -1Q TR bD proper vector b of Q eigenwert maximum 1... b d, be worth descending arrangement by characteristic of correspondence, the dimensionality reduction subspace dimension of d for setting;
3) calculate optimum projection matrix A=QB ∈ R D * d, B=[b wherein 1... b d];
R D * dBe the real number matrix of the capable d row of D, D is the dimension of training data.
Described training data projection module with optimum projection matrix with training data x i(i=1 ... n) project to obtain in the subspace training data in the subspace the picture y i, y i=A Tx i∈ R d(i=1 ... n).
The present invention also provides a kind of Industrial Process Monitoring method based on local discriminatory analysis, comprises the steps:
1) data and the corresponding class of extracting under industrial process object normal condition and all kinds of malfunction identifies as training data; Training data is done standardization, obtain the training data data x after the standardization i(i=1 ... n); To the training data x after the standardization iCarry out the partial structurtes modeling, find the solution and obtain optimum projection matrix, with this optimum projection matrix training data is projected to the subspace and obtain the picture data y of training data in the subspace i(i=1 ... n);
2) gather the industrial process object real time data; Utilize the average of training data in the step (1) and variance to collecting to such an extent that real time data does that translation is flexible to be handled, obtain the real time data y of translation after flexible i(i=1 ... n); Utilize the real time data x after the optimum projection matrix that obtains in the step 1) stretches translation to project in the subspace, obtain the picture y of real time data in the subspace, y=A TX uses based on the nearest neighbor search of Euclidean distance definition and seeks the picture y of training data in the subspace i(i=1 ... n) from picture y nearest some y p, according to y pStatus categories judge the status categories of real time data;
3) status categories to real time data shows.
Step (1) is described, and that training data is made the step of standardization is as follows:
1) computation of mean values: Wherein Be the normal condition of industrial process object or the data point under certain class malfunction, n is the number of data point;
2) calculate variance vectors:
3) translation is flexible: Wherein ./be that vectorial corresponding element is divided by σ xBe the standard deviation vector, The average that obtains for step 1).
In the step (1) to the training data x after the standardization iCarry out the partial structurtes modeling, find the solution and obtain optimum projection matrix, with this optimum projection matrix training data is projected to the subspace and obtain the picture y of training data in the subspace iStep comprise:
Calculate the k-arest neighbors set of data points of each the training data point under the Euclidean distance definition With M (x i) separated into two parts M w(x i) and M b(x i), M wherein w(x i) be and x iBelong to of a sort point, M b(x i) be and x iBelong to inhomogeneous point; K is the size of regional area;
Calculate adjacency matrix W in the local class wAdjacency matrix W between the drawn game category b
Calculate local within class scatter matrix R w, dispersion matrix R between local class bWith local population variance degree matrix R t
R w = 1 2 Σ i = 1 n Σ j = 1 n W w ; i , j ( x i - x j ) ( x i - x j ) T ,
R b = 1 2 Σ i = 1 n Σ j = 1 n W b ; i , j ( x i - x j ) ( x i - x j ) T ,
R t=R w+R b
To R tImplement characteristic value decomposition R t=Q Λ Q T, Q=[q wherein 1... q m] form by nonzero eigenvalue characteristic of correspondence vector; Compute matrix Λ -1Q TR bD proper vector b of Q eigenwert maximum 1... b d, be worth descending arrangement by characteristic of correspondence;
Calculate optimum projection matrix A=QB ∈ R D * d, B=[b wherein 1... b d];
With optimum projection matrix with training data project to obtain in the subspace training data in the subspace the picture y i, y i=A Tx i∈ R d(i=1 ... n).
Between the local class that the present invention adopts in the dispersion maximization drawn game category dispersion minimize criterion can make data after through the judgement projection on local granularity the inhomogeneity data separate as far as possible, homogeneous data is polymerization as far as possible, thereby makes the overlapping minimum between inhomogeneity.Cooperate local classifiers (the present invention has adopted the arest neighbors classification) can reach best classification accuracy to implementing classification through the data after the judgement projection.
Beneficial effect of the present invention mainly shows: realize judging the whether specific category of fault and failure judgement of system simultaneously 1..2. the constraint that not distributed by failure classes number and class, applied widely.3. judging nicety rate is higher than existing system.
Description of drawings
Fig. 1 is the structured flowchart of process monitoring of the present invention system;
Fig. 2 is the block diagram that has the host computer concrete structure in the process monitoring of the present invention system.
Embodiment
Below in conjunction with accompanying drawing the present invention is further described.
With reference to Fig. 1, Fig. 2, the present invention is based on the Industrial Process Monitoring system of local discriminatory analysis, comprise the in-situs tester table 2 that is used to gather industrial process object 1 status data, be used for Distributed Control System (DCS) (DCS), server 4 and host computer 7 that industrial process object is controlled, be useful in the server to the host computer transmission and detect data in real time database 5 and the relational database 6 that is used for storing history data in real time, host computer 7 comprises off-line modeling module 8 and on-line monitoring module 15.
Off-line modeling module 8 comprises: training data extraction module 9, first standardized module 10 (being shown as standardized module among the figure) and model training module 11.
And comprise in the model training module 11 that partial structurtes MBM 12, optimum decision projection find the solution module 13 and training data projection module 14.
On-line monitoring module 15 comprises signal acquisition module 16, second standardized module 17 (being shown as standardized module among the figure), fault detection and diagnosis module 18 and display module 19 as a result.
Workflow below in conjunction with each module declaration Industrial Process Monitoring of the present invention system.,
At first the training data extraction module 9, are used for extracting the historical data of industrial process object 1 under data under the normal condition and all kinds of malfunction from the relational database 6 that saves historical data and are used for model training.Set the data variable that needs extraction before extracting, initial sum termination time, the sample number of sampling interval and each data class.Obtain data set And corresponding class sign l 1... l n∈ 0,1,2...c}.With the normal class of class 0 mark, 1,2...c mark failure classes.The data of extracting are reached first standardized module 10.
First standardized module 10, the training data that imports into is done standardization, and making each variable average is zero, and variance is 1, and the data after the standardization are reached model training module 11, it is stand-by that the average of training data and each variable variance reach second standardized module 17.Wherein standardized step is as follows in first standardized module 10:
1) computation of mean values,
2) calculate variance vectors, each element is the variance of each variable.
3) translation is flexible, calculates Wherein " ./" be that vectorial corresponding element is divided by σ xBe the standard deviation vector.
Training data after the standardization that local MBM 12 receptions first standardized module 10 imports into, local geometry and the modeling respectively of local differentiation structure with dispersion logm certificate between local within class scatter matrix drawn game category, and dispersion matrix, population variance degree matrix reach the judgement projection and find the solution module 13 between the local within class scatter matrix that will calculate, local class, and training data is reached fault detection and diagnosis module 18 in training data projection module 14 and the on-line monitoring module 15.
It is as follows to calculate between local within class scatter matrix, local class the step of dispersion matrix and population variance degree matrix:
1) the k-arest neighbors set of data points of each the training data point under the definition of calculating Euclidean distance With M (x i) separated into two parts M w(x i) and M b(x i), M wherein w(x i) be and x iBelong to of a sort point, M b(x i) be and x iBelong to inhomogeneous point.
2) calculate adjacency matrix W in the local class wAdjacency matrix W between the drawn game category b
3) calculate local within class scatter matrix R w, dispersion matrix R between local class bWith local population variance degree matrix R t
R w = 1 2 Σ i = 1 n Σ j = 1 n W w ; i , j ( x i - x j ) ( x i - x j ) T ,
R b = 1 2 Σ i = 1 n Σ j = 1 n W b ; i , j ( x i - x j ) ( x i - x j ) T ,
R t=R w+R b
The dimension d after module 13 is set projection is found the solution in the optimum decision projection, according to the local within class scatter matrix R of the training data that imports into w, dispersion matrix R between local class bWith local population variance degree matrix R tAsk between the local class of the data that make after the projection maximum and optimum decision projection A dispersion minimum in the local class of dispersion, and reach training data projection module 14.Step is as follows:
1) to R tImplement characteristic value decomposition R t=Q Λ Q T, Q=[q wherein 1... q m] form by nonzero eigenvalue characteristic of correspondence vector.
2) compute matrix Λ -1Q TR bThe d of Q eigenwert maximum vectorial b 1... b d, big to minispread by characteristic of correspondence value size.
3) calculate optimum projection matrix A=QB ∈ R D * d, B=[b wherein 1... b d].
Training data projection module 14, the training data that will import into according to the optimum projection matrix A that imports into projects to and obtains y in the subspace i=A Tx i∈ R d(i=1 ... n), and reach fault detection and diagnosis module 18 in the on-line monitoring module 15.
Signal acquisition module 16 is set the time interval of each data sampling, gathers real time data and imports second standardized module 17 into; The data point that training data average that first standardized module 10 imports in second standardized module, the 17 use off-line modeling modules and variance are imported into each signal acquisition module 16 is implemented the flexible processing of translation, and the result is reached fault detection and diagnosis module 18.The data point that standardisation process imported into any time Calculate Wherein " ./" be that vectorial corresponding element is divided by σ xBe the standard deviation vector.
The optimum projection matrix that fault detection and diagnosis module 18 usefulness model training modules 14 are imported into projects in the subspace y=A with the data x that second standardized module 17 imports into TX; Use is sought y based on the nearest neighbor search of Euclidean distance definition 1... y nIn from the nearest some y of y p, judge that the status categories of active procedure is y pStatus categories l pAnd the result reached display module 19 as a result.
The process status classification imported into according to fault detection and diagnosis module 18 of display module 19 shows the state of active procedure on man-machine interface as a result, and display result is that current system is in normal condition or certain malfunction.

Claims (9)

1. Industrial Process Monitoring system based on local discriminatory analysis, the host computer that comprises the measuring instrument, the Distributed Control System (DCS) that are connected with industrial process object, is used for storing the server of Distributed Control System (DCS) institute image data and is used for the image data that processing server stores is characterized in that:
Described server comprises the real-time data base of the real time data that is used for storage industry process object and is used to store the relational database of the data of described industrial process object under normal condition and all kinds of malfunction;
Described host computer comprises off-line modeling module and on-line monitoring module;
Wherein said off-line modeling module comprises:
A) training data extraction module is used for extracting the class sign of data under normal condition and all kinds of malfunction and each data point correspondence as training data from relational database;
B) first standardized module is used for training data is done standardization, obtains the training data x after the standardization i
C) model training module is used for the training data x after the standardization iCarry out the partial structurtes modeling, find the solution again and obtain optimum projection matrix A, with this optimum projection matrix training data is projected to the subspace and obtain the picture y of training data in the subspace i
Described on-line monitoring module comprises:
A) signal acquisition module; Be used for gathering real time data from real-time data base;
B) second standardized module; Be used to utilize the average of training data of first standardized module and variance to collecting to such an extent that real time data does that translation is flexible to be handled, obtain the real time data of translation after flexible;
C) fault detection and diagnosis module; Real time data after flexible projects in the subspace the optimum projection matrix A that obtains with the model training module with translation, obtains the picture y of real time data in the subspace, uses based on the nearest neighbor search of Euclidean distance definition and seeks the picture y of training data in the subspace iIn from picture y nearest some y p, according to y pStatus categories judge the status categories of real time data;
D) display module as a result; Be used to show the status categories of fault detection and diagnosis module judgement.
2. method for diagnosing faults as claimed in claim 1 is characterized in that: the step that described first standardized module is done standardization to training data is as follows:
1) computation of mean values: Wherein Be the data point of extracting from relational database, n is the data point number, and i is the data point call number;
2) calculate variance vectors: σ x 2 = 1 n - 1 Σ i = 1 n ( x ~ i - x ‾ ) ;
3) translation is flexible: Wherein ./be that vectorial corresponding element is divided by σ xBe the standard deviation vector, The average that obtains for step 1).
3. method for diagnosing faults as claimed in claim 1 is characterized in that: described model training module comprises:
The partial structurtes MBM; With the training data x of dispersion matrix between local within class scatter matrix drawn game category after to standardization iLocal geometry and local differentiate structure modeling respectively;
Module is found the solution in the optimum decision projection; Be used for local within class scatter matrix R according to training data w, dispersion matrix R between local class bWith local population variance degree matrix R tTry to achieve between the local class of data after the projection maximum and optimum projection A dispersion minimum in the local class of dispersion;
The training data projection module; With optimum projection matrix A with training data project to obtain in the subspace training data in the subspace the picture y i
4. method for diagnosing faults as claimed in claim 3 is characterized in that: described partial structurtes MBM is with the training data x of dispersion matrix between local within class scatter matrix drawn game category after to standardization iLocal geometry and local differentiate structure modeling respectively, step is as follows:
1) each training data point x under the definition of calculating Euclidean distance iK-arest neighbors set of data points With M (x i) separated into two parts M w(x i) and M b(x i), M wherein w(x i) be and X iBelong to of a sort point, M b(x i) be and X iBelong to inhomogeneous point; K is the size of regional area;
2) calculate adjacency matrix W in the local class wAdjacency matrix W between the drawn game category b
3) calculate local within class scatter matrix R w, dispersion matrix R between local class bWith local population variance degree matrix R t
R w = 1 2 Σ i = 1 n Σ j = 1 n W w ′ , i , j ( x i - x j ) ( x i - x j ) T ,
R b = 1 2 Σ i = 1 n Σ j = 1 n W b ′ , i , j ( x i - x j ) ( x i - x j ) T ,
R t=R w+R b
5. method for diagnosing faults as claimed in claim 4 is characterized in that: module is found the solution in described optimum decision projection, according to the local within class scatter matrix R of training data w, dispersion matrix R between local class bWith local population variance degree matrix R tAsk between the local class of the data that make after the projection maximum and optimum projection matrix A dispersion minimum in the local class of dispersion, step is as follows:
1) to R tImplement characteristic value decomposition R t=Q Λ Q T, Q=[q wherein 1... q m] form by nonzero eigenvalue characteristic of correspondence vector, m is the nonzero eigenvalue number;
2) compute matrix Λ -1Q TR bD proper vector b of Q eigenwert maximum 1... b d, be worth descending arrangement by characteristic of correspondence, the dimensionality reduction subspace dimension of d for setting;
3) calculate optimum projection matrix A=QB ∈ R D * d, B=[b wherein 1... b d], D is the training data dimension;
6. method for diagnosing faults as claimed in claim 5 is characterized in that: described training data projection module projects to training data with optimum projection matrix and obtains the picture y of training data in the subspace in the subspace i, y i=A Tx i∈ R d(i=1 ... n).
7. the Industrial Process Monitoring method based on local discriminatory analysis is characterized in that, comprises the steps:
1) the class sign of the data point under extraction industrial process object normal condition and all kinds of malfunction and each data point correspondence is as training data; Training data is done standardization, obtain the training data x after the standardization iTo the training data x after the standardization iCarry out the partial structurtes modeling, find the solution and obtain optimum projection matrix, with this optimum projection matrix training data is projected to the subspace and obtain the picture y of training data in the subspace i
2) gather the industrial process object real time data; Utilize the average of training data in the step (1) and variance to collecting to such an extent that real time data does that translation is flexible to be handled, obtain the real time data of translation after flexible; The real time data after flexible projects in the subspace with translation to utilize the optimum projection matrix that obtains in the step 1), obtains the picture y of real time data in the subspace, uses based on the nearest neighbor search of Euclidean distance definition and seeks the picture y of training data in the subspace iIn from picture y nearest some y p, according to y pStatus categories judge the status categories of real time data;
3) status categories to real time data shows.
8. method for diagnosing faults as claimed in claim 7 is characterized in that: step (1) is described, and that training data is made the step of standardization is as follows:
1) computation of mean values: Wherein Be the normal condition of industrial process object or the data point under certain class malfunction, n is the number of data point, and i is the data point call number;
2) calculate variance vectors: σ x 2 = 1 n - 1 Σ i = 1 n ( x ~ i - x ‾ ) ;
3) translation is flexible: Wherein ./be that vectorial corresponding element is divided by σ xBe the standard deviation vector, The average that obtains for step 1).
9. method for diagnosing faults as claimed in claim 8 is characterized in that: in the step (1) to the training data x after the standardization iCarry out the partial structurtes modeling, find the solution and obtain optimum projection matrix, with this optimum projection matrix training data is projected to the subspace and obtain the picture y of training data in the subspace iStep comprise:
Calculate each the training data point x under the Euclidean distance definition iK-arest neighbors set of data points With M (x i) separated into two parts M w(x i) and M b(x i), M wherein w(x i) be and X iBelong to of a sort point, M b(x i) be and x iBelong to inhomogeneous point; K is the size of regional area;
Calculate adjacency matrix W in the local class wAdjacency matrix W between the drawn game category b
Calculate local within class scatter matrix R w, dispersion matrix R between local class bWith local population variance degree matrix R t
R w = 1 2 Σ i = 1 n Σ j = 1 n W w ′ , i , j ( x i - x j ) ( x i - x j ) T ,
R b = 1 2 Σ i = 1 n Σ j = i n W b ′ , i , j ( x i - x j ) ( x i - x j ) T ,
R t=R w+R b
To R tImplement characteristic value decomposition R t=Q Λ Q T, Q=[q wherein 1... q m] form by nonzero eigenvalue characteristic of correspondence vector, m is the nonzero eigenvalue number; Compute matrix Λ -1Q TR bD proper vector b of Q eigenwert maximum 1... b d, the dimensionality reduction subspace dimension of d for setting; Be worth descending arrangement by characteristic of correspondence;
Calculate optimum projection matrix A=QB ∈ R D * d, B=[b wherein 1... b d], D is the training data dimension;
With optimum projection matrix with training data project to obtain in the subspace training data in the subspace the picture y i, y i=A Tx i∈ R d(i=1 ... n).
CN2009101550099A 2009-12-10 2009-12-10 System and method for monitoring industrial process based on local discriminatory analysis Expired - Fee Related CN101738998B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN2009101550099A CN101738998B (en) 2009-12-10 2009-12-10 System and method for monitoring industrial process based on local discriminatory analysis

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2009101550099A CN101738998B (en) 2009-12-10 2009-12-10 System and method for monitoring industrial process based on local discriminatory analysis

Publications (2)

Publication Number Publication Date
CN101738998A true CN101738998A (en) 2010-06-16
CN101738998B CN101738998B (en) 2012-05-30

Family

ID=42462568

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2009101550099A Expired - Fee Related CN101738998B (en) 2009-12-10 2009-12-10 System and method for monitoring industrial process based on local discriminatory analysis

Country Status (1)

Country Link
CN (1) CN101738998B (en)

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102591967A (en) * 2011-12-31 2012-07-18 上海昊沧系统控制技术有限责任公司 Method and device for displaying industrial real-time data
CN104503436A (en) * 2014-12-08 2015-04-08 浙江大学 Quick fault detection method based on random projection and k-nearest neighbor method
CN104699050A (en) * 2015-02-13 2015-06-10 浙江中烟工业有限责任公司 Leaf-shred preparation segment on-line monitoring and fault diagnosing method for cigarette filament treatment driven by data
CN104793604A (en) * 2015-04-10 2015-07-22 浙江大学 Principal component tracking based industrial fault monitoring method and application thereof
CN105004542A (en) * 2015-07-15 2015-10-28 浙江中烟工业有限责任公司 Online monitoring and fault diagnosing method for mixing and flavouring process of cigarette filament production based on principal component analysis
CN106384130A (en) * 2016-09-22 2017-02-08 宁波大学 Fault detection method based on data multi-neighbor-local-feature embedding
CN106610345A (en) * 2015-10-27 2017-05-03 北京卫星环境工程研究所 Health state monitoring system for spatial environment simulation test equipment
CN112413814A (en) * 2020-11-04 2021-02-26 武汉科技大学 Online renewable heating ventilation air conditioner sensor fault detection method based on comprehensive distance

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN201017233Y (en) * 2006-11-23 2008-02-06 浙江大学 Manufacturing production process failure diagnosis device based on wavelet analyzing
CN201017232Y (en) * 2006-11-23 2008-02-06 浙江大学 Industry process non-linearity failure diagnosis device based on fisher
CN100470417C (en) * 2006-12-22 2009-03-18 浙江大学 Fault diagnostic system and method for under industrial producing process small sample condition
CN101446831B (en) * 2008-12-30 2011-05-25 东北大学 Decentralized process monitoring method

Cited By (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102591967A (en) * 2011-12-31 2012-07-18 上海昊沧系统控制技术有限责任公司 Method and device for displaying industrial real-time data
CN102591967B (en) * 2011-12-31 2013-07-17 上海昊沧系统控制技术有限责任公司 Method and device for displaying industrial real-time data
CN104503436A (en) * 2014-12-08 2015-04-08 浙江大学 Quick fault detection method based on random projection and k-nearest neighbor method
CN104699050A (en) * 2015-02-13 2015-06-10 浙江中烟工业有限责任公司 Leaf-shred preparation segment on-line monitoring and fault diagnosing method for cigarette filament treatment driven by data
CN104793604A (en) * 2015-04-10 2015-07-22 浙江大学 Principal component tracking based industrial fault monitoring method and application thereof
CN104793604B (en) * 2015-04-10 2017-05-17 浙江大学 Principal component tracking based industrial fault monitoring method and application thereof
CN105004542A (en) * 2015-07-15 2015-10-28 浙江中烟工业有限责任公司 Online monitoring and fault diagnosing method for mixing and flavouring process of cigarette filament production based on principal component analysis
CN106610345A (en) * 2015-10-27 2017-05-03 北京卫星环境工程研究所 Health state monitoring system for spatial environment simulation test equipment
CN106610345B (en) * 2015-10-27 2019-10-25 北京卫星环境工程研究所 A set of health status monitoring system for space-environment facility
CN106384130A (en) * 2016-09-22 2017-02-08 宁波大学 Fault detection method based on data multi-neighbor-local-feature embedding
CN112413814A (en) * 2020-11-04 2021-02-26 武汉科技大学 Online renewable heating ventilation air conditioner sensor fault detection method based on comprehensive distance
CN112413814B (en) * 2020-11-04 2021-11-19 武汉科技大学 Online renewable heating ventilation air conditioner sensor fault detection method based on comprehensive distance

Also Published As

Publication number Publication date
CN101738998B (en) 2012-05-30

Similar Documents

Publication Publication Date Title
CN101738998B (en) System and method for monitoring industrial process based on local discriminatory analysis
US8868985B2 (en) Supervised fault learning using rule-generated samples for machine condition monitoring
CN103914064B (en) Based on the commercial run method for diagnosing faults that multi-categorizer and D-S evidence merge
CN100470417C (en) Fault diagnostic system and method for under industrial producing process small sample condition
CN101464964B (en) Pattern recognition method capable of holding vectorial machine for equipment fault diagnosis
CN107153409B (en) A kind of nongausian process monitoring method based on missing variable modeling thinking
CN106682303B (en) A kind of three-level inverter method for diagnosing faults based on empirical mode decomposition and decision tree RVM
CN107436597B (en) A kind of chemical process fault detection method based on sparse filtering and logistic regression
CN102291392B (en) Hybrid intrusion detection method based on Bagging algorithm
CN107590506A (en) A kind of complex device method for diagnosing faults of feature based processing
CN106843195B (en) The Fault Classification differentiated based on adaptive set at semi-supervised Fei Sheer
CN106201871A (en) Based on the Software Defects Predict Methods that cost-sensitive is semi-supervised
CN106649789B (en) It is a kind of based on the industrial process Fault Classification for integrating semi-supervised Fei Sheer and differentiating
CN105629958B (en) A kind of batch process method for diagnosing faults based on sub-period MPCA SVM
CN104966161A (en) Electric energy quality recording data calculating analysis method based on Gaussian mixture model
CN106482967A (en) A kind of Cost Sensitive Support Vector Machines locomotive wheel detecting system and method
CN109165604A (en) The recognition methods of non-intrusion type load and its test macro based on coorinated training
CN108304567B (en) Method and system for identifying working condition mode and classifying data of high-voltage transformer
CN109298633A (en) Chemical production process fault monitoring method based on adaptive piecemeal Non-negative Matrix Factorization
CN109389325B (en) Method for evaluating state of electronic transformer of transformer substation based on wavelet neural network
CN201017232Y (en) Industry process non-linearity failure diagnosis device based on fisher
CN110837953A (en) Automatic abnormal entity positioning analysis method
CN201035376Y (en) Failure diagnosis device under small sample conditional in the process of manufacturing production
CN110223193A (en) The method of discrimination and system of operation of power networks state are used for based on fuzzy clustering and RS-KNN model
CN109596912A (en) A kind of decomposition method of non-intrusion type power load

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
C17 Cessation of patent right
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20120530

Termination date: 20121210