CN116611609A - Equipment stability state prediction and evaluation method based on monitoring parameters - Google Patents
Equipment stability state prediction and evaluation method based on monitoring parameters Download PDFInfo
- Publication number
- CN116611609A CN116611609A CN202310415583.3A CN202310415583A CN116611609A CN 116611609 A CN116611609 A CN 116611609A CN 202310415583 A CN202310415583 A CN 202310415583A CN 116611609 A CN116611609 A CN 116611609A
- Authority
- CN
- China
- Prior art keywords
- value
- sequence
- evaluation
- evaluation index
- equipment
- 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
- 238000011156 evaluation Methods 0.000 title claims abstract description 218
- 238000012544 monitoring process Methods 0.000 title claims abstract description 26
- 238000000034 method Methods 0.000 claims abstract description 98
- 238000010248 power generation Methods 0.000 claims abstract description 60
- 230000009467 reduction Effects 0.000 claims abstract description 30
- 238000000354 decomposition reaction Methods 0.000 claims abstract description 19
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 11
- 238000010219 correlation analysis Methods 0.000 claims abstract description 7
- 239000011159 matrix material Substances 0.000 claims description 52
- 238000012216 screening Methods 0.000 claims description 45
- 230000014509 gene expression Effects 0.000 claims description 41
- 230000008569 process Effects 0.000 claims description 20
- 238000004364 calculation method Methods 0.000 claims description 12
- YHXISWVBGDMDLQ-UHFFFAOYSA-N moclobemide Chemical compound C1=CC(Cl)=CC=C1C(=O)NCCN1CCOCC1 YHXISWVBGDMDLQ-UHFFFAOYSA-N 0.000 claims description 12
- 238000012545 processing Methods 0.000 claims description 11
- 230000006641 stabilisation Effects 0.000 claims description 9
- 238000011105 stabilization Methods 0.000 claims description 9
- 238000000605 extraction Methods 0.000 claims description 8
- 239000002131 composite material Substances 0.000 claims description 6
- 238000003064 k means clustering Methods 0.000 claims description 5
- 238000009825 accumulation Methods 0.000 claims description 3
- 239000010687 lubricating oil Substances 0.000 description 16
- 238000012360 testing method Methods 0.000 description 10
- 230000008859 change Effects 0.000 description 9
- 238000010586 diagram Methods 0.000 description 7
- 239000013598 vector Substances 0.000 description 7
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 6
- 230000033228 biological regulation Effects 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 2
- 230000008878 coupling Effects 0.000 description 2
- 238000010168 coupling process Methods 0.000 description 2
- 238000005859 coupling reaction Methods 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 239000003921 oil Substances 0.000 description 2
- 229920006395 saturated elastomer Polymers 0.000 description 2
- 238000013112 stability test Methods 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 1
- 238000007476 Maximum Likelihood Methods 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000001174 ascending effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000003638 chemical reducing agent Substances 0.000 description 1
- 238000007621 cluster analysis Methods 0.000 description 1
- 230000003750 conditioning effect Effects 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 238000007499 fusion processing Methods 0.000 description 1
- 210000004907 gland Anatomy 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 238000002156 mixing Methods 0.000 description 1
- 239000002243 precursor Substances 0.000 description 1
- 238000004445 quantitative analysis Methods 0.000 description 1
- 238000012552 review Methods 0.000 description 1
- 238000007789 sealing Methods 0.000 description 1
- 230000036962 time dependent Effects 0.000 description 1
- 238000012549 training Methods 0.000 description 1
- 238000011426 transformation method Methods 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
Classifications
-
- 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/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
-
- 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
- G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/10—Pre-processing; Data cleansing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/21—Design or setup of recognition systems or techniques; Extraction of features in feature space; Blind source separation
- G06F18/211—Selection of the most significant subset of features
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/21—Design or setup of recognition systems or techniques; Extraction of features in feature space; Blind source separation
- G06F18/213—Feature extraction, e.g. by transforming the feature space; Summarisation; Mappings, e.g. subspace methods
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/23—Clustering techniques
- G06F18/232—Non-hierarchical techniques
- G06F18/2321—Non-hierarchical techniques using statistics or function optimisation, e.g. modelling of probability density functions
- G06F18/23213—Non-hierarchical techniques using statistics or function optimisation, e.g. modelling of probability density functions with fixed number of clusters, e.g. K-means clustering
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/27—Regression, e.g. linear or logistic regression
-
- 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
-
- 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
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2123/00—Data types
- G06F2123/02—Data types in the time domain, e.g. time-series data
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
- Y04S—SYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
- Y04S10/00—Systems supporting electrical power generation, transmission or distribution
- Y04S10/50—Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Business, Economics & Management (AREA)
- General Physics & Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Human Resources & Organizations (AREA)
- Economics (AREA)
- Strategic Management (AREA)
- Computer Vision & Pattern Recognition (AREA)
- General Engineering & Computer Science (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Evolutionary Biology (AREA)
- Artificial Intelligence (AREA)
- Evolutionary Computation (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Business, Economics & Management (AREA)
- Tourism & Hospitality (AREA)
- Entrepreneurship & Innovation (AREA)
- Marketing (AREA)
- Mathematical Physics (AREA)
- Quality & Reliability (AREA)
- Operations Research (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Health & Medical Sciences (AREA)
- Pure & Applied Mathematics (AREA)
- Game Theory and Decision Science (AREA)
- Computational Mathematics (AREA)
- Development Economics (AREA)
- General Health & Medical Sciences (AREA)
- Software Systems (AREA)
- Databases & Information Systems (AREA)
- Algebra (AREA)
- Computing Systems (AREA)
- Primary Health Care (AREA)
- Water Supply & Treatment (AREA)
- Public Health (AREA)
- Educational Administration (AREA)
Abstract
The invention discloses a method for predicting and evaluating equipment stability state based on monitoring parameters, which is specially used for a power generation turbine of a steam power ship. Decomposing the initial time sequence by adopting an MREMD decomposition method, and providing a non-stationary time sequence single-parameter prediction model based on MREMD decomposition and wavelet threshold noise reduction, so that the prediction precision of equipment is improved; the parameter time sequence of each evaluation index is converted into a fluctuation value sequence, so that positive and negative ideal solutions of each evaluation index of the equipment can be determined, the correlation weight of the weighted evaluation parameter sequence of each evaluation index is obtained by combining an improved TOPSIS algorithm and the existing B-type gray correlation analysis method, the final running stability state score of each prediction trend item is obtained, each evaluation index is adjusted according to the stability state score, and the overall stability of the equipment in a future period can be ensured by timely adjusting the correlation parameters.
Description
Technical Field
The invention belongs to the technical field of system stability evaluation, and particularly relates to a device stability state prediction and evaluation method based on monitoring parameters.
Background
The coupling between the parameters related to the running state of the current system is strong, so that equipment is difficult to eliminate, and the cost is high. At present, no mature method for predicting future parameter changes based on the current running state of the system exists, so that state evaluation and maintenance of the system are mainly based on state parameters recorded in the moment of system instability or instability, and in fact, before the instability phenomenon occurs, relevant parameters of the system often have precursors such as state parameter drift, running characteristic change and the like, if the state of equipment cannot be found and adjusted in time, the state of the equipment may be further deteriorated, and finally various faults occur.
For example, the unstable operation of power turbines is particularly evident in extreme variable conditions, cross steam supply and special operating environments of steam powered and nuclear powered vessels. Although the monitoring system ensures the stability of the operation state of the power generation turbine by strictly monitoring the operation state parameters of all the equipment in the system in real time, the problems that the parameter state transition time is longer, the preset control parameters and logic cannot adapt to the technical state change of the equipment and the like still exist due to factors such as thermal inertia and the like. For a power generation turbine in a steam power generation system, the parameters influencing the running condition of the power generation turbine are numerous, the coupling relation is complex, the current monitoring system is difficult to ensure that the power generation turbine can maintain the stable state under all special conditions and disturbance, the power generation turbine lacks corresponding early warning and pre-adjustment capability when running is about to be unstable, and the power generation turbine can only rely on a large number of monitoring personnel to make local adjustment for instability alarms generated when the parameter running in the system is unstable, so that the labor and equipment running cost is huge.
In the aspect of multi-parameter correlation analysis, the entropy weight-ideal solution describes the proximity degree of two evaluation objects through distance measurement, and the method can only show the position relation and can not show the situation change situation among the data sequences of the objects to be evaluated; the gray correlation algorithm describes the correlation by the geometric similarity degree of curves of two study objects, and the method can only reflect the situation change condition of a data curve and cannot reflect the position relation; by adopting the parameter correlation calculation model based on the entropy weight method, the ideal solution and the gray correlation analysis, the similarity of the change trend and the position relation between the parameters can be considered on the basis of determining the importance degree of the parameters.
The prior TOPSIS method constructs a positive ideal solution and a negative ideal solution according to main characteristics of a study object, calculates distances between each evaluation index of the equipment and the positive ideal solution and the negative ideal solution, and then takes the distance degree close to the positive ideal solution and the negative ideal solution as a basis for evaluating the importance degree of the influence factors. The approach degree of two evaluation objects is described through distance measurement, the method can only show the position relation, the situation change situation among the data sequences of the objects to be evaluated can not be shown, and for the data containing non-stationary parameters, the negative ideal solution can not be obtained by using the existing TOPSIS method, and the method can not be applied to the stability evaluation of equipment containing the non-stationary parameters.
It is therefore desirable to provide an improved TOPSIS algorithm by which to combine the two methods to stabilize the operation of a device for a future period of time with device parameters that contain non-stationary parameters, in combination with existing type B gray correlation analysis.
Disclosure of Invention
The invention provides a device stability state prediction and evaluation method based on monitoring parameters, which is used for realizing the prediction and evaluation of the device running stability state containing non-stable parameters, and timely adjusting the running operation of related devices through the stability state score, so that the running stability of the devices in a future period of time is improved.
The invention discloses a device stability state prediction and evaluation method based on monitoring parameters, which comprises the following steps:
s1, analyzing an instability factor affecting the running state of equipment, selecting a plurality of evaluation indexes related to the running stability of the equipment, collecting data sets with the length of N of each evaluation index in a period of time, and forming each data set into an initial time sequence of the running of the evaluation index;
s2, decomposing the initial time sequence by adopting an MREMD decomposition method to obtain n IMF components and a final residual component ;
S3, setting a screening threshold value to screen the IMF components, carrying out wavelet threshold value noise reduction treatment on the IMF components meeting the screening threshold value, recombining to obtain a new IMF component matrix, and processing the new IMF component matrix to obtain a final prediction trend item;
s4, converting the new IMF component matrix to obtain a parameter time sequence of each evaluation index, obtaining a weighted evaluation parameter sequence of each evaluation index by using an entropy weight method, and obtaining positive and negative ideal solutions of each moment point by the weighted evaluation parameter sequence;
s5, respectively calculating the distance relation value between the fluctuation value of each evaluation index and the positive and negative ideal solutions;
s6, respectively calculating gray correlation degrees of fluctuation values and positive and negative ideal solutions of all evaluation indexes by adopting a B-type gray correlation analysis method;
s7, carrying out dimensionless treatment on the distance relation value and the gray correlation value respectively, and fusing the gray correlation value and the distance relation value after the dimensionless treatment into a correlation weight of each evaluation index and the integral running stability state of the equipment;
s8, establishing an operation stability state scoring formula of each evaluation index on the basis of the predicted trend item, obtaining the scoring value of each evaluation index at each moment point in a predicted time period, and establishing an operation scoring formula of the predicted trend item of the whole equipment according to the scoring value and the correlation weight of each evaluation index to obtain the operation scoring of the predicted trend item of the whole equipment;
S9, setting a total reference index, and comparing the operation evaluation score of the overall predicted trend item of the equipment with the total reference index to judge whether adjustment of each evaluation index is needed, so that the operation evaluation score of the overall predicted trend item of the equipment is higher than the total reference index finally.
Further, the step S2 includes the steps of:
s21, utilizing an AR prediction model to perform initial time sequence on each evaluation indexPredicting and extending the left and right endpoints to obtain a new sequence so as to enable +.>Is between two extreme points of the new sequence after extension, wherein the AR prediction model expression is as follows:
(1)
in the above-mentioned method, the step of,the predicted point is the t moment; />The actual values of the parameters at the times t-1, t-2, & gtt-p, & lt/L, & gt>For p+1 real numbers, ">(t=p+1, p+2, & gtis, N) white noise sequence with zero mean;
s22, calculating the average value between adjacent extreme points in the new sequence, wherein the average value point obtained in the new sequence is the initial time sequenceMean point sequence>Wherein k is->The number of the mean value points is k+1, and the +.f is obtained through 3 times spline interpolation>Signal mean sequence>Make->Minus- >Obtain 1 order signal component->The expression is shown as follows:
(2)
s23, calculating the obtained signal componentAs an original signal, the above steps S21 and S22 are repeated for iterative calculation until the signal component after r iterations +.>And the termination condition is satisfied, and the judgment formula of the calculation termination condition is as follows:
(3)
in the above-mentioned method, the step of,for normalizing standard deviation, ++>And->For the signal component h after the r-1 th iteration respectively 1,r-1 The mean point sequence of (t) and the signal component after the r-th iteration +.>The standard deviation of the mean point sequence of (2), P is the conditional probability; definition according to sigma principle->For the initial time sequence->A sense value of the absolute value of the middle extreme point; />And->The mean point of the time series after the initial time series and the r iteration and +.>Ratio sequence of (2); said->、/>And->The formulas of (a) are respectively shown in the following formulas;
(4)
wherein:the z extreme point of the new sequence after extension; />For the initial time sequence->Is a k mean points of (2); />For time series component after the r-1 th iteration +.>Is a k mean points of (2);
s24, willAs an eigenmode function IMF of order 1 1 Output of->Minus->Get component->And willThe above steps S21 to S24 are repeated as an original signal until the residual signal becomes a monotone function or no new IMF component is separated, and the residual signal extraction process is as follows:
(5)
Wherein n is the number of resolvable IMF components;as the final residual component;
s25, after the decomposition, the initial time sequenceCan be expressed as the n IMF components and the final residual component->The sum is as follows:
(6) 。
further, the step S3 includes the steps of:
s31, according to the IMF components and the initial time sequenceCorrelation coefficient of->And root mean square errorObtaining the screening threshold value;
s32, screening the n IMF components by using the screening threshold, and dividing the IMF components into IMF components which do not meet the screening threshold and IMF components which meet the screening threshold;
s33, carrying out wavelet threshold denoising treatment on the components meeting the screening threshold, and recombining the residual components, the IMF components which do not meet the screening threshold and the IMF components subjected to denoising treatment into a new IMF component matrix;
s34, performing feature extraction and dimension reduction on the new IMF component matrix by using a singular value decomposition method to obtain a singular value component matrix;
s35, calculating and obtaining the arrangement entropy value of each singular value component matrix by using an arrangement entropy algorithm;
s36, dividing the new IMF component matrix into two types according to the arrangement entropy value by utilizing a K-means clustering algorithm, taking a type of component with a lower entropy value, reconstructing the type of component into an original trend item, carrying out differential stabilization treatment on the original trend item, obtaining an ARIMA model by fitting, and predicting the original trend item by the ARIMA model to obtain a final predicted trend item.
Further, the step S31 includes the steps of:
s311, calculating the IMF components of each order and an initial time sequenceCorrelation coefficient of->And root mean square errorThe expression is shown as follows:
(7)
in the middle ofAnd->Initial time series>And the average value of IMF components of each order, N is the sequence length of each parameter, and N is the number of IMF components;
s312, calculating a screening threshold value, wherein the screening threshold value comprises a correlation coefficient threshold valueAnd mean square error threshold->The expressions are shown below:
(8)
s313, selecting IMF components meeting the following conditions, and carrying out wavelet threshold denoising to obtain denoised IMF components, wherein the judgment formula of the screening conditions is as follows:
(9)。
further, the step S32 includes the steps of:
s321, estimating the noise of each layer in the IMF components meeting the screening conditions through wavelet detail coefficients of each layer, and simultaneously utilizing the coefficientsThe threshold value of the detail coefficient is lowered layer by layer, and the selection criterion formula of the threshold value is as follows:
(10)
in the method, in the process of the invention,a noise threshold that is a layer j wavelet detail coefficient; />The detail coefficient of the wavelet of the j-th layer; media is an intermediate function, i.e. each coefficient is arranged in descending order, and takes the value of the intermediate number when the coefficient number is odd or takes the value of the intermediate number when the coefficient number is even The mean value of the two numbers; n is the number of data points collected by each evaluation index;
s322, performing threshold processing by using a composite threshold function to obtain an IMF component subjected to wavelet noise reduction, wherein the expression of the composite threshold function is as follows:
(11)
in the method, in the process of the invention,the j-th layer detail coefficient after the noise reduction treatment; sign is a sign function, a e [0,1] is an adjustment coefficient, a=0 is a hard threshold function, and a=1 is a soft threshold function.
Further, the step S4 includes the steps of:
s41, sequentially carrying out accumulation reconstruction on an IMF component subjected to noise reduction of each parameter, a component not subjected to noise reduction treatment and a residual component to form a noise-reduced parameter time sequence, converting the noise-reduced parameter time sequence into a fluctuation value sequence according to a setting value of each evaluation index, obtaining an entropy weight value of the fluctuation value sequence by using an entropy weight method, and obtaining a weighted evaluation parameter sequence of each evaluation index by using the entropy weight value and the fluctuation value sequence;
s42, respectively selecting the point with the largest positive fluctuation value of the weighted evaluation parameter sequence of each evaluation index at the same moment as a positive ideal solution, and taking the point with the largest negative fluctuation value as a negative ideal solution, wherein the calculated formulas of the positive ideal solution and the negative ideal solution are as follows:
(12)
In the method, in the process of the invention,for positive understanding, add>Is a negative ideal solution;/>,/>n represents the number of evaluation indexes.
Further, the expression of the distance relation between the fluctuation value of each evaluation index sequence and the positive and negative ideal solutions is as follows:
(13)
wherein i=1, 2, …, N; j=1, 2, …, n, n represents the number of evaluation indexes,for the distance relation value between the fluctuation value of the j-th evaluation index and the positive ideal solution,/for the j-th evaluation index>And the relationship value of the distance between the fluctuation value of the j-th evaluation index and the negative ideal solution.
Further, the expression of the operation stability state scoring formula of each evaluation index is as follows:
(14)
wherein j=1, 2, …, n;for the weighted evaluation parameter of the jth evaluation index in the ith point in time, +.>Andrespectively representing the positive maximum fluctuation value and the negative maximum fluctuation value of each evaluation index in the running range of the equipment.
Further, the step S9 includes the steps of:
calculating the running evaluation score of the overall predicted trend item of the equipment, and setting a total reference index;
comparing the operation evaluation score of the overall predicted trend item of the equipment with the total reference index;
when the operation evaluation score of the overall prediction trend item of the equipment is larger than or equal to the total reference index, the equipment is indicated to operate well in a future prediction time period;
When the operation evaluation score of the overall prediction trend item of the equipment is smaller than the total reference index, the operation instability of the equipment in a future prediction time period is indicated, and the evaluation index with the lower score is timely adjusted, so that the operation evaluation score of the overall prediction trend item of the equipment is greater than or equal to the total reference index.
Further, the equipment stability state prediction and evaluation method based on the monitoring parameters is specially used for the power generation steam turbine of the steam power ship.
Compared with the prior art, the invention has the beneficial effects that:
1. the invention provides a non-stationary time sequence single parameter prediction model based on MREMD decomposition and wavelet threshold noise reduction, which improves the prediction precision of the model; in the definition of noise IMF components of each IMF component after MREMD decomposition, a subjective IMF component definition method for removing only a first-order IMF component to achieve the purpose of noise reduction is abandoned, and two indexes of a correlation coefficient and a root mean square error are introduced to realize quantitative definition of the IMF component, so that the calculated error is reduced;
2. the invention combines the improved TOPSIS method with the existing B-type gray correlation degree, and converts the noise-reduced parameter time sequence into a fluctuation value sequence according to the setting value of each evaluation index, so that the positive and negative ideal solutions of each evaluation index of the equipment can be determined, and the problem that the negative ideal solutions cannot be obtained when the existing TOPSIS method is applied to the stability evaluation of the equipment containing non-stable parameters is solved;
3. In the method, an improved distance relation value formula is provided when the positive and negative ideal solutions in the improved TOPSIS method are calculated, and the correlation between the fluctuation value of each evaluation index and the positive and negative ideal solutions is calculated through the formula; in addition, as the distance between the two is smaller, the evaluation index is unstable, the weight occupied by the evaluation index in the stability of the whole equipment is larger, namely the weight of the evaluation index and the corresponding distance relation of the evaluation index are in inverse proportion relation, the negative sign is added in the formula of the distance relation value, so that the real situation can be better reflected, and the obtained result is more accurate;
4. according to the invention, through establishing the equipment operation evaluation score and evaluating the whole predicted trend item of the equipment, the state stability evaluation score curve is generated through software, so that the operation condition of the power generation turbine in a period of time in the future can be intuitively observed, the operation state of the equipment can be observed in time by operation operators, and the labor and the equipment operation cost are saved.
Drawings
FIG. 1 is a flow chart of a single evaluation index timing prediction for a power generation turbine of a steam powered vessel of the present invention;
FIG. 2 shows 7 IMF components and 1 residual component R obtained by decomposing the rotational speed of the power turbine of the present embodiment n ;
FIG. 3 is an initial time series x of 7 IMF components of the turbine speed for power generation according to the present embodiment 0 A correlation coefficient threshold screening map of (t);
FIG. 4 is an initial time series x of 7 IMF components of the turbine speed of the power generation turbine of the present embodiment 0 A mean square error threshold screening graph of (t);
FIG. 5 is a schematic diagram of singular value components of a new IMF component matrix of the power turbine speed of the present embodiment;
FIG. 6 is a graph showing the variation of mutual information values of components of the rotational speed of the power generation turbine with delay time according to the present embodiment;
FIG. 7 is a graph showing the variation of pseudo-neighbor ratios of the various components of the speed of the power turbine according to the present embodiment with the embedding dimension;
FIG. 8 is a schematic diagram of the rotational speed of the power turbine and the original trend thereof according to the present embodiment;
FIG. 9 is a schematic diagram of the rotational speed of the power turbine, the raw trend term and the predicted trend term according to the present embodiment;
FIG. 10 is a diagram showing a positive ideal understanding of the weighted evaluation parameter sequence of the evaluation index according to the present embodiment;
FIG. 11 is a diagram showing a negative idealized solution of a weighted evaluation parameter sequence of the evaluation index according to the present embodiment;
Fig. 12 is a schematic diagram showing evaluation scores of the stability state of the power turbine of the present embodiment over a future period of time.
Detailed Description
The invention is described in further detail below with reference to figures 1 to 12 and the specific examples.
The equipment stability state prediction and evaluation method based on the monitoring parameters is specially used for the power generation steam turbine of the steam power ship.
As shown in fig. 1, a device stability state prediction and evaluation method based on monitoring parameters in this embodiment includes the following steps:
s1, analyzing instability factors affecting the running state of equipment, selecting a plurality of evaluation indexes related to the running stability of the equipment, collecting data sets with the length N of each evaluation index in a monitoring system in a period of time, and forming an initial time sequence x of the running of the evaluation indexes by each data set 0 (t)。
The power generation turbine unit of the embodiment comprises a turbine body, a steam stop valve, a speed reducer, a speed regulation system, a lubricating oil system, a steam seal air extraction system and a condensate system.
The fluctuation phenomenon of the monitoring parameters such as condenser vacuum and water level, saturated steam pressure, temperature, lubricating oil pressure, lubricating oil temperature and the like is frequently encountered in the use process of the marine power generation turbine, and equipment damage or valve speed closing can be caused when the fluctuation phenomenon is serious.
Typical instabilities in the operation of a power turbine include unstable condenser vacuum, unstable condenser water level, unstable lubricating oil pressure, unstable lubricating oil temperature, unstable steam parameters of the turbine entering and exiting steam, unstable seal pressure and unstable turbine rotational speed.
The power generation steam turbine is provided with a monitoring system, and in the monitoring system, monitoring parameters directly related to the running state of the power generation steam turbine mainly comprise condenser vacuum, front-nozzle steam pressure, front-nozzle steam temperature, rear-nozzle steam pressure, rear-nozzle steam temperature, rear-turbine bypass valve steam pressure, rear-turbine bypass valve steam temperature, turbine gear unit lubricating oil pressure, turbine gear unit lubricating oil temperature, front-rear pressure difference of a turbine lubricating oil filter, power generation steam turbine rotating speed, condenser water level, condenser vacuum degree, 1# to 3# nozzle rear-steam pressure (the power generation steam turbine unit in the embodiment comprises three nozzles), speed valve front-closing steam temperature, nozzle front-steam pressure, speed valve rear-closing steam temperature, turbine unit adjusting stage steam pressure, turbine unit adjusting stage steam temperature, gland sealing pressure and the like. The parameters are the evaluation indexes.
The above parameters are analyzed and a plurality of evaluation indexes related to the running stability of the equipment are selected from the above parameters:
(1) Because the certain type of power generation steam turbine unit adopts saturated steam, the invention has corresponding relation with the steam pressure of corresponding measuring points, and in quantitative analysis, the influence on the running stability state of the power generation steam turbine can be completely converted into the steam pressure, thus the invention does not include the selection range of parameters in analysis.
(2) The parameters such as the front steam pressure of the nozzle of the turbine unit, the rear steam pressure of the bypass valve of the turbine unit and the like have a fixed correlation, if the parameters are selected too much, the redundancy of the parameters is easy to be caused, the calculated amount and the calculated time are increased, and the front steam pressure of the nozzle is an important regulation parameter in the operation process of the power generation turbine, so that only the front steam pressure of the nozzle is reserved in the parameters.
(3) The lubricating oil temperature of the turbine gear unit and the front-rear pressure difference of the turbine lubricating oil filter are fixed relative parameters, and in the operation process of the steam power generation system, the fluctuation range of the lubricating oil temperature of the turbine gear unit is larger than that of the front-rear pressure difference of the turbine lubricating oil filter, so that the lubricating oil temperature of the turbine gear unit is only reserved, and the lubricating oil temperature is called as the lubricating oil temperature in the subsequent calculation.
(4) In the system operation process, the consistency of the steam pressure fluctuation range behind the nozzles 1# to 3# is weaker than that of the operation change trend, so that the three evaluation indexes in the evaluation of the operation state evaluation indexes of the power generation steam turbine are all required to be reserved.
The final evaluation indexes which can be used as the evaluation indexes of the running stability state of the power generation turbine in the embodiment are 10 evaluation indexes, namely the turbine rotating speed, the condenser vacuum, the steam pressure after the No. 1 nozzle to the No. 3 nozzle, the lubricating oil temperature, the steam pressure after the adjusting stage, the steam pressure before the quick valve closing, the condenser water level and the steam seal pressure.
In this embodiment, the data sets with the length N of the 10 evaluation indexes collected by the monitoring system within a period of time are recorded respectively, and each data set is formed into an initial time sequence of the evaluation index operationIn this embodiment, the number of data points collected by each evaluation index is the same, and the number of data points collected by each evaluation index is 940, and s is used as a unit.
S2, decomposing the initial time sequence by adopting an MREMD decomposition method to obtain 7 IMF components and a final residual component。
In this embodiment, step S2 includes the following steps:
s21, utilizing an AR prediction model to perform initial time sequence on each evaluation indexPredicting and extending the left and right endpoints to obtain a new sequence so as to enable +.>Is between two extreme points of the new sequence after extension, wherein the AR prediction model expression is as follows:
(1)
in the above-mentioned method, the step of,the predicted point is the t moment; />The actual values of the parameters at the times t-1, t-2, & gtt-p, & lt/L, & gt>For p+1 real numbers, ">(t=p+1, p+2, & gtis, N) white noise sequence with zero mean;
s22, calculating the average value between adjacent extreme points in the new sequence, wherein the average value point obtained in the new sequence is the initial time sequenceMean point sequence>Wherein k is->The number of the mean value points is k+1, and the +.f is obtained through 3 times spline interpolation>Signal mean sequence>Make->Minus->Obtain 1 order signal component->The expression is shown as follows:
(2)
s23, calculating the obtained signal componentAs an original signal, the above steps S21 and S22 are repeated for iterative calculation until the signal component after r iterations +. >And the termination condition is satisfied, and the judgment formula of the calculation termination condition is as follows:
(3)
in the above-mentioned method, the step of,for normalizing standard deviation, ++>And->For the signal component h after the r-1 th iteration respectively 1,r-1 The mean point sequence of (t) and the signal component after the r-th iteration +.>The standard deviation of the mean point sequence of (2), P is the conditional probability; definition according to sigma principle->For the initial time sequence->A sense value of the absolute value of the middle extreme point; />And->The mean point of the time series after the initial time series and the r iteration and +.>Ratio sequence of (2); said->、/>And->The formulas of (a) are respectively shown in the following formulas;
(4)
wherein:the z extreme point of the new sequence after extension; />For the initial time sequence->Is a k mean points of (2); />For time series component after the r-1 th iteration +.>Is a k mean points of (2);
s24, willAs an eigenmode function IMF of order 1 1 Output of->Minus->Get component->And willThe above steps S21 to S24 are repeated as an original signal until the residual signal becomes a monotone function or no new IMF component is separated, and the residual signal extraction process is as follows:
(5)
wherein n is the number of resolvable IMF components, and n is 7 in the embodiment as shown in FIG. 2;as the final residual component;
S25, after the decomposition, the initial time sequenceCan be expressed as the 7 IMF components and the final residual component->The sum is as follows:
(6) 。
in this embodiment, adoptDecomposing the time sequence of the rotating speed of the power generation turbine within 0-940S by using the MREMD decomposition method of formulas (1) - (6) to obtain 7 IMF components and one residual componentThe results are shown in FIG. 2.
S3, setting a screening threshold value to screen the IMF components, carrying out wavelet threshold value noise reduction processing on the IMF components meeting the screening threshold value to obtain wavelet noise-reduced IMF components, and mixing the wavelet noise-reduced IMF components with IMF components which do not meet the screening threshold value and residual componentsAnd (3) recombining to obtain a new IMF component matrix, and processing the new IMF component matrix to obtain a final predicted trend item.
In this embodiment, the step S3 includes the following steps:
s31, according to the IMF components and the initial time sequenceCorrelation coefficient of->And root mean square errorObtaining the screening threshold value;
s32, screening 7 IMF components by using the screening threshold, and dividing the IMF components into IMF components which do not meet the screening threshold and IMF components which meet the screening threshold;
s33, carrying out wavelet threshold denoising treatment on the components meeting the screening threshold, and recombining the residual components, the IMF components which do not meet the screening threshold and the IMF components subjected to denoising treatment into a new IMF component matrix;
S34, performing feature extraction and dimension reduction on the new IMF component matrix by using a singular value decomposition method to obtain a singular value component matrix;
s35, calculating and obtaining the arrangement entropy value of each singular value component matrix by using an arrangement entropy algorithm;
s36, dividing the new IMF component matrix into two types according to the arrangement entropy value by utilizing a K-means clustering algorithm, taking a type of component with a lower entropy value, reconstructing the type of component into an original trend item, carrying out differential stabilization treatment on the original trend item, obtaining an ARIMA model by fitting, and predicting the original trend item by the ARIMA model to obtain a final predicted trend item.
In the embodiment, the above formulas (7) - (8) are adopted to calculate IMF components and initial signals of each order after the rotation speed of the power generation turbine is decomposedCorrelation coefficient of->Root mean square error->Correlation coefficient screening threshold ∈ ->And root mean square error screening threshold->Wherein the calculated IMF components of each order are combined with the initial signal +.>Correlation coefficient of->Correlation coefficient screening threshold ∈ ->As shown in fig. 3. IMF components of each order and initial signals after the rotation speed of the power generation turbine is decomposedRoot mean square error>Root mean square error screening threshold- >As shown in fig. 4.
From which it can be derived that the formula (9) satisfies the requirementComponent sum->The remaining 5 IMF components do not meet the requirements.
For a pair ofAnd->The components are subjected to wavelet threshold noise reduction processing, the original signals are decomposed into a plurality of layers of approximate coefficients and detail coefficients by utilizing wavelet transformation, and noise information after wavelet transformation is mainly concentrated in the detail coefficients with smaller absolute values, so that the purpose of removing noise can be achieved by setting the detail coefficients with absolute values smaller than a specified threshold to be 0 and reconstructing the residual wavelet coefficients (namely the approximate coefficients obtained by decomposition and the reserved detail coefficients) back to the original signals by utilizing a wavelet inverse transformation method.
The threshold selection and noise reduction of the present embodiment includes the following steps:
estimating noise of each layer in IMF components meeting screening conditions by wavelet detail coefficients of each layer while utilizing coefficientsReducing the threshold value of the detail coefficient layer by layer, thereby preserving as much of the real signal contained in the high frequency component as possible; and carrying out threshold processing by using the composite threshold function to obtain the IMF component after wavelet noise reduction.
And (3) setting a noise threshold value of the j-th layer wavelet detail coefficient and a composite threshold function in the formula (11) by using the formula (10) to perform threshold processing. The signal to noise ratio and root mean square error of the denoised and non-denoised components are shown in table 1.
In the present embodiment, the noise reduction processing will be performedAnd->Component replacement +.>And->And after the components are accumulated with the remaining 5 IMF components which are not denoised, obtaining a historical operation time sequence of the rotating speed of the power generation turbine after denoise treatment. In this embodiment, the signal-to-noise ratio of the time series of the power turbine rotational speed after noise reduction relative to the time series before noise reduction is 81.4374, and the root mean square error is 0.0025. The residual component, the IMF component which does not meet the screening threshold and the IMF component which is subjected to noise reduction treatment are processed>Component sum->The components are reconstructed into a new IMF component matrix.
The same calculation method is used for calculating other 9 evaluation indexes to respectively obtain n IMF components subjected to wavelet threshold noise reduction treatmentAnd residual component->The new IMF component matrix is formed, and will not be described in detail herein.
The invention provides a non-stationary time sequence single parameter prediction model based on MREMD decomposition and wavelet threshold noise reduction, which improves the prediction precision of the model; in the definition of noise IMF components of each IMF component after MREMD decomposition, a subjective IMF component definition method for removing only the first-order IMF component to achieve the purpose of noise reduction is abandoned, and two indexes of a correlation coefficient and a root mean square error are introduced to achieve quantitative definition of the IMF component, so that the calculated error is reduced.
In this embodiment, feature extraction and dimension reduction are performed on the new IMF component matrix by using a singular value decomposition method to obtain a singular value component matrix, including the following steps:
s341, n pieces of processed wavelet threshold noise reductionComponent and residual component->Removing the mean value and then composing the matrix in the form of column vector>And calculates its covariance matrix +.>The expressions of the two matrices are respectively as follows:
(15)
in the method, in the process of the invention,a diagonal matrix consisting of singular values of a covariance matrix B;for n signals->In this embodiment n is 7; wherein->The characteristic values corresponding to covariance matrix B are arranged in descending order; the matrix U is a matrix formed by eigenvectors corresponding to eigenvalues of the covariance matrix B;
s342, combining the matrix U with the matrixVector corresponding to the zero eigenvalue is removed and reconstructed into K singular value components +.>The expression is as follows:
(16)
in the method, in the process of the invention,is the kth singular value component.
In the embodiment, a new IMF component matrix is formed by the rotating speed operation time sequence of the power generation turbine, singular value decomposition is carried out on the IMF component matrix by utilizing formulas (15) - (16), and 8 singular value components are obtained, namelyAs shown in fig. 5.
In this embodiment, the permutation entropy algorithm is used to calculate and obtain permutation entropy values of each singular value component matrix, including the following steps:
S351, reconstructing the phase by using a delay phase space reconstruction methodReconstructing the 1 st singular value component of the plurality of singular value components>Reconstruction into matrix->The expression is:
(17)
in the method, in the process of the invention,the number of sequences after reconstruction, τ is the delay time and m is the embedding dimension.
In the embodiment, during phase space reconstruction, the mutual information method is used to select the optimal delay time, and the mutual information value of each singular value component of the rotation speed of the power generation turbine is within the delay time range of 1-100 secondsThe variation is shown in fig. 6.
Calculating the optimal embedding dimension by adopting a pseudo-neighbor method, wherein the embedding dimension m is within the dimension variation range of 1-10, and the pseudo-neighbor rate of each singular value componentThe variation with embedding dimension m is shown in fig. 7.
S352, arranging each row of elements of the reconstruction matrix in an ascending order to obtain index values of each row of elements, and then forming index vectors by the index values of each row of elements to form an index matrix, wherein the arrangement modes of each row of elements in the index matrix are commonThe possibility that the number of occurrences of the b-th arrangement in the index matrix is counted as +.>Calculate the probability of its occurrence +.>The expression is:
(18)
wherein, the value range of b is (1, K), K is the number of singular value components.
S353, calculating the arrangement entropy of each singular value component, and arranging the entropyThe expression of (2) is as follows:
(19)
in this embodiment, the first extreme point of the time-dependent curve of the mutual information value in fig. 5 is taken as the optimal delay time of each singular value component of the rotation speed of the power generation turbineThe method comprises the steps of carrying out a first treatment on the surface of the Taking the embedding dimension when the pseudo-neighbor rate is less than 5% in FIG. 6 as the optimal embedding dimension of each singular value component +.>The method comprises the steps of carrying out a first treatment on the surface of the Then substituting it into formula (17), and calculating the arrangement entropy value PE of each component by formula (18) and formula (19), the specific values of which are shown in table 2 below:
in this embodiment, the step S36 includes the following steps:
s361, dividing the new IMF component matrix into two types according to the arrangement entropy value by using a K-means clustering algorithm, and randomly selecting the midpoint of the arrangement entropy value sequenceAs a category center, and select and +.>Furthest>The point is another category center, and the obtained average value of the two categories of points is used as the new category center to carry outClassifying, and randomly selecting a certain point in the sequence of permutation entropy values again>The calculation is performed as a class center until no more changes occur in both class centers. Finally, the lower class and the higher class of average entropy values are obtained. Wherein the remaining object points x in the sequence of calculated permutation entropy values are respectively associated with +. >Andthe distance of (2) can be +.>Instead, the expression is as follows:
(20)
in the above formula, i=1, …, K is the number of singular value components.
S362, screening out singular value components corresponding to the class with lower average entropy value in the two classes of entropy value points, and reconstructing the singular value components into trend items of the original signalsThe expression is as follows:
(21)
in the method, in the process of the invention, taking a matrix with the mean value expanded for IMF components corresponding to entropy points selected through cluster analysis, and replacing the mean value corresponding to unselected components with a zero vector; />For the screened singular value component matrix, the unselected components are replaced by zero vectors; />To->The matrix after the feature vector corresponding to the unselected component is replaced by the zero vector.
In the embodiment, K-means clustering is carried out on the arrangement entropy values of each component of the rotation speed of the power generation turbine by utilizing formulas (20) and (21), and finally selection is carried outAnd->And carrying out trend item reconstruction on the three components to obtain an original trend item, wherein an operation trend chart of the original trend item of the rotation speed of the power generation turbine is shown in fig. 8.
S363, for the trend itemPerforming stability test, and performing differential stabilization treatment on the time sequence with non-stable test result for f times to obtain trend item ++after differential stabilization treatment >The expression of which is shown below,
(22)
in the method, in the process of the invention,the trend item after the differential stabilization treatment is adopted; l is a delay operator; d is the delay order; />Representing the f-time difference.
S364, pairTrend item after stabilization treatmentPerforming AIC criterion grading, wherein the expression of the AIC criterion grading is as follows:
(23)
in the method, in the process of the invention,the maximum likelihood estimation value is the standard deviation of the model sequence residual error, N is the length of the parameter sequence, r is the number of independent parameters of the model, s is the order of the autoregressive coefficient polynomial, and q is the order of the moving average coefficient polynomial.
S365, obtaining an autoregressive coefficient polynomial order S and a moving average coefficient polynomial order q of the model through AIC test, and fitting an ARIMA (S, f, q) model, wherein the expression is as follows:
(24)
in the method, in the process of the invention,is an AR model coefficient; />Coefficients for an ARIMA model; />Is a white noise sequence; c is a constant; />Is->Is a component of the group.
In this example, the trend term shown in FIG. 8 was checked for stationarity, and Phillips-Perron stationarity was used in this example. The test statistic is suitable for the stationarity test of the heteroscedastic occasion and obeys the limit distribution of the corresponding ADF test statistic. The result of the stability test is a non-stable sequence, and the non-stable sequence is converted into a stable sequence after 1 differential calculation; then, performing a red pool information criterion (AIC) on the sequence subjected to differential stabilization to obtain an autoregressive coefficient polynomial order s and a moving average coefficient polynomial order q which are respectively 3 and 1; after 1 contrast score calculation and LB statistic test, an ARIMA (3, 1, 1) model as shown in formula (24) was fitted.
S366, carrying out LB (Ljung-Box) statistic test on the residual sequence of the ARIMA model, wherein the LB statistic test is represented by the following formula:
(25)
wherein N is an integer and is [1, c ]; c is the residual number;Ljung-Box statistics values for residual sequences; />An estimated value of the model residual error; />Is distributed in a chi-square mode.
In this embodiment, the operation trend of the power generation turbine is predicted at 550s, the operation trends from 400s to 550s are taken as training sets, the difference stabilization processing is performed on the original trend item by using the formula (22), the ARIMA (3, 1) model is obtained by fitting the formulas (23) and (24), the operation trend of 60 seconds in the future between 551-610 s is predicted, and the prediction result of the operation trend of the power generation turbine rotating speed is shown in fig. 9. And (3) performing LB (Ljung-Box) statistic test on the residual sequence of the ARIMA model by using a formula of LB statistic test of the formula (25) so as to meet the requirement.
The operation trend of the other 9 evaluation indexes can refer to the above steps, which are not described herein again, and finally the operation trend prediction sequences of the 10 evaluation indexes are obtained.
S4, converting the new IMF component matrix to obtain a parameter time sequence of each evaluation index, obtaining a weighted evaluation parameter sequence of each evaluation index by using an entropy weight method, and obtaining positive and negative ideal solutions of each moment point by the weighted evaluation parameter sequence.
In this embodiment, the step S4 includes the following steps:
s41, sequentially carrying out accumulation reconstruction on the IMF component subjected to noise reduction of each parameter, the component not subjected to noise reduction treatment and the residual component to form a noise-reduced parameter time sequence, converting the noise-reduced parameter time sequence into a fluctuation value sequence according to the setting value of each evaluation index, obtaining the entropy weight of the fluctuation value sequence by using an entropy weight method, and obtaining the weighted evaluation parameter sequence of each evaluation index by using the entropy weight and the fluctuation value sequence. The entropy weights of the evaluation indexes in this embodiment are respectively: the rotation speed of the power generation turbine is 0.0991, and the steam pressure after the nozzles 0.0906,1# to 3# of the vacuum degree of the condenser is 0.0860, 0.0877 and 0.1200 respectively; 0.1106 for the oil temperature, 0.1042 for the vapor pressure after the conditioning stage, 0.0857 for the vapor pressure after the nozzle, 0.1184 for the condenser water level, 0.0977 for the vapor seal pressure.
The invention combines the improved TOPSIS method with the existing B-type gray correlation degree, and converts the noise-reduced parameter time sequence into the fluctuation value sequence according to the setting value of each evaluation index, so that the positive and negative ideal solutions of each evaluation index of the equipment can be determined, and the problem that the negative ideal solutions cannot be obtained when the existing TOPSIS method is applied to the stability evaluation of the equipment containing non-stable parameters is solved.
S42, respectively selecting the point with the largest positive fluctuation value of the weighted evaluation parameter sequence of each evaluation index at the same moment as a positive ideal solution, and taking the point with the largest negative fluctuation value as a negative ideal solution, wherein the calculated formulas of the positive ideal solution and the negative ideal solution are as follows:
(12)
in the method, in the process of the invention,for positive understanding, add>Is a negative ideal solution; />,/>N represents the number of evaluation indexes, in this embodiment, N is 10, and the values of positive and negative ideal solutions of the weighted evaluation parameter sequences of the evaluation indexes are finally obtained, and the values of N are 940, and the results are shown in fig. 10 and 11, respectively.
S5, respectively calculating the distance relation value between the fluctuation value of each evaluation index and the positive and negative ideal solutions, wherein the expression is as follows:
(13)
wherein i=1, 2, …, N, the number of N in this embodiment is 940, and the initial time sequenceThe number of N is the same; j=1, 2, …, n, n represents the number of evaluation indexes, in this embodiment, n has a value of 10, < >>For the distance relation value between the fluctuation value of the j-th evaluation index and the positive ideal solution,/for the j-th evaluation index>And the relationship value of the distance between the fluctuation value of the j-th evaluation index and the negative ideal solution.
In the method, an improved distance relation value formula is provided when the positive and negative ideal solutions in the improved TOPSIS method are calculated, and the correlation between the fluctuation value of each evaluation index and the positive and negative ideal solutions is calculated through the formula; in addition, as the distance between the two is smaller, the evaluation index is unstable, the weight occupied by the evaluation index in the stability of the whole equipment is larger, namely the weight of the evaluation index and the corresponding distance relation of the evaluation index are in inverse proportion relation, the negative sign is added in the formula of the distance relation value, so that the real situation can be better reflected, and the obtained result is more accurate.
S6, respectively calculating gray correlation degrees of fluctuation values and positive and negative ideal solutions of all evaluation indexes by adopting a B-type gray correlation analysis method, wherein the method comprises the following steps of:
s61, calculating the overall displacement difference between the fluctuation value of the jth evaluation index at the ith moment and the positive ideal solutionThe expression is as follows:
(26);
n is the number of data points collected by each evaluation index, and in the embodiment, 940 data points are taken.
S62, respectively calculating the volatility of the j-th evaluation index and the first-order slope difference of the positive ideal solution, wherein the first-order slope difference of the j-th evaluation indexThe expression of (2) is as follows:
(27)
in the above-mentioned method, the step of,a weighted evaluation parameter indicating the (i+1) th time point in the jth evaluation index,/for the (j) th time point>A weighted evaluation parameter indicating an i-th point in time in the j-th evaluation index;
first order slope difference of positive ideal solutionThe expression of (2) is as follows:
(28)
in the above-mentioned method, the step of,representing the maximum forward fluctuation value corresponding to the (i+1) th moment point of n different evaluation indexes; />Representing the maximum forward fluctuation value corresponding to the n different evaluation indexes at the ith moment;
s63, then according to the first-order slope differenceComputing a first order overall slope difference from a first order slope difference of a positive ideal solutionThe expression is as follows:
(29)
Where i=2, 3, …, N, j=1, 2, …, n+1.
S64, calculating a second-order slope difference between the weighted evaluation parameter and the positive ideal solution. Wherein the second order slope difference of the evaluation parameter is weightedThe expression is as follows:
(30);
second order slope difference of positive ideal solutionThe expression of (2) is as follows:
(31)
s65, according to the second-order slope differenceAnd second order slope difference of ideal>Calculating the second order overall slope difference +.>The expression is as follows:
(32);
s66, calculating the change rate difference between the weighted evaluation parameter matrix sequence and the positive ideal solutionThe expression is as follows:
(33);
wherein i=3, 4, …, N; j=1, 2, …, n+1;is->Average value of N points in the sample evaluation indexes; />Is the average value of N points in the positive ideal solution;
s67, calculating the association degree of the volatility of the j-th evaluation index and the similarity degree of the positive ideal solution slopeThe expression is shown as follows:
(34);
s68, calculating the association degree of the volatility of the j-th evaluation index and the similarity of the positive ideal solution change rateThe expression is shown as follows:
(35)
s69, at the association degreeAnd->On the basis of (a) calculating gray correlation degree +_of the j-th evaluation index with respect to the positive ideal solution>The expression is shown as follows:
(36)
Wherein j=1, 2, …, n; omega is an adjusting coefficient, and the value range is (0, 1);
according to the method, the gray correlation degree between the fluctuation value of each evaluation index in the weighted evaluation parameter matrix and the negative ideal solution is calculated。
And S7, carrying out dimensionless treatment on the distance relation value and the gray association degree respectively, and fusing the gray association degree subjected to the dimensionless treatment and the distance relation value into a correlation weight of each evaluation index and the integral running stability state of the equipment.
In this embodiment, the distance relation value is subjected to dimensionless treatmentAnd->As shown in table 3 below. />
In this embodiment, the gray correlation degree is obtained after dimensionless treatmentAnd->As shown in table 4 below.
In this embodiment, after the dimensionless number is obtainedAnd->On the basis of the above, data fusion processing is carried out to obtain a correlation weight value of the volatility of the jth evaluation index and the overall operation stability of the equipment>The expression is as follows:
(37)
wherein, the value range of gamma is [0,1], and the magnitude of gamma determines the curve position relation and geometric shape relation pairIn the present embodiment, γ has a value of 0.5, and the distance relationship value and the weight value of the type B gray correlation degree are the same.
In this embodiment, the final weight of each evaluation index is shown below, where the correlation weight of the rotational speed of the power generation turbine is 0.3357, the correlation weight of the vacuum degree of the condenser is 0.2893,1# to 3# of the correlation weights of the steam pressure after the nozzle is 0.3174, 0.3186 and 0.3449, the correlation weight of the oil temperature is 0.5499, the correlation weight of the steam pressure after the adjustment stage is 0.3296, the correlation weight of the steam pressure after the nozzle is 0.3173, the correlation weight of the water level of the condenser is 0.3405, and the correlation weight of the seal pressure is 0.3199.
S8, on the basis of the predicted trend item, establishing an operation stability state scoring formula of each evaluation index, obtaining the scoring value of each evaluation index at each moment point in a predicted time period, and establishing an operation scoring formula of the predicted trend item of the power generation turbine according to the scoring value and the correlation weight of each evaluation index, so as to obtain the operation scoring of the predicted trend item of the whole equipment.
In this embodiment, the expressions of the 10 evaluation index running stability state scoring formulas can be expressed as follows:
(14)
wherein j=1, 2, …, n;for the weighted evaluation parameter of the jth evaluation index in the ith time point of the predicted interval, the predicted interval in this embodiment is 551-610 s, < - >And->And respectively representing the positive maximum fluctuation value and the negative maximum fluctuation value of each evaluation index in the operation range of the power generation turbine.
In the present embodiment, each evaluation index corresponds toAnd->The specific numerical value of the current operation evaluation index is determined according to actual operation experience, and is obtained through percentage conversion through setting values of all evaluation indexes, wherein the specific process is the same as that of the current operation evaluation index through percentage conversion into fluctuation value data.
As can be seen from the above-mentioned piecewise function (14), whenWhen the error of (2) is within a maximum fluctuation value of 0.5 times, the corresponding stability state score +.>100 points indicates that the evaluation index is very stable at the ith moment of the prediction interval; when->When the error of (2) is within the range of 0.5-1 times of the maximum fluctuation value, the corresponding stability state score +.>60-100 minutes, which indicates that the evaluation index is relatively stable at the ith moment point of the prediction interval; when->When the error of (2) is outside the maximum fluctuation range of 1 times, the corresponding stability state score +.>And the score is less than 60 minutes, which means that the evaluation index is unstable at the ith moment of the prediction interval, at the moment, the software system can give an alarm in time, and related technicians can adjust related parameters in time according to experience, so that the score is restored to more than 60 minutes. Stability status score of the individual evaluation index +. >The value of (2) is not less than 0.
In this embodiment, the stability status score using 10 evaluation indexesAnd the corresponding correlation weight +.>Establishing a real-time operation evaluation score formula of the power generation steam turbine, wherein the expression is as follows:
(38)
wherein Score is an operation evaluation Score of a predicted trend item of the power generation turbine,is the correlation weight of the fluctuation of the j-th evaluation index and the running stability of the power generation turbine, < ->And (3) scoring the running stability state of the jth evaluation index at the ith moment point of the prediction interval, wherein n is 10.
And S9, setting a total reference index, and comparing the operation evaluation score of the predicted trend item of the power generation turbine with the total reference index to judge whether adjustment of each evaluation index is needed, so that the operation evaluation score of the predicted trend item of the power generation turbine is higher than the total reference index finally.
In this embodiment, the step S9 includes the following steps:
s91, calculating operation evaluation scores of predicted trend items of the power generation steam turbine, and setting total reference indexes, wherein the total reference indexes are 80 minutes in the embodiment;
s92, comparing the operation evaluation score of the predicted trend item of the power generation turbine with the total standard index 80;
And S93, when the operation evaluation score of the predicted trend item of the power generation turbine is greater than or equal to the total reference index, the power generation turbine is well operated in a future predicted time period. In the present embodiment, as shown in fig. 12, the score of the operation evaluation score of the power turbine predicting the trend term between 551s to 610s is greater than 80 points, showing that the power turbine is operating well in the future predicted time period.
When the operation evaluation score of the predicted trend item of the power generation turbine is less than 80 minutes, the power generation turbine is judged to be unstable in operation, at the moment, the software system can give an alarm in time, and review of lower evaluation index scores is returned, and related technicians adjust related parameters in time according to experience, so that the overall real-time operation evaluation score of the power generation turbine is not less than 80 minutes, and the normal operation of the power generation turbine is ensured.
In this embodiment, the software system is Matlab, and each calculation and alarm function is implemented based on the software. Of course other similar software may be used. The operation condition of the power generation turbine in the operation process of the thermodynamic system can be intuitively observed according to the stability state evaluation score diagram of the power generation turbine in a period of time in the future, which is helpful for operation operators to observe the operation state of the condenser in time.
And the software can pre-warn the score of any evaluation index at a certain moment in a future prediction interval, and when the score is less than 60 minutes, relevant technicians timely adjust relevant parameters according to experience to restore the score to more than 60 minutes, so that the safe operation of the software is facilitated. Meanwhile, when the operation evaluation score of the predicted trend item of the power generation turbine is lower than 80, the software system can also give an alarm in time, so that the operation evaluation score of the final predicted trend item is not lower than 80, and the normal operation of the power generation turbine is ensured.
Claims (10)
1. The equipment stability state prediction and evaluation method based on the monitoring parameters is characterized by comprising the following steps of:
s1, analyzing an instability factor affecting the running state of equipment, selecting a plurality of evaluation indexes related to the running stability of the equipment, collecting data sets with the length of N of each evaluation index in a period of time, and forming each data set into an initial time sequence of the running of the evaluation index;
s2, decomposing the initial time sequence by adopting an MREMD decomposition method to obtain n IMF components and a final residual component;
S3, setting a screening threshold value to screen the IMF components, carrying out wavelet threshold value noise reduction treatment on the IMF components meeting the screening threshold value, recombining to obtain a new IMF component matrix, and processing the new IMF component matrix to obtain a final prediction trend item;
S4, converting the new IMF component matrix to obtain a parameter time sequence of each evaluation index, obtaining a weighted evaluation parameter sequence of each evaluation index by using an entropy weight method, and obtaining positive and negative ideal solutions of each moment point by the weighted evaluation parameter sequence;
s5, respectively calculating the distance relation value between the fluctuation value of each evaluation index and the positive and negative ideal solutions;
s6, respectively calculating gray correlation degrees of fluctuation values and positive and negative ideal solutions of all evaluation indexes by adopting a B-type gray correlation analysis method;
s7, carrying out dimensionless treatment on the distance relation value and the gray correlation value respectively, and fusing the gray correlation value and the distance relation value after the dimensionless treatment into a correlation weight of each evaluation index and the integral running stability state of the equipment;
s8, on the basis of the predicted trend item, establishing an operation stability state scoring formula of each evaluation index, obtaining scoring values of each evaluation index at each moment point in a predicted time period, and establishing an operation evaluation scoring formula of the predicted trend item of the whole equipment according to the scoring values and the correlation weights of each evaluation index;
s9, setting a total reference index, and comparing the operation evaluation score of the overall predicted trend item of the equipment with the total reference index to judge whether adjustment of each evaluation index is needed, so that the operation evaluation score of the overall predicted trend item of the equipment is higher than the total reference index finally.
2. The method for predicting and evaluating the stability status of a device based on monitoring parameters according to claim 1, wherein the step S2 comprises the following steps:
s21, utilizing an AR prediction model to perform initial time sequence on each evaluation indexPredicting and extending the left and right endpoints to obtain a new sequence so as to enable +.>Is between two extreme points of the new sequence after extension, wherein the AR prediction model expression is as follows:
;
in the above-mentioned method, the step of,the predicted point is the t moment; />The actual values of the parameters at the times t-1, t-2, & gtt-p, & lt/L, & gt>For p+1 real numbers, ">(t=p+1, p+2, & gtis, N) white noise sequence with zero mean;
s22, calculating the average value between adjacent extreme points in the new sequence, wherein the average value point obtained in the new sequence is the initial time sequenceMean point sequence>Wherein k is->The number of the mean value points is k+1, and the +.f is obtained through 3 times spline interpolation>Signal mean sequence>Make->Minus->Obtaining 1-order signal componentsThe expression is shown as follows:
;
s23, calculating the obtained signal componentAs the original signal, repeating the above steps S21 and S22 to perform iterative calculationTo the signal component after r iterations +.>And the termination condition is satisfied, and the judgment formula of the calculation termination condition is as follows:
;
In the above-mentioned method, the step of,for normalizing standard deviation, ++>And->For the signal component h after the r-1 th iteration respectively 1,r-1 The mean point sequence of (t) and the signal component after the r-th iteration +.>The standard deviation of the mean point sequence of (2), P is the conditional probability; definition according to sigma principle->For the initial time sequence->A sense value of the absolute value of the middle extreme point; />And->The mean point of the time series after the initial time series and the r iteration and +.>Ratio sequence of (2);said->、/>And->The formulas of (a) are respectively shown in the following formulas;
;
;
wherein:the z extreme point of the new sequence after extension; />For the initial time sequence->Is a k mean points of (2); />For time series component after the r-1 th iteration +.>Is a k mean points of (2);
s24, willAs an eigenmode function IMF of order 1 1 Output of->Minus->Get component->And will->The above steps S21 to S24 are repeated as an original signal until the residual signal becomes a monotone function or no new IMF component is separated, and the residual signal extraction process is as follows:
;
wherein n is the number of resolvable IMF components;as the final residual component;
s25, after the decomposition, the initial time sequenceCan be expressed as the n IMF components and the final residual component- >The sum is as follows:
。
3. the method for predicting and evaluating the stability state of equipment based on monitoring parameters according to claim 1, wherein the step S3 comprises the following steps:
s31, according to the IMF components and the initial time sequencePhase of (2)Coefficient of closure->Root mean square error->Obtaining the screening threshold value;
s32, screening the n IMF components by using the screening threshold, and dividing the IMF components into IMF components which do not meet the screening threshold and IMF components which meet the screening threshold;
s33, carrying out wavelet threshold denoising treatment on the components meeting the screening threshold, and recombining the residual components, the IMF components which do not meet the screening threshold and the IMF components subjected to denoising treatment into a new IMF component matrix;
s34, performing feature extraction and dimension reduction on the new IMF component matrix by using a singular value decomposition method to obtain a singular value component matrix;
s35, calculating and obtaining the arrangement entropy value of each singular value component matrix by using an arrangement entropy algorithm;
s36, dividing the new IMF component matrix into two types according to the arrangement entropy value by utilizing a K-means clustering algorithm, taking a type of component with a lower entropy value, reconstructing the type of component into an original trend item, carrying out differential stabilization treatment on the original trend item, obtaining an ARIMA model by fitting, and predicting the original trend item by the ARIMA model to obtain a final predicted trend item.
4. A method for predicting and evaluating the stability status of a device based on monitored parameters as claimed in claim 3, wherein said step S31 comprises the steps of:
s311, calculating the IMF components of each order and an initial time sequenceCorrelation coefficient of->Root mean square error->The expression is shown as follows:
;
in the middle ofAnd->Initial time series>And the average value of IMF components of each order, N is the sequence length of each parameter, and N is the number of IMF components;
s312, calculating a screening threshold value, wherein the screening threshold value comprises a correlation coefficient threshold valueAnd mean square error threshold->The expressions are shown below:
;
s313, selecting IMF components meeting the following conditions, and carrying out wavelet threshold denoising to obtain denoised IMF components, wherein the judgment formula of the screening conditions is as follows:
。
5. a method for predicting and evaluating the stability status of a device based on monitored parameters as claimed in claim 3, wherein said step S32 comprises the steps of:
s321, estimating the noise of each layer in the IMF components meeting the screening conditions through wavelet detail coefficients of each layer, and simultaneously utilizing the coefficientsThe threshold value of the detail coefficient is lowered layer by layer, and the selection criterion formula of the threshold value is as follows:
;
in the method, in the process of the invention,a noise threshold that is a layer j wavelet detail coefficient; / >The detail coefficient of the wavelet of the j-th layer; medium is an intermediate value function, and N is the number of data points collected by each evaluation index;
s322, performing threshold processing by using a composite threshold function to obtain an IMF component subjected to wavelet noise reduction, wherein the expression of the composite threshold function is as follows:
;
in the method, in the process of the invention,the j-th layer detail coefficient after the noise reduction treatment; sign is a sign function, a e [0,1] is an adjustment coefficient, a=0 is a hard threshold function, and a=1 is a soft threshold function.
6. The method for predicting and evaluating the stability status of a device based on monitoring parameters according to claim 1, wherein the step S4 comprises the following steps:
s41, sequentially carrying out accumulation reconstruction on an IMF component subjected to noise reduction of each parameter, a component not subjected to noise reduction treatment and a residual component to form a noise-reduced parameter time sequence, converting the noise-reduced parameter time sequence into a fluctuation value sequence according to a setting value of each evaluation index, obtaining an entropy weight value of the fluctuation value sequence by using an entropy weight method, and obtaining a weighted evaluation parameter sequence of each evaluation index by using the entropy weight value and the fluctuation value sequence;
s42, respectively selecting the point with the largest positive fluctuation value of the weighted evaluation parameter sequence of each evaluation index at the same moment as a positive ideal solution, and taking the point with the largest negative fluctuation value as a negative ideal solution, wherein the calculated formulas of the positive ideal solution and the negative ideal solution are as follows:
;
In the method, in the process of the invention,for positive ideal sequence,/o>Is a negative ideal solution sequence; />,/>N represents the number of evaluation indexes.
7. The method for predicting and evaluating the stability state of equipment based on monitoring parameters according to claim 6, wherein the expression of the distance relation between the fluctuation value of each evaluation index sequence and positive and negative ideal solutions is:
;
wherein i=1, 2, …, N; j=1, 2, …, n, n represents the number of evaluation indexes,for the distance relation value between the fluctuation value of the j-th evaluation index and the positive ideal solution,/for the j-th evaluation index>And the relationship value of the distance between the fluctuation value of the j-th evaluation index and the negative ideal solution.
8. The method for predicting and evaluating the stability state of equipment based on monitoring parameters according to claim 1, wherein the expression of the operation stability state scoring formula of each evaluation index is:
;
wherein j=1, 2, …, n;for the weighted evaluation parameter of the jth evaluation index in the ith point in time, +.>And->Respectively representing the positive maximum fluctuation value and the negative maximum fluctuation value of each evaluation index in the running range of the equipment.
9. The method for predicting and evaluating the stability status of a device based on monitored parameters as set forth in claim 8, wherein said step S9 comprises the steps of:
Calculating the running evaluation score of the overall predicted trend item of the equipment, and setting a total reference index;
comparing the operation evaluation score of the overall predicted trend item of the equipment with the total reference index;
when the operation evaluation score of the overall prediction trend item of the equipment is larger than or equal to the total reference index, the equipment is indicated to operate well in a future prediction time period;
when the operation evaluation score of the overall prediction trend item of the equipment is smaller than the total reference index, the operation instability of the equipment in a future prediction time period is indicated, and the evaluation index with the lower score is timely adjusted, so that the operation evaluation score of the overall prediction trend item of the equipment is greater than or equal to the total reference index.
10. The method for predicting and evaluating the stability state of equipment based on monitoring parameters according to any one of claims 1 to 9, wherein the method is dedicated to a power generation turbine of a steam power ship.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310415583.3A CN116611609B (en) | 2023-04-18 | 2023-04-18 | Equipment stability state prediction and evaluation method based on monitoring parameters |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310415583.3A CN116611609B (en) | 2023-04-18 | 2023-04-18 | Equipment stability state prediction and evaluation method based on monitoring parameters |
Publications (2)
Publication Number | Publication Date |
---|---|
CN116611609A true CN116611609A (en) | 2023-08-18 |
CN116611609B CN116611609B (en) | 2024-01-09 |
Family
ID=87677163
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310415583.3A Active CN116611609B (en) | 2023-04-18 | 2023-04-18 | Equipment stability state prediction and evaluation method based on monitoring parameters |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116611609B (en) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117150932A (en) * | 2023-10-30 | 2023-12-01 | 浙江威格泵业有限公司 | Energy efficiency evaluation method of variable-frequency booster shielding pump |
CN117311244A (en) * | 2023-11-28 | 2023-12-29 | 广州宝云信息科技有限公司 | Energy-saving regulation and control method and system based on equipment working condition prediction |
CN117368141A (en) * | 2023-12-07 | 2024-01-09 | 国检测试控股集团湖南华科科技有限公司 | Perchlorate wastewater concentration intelligent detection method based on artificial intelligence |
CN117871948A (en) * | 2024-03-11 | 2024-04-12 | 南方电网科学研究院有限责任公司 | Inter-harmonic detection method and system based on improved VMD algorithm |
CN118094114A (en) * | 2024-04-19 | 2024-05-28 | 四川轻化工大学 | Rotary machine fault diagnosis method for improving BWO-VMD-TOPSIS noise reduction |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104407601A (en) * | 2014-10-09 | 2015-03-11 | 讯腾数码科技(北京)有限公司 | Urban gas pipeline network pressure regulating equipment operation fault diagnosis method |
CN107908863A (en) * | 2017-11-14 | 2018-04-13 | 哈尔滨理工大学 | A kind of hydraulic turbine operating condition decision method based on EMD theories with HHT conversion |
CN109241849A (en) * | 2018-08-07 | 2019-01-18 | 浙江大学 | The method for diagnosing faults of empirical modal and wavelet decomposition and feature selecting towards intelligent power plant steam turbine host |
CN112949986A (en) * | 2021-02-01 | 2021-06-11 | 中国人民解放军海军工程大学 | Steam power system operation stability evaluation method and system, electronic equipment and storage medium |
US20230023931A1 (en) * | 2021-07-08 | 2023-01-26 | Zhejiang University | Hydraulic turbine cavitation acoustic signal identification method based on big data machine learning |
-
2023
- 2023-04-18 CN CN202310415583.3A patent/CN116611609B/en active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104407601A (en) * | 2014-10-09 | 2015-03-11 | 讯腾数码科技(北京)有限公司 | Urban gas pipeline network pressure regulating equipment operation fault diagnosis method |
CN107908863A (en) * | 2017-11-14 | 2018-04-13 | 哈尔滨理工大学 | A kind of hydraulic turbine operating condition decision method based on EMD theories with HHT conversion |
CN109241849A (en) * | 2018-08-07 | 2019-01-18 | 浙江大学 | The method for diagnosing faults of empirical modal and wavelet decomposition and feature selecting towards intelligent power plant steam turbine host |
CN112949986A (en) * | 2021-02-01 | 2021-06-11 | 中国人民解放军海军工程大学 | Steam power system operation stability evaluation method and system, electronic equipment and storage medium |
US20230023931A1 (en) * | 2021-07-08 | 2023-01-26 | Zhejiang University | Hydraulic turbine cavitation acoustic signal identification method based on big data machine learning |
Non-Patent Citations (2)
Title |
---|
卢绪祥;苏一鸣;吴家腾;李录平;: "基于EMD及灰色关联度的滑动轴承润滑状态故障诊断研究", 动力工程学报, vol. 36, no. 1, pages 42 - 47 * |
李济沅;邱玉婷;于竞哲;李莎;周浩;: "基于灰色关联度和TOPSIS的真空断路器运行状态综合评估", 高压电器, vol. 53, no. 3, pages 242 - 247 * |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117150932A (en) * | 2023-10-30 | 2023-12-01 | 浙江威格泵业有限公司 | Energy efficiency evaluation method of variable-frequency booster shielding pump |
CN117150932B (en) * | 2023-10-30 | 2024-01-26 | 浙江威格泵业有限公司 | Energy efficiency evaluation method of variable-frequency booster shielding pump |
CN117311244A (en) * | 2023-11-28 | 2023-12-29 | 广州宝云信息科技有限公司 | Energy-saving regulation and control method and system based on equipment working condition prediction |
CN117311244B (en) * | 2023-11-28 | 2024-02-13 | 广州宝云信息科技有限公司 | Energy-saving regulation and control method and system based on equipment working condition prediction |
CN117368141A (en) * | 2023-12-07 | 2024-01-09 | 国检测试控股集团湖南华科科技有限公司 | Perchlorate wastewater concentration intelligent detection method based on artificial intelligence |
CN117368141B (en) * | 2023-12-07 | 2024-03-01 | 国检测试控股集团湖南华科科技有限公司 | Perchlorate wastewater concentration intelligent detection method based on artificial intelligence |
CN117871948A (en) * | 2024-03-11 | 2024-04-12 | 南方电网科学研究院有限责任公司 | Inter-harmonic detection method and system based on improved VMD algorithm |
CN117871948B (en) * | 2024-03-11 | 2024-05-31 | 南方电网科学研究院有限责任公司 | Inter-harmonic detection method and system based on improved VMD algorithm |
CN118094114A (en) * | 2024-04-19 | 2024-05-28 | 四川轻化工大学 | Rotary machine fault diagnosis method for improving BWO-VMD-TOPSIS noise reduction |
Also Published As
Publication number | Publication date |
---|---|
CN116611609B (en) | 2024-01-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN116611609B (en) | Equipment stability state prediction and evaluation method based on monitoring parameters | |
CN104166787B (en) | A kind of aero-engine method for predicting residual useful life based on multistage information fusion | |
CN108334059B (en) | Particle swarm optimization-based BP neural network model fault diagnosis method | |
CN105467975B (en) | A kind of equipment fault diagnosis method | |
CN109740687B (en) | Fermentation process fault monitoring method based on DLAE | |
CN107895224B (en) | MKECA fermentation process fault monitoring method based on extended nuclear entropy load matrix | |
CN116226646B (en) | Method, system, equipment and medium for predicting health state and residual life of bearing | |
CN111324110A (en) | Fermentation process fault monitoring method based on multiple shrinkage automatic encoders | |
CN112884243A (en) | Air quality analysis and prediction method based on deep learning and Bayesian model | |
CN109948194B (en) | High-voltage circuit breaker mechanical defect integrated learning diagnosis method | |
CN113780420B (en) | GRU-GCN-based method for predicting concentration of dissolved gas in transformer oil | |
CN115421216A (en) | STL-ARIMA-NAR mixed model-based medium-and-long-term monthly rainfall forecasting method | |
CN114707577A (en) | Anomaly detection method and system based on self-confrontation variational self-encoder | |
CN113359435A (en) | Correction method for dynamic working condition data of thermal power generating unit | |
CN116797204A (en) | Primary air quantity fault early warning method for coal mill based on wavelet decomposition and reconstruction and TCN-GRU-Self-Attention | |
CN116933643A (en) | Intelligent data monitoring method based on partial robust M regression and multiple interpolation | |
Wu et al. | Comparing Deep Learning with Statistical Control Methods for Anomaly Detection:(072) | |
CN116306217A (en) | Track similarity residual life prediction method based on slow characteristic information gain ratio | |
CN115587541A (en) | Method for describing turbofan engine process characteristics by utilizing multi-feature fusion residual error network | |
CN115330085A (en) | Wind speed prediction method based on deep neural network and without future information leakage | |
CN107392379B (en) | Lorenz disturbance-based time series wind speed prediction method | |
CN111027612A (en) | Energy metering data feature reduction method and device based on weighted entropy FCM | |
CN116090800B (en) | Equipment stability real-time evaluation method based on monitoring parameters | |
CN114841000B (en) | Soft measurement modeling method based on modal common feature separation | |
CN109740111A (en) | Electric field value prediction algorithm over the ground |
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 |