CN106838979B - Based on boiler furnace grid detection CO and O2The distribution field construction method of concentration - Google Patents
Based on boiler furnace grid detection CO and O2The distribution field construction method of concentration Download PDFInfo
- Publication number
- CN106838979B CN106838979B CN201610981224.4A CN201610981224A CN106838979B CN 106838979 B CN106838979 B CN 106838979B CN 201610981224 A CN201610981224 A CN 201610981224A CN 106838979 B CN106838979 B CN 106838979B
- Authority
- CN
- China
- Prior art keywords
- boiler
- concentration
- flue gas
- model
- furnace
- 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
- 238000001514 detection method Methods 0.000 title claims abstract description 18
- 238000010276 construction Methods 0.000 title claims abstract description 14
- 238000000034 method Methods 0.000 claims abstract description 55
- UGFAIRIUMAVXCW-UHFFFAOYSA-N Carbon monoxide Chemical compound [O+]#[C-] UGFAIRIUMAVXCW-UHFFFAOYSA-N 0.000 claims abstract description 39
- 239000003546 flue gas Substances 0.000 claims abstract description 39
- 238000002485 combustion reaction Methods 0.000 claims abstract description 36
- 239000002817 coal dust Substances 0.000 claims abstract description 24
- 238000012546 transfer Methods 0.000 claims abstract description 17
- 239000012530 fluid Substances 0.000 claims abstract description 9
- 239000002245 particle Substances 0.000 claims description 32
- 238000005070 sampling Methods 0.000 claims description 19
- 239000003245 coal Substances 0.000 claims description 18
- 229910052760 oxygen Inorganic materials 0.000 claims description 16
- 238000009434 installation Methods 0.000 claims description 12
- 238000006243 chemical reaction Methods 0.000 claims description 11
- 239000007787 solid Substances 0.000 claims description 11
- 238000009792 diffusion process Methods 0.000 claims description 10
- 230000005855 radiation Effects 0.000 claims description 10
- 239000007789 gas Substances 0.000 claims description 7
- 208000011580 syndromic disease Diseases 0.000 claims description 6
- 239000000779 smoke Substances 0.000 claims description 5
- 239000000571 coke Substances 0.000 claims description 4
- 238000012545 processing Methods 0.000 claims description 4
- 230000004913 activation Effects 0.000 claims description 3
- 230000015572 biosynthetic process Effects 0.000 claims description 3
- 238000004422 calculation algorithm Methods 0.000 claims description 3
- 210000000038 chest Anatomy 0.000 claims description 3
- 239000003517 fume Substances 0.000 claims description 3
- 239000002594 sorbent Substances 0.000 claims description 3
- 239000000523 sample Substances 0.000 description 18
- 230000008569 process Effects 0.000 description 8
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 description 7
- 239000001301 oxygen Substances 0.000 description 7
- 238000010586 diagram Methods 0.000 description 6
- 238000005259 measurement Methods 0.000 description 5
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 4
- 239000000843 powder Substances 0.000 description 3
- 238000004088 simulation Methods 0.000 description 3
- 230000008901 benefit Effects 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000001914 filtration Methods 0.000 description 2
- 239000008187 granular material Substances 0.000 description 2
- 238000005457 optimization Methods 0.000 description 2
- 238000001556 precipitation Methods 0.000 description 2
- 230000009467 reduction Effects 0.000 description 2
- 239000007921 spray Substances 0.000 description 2
- 239000000126 substance Substances 0.000 description 2
- 230000001052 transient effect Effects 0.000 description 2
- 238000000342 Monte Carlo simulation Methods 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 238000013528 artificial neural network Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004939 coking Methods 0.000 description 1
- 239000000470 constituent Substances 0.000 description 1
- 238000005260 corrosion Methods 0.000 description 1
- 230000007797 corrosion Effects 0.000 description 1
- 230000002950 deficient Effects 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 230000004069 differentiation Effects 0.000 description 1
- 239000000428 dust Substances 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- 238000004134 energy conservation Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 238000001704 evaporation Methods 0.000 description 1
- 230000008020 evaporation Effects 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
- 238000010304 firing Methods 0.000 description 1
- 239000000446 fuel Substances 0.000 description 1
- 238000000227 grinding Methods 0.000 description 1
- 238000010438 heat treatment Methods 0.000 description 1
- 239000013067 intermediate product Substances 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 238000005461 lubrication Methods 0.000 description 1
- 230000000149 penetrating effect Effects 0.000 description 1
- 238000010248 power generation Methods 0.000 description 1
- 239000000047 product Substances 0.000 description 1
- 238000000197 pyrolysis Methods 0.000 description 1
- 230000001869 rapid Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 238000012549 training Methods 0.000 description 1
- 230000005514 two-phase flow Effects 0.000 description 1
- 238000010977 unit operation Methods 0.000 description 1
- 238000009834 vaporization Methods 0.000 description 1
- 230000008016 vaporization Effects 0.000 description 1
Classifications
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F23—COMBUSTION APPARATUS; COMBUSTION PROCESSES
- F23N—REGULATING OR CONTROLLING COMBUSTION
- F23N5/00—Systems for controlling combustion
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F23—COMBUSTION APPARATUS; COMBUSTION PROCESSES
- F23N—REGULATING OR CONTROLLING COMBUSTION
- F23N2223/00—Signal processing; Details thereof
- F23N2223/06—Sampling
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F23—COMBUSTION APPARATUS; COMBUSTION PROCESSES
- F23N—REGULATING OR CONTROLLING COMBUSTION
- F23N2223/00—Signal processing; Details thereof
- F23N2223/10—Correlation
Landscapes
- Engineering & Computer Science (AREA)
- Chemical & Material Sciences (AREA)
- Combustion & Propulsion (AREA)
- Mechanical Engineering (AREA)
- General Engineering & Computer Science (AREA)
- Regulation And Control Of Combustion (AREA)
Abstract
The invention discloses be based on boiler furnace grid detection CO and O2The distribution field construction method of concentration belongs to boiler combustion efficiency detection method technical field, by boiler furnace region wall surface high-temperature flue gas CO, O2The grid type multiple spot on-line checking of concentration and flue gas pressures establishes the flowing, burning of coal dust and heat transfer model in burner hearth using the method for Fluid Mechanics Computation, to predict flue gas CO and O2Motion profile in burner hearth, flue, and then construct flue gas CO and O2Concentration field, provide effectively reliable crucial foundation for boiler combustion adjustment.
Description
Technical field
The invention belongs to boiler combustion efficiency detection method technical fields, are related to CO and O2The distribution field building side of concentration
Method, and in particular to based on boiler furnace grid detection CO and O2The distribution field construction method of concentration.
Background technique
With the development and progress of society, thermal power plant faces the increasing pressure, and electricity needs is continuously increased, to power generation
The requirement of quality is increasingly harsh, and traditional voltage, frequency have not been unique performance assessment criteria, and environmental factor is gradually included,
Form new measurement standard.Coal quality fluctuation and Large-scale machine set increasingly frequently participate in frequency modulation peak regulation, full working scope varying duty fortune
All to unit operation level, more stringent requirements are proposed for capable and long-term middle-low load operation.Currently, power plant is more and more
It faces from power grid, environment and the pressure of itself, energy conservation, consumption reduction, the cry of emission reduction are surging, and burning optimization technology increasingly obtains weight
Depending on.The height of boiler combustion efficiency, into burner hearth fuel can full combustion, the whole efficiency and warp of boiler will be directly affected
Ji property.Therefore realizing that boiler most effectively burns by Combustion System is the effective way for improving boiler whole efficiency.
Generating plant pulverized coal boiler Combustion System design basic ideas are to realize boiler using suitable excess air coefficient as principle
Optimum control of burning is to guarantee suitable excess air coefficient.When the quantity combusted variation of boiler, correspondingly change air output,
To guarantee suitable excess air coefficient, the imperfect combustion loss of boiler and heat loss due to exhaust gas are reduced, makes boiler operatiopn in best
Combustion conditions.Simple measurement O2It can not reflect the quality of coal dust and air troubled water in boiler, even if O2Enough, if it exists
The reasons such as bad are mixed, can also incomplete combustion loss be made to increase.In the actual motion of boiler, only can to economizer exit with
The O of flue gas afterwards2The on-line checking of finite point is carried out, sample point is very limited.For heavy duty boiler, smoke components
It is unevenly distributed and is easy to make oxygen amount detected value not have representativeness.In addition, the originals such as oxygen amount detected value itself is due to leaking out, zero drift
Because that can not reflect the true oxygen amount value in flue gas.For these reasons, in boiler actual moving process, operations staff can only
It is controlled and is burnt with biggish excess air coefficient, it is very uneconomical, cause boiler operatiopn to be easy to happen burner hearth high temperature corrosion, flame
The operations such as deflection, combustion lag, coking and slagging, efficiency of combustion be relatively low puzzlement.A perfect detection means how is found to come pair
Running boiler is adjusted, and is very must with the actual operating efficiency for solving the problems, such as that boiler operatiopn has and improve boiler
It wants.
Current existing related realization scheme is flue gas oxygen content soft-measuring technique: the technology selects steam flow, coal-grinding
Machine coal-supplying amount, total blast volume, total primary wind pressure, pressure fan electric current, air-introduced machine electric current, bellows and burner hearth differential pressure, water supply and smoke evacuation temperature
11 parameters such as degree are established flue gas using feedforward neural network and are contained as input signal using flue gas oxygen content as output parameter
Oxygen amount soft-sensing model.Under various different load operating conditions, multiple groups actual operating data is measured, extracts 1~2 group of data as mould
Type detection uses, remaining is as training sample.Such method solves that oxygen content measurement lag, error are big to ask to a certain extent
Topic.But the accuracy of this hard measurement result is dependent on the integrality of Boundary Test and comprehensive, often due to as-fired coal coal
The fluctuation of matter, the variation of operating condition cause error to increase, and during firing optimization, the application of soft-measuring technique is difficult large-scale
Effectively promote.
Summary of the invention
The present invention in order to overcome the drawbacks of the prior art, is devised based on boiler furnace grid detection CO and O2Concentration
Distribution field construction method establishes the flowing, burning of coal dust and heat transfer model in burner hearth using the method for Fluid Mechanics Computation, with pre-
Survey flue gas CO and O2Motion profile in burner hearth, flue, and then construct flue gas CO and O2Concentration field, for boiler combustion adjust
Effectively reliable crucial foundation is provided.
Specific technical solution adopted by the present invention is: being based on boiler furnace grid detection CO and O2The distribution field structure of concentration
Construction method, it is important to: the method includes:
A, structure, the operation data of station boiler are obtained；
B, to high-temperature flue gas CO, O of boiler furnace region wall surface2Concentration and flue gas pressures carry out grid type multiple spot and examine online
It surveys；
C, the flowing, burning of coal dust and heat transfer model in burner hearth are established using the method for Fluid Mechanics Computation, predicts flue gas
CO and O2Motion profile in burner hearth, flue constructs flue gas CO and O2Concentration field.
The method specifically includes:
A, structure, the operation data of station boiler are obtained: including station boiler dry bottom hopper, burner hearth, furnace arch, furnace nose, burner spray
The detailed dimensions structure of mouth, economizer；Operation data includes at least: the total blast volume of every boiler, total coal amount, each burner
Coal-fired capacity, air quantity, wind-warm syndrome, secondary air temperature, primary air ratio, a secondary air ratio；
B, sampling method is arranged based on burner hearth grid multiple spot, the representative sample of flue gas is taken from burner hearth, i.e., it will probe list
Member is mounted on predetermined position, and the fume sample of hearth region wall surface is obtained by contact unit；
C, the calculating grid of burner hearth and flue is established, wherein the computational domain of full boiler is selected to boiler economizer exit, meter
It calculates domain and is divided into four parts, respectively dry bottom hopper region I, burner region II, burner are up to furnace outlet region
III, furnace outlet to economizer exit region IV；Dry bottom hopper I, burner with up to furnace outlet region III, furnace outlet extremely
Economizer exit region IV uses structured grid, and burner region II uses unstructured grid；
Each region is calculated using three-dimensional steady state, using the non-coupled solution of flowing, heat transfer and combustion model, is used in calculating
Finite Volume Method for Air is discrete to differential governing equation and transport equation to carry out numerical solution for Algebraic Equation set, under each operating condition
CO、O2It is Coal-fired capacity, total blast volume, a wind-warm syndrome, secondary air temperature, primary with the actual measured values of flue gas pressures as input condition
Wind rate, secondary air ratio and air leak rate of air curtain are restrained as boundary condition successive ignition until meeting, and predict flue gas CO, O2In burner hearth to province
Motion profile in the outlet of coal device, constructs CO and O in boiler combustion process2The distributed in three dimensions field of concentration.
Step b is the regulation with American Society of Mechanical Engineers (ASME) standard for smoke sampling: the quantity of measuring point and
Arrangement depends on the size and shape of flue, and sample point interval should be no more than 915mm, and the sum of point used should be no less than at 4 points；
To the boiler of different types, arrange that the sample point number of installation is in strict accordance with ASME in the primary combustion zone addressing of boiler furnace
Standard executes, and points selection is as follows: unit capacity 300MW, 16~24 points, analytical unit 4 of installation points；Unit capacity
600MW, 16~24 points, analytical unit 4 of installation points；Unit capacity 1000MW, 16~24 points of installation points, analytical unit 4
It is a.
The combustion model of coal dust uses diffusion/reaction combined control model, and main includes two, respectively diffusion term and biography
Hot-activity, diffusion term are as follows:
Wherein q is coke consumption rate, and A pre-exponential factor, E is apparent activation energy, POx, sFor the O at particle surface2Partial pressure,
noxTo react order；Conduct heat item are as follows:
Wherein mcFor granular mass, Tg, Tp, TrRespectively gas phase, particle and background radiation temperature, cpFor solid specific heat, fh
For by the reaction heat share of solid sorbent, HreacFor gas-solid reaction heat, εpFor particle blackness, σ is Stefan-Boltzmann constant.
The flow trace model of pulverized coal particle, using stochastic particle tracking model, according to the number of jet stream determine coal dust with
Track granule number；It is assumed that jet stream number is 20, it is meant that and the particle tracked before is divided into 20, the flow that each particle occupies
It is exactly original 1/20, then assigns different velocity fluctuation values to this 20 particles, acquire respective particle trajectories.
It selects when being simulated the heat transfer model to obtain coal dust in the radiant heat transfer to coal dust by based on heat flow method
P1 radiation patterns, compared with discrete ordinates method, Pl radiation patterns take into account radiating scattering.
The PDF model of the method selection standard k- ε two-equation turbulence model, non-premixed combustion Passive Scalar and
The pollution formation model that FLUENT software provides；For the coupled problem of speed and pressure, using the algorithm of SIMPLE；To furnace
The processing of thorax region wall surface uses wall-function method；Simultaneously using the method for owing relaxation.
The beneficial effects of the present invention are: by boiler furnace region wall surface high-temperature flue gas CO, O2Concentration and flue gas pressures
Grid type multiple spot on-line checking, using the method for Fluid Mechanics Computation establish coal dust in burner hearth flowing, burning and heat transfer mould
Type, to predict flue gas CO and O2Motion profile in burner hearth, flue, and then construct flue gas CO and O2Concentration field, for boiler fire
It burns adjustment and effectively reliable crucial foundation is provided.
Using the method based on Computational Fluid Dynamics, it is excellent to can use existing many model library progress for high efficiency and time conservation
Change processing；Intuitive reaction flue gas CO, O2It is mixed can clearly to grasp wind powder in burner hearth for motion profile and the regularity of distribution in burner hearth
Conjunction situation is conducive to the wind powder proportion for targetedly individually adjusting single burner；It can by grid type multiple spot detection technique
Effectively promote the accuracy and validity of concentration field；The actual speed field of the grid type multipoint pressure construction of measurement simultaneously, can be with
The velocity field obtained based on CFD is effectively verified；Burning, which is adjusted, for operations staff provides effective, reliable, comprehensive pass
Key foundation.
Detailed description of the invention
Fig. 1 is boiler structure schematic diagram of the invention；
Fig. 2 is burner hearth and flue structure schematic diagram of the present invention；
Fig. 3 is burner hearth and flue mesh division structure schematic diagram of the present invention；
Fig. 4 is burner region gridding structural schematic diagram of the present invention；
Fig. 5 is the independent method for solving flow chart of the present invention；
Fig. 6 is couple solution method flow diagram of the present invention.
Fig. 7 is sampling schematic diagram of the invention.
Specific embodiment
It elaborates in the following with reference to the drawings and specific embodiments to the present invention:
Specific embodiment is based on boiler furnace grid detection CO and O as shown in Figures 1 to 72The distribution field of concentration constructs
Method, the method include:
A, structure, the operation data of station boiler are obtained: including station boiler dry bottom hopper, burner hearth, furnace arch, furnace nose, burner spray
The detailed dimensions structure of mouth, economizer；Operation data includes at least: the total blast volume of every boiler, total coal amount, each burner
Coal-fired capacity, air quantity, wind-warm syndrome, secondary air temperature, primary air ratio, a secondary air ratio；
B, to high-temperature flue gas CO, O of boiler furnace region wall surface2Concentration and flue gas pressures carry out grid type multiple spot and examine online
It surveys: sampling method being arranged based on burner hearth grid multiple spot, the representative sample of flue gas is taken from burner hearth, i.e., is mounted on contact unit
Predetermined position obtains the fume sample of hearth region wall surface by contact unit, with American Society of Mechanical Engineers (ASME) mark
Quasi- regulation is for smoke sampling: the quantity and arrangement of measuring point depend on the size and shape of flue, sample point interval Ying Buchao
915mm (3 feet) are crossed, and the sum of point used should be no less than at 4 points；To the boiler of different types, in the main burning of boiler furnace
The sample point number of region addressing arrangement, installation is executed in strict accordance with ASME standard, and points selection is as follows: unit capacity
300MW, 16~24 points, analytical unit 4 of installation points；Unit capacity 600MW, 16~24 points of installation points, analytical unit 4
It is a；Unit capacity 1000MW, installation points 16~24 points, analytical unit 4, using analytical unit to the data detected into
Row analytical calculation obtains required data；Boiler furnace region wall surface refers to close to water screen tube to the region of fiery side outer wall, pops one's head in
Unit is made of sampling probe and sampling diversion pipe, and sampling probe is connected by screw thread between the two in the rear end of sampling diversion pipe
It connects, sample gas is exactly drained to sampling probe by the main function for sampling diversion pipe, is equipped with inside sampling probe cylindrical high-precision
Filter is spent, 1 μm of filtering accuracy <, filter efficiency reaches 99%.By the flue gas of filtering, it is single that analysis is entered by heat traced pipeline
Member；It is primary every blowback in four hours in addition, the filter in sampling probe is also connected with blowback compressed air, when blowback,
It will cover in filter surfaces and the dust being deposited in sampling diversion pipe, purged away from sampling diversion pipe, return to furnace
In thorax.In order to obtain more accurate data, the front end for sampling diversion pipe needs to be welded on the fin between two water screen tubes
And flush the sample point for sampling diversion pipe front end face with the front end face of fin, as shown in Figure 7；
C, the flowing, burning of coal dust and heat transfer model in burner hearth are established using the method for Fluid Mechanics Computation, predicts flue gas
CO and O2Motion profile in burner hearth, flue constructs flue gas CO and O2Concentration field:
The calculating grid of burner hearth and flue is established, wherein the computational domain of full boiler is selected to boiler economizer exit, is calculated
Domain is divided into four parts, respectively dry bottom hopper region I, burner region II, burner with up to furnace outlet region III,
Furnace outlet is to economizer exit region IV；Dry bottom hopper I, burner with up to furnace outlet region III, furnace outlet to saving coal
Device exit region IV uses structured grid, and burner region II uses unstructured grid, and wherein structured grid refers to net
All internal points adjacent cells all having the same in lattice region, it can be easily carried out the edge fitting in region, fit
In the calculating that fluid and surface stress are concentrated etc., its major advantage is that the speed that grid generates is fast, high-quality, data knot
Structure is simple, it most typical the disadvantage is that applicable range is narrow, especially as in recent years computer and numerical method
Fast development, requirement of the people to the complexity of domain is higher and higher, in this case, structured grid generation technique
Just seem unable to do what one wishes；Unstructured grid refer to the internal point in net region do not have identical adjacent cells, i.e., with
The connected lattice number of point is different in difference in mesh generation region, as can be seen that structured grid and non-knot from definition
Structure grid has overlapped part, i.e., may include the part of structured grid in unstructured grid, unstructured
Grid can solve the mesh generation of arbitrary shape and any connected region；
Each region is calculated using three-dimensional steady state, using the non-coupled solution of flowing, heat transfer and combustion model, is used in calculating
Finite Volume Method for Air is discrete to differential governing equation and transport equation to carry out numerical solution for Algebraic Equation set, under each operating condition
CO、O2It is Coal-fired capacity, total blast volume, a wind-warm syndrome, secondary air temperature, primary with the actual measured values of flue gas pressures as input condition
Wind rate, secondary air ratio and air leak rate of air curtain are restrained as boundary condition successive ignition until meeting, and predict flue gas CO, O2In burner hearth to province
Motion profile in the outlet of coal device, constructs CO and O in boiler combustion process2The distributed in three dimensions field of concentration；
Wherein, the physical structure and chemical constituent of coal are all sufficiently complex, and burning of coal process can also undergo many complexity
Process, again without obvious boundary between each section, coal dust may undergo two or more processes simultaneously, so the combustion to coal dust
Burning process carries out reasonable description and is not easy to, and it is particularly significant to choose appropriate model to each stage.When numerical value calculates, establish
Volatile matter devolatilization model and coke combustion model make to control to provide the source item in energy-balance equation and component equilibrium equation
Equation group is closed.The precipitation process of volatile matter is also referred to as the pyrolysis of coal in coal, is a sufficiently complex process, is related to object
Manage many processes such as vaporization, chemical bond rupture, the generation of unstable intermediate product and differentiation.Reasonable volatile matter devolatilization model is answered
This can be than more completely describing coal dust in the heating period, product species, quantity and its precipitation rate being therefrom precipitated, so coal
The combustion model of powder uses diffusion/reaction combined control model, and main includes two, respectively diffusion term and heat transfer item, diffusion
Are as follows:
Wherein q is coke consumption rate, and A pre-exponential factor, E is apparent activation energy, POx, sFor the O at particle surface2Partial pressure,
noxTo react order；Conduct heat item are as follows:
Wherein mcFor granular mass, Tg, Tp, TrRespectively gas phase, particle and background radiation temperature, cpFor solid specific heat, fh
For by the reaction heat share of solid sorbent, HreacFor gas-solid reaction heat, εpFor particle blackness, σ is Stefan-Boltzmann constant.
Combustion process of the coal dust in boiler generally can all be related to typical Gas-solid Two-phase Flow, coal dust in the present invention
The flow trace model of grain determines that coal dust tracks granule number according to the number of jet stream using stochastic particle tracking model；It is assumed that penetrating
Flow amount is 20, is meant that and the particle tracked before is divided into 20, and the flow that each particle occupies is exactly original 1/20,
Then different velocity fluctuation values is assigned to this 20 particles, acquires respective particle trajectories.Stochastic Separated Flow Model is with gas phase rapids
Based on the overall view model of stream, to consider gas phase turbulance basic structural feature；It is given at random by the transient state equation of momentum of particle phase
Fluctuation velocity, to consider influence of the DIFFUSION IN TURBULENCE to particle；The tracks of particle are calculated with Monte Carlo method and its along journey
Complexity variation.To the transient state momentum of particle phase on the basis of the way of this method is substantially the equal field when solving gaseous fluid
Equation solution.Stochastic Separated Flow Model is being flowed in research with the particle of the complicated experience such as evaporation, volatilization and out-phase chemical reaction
In variation when have good effect, be at present two-phase lubrication simulation and solid particle combustion simulation in be the most widely used
Model.
It selects when the heat transfer model that the present invention is simulated in the radiant heat transfer to coal dust to obtain coal dust with heat flow method and is
The P1 radiation patterns on basis, compared with discrete ordinates method, Pl radiation patterns take into account radiating scattering, more complex several simulating
There is some superiority when radiation problem in what structure.This method is relatively mature, is used for furnace by many scholars at home and abroad
The simulation of interior radiation.
The PDF model of method selection standard k- ε two-equation turbulence model, non-premixed combustion Passive Scalar of the invention with
And the pollution formation model that FLUENT software provides；For the coupled problem of speed and pressure, using the algorithm of SIMPLE；It is right
The processing of hearth region wall surface uses wall-function method；Simultaneously in order to keep the stability calculated, using the method for deficient relaxation, benefit
With the combination for concentrating solution mode, the flowing computational problem of coal dust inside entire burner hearth is solved, keeps the numerical value being calculated more smart
Really, the concentration field and velocity field constructed is more acurrate, keeps the efficiency of combustion of boiler higher.
One of Fig. 5, Fig. 6 method can be used in specific method for solving, by calculating each items calculated in grid
Parameter obtains concentration field and velocity field in boiler calculation domain.
This construction method has preferable applicability and generalization, can be used for the boiler of various types, to each in boiler
The flue gas CO and O of combustion zone2Concentration field building be a kind of effective method.
Claims (6)
1. being based on boiler furnace grid detection CO and O2The distribution field construction method of concentration, it is characterised in that: the method packet
It includes:
A, obtain station boiler structure, operation data: including station boiler dry bottom hopper, burner hearth, furnace arch, furnace nose, burner nozzle,
The detailed dimensions structure of economizer；Operation data includes at least: the total blast volume of every boiler, total coal amount, the combustion of each burner
Coal amount, air quantity, wind-warm syndrome, secondary air temperature, primary air ratio, a secondary air ratio；
B, to high-temperature flue gas CO, O of boiler furnace region wall surface2Concentration and flue gas pressures carry out grid type multiple spot on-line checking, base
Sampling method is arranged in burner hearth grid multiple spot, and the representative sample of flue gas is taken from burner hearth, i.e., is mounted on contact unit predetermined
Position obtains the fume sample of hearth region wall surface by contact unit；
C, the flowing, burning of coal dust and heat transfer model in burner hearth are established using the method for Fluid Mechanics Computation, predicts flue gas CO and O2
Motion profile in burner hearth, flue constructs flue gas CO and O2Concentration field, establish the calculating grid of burner hearth and flue, wherein
The computational domain of full boiler is selected to boiler economizer exit, and computational domain is divided into four parts, respectively dry bottom hopper region I, combustion
Burner region II, burner are up to furnace outlet region III, furnace outlet to economizer exit region IV；Dry bottom hopper I, combustion
Burner is to use structured grid, burner region up to furnace outlet region III, furnace outlet to economizer exit region IV
II uses unstructured grid；
Each region is calculated using three-dimensional steady state, using the non-coupled solution of flowing, heat transfer and combustion model, using limited in calculating
Volumetric method is discrete to differential governing equation and transport equation to carry out numerical solution for Algebraic Equation set, with CO, O under each operating condition2
With the actual measured values of flue gas pressures as input condition, Coal-fired capacity, total blast volume, a wind-warm syndrome, secondary air temperature, primary air ratio,
Secondary air ratio and air leak rate of air curtain are restrained as boundary condition successive ignition until meeting, and predict flue gas CO, O2In burner hearth to economizer
Motion profile in outlet constructs CO and O in boiler combustion process2The distributed in three dimensions field of concentration.
2. according to claim 1 be based on boiler furnace grid detection CO and O2The distribution field construction method of concentration, feature
Be: step b is with the regulation of American Society of Mechanical Engineers (ASME) standard for smoke sampling: the quantity and arrangement of measuring point take
Certainly in the size and shape of flue, sample point interval should be no more than 915mm, and the sum of point used should be no less than at 4 points；To difference
The boiler of pattern arranges that the sample point number of installation is held in strict accordance with ASME standard in the primary combustion zone addressing of boiler furnace
Row, points selection are as follows: unit capacity 300MW, 16~24 points, analytical unit 4 of installation points；Unit capacity 600MW,
16~24 points, analytical unit 4 of installation points；Unit capacity 1000MW, 16~24 points, analytical unit 4 of installation points.
3. according to claim 2 be based on boiler furnace grid detection CO and O2The distribution field construction method of concentration, feature
Be: the combustion model of coal dust uses diffusion/reaction combined control model, and main includes two, respectively diffusion term and heat transfer
, diffusion term are as follows:
Wherein q is coke consumption rate, and A pre-exponential factor, E is apparent activation energy, POx, sFor the O at particle surface2Partial pressure, noxFor
React order；Conduct heat item are as follows:
Wherein mcFor granular mass, Tg, Tp, TrRespectively gas phase, particle and background radiation temperature, cpFor solid specific heat, fhFor quilt
The reaction heat share of solid sorbent, HreacFor gas-solid reaction heat, εpFor particle blackness, σ is Stefan-Boltzmann constant.
4. according to claim 3 be based on boiler furnace grid detection CO and O2The distribution field construction method of concentration, feature
Be: the flow trace model of pulverized coal particle determines coal dust tracking according to the number of jet stream using stochastic particle tracking model
Grain number；It is assumed that jet stream number is 20, it is meant that and the particle tracked before is divided into 20, the flow that each particle occupies is exactly
Originally 1/20, different velocity fluctuation values then is assigned to this 20 particles, acquires respective particle trajectories.
5. according to claim 4 be based on boiler furnace grid detection CO and O2The distribution field construction method of concentration, feature
It is: selects when being simulated the heat transfer model to obtain coal dust in the radiant heat transfer to coal dust by P1 spoke based on heat flow method
Model is penetrated, compared with discrete ordinates method, Pl radiation patterns take into account radiating scattering.
6. according to claim 5 be based on boiler furnace grid detection CO and O2The distribution field construction method of concentration, feature
Be: the PDF model of the method selection standard k- ε two-equation turbulence model, non-premixed combustion Passive Scalar and
The pollution formation model that FLUENT software provides；For the coupled problem of speed and pressure, using the algorithm of SIMPLE；To furnace
The processing of thorax region wall surface uses wall-function method；Simultaneously using the method for owing relaxation.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610981224.4A CN106838979B (en) | 2016-11-08 | 2016-11-08 | Based on boiler furnace grid detection CO and O2The distribution field construction method of concentration |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610981224.4A CN106838979B (en) | 2016-11-08 | 2016-11-08 | Based on boiler furnace grid detection CO and O2The distribution field construction method of concentration |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106838979A CN106838979A (en) | 2017-06-13 |
| CN106838979B true CN106838979B (en) | 2019-03-12 |
Family
ID=59146330
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610981224.4A Active CN106838979B (en) | 2016-11-08 | 2016-11-08 | Based on boiler furnace grid detection CO and O2The distribution field construction method of concentration |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106838979B (en) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110645592A (en) * | 2019-10-12 | 2020-01-03 | 西安康桥能源技术有限责任公司 | Combustion system improvement method based on multiphase partition coupling |
| CN111396920A (en) * | 2020-03-26 | 2020-07-10 | 华润电力技术研究院有限公司 | Thermal power generating unit and boiler combustion monitoring method and system based on CO measurement |
| CN112240566B (en) * | 2020-11-09 | 2023-10-24 | 国电铜陵发电有限公司 | Online adjustment system and method for offset firing of boiler |
Family Cites Families (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7010461B2 (en) * | 2004-02-09 | 2006-03-07 | General Electric Company | Method and system for real time reporting of boiler adjustment using emission sensor data mapping |
| US9310347B2 (en) * | 2010-11-16 | 2016-04-12 | General Electric Company | Methods and systems for analyzing combustion system operation |
| CN103454105B (en) * | 2013-08-22 | 2016-06-08 | 广东电网公司电力科学研究院 | Biomass recirculating fluidized bed boiler thermal characteristic measuring method and system |
| CN103955599B (en) * | 2014-03-28 | 2017-02-01 | 广东电网公司电力科学研究院 | Acquisition method of combustion characteristic under different circumference air quantities after boiler improvement |
-
2016
- 2016-11-08 CN CN201610981224.4A patent/CN106838979B/en active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN106838979A (en) | 2017-06-13 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Guo et al. | Experimental and numerical investigations on oxy-coal combustion in a 35 MW large pilot boiler | |
| CN106548032B (en) | Based on boiler flue grid detection flue gas CO and O2The distribution field method for reconstructing of concentration | |
| Weber et al. | Fly ash deposition modelling: Requirements for accurate predictions of particle impaction on tubes using RANS-based computational fluid dynamics | |
| Chen et al. | Numerical investigations on different tangential arrangements of burners for a 600 MW utility boiler | |
| Yin et al. | Characterizing and modeling of an 88 MW grate-fired boiler burning wheat straw: Experience and lessons | |
| CN111695249B (en) | Prediction method for heat efficiency of gas boiler | |
| CN106838979B (en) | Based on boiler furnace grid detection CO and O2The distribution field construction method of concentration | |
| Zhang et al. | Analysis and neural network prediction of combustion stability for industrial gases | |
| CN103148506A (en) | Secondary air distribution method for pulverized coal boiler with swirling combustion of front-back hedging | |
| CN107543199A (en) | A kind of pulverized-coal fired boiler online detection of primary air and burning optimization management system | |
| Guo et al. | Effects of gas and particle radiation on IFRF 2.5 MW swirling flame under oxy-fuel combustion | |
| Guo et al. | Experimental and numerical investigations on heat transfer characteristics of a 35MW oxy-fuel combustion boiler | |
| Zhang et al. | Flow behavior of high-temperature flue gas in the heat transfer chamber of a pilot-scale coal-water slurry combustion furnace | |
| Li et al. | Particle sticking behavior near the throat of a low-NOx axial-swirl coal burner | |
| CN103148508A (en) | Secondary air distribution method and system for pulverized coal boiler with swirling combustion of front-back hedging | |
| Darbandi et al. | Numerical study on NOx reduction in a large-scale heavy fuel oil-fired boiler using suitable burner adjustments | |
| Jovanović et al. | Mathematical modelling of swirl oxy-fuel burner flame characteristics | |
| Wang et al. | Mechanism modeling of optimal excess air coefficient for operating in coal fired boiler | |
| CN103343978A (en) | Method for predicting and diagnosing pulverized coal flame stability | |
| Ben-Mansour et al. | Oxy-fuel combustion in a 600 MW gaseous fuel tangentially fired boiler | |
| Chen et al. | Numerical study on the influence of platen super-heater on combustion characteristics in a 600 MW tangentially fired pulverized-coal boiler | |
| Stewart et al. | Creating an Advanced Sensor Network to calculate real-time, mass-weighted flue gas composition and air heater leakage of a coal-fired utility boiler under dynamic operating conditions | |
| Fomenko et al. | Numerical study of the methane-air combustion in the direct-flow burner of the boiler 300 MW TGMP-314 boiler | |
| Gil et al. | Application of numerical analysis for physical and chemical combustion processes in design of boiler units | |
| Shen et al. | Industrial Experimental and Numerical Simulation Study on the Combustion and NO x Formation Characteristics in a 600MWe Utility Boiler with a Novel Swirl Burner Burning Bituminous Coal |
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 | ||
| TR01 | Transfer of patent right | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20190711 Address after: 065000 No. 1 Building, No. 20 Qingxiang Road, Guangyang Economic Development Zone, Langfang City, Hebei Province Patentee after: WATSON ENERGY TECHNOLOGY (LANGFANG) Co.,Ltd. Address before: 065000 No. 10 Changxiang Road, Guangyang Economic Development Zone, Langfang City, Hebei Province Patentee before: WATSON MEASUREMENT & CONTROL TECHNOLOGY (HEBEI) Co.,Ltd. |
|
| PE01 | Entry into force of the registration of the contract for pledge of patent right |
Denomination of invention: Method of constructing distribution field based on grid detection of CO and O2concentrations in boiler furnace Effective date of registration: 20221223 Granted publication date: 20190312 Pledgee: China Everbright Bank Limited by Share Ltd. Langfang branch Pledgor: WATSON ENERGY TECHNOLOGY (LANGFANG) Co.,Ltd. Registration number: Y2022980028822 |
|
| PE01 | Entry into force of the registration of the contract for pledge of patent right | ||
| PC01 | Cancellation of the registration of the contract for pledge of patent right |
Date of cancellation: 20231212 Granted publication date: 20190312 Pledgee: China Everbright Bank Limited by Share Ltd. Langfang branch Pledgor: WATSON ENERGY TECHNOLOGY (LANGFANG) Co.,Ltd. Registration number: Y2022980028822 |
|
| PC01 | Cancellation of the registration of the contract for pledge of patent right | ||
| PE01 | Entry into force of the registration of the contract for pledge of patent right |
Denomination of invention: A Method for Constructing Distribution Fields of CO and O2Concentrations Based on Boiler Furnace Grid Detection Effective date of registration: 20231215 Granted publication date: 20190312 Pledgee: China Everbright Bank Limited by Share Ltd. Langfang branch Pledgor: WATSON ENERGY TECHNOLOGY (LANGFANG) Co.,Ltd. Registration number: Y2023110000529 |
|
| PE01 | Entry into force of the registration of the contract for pledge of patent right |