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

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

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

По физическому смыслу не все переменные потоки х(Л), х(°), у(°> и режимные параметры u(ft) являются независимыми. Во-первых, некоторые входные х<°) или выходные у(°) потоки комплекса фиксированы в соответ­ ствии с технологическим содержанием задачи или пла­ новыми ограничениями. Во-вторых, уравнения матери­ ального баланса (1-2) в некоторых случаях позволяют исключить зависимые переменные и тем самым упро­ стить всю задачу. В частности, для этого можно исполь­

зовать простейшие балансные уравнения вида

(1-3).

В общем случае все ограничения,’ на переменные

определяются соотношениями (1-1) — (1-4).

Однако

в каждой конкретной задаче можно исключить некото­ рые переменные и тем самым упростить задачу. Далее будем предполагать, что такие преобразования выпол­ нены и ограничения (1-3) исключены.

При формулировке задач (1-1) — (1-5) не принима­ лись во внимание ограничения, связанные с наличием входных, промежуточных и выходных конечных запасов продуктов в емкостях (на складах). Эти ограничения могут быть учтены. Так, ограничения, связанные с вход­ ными и выходными запасами комплекса, могут быть включены в множества Х<°) и Y(°), для учета промежуточ­ ных запасов можно каждую из имеющихся промежуточ­ ных емкостей представить как самостоятельный; 5-й агрегат с суммарным за интервал управления входом

x(s,/и выходом у я . Тогда формулы (1-2) сохраняют свой вид, а к системе ограничений при наличии рассматри­ ваемой 5 емкости добавляется ограничение на запасы вида

 

О < х (а)' - у

(8)' + -№ < !•'.

( 1-6)

где

I{QS) — исходный запас;

I*s— верхний возможный пре­

дел

запаса.

 

 

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

у<Л)

(1-7)

пых

 

где

— вектор-столбец показателей качества выход­

ных потоков

£-го

агрегата,

V(B*JX= { Ü^ )xs}, 5 = 1 , .... lk\

Лк)

 

 

 

 

 

входных потоков k-vo

VI';.' — вектор показателей качества

гвх

 

 

 

 

 

tk;

 

известные

агрегата,

= { 0^ } ,

г = 1 ..........

g(ft)

вектор-функции, g{k)=

{g'sk)},

s =

l,

....

lk-

 

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

потоков не изменяется

при переходах между агрегатами.

 

 

 

некоторым

Значения

vBB|x, vBX'

должны принадлежать

 

Г(А )

„ < * >

 

 

 

 

 

 

 

ВЫХ*

V BX

 

 

 

 

 

 

заданным ограниченным множествам значений

 

 

чшШ г -

y(«ï

v(k) ç -

у(Л)

 

( 1-8)

 

 

пых

Vnbix’

вх ^

пх

 

Введение условий (1-6) — (1-8)

сужает множество до­

пустимых управлений

D.

 

 

 

 

 

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

кислород,

метан и др.), входящие в поток. При этом

в одном

потоке1 содержатся различные компоненты,

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

Задача (1-1)— (1-8) значительно упрощается, если рассматривать модель (1-1) —(1-3), (1-7) как двухуров­ невую (рис. 1-3). На верхнем уровне / формулируется модель относительно только переменных x<h\ v(ft) не­ зависимо от переменных и<4 При этом происходит рас-

йределение материальных потоков между агрегатами. Определение переменных u<ft) производится на нижнем уровне //, исходя из условия реализации полученных на верхнем уровне значений х^, y(ft\ v<4 Подобный подход возможен в случае, когда переменные u<fe) не входят или несущественно влияют на критерий (1-5). При этом области; значений переменных на входах и выходах [см. (1-4), (1-8)] такие, что для каждого значения пе-

хму(и]0т

Г*)

S

ТТ“Of)

Рис. 1-3. Структура двухуровневой модели комплекса.

/ — модель верхнего уровня; II — то же нижнего.

Рис. 1-4. Схема комплекса с агрегатами, имеющими один вход и выход.

ременной, например х<й), в Допустимой области при лю­ бых допустимых значениях других рассматриваемых пе­ ременных на входах и выходах найдутся такие значения режимных переменных, которые позволят получить это значение х<4

В ряде работ [21, 22] вместо общей модели (1-1) — (1-3) рассматривается более частная модель с агрега­ тами, имеющими один вход и выход, вида

 

$<*)=/№>(*(*>, и<*>),

 

где k — номер

агрегата;

y(h\

Ф \ х№— скаляры.

Выходной

поток каждого

агрегата

распределяется

между входами многих

агрегатов (рис.

1-4), так что

 

Х & =

2

ак1у {1',

 

 

 

fe* *

 

где ctfcj — доля

потока с выхода /-го агрегата, поступаю­

щая на вход

Æ-ro агрегата;

Rk — множество агрегатов,

с которых потоки поступают на вход k-ro агрегата.

Естественно, О

^ а 1 ;

 

 

2

а « = 1>

 

ke S j

 

где Pj — множество

агрегатов

(включая и выходы ком­

плекса), на которые поступает выходной поток данного /-го агрегата; a.kj могут быть как заданными величина­ ми, так и искомыми переменными задачи.

Вместо введения переменных отдельно для входных хМ и выходных у(г) потоков агрегатов, как это сделано в уравнениях (1-1) — (1-3), можно вводить другие пере­

м енны е^,

относящиеся к расходным параметрам по­

токов между

агрегатами. В данном случае x Jtr —раз.

ходный параметр потока между i-м выходом â-ro агре­ гата и /-м входом г-го агрегата. В случае отсутствия потоков между определенными входом и выходом дан­ ных двух агрегатов соответствующая переменная не вво­ дится. Такая система переменных позволяет сократить общее число переменных в модели. Однако в таких обозначениях сложно записать уравнения материального баланса.

Модель комплекса можно преобразовать в эквива­ лентную, не обращая внимания на ее внутреннюю струк­ туру. Предположим, для общности, что ограничения (1-4) заданы в виде некоторой системы равенств, на­ кладываемых на соответствующие значения входных и выходных переменных. Если средщ ограничений есть неравенства, то добавлением новых переменных их мож­ но превратить в равенства. Тогда для конкретного ком­ плекса можно, подставив (1-1) в уравнения (1-2) и кри^ терий (1-5), исключив ряд переменных и введя новые переменные для расходных параметров материальных

потоков хj (/=1,

. .

R)

и для

управлений

по режим­

ным параметрам

us (s= l,

5), записать ограничения

(1-1)— (1-4) в виде

 

 

 

 

(^î» • • ■»

 

Wj)—0, i— 1, ...,

/г,

и критерий (1-5)

* ^ 0 ,/ = 1 ,

 

(1-9)

 

 

 

 

 

F(xi,

...,

xR, ut, . .

ив) - н т п .

(1-10)

Преобразование задачи к виду (1-9), (1-10) не упро­ щает решение, а служит лишь для стандартизации за<

писи и применения общих методов нелинейного програм­ мирования. Специфика преобразованной модели (1-9), (1-10) будет рассмотрена в гл. 2.

Рассмотрим еще одну модель для приближенного описания производственных возможностей агрегатов комплекса через режимы работы каждого агрегата [9] на примере к-то агрегата с одним входом и многими выходами (рис. 1-5). Предположим, что k-й агрегат мо­ жет работать в 7\ различных режимах и для каждого г-го режима известна производительность агрегата по

выпуску /-го продукта и[р 1\отнесенная к единице входно­

го потока.

Рис. 1-5. Пример агрегата с одним

х (f<L

.и * >

входом и т к выходами.

•Ут

 

Следовательно, если k-\\ агрегат переработает х[к)

входного продукта на г-м режиме, то выход /-го продук­ та составит:

 

,№) _гАгк)у{к)

/ __1

т -

( 1- 11)

 

Уjr

llj

Х г

' / — 1 > • • • •

m k1

общий выход /-го продукта в течение всех Th режи­

мов работы агрегата будет:

 

 

 

#//г)= 2

^-Г/г,^ /г)* / = ) . •••>

Щ.

( 1- 12)

 

Г =-1

 

 

 

 

 

причем

 

 

 

 

 

 

 

 

2

4 к)= х (к).

 

 

 

 

г=I

 

 

 

Таким

образом,

в

задаче оперативного

управления

с моделью

(1-11),

(1-12)

необходимо

найти значения

х (к), т. е.

интенсивности

использования отдельных ре­

жимов в работе /г-го агрегата.

В дальнейшем основное внимание будет уделено по­ становке задачи оперативного управления в форме

(1-1)— (1-5) или (1-9), (1-10), при этом, если это специ­ ально не оговаривается, будут исследоваться задачи опе­ ративного управления в условиях полной информации о параметрах модели комплекса и в предположении, что полученное решение может быть точно реализовано, хо­ тя это и является известной идеализацией.

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

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

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

ративного

управления. Некоторые

примеры

подобных

постановок даны в гл. 7.

 

 

Существенное предположение при формулировке за­

дач (1-1) —(1-5) или (1-9), (1-10)

состоит

в том, что

агрегаты

описываются статическими зависимостями.

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

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

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

Г л а в а в т о р а я

МОДЕЛИ ОПЕРАТИВНОГО УПРАВЛЕНИЯ

Алгоритмы решения сформулированных в гл. 1 задач оперативного управления в значительной мере опреде­ ляются тем, к какому классу задач математического программирования будут сводиться постановки (1-1) — (1-5) или (1-9), (1-10).

Делая различные предположения относительно свойств функций F и G в задаче (1-9), (1-Ю), можно получить различные преобразования модели оператив­ ного управления работой комплекса. Рассмотрим неко­ торые классы этих моделей и параллельно связь свойств преобразованной модели с видом характеристик отдель­ ных агрегатов j<hl

2-1. Линейная модель

В случае линейной модели комплекса формулы (!-§)< (1-10) запишутся следующим образом:

 

aijXj~‘bt==0'

........ пi

 

 

у-1

 

 

 

/?о

j — 1, .... R,

(2-2)

 

C/JC/—-min, Xj>Q,

 

/=i

 

 

переменные

заменены на переменные х с индексами

/= t/?+ 1,

/?+5, i?o=^+*S, a{j, 6г-, Cj — заданные числа.

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

Подобные линейные задачи применительно к задачам распределения материальных и энергетических потоков рассматривались во многих работах [5, 8, 9, 12, 23—32]. При этом существенное внимание было уделено методам расчетов линейных моделей для систем большой раз­ мерности [20, 33].

Предположим, что работа отдельных агрегатов опи­

сывается линейной моделью вида

 

 

уГ = 2 aVixT ' ‘= 1 ........ Пк'

(2-3)

 

t=I

 

 

FM = 2 c f x f \

(2-4)

 

1= 1

 

где

— затраты на функционирование

&-го агрегата;

а{!9, с\к) — заданные числа; множества Х<0), Y(e), Х(/г\ Ÿ<&)

4

1

»

V определены системами линейных равенств и нера­

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

(1-5) затраты FW вхо­

дят линейно.

Тогда из (1-2), (1-3) следует, что преобразованная модель (1-9), (1-10) будет моделью линейного програм­ мирования. Это справедливо и для систем с рециркуля­ цией [5].

В (2-4) предполагается, что затраты на функциони­ рование /г-го агрегата представляются в виде суммы

йЬстбянных затрат, которые нё учитываются в задачё оптимизации, и затрат, линейно зависящих от объема перерабатываемого сырья или объема готового продук­

та. Таким о б р а з о м , в (2-4) — условно-переменные

затраты, отнесенные к единице перерабатываемого сырья.

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

Были разработаны специальные алгоритмы примени­ тельно к оптимизации материальных потоков для линей­ ных задач с двусторонними ограничениями, называемых задачами интервального программирования [31].

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

Обозначим через x skl)L{ расход 5-го компонента, посту­

пающего с выхода /г-го агрегата на вход /-го агрегата. Каждый из компонентов может подаваться на вход лю­ бого из имеющихся N агрегатов. Если предварительно известно, что менаду k-ы и >/-м агрегатами нет связи по

5-му компоненту, то соответствующее x f r) приравнива­

ется нулю.

Суммарный поток s-ro компонента, поступающий на вход /-го агрегата с остальных N агрегатов, можно за­ писать:

(2-5)

а суммарный поток s-ro компонента, уходящего с /-го агрегата, запишется как

(2-6)

к

где x:(0l)—входной ноток комплекса по компоненту 5.

S

Таким оёразом, при наличии iia входе и выходе /-го агрегата г компонентов линейная модель для i-ro агре­

гата

будет иметь вид:

 

 

 

у(П=А(Ох(«,

(2-7)

где А(^ = {я^}} — матрица

технологических

коэффициен-

тов;

у& = {у\,}........ у].0};

хС1>= {лг{°, ....

JCJ0} — векто­

ры.

Подставляя в уравнения (2-7) выражения (2-5) и (2-6), получаем систему уравнений материального ба­

ланса, выраженную через переменные х[к1К

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

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

Для первого агрегата имеем на входе компоненты 1 и 2 и компонент 3 на выходе. Система (2-7) перепишется следующим образом:

1,0) _«(1)

хО)

I

fl(i) х0)

Уз

— “31

Х 1

^

а32 х2

Для второго агрегата имеем на входе компоненты 1 и 3 и ком­ поненты /, 2 и 4 на выходе, из которых 1 и 2 подаются обратно на вход первого агрегата. Система (2-7) примет вид:

,,(2) _ (2) (2) , (2) (2).

У1 — “п x i + “ 1з хз >

*!2)+ 4 1 ’ 4 2>;

4 2>= «‘?> *12’ + «<§>

Принимая во внимание, что согласно рис. 2-1

x|,> = x p

+ *i2l>;

=

хр

=

+

у ^ =

х|М )

*S2' =

 

 

(2-8)

 

4 2,=*Г>

Ч2,“ *3

 

(2) _

v(20)

»

УХ

/

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