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

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

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

дукта 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>42032^ 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

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