CN107563078B - Flame center height coefficient formula correction method and device - Google Patents
Flame center height coefficient formula correction method and device Download PDFInfo
- Publication number
- CN107563078B CN107563078B CN201710807610.6A CN201710807610A CN107563078B CN 107563078 B CN107563078 B CN 107563078B CN 201710807610 A CN201710807610 A CN 201710807610A CN 107563078 B CN107563078 B CN 107563078B
- Authority
- CN
- China
- Prior art keywords
- center height
- flame center
- height coefficient
- coefficient
- formula
- 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
Images
Landscapes
- Regulation And Control Of Combustion (AREA)
Abstract
The embodiment of the invention discloses a flame center height coefficient formula correction method, which can obtain a more accurate hearth outlet flue gas temperature and a flame center height coefficient than the traditional thermodynamic calculation through numerical simulation, and can obtain a flame center height coefficient calculation formula which is similar to the traditional empirical formula form but more fit with the traditional oxygen-enriched combustion boiler by taking an empirical formula of the flame center height coefficient in the traditional thermodynamic calculation as a template and carrying out numerical fitting on a first flame center height coefficient obtained through the numerical simulation, thereby having important guiding significance for the thermodynamic calculation of the boiler under the condition of oxygen-enriched combustion and solving the technical problem that the accurate central height coefficient of the traditional oxygen-enriched combustion boiler cannot be obtained through the empirical formula of the flame center height coefficient of the traditional boiler.
Description
Technical Field
The invention relates to the field of boilers, in particular to a flame center height coefficient formula correction method and device.
Background
The traditional calculation of the height coefficient of the flame center of the boiler is carried out by an empirical formula, but the empirical formula is derived from the introduction of foreign technologies in the last century, and the empirical formula is summarized from the data of the boiler engineering at that time. Due to the development of the current boiler technology, various design parameters are greatly different from those of the current boiler, and the traditional empirical formula is not suitable for calculating the flame center height coefficient of the current oxygen-enriched combustion boiler. Therefore, the technical problem that the accurate central height coefficient of the existing oxyfuel combustion boiler cannot be obtained through the empirical formula of the traditional boiler flame central height coefficient is caused.
Disclosure of Invention
The invention provides a flame center height coefficient formula correction method and a flame center height coefficient formula correction device, and solves the technical problem that the accurate center height coefficient of the conventional oxygen-enriched combustion boiler cannot be obtained through the traditional empirical formula of the flame center height coefficient of the boiler.
The invention provides a flame center height coefficient formula correction method, which comprises the following steps:
s1: calculating the smoke temperature of the outlet of the first furnace corresponding to at least two different first over-fire air distribution rates through numerical simulation;
s2: respectively calculating the flue gas temperature of the outlet of each first furnace to obtain a corresponding first flame center height coefficient;
s3: and (4) performing numerical fitting on the first flame center height coefficient obtained in the step (S2) by taking an empirical calculation formula of the flame center height coefficient in thermal calculation as a template to obtain a corrected flame center height coefficient calculation formula.
Preferably, the step S3 specifically includes:
s301: acquiring an empirical calculation formula of a flame center height coefficient in thermal calculation, and removing the empirical coefficient in the empirical calculation formula to obtain a formula framework;
s302: and (4) performing numerical fitting on the first flame center height coefficient obtained in the step (S2) by taking the formula skeleton as a template to obtain a corrected flame center height coefficient calculation formula.
Preferably, the step S1 is preceded by: step S0;
s0: at least two different first overfire air distribution rates are set.
Preferably, the empirical calculation formula in step S3 is specifically:
wherein M is the flame center height coefficient, M0Is an initial value of the flame center height coefficient, xrFor relative elevation of burner arrangement in boiler, rvIs the inert component ratio of the flue gas in the boiler.
Preferably, the method further comprises the following steps: step S4, step S5, step S6, and step S7;
s4: calculating the temperature of the flue gas at the outlet of the second hearth corresponding to at least one second burn-out air distribution rate different from the first burn-out air distribution rate value through numerical simulation;
s5: calculating the temperature of the flue gas at the outlet of at least one second hearth to obtain a corresponding second flame center height coefficient;
s6: performing error calculation on the second flame center height coefficient and the corrected flame center height coefficient calculation formula through a least square method, and executing the step S7 when the error value is greater than a preset limit value;
s7: and increasing the number of the first overfire air distribution rate values, and returning to the step S1.
The invention provides a flame center height coefficient formula correcting device, which comprises:
the smoke temperature calculation module is used for calculating smoke temperatures of the outlets of the first hearths corresponding to at least two different first over-fire air distribution rates through numerical simulation;
the coefficient calculation module is used for calculating the flue gas temperature of the outlet of each first hearth respectively to obtain the corresponding first flame center height coefficient;
and the fitting correction module is used for performing numerical fitting on the first flame center height coefficient obtained in the coefficient calculation module by taking an empirical calculation formula of the flame center height coefficient in thermal calculation as a template to obtain a corrected flame center height coefficient calculation formula.
Preferably, the fitting correction module specifically includes:
the extraction submodule is used for obtaining an empirical calculation formula of the flame center height coefficient in the thermal calculation, and eliminating the empirical coefficient in the empirical calculation formula to obtain a formula framework;
and the fitting submodule is used for performing numerical fitting on the first flame center height coefficient obtained in the coefficient calculation module by taking the formula skeleton as a template to obtain a corrected flame center height coefficient calculation formula.
Preferably, the method further comprises the following steps: presetting an air distribution module;
and the preset air distribution module is used for setting at least two different first over-fire air distribution rates.
Preferably, the empirical calculation formula in the fitting correction module is specifically:
wherein M is the flame center height coefficient, M0Is an initial value of the flame center height coefficient, xrFor relative elevation of burner arrangement in boiler, rvIs the inert component ratio of the flue gas in the boiler.
Preferably, the method further comprises the following steps: the system comprises a first smoke temperature module, a first coefficient module, a formula verification module and an increase return module;
the second smoke temperature module is used for calculating the smoke temperature of a second hearth outlet corresponding to at least one second burn-out air distribution rate different from the first burn-out air distribution rate value through numerical simulation;
the second coefficient module is used for calculating the temperature of the flue gas at the outlet of at least one second hearth to obtain a corresponding second flame center height coefficient;
the formula verification module is used for carrying out error calculation on the second flame center height coefficient and the corrected flame center height coefficient calculation formula through a least square method, and when the error value is larger than a preset limit value, the addition returning module is triggered;
and the adding and returning module is used for increasing the quantity of the first over-fire air distribution rate values and triggering the smoke temperature calculating module.
According to the technical scheme, the invention has the following advantages:
the invention provides a flame center height coefficient formula correction method, which comprises the following steps: s1: calculating the smoke temperature of the outlet of the first furnace corresponding to at least two different first over-fire air distribution rates through numerical simulation; s2: respectively calculating the flue gas temperature of the outlet of each first furnace to obtain a corresponding first flame center height coefficient; s3: and (4) performing numerical fitting on the first flame center height coefficient obtained in the step (S2) by taking an empirical calculation formula of the flame center height coefficient in thermal calculation as a template to obtain a corrected flame center height coefficient calculation formula.
According to the invention, the more accurate furnace outlet flue gas temperature and flame center height coefficient can be obtained through numerical simulation compared with the traditional thermodynamic calculation, an empirical formula of the flame center height coefficient in the traditional thermodynamic calculation is used as a template, numerical fitting is carried out on the flame center height coefficient obtained through numerical simulation, the approximate form of the traditional empirical formula can be obtained, but the empirical formula is more fit with the flame center height coefficient calculation formula of the existing oxygen-enriched combustion boiler, the method has an important guiding significance for the thermodynamic calculation of the boiler under the oxygen-enriched combustion condition, and the technical problem that the accurate center height coefficient of the existing oxygen-enriched combustion boiler cannot be obtained through the empirical formula of the flame center height coefficient of the traditional boiler is solved.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings used in the description of the embodiments or the prior art will be briefly described below, and it is obvious that the drawings in the following description are only some embodiments of the present invention, and for those skilled in the art, other drawings can be obtained according to these drawings without inventive exercise.
FIG. 1 is a schematic flow chart illustrating an embodiment of a method for correcting a flame center height coefficient formula according to an embodiment of the present invention;
FIG. 2 is a schematic flow chart of another embodiment of a method for correcting a flame center height coefficient formula according to an embodiment of the present invention;
FIG. 3 is a schematic connection diagram of an embodiment of a flame center height coefficient formula modification apparatus according to an embodiment of the present invention;
fig. 4 is a schematic connection diagram of another embodiment of a flame center height coefficient formula modification apparatus according to an embodiment of the present invention.
Detailed Description
The embodiment of the invention provides a flame center height coefficient formula correction method and device, and solves the technical problem that the accurate center height coefficient of the conventional oxygen-enriched combustion boiler cannot be obtained through the traditional empirical formula of the flame center height coefficient of the boiler.
In order to make the objects, features and advantages of the present invention more obvious and understandable, the technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings in the embodiments of the present invention, and it is obvious that the embodiments described below are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
Referring to fig. 1, the present invention provides an embodiment of a method for correcting a flame center height coefficient equation, comprising:
step 101: calculating the smoke temperature of the outlet of the first furnace corresponding to at least two different first over-fire air distribution rates through numerical simulation;
it should be noted that, when the fuel and the intrinsic parameters of the furnace are fixed, the adiabatic combustion temperature of the fuel, the fuel consumption, the furnace wall area, the effective coefficient of the water wall, the heat retention coefficient and the boolean coefficient are all kept unchanged, and the change of the air distribution rate of the over-fire air causes the change of the combustion center, so that the height coefficient of the flame center is changed, and the temperature of the flue gas at the outlet of the furnace is affected.
Step 102: respectively calculating the flue gas temperature of the outlet of each first furnace to obtain a corresponding first flame center height coefficient;
step 103: and (3) performing numerical fitting on the first flame center height coefficient obtained in the step (102) by taking an empirical calculation formula of the flame center height coefficient in thermal calculation as a template to obtain a corrected flame center height coefficient calculation formula.
It should be noted that the empirical calculation formula in step 103 is specifically:
wherein M is the flame center height coefficient, M0Is an initial value of the flame center height coefficient, xrFor relative elevation of burner arrangement in boiler, rvThe inert component ratio of the flue gas in the boiler;
obtaining the flue gas temperature at the outlet of a hearth under different over-fire air rates through numerical simulation, further calculating the flue gas temperature at the outlet of the hearth to obtain a flame center height coefficient M, and then obtaining a calculation expression of the corrected flame center height coefficient M through fitting, wherein the basic form of the calculation expression of the corrected flame center height coefficient M is still the form of an empirical calculation formula, but the empirical coefficient in the empirical calculation formula can be changed or other additional terms can be added;
a boiler thermodynamic calculation method and a boiler numerical simulation method are two most important means in the design and optimization process of a large pulverized coal combustion boiler. Mature thermodynamic calculation methods are derived from a large amount of engineering test data, so that the results of the thermodynamic calculation can be used for checking a numerical simulation model; similarly, the accurate numerical simulation calculation can also be used for correcting the method of thermodynamic calculation, namely the two methods have a coupling relation in predicting the heat transfer in the furnace;
the method has the advantages that the furnace outlet flue gas temperature and the flame center height coefficient which are more accurate than those of the traditional thermodynamic calculation can be obtained through numerical simulation, the empirical formula of the flame center height coefficient in the traditional thermodynamic calculation is used as a template, numerical fitting is carried out on the flame center height coefficient obtained through the numerical simulation, the method can be similar to the traditional empirical formula, the flame center height coefficient calculation formula more fitting the existing oxygen-enriched combustion boiler is obtained, the method has important guiding significance on the thermodynamic calculation of the boiler under the oxygen-enriched combustion condition, and the technical problem that the accurate existing oxygen-enriched combustion boiler center height coefficient cannot be obtained through the empirical formula of the traditional boiler flame center height coefficient is solved.
The above is an embodiment of a method for correcting a flame center height coefficient formula provided by the present invention, and the following is another embodiment of a method for correcting a flame center height coefficient formula provided by the present invention.
Referring to fig. 2, another embodiment of the present invention provides a method for correcting a flame center height coefficient equation, including:
step 201: setting at least two different first over-fire air distribution rates;
it should be noted that the flame center height coefficient obtained by a single overfire air distribution ratio cannot be numerically fitted.
Step 202: calculating the smoke temperature of the outlet of the first furnace corresponding to at least two different first over-fire air distribution rates through numerical simulation;
step 203: respectively calculating the flue gas temperature of the outlet of each first furnace to obtain a corresponding first flame center height coefficient;
it should be noted that the numerical calculation process is as follows:
A. determining a calculation relation of the average temperature of flame in the oxygen-enriched combustion boiler according to the temperature field distribution in the boiler obtained by numerical simulation research of the oxygen-enriched combustion boiler hearth by adopting a similar theoretical method of hearth heat transfer analysis calculation;
wherein, the similar theoretical method of the hearth heat transfer calculation is as follows:
flame temperature and flue gas temperature vary dramatically over their course, but there are many factors that affect the variation of the furnace temperature field along the furnace height, such as: burner placement, fuel properties, structural properties of the heating surface, etc. However, tests show that, for a furnace chamber with a certain height and a water wall distributed all around, the temperature field in the furnace has certain similarity and can be expressed as the following relation:
Θ4＝e-αX-e-βX(1)
in the formulae (1) and (2): t isllTheoretical combustion temperature, in K; x is the relative flame height,
l is the total height of the flame (burner centre to outlet centre), x is the height of the flame from the burner centre, α are empirical coefficients taking into account the effect of heat transfer, combustion, respectively, on flame temperature.
In the formula (1), let X be 1, obtain the fourth power calculation formula of the dimensionless temperature at the furnace outlet:
integrating the formula (1) from 0 to 1 to obtain the average value of the flame temperature of the hearth to the fourth power:
the location of the peak temperature is determined by equation (5):
the position of the highest temperature point is:
due to the fact that
And XmBoth α and β, elimination of α and β after combination (3), (4) and (6), yield:
when X is presentmWhen the time is not changed, the user can select the time,
and
and has a linear relationship. Therefore, the method comprises the following steps:
m and n are both XmI.e.:
plotting this functional relationship will find the intercept of the line to be approximately zero, i.e.
I.e. m ≈ 1. n is actually XmFor different values of the slope of the line, the range of n can be found as follows: n is more than 0.4 and less than or equal to 1.0;
the final flame average temperature was:
namely, it is
In the formula: t ishyIs the average flame temperature in K; t isllTheoretical combustion temperature, in K; t is1"furnace exit smoke temperature, unit K.
B. Based on a basic radiation heat transfer formula, deducing and sorting out a calculation criterion correlation formula of the smoke temperature of the hearth outlet of the oxygen-enriched combustion boiler according to the obtained average temperature of flame in the boiler and by combining numerical simulation data of the hearth of the oxygen-enriched combustion boiler;
C. based on the provided correlation formula of the calculation criterion of the smoke temperature at the outlet of the hearth of the oxygen-enriched combustion boiler, the center height of the flame of the hearth is determined by combining the numerical simulation of the hearth of the oxygen-enriched combustion boiler, and the center height coefficient of the flame is calculated and corrected.
Step 204: acquiring an empirical calculation formula of a flame center height coefficient in thermal calculation, and removing the empirical coefficient in the empirical calculation formula to obtain a formula framework;
it should be noted that the formula skeleton refers to a formula form in which a coefficient obtained according to engineering experience in an empirical calculation formula is removed and changed into an unknown undetermined term, so that the undetermined term exists.
Step 205: performing numerical fitting on the first flame center height coefficient obtained in the step 203 by taking a formula framework as a template to obtain a corrected flame center height coefficient calculation formula;
step 206: calculating the temperature of the flue gas at the outlet of the second hearth corresponding to at least one second burn-out air distribution rate different from the first burn-out air distribution rate value through numerical simulation;
step 207: calculating the temperature of the flue gas at the outlet of at least one second hearth to obtain a corresponding second flame center height coefficient;
it should be noted that, steps 206 and 207 are intended to perform error calculation on the corrected flame center height coefficient according to obtaining a second flame center height coefficient different from the first flame center height coefficient, so as to obtain more accurate effect, and if the first flame center height coefficient is used for error analysis, the error analysis result obtained may have deviation because the corrected flame center height coefficient is numerically fitted by the first flame center height coefficient.
Step 208: calculating the error of the second flame center height coefficient and the corrected flame center height coefficient by a least square method, and executing step 209 when the error is greater than a preset limit value;
it should be noted that, the least square method is an example of the error calculation method in the present invention, and there are various implementation manners for error calculation, which are not listed in the present invention.
Step 209: and increasing the number of the first overfire air distribution rate values, and returning to the step 202.
It should be noted that, if the error value is greater than the preset limit value, the value of the first over-fire air distribution rate is increased to obtain more first flame center height coefficients, and a more accurate correction formula can be obtained by performing numerical fitting through more data;
the invention analyzes the difference of the heat transfer characteristics of the pulverized coal boiler and the traditional pulverized coal boiler by researching the flowing, burning and heat transfer characteristics of the pulverized coal boiler under the condition of oxygen-enriched combustion, corrects and improves the pulverized coal boiler on the basis of the traditional pulverized coal boiler thermodynamic calculation method and numerical simulation method, researches a thermodynamic calculation method suitable for the oxygen-enriched combustion boiler and a relevant submodel for the numerical simulation of the oxygen-enriched combustion boiler, and provides a coupling method of the thermodynamic calculation and the numerical simulation of the oxygen-enriched combustion boiler;
the method has the advantages that the furnace outlet flue gas temperature and the flame center height coefficient which are more accurate than those of the traditional thermodynamic calculation can be obtained through numerical simulation, the empirical formula of the flame center height coefficient in the traditional thermodynamic calculation is used as a template, numerical fitting is carried out on the flame center height coefficient obtained through the numerical simulation, the method can be similar to the traditional empirical formula, the flame center height coefficient calculation formula more fitting the existing oxygen-enriched combustion boiler is obtained, the method has important guiding significance on the thermodynamic calculation of the boiler under the oxygen-enriched combustion condition, and the technical problem that the accurate existing oxygen-enriched combustion boiler center height coefficient cannot be obtained through the empirical formula of the traditional boiler flame center height coefficient is solved.
The above is another embodiment of the method for correcting the flame center height coefficient formula provided by the present invention, and the following is an embodiment of the device for correcting the flame center height coefficient formula provided by the present invention.
Referring to fig. 3, the present invention provides an embodiment of a flame center height coefficient formula correction apparatus, including:
the smoke temperature calculation module 301 is configured to calculate, through numerical simulation, smoke temperatures of first hearth outlets corresponding to at least two different first over-fire air distribution rates;
a coefficient calculation module 302, configured to calculate a flue gas temperature at an outlet of each first furnace, respectively, to obtain a corresponding first flame center height coefficient;
and the fitting correction module 303 is configured to perform numerical fitting on the first flame center height coefficient obtained in the coefficient calculation module by using an empirical calculation formula of the flame center height coefficient in thermal calculation as a template, so as to obtain a corrected flame center height coefficient calculation formula.
The above is an embodiment of a flame center height coefficient formula correction device provided by the present invention, and the following is another embodiment of a flame center height coefficient formula correction device provided by the present invention.
Referring to fig. 4, another embodiment of the present invention provides a flame center height coefficient formula correction apparatus, including:
the smoke temperature calculation module 402 is configured to calculate, through numerical simulation, smoke temperatures of first hearth outlets corresponding to at least two different first over-fire air distribution rates;
the coefficient calculation module 403 is configured to calculate a flue gas temperature at an outlet of each first furnace, and obtain a corresponding first flame center height coefficient;
a fitting correction module 404, configured to perform numerical fitting on the first flame center height coefficient obtained in the coefficient calculation module 403 by using an empirical calculation formula of the flame center height coefficient in thermal calculation as a template, so as to obtain a corrected calculation formula of the flame center height coefficient.
Further, the fitting correction module 404 specifically includes:
the extraction submodule 4041 is used for obtaining an empirical calculation formula of the flame center height coefficient in the thermodynamic calculation, and eliminating the empirical coefficient in the empirical calculation formula to obtain a formula skeleton;
the fitting submodule 4042 is configured to perform numerical fitting on the first flame center height coefficient obtained in the coefficient calculation module 403 with a formula skeleton as a template to obtain a corrected flame center height coefficient calculation formula.
Further, still include: a preset air distribution module 401;
the preset air distribution module 401 is configured to set at least two different first over-fire air distribution rates.
Further, the empirical calculation formula in the fitting correction module 404 is specifically:
wherein M is the flame center height coefficient, M0Is an initial value of the flame center height coefficient, xrFor relative elevation of burner arrangement in boiler, rvIs the inert component ratio of the flue gas in the boiler.
Further, still include: a second smoke temperature module 405, a second coefficient module 406, a formula verification module 407, and an increase return module 408;
a second flue gas temperature module 405, configured to calculate, through numerical simulation, a second furnace outlet flue gas temperature corresponding to at least one second burn-out air distribution rate that is different from the first burn-out air distribution rate value;
a second coefficient module 406, configured to calculate a temperature of flue gas at an outlet of the at least one second furnace, to obtain a corresponding second flame center height coefficient;
the formula verification module 407 is configured to perform error calculation on the second flame center height coefficient and the corrected flame center height coefficient calculation formula by using a least square method, and trigger the increase and return module 408 when the error is greater than a preset limit value;
and a returning module 408 is added for increasing the number of the first over-fire air distribution rate values and triggering the smoke temperature calculation module 402.
It is clear to those skilled in the art that, for convenience and brevity of description, the specific working processes of the above-described systems, apparatuses and units may refer to the corresponding processes in the foregoing method embodiments, and are not described herein again.
In the several embodiments provided in the present application, it should be understood that the disclosed apparatus and method may be implemented in other ways. For example, the above-described apparatus embodiments are merely illustrative, and for example, the division of the modules is merely a logical division, and in actual implementation, there may be other divisions, for example, multiple modules or components may be combined or integrated into another system, or some features may be omitted, or not implemented. In addition, the shown or discussed mutual coupling or direct coupling or communication connection may be an indirect coupling or communication connection through some interfaces, devices or modules, and may be in an electrical, mechanical or other form.
The modules described as separate parts may or may not be physically separate, and parts displayed as modules may or may not be physical modules, may be located in one place, or may be distributed on a plurality of network modules. Some or all of the modules may be selected according to actual needs to achieve the purpose of the solution of the present embodiment.
In addition, functional modules in the embodiments of the present invention may be integrated into one processing module, or each of the modules may exist alone physically, or two or more modules are integrated into one module. The integrated module can be realized in a hardware mode, and can also be realized in a software functional module mode.
The integrated module, if implemented in the form of a software functional module and sold or used as a stand-alone product, may be stored in a computer readable storage medium. Based on such understanding, the technical solution of the present invention may be embodied in the form of a software product, which is stored in a storage medium and includes instructions for causing a computer device (which may be a personal computer, a server, or a network device) to execute all or part of the steps of the method according to the embodiments of the present invention. And the aforementioned storage medium includes: a U-disk, a removable hard disk, a Read-only Memory (ROM), a Random Access Memory (RAM), a magnetic disk or an optical disk, and other various media capable of storing program codes.
The above-mentioned embodiments are only used for illustrating the technical solutions of the present invention, and not for limiting the same; although the present invention has been described in detail with reference to the foregoing embodiments, it will be understood by those of ordinary skill in the art that: the technical solutions described in the foregoing embodiments may still be modified, or some technical features may be equivalently replaced; and such modifications or substitutions do not depart from the spirit and scope of the corresponding technical solutions of the embodiments of the present invention.
Claims (8)
1. A flame center height coefficient formula correction method is characterized by comprising the following steps:
s1: calculating the smoke temperature of the outlet of the first furnace corresponding to at least two different first over-fire air distribution rates through numerical simulation;
s2: respectively calculating the flue gas temperature of the outlet of each first furnace to obtain a corresponding first flame center height coefficient;
s3: performing numerical fitting on the first flame center height coefficient obtained in the step S2 by using an empirical calculation formula of the flame center height coefficient in thermal calculation as a template to obtain a corrected flame center height coefficient calculation formula;
the step S3 specifically includes:
s301: acquiring an empirical calculation formula of a flame center height coefficient in thermal calculation, and removing the empirical coefficient in the empirical calculation formula to obtain a formula framework;
s302: performing numerical fitting on the first flame center height coefficient obtained in the step S2 by taking a formula skeleton as a template to obtain a corrected flame center height coefficient calculation formula;
step S301 specifically includes: the method comprises the steps of obtaining an empirical calculation formula of a flame center height coefficient in thermal calculation, removing the coefficient obtained according to engineering experience in the empirical calculation formula, and changing the removed coefficient into an unknown undetermined term to obtain a formula framework.
2. The method for correcting the flame center height coefficient formula as claimed in claim 1, wherein the step S1 is preceded by the steps of: step S0;
s0: at least two different first overfire air distribution rates are set.
3. The method of claim 1, wherein the empirical formula in step S3 is specifically:
wherein M is the first flame center height coefficient, M0Is an initial value of the flame center height coefficient, xrFor relative elevation of burner arrangement in boiler, rvIs the inert component ratio of the flue gas in the boiler.
4. The method of claim 1, further comprising: step S4, step S5, step S6, and step S7;
s4: calculating the temperature of the flue gas at the outlet of the second hearth corresponding to at least one second burn-out air distribution rate different from the first burn-out air distribution rate value through numerical simulation;
s5: calculating the temperature of the flue gas at the outlet of at least one second hearth to obtain a corresponding second flame center height coefficient;
s6: performing error calculation on the second flame center height coefficient and the corrected flame center height coefficient calculation formula through a least square method, and executing the step S7 when the error value is greater than a preset limit value;
s7: and increasing the number of the first overfire air distribution rate values, and returning to the step S1.
5. A flame center height coefficient formula correction apparatus, comprising:
the smoke temperature calculation module is used for calculating smoke temperatures of the outlets of the first hearths corresponding to at least two different first over-fire air distribution rates through numerical simulation;
the coefficient calculation module is used for calculating the flue gas temperature of the outlet of each first hearth respectively to obtain the corresponding first flame center height coefficient;
the fitting correction module is used for performing numerical fitting on the first flame center height coefficient obtained in the coefficient calculation module by taking an empirical calculation formula of the flame center height coefficient in thermal calculation as a template to obtain a corrected flame center height coefficient calculation formula;
the fitting correction module specifically comprises:
the extraction submodule is used for obtaining an empirical calculation formula of the flame center height coefficient in the thermal calculation, and eliminating the empirical coefficient in the empirical calculation formula to obtain a formula framework;
the fitting submodule is used for performing numerical fitting on the first flame center height coefficient obtained in the coefficient calculation module by taking a formula framework as a template to obtain a corrected flame center height coefficient calculation formula;
the extraction submodule is specifically used for obtaining an empirical calculation formula of the flame center height coefficient in thermodynamic calculation, eliminating the coefficient obtained according to engineering experience in the empirical calculation formula, and changing the eliminated coefficient into an unknown undetermined term to obtain a formula framework.
6. The flame center height coefficient formula correction apparatus of claim 5, further comprising: presetting an air distribution module;
and the preset air distribution module is used for setting at least two different first over-fire air distribution rates.
7. The flame center height coefficient formula correction device of claim 5, wherein the empirical calculation formula in the fitting correction module is specifically:
wherein M is the first flame center height coefficient, M0Is an initial value of the flame center height coefficient, xrFor relative elevation of burner arrangement in boiler, rvIs the inert component ratio of the flue gas in the boiler.
8. The flame center height coefficient formula correction apparatus of claim 5, further comprising: the system comprises a first smoke temperature module, a first coefficient module, a formula verification module and an increase return module;
the second smoke temperature module is used for calculating the smoke temperature of a second hearth outlet corresponding to at least one second burn-out air distribution rate different from the first burn-out air distribution rate value through numerical simulation;
the second coefficient module is used for calculating the temperature of the flue gas at the outlet of at least one second hearth to obtain a corresponding second flame center height coefficient;
the formula verification module is used for carrying out error calculation on the second flame center height coefficient and the corrected flame center height coefficient calculation formula through a least square method, and when the error value is larger than a preset limit value, the addition returning module is triggered;
and the adding and returning module is used for increasing the quantity of the first over-fire air distribution rate values and triggering the smoke temperature calculating module.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710807610.6A CN107563078B (en) | 2017-09-08 | 2017-09-08 | Flame center height coefficient formula correction method and device |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710807610.6A CN107563078B (en) | 2017-09-08 | 2017-09-08 | Flame center height coefficient formula correction method and device |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107563078A CN107563078A (en) | 2018-01-09 |
| CN107563078B true CN107563078B (en) | 2020-09-01 |
Family
ID=60980305
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710807610.6A Active CN107563078B (en) | 2017-09-08 | 2017-09-08 | Flame center height coefficient formula correction method and device |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107563078B (en) |
Families Citing this family (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108870384B (en) * | 2018-04-28 | 2019-07-30 | 东南大学 | The burning of low nitrogen burning circulating fluidized bed boiler and SNCR denitration cooperative optimization method |
| CN111292810A (en) * | 2018-12-06 | 2020-06-16 | 日立汽车系统（中国）有限公司 | Method and device for constructing skeleton mechanism of combustion chemical reaction |
| CN111998383B (en) * | 2020-07-28 | 2024-06-07 | 中国大唐集团科学技术研究院有限公司火力发电技术研究院 | Over-fire air control method based on boiler load and flame central point coordinate quantification |
| CN111947181B (en) * | 2020-08-07 | 2022-04-12 | 中国大唐集团科学技术研究院有限公司华东电力试验研究院 | Method for monitoring central position of flame of power station boiler |
| CN114840898B (en) * | 2022-05-13 | 2024-06-07 | 合肥工业大学智能制造技术研究院 | Method for calculating optimal distance between cooking bench and pan bottom and cooking bench |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102425807A (en) * | 2011-11-23 | 2012-04-25 | 华北电力大学（保定） | Feedforward-feedback composite optimization control method for pulverized coal fired boiler combustion |
| US20120152158A1 (en) * | 2009-12-17 | 2012-06-21 | Mitsubishi Heavy Industries, Ltd. | Solid-fuel-fired burner and solid-fuel-fired boiler |
| CN103148506A (en) * | 2013-03-29 | 2013-06-12 | 广东电网公司电力科学研究院 | Secondary air distribution method for pulverized coal boiler with swirling combustion of front-back hedging |
-
2017
- 2017-09-08 CN CN201710807610.6A patent/CN107563078B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20120152158A1 (en) * | 2009-12-17 | 2012-06-21 | Mitsubishi Heavy Industries, Ltd. | Solid-fuel-fired burner and solid-fuel-fired boiler |
| CN102425807A (en) * | 2011-11-23 | 2012-04-25 | 华北电力大学（保定） | Feedforward-feedback composite optimization control method for pulverized coal fired boiler combustion |
| CN103148506A (en) * | 2013-03-29 | 2013-06-12 | 广东电网公司电力科学研究院 | Secondary air distribution method for pulverized coal boiler with swirling combustion of front-back hedging |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107563078A (en) | 2018-01-09 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107563078B (en) | Flame center height coefficient formula correction method and device | |
| Díez et al. | Modelling of pulverized coal boilers: review and validation of on-line simulation techniques | |
| CN109959283B (en) | thermal checking method and system for high-temperature cooler | |
| CN103267684B (en) | A kind of station boiler pressure restraining element life consumption acquisition methods and system | |
| CN107292006B (en) | Start calculation method and system of supercritical boiler | |
| CN110941186B (en) | Steam temperature control optimization method based on neural network and universal gravitation search algorithm | |
| KR102216820B1 (en) | Test planning device and test planning method | |
| CN103243190B (en) | Method for predicting coal gas consumption of hot blast stove | |
| CN107391894A (en) | A kind of labyrinth calculation of radiation heat transferring method | |
| CN115374636B (en) | Wall temperature prediction method for coupling numerical simulation and performance calculation of boiler | |
| CN103778334A (en) | Method for measuring boiler thermal efficiency of coal fired power plant in real time | |
| CN110442974B (en) | Method and device for optimizing performance of regenerative chamber of horseshoe flame glass kiln | |
| Landfahrer et al. | Numerical model incorporating different oxidizer in a reheating furnace fired with natural gas | |
| Drosatos et al. | Decoupled CFD simulation of furnace and heat exchangers in a lignite utility boiler | |
| Dong et al. | Study on the effect of flame offset on water wall tube temperature in 600° C and 700° C ultra-supercritical boiler | |
| CN107091700A (en) | Temperature Distribution flexible measurement method in burner hearth based on burner hearth Multi sectional slagging situation | |
| CN116359275B (en) | Method, system, equipment and medium for determining hydrodynamic phase transition point of boiler | |
| CN107451304B (en) | Mechanism modeling calculation method of reheater transfer function model | |
| CN110747304B (en) | Method and equipment for determining gas consumption of front-end combustion furnace in blast furnace hot blast stove | |
| CN113359895B (en) | Wall pipe wall overtemperature protection system and control method thereof | |
| CN109583788A (en) | A kind of accounting method and system of fossil-fired unit CO2 emissions | |
| CN112365564A (en) | Temperature monitoring method and device, electronic equipment and readable storage medium | |
| CN114370648A (en) | Automatic control method and device for boiler | |
| CN107977464B (en) | Method and device for accessing benchmarking value data | |
| CN113266844A (en) | Wall temperature deviation control method and device based on opposed firing boiler |
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 | ||
| TA01 | Transfer of patent application right |
Effective date of registration: 20200110 Address after: 510000 room 501-503, annex building, Huaye building, No.1-3, Chumo new street, Xihua Road, Yuexiu District, Guangzhou City, Guangdong Province Applicant after: Guangdong Electric Power Research Institute of energy technology limited liability company Address before: 510080 Dongfeng East Road, Dongfeng, Guangdong, Guangzhou, Zhejiang Province, No. 8 Applicant before: ELECTRIC POWER RESEARCH INSTITUTE, GUANGDONG POWER GRID CO., LTD. Applicant before: Guangdong Electric Power Research Institute of energy technology limited liability company |
|
| TA01 | Transfer of patent application right | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CP01 | Change in the name or title of a patent holder |
Address after: Room 501-503, annex building, Huaye building, No.1-3 Chuimao new street, Xihua Road, Yuexiu District, Guangzhou City, Guangdong Province 510000 Patentee after: China Southern Power Grid Power Technology Co.,Ltd. Address before: Room 501-503, annex building, Huaye building, No.1-3 Chuimao new street, Xihua Road, Yuexiu District, Guangzhou City, Guangdong Province 510000 Patentee before: GUANGDONG ELECTRIC POWER SCIENCE RESEARCH INSTITUTE ENERGY TECHNOLOGY Co.,Ltd. |
|
| CP01 | Change in the name or title of a patent holder |