книги / Типовые задачи оперативного управления непрерывным производством
..pdfдукта j). Поэтому ограничения на допустимую четкость выписы* ваются в виде
|
|
Vs+\,j— |
|
где v*j — заданное значение. |
|
||
по |
Каждый показатель качества целевого продукта определяется |
||
фракционному составу или по значениям выхода (отгона) v,j |
|||
при |
заданных температурах ts. Поэтому ограничения на показате |
||
ли |
качества (как |
правило, односторонние) можно записать: |
|
|
|
»}‘ >(Vs/) |
. |
где |
ÿJ.Q5 — заданное |
значение. |
|
|
При переработке нефти обычно |
задаются ограничения на сле |
дующие показатели качества целевых продуктов: температура кон ца кипения бензина, температура перегонки 90%-ного керосина, температуры перегонки 50, 90 и 96%-ного ДТ, температура вспыш
ки керосина, |
температура вспышки ДТ, температура застыва |
|
ния ДТ. |
|
|
В |
качестве критерия оптимизации процесса разделения можно |
|
рассмотреть |
прибыль [см. (4-80)], однако при первичной перера |
|
ботке |
нефти |
стоимости целевых продуктов — бензина, керосина и |
ДТ, как правило, равны, а стоимость мазута существенно меньше. Поэтому в качестве критерия оптимизации выбирается суммарный выход светлых нефтепродуктов
|
|
|
3 |
|
|
|
|
|
|
F= S р /К 'К ’ |
|
|
<4*98> |
||
|
|
|
/=1 |
|
|
|
|
KOTODbifi максимизируется. |
|
целевого |
продукта, |
зависящая от |
|||
В |
(4-98) |
pj — плотность /-го |
|||||
фракционного состава. |
|
|
|
|
|
||
Будем полагать, что кривые непрерывного фракционного со |
|||||||
става нефти можно аппроксимировать значениями |
при |
темпера |
|||||
турах, приведенных в табл. |
4-1. |
|
|
Таблица 4-1 |
|||
|
|
|
|
|
|
||
Г |
I |
2 |
3 |
4 |
5 |
G |
7 |
ts |
145 |
165 |
206 |
225 |
285 |
330 |
360 |
Ограничения на показатели качества целевых продуктов выра зим через объемные относительные количества фракций, выкипаю щих до температуры (это значения по характеристике НТК), и через количества тех же фракций по разгонке без флегмы (обо значим их ?nj). Теперь ограничены следующие параметры:
конец кипения бензина г~)2^0,98; перегона 90% керосина Ргз^О.ЭО; перегонка 50% ДТ, z>3s^0,5; перегонка 90% ДТ г>зб^0,9; перегонка 96% ДТ 037^0,96;
температура вспышки керосина
р!= / ко—jМ 22^28°С; температура вспышки ДТ
/д= ;Д о+ * У з4^65°С;
температура застывания ДТ ^з_-/з0_|_^3у34— 11°С.
Коэффициенты k2, к, кз, k'a определяются на основе стати
стических данных.
Зависимость vsj(vsj) нелинейна, однако в требуемом диапа зоне изменений vSj ее можно линеаризовать
V s J = P SJ—J—( l-f-G sj) Vs j .
Ограничения (4-97) по фракциям для всех температур, приве денных в табл. 4-1, запишутся следующим образом:
|
|
|
V u U t + |
У12Ц2 = |
»1 ; |
|
|
|
|
|
|
|
|
|
о21ц, - f |
v Z2i h |
— 02Î |
|
|
|
|
|
|
|
|
|
~Ь 032ü 2 “Ь 03 3U3 = |
|
0 3 Î |
|
|
(4-99) |
|||
|
|
|
11л -f- V,xZllz -f- V,i3U3 = |
O.j î |
|
|
|||||
|
|
|
I |
|
|
||||||
|
|
|
ил -J- и2 -f* v 53t i 3 ^ |
o5; |
|
|
|
||||
|
|
|
«1 4 “ ^2 |
0G3«3 ^ |
00 > |
|
I |
|
|
||
|
|
|
U l 4 “ 0 2 |
*4* 0 7 3 ^ 3 |
t>7* |
|
) |
|
|
||
как |
Последние три условия в (4-99) |
имеют вид неравенств, так |
|||||||||
часть фракций, выкипающих |
от 285 до 360°С, содержится |
||||||||||
в мазуте, который не входит в |
целевую |
функцию |
[см. |
(4-88)]. |
|||||||
Естественно, что |
|
|
|
|
|
|
|
|
|
||
|
|
|
|
V s / ^ |
Vs j |
^ |
O.ç/ » |
|
|
|
|
где |
vSj, |
Vsj — заданные |
пределы |
изменения |
переменного |
uej. |
|||||
|
Ограничения на четкость разделения |
|
|
|
|||||||
|
бензина v2\— |
0 1 1 ^ 0 * 1 ; |
|
|
|
|
|
|
|
|
|
|
керосина t>42—032^ 0 *2; |
|
|
|
|
|
|
|
|||
|
ДТ |
O73—Оез<о*з. |
рДовД, |
получаем для |
функционала |
||||||
|
Линеаризуя |
зависимости |
|||||||||
(4-98) |
Р— (рю—Ар10ц)-Ь(р2о—АргО^-НДр^Озг) X |
|
|
||||||||
|
|
|
|
||||||||
|
|
X «2+ (рзо—АраП^-Др^бз) «з— ►max, |
|
(4-100) |
где р, Ар — коэффициенты.
Задача оптимизации разделения нефти состоит в максимизации функционала (4-100) при вышевыписанных ограничениях и являет ся аналогом задачи разделения в относительных единицах [см. формулу (4-81)] при фиксированном значении х.
Полученная линейная задача с переменными коэффициентами стандартным образом сводится к задаче линейного программиро вания.
Vl |
v % |
V 3 |
Vi |
V s |
V Q |
v 7 |
t K |
*5 |
k 2 |
|
0 % |
||||||
|
0 , 1 6 3 |
0 , 2 0 8 0 , 2 8 8 0 , 3 3 6 |
0 , 4 1 4 0 , 4 6 6 |
0 , 5 0 2 |
4 4 |
1 1 8 |
9 , 3 4 |
|
1 0 0 |
||||||
|
|
|
|
|
|
|
Продолжение |
т а 6 л. |
4-2 |
|||||
k |
* 3 |
/г' з |
PlO |
|
P20 |
Рзо |
|
д рг |
Д р а |
Д р ' а |
Д Рз |
д р ' з |
||
|
|
|
||||||||||||
2 0 0 |
1 8 |
1 8 |
0 , 7 2 2 |
0 , 7 9 9 |
0 , 8 6 3 |
0 , 0 3 5 |
0 , 0 2 5 |
0 , 0 2 |
0 , 0 2 |
0 , 0 3 5 |
||||
|
|
|
|
|
|
|
Продолжение |
т а б л. |
4-2 |
|||||
J J » |
V n |
|
|
ÜQI |
V 2 i |
V 2 2 |
v 2 2 |
|
^ £ 3 |
V 23 |
_ V j t |
|||
0 , 7 5 6 |
1 |
0 , 9 2 7 |
1 |
0 , 0 4 8 |
0 , 1 6 |
0 , 1 7 |
0 , 4 4 |
0 , 7 4 3 0 , 9 6 5 |
0 , 7 7 8 |
|||||
|
|
|
|
|
|
|
П р о д о л ж е н и е |
т абл . |
4-2 |
|||||
V 21 |
J h t |
Ü33 |
V p i |
|
V u |
V_35 |
|
УЗС |
V z o |
|
V3G |
V 37 |
|
^37 |
1 |
0 , 0 3 6 |
0 , 2 2 |
0 , 7 9 5 |
0 , 2 6 5 |
0 , 4 1 |
0 , 5 3 3 |
0 , 7 3 5 |
0 , 9 5 5 |
0 , 8 6 3 |
|
1 |
|||
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Таблица 4-3 |
|||
z |
|
I I I |
«2 |
|
|
«3 |
V |
n |
V 12 |
|
Uai |
|
V 2 2 |
|
|
|
|
|
|
|
|
||||||||
0 , 4 7 9 |
0 , 1 5 4 |
0 , 1 3 4 |
0 , 2 4 3 |
0 , 9 2 |
0 , 9 8 5 |
0 , 1 6 |
0 , 4 1 6 |
Продолжение табл. 4-3
V 2 3 |
V n |
V n |
V 3 i |
^33 |
ü 38 |
V 3 I |
|
||||||
|
|
|
|
|
|
0 , 7 4 3 |
0 , 9 1 0 |
0 , 1 4 1 |
0 , 2 4 7 |
0 , 5 2 0 |
0 , 7 3 5 |
0 , 8 8 5 |
|
|
|
Исходные значения коэффициентов для расчетного примера приьецены в табл. 4-2, результаты решения ’— в табл. 4-3. Зна чения v i—vj определяются составом конкретной нефти по харак теристике ИТК. Расчет характеристики ИТК нефти проводится по результатам предыдущей переюнки и анализа полученных светлых нефтепродуктов.
Если задача линейного программирования, к которой сводится сформулированная задача оптимального разделения, не имеет допустимого решения, это означает, что из нефти данного состава при достигнутой на данной установке четкости разделения невоз можно получить нефтепродукты заданного качества. В этом слу чае для получения нефтепродуктов заданного качества необходимо изменить состав нефти, например, путем смешения нефти из двух месторождений.
Глава пят ая
СЛОЖНЫЕ СХЕМЫ СОЕДИНЕНИЯ ТИПОВЫХ ТЕХНОЛОГИЧЕСКИХ ОПЕРАЦИЙ
5-1. Сложная операция
Будем рассматривать следующую модель сложной операции:
nk |
|
|
|
|
№) _ v i |
.,(*) |
* |
|
(5-1) |
# г = 2 |
uîiх! |
|
||
/=i |
|
|
|
|
где k —номер сложной операции в комплексе; |
— |
|||
входной и выходной |
потоки |
операции; |
—t-я состав |
ляющая вектора управления ujA), связанная с перераспре
делением /-го входного потока между mk выходными по токами k-vL операции.
Предположим, что вектор ujA) выбирается произвольно
из выпуклого ограниченного множества üjA);
G II!*’. |
(5-2) |
Это множество определяется ограничениями на управления, например, показателями качества выходных потоков.
Рассмотрим заДЦчу минимизации себестоимости. Тре буется найти такие значения переменных и[®, x{.k), кото
рые удовлетворяют исходным ограничениям (5-1), (5-2) и ограничениям по выходным потокам
ym = yW |
(5-3) |
и обеспечивают минимум функционала
F = S a\ki x? ' |
(5-4) |
/-* |
|
где и!®-- себестоимость переработки |
единицы /-го вход |
ного потока — величина переменная, |
иУ? U*é). |
Сформулированная задача является задачей обоб щенного линейного программирования и может быть сведена к эквивалентной задаче линейного программи рования.
С учетом критерия (5-4) возможна постановка задачи минимизации себестоимости при фиксированных входных
потоках |
и ограничениях на входные потоки |
вида 0 < y{k) < y^\
В этом случае задача минимизации себестоимости превращается в задачу обобщенного линейного програм мирования с переменными
Уш, |
i — 1......mk> s = 0, 1....... |
mk, j= 1....... |
nk. |
L |
SJ |
|
|
Однако, если вместо ограничений (5-2) имеются дру гие ограничения на коэффициенты щ например линей ные зависимости для коэффициентов в разных столбцах, задача становится существенно нелинейной.
Сложная операция является общим случаем по от ношению к другим типовым операциям и технологиче ским комплексам соединения этих операций. Однако структуры комплексов накладывают определенные осо бенности на их модели, а следовательно, на задачи опе ративного управления для таких комплексов. Так, зада чи, рассмотренные в гл. 4 для последовательных и параллельных схем соединения типовых операций,
в большинстве случаев сводились к Задачам линейного программирования. В этой главе рассмотрим два класса часто встречающихся структур комплексов: последова тельно-параллельные схемы и схемы с рециклом, кото рые приводят к нелинейным задачам.
5-2. Последовательно-параллельное соединение типовых операций
Предварительно кратко остановимся на последова тельно-параллельном соединении простых операций. В общем случае задача минимизации издержек для по следовательно-параллельного соединения простых опера ций вида' (4-1) сводится к задаче линейного программи рования с переменными коэффициентами. В столбце переменной у№ будут переменные коэффициенты //i(0(7t)),
х} |
Рис. 5-1. Пример схемы |
ком- |
Jz) |
плекса, состоящего из последо- |
|
" |
вательно-параллельного |
соеди |
|
нения типовых операций. |
|
(3) / —.простая операция; 21 3 — сме-
Цсительные операции.
-зависящие- только от одного режимного пара метра 0(Ч Подобным свойством обладают задачи, в ко торых нет влияния показателей качества вы ходного потока k-й операции на процессы в последую щих операциях. ’ В этом случае применение метода аппроксимации, рассмотренного в § 4-1,а, строго обо сновано.
• При наличии влияния показателей качества одной операции на последующие переменные коэффициенты
/>.(0(ft))» .£л(0(*?), p(ft)(0^fe)) будут одновременно входить в: несколько столбцов и вышеуказанная кусочно-линей
ная аппроксимация не применима.
Аналогичная ситуация возникает при исследовании последовательно-параллельных схем с другими типовы ми операциями.
Рассмотрим комплекс, состоящий из простой onepajции и следующих за ней двух параллельных операций смешения (рис. 5-Î).
Простая операция будет определяться следующей линейной моделью с переменным коэффициентом:
</<’>= foç<4; о<;> = 0Р!1».
гдей,'1, oJ| —качественный показатель входного и выходного потоков; 6 —параметр управления,
0|( 62 —заданные^числа.
Модели смесительных операций будут представлены линейными моделями:
х |
у( 2). |
.(3)
X + 4 3,= ! Л
«4” = о\?х\2>/ут 4-o.Vvf’/ÿ(2,; и*3*= v«'x?'ly«' +Ъ!гх?1ут,
где vl2, Ü22— заданные качественные показатели.
Рассмотрим задачу оперативного управления комп лексом, состоящую в максимизации прибыли
F=c™y{2)+ |
,(3>„(3> |
(I) Д1) |
- |
,(2) Д2) |
с™у™ - |
с\"х™ |
|
||
. Д2) |
(2) _ |
(3) (3) |
|
(3) (3) |
1*0 |
Ой Ai |
|
AQ |
при ограничениях
у п > у ( * ) , y m s z y W ;
( 1 ) |
“ (I) |
v <2) |
(2) |
~(3). |
X |
|
г2 |
X |
»V(3)( |
|
|
|
о\(3) |
V (3) |
Следует учесть также ограничения связи
У(0. |
«Г + д ? ’; |
( 2) |
:Ц(3) = 0<1) |
II |
ип —иП |
Окончательно ограничения |
задачи |
можно |
записать |
|||||||
в следующей форме: |
|
|
|
|
4 3>< x m; |
|||||
|
|
|
|
|
|
|
||||
jt<‘> |
|
|
|
|
|
|
|
|
< x (,>; |
|
0*(ï) |
—*<2) |
|
|
|
|
|
|
= 0 ; |
||
|
|
|
|
+ |
*Г |
|
|
> Р > ; |
||
|
а$)'х™+ |
|
|
4 2> |
|
|
< P > ; |
|||
|
|
|
a2x (22) |
|
|
> 0 ; |
||||
|
|
|
|
aa(d) x<? -j- 04 X2^ |
|
0. |
||||
|
|
|
|
|
|
|
|
_ |
|
(5-5) |
я» (0)=(—о!2) + » i 10); |
« 2 = (— y!i2)+ |
у12); |
|
|||||||
«3 (0)= |
(—v f + |
oi0 0); я*= |
(—ÿ f + âaa). |
|
||||||
Переменный параметр 0 входит в различные столбцы |
||||||||||
матрицы |
ограничений |
(при неизвестных л;(1\ х ^ |
и л:|3)). |
|||||||
Введем новые переменные Ai и Аг, тогда |
|
|
|
|||||||
0==01A1-f-(j2A2, |
*, + |
*, = |
1, |
Я ^О , Я2^=0. |
(5-6) |
|||||
Если подставить (5-6) в (5-5) и провести замену пе |
||||||||||
ременных |
Я1х11) = хп, Xtx W = x lit |
|
|
|
|
|
||||
|
A,j:j2) = jc{f, |
|
||||||||
3 |
r (2 ) ___ r (2) |
j (3) |
_ _ |
(3) |
2 |
(3) ___ |
(3) |
|
|
|
Л 2Л 1 |
---- ^ 1 2 » Л 1Л ] |
----- ^ J J » |
Л 2Л 1 |
----- Л |
12 » |
|
то в (5-5) можно использовать следующие подстановки:
0л:(,) = Ojjc,, -|- в2л:12;
Дополнительно в (5-5) появятся следующие ограни чения:
X tI |
|
= |
xW; |
|
r (2) |
|
|
|
|
л п + * ! ? ..— ,Л 1 |
> |
|||
x (3) |
|
= |
:л:(3) |
|
Л и |
+ |
* I2 ) |
L |
|
|
A 1 |
|
и в силу (5-6)
(5-7)
Наличие ограничений (5-7) делает поставленную за дачу нелинейной.
В исследуемом комплексе первая операция описыва ется линейной моделью с переменными коэффициентами. Если эту простую операцию заменить на смесительную с линейной моделью, то и в этом случае аналогичная за дача распределения материальных потоков для всего комплекса будет нелинейной.
Таким образом, несмотря на относительно простые модели отдельных операций, общая модель последова тельно-параллельного комплекса при учете влияния ка чественных показателей первой операции на последую щие является существенно нелинейной.
Вышеприведенную частную задачу можно решить как линейную параметрическую (0-параметр). Однако в бо лее общем случае, при нескольких параметрах 0; реали зация подобного метода существенно затруднена.
Можно использовать и другой подход, связанный с представлением в исходной задаче произведения пере менных в виде сепарабельной функции с последующей аппроксимацией этой функции [63], например:
(Ы1*
где 6. |
JC<‘>—максимально возможные |
значения перемен |
ных 6 |
и |
|
Каждая из новых переменных 2г- |
аппроксимируется |
в г точках, для этого вводится группа из г новых пере менных
гг
k= 1 |
It~\ |
Для каждого произведения переменных в исходной задаче необходимо ввести 2г новых переменных и четы-
ре ограничения вида (5-8). Таким образом, при указан ной аппроксимации происходит существенное увеличе ние размерности задачи.
5-3. Технологический комплекс с рециклом
Рассмотрим простейшую схему технологического ком плекса с рециклом (рис. 5-2). Такой комплекс включает смесительную, простую и разделительную операции.
Рис. 5-2. Пример схемы комплекса с рециклом.
1 — смесительная операция; 2 — простая операция; 3 — раздели тельная операция.
Смесительная операция описывается моделью x0+ yf]= y 0);
ü!W 3) + ^*o ==<'*/
—„(‘U 1». 02lV!3)+ £,e2JC* = t,2 '^
V ( 1 ) |
V |
ü(1). v |
!) |
V,( 1) |
V ( I ) |
||
|
|
u\ |
» u2 |
|
|
||
,(D |
“ (И |
„(И |
, |
“ (1) |
— заданные величины, ха |
||
где и t, vv2, v\', |
t»j . |
|
v2 |
рактеризующие качественные показатели процесса или концентрации веществ в потоке.
Простая операция описывается моделью ут = e<*>x<*>; 0<2) < 0(2) С в ;2',
где в,21, б[21—заданные числа.
.Для разделительной операции выбрана следующая модель:
в ? > 7 ? ;
"aVi3’ + v'£ y? = < £ V ;
уТ + У ? = хт\
v<3> |
=ц(^ < о (3), s ,i= |
1, 2. |
|
s i . |
si |
si ’ |
’ |