The content of the invention
The technical problems to be solved by the invention are to provide temperature in a kind of burner hearth based on burner hearth Multi sectional slagging situation
Flexible measurement method is distributed, on the basis of station boiler measuring point is not increased, burner hearth is divided into multiple sections according to ignition quality,
The section of burner hearth cigarette temperature for carrying out different sections by the burner hearth inside smoke soft-sensing model established is calculated, and provides stove simultaneously
The real-time slagging situation of different sections inside thorax.
In order to solve the above technical problems, the technical solution adopted in the present invention is:
Temperature Distribution flexible measurement method in a kind of burner hearth based on burner hearth Multi sectional slagging situation, it is characterised in that comprising with
Lower step:
Step one:Burner hearth is divided into main burning area I, burning-out zone II and heat transfer zone III according to ignition quality, wherein main combustion
Burn area I and be divided into x sections according to the burner number of plies, burning-out zone II is divided into y sections according to the burnout degree number of plies, and heat transfer zone III is not segmented, to stove
Thorax carries out relevant data acquisition, main collection boiler real time execution parameter, as-fired coal prime number evidence and boiler furnace structure and design
Parameter;
Wherein, x refers to the burner number of plies, and y refers to the burnout degree number of plies;
Step 2:Assuming that the exit gas temperature of each section in main burning area I, according to each section Calculation of Heat Transfer mould established
Type, calculates the respective water-cooling wall thermal effective coefficient of x section of primary combustion zone;
Step 3:Assuming that the exit gas temperature of each section of burning-out zone II, according to each section Calculation of Heat Transfer mould established
Type, calculates the respective water-cooling wall thermal effective coefficient of y section of burning-out zone;
Step 4:According to the Calculation of Heat Transfer model of heat transfer zone III, the water-cooling wall thermal effective coefficient ψ of heat exchange area is calculatedⅢ;
Step 5:According to the water-cooling wall thermal effective coefficient of the multiple sections calculated, with reference to radiation and conduction heat transfer model
Calculate the slagging thermal resistance of multiple sections;
Step 6:The outlet cigarette temperature for each section assumed in previous steps is checked, step 7 is carried out if all meeting,
If incongruent repeat step two arrives step 6;
Step 7:The real-time exit gas temperature of the multiple sections calculated is exported as the Temperature Distribution in burner hearth, it is defeated
The real-time water-cooling wall thermal effective coefficient for going out the multiple sections calculated is used as the real-time slagging visual data of each section.
Further, in the step one, collection boiler real time execution parameter includes boiler fired coal amount, furnace outlet oxygen
Amount, First air account for total blast volume ratio, Secondary Air account for total blast volume ratio, First air import and export wind-warm syndrome, Secondary Air import and export wind-warm syndrome, it is each
Water-cooling wall wall temperature, the flue gas temperature of hearth outlet of section, real time data is gathered by Power Plant DCS System.
Further, in the step one, as-fired coal prime number divides according to elementary analysis, Industrial Analysis and the calorific value for including coal
Analysis, also needs to the proportioning of different coal samples if burnt coal sample is blending coal, and as-fired coal prime number evidence is obtained by coal analysis.
Further, in the step one, boiler furnace structure and design parameter include the overall heat transfer area of burner hearth, difference
The heat transfer area of section, dischargeable capacity, computed altitude, row's burner arrangement difference in height, burner are averagely arranged height, gone out up and down
Mouth smokestack area, burner hearth air leakage coefficient, the air leakage coefficient of pulverized coal preparation system, chamber structure and design parameter are used and set by boiler
Specification is counted to obtain.
Further, the step 2 detailed process is,
2.1 assume the exit gas temperature T of the section of primary zone the 1stⅠ1", calculate floor burner (the Ith area the 1st of primary zone the 1st
Section) water-cooling wall thermal effective coefficient ψⅠ1:
According to the equation of heat balance of the section of primary zone the 1st
Calculate water-cooling wall thermal effective coefficient ψⅠ1;
2.2 assume the exit gas temperature T of the section of primary zone i-thⅠi", calculate primary zone the i-th floor burner (the Ith area i-th
Section, 1<I≤x) water-cooling wall thermal effective coefficient ψⅠi:
According to the equation of heat balance of the section of primary zone i-th:
Calculate water-cooling wall thermal effective coefficient ψⅠi, herein
Wherein, the i of subscript I represents i-th section in primary zone I, a certain section of Ith area that i representation modules are currently calculated, and i-1 is represented
The previous section of the current calculation of sector of module, 1<Because the heat transfer model of the 1st section is different therefore individually arranges in i≤x, this step
Go out, symbolic interpretation is identical with the i-th section;The n occurred in formula is used as referring to function in algebraically sum formula, without actual meaning
Justice;QkFor the heat for the air (containing leaking out) brought into unit mass fuel in stove, kJ/kg.QrStove is brought into for unit quality fuels
Interior heat, generally equal to fuel net calorific value as received basis, kJ/kg;Q6For other overall heat loss of burner hearth, kJ/kg can root
Chosen according to boiler design book design load;σ0For Boltzmann constant, 5.67 × 10 are generally taken-11kW/(m2·K4);For burner hearth
Overall blackness, is calculated by fuel and obtained, and this is thermodynamic computing general knowledge, is repeated no more;BjiFor Ith area, i-th section calculates quantity combusted,
Kg/s, it is believed thatTⅠi" it is i-th section of Ith area exit gas temperature, K;IⅠi" it is i-th section of Ith area exiting flue gas enthalpy, kJ/
Kg, according to TⅠi" look into and take flue gas enthalpy temperature table to obtain;TⅠiFor Ith area, i-th section of flue gas mean temperature, K;βcrFor the burn-off rate of fuel,
Boiler handbook can be consulted;ψ " be lower curtate to the radiation thermal effective coefficient of upper curtate, typically take 0.1;FⅠiFor the outlet of i-th section of Ith area
Furnace cross-sectional area, m2;HⅠiFor Ith area, i-th section of water-cooling wall heat transfer area, m2;ψⅠiFor i-th section of Ith area water-cooling wall thermal effective coefficient.
Further, the step 3 detailed process is,
3.1 assume the exit gas temperature T of the section of burning-out zone the 1stⅡ1", calculate burning-out zone the 1st floor burnout degree (the IIth area the
1 section) water-cooling wall thermal effective coefficient ψⅡ1:
According to the equation of heat balance of this section
Calculate water-cooling wall thermal effective coefficient ψⅡ1, herein
3.2 assume the exit gas temperature T of burning-out zone kth sectionⅡk", calculate burning-out zone kth floor burnout degree (the IIth area the
K sections, 1<K≤y) water-cooling wall thermal effective coefficient ψⅡk:
According to the equation of heat balance of this section
Calculate water-cooling wall thermal effective coefficient ψⅡk, herein
Wherein, the k of subscript II represents the kth section of burning-out zone II, a certain section of IIth area that k representation modules are currently calculated, k-1 tables
Show the previous section of the current calculation of sector of module, 1<Because the heat transfer model of the 1st section is different therefore independent in k≤y, this step
List, symbolic interpretation is identical with kth section;ΔβcrFor the uncombusted rate of main combustion zone fuel;TⅡk" for II area kth section outlet cigarette
Temperature degree, K;IⅡk" for II area's kth section exiting flue gas enthalpy, kJ/kg, according to TⅡk" look into and take flue gas enthalpy temperature table to obtain;TⅡkFor II
The flue gas mean temperature of area's kth section, K;ψ " be lower curtate to the radiation thermal effective coefficient of upper curtate, typically take 0.1;FⅡkFor II
Area kth section outlet furnace cross-sectional area, m2;HⅡkFor the water-cooling wall heat transfer area of II area's kth section, m2;ψⅡkFor II area's kth section water cooling
Wall thermal effective coefficient.
Further, the step 4 detailed process is,
The equation of heat balance of this section
Calculate water-cooling wall thermal effective coefficient ψⅢ, herein
Wherein, subscript III represents heat transfer zone III;Tf" it is flue gas temperature of hearth outlet, K, i.e. heat transfer zone exit gas temperature;
If" it is furnace outlet flue gas enthalpy, kJ/kg, according to Tf" look into and take flue gas enthalpy temperature table to obtain;TⅢFor the average temperature of flue gas of heat transfer zone
Degree, K;ψpFor the radiation thermal effective coefficient of double of radiation heating-surface of furnace outlet, it can be chosen by boiler design value.FⅢFor heat transfer zone
Export furnace cross-sectional area, m2;HⅢFor the water-cooling wall heat transfer area of heat transfer zone, m2;ψⅢFor heat transfer zone water-cooling wall thermal effective coefficient.
Further, the step 5 detailed process is,
5.1 according to the section coherent radiation of primary zone I the 1st and conduction heat transfer model, calculates the 1st layer of primary zone burner the (the Ith
The 1st section of area) slagging thermal resistance RⅠ1:
By equation group Solve
Slagging thermal resistance RⅠ1;
5.2 according to the section coherent radiation of primary zone I i-th and conduction heat transfer model, calculates i-th layer of primary zone burner the (the Ith
I-th section of area, 1<I≤x) slagging thermal resistance RⅠi:
By equation groupSolve slagging
Thermal resistance RⅠi;
5.3 according to the section coherent radiation of burning-out zone II the 1st and conduction heat transfer model, calculates the 1st layer of burnout degree of burning-out zone (the
The 1st section of IIth area) slagging thermal resistance RⅡ1:
By equation groupSolve slagging
Thermal resistance RⅡ1;
5.4 according to the kth section coherent radiation of burning-out zone II and conduction heat transfer model, calculates burning-out zone kth layer burnout degree (the
II area's kth section, 1<K≤y) slagging thermal resistance RⅡk:
By equation groupSolve knot
Slag thermal resistance RⅡk;
5.5, according to the coherent radiation of heat transfer zone III and conduction heat transfer model, calculate the slagging thermal resistance R in this regionⅢ:
By equation groupSolve slagging
Thermal resistance RⅢ;
Wherein, εzFor the blackness on slagging surface, 0.8~0.9 can use;Tz-Ⅰi、Tz-Ⅱk、Tz-ⅢRespective segments are represented respectively
Slagging surface temperature, K;Tb-Ⅰi、Tb-Ⅱk、Tb-ⅢThe water cooling wall surface temperature of respective segments, K are represented respectively;RⅠi、RⅡk、RⅢRespectively
Represent the slagging thermal resistance of respective segments, m2·K/kW。
Further, the step 6 detailed process is,
6.1 are checked according to check formula:
|RⅠ1(1-ψⅠ2/XⅠ2)-RⅠ2(1-ψⅠ1/XⅠ1)|/RⅠ1(1-ψⅠ2/XⅠ2)≤0.5%, | RⅠ(i-1)(1-ψⅠi/XⅠi)-RⅠi
(1-ψⅠ(i-1)/XⅠ(i-1))|/RⅠ(i-1)(1-ψⅠi/XⅠi)≤0.5% (1<i≤x)、|RⅠx(1-ψⅡ1/XⅡ1)-RⅡ1(1-ψⅠx/XⅠx)|/
RⅠx(1-ψⅡ1/XⅡ1)≤0.5%, | RⅡ(k-1)(1-ψⅡk/XⅡk)-RⅡk(1-ψⅡ(k-1)/XⅡ(k-1))|/RⅡ(k-1)(1-ψⅡk/XⅡk)≤
0.5% (1<k≤y)、|RⅡy(1-ψⅢ/XⅢ)-RⅢ(1-ψⅡy/XⅡy)|/RⅡy(1-ψⅢ/XⅢThe common x+y check formula pair of)≤0.5%
The outlet cigarette temperature T for each section assumed in previous stepsⅠ1"~TⅠx″、TⅡ1"~TⅡy" common x+y cigarette temperature is checked;
6.2 carry out step 7 if all meeting, if there is incongruent repeat step two to arrive step 6;
Wherein, XⅠi、XⅡk、XⅢThe Angle Factor of Waterwall of respective segments is represented respectively, can be according to Structure Calculation.
The present invention compared with prior art, with advantages below and effect:
1st, the invention provides Temperature Distribution flexible measurement method in a kind of burner hearth based on burner hearth Multi sectional slagging situation,
Do not increase on the basis of station boiler measuring point, burner hearth is divided into multiple sections according to ignition quality, passes through the burner hearth established
The section of burner hearth cigarette temperature that inside smoke soft-sensing model carries out different sections is calculated, and provides different sections inside burner hearth simultaneously
Real-time slagging situation;
2nd, the present invention can provide flue-gas temperature distribution and the slagging situation of multiple sections in burner hearth, be hearth combustion adjustment
Data reference is provided with optimization;
3rd, the present invention can be applied to the station boiler of various structures type.
Embodiment
Below in conjunction with the accompanying drawings and the present invention is described in further detail by embodiment, following examples are to this hair
Bright explanation and the invention is not limited in following examples.
The boiler that the embodiment of the present invention is chosen is certain 600MW supercritical once-through boiler, boiler model HG-1956/25.4-
YM5 types, are resuperheat, a direct current cooker for the built-in recirculation pump activation system of supercritical pressure variable-pressure operation band.This
Boiler arranges that single burner hearth, balanced draft, dry ash extraction, turbulent burner are using front-back wall arrangement, opposed firing using Π types.
Boiler front-back wall respectively arranges 3 layers of turbulent burner (LNASB), and above the superiors' coal burner, front-back wall respectively arranges 1 layer of combustion
Most air port.
As shown in Fig. 2 Temperature Distribution flexible measurement method in a kind of burner hearth based on burner hearth Multi sectional slagging situation, comprising with
Lower step:
Step one:Burner hearth is divided into main burning area I, burning-out zone II and heat transfer zone III according to ignition quality, wherein main combustion
Burn area I and be divided into 3 sections according to the burner number of plies, burning-out zone II is divided into 1 section according to the burnout degree number of plies, and heat transfer zone III is typically no longer divided
Section.Relevant data acquisition, main collection boiler real time execution parameter, as-fired coal prime number evidence and boiler furnace structure are carried out to burner hearth
And design parameter.Wherein, boiler real time execution parameter accounts for total blast volume ratio including boiler fired coal amount, oxygen at furnace exit, First air
Example, Secondary Air account for total blast volume ratio, First air and import and export wind-warm syndrome, Secondary Air import and export wind-warm syndrome, water-cooling wall wall temperature, the stove of each section
Thorax exit gas temperature (if can be calculated without measuring point along inverse flue gas flow) etc., can gather real time data by Power Plant DCS System;Enter
Stove coal data is obtained by coal analysis, main elementary analysis, Industrial Analysis and calorimetry including coal etc., such as burning coal
Sample then also needs to the proportioning of different coal samples for blending coal;Chamber structure and design parameter can be used and specification by boiler
Book is obtained, it is necessary to burner hearth entirety heat transfer area, the heat transfer area of different section, dischargeable capacity, computed altitude, arrange burner up and down
Arrangement difference in height, burner averagely arrange height, outlet smokestack area, burner hearth air leakage coefficient, the air leakage coefficient of pulverized coal preparation system.
Step 2:Assuming that the exit gas temperature of each section in main burning area I, according to each section Calculation of Heat Transfer mould established
Type, calculates the respective water-cooling wall thermal effective coefficient of 3 sections of primary combustion zone:
2.1 assume the exit gas temperature T of the section of primary zone the 1stⅠ1", calculate floor burner (the Ith area the 1st of primary zone the 1st
Section) water-cooling wall thermal effective coefficient ψⅠ1:
According to the equation of heat balance of the section of primary zone the 1st
Calculate water-cooling wall thermal effective coefficient ψⅠ1;
2.2 assume the exit gas temperature T of the section of primary zone the 2ndⅠ2", calculate floor burner (the Ith area the 2nd of primary zone the 2nd
Section) water-cooling wall thermal effective coefficient ψⅠ2:
According to the equation of heat balance of the section of primary zone the 2nd:
Calculate water-cooling wall thermal effective coefficient ψⅠ2, herein
2.3 assume the exit gas temperature T of the section of primary zone the 3rdⅠ3", calculate floor burner (the Ith area the 3rd of primary zone the 3rd
Section) water-cooling wall thermal effective coefficient ψⅠ3:
According to the equation of heat balance of the section of primary zone the 3rd:
Calculate water-cooling wall thermal effective coefficient ψⅠ3, herein
Step 3:Assuming that the exit gas temperature of burning-out zone II, according to each section Calculation of Heat Transfer model established, is calculated
The water-cooling wall thermal effective coefficient of 1 section of burning-out zone:
Assuming that the exit gas temperature T of the section of burning-out zone the 1stⅡ1", calculate floor burnout degree (the IIth area the 1st of burning-out zone the 1st
Section) water-cooling wall thermal effective coefficient ψⅡ1:
According to the equation of heat balance of this section
Calculate water-cooling wall thermal effective coefficient ψⅡ1, herein
Step 4:According to the Calculation of Heat Transfer model of heat transfer zone III, the water-cooling wall thermal effective coefficient ψ of heat exchange area is calculatedⅢ:
The equation of heat balance of this section
Calculate water-cooling wall thermal effective coefficient ψⅢ, herein
Step 5:According to the water-cooling wall thermal effective coefficient of the multiple sections calculated, with reference to radiation and conduction heat transfer model
Calculate the slagging thermal resistance of multiple sections:
5.1 according to the section coherent radiation of primary zone I the 1st and conduction heat transfer model, calculates the 1st layer of primary zone burner the (the Ith
The 1st section of area) slagging thermal resistance RⅠ1:
By equation group Solve
Slagging thermal resistance RⅠ1。
5.2 according to the section coherent radiation of primary zone I the 2nd and conduction heat transfer model, calculates the 2nd layer of primary zone burner the (the Ith
The 2nd section of area) slagging thermal resistance RⅠ2:
By equation groupSolve slagging
Thermal resistance RⅠ2。
5.3 according to the section coherent radiation of primary zone I the 3rd and conduction heat transfer model, calculates the 3rd layer of primary zone burner the (the Ith
The 3rd section of area) slagging thermal resistance RⅠ3:
By equation groupSolve slagging heat
Hinder RⅠ3。
5.4 according to the section coherent radiation of burning-out zone II the 1st and conduction heat transfer model, calculates the 1st layer of burnout degree of burning-out zone (the
The 1st section of IIth area) slagging thermal resistance RⅡ1:
By equation groupSolve slagging
Thermal resistance RⅡ1。
5.5, according to the coherent radiation of heat transfer zone III and conduction heat transfer model, calculate the slagging thermal resistance R in this regionⅢ:
By equation groupSolve slagging
Thermal resistance RⅢ。
Step 6:The outlet cigarette temperature for each section assumed in previous steps is checked, step 7 is carried out if all meeting,
If incongruent repeat step two arrives step 6:
6.1 are checked according to check formula:
|RⅠ1(1-ψⅠ2/XⅠ2)-RⅠ2(1-ψⅠ1/XⅠ1)|/RⅠ1(1-ψⅠ2/XⅠ2)≤0.5%, | RⅠ2(1-ψⅠ3/XⅠ3)-RⅠ3(1-
ψⅠ2/XⅠ2)|/RⅠ2(1-ψⅠ3/XⅠ3)≤0.5%, | RⅠ3(1-ψⅡ1/XⅡ1)-RⅡ1(1-ψⅠ3/XⅠ3)|/RⅠ3(1-ψⅡ1/XⅡ1)≤
0.5%th, | RⅡ1(1-ψⅢ/XⅢ)-RⅢ(1-ψⅡ1/XⅡ1)|/RⅡ1(1-ψⅢ/XⅢ)≤0.5% checks formula in previous steps for totally 4
The outlet cigarette temperature T for each section assumedⅠ1″、TⅠ2″、TⅠ3″、TⅡ1" totally 4 cigarette temperature are checked;
6.2 carry out step 7 if all meeting, if incongruent repeat step two arrives step 6.
Step 7:Export the real-time exit gas temperature T of the multiple sections calculatedⅠ1″、TⅠ2″、TⅠ3″、TⅡ1" it is used as stove
Section temperature distribution in thorax, exports the real-time water-cooling wall thermal effective coefficient ψ of the multiple sections calculatedⅠ1、ψⅠ2、ψⅠ3、ψⅠx、ψⅢ
As the real-time slagging visual data of each section, (it shows that more greatly the stronger slagging of water wall absorption radianting capacity is less, and its is smaller
Show that the poorer slagging of water wall absorption radianting capacity is more serious), operations staff is presented to, the ginseng for carrying out boiler operatiopn optimization is used as
Examine.
Above content described in this specification is only illustration made for the present invention.Technology belonging to of the invention
The technical staff in field can be made various modifications or supplement to described specific embodiment or be substituted using similar mode, only
Will without departing from description of the invention content or surmount scope defined in the claims, all should belong to the present invention guarantor
Protect scope.