当前位置:首页 >> 能源/化工 >>

su简单中文使用手册


CWP 软件的安装与简单使用手册

CWP 软件的安装

一.在 LINUX 下建立用户 CWP,在 CWP 下建立目录 path,将源文件 cwp.su.all.37.tar.Z 放 二.到 path 目录下,并建立 bin 文件夹 三.在 CWP 用户主目录下显示隐藏文件,修改.bash_profile 文件,在已有的 export

之后另 起一行,分别添加 export CWPROOT=/home/CWP/path,再于 PATH=$PATH:$HOME/bin 后添加 :/home/CWP/path/bin:/home/CWP 退出保存 四.从终端中分别输入 cd path zcat cwp.su.all.37.tar.Z | tar –xvf… 待终端中反映完毕,分别输入 cd src make install make xtinstall make mglinstall make utils make xminstall make sfinstall 这期间可能有系统安装所等待的时间,不用急,但凡遇到 yes/no,一路 y 下来即可。 四.为了检查是否安装完毕,在终端中输入 Suplane > data.su Suxwigb < data.su & 若出现一个简单的图像,则成功!

CWP 软件的简单说明

一、文中涉及的命令全部以小写形式,均可在终端窗口下输入,以次获取自述帮助。 先说几个命令:suplane 和 suxwigb,more。 suplane 作用是产生一个简单的零偏移距 su 文件,suxwigb 是一个典型的 X—windows 绘制图形工具,如例子: suplane > data.su suxwigb < data.su more < data.su 比较全面的了解它们,请在终端中输入 suplane , suxwigb ,more 。 二、关于 DEMOS 的应用 所有 DEMOS 必须把文件拷到用户根目录下,而后依照 readme 文件中的执行顺序,在终端中输入文件名。注意目录下的文件变动。 三、在执行 DEMOS 文件时,如果想清楚了解程序执行过程,请输入 more programname

由于水平有限,这里的谬误很多,希望大家能在偶尔翻看时,多多留心,发现并改正,衷心希望能和 大家一起学习。谢谢

第一节 两种数据的转换 在 CWP 应用中免不了和两种数据打交道,su 和 segy 格式。它们有联系也有区别。 一.数据的输入输出 1 读写编辑数据 常用命令如下: segyread---read an SEG-Y tape segywite---write an SEG-Y tape segyclean—zeor out unsigned portion of header suaddhead---put headers on bare traces and set the track and ns fields sustrip—remove the SEGYheaders from the traces supaste---paste exitsting SEGY headers on existing data segyhdrs—make seg-yascii and inary headers for segywrite bhedtopar,setbhed---editing the binary header file surange---get max and min values for non-zero header entries suchw---change header word using one or two header word fields sugethw---get the header words in su data suedit---examine segy diskfiles and edit headers suxedit---examine segy diskfiles and edit headers 2 常规的数据转换命令 a2b—convert ascii floats to binary b2a---convert binary floats to ascii ftnstrip---convert fortran floats to c-systels floats recast---recast data type (convert from one data type to another) h2b—convert 8 bit hexadecimal floats to binary transp—transpose an n1 by n2 element matrix 二.SEGY 与 SU 的转换 SEGY 文件格式由三部分组成,镜象头文件,十进制头文件和实际的地震记录,而 SU 格式的文件只含有其第三部分。命令 segyread 可以实现两种格式的转换,使用如下: 在“big-endian”平台上, Segyread tape=/dev/rmt0 verbose=1 endian=1 >data.su 或者更常用一点的 Segyread tape=/dev/rmt0 verbose=1 endian=1 | segyclean > data.su 命令行中斜线部分是 SEGY 磁带的位置(或者是磁盘文件的位置) ,其它应用看自述文档。 三.典型命令的使用 如下命令行输入,注意主目录中文件的变化和终端中的屏幕输出。 1. suplane > data.su suxwigb < data.su & segyhdrs < data.su bhedtopar < binary outpar=binary.par setbhed bfile=binary par=binary.par [parameter1=value paramenter2=value……] segywrite b2a < data.binary n1=5 > data.ascii a2b <data.ascii n1=5 > data.binary

transp < data.binary n1=5 >data.trnsp 2.Su 数据简单操作 suplane >data.su sustrip < data.su head=data.headers > data.bin suaddhead <data.bin ns=1024 > data.su supaste < data.bin ns=1024 head=data.headers >data.su surange < data.su suplane | surange sugethw < data.su key=tracl,tracr.offset,dt,ns | more suplane |sugethw key=tracl,tracr,offset,dt ns output=binary > file.bin more file.bin sushw < data.su key=sx a=6400 c=-100 j=5 > data.new.su suplane |sushw key=sx a=6400 c=-100 j=5 | sugethw key=sx | more suplane | sushw key=offset a=200 b=200 j=5 | sugethw key=offset | more suplane | sushw key=dt,sx ,offset,a=2000,6400,200 b=0,0,200 c=0,-100,0j=0,5,5 | sugethw key=dt,sx,offset | more suedit data.su

第二节 CWP 中图形打印显示工具 一.X-WINDOWS 下的图形打印 1.浮点数据绘图 常用的命令有 xcontour—x contour plot of f(x1,x2) via vector plot call ximage---x imgae plot of a uniformly-sampled function f(x1,x2) xwigb---x wiggle-trace plot of f(x1,x2) via bitmap xgraph---x grapher graphs n[i] pairs of (x,y) coordinates xmovie ---image one or more frames of a uniformly sampled function f(x1,x2) 相应的,在 SU 平台上有以下程序 suxcontour suximage suxwigb suxgraph suxmovie suxmax 例子: suplane | suxcontour title=‖contuour‖ & suplane | suximgae title=‖image‖ & suplane | suxwigb title=‖wigb‖ & suplane | suxgraph title=‖graph‖ & suplane | suxmovie title=‖movie‖ & suplane | suxmax tilte=‖max‖ & suplane | junk1.su suplane | suaddnoise sn=20 >>junk1.su suplane | suaddnoise sn=15 >>junk1.su suplane | suaddnoise sn=10 >>junk1.su suplane | suaddnoise sn=5 >>junk1.su suplane | suaddnoise sn=3 >>junk1.su suplane | suaddnoise sn=2 >>junk1.su suplane | suaddnoise sn=1 >>junk1.su suxmovie < junk1.su n2=32 title=‖frame=%g‖ loop=1 & 2.PS 下的图形打印 pscontour psimage pswigb psgraph psmovie pscube pswigp 在 SU 平台上相应的程序有: supscontour

supsimage supswigb supsgraph supsmovie supscube supswigp 例子: suplane > junk.su supscontour <junk.su title=‖contour‖ > data1.eps & supsimage <junk.su title=‖image‖ label1=―sec‖ label2=‖trace number‖ > data2.eps & supswigb <junk.su title=‖bitmap wigb traces‖ > data3.eps & supsgraph <junk.su title=‖graph‖ > data4.eps & supsmovie <junk.su title=‖movie‖ > data5.eps & supscube <junk.su title=‖coue plot‖ > data6.eps & supswigp <junk.su title=‖wiggle trace‖ > data7.eps & 运行以上程序之后,进入主目录,在所生成的*.eps 文件上点击右键,打开方式选择 ggv, 查看图像。

第三节 常见处理 以下几乎每部分都有 CWP 自带的 DEMOS,参见$CWPROOT/path/src/demos 目录。对 于一般的情况来说,DEMOS 的文件执行具有一定的顺序性。 一.动校正 这一部分有两个命令,sunmo 和 suvelan。使用分别如下:
SUNMO - NMO for an arbitrary velocity function of time and CDP 命令格式为 sunmo <stdin >stdout [optional parameters] 可选参数: tnmo=0 vnmo=2000 anis1=0 anis2=0 cdp= smute=1.5 lmute=25 sscale=1 invert=0 ixoffset=0 upward=0 相应于NMO中速度的NMO时间 相应于NMO中时间的NMO速度 产生四次节点的两各向异性系数 时距曲线中,相应于tnmo的时间 由tnmo和vnmo定义的CDP samples with NMO stretch exceeding smute are zeroed length (in samples) of linear ramp for stretch mute =1 to divide output samples by NMO stretch factor =1 to perform (approximate) inverse NMO do not consider cross-line offset =1 read cross-line offset from trace header =1 to scan upward to find first sample to kill 注意:在处理常速度反射层时, vnmo取单一常数就可以了; 在处理多层变速度反射层时, tnmo 和vnmo总是成对出现,且数据为相对应的数组,如‖tnmo=1,2… vnmo=1500,2000…‖。

关于 NMO 的 DEMO,进入$CWPROOT/path/src/demos/NMO,操作如下: 1. 拷贝 NMO 下的所有文件到主目录下 2. 依次输入 xmodel-make data-xdata-xvclaity-xnmo-clean-xall; 3. 依次输入 psmodel-make data-psdata-psvclaity-psnmo-clean-psall。 请注意:在源码中可能有错误,执行之中有些内容无法实现。 二.叠加 有以下命令:sustack、sudivstack,说明如下:
SUSTACK - stack adjacent traces having the same key header word sustack <stdin >stdout [Optional parameters] Required parameters: none Optional parameters: key=cdp normpow=1.0 header key word to stack on each sample is divided by the normpow'th number of non-zero values stacked (normpow=0 selects no division) verbose=0 verbose = 1 echos information

SUDIVSTACK -

Diversity Stacking using either average power or peak power within windows

Required parameters: None Optional parameters: key=tracf winlen=0.064 key header word to stack on window length in seconds. typical choices: 0.064, 0.128, 0.256, 0.512, 1.024, 2.048, 4.096 if not specified the entire trace is used peak=1 examples: For duplicate field shot records: susort < field.data tracf | sudivstack > stack.data For CDP ordered data: sudivstack < cdp.data key=cdp > stack.data 关于DEMOS,进入$CWPROOT/path/src/demos/Stacking?Phase_Weighted_Stacking,输入 peak power option default is average power

supws文件,观察主目录下的文件变更情况。 三.抽道集 常用命令,susort,使用说明:
SUSORT - sort on any segy header keywords susort <stdin >stdout [[+-]key1 [+-]key2 ...] Susort supports any number of (secondary) keys with either ascending (+, the default) or descending (-) directions for each. Note: Disk input --> any output Pipe input --> Disk output Note: If the the CWP_TMPDIR environment variable is set use its value for the path; else use tmpfile() Example: To sort traces by cdp gather and within each gather by offset with both sorts in ascending order: susort <INDATA >OUTDATA cdp offset 关于DEMOS,进入$CWPROOT/path/src/demos/sorting_trace/tutorial 输入Xsort,PSsort文件 The default sort key is cdp. Only the following types of input/output are supported

即可。 由于以上三部分密不可分,特以速度分析一部分为例, 1、 产生数据

nx=100 fx=0 dx=80 nz=110 fz=0 dz=50 unif2 <model ninf=4 nx=$nx nz=$nz dx=$dx dz=$dz \ v00=5000.0,6000.0,8000.0,10000.0,15000.0 | 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 &

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 dvdz=2.0 dvdx=0.0 v00=5.0 ob=0 # should be half nyquist to avoid all spatial aliasing # velocity gradient (dv/dz) # velocity gradient (dv/dx) # P-wave velocity at surface in kft/s # to eliminate the turned ray contribution # offset information: number, first, spacing # shot information: number, first, spacing # time information: number, spacing

nxo=64 fxo=0.1 dxo=0.1 nxs=12 fxs=1.4 dxs=-0.1 nt=501 dt=0.004

echo "Making synthetic data for demo" susynlv v00=$v00 fpeak=$fpeak ndpfz=10 dvdz=$dvdz ob=$ob \ dt=$dt nt=$nt verbose=0 \ ref=$REF1 ref=$REF2 ref=$REF3 ref=$REF4 \ nxo=$nxo fxo=$fxo dxo=$dxo \ nxs=$nxs fxs=$fxs dxs=$dxs |

suchw key1=cdp key2=gx key3=sx b=1 c=1 d=2 | suaddnoise sn=50 f1=4.0 f2=8.0 f3=20.0 f4=25.0 >modeldata

2、 动校正
sunmo < modeldata cdp=1500,2000,2500,3000 tnmo=0,0.34,0.625,0.875,1.13 vnmo=5000,6000,7050,8410,10000 tnmo=0,0.335,0.630,0.875,1.1 vnmo=5000,6000,7180,8350,9760 tnmo=0,0.35,0.64,0.88,1.1 vnmo=5000,6000,7180,8380,9590 tnmo=0,0.64,0.88,1.1 vnmo=5000,7200,8440,9740 > nmodata & 3、 NMO图 suwind <nmodata key=sx s=1400 j=250 | sugain tpow=2 gpow=0.5 | suximage label1="Time" label2="Gather" \ title="Every 5th Shot Record After NMO" \ windowtitle="NMO" legend=1 units="amplitude" \ f2=1 d2=0.078125 f2num=1 d2num=5 \ wbox=$WIDTH hbox=$HEIGHT xbox=$WIDTHOFF3 ybox=$HEIGHTOFF &

4、抽道集,叠加 susort <nmodata | sustack normpow=1.0 >stackdata 5、 抽道集,叠加绘图 sugain <stackdata tpow=2 gpow=0.5 | suximage label1="Time" label2="CMP" title="Stack" windowtitle="Stack" \ f2=350 d2=50 legend=1 units="amplitude" \ wbox=$WIDTH hbox=$HEIGHT xbox=$WIDTHOFF4 ybox=$HEIGHTOFF &

四.反褶积 常用命令 sufxdecon,说明:
SUFXDECON - random noise attenuation by FX-DECONvolution sufxdecon <stdin >stdout [Optional Parameters] Required Parameters: None Optional Parameters: taper=.1 fmin=6. fmax=.6/(2*dt) ntrw=10 ntrf=4 length of taper minimum frequency to process in Hz maximum frequency to process in Hz time window length (minimum 300ms for lower freqs) number of traces in window number of traces for filter (smaller than ntrw) (accord to twlen)

twlen=entire trace

verbose=0 tmpdir=

=1 for diagnostic print if non-empty, use the value as a directory path prefix for storing temporary files; else, if the CWP_TMPDIR environment variable is set, use its value for the path; else use tmpfile()

关于 DEMOS 的使用,如下: 进入 $CWPROOT/src/demos/……/wiener—levinson,拷贝所有文件到主目录下,在中断中输入
make,依次输入 XDecon1- PSDecon1;XDecon2-PSDecon2;XDecon3-PSDecon3;XDecon4- PSDecon4;Xshape-PSshape。 例 1、 1、产生数据 I=${CWPROOT}/include L=${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 2、显示模型 suxwigb <modeldata title="Data: 64ms reverbs" \ windowtitle="Data" \ label1="Time (sec)" label2="Trace" \ wbox=$WIDTH hbox=$HEIGHT xbox=$WIDTHOFF1 ybox=$HEIGHTOFF1 &

3、自相关 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 &

4、脉冲化信号 supef <modeldata maxlag=.04 | suxwigb label1="Time" windowtitle="Spike" \ wbox=$WIDTH hbox=$HEIGHT xbox=$WIDTHOFF3 ybox=$HEIGHTOFF1 & label2="Trace" title="Spiking Decon: 4,40ms" \

5、脉冲化后的自相关 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 &

6、追踪反射,反射注意反射层 supef <modeldata maxlag=.04 | supef minlag=.05 maxlag=.16 | suxwigb label1="Time" windowtitle="PEF" \ wbox=$WIDTH hbox=$HEIGHT xbox=$WIDTHOFF2 ybox=$HEIGHTOFF2 & label2="Trace" \ title="PEF: 50,160ms (Note dipping reflector!)" \

7、带通滤波 supef <modeldata maxlag=.04 | supef minlag=.05 maxlag=.16 | sufilter f=5,15,80,100 | suxwigb label1="Time" windowtitle="Final" \ wbox=$WIDTH hbox=$HEIGHT xbox=$WIDTHOFF3 ybox=$HEIGHTOFF2 & label2="Trace" title="Spike, PEF, BandPass" \

例2、 1、数据 susynlv nt=201 dt=0.004 ft=0.0 nxo=1 \ nxs=100 dxs=.025 fxs=0.0 er=0 ob=0 \ v00=2.0 dvdz=0 dvdx=0 fpeak=30 smooth=1 \ ref="1:-1.,.1;7.,.5" \ ref="1:-1.,.5;7.,.2" \ ref="2:-1.,.5;7.,.9" |

sushw key=d2 a=.05 >susyn1.out sufilter <susyn1.out >susyn1.flt 2、加噪声 suaddnoise < susyn1.out sufilter >susyn1.sn=4.flt 3、绘图 suxmovie <susyn1.flt \ 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=99.5 & sn=8 f=4,8,100,125 amps=0,1,1,0 |

suxmovie <susyn1.sn=4.flt \ 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=99.5 &

5、反褶积 sufxdecon <susyn1.sn=4.flt >susyn1.sn=4.fx ntrw=100 suxmovie <susyn1.sn=4.fx \ d1num=.1 d2num=10 \ label1="Time [s]" label2="Trace #"

n1=201 n2=100 d1=.004 d2=1 f1=0 f2=1 clip=30. \ -geometry ${WIDTH}x${HEIGHT}+${WIDTHOFF2}+${HEIGHTOFF1} \ title="Fxdcn over Data" perc=99.5 &

sufxdecon <susyn1.flt >susyn1.flt.fx ntrw=100 suxmovie <susyn1.flt.fx \ d1num=.1 d2num=10 \ label1="Time [s]" label2="Trace #"

n1=201 n2=100 d1=.004 d2=1 f1=0 f2=1 clip=30. \ -geometry ${WIDTH}x${HEIGHT}+${WIDTHOFF2}+${HEIGHTOFF2} \ title="Fxdcn over Signal" perc=99.5 &

7、绘图

suop2 susyn1.flt susyn1.flt.fx op=diff | suxmovie \ n1=201 n2=100 d1=.004 d2=1 f1=0 f2=1 clip=30. \ label1="Time [s]" label2="Trace #" title="Diff: signal - fx(signal)" d1num=.1 d2num=10 \ perc=99.5 & -geometry ${WIDTH}x${HEIGHT}+${WIDTHOFF3}+${HEIGHTOFF1} \

suop2 susyn1.flt susyn1.sn=4.fx op=diff | 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}+${WIDTHOFF3}+${HEIGHTOFF2} \ title="Diff: signal - fx(data)" perc=99.5 &

五.带通滤波 常用命令为 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 locut=1 hicut=1 fstoplo=0.10*(nyq) astoplo=0.05 fpasslo=0.15*(nyq) apasslo=0.95 fpasshi=0.40*(nyq) apasshi=0.95 fstophi=0.55*(nyq) astophi=0.05 verbose=0 dt = (from header) 简单的举个例子, Suplane > data.su Subfilt < data.su > dataout Subxwigb < data.su Subxwigb < dataout 关于 DEMOS 的使用,如下: 进入目录,双击 xbfilt1,xbfilt2. =0 for minimum phase filter =0 for no low cut filter =0 for no high cut filter freq(Hz) in low cut stop band upper bound on amp at fstoplo freq(Hz) in low cut pass band lower bound on amp at fpasslo freq(Hz) in high cut pass band lower bound on amp at fpasshi freq(Hz) in high cut stop band upper bound on amp at fstophi =1 for filter design info time sampling interval (sec)

六.偏移 常用的命令为 sugazmig、sumigps、sumigpspi。使用说明如下:
SUGAZMIG - SU version of Jeno GAZDAG's phase-shift migration for zero-offset data. sugazmig <infile >outfile vfile= [optional parameters] Optional Parameters: dt=from header(dt) or ft=0.0 ntau=nt(from data) dtau=dt(from header) ftau=ft .004 time sampling interval midpoint sampling interval dx=from header(d2) or 1.0

first time sample number of migrated time samples migrated time sampling interval first migrated time sample

tmig=0.0 vmig=1500.0 vfile= verbose=0 tmpdir=

times corresponding to interval velocities in vmig interval velocities corresponding to times in tmig name of file containing velocities verbose = 1 echoes information if non-empty, use the value as a directory path prefix for storing temporary files; else if the the CWP_TMPDIR environment variable is set use its value for the path; else use tmpfile()

SUMIGPS - MIGration by Phase Shift with turning rays sumigps <stdin >stdout [optional parms] Required Parameters: None Optional Parameters: dt=from header(dt) or .004 dx=from header(d2) or 1.0 ffil=0,0,0.5/dt,0.5/dt tmig=0.0 vmig=1500.0 vfile= nxpad=0 ltaper=0 verbose=0 tmpdir= time sampling interval distance between sucessive cdp's

trapezoidal window of frequencies to migrate times corresponding to interval velocities in vmig interval velocities corresponding to times in tmig

binary (non-ascii) file containing velocities v(t) number of cdps to pad with zeros before FFT length of linear taper for left and right edges =1 for diagnostic print if non-empty, use the value as a directory path prefix for storing temporary files; else if the the CWP_TMPDIR environment variable is set use its value for the path; else use tmpfile()

SUMIGPSPI - Gazdag's phase-shift plus interpolation migration for zero-offset data, which can handle the lateral velocity variation. sumigpspi <infile >outfile vfile= [optional parameters] Required Parameters: nz= dz= number of depth sapmles depth sampling interval (Please see Notes below concerning the format of vfile) Optional Parameters: dt=from header(dt) or .004 dx=from header(d2) or 1.0 tmpdir= time sampling interval midpoint sampling interval

vfile= name of file containing velocities

if non-empty, use the value as a directory path

prefix for storing temporary files; else if the the CWP_TMPDIR environment variable is set use its value for the path; else use tmpfile()

关于 DEMOS 的使用,如下: 进入目录$CWPROOT/src/demos/Migration,可见
1.Gazdag(零偏移距数据的相移):PSGazmig—Xgazmig 2.Split—step :Xdemo—Pademo 3.Phase—shift-plus-interpolation: Xdemo—Psdemo 依此顺序,按手册说明操作。 叠后偏移举以下例子 例 1、 1、 尖脉冲数据 WIDTH=420 HEIGHT=300 WIDTHOFF1=0 WIDTHOFF2=430 WIDTHOFF3=860 HEIGHTOFF1=5 HEIGHTOFF2=335 HEIGHTOFF3=665 perc=99 n1=64 n2=32 suspike nspk=3 it1=16 ix1=16 it2=32 ix2=16 it3=48 ix3=16 | sushw key=dt,d2 a=50000,.05 >syndata1 suxmovie <syndata1 perc=$perc \ title="Suspike data"\ label1="Time (sec)" label2="Midpoint (km)" \ -geometry ${WIDTH}x${HEIGHT}+${WIDTHOFF1}+${HEIGHTOFF1} &

2、以sugazmig作偏移 sugazmig <syndata1 tmig=0 vmig=1 | suxmovie perc=$perc \ title="SUGAZMIG Migration"\ label1="Migrated Time (sec)" label2="Midpoint (km)" \ -geometry ${WIDTH}x${HEIGHT}+${WIDTHOFF2}+${HEIGHTOFF1} &

3、以sumigps作偏移 sumigps <syndata1 tmig=0 vmig=1 | suxmovie perc=$perc \ title="SUMIGPS Migration"\ label1="Migrated Time (sec)" label2="Midpoint (km)" \ -geometry ${WIDTH}x${HEIGHT}+${WIDTHOFF3}+${HEIGHTOFF1} &

例 2、 1、 人工数据
susynlv nt=$n1 dt=0.04 ft=0.0 nxo=1 \ nxm=$n2 dxm=.05 fxm=0.0 er=0 ob=1 v00=1.0 dvdz=0 dvdx=0 smooth=1 \ ref="0,.5;1.0,.5;2.,1.0;2.5,1.5;3.0,1.0;4.0,.5;5.0,.5" | sushw key=d2 a=.05 >syndata2 suxmovie <syndata2 perc=$perc \ title="Synthetic Data"\ label1="Time (sec)" label2="Distance (km)" \ -geometry ${WIDTH}x${HEIGHT}+${WIDTHOFF1}+${HEIGHTOFF2} & \

2、以sugazmig作偏移 sugazmig <syndata2 tmig=0 vmig=1 | suxmovie perc=$perc \ title="SUGAZMIG Migration"\ label1="Migrated Time (sec)" label2="Midpoint (km)" \ -geometry ${WIDTH}x${HEIGHT}+${WIDTHOFF2}+${HEIGHTOFF2} &

3、以sumigps作偏移 sumigps <syndata2 tmig=0 vmig=1 | suxmovie perc=$perc \ title="SUMIGPS Migration"\ label1="Migrated Time (sec)" label2="Midpoint (km)" \ -geometry ${WIDTH}x${HEIGHT}+${WIDTHOFF3}+${HEIGHTOFF2} &

叠前偏移的例子 1、产生单炮信号单炮2D叠前偏移 suimp2d nt=512 nshot=1 nrec=200 dgx=8 dsx=8 x0=640 z0=768 c=2000 | sukfilter k=.001,.002,.3,.4 amps=0,1,1,0 > 1shot_spike.su

2、信号图 suximage < 1shot_spike.su title="1 shotgather data " \ xbox=10 ybox=10 wbox=350 hbox=350 &

3、K域重建
suspeck1k2 < 1shot_spike.su | suximage x1beg=-.8 x2beg=-.8 x2end=.8 x1end=.8 \ xbox=400 ybox=10 wbox=350 hbox=350 \ legend=1 cmap=hsv4 \ title="K-domain representation 1 shotgather data" &

4、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 > 10shots_spike.su

5、 绘图
suximage < 10shots_spike.su title="10 shotgathers" \ xbox=10 ybox=400 wbox=350 hbox=350 perc=99&

6、 K域时窗
suspeck1k2 < 10shots_spike.su | suximage x1beg=-2.0 x2beg=-2.0 x2end=2.0 x1end=2.0 \ xbox=400 ybox=400 wbox=350 hbox=350 \ legend=1 cmap=hsv4 \ title="K-domain aperture 10 shotgathers, PSPI" &

7、 背景场速度
makevel nx=200 nz=1024 v000=2000.0 dx=8 dz=1.0 dvdx=0.0 dvdz=0.0 >vfile1

transp <vfile1 n1=200>vfile2

8、 分布傅立叶偏移
sumigpresp < 1shot_spike.su nxo=200 nxshot=1 nz=1024 \ dx=8 dz=1.0 fmax=50 vfile=vfile2 > prestack_mig_sp.su 9、 输出图像 suximage < prestack_mig_sp.su title="Prestack migration with sumigpresp" \ xbox=10 ybox=10 wbox=350 hbox=350 &

10、 K域时窗 suspeck1k2 < prestack_mig_sp.su | suximage x1beg=-.2 x2beg=-.2 x2end=.2 x1end=.2 \ xbox=350 ybox=10 wbox=350 hbox=350 \ legend=1 cmap=hsv4 \ title="K-domain aperture 1 shotgather, sumigpresp" &

11、 理论上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 K-domain aperture, 1 shotgather data" &

12、 背景速度 makevel nx=200 nz=1024 v000=2000.0 dx=8 dz=1.0 dvdx=0.0 dvdz=0.0 >vfile1 transp <vfile1 n1=200>vfile2 13、 有限查分偏移 sumigprefd < 1shot_spike.su nxo=200 nxshot=1 nz=1024 \ dx=8 dz=1.0 fmax=50 vfile=vfile2 > prestack_mig_fd.su 14、 图形输出 suximage < prestack_mig_fd.su title="Prestack migration with sumigprfd" \ xbox=10 ybox=10 wbox=350 hbox=350 &

15、K域时窗 suspeck1k2 < prestack_mig_fd.su | suximage x1beg=-.2 x2beg=-.2 x2end=.2 x1end=.2 \ xbox=350 ybox=10 wbox=350 hbox=350 \ legend=1 cmap=hsv4 \ title="K-domain aperture 1 shotgather, sumigprefd" &

16、理论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 K-domain aperture, 1 shotgather data" &

17、背景速度 makevel nx=200 nz=1024 v000=2000.0 dx=8 dz=1.0 dvdx=0.0 dvdz=0.0 >vfile1 transp <vfile1 n1=200>vfile2 18、傅立叶有限差分 sumigpreffd < 1shot_spike.su nxo=200 nxshot=1 nz=1024 \ dx=8 dz=1.0 fmax=50 vfile=vfile2 > prestack_mig_ffd.su 19、偏移输出 suximage < prestack_mig_ffd.su title="Prestack migration with sumigrpeffd" \ xbox=10 ybox=400 wbox=350 hbox=350 &

20、K域时窗 suspeck1k2 < prestack_mig_ffd.su | suximage x1beg=-.2 x2beg=-.2 x2end=.2 x1end=.2 \ xbox=350 ybox=400 wbox=350 hbox=350 \ legend=1 cmap=hsv4 \ title="K-domain aperture 1 shotgather, sumigpreffd" &

21、理论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 K-domain aperture, 1 shotgather data" &

七.典型数据的产生 主要有以下的命令:susynlv,suplane, suspike。下面举susynlv为例,说明一下。 1.单炮数据
susynlv nt=200 dt=0.04 ft=1 nxo=32 dxo=0.5 fxo=0 nxs=1 dxs=0.5 fxs=0.0 > data1.su & suxwigb < data1.su & suaddnoise < data1.su > datanoise1.su sn=20 noise=gauss & suxwigb < datanoise1.su title="datanoise1" & sunmo < data1.su > nmo1.su vnmo=2000 suxwigb < nmo.su title="nmo1.su" & susort < data1.su > sort1.su & suxwigb < sort1.su title="sort1.su" & sustack < sort1.su > sustack1.su & suxwigb < sustack1.su title="sustack1.su" & sufilter < data1.su > filter1.su | suxwigb <filter1.su & sufilter < data1.su > filter1.su f=10,20 amps=1,2 & suxwigb < filter1.su & sufilter < data1.su > filter1.su f=10,20 amps=1,2 dt=0.008 & &

suxwigb < filter1.su & suacor < data1.su > auto1.su & suxwigb < auto1.su & 2.共中心点数据 susynlv nt=200 dt=0.04 ft=1 nxo=32 dxo=0.5 fxo=0 nxm=1 dxm=0.5 fxm=0 > data2.su & suxwigb < data2.su & suaddnoise < data2.su > datanoise2.su sn=20 noise=gauss & suxwigb < datanoise2.su title="datanoise2" & sunmo < data2.su > nmo2.su vnmo=2000 suxwigb <nmo2.su title="nmo2.su" & susort < data2.su > sort2.su & suxwigb < sort2.su title="sort2.su" & sustack < sort2.su > sustack2.su & suxwigb < sustack2.su title="sustack2.su" & suxwigb < sustack.su title="sustack1.su" & sufilter < data2.su > filter2.su | suxwigb <filter2.su & sufilter < data2.su > filter2.su f=10,20 amps=1,2 & suxwigb < filter2.su & sufilter < data2.su > filter2.su f=10,20 amps=1,2 dt=0.08 & suxwigb < filter2.su & suacor < data2.su > auto2.su & suxwigb < auto2.su & 2.产生脉冲串 suspike

SUSPIKE - make a small spike data set suspike [optional parameters] > out_data_file Creates a common offset su data file with up to four spikes 可选参数: nt=64 ntr=32 dt=0.004 offset=400 nspk=4 ix1= ntr/4 it1= nt/4 ix2 = ntr/4 it2 = 3*nt/4 it3 = nt/4; it4 = 3*nt/4; number of time samples时间采样数 number of traces道集数 time sample rate in seconds时间采样率 offset偏移距 number of spikes脉冲数 trace number (from left) for spike #1第一脉冲的道集数(从左 数) time sample to spike #1 第一脉冲的时间采样 trace for spike #2 time for spike #2 time for spike #3 time for spike #4 。 。 。 。 。 。 for impulse response studies

ix3 = 3*ntr/4; trace for spike #3 ix4 = 3*ntr/4; trace for spike #4


相关文章:
su简单中文使用手册
su简单中文使用手册_计算机软件及应用_IT/计算机_专业资料 暂无评价|0人阅读|0次下载|举报文档 su简单中文使用手册_计算机软件及应用_IT/计算机_专业资料。CWP ...
SU系列使用手册
SU系列使用手册_信息与通信_工程科技_专业资料。UPS产品使用手册目录...su简单中文使用手册 33页 免费 SU使用笔记_刘卫华_v2 24页 2下载券 Vray__For...
RClimDex中文用户手册
RClimDex中文用户手册_天文/地理_自然科学_专业资料。...SU25 夏日日数 日最高气温>25℃的日数 天 GSL ...当前的 RClimDex 版 只包括一个简单的数据质量控制...
HDevelop简体中文使用手册
HDevelop简体中文使用手册_计算机软件及应用_IT/计算机_专业资料。halcon 指导 文档...suaJization>5et Par ameters...>Zoom Visualization > Set Parameters...> ...
SR868C8中文说明书
SR868C8中文说明书 隐藏>> 太阳能分体承压控制仪 SR868C8 操作手册 目录 目录--...星期 MO TU WE TH FR SA SU 星期一 星期二 星期三 星期四 星期五 星期六...
DCS SU-25 操作手册
DCS SU-25 操作手册_计算机软件及应用_IT/计算机_专业资料。LOCK ON 中文讲解...!在模拟飞行游戏中 LOCK ON 的起飞步骤 算是比较简单的, 玩家在实际飞行中...
SU560中文说明书
SU560中文说明书_计算机软件及应用_IT/计算机_专业资料。使用说明书 冰淇淋机 SU-560 ● 承蒙您购买本公司冰淇淋机,在此表示衷心的谢意。 ● 请你在使用之前仔细...
TWS用户手册中文
TWS用户手册中文_交通运输_工程科技_专业资料。交易...版权所有。 Sun、Sun Microsystems、 Su 交易平台 ...445 应用简单移动平均指标 ......
VrayForSketchup渲染插件中文使用手册
简单的说这种渲染模式就是把 Primary GI engine 和 Secondary GI engine 都设置...SU+Podium+渲染插件的详... 4页 免费 中文使用手册 10页 免费 Eclipse 插件...
su8_快捷键大全
su8_快捷键大全_建筑/土木_工程科技_专业资料。如何加快设计速度 SU8 官方中文版快捷键 【编辑】 群组 N 组件 G 隐藏 H 隐藏其他 Shift+H 显示上一次 V ...
更多相关标签:
dhtmlx中文使用手册 | gitlab 中文使用手册 | kibana4 中文使用手册 | jira中文使用手册 | jmeter中文使用手册 | mpdf 中文使用手册 | jenkins中文使用手册 | poco c 中文使用手册 |