How to make a seismic wedge model?
3 views (last 30 days)
Show older comments
Maria Amr
on 28 Jan 2021
Answered: Daniel Juliandro Lumban Gaol
on 11 Jun 2021
I am trying to exam your code "Wedge Model" and it seems there is an error and it is not working. I would truly appreciate if you direct me how it is working for an arbitrary reflectivity series. Thank you!
% Copyright: 2018 - Teknik Geofisika, Universitas Pertamina
% URL: https://sites.google.com/site/metkomup/programming
% Cite: Anugerah, Nisfu; Ginting, Gamaliel Rhema; Wicaksono, Gigih Aji; Salsabila, Alda;
% Subakti, Puguh Ari; Syahputra, Loris Alif (2018): Membuat Model Sintetik untuk Model Pembajian.
% figshare. https://doi.org/10.6084/m9.figshare.5946691.v1
% Ricker Wavelet
dt = 0.004;
f = 15;
phi = 3.14;
t0 = 0;
t = [-16:dt:16];
r = (1-2*(phi*f*dt.*(t-t0)).^2)./exp((phi*f*dt.*(t-t0)).^2);
y = decimate(r,120);
% Plotting Hasil Convolution Hasil Membaji
figure(1)
plot(y)
title('Wavelet');
for i = 1:5:1000
c = 1:315;
ff = 1:249;
b = conv(Reflect(:,i),y);
amax(i) = max(b);
figure(1)
plot(10*b+(i*0.25),c);
figure(1)
hold on
title('Wedge Model f=15Hz')
set(gca,'ydir','reverse')
axis([100 250 0 250])
end
% Tunning Thickness
figure(2)
for i = 1:5:1000
plot(i/4,amax(i),'*r');hold on
axis([100 250 0 0.6])
end
0 Comments
Accepted Answer
Daniel Juliandro Lumban Gaol
on 11 Jun 2021
same, your code isnt working. do you guys have any suggestion?
0 Comments
More Answers (0)
See Also
Categories
Find more on Seismology in Help Center and File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!