Calculation device for in-furnace reaction and calculation method for in-furnace reaction

Abstract

To provide a calculation device of an in-furnace reaction capable of calculating with high accuracy.SOLUTION: The calculation device of an in-furnace reaction in which raw material ores are brought into contact with combustion gas and dried and reduced, comprises: a first gas phase physical quantity correcting part for correcting a physical quantity of first gas phase flowing in a first gas region distributed according to a gas region affecting factors out of the gas region including combustion gas, and acquiring corrected first gas phase; a second gas phase physical quantity correcting part for correcting a physical quantity of second gas phase flowing in a second gas region and acquiring corrected second gas phase; a gas mixed amount calculation part for calculating a flow rate of inflow gas moving from the second gas region to a bed layer; a first gas region equilibrium reaction calculation part for calculating an equilibrium reaction of the corrected first gas phase with a first substance existing in the first gas region; a second gas region equilibrium reaction calculating part for calculating an equilibrium reaction of the corrected second gas phase with a second substance existing in the second gas region; and a bed layer equilibrium reaction calculating part for calculating an equilibrium reaction of the bed layer with the inflow gas.SELECTED DRAWING: Figure 2

Description

Technical Field

[0001] The present invention relates to a calculation device for in-furnace reaction and a calculation method for in-furnace reaction.

Background Art

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

[0003] In the dry smelting method using a rotary kiln, raw ore is dried in a rotary dryer so that the adhering moisture is, for example, 15% to 25%, and then the dried ore with reduced adhering water is charged from the charging end of the rotary kiln. Thereafter, the dried ore is heated by the combustion heat of coal supplied from the charging end of the rotary kiln or by a pulverized coal dedicated combustion burner or a pulverized coal and heavy oil co-combustion burner provided at the discharge end of the rotary kiln to dry the dried ore and perform firing.

[0004] As such a dry smelting method using a rotary kiln, for example, in addition to the combustion heat of coal supplied from the charging end and the combustion heat generated by the combustion of pulverized coal, heavy oil, etc. by a burner, there is a method of supplying the combustion heat generated by the combustion of coal charged from the middle of the rotary kiln to provide the heat necessary for drying and partial reduction of the dried ore (see, for example, Patent Document 1).

[0005] Patent Document 1 discloses an operating method of a rotary kiln in which coal is charged into the rotary kiln from a scoop feeder provided in the middle of the rotary kiln, and the dried ore of nickel oxide ore charged from the charging end of the rotary kiln is fired with the combustion heat generated by the combustion of fossil fuel by a burner and subjected to a partial reduction treatment. In this operating method of the rotary kiln, the coal charged from the scoop feeder is thermally decomposed into volatile matter and fixed carbon. The volatile matter is discharged from the charging end together with the combustion gas in the furnace, and the fixed carbon is discharged from the discharge end together with the sintered ore produced by drying and reducing the dried ore.

Prior Art Documents

Patent Documents

[0006]

Patent Document 1

Summary of the Invention

Problems to be Solved by the Invention

[0007] Here, in a well-known operating method such as the operating method of the rotary kiln described in Patent Document 1, the inside of the rotary kiln is treated as one type of gas phase, and the formation of a plurality of different types of gas phases affected by the bed layer has not been considered. Since the gas phase supplied into the rotary kiln is mixed with the generated gas generated from the raw ore, two or more layers with different gas components are generated, such as a layer that is easily affected by the bed layer and a layer that is hardly affected by the bed layer. Therefore, when simulating the behavior of raw ore in a reactor such as a rotary kiln, it is necessary to accurately calculate the in-furnace reaction in consideration of the generation of two or more phases with different gas components in the gas phase due to the generated gas generated from the raw ore.

[0008] One aspect of the present invention aims to provide a calculation device for in-furnace reactions that can calculate the reactions in a reactor with high precision.

Means for Solving the Problems

[0009] One aspect of the in-furnace reaction calculation device according to the present invention is an in-furnace reaction calculation device that moves raw ore supplied from one end side of a reactor toward the other end side, contacts the raw ore with combustion gas supplied from the other end side, and dries and reduces the raw ore, and in one of two gas regions distributed according to a gas region influence factor including product gas generated from a bed layer containing the raw ore among the gas regions containing the combustion gas, a first gas-phase material amount correction unit that corrects the material amount of a first gas phase flowing through a first gas region, and obtains a corrected first gas phase; a second gas-phase material amount correction unit that corrects the material amount of a second gas phase flowing through a second gas region, which is the other of the two gas regions, and obtains a corrected second gas phase; a gas mixing amount calculation unit that calculates the flow rate of inflow gas moving from the second gas region to the bed layer; a first gas region equilibrium reaction calculation unit that calculates an equilibrium reaction between the corrected first gas phase and a first substance present in the first gas region among substances including at least one of solid substances and liquid substances present in the gas region; a second gas region equilibrium reaction calculation unit that calculates an equilibrium reaction between the corrected second gas phase and a second substance present in the second gas region among the substances; a bed layer equilibrium reaction calculation unit that calculates an equilibrium reaction between the bed layer and the inflow gas; is provided.

[0010] One aspect of the in-furnace reaction calculation method according to the present invention is an in-furnace reaction calculation method that moves raw ore supplied from one end side of a reactor toward the other end side, contacts the raw ore with combustion gas supplied from the other end side, and dries and reduces the raw ore, and a first gas-phase material amount correction step of correcting the material amount of a first gas phase flowing through a first gas region, which is one of two gas regions distributed according to a gas region influence factor including product gas generated from a bed layer containing the raw ore among the gas regions containing the combustion gas, and obtaining a corrected first gas phase; a second gas-phase material amount correction step of correcting the material amount of a second gas phase flowing through a second gas region, which is the other of the two gas regions, and obtaining a corrected second gas phase; A gas mixing amount calculation step of calculating the flow rate of the inflow gas moving from the second gas region to the bed layer A first gas region equilibrium reaction calculation step of calculating an equilibrium reaction between the modified first gas phase and a first substance present in the first gas region among substances including at least one of a solid substance and a liquid substance present in the gas region A second gas region equilibrium reaction calculation step of calculating an equilibrium reaction between the modified second gas phase and a second substance present in the second gas region among the substances A bed layer equilibrium reaction calculation step of calculating an equilibrium reaction between the bed layer and the inflow gas and the like

Advantages of the Invention

[0011] One aspect of the in-furnace reaction calculation device according to the present invention can accurately calculate the reaction in the reactor

Brief Description of the Drawings

[0012]

Figure 1

Figure 2

Figure 3

Figure 4

Figure 5

Figure 6

Figure 7

Figure 8

Figure 9

Figure 10

Figure 11

Figure 12

Embodiments for Carrying Out the Invention

[0013] Hereinafter, embodiments of the present invention will be described in detail. For ease of understanding the description, the same reference numerals are given to the same components in each drawing, and duplicate descriptions are omitted. Also, the scales of the members in the drawings may be different from the actual ones.

[0014] When explaining the calculation device for in-furnace reaction according to the embodiment of the present invention, the configuration of the rotary kiln to which the calculation device for in-furnace reaction according to the present embodiment is applied will be explained.

[0015] <Rotary Kiln> FIG. 1 shows a schematic configuration of a rotary kiln to which the calculation device for in-furnace reaction according to the present embodiment is applied. As shown in FIG. 1, the rotary kiln 1 has a kiln body 11 that is rotatable and substantially cylindrical, and a combustion material supply pipe 12 provided in the middle of the kiln body 11.

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

[0017] As for the size of the kiln body 11, for example, it is preferable to use one with an inner diameter of 4.5 m to 5.5 m and a major axis direction length (total length) of 100 m to 110 m.

[0018] One end (the left side in FIG. 1) of the kiln main body 11 has an opening end 11a inserted and closed into the rotary kiln charging end (hereinafter also simply referred to as the "charging end") 14A, and the opening end 11b at the other end (the right side in FIG. 1) is inserted and closed into the rotary kiln discharging end (hereinafter also referred to as the "discharging end") 14B. The kiln main body 11 is disposed in a slightly inclined state from the charging end 14A toward the discharging end 14B and is rotatably supported about an axis.

[0019] A raw material supply pipe 15 for introducing raw ore into the kiln main body 11 is provided through the charging end 14A. A burner 16 is provided at the discharging end 14B, which penetrates the opening end 11b and is introduced into the kiln main body 11.

[0020] As the raw ore, nickel oxide ore (nickel oxide ore) or the like can be used. For example, the raw ore can be pre-dried by a dryer (rotary dryer) such as nickel oxide ore or the like to remove a part of the adhering moisture, and dried ore or the like can be used. The moisture content in the dried ore is about 15% to 25% by mass.

[0021] The nickel oxide ore as the raw ore is not particularly limited, but in the smelting of ferronickel, which is an alloy mainly composed of iron and nickel, garnierite ore or the like is preferably used. As a typical composition of garnierite ore, in terms of dry ore conversion, the Ni grade is 2.1% to 2.5% by mass, the Fe grade is 11% to 23% by mass, the MgO grade is 20% to 28% by mass, the SiO2 grade is 29% to 39% by mass, the CaO grade is less than 0.5% by mass, and the ignition loss is 10% to 15% by mass.

[0022] As the burner 16, a pulverized coal dedicated combustion burner or a mixed combustion burner of pulverized coal and heavy oil can be used. The burner 16 burns a fuel containing pulverized coal or pulverized coal and heavy oil or the like to generate combustion heat in the rotary kiln 1.

[0023] The combustion material supply pipe 12 is provided in the middle of the outer peripheral surface of the kiln body 11 and can supply combustion materials into the kiln body 11. The combustion materials include at least one of substances mainly entering the gas phase such as volatile components and solid substances mainly entering the bed layer such as fixed carbon, and for example, coal materials such as coal can be used.

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

[0025] Fixed carbon is the combustible component mainly composed of carbon excluding ash among char particles (mainly fixed carbon and ash) which are the residues after thermal decomposition after moisture and volatile components have escaped from coal.

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

[0027] The combustion materials are often not single-species coal materials, but a mixture of multiple different species of coal materials, etc., and furthermore, the particle sizes of the coal materials to be input often have a distribution. Note that the particle size refers to the volume-average particle size based on the effective diameter, and the particle size is measured by, for example, the laser diffraction / scattering method, the dynamic light scattering method, or the classification method. When using the laser diffraction / scattering method, the particle size (D 50 ) at which the cumulative particle size distribution from the small particle size side becomes 50% in the volume-based particle size distribution measured by the laser diffraction / scattering method can be used as the average particle size.

[0028] The raw ore is charged into the kiln body 11 from the raw material supply pipe 15 provided at the charging end 14A, and the combustion materials are input into the kiln body 11 from the combustion material supply pipe 12. From the discharge end 14B side, high-temperature combustion gas generated by burning pulverized coal, heavy oil, etc. by the burner 16 installed at the discharge end 14B is blown from the discharge end 14B side toward the charging end 14A side, that is, in the direction opposite to the flow direction of the raw ore.

[0029] Inside the kiln body 11, the raw ore is charged from the charging end 14A. As the kiln body 11 rotates at a predetermined speed, the raw ore charged into the kiln body 11 from the charging end 14A through the raw material supply pipe 15 is conveyed from the opening end 11a side, which is one end side, toward the discharge end 14B, which is the other end side. At this time, while the raw ore moves inside the kiln body 11, it comes into countercurrent contact with the combustion gas flowing from the discharge end 14B toward the charging end 14A side, and is heated by the combustion heat and flame of the high-temperature combustion gas generated by burning fuels such as pulverized coal and heavy oil with the burner 16. Also, the combustion material introduced into the kiln body 11 from the combustion material supply pipe 12 burns due to the combustion gas inside the kiln body 11. The raw ore is also heated by the combustion heat generated by the combustion of the combustion material introduced from the combustion material supply pipe 12. Therefore, as the raw ore moves from the charging end 14A to the discharge end 14B of the kiln body 11 along with the rotation of the kiln body 11, it is heated by the combustion heat and flame of the combustion gas generated by the combustion of the fuel with the burner 16 and the combustion heat generated by the combustion of the combustion material, and the temperature gradually rises.

[0030] Inside the kiln body 11, between the raw ore and the combustion gas, there is a transfer of substances due to the evaporation of moisture contained in the raw ore and the combustion material, the volatilization of volatile components contained in the combustion material, aggregation, and the scattering and falling of ash contained in the burner fuel and the combustion material. The moisture and volatile components generated by the thermal decomposition of the combustion material supplied from the middle of the rotary kiln 1 move together with the combustion gas to the charging end 14A side, and the char particles move to the discharge end 14B together with the raw ore.

[0031] Before the raw ore inside the kiln body 11 reaches the discharge end 14B, the moisture contained in the raw ore is almost completely removed, and it is fired and partially reduced to become sinter. The sinter is discharged from the discharge end 14B.

[0032] The sinter consists of, for example, a temperature of 800 - 900°C and a particle size of about 10 mm - 100 mm.

[0033] At the discharge port of the discharge end 14B, a loster (screening device) 17 is provided to separate sintered ore with a particle size of about 10 mm to 100 mm and sintered lumps (with a particle size of about 100 mm to 500 mm) generated in the rotary kiln 1. The loster 17 is composed of, for example, an iron lattice with an opening of about 100 mm. The sintered ore discharged from the discharge end 14B passes through the loster 17 and then is conveyed to the next process through the sintered ore discharge chute 18.

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

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

[0036] The gas region is the region where the gas phase in the rotary kiln 1 flows, and is distributed to the first gas region and the second gas region so as to have a predetermined material quantity ratio in each region divided according to the gas region influence factor including the generated gas and the like generated from the bed layer. The gas region influence factor refers to the positions of fuel supply facilities such as the raw material supply pipe 15 and the combustion material supply pipe 12, the positions of gas supply facilities such as the burner 16, the generated gas generated in the bed layer, dust, and the in-furnace gas concentration. The first gas region refers to the gas region with little influence of the gas region influence factor. The second gas region refers to the gas region with a large influence of the gas region influence factor.

[0037] The gas phase supplied into the rotary kiln 1 is described as the gas phase G1. The "gas phase G1" is the combustion gas flowing from one adjacent region (region (A + 1)) into region A. The combustion gas includes gases contributing to the equilibrium reaction such as volatile components, oxygen, carbon dioxide, carbon monoxide, hydrogen, etc., and may include inert substances (e.g., nitrogen, etc.) that do not contribute to the reaction in the first gas region equilibrium reaction calculation unit 214-1 and the second gas region equilibrium reaction calculation unit 214-2 described later.

[0038] The gas phase existing in region A is described as the first gas phase G11, the first gas phase containing additives G12, the modified first gas phase G13, the first mixed gas phase G14, the second gas phase G21, the modified second gas phase G22A, the adjusted second gas phase G22B, the second mixed gas phase G23, or the second total gas phase G24. Each gas phase in region A is defined as follows. The "first gas phase G11" refers to the gas phase flowing in the first gas region among the gas phase G1 flowing from one adjacent region (region (A + 1)) into region A. The first gas phase G11 may include inert substances and the like that do not contribute to the reaction in the first gas region equilibrium reaction calculation unit 214-1 described later, similar to the gas phase G1. The "first gas phase containing additives G12" is the gas phase obtained by combining the additive gas phase distributed from the combustion material and the first gas phase. The "first diffused gas G121" is the gas in the first gas phase containing additives G12 existing in the first gas region that moves to the second gas region by convection or diffusion. The "second diffused gas G122" is the gas in the second gas phase G21 existing in the second gas region that moves to the first gas region by convection or diffusion. The "modified first gas phase G13" is the first gas phase in which the amount of substance of the first gas phase containing additives G12 is modified in consideration of the amount of gas in the first gas phase containing additives G12 moving to the second gas region and the amount of gas in the second gas phase G21 existing in the second gas region moving to the first gas phase containing additives G12. The "first mixed gas phase G14" is the first gas phase generated in the first mixed gas phase calculation unit 216-1 described later. The "second gas phase G21" refers to the gas phase of the gas phase G1 flowing from one adjacent region (region (A + 1)) into region A that passes through the second gas region. Similar to the gas phase G1, the second gas phase G21 may contain inert substances and the like that do not contribute to the reaction in the second gas region equilibrium reaction calculation unit 214-2 described later. The "modified second gas phase G22A" is the second gas phase obtained by correcting the amount of substance of the additive-containing first gas phase G12 by considering the amount of the second substance in the additive-containing first gas phase G12 existing in the first gas region that moves to the second gas phase G21 and the amount of the third substance contained in the second gas phase G21 that moves to the additive-containing first gas phase G12. The "adjusted second gas phase G22B" is the second gas phase obtained by removing the flow rate of the inflowing gas that has moved from the modified second gas phase G22A to the bed layer in the gas mixing amount calculation unit 204 described later. The "second mixed gas phase G23" is the second gas phase generated in the second mixed gas phase calculation unit 216-2 described later. The "second total gas phase G24" is the gas phase generated in the gas amount total calculation unit 222 described later.

[0039] The bed layer supplied into the rotary kiln 1 is described as the bed layer S1. The "bed layer S1" refers to the raw ore supplied from the charging end 14A in the kiln main body 11. The raw ore is a solid substance containing ore compounds such as nickel oxide, iron oxide, and magnesium oxide, fixed carbon, etc., and a liquid substance, and may contain inert substances and the like that do not contribute to the reaction in the first bed layer equilibrium reaction calculation unit 215-1 and the second bed layer equilibrium reaction calculation unit 215-2 described later.

[0040] The bed layer existing in region A is described as the first bed layer S11, the additive-containing first bed layer S12, the preliminary modified first bed layer S13A, the modified first bed layer S13B, the first mixed bed layer S14, the second bed layer S21, the preliminary modified second bed layer S21A, the modified second bed layer S21B, or the second mixed bed layer S22. Each bed layer in region A is defined as follows. The "first bed layer S11" refers to the bed layer that exists in the part of the bed layer S1 flowing into region A from one adjacent region (region (A-1)) and has good contact with the gas region on the bed layer surface. Similar to the bed layer S1, the first bed layer S11 may contain inert substances and the like that do not contribute to the reaction in the first bed layer equilibrium reaction calculation unit 215-1 described later. The "first bed layer S12 containing additives" is a bed layer obtained by adding up the added bed layer distributed from the combustion material and the first bed layer S11. The "preliminary corrected first bed layer S13A" refers to the first bed layer in which the amount of the first bed layer S12 containing additives is corrected in consideration of the amount of the seventh substance containing at least one of a solid substance or a liquid substance moving into the gas region in the first bed layer S12 containing additives and the amount of the eighth substance containing at least one of a solid substance or a liquid substance existing in the gas region and moving into the first bed layer S12 containing additives in the second mixture primary amount correction unit 223 described later and contained in the second mixture M22. The "corrected first bed layer S13B" refers to the first bed layer in which the amount of the preliminary corrected first bed layer S13A is corrected in consideration of the amount of the sixth substance containing at least one of a solid substance or a liquid substance in the preliminary corrected first bed layer S13A moving into the preliminary corrected second bed layer S21A and the amount of the fifth substance containing at least one of a solid substance or a liquid substance in the preliminary corrected second bed layer S21A moving into the preliminary corrected first bed layer S13A. The "first mixed bed layer S14" refers to the bed layer generated in the first mixed bed layer calculation unit 218-1 described later. The "second bed layer S21" refers to the bed layer other than the first bed layer S11 among the bed layers S1 flowing into region A from one adjacent region (region A-1), and is the bed layer that exists in the part of the bed layer S1 that is less affected by the gas region inside the bed layer. Similar to the bed layer S1, the second bed layer S21 may contain inert substances and the like that do not contribute to the reaction in the second bed layer equilibrium reaction calculation unit 215-2 described later. The "preliminary corrected second bed layer S21A" is obtained by correcting the amount of material in the second bed layer S21 excluding the amount of the ninth substance containing at least one of the solid component or the liquid component that has moved into the gas region of the second bed layer S21. The "corrected second bed layer S21B" is the second bed layer obtained by correcting the amount of material in the preliminary corrected second bed layer S21A in consideration of the amount of the sixth substance in the preliminary corrected first bed layer S13A that has moved into the preliminary corrected second bed layer S21A and the amount of the fifth substance in the preliminary corrected second bed layer S21A that has moved into the preliminary corrected first bed layer S13A. The "second mixed bed layer S22" is the bed layer generated by the second mixed bed layer calculation unit 218-2 described later.

[0041] Also, the gas flowing from the gas region in region A into the bed layer is described as the inflow gas G221, the first inflow gas G221-1, the second inflow gas G221-2, the first mixed inflow gas G222A, the second mixed inflow gas G222B, or the mixed inflow gas G223. In the following description, each inflow gas in region A is defined as follows. The "inflow gas G221" is the gas flowing from the gas region in region A into the bed layer. The "first inflow gas G221-1" is the inflow gas among the inflow gas G221 that has flowed into the bed layer and flows into the first bed layer S11. The "second inflow gas G221-2" is the inflow gas among the inflow gas G221 that has flowed into the bed layer and flows into the second bed layer S21. The "first mixed inflow gas G222A" is the gas obtained by mixing the unreacted portion of the first bed layer and the first inflow gas generated component generated by the first bed layer equilibrium reaction calculation unit 215-1. The "second mixed inflow gas G222B" is the gas obtained by mixing the unreacted portion of the corrected second bed layer and the second inflow gas generated component generated by the second bed layer equilibrium reaction calculation unit 215-2. The "mixed inflow gas G223" is the gas obtained by adding up the first mixed inflow gas G222A and the second mixed inflow gas G222B.

[0042] The substance supplied into the rotary kiln 1 is described as substance M1. "Substance M1" includes at least one of a solid substance or a liquid substance existing in the gas region.

[0043] The substances existing in region A are described as the first substance M11, the second substance M21, the third substance M3, the fourth substance M4, the fifth substance M5, the sixth substance M6, the seventh substance M7, the eighth substance M8, the ninth substance M9, the first mixed substance M12, the first modified mixed substance M13, the second mixed substance M22, the second modified mixed substance M23, or the second modified mixed substance M24. In the following description, each substance in region A is defined as follows. "The first substance M11" is the substance among the substances M1 flowing into the gas region that flows into the first gas region. "The second substance M21" is the substance among the substances M1 flowing into the gas region that flows into the second gas region. "The third substance M3" is a substance containing at least one of a solid substance or a liquid substance, which is contained in the second modified mixed substance M23 existing in the second gas region. "The fourth substance M4" is a substance containing at least one of a solid substance or a liquid substance, which is contained in the first mixed substance M12 existing in the first gas region. "The fifth substance M5" is one containing at least one of a solid substance or a liquid substance, which is contained in the second mixed bed layer S21A existing in the second bed layer. "The sixth substance M6" is one containing at least one of a solid substance or a liquid substance, which is contained in the preliminary modified first bed layer S13A existing in the first bed layer. "The seventh substance M7" is one containing at least one of a solid substance or a liquid substance, which is contained in the additive-containing first bed layer S12 existing in the first bed layer. "The eighth substance M8" is a substance containing at least one of a solid substance or a liquid substance, which is contained in the second mixed substance M22 existing in the second gas region. "The ninth substance M9" is one containing at least one of a solid substance or a liquid substance, which is contained in the second bed layer S21 existing in the second bed layer. "The "first mixture substance M12" is a mixture of the first substance M1 generated by the first mixture substance calculation unit 217-1 described later." "The "first corrected mixture substance M13" is a mixture of the first substance M1 generated by the first mixture substance amount correction unit 220 described later." "The "second mixture substance M22" is the second mixture substance M21 obtained by mixing the unreacted portion Mr21 of the second substance and the generated component MP21 of the second substance in the second mixture substance calculation unit 217-2 described later." "The "second corrected mixture substance M23" is obtained by adding the amount of movement of the seventh substance M7 moving from the first bed layer and the amount of movement of the ninth substance M9 moving from the second bed layer to the second mixture substance M22, excluding the amount of movement of the eighth substance M8 contained in the second mixture substance M22, in the second mixture substance first amount correction unit 223 described later, and correcting the amount of the second mixture substance M22." "The "second corrected mixture substance M24" is obtained by adding the amount of movement of the fourth substance M4 contained in the first mixture substance M12 to the second corrected mixture substance M23 and correcting the amount of the second corrected mixture substance M23, excluding the third substance M3 moving from the second corrected mixture substance M23 to the first gas region, in the second mixture substance second amount correction unit 224 described later."

[0044] In addition, the reacted portion, unreacted portion, and generated component calculated within region A are described as the first gas reacted portion GR1, the first gas unreacted portion Gr1, the first gas generated component GP1, the second gas reacted portion GR2, the second gas unreacted portion Gr2, the second gas generated component GP2, the first substance reacted portion MR11, the first substance unreacted portion Mr11, the first substance generated component MP11, the second substance reacted portion MR21, the second substance unreacted portion Mr21, the second substance generated component MP21, the first inflow gas reacted portion GR221, the first inflow gas unreacted portion Gr221, the first inflow gas generated component GP221, the second inflow gas reacted portion GR222, the second inflow gas unreacted portion Gr222, the second inflow gas generated component GP222, the first bed layer reacted portion SR1, the first bed layer unreacted portion Sr1, the first bed layer generated component SP1, the second bed layer reacted portion SR2, the second bed layer unreacted portion Sr2, or the second bed layer generated component SP2. These are calculated by any of the respective parts constituting the in-furnace reaction calculation device shown in FIG. 2."

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

[0046] In this embodiment, the bed layer equilibrium reaction calculation unit 215 is composed of two bed layers (a first bed layer equilibrium reaction calculation unit 215-1 and a second bed layer equilibrium reaction calculation unit 215-2).

[0047] Further, the in-furnace reaction calculation device 20 can include a gas region distribution unit 231 that distributes the gas phase G1 into two types of distributed gas phases (a first gas phase G11 and a second gas phase G21) in advance according to the gas region influence factor including the generated gas generated from the bed layer. In this embodiment, the in-furnace reaction calculation device 20 supplies the first gas phase G11 generated by the gas region distribution unit 231 to the first gas-phase material quantity correction unit 202, and supplies the second gas phase G21 generated by the gas region distribution unit 231 to the second gas-phase material quantity correction unit 203.

[0048] The in-furnace reaction calculation device 20 can include a first substance distribution unit 232 that distributes the first substance M1 in advance to two types of first substances (the first substance M11 and the second substance M21). In the present embodiment, the in-furnace reaction calculation device 20 supplies the first substance M11 generated in the first substance distribution unit 232 to the first substance reaction amount calculation unit 211-1, and supplies the second substance M21 generated in the first substance distribution unit 232 to the second substance reaction amount calculation unit 211-2.

[0049] The in-furnace reaction calculation device 20 can include a bed layer distribution unit 233 that distributes the bed layer S1 in advance to two types of distribution bed layers (the first bed layer S11 and the second bed layer S21) according to the reaction rate causative conditions of the solid substances and liquid substances contained in the bed layer S1. In the present embodiment, the in-furnace reaction calculation device 20 supplies the first bed layer S11 generated in the bed layer distribution unit 233 to the first bed layer preliminary material amount correction unit 205, and supplies the second bed layer S21 generated in the bed layer distribution unit 233 to the second bed layer preliminary material amount correction unit 207.

[0050] In the rotary kiln 1, reaction gases after reaction such as water vapor and gas (product gas) are generated from the bed layer due to thermal decomposition and reaction of the raw materials. For example, in a rotary kiln that bakes nickel oxide ore dry ore and performs a partial reduction treatment, water vapor from the bed layer due to dehydration of the adhering water and crystal water remaining in the dry ore, or CO, CO2 gas, etc. generated by combustion of coal input from the middle of the rotary kiln and partial reduction of the ore. When these reaction gases after reaction are released from the bed layer to the gas region, they diffuse throughout the gas region by convection and diffusion, but a certain residence time of the in-furnace gas is required for the concentration in the gas region to become uniform. Usually, in the gas region inside the rotary kiln 1, it is considered that there is a distribution in gas concentration between the vicinity of the bed layer and other regions.

[0051] For example, as shown in FIG. 3, in the bed layer, a first bed layer S11 that is less affected by gas and a second bed layer S21 that has good contact with the gas region are likely to occur. As shown in FIG. 4, along the longitudinal direction of the rotary kiln 1, the first bed layer S11 and the second bed layer S21 are formed. Also, in the gas region, after the reaction gas and dust generated in the bed layer are scattered, in the gas region, a first gas phase G11 that contains little product gas, dust, etc. and a second gas phase G21 that contains a lot of product gas, dust, etc. are likely to occur. As shown in FIG. 4, along the longitudinal direction of the rotary kiln 1, the first gas phase G11 and the second gas phase G21 are formed. Therefore, the reaction in the gas region is easily affected by the product gas, dust, etc. contained in the gas phase.

[0052] Also, the reaction in the bed layer is greatly affected by the gas composition. For example, in the case of the Fe-C-O system and when the CO gas partial pressure is low, the reaction product is Fe3O4, but as the CO partial pressure increases, the product changes to FeO or M-Fe. Therefore, in order to simulate the in-furnace reaction of the rotary kiln 1 and accurately grasp the reaction behavior of the raw materials in the rotary kiln 1, it is necessary to express the gas concentration (gas partial pressure) used in the equilibrium reaction calculation of the bed layer in more detail.

[0053] As methods for simulating a rotary kiln, roughly classified, there are a method in which it is assumed that the gas generated in the ore layer spreads instantaneously throughout the gas region and the gas components in the gas region at the same position in the length direction are treated with a single concentration, and a method in which the calculation of the gas phase is performed by three-dimensional thermal fluid calculation and coupled calculation. In the former method, the gas atmosphere with a high product gas concentration from the bed layer in the gas region near the bed layer cannot be expressed. In the latter method, enormous calculation resources and time are required.

[0054] The in-furnace reaction calculation device 20 divides the gas region into two gas regions (the first gas region and the second gas region) in the rotary kiln 1, and performs equilibrium reaction calculations for the respective bed layers of the first gas phase G11 flowing through the first gas region and the second gas phase G21 flowing through the second gas region. That is, when performing a simulation of the in-furnace reaction based on the Gibbs energy minimization method, the in-furnace reaction calculation device 20 takes into account the post-reaction gas, dust, etc. from the bed layer contained in the gas region in the rotary kiln 1. Then, the in-furnace reaction calculation device 20 independently calculates the reaction amounts required for the equilibrium reaction for the first gas phase G11, which contains almost no post-reaction gas, dust, etc., and the second gas phase G21, which contains a large amount of generated gas, dust, etc., in the gas phase G1. Thereby, the in-furnace reaction calculation device 20 can accurately calculate the equilibrium reaction based on the Gibbs energy minimization method in consideration of the influence of the post-reaction gas, dust, etc. contained in the first gas phase G11 and the second gas phase G21. Since the in-furnace reaction calculation device 20 can represent the gas atmosphere with a high concentration of generated gas from the bed layer in the gas region near the bed layer and does not require a huge amount of computing resources, the in-furnace reaction calculation device 20 can analyze the in-furnace reaction with high accuracy without imposing a large burden.

[0055] In addition, the in-furnace reaction calculation device 20 performs equilibrium reaction calculations for the bed layer and the inflowing gas mixed in the bed layer based on the Gibbs energy minimization method. Thereby, the equilibrium reaction in the bed layer in the rotary kiln 1 can be calculated.

[0056] Therefore, by considering that the gas phase in the rotary kiln 1 exhibits two different reaction behaviors according to the generated gas, dust, etc. from the bed layer contained in the gas phase, the in-furnace reaction calculation device 20 can accurately analyze the reactions occurring in the rotary kiln 1 and improve the calculation accuracy of the reactions. Therefore, the in-furnace reaction calculation device 20 can calculate the reactions in the rotary kiln 1 with higher accuracy.

[0057] As shown in FIG. 2, the in-furnace reaction calculation device 20 inputs the material quantity, temperature, pressure, etc. of the first gas phase G11 and the second gas phase G21 flowing into the gas region of region A and the first bed layer S11 and the second bed layer S21 flowing into the bed layer of region A. The in-furnace reaction calculation device 20 uses these input values to perform material transfer and reaction calculations in each divided region A.

[0058] The first gas phase G11 and the second gas phase G21 are determined according to the material quantity, etc. when the gas phase G1 supplied from the discharge end 14B in the kiln main body 11 is distributed in the region N including the discharge end 14B and moves to region A.

[0059] The first bed layer S11 and the second bed layer S21 are determined according to the material quantity, etc. when the bed layer S1 supplied from the charging end 14A in the kiln main body 11 is distributed in the region 1 including the charging end 14A and moves to region A.

[0060] When performing the material transfer and reaction calculations in region A, the in-furnace reaction calculation device 20 defines the amount of gas (inflow gas G221) mixed into the bed layer among the second gas phase G21 flowing in the gas region. As a method for defining this gas mixing amount, for example, there is a method of giving a fixed value for the ratio of the amount mixed into the bed layer in the total gas amount.

[0061] Further, the in-furnace reaction calculation device 20 defines the amount of material transfer from the solid or liquid material moving in the bed layer to the gas region due to scattering or volatilization and the amount of material transfer from the solid or liquid material moving in the gas region to the bed layer due to falling (sedimentation) or aggregation. As a method for defining this material transfer amount, for example, there is a method of defining for each calculated material based on the sedimentation rate calculated from the specific gravity, particle size, etc. of the solid or liquid material.

[0062] The substances containing at least one of the solid and liquid substances generated after the reaction are combined with the unreacted components in their respective bed layers (the first bed layer and the second bed layer) and used as the input values for the respective bed layers in the next adjacent region. Further, the post-reaction gas generated in each bed layer is combined with the unreacted gas, and further combined with the post-reaction gas of each bed layer and the second gas region to serve as the input value for the adjacent gas region.

[0063] As shown in FIG. 2, the distribution unit 201 distributes the combustion material 240 in terms of mass flow rate to the gas phase and the bed layer contained in the combustion material 240 by presetting the type, components, etc. of the combustion material 240.

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

[0065] The added component-containing first gas phase G12, which includes the first gas phase G11, which is combustion gas flowing from one adjacent region (region (A + 1)) into the first gas region of region A in the furnace, and the added gas phase AG11 of the combustion material distributed by the distribution unit 201, flows in the first gas region.

[0066] The material quantity correction unit 202 of the first gas phase corrects the material quantity of the added component-containing first gas phase G12 flowing through the first gas region, and calculates and obtains the corrected first gas phase G13.

[0067] The material quantity correction unit 203 of the second gas phase corrects the material quantity of the second gas phase G21 flowing through the second gas region, and calculates and obtains the corrected second gas phase G22A.

[0068] The gas mixture quantity calculation unit 204 calculates the mass flow rate of the inflow gas G221 moving from the corrected second gas phase G22A in the second gas region to the bed layer.

[0069] The gas mixture quantity calculation unit 204 calculates a gas phase obtained by subtracting the inflow gas G221 from the corrected second gas phase G22A as the adjusted second gas phase G22B.

[0070] The first bed layer preliminary quantity correction unit 205 corrects the quantity of the first bed layer S12 containing additives, which is a mixture of the first bed layer S11 and the additional bed layer AS11, to calculate a preliminarily corrected first bed layer S13A.

[0071] The quantity of the first bed layer S12 containing additives is corrected considering the transfer amount to the second mixed substance M22 existing in the gas region of the seventh substance M7 containing at least one of a solid substance and a liquid substance in the first bed layer S12 containing additives, and the transfer amount of the eighth substance M8 containing at least one of a solid substance and a liquid substance in the second mixed substance M22 to the first bed layer S12 containing additives.

[0072] The first bed layer quantity correction unit 206 corrects the quantity of the preliminarily corrected first bed layer S13A considering the transfer amount of the sixth substance M6 containing at least one of a solid substance and a liquid substance in the preliminarily corrected first bed layer S13A to the second bed layer S21, and the transfer amount of the fifth substance M5 containing at least one of a solid substance and a liquid substance in the second bed layer S21 to the preliminarily corrected first bed layer S13A, and calculates a corrected first bed layer S13B.

[0073] In this embodiment, the first bed layer quantity correction unit 206 corrects the quantity of the preliminarily corrected first bed layer S13A by adding the transfer amount of the fifth substance M5 contained in the second bed layer S21 to the preliminarily corrected first bed layer S13A, excluding the transfer amount of the sixth substance M6 to the second bed layer S21, and calculates a corrected first bed layer S13B.

[0074] The second bed layer preliminary quantity correction unit 207 corrects the quantity of the second bed layer S21 considering the transfer amount to the second mixed substance M22 existing in the gas region of the ninth substance M9 containing at least one of a solid substance and a liquid substance in the second bed layer S21, and calculates a preliminarily corrected second bed layer S21A.

[0075] The second bed layer material quantity correction unit 208 corrects the material quantity of the pre-corrected second bed layer S21A in consideration of the amount of the sixth substance M6 in the pre-corrected 21st bed layer S13A moving to the second bed layer S21 and the amount of the fifth substance M5 in the pre-corrected second bed layer S21A moving to the pre-corrected first bed layer S13A, and calculates the corrected second bed layer S21B.

[0076] In this embodiment, the second bed layer material quantity correction unit 208 corrects the material quantity of the second bed layer S21 by adding the amount of the sixth substance M6 in the pre-corrected 21st bed layer S13A moving to the second bed layer S21 and subtracting the amount of the fifth substance M5 in the pre-corrected second bed layer S21A moving to the pre-corrected first bed layer S13A, and calculates the corrected second bed layer S21B.

[0077] The inflow gas distribution unit 209 has a function of distributing the inflow gas G221 by mass flow rate into a first inflow gas G221-1 and a second inflow gas G221-2 by presetting the type, components, etc. of the inflow gas G221 calculated by the gas mixing amount calculation unit 204.

[0078] The first gas reaction amount calculation unit 210-1 calculates the first gas reaction amount that contributes to the equilibrium state with the first substance M11 in the corrected first gas phase G13. That is, the first gas reaction amount calculation unit 210-1 calculates the first gas reaction amount that is expected to undergo an equilibrium reaction between the corrected first gas phase G13 and the first substance M11 during the residence time of the corrected first gas phase G13 existing in the first gas region in region A. The reaction amount at the residence time of the corrected first gas phase G13 in region A is calculated using the reaction rate parameters set for each individual substance and the residence time within region A.

[0079] The first gas reaction amount calculation unit 210-1 can calculate the first gas reaction amount that reacts at the residence time within region A using the reaction rate formula.

[0080] The first gas reaction amount calculation unit 210-1 obtains, from the first gas reaction amount, the first gas reaction component GR1 that contributes to the equilibrium state with the first substance M11 in the corrected first gas phase G13.

[0081] In this specification, the reaction component refers to the flow rate that contributes to the reaction due to contact with the modified first gas phase G13 or the first substance M11 existing in the gas region during the residence time when the modified first gas phase G13 or the first substance M11 in region A passes through. Hereinafter, other reaction components also represent the meaning of the flow rate that contributes to the reaction due to contact between any of the gas phase, solid phase, and substance and any of the gas phase, solid phase, and substance.

[0082] That is, the first gas reaction component GR1 is the reaction component that contributes to the equilibrium reaction with the first substance M11 existing in the first gas region among the modified first gas phase G13. The first substance reaction component MR11 is the reaction component that contributes to the equilibrium state with the modified first gas phase G13 among the first substances M11 existing in the first gas region.

[0083] The first gas reaction amount calculation unit 210-1 can obtain the amount that reaches equilibrium between the modified first gas phase G13 and the solid or liquid substance existing in the first substance M11 by adopting, for example, a model that can be calculated from the reaction rate according to the Arrhenius-type reaction rate equation and the residence time of the modified first gas phase G13 in region A. Then, the first gas reaction amount calculation unit 210-1 divides it into the first gas reaction component GR1 corresponding to the gas reaction amount and the first gas unreacted component Gr1 corresponding to the remaining unreacted mass flow rate.

[0084] In the rotary kiln 1, the flow of the modified first gas phase G13 and the first substance M11 existing in the first gas region is in parallel flow. The modified first gas phase G13 has a high flow rate and a short residence time in region A, and the modified first gas phase G13 and the first substance M11 existing in the first gas region do not reach an equilibrium state. In this embodiment, the first gas reaction amount calculation unit 210-1 assumes that only a part of the modified first gas phase G13 in a certain region (region A) in the rotary kiln 1 reaches an equilibrium state with the first substance M11 existing in the first gas region, and divides the modified first gas phase G13 into the gas reaction component GR1 and the gas unreacted component Gr1 other than the gas reaction component GR1 according to the corresponding first gas reaction amount.

[0085] The second gas reaction amount calculation unit 210-2 calculates the second gas reaction amount that contributes to the equilibrium state with the second substance M21 in the adjusted second gas phase G22B obtained by subtracting the inflow gas G221 from the corrected second gas phase G22A. That is, the second gas reaction amount calculation unit 210-2 calculates the second gas reaction amount that is expected to undergo an equilibrium reaction between the adjusted second gas phase G22B and the second substance M21 during the residence time of the adjusted second gas phase G22B present in the first gas region in region A. The reaction amount at the residence time of the adjusted second gas phase G22B in region A is calculated using the reaction rate parameters set for each individual substance and the residence time within region A.

[0086] The second gas reaction amount calculation unit 210-2 can calculate the second gas reaction amount that reacts at the residence time within region A using a reaction rate formula.

[0087] The second gas reaction amount calculation unit 210-2 obtains, from the second gas reaction amount, the second gas reaction component GR2 that contributes to the equilibrium state with the second substance M21 in the adjusted second gas phase G22B.

[0088] That is, the first gas reaction component GR2 is the reaction component that contributes to the equilibrium reaction between the adjusted second gas phase G22B and the second substance M21 present in the second gas region. The second substance reaction component MR21 is the reaction component that contributes to the equilibrium state between the second substance M21 present in the second gas region and the adjusted second gas phase G22B.

[0089] The second gas reaction amount calculation unit 210-2 can obtain the amount that reaches the equilibrium state between the adjusted second gas phase G22B and the solid or liquid substance present in the second substance M21 by adopting, for example, a model that can be calculated from the reaction rate according to the Arrhenius-type reaction rate formula and the residence time of the adjusted second gas phase G22B within region A. Then, the second gas reaction amount calculation unit 210-2 divides it into the second gas reaction component GR2 corresponding to the gas reaction amount and the second gas unreacted component Gr2 corresponding to the remaining unreacted mass flow rate.

[0090] In the rotary kiln 1, the adjusted second gas phase G22B and the flow of the second substance M21 present in the second gas region are in co-current flow. The adjusted second gas phase G22B has a high flow velocity and a short residence time in region A, and the adjusted second gas phase G22B and the second substance M21 present in the second gas region have not reached an equilibrium state. In the present embodiment, the second gas reaction amount calculation unit 210-2 assumes that only a part of the adjusted second gas phase G22B in region A reaches an equilibrium state with the second substance M21 present in the second gas region, and according to the corresponding second gas reaction amount, divides the adjusted second gas phase G22B into a second gas reaction component GR2 and a second gas unreacted component Gr2 other than the second gas reaction component GR2.

[0091] The first substance reaction amount calculation unit 211-1 calculates the first substance reaction amount that contributes to the equilibrium state of the first substance M11 with the modified first gas phase G13. That is, the first substance reaction amount calculation unit 211-1 calculates, as the first substance reaction amount, the reaction amount that the first substance M11 is expected to undergo an equilibrium reaction with the modified first gas phase G13 during the residence time of the first substance M11 present in the gas region in region A. The first substance reaction amount can be calculated from the reaction rate equation. The first substance reaction amount calculation unit 211-1 obtains, from the first substance reaction amount, the first substance reaction component MR11 that contributes to the equilibrium state of the first substance M11 with the modified first gas phase G13.

[0092] Similar to the first gas reaction amount calculation unit 210-1, the first substance reaction amount calculation unit 211-1 can calculate the amount that reaches an equilibrium state between the modified first gas phase G13 and the first substance M11 present in the first gas region using the reaction rate equation.

[0093] Similar to the first gas reaction amount calculation unit 210-1, the first substance reaction amount calculation unit 211-1 can obtain the amount that reaches an equilibrium state between the modified first gas phase G13 and the first substance M11 present in the first gas region by adopting, for example, a model that can be calculated from the reaction rate according to the Arrhenius-type reaction rate equation and the residence time of the first substance M11 in region A. Then, the first substance reaction amount calculation unit 211-1 divides it into a first substance reaction component MR11 corresponding to the first substance reaction amount and a first substance unreacted component Mr11 corresponding to the remaining unreacted mass flow rate.

[0094] In the rotary kiln 1, as described above, the flow of the modified first gas phase G13 and the first substance M11 is in parallel flow. The modified first gas phase G13 has a high flow velocity and a short residence time in region A, and the modified first gas phase G13 and the first substance M11 have not reached an equilibrium state. In the present embodiment, the first substance reaction amount calculation unit 211-1 assumes that only a part of the first substance M11 in a certain region (region A) in the rotary kiln 1 reaches an equilibrium state with the modified first gas phase G13, and according to the corresponding first substance reaction amount, divides the first substance M11 into a first substance reaction component MR11 and a first substance unreacted component Mr11 other than the first substance reaction component MR11.

[0095] The second substance reaction amount calculation unit 211-2 calculates the reaction amount of the second substance M21 that contributes to the equilibrium state with the adjusted second gas phase G22B. That is, the second substance reaction amount calculation unit 211-2 calculates, as the second substance reaction amount, the reaction amount that the second substance M21 is expected to react in an equilibrium reaction with the adjusted second gas phase G22B during the residence time of the second substance M21 present in the second gas region in region A. The second substance reaction amount can be calculated from the reaction rate formula. The second substance reaction amount calculation unit 211-2 obtains, from the second substance reaction amount, a second substance reaction component MR21 of the second substance M21 that contributes to the equilibrium state with the adjusted second gas phase G22B.

[0096] Similar to the first gas reaction amount calculation unit 210-1, the second substance reaction amount calculation unit 211-2 can obtain the amount that reaches an equilibrium state between the adjusted second gas phase G22B and the second substance M21 present in the second gas region by adopting, for example, a model that can be calculated from the reaction rate according to the Arrhenius-type reaction rate formula and the residence time of the second substance M21 in region A. Then, the second substance reaction amount calculation unit 211-2 divides it into a second substance reaction component MR21 corresponding to the second substance reaction amount and a second substance unreacted component Mr21 corresponding to the remaining unreacted mass flow rate.

[0097] In the rotary kiln 1, as described above, the flow of the adjusted second gas phase G22B and the second substance M21 is in parallel flow. The adjusted second gas phase G22B has a high flow velocity and a short residence time in region A, and the adjusted second gas phase G22B and the second substance M21 have not reached an equilibrium state. In the present embodiment, the second substance reaction amount calculation unit 211-2 assumes that only a part of the adjusted second gas phase G22B in region A reaches an equilibrium state with the second substance M21, and according to the corresponding second substance reaction amount, divides the second substance M21 into a second substance reaction component MR21 and a second substance unreacted component Mr21 other than the second substance reaction component MR21.

[0098] The first bed layer reaction amount calculation unit 212-1 calculates the first bed layer reaction amount contributing to the equilibrium state with the first inflow gas G221-1 of the modified first bed layer S13B using a reaction rate equation. That is, the first bed layer reaction amount calculation unit 212-1 calculates the first bed layer reaction amount that the modified first bed layer S13B is expected to react in equilibrium with the first inflow gas G221-1 during the movement time of the modified first bed layer S13B in region A using the first reaction rate equation.

[0099] The first bed layer reaction amount calculation unit 212-1 can calculate the first bed layer reaction amount by using, for example, the reaction rate parameters set for each individual substance included in the modified first bed layer S13B in region A, the residence time in region A, etc. in the first reaction rate equation. Then, the first bed layer reaction amount calculation unit 212-1 calculates the first bed layer reaction component SR1 that contributes to the equilibrium reaction with the first inflow gas G221-1 from the first bed layer reaction amount in the modified first bed layer S13B.

[0100] That is, the first bed layer reaction component SR1 is the reaction component that contributes to the equilibrium reaction with the first inflow gas G221-1 in the modified first bed layer S13B.

[0101] The first bed layer reaction amount calculation unit 212-1 can calculate the first bed layer reaction amount by adopting, for example, a reaction rate conforming to an Arrhenius-type reaction rate equation, a model that can be calculated from the moving time of the modified first bed layer S13B in region A, etc. The first bed layer reaction amount calculation unit 212-1 calculates, for the component estimated to have reached an equilibrium state between the modified first bed layer S13B and the first inflow gas G221-1 only by the residence time of the modified first bed layer S13B in the Arrhenius-type reaction rate equation, the first bed layer reaction amount of the modified first bed layer S13B contributing to the equilibrium reaction with the first inflow gas G221-1. Then, the first bed layer reaction amount calculation unit 212-1 divides it into the first bed layer reaction component SR1 corresponding to the first bed layer reaction amount and the first bed layer unreacted component Sr1 corresponding to the remaining unreacted mass flow rate.

[0102] In the rotary kiln 1, the flows of the modified first bed layer S13B and the first inflow gas G221-1 are countercurrent. The flow rate of the first inflow gas G221-1 is large and the residence time in region A is short, and the modified first bed layer S13B and the first inflow gas G221-1 have not reached an equilibrium state. In this embodiment, the first bed layer reaction amount calculation unit 212-1 assumes that only a part of the modified first bed layer S13B in a certain region (region A) in the rotary kiln 1 reaches an equilibrium state with the first inflow gas G221-1, and divides the modified first bed layer S13B into the first bed layer reaction component SR1 and the first bed layer unreacted component Sr1 other than the first bed layer reaction component according to the corresponding first bed layer reaction amount.

[0103] The second bed layer reaction amount calculation unit 212-2 calculates the bed reaction amount contributing to the equilibrium state between the modified second bed layer S21B and the second inflow gas G221-2 by using the reaction rate equation. That is, the second bed layer reaction amount calculation unit 212-2 calculates the second bed layer reaction amount expected to reach an equilibrium reaction between the modified second bed layer S21B and the second inflow gas G221-2 during the moving time of the modified second bed layer S21B in region A.

[0104] The second bed layer reaction amount calculation unit 212-2 calculates a second bed layer reaction component SR2 that contributes to the equilibrium state with the second inflow gas G221-2 in the corrected second bed layer S21B from the obtained second bed layer reaction amount.

[0105] Similar to the first bed layer reaction amount calculation unit 212-1, for the calculation of the reaction amount reaching the equilibrium state between the corrected second bed layer S21B and the second inflow gas G221-2, the second bed layer reaction amount calculation unit 212-2 can adopt, for example, a reaction rate according to an Arrhenius-type reaction rate equation and a model that can be calculated from the residence time of the preliminary corrected second bed layer S21A in region A. The second bed layer reaction amount calculation unit 212-2 calculates the second bed layer reaction amount of the corrected second bed layer S21B that contributes to the reaction with the second inflow gas G2211-2 only for the component estimated to have reached the equilibrium state between the corrected second bed layer S21B and the second inflow gas G221-2 based on the residence time of the corrected second bed layer S21B using an Arrhenius-type reaction rate equation. Then, the second bed layer reaction amount calculation unit 212-2 divides it into a second bed layer reaction component SR2 corresponding to the second bed layer reaction amount and a second bed layer unreacted component Sr2 corresponding to the remaining unreacted mass flow rate.

[0106] The second bed layer reaction component SR2 is the reaction component that contributes to the equilibrium reaction with the second inflow gas G221-2 in the corrected second bed layer S21B.

[0107] In the rotary kiln 1, the flows of the corrected second bed layer S21B and the second inflow gas G221-2 are countercurrent. The flow rate of the second inflow gas G221-2 is large and the residence time in region A is short, and the corrected second bed layer S21B and the second inflow gas G221-2 have not reached the equilibrium state. In this embodiment, the second bed layer reaction amount calculation unit 212-2 assumes that only a part of the corrected second bed layer S21B in region A reaches the equilibrium state with the second inflow gas G221-2, and divides the corrected second bed layer S21B into a second bed layer reaction component SR2 and a second bed layer unreacted component Sr2 other than that according to the corresponding second bed layer reaction amount.

[0108] The first incoming gas reaction amount calculation unit 213-1 calculates the first incoming gas reaction amount that contributes to the equilibrium state with the modified first bed layer S13B of the first incoming gas G221-1. That is, the first incoming gas reaction amount calculation unit 213-1 calculates the first incoming gas reaction amount that is expected to react in equilibrium with the modified first bed layer S13B during the residence time of the first incoming gas G221-1 present in the bed layer in region A. The first incoming reaction amount at the residence time of the first incoming gas G221-1 in region A can be calculated using the reaction rate parameters set for each individual substance and the residence time within region A, and can be calculated using, for example, a reaction rate formula. Then, the first incoming gas reaction amount calculation unit 213-1 calculates the first incoming gas reaction component GR221 that contributes to the equilibrium reaction with the modified first bed layer S13B from the first incoming gas reaction amount among the first incoming gas G221-1.

[0109] That is, the first incoming gas reaction component GR221 is the reaction component of the first incoming gas G221-1 that contributes to the equilibrium state of the first modified bed layer S13B among the first incoming gas G221-1.

[0110] The first incoming gas reaction amount calculation unit 213-1 can calculate the reaction amount for the first incoming gas G221-1 to reach an equilibrium state with the modified first bed layer S13B using, for example, a model similar to the model used by the first bed layer reaction amount calculation unit 212-1.

[0111] Then, the first incoming gas reaction amount calculation unit 213-1 divides it into the first incoming gas reaction component GR221 corresponding to the first incoming gas reaction amount and the first incoming gas unreacted component Gr221 corresponding to the mass flow rate of the unreacted part other than the remaining first incoming gas reaction component GR221.

[0112] In the rotary kiln 1, as described above, the flow of the first inflowing gas G221-1 and the modified first bed layer S13B is countercurrent. The first inflowing gas G221-1 has a high flow velocity and a short residence time in region A, and the first inflowing gas G221-1 and the modified first bed layer S13B have not reached an equilibrium state. In the present embodiment, the first inflowing gas reaction amount calculation unit 213-1 assumes that only a part of the first inflowing gas G221-1 in region A of the rotary kiln 1 reaches an equilibrium state with the modified first bed layer S13B, and according to the corresponding gas reaction amount, divides the first inflowing gas G221-1 into a first inflowing gas reaction component GR221 and a non-reacted component Gr221 of the other first inflowing gas.

[0113] The second inflowing gas reaction amount calculation unit 213-2 calculates the second inflowing gas reaction amount that contributes to the equilibrium state of the second inflowing gas G221-2 with the modified second bed layer S21B. That is, the second inflowing gas reaction amount calculation unit 213-2 calculates the second inflowing gas reaction amount that is expected to reach an equilibrium reaction with the modified second bed layer S21B during the residence time of the second inflowing gas G221-2 present in the bed layer in region A. The second inflowing reaction amount at the residence time of the second inflowing gas G221-2 in region A can be calculated using the reaction rate parameters set for each individual substance and the residence time in region A, and can be calculated, for example, from the reaction rate formula. Then, the second inflowing gas reaction amount calculation unit 213-2 calculates a second inflowing gas reaction component GR222 that contributes to the equilibrium reaction of the second inflowing gas G221-2 with the modified second bed layer S21B.

[0114] That is, the second inflowing gas reaction component GR222 is the reaction component of the second inflowing gas G221-2 that contributes to the equilibrium state of the second modified bed layer S21B among the second inflowing gas G221-2.

[0115] The second inflowing gas reaction amount calculation unit 213-2 can calculate the reaction amount at which the second inflowing gas G221-2 reaches an equilibrium state with the modified second bed layer S21B using, for example, a model similar to the model used in the second bed layer reaction amount calculation unit 212-2.

[0116] The second inflow gas reaction amount calculation unit 213-2 divides the second inflow gas reaction component GR222 corresponding to the second inflow gas reaction amount and the second inflow gas unreacted component Gr222 corresponding to the mass flow rate of the unreacted part other than the remaining second inflow gas reaction component.

[0117] In the rotary kiln 1, as described above, the flow of the second inflow gas G221-2 and the modified second bed layer S21B is countercurrent. The flow velocity of the second inflow gas G221-2 is high and the residence time in region A is short, and the second inflow gas G221-2 and the modified second bed layer S21B have not reached an equilibrium state. In the present embodiment, the second inflow gas reaction amount calculation unit 213-2 assumes that only a part of the modified second bed layer S21B in region A in the rotary kiln 1 reaches an equilibrium state with the second inflow gas G221-2, and according to the corresponding second inflow gas reaction amount, divides the second inflow gas G221-2 into the second inflow gas reaction component GR222 and the other second inflow gas unreacted component Gr222.

[0118] The first gas region equilibrium reaction calculation unit 214-1 calculates the equilibrium reaction between the modified first gas phase G13 flowing through the first gas region in region A and the first substance M11 existing in the first gas region. That is, the first gas region equilibrium reaction calculation unit 214-1 performs the equilibrium reaction calculation of the first gas region using the gas reaction component GR1 contributing to the equilibrium reaction with the first substance M11 existing in the first gas region among the modified first gas phase G13 and the first substance reaction component MR11 contributing to the equilibrium reaction with the first gas reaction component GR1 among the first substances M11, and calculates at least the composition, amount, etc. of the first gas reactant GR1 and the first substance reactant MR11, and the change and flow rate of each heat quantity when the first gas reactant GR1 and the first substance reactant MR11 reach an equilibrium state.

[0119] The first gas region equilibrium reaction calculation unit 214-1 has the function of performing equilibrium reaction calculations to determine the type, phase, flow rate, etc. of substances (product substances) that may be generated by the reaction between the first gas reaction component GR1 and the first substance reaction component MR11, by presetting the type, phase, flow rate, etc. of the substances (product substances) in advance, so that the free energy of the product substances is minimized. That is, the first gas region equilibrium reaction calculation unit 214-1 can perform equilibrium reaction calculations based on the Gibbs energy minimization method. The first gas region equilibrium reaction calculation unit 214-1 calculates and outputs the changes in heat quantity, flow rate, type, composition, amount, and phase of each when the first gas reaction component GR1 and the first substance reaction component MR11 reach an equilibrium state.

[0120] The first gas region equilibrium reaction calculation unit 214-1 calculates the generated gas resulting from the reaction between the first gas reaction component GR1 and the first substance reaction component MR11, and the unused portion of the first gas reaction component GR1, as the first gas generation component GP1.

[0121] Also, the first gas region equilibrium reaction calculation unit 214-1 calculates the product substance resulting from the reaction between the first gas reaction component GR1 and the first substance reaction component MR11, and the unused portion of the first substance reaction component MR11, as the first substance generation component MP11.

[0122] That is, the first gas generation component GP1 is the gas phase obtained by summing the gas phase generated by the reaction between the first gas reaction component GR1 and the first substance reaction component MR11 and the unused portion of the first gas reaction component GR1. The first substance generation component MP11 is the sum of the product substance generated by the reaction between the first gas reaction component GR1 and the first substance reaction component MR11 and the unused portion of the first substance reaction component MR11.

[0123] The second gas region equilibrium reaction calculation unit 214-2 calculates the equilibrium reaction between the second gas phase G22B flowing through the second gas region in region A and the second substance M21 present in the second gas region. That is, the second gas region equilibrium reaction calculation unit 214-2 uses the gas reaction component GR2 contributing to the equilibrium reaction with the second substance M21 present in the second gas region of the second gas phase G22B and the second substance reaction component MR21 contributing to the equilibrium reaction with the second gas reaction component GR2 among the second substances M21 to perform the equilibrium reaction calculation of the second gas region, and calculates at least the composition, amount, etc. of the second gas reactant GR2 and the first substance reactant MR21, and the change in heat quantity and flow rate of each when the second gas reactant GR2 and the second substance reactant MR21 reach an equilibrium state.

[0124] The second gas region equilibrium reaction calculation unit 214-2 has a function of presetting the type, phase, flow rate, etc. of substances (product substances) that may be generated by the reaction between the second gas reaction component GR2 and the second substance reaction component MR21, and performing equilibrium reaction calculation so that the free energy of the product substances is minimized to determine the type, phase, flow rate, etc. of the product substances. That is, the second gas region equilibrium reaction calculation unit 214-2 can perform equilibrium reaction calculation based on the Gibbs energy minimization method. The second gas region equilibrium reaction calculation unit 214-2 calculates and outputs the change in heat quantity, flow rate, type, composition, amount, and phase of each when the second gas reaction component GR2 and the second substance reaction component MR21 reach an equilibrium state.

[0125] The second gas region equilibrium reaction calculation unit 214-2 calculates the generated gas generated by the reaction between the second gas reaction component GR2 and the second substance reaction component MR21 and the unused portion of the second gas reaction component GR2 as the second gas generation component GP2.

[0126] In addition, the second gas region equilibrium reaction calculation unit 214-2 calculates the product substance generated by the reaction between the second gas reaction component GR2 and the second substance reaction component MR21 and the unused portion of the second substance reaction component MR21 as the second substance generation component MP21.

[0127] That is, the second gas generation component GP2 is a gas phase obtained by summing a gas phase generated by the reaction of the second gas reaction component GR2 and the second substance reaction component MR21, and an unused portion of the second gas reaction component GR2. The first substance generation component MP21 is obtained by summing a product substance generated by the reaction of the gas reaction component GR2 and the second substance reaction component MR21, and an unused portion of the second substance reaction component MR21.

[0128] The bed layer equilibrium reaction calculation unit 215 calculates the equilibrium reaction between the bed layer in region A and the inflowing gas G221.

[0129] The bed layer equilibrium reaction calculation unit 215 includes a first bed layer equilibrium reaction calculation unit 215-1 and a second bed layer equilibrium reaction calculation unit 215-2.

[0130] The first bed layer equilibrium reaction calculation unit 215-1 calculates the equilibrium reaction between the modified first bed layer S13B moving in the bed layer in region A and the first inflowing gas G221 existing in the modified first bed layer S13B. That is, the first bed layer equilibrium reaction calculation unit 215-1 uses the first bed layer reaction component SR1 that contributes to the equilibrium reaction with the first inflowing gas G221 in the modified first bed layer S13B, and the first inflowing gas reaction component GR221 that contributes to the equilibrium state with the modified first bed layer S13B in the first inflowing gas G221 to perform the equilibrium reaction calculation of the first bed layer. Thereby, the first bed layer equilibrium reaction calculation unit 215 calculates at least the changes in the respective heat amounts and the flow rates when the first bed layer reactant SR1 and the first inflowing gas reactant GR221 reach an equilibrium state.

[0131] The first bed layer equilibrium reaction calculation unit 215-1 has the function of performing equilibrium reaction calculations to determine the types, phases, flow rates, etc. of substances (product substances) that may be generated by the reaction between the first bed layer reaction component SR1 and the first inflow gas reaction component GR221, by presetting the types, phases, flow rates, etc. of the product substances in advance, so that the free energy of the product substances is minimized. That is, similar to the gas region equilibrium reaction calculation unit 213, the first bed layer equilibrium reaction calculation unit 215-1 can perform equilibrium reaction calculations based on the Gibbs energy minimization method. The first bed layer equilibrium reaction calculation unit 215-1 calculates and outputs the changes in the amount of heat, flow rates, types, compositions, amounts, and phases of each when the first bed layer reaction component SR1 and the first inflow gas reaction component GR221 reach an equilibrium state.

[0132] The first bed layer equilibrium reaction calculation unit 215-1 calculates the first bed layer generated component SP1 by combining the product substances generated by the reaction between the first bed layer reaction component SR1 and the first inflow gas reaction component GR221, and the unused portion of the first bed layer reaction component SR1. Also, the first bed layer equilibrium reaction calculation unit 215-1 calculates the first inflow gas generated component GP221 by combining the product gas generated by the reaction between the first bed layer reaction component SR1 and the first inflow gas reaction component GR221, and the unused portion of the first inflow gas reaction component GR221.

[0133] That is, the first bed layer generated component SP1 is the sum of the product substances generated by the reaction between the first bed layer reaction component SR1 and the first inflow gas reaction component GR221, and the unused portion of the first bed layer reaction component SR1. The first inflow gas generated component GP221 is the sum of the product gas generated by the reaction between the first bed layer reaction component SR1 and the first inflow gas reaction component GR221, and the unused portion of the first inflow gas reaction component GR221.

[0134] The second bed layer equilibrium reaction calculation unit 215-2 calculates the equilibrium reaction between the modified second bed layer S21B moving within the region A and the second inflow gas G221-2 existing in the bed layer. That is, the second bed layer equilibrium reaction calculation unit 215-2 uses the second bed layer reaction component SR2 that contributes to the equilibrium reaction with the second inflow gas G221-2 in the modified second bed layer S21B and the second inflow gas reaction component GR222 that contributes to the equilibrium state with the modified second bed layer S21B in the second inflow gas G221-2 to calculate the equilibrium reaction of the modified second bed layer S21B. Thereby, at least the composition, amount, etc. of the second bed layer reactant SR2 and the second inflow gas reaction component GR222, and the changes in the respective heat quantities and flow rates when the second bed layer reactant SR2 and the second inflow gas reaction component GR222 reach an equilibrium state are calculated.

[0135] The second bed layer equilibrium reaction calculation unit 215-2 has a function of performing equilibrium reaction calculation so that the free energy of the product substance is minimized by presetting the type, phase, flow rate, etc. of the substance (product substance) that may be generated by the reaction between the second bed layer reaction component SR2 and the second inflow gas reaction component GR222, and determining the type, phase, flow rate, etc. of the product substance. That is, similar to the first bed layer equilibrium reaction calculation unit 215-1, the second bed layer equilibrium reaction calculation unit 215-2 can perform equilibrium reaction calculation based on the Gibbs energy minimization method. The second bed layer equilibrium reaction calculation unit 215-2 calculates and outputs the changes in the respective heat quantities, flow rates, types, compositions, amounts, and phases of the second bed layer reaction component SR2 and the second inflow gas reaction component GR222 when they reach an equilibrium state.

[0136] The second bed layer equilibrium reaction calculation unit 215-2 calculates the second bed layer generated component SP2 as the product substance generated by the reaction between the second bed layer reaction component SR2 and the second inflow gas reaction component GR222 and the unused portion of the second bed layer reaction component SR2. Also, the second bed layer equilibrium reaction calculation unit 215-2 calculates the second inflow gas generated component GP222 as the product gas generated by the reaction between the second bed layer reaction component SR2 and the second inflow gas reaction component GR222 and the unused portion of the second inflow gas reaction component GR222.

[0137] That is, the second bed layer generation component SP2 is the sum of the product substance generated by the reaction of the second bed layer reaction component SR2 and the second inflow gas reaction component GR222, and the unused portion of the second bed layer reaction component SR2. Also, the second inflow gas generation component GP222 is the sum of the generated gas generated by the reaction of the second bed layer reaction component SR2 and the second inflow gas reaction component GR222, and the unused portion of the first inflow gas reaction component GR222.

[0138] The first mixed gas phase calculation unit 216-1 has a function of mixing a plurality of flows, and calculates the flow rate, composition data, etc. of the first mixed gas phase G14, which is a mixed gas phase obtained by mixing the unreacted portion Gr1 of the first gas divided by the first gas reaction amount calculation unit 210-1 and the first gas generation component GP1 generated by the first gas region equilibrium reaction calculation unit 214-1.

[0139] The second mixed gas phase calculation unit 216-2 has a function of mixing a plurality of flows, and calculates the flow rate, composition data, etc. of the second mixed gas phase G23, which is a mixed gas phase obtained by mixing the unreacted portion Gr2 of the second gas divided by the second gas reaction amount calculation unit 210-2 and the second gas generation component GP2 generated by the second gas region equilibrium reaction calculation unit 214-2.

[0140] The first mixture substance calculation unit 217-1 has a function of mixing a plurality of flows, and calculates the flow rate, composition data, etc. of the first mixture substance M12, which is obtained by mixing the unreacted portion Mr11 of the first substance divided by the first substance reaction amount calculation unit 211-1 and the first substance generation component MP11 generated by the first gas region equilibrium reaction calculation unit 214-1.

[0141] The second mixture substance calculation unit 217-2 has a function of mixing a plurality of flows, and calculates the flow rate, composition data, etc. of the second mixture substance M22, which is obtained by mixing the unreacted portion Mr21 of the second substance divided by the second substance reaction amount calculation unit 211-2 and the second substance generation component MP21 generated by the second gas region equilibrium reaction calculation unit 214-2.

[0142] The first mixed bed layer calculation unit 218-1 has a function of mixing a plurality of flows, and calculates the flow rate, composition data, etc. of the first mixed bed layer S14 obtained by mixing the unreacted portion Sr1 of the first bed layer divided by the first bed layer reaction amount calculation unit 212-1 and the first bed layer generated component SP1 generated by the first bed layer equilibrium reaction calculation unit 215-1.

[0143] The second mixed bed layer calculation unit 218-2 has a function of mixing a plurality of flows, and calculates the flow rate, composition data, etc. of the second mixed bed layer S22 obtained by mixing the unreacted portion Sr2 of the second bed layer divided by the second bed layer reaction amount calculation unit 212-2 and the second bed layer generated component SP2 generated by the second bed layer equilibrium reaction calculation unit 215-2.

[0144] The first mixed inflow gas calculation unit 219-1 has a function of mixing a plurality of flows, and calculates the flow rate, composition data, etc. of the first mixed inflow gas G222A obtained by mixing the unreacted portion Gr221 of the first inflow gas divided by the first inflow gas reaction amount calculation unit 213-1 and the first inflow gas generated component GP221 generated by the first bed layer equilibrium reaction calculation unit 215-1.

[0145] The second mixed inflow gas calculation unit 219-2 has a function of mixing a plurality of flows, and calculates the flow rate, composition data, etc. of the second mixed inflow gas G222B obtained by mixing the unreacted portion Gr222 of the second inflow gas divided by the second inflow gas reaction amount calculation unit 213-2 and the second inflow gas generated component GP222 generated by the second bed layer equilibrium reaction calculation unit 215-2.

[0146] The first mixture amount correction unit 220 corrects the amount of the first mixture M12 in consideration of the amount of movement of the third substance M3 in the first gas region of the second mixture M23 to the first mixture M12 and the amount of movement of the fourth substance M4 contained in the first mixture M12 in the first gas region to the second mixture M23, and calculates the first corrected mixture M13.

[0147] The first mixture amount correction unit 220 corrects the amount of the first mixture M12 by adding the amount of movement of the first mixture M12 to the first gas region of the third substance M3 contained in the second mixture M23 and dividing by the amount of movement of the fourth substance M4 contained in the first mixture M12 present in the first gas region to the second mixture M23.

[0148] The inflow gas summing unit 221 has a function of mixing a plurality of flows, sums the first mixed inflow gas G222A generated by the first mixed inflow gas calculation unit 219-1 and the second mixed inflow gas G222B generated by the second mixed inflow gas calculation unit 219-2, and calculates the mixed inflow gas G223.

[0149] The gas amount summing unit 222 has a function of mixing a plurality of flows, sums the second mixed gas phase G23 generated by the second mixed gas phase calculation unit 216-2 and the mixed inflow gas G223 generated by the inflow gas summing unit 221, and calculates the second summed gas phase G24.

[0150] The first-order amount correction unit 223 of the second mixture considers the amount of movement of the seventh substance M7 containing at least one of a solid substance or a liquid substance contained in the first bed layer S12 with additives generated by the first bed layer preliminary amount correction unit 205 to the second mixture M22, the amount of movement of the eighth substance M8 containing at least one of a solid substance or a liquid substance contained in the second mixture M22 generated by the second mixture calculation unit 217-2 to the first bed layer S12 with additives, and the amount of movement of the ninth substance M9 contained in the second mixed bed layer S21 generated by the second bed layer preliminary amount correction unit 207 to the second mixture M22, corrects the amount of the second mixture M22, and calculates the second corrected mixture M23.

[0151] The first-order amount correction unit 223 of the second mixture corrects the amount of the second mixture M22 by adding the amount of movement of the seventh substance M7 contained in the first bed layer S12 with additives to the second mixture M22 and the amount of movement of the ninth substance M9 contained in the second mixed bed layer S21 to the second mixture M22, and dividing by the amount of movement of the eighth substance M8 contained in the second mixture M22 to the first bed layer S12 with additives.

[0152] The second mixture material secondary material quantity correction unit 224 corrects the material quantity of the second corrected mixture material M23 in consideration of the transfer amount of the third substance M3 existing in the first gas region of the second corrected mixture material M23 generated by the second mixture material primary material quantity correction unit 223 to the first mixture material M12 and the transfer amount of the fourth substance M4 contained in the first mixture material M12 to the second corrected mixture material M23, and calculates the second corrected mixture material M24.

[0153] The second mixture material secondary material quantity correction unit 224 excludes the transfer amount of the third substance M3 contained in the second corrected mixture material M23 to the first mixture material M12 and adds the transfer amount of the fourth substance M4 contained in the first mixture material M12 to the second corrected mixture material M23, thereby correcting the material quantity of the second corrected mixture material M23.

[0154] In this embodiment, at least one of the first gas region equilibrium reaction calculation unit 214-1 and the second gas region equilibrium reaction calculation unit 214-2 may include two or more. By separating at least one of the first gas phase G11 and the second gas phase G21 into a plurality of gas phases, the gas regions in the rotary kiln 1 can be divided into three or more types of multiple systems.

[0155] In this embodiment, at least one of the first bed layer equilibrium reaction calculation unit 215-1 and the second bed layer equilibrium reaction calculation unit 215-2 may include two or more. By separating at least one of the first bed layer S11 and the first bed layer S21 into a plurality of bed layers, the bed layers in the rotary kiln 1 can be divided into three or more types of multiple systems.

[0156] In this embodiment, the in-furnace reaction calculation device 20 may not include the distribution unit 201 when the combustion material 240 is not dropped into the rotary kiln 1.

[0157] In this embodiment, the in-furnace reaction calculation device 20 may not include one or more of a first bed layer preliminary material quantity correction unit 205, a first bed layer material quantity correction unit 206, a second bed layer preliminary material quantity correction unit 207, and a second bed layer material quantity correction unit 208. For example, when the in-furnace reaction calculation device 20 does not include the first bed layer preliminary material quantity correction unit 205 and the first bed layer material quantity correction unit 206, the first bed layer reaction quantity calculation unit 212-1 uses the first bed layer S12 instead of the corrected first bed layer S13B. Further, when the in-furnace reaction calculation device 20 does not include the second bed layer preliminary material quantity correction unit 207 and the second bed layer material quantity correction unit 208, the second bed layer reaction quantity calculation unit 212-2 uses the second bed layer S21 instead of the preliminary corrected second bed layer S21A.

[0158] In this embodiment, the in-furnace reaction calculation device 20 may calculate the heat transfer in the furnace using models such as radiation, conduction, and convection as necessary.

[0159] <Calculation method of in-furnace reaction> Next, the in-furnace reaction calculation method according to this embodiment will be described using the in-furnace reaction calculation device according to this embodiment. The in-furnace reaction calculation method according to this embodiment is to move the raw ore supplied from the charging end 14A side of the rotary kiln 1 having the configuration shown in FIG. 1 toward the discharging end 14B side, and while moving, introduce the combustion material 240 from the middle of the movement, and bring it into contact with the combustion gas supplied from the burner 16 provided on the discharging end 14B side of the raw ore to dry and perform reduction.

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

[0161] When the combustion material 240 is dropped into the rotary kiln 1 (step S11: Yes), the in-furnace reaction calculation device 20 provides, as input substances, substances such as volatile components assuming the gas phase and solid substances assuming the bed layer, and uses the distribution unit 201 to distribute the combustion material 240 added into the rotary kiln 1 to the gas phase and the bed layer in terms of mass flow rate (distribution step: step S12).

[0162] Next, the in-furnace reaction calculation device 20 provides, as input substances, the first gas phase G11 which is the combustion gas flowing from one adjacent region (region (A + 1)) into region A and the added gas phase AG11 which is the gas phase of the combustion material 240 distributed by the distribution unit 201. The in-furnace reaction calculation device 20 mixes the first gas phase G11 and the added gas phase AG11 to calculate the added-component-containing first gas phase G12 including the first gas phase G11 and the added gas phase AG11 (gas-phase mixing step of the first gas phase and the added gas phase: step S13).

[0163] Next, the in-furnace reaction calculation device 20 uses the material quantity correction unit 202 of the first gas phase to correct the material quantity of the added-component-containing first gas phase G12 flowing through the first gas region and calculates the corrected first gas phase G13 (correction step of the material quantity of the second gas phase: step S14).

[0164] Next, the in-furnace reaction calculation device 20 uses the material quantity correction unit 203 of the second gas phase to correct the material quantity of the second gas phase G21 flowing through the second gas region and calculates and obtains the corrected second gas phase G22A (correction step of the material quantity of the second gas phase: step S15).

[0165] Next, the in-furnace reaction calculation device 20 uses the gas mixing quantity calculation unit 204 to calculate the flow rate of the inflow gas G221 moving from the corrected second gas phase G22A in the second gas region to the bed layer (gas mixing quantity calculation step: step S16).

[0166] Next, the in-furnace reaction calculation device 20 supplies, as input substances, a first bed layer S11 in which the bed layer S1 flowing into region A from the other adjacent region (region (A-1)) is previously distributed, and an added bed layer AS11 that is a bed layer of the combustion material distributed by the distribution unit 201. The in-furnace reaction calculation device 20 mixes the first bed layer S11 and the added bed layer AS11 to calculate an added-content-containing first bed layer S12 (bed layer mixing step of the first bed layer and the added bed layer: step S17).

[0167] Next, the in-furnace reaction calculation device 20 supplies, as input substances, the added-content-containing first bed layer S12, the amount of movement of the seventh substance M7 contained in the added-content-containing first bed layer S12 to the second mixed substance M22 present in the second gas region, and the amount of movement of the eighth substance M8 contained in the second mixed substance M22 to the added-content-containing first bed layer S12.

[0168] The in-furnace reaction calculation device 20 uses the first bed layer preliminary material amount correction unit 205 to divide by the amount of movement of the seventh substance M7 contained in the added-content-containing first bed layer S12 to the second mixed substance M22, and add the amount of movement of the eighth substance M8 contained in the second mixed substance M22 to the added-content-containing first bed layer S12, thereby correcting the material amount of the added-content-containing first bed layer S12 and calculating a preliminarily corrected first bed layer S13A (first bed layer preliminary material amount correction step: step S18).

[0169] Next, the in-furnace reaction calculation device 20 supplies, as input substances, the preliminarily corrected first bed layer S13A, the amount of movement of the sixth substance M6 contained in the preliminarily corrected first bed layer S13A to the second bed layer S21A, and the fifth substance M5 in the second bed layer S21A that is removed by the second bed layer material amount correction unit 208.

[0170] The in-furnace reaction calculation device 20 uses the first bed layer material quantity correction unit 206 to remove the transfer amount of the sixth substance M6 contained in the preliminary correction first bed layer S13A from the preliminary correction first bed layer S13A, and adds the transfer amount of the fifth substance M5 in the second bed layer S21A removed by the second bed layer material quantity correction unit 208 to the preliminary correction first bed layer S13A. Thereby, the in-furnace reaction calculation device 20 corrects the material quantity of the preliminary correction first bed layer S13A and calculates the corrected first bed layer S13B (preliminary correction first bed layer material quantity correction step: step S19).

[0171] Next, the in-furnace reaction calculation device 20 gives the second bed layer S21 and the transfer amount of the ninth substance M9 contained in the second bed layer S21 to the second mixed substance M22 existing in the second gas region as input substances.

[0172] The in-furnace reaction calculation device 20 uses the second bed layer preliminary material quantity correction unit 207 to correct the material quantity of the second bed layer S21 by removing the transfer amount of the ninth substance M9 contained in the second bed layer S21 to the second mixed substance M22, and calculates the preliminary correction second bed layer S21A.

[0173] Next, the in-furnace reaction calculation device 20 gives the preliminary correction second bed layer S21A, the transfer amount of the sixth substance M6 in the preliminary correction first bed layer S13A removed by the first bed layer material quantity correction unit 206, and the transfer amount of the fifth substance M5 in the preliminary correction second bed layer S21A removed by the second bed layer material quantity correction unit 208 as input substances.

[0174] The in-furnace reaction calculation device 20 uses the second bed layer material quantity correction unit 208 to remove the transfer amount of the fifth substance M5 contained in the preliminary correction second bed layer S21A from the preliminary correction second bed layer S21A, and adds the transfer amount of the sixth substance M6 removed by the first bed layer material quantity correction unit 206 to the preliminary correction second bed layer S21A. Thereby, the in-furnace reaction calculation device 20 corrects the material quantity of the preliminary correction second bed layer S21A and calculates the corrected second bed layer S21B (second bed layer material quantity correction step: step S21).

[0175] Next, the in-furnace reaction calculation device 20 performs a first gas region equilibrium reaction (step of performing the first gas region equilibrium reaction: step S22).

[0176] Next, the in-furnace reaction calculation device 20 performs a second gas region equilibrium reaction (step of performing the second gas region equilibrium reaction: step S23).

[0177] Next, the in-furnace reaction calculation device 20 performs a first bed layer equilibrium reaction (step of performing the first bed layer equilibrium reaction: step S24).

[0178] Next, the in-furnace reaction calculation device 20 performs a second bed layer equilibrium reaction (step of performing the second bed layer equilibrium reaction: step S25).

[0179] Next, the in-furnace reaction calculation device 20 uses the inflow gas summation unit 221 to sum the second mixed inflow gas G222B generated in the step of performing the second bed layer equilibrium reaction (step S25) to the first mixed inflow gas G222A generated in the step of performing the first bed layer equilibrium reaction (step S24) to obtain a mixed inflow gas G223 (inflow gas summation step: step S26).

[0180] Next, the in-furnace reaction calculation device 20 uses the gas amount summation unit 222 to sum the mixed inflow gas G223 generated in the inflow gas summation step (step S26) to the second mixed gas phase G23 generated in the step of performing the second gas region equilibrium reaction (step S19) to obtain a second summed gas phase G24 (gas amount summation step: step S27).

[0181] The in-furnace reaction calculation device 20 moves the second summed gas phase G24 to a region (region (A - 1)) on the charging end 14A side rather than region A.

[0182] Next, the in-furnace reaction calculation device 20 uses the first mixture material quantity correction unit 220 to remove the amount of the fourth substance M4 in the first mixture material M12 that moves to the second gas region during the implementation process of the first gas region equilibrium reaction (step S19), and adds the amount of the third substance M3 in the second mixture material M23 that moves to the first mixture material M12, thereby correcting the material quantity of the first mixture material M12 and calculating the first corrected mixture material M13 (first mixture material quantity correction step: step S28).

[0183] The in-furnace reaction calculation device 20 moves the first corrected mixture material M13 to the region (region (A - 1)) on the charging end 14A side rather than region A.

[0184] Next, the in-furnace reaction calculation device 20 uses the first secondary material quantity correction unit 223 of the second mixture material to add the amount of the seventh substance M7 in the additive-containing first bed layer S12 that moves to the second mixture material M22 during the preliminary material quantity correction step of the first bed layer (step S16) and the amount of the ninth substance M9 in the second mixed bed layer S21 that moves to the second mixture material M22 during the second bed layer preliminary material quantity correction unit 207, and at the same time removes the amount of the eighth substance M8 in the second mixture material M22 that moves to the additive-containing first bed layer S12 during the implementation process of the second bed layer equilibrium reaction (step S25). Thereby, the in-furnace reaction calculation device 20 corrects the material quantity of the second mixture material M22 and calculates the second corrected mixture material M23 (first secondary material quantity preliminary correction step of the second mixture material: step S29).

[0185] Next, the in-furnace reaction calculation device 20 uses the second secondary material quantity correction unit 224 of the second mixture material to remove the amount of the third substance M3 in the second corrected mixture material M23 that moves to the first mixture material M12 during the second mixture material quantity preliminary correction step (step S29), and adds the amount of the fourth substance M4 in the first mixture material M12 that moves to the second corrected mixture material M23 during the first mixture material quantity correction step (step S28), thereby correcting the material quantity of the second corrected mixture material M23 and calculating the second corrected mixture material M24 (second secondary material quantity correction step of the second mixture material: step S30).

[0186] The in-furnace reaction calculation device 20 moves the second modified mixture M24 to a region (region (A-1)) on the charging end 14A side rather than region A.

[0187] Next, the process of performing the first gas region equilibrium reaction (step S22) will be described. FIG. 6 is a flowchart showing the operation of the process of performing the first gas region equilibrium reaction (step S22) in FIG. 5. As shown in FIG. 6, in the process of performing the first gas region equilibrium reaction (step S22), the in-furnace reaction calculation device 20 uses the first gas reaction amount calculation unit 210-1 to calculate the first gas reaction amount that contributes to the equilibrium state with the first substance M11 in the modified first gas phase G13 (first gas reaction amount calculation step: step S221).

[0188] That is, when it is assumed that the modified first gas phase G13 reacts with the first substance M11 existing in the first gas region to reach an equilibrium state, the in-furnace reaction calculation device 20 calculates the gas reaction amount that contributes to the reaction when the modified first gas phase G13 reacts with the first substance M11 existing in the first gas region to reach an equilibrium state.

[0189] Then, based on the calculation result of the first gas reaction amount, the in-furnace reaction calculation device 20 divides it into a first gas reaction component GR1 corresponding to the first gas reaction amount and a first gas unreacted component Gr1 corresponding to the remaining unreacted amount.

[0190] Next, the in-furnace reaction calculation device 20 uses the first substance reaction amount calculation unit 211-1 to calculate the first substance reaction amount that contributes to the equilibrium state with the modified first gas phase G13 among the first substances M11 existing in the first gas region (first substance reaction amount calculation step: step S222).

[0191] That is, when it is assumed that the first substance M11 existing in the first gas region reacts with the modified first gas phase G13 to reach an equilibrium state, the in-furnace reaction calculation device 20 calculates the first substance reaction amount that contributes to the reaction when the first substance M11 existing in the first gas region reacts with the modified first gas phase G13 to reach an equilibrium state.

[0192] Then, based on the calculation result of the first substance reaction amount, the in-furnace reaction calculation device 20 divides it into a first reaction product substance reaction portion MR11 corresponding to the first substance reaction amount and a first substance unreacted portion Mr11 corresponding to the remaining unreacted amount.

[0193] Next, the in-furnace reaction calculation device 20 provides the first gas reaction portion GR1 obtained in the first gas reaction amount calculation step (step S221) and the first substance reaction portion MR11 obtained in the first substance reaction amount calculation step (step S222) as input substances. The in-furnace reaction calculation device 20 uses the first gas region equilibrium reaction calculation unit 214-1 to calculate the equilibrium reaction between the first gas reaction portion GR1 and the first substance reaction portion MR11, and calculates at least the change in heat quantity and the flow rate of each when the first gas reaction portion GR1 and the first substance reaction portion MR11 reach an equilibrium state (first gas region equilibrium reaction calculation step: step S223).

[0194] For the calculation of the equilibrium reaction, for example, the Gibbs energy minimization method or the like can be used.

[0195] The in-furnace reaction calculation device 20 calculates the gas phase generated by the reaction of the first gas reaction portion GR1 and the first substance reaction portion MR11 and the unused portion of the first gas reaction portion GR1 as the first gas generated component GP1, and calculates the product substance generated by the reaction of the first gas reaction portion GR1 and the first substance reaction portion MR11 and the unused portion of the first substance reaction portion MR11 as the first substance generated component MP11.

[0196] Next, the in-furnace reaction calculation device 20 provides the first gas unreacted portion Gr1 obtained in the first gas reaction amount calculation step (step S221) and the first gas generated component GP1 obtained in the first gas region equilibrium reaction calculation step (step S223) as input substances. The in-furnace reaction calculation device 20 uses the first mixed gas phase calculation unit 216-1 to calculate the flow rate, composition data, etc. of the first mixed gas phase G14 including the first gas unreacted portion Gr1 and the first gas generated component GP1 (first mixed gas phase calculation step: step S224).

[0197] The first gas reaction amount GR1 calculated in the first gas reaction amount calculation step (step S221) is the reaction amount when it is assumed that the first mixed gas phase G13 and the first substance M11 have reached an equilibrium state. Therefore, it is usually all used. However, in the gas-phase first gas reaction amount GR1, there may be substances present in an amount greater than the stoichiometric ratio of the reaction or inert substances that do not contribute to the reaction. Inert substances are, for example, nitrogen and the like. In the first mixed gas phase calculation step (step S223), there are substances present in an amount greater than the stoichiometric ratio of the reaction, and as a result, substances remaining after the equilibrium reaction and inert substances are present. The unused portion remaining without being used in the first gas region equilibrium reaction calculation unit 214-1 is calculated as the gas generation component GP1.

[0198] Next, the in-furnace reaction calculation device 20 provides the unreacted portion Mr11 of the first substance obtained in the first substance reaction amount calculation step (step S222) and the first substance generation component MP11 obtained in the first gas region equilibrium reaction calculation step (step S223) as input substances. The in-furnace reaction calculation device 20 uses the first mixture substance calculation unit 217-1 to calculate the flow rate, composition data, etc. of the first mixture substance M12 obtained by mixing the unreacted portion Mr11 of the first substance and the first substance generation component MP11 (first mixture substance calculation step: step S225).

[0199] Next, the implementation step of the second gas region equilibrium reaction (step S23) will be described. FIG. 7 is a flowchart showing the operation of the implementation step of the second gas region equilibrium reaction (step S23) in FIG. 5. As shown in FIG. 7, in the implementation step of the second gas region equilibrium reaction (step S23), the in-furnace reaction calculation device 20 uses the second gas reaction amount calculation unit 210-2 to calculate the second gas reaction amount that contributes to the equilibrium state with the second substance M21 in the adjusted second gas phase G22B (second gas reaction amount calculation step: step S231).

[0200] That is, the in-furnace reaction calculation device 20 calculates the second gas reaction amount that contributes to the reaction when the adjusted second gas phase G22B reacts with the second substance M21 present in the second gas region and reaches an equilibrium state, assuming that the adjusted second gas phase G22B reacts with the second substance M21 present in the second gas region and reaches an equilibrium state.

[0201] Then, based on the calculation result of the second gas reaction amount, the in-furnace reaction calculation device 20 divides it into a second gas reaction component GR2 corresponding to the second gas reaction amount and a second gas unreacted component Gr2 corresponding to the remaining unreacted amount.

[0202] Next, the in-furnace reaction calculation device 20 uses the second substance reaction amount calculation unit 211-2 to calculate the second substance reaction amount that contributes to the equilibrium state with the adjusted second gas phase G22B among the second substances M21 present in the second gas region (second substance reaction amount calculation step: step S232).

[0203] That is, when assuming that the second substance M21 present in the second gas region reacts with the adjusted second gas phase G22B to reach an equilibrium state, the in-furnace reaction calculation device 20 calculates the second substance reaction amount that contributes to the reaction when the second substance M21 present in the second gas region reacts with the adjusted second gas phase G22B to reach an equilibrium state.

[0204] Then, based on the calculation result of the second substance reaction amount, the in-furnace reaction calculation device 20 divides it into a second reaction product substance reaction component MR21 corresponding to the second substance reaction amount and a second substance unreacted component Mr21 corresponding to the remaining unreacted amount.

[0205] Next, the in-furnace reaction calculation device 20 provides the second gas reaction component GR2 obtained in the second gas reaction amount calculation step (step S231) and the second substance reaction component MR21 obtained in the second substance reaction amount calculation step (step S232) as input substances. The in-furnace reaction calculation device 20 uses the second gas region equilibrium reaction calculation unit 214-2 to calculate the equilibrium reaction between the second gas reaction component GR2 and the second substance reaction component MR21, and calculates at least the changes in the respective heat amounts and the flow rates when the second gas reaction component GR2 and the second substance reaction component MR21 reach an equilibrium state (second gas region equilibrium reaction calculation step: step S233).

[0206] For the calculation of the equilibrium reaction, for example, the Gibbs energy minimization method or the like can be used.

[0207] The in-furnace reaction calculation device 20 calculates the gas phase generated by the reaction of the second gas reaction component GR2 and the second substance reaction component MR21, and the unused portion of the second gas reaction component GR2 as the second gas generation component GP2, and calculates the product substance generated by the reaction of the second gas reaction component GR2 and the second substance reaction component MR21, and the unused portion of the second substance reaction component MR21 as the second substance generation component MP21.

[0208] Next, the in-furnace reaction calculation device 20 provides the second gas unreacted component Gr2 obtained in the second gas reaction amount calculation step (step S231) and the second gas generation component GP2 obtained in the second gas region equilibrium reaction calculation step (step S233) as input substances. The in-furnace reaction calculation device 20 uses the second mixed gas phase calculation unit 216-2 to calculate the flow rate, composition data, etc. of the second mixed gas phase G23 including the second gas unreacted component Gr2 and the second gas generation component GP2 (second mixed gas phase calculation step: step S234).

[0209] The second gas reaction component GR2 calculated in the second gas reaction amount calculation step (step S231) is the reaction amount when it is assumed that the second mixed gas phase G21 and the second substance M21 reach an equilibrium state, so usually all of it is used. However, in the gas-phase second gas reaction component GR2, there may be substances present in an amount greater than the stoichiometric ratio of the reaction or inert substances that do not contribute to the reaction. Inert substances are, for example, nitrogen and the like. In the second mixed gas phase calculation step (step S233), there are substances present in an amount greater than the stoichiometric ratio of the reaction and as a result, substances remaining after the equilibrium reaction and inert substances. The unused portion remaining without being used in the second gas region equilibrium reaction calculation unit 214-2 is calculated as the gas generation component GP2.

[0210] Next, the in-furnace reaction calculation device 20 uses the unreacted amount Mr21 of the second substance obtained in the second substance reaction amount calculation step (step S232) and the generated amount MP21 of the second substance obtained in the second gas region equilibrium reaction calculation step (step S233) as input substances. The in-furnace reaction calculation device 20 uses the second mixed substance calculation unit 217-2 to calculate the flow rate, composition data, etc. of the second mixed substance M22 obtained by mixing the unreacted amount Mr21 of the second substance and the generated amount MP21 of the second substance (second mixed substance calculation step: step S235).

[0211] Next, the implementation step (step S24) of the first bed layer equilibrium reaction will be described. FIG. 8 is a flowchart showing the operation of the implementation step (step S24) of the first bed layer equilibrium reaction in FIG. 5. As shown in FIG. 8, in the implementation step (step S24) of the first bed layer equilibrium reaction, the in-furnace reaction calculation device 20 uses the first bed layer reaction amount calculation unit 212-1 to calculate the first bed layer reaction amount that contributes to the equilibrium state with the first inflowing gas reaction amount GR221 in the modified first bed layer S13B (first bed layer reaction amount calculation step: step S241).

[0212] That is, the in-furnace reaction calculation device 20 calculates the first bed layer reaction amount that contributes to the reaction when the modified first bed layer S13B reacts with the first inflowing gas G221-1 to reach an equilibrium state, assuming that the modified first bed layer S13B and the first inflowing gas G221-1 react to reach an equilibrium state.

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

[0214] Next, the in-furnace reaction calculation device 20 uses the first inflowing gas reaction amount calculation unit 213-1 to calculate the first inflowing gas reaction amount that contributes to the equilibrium state with the modified first bed layer S13B in the first inflowing gas G221-1 (first inflowing gas reaction amount calculation step: step S242).

[0215] That is, when assuming that the first inflow gas G221-1 and the modified first bed layer S13B react to reach an equilibrium state, the in-furnace reaction calculation device 20 calculates the first inflow gas reaction amount contributing to the reaction when the first inflow gas G221-1 reacts with the modified first bed layer S13B to reach an equilibrium state.

[0216] Then, based on the calculation result of the first inflow gas reaction amount, the in-furnace reaction calculation device 20 divides it into a first inflow gas reaction component GR221 corresponding to the first inflow gas reaction amount and a first inflow gas unreacted component Gr221 corresponding to the remaining unreacted amount.

[0217] Next, the in-furnace reaction calculation device 20 provides the first bed layer reaction component SR1 obtained in the first bed layer reaction amount calculation step (step S241) and the first inflow gas reaction component GR221 obtained in the first inflow gas reaction amount calculation step (step S202) as input substances. Using the first bed layer equilibrium reaction calculation unit 215-1, the in-furnace reaction calculation device 20 calculates the equilibrium reaction between the first bed layer reaction component SR1 and the first inflow gas reaction component GR221, and calculates at least the change in heat quantity and the flow rate at the time when the first bed layer reaction component SR1 and the first inflow gas reaction component GR221 reach an equilibrium state (first bed layer equilibrium reaction calculation step S203).

[0218] The in-furnace reaction calculation device 20 calculates the bed layer generated by the reaction between the first bed layer reaction component SR1 and the first inflow gas reaction component GR221 and the unused portion of the first bed layer reaction component SR1 as the first bed layer generated component SP1. Also, the in-furnace reaction calculation device 20 calculates the unused portion of the first inflow gas reaction component GR221 generated by the reaction between the first bed layer reaction component SR1 and the first inflow gas reaction component GR221 as the first inflow gas generated component GP221.

[0219] Next, the in-furnace reaction calculation device 20 provides, as input substances, the unreacted portion Sr1 of the first bed layer obtained in the first bed layer reaction amount calculation step (step S241) and the generated component SP1 of the first bed layer obtained in the first bed layer equilibrium reaction calculation step (step S243). The in-furnace reaction calculation device 20 calculates the flow rate, composition data, etc. of the first bed layer S14 containing the additive component including the unreacted portion Sr1 of the first bed layer and the generated component SP1 of the first bed layer by using the first bed layer equilibrium reaction calculation unit 215-1 (first mixed bed layer calculation step: step S244).

[0220] In the first bed layer reaction amount SR1 of the first bed layer calculated in the first bed layer reaction amount calculation step (step S241), there may be substances present in a stoichiometric ratio or more of the reaction or inert substances that do not contribute to the reaction. In the first mixed bed layer calculation step (step S244), among the first bed layer reaction components SR1, the portions that are present in a stoichiometric ratio or more of the reaction and as a result are not used in the reaction and the inert portions are included in the first bed layer generated component SP1 as the unused portions remaining without being used in the first bed layer equilibrium reaction calculation step (step S243) and calculated together.

[0221] The in-furnace reaction calculation device 20 moves the first mixed bed layer S14 to a region on the discharge end 14B side of region A (region (A + 1)).

[0222] Next, the in-furnace reaction calculation device 20 provides, as input substances, the unreacted portion Gr221 of the first inflow gas obtained in the first inflow gas reaction amount calculation step (step S222) and the generated component GP221 of the first inflow gas obtained in the first bed layer equilibrium reaction calculation step (step S233). The in-furnace reaction calculation device 20 calculates the flow rate, composition data, etc. of the first mixed inflow gas G222A including the unreacted portion Gr221 of the first inflow gas and the generated component GP221 of the first inflow gas by using the first mixed inflow gas calculation unit 219-1 (first mixed inflow gas calculation step: step S205).

[0223] Next, the implementation process of the second bed layer equilibrium reaction (step S21) will be described. FIG. 9 is a flowchart showing the operation of the implementation process of the second bed layer equilibrium reaction (step S21) in FIG. 5. As shown in FIG. 9, in the implementation process of the second bed layer equilibrium reaction (step S25), the in-furnace reaction calculation device 20 uses the second bed layer reaction amount calculation unit 212-2 to calculate the second bed layer reaction amount that contributes to the equilibrium state with the second inflow gas reaction component GR222 in the modified second bed layer S21B (second bed layer reaction amount calculation step: step S251).

[0224] That is, when it is assumed that the modified second bed layer S21B and the second inflow gas G221-2 react to reach an equilibrium state, the in-furnace reaction calculation device 20 calculates the second bed layer reaction amount that contributes to the reaction when the modified second bed layer S21B reacts with the second inflow gas G221-2 to reach an equilibrium state.

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

[0226] Next, the in-furnace reaction calculation device 20 uses the second inflow gas reaction amount calculation unit 213-2 to calculate the second inflow gas reaction amount that contributes to the equilibrium state with the modified second bed layer S21B in the second inflow gas G221-2 (second inflow gas reaction amount calculation step: step S252).

[0227] That is, when it is assumed that the second inflow gas G221-2 and the modified second bed layer S21B react to reach an equilibrium state, the in-furnace reaction calculation device 20 calculates the second inflow gas reaction amount that contributes to the reaction when the second inflow gas G221-2 reacts with the modified second bed layer S21B to reach an equilibrium state.

[0228] Then, based on the calculation result of the second inflow gas reaction amount, the in-furnace reaction calculation device 20 divides it into a second inflow gas reaction component GR222 corresponding to the second inflow gas reaction amount and a second inflow gas unreacted component Gr222 corresponding to the remaining unreacted amount.

[0229] Next, the in-furnace reaction calculation device 20 provides, as input substances, the second bed layer reaction amount SR2 obtained in the second bed layer reaction amount calculation step (step S251) and the second inflow gas reaction amount GR222 obtained in the second inflow gas reaction amount calculation step (step S252). Using the second bed layer equilibrium reaction calculation unit 215-2, the in-furnace reaction calculation device 20 calculates the equilibrium reaction between the second bed layer reaction amount SR2 and the second inflow gas reaction amount GR222, and calculates at least the change in heat quantity and the flow rate at the time when the second bed layer reaction amount SR2 and the second inflow gas reaction amount GR222 reach an equilibrium state (second bed layer equilibrium reaction calculation step: step S253).

[0230] The in-furnace reaction calculation device 20 calculates, as the second bed layer generated component SP2, the bed layer generated by the reaction between the second bed layer reaction amount SR2 and the second inflow gas reaction amount GR222, and the unused portion of the second bed layer reaction amount SR2. Further, the in-furnace reaction calculation device 20 calculates, as the second inflow gas generated component GP222, the unused portion of the second inflow gas reaction amount GR222 generated by the reaction between the second bed layer reaction amount SR2 and the second inflow gas reaction amount GR222.

[0231] Next, the in-furnace reaction calculation device 20 provides, as input substances, the second bed layer unreacted amount Sr2 obtained in the second bed layer reaction amount calculation step (step S251) and the second bed layer generated component SP2 obtained in the second bed layer equilibrium reaction calculation step (step S253). Using the second bed layer equilibrium reaction calculation unit 215-2, the in-furnace reaction calculation device 20 calculates the flow rate, composition data, etc. of the second mixed bed layer S22 including the second bed layer unreacted amount Sr2 and the second bed layer generated component SP2 (second mixed bed layer calculation step: step S254).

[0232] In the second bed layer reaction amount calculation step (step S251), in the second bed layer reaction component SR2 of the second bed layer calculated, there may be substances present in a stoichiometric ratio or more of the reaction and inert substances that do not contribute to the reaction. In the second mixed bed layer calculation step (step S254), among the second bed layer reaction component SR2, there are portions that are present in a stoichiometric ratio or more of the reaction and as a result are not used in the reaction and inert portions. The unused portion remaining without being used in the second bed layer equilibrium reaction calculation step (step S253) is calculated as the second bed layer generated component SP2.

[0233] The in-furnace reaction calculation device 20 moves the second mixed bed layer S22 to a region on the discharge end 14B side of region A (region (A + 1)).

[0234] Next, the in-furnace reaction calculation device 20 gives the unreacted portion Gr222 of the second inflowing gas obtained in the second inflowing gas reaction amount calculation step (step S252) and the generated component GP222 of the second inflowing gas obtained in the second bed layer equilibrium reaction calculation step (step S253) as input substances. The in-furnace reaction calculation device 20 calculates the flow rate, composition data, etc. of the second mixed inflowing gas G222B including the unreacted portion Gr222 of the second inflowing gas and the generated component GP222 of the second inflowing gas using the second mixed inflowing gas calculation unit 219-2 (second mixed inflowing gas calculation step: step S255).

[0235] Note that in the in-furnace reaction calculation method according to the present embodiment, each step may be performed in parallel as appropriate, or the order of the steps may be changed.

[0236] For example, in the in-furnace reaction calculation method according to the present embodiment, at least one or more of the steps of the first gas region equilibrium reaction implementation step (step S22), the second gas region equilibrium reaction implementation step (step S23), the first bed layer equilibrium reaction implementation step (step S24), and the second bed layer equilibrium reaction implementation step (step S25) shown in FIG. 5 may be performed in parallel.

[0237] The calculation method of the in-furnace reaction according to this embodiment may perform the inflow gas summation step (step S26) and the first mixed substance amount correction step (step S28) shown in FIG. 5 in parallel, or may perform the inflow gas summation step (step S26) after the first mixed substance amount correction step (step S28).

[0238] The calculation method of the in-furnace reaction according to this embodiment may perform the gas amount summation step (step S27) and the first mixed substance amount correction step (step S28) shown in FIG. 5 in parallel, or may perform the gas amount summation step (step S27) after the first mixed substance amount correction step (step S28).

[0239] The calculation method of the in-furnace reaction according to this embodiment may perform the first gas reaction amount calculation step (step S221) and the first substance reaction amount calculation step (step S222) shown in FIG. 6 in parallel, or may perform the first gas reaction amount calculation step (step S221) after the first substance reaction amount calculation step (step S222).

[0240] The calculation method of the in-furnace reaction according to this embodiment may perform the first mixed gas phase calculation step (step S224) and the first mixed substance calculation step (step S225) shown in FIG. 6 in parallel, or may perform the first mixed gas phase calculation step (step S224) after the first mixed substance calculation step (step S225).

[0241] The calculation method of the in-furnace reaction according to this embodiment may perform the second gas reaction amount calculation step (step S231) and the second substance reaction amount calculation step (step S232) shown in FIG. 7 in parallel, or may perform the second gas reaction amount calculation step (step S231) after the second substance reaction amount calculation step (step S232).

[0242] The calculation method of the in-furnace reaction according to this embodiment may perform the second mixed gas phase calculation step (step S234) and the second mixed substance calculation step (step S235) shown in FIG. 7 in parallel, or may perform the second mixed gas phase calculation step (step S234) after the second mixed substance calculation step (step S235).

[0243] The calculation method of the in-furnace reaction according to this embodiment may perform the first bed layer reaction amount calculation step (step S241) and the first inflowing gas reaction amount calculation step (step S242) shown in FIG. 8 in parallel, or may perform the first bed layer reaction amount calculation step (step S241) after the first inflowing gas reaction amount calculation step (step S242).

[0244] The calculation method of the in-furnace reaction according to this embodiment may perform the first mixed bed layer calculation step (step S244) and the first mixed inflowing gas calculation step (step S245) shown in FIG. 8 in parallel, or may perform the first mixed bed layer calculation step (step S244) after the first mixed inflowing gas calculation step (step S245).

[0245] The calculation method of the in-furnace reaction according to this embodiment may perform the second bed layer reaction amount calculation step (step S251) and the second inflowing gas reaction amount calculation step (step S252) shown in FIG. 9 in parallel, or may perform the second bed layer reaction amount calculation step (step S251) after the second inflowing gas reaction amount calculation step (step S252).

[0246] The calculation method of the in-furnace reaction according to this embodiment may perform the second mixed bed layer calculation step (step S254) and the second mixed inflowing gas calculation step (step S255) shown in FIG. 9 in parallel, or may perform the second mixed bed layer calculation step (step S254) after the second mixed inflowing gas calculation step (step S255).

[0247] The calculation method of the in-furnace reaction according to this embodiment may calculate the heat conduction in the furnace using models such as radiation, conduction, and convection as necessary.

[0248] <Hardware Configuration of In-Furnace Reaction Calculation Device> Next, an example of the hardware configuration of the in-furnace reaction calculation device will be described. FIG. 10 is a hardware configuration diagram of the in-furnace reaction calculation device. As shown in FIG. 10, the in-furnace reaction calculation device 20 is composed of, for example, an information processing device (computer). Physically, it is 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, etc. These are interconnected by a bus 27. Note that the auxiliary storage device 24 and the display device 26 may be provided externally.

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

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

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

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

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

[0254] The display device 26 is a monitor display or the like. On the display device 26, a measurement recording screen and an analysis screen are displayed, and the screen is updated according to input / output operations via the input / output interface 25.

[0255] Each function of the in-furnace reaction calculation device 20 shown in FIG. 10 is realized by causing a main storage device or an auxiliary storage device 24 such as a RAM 22 or a ROM 23 to load simulation software (including an in-furnace reaction calculation program), etc., and causing the CPU 21 to execute a reaction calculation program for raw ore stored in the RAM 22, ROM 23, or auxiliary storage device 24, thereby reading and writing data in the RAM 22, etc., and operating the input / output interface 25 and the display device 26.

[0256] As the in-furnace reaction calculation program, a program having the following configuration can be used. That is, the in-furnace reaction calculation program is a program for causing at least a computer to execute a calculation of an in-furnace reaction in which raw ore supplied from one end side of a reactor is moved toward the other end side while being brought into contact with combustion gas supplied from the other end side to be dried and reduced, a first gas-phase material amount correction step of correcting the material amount of a first gas phase flowing through a first gas region, which is one of two gas regions distributed according to a gas region influence factor including product gas generated from a bed layer containing the raw ore, among the gas regions containing the combustion gas, and obtaining a corrected first gas phase; a second gas-phase material amount correction step of correcting the material amount of a second gas phase flowing through a second gas region, which is the other of the two gas regions, and obtaining a corrected second gas phase; a gas mixing amount calculation step of calculating the flow rate of the inflow gas moving from the second gas region to the bed layer; a first gas region equilibrium reaction calculation step of calculating an equilibrium reaction between the corrected first gas phase and a first substance existing in the first gas region among substances including at least one of a solid substance and a liquid substance existing in the gas region; A second gas region equilibrium reaction calculation step of calculating an equilibrium reaction between the corrected second gas phase and a second substance among the substances present in the second gas region; A bed layer equilibrium reaction calculation step of calculating an equilibrium reaction between the bed layer and the inflowing gas; A program for causing at least a computer to execute the above steps can be used.

[0257] The calculation program for the in-furnace reaction is stored, for example, in a storage device included in a computer such as the main storage device or auxiliary storage device 24 of RAM 22 or ROM 23. Note that the reaction calculation program for the raw ore may be partially or entirely transmitted via a transmission medium such as a communication line and received and recorded (including installation) by a communication module or the like included in the computer. Further, the reaction calculation program for the raw ore may be partially or entirely recorded (including installation) in the computer from a state stored in a portable storage medium such as a CD-ROM, DVD-ROM, or flash memory.

[0258] As described above, the in-furnace reaction calculation device 20 according to the present embodiment includes a first gas phase substance amount correction unit 202, a second gas phase substance amount correction unit 203, a gas mixing amount calculation unit 204, a first gas region equilibrium reaction calculation unit 214-1, and a second gas region equilibrium reaction calculation unit 214-2. The in-furnace reaction calculation device 20 corrects the substance amounts of the first gas phase G12 and the second gas phase G21 containing additives by the first gas phase substance amount correction unit 202 and the second gas phase substance amount correction unit 203, respectively, and then the gas mixing amount calculation unit 204 calculates the flow rate of the inflowing gas G221 that moves from the corrected second gas phase G22A in the second gas region to the bed layer. Then, the in-furnace reaction calculation device 20 calculates an equilibrium reaction between the corrected first gas phase G13 and the first substance M11 in the first gas region equilibrium reaction calculation unit 214-1, and calculates an equilibrium reaction between the adjusted second gas phase G22B and the second substance M21 in the second gas region equilibrium reaction calculation unit 214-2.

[0259] The first gas phase G11 present in the first gas region is a gas phase that contains almost no reactive gas, dust, etc. generated from the bed layer even within the gas region. The second gas phase G21 present in the second gas region is a gas phase that contains a large amount of reactive gas, dust, etc. even within the gas region. The in-furnace reaction calculation device 20 divides the gas region into two regions (the first gas region and the second gas region) according to the proportion affected by reactive gas, dust, etc. generated from the bed layer, and performs equilibrium reaction calculations for each region, so that the equilibrium reaction can be calculated with higher accuracy according to the characteristics of each gas region.

[0260] Therefore, the in-furnace reaction calculation device 20 divides the gas region in the reactor into two series, namely the first gas region and the second gas region, and calculates the equilibrium reaction between the modified first gas phase G13 flowing in each region and the first substance M11, and the equilibrium reaction between the adjusted second gas phase G22B and the second substance M21 respectively based on the Gibbs energy minimization method. Thereby, the in-furnace reaction calculation device 20 can calculate the equilibrium reaction occurring in the gas region in consideration of the influence of the post-reaction gas, dust, etc. generated from the bed layer in the two types of gas regions in the rotary kiln 1. Also, without requiring huge computing resources, it can be considered when calculating the influence of the post-reaction gas etc. from the bed layer that exists at a high concentration near the bed layer. Therefore, the in-furnace reaction calculation device 20 can calculate the equilibrium reaction in the gas region with high accuracy without imposing a large burden.

[0261] In addition, the in-furnace reaction calculation device 20 can calculate the equilibrium reaction between the bed layer and the inflow gas G221 from the second gas region in the bed layer equilibrium reaction calculation unit 215, thereby calculating the equilibrium reaction between the bed layer and the inflow gas G221 mixed from the second gas region. Thereby, the in-furnace reaction calculation device 20 can calculate the equilibrium reaction between the bed layer and the inflow gas G221 from the second gas region based on the Gibbs energy minimization method. Therefore, the in-furnace reaction calculation device 20 can calculate the equilibrium reaction occurring in the bed layer with high accuracy in consideration of the influence of the post-reaction gas in the bed layer (for example, the gas that has flowed into the bed layer from the gas region in the bed layer or the reaction gas generated when the bed layer reacts in equilibrium) etc. that exists at a high concentration near the bed layer in the rotary kiln 1.

[0262] Therefore, the in-furnace reaction calculation device 20 can accurately calculate the reactions occurring in the rotary kiln 1. Since the in-furnace reaction calculation device 20 can calculate the in-furnace reactions occurring in the rotary kiln 1 in a state close to the actual furnace while reducing the burden required for the calculation, it is possible to perform an analysis closer to the actual furnace state with a low load.

[0263] The in-furnace reaction calculation device 20 can include a first gas reaction amount calculation unit 210-1 and a first substance reaction amount calculation unit 211-1. Thereby, the in-furnace reaction calculation device 20 can obtain the first gas reaction component GR1 contributing to the equilibrium reaction of the modified first gas phase G13 in the first gas reaction amount calculation unit 210-1, and can obtain the first substance reaction component MR11 contributing to the equilibrium reaction of the first substance M11 in the first substance reaction amount calculation unit 211-1. Therefore, the in-furnace reaction calculation device 20 can surely perform the equilibrium reaction calculation between the modified first gas phase G13 and the first substance M11 using the modified first gas phase G13 and the first substance M11 required for the reaction in the first gas region equilibrium reaction calculation unit 214-1.

[0264] The in-furnace reaction calculation device 20 can include a first mixed gas phase calculation unit 216-1 and a first mixture substance calculation unit 217-1. Thereby, the in-furnace reaction calculation device 20 can calculate at least the flow rate of the first mixed gas phase G14 obtained by mixing the first unreacted gas component Gr1 and the first generated gas component GP1 in the first mixed gas phase calculation unit 216-1. The in-furnace reaction calculation device 20 can calculate at least the flow rate of the first mixture substance M12 obtained by mixing the first unreacted substance component Mr11 and the first generated substance component MP11 in the first mixture substance calculation unit 217-1. Thereby, the in-furnace reaction calculation device 20 can more accurately calculate the amount of movement of the first gas phase G11 and the first substance M11 to the charging end 14A.

[0265] The in-furnace reaction calculation device 20 can include a second gas reaction amount calculation unit 210-2 and a second substance reaction amount calculation unit 211-2. Thereby, in the second gas reaction amount calculation unit 210-2 of the in-furnace reaction calculation device 20, the second gas reaction component GR2 contributing to the equilibrium reaction of the adjusted second gas phase G22B can be obtained, and in the second substance reaction amount calculation unit 211-2, the first substance reaction component MR21 contributing to the equilibrium reaction of the second substance M21 can be obtained. Therefore, the in-furnace reaction calculation device 20 can surely perform the equilibrium reaction calculation between the adjusted second gas phase G22B and the second substance M21 using the adjusted second gas phase G22B and the second substance M21 required for the reaction in the second gas region equilibrium reaction calculation unit 214-2.

[0266] The in-furnace reaction calculation device 20 can include a second mixed gas phase calculation unit 216-2 and a second mixture substance calculation unit 217-2. Thereby, the in-furnace reaction calculation device 20 can at least calculate the flow rate of the second mixed gas phase G23 obtained by mixing the unreacted gas component Gr2 and the generated gas component GP2 in the second mixed gas phase calculation unit 216-2. The in-furnace reaction calculation device 20 can at least calculate the flow rate of the second mixture substance M22 obtained by mixing the unreacted second substance component Mr21 and the generated second substance component MP21 in the second mixture substance calculation unit 217-2. Thereby, the in-furnace reaction calculation device 20 can more accurately calculate the movement amount of the second gas phase G21 and the second substance M21 to the charging end 14A.

[0267] The in-furnace reaction calculation device 20 can include a first mixture substance amount correction unit 220. Thereby, the in-furnace reaction calculation device 20 excludes the movement amount of the fourth substance M4 contained in the first mixture substance M12 present in the gas region to the second gas region, and adds the movement amount of the third substance M3 present in the second corrected mixture substance M23 in the second gas region to the first mixture substance M12 to calculate the first corrected mixture substance M13. Therefore, the in-furnace reaction calculation device 20 can accurately calculate the movement amount of the first substance M11 that has entered region A to region (A-1).

[0268] The in-furnace reaction calculation device 20 includes an inflow gas distribution unit 209, and the bed layer equilibrium reaction calculation unit 215 can have a first bed layer equilibrium reaction calculation unit 215-1 and a second bed layer equilibrium reaction calculation unit 215-2. The in-furnace reaction calculation device 20 distributes the inflow gas G221 into a first inflow gas G221-1 and a second inflow gas G221-2 in the inflow gas distribution unit 209. Then, the in-furnace reaction calculation device 20 calculates the equilibrium reaction between the modified first bed layer S13B and the first inflow gas G221-1 in the first bed layer equilibrium reaction calculation unit 215-1, and can calculate the equilibrium reaction between the modified second bed layer S21B and the second inflow gas G221-2 in the second bed layer equilibrium reaction calculation unit 215-2. Thereby, since the in-furnace reaction calculation device 20 can calculate the equilibrium reaction between the first bed layer and the second bed layer generated in the bed layer and the inflow gas G221, the equilibrium reaction between the bed layer in the rotary kiln 1 and the inflow gas G221 can be calculated with higher accuracy.

[0269] The in-furnace reaction calculation device 20 can include a first bed layer reaction amount calculation unit 212-1 and a first inflow gas reaction amount calculation unit 213-1. Thereby, the in-furnace reaction calculation device 20 obtains a first bed layer reaction component SR1 contributing to the equilibrium reaction from the modified first bed layer S13B in the first bed layer reaction amount calculation unit 212-1, and can obtain a first inflow gas reaction component GR221 contributing to the equilibrium reaction from the first inflow gas G221-1 in the first inflow gas reaction amount calculation unit 213-1. Therefore, the in-furnace reaction calculation device 20 can surely perform the equilibrium reaction calculation between the modified first bed layer S13B and the first inflow gas G221-1 using the first bed layer reaction component SR1 and the first inflow gas reaction component GR221 required for the reaction in the first bed layer equilibrium reaction calculation unit 215-1.

[0270] The in-furnace reaction calculation device 20 can include a first mixed bed layer calculation unit 218-1 and a first mixed inflow gas calculation unit 219-1. Thereby, the in-furnace reaction calculation device 20 can calculate at least the flow rate of the first mixed bed layer S14 obtained by mixing the unreacted Sr1 in the first bed layer and the first bed layer generated component SP1 in the first mixed bed layer calculation unit 218-1. The in-furnace reaction calculation device 20 can calculate at least the flow rate of the first mixed inflow gas G222A obtained by mixing the unreacted Gr221 in the first inflow gas and the first inflow gas generated component GP221 in the first mixed inflow gas calculation unit 219-1. Thereby, the in-furnace reaction calculation device 20 can accurately calculate the movement amount to the discharge end 14B of the first mixed bed layer S14 and can accurately calculate the movement amount of the first inflow gas G221-1.

[0271] The in-furnace reaction calculation device 20 can include a second bed layer reaction amount calculation unit 212-2 and a second inflow gas reaction amount calculation unit 213-2. Thereby, the in-furnace reaction calculation device 20 can calculate the second bed layer reaction component SR2 contributing to the equilibrium reaction from the modified second bed layer S21B in the second bed layer reaction amount calculation unit 212-2. Further, the in-furnace reaction calculation device 20 can obtain the second inflow gas reaction component GR222 contributing to the equilibrium reaction from the second inflow gas G221-2 in the second inflow gas reaction amount calculation unit 213-2. Therefore, the in-furnace reaction calculation device 20 can surely perform the equilibrium reaction calculation between the modified second bed layer S21B and the inflow gas G221-2 using the second bed layer reaction component SR2 and the second inflow gas reaction component GR222 required for the reaction in the second bed layer equilibrium reaction calculation unit 215-2.

[0272] The in-furnace reaction calculation device 20 can include a second mixed bed layer calculation unit 218-2 and a second mixed inflow gas calculation unit 219-2. Thereby, the in-furnace reaction calculation device 20 can at least calculate the flow rate of the second mixed bed layer S22 obtained by mixing the unreacted Sr2 in the second bed layer and the second bed layer generated component SP2 in the second mixed bed layer calculation unit 218-2. Further, the in-furnace reaction calculation device 20 can at least calculate the flow rate of the second mixed inflow gas G222B obtained by mixing the unreacted Gr222 in the second inflow gas and the second inflow gas generated component GP222 in the second mixed inflow gas calculation unit 219-2. Thereby, the in-furnace reaction calculation device 20 can accurately calculate the amount of movement to the discharge end 14B of the second mixed bed layer S22 and can accurately calculate the amount of movement of the second inflow gas G221-2.

[0273] The in-furnace reaction calculation device 20 can include a first bed layer preliminary material quantity correction unit 205. Thereby, the in-furnace reaction calculation device 20 reduces the amount of movement of the seventh substance M7 in the second gas region of the first bed layer S12 containing the additive component to the second mixed substance M22 and increases the amount of movement of the eighth substance M8 contained in the second mixed substance M22 to the first bed layer S12 containing the additive component in the first bed layer preliminary material quantity correction unit 205, so that the material quantity of the first bed layer S12 containing the additive component can be corrected to calculate the preliminary corrected first bed layer S13A. Therefore, the in-furnace reaction calculation device 20 can accurately obtain the first bed layer reaction component SR1 used in the first bed layer equilibrium reaction calculation unit 215-1. Thus, the in-furnace reaction calculation device 20 can perform the first bed layer equilibrium reaction calculation in the first bed layer equilibrium reaction calculation unit 215-1 more accurately.

[0274] The in-furnace reaction calculation device 20 can include a second bed layer preliminary material quantity correction unit 207. Thereby, in the second bed layer preliminary material quantity correction unit 207 of the in-furnace reaction calculation device 20, the movement amount of the ninth substance M9 contained in the second bed layer S21 to the second mixed substance M22 in the second gas region is reduced, the material quantity of the second bed layer S21 is corrected, and a preliminary corrected second bed layer S21A can be calculated. Therefore, the in-furnace reaction calculation device 20 can accurately obtain the flow rate, etc. of the preliminary corrected second bed layer S21A used in the calculation of the equilibrium reaction in the second bed layer equilibrium reaction calculation unit 215-2. Therefore, the in-furnace reaction calculation device 20 can further improve the calculation accuracy of the equilibrium reaction of the preliminary corrected second bed layer S21A in the second bed layer equilibrium reaction calculation unit 215-2.

[0275] The in-furnace reaction calculation device 20 can include an inflow gas summation unit 221. Thereby, the in-furnace reaction calculation device 20 can calculate a mixed inflow gas G223 obtained by summing the first mixed inflow gas G222A and the second mixed inflow gas G222B. Therefore, the in-furnace reaction calculation device 20 can move the mixed inflow gas G223 present in the bed layer after the equilibrium reaction to the second gas region. Therefore, the in-furnace reaction calculation device 20 can calculate the movement amount of the mixed inflow gas G223 present in the bed layer after the equilibrium reaction to the second mixed gas phase G23 present in the second gas region.

[0276] The in-furnace reaction calculation device 20 can include a gas quantity summation unit 222. Thereby, the in-furnace reaction calculation device 20 can mix the mixed inflow gas G223 generated in the bed layer with the second gas phase G23 to obtain a second summed gas phase G24. Therefore, the in-furnace reaction calculation device 20 can move the mixed inflow gas G223 that has moved from the bed layer in region A to the second gas region and include it in the second gas phase G23 and move it to region (A-1). Therefore, the in-furnace reaction calculation device 20 can calculate the movement amount of the second summed gas phase G24 including the mixed inflow gas G223 and the second mixed gas phase G23 to region (A-1).

[0277] The in-furnace reaction calculation device 20 can include a second mixture primary material quantity correction unit 223. Thereby, in the second mixture primary material quantity correction unit 223 of the in-furnace reaction calculation device 20, the movement amount of the added amount of the eighth substance M8 contained in the second mixture M22 generated in the second mixture calculation unit 217-2 to the first bed layer S12 containing the added amount is reduced, and the movement amount of the seventh substance M7 contained in the first bed layer S12 containing the added amount generated in the first bed layer preliminary material quantity correction unit 205 to the second mixture M22 and the movement amount of the ninth substance M9 contained in the second mixed bed layer S21 generated in the second bed layer preliminary material quantity correction unit 207 to the second mixture M22 can be added. As a result, the in-furnace reaction calculation device 20 can correct the material quantity of the second mixture M22 and calculate the second corrected mixture M23. Therefore, the in-furnace reaction calculation device 20 can calculate more accurately the movement amount of the second substance M21 in the second gas region after the equilibrium reaction to the charging end 14A.

[0278] The in-furnace reaction calculation device 20 can include a second mixture secondary material quantity correction unit 224. Thereby, in the second mixture secondary material quantity correction unit 224 of the in-furnace reaction calculation device 20, the movement amount of the third substance M3 contained in the second corrected mixture M23 to the first mixture M12 is reduced, and the movement amount of the fourth substance M4 contained in the first mixture M12 can be added. As a result, the in-furnace reaction calculation device 20 can correct the material quantity of the second corrected mixture M23 and calculate the second corrected mixture M24. Therefore, the in-furnace reaction calculation device 20 can calculate even more accurately the movement amount of the second substance M21 in the second gas region after the equilibrium reaction to the charging end 14A.

[0279] The in-furnace reaction calculation device 20 can include a gas region distribution unit 231. Thereby, the in-furnace reaction calculation device 20 can distribute the gas region into a first gas region and a second gas region in advance according to the gas region influence factors contained in the gas region immediately after the combustion gas is supplied from the discharge end 14B of the rotary kiln 1. The in-furnace reaction calculation device 20 can accurately calculate the masses of the first gas phase G11 flowing through the first gas region and the second gas phase G21 flowing through the second gas region from the region N including the discharge end 14B until reaching the region A. Therefore, the in-furnace reaction calculation device 20 can more accurately calculate the mass flow rates at which the first gas phase G11 and the second gas phase G21 are used in the equilibrium reaction in the region A. Thus, the in-furnace reaction calculation device 20 can more accurately calculate the respective equilibrium reactions in the first gas region and the second gas region.

[0280] The in-furnace reaction calculation device 20 can include a first substance distribution unit 232. Thereby, the in-furnace reaction calculation device 20 can distribute the first substance contained in the gas region into a first substance M11 and a second substance M21 in advance from immediately after the combustion gas is supplied from the discharge end 14B of the rotary kiln 1. The in-furnace reaction calculation device 20 can accurately calculate the masses of the first substance M11 and the second substance M21 flowing through the gas region from the region N including the discharge end 14B until reaching the region A. Therefore, the in-furnace reaction calculation device 20 can more accurately calculate the mass flow rates at which the first substance M11 and the second substance M21 are used in the equilibrium reaction in the region A. Thus, the in-furnace reaction calculation device 20 can more accurately calculate the respective equilibrium reactions in the first gas region and the second gas region.

[0281] The in-furnace reaction calculation device 20 includes a distribution unit 201. In the distribution unit 201, the gas phase and the bed layer contained in the combustion material 240 introduced during the movement of the raw ore can be distributed by mass flow into an added gas phase AG11 and an added bed layer AS11. The in-furnace reaction calculation device 20 distributes the combustion material 240 introduced from the middle of the rotary kiln 1 into a first gas phase G1 and a first bed layer S1 and treats them as separate fuels respectively. The added gas phase AG11 can be a component constituting the first gas phase G11, and the added bed layer AS11 can be a component constituting the first bed layer S11. The in-furnace reaction calculation device 20 can obtain a first gas reaction component GR1 from the first gas phase G11 and a first bed layer reaction component SR1 from the first bed layer S11. Thereby, the in-furnace reaction calculation device 20 can calculate an equilibrium reaction based on the reaction components required for the equilibrium reaction from the modified first gas phase G13 obtained from the first gas phase G12 containing the added component and the first substance M11 in the first gas region equilibrium reaction calculation unit 214-1. Also, the in-furnace reaction calculation device 20 can calculate an equilibrium reaction based on the reaction components required for the equilibrium reaction from the first inflow gas G221-1 and the modified first bed layer S13B in the first bed layer equilibrium reaction calculation unit 215-1. Therefore, since the in-furnace reaction calculation device 20 can analyze the behavior of substances in the rotary kiln 1 in consideration of the combustion material 240 introduced from the middle of the rotary kiln 1, even when the combustion material 240 is used, the in-furnace reaction can be calculated with high accuracy.

[0282] The in-furnace reaction calculation device 20 can include a first bed layer material quantity correction unit 206. That is, the in-furnace reaction calculation device 20 considers the transfer amount of the sixth substance M6 in the preliminary modified first bed layer S13A, whose material quantity has been corrected by the first bed layer preliminary material quantity correction unit 205, to the preliminary modified second bed layer S21A, and the transfer amount of the fifth substance M5 in the preliminary modified second bed layer S21A to the preliminary modified first bed layer S13A, corrects the material quantity of the preliminary modified first bed layer S13A, and can calculate the modified first bed layer S13B. At this time, the first bed layer equilibrium reaction calculation unit 215-1 uses the modified first bed layer S13B as the first bed layer S11.

[0283] The in-furnace reaction calculation device 20 can, in the first bed layer material quantity correction unit 206, exclude the amount of the sixth substance M6 present in the preliminary correction first bed layer S13A that moves to the second bed layer S21A, and add the amount of the fifth substance M5 present in the second bed S21A that moves to the preliminary correction first bed layer S13A. Therefore, the in-furnace reaction calculation device 20 can more accurately calculate the flow rate and the like of the corrected first bed layer S13B used in the calculation of the equilibrium reaction in the first bed layer equilibrium reaction calculation unit 215-1, and thus can improve the calculation accuracy of the equilibrium reaction of the corrected first bed layer S13B in the first bed layer.

[0284] The in-furnace reaction calculation device 20 can include a second bed layer material quantity correction unit 208. That is, the in-furnace reaction calculation device 20 considers the amount of the fifth substance M5 included in the preliminary correction second bed layer S21A whose material quantity has been corrected by the second bed layer preliminary material quantity correction unit 207 that moves to the preliminary correction first bed layer S13A, and the amount of the sixth substance M6 included in the preliminary correction first bed layer S13A that moves to the preliminary correction second bed layer S21A, corrects the material quantity of the preliminary correction second bed layer S21A, and can calculate the corrected second bed layer S21B. At this time, the second bed layer equilibrium reaction calculation unit 215-2 uses the corrected second bed layer S21B as the second bed layer S21.

[0285] The in-furnace reaction calculation device 20 can, in the second bed layer material quantity correction unit 208, exclude the amount of the fifth substance M5 present in the second bed layer S21A that moves to the preliminary correction first bed layer S13A, and add the amount of the sixth substance M6 present in the preliminary correction first bed layer S13A that moves to the second bed layer S21A. Therefore, the in-furnace reaction calculation device 20 can accurately obtain the flow rate and the like of the preliminary correction second bed layer S21A used in the calculation of the equilibrium reaction in the second bed layer equilibrium reaction calculation unit 215-2, and thus can improve the calculation accuracy of the equilibrium reaction of the corrected second bed layer S21B in the second bed layer.

[0286] The in-furnace reaction calculation device 20 can include the above-described first bed layer material quantity correction unit 206 and second bed layer material quantity correction unit 208. Thereby, the in-furnace reaction calculation device 20 can more accurately calculate the flow rate and the like of the corrected first bed layer S13B used in the calculation of the equilibrium reaction in the first bed layer equilibrium reaction calculation unit 215-1, and can accurately obtain the flow rate and the like of the corrected second bed layer S21B used in the calculation of the equilibrium reaction in the second bed layer equilibrium reaction calculation unit 215-2. Therefore, the in-furnace reaction calculation device 20 can further improve the calculation accuracy of the equilibrium reaction of the corrected first bed layer S13B in the first bed layer S11 and the equilibrium reaction of the corrected second bed layer S21B in the second bed layer S21, respectively.

[0287] The in-furnace reaction calculation device 20 can include a bed layer distribution unit 233. Thereby, the in-furnace reaction calculation device 20 can divide the bed layer S1 into a first bed layer S11 and a second bed layer S21 immediately after being supplied from the charging end 14A of the rotary kiln 1. Since the in-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 until reaching the region A, it can more accurately calculate the mass flow rates at which the first bed layer S11 and the second bed layer S21 are used in the equilibrium reaction in the region A. Therefore, the in-furnace reaction calculation device 20 can more accurately calculate the first bed layer equilibrium reaction calculation and the second bed layer equilibrium reaction calculation.

[0288] In this way, considering that the gas region in the rotary kiln 1 exhibits a plurality of different reaction behaviors, the in-furnace reaction calculation device 20 can calculate the equilibrium reaction occurring in the furnace with high precision. Therefore, by applying the in-furnace reaction calculation device 20 to all regions in the rotary kiln 1, it becomes possible to calculate the entire reaction process in the rotary kiln 1 with high precision.

[0289] A case where the in-furnace reaction calculation device 20 is applied to the entire rotary kiln 1 will be described. As shown in FIG. 11, for example, when the in-furnace reaction calculation device 20 assumes the reaction process in one region A when the inside of the rotary kiln 1 is divided into a plurality of regions as a unit operation model, the reaction process in the rotary kiln 1 can be modeled by a combination of unit operation models. Then, the in-furnace reaction calculation device 20 repeatedly performs (see the arrows in FIG. 11) along the flow of the raw ore or the flow of the combustion gas (including the assumed flow) in a plurality of regions in the rotary kiln 1. And the in-furnace reaction calculation device 20 repeatedly performs until the difference between the calculated value in a predetermined region and the previous calculated value in that region falls within a predetermined range.

[0290] A flowchart of the case where the in-furnace reaction calculation device 20 is applied to the entire rotary kiln 1 is shown in FIG. 12. As shown in FIG. 12, the in-furnace reaction calculation device 20 models the reaction process occurring in the rotary kiln 1 by a combination of a plurality of unit operation models (modeling step: step S21).

[0291] The unit operation model is applied to the in-furnace reaction calculation device 20 shown in FIG. 2 above. For each unit operation model, a gas mixing amount calculation unit 204, an inflow gas distribution unit 209, a first gas region equilibrium reaction calculation unit 214-1, a second gas region equilibrium reaction calculation unit 214-2, a first bed layer equilibrium reaction calculation unit 215-1, a second bed layer equilibrium reaction calculation unit 215-2, etc., which constitute the unit operation model, are prepared in advance.

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

[0293] Next, the in-furnace reaction calculation device 20 performs calculations on the unit operation models modeled in step S21 (calculation step: step S22). Information on the flow is input to the unit operation model.

[0294] The flow information is data such as the components, flow rates, temperatures, rotational speeds, etc. of the raw ore, combustion gas, and combustion materials. When the flow information is input, the unit operation model performs predetermined calculations, and the calculated values (components, flow rates, temperatures, etc. of the raw ore, combustion gas, and combustion materials) are output. From these calculation results, the values used for each component of the in-furnace reaction calculation device 20 shown in FIG. 2 are calculated. The values used for each part of the in-furnace reaction calculation device 20 include the flow rates of various gas phases in region A such as the bed layer S1 and the first gas phase G1, the flow rates of various bed layers in region A such as the first bed layer S1, the flow rates of various substances such as the flow rate of the first substance M1, the gas reaction fraction GR1, the unreacted gas fraction Gr1, the generated gas fraction GP1, etc. in region A, and the flow rates of various reaction fractions, unreacted fractions, and generated fractions.

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

[0296] Next, the in-furnace reaction calculation device 20 determines whether it has calculated up to the final unit operation model (step S23).

[0297] If it has calculated up to the final unit operation model (step 23: Yes), the in-furnace reaction calculation device 20 determines whether there is a previous calculated value of the unit operation model (step S24).

[0298] If there is a previous calculated value (step S24: Yes), the in-furnace reaction calculation device 20 compares the calculated value calculated in the calculation process (step S22) with the previous calculated value (step S25).

[0299] Next, the in-furnace reaction calculation device 20 determines whether the difference between the calculated value and the previous calculated value satisfies the convergence condition (comparison step: step S26).

[0300] As the convergence condition, for example, the difference between the calculated value and the previous calculated value is within a range of several degrees Celsius (for example, 1°C) or less.

[0301] When the difference between the calculated value and the previous calculated value satisfies the convergence condition (step S26: Yes), the in-furnace reaction calculation device 20 ends the calculation. As a result, the reaction process occurring throughout the kiln body 11 of the rotary kiln 1 is analyzed.

[0302] On the other hand, in step S23, when the final unit operation model has not been calculated (step S23: No), the in-furnace reaction calculation device 20 shifts to the unit operation model located in the region (A + 1) or region (A - 1), which is another adjacent unit operation model (step S27). Then, the in-furnace reaction calculation device 20 calculates the unit operation model located in the region (A + 1) or region (A - 1) (step S22).

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

[0304] Therefore, the in-furnace reaction calculation device 20 models the reaction process occurring in the rotary kiln 1 by a combination of a plurality of unit operation models, and performs calculations based on the values set for each of the unit operation models in accordance with the connection order of the unit operation models. In the present embodiment, after calculating each unit operation model in order from the charging end 14A side to the discharging end 14B side of the rotary kiln 1, it is performed from the discharging end 14B side to the charging end 14A side (see FIG. 6). Then, a series of operations are repeated until the calculated value in the predetermined region A satisfies the predetermined convergence condition. As a result, the calculation results such as the flow rates of the combustion gas and the raw ore in each region within the rotary kiln 1 are derived. Thereby, it becomes possible to more accurately analyze the behavior of each substance constituting the combustion gas and the raw ore throughout the kiln body 11 of the rotary kiln 1.

[0305] In this way, the in-furnace reaction calculation device 20 can accurately calculate the reaction in the rotary kiln 1 while taking into account at least that the gas region in the reactor shows two different reaction behaviors in the overall reaction process in the rotary kiln 1, and can reduce the burden involved in the calculation. Therefore, the in-furnace reaction calculation device 20 can more accurately analyze the behavior of each substance in the rotary kiln 1 with a low load while changing the operating conditions of the rotary kiln 1 (for example, the size and rotation speed of the rotary kiln 1, the type and supply amount of the raw ore). Thus, the in-furnace reaction calculation device 20 can be effectively utilized for changing the type and supply amount of the raw ore, pre-examining the improvement of the facilities of the rotary kiln 1, investigating the influence by operating conditions, etc.

[0306] Note that in this embodiment, the raw material supplied into the rotary kiln 1 may be a raw material other than the raw ore.

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

[0308] In this embodiment, in addition to the rotary kiln 1 shown in FIG. 1, any reactor that heats while moving the raw ore from the charging end 14A side toward the discharging end 14B side may be used.

[0309] As described above, the embodiments have been explained. However, the above embodiments are presented as examples, and the present invention is not limited by the above embodiments. The above embodiments can be implemented in various other forms, and various combinations, omissions, replacements, changes, etc. can be made without departing from the gist of the invention. These embodiments and their modifications are included in the scope and gist of the invention, and are included in the invention described in the claims and the equivalent scope thereof.

Explanation of Reference Numerals

[0310] 1 Rotary kiln 11 Kiln body 16 Burner 20. Calculation device for reactor reactions 201 Distribution section of combustion materials (distribution section) 202 1st gas phase quantity correction section 203 Second gas phase quantity correction section 204 Gas Mixture Calculation Unit 205 First Bed Layer Reserve Quantity Correction Section 206 First Bed Layer Mass Correction Section 207 Second Bed Layer Reserve Quantity Correction Section 208 Second Bed Layer Mass Correction Section 209 Inlet gas distribution section 210-1 First gas reaction amount calculation section 210-2 Second gas reaction amount calculation section 211-1 First substance reaction amount calculation section 211-2 Second substance reaction amount calculation section 212-1 First bed layer reaction volume calculation section 212-2 Second bed layer reaction volume calculation section 213-1 First inflow gas reaction amount calculation section 213-2 Second inflow gas reaction amount calculation section 214-1 First gas region equilibrium reaction calculation section 214-2 Second gas region equilibrium reaction calculation section 215-1 First bed layer equilibrium reaction calculation section 215-2 Second bed layer equilibrium reaction calculation section 216-1 First gas mixture calculation section 216-2 Second mixture calculation section 217-1 1st Mixed Material Calculation Department 217-2 2nd Mixed Material Calculation Department 218-1 First mixed bed layer calculation section 218-2 Second mixed bed layer calculation section 219-1 First mixed inflow gas calculation section 219-2 Second mixed inflow gas calculation section 220 First Mixed Material Quantity Correction Section 221 Inflow gas combining section 222 Gas volume summation unit 223 Second Mixture Primary Quantity Correction Department 224 Second mixture material secondary quantity correction unit G1 Gas phase G11 First gas phase G12 First gas phase containing additives G121 First diffusion gas G122 Second diffusion gas G13 Modified first gas phase G14 First mixed gas phase G21 Second gas phase G22A Modified second gas phase G22B Adjusted second gas phase G23 Second mixed gas phase G24 Second total gas phase G221 Inflow gas G221-1 First inflow gas G221-2 Second inflow gas G222A First mixed inflow gas G222B Second mixed inflow gas G223 Mixed inflow gas S1 Bed layer S11 First bed layer S12 First bed layer containing additives S13A Preliminary modified first bed layer S13B Modified first bed layer S14 First mixed bed layer S21 Second bed layer S21A Modified second bed layer S21B Modified second bed layer S22 Second mixed bed layer M1 Substance M11 First substance M12 First modified mixture substance M21 Second substance M3 Third substance M4 Fourth substance M5 Fifth substance M6 Sixth substance M7 Seventh substance M8 Eighth substance M9 Ninth substance GR1 First gas reaction component Gr1 First unreacted gas component GP1 First gas generated component GR2 Second gas reaction component Gr2 Second unreacted gas component GP2 Second gas generated component MR11 First substance reaction component Mr11 First substance unreacted component MP11 First substance generated component MR21 Second substance reaction component Mr21 Second substance unreacted component MP21 Second substance generated component GR221 First inflow gas reaction component Gr221 First inflow gas unreacted component GP221 First inflow gas generated component GR222 Second inflow gas reaction component Gr222 Second inflow gas unreacted component GP222 Second inflow gas generated component SR1 First bed layer reaction component Sr1 First bed layer unreacted component SP1 First bed layer generated component SR2 Second bed layer reaction component Sr2 Second bed layer unreacted component SP2 Second bed layer generated component

Claims

1. A calculation device for in-furnace reactions that dries and reduces raw ore by moving the raw ore supplied from one end side of a reactor toward the other end side and bringing the raw ore into contact with combustion gas supplied from the other end side, comprising: a first gas-phase material quantity correction unit that corrects the material quantity of a first gas phase flowing through a first gas region, which is one of two gas regions distributed according to a gas region influence factor including product gas generated from a bed layer containing the raw ore, among the gas regions containing the combustion gas, to obtain a corrected first gas phase; a second gas-phase material quantity correction unit that corrects the material quantity of a second gas phase flowing through a second gas region, which is the other of the two gas regions, to obtain a corrected second gas phase; a gas mixing quantity calculation unit that calculates the flow rate of the inflow gas moving from the second gas region to the bed layer; a first gas region equilibrium reaction calculation unit that calculates an equilibrium reaction between the corrected first gas phase and a first substance present in the first gas region, among substances including at least one of a solid substance and a liquid substance present in the gas region; a second gas region equilibrium reaction calculation unit that calculates an equilibrium reaction between the corrected second gas phase and a second substance present in the second gas region, among the substances; a bed layer equilibrium reaction calculation unit that calculates an equilibrium reaction between the bed layer and the inflow gas; A calculation device for in-furnace reactions comprising the above components.

2. a first gas reaction quantity calculation unit that calculates a first gas reaction quantity contributing to the equilibrium state between the first substance in the first gas phase and obtains a first gas reaction component contributing to the equilibrium reaction of the first gas phase; a first substance reaction quantity calculation unit that calculates a first substance reaction quantity contributing to the equilibrium state between the first substance and the first gas phase and obtains a first substance reaction component contributing to the equilibrium reaction of the first substance; The calculation device for in-furnace reactions according to Claim 1, comprising the above components.

3. a first mixed gas phase calculation unit that calculates at least the flow rate of a first mixed gas phase obtained by mixing a first unreacted gas component other than the first gas reaction component generated by the first gas reaction quantity calculation unit, and a first gas generation component generated by the reaction between the first gas reaction component and the first substance reaction component generated by the first substance reaction quantity calculation unit in the first gas region equilibrium reaction calculation unit; a first mixture substance calculation unit that calculates at least the flow rate of a first mixture substance obtained by mixing a first unreacted substance component and a first substance generation component generated by the reaction between the first gas reaction component and the first substance reaction component in the first gas region equilibrium reaction calculation unit; The calculation device for in-furnace reactions according to Claim 2, comprising the above components.

4. A second gas reaction amount calculation unit that calculates a second gas reaction amount contributing to the equilibrium state of the second substance in the second gas phase and obtains a second gas reaction component contributing to the equilibrium reaction of the second gas phase; A second substance reaction amount calculation unit that calculates a second substance reaction amount contributing to the equilibrium state of the second substance with the second gas phase and obtains a second substance reaction component contributing to the equilibrium reaction of the second substance; The in-furnace reaction calculation device according to any one of claims 1 to 3, comprising:

5. A second mixed gas phase calculation unit that calculates at least the flow rate of a second mixed gas phase obtained by mixing a second unreacted gas component other than the second gas reaction component generated by the second gas reaction amount calculation unit and a second gas generated component generated by the reaction of the second unreacted gas component and the second gas reaction component generated by the second gas region equilibrium reaction calculation unit and the second substance reaction component generated by the second substance reaction amount calculation unit; A second mixture substance calculation unit that calculates at least the flow rate of a second mixture substance obtained by mixing a second unreacted substance component and a second substance generated component generated by the reaction of the second gas reaction component and the second substance reaction component in the second gas region equilibrium reaction calculation unit; The in-furnace reaction calculation device according to claim 4, comprising:

6. A first mixture substance amount correction unit that corrects the amount of the first mixture substance in consideration of the amount of transfer of a fourth substance containing at least one of a solid substance or a liquid substance included in the first mixture substance containing the first substance generated component generated after the equilibrium reaction of the first substance and the unreacted component of the first substance in the first gas region equilibrium reaction calculation unit to the second gas region, and the amount of transfer of a third substance containing at least one of a solid substance or a liquid substance included in the second mixture substance containing the second substance generated component generated after the equilibrium reaction of the second substance and the unreacted component of the second substance in the second gas region equilibrium reaction calculation unit to the first mixture substance. The in-furnace reaction calculation device according to any one of claims 1 to 5.

7. Comprising an inlet gas distribution unit that distributes the inlet gas into a first inlet gas and a second inlet gas; The bed layer equilibrium reaction calculation unit is A first bed layer equilibrium reaction calculation unit that calculates an equilibrium reaction between a first bed layer, which is one of two types of distributed bed layers distributed according to the segregation and concentration of solid components and liquid components contained in the bed layer, and the first inlet gas; A second bed layer equilibrium reaction calculation unit that calculates an equilibrium reaction between a second bed layer, which is the other of the two types of distributed bed layers, and the second inlet gas; The in-furnace reaction calculation device according to any one of claims 1 to 6 having

8. A first bed layer reaction amount calculation unit that calculates a first bed layer reaction amount contributing to an equilibrium state of the first bed layer with the first inflowing gas and obtains a first bed layer reaction component contributing to the equilibrium reaction of the first bed layer; A first inflowing gas reaction amount calculation unit that calculates a gas reaction amount contributing to an equilibrium state of the first inflowing gas with the first bed layer and obtains a first inflowing gas reaction component contributing to the equilibrium reaction of the first inflowing gas; The in-furnace reaction calculation device according to claim 7, comprising

9. 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 unreacted component of the first bed layer other than the first bed layer reaction component generated by the first bed layer reaction amount calculation unit and a first bed layer generated component generated by the reaction of the first bed layer reaction component and the first inflowing gas reaction component in the first bed layer equilibrium reaction calculation unit; A first mixed inflowing gas calculation unit that calculates at least a flow rate of a first mixed inflowing gas obtained by mixing a first unreacted component of the first inflowing gas other than the first inflowing gas reaction component generated by the first inflowing gas reaction amount calculation unit and a first inflowing gas generated component generated by the reaction of the first bed layer reaction component and the first inflowing gas reaction component in the first bed layer equilibrium reaction calculation unit; The in-furnace reaction calculation device according to claim 8, comprising

10. A second bed layer reaction amount calculation unit that calculates a second bed layer reaction amount contributing to an equilibrium state of the second bed layer with the second inflowing gas and obtains a second bed layer reaction component contributing to the equilibrium reaction of the second bed layer; A second inflowing gas reaction amount calculation unit that calculates a gas reaction amount contributing to an equilibrium state of the second inflowing gas with the second bed layer and obtains a second inflowing gas reaction component contributing to the equilibrium reaction of the second inflowing gas; The in-furnace reaction calculation device according to any one of claims 7 to 9, comprising

11. 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 unreacted component of the second bed layer other than the second bed layer reaction component generated by the second bed layer reaction amount calculation unit and a second bed layer generated component generated by the reaction of the second bed layer reaction component and the second inflowing gas reaction component in the second bed layer equilibrium reaction calculation unit; A second mixed inflow gas calculation unit that calculates at least the flow rate of a second mixed inflow gas obtained by mixing an unreacted portion of the second inflow gas other than the reacted portion of the second inflow gas generated by the second inflow gas reaction amount calculation unit, a second inflow gas generated component generated by the reaction of the reacted portion of the second bed layer and the reacted portion of the second inflow gas in the second bed layer equilibrium reaction calculation unit, The in-furnace reaction calculation device according to claim 10, comprising

12. Considering the amount of movement of a seventh substance containing at least one of a solid substance and a liquid substance in the first bed layer to the second gas region, a second substance generated component generated after the equilibrium reaction of the second substance in the second gas region equilibrium reaction calculation unit, and an unreacted portion of the second substance, the in-furnace reaction calculation device according to any one of claims 7 to 11, comprising a first bed layer preliminary material amount correction unit that corrects the material amount of the first bed layer considering the amount of movement of an eighth substance containing at least one of a solid substance and a liquid substance in a second mixture substance containing the second substance generated component and the unreacted portion of the second substance to the first bed layer.

13. The in-furnace reaction calculation device according to any one of claims 7 to 12, comprising a second bed layer preliminary material amount correction unit that corrects the material amount of the second bed layer considering the amount of movement of a ninth substance containing at least one of a solid substance and a liquid substance in the second bed layer to the second gas region.

14. An inflow gas summation unit that mixes a first inflow gas generated component generated by the reaction of a first bed layer reaction portion contributing to the equilibrium reaction of the first bed layer and the first inflow gas in the first bed layer equilibrium reaction calculation unit, and a second inflow gas generated component generated by the reaction of a second bed layer reaction portion contributing to the equilibrium reaction of the second bed layer and the second inflow gas in the second bed layer equilibrium reaction calculation unit. The in-furnace reaction calculation device according to any one of claims 7 to 13.

15. The in-furnace reaction calculation device according to claim 14, comprising a gas amount summation unit that mixes a mixed inflow gas generated by mixing the first inflow gas generated component and the second inflow gas generated component in the inflow gas summation unit, and a second gas generated component generated by the reaction of the second gas phase and the second substance in the second gas region equilibrium reaction calculation unit.

16. In the second gas region equilibrium reaction calculation unit, considering the amount of movement of an eighth substance containing at least one of a solid substance and a liquid substance in a second mixed substance containing a second substance generation component generated after the equilibrium reaction of the second substance and an unreacted portion of the second substance to the bed layer, the amount of movement of a seventh substance containing at least one of a solid substance and a liquid substance in the bed layer to the second mixed substance, and the amount of movement of a ninth substance containing at least one of a solid substance and a liquid substance in the second bed layer to the second gas region, the amount of substance of the second mixed substance is corrected to calculate a second corrected mixed substance, and a furnace internal reaction calculation device according to any one of claims 7 to 15, comprising a first-order amount-of-substance correction unit for the second mixed substance.

17. In the second-order amount-of-substance correction unit for the second mixed substance, considering the amount of movement of a third substance containing at least one of a solid substance and a liquid substance in the second corrected mixed substance generated by the first-order amount-of-substance correction unit for the second mixed substance to the first gas region, and the amount of movement of a fourth substance containing at least one of a solid substance and a liquid substance in a first mixed substance containing a first substance generation component generated after the equilibrium reaction of the first substance and an unreacted portion of the first substance to the second corrected mixed substance in the first gas region equilibrium reaction calculation unit, the amount of substance of the second corrected mixed substance is corrected, and a furnace internal reaction calculation device according to claim 16, comprising a second-order amount-of-substance correction unit for the second mixed substance.

18. A furnace internal reaction calculation device according to any one of claims 1 to 17, comprising a gas region distribution unit that distributes the gas region to two types of the gas regions in advance according to a gas region influence factor.

19. A furnace internal reaction calculation device according to any one of claims 1 to 18, comprising a first substance distribution unit that distributes the first substance to two types of the first substances in advance according to a gas region influence factor contained in the first substance.

20. Comprising a distribution unit that distributes a combustion material containing at least one of an added gas phase and an added bed layer, which is introduced during the movement of the raw ore, to the added gas phase and the added bed layer in terms of mass flow rate. The first gas region equilibrium reaction calculation unit calculates an equilibrium reaction between a first mixed gas phase containing the first gas phase and the added gas phase and the first substance. The second gas region equilibrium reaction calculation unit calculates an equilibrium reaction between a second mixed gas phase containing the first gas phase and the added gas phase and the second substance. The bed layer equilibrium reaction calculation unit is the in-furnace reaction calculation device according to any one of claims 1 to 19, which calculates the equilibrium reaction between the first mixed bed layer including the bed layer and the added bed layer and the inflowing gas.

21. When it is assumed that the inside of the reactor is divided into a plurality of regions along its major axis direction, The method is repeatedly performed along the flow of the raw ore or the combustion gas in the plurality of regions, and is repeatedly performed until the difference between the calculated value in a predetermined region and the previous calculated value in that region falls within a predetermined range. The in-furnace reaction calculation device according to any one of claims 1 to 20.

22. A method for calculating an in-furnace reaction in which raw ore supplied from one end side of a reactor is moved toward the other end side, and the raw ore is brought into contact with combustion gas supplied from the other end side to be dried and reduced, A first gas phase material amount correction step of correcting the material amount of a first gas phase flowing through a first gas region, which is one of two gas regions distributed according to a gas region influence factor including a product gas generated from a bed layer containing the raw ore, among the gas regions containing the combustion gas, to obtain a corrected first gas phase; A second gas phase material amount correction step of correcting the material amount of a second gas phase flowing through a second gas region, which is the other of the two gas regions, to obtain a corrected second gas phase; A gas mixing amount calculation step of calculating the flow rate of the inflowing gas moving from the second gas region to the bed layer; A first gas region equilibrium reaction calculation step of calculating an equilibrium reaction between the corrected first gas phase and a first substance present in the first gas region among substances including at least one of a solid substance and a liquid substance present in the gas region; A second gas region equilibrium reaction calculation step of calculating an equilibrium reaction between the corrected second gas phase and a second substance present in the second gas region among the substances; A bed layer equilibrium reaction calculation step of calculating an equilibrium reaction between the bed layer and the inflowing gas; A method for calculating an in-furnace reaction including the above steps.

Images