CWP软件操作手册版35p

CW啾件的安装在LINUX下建立用户CWP在CW吓建立目录path,将源文件放到path目录下,并建立bin文件夹在CWFffl户主目录下显示BI藏文件,修改.bash_profile文件,在已有的export之后另起一行,分别添加exportCWPROOT=/home/CWP/path再于PATH=$PATH:$HOME/bin添力口:/home/CWP/path/bin:/home/CWP退出保存从终端中分别输入cdpath…待终端中反映完毕,分别输入cdsrcmakeinstallmakextinstallmakemglinstallmakeutilsmakexminstallmakesfinstall这期间可能有系统安装所等待的时间,不用急,但凡遇到yes/no,一路y下来即可四.为了检查是否安装完毕,在终端中输入Suplane>Suxwigb<&假设出现一个简单的图像,那么成功！CW收件的简单说明一、文中涉及的命令全部以小写形式,均可在终端窗口下输入,以次获取自述帮助.先说几个命令：suplane和suxwigb,more.suplane作用是产生一个简单的零偏移距su文件,suxwigb是一个典型的X-windows绘制图形工具,如例子：suplane>suxwigb<more<比拟全面的了解它们,请在终端中输入suplane,suxwigb,more.二、关于DEMOS勺应用所有DEMOS、须把文件拷到用户根目录下,而后依,照readme文件中的执行顺序,在终端中输入文件名.注意目录下的文件变动.三、在执行DEMOSC件时,如果想清楚了解程序执行过程,请输入moreprogramname由于水平有限,这里的谬误很多,希望大家能在偶尔翻看时,多多留心,发现并改正,衷心希望能和大家一起学习.谢谢第一节两种数据的转换在CW应用中免不了和两种数据打交道,su和segy格式.它们有联系也有区别..数据的输入输出1读写编辑数据常用命令如下：segyread---readanSEG-Ytapesegywite---writeanSEG-Ytapesegyclean—zeoroutunsignedportion〔局部〕ofheadersuaddhead---putheadersonbaretracesandsetthetrackandnsfieldssustrip—removetheSEGYheadersfromthetracessupaste---pasteexisting〔现存的〕SEGYheadersonexistingdatasegyhdrs—makesegyasciiandinnards〔内部〕headersforsegywritebhedtopar,setbhed---editingthebinaryheaderfilesurange---getmaxandminvaluesfornon-zeroheaderentriessuchw---changeheaderwordusingoneortwoheaderwordfieldssugethw---gettheheaderwordsinsudatasuedit---examinesegydiskfilesandeditheaderssuxedit---examinesegydiskfilesandeditheaders2常规的数据转换命令a2b—convertasciifloatstobinaryb2a---convertbinaryfloatstoasciiftnstrip---convertfortranfloatstoc-systelsfloatsrecast-recastdatatype〔convertfromonedatatypetoanother〕h2b—convert8bithexadecimalfloatstobinarytransp—transposeann1byn2elementmatrix.SEGYfSU的转换SEGY文件格式由三局部组成,镜象头文件,十进制头文件和实际的地震记录,而SU格式的文件只含有其第三局部.命令segyread可以实现两种格式的转换,使用如下：在“big-endian〞平台上,Segyreadtape=/dev/rmt0verbose=1endian=1>或者更常用一点的Segyreadtape=/dev/rmt0verbose=1endian=1|segyclean>命令行中斜线局部是SEGY磁带的位置〔或者是磁盘文件的位置〕,其它应用看自述文档..典型命令的使用如下命令行输入,注意主目录中文件的变化和终端中的屏幕输出.suplane>suxwigb<&segyhdrs<bhedtopar<binaryoutpar=setbhedbfile=binarypar=[parameter1=valueparamenter2=value]segywriteb2a<n1=5>a2b<n1=5>transp<n1=5>Su数据简单操作suplane>sustrip<head=>suaddhead<ns=1024>supaste<ns=1024head=>surange<suplane|surangesugethw<key=tracl,,dt,ns|moresuplane|sugethwkey=tracl,tracr,offset,dtnsoutput=binary>moresushw<key=sxa=6400csuplane|sushwkey=sxa=6400c=-100j=5|sugethwkey=sx|moresuplane|sushwkey=offseta=200b=200j=5|sugethwkey=offset|moresuplane|sushwkey=dt,sx,offset,a=2000,6400,200b=0,0,200c=0,-100,0j=0,5,5|sugethwkey=dt,sx,offset|moresuedit第二节CWP中图形打印显示工具一.X-WINDOWS的图形打印.浮点数据绘图常用的命令有

xcontour—xcontourplotoff(x1,x2)viavectorplotcallximage---ximgaeplotofauniformly-sampledfunctionf(x1,x2)xwigb---xwiggle-traceplotoff(x1,x2)viabitmapxgraph---xgraphergraphsn[i]pairsof(x,y)coordinatesxmovie---imageoneormoreframesofauniformlysampledfunctionf(x1,x2)相应的,在SU平台上有以下程序suxcontoursuximagesuxwigbsuxgraphsuxmoviesuxmax例子：suplane|suxcontourtitle=suplane|suximgaetitle=suplane|suxwigbtitle=suplane|suxgraphtitle=suplane|suxmovietitle=suplane|suxmaxtilte=suplane|"contuour"&"image"&“w"contuour"&"image"&“wigb"&〞graph"&"movie"&max&>>suplane|suaddnoisesn=15>>suplane|suaddnoisesn=10>>suplane|suaddnoisesn=5>>suplane|suaddnoisesn=3>>suplane|suaddnoisesn=2>>suplane|suaddnoisesn=1>>suxmovie<n2=32title="frame=%gloop=1&.PS下的图形打印pscontourpsimagepswigbpsgraphpsmoviepscubepswigp在SU平台上相应的程序有：supscontoursupsimagesupswigbsupsgraphsupsmoviesupscube

supswigp例子:suplane>supscontour<title="contour〞>&supswigb<title=supsgraph<title=supsmovie<title=supscube<titlsupswigb<title=supsgraph<title=supsmovie<title=supscube<title=supswigp<title=bitmapwigbtraces〞>&graph〞>&movie〞>&coueplot〞>&wiggletrace〞>&运行以上程序之后,进入主目录,在所生成的*.eps文件上点击右键,翻开方式选择ggv,查看图像.第三节常见处理以下几乎每局部都有CW西带的DEMOS参见$CWPROOT/path/src/demos目录.对于一般的情况来说,DEMO的文件执行具有一定的顺序性.一.动校正这一局部有两个命令,sunmo和suvelan.使用分别如下：SUNMO-NMOforanarbitraryvelocityfunctionoftimeandCDP命令格式为sunmo<stdin>stdout[optionalparameters]可选参数：tnmo=0相应于NMG3速度的NMO]•间vnmo=2000相应于NMG3时间的NM述度anis1=0产生四次节点的两各向异性系数anis2=0时距曲线中,相应于tanis2=0cdp=由tnmo和vnm彘义的CDPsmute=sampleswithNMOstretchexceedingsmutearezeroedlmute=25length(insamples)oflinearrampforstretchmutesscale=1=1todivideoutputsamplesbyNMOstretchfactorinvert=0=1toperform(approximate)inverseNMOixoffset=0donotconsidercross-lineoffset=1readcross-lineoffsetfromtraceheaderupward=0=1toscanupwardtofindfirstsampletokill注意：在处理常速度反射层时,vnmcM单一常数就可以了；在处理多层变速度反射层时,tnmo和vnmo总是成对出现,且数据为相对应的数组,如"tnmo=1,2vnmo=1500,2000…〞.关于NMO勺DEMO进入$CWPROOT/path/src/demos/NMO操作如下：.拷贝NMOF的所有文件到主目录下依次输入xmodel—makedata—xdata—xvclaity—xnmo—clean—xall;依次输入psmodel—makedata—psdata—psvclaity—psnmo—clean一psallo请注意：在源码中可能有错误,执行之中有些内容无法实现.二.叠加有以下命令：sustack、sudivstack,说明如下：SUSTACKstackadjacenttraceshavingthesamekeyheaderwordsustack<stdin>stdout[Optionalparameters]Requiredparameters:noneOptionalparameters:key=cdpheaderkeywordtostackonnormpow=eachsampleisdividedbythenormpow'thnumberofnon-zerovaluesstacked(normpow=0selectsnodivision)verbose=0verbose=1echosinformationSUDIVSTACKDiversityStackingusingeitheraveragepowerorpeakpowerwithinwindowsRequiredparameters:NoneOptionalparameters:key=tracfkeyheaderwordtostackonwinlen=windowlengthinseconds.typicalchoices:,,,ifnotspecifiedtheentiretraceisusedpeak=1peakpoweroptiondefaultisaveragepowerexamples:Forduplicatefieldshotrecords:susort<tracf|sudivstack>ForCDPordereddata:sudivstack<key=cdp>关于DEMOS进入$CWPROOT/path/src/demos/Stacking?Phase_Weighted_Stacking,输入supws文件,观察主目录下的文件变更情况.三.抽道集常用命令,susort,使用说明：SUSORTsortonanysegyheaderkeywordssusort<stdin>stdout[[+-]key1[+-]key2…]susortsupportsanynumberof(secondary)keyswitheitherascending(递升)(+,thedefault)ordescending(-)directionsforeach.Thedefaultsortkeyiscdp.Note:Onlythefollowingtypesofinput/outputaresupportedDiskinput-->anyoutputPipeinput-->DiskoutputNote:IfthetheCWP_TMPDIRenvironmentvariableissetuseitsvalueforthepath;elseusetmpfile()Example:Tosorttracesbycdpgatherandwithineachgatherbyoffsetwithbothsortsinascendingorder:susort<INDATA>OUTDATAcdpoffset关于demos进入$CWPROOT/path/src/demos/sorting_trace/tutorial输入Xsort,PSsort文件即可.由于以上三局部密不可分,特以速度分析一局部为例,1、产生数据nx=100fx=0dx=80nz=110fz=0dz=50unif2<modelninf=4nx=$nxnz=$nzdx=$dxdz=$dz\v00=,,,,|ximagen1=$nzf1=$fzd1=$dzn2=$nxf2=$fxd2=$dxwindowtitle="Model"\legend=1units="wavespeed"\wbox=$WIDTHhbox=$HEIGHTxbox=$WIDTHOFFybox=$HEIGHTOFF&2、产生数据REF1=":,;,〞REF2=":,;,〞REF3=":,;,〞REF4=":,;,〞fpeak=#shouldbehalfnyquisttoavoidallspatialaliasingdvdz=#velocitygradient(dv/dz)dvdx=#velocitygradient(dv/dx)v00=#P-wavevelocityatsurfaceinkft/sob=0#toeliminatetheturnedraycontributionnxo=64fxo=dxo=#offsetinformation:number,first,spacingnxs=12fxs=dxs=#shotinformation:number,first,spacingnt=501dt=#timeinformation:number,spacingecho"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=>modeldata2、动校正sunmo<modeldatacdp=1500,2000,2500,3000tnmo=0,,,,vnmo=5000,6000,7050,8410,10000tnmo=0,,,,vnmo=5000,6000,7180,8350,9760tnmo=0,,,,vnmo=5000,6000,7180,8380,9590tnmo=0,,,vnmo=5000,7200,8440,9740>nmodata&3、NMOSsuwind<nmodatakey=sxs=1400j=250|sugaintpow=2gpow=|suximagelabel1="Time"label2="Gather"\title="Every5thShotRecordAfterNMO"\windowtitle="NMO"legend=1units="amplitude"\f2=1d2=0.078125f2num=1d2num=5\wbox=$WIDTHhbox=$HEIGHTxbox=$WIDTHOFF3ybox=$HEIGHTOFF&4、抽道集,叠加susort<nmodata|sustacknormpow=>stackdata5、抽道集,叠加绘图sugain<stackdatatpow=2gpow=|suximagelabel1="Time"label2="CMP"title="Stack"windowtitle="Stack"\f2=350d2=50legend=1units="amplitude"\wbox=$WIDTHhbox=$HEIGHTxbox=$WIDTHOFF4ybox=$HEIGHTOFF&四.反褶积常用命令sufxdecon,说明：SUFXDECONTandomnoiseattenuationbyFX-DECONvolutionsufxdecon<stdin>stdout[OptionalParameters]RequiredParameters:NoneOptionalParameters:taper=.1lengthoftaperfmin=6.minimumfrequencytoprocessinHz(accordtotwlen)fmax=.6/(2*dt)maximumfrequencytoprocessinHztwlen=entiretracetimewindowlength(minimum300msforlowerfreqs)ntrw=10numberoftracesinwindowntrf=4numberoftracesforfilter(smallerthanntrw)verbose=0=1fordiagnosticprinttmpdir=ifnon-empty,usethevalueasadirectorypathprefixforstoringtemporaryfiles;else,iftheCWP_TMPDIRenvironmentvariableisset,useitsvalueforthepath;elseusetmpfile()关于DEMOS使用,如下：进入$CWPROOT/src/demos//wienerlevinson,拷贝所有文件到主目录下,在中断中输入make,依次输入XDecon1—PSDecon1；XDecon2—PSDecon2；XDecon3—PSDecon3；XDecon4—PSDecon4；Xshape一PSshape.例1、1、产生数据I=${CWPROOT}/includeL=${CWPROOT}/libmake./traces|suaddheadns=512|sushwkey=dta=4000|suaddnoisesn=200|suconvfilter=1152,-384,-904,288,174,-34,-12>modeldatarmtraces2、显示模型suxwigb<modeldatatitle="Data:64msreverbs"\windowtitle="Data"\label1="Time(sec)"label2="Trace"\wbox=$WIDTHhbox=$HEIGHTxbox=$WIDTHOFF1ybox=$HEIGHTOFF1&3、自相关sugain<modeldatatpow=$tpow|suacornt=31|suxwigbtitle="Autocorrelation"\windowtitle="AutoCorr"\label1="Time(sec)"label2="Trace"\wbox=$WIDTHhbox=$HEIGHTxbox=$WIDTHOFF2ybox=$HEIGHTOFF1&4、脉冲化信号supef<modeldatamaxlag=.04|suxwigblabel1="Time"label2="Trace"title="SpikingDecon:4,40ms"\windowtitle="Spike"\wbox=$WIDTHhbox=$HEIGHTxbox=$WIDTHOFF3ybox=$HEIGHTOFF1&5、脉冲化后的自相关supef<modeldatamaxlag=.04|sugaintpow=$tpow|suacornt=31|suxwigbtitle="Autocorafterspike"\windowtitle="AutoCorr/Spike"\label1="Time(sec)"label2="Trace"\wbox=$WIDTHhbox=$HEIGHTxbox=$WIDTHOFF1ybox=$HEIGHTOFF2&6、追踪反射,反射注意反射层supef<modeldatamaxlag=.04|supefminlag=.05maxlag=.16|suxwigblabel1="Time"label2="Trace"\title="PEF:50,160ms(Notedippingreflector!)"\windowtitle="PEF"\wbox=$WIDTHhbox=$HEIGHTxbox=$WIDTHOFF2ybox=$HEIGHTOFF2&7、带通滤波supef<modeldatamaxlag=.04|supefminlag=.05maxlag=.16|sufilterf=5,15,80,100|suxwigblabel1="Time"label2="Trace"title="Spike,PEF,BandPass"\windowtitle="Final"\wbox=$WIDTHhbox=$HEIGHTxbox=$WIDTHOFF3ybox=$HEIGHTOFF2&例2、1、数据susynlvnt=201dt=0.004ft=nxo=1\nxs=100dxs=.025fxs=er=0ob=0\v00=dvdz=0dvdx=0fpeak=30smooth=1\ref="1:-1,1;7,5"\ref="1:-1,5;7,2"\ref="2:-1,5;7,9"|sushwkey=d2a=.05>sufilter<>2、加噪声suaddnoise<sn=8f=4,8,100,125amps=0,1,1,0|sufilter>=3、绘图suxmovie<\label1="Time[s]"label2="Trace#"d1num=.1d2num=10\n1=201n2=100d1=.004d2=1f1=0f2=1clip=30.\geometry${WIDTH}x${HEIGHT}+${WIDTHOFF1}+${HEIGHTOFF1}\title="Signal"perc=&suxmovie<=\label1="Time[s]"label2="Trace#"d1num=.1d2num=10\n1=201n2=100d1=.004d2=1f1=0f2=1clip=30.\geometry${WIDTH}x${HEIGHT}+${WIDTHOFF1}+${HEIGHTOFF2}\title="Data"perc=&5、反褶积sufxdecon<=>=ntrw=100suxmovie<=\label1="Time[s]"label2="Trace#"d1num=.1d2num=10\n1=201n2=100d1=.004d2=1f1=0f2=1clip=30.\geometry${WIDTH}x${HEIGHT}+${WIDTHOFF2}+${HEIGHTOFF1}\title="FxdcnoverData"perc=&label1="Time[s]"label2="Trace#"d1num=.1d2num=10\n1=201n2=100d1=.004d2=1f1=0f2=1clip=30.\geometry${WIDTH}x${HEIGHT}+${WIDTHOFF2}+${HEIGHTOFF2}\title="FxdcnoverSignal"perc=&7、绘图suxmovie\n1=201n2=100d1=.004d2=1f1=0f2=1clip=30.\label1="Time[s]"label2="Trace#"d1num=.1d2num=10\geometry${WIDTH}x${HEIGHT}+${WIDTHOFF3}+${HEIGHTOFF1}\title="Diff:signal-fx(signal)"perc=&suop2=op=diff|suxmovie\label1="Time[s]"label2="Trace#"d1num=.1d2num=10\n1=201n2=100d1=.004d2=1f1=0f2=1clip=30.\geometry${WIDTH}x${HEIGHT}+${WIDTHOFF3}+${HEIGHTOFF2}\title="Diff:signal-fx(data)"perc=&五.带通滤波常用命令为subfilt,说明：SUBFILT-applyButterworthbandpassfiltersubfilt<stdin>stdout[optionalparameters]Requiredparameters:ifdtisnotsetinheader,thendtismandatoryOptionalparameters:(nyquistcalculatedinternally)zerophase=1=0forminimumphasefilterlocut=1=0fornolowcutfilterhicut=1=0fornohighcutfilterfstoplo=*(nyq)freq(Hz)inlowcutstopbandastoplo=upperboundonampatfstoplofpasslo=*(nyq)freq(Hz)inlowcutpassbandapasslo=lowerboundonampatfpasslofpasshi=*(nyq)freq(Hz)inhighcutpassbandapasshi=lowerboundonampatfpasshifstophi=*(nyq)freq(Hz)inhighcutstopbandastophi=upperboundonampatfstophiverbose=0=1forfilterdesigninfodt=(fromheader)timesamplinginterval(sec)简单的举个例子,Suplane>Subfilt<>dataoutSubxwigb<Subxwigb<dataout关于DEMOS勺使用,如下：进入目录,双击xbfilt1,xbfilt2.六.偏移常用的命令为sugazmig、sumigps、sumigpspi.使用说明如下:SUGAZMIGSUversionofJenoGAZDAG'sphase-shiftmigrationforzero-offsetdata.sugazmig<infile>outfilevfile=[optionalparameters]OptionalParameters:dt=fromheader(dt)or.004timesamplingintervaldx=fromheader(d2)ormidpointsamplingintervalft=firsttimesamplentau=nt(fromdata)numberofmigratedtimesamplesdtau=dt(fromheader)migratedtimesamplingintervalftau=ftfirstmigratedtimesampletmig=timescorrespondingtointervalvelocitiesinvmigvmig=intervalvelocitiescorrespondingtotimesintmigvfile=nameoffilecontainingvelocitiesverbose=0verbose=1echoesinformationtmpdir=ifnon-empty,usethevalueasadirectorypathprefixforstoringtemporaryfiles;elseifthetheCWP_TMPDIRenvironmentvariableissetuseitsvalueforthepath;elseusetmpfile()SUMIGPSMIGrationbyPhaseShiftwithturningrayssumigps<stdin>stdout[optionalparms]RequiredParameters:NoneOptionalParameters:dt=fromheader(dt)or.004timesamplingintervaldx=fromheader(d2)ordistancebetweensucessivecdp'sffil=0,0,dt,dttrapezoidalwindowoffrequenciestomigratetmig=timescorrespondingtointervalvelocitiesinvmigvmig=intervalvelocitiescorrespondingtotimesintmigvfile=binary(non-ascii)filecontainingvelocitiesv(t)nxpad=0numberofcdpstopadwithzerosbeforeFFTltaper=0lengthoflineartaperforleftandrightedgesverbose=0=1fordiagnosticprinttmpdir=ifnon-empty,usethevalueasadirectorypathprefixforstoringtemporaryfiles;elseifthetheCWP_TMPDIRenvironmentvariableissetuseitsvalueforthepath;elseusetmpfile()SUMIGPSP-Gazdag'sphase-shiftplusinterpolationmigrationforzero-offsetdata,whichcanhandlethelateralvelocityvariation.sumigpspi<infile>outfilevfile=[optionalparameters]RequiredParameters:nz=numberofdepthsapmlesdz=depthsamplingintervalvfile=nameoffilecontainingvelocities(PleaseseeNotesbelowconcerningtheformatofvfile)OptionalParameters:dt=fromheader(dt)or.004timesamplingintervaldx=fromheader(d2)ormidpointsamplingintervaltmpdir=ifnon-empty,usethevalueasadirectorypathprefixforstoringtemporaryfiles;elseifthetheCWP_TMPDIRenvironmentvariableissetuseitsvalueforthepath;elseusetmpfile()关于DEMOS使用,如下：进入目录$CWPROOT/src/demos/Migration,可见(零偏移距数据的相移)：PSGazmig-Xgazmigstep:XdemoPademoshift-plus-interpolation:XdemoPsdemo依此顺序,按手册说明操作.叠后偏移举以下例子例1、1、尖脉冲数据WIDTH=420HEIGHT=300WIDTHOFF1=0WIDTHOFF2=430WIDTHOFF3=860HEIGHTOFF1=5HEIGHTOFF2=335HEIGHTOFF3=665perc=99n1=64n2=32suspikenspk=3it1=16ix1=16it2=32ix2=16it3=48ix3=16|sushwkey=dt,d2a=50000,.05>syndata1suxmovie<syndata1perc=$perc\title="Suspikedata"\label1="Time(sec)"label2="Midpoint(km)"\geometry${WIDTH}x${HEIGHT}+${WIDTHOFF1}+${HEIGHTOFF1}&2、以sugazmig作偏移sugazmig<syndata1tmig=0vmig=1|suxmovieperc=$perc\title="SUGAZMIGMigration"'label1="MigratedTime(sec)"label2="Midpoint(km)"\geometry${WIDTH}x${HEIGHT}+${WIDTHOFF2}+${HEIGHTOFF1}&3、以sumigps作偏移sumigps<syndata1tmig=0vmig=1|suxmovieperc=$perc\title="SUMIGPSMigration"'label1="MigratedTime(sec)"label2="Midpoint(km)"\geometry${WIDTH}x${HEIGHT}+${WIDTHOFF3}+${HEIGHTOFF1}&例2、1、人工数据susynlvnt=$n1dt=0.04ft=nxo=1\nxm=$n2dxm=.05fxm=er=0ob=1\v00=dvdz=0dvdx=0smooth=1\ref="0,.5;,.5;2.,;,;,;,.5;,.5"|sushwkey=d2a=.05>syndata2suxmovie<syndata2perc=$perc\title="SyntheticData"'label1="Time(sec)"label2="Distance(km)"\geometry${WIDTH}x${HEIGHT}+${WIDTHOFF1}+${HEIGHTOFF2}&2、以sugazmig作偏移sugazmig<syndata2tmig=0vmig=1|suxmovieperc=$perc\title="SUGAZMIGMigration"\label1="MigratedTime(sec)"label2="Midpoint(km)"\geometry${WIDTH}x${HEIGHT}+${WIDTHOFF2}+${HEIGHTOFF2}&3、以sumigps作偏移sumigps<syndata2tmig=0vmig=1|suxmovieperc=$perc\title="SUMIGPSMigration"\label1="MigratedTime(sec)"label2="Midpoint(km)"\geometry${WIDTH}x${HEIGHT}+${WIDTHOFF3}+${HEIGHTOFF2}&叠前偏移的例子1、产生单炮信号单炮2D叠前偏移suimp2dnt=512nshot=1nrec=200dgx=8dsx=8x0=640z0=768c=2000|sukfilterk=.001,.002,.3,.4amps=0,1,1,0>2、信号图suximage<title="1shotgatherdata"\xbox=10ybox=10wbox=350hbox=350&3、K域重建suspeck1k2<|suximagex1beg=x2beg=x2end=.8x1end=.8\xbox=400ybox=10wbox=350hbox=350\legend=1cmap=hsv4\title="K-domainrepresentation1shotgatherdata"&4、10炮记录suimp2dnt=512nshot=10nrec=200dgx=8dsx=160x0=640z0=768c=2000|sukfilterk=.001,.002,.3,.4amps=0,1,1,0>5、绘图suximage<title="10shotgathers"\xbox=10ybox=400wbox=350hbox=350perc=99&6、K域时窗suspeck1k2<|suximagex1beg=x2beg=x2end=x1end=\xbox=400ybox=400wbox=350hbox=350\legend=1cmap=hsv4\title="K-domainaperture10shotgathers,PSPI"&7、背景场速度makevelnx=200nz=1024v000=dx=8dz=dvdx=dvdz=>vfile1transp<vfile1n1=200>vfile28、分布傅立叶偏移sumigpresp<nxo=200nxshot=1nz=1024\dx=8dz=fmax=50vfile=vfile2>9、输出图像suximage<title="Prestackmigrationwithsumigpresp"\xbox=10ybox=10wbox=350hbox=350&10、K域时窗suspeck1k2<|suximagex1beg=x2beg=x2end=.2x1end=.2\xbox=350ybox=10wbox=350hbox=350\legend=1cmap=hsv4\title="K-domainaperture1shotgather,sumigpresp"&11、理论上KM时窗kaperturex0=640z0=768nshot=1sxmin=0szmin=0both=1\ngeo=200dgx=8dsx=1c=2000nstep=0nfreq=100fmin=20fmax=50|xgraphn=400000,0mark=1,0marksize=1,0linewidth=0,1style=normal\width=350height=350\title="TheoreticalK-domainaperture,1shotgatherdata"&12、背景速度makevelnx=200nz=1024v000=dx=8dz=dvdx=dvdz=>vfile1transp<vfile1n1=200>vfile213、有限查分偏移sumigprefd<nxo=200nxshot=1nz=1024\dx=8dz=fmax=50vfile=vfile2>14、图形输出suximage<title="Prestackmigrationwithsumigprfd"\xbox=10ybox=10wbox=350hbox=350&15、K域时窗suspeck1k2<|suximagex1beg=x2beg=x2end=.2x1end=.2\xbox=350ybox=10wbox=350hbox=350\legend=1cmap=hsv4\title="K-domainaperture1shotgather,sumigprefd"&16、理论K域时窗kaperturex0=640z0=768nshot=1sxmin=0szmin=0both=1\ngeo=200dgx=8dsx=1c=2000nstep=0nfreq=100fmin=20fmax=50|xgraphn=400000,0mark=1,0marksize=1,0linewidth=0,1style=normal\width=350height=350\title="TheoreticalK-domainaperture,1shotgatherdata"&17、背景速度makevelnx=200nz=1024v000=dx=8dz=dvdx=dvdz=>vfile1transp