1. 摘要

本文作者为视觉分析构建了一个真实海底的“虚拟世界”。在3D环境中导入底部轮廓。在该模型中，通过地震反射获得的海床地层剖面被数字化为离散点，并用克里金算法进行插值，以在每个地层中产生均匀的网格。然后在每一层构建 Delaunay三角模型，并使用钻井数据进行校准，以校正缓冲区内数据集的深度值。最后，构建的3D海底地层模型通过GPU着色引擎在每一层进行渲染。生成的“3D海底”通过一组关于3D海底及其分析结果的相互关联的实时信息层用于模拟、可视化和分析。

2. 引言

“数字海底”已成为海洋信息技术的热点。虚拟现实在计算机模拟中非常流行，因为它为参与者提供了真实的真实感。可视化技术利用计算机图形学和图像处理理论来表达隐藏在海量数据中的规律和关系。因此，向参与者展示数字海底信息最有效的方式是通过直观的方式表达和渲染海洋场景。

虚拟现实具有交互性、想象性、沉浸感、想象性的特点，具有动态虚拟环境实时展示的优势。而可视化具有直观性、分时性和规律性的特点，具有揭示数据内在规律和各种关系的优势。在将海量海洋数据的可视化集成到虚拟环境中时，很难将它们有机地结合起来。所以3D海底系统一方面完成了海底剖面的直观展示，另一方面保留了交互式漫游分析的实时功能。

由于海底是一种特殊的环境，人们无法直接观察和分析地质特征和现象。因此需要多学科的海洋调查来获取海床下的综合信息，不仅包括直接采样数据（如地下采样和钻井采样），还包括间接调查数据（如地球物理数据和遥感图像）。 此外，将这些多维数据可视化以全面了解海底一直是众多专家和科学家提出的热点问题。传统的海洋地球科学可视化主要集中在平面图、水平剖面和垂直剖面等，在这些二维平面中，很难详细读取海洋地质结构和特征。因此，随着可视化技术的发展，海底高程可以显示为2. 5 d由若干算术模型，如数字高程模型（图DEM），数字地形模型（DTM），数字斜率模型（DSM）。但此类模型只能显示海底表面的地形地貌特征，不能揭示海底内部的地质结构。海洋科学家需要共同组织、访问和查询大量累积和多样化的海洋数据集。3D 可视化可以一次包含更多信息，因此更适合呈现大型数据集。

GIS相关的可视化工作可以分为三类：可视化工具、地理模型的可视化和可视化辅助的3D GIS系统。在过去的几十年里，一些专家尝试了多种方法来研究基于虚拟现实的模拟。应用虚拟现实，开发了3D扩散场景管理子系统，并提出了一种在虚拟环境中模拟岩石破碎的新方法。为了满足不同的应用需求，发布了许多可以对3D对象的内部结构和属性进行建模的3D数据模型。此外，在3D方面的研究可视化、3D拓扑、3D分析和3D web服务也取得了长足的进步。本文作者在现有3D可视化研究成果的基础上，对3D建模方法，提出了一种利用钻井数据和地震剖面对海底地质对象进行交互式分析的应用，旨在在3D环境中直接、完整地显示海底地质结构和特征。此外，随着海洋信息技术的发展，水环境数据的多维动态可视化成为海洋研究的热点，主要体现在对海温、洋流盐度、海水密度等水环境要素信息的模拟仿真。

本文作者的目标不是简单地构建一个复制海床的虚拟环境，而是创建一个平台，该平台可以从海床下地层的变化中捕获数据，然后找到有用的方法（即用于分析和可视化）来使用和显示不同现象的信息，如图1所示. 到目前为止，作者已经基于该基础架构开发了多个应用程序，针对数据建模、仿真、数据可视化和轻量级分析进行了优化。

图1 海底三维视觉分析

由于原始格式和空间参考不一致，这些数据需要预先以统一格式和空间参考进行处理和转换，以便在3D环境中进行集成和可视化。测深数据和海底剖面数据被转换/写入( x, y, z )文本格式，可用于插值和表面建模。钻井数据转换成文本文件格式，记录经度、纬度、水深、孔径、长度、地层深度、地层描述等。这些钻井数据用于校准由测深数据和海底剖面数据构建的空间模型。地貌图，存储在ESRI中shape文件格式，可以通过处理和分析地下采样数据产生。侧扫图像数据以Geo-TIFF格式保存。所有空间数据都使用统一的空间参考（北京-54 坐标系，高斯-克鲁格投影）进行转换。

3. 轻量级3D地质建模

钻井数据可以直接揭示钻井位置的地层结构，进一步对钻井样品进行实验室实验，可以详细显示海底的物理和化学特征。由于海床地理环境复杂，钻井数据获取难度大、成本高，钻井样品的空间覆盖密度很低。如果仅根据有限的钻井资料构建海底地质空间模型，则无法准确、完整地揭示海底的内部结构和特征，尤其是在钻井样品之间的复杂区域。为了获得更多的海底地质结构和特征数据，海洋地球物理勘探获得了大量的海底剖面和地震剖面。这些地球物理勘探数据是通过间接方式获得的，需要科学家解释以反映钻井位置周围的地层信息，存在很大的不确定性。我们提出的平台利用直接钻井数据和间接地球物理数据生成 3通过低成本（轻量级）3D 地质建模方法建立海底地质对象的D模型，以提高3D海底可视化及其对无处不在环境的可视化分析的准确性和性能。

3.1 建模过程

图2 建模过程

离散点准备：对于海底剖面和地震剖面，所有解释地层都被数字化为( x，y，z ) ASCII格式并存储在文本文件中。

3D数据插值：合理的插值可以弥补原始数据的不足，获得额外的平滑数据。在输入的[Xmin,Xmax ]和[ Ymin,Ymax ]等区域空间参数中，将空间范围内的离散点读入内存，并通过预先定义的区间值进行克里金插值，得到规则的网格点。该方法主要考虑空间区域内属性值的可变性，确定有影响的距离范围，在该范围内等待插值的点的值将由采样点的值计算。克里金法是一种计算插值最佳线性无偏估计值的方法。该方法综合考虑样本的几何特征，如形状、大小、待估计区段间的空间分布等，对每个样本赋予一定的系数，进行线性、无偏和最小方差估计。

Delaunay三角剖分生成：基于Bowyer/Watson算法，可以使用网格点数据集生成Delaunay三角剖分，网格点数据集写入几何对象并插入到场景中。Delaunay三角剖分可以保证每个三角形的角度都接近正三角形，遵循三角剖分中最小内角最好的图形优化原则。此外，由不规则分布的离散点数据集生成的 Delaunay 三角剖分网络是独一无二的。

模型标定：控制点数据集由钻孔数据和通过循环缓冲算法计算得到的每个控制点的辐射面积组成。数据点的内插高度值由控制点数据集校准。假设控制点V 0( x 0 , y 0 , z 0) 和来自海底剖面V 1( x 0 , y 0 , z 1) 的插值点，则修正误差H在( 1 )中给出。

侧面构建：首先选择每一层的所有边界点。相邻层i和i +1同侧的点将逆时针写在向量V<i,i +1 >中，如图3所示，这个方向适合光渲染基于GPU，否则生成的网格是不可见的。逆时针矢量生成方法描述如下。以最小X的边界为例，将第i层边界点的数据集写入向量V<i,i +1 >按降序排列。然后将第i +1s层边界点的数据集按升序写入向量V<i,i +1 >。之后，使用V < i, i + 1 > 中的顶点生成多边形。对于地形起伏引起的凹多边形侧面，将进行曲面细分处理，以便完整显示。生成的侧面存储在几何数据结构中，并插入到场景根中进行可视化。

图3 侧面结构图

模型优化：为 Delaunay 三角剖分模型生成的每一层构建三角形条带，以管理和索引三角形单元格，可以减少I/O的频率，提高模型的效率。Strata 几何图形通过边缘折叠方法进行了简化，并且在渲染时会占用更少的系统资源。

属性赋值：属性赋值是多源数据融合和划分的一个非常重要的步骤。首先，根据地质类型在每一层的表面呈现特殊的颜色或纹理。其次，根据沉积物数据生成的沉积物分类，可以通过GPU着色器渲染特殊效果，如镜面反射效果、大理石效果等。第三，动态波浪效果可以通过GPU着色器渲染。

数据融合：为了详细观察地形地貌，侧扫图像以Geo-TIFF格式存储并附加到海底浮雕，沉积物数据以ESRI形状文件格式存储并导入模型进行渲染海底。

剖面渲染：地层剖面渲染由地震剖面完成，其中每个地层根据其岩性类型用定义的颜色或纹理进行渲染。使用这种渲染模式，可以可视化地震剖面，特别是在相交线中，以揭示海底模型的内部现象。这也有助于科学家纠正地震剖面的解释。

剖面分析：剖面分析主要是从用户定义的随机线的3D海床模型中提取剖面的空间信息。用户在图层中单击两次以定义切割线的起点和终点。然后平台将延伸切割线并使其切割底部以获得交点。设起点和终点分别为( x0 , y0 ) 和( x1 , y1 )。

3.2 可视化

本文采用Visual Studio 2010和OpenGL实现“剖面 - 钻孔”建模，开发了虚拟数字海床的原型应用。区域插值为中国山东省东营市附近典型海域的北纬38-39，东经118-119，包括多波束测深数据、侧扫声纳图像、海底剖面、钻井数据和地下采样数据。图5右侧，显示了由初始数据生成的地层剖面图。为了清楚地展示数据，效果由GPU渲染。图5左侧，显示其中钻写入数据后修正地层数据的效果。x轴线坐标是经度，y轴坐标是纬度，z轴坐标是距海底的深度。图6该平台在eclipse和OpenGLES开发的网站浏览器和智能手机应用程序上展示。

图5 可视化，左：钻孔修正；右：地层剖面渲染

图6 3D 海底UI，左图：网站浏览器展示；右：手机展示

4结论

本文作者发布了一种基于“剖面-钻探”模型的泛在3D海底水文平台，可以在3D虚拟环境下实现海底综合数据的建模、可视化和分析。该平台可用于通过多模态交互构建许多有趣的应用程序。3D或空间信息的多维可视化是一个热点和难点的话题。本文作者研究的空间模型和原型应用只是初步的研究和讨论。动态数据的可视化、断层褶皱特殊地质现象的表达等诸多关键问题需要今后研究解决。作者相信这个空间对于更广泛的领域是有价值的，特别是在一个多功能的海底监测设备中。 

原文信息 

[1] Lv, Z., Su, T. An ubiquitous 3D visual analysis platform of seabed. Multimed Tools Appl 78, 34613–34625 (2019). https://doi.org/10.1007/s11042-019-08180-1 译制