SeismicUnix简介学习教程.pptx

《SeismicUnix简介学习教程.pptx》由会员分享，可在线阅读，更多相关《SeismicUnix简介学习教程.pptx（92页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、目录3帮助工具4核心SU程序5绘图6常见处理2关于SU7SU与SHELL语言1SU安装8附录第1页/共92页帮助工具 一、前期准备工作：刚装上ubuntu后，首先需要将网络设置好，并选择合适的软件源。对于教育网来说，推荐用中国科大的软件源。在百度上均可找到软件源的代码。1)首先备份Ubuntu12.04源列表，终端中输入一下命令：(备份下当前的源列表)2)修改更新源 (打开Ubuntu12.04源列表文件)注：软件源的设置应该根据电脑上所装ubuntu版本而定。这里以ubuntu12.04版本的安装为准。第2页/共92页帮助工具 3)将中科大源的代码粘贴进去(可直接放在文件中的前面)：4)通知

2、ubuntu启用新的更新源：终端中键入 至此，源的更新工作已经完毕。注：软件源的设置应该根据电脑上所装ubuntu版本而定。这里以ubuntu12.04版本的安装为准。第3页/共92页帮助工具 在终端中直接用 sudo apt-get install XXX-XXX命令安装相应的软件包。我主要安装了如下的包：build-essential gcc g77 lesstif2 lesstif2-dev libx11-dev libglu1-mesa ligblu1-mesa-dev xlibmesa-glu freeglut3 freeglut3.dev libxmu6 libxmu-dev li

3、bxi6 libxi-dev libxt-dev libglut3-dev gfortran 以安装build-essential为例:在终端中键入sudo apt-get install build-essetial命令，若出现如下画面，表示安装正常启动，稍后的工作就是等待安装结束，进行下一个软件包的安装:后面的软件包的安装也依旧如此，只需在终端中输入相应的安装命令即可(见上面绿体字)。第4页/共92页帮助工具 1)创建su目录：2)在home目录下，用vi或gedit打开.bashrc文件(终端下命令格式 gedit.bashrc):3)在文件最后加上以下三行(这里的yourname=hu

4、angxr)：export CWPROOT=/home/yourname/su#设置seismic unix的根目录。export PATH=$PATH:/home/yourname/su/bin#设置seismic unix的路径。export XAPPLRESDIR=/home/yourname/su/lib/X11/app-defaults#这是README_X中提示要设置的环境变量。第5页/共92页帮助工具 4)保存退出。在终端下运行：source.bashrc 完成环境设置生效。echo$CWPROOT 检查CWPROOT设置是否正确。echo$PATH 检查CWPROOT设置是否正

5、确。如果环境变量设置没问题，那么终端中的界面同下面的一样：至此，系统环境设置完毕。第6页/共92页帮助工具 二、安装Seismic Unix 1.下载Seismic Unix 从下载当前的SU版本。我下载的是cwp_su_all_43R3.tgz(进入上面的网站，点击下面界面红色方框部分的内容即可下载)。具体存放在是/home/yourname/su下：第7页/共92页帮助工具 在终端中进入/home/yourname/su目录，然后键入tar-zxvf cwp_su_all_41_tar.Z 解压文件包：在/home/yourname/su/下生成/src目录。SU软件的相关文件都在其下。安

6、装说明Installation_Instructions和安装配置文件，需要认真阅读。第8页/共92页帮助工具 *3.修改makefile.config文件(这是以前的SU版本需要进行的步骤)进入/home/yourname/su/src 这个目录，可以看到makefile.config 这个文件，要修改里面的三项,输入命令：gedit makefile.config,然后修改所打开的文件：ENDIANFLAG=-DCWP_LITTLE_ENDIAN 句激活（即去掉#）include$(CWPROOT)/src/Rules/gnumake.rules 句激活，并屏蔽掉 include$(CWP

7、ROOT)/src/Rules/oldmake.rules 句,修改后 保存。XDRFLAG=-DSUXDR 句激活（即去掉#）LARGE_FILE_FLAG=-D_FILE_OFFSET_BITS=64-D_LARGEFILE_SOURCE 句激活(即去掉#，否则seismic unix处理数据有大小限制)OPTC=-O3-Wall-Wno-long-long-ansi-pedantic-D_POSIX_SOURCE 句激活(即去掉#，这一项要与cc编译系统相配)CC=cc OPTF=-O-static-march=i686-fno-automatic-fno-second-undersco

8、re 这一项要与fortran编译系统相配 FC=g77,FC=f77,FC=ifort 屏蔽掉，激活FC=gfortran#现在的新的SU版本中的makefile.config该文件已经修改过了，因此不需要进行*3.修改makefile.config文件这一步。第9页/共92页帮助工具 此时在进行正式安装SU之前，建议先看一下/home/yourname/su/src目录下的安装说明Installation_Instructions文件，确保前面的所有工作都已经正确完成(因为新版本的SU的发布，可能会有一些变化，导致它在安装上也会出现有一些变化，如果此时完全照着本文的安装步骤进行可能会出错)

9、。第10页/共92页帮助工具 根据Installation_Instructions 提示，在/home/yourname/su/src目录下分步进行安装：make install (to install the basic set of codes)make xtinstall(to install the X-toolkit applications)make finstall(to install the Fortran codes)make mglinstall(to install the Mesa/Open GL items)先以make install为例，展示其安装时的部分界面

10、：现在终端键入上述的命令:。随之安装开始。若安装结束后出现下边的画面，则表示这一步的安装成功：其它几步的安装也同上，在此不必一一说明。第11页/共92页帮助工具 为了检测安装SU是否成功，在终端中键入下面的管程命令：suplane|suxwigb&如果出现下图，就表示安装成功！第12页/共92页目录3帮助工具4核心SU程序5绘图6常见处理2关于SU7SU与SHELL语言1SU安装8附录第13页/共92页关于SU Seismic Unix地震数据处理系统（简称SU软件包）由科罗拉多矿业学院(Colorado School of Mines,CSM)波场处理中心(the Center for Wa

11、ve Phenomena,CWP)开发的处理程序。此软件是在基于Unix的操作系统上创造一个地震处理环境(用C语言编写)来扩展Unix操作系统来进行地震处理和研究。而且，该软件的所有源码免费提供给想用的人。SU软件包不仅限于进行地震处理。用SU软件可以作与波场相关的一切大套处理，在一定意义上这是一个含有比“Seismic”(地震)一词更广义的软件。SU目的是作为Unix操作系统的一个扩充，因此它兼有Unix系统的很多特性，包括Unix的灵活性和扩展性。就像Unix系统操作一样，SU可以堪称是一种语言(或者元语言)SU是什么第14页/共92页关于SU SU软件包目前没有使用功能驱动的图形用户界面

12、。它是Unix操作系统的一个延伸，一些SU的图形界面只能是有限地访问软件包中的大部分功能。SU不是商业地震软件包的替代品，但是SU能够作为使用的任何一种商业软件包的重要助手。商业软件包用于研究工作，SU常用作一个原型软件包(a prototyping package)。还有如果需要写进新的编码，SU能够作为新的软件应用的开始的基础。SU并不限于地震应用。它可能用于勘探地球物理和通常的信号处理应用中。它也一定能够用于向学生介绍关于“与波有关”的信号处理(特别是傅立叶变换)教学中。这包括雷达、非地震声波和图像资料处理等。再有，SU既不能说是一个3D软件包，也不能说是一个2D软件包，因为其中有很多的

13、滤波和到操作任务与2D和3D中一样。SU不是什么第15页/共92页关于SU 通过访问CWP主页：可以下载最新版本的SU软件包，也可以通过该网站了解更多的CWP/SU的内容。因为SU包含有很多程序，刚开始接触没有那么容易。这里推荐读者可以查阅书籍Seismic Unix地震数据处理系统-SU3.3处理操作系统指南。这本书好比一本“字典”，可以很好地回答SU中一些不可避免的问题。参考资料第16页/共92页目录3帮助工具4核心SU程序5绘图6常见处理2关于SU7SU与SHELL语言1SU安装8附录第17页/共92页帮助工具 有很多工具按不同的详细程度对软件包中的主程序、shell文本以及库函数提供了

14、内部文档：SUHELP-列出CWP/SU程序shell文本。SUNAME-自述文档中各项的名字行以及源码的位置。The sefdoc-是关于大多数可执行主程序和shell文本的一个内部文档工具。通过在命令行键入不带参数、没有用管城（|）或Unix重定向符()进行输入和输出重定向的程序名和shell文本命令，就可以看到自述文档内容。SUDOC-得到编码的DOC列表。SUFIND-自述文档里得到的信息。SUKEYWORD-segy.h 文件中SU关键字说明。Demos-SU中包含了一套demos，位于目录$CWPROOT/src/demos下，主要是shell脚本。$CWPROOT/src/su/

15、main下，SU程序源码。$CWPROOT/src/su/examples下，相关的shell程序。在此，只介绍一下sukeyword、The selfdoc、Demos。其它帮助工具请见附录-帮助工具。第18页/共92页帮助工具 使用方法：在Linux终端窗口命令行中输入不带参数或重定向到文件的程序命令。用途：简单说明程序的用途。举例：sustack#sustack的自述文档段落就打印（显示）到屏幕中 The selfdoc对该程序的用途进行了简短说明 指出如何在命令行键入程序命令，关键字“stdin”和“stdout”分别表示输入取自标准输入和标准输出。Unix这些术语的意思是用户能够从磁

16、盘文件输入和输出参数，或者用Unix重定向输入命令符“”，或管程“|”输入和输出参数。第19页/共92页帮助工具 说明：查看SU道头域中的关键字。使用方法：sukeyword name#name-道头关键字名称 举例：若要查看关键字“tracl”，在终端中键入:sukeyword tracl#tracl-侧线中道顺序号#当从模拟程序或从野外数据中提取数据时，了解上面这个集合的内容是有好处的。SUKEYWORD第20页/共92页帮助工具 说明：SU软件包中有一套演示程序，它们是shell脚本，放在目录$CWPROOT/src/demos中。它们可以用来验证SU程序。建议：如果是刚学SU软件包，可

17、以按以下的顺序运行演示程序：Making_Data目录中有用susynlv程序制作合成记录炮点集和共偏移距剖面的基础内容演示流程。Filtering/Sufilter目录中有一些实际数据处理中消除地滚波和初至演示程序。Deconvolution目录中演示了使用supef和其它工具简单合成脉冲道集，来举例说明去混响和脉冲反褶积处理。演示程序包括使用loops系统检验滤波参数影响的命令。Demos-演示程序第21页/共92页帮助工具 Sorting_Traces指南是一个交互文本，加强了文献中讨论的一些UNIX和SU基本知识。下一步是激活Selecting_Traces演示程序。再接着是NMO演示

18、程序。除此之外，像SYNLV、Migaration等里面也包含重要的演示程序，也应该值得演示。当更加熟习SU的逻辑知识后，就能更自如地按源码运行任何其它演示程序，或者用户按自己的目的修改shell文本来演示程序！第22页/共92页目录3帮助工具4核心SU程序5绘图6常见处理2关于SU7SU与SHELL语言1SU安装8附录第23页/共92页核心SU程序 Seismic Unix的核心程序可以处理多种多样的任务，这些任务可以看作是关于数据研究和处理的。这些任务包括：输入/输出。数据格式转换。设置、查看和编辑道头。浏览查看SU数据。开时窗显示、抽道集和编辑数据。一般操作。变换和滤波操作。对SU数据的

19、地震处理操作。在这一章节中主要介绍前四部分的内容。第24页/共92页核心SU程序 输入/输出输入/输出 1.SEGY格式和SU数据格式：在20实际80年代早期，最通用的数据存储格式是SEG-Y。这是美国勘探地球物理学家协会(Society of Exploration Geophysicists,SEG)的Y格式。这种格式现在仍然广泛使用。SEGY数据格式由三部分组成。第一部分是3200字节的EBCDIC卡片头，包括40个卡片对该磁带进行描述的文本数据。第二部分是400个字节的二进制（卷）头，描述了该磁带卷上内容的相关信息。第三部是实际的地震道数据。每道含有240个字节的道头，数据记录在道头后

20、面。SU数据格式基于SEGY格式的道部分。SEGY和SU道的主要不同在于SU格式的道数据是浮点型，使用运行SU程序的本机上的二进制浮点格式写进。SU数据只由SEGY的道数据组成！SU格式中不保存EBCDIC和二进制卷头。用segyread程序可以把SEGY数据转成SU程序能用的格式。第25页/共92页核心SU程序 SEGY格式SU格式3200个字节ebcid编码区：40条记录80字节=3200字节400字节binary二进制数区域卷头文件道头文件：道头文件：240个字节道头信息；每道数据第26页/共92页核心SU程序 2.SEGYREAD-将SEG-Y数据读入SU 说明：程序segyread用

21、来把数据从SEGY格式转换成SU格式。用法：segyread tape=*segy verbose=1 endian=0 conv=0|segyclean*.su tape后面可以跟文件路径，如 tape=/home/zc/*.sgy 参数：endian=0 or 1#=0 表示“小头”；=1 表示“大头”注意事项：1）当读SEGY磁带或数据文件时，需要使用机器的位序(byter-order)(即头，endian)。SGI、SUN、IBM和所有Motorola芯片机使用的是“大头”(big-endian)或高位(high-byte)IEEE格式，而基于Intel和Dec芯片的系统用的是“小头”

22、(little-endian)或低位(low-byte)；2）在SEGY道头里有可选的道头字(181-240字节)。这些道头字的使用没有标准。SU也不例外。有些SU图形程序使用的参数就存储在这些道头字节中。segyclean程序用于把这些可选的道头字节信息清零，使得SU图形程序不至于被这些信息混淆。#更多关于segyread的参数信息请键入自述文档查看。第27页/共92页核心SU程序 3.SEGYHDRS-为segywrite创建SEGY ascii码和二进制道头 说明：segyhdrs命令创建“header”(道头)和“binary”(二进制)文件。用法：举例：利用suplane命令制作一些

23、测试数据，然后运行segyhdrs命令。注意：“header”（道头）文件是一个ascii码文件。4.SEGYWRITE-写SEGY磁带或磁盘文件 说明：与segyread相联系的命令segywrite程序，可以将SEGY格式按多种不同的方式把数据写到磁带或磁盘文件。segywrite创建文件时需要有header部分和binary部分。如果你没有binary和header文件，你必须用上面提到的程序sgyhdrs来创建它们。最后通过segywrite命令来写磁带，总的命令是：segyhdrs *.su|segywrite tape=*.segy verbose=1 第28页/共92页核心SU程

24、序 输入/输出数据格式转换 常常需要对其它系统的数据进行转换，或需要输入不同格式的数据。在SU里有很多工具和方法处理这些问题。下面列出几个相关的程序：A2B-把ascii码浮点型转换换成二进制码。B2A-把二进制浮点型转换换成ascii码。FTNSTRIP-把Fortran浮点型转换成C浮点型。FTNUNSTRIP-把C浮点型转换成Fortran浮点型。H2B-把8位16进制浮点型转换成二进制。RECAST-改变数据类型(从一种数据类型转到另一种)。SWAPBYTES-不同的数据类型(高、低位)字节(BYTE)交换。SUSWAPBYTES-SU数据字节序从“大头”(big endian)交换为

25、“小头”(little endian)，或者是反过来。下面列举几个常用的程序。第29页/共92页核心SU程序 1.A2B和B2A-ASCII到Binary、Binary到ASCII 说明：在所有的数据格式中，ASCII码最具移植形(也最耗空间)。任何系统均能输入或输出ASCII码数据，还有文本编辑器支持ASCII码。而二进制文件却不能访问。当需编辑二进制码文件时，需要用到A2B程序，将ASCII转换成Binary。反之，需用到b2a，可查看或编辑ASCII码文件。用法：a2b *.bianry 相反的操作是：b2a *.ascii 注意：ASCII码文件可以用通常的文本编辑器进行编辑。只要满足

26、总共个40行、每行80个字符的格式要求，可以在文件中放任何内容。因此上面的n1=number 中的number最好别太小(代表列数)。否则超过40行，ASCII码文件不能正常显示。第30页/共92页核心SU程序 2.FTNSTRIP-Fortran数据输入到C 说明：Fortran语言中二进制数据用“记录头”(BOR)和“记录尾”(EOR)分隔符分开。而C程序语言创建的二进制数据没有这些分隔符。要在C程序中使用Fortran数据需要去掉这些Fortran标识符。用法：3.RECAST-改变二进制数据类型 说明：C支持很多种数据类型，recast程序能对大量的数据类型集进行转换处理。举例：把整型

27、转换成浮点型 用法：可以改变的数据类型：float-浮点型;double-双精度;int-(带符号)整型;short-短整型;char-字符;uchar-无符号字符;long-长整型;ulong-无符号长整型.第31页/共92页核心SU程序 输入/输出设置、查看和编辑道头 在SU数据(SEGY格式例子后面)这些参数的值储存在道头中。有很多程序可以访问这些道头，因此可以根据不同目的设置、浏览和修改这些道头字。SUADDHEAD-把道头加到空道(bare traces)上，设置tracl和ns道头字。SUSTRIP-从道中去掉SEGY道头。SUPASTE-已存在的数据粘贴上已存在的SEGY道头。S

28、URANGE-获取非零道头输入值中最大值和最小值。SUSHW-设置一个或多个道头字，利用道数、取余数(mod)或整除(integer divide)计算道头值或从一个文件中输入道头字。SUCHW-利用一个或两个道头关键字改变道头关键字。SUGETHEW-获取SU数据中道头关键字。SUEDIT-检查segy磁盘文件并编辑道头。SUXEDIT-检查segy磁盘文件并编辑道头第32页/共92页核心SU程序 1.SUADDHEAD-加SU(SEGY风格)道头 说明：如果数据格式正确，也就是C浮点型排列中快速空间方向为每道采样点数增加方向，就可以给这些数据加上道头，以便其它SU程序访问。用法：suadd

29、head *.su 例1：加道头到含有C二进制浮点型数据文件中,用surange检查数据。从上图可以看到,shot.dat文件已经成功从一般的地震数据格式转换成SU数据格式的shot.su。相关的SU程序surange、suedit可以访问该文件.第33页/共92页核心SU程序 例2：对于其它的类型(如整形)使用命令recast 例3：如果数据首先是从Fortran转换而来的整型，那么处理流程为第34页/共92页核心SU程序 2.SUSTRIP-将SU道头与SU数据分离 说明：命令sustrip是命令suaddhead的逆操作。用法：上面的命令系列将去掉文件SU道头，并把SU道头存在文件dat

30、a.headers 3.SUPASTE-将SU道头粘贴到二进制数据 说明：对二进制数据完成处理操作后，可能要把道头粘贴回去，这时用命令supaste。用法：or suspate data.su heade=headers ns=#注意：在用supaste命令时，千万不要把“ns=”参数忘记，否则“粘贴”不成功。第35页/共92页核心SU程序 举例：用之前由正演模拟转换来的数据“shot.dat”上联合使用sustrip和supaste命令。第36页/共92页核心SU程序 4.SURANGE-获得道头字值的范围 说明：对于道头字一条有用的信息是看指定数据集中道头字数值范围，此时可用surange

31、程序。用法：susrange *.su 参数：traces=201#共200道 ns=4000#每道采样点为4000#注意：对于破坏的数据来说，大量的道头中可能显示出奇怪的数值。检测这种问题是命令surange的主要用法之一。第37页/共92页核心SU程序 5.SUGETHW-在SU数据中获取道头关键字的值 说明：前面的surange命令可以查看整个数据道头字表示的最大和最小值。但是，我们尝尝要一道一道查看道头字的数值，并按我们选择的顺序。命令sugethw就是这样的一个工具。用法：surange data.su key=keyword1,keyword2,|more 举例：第38页/共92页

32、核心SU程序 对于“观测系统设置”(geometry setting)，可以按照下面的命令操作：可以看到，此时屏幕上没有显示道头字，是因为 的缘故。第39页/共92页核心SU程序 6.SUSHW-在SU数据中设定道头字 说明：命令sushw可根据各种目的来设置地震道头值。该命令可让用户设置一个或更多的道头字。它常用来对一个特定道头字设置一个值。用法：sushw data.su key=key1,key2,key3 a=a1,a2,a3 参数：key=cdp，设置的道头关键字 a=0，第一道头的值 b=0，组内递增量 c=0，组增量 d=0，道头序号增减 j=ULONG_MAX,ULONG_MA

33、X，组中的元素的数量 第40页/共92页核心SU程序 例1：地震数据“shot.su”中没有表示“dt”(采样间隔)道头字。通过sushw命令给“shot.su”加上道头字“dt”。并键入surange检查。第41页/共92页核心SU程序 对于a,b,c,d,j这些额外的可选参数可用来做更复杂的操作。道头字的值用以下公式计算：1）i=itr+d 2）val(key)=a+b*(i%j)+c*(i/j)这里itr是道号(第一道是itr=0,而不是1)。上面百分号(%)表示取余函数(modulo)，/表示除法。例2：1)设定道头字dt=2000(按2ms采样；2)设定5道中的道头字sx=6400，

34、第二个5道中sx=6300，后面是每5道递减(-100；3)设定每5道数据集道头字offset的值为200、400、6400 步骤1）：或者 第42页/共92页核心SU程序 步骤2）：或者 第43页/共92页核心SU程序 步骤3）：或者 第44页/共92页核心SU程序 当然，可以值引用一次sushw命令就完成所有3中操作：或者 第45页/共92页核心SU程序 补充：补充：在前面的key=dt中，为什么大小为2毫秒的dt其值要用“a=2000”表示，而不是a=2或a=0.002。此时可以键入sukeyword查看。通过上面补充信息，我们可以进一步了解sukeyword的作用。可见，熟习地运用su

35、keyword程序，对于准确地查找相应关键字的信息具有很大的帮助作用！dtdt用微妙用微妙做单位做单位第46页/共92页核心SU程序 7.SUCHW-在SU数据中改变(或计算)道头关键字 说明：有些道头字(如cdp)可以用已有的道头字计算而来，程序suchw就提供这种功能。用法：suchw data.su key=key1,key2,key3 a=a1,a2,a3 参数：key1=cdp，输出的关键字 key2=cdp，输入的关键字 key3=cdp，输入的关键字 a=0，整体偏移一常量 b=1，对第一个输入道头字进行乘除运算 c=0，对第二个输入道头字进行乘除运算 d=1，整体进行乘除运算

36、第47页/共92页核心SU程序 这个程序使用两个 道头字的值(即key2和key3)，利用下面的方程式来计算第三个道头字(key1)：val(key1)=(a+b*val(key2)+c*val(key3)/d 例1）道头字cdp的值偏移一常量(比如-1)suchw outdata a=-1 例2）增加一个常量(比如1000)道头字(比如说ns)suchw key1=ns key2=ns a=100 outfile 例3）设置道头字gx，用sushw命令将偏移距(offset)和sx(shot point)道头字的值相加，然后对sx和gx取平均来计算cdp道头字值。第48页/共92页核心SU程

37、序 8.SUEDIT和SUXEDIT-编辑SU数据中的道头关键字 说明：suedit和suxedit可以检查或改变某个道头。用法：suedit data.su（suxedit区别不大）举例：第49页/共92页目录3帮助工具4核心SU程序5绘图6常见处理2关于SU7SU与SHELL语言1SU安装8附录第50页/共92页绘图 Seismic Unix软件包有小部分(普通C浮点型格式或SU格式)数据绘图的图形工具，可在X-Windows环境下在屏幕显示或PostScript格式硬拷贝。在SU中可能的绘图类型有：contour plots(等值线绘图)gray or colorscale image

38、plots(灰色或彩色图像绘图)wiggle trace plots(地震道波形绘图)line or symbol graph(侧线或符号图)movies(电影)3D cube plots(PostScript only)(3D立体图绘图)(只限于PostScript格式)这些程序自述文档较长，有大量的可选项用来选择图形的外观和标注。然而，取视窗选择数据的功能，应该用程序subset或suwind进行数据处理后才送到绘图程序。第51页/共92页绘图 输入/输出X-Windows绘图程序 在X-Windows环境下，用于浏览通用浮点型数据(没有SU道头)的程序：XCONTOUR-用矢量绘图程序对

39、函数f(x1,x2)的X等值线绘图。XIMAGE-相同采样点函数f(x1,x2)的X图像绘图。XWIGB-用位图对函数f(x1,x2)的X波形地震道绘图。XGRAPH-ni对坐标(x,y)的X GRAPHer图形。XMOVIE-相同采样点函数f(x1,x2)的一幅(帧)或几幅(帧)图像。第52页/共92页绘图 举例：用suplane生成二进制数据，数据大小为n1=64n2=32。然后从SU数据中分离处道头。接着用上面的绘图程序来查看数据(这里只列举出了xcontour，其它的绘图程序的写法均可仿照xcontour，这里就不再逐一列出)说明：符号&是一个Unix命令，告诉工作的shell在后台运

40、行该程序。当程序执行完后，可键入字幕q退出。第53页/共92页绘图 为了SU格式数据绘图，创建了很多程序。它们是：SUXCONTOUR-用矢量绘图程序对地震Unix到文件X等值线绘图。SUXIMAGE-SU数据X-Windows图像绘图。SUXWIGB-SU数据X-Windows位图地震波形绘图。SUXGRAPH-SU数据X-Windows图形绘图。SUXMOVIE-SU数据X电影绘图。SUXMAX-SU数据每道最大值、最小值和最大绝对值X-Windows图。举例：用suplane命令生成的数据对上述各个程序进行测试 suplane|suxcontou title=contour&suplan

41、e|suximage title=image&其它的SU绘图程序均可仿照上面的suxcontour和suximage，这里不再逐一列出。第54页/共92页绘图 输入/输出PostScript绘图程序 可用与通常浮点型数据(数据没有SU道头)PostScript绘图程序有：PSCONTOUR-两个空间函数f(x1,x2)PostScript等值线绘图。PSIMAGE-等间隔函数f(x1,x2)PostScript图像绘图。PSWIGB-用位图对函数f(x1,x2)的PostScript地震道绘图。PSGRAPH-ni对坐标(x,y)PostScript图形。PSMOVIE-等间隔采样函数f(x1

42、,x2,x3)PostScript电影绘图。PSCUBE-数据三维数据体PostScript图像绘图。PSWIGP-用多边形对函数f(x1,x2)的PostScript地震道绘图。举例：用前面去掉道头的data.bin(n1=64,n2=32)数据 在终端中键入gs data.eps可查看该图片。第55页/共92页绘图 这些程序有：SUPSCONTOUR-SU数据PostScript等值线绘图。SUPSIMAGE-SU数据PostScript图像绘图。SUPSCUBE-SU数据PostScrip立体绘图。SUPSWIGB-SU数据PostScrip位图地震波形绘图。SUPSGRAPH-SU数据

43、PostScrip图形绘图。SUPSWIGP-SU数据PostScrip多边形填充波形绘图。SUPSMAX-SU数据每道最大值、最小值和最大绝对值PostScrip绘图。举例：类同于前面的通常浮点型数据PostScript绘图，这里不再叙述。第56页/共92页绘图 在SU中还有几种工具可用于PostScript操作，它们是：PSBBOX-改变存在的PostScript文件的四周的框架(BoudingBOX)。PSMERGE-合并PostScript文件。MERGE2-把多张PostScrip图合并到一页上。MERGE4-把多张图合并到一页上MERGE4 figures onto one pag

44、e。这些操作在这里就不详细阐述。可以键入下面的命令获得更多的信息：more$CWPROOT/bin/merge2 或者下面的命令：more$CWPROOT/bin/merge4 另外，在终端中键入gs*.eps,可查看所有的“*.eps”类型的图片。第57页/共92页绘图 最后展示2张由有限差分正演模拟得来的地震记录的结果。左边一张是一般地震数据（无道头）格式的shot.dat文件，由ximage得到的结果。右边一张是转换得来的加上道头的SU数据格式的shot.su文件，由suximage的得到结果。第58页/共92页目录3帮助工具4核心SU程序5绘图6常见处理2关于SU7SU与SHELL语言

45、1SU安装8附录第59页/共92页常见处理 在第三章的“核心程序”中，主要介绍了一些道头操作的命令。这一章中，将对一些常用的SU数据编辑操作、通用操作、地震模拟工具、地震处理工具做简单介绍。这些部分几乎都有CWP自带的DEMOS,参见$CWPROOT/src/demos目录。对于一般情况，DEMOS的文件执行具有一定的顺序性。第60页/共92页常见处理 输入/输出编辑SU数据 一旦数据读入和道头设置正确后，常常要对数据进行操作和编辑。这些操作任务常常不同，可以进行下面数据处理：Windowed(取时窗)Sorted(道分选)Truncated(道截取)Tapered(斜坡化处理)Zeroed(

46、道充零)作等间隔采样 Concatenated(数据联合)在这里，只简单介绍取视窗的操作suwind。其它处理的相关程序见附录 编辑SU数据。第61页/共92页常见处理 SUWIND-用关键字选定一定时窗地震道 说明：很常见的操作是浏览和处理地震数据中的一小块。suwind命令可以用于选定地震数据窗口大小。参数：key=tracl#设置窗口的关键道头字 min=LONG_MIN#数据关键道头字最小值 max=LONG_MAX#数据关键道头字最大值 count=number#计数参数 j=1#每个数据的第j道不选 s=0#以s为基数(如果(key-s)%j)=0)reject=none#按指定的

47、关键字跳过的地震道 accept=none#用指定的关键字选择地震道 第62页/共92页常见处理 例1：用道数关键字对命令suplane产生的数据进行窗口大小的设置 第63页/共92页常见处理 例2：每个2道选取数据，基数为1 第64页/共92页常见处理 例3：用reject拒绝第3，8，9道。用accept接受第8道。第65页/共92页常见处理 输入/输出SU数据通用操作 SU数据编辑以外的操作任务是进行下面的处理：gaining(增益)resampling(重采样)unary operaions 一元操作(对单个文件的算术操作)。binary operations 二元操作(对两个文件的算

48、术操作)。进行这些处理的程序有：SUADDNOISE-对地震道加噪声 SUGAIN-按不同的增益方式显示地震道 SUOP-在segy数据做一元算术操作 SUOP2-对两个数据做二元操作 在这里，不做详细介绍。具体可参考其自述文档“the selfdoc”。第66页/共92页常见处理 输入/输出变换和滤波操作 地震研究和地震处理的很多方面与数学变换方法有关。特别地，很多地震处理离不开数字Fourier变换。滤波是一个与此相关的主题，因为大多数滤波在频率域进行，或者至少是可以用数学上进行频率域操作来表示。1.1D Fourier变换 SUFFT-从实数时间道变换到复数频率道(正变换)。SUIFFT

49、-从复数频率道变换到实数时间道(逆变换)。SUAMP-从(频率，x)数据输出振幅(amp)、相位(amp)、道实部和虚部。SUFFT-从实数时间道变换到复数频率道(正变换)。SUSPECFX-地震数据(时间T到频率F)Fourier频谱分析。第67页/共92页常见处理 例1：suplane|suxigb title=Original Data&#左图 suplane|sufft|suifft|sushw key=d1,dt a=0,4000|suxwigb&#右图#结果与输入完全一样，除了地震道结果中有更多的采样点，因为变换是要补零(zero-padding)第68页/共92页常见处理 例2：

50、浏览sufit程序输出的振幅谱和相位谱、及实部和虚部 suplane|sufft|suamp mode=amp|suxwigb title=amplitude&suplane|sufft|suamp mode=phase|suxwigb title=phases&suplane|sufft|suamp mode=real|suxwigb title=real&suplane|sufft|suamp mode=imag|suxwigb title=imaginary&例3：SU数据格式可以存储复数数据的实部和虚部。键入下面的命令查看：例4：用suspecfx命令快速查看一个地震道或一块地震数据的