Brevet FR2562908A1 — publié le 19851018

La présente invention concerne les appareils de cris- tallisation par évaporation à marche continue constitués par une cuve fermée disposée horizontalement et divisée par des cloisons verticales en plusieurs compartiments munis d'éléments, tels que des plaques creuses ou des tubes, chauffés par de la vapeur, et qui sont traversés successi- vement par le mélange de cristaux et de liqueur-mère, Une solution sucrée sous-saturée, appelée liqueur standard ou égout, est introduite dans chacun des compartiments et un magma formé d'un mélange de petits cristaux et de liqueur- mère est introduit dans le premier compartiment; la masse cuite produite est extraite du dernier compartiment. Pour automatiser la marche de ces appareils, on a prévu de régler l'alimentation en solution sous-saturée ou la pression de la vapeur de chauffage de chaque compar- timent ou groupe de compartiments de façon à obtenir une masse cuite ayant la composition voulue; on a également proposé de controler la sursaturation de la liqueur-mère dans un ou plusieurs compartiments de cristallisation. Mais, compte tenu des temps de séjour très longs et de la mauvaise représentativité des paramètres mesurés, les systèmes de régulation classiques ne permettent pas une conduite fiable en régime transitoire. Le but de la présente invention est d'optimiser la conduite de l'appareil de cristallisation en prenant en compte, pour la détermination des valenrs de consigne des grandeurs réglées, les caractéristiques de la solution sucrée et du magma alimentant l'appareil de cristallisation et les caractéristiques désirées pour la masse cuite pro- duite, ainsi que les paramètres caractérisant le fonction- nement de l'appareil (pression de calandre, débit de va- peur d'agitation, niveau de masse cuite, etc...). Cette optimisation permet d'augmenter la production pour une qualité donnée de sucre, de réduire les pertes dans les mélasses et de diminuer la consommation d'énergie. Le procédé objet de l'invention consiste à effectuer périodiquement les opérations suivantes: on détermine pour chaque compartiment d'indice 1 la valeur optimale -2- d'un paramètre choisi dans un premier groupe de cinq paramè- tres constitué par le débit d'alimentation en solution su- crée des compartiments rapporté au débit total de solution sucrée à traiter (YEG1), le débit d'eau évaporée dans les compartiments rapporté au débit total d'eau évaporée dans l'appareil de cristallisation (YEV1), le débit de masse cuite sortant des compartiments rapporté au débit total de masse cuite produite (YMC1), le brix (teneur en matières sèches) de la masse cuite dans les compartiments (BMC1) et la teneur en cristaux de sucre de la masse cuite dans les compartiments (RDT1) en utilisant une loi préétablie donnant la valeur dudit paramètre dans les dits compartiments en fonction du brix de la solution sucrée à traiter (BEGT), de la pression de vapeur dans la calandre de l'appareil de cris- tallisation (PCAL)et du rapport R = (QMAG x RDTMAG) / (QEGT x BEGT) avec: QMAG: débit de magma RDTMAG: teneur en cristaux du magma QEGT: débit total de solution sucrée à traiter BEGT: brix de la solution sucrée à traiter, à calculer les valeurs optimales de quatre paramètres d'un second groupe de six paramètres constitué par le débit total de solution sucrée à traiter (QEET), le débit total d'eau évaporée (QEV), le débit total de masse cuite produite (QMC), le brix de la masse cuite produite (BMC), la pression de la vapeur de chauffage (PVAP) et la teneur en cristaux de sucre de la masse cuite produite (RDT), ce second groupe comprenant deux paramètres de conduite (QEGT et PVAP) et quatre paramètres d'état résultants (QEV, QMC, BMC et RDT), en utilisant pour ce calcul les valeurs optimales détermi- nées pour chaque compartiment du paramètre du premier groupe, et des valeurs imposées de deux autres paramètres du second groupe> dont l'un est un paramètre de conduite et l'autre un paramètre d'état résultant et éventuellement un terme correctif dont la valeur est fonction de l'écart entre la valeur imposée dudit paramètre d'état résultant et une valeur mesurée de celui-ci-et on règle l'autre paramètre de conduite de ce second groupe pour le maintenir -3- à une valeur de consigne égale à la valeur optimale calculée. Les paramètres de conduite sont des paramètres sur lesquels on peut agir directement par opposition aux para- mètres d'état résultantsdont les variations résultent d'une modification des paramètres de conduite. Les lois donnant les variations des paramètres du premier groupe peuvent aussi être fonction du débit total de solution à traiter (QEGT) ou du débit total d'eau éva- porée (QEV) ou du débit total de masse cuite (QMC) suivant le choix des paramètres à optimiser. Elles peuvent être présentées sous forme de relations mathématiques ou d'un fichier de courbes ou de tableaux. Pour le calcul des valeurs optimales des paramètres du second groupe, on pourra utiliser, à la place du para- mètre d'état résultant dont les valeurs imposées et mesurées doivent être comparées, une image physique ou chimique dudit paramètre; par exemple le brix de la masse cuite produite (BMC) pourra être remplacé par la densité ou la compacité ou la température de la masse cuite, la teneur en cristaux de sucre de la masse cuite produite (RDT) pour- ra être remplacé par la compacité de la masse cuite et le débit total d'eau évaporée (QEV) pourra être remplacé par le débit d'eau condensée dans le faisceau de chauffage de l'appareil de cristallisation. Il est bien entendu que dans la description et les revendications qui suivent toute désignation de l'un quelconque de ces paramètres comprendra le paramètre lui-même et ses grandeurs de remplacement. Pour calculer les valeurs optimales des quatre para- mètres du second groupe dont un paramètre de conduite et trois paramètres d'état résultant on effectue un bilan mas- sique global, un bilan massique matières sèches, un bilan thermique d'eau évaporée et un bilan de cristallisation. Ces calculs sont effectués à intervalles de temps réguliers, pouvant varier de 1 à 5 minutes, par un système de conduite automatisée comportant un ordinateur de supervision qui re- çoit des informations, en provenance de différents appareils de mesure, sur certaines caractéristiques de la solution sucrée, du magma et de la masse cuite et sur les paramètres de marche de l'appareil de cristallisation ainsi que des données extérieures introduites par exemple au moyen d'un clavier, notamment les valeurs imposées de certains para- mètres. Comme paramètres du second groupe dont les valeurs sont imposées on pourra, par exemple, retenir le débit total de solution sucrée à traiter (QEGT) ou la pression de vapeur de chauffage (PVAP) et l'un quelconque des para- mètres d'état résultants (QEV, QMC, BMC et RDT). Des résultats particulièrement intéressants peuvent être obtenus en choisissant comme paramètre du premier groupe, dont la valeur optimale est déterminée à partir d'un relation mathématique ou d'un fichier de courbes ou de tableaux, le débit d'alimentation en solution sucrée des compartiments (YEG1) et comme paramètre de conduite du second groupe, dont la valeur optimale est calculée, la pression de vapeur de chauffage (PVAP), les deux paramètres dont les valeurs sont imposées, en fonction de la gestion des jets de cristallisation en amont et en aval, étant le débit total de solution sucrée (QEGT) et le brix de la masse cuite produite (BMC). C'est ce mode de mise en oeuvre qui est décrit ci-après, A titre d'exemple non limitatif, en se référant aux dessins l'accompagnant et sur lesquels: La figure 1 est le schéma d'un système de conduite automa- tisée d'un appareil de cristallisation pour la mise en oeuvre du procédé de l'invention, La figure 2 est l'organigramme du procédé, La figure 3 est une courbe de répartition du débit total de solution sucrée entre les compartiments, et La figure 4 est une courbe des variations du débit de vapeur d'agitation en fonction du débit total de solution sucrée à traiter. Sur la figure 1, le numéro de référence 10 désigne un appareil de cristallisation par évaporation A marche continue constitué par une cuve fermée disposée horizon- talement et divisée par des cloisons verticales en plusieurs compartiments 1, 2, 3, 4,..., n-2, n-1, n. Ces cloisons n'atteignent pas le sommet de la cuve de sorte que tous les compartiments communiquent entre eux à la partie supérieure de la cuve. Des ouvertures sont également prévues dans les cloisons pour permettre le passage de la masse cuite d'un compartiment à l'autre. Dans la partie inférieure de la cuve est placé un faisceau de chauffage 12 commun à tous les compartiments et formé de tubes ou de plaques dans lesquels on fait circuler de la vapeur amenée par une conduite 14. Une solution sucrée, qui sera appelée égout dans ce qui suit, est distribuée dans les différents compartiments par des tuyauteries 16; cet égout provient d'un bac 18 auquel les tuyauteries 16 sont reliées par une conduite 20. On peut prévoir une tuyauterie d'alimentation par compar- timent ou unetuyauterie commune à un groupe de compartiments. Dans ce cas, la répartition de débit entre les compartiments du groupe pourra être faite par exemple au moyen d'un dis- positif 52 du type faisant l'objet du brevet 84.00101. Un magma constitué par un mélange de petits cris- taux constituant des germes de cristallisation et de li- queur-mère est introduit dans le premier compartiment au moyen d'une pompe à débit variable 22 prélevant le magma sur une capacité 24. La masse cuite produite est extraite du dernier com- partiment au moyen d'une pompe à débit variable 26. De la vapeur est injectée à la partie inférieure de chaque compartiment, pour agiter la masse cuite et favoriser les échanges thermiques; cette vapeur amenée par une con- duite 28 est distribuée à des injecteurs 30 placés dans le bas des compartiments. Le numéro de référence 32 désigne un système de con- duite automatisée se composant d'un ordinateur de super- vision couplé à un système de conduite de processus com- prenant les régulateurs numériques et l'instrumentation ou à des chaines de régulation analogiques classiques. En variante l'ordinateur peut assurer les fonctions de régu- lation et communiquer avec l'instrumentation par des conver- tisseurs digital-analogique et analogique-digital. Le système reçoit des informations numériques et analogiques, - 6 - les traite et émet les consignes qu'il a optimisées vers les actionneurs; vannes, moteurs, etc... Pour simplifier la figure 1, les différents éléments constitutifs de ce système ont été répartis sur le dessin. Le système 32 reçoit en permanence des informations sur la température, le débit total et le brix de l'égout qui sont mesurées par des instruments 34, 36 et 38, res- pectivement et des informations sur le brix de la masse cuite fournies par un instrument 54. Il reçoit en outre des données externes: résultats d'analyses de laboratoire, caractéristiques imposées par l'exploitant, etc... qui sont introduites par exemple au moyen d'un clavier 39. Il a également en mémoire les caractéristiques technologiques de l'appareil 10: dimensions et volume utile de la cuve, structure du faisceau de chauffage, dimension des tubes du faisceau, etc... Le système comprend en outre un régulateur de la pression de la vapeur de chauffage 40, un régulateur 42 de la pression de vapeur dans la calandre de l'appareil, des régulateurs de débit partiel d'égout 44, un régulateur de débit de vapeur d'agitation 46, un régulateur de débit de magma 48 commandant la pompe 22 et un régulateur de niveau de masse cuite 50 commandant la pompe 26. Conformément à la présente invention, on détermine à intervalles de temps réguliers les valeurs optimales des débits d'alimentation en égout de chaque compartiment ou groupe de compartiments, du débit de magma admis dans le premier compartiment, de la pression de la vapeur de chauf- fage et du débit de vapeur d'agitation injecté dans les compartiments et on les utilise comme valeurs de consigne dans les régulateurs correspondants. Ces intervalles de temps pourront être compris, par exemple, entre une et cinq minutes. Les valeurs optimales des débits d'alimenta- tion QEG1 des compartiments sont déterminées au moyen d'un fichier donnant sous forme de courbes ou de tableaux la répartition en pourcentage du débit total QEGT entre les différents compartiments. Cette répartition dépend de différents paramètres, notamment de la pression de vapeur -7dans la calandre de l'appareil PCAL, du débit total d'égout à traiter QEGT, de son brix BEGT et du rapport R = (QMAG x RDTMAG)-/ (QEGT x BEGT), QMAG étant le débit de magma introduit dans le ler compar- timent et RDTMAG étant la teneur en cristaux du magma. Le fichier cnntient plusieurs courbes ou tableaux correspondant à différentes valeurs de ces paramètres et le système de conduite automatisée sélectionne la courbe ou le tableau à utiliser en fonction des valeurs mesurées ou im- posées de ces paramètres. La fig. 3 montre une telle courbe de répartition du débit total d'égout.On a porté en abscisse les numéros des compartiments et en ordonnées les débits cumulés en pourcentage du débit total. La valeur optimale du débit de magma est calculée à l'aide de la formule R = (QMAG x RDTMAG) / (QEGT x BEGT) à partir des valeurs mesurées de QEGT et BEGT, de la valeur de RDTMAG, déterminée à partir d'analyses en laboratoire du magma, et de la valeur de R fixée par l'exploitant. La valeur optimale du débit total de vapeur d'agita- tion injecté dans les compartiments est déterminée au moyen d'un second fichier donnant sous forme de courbes ou de tableaux le débit de vapeur QINJ en fonction du débit total d'égout à traiter QEGT ou du débit de masse cuite produite QMC ou du débit total d'eau évaporée QEV. La relation entre QINJ et QEGT, QMC ou QEV, dépend de différents paramètres, notamment du brix de l'égout à traiter, de la pression de la calandre et du rapport R et ce fichier contient plusieurs courbes ou tableaux correspondant à différentes valeurs de ces paramètres. La figure 4 montre une courbe des variations du débit de vapeur d'agitation QINJ en fonction du débit d'égout QEGT. On pourrait aussi envisager de régler les débits de vapeur injectés compartiment par compartiment; dans ce cas, le fichier contiendrait autant de jeux de courbes ou de tableaux qu'il y a de compartiments. La valeur optimale de la pression de la vapeur de chauffage PVAP est calculée par le système 32 suivant un procédé schématisé sur la figure 2. Le système de conduite automatisée effectue pour chaque compartiment un bilan de _ 8 _ cristallisation et un bilan thermique d'eau évaporée en utilisant les valeurs mesurées ou déterminées comme indiqué plus haut des différentes grandeurs caractérisant l'égout à traiter et le magma (débit, brix, température, etc...), ainsi que celles des différents paramètres caractérisant le fonctionnement de l'appareil (débits d'alimentation en égout, débit de vapeur d'agitation, pression de vapeur dans la calandre de l'appareil, niveau de masse cuite dans l'ap- pareil, etc...). Pour effectuer les bilans massiques global, de matières sèches et de cristallisation le système 32 utilise les relations suivantes: (1) QMCS = QMCE + yEG1 x QEGT - YEV1 x QEV (2) BMCS = (BMCE x QMCE + BEGT x YEG1 x QEGT) / QMCS (3) PMCS (PMCE x BMCE x QMCE + PEGT x BEGT x YEG1 x QEGT) / (QMCS x BMCS) (4) SD = K x Sp x M x C x ( 1-- l) x (5) RDTS SD + M MES -OM et (6) OMS =(QMCS x RDTS UQMCE x RDTE) QMCE, BMCE, PMCE, TMCE, RDTE et OME étant, respectivement le débit, le brix, la pureté, la température, la teneur en cristaux et la dimension moyenne des cristaux de la masse cuite entrant dans un compartiment, QMCS, BMCS, PMCS, TMCS, RDTS et OMS les grandeurs correspondantes pour la masse cuite sortant d'un compartiment, PEGT et TEGT la pureté et la température de l'égout, SD la masse de sucré déposée par cristallisation dans le compartiment, K le coefficient d'échange massique, Sp la surface spécifique des cristaux, M le débit massique des cristaux (M = QMCE x RDTE), C la concentration en sucre de l'eau-mère dans le compartiment, r la sursaturation de l'eau- mère dans le compartiment, et le temps de séjour dans le compartiment. Pour les bilans thermiques d'eau évaporée, le système 32 utilise la formule suivante: - (7) QEV1 (QMCE x CpMCE x TMCE + YEG1 x QEGT x CpEGT x TEGT - QMCS x CpMCS x TMCS - Pertes thermiques + h x &S x tu) / H -9- CpMCE, CpMCS et CpEGT étant, respectivement, la chaleur spécifique de la masse-cuite entrant dans le compartiment, de la masse cuite sortant du compartiment et de l'égout, h étant le coefficient d'échange thermique, &S la surface d'échange thermique du compartiment, 6tu l'écart de tempé- rature utile, et H l'enthalpie de la vapeur de chauffage. Pour le premier compartiment les valeurs de QMCE, BMCE, PMCE, RDTE, TMCE et CpMCE sont celles des grandeurs corres- pondantes du magma. Pour chaque compartiment, le système 32 effectue d'abord un bilan de cristallisation pour calculer les valeurs des grandeurs caractéristiques de la masse cuite sortant du compartiment, notamment de sa teneur en cristaux RDTS, en fonction des valeurs des grandeurs caractéristiques de la masse cuite et de l'égout entrant dans le compartiment et des paramètres caractérisant la marche de l'appareil. Pour ce calcul, le système suppose connus la pression de la vapeur de chauffage PVAP, le débit d'eau évaporée dans le comparti- ment QEV1 et la teneur en cristaux de la masse cuite sortant du compartiment RDTS. Au début de chaque pas de calcul d'indice i, le système attribue à ces grandeurs les valeurs PVAPi_1i QEVi_1 et RDTS i-1 déterminée au pas de calcul précédent d'indice i-1. Si la valeur calculée RDTSi est égale, à la précision recherchée près, à RDTSi.i, le système retient cette valeur et passe à l'étape suivante du calcul. Si non, un processus itératif permet de calculer une valeur RDTSj, j étant l'indice de cette boucle d'itération, à partir d'une valeur RDTAj attribuée par le système et fonction des valeurs précédentes RDTSj-1 et RDTAj-1 suivant une loi prédéterminée. Le calcul sera recommencé jusqu'à ce que la valeur calculée RDTSj soit égale, avec l'approximation désirée, à la valeur RDTAj attribuée par le système; la valeur retenue de RDTSi sera alors égale à RDTSj. Le système effectue alors le bilan thermique d'eau évaporée du compartiment. Les calculs de bilan de cristalli- sation et de bilan thermique d'eau évaporée sont recommencés jusqu'à ce que la valeur calculée QEVk soit égale, avec l'approximation désirée, à la valeur QEVAk attribuée par le - 10 - système. Le système effectue ainsi un bilan de cristallisation et un bilan thermique d'eau évaporée pour chacun des compar- timents, du premier au dernier, en utilisant comme valeurs des grandeurs caractéristiques de la masse cuite entrant dans un compartiment les valeurs calculées pour la masse cuite sortant du compartiment précédent. Lorsque les bilans ont été bouclés pour tous les compartiments, on compare la valeur calculée de l'une des grandeurs caractéristiques de la masse cuite sortant du dernier compartiment, c'est-à-dire, la masse cuite extraite de l'appareil de cristallisation,à la valeur imposée pour cette grandeur. Dans l'exemple décrit on a choisi le brix de la masse cuite. La valeur calculée du brix BMCC est en fait comparée à la valeur imposée BMCI augmentée d'un terme correctif BMC, qui sera défini ci-après et qui peut être positif, négatif ou nul. Si ces deux valeurs sont égales, à la précision recherchée près, on arrête le calcul et on retient pour valeur optimale de la pression de vapeur de chauffage la valeur PVAPi1 attribuée par le système au début du calcul. Si-non, on remplace cette valeur par une nouvelle valeur calculée suivant une loi déterminée en fonction de l'écart entre BMCI + 4BMC et BMCC, et on recommence tous les calculs en utilisant cette nouvelle valeur. Si le brix BMCC calcul6 à partir de cette nouvelle valeur ne satisfait pas l'égalité BMCC = BMCI + ABMC on recommence l'ensemble des calculs après avoir à nouveau modifié la valeur de PVAP. Lorsque cette égalité est véri- fiée avec l'approximation désirée, les calculs sont arrêtés et la dernière valeur de PVAP retenue comme valeur optimale de la pression de la vapeur de chauffage. La valeur calculée du brix de la masse cuite BMCC est alors comparée à la valeur BMCM mesurée au moyen d'un instrument 54. Si ces deux valeurs diffèrent, le système calcule le terme correctif &BMC qui est fonction de la différence BMCI - BMCM. Ce terme correctif est mis en mé- moire pour être utilisé au pas de calcul suivant comme décrit ci-dessus. 2562 908 - l1 - Les débits d'alimentation en égout des compartiments ou groupes de compartiments sont maintenus égaux aux va- leurs optimales détermin6es comme expliqué plus haut par les régulateurs 44 commandant des vannes placées sur les tuyauteries d'alimentation, ces valeurs optimales étant utilisées comme valeurs de consigne des régulateurs. En variante, on pourrait utiliser des vannes tout ou rien commandées par le système de conduite automatisée, le con- trole des débits s'effectuant alors en boucle ouverte à par- tir de la mesure du débit en position vanne ouverte et en jouant sur les temps d'ouverture et de fermeture des van- nes. La pression de la vapeur de chauffage est maintenue égale à la valeur optimale calculée comme décrit ci-dessus et utilisée comme valeur de consigne du régulateur 40. La pression dans la calandre de l'appareil de cristallisation est maintenue par le régulateur 42 à une valeur constante fixée par l'exploitant. Les débits de vapeur d'agitation et de magma sont maintenus à des valeurs de consigne égales aux valeurs optimales de ces grandeurs par les r6gulateurs 46 et 48 commandant, respectivement, une vanne placée sur la tuyau- terie 28 et le moteur d'entraînement de la pompe 22. Le niveau de masse cuite dans l'appareil est maintenu à une valeur de consigne par le régulateur 50. On pourra maintenir ce niveau constant ou, en variante, admettre des variations de niveau, sur une plage étroite, par exemple pour maintenir le malaxeur aval plein en permanence. De même, en cas de variation importante du débit aval, imposée par la cadence d'essorage par exemple, le débit de la pompe dtex- traction de masse cuite sera modifié pour maintenir le niveau du malaxeur aval et les débits d'égout et de magma seront modifiés en conséquence, dans des limites imposées en valeur haute par la pression maximum de vapeur de chauffage et en valeur basse par le niveau haut des capacités tam- pons placées en a.iont du cristalliseur. En cas de dé- pa$sseent de ces limites, une limitation des d6bits des aeLiers dmont ou aval, suivant le cas, sera nécessaire. - 12 -2562908 Les fichiers utilisés pour déterminer les valeurs optimales des débits d'alimentation des compartiments et du débit de la vapeur d'agitation, pourraient être rempla- cés par des relations mathématiques traduisant les lois YEG1 = f(QEGT, BEGT, R, PCAL) et QINJ = h(QEGT, BEGT, R, PCAL) respectivement. Les grandeurs et paramètres non imposés par l'exploi- tant qui interviennent dans le procédé de conduite automa- tisée objet de la présente invention peuvent être soit mesurés en continu, soit déterminés et introduits dans le système périodiquement, par exemple au moyen d'un cla- vier. Comme cela résulte de la définition générale de l'invention donnée dans le préambule, de nombreuses autres formes de mise en oeuvre sont possibles en choisissant différemment le paramètre dont la valeur optimale est déterminée au moyen d'un fichier ou d'unerelation mathé- matique et les deux paramètres dont les valeurs sont impo- sées. Ce choix dépendra du but visé: production maximale de sucre, réduction de la consommation d'énergie, etc... - 13 - R E V E ND I C AT I 0 N S 1. Procédé de conduite automatisée d'un appareil de cris- tallisation à marche continue pour la production de sucre, constitué par une cuve fermée disposée horizontalement et divisée par des cloisons verticales en plusieurs compar- timents munis d'un faisceau de chauffage, une solution sucrée étant introduite dans chacun des compartiments, un magma formé d'un mélange de petits cristaux et de li- queur mère étant introduit dans le premier compartiment et la masse cuite produite étant extraite du dernier compar- timent, caractérisée en ce qu'on effectue périodiquement les opérations suivantes: on détermine, pour chaque compar- timent, la valeur optimale d'un paramètre d'un premier groupe de cinq paramètres constitué par le débit d'alimen- tation en solution sucrée des compartiments rapporté au débit total de solution sucrée à traiter (YEG1), le débit d'eau évaporée dans les compartiments rapporté au débit total d'eau évaporée dans l'appareil de cristallisation (YEV1), le débit de masse cuite sortant des compartiments rapporté au débit total de masse cuite produite (yMC1), le brix de la masse cuite dans les compartiments (BMC1) et la teneur en cristaux de sucre de la masse cuite dans les compartiments (RDT1), en utilisant une loi préétablie donnant la valeur dudit paramètre dans les compartiments en fonction du brix de la solution sucrée à traiter (BEGT), de la pression de vapeur dans la calandre de l'appareil de cristallisation (PCAL)etdu rapport R = (QMAG x RDTMAG) / (QEGT x BEGT) avec QMAG: débit de magma RDT MAG: teneur n cristaux du magma QEGT: débit total de solution sucrée à traiter BEGT: brix de la solution sucrée à traiter, on calcule les valeurs optimales de quatre paramètres, d'un second groupe de six paramètres constitué par le débit total de solution sucrée à traiter (QEGT), le débit total d'eau évaporée (QEV), le débit total de masse cuite produite (QMC), le brix de la masse cuite produite (BMC), la pression de la _- 14 - 2562908 vapeur de chauffage (PVAP) et la teneur en cristaux de sucre de la masse cuite produite (RDT), ce second groupe comprenant deux paramètres de conduite (QEGT et PVAP) et quatre para- mètres d'état résultants (QEV, QMC, BMC et RDT), en utili- sant pour ce calcul la valeur optimale déterminée pour le paramètre du premier groupe) les valeurs imposées des deux autres paramètres du second groupe dont l'un est un para- mètre de conduite et l'autre un paramètre d'état résultant et éventuellement un terme correctif dont la valeur est fonction de l'écart entre la valeur imposée dudit paramètre d'état résultant et la valeurmesurée de celui-ci et on règle l'autre paramètre de conduite de ce second groupe pour le maintenir à une valeur de consigne égale à la valeur optimale calculée. 2. Procédé selon la revendication 1, caractérisé en ce que les lois donnant les variations des paramètres du premier groupe sont fonction du débit total de solution à traiter (QEGT) ou du débit total d'eau évaporée (QEV) ou du débit total de masse cuite produite (QMC), suivant le choix des paramètres du second groupe à optimiser. 3. Procédé selon la revendication 1 ou 2, caractérisé en ce que la valeur optimale pour chaque compartiment du para- mètre du premier groupe est déterminée au moyen d'une rela- tion mathématique ou d'un fichier de courbes ou de tableaux donnant les valeurs de ce paramètre pour les différents compartiments pour différentes valeurs du brix de la solution sucrée à traiter (BEGT), dudit rapport (R), de la pression de vapeur dans la calandre de l'appareil de cristallisation et éventuellement du débit total de la solution sucrée à traiter (QEGT) ou du débit total de masse cuite produite (QMC) ou du débit total d'eau évaporée (QEV). 4. Procédé selon la revendication 1 ou 2, caractérisé en ce qu'on utilise, à la place du paramètre d'état résultant dont on compare les valeurs imposées et mesurées, une gran- deur représentative dudit paramètre. - 15 - 2562908 - 15 - 5. Procédé selon la revendication 1, 2, 3 ou 4, caractérisé en ce que, pour calculer les valeurs optimales des quatre paramètres du second groupe, on effectue les bilans de cristallisation et thermique d'eau évaporée pour chaque compartiment. 6. Procédé selon l'une quelconque des revendications pré- cedentes, caractérisé en ce qu'on détermine périodiquement les valeurs optimales des-débits d'alimentation en solution sucrée de chaque compartiment en utilisant une loi de répartition du débit dans les compartiments qui est fonc- tion du débit total de solution à traiter (QEGT), du brix de la solution (BEGT), dudit rapport (R) et de la pression de vapeur dans la calandre de l'appareil (PCAL), on calcule la valeur optimale de la pression de vapeur de chauffage (PVAP) en utilisant comme paramètres imposés le débit total de solution sucrée (QEGT), le brix de la masse cuite produite (BMC) et éventuellement un terme correctif ( BMC) dont la valeur est fonction de l'écart entre la valeur impo- sée du brix de la masse cuite (BMCI) et une valeur mesurée de ce brix (BMCM), et on règle la pression de vapeur de chauffage pour le maintenir à une valeur de consigne égale à la valeur optimale calculée. 7. Procédé selon la revendication 6, caractérisé en ce que la pression d vapeur dans la calandre de l'appareil de cris- tallisation (PCAL) est maintenue constante. 8. Procédé selon la revendication 6 ou 7, caractérisé en ce que le débit de magma introduit dans le premier compartiment (QMAG) est réglé de façon à maintenir ledit rapport (R) constant. 9. Procédé selon la revendication 6, 7 ou 8, caractérisé en ce qu'on injecte de la vapeur dans chaque compartiment, à sa partie infèrieure, on détermine la valeur optimale d: débit de vapeur injectée (QINJ) au moyen d'une ou plu- - 16 - sieurs relations mathématiques ou d'un jeu de courbes ou de tableaux donnant les variations de ce débit en fonction du débit total de solution à traiter (QEGT) pour diffé- rentes valeurs du brix de la solution (BEGT), dudit rap- port (R) et de la pression de vapeur dans la calandre de l'appareil de cristallisation (PCAL) et on règle le débit de vapeur injectée pour le maintenir à une valeur de con- signe égale A cette valeur optimale.