su简单中文使用手册.doc

CWP软件的安装与简单使用手册CWP软件的安装一 在LINUX下建立用户CWP，在CWP下建立目录path，将源文件cwp.su.all.37.tar.Z放二 到path目录下，并建立bin文件夹三 在CWP用户主目录下显示隐藏文件，修改.bash_profile 文件，在已有的export之后另起一行，分别添加export CWPROOT=/home/CWP/path，再于PATH=$PATH:$HOME/bin后添加：/home/CWP/path/bin:/home/CWP退出保存四 从终端中分别输入cd pathzcat cwp.su.all.37.tar.Z | tar xvf-待终端中反映完毕，分别输入cd srcmake installmake xtinstallmake mglinstallmake utilsmake xminstallmake sfinstall这期间可能有系统安装所等待的时间，不用急，但凡遇到yes/no，一路y下来即可。四为了检查是否安装完毕，在终端中输入Suplane data.suSuxwigb data.susuxwigb data.su more data.su或者更常用一点的Segyread tape=/dev/rmt0 verbose=1 endian=1 | segyclean data.su命令行中斜线部分是SEGY磁带的位置（或者是磁盘文件的位置），其它应用看自述文档。三典型命令的使用如下命令行输入，注意主目录中文件的变化和终端中的屏幕输出。1suplane data.susuxwigb data.su &segyhdrs data.subhedtopar binary outpar=binary.parsetbhed bfile=binary par=binary.par parameter1=value paramenter2=valuesegywrite b2a data.asciia2b data.binarytransp data.trnsp2Su数据简单操作suplane data.susustrip data.binsuaddhead data.susupaste data.susurange data.susuplane | surangesugethw file.binmore file.binsushw data.new.susuplane |sushw key=sx a=6400 c=-100 j=5 | sugethw key=sx | moresuplane | sushw key=offset a=200 b=200 j=5 | sugethw key=offset | moresuplane | sushw key=dt,sx ,offset,a=2000,6400,200 b=0,0,200 c=0,-100,0j=0,5,5 | sugethw key=dt,sx,offset | moresuedit data.su第二节 CWP中图形打印显示工具一X-WINDOWS下的图形打印1浮点数据绘图常用的命令有xcontourx contour plot of f(x1,x2) via vector plot callximage-x imgae plot of a uniformly-sampled function f(x1,x2)xwigb-x wiggle-trace plot of f(x1,x2) via bitmapxgraph-x grapher graphs ni pairs of (x,y) coordinatesxmovie -image one or more frames of a uniformly sampled function f(x1,x2)相应的，在SU平台上有以下程序suxcontoursuximagesuxwigbsuxgraphsuxmoviesuxmax例子：suplane | suxcontour title=”contuour” &suplane | suximgae title=”image” &suplane | suxwigb title=”wigb” &suplane | suxgraph title=”graph” &suplane | suxmovie title=”movie” &suplane | suxmax tilte=”max” &suplane | junk1.susuplane | suaddnoise sn=20 junk1.susuplane | suaddnoise sn=15 junk1.susuplane | suaddnoise sn=10 junk1.susuplane | suaddnoise sn=5 junk1.susuplane | suaddnoise sn=3 junk1.susuplane | suaddnoise sn=2 junk1.susuplane | suaddnoise sn=1 junk1.susuxmovie junk.susupscontour data1.eps &supsimage data2.eps &supswigb data3.eps &supsgraph data4.eps &supsmovie data5.eps &supscube data6.eps &supswigp data7.eps &运行以上程序之后，进入主目录，在所生成的*.eps文件上点击右键，打开方式选择ggv，查看图像。第三节 常见处理以下几乎每部分都有CWP自带的DEMOS，参见$CWPROOT/path/src/demos目录。对于一般的情况来说，DEMOS的文件执行具有一定的顺序性。一动校正这一部分有两个命令，sunmo和suvelan。使用分别如下：SUNMO - NMO for an arbitrary velocity function of time and CDP 命令格式为 sunmo stdout optional parameters 可选参数： tnmo=0 相应于NMO中速度的NMO时间 vnmo=2000 相应于NMO中时间的NMO速度 anis1=0 产生四次节点的两各向异性系数 anis2=0 时距曲线中，相应于tnmo的时间 cdp= 由tnmo和vnmo定义的CDP smute=1.5 samples with NMO stretch exceeding smute are zeroed lmute=25 length (in samples) of linear ramp for stretch mute sscale=1 =1 to divide output samples by NMO stretch factor invert=0 =1 to perform (approximate) inverse NMO ixoffset=0 do not consider cross-line offset =1 read cross-line offset from trace headerupward=0 =1 to scan upward to find first sample to kill注意:在处理常速度反射层时，vnmo取单一常数就可以了；在处理多层变速度反射层时，tnmo和vnmo总是成对出现,且数据为相对应的数组,如”tnmo=1,2 vnmo=1500,2000”。关于NMO的DEMO，进入$CWPROOT/path/src/demos/NMO，操作如下：1 拷贝NMO下的所有文件到主目录下2 依次输入 xmodelmake dataxdataxvclaityxnmocleanxall；3 依次输入 psmodelmake datapsdatapsvclaitypsnmocleanpsall。请注意：在源码中可能有错误，执行之中有些内容无法实现。二叠加有以下命令：sustack、sudivstack，说明如下:SUSTACK - stack adjacent traces having the same key header word sustack stdout Optional parametersRequired parameters: noneOptional parameters: key=cdp header key word to stack on normpow=1.0 each sample is divided by the normpowth number of non-zero values stacked (normpow=0 selects no division) verbose=0 verbose = 1 echos informationSUDIVSTACK - Diversity Stacking using either average power or peak power within windows Required parameters:None Optional parameters: key=tracf key header word to stack on winlen=0.064 window length in seconds. typical choices: 0.064, 0.128, 0.256, 0.512, 1.024, 2.048, 4.096 if not specified the entire trace is used peak=1 peak power option default is average power examples： For duplicate field shot records: susort stack.data For CDP ordered data: sudivstack stack.data关于DEMOS，进入$CWPROOT/path/src/demos/Stacking?Phase_Weighted_Stacking，输入supws文件，观察主目录下的文件变更情况。三抽道集常用命令，susort，使用说明：SUSORT - sort on any segy header keywordssusort stdout +-key1 +-key2 . Susort supports any number of (secondary) keys with either ascending (+, the default) or descending (-) directions for each. The default sort key is cdp. Note: Only the following types of input/output are supported Disk input - any output Pipe input - Disk output Note: If the the CWP_TMPDIR environment variable is set use its value for the path; else use tmpfile() Example: To sort traces by cdp gather and within each gather by offset with both sorts in ascending order: susort OUTDATA cdp offset关于DEMOS，进入$CWPROOT/path/src/demos/sorting_trace/tutorial 输入Xsort，PSsort文件即可。由于以上三部分密不可分，特以速度分析一部分为例， 1、 产生数据 nx=100 fx=0 dx=80nz=110 fz=0 dz=50unif2 modeldata 2、 动校正 sunmo nmodata &3、 NMO图 suwind nmodata key=sx s=1400 j=250 |sugain tpow=2 gpow=0.5 |suximage label1=Time label2=Gather title=Every 5th Shot Record After NMO windowtitle=NMO legend=1 units=amplitude f2=1 d2=0.078125 f2num=1 d2num=5 wbox=$WIDTH hbox=$HEIGHT xbox=$WIDTHOFF3 ybox=$HEIGHTOFF &4、抽道集，叠加susort stackdata5、 抽道集，叠加绘图sugain stackdata tpow=2 gpow=0.5 |suximage label1=Time label2=CMP title=Stack windowtitle=Stack f2=350 d2=50 legend=1 units=amplitude wbox=$WIDTH hbox=$HEIGHT xbox=$WIDTHOFF4 ybox=$HEIGHTOFF &四反褶积常用命令sufxdecon，说明：SUFXDECON - random noise attenuation by FX-DECONvolution sufxdecon stdout Optional Parameters Required Parameters: None Optional Parameters: taper=.1 length of taper fmin=6. minimum frequency to process in Hz (accord to twlen) fmax=.6/(2*dt) maximum frequency to process in Hz twlen=entire trace time window length (minimum 300ms for lower freqs) ntrw=10 number of traces in window ntrf=4 number of traces for filter (smaller than ntrw) verbose=0 =1 for diagnostic print tmpdir= if non-empty, use the value as a directory path prefix for storing temporary files; else, if the CWP_TMPDIR environment variable is set, use its value for the path; else use tmpfile()关于DEMOS的使用，如下：进入$CWPROOT/src/demos/wienerlevinson,拷贝所有文件到主目录下，在中断中输入make,依次输入XDecon1 PSDecon1；XDecon2PSDecon2；XDecon3PSDecon3；XDecon4PSDecon4；XshapePSshape。例1、1、产生数据I=$CWPROOT/includeL=$CWPROOT/lib make ./traces |suaddhead ns=512 |sushw key=dt a=4000 |suaddnoise sn=200 |suconv filter=1152,-384,-904,288,174,-34,-12 modeldata rm traces2、显示模型suxwigb modeldata title=Data: 64ms reverbs windowtitle=Data label1=Time (sec) label2=Trace wbox=$WIDTH hbox=$HEIGHT xbox=$WIDTHOFF1 ybox=$HEIGHTOFF1 &3、自相关sugain modeldata tpow=$tpow |suacor nt=31 |suxwigb title=Autocorrelation windowtitle=AutoCorr label1=Time (sec) label2=Trace wbox=$WIDTH hbox=$HEIGHT xbox=$WIDTHOFF2 ybox=$HEIGHTOFF1 &4、脉冲化信号supef modeldata maxlag=.04 |suxwigb label1=Time label2=Trace title=Spiking Decon: 4,40ms windowtitle=Spike wbox=$WIDTH hbox=$HEIGHT xbox=$WIDTHOFF3 ybox=$HEIGHTOFF1 &5、脉冲化后的自相关supef modeldata maxlag=.04 |sugain tpow=$tpow |suacor nt=31 |suxwigb title=Autocor after spike windowtitle=AutoCorr/Spike label1=Time (sec) label2=Trace wbox=$WIDTH hbox=$HEIGHT xbox=$WIDTHOFF1 ybox=$HEIGHTOFF2 &6、追踪反射，反射注意反射层supef modeldata maxlag=.04 |supef minlag=.05 maxlag=.16 |suxwigb label1=Time label2=Trace title=PEF: 50,160ms (Note dipping reflector!) windowtitle=PEF wbox=$WIDTH hbox=$HEIGHT xbox=$WIDTHOFF2 ybox=$HEIGHTOFF2 &7、带通滤波supef susyn1.outsufilter susyn1.flt2、加噪声 suaddnoise susyn1.sn=4.flt3、绘图 suxmovie susyn1.flt label1=Time s label2=Trace # d1num=.1 d2num=10 n1=201 n2=100 d1=.004 d2=1 f1=0 f2=1 clip=30. -geometry $WIDTHx$HEIGHT+$WIDTHOFF1+$HEIGHTOFF1 title=Signal perc=99.5 &suxmovie susyn1.sn=4.flt label1=Time s label2=Trace # d1num=.1 d2num=10 n1=201 n2=100 d1=.004 d2=1 f1=0 f2=1 clip=30. -geometry $WIDTHx$HEIGHT+$WIDTHOFF1+$HEIGHTOFF2 title=Data perc=99.5 &5、反褶积sufxdecon susyn1.sn=4.fx ntrw=100suxmovie susyn1.sn=4.fx label1=Time s label2=Trace # d1num=.1 d2num=10 n1=201 n2=100 d1=.004 d2=1 f1=0 f2=1 clip=30. -geometry $WIDTHx$HEIGHT+$WIDTHOFF2+$HEIGHTOFF1 title=Fxdcn over Data perc=99.5 & sufxdecon susyn1.flt.fx ntrw=100suxmovie susyn1.flt.fx label1=Time s label2=Trace # d1num=.1 d2num=10 n1=201 n2=100 d1=.004 d2=1 f1=0 f2=1 clip=30. -geometry $WIDTHx$HEIGHT+$WIDTHOFF2+$HEIGHTOFF2 title=Fxdcn over Signal perc=99.5 &7、绘图suop2 susyn1.flt susyn1.flt.fx op=diff |suxmovie n1=201 n2=100 d1=.004 d2=1 f1=0 f2=1 clip=30. label1=Time s label2=Trace # d1num=.1 d2num=10 -geometry $WIDTHx$HEIGHT+$WIDTHOFF3+$HEIGHTOFF1 title=Diff: signal - fx(signal) perc=99.5 &suop2 susyn1.flt susyn1.sn=4.fx op=diff |suxmovie label1=Time s label2=Trace # d1num=.1 d2num=10 n1=201 n2=100 d1=.004 d2=1 f1=0 f2=1 clip=30. -geometry $WIDTHx$HEIGHT+$WIDTHOFF3+$HEIGHTOFF2 title=Diff: signal - fx(data) perc=99.5 &五带通滤波常用命令为subfilt，说明：SUBFILT - apply Butterworth bandpass filter subfilt stdout optional parameters Required parameters: if dt is not set in header, then dt is mandatory Optional parameters: (nyquist calculated internally) zerophase=1 =0 for minimum phase filter locut=1 =0 for no low cut filter hicut=1 =0 for no high cut filter fstoplo=0.10*(nyq) freq(Hz) in low cut stop band astoplo=0.05 upper bound on amp at fstoplo fpasslo=0.15*(nyq) freq(Hz) in low cut pass band apasslo=0.95 lower bound on amp at fpasslo fpasshi=0.40*(nyq) freq(Hz) in high cut pass band apasshi=0.95 lower bound on amp at fpasshi fstophi=0.55*(nyq) freq(Hz) in high cut stop band astophi=0.05 upper bound on amp at fstophi verbose=0 =1 for filter design info dt = (from header) time sampling interval (sec)简单的举个例子，Suplane data.suSubfilt dataoutSubxwigb data.suSubxwigb dataout关于DEMOS的使用，如下：进入目录，双击xbfilt1，xbfilt2.六偏移常用的命令为sugazmig、sumigps、sumigpspi。使用说明如下：SUGAZMIG - SU version of Jeno GAZDAGs phase-shift migration for zero-offset data. sugazmig outfile vfile= optional parameters Optional Parameters: dt=from header(dt) or .004 time sampling interval dx=from header(d2) or 1.0 midpoint sampling interval ft=0.0 first time sample ntau=nt(from data) number of migrated time samples dtau=dt(from header) migrated time sampling interval ftau=ft first migrated time sample tmig=0.0 times corresponding to interval velocities in vmig vmig=1500.0 interval velocities corresponding to times in tmig vfile= name of file containing velocities verbose=0 verbose = 1 echoes information tmpdir= if