NO315563B1 - Fremgangsmåte for styring av kjemiske synteseprosesser - Google Patents
Fremgangsmåte for styring av kjemiske synteseprosesser Download PDFInfo
- Publication number
- NO315563B1 NO315563B1 NO19975780A NO975780A NO315563B1 NO 315563 B1 NO315563 B1 NO 315563B1 NO 19975780 A NO19975780 A NO 19975780A NO 975780 A NO975780 A NO 975780A NO 315563 B1 NO315563 B1 NO 315563B1
- Authority
- NO
- Norway
- Prior art keywords
- reactor
- measurements
- algorithm
- properties
- predictions
- Prior art date
Links
- 238000000034 method Methods 0.000 title claims abstract description 201
- 238000003786 synthesis reaction Methods 0.000 title claims abstract description 52
- 230000008569 process Effects 0.000 claims abstract description 133
- 230000015572 biosynthetic process Effects 0.000 claims abstract description 49
- 239000000126 substance Substances 0.000 claims abstract description 16
- 238000005259 measurement Methods 0.000 claims description 108
- 238000004422 calculation algorithm Methods 0.000 claims description 104
- 238000004364 calculation method Methods 0.000 claims description 66
- 230000006978 adaptation Effects 0.000 claims description 49
- 239000003153 chemical reaction reagent Substances 0.000 claims description 37
- 229920000642 polymer Polymers 0.000 claims description 37
- 238000004519 manufacturing process Methods 0.000 claims description 36
- 238000006116 polymerization reaction Methods 0.000 claims description 33
- 239000004698 Polyethylene Substances 0.000 claims description 31
- 229920000573 polyethylene Polymers 0.000 claims description 31
- 239000000470 constituent Substances 0.000 claims description 30
- VGGSQFUCUMXWEO-UHFFFAOYSA-N Ethene Chemical compound C=C VGGSQFUCUMXWEO-UHFFFAOYSA-N 0.000 claims description 27
- 239000005977 Ethylene Substances 0.000 claims description 26
- 239000003054 catalyst Substances 0.000 claims description 26
- 239000001257 hydrogen Substances 0.000 claims description 20
- 229910052739 hydrogen Inorganic materials 0.000 claims description 20
- -1 polyethylene Polymers 0.000 claims description 18
- UFHFLCQGNIYNRP-UHFFFAOYSA-N Hydrogen Chemical compound [H][H] UFHFLCQGNIYNRP-UHFFFAOYSA-N 0.000 claims description 16
- 238000006243 chemical reaction Methods 0.000 claims description 15
- 238000004497 NIR spectroscopy Methods 0.000 claims description 8
- 238000013178 mathematical model Methods 0.000 claims description 8
- 230000000694 effects Effects 0.000 claims description 7
- 150000002431 hydrogen Chemical class 0.000 claims description 7
- 239000000463 material Substances 0.000 claims description 7
- 239000000155 melt Substances 0.000 claims description 7
- 238000002156 mixing Methods 0.000 claims description 7
- QQONPFPTGQHPMA-UHFFFAOYSA-N propylene Natural products CC=C QQONPFPTGQHPMA-UHFFFAOYSA-N 0.000 claims description 7
- 125000004805 propylene group Chemical group [H]C([H])([H])C([H])([*:1])C([H])([H])[*:2] 0.000 claims description 7
- 238000005033 Fourier transform infrared spectroscopy Methods 0.000 claims description 4
- 238000005481 NMR spectroscopy Methods 0.000 claims description 4
- 239000004743 Polypropylene Substances 0.000 claims description 4
- 239000000654 additive Substances 0.000 claims description 4
- 230000000996 additive effect Effects 0.000 claims description 4
- 229920001155 polypropylene Polymers 0.000 claims description 4
- 239000000178 monomer Substances 0.000 claims description 3
- 238000012546 transfer Methods 0.000 claims description 3
- 238000009529 body temperature measurement Methods 0.000 claims description 2
- VXNZUUAINFGPBY-UHFFFAOYSA-N 1-Butene Chemical compound CCC=C VXNZUUAINFGPBY-UHFFFAOYSA-N 0.000 description 13
- IAQRGUVFOMOMEM-UHFFFAOYSA-N butene Natural products CC=CC IAQRGUVFOMOMEM-UHFFFAOYSA-N 0.000 description 12
- 238000012937 correction Methods 0.000 description 12
- 239000002904 solvent Substances 0.000 description 12
- 238000010586 diagram Methods 0.000 description 10
- 230000001276 controlling effect Effects 0.000 description 9
- 239000000203 mixture Substances 0.000 description 9
- 239000000243 solution Substances 0.000 description 9
- 239000007789 gas Substances 0.000 description 8
- 239000012071 phase Substances 0.000 description 7
- 230000003068 static effect Effects 0.000 description 7
- 230000008901 benefit Effects 0.000 description 6
- VLKZOEOYAKHREP-UHFFFAOYSA-N n-Hexane Chemical compound CCCCCC VLKZOEOYAKHREP-UHFFFAOYSA-N 0.000 description 6
- 230000008859 change Effects 0.000 description 5
- 239000011347 resin Substances 0.000 description 5
- 229920005989 resin Polymers 0.000 description 5
- 230000001934 delay Effects 0.000 description 4
- 238000007726 management method Methods 0.000 description 4
- 230000003044 adaptive effect Effects 0.000 description 3
- 230000003197 catalytic effect Effects 0.000 description 3
- 239000012530 fluid Substances 0.000 description 3
- 238000010438 heat treatment Methods 0.000 description 3
- 239000007791 liquid phase Substances 0.000 description 3
- 230000009466 transformation Effects 0.000 description 3
- ATUOYWHBWRKTHZ-UHFFFAOYSA-N Propane Chemical compound CCC ATUOYWHBWRKTHZ-UHFFFAOYSA-N 0.000 description 2
- 238000013459 approach Methods 0.000 description 2
- 238000001816 cooling Methods 0.000 description 2
- 239000012809 cooling fluid Substances 0.000 description 2
- 238000000786 liquid-assisted grinding Methods 0.000 description 2
- 238000002844 melting Methods 0.000 description 2
- 230000008018 melting Effects 0.000 description 2
- 239000002574 poison Substances 0.000 description 2
- 231100000614 poison Toxicity 0.000 description 2
- 230000000379 polymerizing effect Effects 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 101001019013 Homo sapiens Mitotic interactor and substrate of PLK1 Proteins 0.000 description 1
- 102100033607 Mitotic interactor and substrate of PLK1 Human genes 0.000 description 1
- 150000001336 alkenes Chemical class 0.000 description 1
- 230000003321 amplification Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000000454 anti-cipatory effect Effects 0.000 description 1
- 239000003426 co-catalyst Substances 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000010924 continuous production Methods 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 239000006185 dispersion Substances 0.000 description 1
- 238000004821 distillation Methods 0.000 description 1
- 238000001035 drying Methods 0.000 description 1
- 238000004836 empirical method Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000001704 evaporation Methods 0.000 description 1
- 230000008020 evaporation Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000012685 gas phase polymerization Methods 0.000 description 1
- 238000005469 granulation Methods 0.000 description 1
- 230000003179 granulation Effects 0.000 description 1
- 239000012535 impurity Substances 0.000 description 1
- 238000002347 injection Methods 0.000 description 1
- 239000007924 injection Substances 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 230000001788 irregular Effects 0.000 description 1
- 239000007788 liquid Substances 0.000 description 1
- 230000007257 malfunction Effects 0.000 description 1
- 238000003199 nucleic acid amplification method Methods 0.000 description 1
- 230000000737 periodic effect Effects 0.000 description 1
- 231100000572 poisoning Toxicity 0.000 description 1
- 230000000607 poisoning effect Effects 0.000 description 1
- 238000005498 polishing Methods 0.000 description 1
- 238000004886 process control Methods 0.000 description 1
- 230000000135 prohibitive effect Effects 0.000 description 1
- 239000001294 propane Substances 0.000 description 1
- 239000002994 raw material Substances 0.000 description 1
- 230000009257 reactivity Effects 0.000 description 1
- 230000001105 regulatory effect Effects 0.000 description 1
- 238000010079 rubber tapping Methods 0.000 description 1
- 238000003860 storage Methods 0.000 description 1
- 239000000725 suspension Substances 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
- 239000004711 α-olefin Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B13/00—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion
- G05B13/02—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric
- G05B13/04—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators
- G05B13/048—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators using a predictor
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B01—PHYSICAL OR CHEMICAL PROCESSES OR APPARATUS IN GENERAL
- B01J—CHEMICAL OR PHYSICAL PROCESSES, e.g. CATALYSIS OR COLLOID CHEMISTRY; THEIR RELEVANT APPARATUS
- B01J19/00—Chemical, physical or physico-chemical processes in general; Their relevant apparatus
- B01J19/0006—Controlling or regulating processes
-
- C—CHEMISTRY; METALLURGY
- C08—ORGANIC MACROMOLECULAR COMPOUNDS; THEIR PREPARATION OR CHEMICAL WORKING-UP; COMPOSITIONS BASED THEREON
- C08F—MACROMOLECULAR COMPOUNDS OBTAINED BY REACTIONS ONLY INVOLVING CARBON-TO-CARBON UNSATURATED BONDS
- C08F10/00—Homopolymers and copolymers of unsaturated aliphatic hydrocarbons having only one carbon-to-carbon double bond
-
- C—CHEMISTRY; METALLURGY
- C08—ORGANIC MACROMOLECULAR COMPOUNDS; THEIR PREPARATION OR CHEMICAL WORKING-UP; COMPOSITIONS BASED THEREON
- C08F—MACROMOLECULAR COMPOUNDS OBTAINED BY REACTIONS ONLY INVOLVING CARBON-TO-CARBON UNSATURATED BONDS
- C08F2/00—Processes of polymerisation
-
- C—CHEMISTRY; METALLURGY
- C08—ORGANIC MACROMOLECULAR COMPOUNDS; THEIR PREPARATION OR CHEMICAL WORKING-UP; COMPOSITIONS BASED THEREON
- C08F—MACROMOLECULAR COMPOUNDS OBTAINED BY REACTIONS ONLY INVOLVING CARBON-TO-CARBON UNSATURATED BONDS
- C08F210/00—Copolymers of unsaturated aliphatic hydrocarbons having only one carbon-to-carbon double bond
- C08F210/16—Copolymers of ethene with alpha-alkenes, e.g. EP rubbers
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B11/00—Automatic controllers
- G05B11/01—Automatic controllers electric
- G05B11/32—Automatic controllers electric with inputs from more than one sensing element; with outputs to more than one correcting element
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B13/00—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion
- G05B13/02—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05D—SYSTEMS FOR CONTROLLING OR REGULATING NON-ELECTRIC VARIABLES
- G05D21/00—Control of chemical or physico-chemical variables, e.g. pH value
- G05D21/02—Control of chemical or physico-chemical variables, e.g. pH value characterised by the use of electric means
-
- C—CHEMISTRY; METALLURGY
- C08—ORGANIC MACROMOLECULAR COMPOUNDS; THEIR PREPARATION OR CHEMICAL WORKING-UP; COMPOSITIONS BASED THEREON
- C08F—MACROMOLECULAR COMPOUNDS OBTAINED BY REACTIONS ONLY INVOLVING CARBON-TO-CARBON UNSATURATED BONDS
- C08F2400/00—Characteristics for processes of polymerization
- C08F2400/02—Control or adjustment of polymerization parameters
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
- Y02P20/00—Technologies relating to chemical industry
- Y02P20/50—Improvements relating to the production of bulk chemicals
- Y02P20/582—Recycling of unreacted starting or intermediate materials
Landscapes
- Chemical & Material Sciences (AREA)
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Organic Chemistry (AREA)
- Chemical Kinetics & Catalysis (AREA)
- General Physics & Mathematics (AREA)
- Automation & Control Theory (AREA)
- Physics & Mathematics (AREA)
- Medicinal Chemistry (AREA)
- Polymers & Plastics (AREA)
- Artificial Intelligence (AREA)
- Medical Informatics (AREA)
- Software Systems (AREA)
- Evolutionary Computation (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Feedback Control In General (AREA)
- Physical Or Chemical Processes And Apparatus (AREA)
- Polymerisation Methods In General (AREA)
- Organic Low-Molecular-Weight Compounds And Preparation Thereof (AREA)
Description
Foreliggende oppfinnelse angår en fremgangsmåte for styring av prosesser for syntese av kjemiske produkter. Den angår også en styringsanordning for å gjøre bruk av denne fremgangsmåten, og en synteseprosess, spesielt syntese av polymer, styrt ved hjelp av denne fremgangsmåten.
I en prosess for syntese av kjemiske produkter som utføres på konvensjo-nell måte, benyttes det kontrollere av PID-type (proporsjonal-integral-differensial) for å styre individuelt et høyt eller lavt antall variabler (temperaturer, strømnings-hastigheter, trykkverdier osv.) som påvirker forløpet av syntesen. For hver temperatur, hver strømningshastighet eller hvert trykk som skal styres, måles med andre ord parameterens faktiske verdi kontinuerlig {eller intermitterende), og en PID-kontroller sammenligner denne faktiske verdien med et instillingspunkt, og virker på den variabel som skal styres for å redusere, dersom dette er korrekt, forskjel-len mellom innstillingspunktet og den målte verdien.
I betraktning av kompleksiteten av de fleste industrielle prosesser med kjemisk syntese, må instillingspunktene for de forskjellige kontrollerne i dag frem-deles justeres empirisk for til slutt å oppnå de ønskede egenskaper ved det syntetiserte produktet. Det benyttes oppskrifter for dette formål, som tilveiebringer kombinasjoner av parametre som er bestemt empirisk, for i en stabil tilstand å oppnå de ønskede egenskaper for det syntetiserte produktet.
Empiriske forhold mellom de kontrollerte variablene og egenskapene ved det syntetiserte produktet kan utledes fra disse oppskriftene ved hjelp av mer eller mindre sofistikerte statistikk-verktøy. Det er imidlertid åpenbart at disse empiriske forholdene neppe kan ta i betraktning de mange innbyrdes avhengigheter som eksisterer mellom de forskjellig variablene som styres separat, eller de ukjente perturbasjoner slik som råmaterialenes innhold av urenheter.
Det er også åpenbart at en tradisjonell styring av typen lukket sløyfe, som anvender målinger av vesentlige egenskaper ved det syntetiserte produktet som tilbakekoplings-korreksjoner, er vanskelig å anvende i de fleste synteseprosesser. Dette er fordi dødtidene som inngår enten i prosessen elle målingene eller analysene som benyttes som tilbakekoplings-korreksjoner, er alt for lange, og uavhen-gighetene mellom de forskjellige variablene som styrer prosessen, er altfor komplekse.
Internasjonal søknad WO 93/24533 beskriver en styringsmetode for en prosess for gassfase-polymerisering av en alfa-olefin i en horisontal reaktor, hvor styringsparametere gjør det mulig å agere i løpet av prosessen for å få polymerens smelte-indeks (MFR) lik det tilsvarende innstillingspunkt, idet metoden innbefatter de følgende trinn: - bestemmelse av forholdene mellom polymerens smelteindeks når den forlater reaktoren, og en første rekke parametere,
- styring av denne første rekken av parametere,
- beregning av polymerens MFR,
- tilpasning av minst en av parameterne for å justere den beregnede MFR til en forutbestemt verdi.
Det har lenge vært kjent at synteseprosesser, spesielt prosesser ved kontinuerlig polymer-syntese (polymeriserings-prosesser), utstyrt med kontrollere med empirisk justerte instillingspunkter, har vesentlige ulemper som kan oppsummeres som følger: - oppstartinger av syntesprosessen tar mye tid, og frembringer store mengder av
produkt utenfor spesifikasjonene;
kvalitetsforandringer foregår langsomt, hvilket også resulterer i produksjon av
store mengder av overgangs-produkter som er utenfor spesifikasjonene; - prosessens produksjonshastighet, dvs. massen av produkt(er) som syntetiseres pr. tidsenhet, er vanskelig å forandre uten å forverre egenskapene ved dette produktet eller disse produktene; - konsistens av det/de syntetiserte produktets/produktenes vesentlige egenskaper er ofte svak, selv i en stabil tilstand.
For å unngå empirisk justering av innstillingspunktene, er det foreslått i faglitteraturen å benytte styringsmetoder for synteseprosesser som gjør bruk av karakteristiske ligninger som benytter modelldannelse av synteseprosessen for å relatere noen egenskaper ved det/de syntetiserte produktet/produktene til reaktorens/reaktorenes arbeidsbetingelser under syntesen. For å begrense kompleksiteten av disse karakteristiske ligningene, er det imidlertid til nå regnet med at det i praksis er nødvendig enten å betrakte utelukkende det statiske tilfellet (stabil tilstand), eller å begrense seg til en høyt forenklet empirisk modelldannelse av prosessens dynamikk. Anvendelsen av en statisk modell er begrenset til styring av en produksjon med temmelig stabil tilstand.
I tilfellet med empirisk modelldannelse, er de karakteristiske ligningene gyldige bare for et smalt gyldighetsområde (nær det punkt hvor modelldannelsen er utført). I begge tilfeller styres oppstartings-fasene og overgangsfasene dårlig. Det er utvilsomt mulig å se for seg «dekning» av et bredere område av driftsbetingelser ved å foreta slik lokal modelldannelse på flere forskjellige punkter i området med driftsparametere, men en slik tilnærmelse blir uoverkommelig så snart det gjøres et forsøk på å styre et antall variabler ved å innvirke på mange parametere.
Det er følgelig ønskelig å ha tilgjengelig en enkel styringsmetode og en enkel styringsanordning som er bedre tilpasset til de spesifikke trekk ved dynamikken i syntesprosesser for kjemiske produkter.
For dette formål angår foreliggende oppfinnelse en fremgangsmåte for
styring av en synteseprosess for minst ett kjemisk produkt i et utstyr som omfatter minst en reaktor (R) som kan opptas i en perfekt blandet reaktor, hvor en eller flere manipulerte variabler (GC) gjør det mulig å innvirke på forløpet av prosessen for å gjøre en eller flere variabler som er relatert til produktets egenskaper og/eller til forløpet av prosessen, og som kalles styrte variabler (GR), lik de tilsvarende innstillingspunkter (Cgr) (eller i det minste så nær de
sistnevnte som mulig), hvor fremgangsmåten innbefatter de følgende trinn:
(a) innmating av innstillingspunkter som vedrører de styrte variablene (CGr); (b) beregning ved hjelp av en prediksjonsenhet (OP), av prediksjoner for de styrte variablene (Pgr), basert på målinger av de prosess-manipulerte variablene M<g>c); (c) anvendelse av en styringsenhet (OC) for å beregne innstillingspunktene for de prosess-manipulerte variablene (Cgc), basert på innstillingspunktene (CGr) og prediksjonene (PGr) for de styrte variablene; (d) overføring av instillingspunktene for de prosess-manipulerte variablene (Cgc) til aktuatorer eller til styringsenheter som styrer aktuatorene, for å innvirke på prosessens forløp;
hvor prediksjonsenheten (OP) er basert på en matematisk modell av prosessen, kalt en direkte modell (M), og er utformet på en slik måte at massen MXr av minst en bestanddel (X) i reaktoren (R) predikteres ved hjelp av ligningen:
hvor:
- FxRin er massestrømmen av bestanddelen X som går inn i reaktoren R;
- xxer oppholdstiden for X i reaktoren (tidskonstant), som er
hvor:
- MXRbetegner den sist beregnede verdi av massen av bestanddelen X som er tilstede i reaktoren R; - 2 Fxdis betegner summen av alle massestrømmer Fxdis som bestanddelen X forsvinner fra reaktoren R for, spesielt ved reaksjon og/eller ved
å forlate reaktoren;
- funksjonen y = LAG (u, x) er løsningen av differensialligningen
beregnet med den øyeblikkelige verdi av u og av x, og med den sist beregnede verdi av y.
Fordelen ved denne fremgangsmåten er at den ovenstående differensialligningen løses ved en enkel algebraisk beregning, f.eks. ved bruk av den følgen-de formel (hvor T betegner det tidsintervall.vanligvis kort i sammenligning med x, som atskiller de suksessive beregningene), eller ved å benytte en formel som er ekvivalent med den sistnevnte:
I det tilfelle hvor massene av flere bestanddeler blir evaluert slik som frem-stilt ovenfor, er fremgangsmåten ifølge oppfinnelsen spesielt fordelaktig, siden disse massene kan beregnes sekvensielt ved slik enkle algebraiske beregninger, ved hyppig om-beregning (generelt er T « x). På den annen side krever tradisjo-nelle metoder samtidig løsning av et system av differensialligninger, og dette krever generelt stor regnekraft og sofistikerte algoritmer for numeriske beregninger (integrasjon); som resultat av dette tar hver iterasjon med beregninger lang tid, og følgelig reagerer en styring av denne type dårlig på hurtige variasjoner.
Den styrte synteseprosessen kan benyttes for syntese av en monomerisk eller polymerisk forbindelse; svært gode resultater er oppnådd i tilfellet med styring av polymeriseringsprosesser. Prosessen strekker seg også til det tilfelle hvor et antall nyttige produkter syntetiseres samtidig i samme prosess. Prosessen kan være kontinuerlig eller ikke-kontinuerlig (sats); styrings-fremgangsmåten ifølge oppfinnelsen gir utmerkede resultater i tilfellet med kontinuerlige prosesser. Den styrte synteseprosessen kan eventuelt utgjøre bare en del av en større prosess, hvor de andre delene av prosessen styres på annerledes måte, eller ikke styres. For at styringsmetoden ifølge oppfinnelsen skal være anvendbar, er det nødven-dig at minst en reaktor kan bli opptatt i en perfekt blandet reaktor, dvs. en reaktor hvor de forskjellige variablene (temperatur, konsentrasjoner av tilstedeværende bestanddeler og lignende) er virtuelt identiske i hvert punkt. Andre mulige reaktorer kan være av stempel-typen (plugg-strøm); slike reaktorer modelldannes matematisk ved bruk av dødtider. Fremgangsmåten er også anvendelig for en prosess som finner sted i et antall reaktorer anordnet i serie og/eller i parallell, og i stand til å frembringe produkter med identiske eller forskjellige egenskaper.
«Bestanddeler» er ment å betegne alle de stoffer/substanser som er tilstede i reaktoren og er ment å inngå i syntesen eller å tillate den: ikke bare start-reagensene samt det/de syntetiserte produktet/produktene, men også eventuelle stoffer som ikke opplever noen omforming, slik som løsemidler, katalysatorer og lignende.
I tillegg til en eller flere reaktorer kan eventuelt anlegget hvor prosessen finner sted, innbefatte andre konvensjonelle anordninger slik som trykkreduseren-de ventiler, strippere, kondensatorer, tørkere, destillasjonskolonner og lignende. Disse hjelpeanordningene kan videre også generelt betraktes som reaktorer (perfekt blandede eller av pluggstrøm-type), selv om det ikke forgår noen kjemisk reaksjon i dem.
I tilfellet med en polymeriserings-prosess kan f.eks. «variablene som angår produktets egenskaper» velges blant molekylmassen, smelte-indeksen, standard densitet, komonomer-innhold når en komonomer er tilstede, og lignende.
Eksempler på «variabler relatert til prosessen» er spesielt temperaturen og trykket som hersker i reaktoren, prosessens produksjonshastighet, konsentrasjonene av de forskjellige reagenser i reaktoren, og lignende. Produksjonshastigheten betegner massen av syntetisert produkt pr. tidsenhet, hvilket imidlertid ikke nødvendigvis er det samme som strømningshastigheten for syntetisert produkt som forlater reaktoren: således er, som eksempel, og spesielt i oppstartings-fasene, massestrømmen av syntetisert produkt som forlater reaktoren, generelt svært liten, eller til og med null, selv om syntesen har begynt, dvs. at denne utløps-strømningshastigheten da er lavere enn produksjonshastigheten. På den annen side kan i en stabil tilstand produksjonshastigheten opptas i massen av syntetisert produkt pr. tidsenhet.
Eksempler på «manipulerte variabler» er strømningshastigheter for reagenser som går inn i reaktoren, effekten som mates til oppvarmingsanordningene, og lignende. Dette er variabler som gjør det mulig å innvirke på forløpet av prosessen, samt generelt på det syntetiserte produktets egenskaper.
Innstillingspunktet eller -punktene for den/de manipulerte variabel/variabler overføres direkte eller indirekte til konvensjonelle aktuatorer, slik som spesielt ventiler, oppvarmingsenheter og lignende. «Indirekte» betyr at de manipulerte variablene kan overføres gjennom en eller flere styringsenheter (som generelt styrer en enkelt variabel, f.eks. PID-kontrollere) som styrer aktuatoren eller aktuatorene («lokal» styring).
For så vidt angår maskinvare, så er generelt prediksjonsenheten og styringsenheten vanlige beregningsanordninger som muliggjør at beregninger kan utføres som en funksjon av deres kabling eller deres programmering; de kan spesielt være datamaskiner eller digitale systemer for ordre-styring (SNCC). En enkelt anordning kan med fordel kombinere prediksjons- og styringsfunksjonene. Beregningsanordningen eller -anordningene som benyttes, er fortrinnsvis av digital type, og leverer resultatene av sine beregninger med regelmessige intervaller (intermitterende). Tidsintervallene som atskiller utmatingen av disse resul- tåtene, kan variere med tiden, og kan også variere i samsvar med det resultat som inngår: det er klart at hurtig varierende variabler må beregnes på nytt hyppi-gere enn langsomt varierende variabler. Skiftregisteret kan anvendes for materiell simulering av dødtidene.
Prediksjonsenheten er basert på en direkte matematisk modell av prosessen (M), hvor reaktoren (R) er opptatt i en perfekt blandet reaktor; en eller flere rene forsinkelser (dødtider) kan eventuelt tas i betraktning for å representere eventuelle reaktorer av pluggstrøm-type, eventuelle transportforsinkelser eller forsinkelser ved oppnåelse av måleresultater osv.
Styringsenheten er fortrinnsvis basert på det omvendte av den direkte modell som anvendes i prediksjonsenheten (omvendt modell).
Generelt innbefatter summen Z Fxdis av alle massestrømmene (Fxdis) for hvilke bestanddelen X forsvinner fra reaktoren R, to ledd: - FRXlsom betegner den massestrøm ved hvilken X forbrukes i en eller flere mulige kjemiske reaksjoner; - Fxout, som betegner den mulige massestrøm av X som forlater reaktoren ved avtapping under reaksjonen, i det (vanlige) tilfellet hvor X ikke forbrukes fullstendig ved reaksjon i denne reaktoren; eller f.eks. ved for-damping, i tilfellet med en åpen reaktor. Fordelen ved fremgangsmåten er at Fxdis-leddene generelt er proporsjo-nale med Mxr,' f.eks. er
(trbetegner tilstedeværelses-tid i reaktor R)
(Rx betegner reaktiviteten for X i reaktoren R).
I dette tilfelle forenkles uttrykket som gir tx, og blir:
Dette uttrykket er uavhengig av MXr, og dette utgjør en spesielt fordelaktig forenk-ling.
En annen fordel ved fremgangsmåten ligger i den periodiske beregning
av tilstedeværelsestiden xx. xx representerer faktisk godt dynamikken for den bestanddel som betraktes i reaktoren. Dette gjør det mulig spesielt å følge forand-
ringen i denne parameteren, hvilket er viktig for å forstå prosessens dynamikk, og følgelig for styringen av den. På den annen side gjør de empiriske metodene av «black box»-type det ikke mulig å oppnå tilgang til denne parameteren.
Beregningen av prediksjonene av de styrte variablene (PGr) kan med fordel i tillegg ta i betraktning en eller flere målinger av styrte variabler (MGr), av manipulerte variabler (MGC) og/eller andre variabler relatert til kjøringen av prosessen (Ma<p>).
På lignende måte kan med fordel beregningen av innstillingspunktene for de prosess-manipulerte variablene (CGc) i tillegg ta i betraktning en eller flere målinger av styrte variabler (MGr), av manipulerte variabler (MGC) og/eller av andre variabler som er relatert til kjøringen av prosessen (Map) som er identisk med eller forskjellige fra de som eventuelt tas i betraktning med beregningen av prediksjonene av de styrte variablene (PGR).
Alle målinger som diskuteres i herværende beskrivelse, er ikke nødvendig-vis direkte målinger, i den forstand at en eller flere av dem kan eventuelt være målinger man kan slutte seg til, dvs. verdier oppnådd ved beregning ut fra en eller et antall andre direkte målinger. Således kan f.eks. ikke produksjonshastigheten for visse eksoterme synteseprosesser måles direkte, men en slutnings-måling av denne kan oppnås ved beregning, f.eks. ut fra (direkte) målinger av strømnings-hastigheten og av inngangs- og utgangs-temperaturene av kjøle-fluidet. 1 det spesielle tilfelle med polymeriserings-prosesser, velges fortrinnsvis den egenskap eller de egenskaper ved polymeren som er involvert i styringen, blant polymerens spesifikke densitet (SD), polymer-smeltens reologiské egenskaper, og dens komonomer-innhold.. Spesielt er den reologiské egenskap eller de reologiské egenskaper som inngår i styrings-fremgangsmåten, med fordel polymerens smelte-indeks og/eller en viskositetsmåling.
En eller flere egenskaper ved polymeren evalueres med fordel ved å anvende en teknikk valgt blant teknikkene nær infrarød-spektroskopi (NIR), Fourier transform infrarød-spektroskopi (FTIR) og kjernemagnetisk resonans (NMR).
Spesielt kan en eller flere egenskaper ved polymeren evalueres ved å anvende et på forhånd etablert korrelasjonsforhold med resultatene av målinger som utføres med nær infrarød spektroskopi (NIR) ved et antall bølgelengder forutbe stemt som en funksjon av polymerens egenskaper, og valgt mellom 0,8 og 2,6 mm.
Ytterligere detaljer vedrørende utførelsen av slike målinger i sammenheng med styring av polymeriserings-prosesser, kan finnes i patentsøknad EP328826 (US 5155184).
For å ta i betraktning mulige avvik mellom målingene og prediksjonene av de styrte variablene, kan det være nyttig å benytte seg av en korreksjon.
En første type korreksjon består i å korrigere innstillingspunktet for minst en styrt variabel (CGr) på grunnlag av det (med fordel filtrerte) avvik mellom målingen (Mgr) og prediksjonen (Pgr) av denne styrte variabelen, for å gjøre styringen effektiv (MGr = CGr) også i nærvær av en feil i prediksjonen av denne styrte variabelen. Denne teknikken kalles vanligvis «intern modell-styring» (IMC).
En andre type korreksjon består i periodisk å tilpasse modellen (M) av prosessen på grunnlag av det (med fordel filtrerte) avvik mellom prediksjonene (PGr) og målingene (MGR) av de styrte variablene, slik at også her modellen av prosessen bør levere prediksjoner av de styrte variablene (PGR) som er så nær som mulig (ideelt like) målingene av disse variablene (MGr), og dette er av avgjørende betydning for en effektiv styring.
Tilpasningen består i å rekalibrere modellen, dvs. i å beregne på nytt en eller flere av modellens parametere; normalt overskrider ikke antallet ombereg-nede parametere antallet styrte variabler som både en prediksjon og en måling er tilgjengelige for. En re-synkronisering (forskyvning i tid) av disse målingene er ofte ønskelig, fremfor alt når de er målinger av egenskaper ved det syntetiserte produktet som det tar lang tid å oppnå. Denne andre typen korreksjon er mer fordelaktig for så vidt den tillater modellen også å bli tilpasset med hensyn på sin dynamikk.
Tilpasningen angår ikke bare den direkte modellen av prosessen (prediksjonsenheten), men også den omvendte modellen (styringsenheten).
I samsvar med en fordelaktig alternativ form, er målingene (Mgr) av de styrte variablene bare involvert i den eventuelle tilpasning av modellen av prosessen, og er ikke direkte involvert i beregningen av innstillingspunktene for prosessens (Cgc) manipulerte variabler.
Dvs. at målingene av de styrte variablene ikke er involvert i den egentlige styringen: fordelen ved dette er at kvaliteten av styringen således ikke påvirkes av den eventuelle tregheten i evalueringen av produktets egenskaper.
Et annet aspekt ved oppfinnelsen er relatert til en styrings-fremgangsmåte slik som beskrevet ovenfor, anvendt på en polymeriseringsprosess, innbefattende ett eller f ler av de følgende tilleggstrinn: - beregning av et innstillingspunkt for temperatur i reaktoren som funksjon av ett eller flere innstillingspunkter for produktets egenskaper; og overføring av temperatur-innstillingspunktet til en eller flere aktuatorer som gjør det mulig å modifisere temperaturen i reaktoren (eventuelt indirekte, dvs. via en eller flere styringsenheter, f.eks. PID-kontrollere, som styrer aktuatoren/aktuatorene); - beregning av en varmebalanse for reaktoren, basert spesielt på temperatur-målinger; anvendelse av denne varmebalansen for å bestemme mengden syntetisert polymer pr. tidsenhet (produksjonshastighet) og/eller katalysator-virkningsgraden og/eller konsentrasjonen av minst en reagens i reaktoren; - beregning av varmemengden som produseres ved polymeriseringen, ved hjelp av en beregning av mengden reagens eller reagenser som polymeriserer;
bestemmelse ved hjelp av dette, av den varmemengde som må tilføyes eller fjernes for å opprettholde reaktor-temperaturen; anvendelse av resultatet av temperaturen denne beregningen (f.eks. ved fremad-mating) for å forbedre temperaturstyringen, for å gi så god tilpasning som mulig til innstillingspunktet for temperatur, spesielt i tilfelle av forandringer i produksjonshastigheten.
Disse alternative formene er basert på den relasjon som eksisterer mellom mengden reagens/reagenser som inngår i reaksjonen, og varmemengden som produseres eller absorberes ved reaksjonen.
I samsvar med en fordelaktig, alternativ form, beregnes egenskapen PxRved en bestanddel «x» i reaktoren R, som opptas i en perfekt blandet reaktor, som følger:
hvor «Px» er en egenskap ved en bestanddel «x», som i hovedsak tilsvarer den lineære blande-loven hvor Wi og w2er masse-fraksjonene av to blandede fraksjoner 1 og 2 av egenskap
Px1+2er egenskapen ved x når den forlater reaktoren etter blanding;
Pxiner egenskapen ved bestanddelen «x» når den går inn i reaktoren R;
Mxr er massen av bestanddelen x i reaktoren R;
Fxin er massestrømmen av bestanddel x som går inn i reaktoren R.
En matematisk transformasjon gjør det noen ganger mulig å gjøre noen variabler som ikke er lineære, lineære (additive): f.eks. adlyder smelte-indeksen for en polymer ikke noen lineær blande-lov, men dens logaritme; derfor utføres den ovenfor nevnte beregningen av Px1+2på denne parameterens logaritme.
I samsvar med en annen fordelaktig og alternativ form, innbefatter styrings-fremgangsmåten ifølge oppfinnelsen de følgende trinn: - innmating av innstillingspunkter som angår en eller flere egenskaper ved det produkt som skal syntetiseres, i en hoved-algoritme; - innmating av innstillingspunktet for prosessens produksjonshastighet i en slave-algoritme; - beregning av konsentrasjons-innstillingspunktene for bestanddelene i reaktoren ved hjelp av hoved-algoritmen, spesielt som en funksjon av innstillingspunktene og av målingene av produktets egenskaper og av målinger eller prediksjoner av de forskjellige bestanddelenes konsentrasjoner i reaktoren; - overføring av de konsentrasjons-innstillingspunkter som beregnes ved hjelp av hoved-algoritmen, som inngangsvariabler i slave-algoritmen; - beregning av strømningshastighets-innstillingspunkter for de bestanddeler som går inn i reaktoren, ved hjelp av slave-algoritmen, spesielt som funksjon av innstillingspunktet for prosessens produksjonshastighet, av konsentrasjons-innstillingspunkter og av strømningshastighets-målinger for de bestanddeler
som går inn i reaktoren, og
- overføring av de strømningshastighets-innstillingspunkter som beregnes ved hjelp av slave-algoritmen, til en eller flere aktuatorer (eventuelt indirekte, dvs. via en eller flere styringsenheter, f.eks. PID-kontrollere som styrer aktuatoren/aktuatorene) for å styre strømninghastighetene for de bestanddeler som går inn i reaktoren,
hvor hoved-algoritmen og/eller slave-algoritmen benyttes slik som beskrevet ovenfor, dvs. ved å anvende LAG-funksjonen til å beregne massen av minst en bestanddel i reaktoren.
Hoved- og slave-algoritmene benyttes også ved hjelp av en eller flere konvensjonelle beregningsanordninger. I samsvar med en fordelaktig, alternativ form, utføres alle beregninger (prediksjon, styring og lignende) med disse to algoritmene ved hjelp av samme beregningsanordning.
Målinger av temperatur (f.eks. temperatur i reaktoren og/eller inngangs-og/eller utgangs-temperaturen for et eventuelt kjølefluidum) inngår med fordel som ytterligere variabler for innmating i prediksjons- og/eller styringsenheten.
Slave-algoritmen tar også fortrinnsvis hensyn til målingene av sammensetningen av bestanddeler som er tilstede i eller forlater reaktoren.
I tillegg innbefatter styrings-fremgangsmåten fortrinnsvis et beregningstrinn med hjelp av slave-algoritmen, som funksjon av strømningshastighets-målingene, av prediksjoner av konsentrasjoner overført til hoved-algoritmen for å beregne prediksjoner av egenskaper som benyttes som ytterligere inngangsvariabler ved beregningen av konsentrasjons-innstillingspunktene.
Hoved- og slavealgoritmene utgjør en styring av kaskade-typen. Det er spesielt fordelaktig at hovedalgoritmen og/eller slave-algoritmen bør være adaptive, dvs. at noen av deres parametere bør om-beregnes periodisk (med regelmessige eller uregelmessige intervaller). En slik tilpasning gjør det mulig spesielt å garantere at den matematiske modellen beskriver prosessen i dens aktuelle tilstand så pålitelig som mulig, selv ved forandringer i noen driftsbetingelser (temperatur, trykk, produksjonshastighet og lignende) og i tilfelle av perturbasjoner (forgifting av katalysatoren etc).
Hovedalgoritmen utfører styringen av produktets egenskaper ved hjelp av en modell basert på karakteristiske ligninger som relaterer produktets egenskaper til de forskjellige bestanddelenes konsentrasjoner i reaktoren, samt eventuelt til temperaturen som hersker i reaktoren. Slave-algoritmen regulerer konsentrasjonene av en eller flere bestanddeler ved å innvirke på mate-strømningshastig-hetene for en eller flere bestanddeler, som kan være forskjellige.
Fordelen ved denne «hoved/slave«-kaskaden ligger i det faktum at hovedmodellen bestemmer med nøyaktighet konsentrasjonene av bestanddeler som er nødvendige for å oppnå egenskaper som er ønskelige for det syntetiserte produktet, og at slavemodellen sikrer at verdiene som påtvinges av hovedmodellen, adlydes. Slaven, som styres av hovedmodellen, er følgelig i stand til: - å bringe konsentrasjonene hurtig til de verdier som er ønsket av hovedmodellen, og holde dem der; og - å styre effektivt prosessens produksjonshastighet uten å forstyrre konsentrasjonene.
Denne hoved/slave-kaskaden er spesielt effektiv fordi både hovedmodellen og slavemodellen tar hensyn til prosessens dynamikk gjennom anvendelsen av LAG-funksjonen i beregningene.
Slavealgoritmen kan i tillegg utformes slik at den forsyner hovedalgoritmen med pålitelige prediksjoner av konsentrasjoner. Fra disse prediksjonene eller målingene av konsentrasjoner utleder hovedalgoritmen pålitelige prediksjoner av egenskaper ved produktet som er i ferd med å syntetiseres i reaktoren. Ved å sammenligne disse prediksjonene av egenskaper med innstillingspunktene for egenskapene, kan hovedalgoritmen, der det passer, gå inn og korrigere konsen-trasjonenes innstillingspunkter. Denne korreksjonen er mulig også før en uover-enstemmelse oppstår mellom en variabel og dens innstil I ingspunkt. Det faktum at prediksjoner av egenskaper, oppnådd fra prediksjoner eller fra målinger av konsentrasjoner, tas hensyn til, gjør det mulig redusere i betraktelig grad de fluktua-sjoner i tid i det syntetiserte produktets egenskaper, og dette resulterer i bedre ensartethet når det gjelder produktets kvalitet.
Hvis egenskapene ved det produkt som skal syntertiseres, avhenger av temperaturen i reaktoren/reaktorene, er det fordelaktig å tilveiebringe en tempe-raturstyring med slavealgoritmen. Denne stiller opp varmebalansen for hver reaktor, og ved hjelp av beregningen av produksjonshastigheten, bestemmer den den varmemengde som må tilføyes eller fjernes for å få tilpasning til de temperatur-innstillingspunkter som beregnes av hovedalgoritmen. Fra disse resultatene utleder den inngangs-innstillingspunkter for varme-styringsenheter i synteseutstyret. Man vil forstå at denne typen prosedyre gjør det mulig å gripe inn på varmestyringsenhetene i synteseutstyret også før temperaturen har foandret seg. Målinger av temperaturer går dessuten med fordel inn som ekstra inngangssignaler i slavealgoritmen.
Hovedalgoritmen innbefatter fortrinnsvis den følgende struktur:
- en prediksjonsenhet basert på en direkte modell av prosessen, som tillater at en prediksjon kan leveres av det syntetiserte produktets egenskaper som funksjon av målinger og/eller av prediksjoner av bestanddelenes konsentrasjoner; - en tilpasningsenhet som sammenligner prediksjonene av egenskaper, beregnet av prediksjonsenheten, med verdier som faktisk er målt må det syntetiserte produktet, og som utleder tilpasningsparametere fra denne sammenligningen, hvor de nevnte tilpasningsparametrene inngår som ytterligere inngangssignaler
til prediksjonsenheten for hovedalgoritmen; og
- en styringsenhet basert på en omvendt modell av prosessen, for å beregne som en funksjon av innstillingspunktene og av prediksjonene av egenskaper ved det produkt som skal syntetiseres, konsentrasjons-innstillingspunkter for slavealgoritmen, idet tilpasningsparametrene også inngår som ytterligere inngangssignaler til styringsenheten.
Slavealgoritmen innbefatter med fordel den følgende struktur:
- en prediksjonsenhet, basert på en direkte modell av prosessen, som tillater levering av en prediksjon av konsentrasjonene av en eller flere av bestanddelene, basert på en materialbalanse i reaktoren; - en tilpasningsenhet som sammenligner prediksjonene av konsentrasjoner, beregnet med den direkte modellen, med målinger av konsentrasjon, og som utleder tilpasningsparametere fra denne sammenligningen, hvor tilpasningsparameterne inngår som ekstra inngangssignaler til prediksjonsenheten i
slavealgoritmen; og
- en styringsenhet, basert på en revers-modell av prosessen, for å beregne, som en funksjon av innstillingspunktet for produksjonshastighet, av de konsentrasjons-innstillingspunkter som beregnes av styringsenheten i hovedalgoritmen, og av de konsentrasjonsprediksjoner som beregnes av prediksjonsenheten i slavealgoritmen, innstillingspunktene for de strømningshastigheter som går inn
j reaktoren, i det tilpasningsparametrene inngår som ytterligere inngangssignaler til styringsenheten i slavealgoritmen.
Prosessens dynamikk beskrives og beregnes med fordel ved hjelp av funksjoner av typen Y = LAG (u, x),
idet denne funksjonen er løsningen av differensialligningen
hvor argumentene u og x varierer med tid. Bruken av denne funksjonen i samsvar med teoremene 1 og 2 som oppstilles nedenfor, gjør det mulig å løse sekvensielt de material-balanser som benyttes av slavealgoritmen, og å beskrive prosessens kinetikk ved å anvende enkle karakteristiske ligninger i hovedalgoritmen. LAG-funksjonen gjør det videre mulig å redusere i betraktelig grad volumet av beregninger som er nødvendige, og gjør det følgelig unødvendig å benytte hurtige og kraftige datamaskiner. I tillegg gjør denne funksjonen det mulig å etablere de direkte og reverserte modellene av prosessen, eller av noen av prosessens deler, på en særlig enkel måte.
De viktigste egenskaper ved den foreslåtte styringen kan oppsummeres som følger: - forventning: styringen begynner å korrigere de målte forstyrrelsene/perturbasjonene også før virkningen av dem er sett på målingene av egenskapene (anvendelse av prediksjoner av konsentrasjoner, av prediksjoner av egenskaper og av prediksjoner av temperaturer i algoritmene); - nøyaktighet også i nærvær av perturbasjoner: den direkte modellen og den omvendte/reverserte modellen omkalibreres kontinuerlig ved anvendelse av
målingene av egenskapene (tilpasning);
- utvidet gyldighet: algoritmen beholder sin gyldighet under forandringene i produksjonshastighet og forandringer i kvalitet, og under oppstartinger og ved stopp (prosessens dynamikk uttrykt som en ligning, anvendelse av prediksjoner
for de variabler ved målinger som innebærer lange dødtider);
- enkelhet: utviklingen og implementeringen forenkles ved hjelp av en
original fremgangsmåte for å uttrykke prosessens dynamikk som en ligning (LAG-funksjonen).
Synteseprosessen som skal styres, modelldannes således på en form som generelt omtales som en «kunnskapsmodell» (første prinsipp-modell), dvs. at dens modell utvikles fra ligninger som reflekterer det detaljerte fysiokjemiske for-løp av prosessen. En slik tilnærmingsmåte gjør det mulig å oppnå, ved hjelp av en matematisk relativt enkel mengde ligninger, resultater som er overlegne i forhold til de som ville ha blitt oppnådd ved hjelp av en empirisk modell av «black box»-type, ved å levere spesielt parametre relatert til virkelige variabler, og bedre gyldighet utenfor identifikasjons-rommet (ekstra polering). De fleste av de empiriske modellene benytter komplekse ligninger, ofte av høy orden hvis det er ønskelig å oppnå en korrekt simulering av prosessens dynamikk, og hvor parametrene (spesielt tidskonstantene) må identifiseres for et nøyaktig driftspunkt; modellen er gyldig bare i det nærmeste området omkring dette driftspunktet. Det er vanskelig å gene-ralisere en slik tilnærmingsmåte til et høyt antall driftspunkter i tilfellet med en faktisk prosess med kjemisk syntese som involverer mange variabler.
I kontrast til dette, og i samsvar med styrings-fremgangsmåten ifølge oppfinnelsen, anvendes en mengde rent statiske, enkle ligninger; prosessens dynamikk simuleres ved hjelp av enkle funksjoner (jamfør LAG-funksjonen ovenfor). Med fordel kan oppholds-tidene (tidskonstanter for ligningene) beregnes på nytt så ofte det er ønskelig, og dette oppstiller ikke noe problem tatt i betraktning lig-ningenes enkelhet. I det endelige resultatet oppnås en mengde med ligninger som er ekstremt enkle og lette å løse i sann tid, også ved en høy frekvens.
Den foreslåtte styrings-fremgangsmåten kan anvendes med fordel på synteseprosesser, spesielt ved kontinuerlig syntese av polymerer (polymeirsering), og spesielt kontinuerlig polymerisering av olefiner slik som f.eks. etylen eller propylen, i flytende fase og i gassfase.
Det beskrives også også en prosess for syntes av ett eller flere kjemiske
produkter, styrt ved hjelp av en styrings-fremgangsmåte i samsvar med oppfinnelsen. Spesielt er det oppnådd svært gode resultater i tilfellet med styring av en prosess for kontinuerlig syntese av polyetylen ved polymerisering av etylen i minst en reaktor, idet reagensene innbefatter etylen, hydrogen og/eller en eventuell komonomer, hvor polymeriserings-reaksjonen finner sted i et løsningsmiddel i nærvær
av en katalysator, og en del av innholdet i reaktoren fjernes kontinuerlig eller intermitterende. Denne prosessen kan like gjeme finne sted i flytende fase og i gassfase; den foregår fortrinnsvis i flytende fase (i et løsningsmiddel).
Fremgangsmåten anvendes på lignende måte ved syntese av polypropylen (idet den viktigste start-monomer i dette tilfellet er propylen istedenfor etylen), idet det også er mulig at propan er tilstede hvis prosessen finner sted i gassfase.
I tilfellet med polypropylen omtales ofte smelte-indeksen som MFI istedenfor Ml.
Oppfinnelsen angår også en styringsanordning, ment for å iverksette styrings-fremgangsmåten ifølge oppfinnelsen, for syntese av ett eller flere kjemiske produkter, innbefattende en slik styringsanordning.
Mer nøyaktig angår oppfinnelsen også en anordning for styring av en synteseprosess for et kjemisk produkt i et synteseutstyr som innbefatter minst en reaktor, hvor anordningen innbefatter: - minst en beregningsenhet; - en anordning for innmating av innstillingspunkter for egenskaper ved det produkt som skal syntetiseres, i beregningsenheten; - en anordning for innmating av innstillingspunkter for produksjonshastighet for det produkt som skal syntetiseres, i beregningsenheten; - enheter for måling av strømningshastigheten for de strømmer som går inn i reaktoren; - enheter for måling av sammensetningen av de strømmer som forlater reaktoren; - enheter for styring av strømningshastigheter (aktuatorer) for å regulere strøm-ningshastighetene for strømmene som går inn i reaktoren; - anordninger for kommunikasjon mellom beregningsenheten, enhetene for måling av strømningshastighet og reguleringsenhetene;
hvor:
- massen av minst en bestanddel beregnes ved hjelp av LAG-funksjonen som er omtalt ovenfor; - beregningsenheten er i stand til å beregne ved hjelp av en hovedalgoritme, og som en funksjon av innstillingspunktene for egenskaper, innstillingspunktene for konsentrasjon av reagenser i reaktoren; - beregningsenheten er i stand til å beregne ved hjelp av en slavealgoritme, og som funksjon av innstillingspunktene for utmating og av innstillingspunkter for konsentrasjon, innstillingspunkter for strømningshastighet for de strømmer som går inn i reaktoren, idet disse innstillingspunktene for strømningshastighet overføres som inngangs-innstillingspunkter til enhetene for styring av strøm-ningshastigheter; - målingene som utføres av enhetene for måling av strømningshastigheter, inngår som .ytterligere inngangssignaler til slavealgoritmen, for å sette den sistnevnte i stand til å beregne prediksjoner av konsentrasjoner som funksjon av disse målingene av strømningshastigheter, idet disse prediksjonene av konsentrasjoner inngår i hovedalgoritmen for å beregne prediksjoner av egenskaper, anvendt som ekstra inngangssignaler i beregningen av innstillingspunkter for konsentrasjoner.
Oppfinnelsen angår også en styringsanordning slik som beskrevet ovenfor, hvor:
- synteseutstyret i tillegg innbefatter:
- enheter for varmestyring, i stand til å styre temperaturen i reaktoren, og
- temperatursensorer; - hovedalgoritmen er i stand til å beregne, som funksjon av innstillingspunktene for egenskaper, innstillingspunkter for temperaturer for reaktoren;
- slavealgoritmen er i stand til:
- å beregne en varmebalanse for reaktoren,
- å løse denne varmebalansen eller disse varmebalansene for å bestemme den varme som må tilføyes til, eller trekkes ut av, syntesen for å tilpasse til innstillingspunktene for temperatur, og - å utlede fra denne varmebalansen eller disse varmebalansene, inngangs-innstillingspunkter for enhetene for varmestyring av reaktoren, og - å motta, som ytterligere inngangsvariabler, de målinger som utføres av temperatursensorene.
Oppfinnelsen angår også en anordning slik som beskrevet ovenfor, hvor målinger som utføres temperatursensorene, inngår som ytterligere inngangssignaler i hovedalgoritmen.
Oppfinnelsen angår også en anordning slik som beskrevet ovenfor, innbefattende: - minst en analysator som er i stand til å levere målinger av de egenskaper som inngår i hovedalgoritmen; og - anordninger for å innmate disse målingene av egenskaper i beregningsenheten;
hvor denne beregningsenheten innbefatter:
- en første prediksjonsenhet basert på en første, direkte modell av prosessen som tillater prediksjon av egenskapene ved det syntetiserte produktet som funksjon av prediksjonene av konsentrasjoner beregnet av slavealgoritmen; - en første tilpasningsenhet som sammenligner prediksjonene av egenskaper, beregnet av den første prediksjonsenheten, med de verdier som faktisk er målt på det syntetiserte produkter, og som fra denne sammenligningen utleder tilpasningsparametre som inngår som ekstra inngangssignaler i den første prediksjonsenheten; og - en første styringsenhet, basert på en første reversert modell, for å beregne, som funksjon av innstillingspunktene og av prediksjonene av egenskaper, innstillingspunkter for konsentrasjoner for slavealgoritmen, idet tilpasningsparametrene også inngår som ekstra inngangssignaler i den første styringsenheten.
Oppfinnelsen angår også en anordning slik som beskrevet ovenfor, som i tillegg innbefatter: - minst en analysator som er i stand til å levere målinger av konsentrasjon av reagensene; og - anordninger for innmatinger av disse målingene av konsentrasjon, i beregningsenheten;
Idet denne beregningsenheten innbefatter
- en andre prediksjonsenhet, basert på en andre direkte modell som tillater prediksjon av konsentrasjonene som funksjon av material-balansen i reaktoren; - en andre tilpasningsenhet som sammenligner prediksjonene av konsentrasjoner, beregnet av den andre prediksjonsenheten, med målingene av konsentrasjon, og som fra denne sammenligningen utleder andre tilpasningsparametre som inngår ytterligere inngangssignaler i den andre prediksjonsenheten; og - en andre styringsenhet, basert på en andre, reversert modell, for å beregne, som funksjon av innstillingspunktene for utmating, innstillingspunkter for konsentrasjon beregnet av hovedalgoritmen og prediksjoner av konsentrasjon fra den andre prediksjonsenheten, innstillingspunkter for de strømningshastigheter som går inn i reaktoren, idet de andre tilpasningsparameterne inngår som ekstra inngangssignaler i den andre styringsenheten.
Oppfinnelsen er angitt i de vedføyde patentkravene.
En konkret utførelsesform av oppfinnelsen er illustrert på basis av en prosess for kontinuerlig syntese av polyetylen (PE), med henvisning til fig. 1 til 10. Disse figurene viser: Fig. 1: et diagram over en krets for fremstilling av polyetylen; Fig. 2: et forenklet diagram over strukturen av en avansert styring i samsvar
med oppfinnelsen; Fig. 3: et grunnleggende diagram over den avanserte styringen anvendt på
fremstillingskretsen ifølge fig. 1; Fig. 4: et grunnleggende diagram over en algoritme for adaptiv styring, an
vendt i det avanserte styringssystemet ifølge fig. 2; Fig. 5: et diagram over strukturen av hovedalgoritmen i det avanserte styrings
systemet ifølge fig. 2; Fig. 6: et diagram over strukturen av slavealgoritmen i det avanserte styrings
systemet ifølge fig. 2; Fig. 7: det generelle diagram over en styrings-fremgangsmåte i samsvar med
oppfinnelsen; og
Fig. 8-10: diagrammet over spesielle, alternative former av fremgangsmåten ifølge oppfinnelsen.
I fig. 7 kan man først se den egentlige synteseprosessen (Pr), som kan styres ved å anvende minst ett innstillingspunkt for en manipulert variabel (CGC)
(f.eks. en eller flere strømningshastigheter for bestanddeler som går inn i reaktoren, en temperatur og lignende) på en egnet aktuator (ventil, varme- eller kjølean-ordning og lignende). Styringen utføres ved hjelp av en styringsenhet (OC) basert på den omvendte matematiske modell av prosessen, og hvis hovedfunksjon det er
å sammenligne innstillingspunktet/innstillingspunktene for de styrte variablene (CGR) (f.eks. en eller flere egenskaper ved det produkt som skal syntertiseres, og/eller en eller flere variabler relatert til forløpet av prosessen) med prediksjonen/prediksjonene for disse variablene (PGR). Prediksjonen eller prediksjonene av de styrte variablene (PGR) beregnes av en prediksjonsenhet (OP) basert på den direkte matematiske modell for prosessen, på grunnlag av målinger av de manipulerte variablene (MGC). Det skal bemerkes at ingen måling av en egenskap eller egenskaper ved det syntetiserte produktet inngår i styringen. Fig. 8 viser en alternativ form av fremgangsmåten ifølge fig. 7, hvor den matematiske modellen for prosessen periodisk tilpasses ved hjelp av en tilpasningsenhet (OA) på grunnlag av avviket (fortrinnsvis filtrert eller behandlet digitalt) mellom prediksjonene (PGR) og målingene (MGR) av de styrte variablene. En resynkronisering (forskyving i tid) av disse målingene og av disse prediksjonene er ofte nødvendig, f.eks. når de angår målinger av egenskaper ved det syntetiserte produktet som det tar lang tid å oppnå. Tilpasningsenheten (OA) sender resultatene av sine beregninger, dvs. dens tilpasnings-instruksjoner, til prediksjonsenheten (for å tilpasse den direkte modell for prosessen) og til styringsenheten (for å tilpasse den reverserte modellen for prosessen). Det skal bemerkes at målingen/målingene av det syntetiserte produktets egenskap/egenskaper bare tas med i betraktningen i tilpasningsprosessen, som generelt finner sted med en mye lavere hyppighet enn den normale styringsprosessen. Disse målingenes eventuelle langsomhet har derfor ingen direkte virkning på kvaliteten av styringen. Fig. 9 viser en annen alternativ form av oppfinnelsen, hvor en eller flere målinger av styrte variabler (MGR) tas i betraktning av styringsenheten (OC), og en eller flere målinger av styrte variabler <MGR') (eventuelt forskjellige) tas i betraktning av prediksjonsenheten (OP). På lignende måte kan i tillegg en eller flere målinger av manipulerte variabler (MGC') tas hensyn til av styringsenheten (OC).
Det er åpenbart at det vil være mulig å skape en annen alternativ form av fremgangsmåten ifølge oppfinnelsen, ved å kombinere de alternative former i fig. 8 og 9, dvs. ved samtidig å benytte en tilpasningsenhet og ved å ta i betraktning en eller flere målinger av manipulerte variabler i styringsenheten, og/eller en eller flere manipulerte variabler i prediksjonsenheten og/eller i styringsenheten.
I fig. 10 er ikke prosessens matematiske modell tilpasset i ordets egentlige forstand, men avviket (fortrinnsvis filtrert) mellom målingene og prediksjonene for de styrte variablene, benyttes til å korrigere innstillingspunktene for de styrte variablene (Cgr). Under disse forhold har denne korreksjonen her blitt vist som en enkel differans: et korrigerende ledd beregnet ved hjelp av tilpasningsenheten OA (som i realiteten her bare er en korreksjonsenhet) subtraheres fra hvert av innstillingspunktene for styrte variabler, og dette gir korrigerte innstillingspunkter Cgr', overført til styringsenheten OC. Det er åpenbart at i noen tilfeller kan korreksjonen involvere operasjoner som er mer komplekse enn en subtraksjon, f.eks. en divi-sjon (i dette tilfellet er det imidlertid mulig å gå tilbake til en subtraksjon ved å betrakte logaritmene av de variabler som behandles). Denne fremgangsmåten kalles vanligvis indre modell-styring (IMC, Internal Model Control).
Hvis det vises til fig. 1, som i form av et diagram viser en krets for kontinue-lig syntese av polyetylen (PE), så foregår polymeriseringen av etylen i en sløyfe-reaktor 10, i suspensjon i et egnet løsemiddel slik som f.eks. heksan. Prosessen er kontinuerlig, dvs. at reagensene injiseres kontinuerlig, og at en del av reaktorens 10 innhold kontinuerlig tappes av. En sirkulasjonspumpe (ikke vist) sikrer homogenitet i reaktorens 10 innhold.
Reagensene som innføres i reaktoren, er etylen «Et», hydrogen «Hy» og buten «Bt» (jamfør henvisningstall 11). En katalysator injiseres også kontinuerlig. Det er viktig å ha god styring over konsentrasjonene av reagensene i reaktoren, fordi egenskapene til PE-harpiksen bestemmes hovedsakelig av konsentrasjonsforholdene Hy/Et og Bt/Et.
Polymeriserings-temperaturen i reaktoren er en ytterligere parameter som påvirker PE-harpiksens egenskaper. Siden polymeriserings-reaksjonen er svært eksotermisk, må reaktorens temperatur styres ved bruk av en kjølekrets 12.
Reaktoren 10 inneholder derfor under drift løsningsmiddel, polymer og reagenser som ikke har reagert ennå, samt katalysator. Innholdet i reaktoren tappes av kontinuerlig gjennom uttappingsrøret 14. Dette uttappings-innholdet går inn i en «stripper» STP 16, som atskiller PE-polymeren og fluidene (løsningsmiddelet og reagensene). Disse fluidene fordampes ved hjelp av injisering av vanndamp, og fjernes til en kondensator CD 18. I denne sistnevnte kondenseres løsnings- middelet på nytt, før det resykles. Reagensene som er lettere, atskilles fra løs-ningsmiddelet, og resykles også. En gassfase-kromatograf GC (20) plassert på utgangen fra kondensatoren 18, tillater at reagens-konsentrajonene Hy/Et og Bt/Et kan bestemmes.
Polymer som fjernes fra stripperen 16, konsentreres i en sentrifuge CFG 22, og tørkes så i en tørker SHLF 24 av fluidisert sjikt-type, før den sendes til av-slutning for granulering der. Prøver tas ved tørkerens 24 utgang for å måle egenskapene til harpiksen i den: krystallinitet (målt via spesifikk densitet «SD») og reologiské egenskaper (smelte-indeks (Ml) eller smelte-strømningstndeks (MFI) og smelte-viskositet «u2», målt i en kapilærstrømning under et skjær på 100 s"<1>).
Denne PE-synteseprosessens dynamikk er langsom og kompleks:
- Sløyfereaktoren 10 oppfører seg som en perfekt blandet reaktor. Følgelig vil enhver endring i mate-strømningshastigheten for en av reagensene, bare bli gradvis reflektert i denne reagensens konsentrasjon i reaktoren. Dette er fordi den nye strømningshastigheten må blandes med hele reaktorens 10 volum for
å bringe den til ny likevekts-konsentrasjon.
Målingen av konsentrasjonen av reagensene foretas av en gassfase-kromatograf 20; dette er et ikke-kontinuerlig instrument som arbeider i suksessive trinn: en gassprøve tas, analyse og deretter utledning av resultatene. Det er derfor en dødtid (fra 5 til 15 min.) mellom forandringene i konsentrasjon og målingen
av dem.
- Egenskapene ved den polymer som fremstilles, avhenger i ethvert øyeblikk av reagensenes konsentrasjoner. Enhver forandring i disse konsentrasjonene påvirker derfor øyeblikkelig egenskapene ved den polymer som fremstilles. På den annen side forandrer gjennomsnitts-egenskapene seg i reaktoren bare gradvis, fordi nylig produsert polymer må blande seg med den polymer som
allerede er tilstede i reaktoren 10 (tilstedeværelses-tid: ± 2 h).
- Når polymeren tappes av fra reaktoren 10, opplever den på nytt en rekke blandinger i de forskjellige utstyrsenheter (STP, CFG og SHLF), i den hensikt å tør-ke den (tilstedeværelses-tid: ±2 h). Prøver av polymer blir så tatt, og analysert av anleggs-laboratoriet. Resultatene av disse analysene vil derfor bare kunne meddeles til produsenten etter en ny dødtid, som kan være betraktelig ± 2 h).
Modelldannelse av prosessen med htelp av en LAG- funksjon
I samsvar med fremgangsmåten ifølge oppfinnelsen utføres den dynamiske modelldannelse av en kontinuerlig synteseprosess ved å ty til antakelser om perfekte blandinger og rene forsinkelser. De perfekte blandingene uttrykkes som en ligning ved hjelp av en funksjon som er velkjent for ingeniører, «LAG»-funksjonen, eller et lavpassfilter (av første orden); denne funksjonen er lineær og enkelt pro-grammerbar.
Den defineres på følgende måte: y = LAG (u, x)
(kjent som «LAG av u i løpet av x») som er løsningen av differensialligningen
hvor argumentene u og x varierer med tid.
Denne ligningen kan løses numerisk (også i sann tid) ved hjelp av en algebraisk ligning av første orden, som har de følgende variabler som argumenter:
- prøvetakingsperioden «T» (eller forløpt tid siden siste iterasjon).
- tilstedeværelsestiden (eller «tidskonstanten») «x» i øyeblikket «T».
- tilstandsvariabelen «y» i det foregående øyeblikk «t-T».
- styringsvariabelen «u» i nåværende øyeblikk «t» (u og x representerer i realiteten de målte eller beregnede verdier i øyeblikket «t» for variablene «u» og «x», som antas å ha vært konstante gjennom det forutgående intervall «T»).
T er fortrinnsvis liten i forhold til x (f.eks. minst 10 ganger mindre), for å øke beregningens nøyaktighet.
Løsningen av den ovenstående ligningen kan f.eks. beregnes ved hjelp av følgende formel: eller enklere (omtrentlig):
Modelldannelse av prosessen ved hjelp av LAG-funksjonen baserer seg på de følgende teoremer: La det være en perfekt blandet reaktor (CSTR) med volum Vr. Denne mates med forskjellige bestanddeler (reagenser eller inerte substanser), innbefattende reagensen «x» (strømningshastighet inn Fxin) som har egenskapen «Px!N» på inngangen. Strømningshastigheten ut «Fout» (uttapping) måles også.
Teorem 1: Anvendelse av LAG-metoden ved beregning av en massebalanse:
I hvert øyeblikk er massen «MxR» av en bestanddel «x» i en perfekt blandet reaktor (CSTR) lik LAG av produktet av inngangs-massestrømmen «Fx)N» og en tid «ix», i løpet av denne samme tiden «tx»:
Tiden «ix» er «tilstedeværelsestiden for «x»; dens verdi er massen av bestanddel MxRdelt på summen av de «utgående» massestrømmer (mengde for brukt ved reaksjonen «FRx», strømningshastighet som forlater reaktoren «Fxout», osv.). Dette teoremet gir således en nøyaktig metode for dynamisk beregning (også i sann tid) av konsentrasjonene i en perfekt blandet reaktor. I realiteten har, hvis «Vr» er reaktorens volum, konsentrasjonen «CxR» av bestanddelen «x» verdien
Siden den totale volumstrømmen av «FV0ut» som forlater reaktoren, er kjent, er reaktorens tilstedeværelsestid «tr» definert:
Massestrømmen «Fx0ut» for bestanddel «x» som forlater reaktoren (avtapping) er således:
Det skal bemerkes at hvis «x» er en inert substans (som ikke gjennomgår noen reaksjon, og forlater reaktoren bare ved avtapping), så gir dette:
I det hyppige tilfelle hvor reaksjonshastigheten for «x» er proporsjonal med dens konsentrasjon CxR, med en proporsjonalitetsfaktor Rx, gir dette videre: og følgelig
Teorem 2: Anvendelse av LAG-metoden ved beregning av en egenskap ved en blanding:
La det være en egenskap «Px» ved en bestanddel, som adlyder den lineære blandingsloven:
hvor Wi og w2er masse-fraksjonene for egenskap Pxiog Px2(med w-i + w2= 1).
I hvert øyeblikk er egenskapen PxRi en perfekt blandet reaktor (CSTR) lik LAG av egenskapen ved inngangen Px1funder en tilstedeværelsestid som har verdien masse-forholdet MxRsom er bestanddelen i reaktoren dividert på den inngående massestrøm (og/eller som fremkommer ved reaksjon) FxiN:
Det er således mulig å ta prosessens dynamikk med i beregningen, og å om-beregne dens tidskonstanter kontinuerlig.
Som forklart ovenfor kan de «egnenskaper» som er involvert her, i noen tilfeller ha gjennomgått en matematisk transformasjon som gjør dem lineære (f.eks. kan logaritmen av smelte-indeksen for en polymer betraktes som å adlyde en lineær blandingslov).
Prinsipp for styringen:
Når modellen for synteseprosessen er etablert, er det nødvendig med
en algoritme som betegner de nødvendige parametere for styringen av denne prosessen.
Fig. 2 viser et forenklet, generelt diagram over den type avansert styring (avansert prosesstyring eller «APC», Advanced Process Controll) som tas i bruk for polymeriserings-prosessen som er beskrevet ovenfor. Man ser at styringssystemet innbefatter en kaskade av to algoritmer, hvor denne kaskaden spesielt styrer PID-kontrollerne for reagensmatingens strømningshastigheter.
De to algoritmene i kaskade også kalt hovedalgoritme og slavealgoritme, er begge adaptive, dynamiske algoritmer basert på modeller med opphav i kunnskap om prosessen (i kontrast til empiriske modeller), som er basert spesielt på den styrte prosessens material-balanser og kinetikk. De anvender den på forhånd definerte LAG-funksjonen.
I fig. 3, som viser prinsippet for styringssystemet i sammenheng med den polymeriseringsprosess som er beskrevet ovenfor, ses det at: - hovedalgoritmen er basert på katalysatorenes karakteristiske ligninger, dvs.
ligninger som angir egenskapene ved PE som funksjon av polymeriserings-temperaturen og av reagensenes konsentrasjoner i reaktoren; den gir slave-algoritmen innstillingspunktene for reagensenes konsentrasjoner for å tilfredsstille innstillingspunktene for egenskapene til PE; - slavealgoritmen er basert på en material-balanse og reaksjonenes kjemiske kinetikk; den gir kontrollerne for mate-strømningshastighet innstillingspunktene for reagensenes strømningshastigheter, som er nødvendige for å tilfredsstille innstillingspunktene for konsentrasjoner, pålagt av hovedalgoritmen, og innstillingspunktet for prosessens produksjonshastighet. Fortrinnsvis beregner den også et antisiperende («feed forward») ledd for innstillingspunktet for temperatur, og forbedrer temperaturstyringen under forandringer i produksjonshastigheten.
Denne typen styring er bare fullstendig nøyaktig dersom modellen er perfekt og tar hensyn til alle mulige perturbasjoner. Generelt er dette ikke tilfelle.
Dette er grunnen til at (se fig. 4) den direkte modellen (og den reverserte modellen) generelt tilpasses kontinuerlig ved å sammenligne prediksjonene med målingene av egenskapene. Denne «tilpasningen» av modellen tillater at dens nøyak-tighet kan opprettholdes i nærvær av perturbasjoner som ikke er modelldannede, og gjør det således mulig å oppnå en nøyaktigere styring under alle forhold.
Slavealgoritmen
Fig. 6 illustrerer prinsippet for slavealgoritmen:
1. en prediksjonsenhet som benytter en direkte modell av den predikterte prosessen, fra målingene av strømningshastighetene som mater reaktoren, og reagens- og polymer-konsentrasjonene; 2. en tilpasningsenhet sammenligner konsentrasjonene av etylen (Et), hydrogen (Hy) og buten (Bt) målt av en analysator (gasskromatograf), med verdiene som er prediktert av den direkte modellen, for å bestemme tre tilpasningsparametere:
- den spesifikke aktivitet av katalysatoren for etylen «AET», i kg/h polyetylen
pr. kg katalysator og pr. kg/cm<3>etylen
- forsterkningsfeilen ved målingen av strømningshastighet av hydrogen «KfHy»
- renheten av buten-matingen «KfBt»; 3. styringsenheten beregner innstillingspunktene for reagensenes mate-strøm-ningshastigheter, fra innstillingspunkter for konsentrasjon, beregnet av hoved-algoritmen, og fra innstillingspunktet for produksjonshastighet; disse innstillingspunktene for strømningshastighet består av en fremovermating (feed forward) basert på den reverserte modellen, og av en tilbakekopling som er proporsjonal med avviket mellom den direkte modellen og innstillingspunktene for konsentrasjon.
For å forstå beregningene som utføres av slavealgoritmen i tilfellet med en prosess for syntese av polyetylen, må det aller først være kjent at det generelt aksepteres at polymeriserings-hastigheten «VitPEt» er proporsjonal:
- med konsentrasjonen av upolymerisert ET cEtR(i kg/m<3>),
- med konsentrasjonen av aktiv katalysator i reaktoren cCAR(i kg/m3), og
- med en proporsjonalitetsfaktor, en funksjon (som er vanskelig å kvantifisere) av temperaturen, av konsentrasjonene av Hy, av Bt og av sam-katalysator, av nærvær av gifter og lignende. Dette leddet benevnes «katalytisk aktivitet» for etylen AEt. I fravær av en større feilfunksjon (gift og lignende) varierer den relativt lite i løpet av en økt.
Hvis Vr er reaktorens volum, er derfor mengden etylen «FpEt» som polymeriserer i reaktoren for hver tidsenhet («polymeriserings-strømningshastighet»):
hvor MEtRer massen av etylen i løsning i reaktoren (i kg).
Det er videre kjent at inntakshastigheten for Hy er omtrent 100 ganger langsommere enn for Et og 10 ganger langsommere for Bt. Dette gir:
hvor MHyRer massen av hydrogen i løsning i reaktoren (i kg), og MBtRer massen av buten i løsning i reaktoren (i kg)
Slavemodellen benytter nå de følgende målinger:
FEtiN= mate-strømingshastighet for etylen (monomer) (kg/h) FSvin= mate-strømningshastighet for løsningsmiddel (heksan) (kg/h) FCAin = mate-strømningshastighet for katalysator , (kg/h) FHyiN= mate-strømningshastighet for hydrogen (transportmiddel) (kg/h) FBtiN= mate-strømningshastighet for buten (komonomer) (kg/h)
Den benytter også de følgende tilpasningsparametere:
AEt = «katalytisk aktivitet» for etylen
KfHy = forsterkningsfeil ved målingen av hydrogen-matingen
KfBt: = renhet av innmatet buten.
De følgende beregningene utføres sekvensvis, i følgende rekkefølge, med høy hyppighet (tidsintervallet som atskiller hver iterasjon er lite sammenlignet med den korteste av tilstedeværelsestidene xx).
Siden reaktorens volum er konstant, er utgående strømningshastighet pr. volumenhet ekvivalent med inngående strømningshastighet pr. volumenhet (ikke-kompresible fluider). Utgående strømningshastighet pr. volumenhet «Fv0ut» kan derfor beregnes som summen av inngående massestrømmer dividert på den densitet som de har reaktoren
(hvor densitetene er som følger: 650 kg/m<3>for løsningsmiddelet, 950 kg/m<3>for polyetylen, 600 kg/m3 for buten). Det antas her at all etylen øyeblikkelig omformes til polyetylen, og strømningshastigheten for hydrogen og for katalysator (noen få kg) ignoreres.
Løsningsmiddelet er kjemisk inert, og forlater reaktoren bare ved å tappes ut fra reaktoren. Løsningsmiddelets masse «MSvR» i reaktoren beregnes ved å anvende teorem 1:
Katalysatoren deaktiveres med en tidskonstant «kd»; massen «MCAR» av den aktive katalysatoren i reaktoren beregnes som følger: og derfor er konsentrasjonen «cCAR» av aktiv katalysator i reaktoren:
Etylen «forlater» reaktoren i uttappings-strømmen og i polymeriserings-reaksjonen. Dens masse i reaktoren «MEtR» beregnes på følgende måte:
Den «rå» (ukalibrerte) masse av hydrogen i reaktoren «MHyRAw» beregnes på lignende måte:
Massen «MHyR», korrigert for å ta hensyn til forsterkningsfeilen ved målingen av hydrogen-matingen, er.
Forholdet Hy: Et i reaktoren er derfor:
HyEtR<=>MHyR/MEtR
Massen av «rå» buten «MBtRAw» beregnes på lignende måte:
Massen «MBtR», korrigert for å ta hensyn til renheten av buten-matingen, er:
Forholdet BtEt i reaktoren er derfor:
BtEtR= MBtR/MEtR
Som allerede påvist, er polymeriserings-strømningshastigheten «FpEt»
(øyeblikkelig produksjonshastighet):
Siden polymeren er inert, og ikke gjennomgår noen reaksjon, er dens masse i reaktoren «MPER»:
Strømningshastigheten for polymer som forlater reaktoren «FPEout» er derfor:
Tilpasning av slavemodellen
Tilpasningsblokken anvender en analysator (f.eks. en gasskromatograf «GC») for å oppnå målingene av konsentrasjonene i reaktoren av etylen «cEtoc», av hydrogen «cHyGc» og av buten «cBtGC» (f.eks. uttrykt i kg/m<3>). Disse målingene sammenlignes med de verdier som er predikter! av den direkte modellen, for å bestemme de følgende tre tilpasningsparametere: - den spesifikke aktivitet av katalysatoren for etylen «AEt», i kg/h polyetylen pr.
kg katalysator og pr. kg/m<3>etylen.
- forsterkningsfeilen ved målingen av hydrogen-strømningshastighet «KfHy».
- renheten av buten-matingen «KfBt».
Gasskromatografen levere prøve-målinger med en forsinkelse på omtrent 6 minutter.
Beregningen av den spesifikke aktivitet for etylen «AEt» følger av likhetene nedenfor:
Siden målingen cEtccer samplet og inneholder støy, er det ikke ønskelig å finne den i en LEAD; det er derfor en foretrukket løsning å anvende:
Det tas hensyn til forsinkelsen på 6 minutter ved målingen av cEtGcved å innføre to LAG-er på tre minutter hver, i modellens verdier, og den endelige formel oppnås:
Beregningen av økningen i hydrogen-strømningshastigheten «KfHy» følger av likhetene nedenfor:
Man tar hensyn til forsinkelsen på 6 minutter ved målingen av cHyGcved å innføre to LAG-er i rekke på 3 minutter hver, i modellens verdi, og den endelige formel oppnås:
På samme måte som for KfHy beregnes korreksjonsparameteren for buten-renhet «KfBt»:
Styringsalgoritmen har som inngangsparametere:
- innstillingspunktene for konsentrasjon som beregnes ved hjelp av hovedalgoritmen; nærmere bestemt innstillingspunktene for konsentrasjonsforholdene
cHYR/cEtR«HyEtsp» og cBtR/cEtR«BtEtSP» (i kg/kg)
- innstillingspunktet for prosessens produksjonshastighet FpEtSp, fastsatt av operatøren - innstillingspunktet for etylen-konsentrasjonen cEtSP, fastsatt av operatøren; og
- konsentrasjonene som beregnes av modellen.
Den beregner innstillingspunktene for mate-strømningshastighetene for reagensene FETSp, FCA Sp, FHYsp og FBtSp. Forskjellige algoritmer kan benyttes, innbefattende MBPC (Model Based Predective Control, modell-basert, prediktiv styring). De kan vanligvis finnes å bestå av en fremovermating basert på den reverserte modellen, og en tilbakekopling som er proporsjonal med avviket mellom den direkte modellen og innstillingspunktene for konsentrasjon.
Styring av etvlen- matinoen:
Fremovermating: verdi for å holde den aktuelle konsentrasjon ved like, basert på inversjon av den stasjonære verdi av følgende ligning: Tilbakekopling: proporsjonal med avviket mellom innstillingspunktet cEtSp og modellen
Innstillingspunkt:
Styring av katalvsator- matingen:
Fremovermating: verdi for å holde den aktuelle konsentrasjonen vedlike, basert på inversjon av den stasjonære verdi av følgende ligning:
Tilbakekopling: proporsjonal med avviket mellom innstillingspunktet FpEtSPog modellen, i samsvar med følgende formel:
Innstillingsverdi:
Styring av hvdrogen- matingen:
Fremovermating: verdi for å holde den aktuelle konsentrasjon vedlike, basert på inversjon av den stasjonære verdi av følgende ligning:
Tilbakekopling: proporsjonal med avviket mellom innstillingspunktet for.forholdet HyEtsp og modellen.
Innstillingsverdi:
Styring av buten- matingen (lignende som for hydrogen):
De foranstående ligningene oppsummerer slavealgoritmens ligninger. De utføres hvert 10 sekund av den digitale styringen og ordre-systemet (SNCC). Hovedalgoritmen
Fig. 5 illustrerer prinsippet ved hovedalgoritmen:
1. dens prediksjonsenhet (basert på en direkte modell) predikterer de viktigste egenskapene ved polymeren (Ml og SD); for dette formål benytter den målingen av polymeriserings-temperaturen, prediksjonene for konsentrasjonene i reaktoren, som leveres av slavemodellen, og tilstedeværelsestidene for PE i de forskjellige enhetene av utstyret; 2. dens tilpasningsenhet sammenligner målingene av Ml og SD, som utføres (også ved utgang fra tørkeren) enten av målelaboratoriet for hver annen time eller ved hjelp av en kontinuerlig analysator, med verdiene som predikteres av den direkte modellen, for å bestemme de to tilpasningsparametrene som er korreksjonsparametere, multiplikative i tilfellet med Ml, og aditive i tilfellet med SD; 3. dens styringsenhet (basert på en reversert modell) beregner innstillingspunktene for konsentrasjonene i reaktoren (forholdene Hy: Et og Bt:Et) ut fra innstillingspunktene for Ml og for SD, som operatøren leverer. Slik som i tilfellet med slavealgoritmen utgjøres denne beregningen av en fremovermating basert på den direkte modellen, og en tilbakekopling som er proporsjonal med avviket
mellom den direkte modellen og operatørens innstillingspunkter.
For en gitt katalysator er harpiksens egenskaper i en stabil tilstand funksjoner av polymeriserings-temperaturen og av reagensenes konsentrasjoner. Blant de forskjellige statiske ligninger som inngår i litteraturen, er de følgende ligninger valgt:
Parameterne a0til a3 og b0til b4 oppnås ved identifikasjon i en stabil tilstand, for et antall harpikser som fremstilles med samme katalysator.
Videre kan alle de forskjellige utstyrsenheter som påmøtes av polyetylen inntil det tidspunkt hvor dets egenskaper blir målt (reaktor, stripper, sentrifuge og så tørker) inngå i, i en første aproksimasjon, perfekt blandede reaktorer.
Hovedalgoritmen benytter følgende målinger som inngangsparametere: TR= temperatur i reaktoren (°C)
V3tp = væskevolum i stripperen (oppnådd ved nivåmåling) (m<3>) MImess = måling av Ml (smelte-indeks)
SDMess = måling av SD (spesifikk densitet)
og de følgende beregninger utføres av slavealgoritmen:
FpEt = øyeblikkelig utmating av polymer (produksjonshastighet) (kg/h) FPEout = strømningshastighet av PE som forlater reaktoren (kg/h) MPEr= masse av PE i reaktoren (kg)
HyEtR= forholdet mellom Hy og Et i reaktoren (kg/kg) BtEtR= forholdet mellom Bt og Et i reaktoren (kg/kg)
De rå, øyeblikkelige verdier (før tilpasning) av SD og av loaritmen til Ml («IMI») beregnes ved hjelp av:
De rå, gjennomsnittlige egenskaper ved utløpet fra reaktoren beregnes ved å anvende teorem 2:
Effektivt gjelder:
- egenskapene IMI og SD tilsvarer ganske godt en lineær lov for blandinger
- sløyfe-reaktoren kan settes lik en perfekt blandet reaktor
- massestrømmen av PE som «går inn i» (fremkommer i) reaktoren, er effektivt FpEt, mengden PE som polymeriserer i et hvilket som helst øyeblikk (produksjonshastighet).
Rå egenskaper ved målingen: siden det er kjent at det er omtrent 500 kg PE pr. m<3>i stripperen, og dersom det antas at stripperen er en perfekt blandet reaktor, beregnes de rå egenskapene ved stripperens utløp som følger:
Siden tilstedeværelses-tiden i sentrifugen er svært kort, kan den ses bort fra.
Tørkeren er en tørker av typen med fluidisert sjikt; den inneholder kontinuerlig omtrent 1400 kg PE. Det kan antas at nivået i stripperen ikke forandrer seg meget, og at strømmen som forlater den, er lik strømmen som går inn i den. Således er strømningshastigheten for PE som går inn i tørkeren, lik FPEout- Ved tørkerens utløp, stedet hvor prøven tas for målingen av egenskapen, er det så de følgende rå verdier:
IMIsh = LAG (IMIstp, 1400/FPEOut)
SDsh = LAG (SDstp, 1400/FPEOUt)
Egenskapene etter tilpasning oppnås ved å ta i bruk tilpasningsparametrene KMI (multiplikativ parameter) og kSD (additiv parameter); reaktorutgangens, stripperens og tørkerens egenskaper etter tilpasning er derfor:
Tilpasning av hovedalgoritmen:
Målinger av egenskaper tar en viss tid for utførelse ± 5 min med en in-line analysator, ± 1 h med utførelse ved hjelp av laboratoriet). For å muliggjøre at tilpasningsparametrene kan beregnes, er det derfor nødvendig å resynkronisere (forskyve i tid) de rå prediksjoner fra modellen med målingene. Dette kan f.eks. gjøres ved hjelp av et skiftregister (her kalt «DELAY function»):
med tmi og tsd= ± 5 min eller ±1 h, i henhold til hvorvidt målingen utføres av en kontinuerlig analysator eller av laboratoriet.
Ved hver ny måling av Ml eller av SD omberegnes den rå tilpasningsparameteren kMI' eller kSD' ved å sammenligne den resynkroniserte, rå modellverdien med den målte verdien:
Disse rå verdiene filtreres for å dempe de hurtige reaksjonene som de eventuelle perturbasjonene (støy) i målingene kan virke inn på prosessen med:
Styringsenhet:
Styringsenheten har som innstillingspunkter verdiene MISp og SDSP, inn-tastet av operatøren. Han beregner innstillingspunktene for forholdene mellom konsentrasjonene i reaktoren, HyEtSPog BtEtSp, som er nødvendige for å oppnå hurtig de ønskede egenskaper MISPog SDSP. Denne beregningen foretas i to trinn: 1. ut fra innstillingspunktene Mlsp og SDSPsom operatøren leverer, og ut fra verdiene etter tilpasning av Ml og av SD i de forskjellige utstyrsenhetene, beregner styringsenheten innstillingspunktene MliSPog SDiSPfor den øyeblikkelige utmatingen. Disse øyeblikkelige innstillingspunktene består av en fremovermating og av en tilbakekopling som er proporsjonal med avviket mellom den direkte modellen og operatørens innstillingspunkter; 2. innstillingspunktene for konsentrasjonsforholdene HyEtSPog BtEtSPberegnes så ved å invertere den statiske ligningen som benyttes ovenfor for beregning av den øyeblikkelige verdi av Ml og av SD.
Innstillingspunkter for de øyeblikkelige egenskapene: tørkerutgangs-egenskapene sammenlignes med egenskapenes innstillingspunkter for å bestemme de ønskelige innstillingspunkter for egenskapene ved stripperéns utløp (sentrifugen ses bort fra):
På tilsvarende måte beregnes de ønskede innstillingspunkter ved reaktorens utløp, ut fra avviket mellom disse stripperutløps-innstillingspunktene og de kalibrerte verdiene ved stripperen:
Endelig beregnes de ønskede innstillingspunkter for den øyeblikkelige utmatingen, ut fra avviket mellom disse reaktorutløps-innstillingspunktene og de tilsvarende kalibrerte verdiene:
Innstillingspunkter for konsentrasjonsforholdene: innstillingspunktene for konsentrasjonsforholdene HyEtSp og BtEtSp oppnås ved å invertere den statiske ligningen som benyttes ovenfor for beregningen av øyeblikkelig verdi av Ml og av SD, ved å substituere i Ml- og SD-leddene de ønskede innstillingspunkter for den øyeblikkelige utmatingen, og ved å anvende tilpasningsparameteren.
Med start fra:
oppnår man: som gir: som gir:
De ovenstående ligninger oppsummerer ligningene i hovedalgoritmen. De utføres for hver 30 sekunder ved hjelp av SNCC.
Med denne prosessen er det mulig å styre polymeriseringen med stor nøy-aktighet. Spesielt: - opprettholdes de styrte egenskapene (Ml og SD) så nær som mulig opptil de ønskede verdier, med minimal spredning - utføres forandringene i kvalitet (og følgelig i Ml- og SD-egenskaper) med hur-tighet og nøyaktighet - utføres polymeriserings-oppstarting og -stopping, samt forandringer i prosessens produksjonshastighet på en hurtigere måte, mens Ml og SD holdes svært nær de ønskede verdier.
Selv om styrings-fremgangsmåten ifølge oppfinnelsen er presentert med hjelp av en prosess for syntese av polyetylen ved kontinuerlig polymerisering av etylen, skal det forstås at styringsmetoden vil være effektiv generelt for andre synteseprosesser, og spesielt for de prosesser som oppviser en eller flere av de følgende egenskapler: - en flervariabel styring er nødvendig fordi et antall variabler påvirker mengden av egenskaper som skal reguleres; - prosessens dynamikk er langsom: blandinger i serier, lange dødtider; - målingene av egenskaper foregår med prøvetaking med lav hyppighet og/eller med støy; - styringen må være dynamisk, dvs. gyldig uansett prosessens produksjonshastighet, samt gjennom forandringer i produksjonshastighet og i kvalitet
(egenskaper) for det produkt som skal syntetiseres;
- det er fordelaktig å estimere noen variabler som ikke måles direkte.
For å kunne implementeres enkelt med de teknikker som er fremvist her, er det tilstrekkelig at: - prosessens statiske ligninger er kjent (de er ofte kjent i det minste i en viss grad, ellers kan ikke prosessen styres); - prosessens dynamikk kan tilnærmes ved bruk av perfekte blandinger og død-tider; - de nødvendige målinger er tilgjengelige og har tilstrekkelig kvalitet (spesielt reagensenes strømningshastigheter og de strømningshastigheter som går gjennom de angjeldende lagringsbeholdere).
Den spesielle anvendelse av LAG-funksjonen som er beskrevet ovenfor, spesielt i teoremene 1 og 2, kan naturligvis utvides til styringsmetoder som er basert på en annen struktur enn den som er gjennomgått, som innbefatter en hovedalgoritme og en slavealgoritme som er forskjellige.
Den kan f.eks. anvendes i en styrt metode som omfatter bare en eneste algoritme.
Eksempler:
8 forsøk med syntese av polyetylen (PE) av 4 forskjellige typer (definert ved hjelp av Ml, SD, osv.) ble utført, henholdsvis ved anvendelse av en konven-sjonell styrings-fremgangsmåte og ved anvendelse av fremgangsmåten ifølge oppfinnelsen. Tabellen nedenfor oppsummerer de resultater som ble funnet på grunnlag av mange målinger av smelte-indeksen for de 8 oppnådde polymerer. Cpk betegner sentrert kapabilitetsverdi (centred capability value) for prosessen.
Det er funnet at kapabilitetsverdien Cpk mer enn dobles ved anvendelse av fremgangsmåten ifølge oppfinnelsen, hvilket viser at egenskapene bare er omtrent halvparten så spredt og/eller bedre sentrert i forhold til de pålagte innstillingspunktene.
Henvisninastall i figurene
Claims (17)
1. Fremgangsmåte for styring av en synteseprosess for minst ett kjemisk produkt i et utstyr som omfatter minst en reaktor (R) som kan tas opp i en perfekt blandet reaktor, hvor en eller flere manipulerte variabler (GC) tillater innvirken på prosessens forløp for å gjøre en eller flere variabler relatert til produktets egenskaper og/eller til prosessens forløp, og som benevnes styrte variabler (GR), lik de tilsvarende innstillingspunkter (Cgr), hvor fremgangsmåten innbefatter de følgen-de trinn: (a) innmating av innstillingspunkter vedrørende de styrte variablene (Cgr); (b) beregning ved hjelp av en prediksjonsenhet (OP), av prediksjoner av de styrte variablene (Pgr), basert på målinger av prosessens manipulerte variabler (Mgc); (c) anvendelse av en styringsenhet (OC) for å beregne innstillingspunktene for prosessens manipulerte variabler (Cgc). basert på innstillingspunktene (Cgr) og prediksjonene (Pgr) for de styrte variablene; (d) overføring av innstillingspunktene for prosessens manipulerte variabler (Cgc) til aktuatorer, eller til styringsenheter som styrer aktuatorene, for å innvirke på prosessens forløp;
hvor prediksjonsenheten (OP) er basert på en matematisk modell av prosessen,
kalt en direkte modell (M),
karakterisert vedat prediksjonsenheten (OP) er utformet på en slik måte at massen (MXr) av minst en bestanddel (X) i reaktoren (R) predikteres ved hjelp av ligningen:
hvor: - FxRiner massestrømmen av bestanddelen X som går inn i reaktoren R; - txer tilstedeværelsestiden for X i reaktoren, hvis verdi er hvor: - MXRbetegner den sist beregnede verdi av massen av bestanddelen X tilstede i reaktoren R; - Z Fxdis betegner summer av alle massestrømmer Fxdis ved hvilke bestanddelen X forsvinner fra reaktoren R, spesielt ved reaksjon og/eller ved å forlate reaktoren; - funksjonen y = LAG (u, x) er løsningen av differensialligningen
beregnet med den øyeblikkelige verdi av u og av x, og med den sist beregnede verdi av y.
2. Fremgangsmåte ifølge krav 1,
karakterisert vedat innstillingspunktet for minst en styrt variabel (Cgr) korrigeres på basis av avviket mellom målingen (Mgr) og prediksjonen (Pgr) for denne styrte variabelen, for å gjøre styringen effektiv også i nærvær av en feil i prediksjonen av denne styrte variabelen (Pgr).
3. Fremgangsmåte ifølge krav 1,
karakterisert vedat modellen (M) av prosessen periodisk tilpasses på grunnlag av avviket mellom prediksjonene (Pgr) og målingene (MGr) av de styrte variablene, slik at modellen av prosessen skal levere prediksjoner av de styrte variablene (Pgr) som er så nær som mulig opptil målingene av disse variablene (MGR).
4. Fremgangsmåte ifølge krav 3,
karakterisert vedat målingene (Mgr) av de styrte variablene bare inngår i den eventuelle tilpasning av modellen av prosessen, og ikke inngår direkte i beregningen av innstillingspunktene for prosessens manipulerte variabler (Cgc)-
5. Fremgangsmåte ifølge et av kravene 1 til 4, anvendt i en polymeriserings-prosess,
karakterisert vedat den innbefatter ett eller flere av de følgende tilleggstrinn: - beregning av et innstiliingspunkt for temperatur i reaktoren som funksjon av ett eller flere innstillingspunkter for produktets egenskaper; og overføring av dette temperatur-innstillingspunktet til en eller flere aktuatorer som gjør det mulig å modifisere temperaturen i reaktoren; - beregning av en varme-balanse for reaktoren, basert særlig på temperatur-målinger; anvendelse av denne varmebalansen for å bestemme mengden av polymer syntetisert pr. tidsenhet og/eller katalysator-effektiviteten og/eller konsentrasjonen av minst en reagens i reaktoren; - beregning av varmemengden som frembringes ved polymeriseringen, ved hjelp av en beregning av mengden av reagens eller reagenser som polymeriserer; bestemmelse ved hjelp av dette, av den varmemengde som må tilføyes eller fjernes for å opprettholde reaktorens temperatur; anvendelse av resultatet av beregningen for å forbedre temperaturstyringen, for å tilpasse så godt som mulig til innstillingspunktet for temperatur, spesielt i tilfelle av forandringer i produksjonshastighet.
6.. Fremgangsmåte ifølge et av de foregående krav,
karakterisert vedat egenskapen PxRved en bestanddel «x» i reaktoren R, som opptas i en perfekt blandet reaktor, beregnes som følger:
hvor «Px» er en egenskap ved en bestanddel «x», som tilsvarer hovedsakelig den lineære blandingsloven
hvor Wi ogW2er masse-fraksjonene for to blandede fraksjoner 1 og 2 av egenskapene Px-i og Px2i
Px1+2 er egenskapen ved x når den forlater reaktoren etter blanding;
PxiNer egenskapen ved bestanddelen «x» når den går inn i reaktoren R;
Mxrer massen av bestanddelen x i reaktoren R;
Fxiner massestrømmen for bestanddel x som går inn i reaktoren R.
7. Fremgangsmåte i et av de foregående krav,
karakterisert vedat den innbefatter følgende trinn: - innmating av innstillingspunkter som angår en eller flere egenskaper ved produktet som skal syntetiseres, i en hovedalgoritme; - innmating av innstillingspunktet for prosessens produksjonshastighet i en slavealgoritme; - beregning av innstillingspunktene for konsentrasjon av bestanddelene i reaktoren, med hovedalgoritmen, spesielt som funksjon av innstillingspunktene og av målingene av produktets egenskaper og målinger eller prediksjoner av konsentrasjonene av de forskjellige bestanddeler i reaktoren; - overføring av innstillingspunktene for konsentrasjon som beregnes av hoved-algoritmen, som inngangsvariabler til slavealgoritmen; - beregning av strømningshastighets-innstillingspunkter for bestanddelene som går inn i reaktoren, med slavealgoritmen, spesielt som funksjon av innstillingspunktet for prosessens produksjonshastighet, av konsentrasjons-innstillingspunktene og av strømningshastighets-målingene av bestanddelene som går inn i reaktoren, og - overføring av strømningshastighets-innstillingspunktene som beregnes med slavealgoritmen, til en eller flere aktuatorer for å styre strømningshastighetene for de bestanddeler som går inn i reaktoren,
hvor hovedalgoritmen og/eller slavealgoritmen benyttes i samsvar med et av de foregående krav.
8. Fremgangsmåte ifølge krav 7,
karakterisert vedat hovedalgoritmen innbefatter: - en prediksjonsenhet basert på en direkte modell av prosessen, som gjør det mulig å levere en prediksjon av egenskapene ved det syntetiserte produktet som en funksjon av målinger og/eller av prediksjoner av konsentrasjonene av bestanddelene; - en tilpasningsenhet som sammenligner prediksjonene av egenskaper beregnet av prediksjonsenheten, med verdier som faktisk er målt på det syntetiserte produktet, og som utleder tilpasningsparametre fra denne sammenligningen, hvor tilpasningsparametrene inngår som tilleggs-innmatingsverdier i prediksjonsenheten i hovedalgoritmen; og - en styringsenhet basert på en reversert modell av prosessen, for å beregne som funksjon av innstillingspunktene og av prediksjonene av egenskaper ved produktet som skal syntetiseres, konsentrasjons-innstillingspunkter for slavealgoritmen, hvor tilpasningsparameterne også inngår som ytterligere inngangs-verdier i styringsenheten.
9. Fremgangsmåte ifølge et av kravene 7 eller 8,
karakterisert vedat slavealgoritmen innbefatter: - en prediksjonsenhet basert på en direkte modell av prosessen, og som gir mulighet for levere en prediksjon av konsentrasjonene av en eller flere av bestanddelene, basert på en material-balanse i reaktoren; - en tilpasningsenhet som sammenligner prediksjonene av konsentrasjoner beregnet ved hjelp av den direkte modellen, med målinger av konsentrasjon, og som utleder tilpasningsparametere fra denne sammenligningen, hvor tilpasningsparametrene inngår som ytterligere innmatingsverdier i prediksjonsenheten i slavealgoritmen; og - en styringsenhet basert på en reversert modell av prosessen, for å beregne som funksjon av produksjonshastighets-innstillingspunktet, av konsentrasjons-innstillingspunktene som er beregnet av styringsenheten i hovedalgoritmen, og av prediksjonene for konsentrasjon som er beregnet av prediksjonsenheten i slavealgoritmen, innstillingspunktene for strømmene som går inn i reaktoren, hvor tilpasningsparametrene inngår som ekstra innmatingsverdier til styringsenheten i slavealgoritmen.
10. Fremgangsmåte ifølge et av de foregående krav,
karakterisert vedat den er anvendt ved styring av den kontinuerlige syntese av polyetylen ved polymerisering av etylen i minst en reaktor, polymeriser-ingsreaksjonen foregår i nærvær av en katalysator, og idét en del av innholdet i reaktoren kontinuerlig eller intermitterende fjernes.
11. Fremgangsmåte ifølge krav 9 eller 10,
karakterisert vedat tilpasningsenheten i slavealgoritmen sammenligner målingene av konsentrasjoner av etylen (Et) hydrogen (Hy) og/eller eventuell ko monomer (Bt), med verdiene som predikteres av prediksjonsenheten i slavealgoritmen, for å bestemme minst en av de følgende tilpasningsparametere: a) den spesifikke aktivitet av katalysatoren for etylen «AEt», i kg/h polyetylen pr.
kg katalysator og pr. kg/m<3>etylen; b) øknings-feilen i målingen av strømningshastighet for hydrogen «KfHy»; c) renheten av komonomer-matingen «KfBt».
12. Fremgangsmåte ifølge et av kravene 1 til 9,
karakterisert vedat den er anvendt ved styring av kontinuerlig syntese av polypropylen ved polymerisering av propylen i minst en reaktor, hvor reagensene innbefatter propylen, hydrogen og/eller en eventuell komonomer, polymeriser-ingsreaksjonen foregår i nærvær av en katalysator, og en del av reaktorens innhold fjernes kontinuerlig eller intermitterende.
13. Fremgangsmåte ifølge krav 9 og 12,
karakterisert vedat tilpasningsenheten i slavealgoritmen sammenligner målingene av konsentrasjonene av propylen (Pe), hydrogen (Hy) og/eller eventuell komonomer (Et), med verdiene som er prediktert av prediksjonsenheten i slavealgoritmen, for å bestemme minst en av de følgende tilpasningsparametre: a) den spesifikke aktivitet av katalysatoren for propylen «APe», i kg/h polypropylen pr. kg katalysator og pr. kg/m<3>propylen; b) øknings-feilen i målingen av hydrogen-strømningshastighet «KfHy»; c) renheten av komonomer-matingen (KfEt».
14. Fremgangsmåte ifølge krav 8,
karakterisert vedat den er anvendt på en polymeriseringsprosess, hvor - smelte-indeksen (Ml) og/eller standard densitet (SD) for polymeren og/eller dens komonomer-innhold måles periodisk; - prediksjonsenheten i hovedalgoritmen beregner rå prediksjoner av (Ml) og av SD som funksjon av temperaturen i reaktoren, av konsentrasjonen i reaktoren og av tilstedeværelses-tiden i de forskjellige utstyrsenheter i polymeriserings-kretsen; - periodisk utfører tilpasningsenheten i hovedalgoritmen følgende: - resynkronisering av de rå prediksjoner av Ml og SD, idet det tas hensyn til tiden som er gått mellom målingene foretatt av Ml og SD og oppnåelse av resultatet av målingene, og sammenligning av de resynkroniserte, rå prediksjoner av Ml og SD, med målingene av Ml og SD, - beregning av en multiplikativ tilpasningsparameter kMI anvendt i den rå prediksjon av Ml for å oppnå en kalibrert prediksjon av Ml, og - beregning av en additiv tilpasningsparameter kSD som anvendes i den rå prediksjon av SD for å oppnå en kalibrert prediksjon av SD.
15. Fremgangsmåte ifølge et av de foregående krav,
karakterisert vedat den er anvendt på en polymeriseirngsprosess, hvor en eller flere egenskaper ved polymeren evalueres ved å anvende en teknikk valgt blant nær infrarød spektroskopi (NIR), Fourier-transform infrarød spektroskopi (FTIR) og kjernemagnetisk resonans (NMR).
16. Fremgangsmåte ifølge et av det foregående krav,
karakterisert vedat den er anvendt på en polymeriseringsprosess, hvor en eller flere egenskaper ved polymeren evalueres ved å anvende et på forhånd etablert korrelasjonsforhold med resultatene av målinger utført med nær in-frarød spektroskopi (NIR) ved et antall forutbestemte bølgelengder, som funksjon av polymerens egenskaper, og valgt mellom 0,8 og 2,6 mm.
17. Anordning for styring av en synteseprosess ifølge ett av kravene 1-16, for et kjemisk produkt i et utstyr for syntese,
karakterisert vedat den innbefatter: - minst en reaktor som kan opptas i en perfekt blandet reaktor; - minst en anordning for innmating av et innstillingspunkt for en egenskap (CGR) ved produktet som skal syntetiseres, i beregningsenheten; - minst en anordning for å innmate et innstillingspunkt for produksjonshastighet for produktet som skal syntetiseres (CGC), i beregningsenheten; - minst en styringsenhet (OC); - minst en prediksjonsenhet (OP); - minst en anordning for å påtrykke en styrings-variabel (CGC) på en egnet aktuator.
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
BE9500506A BE1009406A3 (fr) | 1995-06-09 | 1995-06-09 | Methode de regulation de procedes de synthese de produits chimiques. |
PCT/EP1996/002408 WO1996041822A1 (fr) | 1995-06-09 | 1996-06-01 | Methode de regulation de procedes de synthese de produits chimiques |
Publications (3)
Publication Number | Publication Date |
---|---|
NO975780D0 NO975780D0 (no) | 1997-12-08 |
NO975780L NO975780L (no) | 1998-02-06 |
NO315563B1 true NO315563B1 (no) | 2003-09-22 |
Family
ID=3889027
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
NO19975780A NO315563B1 (no) | 1995-06-09 | 1997-12-08 | Fremgangsmåte for styring av kjemiske synteseprosesser |
Country Status (20)
Country | Link |
---|---|
US (1) | US6144897A (no) |
EP (1) | EP0830394B1 (no) |
JP (1) | JP2000500060A (no) |
KR (1) | KR100432252B1 (no) |
CN (1) | CN1155626C (no) |
AT (1) | ATE177761T1 (no) |
AU (1) | AU704563B2 (no) |
BE (1) | BE1009406A3 (no) |
BR (1) | BR9609240A (no) |
CA (1) | CA2220316C (no) |
CZ (1) | CZ293002B6 (no) |
DE (1) | DE69601794T2 (no) |
EA (1) | EA000343B1 (no) |
ES (1) | ES2131945T3 (no) |
HU (1) | HUP9803021A3 (no) |
NO (1) | NO315563B1 (no) |
PL (1) | PL186591B1 (no) |
RO (1) | RO120973B1 (no) |
UA (1) | UA37277C2 (no) |
WO (1) | WO1996041822A1 (no) |
Families Citing this family (79)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6106785A (en) * | 1997-06-30 | 2000-08-22 | Honeywell Inc. | Polymerization process controller |
US6122557A (en) * | 1997-12-23 | 2000-09-19 | Montell North America Inc. | Non-linear model predictive control method for controlling a gas-phase reactor including a rapid noise filter and method therefor |
GB2334958B (en) * | 1998-02-25 | 2001-11-07 | Porpoise Viscometers Ltd | Melt flow index determination in polymer process control |
US6093211A (en) * | 1998-04-09 | 2000-07-25 | Aspen Technology, Inc. | Polymer property distribution functions methodology and simulators |
BR9803848A (pt) * | 1998-10-08 | 2000-10-31 | Opp Petroquimica S A | Sistema para inferência em linha de propriedades fìsicas e quìmicas, sistema para inferência em linha de variáveis de processo, e, sistema de controle em linha |
DE69920954T2 (de) | 1998-11-30 | 2005-02-03 | Bp Chemicals Ltd. | Verfahren zur polymerisationskontrolle |
US6311134B1 (en) * | 1999-02-09 | 2001-10-30 | Mallinckrodt Inc. | Process and apparatus for comparing chemical products |
FR2790760A1 (fr) * | 1999-03-12 | 2000-09-15 | Bp Chemicals Snc | Procede de polymerisation de l'isobutene |
FR2794757B1 (fr) | 1999-06-11 | 2002-06-14 | Bp Chemicals Snc | Procede de polymerisation de l'isobutene |
EP1214362A1 (en) * | 1999-07-30 | 2002-06-19 | ExxonMobil Chemical Patents Inc. | Raman analysis system for olefin polymerization control |
US6479597B1 (en) | 1999-07-30 | 2002-11-12 | Exxonmobil Chemical Patents Inc. | Raman analysis system for olefin polymerization control |
US6983233B1 (en) | 2000-04-19 | 2006-01-03 | Symyx Technologies, Inc. | Combinatorial parameter space experiment design |
RU2270468C2 (ru) * | 2000-04-27 | 2006-02-20 | Бп Корпорейшн Норт Америка Инк. | Способ управления химическим технологическим процессом |
FR2810325A1 (fr) | 2000-06-16 | 2001-12-21 | Bp Chemicals Snc | Procede de polymerisation de l'isobutene |
WO2002008839A1 (en) * | 2000-07-21 | 2002-01-31 | Unilever N.V. | Method for simulating chemical reactions |
US6625512B1 (en) * | 2000-07-25 | 2003-09-23 | Advanced Micro Devices, Inc. | Method and apparatus for performing final critical dimension control |
DE10037153C2 (de) * | 2000-07-31 | 2003-04-10 | Franz Willeke | Vorrichtung zur Polymerisation in einem Schlaufenreaktor |
DE60141586D1 (de) * | 2000-08-07 | 2010-04-29 | Mitsui Chemicals Inc | Verfahren zur fertigungssteuerung |
US6760631B1 (en) | 2000-10-04 | 2004-07-06 | General Electric Company | Multivariable control method and system without detailed prediction model |
BR0116461A (pt) * | 2000-12-04 | 2003-09-23 | Univation Tech Llc | Processo de polimerização |
US6605675B2 (en) | 2000-12-04 | 2003-08-12 | Univation Technologies, Llc | Polymerization process |
US6996550B2 (en) * | 2000-12-15 | 2006-02-07 | Symyx Technologies, Inc. | Methods and apparatus for preparing high-dimensional combinatorial experiments |
FR2821175A1 (fr) * | 2001-02-19 | 2002-08-23 | Solvay | Methode de regulation d'une propriete d'un produit resultant d'une transformation chimique |
US6458916B1 (en) * | 2001-08-29 | 2002-10-01 | Hitachi, Ltd. | Production process and production apparatus for polybutylene terephthalate |
JP2003076934A (ja) * | 2001-09-03 | 2003-03-14 | Tosoh Corp | ポリマーの物性予測方法及びそれを用いたプラントの運転制御方法 |
CN1329110C (zh) | 2001-09-26 | 2007-08-01 | 伊内奥斯美国公司 | 一体化的化工工艺控制 |
CN100443868C (zh) * | 2001-11-09 | 2008-12-17 | 埃克森美孚化学专利公司 | 利用拉曼光谱分析的聚合物性能的在线测量和控制 |
US7799877B2 (en) | 2001-11-15 | 2010-09-21 | Univation Technologies, Llc | Polymerization monitoring and control using leading indicators |
US6657019B2 (en) * | 2001-11-20 | 2003-12-02 | Basf Corporation | Method and apparatus for predicting polymer latex properties in an emulsion polymerization process to improve the quality and productivity of the polymer latex |
US7846736B2 (en) * | 2001-12-17 | 2010-12-07 | Univation Technologies, Llc | Method for polymerization reaction monitoring with determination of entropy of monitored data |
US7226789B2 (en) * | 2001-12-17 | 2007-06-05 | Unication Technolofies, Llc | Method of applying non-linear dynamics to control a gas-phase polyethylene reactor operability |
DE10216558A1 (de) * | 2002-04-15 | 2003-10-30 | Bayer Ag | Verfahren und Computersystem zur Planung von Versuchen |
EP1359168A1 (en) * | 2002-05-03 | 2003-11-05 | BP Lavéra SNC | Rotomoulding polyethylene and method for producing same |
JP4112900B2 (ja) * | 2002-05-21 | 2008-07-02 | 株式会社山武 | 制御方法及び制御装置 |
EP1578814A2 (en) * | 2002-10-15 | 2005-09-28 | ExxonMobil Chemical Patents Inc. | On-line measurement and control of polymer properties by raman spectroscopy |
CA2409600A1 (en) * | 2002-10-24 | 2004-04-24 | Bayer Inc. | Process for the preparation of isoolefin polymers with controlled free-radical curability |
WO2004040283A2 (en) * | 2002-10-28 | 2004-05-13 | Bp Corporation North America Inc. | Control of a polymerization process |
US7319040B2 (en) * | 2002-10-28 | 2008-01-15 | Ineos Usa Llc | Process control using on-line instrumentation and process models |
US7106437B2 (en) | 2003-01-06 | 2006-09-12 | Exxonmobil Chemical Patents Inc. | On-line measurement and control of polymer product properties by Raman spectroscopy |
WO2005007711A1 (en) * | 2003-07-11 | 2005-01-27 | Innovene Europe Limited | Process for the (co-)polymerisation of ethylene in the gas phase |
EP1644423B1 (en) * | 2003-07-11 | 2006-12-13 | Ineos Europe Limited | Process for the (co-)polymerisation of ethylene in the gas phase |
US7402635B2 (en) * | 2003-07-22 | 2008-07-22 | Fina Technology, Inc. | Process for preparing polyethylene |
US7838605B2 (en) * | 2003-10-17 | 2010-11-23 | Univation Technologies, Llc | Polymerization monitoring and control using improved leading indicators |
US8058366B2 (en) * | 2003-10-17 | 2011-11-15 | Univation Technologies, Llc | Polymerization monitoring and method of selecting leading indicators |
US7400941B2 (en) | 2004-01-14 | 2008-07-15 | Chrevron Phillips Chemical Company Lp | Method and apparatus for monitoring polyolefin production |
WO2005100417A1 (en) * | 2004-04-13 | 2005-10-27 | Union Carbide Chemicals & Plastics Technology Corporation | Use of instantaneous split to improve reactor control |
US7483129B2 (en) | 2005-07-22 | 2009-01-27 | Exxonmobil Chemical Patents Inc. | On-line properties analysis of a molten polymer by raman spectroscopy for control of a mixing device |
US7505129B2 (en) | 2005-07-22 | 2009-03-17 | Exxonmobil Chemical Patents Inc. | On-line analysis of polymer properties for control of a solution phase reaction system |
US7505127B2 (en) * | 2005-07-22 | 2009-03-17 | Exxonmobil Chemical Patents Inc. | On-line raman analysis and control of a high pressure reaction system |
DE602005015220D1 (de) * | 2005-09-16 | 2009-08-13 | Mettler Toledo Ag | Verfahren zur Simulierung eines Prozesses auf Labormassstab |
US7738975B2 (en) | 2005-10-04 | 2010-06-15 | Fisher-Rosemount Systems, Inc. | Analytical server integrated in a process control network |
US8036760B2 (en) * | 2005-10-04 | 2011-10-11 | Fisher-Rosemount Systems, Inc. | Method and apparatus for intelligent control and monitoring in a process control system |
US7444191B2 (en) | 2005-10-04 | 2008-10-28 | Fisher-Rosemount Systems, Inc. | Process model identification in a process control system |
US20070255022A1 (en) * | 2006-04-28 | 2007-11-01 | Fina Technology, Inc. | Fluorinated transition metal catalysts and formation thereof |
US8110518B2 (en) * | 2006-04-28 | 2012-02-07 | Fina Technology, Inc. | Fluorinated transition metal catalysts and formation thereof |
US8138285B2 (en) * | 2007-10-26 | 2012-03-20 | Fina Technology, Inc. | Fluorinated impregnated catalyst systems and methods of forming the same |
ATE528326T1 (de) * | 2007-10-31 | 2011-10-15 | Basell Poliolefine Srl | Verfahren zur steuerung eines lösungsverfahrens zur olefinpolymerisation |
US8114353B2 (en) * | 2007-12-20 | 2012-02-14 | Chevron Phillips Chemical Company Lp | Continuous catalyst activator |
KR101737308B1 (ko) | 2008-12-29 | 2017-05-18 | 바셀 폴리올레핀 이탈리아 에스.알.엘 | 기상 중합 반응기 제어 |
CN101876822B (zh) * | 2009-04-30 | 2012-07-25 | 上海天坛助剂有限公司 | 聚氧乙烯醚生产装置的自动化控制方法 |
US9289739B2 (en) | 2009-06-23 | 2016-03-22 | Chevron Philips Chemical Company Lp | Continuous preparation of calcined chemically-treated solid oxides |
JP2013505489A (ja) | 2009-09-17 | 2013-02-14 | ビーエーエスエフ ソシエタス・ヨーロピア | 化学的エンジニアリング工程を制御する明確なスイッチを有する二自由度制御方法 |
EP2516488B1 (en) | 2009-12-22 | 2013-11-20 | Basell Polyolefine GmbH | Process for the preparation of ethylene homopolymers or copolymers in a high-pressure reactor controlled by a model based predictive controller |
CN103261241B (zh) | 2010-12-22 | 2016-05-11 | 巴塞尔聚烯烃股份有限公司 | 监测乙烯或乙烯与共聚单体在管状反应器中在高压下的聚合的方法 |
US9016075B1 (en) * | 2012-01-04 | 2015-04-28 | The Boeing Company | Aircraft environmental control system and method |
EP2874029A1 (de) * | 2013-11-15 | 2015-05-20 | Bayer Technology Services GmbH | Verfahren zum Betreiben einer zur Durchführung von wenigstens einer chemischen Reaktion eingerichteten Anlage |
CN104656457A (zh) * | 2013-11-19 | 2015-05-27 | 中国石油天然气股份有限公司 | 一种丙烯聚合环管反应器内部压力计算方法 |
US9389161B2 (en) | 2014-04-09 | 2016-07-12 | Exxonmobil Chemical Patents Inc. | On-line FT-NIR method to determine particle size and distribution |
CN105116931A (zh) * | 2015-07-21 | 2015-12-02 | 首钢京唐钢铁联合有限责任公司 | 一种氢氮混合站的氢气控制方法及*** |
IT201600081868A1 (it) | 2016-08-03 | 2018-02-03 | Versalis Spa | Metodo e sistema per il controllo di un impianto di produzione in continuo di un polimero. |
CN106647285B (zh) * | 2017-02-10 | 2019-08-09 | 泉州装备制造研究所 | 一种基于软测量技术的催化剂活性检测方法 |
EP3605249A1 (de) | 2018-08-02 | 2020-02-05 | Siemens Aktiengesellschaft | Verfahren zur synchronisierung, verfahren zum betreiben einer industriellen anlage, vorrichtung, computerprogrammprodukt und computerlesbares medium |
JPWO2020066309A1 (ja) * | 2018-09-28 | 2021-09-24 | 富士フイルム株式会社 | フロー反応設備及び方法 |
CN113260927B (zh) * | 2019-01-31 | 2024-06-18 | 陶氏环球技术有限责任公司 | 过程控制 |
ES2953941T3 (es) | 2019-04-01 | 2023-11-17 | Basf Se | Fabricación continua de prepolímeros de poliuretano |
CN111240284B (zh) * | 2020-01-17 | 2023-03-10 | 陈欣 | 一种模拟生产过程中液体浆料产品产量的品质控制方法 |
KR20220149062A (ko) * | 2021-04-30 | 2022-11-08 | 에스케이인천석유화학 주식회사 | 개질설비의 반응기 제어 방법 및 장치 |
CN113176767B (zh) * | 2021-04-30 | 2022-09-27 | 北京瑞飞华亿科技有限公司 | 一种丙烯回收控制*** |
EP4148035A1 (en) * | 2021-09-14 | 2023-03-15 | Air Liquide Societe Anonyme pour l'Etude et L'Exploitation des procedes Georges Claude | Methanol synthesis based on a mathematical model |
Family Cites Families (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4725654A (en) * | 1982-04-02 | 1988-02-16 | The Dow Chemical Company | Process for anionic polymerization |
IN160886B (no) * | 1983-04-25 | 1987-08-15 | Babcock & Wilcox Co | |
US4668473A (en) * | 1983-04-25 | 1987-05-26 | The Babcock & Wilcox Company | Control system for ethylene polymerization reactor |
US4572819A (en) * | 1983-08-24 | 1986-02-25 | The Dow Chemical Company | Apparatus for anionic polymerization wherein the molecular weight of the polymer is closely controlled |
US4921919A (en) * | 1985-12-10 | 1990-05-01 | Amoco Corporation | Method and apparatus for minimizing polymer agglomerate or lump formation in a gas-phase polypropylene polymerization reactor |
US5065336A (en) * | 1989-05-18 | 1991-11-12 | E. I. Du Pont De Nemours And Company | On-line determination of polymer properties in a continuous polymerization reactor |
EP0486262A1 (en) * | 1990-11-13 | 1992-05-20 | Vinamul Ltd. | Chemical processes |
ES2097520T3 (es) * | 1992-05-29 | 1997-04-01 | Amoco Corp | Polimerizacion de alfa-olefinas. |
JPH063243A (ja) * | 1992-06-23 | 1994-01-11 | Asahi Chem Ind Co Ltd | ゴム状重合体の物性推定方法 |
JPH0628009A (ja) * | 1992-07-07 | 1994-02-04 | Asahi Chem Ind Co Ltd | 重合プロセスの制御方法 |
JP3331031B2 (ja) * | 1993-12-15 | 2002-10-07 | 昭和電工株式会社 | ポリオレフィン重合反応器のスタートアップ制御方法および装置 |
US5570282A (en) * | 1994-11-01 | 1996-10-29 | The Foxboro Company | Multivariable nonlinear process controller |
-
1995
- 1995-06-09 BE BE9500506A patent/BE1009406A3/fr not_active IP Right Cessation
-
1996
- 1996-06-01 EA EA199800049A patent/EA000343B1/ru not_active IP Right Cessation
- 1996-06-01 AT AT96918665T patent/ATE177761T1/de not_active IP Right Cessation
- 1996-06-01 WO PCT/EP1996/002408 patent/WO1996041822A1/fr active IP Right Grant
- 1996-06-01 BR BR9609240A patent/BR9609240A/pt not_active IP Right Cessation
- 1996-06-01 AU AU61246/96A patent/AU704563B2/en not_active Ceased
- 1996-06-01 CZ CZ19973985A patent/CZ293002B6/cs not_active IP Right Cessation
- 1996-06-01 EP EP96918665A patent/EP0830394B1/fr not_active Revoked
- 1996-06-01 HU HU9803021A patent/HUP9803021A3/hu unknown
- 1996-06-01 CA CA002220316A patent/CA2220316C/fr not_active Expired - Fee Related
- 1996-06-01 JP JP9502575A patent/JP2000500060A/ja active Pending
- 1996-06-01 US US08/973,101 patent/US6144897A/en not_active Expired - Lifetime
- 1996-06-01 UA UA97126234A patent/UA37277C2/uk unknown
- 1996-06-01 ES ES96918665T patent/ES2131945T3/es not_active Expired - Lifetime
- 1996-06-01 DE DE69601794T patent/DE69601794T2/de not_active Revoked
- 1996-06-01 KR KR1019970709251A patent/KR100432252B1/ko not_active IP Right Cessation
- 1996-06-01 PL PL96323911A patent/PL186591B1/pl not_active IP Right Cessation
- 1996-06-01 RO RO97-02313A patent/RO120973B1/ro unknown
- 1996-06-01 CN CNB961959509A patent/CN1155626C/zh not_active Expired - Fee Related
-
1997
- 1997-12-08 NO NO19975780A patent/NO315563B1/no unknown
Also Published As
Publication number | Publication date |
---|---|
CZ293002B6 (cs) | 2004-01-14 |
ES2131945T3 (es) | 1999-08-01 |
AU6124696A (en) | 1997-01-09 |
WO1996041822A1 (fr) | 1996-12-27 |
NO975780D0 (no) | 1997-12-08 |
CA2220316A1 (fr) | 1996-12-27 |
DE69601794T2 (de) | 1999-10-07 |
BR9609240A (pt) | 1999-05-11 |
MX9709813A (es) | 1998-03-31 |
JP2000500060A (ja) | 2000-01-11 |
UA37277C2 (uk) | 2001-05-15 |
EA000343B1 (ru) | 1999-04-29 |
CA2220316C (fr) | 2008-11-18 |
AU704563B2 (en) | 1999-04-29 |
EP0830394B1 (fr) | 1999-03-17 |
EA199800049A1 (ru) | 1998-08-27 |
HUP9803021A3 (en) | 1999-05-28 |
PL186591B1 (pl) | 2004-01-30 |
US6144897A (en) | 2000-11-07 |
RO120973B1 (ro) | 2006-10-30 |
NO975780L (no) | 1998-02-06 |
EP0830394A1 (fr) | 1998-03-25 |
CN1192222A (zh) | 1998-09-02 |
CZ398597A3 (cs) | 1998-03-18 |
CN1155626C (zh) | 2004-06-30 |
ATE177761T1 (de) | 1999-04-15 |
KR100432252B1 (ko) | 2004-09-18 |
HUP9803021A2 (hu) | 1999-04-28 |
PL323911A1 (en) | 1998-04-27 |
DE69601794D1 (de) | 1999-04-22 |
KR19990022784A (ko) | 1999-03-25 |
BE1009406A3 (fr) | 1997-03-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
NO315563B1 (no) | Fremgangsmåte for styring av kjemiske synteseprosesser | |
Richards et al. | Measurement and control of polymerization reactors | |
Seki et al. | Industrial application of a nonlinear model predictive control to polymerization reactors | |
EP1259554B1 (en) | Process for the controlled production of polyethylene and its copolymers | |
NO322812B1 (no) | Fremgangsmate og anordning for regulering av en prosess | |
US4956426A (en) | Process for controlled polymerization of stereospecific alpha-olefins having preselected isotacticity | |
US10894840B2 (en) | Unified cooling for multiple polyolefin polymerization reactors | |
EP2167211A2 (en) | Process for the separation of a hydrocarbon-containing feed stream | |
JP2013505489A (ja) | 化学的エンジニアリング工程を制御する明確なスイッチを有する二自由度制御方法 | |
Finkler et al. | Realization of online optimizing control in an industrial semi-batch polymerization | |
NO178790B (no) | Fremgangsmåte ved fremstilling av olefinpolymerer i en autoklavrektor | |
CA3220830A1 (en) | Method for generating artificial intelligence model for process control, process control system based on artificial intelligence model, and reactor comprising same | |
FI88343B (fi) | Foerfarande och anordning foer beaktande av varierande volym och floede vid reglering av genomstroemningsprocesser | |
Bonvin et al. | Optimal grade transition for polyethylene reactors via NCO tracking | |
US10072103B2 (en) | Unified cooling for multiple polyolefin polymerization reactors | |
WO2002016932A1 (en) | Method of estimating the properties of a polymer product | |
US4526657A (en) | Control of a fractional distillation process | |
MXPA97009813A (en) | Method of regulation of quimi product synthesis procedures | |
Morales-Rodelo et al. | Determination and use of feasible operation region in flash distillation control | |
Said et al. | Tuning PID Controller for two-Continuous Stirred Tank Reactors in series with Time Delay | |
JPS6264808A (ja) | ポリオレフインの製造方法および装置 | |
Ogunnaike et al. | Control of polymerization processes | |
Othman et al. | Feedback Control of Industrial Solution Polymerization of Acrylic Acid Using NIR Measurements | |
JPH0348209B2 (no) | ||
Silva et al. | Optimization and nonlinear model predictive control of batch polymerization systems |