CN106524118B - The method for building up of boiler wear resistant explosion-proof temperature field model - Google Patents
The method for building up of boiler wear resistant explosion-proof temperature field model Download PDFInfo
- Publication number
- CN106524118B CN106524118B CN201610866742.1A CN201610866742A CN106524118B CN 106524118 B CN106524118 B CN 106524118B CN 201610866742 A CN201610866742 A CN 201610866742A CN 106524118 B CN106524118 B CN 106524118B
- Authority
- CN
- China
- Prior art keywords
- temperature
- point
- steam
- boiler
- pipe
- 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
Classifications
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F22—STEAM GENERATION
- F22B—METHODS OF STEAM GENERATION; STEAM BOILERS
- F22B35/00—Control systems for steam boilers
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F22—STEAM GENERATION
- F22B—METHODS OF STEAM GENERATION; STEAM BOILERS
- F22B37/00—Component parts or details of steam boilers
- F22B37/02—Component parts or details of steam boilers applicable to more than one kind or type of steam boiler
- F22B37/42—Applications, arrangements, or dispositions of alarm or automatic safety devices
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F22—STEAM GENERATION
- F22B—METHODS OF STEAM GENERATION; STEAM BOILERS
- F22B37/00—Component parts or details of steam boilers
- F22B37/02—Component parts or details of steam boilers applicable to more than one kind or type of steam boiler
- F22B37/42—Applications, arrangements, or dispositions of alarm or automatic safety devices
- F22B37/421—Arrangements for detecting leaks
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F23—COMBUSTION APPARATUS; COMBUSTION PROCESSES
- F23M—CASINGS, LININGS, WALLS OR DOORS SPECIALLY ADAPTED FOR COMBUSTION CHAMBERS, e.g. FIREBRIDGES; DEVICES FOR DEFLECTING AIR, FLAMES OR COMBUSTION PRODUCTS IN COMBUSTION CHAMBERS; SAFETY ARRANGEMENTS SPECIALLY ADAPTED FOR COMBUSTION APPARATUS; DETAILS OF COMBUSTION CHAMBERS, NOT OTHERWISE PROVIDED FOR
- F23M11/00—Safety arrangements
Landscapes
- Engineering & Computer Science (AREA)
- Mechanical Engineering (AREA)
- General Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Thermal Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Combustion & Propulsion (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Control Of Steam Boilers And Waste-Gas Boilers (AREA)
Abstract
The invention discloses a kind of method for building up of boiler wear resistant explosion-proof temperature field simulation model, profiling temperatures in boiler are simulated in real time according to boiler load, flue-gas temperature, vapor (steam) temperature, the outlet practical measuring point of wall temperature, it is tested by the compacted swollen data that all previous maintenance measures, pass through the invalid coefficient of each factor of neural network fashion self study high-temperature material leakage failure, simultaneously according to same unit capacity, with boiler manufacturer, with the data of the same trade or business of the boiler of combustion system, intersect for lateral model and learn, forms temperature field in furnaces simulation model.The present invention forms the Temperature Distribution simulation model of the explosion-proof self study self refresh of boiler wear resistant on the basis of big data, can 360 degree rotation inquiry, the spatial position in overtemperature region in boiler is positioned in real time, calculate the compacted swollen trend of boiler heating surface and Residual Life, emphasis checkpoint guide service is provided to formulate repair schedule, effectively prevent boiler since overtemperature causes boiler venting of dust explosion, improves economy and the safety of boiler operatiopn.
Description
Technical field
The present invention relates to the method for building up of boiler wear resistant explosion-proof temperature field simulation model, belong to boiler wear resistant explosion-proof technology neck
Domain.
Background technology
Boiler is one of important equipment of heat power plant boiler, superheater, reheater, economizer, four class pipe of water-cooling wall
(referred to as " four pipes ") leakage happens occasionally with demolition accident, and boiler only takes forced outage to stop once " four pipes " leakage occurs
The method of stove is repaired, and has been seriously affected the normal production in thermal power plant, has been caused huge economic loss.China's firepower hair
Boiler breakdowns account for about 70% or more of the total accident in power plant in the various accidents of power plant, and " four pipes " quick-fried leakage becomes long-standing problem firepower electricity
The a great problem of factory's safety in production.According to statistics：Boiler is non-every time stop caused by direct economic loss about 2,500,000 or so, and it is general
Unit that electricity power enterprise causes due to pipe explosion accident every year on average is non-to stop 2-3 times, and enterprise direct losses 500-750's ten thousand is direct
Economic benefit, while the casualties that safety accident may be caused to cause.
On the one hand prevention for " four pipes " quick-fried leakage problem is to reinforce the specific aim maintenance of equipment and safeguard, is on the one hand
Reinforce equipment running monitoring and adjustment in time.Apparatus monitoring method is various, generally uses following methods：
(1) according to the compacted swollen record of all previous maintenance of power plant, compacted swollen rate is calculated according to time factor；
Such as：The compacted swollen value 51.4mm of part is checked in 5 monthly tests in 2013, and in August, 2014 checks the compacted swollen value of same position pipe fitting
52.3mm；Predicting in July, 2015, compacted swollen value is 53.02mm.
Calculating process：(52.3-51.4)/15 month=X/12 months；X=0.72mm
In June, 2015, compacted swollen value was：52.3+0.72=53.02mm
Wherein -2014 years in Mays, 2,013 7 monthly mortgage 15 months；- 2015 years 12 months July of in August, 2014.
Such product disadvantage：Compacted swollen prediction accuracy is very low, can make qualitative reference in the case of can only being directed to variation less,
It is quantitative very inaccurate
Boiler during normal operation, required to follow load curve lifting load according to dispatching of power netwoks, different load is coal-fired
Amount, coal-fired mixture ratio, air output, feedwater flow is different, therefore causes exhaust gas volumn, flue gas flow rate, each section of steam flow, respectively
Heating surface exit wall temperature value follows variation, such such multivariable of product not to consider that accuracy is very directly according to time prediction
It is low.
There are environment for such product：Change less for a long time for unit load, unit coal is single, and combustion adjustment is stablized
Under the premise of can do qualitative reference, qualitative reference it is not recommended that.
(2) limited wall temperature measurement point is exported according to heating surface, observes each measuring point overtemperature record, qualitative analysis overtemperature situation.
Such as:Continue overtemperature 20 minutes or more, checks that overtemperature records, pay close attention to；Overtemperature duration exceeds 10 times in period,
It checks that overtemperature records, pays close attention to.
Such product disadvantage：
Wall temperature measurement point distribution is limited, simultaneously because wall temperature measurement point is mounted on, big packet is outer, and the data of measurement are less than the big packet of boiler
Interior actual temperature value, there are practical overtemperature, and measuring point not overtemperature situation, increase overtemperature hidden danger.
Such as certain factory's two-stage superheater point layout, 96 pipe fittings of comb part 1152, only 20 outlet measuring points.By power plant's reality
Between border heating surface has 7-12, it is more than 10000 pipe fittings of mean value computation, measuring point 200-500 or so.Measuring point is well below reality
Demand has 95% or more to lack temperature monitoring, and there are overtemperature hidden danger.
For another example point position is outer in big packet, is not inconsistent with practical pipe fitting position, cannot express pipe fitting each position temperature in burner hearth
Distribution situation is practical to be far below actual demand value to temperature monitoring point.
Different concern standards can be formulated for overtemperature experience difference according to different people, such as：Based on lasting overtemperature 20 minutes with
On, overtemperature duration exceeds 10 times in the period, can only be with formulation is felt, due to not having analysis model, practical operation situation to exist again
Constantly variation cannot effectively analyze situation of change and be executed with feeling.Therefore concern accuracy in overtemperature region is low, and promptness is poor.
There are environment for such product：The special tool of boiler has to this boiler of our factory 5 years or more practical experience, can do qualitative ginseng
It examines, does not accomplish qualitative reference.
(3) traditional material lifetime estimation method：Creep rupture strength method, La-rice (L-M) equation, Robinson life consumptions
Method.
Creep rupture strength method：
Creep rupture strength method calculates remaining life and is built upon on the basis of the creep rupture strength curve on log-log coordinate at present,
At a certain temperature, the stress of creep rupture strength curve is expressed with rupture time relationship by formula, i.e. t=A σ-B
σ in formula -- stress, MPa
T-rupture time, h；
A, B-material constant.
By t=A σ-BBoth sides take logarithm, then can be obtained
Lgt=lgA-Blg σ, thus formula is it is found that stress and the creep rupture strength curve of time are linear relations.Currently, according to
Creep rupture strength curve disperses lower limit and middle limit, and two kinds of shapes of minimum theoretical service life and average theory service life are used in the service life calculates
Formula, i.e.,
(1) the minimum theoretical service life
(2) the average theory service life
T in formulamin-- minimum theoretical service life, h；
T---- average theory service life, h；
σzs-- working stress (pipeline internal pressure reduced stress)
-- 10 under certain temperature4H creep rupture strength temperature values, Mpa；
-- 10 under certain temperature5H creep rupture strength temperature values, Mpa；
-- 2*10 under certain temperature5H creep rupture strength temperature values, Mpa；
The method that industry provides two sets of calculating remaining life, is selected according to different situations, even if under same stress,
The result calculated with average theory Life method using minimum theoretical Life method is differed greatly.Someone is multiplied by the value of coefficient n, n value
It varies with each individual, mixes many artificial factors, the service life, all coefficients lacked data based on experience using theoretical value as foundation
Support.
Often increase 9.8Mpa according to working stress, remaining life reduces an order of magnitude, and remaining life is to working stress wave
Dynamic sensibility is especially big, so slightly deviation remaining life result is mutually far short of what is expected.And the method mainly with people's empirical coefficient be according to
According to the service life accuracy of prediction is very low.
When calculating high-temperature component remaining life using creep rupture strength formula it may first have to determinePersistently
Intensity value.It is the creep rupture strength curve linear extrapolation institute by several experimental points of thousands of hours in short-term currently, determining creep rupture strength value
, but practical proof, curve linear extrapolation method have a long way to go.
La-rice (L-M) equation：
Time-temperature parametric method, basic conception are in creep rupture test, when experiment can be shortened by improving experimental temperature
Between.Attempt by the short-term experiment data under higher temperature come long term data when lower temperature under certain stress of extrapolating.
V=Ae-Q/RT
V-steady creep speed, mm/h；
A-constant；
R-gas constant；
T- absolute temperature, K；
Q-activation energy, 4.1816J/mol
It is assumed that creep fracture time trIt is inversely proportional with secondary creep rates v, tr=AeQ/RT, wherein Q is stress function, can be pushed away
To going out：P (σ)=T (C+lgtr)
Have too many in rummy equation it is assumed that this hypothesis lacks lot of experimental data, also without sufficient creep rupture machine
The physical theory foundation of system, proves according to experimental data, material structure variation difference and creep at a temperature of different creep tests
Fracture mechanism is different, therefore temperature extrapolation method is infeasible.
Robinson life consumption methods：
Assuming that：In creep process, when material given stress and at a temperature of after the t times run, life consumption
Score is t/tr,trIt is the creep fracture time under same stress and temperature condition.Under conditions of stress or temperature change,
Robinson assumes that every part of service life trumpet score is independent for other each parts, and assume when its life consumption score it
It is with for 1, with regard to destroying, i.e. ∑ t/tr=1t certain stress and at a temperature of run time, h；trIn certain stress and temperature
Lower creep fracture time, h；The number of n temperature or stress variation.
In Robinson rules, trValue determination be by original material under certain temperature and stress creep rupture strength experiment come
Estimation.Its correctness depends on creep rupture strength curve linear extrapolation and La-rice Time-temperature parameter is external, so accurate
Exactness is relatively low, unreliable.
Invention content
The present invention provides a kind of boiler wear resistant explosion-proof temperature field stimulation method, according to load, flue-gas temperature, vapor (steam) temperature,
The outlet practical measuring point of wall temperature simulates the overtemperature distribution situation of each heating surface, pipe, inflection point and pipeline section in real time, on the basis of big data
Formed the explosion-proof self study self refresh of boiler wear resistant simulation model, can 360 degree rotation inquiry, in real time position boiler in overtemperature
The spatial position in region calculates the compacted swollen trend of boiler heating surface and Residual Life.
The technical solution used in the present invention is：
A kind of method for building up of boiler wear resistant explosion-proof temperature field simulation model, according to load, flue-gas temperature, vapor (steam) temperature,
The outlet practical measuring point of wall temperature simulates the profiling temperatures of each heating surface, pipe, inflection point and pipeline section in real time, is surveyed by all previous maintenance
The compacted swollen data of amount are tested, and the failure system of each factor of neural network fashion self study high-temperature material leakage failure is passed through
Number, while according to same unit capacity, intersecting for lateral model with the data of the same trade or business of the boiler of combustion system with boiler manufacturer
Study forms and improves Temperature Field Simulation model certainly；Then according to real-time load, flue-gas temperature, vapor (steam) temperature, outlet wall temperature measurement point
Data value binding analysis model, carry out the three-dimensional visualization early warning and analysis in overtemperature region.
The method for building up of boiler wear resistant explosion-proof temperature field simulation model, includes the following steps：
(1), Temperature Field Simulation model foundation：It establishes smoke temperature field, establish vapor (steam) temperature field, computer tube tube internal heat exchange coefficient, meter
Calculate heat transfer outside a tube coefficient, simulation wall temperature temperature field；
(2), models for temperature field is verified：The verification of history big data, maintenance assessment verification.
Smoke temperature field is established to include the following steps：
(1) boiler model coordinate system is determined:Using burner hearth center line and the intersection point of 0 meter of elevation plane of boiler as origin, vertically
Direction is y-axis, is x-axis with screen vertical direction, z-axis is determined according to x-axis and y-axis；
(2) smoke temperature measuring point is chosen:The flue-gas temperature measuring point that each heating surface closes on is chosen, determines measuring point probe co-ordinate, measuring point
Quantity is no less than three；
(3) according to combustion field Temperature Distribution feature, unknown flue-gas temperature point is simulated in conjunction with known smoke temperature measuring point.
Preferably, when the known smoke temperature measuring point of selection be four when, by flue gas direction by the first two point be respectively labeled as A,
B, latter two point are respectively labeled as C, D, consider that accuracy in computation and temperature spot position, smoke temperature field form 4 kinds of calculating moulds
Formula：
Pattern 1 (such as Fig. 1)：The point centered on A points is simulated 1 point in the center line of A and B points, is distinguished in the outside of A, B
Simulate 1 point；1 point is simulated to the directions C, vertically respectively simulates 1 point；Matched curve is formed in X, Y, Z axis respectively；
Pattern 2 (such as Fig. 2)：The point centered on B points simulates 1 point, the difference mould on the outside of A, B in the center line of A and B points
Intend 1 point；1 point is simulated to the directions D, vertically respectively simulates 1 point；Matched curve is formed in X, Y, Z axis respectively；
Pattern 3 (such as Fig. 3)：The point centered on C points simulates 1 point, the difference mould on the outside of C, D in the center line of C and D points
Intend 1 point；1 point is simulated to the directions A, vertically respectively simulates 1 point；Matched curve is formed in X, Y, Z axis respectively；
Pattern 4 (such as Fig. 4)：The point centered on D points simulates 1 point, the difference mould on the outside of C, D in the center line of C and D points
Intend 1 point；1 point is simulated to the directions B, vertically respectively simulates 1 point；Matched curve is formed in X, Y, Z axis respectively；
According to combustion field Temperature Distribution feature (burner hearth central temperature is high, and both sides temperature is low) and the design of flue-gas temperature measuring point
Position；The temperature design of simulation point between two measuring point of same type is that the temperature averages of practical measuring point are multiplied by 1.01-1.05
Coefficient in range, both sides simulation point are multiplied by coefficient within the scope of 0.95-0.99 according to corresponding measuring point temperature point value；Specific choice
Coefficient is related to point layout position according to boiler type, is verified by results of calculation and measured value when maintenance.
9 palace grid space regions are formed by four points, computation model is selected according to actual prediction point region；Region
The case where beyond a computation model, calculates future position temperature value according to different models respectively, is then averaged as pre-
The temperature value of measuring point.As shown in figure 5, selecting computation model according to actual prediction point position region, being shown in Table 1, (* is labeled as choosing
The calculating pattern taken)；The case where one computation model is exceeded for region, such as：Point is calculated to fall in the region of Fig. 25, by
Table 1 predicts point value=(M1+M2+M3+M4)/4 it is found that 4 patterns should be calculated；M1, M2, M3, M4 be respectively according to pattern 1,
4 calculated future position temperature value of pattern 2, pattern 3 and pattern.
1 actual prediction point position region computation model of table
The step of establishing vapor (steam) temperature field be：
It is the steaming of influent header that each heating surface, which has the vapor (steam) temperature measuring point of influent header, every pipe fitting inlet steam temperature,
Stripping temperature；Due to enter each pipe fitting steam flow rate variation range it is small, be considered as it is identical, according to outlet wall temperature measurement point temperature value and
The temperature difference of steam, calculates the outlet steam temperature value of every pipe fitting in pipe；Then according to every pipe fitting inlet steam temperature,
The length of outlet steam temperature and pipe fitting calculates the vapor (steam) temperature each calculated on pipe fitting a little：
(1) each pipe fitting and influent header are established, the relationship of outlet header configures influent header vapor (steam) temperature measuring point；
(2) pipe fitting outlet steam temperature is calculated；
1) selection outlet wall temperature measurement point, it is vertical that measuring point distribution is selected as a cross one；
2) unknown exit wall warm spot is calculated according to known measuring point；
Exit wall temperature value=Wb/Wc × Wd；Wb, Wc are respectively the temperature value of longitudinal known temperature measuring point；Wd is laterally
The temperature value of known temperature measuring point；
3) the vapor (steam) temperature value put on pipe fitting is calculated；
Low fast, the highly endothermic slow feature of outlet steam temperature of absorbing heat of consideration inlet steam temperature, formation temperature and distance
Matched curve；
In T=T goes out-(T go out-T into) × coefficient, according to fair current 0.3, the coefficient of adverse current 0.5 calculates the temperature value of central point；
T goes out：Outlet steam temperature value；T into：Inlet steam temperature value；In T：Central point vapor (steam) temperature value；According to inlet temperature, center
It puts temperature and outlet temperature and generates matched curve along caliber direction distance；Distance according to actual point along caliber direction is brought into
Matched curve generates the vapor (steam) temperature value of this point；
(3) three dimensional space coordinate is calculated:
Each pipe fitting entrance coordinate is configured, according to entrance coordinate, pipe fitting pipe range, corner position, is converted into each calculating point
Spatial position.
The step of computer tube tube internal heat exchange coefficient is：
Computer tube tube internal heat exchange coefficient formula：
Since each pipeline section coefficient of heat transfer gap of heating surface is smaller, using average intraductal heat exchange coefficient as every on pipeline section
Coefficient of heat transfer value；
Wherein, λ is thermal coefficient, dnFor pipe fitting inside radius, ωqIt is steam flow rate, Pr Pu Lang for kinematic coefficient of viscosity, V
Special number；
Heating surface temperature averages=(inlet steam temperature+outlet steam temperature)/2
Heating surface pressure mean values=(inlet steam pressure+outlet vapor pressure)/2
According to average vapor pressure and temperature computation λ, ωq,V,Pr；
The cross-sectional flow area of steam flow rate V=evaporation capacity (t/h) × 1000 × incidence coefficient/averag density/3600/
Evaporation capacity：Measuring point is evaporated by drum side to measure；
Incidence coefficient：Due to evaporation capacity drum measuring point, there are desuperheating water penetratings, total quantity of steam to be higher than before entering superheater
Drum measuring point；Incidence coefficient=(evaporation capacity+attemperation water flow)/evaporation capacity；
Averag density：Saturated vapor calculation formula is brought into according to pressure and temp and calculates specific volume, and specific volume inverse is density；
Cross-sectional flow area=all pipe fittings of this heating surfaceEvery cross-sectional area is π R × R, and R is pipe fitting internal diameter,
XiThus on heating surface single tubular cross-sectional area.
λ is thermal coefficient：According to temperature, pressure, in coefficient of heat transfer table, computation of table lookup thermal coefficient λ.
ωqFor kinematic coefficient of viscosity：According to temperature, pressure, in kinematic coefficient of viscosity table, computation of table lookup ωq。
Pr is Prandtl number:According to temperature, pressure, in Prandtl coefficient table, computation of table lookup Pr.
Coefficient of heat transfer step is outside computer tube：
Heat transfer outside a tube coefficient a1=1/ (R1 × 3.14 × pipe outside diameter × tubular length)；
Wherein, the outer convective heat transfer resistance of R1=pipes external radiation heat exchanged thermoresistance+pipe,
Tubular length=entire heating surface total pipe range
According to heat transfer principle：
Import and export enthalpy difference × evaporation capacity × incidence coefficient=heat flow × tubular length 1.
Heat flow=(average flue-gas temperature-average steam temperature)/thermal resistance is 2.
Dirt heat in the outer convective heat transfer resistance+dust stratification thermal resistance+wall resistance+pipe of thermal resistance=pipe external radiation heat exchanged thermoresistance+pipe
Convective heat transfer resistance is 3. in resistance+pipe
According to 1., 2., R1 is 3. calculated：
(average flue-gas temperature-average steam temperature)/(inlet and outlet enthalpy difference * evaporation capacity * incidence coefficients/tubular length)-
(convective heat transfer resistance in dirtiness resistance+pipe in dust stratification thermal resistance+wall resistance+pipe)=pipe external radiation heat exchanged thermoresistance+heat convection
Thermal resistance=R1
Average flue-gas temperature=heating surface flue-gas temperature central point calculates flue-gas temperature
Average steam temperature=(inlet header box temperature degree mean value+outlet header temperature mean value)/2
Import and export enthalpy difference=H (outlet pressure measuring point, temperature point)-H (inlet pressure measuring point, temperature point)
Evaporation capacity：Measuring point is evaporated by drum side to measure；
Incidence coefficient：Due to evaporation capacity drum measuring point, there are desuperheating water penetratings, total quantity of steam to be higher than before entering superheater
Drum measuring point, incidence coefficient=(evaporation capacity+attemperation water flow)/evaporation capacity；
Tubular length=entire heating surface total pipe range
Dust stratification thermal resistance value, the reference value of common fuel dust stratification thermal resistance is selected to be shown in Table 2 according to different fuel, boiler generally selects
With coal-fired flue-gas dust stratification thermal resistivity=0.00172m2℃/w。
The reference value of 2 dust stratification thermal resistance of table
Dust stratification thermal resistance (DEG C/w)=dust stratification thermal resistivity/π d, d=bores；
Dirtiness resistance (DEG C/w)=fouling resistance coefficient/π d, d=bores.
Tube wall heat conduction coefficient determines according to tubing material and temperature, such as when tubing is 12CrIMoV, and tube wall is led at 500 DEG C
Hot coefficient=33w/ (m DEG C)
3 tube wall heat conduction coefficient of table
Thermal coefficient λ：Under the conditions of steady heat transfer, two 1m apart, area are taken perpendicular to heat conduction direction in interior of articles
For 1 square metre of parallel surface, and 1 degree of the temperature difference of the two panels, then it is flat from a planar conductive to another in 1 second
The heat in face is just defined as the thermal conductivity of the substance, unit：Watt/( meter Du)
Thermal resistance θ：Heat flow is when passing through object, in the temperature difference that object both ends are formed, i.e.,
θ=ln (outer diameter/internal diameter)/(2* π λ * L) λ：Thermal coefficient；L length of pipe sections
Wall resistance=ln (outer diameter/internal diameter)/(2*3.14*L* thermal coefficients)
Convective heat transfer resistance=1/ (π * bore * length * intraductal heat exchange coefficients in pipe)
Simulation wall temperature temperature field step is：
Outside wall temperature=T flue gases-heat flow × (dust stratification thermal resistance+heat transfer outside a tube thermal resistance) or
Outside wall temperature=T steam+heat flow × (dirtiness resistance+intraductal heat exchange thermal resistance)
Heat flow=(T flue gas-T steam)/thermal resistance
Dirt heat in the outer convective heat transfer resistance+dust stratification thermal resistance+wall resistance+pipe of thermal resistance=pipe external radiation heat exchanged thermoresistance+pipe
Convective heat transfer resistance in resistance+pipe.
History big data is verified as：According to power plant's history real time data flue-gas temperature/vapor (steam) temperature/exit wall warm spot, with
And inspection record overhauls in power plant, forms overtemperature factor and compacted swollen matched curve, according to the record sample of maintenance next time, verification is super
Temperature model correctness.
Maintenance assessment is verified as：Operation each pipeline section overtemperature data of boiler, Awaiting Overhaul boiler shutdown are calculated according to overtemperature model
Afterwards, it is verified by macro -graph and compacted swollen inspection record in the emphasis overtemperature region calculated according to overtemperature model；With number
Increase according to amount, sample size increases, automatic to correct overtemperature model, forms self learning model, a newest record of examination is as maintenance
The foundation run afterwards executes；Simultaneously according to same unit capacity, used with the data of the same trade or business of the boiler of combustion system with boiler manufacturer
Intersect in lateral model and learn, forms temperature field in furnaces simulation model.
Temperature field is reminded applied to power plant operations staff, reminds overtemperature key area in real time, is found security risk, is passed through wind
Coal matches, and the adjustment means such as desuperheating water adjust in time, prevents boiler heating surface from continuing overtemperature, improves the boiler heating surface remaining longevity
Life.
Temperature field is applied to power plant's repair schedule and formulates the prompting for checking position, for there are the overtemperature areas of security risk
Domain reminds service personnel to specify corresponding repair schedule, prevents impaired heating surface maintenance from omitting.
For the present invention according to real-time flue-gas temperature/vapor (steam) temperature/exit wall warm spot, combination temperature field computation model is per minute
Thermo parameters method situation is calculated, is reminded in real time for overtemperature region, illustrates that overtemperature region is most worth situation with overtemperature.
Boiler wear resistant explosion-proof temperature of the present invention field analogue simulation is capable of providing 360 degree rotation inquiry, can inquire every pipe
Part, due to pipe fitting 1：1 high-fidelity is drawn, and spatial position in boiler can intuitively be positioned by positioning pipe fitting in graphics；Overtemperature is provided
Basic operating condition, main vapour pressure, Stream temperature degree, main steam flow；There is provided overtemperature areas case, each overtemperature region it is specific
Spatial position, overtemperature are most worth；The spatial position that overtemperature is most worth is provided.
The present invention can also provide historical temperature field:
Historical temperature field provides profiling temperatures of each heating surface within the period.Time shaft can be dragged, checks temperature
Spend distribution situation, can independently select heating surface according to analysis demand, can also global analysis, while when current selection is provided
The overtemperature situation for carving each heating surface provides each heating surface X/Y/Z axis depth profiles.
It is superimposed temperature field:
According to each heating surface situation, overlay analysis heating surface overtemperature.Intuitively show overtemperature distribution situation by three-dimensional, clicks
Each specific each point that calculates of heating surface prompt is most worth temperature, mean temperature and overtemperature duration, and can check details, provides temperature
Spend historical analysis curve and historical temperature data
The present invention is tested by the compacted swollen data that all previous maintenance measures, and passes through neural network fashion self study high temperature material
The invalid coefficient of each factor of leakage failure is expected, while according to same unit capacity, with boiler manufacturer, with the boiler of combustion system
Data of the same trade or business intersect study for lateral model；Form temperature field in furnaces simulation model；Then according to real-time load, flue gas temperature
Degree, vapor (steam) temperature, the data value combination simulation model for exporting wall temperature measurement point, carry out the early warning and analysis in overtemperature region.
After being managed by boiler wear resistant explosion protection system, reduces boiler and booster number occurs, can be enterprise's band every year
Carry out the direct economic benefit of 100-200 ten thousand.Therefore boiler wear resistant it is explosion-proof be ensure unit safety, stabilization, economical operation basis.
It is using advantageous effect caused by above-mentioned technical proposal：
The present invention according to load, flue-gas temperature, vapor (steam) temperature, outlet the practical measuring point of wall temperature simulate in real time each heating surface,
The overtemperature distribution situation of pipe, inflection point and pipeline section forms the imitative of the explosion-proof self study self refresh of boiler wear resistant on the basis of big data
True mode, can 360 degree rotation inquiry, in real time position boiler in overtemperature region spatial position, calculate boiler heating surface it is compacted swollen
Trend and Residual Life provide " personal " guide service for repair schedule, effectively prevent boiler since overtemperature causes boiler
Venting of dust explosion improves economy and the safety of boiler operatiopn.
Description of the drawings
The present invention will be further described in detail below with reference to the accompanying drawings and specific embodiments.
Fig. 1 is that flue-gas temperature point calculates ideograph；
Fig. 2 is the nine grids figure that flue-gas temperature point is formed.
Specific implementation mode
The method for building up of boiler wear resistant explosion-proof temperature of the present invention field simulation model, including：Flue-gas temperature/steam the temperature at scene
Degree/outlet wall temperature measurement point is searched and configuration；The foundation of basic three-dimensional boiler model, parameter configuration；It is mended according to power plant's historical data
Temperature field overtemperature situation is calculated, determines expert model parameter；The every of pipe fitting is calculated in conjunction with expert model according to measuring point data in real time
The real-time temperature values of a point, and tubing temperature standard is exceeded according to pipe fitting material, it is recorded as overtemperature, reminds overtemperature region；Pass through
It is superimposed temperature field and carries out pipe fitting health status.Specifically include the following steps：
Step 1, Temperature Field Simulation model foundation；
1. the flue-gas temperature measuring point at scene is searched and configuration, smoke temperature field is established：
(1) boiler model coordinate system is determined:Using burner hearth center line and the intersection point of 0 meter of elevation plane of boiler as origin, vertically
It is positive direction of the y-axis upwards, screen outward direction is positive direction of the x-axis, and z-axis is determined according to x-axis and y-axis；
(2) smoke temperature measuring point is chosen:The flue-gas temperature measuring point that each heating surface closes on is chosen, determines measuring point probe co-ordinate, measuring point
Quantity is no less than three；
(3) according to combustion field Temperature Distribution feature, unknown flue-gas temperature point is simulated in conjunction with known smoke temperature measuring point.
When choosing four flue-gas temperature measuring points, using following methods：
The first two point is respectively labeled as A, B by flue gas direction, latter two point is respectively labeled as C, D, and it is accurate to consider to calculate
Degree and temperature spot position, smoke temperature field form 4 kinds of calculating patterns：
Pattern 1：The point centered on A points simulates 1 point in the center line of A and B points, 1 is simulated respectively in the outside of A, B
Point；1 point is simulated to the directions C, vertically respectively simulates 1 point；Matched curve is formed in X, Y, Z axis respectively；
Pattern 2：The point centered on B points simulates 1 point in the center line of A and B points, simulates 1 point respectively on the outside of A, B；
1 point is simulated to the directions D, vertically respectively simulates 1 point；Matched curve is formed in X, Y, Z axis respectively；
Pattern 3：The point centered on C points simulates 1 point in the center line of C and D points, simulates 1 point respectively on the outside of C, D；
1 point is simulated to the directions A, vertically respectively simulates 1 point；Matched curve is formed in X, Y, Z axis respectively；
Pattern 4：The point centered on D points simulates 1 point in the center line of C and D points, simulates 1 point respectively on the outside of C, D；
1 point is simulated to the directions B, vertically respectively simulates 1 point；Matched curve is formed in X, Y, Z axis respectively；
9 palace grid space regions are formed by four points, computation model is selected according to actual prediction point region；Region
The case where beyond a computation model, calculates future position temperature value according to different models respectively, is then averaged as pre-
The temperature value of measuring point.As shown in Fig. 2, selecting computation model according to actual prediction point position region, being shown in Table 1, (* is labeled as choosing
The calculating pattern taken)；The case where one computation model is exceeded for region, such as：Point is calculated to fall in the region of Fig. 55, by
Table 1 predicts point value=(M1+M2+M3+M4)/4 it is found that 4 patterns should be calculated；M1, M2, M3, M4 be respectively according to pattern 1,
4 calculated future position temperature value of pattern 2, pattern 3 and pattern.
2. the vapor (steam) temperature at scene/outlet wall temperature measurement point is searched and configuration, vapor (steam) temperature field is established：
It is the steaming of influent header that each heating surface, which has the vapor (steam) temperature measuring point of influent header, every pipe fitting inlet steam temperature,
Stripping temperature；Due to enter each pipe fitting steam flow rate influence it is small, be considered as it is identical, according to outlet wall temperature measurement point temperature value and pipe in
The temperature difference of steam calculates the outlet steam temperature value of every pipe fitting；Then according to every pipe fitting inlet steam temperature, outlet
The length of vapor (steam) temperature and pipe fitting calculates the vapor (steam) temperature each calculated on pipe fitting a little：
(1) each pipe fitting and influent header are established, the relationship of outlet header configures influent header vapor (steam) temperature measuring point；
(2) pipe fitting outlet steam temperature is calculated；
1) selection outlet wall temperature measurement point, it is vertical that measuring point distribution is selected as a cross one；
2) unknown exit wall warm spot is calculated according to known measuring point；
Exit wall temperature value=Wb/Wc × Wd；Wb, Wc are respectively the temperature value of longitudinal known temperature measuring point；Wd is laterally
The temperature value of known temperature measuring point；
3) the vapor (steam) temperature value put on pipe fitting is calculated；
Low fast, the highly endothermic slow feature of outlet steam temperature of absorbing heat of consideration inlet steam temperature, formation temperature and distance
Matched curve；
In T=T goes out-(T go out-T into) × coefficient, according to fair current 0.3, the coefficient of adverse current 0.5 calculates the temperature value of central point；
T goes out：Outlet steam temperature value；T into：Inlet steam temperature value；In T：Central point vapor (steam) temperature value；According to inlet temperature, center
It puts temperature and outlet temperature and generates matched curve along caliber direction distance；Distance according to actual point along caliber direction is brought into
Matched curve generates the vapor (steam) temperature value of this point；
(3) three dimensional space coordinate is calculated:
Each pipe fitting entrance coordinate is configured, according to entrance coordinate, pipe fitting pipe range, corner position, is converted into each calculating point
Spatial position.
3. computer tube tube internal heat exchange coefficient：
Computer tube tube internal heat exchange coefficient formula：
Since each pipeline section coefficient of heat transfer gap of heating surface is smaller, using average intraductal heat exchange coefficient as every on pipeline section
Coefficient of heat transfer value；
Wherein, λ is thermal coefficient, dnFor pipe fitting inside radius, ωqIt is steam flow rate, Pr Pu Lang for kinematic coefficient of viscosity, V
Special number；
Heating surface temperature averages=(inlet steam temperature+outlet steam temperature)/2
Heating surface pressure mean values=(inlet steam pressure+outlet vapor pressure)/2
According to average vapor pressure and temperature computation λ, ωq,V,Pr；
The cross-sectional flow area of steam flow rate V=evaporation capacity (t/h) × 1000 × incidence coefficient/averag density/3600/
Evaporation capacity：Measuring point is evaporated by drum side to measure；
Incidence coefficient：Due to evaporation capacity drum measuring point, there are desuperheating water penetratings, total quantity of steam to be higher than before entering superheater
Drum measuring point；Incidence coefficient=(evaporation capacity+attemperation water flow)/evaporation capacity；
Averag density：Saturated vapor calculation formula is brought into according to pressure and temp and calculates specific volume, and specific volume inverse is density；
Cross-sectional flow area=all pipe fittings of this heating surfaceEvery cross-sectional area is π R × R, and R is pipe fitting internal diameter, Xi
Thus on heating surface single tubular cross-sectional area.
4. the outer coefficient of heat transfer of computer tube：
Heat transfer outside a tube coefficient a1=1/ (R1 × 3.14 × pipe outside diameter × tubular length)；
Wherein, the outer convective heat transfer resistance of R1=pipes external radiation heat exchanged thermoresistance+pipe,
Tubular length=entire heating surface total pipe range
According to heat transfer principle：
Import and export enthalpy difference × evaporation capacity × incidence coefficient=heat flow × tubular length 1.
Heat flow=(average flue-gas temperature-average steam temperature)/thermal resistance is 2.
Dirt heat in the outer convective heat transfer resistance+dust stratification thermal resistance+wall resistance+pipe of thermal resistance=pipe external radiation heat exchanged thermoresistance+pipe
Convective heat transfer resistance is 3. in resistance+pipe
According to 1., 2., R1 is 3. calculated：
(average flue-gas temperature-average steam temperature)/(inlet and outlet enthalpy difference * evaporation capacity * incidence coefficients/tubular length)-
(convective heat transfer resistance in dirtiness resistance+pipe in dust stratification thermal resistance+wall resistance+pipe)=pipe external radiation heat exchanged thermoresistance+heat convection
Thermal resistance=R1
Average flue-gas temperature=heating surface flue-gas temperature central point calculates flue-gas temperature
Average steam temperature=(inlet header box temperature degree mean value+outlet header temperature mean value)/2
Import and export enthalpy difference=H (outlet pressure measuring point, temperature point)-H (inlet pressure measuring point, temperature point)
Evaporation capacity：Measuring point is evaporated by drum side to measure；
Incidence coefficient：Due to evaporation capacity drum measuring point, there are desuperheating water penetratings, total quantity of steam to be higher than before entering superheater
Drum measuring point, incidence coefficient=(evaporation capacity+attemperation water flow)/evaporation capacity；
Tubular length=entire heating surface total pipe range
It tables look-up and obtains dust stratification thermal resistance, coal-fired flue-gas dust stratification thermal resistivity=0.00172m2℃/w。
Dust stratification thermal resistance (DEG C/w)=dust stratification thermal resistivity/π d, d=bores,
Dirtiness resistance (DEG C/w)=fouling resistance coefficient/π d, d=bores, fouling resistance coefficient take 0.0001m2℃/w。
Table look-up to obtain tubing 12CrIMoV, tube wall heat conduction coefficient=33w/ (m DEG C) at 500 DEG C.
Thermal coefficient λ：Under the conditions of steady heat transfer, two 1m apart, area are taken perpendicular to heat conduction direction in interior of articles
For 1 square metre of parallel surface, and 1 degree of the temperature difference of the two panels, then it is flat from a planar conductive to another in 1 second
The heat in face is just defined as the thermal conductivity of the substance, unit：Watt/( meter Du)
Thermal resistance θ=ln (outer diameter/internal diameter)/(2* π λ * L) λ：Tube wall heat conduction coefficient；L length of pipe sections
Wall resistance=ln (outer diameter/internal diameter)/(2*3.14*L* thermal coefficients), L=1m；
Convective heat transfer resistance=1/ (π * bore * length * intraductal heat exchange coefficients in pipe), length=1M
5. simulating wall temperature temperature field：
Outside wall temperature=T flue gases-heat flow × (dust stratification thermal resistance+heat transfer outside a tube thermal resistance) or
Outside wall temperature=T steam+heat flow × (dirtiness resistance+intraductal heat exchange thermal resistance)
Heat flow=(T flue gas-T steam)/thermal resistance
Dirt heat in the outer convective heat transfer resistance+dust stratification thermal resistance+wall resistance+pipe of thermal resistance=pipe external radiation heat exchanged thermoresistance+pipe
Convective heat transfer resistance in resistance+pipe.
By above step, basic three-dimensional boiler model is completed to establish.
Step 2, models for temperature field verification：
1. history big data is verified：According to power plant's history real time data flue-gas temperature/vapor (steam) temperature/exit wall warm spot, with
And power plant's maintenance overhaul record, overtemperature factor and compacted swollen matched curve are formed, according to the record sample of maintenance next time, verification is super
Temperature model correctness.It is mended according to power plant's historical data and calculates temperature field overtemperature situation, determine expert model parameter.
2. maintenance assessment verification：Operation each pipeline section overtemperature data of boiler, Awaiting Overhaul boiler shutdown are calculated according to overtemperature model
Afterwards, it is verified by macro -graph and compacted swollen inspection record in the emphasis overtemperature region calculated according to overtemperature model；With number
Increase according to amount, sample size increases, automatic to correct overtemperature model, forms self learning model, a newest record of examination is as maintenance
The foundation run afterwards executes；Simultaneously according to same unit capacity, used with the data of the same trade or business of the boiler of combustion system with boiler manufacturer
Intersect in lateral model and learn, forms temperature field in furnaces simulation model.
Using above-mentioned simulation model in real time according to measuring point data, in conjunction with expert model, calculate each of pipe fitting put it is real-time
Temperature value, and tubing temperature standard is exceeded according to pipe fitting material, it is recorded as overtemperature, reminds overtemperature region.
Then pipe fitting health status is carried out by being superimposed temperature field.
Claims (8)
1. a kind of method for building up of boiler wear resistant explosion-proof temperature field simulation model, it is characterised in that：According to load, flue-gas temperature,
Vapor (steam) temperature, the outlet practical measuring point of wall temperature simulate the profiling temperatures of each heating surface, pipe, inflection point and pipeline section in real time, pass through
The compacted swollen data that all previous maintenance measures are tested, and each factor of neural network fashion self study high-temperature material leakage failure is passed through
Invalid coefficient, while according to same unit capacity, with boiler manufacturer, with the data of the same trade or business of the boiler of combustion system, for laterally
Model intersects study, is formed and improves Temperature Field Simulation model certainly；Then according to real-time load, flue-gas temperature, vapor (steam) temperature, outlet
The data value binding analysis model of wall temperature measurement point carries out the three-dimensional visualization early warning and analysis in overtemperature region；It includes following steps
Suddenly：
(1), Temperature Field Simulation model foundation：It establishes smoke temperature field, establish vapor (steam) temperature field, computer tube tube internal heat exchange coefficient, computer tube
The outer coefficient of heat transfer, simulation wall temperature temperature field；
(2), models for temperature field is verified：The verification of history big data, maintenance assessment verification；
Described the step of establishing smoke temperature field is：
(1) boiler model coordinate system is determined:Using burner hearth center line and the intersection point of 0 meter of elevation plane of boiler as origin, vertical direction
For y-axis, it is x-axis with screen vertical direction, z-axis is determined according to x-axis and y-axis；
(2) smoke temperature measuring point is chosen:The flue-gas temperature measuring point that each heating surface closes on is chosen, determines measuring point probe co-ordinate, measuring point quantity
No less than three；
(3) according to combustion field Temperature Distribution feature, unknown flue-gas temperature point is simulated in conjunction with known smoke temperature measuring point.
2. the method for building up of boiler wear resistant explosion-proof temperature according to claim 1 field simulation model, it is characterised in that described
Step (3) is：When the known smoke temperature measuring point of selection is four, the first two point is respectively labeled as A, B by flue gas direction, rear two
A point is respectively labeled as C, D, considers that accuracy in computation and temperature spot position, smoke temperature field form 4 kinds of calculating patterns：
Pattern 1：The point centered on A points simulates 1 point in the center line of A and B points, 1 point is simulated respectively in the outside of A, B；To
1 point is simulated in the directions C, vertically respectively simulates 1 point；Matched curve is formed in X, Y, Z axis respectively；
Pattern 2：The point centered on B points simulates 1 point in the center line of A and B points, simulates 1 point respectively on the outside of A, B；To D
1 point is simulated in direction, vertically respectively simulates 1 point；Matched curve is formed in X, Y, Z axis respectively；
Pattern 3：The point centered on C points simulates 1 point in the center line of C and D points, simulates 1 point respectively on the outside of C, D；To A
1 point is simulated in direction, vertically respectively simulates 1 point；Matched curve is formed in X, Y, Z axis respectively；
Pattern 4：The point centered on D points simulates 1 point in the center line of C and D points, simulates 1 point respectively on the outside of C, D；To B
1 point is simulated in direction, vertically respectively simulates 1 point；Matched curve is formed in X, Y, Z axis respectively；
9 palace grid space regions are formed by four points, computation model is selected according to actual prediction point region；Region exceeds
The case where one computation model, calculates future position temperature value according to different models respectively, is then averaged as future position
Temperature value.
3. the method for building up of boiler wear resistant explosion-proof temperature according to claim 1 field simulation model, it is characterised in that described
Establishing vapor (steam) temperature field is：
It is the steam temperature of influent header that each heating surface, which has the vapor (steam) temperature measuring point of influent header, every pipe fitting inlet steam temperature,
Degree；According to the temperature difference of steam in the temperature value and pipe of outlet wall temperature measurement point, the outlet steam temperature value of every pipe fitting is calculated；
Then it according to the length of every pipe fitting inlet steam temperature, outlet steam temperature and pipe fitting, calculates and is each calculated on pipe fitting
The vapor (steam) temperature of point：
(1) each pipe fitting and influent header are established, the relationship of outlet header configures influent header vapor (steam) temperature measuring point；
(2) pipe fitting outlet steam temperature is calculated；
1) selection outlet wall temperature measurement point, it is vertical that measuring point distribution is selected as a cross one；
2) unknown exit wall warm spot is calculated according to known measuring point；
Exit wall temperature value=Wb/Wc × Wd；Wb, Wc are respectively the temperature value of longitudinal known temperature measuring point；Wd is known to transverse direction
The temperature value of temperature point；
3) the vapor (steam) temperature value put on pipe fitting is calculated；
Consider that the low heat absorption of inlet steam temperature is fast, the highly endothermic slow feature of outlet steam temperature forms the fitting of temperature and distance
Curve；
In T=T goes out-(T go out-T into) × coefficient, according to fair current 0.3, the coefficient of adverse current 0.5 calculates the temperature value of central point；T goes out：
Outlet steam temperature value；T into：Inlet steam temperature value；In T：Central point vapor (steam) temperature value；According to inlet temperature, central point temperature
Degree generates matched curve with outlet temperature and along caliber direction distance；Distance according to actual point along caliber direction brings fitting into
Curve generates the vapor (steam) temperature value of this point；
(3) three dimensional space coordinate is calculated:
Each pipe fitting entrance coordinate is configured, according to entrance coordinate, pipe fitting pipe range, corner position, is converted into each sky for calculating point
Between position.
4. the method for building up of boiler wear resistant explosion-proof temperature according to claim 1 field simulation model, it is characterised in that described
The step of computer tube tube internal heat exchange coefficient is：
Computer tube tube internal heat exchange coefficient formula：
Since each pipeline section coefficient of heat transfer gap of heating surface is smaller, using average intraductal heat exchange coefficient as every on pipeline section heat exchange
Coefficient value；
Wherein, λ is thermal coefficient, dnFor pipe fitting inside radius, ωqFor kinematic coefficient of viscosity, V be steam flow rate, Pr is Prandtl
Number；
Heating surface temperature averages=(inlet steam temperature+outlet steam temperature)/2
Heating surface pressure mean values=(inlet steam pressure+outlet vapor pressure)/2
According to average vapor pressure and temperature computation λ, ωq,V,Pr；
The cross-sectional flow area of steam flow rate V=evaporation capacity (t/h) × 1000 × incidence coefficient/averag density/3600/
Evaporation capacity：Measuring point is evaporated by drum side to measure；
Incidence coefficient：Due to evaporation capacity drum measuring point, there are desuperheating water penetratings, total quantity of steam to be higher than drum before entering superheater
Measuring point；Incidence coefficient=(evaporation capacity+attemperation water flow)/evaporation capacity；
Averag density：Saturated vapor calculation formula is brought into according to pressure and temp and calculates specific volume, and specific volume inverse is density；
Cross-sectional flow area=all pipe fittings of this heating surfaceEvery cross-sectional area is π R × R, and R is pipe fitting internal diameter, XiThus
The cross-sectional area of single tubular on heating surface.
5. the method for building up of boiler wear resistant explosion-proof temperature according to claim 1 field simulation model, it is characterised in that described
Coefficient of heat transfer step is outside computer tube：
Heat transfer outside a tube coefficient a1=1/ (R1 × 3.14 × pipe outside diameter × tubular length)；
Wherein, the outer convective heat transfer resistance of R1=pipes external radiation heat exchanged thermoresistance+pipe,
Tubular length=entire heating surface total pipe range
According to heat transfer principle：
Import and export enthalpy difference × evaporation capacity × incidence coefficient=heat flow × tubular length 1.
Heat flow=(average flue-gas temperature-average steam temperature)/thermal resistance is 2.
Dirtiness resistance+pipe in the outer convective heat transfer resistance+dust stratification thermal resistance+wall resistance+pipe of thermal resistance=pipe external radiation heat exchanged thermoresistance+pipe
Interior convective heat transfer resistance is 3.
According to 1., 2., R1 is 3. calculated：
(average flue-gas temperature-average steam temperature)/(inlet and outlet enthalpy difference * evaporation capacity * incidence coefficients/tubular length)-(dust stratification
Convective heat transfer resistance in dirtiness resistance+pipe in thermal resistance+wall resistance+pipe)=pipe external radiation heat exchanged thermoresistance+convective heat transfer resistance=
R1
Average flue-gas temperature=heating surface flue-gas temperature central point calculates flue-gas temperature
Average steam temperature=(inlet header box temperature degree mean value+outlet header temperature mean value)/2
Import and export enthalpy difference=H (outlet pressure measuring point, temperature point)-H (inlet pressure measuring point, temperature point)
Evaporation capacity：Measuring point is evaporated by drum side to measure；
Incidence coefficient：Due to evaporation capacity drum measuring point, there are desuperheating water penetratings, total quantity of steam to be higher than drum before entering superheater
Measuring point, incidence coefficient=(evaporation capacity+attemperation water flow)/evaporation capacity；
Tubular length=entire heating surface total pipe range
Dust stratification thermal resistance=dust stratification thermal resistivity/π d, d is bore；
Dirtiness resistance (DEG C/w)=fouling resistance coefficient/π d, d is bore；
Thermal resistance θ=ln (outer diameter/internal diameter)/(2 π λ L) λ：Thermal coefficient；L length of pipe sections
Wall resistance=ln (outer diameter/internal diameter)/(2 × 3.14 × L × thermal coefficient)；
6. the method for building up of boiler wear resistant explosion-proof temperature according to claim 1 field simulation model, it is characterised in that described
Simulating wall temperature temperature field step is：
Outside wall temperature=T flue gases-heat flow × (dust stratification thermal resistance+heat transfer outside a tube thermal resistance) or
Outside wall temperature=T steam+heat flow × (dirtiness resistance+intraductal heat exchange thermal resistance)
Heat flow=(T flue gas-T steam)/thermal resistance
Dirtiness resistance+pipe in the outer convective heat transfer resistance+dust stratification thermal resistance+wall resistance+pipe of thermal resistance=pipe external radiation heat exchanged thermoresistance+pipe
Interior convective heat transfer resistance.
7. the method for building up of boiler wear resistant explosion-proof temperature according to claim 1 field simulation model, it is characterised in that described
History big data is verified as：It is tieed up according to power plant's history real time data flue-gas temperature/vapor (steam) temperature/exit wall warm spot and power plant
Inspection record is repaiied, overtemperature factor and compacted swollen matched curve are formed, according to the record sample of maintenance next time, verification overtemperature model is just
True property.
8. the method for building up of boiler wear resistant explosion-proof temperature according to claim 1 field simulation model, it is characterised in that described
Maintenance assessment is verified as：It is calculated according to overtemperature model and runs each pipeline section overtemperature data of boiler, after Awaiting Overhaul boiler shutdown, according to super
It is verified by macro -graph and compacted swollen inspection record in the emphasis overtemperature region that temperature model calculates；As data volume increases,
Sample size increases, automatic to correct overtemperature model, forms self learning model, a newest record of examination as run after overhauling according to
According to execution；Simultaneously according to same unit capacity, with boiler manufacturer, with the data of the same trade or business of the boiler of combustion system, for lateral model
Intersect study, forms temperature field in furnaces simulation model.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610866742.1A CN106524118B (en) | 2016-09-30 | 2016-09-30 | The method for building up of boiler wear resistant explosion-proof temperature field model |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610866742.1A CN106524118B (en) | 2016-09-30 | 2016-09-30 | The method for building up of boiler wear resistant explosion-proof temperature field model |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106524118A CN106524118A (en) | 2017-03-22 |
| CN106524118B true CN106524118B (en) | 2018-11-13 |
Family
ID=58344690
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610866742.1A Active CN106524118B (en) | 2016-09-30 | 2016-09-30 | The method for building up of boiler wear resistant explosion-proof temperature field model |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106524118B (en) |
Families Citing this family (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109635341A (en) * | 2018-11-14 | 2019-04-16 | 东方电气集团东方锅炉股份有限公司 | A kind of pressure containing part life-span prediction method of three flue double reheat boiler of tail portion |
| CN109966996A (en) * | 2019-02-26 | 2019-07-05 | 武汉恒力华振科技有限公司 | A kind of system using the melting state in big data analysis prediction hot melt adhesive production process |
| CN111120988A (en) * | 2019-12-11 | 2020-05-08 | 中国大唐集团科学技术研究院有限公司火力发电技术研究院 | Boiler heating surface pipe wall overtemperature early warning method based on hearth temperature field distribution |
| CN111207381A (en) * | 2020-01-13 | 2020-05-29 | 华能莱芜发电有限公司 | Anti-abrasion and anti-explosion monitoring method and system for heating surface of power plant boiler |
| CN111578254B (en) * | 2020-05-29 | 2022-07-12 | 华润电力技术研究院有限公司 | Boiler temperature field distribution prediction method, device and system |
| CN114165777B (en) * | 2020-09-10 | 2023-10-24 | 河北云酷科技有限公司 | Intelligent recognition model for four-pipe leakage of power plant boiler |
| CN112382158A (en) * | 2020-10-17 | 2021-02-19 | 大唐国际发电股份有限公司张家口发电厂 | Boiler anti-abrasion and anti-explosion training system based on three-dimensional simulation system |
| CN112949226B (en) * | 2021-03-25 | 2022-10-04 | 南方电网电力科技股份有限公司 | Prediction method, device and equipment for tube panel wall temperature deviation of high-temperature heating surface of boiler |
| CN114110550B (en) * | 2021-11-12 | 2023-08-29 | 国电内蒙古东胜热电有限公司 | Boiler temperature prediction method and system |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101476715A (en) * | 2009-01-20 | 2009-07-08 | 西安交通大学 | Early warning method for failure of water-cooling wall of power boiler |
| CN105674237A (en) * | 2016-02-22 | 2016-06-15 | 西安西热电站信息技术有限公司 | Real-time measuring method for smoke temperature fields of hearth outlet and horizontal flue of power station boiler |
| CN105698158A (en) * | 2015-11-23 | 2016-06-22 | 大唐鸡西第二热电有限公司 | Method for detecting four-tube leak of circulating fluidized bed boiler |
| CN105758661A (en) * | 2016-05-05 | 2016-07-13 | 华电国际电力股份有限公司技术服务中心 | System and method for evaluating service life of boiler heating surface |
-
2016
- 2016-09-30 CN CN201610866742.1A patent/CN106524118B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101476715A (en) * | 2009-01-20 | 2009-07-08 | 西安交通大学 | Early warning method for failure of water-cooling wall of power boiler |
| CN105698158A (en) * | 2015-11-23 | 2016-06-22 | 大唐鸡西第二热电有限公司 | Method for detecting four-tube leak of circulating fluidized bed boiler |
| CN105674237A (en) * | 2016-02-22 | 2016-06-15 | 西安西热电站信息技术有限公司 | Real-time measuring method for smoke temperature fields of hearth outlet and horizontal flue of power station boiler |
| CN105758661A (en) * | 2016-05-05 | 2016-07-13 | 华电国际电力股份有限公司技术服务中心 | System and method for evaluating service life of boiler heating surface |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106524118A (en) | 2017-03-22 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106524118B (en) | The method for building up of boiler wear resistant explosion-proof temperature field model | |
| CN102494714B (en) | Synchronous reckoning method of utility boiler efficiency and coal heat value as well as ash content and moisture content | |
| CN103646176B (en) | The comprehensive calculation method of energy-saving effect after turbine steam seal transformation | |
| Zhao et al. | Productivity model in hot and humid environment based on heat tolerance time analysis | |
| CN102494325B (en) | Method for monitoring intra-furnace dynamic wall temperature in high-temperature tube system of power station boiler | |
| US20130060471A1 (en) | Estimating energy savings and carbon offsets for buildings in real-time | |
| CN104808634B (en) | A kind of thermal power plant's Visualized management system based on virtual reality | |
| CN102444885B (en) | Method for avoiding overheat tube burst in tubular boiler of power station boiler | |
| Kim et al. | Heating and cooling energy consumption prediction model for high-rise apartment buildings considering design parameters | |
| CN101319256A (en) | Intelligent monitoring method for cooling wall of blast furnace | |
| CN113361192A (en) | Boiler internal heating surface wall temperature safety monitoring and evaluating system | |
| CN106547945A (en) | A kind of energy efficiency test method for being applied to trilogy supply region energy supplying system | |
| CN103390211A (en) | Thermal generator set overall process energy management and cost analysis system | |
| Buratti et al. | Comparison of the energy performance of existing buildings by means of dynamic simulations and artificial neural networks | |
| Wang | Accuracy, validity and relevance of probabilistic building energy models | |
| Zavrl et al. | Comparison of simulated and monitored energy performance indicators on NZEB case study eco silver house | |
| Mayr et al. | Computational analysis of a semi-industrial furnace fired by a flat flame burner under different O2/N2 ratios using the steady laminar flamelet approach | |
| Rospi et al. | Analysis of energy consumption of different typologies of school buildings in the city of Matera (Southern Italy) | |
| Liu et al. | Research on heating performance of heating radiator at low temperature | |
| CN108647478B (en) | BP neural network-based SCR inlet smoke temperature online calculation method for cogeneration unit | |
| KR20130089111A (en) | Estimating method for heat loss of secondary pipeline for district heating of apartment complex | |
| Bisetto et al. | Fire tube heat generators: Experimental analysis and modeling | |
| Duanmu et al. | Heat transfer model of hot-wall Kang based on the non-uniform Kang surface temperature in Chinese rural residences | |
| CN108073763A (en) | A kind of measuring method of station boiler unburned carbon in flue dust | |
| CN103697958A (en) | Real-time measurement method of mass flow of saturated steam of drum outlet of coal-fired unit |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | 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 |