max/iphreeqc
1
0
Fork 0
mirror of https://git.gfz-potsdam.de/naaice/iphreeqc.git synced 2025-12-15 16:18:22 +01:00
Code Issues Packages Projects Releases Wiki Activity
iphreeqc/phreeqc3-examples/examples_pc/ex16.out
Scott Charlton d5345a7aee Add 'phreeqc3-examples/' from commit 'e977363da4c1ae1708ebddca38d5312b6813a337' 
git-subtree-dir: phreeqc3-examples
git-subtree-mainline: 8de0e98ed90bb2b639659c820168ab18c1062992
git-subtree-split: e977363da4c1ae1708ebddca38d5312b6813a337
2018-07-31 16:55:50 -06:00

444 lines
22 KiB
Plaintext
Raw Blame History

Input file: ..\examples\ex16
Output file: ex16.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
RATES
END
------------------------------------
Reading input data for simulation 1.
------------------------------------
TITLE Example 16.--Inverse modeling of Sierra springs
SOLUTION_SPREAD
units mmol/L
Number pH Si Ca Mg Na K Alkalinity S(6) Cl
1 6.2 0.273 0.078 0.029 0.134 0.028 0.328 0.01 0.014
2 6.8 0.41 0.26 0.071 0.259 0.04 0.895 0.025 0.03
INVERSE_MODELING 1
solutions 1 2
uncertainty 0.025
range
phases
Halite
Gypsum
Kaolinite precip
Ca-montmorillonite precip
CO2(g)
Calcite
Chalcedony precip
Biotite dissolve
Plagioclase dissolve
balances
Ca 0.05 0.025
PHASES
Biotite
KMg3AlSi3O10(OH)2 + 6H+ + 4H2O = K+ + 3Mg+2 + Al(OH)4- + 3H4SiO4
log_k 0.0 # No log_k, Inverse modeling only
Plagioclase
Na0.62Ca0.38Al1.38Si2.62O8 + 5.52 H+ + 2.48H2O = 0.62Na+ + 0.38Ca+2 + 1.38Al+3 + 2.62H4SiO4
log_k 0.0 # No log_k, inverse modeling only
END
-----
TITLE
-----
Example 16.--Inverse modeling of Sierra springs
-------------------------------------------
Beginning of initial solution calculations.
-------------------------------------------
Initial solution 1.
-----------------------------Solution composition------------------------------
Elements Molality Moles
Alkalinity 3.280e-004 3.280e-004
Ca 7.800e-005 7.800e-005
Cl 1.400e-005 1.400e-005
K 2.800e-005 2.800e-005
Mg 2.900e-005 2.900e-005
Na 1.340e-004 1.340e-004
S(6) 1.000e-005 1.000e-005
Si 2.730e-004 2.730e-004
----------------------------Description of solution----------------------------
pH = 6.200
pe = 4.000
Specific Conductance (uS/cm, 25 oC) = 37
Density (g/cm3) = 0.99708
Volume (L) = 1.00305
Activity of water = 1.000
Ionic strength = 4.851e-004
Mass of water (kg) = 1.000e+000
Total carbon (mol/kg) = 7.824e-004
Total CO2 (mol/kg) = 7.824e-004
Temperature (deg C) = 25.00
Electrical balance (eq) = 1.400e-005
Percent error, 100*(Cat-|An|)/(Cat+|An|) = 1.90
Iterations = 8
Total H = 1.110139e+002
Total O = 5.550924e+001
----------------------------Distribution of species----------------------------
Log Log Log mole V
Species Molality Activity Molality Activity Gamma cm3/mol
H+ 6.465e-007 6.310e-007 -6.189 -6.200 -0.011 0.00
OH- 1.645e-008 1.604e-008 -7.784 -7.795 -0.011 -4.12
H2O 5.551e+001 1.000e+000 1.744 -0.000 0.000 18.07
C(4) 7.824e-004
CO2 4.539e-004 4.539e-004 -3.343 -3.343 0.000 30.26
HCO3- 3.280e-004 3.200e-004 -3.484 -3.495 -0.011 24.58
CaHCO3+ 2.939e-007 2.867e-007 -6.532 -6.543 -0.011 9.66
NaHCO3 1.049e-007 1.049e-007 -6.979 -6.979 0.000 19.41
MgHCO3+ 1.003e-007 9.778e-008 -6.999 -7.010 -0.011 5.47
CO3-2 2.627e-008 2.378e-008 -7.581 -7.624 -0.043 -4.41
CaCO3 2.806e-009 2.806e-009 -8.552 -8.552 0.000 -14.60
MgCO3 5.927e-010 5.927e-010 -9.227 -9.227 0.000 -17.09
NaCO3- 1.422e-011 1.387e-011 -10.847 -10.858 -0.011 -0.67
Ca 7.800e-005
Ca+2 7.760e-005 7.023e-005 -4.110 -4.153 -0.043 -18.17
CaHCO3+ 2.939e-007 2.867e-007 -6.532 -6.543 -0.011 9.66
CaSO4 1.111e-007 1.111e-007 -6.954 -6.954 0.000 7.50
CaCO3 2.806e-009 2.806e-009 -8.552 -8.552 0.000 -14.60
CaOH+ 1.894e-011 1.847e-011 -10.723 -10.733 -0.011 (0)
CaHSO4+ 4.724e-013 4.606e-013 -12.326 -12.337 -0.011 (0)
Cl 1.400e-005
Cl- 1.400e-005 1.365e-005 -4.854 -4.865 -0.011 18.07
H(0) 5.636e-024
H2 2.818e-024 2.818e-024 -23.550 -23.550 0.000 28.61
K 2.800e-005
K+ 2.800e-005 2.730e-005 -4.553 -4.564 -0.011 9.00
KSO4- 1.749e-009 1.706e-009 -8.757 -8.768 -0.011 (0)
Mg 2.900e-005
Mg+2 2.885e-005 2.612e-005 -4.540 -4.583 -0.043 -21.86
MgHCO3+ 1.003e-007 9.778e-008 -6.999 -7.010 -0.011 5.47
MgSO4 5.444e-008 5.445e-008 -7.264 -7.264 0.000 5.84
MgCO3 5.927e-010 5.927e-010 -9.227 -9.227 0.000 -17.09
MgOH+ 1.541e-010 1.503e-010 -9.812 -9.823 -0.011 (0)
Na 1.340e-004
Na+ 1.339e-004 1.306e-004 -3.873 -3.884 -0.011 -1.40
NaHCO3 1.049e-007 1.049e-007 -6.979 -6.979 0.000 19.41
NaSO4- 5.967e-009 5.820e-009 -8.224 -8.235 -0.011 18.41
NaCO3- 1.422e-011 1.387e-011 -10.847 -10.858 -0.011 -0.67
NaOH 2.095e-022 2.095e-022 -21.679 -21.679 0.000 (0)
O(0) 0.000e+000
O2 0.000e+000 0.000e+000 -45.280 -45.280 0.000 30.40
S(6) 1.000e-005
SO4-2 9.827e-006 8.892e-006 -5.008 -5.051 -0.043 13.98
CaSO4 1.111e-007 1.111e-007 -6.954 -6.954 0.000 7.50
MgSO4 5.444e-008 5.445e-008 -7.264 -7.264 0.000 5.84
NaSO4- 5.967e-009 5.820e-009 -8.224 -8.235 -0.011 18.41
KSO4- 1.749e-009 1.706e-009 -8.757 -8.768 -0.011 (0)
HSO4- 5.594e-010 5.455e-010 -9.252 -9.263 -0.011 40.27
CaHSO4+ 4.724e-013 4.606e-013 -12.326 -12.337 -0.011 (0)
Si 2.730e-004
H4SiO4 2.729e-004 2.730e-004 -3.564 -3.564 0.000 52.08
H3SiO4- 6.542e-008 6.379e-008 -7.184 -7.195 -0.011 27.96
H2SiO4-2 7.598e-015 6.877e-015 -14.119 -14.163 -0.043 (0)
------------------------------Saturation indices-------------------------------
Phase SI log IAP log K(298 K, 1 atm)
Anhydrite -4.93 -9.20 -4.28 CaSO4
Aragonite -3.44 -11.78 -8.34 CaCO3
Calcite -3.30 -11.78 -8.48 CaCO3
Chalcedony -0.01 -3.56 -3.55 SiO2
Chrysotile -15.88 16.32 32.20 Mg3Si2O5(OH)4
CO2(g) -1.88 -3.34 -1.46 CO2
Dolomite -6.89 -23.98 -17.09 CaMg(CO3)2
Gypsum -4.62 -9.20 -4.58 CaSO4:2H2O
H2(g) -20.45 -23.55 -3.10 H2
H2O(g) -1.50 -0.00 1.50 H2O
Halite -10.32 -8.75 1.57 NaCl
O2(g) -42.39 -45.28 -2.89 O2
Quartz 0.42 -3.56 -3.98 SiO2
Sepiolite -10.82 4.94 15.76 Mg2Si3O7.5OH:3H2O
Sepiolite(d) -13.72 4.94 18.66 Mg2Si3O7.5OH:3H2O
SiO2(a) -0.85 -3.56 -2.71 SiO2
Sylvite -10.33 -9.43 0.90 KCl
Talc -12.20 9.20 21.40 Mg3Si4O10(OH)2
Initial solution 2.
-----------------------------Solution composition------------------------------
Elements Molality Moles
Alkalinity 8.951e-004 8.951e-004
Ca 2.600e-004 2.600e-004
Cl 3.000e-005 3.000e-005
K 4.000e-005 4.000e-005
Mg 7.101e-005 7.101e-005
Na 2.590e-004 2.590e-004
S(6) 2.500e-005 2.500e-005
Si 4.100e-004 4.100e-004
----------------------------Description of solution----------------------------
pH = 6.800
pe = 4.000
Specific Conductance (uS/cm, 25 oC) = 95
Density (g/cm3) = 0.99713
Volume (L) = 1.00309
Activity of water = 1.000
Ionic strength = 1.313e-003
Mass of water (kg) = 1.000e+000
Total carbon (mol/kg) = 1.199e-003
Total CO2 (mol/kg) = 1.199e-003
Temperature (deg C) = 25.00
Electrical balance (eq) = -1.400e-005
Percent error, 100*(Cat-|An|)/(Cat+|An|) = -0.73
Iterations = 7
Total H = 1.110150e+002
Total O = 5.551125e+001
----------------------------Distribution of species----------------------------
Log Log Log mole V
Species Molality Activity Molality Activity Gamma cm3/mol
H+ 1.647e-007 1.585e-007 -6.783 -6.800 -0.017 0.00
OH- 6.652e-008 6.386e-008 -7.177 -7.195 -0.018 -4.10
H2O 5.551e+001 1.000e+000 1.744 -0.000 0.000 18.07
C(4) 1.199e-003
HCO3- 8.903e-004 8.554e-004 -3.050 -3.068 -0.017 24.60
CO2 3.047e-004 3.048e-004 -3.516 -3.516 0.000 30.26
CaHCO3+ 2.483e-006 2.387e-006 -5.605 -5.622 -0.017 9.68
MgHCO3+ 6.231e-007 5.983e-007 -6.205 -6.223 -0.018 5.48
NaHCO3 5.333e-007 5.334e-007 -6.273 -6.273 0.000 19.41
CO3-2 2.970e-007 2.531e-007 -6.527 -6.597 -0.069 -4.35
CaCO3 9.296e-008 9.298e-008 -7.032 -7.032 0.000 -14.60
MgCO3 1.443e-008 1.444e-008 -7.841 -7.841 0.000 -17.09
NaCO3- 2.921e-010 2.807e-010 -9.534 -9.552 -0.017 -0.66
Ca 2.600e-004
Ca+2 2.567e-004 2.187e-004 -3.591 -3.660 -0.070 -18.13
CaHCO3+ 2.483e-006 2.387e-006 -5.605 -5.622 -0.017 9.68
CaSO4 7.909e-007 7.911e-007 -6.102 -6.102 0.000 7.50
CaCO3 9.296e-008 9.298e-008 -7.032 -7.032 0.000 -14.60
CaOH+ 2.385e-010 2.290e-010 -9.623 -9.640 -0.018 (0)
CaHSO4+ 8.584e-013 8.242e-013 -12.066 -12.084 -0.018 (0)
Cl 3.000e-005
Cl- 3.000e-005 2.880e-005 -4.523 -4.541 -0.018 18.08
H(0) 3.555e-025
H2 1.778e-025 1.778e-025 -24.750 -24.750 0.000 28.61
K 4.000e-005
K+ 4.000e-005 3.840e-005 -4.398 -4.416 -0.018 9.01
KSO4- 5.713e-009 5.489e-009 -8.243 -8.260 -0.017 (0)
Mg 7.101e-005
Mg+2 7.008e-005 5.977e-005 -4.154 -4.224 -0.069 -21.81
MgHCO3+ 6.231e-007 5.983e-007 -6.205 -6.223 -0.018 5.48
MgSO4 2.850e-007 2.851e-007 -6.545 -6.545 0.000 5.84
MgCO3 1.443e-008 1.444e-008 -7.841 -7.841 0.000 -17.09
MgOH+ 1.424e-009 1.369e-009 -8.846 -8.864 -0.017 (0)
Na 2.590e-004
Na+ 2.585e-004 2.482e-004 -3.588 -3.605 -0.018 -1.38
NaHCO3 5.333e-007 5.334e-007 -6.273 -6.273 0.000 19.41
NaSO4- 2.635e-008 2.531e-008 -7.579 -7.597 -0.017 18.44
NaCO3- 2.921e-010 2.807e-010 -9.534 -9.552 -0.017 -0.66
NaOH 1.585e-021 1.585e-021 -20.800 -20.800 0.000 (0)
O(0) 0.000e+000
O2 0.000e+000 0.000e+000 -42.880 -42.880 0.000 30.40
S(6) 2.500e-005
SO4-2 2.389e-005 2.035e-005 -4.622 -4.692 -0.070 14.03
CaSO4 7.909e-007 7.911e-007 -6.102 -6.102 0.000 7.50
MgSO4 2.850e-007 2.851e-007 -6.545 -6.545 0.000 5.84
NaSO4- 2.635e-008 2.531e-008 -7.579 -7.597 -0.017 18.44
KSO4- 5.713e-009 5.489e-009 -8.243 -8.260 -0.017 (0)
HSO4- 3.265e-010 3.135e-010 -9.486 -9.504 -0.018 40.28
CaHSO4+ 8.584e-013 8.242e-013 -12.066 -12.084 -0.018 (0)
Si 4.100e-004
H4SiO4 4.096e-004 4.098e-004 -3.388 -3.387 0.000 52.08
H3SiO4- 3.970e-007 3.812e-007 -6.401 -6.419 -0.018 27.98
H2SiO4-2 1.920e-013 1.636e-013 -12.717 -12.786 -0.069 (0)
------------------------------Saturation indices-------------------------------
Phase SI log IAP log K(298 K, 1 atm)
Anhydrite -4.07 -8.35 -4.28 CaSO4
Aragonite -1.92 -10.26 -8.34 CaCO3
Calcite -1.78 -10.26 -8.48 CaCO3
Chalcedony 0.16 -3.39 -3.55 SiO2
Chrysotile -10.85 21.35 32.20 Mg3Si2O5(OH)4
CO2(g) -2.06 -3.52 -1.46 CO2
Dolomite -3.99 -21.08 -17.09 CaMg(CO3)2
Gypsum -3.77 -8.35 -4.58 CaSO4:2H2O
H2(g) -21.65 -24.75 -3.10 H2
H2O(g) -1.50 -0.00 1.50 H2O
Halite -9.72 -8.15 1.57 NaCl
O2(g) -39.99 -42.88 -2.89 O2
Quartz 0.59 -3.39 -3.98 SiO2
Sepiolite -7.17 8.59 15.76 Mg2Si3O7.5OH:3H2O
Sepiolite(d) -10.07 8.59 18.66 Mg2Si3O7.5OH:3H2O
SiO2(a) -0.68 -3.39 -2.71 SiO2
Sylvite -9.86 -8.96 0.90 KCl
Talc -6.82 14.58 21.40 Mg3Si4O10(OH)2
---------------------------------------------
Beginning of inverse modeling 1 calculations.
---------------------------------------------
Using Cl1 standard precision optimization routine.
Solution 1:
Input Delta Input+Delta
pH 6.200e+000 + 0.000e+000 = 6.200e+000
Al 0.000e+000 + 0.000e+000 = 0.000e+000
Alkalinity 3.280e-004 + 8.200e-006 = 3.362e-004
C(-4) 0.000e+000 + 0.000e+000 = 0.000e+000
C(4) 7.824e-004 + 0.000e+000 = 7.824e-004
Ca 7.800e-005 + -2.550e-006 = 7.545e-005
Cl 1.400e-005 + 0.000e+000 = 1.400e-005
H(0) 0.000e+000 + 0.000e+000 = 0.000e+000
K 2.800e-005 + -7.000e-007 = 2.730e-005
Mg 2.900e-005 + 0.000e+000 = 2.900e-005
Na 1.340e-004 + 0.000e+000 = 1.340e-004
O(0) 0.000e+000 + 0.000e+000 = 0.000e+000
S(-2) 0.000e+000 + 0.000e+000 = 0.000e+000
S(6) 1.000e-005 + 0.000e+000 = 1.000e-005
Si 2.730e-004 + 0.000e+000 = 2.730e-004
Solution 2:
Input Delta Input+Delta
pH 6.800e+000 + 0.000e+000 = 6.800e+000
Al 0.000e+000 + 0.000e+000 = 0.000e+000
Alkalinity 8.951e-004 + -1.480e-005 = 8.803e-004
C(-4) 0.000e+000 + 0.000e+000 = 0.000e+000
C(4) 1.199e-003 + 0.000e+000 = 1.199e-003
Ca 2.600e-004 + 0.000e+000 = 2.600e-004
Cl 3.000e-005 + 0.000e+000 = 3.000e-005
H(0) 0.000e+000 + 0.000e+000 = 0.000e+000
K 4.000e-005 + 1.000e-006 = 4.100e-005
Mg 7.101e-005 + -8.979e-007 = 7.011e-005
Na 2.590e-004 + 0.000e+000 = 2.590e-004
O(0) 0.000e+000 + 0.000e+000 = 0.000e+000
S(-2) 0.000e+000 + 0.000e+000 = 0.000e+000
S(6) 2.500e-005 + 0.000e+000 = 2.500e-005
Si 4.100e-004 + 0.000e+000 = 4.100e-004
Solution fractions: Minimum Maximum
Solution 1 1.000e+000 1.000e+000 1.000e+000
Solution 2 1.000e+000 1.000e+000 1.000e+000
Phase mole transfers: Minimum Maximum
Halite 1.600e-005 1.490e-005 1.710e-005 NaCl
Gypsum 1.500e-005 1.413e-005 1.588e-005 CaSO4:2H2O
Kaolinite -3.392e-005 -5.587e-005 -1.224e-005 Al2Si2O5(OH)4
Ca-Montmorillon -8.090e-005 -1.100e-004 -5.154e-005 Ca0.165Al2.33Si3.67O10(OH)2
CO2(g) 3.006e-004 2.362e-004 3.655e-004 CO2
Calcite 1.161e-004 1.007e-004 1.309e-004 CaCO3
Biotite 1.370e-005 1.317e-005 1.370e-005 KMg3AlSi3O10(OH)2
Plagioclase 1.758e-004 1.582e-004 1.935e-004 Na0.62Ca0.38Al1.38Si2.62O8
Redox mole transfers:
Sum of residuals (epsilons in documentation): 4.821e+000
Sum of delta/uncertainty limit: 4.821e+000
Maximum fractional error in element concentration: 3.269e-002
Model contains minimum number of phases.
===============================================================================
Solution 1:
Input Delta Input+Delta
pH 6.200e+000 + 0.000e+000 = 6.200e+000
Al 0.000e+000 + 0.000e+000 = 0.000e+000
Alkalinity 3.280e-004 + 8.200e-006 = 3.362e-004
C(-4) 0.000e+000 + 0.000e+000 = 0.000e+000
C(4) 7.824e-004 + 0.000e+000 = 7.824e-004
Ca 7.800e-005 + -2.550e-006 = 7.545e-005
Cl 1.400e-005 + 0.000e+000 = 1.400e-005
H(0) 0.000e+000 + 0.000e+000 = 0.000e+000
K 2.800e-005 + -7.000e-007 = 2.730e-005
Mg 2.900e-005 + 0.000e+000 = 2.900e-005
Na 1.340e-004 + 0.000e+000 = 1.340e-004
O(0) 0.000e+000 + 0.000e+000 = 0.000e+000
S(-2) 0.000e+000 + 0.000e+000 = 0.000e+000
S(6) 1.000e-005 + 0.000e+000 = 1.000e-005
Si 2.730e-004 + 0.000e+000 = 2.730e-004
Solution 2:
Input Delta Input+Delta
pH 6.800e+000 + 0.000e+000 = 6.800e+000
Al 0.000e+000 + 0.000e+000 = 0.000e+000
Alkalinity 8.951e-004 + -1.480e-005 = 8.803e-004
C(-4) 0.000e+000 + 0.000e+000 = 0.000e+000
C(4) 1.199e-003 + 0.000e+000 = 1.199e-003
Ca 2.600e-004 + 0.000e+000 = 2.600e-004
Cl 3.000e-005 + 0.000e+000 = 3.000e-005
H(0) 0.000e+000 + 0.000e+000 = 0.000e+000
K 4.000e-005 + 1.000e-006 = 4.100e-005
Mg 7.101e-005 + -8.980e-007 = 7.011e-005
Na 2.590e-004 + 0.000e+000 = 2.590e-004
O(0) 0.000e+000 + 0.000e+000 = 0.000e+000
S(-2) 0.000e+000 + 0.000e+000 = 0.000e+000
S(6) 2.500e-005 + 0.000e+000 = 2.500e-005
Si 4.100e-004 + 0.000e+000 = 4.100e-004
Solution fractions: Minimum Maximum
Solution 1 1.000e+000 1.000e+000 1.000e+000
Solution 2 1.000e+000 1.000e+000 1.000e+000
Phase mole transfers: Minimum Maximum
Halite 1.600e-005 1.490e-005 1.710e-005 NaCl
Gypsum 1.500e-005 1.413e-005 1.588e-005 CaSO4:2H2O
Kaolinite -1.282e-004 -1.403e-004 -1.159e-004 Al2Si2O5(OH)4
CO2(g) 3.140e-004 2.490e-004 3.794e-004 CO2
Calcite 1.028e-004 8.680e-005 1.182e-004 CaCO3
Chalcedony -1.084e-004 -1.473e-004 -6.906e-005 SiO2
Biotite 1.370e-005 1.317e-005 1.370e-005 KMg3AlSi3O10(OH)2
Plagioclase 1.758e-004 1.582e-004 1.935e-004 Na0.62Ca0.38Al1.38Si2.62O8
Redox mole transfers:
Sum of residuals (epsilons in documentation): 4.821e+000
Sum of delta/uncertainty limit: 4.821e+000
Maximum fractional error in element concentration: 3.269e-002
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: 20
Number of calls to cl1: 62
------------------
End of simulation.
------------------
------------------------------------
Reading input data for simulation 2.
------------------------------------
-------------------------------
End of Run after 0.589 Seconds.
-------------------------------
Reference in New Issue View Git Blame Copy Permalink