Файл: Основы автоматизации для металлургов..pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 14.10.2024

Просмотров: 83

Скачиваний: 0

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

12. Тепло, поглощаемое свободной известью

т

f ^ C a O _

/ j ^ C a o

\ G z

(

dSl)

dt

\ A s l

I 100

V dt

j

 

 

 

 

 

 

M C t o \

G,

I -

l

n

dS

 

2Ap

J 100

V

dt J y

As

 

I 100

 

—0,44 -

 

 

\

 

 

dt

После подстановки получим

 

 

 

 

 

 

 

 

 

д я ( 1

2 ) = [(»o 25) 0,2245 + 338,8017]

 

 

_ 0

0 3

9 9 7

О г

^

 

0.0362G,

dP

\

Л

/

Л ,„

rfS

 

 

 

 

 

 

J".

)

0,017489G,

—0,44

dt

 

 

 

 

 

 

 

Ä

 

z \

 

 

 

 

 

 

 

13. Потери тепла через футеровку:

 

 

 

 

 

 

 

 

[см. уравнение (217)].

 

 

 

 

 

 

 

 

 

 

 

 

14. Тепло, излучаемое из центра реакции

в ванну:

 

 

 

Д Я ( 1

4 ) = *Sts

+ 273)4 _ ( » +

273)4]

 

 

 

 

 

 

 

 

[см. уравнение (216)].

15.Тепло, теряемое с продуктами сгорания:

(15)sp sp Psp

Ч , = ^ с о а 9 5 + ^соЛ 0 5 -

После преобразования и вычисления получим

Д Я ( 1 5 ) = - 0 , 0 0 5 2 6 0 3 О Л / ) ( - - g - ) .

(232)

( 2 3 3 )

(234)

Таким образом, при процессе обезуглероживания прирост тепла за время dQ равняется сумме

 

 

 

15

 

dQ = Д Я ( І ) dt + Д Я ( 2 ) dt+

. . . + Д Я ( 1

5 ) Л =

2 А и і d t -

(235>

 

 

 

î

 

Путем интегрирования

уравнения

(235)

по / от начального

времени 4 до ко­

нечного времени / к получим уравнение для QK Qo, т. е. уравнение для прироста тепла Д<2 всей плавки.

В этом уравнении нужно путем дальнейших преобразований выразить все

связи в зависимости

от количества чугуна

(так как он подвергается

процессу

обезуглероживания)

или от количества FeO,

переходящего

в шлак. Дл я

расчета

требуемого расхода

охлаждающего скрапа (охлаждающих

присадок) нужно в ка­

честве зависимой переменной определить требуемое количество скрапа Gs во всей системе. Так как общее количество металлической шихты определяется требуемой массой плавки Gt, а носителями этой шихты являются чугун и скрап (окалину в этом случае нет необходимости учитывать, так как ее доля в металлической шихте составляет всего 1,5%), справедливо выражение

G = G - + G ,

сz s

откуда

Gz = Gc Gs

где Ce общая масса металлической шихты.

299


Подставив выражение, полученное для

Gz, в предыдущие формулы, решим

всю систему относительно Gs как зависимой

переменной.

 

Однако, поскольку между количествами

чугуна и скрапа существует комута-

тивная связь только по материальному балансу

(с точки

зрения металлической

шихты), а не по энергетическому балансу,

нужно

выразить

изменения массы ком­

понентов (в данном случае прежде всего

Gs)

в функции энергетических разностей

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

Приведенные связи следует дополнить несколькими условиями:

а)

принимают, что углерод, кремний и фосфор содержатся

в жидком чугуне

в виде РезС, FeSi и FeeP;

 

б)

на основании литературных данных [34] предполагают,

что при обезугле­

роживании чугуна кислородом образуется максимально 5% СО2 (при условии, что процесс находится под постоянным контролем) :

в) в уравнениях принята 50%-ная десульфурация. Из этого количества 12% серы окисляется в элементарную серу, которая удаляется с газами, а 88% серы реагирует с СаО и CaS;

г) количество FeO в шлаке (Greo) выражено в виде функции массы скрапа, массы окалины, расстояния от фурмы до зеркала ванны и основности. На основе

экспериментальных результатов

[45], можно получить конечную форму

выра­

жения:

 

 

 

 

 

<3р е 0 =

3175,19 — (34019,9 — GS) (—354 685 • \0~WT + 0,162162)

+

 

+

635,05 ( ß — 2,8) + 0,5 ( G 0

— 635,05),

(236)

где Gs

— масса скрапа;

 

 

 

 

Ht

— расстояние от фурмы до зеркала ванны;

 

 

В — основность;

 

 

 

 

G0

— масса окалины;

 

 

 

д) при расчете тепла, израсходованного на растворение извести, исходят из

реакции CaOs —>СаО(е)

(растворенной в шлаке).

 

 

Хотя температура шлака при кислородно-конвертерном процессе не достигает

температуры плавления

СаО, следует предположить, что для

растворения

СаО

в шлаке необходимо такое же количество тепла, какое требуется для расплав­ ления;

е) обычно температуру чугуна измеряют в миксере. Снижение этой темпе­ ратуры к моменту заливки чугуна в конвертер может быть определено из выра­ жения [45]:

9- = 0,958-т + 37,778,

где т ) т — температура металла в миксере; ж) при расчете принимают, что все реакции протекают только в чугуне, а

скрап добавляют в качестве охлаждающей присадки, которая не участвует в ре­

акции. Поэтому количество чугуна в конце плавки

равно:

 

 

GZK

=

GZQ

-

i^L [(С0 -

С*) +

( S i 0 -

Sis) + (P0

-

Р й ) +

( M n 0 - Mnf t ) +

0 , 5 S 0 ] -

 

 

 

' F e O

 

O F E O

- G V

+

^

 

?

G

 

0

,

 

 

 

 

(237)

 

 

 

 

"

"

 

 

™ F e , 0 ,

 

 

 

 

 

 

 

где

p — содержание

Рез04

в

окалине

(чистота

окалины

колеблется

в

пределах

 

 

 

0,8—1,0);

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

з) количество извести, необходимой для достижения заданной

основности,

можно

выразить

 

 

 

 

 

 

 

 

 

 

 

 

 

 

G c a 0

=

і

^

.

. J _ в

S i o - S l »

 

 

 

 

 

 

 

 

( 238)

^ а о

 

 

A s l

 

<j>

100

 

 

 

z

 

 

 

 

 

 

где

ф — доля

чистой СаО в извести;

 

 

 

 

 

 

 

В — основность,

выраженная

 

как соотношение

СаО/БіОг-

 

 

300


Конечное уравнение для определения количества охлаждающего скрапа (охлаждающих присадок) или температуры выпуска выводится из уравнения (235) с учетом условий, приведенных выше, и здесь не приводится в силу своей громоздкости.

Надежность

физико-химической

модели

 

 

За критерий надежности модели в мировой литературе принимается доля

плавок (в процентах), температура выпуска которых лежит

в пределах

до­

пусков ( ± 1 5 град, от

заданного

значения). Модель,

имеющая

надежность

75%,

считается хорошей, а 80% очень хорошей.

 

 

 

Были проведены

проверочные испытания математической

модели, которая

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

Первые

результаты

оказались

неудовлетворительными. Только 43,2%

плавок

имели температуру,

соответствующую фактической с точностью ± 1 5

град.

Стан­

дартное

отклонение

о* = 30 град.

Гистограмма распределения разности

температур

приведена на рис. 204.

 

 

 

 

25

 

 

 

 

 

го

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

а

20

 

 

 

 

 

15

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

15

 

 

і

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

10

 

 

 

 

 

g

10

 

 

 

 

 

 

 

 

 

 

 

 

-30 -20

-10

0

+10 +20

+30

 

-20

-10

О

+10

+20

 

 

Разность между замеренным

и

 

Разность

меікду замеренной и

 

расчетным

конечным

значением, град

 

расчетной

температуройб,град

 

 

 

 

 

 

 

 

Рис. 204. Гистограмма

р а с п р е д е л е н и я

\xh

Рис.

205. Гистограмма

распределения

fih

 

при

# . „ = 1550° С

 

для

выбранной

совокупности

при

0 в р

=

 

 

 

 

 

 

 

 

= 1550° С

 

 

 

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

 

Таким образом из нескольких кампаний были

выбраны только эти плавки

и

они были вновь промоделированы. Результаты подтвердили ожидания. В пре­

делах

допусков ± 1 5 град

находилось

76% плавок,

стандартное отклонение г> =

=

8,7

град. Гистограмма

распределения

показана на

рис. 205.

Этот факт подтвердил известный принцип, что автоматическая система упра­ вления работает тем точнее, чем более точную информацию она получает, т. е. чем точнее и надежнее работают датчики. Ясное представление о влиянии этого фактора было получено при определении влияния ошибок измерения на конечную температуру стали, как это видно из табл. 34.

Испытания на модели позволили выявить некоторые зависимости, например, необходимую долю скрапа в шихте (%) в зависимости от содержания угле­ рода в чугуне и зависимость этой доли от температуры чугуна. Графическое изо­ бражение этих зависимостей приведено на рис. 206 и 207.

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

301


Т А Б Л И Ц А 34

Влияние ошибок измерения на конечную температуру

Параметр

Масса плавки, кг . . .

Масса чугуна, кг . . . .

Масса скрапа, кг . . .

Масса руды, кг . . . .

Температура чугуна, °С

Чистота извести, % . . Состав чугуна, % '•

Со

Si0 Мпо

Ро So

Состав стали, %: С*

Sis Мп*

Р а

Sa

Доля углерода, окисля­ ющегося до СО, % • • •

Расстояние от фурмы до зеркала ванны, мм . Выбросы, %

Точность

Влияние

ошибки на

измерений

температуру стали,

 

град

 

 

+ 3 0 0

 

+ 2 , 8

+ 3 5 0

 

+ 0 , 8

+ 2 0 0

 

+

1,5

+ 2 0 0

 

+ 4 , 4

+ 5

 

+ 3 , 0

+ 0 , 2

 

+

1,8

± 2

 

+ 2 , 8

+ 0 , 0 5

 

+ 3 , 4

+ 0 , 0 4

 

+ 4 , 1

+ 0 , 0 5

 

+ 3 , 0

+ 0 , 0 0 6

 

+

1,3

± 0 , 0 2

 

+ 0 , 3

+ 0 , 0 3

 

+ 2 , 1

+ 0 , 0 1

 

+ 2 , 5

+ 0 , 0 5

 

+ 3 , 0

+ 0 , 0 0 6

 

+

1,3

+ 0 , 0 2

 

+ 0 , 2

± 1

 

+ 9 , 5

+ 100

 

+

1,8

± 1

 

± 2 , 0

1250

1*А = 1580 °С

 

С0

= 6,23%

 

Sig=

0,65%

1200

10°C=OJ% скрапа

 

 

 

 

 

 

 

 

 

 

 

 

1150

 

 

 

 

 

 

 

 

16

ig

га

 

 

гг

 

 

1100

 

 

 

 

 

 

 

 

 

 

 

 

15

17

19

21

23

25

26

 

 

Доля

скрапа

в

шихте,%

 

 

 

 

Доля

скрапа

б шихте, %

 

Рис.

206.

Зависимость

доли

скрапа

в

ших­

Рис.

207. Зависимость

доли

скрапа

в

ших­

те от

с о д е р ж а н и я

углерода

в чугуне

(# 2 =

 

те

от температуры

чугуна

 

 

 

 

= 1200° С; f>0 =1580° С)

 

 

 

 

 

 

 

 

 

 

 

что

она

является

в

некоторой

степени

общей,

соответствующей

любой

плавке,

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

3 0 2


риальных и энергетических констант, а также максимальной точности измерений. Так как наши знания и возможности в этом смысле являются довольно огра­ ниченными, работа с этой моделью имеет часто очень низкий процент на­ дежности.

2. СТАТИСТИЧЕСКАЯ М О Д Е Л Ь

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

Если параметры процесса изменяются лишь в небольшом диапазоне, можно

сформулировать довольно точное линейное выражение

м і

у я » у* = 2

сцхеі +

Ьха + с,

 

 

 

 

(239)

 

;' =

і

 

 

 

 

 

 

где

у— фактическое значение входной величины

(например,

массы

скрапа);

хеа

у* — расчетное значение этой

величины;

 

 

 

 

— соответственно, входные

и выходные

параметры процесса;

 

а,Ь,

с — постоянные коэффициенты модели.

 

 

 

 

Д л я определения коэффициентов

а, Ь, с используем метод наименьших квад­

ратов таким образом, чтобы сумма

квадратичных отклонений

между

измерен­

ными и расчетными значениями величины у была

минимальной:

 

 

N

 

 

 

 

 

 

 

 

2 ( У ; ~ У - ) 2 = т і п .

 

 

 

 

 

(240)

і = і

 

 

 

 

 

 

 

 

Хотя этот метод и требует минимальных вычислительных мощностей, он не

всегда

дает

хорошие

результаты. В

последнее

время

проводятся опытные ра­

боты, направленные на минимизацию абсолютных значений ошибок. Каким ме­

тодом в конце концов будут определены коэффициенты, будет зависеть

от

ем­

кости

запоминающего устройства

ЭВМ, надежности сбора

измеренных

данных

и других

факторов.

 

 

 

 

Число

N измерений величин,

с которыми модель оперирует, зависит

как

от требуемой статистической надежности определения коэффициентов (с

увели­

чением

N

при одинаковой вероятности ошибки надежность возрастает),

так и

от самого

процесса. Например, если длительность кампании

ограничена

сроком

службы футеровки, то в качестве максимального

числа

измерений можно выбрать

число тех случаев, которые можно охватить в рамках одной

кампании.

 

Уравнение

(239) является очень

простым

выражением

модели

и,

если изве­

стны значения

входных и выходных

величин

е,

х„),

то,

решив

это

уравнение,

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

устройства

(например,

аналоговой

ЭВМ)

Знание

этих коэффициентов

на практике является, однако, препятствием

для широкого применения такой

модели, поскольку в цехе модель должна слу­

жить для

определения

величины

регулирующего воздействия у еще до начала

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

соответствовать конкретному

процессу, так как модель

не

может

приспособиться

к

измененным условиям

хода

процесса. Точно так же

ошибка в

измерении или

в

передаче информации

при

определении коэффициентов

вызовет

их искажение

303