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

1、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-待终端中反映完毕，分别

2、输入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磁带的位置（或者是磁盘文件的

3、位置），其它应用看自述文档。三典型命令的使用如下命令行输入，注意主目录中文件的变化和终端中的屏幕输出。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.bin

4、suaddhead 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,2

5、00 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

6、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=”gr

7、aph” &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.susupla

8、ne | 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目录。对于一般

9、的情况来说，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 sam

10、ples 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

11、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，操作如下：拷贝NMO下的所有文件到主目录下依次输入 xmodelmake dataxdataxvclaityxnmocleanxall；依次输入 psmodelmake datapsdatapsvclait

12、ypsnmocleanpsall。请注意：在源码中可能有错误，执行之中有些内容无法实现。二叠加有以下命令：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 n

13、ormpowth 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 w

14、indow 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，进入$CW

15、PROOT/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 d

16、efault 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

17、 offset with both sorts in ascending order: susort OUTDATA cdp offset关于DEMOS，进入$CWPROOT/path/src/demos/sorting_trace/tutorial 输入Xsort，PSsort文件即可。由于以上三部分密不可分，特以速度分析一部分为例， 产生数据 nx=100 fx=0 dx=80nz=110 fz=0 dz=50unif2 modeldata 动校正 sunmo nmodata &NMO图 suwind nmodata key=sx s=1400 j=250 |sugain tpow=2 g

18、pow=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 stackdata抽道集，叠加绘图sugain stackdata tpow=2 gpow=0.5 |suximage label1=Time

19、 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:

20、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) verb

21、ose=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,依次输入XD

22、econ1 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 rev

23、erbs 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

24、&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/

25、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=$WIDTHOF

26、F2 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.

27、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=.0

28、04 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+$HEIGHT

29、OFF2 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

30、.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 pa

31、rameters 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

32、 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

33、 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 migratio

34、n 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 ti

35、me 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 non-empty, use the value as

36、 a directory path prefix for storing temporary files; else if the the CWP_TMPDIR environment variable is set use its value for the path; else use tmpfile()SUMIGPS - MIGration by Phase Shift with turning rays sumigps stdout optional parms Required Parameters: None Optional Parameters: dt=from header(

37、dt) or .004 time sampling interval dx=from header(d2) or 1.0 distance between sucessive cdps ffil=0,0,0.5/dt,0.5/dt trapezoidal window of frequencies to migrate tmig=0.0 times corresponding to interval velocities in vmig vmig=1500.0 interval velocities corresponding to times in tmig vfile= binary (n

38、on-ascii) file containing velocities v(t) nxpad=0 number of cdps to pad with zeros before FFT ltaper=0 length of linear taper for left and right edges 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 the CWP_TMP

39、DIR environment variable is set use its value for the path; else use tmpfile()SUMIGPSPI - Gazdags phase-shift plus interpolation migration for zero-offset data, which can handle the lateral velocity variation. sumigpspi outfile vfile= optional parameters Required Parameters: nz= number of depth sapm

40、les dz= depth sampling interval vfile= name of file containing velocities (Please see Notes below concerning the format of vfile) Optional Parameters: dt=from header(dt) or .004 time sampling interval dx=from header(d2) or 1.0 midpoint sampling interval tmpdir= if non-empty, use the value as a direc

41、tory path prefix for storing temporary files; else if the the CWP_TMPDIR environment variable is set use its value for the path; else use tmpfile()关于DEMOS的使用，如下：进入目录$CWPROOT/src/demos/Migration，可见1.Gazdag(零偏移距数据的相移)：PSGazmigXgazmig2.Splitstep :XdemoPademo3.Phaseshift-plus-interpolation: XdemoPsdemo依

42、此顺序，按手册说明操作。叠后偏移举以下例子例1、尖脉冲数据 WIDTH=420HEIGHT=300WIDTHOFF1=0WIDTHOFF2=430WIDTHOFF3=860HEIGHTOFF1=5HEIGHTOFF2=335HEIGHTOFF3=665 perc=99n1=64n2=32suspike nspk=3 it1=16 ix1=16 it2=32 ix2=16 it3=48 ix3=16 |sushw key=dt,d2 a=50000,.05 syndata1suxmovie syndata1 perc=$perc title=Suspike data label1=Time (s

43、ec) label2=Midpoint (km) -geometry $WIDTHx$HEIGHT+$WIDTHOFF1+$HEIGHTOFF1 &2、以sugazmig作偏移sugazmig syndata1 tmig=0 vmig=1 |suxmovie perc=$perc title=SUGAZMIG Migration label1=Migrated Time (sec) label2=Midpoint (km) -geometry $WIDTHx$HEIGHT+$WIDTHOFF2+$HEIGHTOFF1 &3、以sumigps作偏移sumigps syndata2suxmovie

44、 syndata2 perc=$perc title=Synthetic Data label1=Time (sec) label2=Distance (km) -geometry $WIDTHx$HEIGHT+$WIDTHOFF1+$HEIGHTOFF2 &2、以sugazmig作偏移sugazmig syndata2 tmig=0 vmig=1 |suxmovie perc=$perc title=SUGAZMIG Migration label1=Migrated Time (sec) label2=Midpoint (km) -geometry $WIDTHx$HEIGHT+$WIDT

45、HOFF2+$HEIGHTOFF2 &3、以sumigps作偏移sumigps 1shot_spike.su2、信号图suximage 1shot_spike.su title=1 shotgather data xbox=10 ybox=10 wbox=350 hbox=350 &3、K域重建suspeck1k2 10shots_spike.su绘图suximage 10shots_spike.su title=10 shotgathers xbox=10 ybox=400 wbox=350 hbox=350 perc=99&K域时窗suspeck1k2 vfile1transp vfile

46、2分布傅立叶偏移sumigpresp prestack_mig_sp.su输出图像suximage prestack_mig_sp.su title=Prestack migration with sumigpresp xbox=10 ybox=10 wbox=350 hbox=350 &K域时窗suspeck1k2 vfile1transp vfile2有限查分偏移sumigprefd prestack_mig_fd.su图形输出suximage prestack_mig_fd.su title=Prestack migration with sumigprfd xbox=10 ybox=10 wbox=350 hbox=350 &15、K域时窗suspeck1k2 vfile1transp vfile2 18、傅立叶有限差分sumigpreffd prestack_mig_ffd.su19、偏移输出suximage prestack_mig_ffd.su title=Prestack migration with sumigrpeffd xbox=10 ybox=400 wbox=350 hbox=350 &20、K域