《phreeqc实例练习》由会员分享，可在线阅读，更多相关《phreeqc实例练习（102页珍藏版）》请在金锄头文库上搜索。
1、wordPHREEQC实例分析例1物种形成分析这个例子计算了海水中矿物质的分布以与一组有关矿物在海水中的饱和程度。为了证明如何在这个模型中应用新的元素，将元素铀添参加由定义的液相模型中wateq.dat是包含于程序分类中的一个数据库文件，它来自于WATEQ4FBall and Nordstrom, 1991，并包含铀。物质形成计算所需要的数据包括温度、Ph、元素的浓度和/或其元素的化合价。海水中的这些数据见表10。这个例子计算中输入的数据组见表11。在模拟中所运用的有关计算的注释包含在TITLE关键字中。SOLUTION数据块定义了海水的成分。注意：元素的化合价用元素化学符号后面圆括号中的数字
2、表示S(6),N(5), N(-3)和O(0)。表10海水的成分未指定浓度时，其浓度的单位为ppm分析的组分PHREEQC符号浓度钙Ca镁Mg钠Na钾K铁Fe.002锰Mn.0002硅石，SiO2Si氯化物Cl碱度，HCO3-Alkalinity硫酸盐，SO42-S(6)硝酸盐，NO3-N(5).29铵，NH4+N(-3).03铀U.0033pH，标准单位pHpe，无单位pe温度，temperature密度，千克/升density用于分配氧化复原元素和计算饱和指数的pe由redox标识符所指定。在这个例子中，用氧化复原电对O(-2)/O(0) 计算的pe值相对应于溶解氧/水，并且这个pe适用于
3、需要pe值的所有的计算。如果redox没有指定，那么缺省的值将会是所输入的pe。缺省的氧化复原标识符可被任何氧化复原元素代替，如输入元素锰时，如此输入的pe被用来表示各种化合价状态的锰；输入铀时，这里是氮/铵电对将会用来计算所形成各种价态铀的pe值。数据组中缺省的单位为ppmunits标识符。这个缺省值可以替换为任何浓度单位，如指定铀的浓度为ppb来代替ppm。因为ppm是一个质量单位，而不是一个摩尔单位，这个程序必须用分子量来将浓度单位转化为摩尔单位。每一种主要物质缺省的分子量在SOLUTION_MASTER_SPECIES输入中指定缺省数据库的值列在表4和附录B中。如果提交的分子量数据不同
4、于其缺省值，必须在输入数据的设置中指定适当的分子量。这可以用gfw标识符来完成，在这里输入真正的分子量，转化硝酸盐的分子量为62.0 g/mol，或是更简便的是以as标识符来完成，在这里输入所使用的化学分子式的单位，正如在这个例子中输入的碱和铵是一样的。注意最后给定的溶解氧O(0)的浓度是1ppm的初始估计值，但它的浓度将会得以调整，直到氧气分压的对数达到。O2(g)的定义是在缺省数据库文件中在PHASES输入(附录B)。当使用相均衡来指定初始浓度正如这个例子中的O(0)，如此仅有一种浓度是得以调整。例如，例如石膏被用来调整钙的浓度，钙的浓度会改变，而硫酸盐的浓度却保持不变。表11 例1的输入
5、数据TITLE Example 1.-Add uranium and speciate seawater.SOLUTION 1 SEAWATER FROM NORDSTROM ET AL. (1979)units ppmredox O(0)/O(-2)Mn 0.0002 peAlkalinity 141.682 as HCO3N(-3) 0.03 as NH4U3.3 ppb N(5)/N(-3)SOLUTION_MASTER_SPECIESSOLUTION_SPECIES#primary master species for U#is also secondary master specie
6、s for U(4)U+4 = U+4U+4 + 4 H2O = U(OH)4 + 4 H+delta_h 24.760 kcalU+4 + 5 H2O = U(OH)5- + 5 H+delta_h 27.580 kcal#secondary master species for U(5)U+4 + 2 H2O = UO2+ + 4 H+ + e-delta_h 31.130 kcal#secondary master species for U(6)U+4 + 2 H2O = UO2+2 + 4 H+ + 2 e-delta_h 34.430 kcalUO2+2 + H2O = UO2OH
7、+ + H+delta_h11.015 kcal2UO2+2 + 2H2O = (UO2)2(OH)2+2 + 2H+delta_h -36.04 kcal3UO2+2 + 5H2O = (UO2)3(OH)5+ + 5H+delta_h -44.27 kcalUO2+2 + CO3-2 = UO2CO3delta_h 0.84 kcalUO2+2 + 2CO3-2 = UO2(CO3)2-2delta_h 3.48 kcalUO2+2 + 3CO3-2 = UO2(CO3)3-4delta_h -8.78 kcalPHASESUraniniteUO2 + 4 H+ = U+4 + 2 H2O
8、delta_h -18.630 kcalEND程序的的数据库文件中不包含铀。这样，当应用这个数据库文件时，输入文件中一定得包括描述热动力学和液相中含铀组分的数据。需要两个关键字来定义铀的形态，即SOLUTION_MASTER_SPECIES和SOLUTION_SPECIES。通过把这两个数据块加到输入文件中，将会在程序运行中确定液相中含铀组分。为把铀稳定地加到列出的元素中，如此这些数据块应参加到数据库文件中。这里铀的数据是说明性的，而不是铀物质的完整描述。使用SOLUTION_MASTER_SPECIES输入来定义含铀的主要物质成分是必要的。因为铀是活泼的氧化复原元素，所以定义具有不同化合价的
9、次要含铀物质也是很有必要的。SOLUTION_MASTER_SPECIES(表11)数据块定义了U+4为主要的含铀物质，同时+4价的铀也是次级主要物质。UO2+是化合价为+5的次级主要含铀物质，UO2+2是化合价为+6的次级主要含铀物质。定义这些液相和其它任何铀络合物的方程必须通过SOLUTION_SPECIES输入来进展。在数据块SOLUTION_SPECIES (表11)中，主要的和次要的物质均附有注释。首要的主要物质总是以恒等反响U+4 = U+4的形式来定义的。次主要物质是在化学反响中仅有的含有电子的液相。另外的氢氧化物和碳酸盐络合物定义为+4和+6价，无+5价。最后，在PHASES输
10、入中定义一种新的含铀矿物。在物质形成模拟中该物质将会被用来计算饱和指数，在计算机运行中的批反响、运移或是反向模拟中，如果没有重新定义，如此不能使用。表12-例1的输出Input file: ex1-Reading data base.-SOLUTION_MASTER_SPECIESSOLUTION_SPECIESPHASESEXCHANGE_MASTER_SPECIESEXCHANGE_SPECIESSURFACE_MASTER_SPECIESSURFACE_SPECIESRATESEND-Reading input data for simulation 1.-SOLUTION 1 SEAW
11、ATER FROM NORDSTROM ET AL. (1979) units ppm redox O(0)/O(-2) Mn 0.0002 pe Alkalinity 141.682 as HCO3 N(5) 0.29 as NO3 N(-3) 0.03 as NH4 U 3.3 ppb N(5)/N(-3)SOLUTION_MASTER_SPECIESSOLUTION_SPECIES U+4 = U+4 U+4 + 4 H2O = U(OH)4 + 4 H+ delta_h 24.760 kcal U+4 + 5 H2O = U(OH)5- + 5 H+ delta_h 27.580 kc
12、al U+4 + 2 H2O = UO2+ + 4 H+ + e- delta_h 31.130 kcal U+4 + 2 H2O = UO2+2 + 4 H+ + 2 e- delta_h 34.430 kcal UO2+2 + H2O = UO2OH+ + H+ delta_h 11.015 kcal 2UO2+2 + 2H2O = (UO2)2(OH)2+2 + 2H+ delta_h -36.04 kcal 3UO2+2 + 5H2O = (UO2)3(OH)5+ + 5H+ delta_h -44.27 kcal UO2+2 + CO3-2 = UO2CO3 delta_h 0.84 kcal UO2+2 + 2CO3-2 = UO2(CO3)2-2 delta_h 3.48 kcal UO2+2 + 3CO3-2 = UO2(CO3)3-4 delta_h -8.78 kcalPHASES Uraninite UO2 + 4 H+ = U+4 + 2 H2O delta_h -18.630 kcalEND-TITLE- Example 1.-Add uranium and speciate seawater.-Beginning of initial solution calculations.
