Petrel解释二维浅地层数据

Petrel是斯伦贝谢开发的一款地质解释和建模软件，有点像地理信息系统的ArcGIS，主要用于数据分析和展示。它不是用来处理原始数据的，而是集成各种处理后的结果数据进行特征分析和目标拾取。当然，它也能读取原始数据，比如二维地层sgy数据，但是图谱较为粗糙，没有那么好看。

Petrel在地质地层解释行业内算是标杆软件，影响力比较大，用的人非常多，在B站有大量教程视频，它的官网也有很多详细的使用技巧。大部分内容讲的是怎么读入三维地震数据和钻井数据，如何解释地层和断层，如何建模等，而讲二维地震数据的内容比较少。

我们所涉及的地层数据主要浅层的和二维的，因此想看看Petrel能做到什么程度。

1、读入Bathy数据，XY坐标为ArcSenconds（经纬度秒）

（1）新建项目，设置项目的坐标系统，只能设置为投影坐标系统，如果不知道设置成什么投影坐标系统，就用项目所在的UTM Zone。

编辑

（2）导入数据，文件类型选择SEG-Y 2D toolbox。

无标题2.jpg

（3）打开Import poststack SEG-Y 2D datasets。

无标题3.jpg

点击Define按钮，弹出Byte location definitions对话框，此时会报一个警告信息，意思就是当前只显示1000个Trace，但是通过改变设置，能加载全部Trace。

点击查看原图

点击左下角的Specify按钮进行设置。总共有3个参数，涉及到Trace header和2D specific两个Tab窗口。在Trace header中，Trace header coordinate unit设置为Arc seconds（经纬度秒），DeadTraces为Set trace live; keep sample value，后者的意思就是加载全部Trace。

无标题5.jpg

在2D specific中，Zero coordinates 设置为Interpolate from neighboring traces。

无标题6.jpg

这三参数设置以后，Byte location definitions对话框的Graph就会显示Trace Header字段包括Line num、SP、CDP、Trace、X Coordinate和YCoordinate等的曲线图。SP和CDP如果不对，可以在File specific settings中重新设置。

点击查看原图

此时，Trace samples、histogram、Seismic、Line locations map窗口分别显示采样数据、振幅直方图、地层图谱、平面位置图。

点击查看原图

（4）指定输入文件的坐标系统，点击Input file CRS的Select按钮。因为输入文件的坐标是经纬度秒，因此，输入文件的坐标系统应为WGS84大地坐标。在第（1）步指定项目坐标系统只能选投影坐标系统，而在这里，坐标系统则可指定为大地坐标系统。这是搞什么鬼？看来，Petrel还是差ArcGIS一着。

无标题12.jpg

最后点击Import poststack SEG-Y 2D datasets窗口的右下角的Run按钮，地层数据就被加载进来了。如果要画地层、断层，在Seismic Interpretation菜单，点击Inert Seismic horizon、Fault，Petrel提供了强大的矢量编辑功能。

点击查看原图

（5）二维地层也可以在3D 窗口中显示，如下图。

点击查看原图

如果你有一个区域的地层测线，那么，就可以形成一个立体地层图。

点击查看原图

同一层的地层数据点就可以生成面。根据面可以生成 true stratigraphic thickness (TST)和 true vertical thickness (TVT)图。

无标题16.jpg

2、读入geospark电火花数据。

如果按照上面的方法读入geosaprk电火花数据，你会发现电火花sgy数据无论如何也不能导入。前面的步骤都是正常的，而且图谱也能预览，但是，一旦点击Run按钮，sgy数据导入不成功，哪怕指定导航文件也不行。我们捣鼓了好几天，尝试各种方法，包括选不同的Seismic文件类型，换其它Petrel版本等，最后都失败了。说明这个问题是Petrel 2018的一个bug。

经过实验，要想成功导入geosaprk数据，我们首先要用geospark后处理软件将经纬度sgy文件转成投影sgy文件，然后在Petrel导入窗口，选择文件类型为SEG-Y seismic data。

无标题17.jpg