#Creates a common offset su data file with up to four spikes
# for impulse response studies
set -x
nt=200 ntr=301 dt=0.001 offset=0 nspk=4 ix1=151 it1=30 ix2=151 it2=60
ix3=151 it3=90 ix4=151 it4=120
suspike nt=$nt ntr=$ntr dt=$dt offset=$offset nspk=$nspk \
sushw key=sx,gx,d2 a=0,0,1 b=1,1,0 >data
数据中一共有四个点,一定要在道头置sx,gx。
用suximage<data &
二、生成速度文件
makevel nx=301 dx=1 nz=201 dz=1 v000=2000 >vfile
ximage< vfile n1=201 legend=1 &
速度图如下:
三、射线追踪计算走时表
先编写参数文件 rayt2d.par
dt=0.001
fz=0 nz=201 dz=1
fx=0 nx=301 dx=1
fxo=0 nxo=301 dxo=1
fzo=0 nzo=201 dzo=1
fxs=0 nxs=101 dxs=3
aperx=150
fa=-90 na=90 da=2 amax=90
fac=0.01 ms=10 ek=1 npv=0
jpfile=jpfile.ray
tfile=tfile
然后
rayt2d <vfile par=rayt2d.par
计算出走时表tfile。这是克希霍夫偏移中最关键的一步,输出的范围一定不要超过速度范围。
四、偏移
编写kdmig.par
fzt=0 nzt=201 dzt=1
fxt=0 nxt=301 dxt=1
fs=0 ns=101 ds=3
fxo=0 nxo=301 dxo=1
fzo=0 nzo=201 dzo=1
aperx=150 dxm=1 ntr=301 ls=0
off0=0 noff=1 doff=20
v0=2000 dvz=0
offmax=500 angmax=90 fmax=250
jpfile=jpfile.kd
ttfile=tfile
然后用下边命令:
sukdmig2d<data par=kdmig.par >kd.data &
suximage <kd.data &
偏移结果如下: