CW 软件操作手册Word下载.docx
《CW 软件操作手册Word下载.docx》由会员分享，可在线阅读，更多相关《CW 软件操作手册Word下载.docx（24页珍藏版）》请在冰点文库上搜索。
suxwigb<
data.su
more<
比较全面的了解它们，请在终端中输入suplane,suxwigb,more。
二、关于DEMOS的应用所有DEMOS必须把文件拷到用户根目录下，而后依照readme文件中的执行顺序，在终端中输入文件名。
注意目录下的文件变动。
三、在执行DEMOS文件时，如果想清楚了解程序执行过程，请输入
moreprogramname
由于水平有限，这里的谬误很多，希望大家能在偶尔翻看时，多多留心，发现并改正，衷心希望能和大家一起学习。
谢谢
第一节两种数据的转换
在CWP应用中免不了和两种数据打交道，su和segy格式。
它们有联系也有区别。
一．数据的输入输出
1读写编辑数据
常用命令如下：
segyread---readanSEG-Ytape
segywite---writeanSEG-Ytape
segyclean—zeoroutunsignedportion（部分）ofheader
suaddhead---putheadersonbaretracesandsetthetrackandnsfields
sustrip—removetheSEGYheadersfromthetraces
supaste---pasteexisting（现存的）SEGYheadersonexistingdata
segyhdrs—makesegyasciiandinnards（内部）headersforsegywritebhedtopar,
setbhed---editingthebinaryheaderfile
surange---getmaxandminvaluesfornon-zeroheaderentries
suchw---changeheaderwordusingoneortwoheaderwordfields
sugethw---gettheheaderwordsinsudata
suedit---examinesegydiskfilesandeditheaders
suxedit---examinesegydiskfilesandeditheaders
2常规的数据转换命令
a2b—convertasciifloatstobinary
b2a---convertbinaryfloatstoascii
ftnstrip---convertfortranfloatstoc-systelsfloats
recast---recastdatatype（convertfromonedatatypetoanother）
h2b—convert8bithexadecimalfloatstobinary
transp—transposeann1byn2elementmatrix
二．SEGY与SU的转换
SEGY文件格式由三部分组成，镜象头文件，十进制头文件和实际的地震记录，而SU格式的文件只含有其第三部分。
命令segyread可以实现两种格式的转换，使用如下：
在“big-endian”平台上，
Segyreadtape=/dev/rmt0verbose=1endian=1>
data.su
或者更常用一点的
Segyreadtape=/dev/rmt0verbose=1endian=1|segyclean>
命令行中斜线部分是SEGY磁带的位置（或者是磁盘文件的位置），其它应用看自述文档。
三．典型命令的使用
如下命令行输入，注意主目录中文件的变化和终端中的屏幕输出。
1．
segyhdrs<
bhedtopar<
binaryoutpar=binary.par
setbhedbfile=binarypar=binary.par[parameter1=valueparamenter2=value……]
segywrite
b2a<
data.binaryn1=5>
data.ascii
a2b<
data.asciin1=5>
data.binary
transp<
data.trnsp
2．Su数据简单操作
sustrip<
data.suhead=data.headers>
data.bin
suaddhead<
data.binns=1024>
supaste<
data.binns=1024head=data.headers>
surange<
suplane|surange
sugethw<
data.sukey=tracl,tracr.offset,dt,ns|more
suplane|sugethwkey=tracl,tracr,offset,dtnsoutput=binary>
file.bin
morefile.bin
sushw<
data.sukey=sxa=6400c
suplane|sushwkey=sxa=6400c=-100j=5|sugethwkey=sx|more
suplane|sushwkey=offseta=200b=200j=5|sugethwkey=offset|more
suplane|sushwkey=dt,sx,offset,a=2000,6400,200b=0,0,200c=0,-100,0j=0,5,5|sugethwkey=dt,sx,offset|more
sueditdata.su
第二节CWP中图形打印显示工具
一．X-WINDOWS下的图形打印
1．浮点数据绘图
常用的命令有
xcontour—xcontourplotoff（x1,x2）viavectorplotcall
ximage---ximgaeplotofauniformly-sampledfunctionf（x1,x2）
xwigb---xwiggle-traceplotoff（x1,x2）viabitmap
xgraph---xgraphergraphsn[i]pairsof（x,y）coordinates
xmovie---imageoneormoreframesofauniformlysampledfunctionf（x1,x2）
相应的，在SU平台上有以下程序
suxcontour
suximage
suxwigb
suxgraph
suxmovie
suxmax
例子：
suplane|suxcontourtitle=”contuour”&
suplane|suximgaetitle=”image”&
suplane|suxwigbtitle=”wigb”&
suplane|suxgraphtitle=”graph”&
suplane|suxmovietitle=”movie”&
suplane|suxmaxtilte=”max”&
suplane|junk1.su
suplane|suaddnoisesn=20>
>
junk1.su
suplane|suaddnoisesn=15>
suplane|suaddnoisesn=10>
suplane|suaddnoisesn=5>
suplane|suaddnoisesn=3>
suplane|suaddnoisesn=2>
suplane|suaddnoisesn=1>
suxmovie<
junk1.sun2=32title=”frame=%g”loop=1&
2．PS下的图形打印
pscontour
psimage
pswigb
psgraph
psmovie
pscube
pswigp
在SU平台上相应的程序有：
supscontour
supsimage
supswigb
supsgraph
supsmovie
supscube
supswigp
junk.su
supscontour<
junk.sutitle=”contour”>
data1.eps&
supsimage<
junk.sutitle=”image”label1=“sec”label2=”tracenumber”>
data2.eps&
supswigb<
junk.sutitle=”bitmapwigbtraces”>
data3.eps&
supsgraph<
junk.sutitle=”graph”>
data4.eps&
supsmovie<
junk.sutitle=”movie”>
data5.eps&
supscube<
junk.sutitle=”coueplot”>
data6.eps&
supswigp<
junk.sutitle=”wiggletrace”>
data7.eps&
运行以上程序之后，进入主目录，在所生成的*.eps文件上点击右键，打开方式选择ggv，查看图像。
第三节常见处理
以下几乎每部分都有CWP自带的DEMOS，参见$CWPROOT/path/src/demos目录。
对于一般的情况来说，DEMOS的文件执行具有一定的顺序性。
一．动校正
这一部分有两个命令，sunmo和suvelan。
使用分别如下：
SUNMO-NMOforanarbitraryvelocityfunctionoftimeandCDP
命令格式为
sunmo<
stdin>
stdout[optionalparameters]
可选参数：
tnmo=0相应于NMO中速度的NMO时间
vnmo=2000相应于NMO中时间的NMO速度
anis1=0产生四次节点的两各向异性系数
anis2=0时距曲线中，相应于tnmo的时间
cdp=由tnmo和vnmo定义的CDP
smute=1.5sampleswithNMOstretchexceedingsmutearezeroed
lmute=25length（insamples）oflinearrampforstretchmute
sscale=1=1todivideoutputsamplesbyNMOstretchfactor
invert=0=1toperform（approximate）inverseNMO
ixoffset=0donotconsidercross-lineoffset
=1readcross-lineoffsetfromtraceheader
upward=0=1toscanupwardtofindfirstsampletokill
注意:
在处理常速度反射层时，vnmo取单一常数就可以了；
在处理多层变速度反射层时，tnmo和vnmo总是成对出现,且数据为相对应的数组,如”tnmo=1,2…vnmo=1500,2000…”。
关于NMO的DEMO，进入$CWPROOT/path/src/demos/NMO，操作如下：
1．拷贝NMO下的所有文件到主目录下
2．依次输入xmodel－makedata－xdata－xvclaity－xnmo－clean－xall；
3．依次输入psmodel－makedata－psdata－psvclaity－psnmo－clean－psall。
请注意：
在源码中可能有错误，执行之中有些内容无法实现。
二．叠加
有以下命令：
sustack、sudivstack，说明如下:
SUSTACK-stackadjacenttraceshavingthesamekeyheaderword
sustack<
stdout[Optionalparameters]
Requiredparameters:
none
Optionalparameters:
key=cdpheaderkeywordtostackon
normpow=1.0eachsampleisdividedbythe
normpow'
thnumberofnon-zerovalues
stacked（normpow=0selectsnodivision）
verbose=0verbose=1echosinformation
SUDIVSTACK-DiversityStackingusingeitheraveragepowerorpeak
powerwithinwindows
Requiredparameters:
None
Optionalparameters:
key=tracfkeyheaderwordtostackon
winlen=0.064windowlengthinseconds.
typicalchoices:
0.064,0.128,0.256,
0.512,1.024,2.048,4.096
ifnotspecifiedtheentiretraceisused
peak=1peakpoweroptiondefaultisaveragepower
examples：
Forduplicatefieldshotrecords:
susort<
field.datatracf|sudivstack>
stack.data
ForCDPordereddata:
sudivstack<
cdp.datakey=cdp>
关于DEMOS，进入$CWPROOT/path/src/demos/Stacking?
Phase_Weighted_Stacking，输入supws文件，观察主目录下的文件变更情况。
三．抽道集
常用命令，susort，使用说明：
SUSORT-sortonanysegyheaderkeywords
susort<
stdout[[+-]key1[+-]key2...]
susortsupportsanynumberof（secondary）keyswitheitherascending（递升）（+,thedefault）ordescending（-）directionsforeach.Thedefaultsortkeyiscdp.
Note:
Onlythefollowingtypesofinput/outputaresupported
Diskinput-->
anyoutput
Pipeinput-->
Diskoutput
IfthetheCWP_TMPDIRenvironmentvariableissetuse
itsvalueforthepath;
elseusetmpfile（）
Example:
Tosorttracesbycdpgatherandwithineachgatherbyoffsetwithbothsortsinascendingorder:
INDATA>
OUTDATAcdpoffset
关于DEMOS，进入$CWPROOT/path/src/demos/sorting_trace/tutorial输入Xsort，PSsort文件即可。
由于以上三部分密不可分，特以速度分析一部分为例，
1、产生数据
nx=100fx=0dx=80
nz=110fz=0dz=50
unif2<
modelninf=4nx=$nxnz=$nzdx=$dxdz=$dz\
v00=5000.0,6000.0,8000.0,10000.0,15000.0|
ximagen1=$nzf1=$fzd1=$dzn2=$nxf2=$fxd2=$dxwindowtitle="
Model"
\
legend=1units="
wavespeed"
wbox=$WIDTHhbox=$HEIGHTxbox=$WIDTHOFFybox=$HEIGHTOFF&
2、产生数据
REF1="
0.0909091:
-4.0,1.0;
12.0,1.0"
REF2="
0.1428570:
-4.0,2.2;
12.0,2.2"
REF3="
0.1111110:
-4.0,3.5;
12.0,3.5"
REF4="
0.2000000:
-4.0,5.0;
12.0,5.0"
fpeak=25.0#shouldbehalfnyquisttoavoidallspatialaliasing
dvdz=2.0#velocitygradient（dv/dz）
dvdx=0.0#velocitygradient（dv/dx）
v00=5.0#P-wavevelocityatsurfaceinkft/s
ob=0#toeliminatetheturnedraycontribution
nxo=64fxo=0.1dxo=0.1#offsetinformation:
number,first,spacing
nxs=12fxs=1.4dxs=-0.1#shotinformation:
nt=501dt=0.004#timeinformation:
number,spacing
echo"
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=25.0>
modeldata
2、动校正
sunmo<
modeldatacdp=1500,2000,2500,3000
tnmo=0,0.34,0.625,0.875,1.13vnmo=5000,6000,7050,8410,10000tnmo=0,0.335,0.630,0.875,1.1vnmo=5000,6000,7180,8350,9760tnmo=0,0.35,0.64,0.88,1.1vnmo=5000,6000,7180,8380,9590
tnmo=0,0.64,0.88,1.1vnmo=5000,7200,8440,9740>
nmodata&
3、NMO图
suwind<
nmodatakey=sxs=1400j=250|
sugaintpow=2gpow=0.5|
suximagelabel1="
Time"
label2="
Gather"
title="
Every5thShotRecordAfterNMO"
windowtitle="
NMO"
amplitude"
f2=1d2=0.078125f2num=1d2num=5\
wbox=$WIDTHhbox=$HEIGHTxbox=$WIDTHOFF3ybox=$HEIGHTOFF&
4、抽道集，叠加
nmodata|
sustacknormpow=1.0>
stackdata
5、抽道集，叠加绘图
sugain<
stackdatatpow=2gpow=0.5|
CMP"
Stack"
f2=350d2=50legend=1units="
wbox=$WIDTHhbox=$HEIGHTxbox=$WIDTHOFF4ybox=$HEIGHTOFF&
四．反褶积
常用命令sufxdecon，说明：
SUFX
升级会员 