第一讲PHREEQC软件简介-地下水污染与防治解析ppt课件.ppt

《第一讲PHREEQC软件简介-地下水污染与防治解析ppt课件.ppt》由会员分享，可在线阅读，更多相关《第一讲PHREEQC软件简介-地下水污染与防治解析ppt课件.ppt（33页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、PHREEQC软件简介 我吓了一跳，蝎子是多么丑恶和恐怖的东西，为什么把它放在这样一个美丽的世界里呢？但是我也感到愉快，证实我的猜测没有错：表里边有一个活的生物PHREEQC的发展过程oPHREEQC之前的版本是PHREEQE。PHREEQE是1980年美国地质调查局的Plummer 和Parkhurst等开发出来的，用于地球化学模拟达15年之久。oPHREEQC于1995年正式推出，它是在PHREEQE的源程序基础上用C语言重写而成的，消除了PHREEQE的缺陷与局限性。PHREEQC发展至今已近10年，其版本不断更新，加入了很多新内容，已发展到版本2中的2.11子版本。由于其现有的模拟能力

2、在同类软件中最强大，所以在国际上被广泛使用。我吓了一跳，蝎子是多么丑恶和恐怖的东西，为什么把它放在这样一个美丽的世界里呢？但是我也感到愉快，证实我的猜测没有错：表里边有一个活的生物PHREEQC的功能PHREEQC以离子联系的水化学模型为基础，可以推算：o（1）生成物和饱和系数；o（2）涉及到可逆反应以及不可逆反应的批反应和一维（1D）运移计算；其中，n涉及到的可逆反应包括水、矿物/无机溶液、气体、固体溶液、表面络合、离子交换平衡；n不可逆反应包括指定成分摩尔转换、动态控制反应、溶液混合和温度变化；o（3）逆向建模，可以发现一系列矿物和气体克分子在水中以不确定组成的转移。我吓了一跳，蝎子是多么

3、丑恶和恐怖的东西，为什么把它放在这样一个美丽的世界里呢？但是我也感到愉快，证实我的猜测没有错：表里边有一个活的生物PHREEQC 第二版的新特点和第一版相比，PHREEQC 第二版的新特点如下：o具有在一维运移计算中模拟弥散（或扩散）和滞流区的能力，使用用户确定的速率表达式模拟分子反应，模拟标准的多种成分或非标准的两种成分的固体溶液的沉淀和溶解，o模拟定体积气相和定压力气相，考虑表面系数或交换位置随着无机物的溶解和沉淀或者分子反应的变化而变化，自动采用多套收敛参数，o打印用户指定量到原始输出文件和（或）适合输入到扩展表格的文件上，以一种与扩展表程序更兼容的形式确定溶解成分等。我吓了一跳，蝎子是

4、多么丑恶和恐怖的东西，为什么把它放在这样一个美丽的世界里呢？但是我也感到愉快，证实我的猜测没有错：表里边有一个活的生物二、PHREEQC的局限性1. 水溶液模型oPHREEQC用相关的离子和Debye-Hckel方程计算非理想溶液，这种类型的水溶液模型在较低的离子强度下比较合适，但在高离子强度下（海水范围或更高）则会失真。o另外，数据库缺乏内部的连续性。因为随程序发布的两个数据库phreeqc.dat和wareq.dat中的lgKs 和反应焓来自于各种文献，未曾系统地确定溶液模型中单个lgK，也未曾确定当前数据库中水溶液模型是否与原始的实验数据相适应。使用时，随程序提供的数据文件可首先考虑，但

5、使用哪种水溶液组分和热力学数据则由用户自己确定。KlgKlg我吓了一跳，蝎子是多么丑恶和恐怖的东西，为什么把它放在这样一个美丽的世界里呢？但是我也感到愉快，证实我的猜测没有错：表里边有一个活的生物2. 离子交换模型o离子交换模型假设交换组分的热力学活度等于它的当量分数，或者用当量分数乘以Debye-Hckel活度系数，以此定义交换组分的活度（Appelo, 1994），更加复杂的交换模型还不能实现。o在许多野外研究中，离子交换模拟需要用研究区的介质进行离子交换试验，通过实验获得的数据来确定合适的模型。我吓了一跳，蝎子是多么丑恶和恐怖的东西，为什么把它放在这样一个美丽的世界里呢？但是我也感到愉快

6、，证实我的猜测没有错：表里边有一个活的生物3. 表面络合oPHREEQC加入了Dzombak 和Morel（1990）发明的双层模型。双层模型是显式地计算扩散层组成（Borkovec and Westall, 1983）和非静电层表面络合（Davis and Kent, 1990）的模型。其它模型，如三层、四层模型在PHREEQC中当前仍无法处理，根据Langmuir或Freundlich的等温吸附理论的吸附模拟可以处理成非静电模型的一种特例。oDavis and Kent（1990）评论表面络合模拟时指出使用物质的浓度作为被吸附组分的标准状态在理论上是有问题的。PHREEQC版本2使用的是表

7、面组分的物质的分数作为其活度而不是物质的浓度。在标准状态下，这一变化对单齿状表面组分是没有影响的，但对多齿状表面则影响很大。表面络合位置数目、表面积、被吸附组分、合适的等的确定也是不确定的，实际模拟时需要用实际的野外介质进行实验研究，通过获取的实验数据来确定合适的表面络合模型，并加入PHREEQC中进行模拟。我吓了一跳，蝎子是多么丑恶和恐怖的东西，为什么把它放在这样一个美丽的世界里呢？但是我也感到愉快，证实我的猜测没有错：表里边有一个活的生物4. 固溶液oPHREEQC使用Guggenheim方法来确定非理想的二元固溶液（Glynn and Reardon, 1990），尚无法计算非理想的三元

8、固溶液。可以通过假设为理想状态来模拟两种或更多种成分的固溶液，但是这种理想状态的假设过于简化，除相同元素的同位素假设处于理想状态之外，其他情况不行。我吓了一跳，蝎子是多么丑恶和恐怖的东西，为什么把它放在这样一个美丽的世界里呢？但是我也感到愉快，证实我的猜测没有错：表里边有一个活的生物5. 溶质运移模拟o模拟一维对流弥散或滞留区中的扩散时，PHREEQC使用的是显式优先差分算法。当网格划分较粗时，此算法出现数值弥散。大量的数值弥散取决于反应模型本身，数值弥散在有些情况下，可能较严重，如线性离子交换、表面络合滞留区的扩散，但当化学反应抵消弥散效应时，数值弥散较小。o因此，推荐的做法是逐步模拟。首先

9、用粗网格尽快获得结果并研究水化学反应过程，最后用细化网格评价数值弥散对反应组分和结果组分的影响。 我吓了一跳，蝎子是多么丑恶和恐怖的东西，为什么把它放在这样一个美丽的世界里呢？但是我也感到愉快，证实我的猜测没有错：表里边有一个活的生物6. 收敛问题oPHREEQC可以发现输入的数据是否存在错误，但不能发现所建立的化学模型的无理性错误。例如，可通过元素的添加或减少来使溶液达到电荷平衡。如果元素无电离组分，或者即使元素的浓度减少到零时仍存在电荷平衡，则数值方法不能收敛。其它的物理性错误有：n（1）当想达到某固定浓度时而加入了碱（事实上应加入酸）；n（2）当非碳酸碱度已经超过了数据输入中的总碱度时数

10、值方法不能收敛。o不过现在数值方法已经非常稳定了，所有已知的收敛问题已解决。偶尔需要使用关键字KNOBS的缩放属性，这种情况只有在总溶解度小于110-15 mol/kg时才需使用。我吓了一跳，蝎子是多么丑恶和恐怖的东西，为什么把它放在这样一个美丽的世界里呢？但是我也感到愉快，证实我的猜测没有错：表里边有一个活的生物7. 逆向模拟o与以前的逆向模拟相比，PHREEQC最大的进步是在确定你想模型的处理时引入了非确定性设置。然而，由于算法处理小数据的方法问题，这种算法显示出的结果中有矛盾的现象。通过修改容许度可使算法能弥补这一问题。o此外，PHREEQC中不能处理矿物沉淀时同位素的Rayleigh分

11、馏效应，这也是一个不足之处。PHREEQCIA Graphical User Interface to the Geochemical Model PHREEQC U.S. Geological Survey April 2002 U.S. Geological Survey April 2002 我吓了一跳，蝎子是多么丑恶和恐怖的东西，为什么把它放在这样一个美丽的世界里呢？但是我也感到愉快，证实我的猜测没有错：表里边有一个活的生物IntroductionoPHREEQCI Version 2 是地球化学计算机程序PHREEQC (Version 2)完全基于窗口的图形用户界面；oPHREEQ

12、CI 可用于交互式地执行 PHREEQC （第二版）所有的模拟计算能力物种形成、批量反应、一维(1D)反应-传输、以及逆向模拟 。我吓了一跳，蝎子是多么丑恶和恐怖的东西，为什么把它放在这样一个美丽的世界里呢？但是我也感到愉快，证实我的猜测没有错：表里边有一个活的生物ApplicationsoMine drainage 矿山排水oRadioactive-waste isolation 放射性废弃物隔离oContaminant migration 污染物迁移oNatural and engineered aquifer remediation 天然或工程含水层修复oAquifer storage

13、and recovery 含水层蓄水和恢复oWater treatment 水处理oNatural systems 天然系统oLaboratory experiments 实验室试验我吓了一跳，蝎子是多么丑恶和恐怖的东西，为什么把它放在这样一个美丽的世界里呢？但是我也感到愉快，证实我的猜测没有错：表里边有一个活的生物Features of the interfaceoCompatible with Windows 95 and 98,Windows NT 4.0, Windows ME, Windows2000, and Windows XP oAutomated installation u

14、sing InstallShield oInteractive access to all of the capabilities of PHREEQCcreate and modify input data files, run simulations, and display results 我吓了一跳，蝎子是多么丑恶和恐怖的东西，为什么把它放在这样一个美丽的世界里呢？但是我也感到愉快，证实我的猜测没有错：表里边有一个活的生物Modeling capabilitiesoAqueous, mineral, gas, surface, ion-exchange, and solid-solut

15、ion equilibria oKinetic reactions o1D diffusion or advection and dispersion with dual-porosity medium oA powerful inverse modeling capability allows identification of reactions that account for the chemical evolution in observed water compositions oExtensive geochemical databases 我吓了一跳，蝎子是多么丑恶和恐怖的东西

16、，为什么把它放在这样一个美丽的世界里呢？但是我也感到愉快，证实我的猜测没有错：表里边有一个活的生物Operation oPHREEQCI Version 2 is a computer program for simulating chemical reactions and transport processes in natural or contaminated water. oPHREEQCI provides all of the capabilities of the geochemical model PHREEQC (Parkhurst and Appelo, 1999), i

17、ncluding：n speciation, nbatch-reaction, n1D reactive-transport, nand inverse modeling. 我吓了一跳，蝎子是多么丑恶和恐怖的东西，为什么把它放在这样一个美丽的世界里呢？但是我也感到愉快，证实我的猜测没有错：表里边有一个活的生物Operation（Continue）oData for PHREEQC are entered through a series of keyword data blocks, each of which provides a specific type of information.

18、nFor example, the SOLUTION keyword data block defines the chemical composition of a solution. PHREEQCI provides tabbed dialog boxes for each PHREEQC keyword. 我吓了一跳，蝎子是多么丑恶和恐怖的东西，为什么把它放在这样一个美丽的世界里呢？但是我也感到愉快，证实我的猜测没有错：表里边有一个活的生物Figure 1. Main screen of PHREEQCI, with the SOLUTION keyword dialog box op

19、en. 我吓了一跳，蝎子是多么丑恶和恐怖的东西，为什么把它放在这样一个美丽的世界里呢？但是我也感到愉快，证实我的猜测没有错：表里边有一个活的生物1. Speciation ModelingoApplicationsSpeciation modeling is useful in situations where the possibility of mineral dissolution or precipitation needs to be known, as in water treatment, aquifer storage and recovery, artificial rech

20、arge, and well injection. oDescriptionSpeciation modeling uses a chemical analysis of a water to calculate the distribution of aqueous species by using an ion-association aqueous model.n The most important results of speciation calculations are saturation indices for minerals, which indicate whether

21、 a mineral should dissolve or precipitate. 我吓了一跳，蝎子是多么丑恶和恐怖的东西，为什么把它放在这样一个美丽的世界里呢？但是我也感到愉快，证实我的猜测没有错：表里边有一个活的生物Keyword data blocks used inspeciation modeling include: SOLUTIONdefine the chemical composition of a solution. SOLUTION_SPREADdefine the chemical composition of multiple solutions in a spre

22、adsheet format. 我吓了一跳，蝎子是多么丑恶和恐怖的东西，为什么把它放在这样一个美丽的世界里呢？但是我也感到愉快，证实我的猜测没有错：表里边有一个活的生物2. Batch-Reaction Modeling oApplicationsBatch-reaction modeling can be applied to problems in laboratory, natural, or contaminated systems. nThe reaction capabilities of PHREEQC have been used frequently in the study

23、 of mine drainage, radioactive-waste disposal, degradation of organic matter, and microbially mediated reactions. 我吓了一跳，蝎子是多么丑恶和恐怖的东西，为什么把它放在这样一个美丽的世界里呢？但是我也感到愉快，证实我的猜测没有错：表里边有一个活的生物Batch-Reaction Modeling (Continue)oDescriptionBatch reactions can be divided into equilibrium and nonequilibrium react

24、ions. nEquilibrium reactions include equilibration of a solution with an assemblage of minerals and specified-pressure gases, ion-exchange sites, surface-complexation sites, a finite gas phase, and (or) solid solutions. nNonequilibrium reactions include kinetic reactions, addition or removal of elem

25、ents from solution, mixing, and changing temperature. 我吓了一跳，蝎子是多么丑恶和恐怖的东西，为什么把它放在这样一个美丽的世界里呢？但是我也感到愉快，证实我的猜测没有错：表里边有一个活的生物Keyword data blocks used in batch-reaction modeling include: oEQUILIBRIUM_PHASESan assemblage of minerals and gases that react to equilibrium (or until exhausted). oEXCHANGEone o

26、r more sets of reaction sites for exchangeable ions, an important reaction for major cations. oSURFACEone or more sets of sites that react by surface complexation, a major reaction for phosphorus, arsenic, and other trace elements. oGAS_PHASEa finite reservoir of gas that reacts with a solution. 我吓了

27、一跳，蝎子是多么丑恶和恐怖的东西，为什么把它放在这样一个美丽的世界里呢？但是我也感到愉快，证实我的猜测没有错：表里边有一个活的生物Keyword data blocks used in batch-reaction modeling (Continue)oSOLID_SOLUTIONSsolids that precipitate as mixtures of minerals, an important reaction for radio nuclides and other trace metals. oKINETICSany nonequilibrium reaction, for w

28、hich a rate expression can be formulated. oRATESBASIC language statements define rate expressions for kinetic reactions. oREACTIONaddition and removal of specified elements from solution. 我吓了一跳，蝎子是多么丑恶和恐怖的东西，为什么把它放在这样一个美丽的世界里呢？但是我也感到愉快，证实我的猜测没有错：表里边有一个活的生物Keyword data blocks used in batch-reaction m

29、odeling (Continue)oMIXmixing together specified fractions of solutions. oREACTION_TEMPERATUREchanging the temperature of the reaction system. oUSEuse the composition of previously defined or saved (SAVE keyword) SOLUTION, EQUILIBRIUM_PHASES, EXCHANGE, SURFACE, GAS_PHASE, SOLID_SOLUTION, or a previou

30、sly defined KINETICS, REACTION, MIX, or REACTION_TEMPERATURE data block. oSAVEsave the composition following the batch reaction of the SOLUTION, EQUILIBRIUM_PHASES, EXCHANGE, SURFACE, GAS_PHASE, or SOLID_SOLUTION for use in subsequent reactions. 我吓了一跳，蝎子是多么丑恶和恐怖的东西，为什么把它放在这样一个美丽的世界里呢？但是我也感到愉快，证实我的猜测

31、没有错：表里边有一个活的生物3. Transport Modeling oApplicationsThe reactive transport capabilities can be used to study ：ncontaminant migration of nutrients, metals, radio nuclides, and organic compounds; nnatural and engineered aquifer remediation; ndiffusion in sediment pore water; nthe chemical evolution of na

32、tural systems; nand laboratory column experiments. 我吓了一跳，蝎子是多么丑恶和恐怖的东西，为什么把它放在这样一个美丽的世界里呢？但是我也感到愉快，证实我的猜测没有错：表里边有一个活的生物Transport Modeling (Continue)oDescriptionReactive-transport modeling simulates advection, dispersion, and chemical reactions as water moves through a 1D column. The column is divide

33、d into a number of cells, and reactant compositions and nonequilibrium reactions can be defined for each cell. All of the reactants and reactions described for batch reactions can be applied to the cells for transport modeling. nThe TRANSPORT keyword data block is used to simulate advection and disp

34、ersion, or pure diffusion, in the column. nThe TRANSPORT data block also can be used to simulate a dual-porosity medium, where water flows through the column but simultaneously allows for diffusion into stagnant side pores. nThe ADVECTION keyword data block is used to simulate purely advective trans

35、port. 我吓了一跳，蝎子是多么丑恶和恐怖的东西，为什么把它放在这样一个美丽的世界里呢？但是我也感到愉快，证实我的猜测没有错：表里边有一个活的生物Keyword data blocks for transport calculations: oTRANSPORT1D advection, dispersion, and reaction modeling. oADVECTION1D advection and reaction modeling. 我吓了一跳，蝎子是多么丑恶和恐怖的东西，为什么把它放在这样一个美丽的世界里呢？但是我也感到愉快，证实我的猜测没有错：表里边有一个活的生物4. In

36、verse Modeling oApplicationsInverse modeling can be used to deduce geochemical reactions and mixing in local and regional aquifer systems, and in aquifer storage and recovery studies. oDescriptionInverse modeling calculates geochemical reactions that account for the change in chemical composition of

37、 water along a flow path. nFor inverse modeling, at least two chemical analyses of water are needed at different points along the flow path, as well as a set of minerals and gases that are potentially reactive. nMole transfers of phases are calculated that account for the change in water composition

38、 along the flow path. nThe numerical method accounts for uncertainties in analytical data.我吓了一跳，蝎子是多么丑恶和恐怖的东西，为什么把它放在这样一个美丽的世界里呢？但是我也感到愉快，证实我的猜测没有错：表里边有一个活的生物Keyword data block for inverse modeling: INVERSE_MODELINGdeduce mixing, mineral, and gas reactions that account for the chemical evolution of

39、waters. 我吓了一跳，蝎子是多么丑恶和恐怖的东西，为什么把它放在这样一个美丽的世界里呢？但是我也感到愉快，证实我的猜测没有错：表里边有一个活的生物How to Obtain PHREEQCI and PHREEQC?oVersions of PHREEQCI and PHREEQC may be obtained at no cost from the USGS at World Wide Web address: http:/water.usgs.gov/software oMany other hydrologic computer models are also available

40、 from this site. 我吓了一跳，蝎子是多么丑恶和恐怖的东西，为什么把它放在这样一个美丽的世界里呢？但是我也感到愉快，证实我的猜测没有错：表里边有一个活的生物ReferencesoParkhurst, D.L. and Appelo, C.A.J., 1999, Users guide to PHREEQC (Version2) A computer program for speciation, batch-reaction, one-dimensional transport, and inverse geochemical calculations: U.S. Geological Survey Water-Resources Investigations Report 99-4259, 310 p.