Pitzer活度系数模型研究与开发
状态方程和活度系数模型
状态方程和活度系数模型1. 引言在化学领域，我们经常需要描述和预测化学反应的行为和性质。
为了达到这一目的，我们需要建立一些数学模型来描述不同物质在不同条件下的行为。
状态方程和活度系数模型是其中两个重要的模型，它们在研究溶液中的物质浓度和反应速率等方面发挥着重要的作用。
2. 状态方程状态方程是一种用数学语言描述物理和化学系统状态的模型。
在化学中，最常见的状态方程是气体状态方程，如理想气体状态方程和范德瓦尔斯状态方程。
理想气体状态方程可以表示为：[ PV = nRT ]其中，P表示气体的压力，V表示气体的体积，n表示气体的物质量，R表示气体常数，T表示气体的温度。
理想气体状态方程假设气体分子之间没有相互作用，体积可忽略不计。
但在一些情况下，气体分子之间会存在相互作用，此时可以使用范德瓦尔斯状态方程：[ (P + )(V - nb) = nRT ]其中，a和b为范德瓦尔斯常数，分别代表相互作用带来的修正。
范德瓦尔斯状态方程适用于高压和低温条件下的气体。
除了气体状态方程，还有液体和固体的状态方程，如液体的状态方程可以由万有气体状态方程经过适当的修正得到：[ P = P^* + ]其中，P^*为理想液体的气体状态方程。
同理，固体的状态方程也可以类似进行修正。
3. 活度系数模型在溶液中，溶质分子和溶剂分子之间会发生相互作用。
这些相互作用会导致溶质分子在溶液中的活动程度与其在纯溶剂中的活动程度不同。
为了描述溶液中物质的活性，引入了活度系数的概念。
活度系数（activity coefficient）描述了溶液中溶质活性和摩尔浓度之间的关系。
活度是指溶液中溶质的实际浓度和理想溶液中相同浓度之间的比值。
活度系数模型是一种用于计算活度系数的数学模型。
常用的活度系数模型有Debye-Huckel模型和Pitzer模型。
Debye-Huckel模型适用于稀溶液中的弱电解质。
该模型通过引入电离度和离子间相互作用来计算活度系数。
海洋物理化学思考题1
海洋物理化学思考题第一章液态水的结构1、液态水有哪两种结构模型一混合型理论模型间隙模型混合模型闪动簇团模型二连续体模型2、“闪动簇团”模型有什么特点？不必假设冰I h结晶的存在，也不必考虑簇团中水分子构成的方式，因而避开了过冷、结冰与流动性等现象所提出的问题“自由”水中的氢键已破裂，分子之间力为偶极-偶极作用和L o n d o n作用簇团具有闪动的特性3、用“闪动簇团”模型解释压力和温度对纯水结构的影响。
压力：簇团理论：簇团由于具有一定的构型而比它周围“自由”的水体积更大，压力增大，将破坏这些簇团使体积减小温度：：簇团的破坏所引起的密度增加和正常热膨胀所引起的密度减小两种因素综合作用的结果。
4、什么是电缩作用？离子与水偶极的局部电荷的相互静电作用把水分子拉近离子周围，使靠近离子附近水的比容变小，即离子周围水的体积缩小了，这种现象称为电缩作用。
5、水的双区模型有几种？6、什么是溶剂化配位水？7、什么是非溶剂化配位水？第二章海水状态方程式1.什么是海水状态方程式2.海水状态方程式可能的用途3.UNESCO高压海水状态方程式的形式4.UNESCO高压海水状态方程式有什么优点？第三章海水热力学基础与盐度、离子强度的关系；离子强度与盐度的关系。
1、ms2、偏摩尔体积与表观摩尔体积的定义，它们之间的关系。
第四章离子水化作用1、根据连续电解质模型，推导水化焓公式2、对连续电解质模型的改进工作有哪些？3、在结构水化模型中，将离子与水的相互作用分成了几个步骤？最后得到的水化焓公式是什么？第五章海水体系中离子－离子相互作用1、离子雰的概念，离子互吸理论的假设。
由于静电作用力的影响，在中心离子（正离子）周围，距离正离子愈近的地方，正电荷密度愈小，负电荷密度愈大，结果在中心周围大部分的正负电荷相互抵消，但总效果犹如在其周围分布着一个大小相等而符号相反的电荷，我们把这一层电荷所构成的球体称为离子雰。
2、离子雰在中心离子J上引出的电能变化是什么，并据此推导D-H活度系数公式。
化工模拟软件介绍
Aspen Plus y 7kt a%vdbwx;O)jI WAspen Plus是一种广泛应用于化工过程的研究开发，设计，生产过程的控制，优化及技术改造等方面的性能优良的软件。
该模拟系统是麻省理工学院于70年代后期研制开发的。
由美国AspenTech公司80年代初推向市场，它用严格和最新的计算方法，进行单元和全过程的计算，为企业提供准确的单元操作模型，还可以评估已有装置的优化操作或新建，改建装置的优化设计。
这套系统功能齐全，规模庞大，可应用于化工，炼油，石油化工，气体加工，煤炭，医药，冶金，环境保护，动力，节能，食品等许多工业领域。
}Q .aU U "q D?IPRO/II w w6A/6 z'&bo !jPRO/II流程模拟程序广泛地应用于化学过程的严格的质量和能量平衡。
从油/气分离到反应精馏，PRO/II提供了最广泛的、最容易使用有效模拟工具。
产品的PROVISION图形用户界面GUI), 提供了一个完全交互的、基于Windows的环境，无论是对于建立简单的，还是复杂的PRO/II模型，它都是理想的环境。
%d.Fn@=bh, n( ewDChemCAD 9o#,P 3 <2 Y 2s' wChemCAD是一个用于对化学和石油工业、炼油、油气加工等领域中的工艺过程进行计算机模拟的应用软件，是工程技术人员用来对连续操作单元进行物料平衡和能量平衡核算的有力工具。
使用它。
可以在计算机上建立和现场装置吻合的数据模型，并通过运算模拟装置的稳态和动态运行，为工艺开发、工程设计以及优化操作提供理论指导。
5 Wl [?6:S2B> WW0KHYSYS D do+dQ tGz)J",Q}?r+oI 3 fPHyprotech出品的Hysys是一个化工流程模拟动态仿真软件，是一款环境模拟设计软件，允许设计者通过概念上的设计而简化制作过程来完成项目工作。
完整的交互性能充分发挥你的创造力。
活度系数计算
电解质溶液活度计算理论进展【摘要】：由于溶液大多数不是理想溶液，需要用活度来代替浓度。
活度系数又是描述活度与浓度的差异程度，因此活度系数的计算对于反应过程相当的重要。
近几年，随着活度系数理论模型的不断发展，活度系数的计算方法也在不断的提高、创新。
本文在回顾电解质溶液热力学经典理论的基础上，对活度系数计算做了综述。
【关键词】：活度系数活度模型热力学模型活度计算Electrolyte solution activity in recent years, progressin computational theoryAbstract：Solution is not ideal because most of the solution need to replace the concentration of activity. Activity coefficient is described differences in degree of activity and concentration, so the calculation of activity coefficients for the reaction process was very important. In recent years, with the activity coefficient of the continuous development of theoretical models, the calculation of activity coefficients are also constantly improving and innovation. In this paper, recalling the classical theory of thermodynamics of electrolyte solution, based on calculations made on the activity coefficient is reviewed.Keywords: Activity coefficient, Activity Model, Thermodynamic model, Activity calculation1、活度与活度系数绝大多数的反应都有溶液（固溶体、冶金熔体及水溶液）参加，而这些溶液经常都不是理想溶液，在进行定量的热力学计算和分析，溶液中各组分的浓度必须代以活度。
K2B.07的Pitzer参数的计算
李丽娟
姚燕
宋彭生 西宁810008)
(中国科学院青海盐湖研究所
摘要:由于硼酸盐在电解质中的特殊性,硼酸盐的热力学问题一直被化学家所重视.本文 根据文献数据,经多种处理获得K:B407不同存在形式的硼氧配阴离子的Pitzer纯盐参数, K2B.O,(OH).电离化学计量平衡常数K..采用所获得的参数计算K,B.0,不同存在形式的 硼氧配阴离子的渗透系数,并与实验值进行比较. 关键词:Pitzer参数,
..(3)
∑扰,--mK+mB
K28407+2 H20§2K++[B405(OHh].
mz=2mB;mKB2mo·mB;123mn
.(4)
4mB3+KmmB—l(nImo=0
∑肋-----Illx+me
l
K28407+2
H20=2K*+【B405(OHh]2'..
.(5)
K++【B40,(0H)'】-§tKS'o,(OH)4】.
1.1215,c(^)0.06399,SD
O.0077.这套数据是采用R.F.Plafford[+l的KzB407溶液25"(2渗透系
数和对应的活度系数一起拟合而来的. 本文采用Pitzer离子相互作用模型,对该文献【61的实验数据进行了理论分析.
2
Pitzer参数计算依据 我们按照Pitzer方程嘲,实验渗透系数计算公式以及不同类型硼氧配阴离子浓度计算公
Abstract
Key WOrds:Pitzer Parameters;K2840,;poly borate
anion
用离子选择性电极测定混合电解质溶液中盐的活度系数实验报告
在诸多计算电解质溶液活度系数的公式中，Pitzer 公式是目前应用最广，精度最高的一
种。对于每一种电解质溶液都有三个参数M������X。对于 KCl，计算公式为：
KCl（分析纯），NaCl（优级纯），去离子水 2） 测量仪器
电子天平，离子计，K+离子选择性电极，Cl-离子选择性电极，电极测量车，磁力搅拌 器，恒温槽。
3） 装置图
4） 5） Fig.1-1 电解质溶液中活度系数的测定装置
6） 1-超级数显恒温槽；2-磁力搅拌器；3-离子活度计；4-电极测量池；5-离子选择性电极
将 KCl 活度系数带入下式：������������������������ = (������±������������������������������������������)2
计算出活度系数������������������������，其中������������������������为 KCl的质量摩尔浓度；������±������������������为 KCl的活度系数。作出
式中条件为 T=298.15K，n=1。 2） KCl-NaCl 混合溶液测定
组成如下无液接电势的电池：������+ ISE | KCl − NaClaq | Cl− ISE。 其电位值为：������������������������ = ������������������������������ + ������������������������������������������������������������ 其中：
Aspen流程模拟基础(入门)-1汇总
ASPEN PLUS模拟固体系统
Aspen Plus在煤的净化和液化、流化床燃烧、高温冶金和湿法冶金， 以及固体废物、聚合物、生物和食品加工业中都得到了应用。
Aspen Plus中固体性质数据有两个来源：一是Solid数据库，它广泛收 集了约3314种纯无机和有机物质的热化学数据；二是和CSIRO数据库 的接口。第10版已合并成新的SOLIDS数据库，这些数据在模拟冶金、 陶瓷矿产业及其它含有固体处理的过程时是必不可少的。还具有一套 通用的处理固体的单元操作模型，包括破碎机、旋风分离器、筛分、 文杜里洗涤器、静电沉淀器、过滤洗涤机和倾析器。此外，Aspen Plus中所有的单元操作都适合于处理固体，例如闪蒸和加热器模型能 计算固体的能量平衡，而反应器模型 RGIBBS可用最小GIBBS自由能来 判断在平衡状态下是否有固相存在。
Aspen Plus还有经过工业考验的能处理反应的分离模型，该模型可在塔 的任意塔板处或所有塔板上处理速率控制反应、化学平衡反应，以及 气、液相反应。反应速率可由置入内部的幂律表示式或由用户提供的 反应动力学程序来计算 。Aspen Plus的简捷算法蒸馏模型需要输入的数 据较少，也具有设计和核算两种型式。在不需要高度精确计算的情况 下可以使用这些模型。
ASPEN PLUห้องสมุดไป่ตู้具有完整的单元操作模型库
Aspen Plus有一套完整的单元操作模型，可以模拟各种操作过程,由单 个原油蒸馏塔的计算到整个合成氨厂的模拟。
由于Aspen Plus系统采用了先进的PLEX数据结构，对于组分数、进 出口物流数、塔的理论板数以及反应数目均无限制，这是Aspen Plus 的一项独特优点，非其它过程模拟软件所能比拟。
第五章 流体状态方程
只包含了纯物质的临界参数，计算方便，在宽广的温度和 压力范围内对非对称烃类体系（C ≤ 16）高压气液平衡具 有足够的精度。 Kordas, A., K. Magoulas, et al. (1995). "Methane-hydrocarbon interaction parameters correlation for the Peng-Robinson and the t-mPR equation of state." Fluid Phase Equilibria 112(1): 33-44.
பைடு நூலகம்
（3）非极性分子对间的总位能关系利用LJ位能关系模型
u(r ) 4 r
6 r 12
LJ位能式中的第一项为排斥作用，第二项为色散能， 对不相同分子和相同分子间的色散能而言，LJ式中的 第二项分布等于式（5-41）、式（5-42）、式（5-43） 的右端，消除极化率得到：
混合规则（Mixing Rule）： 混合物的虚拟参数与混合物的组成和纯物质的参数 之间的关系式； 限于对分子间作用力的认识尚不充分，混合规则很 难从统计热力学严格导出，因此目前的状态方程混 合规则本质上还是经验的； 混合规则的研究非常活跃！
一、经典混合规则（Classical Mixing Rule） 经典混合规则属于van der waals单流体模型， 单流体理论是把由两种或两种以上的物质组成的混 合物作为一种虚拟的纯物质，用一套特征参数来表 征，对任一热力学性质M为：
van氏或rk方程种的a和b随着计算精度的提高方程式中的参数也随之而增加在实际应用时往往会遇到查不到物质特性参数的困难为此需要研究一种基本不包含物质特性参数的新方程即普遍化的状态方程而对比态原理是普遍化状态方程的理论依据
实验3_离子选择性电极测定混合电解质溶液中离子活度系数
化工实验报告实验名称 离子选择性电极测定混合电解质溶液中离子活度系数 班级 化21 姓名 张腾 学号 2012011864 成绩 实验时间 2014.12.3 同组成员 陆言 一、实验预习 1、实验目的（1）学习用离子选择性电极进行电解质溶液中离子活度系数测定的基本原理和方法； （2）使用Pitzer 公式计算混合电解质溶液中盐离子活度系数； （3）实验测定KCl-NaCl 混合溶液中KCl 的活度系数并与理论值比较。
2、实验原理离子选择电极Ionic selective electrode ISE （，）是一种电化学敏感器，具有将溶液中某种特定离子的活度转换成一定电势的能力，其电势与特定离子活度的对数值之间具有线性关系，此关系符合能斯特Nernst （）方程：ln E E k a θ=+式中：E 为电池的电动势，mV ；E θ为标准状态下电池的电动势；k 为电极响应斜率；a 为电解质的活度。
可以使用离子选择性电极来测量电解质溶液的活度系数。
本实验采用氯离子选择电极和钾离子选择电极先测定KCl 单一溶液的电位值，得到电极参数后，再测定KCl 在KCl-NaCl 混合溶液中的电位值，得到活度系数，并与Pitzer 公式计算的理论值比较（1）单一KCl 溶液的测定用单一电解质2KCl-H O 体系组成如下无液接电池：K +ISE KCl 水溶液 Cl −ISE其电位值为：KCl KCl KCl KCl ln E E k a θ=+将KCl 的活度系数KCl γ±代入下式：2KCl KCl KCl =()a m γ±计算出活度a KCl 。
其中m KCl 为KCl 的质量摩尔浓度；γ±KCl 为KCl 的活度系数。
作E KCl −ln a KCl 图，应得到一条直线，从图中求电极参数E KCl ⊝和k KCl 。
理论上讲，电极响应斜率k 可计算得到：/k RT nF =。
ASPENPLUS入门教程(煤气化)
回答“如果…那会怎么样”的问题 在给定的限制内优化工艺条件 辅助确定一个工艺约束部位
第5页
ASPEN PLUS 基本概念
流程图（Flowsheet）
直观形象地表示所模拟系统的流程
模型库（Model Library） 存放可用单元操作模型的库
数据浏览器（Data Browser） 流股（Stream）
第18页
主要内容
ASPEN PLUS简介 ASPEN PLUS 安装方法及界面介绍 通过实例介绍如何建立模拟模型 模型分析工具的使用
第19页
ASPEN PLUS 2006 安装方法
一．运行虚拟光驱DAEMON TOOLS，载入ASPEN 2006的ISO文件 （2.61G那个）
二．运行虚拟盘上的Setup，到添加license这一步
第14页
ASPEN PLUS的优势
可以模拟电解质系统
许多公司已经用Aspen Plus模拟电解质过程，如酸水汽提、苛性 盐水结晶与蒸发、硝酸生产、湿法冶金、胺净化气体和盐酸回收 等。
Aspen Plus提供Pitzer活度系数模型和陈氏模型计算物质的活度系 数，包括强弱电解质、盐类和含有机化合物的电解质系统。这些 模型已广泛地在工业中应用，计算结果准确可靠。
压力改变
•PUMP 泵/料浆泵
•COMPR 单级压缩/膨胀 机 •MCOMPR多级压缩/膨 胀机
•PIPELINE多段管线压 降
•PIPE 单段管线压降
•VALVE 阀压降
第8页
单元操作模型及其主要功能
流控制器
•MULT 乘法器 •DUPL复制器 •CLCHNG 流股复类器 •SELECT 物流选择器 •ANALYZER 物流分析器 •QTVEC 热负荷控制器 •MEASUREMENT测量器
第12讲之1活度和活度系数
离子缔合体；离子对；离子氛
3
用离子强度 I 表示离子与“离子氛”之间作用 的强弱，zi表示溶液中 i 种离子的电荷数，Ci表 示 i 种离子的质量摩尔浓度，则
1 2 I ci zi 2
Debye和Hü ckle 采用静电学Poisson公式和分子运 动论的Boltzmann公式来求活度系数。用 ±式表 示正负离子的平均活度，其最简式为：
2
强电解质溶液理论
1923年，Debye和Hü ckle 提出了强电解质溶 液理论 （1）离子氛和离子强度
强电解质在水溶液中 是完全电离的，但离子并 不是自由的，存在着“离 子氛”(ionic atmosphere) 。 Debye和Hü ckle 将中心 离子周围的异性离子群叫 做离子氛。导致溶液导电 能力下降、电离度下降、 依数性异常。
活度和活度系数
1907年Lewis提出有效浓度的概念。他认为非理想溶 液不符合Raoult定律是因为溶质与溶剂之间有相当复 杂的作用。在没有弄清这些作用之前，可利用实验 数据对实际浓度（c）加以校正，即为有效浓度。 Lewis命名它为活度(activity)，用符号“a” 表示；校 正因子叫活度系数(activity coefficient)，常用符号 表示。
a=c
1
a=c
c ：真实浓度 a ：活度或有效浓度
：活度系数
（亦用 f 表示）
c 越大，离子电荷数越高，的数值越小，a 与
c 的偏离越大。
c 很低， 趋近于1.0 , 可用 c 代替 a 。
Debye、Huckel、Pitzer等科学家对活度系数的理 论根据做了深入研究。提出了强电解质理论。
5
lg 0.509z1 z2 I / m
活度系数计算
电解质溶液活度计算理论进展【摘要】：由于溶液大多数不是理想溶液，需要用活度来代替浓度。
活度系数又是描述活度与浓度的差异程度，因此活度系数的计算对于反应过程相当的重要。
近几年，随着活度系数理论模型的不断发展，活度系数的计算方法也在不断的提高、创新。
本文在回顾电解质溶液热力学经典理论的基础上，对活度系数计算做了综述。
【关键词】：活度系数活度模型热力学模型活度计算Electrolyte solution activity in recent years, progressin computational theoryAbstract：Solution is not ideal because most of the solution need to replace the concentration of activity. Activity coefficient is described differences in degree of activity and concentration, so the calculation of activity coefficients for the reaction process was very important. In recent years, with the activity coefficient of the continuous development of theoretical models, the calculation of activity coefficients are also constantly improving and innovation. In this paper, recalling the classical theory of thermodynamics of electrolyte solution, based on calculations made on the activity coefficient is reviewed.Keywords: Activity coefficient, Activity Model, Thermodynamic model, Activity calculation1、活度与活度系数绝大多数的反应都有溶液（固溶体、冶金熔体及水溶液）参加，而这些溶液经常都不是理想溶液，在进行定量的热力学计算和分析，溶液中各组分的浓度必须代以活度。
【国家自然科学基金】_pitzer理论_基金支持热词逐年推荐_【万方软件创新助手】_20140729
推荐指数 2 2 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
2010年 序号 பைடு நூலகம் 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20
科研热词 相平衡 高寒地区 过碳酸钠 过氧化氢 表面张力 节能技术 盐水体系 环戊烷-甲烷水合物 溶解度计算 溶解度 湖盐开发 活度系数 气体水合物 模型 新型线性预测方程 半理想溶液理论 so42-平衡浓度 pitzer理论 pitzer butler方程
科研热词 硫酸镁 硫酸铵 硫酸钠 相平衡 pitzer模型
推荐指数 2 2 2 2 2
2014年 序号 1 2 3 4 5
2014年 科研热词 湿度计 渗透压 活度 冰点渗透压仪 pitzer模型 推荐指数 1 1 1 1 1
2009年 序号 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18
科研热词 热力学 pitzer电解质溶液理论 隐钾锰型水合二氧化锰 锂离子筛 钾离子筛 量子力学模拟 相平衡 盐度 甲烷水合物 甘氨酸 活度系数 氯化钙 形成条件 多孔介质 从头算势能 丙氨酸 lini0.05mn1.95o4 gibbs自由能
科研热词 推荐指数 相图 2 溶解度 2 水盐体系 2 pitzer参数 2 镁 1 重结晶 1 相平衡 1 氟化 1 晶种 1 废水 1 分配比 1 pitzer理论 1 pitzer模型 1 pitzer方程 1 mnso4溶液 1 caso4-ca(oh)2-h2o体系 1
2013年 序号 1 2 3 4 5
2008年 序号 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18
E-07012-电解质水溶液密度和标准沸点计算
2
电解质水溶液的标准沸点
由标准沸点的基本理论[7]可以推导出计算电解质水溶液标准沸点的公式：
Tb =
* Tb*Δ vap H m * Δ vap H m + RTb* ln a
（14）
在方程（14）中， Tb* 是溶剂水的沸点， Δ vap H m 是水的摩尔蒸发焓， R 是气体常数， a 是
*
溶剂水的活度，其计算公式为：
（4）
在方程（4）中，NA为阿伏加得罗常数， d w 为溶剂水的密度，其单位是g/ml， D 是水的介 电常数； e为电荷单位，其值为 4.802×10
16
-1 －10
esu， κ 为Boltzmann常数，其值为 1.38×10
－
erg．K ；υ = υ M + υ X ，υ M 和 υ X 分别是 1mol电解质电离出的正离子和负离子的物质的量
点的方程式。方程式对样品 1 密度和渗透系数计算值的平均偏差分别为 0.08％、0.84％，利用方程式计 算了样品 1 的标准沸点。对于氯化钠水溶液渗透系数和标准沸点计算值的偏差分别为 0.52％和 1.62 K。 关键词：电解质水溶液 表观摩尔体积 标准沸点 渗透系数
在能源、化工等领域，电解质水溶液有着广泛的应用，必须有相应的热物性数据与之相 对应。因此，本文利用电解质溶液的理论[2]，采用了Pitzer（1973）[3-4]模型，计算电解质水 溶液的密度、渗透系数和标准沸点。
4m (0) 4m 2 φ (1) φ −1 = − + β MX + β MX exp(−2 3m ) + C MX 3 3 1 + 1.2 3m
2 Aφ 3m
（16）
(0) (1) φ 方程（16）中的系数 β MX ， β MX 和 C MX 是由实验数据[9]拟合出来的，本文拟合出来的数据
化工热力学题集
试题一一、 单项选择题（每题1分）1. 虚拟临界常数法是将混合物看成一个虚拟的纯物质，从而将纯物质对比态原理的计算方法用到混合物上。
.A ．正确B ．错误2. 对理想气体有A ．()0/<∂∂T P H B. ()0/>∂∂T P H C. ()0/=∂∂T P H D. ()0/=∂∂P T H3. 熵产生g S ∆是由于 而引起的。
A. 体系与环境之间的热量交换B. 体系与外界功的交换C. 体系内部的不可逆性D. 体系与外界的物质交换4. 以下4个偏导数中只有 是偏摩尔性质。
A ．()j n nv ns i n nU ,,⎪⎪⎭⎫ ⎝⎛∂∂ B. ()j n p ns i n nH ,,⎪⎪⎭⎫ ⎝⎛∂∂ C. ()j n P T i n nG ,,⎪⎪⎭⎫ ⎝⎛∂∂ D. ()j n T nv i n nA ,,⎪⎪⎭⎫ ⎝⎛∂∂ 5. 分离过程的有效能损失,x L E ∆ 。
A. > 0B. < 0C. = 0D.可正可负6. 超临界流体是下列 条件下存在的物质A. 高于T C 和高于P C B ．高于T C 和低于P CC ．低于T C 和高于P CD ．位于T C 和P C 点7. 纯流体在一定温度下，如压力低于该温度下的饱和蒸汽压，则此物质的状态为 。
A ．饱和蒸汽 B.饱和液体 C ．过冷液体 D.过热蒸汽8. 对理想溶液具有负偏差的体系中，各组分活度系数γi 。
A ． >1 B. = 0 C. = 1 D. < 19. 气体经过稳流绝热过程，对外作功，如忽略动能和位能变化，无摩擦损失，则此过程气体焓值 。
A. 增加 B ． 减少 C ．不变 D. 不能确定10. 一气体符合bV RT P -=的状态方程，在等温下当气体从V 1可逆膨胀至V 2时，则体系的熵变 ΔS 为 。
A 、b V b V RT --12lnB 、b V b V R --12lnC 、12ln V V R11. 蒸汽压缩制冷循环中，制冷剂的蒸发温度为200K,冷凝温度为250K ，则理论上该制冷装置的制冷系数ε为 。
利用Pitzer理论计算传统的凝固点降低法求活度系数的方法[发明专利]
![利用Pitzer理论计算传统的凝固点降低法求活度系数的方法[发明专利]](https://img.taocdn.com/s3/m/45feab579a6648d7c1c708a1284ac850ad020433.png)
(19)中华人民共和国国家知识产权局(12)发明专利申请(10)申请公布号 (43)申请公布日 (21)申请号 201611078457.X(22)申请日 2016.11.30(71)申请人 辽宁大学地址 110000 辽宁省沈阳市沈北新区道义南大街58号(72)发明人 佟静 李慧 郑煦 关伟 (74)专利代理机构 沈阳杰克知识产权代理有限公司 21207代理人 罗莹(51)Int.Cl.G01N 25/04(2006.01)(54)发明名称利用Pitzer理论计算传统的凝固点降低法求活度系数的方法(57)摘要本发明涉及利用Pitzer理论计算传统的凝固点降低法求活度系数的方法，包括以下步骤：采用自冷式凝固点测量仪测得溶液凝固点温度T f 与纯溶剂的凝固点温度T 0；通过公式θ＝T 0-T f得到凝固点下降θ值；通过公式计算得到的T f 时的a 1值；再就算得到T f 温度下的渗透系数ф，将渗透系数ф引入Pitzer方程，最终得出溶液的活度系数。
利用该方法测量的活度系数理论上清晰简洁，计算上简单高效，且测得的数据更加精准有效，应用的溶液浓度范围得以扩大，是一种新型的计算活度系数的有效方法。
是Pitzer与凝固点降低法的有效结合。
权利要求书1页 说明书4页 附图2页CN 106645267 A 2017.05.10C N 106645267A1.利用Pitzer理论计算传统的凝固点降低法求活度系数的方法，其特征在于，包括以下步骤：1)采用自冷式凝固点测量仪测得溶液凝固点温度T f 与纯溶剂的凝固点温度T 0；2)将步骤1)测得的溶液凝固点温度T f 与纯溶剂的凝固点温度T 0带入公式：θ＝T 0-T f 得到凝固点下降θ值；3)将步骤2)计算得到的θ值带入到公式(1)中，求得溶液凝固点T f 时的溶剂活度值a 1，-ln a 1＝1.612×10-6(6010θ+2.947θ2-0.01245θ3+…)＝9.688×10-3θ+4.75×10-6θ2-2×10-8θ3+ (1)4)将步骤3)计算得到的T f 时的a 1值带入公式(2)中求得T f 温度下的渗透系数ф，公式(2)如下：φ＝-1000ln a 1/2mM 1 (2)其中m为重量克分子浓度，M 1为溶质分子质量5)将步骤4)计算得到的渗透系数φ，引入Pitzer方程(3)其中，β(o)MX ,β(1)MX 表达离子M和X之间不同类型的短程相互作用以及从溶剂引起的间接力的作用；是三元相互作用的系数；参数b对所有电解质都是1.2；α对所有1:1价电解质是2.0；为渗透系数的Debye-H ückel参数，其计算过程如下：I为离子强度，计算过程如下：公式(4)中，N为阿伏伽德罗常数，ρw 为水的密度，ε为质子电荷，D为纯水的静介电常数；k 为常数，T是温度；公式(5)中，m i 为溶质i的重量克分子浓度，Z i 是i离子的价；6)将步骤5)经多次数据拟合求得的β(0)MX ，β(1)MX 的值，带入公式(6)，计算得到活度系数γ；权 利 要 求 书1/1页CN 106645267 A利用Pitzer理论计算传统的凝固点降低法求活度系数的方法技术领域[0001]本发明涉及一种利用Pitzer理论计算传统的凝固点降低法求活度系数的方法，属于物化计算领域。
第05章 Pitzer电解质溶液理论
2) 上式中第一项f(I)为长程静电能，它是离子强度I的函数； 第二项表示两个离子间的短程作用项（这包括离子间的硬
球排斥能和离子间色散能）。λij为两离子作用系数，它是
离子强度的函数；第三项为三离子间的作用项，μijk为三 离子作用系数。
3) 上式若保留第一项和第二项，则适用于低浓度下电解 散质溶液，加入第三项，适用于高浓度电解质溶液。
mi Zi2
1 2nw
ni Zi2
i
I
1
nw 2nw2
i
ni Zi2
I nw
所以 nw f (I ) f (I ) If '(I ) nw
School of Metallurgical Science and Engineering， Central South University
第二项微分结果为：
结果 与DH结果相比，拟合 结果有明显进步。 结论 基于长程作用和短程作 用共同考虑更能反映微 观粒子间的相互作用。
School of Metallurgical Science and Engineering， Central South University
Pitzer理论基点是统计力学，但严格地用统计力学 原理处理电解质溶液太复杂，因此Pitzer建立了一 个半经验式的统计力学电解质溶液理论。
对于长程相互作用贡献
d
1 1
6kT
i 1
i j
j 1
0
zi z je2
r
dr
gij (r)4r3dr
d
1
6kT
i 1
i j
j 1
0
zi z je2
r
dr
[1
铀成矿溶液中离子活度系数计算的几个问题
铀成矿溶液中离子活度系数计算的几个问题铀成矿溶液中离子活度系数的计算对于研究铀成矿机理和预测铀矿床富集的形成条件具有重要意义。
在铀成矿过程中，铀通常以离子形式存在于溶液中，其溶解度与离子活度紧密相关。
离子活度系数的准确计算有助于更好地理解铀成矿的物理化学过程。
离子活度系数的计算方法主要有Debye-Hückel近似法、Pitzer模型和位相关理论等。
下面将从这几个方面介绍铀成矿溶液中离子活度系数的计算的相关问题。
首先，Debye-Hückel近似法是最早被广泛应用的一种计算离子活度系数的方法。
该方法基于电解质在溶液中的离子作用力场，通过引入离子强度和离子的有效半径来计算离子活度系数。
Debye-Hückel近似法适用于低浓度溶液和相对较小的离子强度，而对于高浓度和高离子强度的溶液，该方法的效果并不好。
Pitzer模型是一种改进的离子活度系数计算方法，它考虑了离子之间的水合、共晶、共聚等效应。
Pitzer模型通过引入参数来修正溶液中离子活度的计算，从而更准确地预测高浓度和高离子强度溶液的离子活度系数。
Pitzer模型在铀成矿溶液中离子活度系数的计算中已经得到了广泛应用。
其次，铀成矿溶液中离子活度系数的计算还需要考虑体系的温度、压力和离子浓度等因素。
随着温度的升高，离子在溶液中的活度系数一般会降低，这是由于离子活动性的热效应所致。
此外，随着溶液中铀离子的浓度增加，离子间相互作用的影响也会增强，离子活度系数可能会下降。
最后，铀成矿溶液中离子活度系数的计算还需要考虑离子水合和配位效应的影响。
铀离子在溶液中通常以多种水合形式存在，不同水合形式的离子活度系数不同。
此外，铀离子还与其他离子形成配位络合物，这也会对离子活度系数产生影响。
因此，在计算铀成矿溶液中离子活度系数时，需要考虑到离子水合和配位效应的综合影响。
总之，铀成矿溶液中离子活度系数的准确计算对于研究铀成矿机理和预测铀矿床富集的形成条件具有重要意义。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
Pitzer活度系数模型研究与开发作者：韩莎莎郑俊强孙晓岩项曙光来源：《当代化工》2020年第01期Research and Development of Pitzer Activity Coefficient ModelHAN;Sha-sha，ZHENG;Jun-qiang，SUN;Xiao-yan，XIANG;Shu-guang（Process Systems Engineering Institute， Qingdao University of Science and Technology，Shandong Qingdao 266042， China）自然界、生命体和工业过程中普遍存在着电解质溶液，是化工行业中的重要组成部分，也是众多过程处理的对象，目前逐渐成为许多有机物和無机物反应的良好媒介，因此对电解质溶液的理论研究、电解质溶液的热力学性质的研究及电解质过程模拟研究具有重要的工业实用价值和理论意义。
其中在电解质溶液理论及含电解质溶液体系的热力学性质方面，Debye[1]、Meissner[2]、Bromley[3]、Chen[4]、陆小华[5]、左有祥[6，7]、Loehe[8]、李以圭[9]和杜艳萍[10]等都做出了很大的贡献。
目前Pitzer是用于计算水电解质溶液体系（尤其是离子强度为6摩尔以下的强电解质体系）的活度系数等热力学性质较为准确的电解质活度系数模型，也是应用最为广泛的电解质溶液理论。
最初1973年，Pitzer修正了D-H理论[1]，得到了经典的半经验Pitzer模型[11]，但适用的浓度较低。
随后为了扩大浓度适用范围，用Margules方程修正了短程项，得到了Pitzer[12]（1980年）模型。
之后，Bromley[3]（1973年）简化的Pitzer模型、Pitzer[13]（1975年）添加的静电非对称混合项、Fürst和Renon[14]（1982年）研究的多种参数对模型用于1-1型电解质固液平衡的影响、李以圭[15，16]（1986年）的Pitzer-Li方程、Simonson等[17]（1986年）的Pitzer-Simonson方程、Kim等[18，19]（1988年）回归的高浓度体系参数、Clegg等[20，21]（1992年）的Clegg-Pitzer模型、李以圭等[22，23]（1994和1997年）的Li-Mather模型、Pitzer[24]（1999年）以及Chen等[25]（2008年）都对Pitzer模型做了相应的修正和完善。
因此，参照Fortran语言编程如Zemaitis[26]中实现含电解质体系的模拟计算过程，也可通过Visual C++编程语言开发Pitzer模型，实现被已有的支持CAPE-OPEN标准的大型通用化工模拟软件所调用，从而对工业中含电解质溶液过程进行设计、模拟、计算和优化，更好地解决较复杂的工程问题。
本文主要是根据Pitzer修正的水电解质溶液体系活度系数计算模型[13]（1975年模型）进行开发并通过对一些应用实例的模拟计算并验证结果对该开发的Pitzer活度系数模型进行分析、讨论和评价。
1 ;Pitzer活度系数模型下文将对Pitzer活度系数模型作详细阐述。
Pitzer活度系数模型以水电解质活度系数模型为基础，区别于其他物性方法的是它与其他活度系数模型无重叠。
它能够准确地计算高达6 mol离子强度以上的电解质水溶液的行为，但不能用于存在除水之外的任何其他溶剂或混合溶剂的体系，任何非水的分子组分都应被认为是溶质，分子溶质存在与否也都可满足，并全被看作亨利组分。
超临界气体的溶解度是采用亨利定律来模拟计算的;汽相逸度系数采用Redlich-Kwong-Soave状态方程，并且所有其他汽相性质也都假定是理想的，但遇到气相中存在缔合行为时，例如羧酸类（RCOOH）或氢氟酸（HF）体系，Redlich-Kwong-Soave则不适用。
对于羧酸类，选择一个非电解质活度系数模型与Hayden-O‘Connell或Nothnagel方法相结合;对于HF体系则选择专用的ENRTL-HF或WILS-HF 物性方法。
本文选取的Pitzer活度系数模型即水电解质活度系数模型计算的系统总过剩吉布斯能的基本方程见公式（1）。
式中：G E;—剩余吉布斯自由能;B、C、θ、ψ—模型中交互参数;;;;;;;R —气体常数;;;;;;;m —质量摩尔浓度，mol/kg;;;;;;;Z —电荷数。
其中长程静电相互作用项为：即可推出活度系数模型：其中： ;n;—摩尔数;*;—非对称。
2 ;模型求解步驟在进料物流的温度T，压力P和质量摩尔浓度m i或组成z i已知条件下，其各个种类（包括所有分子和离子组分）的活度系数的详细求解步骤如下：Pitzer活度系数模型用于计算电解质体系所需要用到的参数主要包括：阳离子-阴离子间的参数β（0）、β（1）、β（2）、β（3）和C;φ;阳离子-阳离子间的参数θcc';阴离子-阴离子间的参数θaa';阳离子1-阳离子2-共同的阴离子三个离子间的参数ψcc'a;阴离子1-阴离子2-共同的阳离子三个离子间的参数ψcaa'以及分子-分子间、分子-离子间的参数β（0）、β（1）和C;φ;还有静电非对称混合效应项Фij。
以上所有参数都可参考Pitzer相关系列文献[26-30]，从而建立一个小型数据库，供开发的模型模拟计算使用。
（1）计算阳-阴离子间的交互参数B ij和C ij因对于不同价型的电解质类型，B ij的计算公式也不相同，其中类似于1-n型（1-1，1-2和2-1等）电解质，B ij和C ij都是针对带不同电荷的离子即阴阳离子间的，其计算公式如下：针对n-m型的电解质并且n大于1、m也大于1（如2-2，2-3和3-4型等），B ij的计算公式如下：其中α1取值2.0;α2取值12.0;α3取值1.4。
相对应的求导得到B ij'（对于1-n型电解质）：求导得到B ij'（对于2-2型电解质）：（2）计算带同种电荷的离子间交互参数θijλ为短程离子间效应描述为离子强度函数。
（3）计算三元交互作用参数ψijkΛ是忽略任何离子强度依赖的三离子交互项。
其中（2）和（3）计算的参数是多电解质组分溶液中可能存在的双离子和三离子间的相互作用。
对于单一的电解质溶液，它们的值都为零。
（4）计算静电非对称混合效应项Фij求导得：对于参数X的计算，根据体系类型不同，其计算方法也不同，当不存在非对称的混合校正项时：当静电非对称混合修正项用多项式表达时为：当静电非对称混合修正项用积分形式表达时为：长程作用项：其中Debye-Hückel参数Aφ、离子强度I x为：水的密度為：以上式子中：k —玻尔兹曼常数;N A—阿伏伽德罗常数;Q e—电子电荷;Ɛ;—介电常数;M;—相对分子质量，kg/kmol。
对f（I）进行求导得：以上各类参数将其代入（3）式就可计算出体系中各个种类的活度系数。
离子种类的活度系数模型：对于水，类似地计算活度系数的对数如下：由于计算出的是以质量摩尔浓度m i为单位标度的活度系数值，可由下列公式将其转化为以摩尔分数x i为单位标度：离子类：水分子：在上一节中求出lnγi之后，即可求：由于计算H m*E需要对活度系数表达式求导，为了避免复杂的求导过程，可以采用以下方式计算：很明显看出，这里采用的是对G m*E进行差分法求导计算：其中T2-T1可以取一个很小的值0.001。
H m*E和G m*E计算完成之后，再按照常规体系的总焓和总吉布斯自由能计算模型得出总焓值和总吉布斯自由能值，既而得到总熵值。
3 ;模型开发通过Visual studio 2005平台采用C++语言进行编程实现，具体开发过程如下。
（1）定義变量//定义组分列表vector<int>compIdList;vector<BOOL>henryList;//定义组分组成vector<double>xm;//定义参数数组集vector<vector<double>>B0ij，B1ij，B2ij， B3ij，Cij，θij，ψijk;（2）函数类型//创建计算活动系数所需的相互作用参数数组。
voidInitialParameters（）;//计算长程静电相互作用项CalcLongRange（double T，vector<double>&x，vector<double>& LR）;//计算短程位能项CalcShortRange（double T，vector<double>&x，vector<double>& SR）;//计算三离子交互项CalcThreeIonInteraction （double T，vector<double>&x，vector<double>& TII）; //计算各个种类的活度系数GetActivityCoefficient（double T，vector<double>&x，vector<double>& ac）;GetActivityCoefficient（CVariant&compIdList， CVariant& composition， double P， double T， VARIANT* value）;//计算Hm*EGetExtraEnthalpyDepature（CVariant&compIdList， CVariant& composition， double P，double T， double * extraEnthalpy）;（3）对亨利组分的处理//初始化一个数组，对应于所有的组分，如果组分为亨利组分，则对应的位置为，否则为0henryList=vector<BOOL>（count，0）;for（inti=0;i<henryComs.GetCount（）;i++） for（int j=0;j<count;j++）if（henryComps.GetLongAt（i）==compIdList[j]）{henryList[j]=1; break;}//以上是亨利组分的识别，识别出的亨利组分和可凝组分的交互作用参数都设为0。
（4）计算纯水的性质CComPtr<IUtilitiesManager>utilitiesManager;//公用工程的接口，计算纯水摩尔体积和焓值utilitiesManager->GetMoleWaterVolumeByTP（T，P，&moleWaterVolume）;HRESULT hr=utilitiesManager->GetMoleEnthalpyByTP（T，P，value）;4 ;模型验证与对比以实际工业生产中的含电解质流程实例对开发的Pitzer模型进行模拟验证，体系中涉及的汽相的性质采用RK-Soave状态方程计算，液相的热力学性质则采用本文开发的Pitzer模型进行模拟。