SU克希霍夫深度偏移脉冲响应试验

 

  SU克希霍夫深度偏移脉冲响应试验

    有不少网友都想要一些实际的例子，作为学习SU的模版，其实SU自带的demos里已经有很多实例，足够大家去钻研学习。

     SU里有一些克希霍夫叠前，叠后深度偏移的例子，但是没有脉冲响应实例，我以前做过一些试验，贡献出来，与大家相互交流。

一、合成单脉冲数据

编写脚本文件 Makedata，内容如下：

#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 \

        ix1=$ix1 it1=$it1 ix2=$ix2 it2=$it2 ix3=$ix3 \

        it3=$it3 ix4=$ix4 it4=$it4 |
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   nt=500
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 &

偏移结果如下：

 

   

    只要走时表计算好，输出参数范围小于等于输入参数范围，就不会有问题。时间原因不细说了，大家自己体验，有什么问题再交流。