基于IDL的地质体真三维重构与可视化研究+——以福州市地热田热储结构为例.pdfVIP

基于IDL的地质体真三维重构与可视化研究 ——以福州市地热田热储结构为例 唐丽玉1 陈崇成i 王钦敏1 缪海岚2 (I．福州大学数据挖掘与信息共享教育部重点实验室，福建省空问信息工程研究中心，福建福州350002；2．联想软件 有限公司，北京100085) 摘要：通过对钻孔线状数据进行有效的分层．直接根据各个层面的标高，应用不规则三角网(11N)表面法 和克里格(kriging)内插算法构建三维体数据，从而实现地质三维重构。并基于IDL开发环境，采用光线投射 法和zbuffer相结合的直接体绘制算法，实现地质三维分析和可视化。同时研发了福州市地热田热储结构系统， 实现了地质真三维交互式切割分析、显示和查询。 关键词：三维重构，可视化，三维交互式切割，IDL 中图分类号： 文献标识类型：A 1 引言 地质三维重构是实现地质可视化的瓶颈[1]。地质体是个复杂的空间实体，很难找到一个统一的三 维重构方法。目前常用的构造方法可归纳为剖面成面法、直接点面法和拓扑分析法【l】。Sirakov等 数据[5]。目前拓扑分析基本采用六面体、四面体模型[1]。地质三维显示是实现地质可视化的关键。目 Volume method)和直接体绘制(Direct 前可视化方法可分为两大类：表面拟合方法(Surface-fitting Cubes)求出子区域的边界多边形并进行排序后，再进行图象合成，其计算速度较传统体元投射法快3—4 Elisa 倍；Ana 型的需求。 本文采用直接点面法，基本思路是，直接将原始的线状数据进行有效的分层，根据各个层面的标高， 应用曲面构造法来生成各个层面，再把各层面组合成三维体数据。可视化采用光线投射法与z_buffer相 (InteractiveData Language)是RESEARCHSYSTEMS公司开发的可视化编程工具、交互式数据语言。 利用它可以灵活方便地分析任何数据，通过交互的2D和3D数据分析及图像处理，可以实现各种分析数 据的可视化，交互的3D分析能力不仅有OFENGL的支持，采用面向对象的技术，它对数据及统计软件包 提供各种科学计算及分析方法[12，13]。 2 地质三维重构 地质三维重构是个非常复杂的问题，因为地质本身是一个三维性、非均质性非常明显的复杂体。实 测资料分布极不均匀和不全面，具有很强的个性特点。但是地质三维重构是解决地质真三维可视化的瓶 颈。下面介绍基于IDL开发的福州市地热田热储结构系统所用的地质三维重构的方法。 2．1数据预处理 本研究区域为福州市地热田集中分布区，有69个钻孔资料。根据热储结构特点，把地层结构基本 概化为：不透水盖层、第四系泥质砂砾卵石层状热水岩组、基岩脉状裂隙热水岩组、不透水基底岩层。 在未揭示到不透水基底岩层的钻孔，采用基岩脉状裂隙热水岩组下推100米作为两者的界面。由数据库 读出的数据整理后格式如表1。 表l钻井结构数据示例 Thec8．$eofborestructuredata Tab．1 表中岩层属性：l表示不透水盖层，2表示第四系泥质砂砾卵石层状热水岩组，3表示基岩脉状裂隙 热水岩组，4表示不透水基底岩层。 地质三维重构 要迸行直接体绘制，首先要进行基于栅格结构的三维重构，即将有限的离散数据在足够重建精度下 方体内的值是变化的，体素内任一点的值可以用八个顶点上采样值的三次线性插值计算出来。 蘑 图1 以体素为构成单位的体结构 Volumestructuremadeofvoxel Fig．1 单凭过于稀少的采样数据无法实现基于栅格结构的三维体绘制，甚至运用插值算法，也会对其内插 获得的数据的客观性有很大的影响。从离散稀疏的采样数据到生成可构造体素的三维矩