2023年8月3日发(作者:)

中南天学硕十学位论文摘要摘要旧版《氧化铝生产热力学数据库》在过去的使用中表现出了强大的生龛力,同时也暴露出了诸多问题。论文针对此,进行了全面的软件优化和版本升级工作,并提出了全新的铝酸钠溶液活度系数估算膜型,实现了数据库对NaAI{OHh-NaOH—H20体系活度系数及其它相关热力学性质的计算。同时,依靠数据库,论文对氧化铝溶出过程中的主要体系进行了热力学分析。新版本数据库按照数据库软件开发的要求,在用户界面的精简性、热力学数据的完整性、计算结果的准确性、反应方程式的配平等方面进行了完善,并对相应的程序代码进行了详细的修改及性能优化,彻底消除了旧版数据库使用中诸如使用不便、配平程序失效、计算失败等同题。此外,为方便用户对数据库的操作及优化升级,特制作了数据库帮助系统,并编写了《氧化铝生产热力学数据库用户手册》和《氧化铝生产热力学数据库程序说明文档》。着眼于氧化铝生产过程中最为重要的NaAI(OH)4一NaOH.H20体系,建立了基于Debye.H0ckel理论的铝酸钠溶液体系活度系数估算馍型,包括:NaAI(OH)。活度系数计算模型以及铝酸钠溶液中Na2COs、Na!S04和Na2H2si04活度系数计算摸型。NaAI(OH)4活度系数计算模型在t=25100。C,aK=1.5~10.5,mN抑HK=O.1~lOmol/kg范围内适用,成功解决平衡溶解度数据在苛性比和苛碱浓度范围上的限制。同时,还对25"C~100"C范囝内Na20.A1203.H20系中三水铝石平衡溶解度估算模型进行了优化,拓宽了其适用范围,提高了其准确度。为检验数据库的准确性及应用价值,通过数据库,论文对拜耳法生产氧l中南人’半’硕十学侍论文摘要化铝溶出过程中的Na20一A1203一H20、Na20一A|203一Si02一H20、CaO.A1203.Si02.H20等体系的主要反应进行了热力学计算和分析,得出了以下结论:(1)高岭石系列物质溶出反应的优先序列为:多水高岭土>迪开石>高岭石:(2)提高温度、Si02和Na20k浓度,均有利于钠硅渣的生成;(3)低饱和系数的水化石馏石容易生成:(4)提高SiO:浓度f『利于水化石榴石的生成;提高温度、Na20k浓度和CtK愈高,均不利丁J水化石榴石的生成;(5)在稳定性方面,低温下水化石榴石稳定,高温下钠硅渣稳定。结果证明,数据库计算结果准确,热力学趋势与工艺中大部分实验事实相吻合。这从一个方面证实了数据库的可靠性和实用性。此次开发优化工作以计算机技术和氧化铝专业知识理论为保障,全面实现了数据库的易用性、准确性、适用性和安全性,确保了新版本数据库合理的生命剧期以及良好的应用前景。关键词氧化铝,热力学数据库,优化,应用n竺苎望!!!竺!塑!!竺呈!!生兰望垒里苎!坠垒篁!TheOptimizationandApplicationofThermodynamicDatabaseforAluminaProductionABSTRACTDuringtheuseofThermod)’namicDatabaseforAluminaProductionoldversion,manyproblemswerefoundinspiteofitsgoodperformance.Inthispaperthesoftwareoptimizationandeditionupgradewerecarriedouttosolvethoseproblems.AndthenewestimationmodelforactivitycoefficientsinNaAI(OH)4-NaOH-HzOsoftwaretOsystemwasestablished,whichwasusedinthedatabasedealwiththerelevantthermodynamiccalculationsimmediately.tOAccordingthedemandofthesoftwareexploration,thecomplementarityuserofthesoftwarewasmadefromthesimplenessofinterface,theintegrityofthermodynamicdata.theaccuracyofcalculationaIresultsandtheequilibriumofreactions.AndalIoldrelevantcodeswerecompletelyamendedtooptimizethewholedatabasecapability.Maindisadvantagesoftheolddatabaseversionwereasovercomedduringthefurtherexploration,suchtheinconvenienceofusing,theinvalidationofequilibriumandthefailureofcalculation,Inaddition,thedatabasehelpsystemwasestablishedand“TheUserManualofThermodynamicDatabaseforAluminaProduction”and“TheProgramInstructionsofThermodynamicDatabaseforAluminaProduction”werewrittenforsupplyingconveniencefortheuserstOoperateandupgradethedarabase.aFocusingontheNaAI(OH)4-NaoH—HzOsystemwhichplaysparamountroleintheindustryofaluminaproduction,acompleteestimationmodelforactivitycoefficientsinNaAI(OH)4一NaOH-H20systembasedontheDebye—Hiickeltheorywasproposedinthispaper.Itincludes:theestimationmodelsforactivtitycoefhcientsofNaAI(OH)4、Na2C03、Na2S04andNa2H2Si04inNaAI(OH)4-NaOH-H20system.Themodelcanbeweltappliedintherangeof25"C~100"(2,q=I.5~10.5,mNnOHK=0.1~10mol/kgtherangelimitsofthevaluesofO.KandsuccessfullysolvedandmNaoHgcalculatedfromsolubilitydata.111璺竺!竺堡!!垒墅!!坠£!Simultaneously,theestimationmodelsforequilibriumsolubilityofaluminaintheNa20.A1203一H20systemwenoptimized.Anditsapplicationandaccuracywerewidenedandimproved.InordertOcheckaccuracyandtheapplicationofthedatabase,themainreactionsintheBayerprocessofaluminaproductionsuchasNa20一A1203一H20、Na20-A1203-Si02一H20andCaO-AIz03-Si02一H20systemwerecomputedandanalysedinthispaper.Andsomeconclusionsweregainedasfollowing:(1)Thepreferentialsequenceofthereactionwithalkaliis:halloysite>dickite>kaolinte;(2)IncreasingthetemperatureandtheconcentrationofSi02andNa20kwillacceleratetheformationofthesodiumhydratehydratewithlowersaturationalumina—silicate;(3)Thegarneteasier;(4)Highconcentration0tKcoefficientwillformofSiozandlowtemperatureandlOWconcentrationofNa20kandfavorablefortheformationofthegarnethydrate;(5)AsgarnethydrateisissteadiersteadieratattOwillbethestability,thesodiumhydratelowtemperatureandthealumina-silicatehightemperature,Theresultsshowthattheareaccuratecalculationresultsfromthedatabaseareandthethermodynamicgendsinaccordancewiththeexperimentalfacts.ThishighlyprovedthereliabilityandpracticabilityoftheThermodynamicDatabaseforAluminaProduction.Thefacility,veracity,applicabilityandsecurityofthedatabasesoftwareweregreatlydevelopedbasedonthecomputetechnologyandaluminatheory.ensuresThenewversionofthesoftware6lmre.thesarisfy.inglifecycleandapplicableKEYWORDS:Alumina,thermodynamicdatabase,optimization,application原创性声明本人声明,所呈交的学位论文是本人在导师指导下进行的研究工作及取得的研究成果。尽我所知,除了论文中特别加以标注和致谢的地方外,论文中不包含其他人已经发表或撰写过的研究成果,也不包含为获得中南大学或其他单位的学位或证书而使用过的材料。与我共同工作的同志对本研究所作的员献均己在论文中作了明确说明。作者签名:Z同期:年月同关于学位论文使用授权说明。本人了解中南大学有关保留、使用学位论文的规定,即:学校有权保留学位论文,允许学位论文被查问和借阅,学校可以公白学位论文的全部或部分内容,可以采用复印、缩印或其它手段保存学位论文:学校可根掘国家或湖南省有关部门规定送交学位论文。艏戤:删刷性g.嗍:年月同±塑叁主竺±芏竺鲨塞一苎二垦奎苎笪垄第一章文献综述1.1热力学概述㈦热力学是关于热、助、能以及在系统状念中由它们产生的变化的研究.是一门研究能量转换以及与转换有关的物质性质(或参数)问相互关系的科学,从广义上讲,热力学是研究系统宏观性喷之『日J的关系。热力学中一个关键的性质是温度,有时热力学也破定义为温度与物质宏观往质之间的关系的bK究。热力学足稳学的宏观理沦.是坂能量转化的观点研宄物质的热胜质。阐明能量从一种形式转换为另一种形式时应遵循的宏观规律,热力学是根据实验结果综合螯理而成的系统理论。它不涉及扬磺的微疆结构和微观粒子的相互作用,也不涉及特殊物质的具体性质,是一种唯象的宏观理沦,具有高度的可嘉眭和普遍睦。热力学的完整理沦体系是由几个基本定律以及相I立的基本状念函数构成的,这些基本定律是以大量实验事实为根据建立起束的。无论多少个物体互帽接牲都能达到热平衡.并且如果A物体同时与B、C两物体处于平衡念,则B、C两物体接触时也一定处于平爱念丽不发生新的变化,这一热平衡规律弥为热力学第零定律。出此可以引入一个状态函故,品宦。温熏是判定一系统是否与其他系统互为热平衡的杯志。热力学第一定律就是能量寸恒定律。是后善在一切:步及热观象的宏哩过程中妁具体友观。描述襄统热运动能量的状念函数是内能。通过作J力、传手丸.系统与外界交换能奄,内能改变。热力学第二定律指出一切涉及热髓象的宏观过程是不可逆的。它阐明了在这些过程中能量转换或传递的方向、条件和限度。相应的念函数是熵.熵的变化指明了热力学过程进行的方向,熵的大小反映了系统所处状念的稳定性。热力学第三定律通常表述为绝对零鏖时,所有纯物质的完美晶体的熵值为零。它意味着。当温度趋近于OK时,所有;疑聚物系的墒值趋近于一个相同的极限值。1.2热力学计算模型物质基本热力学性质的基准叁约定如下【¨51:(1)焓:对纯物质束I兑,规定标准状态下(298K,l大气压下)稳定单质的生成焓H。298为零;对离子来说,规定H+的生成焓拶298.H-为零。(2)熵:在绝对零度时,完整晶体的纯物质,其熵值规定为零:对离子来说,规定任中南人学硕十学付论文第一章文献综述何温度下H+的标准熵为零。1.2l纯物质热力学性质的计算模型【5·6‘7】绝物质热力学性质的计算主要包括:计算任一温嚏T(K)下的热容(Cp)、焓(H-r)、熵(ST)和吉南斯自出能(GT)。计算表达式分别为:Cp=a+b×10‘3T+c×106T-2+d×lO’叶2或Cp=a+bxl013T+c×lO‘T2其中a、b,C,d为热容系故H-r=H。29s+『c,。dT+JCp.dT+.’+』c,dr+∑叫。(1·2)ST=S。298+狰刀+净T仉…+辟JT仉莩半GT=HT·TxST(1-3)(1—4)式中的下标0,l,2…。rl表示相念,AHT.表示相变热。1.2.2水溶液中离子热力学性质的计算模型f%1,1)计尊水溶掖中离子在任一温度T(K)下的焓(H。)、墒(Si)和吉南斯自由能(G.)。计算表达式分别勾:r△H.=△H。299十解刀29fl,‘(1—5)/kS,=ZXS。29s+:。t争△G.=△H。-Tx△S。小,可以忽略不计。故其均未考虑压力的影响。(1·6)(1·7)注:以上三式同样适合纯物质相关热力学性质的计算。在300"C以下压力的影响很(!)计算水溶液中离子的热容(C。)值由于离子的热容值还很缺乏,一股采用Criss提出的“离子熵相对原理”‘42“■束进行近似计算。其汁算方法如下:Cp。=BxTxl04(1.8)在这里,Cpo为离子在温度298和T之闻的平均热容,B是热銮系数。从式(1.8)可以看出,离子的平均热容与温度成线性关系。B54.184×(a+bxS*29&4bs】)(1-9)中奄人学母!十中伊论文第一摩文献综述式中,a.b是仅同离子类型有关的常数。表l一1列出了不同炎型离子的a.b值。S。29sIab。是离子在298K下的标准绝对墒。通过(4—10)式可以计算得出。S。298lahs}2S。208—2092、Z(1.10)其中,Z为离子电荷数(带符号),·20,92(J—mol-*.K。)是H+在298K下的绝对熵。S。!。B为离子的相对墒。可从数据库中查到。表I-t各美耋}子的a,b值嘲Table1-lvaluesofaandbofvariousiOnS离子类掣打~叫a竹…由,简单田f离f简单胡离f及OH‘禽氯珞胡离f’’‘…“~禽氢氧珞月子………。。i页再F—1砭孺r——_j丽矿————i而=广一(XO”。)(XO。(OH)”1)1.65x10。。,I.12x10。5.83x10。111.2x10‘4bI3)H‘热力学性质的计算H+不楫f四种类型的离f之一.其热容与,甚度之问的关系为COpH‘=0.3912T(I一11)从而可以得到如下关系式:H。TH+=0.5x0.3912。(T2-298,152)(I—12)f1.131(1.14)s。rH+(absl=S。298HjabsI+o.3912(%298.15)G。TH+=H。TH+一TxS。r14i曲s1.2.3水溶液中分子热力学性质的计算模型一股情;兄F,弱酸弱碱在水溶液中不能完全离解,呈中性分子的状叁。水溶液中分子的媳力学计算通常采用Helgesonl81给出的近似公式.浚公式主要适用于段低的温度(200℃以下)。假设水溶液中一未离解的中性分子A。B。.其离解方程式为AmBn专mAY++nBz.(1-15)式中m与11一电离时产生的阳离子数和阴离子数Y+与z-——阳离子电荷数和阴离子电荷数离解反应(1一15)的标准反应自由能为:4△G。T(d。s}=一△S629815(d6,[298.15--二fl--exp(exp(b—aT)--c+(T-298.15)/e)}】+△H629815(dls}珊(1·16)式中△S。298.15(dis)一离解反应(1一15)在298.15K时的墒变生童叁竺堡±兰生堡皇AH。298l蛐。)——离解反应(1一15)在298.15K时的焓变.坚一奎苎笙堕a=0.0187:b=.12.74l:c=exp(b+298.15a):0=219:∞5l+ac0离解反应(1—15)的墒变为:△S。T㈨J-一AS。29815m1{exp[exp(b+aT)一e+(T-298.15)/0]}[1+qoexp(b+aT)I/∞其中0=4.106(1·17)离解反应(I一15)的焓变为:AHoT(dⅢ=AO。T【dⅢ一TAS。T(dIs)(1一18)另一方面。离解反应(卜15)在任一温度下的反应吉钿斯自由能也可以写成:△G。Tldtsl=从而mGoTIA,+nO。TtBl—G。T【AmBnl(1·19)(1-20)(1-21)G。TiAmBnl=mG。TIA)+nG。nBl—AG。Tlald¨同理:HoT【AmBn)=mH。T‘A1+nHoT(Bl-△H。T(dIs)在式(1.20)(1.21)中A,B均为离子,其高温下的热力学性质可通过离子相对熵原理求出。1.2.4化学反应的热力学计算模型【35·79I对化学反应来晚。任一温度T(K)下反应焓(△HT)、反应熵(△sT)、标准反应吉白斯自由能(AGr}以及平衡常数K的计算可通过下列公式进行:AH。r=∑bⅣI)一∑0ⅣIJ’I啦拐(I一22)址’硼△s。t=∑b,s-)’}墟锡一∑h黑j(1-23)衄J4罚AO。t=∑b60)一∑bG0)(1·24)厘pt研1成研△GoT=.RTLnK即LnK=.AG01嬲(1-25)式中,H。i’T、s。i,T、G。i,T分别表示i物质的焓、熵及标准吉南斯自由能值;vi表示化学反应计量数:K表示平衡常数。必须指出的是,现有的热力学数据手册,有的列出的是物质的生成自由能△G。f=△H。r—T△S。f有的列出的是章尔吉如斯自由能GoT=H。T--TS。T无论使用哪一种,只要不交叉使用,计算出来的同一反应的吉白斯自由能是一样的·在4中南大中硕十中何论文第一章文献综述本论文中,使用的是生成吉如斯自由能和尘成焓。在以上的讨论中,△GoT代表的是化学反应的标准反应自由能。杯准的含义指的是每一种参加反应的气相物质分_lI部处在】大气压下,参加反应的溶r瘦浓度为l(mol·L“)。而实际7芷产中的情况往往比标准状念条件复杂得多,所以,需要讨论在任意指定的状念下,体系能否发生化学反应。这就要用到下面的等;盛方.智式(1--26)束作为判断化学反应方向的补充。假定反应为aA+bB=IL+mM△GT=△G。T+RTInJ=一RTInK+RTlnJ(1-26)舯,j-籍或籍虹绷混色PL,PⅥ,PA,PB,aL,a、I.aA,aa分别为反应物和生成物的指定分压和反应活变。从(I.26)可以看出,当J<P时.△GT<0,反应能够自发地向右进行,反之.则不能。严格.兑,判断反宦方向应以(1.26)为依掘,但因为J在对数二女内,所以△GoT对△GT的影叫更大。也就是:蔸.当△GoT的绝对值很大时.△GT的IF负号由△GoT束决定。一殴地l并.△(jor<一4I.84f焦时.反应不仅能在标准‘帙仑F进行.在一收惰:兄下也能进,-:△G。T>41.84干焦时.反应则不能自发进行㈣。1.2.5反应的AGoT—T二项式的回归模型在实际计算中,AG。T与T的多项式用起柬比较寐颐,人们喜欢使用G。T的二项式:△G。r=A+BxT。通常可将化学反应的△G。T=f(T)多项式整理成二项式或根据一系列比较准确的/XG。T值利用回归分折法(即最小二乘去)求得AGoT的二项式。按指定温度计算的值越多,用回归分折法计算出柬的AG。T二项式就越准确。最小二乘法的基本原理如下:实验值G,、T,(i-l,2t…,n)误差6i=GI—A—B×Ti目杯函数J_62.=(G.--A--BxT.)2根据极值原理,分别对A、B求偏微商并令其为0,从而得出A、B值:A=否一B×于(1·28)(1·29)B=∑(7=一于)(q一否)/Err,一_2中南人中硕+一中何论艾第一章文献综述瓦甲于2吉喜7:否2吉萋GI相关系数,:圣!三二!!呈三童!1.2.6复杂化合物热力学数据的估算方法1.1.6.1组成相似的无饥化合物热力学陛质的线性估算模型以水合铝酸钙为例加以说明:水合铝酸钙主要存在的分子式形式有:3CaO·A1203·6H20,4CaO·A1203·13H20,4CaO·A1203-19H20.CaO’A1203‘10H20以及2CaO·A1203"8H20等。从氧化物裹示的组成可以看出CaO和H20的分子比在变化,A1203的分子比始终为l,可设汁其,上成焓线性方程式如下:△Ho=△Ho。+m(CaO)△HI。十n(H20)△H二。(1-30)其中△Ho为各水合铝酸钙的生成焓:△Ho。、△HI3、△H2。是和生成焓有关的参数:111、n分别是复杂化合物以氧化物襄示的CaO和H:O的分子比。在标准状念下,各种矸三念的水合铝酸钙的生成焓均可以在数据手册中查副,以多元线眭最小二乘法进行拟合,得到如下结果:△Ho(1d·mol“)=一1633.23.704.2m(CaO)一29498n(H20)同理,可以对水合铝酸的自由能、熵和热容进行拟合,得到的结果如下:AGor(kJ·molll=.1574.76-672.88m(CaO)-237.08n(H20)So(1d·mol“·K“)=57.0l+39.14m(CaO)+38.92n(H20)组成相似的化合物在相同温度下的恒压热容之问也存在一定的线眭关系。即:Cp(T)=Cpn+m(CaO)Cpl+n(H20)Cp2(1-31)通过上式。可以求出复杂无机化合物在不同温叟下的热奋值,然后再利用线性回归求出热容系数,从而可求出其它温度下的热力学性质。1.1.6.2复杂无机化合物(至少包含两种令属离子)的拆分估算模型根据温元凯算法的主体思想,把任何复杂硅酸盐矿物看作各种酸性氧化物和碱性氧化物的复合物。复合时碱性最强的氧化物首先和酸性最强的氧化物反应,然后再和酸性次强的氧化物化合【lo¨l。水在大部分情况下以水合£F的形式存在,有时也作为结构水参与反应。这样,复杂硅酸盐的吉钿斯自由能(或焓)可看成由各氧化物的吉白斯自由能(或焓)之和(EGoOxideorEHoOxide)与氧化物间的反应自由能(或反应焓)之和(EGoRorZHoR)6中南人学硕十学伊论丈第一章文献综述两部分组成。1.3建立氧化铝生产热力学数据库的目的和意义目前,我固的氧化铝尘产还没有一套完整的仿真与自控体系。计算机技术在氧化铝生产中的应用同国际上相比,还存在着很大的差距,特别是蘑f汁算机模拟与控制的软件和数掘库1l常缺乏.已成为发喂爱化铝生产自动化的}要贿碍。今天,任何计算机过程摸拟与自动控制的应啊软件均是建立在大型的热力学数据库的基础之上的。因此,要提高我国计算机在氧化铝生产中的应用水平.加强工业生产的自动化。除了吸收引用国外的先进技术柬结合我固的生产实际外,加强我囤自己的软件和数扼库丌发能力,大力丌展数据库、控制软件等整础性的研究工作十分迫切。我【川的氯化铝生产具有很强的“中因特色“,这是由我固的铝土矿资源决定的,其特点是:(1)铝上矿资源#富,遍及19个省,居世界第四位;<2)铝t矿矿石类型以一水硬铝石为主,占全国总诸量的99%,三水铝石型铝上矿仅约占l%;(3)绝大部分矿石品位冒中等,铝硅比多分1}】在4 ̄10之1日J.约占全国总锗星的85.6I%:(4)矿石具有高铝、高硅、低铁、中等};吕硅比,难藩的特点,且伴生矿多.成份复杂。铝土矿矿右_的这些{寺点决定了我雷的氯化锅r业主要采用混联法和烧结法.目前.除了乎果铝厂采用纯玎尔法外,其它五个氧化铝厂均采喇烧结法或联合法。同玎尔法相比,烧结法工艺流程长、连续性差、能耗高、产品质量差,成本高。我国氧化铝工业的发甏严重滞岳于电解铝工业,每年复化铝缺口部比较大【‘2-’314”】。如何提高产品质量,进一步缩小我国与世界发达国家在氧化铝产能、成本等方面的差距,提高我国氧化铝生产在国际上的竞争力,己成为我国蓑化铝生产厂家主要的攻关方向,因而必须采用新的生产技术去改善(或完善)我固的氧化铝生产工艺,这就需要对氧化铝,土产过程的基础理论有深一步的认识。备矿物及其在体系间的相互反应热力学是氧化铝生产基础理论研究的一个非常重要的方面,其可为新工艺的提出、实际生产控制等提供理论支持。然而,氧化铝生产中的热力学研究目前存在蓄许多困难:l、出于铝土矿及其它含铝的矿物资源如霞石矿、明矾石矿以及各种铝硅酸特岩石等均是组成复杂、化学成分变化很大的矿石,因而在氧化铝,上产中存在着物相多、物理化学行为极其复杂等特点。特别是所涉及的物质大多为二元、三元或三元以上的复氧化物,这一类物质可供参考的热力学数据菲常缺乏;2、现有的热力学数据分散、零乱,没有公认、统~的数据标准:3、计算复杂,手工计算难于完成,且准确性较低:4、缺少先进、现代化的工具进行辅助分析。7中奄人中硕卜中伊论文第一章文献综述针对以上问题,中南大学氧化铝及陶瓷研究所将热力学、计算化学以及数锯库技术的结台起束,研究和丌发氧化铝生产热力学数据库。它的建设成功将解决我国氧化铝工业界广大科技工作苦长期面}临的某些工作难题,可对氧化铝生产辐关工艺过程进行模拟i同时.目前在我国还很少有关于氧化铝生产热力学数掘库方面的报道,它的建成也将填补我国在浚方面的一项空白。I.4本论文的主要研究内容I、氧化铝生产热力学数据库的优化和完善.包括耐用户界面的精简性、热力学数掘的完整性、计算结果的准确性,反瘦方程式的配平等方面进行完善,并对相『萱的程卑代码逆行详细的修改及性能优化:2、编写数堀库源代码程序说明文档以便于今后优化和升级工作的进行,编制数掘库用白使田手册和帮助系统以方便用户使用;3、铝酸钠溶液中活度系数计算模型的完善,含铝酸钠溶液体系中NaAI(OH)。Na_-C03、Na2S04和Na2H2Si04活壹系数计算馍型的建立以及三水铝石平衡溶解度估算漠型的陶造:4、氧化铝尘产热力学数据的应用,包括对玎耳法,圭产氧化铝溶出过程中的Na:O—A1203一H20、Na20·A1203·Si02一HrO、CaO-A1203一Si02一H20等体系中的主要反应进行热力学计算和分析。中啕人学面十学伊论文第.章氧化锅J上产执力学数据的基础俨究第二章氧化铝生产热力学数据的基础研究2.1与热力学相关的基本概念1.状态函数mI(StateFunction)热力学中热力系统在某一瞬问所处的某种宏观物理状;兄称为状念。用宋规定状念的性喷叫做热力学变星或状念性质.一旦规定了状念的{丸力学,殳璧.其他平衡状念性质就是这些变蕈的函数,称为{丸力学函数或状念函数。如压力、温室和比容等。状念参数的数煎仅仅取决于热力系统的状态,丽与达到这种状念昕经历的热力过程无关。2.标准状态㈤(StandardStateS)气体:在p_101.32kPa下的理想气体液体:在p=101.32kPa下的纯液体固体:在p=101.32kPa下固体所处的最稳定的晶体状态3.热容…920I(HeatCapacity,)一定量的物喷,在不发生化学反应和相变的情;兄下。温度升高(或降低)l℃(1K)所r受收(或放出)的热量称为该物质的热容。当系统吸收或放出热量时,它的i昱度可能变化,也可能不变化,这取决于所经历的过程。在传递热量Q的过程中,系统的温室从TI变化到T二.则平均热客定义为:平均媳奋=—羔一7:一Il当Q和T2一TI部趋于零时,这个比值即为热容C:c:1im旦:丝,j-^7j一7;dT:!生z-Ij?j一7;.dT堡c:Iiml每单位质量的热容称为比热容:单位为J,(kg·K)。每单位气体体积的热容为容积热容:c,:.!一塑:风塑22.414珂d7’”d7’单位为kJ/(n·m3.K)。,00为标准状态下气体的密度如果物质的量是l摩尔,它的热容则称为摩尔热容:G=吉等=M拶dq单位为J,(m01.K)。9!塑叁之堡±竺竺堡奎墨三兰墨竺堡竺兰!堕!墼堑!!蔓壁竺茎等匝过程中的热容量成为定压热容:Cp=(等),等客过程中的热容量成为定容热容;c,=【掣),口1热容与温度的关系式一般有两种方式:Cp=a一-b×10‘3T+cxl0吁2+dxl04T2和Cp=a+b,10.jT+cx!05r2。热吝增晕-ACp:化学反应中.产物的总热容减去反壹物的总热容。4.焓‘212二J(EnthalpY)对任何系统来说,焓的定义为:At=U+Pv式中u为系统的内能,尸为系统的压力,V为系统的体积。焓是状念函数。当系统进行某一过程时。焓的变量£XH只与始叁、终奁有关,而与中jt日J过程无关。单位为J/mol或kJ7tool。杯准焓日o:25*6"时处在杯准叁的l摩尔物质的焓值。标准生成焓△H:,:当觇定每种电质在lol325Pa和25。C时的焓值为零,则由有关草质左101325Pa和25"C时生成lmol化合物所吸收的热量为改化合物的杯准,土成焓。化学反应的焓变:等于产物的总生成冶减去反应物的总生成焓:△,H。=∑v。A,H,。可以从一个温嚏的化学反应的焓变求得另一个温度的;教化学反I萱H的嗡变:A,H.。=△Ho+I△cPd7’。5.熵(Entropy):由热力学第二定俸所确定的{丸力状念参数,是系统中分子运动混乱程噎的寰萎.以S表示,单位J/K。熵是弧立系统自发改变状念的可能性的量度,改变看蛹增越大的状念,改变的可能眭越大。熵的微分定义是;系统在可逆过程中同外界交换的热量与系统热力学温受的比值豳:f墼1。Lr几堂标准熵So:25℃时处在杯准念的l摩尔物质的熵疽。杯准生成墒岱?:指在标准状态下,由处在标准奎的稳定单质在25℃时恒温生成处在标准仑的l摩尔化含物所引起的熵变,称为孛,下准,土成膊。6.吉布斯自由能[23l(GibbsFreeEnergyFunction):吉南斯自由能是状念函数。其定义表达式为G=H—TS。衷示体系总能量中具有作非体积功能力的那一部分能量。在25℃时处在丰,下准念的l肇尔物质的占伟斯自由能称为f,下准吉如斯自出能(Go)。标准生成自由能AG?:化台物的标准生成自由能是指在标准状态下,由处在标准态的稳定单质生成该化合物时自由能的变化,称为标准生成自由能。规定平衡状态下的稳定单质在任意温度下的标准生成自由能为零。7.化学反应平衡常数【24】(ChemicalReactionEquilibriumCcnstant)化学反应体系在~定的温度和压力下处于化学平衡状态时,各组分的活度不随时问lO中南人学硕十孚伊论文第一章氧化锅t产执力学数据昀鐾础俨究而改变。化学平衡常数表示处于平衡状念的子组分活度之『日J的数量关系。K越大,反应物的转化事越高。反应越完全。8.活度(Activity)及活度系数(ActivityCoefficients)活度,又称有效浓度,即单位体积电解质溶液中表现出柬的表观离子有敏浓度,即扣除了离子问相互作用的浓度。通常以口(Activity)表示。活度系数表示实际溶液与理想溶液的偏差,通常用r表示。2.2氧化铝生产中相关物质热力学数据的收集与甄别本氧化铝生产热力学数据库的基本数掘包括纯物喷热力学数扼和水溶液热力学数掘。纯物质热力学数掘主要柬自:1986年版B.H.6a6、IXIKHH等编著的{TEPMO丑IIHAMHKACH3HKATOB))1251,1991年版OKnacke编著的{ThermodynamicalPropertiesofInorganicSubstancesII》131和1993年版Barin编苫的{ThermodynamicalDataofPureSubstance))【26l。另外,以1985年版林传仙等编著的《矿物及有关化合物热力学数据手Iil}》f’71作为补充。另外,数据库中还输入了50余种复杂硅酸舱的估算数据,以供研究蕾参考。水溶液热力学数据主要综合了1982年皈DonnaldD,Wagman和WilliamH.Evcans编著的《111eNBSTablesofChemicalThermo#namicProperties,SelectedValuesforInorganicandCtandC2OrganicSubstancesinSIUnits))14j和((TEPMOjIHHAMHKACHnHKATOB》。日i昔是美国国家标准局于1965年丌始编辑的NBSTechnicalNote系列的修订和综合,可靠性较高。在甄别数掘时。原苏联《TEPMOⅡHHAMHKACH。qHKATOB))中的数据权威性较高.国内外大多采用其用于氧化宅吕生产热力学计算.故优先考虑。若某手册数掘与其它手册中的相差较大时,以多数情况为准。同时,把{ThermodynamicalDataofPureSubstance))和《矿物及有关化合物热力学数据手册》的数据作为补充,热容公式参考(《ThermodynamicalPropertiesofInorganicSubstancesII》。对于数据手册中没有热容公式的物质,数据库部进行了一一拟合,并通过热力学函数对拟合值进行返算,证明拟合结果正确后,将拟合值输入数据库中。从数据库中可查到的物质热力学性质(即字段)有:编号(NO)、分子式(formula)、相念(phase)、热容系数(a,b,c,d)、生成焓(/XHO)、生成熵(ASO)、标准焓(Ho)杯准熵(SO)、温度范围(TempUp,TempUown)、晶体转变温度(Ttr)、相变焓(Htr)、相变熵(Str)、熔化温度(Tm)、熔化焓(Hm)、熔化墒(Sm)、蒸发或升华温度(T_b)、蒸发焓(Hb)、蒸发熵(Sb)等。中南大学硕十学伊论文第一二章氧化锅生产热力学数据的蘑础研究目前,氧化铝,土产热力学数掘库已经收集了1199条记录。其中纯物璜为839条,水溶液为360条。2.3铝酸钠溶液活度系数计算模型的研究2.3,1概述氧化铝工业中铅酸钠溶液浓度普遍较高,而且由于铝酸钠溶液本身所具有一系列特殊的性质,如饱和溶液的超常稳定性,使之无法等同于其他电解喷溶液(如:纯NaOH溶液),活度系数研究的难度很大。同时,在铝酸钠溶液分解机理和溶液结构方面。学术界还存在分歧128羽】,也加大了改体系研究的难度。因此,尽管多年柬人们对铝酸钠溶腋体系热力学性质进行了大量的研究工作f30.31.3233】,获得了很多溶解寰.实验数掘。但是直接关于渗透系数和活度系数的研究非常少。1992年.Wesolowski应用Pitzer模型处理三水铝石溶锯度数厄,近似得到了一组有关NaAI(OHh的Pitzer参数和AI(OH)4一,OH‘和Na+混合的参数川。1999年.Park和Englezos采用等压法测定了NaOH.NaAI(OH)4·NaCI—H20体系从mNaOH=0.85,mmNfl=O.15(mol/kg)至0mNaOH=2.6,mNaA“OH)4=0.05。NaA]iOH}4一--0.35,mNacl=1.04(mol,7kg)在298,15K的渗透系数.并用Pitzer模型回归出Pitzer参数PSl。2002年.中南大学化学化工学院的周俊采用等压法测定了总苛胜碱从O.05mol/kg到12mol/kg,苛性比从1.64到5.53的铝酸钠溶液在313.2K的渗透系数,应用Pizter摸型回归得到NaOH.NaAI(OH)4一H20体系的Pitzer参数。其中一组313.2K时的Pitzer参数为:∥…haoH:13㈣NaoH:o.0240o.08669.13”lN。oH:O.31446,.o.00007367,c气aOH:o.003180.∥mNnl{0HH:0.03507,13“’NaAI(OHl4:O.081771,C~aAl(OHw:.O.001066,00H.Al(oH)4‘:t、|,Na+OH’AI(OH)4’:.O.01162t36l。这些实验研究对于铝酸钠溶液活宣系数的计算部有熏要的意义,特别是周俊获得的这组Pitzer参数使得准确计算40"C的NaOH.NaAI(OH)4一H20体系活度系数变得可行,所有相关结论也为以后的研究提供了极有价值的参考。但是,仅仅这些结果显然不足以支持使用Pitzer模型计算一定温度和浓度范围的铝酸钠溶液活度系数,大量其它温度的实验研究工作被期待。另外,有研究表明,现有的Pitzer模型并不适应于处理100"C以上的电解质溶液体系【371。很遗憾,由于以上原因,无法将Pitzer模型作为准备建立估算模型的基础。除直接的实验研究外,Debye.Hiickel模型在涉及含有NaAI(OH)4的溶液中使用比较常见。有人尝试用早期的扩展方程柬处理铝酸钠溶液13。j,但因参数值选取不准确以及数学模型上的缺陷,并没有给铝酸钠溶液活度系数的计算提供太大的帮助。我们注意到。中南人学硕十学伊论文第一章氧化锅生产执力学数据的垦础研究近年比较热门的Helgeson等人提出的Debye.HUckel扩展方程同样玻用于含有NaAI(OH)4的体系。例如:1996年.Diakonov等人在汁算地壳高温流体Na-AI—CI—OH中使用浚扩展方程计算了包含有NaAI(OH)t的混合溶液的平均活噎系数,忽略溶液中备离子a参数的不同,取均值为4A1391。另外,在1995年Pokrovskii和Helgenson关于A1,01.H、O.NaCI体系热力学性质的完整报告中,也旋刘应用该扩展镆型处理铝酸钠溶液。但荏随后讨论中假设OH’和A102。活噎系数相同,并没肓具体i步及活宦系数的计尊ItSl。综上,Debye—Hfickel模型在铝酸钠溶液活度系数的使用上也没有完整的解决方案。由于该模型简单,对于估算模型的建立有较强的吸引力。但是,如果要使用Debye.Hfickel模型作为估算馍型的基础,首先要给出针对NaOH.NaAI(OH)4.H20体系的Debye.H0ckel模型应用假设,同时给出使用Debye·Hilckel模型涉及的物理化学参数的具体值:再则,验证Deb、e—Hfickel模型对NaOH-NaAI(OH)4.H20体系的适用性,确定Deb、e.Hfickel模型町以应用于陔体系;完成上述工作,/j能基于此讨论具体建摸过程。2.3.2先前的计算模型及存在的不足李小破等曾基于一定范围内三水铝石在铝酸钠溶:夜中平衡溶解宦故掘,借助Debye.Hfickel模型,通过选取NaAI(OH)4合理的、#}£,引入”铝酸钠溶液的表观介电常数∥“的概念,将各种影响囚袤归结到£’,提出了铝酸钠溶液备组份活度系数的计算方:去ml。其所建立的活度系数和表观介电常数计算摸型为:’’[NaAIIOH)4=pl+p2’mh,aOHK+p3’mNaOHK‘+p4‘U.K+p5‘O.K2+p6‘lnmNaOHK+p7-lnaK(2-1)F’=pl+p2mNaOHT+p3‘mNaOHT‘十p4‘QK+p5‘aK2+p6。mNaOHT。aK+p7’mNaOHTz"QK+p8。mNaOHTaK。+p9·m~aOHTz·aK‘(2-2)通过上述公式可计算出铝酸钠溶液相关活度系数.但其中也存在一些问题:(1)根据平衡溶解度数据拟合的¨a^IIOHH与CtK、mNaoHT、T的关系模型对40"C时的覆盖性偏差较大:(2)苛性比(Ⅱk)在1.5~5.5和5.5~10.5范围内时,其,舌度系数的变化规律不一致:(3)在通过利用铝酸钠溶掖的表观介电常数占’建立的naA“oH,4计尊馍型,其方程较为繁琐,用其计算活度系数时,计算结果与实际值偏差较大:(4)由Debye.Hnckel模型可知,原s7的计算模型不能解决活度系数大于l的铝酸钠溶液的表观介电常数问题。因此,本文针对NaOH-NaAI(OH)4.H20体系,需解决:(1)探索适用于所有温度下统一、较精确的计算模型;(2)拓宽计算模型的应用范围;(3)提高计算模型计算的准确性,解决活度系数大于l不能计算的难题。!塑叁竺堕±!堡垒塞2.3.3铝酸钠溶液活度系数计算模型的建立兰=三童墼坐堡竺兰竺垄竺垫堡塑堡!|!!竺塑建立NaOH.NaAI(OHh.H20体系活度系数}十算模型豹主要思路妞图2-1所示:l平一婶餐¨卜Yk…Ⅲ..的计算l与文献Y”。..债比牧l1r平衡篆竹FYu。。计算礁掣的拟台I拟台占,与a。、m,,的函数天系Ity。。。汁劳硅掣图2-1Fig.2-1确立孙。n‘)HH估算接型曲虎程堡Flowsheetofbuildingmodelfor丫N“I.oH“2.3.3.1平衡条件下NaOH.NaAI(OH)4.H20体系活度系数计算模型根据文献报道的三水铝石溶解赛数据【4142’,拟合了在乎衡状叁下。铝酸钠溶液活宣系数与温宦、a。和m的函数关系式:丫N直A“oHw=pl+p2’肌+p3’脚二+p4’吼+p5’也2+p6’lnm+p7’lna,(2·3)式(2-3)中不同温度的调整参数值,详见文献¨引。其中40。C时的调骜参数为:pl—1.2557,p2=0.OJ7t.p320,0045,p4=一1.1159.p5=0,0377,p6=一00563-p7--4.0364。由于二氧化铝生产过程中备组份一般处在非平衡、非杯准状念下.因此需将公式(2-3)的适嗣渔圈进行拓展。2.3.3.2非平衡条件下NaOH.NaAl(Olt)4.H20体系活度系数计算模型文献曾报道40‘C时用Pitzer方程研究的铝酸钠溶液活度系数数据㈣,溶液瓯<5.m=0.01~10mol/kg。根据各温度下、备种吼下的平衡溶解度数据,利用公式(2.3)计算出相应的活度系数。将其与40’C时对应的活度系数迸行对比,然后计算出对应活度系致的差值,将这些差值进行拟合回归。得到活度系数差值(△Y)关于温度(t)的函数。所得到的函数形式为:中南人学硕十学伊论文第一章氧化铝生产热力等’数据的基础研究△YI=一O.0682+O.0038+(t-40)△托=O.0888一O.0045+(t-40)(口。S5.5)(吼>5.5)(2-4)(2.5)于是可得到非平衡条件下活度系数(T)关于温度(t)的函数:丫l=丫40一△YI(或△Y2)(2.6)100℃(t≠40℃),在公式(2.4)、(2-5)中,t为摄氏温度(℃)。其适用范围为:v=25吼=1.5~lO.5,伪20.Ol~lOmol/kg。根掘式(2-4)、(2-5)中相关参数盯知:当在低苛性比(口。S5.5)时.i舌度系故碹温宦升高而降低:当在高苛性比(吼>5.5)时.活度系数随.盛宦升高而增々Ⅱ。这可能与铝酸钠溶液中铝酸根离子的赋存状念有关,同时这一变化觇律与文献【蚰l的结论基本一致。2.3.3.3补,AIIOHH的计算模型基于非平衡条件下溶液中备组份平均活度系数的值(公式(2-6)),借助于Debye·Htickel活宦系故修正公式,假设NaAI(OH)4溶液为连续相,引入铝酸钠溶液的表现介电常数f’这一概含。将各种影唰归结于占’的变化f”I。计算对应的占’值。计算£’值的模型公式如下式:!!;!!+兰l!:!旦竺!。:+^.1¨一五I18高0706箸I丢丁s(2,7)×f’:“1‘二1式(:一7)中-P表示锥j酸钠溶液的的密室(g-cm。);肌为Na‘喷量肇尔浓霞(mol-kg“)。m。m~“Ⅲw,.+卅Ⅶw:T为绝对温度(K)。不同口。范围F’与7_和聊的函数关系,得到函数(2-8)、(2-9)。(1)、当瓴=1.5~5.5(含5.5)时,占’=bl+b2+ln(,)+b3*m+b4·朋2+b5*lna,+b6*m2*lnak/t15十b7‘聊3*In口。(2)、当口。=5.5~10,5(不含5.5)时,(2-8)一=bl+b2+t+b3+f2+b4+m+b5+肌2+b6+d^+b7+《+b8+用+a,h/t’5+b0·肌’吼不同吼下的调整参数值,列于表2-1。表2-1不同吼下方程(2-8)和(2-9)中的调整参数值Table2-1Values(2-9)ofadjustableparametersatdifferentcausticity.inEqnst2-8landt2-9、其中,式(2—8)、(2·9)的适用范围为:t=25~100"C,m---0.1~lOmol/kg。研究证明,当式(2—8)、(2—9)中的,分导U5829.47和5360.30时,可不考虑修正项系数打,此时br--O。中南人学硕十学伊论文第_二章氧化铝生产执力学数据的基础研究由式(2.7)可知,当活度系数大于1时,s’无实数锯。另外。当由(2-8)、(2-9)计算出的占’大于一定位时,修正项系数孵O。此时需汁尊打的值。经拟合,打足温室(,)、苛性比l口.)和苛碱浓宦(m)的函数,其函数形式为:(1)、当口。=1.5~5,5(合5.5),由式(2-8)计算出的s7>829.47时,br=cl+c2+,+c3+,+c4*m+c5*m2+c6+口。+c7+口:+c8+m’口。(2)、当口。-5.5~lO.5(不含5.5),由式(2-9)计算出的一>36030时.br=cl+c2+f+c34r,:/+c4+口。+c5+口:+c6+,以2/t+c7+mz*t(2-11)(2-10)不同a.下的调整参数值,列于表2.2。表2-2不同口。下方程(2-10)和I2-lI)中的调鏊参数值Table2-2Vahzesofadjuszableparametersatdifferen!ca).)slzcil"yinEqnsr2一JO)andt2-11)由此,结合Debye,HOcke[公式榴应参数的取值.我们可以借助于段模型得到NaOH.NaAI(OFI)4.H20体系活度系数计算摸型,如公式【2-7)所示。其中,当∥1829.47(瓯=1.5~5.5)或一<360.30(q=5.5~t0.5)时,6产O,占’可由式(2-8)、(2-9)计算而得;当s’>829.47(口。=1.5~5,5)或s’>360.30(口。=5.5~10.5)时,6辟0,b,可由式(2-10)、(2-II)计算雨得。对于不同的占’值,可进行如下解释l“l:1)口。=}5~5.5:当tr/和a.一定时,f’篮随,显度的升高两减小。这可能是出于随蕾温童升高,溶液体积嘭张,每单位体积分子数.成少引起的:当研和温叟相同时,F’值随q的增大而增大,这可能是出于随着口。的增大,溶液中体积较小的OH一离子增多、体积较大的AI(OHh"离子减少导致溶液单位体积所含溶质数增加昕引起的:当温宦和吼相同对,占’值短m的增大蔚增大。这可能是出于/.17趵增大导致溶液草位体积中所含溶质数增加所引起的。.2)瓯--5。5~10.5:当rt'/租aj一定时,s’值随温度的升高而增大。这可能是在高∞溶液中,体系中主要存在AI(OH)4.、AI(OH)6}、OH’等离子。当温宦升高时,溶液中体积较小的OH一离子增多、体积较大的AI(OH)4.离子减少导致溶液单位体积所含溶质数增加所引起的:当m和温度相同时,s’值随瓯的增大而增大;当温度和吼相同时,s’值随加的增大而增大。2.3.4适用性分析由公式(2-8)、(2-9)可知:运用上述NaAI(OH)4,NaOH—H20体系活度系数计算方中南人中项}学伊论文第一帝氧化铝t产执力学数据的基础研究法,计算出t=40"C,口。=l,5、3、5时不同l浓噎条件下的活度系数。将计算结累与文献p6l数据比较,比较情况见图2.2。由图2,2可以看出,在吼分剐为1.5、3.5时,经拟合馍型汁算出m=O.1~10mol/kg的活度系数与文献所列活度系数的相对误差clO%。计算值和文献值均表明:在低浓度时,活噎系数陋浓度的升高而降低:在高浓度时,活度系数谪浓度的升高而升高。计算数掂与文献数掘的吻合情况较好,说明浚模型对于计算铝酸钠溶液的活度系数育很好的适用性,这为热力学上分析铝酸钠体系各组份的行为提供了理论依掘。lOn∞驰叼l00090080>O000孙∞∞加070060O50二4¨,I二m,(tool·k91)/17/(mol·kg。‘)Ill0^_0000加m∞帅舳加的∞02468l012m/(mot·kg‘)图2-2估算模型丫Nn¨oHw计算值与文献值的比较Fig.2-2Comparisonof"fMHoH*calculatedbythecalculationmode!hthisstudy(I)∞d抽ereference(II12.4铝酸钠溶液中Na2C03、Na2S04、Na2H2Si04活度系数计算模型工业铝酸钠溶液中的Na20,除结合成铝酸钠(NaAl02)和NaOH的形态存在外,还以Na2C03,Na2S04,Na2HzSi04等形态存在,其活度的大小对于体系的热力学性质有较大的影响。研究其活度系数计算模型对于体系的热力学计算育一定的理论指导意义。±塑叁竺堡±芏芝堡壅2.4.I计算模型的建立2.4.1.I.笙三兰一兰些!堕羔!笠!兰鳖塑塑墨型竺生Na2C03活度系数计算模型的建立Na2C03—2Na++C012‘(2-1:)在溶液中,NazC03存在如下电离平衡:反应(2.12)的平衡常数为:足:华;‰。×席乞.一Ⅲt2"Ⅶ2m(2-13)式(2.13)中,哎,#.、口№.、口Ⅶ∥^分别为C032",Na+、Na2C03的活度;r№:cfJt为N永03的活度系数;棚:。.、m。,}分别表示Na+、C032-的重量哮尔浓度,单位tool出式(2.13)可得:kg"‘。‰:。,=(当反应(2-12)达到平衡时,(2—14)聊j口+。m‘rl△GT=一2.303RTIogK即:logK=△Gr/·2303RT(2-15)(2-16)经计算,反应(2-12)的△GT为:14934T二+I1.004‘In(T)4T(2_j7)/NGT;.2406373127-1224000/T-26.4213408’T+0.074至此,通过式(2.14)、(2-16)和(2.17)便可计算出某~温度平衡条件下Na2C03的活度系数。根据A.A.阿硌拉诺夫斯基提供的各种温嚏下纯喊在铝酸钠溶液种的溶解度数据∽1·用上述方法计算出了Na2C03相应的活度系数。如表2-3听示:表2-3不同条件下Na2cO,的,舌度系数里坐!::垒垡:i!:竺笪!塑!罂!兰g∑型壁竺竺生堕———一959595959548S411.262.051.73I.511.391.361.251.49I.491,49I.491.491,491.49226001430077.305300424022.4519.7584.4674.6665.9158.2251.614741O.1J250.1533020030.12880.11510.13070.12230.0385O.03920.03950.039600393O,03800.03652502j96jl4V,,·Vqu.u“7003080.02880.0296O.03030.0308O.03160.03t6003180.0312O·0305O.03010,0293O.02884025025030030030030030030030030030030030025747275.781083712004132.0414a,13156.9116994183.68197.85212.21227.15242.28I.49IA.93293.293.293.293293.293.293.293.293.293.2935.9l33788615774869.8761,53560450.4147.0744.2340.6738.2635,1383.50206202425026800278.00118.13131.41145.13159.34174.0918965205.75959525025025025025025025043.72翌!!!::!!.!兰—2羔——旦旦翌111:!竺:!!‘●_■■_■●_●■■■_■______■’__●____■__-●■■■■●■■●■-●_■■■■■■■■■■●_____●■■■-●■__■■____。‘■●。。。’———一一8中商人学硕十学伊论丈第一:章氧化铝t^执力学数据的基础枰窀然后通过数学回p。1的方法。拟合活度系数(TNa2C03)关于温度(T)苛性比(ak)、全碱浓度(mNa*IT))和碳碱浓度(mc0,2。)的函数。所得到的函数彤式为:logr、。,fJ,=pl+p2+InT+p3/T+p‘18logmNa+I.1+p5+lognk-032‘(2.18)其中,T为丌氏温宦(K),mN。+㈤为全碱浓度(moI·kg。’),mCO;2-为碳碱浓变(moI·kg‘’)pl~p5为调整参数,其值见表2—4。表2.4方徉2-is中的阔帮参数值Table2-4ValuesofadjustableparametersinEqnf:一f8】2.4.I,2Na:S04活度系数计算模型的建立运用蹩似fNa2COsi舌宦系数汁尊模犁逭芷的思路.采fflA.A.阿格拉若夫斯基堤哄的各种:益度下Na:SO{在铝酸钠溶液种的溶解度数掘l刈.计算出不同温度、不同碱,衣熏下Na2S04的活度系数,然后通过数学回归,拟合得到Na2SO。活度系数(ha2s04)关于温室‘T)苛性比{瓢)、全碱浓度{m。。+m)和硫酸碱浓度(11"15042。)的函数。所得到的函数形式为:logr~。:.v‘=kl+k2’In(T)十k3’logmha+I十寸k4+logms04z’(2.19)其中,T为丌氏温嚏(K),m、。+(tl为全碱浓度(mol·kg。),mⅫ2‘为硫酸碱浓度(mol-kg’‘),kl~k4为调整参数,其值见表2-5。表2-5方程2-19中的调整参数值Table2-5ValuesofadiustableparametersinEqnf:·19)2.4,1.3Na2H2Si04活度系数计算模型的建立在铝酸钠溶液中,sio:的溶解度比较小(~5叽),对于Si02来说,为稀溶液体系。故可直接运用Debye-Htickel模型计算其活度系数。log元:一掣(2-20)1+dB.1l“式(2-20)9tI表示电解质溶液的离子强度,单位为tool·kg一:乙、Z.分别表示f负9离子的价数:a是关于外离子与中心离子极限距离的参数,单位A;Ar、B,表示Debye.H0ckel参数,单位分别为kgln.mol地、kgV2.mol-U2.cm~,定义为:铲警岛=警f2—21)(2-22)其中,T裹,--,一.:…19r@K),P,p分别表示水的介电常数及畜宦(g·cm’3)。参数a是Debye-Hfickel公式应用的关淀参数。一股而言,a值大于晶体的离子半径之和,小于其水化半径之和。目4≯通常的做法。使用FB离子水化半径和阴离子晶体半径的加和。即Na2H2Si04的a值为:Na+离子水化半径(1.91A)与H2Si042.离子晶体半径(1.61A}相拓昕得I.91A+I.6lA=3.52A。2.4.2适用性分析运用式(2一l8)所述Na2C03活熏系数计算方法。计算不同温变、不同苛碱浓度时Na2C03在铝酸钠溶液中的平衡溶船度。将之与文献数扼州比较,比较隋况见表2-6。Fig.2-6ThecomparisonofsolubtlityofNa2C01calculatedbythemodelin廿IlsstudyandthereferenceT(℃)Na20kL∥L)ak—T(℃)表2-6Na,Oc(g/LlNa2COt平钎毒祥度计算埴与走欲值的比较Na:Ot,(g/L)文献。直22600143.0077.30530042.40224519.7584.4674计算值226.05142.6277.2952.7942,2022.39196994,6774,3766.11Sg4I5I7847.5743,85aL———————二—一文献一直计雉值1.491.49lNa,Ocfg/L)959595959595952502502502502502502504851.262.051.73I.5I1.391.36I.25l49I4925025025030030030030030023963257.47275,781083712004132.04144.13I56.9l1699418368197.85212.2I227.15242.283794359133,78861577.4869876I.5356.04504f47.0744.23406738.2635.1338.0335.97338285.6277.026吼476I.2055.7450.16468243.9940.4538.0434.934l4083.50206.20242.502680027800118.1313l4f1451315934174.0918965205.75493293.293,293.293293293293.293.29329:29663∞3003003003003001.491.49I.4965.9I58.225I6147.4l43.72I.491.49中南人学硕十学伊论文第_二蕈氧化锅士一L媳力学数据的墓础研究由麦2-6可以看出,数据的吻合情况较好,说明了该计算模型对于铝酸钠溶液中Na2C03活度系数的计算具有较好的适用性。此外,运用式(2.19)所述Na2s04活度系数计算方法,计算出在不同温度和不同苛碱浓嚏条件下铝酸钠溶i夜中Na2S04的活度系数,其数据与RobenoNa!SO。的活度系数数据m1摹本吻合,证明了本模犁的准确性。tPabalan等汁算的2.5本章小结l、与热力学相关的基本概念状态函数(StateFunction)、标准状奁(StandardStates),热容(HeatCapacity)、焓(Enthalpy)、熵(Entropy)、吉1{】斯自由能.(Gibbs常故K【ChemicalCoefncients)。ReactionEquilibriumFreeEnergyFunction}、化学反应平衡Constant)、活度(Activit)’)及活度系gt(Activity2、氧化铝生产中相关物质热力学数据的收集与甄别目日F.氧化铝生产媳力学数据唪已经收集了1199条记录。其中纯物质为839条,水溶液为360条。从数据库中可查到的物质热力学性质(即字段)有;编号(NO)、分子式(formula)、相念‘phase)、热容系数(a,b,c,d)、生成焓(AHO)、生成墒(△SO)、标准焓(HO)标准墒(SO)、温度范围(TempUp,TempDown)、晶体转变温度(T仃)、相变焓(Htr)、相交畸(Str),熔化温宣(Tm)、熔化焓(Hm)、熔化墒(Sm)、蒸发或升华温度(Tb)、蒸发冶(Hb)、蒸发熵(Sb)等。3、铝酸钠溶液活宦系数计算模型的研究1.8248lgy2一i—.4—-i—8—:—0719三6jI苓+6,,×106(pm)‘2l1.×(巴)…:p表示铝酸钠溶液的的密度(g-cm。);m为Na+质量唪尔浓度(mo|·kg"1).m=肘、。Ⅲwh+肌mw;r为绝对温度(K),占’为铝酸钠溶液表观介电常数,其与r和m的函数关系洋见2.3节。4,铝酸钠溶液中Na2C03、Na2S04和^fa2H2Si04活度系数计算模型的建立。主塑叁兰堡±堂竺堡奎墨三垦—墼些塑兰兰兰苎堕垡墅型王塑些些三塞.生第三章氧化铝生产热力学数据库的优化与完善3.1概述计算机软件的丌发,属于软件工程学的麓畴,人仃】已经建立了一整套t见范化的方法柬指导软件的丌发,而这些方法同样适胃于软件的优化或再丌发过程14”。软件性能的优化,在理论上是无穷函数,任何软件都有进行优化的空JlBJ。软件工程与数掘库理论的良好结合,是数据库软件丌发及优化的关键。数据库应用软件的优化工作涉及众多方面的内容,需要根掘具体情况而定。数掘库性能优化的一些主要内容,列于表3.1m'491。袁3-I敷据库陛能优化方向Table3·1Optimizingaspectsfordatabasecapabdity执化方向描述数据咋没计是堤毒鲁体性能的最有效建释。数据咩设计包括逻辑数据『车絮卡勾(船丧数据『芊设汁萃|f约柬}以及扬理r}能(缸衄嚣彖冼,对象伊萏和袁引'齐向优化府H{榨序殴汁优化实心I冗和l府H{数据『芊系统提供的实心f只和1.只的一些可_【f{选项提高数据辜性能I冗性能为了提高帮体性能.更改操仃系统(MicrosoftWindowsNT、MicrosoftWindows优化搬务器性能或MicrosoftWindows20001霜『数据唪(SQLSever或Access)98使圳lL确:殳计的裔,旬能够显#堤扁故据7莘眭能止确c殳计Ⅲ户样序,事务边界,锁定_-批处理的使川部雷丁府Hj群序殴计7旺镉我们注意到,除使用结构化,规范化的方法外,很多优化工作最终将归结到代码的编写上,属于编译时优化的7苞畴。对此,编程人员往,芏有一些经验性的技巧,这些编写代码的技巧对应用软件优化有很大的帮助。例如:对VB应用程序而言,可采用减少控件数目、用杯签控件代替文本框、共享图片等方法实现程序大小的优化:选用最佳的变量、数据类型实现速度的优化等。中南入学谚十学竹论文镣二章氧化锅t产扎力学数据牢的优化与完善3.2旧版本数据库存在的问题及相应优化解决方案旧版《氧化铝生产热力学数花库》在过去的使用中表现出了强大的生命力,刚时也暴露出了~些问题。新版的《氧化铝生产热力学数拢库》梅疰改善整个数嬲摩性能上做出尽可能多的努力.并勺口入新的热力学计算模型。藏版设计目标旧坂相同,即:桌面型、}:}研用专业数据库应用软件。这砷产品定也}要是由实际需求所决定的。眩软件属Fj史奇软件.主要用F科研工作。用户群为个人或,J、闭体.FJ时具何专业性强.数据表简单、数瓤晕不欠≮特点。由于个人计算机性能价格比的迅速提高.给了桌面犁数据库以极大的发展空删,也为我们的软件丌发和使用提供了充分的技术支持。下文将对旧版数据库存在的主要问题及对应的优化完善方案和新增内容作分类介绍,脚时为了增加本文的可溃性弓简洁性,省略了数据『聿.程序源代码优化过程的的jf细介绍。1.1关于精简性1.存在的问题在=致掘章的某些界面和J力能中.存在苦:吾三冗长、J力能节复的现象。如数掘计算界图3-1旧版数据库敷据计算界面F|g.3-1ThecaJcuJaZzonmoduleofformeredition3面中胡“纯扬;专(离产)在温瞧T下的焓Ht礴st厦生成}1出自}已经包含了”离子在温熏T下的焓Ht墒St硬屯咙自由能”白勺计雉』由能.如罔3.I昕示:中南大学硕十学1奇论文第二章氧化锅生产拽力学数据车的优化与完善2解决方案针对旧版本数据库在界面方面存在的问题.在不减少功能和影响操作的前提下,对数据库界面进行了简化,对已有的重复功能进行了整合,对其中的浯言进行了精简。以数据计算模块为例,优化后的界面如图3-2所示:图3-2新版数据库数据计簟界面Fig,3·2ThecalculationmoduleofneweditJ(}n3.2.2关于完整性1.先酊数据库基础数据覆盖的内容体系帽图:AlzO,·si02系、CaO-Ah03系、CaO-Si02系、Na20.A1203系,Na20.Si02系。数据库物质记录数:Data表953条,IblAqueous表233条,tblPm'ea表720条。铝酸钠活度系数估算模型适用范围:温度t--40"C~100℃,苛性比ctK=I.5~5,苛碱浓度mHaoHg--05~10mol/kg。2.目前数据库基础数据包含的内容体系帽图:Al∞3一SiOz系、CaO-A120j系、CaO.gi02系、caO.Fe203系、CaO-FeO系、FeO—SiOzg,Na20-A1203-H20系、CaO—A120rSi02系、CaO—FeO-Si02系、NazO.A1203-SiO:一H20系、Na20一CaO—A120rSi02-H20系、Na20一CaO.A1203·C02.H20系。数据库物质记是数:目前,Data表1199条记录。其中纯物质为839条,水溶液为360象。铝酸钠活度系数估算模型适用范围:温度t=25"C~300"C,苛性比aK=1.5~10.5,苛碱浓度rnH目赇=0.1~10mol-kg.‘。±妻叁堂堡±兰生堡壅苎三!墨!兰塑竺兰兰塑!!堕!竺垡!!生圭堂3.2.3关于反应方程式的配平1.存在的问题配平后的反应方程式的系数采用小数表示,这可能给后续计算带束一定误差。如式(3-1)所示:(O.67)Fe+(1)C12=(067)FeCl3(3一1)另外:在某些方程式的配平请求中,程序无法完成配平且弹出出错提示。2优化方案为解决方程式配平中存在的问题,在方程式热力学数掘计算模块中添加了“手动输入系数”和“自动配平系数”功能模块供用户选择,以解决可能出现的问题。其界面如图3—3所示:图3.3反应的AGt计算界面图F嘻3-3TheformsofAGtofreactions3.2.4关于准确性1.存在的问题经原版数据库估算出来的数据,部分值与文献值、手动计算值误差较大,这影响到了数据库软件的应用推广价值。2.优化方案一是对原有的基础数据进行了一一核对,一保证原始数掘的准确性:二是增加了物质的△H298、AS298、AG298等热力学数据.以扩大物质热力学数据的查询范围及增加其计算方式:三是对于有离子和纯物质参与的化学反应,所有物质均统一1采用了式(1-5)、(1石)和(1-7)的计算模型,保证热力学数据的一致性,增加计算结果的准确性。中南人宁硕十学忙论文第二章氧化锅士产热力'-7数据『芊的优化与完簿3.2.5关于程序非正常出错原版数据库程序在运行过程中,有时会突然中断运行并非J下常退出主程序。若需继续操作,需重新启动程序。如图3-4所示:囿3_4数据库程序出错界面图Fig.3-4Theerrorsofformerediuon2|优化方褰通过记录程序出错时的错误代码,找出错误原因,修改相应的原程序并调试,防止程序出错非一常退出现象的缆生。优化思路如图3-5所示:困3-5敷据库程序出错优化思路图Fig.3-5Theoptimizedmca¥ureforthecrror¥ofthedatabase3.2.6关于数据库操作方法1存在的不足数据库包括了3个查询模块(分子式、元素、体系查询),2个类别(水溶液、纯物质)的3种数据编辑方式(更新、添加、删除)和9种数据计算类型。但未对其中的查询方式、输入规则、操作流程等进行必要的说明,对初次使用者来浣有一定难度。中南夫学硕十学佛论文第二章氧化锅士产媳力学数据『芊的优化与完善2解决方案为了使耳】户可以全面系统地了解系统的应用』力能.并学习到其基本使用厅法,特编制了~月余字的《氧化铝生产热力学数据庠(TDAP)的用户使用手册》。从数掂库系统的安装建立到各部分、备模块的使用方法和具体操作步僳都进行了详细说明。如复杂硅酸盟热力学性质拆分估算模块的使用说明如下:图..1-6复杂硅蝮盐热力学性质的折分估霉界面固Fig.玉6Theinterfaceofthermodynamieestimalionofcomplicatesilicates操作步骤:1.输入待估算物质的分子式;2.在酸碱比合物化合产物框中,按照复杂硅酸娥中各氧化物的酸碱性强弱,由强至弱输入各个氧化物的两两化合产物及其系数(系数输入到稼色栏);例如:求CaO·2Si02·A1203·2H20在酸碱化合物化含产物框中应输入:l1CaO*SiO:AI203+sio:注意:此处元素应区分大小写。3.输入计算温度或温度范围;4,点击“计算”按钮进行计算。计算过程中,程序会询问用户某一确定物质的具体晶型或状态,用户根据需要选择即可。选择完毕后将显示计算结果;5.查阅计算结果后,用户点击“完毕”按钮可返回计算窗体:6.若要继续计算.单击“计算下一个”即可。中南大学硕十学伊论文第二章氧化锅t产执力学数据庠的优化与完善3.3其它方面的完善3.3.1元素精确查询原版数据库在进行元素查询时,不是按照元素进行查找。面是按照字母查找,即将包含该元素字母的所有物质均查询出来,即模糊查询。如查询元素“S”时,物质“Si02"也随之被捡索出来。优化时运用“If,…一Then”语句对用户所要查询的元素进行了筛选,达到精确查询的目的,从而解决了类似的问题。3.3.2安全性为保.匹数据库的安全性,设冒了数据库安全密玛.包括数据库登陆啻码、数掘库修改添加密码、Access表打丌密码等。3.3,3程序说明文挡为了给以后的程序完善和升级丌发提供方便,制定了6000余字的数据库源代码说明文档,对在丌发程序中采用的算法逻辑、层次结构、模块『日J的连接、重要变量的使用等情况进行详细的说明。3.3,4数据库用户手{j{}为了方便用户使用数据库.以图文的形式编制了10000余字的《氧化铝生产热力学数据库用户手册》,对数掘库的基本框架、主要功能、使用操作方法、注意事项等情况进行了详细的说明。3.3.5帮助氧化铝生产热力学数据库提供了标准Windows风格的帮助系统,以帮助用户进一步了解和使用本数据库系统。帮助系统的目录是按照数据库的内容层层展开的。当用户进入某个功能模块后,悠获得有关帮助,请直接点击窗口左下角的“帮助”按钮,获得相关帮助内容;也可根据需要在目录中进行选择,或者直接输入关键字查找到所需主题,如图3.7(a、b)所示:中商夫学硕十学伊论文第三章氧化锅生产热力学数据痒的优化与完善Ca)图3-7数据库的帮助系统(a_目录式帮助F嘻3-7b.捡索式帮助)Thehelpsystemofthedatabaseb·Searchesmode)(a-Direclorymode本数据晖的帮助文件是用VisualWorkshop编写的。编写帮助文件的步骤为:Basic6.0所携带的号门编写帮助系统的工具HeI口1、创建帮助主题文件(+.rtf)。帮助主题文件保存了所有的帮助信息,它就象一本书,包括目录和~个个章节,各章节之间通过热区束链援。热区可以是文字或图标,实现主题J日J的跳转。该文件可用Microsoftnf格式。2、创建帮助目录文件(+.cnt)。帮助目录文件是ASCII文件。在帮助主题显示Wbrd编写。保存成框中以{对状目录的形式向用户展示帮助内容,同时在-索引’条中用户可通过关键字进行内容查询。启动应用程序hcw.exc即HelpWorkshop可编写帮助目录文件。3、创建帮助项目文件(+.hpj)。项目文件是用束编译帮助文件的ASCII文件,它将帮助主题文件、目录文件及相关参数传递给帮助系统4、编译和生成帮助文件(·.hip)。直接点击界面上的‘编泽’按钮,该应用程序可自动生成帮助文件。在VB6中通过API函数可对生成的帮助文件实行调用。语法如下:DeclareFunctionWinHelpLib‘‘User32”Alias“WinHelpA”(ByvalhwrldAslong.ByvallpHelpFileAsString,ByvalwCommandAsLong,ByvaldwDataAslong)Aslong29中南大学硕十学付论文第二章氧化铝生产热力学数据『芊的优化与完善3.4本章小结本章从用户界面的精简性、计算结果的准确性、热力学数掘的完整性、反应配平程序,铝酸钠活度系数估算模型的建立等方面角度分析了旧版数据库存在的主要问题,并针对相应问题提出具体的解决方案,论述了具体的优化过程以及对优化的评述,同时进行了新旧数据库的使用对比,晚明新版数掘库的优化成果及明显技术优势。为方便数据库的使用和升级,编制了《氧化铝生产热力学数据唪用户手册》、《氧化铝生产热力学数据库程『芋说明文档》和帮助系统。中南大学硕十学竹论文第四章Na!O.A1203-H:O系热力学分析第四章Na20.A1203.H20系热力学分析不言’而喻,要深入研究氧化铝生产过程中内在的化学规律,用理论柬指导生产非常重要。热力学是研究各种变化过程能量转化关系,以及过程进行的方向和限宦的科学。热力学能为实验研究指明方向..避免实验工作中的盲目性,对优化生产条件、提高生产效率、降低能耗及原材料消耗育重要的指导意义【5051·525”,从而为发展新工艺和改善老工艺以及为有预见性地控制现育生产提供理沦依掘I“l。玎耳法是碱法生产氧化铝中的其中一种i虎程。其是直接利用含有大量游离苛性碱的循环母液处理铝土矿,溶出其中氧化铝得到镭酸钠溶液的,往铝酸钠溶液中添加氢氧化铝(晶种)经长时问搅拌便可分解折出氢氧化铝结晶【5”。拜耳法用在处理低硅铝土矿,特别是用在处理三水铝石型铝土矿时,流程简单能牦低,作业方便,产品质量高,其经济效果远非其它方法所能媲荚。目前全世界生产的氧化铝和氢氧化铝,有90%以上是用玎耳法生产的i”;。4.1三水铝石溶出过程热力学分析玎耳法的实质就是使下一反应在不同的餐件下朝不同的方向交替进行f55】:A1203(1或3)H20+2NaOH+aq一2NaAl(OHh+aq(4一1)铝酸钠溶液在氧化铝工业生产中扮演着非常重要的角色。在通常的热力学讨论中均为考虑活度对反应吉南额自出能的影响。然而铝酸钠溶液活度系数的计算不仅有助于氧化铝生产理论的完善,而且将大大促动氧化铝生产新工艺的研究和丌=发,达到降低成本、节能降耗的目的。现碾据已建立的铝酸钠活度系数计算模型15”,以三水铝石在NaOH溶液中的溶出反应为例,对其进行热力学分析。三水铝石在NaOH溶液中的反应可表示为:该反应的反应吉却斯自由能为:△嘶一f+器,s专≯2△旰+罴ts警.4t(OH),4-NaOH=NaAI(OH)4(4.2)c㈣式(3-3)中,蛳aAI{oHH,(]NaOH分别为NaAI(OH)4、NaOH的活度;YNaAI(oH)4、¨a0H分别为NaAI(OH)4、NaOH的活度系数;mAI(0H斗、moH.分别表示AI(OH)4.、OH’的重量擘±堕盔兰堕±兰堡垒壅尔浓度,单位molkg"‘。笪堕兰型!:竺:垒k旦!:旦!Q墨垫垄竺坌塑根据式(3.3)计算出三水铝石不同温度时在NaOH溶液中的反应吉钿斯自由能值(Na20k=2309/L,ak-2.3),同时将其与不考虑活度时的吉却斯自由能值进行了比较·如图4-1所示:80008000{4000昌20000-2000250270290310330丁,’K350370390410图4-l三水铝石瘩出反应的AG-T图(Na20k=2309/L.OK=2.3)Fig.4-1ThefiguresofAG·Tofthedigestionofalumina(NazOt22309/L.ak22.3)由图4.1计算可以看出,三水铝石在其本身热力学悬度范围内和NaOH反应的吉11】斯自由能值随着温度的升高而减小:到达一定温度后,反应的吉南斯自由能值为负值,说明反应在一定温度时能自发进行。同时,对比考虑活童和不考虑活度时反应吉布斯自由能值的大小可知,活度的引入使反应的吉布斯自由能值有所降低。在相同aK(aK=2.3),不同苛碱浓度下三水铝石溶出反应的△G-T图如图4-2所示:趸2≥t昌042—4280300320图4-2Fig.4-2Na:ok对三水铝石溶出反应AG-T的影响(O.K=2.3)ollEffectofNa20kAG-Tofthedigestionofalumina(aK---2.3)中商人学够十’宇侍论文第四章Na:O.A1201·H10系热力学分析从图4—2可知,在aK不变的条件下,增加苛碱浓度,三水铝石溶出反应的吉南斯自由能值减小。即是说,当aK恒定时,苛碱浓度的增加将有利于溶出反应的进行。在相同苛碱浓度(NazOk=2309/L),不同O.K下三水铝石溶出反应的AG.T图如图4-3所示:星≥t<图4-3Fig.4_3aK对三水铝石瘩出反应△G—T的影响(Na2Cll,=2309/L)onEffectofaKAG—Tofthedigestionofalumina(Na:Ol,=2309/L)从图4.3可知,在苛碱浓度不变的条件下,随蕾旺K值的增加.三水铝石溶出反应的吉白斯自由能值碱,j、。即是浣.当苛碱浓度恒定时,口。值的增加将更有利于溶出反应的进行,进入溶液中的A1203的数量随之增加,文献报道也证实了这一结论㈣。4.2三水铝石平衡溶解度估算铝酸钠溶液是氧化铝生产的基础,有许多研究者曾对在各种温度下的NazO.AIz03.H20系的平衡情况进行过研究,得到的结果大同小异。由于铝酸钠溶液活度系数方面研究的限制,过去人们在计算NazO.鲇20,.H20体系A1203溶解度时,往往以浓度代替活度,忽略活度系数的影响。数据库提供了加入活度系数的平衡溶解度计算模块,束解决这一问题。这些平衡溶解度计算结果对于工业生产有一定指导意义,同时也可以对活度系数计算模型进行检验。数据库计算出的在30~100"C时,A120)在NaOH溶液中的平衡溶解度数据,见表4.I:状念图,如图4.4。表4-l计算出的三水铝石瘩解度数据N赴o(叽’可鬲—丽百—面丁1秽地告下、丽厂丽百面Table4-1GibbsitesolubilitydatacalculatedbythemethodjnthisstudyNaz0.¥/L图“不同温度下的Na20.A1203-H:0系状态图Fig.4-4ThesolubilityCurvesofNa20一A1203-H20systematdifferenttemperatures从图中可以清楚地看到,随温度的升高,等温曲线的曲率逐渐减小,即越来越拉直。这一规律符合一般常识,从而在一定程度上证明了论文所提出的活度计算模型的正确性。34中南人学硕十宁伊论文第四章Na20-AI:03-H20系热力学分析进一步,将不同温度叮A1203在苛碱中的平衡溶解度数掘与文献数据比较,文献选取Russell的数掘【”I,以40"C、70"C、100"C为例.比较隋况见图4-5~4.7。160140120≤100∞占80;604020005000l50Na20.gL200230300图4.540'C平衡露解度计算结襄与丈献IRussell)的比较Fig.4-5ComparisonofequilibriumsolubilitydatacalculatedwiththeestimationmodelinthisstudyandreportedbvRussellat40。C250200d∞150宅;100500O50100150200250300N赴0tglL图4.670"C平衡客解度计算结果与文献(Russell)的比较Fig.4-6ComparisonofequilibriumsolubilitydatacalculatedwiththeestimationmodelinthisstudyandrepottedbyRussellat70"6中南人学硕十学何论文第四蕈Na20.AI:O,-H:O系热力学分忻300250200芒呻苦150o《100500050100150Na?0.g/L200250300图4-7100℃平衡溶解度计算结果与丈献(Russell)的比较Fig.4.7ComparisonofequihbriumsolubilitydatacalculatedwiththeestimationmodelinthisstudyandreportedbyRussellat100"C由图4-5~4.7可以看出,数据的吻合情况较好。溶液离子强度为O.5~11mol·kg“(已做单位换算)左右,计算数据仍然可以很好的与文献值吻合,这说明估箅模型对于低浓度和高浓度铝酸钠溶液都有很好的适用性,完全突破了Debye—Hfickel原始公式对溶液的限制。这也充分证明了该计算模块的实用性和准确性,并为相关研究提供了新的参考。4.3本章小结本章以三水铝石的溶出为例对Na20.A1203.H20体系进行了热力学分析,包括:1、比较了考虑活度和不考虑活童时三水铝石溶出反应的吉伟斯自由能值。结果表明,活度的引入将使溶出反应的吉却斯自由能值降低;2、考查了Na20k和毗对其溶出过程热力学的影响。结果表明,当其它条件不变时,苛碱浓度的增加或提高咏值都将使反应的吉布斯自由能值降低,即有利于溶出反应的进行;3、依靠热力学计算模块——‘‘铝酸钠溶液体系平衡溶解度及活度系数的估算”功能对三水铝石30~100"(2的平衡溶解度数据进行了估算,并与文献值进行了比较,结果基本吻合。这进一步说明了该计算模块的实用性和准确性。中南人学硕十学何论文第百蕈Na20-CaO—A120t-Si02·H20系热力学分析第五章Na20.CaO—A1203.Si02一H20系热力学分析5.INa20—A1203一Si02.H20系热力学分析5.1.1A1203·2Si02"2H20溶出反应热力学分析A1203·2S102·2H20在铝土矿中存在多种形冬,其体系热力学数握温度范围及相变如下:A1203+2Si02’2H20:S(高岭石):298-861K:S(逾丌石):298~700K:S(多水高岭土):298—700K:其与铝酸钠母液发尘如下反应:AI_,03’2S102。2H20+6NaOH+H20--*2NaAI(OH)4+2Na2Si02(OH)2(5-1)哪铷咖哪一三.一瑚脚Ⅲ㈣250300350400T/K450500550600图5-1高岭石溶出反应的△G_T图(Na201,=2309/L,ak-2,3)Fig.5-lAG·Tfigureofthedigestionofkaolinite(Na20k=2309/L.0.K=2.3)ⅡⅡ5-1可知,随温度的升高,反应的吉白斯自由能减小,但相较而言。,其反应的吉白斯自由能值的大小为:高岭石>迪开石>多水高岭土,说明这三种分子式相同的物质在此温度条件下和碱液反应能力大小顺序为:多水高岭土>迪开石>高岭石。.高岭石是Si02在铝土矿中存在的主要形态。现以高岭石为例,考查Na20k和aK对其溶出反应吉布斯自由能值大小的影响,结果如图5·2、5·3所示:中南人学硕十学何论文第百章Na20·CaO-A1201-si02-H20系热力学分析图5--2Fig.5-2aK对高岭石瘩出反应的△G-T的影响(Na20k=2309/L)onEffectofaKAG—Tofthedigestionofkaolinite(Na:Ok=2309/L)图5-3Fig.5—3Na20k对高岭石容出反应AG·T的影响(aK=2.3)onEffectofNa209AG·Tofthedigestionofkaolinite(ctK--2.3)Ⅱ口5-2、5-3可知,随着Na20k浓度和aK的升高,反应的吉向斯自由能减小,即增加苛碱浓度和aK均有利于高岭石溶出反应的进行。5:1.2钠硅渣生成反应热力学分析以碱法生产氧化铝关键在于铝硅分离,在氧化铝生产中硅主要是以钠硅渣(DSP)的形式除去的。生成钠硅渣的反应式如下:2NaAI(OH)4+1.7Na2Si02(OHh。Na20‘A1203’1.7S102‘2.3H20+3.4NaOH+1.7H20(5-2)中南大学硕十学侍论文第五章Na20.Ca(9.A1201.Si02.H20系热力学分折现考查在不同情况下钠硅渣,土成反应的热力学性质。1.碱浓度(Na20k)及苛性比(nK)恒定,二氧化硅浓度(Si02)变化取Na20k=2309/L,aK=2.3,考察当SiO:在O.1~39CL之f’日J变化时,反应(5-2)的吉和斯自由能在T=300~570K之问的变化情况。结果如图5-4所示:4o_:叫薹一』喝司啦Ⅲ啦280330380430480530580T/K图54不同SIO,.浓虚下反应(5-2)的△G·T图(Na:Ok22309'L.O,K32.3)Fig.5.4AG-Tcurvesofthereaction(5·2)atdifferentconcentrationofSiO:(Na:Ok=230∥L.QK=2-3)从图5.4可以看出.Na20k及aK一定时,Si02浓度的愈高,反应(5-2)的吉角斯自由能愈负,即更有利于钠硅渣的生成。2.苛性比(ⅡK)及二氧化硅浓度(Si02)恒定,碱浓度(Na20k)变化取U.K=2.3,Si02=1.5∥L,考察当Na20k在150~3209/L之『Bj变化时,反应(5—2)的星≥t司圈5-¥不同碱农度下反应(5-2)AG-T图(aK----'2.3,SiOr=1.59,L)Fig·5-S△G-Tcurvesofthereaction(5·2)atdifferentconcentrationofNa20R(CtK----2.3,sin2=1.59,L)吉卸斯自由能在T=300~570K之l’日J的变化情况。结果如图5—5所示:±塑盔兰堡±兰旦堡壅兰互兰堕丝旦:壁唑垒兰翌型型生墨垫查:兰!:里从图5-5可以看出,Si02浓度及呶一定时,Na20k浓度的愈高,反应(5-2)的吉南斯自由能愈负,但影响不显著。3.碱浓度(Na20k)及二氧化硅浓度(Si02)恒定,苛性比(aK)变化取Na20k=2309几,sio:=1.59/L,考察当IIK在1.5~3,5之『Bj变化时,反应(5-2)的吉却斯自由能在T=300~570K之『日J的变化情;兄。结果如图5-6所示:o巧E哪≥』dq啪哪嗡280330380430480530580图蹦不同0tK下反应(5-2)的AG-T图(Na20k=2309/L,SiOz=1,5∥L)Fig.5-6AG-TCUIⅣe¥ofthereaction(5-2)atdifferentIlK{Na20,2230e,/L·Si眈=1.Sg/L)从图5-6可以看出,NazOk和si02的浓度一定时,aK愈高,反应(5-2)的吉南斯自由能愈正,即愈难生成钠娃碴。4,生成钠硅渣时Si02的乎衡浓度在铝酸钠溶液中,碱浓度、氧化铝浓度、二氧化硅浓度、温度、品种及杂质等均会影响钠硅渣的生成速度和组成,继而影响钠硅渣的稳定性和二氧化硅在铝酸钠溶液中的平衡浓度。尽管国内外对形成钠硅渣后二氧化硅的平衡浓度避行了很多的实验研究,并提出了较多的二氧化硅平衡浓度经验公式,但由于实验者采用的研究条件不同,得出的实验结果相差较大。从理论上计算或研究形成钠硅渣后二氧化硅的平衡浓度至今未见报道,主要原因是难以制备具有确定化学式【或分子式)的钠硅座,进而无法在观有的实验技术条件下铡得各种钠硅渣的热力学数据。基于此,利用“氧化铝生产热力学数据库软件”中提出的“组成相似的一类化合物其组成与热力学数据存在良好的线性关系”来计算钠硅渣的热力学数据,以研究形成钠硅渣后铝酸钠溶液中二氧化硅的平衡浓度。据钠硅渣(设为Na20·A120y2Si02·2H20)热力学数据的估算值,结合生成钠硅渣的反应方程式,可计算生成钠硅揸时二氧化硅的平衡浓度,即2Na++2AI(0H)4’4-2Si02(OHk2r_Na20‘A1203·2S102’2H20+40H’+2H20(5-3)中寓人学面十中付论文第百章Na20.CaO.AI:O,.Si02-H20系热力学分析Gibbs自由能值为:△Gr。=AG‰一(T-298)△s瑚0+gssAC;万一丁E。垒争7'(5—4)在平衡时,△G≯=一RTinK(5-5)根掘苛性比的定义有:aK一1=x【OH‘】/x【Al(OH)41(5-6)则x[Si02(0H)22.】=(QK—I)x[OH‘]exp(AGro/2RT)/x[Na+】(5—7)掘式(5.7)可计算二氧化硅的平衡浓度,如图5.7。图s-7在不同温度下生成纳硅查时二氯化硅的平衡浓度Fig,5.7TheequilibriumconcentrationofSiO:inthesolutionofsodiumhydratealumino-silicateatdifferenttemperature图5.7说明了对于具有固定结构的钠硅渣,其二氧化硅平衡浓度是随着温度的升高而增大。据此计算出的二氧化硅的平衡浓度与文献的实验测定值【59删进行对比得知,在不同条件下生成钠硅渣时二氧化硅的平衡浓度的基本规律一致,计算结果与文献测定结果基本吻合。这从另一个方面也说明了钠硅渣热力学数据计算的可信。5.2Na20—A1203.CaO.Si02一H:O系热力学分析5.2.1水合硅酸钙在铝酸钠溶液中反应行为的热力学分析水合硅酸钙(以C.S.H表示,下同)是碱法氧化铝生产中的一种重要的含硅物质,在烧结法生产氧化铝过程中,二次反应导致铝和碱的损失主要是由于熟料浸出时水合硅41中南人学硕十中伊论文第五章Na20一CaO-A1201.SiOz—H20系填力学分析酸钙和铝酸钠溶液问的反应造成的,但其与铝酸钠溶液的反应行为不甚清楚。因此,研究水合娃酸钙翔铝酸纳溶液『日J的反应行为具有十分重要的理论和使用价值。C.S,H与铝酸钠溶液问通常被认为发生如下反应:mCaO·nSiOz·xHzO+mC032.+mH20=nH2Si042"+mCaC03+2(m-nlOH‘(5.g)(5-9)mCaO。nSi02’xH20+2nOH。mCa(OH)2+nH2Si042叫x-m)H20利用数掘庠分别计算了反应式(5-g)、(5-9)的反应吉伟斯自由能随反应温度的变化,如图孓8及图5-9所示:图5-8水合硅酸钙与Na2C03反应的△Go-T图Fjg.5.3AGO.TfigureofreactionofhydratecalciumsilicatewithNa2C03图5-9水合硅酸钙与NaOH反应的△Go-T图F噜.S-9AGo.TfigureofreactionofhydratecalciumsilicatewithNaOH中南人学硕十学伊论文第百蕈Na20—CaO—A1201-Si02-H20系孰力学分忻从图5.8可知,在热力学上6CaO·6S102.H20、CaO·Si02·H20、2CaO·Si02·1.17H20、2CaO·Si02*I.6H20和CaO·2S102-2H20均能按式(5-8)与Na2C03发生反应。并随着温度的升高水合砬酸钙越易分解:图5—8与图5-9相比,在同等条件下,Na2C03更易分解水合硅酸钙:从热力学束看,钙硅比(窜尔比:CaO/Si02)高的水合硅酸钙比钙磕比低的更易分解,其反应稳定性顺『芋大体为:CaO·Si02-H20>2CaO·Si02-1.17H20>2CaO'8i02·1.6H20。5.2.2由Ca(OH)2生成水化石榴石的热力学分析Ca(OH)2形成水化石溜石的反应方程式为:3C“OH)2+2AI(OH)4。+xSi02(OH)2。—3CaO·A1203‘xSi02·(6·2x)H20十2(工+1)OH’+2xH20式中x为Si02的饱和系数。x=0.1~2。饭鳓肫鲋舶醣沩:AG=&GT.+罴lg麦【5-1O)在式(5-10)中,当Si02饱和系数X的值不同时,参与反应的Si02的量也会不同。这样计算出柬的热力学数据孰无法进行相关比较和分析。为此,必须确定其反应的统一标准。这罨,统一设Si02的量均为“l”。以此为基准,讨论生成不同饱和系数(x=O.5,1.1.5.2)水化石榴石的反应热力学性质(kJ.molo):6Ca(OHh+4AI(OH)4.+Si02(OH)2。—·2(3CAO·A1203-0.5S102‘5H20)+60H。+2H:O(5—11)3Ca(OH)2+2AI(OH)4‘+Si02(OH)2‘。一3CaO‘A1203·Si02·4H20+40H’+2H20(5-t2)2Ca(OH)2+4/3AI(OH),"+Si02(OHF’一2/3(3CAO·A1203·1.5S102·3H20)+10/30H’+2H20(5-13)1.5Ca(OH)2+AI(OH)4‘+Si02(OHh。一O.5(3Ca0。A1203·2S102。2H20)+30H.+2H20现考察在不同情况下反应钙硅渣生成反应的热力学性质。(5-14)1.碱浓度(Na20k)及苛性比(ctK)恒定,二氧化硅浓度(Si02)变化取Na20k一2309/L,CtK=2.2,考察当Si02在O.1~3蛆之11日】变化时,反应(5-“)~(5-14)的吉布斯自由能在T=300~570K之『日J的变化情况。结果如图5.10(a~b)和图5.1I(a~b)所示:柚2020—020≥昌-20一10毒一20宅昌’10—60—60280330380130T/K'180530580.802803303804301fX480530580(a)2020(b)00一i。二0i。20主三乏‘10-oO昌+40—60.30280330380130T/K180530580—80280330380啪啪啪(c)‰∽(Si02(g/L):a-O.1.b-O.5,c-2.5,d-3.0)same图5-10相同Si02浓度下不同x值反应(5-10)的△G·T图(Na20k--2209/L.OtK=2.3Jofthereaction(5·lO)attheSi02concentrationanddifferentvalueofXrig.s.10△G.Tcurves

2023年8月3日发(作者:)

中南天学硕十学位论文摘要摘要旧版《氧化铝生产热力学数据库》在过去的使用中表现出了强大的生龛力,同时也暴露出了诸多问题。论文针对此,进行了全面的软件优化和版本升级工作,并提出了全新的铝酸钠溶液活度系数估算膜型,实现了数据库对NaAI{OHh-NaOH—H20体系活度系数及其它相关热力学性质的计算。同时,依靠数据库,论文对氧化铝溶出过程中的主要体系进行了热力学分析。新版本数据库按照数据库软件开发的要求,在用户界面的精简性、热力学数据的完整性、计算结果的准确性、反应方程式的配平等方面进行了完善,并对相应的程序代码进行了详细的修改及性能优化,彻底消除了旧版数据库使用中诸如使用不便、配平程序失效、计算失败等同题。此外,为方便用户对数据库的操作及优化升级,特制作了数据库帮助系统,并编写了《氧化铝生产热力学数据库用户手册》和《氧化铝生产热力学数据库程序说明文档》。着眼于氧化铝生产过程中最为重要的NaAI(OH)4一NaOH.H20体系,建立了基于Debye.H0ckel理论的铝酸钠溶液体系活度系数估算馍型,包括:NaAI(OH)。活度系数计算模型以及铝酸钠溶液中Na2COs、Na!S04和Na2H2si04活度系数计算摸型。NaAI(OH)4活度系数计算模型在t=25100。C,aK=1.5~10.5,mN抑HK=O.1~lOmol/kg范围内适用,成功解决平衡溶解度数据在苛性比和苛碱浓度范围上的限制。同时,还对25"C~100"C范囝内Na20.A1203.H20系中三水铝石平衡溶解度估算模型进行了优化,拓宽了其适用范围,提高了其准确度。为检验数据库的准确性及应用价值,通过数据库,论文对拜耳法生产氧l中南人’半’硕十学侍论文摘要化铝溶出过程中的Na20一A1203一H20、Na20一A|203一Si02一H20、CaO.A1203.Si02.H20等体系的主要反应进行了热力学计算和分析,得出了以下结论:(1)高岭石系列物质溶出反应的优先序列为:多水高岭土>迪开石>高岭石:(2)提高温度、Si02和Na20k浓度,均有利于钠硅渣的生成;(3)低饱和系数的水化石馏石容易生成:(4)提高SiO:浓度f『利于水化石榴石的生成;提高温度、Na20k浓度和CtK愈高,均不利丁J水化石榴石的生成;(5)在稳定性方面,低温下水化石榴石稳定,高温下钠硅渣稳定。结果证明,数据库计算结果准确,热力学趋势与工艺中大部分实验事实相吻合。这从一个方面证实了数据库的可靠性和实用性。此次开发优化工作以计算机技术和氧化铝专业知识理论为保障,全面实现了数据库的易用性、准确性、适用性和安全性,确保了新版本数据库合理的生命剧期以及良好的应用前景。关键词氧化铝,热力学数据库,优化,应用n竺苎望!!!竺!塑!!竺呈!!生兰望垒里苎!坠垒篁!TheOptimizationandApplicationofThermodynamicDatabaseforAluminaProductionABSTRACTDuringtheuseofThermod)’namicDatabaseforAluminaProductionoldversion,manyproblemswerefoundinspiteofitsgoodperformance.Inthispaperthesoftwareoptimizationandeditionupgradewerecarriedouttosolvethoseproblems.AndthenewestimationmodelforactivitycoefficientsinNaAI(OH)4-NaOH-HzOsoftwaretOsystemwasestablished,whichwasusedinthedatabasedealwiththerelevantthermodynamiccalculationsimmediately.tOAccordingthedemandofthesoftwareexploration,thecomplementarityuserofthesoftwarewasmadefromthesimplenessofinterface,theintegrityofthermodynamicdata.theaccuracyofcalculationaIresultsandtheequilibriumofreactions.AndalIoldrelevantcodeswerecompletelyamendedtooptimizethewholedatabasecapability.Maindisadvantagesoftheolddatabaseversionwereasovercomedduringthefurtherexploration,suchtheinconvenienceofusing,theinvalidationofequilibriumandthefailureofcalculation,Inaddition,thedatabasehelpsystemwasestablishedand“TheUserManualofThermodynamicDatabaseforAluminaProduction”and“TheProgramInstructionsofThermodynamicDatabaseforAluminaProduction”werewrittenforsupplyingconveniencefortheuserstOoperateandupgradethedarabase.aFocusingontheNaAI(OH)4-NaoH—HzOsystemwhichplaysparamountroleintheindustryofaluminaproduction,acompleteestimationmodelforactivitycoefficientsinNaAI(OH)4一NaOH-H20systembasedontheDebye—Hiickeltheorywasproposedinthispaper.Itincludes:theestimationmodelsforactivtitycoefhcientsofNaAI(OH)4、Na2C03、Na2S04andNa2H2Si04inNaAI(OH)4-NaOH-H20system.Themodelcanbeweltappliedintherangeof25"C~100"(2,q=I.5~10.5,mNnOHK=0.1~10mol/kgtherangelimitsofthevaluesofO.KandsuccessfullysolvedandmNaoHgcalculatedfromsolubilitydata.111璺竺!竺堡!!垒墅!!坠£!Simultaneously,theestimationmodelsforequilibriumsolubilityofaluminaintheNa20.A1203一H20systemwenoptimized.Anditsapplicationandaccuracywerewidenedandimproved.InordertOcheckaccuracyandtheapplicationofthedatabase,themainreactionsintheBayerprocessofaluminaproductionsuchasNa20一A1203一H20、Na20-A1203-Si02一H20andCaO-AIz03-Si02一H20systemwerecomputedandanalysedinthispaper.Andsomeconclusionsweregainedasfollowing:(1)Thepreferentialsequenceofthereactionwithalkaliis:halloysite>dickite>kaolinte;(2)IncreasingthetemperatureandtheconcentrationofSi02andNa20kwillacceleratetheformationofthesodiumhydratehydratewithlowersaturationalumina—silicate;(3)Thegarneteasier;(4)Highconcentration0tKcoefficientwillformofSiozandlowtemperatureandlOWconcentrationofNa20kandfavorablefortheformationofthegarnethydrate;(5)AsgarnethydrateisissteadiersteadieratattOwillbethestability,thesodiumhydratelowtemperatureandthealumina-silicatehightemperature,Theresultsshowthattheareaccuratecalculationresultsfromthedatabaseareandthethermodynamicgendsinaccordancewiththeexperimentalfacts.ThishighlyprovedthereliabilityandpracticabilityoftheThermodynamicDatabaseforAluminaProduction.Thefacility,veracity,applicabilityandsecurityofthedatabasesoftwareweregreatlydevelopedbasedonthecomputetechnologyandaluminatheory.ensuresThenewversionofthesoftware6lmre.thesarisfy.inglifecycleandapplicableKEYWORDS:Alumina,thermodynamicdatabase,optimization,application原创性声明本人声明,所呈交的学位论文是本人在导师指导下进行的研究工作及取得的研究成果。尽我所知,除了论文中特别加以标注和致谢的地方外,论文中不包含其他人已经发表或撰写过的研究成果,也不包含为获得中南大学或其他单位的学位或证书而使用过的材料。与我共同工作的同志对本研究所作的员献均己在论文中作了明确说明。作者签名:Z同期:年月同关于学位论文使用授权说明。本人了解中南大学有关保留、使用学位论文的规定,即:学校有权保留学位论文,允许学位论文被查问和借阅,学校可以公白学位论文的全部或部分内容,可以采用复印、缩印或其它手段保存学位论文:学校可根掘国家或湖南省有关部门规定送交学位论文。艏戤:删刷性g.嗍:年月同±塑叁主竺±芏竺鲨塞一苎二垦奎苎笪垄第一章文献综述1.1热力学概述㈦热力学是关于热、助、能以及在系统状念中由它们产生的变化的研究.是一门研究能量转换以及与转换有关的物质性质(或参数)问相互关系的科学,从广义上讲,热力学是研究系统宏观性喷之『日J的关系。热力学中一个关键的性质是温度,有时热力学也破定义为温度与物质宏观往质之间的关系的bK究。热力学足稳学的宏观理沦.是坂能量转化的观点研宄物质的热胜质。阐明能量从一种形式转换为另一种形式时应遵循的宏观规律,热力学是根据实验结果综合螯理而成的系统理论。它不涉及扬磺的微疆结构和微观粒子的相互作用,也不涉及特殊物质的具体性质,是一种唯象的宏观理沦,具有高度的可嘉眭和普遍睦。热力学的完整理沦体系是由几个基本定律以及相I立的基本状念函数构成的,这些基本定律是以大量实验事实为根据建立起束的。无论多少个物体互帽接牲都能达到热平衡.并且如果A物体同时与B、C两物体处于平衡念,则B、C两物体接触时也一定处于平爱念丽不发生新的变化,这一热平衡规律弥为热力学第零定律。出此可以引入一个状态函故,品宦。温熏是判定一系统是否与其他系统互为热平衡的杯志。热力学第一定律就是能量寸恒定律。是后善在一切:步及热观象的宏哩过程中妁具体友观。描述襄统热运动能量的状念函数是内能。通过作J力、传手丸.系统与外界交换能奄,内能改变。热力学第二定律指出一切涉及热髓象的宏观过程是不可逆的。它阐明了在这些过程中能量转换或传递的方向、条件和限度。相应的念函数是熵.熵的变化指明了热力学过程进行的方向,熵的大小反映了系统所处状念的稳定性。热力学第三定律通常表述为绝对零鏖时,所有纯物质的完美晶体的熵值为零。它意味着。当温度趋近于OK时,所有;疑聚物系的墒值趋近于一个相同的极限值。1.2热力学计算模型物质基本热力学性质的基准叁约定如下【¨51:(1)焓:对纯物质束I兑,规定标准状态下(298K,l大气压下)稳定单质的生成焓H。298为零;对离子来说,规定H+的生成焓拶298.H-为零。(2)熵:在绝对零度时,完整晶体的纯物质,其熵值规定为零:对离子来说,规定任中南人学硕十学付论文第一章文献综述何温度下H+的标准熵为零。1.2l纯物质热力学性质的计算模型【5·6‘7】绝物质热力学性质的计算主要包括:计算任一温嚏T(K)下的热容(Cp)、焓(H-r)、熵(ST)和吉南斯自出能(GT)。计算表达式分别为:Cp=a+b×10‘3T+c×106T-2+d×lO’叶2或Cp=a+bxl013T+c×lO‘T2其中a、b,C,d为热容系故H-r=H。29s+『c,。dT+JCp.dT+.’+』c,dr+∑叫。(1·2)ST=S。298+狰刀+净T仉…+辟JT仉莩半GT=HT·TxST(1-3)(1—4)式中的下标0,l,2…。rl表示相念,AHT.表示相变热。1.2.2水溶液中离子热力学性质的计算模型f%1,1)计尊水溶掖中离子在任一温度T(K)下的焓(H。)、墒(Si)和吉南斯自由能(G.)。计算表达式分别勾:r△H.=△H。299十解刀29fl,‘(1—5)/kS,=ZXS。29s+:。t争△G.=△H。-Tx△S。小,可以忽略不计。故其均未考虑压力的影响。(1·6)(1·7)注:以上三式同样适合纯物质相关热力学性质的计算。在300"C以下压力的影响很(!)计算水溶液中离子的热容(C。)值由于离子的热容值还很缺乏,一股采用Criss提出的“离子熵相对原理”‘42“■束进行近似计算。其汁算方法如下:Cp。=BxTxl04(1.8)在这里,Cpo为离子在温度298和T之闻的平均热容,B是热銮系数。从式(1.8)可以看出,离子的平均热容与温度成线性关系。B54.184×(a+bxS*29&4bs】)(1-9)中奄人学母!十中伊论文第一摩文献综述式中,a.b是仅同离子类型有关的常数。表l一1列出了不同炎型离子的a.b值。S。29sIab。是离子在298K下的标准绝对墒。通过(4—10)式可以计算得出。S。298lahs}2S。208—2092、Z(1.10)其中,Z为离子电荷数(带符号),·20,92(J—mol-*.K。)是H+在298K下的绝对熵。S。!。B为离子的相对墒。可从数据库中查到。表I-t各美耋}子的a,b值嘲Table1-lvaluesofaandbofvariousiOnS离子类掣打~叫a竹…由,简单田f离f简单胡离f及OH‘禽氯珞胡离f’’‘…“~禽氢氧珞月子………。。i页再F—1砭孺r——_j丽矿————i而=广一(XO”。)(XO。(OH)”1)1.65x10。。,I.12x10。5.83x10。111.2x10‘4bI3)H‘热力学性质的计算H+不楫f四种类型的离f之一.其热容与,甚度之问的关系为COpH‘=0.3912T(I一11)从而可以得到如下关系式:H。TH+=0.5x0.3912。(T2-298,152)(I—12)f1.131(1.14)s。rH+(absl=S。298HjabsI+o.3912(%298.15)G。TH+=H。TH+一TxS。r14i曲s1.2.3水溶液中分子热力学性质的计算模型一股情;兄F,弱酸弱碱在水溶液中不能完全离解,呈中性分子的状叁。水溶液中分子的媳力学计算通常采用Helgesonl81给出的近似公式.浚公式主要适用于段低的温度(200℃以下)。假设水溶液中一未离解的中性分子A。B。.其离解方程式为AmBn专mAY++nBz.(1-15)式中m与11一电离时产生的阳离子数和阴离子数Y+与z-——阳离子电荷数和阴离子电荷数离解反应(1一15)的标准反应自由能为:4△G。T(d。s}=一△S629815(d6,[298.15--二fl--exp(exp(b—aT)--c+(T-298.15)/e)}】+△H629815(dls}珊(1·16)式中△S。298.15(dis)一离解反应(1一15)在298.15K时的墒变生童叁竺堡±兰生堡皇AH。298l蛐。)——离解反应(1一15)在298.15K时的焓变.坚一奎苎笙堕a=0.0187:b=.12.74l:c=exp(b+298.15a):0=219:∞5l+ac0离解反应(1—15)的墒变为:△S。T㈨J-一AS。29815m1{exp[exp(b+aT)一e+(T-298.15)/0]}[1+qoexp(b+aT)I/∞其中0=4.106(1·17)离解反应(I一15)的焓变为:AHoT(dⅢ=AO。T【dⅢ一TAS。T(dIs)(1一18)另一方面。离解反应(卜15)在任一温度下的反应吉钿斯自由能也可以写成:△G。Tldtsl=从而mGoTIA,+nO。TtBl—G。T【AmBnl(1·19)(1-20)(1-21)G。TiAmBnl=mG。TIA)+nG。nBl—AG。Tlald¨同理:HoT【AmBn)=mH。T‘A1+nHoT(Bl-△H。T(dIs)在式(1.20)(1.21)中A,B均为离子,其高温下的热力学性质可通过离子相对熵原理求出。1.2.4化学反应的热力学计算模型【35·79I对化学反应来晚。任一温度T(K)下反应焓(△HT)、反应熵(△sT)、标准反应吉白斯自由能(AGr}以及平衡常数K的计算可通过下列公式进行:AH。r=∑bⅣI)一∑0ⅣIJ’I啦拐(I一22)址’硼△s。t=∑b,s-)’}墟锡一∑h黑j(1-23)衄J4罚AO。t=∑b60)一∑bG0)(1·24)厘pt研1成研△GoT=.RTLnK即LnK=.AG01嬲(1-25)式中,H。i’T、s。i,T、G。i,T分别表示i物质的焓、熵及标准吉南斯自由能值;vi表示化学反应计量数:K表示平衡常数。必须指出的是,现有的热力学数据手册,有的列出的是物质的生成自由能△G。f=△H。r—T△S。f有的列出的是章尔吉如斯自由能GoT=H。T--TS。T无论使用哪一种,只要不交叉使用,计算出来的同一反应的吉白斯自由能是一样的·在4中南大中硕十中何论文第一章文献综述本论文中,使用的是生成吉如斯自由能和尘成焓。在以上的讨论中,△GoT代表的是化学反应的标准反应自由能。杯准的含义指的是每一种参加反应的气相物质分_lI部处在】大气压下,参加反应的溶r瘦浓度为l(mol·L“)。而实际7芷产中的情况往往比标准状念条件复杂得多,所以,需要讨论在任意指定的状念下,体系能否发生化学反应。这就要用到下面的等;盛方.智式(1--26)束作为判断化学反应方向的补充。假定反应为aA+bB=IL+mM△GT=△G。T+RTInJ=一RTInK+RTlnJ(1-26)舯,j-籍或籍虹绷混色PL,PⅥ,PA,PB,aL,a、I.aA,aa分别为反应物和生成物的指定分压和反应活变。从(I.26)可以看出,当J<P时.△GT<0,反应能够自发地向右进行,反之.则不能。严格.兑,判断反宦方向应以(1.26)为依掘,但因为J在对数二女内,所以△GoT对△GT的影叫更大。也就是:蔸.当△GoT的绝对值很大时.△GT的IF负号由△GoT束决定。一殴地l并.△(jor<一4I.84f焦时.反应不仅能在标准‘帙仑F进行.在一收惰:兄下也能进,-:△G。T>41.84干焦时.反应则不能自发进行㈣。1.2.5反应的AGoT—T二项式的回归模型在实际计算中,AG。T与T的多项式用起柬比较寐颐,人们喜欢使用G。T的二项式:△G。r=A+BxT。通常可将化学反应的△G。T=f(T)多项式整理成二项式或根据一系列比较准确的/XG。T值利用回归分折法(即最小二乘去)求得AGoT的二项式。按指定温度计算的值越多,用回归分折法计算出柬的AG。T二项式就越准确。最小二乘法的基本原理如下:实验值G,、T,(i-l,2t…,n)误差6i=GI—A—B×Ti目杯函数J_62.=(G.--A--BxT.)2根据极值原理,分别对A、B求偏微商并令其为0,从而得出A、B值:A=否一B×于(1·28)(1·29)B=∑(7=一于)(q一否)/Err,一_2中南人中硕+一中何论艾第一章文献综述瓦甲于2吉喜7:否2吉萋GI相关系数,:圣!三二!!呈三童!1.2.6复杂化合物热力学数据的估算方法1.1.6.1组成相似的无饥化合物热力学陛质的线性估算模型以水合铝酸钙为例加以说明:水合铝酸钙主要存在的分子式形式有:3CaO·A1203·6H20,4CaO·A1203·13H20,4CaO·A1203-19H20.CaO’A1203‘10H20以及2CaO·A1203"8H20等。从氧化物裹示的组成可以看出CaO和H20的分子比在变化,A1203的分子比始终为l,可设汁其,上成焓线性方程式如下:△Ho=△Ho。+m(CaO)△HI。十n(H20)△H二。(1-30)其中△Ho为各水合铝酸钙的生成焓:△Ho。、△HI3、△H2。是和生成焓有关的参数:111、n分别是复杂化合物以氧化物襄示的CaO和H:O的分子比。在标准状念下,各种矸三念的水合铝酸钙的生成焓均可以在数据手册中查副,以多元线眭最小二乘法进行拟合,得到如下结果:△Ho(1d·mol“)=一1633.23.704.2m(CaO)一29498n(H20)同理,可以对水合铝酸的自由能、熵和热容进行拟合,得到的结果如下:AGor(kJ·molll=.1574.76-672.88m(CaO)-237.08n(H20)So(1d·mol“·K“)=57.0l+39.14m(CaO)+38.92n(H20)组成相似的化合物在相同温度下的恒压热容之问也存在一定的线眭关系。即:Cp(T)=Cpn+m(CaO)Cpl+n(H20)Cp2(1-31)通过上式。可以求出复杂无机化合物在不同温叟下的热奋值,然后再利用线性回归求出热容系数,从而可求出其它温度下的热力学性质。1.1.6.2复杂无机化合物(至少包含两种令属离子)的拆分估算模型根据温元凯算法的主体思想,把任何复杂硅酸盐矿物看作各种酸性氧化物和碱性氧化物的复合物。复合时碱性最强的氧化物首先和酸性最强的氧化物反应,然后再和酸性次强的氧化物化合【lo¨l。水在大部分情况下以水合£F的形式存在,有时也作为结构水参与反应。这样,复杂硅酸盐的吉钿斯自由能(或焓)可看成由各氧化物的吉白斯自由能(或焓)之和(EGoOxideorEHoOxide)与氧化物间的反应自由能(或反应焓)之和(EGoRorZHoR)6中南人学硕十学伊论丈第一章文献综述两部分组成。1.3建立氧化铝生产热力学数据库的目的和意义目前,我固的氧化铝尘产还没有一套完整的仿真与自控体系。计算机技术在氧化铝生产中的应用同国际上相比,还存在着很大的差距,特别是蘑f汁算机模拟与控制的软件和数掘库1l常缺乏.已成为发喂爱化铝生产自动化的}要贿碍。今天,任何计算机过程摸拟与自动控制的应啊软件均是建立在大型的热力学数据库的基础之上的。因此,要提高我国计算机在氧化铝生产中的应用水平.加强工业生产的自动化。除了吸收引用国外的先进技术柬结合我固的生产实际外,加强我囤自己的软件和数扼库丌发能力,大力丌展数据库、控制软件等整础性的研究工作十分迫切。我【川的氯化铝生产具有很强的“中因特色“,这是由我固的铝土矿资源决定的,其特点是:(1)铝上矿资源#富,遍及19个省,居世界第四位;<2)铝t矿矿石类型以一水硬铝石为主,占全国总诸量的99%,三水铝石型铝上矿仅约占l%;(3)绝大部分矿石品位冒中等,铝硅比多分1}】在4 ̄10之1日J.约占全国总锗星的85.6I%:(4)矿石具有高铝、高硅、低铁、中等};吕硅比,难藩的特点,且伴生矿多.成份复杂。铝土矿矿右_的这些{寺点决定了我雷的氯化锅r业主要采用混联法和烧结法.目前.除了乎果铝厂采用纯玎尔法外,其它五个氧化铝厂均采喇烧结法或联合法。同玎尔法相比,烧结法工艺流程长、连续性差、能耗高、产品质量差,成本高。我国氧化铝工业的发甏严重滞岳于电解铝工业,每年复化铝缺口部比较大【‘2-’314”】。如何提高产品质量,进一步缩小我国与世界发达国家在氧化铝产能、成本等方面的差距,提高我国氧化铝生产在国际上的竞争力,己成为我国蓑化铝生产厂家主要的攻关方向,因而必须采用新的生产技术去改善(或完善)我固的氧化铝生产工艺,这就需要对氧化铝,土产过程的基础理论有深一步的认识。备矿物及其在体系间的相互反应热力学是氧化铝生产基础理论研究的一个非常重要的方面,其可为新工艺的提出、实际生产控制等提供理论支持。然而,氧化铝生产中的热力学研究目前存在蓄许多困难:l、出于铝土矿及其它含铝的矿物资源如霞石矿、明矾石矿以及各种铝硅酸特岩石等均是组成复杂、化学成分变化很大的矿石,因而在氧化铝,上产中存在着物相多、物理化学行为极其复杂等特点。特别是所涉及的物质大多为二元、三元或三元以上的复氧化物,这一类物质可供参考的热力学数据菲常缺乏;2、现有的热力学数据分散、零乱,没有公认、统~的数据标准:3、计算复杂,手工计算难于完成,且准确性较低:4、缺少先进、现代化的工具进行辅助分析。7中奄人中硕卜中伊论文第一章文献综述针对以上问题,中南大学氧化铝及陶瓷研究所将热力学、计算化学以及数锯库技术的结台起束,研究和丌发氧化铝生产热力学数据库。它的建设成功将解决我国氧化铝工业界广大科技工作苦长期面}临的某些工作难题,可对氧化铝生产辐关工艺过程进行模拟i同时.目前在我国还很少有关于氧化铝生产热力学数掘库方面的报道,它的建成也将填补我国在浚方面的一项空白。I.4本论文的主要研究内容I、氧化铝生产热力学数据库的优化和完善.包括耐用户界面的精简性、热力学数掘的完整性、计算结果的准确性,反瘦方程式的配平等方面进行完善,并对相『萱的程卑代码逆行详细的修改及性能优化:2、编写数堀库源代码程序说明文档以便于今后优化和升级工作的进行,编制数掘库用白使田手册和帮助系统以方便用户使用;3、铝酸钠溶液中活度系数计算模型的完善,含铝酸钠溶液体系中NaAI(OH)。Na_-C03、Na2S04和Na2H2Si04活壹系数计算馍型的建立以及三水铝石平衡溶解度估算漠型的陶造:4、氧化铝尘产热力学数据的应用,包括对玎耳法,圭产氧化铝溶出过程中的Na:O—A1203一H20、Na20·A1203·Si02一HrO、CaO-A1203一Si02一H20等体系中的主要反应进行热力学计算和分析。中啕人学面十学伊论文第.章氧化锅J上产执力学数据的基础俨究第二章氧化铝生产热力学数据的基础研究2.1与热力学相关的基本概念1.状态函数mI(StateFunction)热力学中热力系统在某一瞬问所处的某种宏观物理状;兄称为状念。用宋规定状念的性喷叫做热力学变星或状念性质.一旦规定了状念的{丸力学,殳璧.其他平衡状念性质就是这些变蕈的函数,称为{丸力学函数或状念函数。如压力、温室和比容等。状念参数的数煎仅仅取决于热力系统的状态,丽与达到这种状念昕经历的热力过程无关。2.标准状态㈤(StandardStateS)气体:在p_101.32kPa下的理想气体液体:在p=101.32kPa下的纯液体固体:在p=101.32kPa下固体所处的最稳定的晶体状态3.热容…920I(HeatCapacity,)一定量的物喷,在不发生化学反应和相变的情;兄下。温度升高(或降低)l℃(1K)所r受收(或放出)的热量称为该物质的热容。当系统吸收或放出热量时,它的i昱度可能变化,也可能不变化,这取决于所经历的过程。在传递热量Q的过程中,系统的温室从TI变化到T二.则平均热客定义为:平均媳奋=—羔一7:一Il当Q和T2一TI部趋于零时,这个比值即为热容C:c:1im旦:丝,j-^7j一7;dT:!生z-Ij?j一7;.dT堡c:Iiml每单位质量的热容称为比热容:单位为J,(kg·K)。每单位气体体积的热容为容积热容:c,:.!一塑:风塑22.414珂d7’”d7’单位为kJ/(n·m3.K)。,00为标准状态下气体的密度如果物质的量是l摩尔,它的热容则称为摩尔热容:G=吉等=M拶dq单位为J,(m01.K)。9!塑叁之堡±竺竺堡奎墨三兰墨竺堡竺兰!堕!墼堑!!蔓壁竺茎等匝过程中的热容量成为定压热容:Cp=(等),等客过程中的热容量成为定容热容;c,=【掣),口1热容与温度的关系式一般有两种方式:Cp=a一-b×10‘3T+cxl0吁2+dxl04T2和Cp=a+b,10.jT+cx!05r2。热吝增晕-ACp:化学反应中.产物的总热容减去反壹物的总热容。4.焓‘212二J(EnthalpY)对任何系统来说,焓的定义为:At=U+Pv式中u为系统的内能,尸为系统的压力,V为系统的体积。焓是状念函数。当系统进行某一过程时。焓的变量£XH只与始叁、终奁有关,而与中jt日J过程无关。单位为J/mol或kJ7tool。杯准焓日o:25*6"时处在杯准叁的l摩尔物质的焓值。标准生成焓△H:,:当觇定每种电质在lol325Pa和25。C时的焓值为零,则由有关草质左101325Pa和25"C时生成lmol化合物所吸收的热量为改化合物的杯准,土成焓。化学反应的焓变:等于产物的总生成冶减去反应物的总生成焓:△,H。=∑v。A,H,。可以从一个温嚏的化学反应的焓变求得另一个温度的;教化学反I萱H的嗡变:A,H.。=△Ho+I△cPd7’。5.熵(Entropy):由热力学第二定俸所确定的{丸力状念参数,是系统中分子运动混乱程噎的寰萎.以S表示,单位J/K。熵是弧立系统自发改变状念的可能性的量度,改变看蛹增越大的状念,改变的可能眭越大。熵的微分定义是;系统在可逆过程中同外界交换的热量与系统热力学温受的比值豳:f墼1。Lr几堂标准熵So:25℃时处在杯准念的l摩尔物质的熵疽。杯准生成墒岱?:指在标准状态下,由处在标准奎的稳定单质在25℃时恒温生成处在标准仑的l摩尔化含物所引起的熵变,称为孛,下准,土成膊。6.吉布斯自由能[23l(GibbsFreeEnergyFunction):吉南斯自由能是状念函数。其定义表达式为G=H—TS。衷示体系总能量中具有作非体积功能力的那一部分能量。在25℃时处在丰,下准念的l肇尔物质的占伟斯自由能称为f,下准吉如斯自出能(Go)。标准生成自由能AG?:化台物的标准生成自由能是指在标准状态下,由处在标准态的稳定单质生成该化合物时自由能的变化,称为标准生成自由能。规定平衡状态下的稳定单质在任意温度下的标准生成自由能为零。7.化学反应平衡常数【24】(ChemicalReactionEquilibriumCcnstant)化学反应体系在~定的温度和压力下处于化学平衡状态时,各组分的活度不随时问lO中南人学硕十孚伊论文第一章氧化锅t产执力学数据昀鐾础俨究而改变。化学平衡常数表示处于平衡状念的子组分活度之『日J的数量关系。K越大,反应物的转化事越高。反应越完全。8.活度(Activity)及活度系数(ActivityCoefficients)活度,又称有效浓度,即单位体积电解质溶液中表现出柬的表观离子有敏浓度,即扣除了离子问相互作用的浓度。通常以口(Activity)表示。活度系数表示实际溶液与理想溶液的偏差,通常用r表示。2.2氧化铝生产中相关物质热力学数据的收集与甄别本氧化铝生产热力学数据库的基本数掘包括纯物喷热力学数扼和水溶液热力学数掘。纯物质热力学数掘主要柬自:1986年版B.H.6a6、IXIKHH等编著的{TEPMO丑IIHAMHKACH3HKATOB))1251,1991年版OKnacke编著的{ThermodynamicalPropertiesofInorganicSubstancesII》131和1993年版Barin编苫的{ThermodynamicalDataofPureSubstance))【26l。另外,以1985年版林传仙等编著的《矿物及有关化合物热力学数据手Iil}》f’71作为补充。另外,数据库中还输入了50余种复杂硅酸舱的估算数据,以供研究蕾参考。水溶液热力学数据主要综合了1982年皈DonnaldD,Wagman和WilliamH.Evcans编著的《111eNBSTablesofChemicalThermo#namicProperties,SelectedValuesforInorganicandCtandC2OrganicSubstancesinSIUnits))14j和((TEPMOjIHHAMHKACHnHKATOB》。日i昔是美国国家标准局于1965年丌始编辑的NBSTechnicalNote系列的修订和综合,可靠性较高。在甄别数掘时。原苏联《TEPMOⅡHHAMHKACH。qHKATOB))中的数据权威性较高.国内外大多采用其用于氧化宅吕生产热力学计算.故优先考虑。若某手册数掘与其它手册中的相差较大时,以多数情况为准。同时,把{ThermodynamicalDataofPureSubstance))和《矿物及有关化合物热力学数据手册》的数据作为补充,热容公式参考(《ThermodynamicalPropertiesofInorganicSubstancesII》。对于数据手册中没有热容公式的物质,数据库部进行了一一拟合,并通过热力学函数对拟合值进行返算,证明拟合结果正确后,将拟合值输入数据库中。从数据库中可查到的物质热力学性质(即字段)有:编号(NO)、分子式(formula)、相念(phase)、热容系数(a,b,c,d)、生成焓(/XHO)、生成熵(ASO)、标准焓(Ho)杯准熵(SO)、温度范围(TempUp,TempUown)、晶体转变温度(Ttr)、相变焓(Htr)、相变熵(Str)、熔化温度(Tm)、熔化焓(Hm)、熔化墒(Sm)、蒸发或升华温度(T_b)、蒸发焓(Hb)、蒸发熵(Sb)等。中南大学硕十学伊论文第一二章氧化锅生产热力学数据的蘑础研究目前,氧化铝,土产热力学数掘库已经收集了1199条记录。其中纯物璜为839条,水溶液为360条。2.3铝酸钠溶液活度系数计算模型的研究2.3,1概述氧化铝工业中铅酸钠溶液浓度普遍较高,而且由于铝酸钠溶液本身所具有一系列特殊的性质,如饱和溶液的超常稳定性,使之无法等同于其他电解喷溶液(如:纯NaOH溶液),活度系数研究的难度很大。同时,在铝酸钠溶液分解机理和溶液结构方面。学术界还存在分歧128羽】,也加大了改体系研究的难度。因此,尽管多年柬人们对铝酸钠溶腋体系热力学性质进行了大量的研究工作f30.31.3233】,获得了很多溶解寰.实验数掘。但是直接关于渗透系数和活度系数的研究非常少。1992年.Wesolowski应用Pitzer模型处理三水铝石溶锯度数厄,近似得到了一组有关NaAI(OHh的Pitzer参数和AI(OH)4一,OH‘和Na+混合的参数川。1999年.Park和Englezos采用等压法测定了NaOH.NaAI(OH)4·NaCI—H20体系从mNaOH=0.85,mmNfl=O.15(mol/kg)至0mNaOH=2.6,mNaA“OH)4=0.05。NaA]iOH}4一--0.35,mNacl=1.04(mol,7kg)在298,15K的渗透系数.并用Pitzer模型回归出Pitzer参数PSl。2002年.中南大学化学化工学院的周俊采用等压法测定了总苛胜碱从O.05mol/kg到12mol/kg,苛性比从1.64到5.53的铝酸钠溶液在313.2K的渗透系数,应用Pizter摸型回归得到NaOH.NaAI(OH)4一H20体系的Pitzer参数。其中一组313.2K时的Pitzer参数为:∥…haoH:13㈣NaoH:o.0240o.08669.13”lN。oH:O.31446,.o.00007367,c气aOH:o.003180.∥mNnl{0HH:0.03507,13“’NaAI(OHl4:O.081771,C~aAl(OHw:.O.001066,00H.Al(oH)4‘:t、|,Na+OH’AI(OH)4’:.O.01162t36l。这些实验研究对于铝酸钠溶液活宣系数的计算部有熏要的意义,特别是周俊获得的这组Pitzer参数使得准确计算40"C的NaOH.NaAI(OH)4一H20体系活度系数变得可行,所有相关结论也为以后的研究提供了极有价值的参考。但是,仅仅这些结果显然不足以支持使用Pitzer模型计算一定温度和浓度范围的铝酸钠溶液活度系数,大量其它温度的实验研究工作被期待。另外,有研究表明,现有的Pitzer模型并不适应于处理100"C以上的电解质溶液体系【371。很遗憾,由于以上原因,无法将Pitzer模型作为准备建立估算模型的基础。除直接的实验研究外,Debye.Hiickel模型在涉及含有NaAI(OH)4的溶液中使用比较常见。有人尝试用早期的扩展方程柬处理铝酸钠溶液13。j,但因参数值选取不准确以及数学模型上的缺陷,并没有给铝酸钠溶液活度系数的计算提供太大的帮助。我们注意到。中南人学硕十学伊论文第一章氧化锅生产执力学数据的垦础研究近年比较热门的Helgeson等人提出的Debye.HUckel扩展方程同样玻用于含有NaAI(OH)4的体系。例如:1996年.Diakonov等人在汁算地壳高温流体Na-AI—CI—OH中使用浚扩展方程计算了包含有NaAI(OH)t的混合溶液的平均活噎系数,忽略溶液中备离子a参数的不同,取均值为4A1391。另外,在1995年Pokrovskii和Helgenson关于A1,01.H、O.NaCI体系热力学性质的完整报告中,也旋刘应用该扩展镆型处理铝酸钠溶液。但荏随后讨论中假设OH’和A102。活噎系数相同,并没肓具体i步及活宦系数的计尊ItSl。综上,Debye—Hfickel模型在铝酸钠溶液活度系数的使用上也没有完整的解决方案。由于该模型简单,对于估算模型的建立有较强的吸引力。但是,如果要使用Debye.Hfickel模型作为估算馍型的基础,首先要给出针对NaOH.NaAI(OH)4.H20体系的Debye.H0ckel模型应用假设,同时给出使用Debye·Hilckel模型涉及的物理化学参数的具体值:再则,验证Deb、e—Hfickel模型对NaOH-NaAI(OH)4.H20体系的适用性,确定Deb、e.Hfickel模型町以应用于陔体系;完成上述工作,/j能基于此讨论具体建摸过程。2.3.2先前的计算模型及存在的不足李小破等曾基于一定范围内三水铝石在铝酸钠溶:夜中平衡溶解宦故掘,借助Debye.Hfickel模型,通过选取NaAI(OH)4合理的、#}£,引入”铝酸钠溶液的表观介电常数∥“的概念,将各种影响囚袤归结到£’,提出了铝酸钠溶液备组份活度系数的计算方:去ml。其所建立的活度系数和表观介电常数计算摸型为:’’[NaAIIOH)4=pl+p2’mh,aOHK+p3’mNaOHK‘+p4‘U.K+p5‘O.K2+p6‘lnmNaOHK+p7-lnaK(2-1)F’=pl+p2mNaOHT+p3‘mNaOHT‘十p4‘QK+p5‘aK2+p6。mNaOHT。aK+p7’mNaOHTz"QK+p8。mNaOHTaK。+p9·m~aOHTz·aK‘(2-2)通过上述公式可计算出铝酸钠溶液相关活度系数.但其中也存在一些问题:(1)根据平衡溶解度数据拟合的¨a^IIOHH与CtK、mNaoHT、T的关系模型对40"C时的覆盖性偏差较大:(2)苛性比(Ⅱk)在1.5~5.5和5.5~10.5范围内时,其,舌度系数的变化规律不一致:(3)在通过利用铝酸钠溶掖的表观介电常数占’建立的naA“oH,4计尊馍型,其方程较为繁琐,用其计算活度系数时,计算结果与实际值偏差较大:(4)由Debye.Hnckel模型可知,原s7的计算模型不能解决活度系数大于l的铝酸钠溶液的表观介电常数问题。因此,本文针对NaOH-NaAI(OH)4.H20体系,需解决:(1)探索适用于所有温度下统一、较精确的计算模型;(2)拓宽计算模型的应用范围;(3)提高计算模型计算的准确性,解决活度系数大于l不能计算的难题。!塑叁竺堕±!堡垒塞2.3.3铝酸钠溶液活度系数计算模型的建立兰=三童墼坐堡竺兰竺垄竺垫堡塑堡!|!!竺塑建立NaOH.NaAI(OHh.H20体系活度系数}十算模型豹主要思路妞图2-1所示:l平一婶餐¨卜Yk…Ⅲ..的计算l与文献Y”。..债比牧l1r平衡篆竹FYu。。计算礁掣的拟台I拟台占,与a。、m,,的函数天系Ity。。。汁劳硅掣图2-1Fig.2-1确立孙。n‘)HH估算接型曲虎程堡Flowsheetofbuildingmodelfor丫N“I.oH“2.3.3.1平衡条件下NaOH.NaAI(OH)4.H20体系活度系数计算模型根据文献报道的三水铝石溶解赛数据【4142’,拟合了在乎衡状叁下。铝酸钠溶液活宣系数与温宦、a。和m的函数关系式:丫N直A“oHw=pl+p2’肌+p3’脚二+p4’吼+p5’也2+p6’lnm+p7’lna,(2·3)式(2-3)中不同温度的调整参数值,详见文献¨引。其中40。C时的调骜参数为:pl—1.2557,p2=0.OJ7t.p320,0045,p4=一1.1159.p5=0,0377,p6=一00563-p7--4.0364。由于二氧化铝生产过程中备组份一般处在非平衡、非杯准状念下.因此需将公式(2-3)的适嗣渔圈进行拓展。2.3.3.2非平衡条件下NaOH.NaAl(Olt)4.H20体系活度系数计算模型文献曾报道40‘C时用Pitzer方程研究的铝酸钠溶液活度系数数据㈣,溶液瓯<5.m=0.01~10mol/kg。根据各温度下、备种吼下的平衡溶解度数据,利用公式(2.3)计算出相应的活度系数。将其与40’C时对应的活度系数迸行对比,然后计算出对应活度系致的差值,将这些差值进行拟合回归。得到活度系数差值(△Y)关于温度(t)的函数。所得到的函数形式为:中南人学硕十学伊论文第一章氧化铝生产热力等’数据的基础研究△YI=一O.0682+O.0038+(t-40)△托=O.0888一O.0045+(t-40)(口。S5.5)(吼>5.5)(2-4)(2.5)于是可得到非平衡条件下活度系数(T)关于温度(t)的函数:丫l=丫40一△YI(或△Y2)(2.6)100℃(t≠40℃),在公式(2.4)、(2-5)中,t为摄氏温度(℃)。其适用范围为:v=25吼=1.5~lO.5,伪20.Ol~lOmol/kg。根掘式(2-4)、(2-5)中相关参数盯知:当在低苛性比(口。S5.5)时.i舌度系故碹温宦升高而降低:当在高苛性比(吼>5.5)时.活度系数随.盛宦升高而增々Ⅱ。这可能与铝酸钠溶液中铝酸根离子的赋存状念有关,同时这一变化觇律与文献【蚰l的结论基本一致。2.3.3.3补,AIIOHH的计算模型基于非平衡条件下溶液中备组份平均活度系数的值(公式(2-6)),借助于Debye·Htickel活宦系故修正公式,假设NaAI(OH)4溶液为连续相,引入铝酸钠溶液的表现介电常数f’这一概含。将各种影唰归结于占’的变化f”I。计算对应的占’值。计算£’值的模型公式如下式:!!;!!+兰l!:!旦竺!。:+^.1¨一五I18高0706箸I丢丁s(2,7)×f’:“1‘二1式(:一7)中-P表示锥j酸钠溶液的的密室(g-cm。);肌为Na‘喷量肇尔浓霞(mol-kg“)。m。m~“Ⅲw,.+卅Ⅶw:T为绝对温度(K)。不同口。范围F’与7_和聊的函数关系,得到函数(2-8)、(2-9)。(1)、当瓴=1.5~5.5(含5.5)时,占’=bl+b2+ln(,)+b3*m+b4·朋2+b5*lna,+b6*m2*lnak/t15十b7‘聊3*In口。(2)、当口。=5.5~10,5(不含5.5)时,(2-8)一=bl+b2+t+b3+f2+b4+m+b5+肌2+b6+d^+b7+《+b8+用+a,h/t’5+b0·肌’吼不同吼下的调整参数值,列于表2-1。表2-1不同吼下方程(2-8)和(2-9)中的调整参数值Table2-1Values(2-9)ofadjustableparametersatdifferentcausticity.inEqnst2-8landt2-9、其中,式(2—8)、(2·9)的适用范围为:t=25~100"C,m---0.1~lOmol/kg。研究证明,当式(2—8)、(2—9)中的,分导U5829.47和5360.30时,可不考虑修正项系数打,此时br--O。中南人学硕十学伊论文第_二章氧化铝生产执力学数据的基础研究由式(2.7)可知,当活度系数大于1时,s’无实数锯。另外。当由(2-8)、(2-9)计算出的占’大于一定位时,修正项系数孵O。此时需汁尊打的值。经拟合,打足温室(,)、苛性比l口.)和苛碱浓宦(m)的函数,其函数形式为:(1)、当口。=1.5~5,5(合5.5),由式(2-8)计算出的s7>829.47时,br=cl+c2+,+c3+,+c4*m+c5*m2+c6+口。+c7+口:+c8+m’口。(2)、当口。-5.5~lO.5(不含5.5),由式(2-9)计算出的一>36030时.br=cl+c2+f+c34r,:/+c4+口。+c5+口:+c6+,以2/t+c7+mz*t(2-11)(2-10)不同a.下的调整参数值,列于表2.2。表2-2不同口。下方程(2-10)和I2-lI)中的调鏊参数值Table2-2Vahzesofadjuszableparametersatdifferen!ca).)slzcil"yinEqnsr2一JO)andt2-11)由此,结合Debye,HOcke[公式榴应参数的取值.我们可以借助于段模型得到NaOH.NaAI(OFI)4.H20体系活度系数计算摸型,如公式【2-7)所示。其中,当∥1829.47(瓯=1.5~5.5)或一<360.30(q=5.5~t0.5)时,6产O,占’可由式(2-8)、(2-9)计算而得;当s’>829.47(口。=1.5~5,5)或s’>360.30(口。=5.5~10.5)时,6辟0,b,可由式(2-10)、(2-II)计算雨得。对于不同的占’值,可进行如下解释l“l:1)口。=}5~5.5:当tr/和a.一定时,f’篮随,显度的升高两减小。这可能是出于随蕾温童升高,溶液体积嘭张,每单位体积分子数.成少引起的:当研和温叟相同时,F’值随q的增大而增大,这可能是出于随着口。的增大,溶液中体积较小的OH一离子增多、体积较大的AI(OHh"离子减少导致溶液单位体积所含溶质数增加昕引起的:当温宦和吼相同对,占’值短m的增大蔚增大。这可能是出于/.17趵增大导致溶液草位体积中所含溶质数增加所引起的。.2)瓯--5。5~10.5:当rt'/租aj一定时,s’值随温度的升高而增大。这可能是在高∞溶液中,体系中主要存在AI(OH)4.、AI(OH)6}、OH’等离子。当温宦升高时,溶液中体积较小的OH一离子增多、体积较大的AI(OH)4.离子减少导致溶液单位体积所含溶质数增加所引起的:当m和温度相同时,s’值随瓯的增大而增大;当温度和吼相同时,s’值随加的增大而增大。2.3.4适用性分析由公式(2-8)、(2-9)可知:运用上述NaAI(OH)4,NaOH—H20体系活度系数计算方中南人中项}学伊论文第一帝氧化铝t产执力学数据的基础研究法,计算出t=40"C,口。=l,5、3、5时不同l浓噎条件下的活度系数。将计算结累与文献p6l数据比较,比较情况见图2.2。由图2,2可以看出,在吼分剐为1.5、3.5时,经拟合馍型汁算出m=O.1~10mol/kg的活度系数与文献所列活度系数的相对误差clO%。计算值和文献值均表明:在低浓度时,活噎系数陋浓度的升高而降低:在高浓度时,活度系数谪浓度的升高而升高。计算数掂与文献数掘的吻合情况较好,说明浚模型对于计算铝酸钠溶液的活度系数育很好的适用性,这为热力学上分析铝酸钠体系各组份的行为提供了理论依掘。lOn∞驰叼l00090080>O000孙∞∞加070060O50二4¨,I二m,(tool·k91)/17/(mol·kg。‘)Ill0^_0000加m∞帅舳加的∞02468l012m/(mot·kg‘)图2-2估算模型丫Nn¨oHw计算值与文献值的比较Fig.2-2Comparisonof"fMHoH*calculatedbythecalculationmode!hthisstudy(I)∞d抽ereference(II12.4铝酸钠溶液中Na2C03、Na2S04、Na2H2Si04活度系数计算模型工业铝酸钠溶液中的Na20,除结合成铝酸钠(NaAl02)和NaOH的形态存在外,还以Na2C03,Na2S04,Na2HzSi04等形态存在,其活度的大小对于体系的热力学性质有较大的影响。研究其活度系数计算模型对于体系的热力学计算育一定的理论指导意义。±塑叁竺堡±芏芝堡壅2.4.I计算模型的建立2.4.1.I.笙三兰一兰些!堕羔!笠!兰鳖塑塑墨型竺生Na2C03活度系数计算模型的建立Na2C03—2Na++C012‘(2-1:)在溶液中,NazC03存在如下电离平衡:反应(2.12)的平衡常数为:足:华;‰。×席乞.一Ⅲt2"Ⅶ2m(2-13)式(2.13)中,哎,#.、口№.、口Ⅶ∥^分别为C032",Na+、Na2C03的活度;r№:cfJt为N永03的活度系数;棚:。.、m。,}分别表示Na+、C032-的重量哮尔浓度,单位tool出式(2.13)可得:kg"‘。‰:。,=(当反应(2-12)达到平衡时,(2—14)聊j口+。m‘rl△GT=一2.303RTIogK即:logK=△Gr/·2303RT(2-15)(2-16)经计算,反应(2-12)的△GT为:14934T二+I1.004‘In(T)4T(2_j7)/NGT;.2406373127-1224000/T-26.4213408’T+0.074至此,通过式(2.14)、(2-16)和(2.17)便可计算出某~温度平衡条件下Na2C03的活度系数。根据A.A.阿硌拉诺夫斯基提供的各种温嚏下纯喊在铝酸钠溶液种的溶解度数据∽1·用上述方法计算出了Na2C03相应的活度系数。如表2-3听示:表2-3不同条件下Na2cO,的,舌度系数里坐!::垒垡:i!:竺笪!塑!罂!兰g∑型壁竺竺生堕———一959595959548S411.262.051.73I.511.391.361.251.49I.491,49I.491.491,491.49226001430077.305300424022.4519.7584.4674.6665.9158.2251.614741O.1J250.1533020030.12880.11510.13070.12230.0385O.03920.03950.039600393O,03800.03652502j96jl4V,,·Vqu.u“7003080.02880.0296O.03030.0308O.03160.03t6003180.0312O·0305O.03010,0293O.02884025025030030030030030030030030030030030025747275.781083712004132.0414a,13156.9116994183.68197.85212.21227.15242.28I.49IA.93293.293.293.293293.293.293.293.293.293.2935.9l33788615774869.8761,53560450.4147.0744.2340.6738.2635,1383.50206202425026800278.00118.13131.41145.13159.34174.0918965205.75959525025025025025025025043.72翌!!!::!!.!兰—2羔——旦旦翌111:!竺:!!‘●_■■_■●_●■■■_■______■’__●____■__-●■■■■●■■●■-●_■■■■■■■■■■●_____●■■■-●■__■■____。‘■●。。。’———一一8中商人学硕十学伊论丈第一:章氧化铝t^执力学数据的基础枰窀然后通过数学回p。1的方法。拟合活度系数(TNa2C03)关于温度(T)苛性比(ak)、全碱浓度(mNa*IT))和碳碱浓度(mc0,2。)的函数。所得到的函数彤式为:logr、。,fJ,=pl+p2+InT+p3/T+p‘18logmNa+I.1+p5+lognk-032‘(2.18)其中,T为丌氏温宦(K),mN。+㈤为全碱浓度(moI·kg。’),mCO;2-为碳碱浓变(moI·kg‘’)pl~p5为调整参数,其值见表2—4。表2.4方徉2-is中的阔帮参数值Table2-4ValuesofadjustableparametersinEqnf:一f8】2.4.I,2Na:S04活度系数计算模型的建立运用蹩似fNa2COsi舌宦系数汁尊模犁逭芷的思路.采fflA.A.阿格拉若夫斯基堤哄的各种:益度下Na:SO{在铝酸钠溶液种的溶解度数掘l刈.计算出不同温度、不同碱,衣熏下Na2S04的活度系数,然后通过数学回归,拟合得到Na2SO。活度系数(ha2s04)关于温室‘T)苛性比{瓢)、全碱浓度{m。。+m)和硫酸碱浓度(11"15042。)的函数。所得到的函数形式为:logr~。:.v‘=kl+k2’In(T)十k3’logmha+I十寸k4+logms04z’(2.19)其中,T为丌氏温嚏(K),m、。+(tl为全碱浓度(mol·kg。),mⅫ2‘为硫酸碱浓度(mol-kg’‘),kl~k4为调整参数,其值见表2-5。表2-5方程2-19中的调整参数值Table2-5ValuesofadiustableparametersinEqnf:·19)2.4,1.3Na2H2Si04活度系数计算模型的建立在铝酸钠溶液中,sio:的溶解度比较小(~5叽),对于Si02来说,为稀溶液体系。故可直接运用Debye-Htickel模型计算其活度系数。log元:一掣(2-20)1+dB.1l“式(2-20)9tI表示电解质溶液的离子强度,单位为tool·kg一:乙、Z.分别表示f负9离子的价数:a是关于外离子与中心离子极限距离的参数,单位A;Ar、B,表示Debye.H0ckel参数,单位分别为kgln.mol地、kgV2.mol-U2.cm~,定义为:铲警岛=警f2—21)(2-22)其中,T裹,--,一.:…19r@K),P,p分别表示水的介电常数及畜宦(g·cm’3)。参数a是Debye-Hfickel公式应用的关淀参数。一股而言,a值大于晶体的离子半径之和,小于其水化半径之和。目4≯通常的做法。使用FB离子水化半径和阴离子晶体半径的加和。即Na2H2Si04的a值为:Na+离子水化半径(1.91A)与H2Si042.离子晶体半径(1.61A}相拓昕得I.91A+I.6lA=3.52A。2.4.2适用性分析运用式(2一l8)所述Na2C03活熏系数计算方法。计算不同温变、不同苛碱浓度时Na2C03在铝酸钠溶液中的平衡溶船度。将之与文献数扼州比较,比较隋况见表2-6。Fig.2-6ThecomparisonofsolubtlityofNa2C01calculatedbythemodelin廿IlsstudyandthereferenceT(℃)Na20kL∥L)ak—T(℃)表2-6Na,Oc(g/LlNa2COt平钎毒祥度计算埴与走欲值的比较Na:Ot,(g/L)文献。直22600143.0077.30530042.40224519.7584.4674计算值226.05142.6277.2952.7942,2022.39196994,6774,3766.11Sg4I5I7847.5743,85aL———————二—一文献一直计雉值1.491.49lNa,Ocfg/L)959595959595952502502502502502502504851.262.051.73I.5I1.391.36I.25l49I4925025025030030030030030023963257.47275,781083712004132.04144.13I56.9l1699418368197.85212.2I227.15242.283794359133,78861577.4869876I.5356.04504f47.0744.23406738.2635.1338.0335.97338285.6277.026吼476I.2055.7450.16468243.9940.4538.0434.934l4083.50206.20242.502680027800118.1313l4f1451315934174.0918965205.75493293.293,293.293293293293.293.29329:29663∞3003003003003001.491.49I.4965.9I58.225I6147.4l43.72I.491.49中南人学硕十学伊论文第_二蕈氧化锅士一L媳力学数据的墓础研究由麦2-6可以看出,数据的吻合情况较好,说明了该计算模型对于铝酸钠溶液中Na2C03活度系数的计算具有较好的适用性。此外,运用式(2.19)所述Na2s04活度系数计算方法,计算出在不同温度和不同苛碱浓嚏条件下铝酸钠溶i夜中Na2S04的活度系数,其数据与RobenoNa!SO。的活度系数数据m1摹本吻合,证明了本模犁的准确性。tPabalan等汁算的2.5本章小结l、与热力学相关的基本概念状态函数(StateFunction)、标准状奁(StandardStates),热容(HeatCapacity)、焓(Enthalpy)、熵(Entropy)、吉1{】斯自由能.(Gibbs常故K【ChemicalCoefncients)。ReactionEquilibriumFreeEnergyFunction}、化学反应平衡Constant)、活度(Activit)’)及活度系gt(Activity2、氧化铝生产中相关物质热力学数据的收集与甄别目日F.氧化铝生产媳力学数据唪已经收集了1199条记录。其中纯物质为839条,水溶液为360条。从数据库中可查到的物质热力学性质(即字段)有;编号(NO)、分子式(formula)、相念‘phase)、热容系数(a,b,c,d)、生成焓(AHO)、生成墒(△SO)、标准焓(HO)标准墒(SO)、温度范围(TempUp,TempDown)、晶体转变温度(T仃)、相变焓(Htr)、相交畸(Str),熔化温宣(Tm)、熔化焓(Hm)、熔化墒(Sm)、蒸发或升华温度(Tb)、蒸发冶(Hb)、蒸发熵(Sb)等。3、铝酸钠溶液活宦系数计算模型的研究1.8248lgy2一i—.4—-i—8—:—0719三6jI苓+6,,×106(pm)‘2l1.×(巴)…:p表示铝酸钠溶液的的密度(g-cm。);m为Na+质量唪尔浓度(mo|·kg"1).m=肘、。Ⅲwh+肌mw;r为绝对温度(K),占’为铝酸钠溶液表观介电常数,其与r和m的函数关系洋见2.3节。4,铝酸钠溶液中Na2C03、Na2S04和^fa2H2Si04活度系数计算模型的建立。主塑叁兰堡±堂竺堡奎墨三垦—墼些塑兰兰兰苎堕垡墅型王塑些些三塞.生第三章氧化铝生产热力学数据库的优化与完善3.1概述计算机软件的丌发,属于软件工程学的麓畴,人仃】已经建立了一整套t见范化的方法柬指导软件的丌发,而这些方法同样适胃于软件的优化或再丌发过程14”。软件性能的优化,在理论上是无穷函数,任何软件都有进行优化的空JlBJ。软件工程与数掘库理论的良好结合,是数据库软件丌发及优化的关键。数据库应用软件的优化工作涉及众多方面的内容,需要根掘具体情况而定。数掘库性能优化的一些主要内容,列于表3.1m'491。袁3-I敷据库陛能优化方向Table3·1Optimizingaspectsfordatabasecapabdity执化方向描述数据咋没计是堤毒鲁体性能的最有效建释。数据咩设计包括逻辑数据『车絮卡勾(船丧数据『芊设汁萃|f约柬}以及扬理r}能(缸衄嚣彖冼,对象伊萏和袁引'齐向优化府H{榨序殴汁优化实心I冗和l府H{数据『芊系统提供的实心f只和1.只的一些可_【f{选项提高数据辜性能I冗性能为了提高帮体性能.更改操仃系统(MicrosoftWindowsNT、MicrosoftWindows优化搬务器性能或MicrosoftWindows20001霜『数据唪(SQLSever或Access)98使圳lL确:殳计的裔,旬能够显#堤扁故据7莘眭能止确c殳计Ⅲ户样序,事务边界,锁定_-批处理的使川部雷丁府Hj群序殴计7旺镉我们注意到,除使用结构化,规范化的方法外,很多优化工作最终将归结到代码的编写上,属于编译时优化的7苞畴。对此,编程人员往,芏有一些经验性的技巧,这些编写代码的技巧对应用软件优化有很大的帮助。例如:对VB应用程序而言,可采用减少控件数目、用杯签控件代替文本框、共享图片等方法实现程序大小的优化:选用最佳的变量、数据类型实现速度的优化等。中南入学谚十学竹论文镣二章氧化锅t产扎力学数据牢的优化与完善3.2旧版本数据库存在的问题及相应优化解决方案旧版《氧化铝生产热力学数花库》在过去的使用中表现出了强大的生命力,刚时也暴露出了~些问题。新版的《氧化铝生产热力学数拢库》梅疰改善整个数嬲摩性能上做出尽可能多的努力.并勺口入新的热力学计算模型。藏版设计目标旧坂相同,即:桌面型、}:}研用专业数据库应用软件。这砷产品定也}要是由实际需求所决定的。眩软件属Fj史奇软件.主要用F科研工作。用户群为个人或,J、闭体.FJ时具何专业性强.数据表简单、数瓤晕不欠≮特点。由于个人计算机性能价格比的迅速提高.给了桌面犁数据库以极大的发展空删,也为我们的软件丌发和使用提供了充分的技术支持。下文将对旧版数据库存在的主要问题及对应的优化完善方案和新增内容作分类介绍,脚时为了增加本文的可溃性弓简洁性,省略了数据『聿.程序源代码优化过程的的jf细介绍。1.1关于精简性1.存在的问题在=致掘章的某些界面和J力能中.存在苦:吾三冗长、J力能节复的现象。如数掘计算界图3-1旧版数据库敷据计算界面F|g.3-1ThecaJcuJaZzonmoduleofformeredition3面中胡“纯扬;专(离产)在温瞧T下的焓Ht礴st厦生成}1出自}已经包含了”离子在温熏T下的焓Ht墒St硬屯咙自由能”白勺计雉』由能.如罔3.I昕示:中南大学硕十学1奇论文第二章氧化锅生产拽力学数据车的优化与完善2解决方案针对旧版本数据库在界面方面存在的问题.在不减少功能和影响操作的前提下,对数据库界面进行了简化,对已有的重复功能进行了整合,对其中的浯言进行了精简。以数据计算模块为例,优化后的界面如图3-2所示:图3-2新版数据库数据计簟界面Fig,3·2ThecalculationmoduleofneweditJ(}n3.2.2关于完整性1.先酊数据库基础数据覆盖的内容体系帽图:AlzO,·si02系、CaO-Ah03系、CaO-Si02系、Na20.A1203系,Na20.Si02系。数据库物质记录数:Data表953条,IblAqueous表233条,tblPm'ea表720条。铝酸钠活度系数估算模型适用范围:温度t--40"C~100℃,苛性比ctK=I.5~5,苛碱浓度mHaoHg--05~10mol/kg。2.目前数据库基础数据包含的内容体系帽图:Al∞3一SiOz系、CaO-A120j系、CaO.gi02系、caO.Fe203系、CaO-FeO系、FeO—SiOzg,Na20-A1203-H20系、CaO—A120rSi02系、CaO—FeO-Si02系、NazO.A1203-SiO:一H20系、Na20一CaO—A120rSi02-H20系、Na20一CaO.A1203·C02.H20系。数据库物质记是数:目前,Data表1199条记录。其中纯物质为839条,水溶液为360象。铝酸钠活度系数估算模型适用范围:温度t=25"C~300"C,苛性比aK=1.5~10.5,苛碱浓度rnH目赇=0.1~10mol-kg.‘。±妻叁堂堡±兰生堡壅苎三!墨!兰塑竺兰兰塑!!堕!竺垡!!生圭堂3.2.3关于反应方程式的配平1.存在的问题配平后的反应方程式的系数采用小数表示,这可能给后续计算带束一定误差。如式(3-1)所示:(O.67)Fe+(1)C12=(067)FeCl3(3一1)另外:在某些方程式的配平请求中,程序无法完成配平且弹出出错提示。2优化方案为解决方程式配平中存在的问题,在方程式热力学数掘计算模块中添加了“手动输入系数”和“自动配平系数”功能模块供用户选择,以解决可能出现的问题。其界面如图3—3所示:图3.3反应的AGt计算界面图F嘻3-3TheformsofAGtofreactions3.2.4关于准确性1.存在的问题经原版数据库估算出来的数据,部分值与文献值、手动计算值误差较大,这影响到了数据库软件的应用推广价值。2.优化方案一是对原有的基础数据进行了一一核对,一保证原始数掘的准确性:二是增加了物质的△H298、AS298、AG298等热力学数据.以扩大物质热力学数据的查询范围及增加其计算方式:三是对于有离子和纯物质参与的化学反应,所有物质均统一1采用了式(1-5)、(1石)和(1-7)的计算模型,保证热力学数据的一致性,增加计算结果的准确性。中南人宁硕十学忙论文第二章氧化锅士产热力'-7数据『芊的优化与完簿3.2.5关于程序非正常出错原版数据库程序在运行过程中,有时会突然中断运行并非J下常退出主程序。若需继续操作,需重新启动程序。如图3-4所示:囿3_4数据库程序出错界面图Fig.3-4Theerrorsofformerediuon2|优化方褰通过记录程序出错时的错误代码,找出错误原因,修改相应的原程序并调试,防止程序出错非一常退出现象的缆生。优化思路如图3-5所示:困3-5敷据库程序出错优化思路图Fig.3-5Theoptimizedmca¥ureforthecrror¥ofthedatabase3.2.6关于数据库操作方法1存在的不足数据库包括了3个查询模块(分子式、元素、体系查询),2个类别(水溶液、纯物质)的3种数据编辑方式(更新、添加、删除)和9种数据计算类型。但未对其中的查询方式、输入规则、操作流程等进行必要的说明,对初次使用者来浣有一定难度。中南夫学硕十学佛论文第二章氧化锅士产媳力学数据『芊的优化与完善2解决方案为了使耳】户可以全面系统地了解系统的应用』力能.并学习到其基本使用厅法,特编制了~月余字的《氧化铝生产热力学数据庠(TDAP)的用户使用手册》。从数掂库系统的安装建立到各部分、备模块的使用方法和具体操作步僳都进行了详细说明。如复杂硅酸盟热力学性质拆分估算模块的使用说明如下:图..1-6复杂硅蝮盐热力学性质的折分估霉界面固Fig.玉6Theinterfaceofthermodynamieestimalionofcomplicatesilicates操作步骤:1.输入待估算物质的分子式;2.在酸碱比合物化合产物框中,按照复杂硅酸娥中各氧化物的酸碱性强弱,由强至弱输入各个氧化物的两两化合产物及其系数(系数输入到稼色栏);例如:求CaO·2Si02·A1203·2H20在酸碱化合物化含产物框中应输入:l1CaO*SiO:AI203+sio:注意:此处元素应区分大小写。3.输入计算温度或温度范围;4,点击“计算”按钮进行计算。计算过程中,程序会询问用户某一确定物质的具体晶型或状态,用户根据需要选择即可。选择完毕后将显示计算结果;5.查阅计算结果后,用户点击“完毕”按钮可返回计算窗体:6.若要继续计算.单击“计算下一个”即可。中南大学硕十学伊论文第二章氧化锅t产执力学数据庠的优化与完善3.3其它方面的完善3.3.1元素精确查询原版数据库在进行元素查询时,不是按照元素进行查找。面是按照字母查找,即将包含该元素字母的所有物质均查询出来,即模糊查询。如查询元素“S”时,物质“Si02"也随之被捡索出来。优化时运用“If,…一Then”语句对用户所要查询的元素进行了筛选,达到精确查询的目的,从而解决了类似的问题。3.3.2安全性为保.匹数据库的安全性,设冒了数据库安全密玛.包括数据库登陆啻码、数掘库修改添加密码、Access表打丌密码等。3.3,3程序说明文挡为了给以后的程序完善和升级丌发提供方便,制定了6000余字的数据库源代码说明文档,对在丌发程序中采用的算法逻辑、层次结构、模块『日J的连接、重要变量的使用等情况进行详细的说明。3.3,4数据库用户手{j{}为了方便用户使用数据库.以图文的形式编制了10000余字的《氧化铝生产热力学数据库用户手册》,对数掘库的基本框架、主要功能、使用操作方法、注意事项等情况进行了详细的说明。3.3.5帮助氧化铝生产热力学数据库提供了标准Windows风格的帮助系统,以帮助用户进一步了解和使用本数据库系统。帮助系统的目录是按照数据库的内容层层展开的。当用户进入某个功能模块后,悠获得有关帮助,请直接点击窗口左下角的“帮助”按钮,获得相关帮助内容;也可根据需要在目录中进行选择,或者直接输入关键字查找到所需主题,如图3.7(a、b)所示:中商夫学硕十学伊论文第三章氧化锅生产热力学数据痒的优化与完善Ca)图3-7数据库的帮助系统(a_目录式帮助F嘻3-7b.捡索式帮助)Thehelpsystemofthedatabaseb·Searchesmode)(a-Direclorymode本数据晖的帮助文件是用VisualWorkshop编写的。编写帮助文件的步骤为:Basic6.0所携带的号门编写帮助系统的工具HeI口1、创建帮助主题文件(+.rtf)。帮助主题文件保存了所有的帮助信息,它就象一本书,包括目录和~个个章节,各章节之间通过热区束链援。热区可以是文字或图标,实现主题J日J的跳转。该文件可用Microsoftnf格式。2、创建帮助目录文件(+.cnt)。帮助目录文件是ASCII文件。在帮助主题显示Wbrd编写。保存成框中以{对状目录的形式向用户展示帮助内容,同时在-索引’条中用户可通过关键字进行内容查询。启动应用程序hcw.exc即HelpWorkshop可编写帮助目录文件。3、创建帮助项目文件(+.hpj)。项目文件是用束编译帮助文件的ASCII文件,它将帮助主题文件、目录文件及相关参数传递给帮助系统4、编译和生成帮助文件(·.hip)。直接点击界面上的‘编泽’按钮,该应用程序可自动生成帮助文件。在VB6中通过API函数可对生成的帮助文件实行调用。语法如下:DeclareFunctionWinHelpLib‘‘User32”Alias“WinHelpA”(ByvalhwrldAslong.ByvallpHelpFileAsString,ByvalwCommandAsLong,ByvaldwDataAslong)Aslong29中南大学硕十学付论文第二章氧化铝生产热力学数据『芊的优化与完善3.4本章小结本章从用户界面的精简性、计算结果的准确性、热力学数掘的完整性、反应配平程序,铝酸钠活度系数估算模型的建立等方面角度分析了旧版数据库存在的主要问题,并针对相应问题提出具体的解决方案,论述了具体的优化过程以及对优化的评述,同时进行了新旧数据库的使用对比,晚明新版数掘库的优化成果及明显技术优势。为方便数据库的使用和升级,编制了《氧化铝生产热力学数据唪用户手册》、《氧化铝生产热力学数据库程『芋说明文档》和帮助系统。中南大学硕十学竹论文第四章Na!O.A1203-H:O系热力学分析第四章Na20.A1203.H20系热力学分析不言’而喻,要深入研究氧化铝生产过程中内在的化学规律,用理论柬指导生产非常重要。热力学是研究各种变化过程能量转化关系,以及过程进行的方向和限宦的科学。热力学能为实验研究指明方向..避免实验工作中的盲目性,对优化生产条件、提高生产效率、降低能耗及原材料消耗育重要的指导意义【5051·525”,从而为发展新工艺和改善老工艺以及为有预见性地控制现育生产提供理沦依掘I“l。玎耳法是碱法生产氧化铝中的其中一种i虎程。其是直接利用含有大量游离苛性碱的循环母液处理铝土矿,溶出其中氧化铝得到镭酸钠溶液的,往铝酸钠溶液中添加氢氧化铝(晶种)经长时问搅拌便可分解折出氢氧化铝结晶【5”。拜耳法用在处理低硅铝土矿,特别是用在处理三水铝石型铝土矿时,流程简单能牦低,作业方便,产品质量高,其经济效果远非其它方法所能媲荚。目前全世界生产的氧化铝和氢氧化铝,有90%以上是用玎耳法生产的i”;。4.1三水铝石溶出过程热力学分析玎耳法的实质就是使下一反应在不同的餐件下朝不同的方向交替进行f55】:A1203(1或3)H20+2NaOH+aq一2NaAl(OHh+aq(4一1)铝酸钠溶液在氧化铝工业生产中扮演着非常重要的角色。在通常的热力学讨论中均为考虑活度对反应吉南额自出能的影响。然而铝酸钠溶液活度系数的计算不仅有助于氧化铝生产理论的完善,而且将大大促动氧化铝生产新工艺的研究和丌=发,达到降低成本、节能降耗的目的。现碾据已建立的铝酸钠活度系数计算模型15”,以三水铝石在NaOH溶液中的溶出反应为例,对其进行热力学分析。三水铝石在NaOH溶液中的反应可表示为:该反应的反应吉却斯自由能为:△嘶一f+器,s专≯2△旰+罴ts警.4t(OH),4-NaOH=NaAI(OH)4(4.2)c㈣式(3-3)中,蛳aAI{oHH,(]NaOH分别为NaAI(OH)4、NaOH的活度;YNaAI(oH)4、¨a0H分别为NaAI(OH)4、NaOH的活度系数;mAI(0H斗、moH.分别表示AI(OH)4.、OH’的重量擘±堕盔兰堕±兰堡垒壅尔浓度,单位molkg"‘。笪堕兰型!:竺:垒k旦!:旦!Q墨垫垄竺坌塑根据式(3.3)计算出三水铝石不同温度时在NaOH溶液中的反应吉钿斯自由能值(Na20k=2309/L,ak-2.3),同时将其与不考虑活度时的吉却斯自由能值进行了比较·如图4-1所示:80008000{4000昌20000-2000250270290310330丁,’K350370390410图4-l三水铝石瘩出反应的AG-T图(Na20k=2309/L.OK=2.3)Fig.4-1ThefiguresofAG·Tofthedigestionofalumina(NazOt22309/L.ak22.3)由图4.1计算可以看出,三水铝石在其本身热力学悬度范围内和NaOH反应的吉11】斯自由能值随着温度的升高而减小:到达一定温度后,反应的吉南斯自由能值为负值,说明反应在一定温度时能自发进行。同时,对比考虑活童和不考虑活度时反应吉布斯自由能值的大小可知,活度的引入使反应的吉布斯自由能值有所降低。在相同aK(aK=2.3),不同苛碱浓度下三水铝石溶出反应的△G-T图如图4-2所示:趸2≥t昌042—4280300320图4-2Fig.4-2Na:ok对三水铝石溶出反应AG-T的影响(O.K=2.3)ollEffectofNa20kAG-Tofthedigestionofalumina(aK---2.3)中商人学够十’宇侍论文第四章Na:O.A1201·H10系热力学分析从图4—2可知,在aK不变的条件下,增加苛碱浓度,三水铝石溶出反应的吉南斯自由能值减小。即是说,当aK恒定时,苛碱浓度的增加将有利于溶出反应的进行。在相同苛碱浓度(NazOk=2309/L),不同O.K下三水铝石溶出反应的AG.T图如图4-3所示:星≥t<图4-3Fig.4_3aK对三水铝石瘩出反应△G—T的影响(Na2Cll,=2309/L)onEffectofaKAG—Tofthedigestionofalumina(Na:Ol,=2309/L)从图4.3可知,在苛碱浓度不变的条件下,随蕾旺K值的增加.三水铝石溶出反应的吉白斯自由能值碱,j、。即是浣.当苛碱浓度恒定时,口。值的增加将更有利于溶出反应的进行,进入溶液中的A1203的数量随之增加,文献报道也证实了这一结论㈣。4.2三水铝石平衡溶解度估算铝酸钠溶液是氧化铝生产的基础,有许多研究者曾对在各种温度下的NazO.AIz03.H20系的平衡情况进行过研究,得到的结果大同小异。由于铝酸钠溶液活度系数方面研究的限制,过去人们在计算NazO.鲇20,.H20体系A1203溶解度时,往往以浓度代替活度,忽略活度系数的影响。数据库提供了加入活度系数的平衡溶解度计算模块,束解决这一问题。这些平衡溶解度计算结果对于工业生产有一定指导意义,同时也可以对活度系数计算模型进行检验。数据库计算出的在30~100"C时,A120)在NaOH溶液中的平衡溶解度数据,见表4.I:状念图,如图4.4。表4-l计算出的三水铝石瘩解度数据N赴o(叽’可鬲—丽百—面丁1秽地告下、丽厂丽百面Table4-1GibbsitesolubilitydatacalculatedbythemethodjnthisstudyNaz0.¥/L图“不同温度下的Na20.A1203-H:0系状态图Fig.4-4ThesolubilityCurvesofNa20一A1203-H20systematdifferenttemperatures从图中可以清楚地看到,随温度的升高,等温曲线的曲率逐渐减小,即越来越拉直。这一规律符合一般常识,从而在一定程度上证明了论文所提出的活度计算模型的正确性。34中南人学硕十宁伊论文第四章Na20-AI:03-H20系热力学分析进一步,将不同温度叮A1203在苛碱中的平衡溶解度数掘与文献数据比较,文献选取Russell的数掘【”I,以40"C、70"C、100"C为例.比较隋况见图4-5~4.7。160140120≤100∞占80;604020005000l50Na20.gL200230300图4.540'C平衡露解度计算结襄与丈献IRussell)的比较Fig.4-5ComparisonofequilibriumsolubilitydatacalculatedwiththeestimationmodelinthisstudyandreportedbvRussellat40。C250200d∞150宅;100500O50100150200250300N赴0tglL图4.670"C平衡客解度计算结果与文献(Russell)的比较Fig.4-6ComparisonofequilibriumsolubilitydatacalculatedwiththeestimationmodelinthisstudyandrepottedbyRussellat70"6中南人学硕十学何论文第四蕈Na20.AI:O,-H:O系热力学分忻300250200芒呻苦150o《100500050100150Na?0.g/L200250300图4-7100℃平衡溶解度计算结果与丈献(Russell)的比较Fig.4.7ComparisonofequihbriumsolubilitydatacalculatedwiththeestimationmodelinthisstudyandreportedbyRussellat100"C由图4-5~4.7可以看出,数据的吻合情况较好。溶液离子强度为O.5~11mol·kg“(已做单位换算)左右,计算数据仍然可以很好的与文献值吻合,这说明估箅模型对于低浓度和高浓度铝酸钠溶液都有很好的适用性,完全突破了Debye—Hfickel原始公式对溶液的限制。这也充分证明了该计算模块的实用性和准确性,并为相关研究提供了新的参考。4.3本章小结本章以三水铝石的溶出为例对Na20.A1203.H20体系进行了热力学分析,包括:1、比较了考虑活度和不考虑活童时三水铝石溶出反应的吉伟斯自由能值。结果表明,活度的引入将使溶出反应的吉却斯自由能值降低;2、考查了Na20k和毗对其溶出过程热力学的影响。结果表明,当其它条件不变时,苛碱浓度的增加或提高咏值都将使反应的吉布斯自由能值降低,即有利于溶出反应的进行;3、依靠热力学计算模块——‘‘铝酸钠溶液体系平衡溶解度及活度系数的估算”功能对三水铝石30~100"(2的平衡溶解度数据进行了估算,并与文献值进行了比较,结果基本吻合。这进一步说明了该计算模块的实用性和准确性。中南人学硕十学何论文第百蕈Na20-CaO—A120t-Si02·H20系热力学分析第五章Na20.CaO—A1203.Si02一H20系热力学分析5.INa20—A1203一Si02.H20系热力学分析5.1.1A1203·2Si02"2H20溶出反应热力学分析A1203·2S102·2H20在铝土矿中存在多种形冬,其体系热力学数握温度范围及相变如下:A1203+2Si02’2H20:S(高岭石):298-861K:S(逾丌石):298~700K:S(多水高岭土):298—700K:其与铝酸钠母液发尘如下反应:AI_,03’2S102。2H20+6NaOH+H20--*2NaAI(OH)4+2Na2Si02(OH)2(5-1)哪铷咖哪一三.一瑚脚Ⅲ㈣250300350400T/K450500550600图5-1高岭石溶出反应的△G_T图(Na201,=2309/L,ak-2,3)Fig.5-lAG·Tfigureofthedigestionofkaolinite(Na20k=2309/L.0.K=2.3)ⅡⅡ5-1可知,随温度的升高,反应的吉白斯自由能减小,但相较而言。,其反应的吉白斯自由能值的大小为:高岭石>迪开石>多水高岭土,说明这三种分子式相同的物质在此温度条件下和碱液反应能力大小顺序为:多水高岭土>迪开石>高岭石。.高岭石是Si02在铝土矿中存在的主要形态。现以高岭石为例,考查Na20k和aK对其溶出反应吉布斯自由能值大小的影响,结果如图5·2、5·3所示:中南人学硕十学何论文第百章Na20·CaO-A1201-si02-H20系热力学分析图5--2Fig.5-2aK对高岭石瘩出反应的△G-T的影响(Na20k=2309/L)onEffectofaKAG—Tofthedigestionofkaolinite(Na:Ok=2309/L)图5-3Fig.5—3Na20k对高岭石容出反应AG·T的影响(aK=2.3)onEffectofNa209AG·Tofthedigestionofkaolinite(ctK--2.3)Ⅱ口5-2、5-3可知,随着Na20k浓度和aK的升高,反应的吉向斯自由能减小,即增加苛碱浓度和aK均有利于高岭石溶出反应的进行。5:1.2钠硅渣生成反应热力学分析以碱法生产氧化铝关键在于铝硅分离,在氧化铝生产中硅主要是以钠硅渣(DSP)的形式除去的。生成钠硅渣的反应式如下:2NaAI(OH)4+1.7Na2Si02(OHh。Na20‘A1203’1.7S102‘2.3H20+3.4NaOH+1.7H20(5-2)中南大学硕十学侍论文第五章Na20.Ca(9.A1201.Si02.H20系热力学分折现考查在不同情况下钠硅渣,土成反应的热力学性质。1.碱浓度(Na20k)及苛性比(nK)恒定,二氧化硅浓度(Si02)变化取Na20k=2309/L,aK=2.3,考察当SiO:在O.1~39CL之f’日J变化时,反应(5-2)的吉和斯自由能在T=300~570K之问的变化情况。结果如图5-4所示:4o_:叫薹一』喝司啦Ⅲ啦280330380430480530580T/K图54不同SIO,.浓虚下反应(5-2)的△G·T图(Na:Ok22309'L.O,K32.3)Fig.5.4AG-Tcurvesofthereaction(5·2)atdifferentconcentrationofSiO:(Na:Ok=230∥L.QK=2-3)从图5.4可以看出.Na20k及aK一定时,Si02浓度的愈高,反应(5-2)的吉角斯自由能愈负,即更有利于钠硅渣的生成。2.苛性比(ⅡK)及二氧化硅浓度(Si02)恒定,碱浓度(Na20k)变化取U.K=2.3,Si02=1.5∥L,考察当Na20k在150~3209/L之『Bj变化时,反应(5—2)的星≥t司圈5-¥不同碱农度下反应(5-2)AG-T图(aK----'2.3,SiOr=1.59,L)Fig·5-S△G-Tcurvesofthereaction(5·2)atdifferentconcentrationofNa20R(CtK----2.3,sin2=1.59,L)吉卸斯自由能在T=300~570K之l’日J的变化情况。结果如图5—5所示:±塑盔兰堡±兰旦堡壅兰互兰堕丝旦:壁唑垒兰翌型型生墨垫查:兰!:里从图5-5可以看出,Si02浓度及呶一定时,Na20k浓度的愈高,反应(5-2)的吉南斯自由能愈负,但影响不显著。3.碱浓度(Na20k)及二氧化硅浓度(Si02)恒定,苛性比(aK)变化取Na20k=2309几,sio:=1.59/L,考察当IIK在1.5~3,5之『Bj变化时,反应(5-2)的吉却斯自由能在T=300~570K之『日J的变化情;兄。结果如图5-6所示:o巧E哪≥』dq啪哪嗡280330380430480530580图蹦不同0tK下反应(5-2)的AG-T图(Na20k=2309/L,SiOz=1,5∥L)Fig.5-6AG-TCUIⅣe¥ofthereaction(5-2)atdifferentIlK{Na20,2230e,/L·Si眈=1.Sg/L)从图5-6可以看出,NazOk和si02的浓度一定时,aK愈高,反应(5-2)的吉南斯自由能愈正,即愈难生成钠娃碴。4,生成钠硅渣时Si02的乎衡浓度在铝酸钠溶液中,碱浓度、氧化铝浓度、二氧化硅浓度、温度、品种及杂质等均会影响钠硅渣的生成速度和组成,继而影响钠硅渣的稳定性和二氧化硅在铝酸钠溶液中的平衡浓度。尽管国内外对形成钠硅渣后二氧化硅的平衡浓度避行了很多的实验研究,并提出了较多的二氧化硅平衡浓度经验公式,但由于实验者采用的研究条件不同,得出的实验结果相差较大。从理论上计算或研究形成钠硅渣后二氧化硅的平衡浓度至今未见报道,主要原因是难以制备具有确定化学式【或分子式)的钠硅座,进而无法在观有的实验技术条件下铡得各种钠硅渣的热力学数据。基于此,利用“氧化铝生产热力学数据库软件”中提出的“组成相似的一类化合物其组成与热力学数据存在良好的线性关系”来计算钠硅渣的热力学数据,以研究形成钠硅渣后铝酸钠溶液中二氧化硅的平衡浓度。据钠硅渣(设为Na20·A120y2Si02·2H20)热力学数据的估算值,结合生成钠硅渣的反应方程式,可计算生成钠硅揸时二氧化硅的平衡浓度,即2Na++2AI(0H)4’4-2Si02(OHk2r_Na20‘A1203·2S102’2H20+40H’+2H20(5-3)中寓人学面十中付论文第百章Na20.CaO.AI:O,.Si02-H20系热力学分析Gibbs自由能值为:△Gr。=AG‰一(T-298)△s瑚0+gssAC;万一丁E。垒争7'(5—4)在平衡时,△G≯=一RTinK(5-5)根掘苛性比的定义有:aK一1=x【OH‘】/x【Al(OH)41(5-6)则x[Si02(0H)22.】=(QK—I)x[OH‘]exp(AGro/2RT)/x[Na+】(5—7)掘式(5.7)可计算二氧化硅的平衡浓度,如图5.7。图s-7在不同温度下生成纳硅查时二氯化硅的平衡浓度Fig,5.7TheequilibriumconcentrationofSiO:inthesolutionofsodiumhydratealumino-silicateatdifferenttemperature图5.7说明了对于具有固定结构的钠硅渣,其二氧化硅平衡浓度是随着温度的升高而增大。据此计算出的二氧化硅的平衡浓度与文献的实验测定值【59删进行对比得知,在不同条件下生成钠硅渣时二氧化硅的平衡浓度的基本规律一致,计算结果与文献测定结果基本吻合。这从另一个方面也说明了钠硅渣热力学数据计算的可信。5.2Na20—A1203.CaO.Si02一H:O系热力学分析5.2.1水合硅酸钙在铝酸钠溶液中反应行为的热力学分析水合硅酸钙(以C.S.H表示,下同)是碱法氧化铝生产中的一种重要的含硅物质,在烧结法生产氧化铝过程中,二次反应导致铝和碱的损失主要是由于熟料浸出时水合硅41中南人学硕十中伊论文第五章Na20一CaO-A1201.SiOz—H20系填力学分析酸钙和铝酸钠溶液问的反应造成的,但其与铝酸钠溶液的反应行为不甚清楚。因此,研究水合娃酸钙翔铝酸纳溶液『日J的反应行为具有十分重要的理论和使用价值。C.S,H与铝酸钠溶液问通常被认为发生如下反应:mCaO·nSiOz·xHzO+mC032.+mH20=nH2Si042"+mCaC03+2(m-nlOH‘(5.g)(5-9)mCaO。nSi02’xH20+2nOH。mCa(OH)2+nH2Si042叫x-m)H20利用数掘庠分别计算了反应式(5-g)、(5-9)的反应吉伟斯自由能随反应温度的变化,如图孓8及图5-9所示:图5-8水合硅酸钙与Na2C03反应的△Go-T图Fjg.5.3AGO.TfigureofreactionofhydratecalciumsilicatewithNa2C03图5-9水合硅酸钙与NaOH反应的△Go-T图F噜.S-9AGo.TfigureofreactionofhydratecalciumsilicatewithNaOH中南人学硕十学伊论文第百蕈Na20—CaO—A1201-Si02-H20系孰力学分忻从图5.8可知,在热力学上6CaO·6S102.H20、CaO·Si02·H20、2CaO·Si02·1.17H20、2CaO·Si02*I.6H20和CaO·2S102-2H20均能按式(5-8)与Na2C03发生反应。并随着温度的升高水合砬酸钙越易分解:图5—8与图5-9相比,在同等条件下,Na2C03更易分解水合硅酸钙:从热力学束看,钙硅比(窜尔比:CaO/Si02)高的水合硅酸钙比钙磕比低的更易分解,其反应稳定性顺『芋大体为:CaO·Si02-H20>2CaO·Si02-1.17H20>2CaO'8i02·1.6H20。5.2.2由Ca(OH)2生成水化石榴石的热力学分析Ca(OH)2形成水化石溜石的反应方程式为:3C“OH)2+2AI(OH)4。+xSi02(OH)2。—3CaO·A1203‘xSi02·(6·2x)H20十2(工+1)OH’+2xH20式中x为Si02的饱和系数。x=0.1~2。饭鳓肫鲋舶醣沩:AG=&GT.+罴lg麦【5-1O)在式(5-10)中,当Si02饱和系数X的值不同时,参与反应的Si02的量也会不同。这样计算出柬的热力学数据孰无法进行相关比较和分析。为此,必须确定其反应的统一标准。这罨,统一设Si02的量均为“l”。以此为基准,讨论生成不同饱和系数(x=O.5,1.1.5.2)水化石榴石的反应热力学性质(kJ.molo):6Ca(OHh+4AI(OH)4.+Si02(OH)2。—·2(3CAO·A1203-0.5S102‘5H20)+60H。+2H:O(5—11)3Ca(OH)2+2AI(OH)4‘+Si02(OH)2‘。一3CaO‘A1203·Si02·4H20+40H’+2H20(5-t2)2Ca(OH)2+4/3AI(OH),"+Si02(OHF’一2/3(3CAO·A1203·1.5S102·3H20)+10/30H’+2H20(5-13)1.5Ca(OH)2+AI(OH)4‘+Si02(OHh。一O.5(3Ca0。A1203·2S102。2H20)+30H.+2H20现考察在不同情况下反应钙硅渣生成反应的热力学性质。(5-14)1.碱浓度(Na20k)及苛性比(ctK)恒定,二氧化硅浓度(Si02)变化取Na20k一2309/L,CtK=2.2,考察当Si02在O.1~3蛆之11日】变化时,反应(5-“)~(5-14)的吉布斯自由能在T=300~570K之『日J的变化情况。结果如图5.10(a~b)和图5.1I(a~b)所示:柚2020—020≥昌-20一10毒一20宅昌’10—60—60280330380130T/K'180530580.802803303804301fX480530580(a)2020(b)00一i。二0i。20主三乏‘10-oO昌+40—60.30280330380130T/K180530580—80280330380啪啪啪(c)‰∽(Si02(g/L):a-O.1.b-O.5,c-2.5,d-3.0)same图5-10相同Si02浓度下不同x值反应(5-10)的△G·T图(Na20k--2209/L.OtK=2.3Jofthereaction(5·lO)attheSi02concentrationanddifferentvalueofXrig.s.10△G.Tcurves