Part 1. Characteristics of Technology and Uneven ProcessesCognition, Calculatio
开篇：润墨网以专业的文秘视角，为您筛选了一篇Part 1. Characteristics of Technology and Uneven ProcessesCognition, Calculatio范文，如需获取更多写作素材，在线客服老师一对一协助。欢迎您的阅读与分享！
[a]Iron and Steel Institute of the Ukraine National Academy of Science， Academician Starodubov Sqv.， Dnepropetrovsk， Ukraine.
*Corresponding auhtor.
Received 1 June 2013； accepted 22 Augutst 2013
Abstract
Blast melting is one of the few industrial technologies that preserve the essence and significance by all technical revolutions. This phenomenon exists due to certain properties of the system that ensures exponentially increasing the productivity and linearly lowering the coke rate， that seeks to 200-250 kg/thm. The solution of problems of blast-furnace smelting involves the solution of two analytical tasks： study of the relationship of real parameters and characteristics of the blast melting； forecast of expected parameters and processes on preset parameters of work of the blast melting.
The first task is solved on the basis of balance equations of conservation of mass and energy， the second-based on the method of numerical modeling of processes in radial annular cross-sections along the height of the furnace： multi-zone model of heat-and mass transfer； physico-chemical transformations and mechanics of material and gases. During the numerical and analytical investigation it was shown that the peripheral part of the blast furnace is characterized by the minimum process of direct reduction and also shown that the uniform distribution of burden load provides the minimum fuel consumption.
Key words： Burden； Blast furnace； Heat and material balance； Modeling； Ring radial zones； Vertical temperature zones； Coke rate； Productivity； Natural gas； Blast temperature
I. G. Tovarovskiy （2013）. Substitution of Coke and Energy Saving in Blast Furnaces.Part 1. characteristics of technology and uneven processescognition， calculation， Forecast. Energy Science and Technology， 6（1）， -0. Available from： URL： http：///index.php/est/article/view/10.3968/j.est.1923847920130601.000
DOI： http：///10.3968/j.est.1923847920130601.000
1. CHARACTERISTIC OF THE BLAST-FURNACE SMELTING
Blast melting is one of the few industrial technologies， including processes and unit， to preserve the essence and significance by all technical revolutions. This phenomenon deserves special consideration in terms of its specific and system properties， providing stability in a dynamic industry environment.
Countercurrent principle of technology， carried out in the closed unit shaft type， allows to ensure maximum utilization of the energy input in the base system， and ease the use of the exported products.
The presence at the bottom of the blast furnace carbon extension provides the unique variant， which is typical only for this technology， a feature of combining in one unit three phase state of charge （solid， liquid and softened）， located in a counter with gas. The mentioned set of features would seem to be technically impossible， if the blast smelting was invented for anew， but not exist in reality. However， the course of blast-furnace smelting in modern aggregates is characterized by a high resistance during long-continuous operation. Last reached in the result of long evolution development of technology of securing the benefits inherent in shaft counterflow. Unique properties of the blast furnace， that supply a steady flow of processes at high efficiency， were formed in the long course of its evolution.
The advance of the blast furnace practice went at constantly growing rate during the XIV-XIX centuries and did not slow down in the XX century. Among the major achievements in the relevant equipment and processes are the following：
Intensification and capacity growth of the blowing means and the related increase of the unit volume of the furnaces while their shape， refractories and equipment were upgraded；
Transition to mineral fuel-coke and continuous improvement of its physico-mechanical properties；
Usage of compound sintered iron-containing raw materials with specified properties comprised of a flux component and other additives；
Blast preheating， its enrichment with oxygen and injection of gaseous， liquid and powdered solid fuel and reducing additives， as well as hot reducing gases；
Use of up-to-date means of control and methods of smelting processes and equipment operation monitoring.
The XX century advance resulted in an increase in the unit capacity up to 10-12 thousand t/day and fuel rate reduction to 450-490 kg/thm while there are only 280-300 kg/thm of coke in the total fuel volume， the rest being powdered coal or other coke replacement materials.
Data showed that by the end of each century the specific productivity （t/m3.d） of the best blast furnaces increased exponentially （Figure 1）， while the coke rate （kg/thm） lowered almost linearly； such a phenomenon is not quite typical for the development of technical systems and is a proof of the high unused reserves of the blast furnace practice. As regards the coke rate， we see a weak， but characteristic for such systems， form of changes expressed by two conjugated exponents dying down at the final portion.
The dying-down curve is distinctly seen for the average indices of the blast furnace practice of the second half of the XX century in the USSR （Figure 2）， and is characteristic for other countries. The integral use of the conventional projects known in the art and those being developed to upgrade the blast furnace operation will allow， during the next decade， to reach a specific production capacity of 4 t/m3.d and coke rate 250 kg/thm （with 200-250 kg/thm pulverized coal injection―PCI） thus approaching the ultimate operation characteristics of a blast furnace within its inherent potential.
K-Coke rate， kg/thm， П-Production capacity， t/m3.d.
2. METHODICAL RECEPTIONS OF ANALYSIS
The design methods of the analysis are the formalized tool of an understanding of processes in the quantitative form basing on the fundamental laws of the nature. The perfection of this tool is determined by a level of knowledge of the fundamental laws and correctness of their use for the analysis of specificity of investigated processes. The base laws for the closed systems are the laws of preservation of weight and energy， which use at the design analysis of blast furnace smelting assumes a solution of two base tasks： 1） study of interrelations of parameters and characteristics of real blast-furnace smelting； 2） forecasts of expected characteristics of blast-furnace smelting on preset parameters of work.
The realization of each of the specified tasks is based on a solution of specific problems. In case of the forecast these problems are connected to a correctness of initial assumptions and expedient degree of detailed elaboration of the description of processes （I. G. Tovarovskiy， 2009； I. Tovarovskiy， 2012）. In case of the analysis of real technology major value have completeness， accuracy， reliability of the initial information and way of its treatment.
Results of analysis of real technology depend to a great extent on reliability of the initial data， on which he bases. However the last， as a rule， have errors， which sizes are unknown. This not only complicates the analysis， but sometimes deforms its results. Use of design complex parameters and balances including sizes gauged with errors of a different nature and size is especially complicated.
Author of present article in a course of analysis of the characteristics of branch blast-furnaces working has shown that the revealing of disparities of iron， slag-forming and gasified elements balances is a major component of the analysis and has developed a technique of definition of equivalent rejections of parameters at entering the corrections into the initial data and appropriate computer system of the analysis （I. G. Tovarovskiy， 2009； I. Tovarovskiy， 2012）. The main role of duration and repeatability accepted for the analysis of the individual periods of work of the furnace is shown， and also method of treatment of results on the ground of the principles of selection and data processing are formulated （I. G. Tovarovskiy， 2009； I. Tovarovskiy， 2012）.
The developed system of analysis includes， alongside with gauged initial parameters of processes， also design smelting characteristics： complex parameters of blast regime； iron， slag-forming and gasified elements balances； a heat balance； reducing-thermal and gas dynamics characteristics. The important component of the analysis is the estimation of influence of errors of the initial data on the design characteristics of processes. Without this component the analysis has not that definiteness， which is necessary for reliable conclusions. So， by not estimating an error of definition of any complex parameter in two compared periods of a blast-furnace operation it is impossible to estimate importance of difference of its values in these periods. For such estimation the balances of iron， slag-forming and gasified elements are used. The disparities of these balances are the integrated characteristics of errors of the initial data. The analysis by computer of ponderability of various errors of the account in total value disparities of balances allows to estimate the most probable variants of errors and close to “true” value of design parameters. At comparison of “true” values of parameters in the various periods the “zone of an error”， allowing is taken into account and to judge importance of difference of “true” values. So， if at difference of “true” values of a parameter in two periods or for two furnaces “the zones of errors” do not coincide or impose against each other in a small measure， it is possible to consider the specified difference essential. If the imposing of “zones of an error” is large， difference is insignificant. In a number of cases it is necessary to ascertain impossibility to make the certain conclusions of the analysis. Such conclusion is not less valuable， than determined， since allows to avoid erroneous judgements.
The further work in this direction has revealed opportunities of perfection of a technique of use disparitys of balances for an estimation of influence of errors of the initial data on “validity” of design smelting characteristics. The considered below new methodical developments differ from known. On a base of a resulted below analytical solution of a task it is possible to replace used earlier iterative procedure of disparities separation by a more simple and evident way of an estimation of influence of possible account errors on smelting characteristics and to give a solution more suitable kind for the substantial analysis. Besides in considered statement a task of overlapping of balances of iron and slag-forming constituents of burden for the first time is decided.
3. FORECASTS OF EXPECTED CHARACTERISTICS BY MULTI-ZONE MODEL OF BLAST-FURNACE PROCESSES
When choosing a method of forecasting performance and processes of blast-furnace author used the previously performed analysis， which showed （I. G. Tovarovskiy， 2009； I. Tovarovskiy； I. G. Tovarovskiy， 1987， p.192） that the adequacy of models of real processes depends mainly on the degree of investigated processes. Because the adjustment to the real conditions on the parameters of the internal state can only be very approximate （rather qualitatively）， its production is carried out on the weekends parameters （coke consumption， performance， parameters of cast iron， slag and furnace gas） that does not allow to give an unequivocal assessment of the adequacy of the model to the real processes. This causes the need to enter empirical coefficients， that are not constant， containing insignificant base of which is not always unambiguous. Despite these difficulties model， which help to better understand the processes and set tasks to further study.
Use of results of experimental researches of blast furnaces， the synthesis of theoretical knowledge about the processes significantly promote the development of a comprehensive model of blast-furnace， the most important results were obtained by Japanese （Blast furnace phenomena and modeling， 1987） and Russian developers （Bolshakov & Tovarovskiy， pp. 207-226）. The results obtained illustrate the possibility of a wide use of models for the analysis of real technologies and development of new technological solutions. To date， however， such a large-scale analysis for any one model has not been conducted. The reasons for this are not the only difficulties rethink of the whole technology as a whole system， but also the fact that it requires a specific building of models for convenience of handling them in the course of analytical research.
Having set himself the task of overcoming these difficulties， the author of the article began by creating its own model for analytical researches of processes of blast-furnace smelting. In the presence of the models created by other experts， the creation of his was due， inter alia， necessity compliance with the requirements of consistency of parametric analysis of the performance and processes， including the adequacy of simultaneous reflection on the possibilities of all the processes and indicators on all parameters. Only when using such a model， it is possible to identify a number of regularities， traditionally falling out of attention of researchers and remaining outside the analysis. Specified by the regularities after checking on real objects served as a basis for the deepening of the conclusions and developing new technological solutions， its compliance with the requirements of consistency of parametric analysis of the performance and processes， include the adequacy of simultaneous reflection on the possibilities of all the processes and indicators on all parameters.
Developed in ISI NASU mathematical model of blast-furnace processes is built on the basis of the structural linkage of multiband height and the radius of the blast furnace and General balance of mass and heat. When modeling the blast furnace smelting， the uneven distribution of materials and gases in 12 vertical temperature zones （VTZ） in height and 10 radial of ring zones （RRZ） on the radius of the blast furnace determines the appropriate uneven flow of the processes and polymorphous temperature-concentration， phase and gas-dynamic fields of the furnace volume.
3.1 Zones
The system of discrete material and thermal balances in the radial annular cross-sections for 12 vertical zones is used to describe the heat and mass transfer processes of iron ore reduction in blast furnace. The temperature drop in vertical zones starts from initial temperature to 400°C for the first zone and further with 100°C intervals for other 11 zones with material descend and transformation from the solid to the liquid phase. The temperature of the melt is the bottom boundary of the last zone. The heat and mass transfer processes for specific temperature intervals in each vertical zone are as follows：
From charging temperature to 400°Cheating and moisture evaporation；
From 400° to 900°Сmaterial heating and indirect reduction of iron， reduction of easy reduced materials and dissociation of carbonates and hydrates；
From 900° to 1100°Сheating and a softening of materials， gaseous reduction of iron， carbonates dissociation， Boudourd reaction of CO2 and carbon and coke de-volatilization；
From 1100° to 1300°Сheating and fusion of materials， formation of slag， iron gaseous reduction， carbonates dissociation， Boudourd reaction and reduction of hard reduced materials；
1300°Сtemperature of melt， reduction of the hard reduced materials and melt superheat.
The gas movement in the tuyere zone of blast furnace is three-dimensional. In the bosh the gas flow becomes two-dimensional and in the bellyone-dimensional （Gordon， Maksimov， & Shvydky， 1989； Spirin & Ovchinnikov， 1995）. The gas velocity arrows are parallel in the zone of one-dimensional movement in the furnace stack. However， the magnitude of the gas velocity and the gas flow in various areas of the furnace cross-section depends on burden permeability. It is assumed that isobars should be perpendicular to the material streamlines， which is confirmed by the results of the numerous experimental investigations. Two- and three-dimensional gas flow in the area of furnace bosh and belly also can be represented as one-dimensional. In this case the isobar linked the center of the furnace cross-section at the tuyeres level （the horizontal velocity in this point equal to zero） and the point at the bosh wall with the same static pressure should be assumed as the bottom boundary of material column. Above this boundary （isobar） the gas flow through each ring zone depends on the permeability of the burden in this zone.
The typical pressure distribution along the height of blast furnace for the normal furnace operation is very well established in many publications （Zhavoronkov， 1944； Gudenau， Kreibich， & Nomia， 1979， pp.7-13； Bonnenkamp， Engel， Fix， & Grebe， 1983， pp.21-26； Gudenau， Kreibich， & Peters， 1981， pp.13-18； Shturman & Gudenau， 1982， pp.9-14； Peters， Pot， & Peters， 1986， pp.10-20）. This quadratic function line could be used to approximate the pressure distribution based on continuously measured blast pressure， top gas pressure and the gas pressure in the middle of the stack. The iteration process allows simultaneous solution of the gas dynamics and heat and mass transfer problems in the volume of blast furnace.
3.2 Heat Distribution
The total volume of the top gas， its temperature， and composition are estimated based on the overall heat and mass balances. The values of the gas flow parameters are assigned for each RRZ based on the previous experience as the first approximation. Knowing the boundary conditions on the furnace top it is possible to estimate heat and material balances for the first vertical temperature zone （VTZ） （tminit 4000C） for each RRZ. Further the heat and material balances are estimated for VTZ from 2 to 12 （with 1000C material temperature interval）. The heat distribution amongst the separate RRZ is proportional to the amount of materials loaded into this zone and amount of the hot metal produced in this zone.
During calculation of the burden physical enthalpy increment， linear dependence of the specific heat on temperature has been assumed with conditional division into the following components （kg/thm）： M―metal-forming； S―slag-forming； Cburden carbon （solid）； O―burden oxygen； CDcarbon dioxide of the burden carbonates； Mo ―burden moisture. At the upper boundary of the first VTZ （charge level） the burden initial enthalpy， E， is （kJ/thm）：
For other zones simple expressions of a similar form for DEiII； DEiIII ； …..DEiXII have been also obtained.
For peripheral RRZ （# 10） the burden enthalpy increment decreases by the value of heat losses through a wall. At known overall magnitude of heat losses Qloss （kJ/m2 of surface） its distribution amongst zones can be assumed based on the linear dependence of losses on the gas temperature.
The heat of reactions is associated with the temperature ranges of their course. Inside these ranges they are distributed at regular intervals between corresponding temperature zones.
The uniform distribution of the heat of iron direct reduction between VTZ # 7-10 can be conditionally accepted only for the first approximation of the iterative process. However， because of significant contribution of the direct reduction into the overall heat balance， this distribution is estimated in the course of following iterations based on the reduction process kinetics.
3.3 Estimation of Iron Direct Reduction
The following kinetic equation （I. G. Tovarovskiy， 2009； I. Tovarovskiy， 2012） is used to estimate the distribution of heat and the volume of oxygen reduced during direct and indirect reduction amongst VTZ from 2 to 10： JR= [（CO*-CO）×VCO×kCO+（H*-H）×VH×kH]×（VFO）ji /[（VR）ji+（VC）ji]， m3/（kg×sec）， where： （CO*- CO）， （H* H） are the differences of concentrations of carbon monoxide and hydrogen-equilibrium （*） and actual， m3/m3， respectively； VCO and VHconsumption of CO and Н2 respectively per unit of reduced material， m3/（kg sec）； kCO， kH―constants of CO and H2 oxidation rate respectively， m3/（kg sec）， determined experimentally for each material with an establishment of their temperature dependences： kCO = fCO （t）； kH = fH （t）； （VFO）ji/[（VR）ji + （VC）ji]―relative volume of iron oxides （VFO）ji in the burden， m3/m3； （VR） and （VC）bulk volumes of iron bearing materials and coke， respectively （m3/thm）.
The amount of oxygen reduced from the burden by gases in each zone is calculated by the following formula： OR = JR×tji×GFe-O， m3/thm， where tji is material residence time in the zone （j-# VTZ， i-# RRZ）； GFe-O is quantity of reduced material （iron oxide）， kg/thm. Based on the above parameters， the gas composition （COСО2； Н2Н2О） is recalculated on the boundary of the following zone.
Gas composition on boundaries of zones and its volume are estimated by adding the oxygen produced by the gaseous reduction to a part of CO and Н2 （with their conversion into СО2 and Н2О） and adding to the total amount of gas the values of СОd and СОL （m3）， formed by oxygen of direct reduction of iron and hard to reduce elements respectively， with relevant gasification of the burden solid carbon Сd and СL （kg）. The calculated balance values， are used for estimation of heat transfer parameters in a specific zone （heat capacity of the burden and gas flows， ratio of heat capacities） and the magnitude of volumetric coefficient of heat transfer av = A× （W/（m3×0K）） （A is a constant）. The residence time of material in a specific zone is calculated based on the following formula： t = q / {av×[（VR）ji + （VC）ji]×（T-t）}， sec. Special analytical investigation was performed to derive the formula for av calculations.
3.4 Analytical Investigation of Heat Transfer Coefficient
Comparing the formulas of Zhavoronkov （1944） for DP： DP = M×w1，8×r0，8×T/（d1，2×e1，8×P）， Pa/m with formula （Kitayev， Yaroshenko， & Lazarev， 1966； Timofeyev， 1949； Heunert & Willems， 1959， pp.1545-1552； Bogdandi， Engel， 1971； Aerov & Todes， 1968） for volumetric heat transfer coefficient av ： av= N×w0，6×r0，6×T0，3/（d0，4×e0，6×P0.3）， W/（m3×K）， it becomes possible to derive a new formula for av： av=A×. Here： M， Nconstants； w， r， T， Р―velocity and density （at normal conditions）， temperature and pressure of gas， respectively； d， ethe linear size of voids （analogue to the hydraulic diameter） and voidage of a layer， respectively. Such unconventional form of the formula for heat transfer coefficient is more universal， than known ones， and flexible for fine-tuning of constants and iterative calculations.
4. Gas-dynamic Characteristics
The bulk volume of the materials （VR）ji + （VC）ji， （m3/thm） in a specific zone， incorporated into the calculation procedure， has been defined with respect to the ironbearing components phase transformation with the temperature increase. The bulk volume of the ironbearing materials is defined by the expression （VR）ji = Rji /（gR）ji， where Rji―consumption， kg/t and （gR）jibulk density， kg/m3.
Density of iron-bearing materials total mass： zones 1-8： gRji； zones 9，10： gR9i = （gm +gRji）/2； gR10i = gm ， whereas melt density gm = （Mji+Sji）/（Mji/7800 + Sji/2，500）， kg/m3； where： Mji/7，800―volume of hot metal； Sji/2，500volume of slag， m3/thm； gRji―apparent density， kg/m3 （3）
Bulk density of iron-bearing materials： zones 1-6 （gR）6i； zones 7-10 （gR）7i = [（gR）6i +gRji]/2； （gR）8i = gRji； （gR）9i = （gm +gRji）/2； （ gR）10i = gm. The density of coke total mass in all zones gCji and coke bulk density in all zones （gC）ji are assumed as identical. Based on equations （3） and （4） it is possible to estimate the volume of materials total mass and bulk volume of material in various zones. The linear size of lumps of coke fCji and iron-ore raw materials fRji are estimated by empirical functions of temperature.
The gas velocity under normal conditions could be calculated from the equation of the pressure drop in a packed bed with respect to above specified parameters： w = [DP×e3×P / （M×r0，8×`A1，2×T）]0，555， m3/m2×sec. The gas volume is used in calculation of the heat balance of each zone. Also it is used for estimation of the material residence time in the zone： t = q / {av×[（VR）ji + （VC）ji]×（T-t）}， sec and subsequent estimation of the height and the cross-section area of each zone.
Material flow （VV） and linear velocity （VL） should be determined prior to estimation of the geometry of each zone： VV = [（VR）ji + （VC）ji]×P， m3/sec； VL = VV / AH， m/sec， where： Pblast furnace productivity， thm/sec.
Now it is possible to estimate the height of the VTZ―“hji” for each RRZ： hji = tji×VLji （m） and the distance of top and bottom borders of each VTZ from the stock level as the summation of their heights： Нji = hji + hji-1 + hji-2 +…+ hj1.
4.1 Boundary Conditions for Material Distribution
The boundary conditions for material distribution on the upper level of material column （upper boundary of RRZIi） are determined as follows：
The iron bearing burden mass RI?i and its fraction DI1， DI2， …DI10 （DI1 + DI2 + … +DI10 = 1） are set for each RRZIi.
Initial relative burden distribution （iron bearing materials mass/coke mass） are also assigned for each RRZIi：
RI1relat. = DI1/dI1； RI2relat. = DI2/dI2 … ….. RI10relat. = DI10/dI10， and RI1relat.* RI2relat.* ….. *RI10relat. = 1，0.
The coke mass in each RRZ is calculated based on coke fraction d and RIirelat.： ?dIi = 1.
Such assignment of the raw materials distribution at the furnace top allows to describe all possible variants of materials distribution by coke-less top， two-bell system， rotary charge system （Totem） etc.
4.2 Algorithm Brief Description
The multi-zone balance includes the system of expressions and formulas connected by feed-forward and feedback links. The solution of the multizone balance is coordinated with solution of the overall balance under the predictable mode of operation. The predictable mode requires defining of the parameters of the reference period. The reference period allows estimate parameters of operation for each RRZ and VTZ at existing conditions of the blast furnace. By virtue of this， an iterative algorithm with multistage iterative cycles is used to solve the problem.
4.3 Analysis of the Results of Calculations
The average annual operating parameters of blast furnace （BF） with the useful volume of 5，000 m3 have been used to illustrate the capability of the developed method. The specific production of BF in this period of time was 2 thm/m3/day （working volume） and the coke rate was 412 kg/thm. BF operated with the following actual distribution of raw materials on the furnace top：
The parabolic and uniform types of material distribution also were used for estimation of BF operating parameters.
The parabolic material distribution could be attributed to the two-bell charging system. The uniform material distribution on the furnace top characterizes the bell-less top charging system. The overall coke rate decrease in the case of bell-less top compare to the two-bell system is about 3-4%. Improvement in the process stability allows the decrease in a “reserve of technological heat” in the furnace hearth and as a result of this reduction in hot metal silicon and coke rate. Because of this it is reasonable to expect 4-5% reduction in overall coke consumption. The extend of coke rate decrease was experimentally confirmed by results of operation of blast furnace #6 （with bell-less top） and #5 （two-bell system） at Novolipetsk I & S Works （Bolshakov， Tovarovski， & Shutilev， 2005）.
5. Influence of Burden Materials Distribution on the Blast-furnace Performance
The improvement in a blast furnace performance as a result of optimization of metallic burden distribution was the main goal of the current study. All parameters were estimated with respect to the operation of Blast furnace #5 of Severstal and Blast furnace #9 of Mittal Steel Krivoi Rog. The operation of these furnaces was analyzed in three different modes： retrospective， current and perspective. The real mode corresponds with current furnace operation. There is no natural gas and oxygen injection and blast temperature is decreased to 700 °С for the retrospective mode. The perspective mode assumes PCI and coke oven gas injection with the blast temperature of 1200 °С.
Blast furnace performance at each mode of operation was estimated for two distributions of metallic burden load （MBL）―actual and conditionally-uniform （further-“uniform”）. The “uniform” distribution can be achieved at different MBL on the blast furnace periphery， therefore， the optimization was done based on the minimum total fuel consumption. Other distributions， for example， parabolic， appears to be less effective and characterize mainly two-bell charging system. Analysis of actual furnaces performance confirms the last statement （Bolshakov， Tovarovski， & Shutilev， 2005）.
Results of blast furnace operation for different modes are presented in Tables 1 and 2 for blast furnaces #5 of Severstal and BF #9 of Mittal Steel， Ukraine.
The following conclusions could be made from the results of blast furnace performance investigation：
The uniform distribution of MBL provides the minimum fuel consumption
It is necessary to increase MBL on the periphery of the furnace with decrease in coke rate to achieve the minimum coke consumption.
It is also necessary to increase the average MBL value using such technological measures as PCI， injection of coke oven gas， increase in blast temperature etc.
The major factors in the reduction of the coke rate with uniform distribution of MBL are as follows：
Improvement in heat transfer which finally leads to decrease in temperature of the top gas
Reduction in heat demand for direct reduction due to decrease its fraction and increase in a degree of gas utilization
Reduction in heat losses through the walls because of increase of MBL at periphery of the blast furnace and an increase in the furnace productivity
The improvement in the uniformity of the radial temperature and gas distribution with the transition to the uniform MBL distribution reduces the number of RRZ with degenerative upper heat exchange zone. It is well known that the degenerative heat exchange increases the total fuel rate because of the increase in heat demand.
The results of investigation showed the direction in improvement of MBL distribution： uniform MBL distribution in RRZ 2-9 at maximum possible MBL at periphery. It is also important to maintain the axial gas flow， which mainly depends on quality of coke （especially coke charged into the axial part of the furnace cross-section）， quality of metallic burden， parameters of blast and other technological factors.
The comparison of Severstal and Mittal Steel， Ukraine blast furnaces performance with three types of mode of operation （retrospective， current and perspective） shows the following：
The extend in the coke rate decrease during the transition from the current mode of operation to the perspective one for Severstal blast furnace is lower120，7 kg/thm compare to 148 kg/thm for Mittal Steel， Ukraine blast furnace
The extend in the coke rate increase during the transition from the current mode of operation to the retrospective one is higher for Severstal blast furnace―133 kg/thm versus 122.6 kg/thm for Mittal Steel， Ukraine blast furnace.
Such difference reflects the fact that operation of Severstal blast furnace is much closer to the limiting conditions compare to operation of Mittal Steel， Ukraine blast furnace.
ACKNOWLEDGEMENTS
With special thanks to my colleagues from Iron and Steel Institute of Ukrainian National Academy of Science and professor Gordon I.M. from Hatch， Canada， also to authors given in the references.
REFERENCES
Aerov， M. E.， & Todes， O. M. （1968）. Hydraulic and thermal fundamentals of operation of apparatus with stationary and fluidized granular bed. Leningrad. （in Russian）.
Blast Furnace phenomena and Modelling. （1987）. London and New York： Elsevier applied Science.
Bogdandi， L. A.， & Engel， G. Y. （1971）. Reduction of iron ores （Translation from German）. Moscow， Metallurgiya.
Bolshakov， V. I.， & Tovarovskiy， I. G （Eds.）. （2006）. Knowledge of the processes of blast-furnace smelting. Dnipropetrovsk， Porogy. （in Russian）.
Bolshakov， V. I.， Tovarovskiy， I.， & Shutilev， F. M. （2005）. Effect of charging devices application for blast furnaces. Stal， （7）， 17-20. （in Russian）.
Bonnenkamp， H.， Engel， K.， Fix， V.， & Grebe， K. （1983）. Investigations of the frozen by nitrogen blast furnace. Chernye Metally， （2）， 21-26.
Gordon Y. M.， Maksimov， E. V.， & Shvydky， V. S. （1989）. Mechanics of material and gas movement in shaft furnaces. Alma-Ata. Nauka Kazakh SSR. （in Russian）.
Gudenau， G. V.， Kreibich， K. & Peters， K. H. （1981）. Optimization of blast furnace fusion zone profile. Chernye Metally， （3）， 13-18.
Gudenau， G.V.， Kreibich， K.， & Nomia， Ye. （1979）. Impact of fusion zone on gas flow distribution in blast furnace. Chernye Metally， （22）， 7-13.
Heunert， G.， & Willems， J. （1959）. Stahl and Eisen ， （79）.
Kitayev， B. I.， Yaroshenko， Y. G.， & Lazarev， B. L. （1966）. Heat exchange in blast furnace. Moscow， Metallurgiya. （in Russian）.
Peters， K. H.， Pot， G.， & Peters， M. （1986）. Investigation of the reasons of gas permeability failure in blast furnace. Chernye Metally， （22）， 10-20.
Shturman， E.， & Gudenau， G. V. （1982）. Investigation of blast furnace fusion zone at Shvelgern I & S Works. Chernye Metally， （6/7）， 9-14.
Spirin， N. A.， & Ovchinnikov， Y. N. （1995）. Heat exchange and blast furnace efficiency improvement. Ekaterinburg， UGTUUPI. （in Russian）.
Timofeyev， V. N. （1949）. Heat exchange in the packed bed. VTI News， （2）. （in Russian）.
Tovarovskiy， I. （2012）. Volume 1， Status analysis. Volume 2， Problems and Prospects. Processes of blast furnace smelting （Vol.1， p.595； Vol.2， p.406 ）. LAP-LAMBERT Academic Publishing. （in Russian）.
Tovarovskiy， I. G. （1987）. Blast furnace process parameters. Upgrading and optimization. Moscow， Metallurgiya. （in Russian）.
Tovarovskiy， I. G. （2009）. Domennaya plavka. Monography （2nd ed.）. Dnepropetrovsk， Porogi. （in Russian）.
Zhavoronkov， N. M. （1944）. Hydraulic basics of processes and heat transfer in scrubbers. Moscow， Sovetskaya Nauka. （in Russian）.