CN110008565A - A kind of industrial process unusual service condition prediction technique based on operating parameter association analysis - Google Patents
A kind of industrial process unusual service condition prediction technique based on operating parameter association analysis Download PDFInfo
- Publication number
- CN110008565A CN110008565A CN201910244872.5A CN201910244872A CN110008565A CN 110008565 A CN110008565 A CN 110008565A CN 201910244872 A CN201910244872 A CN 201910244872A CN 110008565 A CN110008565 A CN 110008565A
- Authority
- CN
- China
- Prior art keywords
- data
- parameter
- relevance
- window
- prediction
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 46
- 238000004519 manufacturing process Methods 0.000 title claims abstract description 23
- 238000012098 association analyses Methods 0.000 title abstract 4
- 238000009499 grossing Methods 0.000 claims abstract description 15
- 230000002159 abnormal effect Effects 0.000 claims description 29
- 230000011218 segmentation Effects 0.000 claims description 29
- 238000010219 correlation analysis Methods 0.000 claims description 13
- 238000005070 sampling Methods 0.000 claims description 12
- 238000004364 calculation method Methods 0.000 claims description 10
- 230000005856 abnormality Effects 0.000 claims description 9
- 239000003102 growth factor Substances 0.000 claims description 6
- 239000011159 matrix material Substances 0.000 claims description 6
- 238000005259 measurement Methods 0.000 claims description 5
- 238000004458 analytical method Methods 0.000 claims description 4
- 238000013461 design Methods 0.000 claims description 3
- 238000012417 linear regression Methods 0.000 claims description 3
- 238000012544 monitoring process Methods 0.000 description 3
- 238000003745 diagnosis Methods 0.000 description 2
- 238000003801 milling Methods 0.000 description 2
- FGUUSXIOTUKUDN-IBGZPJMESA-N C1(=CC=CC=C1)N1C2=C(NC([C@H](C1)NC=1OC(=NN=1)C1=CC=CC=C1)=O)C=CC=C2 Chemical compound C1(=CC=CC=C1)N1C2=C(NC([C@H](C1)NC=1OC(=NN=1)C1=CC=CC=C1)=O)C=CC=C2 FGUUSXIOTUKUDN-IBGZPJMESA-N 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000006731 degradation reaction Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 239000007769 metal material Substances 0.000 description 1
- 208000024891 symptom Diseases 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/20—Administration of product repair or maintenance
Landscapes
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Physics & Mathematics (AREA)
- Human Resources & Organizations (AREA)
- Theoretical Computer Science (AREA)
- Strategic Management (AREA)
- Economics (AREA)
- General Physics & Mathematics (AREA)
- Marketing (AREA)
- Entrepreneurship & Innovation (AREA)
- General Business, Economics & Management (AREA)
- Tourism & Hospitality (AREA)
- Quality & Reliability (AREA)
- Operations Research (AREA)
- Geometry (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Development Economics (AREA)
- Game Theory and Decision Science (AREA)
- General Engineering & Computer Science (AREA)
- Testing And Monitoring For Control Systems (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The industrial process unusual service condition prediction technique based on operating parameter association analysis that the invention discloses a kind of, can be applied to the prognostic and health management of industrial process.The present invention starts with from the relevance between industrial process operating parameter, and the association analysis based on operating parameter carries out unusual service condition prediction.In one-parameter forecast period, the present invention predicts each operating parameter by exponential smoothing method according to existing sensing data.In the association analysis stage, the present invention calculates operating parameter relevance by known each operational parameter value and parameter prediction value, and wherein parameter association is indicated by a series of similitude of indexs of representation parameter curve.In the relevance trend prediction stage, the present invention constructs multivariate autoregressive model and predicts parameter association.Method proposed by the present invention accounts for operating parameter relevance, and the prediction result that can be obtained more complete device exception information and more shift to an earlier date is of practical significance for the failure predication of industrial equipment.
Description
Technical Field
The invention belongs to the technical field of reliability engineering, and relates to an industrial process abnormal working condition prediction method based on operation parameter correlation analysis.
Background
With the continuous emergence of complex systems and the increasing demand of real-time monitoring of industrial processes, modern industrial equipment is often equipped with a plurality of sensors to monitor the operation state of the industrial equipment in the operation process. Meanwhile, multiple fault modes may occur in the operation process of the equipment, a certain fault may correspond to a plurality of symptoms, and under the condition, the single sensor information cannot completely reflect the operation state of the equipment, so that fault prediction based on multi-sensor information is generated at the right moment. The failure prediction based on multi-sensor information aims to analyze the operation state of the equipment using comprehensive sensor information, thereby making more reliable equipment diagnosis and prediction. With the continuous development of sensing technology, the use of multiple sensors for condition monitoring, fault diagnosis and prediction of equipment has become a trend.
For a plurality of sensors during the operation of the equipment, the represented operation parameters do not exist independently, and the change of each operation parameter during the operation of the equipment is actually the reaction to the current operation state of the equipment. Under the normal operation state of the industrial equipment, the operation parameters are usually relatively stable and maintained at a relatively stable level, so that the relevance of the operation parameters is relatively stable. However, when an abnormal working condition occurs, the response of each operation parameter to the abnormality is different, so that the relevance of the operation parameters changes, and the relevance and the trend of the relevance change necessarily imply the information of the abnormality and even the fault of the equipment.
Disclosure of Invention
Aiming at the technical situation in the prior art, the invention aims to judge the state of the equipment in the operation process through the relevance of the operation parameters and predict the abnormal working condition of the equipment through predicting the relevance change trend of the operation parameters aiming at the problem that the operation parameters of the equipment have relevance in the operation process.
The concept of the present invention will now be explained as follows:
the invention provides an industrial process abnormal working condition prediction method based on operation parameter relevance change trend. In the single-parameter prediction stage, the method predicts each operation parameter through an exponential smoothing method according to the existing sensor data. In the correlation analysis stage, the invention calculates the correlation of the operation parameters by the known values of the operation parameters and the predicted values of the parameters, wherein the correlation of the parameters is represented by the similarity of a series of indexes representing the parameter curves. In the relevance trend prediction stage, the method constructs a multiple autoregressive model to predict the parameter relevance. The method provided by the invention takes the relevance of the operation parameters into consideration, and can obtain more complete equipment abnormal information and more advanced prediction results.
According to the invention concept, the invention provides an industrial process abnormal working condition prediction method based on operation parameter correlation analysis, which comprises the following steps:
step 1: carrying out fixed step length prediction on a measured value sequence acquired by each sensor in the industrial process by adopting a Holt exponential smoothing model;
step 2: constructing a triple characterized by the variation trend of the operation parameters of the industrial process in a sliding window mode, expressing a measured value sequence in a data window, and calculating the relevance of any two operation parameters in the window through a relevance index based on Euclidean distance;
and step 3: constructing a relevance prediction model according to the calculated relevance data, wherein the relevance prediction model is a multiple autoregressive model, and estimating model parameters through a partial least square algorithm;
and 4, step 4: and for the newly obtained data, predicting the abnormal working condition of the equipment according to the relevance prediction model, and continuously updating the relevance prediction model and the model parameters before predicting the abnormal working condition.
Based on the above scheme, the following implementation manner can be specifically adopted for each step:
preferably, step 1 is specifically as follows:
step 1.1: for industrial equipment with a plurality of sensors, the number of the sensors is recorded as N, when the equipment is in the operation process, operation parameter values, namely sensor data, representing the operation state of the equipment are continuously acquired, and measurement sequences of the sensors are recorded asWherein K represents the length of the sequence,represents the measured value of the sensor i at the kth sampling time point;
step 1.2: for the sensor i, a Holt exponential smoothing model is adopted to predict the measured value sequence of the sensor i, and the measured value is givenIts smoothed value can be calculated according to the following equation:
wherein,to representA smoothed value of (d);the model is a linear growth factor which represents the trend after smoothing, α and β are smoothing coefficients, the value ranges are (0,1), and the initial conditions of the Holt model are as follows:
step 1.3: after obtaining the smooth value and the linear growth factor of the measured data, the result is used for prediction and prediction valueComprises the following steps:
where l represents the prediction step size and τ is the sampling interval of the device sensor signal.
Preferably, step 2 comprises the following substeps:
step 2.1: the operation parameter correlation analysis comprises three stages of time sequence segmentation fitting, triple representation and correlation calculation; in the stage of time series piecewise fitting, for a data window with a fixed length L, it is recorded as a window WjOperating parameter XiI is 1,2, …, N, and L data in the window areWhere j is the start of the window and the corresponding sampling time is tj(ii) a If it is assumed that m data in one of the time slots in the windowCan be fitted by a line segment for the measured values thereinWhich corresponds to a value in the fitted line segment ofThe fitting error ERR of the line segment is calculated as:
when the window data with j as the starting point is subjected to piecewise linearization, the method is characterized in thatStarting to perform line segment fitting on the obtained object; the process of data segmentation fitting is carried out according to the following steps, wherein the set fitting error threshold is recorded as omegaE：
Step (1): set fitting starting point asFitting terminalIs characterized byh is 2; wherein, for the window W with j as the starting pointjThe initial fitting starting point is
Step (2): for dataAdopting a linear regression mode to carry out line segment fitting, thereby obtaining corresponding line segment dataCalculating the fitting error ERR according to the fitting error ERR calculation formula;
and (3): if ERR is less than or equal to omegaEIf so, changing h to h +1, and repeating the step (2); if ERR>ωEIf yes, saving the current fitting end point (namely the data segmentation point), resetting h to be 2, returning to the step 1, and fitting the next part of data by taking the current fitting end point as a new fitting starting point;
repeating the steps (1) to (3) until the window WjAll the data in the linear block are segmented, namely the segmented and linearized data is obtained
Step 2.2: in the segment triplet representation stage, a segment s is described in the form of the following tripletj，
Wherein k isjWhich represents the slope of the line segment,indicates the length of the line segment on the time axis, rjIndicating the rate of increase of the value of the line, i.e. for line dataFor window WjNew sequence after data line segmentationGet its three-tuple sequence representation as s1,s2,…,snN represents the number of line segments after data segmentation in the window;
step 2.3: in the correlation calculation stage, two operation parameters V in the equipment operation processAAnd an operating parameter VBFirstly, segmenting data after line segmentation: at the window WjIn, record parameter VAIs a segmentation point ofParameter VBIs a segmentation point ofWherein n isAAnd nBThe number of line segments after the measured data of the operating parameters A and B are segmented in the window is respectively represented; for parameter VAAnd VBThe segmentation points are merged, repeated items are removed, the segmentation points are arranged from small to large, and the segmentation point sequence of the two parameters is obtained asSubsequently, a parameter V is obtained from the sequence of segmentation points and the representation of tripletsAAnd VBThe new triplet sequence isAnd
to obtainWindow WjTwo internal parameters VAAnd VBAfter the triple sequence, the relevance index d based on Euclidean distance is usedABTo calculate two parameters V in the windowAAnd VBThe relevance of (A):
in the formula:as a parameter VAThe slope of the ith line segment,as a parameter VBThe slope of the ith line segment,as a parameter VAThe numerical growth rate of the ith line segment,as a parameter VBThe numerical growth rate of the ith line segment.
Preferably, step 3 comprises the following substeps:
step 3.1: constructing a relevance prediction model: setting the prediction step length as f, determining model design parameters U and M according to the known data length to enable U + f + M-1 to be equal to a window WjAnd constructing the following matrix:
wherein, { d1,d2,…,dU+f+M-1Denotes a parameter VAAnd VBA correlation sequence of djPresentation Window WjTwo internal parameters VAAnd VBCorrelation index d ofAB；
Step 3.2: for each relevance sequence, using the relevance values of U relevance values to predict the relevance value after f steps, constructing a relevance prediction model:
Fp＝Dpθ
wherein the parameter θ ═ θ1,θ2,…,θU]TAnd can be obtained by a partial least squares algorithm.
Preferably, step 4 comprises the following substeps:
in prediction, a new matrix is constructed using the newly acquired data from the sensors:
subsequently, the relevance values are f-step predicted using the constructed relevance prediction model, i.e.
When in useJudging the abnormality of the equipment, wherein dnormalIs a parameter V when the equipment is in a normal operation state at the initial operation stageAAnd VBCorrelation value of ωpA drift amount threshold for the correlation value relative to the initial normal; if it isReconstructing the model by using new data obtained by the sensor to update the model parameter theta;
and continuously predicting the relevance value of a certain prediction step length along with the updating of the data, thereby predicting the occurrence time of the abnormal working condition of the equipment.
The method for predicting the abnormal working condition of the industrial process based on the correlation analysis of the operation parameters can be used for a complex industrial system with a plurality of sensors. The method starts from the relevance among the operation parameters of the industrial process, and carries out the prediction of the abnormal working condition based on the relevance analysis of the operation parameters. In the single-parameter prediction stage, the method predicts each operation parameter through an exponential smoothing method according to the existing sensor data. In the correlation analysis stage, the invention calculates the correlation of the operation parameters by the known values of the operation parameters and the predicted values of the parameters, wherein the correlation of the parameters is represented by the similarity of a series of indexes representing the parameter curves. In the relevance trend prediction stage, the method constructs a multiple autoregressive model to predict the parameter relevance. The method provided by the invention takes the relevance of the operation parameters into consideration, and can obtain more complete equipment abnormal information and more advanced prediction results. The method provides powerful data support for subsequent equipment health management, is particularly valuable for high-reliability equipment maintenance management, and has wide prospects in the aspect of practical engineering application.
Drawings
FIG. 1 shows measured values and predicted results of turbine operating parameters;
FIG. 2 is a comparison of the predicted trend of correlation between turbine vacuum A and other parameters with the true value;
FIG. 3 shows the prediction result of the time of occurrence of an abnormality in a steam turbine.
Detailed Description
The present invention will now be further described with reference to the accompanying drawings, and some of the principles have been described in detail above, and will not be described again here. The following example illustrates the specific operation steps and verifies the effectiveness of the proposed method using a real case based on the turbine rough protection trip data.
The milling machine data records the operational degradation process of cutting metal material with the milling cutter. The initial working condition of the operation of the steam turbine is load 250MW and condenser vacuum 93kPa, the vacuum value of the condenser is used as an indication parameter, vacuum A starts to indicate abnormity from the 762 th sampling point, and when the vacuum value is reduced to 81kPa, the steam turbine trips. The method for predicting the abnormal working condition of the industrial process comprises the following steps:
step 1: and (3) predicting the measured value sequence acquired by each sensor in the industrial process by using a Holt exponential smoothing model in a fixed step length mode. The method specifically comprises the following substeps:
step 1.1: for industrial equipment with a plurality of sensors, the number of the sensors is recorded as N, when the equipment is in the operation process, operation parameter values, namely sensor data, representing the operation state of the equipment are continuously acquired, and measurement sequences of the sensors are recorded asWherein K represents the length of the sequence,represents the measured value of the sensor i at the kth sampling time point;
step 1.2: for the sensor i, a Holt exponential smoothing model is adopted to predict the measured value sequence of the sensor i, and the measured value is givenIts smoothed value can be calculated according to the following equation:
wherein,to representA smoothed value of (d);the model is a linear growth factor which represents the trend after smoothing, α and β are smoothing coefficients, the value ranges are (0,1), and the initial conditions of the Holt model are as follows:
step 1.3: after obtaining the smooth value and the linear growth factor of the measured data, the result is used for prediction and prediction valueComprises the following steps:
where l represents the prediction step size and τ is the sampling interval of the device sensor signal. In this example, the preset prediction step size is l-15.
According to step 1, for each operational parameter measurement sequence, a fixed-step prediction is performed, the result is given in fig. 1, and at the same time, an actual state monitoring measurement sequence is also given.
Step 2: and constructing a triple group characterized by the variation trend of the operation parameters of the industrial process in a sliding window mode, expressing a measured value sequence in a data window, and calculating the relevance of any two operation parameters in the window through a relevance index based on the Euclidean distance. The method specifically comprises the following substeps:
step 2.1: the operation parameter correlation analysis comprises three stages of time sequence segmentation fitting, triple representation and correlation calculation; in the stage of time series piecewise fitting, for a data window with a fixed length L, it is recorded as a window WjIn this example, the sliding window has a length of 100. Operating parameter XiI is 1,2, …, N, and L data in the window areWhere j is the start of the window and the corresponding sampling time is tj(ii) a If it is assumed that m data in one of the time slots in the windowIt is possible to fit exactly one line segment for the measured values thereinWhich corresponds to a value in the fitted line segment ofThe fitting error ERR of the line segment is calculated as:
when the window data with j as the starting point is subjected to piecewise linearization, the method is characterized in thatStarting to perform line segment fitting on the obtained object; the process of data segmentation fitting is carried out according to the following steps (wherein the set fitting error threshold is recorded as omega)EIn this example ωE＝0.025)：
Step (1): set fitting starting point asFitted endpoint ofh is 2; wherein, for the window W with j as the starting pointjThe initial fitting starting point is
Step (2): for dataAdopting a linear regression mode to carry out line segment fitting, thereby obtaining corresponding line segment dataCalculating the fitting error ERR according to the fitting error ERR calculation formula;
and (3): if ERR is less than or equal to omegaEIf so, changing h to h +1, and repeating the step (2); if ERR>ωEIf yes, saving the current fitting end point (namely the data segmentation point), resetting h to be 2, returning to the step 1, and fitting the next part of data by taking the current fitting end point as a new fitting starting point;
repeating the steps (1) to (3) until the window WjAll the data in the linear block are segmented, namely the segmented and linearized data is obtained
Step 2.2: in the segment triplet representation stage, a segment s is described in the form of the following tripletj，
Wherein k isjWhich represents the slope of the line segment,indicates the length of the line segment on the time axis, rjIndicating the rate of increase of the value of the line, i.e. for line dataFor window WjNew sequence after data line segmentationGet its three-tuple sequence representation as s1,s2,…,snN represents the number of line segments after data segmentation in the window;
step 2.3: in the correlation calculation stage, two operation parameters V in the equipment operation processAAnd an operating parameter VBFirstly, segmenting data after line segmentation: with an operating parameter VAAnd an operating parameter VBFor example, in the window WjIn, record parameter VAIs a segmentation point ofParameter VBIs a segmentation point ofWherein n isAAnd nBThe number of line segments after the measured data of the operating parameters A and B are segmented in the window is respectively represented; for parameter VAAnd VBIs segmented point ofMerging, removing repeated items, arranging from small to large to obtain a division point sequence of two parametersSubsequently, a parameter V is obtained from the sequence of segmentation points and the representation of tripletsAAnd VBThe new triplet sequence isAnd
obtaining a window WjTwo internal parameters VAAnd VBAfter the triple sequence, the relevance index d based on Euclidean distance is usedABTo calculate two parameters V in the windowAAnd VBThe relevance of (A):
in the formula:as a parameter VAThe slope of the ith line segment,as a parameter VBThe slope of the ith line segment,as a parameter VAThe numerical growth rate of the ith line segment,as a parameter VBThe numerical growth rate of the ith line segment.
And step 3: and constructing a relevance prediction model according to the calculated relevance data, wherein the relevance prediction model is a multiple autoregressive model, and estimating model parameters by a partial least square algorithm. The method specifically comprises the following substeps:
step 3.1: constructing a relevance prediction model: setting the prediction step length as f, determining model design parameters U and M according to the known data length to enable U + f + M-1 to be equal to a window WjAnd constructing the following matrix:
wherein, { d1,d2,…,dU+f+M-1Denotes a parameter VAAnd VBA correlation sequence of djPresentation Window WjTwo internal parameters VAAnd VBCorrelation index d ofAB；
Step 3.2: for each relevance sequence, using the relevance values of U relevance values to predict the relevance value after f steps, constructing a relevance prediction model:
Fp＝Dpθ
wherein the parameter θ ═ θ1,θ2,…,θU]TAnd can be obtained by a partial least squares algorithm.
And 4, step 4: and for the newly obtained data, predicting the abnormal working condition of the equipment according to the relevance prediction model, and continuously updating the relevance prediction model and the model parameters before predicting the abnormal working condition. The method specifically comprises the following substeps:
in prediction, a new matrix is constructed using the newly acquired data from the sensors:
subsequently, the relevance values are f-step predicted using the constructed relevance prediction model, i.e.
When in useJudging the abnormality of the equipment, wherein dnormalIs a parameter V when the equipment is in a normal operation state at the initial operation stageAAnd VBCorrelation value of ωpA drift amount threshold for the correlation value relative to the initial normal; if it isReconstructing the model by using new data obtained by the sensor to update the model parameter theta;
and continuously predicting the relevance value of a certain prediction step length along with the updating of the data, thereby predicting the occurrence time of the abnormal working condition of the equipment. In this example, the prediction step f is 10, and the threshold ω is setp＝0.1。
According to the steps 2 to 4, relevance calculation and relevance change trend prediction are carried out, and the result is shown in fig. 2. Subsequently, according to the set failure threshold ωpThe abnormal condition prediction is performed, and the result is shown in fig. 3, and meanwhile, the failure time obtained by detecting the relevance sequence obtained by using the actual data according to the set threshold is also given.
FIG. 1 shows measured values and predicted results of turbine operating parameters. As can be seen from the figure, the single parameter can be predicted well by using exponential smooth prediction. FIG. 2 shows the predicted trend of the correlation between turbine vacuum A and other parameters compared with the true value. As can be seen, a good prediction effect is obtained. Fig. 3 shows the prediction result of the turbine abnormality occurrence time. Fig. 3 shows the correlation between vacuum a and the remaining 6 parameters on the abscissa 1-6, respectively, and the ordinate indicates the predicted occurrence time of the abnormal condition. As can be seen from FIG. 3, a relatively accurate prediction result is obtained by using the prediction of the relevance of several sets of parameters, i.e., the occurrence of abnormal conditions can be predicted in advance by a certain step length. In addition, through the analysis of the raw data, the vacuum a is known to indicate an abnormality from the 762 th sampling point, and the decrease is accelerated, and through the analysis of the correlation change trend, as can be seen from fig. 3, the time for predicting the occurrence of the abnormal condition is the 590 th sampling point (ID6) at the earliest, that is, the occurrence of the abnormality is captured earlier. More importantly, through the relevance change trend, the parameter corresponding to the earliest abnormity generation is found, the judgment of the position of the abnormity generation is facilitated, and the important function of eliminating the abnormity is achieved.
Claims (5)
1. An industrial process abnormal working condition prediction method based on operation parameter correlation analysis is characterized by comprising the following steps:
step 1: carrying out fixed step length prediction on a measured value sequence acquired by each sensor in the industrial process by adopting a Holt exponential smoothing model;
step 2: constructing a triple characterized by the variation trend of the operation parameters of the industrial process in a sliding window mode, expressing a measured value sequence in a data window, and calculating the relevance of any two operation parameters in the window through a relevance index based on Euclidean distance;
and step 3: constructing a relevance prediction model according to the calculated relevance data, wherein the relevance prediction model is a multiple autoregressive model, and estimating model parameters through a partial least square algorithm;
and 4, step 4: and for the newly obtained data, predicting the abnormal working condition of the equipment according to the relevance prediction model, and continuously updating the relevance prediction model and the model parameters before predicting the abnormal working condition.
2. The method for predicting the abnormal working condition of the industrial process based on the correlation analysis of the operation parameters, as recited in claim 1, wherein: step 1 comprises the following substeps:
step 1.1: for industrial equipment with a plurality of sensors, the number of the sensors is recorded as N, when the equipment is in the operation process, operation parameter values, namely sensor data, representing the operation state of the equipment are continuously acquired, and measurement sequences of the sensors are recorded asWherein K represents the length of the sequence,represents the measured value of the sensor i at the kth sampling time point;
step 1.2: for the sensor i, a Holt exponential smoothing model is adopted to predict the measured value sequence of the sensor i, and the measured value is givenIts smoothed value can be calculated according to the following equation:
wherein,to representA smoothed value of (d);the model is a linear growth factor which represents the trend after smoothing, α and β are smoothing coefficients, the value ranges are (0,1), and the initial conditions of the Holt model are as follows:
step 1.3: after obtaining the smooth value and the linear growth factor of the measured data, the result is used for prediction and prediction valueComprises the following steps:
where l represents the prediction step size and τ is the sampling interval of the device sensor signal.
3. The method for predicting the abnormal working condition of the industrial process based on the correlation analysis of the operation parameters, as recited in claim 2, wherein: step 2 comprises the following substeps:
step 2.1: analysis of correlation of operating parameters is divided intoThree stages of time sequence segmentation fitting, triple representation and relevance calculation; in the stage of time series piecewise fitting, for a data window with a fixed length L, it is recorded as a window WjOperating parameter XiI is 1,2, …, N, and L data in the window areWhere j is the start of the window and the corresponding sampling time is tj(ii) a If it is assumed that m data in one of the time slots in the windowCan be fitted by a line segment for the measured values thereinWhich corresponds to a value in the fitted line segment ofThe fitting error ERR of the line segment is calculated as:
when the window data with j as the starting point is subjected to piecewise linearization, the method is characterized in thatStarting to perform line segment fitting on the obtained object; the process of data segmentation fitting is carried out according to the following steps, wherein the set fitting error threshold is recorded as omegaE：
Step (1): set fitting starting point asFitted endpoint ofWherein for j isWindow W of originjThe initial fitting starting point is
Step (2): for dataAdopting a linear regression mode to carry out line segment fitting, thereby obtaining corresponding line segment dataCalculating the fitting error ERR according to the fitting error ERR calculation formula;
and (3): if ERR is less than or equal to omegaEIf so, changing h to h +1, and repeating the step (2); if ERR>ωEIf yes, saving the current fitting end point, resetting h to be 2, returning to the step 1, and fitting the next part of data by taking the current fitting end point as a new fitting starting point;
repeating the steps (1) to (3) until the window WjAll the data in the linear block are segmented, namely the segmented and linearized data is obtained
Step 2.2: in the segment triplet representation stage, a segment s is described in the form of the following tripletj，
Wherein k isjWhich represents the slope of the line segment,indicates the length of the line segment on the time axis, rjIndicating the rate of increase of the value of the line, i.e. for line dataFor window WjNew sequence after data line segmentationGet its three-tuple sequence representation as s1,s2,…,snN represents the number of line segments after data segmentation in the window;
step 2.3: in the correlation calculation stage, two operation parameters V in the equipment operation processAAnd an operating parameter VBFirstly, segmenting data after line segmentation: at the window WjIn, record parameter VAIs a segmentation point ofParameter VBIs a segmentation point ofWherein n isAAnd nBThe number of line segments after the measured data of the operating parameters A and B are segmented in the window is respectively represented; for parameter VAAnd VBThe segmentation points are merged, repeated items are removed, the segmentation points are arranged from small to large, and the segmentation point sequence of the two parameters is obtained asSubsequently, a parameter V is obtained from the sequence of segmentation points and the representation of tripletsAAnd VBThe new triplet sequence isAnd
obtaining a window WjTwo internal parameters VAAnd VBAfter the triple sequence, the relevance index d based on Euclidean distance is usedABTo calculate two parameters V in the windowAAnd VBCorrelation of (2)：
In the formula:as a parameter VAThe slope of the ith line segment,as a parameter VBThe slope of the ith line segment,as a parameter VAThe numerical growth rate of the ith line segment,as a parameter VBThe numerical growth rate of the ith line segment.
4. The method for predicting the abnormal working condition of the industrial process based on the correlation analysis of the operation parameters, as recited in claim 3, wherein: step 3 comprises the following substeps:
step 3.1: constructing a relevance prediction model: setting the prediction step length as f, determining model design parameters U and M according to the known data length to enable U + f + M-1 to be equal to a window WjAnd constructing the following matrix:
wherein, { d1,d2,…,dU+f+M-1Denotes a parameter VAAnd VBA correlation sequence of djPresentation Window WjTwo internal parameters VAAnd VBCorrelation index d ofAB；
Step 3.2: for each relevance sequence, using the relevance values of U relevance values to predict the relevance value after f steps, constructing a relevance prediction model:
Fp＝Dpθ
wherein the parameter θ ═ θ1,θ2,…,θU]TAnd can be obtained by a partial least squares algorithm.
5. The method for predicting the abnormal working condition of the industrial process based on the correlation analysis of the operation parameters, as recited in claim 4, wherein: step 4 comprises the following substeps:
in prediction, a new matrix is constructed using the newly acquired data from the sensors:
subsequently, the relevance values are f-step predicted using the constructed relevance prediction model, i.e.
When in useJudging the abnormality of the equipment, wherein dnormalIs a parameter V when the equipment is in a normal operation state at the initial operation stageAAnd VBCorrelation value of ωpA drift amount threshold for the correlation value relative to the initial normal; if it isThe model is reconstructed using the new data obtained by the sensorsTo update the model parameters θ;
and continuously predicting the relevance value of a certain prediction step length along with the updating of the data, thereby predicting the occurrence time of the abnormal working condition of the equipment.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910244872.5A CN110008565B (en) | 2019-03-28 | 2019-03-28 | Industrial process abnormal working condition prediction method based on operation parameter correlation analysis |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910244872.5A CN110008565B (en) | 2019-03-28 | 2019-03-28 | Industrial process abnormal working condition prediction method based on operation parameter correlation analysis |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110008565A true CN110008565A (en) | 2019-07-12 |
| CN110008565B CN110008565B (en) | 2020-12-11 |
Family
ID=67168700
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910244872.5A Active CN110008565B (en) | 2019-03-28 | 2019-03-28 | Industrial process abnormal working condition prediction method based on operation parameter correlation analysis |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110008565B (en) |
Cited By (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111553606A (en) * | 2020-05-06 | 2020-08-18 | 上海核工程研究设计院有限公司 | Evaluation method for nuclear power equipment component design and equipment index relevance |
| CN111913444A (en) * | 2019-09-07 | 2020-11-10 | 宁波大学 | Chemical process monitoring method based on time sequence multi-block modeling strategy |
| CN111983478A (en) * | 2020-07-07 | 2020-11-24 | 江苏大学 | Electrochemical energy storage power station SOC anomaly detection method based on Holt linear trend model |
| CN112464146A (en) * | 2020-12-03 | 2021-03-09 | 北京航空航天大学 | Key subsystem based on historical telemetering data and single-machine correlation health baseline construction method |
| CN112801426A (en) * | 2021-04-06 | 2021-05-14 | 浙江浙能技术研究院有限公司 | Industrial process fault fusion prediction method based on correlation parameter mining |
| CN113624025A (en) * | 2021-08-18 | 2021-11-09 | 浙江大学 | Condenser vacuum low-jump machine symptom capturing method based on correlation of operation parameters |
| CN113792988A (en) * | 2021-08-24 | 2021-12-14 | 河北先河环保科技股份有限公司 | Online monitoring data anomaly identification method for enterprise |
| CN113867262A (en) * | 2020-06-30 | 2021-12-31 | 华晨宝马汽车有限公司 | Apparatus, method and medium for monitoring operation state of spindle device |
| CN115034094A (en) * | 2022-08-10 | 2022-09-09 | 南通恒强轧辊有限公司 | Prediction method and system for operation state of metal processing machine tool |
| CN115509187A (en) * | 2022-09-20 | 2022-12-23 | 北京中佳瑞通科技有限公司 | Industrial big data processing method and system |
| CN115689071A (en) * | 2023-01-03 | 2023-02-03 | 南京工大金泓能源科技有限公司 | Equipment fault fusion prediction method and system based on correlation parameter mining |
| CN115829422A (en) * | 2023-02-21 | 2023-03-21 | 创银科技(南通)有限公司 | Industrial equipment operation abnormal state identification method based on big data |
| WO2023094916A1 (en) * | 2021-11-23 | 2023-06-01 | International Business Machines Corporation | Identifying persistent anomalies for failure prediction |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20090265155A1 (en) * | 2008-04-21 | 2009-10-22 | Nec Electronics Corporation | Method of predicting reliability of semiconductor device, reliability prediction system using the same and storage medium storing program causing computer to execute the same |
| CN108037387A (en) * | 2017-12-05 | 2018-05-15 | 北京能源集团有限责任公司 | The equipment fault analysis method and device collected based on cluster |
| CN108873859A (en) * | 2018-05-31 | 2018-11-23 | 浙江工业大学 | Based on the bridge-type grab ship unloader fault prediction model method for improving correlation rule |
| CN109086804A (en) * | 2018-07-12 | 2018-12-25 | 中石化石油机械股份有限公司 | A kind of hydraulic device fault forecast method merged based on multi source status monitoring information and reliability characteristic |
| CN109110608A (en) * | 2018-10-25 | 2019-01-01 | 歌拉瑞电梯股份有限公司 | A kind of elevator faults prediction technique based on big data study |
-
2019
- 2019-03-28 CN CN201910244872.5A patent/CN110008565B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20090265155A1 (en) * | 2008-04-21 | 2009-10-22 | Nec Electronics Corporation | Method of predicting reliability of semiconductor device, reliability prediction system using the same and storage medium storing program causing computer to execute the same |
| CN108037387A (en) * | 2017-12-05 | 2018-05-15 | 北京能源集团有限责任公司 | The equipment fault analysis method and device collected based on cluster |
| CN108873859A (en) * | 2018-05-31 | 2018-11-23 | 浙江工业大学 | Based on the bridge-type grab ship unloader fault prediction model method for improving correlation rule |
| CN109086804A (en) * | 2018-07-12 | 2018-12-25 | 中石化石油机械股份有限公司 | A kind of hydraulic device fault forecast method merged based on multi source status monitoring information and reliability characteristic |
| CN109110608A (en) * | 2018-10-25 | 2019-01-01 | 歌拉瑞电梯股份有限公司 | A kind of elevator faults prediction technique based on big data study |
Non-Patent Citations (3)
| Title |
|---|
| CHUNHUI ZHAO 等: "Critical-to-Fault-Degradation Variable Analysis and Direction Extraction for Online Fault Prognostic", 《IEEE TRANSACTIONS ON CONTROL SYSTEMS TECHNOLOGY》 * |
| 卢会国 等: "基于关联规则的故障诊断方法研究", 《成都信息工程学院学报》 * |
| 方名菊 等: "基于形式概念的汽轮机相似故障案例检索", 《山东科技大学学报(自然科学版)》 * |
Cited By (17)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111913444A (en) * | 2019-09-07 | 2020-11-10 | 宁波大学 | Chemical process monitoring method based on time sequence multi-block modeling strategy |
| CN111913444B (en) * | 2019-09-07 | 2022-03-18 | 宁波大学 | Chemical process monitoring method based on time sequence multi-block modeling strategy |
| CN111553606A (en) * | 2020-05-06 | 2020-08-18 | 上海核工程研究设计院有限公司 | Evaluation method for nuclear power equipment component design and equipment index relevance |
| CN113867262A (en) * | 2020-06-30 | 2021-12-31 | 华晨宝马汽车有限公司 | Apparatus, method and medium for monitoring operation state of spindle device |
| CN111983478A (en) * | 2020-07-07 | 2020-11-24 | 江苏大学 | Electrochemical energy storage power station SOC anomaly detection method based on Holt linear trend model |
| CN112464146A (en) * | 2020-12-03 | 2021-03-09 | 北京航空航天大学 | Key subsystem based on historical telemetering data and single-machine correlation health baseline construction method |
| CN112801426A (en) * | 2021-04-06 | 2021-05-14 | 浙江浙能技术研究院有限公司 | Industrial process fault fusion prediction method based on correlation parameter mining |
| CN113624025A (en) * | 2021-08-18 | 2021-11-09 | 浙江大学 | Condenser vacuum low-jump machine symptom capturing method based on correlation of operation parameters |
| CN113792988A (en) * | 2021-08-24 | 2021-12-14 | 河北先河环保科技股份有限公司 | Online monitoring data anomaly identification method for enterprise |
| WO2023094916A1 (en) * | 2021-11-23 | 2023-06-01 | International Business Machines Corporation | Identifying persistent anomalies for failure prediction |
| CN115034094A (en) * | 2022-08-10 | 2022-09-09 | 南通恒强轧辊有限公司 | Prediction method and system for operation state of metal processing machine tool |
| CN115509187A (en) * | 2022-09-20 | 2022-12-23 | 北京中佳瑞通科技有限公司 | Industrial big data processing method and system |
| CN115509187B (en) * | 2022-09-20 | 2023-04-18 | 北京中佳瑞通科技有限公司 | Industrial big data processing method and system |
| CN115689071A (en) * | 2023-01-03 | 2023-02-03 | 南京工大金泓能源科技有限公司 | Equipment fault fusion prediction method and system based on correlation parameter mining |
| CN115689071B (en) * | 2023-01-03 | 2023-05-02 | 南京工大金泓能源科技有限公司 | Equipment fault fusion prediction method and system based on associated parameter mining |
| CN115829422A (en) * | 2023-02-21 | 2023-03-21 | 创银科技(南通)有限公司 | Industrial equipment operation abnormal state identification method based on big data |
| CN115829422B (en) * | 2023-02-21 | 2024-01-02 | 北京瀚海蓝山智能科技有限公司 | Industrial equipment operation abnormal state identification method based on big data |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110008565B (en) | 2020-12-11 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110008565B (en) | Industrial process abnormal working condition prediction method based on operation parameter correlation analysis | |
| Ahmad et al. | A hybrid prognostics technique for rolling element bearings using adaptive predictive models | |
| CN110008253B (en) | Industrial data association rule mining and abnormal working condition prediction method | |
| CA2634328C (en) | Method and system for trend detection and analysis | |
| CN110414154B (en) | Fan component temperature abnormity detection and alarm method with double measuring points | |
| KR102067344B1 (en) | Apparatus and Method for Detecting Abnormal Vibration Data | |
| CN108038044A (en) | A kind of method for detecting abnormality towards continuous monitored target | |
| KR102119661B1 (en) | A method to predict health index transition and residual life for turbomachinery | |
| CN109960232B (en) | Method for selecting leading auxiliary parameter and method for equipment maintenance pre-diagnosis | |
| JP5827426B1 (en) | Predictive diagnosis system and predictive diagnosis method | |
| CN117972614B (en) | Real-time temperature monitoring method in charging and discharging process of mobile power supply | |
| US7949497B2 (en) | Machine condition monitoring using discontinuity detection | |
| CN112700039A (en) | Steady state detection and extraction method for load operation data of thermal power plant | |
| CN117471346A (en) | Method and system for determining remaining life and health status of retired battery module | |
| CN117961643A (en) | Numerical control machine tool online fault monitoring method based on multi-axis current signals | |
| CN118309644A (en) | Pipeline pump operation flow monitoring method and system based on digital twin | |
| CN117150244B (en) | Intelligent power distribution cabinet state monitoring method and system based on electrical parameter analysis | |
| CN111913463B (en) | State monitoring method for chemical volume control system of nuclear power plant | |
| JP2002090266A (en) | Remaining life-predicting device | |
| CN116255555A (en) | Lubricating oil on-line monitoring method and system based on viscosity analysis | |
| EP3869286A1 (en) | System and method for detecting turbine underperformance and operation anomaly | |
| JP7441811B2 (en) | diagnostic equipment | |
| CN118525992B (en) | Perfuming and feeding fault diagnosis method based on time sequence signal analysis | |
| WO2022071189A1 (en) | State-monitoring device and state-monitoring method | |
| KR20240146740A (en) | Predictive diagnosis method and system of nuclear power plant equipment using rule base model |
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 |