LAB 02 - MUTE SHOT NO 16

MATLAB SCRIPT LAB 2.2

clc, clear, close all;
load SeismicData_C;

%% Extract shot number 16
shot_num=16;
p=0;
[Dshot,dt,dx,t,offset]=extracting_shots(D,H,shot_num,p);
Dg=D;

%% Display before muting
scale=1;
figure(5);
mwigb(Dshot,scale,offset,t)
xlabel('Offset(ft)','FontSize',14)
ylabel('Time(s)','FontSize',14)
title('Before muting','FontSize',16)

Figure 1: Shot gather number 16, before muting

%% Mute trace
[i,j]=find(Dg==max(max(Dg)));
Dg(:,j)=0;
Dgz=Dg;

%% Display After Muting
p=0
[Dshot,dt,dx,t,offset]=extracting_shots(Dgz,H,shot_num,p);
scale=1;
figure(6);
mwigb(Dshot,scale,offset,t)
xlabel('Offset(ft)','FontSize',14)
ylabel('Time(s)','FontSize',14)
title('After muting','FontSize',16)

Figure 2: Shot gather number 16, after muting

You already have notice that some high noisy amplitude traces particularly in trace number 31 shot number 16. This is due to the use of bad geophone. Such recorded trace needed to be replace by a zero trace. In MATLAB we can simply use the following step to mute the trace on the complete dataset:-
​
[i,j]=find(Dg==max(max(Dg)));
Dg(:,j)=0;
​