B5 · 出版卷 10
多源地球化学证据融合
Whole-rock、mineral chemistry、isotopes、context
学习目标
完成本课后,你应能够融合野外关系、岩相学、矿物化学、全岩组成、水溶液或稳定同位素资料和地质年代学,而不把它们当成彼此独立的投票;能够建立证据图、审核组成变换、分离物源—过程—时间假设,并为地质结论定义适用边界。
融合不是把很多图摆在一起,而是构建因果与时间模型:每项观测都有采样支持,每项变换都有声明方程,每项推断都有替代模型。矛盾必须保留,因为它常能指出遗漏过程或混合微区。
从对象到事件的证据图
把项目表示为图结构。对象节点 包括露头、区间、颗粒、生长带、脉体和水样;测量节点 包括观测、图像、浓度和同位素比;变换节点 包括空白校正、归一化、对数比、质量平衡、形态计算和日期计算;模型节点 包括混合、反应、分馏、扩散和衰变;事件节点 包括形成、蚀变、变形、冷却和重置。
边必须有动词:采自、测得、经何校正、由何计算、支持、反驳或穿切。直接野外关系不能与区域类比具有相同逻辑地位。每项事件主张都应向后追溯到对象,向前产生可检验预测。
不同数据列并不自动独立。矿物化学和全岩化学可能来自同一样品,年龄和微量元素可能来自同一点位,多个同位素比值可能共享校正。必须建模共享不确定性与来源,避免重复计数制造虚假信心。
组成变换与可比性
地球化学表可同时包含氧化物质量百分数、元素 mg/kg、摩尔比、同位素比和删失值。转换必须记录式量、价态和基准,且保留原始字段。任何归一化表都需保存分母和缺失值规则。
封闭组成携带相对信息。对正组分 x_i,相对于分母 x_D 的加性对数比为
$alr_i=\ln\left(\frac{x_i}{x_D}\right).$
中心对数比把每个组分与几何平均比较,等距对数比使用正交平衡。所选坐标应表达地质对比,并明确处理零值或检出限。把所有删失值替换成同一检出限分数会制造人工聚类。
分类与判别图具有组成、标定和地质适用域。它们不能把蚀变岩变回原始熔体,也不能把混合沉积物自动变成物源指纹。绘图前要检验蚀变、堆晶、粒级和分析可比性,并把图解作为岩相和质量平衡旁的条件证据。
空间与时间融合
位置、深度和微区决定数据能否组合。米尺度蚀变前缘、千米尺度物源区和微米尺度定年域回答不同问题。地图与剖面应携带采样支持和位置不确定性;插值不能仅因欧氏距离近而穿越地质边界。
相对时间约束数值时间。穿切、加生和交代先建立顺序,再解释年龄。若数值日期违背可靠关系,可能表示微区混合、重置、事件归属错误或野外模型错误;不能通过悄悄重命名早期资料消除冲突。
事件表应包含观测关系、定年微区、日期、不确定性、解释和替代模型。区域对比只有在物质连续性、地层、化学或年龄支持时才成立;相似组成可以重复出现,相似年龄也可以对应不同过程。
因果模型与竞争假设
至少写出两个过程模型。例如元素富集可能来自物源加入、流体迁移、残余浓集、矿物分选或分母损失。每个模型应预测矿物载体、互补亏损、空间梯度、同位素效应和时序,并围绕不同预测设计检验。
贝叶斯框架可以显式表达证据更新,但先验和似然资料不足时不必强行给数值概率。定性证据矩阵可以记录每项观测在各模型下是预期、可能还是矛盾。相关观测不能被当成独立似然重复相乘。
模型复杂度必须由证据赢得。在残差没有要求前,两端元混合优于五储库叙述;反过来,把明显多阶段微区强行拟合为一条直线也不是简约,而是证据损失。
结论必须包含适用边界:在哪些区域和条件下可使用模型,哪些样品被排除及原因,哪项不确定性阻止外推。通用教程的成果是方法与边界,而不是对具名地点或机构的主张。
可复现性、参考数据与审核
为对象、样品、制备、分析、参考物质和派生数据分配稳定标识符,使用受控单位与限定符,记录热力学数据库、衰变常数、参考组成、分类边界和代码版本。所有派生图都应能从不可变输入重建。
参考物质具有用途和基体,参考值并不是脱离方法的普适真值;适合某浓度范围的物质未必验证另一个范围。空白、重复样和参考结果应与接受标准一起报告,质量控制失败时必须停止解释并保留记录。
独立审核者应能遍历证据图、重算关键结果并挑战事件归属。网页上线不等同于科学审核。某项真实分析的作者或审核人只能出现在其真实项目来源记录中,不能据此给任何公司或个人赋予与本通用教程的关系。
合成综合案例
一个虚构样线穿过未蚀变岩浆岩、裂隙中心蚀变带和晚期脉体。合成数据包括十二个全岩样、四个样品的矿物图、成对稳定同位素分析,以及核部、边部和脉矿物日期。
全岩 A、B、C 三组分由带外 50:30:20 变为带内 40:40:20。原始百分比看似 A 损失、B 增加。以 C 为对数比分母,
$\ln(A/C):\ \ln(50/20)=0.916\rightarrow\ln(40/20)=0.693,$
$\ln(B/C):\ \ln(30/20)=0.405\rightarrow\ln(40/20)=0.693.$
矿物图显示含 A 长石被破坏、含 B 蚀变相生长,C 赋存在抗风化副矿物。两个候选参照的质量平衡支持中等 A 损失和 B 加入,但副矿物机械富集处范围更宽。
稳定同位素只在交代微区改变,未蚀变核部保持不变。年龄定义一个老核群、相干岩浆边群和较年轻脉体群;两个临近裂隙的边部部分重置。一个露头中脉体晚于蚀变,另一处两者使用同一通道。
首选模型为:岩浆边记录结晶;后来流体沿裂隙引起交代;脉体在流体活动后期或之后沉淀;渗透裂隙附近发生局部重置。替代模型是两次流体事件,因为穿切关系和同位素变化并非处处一致。适用边界排除机械分选地表样与混合烧蚀微区。本案例不作任何经济、组织或具体地点主张。
完成流程与复核
- 盘点对象、支持、单位、制备和质量控制状态。
- 分离原始、校正、归一化和模型字段。
- 数值融合前建立野外与岩相相对时序。
- 核查电荷、质量和同位素承载元素平衡。
- 使用与地质对比相适应的组成变换。
- 把年龄连接到成像微区,并比较系统不确定性。
- 建立至少两个因果—事件模型和判别预测。
- 制图适用范围、排除项和未决边界。
- 发布证据账本、可复现计算和简明结论。
完成复核应回答:每项结果能否追溯到真实观测或明确合成算例?公司和个人名称是否只在书目来源标识中出现?参考尺度、常数、数据库版本与校正是否记录?分析日期是否与解释事件分开?矛盾和剔除分析是否可见?审核者是否知道哪项证据能改变结论?
最终成果是一份年龄—地球化学数据审核,清楚区分测量、日期、过程、事件和不确定性来源。只有在首选地质叙述改变后证据图仍可审核,成果才算通过。
来源与延伸阅读
- Aitchison,组成数据统计分析,封闭地球化学组成对数比分析的基础。
- Winchester 与 Floyd,利用不活动元素进行地球化学判别90057-2),有影响力但必须保留适用域假设的分类方法。
- GeoReM 地质与环境参考物质数据库,参考物质值的可溯源汇编。
- Jochum 与 Nohl,地质参考物质的参考值,经评估的浓度资料与不确定性语境。
- 科学数据 FAIR 指导原则,可复用、可互操作和可审核证据框架。