1、PHREEQC 实例分析例 1物种形成分析这个例子计算了海水中矿物质的分布以及一组有关矿物在海水中的饱和程度。为了证明如何在这个模型中应用新的元素，将元素铀添加入由 phreeqc.dat 定义的液相模型中wateq.dat 是包含于程序分类中的一个数据库文件，它来自于 WATEQ4F（Ball and Nordstrom, 1991） ，并包含铀 。物质形成计算所需要的数据包括温度、Ph、元素的浓度和 /或其元素的化合价。海水中的这些数据见表10。这个例子计算中输入的数据组见表11。在模拟中所运用的有关计算的注释包含在TITLE 关键字中。SOLUTION 数据块定义了海水的成分。注意：元素
2、的化合价用元素化学符号后面圆括号中的数字表示S(6), N(5), N(-3)和O(0)。表 10海水的成分未指定浓度时，其浓度的单位为 ppm分析的组分 PHREEQC 符号 浓度钙 Ca 412.3镁 Mg 1291.8钠 Na 10768.0钾 K 399.1铁 Fe .002锰 Mn .0002硅石，SiO 2 Si 4.28氯化物 Cl 19353.0碱度，HCO 3- Alkalinity 141.682硫酸盐，SO 42- S(6) 2712.0硝酸盐，NO 3- N(5) .29铵，NH 4+ N(-3) .03铀 U .0033pH，标准单位 pH 8.22pe，无单位 pe
3、 8.451温度， temperature 25.0密度，千克/升 density 1.023用于分配氧化还原元素和计算饱和指数的pe由redox标识符所指定。在这个例子中，用氧化还原电对O(-2)/O(0) 计算的pe值相对应于溶解氧 /水，并且这个pe 适用于需要pe值的所有的计算。如果redox没有指定，那么缺省的值将会是所输入的 pe。缺省的氧化还原标识符可被任何氧化还原元素代替，如输入元素锰时，则输入的pe被用来表示各种化合价状态的锰；输入铀时，这里是氮/铵电对将会用来计算所形成各种价态铀的pe值。数据组中缺省的单位为ppm（ units标识符）。这个缺省值可以替换为任何浓度单位，如
4、指定铀的浓度为ppb来代替ppm 。因为ppm 是一个质量单位，而不是一个摩尔单位，这个程序必须用分子量来将浓度单位转化为摩尔单位。每一种主要物质缺省的分子量在SOLUTION_MASTER_SPECIES输入中指定（缺省数据库phreeqc.dat的值列在表4和附录B 中）。如果提交的分子量数据不同于其缺省值，必须在输入数据的设置中指定适当的分子量。这可以用gfw标识符来完成，在这里输入真正的分子量，转化硝酸盐的分子量为62.0 g/mol，或是更简便的是以as标识符来完成，在这里输入所使用的化学分子式的单位，正如在这个例子中输入的碱和铵是一样的。注意最后给定的溶解氧O(0)的浓度是1ppm
5、的初始估计值，但它的浓度将会得以调整，直到氧气分压的对数达到-0.7。O 2(g)的定义是在缺省数据库文件中在PHASES 输入(附录B) 。当使用相均衡来指定初始浓度正如这个例子中的O(0)，则仅有一种浓度是得以调整。例如，例如石膏被用来调整钙的浓度，钙的浓度会改变，而硫酸盐的浓度却保持不变。表 11 例 1 的输入数据TITLE Example 1.-Add uranium and speciate seawater.SOLUTION 1 SEAWATER FROM NORDSTROM ET AL. (1979)units ppmpH 8.22pe 8.451density 1.023te
6、mp 25.0redox O(0)/O(-2)Ca 412.3Mg 1291.8Na 10768.0K 399.1Fe 0.002Mn 0.0002 peSi 4.28Cl 19353.0Alkalinity 141.682 as HCO3S(6) 2712.0N(5) 0.29 gfw 62.0N(-3) 0.03 as NH4U 3.3 ppb N(5)/N(-3)O(0) 1.0 O2(g) -0.7SOLUTION_MASTER_SPECIESU U+4 0.0 238.0290 238.0290U(4) U+4 0.0 238.0290U(5) UO2+ 0.0 238.0290U(
7、6) UO2+2 0.0 238.0290SOLUTION_SPECIES# primary master species for U# is also secondary master species for U(4)U+4 = U+4log_k 0.0U+4 + 4 H2O = U(OH)4 + 4 H+log_k -8.538delta_h 24.760 kcalU+4 + 5 H2O = U(OH)5- + 5 H+log_k -13.147delta_h 27.580 kcal# secondary master species for U(5)U+4 + 2 H2O = UO2+
8、+ 4 H+ + e-log_k -6.432delta_h 31.130 kcal# secondary master species for U(6)U+4 + 2 H2O = UO2+2 + 4 H+ + 2 e-log_k -9.217delta_h 34.430 kcalUO2+2 + H2O = UO2OH+ + H+log_k -5.782delta_h 11.015 kcal2UO2+2 + 2H2O = (UO2)2(OH)2+2 + 2H+log_k -5.626delta_h -36.04 kcal3UO2+2 + 5H2O = (UO2)3(OH)5+ + 5H+log
9、_k -15.641delta_h -44.27 kcalUO2+2 + CO3-2 = UO2CO3log_k 10.064delta_h 0.84 kcalUO2+2 + 2CO3-2 = UO2(CO3)2-2log_k 16.977delta_h 3.48 kcalUO2+2 + 3CO3-2 = UO2(CO3)3-4log_k 21.397delta_h -8.78 kcalPHASESUraniniteUO2 + 4 H+ = U+4 + 2 H2Olog_k -3.490delta_h -18.630 kcalEND程序的的数据库文件phreeqc.dat中不包含铀。这样，当应
10、用这个数据库文件时，输入文件中一定得包括描述热动力学和液相中含铀组分的数据。需要两个关键字来定义铀的形态，即SOLUTION_MASTER_SPECIES和SOLUTION_SPECIES。通过把这两个数据块加到输入文件中，将会在程序运行中确定液相中含铀组分。为把铀稳定地加到列出的元素中，则这些数据块应加入到数据库文件中。这里铀的数据是说明性的，而不是铀物质的完整描述。使用SOLUTION_MASTER_SPECIES输入来定义含铀的主要物质成分是必要的。因为铀是活泼的氧化还原元素，所以定义具有不同化合价的次要含铀物质也是很有必要的。SOLUTION_MASTER_SPECIES(表11)数据
11、块定义了 U+4为主要的含铀物质，同时+4价的铀也是次级主要物质。UO 2 +是化合价为+5 的次级主要含铀物质，UO 2 +2是化合价为+6 的次级主要含铀物质。定义这些液相和其它任何铀络合物的方程必须通过SOLUTION_SPECIES输入来进行。在数据块SOLUTION_SPECIES (表11) 中，主要的和次要的物质均附有注释。首要的主要物质总是以恒等反应（U+4 = U+4）的形式来定义的。次主要物质是在化学反应中仅有的含有电子的液相。另外的氢氧化物和碳酸盐络合物定义为+4和+6 价，无+5价。最后，在PHASES输入中定义一种新的含铀矿物。在物质形成模拟中该物质将会被用来计算饱和
12、指数，在计算机运行中的批反应、运移或是反向模拟中，如果没有重新定义，则不能使用。表 12-例 1 的输出Input file: ex1Output file: ex1.outDatabase file: ./phreeqc.dat-Reading data base.-SOLUTION_MASTER_SPECIESSOLUTION_SPECIESPHASESEXCHANGE_MASTER_SPECIESEXCHANGE_SPECIESSURFACE_MASTER_SPECIESSURFACE_SPECIESRATESEND-Reading input data for simulation 1
13、.-SOLUTION 1 SEAWATER FROM NORDSTROM ET AL. (1979)units ppmpH 8.22pe 8.451density 1.023temp 25.0redox O(0)/O(-2)Ca 412.3Mg 1291.8Na 10768.0K 399.1Fe 0.002Mn 0.0002 peSi 4.28Cl 19353.0Alkalinity 141.682 as HCO3S(6) 2712.0N(5) 0.29 as NO3N(-3) 0.03 as NH4U 3.3 ppb N(5)/N(-3)O(0) 1.0 O2(g) -0.7SOLUTION
14、_MASTER_SPECIESU U+4 0.0 238.0290 238.0290U(4) U+4 0.0 238.0290U(5) UO2+ 0.0 238.0290U(6) UO2+2 0.0 238.0290SOLUTION_SPECIESU+4 = U+4log_k 0.0U+4 + 4 H2O = U(OH)4 + 4 H+log_k -8.538delta_h 24.760 kcalU+4 + 5 H2O = U(OH)5- + 5 H+log_k -13.147delta_h 27.580 kcalU+4 + 2 H2O = UO2+ + 4 H+ + e-log_k -6.4
15、32delta_h 31.130 kcalU+4 + 2 H2O = UO2+2 + 4 H+ + 2 e-log_k -9.217delta_h 34.430 kcalUO2+2 + H2O = UO2OH+ + H+log_k -5.782delta_h 11.015 kcal2UO2+2 + 2H2O = (UO2)2(OH)2+2 + 2H+log_k -5.626delta_h -36.04 kcal3UO2+2 + 5H2O = (UO2)3(OH)5+ + 5H+log_k -15.641delta_h -44.27 kcalUO2+2 + CO3-2 = UO2CO3log_k
16、 10.064delta_h 0.84 kcalUO2+2 + 2CO3-2 = UO2(CO3)2-2log_k 16.977delta_h 3.48 kcalUO2+2 + 3CO3-2 = UO2(CO3)3-4log_k 21.397delta_h -8.78 kcalPHASESUraniniteUO2 + 4 H+ = U+4 + 2 H2Olog_k -3.490delta_h -18.630 kcalEND-TITLE-Example 1.-Add uranium and speciate seawater.-Beginning of initial solution calc
17、ulations.-Initial solution 1. SEAWATER FROM NORDSTROM ET AL. (1979)-Solution composition-Elements Molality MolesAlkalinity 2.406e-03 2.406e-03Ca 1.066e-02 1.066e-02Cl 5.657e-01 5.657e-01Fe 3.711e-08 3.711e-08K 1.058e-02 1.058e-02Mg 5.507e-02 5.507e-02Mn 3.773e-09 3.773e-09N(-3) 1.724e-06 1.724e-06N(
18、5) 4.847e-06 4.847e-06Na 4.854e-01 4.854e-01O(0) 3.746e-04 3.746e-04 Equilibrium with O2(g)S(6) 2.926e-02 2.926e-02Si 7.382e-05 7.382e-05U 1.437e-08 1.437e-08-Description of solution-pH = 8.220 pe = 8.451 Activity of water = 0.981Ionic strength = 6.748e-01Mass of water (kg) = 1.000e+00Total carbon (
19、mol/kg) = 2.180e-03Total CO2 (mol/kg) = 2.180e-03Temperature (deg C) = 25.000Electrical balance (eq) = 7.936e-04Percent error, 100*(Cat-|An|)/(Cat+|An|) = 0.07Iterations = 7Total H = 1.110147e+02Total O = 5.563047e+01-Redox couples-Redox couple pe Eh (volts)N(-3)/N(5) 4.6737 0.2765O(-2)/O(0) 12.3893
20、 0.7329-Distribution of species-Log Log Log Species Molality Activity Molality Activity GammaOH- 2.674e-06 1.629e-06 -5.573 -5.788 -0.215H+ 7.981e-09 6.026e-09 -8.098 -8.220 -0.122H2O 5.551e+01 9.806e-01 -0.009 -0.009 0.000C(4) 2.180e-03HCO3- 1.514e-03 1.023e-03 -2.820 -2.990 -0.170MgHCO3+ 2.195e-04
21、 1.640e-04 -3.658 -3.785 -0.127NaHCO3 1.667e-04 1.948e-04 -3.778 -3.710 0.067MgCO3 8.913e-05 1.041e-04 -4.050 -3.982 0.067NaCO3- 6.718e-05 5.020e-05 -4.173 -4.299 -0.127CaHCO3+ 4.597e-05 3.106e-05 -4.337 -4.508 -0.170CO3-2 3.821e-05 7.959e-06 -4.418 -5.099 -0.681CaCO3 2.725e-05 3.183e-05 -4.565 -4.4
22、97 0.067CO2 1.210e-05 1.413e-05 -4.917 -4.850 0.067UO2(CO3)3-4 1.255e-08 1.183e-10 -7.901 -9.927 -2.025UO2(CO3)2-2 1.814e-09 5.653e-10 -8.741 -9.248 -0.506MnCO3 2.696e-10 3.150e-10 -9.569 -9.502 0.067MnHCO3+ 6.077e-11 4.541e-11 -10.216 -10.343 -0.127UO2CO3 7.429e-12 8.678e-12 -11.129 -11.062 0.067Fe
23、CO3 1.952e-20 2.281e-20 -19.709 -19.642 0.067FeHCO3+ 1.635e-20 1.222e-20 -19.786 -19.913 -0.127Ca 1.066e-02Ca+2 9.504e-03 2.380e-03 -2.022 -2.623 -0.601CaSO4 1.083e-03 1.265e-03 -2.965 -2.898 0.067CaHCO3+ 4.597e-05 3.106e-05 -4.337 -4.508 -0.170CaCO3 2.725e-05 3.183e-05 -4.565 -4.497 0.067CaOH+ 8.60
24、4e-08 6.429e-08 -7.065 -7.192 -0.127Cl 5.657e-01Cl- 5.657e-01 3.528e-01 -0.247 -0.452 -0.205MnCl+ 9.582e-10 7.160e-10 -9.019 -9.145 -0.127MnCl2 9.439e-11 1.103e-10 -10.025 -9.958 0.067MnCl3- 1.434e-11 1.071e-11 -10.844 -10.970 -0.127FeCl+2 9.557e-19 2.978e-19 -18.020 -18.526 -0.506FeCl2+ 6.281e-19 4
25、.693e-19 -18.202 -18.329 -0.127FeCl+ 7.786e-20 5.817e-20 -19.109 -19.235 -0.127FeCl3 1.417e-20 1.656e-20 -19.849 -19.781 0.067Fe(2) 6.909e-19Fe+2 5.205e-19 1.195e-19 -18.284 -18.923 -0.639FeCl+ 7.786e-20 5.817e-20 -19.109 -19.235 -0.127FeSO4 4.845e-20 5.660e-20 -19.315 -19.247 0.067FeCO3 1.952e-20 2
26、.281e-20 -19.709 -19.642 0.067FeHCO3+ 1.635e-20 1.222e-20 -19.786 -19.913 -0.127FeOH+ 8.227e-21 6.147e-21 -20.085 -20.211 -0.127FeHSO4+ 3.000e-27 2.242e-27 -26.523 -26.649 -0.127Fe(3) 3.711e-08Fe(OH)3 2.841e-08 3.318e-08 -7.547 -7.479 0.067Fe(OH)4- 6.591e-09 4.924e-09 -8.181 -8.308 -0.127Fe(OH)2+ 2.
27、118e-09 1.583e-09 -8.674 -8.801 -0.127FeOH+2 9.425e-14 2.937e-14 -13.026 -13.532 -0.506FeSO4+ 1.093e-18 8.167e-19 -17.961 -18.088 -0.127FeCl+2 9.557e-19 2.978e-19 -18.020 -18.526 -0.506FeCl2+ 6.281e-19 4.693e-19 -18.202 -18.329 -0.127Fe+3 3.509e-19 2.796e-20 -18.455 -19.554 -1.099Fe(SO4)2- 6.371e-20
28、 4.760e-20 -19.196 -19.322 -0.127FeCl3 1.417e-20 1.656e-20 -19.849 -19.781 0.067Fe2(OH)2+4 2.462e-24 2.322e-26 -23.609 -25.634 -2.025FeHSO4+2 4.228e-26 1.318e-26 -25.374 -25.880 -0.506Fe3(OH)4+5 1.122e-29 7.679e-33 -28.950 -32.115 -3.165H(0) 0.000e+00H2 0.000e+00 0.000e+00 -44.436 -44.369 0.067K 1.0
29、58e-02K+ 1.041e-02 6.494e-03 -1.982 -2.187 -0.205KSO4- 1.639e-04 1.225e-04 -3.785 -3.912 -0.127KOH 3.137e-09 3.664e-09 -8.504 -8.436 0.067Mg 5.507e-02Mg+2 4.742e-02 1.371e-02 -1.324 -1.863 -0.539MgSO4 7.330e-03 8.562e-03 -2.135 -2.067 0.067MgHCO3+ 2.195e-04 1.640e-04 -3.658 -3.785 -0.127MgCO3 8.913e
30、-05 1.041e-04 -4.050 -3.982 0.067MgOH+ 1.084e-05 8.100e-06 -4.965 -5.092 -0.127Mn(2) 3.773e-09Mn+2 2.171e-09 4.982e-10 -8.663 -9.303 -0.639MnCl+ 9.582e-10 7.160e-10 -9.019 -9.145 -0.127MnCO3 2.696e-10 3.150e-10 -9.569 -9.502 0.067MnSO4 2.021e-10 2.360e-10 -9.695 -9.627 0.067MnCl2 9.439e-11 1.103e-10
31、 -10.025 -9.958 0.067MnHCO3+ 6.077e-11 4.541e-11 -10.216 -10.343 -0.127MnCl3- 1.434e-11 1.071e-11 -10.844 -10.970 -0.127MnOH+ 2.789e-12 2.084e-12 -11.555 -11.681 -0.127Mn(NO3)2 1.375e-20 1.606e-20 -19.862 -19.794 0.067Mn(3) 5.993e-26Mn+3 5.993e-26 4.349e-27 -25.222 -26.362 -1.139N(-3) 1.724e-06NH4+
32、1.648e-06 9.272e-07 -5.783 -6.033 -0.250NH3 7.507e-08 8.769e-08 -7.125 -7.057 0.067N(5) 4.847e-06NO3- 4.847e-06 2.846e-06 -5.315 -5.546 -0.231Mn(NO3)2 1.375e-20 1.606e-20 -19.862 -19.794 0.067Na 4.854e-01Na+ 4.791e-01 3.387e-01 -0.320 -0.470 -0.151NaSO4- 6.053e-03 4.522e-03 -2.218 -2.345 -0.127NaHCO
33、3 1.667e-04 1.948e-04 -3.778 -3.710 0.067NaCO3- 6.718e-05 5.020e-05 -4.173 -4.299 -0.127NaOH 3.117e-07 3.641e-07 -6.506 -6.439 0.067O(0) 3.746e-04O2 1.873e-04 2.188e-04 -3.727 -3.660 0.067S(6) 2.926e-02SO4-2 1.463e-02 2.664e-03 -1.835 -2.574 -0.740MgSO4 7.330e-03 8.562e-03 -2.135 -2.067 0.067NaSO4-
34、6.053e-03 4.522e-03 -2.218 -2.345 -0.127CaSO4 1.083e-03 1.265e-03 -2.965 -2.898 0.067KSO4- 1.639e-04 1.225e-04 -3.785 -3.912 -0.127HSO4- 2.089e-09 1.561e-09 -8.680 -8.807 -0.127MnSO4 2.021e-10 2.360e-10 -9.695 -9.627 0.067FeSO4+ 1.093e-18 8.167e-19 -17.961 -18.088 -0.127Fe(SO4)2- 6.371e-20 4.760e-20
35、 -19.196 -19.322 -0.127FeSO4 4.845e-20 5.660e-20 -19.315 -19.247 0.067FeHSO4+2 4.228e-26 1.318e-26 -25.374 -25.880 -0.506FeHSO4+ 3.000e-27 2.242e-27 -26.523 -26.649 -0.127Si 7.382e-05H4SiO4 7.110e-05 8.306e-05 -4.148 -4.081 0.067H3SiO4- 2.720e-06 2.032e-06 -5.565 -5.692 -0.127H2SiO4-2 7.362e-11 2.29
36、4e-11 -10.133 -10.639 -0.506U(4) 1.040e-21U(OH)5- 1.040e-21 7.773e-22 -20.983 -21.109 -0.127U(OH)4 1.662e-25 1.941e-25 -24.779 -24.712 0.067U+4 0.000e+00 0.000e+00 -46.994 -49.020 -2.025U(5) 1.627e-18UO2+ 1.627e-18 1.216e-18 -17.789 -17.915 -0.127U(6) 1.437e-08UO2(CO3)3-4 1.255e-08 1.183e-10 -7.901
37、-9.927 -2.025UO2(CO3)2-2 1.814e-09 5.653e-10 -8.741 -9.248 -0.506UO2CO3 7.429e-12 8.678e-12 -11.129 -11.062 0.067UO2OH+ 3.386e-14 2.530e-14 -13.470 -13.597 -0.127UO2+2 3.019e-16 9.410e-17 -15.520 -16.026 -0.506(UO2)2(OH)2+2 1.780e-21 5.547e-22 -20.750 -21.256 -0.506(UO2)3(OH)5+ 2.908e-23 2.173e-23 -22.536 -22.663 -0.127-Saturation
