B3 · 出版卷 8
三维构造解释与不确定性
Fault framework、branching、termination、alternative models
学习目标
完成本课后,学习者应能分别定义观测、地质规则和插值;建立拓扑一致的断层面和褶皱面;整合地图、剖面、产状和深部控制;比较显式与隐式表达;传播测量、概念和插值不确定性;评价替代事件模型;设计一项最能区分模型的新观测。
从野外问题开始
一条地图断层迹线绕过褶皱,并被两个钻孔相交。钻孔一在预测深度发现窄断层,钻孔二却在首个模型预测位置上方 140 m 处遇到重复地层。断层是弯曲、重复、被较新断层错开,还是岩芯对比或地表外推错误?
稀疏数据之间几乎总能生成平滑曲面。地质解释增加连续、终止、位移、地层和事件时间规则,但这些规则本身是可检验假设。模型必须显示哪些几何是观测、哪些由插值计算,以及何处存在满足同一证据的替代概念模型。
核心过程模型
三维构造模型包含面、体、线和拓扑。断层面切割或限定岩体;层面终止于断层或侵蚀面;褶皱改变曲面;交汇产生线;位移改变断层两侧的对应关系。只有几何而没有拓扑,会产生不可能的穿插、空隙或重复岩层。
显式建模用三角网、样条或参数面表示曲面;隐式建模用标量场的等值面定义接触或面理。两者都需要输入和约束。增加网格密度或生成更平滑的场不会增加观测。
数据可包括地图迹线、产状、钻孔接触、定向岩芯、地震或其他地球物理解译、地形和相对时序。每类数据都有位置、方向、分辨率和对比不确定性。应统一坐标,同时保留原始来源。
插值依据选择的协方差、平滑度、曲率或构造规则估计观测之间的几何。超出数据包络的外推约束更弱。断层位移和褶皱形态有时需要恢复或运动学算法,不能把所有面独立平滑。
不确定性至少包含:
- 测量不确定性:位置、产状、深度、拾取和分析误差;
- 采样不确定性:不可见体积和覆盖不均;
- 插值不确定性:参数和数学表示;
- 概念不确定性:替代断层、地层对比、事件序列和变形机制;
- 实现不确定性:坐标转换、计算设置、网格分辨率和版本。
同一概念内的模型集合可以探索参数变化,但不同概念必须分别建立。只有在说明采样和先验假设后,概率体才具有可解释含义。
证据与测量
先建立数据登记表,为每项数据赋予稳定标识、坐标参考系、单位、采集方法、质量、不确定性和与机构无关的来源记录。解译接触与原始日志或图像分开保存。使用已知控制点和往返转换检验坐标。
剖面应作为同一模型的查询结果,而不是独立装饰图。地图、剖面和三维视图中的同一曲面必须是同一对象或经过协调。检查断层交汇、层面截断、地层顺序、厚度和位移一致性,并标记无约束终止。
用保留数据检验预测:暂时移除部分钻孔接触或产状,建立模型后再恢复并计算误差。空间聚集的验证点不能检验远处外推。拟合与失败都应报告。
模型发布还需要最小治理记录:输入数据版本、坐标参考系、单位、过滤规则、解释者决策、算法与参数、生成日期、已知冲突和适用尺度。任何重新插值都应生成新版本而不是覆盖旧模型。这样可以区分“数据改变”“规则改变”和“软件实现改变”,也能在结果差异出现时定位原因,而不把责任归给不可追踪的黑箱。
计算示例
一个合成地区记录两期变形:
第一期:层状岩石在低角度逆断层上方缩短,形成向东北倾伏的褶皱。S1 近轴面分布,早期 V1 脉切过层理后与层理一起被褶皱。
第二期:陡倾、走向西北的左行断层切过褶皱,错断逆断层,并在局部把 S1 再活化为脆性组构。晚期 V2 脉充填张性分支并切过 V1。
模型 A 把逆断层延续为一个缓弯曲面穿过两个钻孔,并给后期走滑断层 220 m 左行位移。它符合地表迹线、多数产状和钻孔一,却比钻孔二的重复地层深 140 m。
模型 B 增加一条第一期冲断分支,该分支被第二期断层切割和错开。它符合两个钻孔并预测东侧剖面还有一个截断,但分支没有地表直接露出,模型复杂度也更高。两种模型都符合高置信地图迹线和 V1/V2 时序。A 更简单却与钻孔二冲突;B 通过一个推断面解决冲突。
不能根据平滑度选择。首先复核岩芯对比和深度,然后在预测东部截断附近获取定向深部约束。若相交分支则支持 B;若地层连续则支持弯曲或位置修正后的 A。数据采集前必须记录预测。
对每个曲面赋予“已观测、紧密插值、宽松插值、外推”置信等级。在测量误差范围内生成几何变体,计算截断位置和岩体体积范围。模型 A 与 B 应保持为两个独立概念集合。
误判与不确定性
不能把精美渲染等同准确性。颜色、光照和密集网格会隐藏弱约束;首选模型必须与数据位置、不确定性和替代曲面同时显示。
不能用虚构的确定接触关闭数据空白。若岩性、年龄或地球物理不能对比单元,就应保留未解决状态。拓扑有效但地质错误的模型仍然错误。
避免重复计算派生证据。若地图迹线和剖面拾取来自同一张地球物理图像,它们不是独立证据,权重和验证必须保留共同来源。
模型导出的体积、深度和距离继承概念不确定性。只报告数值插值方差,会在断层数量或地层对比可能改变时严重低估不确定性。
实践调查
使用一条地图迹线、12 个产状、两条剖面和三个合成深部接触建立小型三维模型,保留一个接触和三个产状用于验证。分别建立最小复杂度模型和增加一条断层或褶皱域的替代模型,检查拓扑、地层顺序、断层截断和保留数据误差。
在记录的不确定性范围内扰动位置和产状,为每个概念重建模型集合。在目标深度绘制包络并计算交汇位置范围。编写模型卡,包含数据版本、坐标系、规则、参数、验证、冲突、替代模型和下一项检验。
掌握检查
- 几何与拓扑有何区别?
- 为什么平滑曲面不代表证据充分?
- 同一概念内参数扰动无法表示哪类不确定性?
- 保留数据如何检验解释?
- 什么使一项新观测具有区分力,而不是只增加数量?
资料与延伸阅读
- 三维地质建模中的不确定性,Wellmann 与 Caumon,2018。
- 开源随机地质建模与反演,de la Varga、Schaaf 与 Wellmann,2019。
- 地质图符号标准,观测与推断构造的一致表达。