DECARBURIZATION END POINT DETERMINATION METHOD, DECARBURIZATION END POINT DETERMINATION DEVICE, SECONDARY REFINING OPERATION METHOD FOR STEEL MAKING, AND METHOD FOR PRODUCING MOLTEN STEEL

Abstract

A decarburization end point determination method includes: estimating the carbon concentration and oxygen concentration of the molten steel and carbon dioxide gas concentration of exhaust gas in the vacuum chamber by using measurement values of the carbon concentration and the oxygen concentration of the molten steel, a measurement value of internal pressure of the vacuum chamber, and a model formula; correcting a parameter included in the model formula to reduce at least one of a difference between an estimate value and a measurement value of the oxygen concentration and a difference between an estimate value and a measurement value of the carbon dioxide gas concentration of the exhaust gas; estimating the carbon concentration of the molten steel by using the model formula in which the parameter is corrected; and determining timing when an estimate value reaches a target value as the completion time point of the vacuum decarburization treatment.

Claims

1-8. (canceled)

9. A decarburization end point determination method for determining a completion time point of vacuum decarburization treatment in a facility that performs the vacuum decarburization treatment of reducing carbon concentration of molten steel by degassing a vacuum chamber, the method comprising: estimating the carbon concentration and oxygen concentration of the molten steel and carbon dioxide gas concentration of exhaust gas in the vacuum chamber during the vacuum decarburization treatment by using measurement values of the carbon concentration and the oxygen concentration of the molten steel before the vacuum decarburization treatment is started, a measurement value of internal pressure of the vacuum chamber, and a model formula of the vacuum decarburization treatment; correcting a parameter included in the model formula in such a manner as to reduce at least one of a difference between an estimate value and a measurement value of the oxygen concentration of the molten steel at timing when the oxygen concentration of the molten steel is measured during the vacuum decarburization treatment and a difference between an estimate value and a measurement value of the carbon dioxide gas concentration of the exhaust gas at timing when the carbon dioxide gas concentration of the exhaust gas is measured during the vacuum decarburization treatment; estimating the carbon concentration of the molten steel by using the model formula in which the parameter is corrected; and determining timing when an estimate value reaches a target value as the completion time point of the vacuum decarburization treatment.

10. The decarburization end point determination method according to claim 9, wherein the correcting includes a step of calculating a probable value of the parameter when at least one of the measurement value of the oxygen concentration of the molten steel and the measurement value of the carbon dioxide gas concentration of the exhaust gas is acquired, by an inverse analysis using Bayesian inference.

11. The decarburization end point determination method according to claim 9, wherein the parameter is a decarburization capacity coefficient.

12. The decarburization end point determination method according to claim 10, wherein the parameter is a decarburization capacity coefficient.

13. A decarburization end point determination device that determines a completion time point of vacuum decarburization treatment in a facility that performs the vacuum decarburization treatment of reducing carbon concentration of molten steel by degassing a vacuum chamber, the device comprising a processor comprising hardware, the processor being configured to: estimate the carbon concentration and oxygen concentration of the molten steel and carbon dioxide gas concentration of exhaust gas in the vacuum chamber during the vacuum decarburization treatment by using measurement values of the carbon concentration and the oxygen concentration of the molten steel before the vacuum decarburization treatment is started, a measurement value of internal pressure of the vacuum chamber, and a model formula of the vacuum decarburization treatment; correct a parameter included in the model formula in such a manner as to reduce at least one of a difference between an estimate value and a measurement value of the oxygen concentration of the molten steel at timing when the oxygen concentration of the molten steel is measured during the vacuum decarburization treatment and a difference between an estimate value and a measurement value of the carbon dioxide gas concentration of the exhaust gas at timing when the carbon dioxide gas concentration of the exhaust gas is measured during the vacuum decarburization treatment; estimate the carbon concentration of the molten steel by using the model formula in which the parameter is corrected; and determine timing when an estimate value reaches a target value as the completion time point of the vacuum decarburization treatment.

14. The decarburization end point determination device according to claim 13, wherein the processor is configured to calculate a probable value of the parameter when at least one of the measurement value of the oxygen concentration of the molten steel and the measurement value of the carbon dioxide gas concentration of the exhaust gas is acquired, by an inverse analysis using Bayesian inference.

15. The decarburization end point determination device according to claim 13, wherein the parameter is a decarburization capacity coefficient.

16. The decarburization end point determination device according to claim 14, wherein the parameter is a decarburization capacity coefficient.

17. A secondary refining operation method for steel making, comprising performing operation end determination of a secondary refining process for steel making, by the decarburization end point determination method according to claim 9.

18. A secondary refining operation method for steel making, comprising performing operation end determination of a secondary refining process for steel making, by the decarburization end point determination method according to claim 10.

19. A secondary refining operation method for steel making, comprising performing operation end determination of a secondary refining process for steel making, by the decarburization end point determination method according to claim 11.

20. A secondary refining operation method for steel making, comprising performing operation end determination of a secondary refining process for steel making, by the decarburization end point determination method according to claim 12.

21. A method for producing molten steel, comprising producing the molten steel by using the secondary refining operation method for steel making according to claim 17.

22. A method for producing molten steel, comprising producing the molten steel by using the secondary refining operation method for steel making according to claim 18.

23. A method for producing molten steel, comprising producing the molten steel by using the secondary refining operation method for steel making according to claim 19.

24. A method for producing molten steel, comprising producing the molten steel by using the secondary refining operation method for steel making according to claim 20.

Description

BRIEF DESCRIPTION OF DRAWINGS

[0021] FIG. 1 is a view illustrating a configuration of a vacuum degassing facility to which a decarburization end point determination method that is one embodiment of the present invention is applied.

[0022] FIG. 2 is a flowchart illustrating a flow of decarburization end point determination processing that is one embodiment of the present invention.

[0023] FIG. 3 is a view illustrating a causal relationship between variables related to vacuum decarburization treatment.

[0024] FIG. 4 is a view illustrating time-series changes of input variables.

[0025] FIG. 5 is a view illustrating time-series changes of output variables according to a conventional method.

[0026] FIG. 6 is a view illustrating time-series changes of output variables according to the present invention.

DESCRIPTION OF EMBODIMENTS

[0027] Hereinafter, a decarburization end point determination method that is one embodiment of the present invention will be described in detail with reference to the drawings.

[0028] [Configuration of a Vacuum Degassing Facility]

[0029] First, with reference to FIG. 1, a configuration of a vacuum degassing facility to which the decarburization end point determination method that is the one embodiment of the present invention is applied will be described.

[0030] FIG. 1 is a view illustrating the configuration of the vacuum degassing facility to which the decarburization end point determination method that is the one embodiment of the present invention is applied.

[0031] A vacuum degassing facility 1 to which the decarburization end point determination method according to the one embodiment of the present invention illustrated in FIG. 1 is applied is one of secondary refining facilities in a steel mill of a steelworks. As illustrated in FIG. 1, the vacuum degassing facility 1 includes a ladle 2 that houses molten steel S, a suction pipe 3 and an exhaust pipe 4 immersed in the molten steel S, a vacuum chamber 5, a vacuum device 6, an oxygen feeding lance 7, and a hopper 8.

[0032] When vacuum decarburization treatment is performed by utilization of the vacuum degassing facility 1, an argon gas (reflux gas) A1 is blown into the suction pipe 3 while pressure in the vacuum chamber 5 is reduced to several tens to one hundred mhPa by utilization of the vacuum device 6. Thus, the molten steel S in the suction pipe 3 is raised into the vacuum chamber 5 by an air lift pump action, and the molten steel S is refluxed into the ladle 2 via the exhaust pipe 4. In this process of pressure reduction and reflux, oxygen and carbon in the molten steel S react (decarburization reaction) and form a carbon dioxide gas, and the carbon dioxide gas is taken into bubbles of the argon gas A1 in the molten steel S. Then, the bubbles of the argon gas A1 that takes in the carbon dioxide gas are exhausted to the outside of the vacuum chamber 5 as exhaust gas A2 via the vacuum chamber 5. This reduces carbon concentration of the molten steel S from about 300 ppm to 10 to 100 ppm. In order to accelerate the decarburization reaction, a pure oxygen gas A3 may be blown from the oxygen feeding lance 7 inserted into the vacuum chamber 5 from an upper part of the vacuum chamber 5 at an initial stage of the vacuum decarburization treatment, and the oxygen may be dissolved in the molten steel S. In addition, in the latter half of the vacuum decarburization treatment, an A1 ore or the like may be fed from the hopper 8 in order to stop the decarburization reaction.

[0033] The vacuum degassing facility 1 includes a measurement value collection device 10, a decarburization end point estimation device 11, and a display device 12 as a control system.

[0034] The measurement value collection device 10 is connected to an oxygen feeding lance flowmeter 21, a reflux gas flowmeter 22, a vacuum chamber internal pressure meter 23, a thermometer 24, a dissolved oxygen meter 25, an auxiliary raw material meter 26, an exhaust gas CO concentration meter 27, an exhaust gas CO.sub.2 concentration meter 28, an exhaust gas 02 concentration meter 29, and an exhaust gas flowmeter 30. The measurement value collection device 10 acquires measurement values from these measurement devices at every predetermined control cycle, and outputs the acquired measurement values to the decarburization end point estimation device 11.

[0035] The oxygen feeding lance flowmeter 21 measures a flow rate of the pure oxygen gas A3 blown into the vacuum chamber 5 from the oxygen feeding lance 7 (oxygen flow rate), and inputs a measurement value of the oxygen flow rate to the measurement value collection device 10. The reflux gas flowmeter 22 measures a flow rate of the argon gas A1 blown into the suction pipe 3, and inputs a measurement value of the flow rate of the argon gas A1 to the measurement value collection device 10. The vacuum chamber internal pressure meter 23 measures internal pressure of the vacuum chamber 5 (vacuum chamber internal pressure) and inputs a measurement value of the vacuum chamber internal pressure to the measurement value collection device 10. The thermometer 24 measures an inner wall temperature of a lower portion of the vacuum chamber 5, and inputs a measurement value of the inner wall temperature to the measurement value collection device 10. The dissolved oxygen meter 25 measures dissolved oxygen concentration of the molten steel S in the ladle 2, and inputs a measurement value of the dissolved oxygen concentration to the measurement value collection device 10.

[0036] The auxiliary raw material meter 26 measures weight of auxiliary raw material fed from the hopper 8, and inputs a measurement value of the weight of the auxiliary raw material to the measurement value collection device 10. The exhaust gas CO concentration meter 27 measures CO concentration (carbon dioxide gas concentration) of the exhaust gas A2 exhausted from the vacuum device 6, and inputs a measurement value of the CO concentration to the measurement value collection device 10. The exhaust gas CO.sub.2 concentration meter 28 measures CO.sub.2 concentration of the exhaust gas A2 exhausted from the vacuum device 6, and inputs a measurement value of the CO.sub.2 concentration to the measurement value collection device 10. The exhaust gas 02 concentration meter 29 measures oxygen concentration of the exhaust gas A2 exhausted from the vacuum device 6, and inputs a measurement value of the oxygen concentration to the measurement value collection device 10. The exhaust gas flowmeter 30 measures a flow rate of the exhaust gas A2 exhausted from the vacuum device 6, and inputs a measurement value of the flow rate to the measurement value collection device 10.

[0037] The decarburization end point estimation device 11 includes an information processing device such as a computer, and includes a random access memory (RAM) 11a, a read only memory (ROM) lib, and an arithmetic processing unit 11c. The RAM 11a, the ROM lib, and the arithmetic processing unit 11c are electrically connected via a bus wiring line lid.

[0038] The RAM 11a functions as a working space of the arithmetic processing unit 11c by temporarily storing a computer program executed by the arithmetic processing unit 11c and various kinds of data necessary for executing the computer program. In the present embodiment, when decarburization end point determination processing (described later) is executed, values of variables and constants illustrated in Table 1 to 7 in the following are stored in the RAM 11a in addition to a decarburization end point estimation program 11b1 (described later), and the values of the variables are sequentially updated in a process of the decarburization end point determination processing.

[0039] State variables included in a model formula of a decarburization reaction (described later) are illustrated in Table 1. An adjustment parameter of the model formula of the decarburization reaction is illustrated in Table 2. Input variables of the model formula of the decarburization reaction are illustrated in Table 3. Output variables calculated from the model formula of the decarburization reaction are illustrated in Table 4. Constants included in the model formula of the decarburization reaction are illustrated in Table 5. Intermediate variables used in the model formula of the decarburization reaction (variables temporarily used to calculate the state variables illustrated in Table 1) are illustrated in Table 6. Measurement values of the output variables calculated from the model formula of the decarburization reaction are illustrated in Table 7.

TABLE-US-00001 TABLE 1 Term Unit Description m.sub.c,v kg Carbon weight in molten steel in vacuum chamber m.sub.c,l kg Carbon weight in molten steel in ladle m.sub.o,v kg Dissolved oxygen weight in molten steel in vacuum chamber m.sub.o,l kg Dissolved oxygen weight in molten steel in ladle C.sub.co,g mol mol.sup.−1 Molar fraction of carbon dioxide gas in gas phase in vacuum chamber

TABLE-US-00002 TABLE 2 Term Unit Description K.sub.v kg kg.sup.−1 Decarburization Pa.sup.−1 m.sup.−1 s.sup.−1 capacity coefficient

TABLE-US-00003 TABLE 3 Term Unit Description p.sub.v pa Vacuum chamber internal pressure V.sub.arg m.sup.3 s.sup.−1 Argon gas flow rate V.sub.off25 kg EXHAUST GAS FLOW RATE f.sub.c,add kg s.sup.−1 Mass flow rate of carbon from fed auxiliary raw material to molten steel which rate is calculated from weight measurement value of fed auxiliary raw material f.sub.a,add kg s.sup.−1 Mass flow rate of aluminum from fed auxiliary raw material to molten steel which rate is calculated from weight measurement value of fed auxiliary raw material

TABLE-US-00004 TABLE 4 Term Unit Description C.sub.c,l kg kg.sup.−1 Carbon mass fraction of molten steel in ladle C.sub.o,l kg kg.sup.−1 Oxygen mass fraction of molten steel in ladle C.sub.co,g mol mol.sup.−1 Molar fraction of carbon dioxide gas in gas phase in vacuum chamber

TABLE-US-00005 TABLE 5 Term Unit Value Description D.sub.s m 6.0e−1 Immersion pipe diameter B m.sup.−4/3 kg s.sup.−4/3 7.0e+3 Mass velocity coefficient of reflux g N kg.sup.−1 9.8 Gravitational acceleration M.sub.Al kg mol.sup.−1 2.7e−2 Molar mass of aluminum M.sub.Ar kg mol.sup.−1 4.0e−2 Molar mass of argon M.sub.C kg mol.sup.−1 1.2e−2 Molar mass of carbon M.sub.O kg mol.sup.−1 1.6e−2 Molar mass of oxygen P.sub.ail Pa 1.0e+5 Pressure of reflux argon gas P.sub.atm Pa 1.0e+5 Atmospheric pressure P.sub.i Pa Pressure of formed bubble in molten steel R JK.sup.−1 mol.sup.−1 8.3 Gas constant T.sub.s K Molten steel temperature T.sub.g K Gas temperature in vacuum chamber p kg m.sup.−3 7.0e+3 Molten steel density V.sub.v m.sup.3 1.5 Vacuum chamber volume

TABLE-US-00006 TABLE 6 Term Unit Description Calculation formula C.sub.c,v kg kg.sup.−1 Carbon mass m.sub.c,v/m.sub.s,v fraction of molten steel in vacuum chamber C.sub.c,l kg kg.sup.−1 Carbon mass m.sub.c,l/m.sub.s,l fraction of molten steel in ladle C.sub.o,v kg kg.sup.−1 Oxygen mass m.sub.o,v/m.sub.s,v fraction of molten steel in vacuum chamber C.sub.o,l kg kg.sup.−1 Oxygen mass m.sub.o,l/m.sub.s,l fraction of molten steel in ladle f.sub.arg kg s.sup.−1 Mass velocity of v.sub.argM.sub.ArP.sub.atm/(RT) reflux argon gas f.sub.dec kg s.sup.−1 Decarburization −K.sub.v (K.sub.Cc, vCo, v−p.sub.v) rate f.sub.cir kg s.sup.−1 Reflux flow rate B .Math. D.sub.s.sup.4/3v.sub.arg.sup.1/3 (log (P.sub.ail/p.sub.v) ) .sup.1/3 K Pa kg.sup.2 Carbon-oxygen 10 .sup.(10447/T.sup.s .sup.+ .sup.18.038) kg.sup.−2 reaction coefficient

TABLE-US-00007 TABLE 7 Term Unit Description ĉ.sub.c,l kg kg.sup.−1 Measurement value of carbon mass fraction of molten steel in ladle ĉ.sub.o,l kg kg.sup.−1 Measurement value of oxygen mass fraction of molten steel in ladle ĉ.sub.co,g mol mol.sup.−1 Measurement value of molar fraction of carbon dioxide gas in gas phase in vacuum chamber

[0040] The ROM 11b includes a non-volatile storage device, and stores a computer program such as the decarburization end point estimation program 11b1, and various kinds of control data.

[0041] The arithmetic processing unit 11c includes an electronic circuit such as a central processing unit (CPU) inside the information processing device. The arithmetic processing unit 11c controls an operation of the entire decarburization end point estimation device 11 by loading the computer program stored in the ROM 11b into the RAM 11a and executing the loaded computer program. In the present embodiment, the arithmetic processing unit 11c functions as an initialization unit 11c1, a measurement value reading unit 11c2, a state estimation unit 11c3, and an end point determination unit 11c4 by executing the decarburization end point estimation program 11b1. Functions of these units will be described later.

[0042] The display device 12 includes a known display device such as a liquid crystal display device, and visibly displays various kinds of information in accordance with a control signal from the decarburization end point estimation device 11. In the present embodiment, the display device 12 visually displays an input variable value of the model formula of the decarburization reaction, an output variable value, and a measurement value of the output variable at each calculation cycle similar to that of the decarburization end point estimation device 11. Furthermore, the display device 12 visually displays information related to an end point of the vacuum decarburization treatment.

[0043] In the vacuum degassing facility 1 having such a configuration, the decarburization end point estimation device 11 estimates the carbon concentration of the molten steel S and determines the end point of the vacuum decarburization treatment by executing the following decarburization end point determination processing. In the following, the operation of the decarburization end point estimation device 11 in execution of the decarburization end point determination processing will be described with reference to a flowchart illustrated in FIG. 2.

[0044] [Decarburization End Point Determination Processing]

[0045] FIG. 2 is a flowchart illustrating a flow of the decarburization end point determination processing that is one embodiment of the present invention. The decarburization end point determination processing illustrated in FIG. 2 is started at timing when the vacuum device 6 starts suction processing in the vacuum chamber 5, and the decarburization end point determination processing proceeds to processing of Step S1. Note that the initialization unit 11c1 resets a value of a program counter i that counts a loop count (time step) of the decarburization end point determination processing to 0 at the timing when the decarburization end point determination processing is started.

[0046] In the processing of Step S1, the initialization unit 11c1 stores the values of the constants listed in Table 5 into a region of the RAM 11a which region corresponds to the constants that are included in the model formula of the decarburization reaction and that are listed in Table 5, and stores an initial value 0.0080 [kg kg.sup.−1 Pa.sup.−1 m.sup.−1s.sup.−1] of a decarburization capacity coefficient K.sub.v into a region of the RAM 11a which region corresponds to the decarburization capacity coefficient K.sub.v listed in Table 2. Note that these values may be appropriately changed according to a facility form or an operation form of the vacuum degassing facility 1. As a result, the processing of Step S1 is completed, and the decarburization end point determination processing proceeds to processing of Step S2.

[0047] In the processing of Step S2, the initialization unit 11c1 stores various measurement values output from the measurement value collection device 10 into a region of the RAM 11a which region corresponds to the input variables of the model formula of the decarburization reaction, which variables are listed in Table 3, and the measurement values of the output variables which values are calculated from the model formula of the decarburization reaction and are listed in Table 7. In addition, the initialization unit 11c1 copies the values in the region of the RAM 11a which region corresponds to the measurement values of the output variables, which measurement values are listed in Table 7, to a region of the RAM 11a which region corresponds to the output variables (initial values of the output variables) listed in Table 4. Note that at the start of the decarburization end point determination processing (program counter i=0), values measured immediately after tapping of the molten steel from a converter in the previous process may be used as measurement values of a carbon mass fraction (carbon concentration) and an oxygen mass fraction (dissolved oxygen concentration) of the molten steel listed in Table 7. As a result, the processing of Step S2 is completed, and the decarburization end point determination processing proceeds to processing of Step S3.

[0048] In the processing of Step S3, the state estimation unit 11c3 assumes the measurement value of the carbon mass fraction of the molten steel, the measurement value of the oxygen mass fraction of the molten steel, and a measurement value of a molar fraction of the carbon dioxide gas (carbon dioxide gas concentration) of the exhaust gas which values are stored in the RAM 11a in the processing of Step S2. That is, the state estimation unit 11c3 assumes that an estimate value of the carbon mass fraction of the molten steel, an estimate value of the oxygen mass fraction of the molten steel, and an estimate value of the molar fraction of the carbon dioxide gas of the exhaust gas respectively match the measurement values. Then, by an inverse analysis using Bayesian inference (statistical causal inference), the state estimation unit 11c3 calculates change histories of the carbon mass fraction of the molten steel, the oxygen mass fraction of the molten steel, and the molar fraction of the carbon dioxide gas of the exhaust gas over a past predetermined period (such as past 60 seconds) which change histories are causes of such a result.

[0049] Here, in the Bayesian inference, the estimate values of the carbon mass fraction of the molten steel, the oxygen mass fraction of the molten steel, and the molar fraction of the carbon dioxide gas of the exhaust gas are expressed as probability distribution. Then, in Bayesian theory, when at least one of the carbon mass fraction of the molten steel, the oxygen mass fraction of the molten steel, and the molar fraction of the carbon dioxide gas of the exhaust gas becomes known, that is, when a measurement value thereof is acquired, a causal relationship between the variables related to the vacuum decarburization treatment illustrated in FIG. 3 is traced inversely. Thus, in the Bayesian theory, the most probable decarburization rate and decarburization capacity coefficient, and change histories of the carbon mass fraction of the molten steel, the oxygen mass fraction of the molten steel, and the molar fraction of the carbon dioxide gas of the exhaust gas are estimated.

[0050] Specifically, the state estimation unit 11c3 solves the following optimization problem. This means to calculate, on the assumption that the estimate values of the output variables at the current time step i match the measurement values, a value that is considered to be the most probable in model calculation with respect to each of the values of the state variables and the decarburization capacity coefficient within a past time horizon period h. Specifically, the decarburization reaction is expressed by model formulas expressed in the following formulas (6) to (11), and values of a state variable x and the decarburization capacity coefficient K.sub.v are given by probability distribution. The probability distribution expresses uncertainty that the values of the state variable x and the decarburization capacity coefficient K.sub.v are around these values. Thus, there is also uncertainty that a value of an output variable y during the progress of the decarburization reaction is around this. The following formula (5) expresses the above, and expresses the probability distribution of the output variable y (Gaussian distribution having a standard deviation σ with an estimate value y(i) being an average value) of when the state variable x and the decarburization capacity coefficient K.sub.v are in a given condition. However, a probability density function expressed in the formula (5) may be any kind. Note that an initial value of the probability distribution of the state variable x is set as the Gaussian distribution in which a value measured immediately before the start of the decarburization treatment is an average value and a value that is the same as a standard deviation of a measurement error quantified in advance is held as a standard deviation. However, even with a probability distribution setting method other than the present method, the function of the present invention is not impaired.

[0051] Then, when the value of the current output variable y is determined by actual measurement, the values of the past state variable x and decarburization capacity coefficient K.sub.v estimated uncertainly can also be estimated again with a certain degree of accuracy. The following formulas (1) to (4) express calculation formulas of going back to the past, and mean that the state variable x and the decarburization capacity coefficient K.sub.v return the most probable values (probability p being maximum) under the condition that the value of the output variable y is determined as a measurement value y-hat. Note that P( ) in the formula (1) means an inverse formula of the formula (5). As a result, the processing of Step S3 is completed, and the decarburization end point determination processing proceeds to processing of Step S4.

[00001]$\begin{array}{cc}\underset{K_v}{\max }P\left(x\left(i-h\right),.Math.,x\left(i-1\right),K_V.Math.y\left(i\right)=\hat{y}\left(i\right)\right)& \left(1\right)\end{array}$$\begin{array}{cc}x\left(i\right)=\begin{array}{c}{m}_{c,v}\left(i\right)\\ m_{c,l}\left(i\right)\\ m_{o,v}\left(i\right)\\ m_{o,l}\left(i\right)\\ .Math.c_{\mathrm{co},g}\left(i\right).Math.\end{array}& \left(2\right)\end{array}$$\begin{array}{cc}y\left(i\right)=\left[\begin{array}{c}{c}_{c,l}\left(i\right)\\ c_{o,l}\left(i\right)\\ c_{\mathrm{co},g}\left(i\right)\end{array}\right]& \left(3\right)\end{array}$$\begin{array}{cc}\hat{y}\left(i\right)=\left[\begin{array}{c}{\hat{c}}_{c,l}\left(i\right)\\ {\hat{c}}_{o,l}\left(i\right)\\ {\hat{c}}_{\mathrm{co},g}\left(i\right)\end{array}\right]& \left(4\right)\end{array}$$\begin{array}{cc}P\left(y\left(i\right).Math.x\left(i-h\right),.Math.,x\left(i-1\right),K_V\right)~N\left(y\left(i\right),\sigma \right)& \left(5\right)\end{array}$$\begin{array}{cc}x\left(i+1\right)=x\left(i\right)+\Delta t.Math.\frac{\mathrm{dx}}{\mathrm{dt}}{.Math.}_t& \left(6\right)\end{array}$$\begin{array}{cc}\frac{d{m}_{c,v}}{\mathrm{dt}}=f_{\mathrm{cir}}\left(c_{c,l}-c_{c,v}\right)+f_{\mathrm{dec}}+f_{c,\mathrm{add}}& \left(7\right)\end{array}$$\begin{array}{cc}\frac{d{m}_{c,l}}{\mathrm{dt}}=f_{\mathrm{cir}}\left(c_{c,v}-c_{c,l}\right)& \left(8\right)\end{array}$$\begin{array}{cc}\frac{d{m}_{o,v}}{\mathrm{dt}}=f_{\mathrm{cir}}\left(c_{o,l}-c_{o,v}\right)+\frac{M_o}{M_c}{f}_{\mathrm{dec}}-\frac{3M_o}{2M_{\mathrm{al}}}{f}_{a,\mathrm{add}}& \left(9\right)\end{array}$$\begin{array}{cc}\frac{d{m}_{o,l}}{\mathrm{dt}}=f_{\mathrm{cir}}\left(c_{o,v}-c_{o,l}\right)& \left(10\right)\end{array}$$\begin{array}{cc}\frac{p_v{V}_v}{{\mathrm{RT}}_g}\frac{d{c}_{\mathrm{co},g}}{\mathrm{dt}}=-\frac{p_v}{{\mathrm{RT}}_g}{f}_{\mathrm{dec}}-\frac{P_{atm}{v}_{\mathrm{off},25}}{{\mathrm{RT}}_{atm}}{c}_{\mathrm{co},g}& \left(11\right)\end{array}$

[0052] Note that the left side of the formula (7) expresses a change in carbon weight of the molten steel in the vacuum chamber, the first term of the three terms on the right side expresses an influence of reflux, the second term expresses an influence of the vacuum decarburization, and the third term expresses an influence of the added auxiliary raw material. In addition, the left side of the formula (8) expresses a change in carbon weight of the molten steel in the ladle, and the right side expresses an influence of reflux. In addition, the left side of the formula (9) expresses a change in dissolved oxygen weight of the molten steel in the vacuum chamber, the first term of the three terms on the right side expresses an influence of reflux, the second term expresses an influence of the vacuum decarburization, and the third term expresses an influence of the added auxiliary raw material. In addition, the left side of the formula (10) expresses a change in dissolved oxygen weight of the molten steel in a molten steel ladle, and the right side expresses an influence of reflux. In addition, the formula (11) expresses molar equilibrium between a carbon dioxide gas in a gas phase in the vacuum chamber and the carbon dioxide gas in the exhaust gas (number of moles of carbon dioxide gas molecules is in equilibrium). The expression on the left side expresses the number of moles of the carbon dioxide gas in the gas phase in the vacuum chamber, the first term on the right side expresses the number of moles of the carbon dioxide gas that flows into the gas phase from the molten steel by the vacuum decarburization, and the second term on the right side expresses the number of moles of the carbon dioxide gas diffused to the atmosphere by the exhaust gas.

[0053] In the processing of Step S4, the end point determination unit 11c4 determines whether a value acquired by addition of a predetermined value as a safety margin to the average value of the probability distribution of the carbon concentration of the molten steel (carbon concentration estimate value) calculated in the processing of Step S3 is equal to or smaller than target carbon concentration. Then, as a result of the determination, in a case where the value acquired by the addition of the predetermined value to the carbon concentration estimate value is equal to or smaller than the target carbon concentration (Step S4: Yes), the end point determination unit 11c4 determines that the vacuum decarburization treatment is completed, and ends the series of decarburization end point determination processing. On the other hand, in a case where the value acquired by the addition of the predetermined value to the carbon concentration estimate value is larger than the target carbon concentration (Step S4: No), the end point determination unit 11c4 determines that the vacuum decarburization treatment is not completed, increments the value of the program counter i by one as processing of Step S5, and then returns the decarburization end point determination processing to the processing of Step S2.

[0054] Note that although the average value of the probability distribution of the molten steel carbon concentration is used for the determination of the completion of the vacuum decarburization treatment in the present embodiment, a maximum value, or a value acquired by addition of a value acquired by multiplication of a standard deviation by 3 to the average value may be used. Furthermore, in a case where it is determined that the vacuum decarburization treatment is not completed, the decarburization end point determination processing returns from the processing of Step S4 to the processing of Step S2 via the processing of Step S5. A loop of returning from the processing of Step S4 to the processing of Step S2 is desirably performed in a cycle time of about 5 seconds.

EXAMPLES

[0055] FIG. 4 is a view illustrating a time-series change in input variables in vacuum decarburization treatment with 250 ton/charge and target carbon concentration of 13 ppm. A time-series change in the internal pressure of the vacuum chamber is illustrated in FIG. 4(a), time-series changes in the flow rates of the exhaust gas (line L1) and the argon gas (line L2) are illustrated in FIG. 4(b), and time-series changes in the mass flow rates of carbon (line L3) and aluminum (line L4) flowing into the molten steel from the added auxiliary raw material are illustrated in FIG. 4(c).

[0056] In addition, FIG. 5 is a view illustrating time-series changes in output variables of when the vacuum decarburization treatment illustrated in FIGS. 4(a) to (c) is performed. Time-series changes in the estimate value (line) and the measurement value (plot) of the molten steel carbon concentration are illustrated in FIG. 5(a), time-series changes in the estimate value (line) and the measurement value (plot) of the molten steel oxygen concentration are illustrated in FIG. 5(b), time-series changes in the estimate value (line L5) and the measurement value (line L6) of the carbon dioxide gas concentration of the exhaust gas are illustrated in FIG. 5(c), and a time-series change in the decarburization capacity coefficient is illustrated in FIG. 5(d). Note that the estimate values herein mean values estimated by a conventional method that does not learn the decarburization capacity coefficient of the model formula from online measurement values.

[0057] As illustrated in FIG. 5(a), the molten steel carbon concentration at the time of completion of the vacuum decarburization treatment is estimated to be 13 ppm in the conventional method, and this is a value smaller by 37 ppm than the measurement value of 50 ppm. In addition, as illustrated in FIG. 5(b), the molten steel oxygen concentration at 570 seconds is estimated to be 358 ppm in the conventional method, and this is a value smaller by 33 ppm than the measurement value of 391 ppm. Furthermore, as illustrated in FIG. 5(c), the carbon dioxide gas concentration of the exhaust gas at 90 seconds is estimated to be 58%, and this is a value 40% larger than the measurement value of 18%. The reason why these results are acquired is that a decarburization amount is estimated to be larger than an actual value in the conventional method. As a result, it is erroneously determined that the molten steel carbon concentration reaches the target value of 13 ppm at a time point of 669 seconds, and a decarburization treatment failure is caused since the molten steel carbon concentration does not actually reach the target value.

[0058] On the other hand, FIG. 6 is a view illustrating a result of application of the present invention to the vacuum decarburization treatment illustrated in FIGS. 4(a) to (c). Time-series changes in the estimate value (line) and the measurement value (plot) of the molten steel carbon concentration are illustrated in FIG. 6(a), time-series changes in the estimate value (line) and the measurement value (plot) of the molten steel oxygen concentration are illustrated in FIG. 6(b), time-series changes in the actual measurement value (line L5) and the measurement value (line L6) of the carbon dioxide gas concentration of the exhaust gas are illustrated in FIG. 6(c), and a time-series change in the decarburization capacity coefficient is illustrated in FIG. 6(d). In the present invention, a desirable decarburization capacity coefficient value is learned online from a difference between an online measurement value and an estimate value by the model formula of the decarburization reaction with respect to each of the carbon dioxide gas concentration of the exhaust gas and the molten steel oxygen concentration, and the learned decarburization capacity coefficient value is reflected in the model formula of the decarburization reaction. Specifically, the value of the decarburization capacity coefficient is corrected as illustrated in FIG. 6(d) in such a manner that the molten steel oxygen concentration at 570 seconds matches the measurement value as illustrated in FIG. 6(b) and the estimate value (line L5) of the carbon dioxide gas concentration of the exhaust gas matches the measurement value (line L6) as illustrated in FIG. 6(c).

[0059] Thus, as illustrated in FIG. 6(a), the molten steel carbon concentration can be estimated in such a manner as to finally match the measurement value of the molten steel carbon concentration without direct correction of the molten steel carbon concentration. This means that the value of the decarburization capacity coefficient can be learned online to be a value close to an actual value and reflected in the model formula by application of the present invention even when the value of the decarburization capacity coefficient included in the model formula of the decarburization reaction is set in a manner of being deviated from the actual value. That is, by applying the present invention, it is possible to estimate the molten steel carbon concentration with higher accuracy than the conventional method. As a result, it is possible to determine the end point of the vacuum decarburization treatment at appropriate timing without redundant treatment time. In addition, the secondary refining operation itself in the steel mill becomes highly accurate and stable operation. As a result, treatment time is shortened. Thus, variations in components of the molten steel in the steel mill are reduced, and it becomes possible to produce the molten steel highly accurately and stably.

[0060] Although an embodiment to which the invention made by the present inventors is applied has been described above, the present invention is not limited by the description and drawings included in a part of the disclosure of the present invention according to the present embodiment. That is, other embodiments, examples, operation techniques, and the like made by those skilled in the art on the basis of the present embodiment are all included in the scope of the present invention.

INDUSTRIAL APPLICABILITY

[0061] According to the present invention, a decarburization end point determination method and a decarburization end point determination device that are capable of accurately estimating carbon concentration of molten steel and accurately determining an end point of vacuum decarburization treatment can be provided. Furthermore, according to the present invention, it is possible to provide a secondary refining operation method for steel making capable of performing secondary refining operation highly accurately and stably. Furthermore, according to the present invention, a method for producing molten steel with which method the molten steel can be produced highly accurately and stably can be provided.

REFERENCE SIGNS LIST

[0062] 1 VACUUM DEGASSING FACILITY [0063] 2 LADLE [0064] 3 SUCTION PIPE [0065] 4 EXHAUST PIPE [0066] 5 VACUUM CHAMBER [0067] 6 VACUUM DEVICE [0068] 7 OXYGEN FEEDING LANCE [0069] 8 HOPPER [0070] 10 MEASUREMENT VALUE COLLECTION DEVICE