稀土元素(包括14种稳定的镧系元素La、Ce、Pr、Nd、Sm、Eu、Gd、Tb、Dy、Ho、Er、Tm、Yb和Lu)及钇(Y),在环境中相对不移动,其在天然地下水、河水和海水中的浓度较低(例如Cravotta, 2008a)。相比之下,酸性矿井水(pH < 4.5)中的稀土元素浓度要高几个数量级:波兰Wi?niówka矿区的平均总稀土元素浓度为17,917 μg/L(Migaszewski等人,2019年);西班牙伊比利亚黄铁矿带的浓度为2,739 μg/L(Soyol-Erdene等人,2018年);美国伊利诺伊州煤矿排放物的浓度为1,091 μg/L(Lefticariu等人,2020年);宾夕法尼亚州烟煤和无烟煤矿排放物的浓度分别为285 μg/L和72 μg/L(Cravotta,2008a, 2008b)。通常,酸性矿井排水(AMD)中的SO4、Fe、Al、Mn等金属离子浓度也较高,可能对接收水体造成污染(Nordstrom, 2011;Cravotta和Brady, 2015)。然而,最近有研究将未经处理的AMD及其含有的金属固体视为稀土元素及相关关键矿物的替代国内来源(例如Stewart等人,2017;Vass等人,2019a, 2019b;Hedin等人,2020;Leon等人,2021;Hermassi等人,2022;Balaram等人,2023a, 2023b)。
理解稀土元素的水化学性质对于(1)预测其在环境中的存在形式和迁移路径以及(2)识别和优化从AMD及其他来源中高效回收稀土元素的处理技术至关重要。在常温条件下,稀土元素主要以三价阳离子形式存在,尽管铈可以同时存在于三价和四价状态,而铕可以存在于二价和三价状态(例如Brookins, 1989;Chakhmouradian和Wall, 2012;Williams-Jones等人,2012)。根据pH值和其他溶液性质,自由金属离子或未形成复合物的离子(如Me+3)会与氢氧化物、硫酸盐、碳酸盐、磷酸盐、氟化物、氯化物、草酸盐等配体结合,形成水合物或固体化合物(Brookins, 1989;Lee和Byrne, 1992;Millero, 1992;Liu和Byrne, 1997;Chung等人,1998;Prapaipong等人,1999;Luo和Byrne, 2001, 2004;Luo和Millero, 2004;Schijf和Byrne, 2004, 2021;Williams-Jones等人,2012;Giffaut等人,2014;Migdisov等人,2016)。同时,这些水合金属离子还会与不同组成和结晶度的“水合金属氧化物”形成表面复合物(即吸附),包括水合氧化铁(HFO;Dzombak和Morel, 1990)、水合氧化铝(HAO;Karamilidis和Dzombak, 2010)以及水合锰(III-IV)氧化物(HMO;Tonkin等人,2004)。这些Fe、Al和Mn的固体通常覆盖在河床沉积物上,并是AMD处理污泥的主要成分,稀土元素和金属离子在这些污泥中会被吸附或沉淀(例如Cravotta, 2008a, 2008b;Stewart等人,2017;Hedin等人,2024)。
地球化学物种模型为评估水相-固相相互作用提供了定量框架,并展示了pH值及相应溶液性质对稀土元素在AMD中吸附和去除的影响(Lozano等人,2019;Cravotta, 2021, 2022)。结合水相物种分布模型,扩散双层(DDL)表面复合模型考虑了“强”和“弱”吸附位点上的质子交换反应,用于描述二价(Me+2)和三价(Me+3)阳离子被HFO(Dzombak和Morel, 1990;Verplanck等人,2004;Liu等人,2017;Lozano等人,2020;Horne等人,2026)和HMO(Tonkin等人,2004;Pourret和Davranche, 2013)吸附的情况;而对于HAO,仅考虑了弱吸附位点(Karamilidis和Dzombak, 2010;Lozano等人,2019;Boothe-Lordon等人,2026)。以往涉及表面物种模型的研究有助于解释单一吸附剂系统(如HFO、HAO或HMO)中稀土元素的吸附行为,但未考虑吸附剂的混合物,且在大多数情况下忽略了可能涉及SO4的表面复合反应(Verplanck等人,2004;Pourret和Davranche, 2013;Liu等人,2017;Lozano等人,2019, 2020;Horne等人,2026)。由于SO4通常是AMD中的主要离子,涉及SO4的物种反应对于预测AMD系统中的稀土元素浓度至关重要。例如,先前的研究表明SO4会增加Cu、Zn(Swedlund和Webster, 2001)、Co、Pb和Cd(Swedlund等人,2003)在HFO(铁氢氧化物和施韦特曼石)上的吸附。Lozano等人(2019)还观察到SO4的存在会增加稀土元素在HAO(碱铝矿)上的吸附。由于缺乏用于模拟稀土元素水相-固相相互作用的数据,许多研究(如Verplanck等人,2004;Liu等人,2017;Ashby, 2017)使用了针对二价阳离子在HFO(Dzombak和Morel, 1990)、HAO(Karamilidis和Dzombak, 2010)或HMO(Tonkin等人,2004)吸附的线性自由能关系来估算稀土元素的平衡常数。然而,正如Pourret和Davranche(2013)所指出,并且在本研究中进一步说明的,使用为二价阳离子开发的LFER计算三价稀土元素的结合常数可能会导致不准确性。因此,尽管现有的表面复合模型有助于理解环境中微量元素的迁移行为,但由于缺乏稀土元素表面复合物的热力学数据以及吸附剂化学成分和物理性质的不确定性,这些模型的预测能力受到限制。
本文描述了实验室实验及相应的计算结果,这些计算得出了稀土元素与合成HFO和HAO之间表面复合反应的平衡常数(即结合常数),这些吸附剂代表了AMD系统中的水合金属氧化物。通过一系列台架实验研究了pH值和SO4对稀土元素在H2SO4或HCl介质中与Fe、Al或Mn共沉淀和/或吸附的潜在影响。实验数据通过参数估算软件PEST(Doherty, 2015)和地球化学建模软件PHREEQC(Parkhurst和Appelo, 2013)进行处理,以计算DDL模型中使用的表面复合反应的平衡常数(例如Dzombak和Morel, 1990)。新的稀土元素表面物种反应被纳入了一个综合的热力学数据库文件中,该文件可用于PHREEQC和PHREEQ-N-AMDTreat + REYs水质建模工具(Cravotta, 2021, 2022),这些工具能够预测环境样品中稀土元素的水相、吸附态和沉淀态浓度分布。