面朝大海,春暖花开

作者:xiaokcehui2023-8-11 13:37分类: subbottom 标签: Petral segy

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

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

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

1、读入Bathy数据,XY坐标为ArcSenconds(经纬度秒)

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

无标题1.jpg

(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

电火花sgy数据的图谱显示如下图,这是未处理的原始数据,经过一些简单的滤波操作,就会变成我们通常看到的清晰地层。我们在Petrel中找了好久,没有发现它有任何滤波功能。

点击查看原图

3、小结

与ArcMap一样,Petrel的优势在于成果数据的集成、分析和展绘。对于一般浅地层工程项目,它的作用不大,但是对于大型的浅地层研究项目,它不仅能将收集的各种矢量和栅格数据集成三维地质系统,而且能分析和展示测区的整体地层特征和特殊地质现象,将显著提高了研究报告的质量和深度。



参考文献

1、https://wiki.iodp.tamu.edu/display/LMUG/Importing+Bathy2010+sub-bottom+profile+into+Petrel

温馨提示如有转载或引用以上内容之必要,敬请将本文链接作为出处标注,谢谢合作!

已有 0/8590 人参与

发表评论:

欢迎使用手机扫描访问本站,还可以关注微信哦~