Afficher cette notice en
$$$$ CHI1 NOTICE CHAT 11/09/12 21:15:24 7124 DATE 11/09/12 Operateur CHI1 Voir aussi : -------------- TAB2 = CHI1 TAB1 'COMP' VAL1 < 'LOGK' VAL2 > <'ENTH' VAL3 > ; Objet ------- Le but est de calculer la speciation d'une eau, en tout point d'un domaine a partir de la donnee des concentrations analytiques de chaque composant chimique du systeme. Le calcul se fait en deux temps a l'aide des operateurs CHI1 et CHI2. CHI1 rassemble toutes les donnees relatives a un systeme chimique, et CHI2 effectue la speciation. La terminologie est celle de Mineql. Toutes les concentrations sont donnees en moles par litre. Commentaires ------------ TAB1 est une TABLE indicee par les mots: 'IDEN' <'CHXMX'> <'BDD'> <'CLIM'> <'NVCOMP'> <'NVESP'> <'ECHANGE'> <'TEMPERATURE'> TAB1.IDEN est un objet de type LISTENTI contenant les identifiants (dans la base de donnees),des composants chimi ques a utiliser. TAB1.CHXMX est un objet de type LISTENTI contenant les identifi ants des mineraux a retenir. A defaut on conserve tous les mineraux dont les composants sont dans TAB1.IDEN. TAB1.BDD contient un mot servant a preciser le format de la base de donnees. 'STRASBG' ou 'MINEQL' . 'MINEQL' correspond a la base de donnees standard de Mineql. 'STRASBG' correspond a la base de donnees issue de Kindis. Les formats sont decrits dans le rapport DMT/94/597. L'option par defaut est 'MINEQL'. TAB1.CLIM est une TABLE servant a definir les contraintes chimi ques. Cette TABLE est indicee par des mots tous facultatifs. <'TYP3'> <'COMP3'> <'TYP4'> <'TYP5'> <'TYP6'> TAB1.CLIM.TYP3 est un objet de type LISTENTI contenant les identifiants des especes dont on veut imposer l'activi te. TAB1.CLIM.COMP3 est un objet de type LISTENTI contenant pour chacune des especes de TAB1.CLIM.TYP3 l'identifi ant du composant immobile . Si TAB1.CLIM.TYP3 ne contient que des especes simples cette donnee est inutile. TAB1.CLIM.TYP4 est un objet de type LISTENTI contenant les identifiants des especes precipitees. TAB1.CLIM.TYP5 est un objet de type LISTENTI contenant les identifiants des especes en solution, pouvant etre precipites. TAB1.CLIM.TYP6 est un objet de type LISTENTI contenant les identifiants des especes non prises en compte. TAB1.NVCOMP est une TABLE permettant de rajouter des composants (ne figurant pas dans la base de donnees). Pour n composants , cette TABLE sera indicee par des nombres de 1 a n. Pour le i ieme composant a rajouter TAB1.NVCOMP.i sera une TABLE indicee par les mots: 'IDEN' 'NOM' 'CHARGE' TAB1.NVCOMP.i.IDEN est un entier identifiant du nouveau composant. TAB1.NVCOMP.i.NOM est un mot nom de ce composant TAB1.NVCOMP.i.CHARGE est un entier charge de l'espece simple associee. TAB1.NVESP est une TABLE permettant de rajouter ou de modifier des especes.Pour n especes, cette TABLE sera indicee par des nombres, de 1 a n. Pour la i ieme espece a rajouter TAB1.NVESP.i sera une TABLE indicee par les mots: 'IDEN' 'LOGK' <'ITYP'> <'COMP'> <'STOECH'> TAB1.NVESP.i.IDEN est un entier identifiant de l'espece TAB1.NVESP.i.LOGK est un reel logk de l'espece TAB1.NVESP.i.ITYP entier type de l'espece 2 complexe en solution 3 activite fixee 4 mineraux precipites 5 mineraux dissous 6 non pris en compte dans le calcul TAB1.NVESP.i.COMP objet LISTENTI contenant les identifiants des composants de l'espece. Le nombre de ces identifiants doit etre inferieur a 4 pour une base de donnee de type MINEQL et inferieur a 8 pour une base de donnee de type STRASBG. TAB1.NVESP.i.STOECH objet LISTREEL coefficient stoechiometrique correspondant a chacun de ces composants. TAB1.NVESP.i.NOMESPECE mot. Nom de cette nouvelle espece chimique. ( au plus 32 caracteres.La valeur par defaut est contituee de 32 blancs) TAB1.NVSOSO est une TABLE permettant d'ajouter des solutions solides (qui ne figurent pas encore dans les bases de donnees). Pour n solutions solides, cette TABLE sera indicee par des nombres, de 1 a n. Pour la i ieme solution solide a ajouter TAB1.NVSOSO.i sera une TABLE indicee par les mots: 'IDEN' <'ITYP'> 'SOLID' <'FRACTIO'> TAB1.NVSOSO.i.IDEN est un entier identifiant de la solution solide TAB1.NVSOSO.i.ITYP entier type de la solution solide 3 activite fixee 4 solutions solides precipitees 5 solutions solides dissoutes 6 non pris en compte dans le calcul Pour les types 3 et 4, il faut obligatoirement donner les fractions molaires des poles des solutions solides; pour les types 5 et 6, ce n'est pas obligatoire. TAB.1.NVSOSO.i.SOLID objet LISTENTI contenant les identifiants des poles mineraux purs de la solution solide. Le nombre de ces poles doit etre inferieur a 36. Ces poles sont mis automatiquement en type 6 (ils servent au calcul, mais n'ont pas d'existance physique). TAB1.NVSOSO.i.FRACTIO objet LISTREEL contenant les fractions molaires correspondant a chacun des poles. (Si la solution solide est mise en type 3 ou 4, l'operateur chi1 a besoin des fractions molaires pour calculer les coefficients stoechiometriques ainsi que le logK de la solution solide. Si la solution solide est mise en type 5 ou 6, l'operateur chi2 calculera lui meme les fractions molaires et le reste). TAB1.ECHANGE objet de type LISTENTI contenant les identifiants des sites de surface par echange ionique. L'espece simple associee sera mise automatiquement en type 6 (sert au calcul mais n'a pas d'existence physique) . Soient les deux epeces echangeables c1b et c2b telles que n1 * c1b + n2 * c2b = ns ( nombre total de sites ) liees par le coefficient de selectivite k12=(c1b/c1)**n2*(c2/c2b)**n1 . On definira deux especes c1b et c2b dont les constantes de formation seront k1=c1b/c1*(cidx)**n1 et k2=c2b/c2*(cidx)**n2 (avec cidx concentration de l'espece simple associee au site) telles que k12=(k1)**n2/(k2)**n1. Le coefficient d'activite des especes de surface est prise egale a 1 TAB1.TEMPERATURE mot 'OUI' 'NON' ou entier (1 ou 2). 'NON' on ne tient pas compte de la temperature. C'est l'option par defaut. - Cas de la base STRASBG. Si 'OUI' on prendra en compte les effets thermiques sur le logk, par interpolation de donnees tabulees. - Cas de la base MINEQL. 1 ou 'OUI' on utilise la premiere approximation d'Ulich K(T)=K0+f(H(T)-H(T0)) 2 on utilise la deuxieme approximation d'Ulich K(T)=K0+f((H(T)-H(T0)),(Cp(T)-Cp(T0))) COMP est un mot-cle VAL1 est l'adresse du fichier de composants de la base de donnees. Cette adresse peut etre, soit une chaine de caracteres (le nom du fichier), soit un entier (l'unite de lecture). Dans le cas ou l'on donne le nom du fichier, le code lui affecte la premiere unite de lecture disponible a partir de 80. Cette unite sera liberee a la sortie de l'operateur. LOGK est un mot-cle VAL2 est l'adresse du fichier de logk de la base de donnees. Cette adresse peut etre,soit une chaine de caracteres (le nom du fichier), soit un entier (l'unite de lecture). Dans le cas ou l'on donne le nom du fichier, le code lui affecte la deuxieme unite de lecture disponible a partir de 80. Cette unite sera liberee a la sortie de l'operateur. Dans le cas ou les composants et les logk sont sur le meme fichier il est inutile de preciser 'LOGK' VAL2, toutefois si ces donnees figurent,le code reconnaitra que les fichiers ont le meme nom et lira en sequence. ENTH est un mot-cle VAL3 est l'adresse du fichier de la base de donnees, permettant de prendre en compte les effets thermiques. Cette adresse peut etre,soit une chaine de caracteres (le nom du fichier),soit un entier (l'unite de lecture). Dans le cas ou l'on donne le nom du fichier, le code lui affecte la troisieme unite de lecture disponible a partir de 80. Cette unite sera liberee a la sortie de l'operateur. Dans le cas ou ces donnees et les logk sont sur le meme fichier il est inutile de preciser 'ENTH' VAL3, toutefois si VAL2 et VAL3 sont identiques,le code lira en sequence. TAB2 est une TABLE contenant le systeme chimique. Cette TABLE a pour sous type CHIMI1, et elle est indicee par les mots: 'DESCHI' 'IDEN' <'REDOX'> <'TEMPE'> TAB2.DESCHI est une TABLE contenant la description du systeme chimique. Elle est indicee par les mots: 'IDX' 'NOM' 'CHARGE' 'IDY' 'NN' 'MATRICEA' 'LOGK' 'IDSURF' <'IDZ' 'MATRICEF' 'IDP'> TAB2.DESCHI.IDX objet LISTENTI contenant les identifiants des composants retenus TAB2.DESCHI.NOM objet LISTMOTS contenant les noms de ces composants (sur 32 caracteres). TAB2.DESCHI.CHARGE objet LISTENTI contenant la charge de ces composants. TAB2.DESCHI.IDY objet LISTENTI contenant les identifiants des especes retenues. TAB2.DESCHI.NOMESPECE objet LISTMOTS contenant les noms de ces especes (sur 32 caracteres). TAB2.DESCHI.NN objet LISTENTI de 6 valeurs contenant le nombre d'especes de chaque type c'est a dire: nombre d'especes simples en solution (type 1) nombre de complexes en solution (type 2) nombre d'especes a activite fixee (type 3) nombre de mineraux precipites (type 4) nombre de mineraux precipitables en solution(type 5) nombre d'especes non prises en compte dans le calcul (type 6) TAB2.DESCHI.MATRICEA objet LISTREEL matrice contenant les coefficients stoechiometriques de chaque espece. TAB2.DESCHI.LOGK objet LISTREEL valeur pour chaque espece de la constante de formation logk. TAB2.DESCHI.IDSURF objet LISTENTI nature des especes de type 2. 0 en solution ,1 echange ionique (fixee) ,2 adsorbee TAB2.DESCHI.IDZ objet LISTENTI contenant les identifiants des solutions solides retenues. TAB2.DESCHI.IDP objet LISTENTI contenant les identifiants des poles de solutions solides retenus. TAB2.DESCHI.MATRICEF objet LISTREEL matricie contenant les fractions molaires de chaque solution solide TAB2.IDEN est une TABLE indicee par les mots: 'COMP' 'SOLU' 'NOMSOLU' <'PRECI' 'NOMPRECI'> <'SURF' 'NOMSURF'> <'TYP3' 'NOMTYP3'> <IMP3> <'TYP6' 'NOMTYP6'> <'PARF'> <'REAC'> <'IMMO'> <'SOSO' 'NOMSOSO' 'POLE' 'NOMPOLE'> TAB2.IDEN.COMP objet de type LISTENTI entiers servant a decrire la nature des composants. 1 conservatifs, 2 reactifs, 3 immobiles. TAB2.IDEN.SOLU objet LISTENTI identifiants des especes en solution TAB2.IDEN.NOMSOLU objet LISTMOTS noms attribues par le code a ces especes en solution. Ces noms ont 4 caracteres dont le premier est W suivi eventuellement de 0 ou 00 et du numero d'ordre dans la liste TAB2.DESCHI.IDY de l'identifiant concerne. TAB2.IDEN.PRECI objet LISTENTI identifiants des especes precipites potentiels. TAB2.IDEN.NOMPRECI objet LISTMOTS noms attribues par le code a ces especes. Ces noms ont 4 caracteres dont le premier est W suivi eventuellement de 0 ou 00 et du numero d'ordre dans la liste TAB2.DESCHI.IDY de l'identifiant concerne. TAB2.IDEN.SURF objet LISTENTI identifiants des especes de surface. TAB2.IDEN.NOMSURF objet LISTMOTS noms attribues par le code a ces especes. TAB2.IDEN.TYP3 objet LISTENTI identifiants des especes dont l activite est imposee(meme liste que TAB1.CLIM.TYP3) TAB2.IDEN.NOMTYP3 objet LISTMOTS noms attribues par le code a ces especes. TAB2.IDEN.IMP3 objet LISTENTI identifiants des especes de type 3. (activite imposee, redox). TAB2.IDEN.TYP6 objet LISTENTI identifiants des especes non prises en compte dans le calcul. TAB2.IDEN.NOMTYP6 objet LISTMOTS noms attribues par le code a ces especes. TAB2.IDEN.PARF objet LISTENTI identifiants des composants non reactifs TAB2.IDEN.REAC objet LISTENTI identifiants des composants reactifs TAB2.IDEN.IMMO objet LISTENTI identifiants des composants non transportes( par exemple:composant de surface, composant associe a une espece a activite imposee ) TAB2.IDEN.SOSO objet LISTENTI identifiants des solutions solides. TAB2.IDEN.NOMSOSO objet LISTMOTS noms attribues par le code a ces solutions solides. TAB2.IDEN.POLE objet LISTENTI identifiant des poles de solutions solides. TAB2.IDEN.NOMPOLE objet LISTMOTS noms attribues par le code a ces poles de solutions solides. TAB2.REDOX objet TABLE contenant les donnees relatives aux redox. Cette TABLE est indicee par les mots: 'I1' 'I2' 'A1' 'A2' TAB2.REDOX.I1 objet LISTENTI identifiant du composant mobile TAB2.REDOX.I2 objet LISTENTI identifiant du composant lie au composant mobile. TAB2.REDOX.A1 objet LISTREEL coefficient stoechiometrique de l'electron e- TAB2.REDOX.A2 objet LISTREEL coefficient stoechiometrique du proton H+ TAB2.TEMPE est une TABLE contenant les donnees relatives a la temperature pour les especes selectionnees. Elle est indicee par les mots: 'DELTAH' 'DELCP' 'APPROX' pour la base MINEQL 'COEF' 'LOGK' 'TMIN' 'TMAX' 'NVT' 'NUMT' pour la base STRASBG TAB2.TEMPE.DELTAH objet LISTREEL. Variation d'enthalpie TAB2.TEMPE.DELCP objet LISTREEL. Variation de capacite calorifique TAB2.TEMPE.APPROX entier (1 ou 2) type d'approximation. TAB2.TEMPE.COEF objet LISTREEL.Coefficients d'interpolation. Pour chaque espece retenue,il y a 4 valeurs pour chacune des 8 temperatures de reference. TAB2.TEMPE.TMINA objet LISTREEL. Pour chaque espece retenue, temperature auquelles on donne le LogK.( cette liste est dimensionnee au (nombre d'especes)*(max de valeurs) TAB2.TEMPE.LOGK objet LISTREEL. Valeurs des Logk aux temperatures de la liste TMIMA. TAB2.TEMPE.NVT objet LISTENTI. Pour chaque espece retenue, nombre de temperatures disponibles. TAB2.TEMPE.NUMT objet LISTENTI. Pour chaque espece retenue, identifiant de l'espece.
© Cast3M 2003 - Tous droits réservés.
Mentions légales