A Graphical Tool to Describe the Operating Point of Direct Reduction Shaft Processes
Abstract
:1. Introduction
2. Modeling
2.1. Description of the Midrex NGTM Process
- The top gas (T) exhausted from the shaft furnace, rich in CO and H2, is completely recycled. A total of 30% of the top gas is recycled as a fuel gas (F) in the reformer and the heat recovery unit. The other fraction, called process gas (P), is recycled as a reducing agent.
- The fuel gas (F) is burned in the reformer burners with an additional amount of natural gas and preheated air from the heat recovery unit. The energy from the hot reformer exhaust gases is recovered in a special device.
- The reducing gas is prepared in several steps. First, the process gas (P) is mixed with injected natural gas. This mixture, called the feed gas, is preheated in the heat recovery unit and then injected into the reformer. Cracking takes place in the reformer tubes between CH4, CO2, and H2O. The resulting reformed gas (R) is composed mainly of CO, CO2, H2, and H2O, with a small amount of CH4 remaining (few %).
- The reformed gas (R) is mixed with additional natural gas, called enrichment gas, and pure oxygen, to form the bustle gas (B), whose temperature can reach 900–950 °C.
- The bustle gas (B) is injected into the transition zone located at the mid-height of the shaft. Additional natural gas is also injected into the shaft, at the bottom of the cooling gas loop, and into the transition zone with the bustle gas.
- The shaft furnace is a vertical counter-current gas–solid reactor with a downward flow of iron oxides and an upward flow of hot reducing gas. The iron pellets are then both reduced and carburized. The direct reduced iron achieves a high degree of metallization (92–96%) and a moderate degree of carburization (2–2.5% of carbon in total mass).
2.2. Definition of the Mass Balance Diagram
2.3. Properties of the Mass Balance Diagram
2.3.1. Gas Composition
2.3.2. Gas Mixing
2.3.3. Gas–Solid Reactions in the Shaft Furnace
- Reduction: This reaction corresponds to the transfer of oxygen from the iron-bearing material to the reducing gas. In a counter-current gas–solid configuration, the oxygen variations in the ferrous burden and in the gas are strictly equal (). Consequently, the variation in the gas-specific oxidation is equal to the variation in the burden oxidation degree y, defined by Rist [11] (1963) as the stoichiometric ratio O/Fe of iron oxides:
- Carburization: This reaction describes the transfer of carbon from the gas to the solid. The loss of carbon in the reducing gas is therefore related to the carburization rate of the DRI and can be calculated using Equation (14).
2.4. Thermodynamics Constraints
- If , metallization is thermodynamically impossible: this defines the pre-reduction zone
- If , metallization can occur: this defines the metallization zone.
3. Results
3.1. Reference Midrex NGTM Operating Point
- The humidity of the fuel gas and process gas is 5% and 10%, respectively;
- The composition of natural gas is based on the work of Farhadi et al. [30] (2003);
- The natural gas injected directly into the vessel is assumed to be evenly distributed between the transition zone and the cooling zone;
- The injection of oxygen upstream of the bustle gas (called suroxygenation), although not specified by Sarkar, has been accounted for to maintain mass balance.
3.2. Overview of the Mass Balance Diagram for a Midrex Unit
- Process Gas (node 1):The flow and dry composition of the process gas are measured in the plant using the top gas analyzer, while the moisture content of the gas is estimated based on the vapor saturation pressure derived from the gas temperature. The first node corresponding to the process gas is positioned and labeled P in the diagram. The fuel gas, represented by the node F, is calculated and plotted using the same methodology.
- Feed Gas and Reformed Gas (node 2):The feed gas is a mixture of the process gas with some natural gas injections. The position of the feed gas is easily determined using the method described in Section 2.3.2 on gas mixing. The reformed gas node R overlaps the feed gas node because gaseous chemical reactions do not affect the node position.
- Bustle Gas (node 3):Bustle gas is a mixture of the reformed gas with natural gas and oxygen. The corresponding node B is drawn according to the mixing methodology developed in the Section 2.3.2. Note again that the internal reactions of the gas (combustion with oxygen and water gas shift equilibrium) does not affect the position of the node.
- Reducing Gas (node 4):The reducing gas, represented by node I, corresponds to the theoretical mixtures of all the gases injected into the furnace shaft below the gas–solid reaction zone. It includes the bustle gas, as well as some additional natural gas in the transition zone and in the cooling zone. Optionally, the additional steam accompanying the cooling gas loop and the seal gas can also be considered. For practical reasons, we have added the node TZ, which corresponds to the theoretical mixture of the bustle gas with the natural gas injected in the transition zone.
- Top Gas (node 5):The methodology described in Figure 5 is used to define the graphical position of the top gas. First, the decarburized gas, shown as C, corresponds to the composition of the reducing gas after complete carburization and before reduction. We simplify the process description by assuming that carburization and reduction take place in separate zones. Carburization is achieved by in situ methane reforming and the Boudouard reaction. It is generally accepted that both reactions occur in the lower part of the metallization zone, in the transition zone, favored by the high local temperature and the presence of metallic iron acting as a catalyst as explained by Shams and Moazeni [25] (2015) and simulated by Hamadeh et al. [26] (2018).Finally, the top gas is deduced from the decarburized gas. Based on the reduction rate from pellets to DRI, we determine the position of the top gas node T, as well as the pure wustite point node W’ defined in Section 2.4.
3.3. Application 1: Calculation of the Gas Moisture
- We start with the diagram describing the operating point of the industrial unit, as shown in Section 3.2 and illustrated in Figure 8. At this point, we recover the positions of all the wet gases on the diagram.
- The gas analyzer gives the composition of the dry part of the gas. From this, we calculate the associated oxidation degree and draw the “dry line” on the diagram, representing all the dry gases potentially associated with this composition:
- In addition, we can perform a vector decomposition of the wet gas using the gas mixing methodology explained in Section 2.3.2. We thus define the “wet line” as the line with slope xH2O equal to 1, passing through the node associated with the wet gas:
3.4. Application 2: Mass Balance of the Top Gas Recycling
3.5. Application 3: Estimation of the Instantaneous Production Rate
- 0.
- First, we draw the operating line of the Midrex NGTM working point from the process gas to the decarburized gas (nodes: F, P, R, B, TZ, I and C). At this point, the top gas is not shown on the diagram.
- 1.
- We plot the node Td associated with the dry top gas from the mass balance on the top gas recycling:
- 2.
- We draw the wet line of the top gas, passing through Td with slope = 1.
- 3.
- We also draw the line associated with the reduction path: it is a vertical line () passing through C.
- 4.
- We place the node associated with the top gas T at the intersection of the two lines drawn in steps 2 and 3.
4. Conclusions
- When the composition of the (dry) gas is measured, it is possible to quantify its moisture.
- A graphical mass balance of the top gas recycling allows diagnosis of the consistency of the measurements of flowmeters and gas analyzers.
- The evaluation of the instantaneous production of DRI, as of the instantaneous production of hot metal for the blast furnace, is obtained from a graphical mass balance on the reducing gas at the inlet and outlet of the shaft.
Funding
Conflicts of Interest
Abbreviations
Greek letters | |
specific consumption | |
specific oxidation | |
molar flux (mol·m−3·s−1) | |
Latin letters | |
volumetric fraction of the molecule j in a mix gas | |
local molar concentration of the atomic element i (mol·m−3) | |
molar mass of the element i (g·mol−1) | |
number of atoms i in the molecule j | |
mass flow rate (kg·s−1) | |
volumetric flow rate (m3·s−1) | |
molar volume of the gas (m3·mol−1) | |
w | mass fraction |
gas oxidation degree (derived definition) | |
gas oxidation degree (original definition) | |
y | burden oxidation degree |
superscript | |
gas | related to gaseous element |
s | related to solid element (iron bearing material) |
subscript | |
i | atomic element (C, H, O or Fe) |
References
- World Steel Association. 2002 World Steel in Figures. Available online: https://worldsteel.org/wp-content/uploads/2002-World-Steel-in-Figures.pdf (accessed on 12 June 2023).
- World Steel Association. 2022 World Steel in Figures. Available online: https://worldsteel.org/publications/bookshop/world-steel-in-figures-2023/ (accessed on 12 June 2023).
- Patisson, F.; Mirgaux, O. Hydrogen Ironmaking: How It Works. Metals 2020, 10, 922. [Google Scholar] [CrossRef]
- Patisson, F.; Mirgaux, O.; Birat, J.P. Hydrogen steelmaking. Part 1: Physical chemistry and process metallurgy. Matér. Tech. 2021, 109, 303. [Google Scholar] [CrossRef]
- Ma, Y.; Souza Filho, I.R.; Bai, Y.; Schenk, J.; Patisson, F.; Beck, A.; van Bokhoven, J.A.; Willinger, M.G.; Li, K.; Xie, D.; et al. Hierarchical nature of hydrogen-based direct reduction of iron oxides. Scr. Mater. 2022, 213, 114571. [Google Scholar] [CrossRef]
- Rechberger, K.; Spanlang, A.; Sasiain Conde, A.; Wolfmeir, H.; Harris, C. Green Hydrogen-Based Direct Reduction for Low-Carbon Steelmaking. Steel Res. Int. 2020, 91, 2000110. Available online: https://onlinelibrary.wiley.com/doi/pdf/10.1002/srin.202000110 (accessed on 12 December 2021). [CrossRef]
- Birat, J.P.; Patisson, F.; Mirgaux, O. Hydrogen steelmaking, part 2: Competition with other net-zero steelmaking solutions–geopolitical issues. Matér. Tech. 2021, 109, 307. [Google Scholar] [CrossRef]
- Raabe, D.; Tasan, C.C.; Olivetti, E.A. Strategies for improving the sustainability of structural metals. Nature 2019, 575, 64–74. [Google Scholar] [CrossRef]
- Ghosh, A.; Chatterjee, A. Ironmaking and Steelmaking: Theory and Practice; PHI Learning Pvt. Ltd.: New Delhi, India, 2008. [Google Scholar]
- Nicolle, R. Blast furnace theory: The contribution of the French School. Revue de Métallurgie 2005, 102, 199–217. [Google Scholar] [CrossRef]
- Rist, A.; Bonnivard, G. Réduction d’un lit d’oxydes de fer par un gaz. Revue de Métallurgie 1963, 60, 23–37. [Google Scholar] [CrossRef]
- Rist, A.; Bonnivard, G. Réduction d’un lit d’oxydes de fer par un gaz. Revue de Métallurgie 1966, 63, 296–312. [Google Scholar] [CrossRef]
- Rist, A.; Meysson, N. Étude graphique de la marche du haut fourneau avec vent humide et injections aux tuyères. Revue de Métallurgie 1965, 62, 995–1040. [Google Scholar] [CrossRef]
- Rist, A.; Meysson, N. Recherche graphique de la mise au mille minimale du haut fourneau a faible température de vent. Revue de Métallurgie 1964, 61, 121–146. [Google Scholar] [CrossRef]
- Nicolle, R.; Capelani, R.; Guillot, J.B.; Rist, A. La réduction dans le haut fourneau et la recherche des très faibles mises au mille de coke. Revue de Métallurgie 1980, 77, 769–779. [Google Scholar] [CrossRef]
- Michard, J. Contribution à l’étude du haut fourneau chargé en aggloméré. Les phénomènes de la cuve. Revue de Métallurgie 1961, 58, 41–46. [Google Scholar]
- Rist, A. Analogue Diagrams for Process Metallurgists. ISIJ Int. 1992, 32, 1034–1043. [Google Scholar] [CrossRef]
- Bailera, M.; Nakagaki, T.; Kataoka, R. Revisiting the Rist diagram for predicting operating conditions in blast furnaces with multiple injections. Open Res. Eur. 2021. [Google Scholar] [CrossRef]
- Bailera, M.; Nakagaki, T.; Kataoka, R. Extending the Operating Line Methodology to Consider Shaft and Preheating Injections in Blast Furnaces. ISIJ Int. 2022, 62, 2454–2465. [Google Scholar] [CrossRef]
- David, S.F.; David, F.F.; Nascimento, R.C. An analysis of pulverized coal injection in blast furnaces by rist diagram/Uma análise sobre a injeção de carvão pulverizado em altos-fornos por meio do diagrama de rist. Braz. J. Dev. 2020, 6, 88911–88929. [Google Scholar] [CrossRef]
- Spanlang, A.; Wukovits, W.; Weiss, B. Development of a Blast Furnace Model with Thermodynamic Process Depiction by Means of the Rist Operating Diagram. Chem. Eng. Trans. 2020, 165, 243–247. [Google Scholar] [CrossRef]
- Orre, J.; Ökvist, L.S.; Bodén, A.; Björkman, B. Understanding of Blast Furnace Performance with Biomass Introduction. Minerals 2021, 11, 157. [Google Scholar] [CrossRef]
- Quatravaux, T.; Barros, J.; Gardin, P.; Lucena, G. Adaptation of the Rist Operating Diagram as a Graphical Tool for the Direct Reduction Shaft. Metals 2021, 11, 1953. [Google Scholar] [CrossRef]
- Parisi, D.R.; Laborde, M.A. Modeling of counter current moving bed gas-solid reactor used in direct reduction of iron ore. Chem. Eng. J. 2004, 104, 35–43. [Google Scholar] [CrossRef]
- Shams, A.; Moazeni, F. Modeling and Simulation of the MIDREX Shaft Furnace: Reduction, Transition and Cooling Zones. JOM 2015, 67, 2681–2689. [Google Scholar] [CrossRef]
- Hamadeh, H.; Mirgaux, O.; Patisson, F. Detailed Modeling of the Direct Reduction of Iron Ore in a Shaft Furnace. Materials 2018, 11, 1865. [Google Scholar] [CrossRef] [PubMed]
- Takenaka, Y.; Kimura, Y.; Narita, K.; Kaneko, D. Mathematical model of direct reduction shaft furnace and its application to actual operations of a model plant. Comput. Chem. Eng. 1986, 10, 67–75. [Google Scholar] [CrossRef]
- Sarkar, S.; Bhattacharya, R.; Roy, G.G.; Sen, P.K. Modeling MIDREX Based Process Configurations for Energy and Emission Analysis. Steel Res. Int. 2018, 89, 1700248. Available online: https://onlinelibrary.wiley.com/doi/pdf/10.1002/srin.201700248 (accessed on 12 December 2021). [CrossRef]
- Mirzajani, A.; Ebrahim, H.A.; Nouri, S.M.M. Simulation of a direct reduction moving bed reactor using a three interface model. Braz. J. Chem. Eng. 2018, 35, 1019–1028. [Google Scholar] [CrossRef]
- Farhadi, F.; Motemed Hashemi, M.Y.; Bahrami Babaheidari, M. Modelling and simulation of syngas unit in large scale direct reduction plant. Ironmak. Steelmak. 2003, 30, 18–24. [Google Scholar] [CrossRef]
Gas | Flow Rate (kNm3·h−1) |
---|---|
Fuel gas | 38.3 |
Process gas | 78.5 |
Process natural gas | 18.0 |
Suroxygenation | 2.84 |
Enrichment natural gas | 0.0 |
Transition zone natural gas | 1.3 |
Cooling zone natural gas | 1.3 |
Composition | Bustle Gas | Dry Top Gas | Natural Gas |
---|---|---|---|
CO (%) | 35.0 | 21.4 | 0.0 |
CO2 (%) | 2.0 | 26.3 | 0.0 |
H2 (%) | 55.0 | 49.4 | 0.0 |
H2O (%) | 6.0 | 0.0 | 0.0 |
CH4 (%) | 1.0 | 1.7 | 87.6 |
C2H6 (%) | 0.0 | 0.0 | 8.2 |
C3H8 (%) | 0.0 | 0.0 | 3.1 |
C4H10 (%) | 0.0 | 0.0 | 1.1 |
N2 (%) | 1.0 | 0.0 | 0.0 |
x | 0.069 | 0.242 | −0.342 |
production | 100 t·h−1 |
metallization rate | 92.7% |
carburization rate | 2.2% |
gangue | 5.8% |
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content. |
© 2023 by the author. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).
Share and Cite
Quatravaux, T. A Graphical Tool to Describe the Operating Point of Direct Reduction Shaft Processes. Metals 2023, 13, 1568. https://doi.org/10.3390/met13091568
Quatravaux T. A Graphical Tool to Describe the Operating Point of Direct Reduction Shaft Processes. Metals. 2023; 13(9):1568. https://doi.org/10.3390/met13091568
Chicago/Turabian StyleQuatravaux, Thibault. 2023. "A Graphical Tool to Describe the Operating Point of Direct Reduction Shaft Processes" Metals 13, no. 9: 1568. https://doi.org/10.3390/met13091568
APA StyleQuatravaux, T. (2023). A Graphical Tool to Describe the Operating Point of Direct Reduction Shaft Processes. Metals, 13(9), 1568. https://doi.org/10.3390/met13091568