近年来, 地球空间信息学和对地观测学迅猛发展^[1], 随着计算机图形技术、三维建模技术和可视化技术的迅猛发展, 地下三维可视化建模技术得到了日益广泛的应用。长期以来在传统的地质现象表达中, 通常是以二维平面图和剖面图来表达野外地质的成果^[2], 这种方式存在着空间信息的损失与失真、制图过程繁杂及信息更新困难。地质体三维可视化建模技术正是针对传统的地质信息模拟与表达方法的缺陷, 借助计算机和科学计算可视化技术, 直接从三维空间的角度去理解和表达地质体与地质环境。它具有形象、直观、准确、动态、丰富等特点^[3], 因此, 三维可视化研究有了越来越重要的现实意义和实用价值。在涌现出的众多地质三维可视化软件中, 处在世界领先地位的Micromine公司开发的Micromine软件在这些方面已取得了可喜的成果。本文利用国际矿业界通用的Micromine软件, 以东昆仑造山带1:5万水泥厂幅为例, 探索应用三维地质实体建模的技术和方法。

1.   研究区地质概况

格尔木南1: 5万水泥厂幅位于东昆仑造山带东段, 大地构造上位于昆中断裂带与昆南断裂带之间的东昆仑南地体^[4]。该地区山势雄伟, 相对高差达1000m。该区出露的岩石地层主要有:中-新元古界万保沟群大理岩和砾状灰岩; 奥陶-志留系纳赤台群水泥厂组灰岩、石灰厂组火山岩和碳酸岩、哈拉巴依沟组千枚岩和砂岩等; 下泥盆统牦牛山组砾岩夹砂岩及流纹岩; 上石炭统浩特洛洼组灰岩、板岩和砂岩; 中下三叠统洪水川组砂岩、板岩和砾岩; 新生代沉积沿雪水河和格尔木河谷分布。由于加里东期和印支期强烈造山运动, 上述地层被近东西向的逆冲断层分割成大小不等的构造岩片, 原始沉积层序已被打乱, 形成了一套总体有序、局部无序的构造地层体(图 1)。

图  1  东昆仑造山带水泥厂幅地质略图

1.新生界; 2.中下三叠统洪水川组; 3.上石炭统浩特洛洼组; 4.下泥盆统牦牛山组; 5.奥陶-志留系纳赤台群哈拉巴依沟组; 6.奥陶-志留系纳赤台群石灰厂组碳酸盐岩段; 7.奥陶-志留系纳赤台群石灰厂组火山岩段; 8.奥陶-志留系纳赤台群水泥厂组; 9.中新元古界万保沟群碎屑岩组; 10.中新元古界万保沟群碳酸盐岩组; 11.新元古代花岗斑岩; 12.逆断层

Figure  1.  Geological map of a Cement Plant in the East Kunlun orogenic belt

下载: 全尺寸图片 幻灯片

2.   Micromine软件简介

Micromine软件是澳大利亚Micromine国际矿业软件有限公司开发的从事地质勘探、资源评估、储量计算及露天和地下矿山开采设计的地理信息系统(GIS)三维立体空间平台软件。

软件已经汉化, 适用于WINDOWS95和WINDOWSNT及以上的平台, 自引进中国以来, 在国内多个省(区)的地勘单位、矿山企业取得了较好的使用效果, 受到矿业系统的好评; 软件的记录格式符合中国的格式要求, 具有和其他软件的接口, 使该系统的数据可被其他数据库管理系统及相关软件查询和编辑, 如可以导入CAD、Mapgis、Excel等格式的数据; 软件能够实现各种工程和矿体的三维立体显示和成图, 并根据地质统计学的原理提供进行矿体品位估值和储量计算的各种方法, 同时还提供了矿山开采设计与境界优化等功能; 软件可以制作演示工作成果的三维电影, 以利于进行工作汇报、成果展示和招投标等工作。

3.   三维建模的基本概念及原理

在三维建模的过程中, 地表模型的构建依托于数字地形模型(DTM)的生成, 地下模型的构建实际上是各地质体三维线框的生成以及相互之间布尔运算的结果。

3.1   DTM

DTM (Digital Terrain Model)-数字地形模型是利用一个任意坐标系中大量选择的已知x、y、z的坐标点对连续地面的一个简单的统计表示, 或者说, DTM就是地形表面形态属性信息的数字表达, 是带有空间位置特征和地形属性特征的数字描述。地形表面形态的属性信息一般包括高程、坡度、坡向等。

DTM系统主要是由计算机程序实现的, 应用于各种类型的计算机系统的DTM已经在许多国家开发成功, 尽管使用的方法不同、用户界面各异, 但主要功能都是从离散数据构建出相互连接的网络结构, 以此作为地形的数字模型基础。目前, DTM已经成为GIS的重要组成部分^[5]。

数字地形模型通过野外测量、航测、地图数字化等途径可得到设计所需的地形资料, 这些资料经过计算机处理后可以形成数字地形模型^[6~7]。通常数字地形模型建立的方法, 一是根据现在测量仪器(全站仪)得到的第一手地形测量数据(通常是点坐标)按照一定格式输入到相关软件后可以直接通过软件功能生成数字地形模型; 另一种是根据该地区提供的地质地形图进行矢量化后导入到软件中进行高程赋值, 然后直接生成数字地形模型。

3.2   三维线框实体基本原理

三维线框模型的构建主要是采用了TIN技术(不规则三角网模型)中的Voronoi图与Delaunay三角形算法。TIN是一种表示数字高程模型(DEM)的方法, 它既减少规则格网方法带来的数据冗余, 同时在计算效率(如坡度)方面又优于纯粹基于等高线的方法。这种表面模型可以进行体积估算、表面渲染、切制剖面、快速三维显示等操作, 比计算机图形学中的表面模型有所扩展, 能满足地质勘探中地质制图的基本要求, 也是建立三维实体模型的基础。地质体的地质形态复杂多变, 很难用规则的几何体来描述, 它需要一种灵活、简便、快速的方法来建立地质体的不规则几何模型, 这种模型正是基于这一需要而提出来的。

TIN不规则三角网技术  不规则三角网TIN属于场的数据模型范畴, 三角形单元遍布整个区域无重叠、无缝隙, 实体间的空间关系通过三角形连接表达。用于支持TIN建立的几何构造-Delaunay三角网因其特有的“外接圆规则”和“最大最小角规则”, 使得其成为空间邻近表达的有力工具。这些模型可以看作是不规则三角网TIN在平面空间场表达上发挥的具体作用。

Voronoi图  基于场论观念, 空白区域可看作是目标的相互“竞争”生成空间, 在群结构中, 每个面目标的“影响区域”通过竞争以中剖分结果达到平衡, 这是Voronoi图几何构造的基本思想。Voronoi图又叫泰森多边形或Dirichlet图, 它是由1组连接2邻点直线的垂直平分线组成的连续多边形组成。N个在平面上有区别的点, 按照最邻近原则划分平面, 每个点与它的最近邻区域相关联。Delaunay三角形是由与相邻Voronoi多边形共享1条边的相关点连接而成的三角形。Delaunay三角形的外接圆圆心是与三角形相关的Voronoi多边形的1个顶点。Voronoi三角形是Delaunay图的偶图^[8]。

Delaunay三角形网的通用算法-逐点插入算法  基于散点建立数字地面模型, 常采用在d维的欧几里得空间中构造Delaunay三角形网的通用算法-逐点插入算法。具体算法过程包括:①遍历所有散点, 求出点集的包容盒, 得到作为点集凸壳的初始三角形并放入三角形链表; ②将点集中的散点依次插入, 在三角形链表中找出其外接圆包含插入点的三角形(称为该点的影响三角形), 删除影响三角形的公共边, 将插入点同影响三角形的全部顶点连接起来, 从而完成一个点在Delaunay三角形链表中的插入; ③根据优化准则对局部新形成的三角形进行优化(如互换对角线等), 将形成的三角形放入Delaunay三角形链表; ④循环执行上述第2步, 直到所有散点插入完毕。

3.3   线框模型布尔运算

线框模型布尔运算在复杂地质体建模中主要是对有相交关系的复杂地质体(实体与实体之间、面与实体之间)进行交、差、并等, 基本的布尔运算功能对实体进行交接拼合以得到建模所需的组合实体模型。

目前大多数数字矿山工程软件都可以实现实体模型的并、交和差运算, 表面模型可以进行差、并运算。在建模实际中, 按照软件提供的布尔运算功能进行组合, 能够满足复杂地质体的建模要求。复杂地质体建模中主要遇到的问题是地质体的分支复合和地质体被断层穿插切割等地质现象很难翔实地反映在可视化模型中。在建模中, 实际上2个实体剖面之间是按照直线的方式来进行三角网连接的, 但是遇到呈曲面的断层或者断层破碎带, 实际上这种地质体与断层实体之间的吻合关系是很难在建模实体中反映的。所以只有通过线框模型布尔运算得到的实体组合来逼近地质体实体间的切割和相交等空间形态和关系。

4.   1: 5万水泥厂幅三维地质模型的建立

对于水泥厂幅三维地质模型的建立主要有3方面内容, 数字地形模型(DTM)的生成以及地下地质体模型的构建, 进而将二者无缝结合一体化显示。

4.1   数字地形模型(DTM)

数字地形模型主要由栅格和TIN技术2种数据格式表示。该地质体地表模型的构建采用了TIN技术, 因为TIN模型可以根据地形的具体特点和复杂程度而确定采样点的密度和位置, 可以避免数据的冗余, 并且可以相对准确的描述出地形的结构和一些局部地区, 同时又便于进行地形计算和分析, 能够较好地反映实际地形信息^[9]。具体步骤包括:①将1: 5万水泥厂幅地质地形图在Mapgis中进行矢量化(*.wl), 并压缩存盘。操作:打开Mapgis-输入编辑-矢量化-压缩存盘; ②把矢量化生成的Mapgis线文件(*.wl)转换输出为Mapgis明码格式文件(*.wal)。操作:图形处理-文件转换-装入线-输出Mapgis明码格式; ③将Mapgis明码格式文件(*.wal)导入到Micromine中, 转换成Micromine中的线文件(*. str)。操作:在Micromine菜单栏中, 文件-导入-Mapgis文件-线串; ④根据地形图等高线值为线串的每条等高线进行高程赋值。操作:视图管理器-线(加载)-属性窗口-Z值; ⑤对等高线线文件进行坐标矫正。操作:测量-坐标系统转换-平面; ⑥生成DTM, 其空间数据插值主要是通过TIN技术来实现的。操作:工具栏-创建DTM-选择线(图 2a); ⑦生成的DTM根据研究的需要可直接披覆该区域的地质图、地形图以及遥感图像等。操作:视图管理器-线框-属性-叠加图像。

图  2  水泥厂三维地质模型的建立过程

(注:左下角方位轴绿色箭头指向北(N)。)

Figure  2.  The establishment of 3D geological modeling around the Cement Plant

下载: 全尺寸图片 幻灯片

4.2   地下地质体模型

剖面图的绘制 地下地质体模型的建立主要有2种数据, 钻孔数据和剖面图数据, 这里采用的是剖面图数据。根据水泥厂地质图以及地表DTM在软件中的坐标数据, 我们将由西到东绘制12幅地质剖面图, 具体步骤包括:①在Micromine中加载已经完成的地表DTM, 进行相应设置。操作:视图管理器-线框(地表DTM)-显示选项-二维切块; 工具栏-向西看; ②设置完成后, 软件中显示出地形线, 输入第一个剖面的东坐标。操作:显示设置-范围-东2200;③新建线文件, 对该剖面地形线进行捕捉矢量化, 选择合适的底面线, 生成闭合的剖面线轮廓。操作:工具栏-新线、创建线、捕捉线上点; ④依据相关资料, 对地质图中的地质体走向、产状、断层和地质界线等详细分析和合理的解释, 完成此幅剖面图的绘制, 并且保证剖面图中相同岩性的的岩石连接成一个封闭的曲线, 这样就形成反映地质体现状的剖面图(图 2b); ⑤同理, 重复1~4步骤, 完成其余11幅地质剖面图, 最后生成一个12幅地质剖面图的空间三维线文件(图 2c)。

生成线框模型  ①加载地质剖面图线文件, 新建线框, 按照地质体的趋势, 利用地质体轮廓线连接成面技术在相邻同种地质体剖面线之间用三角网连接三维地质体表面, 直到连接完毕, 最后根据具体情况选择将地质体两端封闭终点或者封闭到点, 这样就形成了地质体的实体。操作:工具栏-创建线框-新三角网-封闭终点(图 2c); ②同理连接完成所有该区域的地质体, 根据地质体颜色对照表分别赋予合适的地质体颜色(图 2d); ③对相交的地质体以及高出地表DTM地质体进行布尔运算, 布尔运算能够使不同实体间的吻合空间关系表现的淋漓尽致。当地质体构造很复杂的情况下, 布尔运算将是一个工程量相当大的工作, 而且要求建模人员格外仔细, 每进行一个运算时必须要对比两块岩体, 判断正确的运算结果。操作:线框-布尔运算(图 2e、2f); ④将地下地质体模型线框(3D)另存为2D的地质体模型线框, 用于截取剖面图时的二维剖面显示。操作:线框属性-显示选项-二维切块。

4.3   数字地形模型(DTM)和地下地质体模型的结合

在Micromine同一个三维视图中加载地表DTM和地下地质体线框(2D和3D均加载), 这样就实现了地表三维数据与地下三维地质体无缝集成三维一体化显示(图 2g)。由此, 按需要可截取任一范围、任一方向、任一比例尺的剖面图, 进一步可以进行矿产勘查分析和地质分析(图 2h)。

5.   基本认识

本文探讨了基于Micromine系统的三维可视化建模技术方法流程, 对于复杂地质体的建模, 地质断层、褶皱等复杂地质结构的表示方面仍然是地质体可视化建模研究的重点之一。本文所建立的1: 5万水泥厂幅三维地质模型, 再现了该地区地质体的三维空间形态和接触关系。但是, 深部地质体的分布与相互关系强烈地依赖于对深部地质的认知程度。尽管如此, 三维可视化建模技术在今后的地质应用中潜力巨大。