WO2010005848A2 - Interpolating seismic data - Google Patents
Interpolating seismic data Download PDFInfo
- Publication number
- WO2010005848A2 WO2010005848A2 PCT/US2009/049426 US2009049426W WO2010005848A2 WO 2010005848 A2 WO2010005848 A2 WO 2010005848A2 US 2009049426 W US2009049426 W US 2009049426W WO 2010005848 A2 WO2010005848 A2 WO 2010005848A2
- Authority
- WO
- WIPO (PCT)
- Prior art keywords
- seismic
- measurements
- seismic measurements
- acquired
- interpolated
- Prior art date
Links
- 238000005259 measurement Methods 0.000 claims abstract description 65
- 238000000034 method Methods 0.000 claims abstract description 56
- 230000008569 process Effects 0.000 claims abstract description 11
- 239000002245 particle Substances 0.000 claims description 24
- 230000033001 locomotion Effects 0.000 claims description 18
- 238000005070 sampling Methods 0.000 claims description 16
- 230000003595 spectral effect Effects 0.000 claims description 9
- 238000001228 spectrum Methods 0.000 description 21
- 238000012545 processing Methods 0.000 description 11
- 239000011159 matrix material Substances 0.000 description 10
- 230000015572 biosynthetic process Effects 0.000 description 8
- 238000005755 formation reaction Methods 0.000 description 8
- 230000001788 irregular Effects 0.000 description 5
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 5
- 238000010586 diagram Methods 0.000 description 4
- 239000004215 Carbon black (E152) Substances 0.000 description 3
- 238000005311 autocorrelation function Methods 0.000 description 3
- 230000008901 benefit Effects 0.000 description 3
- 230000006870 function Effects 0.000 description 3
- 229930195733 hydrocarbon Natural products 0.000 description 3
- 150000002430 hydrocarbons Chemical class 0.000 description 3
- 230000001133 acceleration Effects 0.000 description 2
- 238000013459 approach Methods 0.000 description 2
- 238000004891 communication Methods 0.000 description 2
- 238000006073 displacement reaction Methods 0.000 description 2
- 238000003384 imaging method Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000001629 suppression Effects 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- 238000012952 Resampling Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000001143 conditioned effect Effects 0.000 description 1
- 230000021615 conjugation Effects 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000001151 other effect Effects 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 239000003643 water by type Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S1/00—Beacons or beacon systems transmitting signals having a characteristic or characteristics capable of being detected by non-directional receivers and defining directions, positions, or position lines fixed relatively to the beacon transmitters; Receivers co-operating therewith
- G01S1/02—Beacons or beacon systems transmitting signals having a characteristic or characteristics capable of being detected by non-directional receivers and defining directions, positions, or position lines fixed relatively to the beacon transmitters; Receivers co-operating therewith using radio waves
- G01S1/08—Systems for determining direction or position line
- G01S1/38—Systems for determining direction or position line using comparison of [1] the phase of the envelope of the change of frequency, due to Doppler effect, of the signal transmitted by an antenna moving, or appearing to move, in a cyclic path with [2] the phase of a reference signal, the frequency of this reference signal being synchronised with that of the cyclic movement, or apparent cyclic movement, of the antenna
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/50—Corrections or adjustments related to wave propagation
- G01V2210/57—Trace interpolation or extrapolation, e.g. for virtual receiver; Anti-aliasing for missing receivers
Definitions
- the invention generally relates to interpolating seismic data.
- 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), others to particle motion (e.g., geophones), and industrial surveys may deploy only one type of sensors or both.
- the sensors 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.
- marine surveys Some surveys are known as “marine” surveys because they are conducted in marine environments. However, “marine” surveys may be conducted not only in saltwater environments, but also in fresh and brackish waters.
- a "towed-array” survey an array of seismic sensor-containing streamers and sources is towed behind a survey vessel.
- a technique includes modeling interpolated seismic measurements as a random process that is characterized by seismic measurements acquired at sensor locations and an interpolation error.
- the technique includes determining the interpolated seismic measurements based at least in part on a minimization of the interpolation error.
- an article in another embodiment, includes instructions stored on a computer accessible storage medium that when executed by a processor-based system cause the processor-based system to model interpolated seismic measurements as a random process that is characterized by seismic measurements acquired at sensor locations and an interpolation error.
- the instructions when executed cause the processor-based system to determine the interpolated seismic measurements based at least in part on a minimization of the interpolation error.
- a system in yet another embodiment, includes an interface to receive seismic measurements acquired at non-uniformly spaced sensor locations.
- the system includes a processor to model interpolated seismic measurements at uniformly spaced locations as a random process that is characterized by the seismic measurements acquired at the non-uniformly spaced sensor locations and an interpolation error.
- the processor determines the interpolated seismic measurements based at least in part on a minimization of the interpolation error.
- FIG. 1 is a schematic diagram of a marine acquisition system according to an embodiment of the invention.
- FIGs. 2 and 3 are flow diagrams depicting techniques to interpolate seismic data according to embodiments of the invention.
- FIGs. 4, 5 and 6 are flow diagrams depicting techniques to perform autocorrelation of measured seismic data according to embodiments of the invention.
- FIG. 7 is a schematic diagram of a seismic data processing system according to an embodiment of the invention.
- Fig. 1 depicts an embodiment 10 of a marine seismic data acquisition system in accordance with some embodiments of the invention.
- a survey vessel 20 tows one or more seismic streamers 30 (two exemplary streamers 30 being depicted in Fig. 1) behind the vessel 20.
- the seismic streamers 30 may be several thousand meters long and may contain various support cables (not shown), as well as wiring and/or circuitry (not shown) that may be used to support communication along the streamers 30.
- each streamer 30 includes a primary cable into which is mounted seismic sensors that record seismic signals.
- the seismic sensors are multi-component seismic sensors 58, each of which is capable of detecting a pressure wavef ⁇ eld and at least one component of a particle motion that is associated with acoustic signals that are proximate to the multi-component seismic sensor 58.
- particle motions include one or more components of a particle displacement, one or more components (inline (x), crossline (y) and vertical (z) components (see axes 59, for example)) of a particle velocity and one or more components of a particle acceleration.
- the multi- component seismic sensor 58 may include one or more hydrophones, geophones, particle displacement sensors, particle velocity sensors, accelerometers, pressure gradient sensors, or combinations thereof.
- a particular multi-component seismic sensor 58 may include a hydrophone 55 for measuring pressure and three orthogonally-aligned accelerometers 50 to measure three corresponding orthogonal components of particle velocity and/or acceleration near the seismic sensor 58. It is noted that the multi-component seismic sensor 58 may be implemented as a single device (as depicted in Fig. 1) or may be implemented as a plurality of devices, depending on the particular embodiment of the invention. A particular multi-component seismic sensor 58 may also include pressure gradient sensors 56, which constitute another type of particle motion sensors. Each pressure gradient sensor measures the change in the pressure wavefield at a particular point with respect to a particular direction.
- one of the pressure gradient sensors 56 may acquire seismic data indicative of, at a particular point, the partial derivative of the pressure wavefield with respect to the crossline direction, and another one of the pressure gradient sensors may acquire, a particular point, seismic data indicative of the pressure data with respect to the inline direction.
- the marine seismic data acquisition system 10 includes one or more seismic sources 40 (one exemplary source 40 being depicted in Fig. 1), such as air guns and the like.
- the seismic sources 40 may be coupled to, or towed by, the survey vessel 20.
- the seismic sources 40 may operate independently of the survey vessel 20, in that the sources 40 may be coupled to other vessels or buoys, as just a few examples.
- acoustic signals 42 (an exemplary acoustic signal 42 being depicted in Fig. 1), often referred to as "shots," are produced by the seismic sources 40 and are directed down through a water column 44 into strata 62 and 68 beneath a water bottom surface 24.
- the acoustic signals 42 are reflected from the various subterranean geological formations, such as an exemplary formation 65 that is depicted in Fig. 1.
- the incident acoustic signals 42 that are acquired by the sources 40 produce corresponding reflected acoustic signals, or pressure waves 60, which are sensed by the multi-component seismic sensors 58.
- the pressure waves that are received and sensed by the multi-component seismic sensors 58 include "up going" pressure waves that propagate to the sensors 58 without reflection, as well as “down going” pressure waves that are produced by reflections of the pressure waves 60 from an air- water boundary 31.
- the multi-component seismic sensors 58 generate signals (digital signals, for example), called "traces," which indicate the acquired measurements of the pressure wavefield and particle motion.
- the traces are recorded and may be at least partially processed by a signal processing unit 23 that is deployed on the survey vessel 20, in accordance with some embodiments of the invention.
- a particular multi- component seismic sensor 58 may provide a trace, which corresponds to a measure of a pressure wavefield by its hydrophone 55; and the sensor 58 may provide one or more traces that correspond to one or more components of particle motion, which are measured by its accelerometers 50.
- the goal of the seismic acquisition is to build up an image of a survey area for purposes of identifying subterranean geological formations, such as the exemplary geological formation 65.
- Subsequent analysis of the representation may reveal probable locations of hydrocarbon deposits in subterranean geological formations.
- portions of the analysis of the representation may be performed on the seismic survey vessel 20, such as by the signal processing unit 23.
- the representation may be processed by a seismic data processing system (such as an exemplary seismic data processing system 320 that is depicted in Fig. 7 and is further described below) that may be, for example, located on land or on the vessel 20.
- a seismic data processing system such as an exemplary seismic data processing system 320 that is depicted in Fig. 7 and is further described below
- the down going pressure waves create an interference known as "ghost" in the art.
- the interference between the up going and down going wavefields creates nulls, or notches, in the recorded spectrum. These notches may reduce the useful bandwidth of the spectrum and may limit the possibility of towing the streamers 30 in relatively deep water (water greater than 20 meters (m), for example).
- the technique of decomposing the recorded wavefield into up and down going components is often referred to as wavefield separation, or "deghosting.”
- the particle motion data that is provided by the multi-component seismic sensors 58 allows the recovery of "ghost" free data, which means data that is indicative of the upgoing wavefield.
- the pressure and particle motion measurements typically are acquired at irregularly, or non-uniformly, spaced sensor locations.
- the actual sensor locations do not resemble idealized locations on a uniform grid. Therefore, before the pressure and particle motion measurements that are acquired by the seismic sensors are processed (for such purposes of survey matching, multiple suppression and/or imaging, as examples), the measurements are interpolated to produce interpolated measurements at idealized, uniformly spaced locations.
- Particle motion and pressure wavefield signals that are continuous in the inline (x), crossline (y) and depth (z) directions may be constructed from the sampled pressure and particle motion measurements.
- the Whittaker-Kotel'nikow-Shannon sampling theorem states that any signal f(x) can be reconstructed from its uniformly spaced samples if the sampling interval is less than half the period of the highest frequency component in that signal.
- ⁇ /2 which is known as the Nyquist wavenumber
- the sampling theorem provides the following formula to interpolate any function value from uniformly spaced values f(m/ ⁇ ):
- the sampling theorem provides a way to reconstruct the continuous signal "exactly" from its uniformly spaced samples.
- the signal is sampled at a rate greater than twice the Nyquist rate, i.e., ⁇ .
- the seismic signal may be represented as follows:
- interpolators that are used in seismic data processing are variations of the Yen's interpolator of Type 1.
- interpolators are described in Duijndam, J.J. W., Schonewille, M.A., Hindriks, C.O.H., "Irregular and Sparse Sampling in Exploration Seismology," which is a chapter of L. Zhang, NONUNIFORM SAMPLING: THEORY AND PRACTICE, Kluwer Academic/Plenum Publishers, New York, USA (2001).
- the regularization in Duijndam is formulated as a spectral domain problem, and the spectrum of the signal is estimated by taking a non-uniform Fourier transform of the irregular samples.
- the regularized samples are constructed by using inverse discrete Fourier transformation. It can be shown that this regularization technique is exactly equivalent to Yen's Interpolator of Type 1.
- Zwartjes P.M. and M.D. Sacchi, 2004, "Fourier Reconstruction of Non-Uniformly Sampled, Aliased Data," 74 th Ann. Internal Mtg., Soc. Of Expl. Geophys. (1997-2000).
- Zwartjes' interpolator involves a least squares inversion of the Fourier transformation instead of using an inverse discrete Fourier transform. To this purpose, a cost function is defined, which also involves a non-quadratic penalty term to obtain a parsimonious model.
- Another interpolator is described in Hale, LD.
- Hale's interpolator is based on a more general version of Yen's Interpolator, where a space limited signal assumption is not used. In that case, the uniform samples f(m/ ⁇ ) may be determined by solving a matrix equation similar to Eq. 6. Hale discloses replacing the entries in the inverse matrix by their locally computed approximations.
- the interpolators based on Yen's first theorem usually provide satisfactory results on non-aliased signals that have little high-wavenumber content. However, their performance typically degrades significantly when the interpolated signal has a substantial amount of high wavenumber spectral content.
- Another shortcoming of the interpolators based on Yen's first theorem is that in order to solve Eq. 4, at least as many irregular sampling positions as regular sampling positions are required. Hence, if some seismic traces are dropped out, traces which reside at further locations typically are used to solve the system of equations given by Eq. 4. Usually this degrades the accuracy of the interpolated sample values.
- Yen's first interpolator is exact for infinite length signals, it is an approximation when only a finite extent of the signal is available for interpolation.
- an interpolator for purposes of obtaining an interpolated set of seismic data at uniformly spaced locations, an interpolator is used which models the interpolated seismic data as a stochastic, or random, process. This random process is characterized by the seismic data measurements that have been acquired at non-uniformly-spaced locations and an interpolation error, as described below:
- Ky) ⁇ a m h(x m ) + ⁇ ⁇ m h ⁇ x m ) + w(y) , Eq. 8 m m
- x m represents to the irregular, or non-uniformly, spaced sensor locations
- h(y) represents the interpolated measurement at an arbitrary position "y”
- h(x m ) represents the pressure data measurement acquired at the x m location
- ⁇ m represents a coefficient to be determined, as described below
- “h'(x m )” represents the pressure gradient data (i.e., the particle motion data measurement) acquired at the x m location
- ⁇ m represents a coefficient to be determined, as described below
- w(y) represents the interpolation error.
- the measurements may be two-dimensional (2-D) or three-dimensional (3-D) measurements, depending on the particular embodiment of the invention.
- the statistically optimal (in the least-squares sense) estimate of h(y) may be obtained by minimizing the variance of the interpolation error, as described below:
- Var[w(y)] represents the variance of the interpolation error
- E[] represents the statistical expectation operation, where the random process is assumed for this example to have a mean of zero.
- a minimum norm solution may be obtained by determining the singular value decomposition of II.
- a regularized solution may be obtained by using Tikhonov regularization.
- the statistics of the noise may be incorporated into the interpolator.
- Other variations are contemplated and are within the scope of the appended claims.
- the processor 350 may be coupled to a communication interface 360 for purposes of receiving seismic data that corresponds to the acquired pressure and/or particle motion measurements.
- the processor 350 when executing instructions stored in a memory of the seismic data processing system 320, may receive multi-component data that is acquired by multi- component seismic sensors while in tow.
- the multi-component data may be data that is directly received from the multi-component seismic sensor as the data is being acquired (for the case in which the processor 350 is part of the survey system, such as part of the vessel or streamer) or may be multi-component data that was previously acquired by the seismic sensors while in tow and stored and communicated to the processor 350, which may be in a land-based facility, for example.
- the interface 360 may be a USB serial bus interface, a network interface, a removable media (such as a flash card, CD-ROM, etc.) interface or a magnetic storage interface (IDE or SCSI interfaces, as examples).
- the interface 360 may take on numerous forms, depending on the particular embodiment of the invention.
- the interface 360 may be coupled to a memory 340 of the seismic data processing system 320 and may store, for example, various input and/or output data sets involved with the techniques 120, 140, 170, 180 and/or 200, as indicated by reference numeral 348.
- the memory 340 may store program instructions 344, which when executed by the processor 350, may cause the processor 350 to perform one or more of the techniques that are disclosed herein, such as the techniques 120, 140, 170, 108 and/or 200 and display results obtained via the technique(s) on a display (not shown in Fig. 7) of the system 320, in accordance with some embodiments of the invention.
- the streamers may be in any configuration, such as an over/under configuration, a three-dimensional spread, etc.
- the seismic sensors may be disposed in a sensor cable, such as an ocean bottom cable.
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geology (AREA)
- Environmental & Geological Engineering (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Geophysics (AREA)
- Computer Networks & Wireless Communication (AREA)
- Radar, Positioning & Navigation (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
Description
Claims
Priority Applications (4)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
BRPI0915705A BRPI0915705A2 (en) | 2008-07-05 | 2009-07-01 | method, article and system |
MX2011000192A MX2011000192A (en) | 2008-07-05 | 2009-07-01 | Interpolating seismic data. |
AU2009268798A AU2009268798B2 (en) | 2008-07-05 | 2009-07-01 | Interpolating seismic data |
EP09794992.9A EP2304470A4 (en) | 2008-07-05 | 2009-07-01 | Interpolating seismic data |
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US12/168,125 | 2008-07-05 | ||
US12/168,125 US8861306B2 (en) | 2008-07-05 | 2008-07-05 | Interpolating seismic data |
Publications (2)
Publication Number | Publication Date |
---|---|
WO2010005848A2 true WO2010005848A2 (en) | 2010-01-14 |
WO2010005848A3 WO2010005848A3 (en) | 2010-04-22 |
Family
ID=41464299
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
PCT/US2009/049426 WO2010005848A2 (en) | 2008-07-05 | 2009-07-01 | Interpolating seismic data |
Country Status (6)
Country | Link |
---|---|
US (2) | US8861306B2 (en) |
EP (1) | EP2304470A4 (en) |
AU (1) | AU2009268798B2 (en) |
BR (1) | BRPI0915705A2 (en) |
MX (1) | MX2011000192A (en) |
WO (1) | WO2010005848A2 (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9405028B2 (en) | 2013-02-22 | 2016-08-02 | Ion Geophysical Corporation | Method and apparatus for multi-component datuming |
Families Citing this family (20)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8136383B2 (en) * | 2007-08-28 | 2012-03-20 | Westerngeco L.L.C. | Calibrating an accelerometer |
US8861306B2 (en) * | 2008-07-05 | 2014-10-14 | Westerngeco L.L.C. | Interpolating seismic data |
US20100211320A1 (en) * | 2009-02-13 | 2010-08-19 | Massimiliano Vassallo | Reconstructing a seismic wavefield |
US8554484B2 (en) * | 2009-02-13 | 2013-10-08 | Westerngeco L.L.C. | Reconstructing seismic wavefields |
US8699297B2 (en) * | 2009-02-13 | 2014-04-15 | Westerngeco L.L.C. | Deghosting and reconstructing a seismic wavefield |
US20110103183A1 (en) * | 2009-11-03 | 2011-05-05 | Ahmet Kemal Ozdemir | System and Technique to Increase the Spacing of Particle Motion Sensors on a Seismic Streamer |
US8265875B2 (en) * | 2010-01-29 | 2012-09-11 | Westerngeco L.L.C. | Interpolation of periodic data |
US9010170B2 (en) | 2010-08-16 | 2015-04-21 | Westerngeco L.L.C. | Method and apparatus to test an accelerometer |
US9217805B2 (en) | 2010-10-01 | 2015-12-22 | Westerngeco L.L.C. | Monitoring the quality of particle motion data during a seismic acquisition |
US9043155B2 (en) | 2010-10-07 | 2015-05-26 | Westerngeco L.L.C. | Matching pursuit-based apparatus and technique to construct a seismic signal using a predicted energy distribution |
US8862408B2 (en) | 2011-09-28 | 2014-10-14 | Westerngeco L.L.C. | Determining one or more target positions in an acquisition domain for processing survey data |
US9182512B2 (en) * | 2011-10-26 | 2015-11-10 | Westerngeco L.L.C. | Processing multi-component seismic data |
US20130343156A1 (en) * | 2012-06-25 | 2013-12-26 | Steve Allan Horne | Devices, Systems and Methods for Measuring Borehole Seismic Wavefield Derivatives |
US10353100B2 (en) | 2013-03-21 | 2019-07-16 | Westerngeco L.L.C. | Processing seismic data to attenuate vibration noise |
US20150186563A1 (en) * | 2013-12-30 | 2015-07-02 | Halliburton Energy Services, Inc. | Preconditioning Distinct Subsystem Models in a Subterranean Region Model |
US20150186562A1 (en) * | 2013-12-30 | 2015-07-02 | Halliburton Energy Services, Inc | Preconditioning a Global Model of a Subterranean Region |
US20160327668A1 (en) * | 2014-01-14 | 2016-11-10 | Westerngeco Llc | Interferometry-bsed imaging and inversion |
US20160018543A1 (en) * | 2014-07-21 | 2016-01-21 | Westerngeco L.L.C. | Quality check of compressed data sampling interpolation for seismic information |
US10393901B2 (en) | 2015-08-31 | 2019-08-27 | Pgs Geophysical As | Wavefield interpolation and regularization in imaging of multiple reflection energy |
US10024757B1 (en) | 2017-05-04 | 2018-07-17 | United Launch Alliance, L.L.C. | Non-uniform sampling in bandwidth constrained data acquisition systems |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5978647A (en) | 1997-07-10 | 1999-11-02 | The United States Of America As Represented By The Secretary Of The Navy | Method and apparatus for simulating autocorrelation coefficients in a multipath sonar system |
US6292755B2 (en) | 1999-06-21 | 2001-09-18 | Westerngeco, L.L.C. | 3-D seismic trace extrapolation and interpolation |
US6859734B2 (en) | 2002-10-04 | 2005-02-22 | Paradigm Geophysical Corporation | Method and system for limited frequency seismic imaging |
US7107188B2 (en) | 2003-01-08 | 2006-09-12 | Schlumberger Technology Corporation | Digital pressure derivative method and program storage device |
WO2008008632A2 (en) | 2006-07-12 | 2008-01-17 | Geco Technology B.V. | Time lapse marine seismic surveying employing interpolated multicomponent streamer pressure data |
Family Cites Families (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4821164A (en) * | 1986-07-25 | 1989-04-11 | Stratamodel, Inc. | Process for three-dimensional mathematical modeling of underground geologic volumes |
US4837723A (en) * | 1987-01-30 | 1989-06-06 | Amoco Corporation | Method and apparatus for continuous time-variant digital filtering |
US5012453A (en) * | 1990-04-27 | 1991-04-30 | Katz Lewis J | Inverse vertical seismic profiling while drilling |
US5321613A (en) * | 1992-11-12 | 1994-06-14 | Coleman Research Corporation | Data fusion workstation |
FR2725035A1 (en) * | 1994-09-23 | 1996-03-29 | Geophysique Cie Gle | SEISMIC PROSPECTION PROCESS WITH APPLICATION OF A SELF-DECONVOLUATED PREDICTION ERROR FILTER |
US5617372A (en) * | 1996-08-14 | 1997-04-01 | Western Atlas International, Inc. | Unaliased spatial trace interpolation in the f-k domain |
US5995906A (en) * | 1997-10-03 | 1999-11-30 | Western Atlas International, Inc. | Method for reconciling data at seismic and well-log scales in 3-D earth modeling |
US6510390B1 (en) * | 1999-10-07 | 2003-01-21 | Westerngeco, L.L.C. | 3-D seismic trace extrapolation and interpolation |
FR2834797B1 (en) * | 2002-01-14 | 2004-10-08 | Geophysique Cie Gle | METHOD FOR FILTERING DATA, IN PARTICULAR OF SEISMIC TYPES, BY A RIGGING ANALYSIS |
US6704244B1 (en) * | 2002-10-08 | 2004-03-09 | Pgs Americas, Inc. | Method for acquiring and processing seismic survey data using ocean bottom cables and streamers |
GB2444953B (en) * | 2006-12-19 | 2009-07-22 | Westerngeco Seismic Holdings | Method for obtaining an image of a subsurface by regularizing irregularly sampled seismic data |
US7639564B2 (en) * | 2008-01-11 | 2009-12-29 | Fairfield Industries Incorporated | 3-D TAU-P interpolation |
US7835224B2 (en) * | 2008-03-31 | 2010-11-16 | Westerngeco L.L.C. | Reconstructing low frequency data recordings using a spread of shallow and deep streamers |
US8861306B2 (en) * | 2008-07-05 | 2014-10-14 | Westerngeco L.L.C. | Interpolating seismic data |
-
2008
- 2008-07-05 US US12/168,125 patent/US8861306B2/en active Active
-
2009
- 2009-07-01 BR BRPI0915705A patent/BRPI0915705A2/en not_active IP Right Cessation
- 2009-07-01 AU AU2009268798A patent/AU2009268798B2/en not_active Ceased
- 2009-07-01 WO PCT/US2009/049426 patent/WO2010005848A2/en active Application Filing
- 2009-07-01 MX MX2011000192A patent/MX2011000192A/en active IP Right Grant
- 2009-07-01 EP EP09794992.9A patent/EP2304470A4/en not_active Withdrawn
-
2014
- 2014-09-17 US US14/489,218 patent/US20150003198A1/en not_active Abandoned
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5978647A (en) | 1997-07-10 | 1999-11-02 | The United States Of America As Represented By The Secretary Of The Navy | Method and apparatus for simulating autocorrelation coefficients in a multipath sonar system |
US6292755B2 (en) | 1999-06-21 | 2001-09-18 | Westerngeco, L.L.C. | 3-D seismic trace extrapolation and interpolation |
US6859734B2 (en) | 2002-10-04 | 2005-02-22 | Paradigm Geophysical Corporation | Method and system for limited frequency seismic imaging |
US7107188B2 (en) | 2003-01-08 | 2006-09-12 | Schlumberger Technology Corporation | Digital pressure derivative method and program storage device |
WO2008008632A2 (en) | 2006-07-12 | 2008-01-17 | Geco Technology B.V. | Time lapse marine seismic surveying employing interpolated multicomponent streamer pressure data |
Non-Patent Citations (2)
Title |
---|
P M ZWARTJES ET AL.: "Fourier reconstruction of non-uniformly sampled, aliased seismic data", GEOPHYSICS, vol. 72, no. 1, 2007, pages V21 - V32 |
See also references of EP2304470A4 |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9405028B2 (en) | 2013-02-22 | 2016-08-02 | Ion Geophysical Corporation | Method and apparatus for multi-component datuming |
Also Published As
Publication number | Publication date |
---|---|
EP2304470A4 (en) | 2013-11-13 |
US20100002541A1 (en) | 2010-01-07 |
EP2304470A2 (en) | 2011-04-06 |
AU2009268798A1 (en) | 2010-01-14 |
AU2009268798B2 (en) | 2014-10-02 |
WO2010005848A3 (en) | 2010-04-22 |
BRPI0915705A2 (en) | 2015-10-27 |
US8861306B2 (en) | 2014-10-14 |
MX2011000192A (en) | 2011-04-05 |
US20150003198A1 (en) | 2015-01-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US8861306B2 (en) | Interpolating seismic data | |
US8379481B2 (en) | 3D deghosting of multicomponent or over/under streamer recordings using cross-line wavenumber spectra of hydrophone data | |
US8699297B2 (en) | Deghosting and reconstructing a seismic wavefield | |
US7835224B2 (en) | Reconstructing low frequency data recordings using a spread of shallow and deep streamers | |
US9030910B2 (en) | Jointly interpolating and deghosting seismic data | |
US20100211320A1 (en) | Reconstructing a seismic wavefield | |
US20130182533A1 (en) | Attentuating noise acquired in an energy measurement | |
US20100211322A1 (en) | Interpolating a pressure wavefield along an undersampled direction | |
EP2374026A2 (en) | Using waveform inversion to determine properties of a subsurface medium | |
WO2008157028A2 (en) | Interpolating and deghosting multi-component seismic sensor data | |
WO2010093652A2 (en) | Reconstructing seismic wavefields | |
US10545252B2 (en) | Deghosting and interpolating seismic data |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
121 | Ep: the epo has been informed by wipo that ep was designated in this application |
Ref document number: 09794992 Country of ref document: EP Kind code of ref document: A2 |
|
WWE | Wipo information: entry into national phase |
Ref document number: 2009268798 Country of ref document: AU |
|
WWE | Wipo information: entry into national phase |
Ref document number: MX/A/2011/000192 Country of ref document: MX |
|
NENP | Non-entry into the national phase |
Ref country code: DE |
|
WWE | Wipo information: entry into national phase |
Ref document number: 2009794992 Country of ref document: EP |
|
ENP | Entry into the national phase |
Ref document number: 2009268798 Country of ref document: AU Date of ref document: 20090701 Kind code of ref document: A |
|
WWE | Wipo information: entry into national phase |
Ref document number: 684/CHENP/2011 Country of ref document: IN |
|
ENP | Entry into the national phase |
Ref document number: PI0915705 Country of ref document: BR Kind code of ref document: A2 Effective date: 20110104 |