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

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

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

где 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

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