- [.,,- _. -
-... ,
.
.
-~~--.,...,....,-.-,....,.......,..."..-"----=,..--...,,.----~~~~~=----~-~~~~
tr' ~
<
u
,
.
.'. ,..
,0 i
0 '0-
/
.
Laboratorium voor Cher:nische Technologie
O' D~
·0· 0
I
.,
~-
!
,
~.
0
Verslag behorende bij het fabrieksvoorontwerp
0 U· 0 0 0
van
\r
...
.
.
....Q.A.v..Q.l.k ....e..tL.....~·····~fL.O'1.\."~.............
-
.......
onderwerp:
. d\.e.hb.~.~4L!.y.~.lb.e.~.e.. ... tnl.b.~v..... ~.. .~\A,.~,o.bh.e.{.~..c.....\I..e..~~.o..J......................:.....
0 0 0
.'
"
. adres:
vJo..te.\1 \0 0 ~~,. (l' cl~~
z~ o.Ch. tI\ St ~ 01.\ ~ \
~ J.~~ .
~/ttc..
t
opdrachtdatum :
"I.f t
verslagdatum :
D\Lt
.73 " '1
l
_
r' I l.
I.
Delft, oktober 1914
A.Valk Waterloostraat 119
,
,
•
l
, 1
Rotterdam R.Spronk Zwaanshals 220 Rotterdam
,. I
,
,
r '
a
c
,I
.
,-
L.
Inhoud: pag.
,.
r '
I
,I
Inleiding
1
Reakties
2
Kinetiek De autotherme methanol-reaktor
3 5
T.elichting bij het stroomschema
9
Verantwoording
12
Litteratuur
12
no. 1
,
1
stroomschema van het proces
2
blokschema,' s van de subroutines OUTPUT en INFO
3
blokschema's van de subroutines FCT en FIXBED
4
verklaring van gebruikte symbolen
5
stroomtabel
, ,
-1-
I
IL.
[' l.
l.
Aan de afdeling chemische technologie van de TH te Delft wordt in samenwerking met ingenieursbureau Co mprimo gewerkt aan een methode om met behulp van computerprogramma 's
r -
flowsh e ets t e berekenen.
Ruwweg is de gevolgde 'methode deze: van elk onderdeel (unit) in het processchema worden met behulp van de ingaande stromen de uitgaande stro Den berekend (variabe l en zijn bijvoorbeeld druk , temperatuur, molfracties et c. ) . Een coördinerend programma (EXEC 2) verzorgt de overdracht van gegevens en koppelt op de juiste wDze de te berekenen units. Units in het processchema die voor komen 1n geslo ten stroomkringen
(t.g.v. recycling) worden, om sluitende warmte- en massabalanzen te verkr Dgen , iteratief berekend. Het genoemde EXEC-programma stuurt o ok deze berekeningen. Daar geen dit
rekenpros=am~a
ontw or~en
voor ecn autotherme reactor b estond is
en v ervolgens getest in een voorontwerp van een deel
van een rnethanol -syn these -fabriek.
,
,
!'
-2-
Reakties
Methanol wordt gevormd volgens de onderstaande reaktiev erge lU king : (H
298
=
- 24 . 45 kcal/mol)
r .
De belancrjkste nevenreakt ies zUn : r
•
1)
2)
(H
3)
(H
29 8 298
- 49 .3 kcal/mol ) = +9.12 kcal/mol)
ad 1) Deze reaktie kan onderdrukt worden door gebruik t e maken van een geschikte katalysator zoals ZnO/Cr 0 • 2 3 ad 2) Door het sterk exotherme karakter van deze reaktie is h e t noo dzakelUk een maxinale temperatuur die in de reaktor berei kt kan worden , vast te l eg~en (4 00
ocl.
Gebeurt dit
niet dan z al de re aktie steeds sneller eaan verlopen door e en toene nende
warmte~roductie
wat kan leiden t ot sintering
van de katalysato r. ad 3) Deze reaktie, o ok wel shiftreaktie genoemd , is belangrijk vo or de warmtebeheersing in de reaktor door zijn endotherme karakter. Bovendien wordt een reaktant voor de methanol reak ti e gevormd •
• 1
•
1
r '
r •
I . L
l,
'-
-3-
•
Kinetiek
========
Voor de snelheid waarmee de methanolreaktie verloopt is een van de ,
volgende stappen bepalend. 1) Diffusie van de reaktanten naar het katalysator oppervlak. 2) Adsorptie aan het kata~sator op pervlak . r
•
3) Chemische reaktie. 4) Desorptie van het produkt. 5) Diffusie naar de bulkfase. volgens Natta (lit.1) is de chemische reaktie snelheidsbepalend terwijl hij aanneemt dat de adsoptie te beschrijven is vol gens Langmuir. Uitgaande van een derde orde oppervlakte-reaktie kwam hij tot de volgende reaktiesnelheidsvergelijking:
Hierin zijn A,B,C,D temperatuurfuncties:
A
O.0727.T + 93.1
B
-O.0425.T + 7.45 -O. 005.T + 2.1 -O.133.T + 56.2
e D
De reaktiesne lheidsvergelijking voor de shiftreaktie luidt:
De fugaciteiten op bepaalde afstand van de top der reaktor (A.) l
hangen op de volgende wijze af van de beginfugaciteiten (F.). l
AH20 AC02
FH20 + Y.F e02
ACO
(1-X)F CO
ACH30H ,
'
AH2
FC02 (1-Y)
X.F eo + Y.F e02 FH2 - Y.F e02 - 2AeH30H
Hierin zijn X en Yconversiegraden die gedefinieerd zijn als:
r'
r '
I
'
-4-
,
x
Y
l .
De beginfugaciteiten worden berekend met de relatie: F.1. = x.f. P 1 1 t
0
t
De temperatuurafhankelijkheid van d e evenwichtsconstanten is als volgt: log Kr K s
(5200/T - 12.66)
(lit.2)
exp(7.767/RT - 0.72/100R)
(lit.5)
Opmerking: de reaktiesnelheidsverge lijking gevonden door Natta is o bruikbaar in het temperatuurgebied van 300-400 C en bij drukken van 200-7 0 0 atm.
r e aktiesnelheid van de methanolreaktie reaktiesnelheid van de shiftreaktie evenwichtsconstant e methanolreaktie evenwichtsconstante shiftreaktie fugaciteit component i na bepaalde conversie
F.
fugac~teit
xi
molfractie component i
f.
fugaciteitscoëfficient
X,Y P
conversiegraden resp. CO, CO
T
temperatuur
A,B,C,D
adsorptieconstanten
R
gasconstante
k
fictieve reaktiesnelheidsconstante die het verloop van
1.
1
tot
I
I
, 1
component i aan begin reaktorbuis
2
druk in reaktor
I
de shiftreaktie over de gehele reaktor regelt
,
.
-5I
.
De autotherme me thanol-re aktor
======= = = = ====================
r'
1)
Algemene beschouwing (lit.3)
l.
In een autotherme reaktor (figuur I) stroomt het reaktantenmengsel langs buizen die gevuld zUn met katalysator en waarin zich de exotherme reakties afspelen. Door warmtewisseling met de buizen wordt het reaktantenmengsel opgewarmd. In de top van de reaktor wordt de stroomrichting van het gasmengsel omgekeerd waarna in de buizen de reaktie verloopt. De reaktiewarmteuitgezet te ge n de toptemperatuur geeft een beeld te zien zoals in figuur 11 is geschetst (kromme a). De vorm van de warmteproduktielUn kan verklaard worden met twee effecten.
1) BU lage temperatuur is de reaktiesnelheid klein. Verhoogt men de temperatuur dan zal de reaktiesnelheid exponentieel toenemen en daarmee de warmteproduktie.
2) BU ho g e temperaturen gaat de terug gaande reaktie een grotere rol spelen en neemt de warmteproduktieaf. De warmteafvoer wordt gegeven door de rechte lUn (b) waarvan de helling afhankelUk is van de warmtewisselcapaciteit; is deze groot dan is de helling groot en omgekeerd. stabi
operatie van een autotherme reaktor is slechts mogelUk in
die snUpunten van d e warmte produkt ie- en warmteafvoerlUn waar d( H d)/dT negatief is. pro Optimale stabiele operatie is mogelUk door het warmtewisselend vermogen zo te kiezen dat in de kromme een snUpunt wordt gevonden dicht r'
bU het maximum in de warmteproduktie-lUn;de gemiddelde reaktortemperatuur is dan relatief laag (gunstig i.v.m. sintering van de kalysator), zodat de conversie hoger kan zUn .
.,
. r
)
1
2)
Warmte- en massabalanzen • Voor de berekening van de reaktor worden balanzen opgesteld voor een volumeëlementje (dV) van een buis (figuur 111). Voor de hoeveelheid CO geldt:
~ (CO m
d
z+ z
- CO z ) = -R 1 ·dV
/
I. l.
,.
I .
~
z=L ------
-~-
"'"
T Ul' t
t// l' I
T, l n
fi g uur I
,
• 1
----+ Tt op
figuur 11
r'
r '
-6-
l.
Daar RNatta de omzetting per kg katalysator is, wordt de vergelUking
De inhoud van het volumeëlement is: dV = -Y-.r
2 0
2r
dz
Met invoering van een dimensieloze conversie:
wordt de vergelUking:
(1- E)~ ka t xr ~ .CO
m
2 figuur 111
0
o
Een identieke vergelUking beschrijft de CO
dX /dz 2
=
BoR /C0 2 2
o
0
0
0
0
2
0
0
0
(
1)
conversie in dV
(2 )
o
Daar bij het verlopen van de shiftreaktie CO gevormd wordt verandert vergelijking (1) :
Voor de berekening van het temperatuurverloop in de reaktor wordt een warmte balans ove r dV opgesteld: ~ oc (T
m
•
1
•
J
p
)=warmteproduktie - warmteafvoer door wand d -T 2 z+ z 2 • z 0
2
warmteproduktie
~ kat (1- E.)1tr . dz (R 1AH 1 + R2~H2)
warmteafvoer
U.~rodz(T2 o
z - T1 • z )
Stel: A = -~roU/c o~ p m Combinatie van de laatste vier vergelUkingen geeft:
•
0
•
0
•
(4)
De temperatuur van het reaktantenmengsel wordt alleen bepaald door de warmteoverdracht via de buiswand.De warmtebalans nodig voor het
-7I
•
I
berekenen van het temperatuurverloop aan de koelzijde van de buizen, luidt:
~ .c (T d) - T m p 1 •z 1 • z+ z
=
2~r.dz.U(T2
•z
- T1 ) •z
Deze vergelijking kan worden herleid tot:
••.•• (5 ) r . I
3)
Gebruikte symbolen.
R ,R 1 2
reaktiesnelheden van resp. methanol en shiftreaktie
X"X 2
conversiegraden van resp. CO en CO
2 reaktiewarmte van resp. methanol en shiftreaktie
AH ,AH 2 1 z
afstand tot top van de reaktor
r
straal van een reaktorbuis
dV
volumeëlement reaktorbuis
é
porositeit katalysatorbed
o ,) ka t
dichtheid katalysator
~m
molenflux per buis
c
soortelijke warmte
p
U T
1 .z
warmteoverdrachtscoëfficient ,T
2.z
temperatuur van resp. koel- en reaktiezijde van een reaktorbuis op plaats z.
Het in 2) besproken stelsel differentiaalvergelijkingen (vgl.2,3,
4,5) dient numeriek te worden opgelost. Een probleem is echter dat de bovengrens van vergelijking
5
de benedengrens is van de
overige vergelijkingen.Dit maakt het onmogelijk om de differentiaalvergelijkingen op te lossen uitgaande van de ingangstemperatuur van de reaktanten.De consequenties hiervan op het stroomschema worden in het volgende hoofdstuk besproken. Men is dus genoodzaakt bij de berekening uit te gaan van de toptemperatuur van de reaktor om de in- en uitgangstemperaturen op te lossen. De differentiaalvergelijkingen worden opgelost met de standaard subroutine DHPCG, welke is gebaseerd op de methode van RungaKutta.De benodigde gegevens voor deze subroutine worden berekend r'
-8-
en/of gegeven in het pro g ramma FCT,waarin ook andere systeem afhankelijke berekeningen worden uit ge voerd (zie bijlage .3 ) Via subroutine OUTP worden de oplossing en overgged ragen ian subroutine FIXBED, waarin met de berekende waarden, de uitgaande stromen worden vast ge legd (zie bijla ge n 2 en 3
).
In het onderdeel van subroutine FIXBED dat op de autotherme reaktor betrekking heeft, , worden bovendi e n het aantal benodigde buizen in de reaktor en de voor DHPCG benodigde gewichtsfaktoren,bepaald. r'
De in de verschillende subroutines gebruikte symbolen worden verklaard in bijlage 4.
r •
,
1
,
1
r '
-9-
I, . l.
Toelichting bU h e t stroomschema ===============================
Voor de produktie van methanol wordt uitgegaan van een lichte nafta of aardgas dat in een reformer wordt omgezet met stoom in het reaktantenmengsel voor de methanolsynthese. Dit deel van het produktieproces kon niet berekend worden omdat een module voor de berekening van de reformer ontbrak. BU de berekeningen is daarom uitgegaan van de produktstroom van de reformer. De zuivering van de methanol door destillatie is eveneens niet in de berekening opgenomen. De geproduceerde hoeveelheid methanol wordt berekend door een volledige scheiding tussen methanol en bUprodukten aan te nemen. Om het resulterende stroomschema (bUlage 1) door te rekenen met behulp van het programma EXEC 2, is het eerst noodzakelUk het programma EXEC 1 te draaien om de rekenvolgorde, de loops en de plaats van iteratieunits in deze loops te bepalen. Kenmerkend voor EXEC 2 is dat het flowsheet stroomafwaards doorgerekend wordt. Aan de hand van de rekenvolgorde zull e n zowel de funkties van de verschillende units als de gebruikte modules besproken worden. 1) Condensor In de produktstroom van de reformer is een aanzienlUke hoeveelheid water aanwezig. Omdat bU de methanolsynthese water de ligging van de shiftreaktie negatief beïnvloed, moet zoveel mogelUk water uit het reaktiemengsel worden verwUderd. Voor de berekening van deze condensor is gebruik gemaakt van de module EVAPO. Er deden zich nogal wat problemen voor bU de berekening van de vapor
ratio maar na enige verbeteringen aan
genoe mde berekening was de module redelUk goed bruikbaar. 2) Compressor De compressor in de voedingsstroom is een drietraps tussenge•
1
koelde compressor die noodzakelUk is om de druk van het reaktiemengsel te verhogen van 35 atm tot de reaktiedruk (340 atm) De gebruikte module is genaamd: COMPR.
3) Menger De menger is alleen in het rekenprogramma een noodzakelUke unit. In de praktUk worden op deze plaats in het proces twee stromen samengevoegd (T-stuk). Naam van de module: MIXER.
,
.
f
1
L
_
-10-
,. l
_
4) Warmtewisselaar, reàktor, informatie-unit De reaktor, warmtewisselaar en de informa tie-unit zullen te ge l~k
worden b e handeld, omdat door de
problemen
r~zen
b~
berekeningsw~ze
van EXEC 2
de berekening van de warmtewisselaar.
De warmtewisselaar moet vol gen s het stroomschema berekend worden vóór de reaktor. Zolang echter de reaktor niet is berekend
z~n
er teveel onbekende grootheden om de warmtewisselaar te kunnen berekenen. Om dit probleem te ondervangen is een informatie module ontworpen (b~lage
2).
De funktie van deze module is het doorgeven van
de gegevens van de stroo m uit de mixer
aan de ingaande stroom
van de reaktor, zodat deze berekend kan worden (uitgaande Tt
op
).
Met de voor de reaktor berekende molfrakties en temperaturen r'
van de in- en uitgaande stromen, kan hiern a de warmtewisselaar berekend worden. De informati e -unit is dus alleen een programma-technische module om een beperking van het EXEC 2 pro gramma te omzeilen. Namen van de modules: FIXBED, HEEX1, INFO.
5) Condensor In deze condensor worden de produkten methanol en water gecondenseerd. Het nog aanwezige CH ,CO,C0
met een kleine hoeveel2 heid methanol en water worden gerecirculeerd.
4
Naam van de module: CONDNS.
6) Splitter Indien er geen gas gespuid zou worden neemt de hoeveelheid methaan steeds toe omdat methaan niet in de reaktor wordt o mgezet maar wel aanwezig is in de voedingsstroom. Het is
noodzakel~k
om juist zoveel gas te spuien dat de hoeveel-
heid methaan in het spuigas
gel~k
is aan de hoeveelheid in de
(loop)voedingsstr6om, zodat de inertspiegel op een aanvaard,
1
bare hoogte
bl~ft.
Het spuigas kan als brandstof worden ge-
bruikt in de reformer. Naam van de op methaan gestuurde module: SPLITTER.
7) Compressor Deze compressor verhoogt de druk in de recyclestroom met 10 atm (drukval in reaktor). Naam van de module: COMPR. 8) Iteratie unit Dit is een programma-technische unit die de totale berekening
r
.
-11-
stuurt en bij voldoende
nauwkeur~gheid
van de iteraties van
warmte- en mas 3abalanzen in de loop, de berekeningen stopt. Naam v a n de module: ITER.
9) Eindmodule Een programma-technische module die de verschillende fasen in het EXEC 2 programma stuurt en de berekeningen beëindigt. Naam van de module: ' LAST.
r'
,, .
Voor de berekening van de reaktor en het stroomschema van een (gedeelte van) methanolsynthesefabriek, is gebruik gemaakt van door Okkes - ·Prins bepaalde (optimale) waarden voor d e grootte en samenstelling van de voedingsstroom,reaktor- en reaktievariabelen etc. Gebleken is dat de
mod~le
FIXBED ,waarin opgenomen de subroutine
voor de berekening van een autotherme reaktor, goed voldoet. De resultaten van de berekening zijn weergegeven in de bijgevoegde stromentabel (bijla ge 5).De overi ge resultaten worden gegeven in de computer-uitvoer. Enkele belangrijke g egevens: met de gebruikte waarden is een methanolproduktie van
934~ton
per jaar (8000 uur) mogelijk (bij volledige
zuivering van het ruwe produkt).Hiervoor is een reaktor nodig met 745 buizen van 6 meter lengte. Hoewel de gevolgde berekeningswijze voor de autotherme reaktor goed voldoet, is er toch een groot bezwaar aan verbonden. Een proces dat op deze wijze wordt uitgevoerd is
nam~lijk
niet regelbaar.Een oplos-
sing hiervoor kan zijn dat met een deel van de voedingsstroom een "cold-shot" in de top van de rea.ktor wordt gegeven. Het lijkt ons nuttig en belangrijk dat dit nader wordt bestudeerd. Litteratuur ===========
1)
Natta G. , Journalof catalysis, ~, 349, (1955)
2)
Vachez F. ,Rev.de l'Inst.Franc. du Petrol, l~, 697, (1963)
3)
van Heerden C. , Ind. Eng. Chem.
4)
Okkes A. Prins W. , Optiaaliseringsberekeningen aan een methanol
,12'
1242, (1953)
produktie-eenheid, Intern Rapport , (1972) 5) . 1 I
,
1
r ,
r' I I
Intern Rapport,Nafta-stoom conversie, T.H Delft
'
,
.
,
,
,.
l.
~I
I~
....
~
'";-
3
~
~
:s
rr-
:J
_.
ç
e
=======
~ijlage
e ...
n-
'"
12_.
IC;; ~ ~
.r.
"0
@3 o
0
f'
e
~
....
-~
~
-f
.
()
...
0 C
lö
,,~ .
....,
....
l~
....
}ol
~
... ~
'~~
~J
I"a
.
3
0
"'"0
~
-.,
I;,
!Tl )(
.."
11"., o o ~ 1-
G
M
:t:
~~
..,
I~
"
1"-' I
_. 3
...,
~
e /
la-.
a
;:)
0
(\
..t:
.
t-o
10'
,~
;I)
rn
)(
e ...
Iw"
1:)
e~ ...
..,
/
0
:l
e ...
n 0
1-, ~;. \P
0
!ft
o
<:
f
'
I
.
l
.
bijlage 2
======= = =
subroutine OUTPUT
module INFO
START
START
overdracht gegevens
ja
nee bereken molstroom per comp. per buis
END
nee
ja
schrijf: buislengte conversiegraden,temperaturen
END
r .
stroomgegevens van inen uitgaande stroom aan elkaar gelijk stellen
I
'
subroutine FCT
modul e FIXBED
START
START
overdracht gegevens
overdracht gegevens
berekening fugaciteit e n
schrijf stroomtabel
berekening constant en U:ït r.s. vgl.
nee
SP(15)<5
berekening aantal buizen, molstroom per buis,samenstellin g gasmengsel
geef differe ntiaal vergelijkingen
bepaal beginvoorwaarden voor in tegra tie
END
f '
.. r'
ja
call DHPCG
omzetting symbolen in EXEC variabelen
schrijf strocimtabel
END
berekening van adiabatische en isoth e rme reaktoren
,.
I '
Verklaring van de gebruikte symbolen (met uitzondering van de EXEC-symbolen), onderverdee ld per subroutine: FIXBED (alleen de voor de autotherme reaktor benodigde symbolen worden verklaard): nummering van de komponenten: i=1
CO
2
CO
r .
I
3 4 5
2 H 2 H 0 (G) 2 CH
6
CH 0H (G)
8
CH 0H (L)
4
3 7 H2 0 (L) 3
SAM( 1)
r •
PR(I)
aantal molen van komp. i in de voedingsstroom van een buis (mol/s) percentage aan It .. "
OUTR( 1)
aantal molen van komp. i in de uitgangsstroom van een buis
AS, BS
konstante uit de vergelUking voor de
R
gaskonstante
NBUIS
aantal buizen in reaktor
TTOP
toptemperatuur reaktor
YMAX(I)
waarden gebruikt bU de bepaling der gewichtsfaktoren
SOMY
hulpvariabele
NDIM
aantal differentiaal vergelUkingen
PAR
hulpparameter
"
.1
(kcal/mol. oK)
(oK)
X1
konversie CO tot methanol
X2
konversie CO
PRMT(1)
evenwichts~~RJ/&dte
tot CO 2 ondergrens integratie interval
(2)
bovengrens integratie interval (is lengte buis in meters)
(3)
stapgrootte bU integratie
(4 )
absolute nauwkeurigheid bU integratie
SOMO
hulpvariabele
Y(I)
afhankelUke variabelen
FCT
r'
X
onafhankelUk variabele van stelsel diff. vergelUkingen
CP
soortelUke warmte (cal/oK.mol)
L .
FCO etc.
beginfug aciteit van CO en andere komponenten (atm.)
ACO etc. PTOT
fugaciteit van CO e.a. komp. op plaats van berekening (atm) druk in reaktor
T1
temperatuur
T2
temperatuur in buis
AA,B
konstanten in diff.
S,T,U,V,W,
termen uit de
(atm)
koelz~de
reaktorbuis
( OK) ( OK)
vergel~kingen
reaktiesnelheidsvergel~king
TEL R1,R2 DERy(r)
r •
r'
voor de metha-
nolreaktie reaktiesnelheden voor resp. methanol- en shiftreaktie numerieke waarde diff.
vergel~king
i
Sl~GC"TAAF.l
---------
ne
VOEDING 0.0 0.161\0C
0.0 0.05';"11
r.o
0.0 0.0
0.0 0.0
0.0 0.0
0.0 0.0
0.0 0.0
2.0000 0.0
0.0 0.0
0.0 0.0
0.0 0.0
C.O 0.0
2. 0\)0 0 0.0
0.0 0.0
0.0 0.0
0.0 0.0
30S. ,OlOO 0. 0
0.0 0.0
0.0 1.00000
0.0 0.0
0.0 0.0
0.0 0.0
339.Q998 0.J t. 591
0.0 J.O
1.0000 1.00C'J0
0.0 0.0
0.0 0.0
0.0 0.0
0.0 (I.00C09
0.0 0.0
0.0 0.0
0.0 0.0
0.0 0.0
0.0 0.0
0.0 0.0
0.0 0.0
0.0
4 ti7 .T7h6
J5. ~:~~
1l71.J Q 9Q
O . ~ T'JO O
O.lI,- Jl'
O. l' l~O'l
402.441~
35.n ::0
l3q. qq 9~
0.6~722
0.( !)t~1
0 • .14 117
40~.4412
340. 0; ~ ~
522.P577
n.~~122
o .n 0,i7
0. ~ 4117
~.O
CO~D1UIT
C.O
2
0.20306
0.0 0.061~d
MIX IN
3
0.0 0.20306
0.0 0 .061~S
0.0
CONDI0U
c.:)
O.C'
a.erNO 0.0 0.00266 MIX
0.0 0.0 'l.0 0.012è6
32044.5625
0.0
0.0
0.0
85.3354
35.0()0
~.r6962
O.q1~76
0.0 5026.7227 0.)3731 C.7:3g6
340.00JO 0.00:~4
332.6d46 0.21480
0.0 5026.1227 0.03737 0.7:3i6
340.0(1 '0
577.~~65
0.00r~4
0.21480
0.00009
0.0 0.0
0.0 4822.5586 0.03434 0.68691
330.0 C)" O.G0529
647.9602
648.3999 1.02112
0.0 0.0
0.0 0.0
0.0 0.0
0.0 0.0
0.0
0.0 0.03434
482~.55S6
330.0( ' 0 0.005~9
540.50 34 0.22389
64~.399g
O.6~691
0.02112
0.0 0.0
0.0 0.0
0.0 0.0
0.0 0.0
0.0
0.0 0.00010
0.0 0.0
0.0 0.0
0.0 0.0
0.0 0.0
0.0
HU
0.0 0.04315
'l.0
I1
REAK 9IN
• • • • • • • • • •
0.0 0.04315 0.0 0.02841 CCN010IN 0.0 0.02e47
10
0.0 0.02923 0.0 0.0
11.
,
,
~
o ~-·~ o
314. 0 000
0.00ClG
0.22991
0.0 0.0
126.3735 0.0
330.00)0 0.0
314.0000 0. 0
0.0 0.0
1.0000 0.19766
0.0 0.80214
0.0 0.0
0.0 0.0
0.0
72. 0664 0.70 539
330.C'C J O
314.0('00 0.22991
0.0 'l.00010
0.0 0.0
0.0 0.0
0.0 0.0
0.0 0.0
0.0
314.00 00 0.22991
0.0 0.00010
0.0 0.0
0.0 0.0
0.0 0.0
0.0 C.O
0.0
SPUI 0.0 0.02923
0.0 0.03526
CCMP13 IN 0.0 0.02923
0.0 4624.1133 0.03526 0.7 0 539
330. 0C~(,
G.O 4624.1133 0.03526 0.7 C539
340.0~10
0 .0 0010
31ó.6057 0.22991
0.0 0 .00010
0.0 0.0
0.0 0.0
0.0 0.0
0.0 0.0
0.0
0.0 4624. 2852 0.035Z5 0.7 C541
34 0 .1)(;0
316.~05 7
O.OQct'J
0.22991
0.0 0.00010
0.0 0.0
0.0 0.0
0.0 0.0
0.0 0.0
0.0
1.00';" 0.0
235.')000 0. 0
G.O
0.0 1.00000
0.0 0.0
0.0 0.0
0.0 0.0
'l.0
1. oe,' ')
308. G0 00
').0 (.0
0.0 1.00000
0.0 0.0
0.0 0.0
Cl.O
0.0
'l.0 C.)
0.0 1.0000')
0.0 0.0
0.0 0.0
0.0 0.0
'J.O
0.0 0.C2923
IS"
0.0 0.02923 IN
0.0 0.0
1C992.5078 0.0
0.0 0.0
10992. 5'1 7 él
KOElU IT
C.O 0.0 0.0 0.0
Co
.
O.OOClO
0.00C10
HULPSTR
q
•
30.
I TER IN
)4
IU
,
0.0 4696.1797 0.03526 0.7~539
P~OOST
Q. KOHIN
,
0.~2389
SPLITIN
0.0 0.0
,
,I
WIoSTRH
KOel
.
C.Cl
0.0 0.0
".')
0.0
1.0(': 0
:: .0
o.c
285.1 n OO
o.e
0.0
0.0
I
0 ·0 ..
'.
-{] ,~ .
1.0 ,
0:" 0
0
0 0 0 0
\
0 0 0 0 0 -0
I
0
\0. ..'
,,'
"'0
,.