книги / Типовые задачи оперативного управления непрерывным производством
..pdfПереходя к |
переменным |
х(к\ |
вместо формулировки |
(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-1 — Gik),
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)