Управление формированием рудной базы металлургического предприятия на основе интеграции моделей доменного процесса и разрушения агломерата

Journal of Siberian Federal University. Engineering & Technologies, 2018, 11(2), 168-180

yflK 669.162.1

Managing Iron and Steel Works Ore Base Formation Through Integration of the Blast-Furnace Process and Agglomerate Failure Models

Mikhail Yu. Ryabchikov, Vera V. Grebennikova and Elena S. Ryabchikova*

Magnitogorsk State Technical University by G.I.Nosov 38 Lenin, Magnitogorsk, 455000, Russia

Received 05.06.2016, received in revised form 09.10.2016, accepted 24.12.2017

The paper addresses the issue of optimizing ore feedstock base to operate iron and steel works. It is shown that fluctuating ore feedstock base assessed on the basis of drum cold mechanical strength index TI can result in falling short in achieving the optimal economic values of agglomeration and blast-furnace process. A method of combined use of the stress model of solid inclusion in matrix and the failure model based on the Rittinger's law is suggested to account for effects of varying feedstock agglomerates on agglomeration and on efficiency of the blast-furnace process. The metallurgic agglomerate failure model also accounts for combined use of the stress model developed for solid inclusions in matrix. The algorithmic issues dealing with complex models integration are addressed. The calculation results are provided to estimate the change in the blast furnace production efficiency when agglomerates made of various ore materials are used. The comparison of agglomerates having different strength expressed through TI index is made.

Keywords: ore feedstock base, quality of metallurgical agglomerate, management, blast-furnace process, production efficiency, coke consumption.

Citation: Ryabchikov M.Yu, Grebennikova V.V., Ryabchikova E.S. Managing iron and steel works ore base formation through integration of the blast-furnace process and agglomerate failure models, J. Sib. Fed. Univ. Eng. technol., 2018, 11(2), 168-180. DOI: 10.17516/1999-494X-0020.

© Siberian Federal University. All rights reserved

Corresponding author E-mail address: mr_mgn@mail.ru, greb1981@mail.ru, mika.elena@mail.ru

Управление формированием рудной базы металлургического предприятия на основе интеграции моделей доменного процесса и разрушения агломерата

М.Ю. Рябчиков, В.В. Гребенникова, Е.С. Рябчикова

Магнитогорский государственный технический университет им. Г.И. Носова Россия, 455000, Магнитогорск, пр. Ленина, 38

В статье рассмотрены проблемы, связанные с формированием рудной базы металлургического предприятия. Показано, что при непостоянстве рудной сырьевой базы предприятия оценка полезности агломерируемого сырья с позиции барабанного индекса холодной механической прочности П может не обеспечить достижение оптимальных экономических показателей аглодоменного производства. Для учета индивидуальных особенностей влияния агломератов из различного сырья на производительность доменного процесса предложен способ интеграция модели доменного процесса с моделью разрушения металлургического агломерата, основанной на комбинированном использовании модели напряженного состояния твердого включения в матрице и модели разрушения на основе закона Риттингера. Рассмотрены алгоритмические сложности организации интеграции и поиска решения. Приведены результаты расчетов изменения производительности доменной печи при разрушении агломератов из различных рудных материалов, обладающих различной прочностью при испытании, на барабанный индекс прочности П.

Ключевые слова: рудная сырьевая база предприятия, качество металлургического агломерата, доменный процесс, расход кокса.

Введение

В настоящий момент при сравнении полезности рудных материалов, используемых в доменном процессе после их агломерации, как правило, оперируют усредненными эффектами влияния качественных характеристик агломерата на среднюю производительность и расход кокса в доменном процессе [1]. Основными показателями оценки качества агломерата являются индексы: механической холодной прочности TI (Tumble Index согласно ISO 3271); вос-становимости RI (Reductive Index согласно ISO 4695); прочности после восстановления RDI (Grounding strength in temperature recovery по ISO 4696 - 1).

На рис. 1 приведены результаты расчетов [2] влияния восстановимости и прочности железорудных материалов на технико-экономические показатели доменной печи №9 «ММК» относительно базового варианта работы печи в 2009 г. Представленные зависимости являются усредненными и не предполагают учета рудной принадлежности материалов, используемых при агломерации.

В работе [3] предложена модель разрушения металлургического агломерата, основанная на комбинированном использовании модели напряженного состояния твердого включения в матрице и модели разрушения на основе закона Риттингера. Определены нагрузки, действующие на куски агломерата в испытательном барабане, и предложен способ адаптации модели

455

450

445

440

435

430

425

420

415

1

/ *2

2 / ^ 1

/ у

/ 1 ,2

/ к 1 1

2

60

70

80

90

4 700

4 650

4 600

4 550

4 500

4 450

4 400

4 350

100

Восстановимость (1); холодная механическая прочность (2),

Рис. 1. Рассчитанноеа[2]влияниевосст ановимости(1)ихолоднойпроч ности(2)н арасходкокс а (—) и производительности—(длопсиовийнименной печи NiBOAO«MMK »

Fig. 1. Effect of the reducibility (1) and cold crushing (2) strength, calculated in [2], on the coke consumption (—) and production rate (—) for conditions of blast furnace No. 9 of OJSC Magnitogorsk Iron & Steel Works

разрушенияагломерата с использованием данных о гр анулометрических характеристиках продуктов испытанию При адаптацил прадло асано испонлзовать вер ояпностные распределения параматзювмоьали разрушгнтя.Полученепзеоисинасзиэтих пасомвелзя озосдавности агломерата,оодержанияуглелода вшихтеи кррчностиитаастнокк-Важлым илеольтатом работы подели npai вавном значении инрлоса И1 лля гелсмератов аз еазник матеаи алое.

При условии стабильной рудной базы использование усредненных эффектов влияния на эффективность доменного процесса оправданно. В то же время в текущих экономических условиях взаимоотношения с поставщиками рудных материалов могут отличаться нестабильностью [4].

Например, теперь аглофабрики «ММК» получают усредненную смесь из рудных материалов более чем ста различных месторождений и, исходя из ее параметров, выполняют дозирование топлива и известняка, что является основными управляющими воздействиями на качество агломерата. Для уточнения индивидуальных эффектов влияния механической прочности агломератов из различных железорудных материалов на эффективность доменного процесса при формировании рудной базы предприятия необходима интеграция модели разрушения агломерата [3] с моделью доменного процесса.

Интеграция модели разрушения агломерата в модель доменного процесса

В доменной печи можно выделить две зоны - косвенного восстановления железа и прямого восстановления. В первой (верхней) зоне большое влияние на перепад давления газов оказывает гранулометрический состав шихты и механическая прочность агломерата. Во второй (нижней) зоне важна горячая прочность компонентов шихты, которая характеризует способ- 170 -

1800 1600 " 1400 j£ 1200 ce 1000

g 800 I" 600

1 400

H

200 0

Рис. 2. Принятое допущеине лросвределении темпераеур повыслте ДП:Г=- 0,0019H5 + 0,1085H4 -- 2,94h3 + 41,9h2 - 299,9H + 1778,°C

Fig. 2. Accepted assumption concerning the temperature distribution along the blast furnace height: t =-0,0019H5 + 0,1085H4 - 2,94h3 + 41,9h2 - 299,9H + 1778, °С

ность матееиалов-епротивляться нтгрузвопрв -ысокихтемпернтури-. шения

агломерота [С] ориеииирсдана на8елхт°ююе-ео=. Поотоиуротлмоту88ный даяотдасчет учитывает грнннлометдическдйнотяде в сорхнейч-тти ннчи>-дит нижнтИчасви ит-нет ведется согласномотодлке,предложеннойв]5].

Модель разрушения [3] определяет статистические вероятностные распределения геометрических и прочностных параметров кусков агломерата, представленных заданным количеством (20 тычн), раорушение кото-ых пв мере ввижснеоо ie шахте доменной печи требуется про-моделив-вать.

Длявьшоинениящшближвнныхравчятовнфинярофиквированное распределение температур повысоте доменной печи (ДП) №9 ОАО ММК, изгбраженное на рис. 2, согласующееся с данными [6, 7]. Принятая высота нижней зоны определяется высотой распара и заплечиков и составляет Нрасп + Нзапл = 5 м, а высота верхней зоны Нш = 15 м при общей рабочей высоте 20 м.

Согласно [8] приняты следующие газодинамические параметры базового режима: избыточное дявлеиие ^тья^вкм4 избыточное давлснке ия колошнике 1,в атм.; верхний перепад давления 0,4ятм.;нижнийбебатач даввения К,9 ктм.;обищИпербптддаяления 1,2 атм.; расход гомччего двтья (>0 ]в]Ис; темпенатдра торячетд дд^тьв ЮИ0 °С. Согласно [д] зивиснм-енъ сопро-нквлтдия елла шиояыы ф нижнеИ ч-сти домскихИ птчн м вяетом псромекров шихты и дутья для базового сдрианда домянноо плтвкл выражается уреонением i^-ода

А Р

Аяз р А^ - ИИ,

где Лниз - коэффициент потери давления в зоне размягчения, Па-с2/м6; АРниз - нижний перепад давления, Па; F„ - объемный расход горячего дутья, м3/о.

Коэффициент Лниз использовали для упрощенного учета влияния нижней зоны на измене-ниеперепададавления приизменениирасхода дутья.

Для описания газодинамических процессов в верхней зоне принято известное уравнение С^гош^Е^^^]:

Высота доменной печи, м

w2 AH (1-ош) T 0

АР = м—2-0 2 ш0-t-г,

^ еЩШ 273 (l + P)' тя

где AP - перепад давления в слое исследуемого элемента, Па; ц - коэффициент сопротивления слоя шихты; иа- сторюсть газа черлз поюй шааты, мГо; рп ^ плеоэость газа (воздухаО, кг/м3; и!МУ- E^:t>icoTa шихты, ч;иЛкв.ш _ эквиталентныолиаметр кусков шихты, м; еш - по-

Лто:з:ност:е> inK-^i^i, отн. (^д.; Г-темпэратура газа в исследуемом слое, °С; P - давление газа в ис-следуемомслое,Па.

Исходя из (2) и ряда работ, где используюоея схежие выроженит [Ы Юслезет аависи-месть перемадн давленипот порозносои, эквивгланиног— диаметр куоков шихты - коэффициента сопрогав/етия.

Эквивалененьш диаметр млжет биооь риссчитан ия ситоопои состава менердалов, по кото-Томр няреднляеодя стаанлилшдиаметр аиавиц [10/, л аохваносто - иэ нлношноп и ктжущейсо плотностиматериала[13]:

0 y 4 m

= 1-J0o = 1--7TT7n--е (и

у a гынж яОДЯ Yrr ti d

где й?экв н- элвивалененый миаметр,мм;о,- - оПлнм/ап доля фрй0фии,пвпавшей Bi-ый интервал, отн. ед.;гЭ -середвна/'-оо гнггрвало диамстра— к^лзв, м—; ан, укаа л-насыпнвниксни^лася плотность материала, кг/мс (вкаж = 3,(4 доя агломерата); 4/(тт • 2ЛПОЯП -л воличина;еТкиалал о,,) саму исследуемогя слоя мктериалов в тгечи, Ум3; т - массаматериала вобъемеслоя, кг. Порос-ослоаглтмсрае-и ияонл ре^дт(^дтс1^.в[л:и пп зависимлсоям /11, 1Н, 14]:

0;403 • а034с = 0,403 •

еа d .

0,04

^ = 0,000 • агк0,

где еа - порозность агломерата, отн. ед.; ^экв.а - эквивалентный диаметр кусков агломерата в объеме исследуемогоэлемента,мм; ек — порозность кокса, отн. «од.; ¿¡Л,. _ эквиваоененый дие-метркусков коксавобъемеисследуемого сегмента, мм.

Очевидно, что при моделировании разрушения маоериаоов вЦЩ еиоКоодиоюучшыоать разрушениевсехкомпонентов.еклюиавнеотъемлемыш иомхогент - кокс, объемная доля кото-ряио в шихые знеиимооьна. В ту еке ндемя модель [3] ориентирована на разрушение агломерата. Н а основании выражения

М°=М°0+0Л44Ь^ь5тМ-10,

где А/к - исходная массовая доля мелочи в коксе, %; М-10, М+25 - показатели прочности кокса по ГОСТ 5953-93, %; а3 - коэффициент учет влияния мелочи, можно найти среднее значение мелочи коксаМк>рпри (М_10= 8,5 % [15, 1(5]; -0%° = 3 %; И,ол = 20Э4 м3) по еысоте верхней зоны доменной печи. Притимая длиущение о зинойкюм величендк доли мелочи кокса по

вы(5те шахты, определим долю мелочи в нижней части зоны косвенного восстановления (6), что позволяет рассчитать долю мелочи в любом слое по высоте:

Ме.нон = 2 • Ме.нр -Ме«рх>

где Мк.низ - массовая доля мелочи кокса в нижней части зоны косвенного восстановления, %; Мк.верх - исходная массовая доля мелочи - 5мм в коксе, %.

Для расчета эквивалентного диаметра кусков кокса в зависимости от доли мелочи выполнили анализ экспериментальных данных работы [15] по распределению крупности кусков кокса на разных ступенях разрушения в барабане Сундгрена (рис. 3а), что позволило получить зависимость среднего диаметра кусковкокс а отдолимелоч и-5мм(рис. Зб).Параметры пороз-ности иэкв ивалентного диаметраокатышей приняты на постоянном уровне.

В доменной печи патски шиоты и газт находятся р протиаетоке и их яраречное движение возможно толтко и]е> в соДчювении определхнных а слоги й, иднио из которьгн явлгется ровный сход шихзовых матсрягооз, гар^ктиг1ззузо^вйс^опр^де1^^1^но11 стзпеныз уравновиг^ивания (СУ) шпхеы газогым пороком. В работе [5] для расчета СУ предасжеигвофзжение

СУ = '3 = КГЛР-.V = К, Ал>р = кур

' ^ннуК ИМа^0Я'

гзе /еа-ио1сосл(акткопогиВскояшихтя1 отврквни воздк-псылфурм ян увввняшихты в печи, м; А/. -настяая нраоаниягытоспнсыыете аныноний части, атм.; - насыпная масса шихты, т/м3; ^-площодьпонхрочнтхо хачозик печ^ г-Т-пхромотр адаптации(%10).

Н аличие информации о размерах кусков материалов в процессе их разрушения позволяет перейти от плoщaдиоoпeвганoеoceчeдмх посхв п О^Ха^ рффхктихнoйплoмади(хкмeних)мyсков 50 в овпвхм гк:гыы(^:^тыг но пытогт аиронтП магистяе плчи:

М у2

£ ^е^нрог.аеыглы,.0

од! Эз]

где SjJ,Rjа,Nlа - эффективная площадь, радиус и число кусков 1-й фракцииу-го компонента шихты; k2 - число интервалов размеров кускови число компонентов шихты.

0,08

0.07

ГГ

s X 0.06

с

=t 0,05

CL

о CD 0,04

JO

о 0.03

X

о с 0,02

l_

0.01

0

а)

у А

У\ А —

-Л

20 40 60 80 Средний размер кусков кокса, мм

70

60

50

40

30

20

О 10

100

б)

у = 60.3С ,8е-0.0438х

10 15 Доля мелочи -5 мм. %

20

25

Рис. 3. Распределение крупности кусков кокса на разных ступенях разрушения в барабане Сундгрена

Fig. 3. а) Distribution of the coke lump size at different fracture stages in a Sundgren drum (1, 2, 3 - beginning, middle and completion test stages); б) calculated dependence of the mean coke lump diameter on the share of smalls -5 mm

При выполнении настройки модели на базовый режим параметр КГ упрощенно приняли одинаковым для каждого сегмента по высоте, фракции и компонентам шихты.

Алгоритм расчетов модели газодинамических процессов доменной печи в верхней зоне представлен на рис. 4. На рис. 5 изображена структурная схема взаимодействия параметров блоков модели.

Сложностью используемой схемы расчетов является то, что для первоначального расчета перепадов давления и степеней уравновешивания нет данных о разрушении кусков агломерата и порозности шихты по высоте доменной печи. Поэтому перед началом расчета принимается фиксированная линейная зависимость изменения порозности шихты и эквивалентного диаметра кусков агломерата и кокса по высоте верхней зоны доменной печи.

Расчет перепада давления происходит снизу вверх по слоям. Рассчитываются текущая плотность газового потока рва8, его объемный расход Fvol.dut и динамическая вязкость д. По формуле (2) вычисляется перепад давления АР в текущем слое высотой ЛН=1м. Для следующего сегмента на высоте давление внизу слоя, с учетом потерь в нижележащем, равно Р = Р - АР, далее производится вычисление перепада давлений с шагом расчета АН до верха доменной печи.

Затем рассчитывают степень уравновешивания шихты СУ восходящим потоком газа. При этом слои перебираются напротив сверху вниз для определения силы Qн давления на материалы текущего слоя верхними слоями, расчета разрушения кусков и определения эквивалентного диаметра кусков агломерата dэкв.а и кокса. Находится новая порозность агломерата еа, насыпная плотность агломерата рнаса для каждого слоя и пересчитываются объемные доли компонентов шихты фа, фк в исследуемом сегменте объемом У8е8.

Расчет повторяется циклически (5-10 циклов) до получения установившегося решения по перепадам давления АР по высоте доменной печи и порозности шихты еш.

На рис. 6 показано изменение рассчитанных по модели параметров по высоте верхней части доменной печи №9 ОАО «ММК» для шихты из смеси Соколовского агломерата (49 % по массе), кокса (28 %) и окатышей (23 %) при индексе прочности агломерата Т1 = 70. Общее среднее значение степени уравновешивания шихты (с учетом СУ нижней зоны) составило 0,56 при коэффициенте КГ = 0,08.

Разработанная модель позволяет оценить влияние принадлежности агломерационного сырья к различным месторождениям на СУ доменного процесса при одинаковом индексе холодной механической прочности. В таблице 1 приведены данные об изменении: СУ; расхода дутья, обеспечивающего СУ = 0,57, стандартное отклонение которого составило 0,96 м3/с, или 1,57 %; производительности и расхода кокса согласно эффектам в [1].

Из данных табл. 1 следует значимое влияние рудной принадлежности агломерационного сырья на производительность доменной печи.

На основе зависимостей настроек модели разрушения агломерата от индекса Т1 в [3] можно оценить влияние на эффективность доменного процесса изменения механической прочности агломерата в достаточно широком диапазоне. На рис. 7а видно, что влияние индекса Т1 на производительность доменной печи неоднозначно и может быть как линейным, так и значимо отличным от линейной зависимости.

На рис. 7б результаты расчета по изменению производительности доменной печи представлены в форме зависимостей от содержания углерода в агломерационном процессе, полученных

Рис. 4. Алгоритмрасчетовмоделигазодинамических процессов доменнойпечив верхнейзоне Fig. 4. Algorithm of model calculations of gas-dynamicprocesses in the upper zone of the blast furnace

Рис. 5. Структурная схема взаимодействия параметров блоков модели при расчетах Fig. 5. Structural scheme of the interaction of parameters of model units in calculations

Рис. 6. Изменение рассчитанных по модели параметров по высоте верхней части доменной печи №9 ОАО «ММК» для шихты из смеси Соколовского агломерата

Fig. 6. Change in parameters calculated as per the model along the height of the upper part of blast furnace No. 9 of OJSC Magnitogorsk Iron & Steel Works for the furnace charge made of mixture of Sokolovskoye agglomerate

Таблица 1.Рассчитанноевлияние руднойпринадлежностиагломерационногосырья на пок азатели доменного процесса при TI = 70

Table 1. Calculated effect of ore belonging to the agglomeration of feedstock on the values of the blast-furnace process at TI = 70

№ Принадлежность рудных материалов Основность СУ ^ (м3/с) при СУ = 0,57 ДП", % ДК, %

1 Оленегорский концентрат 0,3 0,563 60,40 -0,93 0,044

2 0,75 0,557 60,ПТ -0,412 0,020

3 1,2 0,558 60,ПТ -0,412 0,020

4 Высоаогорская аглоруда 0,1 0,524 62,60 1,89 - 0,091

5 1,2 0,528 62,П1 1,64 - 0,079

6 Соколовская аглоруда 0,1 0,565 60,27 -1,09 0,052

7 1,2 0,560 60,ПТ - 0,(577 0,032

* Од, ДП, ДК- расходдопьт, изменение производительности печи.

30 50 70 90 1 3 5 7

Индекс TI Содержание углерода. %

Рис. 7. Рассчитанные зависимости снижения производительности доменной печи от индекса TI (а) и содержанаяугиерюда ори агломерчзы а Об)нз оыыошеии юпрежыму знпользования тгломерата из Высокогорской руды при TI = 79 для агломератов из: 1 - Высокогорской руды, В = 1,2; 2 - Высокогорской руды, В = с,1;3д Соколовскойзуды,В м 1,и;° -Соколовсао^удь^В ы0Д

Fig. 7. Calculated dependences of the blast furnace production rate decrease on TI index (a) and carbon content in agglomeration (b) in relation to regime of the 33elomerateu se from Vysokogorskoye ore at TI = 79 for agglomerates of: 1 - Vysokogorskoye ore, B = 1.2; 2 - Vysokogorskoye ore, B = 0.1; 3 - Sokolovskoye ore, B = 1.2; 4 - Sokolovskoye ore, B = 0.1

с использованием моделей [16, 17]. Из формы зависимостей видно, что при нестабильной сырьевой базе поддержание на постоянном уровне содержания углерода топлива в агломерационной шихте может приводить к значимым колебаниям производительности доменной печи, что снижает стабильность доменного процесса.

Заключение

Основным результатом интеграции моделей доменного процесса и модели разрушения агломерата является возможность решения задач формирования рудной базы предприятия с

учетом индивидуальной эффективности использования агломератов из различного сырья в доменном процессе, что актуально в условиях непостоянной сырьевой базы.

Из полученных результатов следует относительно слабое влияние непостоянства газодинамических режимов доменной печи на удельный расход кокса в доменном процессе, что объясняется низким значением используемого усредненного коэффициента из работы И.Г. Това-ровского [1], предполагающего изменение удельного расхода кокса 0,2 % на 10 кПа изменения общего перепада давления. Однако нужно учитывать, что рост механической прочности может сопровождаться в отдельных случаях ростом горячей прочности и изменением восстановимо-сти. Учет данных явлений наряду с вероятностными характеристиками распределений показателей качества агломерата позволяет в перспективе повысить эффективность систем управления формированием рудной базы металлургического предприятия и планирования режимов работы доменных печей.

Список литературы

[1] Товаровский И.Г. Нормативная оценка влияния параметров доменной плавки на расход кокса и производительность. Сталь, 2014, 5, 4-11 [Tovarovskii I.G. Regulatory impact assessment parameters on blast furnace coke consumption and productivity. Steel, 2014, 5, 4-11 (in Russian)]

[2] Дмитриев А.Н., Витькина Г.Ю., Чесноков Ю.А. Создание методических основ анализа влияния качественных характеристик железорудного сырья и кокса на технико-экономические показатели доменного процесса. Физическая химия и технология в металлургии. Сборник трудов, посвященный 60-летию ИМЕТ УрО РАН. Институт металлургии Уральского отделения Российской академии наук, Екатеринбург, 2015, 309-314 [Dmitriev A.N., Vit'kina G.Yu., Chesnokov Yu.A. Creation of methodical bases of the analysis of the influence of qualitative characteristics of iron ore raw materials and coke on technical and economic indicators of the blast furnace process. Proceedings: Physical chemistry and technology in metallurgy. The manifold published for the 60th anniversary of Metallurgy Institute, Ural Branch of the Russian Academy of Sciences, Ekaterinburg, 2015, 309-314 (in Russian)]

[3] Рябчиков М.Ю., Гребенникова В.В., Рябчикова Е.С., Богданов Н.В. Модель разрушения металлургического агломерата. Известия ВУЗов. Черная металлургия, 2016, 3, 159-166 [Ryabchikov M.Yu., Grebennikova V.V., Ryabchikova E.S., Bogdanov N.V. Model of the metallurgical agglomerate failure. News of Higher Educational Institutions. Ferrous metallurgy, 2016, 3, 159-166 (in Russian)]

[4] Волков Ю.В., Соколов И.В., Смирнов А.А. Стратегия освоения сырьевых ресурсов Урала. Горная промышленность, 2006, 4, 57-62 [Volkov Yu.V., Sokolov I.V., Smirnov A.A. Strategy for Ural feedstock supply. Mining industry, 2006, 4, 57-62 (in Russian)]

[5] Спирин Н.А., Лавров В.В., Рыболовлев В.Ю., Гилева Л.Ю., Краснобаев А.В., Швыд-кий В.С., Онорин О.П., Щипанов К.А., Бурыкин А.А. Математическое моделирование металлургических процессов в АСУ ТП. Екатеринбург: ООО «УИПЦ», 2014. 558 с. [Spirin N.A., Lavrov V.V., Rybolovlev V.Yu., Gileva L.Yu., Krasnobaev A.V., Shvydkii V.S., Onorin O.P., Shchipanov K.A., Burykin A.A. Mathematic modeling of metallurgic processes in Automatic Process Control System. Ekaterinburg, OOO «UIPTs», 2014. 558 p. (in Russian)]

[6] Большаков В.И. Познание процессов доменной плавки. Днепропетровск: Пороги, 2006. 439 с. [Bol'shakov V.I. Study of blast furnace smelting processes. Dnepropetrovsk, Porogi, 2006. p. (in Russian)]

[7] Паспорт доменной печи №9 ОАО «ММК». 1999. 28 с. [MDS for blast furnace No. 9 of OJSC Magnitogorsk Iron and Steel Works. 1999. 28 p. (in Russian)]

[8] Спирин Н.А. Компьютерные методы моделирования доменного процесса. Екатеринбург: УГТУ-УПИ, 2005. 301 с. [Spirin N.A. Computer methods for blast furnace smelting process modeling. Ekaterinburg, UGTU-UPI, 2005. 301 p. (in Russian)]

[9] Мухин Р.А. Разработка теоретических основ поведения шихтовых материалов доменной плавки для информационного обеспечения математической модели загрузки доменной печи. Москва: НИТУ МИСиС, 2010. 28 с. [Mukhin R.A. Developing theoretical basis to optimize blast furnace smelting blended materials to provide the information for the mathematical model of the blastfurnace charging. Moskva, NITU MISiS, 2010. 28 p. (in Russian)]

[10] Большаков В.И., Муравьева И.Г., Семенов Ю.С. Применение радиолокационных систем измерения поверхности засыпи шихты для контроля и управления доменной плавкой. Днепропетровск: Пороги, 2013. 364 с. [Bolshakov V.I., Murav'eva I.G., Semenov Yu.S. Application of radar systems to analyze blended mixture stock-line for blast furnace smelting for control and management purpose. Dnepropetrovsk, Porogi, 2013. 364 p. (in Russian)]

[11] Ковшов В.Н., Петренко В.А. Экспериментальные исследования движения шихты и газа в доменной печи. Днепропетровск: Институт технологии, 1996. 124 с. [Kovshov V.N. , Petrenko V.A. Experimental study of the blended mixture and gas motion in blast furnace. Dnepropetrovsk, Institut tekhnologii, 1996. 124 p. (in Russian)]

[12] Тарасов В.П. Газодинамика доменного процесса. М.: Металлургия, 1982. 224 с. [Tarasov V.P. Blast-furnace process gas dynamics. Moskva: Metallurgiya, 1982. 224 p. (in Russian)]

[13] Ковшов В.Н., Галаганов А.И., Чистяков В.Т. Исследование порозности многокомпонентной шихты, Металлургия и коксохимия. Днепропетровск: Техника, 1997, вып. 53, 34-37 [Kovshov V.N., Galaganov A.I., Chistyakov V.T. Study of the multi-component blended mixture void volume. Proceedings: Cast iron metallurgy and coke chemistry. Dnepropetrovsk: Tekhnika 1997, 53, 34-37 (in Russian)]

[14] Kajiwara J. Development of Simulation Model for Burden Distribution of Blast Furnace Top, Tr. ISS Jap, 1983,V. 23, 12,1045-1052.

[15] Иванов Е.Б., Мучник Д.А. Технология производства кокса. Киев: Вища школа, 1976. 232 с. [Ivanov E.B., Muchnik D.A. Coke production technology. Kiev, Vishcha shkola, 1976. 232 p. (in Russian)]

[16] Ryabchikov M.Y., Grebennikova V.V. Simulation of the combined effect of production factors on metallurgical sinter mechanical strength. Metallurgist. 2013, 57, 3-4, 274-283

[17] Гребенникова В. В., Рябчиков М. Ю., Бурдина О. В. Комплексный критерий оценки качества металлургического агломерата. Автоматизированные технологии и производства. 2013, 5, 163-169 [Grebennikova V. V., Ryabchikov M. Yu., Burdina O. V. Metallurgical agglomerate quality assessment complex criterion. Automated technology and production, 2013, 5, 163-169 (in Russian)]