基于Kriging和Marching cube算法的地学3维形态模拟  

3D Geological Modelling Based on Kriging and Marching Cube Algorithm

作　　者：杨鸿翼[1] 刘亮明[1] 赵义来[1] 

机构地区：[1]中南大学地学与环境工程学院计算地学研究中心,长沙410083

出　　处：《中国图象图形学报》2008年第3期531-535,共5页Journal of Image and Graphics

基　　金：教育部科学技术重点研究项目(104261);教育部博士学科点基金项目(20050533038)

摘　　要：地质构造的复杂性和勘查数据的不完备性,决定了对地质要素(地质体和场)的空间形态进行3维模拟是一个面临众多技术难题的艰巨任务。针对此问题,探讨了一种利用空间插值和等值面勾绘来建立地质要素的接近真实的3维可视化模型的方法。Kriging算法,由于其在揭示空间数据非均匀、各向异性分布规律方面具有良好表现,因而被用于对不完备地学数据进行空间插值,即先将分散的不规则分布的数据转换成规则分布的网格数据;再通过Marching cube算法实现由网格离散点构建光滑3维等值曲面;最后通过OpenGL接口实现等值曲面的计算机3维重绘。整个过程采用VC++.net结合OpenGL编程完成。应用这种方法,可实现有限钻孔取样分析数据的矿体品位3维等值面模拟和电磁测深2维反演数据的电阻率3维等值面模拟,其模拟结果接近于传统勘探方法解释的结果,而且其效率高得多,也更加直观,证明此方法是有效的和实用的。It is a tough challenge for the geologists to model 3D shape of geological factors because of the complexity of geologic architecture and the insufficiency of exploration data. This paper presents research on modeling 3D shape of geo-bodies and geo-fields through 3D interpolation and iso-surface constructing from limited and scatted exploration data. The Kriging algorithm,as its favorable for revealing asymmetrical and anisotropic distribution of data in the geo-field, is used as the 3D interpolation method to convert the insufficient and random-arrayed data set into sufficient grid data set. The iso-surfaces are auto-created from the scatter data through Marching cube algorithm. The above processes are achieved through VC ＋ ＋. net and OpenGL programming. By application of the above methods, the 3D iso-surfaces of ore grade from drill data and inversely deduced resistance from magnetotellurics surveying data are simulated. Compared with the traditional methods, this produces for us good results and can be used as more convenient and effective tools for 3D visualization in geo-exploration.

关 键 词：3维形态模拟 KRIGING MARCHING CUBE 等值面 OpenGL VC++.NET 

分 类 号：TP391.9[自动化与计算机技术—计算机应用技术] P208[自动化与计算机技术—计算机科学与技术]