【正文】 he 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。由于以上三部分密不可分，特以速度分析一部分為例， 產(chǎn)生數(shù)據(jù) nx=100 fx=0 dx=80nz=110 fz=0 dz=50unif2 model ninf=4 nx=$nx nz=$nz dx=$dx dz=$dz \ v00=, |ximage n1=$nz f1=$fz d1=$dz n2=$nx f2=$fx d2=$dx windowtitle=Model \ legend=1 units=wavespeed \ wbox=$WIDTH hbox=$HEIGHT xbox=$WIDTHOFF ybox=$HEIGHTOFF amp。,REF2=:,。,REF4=:,。 NMO圖 suwind nmodata key=sx s=1400 j=250 |sugain tpow=2 gpow= |suximage label1=Time label2=Gather \ title=Every 5th Shot Record After NMO \ windowtitle=NMO legend=1 units=amplitude \ f2=1 d2= f2num=1 d2num=5 \ wbox=$WIDTH hbox=$HEIGHT xbox=$WIDTHOFF3 ybox=$HEIGHTOFF amp。四．反褶積常用命令sufxdecon，說明：SUFXDECON random noise attenuation by FXDECONvolution sufxdecon stdin 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 nonempty, use the value as a directory path prefix for storing temporary files。 else use tmpfile()關(guān)于DEMOS的使用，如下：進入$CWPROOT/src/demos/……/wiener—levinson,拷貝所有文件到主目錄下，在中斷中輸入make,依次輸入XDecon1－ PSDecon1；XDecon2－PSDecon2；XDecon3－PSDecon3；XDecon4－PSDecon4；Xshape－PSshape。自相關(guān)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 amp。脈沖化后的自相關(guān)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 amp。帶通濾波supef modeldata maxlag=.04 |supef minlag=.05 maxlag=.16 |sufilter f=5,15,80,100 |suxwigb label1=Time label2=Trace title=Spike, PEF, BandPass \ windowtitle=Final \ wbox=$WIDTH hbox=$HEIGHT xbox=$WIDTHOFF3 ybox=$HEIGHTOFF2 amp。7.,.5 \ ref=1:1.,.5。7.,.9 |sushw key=d2 a=.05 sufilter 加噪聲 suaddnoise sn=8 f=4,8,100,125 amps=0,1,1,0 |sufilter =繪圖 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 ${WIDTH}x${HEIGHT}+${WIDTHOFF1}+${HEIGHTOFF1} \ title=Signal perc= amp。反褶積sufxdecon = = ntrw=100suxmovie = \ label1=Time [s] label2=Trace d1num=.1 d2num=10 \ n1=201 n2=100 d1=.004 d2=1 f1=0 f2=1 clip=30. \ geometry ${WIDTH}x${HEIGHT}+${WIDTHOFF2}+${HEIGHTOFF1} \ title=Fxd over Data perc= amp。繪圖suop2 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 ${WIDTH}x${HEIGHT}+${WIDTHOFF3}+${HEIGHTOFF1} \ title=Diff: signal fx(signal) perc= amp。五．帶通濾波常用命令為subfilt，說明：SUBFILT apply Butterworth bandpass filter subfilt stdin 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=*(nyq) freq(Hz) in low cut stop band astoplo= upper bound on amp at fstoplo fpasslo=*(nyq) freq(Hz) in low cut pass band apasslo= lower bound on amp at fpasslo fpasshi=*(nyq) freq(Hz) in high cut pass band apasshi= lower bound on amp at fpasshi fstophi=*(nyq) freq(Hz) in high cut stop band astophi= upper bound on amp at fstophi verbose=0 =1 for filter design info dt = (from header) time sampling interval (sec)簡單的舉個例子，Suplane Subfilt dataoutSubxwigb Subxwigb dataout關(guān)于DEMOS的使用，如下：進入目錄，雙擊xbfilt1，xbfilt2.六．偏移常用的命令為sugazmig、sumigps、sumigpspi。s phaseshift migration for zerooffset data. sugazmig infile outfile vfile= [optional parameters] Optional Parameters: dt=from header(dt) or .004 time sampling interval dx=from header(d2) or midpoint sampling interval ft= first time sample ntau=nt(from data) number of migrated time samples dtau=dt(from header) migrated time sampling interval ftau=ft f