晶体学与矿物检测

每种矿物都是晶体。它的原子按照具有特定对称性的晶格重复排列,这种对称性体现在矿物如何偏折 X 射线和偏振光上。我们构建机器学习模型,根据这类测量识别矿物,并确保能看清它们是如何得出答案的。

  • Na
  • Cl
  • a = 5.640 Å
拖动可以旋转
图 1 岩盐型结构,石盐 (NaCl)、钾盐 (KCl) 和方铅矿 (PbS) 都具有这种结构。两种离子沿立方体的每条棱交替排列,立方体在各个方向上不断重复。每种矿物的离子都按比例绘制。

用 X 射线读出晶格

X 射线穿过晶体时,每一层原子平面都会反射一小部分光束。来自相邻平面的反射所走的距离略有不同。在大多数角度上,它们的相位错开而相互抵消。只有在少数几个角度上,多走的距离恰好是波长的整数倍,反射相互叠加,探测器记录到一个峰。这就是布拉格定律,2d sin θ = nλ。

衍射仪扫过各个角度,记录下这些峰。峰的位置给出晶格的间距,峰的高度则取决于哪些原子位于何处。

d θ 2d sin θ = nλ

强度

2θ (°)

  • 上一种矿物
图 2 用铜靶 X 射线计算得到的粉末衍射图谱。左图中,第二条射线多走了 2d sin θ 的距离,用彩色标出。石盐和钾盐的结构相同,所以它们的峰按同样的顺序出现。在钾盐中,钾和氯对 X 射线的散射几乎一样强,因此有几个峰几乎消失,第一个峰就是其中之一。峰的位置符合布拉格定律,峰的高度是近似值。

显微镜下的矿物

X 射线衍射读的是粉末,显微镜读的是岩石本来的样子,一个颗粒一个颗粒地看。薄片是把岩石磨到大约 30 微米厚的切片,薄到光可以透过。在偏光显微镜下,每种矿物都有自己的样子:颜色、突起,以及在正交偏光下转动载物台时亮度如何变化。矿物学家凭这些线索用肉眼识别颗粒,图像模型要学的也正是这些线索。

  • 石英
  • 斜长石
  • 黑云母
  • 石榴子石
  • 不透明矿物(硫化物)

点击一个颗粒,查看识别它的依据。

图 3 一块模拟薄片。可以在单偏光和正交偏光之间切换,并转动载物台。在正交偏光下,大多数颗粒每转一圈消光四次,石榴子石在任何角度都保持黑暗,不透明颗粒在两种光下都是黑色。颗粒标注显示的是图像模型学着输出的结果。

用机器学习识别薄片图像是一个活跃的领域。最近一个公开数据集包含 25 种矿物、超过 10 万个由专家标注的颗粒,同时看到单偏光和正交偏光图像的模型,比只看到一种图像的模型表现更好。我们想知道的是,这些模型依靠的是矿物学家使用的线索,还是一些捷径,而这些捷径换到下一座矿山的岩石上就会失效。

模型做什么

来自矿山的样品很少只含一种矿物。它的衍射图谱混合了几种矿物,峰相互重叠,上面还叠加着噪声,它的薄片里则有几百个颗粒。靠人工识别很慢。我们训练模型预测晶体的晶格和空间群,空间群就是晶体全部对称性的完整列表,同时训练模型估计样品中含有哪些矿物。显微镜图像是自然而然的下一步。

这些模型在设计上遵守对称性。晶体在仪器中转了一个方向,仍然是同一个晶体,所以模型的答案不应该因为旋转而改变。

可解释性为什么在这里重要

空间群正是我们在可解释性研究中所说的那种群。我们用来找出群乘法网络内部算法的方法,也可以用来检验晶体学模型学到了什么。地质学家应该能够确认,模型之所以报告黄铁矿,是因为看到了黄铁矿的峰,或是颗粒在显微镜下的样子,而不是因为噪声中的某种图案。

矿物组成也是地质冶金模型用来预测矿石在选矿厂中表现的性质之一。这正是这项工作与我们的深部矿山系统研究的交汇之处。