CN109472063B - Modeling method of energy efficiency model of hot galvanizing unit - Google Patents

Modeling method of energy efficiency model of hot galvanizing unit Download PDF

Info

Publication number
CN109472063B
CN109472063B CN201811223722.8A CN201811223722A CN109472063B CN 109472063 B CN109472063 B CN 109472063B CN 201811223722 A CN201811223722 A CN 201811223722A CN 109472063 B CN109472063 B CN 109472063B
Authority
CN
China
Prior art keywords
heat
zinc
temperature
strip steel
unit
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
Application number
CN201811223722.8A
Other languages
Chinese (zh)
Other versions
CN109472063A (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.)
China University of Geosciences
Original Assignee
China University of Geosciences
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 China University of Geosciences filed Critical China University of Geosciences
Priority to CN201811223722.8A priority Critical patent/CN109472063B/en
Publication of CN109472063A publication Critical patent/CN109472063A/en
Application granted granted Critical
Publication of CN109472063B publication Critical patent/CN109472063B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • CCHEMISTRY; METALLURGY
    • C23COATING METALLIC MATERIAL; COATING MATERIAL WITH METALLIC MATERIAL; CHEMICAL SURFACE TREATMENT; DIFFUSION TREATMENT OF METALLIC MATERIAL; COATING BY VACUUM EVAPORATION, BY SPUTTERING, BY ION IMPLANTATION OR BY CHEMICAL VAPOUR DEPOSITION, IN GENERAL; INHIBITING CORROSION OF METALLIC MATERIAL OR INCRUSTATION IN GENERAL
    • C23CCOATING METALLIC MATERIAL; COATING MATERIAL WITH METALLIC MATERIAL; SURFACE TREATMENT OF METALLIC MATERIAL BY DIFFUSION INTO THE SURFACE, BY CHEMICAL CONVERSION OR SUBSTITUTION; COATING BY VACUUM EVAPORATION, BY SPUTTERING, BY ION IMPLANTATION OR BY CHEMICAL VAPOUR DEPOSITION, IN GENERAL
    • C23C2/00Hot-dipping or immersion processes for applying the coating material in the molten state without affecting the shape; Apparatus therefor
    • C23C2/04Hot-dipping or immersion processes for applying the coating material in the molten state without affecting the shape; Apparatus therefor characterised by the coating material
    • C23C2/06Zinc or cadmium or alloys based thereon
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/08Thermal analysis or thermal optimisation

Abstract

The invention discloses a modeling method of an energy efficiency model of a hot galvanizing unit, which utilizes an energy flow theory, and based on energy conservation, analyzes the energy consumption composition of a heating furnace and a zinc pot in a production line of the hot galvanizing unit, and focuses on the energy flow composition and characteristics of the hot galvanizing unit; and then, respectively establishing a heat balance equation according to the internal relation between the output power and each technological parameter and a mapping mechanism, and summarizing the two parts of energy consumption according to a carbon emission principle to obtain a total energy efficiency model. The energy efficiency model established by the invention provides a basis for optimizing process parameters for later energy efficiency discussion and completing multi-objective sequencing optimization research. The establishment of the model is a theoretical basis for the deep follow-up study, and has definite and reasonable meaning and is indispensable.

Description

Modeling method of energy efficiency model of hot galvanizing unit
Technical Field
The invention belongs to the technical field of energy optimization of a galvanizing unit, and particularly relates to a modeling method of an energy efficiency model of a hot galvanizing unit.
Background
The steel industry is a typical flow industry, and the production process has the characteristics of complex process, strict production conditions, more production equipment, high automation degree and the like. Based on the analysis of cold rolling production line of a certain steel plant, the cold rolling plant is provided with three units of pickling, continuous annealing and galvanization in total. The galvanizing unit in the cold rolling process has huge energy saving potential, and the energy efficiency modeling of the whole production line is finished aiming at the galvanizing production line at present. Ventilation, speed, temperature are the core process parameters.
The area with the greatest energy consumption of the galvanizing unit is as follows: a heating furnace area and a zinc pot area. In order to perform effective energy-saving arrangement, the energy source and consumption of a working area must be analyzed, the mapping relation is found out, and an energy efficiency balance equation of two parts of the continuous annealing furnace and the zinc pot is established.
Disclosure of Invention
In view of the above, the embodiment of the invention provides a modeling method of an energy efficiency model capable of analyzing the relation between the energy source and the consumption of a galvanizing unit.
In order to solve the technical problems, the technical scheme adopted by the embodiment of the invention is that a modeling method of an energy efficiency model of a hot galvanizing unit comprises the following steps:
(1) According to the heat input of the heating furnace, the heat taken away by the strip steel, the heat dissipated by convection of the furnace wall, the heat dissipated by radiation of the furnace wall, the heat taken away by waste gas and protective gas respectively, calculating to obtain a heat balance equation of the heating furnace in the continuous annealing process;
(2) According to the heat input by the steel belt in the zinc pot, the heat supplied by induction heating, the heat dissipated by convection of the surface of the zinc liquid, the heat radiated by the surface of the zinc liquid, the heat dissipated by convection around the zinc pot, the heat radiated by the periphery of the zinc pot, the heat dissipated by convection at the bottom of the zinc pot, the heat radiated by the bottom of the zinc pot and the heat expended by melting zinc ingots, calculating to obtain a heat balance equation of the zinc pot;
(3) And correcting and adding carbon emission factors through the heat balance equation of the heating furnace and the heat balance equation of the zinc pot to obtain an energy efficiency model of the hot galvanizing unit.
The technical scheme provided by the embodiment of the invention has the beneficial effects that: according to the modeling method of the energy efficiency model of the hot galvanizing unit, the embodiment of the invention researches the energy consumption composition of the heating furnace and the zinc pot in the production line of the hot galvanizing unit based on energy conservation according to the energy flow theory, respectively establishes a heat balance equation according to the relation between the output power and each technological parameter, and obtains a total energy efficiency model through correction of carbon emission factors; the method can be used for analyzing a cold rolling production line, and can be used for carrying out ventilation, speed and temperature experimental setting on core technological parameters according to the model, so as to achieve the effects of saving energy, reducing emission and improving yield.
Drawings
FIG. 1 is a schematic diagram of the energy efficiency model composition of a continuous hot dip galvanizing unit module according to an embodiment of the present invention;
FIG. 2 is a schematic diagram showing the comparison of calculated and actual parameters in annealing furnace strip steel according to the method of the embodiment of the invention;
FIG. 3 is a schematic diagram showing the comparison of calculated parameters and actual parameters in a zinc-boiler strip steel according to the method of the embodiment of the invention.
Detailed Description
For the purpose of making the objects, technical solutions and advantages of the present invention more apparent, embodiments of the present invention will be further described with reference to the accompanying drawings.
Example 1
Referring to fig. 1, an embodiment of the present invention provides a modeling method for an energy efficiency model of a hot galvanizing unit, including the following steps:
(1) According to the heat input Q of the heating furnace rq Heat quantity Q taken away by strip steel dg2 Heat loss by convection Q of furnace wall ld Radiation heat dissipation Q of furnace wall lf Heat Q taken away by waste gas and protective gas respectively fq 、Q bh Calculating to obtain a heat balance equation of the heating furnace in the continuous annealing process;
specifically, the heat input Q of the heating furnace rq For gas supply quantity and heat quantity Q taken away by strip steel dg2 Heat loss by convection Q of furnace wall ld Radiation heat dissipation Q of furnace wall lf Heat Q taken away by waste gas and protective gas respectively fq 、Q bh Obtained by the formulas (1) to (6), respectively,
Q rq =V rq h rq (1)
Figure SMS_1
Q ld =α l (T l -T h )S l (3)
Figure SMS_2
Q fq =V fq C fq (T f -T h ) (5)
Q bh =V bh C bh (T q2 -T q1 ) (6)
wherein Q is rq Heat quantity Q obtained for heating furnace dg2 Heat quantity Q taken away by strip steel ld Heat, Q, dissipated for convection to the surface of the furnace wall lf Heat exchange quantity Q for furnace wall radiation fq Heat taken away by waste gas, Q bh The heat quantity taken away by the protective gas is kJ/h; v (V) rq The unit is m, the gas quantity introduced in the standard state and the V is the total quantity of the waste gas in the standard state 3 /h;h rq Is the heat value of fuel gas, and has the unit of kJ/m 3 ;ρ 2 The density of the strip steel is 7.85 multiplied by 10 3 kg/m 3 The method comprises the steps of carrying out a first treatment on the surface of the b is the width of the strip steel, h g The thickness of the strip steel is m; v is the speed of the strip steel, and the unit is m/min; c (C) s The unit is kJ/(kg. Deg.C) which is the mass heat capacity of steel; s is S l Is the heat dissipation area of the furnace wall, and the unit is m 2 ;T ck Is the outlet temperature T of the strip steel i Is the instantaneous temperature T of the strip steel rk Is the inlet temperature T of the strip steel l Is the temperature of furnace wall, T h Is the ambient temperature, T f Is the temperature of the exhaust gas, T q2 Temperature, T, of the shielding gas heated by the radiant tube q1 The inlet temperature of the protective gas is in the unit of DEG C; alpha l The unit of the convection heat transfer coefficient of the furnace wall surface is kJ/(m) 2 ·h·℃);ε lb Is the blackness of the system;
Figure SMS_3
is an angle coefficient; c (C) 1 Is a blackbody radiation coefficient of 4.96×10 3 kJ/(m 2 ·h·K 4 );T l Is the surface temperature of the furnace wall, T h The unit is K for the ambient temperature; c is the heat capacity of the exhaust gas, kJ/(m) 3 ·℃);
The heat capacity of the strip steel is a variable quantity which continuously changes along with the temperature, so that the strip steel cannot be regarded as a constant, and the calculation formula is (7):
C s =1.34×10 -11 t 5 -3.7×10 -8 t 4 +4.007×10 -5 t 3 -0.02101t 2 +5.672t-179.6 (7);
further, the total heat income Q of the heating furnace in the continuous annealing process ly Mainly the gas supply quantity, namely the heat Q obtained by a heating furnace rq The total heat income Q of the heating furnace in the continuous annealing process is obtained by the formula (8) lr Total heat input Q during continuous annealing lr And total heat expenditure Q lc Equilibrium, obtaining a thermal equilibrium equation (12) according to equations (8), (9), (10), (11);
Q lr =Q rq (8)
Q lc =Q dg2 +Q fq +Q ld +Q lf +Q bh (9)
Q lr =Q lc (10)
Q rq =Q dg2 +Q fq +Q ld +Q lf +Q bh (11)
Figure SMS_4
(2) According to the heat Q brought in by the steel belt in the zinc pot dg1 Induction heating to supply heat Q gy Convective heat dissipation Q on surface of zinc liquid bd Surface radiation heat exchange quantity Q of zinc liquid bf Convection heat dissipation Q around zinc pot zd Heat exchanging quantity Q of radiation around zinc pot zf Convection heat dissipation Q at bottom of zinc pot dd Heat exchange quantity Q of zinc pot bottom radiation df Heat of melting zinc ingot Q r Calculating to obtain a heat balance equation of the zinc pot;
specifically, the strip steel in the zinc pot enters heat Q dg1 Induction heating to supply heat Q gy Convective heat dissipation Q on surface of zinc liquid bd Surface radiation heat exchange quantity Q of zinc liquid bf Convection heat dissipation Q around zinc pot zd Heat exchanging quantity Q of radiation around zinc pot zf Convection heat dissipation Q at bottom of zinc pot dd Heat exchange quantity Q of zinc pot bottom radiation df Heat of melting zinc ingot Q r Respectively lead toObtained by the formulas (13) to (21),
Q dg1 =60vbh g ρ 2 ·C s (T dg -T xy ) (13)
Q gy =3600N (14)
Q bd =α b (T b -T h )S b (15)
Figure SMS_5
Q zd =α z (T z -T h )S z (17)
Figure SMS_6
Q dd =α d (T d -T h )S d (19)
Figure SMS_7
Figure SMS_8
wherein Q is dg1 Heat quantity Q of strip steel brought into zinc pot gy Supplying heat, Q, to a zinc pot for induction heating bd Heat quantity Q dissipated by convection of the surface of the zinc liquid bf Heat exchange capacity Q of the surface radiation of the zinc liquid r Heat, Q, expended for melting zinc ingots zd For convection heat dissipation around zinc pot, Q zf For the radiation heat exchange quantity, Q of the periphery of the zinc pot dd Heat dissipation and Q are given off by convection at the bottom of the zinc pot df The heat exchange quantity is radiated to the bottom of the zinc pot, and the units are kJ/h; ρ 2 The density of the strip steel is 7.85 multiplied by 10 3 kg/m 3 ;ρ 1 Is zinc with a density of 7.14X10 3 kg/m 3 The method comprises the steps of carrying out a first treatment on the surface of the b is the width of the strip steel, h g Is the thickness of the strip steel, h x The thickness of the zinc layer is m; v is a beltThe speed of the steel, m/min; c (C) s Is the mass heat capacity of steel, C z The unit is kJ/(kg. DEG C) for the mass heat capacity of zinc; t (T) dg Is the temperature T of the strip steel entering the zinc pot xy Is the working temperature T of zinc liquid b Is the temperature of zinc liquid, T h Is the ambient temperature, T g Is the working temperature T of zinc liquid y The original temperature of the zinc ingot is set at the temperature of DEG C; n is active power, kW; alpha b Is the surface convection heat transfer coefficient of the zinc liquid, kJ/(m) 2 ·h·℃);ε xy The blackness of the zinc liquid surface with zinc oxide ash to the ambient air is the system blackness;
Figure SMS_9
the angular coefficient of the zinc liquid level to the environment heat absorption surface; c (C) 0 For blackbody emissivity, kJ/(m) 2 ·h·K 4 ) The method comprises the steps of carrying out a first treatment on the surface of the L is latent heat of zinc ingot melting, kJ/kg;
further, the total heat expenditure and the heat income of the induction heating zinc pot are equal during production, a heat balance equation (23) of the zinc pot is obtained through formulas (13) to (22),
Q dg1 +Q gy =Q bd +Q bf +Q zd +Q zf +Q dd +Q df +Q r (22)
Figure SMS_10
(3) And correcting and adding carbon emission factors through the heat balance equation of the heating furnace and the heat balance equation of the zinc pot to obtain an energy efficiency model of the hot galvanizing unit.
Specifically, according to the principle of energy conservation, the energy of the two can be simply fused from the angle of energy, the electric energy consumed by the zinc pot and the heat energy consumed by the heating furnace are taken as two independent and related energies, and the two energies are unified by adopting carbon emission factor correction, so that the energy of different sources is measured; according to the heat balance equations (12) and (22) obtained in the steps (1) and (2), the carbon emission factors, the electricity consumption consumed by the zinc pot and the heat consumption consumed by the heating furnace, obtaining a carbon emission energy efficiency model (27) of unit mass of the hot galvanizing unit through formulas (24) to (26),
Figure SMS_11
Figure SMS_12
Figure 10000263555558
Figure SMS_14
Figure SMS_15
wherein, the liquid crystal display device comprises a liquid crystal display device,
Figure SMS_16
c as total carbon emission iEF,elc ,C iEF,gas Carbon emission factors expressed as electric energy and thermal energy, respectively; Σe Ci ,∑G Ci Respectively expressed as the electricity consumption of the zinc pot and the heat consumption of the heating furnace.
Example two
According to the method provided by the embodiment of the invention, the energy efficiency calculation of the hot galvanizing unit is compared with the energy efficiency of M1A1 steel in the actual hot galvanizing unit.
Before heating the strip steel, firstly preheating the strip steel: the exhaust gas is passed through a heat exchanger to transfer heat to the shielding gas, the strip is preheated by a preheating zone having a total path of 44m without any corrosion and maintaining a speed, and the outlet temperature is about 180-200 ℃.
(1) After the strip steel is heated, a part of heat is taken away, the thickness of the steel grade is 0.69mm, the width is 1479mm, the running speed is 110m/min, and the heat Q taken away by the strip steel can be obtained by the formulas (2) and (7) dg2
(2) Part of heat is taken away by the protective gas and the waste gas, so that heat loss is caused, and Q is obtained by the formula (5) fq From formula (6), Q bh Is a measure of (2);
(3) Heat Q of convection and radiation heat exchange of heating furnace surface ld 、Q lf The jetness of the engineering materials used in general can be obtained from the formula (3) and the formula (4), respectively, as shown in Table 1:
TABLE 1 surface jetness of commonly used engineering materials
Figure SMS_17
(4) The total heat expenditure is 46500MJ/h, the heat income is the heat energy brought by gas heating, and the heat value of the burnt mixed gas is 7530+/-418KJ/Nm 3 From this, it can be calculated according to the formula (12) that the ventilation which should be required for heating the strip is 6175.299m 3 /h。
According to the actual monitoring data collected on site, the steel grade passes through seven heating zones in the heating furnace section, and the specific ventilation is shown in table 2:
table 2M2A1 strip heating segment ventilation data
Figure SMS_18
At the moment, the ventilation rate of the actual demand of the strip steel is 6373.759m 3 And/h, obtaining an error by comparing the actual quantity with the calculated quantity.
The ventilation of the same type of strip steel with different specifications and the ventilation of the different types of strip steel are calculated respectively and compared with the actual ventilation, and specific data are obtained as shown in table 3:
TABLE 3 comparison table of calculated parameters and actual parameters of strip steel of annealing furnace
Figure SMS_19
Referring to fig. 2, it is shown that the heat balance equation and the energy efficiency model established by the embodiment of the present invention can effectively predict the ventilation in practical application.
The ceramic induction zinc pot is used in the hot galvanizing unit, the energy sources are mainly the heat of the cooled steel belt and the heating of electric energy, and the output is mainly the heat expenditure of the induction heating zinc pot including the heat dissipation of the surface of the zinc liquid (convection heat exchange Q bd And radiation heat exchange Q bf ) Heat dissipation around the zinc pot (convection heat exchange Q) zd And radiation heat exchange Q zf ) Heat dissipation at the bottom of the zinc pot (convection heat exchange Q) dd And radiation heat exchange Q df ) And heat Q expended in melting zinc ingot r
(1) The strip steel is provided with a temperature before entering the zinc pot, if the temperature is higher than the temperature of the zinc liquid, energy income is brought to the whole, the same strip steel is selected as an annealing furnace, and the heat Q of entering the zinc pot from the strip steel is obtained by taking in (13) dg1
(2) The surface area of the zinc liquid is 12m 2 Convection and radiation heat dissipation quantity Q of zinc liquid surface bd 、Q bf Calculated from the formula (15) and the formula (16); the convection heat exchange quantity Q of the zinc pot can be obtained by the formula (17) and the formula (18) zd And radiation heat exchange Q zf The method comprises the steps of carrying out a first treatment on the surface of the The convection heat transfer Q of the bottom of the zinc pot is obtained by the formula (19) and the formula (20) dd And radiation heat exchange Q df
(3) The zinc ingot absorbs energy when melting, the strip steel galvanization is double-sided galvanization, and each surface is 80g/m 2 The required energy Q at this time is calculated by the formula (21) r
(4) The sum of the energy and the electric energy carried in by the strip is equal to the total heat expenditure, and the electric power required at this time is 201kW, which is obtainable by the formula (23).
Respectively calculating the electric power of the strip steel with the same type and different specifications and the electric power of the strip steel with different types according to the verification data in the annealing furnace, and comparing the electric power with the actual electric power to obtain specific data as shown in table 4
Table 4 comparison table of calculated and actual parameters of zinc pot strip steel
Figure SMS_20
Referring to fig. 3, the thermal equilibrium equation, the energy efficiency model and the data error monitored in the practical application established by the embodiment of the invention are shown to be small, the model precision is high, and the accuracy is good.
In this document, terms such as front, rear, upper, lower, etc. are defined with respect to the positions of the components in the drawings and with respect to each other, for clarity and convenience in expressing the technical solution. It should be understood that the use of such orientation terms should not limit the scope of the protection sought herein.
The embodiments described above and features of the embodiments herein may be combined with each other without conflict.
The foregoing description of the preferred embodiments of the invention is not intended to limit the invention to the precise form disclosed, and any such modifications, equivalents, and alternatives falling within the spirit and scope of the invention are intended to be included within the scope of the invention.

Claims (5)

1. The modeling method of the energy efficiency model of the hot galvanizing unit is characterized by comprising the following steps of:
(1) According to the heat input Q of the heating furnace rq Heat quantity Q taken away by strip steel dg2 Heat loss by convection Q of furnace wall ld Radiation heat dissipation Q of furnace wall lf Heat Q taken away by waste gas and protective gas respectively fq 、Q bh Calculating to obtain a heat balance equation of the heating furnace in the continuous annealing process;
(2) According to the heat Q brought in by the steel belt in the zinc pot dg1 Induction heating to supply heat Q gy Convective heat dissipation Q on surface of zinc liquid bd Surface radiation heat exchange quantity Q of zinc liquid bf Convection heat dissipation Q around zinc pot zd Heat exchanging quantity Q of radiation around zinc pot zf Convection heat dissipation Q at bottom of zinc pot dd Heat exchange quantity Q of zinc pot bottom radiation df Heat of melting zinc ingot Q r Calculating to obtain a heat balance equation of the zinc pot;
(3) The energy efficiency model of the hot galvanizing unit is obtained by adding corrected carbon emission factors through the heat balance equation of the heating furnace and the heat balance equation of the zinc pot;
wherein, the liquid crystal display device comprises a liquid crystal display device,
Figure FDA0004234463470000011
Figure FDA0004234463470000012
Figure FDA0004234463470000013
Figure FDA0004234463470000014
Figure FDA0004234463470000021
wherein, the parameters expressed when i=1, 2 and 3 are parameters related to the calculation of the convection heat dissipation or radiation heat exchange quantity of the surface of the zinc liquid, the periphery of the zinc pot and the bottom of the zinc pot, and V rq The unit is m, the gas quantity introduced in the standard state and the V is the total quantity of the waste gas in the standard state 3 /h;h rq Is the heat value of fuel gas, and has the unit of kJ/m 3 ;ρ 2 The density of the strip steel is 7.85 multiplied by 10 3 kg/m 3 The method comprises the steps of carrying out a first treatment on the surface of the b is the width of the strip steel, h g The thickness of the strip steel is m; v is the speed of the strip steel, and the unit is m/min; c (C) s The unit is kJ/(kg. Deg.C) which is the mass heat capacity of steel; s is S l Is the heat dissipation area of the furnace wall, and the unit is m 2 ;T ck Is the outlet temperature T of the strip steel i Is the instantaneous temperature T of the strip steel rk Is the inlet temperature T of the strip steel l Is the temperature of furnace wall, T h Is the ambient temperature, T f Is the temperature of the exhaust gas, T q2 Temperature, T, of the shielding gas heated by the radiant tube q1 The inlet temperature of the protective gas is in the unit of DEG C; alpha l The unit of the convection heat transfer coefficient of the furnace wall surface is kJ/(m) 2 ·h·℃);ε lb Is of a system blackA degree;
Figure FDA0004234463470000022
is an angle coefficient; c (C) 1 Is a blackbody radiation coefficient of 4.96×10 3 kJ/(m 2 ·h·K 4 ) The method comprises the steps of carrying out a first treatment on the surface of the C is the heat capacity of the exhaust gas, kJ/(m) 3 ·℃);ρ 1 Is zinc with a density of 7.14X10 3 kg/m 3 ;h x The thickness of the zinc layer is m; c (C) z The unit is kJ/(kg. DEG C) for the mass heat capacity of zinc; t (T) dg Is the temperature T of the strip steel entering the zinc pot xy Is the working temperature T of zinc liquid b Is the temperature of zinc liquid, T g Is the working temperature T of zinc liquid y The original temperature of the zinc ingot is set at the temperature of DEG C; n is active power, kW; alpha b Is the surface convection heat transfer coefficient of the zinc liquid, kJ/(m) 2 ·h·℃);ε xy The blackness of the zinc liquid surface with zinc oxide ash to the ambient air is the system blackness; />
Figure FDA0004234463470000024
The angular coefficient of the zinc liquid level to the environment heat absorption surface; c (C) 0 For blackbody emissivity, kJ/(m) 2 ·h·K 4 ) The method comprises the steps of carrying out a first treatment on the surface of the L is latent heat of zinc ingot melting, kJ/kg; />
Figure FDA0004234463470000023
C as total carbon emission iEF,elc ,C iEF,gas Carbon emission factors expressed as electric energy and thermal energy, respectively; Σe Ci ,∑G Ci Respectively expressed as the electricity consumption of the zinc pot and the heat consumption of the heating furnace.
2. The modeling method of an energy efficiency model of a hot galvanizing unit according to claim 1, wherein the heat input quantity Q of the heating furnace rq For gas supply quantity and heat quantity Q taken away by strip steel dg2 Heat loss by convection Q of furnace wall ld Radiation heat dissipation Q of furnace wall lf Heat Q taken away by waste gas and protective gas respectively fq 、Q bh Obtained by the formulas (5) to (10), respectively,
Q rq =V rq h rq (5)
Figure FDA0004234463470000031
Q ld =α l (T l -T h )S l (7)
Figure FDA0004234463470000032
Q fq =V fq C fq (T f -T h ) (9)
Q bh =V bh C bh (T q2 -T q1 ) (10)
wherein Q is rq Heat quantity Q obtained for heating furnace dg2 Heat quantity Q taken away by strip steel ld Heat, Q, dissipated for convection to the surface of the furnace wall lf Heat exchange quantity Q for furnace wall radiation fq Heat taken away by waste gas, Q bh The heat quantity taken away by the protective gas is kJ/h; v (V) rq The unit is m, the gas quantity introduced in the standard state and the V is the total quantity of the waste gas in the standard state 3 /h;h rq Is the heat value of fuel gas, and has the unit of kJ/m 3 ;ρ 2 The density of the strip steel is 7.85 multiplied by 10 3 kg/m 3 The method comprises the steps of carrying out a first treatment on the surface of the b is the width of the strip steel, h g The thickness of the strip steel is m; v is the speed of the strip steel, and the unit is m/min; c (C) s The unit is kJ/(kg. Deg.C) which is the mass heat capacity of steel; s is S l Is the heat dissipation area of the furnace wall, and the unit is m 2 ;T ck Is the outlet temperature T of the strip steel i Is the instantaneous temperature T of the strip steel rk Is the inlet temperature T of the strip steel l Is the temperature of furnace wall, T f Is the temperature of the exhaust gas, T q2 Temperature, T, of the shielding gas heated by the radiant tube q1 The inlet temperature of the protective gas is in the unit of DEG C; alpha l Is a furnace wall surface pairThe flow heat exchange coefficient is in kJ/(m) 2 ·h·℃);ε lb Is the blackness of the system;
Figure FDA0004234463470000033
is an angle coefficient; c (C) 1 Is a blackbody radiation coefficient of 4.96×10 3 kJ/(m 2 ·h·K 4 );T l Is the surface temperature of the furnace wall, T h The unit is K for the ambient temperature; c is the heat capacity of the exhaust gas, kJ/(m) 3 ·℃)。
3. The modeling method of energy efficiency model of hot galvanizing unit according to claim 2, wherein the total heat income Q of the heating furnace in the continuous annealing process is obtained by the formula (11) lr Total heat input Q during continuous annealing lr And total heat expenditure Q lc Equilibrium, obtaining a thermal equilibrium equation (15) according to equations (11), (12), (13), (14);
Q lr =Q rq (11)
Q lc =Q dg2 +Q fq +Q ld +Q lf +Q bh (12)
Q lr =Q lc (13)
Q rq =Q dg2 +Q fq +Q ld +Q lf +Q bh (14)
Figure FDA0004234463470000041
4. the modeling method of an energy efficiency model of a hot galvanizing unit according to claim 1, wherein the heat Q is introduced into the steel strip in the zinc pot dg1 Induction heating to supply heat Q gy Convective heat dissipation Q on surface of zinc liquid bd Surface radiation heat exchange quantity Q of zinc liquid bf Convection heat dissipation Q around zinc pot zd Heat exchanging quantity Q of radiation around zinc pot zf Convection heat dissipation Q at bottom of zinc pot dd Heat exchange quantity Q of zinc pot bottom radiation df Heat of melting zinc ingot Q r Obtained by the formulas (16) to (24), respectively,
Q dg1 =60vbh g ρ 2 ·C s (T dg -T xy ) (16)
Q gy =3600N (17)
Q bd =α b (T b -T h )S b (18)
Figure FDA0004234463470000042
Q zd =α z (T z -T h )S z (20)
Figure FDA0004234463470000043
Q dd =α d (T d -T h )S d (22)
Figure FDA0004234463470000044
Figure FDA0004234463470000045
wherein Q is dg1 Heat quantity Q of strip steel brought into zinc pot gy Supplying heat, Q, to a zinc pot for induction heating bd Heat quantity Q dissipated by convection of the surface of the zinc liquid bf Heat exchange capacity Q of the surface radiation of the zinc liquid r Heat, Q, expended for melting zinc ingots zd For convection heat dissipation around zinc pot, Q zf For the radiation heat exchange quantity, Q of the periphery of the zinc pot dd Heat dissipation and Q are given off by convection at the bottom of the zinc pot df The heat exchange quantity is radiated to the bottom of the zinc pot, and the units are kJ/h; ρ 2 The density of the strip steel is 7.85 multiplied by 10 3 kg/m 3 ;ρ 1 Is zinc with a density of 7.14X10 3 kg/m 3 The method comprises the steps of carrying out a first treatment on the surface of the b is the width of the strip steel, h g Is the thickness of the strip steel, h x The thickness of the zinc layer is m; v is the speed of the strip steel, m/min; c (C) s Is the mass heat capacity of steel, C z The unit is kJ/(kg. DEG C) for the mass heat capacity of zinc; t (T) dg Is the temperature T of the strip steel entering the zinc pot xy Is the working temperature T of zinc liquid b Is the temperature of zinc liquid, T h Is the ambient temperature, T g Is the working temperature T of zinc liquid y The original temperature of the zinc ingot is set at the temperature of DEG C; n is active power, kW; alpha b Is the surface convection heat transfer coefficient of the zinc liquid, kJ/(m) 2 ·h·℃);ε xy The blackness of the zinc liquid surface with zinc oxide ash to the ambient air is the system blackness;
Figure FDA0004234463470000052
the angular coefficient of the zinc liquid level to the environment heat absorption surface; c (C) 0 For blackbody emissivity, kJ/(m) 2 ·h·K 4 ) The method comprises the steps of carrying out a first treatment on the surface of the L is the latent heat of fusion of the zinc ingot and kJ/kg.
5. The method for modeling an energy efficiency model of a hot dip galvanizing unit according to claim 4, wherein the thermal balance equation (26) of the zinc pot is obtained by formulas (16) to (25),
Q dg1 +Q gy =Q bd +Q bf +Q zd +Q zf +Q dd +Q df +Q r (25)
Figure FDA0004234463470000051
the parameters represented by i=1, 2 and 3 are parameters related to the calculation of the convection heat dissipation or radiation heat exchange quantity of the surface of the zinc liquid, the periphery of the zinc pot and the bottom of the zinc pot respectively.
CN201811223722.8A 2018-10-19 2018-10-19 Modeling method of energy efficiency model of hot galvanizing unit Active CN109472063B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811223722.8A CN109472063B (en) 2018-10-19 2018-10-19 Modeling method of energy efficiency model of hot galvanizing unit

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811223722.8A CN109472063B (en) 2018-10-19 2018-10-19 Modeling method of energy efficiency model of hot galvanizing unit

Publications (2)

Publication Number Publication Date
CN109472063A CN109472063A (en) 2019-03-15
CN109472063B true CN109472063B (en) 2023-06-30

Family

ID=65663963

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811223722.8A Active CN109472063B (en) 2018-10-19 2018-10-19 Modeling method of energy efficiency model of hot galvanizing unit

Country Status (1)

Country Link
CN (1) CN109472063B (en)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112464426B (en) * 2019-09-09 2022-11-15 宝山钢铁股份有限公司 Method for calculating temperature of strip steel entering zinc pot
CN112446134B (en) * 2020-10-23 2022-10-14 中南大学 Method for calculating heat loss of furnace body region in electric arc furnace steelmaking process
CN113553707B (en) * 2021-07-21 2022-09-23 中国人民解放军63798部队 Fairing temperature and humidity modeling method based on energy and mass conservation

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2005248208A (en) * 2004-03-01 2005-09-15 Nippon Steel Corp Method for controlling sheet temperature of galvannealed steel sheet during induction heating it
GB201519235D0 (en) * 2015-10-30 2015-12-16 Sharp Kabushiki Kaisha And Faradion Ltd Formation method for sodium ion cell or battery
WO2015192273A1 (en) * 2014-06-19 2015-12-23 国家电网公司 Self-adaptive modeling apparatus for three-phase alternating-current electric arc furnace, and simulation algorithm therefor
CN105219929A (en) * 2014-05-30 2016-01-06 宝山钢铁股份有限公司 The sectional type method of reducing of hot-rolled sheet scale on surface
CN108038319A (en) * 2017-12-20 2018-05-15 中国地质大学(武汉) A kind of method and system of galvanizing unit energy stream energy optimization modeling
CN108034804A (en) * 2017-12-08 2018-05-15 中国地质大学(武汉) A kind of method and system of continuous annealing unit stove area energy consumption modeling

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20050051208A1 (en) * 2003-06-17 2005-03-10 Mount Robert L. System for transferring heat in a thermoelectric generator system
US10127342B2 (en) * 2016-04-08 2018-11-13 Nano And Advanced Materials Institute Limited Method of designing and modifying lithium ion battery cathode materials

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2005248208A (en) * 2004-03-01 2005-09-15 Nippon Steel Corp Method for controlling sheet temperature of galvannealed steel sheet during induction heating it
CN105219929A (en) * 2014-05-30 2016-01-06 宝山钢铁股份有限公司 The sectional type method of reducing of hot-rolled sheet scale on surface
WO2015192273A1 (en) * 2014-06-19 2015-12-23 国家电网公司 Self-adaptive modeling apparatus for three-phase alternating-current electric arc furnace, and simulation algorithm therefor
GB201519235D0 (en) * 2015-10-30 2015-12-16 Sharp Kabushiki Kaisha And Faradion Ltd Formation method for sodium ion cell or battery
CN108034804A (en) * 2017-12-08 2018-05-15 中国地质大学(武汉) A kind of method and system of continuous annealing unit stove area energy consumption modeling
CN108038319A (en) * 2017-12-20 2018-05-15 中国地质大学(武汉) A kind of method and system of galvanizing unit energy stream energy optimization modeling

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
基于碳排放定量分析的低能耗切削加工研究进展;卢晓红等;《组合机床与自动化加工技术》;20151027(第10期);第1-4页 *
酸洗连轧生产线机组轧制功率建模及工艺参数优化;杨杰等;《中国机械工程》;20180929(第19期);第2371-2376页 *
锌锅热平衡测试应用的深入;董建军等;《宝钢技术》;20110815(第04期);第68-72页 *

Also Published As

Publication number Publication date
CN109472063A (en) 2019-03-15

Similar Documents

Publication Publication Date Title
CN109472063B (en) Modeling method of energy efficiency model of hot galvanizing unit
CN106636606B (en) A kind of method for controlling furnace temperature of heating furnace based on simulation model
CN106906351B (en) A kind of board briquette forecasting model and optimum furnace method
CN105423334B (en) Combustion of Hot Air Furnace Intelligent Process Control system and method
CN101736111A (en) Intelligent control method for automatic optimization combustion of hot-blast stove
CN106906350B (en) A kind of H profile steel base heating process temperature distribution calculation method
CN109248928A (en) A kind of hot-rolling heating furnace dynamic method for controlling furnace temperature
CN109182731A (en) A kind of high-carbon-chromium bearing steel continuous casting billet method for heating and controlling based on temperature-time control
CN103243190B (en) Method for predicting coal gas consumption of hot blast stove
CN107391894A (en) A kind of labyrinth calculation of radiation heat transferring method
CN104894362A (en) Method for setting temperature of heating furnace in cold and hot steel billet mixed loading
CN113343514A (en) Method for optimizing heating system of walking beam furnace
CN109022757A (en) The dynamic segmentation control method of tin-plated product annealing temperature
CN108038319B (en) Method and system for optimizing and modeling energy flow energy consumption of galvanizing unit
CN106766943A (en) A kind of vertical cooling furnace apparatus and cooling means of sintering deposit
Díaz-Ibarra et al. Design of a day tank glass furnace using a transient model and steady-state computation fluid dynamics
CN104049649A (en) Model control method of heating furnace temperature
CN103471393A (en) Blast furnace gas double preheating and ignition furnace temperature control method and device
CN207176005U (en) A kind of S types multilayer aluminium bar solution heat stove
CN108984943A (en) Heating furnace steel billet temperature trace model modification method
CN102994730B (en) Method for tracking temperature of steel billet in heating furnace
CN106282533A (en) A kind of temperature-controlled process to be rolled of heating furnace
CN106987704B (en) A kind of cold and hot process temperature control method for loading in mixture slab of pulse type heating furnace
CN102184271B (en) Method for designing vulcanization furnace based on heat flow coupling analysis technique
CN101638717B (en) Method for determining master-slave control proportional factor of hot rolling heating furnace

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