Annual Report 2000 - WIT
Annual Report 2000 - WIT
Annual Report 2000 - WIT
You also want an ePaper? Increase the reach of your titles
YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.
22<br />
For synthetic data, it is easy to pick these events, see Figure 2 (a). As there are<br />
no conflicting events and the signal-to-noise (S/N) ratio is high, the picker follows the<br />
same event over the entire ZO section. The synthetic data have been generated using<br />
a zero-phase Ricker wavelet with a peak frequency of 15 Hz. Picking the maximum<br />
amplitude yields the two-way traveltime from the source to the interface and back to<br />
the receiver at the surface. If seismic events in a real data time section, Figure 2 (b),<br />
have to be followed continuously, problems occur. Usually, the S/N ratio decreases<br />
with increasing traveltime. Events at small traveltimes of the simulated ZO section<br />
might be disturbed by mute zones. Regions with conflicting dips or where scattering<br />
prevails, interrupt the automatic picking, see Figure 2 (b) in the vicinity of CMP no.<br />
550. Thus, half-automatic picking is requested in order to overcome these problems.<br />
A little error is included while picking the maximum amplitude. In the synthetic<br />
and real data-set the actual local maximum cannot always be picked because of the<br />
discrete time sampling. But the picked maximum is very close to it, hence, the error<br />
is small. The traveltime that corresponds to the maximum amplitude of the wavelet<br />
differs slightly from the traveltime that actually represents the spatial location of the<br />
impedance contrast in the time domain. To prevent those errors, one suggestion is to<br />
pick the time sample where the wavelet starts to separate from the noise but that is<br />
even more difficult because one has to define a threshold for the noise.<br />
It is advantageous to pick as much events as possible because the more events are<br />
picked the more interfaces can be inverted and the more precise the model will be.<br />
Another reason is that the inversion algorithm provides only constant velocities within<br />
the layers along the traced rays. Hence, it is preferable to have more picked events<br />
with smaller time intervals in order to construct a more complex model.<br />
SMOOTHING BY MEANS OF STATISTICAL METHODS<br />
The inversion algorithm requires smooth CRS stack attributes. The deviation of the<br />
angle ¡ from trace to trace along an event is small and, thus, the data of the emergence<br />
angles do not require much smoothing. In contrast, the radius of the NIP wavefront<br />
¢¥£3¦6¨ often varies strongly from trace to trace, see Figure 3 (a). To perform a more<br />
stable inversion, the data have to be smoothed.<br />
There are several filters available for smoothing purposes. The smoothing is done<br />
within a pre-defined window of wv G length . This yields a symmetrical window<br />
around the picked time sample that can be set independently in time and/or trace direction.<br />
Here, possible filter choices are:<br />
1. The arithmetic mean, x<br />
<br />
, is the normalised sum of all values 1y/z {| within the