CN106838979A - 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
- CN106838979A CN106838979A CN201610981224.4A CN201610981224A CN106838979A CN 106838979 A CN106838979 A CN 106838979A CN 201610981224 A CN201610981224 A CN 201610981224A CN 106838979 A CN106838979 A CN 106838979A
- Authority
- CN
- China
- Prior art keywords
- boiler
- concentration
- flue gas
- furnace
- region
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000001514 detection method Methods 0.000 title claims abstract description 19
- 238000010276 construction Methods 0.000 title claims abstract description 15
- 238000009826 distribution Methods 0.000 title claims abstract description 15
- 238000000034 method Methods 0.000 claims abstract description 54
- 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 34
- 239000002817 coal dust Substances 0.000 claims abstract description 24
- 238000012546 transfer Methods 0.000 claims abstract description 20
- 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 17
- 238000006243 chemical reaction Methods 0.000 claims description 14
- 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
- 230000033001 locomotion Effects 0.000 claims description 5
- 239000000779 smoke Substances 0.000 claims description 5
- 210000000038 chest Anatomy 0.000 claims description 4
- 239000000571 coke Substances 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
- 239000003517 fume Substances 0.000 claims description 3
- 238000009434 installation Methods 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
- 230000008859 change Effects 0.000 description 5
- 238000010586 diagram Methods 0.000 description 5
- 238000005259 measurement Methods 0.000 description 4
- 230000000694 effects Effects 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 239000000843 powder Substances 0.000 description 3
- 238000004088 simulation Methods 0.000 description 3
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 230000005611 electricity Effects 0.000 description 2
- 238000001914 filtration Methods 0.000 description 2
- 239000008187 granular material Substances 0.000 description 2
- 239000000203 mixture 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
- 238000013528 artificial neural network Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 238000004939 coking Methods 0.000 description 1
- 239000000470 constituent Substances 0.000 description 1
- 230000007797 corrosion Effects 0.000 description 1
- 238000005260 corrosion Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000003111 delayed 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
- 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
- 239000004744 fabric 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
- 238000003825 pressing Methods 0.000 description 1
- 239000000047 product Substances 0.000 description 1
- 238000000197 pyrolysis Methods 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
Abstract
The invention discloses 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 high-temperature flue gas CO, O2The grid type multiple spot on-line checking of concentration and flue gas pressures, flowing, burning and the heat transfer model of coal dust in burner hearth is set up using the method for Fluid Mechanics Computation, to predict flue gas CO and O2Movement locus in burner hearth, flue, and then build flue gas CO and O2Concentration field, for boiler combustion adjustment provides effectively reliable crucial foundation.
Description
Technical field
The invention belongs to boiler combustion efficiency detection method technical field, it is related to CO and O2The distribution field structure side of concentration
Method, and in particular to based on boiler furnace grid detection CO and O2The distribution field construction method of concentration.
Background technology
As the development and progress of society, thermal power plant face the increasing pressure, electricity needs is continuously increased, to generating electricity
The requirement of quality is increasingly harsh, and traditional voltage, frequency has not been unique performance assessment criteria, and environmental factor is gradually included,
Form new criterion.Ature of coal fluctuates and Large-scale machine set increasingly frequently participates in frequency modulation peak regulation, full working scope varying duty fortune
Row and long-term middle-low load operation all propose requirement higher to unit operation level.At present, power plant is increasing
Face from power network, environment and the pressure of itself, energy-conservation, consumption reduction, emission reduction cry it is surging, burning optimization technology increasingly obtains weight
Depending on.Can the height of boiler combustion efficiency, the fuel into burner hearth fully burn, and will directly affect the whole efficiency and warp of boiler
Ji property.Therefore realize that most efficiently burning is the effective way of raising boiler whole efficiency to boiler by Combustion System.
Generating plant pulverized coal boiler Combustion System design basic ideas are, with suitable excess air coefficient as principle, to realize boiler
Burning optimum control is to ensure suitable excess air coefficient.When the quantity combusted of boiler changes, correspondingly change air output,
To ensure suitable excess air coefficient, the imperfect combustion loss of boiler and heat loss due to exhaust gas are reduced, make boiler operatiopn in optimal
Combustion conditions.Simple measurement O2The quality of coal dust and air troubled water in boiler can not be reflected, even if O2Enough, if in the presence of
Mix the reason such as bad, incomplete combustion loss can also increased.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
Skewness easily makes oxygen amount detected value not possess representativeness.In addition, oxygen amount detected value in itself due to leaking out, the original such as zero drift
Because of the true oxygen amount value that can not reflect in flue gas.For these reasons, in boiler actual moving process, operations staff can only
Burning is controlled with larger excess air coefficient, it is very uneconomical, cause boiler operatiopn to be susceptible to burner hearth high temperature corrosion, flame
The operation puzzlement such as relatively low of deflection, combustion lag, coking and slagging, efficiency of combustion.How to find a perfect detection means and come right
Operating boiler is adjusted, and is very must to solve the problems, such as boiler operatiopn and improve the actual operating efficiency of boiler
Want.
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 air pressure, pressure fan electric current, air-introduced machine electric current, bellows and burner hearth differential pressure, feedwater and smoke evacuation temperature
11 parameters such as degree, using flue gas oxygen content as output parameter, are set up flue gas and are contained as input signal using feedforward neural network
Oxygen amount soft-sensing model.Under various different load operating modes, multigroup actual operating data is measured, extract 1~2 group of data as mould
Type detection is used, and remaining is used as training sample.This kind of method solves that oxygen content measurement is delayed, error is big asks to a certain extent
Topic.But, the accuracy of this hard measurement result depends on the integrality of Boundary Test and comprehensive, often due to as-fired coal coal
The fluctuation of matter, the change of operating mode cause error to increase, and during firing optimization, the application of soft-measuring technique is difficult on a large scale
Effectively promote.
The content of the invention
Defect in order to overcome prior art of the invention, devises based on boiler furnace grid detection CO and O2Concentration point
Cloth construction method, flowing, burning and the heat transfer model of coal dust in burner hearth is set up using the method for Fluid Mechanics Computation, to predict
Flue gas CO and O2Movement locus in burner hearth, flue, and then build flue gas CO and O2Concentration field, be boiler combustion adjustment carry
It is provided with the reliable key foundation of effect.
The concrete technical scheme taken of the present invention is：Based on boiler furnace grid detection CO and O2The distribution field structure of concentration
Construction method, it is important to：Described method includes：
A, the structure, the peration data that obtain station boiler；
B, high-temperature flue gas CO, O to boiler furnace region wall2Concentration and flue gas pressures carry out grid type multiple spot and examine online
Survey；
C, the flowing that coal dust in burner hearth is set up using the method for Fluid Mechanics Computation, burning and heat transfer model, predict flue gas
CO and O2Movement locus in burner hearth, flue, builds flue gas CO and O2Concentration field.
Methods described is specifically included：
A, the structure, the peration data that obtain station boiler：Including station boiler furnace hopper, burner hearth, furnace arch, furnace nose, burner spray
Mouth, the detailed dimensions structure of economizer；Peration data at least includes：Every the total blast volume of 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, based on burner hearth grid multiple spot arrange sampling method, the representative sample of flue gas is taken from burner hearth, will pop one's head in list
Unit is arranged on precalculated position, and the fume sample of hearth region wall is obtained by contact unit；
C, the calculating grid for setting up burner hearth and flue, wherein the computational fields of full boiler are selected to boiler economizer exit, meter
Calculate domain and be divided into four parts, respectively furnace hopper region I, burner region II, burner are with up to furnace outlet region
III, furnace outlet are to economizer exit region IV；Furnace hopper I, burner with up to furnace outlet region III, furnace outlet extremely
Economizer exit region IV uses structured grid, burner region II to use unstructured grid；
Each region is calculated using three-dimensional steady state, using flowing, heat transfer and the non-coupled solution of combustion model, is used in calculating
Finite Volume Method for Air is discrete to differential governing equation and transport equation for Algebraic Equation set carries out numerical solution, with each operating mode
CO、O2With the actual measured value of flue gas pressures as input condition, Coal-fired capacity, total blast volume, wind-warm syndrome, secondary air temperature, once
Wind rate, secondary air ratio and air leak rate of air curtain restrain as boundary condition successive ignition until meeting, prediction flue gas CO, O2In burner hearth to province
Movement locus in the outlet of coal device, builds CO and O in boiler combustion process2The distributed in three dimensions of concentration.
Step b is the regulation with ASME (ASME) standard for smoke sampling：The quantity of measuring point and
Arrangement depends on the size and dimension 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 is performed, and points selection is as follows：Unit capacity 300MW, installs 16~24 points of points, analytic unit 4；Unit capacity
600MW, installs 16~24 points of points, analytic unit 4；Unit capacity 1000MW, installs 16~24 points of points, analytic unit 4
It is individual.
The combustion model of coal dust uses diffusion/reaction combined control model, main to include two, respectively diffusion term and biography
Hot-activity, diffusion term is：
Wherein q is coke consumption speed, and A pre-exponential factors, E is apparent activation energy, POx, sIt is the O at particle surface2Partial pressure,
noxIt is reaction exponent number；Heat transfer item be：
Wherein mcIt is granular mass, Tg, Tp, TrRespectively gas phase, particle and background radiation temperature, cpIt is solid specific heat, fh
It is by the reaction heat share of solid sorbent, HreacIt is gas-solid reaction heat, εpIt is particle blackness, σ is Stefan-Boltzmann constant.
The flow trace model of pulverized coal particle, using stochastic particle tracking model, the number according to jet determine coal dust with
Track granule number；It is assumed that jet number is 20, it is meant that and the particle for tracking 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, tries to achieve respective particle trajectories.
When the radiant heat transfer to coal dust is simulated the heat transfer model to draw coal dust from based on heat flow method
P1 radiation patterns, compared with discrete ordinates method, Pl radiation patterns take radiating scattering into account.
Described method selection standard k- ε two-equation turbulence models, the PDF models of non-premixed combustion Passive Scalar and
The pollution formation model that FLUENT softwares are provided；For speed and the coupled problem of pressure, using the algorithm of SIMPLE；To stove
The treatment of thorax region wall uses wall-function method；Simultaneously using the method for owing to relax.
The beneficial effects of the invention are as follows：By to boiler furnace region wall high-temperature flue gas CO, O2Concentration and flue gas pressures
Grid type multiple spot on-line checking, using the method for Fluid Mechanics Computation set up coal dust in burner hearth flowing, burning and heat transfer mould
Type, to predict flue gas CO and O2Movement locus in burner hearth, flue, and then build flue gas CO and O2Concentration field, be boiler combustion
Burn adjustment and effectively reliable key foundation is provided.
Use the method based on Computational Fluid Dynamics, high efficiency and time conservation, it is possible to use existing many model libraries carry out excellent
Change is processed；Reaction flue gas CO, O directly perceived2Movement locus and the regularity of distribution in burner hearth, can clearly grasp wind powder in burner hearth and mix
Conjunction situation, is conducive to the wind powder for targetedly individually adjusting single burner to match；Can by grid type multiple spot detection technique
Effectively lift the accuracy and validity of concentration field；The actual speed of the grid type multipoint pressure construction of measurement simultaneously, can be with
Velocity field to being drawn based on CFD is effectively verified；For operations staff's adjustment burning provides effective, reliable, comprehensive closing
Key foundation.
Brief description of the drawings
Fig. 1 is boiler structure schematic diagram of the invention；
Fig. 2 is burner hearth of the present invention and flue structure schematic diagram；
Fig. 3 is burner hearth of the present invention and flue mesh division structure schematic diagram；
Fig. 4 is burner region gridding structural representation of the present invention；
Fig. 5 is independent method for solving flow chart of the invention；
Fig. 6 is couple solution method flow diagram of the present invention.
Fig. 7 is sampling schematic diagram of the invention.
Specific embodiment
The present invention is elaborated with specific embodiment below in conjunction with the accompanying drawings：
Specific embodiment, as shown in Figures 1 to 7, based on boiler furnace grid detection CO and O2The distribution field of concentration builds
Method, described method includes：
A, the structure, the peration data that obtain station boiler：Including station boiler furnace hopper, burner hearth, furnace arch, furnace nose, burner spray
Mouth, the detailed dimensions structure of economizer；Peration data at least includes：Every the total blast volume of 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, high-temperature flue gas CO, O to boiler furnace region wall2Concentration and flue gas pressures carry out grid type multiple spot and examine online
Survey：Based on burner hearth grid multiple spot arrange sampling method, the representative sample of flue gas is taken from burner hearth, will contact unit be arranged on
Precalculated position, the fume sample of hearth region wall is obtained by contact unit, is marked with ASME (ASME)
Accurate regulation is for smoke sampling：The quantity and arrangement of measuring point depend on the size and dimension of flue, and sample point interval should not surpass
915mm (3 feet) is 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
Region addressing arrangement, the sample point number of installation is performed in strict accordance with ASME standards, and points selection is as follows：Unit capacity
300MW, installs 16~24 points of points, analytic unit 4；Unit capacity 600MW, installs 16~24 points of points, analytic unit 4
It is individual；Unit capacity 1000MW, installs 16~24 points of points, and analytic unit 4 is carried out using analytic unit to the data for detecting
Analytical calculation, draws required data；Boiler furnace region wall refers to the region for pressing close to water screen tube to fiery side outer wall, and probe is single
Unit is made up of sampling probe and sampling mozzle, and sampling probe is connected through a screw thread between the two in the rear end of sampling mozzle,
Sample gas is exactly drained to sampling probe by the Main Function for sampling mozzle, cylinder is provided with inside sampling probe and high-precision is spent
Filter, 1 μm of filtering accuracy ＜, filter efficiency reaches 99%.By the flue gas for filtering, analytic unit is entered by heat traced pipeline；
In addition, the filter in sampling probe is also connected with blowback compressed air, and every blowback in four hours once, when blowback, will
Cover in filter surfaces and the dust being deposited in sampling mozzle, purged away from sampling mozzle, return to burner hearth
It is interior.In order to obtain more accurately data, sampling the front end of mozzle needs on the fin that is welded between two water screen tubes simultaneously
The sample point of sampling mozzle front end face is set to be flushed with the front end face of fin, as shown in Figure 7；
C, the flowing that coal dust in burner hearth is set up using the method for Fluid Mechanics Computation, burning and heat transfer model, predict flue gas
CO and O2Movement locus in burner hearth, flue, builds flue gas CO and O2Concentration field：
The calculating grid of burner hearth and flue is set up, wherein the computational fields of full boiler are selected to boiler economizer exit, is calculated
Domain is divided into four parts, respectively furnace hopper region I, burner region II, burner with up to furnace outlet region III,
Furnace outlet is to economizer exit region IV；Furnace hopper I, burner are with up to furnace outlet region III, furnace outlet to saving coal
Device exit region IV uses structured grid, burner region II to use unstructured grid, and wherein structured grid refers to net
All of internal point all has identical adjacent cells in lattice region, and it can be easily carried out the edge fitting in region, fits
In the calculating of fluid and the medium aspect of surface stress collection, its major advantage be mess generation speed is fast, quality is good, data knot
Structure is simple, and the most typical of it has the disadvantage that applicable scope is narrow, especially as computer in recent years and numerical method
Fast development, requirement more and more higher of the people to the complexity of domain, in this case, structured grid generation technique
Just seem unable to do what one wishes；Unstructured grid refers to internal point in net region without 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 may the part comprising structured grid, destructuring in having overlapped part, i.e. unstructured grid
Grid can solve the mesh generation of arbitrary shape and any connected region；
Each region is calculated using three-dimensional steady state, using flowing, heat transfer and the non-coupled solution of combustion model, is used in calculating
Finite Volume Method for Air is discrete to differential governing equation and transport equation for Algebraic Equation set carries out numerical solution, with each operating mode
CO、O2With the actual measured value of flue gas pressures as input condition, Coal-fired capacity, total blast volume, wind-warm syndrome, secondary air temperature, once
Wind rate, secondary air ratio and air leak rate of air curtain restrain as boundary condition successive ignition until meeting, prediction flue gas CO, O2In burner hearth to province
Movement locus in the outlet of coal device, builds CO and O in boiler combustion process2The distributed in three dimensions of concentration；
Wherein, the physical arrangement and chemical constituent of coal are all sufficiently complex, and burning of coal process can also experience many complexity
Process, again without obvious boundary between each several part, coal dust may simultaneously experience two or more processes, so the combustion to coal dust
Burning process carries out rational description and is not easy to, and it is particularly significant to choose appropriate model to each stage.In numerical computations, set up
Volatile matter devolatilization model and coke combustion model, to provide the source item in energy-balance equation and component equilibrium equation, make control
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 thing
Many processes such as reason vaporization, chemical bond rupture, the generation of unstable intermediate product and differentiation.Rational volatile matter devolatilization model should
This can be than more completely description coal dust in heating period, the product species for therefrom separating out, quantity and its precipitation rate, so coal
The combustion model of powder uses diffusion/reaction combined control model, main to include two, respectively diffusion term and heat transfer item, diffusion
Xiang Wei：
Wherein q is coke consumption speed, and A pre-exponential factors, E is apparent activation energy, POx, sIt is the O at particle surface2Partial pressure,
noxIt is reaction exponent number；Heat transfer item be：
Wherein mcIt is granular mass, Tg, Tp, TrRespectively gas phase, particle and background radiation temperature, cpIt is solid specific heat, fh
It is by the reaction heat share of solid sorbent, HreacIt is gas-solid reaction heat, εpIt is particle blackness, σ is Stefan-Boltzmann constant.
Combustion process of the coal dust in boiler typically can all be related to typical Gas-solid Two-phase Flow, coal dust in the present invention
The flow trace model of grain, using stochastic particle tracking model, the number according to jet determines coal dust tracking granule number；It is assumed that penetrating
Flow amount is 20, and it is meant that and the particle for tracking before is divided into 20, and the flow that each particle occupies is exactly original 1/20,
Then different velocity fluctuation values are assigned to this 20 particles, tries to achieve respective particle trajectories.Stochastic Separated Flow Model is rapid with gas phase
Based on the overall view model of stream, to consider gas phase turbulance basic structural feature；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 change.The way of this method is substantially when gaseous fluid is solved on the basis of equal field to the transient state momentum of particle phase
Equation solution.Stochastic Separated Flow Model is accompanied by the complicated particles for experiencing such as evaporation, volatilization and out-phase chemical reaction in stream in research
In change when have good effect, be at present two-phase lubrication simulation and solid particle combustion simulation in be the most widely used
Model.
The present invention is selected when the radiant heat transfer to coal dust is simulated with the heat transfer model for drawing coal dust
The P1 radiation patterns on basis, compared with discrete ordinates method, Pl radiation patterns take radiating scattering into account, more complicated several simulating
There is some superiority during radiation problem in what structure.This method is relatively ripe, is used for stove by domestic and international many scholars
The simulation of interior radiation.
Method of the present invention selection standard k- ε two-equation turbulence models, the PDF models of non-premixed combustion Passive Scalar and
The pollution formation model that FLUENT softwares are provided；For speed and the coupled problem of pressure, using the algorithm of SIMPLE；To stove
The treatment of thorax region wall uses wall-function method；Simultaneously for the stability for keeping calculating, using lax method is owed, utilize
The combination of solution mode is concentrated, the flowing computational problem of whole burner hearth inside coal dust is solved, makes the numerical value for calculating more accurate,
The concentration field and velocity field for constructing are more accurate, make the efficiency of combustion of boiler higher.
Specific method for solving can be using the one kind in Fig. 5, Fig. 6 method, by calculating each item in calculating grid
Parameter, obtains the 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 build be a kind of effective method.
Claims (7)
1. boiler furnace grid detection CO and O are based on2The distribution field construction method of concentration, it is characterised in that：Described method bag
Include：
A, the structure, the peration data that obtain station boiler；
B, high-temperature flue gas CO, O to boiler furnace region wall2Concentration and flue gas pressures carry out grid type multiple spot on-line checking；
C, the flowing that coal dust in burner hearth is set up using the method for Fluid Mechanics Computation, burning and heat transfer model, prediction flue gas CO and O2
Movement locus in burner hearth, flue, builds flue gas CO and O2Concentration field.
2. according to claim 1 based on boiler furnace grid detection CO and O2The distribution field construction method of concentration, its feature
It is：Methods described is specifically included：
A, the structure, the peration data that obtain station boiler：Including station boiler furnace hopper, burner hearth, furnace arch, furnace nose, burner nozzle,
The detailed dimensions structure of economizer；Peration data at least includes：Every the total blast volume of 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, based on burner hearth grid multiple spot arrange sampling method, the representative sample of flue gas is taken from burner hearth, will contact unit peace
Mounted in precalculated position, the fume sample of hearth region wall is obtained by contact unit；
C, the calculating grid for setting up burner hearth and flue, wherein the computational fields of full boiler are selected to boiler economizer exit, computational fields
Four parts are divided into, respectively furnace hopper region I, burner region II, burner are with up to furnace outlet region III, stove
Thorax is exported to economizer exit region IV；Furnace hopper I, burner are with up to furnace outlet region III, furnace outlet to economizer
Exit region IV uses structured grid, burner region II to use unstructured grid；
Each region is calculated using three-dimensional steady state, using flowing, heat transfer and the non-coupled solution of combustion model, using limited in calculating
Volumetric method is discrete to differential governing equation and transport equation for Algebraic Equation set carries out numerical solution, with CO, O under each operating mode2
With the actual measured value of flue gas pressures as input condition, Coal-fired capacity, total blast volume, wind-warm syndrome, secondary air temperature, primary air ratio,
Secondary air ratio and air leak rate of air curtain restrain as boundary condition successive ignition until meeting, prediction flue gas CO, O2In burner hearth to economizer
Movement locus in outlet, builds CO and O in boiler combustion process2The distributed in three dimensions of concentration.
3. according to claim 2 based on boiler furnace grid detection CO and O2The distribution field construction method of concentration, its feature
It is：Step b is with the regulation of ASME (ASME) standard for smoke sampling：The quantity and arrangement of measuring point take
Certainly in the size and dimension 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 standards in the primary combustion zone addressing of boiler furnace
OK, points selection is as follows：Unit capacity 300MW, installs 16~24 points of points, analytic unit 4；Unit capacity 600MW,
16~24 points of points, analytic unit 4 are installed；Unit capacity 1000MW, installs 16~24 points of points, analytic unit 4.
4. according to claim 3 based on boiler furnace grid detection CO and O2The distribution field construction method of concentration, its feature
It is：The combustion model of coal dust uses diffusion/reaction combined control model, main to include two, respectively diffusion term and heat transfer
, diffusion term is：
Wherein q is coke consumption speed, and A pre-exponential factors, E is apparent activation energy, POx, sIt is the O at particle surface2Partial pressure, noxFor
Reaction exponent number；Heat transfer item be：
Wherein mcIt is granular mass, Tg, Tp, TrRespectively gas phase, particle and background radiation temperature, cpIt is solid specific heat, fhBe by
The reaction heat share of solid sorbent, HreacIt is gas-solid reaction heat, εpIt is particle blackness, σ is Stefan-Boltzmann constant.
5. according to claim 4 based on boiler furnace grid detection CO and O2The distribution field construction method of concentration, its feature
It is：The flow trace model of pulverized coal particle, using stochastic particle tracking model, the number according to jet determines coal dust tracking
Grain number；It is assumed that jet number is 20, it is meant that and the particle for tracking before is divided into 20, and the flow that each particle occupies is exactly
Originally 1/20, different velocity fluctuation values then are assigned to this 20 particles, try to achieve respective particle trajectories.
6. according to claim 5 based on boiler furnace grid detection CO and O2The distribution field construction method of concentration, its feature
It is：When the radiant heat transfer to coal dust is simulated the heat transfer model to draw coal dust from the P1 spokes based on heat flow method
Model is penetrated, compared with discrete ordinates method, Pl radiation patterns take radiating scattering into account.
7. according to claim 6 based on boiler furnace grid detection CO and O2The distribution field construction method of concentration, its feature
It is：Described method selection standard k- ε two-equation turbulence models, the PDF models of non-premixed combustion Passive Scalar and
The pollution formation model that FLUENT softwares are provided；For speed and the coupled problem of pressure, using the algorithm of SIMPLE；To stove
The treatment of thorax region wall uses wall-function method；Simultaneously using the method for owing to relax.
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 true CN106838979A (en) | 2017-06-13 |
| CN106838979B 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) |
Cited By (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 |
| CN112240566A (en) * | 2020-11-09 | 2021-01-19 | 国电铜陵发电有限公司 | Boiler partial combustion online adjusting system and method |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20060136161A1 (en) * | 2004-02-09 | 2006-06-22 | General Electric Company | Method and system for real time reporting of boiler adjustment using emission sensor data mapping |
| US20120122040A1 (en) * | 2010-11-16 | 2012-05-17 | General Electric Company | Methods and Systems for Analyzing Combustion System Operation |
| CN103454105A (en) * | 2013-08-22 | 2013-12-18 | 广东电网公司电力科学研究院 | Method and system for measuring thermal state characteristics of biomass circulating fluidized bedboiler |
| CN103955599A (en) * | 2014-03-28 | 2014-07-30 | 广东电网公司电力科学研究院 | 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
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20060136161A1 (en) * | 2004-02-09 | 2006-06-22 | General Electric Company | Method and system for real time reporting of boiler adjustment using emission sensor data mapping |
| US20120122040A1 (en) * | 2010-11-16 | 2012-05-17 | General Electric Company | Methods and Systems for Analyzing Combustion System Operation |
| CN103454105A (en) * | 2013-08-22 | 2013-12-18 | 广东电网公司电力科学研究院 | Method and system for measuring thermal state characteristics of biomass circulating fluidized bedboiler |
| CN103955599A (en) * | 2014-03-28 | 2014-07-30 | 广东电网公司电力科学研究院 | Acquisition method of combustion characteristic under different circumference air quantities after boiler improvement |
Cited By (4)
| 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 |
| CN112240566A (en) * | 2020-11-09 | 2021-01-19 | 国电铜陵发电有限公司 | Boiler partial combustion online adjusting system and method |
| CN112240566B (en) * | 2020-11-09 | 2023-10-24 | 国电铜陵发电有限公司 | Online adjustment system and method for offset firing of boiler |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106838979B (en) | 2019-03-12 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Chen et al. | Numerical investigations on different tangential arrangements of burners for a 600 MW utility boiler | |
| CN111695249B (en) | Prediction method for heat efficiency of gas boiler | |
| CN106548032A (en) | Based on boiler flue grid detection flue gas CO and O2The distribution field method for reconstructing of concentration | |
| Yu et al. | Numerical investigation of combustion optimization in a tangential firing boiler considering steam tube overheating | |
| CN111339716B (en) | Boiler high-temperature flue gas flow field online agent model construction method | |
| 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. | Experimental and numerical investigations on heat transfer characteristics of a 35MW oxy-fuel combustion boiler | |
| CN106838979A (en) | Based on boiler furnace grid detection CO and O2The distribution field construction method of concentration | |
| CN103148508A (en) | Secondary air distribution method and system for pulverized coal boiler with swirling combustion of front-back hedging | |
| Deng et al. | Effects of flue gas recirculation on combustion and heat flux distribution in 660 MW double-reheat tower-type boiler | |
| Bianco et al. | CFD analysis and design optimization of an air manifold for a biomass boiler | |
| CN116776770B (en) | CFD numerical simulation coupling BP neural network based method for relieving high-temperature corrosion of boiler | |
| CN113074459A (en) | Boiler control optimization method and system | |
| Holkar et al. | CFD Anlysis of Pulverised-Coal Combustion of Burner Used In Furnace with Different Radiation Models | |
| Vondál et al. | Wall heat transfer in gas-fired furnaces: Effect of radiation modelling | |
| CN108875165A (en) | Boiler characteristics scaling method based on operation data | |
| Zheng et al. | Numerical and experimental investigations on the performance of a 300 MW pulverized coal furnace | |
| Xu et al. | An advanced model to assess fouling and slagging in coal fired boilers | |
| Zhu et al. | A numerical study of ash deposition characteristics in a 660mw supercritical tangential 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 | |
| CN105043798A (en) | Method for measuring thermal state characteristics of alcohol-based fuel combustor | |
| Filkoski et al. | Computational fluid dynamics technique as a tool for description of the phenomena occurring in pulverized coal combustion systems | |
| Purnomo et al. | Numerical simulation of coal particle size (fineness) effect to combustion characteristics of sub-critical pulverized coal boiler 600 MW capacity | |
| Gil et al. | Application of numerical analysis for physical and chemical combustion processes in design of boiler units |
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 |
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. |
|
| TR01 | Transfer of patent right | ||
| PE01 | Entry into force 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: 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 |
|
| PC01 | Cancellation 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 |
|
| PE01 | Entry into force 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 |