Digital Processing Of Shallow Seismic Refraction Data With The Refraction Convolution Section by Derecke Palmer M Sc A Thesis Submitted in Fulfillment of the Requirements for the Degree of Doctor of Philosophy School of Geology, The University of New South Wales, Sydney, Australia. September, 2001
256
Embed
Digital processing of shallow seismic refraction data with the refraction convolution section
This document is posted to help you gain knowledge. Please leave a comment to let me know what you think about it! Share it to your friends and learn new things together.
Transcript
Digital Processing
Of
Shallow Seismic Refraction Data
With
The Refraction Convolution Section
by
Derecke Palmer M Sc
A Thesis Submitted in Fulfillmentof the Requirements for the Degree of
Doctor of Philosophy
School of Geology,The University of New South Wales,
Sydney, Australia.
September, 2001
2
Declaration of Originality
I hereby declare that this submission is my own work and to the best of my
knowledge it contains no materials previously published or written by another
person, nor material which to a substantial extent has been accepted for the
award of any other degree or diploma at UNSW or any other educational
institution, except where due acknowledgement is made in the thesis. Any
contribution made to the research by others, with whom I have worked at UNSW
or elsewhere, is explicitly acknowledged in the thesis.
I declare that the intellectual content of this thesis is the product of my own work,
except to the extent that assistance from others in the project’s design and
conception or in style, presentation and linguistic expression is acknowledged.
Derecke Palmer
26 September, 2001
3
Abstract
The refraction convolution section (RCS) is a new method for imaging shallow
seismic refraction data. It is a simple and efficient approach to full trace
processing which generates a time cross-section similar to the familiar reflection
cross-section. The RCS advances the interpretation of shallow seismic refraction
data through the inclusion of time structure and amplitudes within a single
presentation.
The RCS is generated by the convolution of forward and reverse shot records.
The convolution operation effectively adds the first arrival traveltimes of each pair
of forward and reverse traces and produces a measure of the depth to the
refracting interface in units of time which is equivalent to the time-depth function
of the generalized reciprocal method (GRM).
Convolution also multiplies the amplitudes of first arrival signals. To a good
approximation, this operation compensates for the large effects of geometric
spreading, with the result that the convolved amplitude is essentially proportional
to the square of the head coefficient. The signal-to-noise (S/N) ratios of the RCS
show much less variation than those on the original shot records.
The head coefficient is approximately proportional to the ratio of the specific
acoustic impedances in the upper layer and in the refractor, where there is a
reasonable contrast between the specific acoustic impedances in the layers. The
convolved amplitudes or the equivalent shot amplitude products can be useful in
resolving ambiguities in the determination of wavespeeds.
4
The RCS can also include a separation between each pair of forward and
reverse traces in order to accommodate the offset distance in a manner similar to
the XY spacing of the GRM. The use of finite XY values improves the resolution
of lateral variations in both amplitudes and time-depths.
Lateral variations in the near-surface soil layers can affect amplitudes thereby
causing “amplitude statics”. Increases in the thickness of the surface soil layer
correlate with increases in refraction amplitudes. These increases are
adequately described and corrected with the transmission coefficients of the
Zoeppritz equations. The minimum amplitudes, rather than an average, should
be used where it is not possible to map the near surface layers in detail.
The use of amplitudes with 3D data effectively improves the spatial resolution of
wavespeeds by almost an order of magnitude. Amplitudes provide a measure of
refractor wavespeeds at each detector, whereas the analysis of traveltimes
provides a measure over several detectors, commonly a minimum of six. The
ratio of amplitudes obtained with different shot azimuths provides a detailed
qualitative measure of azimuthal anisotropy.
Dip filtering of the RCS removes “cross-convolution” artifacts and provides a
convenient approach to the study of later events.
The RCS facilitates the stacking of refraction data in a manner similar to the CMP
methods of reflection seismology. It can significantly improve S/N ratios.
The RCS is a simple extension of the GRM, which in turn is a generalization from
which most of the standard refraction inversion methods can be derived. The
RCS advances refraction interpretation through the inclusion of time structure
and amplitudes within a single presentation, which is similar to seismic reflection
data. Accordingly, the RCS facilitates the application of current seismic reflection
acquisition, processing and interpretation technology to refraction seismology.
5
Acknowledgements
This work would not have been possible without the support and encouragement
of my supervisor Geoff Taylor, and our head of school, Colin Ward. My focus on
the thesis in the last few years has resulted in some of my academic duties
receiving less than my full attention.
Much of the work for this thesis was carried out between 4:00 am and 6:00 am in
the morning, and it resulted in a number of innocent victims. My wife Coori, and
our two sons, Evan and Heath have had to accommodate an often sleep-
deprived out-of-sorts partner or parent on more than one occasion.
The processing of this and other refraction data has been made possible by
Seismic Un*x developed by the Centre for Wave Propagation Studies at the
Colorado School of Mines. My sincere appreciation to John Stockwell and the
late Jack Cohen for its development, and to Ken Larner for introducing me to SU.
Jacques Jenny of W_Geosoft has generously provided a copy of Visual_SUNT.
Much of the data were acquired when I was an employee of the Geological
Survey of New South Wales. The data for the Mt Bulga 3D survey were acquired
with the assistance of Ross Spencer during a week in the spring of 1986 which
rapidly turned cold and damp. My memory of the survey is of two bedraggled
geophysicists who had forgotten their wet weather clothing wallowing in ankle
deep mud and becoming increasingly frustrated with a temperamental drill rig.
Ian Grierson of Encom Technologies demuxed many of the older field tapes.
6
Contents
Declaration of Originality ________________________________________ 2Abstract ______________________________________________________ 3Acknowledgements_____________________________________________ 5Contents______________________________________________________ 6
1.1 - Recent Innovations in Reflection Seismology ___________________ 101.2 - Recent Innovations in Shallow Refraction Seismology ____________ 111.3 - Digital Processing with the Refraction Convolution Section_________ 141.4 – Outline of Thesis _________________________________________ 191.5 - References______________________________________________ 21
Chapter 2______________________________________________________ 24Inversion of Shallow Seismic Refraction Data – A Review ____________ 24
2.1 - Summary _______________________________________________ 242.2 - Introduction _____________________________________________ 252.3 - Field Data Requirements ___________________________________ 262.4 - Undetected Layers ________________________________________ 272.5 - Incomplete Sampling of Each Layer __________________________ 272.6 - Implications for Model-Based Methods of Inversion ______________ 282.7 - Anisotropy ______________________________________________ 302.8 - The Need to Employ Realistic Models for Refraction Inversion______ 302.9 - The Large Number of Refraction Inversion Methods ______________ 312.10 - Wavefront Reconstruction Methods__________________________ 312.11 - The Intercept Time Method ________________________________ 322.12 - The Reciprocal Methods __________________________________ 332.13 - Data Processing in the Time Domain_________________________ 332.14 - Accommodation of the Offset Distance with Refraction Migration ___ 352.15 - Using Refraction Migration to Recognize Artifacts_______________ 362.16 - Non-uniqueness in Determining Refractor Wavespeeds __________ 372.17 - Fundamental Requirements for Refraction Inversion_____________ 38References __________________________________________________ 39
Chapter 3______________________________________________________ 47Imaging Refractors with the Convolution Section___________________ 47
7
3.1 - Summary _______________________________________________ 473.2 - Introduction _____________________________________________ 483.3 - The Large Variations in Signal-to-Noise Ratios with Refraction Data _ 503.4 - Full Trace Processing Of Refraction Data ______________________ 553.5 - Imaging The Refractor Interface Through The Addition of Forward AndReverse Traveltimes __________________________________________ 583.6 - The Addition of Traveltimes With Convolution ___________________ 613.7 - The Effects of Geometrical Spreading on the Convolution SectionAmplitudes __________________________________________________ 653.8 - Effects Of Refractor Dip On Convolution Amplitudes______________ 693.9 - Conclusions _____________________________________________ 703.10 - References_____________________________________________ 72
Chapter 4______________________________________________________ 75Starting Models For Refraction Inversion__________________________ 75
4.1 - Summary _______________________________________________ 754.2 - Introduction _____________________________________________ 764.3 - Inversion Of A Two Layer Model With The GRM Algorithms________ 784.4 - Time Differences Between Starting Models_____________________ 834.5 - Agreement Between Starting Models And Traveltime Data_________ 864.6 - Discussion ______________________________________________ 874.7 - Conclusions _____________________________________________ 894.8 - References______________________________________________ 90
Chapter 5______________________________________________________ 93Resolving Refractor Ambiguities With Amplitudes __________________ 93
5.1 - Summary _______________________________________________ 935.2 - Introduction _____________________________________________ 945.3 - Amplitude and Wavespeed Relationships ______________________ 955.5 - Mt Bulga Case History _____________________________________ 975.5 - Conclusions ____________________________________________ 1045.6 - References_____________________________________________ 106
Chapter 6_____________________________________________________ 107Efficient Mapping Of Structure And Azimuthal Anisotropy With ThreeDimensional Shallow Seismic Refraction Methods _________________ 107
6.1 - Summary ______________________________________________ 1076.2 - Introduction ____________________________________________ 1086.3 - Data Processing With The GRM ____________________________ 1106.4 - Survey Details __________________________________________ 1116.5 - Analysis of the In-line Traveltime Data________________________ 1136.6 - Analysis of the In-line Amplitude Data ________________________ 1216.7 - Analysis of the Cross-line Traveltime Data ____________________ 1246.8 - The Cross-line Amplitude Data _____________________________ 1286.9 - Discussion and Conclusions _______________________________ 1326.10 - References____________________________________________ 134
8
Chapter 7_____________________________________________________ 137Effects Of Near-Surface Lateral Variations On Refraction Amplitudes _ 137
7.1 - Summary ______________________________________________ 1377.2 - Introduction ____________________________________________ 1387.3 - Traveltime Results _______________________________________ 1397.4 - Effects of Near-surface Lateral Variations on Amplitudes _________ 1447.5 - Relationships Between Amplitudes and Refractor Wavespeeds ____ 1517.6 - Discussion and Conclusions _______________________________ 1537.7 - References_____________________________________________ 155
Chapter 8_____________________________________________________ 157Enhancement of Later Events in the RCS with Dip Filtering _________ 157
8.1 - Summary ______________________________________________ 1578.2 - Introduction ____________________________________________ 1588.3 - Generation of Useful Events and Artifacts in the RCS____________ 1598.4 - Removal of Cross-convolution Artifacts with Dip Filtering _________ 1638.5 - Times for Near-surface Events in the Uncorrected RCS __________ 1668.6 - Near-surface Wavespeeds from the Uncorrected RCS ___________ 1688.7 - Conclusions ____________________________________________ 1728.8 - References_____________________________________________ 172
Chapter 9_____________________________________________________ 173Stacking Seismic Refraction Data in the Convolution Section________ 173
9.1 - Summary ______________________________________________ 1739.2 - Introduction ____________________________________________ 1749.3 – The Cobar Stacked RCS Section ___________________________ 1769.4 - The Static Geophone Spread_______________________________ 1829.4 - Continuous Acquisition of Shallow Seismic Refraction Data _______ 1839.5 – Determination of Fold with RCS Data ________________________ 1859.6 - Discussion and Conclusions _______________________________ 1869.7 - References_____________________________________________ 188
Chapter 10____________________________________________________ 190Discussion and Conclusions ___________________________________ 190
10.1 - Shallow Refraction Seismology for the New Millenium: A PersonalPerspective_________________________________________________ 19010.2 - Conclusions ___________________________________________ 193
Appendix 1 ___________________________________________________ 198Comments on “A brief study of the generalized reciprocal method andsome of the limitations of the method” by Bengt Sjögren.___________ 198
A.1 - Introduction ____________________________________________ 198A.2 - The Use of Average Wavespeeds___________________________ 199A.3 - The Similarities Between The GRM and Sjogren’s Approach ______ 201A.4 - Recognizing And Defining Narrow Zones With Low Wavespeeds InRefractors__________________________________________________ 203
9
A.5 - Use Of Alternative Presentations And Amplitudes For DeterminingWavespeeds In Refractors _____________________________________ 205A.6 - A Systematic Approach With The GRM_______________________ 211A.7 - The Need To Promote Innovation In Shallow Refraction Seismology 212A.8 - References ____________________________________________ 213
Appendix 2 ___________________________________________________ 216Model Determination For Refraction Inversion ____________________ 216
A.1 - Summary ______________________________________________ 216A.2 - Introduction ____________________________________________ 217A.3 - Model and Inversion Strategies _____________________________ 219A.4 - Single Layer Constant Wavespeed Inversion Model_____________ 226A.5 - Two Layer Constant Wavespeed Inversion Model ______________ 229A.6 - Two Layer Wavespeed Reversal Inversion Model ______________ 231A.7 - The Evjen Inversion Model ________________________________ 232A.8 - Transverse Isotropy Inversion Model_________________________ 238A.9 - Errors Related to the Optimum XY Value _____________________ 241A.10 - Discussion and Conclusions ______________________________ 244A.11 - References ___________________________________________ 247A.12 - Appendix: Definition of Variable Wavespeed Media with the GRM 250
Table 1: Summary of amplitude products and wavesppeds.
7.6 - Discussion and Conclusions
This case history provides another good example of the correlation between
head wave amplitude products and the ratios of the wavespeeds. It is complex
with many large variations in depths as well as wavespeeds in both the refractor
and the layer above. Furthermore, it qualitatively confirms the importance of
densities on head wave amplitudes.
154
The case history also provides a valuable insight into the importance of near-
surface variations and geophone coupling on the measured refraction
amplitudes.
Between stations 50 and 72 where there is no outcrop, the amplitudes and the
amplitude products are essentially a function of the wavespeeds in the refractor
and the layer above. However, at station 56, there is an increase in amplitudes
which correlates with an increase in traveltimes and time-depths. These results
indicate that the amplitude variation is related to the near-surface layering, rather
than to the coupling of the geophones with the ground.
The results for the region between stations 26 and 50, where there was
extensive outcrop, support this interpretation, because all of the amplitude
anomalies can be associated with traveltime anomalies.
The presentation of both amplitude products and time-depths for a range of XY
values from zero to more than the optimum value, provides a convenient and
effective method for recognizing near-surface anomalous zones of limited lateral
extent.
The increases in amplitudes are compatible with the transmission coefficients of
the Zoeppritz equations. As the seismic signal approaches the surface from the
refractor, there is an increase in seismic amplitude where there is another layer
with lower wavespeed and or density.
In general, this change in amplitude can be ignored when there are several
continuous layers above the refractor, because the same increase in amplitudes
occurs at each detector. In these situations, the amplitudes are adequately
described with the head coefficients, together with a geometric spreading factor.
155
Where there are lateral changes in the surface layers, such as the irregular
development of a surface soil layer, there can be large variations in amplitudes
by a factor of between 1 and 2. If these layers have sufficient lateral extent so
that they can be mapped, such as the region between stations 26 and 28, then
an approximate correction factor can be computed with the transmission
coefficients of the Zoeppritz equations.
However, this is not always possible. Under these circumstances, the minimum
amplitudes are probably the most representative.
7.7 - References
Bycroft, G. N., 1956, Forced vibrations of a rigid plate on a semi-infinite elastic
space: Roy. Soc. London, 248, 327-368.
Fail, J. P., Grau, G., and Lavergne, M., 1962, Couplage des sismographes avec
le sol: Geophys. Prosp., 10, 128-147.
Lamer, A., 1970, Couplage sol-geophone: Geophys. Prosp., 18, 300-319.
Palmer, D., 1980, The generalized reciprocal method of seismic refraction
interpretation: Society of Exploration Geophysicists.
Palmer, D., 1986, Refraction seismics - the lateral resolution of structure and
seismic: Geophysical Press.
Palmer, D., 2001a, Imaging refractors with the convolution section: Geophysics
66, 1582-1589.
156
Palmer, D., 2001b, Resolving refractor ambiguities with amplitudes: Geophysics
66, 1590-1593.
Pieuchot, M., 1984, Seismic instrumentation: Geophysical Press.
Polak, E. J., 1969, Attenuation of seismic energy and its relation to the properties
of rocks: Ph D thesis, University of Melbourne, p4.7-4.9.
157
Chapter 8
Enhancement of Later Events in theRCS with Dip Filtering
8.1 - Summary
Later events, which occur in the shot records, are also treated in the same
manner as first events with the convolution process. Both the addition of the
traveltimes and the multiplication of amplitudes take place. However, there can
be additional features in which cross-convolution artifacts are also generated.
These artifacts which are formed by the convolution of events from different
refractors, occur as relatively steeply dipping events in the refraction convolution
section (RCS) and therefore, they can be removed by dip filtering. The filtered
RCS shows better continuity of events than is the case with the unfiltered
section.
For events which have traveled through the surface layer, the filtered RCS shows
a series of events which occur at a time which is a function of the distance
between the two shot points and the wavespeed in the surface layer. The time of
this event can be used to improve the estimates of the wavespeed in the surface
layer.
158
8.2 - Introduction
The generation of the refraction convolution section (RCS) (Palmer, 2000)
produces a set of traces with the superficial appearance of a seismic reflection
section. It has been demonstrated that the RCS reproduces the time structure of
the refractor interface with first arrivals, while the amplitudes which are largely
corrected for geometric spreading, are essentially a function of the square of the
head coefficient. It has also been demonstrated that the head coefficient is given
approximately by the ratio of the wavespeeds in the upper layer and refractor.
The RCS amplitudes can be employed to image the refractor, to resolve some of
the ambiguities in the determination of wavespeeds in the refractor, and to obtain
a measure of azimuthal anisotropy with three dimensional methods.
To date, research has focused primarily on the portion of the RCS which
corresponds with the first arrivals, and little attention has been directed at later
events. However, the convolution process performs the same operations on later
arrivals as it does with the first events. These operations are the addition of the
traveltimes in the forward and reverse directions, which replaces of moveout from
trace to trace with a constant amount equal to the reciprocal time, the time from
the forward shot point to the reverse shot point, and the multiplication of the
amplitudes. The addition of the traveltimes produces the relative time structure
on the refracting interface, while the amplitude product effectively compensates
for the large geometric spreading which is characteristic of refraction data. The
true time structure on the interface can be obtained by subtracting the reciprocal
time. As the reciprocal time generally decreases with deeper layers, the
shallower layers occur at later times in the RCS.
One feature of the RCS is the generation of what will be termed cross-
convolution events with later arrivals. In this case, the convolution operation
adds arrivals from different refractors, and therefore generates artifacts which
have no geophysical significance. For example, it is possible to produce an
159
event which is the addition of the traveltimes from the refractor in the forward
direction, with the traveltimes from the surface layer in the reverse direction. In
practice, these artifacts occur as pairs, that is there is also an event produced by
the addition of the traveltimes from the surface layer in the forward direction, with
the traveltimes from the refractor in the reverse direction. Because of the
different moveouts or wavespeeds, these events appear as relatively steeply
dipping features in the RCS.
This study describes the use of dip filtering in the f-k domain (Sheriff and Geldart,
1995), to remove the cross-convolution events, with the aim of enhancing those
later events which may have geological significance.
8.3 - Generation of Useful Events and Artifacts in the RCS
The generation of useful later events in the RCS can best be demonstrated with
the ground-coupled air wave. While it is recognized that the imaging of the air
wave has minimal geological significance, it is employed in this study because its
high amplitude improves its clarity in the RCS.
Figure 8.1 is a shot record from a shallow seismic refraction survey at Mt Bulga,
which has been described previously (Palmer, 2001). The record shows the first
arrivals between about 70 ms and 130 ms as very low amplitude signals, and a
very high amplitude event between 350 ms and 1000 ms. The first arrivals are
refracted from the base of the weathering, while the second arrivals are the
ground-coupled air wave.
Figure 8.2 shows the shot record in the reverse direction. The same two events
can be clearly identified, but in this case, the relative amplitude of the ground-
coupled air wave is lower.
160
Figure 8.1: A shot record showing low amplitude first arrivals between about 70
ms and 130 ms refracted from the base of the weathering, and the high
amplitude ground-coupled air wave between 350 ms and 1000 ms.
161
Figure 8.2: Reverse shot record in the reverse direction. The relative amplitude
of the ground-coupled air wave between 350 ms and 1000 ms is lower.
162
Figure 8.3: The RCS generated with the two shots shown in Figures 8.1 and
8.2. The sampling interval has been halved but there has been no subtraction of
a reciprocal time.
163
The RCS generated with these two shots is shown in Figure 8.3. The sampling
interval has been halved (Palmer, 2001), but there has been no subtraction of a
reciprocal time. The presentation gain is low so that the portion of the RCS
which corresponds with the convolved events from the refractor around 150 ms,
becomes essentially featureless. However, the gain facilitates the recognition of
the strong event at approximately 700 ms between stations 29 and 68. The
limited lateral extent of this event occurs because the recording time of one
second was insufficient to record the air wave at the distant detectors.
In addition, the presentation gain highlights the cross-convolution events which
start a few traces from the left side of the section at 300 ms and continue to
about 550 ms near the right side of the section. The recognition of the
companion artifact which starts on the right hand side and finishes on the left is
not as clearly evident in Figure 8.3 and requires more careful inspection.
8.4 - Removal of Cross-convolution Artifacts with Dip Filtering
The transformation of the in RCS in Figure 8.3 from the time-distance domain to
the frequency-wavenumber (fk) domain with the double Fourier transform, is
shown in Figure 8.4. It shows signal centered on the frequency axis, which
corresponds with the horizontal events, and signal spread out parallel to the
wavenumber axis, which is inferred to correspond with the cross-convolution
events.
Figure 8.6 shows the fk domain after the application of a filter to remove all signal
other than that centered on the frequency axis, while Figure 8.5 shows the RCS
after the application of the filter. The event which corresponds with the time-
depth of the ground-coupled air wave can be clearly seen at about 0.710 s.
164
Figure 8.4: The transformation of the in RCS in Figure 8.3 from the time-distancedomain to the frequency-wavenumber (fk) domain with the double Fouriertransform.
Figure 8.6: The transformation of the in RCS in Figure 8.5 from the time-distance domain to the frequency-wavenumber (fk) domain with the doubleFourier transform.
165
Figure 8.5: Refraction convolution section in Figure 8.3 after dip filtering to
remove cross-convolution events.
166
8.5 - Times for Near-surface Events in the Uncorrected RCS
The RCS in Figure 8.5 which has not been corrected by the subtraction of a
reciprocal time, facilitates the computation of wavespeeds for the near surface
layers.
The traveltime in the forward direction t forward, of a seismic signal travelling
through a near-surface layer, that is, for which the depth can be ignored is
tforward = x / V1 (8.1)
Similarly, the traveltime in the reverse direction treverse, at the same detector is
treverse = (d – x) / V1 (8.2)
where x is the forward shot-to-detector distance and d is the separation between
the forward and reverse shot points, and V1 is the wavespeed in the near-surface
layer.
In the RCS in Figure 8.5, these times are firstly summed, then halved, and they
occur at a time tRCS, where:
tRCS = d / 2 V1 (8.3)
It can be readily shown that the ground-coupled airwaves in Figures 1 and 2 has
wavespeeds of about 335 m/s. Using a value of d, the shot point to shot point
distance, of 480 m, the value of tRCS computed with equation is 0.716 s. This
value is similar to that measured above in Figure 8.5.
167
Figure 8.7: Shot record with shot point at station 26.
168
8.6 - Near-surface Wavespeeds from the Uncorrected RCS
Figures 8.7 and 8.8 are two shot records with shot points at the ends of the
geophone spread at stations 26 and 73. The presentation gains are low to
facilitate recognition of a series of events with a wavespeed of approximately 400
m/s, that is, they arrive at the geophones most distant from the shot points after
about 0.600 s. These events occur over an interval of about 0.15 s, and they are
interpreted to be arrivals from the near-surface layer, rather than the ground-
coupled airwave, because of their lower frequency and inferior continuity
compared to the ground-coupled airwaves in Figures 8.1 and 8.2.
A comparison of the unfiltered and filtered RCS in Figures 8.9 and 8.10, shows
that the dip filtering has removed the cross-convolution events, and that the
horizontal and near-horizontal events are emphasized.
Using equation 8.3, it is readily demonstrated that the group of events with the
wavespeeds of 400 m/s should occur at a time tRCS, of 0.3 s. Figure 8.10 shows
a series of events from about 0.33 s to about 0.47 s with higher amplitudes than
the adjacent events. These events are interpreted to represent signals which
have traveled in the surface soil layer. Using the minimum time of 0.33 s and
equation 8.3, the revised wavespeed for this layer is 360 m/s.
It is also possible to recognize a series of events from about 0.18 s with higher
amplitudes than the adjacent events. These events may correspond with arrivals
which travel through the second layer with a wavespeed of approximately 700
m/s.
While the event associated with the ground-coupled air wave convincingly
demonstrates the generation of meaningful later events in the RCS, the
application to events from the near-surface layers is not as clear. It is likely that
further processing, such as deconvolution may be useful. However initial
169
attempts at simple deconvolution methods were not successful, suggesting that
considerably more research may be required to develop suitable techniques.
Figure 8.8: Reverse shot record with shot point at station 73.
170
Figure 8.9: Unfiltered convolution section.
171
Figure 8.10: Dip filtered convolution section.
172
8.7 - Conclusions
Later events, which occur in the shot records, are also treated in the same
manner as first events with the convolution process. Both the addition of the
traveltimes and the multiplication of amplitudes take place. However, there is an
additional feature in which cross-convolution artifacts are also generated. These
artifacts occur as relatively steeply dipping events in the RCS and therefore, they
can be removed by dip filtering. The filtered RCS shows better continuity of
events than is the case with the unfiltered section.
For events which have traveled through the surface layer, the filtered RCS shows
a series of events which occur at a time which is a function of the distance
between the two shot points and the wavespeed in the surface layer. The time of
this event can be used to improve the estimates of the wavespeed in the surface
layer.
8.8 - References
Palmer, D., 2001, Imaging refractors with the convolution section: Geophysics
66, 1582-1589.
Sheriff, R. E., and Geldart, L. P., 1995, Exploration Seismology, 2nd edition:
Cambridge University Press.
173
Chapter 9
Stacking Seismic Refraction Data inthe Convolution Section
9.1 - Summary
The refraction convolution section (RCS) is an effective domain to vertically stack
shallow seismic refraction data, in order to improve signal-to-noise ratios (S/N).
The convolution operation essentially compensates for the effects of geometrical
spreading, and generates traces with much the same S/N ratios. Such traces
are optimum for stacking, unlike the traces on the original shot records.
A major benefit of stacking in the RCS domain is that it takes places before the
measurement of times or amplitudes. With other approaches which do not
routinely employ stacking, such as tomography, any variations in data quality are
addressed with the application of statistical methods to the traveltimes
determined on the original field data.
An essential requirement for stacking in the RCS domain are data which have
been acquired with a continuous roll along approach typical of reflection
methods, rather than with the more common static spread. Such operations are
more efficient and produce more data from the critical near-surface layers, but
they would require significant re-capitalization of most shallow seismic field
operations.
174
9.2 - Introduction
It is well known that the source energy requirements for seismic refraction
surveys are considerably greater than those for seismic reflection surveys for the
same target. The maximum source-to-detector distance for seismic reflection
surveys is generally less than the target depth, whereas the minimum source-to-
detector distance with refraction surveys is usually greater than four times the
target depth. In addition, the geometric spreading component is the reciprocal of
the distance traveled for reflected signals, while the corresponding function for
refracted signals is the reciprocal of the distance squared. Both the longer path
lengths and the more rapid spreading factors result in low refraction amplitudes
and therefore higher source energy requirements. Commonly, the refraction
source is more than ten times the size of the reflection source for the same
target.
Explosives are the standard energy source in most shallow seismic refraction
surveys, and adequate signal-to-noise (S/N) ratios are readily achieved by
increasing the size of the charge. However, this is not always practical in many
environmentally sensitive or urban areas, and it normally results in either poor
quality data due to insufficient charge sizes, or more commonly, no acquisition of
data at all.
In some cases, it is possible to use vertical signal stacking with repetitive
sources, such as hammers and dropping weights. Nevertheless, this approach
can be of limited usefulness, because many repetitions can be required to obtain
reasonable S/N ratios, especially where urbanization is the major source of
noise.
In addition, vertical stacking can result in slow rates of progress where there are
many source points. Walker et al, (1991) demonstrate that one of the most
important factors in improving the reliability of shallow seismic refraction
175
interpretation is a detailed mapping of the wavespeeds in the layers above the
target refractor using a high density of source points. It is not uncommon to
employ a source point between every other pair of geophones.
This study demonstrates the use of vertical stacking with a CMP-like method
using the refraction convolution section (RCS). Redundant or multi-fold
refraction data are acquired with a continuous roll along approach, which is
standard with reflection acquisition. Multiple overlapping RCS are generated with
pairs of shots with the same shot point to shot point separation. The ensemble of
RCS are then sorted and gathered, in much the same way as reflection shot
records, and the gathers are then stacked.
The RCS is a suitable domain in which to stack refraction data when the shot
size, depth and separation are uniform, because the events have approximately
the same S/N ratios. With effective vertical stacking, the S/N ratio improves as
the square root of the number of traces in the stack, but only when the S/N ratios
of the original traces are much the same. Excessively noisy traces, that is traces
with an anomalously low S/N ratio, can significantly degrade the stack and
reduce the benefits of stacking. This situation occurs with stacking refraction
shot records, because there can be large variations in S/N ratios related to the
effects of geometric spreading. Traces at a given station with nearby source
points will have high S/N ratios while traces at the same station with more distant
source points will have lower S/N ratios. The large range in S/N ratios with
refraction shot records significantly reduces the effectiveness of stacking traces
from various shot records with a surface consistent approach.
Shearer (1991) demonstrates stacking shot records in which the shot-receiver
distance is preserved, but not the individual station locations, in order to improve
S/N ratios with earthquake data (see also Lay and Wallace, 1995, p215-216).
However, this approach is not a viable option with shallow refraction data,
176
because it does not accommodate lateral variations in either the depths to or
wavespeeds in the target refractor.
A major advantage of stacking in the RCS domain is that S/N ratios can be
enhanced prior to the measurement of parameters, such as times or amplitudes.
By contrast, tomographic methods measure traveltimes from the shot records
which have varying S/N ratios, and then seek to minimize any errors with a
statistical approach.
9.3 – The Cobar Stacked RCS Section
High resolution data were recorded with a 48 trace recorder, and single 40 Hz
geophones with a 10 m spacing, as part of a regional seismic reflection survey
(Drummond et al, 1992) across the Cobar Basin (Glen et al, 1994), in the central
west of NSW, Australia. The aim of these high resolution lines was to image the
near surface layers, which in these areas were dipping predominantly in the
vertical direction. The seismic source was a 10 kg charge of a high velocity
seismic explosives at a depth of 40 m, and the shot point interval was 30 m.
The data were recorded with off-end shots in both the forward and reverse
directions in order the obtain large shot-to-detector distances for the vertically
dipping reflection targets. For example, the first shot was at station 96, and the
geophones were from station 95 to station 48. The next shot was at station 90
and the geophones from station 89 to 42, a shift of 60 m. Subsequent shots
continued through to station 48 with geophones between stations 47 and 0. The
geophone array then remained static while the shot points at 60 m intervals
within the array at stations 42, 36, etc., were recorded. The recording process
was then reversed. The shot point at station 3 was recorded with geophones
between stations 4 and 51. Subsequent shots were at 60 m intervals (stations 9,
15, etc.), and the geophone array was moved up by the same amount in each
177
case (stations 10 to 57, 16 to 63, etc.). The resulting maximum refraction fold is
six, which is comparatively low.
Figure 9.1: Forward shot record number 65. Shot point is at station 33
178
Figure 9.2: Reverse shot record number 43. Shot point is at station 90.
The first six to twelve traces near the shot point were arrivals from the surface
layer. In order to generate as many useful convolution traces as possible, pairs
179
of shots spaced 63 stations apart were used. This reduced the fold to between
one and four.
Figures 9.1 and 9.2 are two shot records which show the familiar rapid decay of
refracted energy with distance, and in turn, the large variation in S/N ratios with
offset.
Figures 9.3 to 9.8 are a series of RCS over intervals of approximately 30
stations. The structure of the refractor can be readily seen in the RCS. A major
feature of the five RCS is the approximately uniform S/N ratios.
Figure 9.9 is the stacked section, obtained from the five sorted and gathered
sections. While the structure of the refracting interface can be recognized, there
is only a modest improvement in the S/N, due mainly to the low fold of between
180
one and three. Nevertheless, it demonstrates that stacking is efficacious.
Clearly, a much higher fold is necessary to obtain the full benefits of stacking.
181
Figure 9.9: Stacked convolution section.
182
9.4 - The Static Geophone Spread
For the effective use of stacking, it is necessary to use RCS with relatively
uniform S/N ratios, and therefore to employ uniform acquisition parameters. The
relevant parameters are consistent charge size and depth, as well as a uniform
separation between the forward and reverse sources. Although these field
parameters are the norm with the refraction data acquired with routine reflection
surveys for petroleum and coal exploration and for regional reflection surveys in
fold belts, they are uncommon with most shallow seismic refraction surveys
carried out for geotechnical, groundwater and environmental applications.
The majority of shallow seismic refraction surveys carried out for geotechnical
and other shallow applications acquire data in discrete units or spreads. With
these surveys, a static spread of detectors is used with a multiplicity of source
points located at several offset positions on one side, through the spread and to
offset positions on the other side. The number of shot points recorded for each
spread has increased substantially in recent years in order to improve the
determination of the wavespeed stratification above the target refractor, and it is
now common to record more than eleven shots for a spread of 12 detectors
(Walker et al, 1991). For the typical survey length of 400 m for a road cutting,
approximately eight spreads of 12 detectors with a 2 detector overlap are
required (Walker et al, 1991), making a total of 88 shot points.
However, it is questionable whether even this considerable number of shot points
achieves the stated objectives of defining the wavespeed stratification within the
weathered layer. An inspection of published data (Walker et al, 1991, Fig. 6),
shows that the vast majority of the traveltime data (~ 90%) are arrivals refracted
from the base of the weathering. Although it is essential to ensure some
redundancy in the traveltime data in order to resolve the fundamental ambiguity
of determining the number of layers detected (Palmer, 1986, p21-29; Lankston,
1992), the majority of the data from the main refractor are generally not used in
183
the subsequent data processing stages. It questions the fundamental
effectiveness of the static spread approach to acquiring shallow seismic
refraction data.
There are also concerns about the efficiency of field operations with static
spreads. There can be a comparatively large number of shot points per unit
distance because of the occupation of many shot point locations on two and
often three occasions, as well as the common practice of using an overlap of
several detectors. The repeated occupation of shot points can be
environmentally damaging as well as time consuming. Furthermore, field
operations do not progress smoothly, because the acquisition of data ceases
while the spread of detectors is retrieved and then re-deployed for the next
adjacent spread.
The static spread approach also results in a wide range of source energy
requirements for the different offsets and the different layers above the target.
Relatively low source energies are required for signals propagating in the shallow
near-surface layers, while considerably greater source energies are required for
the deeper target refractors. Since the majority of traveltimes are from the main
refractor, there can be a large source energy requirement. As mentioned
previously, many of these times are not used in the data processing, which
suggests that more efficient approaches may be possible.
This study proposes the continuous acquisition of shallow seismic refraction data
be employed routinely for geotechnical, groundwater and environmental
applications.
9.5- Continuous Acquisition of Shallow Seismic Refraction Data
Continuous acquisition of redundant or multi-fold data is the norm with seismic
reflection methods. With this technique, the source point maintains a fixed
184
position in relation to the detector spread, which for land operations is usually a
split spread with the source in the centre of the detectors. The constant
geometry is obtained by laying out more detectors than there are channels in the
recording instruments and then selecting the required channels with a roll along
switch. Continuous and efficient operations are achieved with a single pass of
the seismic source along the line in conjunction with the continual removal of
detectors from the start of the line after they are no longer required, and their
placement at the other end to which the source is progressing.
Better and more uniform coverage of all refractors commonly occurs.
Roll along acquisition methods can provide better data for either the conventional
or convolution approaches where two or more adjacent static spreads would
normally be employed. Comparisons of field operations show that in fact, there
can be a reduction in the number of shot points per unit distance of coverage.
For example, for a 400 m long survey for a road cutting using a 15 m shot
spacing and a 5 m detector interval, a total of only about 55 shots would be
required. A 12 channel seismic recorder would not be suitable for roll along
operations, because the maximum shot to detector distance of 30 m would
generally be insufficient to record enough arrivals from the base of the
weathering. However, a 24 channel seismic recorder, which is widely used in
shallow refraction surveys, would be suitable, and in many cases might even
permit a reduction in the trace spacing to 3m to further improve the resolution of
the wavespeed stratification in the weathered layer. A 48 channel system would
provide further improvements in data quality through additional reductions in
trace spacing, as well as enhanced capabilities with swath or partial three
dimensional profiling.
The comparatively large number of shot points per unit distance with adjacent
static spreads is a result of the occupation of shot point locations on two and
often three occasions, as well as the common practice of using an overlap of
185
several detectors. In contrast, the continuous recording of refraction data with
the roll along method involves the once-only occupation of each source point,
and incorporates a uniform overlap as an integral part of the method.
Accordingly, it represents a more efficient use of equipment and field personnel
with a lower environmental impact.
For surveys in which only limited coverage is required, it is still desirable to
replace the single static spread with a quasi roll along approach. However, the
number of source points may in fact increase with this approach, in order to
obtain sufficient data redundancy for stacking.
9.6 – Determination of Fold with RCS Data
The maximum fold obtained with continuous refraction acquisition using a split
spread shooting method is similar to that obtained with reflection data, viz.
Maximum fold = Number of detectors / (Shot spacing x 2) (9.1)
Note that the shot spacing is given as the number of detector intervals.
The validity of equation 9.1 can be demonstrated with a simple example.
Suppose that the recording system has 48 channels, the shot is at station 25 and
that the live geophones are from stations 1 to 24 and 26 to 49. For the same
split spread recording pattern, reversed shots at stations 1 and 49, which
represent a shot spacing of 24 stations, are the minimum necessary to compute
a time-depth at each detector, provided all arrivals are from the target refractor.
The maximum resulting fold is therefore one.
186
At the other extreme, if the shot spacing is reduced to a single detector interval,
then the maximum fold is 24. For the more common shot spacing of two detector
intervals, the maximum fold is 12.
In general, not all detectors will record arrivals from the target refractor. Suppose
that the first six arrivals on either side of the shot point are from layers other than
the target. For a shot spacing of two detector intervals, the fold will be 9. This
represents a substantial improvement in efficiency over the static spread
approach. The fold or redundancy of nine would be adequate to resolve most
ambiguities in layer recognition, as well as providing moderate improvements in
S/N through stacking. In addition, 25% of the arrivals are from the shallow
surface layers, which represents a significant improvement over the 10% for
typical static spreads, while the overall shot density per unit distance has been
decreased by as much as 40%. Further increases in the proportion of arrivals
from the near surface layers could be achieved by reducing the station spacing.
This analysis has used shot points at the stations themselves, rather than
midway between, which is also common. The benefit of shots at the detectors is
that reciprocal times, the times between the forward shot point and the reverse
shot point can be readily measured in both directions and then averaged. Any
traveltime delays caused by disturbed ground caused by previous shots, can be
avoided by offsetting the shot points by a few metres at right angles to the line.
9.7 - Discussion and Conclusions
The refraction convolution section (RCS) is an effective domain to vertically stack
shallow seismic refraction data, in order to improve signal-to-noise ratios (S/N).
The convolution operation essentially compensates for the effects of geometrical
spreading, and generates traces with much the same S/N ratios. Such traces
are optimum for stacking, unlike the traces on the original shot records.
187
A major benefit of stacking in the RCS domain is that it takes place before the
measurement of times or amplitudes. With other approaches which do not
routinely employ stacking, such as tomography, any variations in data quality are
addressed with the application of statistical methods to the traveltimes
determined on the original field data.
Stacking data in the RCS requires data with uniform acquisition parameters, such
as are typical of seismic reflection surveys for petroleum exploration on land
using split spread and CMP methods.
This would require a major change in field methods with most shallow seismic
refraction operations. In particular, it would require upgrades in acquisition
systems from 12 or 24 channels to at least 48 and preferably 60 channels for 2D
operations, together with roll along and radio shot firing system capabilities.
However, the costs of re-capitalization would be quickly recovered through
improved operational efficiencies.
The use of a roll along acquisition program would result in a reduction in shot
points by up to 40% and in turn a reduction in field time by at least the same
amount. The daily costs for the average three man field crew with 24 trace
equipment are about $A3000. Accordingly, the cost of a new 60 trace field
system at $A90,000 would be equivalent to thirty days of saved field time.
One obvious application of stacking with the RCS is with the computation of
statics, the corrections for variations in the elevations of source and detectors
and for the weathered layer, for regional seismic reflection surveys in fold belts
(Palmer et al., 2000). Accurate weathering corrections are especially important
with regional reflection studies, because continuous reflectors, which are
common in seismic reflection data in sedimentary basins, are very rare in data
recorded across fold belts. As a result, residual statics routines are not effective,
188
and therefore detailed refraction statics analyses are necessary. Frequently, the
arrivals from the base of the weathering can be poor quality and some form of
signal enhancement prior to the measurement of traveltimes might be beneficial.
9.8 - References
Drummond, B, Goleby, B, Wake-Dyster, K, Glen, R and Palmer, D, 1992, New
tectonic model for the Cobar Basin, NSW points to new exploration models for
targets in the Lachlan Fold Belt: BMR Research Newsletter 16, 16-17.
Glen, R.A., Drummond, B.J., Goleby, B.R., Palmer, D. and Wake-Dyster, K.D.,
1994. Structure of the Cobar Basin New South Wales based on seismic
reflection profiling: Australian Journal of Earth Sciences 41, 341-352.
Lankston, R. W., 1990, High-resolution refraction seismic data acquisition and
processing, in Ward, S. H., ed. Geotechnical and environmental geophysics, vol.
1, Investigations in geophysics no. 5: Society of Exploration Geophysicists, 45-
74.
Lay, T., and Wallace, T. C., 1995, Modern global seismology: Academic Press.
Palmer, D., 1986, Refraction seismics - the lateral resolution of structure and
seismic velocity: Geophysical Press.
Palmer, D., Goleby, B., and Drummond, B., 2000, The effects of spatial sampling
on refraction statics: Explor. Geophys., 31, 270-274.
Shearer, P., 1991, Imaging global body wave phases by stacking long-period
seismograms: J. Geophys. Res., 96, 20,353-20,364.
189
Walker, C., Leung, T. M., Win, M. A., and Whiteley, R. J., 1991, Engineering
seismic refraction: an improved field practice and a new interpretation program,
REFRACT: Explor. Geophys., 22, 423-428.
190
Chapter 10
Discussion and Conclusions
10.1 - Shallow Refraction Seismology for the New Millennium: APersonal Perspective
The point of departure for this study was that most current shallow seismic
refraction operations have not taken advantage of advances in technology for
acquisition, processing or interpretation, they are under-capitalized, they are
relatively inefficient, and that the current seismic reflection technology provides
compelling models for the advancement of shallow refraction seismology. Based
on the results of this work, what then are the major features of seismic refraction
operations which might be appropriate to the requirements and the technology of
the new millennium?
A major achievement of this study is a demonstration of the superiority of 3D
results over 2D. There is simply no substitute for the improved quality and
quantity of information which can be obtained from even simple cost-effective 3D
surveys such as that described in this study. It is essential that 3D refraction
methods be adopted as a matter of some priority.
It is likely that the acceptance of 3D shallow refraction methods will parallel the
acceptance of 3D reflection methods by the petroleum exploration and
production industries and the acceptance of high resolution airborne magnetic
and radiometric data by the mineral exploration industries. Initially, cost was
191
considered to be the major reason for the relatively low levels of acceptance of
these methods. However, this situation changed rapidly when it was widely
demonstrated that high spatial sampling densities in all directions, is one, if not
the most important factor, in reducing risk through improved geological
interpretations. This conclusion is supported by the 3D results described in this
study.
The development of a 3D oriented approach implies the use of specialist seismic
contractors for acquisition in order to employ field systems with greatly increased
capabilities, as well as to promote efficient field operations. It is difficult to justify
the use of relatively expensive professional expertise to carry out routine
unskilled field duties with under-capitalized systems and inefficient operations.
Increasing channel capacity to at least 150 and doubling the number of shot
points could achieve efficient 3D field operations. This would result in an
increase of at least an order of magnitude in the amount of data, and in turn it
would dictate the use of efficient methods of data processing and interpretation.
Full trace processing with the RCS is a simple and efficient approach for
processing any volume of seismic refraction data.
It is likely that the increased quantity and quality of data obtained with 3D surveys
might stimulate a change in the roles of the geophysicist from acquisition and
processing towards interpretation. It also implies inclusion of other geoscientists
at earlier stages of the interpretation process, in order to generate more complex
and more geologically meaningful interpretation models.
The format of data processed with the RCS facilitates the convenient application
of current reflection processing and interpretation technology to shallow seismic
refraction data. Although the existing software developed specifically for
refraction seismology represents many man-years of effort, it is relatively
insignificant when compared with the software developed for reflection
192
seismology. Just as the use of imaging processing software, which was
developed originally for remotely sensed data, has increased the detail of the
geological interpretation of magnetic, radiometric and gravity data, so seismic
reflection software is a vast resource which has the potential to extract even
greater information from refraction data. In particular, the data processed with
the RCS is suitable for analysis with software used for the interpretation of
processed seismic reflection data. Such software includes basic functions for
picking times and amplitudes of horizons, as well as post-processing functions,
such as attribute analysis. Attribute processing of RCS data may have as large
an impact on increasing the detail of the interpreted geological model as it has
with reflection data.
The author’s preference for an approach which is essentially an extension of the
GRM, is hardly surprising. However, other approaches, such as tomography are
currently not viable alternatives. The major shortcoming of tomography is that
the large increase in the number of shot points, commonly by at least an order of
magnitude over a simple GRM approach suggested here, would result in high
and possibly prohibitive costs of acquisition. Furthermore, tomography has yet to
satisfactorily address either the issues of non-uniqueness, large variations in
wavespeeds in the refractor, or anisotropy.
The RCS offers a new approach to generating more complex geological models
from shallow seismic refraction data through the use of the complete seismic
refraction trace and therefore, the use of amplitudes as well as traveltimes. In
time, it may stimulate the development of routine refraction methods which are
comparable in sophistication to current 3D reflection methods.
193
10.2 - Conclusions
This study demonstrates that the refraction convolution section (RCS), generated
by the convolution of forward and reverse shot records, is an efficacious
approach to full trace processing of shallow seismic refraction data.
The convolution operation effectively adds the first arrival traveltimes of each pair
of forward and reverse traces. Like the many standard methods for processing
refraction data which use addition to obtain a measure of the depth to the
refracting interface in units of time, the RCS also produces a similar time image
of the refractor. In this study, the reciprocal time, the traveltime from the forward
shot point to the reverse shot point, is subtracted and the result is then halved
(by halving the sample interval of the trace headers) to form the equivalent of the
time-depth function of the generalized reciprocal method (GRM).
The generation of the RCS requires no estimates of, or assumptions about the
wavespeeds in either the refractor or the overlying layer. Any lateral changes in
refractor wavespeeds are accommodated through the use of forward and reverse
data.
The convolution operation also multiplies the amplitudes of first arrival signals.
This operation compensates for the large effects of geometric spreading to a very
good first approximation, with the result that the convolved amplitude is
essentially proportional to the square of the head coefficient. The signal-to-noise
(S/N) ratio of the RCS shows much less variation than those on the original shot
records.
A significant achievement of this study is the demonstration that the head
coefficient is approximately proportional to the ratio of the specific acoustic
impedances in the upper layer and in the refractor, under the conditions
encountered in most shallow seismic refraction surveys. These conditions are
194
that there is a reasonable contrast between the specific acoustic impedances in
the layers. Although the original theoretical formulations of the head coefficient
were published almost half a century ago, the very limited use of refraction
amplitudes since that time has not resulted in convenient approximations which
might facilitate practical quantitative analysis with routine surveys. It is likely that
the simplification proposed in this work will promote greater use of amplitudes in
routine shallow seismic refraction surveys.
A major part of this study has been the demonstration of the usefulness of either
the convolved amplitudes or the equivalent shot amplitude products. The two 2D
case histories at Mt Bulga demonstrate both the correlation between amplitudes
and wavespeeds, and the use of amplitudes in addressing any ambiguities in the
determination of wavespeeds.
Non-uniqueness in determining wavespeeds in the refractor is an important
issue. Although most geophysicists tacitly accept that the inversion of seismic
refraction data need not necessarily produces a unique solution, the results of
most inversion routines still do not adequately reflect this reality. The non-
uniqueness can occur in the determinations of wavespeeds in both the upper
layer and the refractor and often, they are inter-related. This study proposes
several solutions to non-uniqueness in the refractor wavespeeds. Firstly, the
GRM can be used to generate a family of acceptable starting models for model-
based inversion or tomography. Secondly, the minimum variance criterion of the
GRM can be employed to determine a most likely starting model. Finally,
amplitudes can provide additional valuable information to constrain any starting
models.
The RCS can also include a separation between each pair of forward and
reverse traces in order to accommodate the offset distance in a manner similar to
the XY spacing of the GRM and to improve lateral resolution. The offset distance
is the horizontal separation between the point of refraction on the interface and
195
the point of detection at the surface. Although the differences between the updip
and downdip offset distances can be large, their sum which is obtained with the
optimum XY value, is relatively insensitive to the dip angle. It facilitates the
application of the RCS to deeper refractors where the offset distances are
significant, as well as to very shallow refractors.
The refraction profile across the Mt Bulga massive sulfide orebody demonstrates
that there can be lateral separations of the amplitudes on the forward and
reverse shot records. In this case it was the distinctive low amplitudes
associated with the mineralization. This separation is similar to the optimum XY
value determined from the refractor wavespeed analysis function. In addition,
this case history demonstrates that the accommodation of the offset distance
with finite XY values is efficacious for improving lateral resolution with shallow
refractors and with detector separations as small as 2.5 m.
Another important achievement of this study is the examination of the effects of
variations in the near-surface soil layers on amplitudes or “amplitude statics”.
The profile across the Mt Bulga orebody demonstrates that the increases in the
thickness of the surface soil layer correlate with increases in refraction
amplitudes, and that these increases are adequately described with the
transmission coefficients of the Zoeppritz equations. Where these surface layers
are laterally continuous, the same increases in amplitudes occur at each
detector, and therefore the relative amplitudes are preserved. However, where
the surface layers are laterally discontinuous, the amplitudes can be quite
variable. If the wavespeeds in these zones can be measured, then corrections
can be applied with the Zoeppritz equations. Where this is not possible, then the
minimum amplitudes, rather than an average should be used.
Perhaps the most exciting aspects of this study are the results of the 3D survey.
Even with the nominally 2D structure of the shear zone, there are important
lateral variations in both refractor depth and wavespeed, which could not be
196
predicted on the basis of the earlier 2D survey. In addition there are important
variations in the direction of the rock fabric as inferred from the qualitative
measures of azimuthal anisotropy. These results are a compelling
demonstration that more useful geological interpretations are possible with
simple 3D sets of data with complete spatial coverage in all directions, rather
than with the most detailed inversion of 2D sets of data.
Just as 3D reflection methods have revolutionized petroleum exploration and
production, so it is anticipated that shallow 3D seismic refraction methods will
eventually be recognized as a cost-effective approach to minimizing risk,
especially with geotechnical and environmental investigations. The results of the
3D survey raise the question of whether the 2D model of the subsurface is a
satisfactory approximation for most seismic refraction targets.
Another significant advantage of the use of 3D amplitudes, is that they provide a
measure of refractor wavespeeds at each detector, whereas the analysis of
traveltimes provides a measure over several detectors, commonly a minimum of
six. Therefore, amplitudes effectively improve the spatial resolution by almost an
order of magnitude. It is likely that amplitudes will facilitate the extraction of even
more detail with, for example, the attribute processing methods currently being
used with the interpretation of 3D seismic reflection data.
The RCS provides another approach to the use of later events. “Cross-
convolution” artifacts can be easily removed with simple dip filtering methods,
thereby highlighting those events from other, generally shallower, layers. It is
likely that the application of standard seismic reflection processing steps such as
dip filtering, deconvolution and migration or imaging, will result in the extraction of
further information from the RCS.
197
The RCS also provides an effective approach to the high source energy
requirements of refraction seismology through stacking in a manner similar to the
CMP methods of reflection seismology.
The RCS is a simple and efficient method for full trace processing of shallow
seismic refraction 2D and 3D data. The convolution process is very quick and
not particularly demanding of computing facilities.
The RCS can be viewed as a simple extension of the GRM. It facilitates
improved interpretation of shallow refraction seismic data through the convenient
use of amplitudes as well as traveltimes.
198
Appendix 1
Comments on “A brief study of thegeneralized reciprocal method and
some of the limitations of themethod” by Bengt Sjögren.
A.1 - Introduction
Despite the implications of the title, Sjögren (2000) is essentially a rejection of the
generalized reciprocal method (GRM). It is clear that there are very few aspects
of the GRM, which Sjögren finds acceptable, if in fact there are any at all.
I consider that most of the substance of his critique of the GRM is either wrong or
ill informed, while other aspects are simply matters of opinion. A thorough
response to his paper would be quite lengthy and somewhat technical, and so I
will restrict my response to three main issues only. They are the following:
1. Do we always need to define all layers above the target refractor, or are there
situations where the use of an average wavespeed is more appropriate? I
accept that Sjogren’s re-interpreted depth sections for the two case histories
provide a better estimate of total depths across the complete profiles. However,
the original depth sections in Palmer (1991) generated with the average
wavespeeds are appropriate to the objectives of each case history and to those
of the paper as discussed in the third issue.
199
2. There are no fundamental mathematical differences between the GRM and
the collection of methods used by Sjögren. From both an assessment of his
descriptions of the various methods and an examination of some of his figures, it
is clear that both approaches generate essentially the same processed data, and
that his mean-minus-T method is identical to the wavespeed analysis function of
the GRM. The main differences are found in how the processed data are
interpreted. The GRM provides a systematic and objective framework, whereas
in my opinion, Sjögren's approach is neither consistent nor objective.
3. Sjögren has not addressed the fundamental aims of Palmer (1991), namely
the demonstration of an objective approach for recognizing and defining narrow
zones with low wavespeeds in the refractor. Sjogren’s approach is not
systematic and it frequently relies on personal judgment, and as a consequence,
it can result in the generation of artifacts, such as those shown in Sjögren (2000,
Fig. 5(a)).
A.2 - The Use of Average Wavespeeds
Sjogren’s assertion of the need to define all layers in geotechnical investigations
rather than use an average wavespeed, is not cognizant of the objectives of each
field study in Palmer (1991). In the case of the collapsed doline, the objective
was to determine its depth. Drilling had not been completely successful because
the loss of circulation in the rubble had stopped progress at depths of about 50 m
and before solid rock was encountered. Therefore, the issue is not whether
Sjogren’s detailed layer by layer approach is more useful than the average
wavespeed approach, but rather why both approaches, which give much the
same maximum depths of about 15 m, are clearly at variance with the drilling.
Irrespective of which interpretation approach is used, it is obvious that the
original objectives have not been achieved, that the refraction method is
200
inappropriate for solving the problem, and that the differences in depth
computations between stations 56 and 71 are peripheral to the survey objectives.
Furthermore, Sjogren’s reasoning for his rejection of my explanation that the
refracted energy propagates around rather than under the doline is convoluted,
not convincing and ignores the implications of a genuine three-dimensional
structure.
The second case history was across a fault. An earlier high-resolution reflection
survey had been carried out in order to test whether the method was efficacious
in detecting known faulting in the underlying coal seams. These results were
poor, possible due to the proximity of the line to a busy road and the use of small
explosive charges. The refraction survey was then carried out in order to
generate a more accurate set of statics corrections. The component of the
statics corrections, which effectively replaces the weathered layer with
unweathered material prior to the application of the elevation component, was
generated with the approach described in Dobrin (1976, p.215), Palmer (1995),
and Palmer et al (2000). This method simply scales the time-depths by a factor
which is a function of the average wavespeed in the weathering and the
wavespeed in the refractor.
Therefore, while I accept that Sjogren’s depth computations may be more
appropriate to many types of geotechnical investigations, the use of average
wavespeeds in Palmer (1991) was entirely compatible with the aims of both the
field studies and the paper. Furthermore, the differences in total depths between
Palmer (1991) and Sjogren (2000) do not alter the major conclusions of the
paper with respect to the study of narrow zones with low wavespeeds.
In Palmer (1981, 1992, 2000a, 2000b) I demonstrate the use of the average
wavespeed in accommodating undetected layers, wavespeed reversals and
transverse isotropy. At present, there are no other published approaches to
solving these problems which are commonly encountered in many parts of the
201
world, especially those with deep regoliths. Furthermore, Sjogren appears to
have overlooked other case studies (Palmer, 1980) in which the GRM has
effectively defined all layers above the target refractor.
A.3 - The Similarities Between The GRM and Sjogren’s Approach
Sjogren's indication that he does not accept the usefulness of this average
wavespeed is surprising, especially since he also uses an average overburden
wavespeed with Hales' method where there are multiple layers (Sjogren, 2000,
p.819). In that same paragraph, he also describes varying the XY distances, in a
manner analogous to the GRM. It is clear that there are many similarities
between the two average wavespeeds, but that the GRM has extended the
concept to include several important benefits such as the ability to accommodate
undetected layers, as well as minor differences in the interpretation of the
traveltime graphs. Even though Sjogren has emphasized such differences, there
are only minor differences in depth computations between Palmer (1991) and
Sjogren (2000) at the points where the average wavespeeds were determined.
While, it is acknowledged that the accuracy of the average wavespeed can be
reduced with distance from the point of computation, this is not necessarily a
problem, as was the case with the two field studies.
Furthermore, the similarities extend beyond the average wavespeeds. In Palmer
(1986), I describe Hales’ method and conclude that fundamentally, it is very
similar to the GRM. The similarities are that both methods obtain a measure of
the depth to the refracting interface in units of time through the addition of
forward and reverse traveltimes, and a measure of the wavespeeds in the
refractor from the differencing of the same forward and reverse traveltimes. The
equations describing these two operations for each method are virtually identical.
Furthermore, both methods employ refraction migration in order to accommodate
the offset distance, which is the horizontal distance between the point of
202
refraction on the interface and the point of detection at the surface. The major
difference between the two methods is that Hales’ method achieves the addition,
subtraction and migration with a graphical approach, while the GRM performs
these same operations arithmetically with the scalar traveltimes prior to their
graphical presentation.
Sjogren’s determination of refractor wavespeed in the main refractor begins with
a version of the ABC method, which is a special case of the GRM with a zero XY
spacing. In the first case history but not the second, he then applies what is
clearly another version of the GRM wavespeed analysis algorithms with finite XY
values. Sjogren (2000, p.825) refers to curves 1 and 2 in his Fig. 4 as having
been computed with migrations of 5 m and 7.5 m. These curves bear a
remarkable resemblance to those computed with the GRM with similar XY values
in Palmer (1991, Fig.16). The minor differences are due to Sjogren’s editing of
the traveltime data. He then applies yet another method, namely Hales’ method,
to further refine his wavespeed determinations.
Accordingly, Sjogren’s approach with a number of methods and my approach
with the GRM are essentially generating the same set of computations for the
determination of refractor wavespeeds. Where Sjogren uses a succession of
different techniques, all of which employ addition and subtraction of forward and
reverse traveltimes, together with accommodation of the offset distance with
migration, the GRM achieves the same results within a single presentation, such
as in Palmer (1991, Fig.16).
Are the differences between Sjogren’s approach and the GRM important? An
examination of the wavespeed analysis function in Palmer (1991, Fig.16), reveals
that the graph for an XY value of 5 m has intervals of steeper gradient which
would correspond with lower wavespeeds, and which occur over the same
intervals with low wavespeeds shown in Sjogren (2000, Fig 5(a)). Therefore,
where Sjogren interprets virtually every change in slope in a single graph as
203
evidence of a lateral change in wavespeed, I have concluded that there is
insufficient evidence for the existence of those intervals. This conclusion is
based on an assessment of all of the graphs in Palmer (1991, Fig. 16), and the
recognition of systematic changes in patterns, rather than the detailed
interpretation of a single graph.
The issue then, is whether the approach of Palmer (1991) under-interprets the
data and therefore overlooks intervals with low wavespeeds, or whether Sjogren
(2000) over-interprets the data and generates artifacts which do not really exist.
The issue is important because such features can be very significant in most
geotechnical, groundwater and environmental applications.
A.4 - Recognizing And Defining Narrow Zones With LowWavespeeds In Refractors
This difference between the two approaches introduces the fundamental
question which Palmer (1991) seeks to address. Is there an objective and
systematic approach, which is independent of individual interpretation styles, for
recognizing and defining narrow zones with low wavespeeds in refractors?
At the present time, I am still largely of the opinion that most narrow zones with
low wavespeeds are simply artifacts of the inversion algorithms and individual
interpretation styles. I have several reasons for holding this view.
1. Numerous model studies have shown that the algorithms which seek to
determine the wavespeeds in the refractor through the differencing of forward
and reverse traveltimes, can readily produce narrow zones with alternating high
and low wavespeeds, where there are significant changes in depths to the
refractor. This pattern can be seen for example, in the vicinity of the doline in
204
Sjogren (2000, Fig. 5(a)), and it immediately raises doubts about the veracity of
the lateral changes in wavespeed.
2. The traveltime differences through these zones are small, and frequently they
are within the acceptable errors of the data, commonly plus or minus one
millisecond. For example, in Palmer (1991, Fig. 16), the differences between the
computed points and the line representing the fitted wavespeed function are
generally less than a millisecond for the optimum XY value of 5m. (Also see Fig.
2 to be discussed later.) I question whether any minor changes in slope are
statistically significant.
3. The issue is not resolved with forward modeling with either ray tracing or the
eikonal equation which are employed, for example, with tomographic and other
model-based methods of inversion. The GRM is able to generate a family of
geologically acceptable starting models in which the wavespeeds range from low
to high in narrow zones (Palmer, 2000c; 2000d) and which essentially satisfy the
original traveltime data (Palmer, 1980, p.49-52; 1986, p.106-107) to better than a
millisecond. This is simply another statement of the fundamental problem of non-
uniqueness with all inversion processes (Oldenburg, 1984; Treitel and Lines,
1988), but it is rarely if ever, addressed satisfactorily with refraction methods.
Therefore all of the refractor wavespeed models generated with different XY
spacings in Palmer (1991, Fig. 16), satisfy the traveltime data. Although some of
these can be rejected on simple geological grounds, such as those with negative
wavespeeds which obviously are not geologically realizable, there still remains a
range of models which fit the data to an acceptable accuracy.
Therefore, while I accept that Sjogren has a methodology for determining
wavespeeds in his Figures 1, 2, 4 and 7, I do not accept that he has addressed
the fundamental issues of non-uniqueness, that is of recognizing the artifact from
the real. I consider his approach relies heavily on his familiarity with Hales’
method and modified versions of the ABC method which include a migration
205
process identical to the GRM, and that there is a significant subjective
component. Furthermore, his approach is poorly explained and the use of
different interpretation methods in a non-systematic manner is confusing.
Accordingly, I do not accept that he has demonstrated he has an objective
approach, or that there are narrow zones with low wavespeeds in the refractor in
Figure 5(a).
Sjogren's non-systematic approach can also be demonstrated with the model
data. Sjogren (2000, Fig. 2) changes the inclinations of the slope lines in the
Hales time loop on the basis of the wavespeeds derived from curve 1 for which
the XY value is zero, in order to improve the resolution of the narrow zone with
the low wavespeed. If that approach were to be employed with the first model
(Palmer, 1991, Fig. 2), using the wavespeed analyses for zero XY shown in
Palmer (1991, Fig. 5), then artifacts with both high and low wavespeeds at the
sloping interface, would be generated.
A.5 - Use Of Alternative Presentations And Amplitudes ForDetermining Wavespeeds In Refractors
In Palmer (1991), I present a systematic and objective criterion, generally known
as minimum variance. It is clear that an important aspect of this approach is to
determine a gross model of the refractor wavespeeds, and then to systematically
fit this model as has been done in Palmer (1991, Fig. 16). This is usually an
iterative process, simply because it is difficult to obtain the correct wavespeeds at
the first attempt. It can be somewhat challenging because of the need to
recognize the pattern of the departure of the computed points from the fitted line
as shown in Palmer (1991, Fig. 5), while at the same time accommodating the
normal errors in field data.
206
Recently, I have been experimenting with averaging the wavespeed analysis
function for a range of XY values which range from less than to greater than the
optimum value, with the range of XY values being symmetrical about the
optimum, in order to derive a gross wavespeed model. This process minimizes
many of the apparent changes in wavespeeds due to structure where there are
no narrow lateral changes in wavespeed such as is shown in Palmer (1991, Fig.
5), and it averages many of the errors in picking traveltimes.
Figure 1. Refractor wavespeed analysis function, averaged for XY values from
zero to 10 m.
Using the traveltime data for the doline field study, Figure 1 shows a graph in
which the points computed with the wavespeed analysis function for XY values
207
between zero and 10 m have been averaged. While Figure 1 risks introducing
yet another model of the refractor wavespeeds, namely about 2150 m/s between
stations 39 and 58, and about 3240 m/s elsewhere, there is no indication of any
narrow zones with low wavespeeds.
Figure 2. Differences between the averaged refractor wavespeed analysis
function in Figure 1, and the individual refractor wavespeed analysis functions for
XY values from zero to 10 m.
208
Figure 2 shows the differences between the average and the computed
wavespeed analysis function for XY values between zero and 10 m. The
patterns of these differences are also consistent with there being no narrow
zones with low wavespeeds. The minimum differences occur for an XY value of
5 m, they are essentially random, and they are generally less than a millisecond.
Figure 3. Amplitudes of the forward and reverse offset shots.
I have also been investigating the use of amplitudes as an additional approach to
addressing this fundamental problem of non-uniqueness. In Palmer (2000e,
209
2000f), I demonstrate that the amplitude of the refracted head wave, after
correction for geometrical spreading and refractor dip with either convolution or
multiplication, is essentially a function of the head coefficient. I further
demonstrate that the head coefficient is approximately proportional to the ratio of
the specific acoustic impedances (which is the product of the wavespeed and
density) in the overburden to that in the refractor. Therefore, arrivals from zones
in the refractor with low wavespeeds should exhibit high amplitudes, and vice
versa.
The amplitudes for the shot records are shown in Figure 3. For the shot at
station 1, the amplitudes shown a strong decay between stations 24 and 40,
which is interpreted as the geometric effect, together with an interval with
extremely low amplitudes between stations 46 and 59. The amplitudes for the
reverse shot at station 97 show a much less pronounced geometric effect, but
again there is an interval with very low amplitudes between stations 38 and 44.
These very low amplitudes were a major limitation on the measurement of
accurate and consistent traveltimes.
There are very few model studies on the effects of structure on the refraction
amplitudes. Nevertheless, it is unlikely that the very low amplitudes are
compatible with the simple refraction of energy from under the survey profile, but
rather with some form of scattering. They are compatible with energy
propagating around the doline and being scattered back to the surface through
the highly attenuating medium of the rubble in the collapsed doline, as has been
proposed in Palmer (1991).
Figure 4 shows the product of the amplitudes computed with an XY separation of
5 m. The results can be broadly separated into three regions which correspond
approximately with those recognized in Figure 1 with the wavespeed analysis
function. They include high values between stations 24 and 37, very low values
between stations 37 and 59, and higher values between stations 59 and 71. The
210
lower values between stations 37 and 59 suggest higher rather than lower
wavespeeds. The slightly higher values between stations 37 and 42 suggest
lower wavespeeds, but do not correlate with those determined by Sjogren (2000,
Fig5(a)). However, these results should be used with considerable caution
because of the severe attenuation of seismic energy within the rubble of the
collapsed structure, and because the feature is three-dimensional.
Figure 4. Product of the forward and reverse shot amplitudes shown in Figure 3
with an XY value of 5m in arbitrary units.
Therefore, it seems probable that there are no narrow zones with low
wavespeeds associated with the collapsed doline and that Sjogren has
generated artifacts through an over-interpretation of the processed data.
211
A.6 - A Systematic Approach With The GRM
In summary, Sjogren concludes that “the GRM must be regarded as being of
limited use for detailed and accurate interpretations of refraction seismics for
engineering purposes.” This conclusion is surprising given the very close
similarities between the GRM, Hales’ and the mean-minus-T methods, the latter
two of which he clearly favours. The fact that key features of his figures bear a
remarkable resemblance to parts of the GRM presentations emphasizes the
essential similarities between the two approaches. Instead, he has sought
differences where none really exist, and as a corollary, he has not recognized
similarities where in reality there are many. On the basis of the fundamental
similarities between the GRM and Sjogren's processed data, I conclude that
Sjogren (2000) is more a demonstration of his interpretation style and experience
using a collection of methods rather than a cohesive assessment of the GRM.
However, his approach is neither systematic nor entirely objective, and as a
result it is prone to the generation of artifacts.
It seems that the aim of his paper is to emphasize minor differences, mainly in
the assignment of layers to the traveltime graphs, and then to illogically imply
fundamental shortcomings of the GRM. As such, his paper lacks balance and
objectivity, and it is more in keeping with seeking a conviction in an adversarial
court system than with a scientific journal.
Sjogren (2000) has produced nothing of substance which requires any
fundamental re-assessment of the main features of the GRM in general, or of
Palmer (1991), in particular. He has not satisfactorily addressed the aims of
Palmer (1991), namely an objective method for the recognition and definition of
narrow zones with low wavespeeds. His conclusion that the GRM is unsuitable
for geotechnical applications is not substantiated, and it is based on an
212
incomplete understanding of the GRM and of Palmer (1991), rather than any
genuine shortcomings with the method.
A.7 - The Need To Promote Innovation In Shallow RefractionSeismology
Sjogren’s rejection of the GRM as a useful method for processing and
interpreting shallow seismic refraction data, does little to encourage others to
present new approaches through fear of biased criticism. It is irresponsible and
does not advance the science through balanced and objective debate. In the last
fifty years, innovation in shallow refraction seismology has been rather modest at
best and it has focused predominantly on the various competing methods for
inverting field data. Until there is widespread consensus through a recognition of
fundamental similarities between these inversion methods, there will be little
advancement in other equally important aspects of the science. By comparison,
reflection seismology has achieved major advances through the development of
common midpoint methods, digital signal processing, three-dimensional methods
and sophisticated computer interpretation programs, over the same period of
time. It is now time to move on to the refraction techniques which will be
appropriate to the requirements and the technology of the new millennium.
In Palmer (2000g), I demonstrate the generation of the refraction convolution
section (RCS) through the convolution of forward and reverse traces. The
addition of the traveltimes with convolution is equivalent to that achieved
graphically with Hales’ method and arithmetically with the GRM. The RCS
facilitates full trace processing of seismic refraction data and in turn, the
examination of many important issues such as signal-to-noise ratios, “amplitude
statics”, 3D refraction methods and azimuthal anisotropy, signal processing to
enhance second and later events and stacking data in a manner similar to CMP
reflection methods. A major advantage of the RCS is that it incorporates
213
amplitudes and time structure within a single presentation, facilitating the
resolution of many of the non-uniqueness issues discussed here. It is extremely
rapid and suitable for use with any volume of data, and therefore it can be readily
included in the processing of refraction data with any method.
The RCS is a new approach to obtaining more and better information from
shallow seismic refraction data, and in time it may even supercede the GRM as
well as Hales’ method. However, that possibility can only occur if there is a shift
in culture from one of conflict and an emphasis on minor differences to one of
consensus and an emphasis on fundamental similarities which traditionally, has
characterized the scientific method.
Sjogren’s critique of the GRM does not seek the consensus essential for the
advancement of the science of shallow refraction seismology. Regrettably, it is
neither balanced nor objective, it shows minimal insight into the fundamental
similarities of most methods of refraction inversion, it does not provide an
alternative systematic approach to refraction interpretation, it does not address
the important issues of non-uniqueness, and it does not provide a vision for
future innovation.
A.8 - References
Dobrin, M. B., 1976, Introduction to geophysical prospecting, 3rd edition:
McGraw-Hill Inc.
Oldenburg, D. W., 1984, An introduction to linear inverse theory: Transactions
IEEE Geoscience and Remote Sensing, GE-22(6), 666.
Palmer, D., 1980, The generalized reciprocal method of seismic refraction
interpretation: Society of Exploration Geophysicists.
214
Palmer, D., 1981, An introduction to the generalized reciprocal method of seismic