Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги / Типовые задачи оперативного управления непрерывным производством

..pdf
Скачиваний:
6
Добавлен:
12.11.2023
Размер:
19.07 Mб
Скачать

Переходя к

переменным

х(к\

вместо формулировки

(4-6) получаем:

х<АХ х {к),

 

.... п

О<

1,

 

 

 

(4-7)

А-(«+>) —щах.

Преобразованная задача (4-7) эквивалентна задаче (4-6) и является задачей линейного программирования.

Задача минимизации издероюек. Если в задачу (4-6)

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

yi'l)=zx(n+')>yW,

(4-8)

гдe_j/t">— заданное число, а в качестве функц'ионала при­

нять издержки производства вида

а

= 2

&k)x(k) —-min.

(4-9)

=1

переработки сырья

на /г-й опе­

где с№)— себестоимость

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

Рассмотрим ту же задачу для комплекса последова­ тельно соединенных простых операций вида (4-2). Пред­ положим, что издержки производства иа каждой опе­ рации можно записать как

2<ft)= g/t(6(ft))x<4

(4-10)

где gu — заданная функция от 0<Ч

П

Тогда задача минимизации издержек Fz =

£ z<k} при

А=1

ограничении на план выпуска (4-8) примет следующий вид:

л '1’ — f, (d<4) л:'*» =

0;

 

*<!>— М 0<!)) •*<’>= 0;

 

 

•*<"> — /„(0<”>) *(«+■) = 0;

g i (^(1)) х 1-2*

g t (0(S>)

“b • • • Ч- gn (0tn)) JC("+ min

причем

_

*(«+>) ^<">;

*<*><*(*>,

 

 

(4-11)

ёЛ0,4))=гЛ О |‘,)//,(0(‘)). * = i , .... n\

Ô этой задаче переменные коэффициенты £h и /л, за* висящие от к-то параметра 0(А), входят только в один (k+ 1)-й столбец.

В [20] предложен способ решения подобных задач пу­

тем кусочно-линейной

аппроксимации

 

функций

gk (0(/г)),

f (Wk)) в R точках ôJA), г=

1,

...,R:

 

 

 

 

г*(6<*>).*(‘+,>=

2

ï(9}w) ^

+n;

 

 

 

 

 

Г=1

 

 

 

 

 

 

/* ( о<|'))^+ ,=2 M O -'f '” ;

 

(4-12)

 

 

Г=1

 

 

 

 

 

 

2

л (А+1)= ^ (й+,).

 

 

 

 

г = 1

 

 

 

 

 

 

 

 

 

Подставив (4-12) в (4-11)

и

обозначив

(0^)) = /W-

â<P<w)= g kr, получим вместо

(4-11) следующую

прибли­

женную задачу:

 

 

 

 

 

 

 

 

 

R

x {k\

k — 2,

п;

 

 

 

2 x rk)<{

 

 

 

г= 1

 

 

 

 

 

 

 

 

 

2

 

х!"+'>>?<">;

 

 

 

 

* ' > - 2

^

 

 

=

 

0;

 

(4-13)

Г=1

 

 

 

 

 

 

 

 

2 4 е - 2 ^

3) = 0;

 

 

Г=1

 

Г5=1

 

 

 

 

 

Л

 

t

R

 

(«+ п _

 

 

 

~(л)

2

U

0;

 

 

2 ^

-

=

 

 

 

Г=1

 

 

Г=1

 

 

 

 

 

 

2 2

/г=1 г= I

с дополнительными условнямй

л.{ft) Лк)

о,

х<*>=

0, /• =

1,

q— 1, q-\-2, .... к.

VH

 

 

 

 

(4-14)

 

 

 

 

 

Согласно

[20],

если

//*(0(,{))— строго монотонные

функции, a £/i(0(/i)) —строго выпуклые, то для оптималь­ ного решения задачи (4-13) условия (4-14) всегда вы­ полняются. В этом случае решение задачи (4-13) явля­ ется приближенным решением задачи (4-11).

Рис.

4-2. Схема замещения про­

Рис. 4-3. Параллельное соеди­

стой

операции при кусочно-ли­

нение простых операций.

нейной аппроксимации ее ха­

 

рактеристики.

 

Значения переменных

0(ft> определяются из ре­

шений (4-13) следующим образом:

*(*+■ >=2

А = 1, .... п;

Г = 1

 

Q(* > = 2 r=l

Кусочно-линейная аппроксимация зависимостей

M 0(,l))> êit(Qw) имеет простой физический смысл. Каж­ дая операция с этими характеристиками заменяется па­

раллельным соединением (рис. 4-2) линейных операций вида

{к)

^ = Ы к)>/ '= 1,

R.

л /t)__ ftrSr

Новые условные операции можно рассматривать как некоторые варианты проведения исходной аппроксими­ руемой операции.

Коэффициенты Х{гк) являются интенсивностями исполь­

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

б) Параллельное соединение простых операций

Рассмотрим параллельное соединение m простых операций (рис. 4-3).

Задачи распределения нагрузок для таких схем в случае достаточно общих моделей рассмотрены в [7]. Здесь будут рассмотрены лишь модели с переменными коэффициентами видов (4-1), (4-23 и соответствующие упрощенные методы решения.

Задача максимизации производительности. При огра-

т

ничениях на суммарный входной поток 2 ;t(é,= 3 t0, где А=1

х° — заданное число, задача максимизации общей про­ изводительности по выходу может быть сформулирована следующим образом:

т

d<*> <((*>;

И-15)

т

2 -тах

)

и формально является нелинейной.

Если воспользоваться заменой переменных [см. (4-1)],

то задачу (4-15) можно свести к следующей эквива­ лентной задаче линейного программирования:

О < JC(*X Зс<*> ;

т

2 х {к) =~л:0;

k =i

(4-16)

m

2 — max. k=i )

В результате решения задачи (4-16) находятся зна­ чения переменных х^\ yW, значения 0W находятся из

(4-1).

Однако задачу (4-16) можно решить значительно проще. Допустим, что для отдельных операций соблю­ даются следующие неравенства:

еа^е^^ОСЗ)^; ...

(4-17)

Перенумеруем операции с тем, чтобы соблюдалось (4-17). Так как с учетом (4-1) функционал задачи мож­ но записать как

т

f _ у ***>

CJ Ык) * Л=1

то порядок решения задачи состоит в следующем. Задается полная нагрузка первой операции, т. е.

х(,)= .JO), если ir(1)>;c0, то решением задачи является х№=х^=... = х ж=0, если х(1)<х°, то первой

операции устанавливается нагрузка х(1), а остаточный поток х°—а?(|) относится ко второй операции и т. д. Рас­

смотренная процедура полностью совпадает с процеду­ рой распределения нагрузок при постоянных коэффи­ циентах 0=0 [7].

Задача минимизации издероюек. Рассмотрим парал­

лельное соединение простых операций вида (4-12). Пред­ положим, что производительность комплекса по входу и выходу ограничена, т. е.

т

т

2,y'k)=~y;

S

А=1

/г=1

У5

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

0 < fh(0<*>)у<%>,<«*<‘>;

Ш

k = \

т

_

_

(4-18)

2

у& = у,

 

/г= 1...... т;

л=1

 

 

 

т

ÎA(0(fe))^ (è,^m in,

 

 

2

 

 

ft=i

 

 

 

где i* (»'*>) = в» (в<‘>)Ь(в«*»).

Задача (4-18) представляет собой задачу линейного программирования с переменными коэффициентами, причем коэффициенты, зависящие от /г-го параметра, входят в столбец к . Задача (4-18) может быть решена с помощью кусочно-линейной аппроксимации функций fb(0W) и gk('0(ft)) (см. § 4-1,а). Рассмотрим теперь прак­ тические примеры построения моделей комплексов и их решения.

4-2. Координация нагрузок основных участков аммиачного завода

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

ную пропускную

способность,

определяемую числом

работающих

в данный момент агрегатов и их характеристиками.

завода

[71].

Рассмотрим

структурную

схему

аммиачного

В цехе разделения воздуха 1 ’(рис. 4-4)

методом глубокого

охлаж­

дения из воздуха получают относительно чистые кислород и азот. Эти компоненты вместе с природным газом поступают в цех кон­ версии 2, где при высокой температуре в присутствии катализа­ тора природный газ превращается в сырье для производства аммиа­ ка, которое далее направляется для очистки от углекислого газа в параллельные цехи. В цехе 3 продукты конверсии предварительно сжимаются и подаются в цех 6 для растворения углекислого газа водой. В цехе 5 производится отмывка моноэтаноламиновым раст­ вором, а затем сжатие в цехе 4. Сжатый газ (продукт конверсии) подается в цех медно-аммиачной очистки 7, где освобождается от

угарного газа

и остатков

углекислого газа,

откуда

направляется

в цех синтеза

аммиака 8

для производства

готового

продукта —

жидкого аммиака.

Для ритмичной работы завода необходимо согласование на­ грузок цехов, а внутри цехов — нагрузок агрегатов. Эти функции выполняются автоматизированной системой управления «Каскад», разработанной для аммиачного производства [71].

В системе решаются следующие задачи:

определение цеха, являющегося «узким местом» среди парал­ лельных и последовательных цехов;

сравнение ожидаемой и фактической нагрузок производства и выдача рекомендации по их выравниванию;

расчет оптимальной нагрузки с учетом резерва и выдача реко­ мендации по включению резервных агрегатов или остановке неко­

торых работающих агрегатов.

 

 

 

 

 

 

Исходной информацией для системы являются:

 

/'-го

агрегата

G>ih — максимально

возможная

производительность

k-\\ группы t-го

цеха; bUh — признак

состояния агрегата

(ôifA= i ,

когда агрегат

работает, и 0*,А= 0,

когда

агрегат

остановлен);

a,h — коэффициент

приведения

производительности

к

нормирован­

ной нагрузке (к расходу газообразных продуктов конверсии).

Предполагается, что каждая /е-я

группа

t'-ro

цеха состоит из

«.л параллельно работающих агрегатов.

 

 

 

 

 

Первоначально определяется узкое место и производительность каждого из цехов. Для ожидаемой максимально возможной про­ изводительности k-й группы агрегатов /-го цеха можно записать:

nlk

0 * « = 2 Giik8iik

(4-19)

/=1

или в нормированном виде

Gih=aihG*iii.

В результате сравнения максимальной возможной нагрузки Gц.

fe-й группы с наименьшей нагрузкой G)-1 из максимально возможных

нагрузок предыдущих в этом цехе —il) групп получим:

V’ik — Sg (G f-1Gik),

I 1, если ¥=<0;

где Sg ï = ^ 0 _ е с л и ? > 0

После /г-го сравнения для использования в (Л-|-1 )-м сравнении остается нагрузка Gfti, которая определяется по формуле

Gki = M-ifeG} 1 + (1 — p.jft) Glk.

Полученное на «-м шаге сравнения значение Gni является нор­ мированной нагрузкой узкого места данного цеха, ограничивающей его ожидаемую производительность.

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

завода.

В схеме завода (рис. 4-4) цехи 3—6 соединены в параллельные цепочки. Поэтому предварительно определяются ожидаемые на-

7—9 0 }

97

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

Разность ожидаемой максимально возможной производитель­ ности узкого места и фактической нагрузки производства служит мерой для выдачи рекомендаций по повышению нагрузки или раз­ грузки производства.

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

Рис. 4-4. Структурная схема аммиачного завода.

/ — цех

разделения воздуха;

2 — цех

конверсии; 3 — цех первой компрессии;

4 — цех

второй

компрессии;

5 — цех

моноэтаноламиновой

очистки;

б — цех

водной

очистки;

7 — цех медно-аммиачной очистки; 8 — цех

синтеза

аммиака.

вышеописанным, однако в формуле (4-18) учитывается производи­ тельность как работающих агрегатов, так и находящихся в ре­ зерве.

В результате оператору выдаются рекомендации по пуску ре­ зервных агрегатов или остановке некоторых работающих.

Рассмотренные алгоритмы в системе «Каскад» начинают дей­ ствовать при включении или отключении агрегата или при изме­ нении нагрузки более чем на 2% по сравнению с предыдущим за­ мером.

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

4-3. Смесительная операция

Смесительная операция характеризуется нескольки­ ми входными и одним выходным материальным потоком (см. рис. 2-5). Примеры задач оперативного смешения описаны в [24, 82, 83, 90]. Будем рассматривать сле­ дующую модель смесительной операции:

где индекс k соответствует номеру смесительной опера­ ции; х(к), ук—г’-й входной и общий выходной потоки опе­ рации; и[к) i-я составляющая вектора управления и<Л),

связанного с перераспределением входных потоков.

Предполагается, что

где U/t — множество

допустимых рецептов смешения.

 

Рассмотрим случай, когда множество задано анали­ тически с помощью системы равенств и неравенств вида

2

и[к)=

1 ;

(4-21)

ит <ит <Ът,

(4-22)

i

i

i

 

где и[к), и[к)— заданные вещественные числа.

Будем учитывать также качественные показателя ма­ териальных потоков. Зависимость качественных показате-

леи смеси от и\ , играющих роль управляющих воздей­

ствий, примем линейной

 

 

 

 

V

w _ пк

С

“Г-

* = > ......«*■

(4-23)

 

 

 

 

=

2

 

,(*)

 

 

 

/=1

 

 

 

 

 

— значение

s-ro

качественного показателя выход-

где V5

 

кого

потока;

сг( Ь )' — значение s-го качественного пока-

зателя i~го

 

st

 

 

 

 

 

 

входного потока; тк— число контролируемых

качественных показателей.

 

 

 

Предполагается,

что v(k} являются заданными

величи­

нами,

а на

значение

ем5 1

наложены ограничения

вида

ir '

 

 

 

 

 

 

oli) < ü w

V( к )

(4-24)

где

 

 

 

 

 

— S

 

S

или с учетом (4-23)

v(k\ v w — заданные

числа,

—5

 

] ( k )

 

 

 

 

 

 

 

 

 

S

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

vwum < -(»

(4-25)

 

 

 

 

 

 

 

/=1 si i

s

 

Таким образом, множество U/t определяется системой линейных ограничений (4-21), (4-22), (4-25).

а) Общая задача оперативного управлений смеше­ нием

Требуется-найти значения переменных и у ^ к), удов­

летворяющих кроме ограничений (4-21), (4-22), (4-25) дополнительным ограничениям на запасы исходных ком­ понентов

 

и{к)у ^ < х к{

 

 

(4-26)

и плановым ограничениям

 

 

 

 

где х ^к)—заданный

yW^yW,

 

 

 

(4-27)

максимальный запас

t-го

входного

потока; у№— план выпуска продукта смешения.

решений,

Если существует

множество

допустимых

удовлетворяющих ограничениям

(4-21), (4-22),

(4-25),

(4-27), то возможна

постановка

оптимальной

задачи

оперативного управления • смешением.

Эта задача за­

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

мизирующего (минимизирующего) значения

критерия

F(и, у). В качестве критерия F рассмотрим прибыль

Ч

(4-28)

F=cfy<i)—2 c f u J V .

te.I

где c{k)—стоимость единицы t-ro входного потока; с(к) цена единицы смеси.

Сформулированная задача смешения с критерием (4-28) представляет собой следующую задачу математи­ ческого программирования:

(4-29)

Соседние файлы в папке книги