Petrel是斯伦贝谢开发的一款地质解释和建模软件，有点像地理信息系统的ArcGIS，主要用于数据分析和展示。它不是用来处理原始数据的，而是集成各种处理后的结果数据进行特征分析和目标拾取。当然，它也能读取原始数据，比如二维地层sgy数据，但是图谱较为粗糙，没有那么好看。
Petrel在地质地层解释行业内算是标杆软件，影响力比较大，用的人非常多，在B站有大量教程视频，它的官网也有很多详细的使用技巧。大部分内容讲的是怎么读入三维地震数据和钻井数据，如何解释地层和断层，如何建模等，而讲二维地震数据的内容比较少。
我们所涉及的地层数据主要浅层的和二维的，因此想看看Petrel能做到什么程度。
1、读入Bathy数据，XY坐标为ArcSenconds（经纬度秒）
（1）新建项目，设置项目的坐标系统，只能设置为投影坐标系统，如果不知道设置成什么投影坐标系统，就用项目所在的UTM Zone。
（2）导入数据，文件类型选择SEG-Y 2D toolbox。
（3）打开Import poststack SEG-Y 2D datasets。
点击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。
在2D specific中，Zero coordinates 设置为Interpolate from neighboring traces。
这三参数设置以后，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一着。
最后点击Import poststack SEG-Y 2D datasets窗口的右下角的Run按钮，地层数据就被加载进来了。如果要画地层、断层，在Seismic Interpretation菜单，点击Inert Seismic horizon、Fault，Petrel提供了强大的矢量编辑功能。
（5）二维地层也可以在3D 窗口中显示，如下图。
如果你有一个区域的地层测线，那么，就可以形成一个立体地层图。
同一层的地层数据点就可以生成面。根据面可以生成 true stratigraphic thickness (TST)和 true vertical thickness (TVT)图。
2、读入geospark电火花数据。
如果按照上面的方法读入geosaprk电火花数据，你会发现电火花sgy数据无论如何也不能导入。前面的步骤都是正常的，而且图谱也能预览，但是，一旦点击Run按钮，sgy数据导入不成功，哪怕指定导航文件也不行。我们捣鼓了好几天，尝试各种方法，包括选不同的Seismic文件类型，换其它Petrel版本等，最后都失败了。说明这个问题是Petrel 2018的一个bug。
经过实验，要想成功导入geosaprk数据，我们首先要用geospark后处理软件将经纬度sgy文件转成投影sgy文件，然后在Petrel导入窗口，选择文件类型为SEG-Y seismic data。
电火花sgy数据的图谱显示如下图，这是未处理的原始数据，经过一些简单的滤波操作，就会变成我们通常看到的清晰地层。我们在Petrel中找了好久，没有发现它有任何滤波功能。
3、小结
与ArcMap一样，Petrel的优势在于成果数据的集成、分析和展绘。对于一般浅地层工程项目，它的作用不大，但是对于大型的浅地层研究项目，它不仅能将收集的各种矢量和栅格数据集成三维地质系统，而且能分析和展示测区的整体地层特征和特殊地质现象，将显著提高了研究报告的质量和深度。
参考文献
1、https://wiki.iodp.tamu.edu/display/LMUG/Importing+Bathy2010+sub-bottom+profile+into+Petrel
发表评论：