CN109670625B - NOx emission concentration prediction method based on unscented Kalman filtering least square support vector machine - Google Patents
NOx emission concentration prediction method based on unscented Kalman filtering least square support vector machine Download PDFInfo
- Publication number
- CN109670625B CN109670625B CN201810992640.3A CN201810992640A CN109670625B CN 109670625 B CN109670625 B CN 109670625B CN 201810992640 A CN201810992640 A CN 201810992640A CN 109670625 B CN109670625 B CN 109670625B
- Authority
- CN
- China
- Prior art keywords
- model
- sample
- prediction
- support vector
- value
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 32
- 238000012843 least square support vector machine Methods 0.000 title claims abstract description 17
- 238000001914 filtration Methods 0.000 title claims abstract description 9
- 238000004364 calculation method Methods 0.000 claims abstract description 6
- 239000013598 vector Substances 0.000 claims abstract description 3
- 230000008859 change Effects 0.000 claims description 16
- 238000012549 training Methods 0.000 claims description 14
- 238000010606 normalization Methods 0.000 claims description 10
- 238000005070 sampling Methods 0.000 claims description 9
- 239000011159 matrix material Substances 0.000 claims description 8
- 230000009466 transformation Effects 0.000 claims description 5
- 150000001875 compounds Chemical class 0.000 claims description 4
- 229910052757 nitrogen Inorganic materials 0.000 claims description 4
- 238000005457 optimization Methods 0.000 claims description 4
- 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 claims description 2
- 239000006185 dispersion Substances 0.000 claims description 2
- 238000013507 mapping Methods 0.000 claims description 2
- 238000012216 screening Methods 0.000 claims description 2
- 238000012706 support-vector machine Methods 0.000 claims 3
- 230000007613 environmental effect Effects 0.000 abstract description 4
- 230000009286 beneficial effect Effects 0.000 abstract description 2
- 238000010531 catalytic reduction reaction Methods 0.000 abstract description 2
- MWUXSHHQAYIFBG-UHFFFAOYSA-N Nitric oxide Chemical compound O=[N] MWUXSHHQAYIFBG-UHFFFAOYSA-N 0.000 description 81
- QGZKDVFQNNGYKY-UHFFFAOYSA-N Ammonia Chemical compound N QGZKDVFQNNGYKY-UHFFFAOYSA-N 0.000 description 6
- 238000010586 diagram Methods 0.000 description 4
- 229910021529 ammonia Inorganic materials 0.000 description 3
- 101001095088 Homo sapiens Melanoma antigen preferentially expressed in tumors Proteins 0.000 description 2
- 102100037020 Melanoma antigen preferentially expressed in tumors Human genes 0.000 description 2
- 238000004422 calculation algorithm Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 239000000243 solution Substances 0.000 description 2
- 230000003044 adaptive effect Effects 0.000 description 1
- 238000013473 artificial intelligence Methods 0.000 description 1
- 239000003245 coal Substances 0.000 description 1
- 238000011217 control strategy Methods 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000002347 injection Methods 0.000 description 1
- 239000007924 injection Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
Images
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/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
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/24—Classification techniques
- G06F18/241—Classification techniques relating to the classification model, e.g. parametric or non-parametric approaches
- G06F18/2411—Classification techniques relating to the classification model, e.g. parametric or non-parametric approaches based on the proximity to a decision surface, e.g. support vector machines
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/10—Services
- G06Q50/26—Government or public services
Landscapes
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Tourism & Hospitality (AREA)
- Strategic Management (AREA)
- Human Resources & Organizations (AREA)
- Economics (AREA)
- General Business, Economics & Management (AREA)
- General Engineering & Computer Science (AREA)
- Evolutionary Computation (AREA)
- Data Mining & Analysis (AREA)
- Marketing (AREA)
- Development Economics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Evolutionary Biology (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Artificial Intelligence (AREA)
- Entrepreneurship & Innovation (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Life Sciences & Earth Sciences (AREA)
- Game Theory and Decision Science (AREA)
- Computer Hardware Design (AREA)
- Geometry (AREA)
- Educational Administration (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Feedback Control In General (AREA)
Abstract
The invention discloses a coal-fired unit NOx emission concentration prediction method based on an unscented Kalman filtering least square support vector machine. Determining input and output variables of a dynamic model of the denitration system through analysis of field data and theory, finishing prediction of NOx emission concentration at the current moment through offline calculation of parameters of the dynamic model of the denitration system, continuously updating a kernel parameter value sigma and model parameters alpha and b by adopting unscented Kalman filtering, updating a support vector sample, and predicting the NOx emission concentration at the next moment. The method provided by the invention can accurately predict the NOx emission concentration, is beneficial to further improving the regulation quality of the selective catalytic reduction denitration control system, can be used for judging whether field data is real and accurate, and provides a basis for supervision and law enforcement of environmental protection departments.
Description
Technical Field
The invention belongs to the technical field of thermal power environmental protection monitoring, and particularly relates to a NOx (nitrogen oxide) emission concentration prediction method based on an Unscented Kalman Filter Least square Support Vector Machine (UKF-LSSVM).
Background
With the strong promotion of energy-saving and emission-reducing policies, the emission of nitrogen oxides is more and more valued by environmental protection departments. A coal-fired boiler of a power plant is one of important sources of nitrogen oxides in the environment, and a Selective Catalytic Reduction (SCR) denitration technology is widely used in a thermal power plant at present. However, the denitration system is complex in interior, and coupling correlation among parameter variables is serious, so that the establishment of a corresponding physical model is very difficult. With the wide application of the distributed control system, a large number of operation parameters of the unit and auxiliary equipment related to the operation state are recorded, and with the development of various artificial intelligence and advanced algorithms and the successful application in industry, the research on system data modeling of the denitration system is developed.
The work of predicting and modeling the emission concentration of NOx emitted by a coal-fired unit is relatively less, the emission concentration of NOx emitted by the coal-fired unit is accurately predicted by using relevant operation parameters of the unit, whether field acquisition data is real and accurate can be judged, a basis is provided for supervision and law enforcement of environmental protection departments, and meanwhile, the regulation quality of an SCR control system can be greatly improved and ammonia escape is reduced by combining advanced control strategies such as prediction control, so that the method has important theoretical significance and practical value. Compared with an off-line least square support vector machine algorithm, the model can be updated in time according to the change of the object characteristics, and the method has adaptivity, but the parameter sigma needs to be set in advance, and the prediction performance of the least square support vector machine is greatly influenced.
Disclosure of Invention
The invention aims to solve the technical problem of providing a method for predicting NOx emission concentration based on an unscented Kalman filter least square support vector machine aiming at the defects of the prior art.
In order to achieve the technical purpose, the technical scheme adopted by the invention is as follows:
the method for predicting the NOx emission concentration based on the unscented Kalman filtering least square support vector machine comprises the following steps:
s1: selecting the concentration of NOx discharged by a coal-fired unit as prediction output, and determining a state parameter related to the concentration as an input variable of a model;
s2: sampling at the same frequency, and respectively removing gross errors and normalizing the original data by adopting an improved observed quantity change rate method and an upper and lower limit normalization method;
s3: selecting initial sample data according to the input and output structure of the dynamic model, calculating initial parameters alpha, b and sigma of the dynamic model of the denitration system in an off-line manner based on a least square support vector machine method, and obtaining a model prediction function at the same time;
s4: according to the model prediction function, the model prediction value of the current moment is input and calculated from the previous moment, and the model prediction error is calculated;
s5: judging whether the error exceeds the allowed maximum error, if not, the dynamic model is not corrected, the model parameters are not changed, and returning to the step S4 to continuously calculate the predicted value of the model at the next moment; if the error threshold value is exceeded, updating the model parameters sigma, alpha and b through the unscented transformation, updating the training sample data at the same time, and returning to the step S4 to continuously calculate the predicted value of the model at the next moment.
In the above step S1, the pearson correlation coefficient r is used to measure the correlation of the variables:
in the formula, xiIth sample value, y, for NOx concentration emitted by coal-fired unitiThe i-th sampled value of the state parameter related to the NOx concentration, N is the total number of samples,are the average values of the two sets of variables, respectively; the larger the value of r, the larger the correlation between the two variables, and a correlation threshold r is set0Selecting r value exceeding r0As input variables.
The above-described removal and normalization processing of the gross error in step S2 includes the steps of:
s2.1: the method for detecting the change rate of the improved observed quantity is adopted, whether the change quantity of the observed quantity at the current moment and the previous moment exceeds 3 times of standard deviation is judged to eliminate a coarse error, the change quantity between two adjacent moments is taken as an object, the 3 times of standard deviation criterion is improved, so that the threshold has self-adaptability, and the calculation formula is as follows:
Δxj=x(n-Nh+j)-x(n-Nh+j-1) (2)
wherein j is 2,3, … Nh;
Whereby an improved rate of change of the observed quantity is obtainedAnd judging a gross error by using a 3-time standard deviation criterion;
s2.2: and (3) adopting an upper and lower limit normalization method to normalize all sample data to be between [0 and 1], wherein the formula is as follows:
wherein x and x' are sample values before and after normalization, respectivelymin、xmaxMinimum and maximum values in the sample data.
In the step S3, the off-line calculation of the initial parameters of the dynamic model of the denitration system based on the least square support vector machine method and the obtaining of the model prediction function include the following steps:
for the training sample set T { (x)1,y1),…,(xN,yN) In which xi∈Rd,yiE, R, i is 1,2, …, N, N is the number of training samples, and an optimal decision function is constructed in a high-dimensional feature spaceThe problem is converted into solving the minimized structure risk, and the formula is as follows:
wherein w and b are model parameters, c is a penalty factor,as a non-line from the input space to the high-dimensional feature spaceSexual mapping, eiIs a prediction error;
the Lagrange function corresponding to this is:
wherein α ═ α1 α2 … αN]For Lagrange multipliers, according to the optimized Karush-Kuhn-Tucker (KKT) condition:
and (3) combining Lagrange function expression (8) and KKT optimization condition expression (9) to obtain a linear equation system of the optimization problem:
in the formula (I), the compound is shown in the specification,
y=[y1 y2 … yN]T (13)
in the formula (I), the compound is shown in the specification,taking the radial basis function as a kernel function:and (3) combining the formula (10) with the formula (14) to obtain an optimal decision function, and obtaining a predicted value at the current moment by inputting the predicted value at the previous moment, wherein the optimal decision function, namely a model prediction function, is as follows:
the above step S5 specifically includes the following steps:
for a random variable x ∈ RdThe expected value and covariance matrix are respectivelyP0WhereinP0Initially, a unit matrix is set, and a nonlinear transformation is performed by y ═ g (x) e RmY is a matrix x consisting of a pointk-1∈R2d+1The result is that the number of the first and second,
wherein psi ═ gamma2(d + k) -d is a screening parameter, gamma determines the dispersion degree of sigma points, the parameter k is 0, theta represents the distribution information of sampling points and is 2, and chiz,k-1Is xk-1Is in the z-th column of (1),is the z-th column of the square root of the matrix, WmkIs the desired weight, Wc, of the kth sample pointkThe weight of the k sampling point variance is obtained;
the method for estimating the parameters by using the unscented Kalman filtering comprises the following steps:
s5.1 initialization:
s5.2 σ Point Generation:
s5.3 time update:
χz,k|k-1=F(χz,k-1) (23)
yz,k|k-1=G(χz,k|k-1) (26)
s5.4 parameter updating:
obtaining next moment by iterative solutionAnd PkValue of, fromObtaining the model parameter values b, alpha and sigma at the next moment;
s5.5, updating a sample:
when the prediction error is larger, the sample needs to be updated online, the corresponding sample with the minimum alpha value is replaced by the current sample, and the current sample is set as { x+,y+In which x is+=[x1,x2,…xN]T,y+=[y1,y2…yN]TIt has been determined to remove { x from the set of support vectorsj,yjAnd f, the new training samples are:
xnew+=[x1,…xj-1,xN,xj+1…xN-1]T,ynew+=[y1,…yj-1,yN,yj+1…yN-1]T。
the invention has the following beneficial effects:
the unscented Kalman filter is used for updating the kernel parameter sigma value, so that the adaptivity of the online least square support vector machine is further improved, the prediction precision and the generalization capability are better, and an ideal effect can be achieved when fewer training samples are available.
Drawings
FIG. 1 is a flow chart of the present invention;
FIG. 2 is a sample set of initial training samples of a prediction model according to an embodiment of the present invention;
FIG. 3 is a Pearson correlation coefficient plot of NOx emission concentration for a coal fired unit in an embodiment of the present invention;
FIG. 4 is a diagram of a model for predicting NOx concentration emitted from a coal-fired unit according to an embodiment of the present invention;
FIG. 5 is a diagram showing the prediction results of the LSSVM model in embodiment 1 of the present invention;
FIG. 6 is a diagram of the predicted result of the UKF-LSSVM model in embodiment 1 of the present invention;
FIG. 7 is a diagram of the prediction result of the UKF-LSSVM model in embodiment 2 of the present invention;
FIG. 8 is a graph of variation of the kernel parameter σ with load in the UKF-LSSVM model in embodiment 2 of the present invention.
Detailed Description
Embodiments of the present invention are described in further detail below with reference to the accompanying drawings.
Example 1: the method for predicting the NOx emission concentration based on the unscented Kalman filtering least square support vector machine comprises the following specific steps as shown in figure 1:
step 1: a certain 300MW coal-fired power plant boiler is taken as a research object, and an SCR denitration system device which is relatively mature in technology and relatively wide in application is adopted. The selected data time span is one month, the sampling time is 5s, and the data comprises steady-state data in a normal and stable running state of the system and unsteady-state data in variable load. And selecting 5000 sample points under the condition of load change for modeling simulation. Through repeated experiments, under the condition of considering both the operation speed and the model precision, the number of training samples is selected to be 10, the first 10 samples of the initial sample set are taken as training samples, and the prediction model initial training sample set is shown in fig. 2.
Wherein: pe(t-1) is the unit load; xin1(t-1)、Xin2(t-1) NOx concentration at inlets of the SCR reactors at two sides; xN1(t-1)、XN2(t-1) the amount of ammonia sprayed on both sides; xou1(t-1)、Xou2(t-1) NOx concentration at inlets of the SCR reactors at two sides; y isNOx(t-2)、YNOx(t-1)、YNOxAnd (t) is the concentration of NOx discharged by the coal-fired unit at different sampling moments.
The degree of closeness of the relationship between the variables is checked according to the Pearson correlation coefficient, and the formula is as follows:
in the formula, xiIth sample value, y, for NOx concentration emitted by coal-fired unitiThe i-th sampled value of the state parameter associated with the NOx concentration, N is the total number of samples,are the average values of the two sets of variables, respectively;
pearson correlation coefficients between the candidate variables and the NOx emission concentration of the coal-fired unit are respectively calculated, and specific data are shown in a graph 3. And selecting a variable with the Pearson correlation coefficient | r | ≧ 0.15 as an auxiliary variable of the NOx emission concentration prediction model. And finally, selecting the concentration of NOx at the outlets of the SCR reactors at two sides, the ammonia injection amount of the mixers at two sides, the concentration of NOx at the inlets of the reactors at two sides and the load of the unit as input variables of a prediction model, and taking the concentration of NOx discharged by the coal-fired unit as output variables of the model. The difference from the steady-state modeling is that the dynamic modeling adds the order of the variables to the input and output, and the structure of the specific denitration system dynamic model is shown in fig. 4.
Step 2: the sampling frequency of all data is the same, the coarse error of the original data is removed by adopting an improved observation quantity change rate method, and the data is normalized;
s2.1: and (3) eliminating a coarse error by adopting an improved observation quantity change rate detection method and judging whether the observation quantity change quantity of the current moment and the previous moment exceeds 3 times of a standard deviation (3 sigma). Aiming at the variation between two adjacent moments, the 3 sigma criterion is improved to make the threshold adaptive, and the calculation formula is as follows:
Δxi=x(n-Nh+i)-x(n-Nh+i-1) (2)
wherein i is 2,3, … Nh;
Whereby an improved observation rate of change is determinedAnd used for the 3 sigma criterion for gross error determination.
S2.2: and (3) adopting an upper and lower limit normalization method to normalize all sample data to be between [0,1 ]. The formula is as follows:
wherein x and x' are sample values before and after normalization, respectivelymin、xmaxMinimum and maximum values in the sample data.
S3: selecting initial sample data according to the input and output structure of the dynamic model, and calculating the initial parameters of the dynamic model of the denitration system in an off-line manner based on a least square support vector machine method;
and step 3: according to an input and output structure of a dynamic model, selecting the first 10 sample sets as initial sample data, calculating initial parameters of the dynamic model of the denitration system in an off-line mode based on a least square support vector machine method, wherein the initial parameters are respectively alpha [ -2.55-2.171.351.421.22-1.002.94-2.22-1.192.21 ] and b ═ 0.411, and the sigma initial value is set to be 3;
and 4, step 4: inputting and calculating a model prediction value of the current moment from the previous moment, and calculating a model prediction error;
and 5: judging whether the error exceeds the allowed maximum error, if not, the dynamic model does not need to be corrected, the model parameters are kept unchanged, and returning to the step S4 to continuously calculate the predicted value of the model at the next moment; and if the error threshold is exceeded, updating the model parameters sigma, alpha and b through unscented Kalman, updating the training sample data, returning to the step S4 to continuously calculate the model predicted value at the next moment until the concentration of NOx discharged by the coal-fired unit at 5000 moments is completely predicted.
The model accuracy evaluation standard selects a root mean square error and an average absolute percentage error to evaluate the accuracy degree of the UKF-LSSVM and LSSVM models, and the calculation formula is as follows:
As shown in fig. 6, it can be seen from a comparison between fig. 5 and 6 that MAPE and RMSE of the UKF-LSSVM method are 0.17% and 0.0707mg/m, respectively-3The prediction model has high precision, and the two groups of indexes are respectively smaller than the indexes of the LSSVM, MAPE is 3.83%, and RMSE is 2.4166mg/m-3The UKF-LSSVM can more accurately predict the concentration value of the NOx discharged by the coal-fired unit compared with the LSSVM, has better prediction precision and stronger self-adaptive capacity,
example 2: the present example is different from example 1 in that the initial value of σ is set to 0.1, and other steps and parameters are the same as those of example 1. When the initial value of sigma is set to 0.1, the prediction result of the UKF-LSSVM model is shown in FIG. 7, the variation of the kernel parameter sigma in the UKF-LSSVM model along with the load is shown in FIG. 8, and as can be seen from FIG. 7 and FIG. 8, the UKF-LSSVM model still obtains a good prediction result for any given initial value of the kernel parameter, the model precision is not influenced by the initial kernel parameter sigma value, so that the difficulty in selecting sigma is greatly reduced, and the method has an important effect on facilitating the practical application of the method.
The above is only a preferred embodiment of the present invention, and the protection scope of the present invention is not limited to the above-mentioned embodiments, and all technical solutions belonging to the idea of the present invention belong to the protection scope of the present invention. It should be noted that modifications and embellishments within the scope of the invention may be made by those skilled in the art without departing from the principle of the invention.
Claims (4)
1. The prediction method of the NOx emission concentration based on the unscented Kalman filtering least square support vector machine is characterized by comprising the following steps: the method comprises the following steps:
step S1: selecting the concentration of NOx discharged by a coal-fired unit as prediction output, and determining a state parameter related to the concentration as an input variable of a model;
step S2: sampling at the same frequency, and respectively removing gross errors and normalizing the original data by adopting an improved observed quantity change rate method and an upper and lower limit normalization method;
step S3: selecting initial sample data according to the input and output structure of the dynamic model, calculating initial parameters alpha, b and sigma of the dynamic model of the denitration system in an off-line manner based on a least square support vector machine method, and obtaining a model prediction function at the same time;
step S4: according to the model prediction function, the model prediction value of the current moment is input and calculated from the previous moment, and the model prediction error is calculated;
step S5: judging whether the error exceeds the allowed maximum error, if not, the dynamic model is not corrected, the model parameters are not changed, and returning to the step S4 to continuously calculate the predicted value of the model at the next moment; if the error threshold is exceeded, updating the model parameters sigma, alpha and b through unscented transformation, updating training sample data at the same time, and returning to the step S4 to continuously calculate the model predicted value at the next moment;
the step S3 of calculating the initial parameters of the dynamic model of the denitration system in an off-line manner based on the least square support vector machine method and obtaining the model prediction function at the same time comprises the following steps:
for the training sample set T { (x)1,y1),…,(xN,yN) In which x isi∈Rd,yiE, R, i is 1,2, …, N, N is the number of training samples, and an optimal decision function is constructed in a high-dimensional feature spaceThe problem is converted into solving the minimized structure risk, and the formula is as follows:
wherein w and b are model parameters, c is a penalty factor,for non-linear mapping from the input space to the high-dimensional feature space, eiIs a prediction error;
the Lagrange function corresponding to this is:
wherein α ═ α1 α2…αN]For Lagrange multipliers, according to the optimized Karush-Kuhn-Tucker (KKT) condition:
and (3) combining Lagrange function expression (8) and KKT optimization condition expression (9) to obtain a linear equation system of the optimization problem:
in the formula (I), the compound is shown in the specification,
y=[y1 y2…yN]T (13)
in the formula (I), the compound is shown in the specification,taking the radial basis function as a kernel function:combining the formula (10) and the formula (14) to obtain an optimal decision function, and obtaining a predicted value of the current time from the input of the previous time, wherein the optimal decision function, namely a model prediction function, is as follows:
2. the unscented kalman filter least squares support vector machine-based NOx emission concentration prediction method according to claim 1, characterized by: in step S1, the pearson correlation coefficient r is used to measure the correlation of the variables:
in the formula, xiIth sample value, y, for NOx concentration emitted by coal-fired unitiThe i-th sampled value of the state parameter related to the NOx concentration, N is the total number of samples,are the average values of the two sets of variables, respectively; the larger the value of r, the larger the correlation between the two variables, and a correlation threshold r is set0Selecting r value exceeding r0As input variables.
3. The method of claim 1 for predicting NOx emission concentration based on unscented Kalman filter least squares support vector machines, characterized by: the coarse error removal and normalization process of step S2 includes the following steps:
s2.1: the method for detecting the change rate of the improved observed quantity is adopted, whether the change quantity of the observed quantity at the current moment and the previous moment exceeds 3 times of standard deviation is judged to eliminate a coarse error, the change quantity between two adjacent moments is taken as an object, the 3 times of standard deviation criterion is improved, so that the threshold has self-adaptability, and the calculation formula is as follows:
Δxj=x(n-Nh+j)-x(n-Nh+j-1) (2)
wherein j is 2,3, … Nh;
Whereby an improved observation rate of change is determinedAnd judging a gross error by using a 3-time standard deviation criterion;
s2.2: and (3) adopting an upper and lower limit normalization method to normalize all sample data to be between [0 and 1], wherein the formula is as follows:
wherein x and x' are respectively the sample values before and after normalization, xmin、xmaxMinimum and maximum values in the sample data.
4. The method of predicting NOx emission concentration based on unscented Kalman filter least squares support vector machine according to claim 1, characterized in that: the specific steps of the traceless transformation in step S5 are:
for a random variable x ∈ RdThe expected value and covariance matrix are respectivelyP0WhereinP0Initially, a unit matrix is set, and a nonlinear transformation is performed by y ═ g (x) e RmY is a matrix x consisting of a pointk-1∈R2d+1The result is that,
wherein psi ═ gamma2(d + k) -d is a screening parameter, gamma determines the dispersion degree of sigma points, the parameter k is 0, theta represents the distribution information of sampling points and is 2, and chiz,k-1Is xk-1Is in the z-th column of (1),is the z-th column of the square root of the matrix, WmkIs the desired weight, Wc, of the kth sample pointkThe weight of the k sampling point variance is obtained;
the method for estimating the parameters by using the unscented Kalman filtering comprises the following steps:
s5.1 initialization:
s5.2. Generation of σ Point:
s5.3 time update:
χz,k|k-1=F(χz,k-1) (23)
yz,k|k-1=G(χz,k|k-1) (26)
s5.4 parameter updating:
obtaining next moment by iterative solutionAnd PkValue of, fromObtaining the model parameter values b, alpha and sigma at the next moment;
s5.5, updating a sample:
when the prediction error is larger, the sample needs to be updated on line, the corresponding sample with the minimum alpha value is replaced by the current sample, and the current sample is set as { x+,y+In which x+=[x1,x2,…xN]T,y+=[y1,y2…yN]TIt has been determined to remove { x from the set of support vectorsj,yjAnd f, the new training samples are:
xnew+=[x1,…xj-1,xN,xj+1…xN-1]T,ynew+=[y1,…yj-1,yN,yj+1…yN-1]T。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810992640.3A CN109670625B (en) | 2018-08-28 | 2018-08-28 | NOx emission concentration prediction method based on unscented Kalman filtering least square support vector machine |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810992640.3A CN109670625B (en) | 2018-08-28 | 2018-08-28 | NOx emission concentration prediction method based on unscented Kalman filtering least square support vector machine |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109670625A CN109670625A (en) | 2019-04-23 |
CN109670625B true CN109670625B (en) | 2022-06-17 |
Family
ID=66141627
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810992640.3A Active CN109670625B (en) | 2018-08-28 | 2018-08-28 | NOx emission concentration prediction method based on unscented Kalman filtering least square support vector machine |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109670625B (en) |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110935312B (en) * | 2019-12-16 | 2022-08-30 | 广州珠江电力有限公司 | Dynamic monitoring device and dynamic monitoring method for SCR flue gas denitration system |
CN112506055A (en) * | 2020-12-03 | 2021-03-16 | 国网山东省电力公司电力科学研究院 | Thermal power plant NOx emission optimization control method and system |
CN112697977A (en) * | 2020-12-24 | 2021-04-23 | 重庆大唐国际石柱发电有限责任公司 | Thermal power station boiler flue gas NOx index prediction method |
CN113009890B (en) * | 2021-01-25 | 2021-12-31 | 南京工业大学 | Thermal power plant SCR flue gas denitration control system optimization control method based on BMO-LSSVM |
CN115047894B (en) * | 2022-04-14 | 2023-09-15 | 中国民用航空总局第二研究所 | Unmanned aerial vehicle track measuring and calculating method, electronic equipment and storage medium |
CN114935627A (en) * | 2022-05-07 | 2022-08-23 | 金川集团信息与自动化工程有限公司 | Method for predicting purity of hydrogen chloride at outlet of synthesis furnace |
CN116884516A (en) * | 2023-09-05 | 2023-10-13 | 河南科技学院 | PM2.5 concentration prediction method based on SVM-UKF data fusion |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104091216A (en) * | 2014-07-29 | 2014-10-08 | 吉林大学 | Traffic information predication method based on fruit fly optimization least-squares support vector machine |
CN106096646A (en) * | 2016-06-07 | 2016-11-09 | 衢州学院 | A kind of support vector regression model selection method |
CN106680428A (en) * | 2016-12-19 | 2017-05-17 | 华北电力大学(保定) | Soft measuring method for denitration control system |
-
2018
- 2018-08-28 CN CN201810992640.3A patent/CN109670625B/en active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104091216A (en) * | 2014-07-29 | 2014-10-08 | 吉林大学 | Traffic information predication method based on fruit fly optimization least-squares support vector machine |
CN106096646A (en) * | 2016-06-07 | 2016-11-09 | 衢州学院 | A kind of support vector regression model selection method |
CN106680428A (en) * | 2016-12-19 | 2017-05-17 | 华北电力大学(保定) | Soft measuring method for denitration control system |
Also Published As
Publication number | Publication date |
---|---|
CN109670625A (en) | 2019-04-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109670625B (en) | NOx emission concentration prediction method based on unscented Kalman filtering least square support vector machine | |
Liu et al. | Short-term prognostics of PEM fuel cells: A comparative and improvement study | |
CN110163429B (en) | Short-term load prediction method based on similarity day optimization screening | |
CN107122526B (en) | Experimental segment Mach number modeling method based on abnormal feature subset integration algorithm | |
CN110782067B (en) | Sintering end point fluctuation range prediction method based on fuzzy information granulation | |
CN107403196B (en) | Method for predicting concentration of butane at bottom of debutanizer by instant learning modeling based on spectral clustering analysis | |
CN112287605B (en) | Power flow checking method based on graph convolution network acceleration | |
CN112149714A (en) | Method for determining energy efficiency characteristic index reference value of coal-electric unit based on data mining | |
CN108492026B (en) | Soft measurement method based on integrated orthogonal component optimization regression analysis | |
CN110929835B (en) | Novel silicon carbide-based aviation power converter fault diagnosis method and system | |
CN116305985A (en) | Local intelligent ventilation method based on multi-sensor data fusion | |
CN111639304A (en) | CSTR fault positioning method based on Xgboost regression model | |
CN117592592A (en) | VMD-SSA-LSTM-based power transmission line icing thickness prediction method | |
CN112800672B (en) | Evaluation method, system, medium and electronic equipment for boiler fouling coefficient | |
CN116522065A (en) | Coal mill health degree assessment method based on deep learning | |
CN112416913B (en) | GWO-BP algorithm-based aircraft fuel system state missing value supplementing method | |
CN108182337B (en) | Natural gas purification process modeling method based on ScMiUKFNN algorithm | |
CN112100902A (en) | Lithium ion battery service life prediction method based on stream data | |
Ding et al. | Real-time boiler control optimization with machine learning | |
CN116595883B (en) | Real-time online system state correction method for numerical reactor | |
CN113191082B (en) | Model parameter obtaining method and system based on machine learning and readable medium | |
Yinghua et al. | Research of Quality Prediction Based on Extreme Learning Machine | |
Xu et al. | Adaptive control of a class of nonlinear discrete-time systems using support vector machine | |
Guan et al. | Interval Optimal Controller Design for Uncertain Systems Based on Interval Neural Network | |
Zhang et al. | Improved GM (1, N) model for equipment development cost prediction |
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 |