Study on cation and electronic potential distribution in montmorillonite double layer considering ion size effects
-
摘要:
离子型稀土矿浸出是一个双电层重新平衡的过程,离子尺寸效应对双电层中阳离子和电势分布具有重要影响。本文基于Gouy-Chapman-Stern双电层模型,考虑了离子大小对体系静电自由能所产生的影响,耦合了考虑离子尺寸效应的非线性Boltzmann方程与泊松方程,建立了sPB方程;利用该耦合方程研究了离子尺寸效应对双电层中离子和电势的分布影响,结合经典Poisson-Boltzmann方程(cPB)的计算结果,与杯浸试验结果进行了对比分析。结果表明:NH4+和La3+在蒙脱土表面的吸附存在差异,蒙脱土对NH4+和La3+的饱和吸附量试验值分别为0.635 mmol/g和0.305 mmol/g。基于cPB方程,NH4+和La3+的计算吸附量分别是实际的17倍和11倍;而基于sPB方程,NH4+和La3+的计算吸附量分别为1.75倍和1.25倍。sPB方程所计算的Stern电势小于cPB方程的计算值,且sPB方程所计算的Stern电势随浓度的变化程度不大。离子水化作用对离子的吸附量和Stern电势的变化几乎没有影响。综上,考虑离子尺寸效应的sPB方程能更好描述双电层内阳离子和电势分布情况。
Abstract:The leaching of ionic rare earth is a process of rebalancing the electric double layer, and the ion size effects have an important influence on the cation and potential distribution in the electric double layer. This article coupled the nonlinear Boltzmann equation considering ion size effects with the Poisson equation (sPB) by taking into account the influence of ion size on the system's electrostatic free energy, based on the Gouy-Chapman-Stern double-layer model. This coupling equation was used to study the effects of ion size on the distribution of ions and electric potential in the double layer. Combined with the results of cup leaching experiments, the calculation results were compared and analyzed with the classical Poisson Boltzmann equation (cPB). The results showed that there were differences in the adsorption of NH4+ and La3+ on the surface of montmorillonite, and the measured saturated adsorption amounts were 0.635 mmol/g and 0.305 mmol/g, respectively. Based on the cPB equation, the calculated adsorption amounts of NH4+ and La3+ were about 17 times and 11 times the actual adsorption amounts, respectively. However, according to the sPB equation, they were about 1.75 times and 1.25 times, respectively. The Stern potential calculated by the sPB equation was significantly lower than that by the cPB equation, and the variation of the Stern potential calculated by the sPB equation with concentration was not significant. Ionic hydration had little effect on ion adsorption capacity and Stern potential. Overall, the sPB equation considering ion size could better describe the distribution of cations and electric potential within the electric double layer.
-
氧化铝生产过程中所产生的强碱性废渣被称为赤泥。因生产方法和铝土矿品位的不同,每生产一吨氧化铝大约要产生0.5~2.0 吨的赤泥[1]。2019年中国氧化铝产量7 230.16万吨,伴随产生的赤泥高达1亿吨,累计堆存量超过13亿吨,占地超过12万亩[2]。目前,世界上对赤泥采用的主要处置方法是对其进行筑坝堆存,长期堆积存在溃坝风险。赤泥堆放占用大量土地的同时其中存在的大量碱会污染地下水体与土壤。此外,干燥赤泥遇风会造成空气污染。如何对赤泥进行规模化、无害化处理,仍是急需解决的问题[3]。
赤泥作为氧化铝行业产生的大宗固废之一[4],各国研究人员对其综合利用进行了大量研究,研究方向主要包括:建材领域、有价金属元素提取、化工领域、环保领域[5-6]。其中在建材方面的应用与从赤泥中提取有价金属是目前应用最为成熟的领域,对化工与环保领域的探索则还在进行中。在建材领域应用研究中,赤泥被用来与粉煤灰、砂石混合后制备墙砖[7],以及通过钙化-碳化法将赤泥中碱含量降低后用于部分替代生料制备水泥熟料[8],还有通过熔融法将赤泥与钢渣混合制成微晶玻璃[9]。但由于赤泥中Na元素的存在,制成的产品会出现“反碱”的现象[10],因此,赤泥在建材领域的应用应重点关注Na元素。高铁赤泥中含有多种有价元素,其中的氧化铁含量更是高达30%以上。因而,从赤泥中提取铁是现今研究的重点方向之一,在国内外学者对高铁赤泥的大量提铁研究中,目前主流的提铁方式共有3种:物理磁选、湿法浸出与火法熔炼[11-13]。物理磁选法操作方便,但是磁选效率低,铁精矿杂质较多,还会产生大量废渣[14]。湿法浸出则主要存在酸耗大、铁的溶出率低和酸浸液的分离提纯困难等问题[15]。火法冶金是目前应用较多的工艺,一般采用煤基还原回收铁,其中先焙烧后磁选法可获得铁精矿,金属收得率较高,操作简便;熔融还原法可直接得到生铁,可直接用于炼钢[16-17]。
基于以往研究,东北大学特殊冶金与过程工程研究所发明了“一种高铁赤泥提铁及直接水泥化的方法”[18],该方法是将赤泥与还原剂和造渣剂混料后进行涡流熔融还原,还原后的熔融渣经高温调质制备水泥熟料,以此实现赤泥规模化资源化消纳。赤泥还原渣的主要成分为铝氧化物、硅氧化物、钙氧化物以及钠铝硅钙共同形成的复杂物质,铝酸盐水泥熟料的化学成分主要为Al2O3、CaO、SiO2、TiO2和Fe2O3,两者的化学成分较为相似。铝酸盐水泥主体部分是以Al2O3和CaO结合而成的矿物,两者的质量分数占整个水泥熟料的80%以上。铝酸盐水泥的水化强度主要取决于水泥熟料的化学组成和关键物相的结晶情况,铝酸盐水泥的要求尽可能多地获得CA(CaO·Al2O3,),避免C2AS(CaO·2Al2O3·SiO2)和其他杂相的结晶。因此,研究还原过程中Al、Si元素的含量变化及物相转化的规律对水泥熟料的成分调配及铝酸盐水泥的水化强度有着直接性、关键性的意义。为了实现赤泥涡流熔融还原渣的水泥化,需对还原过程中赤泥中主要造渣元素Al、Si以及有害元素Na的行为进行研究[19]。本文通过对还原渣进行成分分析、物相分析以及微观形貌分析,研究了赤泥熔融还原过程中Na、Al、Si 3种元素的行为变化,为赤泥涡流熔融渣水泥化提供理论依据。
1 实验部分
1.1 实验原料
原料来自于我国某氧化铝厂的高铁赤泥,其化学成分由XRF测定,结果如表1所列。图1所示为原料XRD图谱,分析可得其主要组成为 Fe2O3、TiO2、CaO、NaAlSiO4、SiO2及CaSiO3。还原剂选用了秸秆炭,将计算出应添加的碳质量定为1,控制配碳比为1.1,使用秸秆炭182.5 g,其成分如表2所列。本实验在熔融还原过程中分别使用分析纯CaO 210 g和分析纯CaF2 8.4 g作为添加剂来调节渣的黏度和流动性。
表 1 高铁赤泥的化学成分组成Table 1. Chemical composition of high-iron red mud成分 Na2O Al2O3 SiO2 CaO TiO2 Fe2O3 烧失量 含量 3.02 22.50 6.29 2.03 6.05 55.30 4.81 表 2 秸秆炭的成分组成Table 2. Composition of stalk carbon组分 固定碳 挥发分 灰分 水分 含量 77.11 10.02 7.44 5.43 1.2 实验方法
先将经过粉碎后的赤泥烘干,待高铁赤泥干燥后,研磨过筛备用,取出1 000 g的高铁赤泥及根据配比已计算好的秸秆炭、CaO和CaF2,混合均匀。将混合原料放入石墨坩埚中,随中频感应炉一同升温,观察到加入的原料出现熔化时,记录此时的温度与现象,待物料完全融化后开始计时。开始时反应较为剧烈,故前3 min每隔1 min取样一次,进行水淬,分别记为样品1、样品2、样品3。待到实验进行到一半时间时,再取样一次后水淬,记为样品4。计时结束后,将坩埚内剩余物料倒出后水淬,得到水淬还原终渣,记为样品5。具体工艺流程如图2所示。
本实验使用的仪器主要有:如图3所示的中频感应炉SPZ-160、红外测温枪AR992、电热鼓风干燥箱101。
中频感应炉SPZ-160的最大输入功率为160 kW,振荡频率为1~20 kHz,输入电流为30~320 A,输出电压为70~550 V。
红外测温枪AR992的测温范围为200~2 500 ℃,精确度为±2 ℃,测量物距比为80∶1,使用9 V电池。
电热鼓风干燥箱101的温度设定范围为50~300 ℃,容积为43 L,温度波动率为±1 ℃,输入电压为220 V,功率为800 W。
本实验使用X射线荧光光谱仪和X射线衍射仪分别对实验原料、还原渣等样品进行化学成分分析和物相分析,所需要的检测仪器如表3所列。
表 3 检测仪器Table 3. Testing instrument检测项目 检测设备 型号 生产厂家 Na2O、Al2O3、SiO2、CaO和TFe的定量分析 X射线荧光光谱仪 ZSXPrimus Ⅱ 日本理学公司 样品物相分析 X射线衍射仪 D8 Advance 德国布鲁克公司 样品微观形貌分析 场发射扫描电镜 JSM-7800F 日本电子公司 2 实验结果与分析
2.1 Na元素在还原过程中的走向及矿物转型
Na元素在铝土矿中含量较低,一般不高于5%[20],在高温下经过拜耳法铝酸钠溶液的强碱溶出作用下,一部分的Na 会进入固相形成复杂、多样的铝硅酸钠矿物。根据图1所示的本实验所使用的高铁赤泥XRD图,可以看出Na元素在本研究使用的原料中多以NaAlSiO4形式存在。根据图4可以看出,Na元素含量在还原过程中不断降低, 这是因为在1 000 ℃时高铁赤泥中含有的Na2O发生了还原反应,产生了金属Na。而金属Na的沸点(883 ℃)很低,因此还原出的金属Na有一部分在高温下迅速挥发为钠蒸气。从图5中可以看出,物料中剩余的Na元素则与CaO和SiO2一同生成了NaCaAlSi2O7留在尾渣中。发生的主要反应如下:
Na2O+C=2Na+CO↑ (1) NaAlSiO4+CaO+SiO2= NaCaAlSi2O7 (2) 从图6中可以看出,在尾渣中大部分的Na元素存在于Al元素与Si元素的重合区域,剩余的则弥散分布在渣中,在此重合区域内,碳热还原反应后剩余的Na元素会与Al、Si、O、Ca元素结合,由图5可知,生成的物相主要为NaCaAlSi2O7。
2.2 Al元素在还原过程中的走向及矿物转型
氧化铝碳热还原反应温度一般为2 000~ 2 300 ℃[21],因而在本研究所设定的温度范围内(1 300~1 600 ℃),氧化铝很难与碳发生还原反应生成金属Al。从图1可以看出Al元素在原料中多以氧化物形态与Na和Si共同组成复杂氧化物(NaAlSiO4)。从图7可以看出,在反应整体进行的过程中,氧化铝的含量变化范围是29.5%~31.0%,并未发生较大的变化。图8说明了在反应进行过程中Al元素则主要以氧化铝的形式存在(3CaO·Al2O3);在反应后段有一部分的Al与剩余的Na元素形成NaCaAlSi2O7,另一部分则以Ca2Al(AlSi)O7的形式留在尾渣中。综上,反应过程中Al元素的含量并未发生较大变化,只是在过程中转化成3CaO·Al2O3、NaCaAlSi2O7和Ca2Al(AlSi)O7等几种不同的物相。发生的主要反应如下:
Al2O3+3CaO =3CaO·Al2O3 (3) Al2O3+2CaO+ 2SiO2 = Ca2Al(AlSi)O7 (4) Al元素与Si元素在元素周期表上处在相邻的位置,且均处于金属与非金属的分界线附近,因而虽然Al是金属元素,Si是非金属元素,但2种元素有着相似的化学性质。由图8可知,Al元素与Si元素会与Ca、O及少量的Na元素结合,形成复杂的氧化物组合相,如3CaO·Al2O3、NaCaAlSi2O7 和Ca2Al(AlSi)O7等。从图9中也可以看出,Al元素和Si元素在还原渣中的分布相似,存在大量的重合分布区域。
2.3 Si元素在还原过程中的走向及矿物转型
Si在原料中存在形式较复杂,多以SiO2、CaSiO3、NaAlSiO4、Fe2.454O4Si0.546等形式存在,由于氧化硅与碳的直接碳热还原温度较高,起始还原温度达到了1 600 ℃以上,这意味着在本实验反应进行中很难出现氧化硅被碳直接还原出硅单质的现象。从图10中可以看出氧化硅含量在整体反应进行中的变化范围是9.0%~9.5%,并未发生大范围的变化,只是存在不同物相间的转化。由图11可知,在反应进行的过程中,Si主要与CaO发生反应生成 Ca2SiO4,而在尾渣中Si则以NaCaAlSi2O7、Ca2Al(AlSi)O7、Ca2Al2SiO7、Ca3Al2Si3O12等形式存在。发生的主要反应如下:
Al2O3+2CaO+ 2SiO2 = Ca2Al(AlSi)O7 (4) SiO2+2CaO=Ca2SiO4 (5) Al2O3+2CaO+SiO2 = Ca2Al2SiO7 (6) Al2O3+3CaO+3SiO2 = Ca3Al2Si3O12 (7) 由图11可知,在Al元素与Si元素的重合区域内,Si以NaCaAlSi2O7、Ca2Al(AlSi)O7、Ca2Al2SiO7、Ca3Al2Si3O12等形式存在。从图12中可以看出,还原渣的微观形貌呈不规则的块状分布,渣中主要含有Ca、O、Al、Si元素,此外还有Ti和少量的Na元素。Ca元素和O元素弥散分布在整个渣中。由于实验过程中发生了氧化钠的碳热还原反应,因此Na元素的含量不高,均匀分布在整个渣中。
3 结 论
1) 由于在反应过程中发生了还原反应,Na元素含量不断降低,其余的Na元素则多从NaAlSiO4转化为了NaCaAlSi2O7留在了尾渣中。
2) Al元素在实验温度范围内很难发生碳热还原反应,因而含量没有太大变化,物相主要从3CaO·Al2O3转化成尾渣中的NaCaAlSi2O7、Ca2Al(AlSi)O7等。
3) 因为此反应最高温度也未能达到Si元素的碳热还原反应起始温度,所以Si元素在反应过程中含量变化不大,过程中主要由Ca2SiO4转化为NaCaAlSi2O7、Ca2Al(AlSi)O7 、Ca2Al2SiO7等。
4) 水淬还原渣的微观形貌呈不规则的块状分布,渣中的Ca、O以及少量的Na元素在整个渣中均匀弥散分布,Al和Si元素存在大量的高度重合区域,重合区域内存在着如3CaO·Al2O3、Ca2Al(AlSi)O7、Ca2Al2SiO7、Ca3Al2Si3O12等的Ca、O、Al、Si 4种元素的结合相。Ti元素则存在一些单独的区域与Ca、O元素结合。
赵中波 -
表 1 矿样表面离子吸附量及物理化学参数
Table 1 Ion adsorption capacity and physicochemical parameters of ore samples
矿样 离子吸附量(εX)/(mmol/g) 比表面积(S)/(m2/g) 阳离子交换容量(CEC)/(mmol/100 g) 铵蒙脱土 0.64 66.54 112.35 镧蒙脱土 0.30 表 2 不同阳离子浓度下矿样表面阳离子总吸附量Qx
Table 2 Total cation adsorption capacity Qx on ore samples surface at different ion concentrations
阳离子浓度/ (mmol/L) 0.001 0.004 0.008 0.010 0.040 0.080 0.100 0.400 0.800 铵蒙脱土上NH4+吸附量/(mmol/g) 0.635 0.635 0.634 0.634 0.635 0.635 0.635 0.635 0.636 镧蒙脱土上La3+吸附量/(mmol/g) 0.301 0.301 0.301 0.301 0.301 0.302 0.303 0.305 0.310 表 3 不同浓度下铵蒙脱土矿样表面上NH4+实际吸附量和计算吸附量
Table 3 Actual and calculated adsorption amount of NH4+ on the surface of ammonium-montmorillonite samples at different concentrations
本体溶液中NH4+浓度/ (mmol/L) 实际吸附量/(mmol/g) cPB 方程计算值/(mmol/g) sPB方程计算值(未考虑离子水化)/(mmol/g) sPB方程计算值(考虑离子水化)/(mmol/g) 0.001 0.64 10.96 1.12 1.12 0.004 0.64 10.93 1.12 1.12 0.008 0.63 10.92 1.12 1.12 0.01 0.63 10.89 1.10 1.12 0.04 0.64 10.89 1.11 1.11 0.08 0.64 10.84 1.11 1.11 0.10 0.64 10.82 1.11 1.11 0.40 0.64 10.57 1.10 1.10 0.80 0.64 10.39 1.10 1.11 表 4 不同浓度下镧蒙脱土矿样表面La3+实际吸附量和计算吸附量
Table 4 Actual and calculated adsorption capacity of La3+on the surface of lanthanum-montmorillonite samples at different concentrations
本体溶液中La3+浓度/ (mmol/L) 实际吸附量/(mmol/g) cPB方程计算值/(mmol/g) sPB方程计算值/(mmol/g) 0.001 0.30 3.53 0.38 0.004 0.30 3.54 0.37 0.008 0.30 3.53 0.37 0.01 0.30 3.53 0.38 0.04 0.30 3.44 0.38 0.08 0.30 3.42 0.38 0.10 0.30 3.40 0.38 0.40 0.31 3.16 0.36 0.80 0.31 3.02 0.36 -
[1] 吴昊, 邱廷省, 严华山. 降低离子型稀土原地浸矿残留稀土的对策探讨[J]. 有色金属科学与工程, 2021, 12(2): 105-112. [2] LI Z Y, LI J Y, XU R K, et al. Streaming potential method for characterizing the overlapping of diffuse layers of the electrical double layers between oppositely charged particles[J]. Colloids and Surfaces A: Physicochemical and Engineering Aspects,2015,478:22-29.
[3] XIAO Y F, GAO G H, HUANG L, et al. A discussion on the leaching process of the ion-adsorption type rare earth ore with the electrical double layer model[J]. Minerals Engineering,2018,120:35-43.
[4] XU Q H, SUN Y Y, YANG L F, et al. Leaching mechanism of ion-adsorption rare earth by mono valence cation electrolytes and the corresponding environmental impact[J]. Journal of Cleaner Production,2019,211:566-573.
[5] 谢芳芳,王观石,罗嗣海,等. 离子型稀土尾矿残留铵缓释性分析[J]. 中国环境科学,2021,41(9) 4333-4340. [6] 李贞. 高岭石表面上铵离子吸附的量子化学研究[J]. 有色金属科学与工程, 2017, 8(4): 36-41. [7] GOUY M. Sur la constitution de la charge électrique à la surface d'un électrolyte[J]. Journal de Physique Théorique et Appliquée,1913,25(148): 457-468.
[8] CHAPMAN D L. A contribution to the theory of electrocapillarity[J].The London,Edinburgh,and Dublin Philosophical Magazine and Journal of Science,1913,25(148):475-481.
[9] BARD A J,FAULKNER L R.Electrochemical methods:fundamentals and applications[M].New York:Wiley,1980.
[10] STERN O.Zur theorie der elektrolytischen doppelschicht[J].Zeitschrift Für Elektrochemie Und Angewandte Physikalische Chemie,1924,30(21/22):508-516.
[11] 尚飞,杨成松,张莲海,等.黏土颗粒扩散双电层影响因素分析[J].冰川冻土,2022,44(2):495-505. [12] CHANG F R C,SPOSITO G.The electrical double layer of a disk-shaped clay mineral particle:effects of electrolyte properties and surface charge density[J].Journal of Colloid and Interface Science,1996,178(2):555-564.
[13] CONWAY B E,BOCKRIS J O,AMMAR I A.The dielectric constant of the solution in the diffuse and Helmholtz double layers at a charged interface in aqueous solution[J].Transactions of the Faraday Society,1951,47:756-766.
[14] DOU W,CHEN M H,ZHOU S G.Fast iterative method for local steric Poisson-Boltzmann theories in biomolecular solvation[J].Computer Physics Communications,2023,291:108808.
[15] LIU B, LIU P, XU Z L, et al. Ionic size effects: generalized Boltzmann distributions, counterion stratification, and modified Debye length[J]. Nonlinearity,2013,26(10):2899-2922.
[16] LI B. Minimization of electrostatic free energy and the Poisson-Boltzmann equation for molecular solvation with implicit solvent[J]. Journal on Mathematical Analysis,2009, 40(6): 2536-2566.
[17] REINER E S,RADKE C J.Variational approach to the electrostatic free energy in charged colloidal suspensions:general theory for open systems[J].Journal of the Chemical Society,Faraday Transactions,1990,86(23):3901-3912.
[18] FOGOLARI F,BRIGGS J M.On the variational approach to Poisson-Boltzmann free energies[J].Chemical Physics Letters,1997,281(1/2/3):135-139.
[19] MOLINA F V. Soil colloids: properties and ion binding[M]. CRC Press, 2016.
[20] PLIMPTON S.Fast parallel algorithms for short-range molecular dynamics[J]. Journal of Computational Physics,1995,117(1):1-19.
[21] CYGAN R T,GREATHOUSE J A,KALINICHEV A G.Advances in clayff molecular simulation of layered and nanoporous materials and their aqueous interfaces[J].The Journal of Physical Chemistry C,2021,125(32):17573-17589.
[22] 中华人民共和国环境保护部.土壤 阳离子交换量的测定 三氯化六氨合钴浸提-分光光度法:HJ 889—2017[S].北京:中国环境科学出版社. [23] MARCUS Y.Thermodynamics of solvation of ions.Part 5.—Gibbs free energy of hydration at 298.15 K[J].J Chem Soc,Faraday Trans,1991,87(18):2995-2999.
[24] 李旺,祖波,李嘉雯.电解质对三峡库区泥沙Zeta电位及絮凝特性影响[J].科学技术与工程,2021,21(12):5071-5075. [25] PÁLINKÁS G,RADNAI T,SZÁSZ G I,et al.The structure of an aqueous ammonium chloride solution[J].The Journal of Chemical Physics,1981,74(6):3522-3526.
[26] 樊会敏,赵艳,刘新敏,等.非对称混合电解质体系中恒电荷颗粒表面与Stern电位的测定[J].土壤学报,2015,52(2):446-452. [27] ZHOU S G,WANG Z M,LI B.Mean-field description of ionic size effects with nonuniform ionic sizes:a numerical approach[J].Physical Review.E,Statistical,Nonlinear,and Soft Matter Physics,2011,84(2Pt 1):021901.
[28] NIGHTINGALE E RJr.Phenomenological theory of ion solvation.effective radii of hydrated ions[J].The Journal of Physical Chemistry,1959,63(9):1381-1387.
[29] 郭钟群, 周可凡, 金解放, 等. 流体理化特性对土体渗流规律影响研究进展[J]. 有色金属科学与工程, 2022, 13(4): 116-125. [30] 严奕博,韩鹏,谢瑞珍,等.pH对粉土电化学特性的影响试验研究[J].科学技术与工程,2017,17(17):297-303. [31] GOLDBERG S,SUAREZ D L.Influence of soil solution cation composition on boron adsorption by soils[J].Soil Science,2011,176(2):80-83.