叠后地震剖面制作

1 地震记录制作

首先得到一个合成地震记录，传送门

1.1 模型绘制

• $1000\times 1000$ 的四层模型，从上到下各地层纵波速度分别为1000m/s，2000m/s，3000m/s，4000m/s。

1.2 观测系统建立

• 点击Framework，可得到这个窗口。
  

• 点击Source，设置震源。
  

  • Horizontal line：水平放炮，Vertical line:：井下放炮。
  • Number：放炮数。
  • Interval：炮检距。

• 点击Signal，设置地震波。
  

  • 这里设置的是30Hz的雷克子波。

• Obeservation，检波器设置。
  

  • 水平地面0到1000米，每10米放置一个检波器，即总共101个检波器；采样时间2s，采样时间间隔0.002，总共采样1000次。

1.3 声波方程模拟

• 点击Run图标。
  
  • Invisible 设置为不可见地表。
  • Wave Eqution 设置声波方程。
  • 然后点击run。
• 得到我们的炮集记录和波场快照。
  
  若图像不明显，可调节增益。
  

2 叠后剖面制作

• 初至切除
• 点击Edit，Muting，Upper Line，对初至波（反射波）进行切除。
  
• 右键结束。
• 再次点击Upper line。
  
• 点击Run, Pre-processing procedures, Gathering(GATHER)。
  
  
  • 参数需要记住。
• 然后对道集进行速度拾取。点击Edit, Velocity, Create Velocity Spectra, 绘制速度谱。
  
  
  • 设置开始速度和结束速度，以及采样步长。
    
• 拾取完毕后，进行保存。
  
• 然后对炮集记录进行NMO动校正。
  -
  • 把刚刚的参数填入。
    -载入速度。
    
• 叠加。
• 

• 回到我们的模型，点击File, Convert to, TGR(Tesseral grid)
  
• 再对我们的叠加剖面进行速度载入。
  
• 最后得到柯希霍夫叠后偏移
  
• 用Fimage打开。
  

---

https://www.bilibili.com/video/BV1f64y1U7gj?spm_id_from=333.1007.top_right_bar_window_custom_collection.content.click