在本文中主要对速度深度模型进行各种处理和显示。如果输入的速度深度模型数据是规则数据，可用规则六面体进行剖分；对不规则散列数据，也可用规则六面体剖分，但精度不高，会出现锯齿状；一般的速度深度模型是不规则数据,所以我们这里用Delaunay剖分形成四面体比较精确。速度深度模型的三维可视化流程如下：

速度深度模型的三维可视化流程

    
    七、地震波波前曲面可视化

    在做叠前深度偏移处理时，其中有一项非常重要的工作就是计算地震波走时。为了检验地震波走时计算是否正确，需要把计算结果显示出来。胜利石油管理局与中国科学院地球物理所共同承担的九.五国家自然科学基金会课题“复杂地质体精细速度分析及三维叠前深度偏移理论和方法研究――以桩西埕岛古潜山为例”中就有多种计算走时的方法。中科院杨长春博士基于三维射线跟踪的微变网格方法和波前重建思路，提出了一种适用于三维复杂介质地震波走时的快速算法，可高效地获得三维叠前深度偏移的地震波理论走时；方正茂高级工程师根据Format原理用简单微积分技术和计算排序在矩形网格上计算初至旅行时，弥补了Schneider（1992）算法在实现方法上的不足。在本文的研究过程中，与复杂地质体项目紧密结合，并根据项目的要求，把计算得到的地震波波前曲面实现了三维可视化。

    由于地震波波前曲面传播的三维数据体是规则数据，其数据格式和三维地震数据体的格式相同，所以，地震波波前曲面三维可视化的方法与三维地震数据体的三维可视化基本相同。

    图8是一个速度为2800米/秒的均匀介质地震波前传播曲面图。从图中可以看出，波前面是均匀的半球面。图9是用ZX地区速度深度模型基于Format原理用简单微积分技术和计算排序在矩形网格上计算初至旅行时的地震波波前走时图，它的计算网格是120*500*120米。从图中可以看出，由于速度的变化，波前曲面有变化，并且随着与震源点的距离的增加波前的间隔也随之增大，这与实际情况一致。图10是在一个屏幕上同时显示连续不同时刻的波前曲面切片，这样便于分析。

图8 一个速度为2800米/秒的均匀介质地震波前传播曲面图

    致 谢

    感谢在中科院地球物理所刘洪研究员的精心指导和关怀；感谢胜利石油管理局的余大祥、李进龙、匡 斌、徐兆涛在技术上提供的宝贵指导和帮助。

    参 考 文 献

    1．陈永府、张华等，任意曲面的三角形网格划分,计算机辅助设计与图形学学报，Vol.9,No.5,sep.,1997

    2. 刘洪、李幼铭、高红伟、杨新民等,复杂地质体三维地质模型建立及显示，中科院地球物理所,1998

    3. 卢朝阳、吴成柯,任意多边形内带特征约束的散列数据的最优三角剖分，计算机辅助设计与图形学学报，Vol.9,No.4,Jul.,1997

    4. 张剑秋、张福炎，地球物理勘探可视化工作的挑战与机遇，石油地球物理勘探， Vol.32（6）,1997.12, 884-888

    5. 张剑秋、张福炎，地震层位信息三维可视化,石油物理勘探，Vol.33（1）,1998,2 119-124

    6. Cline A.K. and Renka,R.J.1990,A Constrained two- Dimensional triangulation and solution of closest problems in the presence of barries: SIAMJ. Numer. Anal., 25,1305—1321

    7．Donald Henrn, and M.Pauline Baker, Computer Graphics C Version , Prentice Hall International Inc.

    8. Watson, D.F., 1981, Computing the n—Dimensional Delaunay tessellation with application to Voronoi Polytopes: The computer Journal, 24（2）,167—172