Method, device, system, equipment and medium for determining axial temperature distribution in blast furnace

Abstract

The present application provides a method, device, system, equipment and medium for determining the axial temperature distribution of a blast furnace. The method includes: obtaining smelting data during the smelting of the blast furnace; calculating the theoretical combustion temperature in the tuyere raceway in the blast furnace, as well as the first gas-solid temperature distribution in the dripping zone, softening-melting zone and lump zone in the blast furnace and the first indirect reduction degree of blast furnace ironmaking based on the smelting data; generating a first blast furnace axial temperature distribution curve; judging the first blast furnace axial temperature distribution curve according to a preset convergence condition; if the first blast furnace axial temperature distribution curve meets the preset convergence condition, judging whether the sum of the first indirect reduction degree and the direct reduction degree is a preset threshold; if the sum of the first indirect reduction degree and the direct reduction degree is the preset threshold, determining the axial temperature distribution of the blast furnace according to the first blast furnace axial temperature distribution curve. By determining the axial temperature distribution during the smelting of the blast furnace, the heat state at different height positions in the blast furnace can be obtained, providing a guarantee for the stable operation of ironmaking.

Description

Technical Field

[0001] The present invention relates to the technical field of blast furnace smelting, and particularly to a method, device, system, equipment and medium for determining the axial temperature distribution of a blast furnace. Background Art

[0002] Blast furnace smelting is the main way of modern ironmaking production, which is a production process of physics, chemistry and kinetics under severe conditions such as high temperature, high pressure and airtightness. As a clean energy, hydrogen produces pollution-free water vapor when reducing iron oxides, and "replacing carbon with hydrogen" has attracted wide attention in the ironmaking industry and has become one of the key directions for the low-carbon development of blast furnaces. However, due to the endothermic characteristics of hydrogen during the reduction process, the axial temperature distribution in the blast furnace has changed significantly compared with that of conventional blast furnaces, and key technical common problems such as "cold at the top and hot at the bottom" have emerged. If the axial temperature in the furnace cannot be monitored, it is impossible to immediately adjust the furnace condition to be cooler or hotter, and even make wrong adjustments, resulting in the deterioration of the furnace condition. Therefore, real-time monitoring of the axial temperature distribution in the blast furnace is crucial for ensuring the stable and smooth operation of low-carbon blast furnace smelting.

[0003] The traditional axial temperature distribution in the blast furnace is determined by real-time monitoring of the basic conditions of the blast furnace temperature. The basic conditions of the blast furnace temperature can reflect the temperature distribution law and its changes on the inner edge furnace wall of the blast furnace, so as to monitor the temperature changes in the furnace. However, due to the high temperature, high pressure and corrosive substances inside the blast furnace, in production, it is impossible to finely and effectively monitor the axial temperature distribution of the blast furnace in actual production according to the basic conditions of the blast furnace temperature, and it is also impossible to better propose specific methods for effectively maintaining the efficient operation of the blast furnace at the theoretical and technical levels, which affects the stable operation of blast furnace smelting. Summary of the Invention

[0004] In view of the above-mentioned disadvantages of the prior art, the purpose of the present invention is to provide a method, device, system, equipment and medium for determining the axial temperature distribution of a blast furnace, which is used to solve the technical problem that the axial temperature distribution in the blast furnace cannot be accurately monitored in real time during blast furnace smelting production.

[0005] In view of the above problems, the present invention provides a method for determining the axial temperature distribution of a blast furnace, and the method includes:

[0006] Obtain the smelting data during the blast furnace smelting process, where the smelting data includes the direct reduction degree of blast furnace ironmaking;

[0007] Based on the smelting data, calculate the theoretical combustion temperature in the tuyere raceway area in the blast furnace, the first gas-solid temperature distribution in the dripping zone, softening-melting zone and lump zone in the blast furnace, and the first indirect reduction degree of blast furnace ironmaking respectively;

[0008] Generate a first blast furnace axial temperature distribution curve based on the theoretical combustion temperature, the first gas-solid temperature distribution in the dripping zone, the first gas-solid temperature distribution in the softening-melting zone, and the first gas-solid temperature distribution in the lump zone;

[0009] Judge the first blast furnace axial temperature distribution curve according to the preset convergence condition;

[0010] If the first blast furnace axial temperature distribution curve meets the preset convergence condition, judge whether the sum of the first indirect reduction degree and the direct reduction degree is a preset threshold;

[0011] If the sum of the first indirect reduction degree and the direct reduction degree is the preset threshold, determine the blast furnace axial temperature distribution according to the first blast furnace axial temperature distribution curve.

[0012] In an embodiment of the present invention, based on the smelting data, the gas volume and gas composition in the tuyere raceway, the gas volume and gas composition in the dripping zone, the gas volume and gas composition in the softening-melting zone, and the gas volume and gas composition in the lump zone of the tuyere raceway are also calculated and obtained.

[0013] In an embodiment of the present invention, the obtaining of the theoretical combustion temperature, the gas volume and gas composition in the tuyere raceway includes:

[0014] Based on the smelting data, obtain the raw fuel composition and consumption, the blast volume, the blast temperature, the blast oxygen enrichment rate, the gas injection volume, the total heat released during the combustion processes of coke, pulverized coal and injected gas in the tuyere raceway, and the specific heat capacity of the gas in the tuyere raceway;

[0015] Calculate the raw fuel composition and consumption, the blast volume, the blast temperature, the blast oxygen enrichment rate, the gas injection volume, the total heat, and the specific heat capacity of the gas in the tuyere raceway to obtain the theoretical combustion temperature, the gas volume and gas composition in the tuyere raceway.

[0016] In an embodiment of the present invention, the obtaining of the first gas-solid temperature distribution, the gas volume and gas composition in the dripping zone includes:

[0017] Based on the smelting data, obtain the hot metal composition, the ore ratio, the mass fractions of various substances in the mixed ore, the specific surface area of the raw fuel, and the direct reduction degree, and obtain the gas volume and gas composition in the tuyere raceway;

[0018] Calculate the molten iron composition, the ore ratio, the mass fractions of various substances, and the direct reduction degree to obtain the heat consumption for direct reduction in the cohesive zone, and calculate the direct reduction degree, the molten iron composition, the gas volume in the tuyere raceway, and the gas composition in the tuyere raceway to obtain the gas volume and the gas composition in the cohesive zone;

[0019] Calculate the specific surface area of the raw fuel and the first gas-solid temperature distribution in the tuyere raceway to obtain the gas-solid heat exchange amount in the cohesive zone;

[0020] Query and obtain the heat load of the cohesive zone, and calculate based on the preset height of the cohesive zone, combining the heat consumption for direct reduction in the cohesive zone, the gas-solid heat exchange amount in the cohesive zone, and the heat load of the cohesive zone to obtain the first gas-solid temperature distribution in the cohesive zone.

[0021] In an embodiment of the present invention, the obtaining of the first gas-solid temperature distribution, the gas volume, and the gas composition in the cohesive zone includes:

[0022] Calculate and obtain the gas flow distribution in the cohesive zone according to the gas volume and the gas composition in the cohesive zone;

[0023] Based on the smelting data, obtain the raw fuel composition and consumption, the molten iron output, and the slag amount, and calculate the heat consumption for slag-iron melting and the chemical reaction heat in the cohesive zone based on the raw fuel composition and consumption, the molten iron output, the slag amount, and the gas flow distribution in the cohesive zone;

[0024] Based on the smelting data, obtain the specific surface area of the raw fuel, and calculate in combination with the first gas-solid temperature distribution in the cohesive zone to obtain the gas-solid heat exchange amount in the cohesive zone;

[0025] Query and obtain the heat load of the cohesive zone, and calculate based on the preset width of the cohesive zone, combining the heat consumption for slag-iron melting in the cohesive zone, the chemical reaction heat in the cohesive zone, the gas-solid heat exchange amount in the cohesive zone, and the heat load of the cohesive zone to obtain the first gas-solid temperature distribution, the gas volume, and the gas composition in the cohesive zone.

[0026] In an embodiment of the present invention, the obtaining of the first gas-solid temperature distribution, the gas volume, the gas composition, and the first indirect reduction degree in the lump zone includes:

[0027] Based on the smelting data, obtain the raw fuel composition and consumption, and obtain the gas volume and the gas composition in the cohesive zone, and calculate the chemical reaction heat and the first indirect reduction degree in the lump zone based on the raw fuel composition and consumption, the gas volume, and the gas composition in the cohesive zone;

[0028] Based on the smelting data, obtain the specific surface area of the raw fuels, and acquire the first gas-solid temperature distribution in the cohesive zone, and calculate and obtain the gas-solid heat transfer amount in the lump zone according to the specific surface area of the raw fuels and the first gas-solid temperature distribution in the cohesive zone;

[0029] Calculate the chemical reaction heat in the lump zone and the gas-solid heat transfer amount in the lump zone to obtain the first gas-solid temperature distribution in the lump zone, the gas volume in the lump zone, and the gas composition in the lump zone.

[0030] In an embodiment of the present invention, the judgment of the first axial temperature distribution curve of the blast furnace according to the preset convergence condition further includes:

[0031] If the first axial temperature distribution curve of the blast furnace does not meet the preset convergence condition, then use the gas volume and gas composition in the dripping zone as the input data of the calculation model, and combine with the smelting data to calculate and obtain the second gas-solid temperature distribution in the dripping zone, the cohesive zone, and the lump zone, and the second indirect reduction degree of blast furnace ironmaking again;

[0032] Obtain the theoretical combustion temperature of the tuyere raceway zone, and generate a second axial temperature distribution curve of the blast furnace in combination with the second gas-solid temperature distributions in the dripping zone, the cohesive zone, and the lump zone;

[0033] Judge the second axial temperature distribution curve of the blast furnace according to the preset convergence condition;

[0034] If the second axial temperature distribution curve of the blast furnace meets the preset convergence condition, then judge whether the sum of the second indirect reduction degree and the direct reduction degree is the preset threshold;

[0035] If the second axial temperature distribution curve of the blast furnace does not meet the preset convergence condition, then calculate the third gas-solid temperature distributions in the dripping zone, the cohesive zone, and the lump zone again, and generate a third axial temperature distribution curve of the blast furnace until the third axial temperature distribution curve meets the preset convergence condition.

[0036] In an embodiment of the present invention, the judgment of whether the sum of the indirect reduction degree and the direct reduction degree is the preset threshold further includes:

[0037] If the sum of the indirect reduction degree and the direct reduction degree is not the preset threshold, then adjust the height of the dripping zone, and calculate the fourth gas-solid temperature distributions in the dripping zone, the cohesive zone, and the lump zone and the third indirect reduction degree again until the sum of the third indirect reduction degree and the direct reduction degree is the preset threshold.

[0038] An embodiment of the present invention further provides a device for determining the axial temperature distribution of a blast furnace, characterized in that the device includes:

[0039] An acquisition module, configured to acquire smelting data during the blast furnace smelting process, where the smelting data includes the direct reduction degree of blast furnace ironmaking;

[0040] A calculation module, configured to calculate, based on the smelting data, the theoretical combustion temperature of the tuyere raceway in the blast furnace, the first gas-solid temperature distribution in the dripping zone, softening-melting zone, and lump zone in the blast furnace, and the first indirect reduction degree of blast furnace ironmaking;

[0041] A generation module, configured to generate a first blast furnace axial temperature distribution curve according to the theoretical combustion temperature, the first gas-solid temperature distribution in the dripping zone, the first gas-solid temperature distribution in the softening-melting zone, and the first gas-solid temperature distribution in the lump zone;

[0042] A first judgment module, configured to judge the first blast furnace axial temperature distribution curve according to a preset convergence condition;

[0043] A second judgment module, configured to judge whether the sum of the first indirect reduction degree and the direct reduction degree is a preset threshold if the first blast furnace axial temperature distribution curve meets the preset convergence condition;

[0044] A determination module, configured to determine the blast furnace axial temperature distribution according to the first blast furnace axial temperature distribution curve if the sum of the first indirect reduction degree and the direct reduction degree is the preset threshold.

[0045] An embodiment of the present invention further provides a system for determining the blast furnace axial temperature distribution, characterized in that the system includes:

[0046] A data acquisition subsystem, configured to acquire smelting data during the blast furnace smelting process;

[0047] A tuyere raceway calculation subsystem, configured to solve the theoretical combustion temperature of the tuyere raceway, the gas volume in the tuyere raceway, and the gas composition in the tuyere raceway;

[0048] A dripping zone calculation subsystem, configured to solve the gas-solid temperature distribution in the dripping zone, the gas volume in the dripping zone, and the gas composition in the dripping zone;

[0049] A softening-melting zone calculation subsystem, configured to solve the gas-solid temperature distribution in the softening-melting zone, the gas volume in the softening-melting zone, and the gas composition in the softening-melting zone;

[0050] A lump zone calculation subsystem, configured to solve the gas-solid temperature distribution in the lump zone, the gas volume in the lump zone, the gas composition in the lump zone, and the indirect reduction degree of blast furnace ironmaking;

[0051] An iterative solution subsystem is used to determine whether the axial temperature distribution curve of the blast furnace meets a preset convergence condition, and to determine whether the sum of the indirect reduction degree and the direct reduction degree is a preset threshold, and to adjust the calculation data in the iterative solution process until the axial temperature distribution of the blast furnace is determined.

[0052] An embodiment of the present invention also provides an electronic device, including a processor, a memory, and a communication bus;

[0053] The communication bus is used to connect the processor and the memory;

[0054] The processor is used to execute the computer program stored in the memory to implement the method described in any one of the above embodiments.

[0055] An embodiment of the present invention also provides a computer-readable storage medium, on which a computer program is stored,

[0056] The computer program is used to cause the computer to execute the method described in any one of the above embodiments.

[0057] As described above, a method, device, system, equipment, and medium for determining the axial temperature distribution of a blast furnace provided by an embodiment of the present invention have the following beneficial effects:

[0058] First, by obtaining the smelting data during the smelting process of the blast furnace, where the smelting data includes the direct reduction degree of blast furnace ironmaking; secondly, based on the smelting data, calculate the theoretical combustion temperature in the tuyere raceway area in the blast furnace, and the first gas-solid temperature distributions in the dripping zone, softening-melting zone, and lump zone in the blast furnace and the first indirect reduction degree of blast furnace ironmaking; further, generate the first blast furnace axial temperature distribution curve according to the theoretical combustion temperature, the first gas-solid temperature distribution in the dripping zone, the first gas-solid temperature distribution in the softening-melting zone, and the first gas-solid temperature distribution in the lump zone; then, judge the first blast furnace axial temperature distribution curve according to the preset convergence condition. If the first blast furnace axial temperature distribution curve meets the preset convergence condition, then judge whether the sum of the first indirect reduction degree and the direct reduction degree is a preset threshold. If the sum of the first indirect reduction degree and the direct reduction degree is the preset threshold, finally, determine the axial temperature distribution of the blast furnace according to the first blast furnace axial temperature distribution curve. By performing progressive calculations on the smelting data during the smelting process, the axial temperature distribution during blast furnace smelting can be determined in real time and accurately, and the heat state at different height positions in the blast furnace can be obtained, providing a guarantee for the stable operation of ironmaking.

[0059] It should be understood that the above general description and the following detailed description are only exemplary and explanatory, and cannot limit this application. Description of the Drawings

[0060] Figure 1It is a flowchart of a method for determining the axial temperature distribution of a blast furnace shown in an exemplary embodiment of the present application;

[0061] Figure 2 It is a flowchart of a specific method for determining the axial temperature distribution of a blast furnace shown in an exemplary embodiment of the present application;

[0062] Figure 3 It is a schematic diagram of an exemplary axial temperature distribution curve of a blast furnace shown in an exemplary embodiment of the present application;

[0063] Figure 4 It is a schematic diagram of another exemplary axial temperature distribution curve of a blast furnace shown in an exemplary embodiment of the present application;

[0064] Figure 5 It is a block diagram of a device for determining the axial temperature distribution of a blast furnace shown in an exemplary embodiment of the present application;

[0065] Figure 6 It is a block diagram of a system for determining the axial temperature distribution of a blast furnace shown in an exemplary embodiment of the present application;

[0066] Figure 7 It is a schematic structural diagram of an electronic device provided in an embodiment of the present application. Detailed implementation manners

[0067] The following uses specific specific examples to illustrate the implementation manners of the present invention. Those skilled in the art can easily understand other advantages and effects of the present invention from the content disclosed in this specification. The present invention can also be implemented or applied through other different specific implementation manners. Various details in this specification can also be modified or changed based on different viewpoints and applications without departing from the spirit of the present invention. It should be noted that, without conflict, the following embodiments and the features in the embodiments can be combined with each other.

[0068] It should be noted that the diagrams provided in the following embodiments only illustrate the basic concept of the present invention in a schematic manner. Therefore, only the components related to the present invention are shown in the diagrams, rather than being drawn according to the number, shape, and size of the components in actual implementation. The types, quantities, and proportions of the components in actual implementation can be arbitrarily changed, and the component layout type may also be more complex.

[0069] Blast furnace smelting is the main method of modern ironmaking production, which is a production process that completes physical, chemical, and kinetic processes under harsh conditions such as high temperature, high pressure, and airtightness. During the blast furnace smelting process, it is necessary to monitor the axial temperature in the furnace in real time and immediately adjust the furnace condition to be cooler or hotter, so as to ensure the stable and smooth operation of low-carbon blast furnace smelting is crucial. The traditional axial temperature distribution in the blast furnace is determined by real-time monitoring of the basic conditions of the blast furnace temperature. The basic conditions of the blast furnace temperature can reflect the temperature distribution law and its changes on the inner edge furnace wall of the blast furnace, so as to monitor the temperature changes in the furnace. However, due to the high temperature, high pressure, and corrosive substances inside the blast furnace, in production, it is impossible to finely and effectively monitor the axial temperature distribution of the blast furnace in actual production based on the basic conditions of the blast furnace temperature, and it is also impossible to better propose specific methods for effectively maintaining the efficient operation of the blast furnace at the theoretical and technical levels, which affects the stable operation of blast furnace smelting.

[0070] To solve the above problems, the embodiments of the present application provide a method for determining the axial temperature distribution of a blast furnace. Please refer to Figure 1 , Figure 1 which is a flowchart of the method for determining the axial temperature distribution of a blast furnace shown in an exemplary embodiment of the present application. It should be understood that this method can also be applied to other exemplary implementation environments and be specifically executed by devices in other implementation environments. The present embodiment does not limit the implementation environment applicable to this method. To solve these problems, the embodiments of the present application respectively propose a method for determining the axial temperature distribution of a blast furnace, a device for determining the axial temperature distribution of a blast furnace, a system for determining the axial temperature distribution of a blast furnace, an electronic device, and a computer-readable storage medium. These embodiments will be described in detail below.

[0071] As Figure 1 shown, in an exemplary embodiment, the method for determining the axial temperature distribution of a blast furnace at least includes steps S101 to S106, which are introduced in detail as follows:

[0072] Step S101, obtain the smelting data during the blast furnace smelting process, and the smelting data includes the direct reduction degree of blast furnace ironmaking.

[0073] In blast furnace smelting, iron ore, coke, etc. are charged from the top of the furnace, and preheated air is blown into the tuyere raceway zone along the furnace circumference at the lower part of the blast furnace. At high temperatures, carbon monoxide and hydrogen generated by the combustion of carbon in the coke with oxygen in the blown air remove the oxygen in the iron ore during the upward movement in the furnace, thereby reducing iron. The molten iron is discharged from the iron notch, and the unreduced impurities in the iron ore, etc. form slag, which is discharged from the slag notch, and the generated gas is discharged from the top of the furnace. During the smelting process, there are certain smelting data at each stage. Therefore, by obtaining the smelting data during the blast furnace smelting process, the subsequent determination of the axial temperature distribution of the blast furnace can be carried out. Among them, the smelting data includes the direct reduction degree of blast furnace ironmaking, which refers to the ratio of the amount of iron reduced by the direct reduction method from ferrous oxide to the total amount of reduced iron. The direct reduction process only occurs in the dripping zone. It should be understood that the smelting data during the blast furnace ironmaking process is the real-time data of blast furnace smelting, which can directly and real-time reflect the smelting situation in the blast furnace. Therefore, the axial temperature distribution of the blast furnace can be accurately determined according to the smelting data.

[0074] Step S102: Based on the smelting data, calculate the theoretical combustion temperature in the tuyere raceway zone in the blast furnace, and the first gas-solid temperature distribution in the dripping zone, softening-melting zone, and lump zone in the blast furnace, and the first indirect reduction degree of blast furnace ironmaking.

[0075] It should be noted that according to the different forms of existence of materials, the blast furnace can be divided into five regions: lump zone, softening-melting zone, dripping zone, tuyere raceway zone, and slag-iron zone, which are distributed at different height positions inside the blast furnace. In the lump zone, the evaporation and thermal decomposition of moisture in the burden, the reduction of iron ore, and the heat exchange between the burden and the gas are carried out; in the softening-melting zone, the burden softens in this area and begins to melt and drip at the lower boundary, and the direct reduction reaction of blast furnace ironmaking occurs to form the primary slag; in the dripping zone, various chemical reactions between the dripping liquid slag and iron and the gas and solid carbon occur; in the tuyere raceway zone, the injected fuel reacts with the hot air to generate high-temperature gas; in the slag-iron zone, the slag-metal reaction occurs at the interface between the slag and iron layers and when the iron droplets pass through the slag layer. For the confirmation of the axial temperature distribution of the blast furnace, only the gas-solid temperature distribution in the lump zone, softening-melting zone, dripping zone, and tuyere raceway zone is involved, and the slag-iron zone is generally not considered.

[0076] Based on the smelting data during the blast furnace ironmaking process, the theoretical combustion temperature in the tuyere raceway zone in the blast furnace, the first gas-solid temperature distribution in the dripping zone in the blast furnace, the first gas-solid temperature distribution in the softening-melting zone, the first gas-solid temperature distribution in the lump zone, and the first indirect reduction degree of blast furnace ironmaking can be calculated. The calculation process satisfies the mass conservation, energy conservation, and material conservation. Since the tuyere raceway zone, dripping zone, softening-melting zone, and lump zone are distributed at different height positions inside the blast furnace, the axial temperature distribution of the blast furnace can be monitored in real time and accurately according to the theoretical combustion temperature in the tuyere raceway zone and the first gas-solid temperature distribution in the dripping zone, softening-melting zone, and lump zone.

[0077] In one embodiment, based on the smelting data, the gas volume and gas composition in the tuyere raceway, the gas volume and gas composition in the dropping zone, the gas volume and gas composition in the cohesive zone, and the gas volume and gas composition in the lump zone of the blast furnace are also calculated and obtained.

[0078] Based on the smelting data, while calculating the theoretical combustion temperature in the tuyere raceway of the blast furnace, the first gas-solid temperature distribution in the dropping zone of the blast furnace, the first gas-solid temperature distribution in the cohesive zone, the first gas-solid temperature distribution in the lump zone, and the first indirect reduction degree of blast furnace ironmaking, since the calculations of each region are in a progressive relationship, it is also necessary to calculate the gas volume and gas composition in the tuyere raceway, the gas volume and gas composition in the dropping zone, the gas volume and gas composition in the cohesive zone, and the gas volume and gas composition in the lump zone, so as to complete the calculations of each region.

[0079] In one embodiment, the obtaining of the theoretical combustion temperature, the gas volume and gas composition in the tuyere raceway includes:

[0080] Based on the smelting data, the raw material and fuel composition and consumption, the blast volume, the blast temperature, the blast oxygen enrichment rate, and the gas injection volume are obtained;

[0081] The raw material and fuel composition and consumption, the blast volume, the blast temperature, the blast oxygen enrichment rate, and the gas injection volume are calculated to obtain the theoretical combustion temperature, the gas volume and gas composition in the tuyere raceway.

[0082] In the smelting data, it also includes the raw material and fuel composition and consumption, the blast volume, the blast temperature, the blast oxygen enrichment rate, the gas injection volume, the total heat released during the combustion processes of coke, pulverized coal, and injected gas in the tuyere raceway, and the specific heat capacity of the tuyere raceway gas. Based on the smelting data, the raw material and fuel composition and consumption, the blast volume, the blast temperature, the blast oxygen enrichment rate, the gas injection volume, the total heat released during the combustion processes of coke, pulverized coal, and injected gas in the tuyere raceway, and the specific heat capacity of the tuyere raceway gas of blast furnace ironmaking can be obtained. According to the local material balance and heat balance in the tuyere raceway, the raw material and fuel composition and consumption, the blast volume, the blast temperature, the blast oxygen enrichment rate, the gas injection volume, the total heat released during the combustion processes of coke, pulverized coal, and injected gas in the tuyere raceway, and the specific heat capacity of the tuyere raceway gas are calculated, and the theoretical combustion temperature, the gas volume and gas composition in the tuyere raceway can be obtained.

[0083] It should be noted that the gas volume in the tuyere raceway is calculated by converting all pulverized coal and injected gas into CO, H2, and N2, all the oxygen enrichment (O2) in the blast into CO, and all the H2O in the pulverized coal and blast into CO and H2. The above conversions are affected by the blast temperature. The calculation formula for the gas volume in the tuyere raceway is as follows: In the formula, V BG represents the gas volume in the tuyere raceway, are the volumes of CO generated by the reactions of coke carbon, pulverized coal, injected gas, and water in the blast in the tuyere raceway respectively. Among them, contains the CO generated by the reactions of the fixed carbon, volatile carbon, and moisture in the pulverized coal; are the volumes of H2 generated by the reactions of pulverized coal, injected gas, and water in the blast in the tuyere raceway respectively. Among them, contains the H2 generated by the reactions of the volatile hydrogen and moisture in the pulverized coal; are the volumes of N2 in pulverized coal, injected gas, and blast respectively. Among them, contains the N2 in the volatile matter of pulverized coal and the carrier gas of pulverized coal. The calculation formula for the theoretical combustion temperature in the tuyere raceway is as follows: In the formula, Q b is the total heat released during the combustion of coke, pulverized coal, and injected gas in the tuyere raceway. V BG represents the gas volume in the tuyere raceway. C BG is the specific heat capacity of the gas in the tuyere raceway; the total heat released during the combustion process Q b The calculation formula is: Q b =H C-CO ·(M Coke ·ω Coke-C +M Coal ·ω Coal-C )+H CH4 ·(M Coke ·ω Coke-CH4 +M Coal ·ω Coal-CH4 +M Gas ·ω Gas-CH4 ) In the formula, M coal and M gas represent the pulverized coal consumption per ton of iron and the gas injection amount per ton of iron respectively. H C-CO is the reaction enthalpy of the reaction 2C + O2 = 2CO. H CH4 is the reaction enthalpy of the reaction 2CH4 + O2 = 2CO + 4H2. ω Coke-C , ω Coal-C , ω Coke-CH4 , ω Coal-CH4 , ω Gas-CH4They are respectively the mass fraction of solid carbon in coke, the mass fraction of solid carbon in pulverized coal, the mass fraction of CH4 in coke, the mass fraction of CH4 in pulverized coal, the mass fraction of CH4 in the blown gas, and M coke is the amount of coke participating in combustion in the tuyere raceway zone.

[0084] In one embodiment, obtaining the first gas-solid temperature distribution, the gas volume, and the gas composition in the dripping zone includes:

[0085] Based on the smelting data, obtain the molten iron composition, ore ratio, mass fractions of various substances in the mixed ore, specific surface area of the raw fuel, and direct reduction degree, and obtain the gas volume and gas composition in the tuyere raceway zone;

[0086] Calculate the molten iron composition, ore ratio, mass fractions of various substances, and direct reduction degree to obtain the heat consumption for direct reduction in the dripping zone, and calculate the direct reduction degree, molten iron composition, gas volume, and gas composition in the tuyere raceway zone to obtain the gas volume and gas composition in the dripping zone;

[0087] Calculate the specific surface area of the raw fuel and the theoretical combustion temperature to obtain the gas-solid heat exchange amount in the dripping zone;

[0088] Query and obtain the heat load in the dripping zone, and calculate based on the preset height of the dripping zone, combined with the heat consumption for direct reduction in the dripping zone, the gas-solid heat exchange amount in the dripping zone, and the heat load in the dripping zone to obtain the first gas-solid temperature distribution in the dripping zone.

[0089] The smelting data also includes the molten iron composition, specific surface area of the raw fuel, ore ratio, and mass fractions of various substances in the mixed ore. First, based on the smelting data, obtain the molten iron composition, specific surface area of the raw fuel, ore ratio, mass fractions of various substances in the mixed ore, and direct reduction degree. At the same time, also obtain the gas volume and gas composition in the tuyere raceway zone. Then calculate the molten iron composition, ore ratio, mass fractions of various substances in the mixed ore, and direct reduction degree to obtain the heat consumption for direct reduction in the dripping zone, which includes the heat consumption for direct reduction of iron oxidation and the heat consumption for direct reduction of other oxides such as silicon, manganese, sulfur, phosphorus, and titanium; further, calculate the direct reduction degree, molten iron composition, gas volume, and gas composition in the tuyere raceway zone to obtain the gas volume and gas composition in the dripping zone; then calculate the specific surface area of the raw fuel and the theoretical combustion temperature to obtain the gas-solid heat exchange amount in the dripping zone; finally, query and obtain the heat load in the dripping zone, and calculate based on the preset height of the dripping zone, combined with the heat consumption for direct reduction in the dripping zone, the gas-solid heat exchange amount in the dripping zone, and the heat load in the dripping zone to obtain the first gas-solid temperature distribution in the dripping zone.

[0090] The heat consumption for direct reduction Q in the dripping zone Rd The calculation formula is: where M 矿 is the ore ratio, ω Fe and ω FeS are the mass fractions of Fe and FeS in the mixed ore respectively, Rd is the direct reduction degree, [Si], [Mn], [P], [Ti], [S] are the contents of each element in the hot metal, i.e., the hot metal composition, and H is the reaction heat effect involved in the direct reduction of each element, as shown in Table 1.

[0091] Table 1: Reaction Heat Effects Involved in the Direct Reduction of Each Element

[0092] Number Reaction formula Thermal effect 1 FeO + C = Fe + CO <![CDATA[H Fe > 2 <![CDATA[SiO2 + 2C = Si + 2CO]]> <![CDATA[H si > 3 MnO + C = Mn + CO <![CDATA[H Mn > 4 <![CDATA[P2O5 + 5C = 2P + 5CO]]> <![CDATA[H P > 5 <![CDATA[TiO2 + 2C = Ti + 2CO]]> <![CDATA[H Ti > 6 FeS + CaO + C = Fe + CaS + CO <![CDATA[H S >

[0093] It should be noted that the blast furnace body is divided into five parts: the furnace throat, the furnace shaft, the furnace waist, the furnace belly, and the furnace hearth. The furnace throat plays a role in burden distribution, the furnace shaft, the furnace waist, and the furnace belly play a role in carrying the burden and heat preservation and preheating of the burden, and the furnace hearth plays a role in discharging slag and hot metal. The heat load varies at different positions of the blast furnace. Among them, the heat load of the furnace hearth accounts for 8% - 12% of the total heat load, the heat load of the furnace belly accounts for 25% - 30% of the total heat load, the heat load of the furnace waist accounts for 20% - 25% of the total heat load, and the heat load of the furnace shaft accounts for 35% - 45% of the total heat load. The dripping zone, softening - melting zone, and lump zone of the blast furnace are specifically determined according to the situation to correspond to the furnace shaft, the furnace waist, the furnace belly, or the furnace hearth respectively, so as to determine the heat load range of each region, and then according to actual production, the heat loads of the dripping zone, softening - melting zone, and lump zone are pre - determined.

[0094] In one embodiment, obtaining the first gas - solid temperature distribution, gas volume, and gas composition in the softening - melting zone includes:

[0095] Calculating and obtaining the gas flow distribution in the dripping zone based on the gas volume and gas composition in the dripping zone;

[0096] Based on the smelting data, obtaining the raw fuel composition and consumption, hot metal output, and slag amount, and calculating the heat consumption for slag - iron melting and chemical reaction heat in the softening - melting zone based on the raw fuel composition and consumption, hot metal output, slag amount, and gas flow distribution in the dripping zone;

[0097] Based on the smelting data, obtaining the specific surface area of the raw fuel, and calculating the gas - solid heat exchange amount in the softening - melting zone in combination with the first gas - solid temperature distribution in the dripping zone;

[0098] Querying and obtaining the heat load in the softening - melting zone, and calculating based on the preset width of the softening - melting zone, in combination with the heat consumption for slag - iron melting, chemical reaction heat, gas - solid heat exchange amount, and heat load in the softening - melting zone to obtain the first gas - solid temperature distribution, gas volume, and gas composition in the softening - melting zone.

[0099] First, based on the gas volume and gas composition in the dripping zone, calculate and obtain the gas flow distribution in the dripping zone. Secondly, the smelting data also includes the molten iron output and slag volume. Therefore, based on the smelting data, obtain the raw fuel composition and consumption, molten iron output, and slag volume, and calculate the raw fuel composition and consumption, molten iron output, slag volume, and gas flow distribution in the dripping zone to obtain the heat consumption for slag-iron melting in the softening-melting zone and the chemical reaction heat in the softening-melting zone. Among them, the chemical reaction heat in the softening-melting zone includes the indirect reduction heat of CO in the softening-melting zone, the indirect reduction heat of H2, the evaporation heat of free water in the ore, the decomposition heat of the flux, etc. Further, based on the smelting data, obtain the specific surface area of the raw fuel, and calculate in combination with the first gas-solid temperature distribution in the dripping zone to obtain the gas-solid heat exchange amount in the softening-melting zone. Finally, query and obtain the heat load of the softening-melting zone, and calculate based on the preset width of the softening-melting zone, in combination with the heat consumption for slag-iron melting in the softening-melting zone, the chemical reaction heat in the softening-melting zone, the gas-solid heat exchange amount in the softening-melting zone, and the heat load of the softening-melting zone to obtain the first gas-solid temperature distribution in the softening-melting zone, the gas volume in the softening-melting zone, and the gas composition in the softening-melting zone.

[0100] In one embodiment, obtaining the first gas-solid temperature distribution in the lump zone, the gas volume in the lump zone, the gas composition in the lump zone, and the first indirect reduction degree includes:

[0101] Based on the smelting data, obtain the raw fuel composition and consumption, and obtain the gas volume in the softening-melting zone and the gas composition in the softening-melting zone, and calculate the raw fuel composition and consumption, the gas volume in the softening-melting zone, and the gas composition in the softening-melting zone to obtain the chemical reaction heat in the lump zone and the first indirect reduction degree;

[0102] Based on the smelting data, obtain the specific surface area of the raw fuel, and obtain the first gas-solid temperature distribution in the softening-melting zone, and calculate and obtain the gas-solid heat exchange amount in the lump zone according to the specific surface area of the raw fuel and the first gas-solid temperature distribution in the softening-melting zone;

[0103] Calculate the chemical reaction heat in the lump zone and the gas-solid heat exchange amount in the lump zone to obtain the first gas-solid temperature distribution in the lump zone, the gas volume in the lump zone, and the gas composition in the lump zone.

[0104] The smelting data includes the composition and consumption of the raw fuels. Based on the smelting data, the composition and consumption of the raw fuels are obtained. At the same time, the gas volume and gas composition in the cohesive zone are acquired, and calculations are performed on the composition and consumption of the raw fuels, the gas volume and gas composition in the cohesive zone to obtain the chemical reaction heat in the lumpy zone and the first indirect reduction degree. Among them, the chemical reaction heat in the lumpy zone includes the indirect reduction heat of CO in the lumpy zone, the indirect reduction heat of H2, the evaporation heat of free water in the ore, the decomposition heat of the flux, etc. The first indirect reduction degree is the reduction degree when the iron ore moves down from the burden surface to the root of the cohesive zone after charging, and is obtained from the reduction process of iron oxides in the lumpy zone and the cohesive zone; the smelting data includes the specific surface area of the raw fuels. Based on the smelting data, the specific surface area of the raw fuels is obtained. At the same time, the first gas-solid temperature distribution in the cohesive zone is acquired, and the gas-solid heat transfer amount in the lumpy zone is calculated and obtained according to the specific surface area of the raw fuels and the first gas-solid temperature distribution in the cohesive zone; finally, calculations are performed on the chemical reaction heat in the lumpy zone and the gas-solid heat transfer amount in the lumpy zone to obtain the first gas-solid temperature distribution in the lumpy zone, the gas volume in the lumpy zone and the gas composition in the lumpy zone.

[0105] It should be noted that the gas-solid heat transfer amount is determined by the gas-solid heat transfer coefficient, the specific surface area of the raw fuels, the gas temperature and the burden temperature, that is, Q = h g-s A(T g -T s ), where Q represents the gas-solid heat transfer amount, h g-s represents the gas-solid heat transfer coefficient, A represents the specific surface area of the raw fuels, T g represents the gas temperature, T s represents the burden temperature. Among them, the gas temperature and the burden temperature are the gas temperature and the burden temperature in the above-mentioned theoretical combustion temperature, the first gas-solid temperature distribution in the dripping zone or the first gas-solid temperature distribution in the cohesive zone. Among them, the gas-solid heat transfer coefficient h g-s is determined by the correction coefficient γ, the gas thermal conductivity k g , the burden particle size d s and the Nusselt number Nu, that is Secondly, the calculation formula for the Nusselt number Nu is Nu = 2.0 + 0.6(9Re p ) 1 / 2 Pr 1 / 3 , where Re p is the Reynolds number, where ρ g is the gas density, d s is the burden particle size, u g is the gas flow velocity, μ g is the gas viscosity of the mixed gas, Pr is the Prandtl number, c g is the gas heat capacity, k g is the gas thermal conductivity, μ g is the gas viscosity of the mixed gas.

[0106] Step S103: Generate the first axial temperature distribution curve of the blast furnace based on the theoretical combustion temperature, the first gas-solid temperature distribution in the dripping zone, the first gas-solid temperature distribution in the softening-melting zone, and the first gas-solid temperature distribution in the lump zone.

[0107] After determining the theoretical combustion temperature in the tuyere raceway, the first gas-solid temperature distribution in the dripping zone, the first gas-solid temperature distribution in the softening-melting zone, and the first gas-solid temperature in the lump zone, the first axial temperature distribution curve of the blast furnace can be generated based on the theoretical combustion temperature, the first gas-solid temperature distribution in the dripping zone, the first gas-solid temperature distribution in the softening-melting zone, and the first gas-solid temperature distribution in the lump zone. The axial temperature distribution curve of the blast furnace includes the gas temperature distribution curve and the burden temperature distribution curve. It should be understood that the tuyere raceway, the dripping zone, the softening-melting zone, and the lump zone are at different height positions in the blast furnace, and the theoretical combustion temperature in the tuyere raceway, the first gas-solid temperature distribution in the dripping zone, the first gas-solid temperature distribution in the softening-melting zone, and the first gas-solid temperature distribution in the lump zone can reflect the gas-solid temperature distribution at different height positions in the blast furnace.

[0108] Step S104: Judge the first axial temperature distribution curve of the blast furnace according to the preset convergence condition.

[0109] It should be noted that the preset convergence condition is that the difference between the axial temperature distribution curves of the blast furnace before and after solving is within 1 °C. Among them, the generated first axial temperature distribution curve of the blast furnace is based on the initialized axial temperature distribution curve of the blast furnace. During the iteration process, the currently generated axial temperature distribution curve of the blast furnace is based on the previously generated axial temperature distribution curve of the blast furnace. For the initialized axial temperature distribution curve of the blast furnace, it is determined based on the initialized burden temperature distribution and gas temperature distribution. Among them, the burden temperature distribution T s = 1273.0 - 975.0x / H, T s is the burden temperature distribution, x represents the height from the tuyere raceway, and H is the total height of the blast furnace (from the tuyere raceway to the furnace top); the gas temperature distribution T g = T f -(T f - 200)x / H, T g is the gas temperature distribution, T f is the theoretical combustion temperature in the tuyere raceway, x represents the height from the tuyere raceway, and H is the total height of the blast furnace.

[0110] In one embodiment, judging the first axial temperature distribution curve of the blast furnace according to the preset convergence condition further includes:

[0111] If the first axial temperature distribution curve of the blast furnace does not meet the preset convergence condition, then the gas volume and gas composition in the dripping zone are used as the input data of the calculation model, and combined with the smelting data, the second gas-solid temperature distribution in the dripping zone, the softening-melting zone, and the lump zone and the second indirect reduction degree of blast furnace ironmaking are calculated and obtained again;

[0112] Obtain the theoretical combustion temperature of the tuyere raceway zone, and combine with the second gas-solid temperature distribution in the dripping zone, softening-melting zone, and lump zone to generate the second axial temperature distribution curve of the blast furnace;

[0113] Judge the second axial temperature distribution curve of the blast furnace according to the preset convergence condition;

[0114] If the second axial temperature distribution curve of the blast furnace meets the preset convergence condition, then judge whether the sum of the second indirect reduction degree and the direct reduction degree is the preset threshold;

[0115] If the second axial temperature distribution curve of the blast furnace does not meet the preset convergence condition, then calculate the third gas-solid temperature distribution in the dripping zone, softening-melting zone, and lump zone again, and generate the third axial temperature distribution curve of the blast furnace until the third axial temperature distribution curve of the blast furnace meets the preset convergence condition.

[0116] When the first axial temperature distribution curve of the blast furnace does not meet the preset convergence condition, then take the gas volume and gas composition in the dripping zone as the input data of the calculation model, combine with the smelting data, calculate again and obtain the second gas-solid temperature distribution in the dripping zone, softening-melting zone, and lump zone and the second indirect reduction degree of blast furnace ironmaking; then obtain the theoretical combustion temperature of the tuyere raceway zone, combine with the second gas-solid temperature distribution in the dripping zone, softening-melting zone, and lump zone to generate the second axial temperature distribution curve of the blast furnace; finally, judge the second axial temperature distribution curve of the blast furnace according to the preset convergence condition. If the second axial temperature distribution curve of the blast furnace meets the preset convergence condition, then judge whether the sum of the second indirect reduction degree and the direct reduction degree is the preset threshold. If the second axial temperature distribution curve of the blast furnace still does not meet the above preset convergence condition, then calculate the third gas-solid temperature distribution in the dripping zone, softening-melting zone, and lump zone again, and generate the third axial temperature distribution curve of the blast furnace until the third axial temperature distribution curve of the blast furnace meets the preset convergence condition, and determine the axial temperature distribution of the blast furnace according to the third axial temperature distribution curve of the blast furnace. The above calculation model is a system for determining the axial temperature distribution of the blast furnace.

[0117] It should be understood that in the iterative process, the currently generated axial temperature distribution curve of the blast furnace is based on the previously generated axial temperature distribution curve of the blast furnace to judge whether it meets the preset convergence condition. For example, the above judgment of the second axial temperature distribution curve of the blast furnace according to the preset convergence condition is to judge whether the difference between the second axial temperature distribution curve of the blast furnace and the first axial temperature distribution curve of the blast furnace is within 1 °C.

[0118] Step S105, if the first axial temperature distribution curve of the blast furnace meets the preset convergence condition, then judge whether the sum of the first indirect reduction degree and the direct reduction degree is the preset threshold.

[0119] When the axial temperature distribution curve of the first blast furnace satisfies the preset convergence condition, it is determined whether the sum of the first indirect reduction degree and the direct reduction degree is a preset threshold, where the preset threshold is 1, that is, it is determined whether the sum of the first indirect reduction degree and the direct reduction degree is 1.

[0120] In one embodiment, determining whether the sum of the indirect reduction degree and the direct reduction degree is a preset threshold further includes:

[0121] If it is determined that the sum of the indirect reduction degree and the direct reduction degree is not the preset threshold, the height of the cohesive zone is adjusted, and the fourth gas-solid temperature distribution and the third indirect reduction degree of the cohesive zone, softening-melting zone, and burden zone are calculated again until the sum of the third indirect reduction degree and the direct reduction degree is the preset threshold.

[0122] When the sum of the indirect reduction degree and the direct reduction degree is not the preset threshold 1, the height of the cohesive zone is adjusted, and the fourth gas-solid temperature distribution and the third indirect reduction degree of the cohesive zone, softening-melting zone, and burden zone are calculated again until the sum of the third indirect reduction degree and the direct reduction degree is the preset threshold 1, and the iteration stops. It should be noted that the above-mentioned third indirect reduction degree is the corresponding indirect reduction degree obtained when calculating the fourth gas-solid temperature distribution of the cohesive zone, softening-melting zone, and burden zone.

[0123] Step S106, if the sum of the first indirect reduction degree and the direct reduction degree is the preset threshold, the axial temperature distribution of the blast furnace is determined according to the axial temperature distribution curve of the first blast furnace.

[0124] When the sum of the first indirect reduction degree and the direct reduction degree is the preset threshold 1, it indicates that the obtained axial temperature distribution curve of the first blast furnace is relatively accurate and can reflect the temperature distribution of the blast furnace at different positions in real time and accurately. Then, the axial temperature distribution inside the blast furnace is determined according to this axial temperature distribution curve of the first blast furnace.

[0125] Please refer to Figure 2 , Figure 2 which is a flowchart of a specific method for determining the axial temperature distribution of a blast furnace shown in an exemplary embodiment of the present application. As Figure 2 shown, a specific method for determining the axial temperature distribution of a blast furnace at least includes steps S201 to S206, which are introduced in detail as follows:

[0126] Step S201, obtain the raw fuel composition and consumption, blast volume, blast temperature, blast oxygen enrichment rate, gas injection volume, total heat released during the combustion process of coke, pulverized coal, and injected gas in the tuyere raceway, specific heat capacity of the gas in the tuyere raceway, ore ratio, mass fraction of each substance in the mixed ore, direct reduction degree, hot metal composition, hot metal production, slag volume, and specific surface area of the raw fuel.

[0127] Step S202: Calculate the original fuel composition and dosage, blast volume, blast temperature, blast oxygen enrichment rate, coal injection volume, total heat released during the combustion process of coke, pulverized coal, and injected gas in the tuyere raceway, and the specific heat capacity of the gas in the tuyere raceway to obtain the theoretical combustion temperature, gas volume in the tuyere raceway, and gas composition in the tuyere raceway.

[0128] Step S203: Calculate the hot metal composition, ore ratio, mass fraction of each substance in the mixed ore, and direct reduction degree to obtain the heat consumption for direct reduction in the dripping zone. Also, calculate the direct reduction degree, hot metal composition, gas volume in the tuyere raceway, and gas composition in the tuyere raceway to obtain the gas volume and gas composition in the dripping zone.

[0129] Step S204: Calculate the specific surface area of the original fuel and the theoretical combustion temperature to obtain the gas-solid heat exchange amount in the dripping zone.

[0130] Step S205: Query and obtain the heat load of the dripping zone, and based on the preset height of the dripping zone, calculate by combining the heat consumption for direct reduction in the dripping zone, the gas-solid heat exchange amount in the dripping zone, and the heat load of the dripping zone to obtain the gas-solid temperature distribution in the dripping zone.

[0131] Step S206: Calculate and obtain the gas flow distribution in the dripping zone based on the gas volume and gas composition in the dripping zone. Also, calculate the original fuel composition and dosage, hot metal production, slag volume, and gas flow distribution in the dripping zone to obtain the heat consumption for slag-iron melting and the chemical reaction heat in the softening-melting zone.

[0132] Step S207: Calculate the specific surface area of the original fuel and the gas-solid temperature distribution in the dripping zone to obtain the gas-solid heat exchange amount in the softening-melting zone.

[0133] Step S208: Query and obtain the heat load of the softening-melting zone, and based on the preset width of the softening-melting zone, calculate by combining the heat consumption for slag-iron melting in the softening-melting zone, the chemical reaction heat in the softening-melting zone, the gas-solid heat exchange amount in the softening-melting zone, and the heat load of the softening-melting zone to obtain the gas-solid temperature distribution, gas volume, and gas composition in the softening-melting zone.

[0134] Step S209: Calculate the original fuel composition and dosage, gas volume, and gas composition in the softening-melting zone to obtain the chemical reaction heat and indirect reduction degree in the lump zone.

[0135] Step S210: Calculate the specific surface area of the original fuel and the gas-solid temperature distribution in the softening-melting zone to obtain the gas-solid heat exchange amount in the lump zone.

[0136] Step S211: Calculate the chemical reaction heat and gas-solid heat exchange amount in the lump zone to obtain the gas-solid temperature distribution, gas volume, and gas composition in the lump zone.

[0137] Step S212: Generate the axial temperature distribution curve of the blast furnace based on the theoretical combustion temperature, the gas-solid temperature distribution in the dripping zone, the gas-solid temperature distribution in the softening-melting zone, and the gas-solid temperature distribution in the lump zone.

[0138] Step S213: Determine whether the axial temperature distribution curve of the blast furnace meets the preset convergence condition.

[0139] If it meets the condition, go to Step S214-1; if it does not meet the condition, return to Step S203, that is, take the gas volume and gas composition in the dripping zone as the input data of the calculation model, and combine with each smelting data to calculate again and obtain the gas-solid temperature distribution in the dripping zone, softening-melting zone, lump zone, and the indirect reduction degree of blast furnace ironmaking.

[0140] Step S214-1: Determine whether the sum of the indirect reduction degree and the direct reduction degree is a preset threshold.

[0141] If yes, go to Step S215; if not, go to Step S214-2.

[0142] Step S214-2: Adjust the height of the dripping zone, and then return to Step S205.

[0143] Step S215: Determine the axial temperature distribution of the blast furnace according to the current axial temperature distribution curve of the blast furnace.

[0144] The following takes a 2300m 3 blast furnace as an example, where the effective volume utilization coefficient of the blast furnace is 3.7t / (d·m3).

[0145] Implementation Case 1 - Conventional blast furnace temperature distribution without injecting hydrogen-rich medium.

[0146] According to Table 2 and Table 3 below, using the local material balance and heat balance in the tuyere raceway, the theoretical combustion temperature is calculated to be 2308°C, the gas volume in the tuyere raceway is 1306m3 / tFe, and the gas composition is 40.62% CO, 6.19% H2, and 53.19% N2. For another example, according to the direct reduction degree and hot metal composition shown in Table 2, the total heat consumption for direct reduction in the dripping zone is calculated to be 1171.29 MJ / tFe; according to the slag amount and other data shown in Table 2 above, the total heat consumption for slag-iron melting in the softening-melting zone is calculated to be 1966.21 MJ / tFe; other data are not elaborated here one by one. Due to the large difference in heat load at different heights in the blast furnace, in the implementation case, the heat loads of the hearth, bosh, waist, lower part of the shaft, middle part of the shaft, and upper part of the shaft are selected to account for 10%, 28%, 22%, 15%, 10%, and 15% of the total heat load respectively, and the width of the softening-melting zone is kept unchanged at 1.8m.

[0147] Table 2: Blast furnace operation parameters and hot metal composition

[0148] Item Value Item Value Ore ratio, kg / tFe 1568.79 Direct reduction degree of iron 0.45 Dry coke ratio, kg / tFe 350.00 Slag volume, kg / tFe 308.00 Dry coal ratio, kg / tFe 148.00 Fe content in hot metal 94.97％ Hot blast temperature, °C 1180.00 C content in hot metal 4.20％ Oxygen enrichment rate of blast 4.49％ Si content in hot metal 0.29％ <![CDATA[Blast volume, m 3 tFe]]> 905.39 Mn content in hot metal 0.07％ <![CDATA[Oxygen volume, m 3 / tFe]]> 54.57 P content in hot metal 0.15％ Total heat load, GJ / h 90.00 S content in hot metal 0.11％ <![CDATA[Coal injection rate, m 3 / tFe]]> 0.00 Ti content in hot metal 0.21％

[0149] Table 3: Composition and dosage of raw fuels

[0150] Coke Pulverized coal Mixed ore Flux Moisture content, % 0.1 1 0.32 2 TFe, % 0.064 0.215 59.058 FeO 8.492 Fe2O3 0.091 0.307 74.049 CaO 0.092 0.096 8.072 47.940 SiO2 6.749 6.623 5.689 1.690 MgO 0.138 0.096 1.642 5.480 Al2O3 4.630 2.217 1.783 1.140 MnO 0.010 FeS 0.008 P2O5 0.150 0.050 TiO2 0.085 CO2 43.700 H2O C_solid 86.500 73.066 S 0.740 0.548 CO 0.225 CO2 0.240 CH4 0.101 9.214 H2 0.267 3.609 N2 0.227 1.689 O2 2.534

[0151] The gas-solid heat exchange amount and chemical reaction heat at different heights of the blast furnace (lump zone, dripping zone, softening-melting zone) are controlled by the mass conservation, energy conservation, and material conservation governing equations, and are generally expressed as In the formula, ρ represents the fluid density, t represents time, u represents the velocity of the fluid in different directions, φ represents the general variable, Γ represents various coefficients, S represents the generalized source term, and φ, Γ, and S have specific meanings in the governing equations corresponding to different conservation laws, as shown in Table 4.

[0152] Table 4: Specific meanings of φ, Γ, and S in the governing equations corresponding to different conservation laws

[0153]

[0154] In the table, p is the pressure, T is the temperature, k is the heat transfer coefficient of the gas, c p is the specific heat capacity, S T is the internal heat source, C s is the volume concentration of component s, D s is the diffusion coefficient of component s, S s is the amount of component s generated by the reaction per unit volume per unit time in the calculation system.

[0155] For the obtained axial temperature distribution curve of the blast furnace, please refer to Figure 3 , Figure 3 which is a schematic diagram of an exemplary axial temperature distribution curve of the blast furnace shown in an exemplary embodiment of the present application.

[0156] Implementation Case 2 - Axial temperature distribution of a low-carbon blast furnace during hydrogen-rich medium injection

[0157] In this implementation case, the temperature distribution inside the low-carbon blast furnace considering the injection of hydrogen-rich medium is considered. Among them, the composition of the injected hydrogen-rich medium is shown in Table 5.

[0158] Table 5: Composition of the injection medium

[0159] Composition CO <![CDATA[CO2]]> <![CDATA[CH4]]> <![CDATA[H2]]> <![CDATA[N2]]> Proportion 7.8％ 2.4％ 25％ 63.3％ 1.5％

[0160] The injection amount of the hydrogen-rich medium is 100 m3 / tFe. The composition and dosage of its raw fuels refer to Table 3. The proportion of the heat load at different heights of the blast furnace to the total heat load is exactly the same as that in Implementation Case 1 and will not be repeated. Only the differences are described in Implementation Case 2. After injecting the hydrogen-rich medium, in order to meet the material balance and heat balance inside the blast furnace, the blast parameters and the dosage of raw fuels are adjusted, as shown in Table 6.

[0161] Table 6: Parameter Adjustment

[0162] Item Value Item Value Ore ratio, kg / tFe 1573.71 Slag volume, kg / tFe 348.53 Dry coke ratio, kg / tFe 339.62 Hot blast temperature, °C 1180.00 Dry coal ratio, kg / tFe 112.48 Oxygen enrichment rate of blast 9.00％ Direct reduction degree of iron 0.32 Blast volume, m3 / tFe 773.37 Total heat load, GJ / h 90.00 Oxygen volume, m3 / tFe 78.43

[0163] Based on the data in Table 3 and Table 6 above, using the local material balance and heat balance in the tuyere raceway, the theoretical combustion temperature is calculated to be 2187°C, the gas volume in the tuyere raceway is 1301 m3 / tFe, and the gas composition in the tuyere raceway is 41.35% CO, 13.01% H2, and 45.64% N2. Comparing with Example 1, it is found that the H2 content in the gas increases significantly. For another example, based on the direct reduction degree shown in Table 6 above and the hot metal composition shown in Table 2, the total heat consumption for direct reduction in the dripping zone is calculated to be 879.48 MJ / tFe; based on the slag amount and other data shown in Table 6 above, the total heat consumption for melting slag and iron in the softening-melting zone is calculated to be 1951.10 MJ / tFe; other data will not be elaborated here one by one

[0164] For the obtained axial temperature distribution curve of the blast furnace, please refer to Figure 4 , Figure 4 which is a schematic diagram of another exemplary axial temperature distribution curve of the blast furnace shown in an exemplary embodiment of the present application

[0165] The method for determining the axial temperature distribution of the blast furnace provided in the above embodiment first obtains the smelting data during the smelting process of the blast furnace, where the smelting data includes the direct reduction degree of blast furnace ironmaking; secondly, based on the smelting data, the theoretical combustion temperature in the tuyere raceway in the blast furnace, and the first gas-solid temperature distributions in the dripping zone, softening-melting zone, and lump zone in the blast furnace and the first indirect reduction degree of blast furnace ironmaking are calculated respectively; further, a first axial temperature distribution curve of the blast furnace is generated according to the theoretical combustion temperature, the first gas-solid temperature distribution in the dripping zone, the first gas-solid temperature distribution in the softening-melting zone, and the first gas-solid temperature distribution in the lump zone; then, the first axial temperature distribution curve is judged according to the preset convergence condition. If the first axial temperature distribution curve meets the preset convergence condition, it is further judged whether the sum of the first indirect reduction degree and the direct reduction degree is a preset threshold. If the sum of the first indirect reduction degree and the direct reduction degree is the preset threshold, finally, the axial temperature distribution of the blast furnace is determined according to the first axial temperature distribution curve. By performing progressive calculations on the smelting data during the smelting process, the axial temperature distribution during blast furnace smelting is determined in real time and accurately, and the heat state at different height positions in the blast furnace can be obtained, providing a guarantee for the stable operation of ironmaking

[0166] Please refer to Figure 5 , Figure 5 which is a block diagram of the device for determining the axial temperature distribution of the blast furnace shown in an exemplary embodiment of the present application, as Figure 5As shown in the figure, this embodiment provides a device 500 for determining the axial temperature distribution of a blast furnace. The device includes:

[0167] An acquisition module 501, configured to acquire smelting data during the smelting process of the blast furnace. The smelting data includes the direct reduction degree of blast furnace ironmaking.

[0168] A calculation module 502, configured to calculate, based on the smelting data, the theoretical combustion temperature of the tuyere raceway in the blast furnace, the first gas-solid temperature distribution in the dripping zone, softening-melting zone, and lump zone in the blast furnace, and the first indirect reduction degree of blast furnace ironmaking.

[0169] A generation module 503, configured to generate a first blast furnace axial temperature distribution curve according to the theoretical combustion temperature, the first gas-solid temperature distribution in the dripping zone, the first gas-solid temperature distribution in the softening-melting zone, and the first gas-solid temperature distribution in the lump zone.

[0170] A first judgment module 504, configured to judge the first blast furnace axial temperature distribution curve according to a preset convergence condition.

[0171] A second judgment module 505, configured to judge whether the sum of the first indirect reduction degree and the direct reduction degree is a preset threshold if the first blast furnace axial temperature distribution curve meets the preset convergence condition.

[0172] A determination module 506, configured to determine the axial temperature distribution of the blast furnace according to the first blast furnace axial temperature distribution curve if the sum of the first indirect reduction degree and the direct reduction degree is the preset threshold.

[0173] In this embodiment, the device essentially sets multiple modules to execute the method in any of the above embodiments. For the specific functions and technical effects, reference may be made to the above embodiments, and details are not described herein again.

[0174] Please refer to Figure 6 , Figure 6 which is a schematic diagram of a system for determining the axial temperature distribution of a blast furnace shown in an exemplary embodiment of the present application. As Figure 6 shown in the figure, this embodiment provides a system 600 for determining the axial temperature distribution of a blast furnace. The system includes:

[0175] A data acquisition subsystem 601, configured to acquire smelting data during the smelting process of the blast furnace.

[0176] A tuyere raceway calculation subsystem 602, configured to solve the theoretical combustion temperature of the tuyere raceway, the gas volume in the tuyere raceway, and the gas composition in the tuyere raceway.

[0177] A dripping zone calculation subsystem 603, configured to solve the gas-solid temperature distribution in the dripping zone, the gas volume in the dripping zone, and the gas composition in the dripping zone.

[0178] The cohesive zone calculation subsystem 604 is used to solve the gas-solid temperature distribution, the gas volume, and the gas composition in the cohesive zone;

[0179] The lump zone calculation subsystem 605 is used to solve the gas-solid temperature distribution, the gas volume, the gas composition in the lump zone, and the indirect reduction degree of blast furnace ironmaking;

[0180] The iterative solution subsystem 606 is used to determine whether the axial temperature distribution curve of the blast furnace meets the preset convergence condition, and to determine whether the sum of the indirect reduction degree and the direct reduction degree is a preset threshold, and to adjust the calculation data in the iterative solution process until the axial temperature distribution of the blast furnace is determined.

[0181] Among them, the calculations in the tuyere raceway calculation subsystem 602, the dripping zone calculation subsystem 603, the cohesive zone calculation subsystem 604, and the lump zone calculation subsystem 605 are in a progressive relationship, and the data calculated by the previous system are used for the calculation of the next system.

[0182] Please refer to Figure 7 , Figure 7 which is a schematic structural diagram of an electronic device provided by an embodiment of the present application. The electronic device 700 includes a processor 701, a memory 702, and a communication bus 703;

[0183] The communication bus 703 is used to connect the processor 701 and the memory 702;

[0184] The processor 701 is used to execute the computer program stored in the memory 702 to implement the method in one or more of the above embodiments.

[0185] An embodiment of the present invention also provides a computer-readable storage medium, which is characterized in that a computer program is stored thereon,

[0186] and the computer program is used to cause a computer to execute the method in any one of the above Embodiment 1.

[0187] An embodiment of the present application also provides a non-volatile readable storage medium, in which one or more modules (programs) are stored, and when the one or more modules are applied to a device, the device can be caused to execute the instructions (instructions) included in the steps of Embodiment 1 of the embodiments of the present application.

[0188] It should be noted that the above-mentioned computer-readable medium in the present disclosure can be a computer-readable signal medium, a computer-readable storage medium, or any combination of the two. A computer-readable storage medium can be, for example, but not limited to, an electrical, magnetic, optical, electromagnetic, infrared, or semiconductor system, apparatus, or device, or any combination of the above. More specific examples of a computer-readable storage medium can include, but are not limited to: an electrical connection with one or more wires, a portable computer disk, a hard disk, a random access memory (RAM), a read-only memory (ROM), an erasable programmable read-only memory (EPROM or flash memory), an optical fiber, a portable compact disk read-only memory (CD-ROM), an optical storage device, a magnetic storage device, or any suitable combination of the above. In the present disclosure, a computer-readable storage medium can be any tangible medium that contains or stores a program, and this program can be used by or in combination with an instruction execution system, apparatus, or device. In the present disclosure, a computer-readable signal medium can include a data signal propagated in a baseband or as part of a carrier wave, which carries computer-readable program code. Such a propagated data signal can take various forms, including but not limited to electromagnetic signals, optical signals, or any suitable combination of the above. A computer-readable signal medium can also be any computer-readable medium other than a computer-readable storage medium, and this computer-readable signal medium can send, propagate, or transmit a program for use by or in combination with an instruction execution system, apparatus, or device. The program code contained on a computer-readable medium can be transmitted using any appropriate medium, including but not limited to: wires, optical cables, RF (radio frequency), etc., or any suitable combination of the above.

[0189] The above-mentioned computer-readable medium can be included in the above-mentioned electronic device; or it can exist separately without being assembled into the electronic device.

[0190] Computer program code for performing the operations of the present disclosure can be written in one or more programming languages or combinations thereof. The above-mentioned programming languages include object-oriented programming languages - such as Java, Smalltalk, C++, and also include conventional procedural programming languages - such as the "C" language or similar programming languages. The program code can be executed entirely on the user's computer, partially on the user's computer, executed as an independent software package, partially on the user's computer and partially on a remote computer, or entirely on a remote computer or server. In the case of a remote computer, the remote computer can be connected to the user's computer through any type of network - including a local area network (LAN) or a wide area network (WAN) - or can be connected to an external computer (for example, by using an Internet service provider to connect through the Internet).

[0191] The flowcharts and block diagrams in the accompanying drawings illustrate the possible architectures, functions, and operations of methods and computer program products according to various embodiments of the present disclosure. In this regard, each block in the flowchart or block diagram may represent a module, a segment of a program, or a portion of code that contains one or more executable instructions for implementing a specified logical function. It should also be noted that in some alternative implementations, the functions noted in the blocks may occur in a different order than that noted in the accompanying drawings. For example, two consecutive blocks shown may actually be executed substantially in parallel, or they may sometimes be executed in the reverse order, depending on the functions involved. It should also be noted that each block in the block diagram and / or flowchart, and combinations of blocks in the block diagram and / or flowchart, can be implemented by a dedicated hardware-based system that performs the specified functions or operations, or can be implemented by a combination of dedicated hardware and computer instructions.

[0192] The above embodiments are only illustrative of the principles and effects of the present invention, and are not intended to limit the present invention. Any person familiar with this technology can modify or change the above embodiments without departing from the spirit and scope of the present invention. Therefore, all equivalent modifications or changes completed by those with ordinary knowledge in the technical field without departing from the spirit and technical ideas disclosed by the present invention should still be covered by the claims of the present invention.

Claims

1. A method for determining the axial temperature distribution in a blast furnace, characterized in that, The method includes: Obtaining smelting data during the blast furnace smelting process, where the smelting data includes the direct reduction degree of blast furnace ironmaking; Based on the smelting data, calculating respectively the theoretical combustion temperature of the tuyere raceway in the blast furnace, the first gas-solid temperature distribution in the dripping zone, softening-melting zone, and lump zone in the blast furnace, and the first indirect reduction degree of blast furnace ironmaking; Generating a first blast furnace axial temperature distribution curve according to the theoretical combustion temperature, the first gas-solid temperature distribution in the dripping zone, the first gas-solid temperature distribution in the softening-melting zone, and the first gas-solid temperature distribution in the lump zone; Judging the first blast furnace axial temperature distribution curve according to a preset convergence condition, where the preset convergence condition is that the difference between the blast furnace axial temperature distribution curves before and after solution is within 1°C; If the first blast furnace axial temperature distribution curve meets the preset convergence condition, judging whether the sum of the first indirect reduction degree and the direct reduction degree is a preset threshold; If the sum of the first indirect reduction degree and the direct reduction degree is the preset threshold, determining the blast furnace axial temperature distribution according to the first blast furnace axial temperature distribution curve.

2. The method for determining the axial temperature distribution of a blast furnace according to claim 1, characterized in that, Based on the smelting data, also calculating and obtaining the gas volume and gas composition in the tuyere raceway of the tuyere raceway, the gas volume and gas composition in the dripping zone of the dripping zone, the gas volume and gas composition in the softening-melting zone of the softening-melting zone, and the gas volume and gas composition in the lump zone of the lump zone.

3. The method for determining the axial temperature distribution of a blast furnace according to claim 2, wherein The obtaining of the theoretical combustion temperature, the gas volume in the tuyere raceway, and the gas composition in the tuyere raceway includes: Based on the smelting data, obtaining the raw fuel composition and consumption, blast volume, blast temperature, blast oxygen enrichment rate, gas injection volume, the total heat released during the combustion processes of coke, pulverized coal, and injected gas in the tuyere raceway, and the specific heat capacity of the tuyere raceway gas; Calculating the raw fuel composition and consumption, the blast volume, the blast temperature, the blast oxygen enrichment rate, the gas injection volume, the total heat, and the specific heat capacity of the tuyere raceway gas to obtain the theoretical combustion temperature, the gas volume in the tuyere raceway, and the gas composition in the tuyere raceway.

4. The method for determining the axial temperature distribution of a blast furnace according to claim 2, wherein The obtaining of the first gas-solid temperature distribution in the dripping zone, the gas volume in the dripping zone, and the gas composition in the dripping zone includes: Based on the smelting data, obtaining the hot metal composition, ore ratio, mass fractions of various substances in the mixed ore, raw fuel specific surface area, and direct reduction degree, and obtaining the gas volume in the tuyere raceway and the gas composition in the tuyere raceway; Calculating the hot metal composition, the ore ratio, the mass fractions of various substances, and the direct reduction degree to obtain the heat consumption for direct reduction in the dripping zone, and calculating the direct reduction degree, the hot metal composition, the gas volume in the tuyere raceway, and the gas composition in the tuyere raceway to obtain the gas volume in the dripping zone and the gas composition in the dripping zone; Calculating the raw fuel specific surface area and the theoretical combustion temperature to obtain the gas-solid heat exchange amount in the dripping zone; Query and obtain the heat load of the dripping zone, and calculate based on the preset height of the dripping zone, in combination with the heat consumption for direct reduction in the dripping zone, the gas-solid heat transfer amount in the dripping zone, and the heat load of the dripping zone, to obtain the first gas-solid temperature distribution in the dripping zone.

5. The method for determining the axial temperature distribution of a blast furnace according to claim 2, characterized in that, The obtaining of the first gas-solid temperature distribution in the softening-melting zone, the gas volume in the softening-melting zone, and the gas composition in the softening-melting zone includes: Calculate and obtain the gas flow distribution in the dripping zone according to the gas volume in the dripping zone and the gas composition in the dripping zone; Based on the smelting data, obtain the composition and consumption of the burden materials, the hot metal production, and the slag amount, and calculate the heat consumption for slag-iron melting in the softening-melting zone and the chemical reaction heat in the softening-melting zone from the composition and consumption of the burden materials, the hot metal production, the slag amount, and the gas flow distribution in the dripping zone; Based on the smelting data, obtain the specific surface area of the burden materials, and calculate the gas-solid heat transfer amount in the softening-melting zone in combination with the first gas-solid temperature distribution in the dripping zone; Query and obtain the heat load of the softening-melting zone, and calculate based on the preset width of the softening-melting zone, in combination with the heat consumption for slag-iron melting in the softening-melting zone, the chemical reaction heat in the softening-melting zone, the gas-solid heat transfer amount in the softening-melting zone, and the heat load of the softening-melting zone, to obtain the first gas-solid temperature distribution in the softening-melting zone, the gas volume in the softening-melting zone, and the gas composition in the softening-melting zone.

6. The method for determining the axial temperature distribution of a blast furnace according to claim 2, characterized in that, The obtaining of the first gas-solid temperature distribution in the lump zone, the gas volume in the lump zone, the gas composition in the lump zone, and the first degree of indirect reduction includes: Based on the smelting data, obtain the composition and consumption of the burden materials, and obtain the gas volume in the softening-melting zone and the gas composition in the softening-melting zone, and calculate the chemical reaction heat in the lump zone and the first degree of indirect reduction from the composition and consumption of the burden materials, the gas volume in the softening-melting zone, and the gas composition in the softening-melting zone; Based on the smelting data, obtain the specific surface area of the burden materials, and obtain the first gas-solid temperature distribution in the softening-melting zone, and calculate and obtain the gas-solid heat transfer amount in the lump zone according to the specific surface area of the burden materials and the first gas-solid temperature distribution in the softening-melting zone; Calculate the chemical reaction heat in the lump zone and the gas-solid heat transfer amount in the lump zone to obtain the first gas-solid temperature distribution in the lump zone, the gas volume in the lump zone, and the gas composition in the lump zone.

7. The method for determining the axial temperature distribution of a blast furnace according to any one of claims 1 to 6, characterized in that, The judgment of the first axial temperature distribution curve of the blast furnace according to the preset convergence condition further includes: If the first axial temperature distribution curve of the blast furnace does not meet the preset convergence condition, then use the gas volume in the dripping zone and the gas composition in the dripping zone as the input data of the calculation model, and in combination with the smelting data, calculate and obtain the second gas-solid temperature distribution in the dripping zone, the softening-melting zone, and the lump zone, and the second degree of indirect reduction in blast furnace ironmaking again; Obtain the theoretical combustion temperature of the tuyere raceway zone, and generate a second axial temperature distribution curve of the blast furnace in combination with the second gas-solid temperature distribution in the dripping zone, the softening-melting zone, and the lump zone; Judge the second axial temperature distribution curve of the blast furnace according to the preset convergence condition; If the second axial temperature distribution curve of the blast furnace meets the preset convergence condition, then judge whether the sum of the second degree of indirect reduction and the direct reduction degree is the preset threshold; If the axial temperature distribution curve of the second blast furnace does not meet the preset convergence condition, calculate the third gas-solid temperature distribution in the dripping zone, the softening-melting zone, and the lump zone again, and generate a third axial temperature distribution curve of the blast furnace until the third axial temperature distribution curve of the blast furnace meets the preset convergence condition.

8. The method for determining the axial temperature distribution of a blast furnace according to any one of claims 1 to 6, characterized in that, Then, determining whether the sum of the indirect reduction degree and the direct reduction degree is a preset threshold further includes: If it is determined that the sum of the indirect reduction degree and the direct reduction degree is not the preset threshold, adjust the height of the dripping zone, and calculate the fourth gas-solid temperature distribution and the third indirect reduction degree in the dripping zone, the softening-melting zone, and the lump zone again until the sum of the third indirect reduction degree and the direct reduction degree is the preset threshold.

9. A device for determining the axial temperature distribution of a blast furnace, characterized in that, The device includes: An acquisition module, configured to acquire smelting data during the blast furnace smelting process, where the smelting data includes the direct reduction degree of blast furnace ironmaking; A calculation module, configured to respectively calculate the theoretical combustion temperature in the tuyere raceway area in the blast furnace, and the first gas-solid temperature distribution in the dripping zone, the softening-melting zone, and the lump zone in the blast furnace, and the first indirect reduction degree of blast furnace ironmaking based on the smelting data; A generation module, configured to generate a first axial temperature distribution curve of the blast furnace according to the theoretical combustion temperature, the first gas-solid temperature distribution in the dripping zone, the first gas-solid temperature distribution in the softening-melting zone, and the first gas-solid temperature distribution in the lump zone; A first judgment module, configured to judge the first axial temperature distribution curve of the blast furnace according to a preset convergence condition, where the preset convergence condition is that the difference between the axial temperature distribution curves of the blast furnace before and after the solution is within 1°C; A second judgment module, configured to, if the first axial temperature distribution curve of the blast furnace meets the preset convergence condition, judge whether the sum of the first indirect reduction degree and the direct reduction degree is a preset threshold; A determination module, configured to, if the sum of the first indirect reduction degree and the direct reduction degree is the preset threshold, determine the axial temperature distribution of the blast furnace according to the first axial temperature distribution curve of the blast furnace.

10. A system for determining the axial temperature distribution of a blast furnace, characterized in that, The system includes: A data acquisition subsystem, configured to acquire smelting data during the blast furnace smelting process; A tuyere raceway area calculation subsystem, configured to solve the theoretical combustion temperature, the gas volume in the tuyere raceway area, and the gas composition in the tuyere raceway area; A dripping zone calculation subsystem, configured to solve the gas-solid temperature distribution, the gas volume in the dripping zone, and the gas composition in the dripping zone; A softening-melting zone calculation subsystem, configured to solve the gas-solid temperature distribution, the gas volume in the softening-melting zone, and the gas composition in the softening-melting zone; A lump zone calculation subsystem, configured to solve the gas-solid temperature distribution, the gas volume in the lump zone, the gas composition in the lump zone, and the indirect reduction degree of blast furnace ironmaking; An iterative solution subsystem, configured to judge whether the axial temperature distribution curve of the blast furnace meets the preset convergence condition, and judge whether the sum of the indirect reduction degree and the direct reduction degree is a preset threshold, and adjust the calculation data in the iterative solution process until the axial temperature distribution of the blast furnace is determined, where the preset convergence condition is that the difference between the axial temperature distribution curves of the blast furnace before and after the solution is within 1°C.

11. An electronic device, characterized in that, It includes a processor, a memory, and a communication bus; The communication bus is used to connect the processor and the memory; The processor is configured to execute the computer program stored in the memory to implement the method according to any one of claims 1-8.

12. A computer-readable storage medium, characterized in that, which stores a computer program, the computer program is configured to cause the computer to execute the method according to any one of claims 1-8.

Images