Technical Field
Embodiments of the subject matter disclosed herein generally relate to methods for seismic data acquisition using clustered seismic vibrators that are placed at different vibration points within a predetermined distance from one another and are activated simultaneously.
Discussion of the Background
In the oil and gas industry sector particularly, seismic surveys are commonly used to search for and evaluate subterranean hydrocarbon deposits. Seismic surveys are exploration methods that gather and record data related to seismic wave reflections from interfaces between geological layers. Seismic data is then used to create detailed models of the underlying geological structure.
In operation, source 110 generates seismic waves that may include a surface waves 140 and body waves 160 that may be partially reflected at an interface 170 between two geological layers inside which the seismic waves propagate with different speeds. Each receiver 120 receives the full wavefield (i.e., both surface and body waves) and converts it into an electrical signal.
Source 110 may be a vibratory source. Vibratory sources are actuated in sweeps lasting generally between 2 and 60 seconds, during which time the vibratory sources generate seismic waves whose frequencies vary over a given range (e.g., 2-200 Hz). Receivers 120 record data for a period during and after the sweep time has ended. This period after the sweep time is known as listening time. During listening time, source 110 may be moved from one location to a next location according to a predetermined survey plan. The locations where seismic sources are activated are known as vibration points.
Trucks carrying seismic sources 230-241 move along shot lines SL, stopping at successive vibration points.
A shooting plan includes a large number of vibration points to acquire the information necessary to arrive at conclusions related to the presence and location of subterranean hydrocarbon deposits. Even if the combined sweep and listening time for each one of the vibration points is less than 30 s, acquiring seismic data for all vibration points in a typical survey plan takes a long time, which is expensive in terms of equipment and personnel. Meanwhile, there is a permanent desire to increase vibration points' density to achieve better images of the explored underground structures. Therefore, increasing productivity is an ongoing effort.
One method of increasing productivity of acquiring seismic data known as “slip-sweep” is described in the publication, Wams, J. and Rozemond J. (1998), “Recent Developments in 3-D acquisition using vibroseis in Oman”, Leading Edge 17 No. 8, pp. 1053-1063, the content of which is incorporated herein by reference. According to the slip-sweep method, a next sweep starts before the end of listening time related to a previous sweep. The time interval between successive sweeps is known as slip time. However, to be able to separate seismic data corresponding to the different sweeps, slip time is limited by the requirement to avoid overlapping harmonic energy generated during successive sweeps. For example, using the slip-sweep method for a slip time of 5 s with the plan in
Some harmonic noise reduction algorithms have been developed, for example, as described in U.S. Pat. No. 6,603,707 and U.S. Pat. No. 8,451,686, the contents of which are incorporated herein by reference. Using these harmonic noise reduction algorithms has allowed productivity to further increase.
Another technique of increasing seismic data acquisition productivity known as “independent simultaneous sweeping” (ISS) is described in U.S. Patent Application Publication No. 2012/0290213, the content of which is incorporated herein by reference. ISS uses statistical methods to remove interference due to unsynchronized overlapping sweeps of sources located at large distances there-between (e.g., 2 km).
However, as ISS separation remains imperfect for a large number of closely spaced sources it remains desirable to enhance the statistical source separation by using acquisition schemes specially designed for separation.
To increase data acquisition productivity, closely placed vibratory sources (e.g., less than 200 m apart) are actuated simultaneously to generate seismic vibrations encoded so as to enable separation of individual seismic data related to seismic vibrations produced by each of the simultaneously actuated sources.
According to one embodiment, there is a method for seismic exploration of a subsurface formation. The method includes deploying seismic receivers along receiver lines over the explored subsurface formation. The method further includes grouping vibratory sources in clusters so that at least one of the clusters includes several vibratory sources that are placed at less than a predetermined distance from one another while vibrating, the predetermined distance being less than an average distance between adjacent receiver lines. The method also includes actuating the clusters cyclically, so that, in each cycle, the vibratory sources within one of the clusters vibrate during a substantially same sweep time, being encoded to enable separation of recorded seismic data into subsets corresponding to the seismic vibrations.
According to another embodiment, there is a method for seismic exploration of a subsurface formation. The method includes deploying seismic receivers along receiver lines over the explored subsurface formation. The method further includes grouping vibratory sources in clusters, so that at least one of the clusters includes several vibratory sources that are placed less than a predetermined distance from one another while vibrating, the predetermined distance being less than an average distance between adjacent receiver lines. The method also includes actuating the clusters cyclically to produce seismic vibrations propagating in the subsurface formation, so that, in each cycle, the several vibratory sources in the at least one of the clusters vibrate during a substantially same sweep time, to generate seismic vibrations encoded to enable, for vibration frequencies higher than a threshold frequency, separation of recorded seismic data into subsets corresponding to each of the vibratory sources.
According to another embodiment, there is a method for seismic exploration of a subsurface formation. The method includes simultaneously actuating vibratory sources placed at less than 30 m from one another and acquiring seismic data. Here, individual seismic vibrations produced by the vibratory sources are encoded to enable separation of individual seismic data corresponding to at least one of the vibratory sources, from the seismic data, and the encoding is achieved by one of short time static differences, different constant phases and a combination of time and phase differences.
The accompanying drawings, which are incorporated in and constitute a part of the specification, illustrate one or more embodiments and, together with the description, explain these embodiments. In the drawings:
The following description of the exemplary embodiments refers to the accompanying drawings. The same reference numbers in different drawings identify the same or similar elements. The following detailed description does not limit the invention. Instead, the scope of the invention is defined by the appended claims. The following embodiments are discussed using terminology of a land data acquisition system. However, similar methods may also be pertinent to marine seismic data acquisition and to data acquisition and processing based on similar approaches, such as measuring electromagnetic or ultrasound responses of an unknown structure to excitations injected into the structure at different locations.
Reference throughout the specification to “one embodiment” or “an embodiment” means that a particular feature, structure or characteristic described in connection with an embodiment is included in at least one embodiment of the subject matter disclosed. Thus, the appearance of phrases “in one embodiment” or “in an embodiment” in various places throughout the specification is not necessarily referring to the same embodiment. Further, the particular features, structures or characteristics may be combined in any suitable manner in one or more embodiments.
To increase productivity, vibratory sources are grouped in clusters of two or more that emit simultaneously (i.e., during substantially the same sweep time). The vibratory sources in a cluster start sweeping at the same time except, in some embodiments, for small static time differences (of 50-1000 ms) used for encoding as explained below (these small static time differences being substantially less than the sweep time and constant along a trace). The seismic vibrations emitted during this sweep time are encoded to enable separation of individual seismic data corresponding to the seismic vibration emitted by each of the vibratory sources in the same cluster. Unlike in ISS methods where unsynchronized overlapping sweeps are produced by distant sources, according to various embodiments, in these embodiments, the vibratory sources are physically close, (e.g., at less than 200 m from one another) and are correlated time-wise.
Trucks carrying seismic sources 330-341 move along shot lines SL to place and activate the sources at successive vibration points, a grid of vibration points results from intersecting shot lines SL (parallel to receiver lines 322, 324, 326) with lines SN perpendicular to the shot lines SL. Since, for the embodiment illustrated in
Vibratory sources 330-341 are clustered (i.e., clusters 350, 352, 354, 356, 358, 360) so that sources in the same cluster are actuated simultaneously (i.e., during substantially the same sweep time). Cluster 350 includes sources 330 and 331, cluster 352 includes sources 332 and 333, etc. In
The clusters are actuated one after another, cyclically, according to a predetermined sequence or the source activation may be coordinated by a controller 370 (in this case, the sources do not have to be activated cyclically). Sources of at least one cluster are ready to be activated before a listening time corresponding to a previously activated cluster has ended. The sources may remain clustered throughout the survey. However, changes in this arrangement may be triggered by relief or occasional events (e.g., a source malfunction or additional sources).
Method 400 further includes grouping vibratory sources in clusters, at 420. In
Method 400 further includes actuating the clusters cyclically, at 430. During a cycle, each cluster is actuated once. The vibratory sources pertaining to a cluster vibrate during a substantially same sweep time (i.e., as previously discussed, short static time differences much smaller than the sweep time are considered insignificant). The seismic vibrations produced by the vibratory sources in the same cluster are encoded to enable separation of recorded seismic data into subsets, each subset corresponding to one or more of the seismic vibrations.
The seismic vibrations may be encoded by one of short time static differences, different constant phases, a combination of time and phase differences, etc. For example, considering that the vibrations are sinusoidal, first vibrations Y1 corresponding to a first source and second vibrations Y2 corresponding to a second source (where Y1 and Y2 may be pressure or transverse displacement) may be written as:
Y1(t)=A1(t)sin(2πf1(t)(t−t01)+φ1) (1)
and
Y2(t)=A2(t)sin(2πf2(t)(t−t02)+φ2) (2)
Where A1(t) and A1(t) are maxim amplitude of the first and second vibrations, respectively, f1(t) and f2(t) represent individual sweeps (i.e., frequency variation with time for the respective source), t is time, t01 and t02 are start times of the first and second vibrations, respectively, φ1 and φ2 are initial phases (which also can be function of time and frequency). If the sweeps of the two vibratory sources are the same then f1(t)=f2(t).
If the first seismic vibrations and the second seismic vibrations are encoded via short time static differences, it means that t01≠t02, while f1(t) may be equal to f2(t) and φ1=φ2=0.
If the first seismic vibrations and the second seismic vibrations are encoded via different constant phases, it means that φ1≠φ2 (but φ1 and φ2 are constant with time or frequency) while f1(t) may be equal to f2(t) and t01=t02. For example,
However, in another embodiment, the sign of the phase difference between the first and second seismic vibrations may regularly alternate between +90° or −90° from one cycle to a next cycle. Provided that two adjacent vibrated points have a close enough Earth response, such a phase rotation scheme leads to a well-conditioned linear system of the form: D=φ*E, where E contains the Earth response for each source, D are the recorded shot gathers and
One skilled in the art would appreciate that the hypothesis on locality (i.e., adjacent shots have a close reflectivity) is likely to hold in the Radon domain (frequency, slopes) for local transforms.
If the first seismic vibrations and the second seismic vibrations are encoded via a combination of time and phase differences, it means that φ1=φ2 and t01≠t02 while f1(t) may be equal to f2(t).
The above description discusses an embodiment characterized by two-source clusters. However, in other embodiments, one or more clusters may have more than two sources. For example, in an embodiment with three-source clusters, source phases may rotate, taking, for example, values of
for three consecutive shot locations (i.e., the constant phases are pre-determined according to geographical shot locations). In this matrix, each line corresponds to a successive vibration of the clustered sources.
Additionally,
Since encoding makes it possible to separate (i.e., extract from the acquired seismic data) the first seismic data corresponding to the first vibrations and the second seismic data corresponding to the second seismic vibrations, each vibratory source is moved from a vibration point to a next vibration point during each cycle (i.e., unlike in U.S. Pat. No. 7,050,356 where data for multiple repeated sweeps with different parameters are acquired to perform such a separation).
The clustered sources may be selected to be as close as possible from one another. Related to this feature, in one embodiment, no vibratory source of a cluster may be intercalated between vibratory sources of any other cluster.
Similar to the method described in U.S. Pat. No. 8,451,686 for individual sources, in each cycle, the clusters may be activated at one of a series of times defined relative to the beginning of the cycle. The time intervals between successive times in the series may be substantially equal. In one embodiment, the clusters are activated in the same order in each cycle. In another embodiment, the clusters are activated in a different order in successive cycles (granted that sources in a cluster have time to move to new vibration points before being re-activated in a new cycle).
The above embodiments of simultaneous activation of clustered sources that are physically close may be combined with the slip-sweep method. Thus in one embodiment, at least one of the clusters starts vibrating before the end of seismic receivers' listening time for a previously actuated cluster. In this case, the acquired seismic data may be processed first to extract individual cluster seismic data. For example, an algorithm similar to the one described in U.S. Pat. No. 8,451,686 (HPVA) may be used to extract the individual cluster seismic data. The resulting individual cluster seismic data is then processed to separate the first data and the second data.
In a test run using the survey plan in
Lately, extending the sweep frequency spectrum toward low frequencies (less than 10 Hz) has been considered advantageous. However, in order to be able to separate those low frequencies, the methods are restricted to the use of large static time shifts, which may be counterproductive. At those frequencies, two closely located vibratory sources (i.e., at a distance less than for example ¼ of low-frequency vibrations' wavelength) can emit simultaneously and in phase such that they have the effect of a single source.
According to another embodiment, a sweep of the clustered vibratory sources includes (A) a first period during which the first vibratory source and the second vibratory source in the same cluster have the same phase and vibrate with low frequencies less than or equal to a threshold frequency, and (B) a second period during which the first vibratory source and the second vibratory source have different phases and vibrate with frequencies higher than the threshold frequency. The seismic vibrations produced by different clustered seismic sources (emitting simultaneously) are then encoded to enable separation of seismic data corresponding to individual seismic vibrations at the frequencies higher than the threshold frequency. The advantage of this embodiment is that more low frequency signal is recovered. The signal to noise ratio of frequencies below the threshold will increase by a factor equal to the square root of the number of sources in the cluster, compared to the previous embodiments where phase and shifts are not equal for the whole sweep.
During the second period, the seismic vibrations emitted by different clustered sources may be encoded by one of short time static differences, different constant phases, and a combination of time and phase differences. In one embodiment, during the second period, seismic vibrations emitted by different clustered sources have different constant phases, but phase differences may change randomly in each shooting cycle.
Returning now at
Processor 702 may communicate with other internal and external components through input/output (I/O) circuitry 708 and bussing 710, to provide control signals and the like. For example, processor 702 may communicate with the sensors, electro-magnetic actuator system and/or the pressure mechanism of each source element. Processor 702 carries out a variety of functions as are known in the art, as dictated by software and/or firmware instructions. Server 701 may also include one or more data storage devices, including disk drives 712, CD-ROM drives 714, and other hardware capable of reading and/or storing information, such as a DVD, etc. In one embodiment, software for carrying out the above-discussed steps may be stored and distributed on a CD-ROM 716, removable media 718 or other form of media capable of storing information. The storage media may be inserted into, and read by, devices such as the CD-ROM drive 714, disk drive 712, etc. Server 701 may be coupled to a display 720, which may be any type of known display or presentation screen, such as LCD, plasma displays, cathode ray tubes (CRT), etc. Server 701 may control display 720 to exhibit images based on the acquired seismic data. A user input interface 722 is provided and may include one or more user interface mechanisms such as a mouse, keyboard, microphone, touch pad, touch screen, voice-recognition system, etc.
Server 701 may be coupled to other computing devices, via a network. The server may be part of a larger network configuration as in a global area network (GAN) such as the Internet 728, which allows ultimate connection to the sources, receivers, data storage unit, etc.
The disclosed embodiments provide methods of seismic data acquisition according to which clustered (i.e., physically close) sources are actuated simultaneously. It should be understood that this description is not intended to limit the invention. On the contrary, the exemplary embodiments are intended to cover alternatives, modifications and equivalents, which are included in the spirit and scope of the invention as defined by the appended claims. Further, in the detailed description of exemplary embodiments, numerous specific details are set forth in order to provide a comprehensive understanding of the claimed invention. However, one skilled in the art would understand that various embodiments may be practiced without such specific details.
Although the features and elements of the present exemplary embodiments are described in particular combinations, each feature or element may be usable alone without the other features and elements of the embodiments or in other various combinations with or without other features and elements disclosed herein.
The written description uses examples of the subject matter disclosed to enable any person skilled in the art to practice the same, including making and using the described devices or systems and performing any of the described methods. The patentable scope of the subject matter is defined by the claims, and may include other examples that occur to those skilled in the art. Such examples are intended to be within the scope of the claims.
This application claims priority and benefit from U.S. Provisional Patent Application No. 61/818,002, filed May 1, 2013, for “Simultaneous Vibroseis Acquisition,” the content of which is incorporated in its entirety herein by reference.
Filing Document | Filing Date | Country | Kind |
---|---|---|---|
PCT/EP2014/058825 | 4/30/2014 | WO | 00 |
Publishing Document | Publishing Date | Country | Kind |
---|---|---|---|
WO2014/177614 | 11/6/2014 | WO | A |
Number | Name | Date | Kind |
---|---|---|---|
4715020 | Landrum, Jr. | Dec 1987 | A |
6603707 | Meunier et al. | Aug 2003 | B1 |
7050356 | Jeffryes | May 2006 | B2 |
8451686 | Meunier | May 2013 | B2 |
20060034153 | Meunier et al. | Feb 2006 | A1 |
20090135671 | Meunier | May 2009 | A1 |
20110272207 | Meunier | Nov 2011 | A1 |
20110305105 | Dean | Dec 2011 | A1 |
20120008462 | Liu et al. | Jan 2012 | A1 |
20120290213 | Huo et al. | Nov 2012 | A1 |
Entry |
---|
International Search Report issued in corresponding International application No. PCT/EP2014/058825, date of completion of the International search Jan. 16, 2015. |
Written Opinion of the International Searching Authority issued in corresponding International application No. PCT/EP2014/058825, date of mailing of the International search report Jan. 26, 2015. |
Wams, Jan et al., “Recent developments in 3-D acquisition techniques using vibroseis in Oman,” The Leading Edge, vol. 17, No. 8, pp. 1053-1063, Tulsa, Oklahoma, USA, DOI: 10.1190/1.438089, Aug. 1998. |
Office Action dated Aug. 8, 2017 in related OM Application No. OM/P/2015/000280 (all references cited in OA were previously cited in IDS filed Oct. 21, 2015). |
Number | Date | Country | |
---|---|---|---|
20160077226 A1 | Mar 2016 | US |
Number | Date | Country | |
---|---|---|---|
61818002 | May 2013 | US |