Télécharger chi1.notice

Retour à la liste thématique

Afficher cette notice en

Numérotation des lignes :
$$$$ 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