CN111174569A - Method and system for online prediction of flue gas temperature of calcining section in rotary kiln - Google Patents

Method and system for online prediction of flue gas temperature of calcining section in rotary kiln Download PDF

Info

Publication number
CN111174569A
CN111174569A CN202010046835.6A CN202010046835A CN111174569A CN 111174569 A CN111174569 A CN 111174569A CN 202010046835 A CN202010046835 A CN 202010046835A CN 111174569 A CN111174569 A CN 111174569A
Authority
CN
China
Prior art keywords
flue gas
temperature
kiln
rotary kiln
secondary air
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202010046835.6A
Other languages
Chinese (zh)
Other versions
CN111174569B (en
Inventor
鲁聪
鄢曙光
欧阳安妮
宋紫欣
陈嘉仪
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Wuhan University of Science and Engineering WUSE
Original Assignee
Wuhan University of Science and Engineering WUSE
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Wuhan University of Science and Engineering WUSE filed Critical Wuhan University of Science and Engineering WUSE
Priority to CN202010046835.6A priority Critical patent/CN111174569B/en
Publication of CN111174569A publication Critical patent/CN111174569A/en
Application granted granted Critical
Publication of CN111174569B publication Critical patent/CN111174569B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F27FURNACES; KILNS; OVENS; RETORTS
    • F27BFURNACES, KILNS, OVENS, OR RETORTS IN GENERAL; OPEN SINTERING OR LIKE APPARATUS
    • F27B7/00Rotary-drum furnaces, i.e. horizontal or slightly inclined
    • F27B7/20Details, accessories, or equipment peculiar to rotary-drum furnaces
    • F27B7/42Arrangement of controlling, monitoring, alarm or like devices

Landscapes

  • Engineering & Computer Science (AREA)
  • Mechanical Engineering (AREA)
  • General Engineering & Computer Science (AREA)
  • Muffle Furnaces And Rotary Kilns (AREA)

Abstract

The invention relates to a method and a system for on-line prediction of flue gas temperature of a calcining section in a rotary kiln, which specifically comprise the following steps: construction of Q in rotary kilnsSecondary air、QBurning device、QFlue gas、QWall(s)And QMaterial(s)The heat balance model of (1); calculating Q in the heat balance model according to the parameters of the rotary kiln in stable operationWall(s)、QMaterial(s)And QSecondary airThe value of (d); q to be obtainedWall(s)、QMaterial(s)And QSecondary airSubstituting the value of (A) into the rotary kiln heat balance model in the step S1 to obtain QBurning deviceAnd QFlue gasThe relation between the two, and then the flue gas temperature t of the calcining section is obtainedrTemperature t of kiln tail flue gasyThe relational expression of (1); measuring kiln tail flue gas temperature tyAnd passing through the flue gas temperature t of the calcination sectionrTemperature t of kiln tail flue gasyObtaining the flue gas temperature t of the calcining section by the relational expressionr(ii) a The heat balance model specifically comprises: qSecondary air+QBurning device=QFlue gas+QWall(s)+QMaterial(s). The method has the advantages of effectively predicting the temperature of the combustion section of the rotary kiln and providing theoretical support for the production of the rotary kiln.

Description

Method and system for online prediction of flue gas temperature of calcining section in rotary kiln
Technical Field
The invention relates to the technical field of rotary kilns, in particular to a method and a system for predicting flue gas temperature of a calcining section in a rotary kiln on line.
Background
The rotary kiln is a key device widely applied to production links of cement, metallurgy and the like, and the control of the internal temperature of the rotary kiln is directly related to the stable operation, the product quality and the production cost of the whole production line. The heat released by the fuel combustion can generate a high temperature zone, generally about 1350 ℃ for the pellet rotary kiln under stable work, the temperature is a very key process parameter in production, and is directly related to the fuel consumption, the pellet quality, the emission of nitrogen oxides and barrel ring formation. However, the high-temperature flue gas temperature at the calcining section of the rotary kiln cannot be directly measured for three reasons:
1) the rotary kiln is high-temperature equipment, and the temperature of a calcining section of the lime kiln is about 1350 ℃; the temperature of the calcining section of the cement kiln is even higher than 1500 ℃; the temperature of the calcining section of the pellet rotary kiln is generally higher than 1200 ℃ and does not exceed 1400 ℃;
2) the rotary kiln is in a rotating state when working, the calcining section is positioned in the kiln, and a measuring tool is difficult to extend into the kiln;
3) the fluid in the rotary kiln is in a fast flowing state, generally about 20m/s, and even more than 50m/s at the nozzle, so that the measurement difficulty is high.
Disclosure of Invention
The invention aims to solve the technical problem of providing a method and a system for on-line prediction of flue gas temperature of a calcining section in a rotary kiln, and aims to solve the technical problem.
The technical scheme for solving the technical problems is as follows:
a method for on-line prediction of flue gas temperature of a calcining section in a rotary kiln specifically comprises the following steps:
s1, constructing the Q in the rotary kiln according to the heat balance principle when the rotary kiln stably worksSecondary air、QBurning device、QFlue gas、QWall(s)And QMaterial(s)The heat balance model of (1);
s2, calculating Q in the heat balance model according to the parameters of the rotary kiln in stable operationWall(s)、QMaterial(s)And QSecondary airThe value of (d);
s3, mixing the Q obtained in S2Wall(s)、QMaterial(s)And QSecondary airSubstituting the value of (A) into the rotary kiln heat balance model in the step S1 to obtain QBurning deviceAnd QFlue gasThe relation between the two, and then the flue gas temperature t of the calcining section is obtainedrTemperature t of kiln tail flue gasyThe relational expression of (1);
s4, measuring the temperature t of the kiln tail flue gasyAnd passing through the flue gas temperature t of the calcination sectionrTemperature t of kiln tail flue gasyObtaining the flue gas temperature t of the calcining section by the relational expressionr
The heat balance model specifically comprises: qSecondary air+QBurning device=QFlue gas+QWall(s)+QMaterial(s)
QSecondary airThe energy of the high-temperature combustion-supporting air which directly enters the rotary kiln from the ring cooling section in unit time is expressed;
Qburning deviceRepresents all heat released after the fuel sprayed by the nozzle is completely combusted in unit time;
Qflue gasThe energy of the kiln tail outlet flue gas increased relative to the normal temperature flue gas temperature in unit time is represented;
Qwall(s)The heat quantity dissipated by the wall surface of the kiln in unit time is represented;
Qmaterial(s)Means the increased energy per unit time of pellets exiting the rotary kiln after the completion of the calcination section relative to the energy entering the kiln.
The invention has the beneficial effects that: the invention provides an on-line prediction method for the flue gas temperature of a calcination section of a rotary kiln, which solves the problem that the flue gas temperature of the calcination section cannot be directly measured when the rotary kiln stably works.
On the basis of the technical scheme, the invention can be further improved as follows.
Further, the specific step of S2 includes:
s21: uniformly dividing the kiln wall of the combustion section of the rotary kiln into n sections;
s22: measuring the temperature t of each sectioniAnd surface area SiAnd the ambient temperature t corresponding to the kiln wall of the combustion section of the rotary kiln
S23: measuring the temperature t of each section according to the step S22iAnd surface area SiAnd the ambient temperature t corresponding to the kiln wall of the combustion section of the rotary kilnObtaining the heat dissipation Q of the kiln wall facing the external environment in unit time when the rotary kiln stably worksWall(s)
The beneficial effect of adopting the above further scheme is that the kiln wall of the rotary kiln is equally divided into a plurality of sections, the temperature and the surface area of each section are measured, the temperature is usually measured by an infrared thermometer, so that the heat dissipation capacity of each section is calculated, the whole rotary kiln is further calculated, and the accuracy is high.
Further, in the S23, QWall(s)The calculation formula of (a) is as follows:
Figure BDA0002369728350000031
in the formula, h represents the convective heat transfer coefficient between the wall surface of the kiln and the external environment.
The beneficial effect of adopting above-mentioned further scheme is that through the summation of integrals calculate the heat dissipation capacity of whole rotary kiln wall when the rotary kiln steady operation, the accuracy is high.
Further, the specific step of S3 includes:
s31: measuring the mass m of finished pellets flowing out of a kiln head in unit time and the temperature t of finished pellets discharged from the kiln when the rotary kiln stably worksoutInitial temperature t of green pellets put into the kiln in unit timein
S32: according to the mass m of finished pellets flowing out of the kiln head in unit time and the temperature t of finished pellets discharged from the kiln in the step S31outAnd the initial temperature t of green pellets entering the kiln per unit timeinObtaining the energy Q increased by the pellets before and after entering the kilnMaterial(s)
The beneficial effect of adopting the above further scheme is that the quality of finished pellet ore flowing out of the kiln head in unit time is measured in a manner imaginable to those skilled in the art when the rotary kiln works stably, and the temperature of finished pellets discharged from the kiln and the initial temperature of green pellets entering the kiln in unit time are measured by the infrared thermometer, so that the measurement is convenient and fast.
Further, in the S32, QMaterial(s)The calculation formula of (a) is as follows:
Qmaterial(s)=m×cm×(tout-tin),
In the formula, cmRepresents the specific heat capacity of the pellets.
The beneficial effect of adopting the above further scheme is that the heat of the material is calculated according to each measured parameter and the heat formula.
Further, the specific steps of the steps include:
s41: measuring primary air temperature t sprayed from kiln head in unit time when rotary kiln stably works1Volume of secondary air V2And secondary air temperature t2
S42: according to the measured primary air temperature t sprayed from the kiln head in the unit time measured in the step S421Volume of secondary air V2And secondary air temperature t2Obtaining the heat Q brought by the secondary air in unit timeSecondary air
The beneficial effect of adopting the above further scheme is that the volume of the secondary air is calculated by the way that can be thought by the technical personnel in the field, and the primary air temperature and the secondary air temperature are measured by the infrared thermometer, so that the measurement is convenient and fast.
Further, in the S42, QSecondary airThe calculation formula of (a) is as follows:
Qsecondary air=V2×ρ2×c2×(t2-t1),
In the formula, ρ2Representing secondary air gas density, c2Represents the specific heat capacity of the secondary air.
The beneficial effect of adopting the further scheme is that the heat of the secondary air is calculated according to the measured parameters and the heat formula.
Further, QBurning deviceAnd QFlue gasThe relation between is (V)y-V2)×ρr×cr×(tr-t1)=Vy×ρy×cy×(ty-t1)-Q,
In the formula, VyIs the volume (m) of the outlet flue gas of the rotary kiln in unit time3);ρyIs the density of the outlet flue gas; c. CyIs the specific heat capacity of the outlet flue gas; t is tyIs the outlet flue gas temperature; rhorThe density of the high-temperature flue gas at the calcining section; c. CrThe specific heat capacity of the high-temperature flue gas at the calcining section; t is trIs the high temperature flue gas temperature of the calcining section.
The beneficial effect of adopting the further scheme is that Q is constant, then the heat of the fuel and the flue gas is expressed by a heat formula, a new formula is obtained, and the derivation is convenient and fast.
Further, the step of simulating the calcination section flue gas temperature t obtained in the step S3 by using a simulation is further included between the step S3 and the step S4rTemperature t of kiln tail flue gasyThe relational expression (2) is verified, and if the verification is passed, the step S4 is executed, and if the verification is not passed, the step S2 is returned to be executed.
The method has the beneficial effect that the obtained calcining section flue gas temperature t is subjected to simulationrTemperature t of kiln tail flue gasyIf the verification is passed, the step S4 is executed, and if the verification is not passed, the step S2 is executed again, so that the accuracy is greatly improved.
A system for on-line prediction of flue gas temperature of a calcining section in a rotary kiln comprises the following modules,
a heat balance model building module for building Q-value in the rotary kiln according to the heat balance principle when the rotary kiln stably worksSecondary air、QBurning device、QFlue gas、QWall(s)And QMaterial(s)The heat balance model of (1);
a model parameter calculation module for calculating Q in the heat balance model according to the parameters of the rotary kiln in stable operationWall(s)、QMaterial(s)And QSecondary airThe value of (d);
a relational expression obtaining module for obtaining QWall(s)、QMaterial(s)And QSecondary airSubstituting the numerical value into a rotary kiln heat balance model to obtain QBurning deviceAnd QFlue gasThe relation between the two, and then the flue gas temperature t of the calcining section is obtainedrTemperature t of kiln tail flue gasyThe relational expression of (1);
a calculation module for measuring the kiln tail flue gas temperature tyAnd passing through the flue gas temperature t of the calcination sectionrTemperature t of kiln tail flue gasyObtaining the flue gas temperature t of the calcining section by the relational expressionr
The heat balance model specifically comprises: qSecondary air+QBurning device=QFlue gas+QWall(s)+QMaterial(s);QSecondary airThe energy of the high-temperature combustion-supporting air which directly enters the rotary kiln from the ring cooling section in unit time is expressed;
Qburning deviceRepresents all heat released after the fuel sprayed by the nozzle is completely combusted in unit time;
Qflue gasThe energy of the kiln tail outlet flue gas increased relative to the normal temperature flue gas temperature in unit time is represented;
Qwall(s)The heat quantity dissipated by the wall surface of the kiln in unit time is represented;
Qmaterial(s)Means the increased energy per unit time of pellets exiting the rotary kiln after the completion of the calcination section relative to the energy entering the kiln.
The beneficial effect of adopting above-mentioned further scheme is that the system of calcining section flue gas temperature in rotary kiln is provided to the prediction on line, can be according to the rotary kiln burning section heat balance formula fast prediction rotary kiln burning section flue gas temperature, convenient and fast.
Drawings
FIG. 1 is a flow chart of the method for on-line prediction of flue gas temperature at the calcining section in a rotary kiln according to the present invention;
FIG. 2 is a schematic view of the construction of a rotary kiln according to the present invention;
FIG. 3 shows the kiln tail flue gas temperature t in the present inventionyAnd the temperature t of the flue gas in the calcination sectionrThe relationship toggles a trend line.
Detailed Description
The principles and features of this invention are described in connection with the drawings and the detailed description of the invention, which are set forth below as examples to illustrate the invention and not to limit the scope of the invention.
As shown in fig. 1 to 3, the invention provides a method for online predicting flue gas temperature of a calcining section in a rotary kiln, which specifically comprises the following steps: s1, constructing the Q in the rotary kiln according to the heat balance principle when the rotary kiln stably worksSecondary air、QBurning device、QFlue gas、QWall(s)And QMaterial(s)The heat balance model of (1);
s2, calculating Q in the heat balance model according to the parameters of the rotary kiln in stable operationWall(s)、QMaterial(s)And QSecondary airThe value of (d);
s3, mixing the Q obtained in S2Wall(s)、QMaterial(s)And QSecondary airSubstituting the value of (A) into the rotary kiln heat balance model in the step S1 to obtain QBurning deviceAnd QFlue gasThe relation between the two, and then the flue gas temperature t of the calcining section is obtainedrTemperature t of kiln tail flue gasyThe relational expression of (1);
s4, measuring the temperature t of the kiln tail flue gasyAnd passing through the flue gas temperature t of the calcination sectionrTemperature t of kiln tail flue gasyObtaining the flue gas temperature t of the calcining section by the relational expressionr
The heat balance model specifically comprises: qSecondary air+QBurning device=QFlue gas+QWall(s)+QMaterial(s);QSecondary airThe energy of the high-temperature combustion-supporting air which directly enters the rotary kiln from the ring cooling section in unit time is expressed;
Qburning deviceRepresents all heat released after the fuel sprayed by the nozzle is completely combusted in unit time;
Qflue gasThe energy of the kiln tail outlet flue gas increased relative to the normal temperature flue gas temperature in unit time is represented;
Qwall(s)The heat quantity dissipated by the wall surface of the kiln in unit time is represented;
Qmaterial(s)Means the increased energy per unit time of pellets exiting the rotary kiln after the completion of the calcination section relative to the energy entering the kiln.
The invention provides an on-line prediction method for the flue gas temperature of a calcination section of a rotary kiln, which solves the problem that the flue gas temperature of the calcination section cannot be directly measured when the rotary kiln stably works.
In addition, the temperature of the rotary kiln is 50-200 ℃ when the rotary kiln works stably.
In the present invention, the specific step of S2 includes:
s21: uniformly dividing the kiln wall of the combustion section of the rotary kiln into n sections;
s22: measuring the temperature t of each sectioni(° c) and surface area Si(m2) And the ambient temperature t corresponding to the kiln wall of the combustion section of the rotary kiln(℃);
S23: measuring the temperature t of each section according to the step S22i(° c) and surface area Si(m2) And the ambient temperature t corresponding to the kiln wall of the combustion section of the rotary kilnThe temperature is measured to obtain the heat dissipation Q of the kiln wall facing the external environment in unit time when the rotary kiln stably worksWall(s)
In S23, QWall(s)The calculation formula of (a) is as follows:
Figure BDA0002369728350000071
in the formula, h represents the convective heat transfer coefficient between the wall surface of the kiln and the external environment.
The kiln wall of the rotary kiln is divided into a plurality of sections equally, the temperature and the surface area of each section are measured, the temperature is usually measured by an infrared thermometer, and therefore the heat dissipation capacity of each section is calculated, the whole rotary kiln is further calculated, and the accuracy is high; the heat dissipation capacity of the kiln wall of the whole rotary kiln is calculated through integral summation when the rotary kiln stably works, and the accuracy is high.
In the present invention, the specific step of S3 includes:
s31: measuring the mass m (kg) of finished pellets flowing out of a kiln head in unit time and the temperature t of finished pellets discharged from the kiln when the rotary kiln stably worksout(DEG C) and the initial temperature t of green pellets entering the kiln in unit timein(℃);
S32: according to the mass m (kg) of the finished pellets flowing out of the kiln head in unit time and the temperature t of the finished pellets discharged from the kiln measured in the step S31out(DEG C) and the initial temperature t of green pellets entering the kiln in unit timein(DEG C) to obtain the energy Q increased by the pellets before and after entering the kilnMaterial(s)
In said S32, QMaterial(s)The calculation formula of (a) is as follows:
Qmaterial(s)=m×cm×(tout-tin),
In the formula, cmRepresents the specific heat capacity of the pellet [ J/(kg. K)]。
The quality of finished pellets flowing out of a kiln head in unit time when the rotary kiln stably works is measured in a mode which can be thought by a person skilled in the art, and meanwhile, the temperature of finished pellets discharged from the kiln and the initial temperature of green pellets entering the kiln in unit time are measured by an infrared thermometer, so that the measurement is convenient and quick; and calculating the heat of the material according to the measured parameters and a heat formula.
In the present invention, the specific step of S4 includes:
s41: measuring primary air temperature t sprayed from kiln head in unit time when rotary kiln stably works1(° c), secondary air volume V2(m3) And secondary air temperature t2(℃);
S42: according to the measured primary air temperature t sprayed from the kiln head in the unit time measured in the step S421(° c), secondary air volume V2(m3) And secondary air temperature t2The temperature is higher than the temperature of the secondary air in unit time to obtain the heat Q brought by the secondary air in unit timeSecondary air
In said S42, QSecondary airThe calculation formula of (a) is as follows:
Qsecondary air=V2×ρ2×c2×(t2-t1),
In the formula, ρ2Represents the secondary air gas density (kg/m)3),c2Represents the specific heat capacity of secondary air gas [ J/(kg. K)]。
The volume of the secondary air is calculated in a manner which can be thought by a person skilled in the art, and the primary air temperature and the secondary air temperature are measured by an infrared thermometer, so that the measurement is convenient and quick; and calculating the heat of the secondary air according to the measured parameters and a heat formula.
In the invention, the heat balance formula of the rotary kiln in S1 can be obtained, wherein Q is QWall(s)+QMaterial(s)-QSecondary air,QBurning device-Q=QFlue gas
By converting the formula in step S1, it can be known that the heat dissipated from the kiln wall of the rotary kiln, the heat of the material and the heat of the secondary air can be measured, and therefore the difference between the heat of the fuel and the heat of the flue gas is the theorem Q.
Derivation of Q from the above equationBurning deviceAnd QFlue gasThe relation between is
(Vy-V2)×ρr×cr×(tr-t1)=Vy×ρy×cy×(ty-t1)-Q,
In the formula, VyIs the volume (m) of the outlet flue gas of the rotary kiln in unit time3);ρyThe density of the outlet flue gas (kg/m 3); c. CySpecific heat capacity of outlet flue gas [ J/(kg. K)];tyThe temperature (DEG C) of the outlet flue gas is shown; rhorIs the density (kg/m) of high-temperature flue gas in a calcining section3);crIs the specific heat capacity of high-temperature flue gas at the calcining stage [ J/(kg. K)];trThe air temperature (DEG C) of the high-temperature flue gas at the calcining section is shown.
And Q is a constant, and then the heat of the fuel and the smoke is expressed by a heat formula to obtain a new formula, so that the derivation is convenient and quick.
In the invention, the flue gas temperature t of the calcination sectionrTemperature t of kiln tail flue gasyHas the relation of
Figure BDA0002369728350000091
The flue gas temperature t of the calcining section is obtained through the change of a formularTemperature t of kiln tail flue gasyThe derivation is convenient and fast.
In the invention, the step of simulating the calcination section flue gas temperature t obtained in the step S3 by using a simulation is further included between the step S3 and the step S4rTemperature t of kiln tail flue gasyThe relational expression (2) is verified, and if the verification is passed, the step S4 is executed, and if the verification is not passed, the step S2 is returned to be executed.
The simulation usually adopts Ansys software, which is large-scale general Finite Element Analysis (FEA) software developed by Ansys corporation in the united states, is Computer Aided Engineering (CAE) software growing fastest worldwide, and can interface with most Computer Aided Design (CAD) software to realize data sharing and exchange.
In the present invention, the heat balance analysis of the rotary kiln is as follows:
for the pellet rotary kiln, the heat balance in the steady operation state can be expressed by the simplified relation in step S1, which is simplified by two points:
1) for gas, the energy brought by normal temperature air (i.e. primary air) entering the rotary kiln is taken as a standard, namely the energy is considered to be 0 relatively;
2) for the material, the temperature of the pellet from the grate (feeding the pellet into the rotary kiln) is taken as a standard, that is, the energy brought by the pellet entering the rotary kiln is relatively 0;
Qsecondary air+QBurning device=QFlue gas+QWall(s)+QMaterial(s)
In the formula, QSecondary airIndicating that high-temperature combustion-supporting air directly entering the rotary kiln from a ring cooling section is carried by the high-temperature combustion-supporting air in unit timeThe energy of (a);
Qburning deviceRepresents all heat released after the fuel sprayed by the nozzle is completely combusted in unit time;
Qflue gasThe energy of the kiln tail outlet flue gas increased relative to the normal temperature flue gas temperature in unit time is represented;
Qwall(s)The heat quantity dissipated from the wall surface of the kiln to the outside in unit time is represented;
Qmaterial(s)Means the increased energy per unit time of pellets exiting the rotary kiln after completion of the calcination section relative to their entry into the kiln.
For a rotary kiln operating steadily, QSecondary air、QBurning device、QFlue gas、QWall(s)、QMaterial(s)All the components are constant, however, according to actual production conditions, in the stable working process of the rotary kiln, the internal temperature field is not absolutely stable, certain fluctuation exists, manual regulation is sometimes needed, the temperature of the flue gas at the calcining section is transient, in this case, the heat balance is broken, and Q is obtainedBurning deviceAt a small instantaneous change, the temperature t of the calcination sectionrAlso can change to a certain extent, because the flue gas flow velocity in the kiln is larger, the kiln tail flue gas tyAn immediate corresponding change is immediately followed. At this time, QBurning device、QFlue gasIs independent variable and dependent variable; qSecondary airIs still constant; because the temperature of the outer wall surface of the kiln is influenced by the temperature of the inner wall surface of the kiln through heat conduction, the heat conductivity coefficient is very low due to the obstruction of refractory materials, and the reaction of the outer wall of the kiln to the temperature in the kiln needs a period of time, therefore QWall(s)Is also a constant; the pellet at the outlet is completely calcined, so that the temperature can not obviously change due to internal chemical reaction, and the temperature change caused by heat conduction due to the temperature change of the wall surface of the kiln can not be immediately receivedMaterial(s)Is also a constant. Therefore, the temperature of the molten metal is controlled,the above relationship can be converted to the following relationship:
Qburning device=QFlue gas+Q
In the formula: q ═ QWall(s)+QMaterial(s)-QSecondary airAre constant and can be further calculated by measurable data,
Figure BDA0002369728350000111
in the formula, h represents the convective heat transfer coefficient between the wall surface of the kiln and the external environment; si represents the wall surface area (m) of the i-th stage2);ti、tRespectively representing the temperature of the kiln wall surface of the i-th section and the ambient temperature (DEG C); because the wall temperature of the kiln is not uniform, the temperatures of different sections are different, the temperatures of different areas need to be measured for many times and are obtained by utilizing formula integration,
Qmaterial(s)=m×cm×(tout-tin)
Wherein m represents the mass (kg) of the outlet pellet per unit time; c. CmRepresents the specific heat capacity of the pellet [ J/(kg. K)];tout、tinRespectively showing the temperature (DEG C) of the pellets at the outlet and the pellets at the inlet of the rotary kiln;
Qsecondary air=V2×ρ2×c2×(t2-t1)
In the formula (I); v2Represents the volume (m) of secondary air entering the rotary kiln in unit time3);ρ2Represents the secondary air gas density (kg/m)3);c2Represents the specific heat capacity of secondary air gas [ J/(kg. K)];t2、t1Respectively representing the temperature (DEG C) of secondary air and primary air;
the heat released by the fuel combustion can generate a high temperature zone, generally about 1350 ℃ for the pellet rotary kiln under stable work, the temperature is a very key process parameter in production, and is directly related to the fuel consumption, the pellet quality, the emission of nitrogen oxides and barrel ring formation. The temperature of the area is very high and is positioned in the kiln cylinder, the specific position is difficult to determine, and the velocity of the flow field in the kiln is also high, so that the measurement is difficult. Therefore, the study on the temperature field problem in the kiln is very important! The above contents also show that the kiln tail flue gas temperature and the temperature of the high-temperature area in the kiln have a definite relation, and the kiln tail flue gas temperature is relatively convenient to measure in production, so that the finding of the relation between the kiln tail flue gas temperature and the temperature of the high-temperature area in the kiln is significant.
In addition, QBurning deviceThe rotary kiln has several types of fuels for fuel combustion, the common types of the fuels comprise coal powder, natural gas or liquefied natural gas and the like, and the main combustion component of the rotary kiln is solid C or gaseous CH4The general combustion equation is as follows:
C(s)+O2(g)→CO2(g)
CH4(g)+2O2(g)→CO2(g)+2H2O(g)
the two relationships are the complete overall reaction of fuel combustion, and it can be seen that either solid C combustion or gaseous CH4And (3) combustion, wherein the molecular weight of the gas state before and after the reaction is unchanged, so that the volume of the gas before and after the combustion is not increased. The CO formation is not considered because the oxygen in the rotary kiln is in excess, and even if there is a trace amount of CO at the exit, it is generally not more than 600ppm and can be ignored. NO consideration is given to NO production because the content is also very low, generally not more than 800ppm, which is negligible, and NO production is mainly of fuel type and thermodynamic type, the mechanism for producing NO of fuel type is too complicated, which is not yet clear, and the molecular weight of NO of thermodynamic type is not changed before and after the reaction. Based on this, there are:
Qflue gas=Vy×ρy×cy×(ty-t1)
In the formula, VyThe volume (m) of the outlet flue gas of the rotary kiln in unit time3);ρyDenotes the density (kg/m) of the outlet flue gas3);cyRepresents the specific heat capacity of the outlet flue gas [ J/(kg. K)];ty、t1Respectively showing the temperature of outlet flue gas and the temperature of primary air (DEG C);
Qburning device=(Vy-V2)×ρr×cr×(tr-t1)
In the formula, ρrShows the density (kg/m) of high-temperature flue gas immediately after the completion of heat release by combustion of fuel3) Namely the density of high-temperature flue gas in the calcining section; c. CrRepresents the specific heat capacity of the high-temperature flue gas in the calcining section [ J/(kg. K)];tr、t1Respectively showing the high-temperature flue gas air temperature and the primary air temperature (DEG C) of the calcining section;
in summary, the following relationships are further derived:
(Vy-V2)×ρr×cr×(tr-t1)=Vy×ρy×cy×(ty-t1)-Q
Figure BDA0002369728350000121
because the secondary air entering the rotary kiln, the flue gas after the fuel is completely combusted and the flue gas at the tail of the rotary kiln are high-temperature gases (the temperature is higher than 1100K), and N is mainly contained2、CO2、O2And a trace amount of NOxCO and the like, the main components are basically the same, and the density rho and the specific heat capacity c can be considered to be the same, so that the relational expression can be finally simplified to obtain tyAnd trThe relation of (1):
Figure BDA0002369728350000122
in the formula, Vy、V2、Q、t1Constant, rho, which can be measured when the rotary kiln works stablyy、cyThe density and specific heat capacity of high-temperature flue gas at the outlet of the rotary kiln are the properties of the substance and are constant. The above relation expresses an important meaning: when the rotary kiln works stably, the kiln tail flue gas temperature and the high-temperature flue gas temperature fluctuation of the calcining section have a linear relation, and the linear relation is a linear function relation, and the slope of the linear relation is only related to the kiln tail flue gas flow and the secondary air inlet flow.
The invention also provides a system for on-line prediction of the flue gas temperature of the calcining section in the rotary kiln, which is characterized by comprising the following modules,
a heat balance model building module for building Q-value in the rotary kiln according to the heat balance principle when the rotary kiln stably worksSecondary air、QBurning device、QFlue gas、QWall(s)And QMaterial(s)The heat balance model of (1);
a model parameter calculation module for calculating Q in the heat balance model according to the parameters of the rotary kiln in stable operationWall(s)、QMaterial(s)And QSecondary airThe value of (d);
a relational expression obtaining module for obtaining QWall(s)、QMaterial(s)And QSecondary airSubstituting the numerical value into a rotary kiln heat balance model to obtain QBurning deviceAnd QFlue gasThe relation between the two, and then the flue gas temperature t of the calcining section is obtainedrTemperature t of kiln tail flue gasyThe relational expression of (1);
a calculation module for measuring the kiln tail flue gas temperature tyAnd passing through the flue gas temperature t of the calcination sectionrTemperature t of kiln tail flue gasyObtaining the flue gas temperature t of the calcining section by the relational expressionr
The heat balance model specifically comprises: qSecondary air+QBurning device=QFlue gas+QWall(s)+QMaterial(s);QSecondary airThe energy of the high-temperature combustion-supporting air which directly enters the rotary kiln from the ring cooling section in unit time is expressed;
Qburning deviceRepresents all heat released after the fuel sprayed by the nozzle is completely combusted in unit time;
Qflue gasThe energy of the kiln tail outlet flue gas increased relative to the normal temperature flue gas temperature in unit time is represented;
Qwall(s)The heat quantity dissipated by the wall surface of the kiln in unit time is represented;
Qmaterial(s)Means the increased energy per unit time of pellets exiting the rotary kiln after the completion of the calcination section relative to the energy entering the kiln.
The above description is only for the purpose of illustrating the preferred embodiments of the present invention and is not to be construed as limiting the invention, and any modifications, equivalents, improvements and the like that fall within the spirit and principle of the present invention are intended to be included therein.

Claims (10)

1. A method for on-line prediction of flue gas temperature of a calcining section in a rotary kiln is characterized by comprising the following steps:
s1, constructing the Q in the rotary kiln according to the heat balance principle when the rotary kiln stably worksSecondary air、QBurning device、QFlue gas、QWall(s)And QMaterial(s)The heat balance model of (1);
s2, calculating Q in the heat balance model according to the parameters of the rotary kiln in stable operationWall(s)、QMaterial(s)And QSecondary airThe value of (d);
s3, mixing the Q obtained in S2Wall(s)、QMaterial(s)And QSecondary airSubstituting the value of (A) into the rotary kiln heat balance model in the step S1 to obtain QBurning deviceAnd QFlue gasThe relation between the two, and then the flue gas temperature t of the calcining section is obtainedrTemperature t of kiln tail flue gasyThe relational expression of (1);
s4, measuring the temperature t of the kiln tail flue gasyAnd passing through the flue gas temperature t of the calcination sectionrTemperature t of kiln tail flue gasyObtaining the flue gas temperature t of the calcining section by the relational expressionr
The heat balance model specifically comprises: qSecondary air+QBurning device=QFlue gas+QWall(s)+QMaterial(s);QSecondary airThe energy of the high-temperature combustion-supporting air which directly enters the rotary kiln from the ring cooling section in unit time is expressed;
Qburning deviceRepresents all heat released after the fuel sprayed by the nozzle is completely combusted in unit time;
Qflue gasThe energy of the kiln tail outlet flue gas increased relative to the normal temperature flue gas temperature in unit time is represented;
Qwall(s)Indicating the heat dissipated by the wall of the kiln in a unit timeAn amount;
Qmaterial(s)Means the increased energy per unit time of pellets exiting the rotary kiln after the completion of the calcination section relative to the energy entering the kiln.
2. The method for online predicting of flue gas temperature at the calcining section in the rotary kiln as claimed in claim 1, wherein the specific step of S2 comprises:
s21: uniformly dividing the kiln wall of the combustion section of the rotary kiln into n sections;
s22: measuring the temperature t of each sectioniAnd surface area SiAnd the ambient temperature t corresponding to the kiln wall of the combustion section of the rotary kiln
S23: measuring the temperature t of each section according to the S22iAnd surface area SiAnd the ambient temperature t corresponding to the kiln wall of the combustion section of the rotary kilnObtaining the heat dissipation Q of the kiln wall facing the external environment in unit time when the rotary kiln stably worksWall(s)
3. The method for online prediction of flue gas temperature at the calcining section in the rotary kiln as claimed in claim 2, wherein in S23, QWall(s)The calculation formula of (a) is as follows:
Figure FDA0002369728340000021
in the formula, h represents the convective heat transfer coefficient between the wall surface of the kiln and the external environment.
4. The method for online predicting of flue gas temperature at the calcining section in the rotary kiln as claimed in claim 1, wherein the specific step of S3 comprises:
s31: measuring the mass m of finished pellets flowing out of a kiln head in unit time and the temperature t of finished pellets discharged from the kiln when the rotary kiln stably worksoutAnd the initial temperature t of green pellets entering the kiln per unit timein
S32: finished pellets flowing out of the kiln head per unit time as measured in said S31Mass m, temperature t of finished ball taken out of kilnoutAnd the initial temperature t of green pellets entering the kiln per unit timeinObtaining the energy Q increased by the pellets before and after entering the kilnMaterial(s)
5. The method for online prediction of flue gas temperature at the calcining section in the rotary kiln as claimed in claim 4, wherein in S32, QMaterial(s)The calculation formula of (a) is as follows:
Qmaterial(s)=m×cm×(tout-tin),
In the formula, cmRepresents the specific heat capacity of the pellets.
6. The method for on-line prediction of flue gas temperature of a calcining section in a rotary kiln as claimed in any one of claims 1-5, wherein the specific step of S4 comprises:
s41: measuring primary air temperature t sprayed from kiln head in unit time when rotary kiln stably works1Volume of secondary air V2And secondary air temperature t2
S42: according to the measured primary air temperature t sprayed from the kiln head in the unit time measured in the S421Volume of secondary air V2And secondary air temperature t2Obtaining the heat Q brought by the secondary air in unit timeSecondary air
7. The method for online prediction of flue gas temperature at the calcining section in the rotary kiln as claimed in claim 6, wherein in S42, QSecondary airThe calculation formula of (a) is as follows:
Qsecondary air=V2×ρ2×c2×(t2-t1),
In the formula, ρ2Representing secondary air gas density, c2Represents the specific heat capacity of the secondary air.
8. The method for on-line prediction of flue gas temperature at the calcining section in the rotary kiln as claimed in claim 7, wherein: qBurning deviceAnd QFlue gasThe relation between is (V)y-V2)×ρr×cr×(tr-t1)=Vy×ρy×cy×(ty-t1)-Q,
In the formula, VyThe volume of the outlet flue gas of the rotary kiln in unit time; rhoyIs the density of the outlet flue gas; c. CyIs the specific heat capacity of the outlet flue gas; t is tyIs the outlet flue gas temperature; rhorThe density of the high-temperature flue gas at the calcining section; c. CrThe specific heat capacity of the high-temperature flue gas at the calcining section; t is trIs the high temperature flue gas temperature of the calcining section.
9. The method for on-line prediction of the flue gas temperature of the calcining section in the rotary kiln according to any one of claims 1 to 5, characterized in that: the step between the step of S3 and the step of S4,
subjecting the calcination section flue gas temperature t obtained in the step S3 to simulationrTemperature t of kiln tail flue gasyThe relational expression (2) is verified, and if the verification is passed, the step S4 is executed, and if the verification is not passed, the step S2 is returned to be executed.
10. A system for on-line prediction of flue gas temperature of a calcining section in a rotary kiln is characterized by comprising the following modules,
a heat balance model building module for building Q-value in the rotary kiln according to the heat balance principle when the rotary kiln stably worksSecondary air、QBurning device、QFlue gas、QWall(s)And QMaterial(s)The heat balance model of (1);
a model parameter calculation module for calculating Q in the heat balance model according to the parameters of the rotary kiln in stable operationWall(s)、QMaterial(s)And QSecondary airThe value of (d);
a relational expression obtaining module for obtaining QWall(s)、QMaterial(s)And QSecondary airSubstituting the numerical value into a rotary kiln heat balance model to obtain QBurning deviceAnd QFlue gasThe relation between them, and then obtain the calcinationFlue gas temperature t of burning sectionrTemperature t of kiln tail flue gasyThe relational expression of (1);
a calculation module for measuring the kiln tail flue gas temperature tyAnd passing through the flue gas temperature t of the calcination sectionrTemperature t of kiln tail flue gasyObtaining the flue gas temperature t of the calcining section by the relational expressionr
The heat balance model specifically comprises: qSecondary air+QBurning device=QFlue gas+QWall(s)+QMaterial(s);QSecondary airThe energy of the high-temperature combustion-supporting air which directly enters the rotary kiln from the ring cooling section in unit time is expressed;
Qburning deviceRepresents all heat released after the fuel sprayed by the nozzle is completely combusted in unit time;
Qflue gasThe energy of the kiln tail outlet flue gas increased relative to the normal temperature flue gas temperature in unit time is represented;
Qwall(s)The heat quantity dissipated by the wall surface of the kiln in unit time is represented;
Qmaterial(s)Means the increased energy per unit time of pellets exiting the rotary kiln after the completion of the calcination section relative to the energy entering the kiln.
CN202010046835.6A 2020-01-16 2020-01-16 Method and system for online prediction of flue gas temperature of calcining section in rotary kiln Expired - Fee Related CN111174569B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010046835.6A CN111174569B (en) 2020-01-16 2020-01-16 Method and system for online prediction of flue gas temperature of calcining section in rotary kiln

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010046835.6A CN111174569B (en) 2020-01-16 2020-01-16 Method and system for online prediction of flue gas temperature of calcining section in rotary kiln

Publications (2)

Publication Number Publication Date
CN111174569A true CN111174569A (en) 2020-05-19
CN111174569B CN111174569B (en) 2021-10-22

Family

ID=70648007

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010046835.6A Expired - Fee Related CN111174569B (en) 2020-01-16 2020-01-16 Method and system for online prediction of flue gas temperature of calcining section in rotary kiln

Country Status (1)

Country Link
CN (1) CN111174569B (en)

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112347624A (en) * 2020-10-27 2021-02-09 华北水利水电大学 Glass energy consumption measuring method and device
CN113177332A (en) * 2021-03-09 2021-07-27 广东工业大学 Rotary kiln sintering temperature prediction method based on combination of mechanism and data
CN114018066A (en) * 2021-11-08 2022-02-08 攀钢集团攀枝花钢铁研究院有限公司 Method for predicting flue gas temperature of titanium slag electric furnace
CN114526827A (en) * 2022-02-28 2022-05-24 中冶长天国际工程有限责任公司 Double-hearth limekiln suspension cylinder temperature monitoring method and system and storage medium
CN116151022A (en) * 2023-03-07 2023-05-23 浙江大学 Real-time cement rotary kiln temperature estimation method based on heat balance calculation
CN116611356A (en) * 2023-04-23 2023-08-18 鞍钢股份有限公司 Method for determining convective heat transfer coefficient of titanium white calcination rotary kiln
CN117034788A (en) * 2023-04-23 2023-11-10 鞍钢股份有限公司 Titanium white calcination rotary kiln solid material temperature field calculation method based on continuous equation
CN117308578A (en) * 2023-10-26 2023-12-29 浙江金泰莱环保科技有限公司 Rotary kiln combustion adjusting method based on incineration experiment

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103307872A (en) * 2013-06-26 2013-09-18 中冶长天国际工程有限责任公司 Kiln flue gas temperature detection method and device
CN103323145A (en) * 2013-06-26 2013-09-25 中冶长天国际工程有限责任公司 Method and device for detecting flue-gas temperature in rotary kiln
CN104480300A (en) * 2014-11-20 2015-04-01 中南大学 Pellet production method based on prediction of compressive strength of pellets in rotary kiln
CN104915567A (en) * 2015-06-18 2015-09-16 南京南瑞继保电气有限公司 Cement clinker sintering thermal consumption, thermal efficiency and fire coal thermal value real-time monitoring method based on thermal balance analysis
CN107037787A (en) * 2016-02-03 2017-08-11 中冶长天国际工程有限责任公司 A kind of grate-kiln pelletizing burnup control method and device
CN109539812A (en) * 2017-09-21 2019-03-29 南京梅山冶金发展有限公司 A kind of sleeve kiln intelligent gas controlling device and control method

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103307872A (en) * 2013-06-26 2013-09-18 中冶长天国际工程有限责任公司 Kiln flue gas temperature detection method and device
CN103323145A (en) * 2013-06-26 2013-09-25 中冶长天国际工程有限责任公司 Method and device for detecting flue-gas temperature in rotary kiln
CN104480300A (en) * 2014-11-20 2015-04-01 中南大学 Pellet production method based on prediction of compressive strength of pellets in rotary kiln
CN104915567A (en) * 2015-06-18 2015-09-16 南京南瑞继保电气有限公司 Cement clinker sintering thermal consumption, thermal efficiency and fire coal thermal value real-time monitoring method based on thermal balance analysis
CN107037787A (en) * 2016-02-03 2017-08-11 中冶长天国际工程有限责任公司 A kind of grate-kiln pelletizing burnup control method and device
CN109539812A (en) * 2017-09-21 2019-03-29 南京梅山冶金发展有限公司 A kind of sleeve kiln intelligent gas controlling device and control method

Cited By (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112347624A (en) * 2020-10-27 2021-02-09 华北水利水电大学 Glass energy consumption measuring method and device
CN113177332A (en) * 2021-03-09 2021-07-27 广东工业大学 Rotary kiln sintering temperature prediction method based on combination of mechanism and data
CN113177332B (en) * 2021-03-09 2023-02-24 广东工业大学 Rotary kiln sintering temperature prediction method based on combination of mechanism and data
CN114018066A (en) * 2021-11-08 2022-02-08 攀钢集团攀枝花钢铁研究院有限公司 Method for predicting flue gas temperature of titanium slag electric furnace
CN114018066B (en) * 2021-11-08 2023-09-05 攀钢集团攀枝花钢铁研究院有限公司 Method for predicting flue gas temperature of titanium slag electric furnace
CN114526827A (en) * 2022-02-28 2022-05-24 中冶长天国际工程有限责任公司 Double-hearth limekiln suspension cylinder temperature monitoring method and system and storage medium
CN116151022A (en) * 2023-03-07 2023-05-23 浙江大学 Real-time cement rotary kiln temperature estimation method based on heat balance calculation
CN116151022B (en) * 2023-03-07 2024-04-23 浙江大学 Real-time cement rotary kiln temperature estimation method based on heat balance calculation
CN116611356A (en) * 2023-04-23 2023-08-18 鞍钢股份有限公司 Method for determining convective heat transfer coefficient of titanium white calcination rotary kiln
CN117034788A (en) * 2023-04-23 2023-11-10 鞍钢股份有限公司 Titanium white calcination rotary kiln solid material temperature field calculation method based on continuous equation
CN117308578A (en) * 2023-10-26 2023-12-29 浙江金泰莱环保科技有限公司 Rotary kiln combustion adjusting method based on incineration experiment
CN117308578B (en) * 2023-10-26 2024-09-03 浙江金泰莱环保科技有限公司 Rotary kiln combustion adjusting method based on incineration experiment

Also Published As

Publication number Publication date
CN111174569B (en) 2021-10-22

Similar Documents

Publication Publication Date Title
CN111174569B (en) Method and system for online prediction of flue gas temperature of calcining section in rotary kiln
CN111665809B (en) Segmentation mechanism modeling method and system for rotary cement kiln
WO2018076403A1 (en) One-dimensional simulation method for predicting quality of clinker in rotary cement kiln
CN108388759B (en) Horseshoe flame glass kiln energy consumption modeling and local energy consumption benchmarking method
Xu et al. Modeling of multiphase reaction and slag flow in single-burner coal water slurry gasifier
CN104316559A (en) Testing method capable of accurately reflecting dynamic thermal equilibrium of steel rolling heating furnace
CN113283079A (en) Method for calculating and monitoring boundary of blast furnace tuyere raceway in real time
CN114239430A (en) Method and system for predicting NOx at furnace outlet based on numerical simulation
CN102690072B (en) On-line monitoring and controlling method for incomplete combustion carbides in cement production process
CN114283895B (en) Method for monitoring temperature field of high-temperature reaction zone of zinc oxide rotary volatilizing kiln based on thermodynamic mechanism
Carrasco-Maldonado et al. Oxy-fuel burner investigations for CO2 capture in cement plants
CN103323145B (en) Method and device for detecting flue-gas temperature in rotary kiln
CN107066659A (en) A kind of method that limit of utilization learning machine predicts cement decomposing furnace temperature
Gerhardter et al. Assessment of a novel numerical model for combustion and in-flight heating of particles in an industrial furnace
CN103559382B (en) Lower calorific value online estimating method for pulverized coal as fired in cement decomposing furnace
Liu et al. Investigation of heat transfer phenomena in blast furnace tuyere/blowpipe region
Alyaser Fluid flow and combustion in rotary kiln
CN107025316A (en) A kind of method for monitoring water wall gasifier temperature
López et al. A numerical analysis of the effect of atmospheric pressure on the performance of a heating system with a self-recuperative burner
Rubekina et al. Calculation of the heat transfer in coke dry quenching by numerical modeling
Li et al. Measurements and modeling of fluid flow and thermal processes in an industrial precalciner
Fu et al. Research on soft measurement method of temperature field in cement rotary kiln
Knoll et al. Effects on numerical calculations of in-flight particle trajectories and temperatures considering multiple particle size and shape
CN113283078B (en) Method for calculating and monitoring depth of blast furnace tuyere convolution zone in real time
Schmidt Numerical Model for Heat Transfer Processes in Dry Ash Conveyors

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
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20211022