CN108896287A - A kind of control valve fault detection method based on multiscale analysis - Google Patents
A kind of control valve fault detection method based on multiscale analysis Download PDFInfo
- Publication number
- CN108896287A CN108896287A CN201810634725.4A CN201810634725A CN108896287A CN 108896287 A CN108896287 A CN 108896287A CN 201810634725 A CN201810634725 A CN 201810634725A CN 108896287 A CN108896287 A CN 108896287A
- Authority
- CN
- China
- Prior art keywords
- section
- control valve
- fitting
- point
- slope
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M13/00—Testing of machine parts
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M3/00—Investigating fluid-tightness of structures
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
Abstract
The present invention relates to a kind of control valve fault detection method more particularly to a kind of control valve fault detection methods based on multiscale analysis, it uses multiple dimensioned straight line fitting, and (mutation) point occurs using transformation fitting scale Step wise approximation failure.First, initial gauges are calculated according to the time series of fault-signal, time series is divided into several segments, the slope over 10 is fitted using least square method to each section, the slope of more adjacent matching line segment simultaneously makees difference, and there are catastrophe points within difference maximum, that is, maximum both ends line segment of slope transformation.Then, in the range of two sections of line segments, the above method is continued to use, gradually reduces the range intervals of fault point, until final fitting Scale Convergence is 1, at this point, the maximum point of slope variation is the fault point of former time series.There is superiority in accuracy of this method when finding control valve fault point and validity and computational efficiency when for sophisticated signal.
Description
Technical field
The invention belongs to actuator fault diagnosis technical fields, are related to a kind of control valve failure inspection based on multiscale analysis
Survey method.
Background technique
In Producing Process of Processing Industry, since often in high temperature and pressure, low temperature is toxic, dielectric ring of crystallization or deep-etching
Border and frequent adjusting movement, cause regulating valve to be easy to appear all kinds of failures.And in the course of work of control fluid flow, it adjusts
Save valve it is sensitive whether, be directly related to the quality of entire control system.Once regulating valve breaks down will likely be to work
Industry production brings huge economic loss, or even causes casualties.Therefore, occur in discovery regulating valve operational process in time
Various failures and to take appropriate measures be to guarantee in industrial process that automatic control system is safe and reliable, basis of stable operation.
Currently, the method for diagnosing faults that regulating valve is commonly used has:1) it repairs after regulating valve breaks down or more
The passivity diagnosis changed;2) when coupling apparatus general overhaul, the expection that all regulating valves are safeguarded and overhauled is regularly scheduled
Diagnosis;However both of these approaches not only result in the interruption of industrial processes, and to still can normal use regulating valve into
It the unified dismounting of row and checks, had not only wasted time and resource, but also delay industrial production, all there is certain for both maintenance modes
Drawback.The angle operated normally from the online i.e. valve of control valve according to original working method, the working signal parameter according to valve
Finding method carries out on-line checking to valve, is a kind of new on-line fault diagnosis thinking.
Summary of the invention
For the demand, it is an object of the invention to propose a kind of control valve fault detection side based on multiscale analysis
Method.
A kind of control valve fault detection method based on multiscale analysis uses more rulers to control valve fault-signal
Straight line fitting is spent, failure is sought and point occurs, it is characterised in that include the following steps:
Step 1:Read the time series x of the position of valve containing control valve fault-signal x1,x2,x3,...,xn, n is fault-signal
The number of sampled point, calculates initial gauges, and initial gauges are denoted as l1,N is natural number set;
Step 2:By entire time series x1,x2,x3,...,xnIt is divided into m sections, each section is denoted as sj, wherein j=1,
2 ..., m, each segment length are l1, m and length l1There are relationships:If the length of final stage is less than l1, that
It incorporates it into previous segment data;
Step 3:To sjThe data sequence of the inside selects least square method to carry out straight line fitting, seeks each sjThe fitting of section
The slope k of straight line, is successively denoted as k1,k2,k3,...,km;
Step 4:The absolute value for successively seeking the slope difference of adjacent segments, is denoted as ej=| kj+1-kj|,1≤j≤m-1;
Step 5:Work as e when seeking failure catastrophe point in section where seeking the catastrophe point of failuremax=ej, original sequence at this time
Arrange s1,s2,s3,...,smIn sjWith sj+1Between variation tendency it is maximum, catastrophe point falls within sj, sj+1I.e. in section (j × l1,j×
l1+l1) within;
Step 6:In the section of step 5), new initial gauges are recalculated, repeat step 2)-step 5), confirmation event
Hinder catastrophe point.
A kind of control valve fault detection method based on multiscale analysis, it is characterised in that the selection in step 3
Least square method carries out straight line fitting, seeks each sjSection straight slope the specific steps are:
Step 3.1:Seek residual sum square;If every section of sjFitting a straight line expression formula be x'i=ki+b, i=1,
2 ..., n, then sequence sjIt is denoted as RSS with the residual sum of squares (RSS) of fitting a straight line, shown in calculation expression such as formula (1):
Step 3.2:The formula of being expressed as (2) will be further spread out on the right of formula (1):
In formula (2):To obtain the expression formula (3) of RSS:
Only k, b are unknown quantity in formula (3), and other parameters are given data;
Step 3.3:Keep residual sum of squares (RSS) RSS minimum, then must and only needAnd Respectively i, xiAverage value;
Step 3.4:Find out every section of sj, j=1, the slope k of 2 ..., m fitting a straight line is successively denoted as k1,k2,k3,...,km。
A kind of control valve fault detection method based on multiscale analysis, it is characterised in that the confirmation in step 6
Failure catastrophe point the specific steps are:
Step 6.1:In sj, sj+1Further trouble-shooting catastrophe point in fault section is research with the sequence of intervals length
Object calculates new initial gauges l again2, l2=floor (l1×0.5);
Step 6.2:It is repeated in the method Step wise approximation catastrophe point of step 2- step 5, until final scale ln=2, this
When the siding-to-siding block length searched be 3, fitting Scale Convergence is 1, and the midpoint in this section is time series x1,x2,x3,...,
xnCatastrophe point.
The present invention is by using above-mentioned technology, with multiple dimensioned straight line fitting, using transformation fitting scale Step wise approximation event
(mutation) point occurs for barrier:Firstly, the time series according to fault-signal calculates initial gauges, time series is divided into several
Section fits the slope over 10 using least square method to each section, and the slope of more adjacent matching line segment simultaneously makees difference, and difference is maximum
I.e. there are catastrophe points within the maximum both ends line segment of slope transformation.Then, it in the range of two sections of line segments, continues to use above-mentioned
Method gradually reduces the range intervals of fault point, until final fitting Scale Convergence is 1, at this point, the maximum point of slope variation is
It is the fault point of former time series, accuracy of this method when finding control valve fault point and when for sophisticated signal
There is superiority in validity and computational efficiency.
Detailed description of the invention
Fig. 1 is the control valve fault detection flow chart of the invention based on multiscale analysis;
Fig. 2 is that film gas chamber gas leakage figure occurs for control valve 330s;
Fig. 3 is the schematic diagram of the x of time series containing fault-signal (n);
Fig. 4 is failure detection result figure.
Specific embodiment
The specific embodiment of the invention is described below now in conjunction with Figure of description and embodiment:
Embodiment:
As shown in Figure 1, a kind of control valve fault detection method based on multiscale analysis of the invention, to control valve failure
Signal uses multiple dimensioned straight line fitting, seeks failure and point occurs, specifically comprise the following steps:
Step 1:Control valve is set, film gas chamber principal fault occurs in 330s, reads control valve valve position fault-signal x's
Time series x1,x2,x3,...,xn, n is the number of the sampled point of fault-signal, calculates initial gauges:
Initial gauges are denoted as l1,N is natural number set;
As shown in Fig. 2, control valve, which is arranged, in step 1 occurs film gas chamber principal fault in 330s;
As shown in figure 3, it is horizontal axis (s) for horizontal axis that step 1, which was obtained using the time, valve position (percentage) is the longitudinal axis, control valve event
Barrier signal sequence sequence is expressed as x (n), n=1,2,600;
Step 2:By entire time series x1,x2,x3,...,xnIt is divided into m sections, each section is denoted as sj, wherein j=1,
2 ..., m, each segment length are l1, m and length l1There are relationships:If the length of final stage is less than l1, that
It incorporates it into previous segment data;
Step 3:To sjThe data sequence of the inside selects least square method to carry out straight line fitting, seeks each sjThe straight line of section
Slope;
Step 3.1:Seek residual sum square;If every section of sjFitting a straight line expression formula be x'i=ki+b, i=1,
2 ..., n, then sequence sjIt is denoted as RSS with the residual sum of squares (RSS) of fitting a straight line, shown in calculation expression such as formula (1):
Step 3.2:The formula of being expressed as (2) will be further spread out on the right of formula (1):
In formula (2):To obtain the expression formula (3) of RSS:
Only k, b are unknown quantity in formula (3), and other parameters are given data;
Step 3.3:Keep residual sum of squares (RSS) RSS minimum, then must and only needAnd Respectively i, xiAverage value;
Step 3.4:Find out every section of sj, j=1, the slope k of 2 ..., m fitting a straight line is successively denoted as k1,k2,k3,...,km;
Step 4:Successively seek the absolute value of the slope difference of adjacent segments:
The absolute value of the slope difference of adjacent segments is denoted as ej=| kj+1-kj|,1≤j≤m-1;
Step 5:Section where seeking the catastrophe point of failure:
It seeks needing e when failure catastrophe pointmax=ej, original series s at this time1,s2,s3,...,smIn sjWith sj+1Between change
Change trend is maximum, and catastrophe point falls within sj, sj+1I.e. in section (j × l1,j×l1+l1) within;
Step 6:Confirm failure catastrophe point:
Step 6.1:In sj, sj+1Further trouble-shooting catastrophe point in fault section is research with the sequence of intervals length
Object calculates new initial gauges l again2, l2=floor (l1×0.5);
Step 6.2:It is repeated in the method Step wise approximation catastrophe point of step 2- step 5, until final scale ln=2, this
When the siding-to-siding block length searched be 3, fitting Scale Convergence is 1, and the midpoint in this section is time series x1,x2,x3,...,
xnCatastrophe point;
As shown in figure 4, step 6 failure detection result is located in the section (330s, 360s), failure detection result is compared with subject to
Really.
A kind of control valve fault detection method based on multiscale analysis proposed by the present invention, it is quasi- with multiple dimensioned straight line
It closes, (mutation) point is occurred using transformation fitting scale Step wise approximation failure.Firstly, the time series according to fault-signal calculates
Time series is divided into several segments by initial gauges, fits the slope over 10 using least square method to each section, more adjacent
The slope of matching line segment simultaneously makees difference, and there are catastrophe points within difference maximum, that is, maximum both ends line segment of slope transformation.Then, at this
In the range of two sections of line segments, the above method is continued to use, gradually reduces the range intervals of fault point, until fitting Scale Convergence is
1, at this point, the maximum point of slope variation is the catastrophe point of former time series.Accuracy of this method in trouble-shooting point, with
And there is superiority in the validity and computational efficiency when for sophisticated signal.
Claims (3)
1. a kind of control valve fault detection method based on multiscale analysis, quasi- using multiple dimensioned straight line to control valve fault-signal
It closes, seeks failure and point occurs, it is characterised in that include the following steps:
Step 1:Read the time series x of the position of valve containing control valve fault-signal x1,x2,x3,...,xn, n is the sampling of fault-signal
The number of point, calculates initial gauges, and initial gauges are denoted as l1,l1∈ N, N are natural number set;
Step 2:By entire time series x1,x2,x3,...,xnIt is divided into m sections, each section is denoted as sj, wherein j=1,2 ..., m,
Each segment length is l1, m and length l1There are relationships:If the length of final stage is less than l1, then being received
Enter into previous segment data;
Step 3:To sjThe data sequence of the inside selects least square method to carry out straight line fitting, seeks each sjThe fitting a straight line of section
Slope k, be successively denoted as k1,k2,k3,...,km;
Step 4:The absolute value for successively seeking the slope difference of adjacent segments, is denoted as ej=| kj+1-kj|,1≤j≤m-1;
Step 5:Work as e when seeking failure catastrophe point in section where seeking the catastrophe point of failuremax=ej, original series s at this time1,
s2,s3,...,smIn sjWith sj+1Between variation tendency it is maximum, catastrophe point falls within sj, sj+1I.e. in section (j × l1,j×l1+l1)
Within;
Step 6:In the section of step 5), new initial gauges are recalculated, repeat step 2)-step 5), confirmation failure is prominent
Height.
2. a kind of control valve fault detection method based on multiscale analysis according to claim 1, it is characterised in that step
Selection least square method in rapid 3 carries out straight line fitting, seeks each sjSection straight slope the specific steps are:
Step 3.1:Seek residual sum square;If every section of sjFitting a straight line expression formula be x 'i=ki+b, i=1,
2 ..., n, then sequence sjIt is denoted as RSS with the residual sum of squares (RSS) of fitting a straight line, shown in calculation expression such as formula (1):
Step 3.2:The formula of being expressed as (2) will be further spread out on the right of formula (1):
In formula (2):To obtain the expression formula (3) of RSS:
Only k, b are unknown quantity in formula (3), and other parameters are given data;
Step 3.3:Keep residual sum of squares (RSS) RSS minimum, then must and only needAnd Respectively i, xiAverage value;
Step 3.4:Find out every section of sj, j=1, the slope k of 2 ..., m fitting a straight line is successively denoted as k1,k2,k3,...,km。
3. a kind of control valve fault detection method based on multiscale analysis according to claim 1, it is characterised in that step
Confirmation failure catastrophe point in rapid 6 the specific steps are:
Step 6.1:In sj, sj+1Further trouble-shooting catastrophe point in fault section, using the sequence of intervals length as research object,
New initial gauges l is calculated again2, l2=floor (l1×0.5);
Step 6.2:It is repeated in the method Step wise approximation catastrophe point of step 2- step 5, until final scale ln=2, it is searched at this time
Target-seeking siding-to-siding block length is 3, and fitting Scale Convergence is 1, and the midpoint in this section is time series x1,x2,x3,...,xnIt is prominent
Height.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810634725.4A CN108896287A (en) | 2018-06-20 | 2018-06-20 | A kind of control valve fault detection method based on multiscale analysis |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810634725.4A CN108896287A (en) | 2018-06-20 | 2018-06-20 | A kind of control valve fault detection method based on multiscale analysis |
Publications (1)
Publication Number | Publication Date |
---|---|
CN108896287A true CN108896287A (en) | 2018-11-27 |
Family
ID=64345185
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810634725.4A Pending CN108896287A (en) | 2018-06-20 | 2018-06-20 | A kind of control valve fault detection method based on multiscale analysis |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108896287A (en) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111884180A (en) * | 2020-07-15 | 2020-11-03 | 华北电力大学 | Direct current system protection method and device based on least square fitting |
CN113285977A (en) * | 2020-08-08 | 2021-08-20 | 詹能勇 | Network maintenance method and system based on block chain and big data |
CN113392378A (en) * | 2021-07-16 | 2021-09-14 | 中南大学 | Surrounding rock deformation multipoint mutation identification method and system based on time sequence |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1843822A (en) * | 2006-04-30 | 2006-10-11 | 西安英卓电子科技有限公司 | Dynamic detection system of train wheel pair and detection method thereof |
CN101493378A (en) * | 2009-03-06 | 2009-07-29 | 湖南大学 | Gear failure diagnosis method based on multi-dimension line frequency modulation base sparse signal decomposition |
WO2013180723A1 (en) * | 2012-05-31 | 2013-12-05 | University Of Connecticut | Methods apparatuses for defect diagnosis in a mechanical system |
CN105113375A (en) * | 2015-05-15 | 2015-12-02 | 南京航空航天大学 | Pavement cracking detection system and method based on line structured light |
CN107064315A (en) * | 2017-03-24 | 2017-08-18 | 北京工业大学 | Blade cracks position and depth recognition method based on wavelet analysis |
-
2018
- 2018-06-20 CN CN201810634725.4A patent/CN108896287A/en active Pending
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1843822A (en) * | 2006-04-30 | 2006-10-11 | 西安英卓电子科技有限公司 | Dynamic detection system of train wheel pair and detection method thereof |
CN101493378A (en) * | 2009-03-06 | 2009-07-29 | 湖南大学 | Gear failure diagnosis method based on multi-dimension line frequency modulation base sparse signal decomposition |
WO2013180723A1 (en) * | 2012-05-31 | 2013-12-05 | University Of Connecticut | Methods apparatuses for defect diagnosis in a mechanical system |
CN105113375A (en) * | 2015-05-15 | 2015-12-02 | 南京航空航天大学 | Pavement cracking detection system and method based on line structured light |
CN107064315A (en) * | 2017-03-24 | 2017-08-18 | 北京工业大学 | Blade cracks position and depth recognition method based on wavelet analysis |
Non-Patent Citations (1)
Title |
---|
黄静: "多尺度直线拟合法在时间序列突变点检测中的应用", 《兵工学报》 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111884180A (en) * | 2020-07-15 | 2020-11-03 | 华北电力大学 | Direct current system protection method and device based on least square fitting |
CN111884180B (en) * | 2020-07-15 | 2021-05-18 | 华北电力大学 | Direct current system protection method and device based on least square fitting |
CN113285977A (en) * | 2020-08-08 | 2021-08-20 | 詹能勇 | Network maintenance method and system based on block chain and big data |
CN113285977B (en) * | 2020-08-08 | 2022-07-05 | 山东鼹鼠人才知果数据科技有限公司 | Network maintenance method and system based on block chain and big data |
CN113392378A (en) * | 2021-07-16 | 2021-09-14 | 中南大学 | Surrounding rock deformation multipoint mutation identification method and system based on time sequence |
CN113392378B (en) * | 2021-07-16 | 2024-04-09 | 中南大学 | Surrounding rock deformation multipoint mutation identification method and system based on time sequence |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US9310288B2 (en) | Systems and methods to monitor operating processes | |
CA2679128C (en) | Method and system of using inferential measurements for abnormal event detection in continuous industrial processes | |
CN108896287A (en) | A kind of control valve fault detection method based on multiscale analysis | |
CN109508818B (en) | Online NOx prediction method based on LSSVM | |
US20120330631A1 (en) | Method and System of Using Inferential Measurements for Abnormal Event Detection in Continuous Industrial Processes | |
CN113330463A (en) | Chamber matching in semiconductor device tools using neural networks | |
CN104390657A (en) | Generator set operating parameter measuring sensor fault diagnosis method and system | |
WO2016208315A1 (en) | Plant diagnosis device and plant diagnosis method | |
WO2023207190A1 (en) | Fault early warning system and method for fossil fuel power plant, electronic device, and storage medium | |
CN110779745A (en) | Heat exchanger early fault diagnosis method based on BP neural network | |
CN106200624B (en) | Based on the Industrial Boiler method for diagnosing faults for intersecting segmentation PCA | |
CN117029968A (en) | Traffic data diagnosis method, system, storage medium and electronic equipment | |
KR102185017B1 (en) | Integrated information management and predictive maintenance system for gas supply facilities | |
CN104062903A (en) | Carbon fiber coagulating bath reconstitution controller based on immunologic mechanism | |
CN112541296A (en) | SO2 prediction method based on PSO-LSSVM | |
CN113836794B (en) | Soft and hard combined fly ash carbon content online monitoring method | |
KR100340967B1 (en) | Method and apparatus for detecting fault using General Regression Neural Network | |
CN112328590B (en) | Deep cleaning method for operation data of thermal equipment | |
JPH06303677A (en) | Pipe network management method and system utilizing the same | |
RU2794126C1 (en) | Device for determining blast furnace failure, method for detecting blast furnace fault and method for operating blast furnace | |
CN114117934B (en) | Power station boiler water-cooled wall temperature on-line prediction method based on gated neural network | |
Martini et al. | Gross error detection based on Serial Elimination: Applications to an industrial Gas Turbine | |
RU2801843C1 (en) | Information-measuring control system for automatic control of the temperature parameters of gas field facilities | |
KR20230095408A (en) | Display manufacturing equipment malfunction real-time inspection system | |
KR102561073B1 (en) | Smart Performance Monitoring System Using Real-time Efficiency Control Limit for Industrial Boilers |
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 | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20181127 |