Afficher cette notice en
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