{'fraction_non_alphanumeric': 0.034102597068655184, 'fraction_numerical': 0.02256364103882746, 'mean_word_length': 4.848308270676692, 'pattern_counts': {'":': 0, '<': 0, '<?xml version=': 0, '>': 0, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 0, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Use of Hydrogen–Rich Gas in Blast Furnace Ironmaking of V–bearing Titanomagnetite: Mass and Energy Balance Calculations

The iron and steel industry is a major CO2 emitter and an important subject for the implementation of carbon emission reduction goals and tasks. Due to the complex ore composition and low iron grade, vanadium–bearing titanomagnetite smelting in a blast furnace consumes more coke and emits more carbon than in an ordinary blast furnace. Injecting hydrogen–rich gas into blast furnace can not only partially replace coke, but also reduce the carbon emission. Based on the whole furnace and zonal energy and mass balance of blast furnace, the operation window of the blast furnace smelting vanadium–bearing titanomagnetite is established in this study on the premise that the thermal state of the blast furnace is basically unchanged (raceway adiabatic flame temperature and top gas temperature). The effects of different injection amounts of hydrogen–rich gases (shale gas, coke oven gas, and hydrogen) on raceway adiabatic flame temperature and top gas temperature, and the influence of blast temperature and preheating temperature of hydrogen–rich gases on operation window are calculated and analyzed. This study provides a certain theoretical reference for the follow–up practice of hydrogen–rich smelting of vanadium–bearing titanomagnetite in blast furnace.


Introduction
In 2020, China set the goals of "carbon peaking" in 2030 and "carbon neutralization" in 2060 [1]. The iron and steel industry is a major CO 2 emitter and an important subject for the implementation of carbon emission reduction goals and tasks. It is an urgent strategic task to vigorously develop new low-carbon or carbon-free metallurgical technologies, especially for ironmaking systems. In recent years, several researches have been carried out around low-carbon ironmaking technology at home and abroad [2][3][4][5][6], the research focuses on the reduction of iron ore by replacement carbon with hydrogen and, thus, reduction of carbon emission. Although non-blast-furnace ironmaking technology can play a certain role in carbon reduction, it is not the current mainstream ironmaking process; so, the research focus should be on blast furnace ironmaking in the short term. Blast furnace ironmaking is a highly carbon-intensive process, and each ton of hot metal (THM) produced emits roughly 1.8 tons of CO 2 , most of which is consumed in the form of coke. The reducing agent and heating agent of coke can be replaced by injecting hydrogen-rich fuel (with low carbon concentration). The carbon content of shale gas and coke oven gas is about 75% and 23%, respectively, both lower than that of coke (usually 84~91%) [7]. In industry, there have been practical cases for blast furnaces to inject hydrogen-rich gas with CH 4 as the main component, a method widely used in North America, Russia, and other regions; for example, in 2015, 13 blast furnaces in North America (with a total output of 16.6 million THM) only injected natural gas, with an average injection volume of 83 kg/THM, and 12 blast furnaces simultaneously injected natural gas and pulverized coal [8][9][10][11].
In the 1960s and 1970s, Chongqing Iron and Steel Corporation had the practical experience of injecting natural gas into the furnace, but it was abandoned because the price of natural gas was much higher than that of pulverized coal at that time and the economic benefit was low. Now, huge reserves of shale gas resources have been explored in Southwest China, and with the continuous development of mining technology, the cost continues to decline, which provides a practical basis for injecting shale gas into the blast furnace. At the same time, China has fully mastered the smelting technology of vanadium-bearing titanomagnetite in blast furnace and formed a unique technical standard and system. However, at present, the research on smelting vanadium-bearing titanomagnetite by hydrogen-rich or full hydrogen blast furnace combined with shale gas has not been explored in-depth. In order to accurately evaluate the feasibility of hydrogen-rich gas injection into a vanadium-bearing titanomagnetite blast furnace, the current study took a 1750 m 3 blast furnace in an iron and steel plant in Panxi area as the research object and established a mathematical model of blast furnace operation window based on the energy and mass balance of blast furnace by calculating and analyzing the reasonable operation window of shale gas, coke oven gas, and hydrogen. It provides a certain theoretical reference for the follow-up practice of hydrogen-rich smelting of vanadium-bearing titanomagnetite in blast furnace.

Mathematical Model of Energy and Mass Balance in Blast Furnace
According to the whole furnace and zonal heat balance [12][13][14], combined with the characteristics of hydrogen-rich gas injection, the mathematical model of energy and mass balance of hydrogen-rich gas injection was established. The preconditions for establishing energy and mass balance are as follows: (1) The types of hydrocarbon combination in shale gas are complex, but the main component is CH 4 , and the cracking process of other alkanes is similar to that of methane. Therefore, in order to simplify the calculation process, the alkane components are converted into CH 4 in this study. (2) CH 4 in the hydrogen-rich gas undergoes cracking and combustion reaction at the tuyeres to generate CO and H 2 . (3) The relative error of material balance is within 0.1%, which ensures the accuracy of solving the gas composition equation in each zone. (4) The heat loss of the whole furnace and area is 5~8%.
The establishment method of energy and mass balance is as follows: (1) Mass balance calculation. The iron ore consumption was obtained according to the Fe balance. Before injecting hydrogen-rich gas, the blast volume and top gas composition were measured, and the direct reduction degree could be obtained according to actual production data. Then, when hydrogen-rich gas injection was simulated, the blast volume was obtained according to the oxygen balance; the furnace hearth and top gas compositions were calculated according to the C, H, and N balances; the direct reduction degree after injecting hydrogen-rich gas could be approximately calculated by the empirical formula as follows [14]: (2) Heat balance calculation. The heat balance of the whole furnace was calculated based on the Gass Law, without considering the actual reaction process in the blast furnace, and with the material feeding state as the starting point and the output state as the end point; the zonal heat balance was based on the actual reaction process, taking into account the main energy sources and heat consumption in the blast furnace. The thermal effects of all chemical reactions are considered in the heat balance calculation of blast furnace, which is based on the database data of FactSage 8.1 software and related literatures. The thermal effects of different chemical reactions are different at different temperatures, among which some major chemical reactions and their thermal effects are as follows:

Thermodynamics Analysis of Reduction of Vanadium and Titanium Oxides
Considering that vanadium-bearing titanomagnetite contains a lot of titanium oxides and vanadium oxides, it is necessary to research the reduction of titanium oxides and vanadium oxides in the process of smelting vanadium-bearing titanomagnetite by hydrogenrich or full hydrogen blast furnace combined with shale gas. Figure 1 depicts the Gibbs' free energy change as a function of temperature in the reduction of V 2 O 5 and TiO 2 with different reducing agents (CO, H 2 , and C). It can be found that the reduction of V 2 O 5 and TiO 2 is carried out in stages: Compared with CO and H 2 , vanadium oxide and titanium oxide can be easily reduced by carbon at high temperature under standard conditions. Affected by the complex conditions of blast furnace hearth, the contents of vanadium and titanium in hot metal are 0.31% and 0.19%, respectively, according to the actual smelting data of vanadium titanium blast furnace, which are mainly obtained by direct reduction.

Raceway Adiabatic Flame Temperature
The raceway adiabatic flame temperature (RAFT) was calculated by assuming that the combustion products are CO, H2, and N2 only [15], considering the reaction of the blast air and injectants with carbon from coke at 1550 °C and the effects of ash−forming oxides, and the physical heat contained in the fuel and blast can be transferred to the combustion products after deducting the heat consumption in the combustion zone. With the in−depth study of physical and chemical processes in blast furnace, the RAFT calculation model [16][17][18][19] has been continuously modified and improved, and the following points were also considered based on the previous RAFT calculation: the influence of temperature on heat capacity of raw fuel, blast, and hydrogen−rich gas; the heat consumption of heating up, melting, and SiO2 gasification of ash after combustion; the sensible heat and reaction heat of pulverized coal and hydrogen−rich gas.
According to the injection characteristics of hydrogen−rich blast furnace, the raceway adiabatic flame temperature (RAFT) was modified from two aspects: firstly, the decomposition heat of methane was calculated by using FactSage 8.1 software; secondly, the combustion rate of shale gas was obtained by simulating the carbon balance of blast furnace.

Raceway Adiabatic Flame Temperature
The raceway adiabatic flame temperature (RAFT) was calculated by assuming that the combustion products are CO, H 2 , and N 2 only [15], considering the reaction of the blast air and injectants with carbon from coke at 1550 • C and the effects of ash-forming oxides, and the physical heat contained in the fuel and blast can be transferred to the combustion products after deducting the heat consumption in the combustion zone. With the indepth study of physical and chemical processes in blast furnace, the RAFT calculation model [16][17][18][19] has been continuously modified and improved, and the following points were also considered based on the previous RAFT calculation: the influence of temperature on heat capacity of raw fuel, blast, and hydrogen-rich gas; the heat consumption of heating up, melting, and SiO 2 gasification of ash after combustion; the sensible heat and reaction heat of pulverized coal and hydrogen-rich gas.
According to the injection characteristics of hydrogen-rich blast furnace, the raceway adiabatic flame temperature (RAFT) was modified from two aspects: firstly, the decomposition heat of methane was calculated by using FactSage 8.1 software; secondly, the combustion rate of shale gas was obtained by simulating the carbon balance of blast furnace.
RAFT calculation formula is

Top Gas Temperature
The top gas temperature (TGT) with certain heat loss conditions was calculated based on the zonal heat balance [14]. The influence of injecting hydrogen-rich gas and direct reduction degree on the top gas temperature was considered in the calculation, which was reflected in the items q hgas , q i , and V tgas . In addition, the effects of temperature and gas composition on heat capacity (C tgas ) were also considered.
The calculation formula is as follows: In the middle and lower part of the blast furnace, the reducing ability of carbon is much greater than that of CO and H 2 due to the higher temperature (gas temperature greater than 1000 • C and burden temperature greater than 950 • C, based on literature data), suggesting that the reduction reaction in this high temperature area is mainly direct reduction [14,20]. Therefore, since the direct reduction degree will directly affect the amount of gas in the high temperature area, the direct reduction enthalpy is considered in q hgas . The calculation formula of each term of the TGT equation is as follows: It is worth noting that the individual heat in these equations varies with the type and injection rate of hydrogen-rich gas. Based on the energy mass balance model and TGT calculation model, a series of calculations (under the condition that the thermal state of the blast furnace is kept stable) can be used to obtain the relevant parameter values with injecting different hydrogen-rich gases, as shown in Tables A1-A3. The variation law of the top gas composition, individual heats, and direct reduction degree with the injection rate of hydrogen-rich gas is basically consistent with the calculation in the reference [14].

Mathematical Model of Operation Window
When hydrogen-rich gas is injected, the amount of hearth gas increases, the RAFT decreases, and the TGT increases. The ranges of RAFT and TGT can be adjusted by increasing the blast oxygen enrichment rate. RAFT is an important parameter for the blast furnace operator to judge the thermal state of the hearth. It determines the initial gas temperature of the hearth, thus affecting the heat transfer, reduction, slagging, desulfurization, hot metal temperature, composition, etc. When hydrogen-rich fuel is injected, the RAFT becomes low, which reduces the replacement ratio, and with the increase in hydrogen-rich fuel consumption, the furnace condition may even deteriorate. Too high RAFT will also cause the initial gas volume expansion in the hearth, and a large amount of SiO volatilizes, which increases the resistance to the charge column, affects the decline of the charge, and even makes the blast furnace difficult to run and suspend. Too-high or too-low TGT is also unfavorable to the production of blast furnace. Too-high TGT will increase the heat loss in the furnace, increase the coke ratio, and shorten the service life of furnace top charging equipment and bag filter. Too-low TGT will cause the gas temperature to be lower than the dew point, further causing dust to block the filter cloth, making it lose its filtering function and greatly reducing the dust removal efficiency. Therefore, in order to discuss the operation constraints of injecting hydrogen-rich gas, the allowable injection amount of hydrogen-rich gas (Abscissa) and blast oxygen enrichment rate (Ordinate) are often used to describe the feasible operation window on the premise of maintaining the thermal state of blast furnace (reasonable RAFT and TGT) [15,21,22]. The calculation process of the operation window is shown in the Figure 2. According to the previous study [23], the RAFT is controlled at 2050~2200 °C during pulverized coal injection and 1900 °C for natural gas injection. The lower limit of RAFT can be as low as 1704 °C when natural gas injection is recommended for some North American blast furnaces [24]. K. S. Abdel Halim et al. [25] studied the variation range of RAFT when injecting natural gas into four blast furnaces in Russia and Egypt. After injecting natural gas, RAFT showed a downward trend. For example, when SSC 1# blast furnace reached the maximum injection volume (158 m 3 /THM), the RAFT decreased by more than 200 °C. The difference in RAFT calculation methods of blast furnaces in different iron and steel plants lead to the different blast parameters, and the difference of RAFT can reach 150 °C [26]. Therefore, according to the on−site production parameters of a 1750 m 3 blast furnace smelting vanadium−bearing titanomagnetite in Panxi, China, the RAFT is calculated to be 2150~2250 °C; in combination with the charge properties, slag characteristics, top equipment, and dust removal restrictions, this study sets the operation windows of TGT in the range 110~180 °C [27] and the RAFT in the ranges of 2000~2100 °C and 2100~2200 °C.

Calculation Conditions
Due to the large amount of calculation data of energy and mass balance, we only list the main parameters used for calculation. According to the 1750 m 3 blast furnace, the main production parameters under normal generation conditions (only PCI) are shown in Table  A4 and the composition of mixed ore is shown in Table A5. The compositions of According to the previous study [23], the RAFT is controlled at 2050~2200 • C during pulverized coal injection and 1900 • C for natural gas injection. The lower limit of RAFT can be as low as 1704 • C when natural gas injection is recommended for some North American blast furnaces [24]. K. S. Abdel Halim et al. [25] studied the variation range of RAFT when injecting natural gas into four blast furnaces in Russia and Egypt. After injecting natural gas, RAFT showed a downward trend. For example, when SSC 1# blast furnace reached the maximum injection volume (158 m 3 /THM), the RAFT decreased by more than 200 • C. The difference in RAFT calculation methods of blast furnaces in different iron and steel plants lead to the different blast parameters, and the difference of RAFT can reach 150 • C [26]. Therefore, according to the on-site production parameters of a 1750 m 3 blast furnace smelting vanadium-bearing titanomagnetite in Panxi, China, the RAFT is calculated to be 2150~2250 • C; in combination with the charge properties, slag characteristics, top equipment, and dust removal restrictions, this study sets the operation windows of TGT in the range 110~180 • C [27] and the RAFT in the ranges of 2000~2100 • C and 2100~2200 • C.

Calculation Conditions
Due to the large amount of calculation data of energy and mass balance, we only list the main parameters used for calculation. According to the 1750 m 3 blast furnace, the main production parameters under normal generation conditions (only PCI) are shown in Table A4 and the composition of mixed ore is shown in Table A5. The compositions of hydrogen-rich gas (shale gas, coke oven gas, and hydrogen), coke, and PCI for calculation are shown in Tables A6-A8, respectively. The composition of hot metal and slag are shown in Tables A9 and A10, respectively. Meanwhile, the calorific value of various fuels can be obtained based on the given composition of coke, coal, and hydrogen-rich gas, as shown in Table A11 [28,29].

Calculation Results and Analysis
Blast furnace smelting is a complex and high-temperature reaction process of gasliquid-solid countercurrent. Replacing coke with hydrogen-rich gas will change the fuel structure of blast furnace smelting, which will have a great impact on the blast furnace smelting. In order to determine the feasibility of injection of hydrogen-rich gas into blast furnace with smelting vanadium-bearing titanomagnetite, the influence of three types of injected hydrogen-rich gases (shale gas, coke oven gas, and hydrogen) on the blast furnace smelting parameters was studied.
Blast furnace smelting is a complex and high-temperature reaction process of gasliquid-solid countercurrent. Replacing coke with hydrogen-rich gas will change the fuel structure of blast furnace smelting, which will have a great impact on blast furnace smelting. In order to determine the feasibility of injection of hydrogen-rich gas into blast furnace with smelting vanadium-bearing titanomagnetite, the influence of three types of injected hydrogen-rich gases (shale gas, coke oven gas, and hydrogen) on the blast furnace smelting parameters was studied. Figure 3 shows the variations of raceway adiabatic flame temperature and coke ratio after injecting hydrogen-rich gas with fixed PCI (108 kg/THM). When only pulverized coal is injected, the RAFT is 2230.3 • C. As each of the three hydrogen-rich gases replaces part of the coke and starts injection, the RAFT decreases at a slower rate with the increase in the injection quantity. However, the effect of different hydrogen-rich gases' injection on the RAFT is obviously different. When the shale gas is injected (Figure 3a), the RAFT decreases at a rate of 4.7~3.4 • C/kg with the increase in injection quantity; the RAFT also drops to 2000 • C with the injection quantity of shale gas reaching 50.4 kg/THM. Figure 3b shows that injecting coke oven gas will decrease the RAFT by 3.7~2.9 • C/kg, while the RAFT drops to 2000 • C with an injection quantity of 62.8 kg/THM. Figure 3c shows that the decrease rate of RAFT is 10.2~7 • C /kg when pure hydrogen is injected into the furnace, and it only needs 22.8 kg/THM of injection for a RAFT of 2000 • C. This indicates that the decrease rate of the RAFT is the fastest when injecting hydrogen, while the decrease rate of the RAFT with injecting coke oven gas is the slowest. In addition, it can be also found that the coke ratio decreases at a gradually slower rate with the uniform increase in three hydrogen-rich gases injection quantity. It means that the substitution capacity of the three hydrogen-rich gases for coke gradually decreases with the increase in injection quantity, which may explain why the decreasing rate of the RAFT gradually slows down.   Figure 4 shows the variations of top gas temperature (TGT) and coke ratio after injecting hydrogen−rich gas with fixed PCI (108 kg/THM). It can be found that the initial TGT is 155.1 °C when only pulverized coal is injected. With the increase in hydrogen−rich gases injection quantity, the TGT shows an upward trend, while the increasing rate of the TGT is different with injecting different hydrogen−rich gases. As shown in Figure 4a, the increasing rate of TGT is about 1.1 °C/kg as shale gas injection quantity increases, while it is 0.8~0.9 °C/kg with coke oven gas injection (Figure 4b). This means that the injection of both hydrogen−rich gases can steadily increase the TGT. However, Figure 4c shows that the increase rate of the TGT with pure hydrogen injection varies from 0.7 to 2.2 °C/kg, which is not as stable as the injection of the other two hydrogen−rich gases. This may be caused by the increase in heat income per ton of iron in the blast furnace, which is due to the rapid decline in the substitution capacity of unit hydrogen injection quantity for coke.  Figure 4 shows the variations of top gas temperature (TGT) and coke ratio after injecting hydrogen-rich gas with fixed PCI (108 kg/THM). It can be found that the initial TGT is 155.1 • C when only pulverized coal is injected. With the increase in hydrogen-rich gases injection quantity, the TGT shows an upward trend, while the increasing rate of the TGT is different with injecting different hydrogen-rich gases. As shown in Figure 4a, the increasing rate of TGT is about 1.1 • C/kg as shale gas injection quantity increases, while it is 0.8~0.9 • C/kg with coke oven gas injection (Figure 4b). This means that the injection of both hydrogen-rich gases can steadily increase the TGT. However, Figure 4c shows that the increase rate of the TGT with pure hydrogen injection varies from 0.7 to 2.2 • C/kg, which is not as stable as the injection of the other two hydrogen-rich gases. This may be caused by the increase in heat income per ton of iron in the blast furnace, which is due to the rapid decline in the substitution capacity of unit hydrogen injection quantity for coke.  Through previous calculations, it has been found that injecting hydrogen−rich gas will decrease the RAFT and increase the TGT. In order to keep the RAFT and TGT within a reasonable range, it is necessary to adjust the blast injection rate and the oxygen enrichment rate in the blast, which is called the operating window (the functional relationship between the oxygen enrichment rate and the blast injection rate). Considering Through previous calculations, it has been found that injecting hydrogen-rich gas will decrease the RAFT and increase the TGT. In order to keep the RAFT and TGT within a reasonable range, it is necessary to adjust the blast injection rate and the oxygen enrichment rate in the blast, which is called the operating window (the functional relationship between the oxygen enrichment rate and the blast injection rate). Considering that the operating window constrained by the RAFT and TGT is also affected by the blast temperature and preheating temperature of injected hydrogen-rich gases in actual production, the influence of different blast temperatures and preheating temperatures on the operating window for hydrogen-rich gas injection is also simulated. Figure 5 shows the calculated operating window for three hydrogen-rich gases at two blast temperatures (at constant preheating temperature of hydrogen-rich gases of 25 • C). The constraint ranges of the TGT and RAFT are 110~180 • C and 2000~2200 • C, respectively. The difference between window 1 and window 2 is the constraint range on RAFT: window 1 is 2100~2200 • C and window 2 is 2000~2100 • C. It can be found that all the operating windows for three hydrogen-rich gases shift to the lower right with the increase in blast temperature, while the shape and area of the operation window have not changed significantly. This indicates that increasing blast temperature at the same injection rate is beneficial to reduce the upper and lower limits of oxygen enrichment rate, while the relative adjustment range of oxygen enrichment rate and injection rate will not change. As the blast temperature increases from 1100 • C to 1200 • C, the injection rate ranges of shale gas, coke oven gas, and hydrogen in their respective integrated operating windows (the constraint range of the RAFT is 2000~2200 • C) change from 0~90 kg/THM, 0~107.2 kg/THM, 0~31.8 kg/THM to 0~90.4 kg/THM, 0~108.3 kg/THM, 0~32.9 kg/THM, respectively, while the oxygen enrichment rates corresponding to their maximum injection rate change from 14.2%, 13.3%, 9% to 12.9%, 12.1%, 8%, respectively. When the injection rate of shale gas, coke oven gas, and hydrogen is 20 kg/THM, the oxygen enrichment rate ranges corresponding to their injection rate change from 3.3~7.2%, 2.9~6.8%, 3.1~7.3% to 2.4~6.1%, 2~5.7%, 2.2~6.2%, respectively. It shows that the increase in blast temperature can slightly increase the injection rate range of the operating window, but can reduce the oxygen enrichment rate of hydrogen-rich gases. In addition, it can be seen from Figure 5 that the calculated operating window for hydrogen injection has the smallest variation range for blast injection rate and oxygen enrichment rate compared with the other two hydrogen-rich gases. This means that injecting hydrogen has the lowest fault tolerance and more stringent operation requirements. Figure 6 shows the calculated operating window for three hydrogen-rich gases at different preheating temperatures (at constant blast temperature of 1100 • C). It can be found that when the gas injected is shale gas (Figure 6a,b) and coke oven gas (Figure 6c,d), the increase in preheating temperature will also cause the operating window to shift slightly to the lower right, demonstrating that increasing preheating temperature at the same injection rate can only slightly reduce the upper and lower limits of oxygen enrichment rate. Preheating hydrogen (Figure 6e,f) has little effect on the operating window, just a slight shift to the right. As the preheating temperature increases from 25 • C to 600 • C, the injection rate ranges of shale gas, coke oven gas, and hydrogen in their respective integrated operating windows change from 0~90 kg/THM, 0~107.2 kg/THM, 0~31.8 kg/THM to 0~94 kg/THM, 0~112.5 kg/THM, 0~32.1 kg/THM, respectively, while the oxygen enrichment rates corresponding to their maximum injection rate change from 14.2%, 13.3%, 9% to 13.7%, 13%, 8.9%, respectively. When the injection rate of shale gas, coke oven gas, and hydrogen is 20 kg/THM, the oxygen enrichment rate ranges corresponding to their injection rate change from 3.3~7.2%, 2.9~6.8%, 3.1~7.3% to 3.1~7%, 2.8~6.7%, 3.1~7.3%, respectively. This shows that the increase in preheating temperature can obviously increase the injection rate range of the operating window, but has little effect on the oxygen enrichment rate. By comparing the operation windows in the literature, it can be found that the calculated results are basically consistent with the literature [14,26]. seen from Figure 5 that the calculated operating window for hydrogen injection has th smallest variation range for blast injection rate and oxygen enrichment rate compare with the other two hydrogen−rich gases. This means that injecting hydrogen has th lowest fault tolerance and more stringent operation requirements.  Figure 6 shows the calculated operating window for three hydrogen−rich gases different preheating temperatures (at constant blast temperature of 1100 °C). It can b found that when the gas injected is shale gas (Figure 6a,b) and coke oven gas (Figure 6c,d the increase in preheating temperature will also cause the operating window to shi operation windows, the operation windows obtained by theoretical calculation can be approximately seen as a parallelogram, which means that too low or too high an injection rate of hydrogen−rich gas will reduce the adjustment range of oxygen enrichment rate, thus reducing the fault tolerance rate of blast furnace operation. However, due to the lack of actual production parameters, it is difficult to make a comparison to determine the best range of injection rates for different hydrogen−rich gases.  In addition, although different hydrogen-rich gas compositions can obtain different operation windows, the operation windows obtained by theoretical calculation can be approximately seen as a parallelogram, which means that too low or too high an injection rate of hydrogen-rich gas will reduce the adjustment range of oxygen enrichment rate, thus reducing the fault tolerance rate of blast furnace operation. However, due to the lack of actual production parameters, it is difficult to make a comparison to determine the best range of injection rates for different hydrogen-rich gases.

Conclusions
The composition of injected hydrogen-rich gas has an impact on RAFT and TGT; different hydrogen-rich gas compositions can obtain different operation windows. Based on the production data of a blast furnace smelting vanadium-bearing titanomagnetite, the energy and mass balance calculation was carried out, and the operation windows of the blast furnace with hydrogen-rich gases injection were determined (the top gas temperature (TGT) was 110-180 • C, and the raceway adiabatic flame temperature (RAFT) was 2000-2200 • C). When the blast temperature increases from 1100 • C to 1200 • C, the theoretical injection quantities of shale gas, coke oven gas, and hydrogen change from 0~90 kg/THM, 0~107.2 kg/THM, 0~31.8 kg/THM to 0~90.4 kg/THM, 0~108.3 kg/THM, 0~32.9 kg/THM, respectively, and the corresponding oxygen enrichment rates change from 1.4~14.2%, 1.4~13.3%, 1.4~9% to 1.3~12.9%, 1.3~12.1%, 1.3~8%, respectively. The theoretical calculation results show that the increase in blast temperature can significantly reduce the oxygen enrichment rate of hydrogen-rich gas, while the increase in preheating temperature can significantly increase the injection quantities. Injection hydrogen-rich gas can greatly reduce the use of coke, thereby reducing greenhouse gas emissions.

Data Availability Statement:
The data presented in this study are available on request from the corresponding author.