PHREEQC软件简介地下水污染与防治

PHREEQC软件简介PHREEQC的开展过程PHREEQC之前的版本是PHREEQE。PHREEQE是1980年美国地质调查局的Plummer和Parkhurst等开发出来的，用于地球化学模拟达15年之久。PHREEQC于1995年正式推出，它是在PHREEQE的源程序根底上用C语言重写而成的，消除了PHREEQE的缺陷与局限性。PHREEQC开展至今已近10年，其版本不断更新，参加了很多新内容，已开展到版本2中的2.11子版本。由于其现有的模拟能力在同类软件中最强大，所以在国际上被广泛使用。PHREEQC的功能PHREEQC以离子联系的水化学模型为根底，可以推算：〔1〕生成物和饱和系数；〔2〕涉及到可逆反响以及不可逆反响的批反响和一维〔1D〕运移计算；其中，涉及到的可逆反响包括水、矿物/无机溶液、气体、固体溶液、外表络合、离子交换平衡；不可逆反响包括指定成分摩尔转换、动态控制反响、溶液混合和温度变化；〔3〕逆向建模，可以发现一系列矿物和气体克分子在水中以不确定组成的转移。PHREEQC第二版的新特点和第一版相比，PHREEQC第二版的新特点如下：具有在一维运移计算中模拟弥散〔或扩散〕和滞流区的能力，使用用户确定的速率表达式模拟分子反响，模拟标准的多种成分或非标准的两种成分的固体溶液的沉淀和溶解，模拟定体积气相和定压力气相，考虑外表系数或交换位置随着无机物的溶解和沉淀或者分子反响的变化而变化，自动采用多套收敛参数，打印用户指定量到原始输出文件和〔或〕适合输入到扩展表格的文件上，以一种与扩展表程序更兼容的形式确定溶解成分等。二、PHREEQC的局限性1.水溶液模型PHREEQC用相关的离子和Debye-Hückel方程计算非理想溶液，这种类型的水溶液模型在较低的离子强度下比较适宜，但在高离子强度下〔海水范围或更高〕那么会失真。另外，数据库缺乏内部的连续性。因为随程序发布的两个数据库phreeqc.dat和wareq.dat中的lgKs’和反响焓来自于各种文献，未曾系统地确定溶液模型中单个lgK，也未曾确定当前数据库中水溶液模型是否与原始的实验数据相适应。使用时，随程序提供的数据文件可首先考虑，但使用哪种水溶液组分和热力学数据那么由用户自己确定。2.离子交换模型离子交换模型假设交换组分的热力学活度等于它的当量分数，或者用当量分数乘以Debye-Hückel活度系数，以此定义交换组分的活度〔Appelo,1994〕，更加复杂的交换模型还不能实现。在许多野外研究中，离子交换模拟需要用研究区的介质进行离子交换试验，通过实验获得的数据来确定适宜的模型。3.外表络合PHREEQC参加了Dzombak和Morel〔1990〕创造的双层模型。双层模型是显式地计算扩散层组成〔BorkovecandWestall,1983〕和非静电层外表络合〔DavisandKent,1990〕的模型。其它模型，如三层、四层模型在PHREEQC中当前仍无法处理，根据Langmuir或Freundlich的等温吸附理论的吸附模拟可以处理成非静电模型的一种特例。DavisandKent〔1990〕评论外表络合模拟时指出使用物质的浓度作为被吸附组分的标准状态在理论上是有问题的。PHREEQC版本2使用的是外表组分的物质的分数作为其活度而不是物质的浓度。在标准状态下，这一变化对单齿状外表组分是没有影响的，但对多齿状外表那么影响很大。外表络合位置数目、外表积、被吸附组分、适宜的等确实定也是不确定的，实际模拟时需要用实际的野外介质进行实验研究，通过获取的实验数据来确定适宜的外表络合模型，并参加PHREEQC中进行模拟。4.固溶液PHREEQC使用Guggenheim方法来确定非理想的二元固溶液〔GlynnandReardon,1990〕，尚无法计算非理想的三元固溶液。可以通过假设为理想状态来模拟两种或更多种成分的固溶液，但是这种理想状态的假设过于简化，除相同元素的同位素假设处于理想状态之外，其他情况不行。5.溶质运移模拟模拟一维对流弥散或滞留区中的扩散时，PHREEQC使用的是显式优先差分算法。当网格划分较粗时，此算法出现数值弥散。大量的数值弥散取决于反响模型本身，数值弥散在有些情况下，可能较严重，如线性离子交换、外表络合滞留区的扩散，但当化学反响抵消弥散效应时，数值弥散较小。因此，推荐的做法是逐步模拟。首先用粗网格尽快获得结果并研究水化学反响过程，最后用细化网格评价数值弥散对反响组分和结果组分的影响。6.收敛问题PHREEQC可以发现输入的数据是否存在错误，但不能发现所建立的化学模型的无理性错误。例如，可通过元素的添加或减少来使溶液到达电荷平衡。如果元素无电离组分，或者即使元素的浓度减少到零时仍存在电荷平衡，那么数值方法不能收敛。其它的物理性错误有：〔1〕当想到达某固定浓度时而参加了碱〔事实上应参加酸〕；〔2〕当非碳酸碱度已经超过了数据输入中的总碱度时数值方法不能收敛。不过现在数值方法已经非常稳定了，所有的收敛问题已解决。偶尔需要使用关键字KNOBS的缩放属性，这种情况只有在总溶解度小于1×10-15mol/kg时才需使用。7.逆向模拟与以前的逆向模拟相比，PHREEQC最大的进步是在确定你想模型的处理时引入了非确定性设置。然而，由于算法处理小数据的方法问题，这种算法显示出的结果中有矛盾的现象。通过修改容许度可使算法能弥补这一问题。此外，PHREEQC中不能处理矿物沉淀时同位素的Rayleigh分馏效应，这也是一个缺乏之处。PHREEQCI—AGraphicalUserInterfacetothe

GeochemicalModelPHREEQCU.S.GeologicalSurveyApril2002IntroductionPHREEQCIVersion2是地球化学计算机程序PHREEQC(Version2)完全基于窗口的图形用户界面；PHREEQCI可用于交互式地执行PHREEQC〔第二版〕所有的模拟计算能力—物种形成、批量反响、一维(1D)反响-传输、以及逆向模拟。ApplicationsMinedrainage矿山排水Radioactive-wasteisolation放射性废弃物隔离Contaminantmigration污染物迁移Naturalandengineeredaquiferremediation天然或工程含水层修复Aquiferstorageandrecovery含水层蓄水和恢复Watertreatment水处理Naturalsystems天然系统Laboratoryexperiments实验室试验FeaturesoftheinterfaceCompatiblewithWindows95and98,WindowsNT4.0,WindowsME,Windows2000,andWindowsXPAutomatedinstallationusingInstallShieldInteractiveaccesstoallofthecapabilitiesofPHREEQC—createandmodifyinputdatafiles,runsimulations,anddisplayresultsModelingcapabilitiesAqueous,mineral,gas,surface,ion-exchange,andsolid-solutionequilibriaKineticreactions1Ddiffusionoradvectionanddispersionwithdual-porositymediumApowerfulinversemodelingcapabilityallowsidentificationofreactionsthataccountforthechemicalevolutioninobservedwatercompositionsExtensivegeochemicaldatabasesOperationPHREEQCIVersion2isacomputerprogramforsimulatingchemicalreactionsandtransportprocessesinnaturalorcontaminatedwater.PHREEQCIprovidesallofthecapabilitiesofthegeochemicalmodelPHREEQC(ParkhurstandAppelo,1999),including：speciation,batch-reaction,1Dreactive-transport,andinversemodeling.Operation〔Continue〕DataforPHREEQCareenteredthroughaseriesofkeyworddatablocks,eachofwhichprovidesaspecifictypeofinformation.Forexample,theSOLUTIONkeyworddatablockdefinesthechemicalcompositionofasolution.PHREEQCIprovidestabbeddialogboxesforeachPHREEQCkeyword.Figure1.MainscreenofPHREEQCI,withtheSOLUTIONkeyworddialogboxopen.1.SpeciationModelingApplications—Speciationmodelingisusefulinsituationswherethepossibilityofmineraldissolutionorprecipitationneedstobeknown,asinwatertreatment,aquiferstorageandrecovery,artificialrecharge,andwellinjection.Description—Speciationmodelingusesachemicalanalysisofawatertocalculatethedistributionofaqueousspeciesbyusinganion-associationaqueousmodel.Themostimportantresultsofspeciationcalculationsaresaturationindicesforminerals,whichindicatewhetheramineralshoulddissolveorprecipitate.Keyworddatablocksusedinspeciationmodelinginclude:SOLUTION—definethechemicalcompositionofasolution.SOLUTION_SPREAD—definethechemicalcompositionofmultiplesolutionsinaspreadsheetformat.2.Batch-ReactionModelingApplications—Batch-reactionmodelingcanbeappliedtoproblemsinlaboratory,natural,orcontaminatedsystems.ThereactioncapabilitiesofPHREEQChavebeenusedfrequentlyinthestudyofminedrainage,radioactive-wastedisposal,degradationoforganicmatter,andmicrobiallymediatedreactions.Batch-ReactionModeling(Continue)Description—Batchreactionscanbedividedintoequilibriumandnonequilibriumreactions.Equilibriumreactionsincludeequilibrationofasolutionwithanassemblageofmineralsandspecified-pressuregases,ion-exchangesites,surface-complexationsites,afinitegasphase,and(or)solidsolutions.Nonequilibriumreactionsincludekineticreactions,additionorremovalofelementsfromsolution,mixing,andchangingtemperature.Keyworddatablocksusedinbatch-reactionmodelinginclude:EQUILIBRIUM_PHASES—anassemblageofmineralsandgasesthatreacttoequilibrium(oruntilexhausted).EXCHANGE—oneormoresetsofreactionsitesforexchangeableions,animportantreactionformajorcations.SURFACE—oneormoresetsofsitesthatreactbysurfacecomplexation,amajorreactionforphosphorus,arsenic,andothertraceelements.GAS_PHASE—afinitereservoirofgasthatreactswithasolution.Keyworddatablocksusedinbatch-reactionmodeling(Continue)SOLID_SOLUTIONS—solidsthatprecipitateasmixturesofminerals,animportantreactionforradionuclidesandothertracemetals.KINETICS—anynonequilibriumreaction,forwhicharateexpressioncanbeformulated.RATES—BASIClanguagestatementsdefinerateexpressionsforkineticreactions.REACTION—additionandremovalofspecifiedelementsfromsolution.Keyworddatablocksusedinbatch-reactionmodeling(Continue)MIX—mixingtogetherspecifiedfractionsofsolutions.REACTION_TEMPERATURE—changingthetemperatureofthereactionsystem.USE—usethecompositionofpreviouslydefinedorsaved(SAVEkeyword)SOLUTION,EQUILIBRIUM_PHASES,EXCHANGE,SURFACE,GAS_PHASE,SOLID_SOLUTION,orapreviouslydefinedKINETICS,REACTION,MIX,orREACTION_TEMPERATUREdatablock.SAVE—savethecompositionfollowingthebatchreactionoftheSOLUTION,EQUILIBRIUM_PHASES,EXCHANGE,SURFACE,GAS_PHASE,orSOLID_SOLUTIONforuseinsubsequentreactions.3.TransportModelingApplications—Thereactivetransportcapabilitiescanbeusedtostudy：contaminantmigrationofnutrients,metals,radionuclides,andorganiccompounds;naturalandengineeredaquiferremediation;diffusioninsedimentporewater;thechemicalevolutionofnaturalsystems;andlaboratorycolumnexperiments.TransportModeling(Continue)Description—Reactive-transportmodelingsimulatesadvection,dispersion,andchemicalreactionsaswatermovesthrougha1Dcolumn.Thecolumnisdividedintoanumberofcells,andreactantcompositionsandnonequilibriumreactionscanbedefinedforeachcell.Allofthereactantsandreactionsdescribedforbatchreactionscanbeappliedtothecellsfortransportmodeling.TheTRANSPORTkeyworddatablockisusedtosimulateadvectionanddispersion,orpurediffusion,inthecolumn.TheTRANSPORTdatablockalsocanbeusedtosimulateadual-porositymedium,wherewaterflowsthroughthecolumnbutsimultaneouslyallowsfordiffusionintostagnantsidepores.TheADVECTIONkeyworddatablockisusedtosimulatepurelyadvectivetransport.Keyworddatablocksfortransportcalculations:TRANSPORT—1Dadvection,dispersion,andreactionmodeling.ADVECTION—1Dadvectionandreactionmodeling.4.InverseModelingApplications—Inversemodelingcanbeusedtodeducegeochemicalreactionsandmixinginlocalandregionalaquifersystems,andinaquiferstorageandrecoverystudies.Description—Inversemodelingcalculatesgeochemicalreactionsthataccountforthechangeinchemicalcompositionofwateralongaflowpath.Forinversemodeling,atleasttwochemicalanalysesofwa