phreeqc实例

合集下载

phreeqc中文使用手册

phreeqc中文使用手册

phreeqc中文使用手册【最新版】目录1.Phreeqc 简介2.安装与配置 Phreeqc3.Phreeqc 的基本语法4.Phreeqc 的应用示例5.总结正文【1.Phreeqc 简介】Phreeqc 是一个用于处理化学问题的开源软件包,特别是在处理矿物溶解度、岩石水化学和环境地球化学方面的问题非常有效。

它基于 Python 编程语言,通过大量的函数和方法,为用户提供了一个强大的工具来解决各种化学问题。

【2.安装与配置 Phreeqc】在使用 Phreeqc 之前,你需要先安装 Python,并使用 pip 安装Phreeqc。

安装过程非常简单,只需要在终端中输入以下命令:```pip install phreeqc```安装完成后,你可以通过编写 Python 脚本来使用 Phreeqc。

需要注意的是,Phreeqc 对于 Python 的版本有一定的要求,如果你使用的Python 版本过低,可能需要升级到新版本才能正常使用 Phreeqc。

【3.Phreeqc 的基本语法】Phreeqc 的基本语法非常简单,主要包括以下几个部分:- 导入模块:在使用 Phreeqc 之前,需要先导入 phreeqc 模块。

- 创建对象:Phreeqc 中的许多功能都需要创建对象,例如溶液对象、离子对象等。

- 调用方法:通过对象调用相应的方法,可以实现各种化学计算。

下面是一个简单的 Phreeqc 使用示例:```pythonimport phreeqc# 创建一个溶液对象solution = phreeqc.Solution("example.iqc")# 添加物质到溶液中solution.add("CaCO3", 1.0)solution.add("H2O", 1000.0)# 计算溶液中的离子浓度solution.calculate_ions()# 输出结果solution.print_ions()```【4.Phreeqc 的应用示例】下面是一个使用 Phreeqc 计算矿物溶解度的示例:```pythonimport phreeqc# 创建一个溶液对象solution = phreeqc.Solution("example.iqc")# 添加物质到溶液中solution.add("CaCO3", 1.0)solution.add("H2O", 1000.0)# 设置溶液的温度和压力solution.set_temperature(25)solution.set_pressure(101325)# 计算矿物的溶解度solution.calculate_dissolution("CaCO3")# 输出结果solution.print_dissolution()```【5.总结】Phreeqc 是一个功能强大的化学计算软件包,可以帮助用户解决各种化学问题。

phreeqc中文使用手册

phreeqc中文使用手册

phreeqc中文使用手册摘要:1.Phreeqc 简介2.安装与配置3.使用方法4.常见问题5.总结正文:1.Phreeqc 简介Phreeqc 是一个用于解决化学问题的计算机程序,特别是在矿物溶解度、离子交换、吸附和沉淀等方面具有很高的应用价值。

该程序采用离子强度罗兰兹法,可以模拟各种复杂的水文地球化学过程。

Phreeqc 具有丰富的功能和灵活性,适用于科研、教学和工程应用等多种场景。

2.安装与配置在使用Phreeqc 之前,首先需要进行程序的安装。

根据用户的操作系统,可以选择相应的安装包进行安装。

安装完成后,需要对Phreeqc 进行配置,包括设置计算参数、定义输入输出文件格式等。

具体配置方法可以参考Phreeqc 的用户手册或在线教程。

3.使用方法Phreeqc 的使用方法主要包括以下几个步骤:(1)准备输入文件:根据需要解决的问题,编写相应的输入文件。

输入文件主要包括化学物质的化学式、初始浓度、溶液的物理化学性质等。

(2)运行Phreeqc:在命令行界面,输入Phreeqc 命令并指定输入文件,启动程序进行计算。

(3)查看结果:Phreeqc 计算完成后,会生成相应的输出文件。

用户可以根据输出文件查看计算结果,包括溶液中各离子的浓度、pH 值、溶解度积等。

4.常见问题在使用Phreeqc 过程中,可能会遇到一些常见问题,如输入文件格式错误、计算结果不符合预期等。

针对这些问题,可以通过以下方法进行解决:(1)检查输入文件:确保输入文件格式正确,化学物质的化学式和初始浓度正确无误。

(2)调整计算参数:根据问题的实际情况,调整Phreeqc 的计算参数,如计算方法、离子强度等。

(3)参考教程和手册:查阅Phreeqc 的用户手册或在线教程,了解程序的使用方法和注意事项。

5.总结Phreeqc 是一个功能强大的化学计算程序,广泛应用于水文地球化学、环境科学等领域。

phreeqc中文使用手册

phreeqc中文使用手册

phreeqc中文使用手册摘要:1.phreeqc 简介2.phreeqc 的安装与配置3.phreeqc 的使用方法4.phreeqc 的常见问题与解决方案5.总结正文:【1.phreeqc 简介】phreeqc 是一款用于计算机视觉和机器学习的开源库,它提供了许多用于处理图像和视频的先进功能。

phreeqc 的命名来源于“快速、高效、易用”三个词的缩写,这也是它最显著的特点。

【2.phreeqc 的安装与配置】phreeqc 可以通过pip 进行安装,安装命令为:`pip install phreeqc`。

安装完成后,你需要配置phreeqc 的环境变量,以便在程序中调用它。

具体配置方法可以参考phreeqc 的官方文档。

【3.phreeqc 的使用方法】phreeqc 的使用方法非常简单。

首先,你需要导入phreeqc 库,然后调用相应的函数来实现你想要的功能。

例如,如果你想要对一张图片进行处理,你可以使用`phreeqc.process_image()`函数;如果你想要对一段视频进行处理,你可以使用`phreeqc.process_video()`函数。

【4.phreeqc 的常见问题与解决方案】在使用phreeqc 的过程中,你可能会遇到一些问题。

例如,如果你在使用过程中遇到了“找不到指定模块”的错误,那么可能是因为你没有正确安装phreeqc 或者没有正确配置环境变量。

解决这个问题的方法是检查你的安装和配置是否正确。

【5.总结】总的来说,phreeqc 是一款功能强大、易于使用的计算机视觉和机器学习库,它为开发者提供了许多方便实用的工具。

phreeqc中文使用手册

phreeqc中文使用手册

phreeqc中文使用手册【实用版】目录1.Phreeqc 简介2.安装与配置 Phreeqc3.Phreeqc 的使用方法4.Phreeqc 的功能与应用领域5.常见问题与解答正文【1.Phreeqc 简介】Phreeqc 是一款开源的中文自然语言处理工具,全称为“现代汉语语义分析工具集”。

它由哈尔滨工业大学的自然语言处理实验室研发,旨在为中文自然语言处理提供一套功能完善、性能高效的工具。

Phreeqc 集成了多种中文自然语言处理技术,如分词、词性标注、命名实体识别、依存句法分析等,广泛应用于文本分类、情感分析、信息抽取、机器翻译等领域。

【2.安装与配置 Phreeqc】Phreeqc 支持 Linux、macOS 和 Windows 操作系统。

安装过程相对简单,用户只需根据官方提供的安装指南进行操作即可。

对于 Linux 系统,可以使用包管理器如 apt、yum 等进行安装;对于 macOS 和 Windows 系统,则需要从源代码编译安装。

配置 Phreeqc 主要包括设置环境变量、数据路径、模型路径等。

环境变量用于指定 Phreeqc 的可执行文件位置;数据路径用于存储Phreeqc 的词库、规则文件等;模型路径则用于存储 Phreeqc 的模型文件。

用户可以根据自己的需求进行相应的配置。

【3.Phreeqc 的使用方法】Phreeqc 提供了丰富的 API 接口,用户可以通过命令行、Python、Java 等方式调用这些接口。

下面以命令行为例,介绍 Phreeqc 的基本使用方法:(1) 分词:使用`phreeqc_cut`命令进行分词,如`phreeqc_cut -i input.txt -o output.txt`。

(2) 词性标注:使用`phreeqc_pos`命令进行词性标注,如`phreeqc_pos -i input.txt -o output.txt`。

(3) 命名实体识别:使用`phreeqc_ner`命令进行命名实体识别,如`phreeqc_ner -i input.txt -o output.txt`。

phreeqc中文使用手册

phreeqc中文使用手册

phreeqc中文使用手册摘要:一、引言- 介绍phreeqc中文使用手册的目的和适用对象二、phreeqc软件概述- 解释phreeqc的含义和作用- 介绍phreeqc软件的发展历程和主要功能三、phreeqc中文使用手册的结构- 概述各个章节的内容和主题四、软件安装与配置- 说明软件的安装流程和注意事项- 介绍如何配置phreeqc以满足用户需求五、模型与方法- 详述phreeqc软件中包含的各种模型和计算方法- 解释如何选择合适的模型和方法进行计算六、输入与输出- 说明如何准备和输入数据- 介绍phreeqc软件的输出结果及其含义七、常见问题与解决方案- 分析在使用过程中可能遇到的问题- 提供相应的解决方法和技巧八、软件更新与维护- 介绍如何获取软件的最新版本- 说明如何进行软件的更新和维护九、结论- 总结phreeqc中文使用手册的主要内容- 强调软件在实际应用中的优势和价值正文:一、引言phreeqc中文使用手册是为了帮助用户更好地理解和使用phreeqc软件而编写的。

本手册适用于phreeqc软件的初学者和有经验的用户,旨在提供全面、详细的软件操作指南和实用技巧。

通过阅读本手册,用户可以掌握phreeqc软件的使用方法,充分发挥其在地球化学计算领域的优势。

二、phreeqc软件概述phreeqc是一款功能强大的地球化学计算软件,全称为“PHREEQC - A Geochemical Modeling Program for Speciation, Batch, and One-Dimensional Transport”。

该软件由美国地质调查局(USGS)开发,基于PHREEQE软件,适用于多种场景下的地球化学计算。

phreeqc软件可以用于解决以下问题:1.分析化学反应在特定条件下的平衡性质;2.计算溶液中各种离子的活度系数;3.模拟物质在环境中的迁移和转化过程;4.为地质调查、环境保护和资源开发等项目提供地球化学数据支持。

Phreeqc 使用教程

Phreeqc 使用教程
3. Now you should have the most
up to date database in your GEMS programme.
2. Finding information in GEMS
1. Open the GEMS project you
2.
wish to find the thermodynamic 3. info you want to use.
1. 定义物相. 2. 选择物相.
3.反应物的初始摩尔数。如果缺省, 默认为1mol 4. 反应物当前摩尔数 1 1
当反应发生时,这个数值会增加或者减少。 如果缺省,默认为1mol
5. 速率表达式中的参数。 6. 积分步骤公差(mols) 对每一个时间间隔的每一步积分而言, 速率表达式的第五次序和第四次序之间 的不同一定会小于这个公差,或时间间 隔自动减少。反应重要性小的一些浓度 差异需要更小一些的公差。动力学积分 数字上的准确性可通过减少公差实现, 以决定结果改变的重要性。缺省为10-8. 2 3 4 6 5
1
2
3
二、Solution溶液
1.选择溶液的温度,加入水质量,pH,密度。 2.选择进入元 素成分,选择 溶液中元素。 3. 注意单位。 4. 选中元素,注 意填入的是浓度, 不是物质的量。 1 3 pH 温度 溶液密度 水质量 1
2
三、 Equilibrium phase平衡相。
1. 选择所需要的反应物与产物, 跟据自己研究的条件,决定是 否选择(因为reaction 模块中 也可以添加反应物)。
2. Go into ‘Database mode’.
3. Click on ‘DComp’. Search for

PHREEQCI模拟软件简介

PHREEQCI模拟软件简介

PHREEQC是由美国地调所开发的水文地球化学模拟软件,是用C语言编写的进行低温水文地球化学计算的计算机程序,可进行正向模拟和反向模拟,几乎能解决水、气、岩土相互作用系统中所有平衡热力学和化学动力学问题,包括水溶物配合、吸附一解吸、离子交换、表面配合、溶解沉淀、氧化一还原。

PHREEqc是一个加入了非确定项控制的水一岩作用模拟软件。

正向模拟能根据给定的反应机理来预测水的组分和质量的迁移,可以进行配分和饱和指数以及一次投药反应和一维运移计算。

反向模拟根据观测的化学和同位素资料来确定水岩反应机理,即计算造成水流途径上初始和最终水组分差异所必须溶解或析出的矿物和气体物质的量。

PHREEQC有一个强大的热力学数据库供输入和运行使用。

该数据库主要由 phreeqc.dat、wateq4f.dat、minteq.dat、LLNL.dat 4个数据库组成。

每一数据库均有水溶液主要组分,水溶液一般组分,相 (气体和矿物 )、表面主要组分和表面一般组分数据块。

PHREEQC在不同化学条件下的地下水化学环境中的应用

PHREEQC在不同化学条件下的地下水化学环境中的应用

PHREEQC在不同化学条件下的地下水化学环境中的应用作者摘要:在高矿化度地下水分布地区实行地浸采矿是一项世界性难题, 因为当地下水矿化度超过5g/ L时就不适合进行地浸。

我国某砂岩型矿含矿含水层地下水的矿化度达到8~12 g/ L ,无论是采用酸法和碱法地浸,都发生堵塞现象。

为了解决这个难题,通过淡化,降低矿区地下水的矿化度,即降低水中Ca2 + 、Mg2 + 及SO42 - 的含量,从而避免了石膏的析出。

在充入CO2 气体的条件下,碳酸盐的饱和指数降为负值。

从而达到地浸顺利进行的目的。

应用PHREEQC模拟软件,对我国某高矿化度地下水的矿地浸溶浸剂进行模拟,通过计算饱和指数,确定溶质水解沉淀的水文地球化学条件的临界值。

研究在该地地浸工艺过程中溶质的存在形式,以及用常规地浸产生沉淀堵塞的原因,为解决地下水高矿化度的难题提供依据和办法。

关键字:高矿化度;地浸;淡化;PHREEQCAbstract:In high salinity groundwater distribution area of in-situ leaching mining is a worldwide problem, because the local water mineralization degree more than 5g/ L is not suitable for in-situ leaching. A sandstone type ore in China ore-bearing aquifer groundwater mineralization degree reached 8 ~ 12 g/ L, either by acid and alkali leaching, are blocked. In order to solve this problem, through desalination, reduce the degree of mineralization of groundwater in the lower water content, Ca2 +, Mg2 + and SO42 -, so as to avoid the precipitation of gypsum. In filling CO2 gas conditions, the saturation index of carbonate fall below zero. So as to achieve the purpose of in-situ leaching smoothly. Application of PHREEQC simulation software, to simulate the high salinity groundwater ore leaching infusion, by calculating the saturation index, the critical value to determine the hydrogeochemical conditions of solute precipitation. In the form of in-situ leaching of solutes in the process, as well as with conventional immersion precipitation blockage reason, provides the basis and the way to solve the problem of groundwater with high mineralization.Keywords:High salinity; leaching; desalination; PHREEQC1.PHREEQC模拟软件简介PHREEQC 是由美国地调所(USGS)在PHREEQE 的基础上开发的用于计算多种低温水文地球化学反应的计算机软件[1] ,是用 C 语言编写的进行低温水文地球化学计算的计算机程序,可进行正向模拟和反向模拟,几乎能解决水、气、岩土相互作用系统中所有平衡热力学和化学动力学问题,包括水溶物配合、吸附-解吸、离子交换、表面配合、溶解-沉淀、氧化-还原。

phreeqc中文使用手册

phreeqc中文使用手册

phreeqc中文使用手册【原创版】目录1.Phreeqc 简介2.安装与配置 Phreeqc3.Phreeqc 的基本语法4.Phreeqc 的应用示例5.总结正文【1.Phreeqc 简介】Phreeqc 是一个用于解决化学问题的计算机程序,特别是在矿物加工和矿物水文学方面。

它是一个基于 Python 的免费软件,可以运行在Windows、Linux 和 Mac OS X 等操作系统上。

Phreeqc 提供了大量的功能,包括溶液模拟、沉淀模拟、酸碱滴定、质量平衡计算等。

【2.安装与配置 Phreeqc】Phreeqc 的安装相对简单。

首先,需要下载 Phreeqc 的源代码,并按照官方提供的安装指南进行安装。

安装完成后,需要配置 Phreeqc 的环境变量,以便在运行时可以找到 Phreeqc 的可执行文件。

此外,Phreeqc 还需要一些依赖库,例如 Python 的 SciPy 和 NumPy 库,需要提前安装。

【3.Phreeqc 的基本语法】Phreeqc 的基本语法包括两个部分:输入文件和命令行参数。

输入文件包含了要解决的问题的详细描述,包括物质的化学式、初始浓度等信息。

命令行参数则指定了要运行的 Phreeqc 程序的选项和参数。

【4.Phreeqc 的应用示例】以下是一个 Phreeqc 的应用示例:假设有一个含有铜离子和氢氧根离子的溶液,我们想要知道当 pH 值为多少时,铜离子会开始沉淀。

我们可以使用 Phreeqc 的沉淀模拟功能来解决这个问题。

首先,我们需要编写一个输入文件,描述溶液的初始条件和要解决的问题。

然后,我们在命令行中指定输入文件和沉淀模拟选项,运行 Phreeqc 程序。

最后,我们得到结果:当 pH 值小于 4.76 时,铜离子开始沉淀。

【5.总结】Phreeqc 是一个功能强大的化学计算软件,特别适用于矿物加工和矿物水文学领域。

它的安装和配置相对简单,使用起来也较为方便。

第一讲PHREEQC软件简介_地下水污染与防治解析

第一讲PHREEQC软件简介_地下水污染与防治解析

Introduction
PHREEQCI Version 2 是地球化学计算机程序 PHREEQC (Version 2)完全基于窗口的图形用 户界面;
PHREEQCI 可用于交互式地执行 PHREEQC (第二版)所有的模拟计算能力—物种形成、 批量反应、一维(1D)反应-传输、以及逆向模拟 。
(3)逆向建模,可以发现一系列矿物和气体克分子在水中以 不确定组成的转移。
PHREEQC 第二版的新特点
和第一版相比,PHREEQC 第二版的新特点如下: 具有在一维运移计算中模拟弥散(或扩散)和滞
流区的能力,使用用户确定的速率表达式模拟分 子反应,模拟标准的多种成分或非标准的两种成 分的固体溶液的沉淀和溶解, 模拟定体积气相和定压力气相,考虑表面系数或 交换位置随着无机物的溶解和沉淀或者分子反应 的变化而变化,自动采用多套收敛参数, 打印用户指定量到原始输出文件和(或)适合输 入到扩展表格的文件上,以一种与扩展表程序更 兼容的形式确定溶解成分等。
Kinetic reactions 1D diffusion or advection and dispersion with
dual-porosity medium A powerful inverse modeling capability allows
identification of reactions that account for the chemical evolution in observed water compositions Extensive geochemical databases
Operation
PHREEQCI Version 2 is a computer program for simulating chemical reactions and transport processes in natural or contaminated water.

phreeqc指数对方解石的饱和状况

phreeqc指数对方解石的饱和状况

phreeqc指数对方解石的饱和状况标题:探寻phreeqc指数对方解石饱和状况的关系介绍:在地球科学领域,方解石是一种常见的矿物,具有广泛的应用领域。

了解方解石的饱和状况对于地质学家和环境科学家来说至关重要。

本文将通过深入研究phreeqc指数及其与方解石饱和状况之间的关系,揭示其背后的科学原理和应用价值。

1. phreeqc指数的概念及其计算方法1.1 什么是phreeqc指数1.2 phreeqc指数的计算方法2. 方解石的饱和状况与phreeqc指数的关系2.1 方解石的定义和饱和状况2.2 方解石饱和度的表征方式2.3 phreeqc指数与方解石饱和状况的影响因素2.3.1 温度对方解石饱和度的影响2.3.2 溶液中的离子浓度对方解石饱和度的影响2.3.3 pH值对方解石饱和度的影响3. phreeqc指数的应用与意义3.1 环境科学领域中的应用3.2 地球科学领域中的应用3.3 工程领域中的应用4. 通过phreeqc指数了解方解石饱和状况的优势4.1 高精度的方解石饱和度计算4.2 丰富的实验数据支持4.3 可迁移性与普适性总结:本文深入研究了phreeqc指数与方解石饱和状况之间的关系。

通过phreeqc指数的计算和分析,我们可以更好地了解方解石的饱和度以及影响因素。

本文还分析了phreeqc指数在环境科学、地球科学和工程领域中的应用和意义。

通过对phreeqc指数的研究,我们能够更全面、深刻地理解方解石饱和状况,并为相关领域的研究和实践提供有力的支持。

观点与理解:phreeqc指数作为一种计算和表征方解石饱和状况的指标,在地球科学和环境科学领域扮演着重要角色。

它提供了高精度的方解石饱和度计算和预测,为我们深入探索方解石行为提供了有力工具。

phreeqc 指数的应用广泛,可在各种实际问题和工程项目中发挥作用,为解决相关问题提供科学依据和技术支持。

文章完。

phreeqc指数(研究方解石饱和状况的强有力指标)在研究方解石饱和状况中,phreeqc指数具有许多优势使得其成为一种重要的工具。

phreeqc实例

phreeqc实例

phreeqc实例例2——矿物相的溶解平衡这个例子测定了最稳定相——石膏或是无水石膏在一定温度范围内的溶解性。

输入数据组见表13。

仅用pH和温度来定义纯水溶液。

缺省单位是millimolal,但没有指定浓度。

缺省状态下,pe为4.0,缺省的氧化还原计算使用pe,且水的密度为1.0(这是没有必要的,因为没有浓度为“每升”)。

在批反应期间,所有允许反应达到指定的饱和指数的反应都列在EQUILIBRIUM_PHASES中,无论开始时它们是否存在。

输入数据包含相的名字(之前通过PHASES输入在数据库或输入文件中定义)、指定的饱和指数和以摩尔数表示的当前存在的相的数量。

如果一种相在当前不存在,则在纯相集合中给定其浓度为0.0 mol。

在这个例子中,石膏和无水石膏允许反应来达到平衡(饱和指数等于0.0),初始相的集合中每种矿物都为1mol。

每一种矿物或是反应达到平衡或是被耗尽。

在大多数的情况下,1 mol的单纯相足够达到反应平衡。

表13--例2输入数据的设置TITLE Example 2.--Temperature dependence of solubilityof gypsum and anhydriteSOLUTION 1 Pure waterpH 7.0temp 25.0EQUILIBRIUM_PHASES 1Gypsum 0.0 1.0Anhydrite 0.0 1.0REACTION_TEMPERATURE 125.0 75.0 in 51 stepsSELECTED_OUTPUTfile ex2.selsi anhydrite gypsumEND在REACTION_TEMPERATURE数据块中,指定反应温度计算步长为1?,从25?开始,75?结束,共计算温度变化51次。

温度的每一度变化都详细指定输入数据,在可能的情况下,在EQUILIBRIUM_PHASES(石膏和无水石膏)中定义的物质将会反应达到平衡,或是直到两种物质都完全溶解。

phreeqc实例练习

phreeqc实例练习

PHREEQC实例分析例1——物种形成分析这个例子计算了海水中矿物质的分布以及一组有关矿物在海水中的饱和程度。

为了证明如何在这个模型中应用新的元素,将元素铀添加入由phreeqc.dat定义的液相模型中[wateq.dat是包含于程序分类中的一个数据库文件,它来自于WATEQ4F(Ball and Nordstrom, 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硅石,SiO2Si 4.28氯化物Cl 19353.0碱度,HCO3-Alkalinity 141.682硫酸盐,SO42-S(6) 2712.0硝酸盐,NO3-N(5) .29铵,NH4+N(-3) .03铀U .0033pH,标准单位pH 8.22pe,无单位pe 8.451温度,℃temperature 25.0密度,千克/升density 1.023用于分配氧化还原元素和计算饱和指数的pe由redox标识符所指定。

在这个例子中,用氧化还原电对O(-2)/O(0) 计算的pe值相对应于溶解氧/水,并且这个pe适用于需要pe值的所有的计算。

如果redox没有指定,那么缺省的值将会是所输入的pe。

缺省的氧化还原标识符可被任何氧化还原元素代替,如输入元素锰时,则输入的pe被用来表示各种化合价状态的锰;输入铀时,这里是氮/铵电对将会用来计算所形成各种价态铀的pe值。

反向模拟在PHREEQC程序中的实现

反向模拟在PHREEQC程序中的实现

反向模拟在PHREEQC程序中的实现反匐模量薹薹RQC程序中的实现马骏中船勘察设计研究院2003331PHREEQC模型及原理简介PHREEQC是由美国地调所开发的水文地球化学模拟软件,它是C语言编写的进行低温水文地球化学计算的计算机程序,其主要是根据用户的输入命令,PHREEQC将选择其中的某些方程来描述相应的化学反应过程,这些方程组成的方程组,采用改进的牛顿拉斐逊(Newton-Raphson)方法进行迭代求解.PHREEQC共提供了四个数据库供用户进行选择应用.输入文件是需要用户编写的文本文件,文件给出命令(反应模式)供模型读入并进行模拟,也可以在此文件中对数据库进行修改和特别选择计算输出结果;标准输出文件是PHREEQC在模拟运算过程中的输出结果;选择性输出文件是根据用户需要选择性输出的计算结果.国内目前对PHREEQC的应用多限于进行化学组分的分析.2反向地球化学模拟反向地球化学模拟是依据观测到的水化学资料来确定系统中所进行的水一岩反应,也就是对观测到的水化学资料做出解释,目的是查找出地下水和不同矿物,气体之间发生复杂的反应状况,在合理条件下给予量化.在反向模拟模型中,从来不存在一个"唯一的"答案,而是一系列的可能路径,并且在此列出的其中之一.反向模拟计算中, 最困难的是反应阶段的选择,如化学恒等式和合成物的选择,在地质运用中,必须根据野外地质条件和室内地质实验确定的矿物类型来选择,包括显微分析,化学分析等.对于反应矿物的含量在本模型中不是要求精确度很高的,但是,由于水一岩作用的复杂性和地质矿物成分的复杂性,潜在的反应类型是必须要认清的,在模型输入时是必须的.为了方便,logK系数设为0, logK不影响本反向模型的过程和结果,因为在本模型中只对化学剂量进行计算,求出监测孔地下水中相应矿物的饱和指数等. 3反向模拟在PHREEQC程序中的实现过程水文地质研究的一个重要方面是确定地下水的来源,正确估计可能的地质的或人为的影响,以及地下水活动对工程的影响.于是人们想通过对地下水化学组成的认识来恢复其流程.例如,从所研究的井水的化学组成以及雨水的化学分析结果,可以推测,在雨水的下渗过程中,必须经过何种地质构造,通过与矿物相的反应(溶解,沉淀),使其化学组成发生改变,最终形成井中的水.同样的道理,我们可以通过对检测孔水质和原始库水的分析比较,大致确定坝基水的来源和演化过程.在输入数据前,需要对矿物组分的含量作变换,由mg/L转换为mol/L.在模拟进行前,一般检查每个分析的分析误差,由电荷守恒概念可知,需要对地下水分析结果中进行电荷平衡.如Ca,Mg等包含地下水指标信息多,在输入模型中,以Ca作调整平衡模型中电荷.取模型不确定系数为0.05,即5%.在程序中,不确定系数大于0.1,则最终计算结果没有大的意义.4实例研究以陈村大坝为例,分别以库水和监测孔水样作为初始条件水和反应后的水,水化学组分的差异归因于水一岩的相互作用.转幕|页表1大坝灌浆廊道地下水水质组分及含量统计表表2大坝排水廊道地下水水质组分及含量统计表水质组分岔量(molL)编号主要纽分灌10-4灌l7-4灌24-2溢25-1la2296e-05l994e-041753e一041.448e—O42K'1999e51440e-049.993e—O51.100e-043Ca'1787e—O33679e一035.087e-048270e~0442.766e-O48859e-048.201e一052941e一04bCl1560e一041.650e-041230e-041320e-046S矾6.708e-054.390e一051.129e?045161e一057CO:.'6377e一068O05e-069918一044,750e—O5编号衍生组分1HCO3021e-035987e-O33.690e一041.544e-032CaHCO5Olie一05l8:tOe一041948e—O61.27Oe—O5 3CaSO;.1261e-05l326e∞57.553e-065164e-064CaCO:"1013e-052.049e-051.250e—O44.0lOe-055CaOH8307e—O99.040e—O97.930e—O76.208e—08 6CaHS0'2.262e—l24.256e一124.256~一l55.926e-14 KSO:6830e-092.838e-086.431e∞83103e—O88KOH2.261e~129.089e-123.60te一091.946e一109Mg}{CO,7.158e一064046e-052885e--074158e一06 10MgSO.2302e一063780e-061433e一062163e一06 l1Mgc0.89~Oe-O72824e—O61147e-058120e-0612MH'2314e—084.799e一082802e—O64842e-07l3NaHC0331Be一085.360e—O73.275e一081.108e-07l4NaSO5604e—092814e—O88052e一0829t6e一08 15Na~O31986e-091919e_085882e—079,993e—O8I6NaOH.4.957e—l22407e-111.204e-084.885e-10l7HS0.1447e-101,493e一1O8565e一137.644e-12承爨组分台量(mol/L)编号主要组铮排7-1排8-2排9—1摊i4-21Na2.118e-041.338e-042.6l5e-042.794e-04K'I.499e-041.1tOe-041.349e-04t.5O9e—O43Ca5976e-049905e-046.271e一046848e一0441857e一044745e—O4l05Oe一041.422e-04Cit.720e一04I68Oe-041220e一041830e一046SO8.247e一∞4.840e-057.025e-05l435e-Og7C0:3.558e-041.910e—O3706e一045,605e-05缡号斯生组分lHCOf2565e-042】79e—O32600e-04】237e-032CaHCO:1.548e一062.085e一05I.653e一068,549e-06 3CaSG~6.146e-065465e一065538e一061,225e-054CaCO:.2.237e一041820e-052464e-044035e-05CaOH2.069e一062.O27e-082.244e-067.740e一086CaHS01.540e—l52+269e-131.345e-159398e-147KSO;6.863e-082.849e-085282e—O8l_∞le一078KOH1.215e—O85426e—l11.128e-083.992e一109HCO'4425e—O79.206e一062545e一07l634e一0610~lgSO{.2.248e-063086e—O6109le-062995e-06l1CO,3.956e-054969e一062347e-054.772e一06l2g侧1.409e一052131e-O?8.23le~063525e-0713NaHCO52713e-081.424e一073.273e-08I.726e-07l4NaSO一6.920e-082.453e—O87.o28e一081.588e—O7 15NaCO31l11e一063.606e-08I.380e-062.31Oe-O7l6NaOH3.272e一081.248e一104.008e-08l_410e-0917HSO[2.709e—l32.516e一1l2246e-131.443e-1l层吸附于介孔内表面;②在0.25到0.50中等压力阶段,吸附量急剧增加,吸附等温线急速向上翘,呈现两个明显的拐点特征,表明合成分子筛介孔孔径分布很窄II③而后在0.50N0.80的较高压力阶段,介孔内的吸附基本趋于平衡,随着相对压力的增加吸附曲线变得平缓,N1分子以单层到多层吸附于介孔外表面;滞后环的出现是由于氮气在粒子问的空穴(次级孔)内引起的毛细凝聚造成的I1.④最后由于氮气在大颗粒间的凝聚使得在高压阶段产生一个突跃.从图2(b)可以看出,分子筛具有较窄的孔径分布特征,介孔结构主要分布在2.3nm 附近.用BET方法计算所合成的高铝介孔分子筛的比表面积为300.13m2/g;由BJH方法对脱附曲线数据计算得到平均孔径为3. 36nm,平均孔容为0.25cm/g.2.3SAXS图3是由SAXS得到的介孔分子筛样品的孔径分布曲线.该结果与前面N,吸附一脱附的结果一致.3结论以天然矿物微斜长石为原料,十六烷基三甲基溴化铵(CTAB)为模板剂,白炭黑作为补充硅源,在碱性条件下水热合成了高铝含量的介孔分子筛.结果显示:所合成样品的晶胞参数为4.17nm,平均孔径为3. 36nm,平均孔容为0.25cm3/g,BET比表面积为300.13m/g.图3由SAXS得到的样品的子L径分布曲线誊上接毒计算结果的整理与分析:依据PHREEQC软件,对代表性孔灌l0—4,灌174,灌24—2,灌25一l,{{E71,{{E8—2,排9一l,排l4—2孔地下水中各组分及其衍生组分的浓度进行计算,得到如下大坝坝址区监N~I4L地下水水质组分及含量统计表l.由表l,表2可见:该大坝坝基地下水中各化学组份除主要是以单一离子或络阴离子的形式存在以外,在一定的水化学环境下也可以以络合物的形式存在.两个计算部位坝基地下水中的K和C1基本上是以单一离子(K和Cl~)的形式存在,特别是Cl,其计算浓度占分析浓度的百分比基本上是100%;而Na,Ca~1:IMg等元素除分别主要以单一离子和络阴离子形式存在外,还常以无机络合物的形式存在;至于碳元素除主要以络阴离子(HCO)形式存在外,还以CaHCO',MgHCO,NaHCO., CaCO.,MgCO.,NaCO3等络合物组份形式存在,而以CO~2-形式存在的数量占其总含量较小,界于0.1%~4.1%之间.对于所述的8个代表性水样监测孔,坝基地下水溶液中自由Caz'相应组分总量中占有绝对优势,其含量百分比均已超过了70%,坝基水中Ca的迁移,主要来自水与帷幕水泥石间的作用,这是坝基化学潜蚀的主要方式之一,也是帷幕遭受侵蚀的标志.8个代表性水样监NII~L中灌24—2孔SO一相对含量最大,达到92.5%;灌l7—4孔的自由so相对含量最小,只有72. O%.造成这种现象的原因是,以CaSQ~和MgSO.形式存在SO在含S组分总量中占了较大的比例.Ca和Mg的情况与SO类似,由于它们与CO,HCO及SO…的络合作用,使得以自由离子形式存在的Ca和Mg合量有一定程度的减少. 同样,对于矿化度较高的水样,以CaCOo 和CaHCO形式存在的CO在含C组分总量中所占的比例也达到了不可忽略的程度.总之,阴阳离子问的络合作用使得以自由离子形式存在的组分含量有不同程度的减少,自由C1,Na'和K含量减少的程度要小一些,而SO,Ca和Mg含量减少的程度则较大.对于离子的吸附特性,由于CaCO,., CaSO.,MgCO.,MgSO.等络合组份的存在,游离的Ca,Mg等离子的浓度降低,使得Ca,Mg等组份被介质吸附的能力降低,从而增强了这些组份的迁移能力,进而使得Ca,Mg碳酸盐的饱和度降低.综上所述,若仅以分析浓度去研究坝基水一岩作用进行坝基水一岩作用模拟或以此为根据分析评价坝基地下水对帷幕的侵蚀性,势必会有很大的误差,甚至得出错误的结论.所以进行坝基地下水水质组分及其络合物形式的定量计算显得有实际意义.5结论本文着重介绍了水一岩作用的反向模拟,运用PHREEQC程序对陈村大坝坝址区地下水进行了反向模拟计算,结果显示了各代表性监N:tL水质络合离子的组分含量和矿物饱和指数,并进而进行了络合离子的含量分析,饱和指数的结果分析.一一置々..,一。

pdca检验科案例

pdca检验科案例

PDCA检验科案例介绍PDCA(Plan-Do-Check-Act)是一种管理方法论,也被称为“循环质量管理”。

它的核心思想是通过循环的过程不断改进和提高工作效率和质量。

本文将通过一个科学实验的案例来探讨PDCA检验的应用。

实验背景在科学研究中,实验是验证假设和推论的重要手段。

然而,在进行实验之前,科学家需要制定实验计划,明确实验目的、方法和预期结果。

而且,在实验过程中,科学家还需要对实验数据进行检验,确保实验结果的准确性和可靠性。

实验目的本实验的目的是通过PDCA检验方法,验证一种新型化肥对植物生长的影响。

具体而言,我们想要确定该化肥是否能够提高植物的生长速度和产量。

实验步骤计划(Plan)1.确定实验目标:验证新型化肥对植物生长的影响。

2.设计实验方案:选择一种适合的植物作为实验对象,并将其分为两组,一组使用新型化肥,另一组使用传统化肥作为对照组。

每组分别设置若干个重复样本。

3.制定实验步骤:准备土壤和种子,按照实验方案进行组织实验。

确保实验条件的一致性。

4.制定数据采集计划:确定需要采集的数据,如植物的高度、叶片数量和果实产量等。

执行(Do)1.准备实验材料:购买新型化肥和传统化肥,并准备好土壤、种子和其他实验用具。

2.实施实验方案:按照计划进行实验,确保每组样本的处理方式和环境条件相同。

3.记录实验过程和数据:记录实验过程中的关键步骤和观察结果,并及时采集实验数据。

检查(Check)1.数据分析:对实验数据进行统计分析,比较使用新型化肥和传统化肥的植物生长情况。

2.结果评估:根据数据分析结果,评估新型化肥对植物生长的影响是否显著。

3.检查实验过程:回顾实验过程中的操作是否符合计划,是否存在误差或偏差。

行动(Act)1.结论和建议:根据实验结果,得出结论并提出相应的建议,如是否推广使用新型化肥。

2.改进措施:如果实验结果不理想,提出改进措施,如调整实验方案、改进实验条件等。

3.实施改进:根据改进措施,进行新一轮的实验,验证改进效果。

phreeqc中文使用手册

phreeqc中文使用手册

phreeqc中文使用手册摘要:一、引言- 介绍phreeqc 中文使用手册的目的和适用对象二、phreeqc 软件简介- 简要介绍phreeqc 软件的背景、功能和应用领域三、phreeqc 中文使用手册的结构- 详述手册中各章节的内容和作用四、phreeqc 软件安装与配置- 说明软件的安装流程和注意事项- 介绍如何进行软件配置以满足不同用户需求五、phreeqc 模型与输入参数- 详述phreeqc 软件中涉及的各种模型及作用- 解释输入参数的意义和设置方法六、phreeqc 软件操作步骤- 给出使用phreeqc 软件进行计算的具体步骤七、phreeqc 结果解读与分析- 说明如何解读和分析软件输出的结果- 提供一些实际应用案例以供参考八、phreeqc 软件的局限性与未来发展- 阐述phreeqc 软件在实际应用中可能存在的局限性- 展望软件未来发展的方向和可能的功能改进九、附录- 提供一些实用的附录,如软件操作中的常见问题解答、术语解释等正文:一、引言《phreeqc 中文使用手册》旨在帮助广大用户更好地理解和使用phreeqc 软件,为相关领域的科研和工程应用提供支持。

本手册适用于具有一定地球化学背景知识的读者,无论是初学者还是有经验的用户,都可以从中学到有关phreeqc 软件的详细信息。

二、phreeqc 软件简介phreeqc 是一款功能强大的地球化学计算软件,主要用于研究地下水、地表水、土壤等介质中的化学物质迁移转化过程。

该软件广泛应用于水资源管理、环境保护、地质勘探、农业等领域,为科研人员和工程师提供了方便的分析工具。

三、phreeqc 中文使用手册的结构本手册共分为九章,分别为:引言、phreeqc 软件简介、phreeqc 中文使用手册的结构、phreeqc 软件安装与配置、phreeqc 模型与输入参数、phreeqc 软件操作步骤、phreeqc 结果解读与分析、phreeqc 软件的局限性与未来发展以及附录。

Phreeqc教程2 -热力学数据库.

Phreeqc教程2 -热力学数据库.
qrxn = - (cs solution J/g-K) (mass of solution g) (DT K)
= - (4.18 J/g-K) [(1.0g/mL)(100 mL)] (6.5 K)
= - 2700 J or 2.7 kJ
DH = 2.7 kJ
Enthalpy change per mole of HCl = (-2.7 kJ)/(0.050 mol)
– HS- + Hg = HgS + H+ + 2e-
– H2S + Hg = HgS + 2 H+ + 2e-
E
E
0
RT nF
ln
HS H
– Where E is Ep from the analytical peak, E0 is the formal potential, n is # e-’s, F is Faraday’s constant
Determining Entropy
• As for H, dS=q/dT can be measured as heat energy (q)
• Another way to think of entropic energy – for any reaction, energy is ‘dispersed’ to/from the surroundings –measured from 0K (actually just close to it), where S0=0 for ANY substance (at 0 K, atoms do not MOVE!)
• There can be significant error in the thermodynamic data used in different databases.

PHREEQC在地下水溶质反应运移模拟中的应用

PHREEQC在地下水溶质反应运移模拟中的应用
考虑 到 动 态 生 物 和 化 学 反 应 的 复 杂 多 样 性 , PHREEQC 准许用户用 Basic 语言编写相应的动态反应 语句直接嵌入输入文件 (或数据库) 。PHREEQC 采用 了龙格 - 库塔 ( Runge2Kutta) 法 ,通过反应速度在时间 上的积分来模拟动态反应过程 。
助项目 (50178040) 作者简介 : 毛晓敏 (19712) ,女 ,博士后 ,主要从事水文水资源及
水环境方面的研究 。 E2mail :xmao @ed. ac. uk
移结合起来 ,发展多学科交叉的地下水溶质反应运移 模型[5] 。
本文利用 PHREEQC 对地下水运动情况下土壤中 溶质交换和运移 、氧化还原反应和生物降解分别进行 了模拟 ,并与有关资料进行了验证分析 。
液进 行 淋 洗 交 换 ( 溶 液 浓 度 分 别 为 012N 、210N 和
210N) ,同时采集渗出液观察含标记钙离子溶液浓度的
变化 。
(2) 主要反应概述 因为在两种钙离子存在的条
件下 ,土壤溶液中有如下离子交换反应 :
Ca2+ + Ca 3 - X2
Ca - X2 + Ca 3 2+
(1)
· 2 0 ·
水文地质工程地质
2004 年第 2 期
PHREEQC 在地下水溶质反应2运移模拟中的应用
毛晓敏1 ,刘 翔2 ,Barry D A1 (11 英国爱丁堡大学土木与环境研究所 ,爱丁堡 EH9 3JN ;
21 清华大学环境科学与工程系 ,北京 100084)
摘要 : 由于地下水污染的加剧 ,对地下水中污染物运移规律的研究日益受到重视 。地下水中的溶质在运移过程中伴随 着溶质组分间的化学反应 ,因此需要建立地下水溶质运移与化学反应的耦合模型 。PHREEQC 是近年来发展起来的描述 局部平衡反应 、动态生物化学反应的水文地球化学模拟软件 。本文利用该模拟软件对一维地下水流动过程中溶质离子 交换反应和动态氧化还原反应进行了模拟 。结果表明 , PHREEQC 能够成功地进行溶质运移情况下复杂水化学反应模 拟 ,但对于复杂地下水流和溶质运动的情况 ,有必要耦合其它的地下水流动和溶质运移软件来共同完成 。 关键词 : 地下水 ; 溶质 ; 反应运移 ; 离子交换 ; 动态氧化还原 ; PHREEQC 中图分类号 : P64112 ;P64113 文献标识码 : A 文章编号 : 100023665 (2004) 0220020205

phreeqc中文使用手册

phreeqc中文使用手册

phreeqc中文使用手册摘要:1.引言2.phreeqc 软件的概述3.phreeqc 软件的使用方法a.安装和设置b.输入数据c.运行模型d.输出结果4.phreeqc 软件的应用领域5.phreeqc 软件的优缺点分析6.结论正文:【引言】phreeqc 是一款中文使用手册,为用户提供了关于phreeqc 软件的使用方法和详细说明。

本文将重点介绍phreeqc 软件的概述、使用方法、应用领域、优缺点分析等内容,帮助用户更好地理解和应用该软件。

【phreeqc 软件的概述】phreeqc 是一款基于化学反应原理的地球化学模型软件,广泛应用于水文地质、环境地质、土壤改良等领域。

通过模拟地下水、地表水、土壤等介质中化学物质的迁移、转化等过程,为用户提供地球化学环境问题的解决方案。

【phreeqc 软件的使用方法】使用phreeqc 软件需要进行以下步骤:a.安装和设置:根据操作系统的不同,下载相应的软件安装包并进行安装。

在安装过程中,用户需要设置一些基本参数,如模型类型、单位制等。

b.输入数据:用户需要输入关于地质、水文、化学物质等方面的数据,包括各类参数、边界条件等。

数据可以来自于实验、观测或文献等。

c.运行模型:在完成数据输入后,用户可以启动模型运行。

软件会根据设定的模型类型和参数进行计算,得出模拟结果。

d.输出结果:软件会生成各种结果文件,包括数据表格、图形等。

用户可以根据需要查看和分析这些结果,以了解地球化学过程和环境问题。

【phreeqc 软件的应用领域】phreeqc 软件广泛应用于以下领域:1.水文地质:研究地下水、地表水等地质水体的化学特征、物质循环等过程。

2.环境地质:评估污染物的迁移、转化、归宿等过程,为污染治理提供依据。

3.土壤改良:研究土壤中化学物质的变化规律,为土壤改良和农业生产提供指导。

4.其他领域:还包括地质资源勘查、地质灾害预测等。

【phreeqc 软件的优缺点分析】优点:1.强大的模拟能力:phreeqc 软件能够模拟复杂的地球化学过程,为用户提供准确的模拟结果。

  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

PHREEQC实例分析例2——矿物相的溶解平衡这个例子测定了最稳定相——石膏或是无水石膏在一定温度围的溶解性。

输入数据组见表13。

仅用pH和温度来定义纯水溶液。

缺省单位是millimolal,但没有指定浓度。

缺省状态下,pe为4.0,缺省的氧化还原计算使用pe,且水的密度为1.0(这是没有必要的,因为没有浓度为“每升”)。

在批反应期间,所有允许反应达到指定的饱和指数的反应都列在EQUILIBRIUM_PHASES中,无论开始时它们是否存在。

输入数据包含相的名字(之前通过PHASES输入在数据库或输入文件中定义)、指定的饱和指数和以摩尔数表示的当前存在的相的数量。

如果一种相在当前不存在,则在纯相集合中给定其浓度为0.0 mol。

在这个例子中,石膏和无水石膏允许反应来达到平衡(饱和指数等于0.0),初始相的集合中每种矿物都为1mol。

每一种矿物或是反应达到平衡或是被耗尽。

在大多数的情况下,1 mol的单纯相足够达到反应平衡。

表13--例2输入数据的设置TITLE Example 2.--Temperature dependence of solubilityof gypsum and anhydriteSOLUTION 1 Pure waterpH 7.0temp 25.0EQUILIBRIUM_PHASES 1Gypsum 0.0 1.0Anhydrite 0.0 1.0REACTION_TEMPERATURE 125.0 75.0 in 51 stepsSELECTED_OUTPUTfile ex2.selsi anhydrite gypsumEND在REACTION_TEMPERATURE数据块中,指定反应温度计算步长为1℃,从25℃开始,75℃结束,共计算温度变化51次。

温度的每一度变化都详细指定输入数据,在可能的情况下,在EQUILIBRIUM_PHASES(石膏和无水石膏)中定义的物质将会反应达到平衡,或是直到两种物质都完全溶解。

最后在每一次计算之后,用SELECTED_OUTPUT来输出石膏和无水石膏的饱和指数到文件ex2.sel中。

然后用这个文件来生成图5。

初始溶液计算和第一次批反应的结果见表14。

纯水中物质的分布列于标题“开始初始溶液计算”下。

在25℃下,与所给定石膏和无水石膏量达到平衡的体系是第一次批反应,反应结果列在标题“开始批反应计算”之后。

紧接着这个标题,定义了批反应阶段,接下来是计算中确定关键字数据的列表。

在这个例子中,计算时所应用的反应组分保存为序号1,纯相集合保存为序号1,反应温度保存为序号1。

概念上,溶液和纯相放在一个烧杯中,并调整温度为25℃,使其反应达到体系平衡。

图5--在25到75℃的温度围,溶液中石膏和硬石膏稳定时的饱和指数,在子标题“Phase assemblage”下,列出了由EQUILIBRIUM_PHASES所定义的每一种相的饱和指数和数量。

在第一次的批反应阶段,最后相的集合不含无水石膏,对溶液而言这并没有达到平衡(饱和指数为-0.22),1.985mol的石膏,与溶液达到平衡(饱和指数等于0.0)。

所有的无水石膏得以溶解,且大多数的钙和硫酸盐是生成石膏沉淀下来。

“溶液组分”表明,15.64 mmol/kgw的钙和硫酸盐留在溶液中,这确定了纯水中石膏的溶解度。

然而,在液相中每一成分的摩尔总数为15.08,这是因为水的质量仅为0.9645kg(“溶液的定义”)。

在形成石膏(CaSO4 .2H2O)时,水从溶液中移除。

因此,在批反应计算中,溶剂水的质量并不是常数;在溶解和沉淀相中,反应和水合作用可以增加或减少溶剂水的质量。

所有批反应阶段的饱和指数绘于图5。

在每一阶段,纯水与相在不同的温度下反应(反应并不累加)。

PHREEQC缺省的数据块表明:在温度低于57℃时,石膏是稳定相(饱和指数等于0.0);在高于这个温度时,无水石膏是稳定相。

表14—例2的选择性输出-------------------------------------------Beginning of initial solution calculations.-------------------------------------------Initial solution 1. Pure water-----------------------------Solutioncomposition------------------------------Elements Molality MolesPure water----------------------------Description of solution----------------------------pH = 7.000pe = 4.000Activity of water = 1.000Ionic strength = 1.001e-07Mass of water (kg) = 1.000e+00Total alkalinity (eq/kg) = 1.082e-10Total carbon (mol/kg) = 0.000e+00Total CO2 (mol/kg) = 0.000e+00Temperature (deg C) = 25.000Electrical balance (eq) = -1.082e-10Percent error, 100*(Cat-|An|)/(Cat+|An|) = -0.05Iterations = 0Total H = 1.110124e+02Total O = 5.550622e+01----------------------------Distribution of species----------------------------Log Log Log Species Molality Activity Molality Activity GammaOH- 1.002e-07 1.001e-07 -6.999 -6.999 -0.000 H+ 1.001e-07 1.000e-07 -7.000 -7.000 -0.000 H2O 5.551e+01 1.000e+00 0.000 0.000 0.000H(0) 1.416e-25H2 7.079e-26 7.079e-26 -25.150 -25.150 0.000O(0) 0.000e+00O2 0.000e+00 0.000e+00 -42.080 -42.080 0.000------------------------------Saturationindices-------------------------------Phase SI log IAP log KTH2(g) -22.00 -22.00 0.00 H2O2(g) -39.12 44.00 83.12 O2-----------------------------------------Beginning of batch-reaction calculations.-----------------------------------------Reaction step 1.Using solution 1. Pure waterUsing pure phase assemblage 1.Using temperature 1.-------------------------------Phaseassemblage--------------------------------Moles in assemblage Phase SI log IAP log KT Initial Final DeltaAnhydrite -0.22 -4.58 -4.36 1.000e+00 -1.000e+00 Gypsum 0.00 -4.58 -4.58 1.000e+00 1.985e+00 9.849e-01-----------------------------Solutioncomposition------------------------------Elements Molality MolesCa 1.564e-02 1.508e-02S 1.564e-02 1.508e-02----------------------------Description of solution----------------------------pH = 7.062 Charge balance pe = 10.691 Adjusted to redox equilibriumActivity of water = 1.000Ionic strength = 4.178e-02Mass of water (kg) = 9.645e-01Total alkalinity (eq/kg) = 1.122e-10Total carbon (mol/kg) = 0.000e+00Total CO2 (mol/kg) = 0.000e+00Temperature (deg C) = 25.000Electrical balance (eq) = -9.766e-11Percent error, 100*(Cat-|An|)/(Cat+|An|) = -0.00Iterations = 23Total H = 1.070728e+02Total O = 5.359671e+01----------------------------Distribution of species----------------------------Log Log Log Species Molality Activity Molality Activity GammaOH- 1.402e-07 1.155e-07 -6.853 -6.937 -0.084 H+ 1.006e-07 8.665e-08 -6.997 -7.062 -0.065 H2O 5.551e+01 9.996e-01 -0.000 -0.000 0.000 Ca 1.564e-02Ca+2 1.045e-02 5.176e-03 -1.981 -2.286 -0.305 CaSO4 5.191e-03 5.242e-03 -2.285 -2.281 0.004 CaOH+ 1.192e-08 9.910e-09 -7.924 -8.004 -0.080H(0) 4.363e-39H2 2.181e-39 2.202e-39 -38.661 -38.657 0.004O(0) 1.701e-15O2 8.504e-16 8.587e-16 -15.070 -15.066 0.004S(-2) 0.000e+00HS- 0.000e+00 0.000e+00 -117.650 -117.734 -0.084 H2S 0.000e+00 0.000e+00 -117.807 -117.803 0.004 S-2 0.000e+00 0.000e+00 -123.278 -123.590 -0.312S(6) 1.564e-02SO4-2 1.045e-02 5.075e-03 -1.981 -2.295 -0.313 CaSO4 5.191e-03 5.242e-03 -2.285 -2.281 0.004 HSO4- 5.141e-08 4.276e-08 -7.289 -7.369 -0.080------------------------------Saturationindices-------------------------------Phase SI log IAP log KTAnhydrite -0.22 -4.58 -4.36 CaSO4Gypsum 0.00 -4.58 -4.58 CaSO4:2H2OH2(g) -35.51 -35.51 0.00 H2H2S(g) -116.81 -158.45 -41.64 H2SO2(g) -12.11 71.01 83.12 O2Sulfur -87.18 -122.94 -35.76 S例3——混合反应这个例子证明了PHREEQC进行一系列地球化学模拟的能力,在同一次运行中,最后的模拟依赖于以前模拟的结果。

相关文档
最新文档