Method for Predicting Carbon Consumption and Carbon Emission of Hydrogen-Rich Blast Furnace (BF) Based on C-rd Theory

Abstract

The present disclosure provides a method for predicting carbon consumption and carbon emission of a hydrogen-rich blast furnace (BF) based on a C-r.sub.d theory. The method includes: inputting raw material and fuel conditions, product parameters, injection parameters and a heat loss verification coefficient, determining a reduction degree of hydrogen, and establishing a carbon consumption calculation model; determining a direct reduction degree, predicted values of various carbon consumption items and a predicted fuel ratio according to a carbon balance relation; calculating, if an absolute value of a difference between the predicted fuel ratio and an initial fuel ratio is less than a preset value and the various indicators satisfy constraints on the material and heat balance, and outputting the parameter condition and the corresponding calculation results of the various indicators, or otherwise adjusting the initial fuel ratio, heat loss verification coefficient, the raw material and fuel conditions, and the injection parameters.

Claims

1. A method for predicting carbon consumption and carbon emission of a hydrogen-rich blast furnace (BF) based on a C-r.sub.d theory, comprising: inputting raw material and fuel conditions, product parameters, injection parameters and a heat loss verification coefficient of a BF, wherein the raw material and fuel conditions comprise a type of iron ore, a burden structure, an initial fuel ratio, chemical components of a raw material and a fuel, and blast parameters; the blast parameters comprise a blast humidity, a blast temperature, an initial blast volume, an oxygen enrichment rate, and an oxygen (O.sub.2) purity; the product parameters comprise components of molten iron, a temperature of the molten iron, components of slag, a carbon monoxide (CO) utilization rate, components of dust, and a dust volume; and the injection parameters comprise a type of an injected hydrogen-rich medium, an injection amount of the hydrogen-rich medium, an injection temperature, an injection position, and a hydrogen (H.sub.2) utilization rate; determining a reduction degree of hydrogen according to the blast parameters, the initial fuel ratio and the injection parameters; respectively establishing relations of various carbon consumption items of the BF with a direct reduction degree according to the reduction degree of the hydrogen, the raw material and fuel conditions, the product parameters, the injection parameters and the heat loss verification coefficient to obtain a carbon consumption calculation model, wherein the various carbon consumption items comprise carburization of BF molten iron, carbon consumption in direct reduction of trace elements, carbon consumption in direct reduction of iron oxide, carbon consumption in chemical reaction, and carbon consumption for heat supply; and the carbon consumption calculation model comprises a function for calculating an amount of the carburization of the BF molten iron, a function for calculating an amount of the carbon consumption in the direct reduction of the trace elements, a function for calculating an amount of the carbon consumption in the direct reduction of the iron oxide, a function for calculating an amount of the carbon consumption in the chemical reaction, and a function for calculating an amount of the carbon consumption for the heat supply; determining the direct reduction degree according to a carbon balance relation, and determining predicted values of the various carbon consumption items and a predicted fuel ratio according to the direct reduction degree and the carbon consumption calculation model; determining whether an absolute value of a difference between the predicted fuel ratio and the initial fuel ratio is less than a preset value to obtain a first determination result; adjusting, if the first determination result indicates no, the initial fuel ratio and the initial blast volume, and going back to the step of determining a reduction degree of hydrogen according to the blast parameters, the initial fuel ratio and the injection parameters; calculating, if the first determination result indicates yes, a material balance and a heat balance of the BF according to the direct reduction degree, the reduction degree of the hydrogen and the predicted fuel ratio to obtain calculation results of various indicators, wherein the calculation results of the various indicators comprise the predicted values of the various carbon consumption items, carbon emission per ton of iron, the predicted fuel ratio, a raceway adiabatic flame temperature (RAFT) in a tuyere area, a heat loss in a high-temperature area, a temperature of top gas in a low-temperature area, and an overall heat loss of the BF; determining whether the calculation results of the various indicators satisfy constraints on a material balance and a heat balance in each of areas of the BF to obtain a second determination result, wherein the constraints on the material balance and the heat balance in each of the areas of the BF comprise a constraint on a material balance in the tuyere area, a constraint on a heat balance in the tuyere area, a constraint on the RAFT, a constraint on a material balance in the high-temperature area, a constraint on a heat balance in the high-temperature area, a constraint on a material balance in the low-temperature area, a constraint on a heat balance in the low-temperature area, a constraint on a material balance in the whole furnace, and a constraint on a heat balance in the whole furnace; adjusting, if the second determination result indicates no, the heat loss verification coefficient, the raw material and fuel conditions, and the injection parameters, and going back to the step of determining a reduction degree of hydrogen according to the blast parameters, the initial fuel ratio and the injection parameters; and outputting, if the second determination result indicates yes, the raw material and fuel conditions, the product parameters, the injection parameters, the heat loss verification coefficient and the corresponding calculation results of the various indicators, wherein the calculation results of the various indicators are used to evaluate a carbon source, a carbon direction, carbon consumption and carbon emission of the BF.

2. The method for predicting carbon consumption and carbon emission of a hydrogen-rich BF based on a C-r.sub.d theory according to claim 1, wherein the function for calculating the amount of the carbon consumption in the chemical reaction is specifically expressed by: w ( C ) i = n i ? 12 56 ( 1 - r d - r H 2 ) ? w [ Fe ] + 12 56 ? n 2 ? 1 3 ( r d + r H 2 ) ? w [ Fe ] = 1 ? CO w ( C ) CO 2 , wherein, w(C).sub.i is the amount of the carbon consumption in the chemical reaction, r.sub.d is the direct reduction degree, r.sub.H.sub.2 is the reduction degree of the hydrogen, w[Fe] is a content of ferrum (Fe) in the molten iron, ?.sub.CO is the CO utilization rate, w(C).sub.CO.sub.2 is an amount of carbon reduced to generate carbon dioxide (CO.sub.2), n.sub.1 is an excess coefficient of a reductant CO for reducing ferrous oxide (FeO) to Fe, and n.sub.2 is an excess coefficient of the reductant CO for reducing ferroferric oxide (Fe.sub.3O.sub.4) to FeO.

3. The method for predicting carbon consumption and carbon emission of a hydrogen-rich BF based on a C-r.sub.d theory according to claim 1, wherein the function for calculating the amount of the carbon consumption for the heat supply is specifically expressed by: w ( C ) = w ( C ) c ombustion + w ( C ) d , wherein, w(C) is the amount of the carbon consumption for the heat supply, w(C).sub.d is the amount of the carbon consumption in the direct reduction of the iron oxide, and w(C).sub.combustion is an amount of carbon combusted in front of a tuyere.

4. The method for predicting carbon consumption and carbon emission of a hydrogen-rich BF based on a C-r.sub.d theory according to claim 2, wherein the amount of the carbon reduced to generate the CO.sub.2 is calculated by: w ( C ) CO 2 = 1 2 5 6 [ Fe ] ( 1 - r d - r H 2 ) + 1 2 ? 1 2 5 6 w [ Fe ] .

5. The method for predicting carbon consumption and carbon emission of a hydrogen-rich BF based on a C-r.sub.d theory according to claim 3, wherein the amount of the carbon combusted in front of the tuyere is calculated by: w ( C ) combustion ( q CO + q hot blast ) + Q injection + Q others = Q reduction + Q molten iron + Q slag + Q top gas + Q H 2 O cracking + Q heat loss + Q cracking of pulverized coal + Q H 2 O evaporation + Q , wherein, q.sub.CO is heat that 1 kg of carbon is combusted in front of the tuyere to generate CO, q.sub.hot blast is heat carried in by blast air, Q.sub.injection is physical heat carried by the injected hydrogen-rich medium, Q.sub.reduction is heat consumed by reduction, Q.sub.molten iron is heat carried off by the molten iron, Q.sub.slag is heat carried off by the slag, Q.sub.top gas is heat carried off by the top gas, Q H 2 O cracking is heat consumed by cracking of water (H.sub.2O) in front of the tuyere, Q.sub.cracking of pulverized coal is heat consumed by cracking of pulverized coal in front of the tuyere, Q H 2 O evaporation is heat carried off by evaporation of H.sub.2O in the raw material and the fuel of the BF, Q.sub.heat loss is a heat loss, Q is the heat loss verification coefficient, and Q.sub.others is other heat of the BF, comprising heat carried in by hot charging of the raw material of the BF; the heat consumed by the reduction is calculated by: Q reduction = .Math. k Q reduction , k = Q direct reduction + Q indirect reduction + Q Si , Mn , P = 152190 56 w [ Fe ] ? r d + 47522 2 ? 72 w silicate + 1549 160 w Fe 2 O 3 + 20883 232 w Fe 3 O 4 + 27718 56 w [ Fe ] ? r H 2 - 13605 56 w [ Fe ] ( 1 - r d - r H 2 ) + 5087 ? w [ Mn ] + 22049 ? w [ Si ] + 15490 ? w [ P ] , wherein, k denotes different reduction items, Q.sub.direct reduction is heat consumed by the direct reduction of the iron oxide, Q.sub.indirect reduction is heat consumed by indirect reduction of the iron oxide, Q.sub.Si,Mn,P is heat consumed by the reduction of the trace elements, w.sub.silicate is a mass of silicon dioxide (SO.sub.2) in an iron bearing material, w.sub.Fe.sub.2.sub.O.sub.3 is a mass of ferric oxide (Fe.sub.2O.sub.3) in the iron bearing material, w.sub.Fe.sub.3.sub.O.sub.4 is a mass of Fe.sub.3O.sub.4 in the iron bearing material, w[Fe] is a content of Fe in the molten iron, w[Mn] is a content of manganese (Mn) in the molten iron, w[Si] is a content of silicon (Si) in the molten iron, w[P] is a content of phosphorus (P) in the molten iron, r.sub.d is the direct reduction degree, and r.sub.H.sub.2 is the reduction degree of the hydrogen; and the heat carried off by the top gas is calculated by: Q top gas = Q CO 2 , top gas + Q CO , top gas + Q N 2 , top gas + Q H 2 , top gas + Q H 2 O , top gas = q CO 2 , top gas 12 w ( C ) CO 2 + q CO , top gas 12 ( w ( C ) i - w ( C ) CO 2 ) + q N 2 , top gas 22.4 ? ( N 2 ) b ? V b + q H 2 , top gas 2 w H 2 ? H 2 + q H 2 O , top gas 2 ? w H 2 ( 1 - ? H 2 ) wherein, Q.sub.CO.sub.2.sub.,top gas is heat carried off by CO.sub.2 in the top gas, Q.sub.CO,top gas is heat carried off by CO in the top gas, Q N 2 , top gas is heat carried off by nitrogen (N.sub.2) in the top gas, Q H 2 , top gas is heat carried off by H.sub.2 in the top gas, Q H 2 O , top gas is heat carried off by H.sub.2O in the top gas, q.sub.CO.sub.2.sub.,top gas is heat carried off by 1 mol of CO.sub.2 in the top gas, q.sub.CO,top gas is heat carried off by 1 mol of CO in the top gas, Q N 2 , top gas is heat carried off by 1 mol of N.sub.2 in the top gas, Q H 2 , top gas is heat carried off by 1 mol of H.sub.2 in the top gas, Q H 2 O , top gas is heat carried off by 1 mol of H.sub.2O in the top gas, w(C).sub.i is the amount of the carbon consumption in the chemical reaction, w(C).sub.CO.sub.2 is the amount of the carbon reduced to generate CO.sub.2, V.sub.b is the initial blast volume, (N.sub.2).sub.b is a volume fraction of N.sub.2 in the blast air, ?.sub.H.sub.2 is the H.sub.2 utilization rate, and w.sub.H.sub.2 is a total volume of H.sub.2 charged to the furnace.

6. The method for predicting carbon consumption and carbon emission of a hydrogen-rich BF based on a C-r.sub.d theory according to claim 1, wherein the function for calculating the amount of the carbon consumption in the direct reduction of the iron oxide is specifically expressed by: w ( C ) d = 1 2 5 6 ? w [ Fe ] ? r d , wherein, w(C).sub.d is the amount of the carbon consumption in the direct reduction of the iron oxide, w[Fe] is a content of Fe in the molten iron, and r.sub.d is the direct reduction degree.

7. The method for predicting carbon consumption and carbon emission of a hydrogen-rich BF based on a C-r.sub.d theory according to claim 1, wherein the carbon balance relation is specifically expressed by: w ( C ) i = C injection + w ( C ) + w ( C ) carburization + w ( C ) S i , Mn , P , wherein, w(C).sub.i is the amount of the carbon consumption in the chemical reaction, C.sub.injection is an equivalent amount of carbon carried in by the injected hydrogen-rich medium, w(C) is the amount of the carbon consumption for the heat supply, w(C).sub.carburization is the amount of the carburization of the BF molten iron, and w(C).sub.Si,Mn,P is the amount of the carbon consumption in the direct reduction of the trace elements.

8. The method for predicting carbon consumption and carbon emission of a hydrogen-rich BF based on a C-r.sub.d theory according to claim 1, wherein the predicted fuel ratio is calculated by: K = ( w ( C ) + w ( C ) Si , Mn , P + w ( C ) carburization ) / w ( C ) k , wherein, K is the predicted fuel ratio, the w(C) is the amount of the carbon consumption for the heat supply, w(C).sub.carburization is the amount of the carburization of the BF molten iron, w(C).sub.Si,Mn,P is the amount of the carbon consumption in the direct reduction of the trace elements, and w(C).sub.k is a mass percent of carbon in the fuel.

9. The method for predicting carbon consumption and carbon emission of a hydrogen-rich BF based on a C-r.sub.d theory according to claim 1, wherein the adjusting the initial fuel ratio and the initial blast volume specifically comprises: updating the initial fuel ratio with a sum of the predicted fuel ratio and the preset value; and determining a predicted blast volume according to a carbon-oxygen balance in front of a tuyere, and updating the initial blast volume with the predicted blast volume, wherein the predicted blast volume is calculated by: V b = w ( C ) c o mbustion ? 22.4 24 ? v ( O 2 ) b , wherein, V.sub.b is the predicted blast volume, w(C).sub.combustion is an amount of carbon combusted in front of the tuyere, and v(O.sub.2).sub.b is a volume fraction of O.sub.2 in blast air.

10. The method for predicting carbon consumption and carbon emission of a hydrogen-rich BF based on a C-r.sub.d theory according to claim 1, after the respectively establishing relations of various carbon consumption items of the BF with a direct reduction degree according to the reduction degree of the hydrogen, the raw material and fuel conditions, the product parameters, the injection parameters and the heat loss verification coefficient to obtain a carbon consumption calculation model, further comprising: drawing and outputting a carbon consumption line distribution diagram according to the carbon consumption calculation model.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

[0048] To describe the technical solutions in the examples of the present disclosure or in the prior art more clearly, the accompanying drawings required for the examples are briefly described below. Apparently, the accompanying drawings in the following description show merely some examples of the present disclosure, and those of ordinary skill in the art may still derive other accompanying drawings from these accompanying drawings without creative efforts.

[0049] FIG. 1 is a method for predicting carbon consumption and carbon emission of a hydrogen-rich BF based on a C-r.sub.d theory according to the present disclosure;

[0050] FIG. 2 illustrates a carbon consumption line distribution diagram of a BF in all-coke smelting and coal injection according to the present disclosure;

[0051] FIG. 3 illustrates a carbon consumption line distribution diagram of a BF in single COG injection according to the present disclosure; and

[0052] FIG. 4 illustrates a carbon consumption line distribution diagram of a BF in injection of coal+COG according to the present disclosure.

DETAILED DESCRIPTION OF THE EMBODIMENTS

[0053] The technical solutions of the embodiments of the present disclosure are clearly and completely described below with reference to the accompanying drawings. Apparently, the described embodiments are merely a part rather than all of the embodiments of the present disclosure. All other embodiments obtained by those of ordinary skill in the art based on the embodiments of the present disclosure without creative efforts shall fall within the protection scope of the present disclosure.

[0054] An objective of the present disclosure is to provide a method for predicting carbon consumption and carbon emission of a hydrogen-rich BF based on a C-r.sub.d theory, so as to obtain relations of operation parameters in different hydrogen enrichment processes with carbon consumption and carbon emission of the BF, predict and evaluate carbon consumption and carbon emission of the different hydrogen enrichment processes, and elaborate a carbon source, a carbon direction, carbon consumption and carbon emission in the operation parameters of the different hydrogen enrichment processes, thereby providing theoretical guidance for hydrogen enrichment of the BF and low-carbon energy transformation of the BF.

[0055] The present disclosure uses the following technical solutions: According to a principle of a C-r.sub.d relation proposed by A. H. Ramm (namely a C-r.sub.d theory), and in view of influences of hydrogen, linear relations of carbon consumption items in a BF with a direct reduction degree r.sub.d are established. The carbon consumption items in the BF include carburization of molten iron, carbon consumption in reduction of trace elements, carbon consumption in direct reduction of iron, carbon consumption in indirect reduction, and carbon consumption for heat supply. The relations of the various carbon consumption items with the r.sub.d are calculated. Loop iteration is carried out according to a carbon balance relation to obtain the r.sub.d, r.sub.H.sub.2 and a predicted fuel ratio. The r.sub.d, the r.sub.H.sub.2 and the predicted fuel ratio are used to calculate a material balance and a heat balance in each of areas of the BF. According to constraints on the material balance and the heat balance in the area, a heat loss verification coefficient, raw material and fuel conditions, and injection parameters are adjusted constantly, thereby predicting various carbon consumption items and carbon emission of the BF in different hydrogen enrichment processes.

[0056] To make the above objectives, features, and advantages of the present disclosure clearer and more comprehensible, the present disclosure will be further described in detail below with reference to the accompanying drawings and the specific implementations.

[0057] As shown in FIG. 1, the present disclosure provides a carbon evaluation method for a hydrogen-rich BF, including the following steps: [0058] Step S1: Input raw material and fuel conditions, product parameters, injection parameters and a heat loss verification coefficient of a BF, where the raw material and fuel conditions include a type of iron ore, a burden structure, an initial fuel ratio, chemical components of a raw material and a fuel, and blast parameters; the blast parameters include a blast humidity, a blast temperature, an initial blast volume, an oxygen enrichment rate, and a O.sub.2 purity; the product parameters include components of molten iron, a temperature of the molten iron, components of slag, a CO utilization rate, components of dust, and a dust volume; and the injection parameters include a type of an injected hydrogen-rich medium, an injection amount of the hydrogen-rich medium, an injection temperature, an injection position, and a H.sub.2 utilization rate. [0059] Step S2: Determine a reduction degree of hydrogen according to the blast parameters, the initial fuel ratio and the injection parameters. [0060] Step S3: Respectively establish relations of various carbon consumption items of the BF with a direct reduction degree according to the reduction degree of the hydrogen, the raw material and fuel conditions, the product parameters, the injection parameters and the heat loss verification coefficient to obtain a carbon consumption calculation model, where the various carbon consumption items include carburization of BF molten iron, carbon consumption in direct reduction of trace elements, carbon consumption in direct reduction of iron oxide, carbon consumption in chemical reaction, and carbon consumption for heat supply; and the carbon consumption calculation model includes a function for calculating an amount of the carburization of the BF molten iron, a function for calculating an amount of the carbon consumption in the direct reduction of the trace elements, a function for calculating an amount of the carbon consumption in the direct reduction of the iron oxide, a function for calculating an amount of the carbon consumption in the chemical reaction, and a function for calculating an amount of the carbon consumption for the heat supply.

[0061] Specifically, according to the raw material and fuel conditions (namely raw material conditions and fuel conditions) of the BF and product requirements, relations of the carburization of the BF molten iron, the carbon consumption in the direct reduction of the trace elements, the carbon consumption in the direct reduction of the iron oxide, the carbon consumption in the chemical reaction and the carbon consumption for the heat supply with the direct reduction degree are established. The relation of each carbon consumption item C and the direct reduction degree r.sub.d is calculated as follows: [0062] (1) The amount w(C).sub.carurtization of the carburization of the molten iron:

[00021] w ( C ) carburization = 1.34 + 2. 54 ? 1 0 - 3 t + 0.04 ? w [ Mn ] - 0.3 w [ Si ] - 0.35 w [ P ] - 0.4 w [ S ] ( 1 )

[0063] where, t is the temperature of the molten iron; w[Mn], w[Si], w[P] and w[S] are respectively contents of elements in the molten iron by a mass percent; 2.54, 0.04, 0.3, 0.35 and 0.4 are respectively relational coefficients of the amount of the carburization with the temperature, w[Mn], w[Si], w[P], and w[S] in the molten iron. [0064] (2) The amount w(C).sub.Si,Mn,P of the carbon consumption in the direct reduction of the trace elements:

[00022] w ( C ) Si , Mn , P = 1 2 2 8 ? 2 ? w [ Si ] + 1 2 5 6 ? w [ Mn ] + 6 0 6 2 ? w [ P ] ( 2 )

[0065] where, 12/28?2, 12/56 and 60/62 are respectively masses of carbon consumed to reduce 1 mol of [Si], 1 mol of [Mn] and 1 mol of [P]. [0066] (3) The amount w(C).sub.d of the carbon consumption in the direct reduction of the iron oxide:

[00023] w ( C ) d = 1 2 5 6 ? w [ Fe ] ? r d ( 3 )

[0067] where, r.sub.d is the direct reduction degree defined by the Pavlov, and specifically refers to a direct reduction degree of [Fe]. [0068] (4) The amount of the carbon consumption in the chemical reaction, namely an amount w(C).sub.i of carbon consumption in indirect reduction of the iron oxide:

[00024] w ( C ) i = n 1 ? 1 2 5 6 ( 1 - r d - r H 2 ) ? w [ Fe ] + 1 2 5 6 ? n 2 ? 1 3 ( r d + r H 2 ) ? w [ Fe ] = 1 ? CO w ( C ) CO 2 ( 4 ) w ( C ) CO 2 = 1 2 5 6 w [ Fe ] ( 1 - r d - r H 2 ) + 1 2 ? 1 2 5 6 w [ Fe ] + w ( C ) Si , Mn , P ( 5 ) r H 2 = w H 2 ? H 2 ? 5 6 2 ? 1 w [ F e ] = [ K ? ( H 2 ) fuel + V b ? ( H 2 ) b + V injection ? ( H 2 ) injection ] ? ? H 2 ? 5 6 2 ? 1 w [ Fe ] ( 6 )

[0069] where, r.sub.H.sub.2 is the reduction degree of the hydrogen, ?.sub.H.sub.2 is the H.sub.2 utilization rate, w.sub.H.sub.2 is a total volume of H.sub.2 charged to the furnace, specifically including H.sub.2 from cracking of H.sub.2O in the raw material, the fuel and blast air as well as injected H.sub.2, K is the initial fuel ratio, V.sub.b is the initial blast volume, V.sub.injection is a mass of an injected hydrogen-rich medium, (H.sub.2).sub.fuel, (H.sub.2).sub.b, and (H.sub.2).sub.injection are respectively equivalent H.sub.2 volumes in the fuel, the blast air and the injected hydrogen-rich medium, ?.sub.CO is a CO utilization rate in coal gas, w(C).sub.CO.sub.2 is an amount of carbon reduced to generate CO.sub.2, n.sub.1 is an excess coefficient of a reductant CO for reducing FeO to Fe, and n.sub.2 is an excess coefficient of the reductant CO for reducing Fe.sub.3O.sub.4 to FeO. [0070] (5) The amount w(C) of the carbon consumption from the heat supply:

[00025] w ( C ) = w ( C ) combustion + w ( C ) d ( 7 )

[0071] where, w(C) is the amount of the carbon consumption for the heat supply, and w(C).sub.combustion is an amount of carbon combusted in front of a tuyere and may be calculated through a second heat balance of the BF, namely ?.sub.i Q.sub.heatsupply.sub.i=?.sub.j Q.sub.heat consumption j, ?.sub.i Q.sub.heat supply i being a total amount of heat supply per ton of iron, and ?.sub.j Q.sub.heat consumption j being a total amount of heat consumption per ton of iron, at a unit of kJ/tFe. Specifically:

[00026] w ( C ) combustion ( q CO + q hot blast ) + Q injection + Q others = Q reduction + Q molten iron + Q slag + Q top gas + Q H 2 O cracking + Q heat loss + Q cracking of pulverized coal + Q H 2 O evaporation + Q ( 8 )

[0072] where, q.sub.CO is heat that 1 kg of carbon is combusted in front of the tuyere to generate CO; q.sub.hot blast heat carried in by hot-blast air required to combust 1 kg of carbon; Q.sub.injection is physical heat carried by the injected hydrogen-rich medium, which is associated with the injection amount of the hydrogen-rich medium and the injection temperature, and calculated by Eq. (9); Q.sub.others is other heat of the BF, such as heat carried in by hot charging of the raw material of the BF; Q.sub.reduction is heat consumed by reduction, which is associated with the direct reduction degree, the reduction degree of the hydrogen and the components of the molten iron, and calculated by Eq. (10); Q.sub.molten iron is heat carried off by the molten iron; Q.sub.slag is heat carried off by the slag; Q.sub.top gas is heat carried off by top gas, which is associated with components and a temperature of the top gas and calculated by Eq. (11), the components of the top gas depending on the direct reduction degree and the reduction degree of the hydrogen;

[00027] Q H 2 0 cracking

is heat consumed by cracking of H.sub.2O in front of the tuyere; Q.sub.heat loss is a heat loss, which is specifically 5-10% of total heat; Q.sub.cracking of pulverized coal is heat consumed by cracking of pulverized coal in front of the tuyere;

[00028] Q H 2 0 evaporation

is heat carried off by evaporation of H.sub.2O in the raw material and the fuel of the BF; and Q is the heat loss verification coefficient. As can be seen, except that Q.sub.reduction and Q.sub.top gas are associated with the direct reduction degree and the reduction degree of the hydrogen, other items can be solved by the raw material and fuel conditions of the BF, the product parameters and the injection parameters to obtain specific values. The reduction degree of the hydrogen can also be obtained by Eq. (6). Hence, a relation between w(C).sub.combustion and the direct reduction degree r.sub.d can be established, and substituted into Eq. (7) to obtain a relation between w(C) and the r.sub.d.

[00029] Q inject.Math.on = .Math. i ? 2 9 8 T injection c P , i dT ( 9 )

[0073] where, T.sub.injection is the injection temperature, i is the injected hydrogen-rich medium, and c.sub.P,i is a molar heat capacity at a constant pressure for the injected hydrogen-rich medium i.

[00030] Q reduction = .Math. k Q reduction , k = Q direct reduction + Q indirect reduction + Q Si , Mn , P = 152190 56 w [ Fe ] ? r d + 47522 2 ? 72 w silicate + 1549 160 w Fe 2 O 3 + 20883 232 w Fe 3 O 4 + 27718 56 w [ Fe ] ? r H 2 - 13605 56 w [ Fe ] ( 1 - r d - r H 2 ) + 5087 ? w [ Mn ] + 22049 ? w [ Si ] + 15490 ? w [ P ] ( 10 )

[0074] where, k denotes different reduction items; Q.sub.direct reduction is heat consumed by the direct reduction of the iron oxide, Q.sub.indirect reduction is heat consumed by the indirect reduction of the iron oxide, and Q.sub.Si,Mn,P is heat consumed by the reduction of the trace elements, at a unit of kJ/tFe; w.sub.silicate is a mass of SiO.sub.2 in an iron bearing material, w.sub.Fe.sub.2.sub.O.sub.3 is a mass of Fe.sub.2O.sub.3 in the iron bearing material, and w.sub.Fe.sub.3.sub.O.sub.4 is a mass of Fe.sub.3O.sub.4 in the iron bearing material, at a unit of kg/tFe; 152190/56 is heat consumed to indirectly reduce 1 mol of FeO, 47522/2?72 is heat consumed to reduce SiO.sub.2, 1549/160 is heat consumed to reduce Fe.sub.2O.sub.3, 20883/232 is heat consumed to reduce Fe.sub.3O.sub.4, 27718/56 is heat consumed to reduce the hydrogen, 5087 is heat consumed to reduce Mn, 22049 is heat consumed to reduce Si, 15490 is heat consumed to reduce P, and ?13605/56 is heat released by reducing 1 mol of FeO with CO.

[00031] Q top gas = Q CO 2 , top gas + Q CO , top gas + Q N 2 , top gas + Q H 2 , top gas + Q H 2 O , top gas = q CO 2 , top gas 12 w ( C ) CO 2 + q CO , top gas 12 ( w ( C ) i - w ( C ) CO 2 ) + q N 2 , top gas 22.4 ? ( N 2 ) b ? V b + q H 2 , top gas 2 w H 2 ? H 2 + q H 2 O , top gas 2 ? w H 2 ( 1 - ? H 2 ) ( 11 )

[0075] where, Q.sub.CO.sub.2.sub.,top gas, Q.sub.CO,top gas,

[00032] Q N 2 , top gas , Q H 2 , top gas , Q H 2 O , top gas

respectively heat carried off by CO.sub.2, CO, N.sub.2, H.sub.2 and H.sub.2O in the top gas, at a unit of kJ/tFe; q.sub.CO.sub.2.sub.,top gas, q.sub.CO,top gas,

[00033] Q N 2 , top gas , Q H 2 , top gas , Q H 2 O , top gas

are respectively heat carried off by 1 mol of CO.sub.2, 1 mol of CO, 1 mol of N.sub.2, 1 mol of H.sub.2 and 1 mol of H.sub.2O in the top gas, at a unit of kJ/mol, and calculated by Eq. (12); V.sub.b is the initial blast volume; and (N.sub.2).sub.b is a volume fraction of N.sub.2 in the blast air.

[00034] q j , top gas = ? 2 9 8 T top gas c P , j d T ( 12 )

[0076] where, j is a component of the top gas, mainly including CO.sub.2, CO, N.sub.2, H.sub.2 and H.sub.2O, c.sub.P,j is a molar heat capacity at a constant pressure for the component j, for example: c.sub.P,O.sub.2 is a molar heat capacity at a constant pressure for O.sub.2, and c.sub.P,N.sub.2 is a molar heat capacity at a constant pressure for N.sub.2; and T.sub.top gas is a temperature of the top gas.

[0077] After the relations of the various carbon consumption items with the direct reduction degree are established, the carbon consumption calculation model is obtained. When the direct reduction degree is within [0, 1?r.sub.H.sub.2], a carbon consumption line distribution diagram is drawn with the direct reduction degree as a horizontal axis and the carbon consumption as a longitudinal coordinate. It is to be noted that the carbon consumption in the direct reduction, the carbon consumption in the chemical reaction and the carbon consumption for the heat supply, except the carburization of the molten iron, are drawn with w(C).sub.carburization+w(C).sub.Si,Mn,P as a starting point, namely after the relation between the carbon consumption and the direct reduction degree is obtained, a constant term (w(C).sub.carbrization+w(C).sub.Si,Mn,P) is added. [0078] Step S4: Determine the direct reduction degree according to a carbon balance relation, and determine predicted values of the various carbon consumption items and a predicted fuel ratio according to the direct reduction degree and the carbon consumption calculation model.

[0079] The carbon balance relation is specifically expressed by:

[00035] w ( C ) i = C injection + w ( C ) + w ( C ) carburiztion + w ( C ) S i , Mn , P ( 13 )

[0080] where, w(C).sub.i is the amount of the carbon consumption in the chemical reaction, C.sub.injection is an equivalent amount of carbon carried in by the injected hydrogen-rich medium, w(C) is the amount of the carbon consumption for the heat supply, w(C).sub.carburization is the amount of the carburization of the BF molten iron, and w(C).sub.Si,Mn,P is the amount of the carbon consumption in the direct reduction of the trace elements.

[0081] Specifically, according to the carbon balance relation, equations for calculating C.sub.injection, w(C), w(C).sub.carburization and w(C).sub.Si,Mn,P are substituted into Eq. (13), namely Eqs. (1)-(12) are substituted into Eq. (13). On two sides of the resulting equation, only r.sub.d is unknown. Solving the equation obtains the direct reduction degree r.sub.d. The r.sub.d is brought back to the carbon consumption calculation model, thereby obtaining the various carbon consumption items of the BF.

[0082] According to the carbon-oxygen balance in front of the tuyere, the predicted fuel ratio is calculated by:

[00036] K = ( w ( C ) + w ( C ) S i , Mn , P + w ( C ) carburization ) / w ( C ) k ( 14 )

[0083] where, K is the predicted fuel ratio, and w(C).sub.k is a mass percent of carbon in the fuel. [0084] Step S5: Determine whether an absolute value of a difference between the predicted fuel ratio and the initial fuel ratio is less than a preset value to obtain a first determination result. Preferably, the preset value is 0.2. [0085] Step S6: Adjust, if the first determination result indicates no, the initial fuel ratio and the initial blast volume, and go back to Step S2.

[0086] Specifically, the initial fuel ratio is updated with a sum of the predicted fuel ratio and the set value. According to the carbon-oxygen balance in front of the tuyere, a predicted blast volume is determined, and the initial blast volume is updated with the predicted blast volume. [0087] Step S7: Calculate, if the first determination result indicates yes, a material balance and a heat balance of the BF according to the direct reduction degree, the reduction degree of the hydrogen and the predicted fuel ratio to obtain calculation results of various indicators, where the calculation results of the various indicators include the predicted values of the various carbon consumption items, carbon emission per ton of iron, the predicted fuel ratio, an RAFT in a tuyere area, a heat loss in a high-temperature area, a temperature of top gas in a low-temperature area, and an overall heat loss of the BF. Generally, the temperature of 900-1,000? C. in the BF is used as a boundary between the high-temperature area and the low-temperature area. The high-temperature area refers to an area with a temperature higher than the boundary, and the low-temperature area refers to an area with a temperature lower than the boundary.

[0088] In addition, according to the actual production requirements, the calculation results of the various indicators may further include a distribution of the reduction degree, a coal generation process, a heat balance sheet in the high-temperature area, a heat balance sheet in the low-temperature area, and an overall heat balance sheet and an overall material balance sheet of the BF. Based on the calculation results of the various indicators, the present disclosure can further be used to analyze trends of the various indicators under different hydrogen enrichment conditions to obtain adjustment directions. [0089] Step S8: Determine whether the calculation results of the various indicators satisfy constraints on a material balance and a heat balance in each of areas of the BF to obtain a second determination result, where the constraints on the material balance and the heat balance in each of the areas of the BF include a constraint on a material balance in the tuyere area, a constraint on a heat balance in the tuyere area, a constraint on the RAFT, a constraint on a material balance in the high-temperature area, a constraint on a heat balance in the high-temperature area, a constraint on a material balance in the low-temperature area, a constraint on a heat balance in the low-temperature area, a constraint on a material balance in the whole furnace, and a constraint on a heat balance in the whole furnace. [0090] Step S9: Adjust, if the second determination result indicates no, the heat loss verification coefficient, the raw material and fuel conditions, and the injection parameters, and go back to Step S2. [0091] Step S10: Output, if the second determination result indicates yes, the raw material and fuel conditions, the product parameters, the injection parameters, the heat loss verification coefficient and the corresponding calculation results of the various indicators, where the calculation results of the various indicators are used to evaluate a carbon source, a carbon direction, carbon consumption and carbon emission of the BF.

[0092] Specifically, the r.sub.d, r.sub.H2 and various carbon consumption values are substituted into calculation of the material balance and the heat balance in each of the areas of the BF to observe whether the constraints on the RAFT, the high-temperature area, the low-temperature area, the heat balance in the whole furnace and the like are satisfied. If yes, the results are output. If no, the heat loss verification coefficient in the function for calculating the amount of the carbon consumption for the heat supply is adjusted, and the input raw material and fuel conditions and injection parameters are adjusted, until the constraints are satisfied. It is to be noted that upon adjustment, the heat loss accounts for 5-10% of total heat.

[0093] Further, after the step of respectively establishing relations of various carbon consumption items of the BF with a direct reduction degree according to the reduction degree of the hydrogen, the raw material and fuel conditions, the product parameters, the injection parameters and the heat loss verification coefficient to obtain a carbon consumption calculation model, the method further includes: Draw and output a carbon consumption line distribution diagram according to the carbon consumption calculation model. From the carbon consumption line distribution diagram, the carbon source and the carbon direction can be directly known.

[0094] Specific examples are provided hereinafter to verify correctness of the method.

[0095] Example 1: Relations between various carbon consumption items C of the hydrogen-rich BF and r.sub.d are established, and intersections are obtained, thereby verifying the calculation method.

[0096] Production parameters from a BF in China are used as an initial condition for calculation. The BF has a CO utilization rate of 50%, thus obtaining the carbon consumption w(C).sub.i=w(C).sub.?CO=0.5=1/0.5*w(C).sub.CO.sub.2 in chemical reaction. The coal ratio of the BF is fixed at 167 kg/tFe, and the coke ratio (fuel ratio of the BF=fuel ratio+coal ratio) in the fuel ratio of the BF is predicted. With the carbon consumption calculation model of the present disclosure, the coke ratio is 358 kg/tFe. However, the coke ratio of the BF in actual production is 355 kg/tFe. Therefore, the error is less than 1%, which verifies the correctness of the model. Through calculation of the model, the carbon consumption line distribution diagram of the BF in all-coke smelting and coal injection is as shown in FIG. 2. As can be seen from FIG. 2, the carbon consumption of the BF in the all-coke smelting is the crease line OPQ. In the section OP, the carbon consumption of the BF depends on the chemical reaction. In the section PQ, the carbon consumption of the BF depends on the heat supply. In the coal injection of the BF, pulverized coal can carry in a part of hydrogen-containing substances, such that the reduction degree of the hydrogen in the BF is increased, the direct reduction degree r.sub.d is lowered, the carbon consumption in the indirect reduction is reduced, the lines w(C).sub.CO2, w(C).sub.?CO=0.5 move down, and the total carbon consumption of the BF is reduced from 437 kg/tFe at the point P to 433 kg/tFe at the point M. On the other hand, the pulverized coal can particulate in combustion and reduction instead of a part of coke to effectively reduce the coke ratio. As can be seen from FIG. 2, the carbon of the BF in the coal injection comes from the following sources: 124 kg/tFe from the coal injection, and 309 kg/tFe from the coke ratio. Directions of the carbon are as follows: 43.40 kg/tFe for the carburization, 5.83 kg/tFe for the direct reduction of the trace elements, 97.10 kg/tFe for the direct reduction, and 286.59 kg/tFe for the combustion. The total carbon consumption is 433 kg/tFe, and the carbon emission in indirect reduction for generating CO.sub.2 is 191.84 kg/tFe.

[0097] Example 2: The coal injection is changed into single COG injection to predict influences of hydrogen enrichment with the single COG injection on carbon consumption of the BF.

[0098] 300 m.sup.3/tHM COG at 25? C. is used as an injection condition. Through calculation of the model, the C-r.sub.d line distribution is as shown in FIG. 3. As can be seen from FIG. 3, due to a high content of hydrogen in the COG, the reduction degree of the hydrogen is increased, the carbon consumption w(C).sub.d in the direct reduction and the carbon consumption w(C).sub.CO2 in the indirect reduction are significantly reduced, and the carbon consumption w(C) for the heat supply rises slightly. From FIG. 3, with the COG injection, the intersection between the carbon consumption in the chemical reaction and the carbon consumption for the heat supply is changed from D to D, the direct reduction degree is changed from b to a, and the total carbon consumption of the BF is reduced from 433 kg/tFe to 405 kg/tFe. 300 Nm.sup.3/tFe of COG are injected, such that about 57 kg of carbon is carried in to participate in the indirect reduction or provide excessive CO for the indirect reduction. In the hydrogen-rich injection, the carbon of the BF comes from the following sources: 57 kg/tFe from the COG, and 348 kg/tFe from the coke ratio. Directions of the carbon are as follows: 43.40 kg/tFe for the carburization, 5.83 kg/tFe for the direct reduction of the trace elements, 52.66 kg/tFe for the direct reduction, 246.47 kg/tFe for the combustion, and 57 kg/tFe as the excessive carbon for the indirect reduction. The total carbon consumption is 405 kg/tFe, and the carbon emission in indirect reduction for generating CO.sub.2 is 178.12 kg/tFe. In the hydrogen-rich injection, the carbon consumption for the heat supply is 348 kg/tFe. Assuming that the carbon consumption for the heat supply is only provided by the coke, the coke ratio is 405 kg/tHM, which is significantly higher than that of the coal injection. Therefore, the single COG injection for the hydrogen enrichment can achieve energy conservation and emission reduction of the BF, but cannot reduce the coke ratio of the BF. This indicates that the single COG injection should be combined with the coal injection.

[0099] Example 3: The coal injection is changed into hybrid injection of coal and COG to predict influences of hydrogen enrichment with the hybrid injection of the coal and the COG on carbon consumption of the BF.

[0100] By fixing the coal injection rate at 167 kg/tFe, carbon consumption line distribution of the BF with the hybrid injection of the coal and the COG is as shown in FIG. 4. As can be seen from FIG. 4, the carbon consumption per ton of iron for the BF with the hybrid injection of the coal and the COG is reduced from 433 kg to 398 kg. Meanwhile, the injected hydrogen-containing carbonaceous substances (pulverized coal and COG) participate in the indirect reduction and heat supply instead of a part of coke. In the hydrogen-rich injection, the carbon of the BF comes from the following sources: 57 kg/tFe from the COG, 124 kg/tFe from the pulverized coal, and 217 kg/tFe from the coke ratio. Directions of the carbon are as follows: 43.40 kg/tFe for the carburization, 5.83 kg/tFe for the direct reduction of the trace elements, 41.3 kg/tFe for the direct reduction, 250.32 kg/tFe for the combustion, and 57 kg/tFe as the excessive carbon for the indirect reduction. The total carbon consumption is 398 kg/tFe, and the carbon emission in indirect reduction for generating CO.sub.2 is 174.37 kg/tFe. In the hydrogen-rich injection, the coke ratio is reduced to 252 kg/tFe, which satisfies requirements of the BF on energy conservation and emission reduction. However, by substituting the r.sub.d, r.sub.H2 and predicted fuel ratio in the condition to calculation of the material balance and the heat balance in each of the areas, the RAFT is reduced from 2134? C. to 1642? C. due to little physical heat in the injected COG and a fact that the CH.sub.4 in the COG is cracked dramatically in the tuyere area to cause volume expansion. This indicates that heat at a lower side of the BF is of importance to the BF with the hybrid injection of the coal and the COG, and the heat loss verification coefficient, the raw material and fuel conditions and the injection parameters are further to be adjusted to satisfy the constraints on the material balance and the heat balance in each of the areas of the BF.

[0101] Example 4: While the constraints on the material balance and the heat balance in each of the areas of the BF are satisfied, energy conservation and carbon emission effects of different hydrogen enrichment processes that are predicted with the method for predicting carbon consumption and carbon emission of a hydrogen-rich BF based on a C-r.sub.d theory are as shown in Table 1. As compared with results of the BF from foreign hydrogen-enrichment test, the present disclosure is reasonable and reliable.

TABLE-US-00001 TABLE 1 Comparison of different hydrogen enrichment processes in energy conservation and carbon emission Injection amount Coke ratio CO.sub.2 emission Comparison at Condition per ton of iron kg/tFe RAFT ? C. reduction % home and abroad Coal injection 167 kg/tFe 358 2134 1% BF in China Single COG 119 2101 3% Japanese COURSE50 injection Nm.sup.3/tFe project Single COG 119 2097 3.17%.sup. Present disclosure injection Nm.sup.3/tFe Single H.sub.2 310 2091 12% Japanese COURSE50 injection Nm.sup.3/tFe project Single H.sub.2 310 2082 11% Present disclosure injection Nm.sup.3/tFe Top Circulation volume 273 1852 22.67% Present disclosure gas-gas-making 775 furnace-BF top Nm.sup.3/tFe gas circulation Injection of coal + 167 kg/tFe of coal + 305-280 1900 0-5% Present disclosure COG 132.5-357.5 Nm.sup.3/tFe of COG Injection of coal + 167 kg/tFe of coal 316-300 1900 9%-13% Present disclosure H.sub.2 167-500 Nm.sup.3/tFe of H.sub.2 Injection of coal + 167 kg/tFe of coal + 272-182 1900 19.5%-39.7% Present disclosure 1200? C. H.sub.2 295-800 Nm.sup.3/tFe of H.sub.2

[0102] Each example of the present disclosure is described in a progressive manner, each example focuses on the difference from other examples, and the same and similar parts between the examples may refer to each other.

[0103] Specific examples are used herein to explain the principles and embodiments of the present disclosure. The foregoing description of the embodiments is merely intended to help understand the method of the present disclosure and its core ideas; besides, various modifications may be made by those of ordinary skill in the art to specific embodiments and the scope of application in accordance with the ideas of the present disclosure. In conclusion, the content of the present specification shall not be construed as limitations to the present disclosure.