• 中文核心期刊
  • 中国科技核心期刊
  • CSCD来源期刊
  • DOAJ 收录
  • Scopus 收录

高温高压下元素配分的原位实验与计算模拟研究进展

黄圣轩, 巫翔, 秦善

黄圣轩, 巫翔, 秦善. 高温高压下元素配分的原位实验与计算模拟研究进展[J]. 岩矿测试, 2016, 35(2): 117-126. DOI: 10.15898/j.cnki.11-2131/td.2016.02.002
引用本文: 黄圣轩, 巫翔, 秦善. 高温高压下元素配分的原位实验与计算模拟研究进展[J]. 岩矿测试, 2016, 35(2): 117-126. DOI: 10.15898/j.cnki.11-2131/td.2016.02.002
Sheng-xuan HUANG, Xiang WU, Shan QIN. Research Progress on in situ Experimental and Theoretical Simulations of Element Partitioning under High Temperature and High Pressure[J]. Rock and Mineral Analysis, 2016, 35(2): 117-126. DOI: 10.15898/j.cnki.11-2131/td.2016.02.002
Citation: Sheng-xuan HUANG, Xiang WU, Shan QIN. Research Progress on in situ Experimental and Theoretical Simulations of Element Partitioning under High Temperature and High Pressure[J]. Rock and Mineral Analysis, 2016, 35(2): 117-126. DOI: 10.15898/j.cnki.11-2131/td.2016.02.002

高温高压下元素配分的原位实验与计算模拟研究进展

基金项目: 

国家自然科学基金资助项目(U1232204, 41473056)

详细信息
    作者简介:

    黄圣轩, 硕士研究生, 地球化学专业.E-mail:sxhuang@pku.edu.cn

    通讯作者:

    巫翔, 研究员, 主要从事实验地球化学研究.E-mail:xiang.wu@pku.edu.cn

  • 中图分类号: P599;P574.2;O657.34

Research Progress on in situ Experimental and Theoretical Simulations of Element Partitioning under High Temperature and High Pressure

  • 摘要: 地球内部元素的配分行为具有重要意义, 有助于揭示地球内部物质循环、元素成矿、超临界流体性质等基础地质问题。同步辐射微聚焦X射线荧光光谱(SR-μXRF)结合金刚石压腔(DAC)技术, 可以在高温高压条件下实现高准确度和高分辨率的原位(in situ)测试元素在不同物相间的浓度, 并获取配分系数。而基于地球化学热力学的理论计算模拟, 可以获得元素在不同物相间的赋存形式, 有助于揭示元素配分机制。本文综述了SR-μXRF结合DAC技术和地球化学热力学计算模拟原位测定高温高压下元素配分的研究方法, 及其在地球内部元素迁移、成矿作用和地球早期形成过程等领域中的应用进展, 认为利用该套方法研究高温高压下元素配分行为的体系将变得更加复杂, 对温度、压力条件则要求更高。本文旨在更新和丰富高温高压下元素配分系数数据, 深化对地球内部物质循环等问题的认知。
  • “成分”和“结构”是研究物质的两个最重要、最根本的维度,其决定了物质的物理、化学性质等。“温度”和“压力”作为两个重要的热力学参数,深刻地影响着物质的组成和结构。地球内部的高温高压环境使地球内部物质的组成和结构丰富多彩,同时也给地球内部物质组成和结构的研究带来了巨大挑战。自同步辐射X射线结合金刚石压腔(DAC)技术的原位(in situ)测试方法得到应用以来,地球内部物质结构的研究取得了突飞猛进的进展[1, 2, 3, 4]。但是,地球内部无时无刻不进行着复杂的物理化学变化、能量交换和物质循环,使得地球内部物质成分的研究更为复杂。

    在地球内部物质成分的研究过程中,学者们着重关注的是元素在不同物相之间的配分行为。地球内部物质的主要存在形式为固相(solid)、熔体相(melt)、流体相(fluid),尤其地壳至上地幔这个区域,是三相相互作用最为明显且重要的位置。在地壳至上地幔区域存在着大量的地质流体,Schmidt等[5]估算在截面积为1平方米、高度为大洋岩石圈厚度的空间范围内,平均约有106千克的水在俯冲带被带进地球内部。同时,流体相与固相、熔体相相比密度小、黏度小,迁移能力强,具有更高的反应活性,可以与围岩或熔体反应达到再平衡。利用元素在三相中的配分行为及基本的地质、物理、化学定律等,可以推算出在某一温度、压力、地质背景条件下元素的含量及其输运过程;也可以在已测得元素含量的情况下,定量地分析出在特定的温度、压力条件下,各物相的相对含量以及对形成某一地质作用的贡献程度[6, 7, 8]。因此,研究元素的配分行为可以有效地反演地质过程,理解区域内元素的迁移、物质的循环。其次,元素的配分行为、流体的成分以及环境的物理化学状态决定了流体富集元素的种类、富集元素的存在形式[9, 10]。研究元素配分行为也可以用于模拟元素富集时的地质环境、迁移形式、富集状态、成矿机制等;此外还有助于了解地球内部超临界流体的组成、性质。超临界流体(supercritical fluid,SCF)是指温度、压力超过临界点处于超临界状态的物质,具有极强的溶解能力、很好的热传导率和较高的反应活性。超临界流体的性质是密度的函数,而密度在超临界区间内又是温度、压力的函数,因此地球内部不同深度的流体有不同的性质,并对地球内部物质产生不同的影响[11, 12, 13, 14, 15]

    目前,元素配分行为研究的样品主要来源于:① 天然样品[16, 17, 18, 19];② 大压机或活塞圆筒等合成的样品[20, 21, 22, 23, 24, 25, 26]。利用电子探针技术(EMP)或电感耦合等离子体质谱法(ICP-MS)测试样品中元素在不同物相之间的浓度,进而确定其配分行为。但是,天然样品在分析之前受到了不同程度的物理、化学作用的改造,合成样品在分析之前会经历泄压、淬火的过程,这些都会给测试结果带来不确定性。同步辐射微聚焦X射线荧光光谱(SR-μXRF)结合DAC技术实现了原位测量元素在各物相之间的浓度[27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37]。首先,这种测试手段可以实现实时、原位直接测试待测元素在高温高压下的配分情况,便于分析非平衡态,有助于研究元素的扩散行为、配分机制,探索地球内部的动力学过程。in situ SR-μXRF可探测的最低浓度在几个μg/g(一般是5 μg/g)以下,具有很高的分辨率,可以满足微量元素的低浓度限制。此外,原位测试可以避免分析天然样品、泄压淬火样品的不确定性[33],测试结果具有很高的可信度。

    尽管in situ SR-μXRF结合DAC技术在元素配分行为的研究方面具有巨大的优势,但截至目前,DAC结合同步辐射微聚焦X射线谱学多用于测试物质在极端条件下的结构,在成分分析上还没有被很好地利用。现阶段利用该实验手段研究的体系较简单,体系的温度、压力条件不高,可靠的高温高压下元素在不同物相之间(尤其是含有流体相)的配分系数数据稀少,限制了对地球内部物质循环的问题的认知。此外,随着计算机技术的飞速发展,基于地球化学热力学的理论计算逐渐兴起,它既可以为实验结果提供热力学层面的支持,也可以为实验的开展提供新思路,但该方法刚刚兴起,学术界对其重视程度不够。本文在国内外最新研究成果基础上,阐述了SR-μXRF主要原理、DAC实验技术、理论计算模拟方法及其在高温高压下元素配分研究中的应用。

    同步辐射是带电粒子运动速度接近光速在电磁场中偏转时发出的一种电磁辐射。新的同步辐射光源的光斑尺寸很小(例如,上海同步辐射光源最小光斑尺寸可达1.6 μm×1.8 μm),可以实现微聚焦测试。同步辐射X射线与物质相互作用时,元素的原子通过内层电子的跃迁产生的荧光光谱是该元素的特征标志,可以定性判定待测物质的元素组成。同时,X射线荧光谱线的峰面积与元素的含量呈正相关。在没有元素含量已知的标准参考物质时,可以利用峰面积的比值分析待测物质元素的相对含量;在有元素含量已知的标准参考物质时,可以定量测定待测物质元素的绝对含量。

    DAC型号和种类繁多,但其基本构成和工作原理相同。DAC大体由三部分组成:金刚石压砧、支撑加压部分、外部机械装置。金刚石压砧由一对金刚石和密封垫组成。金刚石一般会粘在碳化钨座垫上,密封垫中间打一个小孔作为样品腔。根据不同的金刚石压腔,外部的支撑加压部分会有不同。有活塞-圆筒式加压、螺栓式加压等。外部机械装置是固定金刚石压砧和加压装置的部件。1993年Bassett等[38]基于Merrill-Basset型DAC发明了水热金刚石压腔(HDAC),它的大小是传统DAC的3~4倍。在原有的DAC基础上,HDAC附加了一套完整的加温、测温装置(图 1)。水热金刚石压腔的温度可以达到1200℃,压力达到10 GPa,可以模拟地壳至上地幔温压条件下的地质作用,尤其是适用于模拟流体与地质物质之间的相互作用[39, 40]

    图  1  HDAC腔外加温、测温装置及HDAC结合SR-μXRF实验光路示意图(修改自文献[29, 39])
    碳化钨座垫包裹有金属丝用以对上下两个金刚石压砧加热。与金属丝相连有可变电阻,可以保证上下两个金刚石压砧被加热到的温度是一致的(一般误差在5℃以下)。与金刚石相连有一组热电偶用以测量温度,由于其不在样品腔内,因此温度的测量不受压力的影响,准确度很高(一般测量精度在0.1℃)。探测器a表示传统测试方法,探测器与入射X射线成~20°;探测器b表示带有凹槽的HDAC测试方法,探测器与入射X射线成~90°,X射线荧光只来自于同一物相中的元素。
    Figure  1.  Schematic diagram of heaters, measuring of the temperature of the sample and the optical path in HDAC combined with SR-μXRF (Modified from Reference [29, 39])

    利用XRF对DAC内的样品进行测定的传统方法是使入射X射线与DAC上平面成~10°,而接收位置偏离入射X射线~20°(图 1)。但是,在产生X射线荧光的同时,入射X射线也会与物质产生瑞利散射作用和康普顿散射作用,其产生的背底会影响X射线荧光光谱的分析。Schmidt等[29]将水热金刚石压腔内的下压砧金刚石中间处形成一个凹槽(图 1),接收由此射出的X射线荧光。这种测试方法使散射产生的背景最小,得到的光谱的信噪比最好。同时,当金刚石压腔内存在多种物相时,此处接收的X射线荧光是同一物相中的元素产生的,不存在共激发现象,提高了测量的准确性。

    除实验模拟外,计算模拟也是研究科学问题的重要方法,尤其是现代计算机技术得到了飞速发展,使得计算模拟成为验证实验或预测新现象不可或缺的手段。应用于地学的计算模拟一般称为地球化学反应模拟。地球化学反应模拟建立在热力学平衡基础上,它作为一种新兴的研究方法,在研究地球化学反应体系中物质的存在形式、转化方式、迁移规律、反应过程、平衡状态等方面都发挥了重要的作用。国际上被广泛应用的地球化学模拟软件主要有EQ3/6[41]、PHREEQC[42]、CHESS[43]、MINEQL+[44]、GWB[45]等。

    地球内部物质组成元素种类繁多,将每一种元素的配分行为、迁移模式研究清楚既不现实,亦无价值。因此,学术界将主要精力集中在具有地球物质循环指示性的元素(如大离子亲石元素、高场强元素和稀土元素)以及对人类生活具有重要作用的稀有矿产元素(如铂族元素)等方面。此外,随着实验技术的不断发展,学者们开始尝试在更高的温度、压力条件下探索下地幔、核幔边界物质循环以及地球形成早期的物质分异的演化规律。

    Rb、Sr是用于指示岩浆、岩石成因的大离子亲石元素。两者的固相-熔体相分配系数Dis/m已被研究得非常清楚。但是在火成岩形成过程中,流体与熔体之间的反应是非常常见的。Dis/m已经不能完整地描述Rb、Sr在火成岩形成过程中的全部配分行为,因而其流体相-熔体相分配系数Dif/m以及温度、压力、组分等对Dif/m产生的影响备受关注。Bai等[23]利用EMP和ICP-MS分析淬火样品测定了Rb、Sr的Dif/m,研究发现在750℃、0.4 GPa条件下DRbf/m和DSrf/m都大于1。Borchert等[36]利用SR-μXRF和HDAC技术并对比淬火实验,系统地研究了Rb、Sr的Dif/m及其影响因素,研究表明在实验条件范围内,DRbf/m和DSrf/m都小于1,即在不同条件下Rb、Sr都倾向于进入熔体相(图 2)。在一定温度下,Dif/m随压力的增加而增大,并且在流体相中不含Cl-时,压力对Dif/m的影响更为显著;同时,Dif/m随流体相中Cl-浓度的增加而增大,尤其是在流体相中含有Cl-时,DRbf/m和DSrf/m较流体相的成分为H2O时增加了近两个数量级。之后Borchert等[46]通过X射线吸收谱和第一性原理计算发现(第一性原理计算是基于量子力学、密度泛函理论的计算模拟,不同于前文提到的基于热力学平衡的地球化学热力学计算模拟),在高温高压下,Sr在含有Cl-的流体中会形成[SrCl(H2O)6]+,增加了Sr在流体相中的稳定性,使得DSrf/m大幅度增大。

    图  2  不同流体相浓度条件下,Rb(a)和Sr(b)的流体相-熔体相分配系数随压力的变化图(修改自文献[36])
    Figure  2.  Fluid-melt partition coefficients of Rb (a) and Sr (b) as a function of pressure for various fluid compositions (Modified from Reference[36])

    但Borchert等仅对HDAC内的流体相进行了原位测试,对熔体相仍然测试的是淬火样品。随后,Kawamoto等[37]改进了Borchert等的实验,对HDAC内的流体相和熔体相都进行原位测试,研究发现DRbf/m约比DSrf/m大了近一个数量级,即含Cl-流体运移Rb的能力比运移Sr的能力要大得多。这意味着在流体中含有一定量的Cl-时,Rb/Sr值将发生很大变化,其指示的地质过程界限也需要重新厘定。此外,对比Borchert等的原位测试结果和淬火实验结果可以发现,当流体中含有Cl-时,两种实验手段获得的分配系数相近,而当流体中不含Cl-时,两种实验手段获得的分配系数差别很大。Borchert等[33]的研究表明Rb、Sr在流体相与熔体相之间达到平衡是非常快速的,在淬火过程中会发生逆反应,实现Rb、Sr在流体相与熔体相之间的再平衡,由此分配系数产生偏差。而Rb、Sr在含有Cl-的流体中会形成稳定的配合物,可以阻止逆向反应的持续发生,因而两种实验方法得到的结果相吻合。这里也表明了淬火实验的局限性和原位测试技术的优越性。

    高场强元素也是研究岩石成因、区分岩浆类型的有效示踪剂。例如,与大洋中脊玄武岩相比,弧岩浆具有明显的富集大离子亲石元素、亏损高场强元素的地球化学特征[47]。在俯冲板片向下俯冲的过程中,高温高压条件下形成了金红石等稳定的矿物相,高场强元素如Nb、Ta等以类质同象形式进入金红石晶格脱离俯冲板片。同时,金红石的高稳定性造成寄宿于金红石中的高场强元素的活动性也非常低。Tropper等[48]通过实验证明了金红石在纯水中的低溶解度(约150 μg/g),间接证明了以上推断。然而,俯冲板片俯冲脱水产生的流体并不是纯水。在高温高压下,这些流体中含有卤素(Cl、F),碱金属(Na、K)和硅铝质物质等[11]。Manning等[32]利用SR-μXRF和HDAC技术研究发现,金红石的溶解量与流体中NaAlSi3O8和Na2Si3O7的溶解量存在正相关性,随着流体中铝硅酸盐溶解量的增加,Ti的溶解量可由~5 μg/g增至~4000 μg/g,尤其是Na2Si3O7对金红石溶解的促进作用更明显。Manning等认为在高温高压条件下,流体中存在的碱性离子(如Na+、K+等)可以与Ti4+构成配合物,增加了金红石的溶解度。在NaAlSi3O8中,Na+/Al3+=1,Al3+全部占据四面体中心位置,为造网元素。而不含Al3+的硅酸盐(硅含量不变)的聚合度更低,便于Ti4+进入。Antignano等[49]利用失重法测定了金红石的溶解量,在实验条件范围内Ti的溶解量达~10 μg/g至~1000 μg/g。同时,Antignano等研究表明Na+/Al3+正相关于压力。因此,这可能意味着金红石的溶解度也将正相关于压力,但这需要更多的实验来验证。这些研究结果为低硅埃达克岩的富Ti结果提供了另一种机制:传统观念认为低硅埃达克岩富集Ti是地幔楔部分熔融造成的[50],而Manning等的研究结果表明富含卤素和铝硅酸盐的流体也可以迁移一定量的Ti,导致富Ti的埃达克岩的形成。

    同时,Tanis等[51]利用SR-μXRF和HDAC技术主要模拟了冷俯冲温压条件下,金红石中类质同象存在的Nb在流体中的迁移行为,研究发现,在含有氯化钠的流体中,Nb的溶解量正相关于温度,而与压力、流体密度呈负相关性;Nb在含有氟化钠的流体中的溶解量更多。Tanis等认为Nb、Ti除了可以与碱金属离子形成配合物,还可以与流体中的卤素形成配合物,提高高场强元素在流体中的稳定性。这意味着在俯冲板片俯冲脱水过程中形成的金红石虽然是稳定的矿物相,但是地球内部地质流体仍然可以输运部分Ti以及寄宿于金红石内的高场强元素。地球内部地质流体输运高场强元素的能力强烈依赖于流体的成分:流体中溶解的铝硅酸盐、卤素等都可以提高流体输运高场强元素的能力。

    以上研究打破了金红石及寄宿于金红石内的高场强元素具有极低活动性的传统观念,在具有一定量的流体参与的岩浆演化过程是更为复杂的。

    稀土元素(REE)的离子半径相近,在各种地质作用中总是同时存在。稀土元素离子的价态一般为+3,Eu、Ce为可变价元素,可分别呈现+2、+4价。离子价态、离子半径的改变使Eu、Ce的地球化学行为发生改变,在稀土元素丰度图谱中表现为异常,这些异常被用来指示特殊的地球化学作用。Allen等[52]研究了热液流体对橄榄岩的交代作用及其对橄榄岩中稀土元素的影响。在400℃、50 MPa下,含有Cl-的热液流体与橄榄岩反应之后溶解的阳离子(如Na+、Ca2+等)显著增加并出现了Eu正异常,而在低Cl-热液流体中没有该异常。热液流体交代辉石,使可以与Eu2+类质同象的Ca2+进入流体,Eu2+也进入流体。GWB软件模拟结果(图 3)表明该温压条件、弱酸性环境中Eu2+主要以EuCl42-形式存在,其他稀土元素氯离子配合物不及EuCl42-稳定,这成为热液中Eu2+正异常的主导因素。而Loges等[53]在德国黑森林的克莱拉萤石矿里的锰氧化物中发现Ce强烈负异常,这与传统的Ce4+活动性更低的结果相悖。GWB软件模拟排除了次生锰氧化物矿物中强烈的Ce负异常来源于热液流体中稀土元素丰度模式的异常。Loges等通过分析样品认为Ce4+的活动性比其他+3价稀土元素的活动性更高,而不是更低。在次生矿物生成、REE再分配的过程中,高活动性的Ce4+倾向于保存在热液流体中,其他稀土元素则倾向于进入次生矿物相。Ce4+的高活动性是由于在热液流体中形成了稳定的配合物。Ce4+不仅可以与Cl-、F-、CO32-等形成配合物,还可以与有机配体形成稳定的配合物,极大地增加了其在流体中的溶解度[54]。稀土元素的活动性以及在不同物相之间的配分行为已经研究了几十年,但仍有诸多影响稀土元素的因素没有被挖掘或是认识不够深入。而理论计算模拟可以在热力学和动力学上提供证据或是思路,以使这方面的研究更加完善。

    图  3  GWB软件计算模拟的400℃、50 MPa、pH=5.2条件下Eu在不同Cl-浓度的热液流体中的存在形式及其浓度(修改自文献[52])
    Figure  3.  Results of fluid speciation calculation at 400℃ and 50 MPa using the GWB software. The fluid's pH was fixed at 5.2 (Modified from Reference [52])

    在20世纪70年代,μg/g量级以下的物质浓度的测量成为现实。铂族元素(PGE)也因此可以作为岩浆、岩浆矿床成因的示踪剂,用以研究硫化物与岩浆的分离时间、物质平衡、贵金属矿床等问题[55, 56]。Augé等[57]研究发现在斑岩型Cu-Au矿床中,Pd、Pt的含量异常高。在其他一些类型的热液矿床中也有类似的发现[58, 59]。这些研究说明热液流体可以迁移岩石中的Pd、Pt,对PGE的地球化学示踪具有明显影响。Barnes等[60]利用GWB软件模拟了Pd、Pt在寄主岩石和热液流体之间的配分行为。计算模拟发现,在300℃、氧化条件下,热液流体中Pd、Pt的主要存在形式为PdCl42-、PtCl42-;在300℃、较还原条件下,热液流体中Pd、Pt的主要存在形式为Pd(HS)2(aq)、Pt(HS)2(aq);在弱氧化、弱还原条件下,两者在热力学上都变得不稳定,Pd、Pt的溶解度达到极小值。酸性热液流体中Pd的溶解度达到几十个ng/g,可以使寄主岩石中的Pd显著改变。在赤铁矿-磁铁矿-菱镁矿(HM)和磁铁矿-黄铁矿-磁黄铁矿-镍黄铁矿(MPP)两个体系中,热液流体中Pd、Pt的溶解量差异很大(表 1)。在HM体系中,Pd、Pt的溶解量为10-8~10-13 ng/g,可以认为热液流体几乎无法从岩石中迁移Pd、Pt;在MPP体系中,Pd、Pt的溶解量在ng/g量级,热液中Pd、Pt的主要存在形式是Pd(HS)2(aq)、Pt(HS)2(aq)。热液流体中硫的含量对Pd、Pt的溶解、迁移至关重要。之后,Mei等[61]利用分子动力学和X射线吸收谱佐证了这一点。同时,计算模拟的结果也表明可以利用原始矿床中Pd、Pt含量的减少量来指示高硫热液及反演固-液相相互作用。

    表  1  GWB软件计算模拟的不同体系下Pd、Pt在热液流体中的溶解度(修改自文献[60])
    Table  1.  Palladium and platinum solubility in hydrothermal fluids buffered by two mineral assemblages simulated by the GWB software (Modified from Reference [60])
    因素赤铁矿-磁铁矿-菱镁矿(HM)体系磁铁矿-黄铁矿-磁黄铁矿-镍黄铁矿(MPP)体系
    初始流体成分3 mol/kg氯化钠;10-5mol/kg稀盐酸3 mol/kg氯化钠;10-5 mol/kg稀盐酸
    固相体系成分赤铁矿-磁铁矿-菱镁矿体系磁铁矿-黄铁矿-磁黄铁矿-镍黄铁矿体系
    Pd终态浓度(ng/g)3.2×10-84.7
    Pt终态浓度(ng/g)9.3×10-130.8
    硫化氢终态溶解量(mol/kg)00.0007
    终态酸碱度(pH)6.65.8
    终态氧逸度(logfO2)-31-34
    流体中Pd、Pt主要存在形式PdCl42-;PtCl42-Pd(HS)2(aq);Pt(HS)2(aq)
    下载: 导出CSV 
    | 显示表格

    地球科学家通过高温高压实验、计算机模拟等手段对地球形成的早期历史进行了大量的研究[62, 63]。Snow等[64]研究发现地幔和地壳中的高度亲铁元素(HSE,如Pd、Pt、Ru等)的浓度比通过分配系数计算得到的浓度要高。后增薄层模型认为核幔分异之后地球受到了陨石撞击,造成地幔中亲铁元素的浓度异常[65]。Righter等[66]研究发现HSE的亲和性随压力的增加而减小,前人实验测得的分配系数存在一定的局限性。为了进一步探索地球形成早期的演化历史和更深处的物理化学过程,Petitgirard等[67]将激光加温金刚石压腔(LHDAC)与X射线衍射(XRD)、X射线荧光光谱(XRF)技术结合起来,实现了在更高温度、压力条件下原位研究元素的配分行为,用于探索地球更深处的物理化学过程。Petitgirard等在激光加温的同时对玻璃态硅酸盐-铁镍合金进行XRF采谱,获得了在不同时间下样品中元素浓度的变化(图 4)。在玻璃态硅酸盐发生熔融时(铁镍合金已经熔融),硅酸盐中的Ru、Pd进入铁镍熔体中,而Zr则保留在硅酸盐熔体中,显示了Ru、Pd的亲铁性和Zr的相对亲石性。由于条件所限,Petitgirard等没有提取这些元素的分配系数,但是这套装置为研究亲铁元素在地核、地幔之间的配分行为以及元素的扩散机制,进而探索地球形成的早期历史提供了强有力的保障。此外,在实际解决这类问题时,样品中一般不会存在流体,因此可以将样品在高温高压环境中淬火,利用EMP、ICP-MS等手段辅助XRF测试元素的浓度。该种实验装置和研究策略最近几年才逐步发展起来,这方面的成果还没有被广泛报道。以欧洲同步辐射光源(ESRF)和德国同步辐射光源(DESY)为依托的一些欧洲学者已经展开了深入的研究,在该领域走在世界前列[68, 69]

    图  4  铁镍熔体中Ru、Pd、Zr元素的X射线荧光强度随时间变化图(修改自文献[67])
    Figure  4.  In situ evolution of the XRF intensity of Ru, Pd, and Zr in the metallic foil as a function of time (Modified from Reference [67])

    地球内部物质的成分和结构研究是地球科学中两个永恒的主题。物质结构研究技术发展已相对成熟、完整,但其成分研究技术尚不完善。DAC结合SR-μXRF实验方法和理论计算模拟,为原位研究地球内部物质的成分打开了一扇新的大门,并且已经在研究地球内部固相、熔体相、流体相三相相互作用,元素赋存方式、迁移规律、成矿作用等方面取得了一定的进展。最近的高温高压下元素配分行为的实验研究充分表明,地球深部富含卤素、铝硅酸盐的流体相可以溶解、运移数量可观的大离子亲石元素以及被认为活动性非常低的高场强元素、稀土元素,并且随着温度压力的增加,这种作用也越来越显著。这些结果更新了我们原有的对一些岩石成因、物质循环模式的认知。

    但是目前研究的体系相对单一,大多数研究都是着重探讨流体中卤素对元素配分行为的影响,还有一些重要的因素没有考虑,如S元素、酸碱度、氧逸度等。同时,多数研究的熔体相都是硅酸盐熔体,熔体相为碳酸盐的研究还没有被广泛报道。其次,目前研究体系的压力较低,一般不超过2 GPa,这对于研究俯冲带地区的元素迁移、物质循环是不够的。尽管近年来对流体包裹体的研究不断深入,一些元素在气相中被运移的报道屡见不鲜,而目前还没有利用前述实验方法进行相关研究。另外,与第一性原理计算模拟不同,基于热力学平衡的地球化学计算模拟需要确定、合理的数据库,由于高温高压下元素在不同物相之间的配分系数数据稀少,这也限制了计算模拟发挥其巨大功能。以上诸多不足使得距离建立一套完整的元素配分理论、元素示踪模型、元素循环模式还很远。因此,未来该领域的研究体系将不断复杂,考虑因素更加全面;温度压力条件将不断升高,以真实反映地球内部高温高压下元素循环状态;通过实验建立准确的高温高压条件下元素配分系数数据库,使实验与计算模拟相辅相成,不断深化人类对于地球内部元素循环、物质组成的认知。

  • 图  1   HDAC腔外加温、测温装置及HDAC结合SR-μXRF实验光路示意图(修改自文献[29, 39])

    碳化钨座垫包裹有金属丝用以对上下两个金刚石压砧加热。与金属丝相连有可变电阻,可以保证上下两个金刚石压砧被加热到的温度是一致的(一般误差在5℃以下)。与金刚石相连有一组热电偶用以测量温度,由于其不在样品腔内,因此温度的测量不受压力的影响,准确度很高(一般测量精度在0.1℃)。探测器a表示传统测试方法,探测器与入射X射线成~20°;探测器b表示带有凹槽的HDAC测试方法,探测器与入射X射线成~90°,X射线荧光只来自于同一物相中的元素。

    Figure  1.   Schematic diagram of heaters, measuring of the temperature of the sample and the optical path in HDAC combined with SR-μXRF (Modified from Reference [29, 39])

    图  2   不同流体相浓度条件下,Rb(a)和Sr(b)的流体相-熔体相分配系数随压力的变化图(修改自文献[36])

    Figure  2.   Fluid-melt partition coefficients of Rb (a) and Sr (b) as a function of pressure for various fluid compositions (Modified from Reference[36])

    图  3   GWB软件计算模拟的400℃、50 MPa、pH=5.2条件下Eu在不同Cl-浓度的热液流体中的存在形式及其浓度(修改自文献[52])

    Figure  3.   Results of fluid speciation calculation at 400℃ and 50 MPa using the GWB software. The fluid's pH was fixed at 5.2 (Modified from Reference [52])

    图  4   铁镍熔体中Ru、Pd、Zr元素的X射线荧光强度随时间变化图(修改自文献[67])

    Figure  4.   In situ evolution of the XRF intensity of Ru, Pd, and Zr in the metallic foil as a function of time (Modified from Reference [67])

    表  1   GWB软件计算模拟的不同体系下Pd、Pt在热液流体中的溶解度(修改自文献[60])

    Table  1   Palladium and platinum solubility in hydrothermal fluids buffered by two mineral assemblages simulated by the GWB software (Modified from Reference [60])

    因素赤铁矿-磁铁矿-菱镁矿(HM)体系磁铁矿-黄铁矿-磁黄铁矿-镍黄铁矿(MPP)体系
    初始流体成分3 mol/kg氯化钠;10-5mol/kg稀盐酸3 mol/kg氯化钠;10-5 mol/kg稀盐酸
    固相体系成分赤铁矿-磁铁矿-菱镁矿体系磁铁矿-黄铁矿-磁黄铁矿-镍黄铁矿体系
    Pd终态浓度(ng/g)3.2×10-84.7
    Pt终态浓度(ng/g)9.3×10-130.8
    硫化氢终态溶解量(mol/kg)00.0007
    终态酸碱度(pH)6.65.8
    终态氧逸度(logfO2)-31-34
    流体中Pd、Pt主要存在形式PdCl42-;PtCl42-Pd(HS)2(aq);Pt(HS)2(aq)
    下载: 导出CSV
  • doi: 10.1126/science.1095932
    doi: 10.1126/science.1250274
    doi: 10.1016/j.epsl.2014.01.019
    doi: 10.1016/j.epsl.2015.02.008
    doi: 10.1029/97JB00788
    doi: 10.2113/gsecongeo.84.2.291
    doi: 10.1038/nature02644
    doi: 10.1038/nature03971
    doi: 10.1093/petrology/egm073
    doi: 10.1073/pnas.1010968108
    doi: 10.1016/S0016-7037(02)01049-9
    doi: 10.1007/s00126-006-0057-7
    doi: 10.1016/j.gca.2008.01.034
    doi: 10.2113/gsecongeo.103.5.877
    doi: 10.1038/348531a0
    doi: 10.1016/0016-7037(95)00244-T
    doi: 10.1127/ejm/9/3/0569
    doi: 10.1127/ejm/11/3/0415
    doi: 10.1016/j.gca.2004.11.005
    doi: 10.2138/am.2007.2066
    doi: 10.1016/S0168-583X(03)01100-5
    doi: 10.2138/am-2003-0705
    doi: 10.2138/am-2003-2-305
    doi: 10.1080/08957950701384572
    doi: 10.1016/j.gca.2010.04.001
    doi: 10.1016/j.gca.2009.10.033
    doi: 10.1186/1880-5981-66-1
    doi: 10.1063/1.1143931

    Parkhurst D L, Appelo C A J. User's Guide to PHREEQC (Version 2): A Computer Program for Speciation, Batch-reaction, One-dimensional Transport, and Inverse Geochemical Calculations, US Geological Survey Water-Resources Investigations Report 99-4259[R]. Colorado, Denver, 1999.

    Van der Lee J. Thermodynamic and Mathematical Concepts of CHESS, Technical Report Nr. LHM/RD/98/39[R]. Fontainebleau, France: Ecole des Mines de Paris, 1998.

    doi: 10.1016/0198-9715(92)90053-T
    doi: 10.1016/j.gca.2014.07.010
    doi: 10.1098/rsta.1991.0054
    doi: 10.2138/am.2005.1806
    doi: 10.2138/am-2015-5031
    doi: 10.1016/j.gca.2004.07.016
    doi: 10.1016/j.gca.2012.03.017
    doi: 10.1016/j.gca.2010.07.010
    doi: 10.2113/gscanmin.43.4.1355
    doi: 10.2113/gsecongeo.89.5.1053
    doi: 10.2113/econgeo.105.6.1047
    doi: 10.1016/j.oregeorev.2011.08.004
    doi: 10.1016/j.gca.2015.04.009
    doi: 10.1038/nature08477
    doi: 10.1038/34396
    doi: 10.1038/35019050
    doi: 10.1038/ngeo180
    doi: 10.1063/1.3680573

    Petitgirard S, Borchert M, Appel K, et al. Partitioning of HSEs between Silicate and Iron Melts at Condition Prevailing for the Early Earth[R]. DESY Photon Science Annual Report-User Experiment Reports, 2012.

    Borchert M, Petitgirard S, Appel K, et al. Sulphur and Selenium Partitioning between Silicate and Metallic Melts at Core-mantle Boundary Conditions[R]. DESY Photon Science Annual Report-User Experiment Reports, 2012.

  • 期刊类型引用(1)

    1. 巫翔,高春晓,王超. 高温高压装置研制和技术创新的发展现状与趋势. 地球科学. 2022(08): 2757-2764 . 百度学术

    其他类型引用(1)

图(4)  /  表(1)
计量
  • 文章访问数:  1621
  • HTML全文浏览量:  289
  • PDF下载量:  16
  • 被引次数: 2
出版历程
  • 收稿日期:  2015-12-15
  • 修回日期:  2016-02-19
  • 录用日期:  2016-03-04
  • 发布日期:  2016-02-24

目录

/

返回文章
返回