CN101498458B - Model updating method for on-line boiler combustion optimization - Google Patents
Model updating method for on-line boiler combustion optimization Download PDFInfo
- Publication number
- CN101498458B CN101498458B CN200910096410XA CN200910096410A CN101498458B CN 101498458 B CN101498458 B CN 101498458B CN 200910096410X A CN200910096410X A CN 200910096410XA CN 200910096410 A CN200910096410 A CN 200910096410A CN 101498458 B CN101498458 B CN 101498458B
- Authority
- CN
- China
- Prior art keywords
- model
- alpha
- centerdot
- sigma
- gamma
- 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.)
- Expired - Fee Related
Links
Classifications
-
- 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
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
- Y02P90/00—Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
- Y02P90/02—Total factory control, e.g. smart factories, flexible manufacturing systems [FMS] or integrated manufacturing systems [IMS]
Landscapes
- Feedback Control In General (AREA)
Abstract
The invention relates to a method for updating a boiler on-line combustion optimization model. The method utilizes the data exceeding the predictive error limit of an exit model to build a new model, then utilizes combustion data of a boiler and an optimization algorithm to find an optimum new model and the weight coefficient of the existing model, and utilizes an optimum weight coefficient to combine the new model with the existing model so as to predict the working conditions of a new boiler for realizing model update. The invention overcomes the defects that the traditional update method gives up the existing model and can not utilize the study result of the existing model to fully utilize the study result of the existing model, thereby greatly reducing the amount and the time of the calculation work of model update.
Description
Technical field
The invention belongs to information control technology field, relate to the incremental learning technology, particularly relate to a kind of model update method of on-line boiler combustion optimization.
Background technology
The burning optimization of station boiler is the important technical of energy-saving and emission-reduction, and its target is by adjusting the running status that operational factors such as boiler wind speed adjustment, coal supply obtain high efficiency or low pollution emission.The collocation of a lot of operational factors of boiler such as air distribution situation, coal supply situation and oxygen amount etc. all has direct influence to the boiler combustion state, and the configuration of operating parameters such as different air distributions, coal supply can directly cause the situation of the discharging of different boiler efficiencies and dusty gas.For given station boiler, under certain loading condiction, there is a kind of operating parameter allocation plan of optimum for the index of different fired states, can make the fired state optimization.Yet, very complicated coupled relation is arranged between the operating parameter of boiler, the configuration that find optimum operating parameter is very difficult.Along with continuous progress in science and technology, the boiler operatiopn automaticity improves constantly, but the boiler combustion optimization problem never well is resolved.
The research focus of boiler combustion optimization is by data mining at present, in a large amount of different actual operation parameters combinations, the method of applied for machines study, excavate the relational model between operational factor and boiler combustion index, utilize again and optimize algorithm carries out boiler in conjunction with the model of excavating online burning optimization.This method uses manpower and material resources sparingly, and can find the parameter configuration more excellent than artificial experiment, but because the characteristic of boiler plant changes to some extent along with the growth meeting of service time, and use coal and also can change sometimes, how guaranteeing therefore that model can upgrade fast and efficiently to adapt to news becomes the bottleneck problem of this method.This problem and modeling method, sample data are chosen and update strategy etc. all has much relations.
Summary of the invention
The objective of the invention is at the bottleneck problem in the boiler combustion optimization, propose a kind of model update method of taking into account historical learning outcome of model and new situation of change.
The present invention utilizes to exceed the data that have the model prediction limit of error and set up new model, and then utilize the boiler combustion data and optimize algorithm, seek the optimum new model and the weight coefficient of existing model, utilize optimum weight coefficient that new model and existing model are combined, jointly new boiler working condition is predicted that implementation model upgrades.This method has overcome and will have model in traditional update method and abandon fully, can not utilize the shortcoming of the learning outcome of existing model, has made full use of the learning outcome of existing model, has shortened the amount of calculation and the time of model modification greatly.
Technical scheme of the present invention is to make sample by the data of utilizing existing model prediction to exceed the limit of error, set up new model, and with new model and means such as existing model combines, establish a kind of model update method of on-line boiler combustion optimization, utilized the renewal of implementation model fast and efficiently of this method.
The step of the inventive method comprises:
Step (1) is set up the prediction error database that has had model.According to concrete boiler combustion situation with to the requirement of model prediction precision, the permission predicated error limit δ of setting model, when image data, error between judgment models predicted value and the actual motion value and the size that allows predicated error limit δ, if predicated error is greater than δ, promptly | V
c-V
s|>δ, wherein V
cBe model predication value, V
sBe actual operating data, the data that then will transfinite deposit in the prediction error database, in order to the usefulness of model modification.
Step (2) is set up new model.When master mould need upgrade, utilize data in the prediction error database as training sample, sample can be expressed as { x
i, y
i}
I=1 N, x wherein
iRepresent the boiler operating parameter vector of i group, y as the input data
iRepresent the parameter of i group as the sign boiler combustion state of output parameter, adopt algorithm of support vector machine to set up new model, kernel function is elected RBF as:
φ (x) is a mapping function, establishes the object function of being asked to be: f (x
i)=w φ (x
i)+b, f (x
i) being the boiler combustion index prediction value of model output, w is the weight coefficient vector, b is an intercept.Introduce relaxation factor ξ
* i〉=0 and ξ
i〉=0 and allow error of fitting ε, model can be by in constraint:
Obtain, wherein constant c>0 is penalty coefficient.This minimization problem is a convex quadratic programming problem,
Introduce Lagrangian:
At the saddle point place, function L is about w, b, ξ
i, ξ
i *Minimal point, also be α
i, α
i *, γ
i, γ
i *Maximal point, minimization problem are converted into the maximization problems of asking its dual problem.
Lagrangian L is about w at the saddle point place, b, ξ
i, ξ
i *Minimal point:
Can get the dual function of Lagrangian:
At this moment,
According to Ku En-Plutarch (KKT) conditional theorem, have following formula to set up at saddle point:
By following formula as seen,
α
iAnd α
i *Can not be non-zero simultaneously, can get:
Can obtain b from following formula, obtain model.Empty the prediction error database, in order to upgrade the usefulness of image data next time.
Step (3) is determined the weight of new model and existing model.Gather the different operating conditions of up-to-date boiler data down as test samples, use the master mould prediction and predict weighted-average method, check data is predicted, i.e. y with new model
i=α D
n+ β D
o, y wherein
iBe i group Optimization of operating target prediction value, D
nBe new model predicted value, D
oBe original model predication value, α is a new model predicted value weight coefficient, β is the prediction weight coefficient of original model, and alpha+beta=1, definite employing particle cluster algorithm iteration optimizing of α and β is determined, definition particle cluster algorithm initial position Z vector respectively tie up component, be respectively new model weight and original model weight beta, object function is: the minJ=∑ | y
i-Y
i|, Y wherein
iIt is the actual motion value of i group operating mode.When J had reached minimum, reaches setting value or finished iterations, training was finished, and obtains the weight coefficient of new model and existing model.
Step (4) combines original model with the weight proportion of new model by optimum, constitute new model, i.e. D=α D
n+ β D
o, wherein D realizes the renewal of boiler combustion optimization model for the model after upgrading.
The model update method that the present invention proposes has made full use of the useful information that existing model comprised, significantly reduced the workload of model modification, improve the efficient of model modification, satisfied the actual requirement of boiler combustion on-line optimization, guaranteed the real-time and the accuracy of boiler combustion optimization.
The specific embodiment
A kind of model update method of on-line boiler combustion optimization, concrete steps are:
(1) sets up the prediction error database that has had model.According to concrete boiler combustion situation with to the requirement of model prediction precision, the permission predicated error limit δ of setting model, when image data, error between judgment models predicted value and the actual motion value and the size that allows predicated error limit δ, if predicated error is greater than δ, promptly | V
c-V
s|>δ, wherein V
cBe model predication value, V
sBe actual operating data, the data that then will transfinite deposit in the prediction error database, in order to the usefulness of model modification.
(2) set up new model.When master mould need upgrade, utilize data in the prediction error database as training sample, sample can be expressed as { x
i, y
i}
I=1 N, x wherein
iRepresent the boiler operating parameter vector of i group, y as the input data
iRepresent the parameter of i group as the sign boiler combustion state of output parameter, adopt algorithm of support vector machine to set up new model, kernel function is elected RBF as:
φ (x) is a mapping function, establishes the object function of being asked to be: f (x
i)=w φ (x
i)+b, f (x
i) being the boiler combustion index prediction value of model output, w is the weight coefficient vector, b is an intercept.Introduce relaxation factor ξ
* i〉=0 and ξ
i〉=0 and allow error of fitting ε, model can be by in constraint:
Obtain, wherein constant c>0 is penalty coefficient.This minimization problem is a convex quadratic programming problem, introduces Lagrangian:
Wherein:
Be Lagrange's multiplier.
At the saddle point place, function L is about w, b, ξ
i, ξ
i *Minimal point, also be α
i, α
i *, γ
i, γ
i *Maximal point, minimization problem are converted into the maximization problems of asking its dual problem.
Lagrangian L is about w at the saddle point place, b, ξ
i, ξ
i *Minimal point:
Can get the dual function of Lagrangian:
At this moment,
According to Ku En-Plutarch (KKT) conditional theorem, have following formula to set up at saddle point:
By following formula as seen,
α
iAnd α
i *Portion can not be a non-zero simultaneously, can get:
Can obtain b from following formula, obtain model.Empty the prediction error database, in order to upgrade the usefulness of image data next time.
(3) determine the weight of new model and existing model.Gather the different operating conditions of up-to-date boiler data down as test samples, use the master mould prediction and predict weighted-average method, check data is predicted, i.e. y with new model
i=α D
n+ β D
o, y wherein
iBe i group Optimization of operating target prediction value, D
nBe new model predicted value, D
oBe original model predication value, α is a new model predicted value weight coefficient, β is the prediction weight coefficient of original model, and alpha+beta=1, definite employing particle cluster algorithm iteration optimizing of α and β is determined, definition particle cluster algorithm initial position Z vector respectively tie up component, be respectively new model weight and original model weight beta, object function is: min J=∑ | y
i-Y
i|, Y wherein
iIt is the actual motion value of i group operating mode.When J had reached minimum, reaches setting value or finished iterations, training was finished, and obtains the weight coefficient of new model and existing model.
(4) original model is combined with the weight proportion of new model by optimum, constitute new model, i.e. D=α D
n+ β D
o, wherein D realizes the renewal of boiler combustion optimization model for the model after upgrading.
Claims (1)
1. the model update method of an on-line boiler combustion optimization is characterized in that the step of this method comprises:
Step (1). set up the prediction error database there has been model: according to concrete boiler combustion situation with to the requirement of model prediction precision, the permission predicated error limit δ of setting model, when image data, error between judgment models predicted value and the actual motion value and the size that allows predicated error limit δ, if predicated error is greater than δ, promptly | V
c-V
s|>δ, the data that then will transfinite deposit in the prediction error database, wherein V
cBe model predication value, V
sBe actual operating data;
Step (2). set up new model: when master mould need upgrade, utilize data in the prediction error database as training sample, schedule of samples is shown { x
i, y
i}
I=1 N, x wherein
iRepresent the boiler operating parameter vector of i group, y as the input data
iRepresent the parameter of i group as the sign boiler combustion state of output parameter, adopt algorithm of support vector machine to set up new model, kernel function is elected RBF as
φ (x) is a mapping function, establishes the object function of being asked to be: f (x
i)=w φ (x
i)+b, f (x
i) being the boiler combustion index prediction value of model output, w is the weight coefficient vector, b is an intercept; W, b by
In condition
Find the solution and get, constant c>0 is penalty coefficient; Introduce relaxation factor ξ in the formula
* iAnd ξ
i, ξ
* i〉=0 and ξ
i〉=0, ε is for allowing error of fitting;
This minimization problem is a convex quadratic programming problem, introduces Lagrangian:
At the saddle point place, function L is about w, b, ξ
i, ξ
i *Minimal point, also be α
i, α
i *, γ
i, γ
i *Maximal point, minimization problem are converted into the maximization problems of asking its dual problem;
Lagrangian L is about w at the saddle point place, b, ξ
i, ξ
i *Minimal point:
Can get the dual function of Lagrangian:
At this moment,
According to Ku En-Plutarch (KKT) conditional theorem, have following formula to set up at saddle point:
Can obtain b from following formula, obtain model; Empty the prediction error database, in order to upgrade the usefulness of image data next time;
Step (3). determine the weight of new model and existing model; Data under the different operating conditions of collection boiler are used the master mould prediction and are predicted weighted-average method with new model as test samples, check data are predicted, i.e. y
i=α D
n+ β D
o, y wherein
iBe i group Optimization of operating target prediction value, D
nBe new model predicted value, D
oBe original model predication value, α is a new model predicted value weight coefficient, and β is the prediction weight coefficient of original model, and alpha+beta=1;
Definite employing particle cluster algorithm iteration optimizing of α and β determines, definition particle cluster algorithm initial position Z vector respectively tie up component, be respectively new model weight and original model weight beta, object function is: the minJ=∑ | y
i-Y
i|, Y wherein
iBe the actual motion value of i group operating mode, when J had reached minimum, reaches setting value or finished iterations, training was finished, and obtains the weight coefficient of new model and existing model;
Step (4). original model is combined with the weight proportion of new model by optimum, constitute new model, i.e. D=α D
n+ β D
o, wherein D realizes the renewal of boiler combustion optimization model for the model after upgrading.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN200910096410XA CN101498458B (en) | 2009-03-02 | 2009-03-02 | Model updating method for on-line boiler combustion optimization |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN200910096410XA CN101498458B (en) | 2009-03-02 | 2009-03-02 | Model updating method for on-line boiler combustion optimization |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101498458A CN101498458A (en) | 2009-08-05 |
CN101498458B true CN101498458B (en) | 2010-08-11 |
Family
ID=40945650
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN200910096410XA Expired - Fee Related CN101498458B (en) | 2009-03-02 | 2009-03-02 | Model updating method for on-line boiler combustion optimization |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101498458B (en) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102252343A (en) * | 2011-05-05 | 2011-11-23 | 浙江宜景环保科技有限公司 | Method for optimizing combustion of porous medium combustor |
CN102252342A (en) * | 2011-05-05 | 2011-11-23 | 衢州远景资源再生科技有限公司 | Model updating method for online combustion optimization of porous medium combustor |
CN105808945A (en) * | 2016-03-07 | 2016-07-27 | 杭州电子科技大学 | Mixed intelligent boiler combustion efficiency optimization method |
Families Citing this family (19)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102184287B (en) * | 2011-05-05 | 2012-08-22 | 杭州电子科技大学 | Modelling method for combustion optimization of waste plastics oil refining |
CN102194057B (en) * | 2011-05-05 | 2013-02-06 | 浙江宜景环保科技有限公司 | Method for updating model of waste plastic oil refining combustion optimization |
CN102222128B (en) * | 2011-05-05 | 2013-01-02 | 杭州电子科技大学 | Method for combustion optimization of waste plastics oil refining |
CN102184450A (en) * | 2011-05-05 | 2011-09-14 | 衢州远景资源再生科技有限公司 | Modeling method for combustion optimization of porous medium combustor |
CN102778538B (en) * | 2012-07-06 | 2014-12-10 | 广东电网公司电力科学研究院 | Soft measuring method based on improved SVM (Support Vector Machine) for measuring boiler unburned carbon content in fly ash |
CN102799939A (en) * | 2012-07-16 | 2012-11-28 | 杭州电子科技大学 | Biomass furnace combustion-optimized model updating method |
CN102981409B (en) * | 2012-12-11 | 2015-05-06 | 富通集团有限公司 | Method for optimizing operation of coaxial cable sheathing machine |
CN103064370B (en) * | 2012-12-17 | 2015-05-06 | 富通集团有限公司 | Cooperating production method of multiple coaxial-cable shield machines |
CN103576553B (en) * | 2013-11-06 | 2016-08-24 | 华北电力大学(保定) | A kind of fractional order self-adjusting control of coal-burning boiler vapor (steam) temperature |
CN104879750B (en) * | 2015-06-15 | 2017-05-10 | 神华集团有限责任公司 | Combustion optimizing device, system and method for circulating fluidized bed boiler |
CN106594794B (en) * | 2016-12-22 | 2019-03-08 | 杭州电子科技大学 | A kind of hybrid intelligent boiler efficiency burning optimization model update method |
CN107016176A (en) * | 2017-03-24 | 2017-08-04 | 杭州电子科技大学 | A kind of hybrid intelligent overall boiler burning optimization method |
CN108257121B (en) * | 2018-01-09 | 2019-01-25 | 北京百度网讯科技有限公司 | Method, apparatus, storage medium and the terminal device that product defects detection model updates |
CN108242054A (en) | 2018-01-09 | 2018-07-03 | 北京百度网讯科技有限公司 | A kind of steel plate defect detection method, device, equipment and server |
CN109165244A (en) * | 2018-09-06 | 2019-01-08 | 北京航天控制仪器研究所 | A kind of stream data processing method of the online multidimensional output based on correlation |
CN111240282B (en) * | 2019-12-31 | 2021-12-24 | 联想(北京)有限公司 | Process optimization method, device, equipment and computer readable storage medium |
CN111459030B (en) * | 2020-04-16 | 2022-03-29 | 南京英璞瑞自动化科技有限公司 | Self-adaptive modeling method for closed-loop combustion optimization of boiler |
CN111783968B (en) * | 2020-06-30 | 2024-05-31 | 山东信通电子股份有限公司 | Power transmission line monitoring method and system based on cloud edge cooperation |
CN114880927A (en) * | 2022-04-29 | 2022-08-09 | 广东大唐国际雷州发电有限责任公司 | Intelligent power plant monitoring method, system, equipment and storage medium |
-
2009
- 2009-03-02 CN CN200910096410XA patent/CN101498458B/en not_active Expired - Fee Related
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102252343A (en) * | 2011-05-05 | 2011-11-23 | 浙江宜景环保科技有限公司 | Method for optimizing combustion of porous medium combustor |
CN102252342A (en) * | 2011-05-05 | 2011-11-23 | 衢州远景资源再生科技有限公司 | Model updating method for online combustion optimization of porous medium combustor |
CN102252343B (en) * | 2011-05-05 | 2013-05-08 | 浙江宜景环保科技有限公司 | Method for optimizing combustion of porous medium combustor |
CN105808945A (en) * | 2016-03-07 | 2016-07-27 | 杭州电子科技大学 | Mixed intelligent boiler combustion efficiency optimization method |
CN105808945B (en) * | 2016-03-07 | 2018-11-13 | 杭州电子科技大学 | A kind of hybrid intelligent boiler efficiency burning optimization method |
Also Published As
Publication number | Publication date |
---|---|
CN101498458A (en) | 2009-08-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101498458B (en) | Model updating method for on-line boiler combustion optimization | |
CN101498457B (en) | Boiler combustion optimizing method | |
CN101498459B (en) | Modeling method for boiler combustion optimization | |
CN113657031B (en) | Heat supply dispatching automation realization method, system and platform based on digital twin | |
CN102252342A (en) | Model updating method for online combustion optimization of porous medium combustor | |
CN102252343B (en) | Method for optimizing combustion of porous medium combustor | |
CN102184287B (en) | Modelling method for combustion optimization of waste plastics oil refining | |
CN106773682B (en) | The glass furnace bottom of pond temperature intelligent forecast Control Algorithm being dynamically determined based on time lag | |
CN106594794B (en) | A kind of hybrid intelligent boiler efficiency burning optimization model update method | |
CN102194057B (en) | Method for updating model of waste plastic oil refining combustion optimization | |
CN103869783B (en) | A kind of concentrate yield on-line prediction method | |
CN102842066A (en) | Modeling method for combustion optimization of biomass furnace | |
CN109798646A (en) | A kind of air quantity variable air conditioner control system and method based on big data platform | |
CN112700066A (en) | Optimal time scale coordination method for scheduling of electric-thermal integrated energy system | |
CN103617470A (en) | Optimization method for mineral separation comprehensive production index under equipment capability changing condition | |
CN102184450A (en) | Modeling method for combustion optimization of porous medium combustor | |
CN113589693A (en) | Cement industry decomposing furnace temperature model prediction control method based on neighborhood optimization | |
CN105808945B (en) | A kind of hybrid intelligent boiler efficiency burning optimization method | |
CN102021254B (en) | Distributed coordination control method of blast furnace system and pressure recovery turbine (TRT) device | |
Zhang et al. | An integrated energy efficiency evaluation method for forging workshop based on IoT and data-driven | |
CN110118382A (en) | General operation regulation strategy identification and evaluation method for heat exchange station | |
CN116440670B (en) | Limestone slurry density stability control method | |
CN102799939A (en) | Biomass furnace combustion-optimized model updating method | |
CN107119159A (en) | A kind of optimization method for the blast furnace material distribution process burden distribution matrix that there is integer programming problem | |
CN104656436B (en) | A kind of decomposition furnace outlet Temperature Modeling method |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
C17 | Cessation of patent right | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20100811 Termination date: 20130302 |