Brevet FR2773243A1 — publié le 19990702

DOMAINE TECHNIQUE La présente invention concerne un procédé de fabrication optimale d'un mélange ou d'une suite de mélanges de produits liquides, pâteux ou pulvérulents à partir d'une pluralité de produits appelés bases. Il trouve son application dans les industries chimiques, les unités de fabrication de lubrifiants, les industries agro-alimentaires, dans les raffineries de pétrole et les dépôts de produits pétroliers notamment pour la fabrication de carburants et de combustibles. ETAT DE LA TECHNIQUE ANTERIEURE La fabrication d'un produit conforme à des spécifications données, par mélange de produits appelés bases de caractéristiques physico-chimiques connues, comporte une étape dite de formulation qui à pour but de déterminer les proportions de chacune des bases. Cette fabrication doit en plus respecter des contraintes opérationnelles connues, telles que, dans le cas de produits liquides, les débits minima et maxima des pompes de transfert des bases, les quantités de bases disponibles, les volumes de creux dans les bacs de stockage du produit fabriqué. Pour un produit donné à fabriquer, on dispose en général d'un simulateur ou d'un modèle qui permet, à partir des caractéristiques connues par mesure des bases et de formulations empiriques, de vérifier que les caractéristiques du futur mélange seront acceptables. Dans le cas o l'on doit fabriquer plusieurs mélanges à la suite les uns des autres, il convient également de vérifier que les quantités des bases disponibles permettront d'effectuer la totalité des mélanges. Ces procédés de simulation et de modélisation sont longs à mettre en oeuvre et n'optimisent aucun critère. Un autre procédé connu utilisant les techniques de programmation linéaire permet d'optimiser le coût de fabrication d'un mélange quand les coûts des bases sont connus. Ce procédé est généralement utilisé pour établir des programmes mensuels de fabrication, il ne prend pas en compte les contraintes journalières d'exploitation des usines. Un autre procédé connu permet de maximiser la flexibilité, c'est à dire la possibilité de réaliser de futurs mélanges dont on ne peut pas prévoir les spécifications ou la demande. Ce procédé présente l'inconvénient de ne pas prendre en compte le coût de réalisation de chaque mélange, ni celui d'une suite de mélanges. En appliquant un autre procédé connu il est possible de prendre en compte plusieurs critères, par exemple un critère de coût Ci et un critère de flexibilité C2. Ce procédé consiste à: - affecter un poids à chaque critère soit a au critère Cl et b au critère C2, - minimiser le critère pondéré aC1 - b C2, - rechercher les valeurs des variables de la formulation qui minimisent le critère pondéré. Ce procédé est difficile à mettre en oeuvre pour les raisons suivantes: il oblige l'utilisateur à fixer la valeur des poids affectés aux critères par tâtonnements. En procédant ainsi on constate la plupart du temps que l'on optimise soit l'un ou soit l'autre critère et il est très difficile de trouver une solution qui réalise un compromis acceptable entre les deux critères, - il est long et fastidieux à mettre en oeuvre et n'est pas automatisable. EXPOSE DE L'INVENTION La présente invention a justement pour objet de remédier à ces inconvénients et notamment de fournir un procédé de fabrication d'un mélange ou d'une suite de mélanges de produits liquides, pâteux ou pulvérulents à partir de bases dont les caractéristiques sont connues, dont les formulations sont optimisées. Il trouve son application dans toutes les industries qui fabriquent des produits par mélange de produits de base, dont notamment l'industrie pétrolière pour la fabrication des carburants et des combustibles. A cette fin la présente invention propose un procédé de fabrication d'un mélange de produits comportant une étape de formulation optimale à partir d'un modèle du mélange, qui consiste à exécuter les opérations suivantes: - choisir au moins un premier et un deuxième critère élémentaires d'optimisation, - déterminer un intervalle de variation pour chacun des premier et deuxième critères élémentaires d'optimisation, chaque intervalle étant défini par les valeurs extrêmes du critère considéré, calculées au moyen du modèle du mélange et d'un algorithme d'optimisation, - définir une pluralité de points dans l'intervalle de variation du premier critère élémentaire d'optimisation ainsi déterminé, - calculer pour lesdits points, au moyen du modèle du mélange et d'un algorithme d'optimisation, les valeurs optimales du deuxième critère élémentaire d'optimisation, qui sont des éléments de l'ensemble des solutions efficientes, ainsi que les sensibilités dudit critère aux dits points, sélectionner parmi l'ensemble des solutions efficientes une solution préférée correspondant à un compromis prédéterminé entre les deux critères élémentaires d'optimisation, - calculer une formulation optimale correspondant à ladite solution préférée. Selon une autre caractéristique du procédé de l'invention la solution préférée est sélectionnée en choisissant un point parmi la pluralité de points définis dans l'intervalle de variation du premier critère d'optimisation. Selon une autre caractéristique du procédé de l'invention, la sélection de la solution préférée dans l'ensemble des solutions efficientes consiste à calculer la valeur du deuxième critère élémentaire d'optimisation à partir d'une valeur du premier critère élémentaire d'optimisation choisie par interpolation entre deux points consécutifs de la pluralité de points définis dans l'intervalle de variation dudit critère. Selon une autre caractéristique du procédé de l'invention, la sélection de la solution préférée dans l'ensemble des solutions efficientes consiste à choisir la solution pour laquelle la variation de sensibilité du deuxième critère d'optimisation est maximale. Selon une autre caractéristique du procédé de l'invention, la sélection de la solution préférée dans l'ensemble des solutions efficientes consiste à élaborer un critère pondéré d'optimisation avec les premier et deuxième critères d'optimisation affectés de poids libres puis à rechercher la solution qui optimise ledit critère pondéré d'optimisation. Selon une autre caractéristique du procédé de l'invention, la sélection de la solution préférée dans l'ensemble des solutions efficientes consiste à élaborer un critère pondéré d'optimisation avec les premier et deuxième critères d'optimisation affectés de poids libres respectivement égaux à a et b puis à fixer des valeurs aux poids libres telles que le rapport b sur a soit égal à la valeur du quotient de l'intervalle de variation du deuxième critère d'optimisation par l'intervalle de variation du premier critère élémentaire d'optimisation puis à rechercher la solution correspondante. Selon une autre caractéristique du procédé de l'invention le premier critère élémentaire d'optimisation est un critère de coût et le deuxième critère élémentaire d'optimisation est un critère de flexibilité. Selon une autre caractéristique le procédé de l'invention est utilisé pour fabriquer un carburant. Selon une dernière caractéristique le procédé de l'invention est utilisé pour fabriquer un suite de mélanges. BREVE DESCRIPTION DES DESSINS L'invention sera mieux comprise à l'aide de la description qui suit, d'un mode de réalisation donné à titre d'exemple et en référence aux dessins annexés dans lesquels: - la figure 1 est un diagramme qui représente le domaine d'évolution des deux critères élémentaires d'optimisation. - la figure 2 représente graphiquement la pluralité de points définis dans l'intervalle de variation du premier critère élémentaire d'optimisation et les éléments correspondants de l'ensemble des solutions efficientes - la figure 3 représente l'ensemble des solutions efficientes dans l'espace (C2,C1) d'un mélange de peintures. - La figure 4 représente l'ensemble chiffré des solutions efficientes dans l'espace (C2,C1) d'un mélange de mélange de peintures. - La figure 5 représente l'ensemble chiffré des solutions efficientes dans l'espace (C2,C1) d'un mélange de mélange de carburants. EXPOSE DETAILLE DE L'INVENTION D'une manière générale le procédé de l'invention est utilisé pour fabriquer des mélanges de produits en prenant en compte au moins deux critères d'optimisation. Le procédé de l'invention permet de trouver une solution au problème de la formulation optimale, qui assure un compromis réel entre les valeurs d'au moins deux critères d'optimisation, donc intermédiaire entre les solutions extrêmes. Pour la formulation d'un mélange le procédé de l'invention consiste dans une première étape à choisir: -un premier critère élémentaire d'optimisation, par exemple un critère de coût de fabrication du mélange, le coût des bases étant connu, -un deuxième critère élémentaire d'optimisation, par exemple un critère de flexibilité qui représente la possibilité de réaliser ultérieurement les mélanges les plus divers à partir des bases disponibles. Puis à déterminer un intervalle de variation pour les critères élémentaires d'optimisation C1 et C2. Dans le cas considéré on cherche la valeur minimale de Cl soit Clmin qui correspond au coût minimal de fabrication du mélange et la valeur maximale de C2 soit C2max qui correspond à la flexibilité maximale, en calculant les formulations qui correspondent à ces extrêmes au moyen du modèle du mélange qui représente la relation mathématique entre les caractéristiques des bases et celles du mélange. A la valeur Clmin correspond une solution qui permet de calculer une valeur C2min du critère C2 et à la valeur C2max correspond une valeur Clmax du critère C1. La figure 1 illustre sous forme graphique les relations entre ces différentes valeurs. La surface circonscrite par la courbe 1 représente les valeurs de C2 en fonction de celles de Cl. La valeur minimale de Cl correspond à l'abscisse du point P. déterminée au moyen du modèle et l'ordonnée du point P correspond à la valeur C2min. La valeur maximale de C2 correspond à l'ordonnée du point R et l'abscisse de ce point R correspond à la valeur Clmax. L'intervalle de variation du critère Cl est l'intervalle 2 compris entre les valeurs Clmin et Clmax. L'intervalle de variation du critère C2 est l'intervalle 3 compris entre les valeurs C2min et C2max. On définit ensuite sur l'intervalle de variation de Cl une pluralité de points, par exemple 7, répartis arbitrairement pour couvrir la totalité de l'intervalle 2. Cette pluralité de points est représentée sur la figure 2 par les points 6a à 6g. Pour chacun de ces points 6a à 6g, qui correspondent à des valeurs de Cl, on calcule les valeurs maximales de C2 et les sensibilités correspondantes au moyen du modèle du mélange et d'un algorithme d'optimisation. Si on porte sur le diagramme de la figure 2 les points d'abscisses égales à celles des points 6a à 6g et d'ordonnées égales aux valeurs maximales correspondantes de C2 on obtient la pluralité de points 7a à 7g, qui sont des éléments de l'ensemble des solutions efficientes. La sensibilité en un des points 7a à 7g est égale par définition à la pente p de la courbe 4 obtenue en joignant les points 7a à 7g. Cette courbe 4 est une approximation de l'ensemble des solutions efficientes. A une variation donnée AC1 correspond une variation AC2, telle que AC2 = p.AC1. Une pente p élevée en un point de la courbe 4 signifie qu'en ce point, pour une variation faible du coût de fabrication on peut obtenir une forte augmentation de la flexibilité. Inversement dans les zones de faible pente de la courbe Cl une augmentation importante du coût de fabrication n'entraîne qu'une faible augmentation de la flexibilité. Pour sélectionner la solution préférée c'est à dire la solution considérée par le fabricant du mélange comme le meilleur compromis entre le coût de fabrication et la flexibilité, on fixe une valeur à Cl en choisissant un des points parmi la suite 6a à 6g. L'ordonnée du point choisi donne la valeur de C2. Un des avantages du procédé de l'invention est son caractère automatisable qui élimine les tâtonnements. De plus après calcul des coordonnées des points 7a à 7g la courbe 4 peut être visualisée sur un écran. Ainsi l'opérateur qui doit réaliser le mélange dispose d'un moyen extrêmement efficace pour choisir le compromis entre les deux critères élémentaires d'optimisation qu'il considère comme le meilleur compte tenu des contraintes de fabrication du moment et des futurs mélanges dont il a connaissance. EXEMPLE 1 Le mélange pour lequel une formulation optimale est recherchée est un mélange de peintures de volume et de couleur donnés, à réaliser à partir de trois bases de couleurs connues disponibles dans des quantités connues. Dans la description qui va suivre les données et variables qui seront utilisées sont définies comme suit: *Données et variables relatives aux bases: - j est une variable indicielle qui varie de 1 à 3, - Bj représente les bases disponibles, - Vjo représente le volume initial disponible de la base Bj, - Dj représente le coût de la base Bj par unité de volume, - Rj représente l'indice de couleur de la base Bj, *Données relatives au mélange à réaliser: - Vm représente le volume du mélange à réaliser, - Rm représente l'indice de couleur du mélange à réaliser, - Vj représente le volume de la base Bj, utilisé pour réaliser le mélange. Le but du procédé de l'invention est de trouver une formulation qui assure un compromis réel entre les valeurs d'un premier et d'un deuxième critère élémentaires d'optimisation qu'on appelle respectivement Cl et C2 et qui sont définis comme suit: - Cl est un critère de coût de fabrication du mélange défini par la formule suivante: Cl = Vj. Dj Ci est égal à la somme des coûts des bases utilisées - C2 est un critère de flexibilité défini par la formules suivante: C2 = [ Z (Vjo - Vj)(Vio - Vi) Ri]1/2 l donnée. Le symbole I signifie: valeur absolue du déterminant. En développant cette formule on obtient: C2 = [(Vlo-Vl)(V2o- V2)IR2-Rl1 +(Vlo-Vl)(V3o-V3) IR3-Rl1 +(V2o-V2)(V3o-V3) IR3-RlI] 1/2 Dans le présent exemple le modèle du mélange est la relation mathématique qui relie l'indice de couleur Rm du mélange aux indices de couleur Rj des bases. On a choisi le modèle linéaire suivant: Vm.Rm = Vj.Rj et Vm = Vj J j Avec ce modèle et ces critères les solutions aux problèmes d'optimisation qui devront être résolus dans les étapes suivantes, sont toutes uniques. L'étape suivante du procédé de l'invention consiste à déterminer les intervalles de variations de Cl et C2. Pour cela, comme on cherche à minimiser le coût du mélange et à maximiser la flexibilité pour se ménager la possibilité de réaliser d'autres mélanges dans les meilleures conditions, on recherche la valeur minimale de Ci qu'on appellera Clmin puis la valeur maximale de C2 qu'on appellera C2max. Pour calculer Clmin on cherche tout d'abord, au moyen d'un programme d'optimisation linéaire mettant en oeuvre la méthode dite du Simplex, les valeurs de Vl, V2, V3 qui résolvent l'équation min Cl (V1,V2,V3) sous les trois contraintes suivantes: 1- L'indice de couleur du mélange Rm doit être compris entre les valeurs maximales et minimales des indices de couleur des bases Rj: Rmin Rm = (Z Vj.Rj)/Vm) base Bj. Le programme linéaire d'optimisation donne les solutions V*i, V*2, V*3. Avec ces valeurs on calcule la valeur Clmin de Ci et la valeur correspondante C2min de C2 soit: Clmin = V*1 D1+ V*2 D2+ V*3 D3 C2min = [(Vlo-V*l)(V2o-V*2) IR2-Rll +(Vlo-V*l)(V3o-V*3) IR3-Rl1 +(V2oV*l)(V3o-V*3) IR3-Rll] 1/2 Pour calculer C2max on cherche tout d'abord, au moyen de la méthode de programmation non linéaire SQP (Sequential Quadratique Programming), les valeurs de V1, V2, V3 qui résolvent l'équation max C2 (V1,V2,V3) sous les trois contraintes exprimées ci-dessus. Soit V01, V 2, V 3 les solutions données par le programme d'optimisation non linéaire avec lesquelles on calcule C2max et Clmax. C2max = [(Vlo-V l)(V2o-V 2)IR2-Rll +(Vlo-V l)(V3o- V 3) IR3-Rll +(V20o-VO2) (V30o-VO3) R3-Rl] 1/2 Clmax = V 01.D1 + V 2 D2 + V 3 D3 Les intervalles de variation de Cl et C2 sont: Intervalle CI: [Clmin, Clmax] Intervalle C2: [C2min, C2max] Clmin, Clmax, C2min, C2max ayant les valeurs calculées précédemment. On définit ensuite 8 points sur l'intervalle de variation de Ci dont les extrémités de manière à déterminer sensiblement 7 intervalles égaux à (Clmax - C2min)/7. Ces points sont repérés 6a, 6b, 6c, 6d, 6e, 6f, 6g et 6h, sur l'axe C1 du graphique de la figure 3, ils ont pour abscisses respectives les valeurs xi avec i=l à 8. A partir de chacune des valeurs xi du critère Ci, on calcule la valeur correspondante de max C2(V1,V2,V3), par programmation non linéaire sous les trois contraintes exprimées ci-dessus plus une quatrième: xi = I Dj.Vj Pour chaque point d'abscisse xi le programme d'optimisation non linéaire donne les solutions V*li, V*2i, V*3i ainsi que la valeur de la sensibilité de C2 au point d'abscisse considéré. A partir de ces solutions on calcule les valeurs correspondantes de C2, soit yi yi = t(Vlo-V*li)(V2o-V*2i) IR2-Rll +(Vlo-V*li)(V3o-V*3i) IR3- R1I +(V2o-V*2i)(V3o-V*3i) IR3-R2 1] /2 Remarque: la valeur de C2 pour x8 soit y8 a été calculée précédemment, elle est égale à C2max. Sur la figure 3 sont représentés les points 7a, 7b, 7c, 7d, 7e, 7f, 7g et 7h, dont les coordonnées sont (xl,yl), (x2,y2), (x3,y3), (x4,y4), (x5,y5), (x6,y6), (x7,y7) et (x8,y8). On a également tracée sur la figure 3 la courbe 4 qui passe par les points 7a à 7h qui représente l'ensemble des solutions efficientes. La pente en un point de la courbe 4 représente la sensibilité du critère C2 par rapport aux variations du critère C1 en ce point. Les valeurs des pentes aux points 7a à 7h sont données par le multiplicateur de Lagrange, lui même calculé par le programme d'optimisation non linéaire en même temps que les solutions. On désigne ces pentes par pl, p2 à p8 pour les points respectivement 7a à 7h. Pour terminer on sélectionne la solution préférée, c'est à dire celle qui réalise le meilleur compromis entre les valeurs du critère de coût C1 et du critère ce flexibilité C2 prédéterminée en fonction des besoins de la production. Pour cela on recherche automatiquement le point de la courbe 4 pour lequel la variation de sensibilité est maximale, au moyen d'un algorithme qui calcule le maximum de la fonction IPi - Pill avec i = 2 à 8 Soit k la valeur de i ainsi déterminée. On calcule p* avec cette valeur de k, par la formule p* = - (pk + pk+l) On résout le problème d'optimisation min p*.C1 - C2 par un programme non linéaire avec les contraintes suivantes: Vm = Z Vj Rm - Rj.Vj Vm j 0 V*3p, qui constituent la solution préférée. Avec ces valeurs on calcule les valeurs optimales des critères au compromis: soit Cl*=Cl(V*lp,V*2p,V*3p) et C2*=C2(V*lp,V*2p,V*3p) La formulation optimale du mélange de volume Vm de couleur Rm est la suivante: Volume de la base Bl: V*lp Volume de la base B2: V*2p Volume de la base B3: V*3p Une autre manière particulièrement intéressante pour sélectionner la solution préférée consiste: - à élaborer un critère z pondéré d'optimisation sous la forme de l'expression suivante: z = aC1 + bC2 dans laquelle: Cl représente le premier critère élémentaire d'optimisation c'est à dire le critère de coût de fabrication, C2 représente le deuxième critère élémentaire d'optimisation c'est à dire le critère de flexibilité. 30. a représente le poids affecté au critère CI dans le critère pondéré z, b représente le poids affecté au critère C2 dans le critère pondéré z. C2 max- C2 min On choisit b = -1 et a = Cl maxCl min a représente la sensibilité moyenne. - à chercher au moyen d'une méthode de programmation non linéaire les valeurs de Vl, V2, V3 qui résolvent l'équation: min z sous les trois contraintes exprimées précédemment et rappelées ci-après: 1) R min précédemment. Les valeurs de V1, V2, V3 ainsi déterminées constituent la solution préférée. Ainsi pour fabriquer un mélange de peinture grise de volume Vm = 500 m3 d'indice de couleur Rm = 50, à partir de trois bases de couleurs respectivement noire, blanche et grise dont les données et caractéristiques sont les suivantes: TABLEAU 1 BASES Couleur Indice: j Volume Coût:Dj Indice de :i disponible: en francs couleur: Rj Vjo en m3 Noire 1 500 95 1 Blanche 2 500 105 99 Grise 3 500 101 50 les intervalles de variation de Ci et de C2 sont les suivants: Intervalle Cl: [Clmin = 50 000 F. Clmax = 50 500 F] Intervalle C2: [C2min = 4,2866, C2max = 4,9497] Les résultats des calculs d'optimisation sont les suivants: TABLEAU 2 Points sur Solutions efficientes Intervalle CI Repères Abscisses Repères Ordonnées Sensibilités xi yi pi 6a 50 000 7a 4,2866 0, 002597 6b 50 071 7b 4,4721 0,002347 6c 50 142 7c 4,6233 0,001892 6d 50 214 7d 4,7434 0,001475 6e 50 285 7e 4,8347 0,001085 6f 50 357 7f 4,8989 0,000714 F! 6g 50 428 7g 4,9371 10, 000354 6h 50 500 7h 4,9497 0,000000 Remarque: dans cet exemple, on choisit Rm = Rmin = Rmax traduisant le fait que l'on veut fabriquer un mélange dont l'indice de couleur est déterminé. Si on voulait fabriquer un mélange dont l'indice de couleur serait compris entre les deux valeurs Rmin et Rmax, il suffirait de remplacer la contrainte d'égalité par deux contraintes d'inégalité. Les valeurs xi et yi du Tableau 2 sont les coordonnées des points 7a à 7h représentés sur la figure 4. Sur cette figure 4 l'axe des abscisses représente les valeurs du critère Ci et l'axe des ordonnées les valeurs du critère C2. La courbe 4 obtenue en joignant les points 7a à 7h représente l'ensemble des solutions efficientes. La solution préférée qui réalise le meilleur compromis entre les valeurs du critère de coût Cl et du critère de flexibilité C2, sélectionnée en recherchant le point de la courbe 4 pour lequel la variation de sensibilité est maximale correspond à: CI* = 50 106 francs et C2* = 4,5498 point repéré S sur la courbe 4. la sensibilité p* pour ces valeurs de Cl et C2 étant égale à 0,002120. La formulation optimale du mélange est: Volume de peinture noire: 197 m3 soit 39,37% de Vm Volume de peinture blanche: 197 m3 soit 39,37% de Vm Volume de peinture grise: 106 m3 soit 21,25% de Vm pour fabriquer 500 m3 d'un mélange dont l'indice de couleur est égal à 50. EXEMPLE 2 Le mélange à fabriquer pour lequel une formulation optimale est recherchée est un carburant dont les caractéristiques doivent être conformes à des spécifications commerciales. Pour réaliser ce mélange on dispose dans des bacs de stockage de sept bases en quantités connues dont les caractéristiques physiques et physico- chimiques ont été mesurées à partir d'échantillons prélevés dans les bacs de stockage et/ou au moyen d'analyseurs en ligne pendant leur fabrication. Dans la description qui va suivre les données et variables qui seront utilisées sont définies comme suit: * Données et variables relatives aux bases: - j est une variable indicielle qui varie de 1 à 7, - Bj représente les bases disponibles, - Vjo représente le volume initial disponible de la base Bj, - Dj représente le coût de la base Bj par unité de volume, - Rij représente la caractéristique i de la base Bj, i variant de 1 à 6 - TK18, TK19, TK21 et TK04 sont des repères de bacs de stockage des bases B1,B2, B3 et B4 - Ref, FCC, Butane caractérisent l'origine des bases B5, B6 et B7 * Données relatives au mélange à réaliser: - Vm représente le volume du mélange à réaliser, - Rim représente la caractéristique i du mélange à réaliser i variant de 1 à 6, - Vj représente le volume de la base Bj, utilisé pour réaliser le mélange. * Caractéristiques des bases et du mélange - MV signifie masse volumique RON signifie indice d'octane recherche - MON signifie indice d'octane moteur - RVP signifie tension de vapeur - %D 70 o signifie pour-cent distillé à 70 C - %D 100 signifie pour-cent distillé à 100 C Le but du procédé de l'invention est de fabriquer un mélange selon une formulation qui assure un compromis réel entre les valeurs optimales de deux critères élémentaires d'optimisation qu'on appelle respectivement Cl et C2 et qui sont définis comme suit: - C1 est un critère de coût de fabrication du mélange défini par la formule suivante: Cl = Z Vj. Dj qui exprime que Ci est égal à la somme des coûts des bases utilisées - C2 est un critère de flexibilité basé sur le volume, le RON et le RVP, défini par la formules suivante: 1/3 F 1 RONi RVPi C2 = E (Vio- Vi)(Vjo-Vj)(Vko-Vk) 1 RONj RVPj 1 varient de 1 à 7. Vio, Vjo et Vko sont les volumes initiaux disponibles respectivement des bases i, j et k. Dans le présent exemple le modèle linéaire du mélange est la relation mathématique qui relie les caractéristiques Rim du mélange aux caractéristiques Rij des bases ou à des indices représentatifs de ces caractéristiques. TABLEAU 3 Bases Caractéristiques des bases Bj Rij Indice: Origine Volume MV: RON: MON: j disponible Rlj R2j R3j Vjo en m3 t/m3 1 TK18 8 000 0,713 95,8 93,1 2 TKl9 2 500 0,769 95,0 85,0 3 TK21 3 000 0,746 118 101 4 TK04 5 000 0,649 70,4 69,3 Ref 21 700 0,815 102 92,6 6 FCC 37 500 0,741 92,7 81,3 7 Butane 3 000 0,584 97,1 88,3 7 Butane 300 |,5497, 1 88,3 TABLEAU 4 Bases Caractéristiques des bases Rij Indice: Origine Volume RVP: D 70: %D 1000: j disponible R4j R5j R6j en Bar 1 TK18 8 000 0,4 10 50 2 TK19 2 500 0,83 20 48 3 TK21 3 000 0,53 100 100 4 TK04 i 5 000 0,94 95 100 I I. I Ref 21 700 0,25 6 31 6 FCC 37 500 0,7 30 45 7 Butane 3 000 4,73 100 100 TABLEAU 5 Bases Indice: Origine Coût j en DM 1 TK18 359 2 TK19 239 3 TK21 499 4 TKO4 303 Ref 331 6 FCC 318 7 Butane 225 TABLEAU 6 Caractéristiques du mélange carburant à formuler Volume MV RON MON en m3 en t/m3 R2m R3m Rlm Valeur 12 541 0,726 91,6 83,1 minimale Valeur 12 541 0,779 93 86 maximale TABLEAU 7 Caractéristiques du mélange à formuler (suite) Volume RVP %D 70 %D 100 en m3 en Bar R5m R6m i R4m Valeur 12 541 0,46 15 40 minimale Valeur 12 541 0,68 42 65 maximale En déroulant les différentes étapes du procédé de l'invention on obtient les résultats suivants: Les intervalles devariation de Ci et C2 sont les suivants: Intervalle Cl: [3 875 430, 4 066 472] Intervalle C2: [41,0, 50,3] La courbe 4 de la figure 5 représente l'ensemble des solutions efficientes dans l'espace Ci, C2. La formulation optimale et les caractéristiques du mélange formulé qui correspondent à la solution préférée, c'est à dire celle qui réalise le meilleur compromis entre les valeurs des deux critères d'optimisation, sont données ci-après dans les tableaux 8 et 9. Sur la figure 5 le point S correspond à la solution préférée pour laquelle la variation de sensibilité du critère C2 par rapport au critère Cl est maximale. En considérant la figure 5 on constate que entre les point 7a et 7b à une faible variation du critère de coût C1 correspond une très forte augmentation de la valeur du critère de flexibilité C2. Inversement entre les points 7b et 7h, à une forte augmentation du coût correspond une très faible augmentation de la valeur du critère de flexibilité. TABLEAU 8 Bases Formulation optimale du mélange Indice: Origine Volume Volume j de la base j de la base j en %Vm en m3 1 TK18 13,65 1711,3 2 TKl9 10,4 1304 3 TK21 0 0 4 TK04 1,62 203,7 Ref 0 0 6 FCC 74,33 9322,1 7 Butane 0 0 TABLEAU 9 Caractéristiques du mélange formulé Volume Coût MV RON MON RVP %D 70 %D 1000 en m3 en MDM en t/m3 12 541 3,952 0,7386 93 83,1 0,676 27 47 Les exemples d'application du procédé de l'invention donnés ci-dessus ne sont pas limitatifs, au contraire le procédé de l'invention peut être utilisé pour optimiser la formulation de mélanges de produits très divers tels que des poudres, des liquides ou des pâtes dans des domaines très variés de l'industrie. Le nombre de critères d'optimisation et le nombre des caractéristiques pris en compte par le procédé de l'invention ne sont limités pas limités à deux comme dans les exemples. Il peuvent éventuellement être limités par les programmes d'optimisation linéaires et non linéaires utilisés. REVENDTCATTIONS 1- Procédé de fabrication d'un mélange de produits comportant une étape de formulation optimale à partir d'un modèle du mélange caractérisé en ce que ladite étape consiste à exécuter les opérations suivantes: - choisir au moins un premier et un deuxième critère élémentaires d'optimisation, déterminer un intervalle de variation pour chacun des premier et deuxième critères élémentaires d'optimisation chaque intervalle étant défini par les valeurs extremes du critère considéré, calculées au moyen du modèle du mélange et d'un algorithme d'optimisation, - définir une pluralité de points dans l'intervalle de variation du premier critère élémentaire d'optimisation ainsi déterminé, - calculer pour lesdits points, au moyen du modèle du mélange et d'un algorithme d'optimisation, les valeurs optimales du deuxième critère élémentaire d'optimisation, qui sont des éléments de l'ensemble des solutions efficientes, ainsi que les sensibilités dudit critère aux dits points, - sélectionner parmi l'ensemble des solutions efficientes une solution préférée correspondant à un compromis prédéterminé entre les deux critères élémentaires d'optimisation, - calculer la formulation optimale correspondant à ladite solution préférée. 2- Procédé selon la revendication 1 caractérisé en ce que la solution préférée est sélectionnée en choisissant un point parmi la pluralité de points définis dans l'intervalle de variation du premier critère d'optimisation. 3- Procédé selon la revendication 1 caractérisé en ce que la sélection de la solution préférée dans l'ensemble des solutions efficientes consiste à calculer la valeur du deuxième critère élémentaire d'optimisation à partir d'une valeur du premier critère élémentaire d'optimisation choisie par interpolation entre deux points consécutifs de la pluralité de points définis dans l'intervalle de variation dudit critère. 4- Procédé selon la revendication 1 caractérisé en ce que la sélection de la solution préférée dans l'ensemble des solutions efficientes consiste à choisir la solution pour laquelle la variation de sensibilité du deuxième critère d'optimisation est maximale. - Procédé selon la revendication 1 caractérisé en ce que la sélection de la solution préférée dans l'ensemble des solutions efficientes consiste à élaborer un critère pondéré d'optimisation avec les premier et deuxième critères d'optimisation affectés de poids libres puis à rechercher la solution qui optimise ledit critère pondéré d'optimisation. 6- Procédé selon la revendication 1 caractérisé en ce que la sélection de la solution préférée dans l'ensemble des solutions efficientes consiste à élaborer un critère pondéré d'optimisation avec les premier et deuxième critères d'optimisation affectés de poids libres respectivement égaux à a et b puis à fixer des valeurs aux poids libres telles que le rapport b sur a soit égal à la valeur du quotient de l'intervalle de variation du deuxième critère d'optimisation par l'intervalle de variation du premier critère élémentaire d'optimisation puis à rechercher la solution correspondante. 7- Procédé selon une des revendications 1 à 6 caractérisé en ce que le premier critère élémentaire d'optimisation est un critère de coût et le deuxième critère élémentaire d'optimisation est un critère de flexibilité. 8- Application du procédé selon une des revendications 1 à 7 à la fabrication d'un carburant. 9- Application du procédé selon une des revendications 1 à 8 à la fabrication d'une suite de mélanges.