最新phreeqc实例练习.docx

1、最新phreeqc实例练习phreeqc实例练习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符号浓度钙Ca412.3镁Mg1291.8钠Na10768.0钾K399.1铁Fe.002锰Mn.0002硅石，SiO2Si4.28氯化物Cl19353.0碱度，HCO3-Alkalinity141.682硫酸盐，SO42-S(6)2712.0硝酸盐，NO3-N(5).29铵，NH4+N(-3).03铀U.0033pH，标准单位pH8.22pe，无单位pe8.451温度，temperature25.0密度，

3、千克/升density1.023用于分配氧化还原元素和计算饱和指数的pe由redox标识符所指定。在这个例子中，用氧化还原电对O(-2)/O(0) 计算的pe值相对应于溶解氧/水，并且这个pe适用于需要pe值的所有的计算。如果redox没有指定，那么缺省的值将会是所输入的pe。缺省的氧化还原标识符可被任何氧化还原元素代替，如输入元素锰时，则输入的pe被用来表示各种化合价状态的锰；输入铀时，这里是氮/铵电对将会用来计算所形成各种价态铀的pe值。数据组中缺省的单位为ppm（units标识符）。这个缺省值可以替换为任何浓度单位，如指定铀的浓度为ppb来代替ppm。因为ppm是一个质量单位，而不是一个

4、摩尔单位，这个程序必须用分子量来将浓度单位转化为摩尔单位。每一种主要物质缺省的分子量在SOLUTION_MASTER_SPECIES输入中指定（缺省数据库phreeqc.dat的值列在表4和附录B中）。如果提交的分子量数据不同于其缺省值，必须在输入数据的设置中指定适当的分子量。这可以用gfw标识符来完成，在这里输入真正的分子量，转化硝酸盐的分子量为62.0 g/mol，或是更简便的是以as标识符来完成，在这里输入所使用的化学分子式的单位，正如在这个例子中输入的碱和铵是一样的。注意最后给定的溶解氧O(0)的浓度是1ppm的初始估计值，但它的浓度将会得以调整，直到氧气分压的对数达到-0.7。O2(

5、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.023temp 25.0redox O(0)/O(-2)Ca 412.3Mg 1291.8Na 1

6、0768.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(6) UO2+2 0.0 238.0290SOLUTION_SPECIES# prima

7、ry 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+ + 4 H+ + e-log_k -6.432delta_h 31.130 kcal#

8、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_k -15.641delta_h -44.27 kcalUO2+2 + CO3-2 =

9、 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)数据块定义了U+4为主要的含铀物质，同时+4价的铀也是次级主要物质。UO2 +是化合价为+5

11、的次级主要含铀物质，UO2 +2是化合价为+6的次级主要含铀物质。定义这些液相和其它任何铀络合物的方程必须通过SOLUTION_SPECIES输入来进行。在数据块SOLUTION_SPECIES (表11)中，主要的和次要的物质均附有注释。首要的主要物质总是以恒等反应（U+4 = U+4）的形式来定义的。次主要物质是在化学反应中仅有的含有电子的液相。另外的氢氧化物和碳酸盐络合物定义为+4和+6价，无+5价。最后，在PHASES输入中定义一种新的含铀矿物。在物质形成模拟中该物质将会被用来计算饱和指数，在计算机运行中的批反应、运移或是反向模拟中，如果没有重新定义，则不能使用。表12-例1的输出In

12、put file: ex1Output file: ex1.outDatabase file: ./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.- SOLUTION 1 SEAWATER FROM NORDSTROM ET AL.

13、(1979) units ppm pH 8.22 pe 8.451 density 1.023 temp 25.0 redox O(0)/O(-2) Ca 412.3 Mg 1291.8 Na 10768.0 K 399.1 Fe 0.002 Mn 0.0002 pe Si 4.28 Cl 19353.0 Alkalinity 141.682 as HCO3 S(6) 2712.0 N(5) 0.29 as NO3 N(-3) 0.03 as NH4 U 3.3 ppb N(5)/N(-3) O(0) 1.0 O2(g) -0.7 SOLUTION_MASTER_SPECIES U U+4 0

14、.0 238.0290 238.0290 U(4) U+4 0.0 238.0290 U(5) UO2+ 0.0 238.0290 U(6) UO2+2 0.0 238.0290 SOLUTION_SPECIES U+4 = U+4 log_k 0.0 U+4 + 4 H2O = U(OH)4 + 4 H+ log_k -8.538 delta_h 24.760 kcal U+4 + 5 H2O = U(OH)5- + 5 H+ log_k -13.147 delta_h 27.580 kcal U+4 + 2 H2O = UO2+ + 4 H+ + e- log_k -6.432 delta

15、_h 31.130 kcal U+4 + 2 H2O = UO2+2 + 4 H+ + 2 e- log_k -9.217 delta_h 34.430 kcal UO2+2 + H2O = UO2OH+ + H+ log_k -5.782 delta_h 11.015 kcal 2UO2+2 + 2H2O = (UO2)2(OH)2+2 + 2H+ log_k -5.626 delta_h -36.04 kcal 3UO2+2 + 5H2O = (UO2)3(OH)5+ + 5H+ log_k -15.641 delta_h -44.27 kcal UO2+2 + CO3-2 = UO2CO

16、3 log_k 10.064 delta_h 0.84 kcal UO2+2 + 2CO3-2 = UO2(CO3)2-2 log_k 16.977 delta_h 3.48 kcal UO2+2 + 3CO3-2 = UO2(CO3)3-4 log_k 21.397 delta_h -8.78 kcal PHASES Uraninite UO2 + 4 H+ = U+4 + 2 H2O log_k -3.490 delta_h -18.630 kcal END-TITLE- Example 1.-Add uranium and speciate seawater.-Beginning of

17、initial solution calculations.-Initial solution 1. SEAWATER FROM NORDSTROM ET AL. (1979)-Solution composition- Elements Molality Moles Alkalinity 2.406e-03 2.406e-03 Ca 1.066e-02 1.066e-02 Cl 5.657e-01 5.657e-01 Fe 3.711e-08 3.711e-08 K 1.058e-02 1.058e-02 Mg 5.507e-02 5.507e-02 Mn 3.773e-09 3.773e-

18、09 N(-3) 1.724e-06 1.724e-06 N(5) 4.847e-06 4.847e-06 Na 4.854e-01 4.854e-01 O(0) 3.746e-04 3.746e-04 Equilibrium with O2(g) S(6) 2.926e-02 2.926e-02 Si 7.382e-05 7.382e-05 U 1.437e-08 1.437e-08-Description of solution- pH = 8.220 pe = 8.451 Activity of water = 0.981 Ionic strength = 6.748e-01 Mass

19、of water (kg) = 1.000e+00 Total carbon (mol/kg) = 2.180e-03 Total CO2 (mol/kg) = 2.180e-03 Temperature (deg C) = 25.000 Electrical balance (eq) = 7.936e-04 Percent error, 100*(Cat-|An|)/(Cat+|An|) = 0.07 Iterations = 7 Total H = 1.110147e+02 Total O = 5.563047e+01-Redox couples- Redox couple pe Eh (

20、volts) N(-3)/N(5) 4.6737 0.2765 O(-2)/O(0) 12.3893 0.7329-Distribution of species- Log Log Log Species Molality Activity Molality Activity Gamma OH- 2.674e-06 1.629e-06 -5.573 -5.788 -0.215 H+ 7.981e-09 6.026e-09 -8.098 -8.220 -0.122 H2O 5.551e+01 9.806e-01 -0.009 -0.009 0.000C(4) 2.180e-03 HCO3- 1.

21、514e-03 1.023e-03 -2.820 -2.990 -0.170 MgHCO3+ 2.195e-04 1.640e-04 -3.658 -3.785 -0.127 NaHCO3 1.667e-04 1.948e-04 -3.778 -3.710 0.067 MgCO3 8.913e-05 1.041e-04 -4.050 -3.982 0.067 NaCO3- 6.718e-05 5.020e-05 -4.173 -4.299 -0.127 CaHCO3+ 4.597e-05 3.106e-05 -4.337 -4.508 -0.170 CO3-2 3.821e-05 7.959e

22、-06 -4.418 -5.099 -0.681 CaCO3 2.725e-05 3.183e-05 -4.565 -4.497 0.067 CO2 1.210e-05 1.413e-05 -4.917 -4.850 0.067 UO2(CO3)3-4 1.255e-08 1.183e-10 -7.901 -9.927 -2.025 UO2(CO3)2-2 1.814e-09 5.653e-10 -8.741 -9.248 -0.506 MnCO3 2.696e-10 3.150e-10 -9.569 -9.502 0.067 MnHCO3+ 6.077e-11 4.541e-11 -10.2

23、16 -10.343 -0.127 UO2CO3 7.429e-12 8.678e-12 -11.129 -11.062 0.067 FeCO3 1.952e-20 2.281e-20 -19.709 -19.642 0.067 FeHCO3+ 1.635e-20 1.222e-20 -19.786 -19.913 -0.127Ca 1.066e-02 Ca+2 9.504e-03 2.380e-03 -2.022 -2.623 -0.601 CaSO4 1.083e-03 1.265e-03 -2.965 -2.898 0.067 CaHCO3+ 4.597e-05 3.106e-05 -4

24、.337 -4.508 -0.170 CaCO3 2.725e-05 3.183e-05 -4.565 -4.497 0.067 CaOH+ 8.604e-08 6.429e-08 -7.065 -7.192 -0.127Cl 5.657e-01 Cl- 5.657e-01 3.528e-01 -0.247 -0.452 -0.205 MnCl+ 9.582e-10 7.160e-10 -9.019 -9.145 -0.127 MnCl2 9.439e-11 1.103e-10 -10.025 -9.958 0.067 MnCl3- 1.434e-11 1.071e-11 -10.844 -1

25、0.970 -0.127 FeCl+2 9.557e-19 2.978e-19 -18.020 -18.526 -0.506 FeCl2+ 6.281e-19 4.693e-19 -18.202 -18.329 -0.127 FeCl+ 7.786e-20 5.817e-20 -19.109 -19.235 -0.127 FeCl3 1.417e-20 1.656e-20 -19.849 -19.781 0.067Fe(2) 6.909e-19 Fe+2 5.205e-19 1.195e-19 -18.284 -18.923 -0.639 FeCl+ 7.786e-20 5.817e-20 -

26、19.109 -19.235 -0.127 FeSO4 4.845e-20 5.660e-20 -19.315 -19.247 0.067 FeCO3 1.952e-20 2.281e-20 -19.709 -19.642 0.067 FeHCO3+ 1.635e-20 1.222e-20 -19.786 -19.913 -0.127 FeOH+ 8.227e-21 6.147e-21 -20.085 -20.211 -0.127 FeHSO4+ 3.000e-27 2.242e-27 -26.523 -26.649 -0.127Fe(3) 3.711e-08 Fe(OH)3 2.841e-0

27、8 3.318e-08 -7.547 -7.479 0.067 Fe(OH)4- 6.591e-09 4.924e-09 -8.181 -8.308 -0.127 Fe(OH)2+ 2.118e-09 1.583e-09 -8.674 -8.801 -0.127 FeOH+2 9.425e-14 2.937e-14 -13.026 -13.532 -0.506 FeSO4+ 1.093e-18 8.167e-19 -17.961 -18.088 -0.127 FeCl+2 9.557e-19 2.978e-19 -18.020 -18.526 -0.506 FeCl2+ 6.281e-19 4

28、.693e-19 -18.202 -18.329 -0.127 Fe+3 3.509e-19 2.796e-20 -18.455 -19.554 -1.099 Fe(SO4)2- 6.371e-20 4.760e-20 -19.196 -19.322 -0.127 FeCl3 1.417e-20 1.656e-20 -19.849 -19.781 0.067 Fe2(OH)2+4 2.462e-24 2.322e-26 -23.609 -25.634 -2.025 FeHSO4+2 4.228e-26 1.318e-26 -25.374 -25.880 -0.506 Fe3(OH)4+5 1.

29、122e-29 7.679e-33 -28.950 -32.115 -3.165H(0) 0.000e+00 H2 0.000e+00 0.000e+00 -44.436 -44.369 0.067K 1.058e-02 K+ 1.041e-02 6.494e-03 -1.982 -2.187 -0.205 KSO4- 1.639e-04 1.225e-04 -3.785 -3.912 -0.127 KOH 3.137e-09 3.664e-09 -8.504 -8.436 0.067Mg 5.507e-02 Mg+2 4.742e-02 1.371e-02 -1.324 -1.863 -0.539 MgSO4 7.330e-03 8.562