su简单中文使用手册.docx

su简单中文使用手册.docx

《su简单中文使用手册.docx》由会员分享，可在线阅读，更多相关《su简单中文使用手册.docx（39页珍藏版）》请在冰豆网上搜索。

su简单中文使用手册

CWP软件的安装与简单使用手册

CWP软件的安装

一．在LINUX下建立用户CWP，在CWP下建立目录path，将源文件cwp.su.all.37.tar.Z放

二．到path目录下，并建立bin文件夹

三．在CWP用户主目录下显示隐藏文件，修改.bash_profile文件，在已有的export之后另起一行，分别添加

exportCWPROOT=/home/CWP/path，再于

PATH=$PATH:

$HOME/bin后添加

：

/home/CWP/path/bin:

/home/CWP

退出保存

四．从终端中分别输入

cdpath

zcatcwp.su.all.37.tar.Z|tar–xvf-

…

待终端中反映完毕，分别输入

cdsrc

makeinstall

makextinstall

makemglinstall

makeutils

makexminstall

makesfinstall

这期间可能有系统安装所等待的时间，不用急，但凡遇到yes/no，一路y下来即可。

四．为了检查是否安装完毕，在终端中输入

Suplane>data.su

Suxwigb

若出现一个简单的图像，则成功！

CWP软件的简单说明

一、文中涉及的命令全部以小写形式，均可在终端窗口下输入，以次获取自述帮助。

先说几个命令：

suplane和suxwigb，more。

suplane作用是产生一个简单的零偏移距su文件，suxwigb是一个典型的X—windows绘制图形工具，如例子：

suplane>data.su

suxwigb

more

比较全面的了解它们，请在终端中输入suplane,suxwigb,more。

二、关于DEMOS的应用所有DEMOS必须把文件拷到用户根目录下，而后依照readme文件中的执行顺序，在终端中输入文件名。

注意目录下的文件变动。

三、在执行DEMOS文件时，如果想清楚了解程序执行过程，请输入

moreprogramname

由于水平有限，这里的谬误很多，希望大家能在偶尔翻看时，多多留心，发现并改正，衷心希望能和大家一起学习。

谢谢

第一节两种数据的转换

在CWP应用中免不了和两种数据打交道，su和segy格式。

它们有联系也有区别。

一．数据的输入输出

1读写编辑数据

常用命令如下：

segyread---readanSEG-Ytape

segywite---writeanSEG-Ytape

segyclean—zeoroutunsignedportionofheader

suaddhead---putheadersonbaretracesandsetthetrackandnsfields

sustrip—removetheSEGYheadersfromthetraces

supaste---pasteexitstingSEGYheadersonexistingdata

segyhdrs—makeseg-yasciiandinaryheadersforsegywrite

bhedtopar,setbhed---editingthebinaryheaderfile

surange---getmaxandminvaluesfornon-zeroheaderentries

suchw---changeheaderwordusingoneortwoheaderwordfields

sugethw---gettheheaderwordsinsudata

suedit---examinesegydiskfilesandeditheaders

suxedit---examinesegydiskfilesandeditheaders

2常规的数据转换命令

a2b—convertasciifloatstobinary

b2a---convertbinaryfloatstoascii

ftnstrip---convertfortranfloatstoc-systelsfloats

recast---recastdatatype（convertfromonedatatypetoanother）

h2b—convert8bithexadecimalfloatstobinary

transp—transposeann1byn2elementmatrix

二．SEGY与SU的转换

SEGY文件格式由三部分组成，镜象头文件，十进制头文件和实际的地震记录，而SU格式的文件只含有其第三部分。

命令segyread可以实现两种格式的转换，使用如下：

在“big-endian”平台上，

Segyreadtape=/dev/rmt0verbose=1endian=1>data.su

或者更常用一点的

Segyreadtape=/dev/rmt0verbose=1endian=1|segyclean>data.su

命令行中斜线部分是SEGY磁带的位置（或者是磁盘文件的位置），其它应用看自述文档。

三．典型命令的使用

如下命令行输入，注意主目录中文件的变化和终端中的屏幕输出。

1．

suplane>data.su

suxwigb

segyhdrs

bhedtopar

setbhedbfile=binarypar=binary.par[parameter1=valueparamenter2=value……]

segywrite

b2adata.ascii

a2bdata.binary

transpdata.trnsp

2．Su数据简单操作

suplane>data.su

sustripdata.bin

suaddheaddata.su

supastedata.su

surange

suplane|surange

sugethw

suplane|sugethwkey=tracl,tracr,offset,dtnsoutput=binary>file.bin

morefile.bin

sushwdata.new.su

suplane|sushwkey=sxa=6400c=-100j=5|sugethwkey=sx|more

suplane|sushwkey=offseta=200b=200j=5|sugethwkey=offset|more

suplane|sushwkey=dt,sx,offset,a=2000,6400,200b=0,0,200c=0,-100,0j=0,5,5|sugethwkey=dt,sx,offset|more

sueditdata.su

第二节CWP中图形打印显示工具

一．X-WINDOWS下的图形打印

1．浮点数据绘图

常用的命令有

xcontour—xcontourplotoff（x1,x2）viavectorplotcall

ximage---ximgaeplotofauniformly-sampledfunctionf（x1,x2）

xwigb---xwiggle-traceplotoff（x1,x2）viabitmap

xgraph---xgraphergraphsn[i]pairsof（x,y）coordinates

xmovie---imageoneormoreframesofauniformlysampledfunctionf（x1,x2）

相应的，在SU平台上有以下程序

suxcontour

suximage

suxwigb

suxgraph

suxmovie

suxmax

例子：

suplane|suxcontourtitle=”contuour”&

suplane|suximgaetitle=”image”&

suplane|suxwigbtitle=”wigb”&

suplane|suxgraphtitle=”graph”&

suplane|suxmovietitle=”movie”&

suplane|suxmaxtilte=”max”&

suplane|junk1.su

suplane|suaddnoisesn=20>>junk1.su

suplane|suaddnoisesn=15>>junk1.su

suplane|suaddnoisesn=10>>junk1.su

suplane|suaddnoisesn=5>>junk1.su

suplane|suaddnoisesn=3>>junk1.su

suplane|suaddnoisesn=2>>junk1.su

suplane|suaddnoisesn=1>>junk1.su

suxmovie

2．PS下的图形打印

pscontour

psimage

pswigb

psgraph

psmovie

pscube

pswigp

在SU平台上相应的程序有：

supscontour

supsimage

supswigb

supsgraph

supsmovie

supscube

supswigp

例子：

suplane>junk.su

supscontourdata1.eps&

supsimagedata2.eps&

supswigbdata3.eps&

supsgraphdata4.eps&

supsmoviedata5.eps&

supscubedata6.eps&

supswigpdata7.eps&

运行以上程序之后，进入主目录，在所生成的*.eps文件上点击右键，打开方式选择ggv，查看图像。

第三节常见处理

以下几乎每部分都有CWP自带的DEMOS，参见$CWPROOT/path/src/demos目录。

对于一般的情况来说，DEMOS的文件执行具有一定的顺序性。

一．动校正

这一部分有两个命令，sunmo和suvelan。

使用分别如下：

SUNMO-NMOforanarbitraryvelocityfunctionoftimeandCDP

命令格式为

sunmostdout[optionalparameters]

可选参数：

tnmo=0相应于NMO中速度的NMO时间

vnmo=2000相应于NMO中时间的NMO速度

anis1=0产生四次节点的两各向异性系数

anis2=0时距曲线中，相应于tnmo的时间

cdp=由tnmo和vnmo定义的CDP

smute=1.5sampleswithNMOstretchexceedingsmutearezeroed

lmute=25length（insamples）oflinearrampforstretchmute

sscale=1=1todivideoutputsamplesbyNMOstretchfactor

invert=0=1toperform（approximate）inverseNMO

ixoffset=0donotconsidercross-lineoffset

=1readcross-lineoffsetfromtraceheader

upward=0=1toscanupwardtofindfirstsampletokill

注意:

在处理常速度反射层时，vnmo取单一常数就可以了；在处理多层变速度反射层时，tnmo和vnmo总是成对出现,且数据为相对应的数组,如”tnmo=1,2…vnmo=1500,2000…”。

关于NMO的DEMO，进入$CWPROOT/path/src/demos/NMO，操作如下：

1．拷贝NMO下的所有文件到主目录下

2．依次输入xmodel－makedata－xdata－xvclaity－xnmo－clean－xall；

3．依次输入psmodel－makedata－psdata－psvclaity－psnmo－clean－psall。

请注意：

在源码中可能有错误，执行之中有些内容无法实现。

二．叠加

有以下命令：

sustack、sudivstack，说明如下:

SUSTACK-stackadjacenttraceshavingthesamekeyheaderword

sustackstdout[Optionalparameters]

Requiredparameters:

none

Optionalparameters:

key=cdpheaderkeywordtostackon

normpow=1.0eachsampleisdividedbythe

normpow'thnumberofnon-zerovalues

stacked（normpow=0selectsnodivision）

verbose=0verbose=1echosinformation

SUDIVSTACK-DiversityStackingusingeitheraveragepowerorpeak

powerwithinwindows

Requiredparameters:

None

Optionalparameters:

key=tracfkeyheaderwordtostackon

winlen=0.064windowlengthinseconds.

typicalchoices:

0.064,0.128,0.256,

0.512,1.024,2.048,4.096

ifnotspecifiedtheentiretraceisused

peak=1peakpoweroptiondefaultisaveragepower

examples：

Forduplicatefieldshotrecords:

susortstack.data

ForCDPordereddata:

sudivstackstack.data

关于DEMOS，进入$CWPROOT/path/src/demos/Stacking?

Phase_Weighted_Stacking，输入supws文件，观察主目录下的文件变更情况。

三．抽道集

常用命令，susort，使用说明：

SUSORT-sortonanysegyheaderkeywords

susortstdout[[+-]key1[+-]key2...]

Susortsupportsanynumberof（secondary）keyswitheitherascending（+,thedefault）ordescending（-）directionsforeach.Thedefaultsortkeyiscdp.

Note:

Onlythefollowingtypesofinput/outputaresupported

Diskinput-->anyoutput

Pipeinput-->Diskoutput

Note:

IfthetheCWP_TMPDIRenvironmentvariableissetuse

itsvalueforthepath;elseusetmpfile（）

Example:

Tosorttracesbycdpgatherandwithineachgatherbyoffsetwithbothsortsinascendingorder:

susortOUTDATAcdpoffset

关于DEMOS，进入$CWPROOT/path/src/demos/sorting_trace/tutorial输入Xsort，PSsort文件即可。

由于以上三部分密不可分，特以速度分析一部分为例，

1、产生数据

nx=100fx=0dx=80

nz=110fz=0dz=50

unif2

v00=5000.0,6000.0,8000.0,10000.0,15000.0|

ximagen1=$nzf1=$fzd1=$dzn2=$nxf2=$fxd2=$dxwindowtitle="Model"\

legend=1units="wavespeed"\

wbox=$WIDTHhbox=$HEIGHTxbox=$WIDTHOFFybox=$HEIGHTOFF&

2、产生数据

REF1="0.0909091:

-4.0,1.0;12.0,1.0"

REF2="0.1428570:

-4.0,2.2;12.0,2.2"

REF3="0.1111110:

-4.0,3.5;12.0,3.5"

REF4="0.2000000:

-4.0,5.0;12.0,5.0"

fpeak=25.0#shouldbehalfnyquisttoavoidallspatialaliasing

dvdz=2.0#velocitygradient（dv/dz）

dvdx=0.0#velocitygradient（dv/dx）

v00=5.0#P-wavevelocityatsurfaceinkft/s

ob=0#toeliminatetheturnedraycontribution

nxo=64fxo=0.1dxo=0.1#offsetinformation:

number,first,spacing

nxs=12fxs=1.4dxs=-0.1#shotinformation:

number,first,spacing

nt=501dt=0.004#timeinformation:

number,spacing

echo"Makingsyntheticdatafordemo"

susynlvv00=$v00fpeak=$fpeakndpfz=10dvdz=$dvdzob=$ob\

dt=$dtnt=$ntverbose=0\

ref=$REF1ref=$REF2ref=$REF3ref=$REF4\

nxo=$nxofxo=$fxodxo=$dxo\

nxs=$nxsfxs=$fxsdxs=$dxs|

suchwkey1=cdpkey2=gxkey3=sxb=1c=1d=2|

suaddnoisesn=50f1=4.0f2=8.0f3=20.0f4=25.0>modeldata

2、动校正

sunmo

tnmo=0,0.34,0.625,0.875,1.13vnmo=5000,6000,7050,8410,10000tnmo=0,0.335,0.630,0.875,1.1vnmo=5000,6000,7180,8350,9760tnmo=0,0.35,0.64,0.88,1.1vnmo=5000,6000,7180,8380,9590

tnmo=0,0.64,0.88,1.1vnmo=5000,7200,8440,9740>nmodata&

3、NMO图

suwind

sugaintpow=2gpow=0.5|

suximagelabel1="Time"label2="Gather"\

title="Every5thShotRecordAfterNMO"\

windowtitle="NMO"legend=1units="amplitude"\

f2=1d2=0.078125f2num=1d2num=5\

wbox=$WIDTHhbox=$HEIGHTxbox=$WIDTHOFF3ybox=$HEIGHTOFF&

4、抽道集，叠加

susort

sustacknormpow=1.0>stackdata

5、抽道集，叠加绘图

sugain

suximagelabel1="Time"label2="CMP"title="Stack"windowtitle="Stack"\

f2=350d2=50legend