How to make a seismic wedge model?

4 visualizaciones (últimos 30 días)
Maria Amr
Maria Amr el 28 de En. de 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

Respuesta aceptada

Daniel Juliandro Lumban Gaol
Daniel Juliandro Lumban Gaol el 11 de Jun. de 2021
same, your code isnt working. do you guys have any suggestion?

Más respuestas (0)

Categorías

Más información sobre Seismology en Help Center y File Exchange.

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by