Reactor reaction calculation device and reactor reaction calculation method

Abstract

To provide a calculation apparatus of a furnace reaction capable of calculating a reaction in a reaction furnace with high accuracy.SOLUTION: A furnace reaction calculation apparatus includes a gas mixing amount calculation unit that calculates a flow rate of an inflow gas that moves from a gas region in the furnace to a bed layer containing a raw material ore in the first gas phase containing a combustion gas in the reaction furnace, an inflow gas distribution unit distributing the inflow gas into a first inflow gas and a second inflow gas, a first bed layer reaction amount calculation unit determining a first bed layer reaction portion by calculating a first bed layer reaction amount contributing to an equilibrium state of the one first bed layer to which a bed layer is distributed with a first reaction rate formula, and a second bed layer reaction amount calculation unit for determining a second bed layer reaction portion by calculating the second bed layer reaction portion that contributes to the equilibrium state of the second bed layer with a second reaction rate formula, a first bed layer equilibrium reaction calculation unit calculating the equilibrium reaction of the first bed layer reaction portion, a second bed layer equilibrium reaction calculation unit calculating an equilibrium reaction of the second bed layer reaction portion, and a gas region equilibrium reaction calculation unit calculating an equilibrium reaction of the first phase.SELECTED DRAWING: Figure 2

Description

[Technical Field]

[0001] The present invention relates to a reactor reaction calculation device and a reactor reaction calculation method. [Background technology]

[0002] As a method for smelting laterite ores (nickel oxide ores) such as limonite ore and saprolite ore, which are types of oxide ores, a dry smelting method is known in which a rotary kiln, a moving hearth furnace, or the like is used to produce ferronickel, which is an alloy mainly composed of iron and nickel.

[0003] In the rotary kiln dry smelting method, raw ore is dried in a rotary dryer to reduce the adhering moisture content to, for example, 15% to 25%, and the dried ore with reduced adhering moisture is then charged into the charging end of the rotary kiln. The dried ore is then heated by the combustion heat of coal supplied from the charging end of the rotary kiln or by a pulverized coal-fired burner or a pulverized coal and heavy oil-fired burner installed at the discharge end of the rotary kiln, thereby drying and firing the dried ore.

[0004] As a dry smelting method using such a rotary kiln, for example, in addition to the combustion heat of coal supplied from the charging end and the combustion heat generated by burning pulverized coal or heavy oil in a burner, there is a method in which the combustion heat generated by burning coal charged midway through the rotary kiln is used to provide the heat necessary for drying and partial reduction of the dried ore.

[0005] For example, Patent Document 1 discloses a method of operating a rotary kiln in which coal is fed into the rotary kiln from a scoop feeder installed midway through the rotary kiln, and dried nickel oxide ore fed into the charging end of the rotary kiln is calcined and partially reduced using the heat of combustion generated by the combustion of fossil fuels in a burner. In this rotary kiln operating method, the coal fed into the scoop feeder is thermally decomposed into volatile matter and fixed carbon, and the volatile matter is discharged from the charging end together with the combustion gases in the furnace, and the fixed carbon is discharged from the discharge end together with the calcined ore produced by drying and reducing the dried ore. [Prior art documents] [Patent documents]

[0006] [Patent Document 1] Patent No. 5967616 Summary of the Invention [Problem to be solved by the invention]

[0007] Here, in known operating methods such as the rotary kiln operating method described in Patent Document 1, the bed layer is treated as a single type of layer, and differences in the reaction rate of the raw ore within the bed layer are not considered. The raw ore supplied to a rotary kiln typically has a particle size distribution and differs in density, composition, etc., resulting in differences in reaction rate due to differences in the particle size, density, composition, etc. of the raw ore. Furthermore, when multiple types of raw ore are supplied, the porosity, composition, etc. of the raw ore vary depending on the type of raw ore, resulting in differences in reaction rate. Therefore, when simulating the behavior of raw ore in a reactor such as a rotary kiln, it is necessary to accurately calculate the reaction within the furnace, taking into account differences in reaction rate due to differences in the physical properties of the supplied raw ore, such as particle size, density, composition, and type.

[0008] An object of one aspect of the present invention is to provide a reactor reaction calculation device that can calculate the reaction in a reactor with high accuracy. [Means for solving the problem]

[0009] One aspect of the apparatus for calculating reactions in a furnace according to the present invention is an apparatus for calculating reactions in a furnace, in which raw material ore is supplied from one end of a reactor and moved toward the other end while the raw material ore is brought into contact with combustion gas supplied from the other end, thereby drying and reducing the raw material ore, a gas mixing amount calculation unit that calculates the flow rate of an inflow gas that moves from the gas region in the reactor to a bed layer containing the raw ore, among a first gas phase containing the combustion gas that flows through a gas region in the reactor; an inflow gas distribution unit that distributes the inflow gas into a first inflow gas and a second inflow gas; a first bed layer reaction amount calculation unit that calculates, using a first reaction rate equation, a first bed layer reaction amount that contributes to an equilibrium state with the first inflow gas in one of two types of distributed bed layers among the bed layers distributed according to reaction rate-dependent conditions of the solid material and the liquid material contained in the bed layer, and obtains a first bed layer reaction amount that contributes to the equilibrium reaction of the first bed layer; a second bed layer reaction amount calculation unit that calculates a second bed layer reaction amount that contributes to an equilibrium state with the second inflow gas in the other second bed layer of the two types of distributed bed layers using a second reaction rate equation different from the first reaction rate equation, and obtains a second bed layer reaction amount that contributes to the equilibrium reaction in the second bed layer; a first bed layer equilibrium reaction calculation unit that calculates an equilibrium reaction between the first bed layer reaction component and the first inflow gas; a second bed layer equilibrium reaction calculation unit that calculates an equilibrium reaction between the second bed layer reaction component and the second inflow gas; a gas region equilibrium reaction calculation unit that calculates an equilibrium reaction between the first gas phase and a first substance including at least one of a solid substance and a liquid substance present in the gas region; Equipped with.

[0010] One aspect of the method for calculating a reaction in a furnace according to the present invention is a method for calculating a reaction in a furnace, in which a raw material ore is supplied from one end of a reactor and moved toward the other end while the raw material ore is brought into contact with a combustion gas supplied from the other end, thereby drying and reducing the raw material ore, a gas mixing amount calculation step of calculating a flow rate of an inflow gas that moves from the gas region in the reactor to a bed layer containing the raw ore, among a first gas phase containing the combustion gas flowing through a gas region in the reactor; an inflow gas distribution step of dividing the inflow gas into a first inflow gas and a second inflow gas; a first bed layer reaction amount calculation step of calculating, by using a first reaction rate equation, a first bed layer reaction amount contributing to an equilibrium state with the first inflow gas in one of two types of distributed bed layers distributed among the bed layers according to reaction rate-dependent conditions of the solid and liquid substances contained in the bed layers, to obtain a first bed layer reaction amount contributing to the equilibrium reaction in the first bed layer; a second bed layer reaction amount calculation step of calculating a second bed layer reaction amount contributing to an equilibrium state with the second inlet gas in the other of the two types of distributed bed layers using a second reaction rate equation different from the first reaction rate equation to obtain a second bed layer reaction amount contributing to the equilibrium reaction in the second bed layer; a first bed layer equilibrium reaction calculation step of calculating an equilibrium reaction between the first bed layer reaction component and the first inflow gas; a second bed layer equilibrium reaction calculation step of calculating an equilibrium reaction between the second bed layer reaction component and the second inflow gas; a gas region equilibrium reaction calculation step of calculating an equilibrium reaction between the first gas phase and a first substance including at least one of a solid substance and a liquid substance present in the gas region; Includes. [Effects of the Invention]

[0011] One aspect of the reactor reaction calculation device according to the present invention can calculate the reaction in the reactor with high accuracy. [Brief explanation of the drawings]

[0012] [Figure 1] 1 shows a schematic configuration of a rotary kiln to which a calculation device for furnace reactions according to an embodiment of the present invention is applied. [Figure 2] 1 is a block diagram showing the functions of a reactor reaction calculation device according to an embodiment of the present invention. [Figure 3] FIG. 1 is an explanatory diagram showing the state of segregation of bed layers in a rotary kiln. [Figure 4] FIG. 2 is an explanatory diagram showing the state of segregation of bed layers in the longitudinal direction in a rotary kiln. [Figure 5] 1 is a flowchart illustrating a method for calculating a reactor reaction according to an embodiment of the present invention. [Figure 6] 6 is a flowchart showing the operation of the gas region equilibrium reaction implementation step (step S19) in FIG. 5. [Figure 7] 6 is a flowchart showing the operation of the first bed layer equilibrium reaction implementation step (step S20) of FIG. 5. [Figure 8] 4 is a flowchart showing the operation of the second bed layer equilibrium reaction implementation step (step S21) of FIG. 3. [Figure 9] FIG. 1 is a hardware configuration diagram of a reactor reaction calculation device. [Figure 10] FIG. 2 is a diagram showing the flow of combustion gas and raw ore in a rotary kiln. [Figure 11] 1 is a flowchart illustrating a case where a method for calculating furnace reactions according to an embodiment of the present invention is applied to the entire rotary kiln. DETAILED DESCRIPTION OF THE INVENTION

[0013] Hereinafter, embodiments of the present invention will be described in detail. To facilitate understanding of the description, the same components in the drawings will be denoted by the same reference numerals, and duplicated descriptions will be omitted. Also, the scale of each component in the drawings may differ from the actual scale.

[0014] Before describing the calculation device for furnace reactions according to the embodiment of the present invention, the configuration of a rotary kiln to which the calculation device for furnace reactions according to the present embodiment is applied will be described.

[0015] <Rotary kiln> Fig. 1 shows the schematic configuration of a rotary kiln to which the calculation device for furnace reactions according to this embodiment is applied. As shown in Fig. 1, the rotary kiln 1 has a rotatable, approximately cylindrical kiln body 11 and a combustion material supply pipe 12 provided midway through the kiln body 11.

[0016] The kiln body 11 is a kiln having a cylindrical hollow structure, and is made of carbon steel with a thickness of 15 to 30 mm. The kiln body 11 is preferably provided with a refractory material on the inner peripheral wall surface thereof to enhance heat resistance.

[0017] The size of the kiln body 11 is preferably, for example, an inner diameter of 4.5 m to 5.5 m and a length in the major axis direction (total length) of 100 m to 110 m.

[0018] The kiln body 11 has an open end 11a at one end (left side in FIG. 1) that is inserted into and closed by the rotary kiln charging end (hereinafter also referred to as the "charging end") 14A, and an open end 11b at the other end (right side in FIG. 1) that is inserted into and closed by the rotary kiln discharge end (hereinafter also referred to as the "discharge end") 14B. The kiln body 11 is disposed at a slight incline from the charging end 14A toward the discharge end 14B, and is supported so as to be rotatable about its axis.

[0019] A raw material supply pipe 15 is installed through the charging end 14A to introduce raw material ore into the kiln body 11. A burner 16 is installed at the discharge end 14B and is introduced into the kiln body 11 through the open end 11b.

[0020] The raw ore may be nickel oxide ore (nickel oxide ore). The raw ore may be, for example, dried ore obtained by pre-drying nickel oxide ore with a dryer (rotary dryer) to remove some of the adhering moisture. The moisture content in the dried ore is about 15% to 25% by mass.

[0021] The nickel oxide ore used as the raw material ore is not particularly limited, but garnierite ore is preferably used in the smelting of ferronickel, an alloy primarily composed of iron and nickel. A typical composition of garnierite ore, calculated on a dried ore basis, is 2.1 mass% to 2.5 mass% Ni, 11 mass% to 23 mass% Fe, 20 mass% to 28 mass% MgO, 29 mass% to 39 mass% SiO2, and less than 0.5 mass% CaO, and the loss on ignition is 10 mass% to 15 mass%.

[0022] The burner 16 may be a pulverized coal-only burner or a pulverized coal and heavy oil mixed-fuel burner. The burner 16 burns pulverized coal or a fuel containing pulverized coal and heavy oil, etc., to generate combustion heat within the rotary kiln 1.

[0023] The combustion material supply pipe 12 is provided midway along the outer circumferential surface of the kiln body 11 and can supply combustion material into the kiln body 11. The combustion material includes at least one of substances that mainly enter the gas phase, such as volatile matter, and solid substances that mainly enter the bed layer, such as fixed carbon, and can be, for example, a carbonaceous material such as coal.

[0024] The volatile components include volatile substances such as hydrocarbon compounds, sulfur, and halogens.

[0025] Fixed carbon is the combustion fraction of char particles (mainly fixed carbon and ash), which are the residue after pyrolysis after moisture and volatile matter have been removed from coal, and is mainly composed of carbon, excluding ash.

[0026] In Figure 1, only one combustion material supply pipe 12 is provided on the outer peripheral surface of the kiln body 11, but multiple pipes may be provided on the outer peripheral surface of the kiln body 11 along the axial direction or around the axis of the kiln body 11.

[0027] Combustion materials are often not a single type of carbonaceous material, but a mixture of several different types of carbonaceous material, and the particle size of the carbonaceous material added often has a distribution. Note that particle size refers to the volume average particle size based on the effective diameter, and particle size is measured by, for example, laser diffraction / scattering method, dynamic light scattering method, or classification method. When using laser diffraction / scattering method, the particle size (D) is the particle size at which the cumulative particle size distribution from the small particle size side reaches 50% in the volume-based particle size distribution measured by laser diffraction / scattering method. 50 ) can be used as the average particle size.

[0028] Raw ore is charged into the kiln body 11 through raw material supply pipe 15 installed at charging end 14A, and combustion material is introduced into the kiln body 11 through combustion material supply pipe 12. High-temperature combustion gas generated by burning pulverized coal, heavy oil, etc. using burner 16 installed at discharge end 14B is blown from discharge end 14B toward charging end 14A, i.e., in the opposite direction to the flow of the raw ore.

[0029] In the kiln body 11, raw ore is charged through the charging end 14A. As the kiln body 11 rotates at a predetermined speed, the raw ore is transported from the opening end 11a of the kiln body 11 through the raw material supply pipe 15 from one end (open end 11a) toward the other end (discharge end 14B). As the raw ore moves through the kiln body 11, it comes into countercurrent contact with combustion gas flowing from the discharge end 14B toward the charging end 14A. The raw ore is heated by the combustion heat and flame of the high-temperature combustion gas generated by burning fuel such as pulverized coal or heavy oil in the burner 16. Furthermore, combustion material introduced into the kiln body 11 through the combustion material supply pipe 12 is combusted by the combustion gas in the kiln body 11. The raw ore is also heated by the combustion heat generated by the combustion of the combustion material introduced through the combustion material supply pipe 12. Therefore, as the kiln body 11 rotates, the raw ore moves from the charging end 14A to the discharge end 14B of the kiln body 11, and is heated by the combustion heat and flame of the combustion gas generated by the combustion of fuel in the burner 16 and the combustion heat generated by the combustion of the combustion material, gradually increasing its temperature.

[0030] Within the kiln body 11, material transfer occurs between the raw ore and combustion gas due to the evaporation of moisture contained in the raw ore and combustion materials, the volatilization and condensation of volatile matter contained in the combustion materials, and the scattering and falling of ash contained in the burner fuel and combustion materials. The moisture and volatile matter produced by the thermal decomposition of the combustion materials supplied midway through the rotary kiln 1 move toward the charging end 14A together with the combustion gas, and char particles move toward the discharge end 14B together with the raw ore.

[0031] By the time the raw ore in the kiln body 11 reaches the discharge end 14B, the moisture contained in the raw ore is almost completely removed, the raw ore is burned, and the raw ore is partially reduced to become cinder, which is then discharged from the discharge end 14B.

[0032] The calcined ore has a temperature of, for example, 800 to 900°C and a particle size of about 10 mm to 100 mm.

[0033] A grate (sieving device) 17 is provided at the discharge opening of the discharge end 14B to separate the cinders with particle diameters of about 10 mm to 100 mm from the sintered lumps (particle diameters of about 100 mm to 500 mm) generated in the rotary kiln 1. The grate 17 is made of, for example, an iron grate with an opening of about 100 mm. The cinders discharged from the discharge end 14B pass through the grate 17 and then pass through a cinder discharge chute 18 to be transported to the next process.

[0034] <Calculation device for reactor reactions> Next, the calculation device for furnace reactions according to this embodiment will be described. FIG. 2 is a block diagram showing the functions of the calculation device for furnace reactions according to this embodiment. In FIG. 2, the calculation device for furnace reactions will be described as a unit operation model in region A, where combustion material is supplied from combustion material supply pipe 12, when the rotary kiln 1 is divided into multiple regions. In FIG. 2, the adjacent region on the side where combustion gas is blown in from region A is referred to as region (A+1), and the adjacent region on the side where raw material ore is charged is referred to as region (A-1).

[0035] In the following description, the gas region means a region in the rotary kiln 1 where a gas phase containing gas and dust flows, and the bed layer (bed region) means a region where the raw ore moves.

[0036] The gas phases in region A are referred to as a first gas phase G1, a second gas phase G2, a third gas phase G3, or a third mixed gas phase G3'. Each gas phase in region A is defined as follows: The "first gas phase G1" is the combustion gas that flows into region A from one of the adjacent regions (region (A+1)). The combustion gas includes gases that contribute to the equilibrium reaction, such as volatile matter, oxygen, carbon dioxide, carbon monoxide, and hydrogen, and may also include inert substances that do not contribute to the reaction (e.g., nitrogen, etc.). The "second gas phase G2" is a gas phase obtained by combining the added gas phase distributed from the combustion material and the first gas phase. The "third gas phase G3" is a gas phase generated in the mixed gas phase calculation unit 215, which will be described later. The "third mixed gas phase G3'" is a gas phase generated in the gas amount summing unit 222, which will be described later.

[0037] The bed layer fed into the rotary kiln 1 is described as bed layer S1. The "bed layer S1" refers to raw ore supplied from the charging end 14A inside the kiln body 11. The raw ore includes solid materials, such as ore compounds such as nickel oxide, iron oxide, and magnesium oxide, and fixed carbon, as well as liquid materials, and may also include inactive materials that do not contribute to the reactions in the first bed layer equilibrium reaction calculation unit 214-1 and the second bed layer equilibrium reaction calculation unit 214-2, which will be described later.

[0038] The bed layers present in Area A are referred to as the first bed layer S11, the preliminary corrected first bed layer S11B, the corrected first bed layer S11C, the first mixed bed layer S12, the second bed layer S21, the corrected second bed layer S21A, or the second mixed bed layer S22. Each bed layer in Area A is defined as follows: The "first bed layer S11" refers to a bed layer containing solid and liquid substances with a fast reaction rate, which are distributed according to the reaction rate-causing conditions of the solid and liquid substances contained in the bed layer S1, among the bed layers S1 flowing into the region A from one adjacent region (region (A-1)). The first bed layer S11 may contain inactive substances that do not contribute to the reaction in the first bed layer equilibrium reaction calculation unit 214-1 described later, similar to the bed layer S1. The reaction rate-attributing conditions are set appropriately depending on the properties of the solid and liquid substances contained in the bed layer. The properties of the solid substances include the particle size, specific gravity, type, segregation, etc. of the solid substance. The properties of the liquid substances include the specific gravity, type, segregation, etc. of the liquid substance. The "first bed layer S11A" is a bed layer obtained by combining the additive bed layer distributed from the combustion material and the first bed layer S11. The "preliminarily corrected first bed layer S11B" is a first bed layer in which the quantity of material in the first bed layer S11 has been corrected taking into consideration the amount of a fourth substance contained in the first bed layer S11, which contains at least one of a solid substance or a liquid substance, moving into the gas region, and the amount of a fifth substance contained in the first mixed substance present in the gas region and which contains at least one of a solid substance or a liquid substance, and which is generated in the first mixed substance calculation unit described below, moving into the first bed layer S11. The "corrected first bed layer S11C" is a first bed layer in which the quantity of material in the preliminary corrected first bed layer S11B has been corrected, taking into account the amount of a second substance contained in the preliminary corrected first bed layer S11B, which contains at least one of a solid substance or a liquid substance, that has moved to the second bed layer S21, and the amount of a third substance contained in the second bed layer S21, which contains at least one of a solid substance or a liquid substance, that has moved to the preliminary corrected first bed layer S11B. The "first mixed bed layer S12" is a bed layer generated in the first mixed bed layer calculation unit 218, which will be described later. The "second bed layer S21" refers to a bed layer other than the first bed layer S11 among the bed layers S1 flowing into the region A from one adjacent region (region (A-1)), and including solid and liquid substances having a slower reaction rate than the first bed layer S11. The second bed layer S21, like the bed layer S1, may include inactive substances that do not contribute to the reaction in the second bed layer equilibrium reaction calculation unit 214-2 described later. The "corrected second bed layer S21A" is a second bed layer in which the quantity of material in the second bed layer S21 has been corrected taking into account the amount of a second substance contained in the preliminary corrected first bed layer S11B, which contains at least one of a solid substance or a liquid substance, that has moved to the second bed layer S21, and the amount of a third substance contained in the second bed layer S21, which contains at least one of a solid substance or a liquid substance, that has moved to the preliminary corrected first bed layer S11B. The "second mixed bed layer S22" is a bed layer generated in the second mixed bed layer calculation unit 220, which will be described later.

[0039] The gases flowing into the bed from the gas region in region A are referred to as inlet gas G11, first inlet gas G111, second inlet gas G112, first mixed inlet gas G111A, second mixed inlet gas G112A, or mixed inlet gas G12. In the following description, the inlet gases in region A are defined as follows: "Inlet gas G11" is the gas that flows into the bed layer from the gas region in region A. The "first inflow gas G111" is the inflow gas that flows into the first bed layer S11 out of the inflow gas G11 that has flowed into the bed layer. The "second inflow gas G112" is the inflow gas that flows into the second bed layer S21 out of the inflow gas G11 that has flowed into the bed layer. The "first mixed inflow gas G111A" is a gas obtained by mixing the unreacted portion of the first bed layer and the product portion of the first inflow gas generated in the first bed layer equilibrium reaction calculation section 214-1. The "second mixed inflow gas G112A" is a gas obtained by mixing the corrected second bed layer unreacted portion and the second inflow gas product portion generated in the second bed layer equilibrium reaction calculation section 214-2. The "mixed inflow gas G12" is a gas obtained by combining the first mixed inflow gas G111A and the second mixed inflow gas G112A.

[0040] The substances present in region A are referred to as first substance M1, second substance M2, third substance M3, fourth substance M4, fifth substance M5, first mixed substance M11 or first modified mixed substance M12. In the following description, each substance in region A is defined as follows: The "first substance M1" is a substance that flows from the gas region of region (A+1) into the gas region of region A, and includes at least one of a solid substance and a liquid substance present in the gas region of region A. The "second substance M2", "third substance M3" and "fourth substance M4" include at least one of a solid substance and a liquid substance present in the bed layer. The "fifth substance M5" includes at least one of a solid substance and a liquid substance that exists in the gas region. The "first mixed material M11" is a mixture of the first material M1 generated in the first mixed material calculation unit 216, which will be described later. The "first modified mixed material M12" is a mixture of the first material M1 produced in the first mixed material quantity modifying unit 223, which will be described later.

[0041] The reactants, unreacted components, and products calculated in region A are referred to as gas reactant GR1, gas unreacted component Gr1, gas product GP1, first substance reactant MR1, first substance unreacted component Mr1, first substance product MP1, first inflow gas reactant GR111, first inflow gas unreacted component Gr111, first inflow gas product GP111, second inflow gas reactant GR112, second inflow gas unreacted component Gr112, second inflow gas product GP112, first bed layer reactant SR1, first bed layer unreacted component Sr1, first bed layer product SP1, second bed layer reactant SR2, second bed layer unreacted component Sr2, or second bed layer product SP2. These are calculated by any of the components constituting the reactor reaction calculation device shown in Figure 2.

[0042] The configuration of the calculation device for in-furnace reactions shown in Fig. 2 will be described. As shown in Fig. 2, the calculation device for in-furnace reactions 20 includes a distribution unit (distribution unit) 201 for combustion materials, a gas mixture amount calculation unit 202, a first bed layer preliminary material amount correction unit 203, a first bed layer material amount correction unit 204, a second bed layer material amount correction unit 205, an inflow gas distribution unit 206, a gas reaction amount calculation unit 207, a first material reaction amount calculation unit 208, a first bed layer reaction amount calculation unit 209, a first inflow gas reaction amount calculation unit 210, a second bed layer reaction amount calculation unit 211, a second inflow gas The system includes a gas reaction amount calculation unit 212, a gas region equilibrium reaction calculation unit 213, a first bed layer equilibrium reaction calculation unit 214-1, a second bed layer equilibrium reaction calculation unit 214-2, a mixed gas phase calculation unit 215, a first mixed substance calculation unit 216, a first mixed inflow gas calculation unit 217, a first mixed bed layer calculation unit 218, a second mixed inflow gas calculation unit 219, a second mixed bed layer calculation unit 220, an inflow gas summing unit 221, a gas amount summing unit 222, and a first mixed substance amount correction unit 223.

[0043] The reactor reaction calculation device 20 may further include a bed layer distribution unit 224 that distributes the bed layer S1 into two types of distributed bed layers (a first bed layer S11 and a second bed layer S21) in advance according to reaction rate-dependent conditions of the solid and liquid materials contained in the bed layer S1. In this embodiment, the reactor reaction calculation device 20 supplies the first bed layer S11 generated in the bed layer distribution unit 224 to the first bed layer preliminary quantity correction unit 203, and supplies the second bed layer S21 generated in the bed layer distribution unit 224 to the second bed layer quantity correction unit 205.

[0044] When raw ores with different particle sizes, densities, compositions, and types of ores are simultaneously fed into the rotary kiln 1, striped segregation, in which the first bed layer S11 and the second bed layer S21 are segregated into stripes, occurs in the bed layer, as shown in Figure 3. When striped segregation occurs, the first bed layer S11 and the second bed layer S21 periodically exhibit uneven composition distribution along the longitudinal direction of the rotary kiln 1, as shown in Figure 4. This may result in significantly different reaction patterns for each stripe. When raw ores with different particle sizes, densities, and other characteristics are fed into the rotary kiln 1, treating the bed layer in the rotary kiln 1 as a single concentration cannot accurately represent the composition distribution within the raw ore due to the different reaction rates of the raw ores within the bed layer. Therefore, even if a simulation and analysis of the furnace reaction are performed, it is difficult to fit the simulation to the actual furnace reaction. Therefore, when considering the phenomena occurring inside the furnace, it is necessary to consider not only the difference in reaction rate of the raw ore but also the composition distribution in the bed layer caused by segregation.

[0045] The furnace reaction calculation device 20 divides the rotary kiln 1 into a first bed layer S11 and a second bed layer S21 and performs equilibrium reaction calculations for each bed layer. That is, when simulating furnace reactions based on the Gibbs energy minimization method, the furnace reaction calculation device 20 takes into account the different reaction rates of the raw ore in the rotary kiln 1 and the composition distribution (unevenness) within the raw ore. The furnace reaction calculation device 20 then calculates the reaction amounts required for each equilibrium reaction using different reaction rate equations independently for the first bed layer S11, which has a fast reaction rate, and the second bed layer S21, which has a slow reaction rate. As a result, the furnace reaction calculation device 20 can accurately calculate equilibrium reactions based on the Gibbs energy minimization method, taking into account the difference in reaction rate of the raw ore contained in the first bed layer S11 and the first bed layer S12 and the composition distribution within the raw ore. Therefore, the furnace reaction calculation device 20 can analyze the furnace reaction in the bed layer with high accuracy. The furnace reaction calculation device 20 can perform fitting with higher accuracy, so it can be effectively used to analyze operational data. For example, since it can calculate two different states occurring in the rotary kiln 1 simultaneously at the same position, it can also represent a state in which opposing reactions, such as an oxidation reaction caused by an oxidizing gas and a reduction reaction caused by a reducing substance in the bed layer, occur at the same position in the longitudinal direction of the rotary kiln 1, which was difficult with conventional methods. Therefore, the furnace reaction calculation device 20 can perform fitting that is closer to the actual furnace state.

[0046] The furnace reaction calculation device 20 also calculates the equilibrium reactions between the gas flowing in the gas region and the solid or liquid substances present in the gas region based on the Gibbs energy minimization method, thereby enabling calculation of the equilibrium reactions in the gas region of the rotary kiln 1.

[0047] Therefore, by taking into consideration that the bed layer in the rotary kiln 1 exhibits two different types of reaction behavior depending on the particle size of the raw ore, etc., the calculation device 20 for in-furnace reactions can analyze the reactions occurring in the rotary kiln 1 with high precision and improve the calculation accuracy of the reactions. Therefore, the calculation device 20 for in-furnace reactions can calculate the reactions in the rotary kiln 1 with higher precision.

[0048] 2, the reactor reaction calculation device 20 inputs the quantities, temperatures, pressures, etc. of the first bed layer S11 and the second bed layer S21 flowing into the bed layers of the region A and the first gas phase G1 flowing into the gas region. The reactor reaction calculation device 20 uses these input values ​​to perform mass transfer and reaction calculations in each divided region A.

[0049] The first bed layer S11 and the second bed layer S21 are determined based on the amount of material when the bed layer S1 supplied from the charging end 14A in the kiln body 11 is distributed in the area including the charging end 14A and moved to area A.

[0050] The reactor reaction calculation device 20 defines the amount of gas (inlet gas G11) mixed into the bed layer in the first gas phase G1 flowing in the gas region when calculating mass transfer and reactions in the region A. As a method for defining this mixed gas amount, for example, there is a method in which the ratio of the amount mixed into the bed layer to the total amount of gas is given as a fixed value.

[0051] Furthermore, the reactor reaction calculation device 20 defines the amount of mass transfer to the gas region due to scattering or volatilization of solid or liquid materials moving in the bed layer, and the amount of mass transfer to the bed layer due to falling (sedimentation) or coagulation of solid or liquid materials moving in the gas region. For example, the amount of mass transfer can be defined for each calculation substance based on the settling velocity calculated from the specific gravity, particle size, etc. of the solid or liquid material.

[0052] As shown in FIG. 2, the distribution unit 201 distributes the combustion material 230 to the gas phase and the bed layer contained in the combustion material 230 at a mass flow rate by setting the type, components, etc. of the combustion material 230 in advance.

[0053] The distributed gas phase and bed layer are supplied into the region A as an added gas phase AG11 and an added bed layer AS11.

[0054] The gas mixture amount calculation unit 202 calculates the mass flow rate of the inflow gas G11 moving from the gas region to the bed layer, out of the second gas phase G2 which includes the first gas phase G1, which is the combustion gas flowing from one adjacent region (region (A+1)) into the gas region of region A in the furnace, and the added gas phase AG11 of the combustion material distributed by the distribution unit 201.

[0055] The first bed layer preliminary quantity correcting unit 203 corrects the quantity of the first bed layer S11A obtained by mixing the first bed layer S11 and the additive bed layer AS11, and calculates the preliminary corrected first bed layer S11B.

[0056] The amount of material in the first bed layer S11A is corrected taking into account the amount of movement of a fourth substance M4, which is contained in the first bed layer S11A and includes at least one of a solid substance or a liquid substance, to the first mixed substance M11 present in the gas region, and the amount of movement of a fifth substance M5, which is contained in the first mixed substance M11 present in the gas region and includes at least one of a solid substance or a liquid substance, to the first bed layer S11A.

[0057] The first bed layer quantity correction unit 204 corrects the quantity of the first bed layer S11 by subtracting the amount of second substance M2, which includes at least one of a solid substance and a liquid substance and is contained in the preliminary corrected first bed layer S11B, from the amount of second substance M2 transferred to the second bed layer S21, to calculate the corrected first bed layer S11C.

[0058] In this embodiment, the first bed layer quantity correcting unit 204 corrects the quantity of the first bed layer S11 in consideration of the amount of the second substance M2 transferred to the second bed layer S21 and the amount of the third substance M3, which is contained in the second bed layer S21 and includes at least one of a solid substance and a liquid substance, transferred to the preliminary corrected first bed layer S11B. As a result, the corrected first bed layer S11C is calculated.

[0059] The second bed layer quantity correcting unit 205 corrects the quantity of the second bed layer S21 by subtracting the amount of the third substance M3 contained in the second bed layer S21 that has moved to the preliminary corrected first bed layer S11B, and calculates the corrected second bed layer S21A.

[0060] In this embodiment, the second bed layer quantity correcting unit 205 corrects the quantity of the first bed layer S11 in consideration of the amount of the second substance M2 contained in the preliminary corrected first bed layer S11B moving to the second bed layer S21 and the amount of the third substance M3 contained in the second bed layer S21 moving to the preliminary corrected first bed layer S11B. This allows the corrected first bed layer S11C to be calculated.

[0061] The inflow gas distribution unit 206 has the function of distributing the inflow gas G11 into a first inflow gas G111 and a second inflow gas G112 by mass flow rate by setting in advance the type and components of the inflow gas G11 calculated by the gas mixture amount calculation unit 202.

[0062] The gas reaction amount calculation unit 207 calculates the amount of gas reaction that contributes to the equilibrium state between the second gas phase G2 and the first substance M1. That is, the gas reaction amount calculation unit 207 calculates the amount of gas reaction that is expected to occur when the second gas phase G2 and the first substance M1 reach equilibrium during the residence time of the second gas phase G2 in the gas region in region A. The reaction amount during the residence time of the second gas phase G2 in region A is calculated using reaction rate parameters set for each substance and the residence time in region A. The gas reaction amount can be calculated using a reaction rate equation to determine the amount of reaction during the residence time in region A. Then, the gas reaction amount calculation unit 207 calculates the amount of reaction that contributes to the equilibrium state between the second gas phase G2 and the first substance M1 from the gas reaction amount as the gas reaction amount GR1.

[0063] In this specification, the reactive component refers to the flow rate that contributes to a reaction by contact with the second gas phase G2 or the first substance during the residence time of the first substance M1 present in the second gas phase G2 or gas region in region A. Hereinafter, other reactive components also refer to the flow rate that contributes to a reaction by contacting any of the gas phase, solid phase, and substance with any of the gas phase, solid phase, and substance.

[0064] That is, the gas reaction portion GR1 is a reaction portion that contributes to an equilibrium reaction with the first substance M1 present in the gas region of the second gas phase G2. The first substance reaction portion MR1 is a reaction portion that contributes to an equilibrium state with the second gas phase G2 of the first substance M1 present in the gas region flowing in from one adjacent region (region (A+1)).

[0065] The gas reaction amount calculation unit 207 can determine the amount of the second gas phase G2 and the solid or liquid substance present in the first substance M1 that reaches an equilibrium state by employing, for example, a model that can calculate the amount from the reaction rate according to the Arrhenius reaction rate equation and the residence time of the second gas phase G2 in the region A. Then, the gas reaction amount calculation unit 207 divides the amount into a gas reaction amount GR1 corresponding to the gas reaction amount and an unreacted gas amount Gr1 corresponding to the remaining unreacted mass flow rate.

[0066] In the rotary kiln 1, the second gas phase G2 and the first substance M1 present in the gas region flow in parallel. Because the second gas phase G2 has a high flow rate and a short residence time in region A, the second gas phase G2 and the first substance M1 present in the gas region do not reach equilibrium. In this embodiment, the gas reaction amount calculation unit 207 assumes that only a portion of the second gas phase G2 in a certain region (region A) within the rotary kiln 1 reaches equilibrium with the first substance M1 present in the gas region, and divides the second gas phase G2 into a gas reaction portion GR1 and an unreacted gas portion Gr1 other than the gas reaction portion GR1 according to the corresponding gas reaction amount.

[0067] The first material reaction amount calculation unit 208 calculates the first material reaction amount that contributes to the equilibrium state of the first material M1 with the second gas phase G2. That is, the first material reaction amount calculation unit 208 calculates the reaction amount that is expected to occur when the first material M1 reacts in equilibrium with the second gas phase G2 during the residence time of the first material M1 in the gas region in region A as the first material reaction amount. The first material reaction amount calculation unit 208 can be calculated using a reaction rate equation. Then, from the first material reaction amount, the first material reaction amount calculation unit 208 calculates the first material reaction amount MR1 of the first material M1 that contributes to the equilibrium state with the second gas phase G2.

[0068] Like the gas reaction amount calculation unit 207, the first substance reaction amount calculation unit 208 can determine the amount of the second gas phase G2 and the first substance M1 present in the gas region that reaches an equilibrium state by employing, for example, a model that can calculate the amount from the reaction rate according to the Arrhenius reaction rate equation and the residence time of the first substance M1 in the region A. Then, the first substance reaction amount calculation unit 208 divides the amount into a first substance reaction amount MR1 corresponding to the first substance reaction amount and an unreacted first substance amount Mr1 corresponding to the remaining unreacted mass flow rate.

[0069] As described above, in the rotary kiln 1, the second gas phase G2 and the first substance M1 flow in parallel, and the second gas phase G2 has a high flow rate and a short residence time in region A, so the second gas phase G2 and the first substance M1 do not reach equilibrium. In this embodiment, the first substance reaction amount calculation unit 208 assumes that only a portion of the second gas phase G2 in a certain region (region A) in the rotary kiln 1 reaches equilibrium with the first substance M1, and divides the first substance M1 into a reacted portion MR1 of the first substance and an unreacted portion Mr1 of the first substance other than the reacted portion MR1 of the first substance according to the corresponding first substance reaction amount.

[0070] The first bed layer reaction amount calculation unit 209 calculates the amount of the first bed layer reaction that contributes to the equilibrium state of the modified first bed layer S11C with the first inflow gas G111 using the first reaction rate equation. That is, the first bed layer reaction amount calculation unit 209 calculates the amount of the first bed layer reaction that is expected to occur in equilibrium with the first inflow gas G111 in the modified first bed layer S11C during the movement time of the modified first bed layer S11C in the region A using the first reaction rate equation.

[0071] The first bed layer reaction amount calculation unit 209 can calculate the first bed layer reaction amount by using, for example, reaction rate parameters set for each substance included in the modified first bed layer S11C in region A, residence time in region A, etc. in the first reaction rate equation. Then, the first bed layer reaction amount calculation unit 209 calculates the first bed layer reaction amount SR1 that contributes to the equilibrium reaction with the first inflow gas G111 in the modified first bed layer S11C from the first bed layer reaction amount.

[0072] That is, the first bed layer reaction portion SR1 is a reaction portion that contributes to the equilibrium reaction with the first inflow gas G111 in the modified first bed layer S11C.

[0073] The first bed layer reaction amount calculation unit 209 can calculate the first bed layer reaction amount by, for example, employing a model that can calculate the reaction rate according to the Arrhenius reaction rate equation and the transit time of the modified first bed layer S11C within region A. The first bed layer reaction amount calculation unit 209 calculates the first bed layer reaction amount of the modified first bed layer S11C that contributes to the equilibrium reaction with the first inflow gas G111 only for the amount estimated by the Arrhenius reaction rate equation that the modified first bed layer S11C and the first inflow gas G111 have reached equilibrium based on the residence time of the modified first bed layer S11C. The first bed layer reaction amount calculation unit 209 then divides the first bed layer reaction amount into a first bed layer reaction amount SR1 corresponding to the first bed layer reaction amount and a first bed layer unreacted amount Sr1 corresponding to the remaining unreacted mass flow rate.

[0074] In the rotary kiln 1, the modified first bed layer S11C and the first inflow gas G111 flow countercurrently. Because the first inflow gas G111 has a high flow velocity and a short residence time in region A, the modified first bed layer S11C and the first inflow gas G111 do not reach equilibrium. In this embodiment, the first bed layer reaction amount calculation unit 209 assumes that only a portion of the modified first bed layer S11C in a certain region (region A) in the rotary kiln 1 reaches equilibrium with the first inflow gas G111, and divides the first bed layer reacted portion SR1 into the first bed layer reacted portion SR1 and the first bed layer unreacted portion Sr1 other than the first bed layer reacted portion according to the corresponding first bed layer reaction amount.

[0075] The first inlet gas reaction amount calculation unit 210 calculates the first inlet gas reaction amount that contributes to the equilibrium state between the first inlet gas G111 and the modified first bed layer S11C. That is, the first inlet gas reaction amount calculation unit 210 calculates the first inlet gas reaction amount that is expected to occur when the first inlet gas G111 is in equilibrium with the modified first bed layer S11C during the residence time of the first inlet gas G111 in the bed layer in region A. The first inlet reaction amount during the residence time of the first inlet gas G111 in region A can be calculated using reaction rate parameters set for each substance and the residence time in region A, for example, using a reaction rate equation. Then, the first inlet gas reaction amount calculation unit 210 calculates the first inlet gas reaction amount GR111 that contributes to the equilibrium reaction with the modified first bed layer S11C from the first inlet gas reaction amount.

[0076] That is, the first inflow gas reaction portion GR111 is the reaction portion of the first inflow gas G111 that contributes to the equilibrium state of the first modified bed layer S11B.

[0077] The first inflow gas reaction amount calculation unit 210 can calculate the reaction amount at which the first inflow gas G111 reaches equilibrium with the modified first bed layer S11C, for example, using a model similar to the model used in the first bed layer reaction amount calculation unit 209.

[0078] Then, the first inflow gas reaction amount calculation unit 210 divides it into a first inflow gas reaction amount GR111 corresponding to the first inflow gas reaction amount and a first inflow gas unreacted amount Gr111 corresponding to the unreacted mass flow rate other than the remaining first inflow gas reaction amount GR111.

[0079] As described above, in the rotary kiln 1, the first inflow gas G111 and the modified first bed layer S11C flow countercurrently, and the first inflow gas G111 has a high flow rate and a short residence time in region A, so the first inflow gas G111 and the modified first bed layer S11C do not reach equilibrium. In this embodiment, the first inflow gas reaction amount calculation unit 210 assumes that only a portion of the first inflow gas G111 in a certain region (region A) in the rotary kiln 1 reaches equilibrium with the modified first bed layer S11C, and divides the first inflow gas G111 into a reacted portion of the first inflow gas GR111 and an unreacted portion of the first inflow gas Gr111 according to the corresponding gas reaction amount.

[0080] The second bed layer reaction amount calculation unit 211 calculates the amount of bed reaction that contributes to the equilibrium state of the modified second bed layer S21A with the second inflow gas G112 using the second reaction rate equation. That is, the second bed layer reaction amount calculation unit 211 calculates the amount of second bed layer reaction that is expected to occur in equilibrium with the second inflow gas G112 in the modified second bed layer S21A during the movement time of the modified second bed layer S21A in the region A (i.e., the movement time of each compound including solid substances and liquid substances present in the modified second bed layer S21A) using the second reaction rate equation.

[0081] The second reaction rate equation is a reaction rate equation different from the first reaction rate equation used in the first bed layer reaction amount calculation unit 209 .

[0082] The second bed layer reaction amount calculation unit 211 calculates the second bed layer reaction amount SR2 of the corrected second bed layer S21A that contributes to the equilibrium state with the second inflow gas G112 from the second bed layer reaction amount obtained using the second reaction rate equation.

[0083] The second bed layer reaction amount calculation unit 211 can calculate the reaction amount that reaches equilibrium between the modified second bed layer S21A and the second inflow gas GR112, similarly to the first bed layer reaction amount calculation unit 209, by employing, for example, a model that can be calculated from a reaction rate according to an Arrhenius reaction rate equation and the transit time of the modified second bed layer S21A within region A. The second bed layer reaction amount calculation unit 211 calculates the second bed reaction amount of the modified second bed layer S21A that contributes to the reaction with the second inflow gas GR112 only for the amount estimated by the Arrhenius reaction rate equation that the modified second bed layer S21A and the second inflow gas GR112 reach equilibrium based on the residence time in the modified second bed layer S21A. The second bed layer reaction amount calculation unit 211 then divides the amount into a second bed layer reaction amount SR2 corresponding to the second bed layer reaction amount and a second bed layer unreacted amount Sr2 corresponding to the remaining unreacted mass flow rate.

[0084] The second bed layer reaction portion SR2 is a reaction portion that contributes to the equilibrium reaction with the second inflow gas G112 in the modified second bed layer S2A.

[0085] In the rotary kiln 1, the modified second bed layer S21A and the second inflow gas G112 flow countercurrently. Because the second inflow gas G112 has a high flow rate and a short residence time in region A, the modified second bed layer S21A and the second inflow gas G112 do not reach equilibrium. In this embodiment, the second bed layer reaction amount calculation unit 211 assumes that only a portion of the modified second bed layer S2A in a certain region (region A) in the rotary kiln 1 reaches equilibrium with the second inflow gas G112, and divides the second bed layer reacted portion SR2 into the second bed layer reacted portion SR2 and the remaining second bed layer unreacted portion Sr2 according to the corresponding second bed layer reaction amount.

[0086] The second inlet gas reaction amount calculation unit 212 calculates the second inlet gas reaction amount that contributes to the equilibrium state of the second inlet gas G112 with the modified second bed layer S21A. That is, the second inlet gas reaction amount calculation unit 212 calculates the second inlet gas reaction amount that the second inlet gas G112 is expected to undergo in equilibrium with the modified second bed layer S21A during the residence time of the second inlet gas G112 in the bed layer in region A. The second inlet reaction amount during the residence time of the second inlet gas G112 in region A can be calculated using reaction rate parameters set for each substance and the residence time in region A, for example, from a reaction rate equation. Then, the second inlet gas reaction amount calculation unit 212 calculates the second inlet gas reaction amount GR112 that contributes to the equilibrium reaction with the second modified bed layer S21A of the second inlet gas G112.

[0087] That is, the second inflow gas reaction portion GR112 is the reaction portion of the second inflow gas G112 that contributes to the equilibrium state of the second amended bed layer S21A.

[0088] The second inflow gas reaction amount calculation unit 212 can calculate the reaction amount at which the second inflow gas G112 reaches equilibrium with the modified second bed layer S21A, for example, using a model similar to the model used in the second bed layer reaction amount calculation unit 211.

[0089] The second inflow gas reaction amount calculation unit 212 divides the second inflow gas into a second inflow gas reaction amount GR112 corresponding to the second inflow gas reaction amount and a second inflow gas unreacted amount Gr112 corresponding to the unreacted mass flow rate other than the remaining second inflow gas reaction amount.

[0090] As described above, in the rotary kiln 1, the second inflow gas G112 and the modified second bed layer S21A flow countercurrently, and the second inflow gas G112 has a high flow rate and a short residence time in region A, so the second inflow gas G112 and the modified second bed layer S21A do not reach equilibrium. In this embodiment, the second inflow gas reaction amount calculation unit 212 assumes that only a portion of the modified second bed layer S21A in a certain region (region A) in the rotary kiln 1 reaches equilibrium with the second inflow gas G112, and divides the second inflow gas G112 into the reacted portion of the second inflow gas GR112 and the remaining unreacted portion of the second inflow gas GR112 according to the corresponding second inflow gas reaction amount.

[0091] The gas region equilibrium reaction calculation unit 213 calculates the equilibrium reaction between the second gas phase G2 flowing in the gas region within region A and the first substance M1 present in the gas region. That is, the gas region equilibrium reaction calculation unit 213 performs equilibrium reaction calculation of the gas region using the gas reaction portion GR1 of the second gas phase G2 that contributes to the equilibrium reaction with the first substance M1 present in the gas region, and the first substance reaction portion MR1 of the first substance M1 that contributes to the equilibrium reaction with the gas reaction portion GR1, and calculates at least the compositions, amounts, etc. of the gas reactant GR1 and the first substance reactant MR1, and the changes in the heat quantities and flow rates of the gas reactant GR1 and the first substance reactant MR1 when the gas reactant GR1 and the first substance reactant MR1 reach an equilibrium state.

[0092] The gas region equilibrium reaction calculation unit 213 has a function of performing equilibrium reaction calculations to determine the type, phase, flow rate, etc. of the product substance that may be generated by the reaction between the gas reaction component GR1 and the first substance reaction component MR1 so that the free energy of the product substance is minimized by setting the type, phase, flow rate, etc. in advance. That is, the gas region equilibrium reaction calculation unit 213 can perform equilibrium reaction calculations based on the Gibbs energy minimization method. The gas region equilibrium reaction calculation unit 213 calculates and outputs the change in heat quantity, flow rate, type, composition, amount, phase, etc. of the product substance when the gas reaction component GR1 and the first substance reaction component MR1 reach equilibrium.

[0093] The gas region equilibrium reaction calculation unit 213 calculates the product gas generated by the reaction between the gas reaction amount GR1 and the first substance reaction amount MR1 and the unused amount of the gas reaction amount GR1 as the gas product amount GP1.

[0094] In addition, the gas region equilibrium reaction calculation unit 213 calculates the product material generated by the reaction between the gas reaction amount GR1 and the first substance reaction amount MR1 and the unused amount of the first substance reaction amount MR1 as the first substance product amount MP1.

[0095] That is, the gas product GP1 is a gas phase obtained by the reaction of the gas reaction GR1 with the first substance reaction MR1, and the unused portion of the gas reaction GR1. The first substance product MP1 is a gas phase obtained by the reaction of the gas reaction GR1 with the first substance reaction MR1, and the unused portion of the first substance reaction MR1.

[0096] The first bed layer equilibrium reaction calculation unit 214-1 calculates the equilibrium reaction between the modified first bed layer S11C moving through the bed layers in region A and the first inflow gas G111 present in the modified first bed layer S11C. That is, the first bed layer equilibrium reaction calculation unit 214-1 performs the equilibrium reaction calculation of the first bed layer using the first bed layer reaction portion SR1 that contributes to the equilibrium reaction with the first inflow gas G111 in the modified first bed layer S11C and the first inflow gas reaction portion GR111 that contributes to the equilibrium state with the modified first bed layer S11C. As a result, the first bed layer equilibrium reaction calculation unit 214-1 calculates at least the change in heat quantity and flow rate of each of the first bed layer reactant SR1 and the first inflow gas reactant GR111 when they reach equilibrium.

[0097] The first bed layer equilibrium reaction calculation unit 214-1 has a function of performing equilibrium reaction calculations to determine the type, phase, flow rate, etc. of the product (product) that may be produced by the reaction between the first bed layer reaction component SR1 and the first inflow gas reaction component GR111 so as to minimize the free energy of the product, by setting in advance the type, phase, flow rate, etc. of the product. That is, the modified first bed layer S11C can perform equilibrium reaction calculations based on the Gibbs energy minimization method, similar to the gas region equilibrium reaction calculation unit 213. The first bed layer equilibrium reaction calculation unit 214-1 calculates and outputs the change in heat quantity, flow rate, type, composition, amount, phase, etc. of the product when the first bed layer reaction component SR1 and the first inflow gas reaction component GR111 reach equilibrium.

[0098] The first bed layer equilibrium reaction calculation unit 214-1 calculates the product gas generated by the reaction of the first bed layer reaction SR1 with the first inflow gas reaction GR11 and the unused portion of the first bed layer reaction SR1 as the first bed layer product SP1. The first bed layer equilibrium reaction calculation unit 214-1 also calculates the product gas generated by the reaction of the first bed layer reaction SR1 with the first inflow gas reaction GR111 and the unused portion of the first inflow gas reaction GR111 as the first inflow gas product GP111.

[0099] That is, the first bed layer product SP1 is the sum of the product material produced by the reaction of the first bed layer reaction SR1 with the first inflow gas reaction GR11 and the unused portion of the first bed layer reaction SR1. The first inflow gas product GP111 is the sum of the product gas produced by the reaction of the first bed layer reaction SR1 with the first inflow gas reaction GR111 and the unused portion of the first inflow gas reaction GR111.

[0100] The second bed layer equilibrium reaction calculation unit 214-2 calculates the equilibrium reaction between the modified second bed layer S21A moving through the bed layer in region A and the second inlet gas G112 present in the bed layer. That is, the second bed layer equilibrium reaction calculation unit 214-2 performs equilibrium reaction calculation of the modified second bed layer S21A using the second bed layer reaction portion SR2 that contributes to the equilibrium reaction with the second inlet gas G112 in the modified second bed layer S21A and the second inlet gas reaction portion GR112 that contributes to the equilibrium state with the modified second bed layer S21A. As a result, at least the compositions and amounts of the second bed layer reactant SR2 and the second inlet gas reactant GR112, and the changes in the heat quantities and flow rates of the second bed layer reactant SR2 and the second inlet gas reactant GR112 when they reach equilibrium are calculated.

[0101] The second bed layer equilibrium reaction calculation unit 214-2 has a function of performing equilibrium reaction calculations to determine the type, phase, flow rate, etc. of the product (product) that may be produced by the reaction between the second bed layer reaction component SR2 and the second inflow gas reaction component GR112 so that the free energy of the product is minimized by setting the type, phase, flow rate, etc. in advance. That is, the second bed layer equilibrium reaction calculation unit 214-2 can perform equilibrium reaction calculations based on the Gibbs energy minimization method, similar to the first bed layer equilibrium reaction calculation unit 214-1. The second bed layer equilibrium reaction calculation unit 214-2 calculates and outputs the change in heat quantity, flow rate, type, composition, amount, and phase of the product when the second bed layer reaction component SR2 and the second inflow gas reaction component GR112 reach equilibrium.

[0102] The second bed layer equilibrium reaction calculation unit 214-2 calculates the product gas generated by the reaction of the second bed layer reaction SR2 with the second inflow gas reaction GR112 and the unused portion of the second bed layer reaction SR2 as the second bed layer product SP2. The second bed layer equilibrium reaction calculation unit 214-2 also calculates the product gas generated by the reaction of the second bed layer reaction SR2 with the second inflow gas reaction GR112 and the unused portion of the second inflow gas reaction GR112 as the second inflow gas product GP112.

[0103] That is, the second bed layer product SP2 is the sum of the product material produced by the reaction of the second bed layer reaction SR2 with the second inflow gas reaction GR112 and the unused portion of the second bed layer reaction SR2. The second inflow gas product GP112 is the sum of the product gas produced by the reaction of the second bed layer reaction SR2 with the second inflow gas reaction GR112 and the unused portion of the second inflow gas reaction GR112.

[0104] The mixed gas phase calculation unit 215 has the function of mixing multiple flows, and calculates the flow rate and composition data of the third gas phase G3, which is a mixed gas phase obtained by mixing the unreacted gas portion Gr1 separated by the gas reaction amount calculation unit 207 and the first gas product portion GP1 generated by the gas region equilibrium reaction calculation unit 213.

[0105] The first mixed substance calculation unit 216 has a function of mixing a plurality of flows, and calculates the flow rate, composition data, etc. of the first mixed substance M11 obtained by mixing the unreacted portion of the first substance Mr1 divided by the first substance reaction amount calculation unit 208 and the produced portion of the first substance MP1 generated by the gas region equilibrium reaction calculation unit 213.

[0106] The first mixed inflow gas calculation unit 217 has the function of mixing multiple flows, and calculates the flow rate and composition data of the first mixed inflow gas G111A, which is a mixture of the first inflow gas unreacted portion Gr111 separated in the first inflow gas reaction amount calculation unit 210 and the first inflow gas product portion GP111 generated in the first bed layer equilibrium reaction calculation unit 214-1.

[0107] The first mixed bed layer calculation unit 218 has a function of mixing multiple flows, and calculates the flow rate, composition data, etc. of the first mixed bed layer S12 obtained by mixing the first bed layer unreacted portion Sr1 separated by the first bed layer reaction amount calculation unit 209 and the first bed layer product portion SP1 generated by the first bed layer equilibrium reaction calculation unit 214-1.

[0108] The second mixed inflow gas calculation unit 219 has the function of mixing multiple flows, and calculates the flow rate and composition data of the second mixed inflow gas G112A, which is a mixture of the second inflow gas unreacted portion Gr112 separated in the second inflow gas reaction amount calculation unit 212 and the second inflow gas product portion GP112 generated in the second bed layer equilibrium reaction calculation unit 214-2.

[0109] The second mixed bed layer calculation unit 220 has a function of mixing multiple flows, and calculates the flow rate, composition data, etc. of the second mixed bed layer S22 obtained by mixing the second bed layer unreacted portion Sr2 separated by the second bed layer reaction amount calculation unit 211 and the second bed layer product portion SP2 generated by the second bed layer equilibrium reaction calculation unit 214-2.

[0110] The inflow gas combining unit 221 has the function of mixing multiple flows, and combines the first mixed inflow gas G111A generated in the first mixed inflow gas calculation unit 217 and the second mixed inflow gas G112A generated in the second mixed inflow gas calculation unit 219 to calculate the mixed inflow gas G12.

[0111] The gas amount summing unit 222 has the function of mixing multiple flows, and sums the third gas phase G3 generated in the mixed gas phase calculation unit 215 and the mixed inflow gas G12 generated in the inflow gas summing unit 221 to calculate the third mixed gas phase G3'.

[0112] The first mixed substance quantity correction unit 223 corrects the quantity of the first mixed substance M11 in consideration of the amount of movement of the first mixed substance M11 to the first bed layer S11 generated in the first mixed substance calculation unit 216, and calculates the first corrected mixed substance M12.

[0113] The amount of the first mixture material M11 is corrected taking into account the amount of the fourth material M4 contained in the first bed layer S11A that moves to the first mixture material M11 present in the gas region and the amount of the fifth material M5 contained in the first mixture material M11 that moves to the first bed layer S11A.

[0114] In this embodiment, at least one of the first bed layer equilibrium reaction calculation unit 214-1 and the second bed layer equilibrium reaction calculation unit 214-2 may be provided in two or more units. At least one of the first bed layer S11 and the first bed layer S21 is further separated into a plurality of bed layers, so that the bed layers in the rotary kiln 1 are divided into three or more types of multiple systems.

[0115] In this embodiment, the calculation device 20 for furnace reactions does not need to include the distribution unit 201 when the combustion material 230 is not dropped into the rotary kiln 1. In this case, the gas reaction portion is the reaction portion of the first gas phase G1 that contributes to the equilibrium reaction with the bed layer S1, and the bed reaction portion is the reaction portion of the bed layer S1 that contributes to the equilibrium state with the first gas phase G1.

[0116] In this embodiment, the reactor reaction calculation device 20 does not need to include one or more of the first bed layer preliminary quantity correcting unit 203, the first bed layer quantity correcting unit 204, and the second bed layer quantity correcting unit 205. For example, when the reactor reaction calculation device 20 does not include the first bed layer preliminary quantity correcting unit 203 and the first bed layer quantity correcting unit 204, the first bed layer reaction quantity calculation unit 209 uses the first bed layer S11 instead of the corrected first bed layer S11C. Furthermore, when the reactor reaction calculation device 20 does not include the second bed layer quantity correcting unit 205, the second bed layer reaction quantity calculation unit 211 uses the second bed layer S21 instead of the corrected second bed layer S21A.

[0117] In this embodiment, the reactor reaction calculation device 20 may calculate the heat conduction in the reactor using a model of radiation, conduction, convection, or the like, as necessary.

[0118] <Calculation method for reactor reactions> Next, a method for calculating furnace reactions according to this embodiment will be described using the furnace reaction calculation device according to this embodiment. In the method for calculating furnace reactions according to this embodiment, raw ore is supplied from the charging end 14A of the rotary kiln 1 having the configuration shown in Figure 1, and while the raw ore is being moved toward the discharge end 14B, combustion material 230 is introduced midway through the movement, and the raw ore is brought into contact with combustion gas supplied from a burner 16 provided on the discharge end 14B side, whereby it is dried and reduced.

[0119] 5 is a flowchart illustrating a method for calculating a furnace reaction according to this embodiment. As shown in FIG. 5, the furnace reaction calculation device 20 checks whether or not the combustion material 230 has been dropped into the rotary kiln 1 (checking step: step S11).

[0120] If combustion material 230 has been dropped into the rotary kiln 1 (step S11: Yes), the furnace reaction calculation device 20 provides materials such as volatile matter assumed to be in the gas phase and solid materials assumed to be in the bed layer as input materials, and uses the distribution unit 201 to distribute the combustion material 230 added into the rotary kiln 1 between the gas phase and the bed layer at a mass flow rate (distribution process: step S12).

[0121] Next, the calculation device 20 of the in-furnace reaction receives as input materials a first gas phase G1, which is a combustion gas flowing into the region A from one adjacent region (region (A+1)), and an added gas phase AG11, which is the gas phase of the combustion material 230 distributed by the distribution unit 201. The calculation device 20 of the in-furnace reaction mixes the first gas phase G1 and the added gas phase AG11, and calculates a second gas phase G2 containing the first gas phase G1 and the added gas phase AG11 (gas phase mixing process of the first gas phase and the added gas phase: step S13).

[0122] Next, the reactor reaction calculation device 20 uses the gas mixing amount calculation unit 202 to calculate the flow rate of the inflow gas G11 moving from the gas region to the bed layer from the second gas phase G2 including the first gas phase G1 and the added gas phase AG11 (gas mixing amount calculation process: step S14).

[0123] Next, the calculation device 20 of the in-furnace reaction provides as input materials the first bed layer S11, into which the bed layer S1 flowing into the region A from the other adjacent region (region (A-1)) has been distributed in advance, and the additional bed layer AS11, which is a bed layer of the combustion material distributed by the distribution unit 201. The calculation device 20 of the in-furnace reaction mixes the first bed layer S11 and the additional bed layer AS11 to calculate the first bed layer S11A (bed layer mixing step of the first bed layer and the additional bed layer: step S15).

[0124] Next, the reactor reaction calculation device 20 provides the first bed layer S11A, the amount of the fourth substance M4 contained in the first bed layer S11A that has moved to the gas region, and the amount of the fifth substance M5 contained in the first mixed substance quantity correction unit 223 that has moved to the first bed layer S11A as input substances.

[0125] The reactor reaction calculation device 20 corrects the quantity of the first bed layer S11A by using the first bed layer preliminary quantity correction unit 203 to subtract the amount of the fourth substance M4 contained in the first bed layer S11A that has moved to the gas region and add the amount of the fifth substance M5 contained in the first mixed substance quantity correction unit 223 that has moved to the first bed layer S11A, thereby calculating the preliminary corrected first bed layer S11B (first bed layer preliminary quantity correction step: step S16).

[0126] Next, the reactor reaction calculation device 20 provides as input materials the preliminary correction first bed layer S11B, the amount of the fourth substance M4 contained in the preliminary correction first bed layer S11B moving to the second bed layer S21, and the fifth substance M5 in the second bed layer S21 divided by the second bed layer material quantity correction unit 205.

[0127] The reactor reaction calculation device 20 subtracts the movement amount of the second substance M2 contained in the preliminary correction first bed layer S11B from the preliminary correction first bed layer S11B using the first bed layer quantity correction unit 204, and adds the movement amount of the third substance M3 in the second bed layer S21 divided by the second bed layer quantity correction unit 205 to the preliminary correction first bed layer S11B. As a result, the reactor reaction calculation device 20 corrects the quantity of the preliminary correction first bed layer S11B and calculates the corrected first bed layer S11C (preliminary correction first bed layer quantity correction step: step S17).

[0128] Next, the reactor reaction calculation device 20 provides the second bed layer S21, the movement amount of the second substance M2 in the preliminary corrected first bed layer S11B divided by the first bed layer quantity correction unit 204, and the third substance M3 in the second bed layer S21 divided by the second bed layer quantity correction unit 205 as input materials.

[0129] The reactor reaction calculation device 20 uses the second bed layer quantity correction unit 205 to subtract the movement amount of the third substance M3 contained in the second bed layer S21 from the second bed layer S21, and adds the movement amount of the second substance M2 divided by the first bed layer quantity correction unit 204 to the second bed layer S21. As a result, the reactor reaction calculation device 20 corrects the quantity of the second bed layer S21 and calculates the corrected second bed layer S21A (second bed layer quantity correction step: step S18).

[0130] Next, the reactor reaction calculation device 20 performs a gas region equilibrium reaction (a gas region equilibrium reaction execution step: step S19).

[0131] Next, the reactor reaction calculation device 20 performs a first bed layer equilibrium reaction (first bed layer equilibrium reaction execution step: step S20).

[0132] Next, the reactor reaction calculation device 20 performs the second bed layer equilibrium reaction (step of performing the second bed layer equilibrium reaction: step S21).

[0133] Next, the reactor reaction calculation device 20 uses the inflow gas combining unit 221 to combine the first mixed inflow gas G111A generated in the first bed layer equilibrium reaction implementation process (step S20) with the second mixed inflow gas G111B generated in the second bed layer equilibrium reaction implementation process (step S21) to obtain the mixed inflow gas G12 (inflow gas combining process: step S22).

[0134] Next, the reactor reaction calculation device 20 uses the gas amount summation unit 222 to sum the third gas phase G3 generated in the gas region equilibrium reaction implementation process (step S19) with the mixed inflow gas G12 generated in the first mixed inflow gas and second mixed inflow gas summation process (step S22) to obtain a third combined gas phase G3' (gas amount summation process: step S23).

[0135] The reactor reaction calculation device 20 moves the third combined gas phase G3' to a region (region (A-1)) closer to the charging end 14A than the region A.

[0136] Next, the reactor reaction calculation device 20 uses the first mixture substance quantity correction unit 223 to correct the quantity of the first mixture substance M11 by dividing the first mixture substance M11 produced in the gas region equilibrium reaction implementation process (step S19) by the amount of movement of the fourth substance M4 produced in the first bed layer preliminary quantity correction process (step S16) and adding the amount of movement of the fifth substance M5 from the first mixture substance M11 to the first bed layer S11, thereby calculating the first corrected mixture substance M12 (first mixture substance quantity correction process: step S24).

[0137] The calculation device 20 for the in-furnace reaction moves the first corrected mixture M12 to a region (region (A-1)) closer to the charging end 14A than the region A.

[0138] Next, the gas region equilibrium reaction implementation step (step S19) will be described. Fig. 6 is a flowchart showing the operation of the gas region equilibrium reaction implementation step (step S19) in Fig. 5. As shown in Fig. 6, in the gas region equilibrium reaction implementation step (step S19), the reactor reaction calculation device 20 uses the gas reaction amount calculation unit 207 to calculate the gas reaction amount that contributes to the equilibrium state with the first substance reaction amount MR1 in the second gas phase G2 (gas reaction amount calculation step: step S191).

[0139] That is, assuming that the second gas phase G2 reacts with the first substance M1 present in the gas region to reach an equilibrium state, the furnace reaction calculation device 20 calculates the amount of gas reaction that contributes to the reaction when the second gas phase G2 reacts with the first substance M1 present in the gas region to reach an equilibrium state.

[0140] Then, based on the calculation result of the gas reaction amount, the reactor reaction calculation device 20 divides the gas reaction amount into a gas reaction amount GR1 corresponding to the gas reaction amount and an unreacted gas amount Gr1 corresponding to the remaining unreacted amount.

[0141] Next, the reactor reaction calculation device 20 uses the first substance reaction amount calculation unit 208 to calculate the first substance reaction amount of the first substance M1 present in the gas region that contributes to the equilibrium state with the second gas phase G2 (first substance reaction amount calculation process: step S192).

[0142] That is, assuming that the first substance M1 present in the gas region reacts with the second gas phase G2 to reach an equilibrium state, the reactor reaction calculation device 20 calculates the amount of reaction of the first substance that contributes to the reaction when the first substance M1 present in the gas region reacts with the second gas phase G2 to reach an equilibrium state.

[0143] Then, based on the calculation result of the first substance reaction amount, the reactor reaction calculation device 20 divides it into a first substance reaction amount MR1 corresponding to the first substance reaction amount and a first substance unreacted amount Mr1 corresponding to the remaining unreacted amount.

[0144] Next, the reactor reaction calculation device 20 receives as input materials the gas reaction amount GR1 obtained in the gas reaction amount calculation step (step S191) and the first substance reaction amount MR1 obtained in the first substance reaction amount calculation step (step S192). The reactor reaction calculation device 20 calculates the equilibrium reaction between the gas reaction amount GR1 and the first substance reaction amount MR1 using the gas region equilibrium reaction calculation unit 213, and calculates at least the change in heat quantity and flow rate of each when the gas reaction amount GR1 and the first substance reaction amount MR1 reach an equilibrium state (gas region equilibrium reaction calculation step: step S193).

[0145] For example, the Gibbs energy minimization method can be used to calculate the equilibrium reaction.

[0146] The reactor reaction calculation device 20 calculates the gas phase produced by the reaction between the gas reaction amount GR1 and the first substance reaction amount MR1 and the unused portion of the gas reaction amount GR1 as the first gas product amount GP1, and calculates the product material produced by the reaction between the gas reaction amount GR1 and the first substance reaction amount MR1 and the unused portion of the first substance reaction amount MR1 as the first substance product amount MP1.

[0147] Next, the reactor reaction calculation device 20 receives as input the unreacted gas Gr1 obtained in the gas reaction amount calculation step (step S191) and the gas product GP1 obtained in the gas region equilibrium reaction calculation step (step S193). The reactor reaction calculation device 20 uses the mixed gas phase calculation unit 215 to calculate the flow rate, composition data, etc. of the third gas phase G3 containing the unreacted gas Gr1 and the first gas product GP1 (mixed gas phase calculation step: step S194).

[0148] The gas reaction amount GR1 calculated in the gas reaction amount calculation step (step S191) is the reaction amount when it is assumed that the second gas phase G2 and the first substance M1 have reached equilibrium, and therefore is normally used in its entirety. However, the gas phase gas reaction amount GR1 may contain substances that exist in excess of the stoichiometric ratio of the reaction or inactive substances that do not contribute to the reaction. Examples of inactive substances include nitrogen. In the mixed gas phase calculation step (step S194), substances that exist in excess of the stoichiometric ratio of the reaction and that result in remaining after the equilibrium reaction or inactive substances exist. The unused portion that remains unused in the gas region equilibrium reaction calculation unit 213 is calculated as the gas production amount GP1.

[0149] Next, the in-furnace reaction calculation device 20 receives as input materials the unreacted first substance Mr1 obtained in the first substance reaction amount calculation step (step S192) and the first substance product MP1 obtained in the gas region equilibrium reaction calculation step (step S193). The in-furnace reaction calculation device 20 uses the first mixture calculation unit 216 to calculate the flow rate, composition data, etc. of the first mixture M11 obtained by mixing the unreacted first substance Mr1 and the first substance product MP1 (first mixture calculation step: step S195).

[0150] Next, the first bed layer equilibrium reaction execution step (step S20) will be described. Fig. 7 is a flowchart showing the operation of the first bed layer equilibrium reaction execution step (step S20) in Fig. 5. As shown in Fig. 7, in the first bed layer equilibrium reaction execution step (step S20), the reactor reaction calculation device 20 uses the first bed layer reaction amount calculation unit 209 to calculate the first bed layer reaction amount that contributes to the equilibrium state with the first inflow gas reaction amount GR111 in the corrected first bed layer S11C using the first reaction rate equation (first bed layer reaction amount calculation step: step S201).

[0151] That is, assuming that the modified first bed layer S11C and the first inflow gas G111 react to reach an equilibrium state, the calculation device 20 for the in-furnace reaction calculates, using the first reaction rate equation, the first bed layer reaction amount contributing to the reaction when the modified first bed layer S11C reacts with the first inflow gas G111 to reach an equilibrium state.

[0152] Then, based on the calculation result of the first bed layer reaction amount, the reactor reaction calculation device 20 divides it into a first bed layer reaction amount SR1 corresponding to the first bed layer reaction amount and a first bed layer unreacted amount Sr1 corresponding to the remaining unreacted amount.

[0153] Next, the reactor reaction calculation device 20 uses the first inflow gas reaction amount calculation unit 210 to calculate the first inflow gas reaction amount of the first inflow gas G111 that contributes to the equilibrium state with the modified first bed layer S11C (first inflow gas reaction amount calculation process: step S202).

[0154] That is, assuming that the first inflow gas G111 and the modified first bed layer S11C react to reach an equilibrium state, the reactor reaction calculation device 20 calculates the reaction amount of the first inflow gas that contributes to the reaction when the first inflow gas G111 reacts with the modified first bed layer S11C to reach an equilibrium state.

[0155] Then, based on the calculation result of the first inflow gas reaction amount, the furnace reaction calculation device 20 divides it into a first inflow gas reaction amount GR111 corresponding to the first inflow gas reaction amount and a first inflow gas unreacted amount Gr111 corresponding to the remaining unreacted amount.

[0156] Next, the reactor reaction calculation device 20 receives as input the first bed layer reaction amount SR1 obtained in the first bed layer reaction amount calculation step (step S201) and the first inflow gas reaction amount GR111 obtained in the first inflow gas reaction amount calculation step (step S202). The reactor reaction calculation device 20 uses the first bed layer equilibrium reaction calculation unit 214-1 to calculate the equilibrium reaction between the first bed layer reaction amount SR1 and the first inflow gas reaction amount GR111, and calculates at least the changes in the heat quantities and the flow rates of the first bed layer reaction amount SR1 and the first inflow gas reaction amount GR111 when they reach equilibrium (first bed layer equilibrium reaction calculation step S203).

[0157] The reactor reaction calculation device 20 calculates the bed layer generated by the reaction of the first bed layer reaction portion SR1 and the first inflow gas reaction portion GR111 and the unused portion of the first bed layer reaction portion SR1 as the first bed layer product portion SP1. The reactor reaction calculation device 20 also calculates the gas phase generated by the reaction of the first bed layer reaction portion SR1 and the first inflow gas reaction portion GR111 and the unused portion of the first inflow gas reaction portion GR111 as the first inflow gas product portion GP111.

[0158] Next, the reactor reaction calculation device 20 receives as input materials the first bed layer unreacted fraction Sr1 obtained in the first bed layer reaction amount calculation step (step S201) and the first bed layer product fraction SP1 obtained in the first bed layer equilibrium reaction calculation step (step S203).The reactor reaction calculation device 20 uses the first bed layer equilibrium reaction calculation unit 214-1 to calculate the flow rate, composition data, etc. of the first mixed bed layer S12 containing the first bed layer unreacted fraction Sr1 and the first bed layer product fraction SP1 (first mixed bed layer calculation step: step S204).

[0159] The first bed layer reaction amount SR1 of the first bed layer calculated in the first bed layer reaction amount calculation step (step S201) may contain substances that exist in an amount greater than the stoichiometric ratio of the reaction and inactive substances that do not contribute to the reaction. In the first mixed bed layer calculation step (step S204), the first bed layer reaction amount SR1, that is present in an amount greater than the stoichiometric ratio of the reaction and is not used in the reaction as a result, and the inactive portion are calculated and included in the first bed layer product amount SP1 as an unused portion that remains unused in the first bed layer equilibrium reaction calculation step (step S203).

[0160] The reactor reaction calculation device 20 moves the first mixed bed layer S12 to a region (region (A+1)) closer to the discharge end 14B than the region A.

[0161] Next, the reactor reaction calculation device 20 receives as input materials the first inlet gas unreacted portion Gr111 obtained in the first inlet gas reaction amount calculation step (step S202) and the first inlet gas product portion GP111 obtained in the first bed layer equilibrium reaction calculation step (step S203). The reactor reaction calculation device 20 uses the first mixed inlet gas calculation unit 217 to calculate the flow rate, composition data, etc. of the first mixed inlet gas G111A containing the first inlet gas unreacted portion Gr111 and the first inlet gas product portion GP111 (first mixed inlet gas calculation step: step S205).

[0162] Next, the second bed layer equilibrium reaction execution step (step S21) will be described. Fig. 8 is a flowchart showing the operation of the second bed layer equilibrium reaction execution step (step S21) in Fig. 3. As shown in Fig. 8, in the second bed layer equilibrium reaction execution step (step S21), the reactor reaction calculation device 20 uses the second bed layer reaction amount calculation unit 211 to calculate the second bed layer reaction amount that contributes to the equilibrium state with the second inflow gas reaction amount GR112 in the corrected second bed layer S21A using the second reaction rate equation (second bed layer reaction amount calculation step: step S211).

[0163] That is, assuming that the modified second bed layer S21A and the second inflow gas G112 react to reach an equilibrium state, the calculation device 20 for the in-furnace reaction calculates, using the second reaction rate equation, the amount of second bed layer reaction that contributes to the reaction when the modified second bed layer S21A reacts with the second inflow gas G112 to reach an equilibrium state.

[0164] Then, based on the calculation result of the second bed layer reaction amount, the reactor reaction calculation device 20 divides it into a second bed layer reaction amount SR2 corresponding to the second bed layer reaction amount and a second bed layer unreacted amount Sr2 corresponding to the remaining unreacted amount.

[0165] Next, the reactor reaction calculation device 20 uses the second inflow gas reaction amount calculation unit 212 to calculate the second inflow gas reaction amount of the second inflow gas G112 that contributes to the equilibrium state with the modified second bed layer S21A (second inflow gas reaction amount calculation process: step S212).

[0166] That is, assuming that the second inflow gas G112 reacts with the modified second bed layer S21A to reach an equilibrium state, the reactor reaction calculation device 20 calculates the amount of second inflow gas reaction that contributes to the reaction when the second inflow gas G112 reacts with the modified second bed layer S21A to reach an equilibrium state.

[0167] Then, based on the calculation result of the second inflow gas reaction amount, the furnace reaction calculation device 20 divides it into a second inflow gas reaction amount GR112 corresponding to the second inflow gas reaction amount and a second inflow gas unreacted amount Gr112 corresponding to the remaining unreacted amount.

[0168] Next, the reactor reaction calculation device 20 receives as input materials the second bed layer reaction amount SR2 obtained in the second bed layer reaction amount calculation step (step S211) and the second inflow gas reaction amount GR112 obtained in the second inflow gas reaction amount calculation step (step S212). The reactor reaction calculation device 20 uses the second bed layer equilibrium reaction calculation unit 214-2 to calculate the equilibrium reaction between the second bed layer reaction amount SR2 and the second inflow gas reaction amount GR112, and calculates at least the changes in the heat quantities and the flow rates of the second bed layer reaction amount SR2 and the second inflow gas reaction amount GR112 when they reach equilibrium (first bed layer equilibrium reaction calculation step: step S213).

[0169] The reactor reaction calculation device 20 calculates the bed layer generated by the reaction of the second bed layer reaction portion SR2 with the second inflow gas reaction portion GR112 and the unused portion of the second bed layer reaction portion SR2 as the second bed layer product portion SP2. The reactor reaction calculation device 20 also calculates the gas phase generated by the reaction of the second bed layer reaction portion SR2 with the second inflow gas reaction portion GR112 and the unused portion of the second inflow gas reaction portion GR112 as the second inflow gas product portion GP112.

[0170] Next, the reactor reaction calculation device 20 receives as input materials the second bed layer unreacted fraction Sr2 obtained in the second bed layer reaction amount calculation step (step S211) and the second bed layer product fraction SP2 obtained in the second bed layer equilibrium reaction calculation step (step S213).The reactor reaction calculation device 20 uses the second bed layer equilibrium reaction calculation unit 214-2 to calculate the flow rate, composition data, etc. of the second mixed bed layer S22 containing the second bed layer unreacted fraction Sr2 and the second bed layer product fraction SP2 (second mixed bed layer calculation step: step S214).

[0171] The second bed layer reaction amount SR2 of the second bed layer calculated in the second bed layer reaction amount calculation step (step S211) may contain substances that exist in an amount greater than the stoichiometric ratio of the reaction and inactive substances that do not contribute to the reaction. In the second mixed bed layer calculation step (step S214), the second bed layer reaction amount SR2 contains a portion that exists in an amount greater than the stoichiometric ratio of the reaction and is not used in the reaction as a result, and an inactive portion. The unused portion that remains unused in the second bed layer equilibrium reaction calculation step (step S213) is calculated as the second bed layer product amount SP2.

[0172] The reactor reaction calculation device 20 moves the second mixed bed layer S22 to a region (region (A+1)) closer to the discharge end 14B than the region A.

[0173] Next, the reactor reaction calculation device 20 receives as input materials the second inlet gas unreacted portion Gr112 obtained in the second inlet gas reaction amount calculation step (step S212) and the second inlet gas product portion GP112 obtained in the second bed layer equilibrium reaction calculation step (step S213). The reactor reaction calculation device 20 uses the second mixed inlet gas calculation unit 219 to calculate the flow rate, composition data, etc. of the second mixed inlet gas G112A containing the second inlet gas unreacted portion Gr112 and the second inlet gas product portion GP112 (second mixed inlet gas calculation step: step S215).

[0174] In the method for calculating the furnace reaction according to this embodiment, the steps may be performed in parallel or the order of the steps may be changed as appropriate.

[0175] For example, the method for calculating the furnace reaction according to this embodiment may perform at least one of the steps of carrying out the gas region equilibrium reaction (step S19), the step of carrying out the first bed layer equilibrium reaction (step S20), and the step of carrying out the second bed layer equilibrium reaction (step S21) in parallel.

[0176] In the method for calculating the furnace reaction according to this embodiment, the inflow gas summing step (step S22) and the first mixture substance amount correction step (step S24) may be performed in parallel, or the inflow gas summing step (step S22) may be performed after the first mixture substance amount correction step (step S24).

[0177] In the method for calculating the furnace reaction according to this embodiment, the gas amount summing step (step S23) and the first mixture substance amount correction step (step S24) may be performed in parallel, or the gas amount summing step (step S23) may be performed after the first mixture substance amount correction step (step S24).

[0178] In the method for calculating the furnace reaction according to this embodiment, the gas reaction amount calculation step (step S191) and the first substance reaction amount calculation step (step S192) may be performed in parallel, or the gas reaction amount calculation step (step S191) may be performed after the first substance reaction amount calculation step (step S192).

[0179] In the method for calculating the furnace reaction according to this embodiment, the mixed gas phase calculation step (step S194) and the first mixed substance calculation step (step S195) may be performed in parallel, or the mixed gas phase calculation step (step S194) may be performed after the first mixed substance calculation step (step S195).

[0180] In the method for calculating the in-furnace reaction according to this embodiment, the first bed layer reaction amount calculation step (step S201) and the first inflow gas reaction amount calculation step (step S202) may be performed in parallel, or the first bed layer reaction amount calculation step (step S201) may be performed after the first inflow gas reaction amount calculation step (step S202).

[0181] In the method for calculating the in-furnace reaction according to this embodiment, the first mixed bed layer calculation step (step S204) and the first mixed inflow gas calculation step (step S205) may be performed in parallel, or the first mixed bed layer calculation step (step S204) may be performed after the first mixed inflow gas calculation step (step S205).

[0182] In the method for calculating the in-furnace reaction according to this embodiment, the second bed layer reaction amount calculation step (step S211) and the second inflow gas reaction amount calculation step (step S212) may be performed in parallel, or the second bed layer reaction amount calculation step (step S211) may be performed after the second inflow gas reaction amount calculation step (step S212).

[0183] In the method for calculating the in-furnace reaction according to this embodiment, the second mixed bed layer calculation step (step S214) and the second mixed inflow gas calculation step (step S215) may be performed in parallel, or the second mixed bed layer calculation step (step S214) may be performed after the second mixed inflow gas calculation step (step S215).

[0184] In the method for calculating the reaction inside the furnace according to this embodiment, the heat conduction inside the furnace may be calculated using a model of radiation, conduction, convection, or the like, as necessary.

[0185] <Hardware configuration of reactor reaction calculation device> Next, an example of the hardware configuration of a reactor reaction calculation device will be described. FIG. 9 is a hardware configuration diagram of the reactor reaction calculation device. As shown in FIG. 9, the reactor reaction calculation device 20 is configured, for example, by an information processing device (computer), and can be physically configured as a computer system including a CPU (Central Processing Unit: processor) 21, which is an arithmetic processing unit, a RAM (Random Access Memory) 22 and a ROM (Read Only Memory) 23, which are main storage devices, an auxiliary storage device 24, an input / output interface 25, a display device 26, which is an output device, and the like. These are connected to each other by a bus 27. The auxiliary storage device 24 and the display device 26 may be provided externally.

[0186] The CPU 21 controls the overall operation of the furnace reaction calculation device 20 and performs various information processing. The CPU 21 executes a raw material ore reaction calculation program stored in the ROM 23 or auxiliary storage device 24, and controls the display operation of the measurement recording screen and the analysis screen.

[0187] The RAM 22 is used as a work area for the CPU 21 and may include a non-volatile RAM for storing main control parameters and information.

[0188] The ROM 23 stores a basic input / output program, etc. A raw material ore reaction calculation program may be stored in the ROM 23.

[0189] The auxiliary storage device 24 is a storage device such as an SSD (Solid State Drive) or an HDD (Hard Disk Drive), and stores, for example, a raw ore reaction calculation program and various data, files, etc. required for the operation of the furnace reaction calculation device 20.

[0190] The input / output interface 25 includes both a user interface such as a touch panel, keyboard, display screen, and operation buttons, and a communication interface that takes in information from an external data recording server or the like and outputs analysis information to other electronic devices.

[0191] The display device 26 is a monitor display, etc. The display device 26 displays a measurement recording screen and an analysis screen, and the screen is updated in response to input / output operations via the input / output interface 25.

[0192] Each function of the reactor reaction calculation device 20 shown in Figure 9 is realized by loading simulation software (including a reactor reaction calculation program) into a main memory device such as RAM 22 or ROM 23 or an auxiliary memory device 24, and executing a raw ore reaction calculation program stored in RAM 22, ROM 23 or the auxiliary memory device 24 by CPU 21, thereby reading and writing data in RAM 22, etc., and operating an input / output interface 25 and a display device 26.

[0193] The calculation program for the reactor reaction can be a program with the following configuration. That is, the calculation program for the in-furnace reaction is a program that causes at least a computer to execute calculations for an in-furnace reaction in which raw material ore supplied from one end of a reactor is brought into contact with combustion gas supplied from the other end while moving the raw material ore toward the other end, thereby drying and reducing the raw material ore, a gas mixing amount calculation step of calculating a flow rate of an inflow gas that moves from the gas region in the reactor to a bed layer containing the raw ore, among a first gas phase containing the combustion gas flowing through a gas region in the reactor; an inflow gas distribution step of dividing the inflow gas into a first inflow gas and a second inflow gas; a first bed layer reaction amount calculation step of calculating, by using a first reaction rate equation, a first bed layer reaction amount contributing to an equilibrium state with the first inflow gas in one of two types of distributed bed layers distributed among the bed layers according to reaction rate-dependent conditions of the solid and liquid substances contained in the bed layers, to obtain a first bed layer reaction amount contributing to the equilibrium reaction in the first bed layer; a second bed layer reaction amount calculation step of calculating a second bed layer reaction amount contributing to an equilibrium state with the second inlet gas in the other of the two types of distributed bed layers using a second reaction rate equation different from the first reaction rate equation to obtain a second bed layer reaction amount contributing to the equilibrium reaction in the second bed layer; a first bed layer equilibrium reaction calculation step of calculating an equilibrium reaction between the first bed layer reaction component and the first inflow gas; a second bed layer equilibrium reaction calculation step of calculating an equilibrium reaction between the second bed layer reaction component and the second inflow gas; a gas region equilibrium reaction calculation step of calculating an equilibrium reaction between the first gas phase and a first substance including at least one of a solid substance and a liquid substance present in the gas region; A program that causes a computer to execute at least the above can be used.

[0194] The calculation program for the furnace reaction is stored in a storage device provided in the computer, such as a main storage device such as RAM 22 or ROM 23, or an auxiliary storage device 24. The raw ore reaction calculation program may be configured so that a part or all of it is transmitted via a transmission medium such as a communication line, and is received and recorded (including installed) by a communication module or the like provided in the computer. The raw ore reaction calculation program may be configured so that a part or all of it is stored in a portable storage medium such as a CD-ROM, DVD-ROM, or flash memory, and then recorded (including installed) in the computer.

[0195] As described above, the reactor reaction calculation device 20 according to this embodiment includes the gas mixture amount calculation unit 202, the inflow gas distribution unit 206, the first bed layer reaction amount calculation unit 209, the second bed layer reaction amount calculation unit 211, the gas region equilibrium reaction calculation unit 213, the first bed layer equilibrium reaction calculation unit 214-1, and the second bed layer equilibrium reaction calculation unit 214-2. The reactor reaction calculation device 20 calculates the first bed layer reaction amount SR1 contributing to the equilibrium reaction from the corrected first bed layer S11C using the first reaction rate equation in the first bed layer reaction amount calculation unit 209, and calculates the second bed layer reaction amount SR2 contributing to the equilibrium reaction from the corrected second bed layer S21A using the second reaction rate equation in the second bed layer reaction amount calculation unit 211. Furthermore, the reactor reaction calculation device 20 calculates the flow rate of inflow gas G11, which moves from the gas region to the bed layer, in the second gas phase G2, in the gas mixture amount calculation unit 202, and distributes the inflow gas G11 into a first inflow gas G111 and a second inflow gas G112 in the inflow gas distribution unit 206. The reactor reaction calculation device 20 calculates the equilibrium reaction between the first bed layer reaction component S11B and the first inflow gas G111 in the first bed layer equilibrium reaction calculation unit 214-1, and calculates the equilibrium reaction between the second bed layer reaction component SR2 and the second inflow gas reaction component GR112 in the second bed layer equilibrium reaction calculation unit 214-2.

[0196] The modified first bed layer S11C is a portion of the bed layer S11 where the reaction rate is fast, and the modified second bed layer S21A is a portion of the bed layer S11 where the reaction rate is slow. The reactor reaction calculation device 20 calculates the first bed layer reaction amount SR1 using the first reaction rate equation for the portion of the bed layer S11 where the reaction rate is fast, such as the modified first bed layer S11C, and calculates the second bed layer reaction amount SR2 using the second reaction rate equation for the portion of the bed layer S11 where the reaction rate is slow, such as the modified second bed layer S21A.

[0197] Therefore, the calculation device 20 for in-furnace reactions divides the bed layers in the reactor into two series, the modified first bed layer S11C and the modified second bed layer S21A, and can calculate the equilibrium reactions between these and the inflow gas G11 based on the Gibbs energy minimization method. Therefore, the calculation device 20 for in-furnace reactions can accurately calculate the equilibrium reactions occurring in the bed layers by taking into account the reaction behavior of two types of bed layers with different reaction rates within the bed layers in the rotary kiln 1.

[0198] This allows the furnace reaction calculation device 20 to calculate the equilibrium reactions that occur in the bed layer, taking into account the areas in the bed layer where the reaction rate differs and the composition distribution in the raw ore present in the bed layer, thereby enabling fitting that is closer to the actual state inside the furnace.

[0199] Furthermore, the furnace reaction calculation device 20 can calculate the equilibrium reaction between the second gas phase G2 and the first substance M1 present in the gas region by calculating the equilibrium reaction using the gas reaction portion GR1 and the first substance reaction portion MR1 in the gas region equilibrium reaction calculation unit 213. This allows the furnace reaction calculation device 20 to calculate the equilibrium reaction between the second gas phase G2 present in the gas region and the first substance M1 based on the Gibbs energy minimization method. Therefore, the furnace reaction calculation device 20 can calculate the equilibrium reaction occurring in the gas region with high accuracy, taking into account the equilibrium reaction between the second gas phase G2 and the first substance M1 in the gas region of the rotary kiln 1.

[0200] Therefore, the furnace reaction calculation device 20 can calculate the reactions occurring in the bed layer and gas region in the rotary kiln 1 with high accuracy.

[0201] The in-furnace reaction calculation device 20 can include a first inflow gas reaction amount calculation unit 210. This allows the in-furnace reaction calculation device 20 to calculate the first inflow gas reaction amount GR111 that contributes to the equilibrium reaction from the first inflow gas G111 in the first inflow gas reaction amount calculation unit 210. Therefore, the in-furnace reaction calculation device 20 can reliably perform equilibrium reaction calculation with the modified first bed layer S11C in the first bed layer equilibrium reaction calculation unit 214-1 using the necessary first inflow gas reaction amount GR111.

[0202] The reactor reaction calculation device 20 may include a first mixed inlet gas calculation unit 217 and a first mixed bed layer calculation unit 218. As a result, the reactor reaction calculation device 20 can use the first mixed inlet gas calculation unit 217 to calculate at least the flow rate of the first mixed inlet gas G111A obtained by mixing the first inlet gas unreacted portion Gr111 and the first inlet gas product portion GP111. The reactor reaction calculation device 20 can use the first mixed bed layer calculation unit 218 to calculate at least the flow rate of the first mixed bed layer S12 obtained by mixing the first bed layer unreacted portion Sr1 and the first bed layer product portion SP1. As a result, the reactor reaction calculation device 20 can accurately calculate the amount of movement of the first inlet gas G111 and the amount of movement of the first mixed bed layer S12 to the discharge end 14B.

[0203] The in-furnace reaction calculation device 20 can include a second inflow gas reaction amount calculation unit 212. This allows the in-furnace reaction calculation device 20 to calculate the second inflow gas reaction amount GR112 contributing to the equilibrium reaction from the second inflow gas G112 in the second inflow gas reaction amount calculation unit 212. Therefore, the in-furnace reaction calculation device 20 can reliably perform equilibrium reaction calculation with the modified second bed layer S21A in the second bed layer equilibrium reaction calculation unit 214-2 using the necessary second inflow gas reaction amount GR112.

[0204] The reactor reaction calculation device 20 may include a second mixed inlet gas calculation unit 219 and a second mixed bed layer calculation unit 220. As a result, the reactor reaction calculation device 20 can use the second mixed inlet gas calculation unit 219 to calculate at least the flow rate of the second mixed inlet gas G112A obtained by mixing the second inlet gas unreacted portion Gr112 and the second inlet gas product portion GP112. The reactor reaction calculation device 20 can use the second mixed bed layer calculation unit 220 to calculate at least the flow rate of the second mixed bed layer S22 obtained by mixing the second bed layer unreacted portion Sr2 and the second bed layer product portion SP2. As a result, the reactor reaction calculation device 20 can accurately calculate the amount of movement of the second inlet gas G112 and the amount of movement of the second mixed bed layer S22 to the discharge end 14B.

[0205] The reactor reaction calculation device 20 can include a first bed layer quantity correction unit 204. This allows the reactor reaction calculation device 20 to subtract the amount of the second substance M2 present in the preliminary correction first bed layer S11B transferred to the second bed layer S21 and add the amount of the third substance M3 present in the second bed layer S21 transferred to the preliminary correction first bed layer S11B. Therefore, the reactor reaction calculation device 20 can more accurately calculate the flow rate of the corrected first bed layer S11C used in the calculation of the equilibrium reaction in the first bed layer equilibrium reaction calculation unit 214-1, and can further improve the calculation accuracy of the equilibrium reaction of the corrected first bed layer S11C in the first bed layer S11.

[0206] The reactor reaction calculation device 20 can include a second bed layer quantity correction unit 205. This allows the reactor reaction calculation device 20 to subtract the amount of the third substance M3 present in the second bed layer S21 transferred to the corrected first bed layer S11C and add the amount of the second substance M2 present in the preliminary corrected first bed layer S11B transferred to the second bed layer S21. Therefore, the reactor reaction calculation device 20 can accurately determine the flow rate and the like of the corrected second bed layer S21A used for calculating the equilibrium reaction in the second bed layer equilibrium reaction calculation unit 214-2, thereby further improving the calculation accuracy of the equilibrium reaction of the corrected second bed layer S21A in the second bed layer S21.

[0207] The reactor reaction calculation device 20 can include a gas reaction amount calculation unit 207 and a first substance reaction amount calculation unit 208. As a result, the reactor reaction calculation device 20 can calculate, in the gas reaction amount calculation unit 207, a gas reaction amount GR1 that contributes to the equilibrium reaction from the second gas phase G2 including the first gas phase G1, and can calculate, in the first substance reaction amount calculation unit 208, a first substance reaction amount MR1 that contributes to the equilibrium reaction from the first substance M1 present in the gas region. Therefore, the reactor reaction calculation device 20 can calculate the reaction amounts to be used in the gas region equilibrium reaction calculation unit 213 from the second gas phase G2 and the first substance M1 present in the gas region, and can reliably perform equilibrium reaction calculation in the gas region using the necessary reaction amounts.

[0208] The in-furnace reaction calculation device 20 can include a mixed gas phase calculation unit 215 and a first mixed substance calculation unit 216. The mixed gas phase calculation unit 215 of the in-furnace reaction calculation device 20 can calculate at least the flow rate of the third gas phase G3 obtained by mixing the unreacted gas portion Gr1 and the gas product portion GP1. The first mixed substance calculation unit 216 of the in-furnace reaction calculation device 20 can calculate at least the flow rate of the first mixed substance M11 obtained by mixing the unreacted first substance portion Mr1 and the first substance product portion MP1. This allows the in-furnace reaction calculation device 20 to more accurately calculate the amount of movement of the first gas phase G1 and the first substance M1 to the charging end 14A.

[0209] The in-furnace reaction calculation device 20 can include a first mixture substance amount correction unit 223. This allows the in-furnace reaction calculation device 20 to include the fifth substance M5 contained in the first mixture substance M11 present in the gas region in the first bed layer S11A, and also to include the fourth substance M4 present in the first bed layer S11A in the first mixture substance M11 present in the gas region. Therefore, the in-furnace reaction calculation device 20 can take into account the amount of the fifth substance M5 contained in the first mixture substance M11 transferred to the bed layer and the amount of the fourth substance M4 contained in the first bed layer S11A transferred to the first mixture substance M11, and can therefore accurately calculate the amount of the first mixture substance M11 present in region A transferred to region (A-1).

[0210] The reactor reaction calculation device 20 can include a first bed layer reserve quantity correction unit 203. This allows the reactor reaction calculation device 20 to include the fifth substance M5 contained in the first mixture substance M11 present in the gas region in the first bed layer S11A, and to include the fourth substance M4 present in the first bed layer S11A in the first mixture substance M11 present in the gas region. Therefore, the reactor reaction calculation device 20 can take into account the amount of the fifth substance M5 contained in the first mixture substance M11 transferred to the bed layer and the amount of the fourth substance M4 contained in the first bed layer S11A transferred to the first mixture substance M11, and can accurately calculate the amount used for the equilibrium reaction of the first bed layer S11A present in region A. Therefore, the reactor reaction calculation device 20 can accurately calculate the first bed layer reaction amount SR1 used in the first bed layer equilibrium reaction calculation unit 214-1, and can perform more accurate first bed layer equilibrium reaction calculation.

[0211] The reactor reaction calculation device 20 can include an inflow gas summing unit 221. This allows the reactor reaction calculation device 20 to calculate the mixed inflow gas G12 obtained by summing the first mixed inflow gas G111A and the second mixed inflow gas G112A. Therefore, the reactor reaction calculation device 20 can move the gas components present in the bed layer to the gas region. Therefore, the reactor reaction calculation device 20 can move the gas components present in the bed layer to the gas region and then to the region (A-1).

[0212] The reactor reaction calculation device 20 can include a gas amount summing unit 222. This allows the reactor reaction calculation device 20 to mix the mixed inflow gas G12 generated in the bed layer into the third gas phase G3, so that the gas components present in the bed layer can also be included in the third gas phase G3. Therefore, the reactor reaction calculation device 20 can reliably move the gas components present in region A to region (A-1).

[0213] The furnace reaction calculation device 20 may include a bed layer distribution unit 224. This allows the furnace reaction calculation device 20 to divide the bed layer S1 into the first bed layer S11 and the second bed layer S21 immediately after it is supplied from the charging end 14A of the rotary kiln 1. The furnace reaction calculation device 20 can accurately calculate the masses of the first bed layer S11 and the second bed layer S21 from the region 1 including the charging end 14A to the region A, and therefore can more accurately calculate the mass flow rates used for the equilibrium reactions of the first bed layer S11 and the second bed layer S21 in the region A. Therefore, the furnace reaction calculation device 20 can more accurately calculate the first bed layer equilibrium reaction calculation and the second bed layer equilibrium reaction calculation.

[0214] The calculation device 20 for the in-furnace reaction includes a distribution unit 201, which can distribute the gas phase and bed layer contained in the combustion material 230 added during the movement of the raw ore into an additional gas phase AS11 and an additional bed layer AG11 at a mass flow rate. The calculation device 20 for the in-furnace reaction distributes the combustion material 230 added midway through the rotary kiln 1 into a first gas phase G1 and a first bed layer S1, treating them as separate fuels, and can treat the additional gas phase AG11 as a component that makes up the second gas phase G2, and the additional bed layer AS11 as a component that makes up the first bed layer S11. As a result, the furnace reaction calculation device 20 calculates the gas reaction amount GR11 from the second gas phase G2 and the first bed layer reaction amount SR11 from the modified first bed layer S11C, and the first bed layer equilibrium reaction calculation unit 214-1 can calculate the equilibrium reaction based on the reaction amount required for the equilibrium reaction from the first inflow gas G111 and the modified first bed layer S11C. Therefore, the furnace reaction calculation device 20 can analyze the behavior of materials in the rotary kiln 1 taking into account the combustion material 230 that is added midway through the rotary kiln 1, so that the furnace reaction can be calculated with high accuracy even when the combustion material 230 is used.

[0215] In this way, the furnace reaction calculation device 20 can accurately calculate the equilibrium reactions that occur in the furnace, taking into account that the bed layers in the rotary kiln 1 exhibit multiple different types of reaction behavior. Therefore, by applying the furnace reaction calculation device 20 to the entire area inside the rotary kiln 1, it is possible to accurately calculate the entire reaction process inside the rotary kiln 1.

[0216] The following describes a case where the in-furnace reaction calculation device 20 is applied to the entire rotary kiln 1. For example, as shown in FIG. 10, when the rotary kiln 1 is divided into multiple regions, and the reaction process in one region A is assumed to be a unit operation model, the in-furnace reaction calculation device 20 can model the reaction process in the rotary kiln 1 by combining unit operation models. The in-furnace reaction calculation device 20 then repeatedly performs calculations in the multiple regions of the rotary kiln 1 along the flow of raw ore or the flow of combustion gas (including the assumed flow) (see arrows in FIG. 6). The in-furnace reaction calculation device 20 then repeatedly performs calculations until the difference between the calculated value in a specified region and the previous calculated value in that region falls within a specified range.

[0217] Fig. 11 shows a flowchart when the in-furnace reaction calculation device 20 is applied to the entire rotary kiln 1. As shown in Fig. 11, the in-furnace reaction calculation device 20 models the reaction process occurring in the rotary kiln 1 by combining multiple unit operation models (modeling step: step S21).

[0218] The unit operation model is applied with the reactor reaction calculation device 20 shown in Figure 2. For each unit operation model, a gas mixture amount calculation unit 202, an inflow gas distribution unit 206, a gas region equilibrium reaction calculation unit 213, a first bed layer equilibrium reaction calculation unit 214-1, a second bed layer equilibrium reaction calculation unit 214-2, etc., which constitute the unit operation model, are prepared in advance.

[0219] Each unit operation model is interconnected along the flow of raw material ore and combustion gas (including assumed flows).

[0220] Next, the reactor reaction calculation device 20 performs calculations on the unit operation model modeled in step S21 (calculation step: step S22). Flow information is input to the unit operation model.

[0221] The flow information includes data such as the composition, flow rate, temperature, and rotation speed of the raw ore, combustion gas, and combustion material. When the flow information is input, the unit operation model performs a predetermined calculation and outputs the calculated values ​​(composition, flow rate, temperature, etc. of the raw ore, combustion gas, and combustion material). From these calculation results, values ​​used for each component of the reactor reaction calculation device 20 shown in Figure 2 are calculated. The values ​​used for each component of the reactor reaction calculation device 20 include the flow rates of various gas phases in region A, such as bed layer S1 and first gas phase G1, the flow rates of various bed layers in region A, such as first bed layer S1, the flow rates of various substances, such as the flow rate of first substance M1, and the flow rates of various reacted, unreacted, and produced components in region A, such as gas reacted component GR1, unreacted gas component Gr1, and gas produced component GP1.

[0222] In this embodiment, calculation of the unit operation model is performed starting from the unit operation model located closest to the charging end 14A.

[0223] Next, the reactor reaction calculation device 20 determines whether or not calculation has been completed up to the final unit operation model (step S23).

[0224] When calculations have been completed up to the final unit operation model (step S23: Yes), the reactor reaction calculation device 20 determines whether or not there is a previous calculation value for the unit operation model (step S24).

[0225] If there is a previous calculation value (step S24: Yes), the reactor reaction calculation device 20 compares the calculation value calculated in the calculation step (step S24) with the previous calculation value (step S25).

[0226] Next, the reactor reaction calculation device 20 determines whether or not the difference between the calculated value and the previous calculated value satisfies the convergence condition (comparison step: step S26).

[0227] The convergence condition is, for example, that the difference between the calculated value and the previous calculated value is within a range of several degrees Celsius (for example, 1 degree Celsius) or less.

[0228] If the difference between the calculated value and the previous calculated value satisfies the convergence condition (step S26: Yes), the furnace reaction calculation device 20 ends the calculation. This allows the reaction process occurring throughout the entire kiln body 11 of the rotary kiln 1 to be analyzed.

[0229] On the other hand, in step S23, if the calculation has not been completed up to the final unit operation model (step S23: No), the reactor reaction calculation device 20 proceeds to the unit operation model located in region (A+1) or region (A-1), which is another adjacent unit operation model (step S27).Then, the reactor reaction calculation device 20 performs calculations for the unit operation model located in region (A+1) or region (A-1) (step S22).

[0230] In step S24, if there is no previous calculated value (step S24: No), or in step S26, if the difference between the calculated value and the previous calculated value does not satisfy the convergence condition (step S26: No), the reactor reaction calculation device 20 moves to the first unit operation model (step S28).

[0231] Therefore, the furnace reaction calculation device 20 models the reaction process occurring in the rotary kiln 1 using a combination of multiple unit operation models, and performs calculations based on the values ​​set for each unit operation model in accordance with the connection order of the unit operation models. In this embodiment, calculations for each unit operation model are performed sequentially from the charging end 14A to the discharge end 14B of the rotary kiln 1, and then from the discharge end 14B to the charging end 14A (see Figure 6). This series of operations is then repeated until the calculated values ​​in a predetermined region A satisfy a predetermined convergence condition. As a result, calculation results such as the flow rates of combustion gas and raw ore in each region of the rotary kiln 1 are derived. This enables more accurate analysis of the behavior of each substance constituting the combustion gas and raw ore throughout the kiln body 11 of the rotary kiln 1.

[0232] In this way, the calculation device 20 for furnace reactions can calculate the reactions in the reactor with high precision, taking into account the fact that the bed layers in the reactor exhibit two different types of reaction behavior throughout the entire reaction process in the rotary kiln 1. This allows for more accurate analysis of the behavior of each substance in the rotary kiln 1 while changing the operating conditions of the rotary kiln 1 (for example, the size and rotation speed of the rotary kiln 1, the type of raw ore, and the feed amount). Therefore, the calculation device 20 for furnace reactions can be effectively used for changing the type and feed amount of raw ore, preliminarily examining equipment improvements for the rotary kiln 1, and investigating the effects of operating conditions, etc.

[0233] In this embodiment, the raw material supplied into the rotary kiln 1 may be a raw material other than raw ore.

[0234] In this embodiment, the combustion material dropped into the rotary kiln 1 shown in FIG. 1 from the middle does not necessarily contain both volatile matter and bed layer, and may contain other substances such as ash in addition to volatile matter and bed layer.

[0235] In this embodiment, any reactor other than the rotary kiln 1 shown in FIG. 1 may be used as long as it heats raw ore while moving it from the charging end 14A toward the discharging end 14B.

[0236] Although the embodiments have been described above, they are presented as examples and the present invention is not limited to the above embodiments. The above embodiments can be implemented in various other forms, and various combinations, omissions, substitutions, modifications, etc. can be made without departing from the spirit of the invention. These embodiments and their modifications are included within the scope and spirit of the invention, and are also included in the inventions and their equivalents as set forth in the claims. [Explanation of symbols]

[0237] 1. Rotary kiln 11 Kiln body 16 burners 20. Furnace reaction calculation device 201 Distribution section of combustion materials (distribution section) 202 Gas mixture calculation unit 203 First Bed Layer Reserve Quantity Correction Section 204 First Bed Layer Material Volume Correction Section 205 Second Bed Layer Mass Correction Section 206 Inlet gas distribution section 207 Gas reaction amount calculation section 208 First substance reaction amount calculation unit 209 First bed layer reaction volume calculation section 210 First inflow gas reaction amount calculation unit 211 Second bed layer reaction volume calculation section 212 Second inflow gas reaction amount calculation unit 213 Gas Region Equilibrium Reaction Calculation Section 214-1 First bed layer equilibrium reaction calculation section 214-2 Second bed layer equilibrium reaction calculation section 215 Gas-phase mixture calculation section 216 1st Mixed Material Calculation Department 217 First mixed inflow gas calculation unit 218 First Mixed Bed Layer Calculation Section 219 Second mixed inflow gas calculation section 220 Second Mixed Bed Layer Calculation Section 221 Inflow gas combining section 222 Gas volume summing unit 223 First Mixed Material Quantity Correction Section 224 Bed Layer Distribution Section G1 First gas phase G2 Second gas phase G3 Third gas phase G3' Third mixture G11 Inlet gas G111 First inlet gas G112 Second inlet gas G111A First mixed inlet gas G112A Second mixed inlet gas G12 Mixed Inlet Gas S1 bed layer S11, S11A 1st bed layer S11B Preliminary revised first bed layer S11C Revised Bed 1 S12 First mixed bed layer S21 Second Bed Layer S21A Revised Bed 2 S22 Second mixed bed layer M1 1st substance M2 second substance M3 3rd substance M4 4th substance M5 5th substance M11 1st mixed substance M12 1st Modified Mixed Substance GR1 Gas Reaction Gr1 Unreacted gas GP1 Gas Generating Components MR1 First substance reaction Mr1 Unreacted first substance MP1 First substance production component GR111 First inflow gas reaction Gr111 1st inflow gas unreacted GP111 First inflow gas generation component GR112 Second inflow gas reaction component Gr112 Unreacted second inflow gas GP112 Second inflow gas generation component SR1 First bed layer reaction Sr1 Unreacted portion of the first bed layer SP1 First bed layer product SR2 Second bed layer reaction Sr2 Unreacted portion of the second bed layer SP2 Second bed layer product

Claims

1. A calculation device for a reaction in a furnace, in which raw material ore is supplied from one end of a reactor and moved toward the other end of the reactor, and the raw material ore is brought into countercurrent contact with combustion gas supplied from the other end, thereby drying and reducing the raw material ore, a gas mixing amount calculation unit that calculates the flow rate of an inflow gas that moves from the gas region in the reactor to a bed layer that is the raw material ore, among the first gas phase that is the combustion gas flowing through the gas region in the reactor; an inflow gas distribution unit that distributes the inflow gas into a first inflow gas and a second inflow gas; a first bed layer reaction amount calculation unit that calculates, using a first reaction rate equation, a first bed layer reaction amount that contributes to an equilibrium state with the first inflow gas in one of two types of distributed bed layers among the bed layers distributed according to reaction rate-dependent conditions of the solid material and the liquid material contained in the bed layer, thereby determining a first bed layer reaction amount that contributes to the equilibrium reaction of the first bed layer; a second bed layer reaction amount calculation unit that calculates a second bed layer reaction amount that contributes to an equilibrium state with the second inflow gas in the other second bed layer of the two types of distributed bed layers using a second reaction rate equation different from the first reaction rate equation, and obtains a second bed layer reaction amount that contributes to the equilibrium reaction in the second bed layer; a first bed layer equilibrium reaction calculation unit that calculates an equilibrium reaction between the first bed layer reaction component and the first inflow gas; a second bed layer equilibrium reaction calculation unit that calculates an equilibrium reaction between the second bed layer reaction component and the second inflow gas; a gas region equilibrium reaction calculation unit that calculates an equilibrium reaction between the first gas phase and a first substance including at least one of a solid substance and a liquid substance present in the gas region based on a Gibbs energy minimization method; a first inflow gas reaction amount calculation unit that calculates a gas reaction amount that contributes to an equilibrium state of the first inflow gas with the first bed layer, and obtains a first inflow gas reaction amount that contributes to the equilibrium reaction of the first inflow gas; a first mixed inflow gas calculation unit that calculates at least a flow rate of a first mixed inflow gas obtained by mixing a first inflow gas unreacted portion, which is an unreacted portion other than the first inflow gas reacted portion generated in the first inflow gas reaction amount calculation unit, and a first inflow gas product portion generated by the reaction of the first bed layer reacted portion and the first inflow gas reacted portion in the first bed layer equilibrium reaction calculation unit; a first mixed bed layer calculation unit that calculates at least a flow rate of a first mixed bed layer obtained by mixing a first bed layer unreacted portion, which is an unreacted portion other than the first bed layer reacted portion generated in the first bed layer reaction amount calculation unit, and a first bed layer product portion generated by a reaction between the first bed layer reacted portion and the first inflow gas reacted portion in the first bed layer equilibrium reaction calculation unit; a second inflow gas reaction amount calculation unit that calculates a gas reaction amount that contributes to an equilibrium state of the second inflow gas with the second bed layer, and obtains a second inflow gas reaction amount that contributes to the equilibrium reaction of the second inflow gas; a second mixed inflow gas calculation unit that calculates at least a flow rate of a second mixed inflow gas obtained by mixing a second inflow gas unreacted portion, which is an unreacted portion other than the second inflow gas reacted portion generated in the second inflow gas reaction amount calculation unit, and a second inflow gas product portion generated by the reaction of the second bed layer reacted portion and the second inflow gas reacted portion in the second bed layer equilibrium reaction calculation unit; a second mixed bed layer calculation unit that calculates at least a flow rate of a second mixed bed layer obtained by mixing a second bed layer unreacted portion, which is an unreacted portion other than the second bed layer reacted portion generated in the second bed layer reaction amount calculation unit, and a second bed layer product portion generated by a reaction between the second bed layer reacted portion and the second inflow gas reacted portion in the second bed layer equilibrium reaction calculation unit; Equipped with the reaction rate-causing condition is both one or more properties of a solid substance contained in the bed layer selected from the group consisting of particle size, specific gravity, type, and segregation, and one or more properties of a liquid substance contained in the bed layer selected from the group consisting of specific gravity, type, and segregation, The reactor is divided into a plurality of regions, and when the reaction process in one region is assumed to be a unit operation model, the reactor calculation device is applied to each of the unit operation models, models the reaction process in the reactor by combining the unit operation models, and repeatedly calculates the unit operation models in the plurality of regions in the reactor along the flow of the raw material ore or the flow of the combustion gas, until the difference between the calculated value of the temperature of the first gas phase and the first bed layer in a predetermined region and the previous calculated value in that region satisfies a predetermined convergence condition.

2. a first bed layer material quantity correcting unit that corrects the material quantity of the first bed layer by subtracting a movement amount of a second material, including at least one of a solid material and a liquid material, contained in the first bed layer from the movement amount of the second material to calculate a corrected first bed layer; 2. The apparatus for calculating a reactor reaction according to claim 1, wherein the first bed layer reaction amount calculation unit uses the modified first bed layer as the first bed layer.

3. a second bed layer material quantity correcting unit that corrects the material quantity of the second bed layer by subtracting the amount of a third substance, including at least one of a solid substance and a liquid substance, contained in the second bed layer from the amount of movement of the third substance to the first bed layer, and calculates a corrected second bed layer; 3. The apparatus for calculating a reactor reaction according to claim 1, wherein the second bed layer reaction amount calculation unit uses the modified second bed layer as the second bed layer.

4. a gas reaction amount calculation unit that calculates a gas reaction amount that contributes to an equilibrium state between the first gas phase and the first substance, and obtains a gas reaction amount that contributes to an equilibrium reaction of the first gas phase; a first substance reaction amount calculation unit that calculates a first substance reaction amount that contributes to the equilibrium state of the first substance with the first gas phase, and obtains a first substance reaction amount that contributes to the equilibrium reaction of the first substance; The reactor reaction calculation device according to any one of claims 1 to 3, comprising:

5. The unreacted amount of gas other than the reacted amount of gas generated in the gas reaction amount calculation unit is defined as the unreacted amount of gas, The unreacted amount of the first substance other than the reacted amount of the first substance generated in the first substance reaction amount calculation unit is defined as the unreacted amount of the first substance, a mixed gas phase calculation unit that calculates at least a flow rate of a mixed gas phase obtained by mixing the unreacted gas and a gas product generated by the reaction of the gas reactant with the first substance reactant in the gas region equilibrium reaction calculation unit; a first mixed substance calculation unit that calculates at least a flow rate of a first mixed substance obtained by mixing the unreacted portion of the first substance and a first substance product portion generated by the reaction of the gas reactant and the first substance reactant in the gas region equilibrium reaction calculation unit; The reactor reaction calculation device according to claim 4, comprising:

6. 6. The apparatus for calculating a reaction in a furnace according to claim 5, further comprising a first mixed substance quantity correction unit that corrects the quantity of the first mixed substance and calculates a first corrected mixed substance, taking into account the amount of movement of a fourth substance, which is contained in the first bed layer and includes at least one of a solid substance or a liquid substance, to the gas region and the amount of movement of a fifth substance, which is contained in the first mixed substance generated by the first mixed substance calculation unit and includes at least one of a solid substance or a liquid substance, to the first bed layer.

7. 7. The apparatus for calculating a reactor reaction according to claim 5, further comprising a first bed layer preliminary quantity correction unit that corrects the quantity of the first bed layer in consideration of the amount of movement of a fourth substance, which is contained in the first bed layer and includes at least one of a solid substance and a liquid substance, to the gas region, and the amount of movement of a fifth substance, which is contained in the first mixture generated by the first mixture calculation unit and includes at least one of a solid substance and a liquid substance, to the first bed layer.

8. 8. The apparatus for calculating a reaction in a furnace according to claim 1, further comprising an inflow gas combining unit that mixes a first inflow gas product generated by a reaction between the first bed layer reaction component and the first inflow gas in the first bed layer equilibrium reaction calculation unit and a second inflow gas product generated by a reaction between the second bed layer reaction component and the second inflow gas in the second bed layer equilibrium reaction calculation unit.

9. 9. The apparatus for calculating a reaction in a furnace according to claim 8, further comprising a gas quantity summing unit that mixes a mixed inflow gas generated by mixing the first inflow gas generation component and the second inflow gas generation component in the inflow gas summing unit with a gas generation component generated by the reaction of the first gas phase with the first substance in the gas region equilibrium reaction calculation unit.

10. The calculation device for a reactor reaction according to any one of claims 1 to 9, further comprising a bed layer distribution unit that distributes the first bed layer into two types of distribution bed layers in advance according to the reaction rate-causing conditions of the solid material and the liquid material contained in the first bed layer.

11. A distribution unit is provided which distributes the combustion material, which is introduced during the movement of the raw ore and includes at least one of a gas phase and the bed layer, to the added gas phase and the added bed layer at a mass flow rate, with the gas phase being an added gas phase and the bed layer being an added bed layer, the gas region equilibrium reaction calculation unit calculates an equilibrium reaction between the first substance and a second gas phase including the first gas phase and the additional gas phase; The first bed layer equilibrium reaction calculation unit calculates the equilibrium reaction between the inflow gas and a first combined distribution bed layer including the first bed layer and the additive bed layer.

12. A calculation device for furnace reactions according to any one of claims 1 to 11, wherein the raw ore is nickel oxide ore whose main components are iron and nickel.

13. A method for calculating a reaction in a reactor, in which raw material ore is supplied from one end of a reactor and moved toward the other end of the reactor, and the raw material ore is brought into countercurrent contact with combustion gas supplied from the other end of the reactor, thereby drying and reducing the raw material ore, a gas mixing amount calculation step of calculating a flow rate of an inflow gas that moves from the gas region in the reactor to a bed layer that is the raw material ore, out of the first gas phase that is the combustion gas flowing through the gas region in the reactor; an inlet gas distribution step of dividing the inlet gas into a first inlet gas and a second inlet gas; a first bed layer reaction amount calculation step of calculating, by using a first reaction rate equation, a first bed layer reaction amount contributing to an equilibrium state with the first inflow gas in one of two types of distributed bed layers distributed among the bed layers according to reaction rate-dependent conditions of the solid and liquid substances contained in the bed layer, thereby obtaining a first bed layer reaction amount contributing to the equilibrium reaction in the first bed layer; a second bed layer reaction amount calculation step of calculating a second bed layer reaction amount contributing to an equilibrium state with the second inlet gas of the other of the two types of distributed bed layers using a second reaction rate equation different from the first reaction rate equation to obtain a second bed layer reaction amount contributing to the equilibrium reaction of the second bed layer; a first bed layer equilibrium reaction calculation step of calculating an equilibrium reaction between the first bed layer reaction component and the first inflow gas; a second bed layer equilibrium reaction calculation step of calculating an equilibrium reaction between the second bed layer reaction component and the second inflow gas; a gas region equilibrium reaction calculation step of calculating an equilibrium reaction between the first gas phase and a first substance including at least one of a solid substance and a liquid substance present in the gas region based on a Gibbs energy minimization method; a first inflow gas reaction amount calculation step of calculating a gas reaction amount contributing to an equilibrium state of the first inflow gas with the first bed layer, thereby obtaining a first inflow gas reaction amount contributing to the equilibrium reaction of the first inflow gas; a first mixed inflow gas calculation step of calculating at least a flow rate of a first mixed inflow gas obtained by mixing a first inflow gas unreacted portion, which is an unreacted portion other than the first inflow gas reacted portion generated in the first inflow gas reaction amount calculation step, and a first inflow gas product portion generated by the reaction of the first bed layer reacted portion and the first inflow gas reacted portion in the first bed layer equilibrium reaction calculation step; a first mixed bed layer calculation step of calculating at least a flow rate of a first mixed bed layer obtained by mixing a first bed layer unreacted portion, which is an unreacted portion other than the first bed layer reacted portion generated in the first bed layer reaction amount calculation step, and a first bed layer product portion generated by the reaction of the first bed layer reacted portion with the first inflow gas reacted portion in the first bed layer equilibrium reaction calculation step; a second inflow gas reaction amount calculation step of calculating a gas reaction amount contributing to an equilibrium state of the second inflow gas with the second bed layer, thereby obtaining a second inflow gas reaction amount contributing to the equilibrium reaction of the second inflow gas; a second mixed inflow gas calculation step of calculating at least the flow rate of a second mixed inflow gas obtained by mixing a second inflow gas unreacted portion, which is an unreacted portion other than the second inflow gas reacted portion generated in the second inflow gas reaction amount calculation step, and a second inflow gas product portion generated by the reaction of the second bed layer reactant and the second inflow gas reactant in the second bed layer equilibrium reaction calculation step; a second mixed bed layer calculation step of calculating at least a flow rate of a second mixed bed layer obtained by mixing a second bed layer unreacted portion, which is an unreacted portion other than the second bed layer reacted portion generated in the second bed layer reaction amount calculation step, and a second bed layer product portion generated by the reaction of the second bed layer reacted portion with the second inflow gas reacted portion in the second bed layer equilibrium reaction calculation step; Including, the reaction rate-causing condition is both one or more properties of a solid substance contained in the bed layer selected from the group consisting of particle size, specific gravity, type, and segregation, and one or more properties of a liquid substance contained in the bed layer selected from the group consisting of specific gravity, type, and segregation, The method for calculating the reaction in the reactor is applied to each unit operation model when the reactor is divided into a plurality of regions and the reaction process in one region is assumed to be a unit operation model, the reaction process in the reactor is modeled by combining the unit operation models, and calculations of the unit operation models in the plurality of regions in the reactor are repeatedly performed along the flow of the raw material ore or the flow of the combustion gas, until the difference between the calculated value of the temperature of the first gas phase and the first bed layer in a predetermined region and the previous calculated value in that region satisfies a predetermined convergence condition.

Images