FR 2602243 A2 19880205 FR 8611191 A 19860801 La présente invention a pour objet un procédé de conduite automatisée d'un appareil de cristallisation a marche continue pour la production de sucre, constitue par une cuve fermée, divisée par des cloisons en plusieurs compartiments 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 liqueur-mère étant introduit dans le premier compartiment et la masse cuite produite étant extraite du dernier compartiment. Selon 1 'une des formes de mise en oeuvre du procédé objet du brevet principal, on détermine périodiquement la valeur optimale d'un premier parametre, par exemple le débit d'alimentation en solution sucrée de chaque compartiment, en utilisant un loi préétablie donnant la valeur dudit paramètre dans les compartiments en fonction du brix de la solution à traiter, de la pression de vapeur dans la calandre de l'appareil et du rapport R = (QMAG x RDTMS) / (BEAT x BEGT) avec QRG . débit de magma EDTMAG : teneur en cristaux du magma QEGT : débit total de solution sucrée à traiter BEGT : brix de la solution sucrée à traiter, on calcule la valeur optimale de la pression de vapeur de chauffage ou du débit total de solution à traiter en utilisant pour ce calcul les valeurs imposées du dépit total de solution sucrée ou de la pression de valeur de chauffage, respectivement, et du brix de la masse cuite produite et éventuellement un terme correctif dont la valeur est fonction de 1 'écart entre la valeur imposée du brix de la masse cuite et la valeur mesurée de ce brix, et on règle la pression de vapeur de chauffage ou le débit total de solution, respectivement, pour le maintenir à une valeur de consigne égale à la valeur optimale calculée. Le but de la présente invention est d'apporter a ce procédé des perfectionnements permettant une autoadaptation du modèle mathématique utilise pour les calculs aux conditions réelles de fonctionnement de l'appareil. Suivant une première caractéristique de l'invention, les valeurs d'un ou plusieurs paramètres utilises dans les calculs de bilan thermique de chaque compartiment sont déterminées, à chaque pas de calcul, en fonction de la valeur calculée au pas de calcul précédent et de la valeur mesurée du débit d'eau condensée dans le faisceau de chauffage. On effectue ainsi, à chaque pas de calcul, une correction des paramètres utilises pourle calcul, ce qui permet notamment de prendre en compte 1 'encrassement progressif des éléments du faisceau de chauffage entre deux nettoyages. Cette correction peut porter sur les coefficients -d'échange thermique etiou les surfaces d'échange thermique utilisés pour le calcul des bilans thermiques.Les coefficients d'échange thermique peuvent, par exemple, être déterminés au moyen de la formule h = ho (a + b .QCOND + D.VÀP) 2QNOM dans laquelle QCOND est le débit mesuré d'eau condensée dans le faisceau de chauffage, QVAP est le débit d'eau condensée calculé au pas de calcul précédent, QNOM est le débit total d'eau condensée correspondant à la valeur ho du coefficient d'échange et a et b sont des constantes propres à chaque compartiment ou groupe de compartiments. Pour le calcul des surfaces d'échange effectives ou utiles on pourra utiliser la formule surface échange effective = surface d'échange théorique . FXj avec FXj = 1-1 - - XENCRA . (1 - FXDj) 1 - FXDn n étant le nombre de compartiments, j l'indice du compartiment considéré, FXDj une constante caractéristique du compartiment et XENCRN une grandeur fonction de 1 'écart entre la valeur mesurée QCOND et la valeur calculée BVAP du débit total d'eau condensée. XENCP peut, par exemple, être él borée par un régulateur de type PID recevant, en entrée, les valeurs mesurée et calculée du débit d'eau condensée. Selon une autre caractéristique de l'invention, on passe automatiquement d'un mode de conduite avec détermination de la valeur optimale de la pression de la vapeur de chauffage, le débit total de solution à traiter étant maintenu à une valeur prédéterminée, à un mode de conduite avec détermination de la valeur optimale du débit total de solution à traiter, la pression de la vapeur de chauffage étant imposée, lorsque la valeur optimale de la pression de vapeur de chauffage est supérieure à la pression de vapeur disponible, et on revient au premier mode de conduite lorsque la pression de vapeur disponible devient supérieure à la valeur imposée de la -pression de la vapeur de chauffage. La description qui suit d'un mode de mise en oeuvre de l'invention, donné à titre d'exemple non limitatif, se réfère aux dessins qui l'accompagnent et sur lesquels La figure 1 est le schéma d'un système de conduite automatique d'un appareil de critallisation pour la mise en oeuvre du procédé de l'invention, et La figure 2 est l'organigramme du procédé. Sur la figure 1, le numéro de référence 1 désigne un appareil de cristallisation par évaporation à marche continue constitué par une cuve fermée disposée horizontalement et divisée par des cloisons verticales en plusieurs compartiments 1, 2, 3, 4, ..., n-2, n-l, 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 1 commun a 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 i8 auquel les tuyauteries 16 sont reliées par une conduite 20. On peut prévoir une tuyauterie d'alimentation par compartiment ou une tuyauterie commune a 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 dispositif 52 du type faisant l'objet du brevet 84.t)0101. Un magma constitué par un mélange de petits cristaux constituant des germes de cristallisation et de liqueurmè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 cc-3artiment 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 S cette vapeur amenée par une conduite 28 est distribuée a des injecteurs 30 placés dans le bas des compartiments. Le numéro de référence 32 désigne un système de conduite automatisée se composant d'un ordinateur de supervision couplé a un système de conduite de processus comprenant les régulateurs numériques et l'instrumentation ou à des chaines de régulation numériques déportées, ou à des chaines de régulation analogiques classiques. En variante l'ordinateur peut assurer les fonctions de régulation et communiquer avec l'instrumentation par des convertisseurs digital-analogique et analogique-dig,ital. Le système reçoit des informations numériques et analogiques, les traite et émet les consignes qu'il a optimisées vers les actionneurs : vannes, moteurs5 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 Z2 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, respectivement, 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 1 '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 4, un régulateur 42 de la pression de vapeur dans la calandre de l'apPareíl. de 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 a l'invention du brevet principal, 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 chauffage 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, par exemple, de 15 secondes. Les valeurs optimales des débits d'alimentation EGj des compartiments, du débit de magma QMG et du débit de vapeur d'agitation QINJ sont déterminées comme décrit dans le brevet principal. La valeur optimale de la pression de la vapeur de chauffage PVAP est calculée par le système 32 suivant le prt-édé présenté sous forme d'organigramme sur la figure 2. Le système 32 effectue pour chaque compartiment un bilan de 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 1 'égout a 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'appareil, etc...) comme décrit au brevet principal Pour effectuer les bilans thermiques, le système 32 utilise la formule suivante D-.EVj = (D-.MCE X EhMCE + YEGj x FIENT x EhEGT - QMCS x EhMCS - Pertes thermiques + h x #S x #tu) / H dans laquelle QEVj est la quantité d'eau évaporée dans le compartiment j, QMCE et D-.MCS étant, respectivement, les débits de la masse cuite entrant dans le compartiment, et sortant du compartiment, YEGj étant le débit d'alimentation en égout du compartiment rapporté au débit total d'égout, EhMCE, EhMCS et EhEGT étant, respectivement, l'enthalpie 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 effective du compartiment, d tu l'écart de température utile, et H l'enthalpie spécifique de la vapeur de chauffage. Pour le premier compartiment les valeurs de QMCE, BMCE, PMCE, RDTE, TMCE et EhMCE sont celles des grandeurs correspondantes du magma. Conformément a la présente invention, les coefficients d'échange thermique sont donnés par la formule h = ho . (a + b QCOND + QVAP) 2 QNOM et leur valeur est déterminée à chaque pas de calcul en fonction du débit d'eau condensée SECOND mesuré au moyen du débitmètre 56 et dudébit d'eau condensée QVAP calcule par le système 32, a et b étant des constantes propres à chaque compartiment ou groupe de compartiments et ho le coefficient d'échange thermique théorique correspondant å des conditions de fonctionnement prises pour référence (débit nominal, appareil propre...) et pour lesquelles le débit d'eau condensée est égal à QNOM. Les valeurs de a, b, ho et ONOM sont stockées dans un fichier du système 32. La surface d'échange thermique effective utilisée pour le calcul est donnée par la formule #S = #S théorique . FXj avec FXj = 1~1 - XENCRA . (1 - FXDj) 1- FXDn Les valeurs théoriques des surfaces d'échange thermique ( d S théorique) des différents compartiments sont calculées initialement à partir des données géométriques (diamètre et épaisseur des tubes, écartement dC tubes du faisceau...) et stockées dans un fichier du système 32. Les constantes FXDj sont également stockées dans un fichier. Le coefficient XENCR permet de tenir compte de l'encrassement des tubes du faisceau. Sa 'valeur est celle du signal de sortie d'un régulateur numérique de type PID, incorporé au système 32 et à l'entrée duquel sont appliqués des signaux représentant les valeurs SECOND (débit d'eau condensée mesurez par le débitmètre 56) et QVNF' (débit d'eau condensée calculée). La valeur de XENCRR est donc fonction de l'écart entre les valeurs de OCOND et OVAP. Le système effectue un bilan de cristallisation et un bilan thermique d'eau évaporée pour chacun des compartiments, 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 du brix BMCC à la valeur imposée BMCI. Si ces deux valeurs sont égales, a 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 PVAPi-1 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 et BMCC, et on recommence tous les calculs en utilisant cette nouvelle valeur.Si le brix BMCC calcule a partir de cette nouvelle valeur ne satisfait pas l'égalité BMCC = BMCI on recommence 1 'ensemble des calculs après avoir à nouveau modifié la valeur de PV4P. Lorsque cette égalité est vérifié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 et utilisée comme valeur de consigne du régulateur 40. Le débit d'eau condensée QVRP calculé par le système Z2 å chaque pas de calcul est utilisÉ, avec le débit d'eau condensée BCOND mesuré par le débitmètre 56, pour déterminer les nouvelles valeurs des coefficients d'échange thermique h et du coefficient XENCRçi, et cellesci sont utilisées au pas de calcul suivant. Le fonctionnement décrit ci-dessus, qui consiste a imposer un débit total d'égout et à régler la pression de la vapeur de chauffage pour la maintenir à la valeur optimale calculée par le système 32 est le fonctionnement normal. Cependant, lorsque la pression sous laquelle la vapeur de chauffage est disponible n'est pas suffisante, la vanne 39 commandée par le régulateur 40 s ouvre en grand sans que l'on puisse amener la pression à la valeur optimale calculée.Conforménent a 1 'invention, le système maintient alors la pression de la vapeur de chauffage constante et inférieure ou au plus égale à la pression maximale sous laquelle la vapeur peut être fournie et calcule la valeur optimale du débit total d'égout alXnentant l'appareil 10 puis, au moyen d'une loi de répartition préétablie, les débits d'alimentation en égout des différents compartiments, les valeurs calculées étant imposées comme consignes a régulateurs de-débit partiels 44, pour maintenir le brix de la masse cuite produite à la valeur imposée. Le passage du premier mode de conduite au second mode s'effectue automatiquement lorsque le degré d'ouverture de la vanne 39 depasse une valeur prédéterminée, par exemple 80% Le retour du premier mooe de conduite s'effectue lorsque le degré d'ouverture de la vanne 39 redevient inférieur a une seconde valeur prédéterminée, inférieure à la première, par exemple 70% REVENDICATIONS 1. Procédé de conduite automatisée d'un appareil de cristallisation à marche continue pour la production de sucre selon les revendications 1 et 5 du brevet principal, caractérisé en ce que les valeurs d'un ou plusieurs paramètres (h, t S) utilisés dans les calculs de bilan thermique des différents compartiments sont déterminées à chaque pas de calcul en fonction de la valeur calculée au pas de calcul précédent (QVAP) et de la valeur mesurée (OCOND) du débit d'eau condensée dans le faisceau de chauffage. 2 Procédé selon la revendication 1, caractérisé en ce que ledit paramètre est le coefficient d'échange thermique a travers les éléments du faisceau de chauffage (h). 3. Procédé selon la revendication 1 ou 2, caractérisé en ce que ledit paramètre est la surface d'échange thermique des éléments du faisceau de chauffage (1 Si. 4, Procédé selon la revendication 2, caractérisé en ce que le coefficient d'échange thermique (h) est déterminé au moyen de la formule : h = ho (a + b . OCOND + VAP) 2 BNOFî dans laquelle ho, a et b sont des constantes propres a chaque compartiment et QNOM est une constante représentant le débit d'eau condensée produit quand h = ho. 5. Procédé selon la revendication 3, caractérisé en ce que la surface d'échange thermique (#S) est déterminée au moyen des formules #S = #S théorique x FXj et FXj = 1~1 - XENCA ( 1 . FXDj) 1 - FXDn dans lesquelles #S théorique et FXDj sont des constante. propres à chaque compartiment et XENCRm est un coefficient correcteur qui est fonction de 1 'écart entre la valeur mesurée (QCOND) et la valeur calculée (QVAP) du débit d'eau condensée. 6. Procédé selon la revendication 5, caractérisé en ce que la valeur dudit coefficient correcteur (XENCP) est celle du signal de sortie d'un régulateur de type PID à l'entrée duquel sont appliqués des signaux proportionnels aux valeurs mesurées et calculées (QCOND, QVNP) du débit d'eau condensée. 7. Procédé selon l'une quelconque des revendications précédentes, caractérisé en ce qu'on calcule la valeur optimale de la pression de la vapeur de chauffage (PVAP) et on maintient la pression de la vapeur de chauffage a la valeur optimale calculée, le débit total de solution å traiter (EGT) étant maintenu a une valeur prédéterminée, tant que ladite valeur optimale est inférieure à la pression de vapeur disponible et on calcule la valeur optimale du débit total de solution à traiter et on maintient ce débit à la valeur optimale calculée, la pression de vapeur de chauffage étant maintenue à une valeur imposée, tant que la pression de vapeur disponible reste inférieure à la valeur imposée.