MagMin_PT——面向岩浆岩的Excel矿物分类与地质温压计程序
- 2026-11-06 22:48:13
一、论文核心信息提炼
核心结论：MagMin_PT是一个免费、开源的Excel工作簿，可将电子探针（EPMA）成分数据快速转化为常见岩浆矿物的结构式、端元组分、分类图解和温压（P–T）结果，并内置矿物–熔体平衡检验，降低了岩浆岩矿物学与温压计算的软件门槛。
标题：MagMin_PT: An Excel-based mineral classification and geothermobarometry program for magmatic rocks
作者：Mesut Gündüz, Kürşad Asan
发表期刊：Mineralogical Magazine
发表年份：2022
二、研究背景与核心科学问题
背景
岩浆岩中的造岩矿物记录了从地幔源区、岩浆房到浅部地壳的物理化学条件（温度T、压力P、氧逸度fO₂、挥发分含量等），这些信息以主量元素、微量元素和同位素组成的形式“密码式”地保存在矿物和全岩成分中（图1）。电子探针（EPMA）是获取矿物成分最常用的手段，但其输出为氧化物重量百分比，必须经过阳离子结构式换算和端元组分计算，才能用于矿物分类和地质温压计（geothermobarometry）。
此前已有大量Excel或独立程序，但多数是矿物专属的计算器：角闪石（Yavuz, 1996, 2007; Locock, 2014; Li et al., 2020a）、云母（Yavuz and Öztas, 1997; Li et al., 2020b）、辉石（Yavuz, 2013）、石榴子石（Yavuz and Yildirim, 2020）、磁铁矿–钛铁矿（Lepage, 2003; Yavuz, 2021）等。少数通用程序（如Afifi and Essene, 1988; Brandelik, 2009）也主要停留在结构式或分类层面，缺乏跨矿物的综合温压计算能力。因此，研究者需要一个用户友好、无需编程、覆盖常见造岩矿物、并集成平衡检验与多种温压计的Excel平台。
图1
核心科学问题
如何将EPMA矿物成分数据标准化、自动化地转换为结构式、端元组分和矿物分类结果？
如何在一个统一界面中集成橄榄石、辉石、角闪石、黑云母、长石、磁铁矿–钛铁矿、磷灰石和锆石的温压计、氧逸度计、岩浆水含量和饱和温度计算，并在计算前进行矿物–熔体平衡检验？
三、软件设计与研究方法
软件架构
MagMin_PT是一个Microsoft Excel工作簿，包含11个工作表：Instructions、Data Input、Olivine、Pyroxene、Amphibole、Biotite、Feldspar、Magnetite-Ilmenite、Apatite、Zircon和Conversions。操作顺序为（图2）：① Data input → ② 选择矿物 → ③ Calculation。由于部分温压计为迭代方程，用户需在Excel中启用“Enable iterative calculation”以避免循环引用错误。
图2
数据输入与结构式计算
输入数据为矿物氧化物的wt.%（EPMA分析），必要时还需输入全岩或玻璃成分，用于矿物–熔体温压计。结构式换算基于氧原子数和电价平衡：
| 矿物 | 计算基准 | 端元/分类 |
|---|---|---|
| 橄榄石 | 4个氧 | Fo–Fa |
| 辉石 | 6个氧 | 辉石 quadrilateral（图3） |
| 角闪石 | IMA 2012方案 | 8个二元分类图（图4） |
| 黑云母 | 11、22或24个氧 | 黑云母分类 |
| 长石 | 8、16或32个氧 | Ab–An–Or三元图（图5） |
| 磁铁矿/钛铁矿 | 4/3个氧 | FeO–Fe₂O₃–TiO₂三元图 |
| 磷灰石 | 25或26个氧 | 磷灰石族 |
| 锆石 | 4个氧 | ZrSiO₄ |
Fe³⁺估算是程序的关键环节：辉石和磁铁矿–钛铁矿采用Droop（1987） 的化学计量法；角闪石按Hawthorne et al.（2012）IMA 2012附录III的经验规则；黑云母除传统Dymek（1983）方法外，还提供Li et al.（2020b）的机器学习方案；角闪石也提供Li et al.（2020a） 的机器学习结构式计算作为IMA方案的替代。
图3
图4
图5
平衡检验
温压计的前提是矿物–熔体或矿物–矿物化学平衡。MagMin_PT提供三类检验：
岩相学观察：自形晶形通常指示平衡；熔蚀、反应边、环带、假象等为不平衡标志。
矿物–熔体分配系数：对橄榄石、斜方辉石、单斜辉石和角闪石计算Fe–Mg交换系数 K_D（Eq. 1–2）。可接受范围：橄榄石 0.30 ± 0.03；斜方辉石 0.29 ± 0.06；单斜辉石 0.27 ± 0.03；角闪石 0.28 ± 0.11。辉石还提供Rhodes图（液相Mg# vs. 矿物Mg#）进行图示检验。
预测–观测组分对比：主要适用于单斜辉石，将预测组分与实测组分投在1:1线上，偏离越小越接近平衡。
内置温压计与计算模块
| 矿物/模块 | 主要计算内容 | 代表性校准 |
|---|---|---|
| 橄榄石 | 橄榄石–熔体温度计 | Helz and Thornber (1987); Beattie (1993); Sisson and Grove (1993); Putirka et al. (2007); Putirka (2008) |
| 辉石 | 两辉石温度计/压力计；辉石–熔体温度计 | Wood and Banno (1973); Wells (1977); Lindsley and Andersen (1983); Brey and Kohler (1990); Putirka (2008) |
| 角闪石 | Al-in-hornblende压力计；角闪石–斜长石、角闪石–熔体温压计；氧缓冲剂；熔体水含量 | Hammerstrom and Zen (1986); Hollister et al. (1987); Johnson and Rutherford (1989); Blundy and Holland (1990); Schmidt (1992); Anderson and Smith (1995); Holland and Blundy (1994); Molina et al. (2015, 2021); Putirka (2016); Ridolfi et al. (2010); Ridolfi (2021) |
| 黑云母 | 黑云母温压计；氧缓冲剂 | Luhr et al. (1984); Uchida et al. (2007); Wones and Eugster (1965); Wones (1989) |
| 长石 | 斜长石–熔体、碱性长石–熔体、二长石温度计；水含量计 | Putirka (2005, 2008); Pichavant and Macdonald (2007) |
| 磁铁矿–钛铁矿 | 温度与fO₂；缓冲曲线（MW、QFM、NNO、HM） | Powell and Powell (1977); Spencer and Lindsley (1981); Andersen and Lindsley (1985); Sauerzapf et al. (2008); Frost (1991) |
| 磷灰石 | 磷灰石饱和温度 | Piccoli et al. (1999); Piccoli and Candela (2002) |
| 锆石 | 锆石饱和温度；Ti-in-zircon温度计 | Hanchar and Watson (2003); Boehnke et al. (2013); Watson et al. (2006) |
此外，程序可计算氧逸度（log fO₂）、岩浆水含量（wt.% H₂O） 和饱和温度，并将结果绘制为分类图、P–T图和fO₂–T图。所有图形和表格可导出为gif/jpeg/tiff或复制到其他软件中。
四、主要结果与软件功能
1. 实现了常见造岩矿物的“一站式”结构式与分类计算
MagMin_PT覆盖了岩浆岩中最常见的9类矿物（橄榄石、单斜辉石、斜方辉石、两辉石、角闪石、黑云母、长石、磁铁矿–钛铁矿、磷灰石、锆石）。用户输入EPMA数据后，程序自动完成阳离子分配、端元计算和分类投图（图3–图5）。例如：
辉石按Morimoto et al.（1988） 方案分配T、M1、M2位，并采用Droop（1987） 估算Fe³⁺，投在Lindsley and Andersen（1983） 的辉石等温线图上（图3）；
角闪石按IMA 2012（Hawthorne et al., 2012）计算，生成8个二元分类图（图4）；
长石按Ab–An–Or三元体系分类（图5）。
2. 将平衡检验前置，减少温压计误用
程序在计算P–T之前自动计算K_D(Fe–Mg) 并给出可接受范围。若K_D偏离合理区间，说明所选矿物–熔体对可能不平衡，温压结果不可靠。这一设计把岩石学判断嵌入软件流程，而不是让用户盲目套用温压计。
3. 集成了多套温压计，便于交叉验证
同一矿物往往有多种校准。MagMin_PT允许用户同时获得多个温压计的结果，比较其一致性。例如角闪石模块同时包含多套Al-in-hornblende压力计，但作者明确提醒：这些压力计适用条件受限，不能无条件用于所有岩浆系统（Erdmann et al., 2014; Putirka, 2016）。这种“多校准并列”的设计有助于识别系统偏差和适用性边界。
4. 支持机器学习辅助的Fe³⁺和结构式计算
对角闪石和黑云母这类含水、复杂固溶体矿物，传统化学计量法难以唯一确定Fe³⁺/Fe²⁺。MagMin_PT引入了Li et al.（2020a, 2020b） 的机器学习方法，与IMA 2012或Dymek（1983）方案并列，用户可比较不同方法的结果。这反映了近年矿物温压计领域的一个趋势：数据驱动方法正在补充传统经验校准。
5. 可计算氧逸度、岩浆水含量和饱和温度
除P–T外，程序还可输出：
氧缓冲曲线（MW、QFM、NNO、HM）及log fO₂（Eq. 5；Frost, 1991）；
角闪石成分计算的熔体水含量（Ridolfi et al., 2010; Ridolfi, 2021）；
斜长石–熔体hygrometer 和Al₂O₃ hygrometer（Pichavant and Macdonald, 2007）的水含量；
磷灰石饱和温度和锆石饱和温度、Ti-in-zircon温度。
这些参数将矿物成分与岩浆挥发分和氧化还原状态联系起来，扩展了程序的应用范围。
五、科学意义与结论
如何回答核心科学问题？
标准化与自动化：MagMin_PT把EPMA数据到结构式、端元、分类和P–T结果的流程封装在Excel中，减少了手工计算出错的风险，也使不同研究之间的结果更具可比性。
综合性与可及性：与矿物专属程序相比，MagMin_PT在一个工作簿中覆盖多类矿物和多种温压计；与需要编程或商业软件的平台相比，Excel界面无需编程基础，便于教学和常规研究使用。
对岩浆岩岩石学研究的启示
从“单矿物计算”走向“系统交叉验证”：同一岩石中不同矿物的温压结果可以相互检验。例如，橄榄石–熔体温度、辉石–熔体温度和锆石Ti温度若一致，可增强对岩浆房条件的信心；若不一致，则提示不平衡、多期岩浆或校准适用性问题。
平衡检验的“门槛”作用：K_D检验和Rhodes图把岩相学观察与定量计算连接起来，避免把不平衡矿物对投入温压计。这对火山岩中常见的环带、熔蚀和反应边晶体尤为重要。
对含水矿物的谨慎处理：角闪石和黑云母的Fe³⁺估算、结构式分配和压力计适用性仍是薄弱环节。MagMin_PT并列多种方案的做法，实际上提醒用户：含水矿物的温压结果应视为半定量约束，需结合岩相学和实验岩石学判断。
局限性与使用建议
程序依赖化学平衡假设；K_D落在“可接受范围”并不证明平衡，只是必要非充分条件。
Al-in-hornblende压力计仅在特定成分、温度和氧逸度范围内有效，不能外推到所有岩浆系统。
Excel迭代计算需要手动启用；大数据量或复杂公式可能导致计算缓慢。
部分温压计的经验校准范围有限，用户应检查其适用P–T–X区间，避免外推。
程序是计算工具，不能替代详细的岩相学观察和地质背景分析。
六、总结
MagMin_PT是一个面向岩浆岩的免费Excel矿物分类与温压计算程序，覆盖橄榄石、辉石、角闪石、黑云母、长石、磁铁矿–钛铁矿、磷灰石和锆石等常见造岩矿物。它将EPMA成分数据自动转换为结构式、端元组分和分类图解，并集成矿物–熔体平衡检验、多种温压计、氧逸度计、岩浆水含量和饱和温度计算。其核心价值在于降低计算门槛、促进标准化、支持多方法交叉验证；同时，作者通过并列多种校准和明确适用条件，提醒使用者温压计结果必须建立在平衡检验和岩相学判断之上。对于岩浆岩岩石学、火山学和中深部地壳过程研究，MagMin_PT提供了一个实用、透明且易于复现的计算平台。
