CN110516323A - A kind of pressure pipeline damage forecast method based on Time-Series analysis - Google Patents

A kind of pressure pipeline damage forecast method based on Time-Series analysis Download PDF

Info

Publication number
CN110516323A
CN110516323A CN201910729109.1A CN201910729109A CN110516323A CN 110516323 A CN110516323 A CN 110516323A CN 201910729109 A CN201910729109 A CN 201910729109A CN 110516323 A CN110516323 A CN 110516323A
Authority
CN
China
Prior art keywords
formula
stress
days
following
sequence
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
CN201910729109.1A
Other languages
Chinese (zh)
Other versions
CN110516323B (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.)
Hefei General Machinery Research Institute Special Equipment Inspection Station Co Ltd
Hefei General Machinery Research Institute Co Ltd
Original Assignee
Hefei General Machinery Research Institute Special Equipment Inspection Station Co Ltd
Hefei General Machinery Research Institute Co Ltd
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 Hefei General Machinery Research Institute Special Equipment Inspection Station Co Ltd, Hefei General Machinery Research Institute Co Ltd filed Critical Hefei General Machinery Research Institute Special Equipment Inspection Station Co Ltd
Priority to CN201910729109.1A priority Critical patent/CN110516323B/en
Publication of CN110516323A publication Critical patent/CN110516323A/en
Application granted granted Critical
Publication of CN110516323B publication Critical patent/CN110516323B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

The pressure pipeline damage forecast method based on Time-Series analysis that the present invention relates to a kind of.The present invention calculates the first-order difference of monitor stress value first, by extracting stationarity information, the stress value difference in the following several days periods is predicted in conjunction with autoregressive moving average Time-Series analysis arma modeling, and then calculate the stress value in the following several days periods, finally by the development trend of the method for weighted moving average (EWMA) analysis predicted stresses value, to carry out damage alarming.The present invention may be implemented to carry out real-time early warning in in-service pipeline, have the characteristics that early warning is timely, accurate, this provides technical support for the safe operation of pipeline long period.

Description

A kind of pressure pipeline damage forecast method based on Time-Series analysis
Technical field
The invention belongs to process industries to be safely operated safeguards technique field, be specifically related to a kind of based on Time-Series analysis Pressure pipeline damage forecast method.
Background technique
Pressure pipeline is the important component in process industry, and safety operation level directly determines petroleum chemical enterprise Safety operation level.With the expansion of petroleum chemical enterprise's scale, the limitation of conventional detection means is also increasingly apparent, such as: day Often consumption manpower is too many, and specific position has certain risk, is easy the presence of detection dead angle;It can only for equipment safety evaluation Based on the test data under shutdown status, it can not reflect the truth under service state, it can not be to in-service device structure state Be measured in real time, twice detect during whether safety can only be estimated according to testing result before, but it is this estimate usually by The influence of many uncertain factors in pipeline performance environment and there is very big error.
In order to ensure its safety operation level, domestic and international researcher proposes to carry out real-time monitoring to pipe stress state, Pipeline damages when thinking monitor value beyond stress threshold.In practical engineering application, each position lesion development trend is different, Once beyond damaging after stress threshold, trend development is too fast to be often led to have little time to prepare emergency preplan.Open source information shows phase It closes researcher and the following stress of monitoring data prediction is based on by autoregressive moving average Time Series Analysis Method (ARMA (p, q)) Value, but the fluctuation of stress value stationarity caused by pipeline operating condition changes, will lead to predicted value distortion.
Summary of the invention
In order to solve the above technical problem, the present invention provides a kind of pressure pipeline damage forecast side based on Time-Series analysis Method.This method provides technical support for the safe operation of pipeline long period.
In order to achieve the object of the present invention, the invention adopts the following technical scheme:
A kind of pressure pipeline damage forecast method based on Time-Series analysis, comprising the following steps:
Step 1, pressure pipeline stress value is obtained by sensor unit, the stress value frequency acquisition is f times/day;
Step 2, the stress value sequence accumulated in time history M days is denoted as:
t}={ τ12,...,τM×f} (1)
Assuming that steady sequency subject to the stress value sequence, i.e., the described stress value is fluctuated around a certain steady state value small range, Its fluctuation amplitude meets normal distribution;
Step 3, the stress value of the following N days pressure pipelines is obtained by Time-Series analysis;
Step 3-1 extracts stationarity information based on calculus of finite differences;
It selects first-order difference to extract the stationarity information of monitor stress value trend, it is poor to obtain the steady random stress of zero-mean Sub-sequence, as follows:
{Δτt}={ Δ τ2,Δτ3,…ΔτM×f} (2)
In formula (2), Δ τttt-1, wherein t=2,3 ..., M × f;
Step 3-2, stress difference sequence { Δ τtAutoregressive moving average Time-Series analysis ARMA (p, q) model parameter estimates Meter;
Stress difference sequence { Δ τtInternal data relationship is expressed as with ARMA (p, q) model:
In formula (3) formulaReferred to as Parameters of Autoregressive Models;θ1, θ2..., θqReferred to as sliding average mould Shape parameter;
A in formula (3)t-1, at-2..., at-qReferred to as preceding q walks distracter, and mutually indepedent and obedience mean value is 0, at-1, at-2..., at-qMeeting variance isNormal distribution, be denoted as:
I=t-q in formula (4) ..., t-2, t-1, t
Step 3-3 determines autoregressive moving average Time-Series analysis ARMA (p, q) model order;
To formula (3) two sides respectively multiplied by Δ τt-kAnd mathematic expectaion is sought, it obtains:
In formula: Rk=E [Δ τtΔτt-k], Ri-k=E [Δ τt-iΔτt-k],And:
Because of RkFor even function, k=q+1, q+2 ..., q+p are taken so working as, using formula (6) property, after formula (5) expansion It obtains:
Following matrix is obtained according to formula (7):
Parameters of Autoregressive Models can be acquired by formula (8)
It brings the Parameters of Autoregressive Models into formula (5), while taking k=0,1,2 ..., q, and obtained using formula (6) property To following equation:
In formula (9):
By the way that moving average model parameter θ can be found out in formula (9)1, θ2..., θq
Step 3-4 determines autoregressive moving average Time-Series analysis ARMA (p, q) model order;
Using information criterion abstraction sequence { Δ τtIn maximum fault information, information criterion function are as follows:
P, q < 5 are selected for pipeline configuration, calculate p, all combined values of q take so that AIC (p, q) value the smallest p, q Group is combined into ARMA (p, q) model order;
Step 3-5 calculates the following N days stress values:
After determining ARMA (p, q) model order, the following N days stress values can be calculated, steps are as follows:
It calculates the following N days stress difference and takes t=M × f+1, M × f+2 ... using formula (3), (M+1) × f obtains future F stress difference is as follows within 1 day:
Using same method, then the available following N days stress difference, convolution (2) construct stress difference sequence
Further acquire stress sequence
M≤100, N≤20 are generally taken in order to ensure solving precision for pipeline configuration;
Step 4, stress development trend is analyzed by the method for weighted moving average (EWMA):
The expectation μ and variances sigma of preceding M days stress values2It is indicated with maximal possibility estimation are as follows:
To stress sequence { τtWeighted moving average analysis is done, the control variable of i-th of stress value is expressed as in sequence:
In formula (17), Z is taken0=μ, λ are smoothing parameter, and 0 < λ < 1, for pipeline architecture, general straight tube take λ= 0.6;
According to EWMA control figure statistics distribution feature, variable Z can be controlled in the hope of its statistical natureiExpectation:
In formula (20):
Control variable ZiVariance are as follows:
According to preceding M days stress values, the control for acquiring EWMA control figure is limited to:
Wherein UCL is upper control limit, and LCL is lower control limit;
Step 5, damage alarming:
The control variable Z of stress value in N days following is calculated according to formula (17)iThe development trend of value, as control variable ZiValue It is gradually deviated from and is more than that control is prescribed a time limit, then it is assumed that will damage within following N days;When development trend is in above-mentioned control limit, it is believed that There is no updating Stressing history data sequence at this time for damage are as follows:
t}={ τf+1,...,τM×f,...,τ(M+1)×f-1(M+1)×f} (23)
Step 3,4 are repeated, the control variable Z to stress value in following N days is continuediDevelopment trend is judged.
The beneficial effects of the present invention are:
The present invention calculates the first-order difference of monitor stress value first, by extracting stationarity information, slides in conjunction with autoregression Average Time-Series analysis arma modeling predicts the stress value difference in the following several days periods, and then calculates the following several days weeks Stress value in phase, finally by the development trend of the method for weighted moving average (EWMA) analysis predicted stresses value, to be damaged Early warning.The present invention may be implemented to carry out real-time early warning in in-service pipeline, have the characteristics that early warning is timely, accurate, this is pipeline Long period safe operation provides technical support.
Detailed description of the invention
Fig. 1 is preceding 70 days stress time-history curves.
Fig. 2 is that preceding 70 days stress differences divide time-history curves.
Fig. 3 is the difference time-history curves of preceding 70 days stress and following 15 days predicted stresses.
Fig. 4 is the time-history curves of preceding 70 days stress and following 15 days predicted stresses.
Fig. 5 is the control variable trends curve of preceding 70 days stress and following 15 days predicted stresses.
Fig. 6 is the control variable trends curve of following 15 days predicted stresses after 1025 days.
Specific embodiment
More specific detail is made to technical solution of the present invention below with reference to embodiment:
A kind of pressure pipeline damage forecast method based on Time-Series analysis, comprising the following steps:
Step 1, pressure pipeline stress value being obtained by sensor unit, the stress value frequency acquisition is f=2 times/ It;
Step 2, the stress value sequence accumulated in time history 70 days is denoted as:
t}={ τ12,…τ140} (1)
Above-mentioned stress value sequence time history curve is as shown in Figure 1;
Assuming that steady sequency subject to the stress value sequence, i.e., the described stress value is fluctuated around a certain steady state value small range, Its fluctuation amplitude meets normal distribution, and above-mentioned stress value sequence time history curve is as shown in Figure 1;
Step 3, the stress value of following 15 days pressure pipelines is obtained by Time-Series analysis;
Step 3-1 extracts stationarity information based on calculus of finite differences;
It selects first-order difference to extract the stationarity information of monitor stress value trend, it is poor to obtain the steady random stress of zero-mean Sub-sequence, as follows:
{Δτt}={ Δ τ2,Δτ3,…,Δτ140} (2)
In formula (2), Δ τttt-1, wherein t=2,3 ..., 140;
Above-mentioned stress Differential time course curve is as shown in Figure 2;
Step 3-2, stress difference sequence { Δ τtAutoregressive moving average Time-Series analysis ARMA (p, q) model parameter estimates Meter;
Stress difference sequence { Δ τtInternal data relationship is expressed as with ARMA (p, q) model:
In formula (3) formulaReferred to as Parameters of Autoregressive Models;θ1, θ2..., θqReferred to as sliding average mould Shape parameter;
A in formula (3)t-1, at-2..., at-qReferred to as preceding q walks distracter, and mutually indepedent and obedience mean value is 0, at-1, at-2..., at-qMeeting variance isNormal distribution, be denoted as:
I=t-q in formula (4) ..., t-2, t-1, t
Step 3-3 determines autoregressive moving average Time-Series analysis ARMA (p, q) model order;
To formula (3) two sides respectively multiplied by Δ τt-kAnd mathematic expectaion is sought, it obtains:
In formula: Rk=E [Δ τtΔτt-k], Ri-k=E [Δ τt-iΔτt-k],And:
Because of RkFor even function, reason is as follows:
Rk=E [Δ τtΔτt-k]=E [Δ τt-kΔτt]=E [Δ τt-kΔτt-k+k]
=E [Δ τsΔτs-(-k)]=E [Δ τtΔτt-(-k)]=R-k
So, using formula (6) property, will be obtained after formula (5) expansion when taking k=q+1, q+2 ..., q+p:
Following matrix is obtained according to formula (7):
Parameters of Autoregressive Models can be acquired by formula (8)
It brings the Parameters of Autoregressive Models into formula (5), while taking k=0,1,2 ..., q, and obtained using formula (6) property To following equation:
That is:
In formula (9):
By the way that moving average model parameter θ can be found out in formula (9)1, θ2..., θqAnd σa
Step 3-4 determines autoregressive moving average Time-Series analysis ARMA (p, q) model order;
Using information criterion abstraction sequence { Δ τtIn maximum fault information, information criterion function are as follows:
P, q < 5 are selected for pipeline configuration, calculate p, all combined values of q take so that AIC (p, q) value the smallest p, q Group is combined into ARMA (p, q) model order;
It acquires
Table 1 difference p, q combine corresponding AIC (p, q) value
p q AIC(p,q)
0 1 3.37
0 2 3.11
3 2 1.95
4 3 2.67
4 4 2.74
Determine ARMA (p, q) model order, p=3, q=2;At this point, Parameters of Autoregressive Models are as follows: σa=0.12;Moving average model parameter θ1=0.53, θ2=-0.79;
Step 3-5 calculates the following 1-15 days stress values:
After determining ARMA (p, q) model order, the following N days stress values can be calculated, steps are as follows:
Following 1 day stress difference is calculated, using formula (3), takes t=141,142 obtain following 1 day 2 stress difference such as Under:
A in formula141, a142To obey NID=(0, σa) random number;Using same method, t=142 is taken, 143 ..., 168,169,170 available following 15 days stress difference values, convolution (2) construct stress difference sequence:
Above-mentioned stress difference sequence time history curve is as shown in Figure 3;
Further acquire stress sequence
Above-mentioned stress value sequence time course curve is as shown in Figure 4;
M≤100, N≤20 are generally taken in order to ensure solving precision for pipeline configuration;
Step 4, stress development trend is analyzed by the method for weighted moving average (EWMA):
The expectation μ and variances sigma of preceding 70 days stress values2It is indicated with maximal possibility estimation are as follows:
To stress sequence { τtWeighted moving average analysis is done, the control variable of i-th of stress value is expressed as in sequence:
In formula (17), Z is taken0=μ, λ are smoothing parameter, and 0 < λ < 1, for pipeline architecture, general straight tube take λ= 0.6;
According to EWMA control figure statistics distribution feature, variable Z can be controlled in the hope of its statistical natureiExpectation:
In formula (20):
Control variable ZiVariance are as follows:
According to preceding M days stress values, the control for acquiring EWMA control figure is limited to:
Wherein UCL is upper control limit, and LCL is lower control limit;
Step 5, damage alarming:
The control variable Z of stress value in N days following is calculated according to formula (17)iThe development trend of value;
As shown in figure 5, providing control limit and control variable change situation, it can be seen that control variable is controlling within first 70 days In limit, 15 days futures can be seen that according to 15 days predicted values of future, in control limit, structure will not damage control variable Wound updates data sequence of the stress value in the 2nd~71 day time history at this time are as follows:
t}={ τ34,…,τ140141142} (23)
Step 3,4 are repeated, the control variable Z to stress value in following N days is continuediDevelopment trend is judged.Work as acquisition After 1025th day data, as shown in Figure 6, it is predicted that following 15 days control variables will exceed control limit, it is believed that following 15 days pipelines It will damage.

Claims (1)

1. a kind of pressure pipeline damage forecast method based on Time-Series analysis, it is characterised in that the following steps are included:
Step 1, pressure pipeline stress value is obtained by sensor unit, the stress value frequency acquisition is f times/day;
Step 2, the stress value sequence accumulated in time history M days is denoted as:
t}={ τ12,...,τM×f} (1)
Assuming that steady sequency subject to the stress value sequence, i.e., the described stress value is fluctuated around a certain steady state value small range, wave Dynamic amplitude meets normal distribution;
Step 3, the stress value of the following N days pressure pipelines is obtained by Time-Series analysis;
Step 3-1 extracts stationarity information based on calculus of finite differences;
It selects first-order difference to extract the stationarity information of monitor stress value trend, obtains the steady random stress difference sequence of zero-mean Column, as follows:
{Δτt}={ Δ τ2,Δτ3,…ΔτM×f} (2)
In formula (2), Δ τttt-1, wherein t=2,3 ..., M × f;
Step 3-2, stress difference sequence { Δ τtAutoregressive moving average Time-Series analysis ARMA (p, q) model parameter estimation;
Stress difference sequence { Δ τtInternal data relationship is expressed as with ARMA (p, q) model:
In formula (3) formulaReferred to as Parameters of Autoregressive Models;θ1, θ2..., θqReferred to as moving average model is joined Number;
A in formula (3)t-1, at-2..., at-qReferred to as preceding q walks distracter, and mutually indepedent and obedience mean value is 0, at-1, at-2..., at-qMeeting variance isNormal distribution, be denoted as:
I=t-q in formula (4) ..., t-2, t-1, t
Step 3-3 determines autoregressive moving average Time-Series analysis ARMA (p, q) model order;
To formula (3) two sides respectively multiplied by Δ τt-kAnd mathematic expectaion is sought, it obtains:
In formula: Rk=E [Δ τtΔτt-k], Ri-k=E [Δ τt-iΔτt-k],And:
Because of RkFor even function, k=q+1, q+2 ..., q+p are taken so working as, using formula (6) property, will be obtained after formula (5) expansion:
Following matrix is obtained according to formula (7):
Parameters of Autoregressive Models can be acquired by formula (8)
It brings the Parameters of Autoregressive Models into formula (5), while taking k=0,1,2 ..., q, and it is following to utilize formula (6) property to obtain Equation:
In formula (9):
By the way that moving average model parameter θ can be found out in formula (9)1, θ2..., θq
Step 3-4 determines autoregressive moving average Time-Series analysis ARMA (p, q) model order;
Using information criterion abstraction sequence { Δ τtIn maximum fault information, information criterion function are as follows:
P, q < 5 are selected for pipeline configuration, calculate p, all combined values of q take so that AIC (p, q) value the smallest p, q combination For ARMA (p, q) model order;
Step 3-5 calculates the following N days stress values:
After determining ARMA (p, q) model order, the following N days stress values can be calculated, steps are as follows:
It calculates the following N days stress difference and takes t=M × f+1, M × f+2 ... using formula (3), (M+1) × f obtains 1 day f of future A stress difference is as follows:
Using same method, then the available following N days stress difference, convolution (2) construct stress difference sequence:
Further acquire stress sequence:
M≤100, N≤20 are generally taken in order to ensure solving precision for pipeline configuration;
Step 4, stress development trend is analyzed by the method for weighted moving average (EWMA):
The expectation μ and variances sigma of preceding M days stress values2It is indicated with maximal possibility estimation are as follows:
To stress sequence { τtWeighted moving average analysis is done, the control variable of i-th of stress value is expressed as in sequence:
In formula (17), Z is taken0=μ, λ are smoothing parameter, and 0 < λ < 1, and for pipeline architecture, general straight tube takes λ=0.6;
According to EWMA control figure statistics distribution feature, variable Z can be controlled in the hope of its statistical natureiExpectation:
In formula (20):
Control variable ZiVariance are as follows:
According to preceding M days stress values, the control for acquiring EWMA control figure is limited to:
Wherein UCL is upper control limit, and LCL is lower control limit;
Step 5, damage alarming:
The control variable Z of stress value in N days following is calculated according to formula (17)iThe development trend of value, as control variable ZiIt is worth gradually inclined It prescribes a time limit from and more than control, then it is assumed that will damage within following N days;When development trend is in above-mentioned control limit, it is believed that damage does not have There is generation, update Stressing history data sequence at this time are as follows:
t}={ τf+1,...,τM×f,...,τ(M+1)×f-1(M+1)×f} (23)
Step 3,4 are repeated, the control variable Z to stress value in following N days is continuediDevelopment trend is judged.
CN201910729109.1A 2019-08-08 2019-08-08 Pressure pipeline damage prediction method based on time sequence analysis Active CN110516323B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910729109.1A CN110516323B (en) 2019-08-08 2019-08-08 Pressure pipeline damage prediction method based on time sequence analysis

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910729109.1A CN110516323B (en) 2019-08-08 2019-08-08 Pressure pipeline damage prediction method based on time sequence analysis

Publications (2)

Publication Number Publication Date
CN110516323A true CN110516323A (en) 2019-11-29
CN110516323B CN110516323B (en) 2023-03-24

Family

ID=68623933

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910729109.1A Active CN110516323B (en) 2019-08-08 2019-08-08 Pressure pipeline damage prediction method based on time sequence analysis

Country Status (1)

Country Link
CN (1) CN110516323B (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113436733A (en) * 2021-08-26 2021-09-24 肾泰网健康科技(南京)有限公司 Characteristic construction method of hemodialysis scheme generation model based on fusion experience

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102042848A (en) * 2010-11-23 2011-05-04 北京航空航天大学 Prediction method of multi-functional parameter accelerated degradation testing product life based on multivariate hybrid time sequence analysis
CN102650527A (en) * 2012-05-25 2012-08-29 北京航空航天大学 Temperature compensation method for denoising fiber-optic gyroscope on basis of time series analysis
CN106059661A (en) * 2015-12-25 2016-10-26 国家电网公司 Time sequence analysis based optical transmission network trend prediction method
CN106529145A (en) * 2016-10-27 2017-03-22 浙江工业大学 Bridge monitoring data prediction method based on ARIMA-BP neural network
JP2018059464A (en) * 2016-10-06 2018-04-12 いすゞ自動車株式会社 Monitoring system of damage state caused by temperature change and monitoring method of damage state caused by temperature change

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102042848A (en) * 2010-11-23 2011-05-04 北京航空航天大学 Prediction method of multi-functional parameter accelerated degradation testing product life based on multivariate hybrid time sequence analysis
CN102650527A (en) * 2012-05-25 2012-08-29 北京航空航天大学 Temperature compensation method for denoising fiber-optic gyroscope on basis of time series analysis
CN106059661A (en) * 2015-12-25 2016-10-26 国家电网公司 Time sequence analysis based optical transmission network trend prediction method
JP2018059464A (en) * 2016-10-06 2018-04-12 いすゞ自動車株式会社 Monitoring system of damage state caused by temperature change and monitoring method of damage state caused by temperature change
CN106529145A (en) * 2016-10-27 2017-03-22 浙江工业大学 Bridge monitoring data prediction method based on ARIMA-BP neural network

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
肖汉斌等: "基于时序分析的桥式起重机结构损伤预警研究", 《武汉理工大学学报》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113436733A (en) * 2021-08-26 2021-09-24 肾泰网健康科技(南京)有限公司 Characteristic construction method of hemodialysis scheme generation model based on fusion experience
CN113436733B (en) * 2021-08-26 2021-11-30 肾泰网健康科技(南京)有限公司 Characteristic construction method of hemodialysis scheme generation model based on fusion experience

Also Published As

Publication number Publication date
CN110516323B (en) 2023-03-24

Similar Documents

Publication Publication Date Title
CN108681633B (en) Condensate pump fault early warning method based on state parameters
CN106873571B (en) early warning method based on data and model fusion
US20230023931A1 (en) Hydraulic turbine cavitation acoustic signal identification method based on big data machine learning
KR100867938B1 (en) Prediction method for watching performance of power plant measuring instrument by dependent variable similarity and kernel feedback
CN108804740B (en) Long-distance pipeline pressure monitoring method based on integrated improved ICA-KRR algorithm
CN110533294B (en) Nuclear power plant operation fault alarm method based on artificial intelligence technology
CN105404280A (en) Industrial process fault detection method based on autoregression dynamic hidden variable model
CN103631145B (en) Multiple operating modes process monitoring method and system based on monitor control index switching
CN106483847B (en) A kind of water cooler fault detection method based on adaptive ICA
CN111354496B (en) Nuclear power plant accident online diagnosis and state tracking prediction method
CA2723695A1 (en) Device and method to monitor a gas turbine
CN116938676A (en) Communication risk combined early warning method based on data source analysis
CN104268649A (en) Water pipe water leakage detecting method based on wavelet singularity analysis and ARMA model
CN116292367A (en) Power plant fan system abnormal condition detection method based on one-dimensional convolution
KR102110319B1 (en) System for generating learning data
CN110516323A (en) A kind of pressure pipeline damage forecast method based on Time-Series analysis
Zhang et al. Distributed quality-related process monitoring framework using parallel DVIB-VAE-mRMR for large-scale processes
KR101956530B1 (en) Method and apparatus for prediction of wind direction, and method for yaw control of wind turbines using the same
CN111443602B (en) Hidden Markov-based valve health degree estimation and life prediction method
CN112668794A (en) Maintenance optimization method for aircraft engine under uncertain life prediction information
US11339763B2 (en) Method for windmill farm monitoring
CN115828744A (en) White light LED fault on-line diagnosis and service life prediction method
CN105279553A (en) Method for identifying fault degree of high-pressure heater water supply system
CN113919224A (en) Main steam flow prediction method and system
Kumar et al. Gas turbine engine operational data analysis for anomaly detection: Statistical vs. neural network approach

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