CN109284581B - Method for analyzing matching characteristic of boiler heating surface working medium flow and heat load distribution - Google Patents
Method for analyzing matching characteristic of boiler heating surface working medium flow and heat load distribution Download PDFInfo
- Publication number
- CN109284581B CN109284581B CN201811388758.1A CN201811388758A CN109284581B CN 109284581 B CN109284581 B CN 109284581B CN 201811388758 A CN201811388758 A CN 201811388758A CN 109284581 B CN109284581 B CN 109284581B
- Authority
- CN
- China
- Prior art keywords
- enthalpy rise
- factors
- heating surface
- working medium
- distribution
- 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
- 238000010438 heat treatment Methods 0.000 title claims abstract description 72
- 238000000034 method Methods 0.000 title claims abstract description 39
- 238000001816 cooling Methods 0.000 claims abstract description 41
- 230000002159 abnormal effect Effects 0.000 claims description 44
- 230000005012 migration Effects 0.000 claims description 36
- 238000013508 migration Methods 0.000 claims description 36
- 239000003245 coal Substances 0.000 claims description 16
- 230000003068 static effect Effects 0.000 claims description 13
- 238000007619 statistical method Methods 0.000 claims description 9
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 9
- 230000008859 change Effects 0.000 claims description 7
- 108010074864 Factor XI Proteins 0.000 claims description 6
- 238000004364 calculation method Methods 0.000 claims description 5
- 238000004939 coking Methods 0.000 claims description 4
- 238000009825 accumulation Methods 0.000 claims description 3
- 238000012360 testing method Methods 0.000 claims description 3
- 238000010521 absorption reaction Methods 0.000 description 4
- UGFAIRIUMAVXCW-UHFFFAOYSA-N Carbon monoxide Chemical compound [O+]#[C-] UGFAIRIUMAVXCW-UHFFFAOYSA-N 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 2
- 239000002826 coolant Substances 0.000 description 2
- 239000003546 flue gas Substances 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 230000000630 rising effect Effects 0.000 description 2
- 230000005856 abnormality Effects 0.000 description 1
- 230000006978 adaptation Effects 0.000 description 1
- 230000003321 amplification Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000003199 nucleic acid amplification method Methods 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/08—Thermal analysis or thermal optimisation
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
- Y04S—SYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
- Y04S10/00—Systems supporting electrical power generation, transmission or distribution
- Y04S10/50—Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Testing Resistance To Weather, Investigating Materials By Mechanical Methods (AREA)
- Control Of Steam Boilers And Waste-Gas Boilers (AREA)
Abstract
The invention discloses a method for analyzing the distribution matching characteristic of a working medium flow and a heat load on a heating surface of a boiler, which comprises the steps of calculating a matching factor of a cooling working medium flow and the heat load on the heating surface, and grading the mismatching condition; and calculating the influence of each influence factor on the matching factor, grading according to the influence degree, analyzing the mismatching degree of the flow of the heating surface cooling working medium and the heat load, and diagnosing the significant factor causing the mismatching of the flow of the heating surface cooling working medium and the heat load so as to determine the mismatching sequence of the distribution of the heating surface working medium and the heat load of the treatment boiler. The invention can diagnose the controllable significant influence factors, the controllable semi-significant influence factors, the semi-controllable significant influence factors and the semi-controllable semi-significant influence factors which are not matched with the heating surface by analyzing the matching characteristics of the flow distribution and the heat load distribution of the cooling working medium on the heating surface of the boiler, and quickly and accurately determine the treatment order of the mismatching of the flow distribution and the heat load distribution of the cooling working medium on the heating surface.
Description
Technical Field
The invention belongs to an analysis method for matching boiler heating surface with heat absorption capacity, namely heat load distribution with positive feedback characteristics, and particularly relates to a method for analyzing matching characteristics of boiler heating surface working medium flow and heat load distribution.
Background
The superheater and the reheater of the boiler and the water wall of the once-through boiler have typical positive feedback characteristics, the heat load of the heat exchange tube is larger, the outlet working medium temperature of the heat exchange tube is higher, the density of the working medium is low, so that the flow of the cooling working medium in the tube is reduced, and the outlet working medium temperature is higher. During design, a method of controlling the flow resistance coefficient of the heat exchange tube to be matched with heat absorption capacity is adopted to enable the temperature of the outlet working medium to be basically consistent, and the working medium often shows larger difference in operation. At present, the distribution deviation of the outlet wall temperature is analyzed more. However, a method for forming a system is not provided for expressing the matching degree of the working medium flow of the heating surface and the heat load, and a special method is not provided for diagnosing the reason of mismatching of the working medium flow of the heating surface and the heat load.
At present, the condition of large temperature distribution deviation of the outlet wall of a heating surface of a boiler is prominent, and the safety of the operation of the heating surface of the boiler is greatly influenced, such as the conditions that part of heat exchange tube materials are aged too fast, part of oxide skins on the inner wall of the heat exchange tube grow too fast, the bonding state with a base material is not good, a water wall tube is pulled to crack, and a large area of transverse cracks are formed on the water wall tube. The reason for causing the large wall temperature distribution deviation is that the flow of the cooling working medium in the heat exchange tube is not matched with the heat absorption capacity of the heat exchange tube, and the factors causing the mismatching are as follows: the flow resistance coefficient of the heat exchange tubes is distributed; deviation of heat absorption capacity of different heated tubes; the temperature and flow distribution of the flue gas in the furnace, namely the heat load distribution; coking and scaling on the surface of the heating surface; the flue gas flow and temperature distribution change caused by the boiler heat load change; the ratio of radiant heat exchange to convective heat exchange caused by heat load is changed; the temperature distribution of the flame of the hearth is changed due to the change of the heat load of the boiler; the temperature distribution of the flame of the hearth is changed due to the switching of the coal mill; static pressure distribution of inlet and outlet header; the distribution of the steam-water system resistance coefficient deviates from the theory; orifice fouling, etc. The expression of the matching characteristic of the flow of the cooling working medium in the heat exchange tube of the heating surface and the heat load, the amplitude and the speed of steam fluctuation in the variable load process, the amplification effect of the wall temperature fluctuation compared with the steam temperature and the like are also provided. The change of some factors is difficult to control and belongs to uncontrollable factors; some factors can be changed and belong to controllable factors; some factors can be partially changed by adjusting control, and belong to semi-controllable factors.
Disclosure of Invention
The technical problem to be solved by the invention is to provide a method for analyzing the matching characteristic of the boiler heating surface working medium flow and the heat load distribution aiming at the defects of the prior art.
In order to achieve the technical purpose, the technical scheme adopted by the invention is as follows:
the method for analyzing the matching characteristic of the distribution of the working medium flow and the heat load of the heating surface of the boiler calculates the matching factor of the cooling working medium flow and the heat load of the heating surface and grades the mismatching condition; calculating the influence of each influence factor on the matching factor, grading according to the influence degree, analyzing the mismatching degree of the heating surface cooling working medium flow and the heat load, and diagnosing the significant factor causing the mismatching of the heating surface cooling working medium flow and the heat load so as to determine the mismatching sequence of the heating surface working medium flow and the heat load distribution of the treatment boiler, wherein the method specifically comprises the following steps:
the method comprises the following steps: testing the wall temperature distribution of the outlet of the heating surface under different operating conditions;
step two: calculating the enthalpy rise distribution of the working medium of each heated pipe according to the operation parameters of the boiler and the wall temperature distribution;
step three: counting the distribution of enthalpy rise, and calculating average enthalpy rise, standard enthalpy rise deviation, maximum enthalpy rise deviation and relative maximum enthalpy rise deviation; abnormal enthalpy rise concentration ratio, maximum enthalpy rise migration and abnormal enthalpy rise distribution migration;
calculating average enthalpy rise, standard enthalpy rise deviation, maximum enthalpy rise deviation and relative maximum enthalpy rise deviation; the formulas of the abnormal enthalpy rise concentration ratio, the maximum enthalpy rise migration and the abnormal enthalpy rise distribution migration are as follows:
mean enthalpy rise:
standard enthalpy rise deviation:
enthalpy rise deviation, i.e., the difference between the heat exchange tube enthalpy rise and the average enthalpy rise: Δ Hij＝Hij-Hp kJ/kg (3)；
Maximum enthalpy rise deviation: Δ Hmax＝maxHij-Hp kJ/kg (4)；
Relative maximum enthalpy rise deviation, i.e. maximum enthalpy rise deviation to average enthalpy rise ratio: h ismax＝Hmax/Hp (5)；
Standard relative enthalpy rise deviation, i.e., the ratio of heat exchange tube enthalpy rise deviation to standard enthalpy rise deviation: h isδ＝Hp/δH (6)；
Maximum standard enthalpy rise deviation, i.e. the ratio of maximum enthalpy rise deviation to standard enthalpy rise deviation: h isδmax＝Hmax/δH (7)；
Continuous abnormal enthalpy rise concentration ratio gammai: counting that standard enthalpy rise deviation of the heat exchange tube is more than 3 times deltaHThe interval between the two statistical points is not more than two, the pipe is regarded as continuous abnormality, the percentage of the continuous abnormal points in the statistical sample is called as continuous abnormal enthalpy rising concentration ratio, and the formula is as follows: gamma rayi＝100×∑Cont(h δ max＞3)/N％ (8)；
Maximum enthalpy rise migration τkn: under the working condition k and the working condition n, the number of statistical points of the position migration of the maximum enthalpy riser tube is called the maximum enthalpy riser migration in percentage of the statistical sample number, and the formula is as follows: tau iskn＝100×(S(Hmax)k-S(Hmax)n)/N％ (9)；
Wherein, S (Hmax)kThe position point of the water wall tube with the largest enthalpy rise under the operation working condition with the number of k is shown;
the central position Cen of the enthalpy rise point, i.e. Cont (h)δmax>3) A symmetrical center position point of continuous abnormal enthalpy rise distribution;
abnormal enthalpy rise distribution migration τγ kn: counting the central position of the abnormal enthalpy rise point which is continuously distributed, wherein the percentage of the number of the statistical points of the central position migration in the statistical sample number under the working condition k and the working condition n is called abnormal enthalpy rise distribution migration, and the formula is as follows: tau isγkn＝100×(Cenk(Cont(hδmax>3))-Cenn(Cont(hδmax>3)))/N％ (10)。
Step four: calculating a matching factor of the flow of the cooling working medium on the heating surface and the heat load by adopting a statistical analysis method of a weighting factor; piece of clothThe matching factor is the matching property of the flow of the cooling working medium on the heating surface and the heat load, and is expressed by the distribution of enthalpy rise of the heat exchange tube and is marked as xi; the matching factor can be divided into a single working condition matching factor and a multi-working condition matching factor, and the multi-working condition matching factor xiht k,nThe functions of average enthalpy rise, standard enthalpy rise deviation, maximum enthalpy rise deviation, relative maximum enthalpy rise deviation, standard maximum enthalpy rise deviation, abnormal enthalpy rise concentration ratio, maximum enthalpy rise migration and abnormal enthalpy rise distribution migration are calculated by adopting a weighted statistical method; the single working condition matching factor xikOnly considering the matching characteristics under a single operating condition, the method is a function of average enthalpy rise, standard enthalpy rise deviation, maximum enthalpy rise deviation, relative maximum enthalpy rise deviation, standard maximum enthalpy rise deviation and abnormal enthalpy rise concentration, and a weighted statistical method is also adopted for calculation, wherein the specific calculation formula is as follows:
ξht k＝f1(Hp,δH,Hmax,hmax,hδmax,γi)
＝α1×Hp+α2×δH+α3×Hmax+α4×hmax+α5×hδmax+α6×γi (11)；
ξht k,n＝f2(Hp,δH,Hmax,hmax,hδmax,γi,τkn,τγkn)＝β1×(Hp-Hp)+β2×(δk H-δn H)+β3×(Hk max -Hn max)+β4×(hk max-hn max)+β5×(hk δmax-hn δmax)+β6×(γk i-γn i)+β7×τkn+β8×τγkn(12) in the formula, α1、α2、α3、α4、α5、α6、β1、β2、β3、β4、β5、β6、β7And beta8Are all weighting coefficients.
Step five: grading the matching characteristics of the flow of the cooling working medium on the heating surface and the heat load according to the calculated matching factors; the grading method comprises the following steps: given standard matching factor xi01、ξ02When the calculated heating surface matching factor xiht>ξ01Time, significant mismatch; when xi02>ξht>ξ01Half significant mismatch; when xiht<ξ02When it is, it is a match.
Step six: grading factors influencing the matching characteristic of the flow of the cooling working medium on the heating surface and the heat load; specifically, the factors affecting the flow distribution of the cooling medium include: the flow resistance coefficient of the heat exchange tube, the aperture of a throttling hole, the static pressure distribution of an inlet header, the static pressure distribution of an outlet header, the output of a boiler and the blockage condition of foreign matters; factors that influence the distribution of the thermal load are: boiler output, hearth coking, heating surface ash accumulation, coal mill operation, primary air uniformity, secondary air door adjustment and coal quality of coal; the grading result is as follows: the uncontrollable factors are: the heat exchange tube flow resistance coefficient, inlet header static pressure distribution, outlet header static pressure distribution and coal quality; the controllable factors are as follows: the diameter of a throttling hole is used when the coal mill is put into use; the remaining factors are semi-controllable factors.
Step seven: calculating the influence of the factors on the matching factors of the flow of the cooling working medium on the heating surface and the heat load, and grading the influence degree of the factors; the method for grading the influence degree of the factors comprises the following steps: calculating the matching factor xi of the heating surface cooling working medium flow and the heat load of the controllable factor and the semi-controllable factor through different working conditions influencing the change of the factorsht k,nMatching the criterion factor xief01、ξef02Comparing, grading the influence of the factors on the matching factors, xiht k,n>ξef01Time, is a significant influencing factor; when xief02>ξht k,n>ξef01Time, is a semi-significant influencing factor; when xiht k,n<ξef02Is not significantAnd (4) influencing factors.
Step eight: grading the influence factors; according to the controllable characteristics of the influence factors and the influence degree on the matching factors, the influence factors are divided into: controllable significant influence factors, semi-controllable significant influence factors and non-controllable significant influence factors; controllable semi-significant influencing factors, semi-controllable semi-significant influencing factors and non-controllable semi-significant influencing factors; no significant influencing factor.
Step nine: the treatment is carried out according to controllable obvious influence factors, controllable semi-obvious influence factors, semi-controllable obvious influence factors and semi-controllable semi-obvious influence factors when the working medium flow of the heating surface of the boiler is not matched with the heat load distribution, and the treatment is not considered for other factors.
The invention has the following beneficial effects:
the invention can diagnose the controllable obvious influence factors, the controllable semi-obvious influence factors, the semi-controllable obvious influence factors and the semi-controllable semi-obvious influence factors which are not matched with the heating surface by analyzing the matching characteristics of the flow distribution and the heat load distribution of the cooling working medium of the heating surface of the boiler, and determine the treatment sequence of the mismatching of the flow distribution and the heat load distribution of the cooling working medium of the heating surface.
Detailed Description
The present invention will be described in further detail with reference to specific examples.
The overall idea of the invention is as follows: calculating a matching factor of the flow of the cooling working medium on the heating surface and the heat load, and grading the mismatching condition; and calculating the influence of each influence factor on the matching factor, and grading according to the influence. Therefore, the mismatching degree of the flow of the cooling working medium on the heating surface and the heat load is analyzed, and the significant factor causing the mismatching of the flow of the cooling working medium on the heating surface and the heat load is diagnosed.
In the embodiment, the method provided by the invention is used for realizing the matching characteristic analysis of the cooling working medium flow of the water cooled wall of the once-through boiler adopting the vertical water cooled wall for the lower radiation and the heat load of the hearth, and specifically comprises the following steps:
the method comprises the following steps: testing the wall temperature distribution of the outlet of the heating surface under different operating conditions;
step two: calculating the enthalpy rise distribution of the working medium of each heated pipe according to the operation parameters of the boiler and the wall temperature distribution;
step three: counting the distribution of enthalpy rise, and calculating average enthalpy rise, standard enthalpy rise deviation, maximum enthalpy rise deviation and relative maximum enthalpy rise deviation; abnormal enthalpy rise concentration ratio, maximum enthalpy rise migration and abnormal enthalpy rise distribution migration;
calculating average enthalpy rise, standard enthalpy rise deviation, maximum enthalpy rise deviation and relative maximum enthalpy rise deviation; the formulas of the abnormal enthalpy rise concentration ratio, the maximum enthalpy rise migration and the abnormal enthalpy rise distribution migration are as follows:
mean enthalpy rise:
standard enthalpy rise deviation:
enthalpy rise deviation, i.e., the difference between the heat exchange tube enthalpy rise and the average enthalpy rise: Δ Hij＝Hij-Hp kJ/kg (3)；
Maximum enthalpy rise deviation: Δ Hmax＝maxHij-Hp kJ/kg (4)；
Relative maximum enthalpy rise deviation, i.e. maximum enthalpy rise deviation to average enthalpy rise ratio: h ismax＝Hmax/Hp (5)；
Standard relative enthalpy rise deviation, i.e., the ratio of enthalpy rise deviation of the heat exchange tube to standard enthalpy rise deviation: h isδ＝Hp/δH (6)；
Maximum standard enthalpy rise deviation, i.e. the ratio of maximum enthalpy rise deviation to standard enthalpy rise deviation: h isδmax＝Hmax/δH (7)；
Continuous abnormal enthalpy rise concentration ratio gammai: counting that standard enthalpy rise deviation of the heat exchange tube is more than 3 times deltaHThe interval between the two statistical points is not more than two, and the continuous abnormal points account for the percentage of the statistical sample and are called continuous abnormal enthalpy rising setMedium, the formula is: gamma rayi＝100×∑Cont(h δ max＞3)/N％ (8)；
Maximum enthalpy rise migration τkn: under the working condition k and the working condition n, the number of statistical points of the position migration of the maximum enthalpy riser tube is called the maximum enthalpy riser migration in percentage of the statistical sample number, and the formula is as follows: tau iskn＝100×(S(Hmax)k-S(Hmax)n)/N％ (9)；
Wherein, S (Hmax)kThe position point of the water wall tube with the largest enthalpy rise under the operation working condition with the number of k is shown;
the central position Cen of the enthalpy rise point, i.e. Cont (h)δmax>3) A symmetrical center position point of continuous abnormal enthalpy rise distribution;
abnormal enthalpy rise distribution migration tauγ kn : counting the central position of the abnormal enthalpy rise point which is continuously distributed, wherein the percentage of the number of the statistical points of the central position migration in the statistical sample number under the working condition k and the working condition n is called abnormal enthalpy rise distribution migration, and the formula is as follows: tau.γkn＝100×(Cenk(Cont(hδmax>3))-Cenn(Cont(hδmax>3)))/N％ (10)。
Step four: calculating a matching factor of the flow of the cooling working medium on the heating surface and the heat load by adopting a statistical analysis method of a weighting factor; the matching factor is the matching performance of the flow of the cooling working medium on the heating surface and the heat load, and is expressed by the distribution of enthalpy rise of the heat exchange tube and is marked as xi; the matching factor can be divided into a single working condition matching factor and a multi-working condition matching factor, and the multi-working condition matching factor xiht k,nThe functions of average enthalpy rise, standard enthalpy rise deviation, maximum enthalpy rise deviation, relative maximum enthalpy rise deviation, standard maximum enthalpy rise deviation, abnormal enthalpy rise concentration ratio, maximum enthalpy rise migration and abnormal enthalpy rise distribution migration are calculated by adopting a weighted statistical method; the single working condition matching factor xikOnly considering the matching characteristic under a single operation condition, the method is a function of average enthalpy rise, standard enthalpy rise deviation, maximum enthalpy rise deviation, relative maximum enthalpy rise deviation, standard maximum enthalpy rise deviation and abnormal enthalpy rise concentration, and a weighted statistical method is also adopted for calculation, wherein the specific calculation formula is as follows:
ξht k＝f1(Hp,δH,Hmax,hmax,hδmax,γi)
＝α1×Hp+α2×δH+α3×Hmax+α4×hmax+α5×hδmax+α6×γi(11)；
ξht k,n＝f2(Hp,δH,Hmax,hmax,hδmax,γi,τkn,τγkn)＝β1×(Hp-Hp)+β2×(δk H-δn H)+β3×(Hk max -Hn max)+β4×(hk max-hn max)+β5×(hk δmax-hn δmax)+β6×(γk i-γn i)+β7×τkn+β8×τγkn(12) in the formula, α1、α2、α3、α4、α5、α6、β1、β2、β3、β4、β5、β6、β7And beta8Are all weighting coefficients.
Step five: grading the matching characteristics of the flow of the cooling working medium on the heating surface and the heat load according to the calculated matching factors; the grading method comprises the following steps: given standard matching factor xi01、ξ02When the calculated heating surface matching factor xiht>ξ01Time, significant mismatch; when xi02>ξht>ξ01When, it is a half significant mismatch; when xiht<ξ02When it is, it is a match.
Step six: grading factors influencing the matching characteristic of the flow of the cooling working medium on the heating surface and the heat load; specifically, the factors affecting the flow distribution of the cooling medium include: the flow resistance coefficient of the heat exchange tube, the aperture of a throttling hole, the static pressure distribution of an inlet header, the static pressure distribution of an outlet header, the output of a boiler and the blockage condition of foreign matters; factors that influence the distribution of the thermal load are: boiler output, hearth coking, heating surface ash accumulation, coal mill operation, primary air uniformity, secondary air door adjustment and coal quality of coal; the grading result is as follows: the uncontrollable factors are: the heat exchange pipe flow resistance coefficient, the inlet header static pressure distribution, the outlet header static pressure distribution and the coal quality; the controllable factors are as follows: the diameter of a throttling hole is used when the coal mill is put into use; the remaining factors are semi-controllable factors.
Step seven: calculating the influence of the factors on the matching factors of the flow of the cooling working medium on the heating surface and the heat load, and grading the influence degree of the factors; the method for grading the influence degree of the factors comprises the following steps: calculating the matching factor xi of the heating surface cooling working medium flow and the heat load of the controllable factor and the semi-controllable factor through different working conditions influencing the change of the factorsht k,nMatching the criterion factor xief01、ξef02Comparing, grading the influence of the factors on the matching factors, xiht k,n>ξef01Time, is a significant influencing factor; when xief02>ξht k,n>ξef01Time, is a semi-significant influencing factor; when xiht k,n<ξef02Time is an insignificant influence.
Step eight: grading the influence factors; according to the controllable characteristics of the influence factors and the influence degree on the matching factors, the influence factors are divided into: controllable significant influence factors, semi-controllable significant influence factors and non-controllable significant influence factors; controllable semi-significant influencing factors, semi-controllable semi-significant influencing factors and non-controllable semi-significant influencing factors; and (4) no significant influence factor.
Step nine: the treatment is carried out according to controllable obvious influence factors, controllable semi-obvious influence factors, semi-controllable obvious influence factors and semi-controllable semi-obvious influence factors when the working medium flow of the heating surface of the boiler is not matched with the heat load distribution, and the treatment is not considered for other factors.
The above is only a preferred embodiment of the present invention, and the protection scope of the present invention is not limited to the above-mentioned embodiments, and all technical solutions belonging to the idea of the present invention belong to the protection scope of the present invention. It should be noted that modifications and adaptations to those skilled in the art without departing from the principles of the present invention may be apparent to those skilled in the relevant art and are intended to be within the scope of the present invention.
Claims (7)
1. The method for analyzing the matching characteristic of the boiler heating surface working medium flow and the heat load distribution is characterized by comprising the following steps of: calculating a matching factor of the flow of the cooling working medium on the heating surface and the heat load, and grading the mismatching condition; calculating the influence of each influence factor on the matching factor, grading according to the influence degree, analyzing the mismatching degree of the heating surface cooling working medium flow and the heat load, and diagnosing the significant factor causing the mismatching of the heating surface cooling working medium flow and the heat load so as to determine the mismatching sequence of the heating surface working medium flow and the heat load distribution of the treatment boiler, wherein the method specifically comprises the following steps:
the method comprises the following steps: testing the wall temperature distribution of the outlet of the heating surface under different operating conditions;
step two: calculating the enthalpy rise distribution of the working medium of each heated pipe according to the operation parameters of the boiler and the wall temperature distribution;
step three: counting the distribution of enthalpy rise, and calculating average enthalpy rise, standard enthalpy rise deviation, maximum enthalpy rise deviation and relative maximum enthalpy rise deviation; abnormal enthalpy rise concentration ratio, maximum enthalpy rise migration and abnormal enthalpy rise distribution migration;
step four: calculating a matching factor of the flow of the cooling working medium on the heating surface and the heat load by adopting a statistical analysis method of a weighting factor; step five: grading the matching characteristics of the flow of the cooling working medium on the heating surface and the heat load according to the calculated matching factors;
step six: grading factors influencing the matching characteristic of the flow of the cooling working medium on the heating surface and the heat load;
step seven: calculating the influence of the factors on the matching factors of the flow of the cooling working medium on the heating surface and the heat load, and grading the influence degree of the factors;
step eight: grading the influence factors;
step nine: the treatment is carried out according to controllable obvious influence factors, controllable semi-obvious influence factors, semi-controllable obvious influence factors and semi-controllable semi-obvious influence factors when the working medium flow of the heating surface of the boiler is not matched with the heat load distribution, and the treatment is not considered for other factors.
2. The method for analyzing the matching characteristic of the working medium flow and the heat load distribution of the heating surface of the boiler according to claim 1, wherein the method comprises the following steps: calculating average enthalpy rise, standard enthalpy rise deviation, maximum enthalpy rise deviation and relative maximum enthalpy rise deviation; the formulas of the abnormal enthalpy rise concentration ratio, the maximum enthalpy rise migration and the abnormal enthalpy rise distribution migration are as follows:
mean enthalpy rise:
standard enthalpy rise deviation:
enthalpy rise deviation, i.e., the difference between the heat exchange tube enthalpy rise and the average enthalpy rise: Δ Hij＝Hij-HpkJ/kg (3)；
Maximum enthalpy rise deviation: Δ Hmax＝maxHij-HpkJ/kg (4)；
Relative maximum enthalpy rise deviation, i.e. maximum enthalpy rise deviation to average enthalpy rise ratio: h ismax＝Hmax/Hp (5)；
Standard relative enthalpy rise deviation, i.e., the ratio of heat exchange tube enthalpy rise deviation to standard enthalpy rise deviation: h isδ＝Hp/δH (6)；
Maximum standard enthalpy rise deviation, i.e. the ratio of maximum enthalpy rise deviation to standard enthalpy rise deviation: h is a total ofδmax＝Hmax/δH (7)；
Continuous abnormal enthalpy rise concentration ratio gammai: counting that standard enthalpy rise deviation of heat exchange tube is larger than 3 times deltaHThe interval between the two statistical points is not more than two, the statistical samples are regarded as continuous abnormal, and the number of the continuous abnormal points accounts for the percentage of the statistical samplesCalled continuous abnormal enthalpy rise concentration ratio, the formula is: gamma rayi＝100×∑Cont(h δ max＞3)/N％ (8)；
Maximum enthalpy rise migration τkn: under the working condition k and the working condition n, the number of statistical points of the position migration of the maximum enthalpy riser tube is called the maximum enthalpy riser migration in percentage of the statistical sample number, and the formula is as follows: tau iskn＝100×(S(Hmax)k-S(Hmax)n)/N％ (9)；
Wherein, S (Hmax)kThe position point of the water wall tube with the largest enthalpy rise under the operation working condition with the number of k is shown;
the central position Cen of the point of rise of the abnormal enthalpy, i.e. Cont (h)δmax>3) A symmetrical center position point of continuous abnormal enthalpy rise distribution;
abnormal enthalpy rise distribution migration τγ kn : counting the central positions of abnormal enthalpy rise points which are continuously distributed, wherein the number of counting points of the central position migration accounts for the percentage of the number of counting samples under the working condition k and the working condition n, and is called abnormal enthalpy rise distribution migration, and the formula is as follows: tau isγkn＝100×(Cenk(Cont(hδmax>3))-Cenn(Cont(hδmax>3)))/N％ (10)。
3. The method for analyzing the matching characteristic of the boiler heating surface working medium flow and the heat load distribution as claimed in claim 2, wherein the method comprises the following steps: fourthly, the matching factor is the matching performance of the flow of the cooling working medium on the heating surface and the heat load, and is expressed by the distribution of the enthalpy rise of the heat exchange tube and is marked as xi; the matching factor can be divided into a single working condition matching factor and a multi-working condition matching factor, and the multi-working condition matching factor xiht k,nThe functions of average enthalpy rise, standard enthalpy rise deviation, maximum enthalpy rise deviation, relative maximum enthalpy rise deviation, standard maximum enthalpy rise deviation, abnormal enthalpy rise concentration ratio, maximum enthalpy rise migration and abnormal enthalpy rise distribution migration are calculated by adopting a weighted statistical method; the single working condition matching factor xikThe matching characteristic under the single operation condition is only considered, and is a function of average enthalpy rise, standard enthalpy rise deviation, maximum enthalpy rise deviation, relative maximum enthalpy rise deviation, standard maximum enthalpy rise deviation and abnormal enthalpy rise concentration, and the sameCalculating by adopting a weighted statistical method, wherein a specific calculation formula is as follows:
ξht k＝f1(Hp,δH,Hmax,hmax,hδmax,γi)
＝α1×Hp+α2×δH+α3×Hmax+α4×hmax+α5×hδmax+α6×γi (11)；
ξht k,n＝f2(Hp,δH,Hmax,hmax,hδmax,γi,τkn,τγkn)＝β1×(Hp-Hp)+β2×(δk H-δn H)+β3×(Hk max-Hn max)+β4×(hk max-hn max)+β5×(hk δmax-hn δmax) +β6×(γk i-γn i)+β7×τkn+β8×τγkn (12)；
in the formula, alpha1、α2、α3、α4、α5、α6、β1、β2、β3、β4、β5、β6、β7And beta8Are all weighting coefficients.
4. The method for analyzing the matching characteristic of the working medium flow and the heat load distribution of the heating surface of the boiler as claimed in claim 3, wherein: step five, the grading method comprises the following steps: given standard matching factor xi01、ξ02When the calculated heating surface matching factor xiht>ξ01Time, significant mismatch; when xi02>ξht>ξ01When, it is a half significant mismatch; when xiht<ξ02When it is, it is a match.
5. The method for analyzing the matching characteristic of the boiler heating surface working medium flow and the heat load distribution as claimed in claim 4, wherein the method comprises the following steps: sixthly, the factors influencing the matching characteristics of the cooling working medium flow of the heating surface and the heat load are divided into factors influencing the flow distribution of the cooling working medium and factors influencing the heat load distribution, specifically, the factors influencing the flow distribution of the cooling working medium comprise: flow resistance coefficient of the heat exchange tube, orifice diameter, static pressure distribution of an inlet header, static pressure distribution of an outlet header, boiler output and foreign matter blockage conditions; factors that influence the distribution of the thermal load are: boiler output, hearth coking, heating surface ash accumulation, coal mill operation, primary air uniformity, secondary air door adjustment and coal quality of coal; the grading result is as follows: the uncontrollable factors are: the heat exchange pipe flow resistance coefficient, the inlet header static pressure distribution, the outlet header static pressure distribution and the coal quality; the controllable factors are as follows: the diameter of a throttling hole is used when the coal mill is put into use; the remaining factors are semi-controllable factors.
6. The method for analyzing the matching characteristic of the working medium flow and the heat load distribution of the heating surface of the boiler as claimed in claim 5, wherein: the seventh step is that the method for grading the influence degree of the factors comprises the following steps: calculating the matching factor xi of the heating surface cooling working medium flow and the heat load of the controllable factor and the semi-controllable factor through different working conditions influencing the change of the factorsht k,nMatching the criterion factor xief01、ξef02Comparing, grading the influence of the factors on the matching factors, xiht k,n>ξef01Time, is a significant influencing factor; when xief02>ξht k,n>ξef01Time, is a semi-significant influencing factor; when xiht k,n<ξef02Time is an insignificant influence.
7. The method for analyzing the matching characteristic of the working medium flow and the heat load distribution of the heating surface of the boiler according to any one of claims 1 to 6, wherein the method comprises the following steps: step eight, the method for grading the influence factors comprises the following steps: according to the controllable characteristics of the influence factors and the influence degree on the matching factors, the influence factors are divided into: controllable significant influence factors, semi-controllable significant influence factors and non-controllable significant influence factors; controllable semi-significant influencing factors, semi-controllable semi-significant influencing factors and non-controllable semi-significant influencing factors; and (4) no significant influence factor.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811388758.1A CN109284581B (en) | 2018-11-20 | 2018-11-20 | Method for analyzing matching characteristic of boiler heating surface working medium flow and heat load distribution |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811388758.1A CN109284581B (en) | 2018-11-20 | 2018-11-20 | Method for analyzing matching characteristic of boiler heating surface working medium flow and heat load distribution |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109284581A CN109284581A (en) | 2019-01-29 |
| CN109284581B true CN109284581B (en) | 2022-06-03 |
Family
ID=65175874
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811388758.1A Active CN109284581B (en) | 2018-11-20 | 2018-11-20 | Method for analyzing matching characteristic of boiler heating surface working medium flow and heat load distribution |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109284581B (en) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111765798B (en) * | 2019-03-31 | 2022-03-11 | 上海梅山钢铁股份有限公司 | Waste heat boiler heat pipe exchanger online intelligent diagnosis and bypass protection device and control method |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101493861A (en) * | 2009-01-15 | 2009-07-29 | 哈尔滨工业大学 | Hydrodynamic numerical implementation for natural circulation boiler |
| CN105202523A (en) * | 2015-10-28 | 2015-12-30 | 贵州电网有限责任公司电力科学研究院 | Method for calculating W flame supercritical power station boiler sewage cooling wall flue gas thermal load distribution online |
| CN108105749A (en) * | 2017-11-30 | 2018-06-01 | 湖南省湘电试验研究院有限公司 | Working medium flow On-line Measuring Method and system in a kind of water screen tube |
-
2018
- 2018-11-20 CN CN201811388758.1A patent/CN109284581B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101493861A (en) * | 2009-01-15 | 2009-07-29 | 哈尔滨工业大学 | Hydrodynamic numerical implementation for natural circulation boiler |
| CN105202523A (en) * | 2015-10-28 | 2015-12-30 | 贵州电网有限责任公司电力科学研究院 | Method for calculating W flame supercritical power station boiler sewage cooling wall flue gas thermal load distribution online |
| CN108105749A (en) * | 2017-11-30 | 2018-06-01 | 湖南省湘电试验研究院有限公司 | Working medium flow On-line Measuring Method and system in a kind of water screen tube |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109284581A (en) | 2019-01-29 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US5181482A (en) | Sootblowing advisor and automation system | |
| CN105160158B (en) | The processing method of boiler secondary air baffle characteristics test data | |
| CN109284581B (en) | Method for analyzing matching characteristic of boiler heating surface working medium flow and heat load distribution | |
| CN109992825B (en) | Boiler four-tube real-time service life assessment method considering wall thickness reduction and overheating influence | |
| WO2021208343A1 (en) | Neural network model-based method and apparatus for water spray adjustment of water spray desuperheater | |
| CN106773955A (en) | A kind of boiler soot-blowing optimizes system and its optimization method | |
| CN107525277A (en) | A kind of low nitrogen wall-hung condensing gas furnace and its application process | |
| WO2023279601A1 (en) | Intelligent soot blowing control method based on multi-dimensional evaluation factors, and system and storage medium | |
| CN110846490B (en) | Optimization calculation method for postweld heat treatment heating rate of 9% Cr hot-strength steel pipeline | |
| CN112325267B (en) | Pressure-stabilizing steam blowing pipe system and method for double reheat boiler | |
| CN111242279B (en) | Slag bonding prediction system and method for ultra-supercritical pulverized coal boiler furnace | |
| CN107143840B (en) | A kind of steam boiler system of intelligent control blowdown reference data | |
| CN107101191B (en) | A kind of cloud computing boiler blow-out system | |
| CN103244942B (en) | Utility boiler heating surface design method applicable to flue-gas temperature deviation fixation | |
| CN107084382B (en) | A kind of component feeds back automatic control boiler drainage | |
| LU500112B1 (en) | Cooling system of a metallurgical furnace | |
| CN112664965B (en) | Coal-fired boiler soot blowing operation optimization method based on experimental quantitative analysis | |
| CN104566413B (en) | A kind of method of fast selecting boiler blow piping parameter | |
| CN107246604B (en) | Method and system for adjusting reheat steam temperature of boiler | |
| CN112797398A (en) | Ultra-supercritical secondary reheating unit boiler temperature regulating system | |
| CN111304393A (en) | Method for checking main technological parameters of plate heat exchanger | |
| CN205208602U (en) | Modular intelligence phase transition heat transfer device | |
| CN107557547A (en) | A kind of post weld heat treatment technique of austenitic stainless steel containing niobium and fluid bed | |
| CN115161445B (en) | Method for optimizing medium-frequency induction heating local postweld heat treatment parameters of 9% Cr hot-strength steel pipeline | |
| CN114413251B (en) | System for adjusting main steam temperature of Pi-type opposite-flow direct-current boiler |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |