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這期間可能有系統(tǒng)安裝所等待的時間,不用急,但凡遇到y(tǒng)es/no,一路y下來即可四.為了檢查是否安裝完畢,在終端中輸入Suplane>Suxwigb<&假設(shè)出現(xiàn)一個簡單的圖像,那么成功！CW收件的簡單說明一、文中涉及的命令全部以小寫形式,均可在終端窗口下輸入,以次獲取自述幫助.先說幾個命令：suplane和suxwigb,more.suplane作用是產(chǎn)生一個簡單的零偏移距su文件,suxwigb是一個典型的X-windows繪制圖形工具,如例子：suplane>suxwigb<more<比擬全面的了解它們,請在終端中輸入suplane,suxwigb,more.二、關(guān)于DEMOS勺應(yīng)用所有DEMOS、須把文件拷到用戶根目錄下,而后依,照readme文件中的執(zhí)行順序,在終端中輸入文件名.注意目錄下的文件變動.三、在執(zhí)行DEMOSC件時,如果想清楚了解程序執(zhí)行過程,請輸入moreprogramname由于水平有限,這里的謬誤很多,希望大家能在偶爾翻看時,多多留心,發(fā)現(xiàn)并改正,衷心希望能和大家一起學習.謝謝第一節(jié)兩種數(shù)據(jù)的轉(zhuǎn)換在CW應(yīng)用中免不了和兩種數(shù)據(jù)打交道,su和segy格式.它們有聯(lián)系也有區(qū)別..數(shù)據(jù)的輸入輸出1讀寫編輯數(shù)據(jù)常用命令如下：segyread---readanSEG-Ytapesegywite---writeanSEG-Ytapesegyclean—zeoroutunsignedportion〔局部〕ofheadersuaddhead---putheadersonbaretracesandsetthetrackandnsfieldssustrip—removetheSEGYheadersfromthetracessupaste---pasteexisting〔現(xiàn)存的〕SEGYheadersonexistingdatasegyhdrs—makesegyasciiandinnards〔內(nèi)部〕headersforsegywritebhedtopar,setbhed---editingthebinaryheaderfilesurange---getmaxandminvaluesfornon-zeroheaderentriessuchw---changeheaderwordusingoneortwoheaderwordfieldssugethw---gettheheaderwordsinsudatasuedit---examinesegydiskfilesandeditheaderssuxedit---examinesegydiskfilesandeditheaders2常規(guī)的數(shù)據(jù)轉(zhuǎn)換命令a2b—convertasciifloatstobinaryb2a---convertbinaryfloatstoasciiftnstrip---convertfortranfloatstoc-systelsfloatsrecast-recastdatatype〔convertfromonedatatypetoanother〕h2b—convert8bithexadecimalfloatstobinarytransp—transposeann1byn2elementmatrix.SEGYfSU的轉(zhuǎn)換SEGY文件格式由三局部組成,鏡象頭文件,十進制頭文件和實際的地震記錄,而SU格式的文件只含有其第三局部.命令segyread可以實現(xiàn)兩種格式的轉(zhuǎn)換,使用如下：在“big-endian〞平臺上,Segyreadtape=/dev/rmt0verbose=1endian=1>或者更常用一點的Segyreadtape=/dev/rmt0verbose=1endian=1|segyclean>命令行中斜線局部是SEGY磁帶的位置〔或者是磁盤文件的位置〕,其它應(yīng)用看自述文檔..典型命令的使用如下命令行輸入,注意主目錄中文件的變化和終端中的屏幕輸出.suplane>suxwigb<&segyhdrs<bhedtopar<binaryoutpar=setbhedbfile=binarypar=[parameter1=valueparamenter2=value]segywriteb2a<n1=5>a2b<n1=5>transp<n1=5>Su數(shù)據(jù)簡單操作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第二節(jié)CWP中圖形打印顯示工具一.X-WINDOWS的圖形打印.浮點數(shù)據(jù)繪圖常用的命令有

xcontour—xcontourplotoff(x1,x2)viavectorplotcallximage---ximgaeplotofauniformly-sampledfunctionf(x1,x2)xwigb---xwiggle-traceplotoff(x1,x2)viabitmapxgraph---xgraphergraphsn[i]pairsof(x,y)coordinatesxmovie---imageoneormoreframesofauniformlysampledfunctionf(x1,x2)相應(yīng)的,在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平臺上相應(yīng)的程序有：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,查看圖像.第三節(jié)常見處理以下幾乎每局部都有CW西帶的DEMOS參見$CWPROOT/path/src/demos目錄.對于一般的情況來說,DEMO的文件執(zhí)行具有一定的順序性.一.動校正這一局部有兩個命令,sunmo和suvelan.使用分別如下：SUNMO-NMOforanarbitraryvelocityfunctionoftimeandCDP命令格式為sunmo<stdin>stdout[optionalparameters]可選參數(shù)：tnmo=0相應(yīng)于NMG3速度的NMO]?間vnmo=2000相應(yīng)于NMG3時間的NM述度anis1=0產(chǎn)生四次節(jié)點的兩各向異性系數(shù)anis2=0時距曲線中,相應(yīng)于tanis2=0cdp=由tnmo和vnm彘義的CDPsmute=sampleswithNMOstretchexceedingsmutearezeroedlmute=25length(insamples)oflinearrampforstretchmutesscale=1=1todivideoutputsamplesbyNMOstretchfactorinvert=0=1toperform(approximate)inverseNMOixoffset=0donotconsidercross-lineoffset=1readcross-lineoffsetfromtraceheaderupward=0=1toscanupwardtofindfirstsampletokill注意：在處理常速度反射層時,vnmcM單一常數(shù)就可以了；在處理多層變速度反射層時,tnmo和vnmo總是成對出現(xiàn),且數(shù)據(jù)為相對應(yīng)的數(shù)組,如"tnmo=1,2vnmo=1500,2000…〞.關(guān)于NMO勺DEMO進入$CWPROOT/path/src/demos/NMO操作如下：.拷貝NMOF的所有文件到主目錄下依次輸入xmodel—makedata—xdata—xvclaity—xnmo—clean—xall;依次輸入psmodel—makedata—psdata—psvclaity—psnmo—clean一psallo請注意：在源碼中可能有錯誤,執(zhí)行之中有些內(nèi)容無法實現(xiàn).二.疊加有以下命令：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>關(guān)于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關(guān)于demos進入$CWPROOT/path/src/demos/sorting_trace/tutorial輸入Xsort,PSsort文件即可.由于以上三局部密不可分,特以速度分析一局部為例,1、產(chǎn)生數(shù)據(jù)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、產(chǎn)生數(shù)據(jù)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()關(guān)于DEMOS使用,如下：進入$CWPROOT/src/demos//wienerlevinson,拷貝所有文件到主目錄下,在中斷中輸入make,依次輸入XDecon1—PSDecon1；XDecon2—PSDecon2；XDecon3—PSDecon3；XDecon4—PSDecon4；Xshape一PSshape.例1、1、產(chǎn)生數(shù)據(jù)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、自相關(guān)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、脈沖化后的自相關(guān)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、數(shù)據(jù)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關(guān)于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()關(guān)于DEMOS使用,如下：進入目錄$CWPROOT/src/demos/Migration,可見(零偏移距數(shù)據(jù)的相移)：PSGazmig-Xgazmigstep:XdemoPademoshift-plus-interpolation:XdemoPsdemo依此順序,按手冊說明操作.疊后偏移舉以下例子例1、1、尖脈沖數(shù)據(jù)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、人工數(shù)據(jù)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、產(chǎn)生單炮信號單炮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