Vibratory signal sweep seismic prospecting method and apparatus

Abstract

Disclosed are a method and apparatus for vibratory seismic prospecting in which a vibratory signal sweep is produced by a vibrating apparatus in accordance with a power spectrum function S(f), which is proportional to a noise component function n(f) and a subsurface signal amplitude attenuation function B(f) for an area to be surveyed, to thereby produce reflection signals having an improved signal-to-noise ratio.

Description

BACKGROUND OF THE INVENTION

1. Field of the Invention

The present invention relates to a seismic prospecting method and apparatus which use an improved vibratory signal sweep, and more particularly to a seismic prospecting method and apparatus which use a vibratory signal sweep which is formulated in consideration of attenuation and signal scattering characteristics of the particular earth formation target with which the vibratory signal sweep is used.

2. Description of the Prior Art

The use of vibratory apparatus in seismic prospecting is well known in the art. Commonly, a vibratory sweep is used in such apparatus to vibrate the earth. The sweep typically lasts for 20 to 30 seconds, during which time the instantaneous frequency of the vibratory oscillating signal varies linearly and monotonically with time, usually from a first lower frequency f.sub.1 to a second higher frequency f.sub.2. The amplitude of the oscillation signal remains substantially constant over the duration of the sweep, but it is preferably linearly tapered near the beginning and ending of the sweep to avoid signal overshoots and to facilitate signal processing of the sweep wave reflected from subsurface formations.

While this type of vibratory sweep is considered good for general purpose seismic prospecting, it has limitations since the various frequency components of the sweep are affected differently by the signal scattering and absorption effects of the earth formation. Commonly, when the vibratory sweep signal is received by a receiver after being reflected by subsurface formation conditions, the higher frequency components thereof are scattered and attenuated to a greater degree than the low frequency components. For data processing purposes, the received reflected signal waveform should have as flat an amplitude characteristic as possible. However, processing the received signal to yield a flat amplitude spectrum by amplifying (equalizing) the higher frequency components of the reflected wave also undesirably increases a background noise component at the higher frequencies, so that the noise component increases in amplitude, with increasing frequency of the vibratory sweep.

Many variations to the linear vibratory sweep discussed above have been proposed to control the amplitude of the vibratory sweep signal throughout the frequency spectrum of interest. See, for example, the paper "Signal Design In The `Vibroseis`.RTM. Technique" by Goupillaud, published in Geophysics, Vol. 41, No. 6, December 1976, pages 1291-1304. However, these variations are based on changing the vibratory signal sweep in accordance with a predetermined mathematical function without regard to the specific characteristics of the subsurface formation with which the vibratory apparatus is used. Accordingly, although one or more of the known vibratory sweep signal patterns may work well with one type of subsurface formation, they will not necessarily work well with other formations having different signal scattering and attenuation characteristics.

SUMMARY OF THE INVENTION

Accordingly, one object of the invention is the provision of a method and apparatus for seisimic prospecting which employ an improved vibratory signal sweep pattern formulated in accordance with subsurface signal attenuation and noise characteristics of an area to be surveyed, so that the sweep pattern is optimized for a particular area being explored and the frequency spectrum of the reflected wave can be amplitude equalized, with an improved signal-to-noise ratio, thus facilitating subsequent processing of the reflected wave.

These objects are obtained in the method of the invention which comprises the steps of:

determining a subsurface amplitude attenuation function B(f) for a frequency spectrum between frequency components f.sub.1 and f.sub.2 for a subsurface area to be surveyed; defining a noise component function n(f) associated with the area for the frequency spectrum; producing a vibratory sweep signal over the frequency spectrum which is a function S(f) of the subsurface amplitude attenuation function B(f) and the noise component function n(f); providing the vibratory sweep signal to an earth vibratory apparatus; and vibrating the earth with the apparatus.

These objects are also obtained with an apparatus of the invention which comprises:

means for generating a vibratory sweep signal over a frequency spectrum existing between a first frequency f.sub.1 and a second frequency f.sub.2, the vibratory sweep signal having a frequency domain power spectrum function S(f) which is a function of a frequency domain amplitude attenuation function B(f) and a frequency domain noise component function n(f) over the frequency spectrum for a subsurface area to be logged; and means responsive to the signal generating means for vibrating an earth formation with said vibratory sweep signal.

The above objects, features and advantages of the invention, as well as others, will be more clearly discerned from the ensuing detailed description of the invention which is provided in connection with the accompanying drawings.

BRIEF DESCRIPTION OF THE DRAWINGS

FIG. 1 shows a vibratory seismic prospecting apparatus employing the method of the invention;

FIG. 2 shows an amplitude spectra of a conventional vibratory signal sweep;

FIG. 3 shows an amplitude spectra of recorded reflection signals and background noise typically associated with the conventional vibratory signal sweep;

FIG. 4 shows the FIG. 3 amplitude spectra after signal correlation;

FIG. 5 shows the FIG. 4 amplitude spectra after amplitude spectrum flattening;

FIG. 6 shows a power spectra after the amplitude spectra of FIG. 4 is flattened;

FIG. 7 shows a vibratory signal sweep power spectra produced in accordance with the method of the invention; FIG. 8 shows a sweep time function for producing the vibratory signal sweep power spectra shown in FIG. 7;

FIG. 9 shows in block diagram form a frequency domain signal generator for producing the FIG. 6 power spectra; and

FIG. 10 shows in block diagram form a time domain signal generator for producing the FIG. 6 power spectra.

DETAILED DESCRIPTION OF THE INVENTION

To fully explain the invention, a brief description of the operation of a typical vibratory signal generating system will be first provided in connection with FIG. 1. In the ensuing description, functions denoted by (f) are represented in the frequency domain, while functions denoted by (t) are in the time domain.

In vibratory sweep seismic prospecting a vibratory apparatus 13 is provided as an acoustic sound source, which generates vibratory acoustic waves into an earth formation. The acoustic waves are typically generated over a time perid of 20-30 seconds, during which the frequency of the wave is varied from a lower frequency f.sub.1 to a higher frequency f.sub.2 with a constant amplitude, as shown in FIG. 2. The amplitude may also be tapered at the beginning and end points of the sweep, as shown by the dotted lines in FIG. 2. The acoustic waves are reflected by boundary conditions B.sub.1, B.sub.2, etc. in the earth formation and are received as reflected waves R.sub.1, R.sub.2, etc. at one or more surface receivers 15. The reflected waves R.sub.1, R.sub.2 are digitally recorded in recording apparatus 17 and later processed in processing apparatus 19. As well known, processing apparatus 19 produces useful information on subsurface lithology from the recorded reflection signals R.sub.1, R.sub.2, etc.

Assuming the acoustic wave signal, i.e., the vibratory signal, has a power spectrum function S(f), the reflection signal amplitude from a specific target depth is .sqroot.S(f) B(f), where B(f) represents a reflection wave amplitude attenuation characteristic. The reflection signal .sqroot.S(f) B(f) also has associated therewith a background noise power spectrum .sqroot.n(f). Both the reflection signal amplitude spectra and the background noise component for the FIG. 2 transmitted signal are illustrated in FIG. 3.

After an autocorrelation function is applied (multiplication of the reflection signal by the square root of the sweep signal S(f)) to the recorded and received reflected waves by processing apparatus 19, the amplitude spectra becomes S(f)B(f), while the noise amplitude spectra becomes .sqroot.n(f)S(f); both are illustrated in FIG. 4. Subsequent filtering, such as by using deconvolution techniques, will flatten the reflection signal spectrum for further processing, resulting in a Klauder-like wavelet. After filtering, the signal amplitude spectra has a uniform magnitude in the frequency band (f.sub.1 -f.sub.2) of interest. For simplicity, it will be assumed this magnitude is unity. However, the noise amplitude spectra is now ##EQU1## which is shown in FIG. 5. The corresponding noise power spectra ##EQU2## and unit magnitude signal power spectra are illustrated in FIG. 6. As shown therein, the amount of noise present is quite significant.

The method of the invention provides a seismic prospecting method using a unique vibratory signal sweep so as to minimize the noise energy content of the reflection signals as much as possible, and increase the signal-to-noise ratio, all within the physical limitations of the vibratory apparatus 13. The unique vibrator signal sweep produced in accordance with the invention can be represented as a power spectrum in the frequency domain S(f), which must meet the following conditions: ##EQU3## where S(f) is the vibratory sweep signal to be designed to optimize the signal-to-noise ratio and the function L(f) and the constant C represent physical limitations of the vibratory apparatus 13. L(f) is a function representing the vibrating time required for each frequency, such that a unit amount of vibratory energy can be supplied to the earth. This function can be experimentally determined for a given vibratory apparatus.

To optimize the signal-to-noise ratio, the noise energy content ##EQU4## must be minimized under the constraints of Equation (1). Using a Lagrange multiple formulation, the following equation can be written from Equation (1) and the noise power spectra ##EQU5## which is to be minimized. In this equation, .lambda. is the Lagrange multiplier, i.e., an abitrary constant scaler factor, and f.sub.1, f.sub.2 define the limits of the frequency band of interest.

Using variational calculus, the first variation is given by ##EQU6## If .delta.I=0 for any variation of .delta.S(f), then ##EQU7## at each frequency within the frequency band (f.sub.1, f.sub.2).

Or, stated otherwise, ##EQU8## in order to minimize noise and maximize the signal-to-noise ratio. This equation specifies that more vibratory energy (power) should be supplied at the frequency where background noise is high and attenuation is severe. Thus, when the functions n(f), B(f) and L(f) are known, the optimum power spectrum of the sweep is given by Equation 5.

Moreover, Equation 5 can be further simplified if certain assumptions are made. First, that the physical limitation function L(f)=1, which, in most cases, is true; and, second, that the noise power spectrum n(f)=n, a constant within the frequency band f.sub.1 -f.sub.2 of interest. The second assumption is also generally true, as the noise component is typically uniform in amplitude across a frequency band of interest. With these assumptions, Equation 5 can be rewritten as: ##EQU9##

The attenuation function B(f) from a vibratory source to a receiver can be a measured function for a particular area of interest or it can be represented by a constant Q characteristic as follows:

B(f)=e.sup.-.alpha.f (7) When the attenuation function B(f) is represented by Equation 7, Equation 6 reduces to ##EQU10## To determine the time sweep function f(t) for the vibratory apparatus 13, which is required to produce the power spectrum function S(f) of Equation 8, it is necessary to make further assumptions about the variations of amplitude as a function of frequency in the time domain. If the assumption of a uniform amplitude is made, more energy can be pumped into the ground at a specific frequency by keeping the vibrator operating at that frequency a proportionally longer time span. For simplicity, amplitude tapering at the beginning and end of the sweep is not considered in the design of the time sweep function. To derive the time sweep function f(t), it is noted that (as illustrated in FIG. 7) the energy at frequency f within a narrow band .DELTA.f is proportional to ##EQU11## This quantity should be equated to S(f).DELTA.f for optimum noise rejection at the correlated signal. Thus, ##EQU12## Integrating, one obtains, ##EQU13## where c.sub.1, and the factor .sqroot.n/.lambda. can be determined from the design parameters T, f.sub.1, and f.sub.2, where T is sweep duration and f.sub.1 and f.sub.2 define the sweep boundaries. To determine C.sub.1 and .sqroot.n/.lambda., the following two conditions must be satisfied: Condition 1. When t=0, f=f.sub.1. Condition 2. When t=T, f=f.sub.2.

From Condition 1: ##EQU14##

Thus, ##EQU15##

Condition 2 leads to: ##EQU16## or, ##EQU17## So, ##EQU18## or ##EQU19##

Simplifying, produces: ##EQU20##

This function is illustrated in FIG. 8.

A SPECIFIC EXAMPLE

A specific example of the invention is the design of a 16 second vibratory sweep from 15 Hz-100 Hz. From examining the amplitude spectrum of a dynamite explosion, which is computed from a time window from 2 seconds to 3 seconds (the target reflection is about 2.5 sec.) about the explosive time, the amplitude attenuation function (B(f)) is estimated to be flat and approximately 2.63 from 15 Hz to 100 Hz. That is,

B(f)=e.sup..alpha.(f.sbsp.2.sup.-f.sbsp.1.sup.) =2.63 (21)

Since f.sub.1 =15 Hz, and f.sub.2 =100 Hz, then

.alpha.=0.0114 sec. Also, e.sup..alpha.f.sbsp.2 =3.119, e.sup..alpha.f.sbsp.1 =1.1865, and T=16 sec. The sweep function f(t) (Equation 20) then becomes: f(t)=87.72[ln (t+9.8235)-2.114] (22) where t is in seconds, f is in Hertz and ln is the natural logarithm.

To further illustrate the invention, two other vibratory sweep functions f(t) can also be designed for the same target depth at the same area (same attenuation B(f) and noise n(f) functions), but with

f.sub.1 =15 Hz, f.sub.2 =90 Hz and T=16 seconds, and f.sub.1 =15 Hz, f.sub.2 =90 Hz and T=21 seconds

1. For

f.sub.1 =15 Hz f.sub.2 =90 Hz T=16 seconds, then, f(t)=87.72[ln (t+11.84)-2.3]. (23)

2. For

f.sub.1 =15 Hz f.sub.2 =90 Hz T=21 seconds, then, f(t)=87.72[ln (t+15.54)-2.572]. (24)

Once an optimal vibratory signal sweep has been designed, the vibratory signal is applied to vibratory apparatus 13 in the same manner as the more conventional sweep signals to generate, at one or more receivers 15, signals which, when processed by processing apparatus 19 and spectrum flattened (equalized), have an improved signal-to-noise ratio. The amplitude of the vibratory sweep signal is substantially constant over substantially the entire frequency spectrum defined by frequencies f.sub.1 and f.sub.2. However, amplitude tapering at the beginning and end of the vibratory signal sweep may also be applied to smooth the frequency response power function S(f) at f.sub.1 and f.sub.2 and to produce a sharper correlated vibratory signal sweep signature.

FIG. 9 illustrates in block diagram format the signal generating apparatus 11 required to produce the sweep power spectrum function S(f) described above with reference to Equation (5). It includes frequency domain function generators n(f), B(t) and L(t), denoted as elements 21, 23 and 25, a divider 27, a (1/.lambda.) multiplier 29, a square root circuit 31 and an inverter/multiplier 33. The various function generators can be configured as ROM stored digitized functions which are read out to a digital-to-analog converter.

In the time domain, once a specific time sweep function f(t) is determined for a specific application, the signal generating apparatus 11 can be formed as shown in FIG. 10. The time sweep function f(t) for a particular application can be digitally stored in ROM 41 and read out by sequential address signals generated by counter 43, counting the output pulses of oscillator 45. The digitized function f(t) is then converted into an analog signal by digital-to-analog converter 47, which controls the sweep rate of sweep signal generator 49 to produce an output signal having the desired frequency domain power spectrum S(f).

Although a preferred embodiment of the method and apparatus of the invention has been described above, it should be apparent that many modifications can be made thereto without departing from the spirit and scope of the invention. Accordingly, the invention is not limited by the foregoing description, but is only limited by the scope of the appended claims.

Claims

1. A method for seismic surveying using a vibratory signal sweep comprising the steps of:

determining a frequency domain subsurface amplitude attenuation function B(f) for a frequency spectrum existing between frequency components f.sub.1 and f.sub.2 for a subsurface area to be surveyed;

determining a frequency domain noise component function n(f) associated with said area for said frequency spectrum; and

producing a vibratory sweep signal over said frequency spectrum which has a frequency domain power spectrum function S(f) which is a function of said subsurface amplitude attenuation function B(f) and said noise component function n(f).

2. The method of claim 1, further comprising the steps of:

providing said vibratory sweep signal to an earth vibratory apparatus; and

vibrating the earth with said apparatus.

3. A method as in claim 1, wherein said vibratory sweep signal power spectrum function S(f) is a function proportional to the square root of said noise component function n(f) and inversely proportional to said amplitude attenuation function B(f).

4. A method according to claim 3, wherein said vibratory sweep signal power spectrum function is defined as follows: ##EQU21## where .lambda. is an arbitrary scale factor.

5. A method according to claim 3, further comprising the step of determining a function L(f) representing the vibrating time required for a particular vibrating apparatus to produce a unit amount of vibratory energy into the earth for each frequency within said spectrum, and wherein said vibratory sweep signal function S(f) is proportional to the square root of a quantity obtained by dividing said noise component function n(f) by said physical limitation function L(f).

6. A method according to claim 5, wherein said vibratory signal sweep power spectrum function is as follows: ##EQU22## where .lambda. is an arbitrary scale factor.

7. A method according to claim 2, further comprising the steps of

receiving and recording a reflection signal produced by vibrating the earth with said sweep signal power spectrum function S(f).

8. A method according to claim 4, wherein said vibratory signal sweep power spectrum function is as follows: ##EQU23## where n(f) is represented as a noise component constant n.

9. A method according to claim 1, wherein said power spectrum function S(f) is produced in accordance with a time function f(t) for operation of a sweep signal generator which is defined as follows: ##EQU24## where t represents the instantaneous time of the sweep, T represents the duration of the sweep, f.sub.1 and f.sub.2 are the starting and ending frequencies of the sweep and .alpha. is derived from the equation B(f)=e.sup..alpha.(f.sbsp.2.sup.-f.sbsp.1.sup.).

10. A method as in claim 1, wherein said vibratory sweep signal has a substantially constant amplitude over at least a substantial portion of the frequency spectrum defined by f.sub.1 and f.sub.2.

11. A method as in claim 10, wherein said vibratory sweep signal is amplitude tapered at the beginning and ending portions thereof and otherwise has said constant amplitude over said frequency spectrum.

12. A method as in claim 1, wherein said amplitude attenuation function B(f) is a constant Q function of the form B(f)=e.sup.-.alpha.f.

13. A method as in claim 1, wherein said noise function n(f) is defined as a constant value over said frequency spectrum.

14. An apparatus for vibratory sweep signal seismic prospecting comprising:

means for generating a vibratory sweep signal over a frequency spectrum existing between a first frequency f.sub.1 and a second frequency f.sub.2, said vibratory sweep signal having a frequency domain power spectrum function S(f) which is a function of a frequency domain amplitude attenuation function B(f) and a frequency domain noise component function n(f) over said frequency spectrum for a subsurface area to be surveyed; and

means responsive to said signal generating means for vibrating an earth formation with said vibratory sweep signal.

15. An apparatus as in claim 14, wherein said vibratory sweep signal power spectrum function S(f) is a function proportional to the square root of said noise component function n(f) and inversely proportional to said amplitude attenuation function B(f).

16. An apparatus according to claim 15, wherein said vibratory sweep signal power spectrum function is defined as follows: ##EQU25## where .lambda. is an arbitrary scale factor.

17. An apparatus as in claim 14, wherein said frequency domain power spectrum function S(f) is further proportional to a frequency domain function L(f) which represents the vibrating time required by said vibrating means for each frequency in said spectrum to produce a unit amount of energy into the earth.

18. An apparatus according to claim 17, wherein said vibratory signal sweep power spectrum function is as follows: ##EQU26## where .lambda. is an arbitrary scale factor.

19. An apparatus as in claim 14, further comprising means for receiving and recording a reflection signal produced by vibrating the earth with said vibratory sweep signal.

20. An apparatus according to claim 16, wherein said vibratory signal sweep power spectrum function is as follows: ##EQU27## where n(f) is represented as a noise component constant n.

21. An apparatus according to claim 14, wherein said power spectrum function S(f) is produced in accordance with a time function f(t) for said generating means which is defined as follows: ##EQU28## where t represents the instantaneous time of the sweep, T represents the duration of the sweep, f.sub.1 and f.sub.2 are the starting and ending frequencies of the sweep and .alpha. is derived from the equation B(f)=e.sup..alpha.(f.sbsp.2.sup.-f.sbsp.1.sup.).

22. An apparatus as in claim 14, wherein said vibratory sweep signal has a substantially constant amplitude over at least a substantial portion of the frequency spectrum defined by f.sub.1 and f.sub.2.

23. An apparatus as in claim 22, wherein said vibratory sweep signal is amplitude tapered at the beginning and ending portions thereof and otherwise has said constant amplitude over said frequency spectrum.

24. An apparatus as in claim 14, wherein said amplitude attenuation function B(f) is a constant Q function of the form B(f)=e.sup.-.alpha.f.

25. A method as in claim 14, wherein said noise function n(f) is defined as a constant value over said frequency spectrum.