Patent Grant
9594180
Embodiments of the subject matter disclosed herein generally relate to methods and systems for seismic data processing and, more particularly, to mechanisms and techniques for removing ghost reflections from seismic data.
Seismic data acquisition and processing techniques are used to generate a profile (image) of a geophysical structure (subsurface) of the strata underlying the land surface or seafloor. Among other things, seismic data acquisition involves the generation of acoustic waves and the collection of reflected/refracted versions of those acoustic waves to generate the image. This image does not necessarily provide an accurate location for oil and gas reservoirs, but it may suggest, to those trained in the field, the presence or absence of oil and/or gas reservoirs. Thus, providing an improved image of the subsurface in a shorter period of time is an ongoing process in the field of seismic surveying.
Mapping subsurface geology during exploration for oil, gas, and other minerals and fluids uses a form of remote sensing to construct two-dimensional or three-dimensional images of the subsurface. The process is known as seismic surveying, wherein an energy source transmits pressure pulses into the earth. These pressure pulses can be reflected by geological interfaces associated with the earth and subsequently recorded at the surface by arrays of detectors.
A contaminant which causes problems in processing and interpreting the subsurface data is a phenomenon known as a ghost reflection. Ghost reflections are caused by the sea-surface operating as a mirror and reflecting up-going pressure-waves. There are source-side ghosts caused by the source energy being reflected back from the sea-surface and there are receiver-side ghosts incident at the detectors as down going surface reflections of the up-going energy from the subsurface.
These ghost reflections do more than complicate the subsurface image with additional waves, the mirror effect which produces the ghost reflections changes the phase of the reflection by 180 degrees such that, in some circumstances, the energy constructively interferes with the desired signal to magnify it and, in other circumstances, destructively interferes with and destroys the desired signal.
Water surface reflections of very long wavelength, low frequency seismic waves, destructively interfere such that there is always a null or notch at zero Hertz. These notches and the associated filtering limit the extent to which subsurface reflections or events can be resolved. This damaging process results in blurred images at best, and at worst, fictitious reflections when the ghost energy lags significantly behind the primary reflection energy. To counter the effect of the notches, conventional seismic surveys are designed to use shallow source and receiver arrays to ensure that the second notch lies at high frequencies. However, the resultant sloped bandpass filter is a pervading problem as it causes seemingly irretrievable damage to low frequency information which is increasingly being sought in the industry to deliver extra value in the interpretation process. Furthermore, a relatively shallow tow of the receivers exposes them to the very noisy-environment that exists close to the sea surface which in itself can be problematic.
Most existing mechanisms to de-ghost seismic data fundamentally rely on recording alternate views of the same data. Data are deliberately acquired with different ghost characteristics so that when combined there is improved signal spectrum coverage and the damaging notches are filled.
There are three commonly used data acquisition variations designed to tackle the ghost problem: slanted streamers; parallel streamers arranged vertically one above another at the same horizontal position, known as “over-under”, and mixtures of different types of detectors (e.g., use of hydrophones with vector-sensors such as accelerometers) at a coincident position, are representative.
The former techniques exploit variations in the recorded ghost effect, which can be processed together to de-ghost the signals. The latter exploits the fact that up-going and down-going energy exhibit different polarity which one type of detector is able to observe, whereas the other type does not. This allows the ghost energy to be removed by careful summation of the two signals based on the fact that the ghost is of opposite polarity in one of the recorded datasets. The slanted streamer array technique deploys (current) standard streamer equipment, whereas the other two data acquisition techniques use an increased number of streamers or sets of duplicate detectors, often called ‘dual sensors’. These therefore increase, doubling at maximum, the amount of data traces recorded.
Once recorded, the data are routinely processed in a computer. The term de-ghosting is used to describe the computer-based step to reduce or remove the ghost effect from the data. Many de-ghosting processes today either involve some adaptive summation to extinguish the ghost by polarity difference, or adaptive summation of differently ghosted waveforms to recombine the primary signal present in both, in essence to infill the spectral notches.
The majority of marine seismic data is acquired with conventionally towed streamers which are equipped with only hydrophones and are all typically towed at the same constant depth for any single survey. Fundamentally, it would be advantageous if methods and systems existed that were able to de-ghost the conventionally acquired marine seismic data. Such a method would not depend on special acquisition geometry or special hardware. Instead, in one example, it would be based on only post-acquisition, specially configured computing hardware and/or software, and therefore be backward compatible with all previously acquired and existing seismic data.
Accordingly, it would be desirable to provide systems and methods that avoid the afore-described problems and drawbacks associated with deghosting seismic data.
Thus embodiments described herein provide techniques and systems for deghosting seismic data in a manner which preserves primary wave data without amplifying noise.
According to an embodiment, a method, stored in a memory and executing on a processor, for attenuating ghost reflections in seismic trace data, includes the steps of generating a time-variant deghosting filter based on a cost function using at least one constraint which preserves primary wave data; and applying the time-variant deghosting filter to the seismic trace data for attenuating the ghost reflections.
According to another embodiment, a method, stored in a memory and executing on a processor, for attenuating ghost reflections from seismic trace data, includes the step of applying a time-variant deghosting filter to the seismic trace data.
According to another embodiment, a method, stored in a memory and executing on a processor, for attenuating ghost reflections in seismic trace data, includes the steps of: selecting a cost function and at least one associated constraint for optimizing removal of the ghost reflections wherein the cost function is based on one or more of a plurality of computations comprising a variance computation, an LP Norm computation, a hybrid norm computation and a Wiener entropy computation and the optimizing is based on a predetermined library of predefined filters, generating a time-variant deghosting filter based on the cost function wherein the time-variant deghosting filter is based on a plurality of parameters comprising a seismic source depth, a seismic receiver depth and an estimated sea-surface reflection coefficient; and applying the time-variant deghosting filter to the seismic trace data for attenuating the ghost reflections wherein the time-variant deghosting filter is recalculated for each of a plurality of time intervals associated with the seismic trace data.
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. Some of the following embodiments are discussed, for simplicity, with regard to the terminology and structure of various marine seismic acquisition systems or configurations. However, the embodiments to be discussed next are not limited to these configurations, but may be extended to other arrangements as discussed later.
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 the 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.
According to various embodiments herein, methods and systems for reducing or removing ghost reflections from marine seismic data are described. The methods and embodiments are based on, for example, constructing a time-variant deghosting filter through a constrained computation optimized by a cost function that is selected to retain the primary data while reducing or removing the effects of the ghost(s). Such methods and systems can, for example, be used to reduce the ghosts associated with seismic images without the use of special acquisition geometry or special hardware.
In order to provide some context for the subsequent embodiments related to such deghosting techniques, consider first a seismic data acquisition process and system as will now be described with respect to
One or more source arrays 108 can also be towed by ship 102, or another ship (not shown), for generating seismic waves. Source arrays 108 can be placed either in front of or behind the receivers 112 (one representative receiver per streamer), or both behind and in front of the receivers 112. The seismic waves generated by the source arrays 108 propagate downward, reflect off of, and penetrate the seafloor, wherein the refracted waves eventually are reflected by one or more reflecting structures (not shown in
Looking to
The embodiments below are discussed without specifying the type of seismic receivers used to record the seismic data or specific towing arrangements, as the embodiments can be used with any such receivers or towing arrangements. However those skilled in the art will recognize that, in general for a marine seismic survey, streamers that are towed by one or more vessels and the streamers include the seismic receivers. The streamers may be horizontal or slanted or have a curved profile as illustrated in
The curved streamer 300 of
With this context in mind, a discussion of ghosts and then deghosting according to the embodiments will now be presented. Considering the wave patterns associated with an offshore seismic survey, up-going seismic waves are reflected downwards by the sea surface. At any given depth below the surface, the up and down going waves interfere constructively or destructively as earlier described. Looking to
The embodiments set forth below describe methods and systems for attenuating the effects of this ghost interference through seismic data processing. It should be noted that a simplified (multiple-free) model of one trace of sampled data, x[n], from a seismic acquisition can be represented as:
x[n]=s[n]*g_s[n]*g_r[n]*q[n]*r[n]+w[n] (1)
where n is the temporal sample number, s[n] is the far-field signature of the seismic source, g_s[n] and g_r[n] models the source and receiver ghosts, respectively, r[n] models the earth's reflectivity series, q[n] models the effect of attenuation, and w[n] is noise.
The ghost operators, g_s[n] and g_r[n], are of the form:
g[n]=δ[n]−γδ[n−τ] (2)
where δ is the Dirac delta function, γ is the sea surface reflection coefficient and τ is the propagation delay between the direct wave and the ghost. Assuming a perfectly flat sea surface, the reflection coefficient, γ, is close to 1. However, in practice, γ depends strongly on the sea-state, and varies in both space and time. A typical range of values can be between [0.6, 0.99].
The implication of this practical observation to the embodiments is that, in most cases, the ghost reflection will be weaker than the up-going signal, and there will not be a total destructive interference at the notch frequency, e.g., as shown in
The delay, τ, in Equation (2) is given by the relation:
τ=2z/(v cos(α)) (3)
where z is the depth of the source or receiver, v denotes the speed of sound in water, which is close to 1500 m/s in this purely illustrative example, and α is the take-off angle or the angle of incidence with which the signal impinges upon the sea surface as measured with respect to the vertical. This incidence angle again depends on the reflector dip and the source-receiver offset.
One way to remove ghost effects from acquired seismic data is to calculate, and then apply to the received seismic data, an inverse filter of the ghost models in equation (2). Generally speaking, the inverse filter of the ghost operator g[n] is a filter g−1[n] for which g[n]*g−1[n]=δ[n] where * denotes convolution and δ[n] is the Dirac delta function. In the frequency domain, the ghost operator G(ω) is of the form: G(ω)=1−γ exp(−iωT). The inverse filter of the ghost operator G(ω) is a filter G−1(ω) for which G(ω)G−1(ω)=1, i.e., G−1(ω)=1/G(ω). Thus, for a given pair (γ, τ), the inverse filter g−1[n], or its frequency domain formulation G−1(ω, is uniquely defined.
Stable application of an inverse ghost filter operator for either the source, the receiver, or both, is difficult because of close-to-zero values at the notch (near singularities) and the variability of the parameters. Thus, a straightforward application of the inverse filter generates an undesirable amplification of the noise at the notch frequencies. Thus, according to these embodiments, when solving for the ghost operator is instead transformed into a constrained optimization task, it becomes possible to derive a stable technique that removes the effect of the sea surface in a time-variant fashion, without significantly affecting the noise term in equation (1).
To illustrate this concept according to an embodiment, a technique to attenuate the effect of the receiver side ghost will now be described. In this embodiment, a time-variant filter, f[n], is constructed which removes or attenuates the ghost reflection, without significant effect on the noise, and is found through a constrained optimization of the form:
minf[n]C(f[n],x[n])
“subject to”|f[n]*g_r[n]*
where C(f[n], x[n]) is a suitable cost function,
Looking to
It will be appreciated that various forms of cost functions can be applied to determine the ghost filter in equation (4), both in the time and frequency domain. Examples include:
Variance: C(f[n],x[n])=E{|f[n]*x[n]|2},
LP norm: C(f[n],x[n])=(Σ|f[n]*x[n]|P)1/P,
(where P≧1 is a real number)
Hybrid Norm, and
Wiener entropy (a measure of spectral flatness):
Numerous variants are contemplated. For example, it is possible to choose a number of different optimization criteria. In another example, a combination of different cost functions is used. In yet another example, two or more cost functions are simultaneously combined. In still another example, two or more cost functions are simultaneously combined and weights are applied to each.
The model and optimization in Equation (1) and Equation (4), respectively, describe a method to remove the sea surface effects. In another example, the method is used to correct for the imperfections in source output and attenuation of the signal as it propagates through the earth's interior, in addition to the ghost effects. In another example, the optimization in Equation (4) is solved by a numerical solver, e.g. through linear programming.
In another example, the solution space of the optimization is limited by using a library of K predefined filters, fk[n], k=1, . . . , K, which fulfill the constraint in Equation (4). For each sample along a seismic trace, each filter fk[n] is applied and the cost function is calculated. In one example, a predetermined criterion is used to select the optimal filter. The optimization is then of the form:
minkC(fk[n],x[n]),k=1, . . . , K. (5)
and f[n]=arg min C(fk[n],x[n]).
In one example, the candidate filters for receiver side deghosting are constructed as inverse ghost operators of the form,
fk[n]=gr−1[n] (6)
by systematically varying γ and τ. Here, the indice k refers to a given value of the pair (γ, τ). The candidate filter for the value γ=0 can also be included. The actual implementation of the inverse filters can be varied, and the filters are not necessarily constrained to be of the form in Equation (6).
The foregoing discussion focuses on receiver-side deghosting, however the present invention is not so limited. In one example, the same technique is applied to attenuate the source side ghost. In a further example, the method is extended to attenuate both the source and receiver ghosts. Considering
In another embodiment, the inverse operator and the optimization scheme described above is applied on a trace-by-trace basis, in one example in a shot or, in another example in a CDP-gather. As can be appreciated, under the influence of noise, the optimal filter of an embodiment is found using a number of neighboring traces.
Continuing with an embodiment, as will be appreciated by those skilled in the art, seismic data consist of a signal component and a noise component. However, as a general definition, it can be said that any recorded energy which interferes with the desired signal can be considered as noise. The noise can be classified as background noise (for instance wind, swell, noise from nearby production, or interference from nearby seismic acquisition), source-generated noise (for instance direct and scattered waves or multiples), and instrument noise. Noise can show up as coherent or incoherent energy in seismic gathers. This diversity of noise-types, with different characteristics, makes separation of signal and noise a challenging process. As mentioned previously, in some examples, the performance of the method is dependent on the signal-to-noise ratio for any given trace. Accordingly, in one example, before applying the de-ghosting procedure described herein, it is beneficial to remove as much noise as possible from the data. In practice, in one example, this means that de-ghosting optimally takes place late in the seismic processing sequence, and, in another example, only after various types of background noise and multiples are removed from the data.
In one example of constructing a time-variant de-ghosting filter, the filter parameters are adjusted along the seismic trace to adapt to time-variant changes in the recorded seismic data. The effect of the sea surface reflection(s) on the seismic data can be described by a filter, g[n]. The inverse of filter g[n] can remove the signal associated with these reflections. However, the inverse of filter g[n] is not a stable filter, but tends to boost both high and low frequency noise. Accordingly, it is desirable to find a de-ghosting filter that removes or otherwise attenuates the sea surface reflections without amplifying the noise. Although not a true inverse filter, the constructed filter removes or otherwise attenuates the signal associated with the sea surface reflections without amplifying the noise. The filter will thus “adapt” to the recorded data, in the sense that it will attempt to perform as good a de-ghosting as possible without unduly amplifying the noise. Hence, the filter characteristics will change along the seismic trace, making it a time-variant filter. In practice, in one example, a new set of filter coefficients is found for each time instance, n. The filter coefficients are applied to a time window (subset) of the seismic trace. The time interval is moved along the trace (in time), and a new set of filter coefficients are computed. The filter, being a time-variant filter, is applied “in a time-variant fashion”.
In this way, an optimal de-ghosting filter for a given time instance, n, is found (constructed) through an optimization. In another example, a set of predefined candidate filters are evaluated using a similar cost function as used in the optimization step for each time instance along the trace.
In one example, a constrained optimization is an optimization where the filter output must fulfill a given criteria. In another example, the constraints force the filter to only attenuate the sea surface ghost, not the “primary” data. In another example, the predetermined constraints include a model of the source signature and a model of the ghost filter.
In one example, a cost function gives a measure of how good a given solution is. In one example, the cost function measures of how well a candidate filter, fk[n], performs de-ghosting when applied to the seismic trace, x[n]. It is desirable to make the “cost” as small as possible, hence the minimization. Alternatively, a cost function can be called an optimization function.
In one example, the constraint imposed on the filter, f[n], is that it leaves the “primary” component undistorted in the seismic trace after application of the filter. In one example, this is achieved by setting the constraint criteria to ensure that if the filter is applied to a seismic trace containing a single event, the amplitude of the output is equal to or larger than the amplitude of the “ghost-free” signal, for the duration of the ghost-free part of the signal. This is illustrated in
In one example, retaining the primary data means that all data associated with the “direct wave” is still present in the seismic trace after de-ghosting. Hence, the seismic events, as they would have been recorded without a reflecting sea surface, are preserved in the seismic trace in the de-ghosting process.
In one example, attenuating the effects of the ghost(s) means removing or attenuating the component of the data associated with the sea surface reflection(s). Removing implies that that we have removed all signals associated with the reflections, whereas attenuating implies that there can still be some residual energy from the reflected wave (ghost) left in the seismic trace.
In one example, the seismic data contains a “direct wave” component and a “reflected wave” component, where the reflected wave comes from reflection(s) off the sea surface. Embodiments described herein want to remove data associated with the “reflected wave” from the seismic trace, and leave the data associated with the direct wave, which are called “primary data”.
In one example, “primary data” is defined as data as recorded without a source- and/or receiver-side sea surface reflection. This means that if there were no reflecting sea surface present when acquiring the seismic data trace, this is what the recorded seismic trace would look like. In some examples, a distinction has been made between the reflection from above the source, and the reflection from above the receiver, as the method, in certain examples, can be used to attenuate either or both.
A system for practicing the embodiments is contemplated where a computer system is configured and programmed for seismic processing. In one embodiment, a library of pre-determined filters is stored and/or otherwise accessible by the computer system. In the course of seismic processing, seismic trace records are read trace-by-trace. Each seismic trace is then subjected to testing by applying two or more of the pre-determined filters to a time window portion of the seismic trace. Each output from each of the two or more pre-determined filters is input into a pre-determined cost function, such as, in one example, the L2-norm. The results from the cost function, one result for each filter tested, are compared based on predetermined selection criteria, such as, in one example, selecting the cost function result having the minimum value. Based on the comparison of the results, one of the pre-determined filters is selected. Then, the selected filter is applied to the trace in that time window portion of the seismic trace. The ghost in the time window is thus attenuated. The process is repeated for another time window portion of the seismic trace. On completion of the filtering of the seismic trace, which is now the de-ghosted trace, the newly filtered seismic trace is output. The process is repeated for the next seismic trace.
In a further example, two or more traces are input. The test filter results from each of the two or more traces are combined or otherwise stacked prior to input into the cost function. In another example, the results from each of the two or more traces are input into the cost function. The matrix of results from the cost function is then compared to select which pre-determined filter to apply. In a further example, the two or more traces have a common receiver offset. In another example, the two or more traces have a common source. In another example, the two or more traces have a common mid-point. In another example, the library of pre-determined filters is substituted with filters that are derived from the input seismic data.
Looking now to
Looking now to
Looking now to
Next at step 1004, the method embodiment 1000 generates a time-variant deghosting filter based on the previously selected cost function. It should be noted in the method embodiment that the time-variant deghosting filter is based on a plurality of parameters comprising a seismic source depth, a seismic receiver depth and an estimated sea-surface reflection coefficient. It should further be noted in the embodiment that the determination can be performed in a frequency domain or in a time domain.
Next at step 1006, the method embodiment 1000 applies the time-variant deghosting filter to the seismic trace data to attenuate the ghost reflections. It should be noted in the method embodiment 1000 that the time-variant deghosting filter is recalculated for each of a plurality of time intervals associated with the seismic trace data.
Looking now to
Next in the embodiment node 1100, the deghosting filter component 1104 provides the capability to compute a time-variant deghosting filter. Continuing with the embodiment node 1100, the engine component 1106 provides the capability for deghosting the seismic trace data 1110. It should be noted in the embodiment that the seismic trace data is deghosted by applying the time-variant deghosting filter to the seismic trace data 1110. Next in the embodiment node 1100, the output component provides the capability to output the deghosted seismic trace data for further processing.
The computing device(s) or other network nodes involved in attenuating ghost reflections in seismic trace data as set forth in the above described embodiments may be any type of computing device capable of processing and communicating seismic data associated with a seismic survey. An example of a representative computing system capable of carrying out operations in accordance with these embodiments is illustrated in
Data storage unit 1216 itself can comprise hard disk drive (HDD) 1218 (these can include conventional magnetic storage media, but, as is becoming increasingly more prevalent, can include flash drive-type mass storage devices 1220, among other types), ROM device(s) 1222 (these can include electrically erasable (EE) programmable ROM (EEPROM) devices, ultra-violet erasable PROM devices (UVPROMs), among other types), and random access memory (RAM) devices 1224. Usable with USB port 1210 is flash drive device 1220, and usable with CD/DVD R/W device 1212 are CD/DVD disks 1226 (which can be both read and write-able). Usable with diskette drive device 1214 are floppy diskettes 1228. Each of the memory storage devices, or the memory storage media (1218, 1220, 1222, 1224, 1226, and 1228, among other types), can contain parts or components, or in its entirety, executable software programming code (software) 1230 that can implement part or all of the portions of the method described herein. Further, processor 1208 itself can contain one or different types of memory storage devices (most probably, but not in a limiting manner, RAM memory storage media 1224) that can store all or some of the components of software 1230.
In addition to the above described components, system 1200 also comprises user console 1232, which can include keyboard 1234, display 1236, and mouse 1238. All of these components are known to those of ordinary skill in the art, and this description includes all known and future variants of these types of devices. Display 1236 can be any type of known display or presentation screen, such as liquid crystal displays (LCDs), light emitting diode displays (LEDs), plasma displays, cathode ray tubes (CRTs), among others. User console 1232 can include one or more user interface mechanisms such as a mouse, keyboard, microphone, touch pad, touch screen, voice-recognition system, among other inter-active inter-communicative devices.
User console 1232, and its components if separately provided, interface with server 1202 via server input/output (I/O) interface 1240, which can be an RS232, Ethernet, USB or other type of communications port, or can include all or some of these, and further includes any other type of communications means, presently known or further developed. System 1200 can further include communications satellite/global positioning system (GPS) transceiver device 1242, to which is electrically connected at least one antenna 1244 (according to an embodiment, there would be at least one GPS receive-only antenna, and at least one separate satellite bi-directional communications antenna). System 1200 can access internet 1246, either through a hard wired connection, via I/O interface 1240 directly, or wirelessly via antenna 1244, and transceiver 1242.
Server 1202 can be coupled to other computing devices, such as those that operate or control the equipment of ship 102 of
According to a further embodiment, system 1200, being designed for use in seismic exploration, will interface with one or more sources 1248, 1250 and one or more receivers 1252, 1254. As further previously discussed, sources 1248, 1250 and receivers 1252, 1254 can communicate with server 1202 either through an electrical cable that is part of streamer 1256, 1258, or via a wireless system that can communicate via antenna 1244 and transceiver 1242 (collectively described as communications conduit 1260).
According to further exemplary embodiments, user console 1232 provides a means for personnel to enter commands and configuration into system 1200 (e.g., via a keyboard, buttons, switches, touch screen and/or joy stick). Display device 1236 can be used to show: source/receiver 1256, 1258 position; visual representations of acquired data; source 1248, 1250 and receiver 1252, 1254 status information; survey information; and other information important to the seismic data acquisition process. Source and receiver interface unit 1204 can receive the seismic data from receiver 1252, 1254 though communication conduit 1260 (discussed above). Source and receiver interface unit 1204 can also communicate bi-directionally with sources 1248, 1250 through the communication conduit 1260. Excitation signals, control signals, output signals and status information related to source 1248, 1250 can be exchanged by communication conduit 1260 between system 1200 and source 1248, 1250.
Bus 1206 allows a data pathway for items such as: the transfer and storage of data that originate from either the source sensors or receivers; for processor 1208 to access stored data contained in data storage unit memory 1216; for processor 1208 to send information for visual display to display 1236; or for the user to send commands to system operating programs/software 1230 that might reside in either the processor 1208 or the source and receiver interface unit 1204.
System 1200 can be used to implement the methods described above associated with attenuating ghost reflections in seismic trace data according to an exemplary embodiment. Hardware, firmware, software or a combination thereof may be used to perform the various steps and operations described herein. According to an exemplary embodiment, software 1230 for carrying out the above discussed steps can be stored and distributed on multi-media storage devices such as devices 1218, 1220, 1222, 1224, 1226, and/or 1228 (described above) or other form of media capable of portably storing information (e.g., universal serial bus (USB) flash drive 1220). These storage media may be inserted into, and read by, devices such as the CD-ROM drive 1212, the disk drive 1214, among other types of software storage devices.
Looking now to
It should be noted in the embodiments described herein that these techniques can be applied in either an “offline”, e.g., at a land-based data processing center or an “online” manner, i.e., in near real time while onboard the seismic vessel. For example, predicting a desired seismic quantity at a desired location of a desired depth can occur as the seismic data is recorded onboard the seismic vessel. In this case, it is possible for the prediction to be generated as a measure of the quality of the sampling run.
The disclosed exemplary embodiments provide a server node, and a method for predicting a desired seismic quantity at a desired location of a desired depth. 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. Further, in the detailed description of the exemplary embodiments, numerous specific details are set forth in order to provide a comprehensive understanding of the 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 the embodiments in particular combinations, each feature or element can be used alone without the other features and elements of the embodiments or in various combinations with or without other features and elements disclosed herein. The methods or flow charts provided in the present application may be implemented in a computer program, software, or firmware tangibly embodied in a computer-readable storage medium for execution by a general purpose computer or a processor.
This 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 any devices or systems and performing any incorporated 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 other examples are intended to be within the scope of the claims.
The present application is related to, and claims priority from U.S. Provisional Patent Application No. 61/721,234, filed Nov. 1, 2012, entitled “Method and System for Removing Ghost Reflections from Marine Seismic Data,” to Johan-Fredrik SYNNEVAAG, the disclosure of which is incorporated herein by reference, as well as U.S. Provisional Patent Application No. 61/722,204, entitled “Seismic Target Illumination”, filed on Nov. 4, 2012, the disclosure of which is incorporated here by reference.
| Number | Name | Date | Kind |
|---|---|---|---|
| 5148406 | Brink et al. | Sep 1992 | A |
| 6493636 | DeKok | Dec 2002 | B1 |
| 20100074049 | Kragh et al. | Mar 2010 | A1 |
| 20110292762 | Ozdemir | Dec 2011 | A1 |
| 20120221075 | Bentwich | Aug 2012 | A1 |
| Entry |
|---|
| Kundur et al., ‘Blind Image Deconvolution: An Algorithmic Approach to Practical Image Restoration’, May 1996, UOT Publicaition, pp. 1-42. |
| M. Van Der Baan, et al., “Robust Wavelet Estimation and Blind Deconvolution of Noisy Surface Seismics”, Geophysics, Sep.-Oct. 2008, pp. V37-V46, vol. 73, No. 5. |
| B.J. Posthumus, “Deghosting Using a Twin Streamer Configuration”, Geophysical Prospecting, Apr. 1993, pp. 267-286, vol. 41, No. 3. |
| J. Robertsson, et al., “Rough-Sea Deghosting Using a Single Streamer and a Pressure Gradient Approximation”, Geophysics, Nov.-Dec. 2002, pp. 2005-2011, vol. 67, No. 6. |
| Y. Zhang et al., “A New Bidirectional Deconvolution Method that Overcomes the Minimum Phase Assumption”, SEP Report No. 142, Jun. 2011, pp. 1-13. |
| Z. Zhou, et al., “Analysis of a Broadband Processing Technology Applicable to Conventional Streamer Data”, First Break, Oct. 2012, pp. 77-82, vol. 30. |
| Number | Date | Country | |
|---|---|---|---|
| 20140121981 A1 | May 2014 | US |
| Number | Date | Country | |
|---|---|---|---|
| 61722204 | Nov 2012 | US | |
| 61721234 | Nov 2012 | US |
