Télécharger chi1.notice

Retour à la liste

Afficher cette notice en

Numérotation des lignes :
   1 : $$$$ CHI1     NOTICE  CHAT      11/09/12    21:15:24     7124           
   2 :                                              DATE     11/09/12
   3 :   
   4 :   Operateur CHI1                                Voir aussi :
   5 :    --------------  
   6 :        TAB2 = CHI1 TAB1 'COMP' VAL1 < 'LOGK' VAL2 > <'ENTH' VAL3 >    ;
   7 : 
   8 : 
   9 :         Objet
  10 :        -------
  11 :        Le but est de calculer la speciation d'une eau, en tout point
  12 :        d'un domaine a partir de la donnee des concentrations analytiques
  13 :        de chaque composant chimique du systeme. Le calcul se fait en
  14 :        deux temps a l'aide des operateurs CHI1 et CHI2.
  15 :        CHI1 rassemble toutes les donnees relatives a un systeme chimique,
  16 :        et CHI2 effectue la speciation.
  17 :         La terminologie est celle de Mineql.
  18 : 
  19 :           Toutes les concentrations sont donnees en moles par litre.
  20 : 
  21 :        Commentaires
  22 :        ------------
  23 :        TAB1 est une TABLE indicee par les mots:
  24 :             'IDEN'  <'CHXMX'>  <'BDD'>  <'CLIM'>  <'NVCOMP'>
  25 :             <'NVESP'> <'ECHANGE'>  <'TEMPERATURE'>
  26 : 
  27 :        TAB1.IDEN est un objet de type LISTENTI contenant les
  28 :              identifiants (dans la base de donnees),des composants chimi
  29 :              ques a utiliser.
  30 : 
  31 :        TAB1.CHXMX est un objet de type LISTENTI contenant les identifi
  32 :              ants des mineraux a retenir. A defaut on conserve tous les
  33 :              mineraux dont les composants sont dans TAB1.IDEN.
  34 : 
  35 :        TAB1.BDD  contient un mot servant a preciser le format de la base
  36 :              de donnees. 'STRASBG' ou 'MINEQL' .
  37 :              'MINEQL' correspond a la base de donnees standard de Mineql.
  38 :              'STRASBG' correspond a la base de donnees issue de Kindis.
  39 :              Les formats sont decrits dans le rapport DMT/94/597.
  40 :              L'option par defaut est 'MINEQL'.
  41 : 
  42 :        TAB1.CLIM  est une TABLE servant a definir les contraintes chimi
  43 :              ques. Cette TABLE est indicee par des mots tous facultatifs.
  44 :              <'TYP3'> <'COMP3'> <'TYP4'>  <'TYP5'> <'TYP6'>
  45 : 
  46 :              TAB1.CLIM.TYP3 est un objet de type LISTENTI contenant les
  47 :                   identifiants des especes dont on veut imposer l'activi
  48 :                   te.
  49 : 
  50 :              TAB1.CLIM.COMP3  est un objet de type LISTENTI contenant
  51 :                   pour chacune des especes de TAB1.CLIM.TYP3 l'identifi
  52 :                   ant du composant immobile . Si TAB1.CLIM.TYP3 ne
  53 :                   contient que des especes simples cette donnee est 
  54 :                   inutile.
  55 : 
  56 :              TAB1.CLIM.TYP4 est un objet de type LISTENTI contenant les
  57 :                   identifiants des especes precipitees.
  58 : 
  59 :              TAB1.CLIM.TYP5 est un objet de type LISTENTI contenant les
  60 :                   identifiants des especes en solution, pouvant etre
  61 :                   precipites.
  62 : 
  63 :              TAB1.CLIM.TYP6 est un objet de type LISTENTI contenant les
  64 :                   identifiants des especes non prises en compte.
  65 : 
  66 :        TAB1.NVCOMP  est une TABLE permettant de rajouter des composants
  67 :             (ne figurant pas dans la base de donnees). Pour n composants
  68 :              , cette TABLE sera indicee par des nombres de 1 a n.
  69 :              Pour le i ieme composant a rajouter TAB1.NVCOMP.i sera une
  70 :              TABLE indicee par les mots:  'IDEN' 'NOM' 'CHARGE'
  71 :              TAB1.NVCOMP.i.IDEN est un entier identifiant du nouveau
  72 :                        composant.
  73 :              TAB1.NVCOMP.i.NOM est un mot nom de ce composant
  74 :              TAB1.NVCOMP.i.CHARGE  est un entier charge de l'espece
  75 :                        simple associee.
  76 : 
  77 :        TAB1.NVESP est une TABLE permettant de rajouter ou de modifier 
  78 :              des especes.Pour n especes, cette TABLE sera indicee par
  79 :              des nombres, de 1 a n.
  80 :              Pour la i ieme espece a rajouter TAB1.NVESP.i sera une
  81 :              TABLE indicee par les mots:
  82 :                   'IDEN' 'LOGK' <'ITYP'> <'COMP'> <'STOECH'>
  83 : 
  84 :              TAB1.NVESP.i.IDEN est un entier identifiant de l'espece
  85 : 
  86 :              TAB1.NVESP.i.LOGK est un reel logk de l'espece
  87 : 
  88 :              TAB1.NVESP.i.ITYP entier type de l'espece
  89 :                                2   complexe en solution
  90 :                                3   activite fixee
  91 :                                4   mineraux precipites
  92 :                                5   mineraux dissous
  93 :                                6   non pris en compte dans le calcul
  94 : 
  95 :              TAB1.NVESP.i.COMP objet LISTENTI contenant les identifiants
  96 :                               des composants de l'espece. Le nombre de
  97 :                               ces identifiants doit etre inferieur a 4
  98 :                               pour une base de donnee de type MINEQL
  99 :                               et inferieur a 8 pour une base de donnee
 100 :                               de type STRASBG.
 101 : 
 102 :              TAB1.NVESP.i.STOECH    objet LISTREEL coefficient
 103 :                               stoechiometrique correspondant a chacun 
 104 :                               de ces composants.
 105 :              TAB1.NVESP.i.NOMESPECE  mot. Nom de cette nouvelle espece
 106 :                                chimique. ( au plus 32 caracteres.La
 107 :                                valeur par defaut est contituee de 32
 108 :                                blancs)
 109 :        
 110 :        TAB1.NVSOSO est une TABLE permettant d'ajouter des solutions
 111 :              solides (qui ne figurent pas encore dans les bases de 
 112 :              donnees). Pour n solutions solides, cette TABLE sera indicee 
 113 :              par des nombres, de 1 a n.
 114 :              Pour la i ieme solution solide a ajouter TAB1.NVSOSO.i sera 
 115 :              une TABLE indicee par les mots:
 116 :                       'IDEN' <'ITYP'> 'SOLID' <'FRACTIO'>
 117 : 
 118 :              TAB1.NVSOSO.i.IDEN est un entier identifiant de la solution 
 119 :                               solide
 120 : 
 121 :              TAB1.NVSOSO.i.ITYP entier type de la solution solide
 122 :                                3  activite fixee
 123 :                                4  solutions solides precipitees
 124 :                                5  solutions solides dissoutes
 125 :                                6  non pris en compte dans le calcul
 126 :                          Pour les types 3 et 4, il faut obligatoirement
 127 :                          donner les fractions molaires des poles des
 128 :                          solutions solides; pour les types 5 et 6, ce 
 129 :                          n'est pas obligatoire.   
 130 : 
 131 :               TAB.1.NVSOSO.i.SOLID objet LISTENTI contenant les 
 132 :                          identifiants des poles mineraux purs de la
 133 :                          solution solide. Le nombre de ces poles doit
 134 :                          etre inferieur a 36. Ces poles sont mis 
 135 :                          automatiquement en type 6 (ils servent au calcul,
 136 :                          mais n'ont pas d'existance physique).
 137 :                
 138 :               TAB1.NVSOSO.i.FRACTIO objet LISTREEL contenant les fractions
 139 :                          molaires correspondant a chacun des poles.
 140 :                          (Si la solution solide est mise en type 3 ou 4, 
 141 :                          l'operateur chi1 a besoin des fractions molaires 
 142 :                          pour calculer les coefficients stoechiometriques
 143 :                          ainsi que le logK de la solution solide. 
 144 :                          Si la solution solide est mise en type 5 ou 6, 
 145 :                          l'operateur chi2 calculera lui meme les fractions 
 146 :                          molaires et le reste).
 147 : 
 148 :        TAB1.ECHANGE objet de type LISTENTI contenant les identifiants 
 149 :              des sites de surface par echange ionique.
 150 :              L'espece simple associee sera mise automatiquement en 
 151 :              type 6 (sert au calcul mais n'a pas d'existence physique) .
 152 :              Soient les deux epeces echangeables c1b et c2b telles que
 153 :              n1 * c1b + n2 * c2b = ns ( nombre total de sites ) liees
 154 :              par le coefficient de selectivite 
 155 :                                 k12=(c1b/c1)**n2*(c2/c2b)**n1 .
 156 :              On definira deux especes c1b et c2b dont les constantes
 157 :              de formation seront  k1=c1b/c1*(cidx)**n1 et
 158 :              k2=c2b/c2*(cidx)**n2  (avec cidx concentration de l'espece
 159 :              simple associee au site) telles que k12=(k1)**n2/(k2)**n1.
 160 :              Le coefficient d'activite des especes de surface est prise
 161 :              egale a 1
 162 : 
 163 :        TAB1.TEMPERATURE mot 'OUI'  'NON' ou entier (1 ou 2).
 164 :               'NON' on ne tient pas compte de la temperature. C'est 
 165 :               l'option par defaut.
 166 :             - Cas de la base STRASBG.
 167 :               Si 'OUI' on prendra en compte les effets thermiques sur 
 168 :               le logk, par interpolation de donnees tabulees.
 169 :             - Cas de la base MINEQL.
 170 :                1 ou 'OUI' on utilise la premiere approximation d'Ulich
 171 :                     K(T)=K0+f(H(T)-H(T0))
 172 : 
 173 :                2 on utilise la deuxieme approximation d'Ulich 
 174 :                     K(T)=K0+f((H(T)-H(T0)),(Cp(T)-Cp(T0)))
 175 : 
 176 :        COMP est un mot-cle
 177 : 
 178 :        VAL1  est l'adresse du fichier de composants de la
 179 :              base de donnees. Cette adresse peut etre, soit une chaine 
 180 :              de caracteres (le nom du fichier), soit un entier (l'unite
 181 :              de lecture). Dans le cas ou l'on donne le nom du fichier, 
 182 :              le code lui affecte la premiere unite de lecture disponible
 183 :              a partir de 80. Cette unite sera liberee a la sortie de
 184 :              l'operateur.
 185 : 
 186 :        LOGK  est un mot-cle
 187 : 
 188 :        VAL2  est l'adresse du fichier de logk de la base de donnees.
 189 :              Cette adresse peut etre,soit une chaine de caracteres
 190 :              (le nom du fichier), soit un entier (l'unite de lecture).
 191 :              Dans le cas ou l'on donne le nom du fichier, le code lui
 192 :              affecte la deuxieme unite de lecture disponible a partir de
 193 :              80. Cette unite sera liberee a la sortie de l'operateur.
 194 :              Dans le cas ou les composants et les logk sont sur le meme
 195 :              fichier il est inutile de preciser 'LOGK' VAL2, toutefois
 196 :              si ces donnees figurent,le code reconnaitra que les
 197 :              fichiers ont le meme nom et lira en sequence.
 198 : 
 199 :        ENTH  est un mot-cle
 200 : 
 201 :        VAL3  est l'adresse du fichier de la base de donnees,
 202 :              permettant de prendre en compte les effets thermiques.
 203 :              Cette adresse peut etre,soit une chaine de caracteres 
 204 :              (le nom du fichier),soit un entier (l'unite de lecture).
 205 :              Dans le cas ou l'on donne le nom du fichier, le code 
 206 :              lui affecte la troisieme unite de lecture disponible a
 207 :              partir de 80. Cette unite sera liberee a la sortie de 
 208 :              l'operateur.
 209 :              Dans le cas ou ces donnees et les logk sont sur le meme
 210 :              fichier il est inutile de preciser 'ENTH' VAL3, toutefois
 211 :              si VAL2 et VAL3 sont identiques,le code lira en sequence.
 212 : 
 213 : 
 214 : 
 215 : 
 216 :        TAB2 est une TABLE contenant le systeme chimique. Cette TABLE a
 217 :             pour sous type CHIMI1, et elle est indicee par les mots:
 218 :                     'DESCHI' 'IDEN' <'REDOX'>  <'TEMPE'>
 219 : 
 220 :        TAB2.DESCHI est une TABLE contenant la description du systeme
 221 :              chimique. Elle est indicee par les mots:
 222 :             'IDX' 'NOM' 'CHARGE' 'IDY' 'NN' 'MATRICEA' 'LOGK' 'IDSURF'
 223 :             <'IDZ' 'MATRICEF' 'IDP'>
 224 :  
 225 :              TAB2.DESCHI.IDX  objet LISTENTI contenant les identifiants
 226 :                    des composants retenus
 227 :       
 228 :              TAB2.DESCHI.NOM  objet LISTMOTS contenant les noms de ces
 229 :                    composants (sur 32 caracteres).
 230 :       
 231 :              TAB2.DESCHI.CHARGE objet LISTENTI contenant la charge de
 232 :                    ces composants.
 233 :       
 234 :              TAB2.DESCHI.IDY objet LISTENTI contenant les identifiants 
 235 :                    des especes retenues.
 236 : 
 237 :              TAB2.DESCHI.NOMESPECE  objet LISTMOTS contenant les noms
 238 :                    de ces especes (sur 32 caracteres).
 239 :       
 240 :              TAB2.DESCHI.NN objet LISTENTI de 6 valeurs contenant le
 241 :                    nombre d'especes de chaque type c'est a dire:
 242 :                     nombre d'especes simples en solution (type 1)
 243 :                     nombre de complexes en solution (type 2)
 244 :                     nombre d'especes a activite fixee (type 3)
 245 :                     nombre de mineraux precipites (type 4)
 246 :                     nombre de mineraux precipitables en solution(type 5)
 247 :                     nombre d'especes non prises en compte dans le calcul
 248 :                     (type 6)
 249 :       
 250 :              TAB2.DESCHI.MATRICEA objet LISTREEL matrice contenant les 
 251 :                    coefficients stoechiometriques de chaque espece.
 252 :       
 253 :              TAB2.DESCHI.LOGK  objet LISTREEL valeur pour chaque espece
 254 :                    de la constante de formation logk.
 255 :       
 256 :              TAB2.DESCHI.IDSURF objet LISTENTI  nature des especes de 
 257 :                    type 2.
 258 :                    0 en solution ,1 echange ionique (fixee) ,2 adsorbee
 259 : 
 260 :              TAB2.DESCHI.IDZ objet LISTENTI contenant les identifiants 
 261 :                    des solutions solides retenues.
 262 : 
 263 :              TAB2.DESCHI.IDP objet LISTENTI contenant les identifiants 
 264 :                    des poles de solutions solides retenus.
 265 : 
 266 :              TAB2.DESCHI.MATRICEF objet LISTREEL matricie contenant les 
 267 :                    fractions molaires de chaque solution solide
 268 : 
 269 :        TAB2.IDEN est une TABLE indicee par les mots:
 270 :              'COMP' 'SOLU' 'NOMSOLU' <'PRECI' 'NOMPRECI'>
 271 :              <'SURF' 'NOMSURF'>  <'TYP3' 'NOMTYP3'> <IMP3>
 272 :              <'TYP6' 'NOMTYP6'>  <'PARF'> <'REAC'>  <'IMMO'>
 273 :              <'SOSO' 'NOMSOSO' 'POLE' 'NOMPOLE'>
 274 : 
 275 :              TAB2.IDEN.COMP objet de type LISTENTI entiers servant a
 276 :                    decrire la nature des composants.
 277 :                    1 conservatifs, 2 reactifs, 3 immobiles.
 278 :       
 279 :              TAB2.IDEN.SOLU objet LISTENTI identifiants des especes en
 280 :                    solution
 281 :       
 282 :              TAB2.IDEN.NOMSOLU objet LISTMOTS noms attribues par le code
 283 :                    a ces especes en solution. Ces noms ont 4 caracteres
 284 :                    dont le premier est W suivi eventuellement de 0 ou 00
 285 :                    et du numero d'ordre dans la liste TAB2.DESCHI.IDY de
 286 :                    l'identifiant concerne.
 287 :       
 288 :              TAB2.IDEN.PRECI objet LISTENTI identifiants des especes 
 289 :                    precipites potentiels.
 290 :       
 291 :              TAB2.IDEN.NOMPRECI objet LISTMOTS noms attribues par le 
 292 :                    code a ces especes. Ces noms ont 4 caracteres dont le
 293 :                    premier est W suivi eventuellement de 0 ou 00 et du
 294 :                    numero d'ordre dans la liste  TAB2.DESCHI.IDY  de
 295 :                    l'identifiant concerne.
 296 :       
 297 :              TAB2.IDEN.SURF objet LISTENTI identifiants des especes de 
 298 :                    surface.
 299 :       
 300 :              TAB2.IDEN.NOMSURF objet LISTMOTS noms attribues par le code 
 301 :                    a ces especes.
 302 :       
 303 :              TAB2.IDEN.TYP3 objet LISTENTI identifiants des especes dont
 304 :                    l activite est imposee(meme liste que TAB1.CLIM.TYP3)
 305 :       
 306 :              TAB2.IDEN.NOMTYP3 objet LISTMOTS noms attribues par le code
 307 :                    a ces especes.
 308 :       
 309 :              TAB2.IDEN.IMP3 objet LISTENTI identifiants des especes de
 310 :                    type 3. (activite imposee, redox).
 311 :       
 312 :              TAB2.IDEN.TYP6  objet LISTENTI identifiants des especes non
 313 :                    prises en compte dans le calcul.
 314 :       
 315 :              TAB2.IDEN.NOMTYP6  objet LISTMOTS noms attribues par le code
 316 :                    a ces especes.
 317 :       
 318 :              TAB2.IDEN.PARF objet LISTENTI identifiants des composants 
 319 :                    non reactifs
 320 :       
 321 :              TAB2.IDEN.REAC objet LISTENTI identifiants des composants
 322 :                    reactifs
 323 :       
 324 :              TAB2.IDEN.IMMO objet LISTENTI identifiants des composants
 325 :                    non transportes( par exemple:composant de surface,
 326 :                    composant associe a une espece a activite imposee )
 327 : 
 328 :              TAB2.IDEN.SOSO objet LISTENTI identifiants des solutions 
 329 :                    solides.
 330 : 
 331 :              TAB2.IDEN.NOMSOSO objet LISTMOTS noms attribues par le code a
 332 :                    ces solutions solides.
 333 : 
 334 :              TAB2.IDEN.POLE objet LISTENTI identifiant des poles de 
 335 :                    solutions solides.
 336 : 
 337 :              TAB2.IDEN.NOMPOLE objet LISTMOTS noms attribues par le code a
 338 :                    ces poles de solutions solides.
 339 : 
 340 :        TAB2.REDOX objet TABLE contenant les donnees relatives aux redox.
 341 :              Cette TABLE est indicee par les mots: 'I1' 'I2' 'A1' 'A2'
 342 : 
 343 :              TAB2.REDOX.I1  objet LISTENTI identifiant du composant
 344 :                    mobile
 345 :       
 346 :              TAB2.REDOX.I2 objet LISTENTI identifiant du composant lie 
 347 :                    au composant mobile.
 348 :       
 349 :              TAB2.REDOX.A1 objet LISTREEL coefficient stoechiometrique
 350 :                    de l'electron e-
 351 :       
 352 :              TAB2.REDOX.A2 objet LISTREEL coefficient stoechiometrique
 353 :                    du proton H+
 354 : 
 355 :        TAB2.TEMPE est une TABLE contenant les donnees relatives a la
 356 :                temperature pour les especes selectionnees. Elle est
 357 :                indicee par les mots:
 358 :                'DELTAH' 'DELCP' 'APPROX'  pour la base MINEQL
 359 :                'COEF' 'LOGK' 'TMIN' 'TMAX' 'NVT' 'NUMT' pour la base
 360 :                  STRASBG
 361 : 
 362 :              TAB2.TEMPE.DELTAH objet LISTREEL. Variation d'enthalpie
 363 : 
 364 :              TAB2.TEMPE.DELCP objet LISTREEL. Variation de capacite
 365 :                    calorifique
 366 : 
 367 :              TAB2.TEMPE.APPROX entier (1 ou 2) type d'approximation.
 368 : 
 369 :              TAB2.TEMPE.COEF objet LISTREEL.Coefficients d'interpolation.
 370 :                    Pour chaque espece retenue,il y a 4 valeurs pour 
 371 :                    chacune des 8 temperatures de reference.
 372 : 
 373 :              TAB2.TEMPE.TMINA objet LISTREEL. Pour chaque espece retenue,
 374 :                    temperature auquelles on donne le LogK.( cette liste
 375 :                    est dimensionnee au (nombre d'especes)*(max de valeurs)
 376 : 
 377 :              TAB2.TEMPE.LOGK objet LISTREEL. Valeurs des Logk aux
 378 :                    temperatures de la liste TMIMA.
 379 : 
 380 :              TAB2.TEMPE.NVT objet LISTENTI. Pour chaque espece retenue,
 381 :                     nombre de temperatures disponibles.
 382 : 
 383 :              TAB2.TEMPE.NUMT objet LISTENTI. Pour chaque espece retenue,
 384 :                     identifiant de l'espece.
 385 : 
 386 : 
 387 : 
 388 : 
 389 :  
 390 : 
 391 :  
 392 : 
 393 :  
 394 :  
 395 :  
 396 :  
 397 :  
 398 :  
 399 :  

© Cast3M 2003 - Tous droits réservés.
Mentions légales