книги / Типовые задачи оперативного управления непрерывным производством
..pdfгде auh=fii(uik), Cih=Fi(uih), Xiu— новые переменные,
T
удовлетворяющие условию 2 AIA= aV
Л=1
Подставляя (2-36) в (2-35), получаем аппроксими рующую задачу
т R 1
2 |
ai^kxki “h 2 aiix i = |
^ ’ Xkl’ x i ^ |
k=\ |
/=2 |
(2-37) |
T |
R |
2 |
+ S c/*/“*mîri* |
|
k=l |
/«2 |
|
Из решения задачи (2-37) определяются значения |
||
переменных x*i/£ и режимного параметра |
|
|
|
« * ,= 2 И**Я*1Л, |
(2*38) |
|
k= \ |
|
где k*ih‘ —x*ihfx*i.
Для того чтобы определение M*I в (2-38) соответст вовало принятой в (2-36) точности аппроксимации, не обходимо, чтобы в решении (2-37) были отличны от ну ля не более двух смежных значений Xi/t:
Xi л, ^1,/i+i^O, xsi = 0 , |
s = 1, . . . , k — 1, |
k + 2 , . |
., 7\ |
(2-39) |
Согласно [20] это условие выполняется, если |
fn(ui) |
|||
являются линейными |
функциями от |
строго |
монотонной |
функции f(ui), a F\{Ui)— строго выпуклая функция. В [40] изложенные приемы аппроксимации примене
ны в случае нескольких управлений. Аппроксимирующая задача (2-37) допускает следующую техническую интер претацию. Процесс преобразования первого входного по тока Xi можно осуществить Т способами, причем Xiu— интенсивность, сц{— стоимость k-то способа, /г=1, ..., Г. Возможен любой промежуточный способ, который опре деляется как средний с весом Xi/£=Xi/i/Xi.
При помощи простой интерпретации можно непо средственно формулировать линейную аппроксимирую щую задачу (2-37), не приводя исходной нелинейной по становки (2-35), что и делается в ряде работ. При этом необходимо иметь в виду ограничения иа вид функции /п(и), F\(u) или нелинейное условие (2-39).
В [20] рассмотрена также задача, когда в одном столбце матрицы задачи переменные коэффициенты за
висят от двух управлений, т. |
|
е. в |
(2-22), (2-23) |
fn(u) = |
|||
“ fti(u\y ti2), |
Fi{u)=F\{tiu |
|
и2). |
Нелинейные |
функции |
||
f«(wi, и2), F i( u i , и2) |
могут быть |
аппроксимированы ку |
|||||
сочно-линейными на сетке |
|
значений щи, и2Р, k= \ , ... |
|||||
..., Ти р—1, |
..., Т2: |
т, |
|
т3 |
|
|
|
|
|
|
|
|
|||
|
^г) |
= S |
S |
|
|
||
|
|
k = l |
p |
= 1 |
|
(2-40) |
|
|
|
|
T, |
Ta |
|
||
|
|
|
|
|
|||
|
Fi faii ^a)*^i==2 |
|
S |
^йр^Лр* |
|
||
|
|
fc=I p=I |
|
|
Здесь переменные *iftp соответствуют интенсивности использования варианта со значениями управлений Щн, и2р и
"2iX,ke=X,.
ь,р
Выражения для }ц{щ, и2), F\{uu и2) из (2-40) под ставляются в (2-22), (2-23). Если в результате решения полученной задачи не более четырех смежных xibp^O, а остальные *iftp=0, то определяются приближенные значения щ, иг.
Если в задаче можно реализовать только одно реше ние, а ответом в аппроксимирующей задаче служат не смежные значения х\ьР, то технологическая трактовка полученных решений вызывает затруднения.
Аппроксимация на сетке практически может быть использована в случае двух режимных параметров для /-го столбца переменных коэффициентов задачи. Увели чение числа режимных параметров приводит к сущест венному увеличению размерности.
Включение агрегатов, описываемых моделями с пере менными коэффициентами вида (2-25), (2-26), в ком плекс агрегатов с такими же линейными характеристи
ками приводит |
к модели (2-22), (2-23) общего вида. |
В этих моделях |
переменные коэффициенты в разных |
столбцах могут быть функциями одних и тех же управ лений,
Приближенное решение задачи (2-22), (2-23) может быть основано, на аппроксимации зависимостей (2-25),
.(2-26) кусочно-линейными для Т фиксированных значе ний управлений, как это было сделано для задачи
(2-35). Аналогично по результатам решения аппрокси мирующей задачи определяется значение управления щ. Трудности трактовки результатов могут возникнуть при нескольких несмежных *«> 0, k= l, . . T.
Рассмотрим более общую модель системы с перемен ными коэффициентами [42]
/?
У/=Я|* + |
3 |
a l l X j , î = l , |
|
nr, |
(2-41) |
||
где |
|
/= « |
|
|
|
|
|
|
|
|
|
|
|
|
|
a i / = £ i / + S b‘!'a'’ ! = |
= |
I ............. ... |
l=°> |
!> • |
•» n; (2-42) |
||
r= 1 |
|
|
|
|
|
|
|
ar— независимые переменные коэффициенты |
(парамет |
||||||
ры), r = l, .... L\ dij, bijr — заданные постоянные коэффи |
|||||||
циенты. |
|
|
|
|
|
|
|
Подставим |
(2-42) |
в (2-41) .и |
введем |
обозначения |
|||
|
|
|
/? |
|
|
|
|
* /г (Х) = |
6/ог + |
2 |
bUrX i> * — |
1. • • •. т • |
|
||
|
|
/= 1 |
|
|
|
|
|
Тогда (2-41) можно записать как |
|
|
|
||||
|
/г |
|
L |
|
|
|
|
Уi= ^ io + S |
+ |
Х агг/г(х), |
t = l , . |
/я, |
|||
/=1 |
|
г=,1 |
|
|
|
|
|
|
|
0 < а Г< 1 , г==1, .... |
L. |
(2-43) |
|||
Рассмотрим свойства этой модели на примере агрегата, имею |
||||||||
щего два входа (/?= 2), три |
выхода |
( т = 3) |
и |
два переменных ко |
||||
эффициента |
(L = 2). |
|
|
|
|
|
|
|
Уравнения в этом случае запишутся так |
|
|
|
|||||
Уг = |
flio + |
(<*!!*, 4 |
о12Хг) 4 |
a,z„ (х) 4- a2z12 (х); |
] |
|||
у 2 |
= |
д20 4 |
(a2iXi 4 |
а22ха) + |
“izsi (х) 4 |
«г^гг (х) ; |
> (2-44) |
|
«/з = |
4 |
(<*31*1 + |
032Х2) 4 |
«1*31 (*) + |
°2г32 (х). |
) |
||
где 0 ^ а |
, |
1, £ = |
1, 2; |
|
|
|
|
|
zi\= b ioi4 6 11t*i4 *121X2;
2|2=&1ог4^112ХI4 ^I22X2Î
Z2i=&20l4 ^2llXl4 ^22lX2;
222=&2024^212Xl4^22lX2;
Z31=&30l4ft«nXi4Ô32lX2;
z32= ^3024^312Xl4^M2X2-
Если каждый из переменных коэффициентов входит лишь в определенный столбец л\-, то, как уже отмечалось, исходная за дача может быть сведена к линейной путем замены переменных. Для нашего примера такой случай возникает, например, когда все Ък02 равны нулю, 6121=Г»221=^Э21=0 И &112= 6212= 6з12= 0. В ЭТОМ
случае
У\— аю+ (ÛI1х\-\-а\2X2)-f-(x1bшх î+ОгЬ122X2
у2= <*20+ |
(a2lX\-\-Q22x 2) —j—Oti Ьц\Х |
222*21 |
|||||
# з = а з г + |
( « 31* 1+ |
|
032* 2) +-<Xi bai\x i+ ct2& 322*2 - |
||||
Вводя новые переменные Mi=ai.vi, «2= 012*2, приходим к следую |
|||||||
щей эквивалентной линейной модели: |
|
|
|
||||
Ух — « ю |
+ |
( « u * i |
+ |
« 12* 2) |
+ |
b n iU i + |
Ь122и 2; |
Уг — « го |
+ |
( « 21*1 |
+ |
« 22* 2) |
+ |
bilx ut + |
b22iaz; |
У3= а з0 - f |
(a3ix j + |
азгх 2) |
+ |
bm a , + |
b322u2; |
||
()<;«, |
х г; 0 |
|
x 2. |
|
|
Однако можно указать более общее условие сведения полученной модели к линейной. Согласно [42] модель может быть -сведена к линейной, если
|
|
bijr-—ëirCrjt |
(2-46) |
где gir, |
Crj— элементы матриц G и С, имеющих размер |
||
ность соответственно |
(mXL) и |
( LX( R4-1)). |
|
Для |
нашего примера |
размерность |
матриц (3X2) и (2X3). |
Используя (2-46), получаем: |
|
2ll=gü (Cio4-Cll*l4-CI2*2) î Zl2=£l2(C2o4-£21*l4-C22*2) î
Z2 i= g 2 i (cio-j-cnXi-j-cisXÿ) ;
Z22~ë22(C2Crj-^2]Xi-j-Cz2X2)
Z 3 l~ g z i («10+ « 1I*1+Cl2*2)î
232—g32(020+C2I*I+c22*2)• Вводя новые переменные
Hi=ai(c,o4-cnX|-fci2X2), «2= a 2(c2o-|- +« 21*1+ 022*2),
можно при условии постоянства знака функций Oio+«n*i+«i2*2 п C20+C2Îх 14 "с22Х2 в области изменения Х\ и х2 перейти к следующей
линейной модели:
Уг — |
«10 + « 1 1 * 1 |
+ |
« 1 2 |
* 2 + £ l l « I + g l 2 « 2 ; |
|
||
Уг — |
«20 + « 2 1 * 1 |
+ |
« 2 2 * |
2 + |
i» 2 1 « l + |
^ 2 2 « 2 F |
(2-47) |
|
, |
, |
«32*2 + |
, |
I |
||
Уг — «30 + «31*1 + |
g3i«i + |
й зги 2, |
|
0 < « , < С , 0 + С П Х, + С12Х 2,
0 < ы 2< с20 + «21*1 «Я2*2* |
s |
М одель |
(2-47) |
получила |
название диапазонной. Действительно, |
||||||||||
при фиксированных входах Х\ и |
х2 значения |
выходов уи Уг> Уз ле |
|||||||||||
ж ат в |
диапазоне, |
определяемом коэффициентами |
g iT и |
граничными |
|||||||||
значениями щ и и2. З адавая параметры щ и |
и2, |
можно |
установить |
||||||||||
требуемые |
значения |
входов |
в |
пределах |
возможностей |
модели |
|||||||
(рис. 2-2). |
|
|
|
|
|
|
|
|
|
|
|
|
|
Диапазонные модели в общем случае можно записать |
|||||||||||||
в азиде |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
y = |
a0 + |
Ax + Gu; |
|
|
|
(2-48) |
||
|
|
|
|
|
0 < u < c 0 + |
Cx, |
|
|
|
||||
|
|
|
|
|
|
|
|
|
|||||
где ао, Со— векторы; A, G, С — матрицы коэффициентов; |
|||||||||||||
у >и х — векторы выходных |
и входных потоков |
объекта; |
|||||||||||
и — вектор переменных, обе |
|
|
|
|
|
||||||||
спечивающих |
диапазонную |
|
|
|
|
|
|||||||
зависимость |
выходов |
|
от |
|
|
|
|
|
|||||
входов. |
|
|
|
|
|
|
|
|
|
|
|
|
|
В (2-48) имеется неодно |
|
|
|
|
|
||||||||
значная зависимость выход |
|
|
|
|
|
||||||||
ных переменных у от вход |
|
|
|
|
|
||||||||
ных х в виде диапазонов |
|
|
|
|
|
|
|||||||
возможных |
значений |
выхо |
|
|
|
|
|
||||||
дов, |
причем |
выходы в пре |
|
|
|
|
|
||||||
делах |
диапазонов |
связаны |
|
Рис. 2-2. |
Характеристики |
диа |
|||||||
между собой, эта связь осу |
пазонной |
модели. |
|
|
|||||||||
ществляется |
через |
перемен |
|
|
|
|
|
ные и. Пример характеристики подобной модели для одномерного случая приведен на рис. 2-2.
Если вектор и имеет такую же размерность, как и у, то выходы в диапазонах между собой не связаны, т. е., варьируя технологические параметры агрегата, можно устанавливать любые значения выходов в преде лах диапазонов. Это варьирование будет сказываться на значениях переменных и. Если размерность и равна единице, то в пределах диапазонов каждому значению одного из выходов с ненулевым g будет соответствовать единственное значение любого другого выхода в своем диапазоне. Естественно, что когда С о = 0 и сг,==0, модель (2-48) превращается в обычную линейную модель.
Диапазонную линейную модель (2-48) можно рас сматривать как частный случай модели с переменными коэффициентами (2-43), поскольку последняя сводится к (2-48) при условии (2-46). Однако класс диапазонных моделей представляет самостоятельный интерес, так как
âTH моДелй ДоСтйтоЧйо хорошо оНиСывйют -работу реаль ных агрегатов -и имеют ряд преимуществ при экспери ментальном 'построении моделей.
Т а б л и ц а 2-1
Модель агрегата
Нелинейная
Линейная
С граничными ва риантами
Спеременными ко эффициентами
Диапазонная
Уравнение агрегата |
(зависимость |
Нрмер форму |
|||||
лы и парагра |
|||||||
вывода от входа н параметров) |
|
фа |
|||||
y(*) = |
f(fc)(X(A), |
и(/г), £<*)) |
и |
( l - о |
|||
|
|
|
|
|
|
в § 2-2 |
|
Ч |
|
|
|
|
|
|
(2-3) |
|
|
|
|
|
............ .... |
|
|
м |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
ч |
|
|
|
|
|
|
|
у [к) = 2 |
t f k h f K |
/ = 1, . . . . mk |
|
(1-12) |
|||
r= l |
|
|
|
|
|
и |
в § 2-4 |
|
|
|
Tk |
|
|
|
|
* ( * ) = |
2 |
xj.** |
|
|
|||
|
|
r=,l |
|
|
|
|
|
|
|
4 |
|
|
|
|
|
|
|
|
|
|
|
|
(2-25) |
|
|
/= i |
|
|
mk |
|
|
i' = |
l , |
. . . » |
|
|
|||
y = a 0 + A x - f G u |
|
(2-48) |
|||||
0 ^ |
U ^ |
Cg -j- Cx |
|
||||
|
|
В табл. 2-1 сведены различные виды моделей агре гатов, рассмотренные в гл. 1, 2, -и даны формулы, при водящие к линейным или нелинейным моделям ком плекса в целом.
2-4. Пример модели нефтеперерабатывающего завода
В качестве примера модели крупного непрерывного комплекса рассмотрим модель распределения материаль ных потоков производства топлив нефтеперерабатываю щего завода (НПЗ).
Подобная модель используется как для целей опе ративного управления, так и для тех-нико-экономическо-
го планирования [9, 12, 13]. Рассмотрим вариант ли нейной модели для задачи технико-экономического пла нирования.
Модель НПЗ строится на основе моделей отдельных технологических операций и процессов, связанных мате риальными потоками, и ограничений на входные и выходные потоки технологической схемы. Обычно произ-
Рис. 2-3. Схема основных материальных потоков НПЗ.
/ — установка первичной переработки нефти |
(ABT); 2 —-четкая ректификация |
|||||
бензина; 3 — газофракционирующая установка |
(ГФУ); |
4 — риформинг-1; |
5 — |
|||
риформинг-2; |
6 — гидроочистка дизельных |
топлив; |
7 — смешение бензина |
|||
АИ-93; 8 — смешение |
бензина А-76; 9 — смешение дизельного топлива |
(ДТ) |
||||
с содержанием |
серы |
0,5%; /0 —- смешение ДТ с |
содержанием серы 0,2%; 11— |
|||
смешение мазута. |
|
|
|
|
|
водство топлив НПЗ 1включает следующие технологиче ские процессы: первичную переработку нефти, при ко торой из -нефти выделяют бензин, керосин, дизельное топливо, мазут; риформинг узких фракций бензина, крекинг, гидроочистку дизельного топлива, фракциони рование газов и др. Основные товарные продукты: бен зин, дизельное топливо и мазут получают путем смеше ния из нескольких компонентов. Керосин получают на установках первичной переработки без последующего смешения. На (рис. 2-3 показана упрощенная технологи ческая ехема производства топлив НПЗ, включающая перечисленные процессы, кроме крекинга. В этой схеме цех первичной переработки нефти, состоящий из не скольких установок, условно показан как одна установ-
ка Î. На схеме показаны не все материальные потоки, а только основные.
а) Установки
Отдельные установки являются достаточно сложны ми технологическими объектами. Модели установок за даны в виде связей вход— выход и предназначены для решения задачи распределения потоков между уста новками.
Установки описываются линейными моделями с пере менными коэффициентами (коэффициентами отбора). Множество допустимых значений этих коэффициентов
|
|
|
|
Таблица 2-2 |
Гервичная переработка |
нефти |
|
|
|
|
|
|
Коэффициент отбора при |
|
|
Получаемый продукт |
режиме работы |
||
|
|
|
||
|
|
|
1 |
2 |
Сухой г а з .................. |
|
0,015 |
0,015 |
|
Головка |
стабилизации |
|
0,01 |
0,01 |
Бензин прямогонный |
|
0,1 |
0,10 |
|
Керосин |
• .................. |
|
0,125 |
0,125 |
Дизельное топливо . |
|
0,23 |
0,16 |
|
Фракция |
350—420°С |
|
0,07 |
0,07 |
Фракция |
420—490°С |
|
0,09 |
0,09 |
Мазут прямогонный , |
|
0,031 |
0,101 |
|
Гудрон |
...................... |
|
0,32 |
0,32 |
Потери |
...................... |
|
0,009 |
0,009 |
Ч еткая ректификация бензина
Получаемый продукт
Сухой газ .................................................................
Фракция, |
выкипающая до 62°С . . |
Фракция |
62— 105° С ............................................ |
Фракция |
105— 180°С . . . . |
Потери |
.................................................................. |
Таблица 2-3
Коэффициент отбора при режиме работы
î |
2 |
0,03 |
0 ,0 3 |
0,13 |
0,13 |
0,32 |
0,38 |
0,508 |
0,448 |
0,012 |
0,012 |
|
|
|
Таблица 2-4 |
Газофракцйонирую щ ая устан овк а (ГФУ) |
|
|
|
|
|
Коэффициент отбора тфи |
|
Получаемый продукт |
|
режиме работы |
|
|
|
|
|
|
|
! |
2 |
Газ c y x o ii ............................... |
|
0,11 |
0,11 |
Пропан-бугаиовая фракция |
. . |
0 ,1 8 |
0 ,1 8 |
Н -б у т а н ...................... |
|
0 ,2 3 |
0 ,2 6 |
И з о б у т а н ...................... |
. . . |
0 ,0 9 |
0 ,1 2 |
Газовый бензин ..................................................... |
|
0 ,3 7 7 |
0 ,3 1 7 |
Потери ...................................................................... |
|
0 ,0 1 3 |
0 ,0 1 3 |
Таблица 2-5
Риформинг-1
Продукт, получаемый при переработке фракции 105—180°С
Газ с у х о й ...................... |
....................... |
|
Головка |
стабилизации |
...................... |
Бензин |
риформинга . |
............................... |
В о д о р о д ........................... |
.................. |
|
Потери |
........................... |
.................. |
Коэффициент отбора при режиме работы
I |
2 |
0,073 |
0,073 |
0,028 |
0,055 |
0,834 |
0,807 |
0,046 |
0,046 |
0,019 |
0,019 |
Таблица 2-6
Риформинг-3
Продукт, поручаемый при |
Коэффн- |I |
Продукт, получаемый при |
Коэффи |
|
переработке фракции 62—105°С |
циент от-1I переработке ({тракции 62—105°С |
циент от |
||
|
|
бора |
|
бора |
Газ сухой |
........................... |
0,05 |
Бензол, толуол и др. . |
0,395 |
Горловка |
стабилизации |
0,05 |
Рафинат .............................. |
0,438 |
В о д о р о д ............................... |
|
0,035 |
Потери ............................... |
0,032 |
|
|
|
Таблица 2-7 |
|
Г идроочистка дизельного топлива |
|
|||
Получаемый продукт |
Коэффи |
Получаемый продукт |
Коэффи |
|
циент |
циент |
|||
|
|
отбора |
|
отбора |
Газ с у х о й ........................... |
топливо гид |
0,0158 |
О т г о н ................................... |
0,0089 |
Дизельное |
|
Сероводород ...................... |
0,0099 |
|
роочищенное .................. |
0,9555 |
Потери ............................... |
0,0099 |
определяется через пряничные режимы. Заданные зна чения (коэффициентов отбора на праничных режимах для рассматриваемых установок на рис. 2-3 приведены в табл. 2-2—2-7. Граничные режимы задаются при фиксированном составе сырья.
При использовании моделей с граничными (режима ми предполагается, что показатели качества получен ных продуктов одинаковы для всех режимов работы установок.
Переменные в модели планирования работы НПЗ — количество сырья Xk, перерабатываемое по &-му режиму. Интенсивность использования каждого режима Хи= =Xh/xо. Ограничения на производительность каждой установки записываются в матрице общей задачи .в виде
к
где суммирование осуществляется по всем режимам (b— максимальная производительность).
В такой форме ограничения представлены в матрице задачи линейного программирования.
б) Смешение
Товарные продукты: бензин, дизельное топливо, ма зут получают путем смешения компонентов, «получаемых на различных установках. Готовый «продукт должен иметь соответствующие ГОСТ или техническим условиям показатели качества, которые записываются в виде ограничений. Переменными модели являются количества Xj каждого компонента, поступающего в топливо данной марки. Для бензинов при планировании наиболее важ ны ограничения на октановое число, упругость паров и содержание серы. При текущем планировании обычно можно считать, что в этилированные марки бензина добавляется 'максимальное количество этиловой жидко сти. Тогда зависимость октанового числа смеси от окта новых чисел компонентов также будет линейной. Таким образом, все ограничения на показатели качества бен зинов в задаче предполагаются линейными и записыва ются в следующем виде:
Пх |
|
по октановому числу ^ (а*ч — агост) |
О, |
/с=1