b1280 Technique de L'ingénieur

Optimisation des processus
nergtiques : gestion optimale

par
Christian CHATELAIN
Ingnieur de lcole Nationale Suprieure dArts et Mtiers
Ingnieur-Chercheur la Direction des tudes et Recherches de lectricit de France (EDF)
Jean-Claude DUCROCQ
Ingnieur Civil des Mines
Benot MIGNARD
Ingnieur Civil des Mines
et
Michel COEYTAUX
Ingnieur de lcole Centrale des Arts et Manufactures
Chef de Produit SECI Manager et salles de contrle de la socit SERETE
1.
1.1
1.2
1.3
Analyse des critres conomiques .....................................................

Critres de choix dun procd...................................................................
Optimisation dun processus de fabrication .............................................
Problme des donnes................................................................................
2.
Constitution et gestion optimises des rseaux dnergie

dans les industries...................................................................................
Principes de la modlisation.......................................................................
Formulation mathmatique ........................................................................
Conclusion....................................................................................................
8
9
10
14
3.1
3.2
3.3
3.4
Optimisation des schmas de transferts dnergie dans les

procds industriels................................................................................
Position du problme ..................................................................................
Une mthode doptimisation......................................................................
Exemple dun procd de la chimie...........................................................
Conclusion....................................................................................................
16
16
16
18
18
4.
4.1
4.2
4.3
4.4
4.5
4.6
Gestion de lnergie en temps rel .....................................................

Dfinitions ....................................................................................................
Fonctions gnrales de gestion et leur hirarchie....................................
Structure des systmes de contrle-commande ......................................
Logiciels........................................................................................................
Rentabilit dune gestion optimale en temps rel ....................................
volutions prvisibles .................................................................................
19
19
20
22
23
23
24
Rfrences bibliographiques .........................................................................
24
2.1
2.2
2.3
B 1 282
2 - 1989
3.
B 1 282 - 2
article Optimisation des processus nergtiques se compose de trois

articles :
Concepts [B 1 280] ;
Applications [B 1 281] ;
Gestion optimale [B 1 282].
Largement mise en uvre dans la plupart des procds industriels, lnergie
est un bien prcieux quil faut conomiser, mais non sans comptabiliser
paralllement les investissements et les frais qui sattachent aux oprations dont
elle est lobjet.
Toute reproduction sans autorisation du Centre franais dexploitation du droit de copie est strictement interdite.
Techniques de lIngnieur, trait Gnie nergtique
B 1 282 1
OPTIMISATION DES PROCESSUS NERGTIQUES : GESTION OPTIMALE ___________________________________________________________________________
Ainsi la gestion de lnergie se prsente-t-elle comme un problme doptimisation technico-conomique qui se pose diffrents niveaux :
dans le choix de lquipement qui permettra dassurer une mme fonction
globale de la faon la moins coteuse ;
dans le mode dexploitation le plus judicieux de cet quipement, souvent
indissociable de son choix ;
dans le meilleur compromis performance-prix attach chaque appareil
participant la fonction ;
enfin, dans le contrle instantan du fonctionnement de lensemble pour le
maintenir en permanence au plus prs de ses conditions optimales.
Les spcialistes qui ont trait, en collaboration, ce vaste sujet lont dcoup
en quatre paragraphes :
analyse des critres conomiques ;
constitution et gestion optimises des rseaux dnergie dans les
industries ;
optimisation des schmas de transferts dnergie dans les procds
industriels ;
gestion de lnergie en temps rel ;
o se trouvent approfondis les diffrents aspects que lon vient dvoquer.
1. Analyse des critres

conomiques
Nota : le lecteur se reportera utilement la rubrique Management de lentreprise du
trait Lentreprise industrielle, et en particulier larticle Choix des investissements
[A 4 450].
1.1 Critres de choix dun procd

Le comportement des entreprises est principalement ax sur le
futur, puisquil concerne les projets dinvestissements et les besoins
de capitaux. La technologie offrant, presque toujours, une diversit
de procds pour une production donne, un choix doit tre fait entre
ces diffrents procds. La plupart des dcisions qui en dcoulent
ncessitent un arbitrage entre le prsent et lavenir.
De deux projets, lindustriel choisira le plus rentable. Encore faut-il
sentendre sur ce terme ! Le cot conomique repose sur deux
facteurs : lun est dfinitif, il sagit de linvestissement et de son
financement pour lequel lchelonnement des dpenses est bien
connu ; lautre, le cot dexploitation, ne peut tre questim selon
les caractristiques techniques probables et les conditions conomiques extrapoles partir de la situation actuelle.
Le problme (rendre un procd le plus conomique possible) bute
sur un premier obstacle : comment concilier cot dinvestissement
et cots dexploitation ? En effet, linvestissement est pay au tout
dbut de la mise en place du procd, tandis que les frais dexploitation sont dbourser pendant toute la dure de vie de linstallation.
Il faudrait mme, en toute logique, intgrer, ces cots, les frais de
dclassement et de remplacement.
Le choix peut se faire par le calcul de la valeur actuelle (ou bilan
actualis) de chacun des projets caractriss par linvestissement et
par les dpenses et recettes prvisibles pour chaque anne venir.
au prix que le secteur productif est prt payer pour disposer de

moyens de financement supplmentaires et, rciproquement, au
prix que la collectivit est prte accepter pour renoncer une
satisfaction immdiate.
La majorit des individus prfrent disposer dun bien immdiatement au lieu den disposer seulement dans un ou deux ans. Cette
prfrence nest pas uniquement fonde sur lincertitude de lavenir
ou sur linflation, mais aussi sur la croissance conomique. Dune
promesse ferme de recevoir 1 000 francs dans cinq ans, on peut
prfrer les recevoir tout de suite, afin de faire travailler cet argent
pour en tirer profit.
Plus la croissance est rapide, plus les besoins dinvestissement
sont grands, et plus le taux dactualisation est lev.
Pour des raisons de simplification, nous supposerons ce taux
dactualisation a constant dans le temps, et les cots seront exprims
en francs de lanne 0. De plus, linvestissement est entirement
imput sur cette premire anne, les cots proportionnels
commenant lanne 1.
De ce fait, on a la relation dquivalence suivante :
1
1 franc de lanne n = ----------------------franc de lanne 0
( 1 + a )n
1.1.2 Valeur du taux dactualisation choisir

Si le dcideur appartient une entreprise prive, lestimation de
ce taux est, en gnral, subjectif.
Dans le cas o le financement se fait uniquement par un emprunt,
le taux sera de lordre de :
a = i (1 b) + s
avec
b
i
1.1.1 Actualisation des cots futurs

Pour exprimer ces dpenses dans une mme unit, il est ncessaire de disposer dun indicateur global, le taux dactualisation ; il
mesure la raret des ressources de financement que la collectivit
peut mettre la disposition du secteur productif [1] : cela correspond
B 1 282 2
t a u x d i m p o s i t i o n d e s b n fi c e s i n d u s t r i e l s e t
commerciaux,
taux dintrt de lemprunt,
s
taux de risque ou coefficient de scurit.
Dans le cas o une part dautofinancement intervient, le calcul de
ce taux dactualisation est encore plus subjectif et arbitraire.
__________________________________________________________________________ OPTIMISATION DES PROCESSUS NERGTIQUES : GESTION OPTIMALE
Entre deux projets dont lun procure un bnfice actualis certain

et lautre le mme bnfice mais moins probable, lentrepreneur
prfre gnralement le premier. Il existe deux possibilits de tenir
compte de cela :
soit alourdir les projets alatoires dun montant de scurit ;
soit augmenter le taux dactualisation de ces projets
(+ 5 20 %).
Si le dcideur appartient un service public, il na aucun choix : il
prend le taux dactualisation fix par le Commissariat au Plan 9 %.
En rgle gnrale, ces bilans sont calculs en francs constants.
Toutefois, si lon dsire constituer un chancier en francs courants, il faut alors intgrer linflation dans le taux dactualisation :
r est donc une fonction monotone dcroissante du temps de rcupration. Minimiser le temps de rcupration revient donc maximiser le taux de rentabilit interne.
Par contre, bnfice actualis et temps de rcupration ne sont
pas toujours quivalents comme le montre lexemple suivant : (0)
dure de vie T ......................(ans)

investissement I .................... (kF)
bnfice brut.....................(kF/an)
temps de rcupration.......(ans)
taux de rentabilit r................(%)
bnfice actualis B............. (kF)
Projet I
Projet II
10
1 000
500
2,0
49,1
2 209
10
2 000
750
2,7
35,7
2 813
choix du projet I
choix du projet I
choix du projet II
1 + a = (1 + a ) (1 + i )
Cependant, la difficult dvaluer linflation future ne plaide
pas en faveur de cette mthode qui, de toute faon, napporte
aucune indication supplmentaire sur le choix faire.
Les calculs ultrieurs sont donc toujours raliss en francs
constants.
Remarquons ici que francs constants ne veut pas dire prix
constants. Des drives de prix peuvent tre incluses pour certains
biens (lnergie, par exemple).
1.1.3 Critre conomique utiliser

Le bilan actualis B (ou bnfice net actualis), le discounted
cashflow des Amricains, est la diffrence entre la somme des
recettes Rn et la somme des dpenses Dn pendant la dure de vie T
de linstallation :
T
B = I+
n=1
( Rn Dn )
VT
--------------------------+ ---------------------( 1 + a )n
( 1 + a )T
(n = 1 T )
(1)
avec
I
investissement,
VT valeur rsiduelle de linvestissement.
Dautres critres de choix existent, mais prsentent des dfauts,
comme nous le verrons au paragraphe 1.2.2.
Le taux de rentabilit interne r est gal au taux dactualisation qui
annule le bnfice actualis. Une opration est rentable si son taux
de rentabilit dpasse le taux dactualisation.
Le temps de rcupration est le temps ncessaire pour annuler le
bilan prcdent (en gnral non actualis) :
Investissement
Temps de rcupration = -----------------------------------------------Bnfice annuel
(2)
Son usage revient privilgier les investissements, cest--dire

rcuprer la mise de fonds le plus rapidement possible.
Le temps de retour sapplique surtout la comparaison de deux
projets entre eux :
Investissement
Temps de retour = ----------------------------------------------------------- Cot dexploitation
(3)
Les deuxime et troisime critres sont quivalents si les

recettes R et les dpenses D sont invariantes au fil des ans et si la
valeur rsiduelle est nglige. En effet, dans ce cas, le taux de rentabilit interne r se dtermine par lquation :
1 1 / ( 1 + r )T
I
------------------------------------ = ---------------r
RD
(4)
Si les dures de vie des quipements sont diffrentes, il faut choisir

comme horizon de calcul le plus petit commun multiple des dures
de vie de chacun des projets, en procdant des renouvellements
lidentique. On comparera, par exemple, trois oprations successives de centrales thermiques (dure de vie de 30 ans) et une opration de centrale hydraulique (dure de vie estime 90 ans).
On peut aussi prendre une dure T quelconque (la plus petite dure
de vie par exemple) et faire les comparaisons en tenant compte de
la valeur rsiduelle des quipements la fin de cette priode.
1.2 Optimisation dun processus

de fabrication
Choisir entre deux technologies diffrentes suppose que chacun
de ces deux projets est optimal linstant de la dcision. Pour que
ce choix soit cohrent, il est indispensable que les critres qui rgissent loptimisation de chaque procd soient identiques non seulement entre eux, mais aussi avec le critre de dcision finale. Bien
que cette condition paraisse vidente, elle nest pas toujours facile
respecter.
Pour des raisons dopportunit commerciale, les fournisseurs
peuvent, dans certains cas, avoir tendance privilgier le court terme
(investissement minimal) au dtriment des dpenses futures
(nergie par exemple). Mais ce nest pas la seule raison dun manque
doptimalit de certaines installations.
Les deux principales causes sont le cot supplmentaire de ltude
doptimisation et la dispersion des connaissances et des responsabilits. Lensemblier na pas toujours la possibilit de matriser la
fonction investissement des lments du procd quil sous-traite ;
et il ne peut multiplier ses demandes de devis linfini. Quant au
sous-traitant, il est rarement inform de lensemble dans lequel son
appareil simbrique. De plus, la fourniture des utilits (nergie)
nentre gnralement pas dans la limite des fournitures de lappel
doffre, ce qui fait que celles-ci sont rarement intgres dans la fonction conomique optimiser. On pourra voir linfluence de cette
dernire constatation dans les exemples traits au paragraphe 1.2.2.
1.2.1 Multitude de variantes conomiques

pour une solution
Tout acte industriel met en uvre des matriels (la machine), de
lnergie (ce qui la fait fonctionner) et du contrle-commande
(llment de dcision : main-duvre ou automatisme).
Avant daboutir cette machine industrielle, il a fallu tudier les
phnomnes mis en uvre, les tester en laboratoire afin de valider
sa faisabilit, puis extrapoler ces rsultats lchelle de la production
industrielle, souvent par lintermdiaire dun pilote semi-industriel.
B 1 282 3
Toutes ces tapes nous confortent sur le bon fonctionnement technique de notre machine ; mais quen est-il de sa performance
conomique ?
Une fois le schma du procd retenu, il nous reste dterminer
ses paramtres techniques ; ce sont, par exemple, une vitesse de
rotation, un pincement dchangeur, une pression de fonctionnement, une paisseur de calorifuge, etc.
Presque tous les processus lmentaires de fabrication ralisent
un transfert et peuvent se rduire une fonction de la forme :
F = kSEspcif
(5)
dans laquelle :
Espcif est la variation du potentiel de transformation (consommation dnergie unitaire) ncessaire pour raliser ce transfert,
F
est le flux, cest--dire la production horaire de la machine,
k
est un coefficient technique qui dpend de la technologie
employe,
S
est une surface de transformation qui dtermine la taille de
lappareil et donc linvestissement consommateur
dnergie.
Les exemples de tels processus sont nombreux. Citons-en deux.
Alimentation dun appareil de puissance lectrique F (kW) sous une
tension continue U travers un cble lectrique de conductivit et de
longueur L. Si S est la section du cble, les pertes en ligne sont :
L
W ( kW ) = 2 ------ ( F/U ) 2
S
Lnergie ncessaire au processus nest gnralement pas de

lnergie marchande comme le fuel lourd ou llectricit haute
tension alimentant lusine. Il faut la transformer par lintermdiaire
de machines (chaudire, poste et transformateur lectrique) en
nergie utilisable (vapeur deau, lectricit MT ou BT, etc.), puis la
transporter dans un rseau interne jusquaux machines utilisatrices.
Cet investissement producteur dnergie ne doit pas tre oubli. Il
varie dans le mme sens que la consommation dnergie, et donc
en sens inverse de linvestissement du procd. Linvestissement
total prsente, dans la plupart des cas, un minimum pour une valeur
donne du paramtre nergtique.
Bien entendu, les quations relles utiliser sont parfois plus
compliques ; mais en fin de compte, il est souvent possible, dans
une plage de variation dtermine, dapproximer celles-ci par des
fonctions de type puissance [2] :
investissement producteur dnergie :
Ipro = A + B E
investissement consommateur dnergie :
Icon = C + D E
cot annuel de lnergie :
C ann = F E
avec = 0,8 1 et = 0,25 1 ;
A, B, C, D, F coefficients fonction du procd envisag ;
E consommation dnergie.
Le cot total sur la priode T est donn par :
F / U tant le courant consomm.

Cette quation peut se mettre sous la forme :
U2
F = ---------- S ( W/F )
2L
avec E = W /F consommation spcifique.
Un module dosmose inverse (par exemple, en concentration de
lactosrum) se compose de membranes de surface S (m2 ), dune
pompe de dbit Ft (m3/ h) crant une surpression de p (bar). Si la pression osmotique est faible vis--vis de la pression p, le dbit de permat
(eau) est F (m3/ h) = hS p, avec h coefficient de dbit de la membrane.
Lnergie dpense la pompe est :
C total = I pro + I con + C ann

avec
a
T
F = Ft (1 C i / C f )
Ces quations scrivent aussi :
Le coefficient technique k dpend souvent du potentiel nergtique et introduit un lment de non-linarit dans les quations. Par
exemple, le coefficient dchange thermique diminue lorsque lcart
de temprature diminue.
B 1 282 4
(6)
Prenons, comme exemple, un procd de concentration :

linvestissement producteur dnergie (chaudire fuel, par
exemple) est gal :
Ipro = 600 E 0,88
linvestissement consommateur dnergie (cas dun vaporateur
multiples effets) est gal :
Icon = 287 700 E 0,75
le cot annuel de lnergie (fuel lourd) est gal :
Cann = 120,3 E (en kF)
E est la consommation de fluide nergtique et est exprime en
tonnes de vapeur deau par heure.
Le cot total sur la priode T est :
h
F = ( 1 C i /C f ) ------------- S ( W/F )
0,04
Dans tous les cas industriels qui nous intressent, la production
est fixe par le cahier des charges.
Daprs la relation (5) et pour une production donne, la
consommation spcifique dnergie du procd est inversement
proportionnelle au paramtre fixant la taille de lquipement. De ce
fait, linvestissement de lappareil consommateur dnergie (facteur
dchelle ngatif) varie en sens inverse de la consommation
dnergie. On retrouve ici le principe bien connu : conomiser
lnergie cote cher en investissement.
taux dactualisation annuel,

dure de vie des quipements.
W (en kWh/ h) 0,04 Ft p

Soit C i la concentration initiale du produit et C f sa concentration finale,
alors les dbits F t et F sont lis par la relation :
1
1 ---------------------( 1 + a )T
---------------------------------a
C total =
600E 0,88
+ 287 700 E
0,75
+ 120,3 E
1
1 ---------------------T
(1 + a)
-------------------------------a
1.2.2 Principe doptimisation

Comme nous avons un degr de libert (lnergie dans notre
exemple) pour notre processus, il est judicieux de choisir la meilleure
valeur de ce paramtre. Ce dernier nest gnralement pas directement lnergie consomme E, mais une caractristique technique S
de lappareil lie cette consommation par la relation (5).
Le premier critre est de rendre minimal le cot total F ( E )

[relation (6)], ce qui revient annuler la drive de la fonction :
F (E )/ E avec 2F (E )/ E 2 positif. La relation (5) conduit une relation E = G ( S ). Le critre F ( E ) devient F ( G ( S ) ) que nous
notons H (S ). La minimisation de ce nouveau critre sexprime par :
F ( E ) G ( S )
H ( S )
--------------------- = 0
soit -------------------- ---------------------- = 0
E
S
S
F ( E ) 2 G ( S ) 2 F ( E )
G ( S )
2H (S)
- > 0 soit ------------------- ------------------------ --------------------et ---------------------+ --------------------E
S
S 2
S 2
E 2
Tableau rcapitulatif des rsultats obtenus

(base 100 = cots au point A)
Critre
>0
Deux cas qui respectent ces quations peuvent se prsenter :

F (E )/ E est nul et le critre est vrifi ;
G (S )/ S est nul avec 2G (S )/ S 2 de mme signe que
F (E )/ E.
Ce deuxime cas a trs peu de chances de se produire, car il signifie
quil existe un extrmum de la consommation dnergie pour une
valeur donne du paramtre technique.
On peut donc dire quoptimiser le critre selon un paramtre
technique ou selon lnergie, comme dans notre exemple, aboutit
raliser un mme quipement.
nergie
Temps
Investisse- Cot
de
ment
total rcupration
(ans)
Cot total
Cot hors utilits
A
B
100
125
100
93
100
102
1,86
1,77
Investissement
200
86
119
1,81
Temps de rcupration
avec utilits
D
sans utilits
E
150
415
89
103
106
192
1,75
2,99
avec utilits
sans utilits
D
E
138
311
90
92
104
155
2,47
3,58
avec utilits
sans utilits
D
E
158
519
88
115
108
230
1,35
2,67
Nous ne dvelopperons pas ici le cas de plusieurs degrs de

libert ; cela peut arriver avec des paramtres techniques
secondaires (par exemple, la pression daspiration au compresseur
pour une vaporation avec compression mcanique de vapeur est
un paramtre secondaire par rapport au taux de compression).
Linfluence de ces paramtres, souvent du second ordre, doit
toutefois tre envisage lors dune tude dun projet rel. Cela peut
aussi se produire si le systme tudi comporte plusieurs processus
lmentaires lis entre eux. Le principe appliqu reste cependant le
mme, bien que sa rsolution soit moins aise.
Dautres critres peuvent galement tre pris en compte.
Minimisation du cot total (figure 1a ) ; lorsque lon optimise un
processus, les recettes ne sont pas lies la technologie employe,
mais la vente du produit fini. Nous pouvons donc remplacer la
maximisation du bnfice actualis par la minimisation du cot total,
somme de linvestissement total et du cot actualis dexploitation
[relation (6)].
Lnergie consomme (point A) est de 20 t /h pour un cot total
de 54 200 kF. Si les utilits (investissement producteur dnergie)
sont oublies (point B), la consommation nergtique passe
25 t /h, soit une augmentation de 25 % avec un cot total rel
(point B) de 55 200 kF (+ 2 %) et un investissement rduit de 7 %.
Minimisation de linvestissement total (figure 1a) : si ce critre
permet de rduire linvestissement de 14 % (point C), il augmente le
cot actualis de 19 % (point C ) et la consommation dnergie
de 100 % !
En comparant ces deux projets (A et C), on peut constater
(figure 1b ) que la solution investissement minimal devient plus
onreuse ds la deuxime anne de fonctionnement ; ce second
critre est donc carter.
Minimisation du temps de rcupration (figure 2) : comme nous
lavons vu au paragraphe 1.1.3, cela revient galement maximiser
le taux de rentabilit interne. Ce critre pose une relle difficult
pour estimer correctement le montant de la recette annuelle, lequel
dpend surtout du march et du circuit commercial du produit
fabriqu.
titre dexemple, on peut calculer les optimums lorsque cette
recette varie de + 25 % (points D et E) ou de 25 % (points D
et E). Si le temps de rcupration optimal est suprieur un an,
alors la consommation dnergie est suprieure celle du point A
(+ 50 % pour le point D). Avec ce critre, ne pas tenir compte de
linvestissement des utilits donne un optimum farfelu (point E)
puisquil cote plus cher en investissement et en exploitation que
la solution correspondant linvestissement minimal (point C).
Figure 1 Optimisation du cot dun processus
Ce critre est donc, lui aussi, carter des processus doptimisation et de dcision finale. Seul le cot total, ou le bnfice actualis,
est utilis par la suite.
(0)
B 1 282 5
Figure 2 Optimisation du temps de rcupration

Figure 3 Influence du taux dactualisation sur le cot total
Sur la relation (6), on constate que les paramtres fixant la

valeur de loptimum sont au nombre de trois : le taux dactualisation, la dure de vie des quipements et la variation relative des
prix de lnergie par rapport au cot dinvestissement.
1.2.3 Influence du taux dactualisation

Entre un processus optimis un taux de 9 %, et celui optimis
pour un industriel peu confiant dans lavenir (a = 25 %), la
consommation dnergie du second est augmente de 23 %
(figure 3).
linverse, dans un pays conomie dirige o le taux
dactualisation serait nul, le processus consommerait 17 % dnergie en moins.
1.2.4 Influence de la dure de vie des quipements

Ce que lon entend ici nest pas la dure de la vie physique de la
machine, mais sa dure de fonctionnement conomique. En effet,
une centrale lectrique, par exemple, est prvue pour fonctionner
de nombreuses annes car lvolution de la demande future dlectricit est bien cerne. Par contre, dans certaines branches industrielles,
lvolution du march et des techniques est telle que lobsolescence
du produit ou de la machine est rapidement atteinte.
Cette dure de fonctionnement est donc trs lie au type
dindustrie qui utilise lquipement. Il nest, ainsi, pas anormal de
trouver des quipements optimiss ayant des consommations trs
diffrentes si les conditions dutilisation sont elles-mmes diffrentes (figure 4).
lextrme, entre le secteur public (a = 9 % et T = 30 ans) et une
industrie en volution trs rapide et ne prenant pas de risques
(a = 25 % et T = 3 ans), la diffrence de consommation nergtique
dun mme processus est de lordre de 80 %.
B 1 282 6
Figure 4 Influence de la dure de vie de linstallation sur le cot total
1.2.5 Influence de lenchrissement de lnergie

Le troisime lment modifiant loptimum conomique est la
variation du cot des nergies vis--vis du cot dinvestissement.
Il peut tre regard de deux manires.
Dune faon statique : pour un mme dcideur et un mme
processus, loptimum actuel est-il identique celui dil y a dix ans ?
Si lon considre que les cots de notre exemple sont ceux de 1982,
la figure 5 montre lvolution du processus optimal entre 1973
et 1984. Entre ces deux dates, le rapport prix du fuel lourd sur cot
dinvestissement a t multipli par prs de 4. De ce fait, le processus
optimis en 1984 consomme 36 % dnergie de moins que celui
davant le premier choc ptrolier.
Ceci montre bien quil ne suffit pas de reprendre des tudes antrieures pour un nouveau projet et que la consommation dnergie
dun procd nest pas une constante universelle dpendant uniquement de la technologie mise en uvre.
Dune faon dynamique, cest--dire en incluant une drive d sur
le prix des nergies dans la relation (6) du cot total :

1+d
C total = I pro + I con + C ann 1 --------------1+a
1+d
-------------ad
T
(7)
Ici, nous navons quun combustible : le fuel lourd. Dans un cas

rel, plusieurs nergies sont ncessaires (lectricit et combustible,
par exemple), ce qui fait autant de coefficients de drive de prix.
Une drive a moins dinfluence sur loptimum que la modification
du prix de lnergie car elle ne sapplique progressivement que sur
les annes futures, elles-mmes pondres par lactualisation. Cela
ne veut pas dire quil ne faut pas en tenir compte, comme le montre
la figure 6. En partant des cots de 1982, on applique une drive
sur le fuel variant entre 12,5 % (retour au prix de 1973) + 13,4 %
(prix multipli par 4 au bout de 11 ans). La consommation optimale
dnergie est alors comprise entre + 23 % et 25 % de la valeur
calcule sans drive du combustible.
Avoir une incertitude sur le prix futur de lnergie est parfaitement comprhensible, il ne convient pas pour autant de majorer le
taux dactualisation. En effet, cette incertitude naffecte quune partie du cot total (le cot annuel de lnergie en question).
Dans le tableau suivant, on calcule, partir de lexemple prcdent,
les taux dactualisation avec une drive nulle [relation (6)] qui
donnent le mme optimum nergtique quune drive de prix avec
un taux dactualisation de 9 % [relation (7)]. Rsultat : si lnergie a
tendance augmenter, alors il faut diminuer le taux dactualisation,
cest--dire prendre un taux de risque ngatif !
(0)
nergie loptimum
(t/h)
a (pour d = 0 %)
(%)
16,6
18,1
20,0
21,9
23,7
25,4
27,1
0
3,8
9,0
14,7
21,1
28,2
36,3
d (pour a = 9 %)
(%)
9
5
0
5
10
15
20
1.2.6 Limites des paramtres

Il peut arriver que la valeur optimale dun paramtre corresponde
une impossibilit technique. Cest le cas par exemple de lcart
de temprature sur un vaporateur flot tombant, qui cre des
instabilits de fonctionnement au-dessous de 1 ou 2 oC. Ces limites
voluent avec le progrs technique.
Il arrive aussi que la limite de scurit soit fixe par une norme,
comme par exemple pour la densit maximale de courant admissible dans un cble lectrique.
1.2.7 Discontinuits de linvestissement

Presque tous les matriels sont choisir dans le catalogue dun
constructeur. Cela veut dire que le paramtre technique varie par
paliers. Cest le cas des moteurs, des pompes, des cbles lectriques,
des diamtres de tuyauterie, des chaudires, etc.
Dans certains cas, il y a mme passage dune technologie bon march une technologie plus sophistique. Au-del dun taux de
compression donn, il faut, par exemple, utiliser un compresseur
deux tages qui cote environ 50 % plus cher.
Figure 5 Influence du cot des nergies
Si tous les appareils du processus sont de ce type, loptimisation

se fait sur une variable discrtise. Sil reste un appareil dont la taille
na pas ce type de contrainte (surface dchange par exemple), la
recherche de loptimum doit se faire dans chaque gamme de matriels (figure 7).
B 1 282 7
1.3 Problme des donnes

1.3.1 Cots dinvestissement
tant donn toutes les approximations faites sur le taux
dactualisation, la dure dutilisation, la drive du prix de lnergie,
mais aussi lallure assez plate du minimum de la courbe de cot total,
il nest pas indispensable de rechercher une trop grande prcision
pour ltablissement du cot dinvestissement. Des formules de type
puissance sont gnralement suffisantes : I = A + BE ( 1.2.1), par
exemple. Nombre de publications fournissent ces donnes. Pour les
matriels du gnie chimique, citons le manuel dvaluation conomique des procds [3] et, pour les matriels lectriques, le recueil
dAmpre 82 [4].
Utiliser ce type de formules apporte galement lassurance de ne
pas trop grever le cot dtude du projet.
1.3.2 Cots proportionnels

Par mesure de simplification, le cot dentretien est souvent pris
gal un certain pourcentage de linvestissement (3 6 % selon les
matriels).
Le cot des nergies est le cot actuel, auquel on peut appliquer
plusieurs scnarios de drives. Pour les nergies dont les prix sont
saisonnaliss (tarifs), il suffit de prendre le cot moyen marginal car
la variation du paramtre technique ne modifie pas la rpartition
dnergie dans les diffrentes tranches de tarification. Dans ce cas,
il faut quand mme vrifier que les contraintes du tarif restent
respectes (puissance souscrite en heures creuses infrieure ou
gale la puissance en heures pleines par exemple).
Figure 6 Influence dune drive d du prix des combustibles (fuel)
2. Constitution et gestion
optimises des rseaux
dnergie dans les industries
Figure 7 Effets de la discontinuit des investissements.

Cas de deux optimums I et II
B 1 282 8
Le choix dun investissement nergtique et la gestion optimale

de lnergie sur un site industriel doivent tre envisags en prenant
en compte deux dimensions importantes du problme :
linsertion de lquipement nergtique au sein dun systme
dquipements interdpendants : cest la dimension systmique du
problme rsoudre ;
les variations au cours dune anne dexploitation des tarifs des
nergies, des besoins des procds de fabrication en utilits (vapeur,
nergie mcanique, etc.), des disponibilits des quipements : cest
la dimension temporelle.
Illustrons qualitativement ces deux points sur le systme
chaleur-force de la figure 8. On pourrait envisager linstallation
dune compression mcanique de vapeur (CMV) substituant une
puissance lectrique 2 MW 20 t/h de vapeur basse pression dans
un atelier du process, induisant ainsi des modifications sensibles
au niveau du systme chaleur-force. Le tonnage de vapeur conomis conduit, en effet, une baisse de lautoproduction qui peut
tre du mme ordre de grandeur que la puissance de la CMV
elle-mme.
valuer la rentabilit de la compression de vapeur lchelle de
latelier concern conduit donc une erreur de diagnostic par
non-prise en compte des interdpendances entre quipements lies
laspect systme. De la mme faon, une modification des
diffrents besoins des procds en vapeur haute, moyenne et basse
Figure 8 Exemple de systme chaleur-force, auquel on envisage dinstaller une compression mcanique de vapeur (CMV)
pression conduit une modification de lallocation de charge entre

les turbines, les dtentes statiques et les chaudires du systme
chaleur-force si lon dsire minimiser le cot de lensemble des
approvisionnements nergtiques du site. Une variation des prix
relatifs des combustibles sous chaudires et de llectricit du rseau
conduit dailleurs la mme remarque.
Lexemple de la figure 8 montre la ncessit dun outil mthodologique et de son support informatique permettant :
loptimisation des choix dinvestissements nergtiques ;
la gestion optimale de lnergie, moyennant la prise en compte
de lensemble du systme nergtique et des diffrentes variations
pouvant intervenir au cours dune anne dexploitation.
2.1 Principes de la modlisation

Un systme nergtique industriel et son exploitation peuvent tre
caractriss par les lments suivants :
un ensemble dquipements producteurs ou transformateurs
dnergie dont le but final est de fournir tout instant les utilits
(vapeur diffrents niveaux thermiques, eau chaude, lectricit, etc.)
ncessaires au fonctionnement du procd de fabrication ;
la topologie du systme chaleur-force, cest--dire la faon dont
les quipements sont relis entre eux ;
le fonctionnement du procd, gnralement fix par des

impratifs de production et par des contraintes techniques ; il nest
donc gnralement pas sujet optimisation au cours du temps, sauf
en cas dinstallation dun nouvel quipement ; pour cette raison, le
procd peut tre reprsent par la description de ses besoins en
utilits au cours du temps ; un investissement au niveau du process
est alors trait au travers de la modification induite des besoins du
procd en utilits ;
les sources dnergie disponibles, ainsi que leurs tarifications
respectives.
Le choix optimal des investissements puis la gestion de lnergie
peuvent alors tre poss dans les termes suivants : minimiser
lensemble des cots nergtiques et des charges lies aux
investissements en respectant :
les caractristiques de fonctionnement des diffrents quipements, cest--dire les relations entre leurs entrants et leurs sortants ;
la topologie du rseau nergtique ; cela peut tre obtenu en
crivant que la somme des productions de chaque produit doit tre
au moins gale la somme de ses utilisations ;
les capacits maximales et les minimums techniques des
diffrents quipements ;
les besoins du procd en utilits ;
dventuelles contraintes techniques propres au systme
tudi.
B 1 282 9
La minimisation de la fonction conomique propose doit tre

effectue sur une anne dexploitation complte et les contraintes
prcdemment listes doivent tre vrifies sur chacune des
priodes reprsentant :
les variations dactivit du procd au cours du temps ;
les variations des cots des nergies en raison de lexistence
de structures tarifaires de certaines nergies (lectricit, par
exemple).
Pour raliser cet objectif, on peut utiliser les programmations
linaire et mixte (article Optimisation des processus nergtiques :
concepts [B 1 280], dans le prsent trait). Les contraintes cites sont
crites sous une forme linaire ou linarise par morceaux, puis
dupliques sur chacune des N priodes de la modlisation. ces
contraintes monopriodiques peuvent venir sajouter dautres
contraintes liant des priodes entre elles (phnomnes de stockage,
prime fixe de certaines tarifications nergtiques). On remarquera
que les contraintes descriptives de la topologie du systme
permettent lintgration de la dimension systmique, tandis que la
duplication des contraintes sur un ensemble de N priodes descriptives dune anne permet la prise en compte de la dimension temporelle du problme.
Le recours la programmation linaire pour traiter ce type de
problme se justifie notamment par :
la certitude de parvenir loptimum ;
le nombre lev de variables mises en jeu, ce qui poserait de
srieux problmes pour les techniques doptimisation non linaire
(plusieurs milliers de variables et de contraintes dans la pratique).
2.2.1 Modlisation du groupe turboalternateur

La courbe caractristique du groupe turboalternateur, relation
entre la puissance Z sur larbre et le dbit de vapeur T2 (figure 10),
peut tre modlise par les relations suivantes :
1
T 2 = b 1 t 1 T 2 + b 2 t 2 T 2
Z = b 1 t 1 Z 1 + b 2 t 2 Z 2
b1 t1
avec
t2 b2
b2 b1 0
t 1 1, t 2 1
o b1 et b 2 sont des variables boolennes et t 1 et t 2 des variables
continues.
1
Les constantes T 2 et T 2 dsignent les longueurs (en t/ h) des

deux intervalles de dbit dont on prend les fractions t1 et t 2 .
Les variables boolennes b1 et b 2 ont pour rle de nautoriser
2
lutilisation de lintervalle T 2 quune fois lintervalle T 2

compltement utilis (t1 = 1). Dans ce cas, en effet, b1 = b 2 = 1, ce
qui implique t1 = 1. Puisque b 1 = 1 t 1 et t 1 1 , il en rsulte alors :
1
T 2 = T 2 + t 2 T 2
avec 0 t 2 1 (on a de mme Z = Z 1 + t 2 Z 2).
Si b1 = 1 et b2 = 0, t 2 b 2 = 0 conduit t 2 = 0, soit :
2.2 Formulation mathmatique

On considre le systme nergtique de la figure 9 compos
initialement dune chaudire haute pression, dun groupe
turboalternateur, dune dtente statique et dun ballon de stockage
de la vapeur basse pression. On envisage linvestissement dune
chaudire de capacit inconnue X 6 . Lanne dexploitation est reprsente par deux jours types comportant chacun trois configurations
dexploitation pour lesquelles les besoins du process en vapeur Q v

et en lectricit
Qe
Par contre, la contrainte b 2 b 1 0 interdit quant elle des solutions du type t1 = 0 et t 2 0. Dans ce cas en effet, t 2 b 2 conduit
b2 = 1, et b 2 b 1 0 induit b 1 = 1 t 1 , soit t 1 = 1 ce qui contredit
t 1 = 0.
ont t indiqus.
La version tarifaire de llectricit comporte cinq postes qui

regroupent les priodes dactivit du systme de la faon suivante :
heures de pointe
: priode 1
heures pleines dhiver : priode 2
heures creuses dhiver : priode 3
heures pleines dt
: priodes 4 et 5
heures creuses dt
: priode 6
Avant de dfinir la fonction conomique minimiser et de lister
lensemble des contraintes du problme, deux points mritent des
explications complmentaires : la modlisation du fonctionnement
du groupe turboalternateur et la modlisation du tarif EDF. Les
autres contraintes numres au paragraphe 2.2.3 sont beaucoup
plus simples et ne ncessitent pas de commentaires particuliers.
B 1 282 10
T 2 = t 1 T 2
2.2.2 Modlisation de la tarification EDF

Nota : le lecteur pourra se reporter la rubrique Rseaux lectriques. Planification et
exploitation du trait Gnie lectrique.
Comme lindique la figure 11, une version tarifaire (ici 5 postes)

est modlise par :
des souscriptions incrmentales de puissance X1 , X 2 , ..., X 5
factures respectivement aux cots unitaires dcroissants P1 , P2 ,
..., P5 , correspondant une prime fixe ;
des utilisations V i , sur la priode , de la puissance souscrite

sur le i me poste tarifaire ; elles sont factures aux prix Pei de
lnergie sur le i me poste tarifaire ; lappel au rseau T 4 sur la
priode est alors la somme des utilisations autorises V i des diffrentes puissances souscrites.
Figure 9 Exemple dun systme nergtique modlis
B 1 282 11
La puissance X1 souscrite en priode 1 sera utilise en priode 2

puisque P1 est dj pay. Tout supplment de puissance de la
priode 2 sur la priode 1 sera souscrit en X2 au cot unitaire P2 et
non en X1 supplmentaire car P2 est moins cher que P1 (le programme minimise les cots).
Les contraintes de prise en compte de la tarification lectrique
sont alors :
1
T 4
2
T 4
3
T 4
4
T 4
5
T 4
= V1
2
= V1+V2
3
V1
4
V1
5
V1
=
=
=
+
+
+
3
V2+
4
V2+
5
V2+
V3
4
V3+V4
V 1 X1
2
V 1 X1
3
V 1 X 1
4
V 1 X 1
5
V 1 X 1
V 3 X 3
4
V 3 X3
5
V 3 X3
V 2 X 2
3
V 2 X2
4
V 2 X2
5
V 2 X2
V3+V4+V5
Figure 10 Courbe caractristique dun groupe turboalternateur
V 4 X4
5
V 4 X4
V 5 X5
Les significations des diffrentes variables ( V i , X i ) sont explicites sur lexemple de la figure 12.
2.2.3 Liste des contraintes

Les diffrentes contraintes sont listes par type, avec entre parenthses gauche leur nombre. Les contraintes dexploitation (autoproduction minimale de 1 MW, fonctionnement de la chaudire) sont
donnes titre dexemple ; dautres types de phnomnes peuvent
tre pris en compte.
Les structures tarifaires du gaz et de llectricit retenues pour cet
expos ne prsentent quune valeur dexemple pour lcriture des
contraintes quelles engendrent. Il est vident que ces structures tarifaires peuvent voluer au cours du temps et que la prsente
modlisation sapplique toute forme dnergie.
Figure 11 Modlisation de la tarification EDF : version 5 postes
Bilans emplois-ressources des produits (pour = 1, 2, 3, 4, 5, 6)
1
T 2
b2
) 1,183
T3
avec
(6) T 1 + ( b 1 t 1 T 2 + b 2 t 2 T 2 ) + T 3 0
(6)
(b 1
t1
t2
2
T 2
2
Capacits des quipements (pour = 1, 2, 3, 4, 5, 6)
T5
Qv
(6) ( b 1 t 1 Z + b 2 t 2 Z ) T 4 + 0,652 T 5 Q e
o 1,183
0,652
correspond la quantit de vapeur basse pression

fournie par la dtente-dsurchauffe, pour 1 tonne de
vapeur haute pression,
correspond la consommation spcifique de la
chaudire lectrique (en MWh/t de vapeur).
Contraintes de stockage - dstockage

(6) E + 1 = E + (S + 1 R + 1) D + 1
(4) E C s ( = 2 5 )
B 1 282 12
D (heures) dures respectives des priodes .
et E 1 = E 6 = 0
(6) T 1 100
(6)
(6)
(6)
(6)
T2
T3
T4
T5
75
100
20
X6 0
Ces bornes suprieures correspondent aux limites de fonctionnement pleine charge des quipements.
Contraintes spcifiques dexploitation
Autoproduction minimale :
(6) ( b 1 t 1 Z
+ b 2 t 2 Z ) 1
Figure 12 Signification des variables V i et Xi dans la modlisation de la tarification EDF
Structure tarifaire du gaz
Chaudire charge non modulable sur la journe :

(4)
1
T1
=T
2
1
1
T1
= T
3
1
4
T1
= T
5
1
4
1
= T
6
1
1
1
2
2
3
3
24
( 1 ) 0,75 [ T 1 D + T 1 D + T 1 D ] ------------------------------------- X7 0
D1 + D2 + D3
Contraintes sur les variables entires (modlisation de la turbine)
4
4
5
5
6
6
24
( 1 ) 0,75 [ T 1 D + T 1 D + T 1 D ] ------------------------------------- X7 0
D4 + D5 + D6
(6) b 1 t 1
(6) t 2 b 2
(6)
(6)
(6)
t11
t21
b2b1
0
B 1 282 13
Structure tarifaire de llectricit

1
1
(1) T 4 V 1 = 0
2
2
2
(1) T 4 V 1 V 2 = 0
3
3
3
3
(1) T 4 V 1 V 2 V 3
=0
4
4
4
4
4
(1) T 4 V 1 V 2 V 3 V 4 = 0
5
5
5
5
5
(1) T 4 V 1 V 2 V 3 V 4 = 0
6
6
6
6
6
6
(1) T 4 V 1 V 2 V 3 V 4 V 5
1
V 1 X1 0
X1 0
5
V 1 X1 0
6
V 1 X1 0
( 5 ) V 2 X2 0
V 2 X2 0
3
V 2 X2 0
4
V 2 X2 0
V 2 X2 0
X1 0
5
6
Nota : on rappelle que 1 th = 4,185 5 10 6 J, la thermie ntant plus une unit lgale.
( 4 ) V 3 X3 0
V 3 X3 0
6
X3 0
V 3 X3 0
( 3 ) V 4 X4 0
V 4 X4 0
5
V4
Lexemple prcdent a t rsolu avec les donnes suivantes :

enthalpie du rseau de vapeur haute pression : 770 th/t ;
enthalpie du rseau de vapeur basse pression : 650 th/t ;
rendement de la chaudire existante : 85 % ;
besoins du process en vapeur basse pression et en lectricit
(en trait fort, figure 13) ;
tarif lectrique retenu correspondant au barme du tarif EDF
5 postes de taille A (article Tarifs de llectricit en France : barme
des prix [D 4 935], dans le trait Gnie lectrique) en vigueur au
01/01/1987 ;
charges annuelles dinvestissement admises pour la chaudire
lectrique : 100 kF/ MW, ce qui correspond un temps de retour
denviron 3 ans.
=0
2
V 1
3
V 1
( 6 ) V 1 X1 0
4
V3
immdiatement que le traitement dun cas rel (plusieurs dizaines

dquipements et de priodes) conduit des problmes de grandes
tailles ncessitant :
le recours oblig aux programmations linaire (algorithme du
simplex) ou mixte ;
un programme gnrateur des contraintes et de la fonction
conomique minimiser partir des informations caractristiques
du cas tudier.
X4 0
( 1 ) V 5 X5 0
2.2.4 Fonction conomique minimiser

Cette fonction comporte les termes suivants :
cot de lnergie lectrique ;
prime fixe de souscription de llectricit ;
cot de lnergie gaz : deux tranches de tarification de cots
Pc 1 , Pc 2 , la premire tranche tant relative une nergie donne A
(valeur contractuelle) ;
prime fixe de souscription du gaz FG et charge annuelle F lie
lventuel investissement de la chaudire lectrique.
Cot total minimiser :
1
Min { [ Pe1 T 4 D + Pe 2 T 4 D + Pe3 T 4 D

4
+ P e 4( T 4 D + T 4 D ) + P e5 T 4 D ]
2.3 Conclusion
+ ( P1 X1 + P2 X2 + P3 X3 + P4 X4 + P5 X5 )
1
+ [ P c 2 0,75 ( T 1 D + T 1 D 2 + T 1 D
4
+ T 1 D 4 + T 1 D + T 1 D A ) + A 0,75 P c1 ]
+ ( FGX 7 + FX 6 ) }
2.2.5 Rsolution
Cet exemple se prsente donc sous la forme dun problme
doptimisation en variables mixtes comportant 125 contraintes,
72 variables continues et 12 variables boolennes pour seulement
4 quipements et une anne compose de 6 priodes. On conoit
B 1 282 14
Pour un prix du combustible fossile sous chaudire gal

0,10 F/ th PCI, les rsultats indiquent notamment que (figure 13) :
en heures creuses dt, lexploitant a conomiquement intrt
dtendre statiquement la vapeur et acheter la totalit de lnergie
lectrique sur le rseau externe EDF ;
sur les autres priodes tarifaires, la stratgie optimale consiste
turbiner les besoins de vapeur du procd et acheter sur le rseau
lectrique externe le complment des besoins lectriques de ce
mme procd ; il convient de remarquer que larrt-dmarrage du
groupe turbo-alternateur entre heures pleines et creuses dt donc
au cours de la mme journe est difficile raliser en pratique,
mme si cela correspond un optimum conomique ; pour une
application relle, il conviendrait donc dadjoindre au problme une
contrainte supplmentaire traduisant le minimum technique du
groupe turboalternateur en heures creuses dt ;
le prix du combustible fossile est trop bas pour justifier linvestissement dune chaudire lectrique.
Lorsque le prix du combustible fossile sous chaudire est
de 0,20 F/ th PCI, on constate (figure 14) que :
linvestissement dune chaudire lectrique de 20 t/h est
conomiquement justifi ; elle fournit alors la totalit de la vapeur
basse pression en heures creuses dt et une partie seulement en
heures pleines, le complment tant assur par la dtente statique ;
lexploitant a conomiquement intrt stopper lautoproduction lectrique dorigine thermique pendant la totalit de lt
tarifaire.
Parmi les rsultats correspondant lexemple donn, on trouve

notamment :
la capacit optimale X6 de lquipement dont on envisage
linvestissement ;
le niveau des souscriptions de puissance X1 , ..., X 5 de lnergie
lectrique ;
les niveaux dactivits T 2 et T 3 du groupe turboalternateur

et de la dtente statique ;
les achats des diffrentes nergies au cours du temps en fonction de leurs prix et des besoins du procd.
La modlisation retenue permet donc de dterminer le choix et
le dimensionnement optimal dun investissement nergtique, le
choix optimal des souscriptions des diffrentes nergies, la gestion
optimale des appareils producteurs ou transformateurs dnergie,
ainsi que les achats dnergie lextrieur du site tudi. Elle rpond
donc bien aux objectifs annoncs dans lintroduction.
Figure 13 Besoins du process en vapeur basse pression

et en lectricit
Figure 14 Besoins du process en vapeur basse pression

et en lectricit pour un prix du combustible fossile du double
par rapport la figure 13
B 1 282 15
3. Optimisation des schmas

de transferts dnergie
dans les procds industriels
3.1 Position du problme
Dans un procd industriel, le schma des transferts dnergie est
conu en essayant de maximiser le rendement nergtique global,
cest--dire en valorisant au mieux lnergie utilise dans le procd.
Cette optimisation nergtique saccompagne dune optimisation
conomique tenant compte des investissements et des cots de
fonctionnement.
Une connaissance prcise du procd et de ses contraintes opratoires est ncessaire : en chimie par exemple, les tempratures et
les pressions de raction devront tre respectes dans toute
proposition de nouveau systme de transferts dnergie.
Les procds industriels comportent frquemment de nombreux
besoins nergtiques. Lnergie peut tre utilise sous forme de force
motrice (pompes, compresseurs), transforme en nergie chimique
(raction endothermique) ou en nergie thermique (fours,
chaudires).
Paralllement ces besoins, les effluents thermiques internes au
procd (condensation, refroidissement de produits, etc.)
constituent des sources dnergie. Cela sexplique par le fait que les
quantits dnergie mises en jeu sont souvent trs suprieures aux
besoins thoriques : dans une colonne de distillation, par exemple,
la sparation de deux corps ncessite un appoint dnergie important
au bouilleur, alors que le travail minimal de sparation est le plus
souvent faible. La diffrence entre lnergie pratiquement ncessaire
et celle thoriquement ncessaire se retrouve au niveau du
condenseur, une temprature plus faible que celle laquelle elle
a t fournie au bouilleur.
Afin damliorer le rendement nergtique du procd, on cherche
rutiliser cette nergie disponible, cest--dire laffecter un
besoin interne au procd ou appartenant un autre procd de la
mme plate-forme industrielle. On conoit ainsi de nouveaux
schmas nergtiques ; ceux-ci peuvent tre nombreux, varis et
complexes car la diversit des sources et des besoins dnergie
sajoute celle des techniques de transferts dnergie (simple changeur si les tempratures des sources et besoins sont compatibles,
pompe chaleur si le niveau de temprature doit tre mont, etc.)
do lintrt de la modlisation de choix dquipements nergtiques (modles daide la dcision).
Un schma nergtique pourra tre retenu si lon dmontre, dans
un premier temps, sa faisabilit technique et, dans un second temps,
son intrt conomique.
En ce qui concerne la faisabilit technique, de nombreux critres
doivent entrer en ligne de compte avant toute proposition de
nouveau schma :
la technologie pour le matriel propos est-elle bien connue ? ;
cette technologie peut-elle adapte tre au process (matriel
antidflagrant, par exemple) ? Est-elle compatible avec les fluides
mis en jeu dans le procd ? ;
les performances souhaites sont-elles ralisables (taux de
compression, densit de transfert de chaleur, tenue des
matriaux) ? ;
lencombrement du matriel envisag nest-il pas rdhibitoire ?
Pour sassurer de lintrt conomique de ces rcuprations
nergtiques, plusieurs points devront galement tre rsolus :
lexploitation dune source dnergie de faible puissance ou
trop faible temprature nest pas forcment conomiquement
justifie ;
la quantit dnergie rcupre augmente si lon diminue le
pincement des changeurs (carts de temprature entre fluide chaud
et fluide froid) ; cela implique une augmentation de la surface
dchange, et donc aussi de linvestissement : la rcupration devra
tre pousse jusqu un optimum quil conviendra de trouver.
B 1 282 16
Chaque schma nergtique devra tre tudi dans son

ensemble : en effet, la valorisation dune source dnergie entrane
une modification du bilan thermique du procd et ne peut donc
tre tudie isolment.
Ces schmas nergtiques devront tre complts, le plus
souvent, par un apport dnergie externe au procd (production de
vapeur en chaudire, fours, nergie mcanique de compression,
etc.).
Lorsque les schmas de transfert dnergie paraissant les plus intressants ont t slectionns, puis modliss, il convient de dterminer le systme daffectation optimal des ressources aux besoins
dnergie. Cette optimisation devra tre ralise sous certaines
contraintes, garantes des faisabilits technique et thermodynamique, et des possibilits dadaptation des techniques
proposes au procd :
contraintes de procd (respect des tempratures et des
pressions de raction, quilibre du bilan nergtique, etc.) ;
contraintes techniques (limitation du taux de compression,
etc.) ;
contraintes thermodynamiques (compatibilit des tempratures et des pressions du fluide dcrivant un cycle dans une pompe
chaleur ou une compression mcanique de vapeur, etc.) ;
contraintes conomiques (limitation de la surface dchange,
etc.).
La fonction objectif (fonction optimiser) pourra tre une fonction
conomique F tenant compte de linvestissement et du cot
dexploitation. On cherchera, par exemple, maximiser la fonction
suivante :
F = E I
avec
reprsentant la diminution du cot dexploitation par

rapport la situation initiale (en francs/an) ; E dpend,
entre autres, du cot des diffrentes formes dnergie en
concurrence (gaz, lectricit, fuel, charbon),
temps de retour marginal admissible, fix au dpart (en
annes) ; en effet loptimum : dF = 0 donc dE dI = 0,
soit = dI/dE, ce qui correspond bien un temps de
retour de linvestissement marginal,
investissement d aux modifications du procd
impliques par le nouveau schma nergtique (en
francs).
Cette fonction conomique, tout comme les contraintes ou les

quations de modlisation, peuvent prsenter un caractre fortement non linaire (une surface dchange, par exemple, fait intervenir une fonction logarithme). La mthode adopte devra donc offrir
la possibilit de travailler en non linaire.
La mthode du gradient projet, dont un algorithme a t dvelopp par EDF (Direction des tudes et Recherches), est prsente
sommairement ( 3.2.2).
Nota : il est possible doptimiser selon un autre critre : minimisation de la
consommation dnergie primaire, par exemple. Il suffit pour cela que la fonction objectif
puisse sexprimer en fonction des variables retenues dans le systme optimiser.
3.2 Une mthode doptimisation

3.2.1 Modlisation - simulation
Partant dun procd existant, il est ncessaire, dans un premier
temps, de recenser et de dcrire ses diffrents flux matire et nergie
laide de donnes intensives (temprature, pression, fraction
vapeur, composition molaire, etc.) et de donnes extensives
(dbit-masse, par exemple). On en dduit lensemble des sources
de chaleur disponibles, leur niveau de temprature, ainsi que les
besoins nergtiques du procd.
Souvent, dans ce schma de flux initial, des rcuprations nergtiques existent dj (rchauffage de ractifs par des produits, production de vapeur dans une chaudire de rcupration, etc.). La mthode
doit permettre de vrifier si le schma actuel est le plus performant,
et, dans le cas contraire, il convient de dterminer ce qui est

remettre en cause et par quoi le schma nergtique peut tre
complt.
La faisabilit des transferts nergtiques tant lie principalement
aux quantits dnergie disponibles et requises et leurs niveaux
de temprature respectifs, une analyse thermodynamique peut
permettre de suggrer de nouveaux systmes de transfert dnergie.
Lemploi de diverses techniques peut tre envisag. Leurs performances doivent guider la conception des nouveaux schmas.
ce niveau, une bonne connaissance du procd et une matrise
des contraintes inhrentes ce procd sont indispensables.
Loprateur, qui pourrait terme tre remplac par un systme
expert, doit choisir les grandes options de transferts dnergie qui
caractriseront les schmas jugs plausibles. Tous les cas daffectation des sources aux besoins avec les techniques existantes
pourraient tre envisags, mais le temps de calcul serait norme ;
de plus, la majorit des cas simuls seraient irralistes.
Les diffrentes techniques intervenant dans le schma de transferts dnergie sont dcrites par un jeu de paramtres connus au
dpart et des variables que la simulation dterminera.
Prenons lexemple dune compression mcanique de vapeur deau
aprs dtente dans un ballon (figure 15).
Neuf quations rgissent ce systme :
(0)
changeur
calcul du bilan thermique
calcul de la surface dchange
ballon
calcul du bilan matires
compresseur
calcul de la temprature de sortie
calcul de la puissance
dsurchauffe
calcul du bilan matires
pompe
calcul de la puissance
Le jeu de variables peut tre le suivant :

S
surface dchange ;
T1 temprature la sortie de lchangeur ;
T3 temprature la sortie du compresseur ;
Wc puissance du compresseur ;
Wp puissance de la pompe ;
q2 dbit de vapeur en sortie du ballon ;
q3 dbit de liquide en sortie du ballon ;
q4 dbit de dsurchauffe ;
qv dbit de vapeur.
Les autres grandeurs caractrisant le systme sont disponibles
dans la description initiale des flux du procd ou sont fixes
a priori. La simulation de lensemble des techniques composant le
schma nergtique fournit ainsi les grandeurs caractristiques du
nouveau systme de flux du procd (tempratures, dbits) et les
grandeurs caractristiques des techniques choisies (surface
dchange, puissance des pompes et compresseurs, etc.).
3.2.2 Optimisation technico-conomique

Les grandeurs calcules lors de la simulation permettent a priori
de dterminer approximativement les investissements et les cots
nergtiques (si lon a prcis, au pralable, un cot des diffrentes
formes dnergie) ; il est donc possible de calculer la fonction objectif
prcdemment dfinie : F = E I.
Figure 15 Schma dune compression mcanique de vapeur deau
Reprenons lexemple du systme prcdent : certains paramtres,

fixs a priori, peuvent tre considrs comme variables ; cest le cas,
par exemple, du dbit q1 et de la pression du ballon p 2 . Si lon ajoute
ces deux variables (dites variables de contrle) aux neuf variables
initiales, le systme dquations comporte deux degrs de libert.
On entre ainsi dans la phase doptimisation, o lon cherche maximiser un critre (ici conomique) en respectant certaines contraintes.
Dans le cas prcdent, le jeu de contraintes pourrait tre le
suivant :
: pression du ballon infrieure la pression
p2 p1
lentre de lchangeur ;
: dbit de dsurchauffe infrieur au dbit de
q4 q3
liquide sortant du ballon ;
q v q v max
: dbit de vapeur infrieur ou gal au besoin
total ;
: limitation de la temprature la sortie du
T 3 T 3 max
compresseur ;
T 1 T fe
compatibilit des tempratures changeurT1 T2

ballon ;
S S max
: limitation de la taille de lchangeur.

Dune manire plus gnrale, si un systme est modlis avec
n e quations et comporte n v variables (n v > n e), il existe n v n e
degrs de libert.
Plusieurs mthodes doptimisation non linaire existent, mais
aucune dentre elles nest universelle : une adaptation de la mthode
au problme trait est souvent ncessaire.
Les principaux types de mthodes existants (article Optimisation
des processus nergtiques : concepts [B 1 280] dans ce trait)
sont les suivants :
mthodes de pnalisation
relativement anciennes ;
mthodes de gradient
mthodes du lagrangien augment (1975) ;
mthodes quadratiques squentielles (1978).
B 1 282 17
Les deux derniers types de mthodes, encore trs jeunes, sont

promis vraisemblablement un dveloppement important, compte
tenu de la facilit de programmation pour le lagrangien augment
et des performances des mthodes quadratiques squentielles.
Des amliorations sont toutefois attendre.
Les mthodes de pnalisation ne sont pas trs efficaces, et lon
prfre gnralement une mthode de gradient qui prsente lavantage dtre particulirement robuste.
Les mthodes de gradient ont souvent une convergence lente mais
certaines dentre elles conviennent bien aux problmes que lon se
pose ici, comportant souvent de nombreuses quations (20 50) et
peu de degrs de libert (2 5).
Prsentons sommairement lune de ces mthodes : le gradient
projet.
Partant dun point (x1 , x2, ..., xn ) [n tant gal au nombre de
degrs de libert (n v ne )] reprsentant les diffrentes variables de
contrle du systme, on calcule le gradient de la fonction objectif ;
celui-ci fournit la direction dans laquelle il est souhaitable de faire
varier le point de dpart.
Un nouveau point est calcul, o le gradient est projet sur un
hyperplan tangent la contrainte limitante ; cela permet de ne
conserver que la composante du gradient qui reprsente la direction
dans laquelle la progression est possible tout en restant dans le
domaine ralisable. Le processus doptimisation sarrte lorsque les
variations du critre conomique sont infrieures une prcision
dfinie au pralable.
Il convient de sassurer dans cette mthode, comme dans bon
nombre des autres, que le rsultat nest pas un optimum local. Il
est conseill de raliser plusieurs optimisations avec des points de
dpart diffrents afin de confirmer loptimum.
3.3 Exemple dun procd de la chimie

Les procds de la chimie sont souvent trs complexes ; les
transferts dnergie y sont nombreux :
condensation en tte de colonne distiller ;
vaporation en fond de colonne distilller ;
refroidissement de produits dune raction exothermique ;
etc.
Lensemble de ces sources et besoins nergtiques peut tre reprsent laide dun diagramme temprature-puissance (figure 16) ;
il est ainsi possible dapprcier la faisabilit thermodynamique des
diverses rcuprations de chaleur envisageables et deffectuer une
premire slection des techniques de transfert.
tant donn la diversit des produits que lon peut rencontrer en
chimie, il est ncessaire, pour la simulation des transferts dnergie,
de disposer dune base de donnes thermodynamiques fournissant
les principales caractristiques des produits intervenant dans le
procd. De plus, le calcul des enthalpies pour lcriture des bilans
thermiques ncessite lemploi de modles thermodynamiques ;
ceux-ci sont plus ou moins prcis, plus ou moins compliqus mettre
en uvre et doivent souvent tre adapts aux produits.
Les schmas de transferts dnergie, une fois simuls puis optimiss comme nous lavons dcrit prcdemment, peuvent tre
compars dun point de vue nergtique et conomique. Divers
corps dhypothses dvolution des prix des diffrentes formes
dnergie peuvent tre envisags : on dtermine ainsi les cots
seuils de ces nergies assurant lintrt conomique des diffrents
schmas.
3.4 Conclusion
La mthodologie dcrite est adapte aux systmes de transferts
dnergie troitement lis aux procds. Elle permet de simuler un
grand nombre de ces systmes et de slectionner ceux qui peuvent
B 1 282 18
Figure 16 Analyse thermodynamique : diagramme temprature puissance
prsenter un intrt. Loptimisation nest possible que lorsque lon

connat lexpression des investissements en fonction des variables
de contrle des techniques simules. Elle est intressante pour la
comparaison de diffrents projets entre eux, dun point de vue
conomique.
Cette dmarche doit tre considre comme un filtre : elle limine
rapidement les schmas nergtiques sans intrt et conserve ceux
qui mritent une attention particulire et une tude dingnierie plus
fine.
4. Gestion de lnergie
en temps rel
Le critre ne jouant pas le rle de critre objectif, production autonome dans le premier cas cit ci-dessus, cot dexploitation dans
le second, peut tre ventuellement born et jouer alors un rle de
contrainte dans le calcul doptimisation.
Depuis le milieu des annes soixante dix, deux dveloppements

ont donn une importance croissante au problme de la gestion de
lnergie. Dune part, les industriels ont investi dans la transformation de leurs installations pour diminuer leurs consommations spcifiques (rcupration dnergie) et pour les rendre aptes utiliser
plusieurs sources dnergie primaire. Dautre part, la tarification des
nergies, dans un march devenu trs concurrentiel, en particulier
depuis 1985, sest considrablement affine.
Dans ce contexte, la tche de lexploitant est devenue trs
complexe. Ses liberts de choix, mais aussi ses contraintes
techniques, se sont accrues en mme temps que les conditions
conomiques devenaient plus variables, les prix relatifs des nergies
changeant toutes les semaines, voire tous les jours (lectricit bien
plaire prix spot des quantits limites de fuel). De stables long
terme et affichs, les prix des nergies sont devenus ngocis et
phmres.
La mise disposition de lexploitant doutils informatiques puissants et performants, laidant effectuer en permanence les
meilleurs choix possibles, est aujourdhui une ncessit dans la
mesure o ils permettent de tirer parti et rentabiliser au mieux les
investissements passs, en adaptant, en temps rel, les conditions
techniques dexploitation des installations un environnement
conomique constamment changeant.
4.1.1.2 Contraintes
De natures varies, on peut distinguer :
les contraintes purement techniques lies la conception et
au dimensionnement des quipements, telles que la capacit maximale dune chaudire, la pression maximale admissible dun rservoir ou le minimum technique dun brleur ;
les contraintes dexploitation imposes par le systme de rgulation, telles que le dmarrage en cascade de chaudires, le rapport
dbit dair/dbit de combustible ou la loi de variation du dbit et
de la temprature de dpart dun rseau de chaleur en fonction de
la temprature extrieure ;
les contraintes de quantits dnergie livres aux consommateurs dont la satisfaction des besoins est en gnral impose ;
les contraintes dapprovisionnement du site en gaz, lectricit,
hydrocarbures, eau, etc., exprimes au travers des contrats conclus
entre les fournisseurs et le site ;
les contraintes de qualit de la fourniture, telles que la siccit
de la vapeur (exemple dune opration de vaporisage en teinturerie),
le niveau de temprature du fluide caloporteur (exemple dune
pasteurisation du lait).
Comme indiqu prcdemment, la contrainte peut tre de nature
conomique lorsque le cot ne constitue pas le critre objectif.
4.1.1.3 Liberts de choix
4.1 Dfinitions
Avant toute chose, il apparat ncessaire de prciser ce que
recouvrent certains termes couramment utiliss dans le domaine de
la gestion optimale en temps rel de systmes nergtiques.
4.1.1 Optimisation
Optimiser un processus de production nergtique ou autre na
de signification que si lon prcise :
le critre dvaluation retenu ;
le cadre des contraintes techniques, dapprovisionnement de
matires de base (combustibles, eau, lectricit, etc.), de la qualit
de la fourniture, etc., dans lequel loptimum est recherch ;
les liberts de choix offertes lexploitant du systme considr
pour atteindre cet optimum.
4.1.1.1 Critre dvaluation
Le critre dvaluation ou encore critre objectif est ncessairement unique. Il peut tre de nature technique ou conomique selon
les buts poursuivis. On pourra chercher, par exemple :
minimiser un cot dexploitation ;
maximiser une production autonome dnergie lectrique ;
ou encore :
minimiser un risque de dfaillance dun matriel ;
minimiser la pollution acide mise par un site.
Chacun de ces objectifs correspond respectivement au souci
dassurer une production moins chre, plus forte, de meilleure
qualit, moins polluante.
Ces critres ne peuvent tre en aucun cas optimiss simultanment : ou bien lon recherche une production autonome dnergie
lectrique maximale au dtriment du cot dexploitation, par
exemple en hiver lors dune dfaillance du rseau lectrique
alimentant lusine, ou bien lon recherche un cot dexploitation
minimal et la production autonome dlectricit est alors un rsultat
fatal du calcul doptimisation.
Par sa nature, un systme ne peut tre optimis que sil existe

un choix alternatif la solution propose initialement. Dans
lexploitation dun systme nergtique global, lallocation de charge
entre plusieurs chaudires, le niveau de temprature du dgazeur,
la pression dun collecteur de vapeur, lentranement dune pompe
ou dun ventilateur soit par moteur lectrique, soit par turbine, le
dmarrage ou larrt de certains quipements, constituent autant de
variables sur lesquelles on peut agir pour optimiser le critre
dvaluation retenu.
Au niveau dun quipement comme une chaudire classique
quipe de son brleur, loptimisation du cot dexploitation passera
essentiellement par le rglage adquat dune seule variable, lexcs
dair.
4.1.1.4 Optimisation locale et optimisation densemble
Dans un systme nergtique donn, on a coutume de distinguer
deux niveaux doptimisation :
le premier dit local concerne les diffrents quipements,
composant le systme, envisags sparment ;
le second dit densemble considre le systme nergtique
comme un tout dont les lments interagissent et ventuellement
sopposent ; ce second niveau se justifie dans la mesure o la somme
doptimums partiels ne conduit pas, en gnral, un optimum
global.
Optimisation locale : le cas typique dune optimisation locale est
constitu par une chaudire. Son cot dexploitation sera dautant
plus faible que la combustion sera mieux rgle, comme lindiquent
les courbes classiques de pertes en fonction de lexcs dair
(figure 17).
Il sagit bien dun processus doptimisation, car il existe :
un critre objectif, qui est le rendement de la chaudire ; ce
critre apparemment technique est en fait conomique puisquil est
une traduction directe du cot dexploitation ;
au moins une contrainte, qui est la quantit de vapeur dlivre
par la chaudire une pression et une temprature donnes ;
une libert de choix, qui est lexcs dair de combustion.
B 1 282 19
La mise en uvre de cette optimisation sera plus ou moins bien

ralise selon que lon dispose dun simple rgulateur analogique
ou dun contrle numrique disposant dalgorithmes de calculs
sophistiqus.
Loptimisation locale est de nature dynamique (article Optimisation des processus nergtiques : concepts [B 1 280] dans ce trait).
Optimisation densemble : dans lexemple prcdent, le processus doptimisation est indpendant de lenvironnement de la
chaudire. Que celle-ci alimente un changeur ou une turbine
vapeur, lenvironnement reste identique.
Au contraire, si lon considre un ensemble chaudire-groupe
turboalternateur, la pression et la temprature de la vapeur dlivre
sont des paramtres importants du point de vue de loptimisation
du systme : plus la pression et la temprature sont leves, plus
la turbine produira dnergie mcanique dbit de vapeur donn,
mais plus le rendement de la chaudire diminuera. Selon le cot de
llectricit un instant donn, on aura intrt sacrifier le rendement de la production de vapeur au profit de la production dlectricit et, inversement, si llectricit est peu chre, on aura intrt
baisser la pression et la temprature de la vapeur jusqu des valeurs
minimales compatibles avec la conception de la chaudire et
arrter le groupe turboalternateur, la vapeur tant dtendue statiquement dans le by-pass.
Un problme doptimisation densemble suppose ncessairement
que lon cre un modle mathmatique de linstallation considre
et que lon dispose dalgorithmes de rsolution et doptimisation de
ce systme dquations reprsentant le fonctionnement de chaque
quipement et de leurs interactions. Loptimisation a alors pour objet
de rechercher les valeurs de consigne des rgulations, les conditions
de marche des quipements et lutilisation des combustibles
disponibles conduisant la meilleure valeur possible du critre
retenu.
4.1.2 Temps rel

La conduite dinstallations en temps rel a souvent t associe
la notion de rapidit de rponse des moyens de gestion mis en
uvre. En ralit, le fait quun systme de contrle-commande
rponde en temps rel une perturbation dpend essentiellement
de linertie de linstallation contrle.
Lorsquil sagit de grer des alarmes ou des scurits, il est en

gnral important de ragir dans des temps trs courts, allant de
quelques millisecondes quelques secondes.
Il en est de mme pour la plupart des fonctions de rgulation telles
que :
rgulation de vitesse dune turbine gaz (une centaine de
millisecondes) ;
rgulation de pression dune chaudire (quelques secondes
une minute).
Dans tous ces cas, la rapidit de lorgane de gestion est importante.
Lorsquil sagit doptimiser le fonctionnement dun quipement,
la rapidit de rponse est moins fondamentale et se situe des temps
de lordre de la minute (exemple dun rglage de lexcs dair dun
brleur).
Une optimisation densemble ne peut tre, quant elle, que
relativement lente. Le temps de rponse ncessaire est alors de
quelques minutes une demi-heure, selon ltendue du systme
nergtique considr.
4.2 Fonctions gnrales de gestion

et leur hirarchie
4.2.1 Fonctions
La gestion dun systme nergtique met en uvre diffrentes
fonctions qui sont classes par niveaux de complexit croissante
(figures 18 et 23).
Niveau 0 - Fonction dacquisition de donnes, de rgulation
simple (boucles PID), dactions rflexes tout ou rien : ce niveau,
les matriels utiliss sont essentiellement des capteurs-transmetteurs- convertisseurs, des rgulateurs et des automates simples ; ils
sont implants localement. Les temps de rponse demands sont
de lordre de 0,3 s.
Nota : boucle PID (proportionnelle, intgtrale et drive), article Principes des chanes
de rgulation [R 7 090] dans le trait Informatique industrielle.
Niveau 1 - Fonction de contrle avanc

Celui-ci porte sur des variables labores (exemple, la puissance)
et comportent des fonctions de calcul plus ou moins complexes telles
que : algorithme de dcouplage, calculs prdictifs, compensation
dynamique (avance, retard, temps mort, etc.).
Figure 18 Hirarchie des fonctions de contrle-commande
Figure 17 Optimisation de lexcs dair de combustion
B 1 282 20
Le contrle avanc repose sur des automates programmables

disposant de quelques kilo-octets ou quelques dizaines de
kilo-octets de mmoire vive implants localement ou, ventuellement, en salle de contrle pour de petites installations et selon la
technologie utilise.
Les temps de rponse demands sont de lordre de quelques
secondes environ 1 min.
Niveau 2 - Fonction doptimisation locale
Celle-ci porte sur des quipements individualiss ou ventuellement un ensemble homogne (exemple dune optimisation de
lallocation de charge entre plusieurs chaudires). Elle met ncessairement en uvre un modle de fonctionnement de lquipement
ou des quipements considrs. Le calculateur supportant cette
fonction tait jusqu prsent implant au niveau central compte tenu
de la puissance ncessaire des appareils du point de vue capacit
de calcul et espace mmoire vive. Leur miniaturisation, la baisse des
cots et laugmentation de leurs performances tendent aujourdhui
dcentraliser la fonction doptimisation locale proximit des quipements concerns.
Les temps de rponse demands sont compris entre 1 min
et 20 min.
Niveau 3 - Fonction doptimisation densemble
Celle-ci est fonde sur un modle labor du systme nergtique ; elle fait appel des algorithmes de rsolution et doptimisation sophistiqus exigeant des ressources informatiques importantes se situant, dbut 1988, la frontire entre les mini- et les
micro-ordinateurs.
Des mmoires vives de plusieurs mgaoctets, des vitesses de
calcul suprieures ou gales au Mips (million dinstructions par
seconde) des microprocesseurs 32 bits sont alors ncessaires pour
mener bien des calculs dont la dure peut tre comprise entre 10
et 60 min, ce qui reste compatible avec linertie des systmes
physiques considrs. La fonction doptimisation densemble est par
nature centralise dans la mesure o elle suppose une concentration
de linformation.
Selon les dimensions et la complexit des installations nergtiques envisages et le degr dlaboration de leur systme de gestion,
on met en uvre un ou plusieurs des niveaux prcdents. Les
niveaux 0 et 1 suffisent pour grer correctement de petites installations nergtiques. Les grandes plates-formes industrielles, raffineries de ptrole, usines chimiques ou chauffages urbains importants
ont intrt enrichir leur outil de gestion jusquau niveau 3.
Les niveaux 4 et 5 font partie, fin 1988, du domaine du futur et
dpassent le simple domaine de lnergie.
Le niveau 4 correspond loptimisation gnrale des conditions
techniques dexploitation (niveaux 1, 2 et 3) en relation avec les
contraintes dapprovisionnement de lusine et surtout de livraison
des produits finis la clientle. ce niveau, il est rpondu la
question : quel est le rglage de mes units de production permettant
de satisfaire mon planning de livraison, un cot minimal
dexploitation ?
Le niveau 5 intgre, en plus, la dimension financire de lentreprise
et correspond donc une optimisation globale et complte de ses
diffrentes fonctions en tenant compte de toutes les interactions et
contraintes. Comme image on peut dire qu ce niveau, le rglage
du point de consigne de fonctionnement dun quipement est
analys en relation avec son impact sur le bilan de lentreprise.
certains consommateurs en fonction de seuils de puissance atteints

et/ou selon des tranches horaires prdfinies.
Ces appareils programmables, dans leurs versions les plus
simples, ne travaillent qu partir dinformations du type tout ou rien
et ne comportent aucune fonction de calcul. Leur action est galement du type tout ou rien.
Des systmes de ce genre permettent dapporter des gains de
lordre de 10 %.
4.2.2.2 Niveau 1 : dmarrage optimis dun four lectrique
La solution de rfrence I (figure 19) consomme 25 % dnergie
en plus de loptimum thermique solution II et cote deux fois plus
cher que loptimum conomique solution III. La solution IV peut
tre considre comme un optimum technico-conomique car elle
vite les pointes de temprature tout en restant conomiquement
acceptable.
Cet exemple illustre bien lavantage quapportent les microprocesseurs reprogrammables, car ils peuvent sadapter trs
rapidement de nouvelles conditions dexploitation, de tarification
ou de contraintes technologiques.
Les gains apports par de tels quipements sont de lordre de 10
20 %.
4.2.2.3 Niveau 1 : amlioration de la rgulation du systme
nergtique dune papeterie
Dans lindustrie de fabrication du papier, il est relativement
frquent quau cours de son schage, le papier casse brutalement,
ce qui conduit un arrt quasi instantan de la production et une
chute brutale des consommations de vapeur.
Il est, de plus, courant que ces industries disposent dun groupe
turboalternateur dtendant la vapeur produite par des chaudires
haute pression jusqu un niveau directement utilisable par les
machines papier.
Les systmes classiques de rgulation taient fonds sur le
principe suivant (figure 20) : peu de temps aprs la casse du papier,
la pression en aval de la turbine a tendance monter, ce qui a pour
effet de fermer sa vanne dadmission entranant son tour une
augmentation de pression au niveau de la chaudire do il rsultera
une baisse des feux. En gnral, la brutalit du phnomne conduit
invariablement faire fonctionner les soupapes de la chaudire.
Le systme de contrle avanc indiqu sur la figure 20 va
permettre danticiper directement la rduction de la charge de la
chaudire dans des conditions dpendant des niveaux de pression,
dbit, etc. existant au moment de la casse et dviter toute perte de
vapeur.
4.2.2 Exemples dapplication

4.2.2.1 Niveau 0 : dlestage/relestage de puissance lectrique
Dans le but de suivre au mieux le contrat dapprovisionnement
en lectricit souscrit auprs dEDF, il est maintenant frquent
dutiliser un petit automate permettant de dlester ou relester
Figure 19 Diagramme de puissance de chauffage

dun four lectrique en heures creuses
B 1 282 21
Les gains nergtiques dans un tel cas peuvent atteindre quelques

diximes de pour-cent de la consommation annuelle.
4.2.2.4 Niveau 2 : suivi optimis du contrat EDF
Limportance des primes de puissance souscrite a conduit certains
gros consommateurs dlectricit disposant de moyens leur
permettant de moduler leur appel de puissance sur le rseau (par
exemple, groupes turboalternateurs, process de fabrication adaptable) squiper doutils sophistiqus de suivi du contrat EDF.
EDF ne comptabilisant la puissance appele qu partir de lnergie
moyenne consomme sur une priode de 10 min (figure 21), il est
alors possible dtablir des programmes de calcul et des stratgies
daction qui, en fonction de lallure de la courbe de consommation
lintrieur de cette priode et de diffrents paramtres dexploitation (niveau de charge des groupes turboalternateurs et temps de
rponse des chelons de puissances par exemple), vont permettre
de moduler la consommation sur le rseau EDF pour tangenter le
contrat.
Le problme peut tre assez complexe rsoudre, puisquil doit

prendre en compte des stratgies daction qui sont en gnral diffrentes selon les tranches horo-saisonnires, le fait que certains matriels puissent tre larrt pour des raisons de maintenance, etc.
4.2.2.5 Niveau 2 : allocation de charge entre chaudires
Dans la mesure o une installation comporte plusieurs chaudires
de caractristiques diffrentes, il existe un optimum de rpartition
de la production entre ces diffrentes units pour satisfaire un besoin
donn de vapeur. Le rglage correspondant nest pas toujours intuitif, comme le montre le cas dtaill ci-aprs.
Considrons deux chaudires fonctionnant en parallle (figure 22)
utilisant le mme combustible dont les caractristiques sont :
chaudire A : dbit compris entre 40 et 100 t/h de vapeur avec
un rendement compris entre 90 et 92 % ;
chaudire B : dbit compris entre 20 et 50 t/h de vapeur avec
un rendement compris entre 88 et 90 %.
La quantit de vapeur produire est de 120 t/h.
Intuitivement, lexploitant charge la chaudire de plus faible rendement son minimum. Dans ces conditions, la chaudire A produit
100 t/h et la chaudire B 20 t/h.
En fait, ce rsultat nest pas optimal comme le montre un calcul
suffisamment fin tenant compte du fait que lvolution des rendements nest pas linaire en fonction de la charge.
En effet, le rsultat obtenu par une optimisation non linaire est
le suivant : la chaudire A doit tre charge 78 t / h et la
chaudire B 42 t/h.
Cette situation correspond au fait quau point considr les cots
marginaux de la vapeur pour les deux chaudires sont identiques
ce qui, mathmatiquement, sexprimera par lgalit des drives
partielles du cot de la vapeur par rapport aux dbits de chaque
chaudire, traduisant bien lexistence dun optimum au point
considr.
4.3 Structure des systmes

de contrle-commande
Les systmes modernes de contrle-commande sont caractriss
par une structure dcentralise et hirarchise (figure 23).
Aux niveaux infrieurs, les fonctions dacquisition et de rgulation, quil sagisse densembles simples (rgulateurs PID) ou relativement labors (rgulation prdictive), sont rparties localement
et ralises au moyen de microprocesseurs situs en gnral
proximit des units.
Figure 20 Contrle amlior du systme nergtique
dune papeterie
Figure 21 Allure des consommations lors dun suivi optimis

du contrat EDF
B 1 282 22
Figure 22 Allocation de charge entre chaudires
On peut cependant considrer que les possibilits des matriels

sont souvent trs fortement sous-employes en raison de labsence
de logiciels dapplication adapts aux cas concrets rencontrs.
4.4 Logiciels
Figure 23 Structure dcentralise et hirarchise

des grands systmes de contrle-commande
Les petites installations, en particulier celles que lon rencontre

dans les petites et moyennes entreprises, ne ncessitent pas en gnral de sophistiquer davantage le systme de contrle-commande.
Seule une centralisation des informations pour savrer utile pour
mieux suivre certains ratios de consommations et donc agir au mieux
moyen terme.
Au contraire, les plates-formes nergtiques complexes ne
peuvent plus tre exploites partir du simple bon sens comme la
montr lexemple du paragraphe 4.2.2.5. Des moyens puissants de
calcul, exigeant des quantits importantes dinformations,
deviennent ncessaires.
Le systme de contrle-commande doit alors sappuyer sur une
interconnexion gnrale des quipements au moyen dun cble
coaxial (bus) qui va permettre :
dune part, un change permanent et intense dinformations
(1 plusieurs Mbit / s) entre les diffrents microprocesseurs qui
auront alors la possibilit dintgrer dans leurs logiques de contrle
des variables externes aux sous-ensembles quils contrlent ;
dautre part, un renvoi dinformations vers un poste central qui,
aprs un traitement appropri de ces donnes, par exemple des
calculs doptimisation, renverra des consignes labores vers les
diffrents microprocesseurs installs sur le site.
Jusquau tout dbut des annes 80, pratiquement toutes les fonctions de calcul doptimisation de fonctionnement dquipements ou
de sous-ensembles dquipements taient effectivement ralises
au niveau central. La baisse du cot des mmoires vives, laccroissement des vitesses de calcul, le durcissement des matriels vis--vis
denvironnements industriels agressifs font que lon tend actuellement dporter ces fonctions un niveau local.
Le poste central tend en contrepartie senrichir de fonctions
doptimisation globales qui ncessitent lutilisation de logiciels
lourds et donc des matriels informatiques puissants.
Une dizaine de constructeurs (CGEE-Alsthom, Hartmann et Braun,
Foxboro, etc.) proposent aujourdhui des systmes numriques de
contrle-commandes capables de traiter toutes les fonctions de
contrle-commande.
Dans de telles installations, lcran vido couleur et des claviers
spcialiss ont entirement limin les traditionnels boutonspoussoirs, lampes tmoins et verrines dalarme quipant les
panneaux de contrle traditionnels. La totalit du dialogue
homme-machine passe par lintermdiaire de moyens informatiques
utilisant des codes alpha-numriques et des symboles graphiques.
Tous les constructeurs disposent de logiciels standards permettant

de configurer facilement les diffrentes boucles de rgulation, aussi
complexes soient-elles, ou de programmer leurs automates.
Les logiciels plus sophistiqus, permettant dtablir des bilans
masse et nergie dtaills et fiables, doptimiser des quipements,
voire doptimiser un systme nergtique complet, sont beaucoup
plus rares.
En pratique, sauf exceptions (contrle de lexcs dair de combustion, par exemple), chaque cas est trait sparment et individuellement : de nombreux lves-ingnieurs au cours de leurs stages de
fin dtudes auront ainsi crit des programmes dapplication spcifiques un site, en langage FORTRAN ou BASIC. Linconvnient de
cette mthode est que, en gnral, eux seuls connaissent les limites
du programme et que, les sites voluant, il devient rapidement
inapplicable.
De grands efforts sont faits actuellement pour crer des logiciels
ou plutt des progiciels facilement configurables en fonction de la
ralit dun site donn et donc permettant de sadapter aux volutions invitables de celui-ci.
La socit SERETE, en collaboration avec les socits Elf,
CGEE-Alsthom et lInstitut Franais du Ptrole, a fait des efforts
importants dans ce domaine pour laborer un logiciel configurable
doptimisation globale de systmes nergtiques complexes,
logiciel dnomm SECI-Manager. Bas sur des principes doptimisation non linaire avec prise en compte intgrale de variables
discrtes, reprsentatives en particulier des tats marche ou arrt
dquipements, celui-ci est oprationnel de manire industrielle
depuis 1986 aprs cinq annes de dveloppement.
Le choix de logiciels industriels, oprationnels et intgrables dans
des systmes de contrle-commande reste encore extrmement
restreint. Il en rsulte un sous-emploi des matriels disponibles dont
la puissance de traitement est souvent bien suprieure ce qui est
rellement utilis.
4.5 Rentabilit dune gestion optimale

en temps rel
La rentabilit que lon peut escompter dune gestion optimale en
temps rel dun systme nergtique a deux origines bien distinctes.
Un premier bnfice peut tre obtenu par la mise en uvre de
fonctions de niveau 0, 1 et 2 dcrites au paragraphe 4.2.1. Celles-ci
ont, en effet, pour consquence de rduire les variations des paramtres de conduite de linstallation de ce qui permet de les maintenir
proximit de valeurs spcifies ou de valeurs limites (exemples
cits au paragraphe 4.2.2.
Linstallation travaillant en permanence dans de bonnes conditions, il en rsulte une premire conomie dexploitation. Pour
quantifier cette conomie, il est bien entendu ncessaire de disposer
dun cas concret et dune situation de rfrence par rapport auxquels
il y a lieu de raisonner. Le temps de retour peut tre compris entre
quelques mois pour une installation neuve et trois ans lors du
remodelage dune installation existante.
Un second bnfice peut tre obtenu par la mise en uvre dune
fonction de niveau 3 ( 4.2.1) qui va permettre de dfinir des points
de consigne dexploitation optimums grce lutilisation dun
modle densemble reliant les performances conomiques de linstallation ses variables de contrle.
B 1 282 23
Dans ce cas, il est en gnral bien difficile de dterminer ou

mme destimer a priori les gains que lon peut attendre dune telle
optimisation, le rsultat tant rarement intuitif mme dans des cas
apparemment simples et parfaitement dlimits (exemple cit au
paragraphe 4.2.2.5).
Lexprience pratique montre cependant que les gains se situent
entre 0,5 et 5 % selon les cas et la complexit du problme pos.
4.6 volutions prvisibles

Dans le domaine de matriels, la gestion optimale en temps rel
des systmes nergtiques a pu prendre son essor au dbut des
annes 80 grce au dveloppement considrable de la microinformatique et surtout lamlioration constante des matriels en
matire de robustesse, fiabilit, prix, performances et rsistance
des environnements industriels.
Des systmes de contrle entirement centraliss du dbut des
annes 70, on est pass dans les annes 80 lutilisation de microprocesseurs locaux possdant des capacits mmoire et des puissances de calcul tendues.
Les annes 90 devraient voir une dissmination de lintelligence
jusquau niveau des capteurs, lensemble des quipements tant
relis entre eux par un bus unique.
Simultanment, le poste central dispose de microprocesseurs de

plus en plus puissants permettant une extension des fonctions de
calcul doptimisation et de traitement des donnes. On tend
connecter ces systmes avec des ordinateurs de gestion, les bases
de donnes devenant communes aux diffrentes fonctions techniques et financires de lentreprise.
Dans le domaine des logiciels, ceux facilitant le dialogue
homme-machine ne sont fortement dvelopps. Ils ont pour consquence de faire disparatre en partie le clavier qui se spcialise de
plus en plus, lcran tactile simposant progressivement avec une
imagerie et une dfinition samliorant constamment.
La configuration des systmes sappuie sur des logiciels plus
souples et plus faciles demploi.
Des logiciels de calcul doptimisation configurables sintressant
tout ou partie du systme nergtique devraient faire progressivement leur apparition.
Enfin, des logiciels plus globaux, nenvisageant le systme nergtique que comme un sous-ensemble du processus global de
fabrication de lentreprise, devraient voir le jour dici quelques
annes.
Il est probable quil sera alors ncessaire de disposer de matriels
informatiques nettement plus performants quactuellement si lon
veut que la gestion optimise que lon pourra raliser avec eux soit
effectivement en temps rel compte tenu de la taille des logiciels
quil y aura alors lieu de mettre en uvre.
Rfrences bibliographiques
[1]
[2]
[3]
BOITEUX (M.). Note sur le taux dactualisation. Revue dconomie politique, Paris (1969).
LEGOFF (P.). Optimisation nergtique des
procds industriels. Rev. gn. therm. (F),
nos 193 et 194 (1978).
CHAUVEL (A.), LEPRINCE (P.), BARTHEL (Y.),
RAIMBAULT (C.) et ARLIE (J.-P.). Manuel
dvaluation conomique des procds.
Technip (1976).
B 1 282 24
[4]
[5]
[6]
Applications de llectricit dans les procds

industriels. Ampre 82 (chapitre IP1) (F),
Librairie Lavoisier, Collection lectra (1986).
GJUMBIR (M.) et OLUJIC (Z.). Effective
ways to solve single nonlinear equations.
Chemical Engineering, juil. 1984.
CHINNECK (J.W.) et CHANDRASHEKAR (M.).
Models of large scale industrial energy
systems. Optimisation and synthesis. Energy,
vol. 9 (1984).
[7]
[8]
CHENG (R.C.H.) et JONES (G.). An automatic

method for modelling and analysing industrial
chemical plants. European journal of operational research, no 14 (1983).
GOURLIA (J.P.). Analyse exergtique de la
distillation. (Energtique industrielle 3.
Applications en gnie chimique, changeurs,
sparateurs, racteurs, de P. LE GOFF), Librairie Lavoisier, Collection Technique et
Documentation (1982).

b1280 Technique de L'ingénieur

Transféré par

Droits d'auteur :

Formats disponibles

b1280 Technique de L'ingénieur

Transféré par

Informations du document

Titre original

Copyright

Formats disponibles

Partager ce document

Partager ou intégrer le document

Options de partage

Avez-vous trouvé ce document utile ?

Ce contenu est-il inapproprié ?

Droits d'auteur :

Formats disponibles

b1280 Technique de L'ingénieur

Transféré par

Droits d'auteur :

Formats disponibles

Optimisation des processus

nergtiques : gestion optimale

Analyse des critres conomiques .....................................................

Constitution et gestion optimises des rseaux dnergie

Optimisation des schmas de transferts dnergie dans les

Gestion de lnergie en temps rel .....................................................

Rfrences bibliographiques .........................................................................

article Optimisation des processus nergtiques se compose de trois

OPTIMISATION DES PROCESSUS NERGTIQUES : GESTION OPTIMALE ___________________________________________________________________________

1. Analyse des critres

1.1 Critres de choix dun procd

au prix que le secteur productif est prt payer pour disposer de

1.1.2 Valeur du taux dactualisation choisir

1.1.1 Actualisation des cots futurs

__________________________________________________________________________ OPTIMISATION DES PROCESSUS NERGTIQUES : GESTION OPTIMALE

Entre deux projets dont lun procure un bnfice actualis certain

dure de vie T ......................(ans)

1.1.3 Critre conomique utiliser

Son usage revient privilgier les investissements, cest--dire

Les deuxime et troisime critres sont quivalents si les

Si les dures de vie des quipements sont diffrentes, il faut choisir

1.2 Optimisation dun processus

1.2.1 Multitude de variantes conomiques

OPTIMISATION DES PROCESSUS NERGTIQUES : GESTION OPTIMALE ___________________________________________________________________________

Lnergie ncessaire au processus nest gnralement pas de

F / U tant le courant consomm.

C total = I pro + I con + C ann

Prenons, comme exemple, un procd de concentration :

taux dactualisation annuel,

W (en kWh/ h) 0,04 Ft p

1.2.2 Principe doptimisation

__________________________________________________________________________ OPTIMISATION DES PROCESSUS NERGTIQUES : GESTION OPTIMALE

Le premier critre est de rendre minimal le cot total F ( E )

Tableau rcapitulatif des rsultats obtenus

Deux cas qui respectent ces quations peuvent se prsenter :

Nous ne dvelopperons pas ici le cas de plusieurs degrs de

Figure 1 Optimisation du cot dun processus

OPTIMISATION DES PROCESSUS NERGTIQUES : GESTION OPTIMALE ___________________________________________________________________________

Figure 2 Optimisation du temps de rcupration

Sur la relation (6), on constate que les paramtres fixant la

1.2.3 Influence du taux dactualisation

1.2.4 Influence de la dure de vie des quipements

Figure 4 Influence de la dure de vie de linstallation sur le cot total

__________________________________________________________________________ OPTIMISATION DES PROCESSUS NERGTIQUES : GESTION OPTIMALE

1.2.5 Influence de lenchrissement de lnergie

Ici, nous navons quun combustible : le fuel lourd. Dans un cas

1.2.6 Limites des paramtres

1.2.7 Discontinuits de linvestissement

Figure 5 Influence du cot des nergies

Si tous les appareils du processus sont de ce type, loptimisation

OPTIMISATION DES PROCESSUS NERGTIQUES : GESTION OPTIMALE ___________________________________________________________________________

1.3 Problme des donnes

1.3.2 Cots proportionnels

Figure 7 Effets de la discontinuit des investissements.

Le choix dun investissement nergtique et la gestion optimale

__________________________________________________________________________ OPTIMISATION DES PROCESSUS NERGTIQUES : GESTION OPTIMALE

pression conduit une modification de lallocation de charge entre

2.1 Principes de la modlisation

le fonctionnement du procd, gnralement fix par des