phreeqc实例练习.docx
《phreeqc实例练习.docx》由会员分享,可在线阅读,更多相关《phreeqc实例练习.docx(144页珍藏版)》请在冰豆网上搜索。
phreeqc实例练习
PHREEQC实例分析
例1——物种形成分析
这个例子计算了海水中矿物质的分布以及一组有关矿物在海水中的饱和程度。
为了证明如何在这个模型中应用新的元素,将元素铀添加入由phreeqc.dat定义的液相模型中[wateq.dat是包含于程序分类中的一个数据库文件,它来自于WATEQ4F(BallandNordstrom,1991),并包含铀]。
物质形成计算所需要的数据包括温度、Ph、元素的浓度和/或其元素的化合价。
海水中的这些数据见表10。
这个例子计算中输入的数据组见表11。
在模拟中所运用的有关计算的注释包含在TITLE关键字中。
SOLUTION数据块定义了海水的成分。
注意:
元素的化合价用元素化学符号后面圆括号中的数字表示[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
硅石,SiO2
Si
4.28
氯化物
Cl
19353.0
碱度,HCO3-
Alkalinity
141.682
硫酸盐,SO42-
S(6)
2712.0
硝酸盐,NO3-
N(5)
.29
铵,NH4+
N(-3)
.03
铀
U
.0033
pH,标准单位
pH
8.22
pe,无单位
pe
8.451
温度,℃
temperature
25.0
密度,千克/升
density
1.023
用于分配氧化还原元素和计算饱和指数的pe由redox标识符所指定。
在这个例子中,用氧化还原电对O(-2)/O(0)计算的pe值相对应于溶解氧/水,并且这个pe适用于需要pe值的所有的计算。
如果redox没有指定,那么缺省的值将会是所输入的pe。
缺省的氧化还原标识符可被任何氧化还原元素代替,如输入元素锰时,则输入的pe被用来表示各种化合价状态的锰;输入铀时,这里是氮/铵电对将会用来计算所形成各种价态铀的pe值。
数据组中缺省的单位为ppm(units标识符)。
这个缺省值可以替换为任何浓度单位,如指定铀的浓度为ppb来代替ppm。
因为ppm是一个质量单位,而不是一个摩尔单位,这个程序必须用分子量来将浓度单位转化为摩尔单位。
每一种主要物质缺省的分子量在SOLUTION_MASTER_SPECIES输入中指定(缺省数据库phreeqc.dat的值列在表4和附录B中)。
如果提交的分子量数据不同于其缺省值,必须在输入数据的设置中指定适当的分子量。
这可以用gfw标识符来完成,在这里输入真正的分子量,转化硝酸盐的分子量为62.0g/mol,或是更简便的是以as标识符来完成,在这里输入所使用的化学分子式的单位,正如在这个例子中输入的碱和铵是一样的。
注意最后给定的溶解氧O(0)的浓度是1ppm的初始估计值,但它的浓度将会得以调整,直到氧气分压的对数达到-0.7。
[O2(g)的定义是在缺省数据库文件中在PHASES输入(附录B)]。
当使用相均衡来指定初始浓度[正如这个例子中的O(0)],则仅有一种浓度是得以调整。
例如,例如石膏被用来调整钙的浓度,钙的浓度会改变,而硫酸盐的浓度却保持不变。
表11例1的输入数据
TITLEExample1.--Adduraniumandspeciateseawater.
SOLUTION1SEAWATERFROMNORDSTROMETAL.(1979)
unitsppm
pH8.22
pe8.451
density1.023
temp25.0
redoxO(0)/O(-2)
Ca412.3
Mg1291.8
Na10768.0
K399.1
Fe0.002
Mn0.0002pe
Si4.28
Cl19353.0
Alkalinity141.682asHCO3
S(6)2712.0
N(5)0.29gfw62.0
N(-3)0.03asNH4
U3.3ppbN(5)/N(-3)
O(0)1.0O2(g)-0.7
SOLUTION_MASTER_SPECIES
UU+40.0238.0290238.0290
U(4)U+40.0238.0290
U(5)UO2+0.0238.0290
U(6)UO2+20.0238.0290
SOLUTION_SPECIES
#primarymasterspeciesforU
#isalsosecondarymasterspeciesforU(4)
U+4=U+4
log_k0.0
U+4+4H2O=U(OH)4+4H+
log_k-8.538
delta_h24.760kcal
U+4+5H2O=U(OH)5-+5H+
log_k-13.147
delta_h27.580kcal
#secondarymasterspeciesforU(5)
U+4+2H2O=UO2++4H++e-
log_k-6.432
delta_h31.130kcal
#secondarymasterspeciesforU(6)
U+4+2H2O=UO2+2+4H++2e-
log_k-9.217
delta_h34.430kcal
UO2+2+H2O=UO2OH++H+
log_k-5.782
delta_h11.015kcal
2UO2+2+2H2O=(UO2)2(OH)2+2+2H+
log_k-5.626
delta_h-36.04kcal
3UO2+2+5H2O=(UO2)3(OH)5++5H+
log_k-15.641
delta_h-44.27kcal
UO2+2+CO3-2=UO2CO3
log_k10.064
delta_h0.84kcal
UO2+2+2CO3-2=UO2(CO3)2-2
log_k16.977
delta_h3.48kcal
UO2+2+3CO3-2=UO2(CO3)3-4
log_k21.397
delta_h-8.78kcal
PHASES
Uraninite
UO2+4H+=U+4+2H2O
log_k-3.490
delta_h-18.630kcal
END
程序的的数据库文件phreeqc.dat中不包含铀。
这样,当应用这个数据库文件时,输入文件中一定得包括描述热动力学和液相中含铀组分的数据。
需要两个关键字来定义铀的形态,即SOLUTION_MASTER_SPECIES和SOLUTION_SPECIES。
通过把这两个数据块加到输入文件中,将会在程序运行中确定液相中含铀组分。
为把铀稳定地加到列出的元素中,则这些数据块应加入到数据库文件中。
这里铀的数据是说明性的,而不是铀物质的完整描述。
使用SOLUTION_MASTER_SPECIES输入来定义含铀的主要物质成分是必要的。
因为铀是活泼的氧化还原元素,所以定义具有不同化合价的次要含铀物质也是很有必要的。
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输入中定义一种新的含铀矿物。
在物质形成模拟中该物质将会被用来计算饱和指数,在计算机运行中的批反应、运移或是反向模拟中,如果没有重新定义,则不能使用。
表12--例1的输出
Inputfile:
ex1
Outputfile:
ex1.out
Databasefile:
../phreeqc.dat
------------------
Readingdatabase.
------------------
SOLUTION_MASTER_SPECIES
SOLUTION_SPECIES
PHASES
EXCHANGE_MASTER_SPECIES
EXCHANGE_SPECIES
SURFACE_MASTER_SPECIES
SURFACE_SPECIES
RATES
END
------------------------------------
Readinginputdataforsimulation1.
------------------------------------
SOLUTION1SEAWATERFROMNORDSTROMETAL.(1979)
unitsppm
pH8.22
pe8.451
density1.023
temp25.0
redoxO(0)/O(-2)
Ca412.3
Mg1291.8
Na10768.0
K399.1
Fe0.002
Mn0.0002pe
Si4.28
Cl19353.0
Alkalinity141.682asHCO3
S(6)2712.0
N(5)0.29asNO3
N(-3)0.03asNH4
U3.3ppbN(5)/N(-3)
O(0)1.0O2(g)-0.7
SOLUTION_MASTER_SPECIES
UU+40.0238.0290238.0290
U(4)U+40.0238.0290
U(5)UO2+0.0238.0290
U(6)UO2+20.0238.0290
SOLUTION_SPECIES
U+4=U+4
log_k0.0
U+4+4H2O=U(OH)4+4H+
log_k-8.538
delta_h24.760kcal
U+4+5H2O=U(OH)5-+5H+
log_k-13.147
delta_h27.580kcal
U+4+2H2O=UO2++4H++e-
log_k-6.432
delta_h31.130kcal
U+4+2H2O=UO2+2+4H++2e-
log_k-9.217
delta_h34.430kcal
UO2+2+H2O=UO2OH++H+
log_k-5.782
delta_h11.015kcal
2UO2+2+2H2O=(UO2)2(OH)2+2+2H+
log_k-5.626
delta_h-36.04kcal
3UO2+2+5H2O=(UO2)3(OH)5++5H+
log_k-15.641
delta_h-44.27kcal
UO2+2+CO3-2=UO2CO3
log_k10.064
delta_h0.84kcal
UO2+2+2CO3-2=UO2(CO3)2-2
log_k16.977
delta_h3.48kcal
UO2+2+3CO3-2=UO2(CO3)3-4
log_k21.397
delta_h-8.78kcal
PHASES
Uraninite
UO2+4H+=U+4+2H2O
log_k-3.490
delta_h-18.630kcal
END
-----
TITLE
-----
Example1.--Adduraniumandspeciateseawater.
-------------------------------------------
Beginningofinitialsolutioncalculations.
-------------------------------------------
Initialsolution1.SEAWATERFROMNORDSTROMETAL.(1979)
-----------------------------Solutioncomposition------------------------------
ElementsMolalityMoles
Alkalinity2.406e-032.406e-03
Ca1.066e-021.066e-02
Cl5.657e-015.657e-01
Fe3.711e-083.711e-08
K1.058e-021.058e-02
Mg5.507e-025.507e-02
Mn3.773e-093.773e-09
N(-3)1.724e-061.724e-06
N(5)4.847e-064.847e-06
Na4.854e-014.854e-01
O(0)3.746e-043.746e-04EquilibriumwithO2(g)
S(6)2.926e-022.926e-02
Si7.382e-057.382e-05
U1.437e-081.437e-08
----------------------------Descriptionofsolution----------------------------
pH=8.220
pe=8.451
Activityofwater=0.981
Ionicstrength=6.748e-01
Massofwater(kg)=1.000e+00
Totalcarbon(mol/kg)=2.180e-03
TotalCO2(mol/kg)=2.180e-03
Temperature(degC)=25.000
Electricalbalance(eq)=7.936e-04
Percenterror,100*(Cat-|An|)/(Cat+|An|)=0.07
Iterations=7
TotalH=1.110147e+02
TotalO=5.563047e+01
---------------------------------Redoxcouples---------------------------------
RedoxcouplepeEh(volts)
N(-3)/N(5)4.67370.2765
O(-2)/O(0)12.38930.7329
----------------------------Distributionofspecies----------------------------
LogLogLog
SpeciesMolalityActivityMolalityActivityGamma
OH-2.674e-061.629e-06-5.573-5.788-0.215
H+7.981e-096.026e-09-8.098-8.220-0.122
H2O5.551e+019.806e-01-0.009-0.0090.000
C(4)2.180e-03
HCO3-1.514e-031.023e-03-2.820-2.990-0.170
MgHCO3+2.195e-041.640e-04-3.658-3.785-0.127
NaHCO31.667e-041.948e-04-3.778-3.7100.067
MgCO38.913e-051.041e-04-4.050-3.9820.067
NaCO3-6.718e-055.020e-05-4.173-4.299-0.127
CaHCO3+4.597e-053.106e-05-4.337-4.508-0.170
CO3-23.821e-057.959e-06-4.418-5.099-0.681
CaCO32.725e-053.183e-05-4.565-4.4970.067
CO21.210e-051.413e-05-4.917-4.8500.067
UO2(CO3)3-41.255e-081.183e-10-7.901-9.927-2.025
UO2(CO3)2-21.814e-095.653e-10-8.741-9.248-0.506
MnCO32.696e-103.150e-10-9.569-9.5020.067
MnHCO3+6.077e-114.541e-11-10.216-10.343-0.127
UO2CO37.429e-128.678e-12-11.129-11.0620.067
FeCO31.952e-202.281e-20-19.709-19.6420.067
FeHCO3+1.635e-201.222e-20-19.786-19.913-0.127
Ca1.066e-02
Ca+29.504e-032.380e-03-2.022-2.623-0.601
CaSO41.083e-031.265e-03-2.965-2.8980.067
CaHCO3+4.597e-053.106e-05-4.337-4.508-0.170
CaCO32.725e-053.183e-05-4.565-4.4970.067
CaOH+8.604e-086.429e-08-7.065-7.192-0.127
Cl5.657e-01
Cl-5.657e-013.528e-01-0.247-0.452-0.205
MnCl+9.582e-107.160e-10-9.019-9.145-0.127
MnCl29.439e-111.103e-10-10.025-9.9580.067
MnCl3-1.434e-111.071e-11-10.844-10.970-0.127
FeCl+29.557e-192.978e-19-18.020-18.526-0.506
FeCl2+6.281e-194.693e-19-18.202-18.329-0.127
FeCl+7.786e-205.817e-20-19.109-19.235-0.127
FeCl31.417e-201.656e-20-19.849-19.7810.067
Fe
(2)6.909e-19
Fe+25.205e-191.195e-19-18.284-18.923-0.639
FeCl+7.786e-205.817e-20-19.109-19.235-0.127
FeSO44.845e-205.660e-20-19.315-19.2470.067
FeCO31.952e-202.281e-20-19.709-19.6420.067
FeHCO3+1.635e-201.222e-20-19.786-19.913-0.127
FeOH+8.227e-216.147e-21-20.085-20.211-0.127
FeHSO4+3.000e-272.242e-27-26.523-26.649-0.127
Fe(3)3.711e-08
Fe(OH)32.841e-083.318e-08-7.547-7.4790.067
Fe(OH)4-6.591e-094.924e-09-8.181-8.308-0.127
Fe(OH)2+2.118e-091.583e-09-8.674-8.801-0.127
FeOH+29.425e-142.937e-14-13.026-13.532-0.506
FeSO4+1.093e-188.167e-19-17.961-18.088-0.127
FeCl+29.557e-192.978e-19-18.020-18.526-0.506
FeCl2+6.281e-194.693e-19-18.202-18.329-0.127
Fe+33.509e-192.796e-20-18.455-19.554-1.099
Fe(SO4)2-6.371e-204.760e-20-19.196-19.322-0.127
FeCl31.417e-201.656e-20-19.849-19.7810.067
Fe2(OH)2+42.462e-242.322e-26-23.609-25.634-2.025
FeHSO4+24.228e-261.318e-26-25.374-25.880-0.506
Fe3(OH)4+51.122e-297.679e-33-28.950-32.115-3.165
H(0)0.000e+00
H20.000e+000.000e+00-44.436-44.3690.067
K1.058e-02
K+1.041e-026.494e-03-1.982-2.187-0.205
KSO4-1.639e-041.225e-04-3.785-3.912-0.127
KOH3.137e-093.664e-09-8.504-8.4360.067
Mg5.507e-02
Mg+24.742e-021.371e-02-1.324-1.863-0.539
MgSO47.330e-038.562e-03-2.135-2.067