CN107272655B - Batch process fault monitoring method based on multistage ICA-SVDD - Google Patents

Batch process fault monitoring method based on multistage ICA-SVDD Download PDF

Info

Publication number
CN107272655B
CN107272655B CN201710599054.8A CN201710599054A CN107272655B CN 107272655 B CN107272655 B CN 107272655B CN 201710599054 A CN201710599054 A CN 201710599054A CN 107272655 B CN107272655 B CN 107272655B
Authority
CN
China
Prior art keywords
data
stage
ica
matrix
svdd
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.)
Active
Application number
CN201710599054.8A
Other languages
Chinese (zh)
Other versions
CN107272655A (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.)
Jiangnan University
Original Assignee
Jiangnan University
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 Jiangnan University filed Critical Jiangnan University
Priority to CN201710599054.8A priority Critical patent/CN107272655B/en
Publication of CN107272655A publication Critical patent/CN107272655A/en
Application granted granted Critical
Publication of CN107272655B publication Critical patent/CN107272655B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B23/00Testing or monitoring of control systems or parts thereof
    • G05B23/02Electric testing or monitoring
    • G05B23/0205Electric testing or monitoring by means of a monitoring system capable of detecting and responding to faults
    • G05B23/0218Electric testing or monitoring by means of a monitoring system capable of detecting and responding to faults characterised by the fault detection method dealing with either existing or incipient faults
    • G05B23/0243Electric testing or monitoring by means of a monitoring system capable of detecting and responding to faults characterised by the fault detection method dealing with either existing or incipient faults model based detection method, e.g. first-principles knowledge model

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Automation & Control Theory (AREA)
  • Complex Calculations (AREA)

Abstract

The invention discloses a kind of batch process fault monitoring methods based on multistage ICA-SVDD.For process mechanism complexity and there are the batch processes of multiple operational phases.The multistage negotiation and data distribution non-Gaussian system problem having for some batch processes, using a kind of improved divided stages and fault monitoring method.Divided stages are carried out according to the similarity of each timeslice and K mean algorithm first, then the characteristic information that Independent Component Analysis extracts non-gaussian is utilized respectively to each stage, it is finally introducing Support Vector data description algorithm and Statistic analysis models is established respectively to independent element and remaining Gauss residual error space, realize the malfunction monitoring to whole process.Malfunction monitoring applied to an actual conductor etching process, the results showed that this method has more preferably monitoring effect to multistage batch process.

Description

Batch process fault monitoring method based on multistage ICA-SVDD
Technical field
The present invention relates to the batch process fault monitoring methods based on multistage ICA-SVDD, belong to industrial process failure and examine Disconnected and hard measurement field.
Background technique
Batch process is a kind of important industrial mode of production, and process mechanism is complicated and there are multiple operation ranks Section, and influence of the product quality vulnerable to uncertain factor.In order to guarantee batch production process safe and reliable operation and The high quality of product is pursued, and is needed to establish effective process monitoring system and is carried out failure monitoring to batch production process.Polynary system Meter course control method for use has been widely applied in Batch process monitoring, such as multidirectional pivot analysis (multiway Principal component analysis, MPCA) and multidirectional partial least squares analysis (multiway partial least Squares, MPLS).But most of data description method relevant to principal component analysis and partial least squares analysis have number It is linear limitation according to the relationship met between Gaussian Profile and different variables.
For the non-Gaussian feature of batch process data, independent component analysis (independent component Analysis, ICA) it is introduced in process monitoring field, in order to adapt to the malfunction monitoring to batch process, some scholars are ICA Method is extended to multidirectional independent component analysis (multiway ICA, MICA) method.Although can handle non-gaussian data, It is to be based on Density Estimator method calculating process complexity in the confidence limit of determination process monitoring statisticss amount and parameter can not It is accurate to obtain, when dimension is larger, the problems such as not can avoid Density Estimator bring " dimension disaster ".In addition, for The monitoring of nonlinear batch process, traditional MPCA/MPLS method also extend to its non-linear form, such as multidirectional core PCA and more To core PLS, however, the course monitoring method based on these nonlinear methods is also required to Gaussian distributed.
The multioperation stage is an inherent characteristic of many batch processes, if realizing entire batch using single modeling pattern It is bad to will lead to monitoring effect of the model in different phase for the monitoring of secondary process.For this characteristic, domestic and foreign scholars are Done a large amount of research, by batch process carry out reasonable divided stages and in sub-stage establishment process monitoring Model, to improve Monitoring Performance.
Process monitoring may be considered a monodrome classification problem, because the task of monitoring is by normal data and number of faults According to separation.Support Vector data description (Support vector data description, SVDD) algorithm is a kind of initial The monodrome classification method proposed by Tax and Duin.Normal data sample space is mapped to high dimensional feature by nonlinear transformation A model is simultaneously established in space, so that normal data be separated with fault data, achievees the purpose that procedure fault monitors.It uses SVDD algorithm carry out malfunction monitoring can handle simultaneously do not meet be between Gaussian Profile and variable non-linear relation data. SVDD has been used for the fields such as damage check, image classification, pattern-recognition, and the application in process monitoring field also starts Paid attention to.
A kind of fault monitoring method of the multistage batch process based on independent component analysis and Support Vector data description, It can effectively improve the malfunction monitoring performance of batch process.Due to most of with principal component analysis and partial least squares analysis phase It is linear limitation that the data description method of pass, which has the relationship between data fit Gaussian Profile and different variables,.The party Method can solve process data non-gaussian and nonlinear monitoring problem simultaneously, have more preferably malfunction monitoring effect.
Summary of the invention
It is directed to multistage, non-Gaussian system and non-linear, the process mechanism complexity of batch process presentation, and product quality Vulnerable to the influence of uncertain factor, in order to improve the malfunction monitoring performance to multistage batch process, the present invention provides one kind Batch process fault monitoring method based on multistage ICA-SVDD.
Three-dimensional data is unfolded along batch direction first, fuzzy divided stages are carried out according to the similarity of each timeslice, Available initial cluster number carries out accurate divided stages will pass through K mean algorithm;After determining the stage, then Three-dimensional data is unfolded along variable direction, is then utilized respectively ICA method to each stage and carries out feature extraction, extracts corresponding Non-gaussian characteristic information and residual information;Finally the non-gaussian space of independent element and remaining residual error space are passed through respectively SVDD algorithm establishes Statistic analysis models, realizes the on-line fault monitoring to whole process.
The purpose of the present invention is what is be achieved through the following technical solutions:
Batch process fault monitoring method based on multistage ICA-SVDD, the method includes following procedure: being directed to The multioperation stage feature of batch process needs to carry out reasonable divided stages to production process, to establish multiple submodels Carry out malfunction monitoring.Fuzzy divided stages are carried out to production process according to the similarity of the mean vector of each timeslice first, Obtain preliminary number of stages.Then it is classified as one kind at the time of by K mean algorithm that data characteristics is similar, and then is obtained more Accurate divided stages.
The characteristic information that Independent Component Analysis extracts non-gaussian is utilized respectively to each stage, when using ICA algorithm It after extracting whole independent elements, is rearranged according to non-gaussian degree size, the stronger independent element of d independence before choosing Obtain homographyDue to having carried out multistage division to batch process, so will carry out ICA to each stage Analysis is to carry out feature extraction.And then corresponding non-gaussian characteristic information of each stage and residual information can be extracted, so as to Monitoring and statistics model is established to carry out malfunction monitoring.
It is finally introducing Support Vector data description algorithm and system is established respectively to independent element and remaining Gauss residual error space Analysis model is counted, realizes the malfunction monitoring to whole process.Normal data sample space is mapped to height by nonlinear transformation Dimensional feature space simultaneously establishes a model, so that normal data be separated with fault data, reaches the mesh of procedure fault monitoring 's.Carrying out malfunction monitoring using SVDD algorithm can handle that not meet be non-linear relation between Gaussian Profile and variable simultaneously Data.
Detailed description of the invention
Fig. 1 is the Batch process monitoring method flow diagram based on multistage ICA-SVDD;
The expansion of Fig. 2 batch process data;
Fig. 3 multistage division result;
The normal batch monitoring result of Fig. 4;
Fig. 5 failure TCP+50 batch monitoring result;
Specific embodiment
Below with reference to shown in Fig. 1, the present invention is further described:
The data used is collected during an actual semiconductor etch processes, is half-and-half led respectively The normal data and fault data of body etching carry out malfunction monitoring.
Step 1: the 3-D data set X (I × J × K) of batch process carries out two-dimensional development, wherein I represents lot number, J generation Table variable number, K represent sampling number.Using along batch direction and the data processing method combined along variable direction, first by three The data X (I × J × K) of dimension form is converted into two-dimensional matrix X (I × KJ) along batch direction, then standardizes two-dimensional matrix;Again It is reconfigured according to variable direction, forms new two-dimensional matrix X (KI × J).Two step data method of deploying are as shown in Figure 2.
Step 2: reasonable divided stages being carried out to production process, carry out malfunction monitoring to establish multiple submodels.It is first Fuzzy divided stages are first carried out to production process according to the similarity of the mean vector of each timeslice, obtain the preliminary stage Number.Then it is classified as one kind at the time of by K mean algorithm that data characteristics is similar, and then obtains more accurate divided stages.
Step1: first it is unfolded three dimensional process data X (I × J × K) to obtain two-dimensional matrix X (I × KJ) by batch direction, so The 2-D data timeslice matrix X that temporally axis direction is cut into batch afterwards and variable formsk(I × J), k=1,2 ..., K.
Step2: each timeslice matrix X is soughtkThe mean vector of (I × J), is denoted asThese mean values Vector represents the characteristic information of each timeslice, carries out initial stage division to timeslice using these characteristic informations, goes forward side by side The identification in row each stage, with first timeslice X1Benchmark X as first stagebase, then according to similarity calculation Formula:
Successively calculate XbaseSubsequent timeslice and its similarity, and similarity threshold α is set, if X2And XbasePhase It is greater than threshold alpha like degree, then it is assumed that X2Present period is also belonged to, then proceedes to calculate next timeslice and XbaseSimilarity; Otherwise it is assumed that X2Belong to next stage, and enables Xbase=X2, continue by above-mentioned steps.
Similar timeslice is connected to form a period according to similarity, obtains preliminary fuzzy division.It can obtain To corresponding stage number P, this provides foundation to choose cluster numbers with clustering algorithm below, still, this fuzzy division method Will appear certain points or few continuity point can not accurately be divided into some stage.
Step3: being clustered by mean vector of the K-means clustering algorithm to timeslice, algorithm input be mean value to Duration setAnd cluster number P, P cluster centre is arbitrarily selected, successive ignition calculating is carried out, when algorithm is full When the sufficient condition of convergence, the cluster centre of available P subclass calculates each mean vectorTo all cluster centres away from From, so that it may it obtainsFor the membership of P subclass, since the input of clustering algorithm arranges sequentially in time Timeslice mean vector, therefore sequentially in time, the point in affiliated stage can will can not be determined in fuzzy division, be divided into In one corresponding stage, so that it may obtain more accurate divided stages.The multistage division result of conductor etching process is shown in Fig. 3.
Step 3: carrying out feature letter using independent component analysis (independent component analysis, ICA) Breath extracts, and data higher-order statistics are more fully utilized in ICA, and can further extract phase from observation data Mutual independent latent variable, these latent variables can more constitutionally extraction reaction process features.
ICA model is defined as
X=AS+E (2)
Wherein X=[x (1), x (2) ..., x (n)] ∈ Rm×nIt is observation data matrix, A=[a1,a2,...,ad]∈Rm×d It is unknown hybrid matrix, S=[s (1), s (2) ..., s (n)] ∈ Rd×nIt is hiding independent element matrix, E ∈ Rm×nIt is Residual matrix.N is the number of samples of acquisition, by d≤m it is found that a kind of ICA and data compression technique similar with PCA in fact, Information as much as possible is described by data as few as possible.
The purpose of ICA is to estimate hybrid matrix W and independent element S from observation data X, therefore, ICA target: is found One solves mixed matrix W, source signal can be isolated from observation signal, i.e.,
As d=m, that is, the inverse time that mixed matrix W is hybrid matrix A is solved,It is exactly the best estimate of independent element S.
It after extracting whole independent elements using ICA algorithm, is rearranged according to non-gaussian degree size, d before choosing The stronger independent element of independence obtains homographyDue to having carried out multistage division to batch process, so wanting To the X in each stagep(KpI × J) ICA analysis is carried out to carry out feature extraction, wherein p indicates corresponding pth stage, XpTable Show the two-dimensional matrix that the pth stage is unfolded along variable direction, KpIt indicates pth stage corresponding number of samples, and then can extract Each stage corresponding non-gaussian characteristic information and residual information carry out malfunction monitoring to establish monitoring and statistics model.
Step 4: being carried out using Support Vector data description (Support vector data description, SVDD) The on-line monitoring of batch process, SVDD are a kind of monodrome sorting algorithms, and basic thought is for training dataset X={ xi,i =1 ..., N }, by Φ: X → F of non-linear transfer by luv space data projection to feature space { Φ (xi), i=1 ..., N }, the suprasphere that may include the minimum volume of all data samples then can be found in feature space, SVDD passes through core The input space is mapped to higher dimensional space to learn to obtain flexible and accurate data descriptive model, obtained hypersphere by function Volume data, which describes boundary, to be indicated by the supporting vector of sub-fraction.
In order to construct such minimal hyper-sphere, SVDD needs to solve following optimization problem:
In formula (4), a be suprasphere the centre of sphere, R be suprasphere radius, penalty coefficient C weighed suprasphere volume and The false segmentation rate of training sample, slack variable ξiIntroducing represent penalty term accidentally point generated to i-th training sample, it is above-mentioned excellent Change problem, which can be converted into, solves corresponding dual problem:
It is to introduce kernel function K (x hereini,xj) replace interior Product function (xi,xj), using quadratic programming, a can be found outi, such as Fruit x0An arbitrary supporting vector is represented, then the centre of sphere of suprasphere and radius may be expressed as:
For new sample xnew, its available distance for arriving the suprasphere centre of sphere:
If DistnewThen the sample is normal by≤R;Conversely, the sample is exceptional sample.
After carrying out feature extraction using each stage of the ICA to batch process, the non-height of independent element can be respectively obtained The data in this space and remaining residual error space, by SVDD method to the independent element extractedIt is established with residual matrix E Statistic analysis models, first against independent elementIt is as follows to establish SVDD model:
In formula:Indicate k-th of sample, αkFor the Lagrange multiplier of corresponding sample, K () indicates Gaussian kernel letter Number, available independent elementThe centre of sphere and radius for being formed by SVDD suprasphere are as shown in (9) formula:
Equally, then to residual matrix E it is as follows to establish SVDD model:
So, the centre of sphere and radius that residual matrix E is formed by SVDD suprasphere are as shown in (11) formula:
For each stage of batch process, the suprasphere radius R in corresponding stage, online process monitoring can be respectively obtained When, as the sampled data x for obtaining current timenewWhen, after data prediction and ICA feature extraction, independence is sought respectively IngredientWith residual matrix EnewTo the distance of the corresponding centre of sphereAnd DistE_new, therefore as Dist≤R, it can recognize Sample for current time is normal, and works as Dist > R, then it represents that the sample at current time is fault data.
After carrying out divided stages and pretreatment to process data, for the failure of relatively simple modeling and multi-stage modeling Monitoring effect establishes the Statistical monitor model of single and multistage PCA, ICA, ICA-SVDD respectively, all statistics Statisti-cal control limit is set as 99%.Fig. 4 gives six kinds of methods to the monitored results of normal batch, it can be seen that six kinds of sides Method can obtain good monitored results to normal batch.Wherein (a), (b), the event that (c) is the normal batch of single modeling method Barrier monitoring figure, (d), (e), (f) be the normal batch of multi-stage modeling method malfunction monitoring figure.Fig. 5 gives TCP+50 failure Malfunction monitoring figure, wherein (a), (b), (c) be single modeling method failure batch monitoring figure, (d), (e), (f) be it is multistage The monitoring figure of section modeling method failure batch.From figure 5 it can be seen that either single model is still for TCP+50 failure The malfunction monitoring of multiphase confinement, ICA-SVDD are to be substantially better than other two methods.

Claims (2)

1. the method for the batch process malfunction monitoring based on multistage ICA-SVDD, which is characterized in that this method step are as follows:
Step 1: the 3-D data set X (I × J × K) of batch process carries out two-dimensional development, wherein I represents lot number, and J, which is represented, to be become Number is measured, K represents sampling number;Using along batch direction and the data processing method combined along variable direction, first by three-dimensional shaped The data X (I × J × K) of formula is converted into two-dimensional matrix X (I × KJ) along batch direction, then standardizes two-dimensional matrix;According still further to Variable direction reconfigures, and forms new two-dimensional matrix X (KI × J);
Step 2: reasonable divided stages being carried out to production process, carry out malfunction monitoring to establish multiple submodels;Root first Fuzzy divided stages are carried out to production process according to the similarity of the mean vector of each timeslice, obtain preliminary number of stages; Then it is classified as one kind at the time of by K mean algorithm that data characteristics is similar, and then obtains more accurate divided stages;
Step1: first it is unfolded three dimensional process data X (I × J × K) to obtain two-dimensional matrix X (I × KJ) by batch direction, then press Time-axis direction is cut into the 2-D data timeslice matrix X of batch and variable compositionk(I × J), k=1,2 ..., K;
Step2: each timeslice matrix X is soughtkThe mean vector of (I × J), is denoted asThese mean vectors The characteristic information of each timeslice is represented, initial stage division is carried out to timeslice using these characteristic informations, and carry out each The identification in a stage, with first timeslice X1Benchmark X as first stagebase, then according to calculating formula of similarity:
Successively calculate XbaseSubsequent timeslice and its similarity, and similarity threshold α is set, if X2And XbaseSimilarity Greater than threshold alpha, then it is assumed that X2Present period is also belonged to, then proceedes to calculate next timeslice and XbaseSimilarity;Otherwise, Think X2Belong to next stage, and enables Xbase=X2, continue by above-mentioned steps;
Similar timeslice is connected to form a period according to similarity, obtains preliminary fuzzy division;It is available right The stage number P answered, this provides foundation to choose cluster numbers with clustering algorithm below, and still, this fuzzy division method can go out Existing certain points or few continuity point can not accurately be divided into some stage;
Step3: being clustered by mean vector of the K-means clustering algorithm to timeslice, and algorithm input is mean vector collection It closesAnd cluster number P, P cluster centre is arbitrarily selected, successive ignition calculating is carried out, is restrained when algorithm meets When condition, the cluster centre of available P subclass calculates each mean vectorTo the distance of all cluster centres, so that it may To obtainFor the membership of P subclass, since the input of clustering algorithm is that the timeslice that arranges sequentially in time is equal Be worth vector, therefore sequentially in time, can will can not be determined in fuzzy division belonging to the stage point, be divided into one it is corresponding In stage, so that it may obtain more accurate divided stages;
Step 3: carrying out characteristic information using independent component analysis (independent component analysis, ICA) and mention It takes, data higher-order statistics are more fully utilized in ICA, and can further extract from observation data mutually indepedent Latent variable, these latent variables can more constitutionally extract reaction process feature;
ICA model is defined as
X=AS+E (2)
Wherein X=[x (1), x (2) ..., x (n)] ∈ Rm×nIt is observation data matrix, A=[a1,a2,...,ad]∈Rm×dIt is Unknown hybrid matrix, S=[s (1), s (2) ..., s (n)] ∈ Rd×nIt is hiding independent element matrix, E ∈ Rm×nIt is residual error Matrix;N be acquisition number of samples, by d≤m it is found that ICA in fact it is similar with PCA be also a kind of data compression technique, by use up May data less information as much as possible described;
The purpose of ICA is to estimate hybrid matrix A and independent element S from observation data X, therefore, ICA target: finds one Mixed matrix W is solved, source signal can be isolated from observation signal, i.e.,
As d=m, that is, the inverse time that mixed matrix W is hybrid matrix A is solved,It is exactly the best estimate of independent element S;
It after extracting whole independent elements using ICA algorithm, is rearranged according to non-gaussian degree size, selection first d independent The stronger independent element of property obtains homographyDue to having carried out multistage division to batch process, so will be to every The X in a stagep(KpI × J) ICA analysis is carried out to carry out feature extraction, wherein p indicates corresponding pth stage, XpIndicate pth The two-dimensional matrix that stage is unfolded along variable direction, KpIt indicates pth stage corresponding number of samples, and then each rank can be extracted The corresponding non-gaussian characteristic information of section and residual information carry out malfunction monitoring to establish monitoring and statistics model;
Step 4: interval is carried out using Support Vector data description (Support vector data description, SVDD) The on-line monitoring of process, SVDD are a kind of monodrome sorting algorithms, and basic thought is for training dataset X={ xi, i= 1 ..., N }, by Φ: X → F of non-linear transfer by luv space data projection to feature space { Φ (xi), i=1 ..., N }, Then the suprasphere that may include the minimum volume of all data samples can be found in feature space, SVDD passes through kernel function The input space is mapped to higher dimensional space to learn to obtain flexible and accurate data descriptive model, obtained hypersphere volume data Describing boundary is indicated by the supporting vector of sub-fraction, and in order to construct such minimal hyper-sphere, SVDD needs Solve following optimization problem:
In formula (4), a is the centre of sphere of suprasphere, and R is the radius of suprasphere, and penalty coefficient C has weighed the volume and training of suprasphere The false segmentation rate of sample, slack variable ξiIntroducing represent penalty term accidentally point, above-mentioned optimization problem generated to i-th training sample It can be converted into and solve corresponding dual problem:
Wherein, αiRepresent the Lagrange multiplier of corresponding i-th of training sample, αjThe glug for representing corresponding j-th of training sample is bright Day multiplier;
It is to introduce kernel function K (x hereini,xj) replace interior Product function (xi,xj), using quadratic programming, a can be found outiIf x0 An arbitrary supporting vector is represented, then the centre of sphere of suprasphere and radius may be expressed as:
For new sample xnew, its available distance for arriving the suprasphere centre of sphere:
If DistnewThen the sample is normal by≤R;Conversely, the sample is exceptional sample;
After carrying out feature extraction using each stage of the ICA to batch process, the non-gaussian that can respectively obtain independent element is empty Between and remaining residual error space data, by SVDD method to the independent element extractedIt establishes and counts with residual matrix E Analysis model, first against independent elementIt is as follows to establish SVDD model:
In formula:Indicate k-th of sample, αkFor the Lagrange multiplier of corresponding sample, K () indicates gaussian kernel function, can To obtain independent elementThe centre of sphere and radius for being formed by SVDD suprasphere are as shown in (9) formula:
Equally, then to residual matrix E it is as follows to establish SVDD model:
Wherein, βkRepresent the Lagrange multiplier of corresponding k-th of training sample;
So, the centre of sphere and radius that residual matrix E is formed by SVDD suprasphere are as shown in (11) formula:
For each stage of batch process, the suprasphere radius R in corresponding stage can be respectively obtained, when online process monitoring, As the sampled data x for obtaining current timenewWhen, after data prediction and ICA feature extraction, independent element is sought respectivelyWith residual matrix EnewTo the distance of the corresponding centre of sphereAnd DistE_new, therefore as Dist≤R, it is believed that when The sample at preceding moment is normal, and works as Dist > R, then it represents that the sample at current time is fault data.
2. the method for the batch process malfunction monitoring according to claim 1 based on multistage ICA-SVDD, feature exist In, the problem of there is multistage negotiation and data distribution non-Gaussian system due to batch process, based on independent component analysis and support to The fault monitoring method for measuring the multistage batch process of data description, can solve process data non-gaussian and nonlinear simultaneously Monitoring problem.
CN201710599054.8A 2017-07-21 2017-07-21 Batch process fault monitoring method based on multistage ICA-SVDD Active CN107272655B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710599054.8A CN107272655B (en) 2017-07-21 2017-07-21 Batch process fault monitoring method based on multistage ICA-SVDD

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710599054.8A CN107272655B (en) 2017-07-21 2017-07-21 Batch process fault monitoring method based on multistage ICA-SVDD

Publications (2)

Publication Number Publication Date
CN107272655A CN107272655A (en) 2017-10-20
CN107272655B true CN107272655B (en) 2019-08-06

Family

ID=60078476

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710599054.8A Active CN107272655B (en) 2017-07-21 2017-07-21 Batch process fault monitoring method based on multistage ICA-SVDD

Country Status (1)

Country Link
CN (1) CN107272655B (en)

Families Citing this family (17)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109177101A (en) * 2018-06-28 2019-01-11 浙江工业大学 A kind of injection molding machine batch process fault detection method
CN108960329B (en) * 2018-07-06 2020-11-06 浙江科技学院 Chemical process fault detection method containing missing data
CN109062196B (en) * 2018-10-31 2020-12-15 东北大学 Blast furnace process monitoring and fault diagnosis method integrated with PCA-ICA
CN109669413B (en) * 2018-12-13 2021-01-08 宁波大学 Dynamic non-Gaussian process monitoring method based on dynamic latent independent variables
CN109754010B (en) * 2018-12-29 2021-04-02 北京化工大学 Intermittent process multi-mode partitioning method for time-series constraint fuzzy clustering
CN110245460B (en) * 2019-06-28 2023-04-25 北京工业大学 Intermittent process fault monitoring method based on multi-stage OICA
CN110701487B (en) * 2019-09-18 2021-08-24 浙江工业大学 KPCA and Cas-SVDD-based multi-working-condition pipeline leakage detection method
CN111122135B (en) * 2019-12-05 2021-05-28 西安交通大学 Method for evaluating looseness degree of flange bolt connection structure
CN111160811B (en) * 2020-01-17 2023-10-03 北京工业大学 Batch process fault monitoring method based on multi-stage FOM-SAE
CN111736567B (en) * 2020-05-12 2021-10-26 江南大学 Multi-block fault monitoring method based on fault sensitivity slow characteristic
CN112418348A (en) * 2020-12-11 2021-02-26 大连理工大学 Image source identification method based on envelope optimization
CN113190537A (en) * 2021-03-22 2021-07-30 广东电网有限责任公司东莞供电局 Data characterization method for emergency repair site in monitoring area
CN113311796B (en) * 2021-06-04 2022-04-22 北京工业大学 Fermentation process stage division method based on joint typical variable matrix
CN114280935A (en) * 2021-12-16 2022-04-05 北京工业大学 Multi-stage fermentation process fault monitoring method based on semi-supervised FCM and SAE of information entropy
CN116028783B (en) * 2023-03-30 2023-06-20 北京科技大学 Strip steel hot continuous rolling micro fault real-time detection method and device based on process data
CN116776258B (en) * 2023-08-24 2023-10-31 北京天恒安科集团有限公司 Power equipment monitoring data processing method and system
CN117350328B (en) * 2023-09-11 2024-06-25 江南大学 LSTM-CVAE-based lactobacillus fermentation process fault detection and diagnosis method

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103777627A (en) * 2014-01-24 2014-05-07 浙江大学 Batch process online-monitoring method based on small number of batches
CN103970092A (en) * 2014-04-13 2014-08-06 北京工业大学 Multi-stage fermentation process fault monitoring method based on self-adaption FCM algorithm
CN104122485A (en) * 2014-07-23 2014-10-29 国网天津市电力公司 Recording file based line fault analysis
CN104699894A (en) * 2015-01-26 2015-06-10 江南大学 JITL (just-in-time learning) based multi-model fusion modeling method adopting GPR (Gaussian process regression)
CN105159247A (en) * 2015-08-05 2015-12-16 北京工业大学 Information transmission based phase affiliation judgment method for real time sampling points in intermittent process
CN105425779A (en) * 2015-12-24 2016-03-23 江南大学 ICA-PCA multi-working condition fault diagnosis method based on local neighborhood standardization and Bayesian inference
CN106300371A (en) * 2016-08-29 2017-01-04 武汉理工大学 A kind of low-voltage ride-through method of wound brushless double-fed wind power generator group
CN106483847A (en) * 2016-09-20 2017-03-08 北京工业大学 A kind of handpiece Water Chilling Units fault detection method based on self adaptation ICA

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8849380B2 (en) * 2007-11-26 2014-09-30 Canfield Scientific Inc. Multi-spectral tissue imaging

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103777627A (en) * 2014-01-24 2014-05-07 浙江大学 Batch process online-monitoring method based on small number of batches
CN103777627B (en) * 2014-01-24 2016-03-30 浙江大学 A kind of batch process on-line monitoring method based on a small amount of batch
CN103970092A (en) * 2014-04-13 2014-08-06 北京工业大学 Multi-stage fermentation process fault monitoring method based on self-adaption FCM algorithm
CN104122485A (en) * 2014-07-23 2014-10-29 国网天津市电力公司 Recording file based line fault analysis
CN104699894A (en) * 2015-01-26 2015-06-10 江南大学 JITL (just-in-time learning) based multi-model fusion modeling method adopting GPR (Gaussian process regression)
CN105159247A (en) * 2015-08-05 2015-12-16 北京工业大学 Information transmission based phase affiliation judgment method for real time sampling points in intermittent process
CN105425779A (en) * 2015-12-24 2016-03-23 江南大学 ICA-PCA multi-working condition fault diagnosis method based on local neighborhood standardization and Bayesian inference
CN106300371A (en) * 2016-08-29 2017-01-04 武汉理工大学 A kind of low-voltage ride-through method of wound brushless double-fed wind power generator group
CN106483847A (en) * 2016-09-20 2017-03-08 北京工业大学 A kind of handpiece Water Chilling Units fault detection method based on self adaptation ICA

Also Published As

Publication number Publication date
CN107272655A (en) 2017-10-20

Similar Documents

Publication Publication Date Title
CN107272655B (en) Batch process fault monitoring method based on multistage ICA-SVDD
Luo et al. Unsupervised learning of long-term motion dynamics for videos
Kantorov et al. Efficient feature extraction, encoding and classification for action recognition
CN109116834B (en) Intermittent process fault detection method based on deep learning
CN109271975A (en) A kind of electrical energy power quality disturbance recognition methods based on big data multi-feature extraction synergetic classification
CN109633369B (en) Power grid fault diagnosis method based on multi-dimensional data similarity matching
CN112200104B (en) Chemical engineering fault diagnosis method based on novel Bayesian framework for enhanced principal component analysis
CN109389325B (en) Method for evaluating state of electronic transformer of transformer substation based on wavelet neural network
CN106952293B (en) Target tracking method based on nonparametric online clustering
CN105955214B (en) Batch process fault detection method based on sample time-series and neighbour's affinity information
Ren et al. Power system event classification and localization using a convolutional neural network
CN106600602A (en) Clustered adaptive window based hyperspectral image abnormality detection method
CN114676742A (en) Power grid abnormal electricity utilization detection method based on attention mechanism and residual error network
Zhu et al. Networked time series shapelet learning for power system transient stability assessment
CN110045227B (en) power distribution network fault diagnosis method based on random matrix and deep learning
CN112507778B (en) Loop detection method of improved bag-of-words model based on line characteristics
CN109829494A (en) A kind of clustering ensemble method based on weighting similarity measurement
CN108375729B (en) Degenerate state appraisal procedure is clustered based on the motor SOM that Fisher differentiates
CN110866439A (en) Hyperspectral image joint classification method based on multi-feature learning and superpixel kernel sparse representation
CN101738998A (en) System and method for monitoring industrial process based on local discriminatory analysis
Song et al. Feature extraction and target recognition of moving image sequences
Gu et al. Identification of concurrent control chart patterns with singular spectrum analysis and learning vector quantization
Yang et al. An intelligent singular value diagnostic method for concrete dam deformation monitoring
CN117041017A (en) Intelligent operation and maintenance management method and system for data center
CN115375921A (en) Two-stage non-intrusive load identification method and terminal

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant