陈桂廷的个人博客分享 http://blog.sciencenet.cn/u/jinshi2015 Reproducible Research Philosophy

博文

简单模型kirchhoff叠前深度偏移实验(SU Shell)

已有 4836 次阅读 2016-7-14 20:55 |系统分类:科研笔记| Shell

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 完成整个实验过程,结果如下:

正演记录:


偏移结果:





http://blog.sciencenet.cn/blog-2834901-990702.html

上一篇:Seismic Unix 安装(2016年)
下一篇:Reproducible research philosophy

0

该博文允许注册用户评论 请点击登录 评论 (0 个评论)

数据加载中...

Archiver|手机版|科学网 ( 京ICP备07017567号-12 )

GMT+8, 2021-6-13 10:08

Powered by ScienceNet.cn

Copyright © 2007- 中国科学报社

返回顶部