The present invention relates to the technical field of blast furnace ironmaking processes, in particular to a calculating and real-time monitoring method for a boundary of a blast furnace tuyere raceway.
During blast furnace ironmaking production, high-temperature high-speed air is blown into a blast furnace through a blast furnace tuyere. Due to the effect of blast wind, a region in which coke circulates is formed near the front of a tuyere, which is called a blast furnace tuyere raceway. The tuyere raceway is located at the front end of the lower tuyere of a vertical furnace body, and is formed by an intense combustion reaction between coke and auxiliary fuel injected into the furnace and oxygen in the blast wind. Mixed airflow of the blast wind and gas performs circulating motion in the region, accompanied by high-speed rotation of small coke granules and unburned coal powder, as well as a burnout process of crushed coke in a swirling process. The size of the raceway is directly related to blast wind parameters, raw fuel conditions and other factors. It is the birthplace of heat energy and a gas reducing agent, providing heat and energy supply for the whole blast furnace production. The depth of the tuyere raceway and complex physical and chemical reactions inside determine the primary distribution of coal gas flow in the blast furnace and the declining state of upper furnace charge, reflect the burning state of the coke, is the basis for smooth operation of furnace conditions, and plays a crucial role in the smelting process.
The study on the characteristics of the blast furnace tuyere raceway can be mainly divided into two aspects: direct method study on the characteristics of the raceway and indirect method study on the characteristics of the raceway. Firstly, for the direct study method on the characteristics of the blast furnace tuyere raceway, study is carried out by directly detecting the relevant parameters that represent the blast furnace raceway, mainly focusing on the direct measurement of parameters such as the size, shape, and temperature of the raceway. However, there are problems of significant fluctuations in measurement results caused by the situation that instrument equipment is easy to affect by the actual environment inside the furnace, high instrument cost, and inability to achieve real-time monitoring, so that the direct study method cannot be fully popularized in small and medium enterprises. The direct study method can be divided into study by an empirical observation method and study by an actual measurement method. Secondly, the indirect study method on the characteristics of the blast furnace raceway, namely a model study method, comprises the following two aspects: firstly, by establishing a physical parameter experimental model of the blast furnace tuyere raceway, experimental detection is carried out on the model based on the characteristics of the blast furnace tuyere raceway. However, due to the internal presence of the raceway but the complex and variable reactions inside the raceway, a cold-state model cannot reflect the actual internal state of the raceway well. A commonly used method is to establish an Euler mathematical model for solving based on the transfer of momentum, mass and heat during the motion of the raceway. However, if existing Euler models are adopted, modelling is complex, a large number of parameters are required, calculation is difficult, a long time is consumed, the purpose of real-time monitoring of reactions is difficult to achieve, and experimental models cannot reflect the actual internal state of the raceway well. Secondly, an established two-dimensional or three-dimensional mathematical model for the blast furnace tuyere raceway is used for numerically simulating the chemical reaction process in the raceway, so as to achieve the purpose of studying the characteristic change law of the raceway.
A technical problem to be solved by the present invention is to provide a calculating and real-time monitoring method for a boundary of a blast furnace tuyere raceway so as to address the shortcomings of a direct measurement method and an experimental model method mentioned above, and to improve a mechanism mathematical model. On the basis of a depth model for a raceway, by establishing an equilibrium equation of two forces at each point at the boundary of the raceway, the changes in the boundary of the raceway can be solved efficiently in real time, the change law of the boundary of the raceway is obtained, the impact of internal parameters of the raceway on the depth and height of the raceway is studied, the changes in the depth and height of the raceway are adjusted by controlling blast wind parameters, and reliable guarantees are provided for the stable operation of the blast furnace.
To solve the above-mentioned technical problems, the technical solution adopted by the present invention is to provide a calculating and real-time monitoring method for the boundary of a blast furnace tuyere raceway, comprising the following steps.
Step 1: establishing a depth calculation model for the raceway according to a formation principle of the blast furnace tuyere raceway, and further obtaining a calculation formula for a depth of the raceway so as to obtain a change law of the depth of the raceway.
When a cavity in the raceway is in a stable motion state, taking a micro element region A at a deepest position inside the raceway as a research object, at this time, enabling the region A to be balanced under a combined action of a blast wind gas impulse and a resistance of a coke layer, and establishing the depth calculation model for the raceway according to an equilibrium of two forces to solve changes of the depth of the raceway, as shown in the following formulas:
Wherein FA represents the blast wind gas impulse at the region A, FB represents the resistance of the coke layer at the region A, ρg0 represents a blast wind density in a standard state, Vg represents the blast wind volume at the region A, ST represents an area of a tuyere, DR represents the depth of the raceway, DT represents a diameter of the tuyere, α is a constant which is used for representing a relationship between the depth and a width of the raceway, P represents the blast wind pressure, Tm represents a temperature of the raceway, Sp represents a total cross-sectional area of coke particles at the region A, pp represents a coke density, Vp represents a volume of the coke particles at the region A, g represents an acceleration of gravity, DPR represents a diameter of coke in front of the boundary of the raceway, DPR=0.6DPC, DPC represents a diameter of the coke before entering a furnace, and Tw represents blast wind temperature.
The calculation formula for the depth of the raceway is obtained through the depth calculation model for the raceway:
Wherein K and β both are coefficients to be determined, and ρg represents a blast wind density under an actual wind temperature and wind pressure.
Step 2: establishing a boundary model for the blast furnace tuyere raceway through the depth calculation model for the blast furnace tuyere raceway, and determining a calculation formula for the boundary of the raceway.
When an internal motion of the raceway is in a stable state, randomly taking a point B at the boundary of the raceway as a research object, at this time, enabling the point B to be balanced under the combined action of the blast wind gas impulse and the resistance of the coke layer, and establishing a mathematical model for any point at the boundary of the raceway according to the equilibrium of two forces to solve changes of the boundary of the raceway.
FD represents a blast wind gas impulse at the point B, Fb represents a resistance of the coke layer at the point B to a blast wind airflow, MB represents a blast wind mass flow rate at the point B, UB represents a wind speed at the point B, SB represents a cross-sectional area of the point B, SA represents a cross-sectional area of the region A, ρg represents the blast wind density under the actual wind temperature and wind pressure, P0 represents 1 standard atmospheric pressure, L represents a curve distance from the deepest position of the raceway to any point on the boundary of the raceway, σ represents a loss rate of blast wind on the boundary of the raceway, H represents a total height of a blast furnace, h represents a vertical distance from the point B to the tuyere, a and b are shape parameters of the boundary of the raceway, I represents a total number of coke granules included in the cross-section at the point B, and DPRi is a diameter of an ith coke granule at the point B of the boundary of the raceway.
Because the boundary from the deepest position of the raceway to any point on the boundary of the raceway is in a curved shape, L=(DR−x)a+hb is set, and wherein coordinates of the point B of the boundary of the raceway relative to the tuyere are (x, h).
The calculation formula for the boundary of the raceway is further obtained through the mathematical model for any point at the boundary of the raceway.
K4 and K4 are coefficients to be determined.
Step 3: obtaining modelling parameters, analyzing an impact of the modelling parameters on the boundary model for the raceway, and determining main parameters that affect the boundary of the raceway.
According to the boundary model for the raceway and relevant parameters involved in modelling, and a change law of the boundary of the raceway, the blast wind pressure P and the blast wind volume Vg in blast wind parameters, as well as the temperature Tm of the raceway and the loss rate σ of the blast wind at the coke of the raceway in the blast wind parameters, have a direct impact on the changes of the boundary of the raceway; according to the calculation formula for the boundary of the raceway, reducing the blast wind pressure P, increasing the blast wind volume Vg and increasing the temperature Tm and reducing the loss rate of the blast wind at the coke of the raceway are all beneficial for developing the raceway towards a center and increasing the depth of the raceway, wherein the blast wind pressure P and the blast wind volume Vg are adjustable, while changes of the temperature Tm in the raceway are influenced by various factors, whereby the temperature is not directly adjustable to cause changes.
Step 4: solving a height of the raceway through the calculation formula for the boundary of the raceway.
Deriving a variable h in the calculation formula for the boundary of the raceway, and making the derived result equal to zero, thereby calculating a maximum value of the raceway in a vertical direction, that is, the height of the raceway, as shown in the following formula:
Wherein G is the height of the raceway, and K5 is a coefficient to be determined.
Step 5: when the height or the depth of the raceway exceeds a set range, adjusting a blast wind pressure P and a blast wind volume Vg to restore the height or the depth of the raceway to a normal range.
Beneficial effects generated by adopting the above technical solution lie in that in the calculating and real-time monitoring method for the boundary of a blast furnace tuyere raceway provided by the present invention, firstly, the depth model for the raceway is determined, based on the depth model, the mathematical model for any point on the boundary of the raceway is deduced through physical principles and chemical principles, and through the mathematical model, the pocket shape of the blast furnace tuyere raceway can be intuitively reflected; the variable of temperature of the raceway is successfully introduced, two blast wind parameters that affect the depth and boundary of the raceway, namely blast wind pressure and blast wind temperature are obtained, the impact of the two parameters on the changes in the depth and boundary of the raceway is studied, and when there are unreasonable changes in the depth and boundary of the raceway, the two parameters are adjusted timely to restore the depth and boundary of the raceway to the normal range. According to the method disclosed by the present invention, the changes of the depth and boundary of the raceway can be monitored in real-time, providing safety guidance for the actual production of blast furnaces so as to ensure the safety of life and property of enterprises.
The specific implementations of the present invention are described in more detail below with reference to the accompanying drawings and embodiments. The following embodiments are intended to illustrate the present invention, rather than to limit the scope of the present invention.
The formation of the blast furnace tuyere raceway is shown in
In the embodiments of the present invention, based on the depth calculation model for the raceway, a research object is any point B on the boundary of the raceway. By analyzing the force at the point B and comparing the force with that in the deepest point, a boundary model of the raceway is established, and monitoring the boundary of the raceway is carried out based on the boundary model of the raceway.
A calculating and real-time monitoring method for each point on the boundary of a blast furnace tuyere raceway comprises the following steps.
Step 1: establishing a depth calculation model for the raceway according to a formation principle of the blast furnace tuyere raceway, as shown in
When an internal motion of the raceway is in a stable state, taking a region A at a deepest position inside the raceway as a research object, at this time, enabling the region A to be balanced under a combined action of a blast wind gas impulse and a resistance of a coke layer, and establishing the depth calculation model for the raceway according to an equilibrium of two forces to solve changes of the depth of the raceway, as shown in the following formula:
Wherein FA represents the blast wind gas impulse at the region A, FB represents the resistance of the coke layer at the region A, ρg0 represents a blast wind density in a standard state, Vg represents the blast wind volume at the region A, ST represents an area of the tuyere, DR represents the depth of the raceway, DT represents a diameter of the tuyere, a is a constant which is used for representing a relationship between the depth and a width of the raceway, P represents the blast wind pressure, Tm represents a temperature of the raceway, measured by a CCD temperature measuring instrument, SP represents the total cross-sectional area of coke particles at the region A, ρp represents a coke density, VP represents a volume of the coke particles at the region A, g represents an acceleration of gravity, DPR represents a diameter of coke in front of the boundary of the raceway, DPR=0.6DPC, DPC represents a diameter of the coke before entering the furnace, ρg represents blast wind density under actual wind temperature and wind pressure, and Tw represents a blast wind temperature.
The calculation formula for the depth of the raceway is obtained through the depth calculation model for the raceway:
Wherein K and β are coefficients to be determined.
After extensive experimental verification, the calculation formula for the depth of the raceway is solved finally in the embodiment.
Step 2: establishing a boundary model for the blast furnace tuyere raceway through the depth calculation model for the blast furnace tuyere raceway, and determining a calculation formula for the boundary of the raceway.
When blast wind in the blast furnace tuyere raceway blows from the tuyere to the deepest position, the blast wind will continue to move up and down until it reaches the highest point in the tuyere raceway; therefore, the model can be derived by using the deepest position in the blast furnace tuyere raceway as a source of upward blast wind; the blast wind starts to blow up from the deepest position in the tuyere raceway and blows out of the cavity in the tuyere raceway until it reaches the highest point in the raceway; on the basis of an established formula of the depth calculation model for the blast furnace tuyere raceway, and the boundary model for each point at the boundary of the blast furnace tuyere raceway can be derived based on the principle of an equilibrium of forces due to the loss of the blast wind at the boundary of the raceway.
After being blown therein from the tuyere, high-temperature blast wind reaches the deepest point A in the raceway in the stroke being approximately a conical pipeline, and then the blast wind moves upwards; a micro element region B is formed at a horizontal distance of x meters from the raceway tuyere and a vertical distance of h meters from the raceway tuyere, with a cross-sectional area of SB. At a point B, under a combined action of a blast wind gas impulse FD and a resistance Fb of the coke layer at the point B, a balance state is achieved; at this time, it is considered that the raceway is in a stable operating state and satisfies an expression:
Wherein the magnitude of the blast wind gas impulse FD is mainly proportional to a blast wind mass and a blast wind speed, namely:
In the expression, MB represents a blast wind mass flow rate at the position, which specifically refers to a mass of fluid passing through a cross-section per unit time, namely:
Wherein ρg represents a blast wind density, Vg represents a blast wind volume, L represents a curve distance from the deepest position of the raceway to any point on the boundary of the raceway, and σ represents a loss rate of the blast wind on the boundary of the raceway.
L is further discussed, because the boundary from the deepest position of the raceway to any point on the boundary of the raceway is in a curved shape, assuming L is equal to (DR−x)a+hb, wherein (x, h) is the coordinate of the point B relative to the tuyere, and a and b are curve parameters.
UB represents a wind speed at the point B, which is proportional to the blast wind pressure and also directly proportional to the actual temperature inside the raceway, namely:
In the expression, MB represents a blast wind mass flow rate at the point B, ρg represents the blast wind density, SB represents a cross-sectional area of the point B, P represents the blast wind pressure, P0 represents 1 standard atmospheric pressure, Tm represents the temperature of the raceway, actually measured using a temperature measuring instrument, and T0 is equal to 298K.
SB is the cross-sectional area of the point B, with the expressions:
In the expression, WR represents a width of the raceway.
Sp is a total cross-sectional area of coke particles at the point B, with the expression:
In the expression, DPR represents a diameter of coke in front of the boundary of the raceway, DP, represents a diameter of the coke before entering the furnace, DPR=0.6DPC.
The expression of the blast wind gas impulse FD at the point B can be further deduced as follows through combining equations (3) and (4):
The experiment proves that the relationship between the depth DR of the raceway and the width WR of the raceway is:
In the equation, α is a constant, and α is equal to 0.2 to 0.6;
ST is set to represent the area of the tuyere, with the expression:
By combining equation (5), equation (8) and equation (9), it can be inferred:
Finally, the final expression of the blast wind gas impulse FD at the point B can be obtained as follows through substituting the equation (10) into the equation (7):
Fb represents the resistance of the coke layer at the point B to a blast wind airflow, which is mainly determined by a compaction degree of the coke, and is approximately proportional to the gravity of coke granules and the blast wind pressure, namely:
In the expression, ρp represents the coke density, VP represents the volume of the coke particles at the point A, g represents the acceleration of gravity, DPR represents the diameter of the coke in front of the boundary of the raceway, P represents the blast wind pressure, H represents a total height of the blast furnace, and h represents a vertical distance from the point B to the tuyere; the blast wind density ρg can be further expressed as:
In the expression, ρg0 represents the blast wind density under the standard state, Tw represents the blast wind temperature, P represents the blast wind pressure, and P0 represents 1 standard atmospheric pressure.
Finally, by substituting the expressions (11), (12) and (13) into the expression (1), the calculation formula for the boundary of the raceway can be obtained as follows:
Wherein K3 and K4 are coefficients to be determined.
By using the calculation formula for the boundary of the raceway mentioned above, a mathematical model of the raceway can be further established to study the trend of changes in the boundary of the raceway.
Step 3: obtaining modelling parameters, analyzing an impact of the modelling parameters on the boundary model for the raceway, and determining main parameters that affect the boundary of the raceway.
According to the boundary model for the raceway as shown in
In the embodiment, the relevant parameters used in the modelling process of the boundary model of the raceway are shown in Table 1:
Step 4: solving a height of the raceway through the calculation formula for the boundary of the raceway.
Deriving a variable h in the calculation formula for the boundary of the raceway, and making the derived result equal to zero, thereby calculating a maximum value of the raceway in a vertical direction, that is, the height of the raceway, as shown in the following formula:
Wherein G is the height of the raceway, and K5 is a coefficient to be determined.
After extensive experimental verification, finally K5 being equal to 2.12 is obtained, and the calculation formula for the height of the raceway is:
Step 5: when the height or the depth of the raceway exceeds a set range, adjusting the blast wind pressure P and the blast wind volume Vg to restore the height or the depth of the raceway to a normal range.
In the embodiment, the changes in the depth and height of the raceway under smooth operation of the blast furnace are shown as
Finally, it should be noted that the above embodiments are merely intended to describe the technical solutions of the present invention, rather than to limit the present invention. Although the present invention is described in detail with reference to the above embodiments, those ordinary skilled in the art should understand that they may still make modifications to the technical solutions described in the above embodiments or make equivalent replacements to some or all technical features thereof. However, these modifications or replacements do not cause the essence of the corresponding technical solutions to depart from the scope of the technical solutions of the embodiments of the present invention.
Number | Date | Country | Kind |
---|---|---|---|
202110566801.4 | May 2021 | CN | national |
Filing Document | Filing Date | Country | Kind |
---|---|---|---|
PCT/CN2021/135430 | 12/3/2021 | WO |