三维地质建模方法
2020-01-20 · 技术研发知识服务融合发展。
自20世纪80年代以来,研究人员提出了许多三维地质模型来模拟地质体,使这方面的研究有了长足的发展。通过对国内外大量的三维地质建模方面的文献和专业软件的研究分析,三维地质建模方法大体可归纳为三类:离散点源法、剖面框架法和多源数据耦合建模法。
1.2.3.1 离散点源法
在地质找矿中,经常需要根据少量的离散点采样数据(如地质测绘或钻孔资料)来获取地质体的形状,从而为进一步指导找矿起指导性的作用。因此,研究如何实现空间散乱点数据场可视化的方法具有一定的意义。
Carlson(1987)从地质学的角度提出了地下空间结构的三维概念模型,并提出用单纯复形模型(Simplicial Complex Model)来建立地质模型。Victor(1993)、Pilout(1994)则具体应用Delaunay四面体的三维矢量数据模型研究离散点地质建模问题。Lattuada(1995)对3DDT(3 Dimensional Delaunay Triangulation)在地质领域内的应用进行了研究,表明四面体格网能很好地用于地质体的三维建模,优点包括:四面体单元易于建立索引;模型易于手工编辑;可通过相邻关系导出拓扑结构;约束三角剖分易于实现面约束;四面体非常便于可视化,同时具有较高的表达精度;易于实现搜索和关系查询等。Courrioux et al.(2001)基于Voronoi图实现了地质对象实体的自动重构。Frank et al.(2007)采用隐函数法(implicit function)表达三维曲面,对离散点集进行三维重构(reconstruction),用来模拟断层和盐丘(salt dome)。杨钦(2001,2005)利用离散点源信息构建地层与断层结构面,依此作为约束条件约束Delaunay剖分建立三维地质模型。
钻孔数据也属于一种点源信息。它实质上是将原始的点、线数据进行有效的分层,根据各层面标高应用曲面构造法来生成各个层面或实体。围绕钻孔数据进行三维地质建模已有许多学者进行了研究,其中较早利用钻孔数据进行三维地质模拟的是加拿大学者Houlding(1994,2000),利用钻孔孔口点位信息进行 Delaunay三角剖分,作为“主 TIN(Primary TIN)”,其他地层面则通过高程映射实现。张煜等(2001)对其建模方法进行了深入研究与发展,在垂直钻孔的理想状态下,采用三棱柱(Tri-prism,TP)数据模型建立三维地质模型,并给出了相关的剖切算法。Lemon et al.(2003)采用“地层层位法”建立三维地层模型,并采用自定义剖面(user-defined cross-sections)的方法对地质模型进行局部交互修正。吴江斌(2003)、朱合华等(2003)提出一种基于钻孔数据的二分拓扑数据结构的建模算法,尝试采用基于钻孔数据的四面体体元模型构建地下三维地质模型;四面体结构在表达复杂结构上则较灵活,但是使用四面体表示空间实体会产生大量的冗余,且生成四面体的算法比较复杂。张芳(2005)采用Delaunay三角构网技术,利用钻孔数据构建三维地层层面模型,同时引入“界面分片”思想,以适应于海量数据模型的可视化表达,但缺少对地质体属性信息的表达。在三棱柱模型的基础上,针对钻孔存在偏斜问题,类三棱柱(Analogical Tri-prism,ATP)(齐安文等,2002)、广义三棱柱(Generalized Tri-prism,GTP)(Wu,2004)方法先后被提出,用来进行三维地质建模,已被证明广泛适应于矿山、石油等深部地质问题建模;同时,似三棱柱(Similar Triprism,STP)(Gong et al.,2004)也被提出用于解决钻孔倾斜问题,如郑蔚等(2005)基于钻孔数据采用STP建立三维地质模型对地下空间进行虚拟漫游。STP与GTP本质上是相同的。基于钻孔数据建立三维地质建模,这一看似简单的数据模型方法,经历了10多年的发展历程:从初期的TP数据模型,适用于钻孔垂直成层、地层等厚的理想情况,发展到STP、GTP适用于钻孔不垂直且地层不等厚的常见情况。
1.2.3.2 剖面框架法
剖面框架法就是在收集整理原始地质勘探资料的基础上,建立分类数据库,人工交互生成大量的二维地质剖面,然后应用曲面构造法生成各层位面表达三维地质模型,或者利用体元表示法直接进行地质体建模(Chae et al.,1999)。
利用地质剖面表达研究区域三维地质现象的初级形式是序列地质剖面法(朱小弟等,2001)。序列地质剖面构模技术实质是传统地质制图方法的计算机实现,即通过平面图或剖面图来描述地质构造,记录地质信息,如图1.2所示。其特点是将3D问题2D化,在空间上采用若干平行或近似平行的地质剖面来表达研究区域的地质分布特征,但它在空间表达上是不完整的,它把剖面之间的地层或构造分布情况留给工程设计人员去“想象”。这种构模方法难以完整表达3D矿床及其内部构造。
基于剖面信息建立真三维模型具有很大的发展空间,对于复杂地质构造区域具有很好的适应性,成为当前地质建模的主要方法之一。然而,基于剖面进行三维重构得到完善发展的是在医学领域,后来迅速扩展到其他领域。在医学领域里,通过电脑断层扫描(CAT)或者核磁共振(MRI)等技术,可以获得一系列相互平行的人体切片图像,通过提取对象的边界,基于轮廓线算法,生成三维人体模型。地质剖面信息同医学切片信息一样,都是反映研究对象的某一特定断面上的构造分布,可以借助医学三维人体建模技术来构造三维地质模型。较早将医学领域的切面三维建模引入地学领域的是在考古学方面的应用(Tipper,1976,1977;Herbert et al.,1995),主要应用在古生物的三维重构方面,而应用在三维地质建模方面的文献并不很多。
图1.2 序列地质剖面构模实例
公认的剖面三维重构的代表之作是Keppel的文章(Meyres et al.,1992;Herbert et al.,1995,2001;Xu et al.,2003;屈红刚等,2003)。在Keppel的研究基础上,Meyres(1992)将剖面建模方法分为4个子问题:对应问题(correspondence problem)、构网问题(tiling problem)、分支问题(branching problem)和光滑问题(fitting problem):①对应问题解决相邻剖面之间的轮廓线匹配问题;②构网问题主要解决轮廓线之间的三角形构网问题,考虑满足某个准则,例如最大体积法(Keppel,1975)、最小面积法(Fuchsetal.,1977)等;③分支问题是解决同一对象在不同剖面上的组成部分的个数不同的问题;④光滑问题主要解决将初始生成的三角网进行插值,从而得到更加光滑的三角网。
屈红刚等(2003)提出基于含拓扑剖面地质建模方法来实现复杂地质的三维建模的对应问题,邓飞等(2007)则对一般意义上的剖面地质建模进行了讨论。
1.2.3.3 多源数据耦合建模法
随着计算机性能的提高,具备了对海量数据的处理能力,人们对建立的地质模型要求也不断提高,希望能够建立高精度和高复杂度的地质模型(Turner,2003,2006;Calcagno et al.,2006;Kaufmann et al.,2008)。提高模型的精度可以通过插值的方法来实现,但更好的方法是通过增加约束信息来对初始地质模型进行细化,这就涉及耦合多源数据来建立地质模型的问题。
早在1993年,Houlding提出三维地学建模概念的时候就强调地质解释信息具备对模型的修正(revision)功能。并且指出矿业工程有地质勘探数据、人工绘制数据及施工数据,还有不确定性的需要通过地质统计学进行估计的数据(Houlding,2000),最终的地质模型需要综合考虑这些种类不同的数据。
McInerney et al.(2005a,b)认为三维地质建模只能部分上是一个数字地质采样过程,更重要的是地质学家的人工解释过程。并且尖锐地指出,不要指望一些计算机软件能够自动并成功地“建模”! 让一个有经验的地质学家输入解释性的信息进行建模,是现实和必要的;而软件只是建模过程中提供便利的一个工具(There is no expectation that some computer software will successfully and automatically“builda model”! The reality is that interpretative input from a skilled geologist is essential to build a model;the software is simply a tool to facilitate the model-building process)。其要求实际上是,地质建模不仅要考虑地质勘探所获取的确定性数据,还应加入地质工程人员对地质构造的解释性数据,这就构成多源地质建模的基本思想。
Mallet(2002)针对地质体建模的特殊性和复杂性,以点、线数据为主要数据源,建立以三角形为基本单元的三维曲面,采用离散光滑插值技术(Discrete Smooth Interploate,DSI)使曲面的粗糙度最小,并作为GOCAD的核心技术,得到了许多地球物理公司和石油公司的支持。
相比较国外以石油、矿业工程为主要应用领域的三维地质建模,钟登华等(2006)则从水利水电工程地质领域,研究多源地质数据建立坝区的三维地质模型。Wu et al.(2005)提出一种逐步细分的多源数据集成地质建模方法,考虑到地质数据大多比较稀疏和低采样率的特征,采用逐步细化的方法对初始地质模型不断修正。
地质构造的复杂性和认识的阶段性,使多源地质建模引起越来越多的研究兴趣。32届国际地质大会(International Geological Conference,IGC)于2004年在意大利佛罗伦萨召开,在“地质的复兴(The Renaissance of Geology)”(Zanchi et al.,2007)议题上,多名国际知名的地学建模专家共同提到了多源地质建模问题。其中,Zanchi et al.(2008)借助商业软件对意大利境内阿尔卑斯山(Alps)利用多源地质建模问题进行研究,并应用于滑坡稳定性分析。西方发达国家主要将地质建模应用于能源与环境领域,这是为数不多的在工程建设领域开辟蹊径的研究。无独有偶,Kaufmann et al.(2008)尝试采用多源地质建模,研究在废弃煤矿巷道内进行天然气储存问题。
总体来看,三维地质建模技术是一个从简单地层模拟到复杂地质构造模拟的发展过程。从最初基于单一数据建立简单层状三维地质模型,到综合利用多源数据建立复杂地质模型,能够反映地质构造的空间特征。