max/iphreeqc
1
0
Fork 0
mirror of https://git.gfz-potsdam.de/naaice/iphreeqc.git synced 2025-12-16 00:28:23 +01:00
Code Issues Packages Projects Releases Wiki Activity
iphreeqc/ex18.out
Charlton, Scott R a9d1de8fa2 Updated output files (from windows)
2024-04-18 13:41:57 -06:00

585 lines
29 KiB
Plaintext
Raw Blame History

Input file: ../examples/ex18
Output file: ex18.out
Database file: ../database/phreeqc.dat
------------------
Reading data base.
------------------
SOLUTION_MASTER_SPECIES
SOLUTION_SPECIES
PHASES
EXCHANGE_MASTER_SPECIES
EXCHANGE_SPECIES
SURFACE_MASTER_SPECIES
SURFACE_SPECIES
CALCULATE_VALUES
RATES
END
------------------------------------
Reading input data for simulation 1.
------------------------------------
TITLE Example 18.--Inverse modeling of Madison aquifer
SOLUTION 1 Recharge number 3
units mmol/kgw
temp 9.9
pe 0.
pH 7.55
Ca 1.2
Mg 1.01
Na 0.02
K 0.02
Fe(2) 0.001
Cl 0.02
S(6) 0.16
S(-2) 0
C(4) 4.30
isotope 13C -7.0 1.4
isotope 34S 9.7 0.9
SOLUTION 2 Mysse
units mmol/kgw
temp 63.
pH 6.61
pe 0.
redox S(6)/S(-2)
Ca 11.28
Mg 4.54
Na 31.89
K 2.54
Fe(2) 0.0004
Cl 17.85
S(6) 19.86
S(-2) 0.26
C(4) 6.87
isotope 13C -2.3 0.2
isotope 34S(6) 16.3 1.5
isotope 34S(-2) -22.1 7
INVERSE_MODELING 1
solutions 1 2
uncertainty 0.05
range
isotopes
13C
34S
balances
Fe(2) 1.0
ph 0.1
phases
Dolomite dis 13C 3.0 2
Calcite pre 13C -1.5 1
Anhydrite dis 34S 13.5 2
CH2O dis 13C -25.0 5
Goethite
Pyrite pre 34S -22. 2
CaX2 pre
Ca.75Mg.25X2 pre
MgX2 pre
NaX
Halite
Sylvite
PHASES
Sylvite
KCl = K+ + Cl-
log_k 0.0
CH2O
CH2O + H2O = CO2 + 4H+ + 4e-
log_k 0.0
EXCHANGE_SPECIES
0.75Ca+2 + 0.25Mg+2 + 2X- = Ca.75Mg.25X2
log_k 0.0
END
-----
TITLE
-----
Example 18.--Inverse modeling of Madison aquifer
-------------------------------------------
Beginning of initial solution calculations.
-------------------------------------------
Initial solution 1. Recharge number 3
-----------------------------Solution composition------------------------------
Elements Molality Moles
C(4) 4.300e-03 4.300e-03
Ca 1.200e-03 1.200e-03
Cl 2.000e-05 2.000e-05
Fe(2) 1.000e-06 1.000e-06
K 2.000e-05 2.000e-05
Mg 1.010e-03 1.010e-03
Na 2.000e-05 2.000e-05
S(6) 1.600e-04 1.600e-04
----------------------------Description of solution----------------------------
pH = 7.550
pe = 0.000
Specific Conductance (µS/cm, 10°C) = 277
Density (g/cm³) = 0.99999
Volume (L) = 1.00035
Viscosity (mPa s) = 1.31356
Activity of water = 1.000
Ionic strength (mol/kgw) = 6.542e-03
Mass of water (kg) = 1.000e+00
Total alkalinity (eq/kg) = 4.016e-03
Total CO2 (mol/kg) = 4.300e-03
Temperature (°C) = 9.90
Electrical balance (eq) = 1.063e-04
Percent error, 100*(Cat-|An|)/(Cat+|An|) = 1.24
Iterations = 8
Total H = 1.110164e+02
Total O = 5.551946e+01
----------------------------Distribution of species----------------------------
Log Log Log mole V
Species Molality Activity Molality Activity Gamma cm³/mol
OH- 1.126e-07 1.034e-07 -6.948 -6.985 -0.037 -4.98
H+ 3.038e-08 2.818e-08 -7.517 -7.550 -0.033 0.00
H2O 5.551e+01 9.999e-01 1.744 -0.000 0.000 18.02
C(4) 4.300e-03
HCO3- 3.929e-03 3.622e-03 -2.406 -2.441 -0.035 22.89
CO2 2.971e-04 2.974e-04 -3.527 -3.527 0.000 33.66
MgHCO3+ 3.094e-05 2.845e-05 -4.510 -4.546 -0.036 4.93
CaHCO3+ 3.027e-05 2.794e-05 -4.519 -4.554 -0.035 8.96
CO3-2 5.765e-06 4.167e-06 -5.239 -5.380 -0.141 -6.10
CaCO3 4.730e-06 4.737e-06 -5.325 -5.324 0.001 -14.66
MgCO3 2.204e-06 2.208e-06 -5.657 -5.656 0.001 -17.07
FeHCO3+ 2.086e-07 1.919e-07 -6.681 -6.717 -0.036 (0)
FeCO3 5.288e-08 5.296e-08 -7.277 -7.276 0.001 (0)
NaHCO3 3.675e-08 3.686e-08 -7.435 -7.433 0.001 31.86
KHCO3 2.291e-08 2.291e-08 -7.640 -7.640 0.000 40.92
(CO2)2 9.408e-10 9.423e-10 -9.026 -9.026 0.001 67.31
Ca 1.200e-03
Ca+2 1.152e-03 8.324e-04 -2.938 -3.080 -0.141 -18.31
CaHCO3+ 3.027e-05 2.794e-05 -4.519 -4.554 -0.035 8.96
CaSO4 1.277e-05 1.279e-05 -4.894 -4.893 0.001 6.78
CaCO3 4.730e-06 4.737e-06 -5.325 -5.324 0.001 -14.66
CaOH+ 5.328e-09 4.901e-09 -8.273 -8.310 -0.036 (0)
CaHSO4+ 2.133e-12 1.962e-12 -11.671 -11.707 -0.036 (0)
Cl 2.000e-05
Cl- 2.000e-05 1.838e-05 -4.699 -4.736 -0.037 17.40
FeCl+ 1.461e-11 1.344e-11 -10.835 -10.872 -0.036 (0)
HCl 1.948e-13 1.960e-13 -12.711 -12.708 0.003 (0)
Fe(2) 1.000e-06
Fe+2 7.297e-07 5.298e-07 -6.137 -6.276 -0.139 -23.20
FeHCO3+ 2.086e-07 1.919e-07 -6.681 -6.717 -0.036 (0)
FeCO3 5.288e-08 5.296e-08 -7.277 -7.276 0.001 (0)
FeSO4 6.844e-09 6.855e-09 -8.165 -8.164 0.001 39.09
FeOH+ 1.965e-09 1.811e-09 -8.707 -8.742 -0.036 (0)
FeCl+ 1.461e-11 1.344e-11 -10.835 -10.872 -0.036 (0)
Fe(OH)2 1.369e-13 1.371e-13 -12.864 -12.863 0.001 (0)
FeHSO4+ 1.357e-15 1.249e-15 -14.867 -14.904 -0.036 (0)
Fe(OH)3- 1.678e-16 1.546e-16 -15.775 -15.811 -0.036 (0)
H(0) 1.316e-18
H2 6.579e-19 6.588e-19 -18.182 -18.181 0.001 28.63
K 2.000e-05
K+ 1.995e-05 1.833e-05 -4.700 -4.737 -0.037 8.43
KSO4- 2.828e-08 2.615e-08 -7.548 -7.582 -0.034 12.38
KHCO3 2.291e-08 2.291e-08 -7.640 -7.640 0.000 40.92
Mg 1.010e-03
Mg+2 9.648e-04 7.000e-04 -3.016 -3.155 -0.139 -21.10
MgHCO3+ 3.094e-05 2.845e-05 -4.510 -4.546 -0.036 4.93
MgSO4 1.203e-05 1.206e-05 -4.920 -4.919 0.001 -14.55
MgCO3 2.204e-06 2.208e-06 -5.657 -5.656 0.001 -17.07
MgOH+ 2.321e-08 2.144e-08 -7.634 -7.669 -0.034 (0)
Mg(SO4)2-2 7.109e-09 5.202e-09 -8.148 -8.284 -0.136 1.52
Na 2.000e-05
Na+ 1.994e-05 1.836e-05 -4.700 -4.736 -0.036 -2.42
NaHCO3 3.675e-08 3.686e-08 -7.435 -7.433 0.001 31.86
NaSO4- 2.687e-08 2.478e-08 -7.571 -7.606 -0.035 -19.98
NaOH 1.896e-22 1.899e-22 -21.722 -21.721 0.001 (0)
O(0) 0.000e+00
O2 0.000e+00 0.000e+00 -61.151 -61.151 0.001 28.94
S(6) 1.600e-04
SO4-2 1.351e-04 9.732e-05 -3.869 -4.012 -0.143 14.48
CaSO4 1.277e-05 1.279e-05 -4.894 -4.893 0.001 6.78
MgSO4 1.203e-05 1.206e-05 -4.920 -4.919 0.001 -14.55
KSO4- 2.828e-08 2.615e-08 -7.548 -7.582 -0.034 12.38
NaSO4- 2.687e-08 2.478e-08 -7.571 -7.606 -0.035 -19.98
Mg(SO4)2-2 7.109e-09 5.202e-09 -8.148 -8.284 -0.136 1.52
FeSO4 6.844e-09 6.855e-09 -8.165 -8.164 0.001 39.09
HSO4- 2.131e-10 1.960e-10 -9.671 -9.708 -0.036 38.92
CaHSO4+ 2.133e-12 1.962e-12 -11.671 -11.707 -0.036 (0)
FeHSO4+ 1.357e-15 1.249e-15 -14.867 -14.904 -0.036 (0)
------------------------------Saturation indices-------------------------------
Phase SI** log IAP log K(283 K, 1 atm)
Anhydrite -2.97 -7.09 -4.12 CaSO4
Aragonite -0.20 -8.46 -8.25 CaCO3
Arcanite -11.37 -13.49 -2.11 K2SO4
Calcite -0.06 -8.46 -8.40 CaCO3
CH2O -33.73 -33.73 0.00 CH2O
CO2(g) -2.26 -3.53 -1.27 CO2
Dolomite -0.11 -16.99 -16.88 CaMg(CO3)2
Epsomite -5.33 -7.17 -1.83 MgSO4:7H2O
Gypsum -2.49 -7.09 -4.60 CaSO4:2H2O
H2(g) -15.13 -18.18 -3.05 H2
H2O(g) -1.91 -0.00 1.91 H2O
Halite -11.03 -9.47 1.56 NaCl
Hexahydrite -5.58 -7.17 -1.59 MgSO4:6H2O
Kieserite -5.85 -7.17 -1.32 MgSO4:H2O
Melanterite -7.88 -10.29 -2.41 FeSO4:7H2O
Mirabilite -11.51 -13.48 -1.97 Na2SO4:10H2O
O2(g) -58.39 -61.15 -2.76 O2
Siderite -0.86 -11.66 -10.79 FeCO3
Sylvite -9.47 -9.47 0.00 KCl
Thenardite -13.26 -13.48 -0.23 Na2SO4
**For a gas, SI = log10(fugacity). Fugacity = pressure * phi / 1 atm.
For ideal gases, phi = 1.
Initial solution 2. Mysse
-----------------------------Solution composition------------------------------
Elements Molality Moles
C(4) 6.870e-03 6.870e-03
Ca 1.128e-02 1.128e-02
Cl 1.785e-02 1.785e-02
Fe(2) 4.000e-07 4.000e-07
K 2.540e-03 2.540e-03
Mg 4.540e-03 4.540e-03
Na 3.189e-02 3.189e-02
S(-2) 2.600e-04 2.600e-04
S(6) 1.986e-02 1.986e-02
----------------------------Description of solution----------------------------
pH = 6.610
pe = 0.000
Specific Conductance (µS/cm, 63°C) = 10510
Density (g/cm³) = 0.98523
Volume (L) = 1.01939
Viscosity (mPa s) = 0.45611
Activity of water = 0.999
Ionic strength (mol/kgw) = 7.104e-02
Mass of water (kg) = 1.000e+00
Total alkalinity (eq/kg) = 5.286e-03
Total CO2 (mol/kg) = 6.870e-03
Temperature (°C) = 63.00
Electrical balance (eq) = 3.215e-03
Percent error, 100*(Cat-|An|)/(Cat+|An|) = 3.23
Iterations = 8 (16 overall)
Total H = 1.110179e+02
Total O = 5.560451e+01
---------------------------------Redox couples---------------------------------
Redox couple pe Eh (volts)
S(-2)/S(6) -3.6540 -0.2437
----------------------------Distribution of species----------------------------
Log Log Log mole V
Species Molality Activity Molality Activity Gamma cm³/mol
OH- 5.957e-07 4.608e-07 -6.225 -6.337 -0.112 -3.63
H+ 2.959e-07 2.455e-07 -6.529 -6.610 -0.081 0.00
H2O 5.551e+01 9.985e-01 1.744 -0.001 0.000 18.35
C(4) 6.870e-03
HCO3- 4.595e-03 3.660e-03 -2.338 -2.436 -0.099 25.86
CO2 1.753e-03 1.772e-03 -2.756 -2.751 0.005 36.36
CaHCO3+ 2.484e-04 1.994e-04 -3.605 -3.700 -0.095 10.66
NaHCO3 1.852e-04 1.914e-04 -3.732 -3.718 0.014 31.20
MgHCO3+ 6.121e-05 4.775e-05 -4.213 -4.321 -0.108 6.10
CaCO3 1.668e-05 1.695e-05 -4.778 -4.771 0.007 -14.51
KHCO3 5.118e-06 5.126e-06 -5.291 -5.290 0.001 41.50
CO3-2 2.699e-06 1.086e-06 -5.569 -5.964 -0.395 -2.58
MgCO3 1.518e-06 1.543e-06 -5.819 -5.812 0.007 -17.09
(CO2)2 1.559e-07 1.584e-07 -6.807 -6.800 0.007 72.72
FeHCO3+ 1.158e-08 9.117e-09 -7.936 -8.040 -0.104 (0)
FeCO3 6.385e-10 6.490e-10 -9.195 -9.188 0.007 (0)
Ca 1.128e-02
Ca+2 7.843e-03 3.171e-03 -2.106 -2.499 -0.393 -17.58
CaSO4 3.172e-03 3.224e-03 -2.499 -2.492 0.007 8.42
CaHCO3+ 2.484e-04 1.994e-04 -3.605 -3.700 -0.095 10.66
CaCO3 1.668e-05 1.695e-05 -4.778 -4.771 0.007 -14.51
CaHSO4+ 1.310e-08 1.031e-08 -7.883 -7.987 -0.104 (0)
CaOH+ 2.719e-09 2.141e-09 -8.566 -8.669 -0.104 (0)
Cl 1.785e-02
Cl- 1.785e-02 1.388e-02 -1.748 -1.858 -0.109 18.29
HCl 8.646e-10 9.269e-10 -9.063 -9.033 0.030 (0)
FeCl+ 6.060e-10 4.771e-10 -9.218 -9.321 -0.104 (0)
Fe(2) 4.000e-07
Fe(HS)2 2.863e-07 2.910e-07 -6.543 -6.536 0.007 (0)
Fe+2 5.995e-08 2.491e-08 -7.222 -7.604 -0.381 -19.95
FeSO4 3.584e-08 3.643e-08 -7.446 -7.439 0.007 -6.81
FeHCO3+ 1.158e-08 9.117e-09 -7.936 -8.040 -0.104 (0)
Fe(HS)3- 4.609e-09 3.628e-09 -8.336 -8.440 -0.104 (0)
FeCO3 6.385e-10 6.490e-10 -9.195 -9.188 0.007 (0)
FeCl+ 6.060e-10 4.771e-10 -9.218 -9.321 -0.104 (0)
FeOH+ 5.020e-10 3.976e-10 -9.299 -9.401 -0.101 (0)
Fe(OH)2 2.540e-13 2.582e-13 -12.595 -12.588 0.007 (0)
FeHSO4+ 1.029e-13 8.100e-14 -12.988 -13.092 -0.104 (0)
Fe(OH)3- 6.860e-17 5.434e-17 -16.164 -16.265 -0.101 (0)
H(0) 1.219e-09
H2 6.097e-10 6.198e-10 -9.215 -9.208 0.007 28.58
K 2.540e-03
K+ 2.340e-03 1.814e-03 -2.631 -2.741 -0.111 10.01
KSO4- 1.954e-04 1.605e-04 -3.709 -3.794 -0.085 15.45
KHCO3 5.118e-06 5.126e-06 -5.291 -5.290 0.001 41.50
Mg 4.540e-03
MgSO4 2.446e-03 2.527e-03 -2.612 -2.597 0.014 -0.04
Mg+2 1.986e-03 8.303e-04 -2.702 -3.081 -0.379 -22.49
MgHCO3+ 6.121e-05 4.775e-05 -4.213 -4.321 -0.108 6.10
Mg(SO4)2-2 4.582e-05 2.013e-05 -4.339 -4.696 -0.357 26.89
MgCO3 1.518e-06 1.543e-06 -5.819 -5.812 0.007 -17.09
MgOH+ 3.185e-07 2.573e-07 -6.497 -6.590 -0.093 (0)
Na 3.189e-02
Na+ 2.911e-02 2.304e-02 -1.536 -1.638 -0.102 0.09
NaSO4- 2.598e-03 2.072e-03 -2.585 -2.684 -0.098 -3.67
NaHCO3 1.852e-04 1.914e-04 -3.732 -3.718 0.014 31.20
NaOH 1.044e-18 1.062e-18 -17.981 -17.974 0.007 (0)
O(0) 0.000e+00
O2 0.000e+00 0.000e+00 -63.096 -63.089 0.007 32.51
S(-2) 2.600e-04
HS- 1.480e-04 1.145e-04 -3.830 -3.941 -0.112 21.38
H2S 1.114e-04 1.132e-04 -3.953 -3.946 0.007 42.59
Fe(HS)2 2.863e-07 2.910e-07 -6.543 -6.536 0.007 (0)
Fe(HS)3- 4.609e-09 3.628e-09 -8.336 -8.440 -0.104 (0)
S-2 1.440e-09 5.667e-10 -8.842 -9.247 -0.405 (0)
(H2S)2 1.351e-09 1.373e-09 -8.869 -8.862 0.007 27.88
S(6) 1.986e-02
SO4-2 1.136e-02 4.441e-03 -1.945 -2.353 -0.408 23.91
CaSO4 3.172e-03 3.224e-03 -2.499 -2.492 0.007 8.42
NaSO4- 2.598e-03 2.072e-03 -2.585 -2.684 -0.098 -3.67
MgSO4 2.446e-03 2.527e-03 -2.612 -2.597 0.014 -0.04
KSO4- 1.954e-04 1.605e-04 -3.709 -3.794 -0.085 15.45
Mg(SO4)2-2 4.582e-05 2.013e-05 -4.339 -4.696 -0.357 26.89
HSO4- 3.436e-07 2.705e-07 -6.464 -6.568 -0.104 41.70
FeSO4 3.584e-08 3.643e-08 -7.446 -7.439 0.007 -6.81
CaHSO4+ 1.310e-08 1.031e-08 -7.883 -7.987 -0.104 (0)
FeHSO4+ 1.029e-13 8.100e-14 -12.988 -13.092 -0.104 (0)
------------------------------Saturation indices-------------------------------
Phase SI** log IAP log K(336 K, 1 atm)
Anhydrite -0.11 -4.85 -4.75 CaSO4
Aragonite 0.21 -8.46 -8.67 CaCO3
Arcanite -6.28 -7.84 -1.55 K2SO4
Calcite 0.46 -8.46 -8.92 CaCO3
CH2O -14.57 -14.57 0.00 CH2O
CO2(g) -0.95 -2.75 -1.80 CO2
Dolomite 0.55 -17.51 -18.06 CaMg(CO3)2
Epsomite -3.94 -5.44 -1.50 MgSO4:7H2O
FeS(ppt) -1.02 -4.93 -3.92 FeS
Gypsum -0.19 -4.85 -4.66 CaSO4:2H2O
H2(g) -6.07 -9.21 -3.14 H2
H2O(g) -0.65 -0.00 0.65 H2O
H2S(g) -2.63 -10.55 -7.92 H2S
Halite -5.09 -3.50 1.60 NaCl
Hexahydrite -3.92 -5.44 -1.51 MgSO4:6H2O
Kieserite -4.11 -5.43 -1.32 MgSO4:H2O
Mackinawite -0.29 -4.93 -4.65 FeS
Melanterite -8.11 -9.96 -1.85 FeSO4:7H2O
Mirabilite -5.58 -5.63 -0.06 Na2SO4:10H2O
O2(g) -60.02 -63.09 -3.07 O2
Pyrite 7.97 -9.57 -17.54 FeS2
Siderite -2.47 -13.57 -11.10 FeCO3
Sulfur -2.13 1.97 4.09 S
Sylvite -4.60 -4.60 0.00 KCl
Thenardite -5.14 -5.63 -0.49 Na2SO4
**For a gas, SI = log10(fugacity). Fugacity = pressure * phi / 1 atm.
For ideal gases, phi = 1.
---------------------------------------------
Beginning of inverse modeling 1 calculations.
---------------------------------------------
Using Cl1 standard precision optimization routine.
Solution 1: Recharge number 3
Input Delta Input+Delta
pH 7.550e+00 + 0.000e+00 = 7.550e+00
Alkalinity 4.016e-03 + 1.063e-04 = 4.122e-03
C(-4) 0.000e+00 + 0.000e+00 = 0.000e+00
C(4) 4.300e-03 + 0.000e+00 = 4.300e-03
Ca 1.200e-03 + 0.000e+00 = 1.200e-03
Cl 2.000e-05 + 0.000e+00 = 2.000e-05
Fe(2) 1.000e-06 + 0.000e+00 = 1.000e-06
Fe(3) 0.000e+00 + 0.000e+00 = 0.000e+00
H(0) 0.000e+00 + 0.000e+00 = 0.000e+00
K 2.000e-05 + 0.000e+00 = 2.000e-05
Mg 1.010e-03 + 0.000e+00 = 1.010e-03
Na 2.000e-05 + 0.000e+00 = 2.000e-05
O(0) 0.000e+00 + 0.000e+00 = 0.000e+00
S(-2) 0.000e+00 + 0.000e+00 = 0.000e+00
S(6) 1.600e-04 + 0.000e+00 = 1.600e-04
X 0.000e+00 + 0.000e+00 = 0.000e+00
13C(-4) -7 + 0 = -7
13C(4) -7 + 0 = -7
34S(-2) 9.7 + 0 = 9.7
34S(6) 9.7 + 0 = 9.7
Solution 2: Mysse
Input Delta Input+Delta
pH 6.610e+00 + 0.000e+00 = 6.610e+00
Alkalinity 5.286e-03 + 0.000e+00 = 5.286e-03
C(-4) 0.000e+00 + 0.000e+00 = 0.000e+00
C(4) 6.870e-03 + 0.000e+00 = 6.870e-03
Ca 1.128e-02 + 0.000e+00 = 1.128e-02
Cl 1.785e-02 + 0.000e+00 = 1.785e-02
Fe(2) 4.000e-07 + 0.000e+00 = 4.000e-07
Fe(3) 0.000e+00 + 0.000e+00 = 0.000e+00
H(0) 1.219e-09 + 0.000e+00 = 1.219e-09
K 2.540e-03 + 0.000e+00 = 2.540e-03
Mg 4.540e-03 + 0.000e+00 = 4.540e-03
Na 3.189e-02 + -1.229e-03 = 3.066e-02
O(0) 0.000e+00 + 0.000e+00 = 0.000e+00
S(-2) 2.600e-04 + 0.000e+00 = 2.600e-04
S(6) 1.986e-02 + 9.930e-04 = 2.085e-02
X 0.000e+00 + 0.000e+00 = 0.000e+00
13C(-4) -2.3 + 0 = -2.3
13C(4) -2.3 + 0 = -2.3
34S(-2) -22.1 + 0 = -22.1
34S(6) 16.3 + 0 = 16.3
Isotopic composition of phases:
13C Dolomite 3 + 0 = 3
13C Calcite -1.5 + 0 = -1.5
34S Anhydrite 13.5 + -0.715402 = 12.7846
13C CH2O -25 + 3.93324 = -21.0668
34S Pyrite -22 + 2 = -20
Solution fractions: Minimum Maximum
Solution 1 1.000e+00 9.999e-01 1.000e+00
Solution 2 1.000e+00 1.000e+00 1.000e+00
Phase mole transfers: Minimum Maximum Formula (Approximate SI in solution 1, 2 at 336 K, 1 atm)
Dolomite 1.120e-02 1.022e-02 1.194e-02 CaMg(CO3)2 ( 1.06, 0.55)
Calcite -2.404e-02 -2.597e-02 -2.113e-02 CaCO3 ( 0.46, 0.46)
Anhydrite 2.293e-02 2.037e-02 2.378e-02 CaSO4 ( -2.35, -0.11)
CH2O 4.222e-03 2.482e-03 5.808e-03 CH2O ( , )
Goethite 9.866e-04 5.145e-04 1.418e-03 FeOOH ( , )
Pyrite -9.872e-04 -1.417e-03 -5.162e-04 FeS2 ( , )
MgX2 -7.666e-03 -8.591e-03 -6.972e-03 MgX2 ( , )
NaX 1.533e-02 1.394e-02 1.718e-02 NaX ( , )
Halite 1.531e-02 1.429e-02 1.633e-02 NaCl (-11.07, -5.09)
Sylvite 2.520e-03 2.392e-03 2.648e-03 KCl ( -9.47, -4.60)
Redox mole transfers:
Fe(3) 9.866e-04
H(0) -1.219e-09
S(-2) -2.234e-03
Sum of residuals (epsilons in documentation): 2.312e+00
Sum of delta/uncertainty limit: 4.444e+00
Maximum fractional error in element concentration: 5.000e-02
Model contains minimum number of phases.
===============================================================================
Solution 1: Recharge number 3
Input Delta Input+Delta
pH 7.550e+00 + 0.000e+00 = 7.550e+00
Alkalinity 4.016e-03 + 1.063e-04 = 4.122e-03
C(-4) 0.000e+00 + 0.000e+00 = 0.000e+00
C(4) 4.300e-03 + 0.000e+00 = 4.300e-03
Ca 1.200e-03 + 0.000e+00 = 1.200e-03
Cl 2.000e-05 + 0.000e+00 = 2.000e-05
Fe(2) 1.000e-06 + 0.000e+00 = 1.000e-06
Fe(3) 0.000e+00 + 0.000e+00 = 0.000e+00
H(0) 0.000e+00 + 0.000e+00 = 0.000e+00
K 2.000e-05 + 0.000e+00 = 2.000e-05
Mg 1.010e-03 + 0.000e+00 = 1.010e-03
Na 2.000e-05 + 0.000e+00 = 2.000e-05
O(0) 0.000e+00 + 0.000e+00 = 0.000e+00
S(-2) 0.000e+00 + 0.000e+00 = 0.000e+00
S(6) 1.600e-04 + 0.000e+00 = 1.600e-04
X 0.000e+00 + 0.000e+00 = 0.000e+00
13C(-4) -7 + 0 = -7
13C(4) -7 + 0 = -7
34S(-2) 9.7 + 0 = 9.7
34S(6) 9.7 + 0 = 9.7
Solution 2: Mysse
Input Delta Input+Delta
pH 6.610e+00 + 0.000e+00 = 6.610e+00
Alkalinity 5.286e-03 + 0.000e+00 = 5.286e-03
C(-4) 0.000e+00 + 0.000e+00 = 0.000e+00
C(4) 6.870e-03 + -3.424e-04 = 6.528e-03
Ca 1.128e-02 + 0.000e+00 = 1.128e-02
Cl 1.785e-02 + 0.000e+00 = 1.785e-02
Fe(2) 4.000e-07 + 0.000e+00 = 4.000e-07
Fe(3) 0.000e+00 + 0.000e+00 = 0.000e+00
H(0) 1.219e-09 + 0.000e+00 = 1.219e-09
K 2.540e-03 + 0.000e+00 = 2.540e-03
Mg 4.540e-03 + 0.000e+00 = 4.540e-03
Na 3.189e-02 + -1.229e-03 = 3.066e-02
O(0) 0.000e+00 + 0.000e+00 = 0.000e+00
S(-2) 2.600e-04 + 0.000e+00 = 2.600e-04
S(6) 1.986e-02 + 9.930e-04 = 2.085e-02
X 0.000e+00 + 0.000e+00 = 0.000e+00
13C(-4) -2.3 + 0 = -2.3
13C(4) -2.3 + 0 = -2.3
34S(-2) -22.1 + 0 = -22.1
34S(6) 16.3 + 0 = 16.3
Isotopic composition of phases:
13C Dolomite 3 + 2 = 5
13C Calcite -1.5 + -1 = -2.5
34S Anhydrite 13.5 + -0.146098 = 13.3539
13C CH2O -25 + 5 = -20
34S Pyrite -22 + 2 = -20
Solution fractions: Minimum Maximum
Solution 1 1.000e+00 1.000e+00 1.000e+00
Solution 2 1.000e+00 1.000e+00 1.000e+00
Phase mole transfers: Minimum Maximum Formula (Approximate SI in solution 1, 2 at 336 K, 1 atm)
Dolomite 5.446e-03 4.995e-03 5.842e-03 CaMg(CO3)2 ( 1.06, 0.55)
Calcite -1.215e-02 -1.337e-02 -1.066e-02 CaCO3 ( 0.46, 0.46)
Anhydrite 2.254e-02 2.037e-02 2.298e-02 CaSO4 ( -2.35, -0.11)
CH2O 3.488e-03 2.482e-03 4.302e-03 CH2O ( , )
Goethite 7.909e-04 5.145e-04 1.016e-03 FeOOH ( , )
Pyrite -7.915e-04 -1.015e-03 -5.162e-04 FeS2 ( , )
Ca.75Mg.25X2 -7.666e-03 -8.591e-03 -6.972e-03 Ca.75Mg.25X2 ( , )
NaX 1.533e-02 1.394e-02 1.718e-02 NaX ( , )
Halite 1.531e-02 1.429e-02 1.633e-02 NaCl (-11.07, -5.09)
Sylvite 2.520e-03 2.392e-03 2.648e-03 KCl ( -9.47, -4.60)
Redox mole transfers:
Fe(3) 7.909e-04
H(0) -1.219e-09
S(-2) -1.843e-03
Sum of residuals (epsilons in documentation): 3.320e+00
Sum of delta/uncertainty limit: 7.370e+00
Maximum fractional error in element concentration: 5.000e-02
Model contains minimum number of phases.
===============================================================================
Summary of inverse modeling:
Number of models found: 2
Number of minimal models found: 2
Number of infeasible sets of phases saved: 28
Number of calls to cl1: 80
------------------
End of simulation.
------------------
------------------------------------
Reading input data for simulation 2.
------------------------------------
Reference in New Issue View Git Blame Copy Permalink