PHREEQC使用手册使用手册使用手册使用手册第一部分第一部分第一部分第一部分关键字关键字关键字关键字SOLUTION这个关键字数据块是用来定义初始溶液的温度和化学组分。所有输入的浓度都将转化成以摩尔为单位,即等于元素的摩尔数、元素的化合价状态和水的质量。若在每一种溶液中离子存在状态的计算是有效的,那每一种溶液对随后的批反应、运移反应或反向模拟计算是可行的。离子存在的能力是用来调整单个元素的浓度以达到电荷的平衡或是与纯相的平衡。数据块的例子数据块的例子数据块的例子数据块的例子Line0:SOLUTION25Testsolutionnumber25Line1:temp25.0Line2:pH7.0chargeLine3:pe4.5Line4:redoxO(-2)/O(0)Line5:unitsppmLine6:density1.02Line7a:Ca80.Line7b:S(6)96.asSO4Line7c:S(-2)1.asSLine7d:N(5)N(3)14.asNLine7e:O(0)8.0Line7f:C61.0asHCO3CO2(g)-3.5Line7g:Fe55.ug/kgsasFeS(6)/S(-2)PyriteLine8a:-isotope13C-12.1.#permilPDBLine8b:-isotope34S15.1.5#permilCDTLine9:-water0.5#kg解释说明解释说明解释说明解释说明Line0:溶液溶液溶液溶液((((SOLUTION))))[序号(number)][描述(scription)]溶液溶液溶液溶液((((SOLUTION))))--是数据块的关键字。[序号(number)--正数是用来指定下列溶液的组分。可以以m到n形式给定的一组数,这里m小于n,并且这两个数是用没有空格间隔的连接号分开的。缺省为1。[描述(scription)]—描述溶液的可选的部分。Line1:温度温度温度温度((((temp))))温度(temperature)温度温度温度温度((((temp))))--这一行中表明的是模拟用到的温度。可选用,temperature,或t[emperature].温度(temperature)--温度(摄氏度),缺省为25摄氏度。Line2:pHpH[(电荷电荷电荷电荷((((charge))))或相的名字(phasename)[饱和指数(saturationindex)])]pH--这一行表明的是溶液的pH值。可选用-ph。pH--pH值,氢离子活度系数的负对数。电荷电荷电荷电荷((((charge))))--表明调整pH值以取得电荷的平衡。如果电荷电荷电荷电荷((((charge))))特定为pH值而指定的,那么它将不能应用于其他的元素。相的名字(phasename)--调整pH值以达到指定相的指定的饱和指数。饱和指数(saturationindex)--pH值将会得到调整以达到特定相中的饱和指数。缺省为0.0。如果第二行没有输入,则缺省的pH值为7.0,同时指定电荷(charge)和相的名字是不允许的。注意,指定相是可行的,但不可能调整pH值以达到指定的饱和指数。Line3:pepe[(电荷电荷电荷电荷((((charge))))或相的名字(phasename)[饱和指数(saturationindex)])]pe--这一行输入pe值。可选用-pe。pe--pe值,电子活度系数的负对数。电荷电荷电荷电荷((((charge))))--(不推荐使用)表明调整pe值以达到电荷的平衡。相的名字(phasename)--在特定的相中,pe值将会得到调整以达到指定的饱和指数。饱和指数(saturationindex)--调整pe值到以达到特定相的指定的饱和指数。缺省为0.0。如果第三行没有输入,则缺省的pe值为7.0,同时指定电荷(charge)和相的名字是不允许的。不推荐调整pe值以达到电荷的平衡。在调整pe值以达到某相指定的饱和指数时,应当是非常小心的。因为通常情况下,这是不可能的。Line4:氧化还原化应氧化还原化应氧化还原化应氧化还原化应((((redox))))氧化还原电对(redoxcouple)氧化还原化应氧化还原化应氧化还原化应氧化还原化应((((redox))))--表明了用来定义计算电势pe的氧化还原电对。这个电势pe用于任何氧化还原元素,通过它,这种氧化还原元素的电势pe需来用来决定元素化合价状态的分配。可选用-r[edox]。氧化还原电对(redoxcouple)--氧化还原电对决定了电势pe。氧化还原电对的确定是以一种元素的两种化合价的状态来指定的,中间以“/”隔开,中间不允许有空格。如果没有包含这样的line4,那么输入的pe值是缺省的。氧化还原电位(-redox)的使用不改变pe输入值。示例数据块用溶解氧的浓度[行7e中由O(0)定义]和从水中形成O2(aq)的氧化还原半反应(缺省数据库中在SOLUTION_SPECIES数据块中定义),来计算缺省的电势pe。Line5:单位单位单位单位((((units))))浓度单位(concentrationunits)单位单位单位单位((((units))))--这一行中表明缺省的浓度单位,可选用,-u[nits]。浓度单位(concentrationunits)--缺省的浓度单位。有三组浓度单位是允许的,浓度(l)每升(“/L”),(2)每千克溶液(“/kgs”),或(3)每千克水(“/kgw”)。溶液中的所有的浓度单位必须在同一个组中,或是以克,或是以摩来表示,并且前缀milli(m)和micro(u)是可接受的。在“每千克溶液”组中的“ppt”每千百万的部分,“ppm”是每亿的部分,和“ppb”都是可以接受的。缺省是mmol/kgw(每千克水的毫摩)。Line6:密度密度密度密度((((density))))密度(density)密度密度密度密度((((density))))--指示在这一行中输入密度。可选用,dens,or-d[ensity]。密度(density)--溶液的密度,kg/L(等于g/cm3),缺省1.0。这个密度仅适用于输入的浓度单位是“每升”。Line7:元素目录(elementlist),浓度(concentration),[单位(units)],([asformula]or[gfwgfw]),[氧化还还电对redoxcouple],[(电荷或相的名字(chargeorphasename)[饱和指数(saturationindex)])]元素目录(elementlist)--元素的名字或化合价的目录,以空格隔开(看line7d)。元素的名字和化合价的状态必须符合于在SOLUTION_MASTER_SPECIES第一栏中的项目。浓度(concentration)--溶液中元素的浓度或是溶液中的元素化合价状态浓度的总和。[单位(units)]—元素的浓度单位(line7g)。如果单位没有被指定,缺省的单位(若line5存在,单位的值是存在的,若不存在line5,则为mmol/kgw)将被设定。asformula–表明了从计算克分子式重量中给定的一种化学分子式,formula。仅在输入浓度是以质量为单位时,克分子式的重量是需要的。在这种溶液中的这种元素的克分子式的重量是用来把质量单位转化成摩尔单位,它不为将来的应用而保存。如果并没有指定一种克分子式,缺省的克分子的定义是在SOLUTION_MASTER_SPECIES中。对碱度而言,公式应当给与克当量。对碳酸钙的碱度,克当量分子式为Ca0.5(CO3)0.5,随着程序的运行,这在phreeqc.dat和wateq4f.dat数据块中是缺省的。gfwgfw--这表明的是将会碰到的克分子量,gfw。当输入的浓度是以质量单位表示的时候,克分子量才是需要的。在这种溶液中的这种元素,特定的克分子量是用来把质量单位转化成摩尔单位,它不为将来的应用来保存。如果没有指定克分子量,缺省的克分子量的定义是在SOLUTION_MASTER_SPECIES中。对碱度而言,克当量是应当使用的,碳酸钙的碱度,克当量大约为50.04g/eq。氧化还还电对redoxcouple—在元素列表中,氧化还原电对是用来定义元素和元素化合价的状态。仅在氧化还原是活泼的或是以下两种情况下,定义氧化还原电对是恰当的。这两种情况是(1)元素的总量是确定的(在元素名字中间没有括号)或是(2)两种或两种以上的化合价的状态被确定(跟在元素名字后面的括号中的化合价的状态被确定)。对于非氧化还原活性的元素或是单个元素的化合价状态定义其氧化还原电对是没有必要的。初始溶液的计算不需要所有氧化还原元素的氧化还原电对的氧化还原平衡。指定一种氧化还原电对(aredoxcouple)将会加强可选择的氧化还原平衡。氧化还原电对的平衡将会与指定的氧化还原元素相平衡。氧化还原元素的定义将会以一种元素以“/”隔开的两种化合价的状态来确定,中间不允许有空格。指定的氧化还原电对与缺省的pe或是缺省的氧化还原电对无关,这个氧化还原电对被用来计算不同化合价状态元素下的溶液中的pe值。如果没有碰到氧化还原电对,则将用到在line4中缺省定义的氧化还原电对,如果没有line4,则将使用pe。电荷电荷电荷电荷charge--表明通过调整这种元素的浓度来达到电荷的平衡,这种元素要具有离子的组分。如果电荷charge是为一种元素所确定,那它不应当指定pH或其他任何的元素。(注释这比通过调整去除所有指定的元素具有更大的电荷不平衡性是可能的,在这种情况下,问题将是不可解决的)。相的名字phasename—调整元素的浓度以达到指定纯相中指定的饱和指数。可以肯定的是与相达到指定的平衡是合理的,应当在相中拥有这种元素的组分。如果已为这种元素指定电荷charge,那么相的名字phasename是不能被应用的。饱和指数(saturationindex)--调整元素的浓度以达到给定纯相中的饱和指数。缺省为0.0。Line8:同位素同位素同位素同位素((((-isotope))))名字(name),值(value),[不确定性限(uncertaintylimit)]同位素同位素同位素同位素((((-isotope))))--在这一行中表明的是一种元素的同位素成分或是碰到的元素的化合价。同位素数据块仅适用于反向模拟反应中。可选用,isotope,或-i[sotope]。名字(name)--同位素的名字。这个名字必须以质量数开头,后面是元素的名字或元素化合价状态的名字,它的定义是在SOLUTION_MASTER_SPECIES中。值(value)--元素或元素化合价态的组分。单位通常是一种比率、千分率或是碳的百分数。不确定性限(uncertaintylimit)--这个不确定性限是用在反向模拟中。这个值在SOLUTION数据块中是可选的,或是应用缺省的不确定性限(INVERSE_MODELING),或是在INVERSE_MODELING数据块中以-isotopes标识符定义的不确定性限Line9:水水水水((((-water))))质量(mass)水水水水((((-water))))--这一行中表明的是水的质量。溶液摩尔数的计算是根据输入的浓度和溶液中水的质量所决定的摩尔数。可选用,water,or-w[ater]。质量(mass)--溶液中水的质量,缺省为1kg。注释在溶液(SOLUTION))))输入行中,输入的顺序是不重要的。在一行中,同时指定“as”和“gfw”是不允许的。在一行中,同时指定“charge”和相的名字也是不允许。指定化合价状态的浓度或多次定义元素的浓度是不允许的。例如,指定总Fe和Fe(+2)的浓度是不允许的,因为铁离子在隐含状态下,被定义了两次。碱度或总碳或者是两者都可以在溶液输入中被确定,如果碱度和总碳都被确定了,那将调整pH来取得指定的碱度。如果碱度的单位是以碳酸钙来表示的,所应用的正确的分子式是“