US10310111B2 - Wave-fields separation for seismic recorders distributed at non-flat recording surfaces - Google Patents
Wave-fields separation for seismic recorders distributed at non-flat recording surfaces Download PDFInfo
- Publication number
- US10310111B2 US10310111B2 US13/670,900 US201213670900A US10310111B2 US 10310111 B2 US10310111 B2 US 10310111B2 US 201213670900 A US201213670900 A US 201213670900A US 10310111 B2 US10310111 B2 US 10310111B2
- Authority
- US
- United States
- Prior art keywords
- going
- datum
- seismic data
- wave fields
- transformed
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active, expires
Links
- 238000000926 separation method Methods 0.000 title description 12
- 238000000034 method Methods 0.000 claims abstract description 50
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 34
- 230000006870 function Effects 0.000 claims description 19
- 239000002245 particle Substances 0.000 claims description 12
- 239000011435 rock Substances 0.000 claims description 8
- 238000013213 extrapolation Methods 0.000 claims description 3
- 230000001131 transforming effect Effects 0.000 claims 2
- 238000004422 calculation algorithm Methods 0.000 description 6
- 230000008569 process Effects 0.000 description 6
- 230000009466 transformation Effects 0.000 description 5
- NJPPVKZQTLUDBO-UHFFFAOYSA-N novaluron Chemical compound C1=C(Cl)C(OC(F)(F)C(OC(F)(F)F)F)=CC=C1NC(=O)NC(=O)C1=C(F)C=CC=C1F NJPPVKZQTLUDBO-UHFFFAOYSA-N 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 230000015654 memory Effects 0.000 description 3
- 230000001133 acceleration Effects 0.000 description 2
- 238000004590 computer program Methods 0.000 description 2
- 238000000354 decomposition reaction Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000007246 mechanism Effects 0.000 description 2
- 230000002123 temporal effect Effects 0.000 description 2
- 230000003044 adaptive effect Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000013500 data storage Methods 0.000 description 1
- 230000003412 degenerative effect Effects 0.000 description 1
- 238000005553 drilling Methods 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 229930195733 hydrocarbon Natural products 0.000 description 1
- 150000002430 hydrocarbons Chemical class 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000010287 polarization Effects 0.000 description 1
- 238000012876 topography Methods 0.000 description 1
- 238000000844 transformation Methods 0.000 description 1
Images
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
-
- 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
- G01V1/364—Seismic filtering
- G01V1/368—Inverse filtering
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/38—Seismology; Seismic or acoustic prospecting or detecting specially adapted for water-covered areas
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/38—Seismology; Seismic or acoustic prospecting or detecting specially adapted for water-covered areas
- G01V1/3843—Deployment of seismic devices, e.g. of streamers
- G01V1/3852—Deployment of seismic devices, e.g. of streamers to the seabed
-
- 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/56—De-ghosting; Reverberation compensation
Definitions
- Embodiments of the subject matter disclosed herein generally relate to methods and systems and, more particularly, to mechanisms and techniques for separating up-going and down-going wave fields from seismic data recorded underwater or under the surface of the earth by a seismic receiver.
- Marine and land seismic data acquisition and processing generate an image of the geophysical structure (subsurface). While this image/profile does not provide a precise location for oil and gas reservoirs, it suggests, to those trained in the field, the presence or absence of oil and/or gas reservoirs. Thus, providing a high-resolution image of the subsurface is an ongoing process for the exploration of natural resources, including, among others, oil and/or gas.
- marine systems for the recording of seismic waves are based on towed streamers or on seafloor-deployed cables or nodes.
- An example of traditional marine system for recording seismic waves at the seafloor is illustrated in FIG. 1 and this system is described in European Patent No. EP 1 217 390, the entire content of which is incorporated herein by reference.
- plural seismic receivers 10 are removably attached to a pedestal 12 together with a memory device 14 . Multiple such receivers are deployed on the bottom 16 of the ocean.
- a source vessel 18 tows a seismic source 20 that is configured to emit seismic waves 22 and 24 . Seismic waves 22 propagate downward, toward the ocean bottom 16 .
- the seismic wave (primary) is recorded (as a trace) by the seismic receiver 10 , while the seismic waves 24 reflected at the water surface 28 are detected by the receivers 10 at a later time. Since the interface between the water and air is well approximated as a quasi-perfect reflector (i.e., the water surface acts as a mirror for the acoustic or seismic waves), the reflected wave 24 travels back toward the receiver 10 .
- This reflected wave is traditionally referred to as a ghost wave because this wave is due to a spurious reflection.
- the ghosts are also recorded by the receivers 10 , but with a different polarization and a time lag relative to the primary wave 22 . As the primary wave 22 moves in an upward direction toward the receiver 10 , this wave is sometimes called an up-going wave-field, and as the ghost 24 moves in a downward direction toward the receiver 10 , this wave is sometimes called a down-going wave-field.
- FIG. 1 also shows the receiver 10 being configured to detach from the pedestal 12 and to rise to the water surface 28 to be retrieved by a collecting boat 30 . Based on the data collected by the receiver 10 , an image of the subsurface is generated by further analyses.
- every arrival of a marine seismic wave at receiver 10 is accompanied by a ghost reflection.
- ghost arrivals trail their primary arrival and are generated when an upward traveling wave is recorded a first time on submerged equipment before being reflected at the surface-air contact.
- Primary and ghost (receiver-side ghost and not the source-side ghost) signals are also commonly referred to as up-going and down-going wave-fields.
- the time delay between an event and its ghost depends entirely upon the depth of the receiver 10 and the wave velocity in water (this can be measured and is considered to be approximately 1500 m/s). It can be only a few milliseconds for towed streamer data (depths of less than 15 meters) or up to hundreds of milliseconds for deep Ocean Bottom Cable (OBC) and Ocean Bottom Node (OBN) acquisitions.
- OBC Ocean Bottom Cable
- OBN Ocean Bottom Node
- Multi-component marine acquisition uses receivers that are capable of measuring at least two different parameters, for example, water pressure (recorded with a hydrophone) and water particle acceleration or velocity (recorded with a geophone or accelerometer).
- water pressure recorded with a hydrophone
- water particle acceleration or velocity recorded with a geophone or accelerometer.
- multi-component marine acquisitions deliver, besides a pressure recording P, at least a vertical particle velocity (or acceleration) component Z.
- Wave-field separation allows the separation of the recorded wave-field into its individual parts: up-going and down-going waves.
- Various techniques are known in the field for wave-field separation, e.g., Amundsen, 1993 , Wavenumber - based filtering of marine point source data , Geophysics; or Ball and Corrigan, 1996 , Dual - sensor summation of noisy ocean - bottom data , SEG Ann.
- the recording surface is a planar surface.
- the ocean bottom is a non-planar acquisition surface.
- the towed-streamer depth may vary along its length, or buried receivers may be deployed at variable depth.
- the planar surface assumption fails, and the collected data may generate spurious effects in the final image unless it is corrected.
- the method includes a step of receiving seismic data (P o , Z o ) recorded in the time-space domain with seismic recorders distributed on a first datum, wherein the first datum is non-flat; a step of establishing a mathematical relation between transformed seismic data (P, Z) and the up-going and down-going wave fields (U, D) on a second planar datum; and a step of solving with an inversion procedure, run on a processor, the mathematical relation to obtain the up-going and down-going wave fields (U, D) for the second datum.
- the second datum is different from the first datum.
- a computing device for separating up-going and down-going wave fields (U, D) in seismic data related to a subsurface of a body of water or to a subsurface of a body of rock.
- the computing device includes an interface configured to receive seismic data (P o , Z o ) recorded in the time-space domain with seismic recorders distributed on a first datum, wherein the first datum is non-flat; and a processor connected to the interface.
- the processor is configured to receive a mathematical relation between transformed seismic data (P, Z) and the up-going and down-going wave fields (U, D) on a second planar datum, and solve with an inversion procedure the mathematical relation to obtain the up-going and down-going wave fields (U, D) for the second datum.
- the second datum is different from the first datum.
- a computer readable medium including computer executable instructions, wherein the instructions, when executed by a processor, implement instructions for separating up-going and down-going wave fields (U, D) in seismic data related to a subsurface of a body of water or to a subsurface of a body of rock.
- the instructions include receiving seismic data (P o , Z o ) recorded in the time-space domain with seismic recorders distributed on a first datum, wherein the first datum is non-flat; establishing a mathematical relation between transformed seismic data (P, Z) and the up-going and down-going wave fields (U, D) on a second planar datum; and solving with an inversion procedure, run on a processor, the mathematical relation to obtain the up-going and down-going wave fields (U, D) for the second datum.
- the second datum is different from the first datum.
- FIG. 1 is a schematic diagram of a conventional seismic data acquisition system having plural seismic receivers provided at the ocean bottom;
- FIG. 2 is a schematic diagram illustrating plural seismic receivers provided on a non-flat datum and a flat datum at which up- and down-going wave-fields are calculated according to an exemplary embodiment
- FIG. 3 is a flowchart illustrating a method for separating up- and down-going wave-fields according to an exemplary embodiment
- FIGS. 4 and 5 are graphs illustrating synthetic P and Z components according to an exemplary embodiment
- FIGS. 6 and 7 are graphs illustrating up- and down-going wave-fields separated according to an exemplary embodiment
- FIG. 8 is a flowchart of another method for separating up- and down-going wave-fields according to an exemplary embodiment.
- FIG. 9 is a schematic diagram of an apparatus configured to run a separation method according to an exemplary embodiment.
- a novel method for separating up- and down-going components includes a first step (i) of determining equations relating the desired separation results on a planar acquisition datum to the available multi-component recordings.
- equations can be formulated in the f-k (frequency wave-number) domain, in the tau-p domain, or other equivalent domains, or a combination of these domains.
- wave-field extrapolation terms to be discussed later
- media properties e.g., sound velocity in water or in rock layers.
- the novel method further includes a step (ii) of inverting the equations from step (i) to find the desired separation results as a function of the available recordings.
- This inversion step can be carried out using a variety of algorithms, for example, analytically or by means of a least-squares inversion. It is noted that the amount of seismic data that is used with the equations and the inversion process require specialized computer software to be implemented on a computing device. As the volume of seismic data that needs to be processed for separating the up-going and down-going components is large, it is impractical, if not impossible, for a human being to do all these calculations in his or her mind.
- FIG. 2 illustrates a non-flat ocean bottom 40 on which plural receivers 42 - 50 have been distributed.
- a body of water 52 is located above the receivers and has a water/air interface 54 .
- a structure 56 is buried below the ocean bottom 40 and it is desired to be imaged with the novel method.
- Another example of a non-flat acquisition surface is the situation when towed-streamer depths vary along their lengths, or when receivers are buried beneath the earth's surface at variable depths.
- each receiver 42 - 50 is configured to record a water pressure P and a particle velocity Z along a z-axis
- equations matrix-to-squares relating (i) the up-going U and down-going D waves at a planar (and also flat) datum 60 and (ii) the recorded and transformed P and Z seismic data on a non-flat datum 62 is given by:
- the original seismic data P o and Z o is a function of the position x i and z i of each receiver 42 - 50 and also a time t at which the data is recorded.
- the y i component is considered to be zero in the example shown in FIG. 2 .
- the method is also applicable to a situation in which the seismic data is recorded with only one-component receivers, i.e., P o or Z o or other component. In this situation, an equation relating the recorded component to the up-going wavefield on a flat datum can be derived and inverted.
- the inversion result might, in this case, contain noise due to the presence of receiver ghost notches. However, this noise can be effectively reduced or removed by the process of stacking or using traditional noise attenuation and signal enhancement techniques.
- the original seismic data P o and Z o is transformed (from the time-space domain), in this example, with a temporal Fast Fourier transformation (FFT) so that the time component t is now an angular frequency component ⁇ .
- FFT temporal Fast Fourier transformation
- the P and Z components in equations (1) and (2) are the temporal FFT of P o and Z o .
- the U and D components in equations (1) and (2) are written in the f-k domain (with f being the frequency corresponding to the angular frequency ⁇ , k being the horizontal wave-number and k z being the vertical wave-number) and these are the up-going and down-going wave-fields desired to be calculated.
- the f-k domain is one possible transformation. Other transformation or transformations may be used.
- the vertical distance ⁇ z i in equations (1) and (2) is the depth difference between the planar datum 60 and the non-planar datum 62 at receiver i.
- the density of the water is represented by ⁇ , and N k is a normalization factor related to the number of receivers at the ocean bottom.
- the terms e j2 ⁇ kx i present in both equations (1) and (2) are related to a spatial inverse FFT that transforms the wave-numbers k to the spatial coordinates of the sensors.
- Wave-field extrapolators for the up-going and down-going wave-fields are also present in equations (1) and (2).
- the wave-field extrapolators can be found in the equations relating the P and Z components to the U and D components irrespective of the transformation domain employed.
- the wave-field extrapolators are given by e ⁇ jk z ⁇ z i .
- the wave-field extrapolators have opposite signs for the U and D components, and they depend from the vertical wave-number and the depth difference between the planar datum 60 and the non-planar datum 62 at receiver i.
- the wavefield extrapolators in this example apply to acoustic propagation with a constant velocity.
- planar datum 60 at which the U and D fields are calculated can be above or below one or more of the receivers 42 - 62 .
- the embodiment shown in FIG. 2 illustrates the planar datum 60 above the receivers. Further, it is possible to have the planar datum 60 to have a flat shape. Furthermore, it is possible that the planar datum 60 is above but close to the receivers 42 - 62 .
- equations (1) and (2) are linear in U and D and, thus, the equations can be inverted using a variety of known algorithms. The details of these algorithms are omitted herein.
- the novel process discussed above may be implemented in a computing device that is provided with dedicated software for separating the up- and down-going components. The computing device is discussed later with regard to FIG. 9 .
- step 300 seismic data (at least two components are recorded, e.g., P and Z) is recorded with corresponding seismic sensors that are provided on the bottom of the ocean.
- the seismic data is transformed in a desired first domain in step 302 .
- the first domain may be the space-frequency domain.
- step 304 equations relating (1) the seismic data transformed in the first domain to (2) up- and down-going wave-fields in a second domain are established.
- the second domain is different from the first domain and may be, for example, the f-k domain.
- the up- and down-going wave-fields correspond to a desired planar acquisition datum, while the transformed seismic data corresponds to a non-flat datum.
- Other domains for the first and second domains are possible.
- step 306 The equations are inverted in step 306 to find the desired separation results as a function of the available recordings. Then, after various processing steps which are known in the art and not repeated herein, an image of the surveyed subsurface is generated in step 308 based on the separated U and/or D.
- the method noted above is now applied to a set of synthetic P and Z data.
- the synthetic P data is illustrated in FIG. 4
- the synthetic Z data is illustrated in FIG. 5 .
- the data is generated as being recorded with a certain offset (distance along X axis) from the source and at a time t (on Y axis) from a non-flat acquisition datum.
- the data shown in FIGS. 4 and 5 is calculated, for example, via acoustic modeling over a half-space and, therefore, it includes only direct arrivals. In this respect, it should be noted that the direct arrival is the most difficult type of event to separate, because its propagation angles are normally wider than for other events.
- FIGS. 6 and 7 show that the down-going wave-field is symmetric, while FIG. 6 shows that the up-going wave-field is complex due to the reflection at the non-flat datum (i.e., ocean bottom).
- the method includes a step 800 of receiving seismic data (P o , Z o ) recorded in the time-space domain with seismic recorders distributed on a first datum, wherein the first datum is non-flat; a step 802 of establishing a mathematical relation between transformed seismic data (P, Z) and the up-going and down-going wave-fields (U, D) on a second plane datum; and a step 804 of solving with an inversion procedure, run on a processor, the mathematical relation to obtain the up-going and down-going wave-fields (U, D) for the second datum.
- the second datum is different from the first datum.
- FIG. 9 An example of a representative computer system capable of carrying out operations in accordance with the exemplary embodiments discussed above is illustrated in FIG. 9 .
- Hardware, firmware, software or a combination thereof may be used to perform the various steps and operations described herein.
- the exemplary computer system 900 suitable for performing the activities described in the exemplary embodiments may include a server 901 .
- a server 901 may include a central processor unit (CPU) 902 coupled to a random access memory (RAM) 904 and to a read-only memory (ROM) 906 .
- the ROM 906 may also be other types of storage media to store programs, such as programmable ROM (PROM), erasable PROM (EPROM), etc.
- the processor 902 may communicate with other internal and external components through input/output (I/O) circuitry 908 and bussing 910 , to provide control signals and the like.
- the processor 902 carries out a variety of functions as are known in the art, as dictated by software and/or firmware instructions.
- the server 901 may also include one or more data storage devices, including hard disk drives 912 , CD-ROM drives 914 , and other hardware capable of reading and/or storing information such as a DVD, etc.
- software for carrying out the above-discussed steps may be stored and distributed on a CD-ROM or DVD 916 , removable media 918 or other form of media capable of portably storing information. These storage media may be inserted into, and read by, devices such as the CD-ROM drive 914 , the drive 912 , etc.
- the server 901 may be coupled to a display 920 , which may be any type of known display or presentation screen, such as LCD or LED displays, plasma displays, cathode ray tubes (CRT), etc.
- a user input interface 922 is provided, including one or more user interface mechanisms such as a mouse, keyboard, microphone, touch pad, touch screen, voice-recognition system, etc.
- the server 901 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 928 .
- GAN global area network
- the exemplary embodiments may be embodied in a wireless communication device, a telecommunication network, as a method or in a computer program product. Accordingly, the exemplary embodiments may take the form of an entirely hardware embodiment or an embodiment combining hardware and software aspects. Further, the exemplary embodiments may take the form of a computer program product stored on a computer-readable storage medium having computer-readable instructions embodied in the medium. Any suitable computer-readable medium may be utilized, including hard disks, CD-ROMs, digital versatile discs (DVD), optical storage devices, or magnetic storage devices such as floppy disk or magnetic tape. Other non-limiting examples of computer-readable media include flash-type memories or other known types of memories.
- the disclosed exemplary embodiments provide an apparatus and a method for seismic data processing. 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 the 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.
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Oceanography (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
Description
Claims (20)
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
FR1161720A FR2984525B1 (en) | 2011-12-15 | 2011-12-15 | WAVE FIELD SEPARATION FOR SEISMIC RECORDINGS DISTRIBUTED ON NON-PLANAR RECORDING SURFACES |
FR1161720 | 2011-12-15 |
Publications (2)
Publication Number | Publication Date |
---|---|
US20130155811A1 US20130155811A1 (en) | 2013-06-20 |
US10310111B2 true US10310111B2 (en) | 2019-06-04 |
Family
ID=47323992
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US13/670,900 Active 2037-12-22 US10310111B2 (en) | 2011-12-15 | 2012-11-07 | Wave-fields separation for seismic recorders distributed at non-flat recording surfaces |
Country Status (5)
Country | Link |
---|---|
US (1) | US10310111B2 (en) |
EP (1) | EP2605047B1 (en) |
BR (1) | BR102012032227A2 (en) |
CA (1) | CA2798794A1 (en) |
FR (1) | FR2984525B1 (en) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US11209564B2 (en) * | 2018-07-18 | 2021-12-28 | Saudi Arabian Oil Company | Deghosting of seismic data through echo- deblending using coincidence filtering |
CN112147691B (en) * | 2019-06-28 | 2024-05-07 | 中国石油化工股份有限公司 | Quick coding ordering-free reference plane correction method and system |
CN114460646B (en) * | 2022-04-13 | 2022-06-28 | 山东省科学院海洋仪器仪表研究所 | Reflected wave travel time inversion method based on wave field excitation approximation |
Citations (28)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1217390A1 (en) | 2000-12-15 | 2002-06-26 | Institut Francais Du Petrole | Method and apparatus for marine seismic exploration, using detectors coupled with the water bottom |
GB2405473A (en) | 2003-08-23 | 2005-03-02 | Westerngeco Ltd | Multiple attenuation in marine seismic survey |
WO2006037815A1 (en) | 2004-10-08 | 2006-04-13 | Compagnie Generale De Geophysique | Improvement to seismic processing for the elimination of multiple reflections |
US7660191B2 (en) * | 2005-07-12 | 2010-02-09 | Westerngeco L.L.C. | Methods and apparatus for acquisition of marine seismic data |
US20100091610A1 (en) | 2008-10-14 | 2010-04-15 | Walter Sollner | Method for imaging a sea-surface reflector from towed dual-sensor streamer data |
US7706211B2 (en) * | 2006-02-06 | 2010-04-27 | Shell Oil Company | Method of determining a seismic velocity profile |
US7751279B2 (en) * | 2006-05-03 | 2010-07-06 | Baker Hughes Incorporated | Sub-salt reflection tomography and imaging by walkaway VSP survey |
US8077544B2 (en) * | 2008-03-28 | 2011-12-13 | Westerngeco L.L.C. | Dual-wavefield multiple attenuation |
US8077543B2 (en) * | 2007-04-17 | 2011-12-13 | Dirk-Jan Van Manen | Mitigation of noise in marine multicomponent seismic data through the relationship between wavefield components at the free surface |
US8275551B2 (en) * | 2006-11-22 | 2012-09-25 | Shell Oil Company | Method of imaging of seismic data involving a virtual source, methods of producing a hydrocarbon fluid, and a computer readable medium |
US8467265B2 (en) * | 2010-01-26 | 2013-06-18 | Westerngeco L.L.C. | Interpolation and deghosting of seismic data acquired in the presence of a rough sea |
US8693281B2 (en) * | 2007-02-24 | 2014-04-08 | Westerngeco L.L.C. | Marine seismic acquisition method and system |
US8937848B2 (en) * | 2010-07-27 | 2015-01-20 | Cggveritas Services Sa | Methods and systems to eliminate undesirable variations in time-lapse seismic surveys |
US8971149B2 (en) * | 2012-05-10 | 2015-03-03 | Cggveritas Services Sa | Method, device and processing algorithm for multiple and noise elimination from marine seismic data |
US9110191B2 (en) * | 2009-03-30 | 2015-08-18 | Westerngeco L.L.C. | Multiple attenuation for ocean-bottom seismic data |
US9134445B2 (en) * | 2009-01-16 | 2015-09-15 | Schlumberger Technology Corporation | Processing seismic data |
US9176249B2 (en) * | 2011-12-21 | 2015-11-03 | Cggveritas Services Sa | Device and method for processing variable depth streamer data |
US9291737B2 (en) * | 2012-05-31 | 2016-03-22 | Pgs Geophysical As | Methods and systems for imaging subterranean formations with primary and multiple reflections |
US9423518B2 (en) * | 2012-02-09 | 2016-08-23 | Pgs Geophysical As | Method for processing dual-sensor streamer data with anti-alias protection |
US9435905B2 (en) * | 2012-04-19 | 2016-09-06 | Cgg Services Sa | Premigration deghosting of seismic data with a bootstrap technique |
US9477001B2 (en) * | 2012-05-11 | 2016-10-25 | Exxonmobil Upstream Research Company | Redatuming seismic data with correct internal multiples |
US9535181B2 (en) * | 2013-04-29 | 2017-01-03 | Cgg Services Sa | Device and method for wave-field reconstruction |
US9684085B2 (en) * | 2013-01-15 | 2017-06-20 | Cgg Services Sas | Wavefield modelling and 4D-binning for seismic surveys from different acquisition datums |
US9784870B2 (en) * | 2013-10-22 | 2017-10-10 | Cgg Services Sas | Demultiple using up/down separation of towed variable-depth streamer data |
US9791581B2 (en) * | 2013-06-07 | 2017-10-17 | Cgg Services Sas | Method and system for simultaneous acquisition of pressure and pressure derivative data with ghost diversity |
US9864084B2 (en) * | 2013-06-07 | 2018-01-09 | Cgg Services Sas | Coherent noise attenuation method |
US10012747B2 (en) * | 2015-03-30 | 2018-07-03 | Cgg Services Sas | Device and method for joint 3-dimensional deghosting of multiple vintages |
US10054704B2 (en) * | 2013-07-01 | 2018-08-21 | Westerngeco L.L.C. | Predicting multiples in survey data |
-
2011
- 2011-12-15 FR FR1161720A patent/FR2984525B1/en not_active Expired - Fee Related
-
2012
- 2012-11-07 US US13/670,900 patent/US10310111B2/en active Active
- 2012-12-13 CA CA2798794A patent/CA2798794A1/en not_active Abandoned
- 2012-12-14 EP EP12197133.7A patent/EP2605047B1/en active Active
- 2012-12-17 BR BR102012032227-7A patent/BR102012032227A2/en not_active IP Right Cessation
Patent Citations (32)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6625083B2 (en) | 2000-12-15 | 2003-09-23 | Institut Francais Du Petrole | Method and device intended for seismic exploration of an underwater subsurface zone using seismic receivers coupled with the water bottom |
EP1217390A1 (en) | 2000-12-15 | 2002-06-26 | Institut Francais Du Petrole | Method and apparatus for marine seismic exploration, using detectors coupled with the water bottom |
GB2405473A (en) | 2003-08-23 | 2005-03-02 | Westerngeco Ltd | Multiple attenuation in marine seismic survey |
US7710821B2 (en) * | 2003-08-23 | 2010-05-04 | Schlumberger Technology Corporation | Multiple attenuation method |
WO2006037815A1 (en) | 2004-10-08 | 2006-04-13 | Compagnie Generale De Geophysique | Improvement to seismic processing for the elimination of multiple reflections |
US7660191B2 (en) * | 2005-07-12 | 2010-02-09 | Westerngeco L.L.C. | Methods and apparatus for acquisition of marine seismic data |
US7706211B2 (en) * | 2006-02-06 | 2010-04-27 | Shell Oil Company | Method of determining a seismic velocity profile |
US7751279B2 (en) * | 2006-05-03 | 2010-07-06 | Baker Hughes Incorporated | Sub-salt reflection tomography and imaging by walkaway VSP survey |
US8275551B2 (en) * | 2006-11-22 | 2012-09-25 | Shell Oil Company | Method of imaging of seismic data involving a virtual source, methods of producing a hydrocarbon fluid, and a computer readable medium |
US8693281B2 (en) * | 2007-02-24 | 2014-04-08 | Westerngeco L.L.C. | Marine seismic acquisition method and system |
US8077543B2 (en) * | 2007-04-17 | 2011-12-13 | Dirk-Jan Van Manen | Mitigation of noise in marine multicomponent seismic data through the relationship between wavefield components at the free surface |
US8077544B2 (en) * | 2008-03-28 | 2011-12-13 | Westerngeco L.L.C. | Dual-wavefield multiple attenuation |
AU2009212893A1 (en) * | 2008-10-14 | 2010-04-29 | Pgs Geophysical As | Method of imaging a sea-surface reflector from towed dual-sensor streamer data |
EP2177933A2 (en) | 2008-10-14 | 2010-04-21 | PGS Geophysical AS | Method for Imaging a Sea-Surface Reflector from Towed Dual-Sensor Streamer Data |
US20100091610A1 (en) | 2008-10-14 | 2010-04-15 | Walter Sollner | Method for imaging a sea-surface reflector from towed dual-sensor streamer data |
US9134445B2 (en) * | 2009-01-16 | 2015-09-15 | Schlumberger Technology Corporation | Processing seismic data |
US9110191B2 (en) * | 2009-03-30 | 2015-08-18 | Westerngeco L.L.C. | Multiple attenuation for ocean-bottom seismic data |
US8467265B2 (en) * | 2010-01-26 | 2013-06-18 | Westerngeco L.L.C. | Interpolation and deghosting of seismic data acquired in the presence of a rough sea |
US8937848B2 (en) * | 2010-07-27 | 2015-01-20 | Cggveritas Services Sa | Methods and systems to eliminate undesirable variations in time-lapse seismic surveys |
US9176249B2 (en) * | 2011-12-21 | 2015-11-03 | Cggveritas Services Sa | Device and method for processing variable depth streamer data |
US9423518B2 (en) * | 2012-02-09 | 2016-08-23 | Pgs Geophysical As | Method for processing dual-sensor streamer data with anti-alias protection |
US9435905B2 (en) * | 2012-04-19 | 2016-09-06 | Cgg Services Sa | Premigration deghosting of seismic data with a bootstrap technique |
US8971149B2 (en) * | 2012-05-10 | 2015-03-03 | Cggveritas Services Sa | Method, device and processing algorithm for multiple and noise elimination from marine seismic data |
US9477001B2 (en) * | 2012-05-11 | 2016-10-25 | Exxonmobil Upstream Research Company | Redatuming seismic data with correct internal multiples |
US9291737B2 (en) * | 2012-05-31 | 2016-03-22 | Pgs Geophysical As | Methods and systems for imaging subterranean formations with primary and multiple reflections |
US9684085B2 (en) * | 2013-01-15 | 2017-06-20 | Cgg Services Sas | Wavefield modelling and 4D-binning for seismic surveys from different acquisition datums |
US9535181B2 (en) * | 2013-04-29 | 2017-01-03 | Cgg Services Sa | Device and method for wave-field reconstruction |
US9791581B2 (en) * | 2013-06-07 | 2017-10-17 | Cgg Services Sas | Method and system for simultaneous acquisition of pressure and pressure derivative data with ghost diversity |
US9864084B2 (en) * | 2013-06-07 | 2018-01-09 | Cgg Services Sas | Coherent noise attenuation method |
US10054704B2 (en) * | 2013-07-01 | 2018-08-21 | Westerngeco L.L.C. | Predicting multiples in survey data |
US9784870B2 (en) * | 2013-10-22 | 2017-10-10 | Cgg Services Sas | Demultiple using up/down separation of towed variable-depth streamer data |
US10012747B2 (en) * | 2015-03-30 | 2018-07-03 | Cgg Services Sas | Device and method for joint 3-dimensional deghosting of multiple vintages |
Non-Patent Citations (21)
Title |
---|
Bancroft-Crewes, John C. "A visualization of the relationship between Kirchhoff migration and seismic inversion.". 2002. * |
C.P.A. Wapenaar et al., "Elastic redatuming of multicomponent seismic data", Geophysical Prospecting 40, presented May-Jun. 1990, pp. 465-482. |
Center for Subsurface Imaging. Chapter 10 Phase Shift Migration. URL: [csim.kaust.edu.sa/files/Short.Courses/BP.2011/PPT/chphasemig.pdf]. 2006. * |
Fortuin, Leonard. "Survey of literature on reflection and scattering of sound waves at the sea surface." The Journal of the Acoustical Society of America 47.5B (1970): 1209-1228. * |
Fred Barr, "Dual sensor OBC technology" The Leading Edge, Houston, Texas, Jan. 1997. |
French Preliminary Search issued in related French Application No. 1161720 dated Sep. 21, 2012. |
Gerard Schuster. Equations for Redatuming Seismic Data by an Asymptotic Inverse. URL: [utam.gg.utah.edu/tomo05/05_ann/PDF/jerry_inter_GRT_single.pdf]. Feb. 6, 2006. * |
Gerard Schuster. Static Corrections to Remove Elevation and Near-Surface Heterogeneities. URL: [http://utam.gg.utah.edu/stanford/node17.html]. 1998. * |
Hocine Tabti et al., "Conventional versus dual sensor streamer deghosting: a case study from a Haltenbanken survey", First Break, vol. 27, Nov. 2009, pp. 137-141. |
IUPAC. Compendium of Chemical Terminology, 2nd ed. (the "Gold Book"). Compiled by A. D. McNaught and A. Wilkinson. Blackwell Scientific Publications, Oxford (1997). XML on-line corrected version: http://goldbook.iupac.org (2006-) created by M. Nic, J. Jirat, B. Kosata; updates compiled by A. Jenkins. ISBN 0-9678550-9-8. doi:10.1351/goldbook. * |
K.M. Schalkwijk et al., "Adaptive decomposition of multi-component ocean-bottom seismic data into downgoing and upgoing P and S waves", Geophysic, vol. 68, No. 3, May-Jun. 2003, pp. 1091-1102. |
Lasse Amundsen, "Wavenumber-Based Filtering of Marine Point Source Data", Geophysics, Sep. 1993, pp. 1104-1107. |
McMechan, George A., and H. W. Chen. "Implicit static corrections in prestack migration of common-source data." Geophysics 55.6 (1990): 757-760. * |
Palacharla, Gopal. "Migration from a non-flat datum via reverse-time extrapolation." Stanford Exploration Project, Report 84. May 9, 2001. pp. 1-50. * |
Robertsson et al., "On the use of multicomponent streamer recordings for reconstruction of pressure wavefields in the crossline direction", Geophysics vol. 73, No. 5, Sep.-Oct. 2008, pp. A45-A49. |
Sabri, M., and Willem Steenaart. "An approach to band-limited signal extrapolation: the extrapolation matrix." Circuits and Systems, IEEE Transactions on 25.2 (1978): 74-78. * |
Stanford Exploration Project. Chapter 2, Why time and space?. URL: [sepwww.stanford.edu/sep/prof/iei2/40_02.pdf]. 2008. * |
Vaughn Ball et al., "Dual-sensor summation of noisy ocean-bottom data", SEG Annual Meeting Nov. 10-15, 1996, Denver, Colorado, pp. 28-31. |
Wapenaar, C. P. A., H. L. H. Cox, and A. J. Berkhout. "Elastic Redatuming of Multicomponent Seismic DATA1." Geophysical prospecting 40.4 (1992): 465-482. * |
Ziolkowski, A. M., D. B. Taylor, and R. G. K. Johnston. "Marine seismic wavefield measurement to remove sea-surface multiples." Geophysical Prospecting 47.6 (1999): 841-870. * |
Ziolkowski, Anton, David Taylor, and Rodney Johnston. "Multiple wavefields: separating incident from scattered, up from down, and primaries from multiples." SEG Technical Program Expanded Abstracts 1998. Society of Exploration Geophysicists, 1998. 1499-1502. * |
Also Published As
Publication number | Publication date |
---|---|
US20130155811A1 (en) | 2013-06-20 |
EP2605047A1 (en) | 2013-06-19 |
EP2605047B1 (en) | 2015-01-21 |
BR102012032227A2 (en) | 2014-03-11 |
FR2984525A1 (en) | 2013-06-21 |
FR2984525B1 (en) | 2014-01-17 |
CA2798794A1 (en) | 2013-06-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US8792296B2 (en) | Method and apparatus for deghosting seismic data | |
US9103941B2 (en) | Device and method for deghosting variable depth streamer data | |
US10345469B2 (en) | Device and method for correcting seismic data for variable air-water interface | |
CA2733699C (en) | Method for separating up and down propagating pressure and vertical velocity fields from pressure and three-axial motion sensors in towed streamers | |
EP2249182B1 (en) | Method for calculation of seismic attributes from seismic signals | |
US9116256B2 (en) | Method and device for wave fields separation in seismic data | |
EP2946232B1 (en) | Wavefield modelling and 4d-binning for seismic surveys with different acquisition datums | |
US9176249B2 (en) | Device and method for processing variable depth streamer data | |
US9448315B2 (en) | Device and method for denoising ocean bottom data | |
US20120026830A1 (en) | Methods and systems to eliminate undesirable variations in time-lapse seismic surveys | |
SG185884A1 (en) | Methods and apparatus for seismic exploration using pressure changes caused by sea-surface variations | |
US9791581B2 (en) | Method and system for simultaneous acquisition of pressure and pressure derivative data with ghost diversity | |
AU2014201497A1 (en) | A system and method for modelling three dimensional shallow water multiples using predictive sea floor reconstruction | |
EP3956696B1 (en) | Attenuation of low-frequency noise in continously recorded wavefields | |
US9201154B2 (en) | Device and method for denoising ocean bottom data | |
US20120213032A1 (en) | Method for pz summation of 3-dimensional wide azimuth receiver gathers and device | |
US10310111B2 (en) | Wave-fields separation for seismic recorders distributed at non-flat recording surfaces |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
AS | Assignment |
Owner name: CGGVERITAS SERVICES SA, FRANCE Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNOR:GRION, SERGIO;REEL/FRAME:029270/0268 Effective date: 20120116 |
|
STPP | Information on status: patent application and granting procedure in general |
Free format text: NOTICE OF ALLOWANCE MAILED -- APPLICATION RECEIVED IN OFFICE OF PUBLICATIONS |
|
AS | Assignment |
Owner name: CGG SERVICES SA, FRANCE Free format text: CHANGE OF NAME;ASSIGNOR:CGGVERITAS SERVICES SA;REEL/FRAME:048391/0168 Effective date: 20130522 |
|
AS | Assignment |
Owner name: CGG SERVICES SAS, FRANCE Free format text: CHANGE OF NAME;ASSIGNOR:CGG SERVICES SA;REEL/FRAME:048405/0506 Effective date: 20161107 |
|
STPP | Information on status: patent application and granting procedure in general |
Free format text: PUBLICATIONS -- ISSUE FEE PAYMENT VERIFIED |
|
STCF | Information on status: patent grant |
Free format text: PATENTED CASE |
|
MAFP | Maintenance fee payment |
Free format text: PAYMENT OF MAINTENANCE FEE, 4TH YEAR, LARGE ENTITY (ORIGINAL EVENT CODE: M1551); ENTITY STATUS OF PATENT OWNER: LARGE ENTITY Year of fee payment: 4 |