A prestack seismic inversion
暂无评分
资源描述
叠前三参数反演横波速度、纵波速度、密度。
%*************************************************************************
% 反射系数、子波矩阵、合成理论记录
%*************************************************************************
ang=[5,15,25].*pi/180; %角度(弧度制)
anglen=length(ang);
r=rflecoe(vp,vs,den,ang,'aki_richard'); %两种反射系数的计算 'zoeppritz'
Ks=60;f0=40;dt=0.001; %子波长度、主频、采样率
wvletMatrix=waveletMatrix(f0,Ks,dt,length(r)); %子波矩阵
syndat=synrec(wvletMatrix,Ks,r); %合成理论记录
%*************************************************************************
% 迭代计算
%*************************************************************************
fori=1:3
ir=rflecoe(ivp,ivs,iden,ang,'aki_richard');
&nb