The invention generally relates to a dithered slip sweep vibroseis acquisition system and technique.
Seismic exploration involves surveying subterranean geological formations for hydrocarbon deposits. A survey typically involves deploying seismic source(s) and seismic sensors at predetermined locations. The sources generate seismic waves, which propagate into the geological formations creating pressure changes and vibrations along their way. Changes in elastic properties of the geological formation scatter the seismic waves, changing their direction of propagation and other properties. Part of the energy emitted by the sources reaches the seismic sensors. Some seismic sensors are sensitive to pressure changes (hydrophones) and others are sensitive to particle motion (e.g., geophones). Industrial surveys may deploy only one type of sensors or both. In response to the detected seismic events, the sensors generate electrical signals to produce seismic data. Analysis of the seismic data can then indicate the presence or absence of probable locations of hydrocarbon deposits.
One type of seismic source is an impulsive energy source, such as dynamite for land surveys or a marine air gun for marine surveys. The impulsive energy source produces a relatively large amount of energy that is injected into the earth in a relatively short period of time. Accordingly, the resulting data generally has a relatively high signal-to-noise ratio, which facilitates subsequent data processing operations. The use of an impulsive energy source for land surveys may pose certain safety and environmental concerns.
Another type of seismic source is a seismic vibrator, which is used in connection with a “vibroseis” survey. For a seismic survey that is conducted on dry land, the seismic vibrator imparts a seismic source signal into the earth, which has a relatively lower energy level than the signal that is generated by an impulsive energy source. However, the energy that is produced by the seismic vibrator's signal lasts for a relatively longer period of time.
In an embodiment of the invention, a technique includes generating vibroseis sweeps for a vibroseis survey to produce seismic data acquired in response to seismic signals produced by the sweeps. The generation of the vibroseis sweeps including temporally arranging the sweeps into time-overlapping groups. The technique includes regulating a timing of the groups relative to each other based on a slip time. The technique also includes regulating a timing of the sweeps of each group such that consecutive sweep firings of each group are spaced apart by a time substantially less than the slip time.
In another embodiment of the invention, a technique includes receiving seismic data acquired in response to vibroseis sweeps arranged in time-overlapping groups. The technique includes processing the seismic data based on a slip time separation technique to separate the data into first datasets that are each associated with a different one of the groups. The technique includes processing the first datasets based on a source dithering separation technique to separate each first dataset into second datasets that are each associated with a different one of the sweeps.
In another embodiment of the invention, a system includes actuators and a control subsystem. The actuators generate vibroseis sweeps for a vibroseis survey to produce seismic data acquired in response to seismic signals produced by the sweeps. The vibroseis sweeps are arranged in time-overlapping groups. The control subsystem regulates a timing of the groups relative to each other based on slip time and regulates a timing of the sweeps of each group such that consecutive sweep firings of each group are spaced apart by a time substantially less than the slip time.
In yet another embodiment of the invention, a system includes an interface and a processor. The interface receives seismic data acquired in response to vibroseis sweeps temporally arranged in time-overlapping groups. The processor processes the seismic data based on a slip time separation technique to separate the data into first datasets that are each associated with a different one of the groups; and the processor processes the first datasets based on a source dithering separation technique to separate each first dataset into second datasets that are each associated with a different one of the sweeps.
Advantages and other features of the invention will become apparent from the following drawing, description and claims.
Referring to
For purposes of generating the signal 15, the seismic vibrator 10 may contain an actuator (a hydraulic or electromagnetic actuator, as examples) that drives a vibrating element 11 in response to a sweep pilot signal (called “DF(t)” in
Among its other features, the seismic vibrator 10 may include a signal measuring apparatus 13, which includes sensors (accelerometers, for example) to measure the signal 15 (i.e., to measure the output ground force of the seismic vibrator 10). As depicted in
The vibrating element 11 contains a reaction mass that oscillates at a frequency and amplitude that is controlled by the DF(t) pilot signal: the frequency of the DF(t) signal sets the frequency of oscillation of the reaction mass; and the amplitude of the oscillation, in general, is controlled by a magnitude of the DF(t) signal. During the generation of the sweep, the frequency of the DF(t) signal transitions (and thus, the oscillation frequency of the reaction mass transitions) over a range of frequencies, one frequency at time. The amplitude of the DF(t) signal may be linearly or non-linearly varied during the generation of the sweep pursuant to a designed amplitude-time envelope.
It is noted that unlike the seismic vibrator 10, a seismic vibrator may alternatively be constructed to be located in a borehole, in accordance with other embodiments of the invention. Thus, seismic sensors, such as geophones, may alternatively be disposed in a borehole to record measurements produced by energy that is injected by borehole-disposed vibrators. Although specific examples of surface-located seismic vibrators and seismic sensors are described herein, it is understood that the seismic sensors and/or the seismic vibrators may be located downhole in accordance with other embodiments of the invention.
Due to the mechanics and movement of the seismic vibrator, the overall time consumed in generating a vibroseis sweep significantly exceeds the sweep length, or duration, which is just one component of the overall time. For example, the overall time involved in generating a particular vibroseis sweep includes a time associated with deploying the base plate (such as the base plate 12 depicted in
For purposes of further increasing the acquisition efficiency when multiple seismic vibrators are used, a “slip sweep” technique may be used, as described in, for example, Rozemond, H. J., 1996, Slip-sweep acquisition, 66th Annual International Meeting, SEG, Expanded Abstracts, 64-67. In the slip sweep technique, a particular sweep begins without waiting for the previous sweep to terminate. In the absence of harmonic noise, if the time interval between the beginning, or firing, of consecutive sweep sequences (called the “slip time”) is greater than the listening time, then the seismic responses to the consecutive sweep sequences do not overlap in the time-frequency domain, which facilitates separation of the measurements.
If, however, harmonic noise is present, the seismic responses do overlap in the time-frequency domain and correspondingly deteriorate the quality of the acquired seismic data. U.S. Pat. No. 7,050,356, entitled “METHOD OF SEISMIC SURVEYING,” which issued on May 23, 2006, and is hereby incorporated by reference discloses a technique to attenuate this harmonic noise when the slip time is above a certain value that does not depend upon the acquisition parameters. PCT Publication No. 2008/110743, entitled “VIBROSEIS ACQUISITION METHOD,” which has a filing date of Mar. 9, 2007, and is hereby incorporated by reference in its entirety, extends this technique to smaller slip times, which are equal to or greater than the listening time.
In accordance with embodiments of the invention described herein, the efficiency of the vibroseis survey is further enhanced by combining the slip sweep technique with a source dithering technique. With source dithering, seismic sources are operated close in time, which causes the corresponding seismic signals to interfere with each other for at least part of the record in the recorded time-frequency space. However, as described in U.S. patent application Ser. No. 11/964,402, entitled “SEPARATION OF SIMULTANEOUS SOURCE DATA,” which was filed on Dec. 26, 2007, and is hereby incorporated by reference in its entirety, the data acquired from simultaneously or near simultaneously fired interfering seismic sources may be processed via a source separation technique (called a “source dithering separation technique” herein) to separate the acquired data according to the producing source. More specifically, in this technique, the acquired seismic data undergoes a preprocessing step to produce corresponding datasets that each correspond to one of the interfering sources and do not correspond to the other interfering sources. In accordance with embodiments of the invention, the source dithering involves slightly varying the times at which the sweep sequences in a dithered group are initiated, or fired, via a predetermined schedule or through a randomly or pseudo randomly generated timing.
The initiation, or firing, times are separated by significantly less durations than the above-described slip times. More specifically, the sources in a dithered group are varied by dither times. Both the slip and dither times are time intervals between the beginning times of consecutive sweeps. The slip time is always larger than the listening time, whereas the dither time is significantly smaller than the listening time.
A slip time separation technique would not sufficiently separate the measurements produced by a particular dithered group of sweeps. However, the measurements acquired from a source dithered group may be source separated as follows. When traces are collected into a domain that includes many firings of each interfering source and are aligned such that time zero corresponds to the firing time for a specific source, then the signal from the source appears coherent while the signal from the other sources appears incoherent. As described in U.S. patent application Ser. No. 11/964,402, although the signal from the specific source appears to be coherent while the signal from the other sources appears to be incoherent, the apparently incoherent signals are not mathematically incoherent due to the known timing delays. Therefore, the recorded data may be decomposed into data associated with each source using, for example, Radon transforms.
In accordance with embodiments of the invention described herein, a vibroseis-based seismic survey uses both the source dithering and slip sweep techniques to enhance the efficiency of the survey. The technique is illustrated using an exemplary system 100 that is schematically depicted in
The vibroseis sweeps are temporally arranged in time-overlapping groups; and the groups are generated successively and overlap each other in time such that the initiation, or firing, of the first sweep of each group is separated by at least a slip time from the initiation, or firing, of the first sweep of the previously generated group. For the example that is depicted in
In the absence of harmonic noise and if the earth response has a finite length, the slip time condition 120 imposes a time interval equal to the earths impulse response. In presence of harmonic noise, a more restrictive slip time condition is applied. Due to the slip time condition 120, a slip time separation technique may be used to separate the acquired seismic data into datasets, where each dataset is associated with the measurements produced by a different pair 102, 104, 106, 108, 110 or 112. The slip sweep separation may involve separating the acquired seismic data according to the pairs based on the time-frequency separation that is due to the slip condition 120. The slip sweep separation may involve additional processing if significant harmonic energy is present, as described in the references mentioned above.
The sweeps of each pair are dithered, which means that the sweeps of each pair are fired at different times, but the firing times of the sweeps are substantially closer together in time than the slip time. Therefore, the seismic signals that are produced by the sweeps of each pair substantially interfere with each other in the time-frequency space, and a source dithering separation process is used to source separate the signals acquired for each group. Thus, the second sweeps 102b, 104b, 106b, 108b, 110b and 112b are constrained by a dithering condition 124.
Due to the dithering condition 124, a source dithering separation technique, such as the one disclosed in U.S. patent application Ser. No. 11/964,402 (as a non-limiting example), may be used to source separate the measurements attributable to each pair. The source dithering separation technique may be carried out in the common receiver domain.
As a more specific example, a dithered slip sweep sequence may be generated in the following manner, in accordance with embodiments of the invention. The sequence is partially defined by a dithering time called “ΔTi” which may be positive or negative and represents the firing time of the second sweep in the dithered pair (where “i” is the dithered pair index), relative to the firing time of the master sweep of the pair. The minimum and maximum dithering times are denoted by “Tmin” and “Tmax,” respectively, and the following conditions for the firing sequence applies. The firing times for both sweeps 104a and 104b of the pair 104 need to be more than the slip time, ST, after the later of the firing times of the sweeps 102a and 102b. Also, the sweep 104b fires at time ΔT2 after the sweep 104a. The earliest firing time for 104a is therefore given by the sum of ST and maximum(0, ΔT1) less minimum(0, ΔT2). The worst case occurs when ΔT1=Tmax and ΔT2=Tmin.
The firing times may be made independent of the actual dither times by assuming the worst case for each dither, in which case the sweep 104a fires at the sum of ST and maximum(0, Tmax) less minimum(0, Tmin) after 102a. There are three cases for this: Tmax>Tmin>0, Tmax>0>Tmin and 0>Tmax>Tmin, with corresponding firing delays of ST+Tmax, ST+Tmax−Tmin and ST-Tmin, respectively. In the special case for which all the ΔTi are positive, the earliest firing time for the sweep 104b becomes ST+ΔT1, and this, along with ST=6 second (s), is the case shown in
The example that is depicted in
Furthermore, for this example, the six vibrators whose sweep start time is dictated by the slip-sweep condition 120 are located on different shot lines.
The acquisition parameters listed above imply a cycle time (called “CT”) per vibrator of CT=36.0+Σi=1.6 ΔTi seconds (s). This cycle time may be kept by each vibrator because the total move-up time (SL+PU+PD+MU=31s) is less than the cycle time CT.
In the above example, the maximum number of successive sweeps that may interfere with a given sweep is five. For example, the sensed seismic signals produced by the sweeps 102b, 104a, 104b, 106a and 106b are recorded before the end of the listening time of 102a, and therefore may interfere with it; and the slip sweep and source dithering separation techniques that are disclosed herein may be applied to remove these interfering signals.
The above example is included for illustrative purposes only, as other embodiments are contemplated and are within the scope of the appended claims. For example, other variants include vibrators whose start times are slip-sweep constrained and located on the same shot line and/or pairs of vibrators that form dithered pairs on different source lines. As another example, the actual vibrators to which either the dithering or slip-sweep condition is applied may change during the survey; and the recording system stores in memory the previous two start times to decide whether a dithering or a slip-sweep condition must be applied to the start time of the next sweep. As yet another example, there may be more than two sweeps in a dithered group.
To summarize,
Additionally, the parallelograms 3042, 3062 and 3082 refer to the second harmonics due to the second, third and fourth shots, respectively.
For this example, the firing time of shot one is at time zero, and the firing time of shot two is the dither time, or ΔT1. The firing time for shot three is MST+ΔT1. The firing time of shot three is therefore linked to the firing time of shot two by a slip-time condition, whereas the firing time of shot two is linked to the firing time of shot one by a dithering time condition.
The fundamental response 3021 in the interval of interest (i.e., the interval of interest between time zero and time SL+LTIM) is contaminated by the fundamental responses 3041, 3061 and 3081 and by the harmonic responses 3022, 3042 3062 and 3082. The contamination that does not overlap in time-frequency domain may be eliminated using such techniques as those described in U.S. Pat. No. 7,050,356 and PCT Publication No. 2008/110743. The contamination that overlaps in the time-frequency domain may be eliminated with a three or higher fold (depending on the dither time) dithering separation approach, as disclosed in U.S. patent application Ser. No. 11/964,402.
In accordance with some embodiments of the invention, the slip sweep separation is followed by the dither separation. The contamination caused by fundamental responses 3061 and 3081 and the harmonic responses 3022, 3042 3062 and 3082 is eliminated using the slip sweep approach. The contamination due to the fundamental response 3041 (the most severe one) is successively attenuated with a two-fold dithering separation approach. The 3041 contamination may be eliminated by applying the dithering separation technique in the common receiver gather domain, such as the technique disclosed in U.S. patent application Ser. No. 11/964,402.
In other embodiments of the invention, the contamination caused by the responses 3061 3081, 3022 and 3042 is eliminated using the slip sweep approach. For this approach, the contamination due to the response 3041 (the most severe one), 3062 and 3082 is successively attenuated with a four-fold dithering separation approach.
In another embodiment of the invention, a dithering separation approach may be first applied in the interval of interest for all of the contaminating sources, and this approach has been followed by a slip sweep approach to eliminate the remaining contamination.
In another variation, the slip sweep and dither approaches may be simultaneously applied.
In yet another variation, all of the coherent sources in the time interval from time zero to time SL+LTIM are simultaneously migrated. The time-frequency representation of all of the coherent sources is in fact known with a high level of accuracy due to the fact that the signal emitted by each vibrator is controlled by a feedback system that minimizes the difference between the transmitted ground-force and a reference signal.
In yet another approach, the dithering time is frequency dependent. Such dependency may be obtained in the simplest with constant phase shift between sweeps.
Referring to
It is noted that
For the processing system 400 that is depicted in
Other embodiments are contemplated and are within the scope of the appended claims. For example, in accordance with some embodiments of the invention, because of logistics, the dither and slip sweep sequences may not be interlaced. For example, because of logistics, a cycling sequence like the following may be applied: dither-slip sweep-slip sweep, etc. This may happen if the third vibrator to which a dither condition is supplied is not ready and therefore, a less strict slip time condition is applied.
As an example of another variation, the vibrators may be fired in a different order than the order described in connection with
As another example, in some embodiments of the invention, a set of different initial phases may be applied to the sweeps. As an example, this set of initial phases may be determined according to the technique disclosed in U.S. Pat. No. 4,982,374, entitled, “METHOD OF SOURCE CODING AND HARMONIC CANCELLATION FOR VIBRATIONAL GEOPHYSICAL SURVEY SOURCES,” which issued on Jan. 1, 1991, and is hereby incorporated by reference in its entirety.
While the present invention has been described with respect to a limited number of embodiments, those skilled in the art, having the benefit of this disclosure, will appreciate numerous modifications and variations therefrom. It is intended that the appended claims cover all such modifications and variations as fall within the true spirit and scope of this present invention.
This is a divisional application of co-pending U.S. patent application Ser. No.: 12/247,284 to Claudio Bagaini, et al., filed on Oct. 8, 2008, and entitled “Dithered Slip Sweep Vibroseis Acquisition System and Technique,” which is hereby incorporated in its entirety for all intents and purposes by this reference.
Number | Name | Date | Kind |
---|---|---|---|
6519533 | Jeffryes | Feb 2003 | B1 |
20050128874 | Herkenhoff | Jun 2005 | A1 |
20070195644 | Marples | Aug 2007 | A1 |
20100085836 | Bagaini et al. | Apr 2010 | A1 |
20120008462 | Liu | Jan 2012 | A1 |
Number | Date | Country | |
---|---|---|---|
20150177397 A1 | Jun 2015 | US |
Number | Date | Country | |
---|---|---|---|
Parent | 12247284 | Oct 2008 | US |
Child | 14637363 | US |