四种蛇纹石化程度定量分析方法的系统比较与适用性评估

《Chemical Geology》：Comparative evaluation of analytical methods for quantifying serpentinization extent

　　

蛇纹石化作用是指超镁铁质岩石（以橄榄石为主要矿物）与水之间发生的水岩反应，因其在从地表到地球内部各圈层挥发分循环中的广泛产出及关键作用而备受关注。蛇纹石化程度是研究蛇纹石化动力学的重要定量指标。然而，目前用于量化蛇纹石化程度的分析方法差异较大，包括重量法、热重（TG）分析、X射线衍射（XRD）分析、拉曼光谱、傅里叶变换红外光谱（FTIR）以及扫描电子显微镜（SEM）图像分析等。由于缺乏可用于实验室间数据比对的蛇纹石标样，难以评估已报道的蛇纹石化速率差异究竟反映真实的实验差异还是方法学假象。本研究系统比较并评估了四种方法——XRD、TG分析、拉曼光谱和FTIR光谱——测定蛇纹石化程度的稳健性。为此，研究人员制备了蛇纹石化反应涉及的标样矿物（橄榄石、利蛇纹石和氢氧镁石）及天然橄榄岩，并将其混合配制成三个二元标样混合物系列：橄榄石–利蛇纹石（Ol–Lz）、橄榄岩–利蛇纹石（Prd–Lz）和橄榄岩–氢氧镁石（Prd–Brc）。由于蛇纹石和氢氧镁石在天然和实验产物中常共存，研究人员进一步制备了利蛇纹石/氢氧镁石质量比固定为4.75:1（即反应式（1）的化学计量产物比例）的三元橄榄石–利蛇纹石–氢氧镁石（Ol–Lz–Brc）系列。采用全部四种方法为各系列建立校准曲线，并定量评估其精密度和准确度。随后将这些方法应用于实验室热液实验的蛇纹石化产物，以评估其在实际条件下的表现。研究结果明确了每种方法的优势、局限性和适用场景，为根据特定矿物组合和反应程度选择方法提供了实用指南。

蛇纹石化是超镁铁质岩石（以橄榄石为主要矿物）与水之间的水岩反应，因广泛产出并驱动地表至地球内部各圈层间的挥发分循环而备受关注（Evans et al., 2013; O'Hanley, 1996）。蛇纹岩主要由蛇纹石族矿物（纤蛇纹石、利蛇纹石和叶蛇纹石）组成，是洋壳和陆壳中H₂O及流体活动性元素的重要储库。蛇纹石俯冲可将C、N、H₂O、B、Cl等挥发分输送至地球深部，对地幔动力学过程至关重要。蛇纹石化通过形成碳酸盐矿物促进碳捕集与封存，有助于缓解大气CO₂升高驱动的全球变暖。蛇纹石化过程中产生的H₂、CH₄、NH₃等还原性气体不仅对可持续能源应用具有意义，也是早期地球生命起源与演化的重要驱动力，并支撑极端条件下的独特微生物生态系统。蛇纹石化程度是研究蛇纹石化动力学的关键定量指标。橄榄石蛇纹石化的主要产物为蛇纹石（Mg₃Si₂O₅(OH)₄）和氢氧镁石（Mg(OH)₂），其相对比例随温度、流体化学性质和时间变化。本文中蛇纹石化程度指蛇纹石在全岩固体中的质量分数（wt%），需与橄榄石转化程度区分：由于镁橄榄石完全蛇纹石化生成蛇纹石与氢氧镁石的质量比为4.75:1，完全反应样品含约82.7 wt%蛇纹石而非100 wt%。然而，动力学研究中采用的分析方法差异很大，包括重量法、热重（TG）分析、X射线衍射（XRD）分析、拉曼光谱、傅里叶变换红外光谱（FTIR）和扫描电子显微镜（SEM）图像分析。由于缺乏蛇纹石标样用于实验室间数据比对，难以评估已报道的蛇纹石化速率差异反映真实实验差异还是方法学假象。早期热液实验通常采用重量法监测蛇纹石化程度，但该方法仅能分辨总体重量差异而非各相，不适用于含多种含水矿物的组合。TG分析可定量矿物组合中的分子水和结构羟基，已证明对仅含单一含水矿物的体系有效，但当蛇纹石和氢氧镁石共存时，二者脱水温度范围重叠可能导致蛇纹石含量高估。XRD方法广泛用于监测金刚石压腔实验中蛇纹石化反应的时间演化，结合Rietveld方法可定量蛇纹石化程度，但热液蛇纹石产物的XRD峰强度受结晶度强烈影响，需进一步验证。拉曼光谱靶向蛇纹石3600–3720 cm^?1区域的羟基振动，Marcaillou et al.（2011）基于橄榄岩与天然蛇纹石的机械混合物建立了蛇纹石含量与积分峰面积的校准曲线，但当存在其他含水矿物（如氢氧镁石）时，其重叠的–OH振动峰难以区分。FTIR原理与拉曼光谱相似，同样受多种含水矿物共存时峰干扰的限制。SEM利用背散射电子图像测量矿物面积比例，可检测微量矿物，但人工选择的主观性可能引入系统偏差。综上，目前常用的四种方法（TG分析、XRD、拉曼光谱和FTIR光谱）各具不同的操作原理、校准方法和分析局限性，尽管在蛇纹石化研究中广泛采用，但尚未在相同实验条件下对这四种方法进行系统交叉比较。以往比较仅限于两两评估，通常涉及这四种技术之一与未纳入本研究的方法（如BSE图像分析或重力分离）。因此，已报道蛇纹石化速率的差异究竟反映真实实验差异还是方法学假象仍不清楚，亟需以同一套实验产物为基准对这些方法进行全面评估，以确定其相对准确度、精密度和适用范围。本研究系统比较和评估了XRD、TG分析、拉曼光谱和FTIR光谱四种方法测定蛇纹石化程度的稳健性。为此，研究人员制备了蛇纹石化涉及的标样矿物（橄榄石、利蛇纹石和氢氧镁石）及天然橄榄岩，混合配制成三个二元标样混合物系列：Ol–Lz、Prd–Lz和Prd–Brc。由于蛇纹石和氢氧镁石在天然和实验产物中常共存，进一步制备了利蛇纹石/氢氧镁石质量比固定为4.75:1的三元Ol–Lz–Brc系列。采用全部四种方法为各系列建立校准曲线，定量评估其精密度和准确度，随后应用于实验室热液实验的蛇纹石化产物以评估实际条件下的表现。研究结果明确了每种方法的优势、局限性和适用场景，为根据特定矿物组合和反应程度选择方法提供了实用指南。

研究人员首先制备了三种纯矿物（橄榄石、利蛇纹石和氢氧镁石）和一种天然橄榄岩作为标样材料，其组成通过X射线荧光（XRF）分析确定。橄榄石采自中国吉林省敦化，利蛇纹石采自山东泰安，氢氧镁石为Macklin氢氧化镁粉末，橄榄岩采自汉诺坝地幔捕虏体。将橄榄石和利蛇纹石粉末按约0–100 wt%的质量比混合，在玛瑙研钵中研磨10分钟，制成粒径5–50 μm的Ol–Lz标样混合物系列，同样程序用于制备Prd–Lz和Prd–Brc系列。三元Ol–Lz–Brc系列中利蛇纹石/氢氧镁石质量比固定为4.75:1，对应镁橄榄石理想蛇纹石化反应产物的化学计量比例。另制备了利蛇纹石–氢氧镁石质量比1:1的混合物，仅用XRD和TG分析。所有混合物均用XRD、TG、拉曼和FTIR分析以建立校准曲线。

XRD分析使用Bruker D8 ADVANCE X射线衍射仪，Cu–Kα辐射，45 kV和200 mA，2θ范围5–70°，步长0.02°，扫描速度2°/min。粉末样品采用常规前装样技术制备，每个样品重复测量四次以评估重现性。XRD数据用Jade软件处理，采用Rietveld方法、峰高比（PHR）方法和峰面积比（PAR）方法三种途径生成校准曲线。Rietveld精修使用Jade软件定量各标样混合物中蛇纹石和/或氢氧镁石的比例，参考标样包括夏威夷Salt Lake Crater橄榄石、马其顿Radusa Mine利蛇纹石和氢氧镁石。PHR和PAR方法测量橄榄石、利蛇纹石和氢氧镁石特征衍射峰的峰高和积分面积，选取橄榄石（2θ = 17.34°）的(020)峰、利蛇纹石（2θ = 12.22°）的(001)峰和氢氧镁石（2θ = 18.58°）的(001)峰用于校准。PHR方法将峰高比与矿物比例关联，PAR方法将峰面积比与矿物比例关联，均采用式（2）拟合，其中A为非线性回归确定的比例常数，W_max为目标矿物在相关系列中的最大质量分数。

TG分析使用METTLER TOLEDO TGA2热重分析仪，约15 mg样品粉末置于铂坩埚中，在N₂气氛（100 mL/min）下防止Fe氧化，50 °C平衡后以30 °C/min升温至1000 °C，每个样品三次重复测量，称量精度0.1 μg。根据DTG曲线和前人的研究，利蛇纹石脱水温度范围定义为400–800 °C，氢氧镁石为300–700 °C。对于Ol–Lz和Prd–Lz混合物，测量400–800 °C之间的重量损失，利用利蛇纹石理论H₂O含量（13.0 wt%）计算参考利蛇纹石含量。对于Prd–Brc混合物，评估了全脱水范围（300–700 °C）和主要脱水范围（300–475 °C）两个温度区间。三元Ol–Lz–Brc混合物中两种含水矿物共存，采用475 °C作为分界温度，将300–475 °C的重量损失归属于氢氧镁石，475–800 °C的归属于利蛇纹石。

拉曼光谱使用Horiba LabRAM HR800 Evolution共聚焦拉曼显微镜，空间分辨率优于1 μm（水平）和2 μm（垂直）。粉末样品压平后先进行单点测量，随后在随机选择的400 × 400 μm²区域进行映射扫描，间隔80 μm，每个映射扫描获得36个点（6 × 6网格）的光谱。使用532 nm激光，单点测量积分时间5 s，映射扫描2 s。利用3600–3720 cm^?1羟基伸缩振动带的积分面积建立校准曲线，利蛇纹石特征峰在3687 cm^?1，氢氧镁石在3655 cm^?1。三元混合物中使用三个积分范围：全带（3600–3720 cm^?1）用于利蛇纹石+氢氧镁石总含量，3667–3720 cm^?1仅用于利蛇纹石，3636–3667 cm^?1仅用于氢氧镁石。

FTIR使用Thermo Nicolet iS50 FTIR光谱仪，400–4000 cm^?1范围，分辨率4 cm^?1，32次扫描累积。纯KBr片测量作为背景。KBr压片由约1 mg样品粉末与200 mg KBr充分混合制成约500 μm厚度，每个样品三次重复测量。利用SiO伸缩带（约700–1140 cm^?1）和–OH伸缩带（3600–3720 cm^?1）的积分面积建立校准曲线，采用归一化比值A_–OH/(A_Si–O + A_–OH) × 100%进行定量，所有光谱吸光度均在比尔–朗伯定律线性范围内（<0.7）。

XRF分析用于确定样品主量元素组成，粉末样品105 °C干燥后1000 °C灼烧计算烧失量（LOI），约1 g样品与Li₂B₄O₇–LiBO₂–LiNO₃熔剂混合后在铂坩埚中1050 °C熔融，玻璃片用Pananalytical XRF光谱仪分析，分析不确定度小于0.1 wt%。电子探针分析使用JEOL JXA-iSP100电子探针显微分析仪测定橄榄岩中橄榄石和斜方辉石的主量元素组成，加速电压15 kV，束流20 nA，束斑直径3 μm。

热液蛇纹石化实验采用与标准橄榄岩相同的橄榄岩粉末（约10.0 g），装入200 mL钛合金热液釜中，加入去离子水（约100.0 g），抽真空除去空气后充入N₂至初始压力1 MPa防止氧化。三组实验分别在200 °C/3.3 MPa、250 °C/5.8 MPa和300 °C/10.3 MPa条件下进行，运行5–15天以获得不同蛇纹石化程度的样品。实验后自然冷却至室温，固体产物回收后在85 °C干燥48 h，玛瑙研钵中研磨确保均匀，全部实验产物用上述四种定量方法分析。

研究结果表明，三种纯矿物和天然橄榄岩的主量元素组成通过XRF确定，橄榄石化学式为(Mg_0.91Fe_0.09)₂SiO₄，利蛇纹石和氢氧镁石的LOI分别为13.47 wt%和31.20 wt%，化学式分别为(Mg_0.9Fe_0.1)₃Si₂O₅(OH)₄和Mg(OH)₂。橄榄岩XRD图谱仅显示橄榄石和斜方辉石的典型峰，归类为方辉橄榄岩，副矿物斜长石和尖晶石低于XRD检测限。基于SiO₂质量平衡计算，橄榄石和斜方辉石比例分别约为85 wt%和15 wt%。

校准曲线方面，XRD Rietveld方法的Ol–Lz校准曲线（y = (0.971 ± 0.025)x，R² = 0.9961）和Prd–Lz校准曲线（y = (1.002 ± 0.028)x，R² = 0.9953）与1:1参考线吻合良好，Prd–Brc校准曲线（y = (0.967 ± 0.014)x，R² = 0.9988）略有偏差。PHR和PAR校准曲线均拟合良好（R² > 0.99），拟合参数A范围为2.351–8.160。三元体系中Rietveld精修对Lz–Brc混合物中氢氧镁石含量重现良好（斜率= 1.024 ± 0.030，R² = 0.9966），但利蛇纹石低估约8%（斜率= 0.919 ± 0.021），三元混合物中氢氧镁石高估约24%（斜率= 1.236 ± 0.047）。TG校准曲线Ol–Lz和Prd–Lz均呈优异线性拟合（R² = 0.999），斜率分别为0.1209 ± 0.0015和0.1215 ± 0.0017，与理论斜率0.13一致。Prd–Brc的300–475 °C和300–700 °C两种校准方法均拟合优异（R² = 0.9996），斜率分别为0.2732 ± 0.0018和0.3056 ± 0.0018，后者与氢氧镁石理论校准曲线（y = 0.309x）接近。三元Ol–Lz–Brc混合物中利蛇纹石（475–800 °C）和氢氧镁石（300–475 °C）校准曲线均优异（R² > 0.997），斜率分别为0.1253 ± 0.0010和0.2930 ± 0.0056，均高于二元校准，表明两种矿物因脱水区间重叠而各自捕获了对方的部分重量损失，导致二者均被高估。Lz–Brc混合物中300–475 °C重量损失随利蛇纹石比例增加而线性降低（R² = 0.9919），从纯利蛇纹石的约0.9 wt%降至纯氢氧镁石的约27.5 wt%。

拉曼校准曲线R²为0.9433–0.9849，但测量分散性大，Ol–Lz和Prd–Lz曲线在低利蛇纹石比例（<40 wt%）处系统性高估，归一化峰面积接近实际混合比的两倍，表明拉曼方法在低蛇纹石化程度下不可靠。三元混合物中全–OH带积分面积与利蛇纹石+氢氧镁石总含量的R² = 0.8934，利蛇纹石子区间（3667–3720 cm^?1）R² = 0.9317，氢氧镁石子区间（3636–3667 cm^?1）R² = 0.5936，相关性差。

FTIR校准曲线基于Si–O和–OH带积分面积，全程拟合优异（R² > 0.99），测量不确定度显著较低。Ol–Lz和Prd–Lz曲线呈线性，斜率分别为0.0741 ± 0.0025和0.0714 ± 0.0014，Prd–Brc曲线向上凹，用二次函数拟合（y = (0.0100 ± 0.0003)x²，R² = 0.9933），因氢氧镁石无SiO带，归一化比值分母随氢氧镁石增加而减小。三元混合物中归一化比值与利蛇纹石+氢氧镁石总含量线性相关（y = (0.0631 ± 0.0022)x，R² = 0.9925）。

精密度和准确度评估表明，TG分析产生最准确的校准曲线，三个二元系列拟合误差均低于±1.5%（Ol–Lz为±1.2%，Prd–Lz为±1.4%，Prd–Brc为±0.8%），Lz–Brc混合物除外（±3.6%）。XRD方法中Rietveld表现最佳（±1.4%–2.9%），其次为PHR（±1.8%–4.9%）和PAR（±2.8%–6.5%）。FTIR光谱误差中等（±2.0%–3.4%），拉曼光谱误差最大（±5.1%–10.0%）。三元体系中TG分析仍最准确（利蛇纹石±0.8%，氢氧镁石±1.9%），FTIR对利蛇纹石+氢氧镁石总含量误差为±3.5%，XRD方法显著恶化（PHR氢氧镁石±14.8%，PAR±12.9%），Rietveld相对稳定（±2.3%和±3.8%），拉曼表现最差（利蛇纹石±11.0%，氢氧镁石±32.0%）。各情况下氢氧镁石误差均超过利蛇纹石，表明三元系列中氢氧镁石的低丰度显著降低所有四种方法的准确度。

讨论部分对四种方法进行了系统评估。XRD具有原位观测和无需物理标样校准的优势，Rietveld精修仅用数据库晶体学参考数据即可实现准确定量，但存在四方面局限：优选取向影响（尤其对片状层状硅酸盐矿物）、结晶度差异（实验产物中蛇纹石衍射峰宽而低，与校准标样差异显著，本研究标准利蛇纹石(001)峰FWHM为0.176°，而热液实验产物5天后FWHM为0.415°、15天后为0.267°）、低丰度相检测能力有限以及第三相存在时峰比方法急剧退化。TG分析的主要优势在于质量损失模式仅由化学组成和矿物种类决定，不受结晶度、粒度或优选取向影响，对仅含单一含水矿物的体系高度有效，但面临蛇纹石和氢氧镁石脱水温度范围不完全分离的关键限制，二者共存时使用单一分界温度引入系统误差，本研究三元实验直接量化了该效应：采用475 °C分界，利蛇纹石校准斜率从二元0.1209升至三元0.1253，氢氧镁石从0.2732升至0.2930，二者分别被高估3.6%和7.3%。含绿泥石和菱镁矿等热分解范围重叠的矿物时TG分析可靠性进一步受限。拉曼光谱因激光光斑微小（<2 μm），每次测量采样面积远小于粉末样品组成变化尺度，定量需通过多点分析或映射重建相比例，在低蛇纹石化程度下系统性高估，三元混合物中氢氧镁石拟合误差达±32.0%，但保留区分蛇纹石多型和快速无损鉴定的优势，最适角色为相和多型鉴定及相间空间关系表征而非总体定量。FTIR光谱为总体定量方法，KBr压片整合整个约1 mg粉末的吸收，校准曲线全程拟合优异，误差±2.0%–3.5%仅次于TG分析，但需注意归一化比值仅在基质和含水相SiO带强度可比时近似线性，且利蛇纹石（3682 cm^?1）和氢氧镁石（3692 cm^?1）的–OH伸缩带仅相隔约10 cm^?1无法分辨，三元混合物中仅能恢复利蛇纹石+氢氧镁石总含量。

将四种方法应用于热液蛇纹石化实验产物，XRD和SEM-EDS表征确认实验产物由蛇纹石（利蛇纹石）及不同量残余橄榄石和辉石组成，未检测到氢氧镁石，蛇纹石含量随反应时间（5–15天）从低程度（10–20 wt%）增至中等程度（约40 wt%）。实验蛇纹石结晶度显著低于校准标样，(001)反射远宽于标准利蛇纹石（FWHM = 0.21–0.55° vs 0.176°），–OH伸缩带也更宽（FWHM = 38.7–47.1 cm^?1 vs 29.2 cm^?1），宽33%–61%，该结晶度差异是方法间差异的主要来源。TG分析因对结晶度和优选取向效应不敏感而作为参考基线。三种XRD方法产生系统性分歧结果，Rietveld方法对低蛇纹石化程度样品尤其不可靠，甚至出现短时间实验蛇纹石含量高于同条件下长时间实验的异常结果，PHR和PAR方法系统性低估蛇纹石含量，且PHR一致低于PAR。FTIR定量结果关键取决于数据处理方法和积分策略，FTIR-T（本研究校准曲线）相对TG基线系统性高估蛇纹石含量，FTIR-H（Huang et al., 2023; Huang et al., 2017b校准曲线）系统性低估，差异源于归一化所用SiO带不同。拉曼光谱变异性最大，同一样品重复测量间蛇纹石含量估计分散显著。这些方法间差异对解释蛇纹石化动力学具有重要意义，方法选择关键影响蛇纹石化程度测定，尤其对低至中等蛇纹石化程度（<50 wt%）样品，差异主要源于结晶度不匹配，但方法一致应用时系统偏差恒定并在相对比较中抵消，不影响相对趋势研究。

方法选择建议方面，应用定量方法前应进行初步XRD分析和/或岩相观察以确定矿物组合并评估样品结晶度。对于仅含单一含水矿物的简单体系，推荐TG分析，测量重量损失结合校准曲线估计矿物比例，精度±0.8%–1.4%，对结晶度差的样品尤具优势。对于蛇纹石和氢氧镁石共存的体系，TG分析仍为推荐主要方法但需谨慎解释，因脱水范围重叠任何单一分界温度均导致两相各自捕获对方部分重量损失而均被高估，偏差大小强烈取决于利蛇纹石/氢氧镁石比例，蛇纹石主导组合中偏差适中而氢氧镁石富集组合中显著，建议用XRD Rietveld精修独立确认矿物比例，不推荐XRD峰高比和峰面积比方法（第三相存在时准确度急剧下降），FTIR可可靠测定蛇纹石+氢氧镁石总含量但无法单独分辨两相，拉曼可用于多型判别但不应依赖其定量次要含水相。对于含多种含水矿物和碳酸盐的复杂体系，本研究讨论的任何单一方法均无法实现可靠定量，必须结合SEM-EDS和/或EPMA等微分析技术进行相鉴定和模式丰度估算，并采用多种互补方法相互印证。

结论部分，本研究系统评估了四种常用蛇纹石化程度定量方法（XRD、TG分析、拉曼光谱和FTIR光谱），首先用三个二元标样混合物系列（Ol–Lz、Prd–Lz和Prd–Brc）和三元Ol–Lz–Brc系列建立全部四种方法的校准曲线以定量其精密度和准确度，随后应用于不同温压条件下实验室蛇纹石化实验产物。主要认识和推荐如下：TG分析对含单一含水矿物的简单体系具有最高精密度和准确度，但蛇纹石和氢氧镁石脱水温度范围重叠，任何分界温度均导致两相各自捕获对方部分重量损失而均被高估，偏差大小和方向取决于利蛇纹石/氢氧镁石比例及重量损失分配方式：采用独立校准曲线和475 °C分界时两相均被高估（本研究），而基于理论H₂O含量定量则蛇纹石高估、氢氧镁石低估（Miller et al., 2016）。XRD方法（Rietveld、PHR和PAR）受优选取向和结晶度变化强烈影响，Rietveld方法对低蛇纹石化程度样品尤其不可靠，峰高比和峰面积比方法在第三相存在时进一步恶化而Rietveld精修相对稳定。FTIR光谱精密度优于拉曼光谱，但定量准确度关键取决于数据处理协议和积分策略的选择，三元混合物中可可靠恢复蛇纹石+氢氧镁石总含量但无法单独分辨两相。拉曼光谱因在非均质样品中采样面积有限而变异性最大，无法定量次要量含水相。将这些方法应用于热液实验产物揭示出大的方法间差异，主要源于结晶度差的反应产物与结晶良好的校准标样之间的不匹配。因此强调实践中方法选择应针对特定矿物组合和蛇纹石化程度量身定制，对于含多种含水矿物和碳酸盐矿物的复杂体系，必须结合SEM-EDS和/或EPMA等微分析技术。这些发现为改进蛇纹石化研究的数据准确度和实验室间可比性提供了实用指南。