OPTIMIZĂRI STATICE ŞI DINAMICE
FUNCŢII OBIECTIV MULTICRITERIALE ŞI POLIOPTIMIZARE
METODE DE DECOMPOZIŢIE
STRUCTURILE MULTISTRAT ŞI MULTINIVEL
1. OPTIMIZǍRI STAŢIONARE
Exemplu: o instalaţie tehnologicǎ din industria chimicǎ
-R –reactorul
-S -schimbǎtorul de cǎldurǎ
-D –decantorul
-C-coloană de distilare
-PA-paleta de amestecare
-Aşi B– reactanţii - cu debitele QAşi QB
-debitul - fracţiunea recirculatǎ prin reactor
din debitul total QPB de produs de bazǎ al coloanei
-debitul - restul produsului de bazǎ,
folosit în calitate de combustibil
1
PB
Q
2
PB
Q
Din reactorul R, în care temperatura are valoarea T,
se obţine debitul total QR,
care însumeazǎ debitele a 6 componente rezultate din reacţie :
-debitele QRA şi QRB
corespunzǎtoare cantitǎţilor din reactanţii Aşi B
care nu au participat la reacţie
-debitele QRE şi QRF ale produselor secundare Eşi F
obţinute din reacţie
-debitul QRG de produs secundar greu, uleios, G, obţinut din reacţie
- debitul QRP de produs principal
Debitul total QRde ieşire din reactor
intrǎ în schimbǎtorul de caldurǎ S
(în care are loc un proces de rǎcire
datoritǎ fluidului de rǎcire FR)
trecând apoi în decantorul D,
care separǎ produsul greu G,
la baza decantorului rezultând debitul QRG;
produsul greu Gnu poate fi utilizat
şi, în plus, necesitǎ o anumitǎ tratare,
deci implicǎ o anumitǎ cheltuialǎ,
înainte de a fi evacuate ca reziduu,
evitându-se astfel efectele de poluare
provocate de produsul respectiv
Produsele intermediare Eşi F
pot fi folosite numai în calitate de combustibili
şi au, deci, o valoare utilǎ corespunzǎtore;
folosirea debitului (o parte din debitul QPB
de produs de bazǎ al coloanei C)
în calitate de combustibil este posibilǎ
datoritǎ acestei proprietǎţi utile
a produselor intermediare Eşi F
1
PB
Q
Dupǎ separarea produsului greu G, în decantor,
restul componentelor debitului QRintrǎ în coloana C,
din care se obţine la vârf produsul principal P,
(cu debitul QP ),
iar la bazǎ rezultǎ debitul de produs de bazǎ PB,
(cu debitul QPB ),
împǎrţit în debitul de combustibil
şi debitul de recirculare în reactor
Funcţionarea instalaţiei tehnologice
este caracterizatǎ de 12 mǎrimi :
debitele QA, QB, QR, QRA, QRB, QRE, QRF,
QRG, QRP, QP, , şi temperatura Tdin reactor
mǎrimi ce intervin
atât în expresia criteriului de optimizare,
cât şi în expresiile restricţiilor
Optimizarea staţionarǎ are ca obiectiv
sǎ determine pentru instalaţia tehnologicǎ
acel regim de funcţionare de lungǎ duratǎ
(acele valori staţionare optime
ale mǎrimilor caracterizând funcţionarea instalaţiei)
care sǎ asigure maximizarea
ratei anuale de recuperare a investiţiilor
1
PB
Q
2
PB
Q
Expresia criteriului de optimizare
(a ratei anuale de recuperare a investiţiilor)
se obţine în funcţie de unele dintre aceste 12 mǎrimi
Prin intermediul diferitelor grupe de mǎrimi
se obţin expresiile restricţiilor de tip egalitate (9 restricţii)
şi ale restricţiilor de tip inegalitate (14 restricţii)
Restricţiile de tip egalitate
rezultǎ scriind ecuaţiile bilanţurilor de material pe ansamblu
şi pentru fiecare componentǎ
De exemplu, pentru ansambul instalaţiei tehnologice
rezultǎ relaţia de egalitate a debitelor de intrare şi de ieşire
1
PBPRGBA QQQQQ
Optimizarea staţionarǎ
îşi propune ca obiectiv
determinarea valorilor staţionare optime ale unor mǎrimi,
care vor fi menţinute timp îndelungat
pentru obţinerea unui maxim al criteriului,
reprezentat de rata anualǎ de recuperare a investiţiilor,
criteriul însuşi fiind formulat
prin intermediul unui interval mare de timp
Menţinerea anumitor mǎrimi la valorile optime stabilite
se va face prin intermediul unor bucle locale
de reglare automatǎ
în care vor interveni regimuri tranzitorii,
eventuala optimizare a acestor regimuri
formeazǎ obiectul unei alte clase de probleme :
problemele de optimizare dinamicǎ
2. OPTIMIZǍRI DINAMICE
Exemplu: cazul cel mai simplu,
al unui sistem de reglare automatǎ în funcţie de ieşire,
monovariabil (cu o singurǎ mǎrime de referinţǎ w
şi o singurǎ mǎrime reglatǎ), liniar şi continuu
-RA - este regulatorul automat
-EE – elementul de execuţie
-IT – instalaţia tehnologicǎ
-Tr –traductorul
-p1, p2, …, pk – perturbǎrile
care acţioneazǎ asupra instalaţiei tehnologice
Grupând elementele EE, IT şi Tr într-un singur bloc echivalent F,
rezultǎ pentru schema de elemente aspectul
-Feste blocul echivalent al pǎrţii fixate a sistemului
-ε–eroarea
-u– mǎrimea de comandǎ
-y– mǎrimea reglatǎ, respectiv mǎrimea de ieşire a sistemului
Regimul tranzitoriu al mǎrmii reglate y,
provocat de o variaţie treaptǎ unitarǎ a mǎrimii de referinţǎ w
În cazul rǎspunsului la variaţia treaptǎ unitarǎ pozitivǎ
a mǎrimii de referinţǎ w
0
I dt
2
1
0
J dt
0
2
2dtuJ
0
22
213 )( dtuJJJ
Polioptimizare
dtdtdtJ)( 2
2
2
1
0 0 0
2
2
2
14
u)f(x,x
.
x
u
Criteriile integrale pǎtratice
se formuleazǎ şi în limbajul intrare-stare-ieşire,
care descrie funcţionarea sistemului
prin intermediul ecuaţiilor de stare
vectorul stǎrilor
(de regulǎ) vectorul comenzilor,
dar se poate defini şi un set de parametri.
0
5)( dtJRuuQxxTT
0
6)( dtJeRuuyQy TT
x, u, y-vectrii stǎrilor, comenzilor şi mǎrimilor de ieşire (reglate)
Q, Qe, R–matrice de ponderare
FOLOSIREA METODELOR DE DESCOMPUNERE
(DECOMPOZIŢIE)
În ultima vreme a cǎpǎtat o extindere din ce în ce mai
largǎ otehnicǎ de optimizare bazatǎ pe descompunerea
unor probleme de optimizare pentru instalaţii sau sisteme
cu structuri mari şi de complexitate ridicatǎ în probleme de
optimizare ale unor subansambe (subsisteme) mai simple
ale instalaţiilor respective; aceastǎ tehnicǎ este numitǎ
“tehnicǎ de descompunere” sau “tehnicǎ de decompoziţie”
O asemenea descompunere poate prezenta avantaje
importante pentru obţinerea soluţiilor optimizǎrii, întrucât
problemele de optim ale subsistemelor se rezolvǎ mai uşor,
datoritǎ numǎrului redus de variabile care intervin în
fiecare subsistem; pentru obţinerea unui optim pe
ansamblu, pe lângǎ optimizarea funcţionǎrii fiecǎrui
subsistem este necesarǎ şi oacţiune de coordonare centralǎ,
pentru compensarea legǎturilor de interacţiune dintre
subsisteme
Se obţine astfel o structurǎ ierahizatǎ,cu douǎ sau mai
multe nivele; în cazul a douǎ nivele, menţionat anterior, la
nivelul inferior are loc optimizarea funcţionǎrii fiecǎrui
subsistem, iar la nivelul superior se efectueazǎ acţiunea de
coordonare pentru obţinerea optimului de ansamblu
Considerentele referitoare la ierahizarea pe douǎ nivele se
pot generaliza uşor pentru cazul mai multor nivele
Structura unui sistem multinivel cu douǎ nivele este
reprezentatǎ în figura de mai jos
La nivelul inferior se gǎsesc subsistemele 1, 2, 3,
care asigurǎ optime locale,
de exemplu prin maximizarea unor criterii I1, I2, I3,
iar nivelul superior este constituit de coordonatorul 4,
care asigurǎ un opitim global
(schimbând informaţii în ambele sensuri cu subsistemele 1, 2, 3,
informaţii reprezentate simbolic prin sǎgeţi)
de exemplu prin minimizarea unui criteriu I4
Exemplu
Cele douǎ subsisteme formate, subsistemul 1şi subsistemul 2,
sunt legate prin douǎ mǎrimi :
debitul de intrare în coloana C, notat cu
la ieşirea din subsistemul 1
la intrarea în subsistemul 2
fracţiunea recirculatǎ din produsul de bazǎ, notata cu
la ieşirea din subsistemul 2
la intrarea în subsistemul 1
C
Q
^
C
Q
2
PB
Q
2
^
PB
Q
La restricţiile anterioare de tip egalitate se adaugǎ alte restricţii.
Una dintre ele rezultǎ din ecuaţia bilanţului de material
pentru debitul care iese din subsistemul 1:
QC = QR –QRG
aceastǎ ecuaţie conducînd la restricţia
QR –QRG –QC = 0
Alte douǎ restricţii de tip egalitate
rezultǎ din necesitatea egalitǎţii debitelor
de ieşire din unul dintre subsisteme
şi de intrare în celǎlalt
C
CQQ
^
2
22
^
PB
PB QQ
0
^ C
CQQ
^0
22 PBPB QQ
21
^
PBPBP
CQQQQ
0
^
21 C
PBPBPQQQQ
Pentru subsistemul 2rezultǎ osingurǎ ecuaţie de bilanţ de material :
din aceasta obţinându-se restricţia
Stratul 1
conţine algoritmii de reglare automatǎ
şi acţioneazǎ în permanenţǎ pentru ca valorile mǎrimilor reglate,
constituind vectorul Y,
sǎ urmǎreascǎ valorile prescrise prin mǎrimile de referinţǎ,
constituind vectorul W;
stratul 1transmite elementelor de execuţie din blocul F
vectorul mǎrimilor de comandǎ C
Asupra blocului Facţioneazǎ şi vectorul perturbǎrilor P
Stratul 2
include un algoritm de optimizare, care,
pe baza informaţiilor Y şi X
(primite de la blocul F) şi a informaţiilor n
(primite de la mediul înconjurǎtor)
elaboreazǎ valorile optime prescise
pentru mǎrimile reglate Y,
transmiţând stratului 1aceste valori
sub forma vectorului mǎrimilor de referinţǎ W
valorile optime sunt calculate
de algoritmul de optimizare
în ipoteza cǎ parametrii instalaţiei tehologice
(factori de amplificare, constante de timp, timp mort,
constituind un vector A –litera mare pentru α)
nu suferǎ variaţii în timp
stratul 3
realizeazǎ un algoritm de adaptare
Pe baza informaţiilor X, Y şi Z,
primite de la blocul F,
şi a informaţiilor q primite de la mediul înconjurǎtor,
algoritmul 3 determinǎ valorile parametrilor A(α)
şi le transmite stratului 2,
pentru ca algoritmul de optimizare din acest strat
sǎ ia în considerare valorile respective,
uneori fiind modificatǎ în mod corespunzǎtor
structura acestui algoritm
Stratul 3asigurǎ astfel o actualizare
a valorilor parametrilor modelului matematic
al instalaţiei tehnologice
Pe lângǎ cele trei straturi menţionate, poate exista
şi un al patrulea strat, cu algoritm de auto-organizare,
care, în funcţie de anumiţi factori exteriori –de
exemplu factori economici –poate asigura eventuale
modificǎri ale structurii sistemului ierarhizat
În sistemele cu structurǎ multistrat, fiecare strat
intervine la intervale de timp diferite, intervalele fiind cu
atât mai mari cu cât stratul are o poziţie mai înaltǎ
De asemenea, fiecare strat foloseşte un grad de
simplificare diferit pentru modelul instalaţiei tehnologice,
modelul fiind cu atât mai simplificat şi mai abstractizat –
prin evitarea considerǎrii unor detalii –cu cât poziţia
stratului este mai înaltǎ
În unele cazuri, structura multistrat are un aspect
preponderent temporal
Astfel, de exemplu, în cazul optimizǎrii utilizǎrii resurselor
de apǎ, stratul superior determinǎ soluţia pe termen lung, de un
an, folosind un model simplificat, în care toate rezervoarele mai
puţin importante sunt înlocuite printr-un rezervor echivalent
Pe baza calculelor de optimizare din stratul superior, se
transmite stratului mijlociu o soluţie, care nu este însǎ detailatǎ,
întrucât rezervoarele de dimensiuni medii şi mici nu au fost
considerate la volumele lor reale
Stratul mijlociu calculeazǎ mai precis soluţia de
optimizare pe termen de olunǎ, considerând datele
exacte ale rezervoarelor medii şi înlocuind toate
rezervoarele mici printr-un rezervor echivalent,
rezultatul optimizǎrii pentru sfârşitul primei zi a lunii
fiind transmis stratului inferior
În cadrul stratului inferior se calculeazǎ soluţia
detailatǎ pe termen de ozi – ţinând seama de
dimensiunile reale ale rezervoarelor mici –şi se asigurǎ
realizarea acestei soluţii prin intermediul
echipamentelor de reglare automatǎ ale instalaţiilor de
distribuţie a apei.
Pe mǎsura realizǎrii soluţiilor calculate (deci dupǎ fiecare
zi, respectiv dupǎ fiecare lunǎ) se transmit semnale de reacţie
de la stratul inferior la cel mediu şi de la cel mediu la cel
superior, aceste informaţii fiind utilizate ca valori iniţiale – de
cǎtre straturile care le primesc –pentru recalcularea soluţiilor
optime pe intervalele aferente, de olunǎ şi respectiv un an
Descompuneri de tip temporal intervin şi în cadrul
metodelor de optimizare bazate pe principiul maximului
discret şi pe programarea dinamicǎ
