【正文】 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。第三節(jié) 常見處理以下幾乎每部分都有CWP自帶的DEMOS，參見$CWPROOT/path/src/demos目錄。2．PS 下的圖形打印pscontourpsimagepswigbpsgraphpsmoviepscubepswigp在SU平臺上相應的程序有：supscontoursupsimagesupswigbsupsgraphsupsmoviesupscubesupswigp例子：suplane supscontour title=”contour” amp。1．suplane suxwigb amp。二、關于DEMOS的應用 所有DEMOS必須把文件拷到用戶根目錄下，而后依照readme文件中的執(zhí)行順序，在終端中輸入文件名。若出現(xiàn)一個簡單的圖像，則成功！CWP軟件的簡單說明一、文中涉及的命令全部以小寫形式，均可在終端窗口下輸入，以次獲取自述幫助。一．數(shù)據(jù)的輸入輸出1 讀寫編輯數(shù)據(jù)常用命令如下：segyreadread an SEGY tapesegywitewrite an SEGY tapesegyclean—zeor out unsigned portion（部分） of headersuaddheadput headers on bare traces and set the track and ns fields sustrip—remove the SEGYheaders from the tracessupastepaste existing(現(xiàn)存的) SEGY headers on existing datasegyhdrs—make segyascii and innards(內(nèi)部) headers for segywrite bhedtopar,setbhedediting the binary header filesurangeget max and min values for nonzero header entriessuchwchange header word using one or two header word fieldssugethwget the header words in su datasueditexamine segy diskfiles and edit headerssuxeditexamine segy diskfiles and edit headers2 常規(guī)的數(shù)據(jù)轉(zhuǎn)換命令a2b—convert ascii floats to binaryb2aconvert binary floats to asciiftnstripconvert fortran floats to csystels floatsrecastrecast data type (convert from one data type to another)h2b—convert 8 bit hexadecimal floats to binarytransp—transpose an n1 by n2 element matrix二．SEGY與SU的轉(zhuǎn)換SEGY文件格式由三部分組成，鏡象頭文件，十進制頭文件和實際的地震記錄，而SU格式的文件只含有其第三部分。suplane | suxmovie title=”movie” amp。supscube title=”coue plot” amp。請注意：在源碼中可能有錯誤，執(zhí)行之中有些內(nèi)容無法實現(xiàn)。,REF3=:,。例產(chǎn)生數(shù)據(jù)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 traces顯示模型suxwigb modeldata title=Data: 64ms reverbs \ windowtitle=Data \ label1=Time (sec) label2=Trace \ wbox=$WIDTH hbox=$HEIGHT xbox=$WIDTHOFF1 ybox=$HEIGHTOFF1 amp。7.,.2 \ ref=2:1.,.5。使用說明如下：SUGAZMIG SU version of Jeno GAZDAG39。 else if the the CWP_TMPDIR environment variable is set use its value for the path。,。10炮記錄suimp2d nt=512 nshot=10 nrec=200 dgx=8 dsx=160 x0=640 z0=768 c=2000 |sukfilter k=.001,.002,.3,.4 amps=0,1,1,0 繪圖suximage title=10 shotgathers \xbox=10 ybox=400 wbox=350 hbox=350 perc=99amp。1背景速度makevel nx=200 nz=1024 v000= dx=8 dz= dvdx= dvdz= vfile1transp vfile1 n1=200vfile2 1傅立葉有限差分sumigpreffd nxo=200 nxshot=1 nz=1024 \dx=8 dz= fmax=50 vfile=vfile2 1偏移輸出suximage title=Prestack migration with sumigrpeffd \xbox=10 ybox=400 wbox=350 hbox=350 amp。suxwigb title=datanoise1 amp。sufilter f=10,20 amps=1,2 amp。suaddnoise sn=20 noise=gauss amp。sufilter | suxwigb amp。 it2 = 3*nt/4 time for spike 2 。 time for spike 4 。 trace for spike 3 。suxwigb amp。sunmo vnmo=2000suxwigb title= amp。sufilter f=10,20 amps=1,2 dt= amp。 suxwigb title= amp。2理論K域時窗 kaperture x0=640 z0=768 nshot=1 sxmin=0 szmin=0 both=1 \ngeo=200 dgx=8 dsx=1 c=2000 nstep=0 nfreq=100 fmin=20 fmax=50 |xgraph n=400000,0 mark=1,0 marksize=1,0 linewidth=0,1 style=normal \width=350 height=350 \title=Theoretical Kdomain aperture, 1 shotgather data amp。 背景場速度makevel nx=200 nz=1024 v000= dx=8 dz= dvdx= dvdz= vfile1transp vfile1 n1=200vfile2 分布傅立葉偏移sumigpresp nxo=200 nxshot=1 nz=1024 \dx=8 dz= fmax=50 vfile=vfile2 輸出圖像suximage title=Prestack migration with sumigpresp \xbox=10 ybox=10 wbox=350 hbox=350 amp。,.5。疊后偏移舉以下例子例 尖脈沖數(shù)據(jù) 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 (sec) label2=Midpoint (km) \ geometry ${WIDTH}x${HEIGHT}+${WIDTHOFF1}+${HEIGHTOFF1} amp。 else if the the CWP_TMPDIR environment variable is set use its value for the path。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}+${HEIGHTOFF2} \ title=Data perc= amp。脈沖化信號supef modeldata maxlag=.04 |suxwigb label1=Time label2=Trace title=Spiking Decon: 4,40ms \ windowtitle=Spike \ wbox=$WIDTH hbox=$HEIGHT xbox=$WIDTHOFF3 ybox=$HEIGHTOFF1 amp。, fpeak= should be half nyquist to avoid all spatial aliasingdvdz= velocity gradient (dv/dz)dvdx= velocity gradient (dv/dx)v00= Pwave velocity at surface in kft/so