Intelligent kiln control method for tracking multi-energy combustion carbon emission
Technical Field
The invention relates to the technical field of kilns, in particular to an intelligent kiln control method for tracking multi-energy combustion carbon emission.
Background
Industrial kilns are equipment for thermal processing or calcining of products by consuming energy to generate heat, and require a large amount of energy to supply heat and burn out a large amount of carbon dioxide during operation. At present, kiln adjustment modes of continuous kilns (including roller kilns, tunnel kilns, push plate kilns and the like) are mostly full manual or semi-automatic, and when a production line changes production, due to the change of products and yields, energy consumption and smoke discharge quantity change, the frequency of a fan and valve elements are required to be manually adjusted, the operation is complex, and the product is unqualified due to human errors easily.
At present, an advanced kiln can adopt an electric and gas multi-energy mixed combustion mode to provide heat, energy source change is carried out according to the economy of different heating modes, and a fan and a valve of the kiln are automatically adjusted. However, in the process of changing energy types, the kiln needs to have stable internal temperature to ensure normal firing of products; therefore, how much input amount of the newly used energy or energy combination can just replace the original energy is important for kiln operation; however, due to the influence of various factors such as internal temperature, smoke amount, environmental temperature and the like in the operation process of the kiln, the energy source change cannot adopt a simple proportional relationship; at present, no method is available for estimating the energy input quantity during energy change; the operation of gradually increasing new energy and reducing original energy, increasing new energy and reducing original energy again after the temperature in the kiln is stable until the new energy is completely changed can only be adopted, the energy changing process is slow, and meanwhile, the accurate control and adjustment of the temperature in the kiln are not facilitated. In addition, the kiln lacks monitoring of carbon emission indexes, and the carbon emission of kiln combustion cannot be ensured to meet the environmental protection requirement.
Disclosure of Invention
In order to overcome the defects and shortcomings in the prior art, the invention aims to provide an intelligent kiln control method for tracking the carbon emission of multi-energy combustion; the method is based on the balance of kiln heat income and expenditure, builds a mathematical model of kiln heat balance, combines carbon emission indexes and economic indexes to obtain the energy source type proportion and corresponding input quantity, can meet the carbon emission requirement, reduce the kiln use cost, enable the combustion to reach the environment-friendly state of optimal economy, can quickly realize energy source type change, stabilize the kiln internal temperature and ensure the product firing quality.
In order to achieve the above purpose, the invention is realized by the following technical scheme: an intelligent kiln control method for tracking multi-energy burning carbon emission is used for kilns which comprehensively use fuel, green electricity and ash electricity as energy sources; the method comprises the following steps:
s1, constructing a kiln heat balance mathematical model with self-adaptive parameters; the adaptive parameters include fuel type, fuel consumption flow m r And electric heating real-time power W d The method comprises the steps of carrying out a first treatment on the surface of the The kiln heat balance mathematical model is as follows:
total kiln heat input Qsr = fuel carry-in sensible heat Q x (m r ) +Fuel Combustion exotherm Q r (m r ) +combustion air brought into sensible heat Q zk (m r ) +electric heating quantity Q d (W d ) +green body brought into sensible heat Q pt +Cooling wind brought into sensible Heat Q lq ;
Total furnace heat expenditure Qjq = flue gas carried heat Q yq (m r ) +product carry-over sensible heat Q cp Heat consumption Q of +firing sz Sensible heat Q is taken away by +discharging wind sc +heat loss Q ss ;
S2, acquiring real-time kiln production parameters; according to kiln production parameters, calculating sensible heat Q carried by a blank in a kiln heat balance mathematical model pt Sensible heat Q is brought in by cooling air lq Sensible heat Q of product cp Heat consumption Q of firing sz Sensible heat Q carried away by furnace outlet wind sc Heat loss Q ss The method comprises the steps of carrying out a first treatment on the surface of the Based on the numerical balance of the total heat income Qsr of the kiln and the total heat expenditure Qjq of the kiln, obtaining the fuel consumption flow m r And electric heating real-time power W d Is a value interval and constraint relation;
s3, matching the green electric power W which can be provided in real time under the condition that the carbon emission Te is less than or equal to the upper limit of carbon emission l The fuel consumption flow rate m obtained in step S2 r And electric heating real-time power W d Further limiting the fuel consumption flow m in the value interval and the constraint relation of (1) r And electric power W of ash h Is a value interval and constraint relation;
s4, fuel consumption flow m obtained according to the step S3 r And electric power W of ash h Matching real-time energy cost data, minimizing total heating energy cost Ae as optimization target, and determining fuel type and fuel consumption flow m r And electric power W of ash h Is a value of (2);
s5, according to the fuel type and the fuel consumption flow m r And electric power W of ash h To determine the sensible heat Q carried by the fuel in the mathematical model of the kiln heat balance x (m r ) Heat release Q of fuel combustion r (m r ) The combustion air is brought into sensible heat Q zk (m r ) Heat quantity Q is taken away by flue gas yq (m r ) To adjust kiln operation;
and after the kiln runs for a set time, repeating the steps S2 to S5 until the control is finished.
Preferably, the fuel is brought into sensible heat Q x (m r ) The method comprises the following steps: q (Q) x (m r )=m r ×C r ×t r ;
Wherein C is r Specific heat of fuel, determined by the type of fuel; t is t r The temperature of the fuel entering the kiln is detected by a temperature measuring module at the fuel inlet;
the fuel burns and releases heat Q r (m r ) The method comprises the following steps: q (Q) r (m r )=m r ×q r ;
Wherein q r Is the lower calorific value of the fuel and is determined by the type of the fuel;
the combustion air brings sensible heat Q zk (m r ) The method comprises the following steps: q (Q) zk (m r )=V zk ×C zk ×t zk ;
Wherein V is zk For combustion air flow, the combustion fan is used for controlling the fuel type and the fuel consumption flow m r Determining; t is t zk The temperature of the combustion air is detected by a temperature measuring module in the combustion fan; c (C) zk Specific heat for combustion supporting;
the electric heating quantity Q d (W d ) The method comprises the following steps: q (Q) d (W d )=W d ×ξ d ;
Wherein, xi d For the electric heating efficiency, set by the user;
the flue gas takes away heat Q yq (m r ) The method comprises the following steps: q (Q) yq (m r ) =V yq ×C yq ×t yq ;
Wherein V is yq The flow of the exhaust gas is regulated and controlled by a smoke exhaust fan, and the flow m of the fuel consumption is regulated and controlled by the smoke exhaust fan r And kiln internal temperature determination; t is t yq The temperature of the exhaust fume is detected by a temperature measuring module in the exhaust fume pipe; c (C) yq Specific heat of discharged smoke.
Preferably, the blank is brought into sensible heat Q pt The method comprises the following steps: q (Q) pt =m pt ×t pt [(1-W pt )C pt +4.18W pt ];
Wherein m is pt The kiln feeding amount of the blank body in unit time is shown; t is t pt The kiln feeding temperature of the green body; w (W) pt The water content of the green body is; c (C) pt To obtain specific heat at kiln inlet temperature, the green body is subjected to kiln inlet temperature t pt The calculation results are that: c (C) pt =0.836+2.63×10 -4 ×t pt ;
The product takes out sensible heat Q cp The method comprises the following steps: q (Q) cp =m cp ×t cp ×C cp ;
Wherein m is cp Yield per unit time; t is t cp The kiln outlet temperature of the product is detected by a temperature measuring module at the kiln outlet; c (C) cp To obtain specific heat at kiln outlet temperature, the kiln outlet temperature t of the product cp The calculation results are that: c (C) cp =0.836+2.63×10 -4 ×t cp 。
Preferably, the cooling wind is brought into sensible heat Q lq The method comprises the following steps:
Q lq =Q jl +Q wl ;
wherein Q is jl Represents the sensible heat brought in by the quenching wind; q (Q) wl Represents sensible heat brought in by tail cold air;
sensible heat Q is taken away by the furnace outlet wind sc The method comprises the following steps:
Q sc =Q yr +Q wc ;
wherein Q is yr Represents that the sensible heat is taken away by the hot air; q (Q) wc And represents the sensible heat taken away by the hot air of the kiln tail chimney.
Preferably, the quench wind brings in sensible heat Q jl Sensible heat Q is brought in by tail cold air wl Sensible heat Q is taken away by hot air pumping yr Sensible heat Q is taken away to kiln tail chimney hot air wc The calculation method of (1) is as follows:
Q i =V i ×C i ×t i ,i=jl,wl,yr,wc;
wherein V is i Is the wind flow; t is t i Is the wind temperature; c (C) i Specific heat of wind, from the temperature t of wind i The calculation results are that: c (C) i =1.3+7.45×10 -6 t i +5.95×10 -8 t i 2 。
Preferably, the firing heat consumption Q sz The method comprises the following steps: q (Q) sz =green body moisture evaporation and heat consumption Q zs Decomposition of Clay Heat consumption Q during the+ firing process h ;
The moisture of the green body evaporates and heats up the heat consumption Q zs The method comprises the following steps: q (Q) zs =m zs ×(2490+1.93t y );
Wherein m is zs As the total water quantity, m zs =adsorbed water mass m xs Water amount of+ crystal m js ;m xs =m pt ×W pt ;m js =m pt ×(1-W pt )-m gp (MgO×44/40+CaO×44/56) ;m gp For dry blank mass, m gp =m pt -m zs ;t y The temperature of the flue gas leaving the kiln is detected by a temperature measuring module in the smoke exhaust pipe;
the firing process decomposes the clay heat consumption Q h The method comprises the following steps: q (Q) h =2.1×m gp ×Al 2 O 3 ;
Wherein Al is 2 O 3 Is oxidized for the formulaAluminum content.
Preferably, the heat loss Q ss The method comprises the following steps:
Q ss =incomplete combustion loss Q bx Surface heat loss Q of kiln bs +other heat loss Q t 。
Preferably, the carbon emission Te is:
Te=m r ×τ r +W h ×τ d ;
wherein τ r Represents the carbon-containing factor of the fuel gas; τ d Represents an gray electrical carbon factor; w (W) h =W d -W l 。
Preferably, the real-time energy cost data includes an ash electricity unit price a d Price aτ of carbon tax and unit price a of fuel gas r ;
The total cost Ae of the heating energy is as follows: ae=a d ×W h +τ d ×W h ×aτ+m r ×τ r ×aτ+m r ×a r ;
Wherein τ r Represents the carbon-containing factor of the fuel gas; τ d Representing the gray electrical carbon factor.
Preferably, the step S5 further comprises, depending on the fuel type, the fuel consumption flow rate m r And electric power W of ash h Calculating the carbon emission Te; the carbon emission of kiln combustion can be tracked.
Compared with the prior art, the invention has the following advantages and beneficial effects:
the above-mentioned carbon emissions define the use interval of fuel gas and ash electricity. The heating energy total cost has the lowest value in the interval, so the method can find the optimal economic energy collocation scheme under the carbon emission standard
1. The method of the invention constructs a mathematical model of kiln heat balance based on kiln heat income expenditure balance, combines carbon emission index and economic index to obtain energy source type proportion and corresponding input quantity, can meet carbon emission requirement, reduce kiln use cost, lead combustion to reach environment-friendly state of optimal economy, can quickly realize energy source type change, stabilize kiln internal temperature and ensure product firing quality;
2. the method can calculate the steady-state data of kiln parameters such as flue gas, combustion air, fuel gas and the like through the mathematical model of kiln heat balance when the kiln is subjected to optimal economy adjustment of combustion, thereby correspondingly adjusting the opening of a valve element of the fan and leading the heat balance in the kiln to be in dynamic balance.
Drawings
FIG. 1 is a flow chart of the intelligent kiln control method for multi-energy combustion carbon emission tracking of the invention.
Detailed Description
The invention is described in further detail below with reference to the drawings and the detailed description.
Examples
The intelligent kiln control method for tracking the carbon emission of multi-energy combustion is used for kilns comprehensively using fuel, green electricity and ash electricity as energy sources.
The intelligent kiln control method for tracking the carbon emission of the multi-energy combustion, as shown in fig. 1, comprises the following steps:
s1, constructing a kiln heat balance mathematical model with self-adaptive parameters; the adaptive parameters include fuel type, fuel consumption flow m r And electric heating real-time power W d The method comprises the steps of carrying out a first treatment on the surface of the The kiln heat balance mathematical model is as follows:
total kiln heat input Qsr = fuel carry-in sensible heat Q x (m r ) +Fuel Combustion exotherm Q r (m r ) +combustion air brought into sensible heat Q zk (m r ) +electric heating quantity Q d (W d ) +green body brought into sensible heat Q pt +Cooling wind brought into sensible Heat Q lq ;
Total furnace heat expenditure Qjq = flue gas carried heat Q yq (m r ) +product carry-over sensible heat Q cp Heat consumption Q of +firing sz Sensible heat Q is taken away by +discharging wind sc +heat loss Q ss ;
The calculation mode of each numerical value is as follows:
(1) Sensible heat Q of fuel x (m r ) The method comprises the following steps: q (Q) x (m r )=m r ×C r ×t r ;
Wherein, the liquid crystal display device comprises a liquid crystal display device,C r specific heat of fuel, determined by the type of fuel; t is t r The temperature of the fuel entering the kiln is detected by a temperature measuring module at the fuel inlet;
(2) Fuel combustion exotherm Q r (m r ) The method comprises the following steps: q (Q) r (m r )=m r ×q r ;
Wherein q r Is the lower calorific value of the fuel and is determined by the type of the fuel;
(3) Combustion air is brought into sensible heat Q zk (m r ) The method comprises the following steps: q (Q) zk (m r )=V zk ×C zk ×t zk ;
Wherein V is zk For combustion air flow, the combustion fan is used for controlling the fuel type and the fuel consumption flow m r Determining; t is t zk The temperature of the combustion air is detected by a temperature measuring module in the combustion fan; c (C) zk To aid specific heat, the temperature t of the aid air zk And (3) calculating to obtain: 1.3+7.45X10 -6 t zk +5.95×10 -8 t zk 2 ;
(4) Electric heating quantity Q d (W d ) The method comprises the following steps: q (Q) d (W d )=W d ×ξ d ;
Wherein, xi d Is the electric heating efficiency;
(5) The green body is brought into sensible heat Q pt The method comprises the following steps: q (Q) pt =m pt ×t pt [(1-W pt )C pt +4.18W pt ];
Wherein m is pt The kiln feeding amount of the blank body in unit time is shown; t is t pt The kiln feeding temperature of the green body; w (W) pt The water content of the green body is; c (C) pt To obtain specific heat at kiln inlet temperature, the green body is subjected to kiln inlet temperature t pt The calculation results are that: c (C) pt =0.836+2.63×10 -4 ×t pt ;
(6) Cooling air is brought into sensible heat Q lq The method comprises the following steps: q (Q) lq =Q jl +Q wl ;
Wherein Q is jl Represents the sensible heat brought in by the quenching wind; q (Q) wl Represents sensible heat brought in by tail cold air;
quenching wind bringing in sensible heat Q jl The method comprises the following steps: q (Q) jl =V jl ×C jl ×t jl ;
Wherein V is jl The flow of the rapid cooling air is detected by a rapid cooling fan flow detection module; t is t jl The temperature is the temperature of the rapid cooling air, and is detected by a temperature measuring module in the rapid cooling fan; c (C) jl To quench the specific heat of wind, from the temperature t of the quench wind jl The calculation results are that: c (C) jl =1.3+7.45×10 -6 t jl +5.95×10 -8 t jl 2 ;
Tail cooling wind brought into sensible heat Q wl The method comprises the following steps: q (Q) wl =V wl ×C wl ×t wl ;
Wherein V is wl The tail cold air flow is detected by a tail air cooler flow detection module; t is t wl The temperature of the tail cold air is detected by a temperature measuring module in the tail air cooler; c (C) wl Specific heat of the tail cold air, from the temperature t of the tail cold air wl The calculation results are that: c (C) wl =1.3+7.45×10 -6 t wl +5.95×10 -8 t wl 2 ;
(7) Heat Q is taken away by flue gas yq (m r ) The method comprises the following steps: q (Q) yq (m r ) =V yq ×C yq ×t yq ;
Wherein V is yq The flow of the exhaust gas is regulated and controlled by a smoke exhaust fan, and the flow m of the fuel consumption is regulated and controlled by the smoke exhaust fan r And kiln internal temperature determination; t is t yq The temperature of the exhaust fume is detected by a temperature measuring module in the exhaust fume pipe; c (C) yq To discharge specific heat, the temperature t of discharged smoke yq And (3) calculating to obtain: c (C) yq =1.3+7.45×10 -6 t yq +5.95×10 -8 t yq 2 ;
(8) Sensible heat Q of product cp The method comprises the following steps: q (Q) cp =m cp ×t cp ×C cp ;
Wherein m is cp Yield per unit time; t is t cp The kiln outlet temperature of the product is detected by a temperature measuring module at the kiln outlet; c (C) cp To obtain specific heat at kiln outlet temperature, the kiln outlet temperature t of the product cp The calculation results are that: c (C) cp =0.836+2.63×10 -4 ×t cp ;
(9) Burning outHeating heat Q sz The method comprises the following steps: q (Q) sz =green body moisture evaporation and heat consumption Q zs Decomposition of Clay Heat consumption Q during the+ firing process h ;
Moisture evaporation and heating heat consumption Q of green body zs The method comprises the following steps: q (Q) zs =m zs ×(2490+1.93t y );
Wherein m is zs As the total water quantity, m zs =adsorbed water mass m xs Water amount of+ crystal m js ;m xs =m pt ×W pt ;m js =m pt ×(1-W pt )-m gp (MgO×44/40+CaO×44/56) ;m gp For dry blank mass, m gp =m pt -m zs ;t y The temperature of the flue gas leaving the kiln is detected by a temperature measuring module in the smoke exhaust pipe;
decomposition of clay heat consumption Q during firing h The method comprises the following steps: q (Q) h =2.1×m gp ×Al 2 O 3 ;
Wherein Al is 2 O 3 Is the content of the formula alumina;
(10) Sensible heat Q is taken away by furnace outlet wind sc The method comprises the following steps: q (Q) sc =Q yr +Q wc ;
Wherein Q is yr Represents that the sensible heat is taken away by the hot air; q (Q) wc Sensible heat is taken away by hot air representing a kiln tail chimney;
sensible heat Q is taken away by hot air pumping yr The method comprises the following steps: q (Q) yr =Vyr×Cyr×tyr;
Wherein Vyr is the flow of the hot air, and is detected by a flow detection module of the hot air pipe; tyr is the temperature of the hot air, and is detected by a temperature measuring module in the hot air pipe; cyr is the specific heat of the hot air, and is calculated by the temperature tyr of the hot air: cyr=1.3+7.45×10 -6 tyr +5.95×10 -8 tyr 2 ;
Sensible heat Q is taken away to kiln tail chimney hot air wc The method comprises the following steps: q (Q) wc =V wc ×C wc ×t wc ;
Wherein V is wc The tail exhaust flow is detected by a tail exhaust fan flow detection module; t is t wc The temperature of the tail exhaust fan is detected by a temperature measuring module in the tail exhaust fan; c (C) wc Specific heat of tail draftDegree t wc The calculation results are that: c (C) wc =1.3+7.45×10 -6 t wc +5.95×10 -8 t wc 2 ;
(11) Heat loss Q ss The method comprises the following steps: q (Q) ss =incomplete combustion loss Q bx Surface heat loss Q of kiln bs +other heat loss Q t The method comprises the steps of carrying out a first treatment on the surface of the Incomplete combustion loss Q bx Surface heat loss Q of kiln bs Other heat losses Q t And a table can be established according to detection by setting by a user, and the value can be determined through table lookup.
S2, acquiring real-time kiln production parameters; according to kiln production parameters, calculating sensible heat Q carried by a blank in a kiln heat balance mathematical model pt Sensible heat Q is brought in by cooling air lq Sensible heat Q of product cp Heat consumption Q of firing sz Sensible heat Q carried away by furnace outlet wind sc Heat loss Q ss The method comprises the steps of carrying out a first treatment on the surface of the Based on the numerical balance of the total heat income Qsr of the kiln and the total heat expenditure Qjq of the kiln, obtaining the fuel consumption flow m r And electric heating real-time power W d Is a value interval and constraint relation;
s3, matching the green electric power W which can be provided in real time under the condition that the carbon emission Te is less than or equal to the upper limit of carbon emission l The fuel consumption flow rate m obtained in step S2 r And electric heating real-time power W d Further limiting the fuel consumption flow m in the value interval and the constraint relation of (1) r And electric power W of ash h Is a value interval and constraint relation; the carbon emission binary function Te is:
Te=m r ×τ r +W h ×τ d ;
wherein τ r Represents the carbon-containing factor of the fuel gas; τ d Represents an gray electrical carbon factor; w (W) h =W d -W l ;
S4, fuel consumption flow m obtained according to the step S3 r And electric power W of ash h Matching real-time energy data, minimizing total heating energy cost Ae as optimization target, and determining fuel type and fuel consumption flow m r And electric power W of ash h Is a value of (2);
the total heating energy cost Ae is: ae=a d ×W h +τ d ×W h ×aτ+m r ×τ r ×aτ+m r ×a r ;
Wherein a is d Is ash electricity unit price; aτ is a carbon tax price; a, a r Is the unit price of fuel gas;
s5, according to the fuel type and the fuel consumption flow m r And electric power W of ash h To determine the sensible heat Q carried by the fuel in the mathematical model of the kiln heat balance x (m r ) Heat release Q of fuel combustion r (m r ) The combustion air is brought into sensible heat Q zk (m r ) Heat quantity Q is taken away by flue gas yq (m r ) To adjust kiln operation;
specifically, the fuel type and the fuel consumption flow rate m are determined r Electric ash power W h Flow V of combustion-supporting air zk Flow V of exhaust gas yq Adjusting a fan and a valve corresponding to the kiln;
step S5 preferably further comprises: according to fuel type, fuel consumption flow m r And electric power W of ash h Calculating the carbon emission Te;
and after the kiln runs for a set time, repeating the steps S2 to S5 until the control is finished.
Because each fan and valve element of the kiln can change according to the running state in the running process of the kiln, after the kiln runs for a set time, the steps S2 to S5 are repeatedly executed, so that the balance of the heat income and the heat expenditure of the kiln can be kept in the running process of the kiln, and the temperature in the kiln is stabilized.
The above examples are preferred embodiments of the present invention, but the embodiments of the present invention are not limited to the above examples, and any other changes, modifications, substitutions, combinations, and simplifications that do not depart from the spirit and principle of the present invention should be made in the equivalent manner, and the embodiments are included in the protection scope of the present invention.