MATLAB SCRIPT LAB 1.2
load SeismicData_C.matwhos
load ('SeismicData_C.mat','H')
[sx,sy,gx,gy,shot_gathers,num_trace_per_sg,sz,gz] = extracting_geometry(H);
%% Display wiggle shot num 16
load SeismicData_C.mat
shot_num=16;
p=0;
[Dshot,dt,dx,t,offset]=extracting_shots(D,H,shot_num ,p);
scale=1;
Figure 1: variable area display shot number 16
%% Display grey scaled variable density display
mwigb(Dshot,scale,offset,t)
xlabel('Offset(ft)','FontSize',14)
ylabel('Time(s)','FontSize',14)
Figure 2: grey scaled density display shot no 16
%% Display color variable density display
figure,simage_display(Dshot,offset,t,0)
xlabel('Offset(ft)','FontSize',14)
ylabel('Time(s)','FontSize',14)
Figure 3: colored variable density display shot no 16
Figure 4: wiggle on top of colored plotting area display shot gather no 16
In shot gather number 16, notice that trace number 31 amplitudes are increasing as time increases. This may require editing by muting this trace since it will affect the subsequent processing step. In the next section, we will extract more shot gathers to enhance the image or to appear the noise trace.
No comments:
Post a Comment