基于matlab的地震正演模型实现

内涵两个matlab程序和论文，模型为三层，含油模型和不含油模型。根据褶积理论 ,结合地质模型 ,在 MATLAB 环境中编写程序实现了地震正演。为了最大可 能的与实际情况相符合 ,地震模型中使用的子波是稳 定可实现的子波。

n = 3000; %采样点数
dx = 50; %道间距 
dt = 0.002; %采样间隔
nl = 3; %界面数 
m = 80; %道数 
v0 = [1000 1500 2000 2500] ; %第 1 层速度 
v1 = [1000 1500 1300 2500] ; %第 2 层速度 
v2 = [1000 1500 2000 2500] ; %第 3 层速度 
h（1 1） = 800;
h（2 1） = 1200; 
h（3 1） = 1800; %该模型为各层平行 无倾角 
xmax = （m/ 2 + 1）*dx ; 
ymax = n * dt ; %坐标范围 

p = 2; %图象显示数字 以下部分为子波采样 
a0= 200;
f = 20;
nw = 60;
b = 30; %子波参数 
tt = 0:dt :（nw - 1）*dt ;
wb =a0 * sin（2 * pi *f * tt） .* exp（ - b * tt） ; %视速 度 
hva0= h（3  :） .* v0（1:nl） ;
hvb0= h（3  :） ./ v0（1 :nl） ; 
hva1= h（3  :） .* v1（1:nl） ;
hvb1= h（3  :） ./ v1（1 :nl） ;
hva2= h（3  :） .* v2（1:nl） ;
hvb2= h（3  :） ./ v2（1 :nl） ;

for i = 2:nl 
    hva0（i） = hva0（i） + hva0（i - 1） ;
    hvb0（i） = hvb0（i） + hvb0（i - 1） ;
    hva1（i） = hva1（i） + hva1（i - 1） ;
    hvb1

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2020-08-20 11:53  地震正演\
     文件        3102  2020-08-20 10:45  地震正演\main.m
     文件        2630  2020-08-20 11:44  地震正演\main1.m
     文件       30182  2020-08-20 11:49  地震正演\不含油模型.png
     文件       30958  2020-08-20 11:50  地震正演\含油模型.png
     文件       73707  2020-08-20 11:52  地震正演\地质模型.png