Lab 01 - Examine Seismic Data 2
Lab 01 - Examine Seismic Data 2
1. INTRODUCTION
When geophysicists meet with their seismic data for the first time, they carefully must look into the
seismic data itself and its header, which contains information and tabulation of parameters used to
acquire the data. Usually, the data is stored in magnetic tapes or hard disks and is saved in various
standard data formats such as the well-known SEG-Y or Seismic Unix formats.
The seismic data and its header information are stored in the MATLAB data format file called
SeismicData.mat. In this lab session, you are require to examine this data along with its header
information, by using MATLAB. This will assist you in further analysis and processing of this real reflection
seismic data, which will be performed in subsequent lab session.
Bytes
Class
7132752
3303008
double
struct
Attributes
Size
Bytes
Class
1501 x594
1x594
1x594
7132752
3303008
4752
double
struct
double
Attributes
You can explore more fields in a similar way. A MATLAB function called extracting_geometry.m (provided
with the manual) can be used to get useful geometrical information such as the source coordinates (sx &
sy), the receiver coordinates (gx & gy), the source elevations (gz), and others. The following MATLAB code
shows an example of how one can extract certain geometrical acquisition parameters:
1 load ('SeismicData.mat','H')
2 [sx ,sy ,gx ,gy ,shot_gathers ,num_trace_per_sg ,sz ,gz] = extracting_geometry (H);
These are different geometrical variables obtained from the seismic header structure H. We then can use
such variables to plot the number of seismic traces versus the number shot gathers as in the following
code:
1 figure ,stem (shot_gathers ,num_trace_per_sg)
2 xlabel ('Shot gather numbers','FontSize' ,14)
3 ylabel ('Number of traces/shot gather','FontSize' ,14)
4 axis ([0,max(shot_gathers)+1,0,max(num_trace_per_sg)+2])
5 set (gca ,'YMinorGrid','on')
Also, here is another example to plot the number of sources x-axis locations versus the number of traces:
1 figure ,plot(sx,'.')
2 xlabel('Number of traces','FontSize' ,14)
3 ylabel('Sources xaxis locations (ft)','FontSize' ,14)
4 axis tight
5 grid
To plot the receivers x-axis locations versus the number of traces, you need to alter the code above by
calling the receiver parameter. Please plot the source and receivers elevation vs the number of traces.
Figure 1.1 shows different useful plots all of which are obtained from the seismic header variable H. A
final example is plotting the so-called seismic stacking chart (see Figure 1.2) using the MATLAB function
stack_chart.m (provided with the manual) where the vertical axis shows the shot number and the
horizontal axis shows the shot and receiver x-coordinate (we will discuss stacking charts in more details
in lab 5).
1 num_shots=length(shot_gathers);
2 stacking_chart(sx ,gx ,num_shots ,num_trace_per_sg);
Figure 1.1: Various geometrical information plots: (a) the number of traces per shot, (b) the source
elevation profile for each trace and (c) the receiver elevation per seismic trace.
Figure 1.3: Various displays for seismic data shot gather number 8: (a) variable area display, (b)
grayscaled variable density display and (c) colored variable density display. The color bars in (b)
and (c) refers to the amplitude dynamic range of the data.
If we display shot gather number 16 (see Figure 1.4), 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 steps. In chapter 2, we shed more light on this issue. You may also interested in
displaying a group of seismic shot gathers concatenated together using the same function
extracting_shots.m. For example, you can extract shot gathers number 4-6 where you must provide
shot_num=4:6; in line 6 of the above code. Figure 1.5 shows these extracted shot gathers in their variable
area display.
Figure 1.4: Various displays for seismic data similar to those in Figure 1.3 but for shot gather
number 16. Clearly, trace number 31 (at offset 0) requires muting (replacing it by zeros) as will be
discussed later in chapter 2.
5. COMPUTER ASSIGNMENTS
1. Load the seismic data SeismicData.mat.
2. Examine the seismic data header information and obtain the following information: common depthpoint (CDP) values, number of time samples (nt), spatial sampling interval (dx) and time sampling
interval (dt).
3. Display shot gathers 1215 (as the case in Figure 1.5) using the:
a.
b.
c.
d.
Wiggle plotting
Gray-scaled image plotting
Colored image plotting
Wiggle on top of colored plotting.
Comment on the shot gather amplitudes .What do you notice from a trace to another?