蒋平(嫩江县国土资源局)
摘要:国外地质三维建模和可视化研究发展较快,我国现有的地理信息系统(GIS)都主要表达二维的地表地物的图形和属性信息,要扩展到真三维包含地下地质结构的地质信息系统还有差距。本文主要针对三维建模与可视化的关键技术进行了论述,希望对同行业起到抛砖引玉的作用。
关键词:工程地质三维建模可视化技术
0引言
纵观国内外几种软件的研究与开发现状,它们为工程地质三维建模与可视化打下了很好的技术基础,提供了很宝贵的开发经验。但是,对于工程地质专业的地质体建模与可视化分析的针对性不强,不能够很好地满足工程地质生产与研究的专业功能需要。因此本文将从分析工程地质的三维建模和可视化的关键技术问题入手,简单描述作者在工程地质三维建模和可视化方面的初步开发研究成果。
一个大型地质工程项目从可行性研究阶段、初步设计阶段到详细设计阶段,乃至到工程施工与运行阶段,往往积累了大量的地质资料,用三维模型图形图像来表达和解释如此庞大的资料,比光靠数据库和图表图纸等传统手段来得有效的多。建立工程地质体的三维模型,处理岩层界面与结构面组合关系,逼真反映地下主要地质结构全貌,将为工程地质工作者分析研究工程地质现象和发现掌握岩土体结构规律,提供一种崭新的研究手段和研究方法。
1离散数据的插值与拟合
工程地质复杂地质体中的各种地质信息,包括地表地形、地下水位、地层界面、断层、节理、风化带分布、侵入体及各种地球物理、地球化学、岩土体的物理力学参数或数据的等值面(线)等,都可以看作是三维空间中的函数,它们的拟合函数要根据实际勘测数据建立,实测数据越丰富,越能够真实描绘出这些信息的空间分布规律。地表地形测量数据、地下水位埋深测量信息等的单值曲面图形生成可归结为双自变量离散数据的插值和拟合,多值曲面如倒转褶皱和空间等值面等,则应采用多参变量插值等其他一些较复杂的方法。空间曲面插值函数有以下构造方法,如与距离成反比的加权方法(Shepard方法),径向基函数插值法(Multiquadric方法),平面弹性理论插值法等,它们同样适用于单个连续地层界面、地球物理勘探数据、地球化学勘探数据以及岩土体物理力学参数在地质体空间的分布。
2三维数据结构
工程地质体一般是不规则形体,在计算机图形学中曲线和曲面总是分别通过很多微小直线段和微小三角面逼近来模拟地层岩性界线和岩层曲面,即岩层界面(和地表曲线、地下水位面等地质层面界线)和岩层曲面都分别是许多微小直线段和微小三角面的集合。地质体三维空间数据结构是工程地质三维建模和可视化的基础,这就要求必须具备有效的分层的三维数据结构,能够确保人机交互和查询的实现。
3曲面求交
地质体中存在大量各种层面,当出现地层不整合、断层错断岩层、地层尖灭和地下水出露于河谷地表等情形时,就自然会遇到曲面间求交的问题;地质体三维模型的上部边界是地表曲面,通过数学方法拟合出的岩层面或地下水位面不应超出地表曲面,即超出部分不应显示。同样的,当显示多层地层时,下面的每一岩层应以其上一岩层为边界。因此,为了可视化地层界面必须要解决地层面与地表、断层面和其他地层面的求交问题。另一方面,在剖面图成图时,地质界线的绘制是通过显示剖面(平面)与各种地质界面(曲面)求交所得出的交线。因此曲面求交包括地质界面(层面)之间的相交,和地质界面与剖面的相交两类问题。
4三维拓扑结构分析
从地质学角度看,拓扑是地质对象间关系的表格,拓扑表存储层位间上覆、下伏和交切(被断层切割后地层的拓扑表达)等的地层学关系及地质空间位置关系。拓扑也可视为允许这些地质关系合理储存的数据结构。例如,考虑多层地层,上一个岩层的底面和与其相邻的下一个岩层的顶面是上下岩层这两个实体的公共部分或共享边界,它们之间的拓扑关系就是相邻和同一的关系,在存储数据时只存储上一个岩层的底面或其相邻的下一个岩层的顶面,即相邻岩层的边界曲面可以存为一个地层曲面,大大减少数据存储量。评价地质模型系统的优缺点往往决定于描述地质对象所用的拓扑结构。
5可视化技术
工程地质复杂地质体可视化,是利用计算机技术将工程勘测获得的数据,转换为形象直观的便于进行交互分析的地下地质结构空间形态的立体图和剖面图形,其基础是工程数据和测量数据的可视化。利用可视化技术可以从庞大的地质勘测数据中构造出地质工程中对于边破稳定性和地下硐室变形破坏等起关键作用的岩层和结构面,并显示其范围、走向和相互交切关系,帮助工程地质人员对原始数据做出正确解释,继而为工程地质分析具体问题提供决策支持。
6三维建模与可视化技术的初步应用
6.1地质工程复杂地质体三维建模与可视化的研究框图
地质工程复杂地质体三维建模与可视化的过程是,从勘探数据库中提取各种地质信息的坐标位置及岩土体的物理力学参数,通过不同的拟合与插值函数得到地质层面(曲面)和地质实体的三维计算机图形显示,表达地质信息在研究区域内的分布规律。生成地质岩层面和地质实体后,实现从任意角度观察建立的模型,实现根据指定的剖面走向、倾向和倾角生成垂直剖面。
6.2初步开发与应用
工程勘测空间数据库在收集整理现场勘测数据后录入各分项数据表,这些数据表不仅包括地质信息的位置数据,更重要的是提供属性数据。以地层岩性数据表为例,要求录入钻孔编号、岩层起始深度、岩层终止深度、层厚、岩性(地层名称)、地层代码(地层年代)、岩层走向、岩层倾向、岩层倾角、接触关系、地质描述等数据。随着工程勘测的进展,能够方便地修改补充和管理勘测数据。