Next Article in Journal
First-Principle Investigation into Mechanical Properties of Al6Mg1Zr1 under Uniaxial Tension Strain on the Basis of Density Functional Theory
Next Article in Special Issue
Numerical Study on Heat Transfer Characteristic of Hot Metal Transportation before EAF Steelmaking Process
Previous Article in Journal
Orientation Dependence of B19’-Martensite Reorientation Stress and Yield Stress in TiNi Single Crystals
Previous Article in Special Issue
Numerical Simulation of Heat and Mass Transfer Behavior during Iron Ore Sintering: A Review
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Graphical Tool to Describe the Operating Point of Direct Reduction Shaft Processes

by
Thibault Quatravaux
Institut Jean Lamour, UMR 7198 Université de Lorraine/CNRS, Campus Artem, 2 allée André Guinier-BP 50840, LabEx DAMAS, CEDEX 01, F-54011 Nancy, France
Formerly ArcelorMittal R&D Maizières Process.
Metals 2023, 13(9), 1568; https://doi.org/10.3390/met13091568
Submission received: 4 August 2023 / Revised: 23 August 2023 / Accepted: 30 August 2023 / Published: 7 September 2023
(This article belongs to the Special Issue Advances in Ironmaking and Steelmaking Processes (2nd Edition))

Abstract

:
This article presents a new graphical tool for direct reduction shaft processes inspired by the Rist diagram developed for blast furnaces. The tool represents gas flows using vectors, with specific consumption and specific oxidation as components to indicate gas/iron ratios. Key features include consideration of gas chemical composition for vector directions, easy visual representation of gas mixtures, as well as reduction and carburization rates of direct reduced iron (DRI). The tool also includes thermodynamic conditions for reduction from the Chaudron diagram, analogous to the Rist diagram. Several practical applications are presented, including quantifying gas moisture, evaluating the measurement consistency of flowmeters and gas analyzers in top gas recycling, and evaluating instantaneous DRI production by analyzing reducing gas at the inlet and outlet of the shaft. This graphical tool could be useful for production teams to monitor and optimize process flows and promote understanding among students, engineers, technicians, and operators. Its potential for online use further enhances its practical value. As a result, the tool is of significant academic and industrial interest in improving process efficiency and optimization.

1. Introduction

The production of iron by direct reduction is the main alternative to traditional blast furnace iron production. Over the past two decades, this process has gained considerable popularity. According to the World Steel Association [1,2], direct reduced iron (DRI) now accounts for 10% of total iron production, a significant increase from 1.5% in the 1980s, and reached 125 Mt in 2022, about three times the annual production in the early 2000s. Many current research efforts, both industrial and academic, are devoted to the study of the kinetics of direct reduction processes in which the reducing agent is wholly or partially composed of hydrogen [3,4,5,6]. This path is emerging as one of the most promising methods considered in the strategy for decarbonization of the steel industry [7,8].
Among the direct reduction processes, the gas-based shaft furnace is the most widely used technology, largely represented by the Midrex NGTM and Tenova HYL processes, which use natural gas as the reducing agent for iron oxides, as explained by Ghosh and Chatterjee [9] (2008).
As direct reduction is a relatively new family of processes in ironmaking, with the first reactors being less than 50 years old, there is little literature devoted to a detailed understanding of its operating point. It is interesting to look at the recent evolution of blast furnace operating points over the last few decades, which is the result of a series of complementary research studies carried out between the 1950s and the 1970s, which led to a paradigm shift in the approach to this process [10].
Among these notable studies, Rist and Bonnivard [11] (1963) introduced a graphical representation of the operating point of a blast furnace, incorporating a model of staggered iron oxide reduction mechanisms in a theoretical counter-current gas–solid reactor setup. Rist and Bonnivard [12] (1966) further refined their approach by considering the specific operating conditions of the blast furnace. Several studies are based on the application of the diagram to the optimization of the operating point, such as the effect of wind injection conditions at the tuyeres by Rist and Meysson [13] (1965), or the determination of the minimum coke yield by the same authors [14] (1964), later also studied by Nicolle et al. [15] (1980).
This graphical tool provided visual support for the blast furnace modeling developed by Michard [16] (1961), which was based on the theoretical work provided by Kitaev in the 1950s. After 30 years of industrial and academic applications, Rist [17] (1992) concluded that this approach demonstrated the value of a synthetic description of the entire reactor system. The pedagogical benefits of this methodology were obvious, facilitating process understanding and the determination of likely operating points through an intuitive visual approach requiring minimum calculations. For these reasons, it remains a topical tool for the study and optimization of blast furnace operation [18,19,20], as well as for the evolution of the process towards decarbonization [21], by hydrogen injection or the use of biomass [22].
Quatravaux et al. [23] (2021) have recently successfully adapted the Rist operating diagram to direct reduction shaft processes, specifically for the Midrex NGTM process. They have shown that this type of approach remains relevant among the many tools currently available for process study, complementing the various numerical solutions, whether via 1D coupled reactor models provided by Parisi and Laborde [24] (2004) or Shams and Moazeni [25] (2015), as well as simulations of coupled phenomena computed by the finite volume method developed by Hamadeh et al. [26] (2018).
However, this diagram is only a partial description of the DRI process, as it deals only with the reduction zone of the shaft furnace. It hides part of the process including top gas recycling and reducing gas preparation.
Therefore, this manuscript presents a novel graphical tool developed specifically for direct reduction shaft processes which is applied to the description of the Midrex NGTM process. It provides a global approach to an industrial plant and incorporates many of the concepts found in the Rist operating diagram.

2. Modeling

2.1. Description of the Midrex NGTM Process

The main differences between direct reduction and blast furnace iron ore reduction are the use of natural gas as a reducing agent and as a heat source with an operating point at temperatures below 1000 °C, without melting the iron-bearing material. The Midrex NGTM process is the most widely used technology in the steel industry for the production of direct reduced iron (DRI), accounting for 60% of total DRI production. The Midrex NGTM process industrial unit is based on the combination of three main components: a shaft furnace, a reformer, and a heat recovery system. This configuration allows for optimum recovery of natural gas by fully recycling the gases from the shaft furnace into the process. Figure 1 shows the operating principle of the Midrex NGTM process.
The description of the operating point of this process is addressed through the main steps identified and described below:
  • 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

We will consider the characteristics of a gas flow (G) at each point of the installation. We will describe this gas in terms of adimensional variables that quantify some stoichiometric relationships between this gas flow and that of the iron-bearing material passing through the shaft furnace.
The first variable μ , called specific consumption by Rist, represents the stoichiometric ratio between the carbon and hydrogen gas flow and the iron flow passing through the shaft furnace, defined by (1):
μ = ϕ C g a s + ϕ H g a s 2 ϕ F e s
where ϕ C g a s and ϕ H g a s are the molar flux of carbon and hydrogen calculated with the following relation:
ϕ i g a s = Q v g a s V m × m o l e c u l e j n i j · % j
Q v g a s is the volumetric flow rate of the gas, V m is the molar volume of the gas, % j is the volumetric fraction of molecule j in the gas, and n i j is the number of atoms i in the molecule j.
The gas is mainly composed of CO, CO2, H2, H2O, and hydrocarbons CmHn; the molar gas fluxes can be approximated by the following relations:
ϕ C g a s = Q v g a s V m × % C O + % C O 2 + m · % C m H n
ϕ H g a s = Q v g a s V m × 2 % H 2 + 2 % H 2 O + n · % C m H n
ϕ O g a s = Q v g a s V m × % C O + 2 % C O 2 + % H 2 O
ϕ F e s is the molar flux of iron from solid s into the furnace shaft. It is obtained by Equation (6) as a function of the mass flow rate Q m s and the mass fraction of iron w F e s in the solid. The reference solid used in the calculation can be pellets or DRI.
ϕ F e s = Q m s × w F e s M F e
Similarly, the specific oxidation ν is defined as the stoichiometric ratio of carbon and oxygen gas flux to iron flux through the shaft:
ν = ϕ O g a s ϕ C g a s ϕ F e s
Rist defined the oxidation degree of a gas as the following stoichiometric ratio:
x = ϕ O g a s ϕ C g a s ϕ C g a s + ϕ H g a s 2
From this definition, we can deduce:
x = ν μ
Finally, the two variables μ and ν are plotted on a diagram on the x-axis and y-axis, respectively. Any gas G flowing through the direct reduction unit can be plotted on such a graph, as shown in Figure 2.

2.3. Properties of the Mass Balance Diagram

From the definitions given in Section 2.2, we derive the following implications:

2.3.1. Gas Composition

In this diagram, the gas flow (G) can be represented by a vector approach, so that the oxidation degree x corresponds to the direction of the vector O G .
The directions of the main gaseous molecules encountered in direct reduction processes are shown in the diagram in Figure 3. In addition, pure carbon is also mentioned, although it is not a gas, in order to introduce the phenomena of carburization with depletion of carbon in the gas.

2.3.2. Gas Mixing

Both specific consumption and oxidation are extensive variables, so a simple law of additivity can be applied in the context of a gas mixture:
μ m i x = g a s m i x e d μ g a s
ν m i x = g a s m i x e d ν g a s
With a vector approach, the gas mixing process is simply described by Equation (12) and easily interpreted graphically, as shown in Figure 4.
O G m i x = g a s m i x e d O G g a s

2.3.3. Gas–Solid Reactions in the Shaft Furnace

Two main types of gas–solid reactions occur in direct reduction processes:
  • 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 ( δ ϕ O s = δ ϕ O g ). 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:
    δ ν = δ y
  • 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).
    δ ν = δ μ = w C D R I w F e D R I · M F e M C
w C D R I and w F e D R I are the mass fraction of carbon and iron in the DRI.
Finally, Figure 5 shows the graphical interpretation of these reactions in the mass balance diagram.

2.4. Thermodynamics Constraints

Similar to the development of the Rist diagram, it is possible to include thermodynamic constraints on the reduction mechanisms that were highlighted by Chaudron.
As a reminder, Chaudron established equilibrium conditions as a function of temperature between different iron oxides and the gas mixtures H2–H2O, CO–CO2. The diagram thus represents gas–wustite and gas–magnetite thermodynamic equilibria as a function of the composition and temperature of the reducing gas, as shown in Figure 6.
We define the metallization line passing through the origin with a slope equal to x W , which is the gas oxidation degree in equilibrium with Wustite and metallic iron. This value, x W , is the oxidation degree of the so-called “Wustite Point” in the Rist diagram. Quatravaux et al. [23] provided the methodology for calculating x W according to the Chaudron diagram:
x C h a u d r o n = % C O 2 + % H 2 O % C O + % C O 2 + % H 2 + % H 2 O
x C h a u d r o n is the oxidation degree of the reducing gas as defined in the Chaudron diagram formalism.
We then derive the corresponding oxidation degree for the mass balance diagram according to the approach presented by Quatravaux et al. [23], based on the following relationships:
x W = ( 1 α ) · x C h a u d r o n α
α = m + n 2 % C m H n η C H
where:
η C H = % C O + % C O 2 + % H 2 + % H 2 O + m + n 2 % C m H n
α is a coefficient to transpose the equilibrium lines from the Chaudron diagram to the corresponding gas oxidation degree defined by Rist, assuming the presence of hydrocarbons C m H n in the reducing gas.
% C m H n is the CmHn content of the gas in the reduction zone. It is assumed to be the same as the top gas.
The metallization line divides the reduction zone into two distinct parts for which the reduction conditions are thermodynamically different, as shown in Figure 7:
  • If s l o p e > x W , metallization is thermodynamically impossible: this defines the pre-reduction zone
  • If s l o p e < x W , metallization can occur: this defines the metallization zone.
A detailed calculation of the metallization line has been presented by Quatravaux et al. [23]. The slope of the metallization line depends on the temperature and composition of the reducing gas. In the present study, we consider the local temperature in the schaft at the start of metallization. This is estimated to be 800 °C, in agreement with in-situ measurements by Takenaka and Kimura [27]. According to the definition of the Rist diagram, it is also possible to draw the deviation to ideality ω , which is the difference between the iron oxidation degree at the beginning of metallization and the equivalent pure wustite point; the point is assumed to correspond to y = 1.056 for such a temperature, represented by the node W′, as described in Figure 7.
ω must be positive to obtain a reduction in the iron oxides up to the metallization, so W′ must be located in the metallization zone. Otherwise, the reducing power of the gas would be insufficient, and the reduction would stop at the wustite stage.
In conclusion, the metallization line—in particular, ω —is an indicator that allows for visualization of the feasibility of a metallization according to the reducing gas injected in the shaft furnace. It could be useful in the context of a prospective search for new breakthrough operating points.

3. Results

3.1. Reference Midrex NGTM Operating Point

For pedagogical purposes, the theoretical operating point of a Midrex unit as proposed by Sarkar et al. [28] (2018) serves as a reference. The extensive details provided in this reference outweigh the descriptions of operating points taken from various industrial plants, including those in Siderca (Argentina), Gilmore (United States) by Parisi and Laborde [24] (2004), and Khorasan (Iran) by Mirzajani et al. [29] (2018). We have completed the dataset under the following assumptions:
  • 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.
Table 1 lists the gas flow rates, while Table 2 lists their compositions. Data related to the DRI can be found in Table 3.

3.2. Overview of the Mass Balance Diagram for a Midrex Unit

Figure 8 presents the corresponding graphical operating point described above. To facilitate the reading of the diagram, the reader is referred in parallel to the process description in Figure 1.
The operating point of the Midrex NGTM process is drawn through the following main steps:
  • 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

A first application of this diagram is the quantification of gas humidity in Midrex NGTM plants. This is possible if the gas composition is measured regularly, which is the case for the reformed gas, bustle gas, and top gas.
Let us first consider the wet gas as a mixture between its dry part and steam. The dry part of the gas is plotted on the diagram using the following method:
  • 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 x a n a l y z e r and draw the “dry line” on the diagram, representing all the dry gases potentially associated with this composition:
    ν = x a n a l y z e r · μ
  • 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:
    ν = ν w e t + x H 2 O · ( μ μ w e t )
The node associated with the dry fraction of the gas is positioned at the intersection of the two lines defined above. The corresponding graphical construction is shown in Figure 9.
Finally, the humidity of the gas can be calculated directly from the coordinates of the dry part of the gas according to the following relationship:
% H 2 O = 1 1 x w e t 1 + x w e t · ( μ d r y 1 ) ν d r y
Figure 10 shows the plot of the dry portions of the fuel gas, process gas, and top gas. These three nodes are located on the “dry top gas” line determined from the gas analyzer measurements. This shows the importance of the condensation phenomenon during the recycling of the top gas to the fuel and process gas.

3.4. Application 2: Mass Balance of the Top Gas Recycling

A second application of this diagram concerns the mass balance of the top gas recycling. This approach can be used to check the consistency of the process operating point and thus diagnose the quality of the instrumentation measurements.
Since the amount of condensed water from the top gas is not measured, a wet gas mass balance is not possible. Therefore, we apply this approach to the dry fractions of the gas, as defined in the previous section.
We also assume that the flare gas is negligible in the mass balance. In other words, the top gas is completely recycled without any losses. The gas mass balance of the dry top gas is described by the following relationship:
Q v d r y t o p g a s = Q v d r y p r o c e s s g a s + Q v d r y f u e l g a s
On the diagram, we identify the dry portions of the fuel gas, process gas, and top gas by their respective nodes designated by Fd, Pd, and Td according to the graphical methodology shown in Figure 10. The mass balance corresponding to the dry top gas recycling is described by the following simple vector relationship:
O T d = O P d + O F d
This relationship is never strictly verified because of all the measurement uncertainties and the instrumentation drifts. Therefore, we calculate the equilibrium deviation on the top gas recycling mass balance ( D O T G ):
D O T G = 1 O P d + O F d O T d
D O T G is a simple and immediate quantitative criterion to diagnose the consistency of the main process measurements made on the plant, since all mass flow rates and the gas analyzer are involved in this calculation.
We illustrate the application of this criterion in the context of the ideal operating point described in Section 3.1. This is shown in Figure 11, where the relation (23) is perfectly validated. In addition, we use this criterion in the context of a 10% overestimation of the natural gas flow rate injected upstream of the reformer, as shown in Figure 12.

3.5. Application 3: Estimation of the Instantaneous Production Rate

The instantaneous hot metal production in the blast furnace is a parameter calculated from the difference in oxygen per unit of time between the inlet gas at the tuyeres and the outlet top gas at the throat. Since this difference is due to the reduction in the ferrous burden, it is then possible to deduce the associated hot metal production.
This parameter is a key piece of information for controlling the operation of the blast furnace, as it makes it possible to monitor the incoming flows (gas, coke, wind, …) in terms of rate per ton of iron and provides an immediate information on any disturbance of the operating point.
The instantaneous production is therefore a direct consequence of the progress of the iron oxide reduction reaction in the process. This is different from hot metal production, which is measured by casting weight. These two concepts are relatively analogous in the context of operating balances. However, they diverge in the realm of process control, where data on the weight of castings are too inaccurate and too late.
In the case of the Midrex NGTM process, production is confused with the extraction rate of the DRI at the bottom of the shaft furnace. However, for similar reasons and problems as with the blast furnace, it is necessary to make a semantic distinction between (instantaneous) production and extraction.
As with the blast furnace, it is possible to estimate an instantaneous DRI production by calculating the oxygen balance on the reducing gas between the bottom and top of the reduction zone.
The mass balance diagram provides access to this information. This application is interesting for online use of the tool in order to detect fluctuations in the operating point and to make the necessary corrections to the process in a short time.
In order to apply this method, it is necessary to have previously established a material balance for the recycling of the top gas, whose methodology is described in the Section 3.4 and to quantify the DOTG drift.
In addition, we assume that the carburization and metallization rates are stable parameters over time with respect to variations in gas characteristics. This assumption is reasonable because the residence time of the material in the shaft furnace (several hours) is much longer than that of the gas (a few seconds). Therefore, we rely on the chemical measurements performed regularly on DRI.
We first determine the instantaneous reduction rate graphically based on the extraction rate. The corresponding methodology is shown in Figure 13.
The graphical methodology includes the following main steps:
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:
O T d = 1 1 D O T G ( O F d + O P d )
2.
We draw the wet line of the top gas, passing through Td with slope x H 2 O = 1.
3.
We also draw the line associated with the reduction path: it is a vertical line ( i s o μ ) 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.
We define the variation in the oxidation degree of iron obtained using this graphical construction based on the DRI extraction rate δ y e x t r a c t i o n :
δ y e x t r a c t i o n = ν T ν C
The actual variation in the oxidation degree of the ferrous burden, based on instantaneous production, is derived from the following relationship:
δ y p r o d u c t i o n = y p e l l e t y D R I
where y p e l l e t 1.5 and:
y D R I = 1.056 × ( 1 metallization )
The instantaneous production can be easily deduced from the assumed extraction and a correcting coefficient α :
production = extraction × α
The graphical construction described in Figure 13, plotted according to the instantaneous production, would allow the Equations (26) and (27) to be equal. The correction coefficient is calculated with the following equation:
α = δ y e x t r a c t i o n + 2 ( ν C ν I ) δ y p r o d u c t i o n + 2 ( ν C ν I )
It should be noted that the vector I C corresponds to the rate of carburization of the DRI. Therefore, its length is independent of the iron molar flux chosen for the graphical representation.

4. Conclusions

In this paper, we have presented a new graphical tool developed for the direct reduction shaft processes, inspired by the Rist diagram originally dedicated to the blast furnace.
Any gas flow in the plant can be represented on this diagram by a vector whose components, called specific consumption and specific oxidation, represent certain stoichiometric gas/iron ratios.
This particular graphical tool has several unique features. First, the chemical composition of each individual gas is responsible for determining the direction of the associated vector. In addition, a simple vector addition can accurately represent successive gas mixtures. Finally, this tool allows direct graphical representation of the reduction and carburization rates associated with direct reduced iron (DRI).
In addition, the thermodynamic conditions for reduction, highlighted in the Chaudron diagram, are plotted in this tool and are used in a manner strictly equivalent to the Rist diagram.
Several applications of this tool are presented:
  • 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.
As a result, this graphical tool can be useful for production teams to monitor and optimize the process flow. In this sense, it would benefit from being used online.
This tool therefore promotes dialog within a community of students, as well as engineers, technicians, and operators, in order to better understand the process and its optimization.
In this context, it would be of great academic and industrial interest.

Funding

This research received no external funding.

Conflicts of Interest

The author declares no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
Greek letters
μ specific consumption
ν specific oxidation
ϕ molar flux (mol·m−3·s−1)
Latin letters
a j volumetric fraction of the molecule j in a mix gas
c i local molar concentration of the atomic element i (mol·m−3)
M i molar mass of the element i (g·mol−1)
n i j number of atoms i in the molecule j
Q m mass flow rate (kg·s−1)
Q v volumetric flow rate (m3·s−1)
V m molar volume of the gas (m3·mol−1)
wmass fraction
x g a s gas oxidation degree (derived definition)
X g a s gas oxidation degree (original definition)
yburden oxidation degree
superscript
gasrelated to gaseous element
srelated to solid element (iron bearing material)
subscript
iatomic element (C, H, O or Fe)

References

  1. 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).
  2. 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).
  3. Patisson, F.; Mirgaux, O. Hydrogen Ironmaking: How It Works. Metals 2020, 10, 922. [Google Scholar] [CrossRef]
  4. 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]
  5. 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]
  6. 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]
  7. 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]
  8. Raabe, D.; Tasan, C.C.; Olivetti, E.A. Strategies for improving the sustainability of structural metals. Nature 2019, 575, 64–74. [Google Scholar] [CrossRef]
  9. Ghosh, A.; Chatterjee, A. Ironmaking and Steelmaking: Theory and Practice; PHI Learning Pvt. Ltd.: New Delhi, India, 2008. [Google Scholar]
  10. Nicolle, R. Blast furnace theory: The contribution of the French School. Revue de Métallurgie 2005, 102, 199–217. [Google Scholar] [CrossRef]
  11. 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]
  12. 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]
  13. 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]
  14. 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]
  15. 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]
  16. 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]
  17. Rist, A. Analogue Diagrams for Process Metallurgists. ISIJ Int. 1992, 32, 1034–1043. [Google Scholar] [CrossRef]
  18. 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]
  19. 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]
  20. 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]
  21. 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]
  22. 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]
  23. 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]
  24. 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]
  25. 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]
  26. 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]
  27. 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]
  28. 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]
  29. 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]
  30. 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]
Figure 1. Operating principle of the Midrex NGTM process.
Figure 1. Operating principle of the Midrex NGTM process.
Metals 13 01568 g001
Figure 2. Graphical representation of a gas G in the mass balance diagram.
Figure 2. Graphical representation of a gas G in the mass balance diagram.
Metals 13 01568 g002
Figure 3. Directions of the main gases.
Figure 3. Directions of the main gases.
Metals 13 01568 g003
Figure 4. Graphical representation of the gas mixing.
Figure 4. Graphical representation of the gas mixing.
Metals 13 01568 g004
Figure 5. Graphical representation of reduction and carburization.
Figure 5. Graphical representation of reduction and carburization.
Metals 13 01568 g005
Figure 6. Graphical representation of the Chaudron (or Bauer–Glaessner) diagram for a typical reducing gas in the Midrex NGTM shaft.
Figure 6. Graphical representation of the Chaudron (or Bauer–Glaessner) diagram for a typical reducing gas in the Midrex NGTM shaft.
Metals 13 01568 g006
Figure 7. Graphical representation of the metallization zone and the deviation to ideality.
Figure 7. Graphical representation of the metallization zone and the deviation to ideality.
Metals 13 01568 g007
Figure 8. Balance diagram for a Midrex NGTM representative theoretical working point.
Figure 8. Balance diagram for a Midrex NGTM representative theoretical working point.
Metals 13 01568 g008
Figure 9. Graphical determination of the dry part of a gas.
Figure 9. Graphical determination of the dry part of a gas.
Metals 13 01568 g009
Figure 10. Graphical determination of the dry parts of the fuel gas, process gas, and top gas.
Figure 10. Graphical determination of the dry parts of the fuel gas, process gas, and top gas.
Metals 13 01568 g010
Figure 11. Ideal mass balance of the top gas recycling.
Figure 11. Ideal mass balance of the top gas recycling.
Metals 13 01568 g011
Figure 12. Influence of an overestimation of the injected natural gas flow measurement for feed gas preparation (+10%) on the deviation of the top gas recycling mass balance.
Figure 12. Influence of an overestimation of the injected natural gas flow measurement for feed gas preparation (+10%) on the deviation of the top gas recycling mass balance.
Metals 13 01568 g012
Figure 13. Determination of the instantaneous reduction rate.
Figure 13. Determination of the instantaneous reduction rate.
Metals 13 01568 g013
Table 1. Gas volume flow rates.
Table 1. Gas volume flow rates.
GasFlow Rate (kNm3·h−1)
Fuel gas38.3
Process gas78.5
Process natural gas18.0
Suroxygenation2.84
Enrichment natural gas0.0
Transition zone natural gas1.3
Cooling zone natural gas1.3
Table 2. Gas composition.
Table 2. Gas composition.
CompositionBustle GasDry Top GasNatural Gas
CO (%)35.021.40.0
CO2 (%)2.026.30.0
H2 (%)55.049.40.0
H2O (%)6.00.00.0
CH4 (%)1.01.787.6
C2H6 (%)0.00.08.2
C3H8 (%)0.00.03.1
C4H10 (%)0.00.01.1
N2 (%)1.00.00.0
x0.0690.242−0.342
Table 3. Production parameters.
Table 3. Production parameters.
production100 t·h−1
metallization rate92.7%
carburization rate2.2%
gangue5.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.

Share and Cite

MDPI and ACS Style

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

AMA Style

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 Style

Quatravaux, 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 Style

Quatravaux, 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

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop