《phreeqc实例练习(1)》由会员分享，可在线阅读，更多相关《phreeqc实例练习(1)（111页珍藏版）》请在金锄头文库上搜索。
1、PHREQC实例分析例-物种形成分析这个例子计算了海水中矿物质的分布以及一组有关矿物在海水中的饱和程度。为了证明如何在这个模型中应用新的元素，将元素铀添加入由phreeqc。dat定义的液相模型中ateq.da是包含于程序分类中的一个数据库文件,它来自于WATEQ4（BallndNodstrom,191），并包含铀。物质形成计算所需要的数据包括温度、元素的浓度和/或其元素的化合价。海水中的这些数据见表。这个例子计算中输入的数据组见表1。在模拟中所运用的有关计算的注释包含在TITLE关键字中。SOLTION数据块定义了海水的成分.注意：元素的化合价用元素化学符号后面圆括号中的数字表示（6）， N
2、（5）, （3）和O(0）.表10海水的成分未指定浓度时，其浓度的单位为ppm分析的组分PHREEQC符号浓度钙a412.3镁Mg1291。8钠a107。0钾K391铁F。锰。00硅石,SiO2Si4.28氯化物Cl19353碱度，HO3Akalinity141.682硫酸盐，SO42(6）272。0硝酸盐，NO3N（5)29铵，NH4（3）。03铀.033p,标准单位H。2pe，无单位pe8。451温度，tmprature5。0密度,千克/升ity1。23用于分配氧化还原元素和计算饱和指数的p由redo标识符所指定。在这个例子中，用氧化还原电对（-2)/O（）计算的pe值相对应于溶解氧/水,
3、并且这个e适用于需要pe值的所有的计算。如果eox没有指定，那么缺省的值将会是所输入的p。缺省的氧化还原标识符可被任何氧化还原元素代替，如输入元素锰时,则输入的e被用来表示各种化合价状态的锰；输入铀时,这里是氮/铵电对将会用来计算所形成各种价态铀的pe值.数据组中缺省的单位为pm（unis标识符).这个缺省值可以替换为任何浓度单位，如指定铀的浓度为来代替pm。因为pm是一个质量单位，而不是一个摩尔单位，这个程序必须用分子量来将浓度单位转化为摩尔单位。每一种主要物质缺省的分子量在OUTION_MAE_PCES输入中指定(缺省数据库hre。dat的值列在表4和附录B中）。如果提交的分子量数据不同于
4、其缺省值,必须在输入数据的设置中指定适当的分子量。这可以用gfw标识符来完成，在这里输入真正的分子量，转化硝酸盐的分子量为62 gmol，或是更简便的是以s标识符来完成，在这里输入所使用的化学分子式的单位，正如在这个例子中输入的碱和铵是一样的。注意最后给定的溶解氧O(0）的浓度是ppm的初始估计值,但它的浓度将会得以调整，直到氧气分压的对数达到-0。O2(）的定义是在缺省数据库文件中在PHSS输入（附录）.当使用相均衡来指定初始浓度正如这个例子中的O()，则仅有一种浓度是得以调整。例如,例如石膏被用来调整钙的浓度，钙的浓度会改变，而硫酸盐的浓度却保持不变。表1 例1的输入数据TITLEExam
5、le .dd uranium anseciae sawatSOLUTIN 1 ATER RMNRSTOM ET L.(979）uni pH 2p 851dnity 。023temp 20redo O（0）O（2）a 12。Mg 1291。a 10768K39.1e000M 00002 pi 48Cl353。Alaliniy 11。682 as O3S（6) 2120N（） 0。29gfw 2。0N(-）.3 as N4U3.3ppb N（5)/N(-3)（0） 1 (g) 0.OLUTINMTER_SECEUU+40. 2380290 2380290(4）U 0.0。0290U(5) UO2
6、.0 238。0290（6） UO2+0.0 238。0290SLUTION_PECI pmay master peiesfor U#i lsoseconary sterpeies fr U（4）U+4 = U4log_k 0.0+ 4 H2O = (O） 4 H+log_k -.538delta_h24.76 kclU+4 + 52= U(OH)5 5 +log_ -13147detah 7。580 kc ecoa aster sece or U(5）U+4+ HO = + 4 H+ -lg_k。4dlta_h31.130 kcl ecodary astr seies for U（6)+4+
7、2 H2 = O22+ 4 H+ elg_9.217elh3430 kalO2+2 H2O = O2OH+ + H+lg_k58delta_h11.0 ka2 + 2H2= （U)2（OH）+ 2Hog_k -6ta_ -36。0 caUO2+ + H2O = （O)3（OH）+ + 5Hlg_-1.641elt_h -4.2alU2+2 + CO32 = UOC3ogk 1006dea_.84 kclU2+2 2CO3-2=UO2（CO3)2-2lo_k 16。977delta_ 3。4 kUO2+ + 3C3-2= UO2(O）3-log_k21.39dlh 8.7 kcalPAESUan
8、eUO+ H+ =U4 + 2H2Og -90deta_h 830 kcalEND程序的的数据库文件preeqcat中不包含铀。这样，当应用这个数据库文件时，输入文件中一定得包括描述热动力学和液相中含铀组分的数据.需要两个关键字来定义铀的形态，即SLINMASTER_CIE和SLTON_PEIES。通过把这两个数据块加到输入文件中,将会在程序运行中确定液相中含铀组分。为把铀稳定地加到列出的元素中，则这些数据块应加入到数据库文件中。这里铀的数据是说明性的，而不是铀物质的完整描述。使用OLUTION_ASTRSPECIS输入来定义含铀的主要物质成分是必要的。因为铀是活泼的氧化还原元素，所以定义具有
9、不同化合价的次要含铀物质也是很有必要的.SLUTION_MSER_SPCIES（表11）数据块定义了+为主要的含铀物质，同时+4价的铀也是次级主要物质.UO2 +是化合价为+5的次级主要含铀物质，UO2 +2是化合价为+6的次级主要含铀物质。定义这些液相和其它任何铀络合物的方程必须通过SOLTIN_PECIES输入来进行。在数据块SLTO_SPCIE（表11)中，主要的和次要的物质均附有注释。首要的主要物质总是以恒等反应(+ U+4)的形式来定义的。次主要物质是在化学反应中仅有的含有电子的液相。另外的氢氧化物和碳酸盐络合物定义为+和+价，无+5价.最后，在HASES输入中定义一种新的含铀矿物。
10、在物质形成模拟中该物质将会被用来计算饱和指数,在计算机运行中的批反应、运移或是反向模拟中，如果没有重新定义，则不能使用.表12-例1的输出Inut ile： ex1Outut fle： x。outDatabsfi： ./hreeqc。dat-Readigdata base.-SOLTONMSTERSISSOUTINSECESHASESEXCANGE_MSTERSEIESEXCHANGE_SPECIESSRFACE_MSTER_SPEISURAEPECIERATEEN-eadig iput data fr siulatio-SLTION 1 SAWTERFROM RSTROM ET L. （197） units m pH 822 pe 8451 densit 10 temp 2.0 redox O（0)/O（2） C 412.3 Mg 12918 Na 1070 399. e 000 Mn 0。002 p S 42 l 19353.0 Allity 141。68asHC3 S（) 272。0 （5） 0.29 asNO3 N（-3） 0.03 as NH4 U . ppb N(5）/N(3） （0) 10 O（g） 0.7SOLUTION_MSTER_PEIS U +4 0. 238020 30290 U（4） U+ .0 2380290 (5)
