次生矿物有序-无序状态对地热水与地层水多组分地热温标的影响

《Geochimica et Cosmochimica Acta》：Impact of order–disorder in secondary minerals on the multicomponent geothermometry of formation waters and geothermal waters

　　

多组分地热温标是估算地下温度的关键工具，但其精度常受矿物晶体有序-无序状态的显著影响，而标准热力学数据库往往未能充分表征这一状态。本研究通过系统性地纳入关键次生矿物的可变有序形式，解决了这一局限。研究人员首先利用Landau理论计算了常见次生矿物（包括长石、碳酸盐、云母、铁氧化物、硫化物、榍石及石英）在有序-无序参数Q值从0（完全无序）到1（完全有序）范围内的无序热力学性质（吉布斯自由能、焓、熵、体积）。这些数据被用于显著扩充SUPCRT92热力学数据库及后续的PHREEQC-LLNL数据库。增强后的地热温标方法结合修改后的数据库，应用于美国加利福尼亚Kettleman North Dome油田Temblor层的两组地层水，以及冰岛Reykjanes与Svartsengi、新西兰Tauhara、哥斯达黎加Miravalles、土耳其Kizildere等地热田的十五组地热水。案例研究表明，选择具有合适Q值的矿物形式，可使矿物饱和指数（log Q/K）在合理的平衡温度下实现更好的收敛，相较于仅使用完全有序或未明确有序状态的矿物计算，其结果往往与独立估算的储层条件或已知的成岩/热液矿物组合更为吻合。这一成果至关重要，因为晶体有序-无序显著影响着许多固相的热力学稳定性与反应性。这些矿物作为油气田地层的成岩相和热液蚀变产物，在决定流体-岩石平衡中起着关键作用。忽视其特定的有序-无序状态可能导致错误的温度估算和对储层条件的误判。结果强调了在稳健地球化学模拟中考虑有序-无序现象的重要性，并凸显了扩充后的热力学数据集对未来地热温标应用的实用价值。

本研究发表于《Geochimica et Cosmochimica Acta》，针对多组分地热温标中因忽略矿物有序-无序状态导致温度估算偏差的问题展开。传统热力学数据库通常仅收录矿物的单一有序态，未考虑实际地质过程中矿物结晶速率、冷却历史对离子占位的影响，使得流体-岩石平衡模拟结果与实测储层温度存在系统性误差。为此，研究人员系统计算了长石、碳酸盐、云母等19种次生矿物在有序-无序参数Q（0至1区间）下的热力学性质，构建了包含可变有序态的SUPCRT92与PHREEQC-LLNL数据库，并将其应用于全球5个地热系统与1个含油气系统的17组水样模拟，证实引入Q参数可显著提升温度反演精度。

研究采用三个关键技术方法：一是基于Landau理论计算矿物无序热力学性质，覆盖Q值0至1、步长0.1的区间；二是通过参考反应法将新矿物与原有数据库关联，确保热力学内部一致性，并将数据录入SUPCRT92后回归得到PHREEQC所需的温度-溶度积系数；三是采用多组分地热温标模拟，针对水样缺失组分（如Al、Fe³⁺）通过假设与特定有序态矿物平衡进行约束，并强制氧化还原非平衡以避免模拟前后化学组成漂移。样本涵盖美国加州Kettleman North Dome油田Temblor层两组地层水，以及冰岛、新西兰、哥斯达黎加、土耳其五地十五组地热水。

研究结果如下：

4.1 美国加州Kettleman North Dome油田287-19J井（Temblor层II带）：该区地层水平衡温度为105°C，模拟显示其与Q=0.76的冰长石、Q=0.72的白云母、Q=0.67的方解石、Q=0.25的白云石、Q=0.77的钛铁矿及Q=0.48的磁铁矿达到平衡，石英则以完全无序态接近饱和。各矿物饱和指数曲线在105°C附近收敛，与实测温度一致。

4.2 美国加州Kettleman North Dome油田272R-32J井（Temblor层V带）：平衡温度为100°C，对应Q=0.78的冰长石、Q=0.57的白云母、Q=0.65的方解石及Q=0.53的磁铁矿，石英、钠云母与钛铁矿以Q=0.9的高有序态接近平衡。

5.1 冰岛Reykjanes 8井：平均井温284±6°C，水样与Q=0.80的钠长石、Q=0.61的冰长石、Q=0.20的白云石及Q=0.58的钛铁矿平衡；Q=0.9的铁白云石与有序/无序绿帘石的饱和温度与井温偏差在7°C以内。

5.2 冰岛Svartsengi 4井：平均井温242.2±0.4°C，与Q=0.55的冰长石、Q=0.90的钠长石及Q=0.61的钛铁矿平衡；绿帘石的有序与无序态饱和温度均略高于井温。

5.3 新西兰Tauhara TH-06井：多矿物饱和温度平均为256.2±2.2°C，与井口测温及放汽焓推算的253°C吻合。

5.4 哥斯达黎加Miravalles PGM系列井：采用“FixAl”法（假设与有序钠长石平衡）约束Al浓度，1984至1994年间四口井的平均平衡温度为235-242°C，与历史生产数据一致；低SiO₂样品的玉髓饱和温度偏低，石英-玉髓混合态则偏高。

5.5 土耳其Kizildere KD系列井：同样采用“FixAl”法，四口井平均平衡温度为202-214°C，与储层测温范围匹配。

讨论部分指出，引入Q参数后，模拟识别的成岩/热液矿物组合与实际钻孔岩心观察高度吻合，仅少数矿物存在差异。研究证实，矿物的有序-无序状态直接影响其热力学稳定性，忽略该因素会导致温度估算偏差。结论强调，扩充后的热力学数据库为多组分地热温标提供了更可靠的工具，未来可推广至更广泛的地质流体模拟，同时呼吁进一步开展矿物有序-无序热力学性质的实验测定工作。