SU shell

分成四个文件：

process  实现流程，包括创建速度模型，正演共炮数据，射线追踪计算走时表，kirchhoff深度偏移

input  速度模型参数

kdmig.par   偏移参数

rayt2d.par 射线追踪参数

把参数单独放在成3个文件中读取，简化process中的代码。

实现步骤：

1.vi process

 输入如下代码并保存

set -v

   nz=51 dz=50 fz=.0  labelz="Depth (m)"

   nx=71 dx=50 fx=0.0  labelx="Distance (m)"

   ninf=0  npmax=201

   unif2 <input >vfile  ninf=$ninf  npmax=$npmax \

   nz=$nz dz=$dz fz=$fz nx=$nx dx=$dx fx=$fx \

   v00=1500  dvdz=0.8

   set -x

   nt=725 dt=0.004 ft=0.0 tmin=0.2 nxs=6 dxs=0.3 fxs=0.05

   fxo=0.1 dxo=0.025 nxo=60

   fpeak=30 er=0 ls=1

   susynlv  nt=$nt dt=$dt ft=$ft nxs=$nxs dxs=$dxs fxs=$fxs  fxo=$fxo\

   dxo=$dxo nxo=$nxo tmin=$tmin v00=1.5 dvdx=0.0 dvdz=0.8 \

   fpeak=$fpeak er=1 ls=1

   ref="1:0,0.75;0.5,0.75;1,1.5;1.5,1.5,1.75,1.2;3,1.2">data

   rayt2d <vfile par=rayt2d.par

   sufrac < data phasefac=.25  > data.tmp

   sukdmig2d < data.tmp par=kdmig.par>kd.data   &

   suximage<kd.data perc=99.5 cmap=hsv5         &

   exit

2. vi input

   1       0
   200     0
   1      -99999
   

3. vi rayt2d.par

  dt=0.004   nt=501
  fz=0 nz=51 dz=50
  fx=0 nx=71 dx=50
  fxo=0 nxo=71 dxo=50
  fzo=0 nzo=51 dzo=50
  fxs=0 nxs=36 dxs=100
  aperx=2500
  fa=-75 na=76 da=2 amax=75
  fac=0.01 ms=10 ek=1 npv=0
  jpfile=jpfile.ray
  tfile=tfile

4. vi kdmig.par

 fzt=0 nzt=51 dzt=50
 fxt=0 nxt=71 dxt=50
 fs=0 ns=36 ds=100
 fxo=0 nxo=141 dxo=25
 fzo=0 nzo=251 dzo=10
 aperx=2000 dxm=12.5 ntr=1000 ls=1
 off0=0 noff=1 doff=200
 v0=1500 dvz=0.0
 offmax=4000 angmax=60 fmax=55
 jpfile=jpfile.kd
 ttfile=tfile

4个文件创建完成后输入：chmod +x process 设置执行权限。

之后执行：./process 完成整个实验过程，结果如下：

正演记录：

偏移结果：