How to make a seismic wedge model?
1 Ansicht (letzte 30 Tage)
Ältere Kommentare anzeigen
Maria Amr
am 28 Jan. 2021
Beantwortet: Daniel Juliandro Lumban Gaol
am 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 Kommentare
Akzeptierte Antwort
Daniel Juliandro Lumban Gaol
am 11 Jun. 2021
same, your code isnt working. do you guys have any suggestion?
0 Kommentare
Weitere Antworten (0)
Siehe auch
Kategorien
Mehr zu Seismology finden Sie in Help Center und File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!