[go: up one dir, main page]

WO2016168280A1 - Génération d'un modèle précis de bruit et sa soustraction de données sismiques - Google Patents

Génération d'un modèle précis de bruit et sa soustraction de données sismiques Download PDF

Info

Publication number
WO2016168280A1
WO2016168280A1 PCT/US2016/027261 US2016027261W WO2016168280A1 WO 2016168280 A1 WO2016168280 A1 WO 2016168280A1 US 2016027261 W US2016027261 W US 2016027261W WO 2016168280 A1 WO2016168280 A1 WO 2016168280A1
Authority
WO
WIPO (PCT)
Prior art keywords
noise
signal
measurements
seismic data
seismic
Prior art date
Application number
PCT/US2016/027261
Other languages
English (en)
Inventor
Massimiliano Vassallo
Paul RAS
Ali Ozbek
Nihed EL ALLOUCHE
Yousif Izzeldin Kamil Amin
Original Assignee
Schlumberger Technology Corporation
Schlumberger Canada Limited
Services Petroliers Schlumberger
Geoquest Systems B.V.
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Schlumberger Technology Corporation, Schlumberger Canada Limited, Services Petroliers Schlumberger, Geoquest Systems B.V. filed Critical Schlumberger Technology Corporation
Publication of WO2016168280A1 publication Critical patent/WO2016168280A1/fr

Links

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
    • G01V1/364—Seismic 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/28—Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/32—Transforming one recording into another or one representation into another
    • G01V1/325—Transforming one representation into another
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00—Details of seismic processing or analysis
    • G01V2210/20—Trace signal pre-filtering to select, remove or transform specific events or signal components, i.e. trace-in/trace-out
    • G01V2210/21—Frequency-domain filtering, e.g. band pass
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00—Details of seismic processing or analysis
    • G01V2210/20—Trace signal pre-filtering to select, remove or transform specific events or signal components, i.e. trace-in/trace-out
    • G01V2210/24—Multi-trace filtering
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00—Details of seismic processing or analysis
    • G01V2210/20—Trace signal pre-filtering to select, remove or transform specific events or signal components, i.e. trace-in/trace-out
    • G01V2210/27—Other pre-filtering
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00—Details of seismic processing or analysis
    • G01V2210/30—Noise handling
    • G01V2210/32—Noise reduction
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00—Details of seismic processing or analysis
    • G01V2210/30—Noise handling
    • G01V2210/34—Noise estimation
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00—Details of seismic processing or analysis
    • G01V2210/40—Transforming data representation
    • G01V2210/45—F-x or F-xy domain
    • 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

  • seismic acquisition technologies based on the measurement of data at discrete locations may choose between either oversampling the signal and sampling the noise, or sampling the signal and under-sampling the noise.
  • the former results in a larger number of receivers to acquire the signal of interest, which results in increased cost.
  • the latter may result in a reduced cost acquisition system but may increase a risk of data being compromised by aliased noise, which may be challenging to attenuate during processing processes.
  • vibration noise propagating along the cable is dominant on the accelerometer measurements. This noise propagates along the cable more slowly than seismic events reflected from the subsurface.
  • Other sources of noise also affect accelerometers, also propagating slowly along the cable. Examples of such noise sources are torsional noise, weather related noise, noise associated to positioning equipment, currents, and noise generated by external sources of perturbation, such as barnacles on the cables.
  • Embodiments of the disclosure may provide a method for noise mitigation.
  • the method may include obtaining seismic data including measurements of a seismic wavefield for at least one subsurface volume.
  • the method may also include obtaining at least one component of a spatial gradient of the measurements of the seismic wavefield.
  • the method may include determining a representation of the seismic data based at partially on the measurements of the seismic wavefield and the at least one component of the spatial gradient.
  • the method may include identifying a signal of interest and noise in the representation of the seismic data based at least partially on different characteristics of the signal of interest and the noise in a domain of the representation.
  • the method may include calculating at least one of a signal model or a noise model from the signal of interest and noise identified in the representation of the seismic data.
  • the method may also include determining a noise attenuated signal for the seismic data based at least partially on the at least one of the signal model or the noise model.
  • obtaining the measurements of the seismic wavefield may include receiving acoustic signals collected for the at least one subsurface volume.
  • Obtaining the at least one component of the spatial gradient of the measurement may include determining the at least one component of the spatial gradient from the acoustic signals.
  • Determining the representation of the seismic data may include processing the measurements and the at least one component of the spatial gradient using an analysis process.
  • the analysis process may include at least one of multichannel interpolation by matching pursuit, generalized matching pursuit, an extended generalized matching pursuit, a finite difference multichannel interpolation by matching pursuit, or a greedy algorithm.
  • obtaining the measurements of the seismic wavefield may include receiving acoustic signals collected for the at least one subsurface volume.
  • Obtaining the at least one component of the spatial gradient of the measurements may include receiving the at least one component of the spatial gradient collected for the at least one subsurface volume.
  • Determining the reconstructed wavefield for the seismic data may include processing the measurements and the at least one component of the spatial gradient using an analysis process.
  • the analysis process may include at least one of multichannel interpolation by matching pursuit, generalized matching pursuit, an extended generalized matching pursuit, a finite difference multichannel interpolation by matching pursuit, or a greedy algorithm.
  • determining the noise attenuated signal for the seismic data may include selecting the signal of interest from the representation of the seismic data as the noise attenuated signal.
  • Embodiments of the disclosure may provide a non-transitory computer-readable medium storing instructions.
  • the instructions when executed by one or more processors of a computing system, may cause the computing system to perform a method.
  • the method may include obtaining seismic data including measurements of a seismic wavefield for at least one subsurface volume.
  • the method may also include obtaining at least one component of a spatial gradient of the measurements of the seismic wavefield.
  • the method may include determining a representation of the seismic data based at partially on the measurements of the seismic wavefield and the at least one component of the spatial gradient.
  • the method may include identifying a signal of interest and noise in the representation of the seismic data based at least partially on different characteristics of the signal of interest and the noise in a domain of the representation.
  • the method may include calculating at least one of a signal model or a noise model from the signal of interest and noise identified in the representation of the seismic data.
  • the method may also include determining a noise attenuated signal for the seismic data based at least partially on the at least one of the signal model or the noise model.
  • Embodiments of the disclosure may provide a computing system.
  • the computing system may include one or more processors and a memory system including one or more non-transitory computer-readable media storing instructions.
  • the instructions when executed by one or more processors, may cause the computing system to perform a method.
  • the method may include obtaining seismic data including measurements of a seismic wavefield for at least one subsurface volume.
  • the method may also include obtaining at least one component of a spatial gradient of the measurements of the seismic wavefield.
  • the method may include determining a representation of the seismic data based at partially on the measurements of the seismic wavefield and the at least one component of the spatial gradient.
  • the method may include identifying a signal of interest and noise in the representation of the seismic data based at least partially on different characteristics of the signal of interest and the noise in a domain of the representation.
  • the method may include calculating at least one of a signal model or a noise model from the signal of interest and noise identified in the representation of the seismic data.
  • the method may also include determining a noise attenuated signal for the seismic data based at least partially on the at least one of the signal model or the noise model.
  • Figures 1 A, IB, 1C, ID, 2, and 3 illustrate simplified, schematic views of an oilfield and its operation, according to an embodiment.
  • Figure 4 illustrates a flowchart of an example of a method for noise attenuation, according to an embodiments.
  • Figures 5 A and 5B illustrate plots of traces used as input to an embodiment of a method of the present disclosure, according to an embodiment.
  • Figures 5C and 5D illustrate plots of traces after employing an embodiment of the method of the present disclosure, according to an embodiment.
  • Figures 6, 7, 8, 9, 10, 11, 12, 13, and 14 illustrate flowcharts of different examples of methods for noise attenuation, according to various embodiments.
  • Figures 15A and 15B illustrate another flowchart of an example of another method for noise attenuation, according to an embodiment.
  • Figure 16 illustrates a schematic view of a computing or processor system for performing the method, according to an embodiment.
  • first, second, etc. may be used herein to describe various elements, these elements should not be limited by these terms. These terms are used to distinguish one element from another. For example, a first object or step could be termed a second object or step, and, similarly, a second object or step could be termed a first object or step, without departing from the scope of the present disclosure.
  • the first object or step, and the second object or step are both, objects or steps, respectively, but they are not to be considered the same object or step.
  • Figures 1 A- ID illustrate simplified, schematic views of oilfield 100 having subterranean one or more geological formations 102 containing reservoir 104 therein in accordance with implementations of various technologies and techniques described herein.
  • Figure 1 A illustrates a survey operation being performed by a survey tool, such as seismic truck 106.1, to measure properties of the subterranean formation.
  • the survey operation may be a seismic survey operation for producing sound vibrations.
  • one such sound vibration e.g., sound vibration 1 12 generated by source 1 10
  • a set of sound vibrations may be received by sensors, such as geophone-receivers 1 18, situated on the earth's surface.
  • the data received 120 may be provided as input data to a computer system, for example, a computer 122.1 of a seismic truck 106.1, and responsive to the input data, computer 122.1 generates seismic data output 124.
  • This seismic data output may be stored, transmitted or further processed as desired, for example, by data analysis, data reduction, and the like.
  • Figure IB illustrates a drilling operation being performed by drilling tools 106.2 suspended by rig 128 and advanced into subterranean formations 102 to form wellbore 136.
  • Mud pit 130 may be used to draw drilling mud into the drilling tools via flow line 132 for circulating drilling mud down through the drilling tools, then up wellbore 136 and back to the surface.
  • the drilling mud may be filtered and returned to the mud pit.
  • a circulating system may be used for storing, controlling, or filtering the flowing drilling mud.
  • the drilling tools may be advanced into subterranean formations 102 to reach reservoir 104. Each well may target one or more reservoirs.
  • the drilling tools may be adapted for measuring downhole properties using logging while drilling tools.
  • the logging while drilling tools may also be adapted for taking core sample 133 as shown.
  • Computer facilities may be positioned at various locations about the oilfield 100 (e.g., the surface unit 134) and/or at remote locations.
  • Surface unit 134 may be used to communicate with the drilling tools and/or offsite operations, as well as with other surface or downhole sensors.
  • Surface unit 134 may be capable of communicating with the drilling tools to send commands to the drilling tools, and to receive data therefrom.
  • Surface unit 134 may also collect data generated during the drilling operation and produce data output 135, which may then be stored or transmitted.
  • Sensors (S), such as gauges, may be positioned about oilfield 100 to collect data relating to various oilfield operations as described previously. As shown, sensor (S) may be positioned in one or more locations in the drilling tools and/or at rig 128 to measure drilling parameters, such as weight on bit, torque on bit, pressures, temperatures, flow rates, compositions, rotary speed, and/or other parameters of the field operation. Sensors (S) may also be positioned in one or more locations in the circulating system.
  • Drilling tools 106.2 may include a bottom hole assembly (BHA) (not shown), generally referenced, near the drill bit (e.g. , within several drill collar lengths from the drill bit).
  • BHA bottom hole assembly
  • the bottom hole assembly may include capabilities for measuring, processing, and storing information, as well as communicating with surface unit 134.
  • the bottom hole assembly may further include drill collars for performing various other measurement functions.
  • the bottom hole assembly may include a communication subassembly that communicates with surface unit 134.
  • the communication subassembly may be adapted to send signals to and receive signals from the surface using a communications channel such as mud pulse telemetry, electro-magnetic telemetry, or wired drill pipe communications.
  • the communication subassembly may include, for example, a transmitter that generates a signal, such as an acoustic or electromagnetic signal, which is representative of the measured drilling parameters. It will be appreciated by one of skill in the art that a variety of telemetry systems may be employed, such as wired drill pipe, electromagnetic or other known telemetry systems.
  • the wellbore may be drilled according to a drilling plan that is established prior to drilling.
  • the drilling plan may set forth equipment, pressures, trajectories and/or other parameters that define the drilling process for the wellsite.
  • the drilling operation may then be performed according to the drilling plan. As information is gathered, the drilling operation may deviate from the drilling plan. Additionally, as drilling or other operations are performed, the subsurface conditions may change.
  • the earth model may also be adjusted as information is collected
  • the data gathered by sensors (S) may be collected by surface unit 134 and/or other data collection sources for analysis or other processing.
  • the data collected by sensors (S) may be used alone or in combination with other data.
  • the data may be collected in one or more databases and/or transmitted on or offsite.
  • the data may be historical data, real time data, or combinations thereof.
  • the real time data may be used in real time, or stored for later use.
  • the data may also be combined with historical data or other inputs for further analysis.
  • the data may be stored in separate databases, or combined into a single database.
  • Surface unit 134 may include transceiver 137 to allow communications between surface unit 134 and various portions of the oilfield 100 or other locations.
  • Surface unit 134 may also be provided with or functionally connected to one or more controllers (not shown) for actuating mechanisms at oilfield 100.
  • Surface unit 134 may then send command signals to oilfield 100 in response to data received.
  • Surface unit 134 may receive commands via transceiver 137 or may itself execute commands to the controller.
  • a processor may be provided to analyze the data (locally or remotely), make the decisions and/or actuate the controller. In this manner, oilfield 100 may be selectively adjusted based on the data collected. This technique may be used to optimize (or improve) portions of the field operation, such as controlling drilling, weight on bit, pump rates, or other parameters. These adjustments may be made automatically based on computer protocol, and/or manually by an operator. In some cases, well plans may be adjusted to select optimum (or improved) operating conditions, or to avoid problems.
  • Figure 1C illustrates a wireline operation being performed by wireline tool 106.3 suspended by rig 128 and into wellbore 136 of Figure IB.
  • Wireline tool 106.3 may be adapted for deployment into wellbore 136 for generating well logs, performing downhole tests and/or collecting samples.
  • Wireline tool 106.3 may be used to provide another method and apparatus for performing a seismic survey operation.
  • Wireline tool 106.3 may, for example, have an explosive, radioactive, electrical, or acoustic energy source 144 that sends and/or receives electrical signals to surrounding subterranean formations 102 and fluids therein.
  • Wireline tool 106.3 may be operatively connected to, for example, geophones 1 18 and a computer 122.1 of a seismic truck 106.1 of Figure 1A. Wireline tool 106.3 may also provide data to surface unit 134. Surface unit 134 may collect data generated during the wireline operation and may produce data output 135 that may be stored or transmitted. Wireline tool 106.3 may be positioned at various depths in the wellbore 136 to provide a survey or other information relating to the subterranean formation 102.
  • Sensors (S) such as gauges, may be positioned about oilfield 100 to collect data relating to various field operations as described previously. As shown, sensor S may be positioned in wireline tool 106.3 to measure downhole parameters which relate to, for example porosity, permeability, fluid composition and/or other parameters of the field operation.
  • Figure ID illustrates a production operation being performed by production tool 106.4 deployed from a production unit or Christmas tree 129 and into completed wellbore 136 for drawing fluid from the downhole reservoirs into surface facilities 142.
  • the fluid flows from reservoir 104 through perforations in the casing (not shown) and into production tool 106.4 in wellbore 136 and to surface facilities 142 via gathering network 146.
  • Sensors (S), such as gauges, may be positioned about oilfield 100 to collect data relating to various field operations as described previously. As shown, the sensor (S) may be positioned in production tool 106.4 or associated equipment, such as Christmas tree 129, gathering network 146, surface facility 142, and/or the production facility, to measure fluid parameters, such as fluid composition, flow rates, pressures, temperatures, and/or other parameters of the production operation.
  • production tool 106.4 or associated equipment, such as Christmas tree 129, gathering network 146, surface facility 142, and/or the production facility, to measure fluid parameters, such as fluid composition, flow rates, pressures, temperatures, and/or other parameters of the production operation.
  • Production may also include injection wells for added recovery.
  • One or more gathering facilities may be operatively connected to one or more of the wellsites for selectively collecting downhole fluids from the wellsite(s).
  • Figures 1 A-1D illustrate tools used to measure properties of an oilfield
  • the tools may be used in connection with non-oilfield operations, such as gas fields, mines, aquifers, storage or other subterranean facilities.
  • non-oilfield operations such as gas fields, mines, aquifers, storage or other subterranean facilities.
  • various measurement tools capable of sensing parameters, such as seismic two-way travel time, density, resistivity, production rate, etc., of the subterranean formation and/or its geological formations may be used.
  • Various sensors (S) may be located at various positions along the wellbore and/or the monitoring tools to collect and/or monitor the desired data. Other sources of data may also be provided from offsite locations.
  • Figures 1 A-1D are intended to provide a brief description of an example of a field usable with oilfield application frameworks. Part of, or the entirety, of oilfield 100 may be on land, water and/or sea, as further discussed below. Also, while a single field measured at a single location is depicted, oilfield applications may be utilized with any combination of one or more oilfields, one or more processing facilities and one or more wellsites.
  • Figure 2 illustrates another schematic view, partially in cross section of oilfield 200 having data acquisition tools 202.1, 202.2, 202.3 and 202.4 positioned at various locations along oilfield 200 for collecting data of subterranean formation 204 in accordance with implementations of various technologies and techniques described herein.
  • Data acquisition tools 202.1-202.4 may be the same as data acquisition tools 106.1-106.4 of Figures 1A-1D, respectively, or others not depicted. As shown, data acquisition tools 202.1-202.4 may generate data plots or measurements 208.1-208.4, respectively. These data plots are depicted along oilfield 200 to demonstrate the data generated by the various operations.
  • Data plots 208.1-208.3 are examples of static data plots that may be generated by data acquisition tools 202.1-202.3, respectively; however, it should be understood that data plots 208.1- 208.3 may also be data plots that are updated in real time. These measurements may be analyzed to better define the properties of the formation(s) and/or determine the accuracy of the measurements and/or for checking for errors. The plots of each of the respective measurements may be aligned and scaled for comparison and verification of the properties.
  • Static data plot 208.1 may be a seismic two-way response over a period of time.
  • Static plot 208.2 may be core sample data measured from a core sample of the formation 204.
  • the core sample may be used to provide data, such as a graph of the density, porosity, permeability, or some other physical property of the core sample over the length of the core. Tests for density and viscosity may be performed on the fluids in the core at varying pressures and temperatures.
  • Static data plot 208.3 may be a logging trace that typically provides a resistivity or other measurement of the formation at various depths.
  • a production decline curve or graph 208.4 may be a dynamic data plot of the fluid flow rate over time.
  • the production decline curve may provide the production rate as a function of time.
  • measurements may be taken of fluid properties, such as flow rates, pressures, composition, etc.
  • the subterranean structure 204 may have a plurality of geological formations 206.1- 206.4. As shown, this structure may have several formations or layers, including a shale layer 206.1, a carbonate layer 206.2, a shale layer 206.3 and a sand layer 206.4. A fault 207 may extend through the shale layer 206.1 and the carbonate layer 206.2.
  • the static data acquisition tools may be adapted to take measurements and detect characteristics of the formations.
  • oilfield 200 may contain a variety of geological structures and/or formations, sometimes being complex. In some locations, typically below the water line, fluid may occupy pore spaces of the formations.
  • Each of the measurement devices may be used to measure properties of the formations and/or its geological features. While each acquisition tool is shown as being in specific locations in oilfield 200, it will be appreciated that one or more types of measurement may be taken at one or more locations across one or more fields or other locations for comparison and/or analysis.
  • the data collected from various sources may then be processed and/or evaluated.
  • seismic data displayed in static data plot 208.1 from data acquisition tool 202.1 may be used by a geophysicist to determine characteristics of the subterranean formations and features.
  • the seismic data may be processed and analyzed to aid in the evaluation of the seismic data, as described herein.
  • the core data shown in static plot 208.2 and/or log data from well log 208.3 may be used by a geologist to determine various characteristics of the subterranean formation.
  • the production data from graph 208.4 may be used by the reservoir engineer to determine fluid flow reservoir characteristics.
  • the data analyzed by the geologist, geophysicist and the reservoir engineer may be analyzed using modeling techniques.
  • FIG. 3 illustrates a side view of a marine-based survey 360 of a subterranean geological subsurface 362 in accordance with one or more implementations of various techniques described herein.
  • Subsurface 362 may include seafloor surface 364.
  • Seismic sources 366 may include marine sources such as vibroseis or airguns, which may propagate seismic waves 368 (e.g., energy signals) into the Earth over an extended period of time or at a nearly instantaneous energy provided by impulsive sources.
  • the seismic waves may be propagated by marine sources as a frequency sweep signal.
  • marine sources of the vibroseis type may initially emit a seismic wave at a low frequency (e.g., 5 Hz) and increase the seismic wave to a high frequency (e.g., 80-90Hz) over time.
  • the component(s) of the seismic waves 368 may be reflected and converted by seafloor surface 364 (i.e., reflector), and seismic wave reflections 370 may be received by a plurality of seismic receivers 372.
  • Seismic receivers 372 may be disposed on a plurality of streamers (i.e., streamer array 374).
  • the seismic receivers 372 may generate electrical signals representative of the received seismic wave reflections 370.
  • the electrical signals may be embedded with information regarding the subsurface 362 and captured as a record of seismic data.
  • each streamer may include streamer steering devices such as a bird, a deflector, a tail buoy and the like, which are not illustrated in this application.
  • the streamer steering devices may be used to control the position of the streamers in accordance with the techniques described herein.
  • seismic wave reflections 370 may travel upward and reach the water/air interface at the water surface 376, a portion of reflections 370 may then reflect downward again (i.e., sea-surface ghost waves 378) and be received by the plurality of seismic receivers 372.
  • the sea-surface ghost waves 378 may be referred to as surface multiples.
  • the point on the water surface 376 at which the wave may be reflected downward is generally referred to as the downward reflection point.
  • the electrical signals may be transmitted to a vessel 380 via transmission cables, wireless communication or the like.
  • the vessel 380 may then transmit the electrical signals to a data processing center.
  • the vessel 380 may include an onboard computer capable of processing the electrical signals (i.e., seismic data).
  • seismic data i.e., seismic data
  • surveys may be of formations deep beneath the surface.
  • the formations may typically include multiple reflectors, some of which may include dipping events, and may generate multiple reflections (including wave conversion) for receipt by the seismic receivers 372.
  • the seismic data may be processed to generate a seismic image of the subsurface 362.
  • Marine seismic acquisition systems tow each streamer in streamer array 374 at the same depth (e.g., 5-10m).
  • marine based survey 360 may tow the streamers in streamer array 374 at different depths such that seismic data may be acquired and processed in a manner that avoids the effects of destructive interference due to sea-surface ghost waves.
  • marine- based survey 360 of Figure 3 illustrates eight streamers towed by vessel 380 at eight different depths. The depth of each streamer may be controlled and maintained using the birds disposed on each streamer.
  • Variable-depth single-streamer marine streamer deghosting decomposes the upgoing wavefield in terms of its wavenumber along a single horizontal coordinate. This may be done on a frequency-by-frequency basis. A linear system may then be built that models the ghost operator, assuming that there is no propagation orthogonal to this coordinate axis. The upgoing wavefield may then be estimated by inverting this linear system. If the out-of-plane propagation is large, the inverse of the linear system may become unstable, and damping may be applied.
  • seismic data including seismic measurements of at least one subsurface volume may be processed to remove or reduce the noise contained in the seismic data.
  • seismic data including seismic measurements of the at least one subsurface volume may be analyzed using a multi -component interpolation technology, or multi-component spectral analysis, to create an accurate model of the noise present on seismic measurements.
  • the estimated model of noise may be accurate enough to allow its direct subtraction from the seismic measurement without an adaptive filter, although such an adaptive filter may be used.
  • the multi-component interpolation technology may be used for attenuation of noise that has a slower propagation velocity than the signal, in the case it is aliased during the sampling process, which may make it difficult to attenuate from the measured data.
  • Noise such as "ground roll” in land acquisition or "vibration noise” in multi-sensor marine measurements may be targets of the present methods.
  • the methods include data sampled by receivers that are able to measure either a signal and its directional gradient along the direction of propagation of the noise, or groups of at least two receivers deployed at close distance, to allow a usable estimate of the gradient with finite different methods (or others).
  • the process of analyzing the seismic data may utilize dynamic range effects, when the noise is much stronger than the signal of interest. In some embodiments, this may enable an effective and signal safe attenuation of the noise that is outside of the theoretical bandwidth of signal.
  • the noise targeted by the analysis can be either coherent or incoherent without particular distinction.
  • the analysis may allow separating the near surface noise from the deeper seismic reflections. Incoherent ambient noise may also be reduced when using the analysis process.
  • the process of analyzing the seismic data may be applied to seismic data sampled by receivers able to measure either, a signal and its directional gradient along the direction of propagation of the noise, or groups of at least two receivers deployed at close distance, to allow a usable estimate of the gradient with finite difference methods.
  • Figure 4 illustrates an example of a method 400 of removing noise from seismic data, according to various embodiments.
  • the illustrated processes of the method 400 are examples and that any of the illustrated processes can be removed, additional processes may be added, and the order of the illustrated processes may be changed.
  • any of the illustrated processes of the method 400 may be performed by one or more computer systems, for example, the computer systems illustrated in Figures 1 A-ID, Figure 2, and Figure 3, described above, and/or the computer systems described below in Figures 15A and 15B.
  • seismic data including measurements of at least one subsurface volume may be received.
  • the seismic data may be received by one or more computer systems, for examples, any of the computer systems described herein.
  • the seismic data may include measurements measured directly by seismic survey tools.
  • the measurements measured directly by seismic survey tools may include direct sensor measurements of acoustic return signal, direct particle motion measurements, direct measurement of sensor locations, and the like.
  • the seismic data may include measurements generated by processing the measurements measured directly by seismic survey tools.
  • the direct sensor measurements may be processed to generate an estimated measurement for a group of two or more sensors.
  • the direct particle motion measurements may be processed to covert the measurements to a spatial gradient.
  • measurement of a group of two or more sensors may be processed to estimate the spatial gradient of the particular motion measurement, e.g., a finite difference operator or other difference operators may be used to estimate the gradients.
  • Figure 5 A illustrates a plot showing seismic data, i.e., input traces sampled at 50m spacing
  • Figure 5B illustrates the representation in the frequency- wavenumber (FK) domain of a portion of these data.
  • the vertical axis is time, measured in seconds, and the horizontal is space, in meters.
  • the vertical axis represents frequency, in Hz, and the horizontal the wavenumber, in 1/m units.
  • the seismic data was obtained by averaging measurements of couples of receivers at short distance. In some embodiments, these may be direct measurements of single sensors, eventually processed. The severe alias of the near surface noise is recognizable, especially in the FK plot at low frequencies.
  • the seismic data may include data that was collected using any type and/or configuration of seismic acquisition system, for example, survey tools as described above.
  • the seismic data may be collected by any kind of seismic acquisition system that allows one to measure or to estimate a spatial gradient of a measured wavefield, either at the source or receiver side, or combination of both.
  • the seismic survey tool may include single or multi-sensor marine acquisition systems that include receivers spaced along the cables.
  • the marine acquisition systems may also include a coarse set of accelerometers along the cable with a layout, for example, as a sparse series of couples.
  • nodes may contain a number of few seismic receivers in a small space, allowing the estimate of the gradient of the wavefield of interest.
  • the noise to be separated from the reflection signal may include near surface events, as well as incoherent ambient noise.
  • spatial gradient may be estimated along the cable direction, by combining couples or groups of consecutive sensors. This may allow the attenuation of noise propagating along the inline direction, such as near surface noise.
  • the marine acquisition systems may include a design for seismic streamers, containing a smaller number of receivers that offer another noise attenuation strategy, not requiring dense sampling along the cable.
  • the spatial gradient may be computed across cables, at low frequencies where the wavefield would not be aliased.
  • the seismic survey tool may include a land acquisition system that includes a single or a group of two or more receivers deployed at sampling positions.
  • the seismic data may include measurements from traditional sensors that positioned in a way that facilitates groups of seismic traces to be processed so that an estimate of the gradient can be computed by applying finite difference operators.
  • the positioning may include couples of receivers at short distance, each couple being at wide distance from the next couple of receivers.
  • the land acquisition systems may include sticks or sensor technologies that measure directly a spatial gradient of the wavefield.
  • the seismic survey tool may include ocean bottom nodes (OBN) acquisition systems.
  • OBN acquisition systems nodes may contain one or more seismic receivers in a small space, allowing the estimate of the gradient of the wavefield of interest.
  • the noise to be separated from the reflection signal may include near surface events, as well as incoherent ambient noise.
  • spatial gradient may be estimated along the cable direction, by combining couples or groups of consecutive sensors. This would allow the attenuation of noise propagating along the inline direction, such as near surface noise. This could also be applied to multi-measurements streamer to estimate the inline gradient of the acceleration measurements, hence allowing attenuation of the slow and strong vibration and torsional noise from the acceleration measurements.
  • the OBN acquisition systems may include a design for seismic streamers, containing a smaller number of receivers that offer a noise attenuation strategy, not requiring dense sampling along the cable.
  • spatial gradient may be computed across cables, at low frequencies where the wavefield would not be aliased.
  • any of the seismic acquisition systems may include seismic source devices that allow spatial gradients to be measured or estimated on the source side.
  • marine acquisition systems or OBN acquisition systems may include marine vibrators that would enable the generation of a source gradient wavefield.
  • the source device may allow the estimate the source gradient. Once the source gradient is available, the seismic data can be used to further attenuate noise outside of the theoretical bandwidth of the signal.
  • the combination of source and receiver gradients, either measured or estimated may allow the use of the seismic data in a variety of domains, such as, the offset-azimuth-midpoint.
  • a reconstructed wavefield may be determined for the seismic data based at least partially on the measurements and spatial gradient data associated with the seismic data.
  • the seismic data may be analyzed using multi-channel interpolation and/or spatial analysis techniques. The analysis may include using the measured wavefield, an estimate of the measured wavefield, and the spatial gradient (estimated or measured), or combination thereof to determine the wavefield components, e.g., a signal of interest and noise signal, of the seismic data.
  • multi-channel interpolation may allow the reconstruction of the measured signal on a fine grid of receivers even when its original spatial sampling may be too coarse to allow such reconstruction on the fine grid in theory.
  • the multi-channel interpolation techniques may process the gradient measurements, which provide additional information.
  • MIMAP Multichannel Interpolation by Matching Pursuit
  • MIMAP is an analysis technique based at least partially on spectral estimation of the seismic data.
  • the spectral estimate may enable the computation of the interpolated wavefield.
  • MIMAP can include an analysis technique that uses pressure and gradients to interpolate the pressure wavefield.
  • MIMAP may allow reconstruction of wavefield affected by severe spatial aliasing.
  • MIMAP may estimate a reconstructed wavefield by iteratively matching the input measurements until these are fully explained by a set of basis functions.
  • MFMAP may match iteratively the input measurements by trying to estimate their spatial spectrum and determine iteratively which basis functions are likely components of the wavefield.
  • MFMAP may achieve reconstruction by estimating the spectrum of the data iteratively until a stop condition, such as the norm of the error being below a predetermined value, is reached.
  • MFMAP may also be utilized to determine the components of the wavefield where spatial gradient is not directly measured, but an estimated spatial gradient is available or determined.
  • an estimated spatial gradient may be generated by a finite difference operator on the measurements of a group of receivers. MFMAP may treat the estimated spatial gradient as if it was an actual spatial gradient for the analysis.
  • GMP Generalized Matching Pursuit
  • GMP may be used to analyze the seismic data and determine the components of the wavefield.
  • GMP is a generalization of MFMAP.
  • GMP may process the pressure wavefield as well as horizontal and/or vertical components of the gradient of the pressure measured by multi-sensor seismic acquisition system.
  • GMP may achieve joint interpolation and deghosting of the pressure wavefield.
  • GMP may estimate a reconstructed wavefield by matching iteratively the input measurements until these are fully explained by a set of basis functions. For example, the GMP may match iteratively the input measurements by trying to estimate their spatial spectrum and determine iteratively which basis functions are likely components of the wavefield.
  • GMP may achieve reconstruction by estimating the spectrum of the data iteratively until a stop condition, such as the norm of the error being below a predetermined value.
  • GMP may also be utilized to determine the components of the wavefield where spatial gradient is not directly measured, but an estimated spatial gradient is available or determined.
  • an estimated spatial gradient may be generated by a finite difference operator on the measurements of a group of receivers. GMP may treat the estimated spatial gradient as if it was an actual spatial gradient for the analysis.
  • E-GMP Extended Generalized Matching Pursuit
  • E-GMP may be used to analyze the seismic data and determine the components of the wavefield.
  • E-GMP may estimate a reconstructed wavefield by matching iteratively the input measurements until these are fully explained by a set of basis functions. For example, the E-GMP may match iteratively the input measurements by trying to estimate their spatial spectrum and determine iteratively which basis functions are likely components of the wavefield.
  • E-GMP may achieve reconstruction by estimating the spectrum of the data iteratively until a stop condition is reached.
  • E-GMP may process a more generic set of data, provided that the relation between the input measurements and the desired output is known.
  • E-GMP may determine the components of the wavefield where spatial gradient is not directly measured, but an estimated spatial gradient is available or determined.
  • an estimated spatial gradient may be generated by a finite difference operator on the measurements of a group of receivers.
  • E-GMP may consider the actual model used to generate the estimated gradient, and may matches the data respecting such model.
  • E-GMP may use this estimated gradient for interpolation purposes without suffering of the approximation error in estimating the gradient.
  • E-GMP may use the finite different operator to model the relation between a signal of interest and its estimated spatial gradient.
  • the use of an actual operator used to estimate the spatial gradient may provide an improved match than the use of the ideal gradient operator that would be used within GMP.
  • finite difference MFMAP may be used to analyze the seismic data and determine the components of the wavefield.
  • FD-MFMAP is a particular realization of E-GMP that generalizes MFMAP, i.e., extends MFMAP to the finite difference gradient.
  • an estimated spatial gradient may be generated by a finite difference operator on the measurements of a group of receivers.
  • FD-MFMAP may treat the estimated spatial gradient as if it were an actual spatial gradient for the analysis.
  • MFMAP, GMP, E-GMP, and FD-MEVIAP are techniques based on spectral estimation of the input data: in these cases, the spectral estimate supports the computation of the interpolated wavefield. In these examples, the spectral estimation performed by these (and other) techniques may produce a signal of interest, without any interpolated data being generated.
  • non-matching point algorithms may be utilized to determine wavefield components.
  • greedy algorithms may be utilized to determine the components of the wavefield.
  • a greedy algorithm may be an algorithm that follows the problem solving heuristic of making the locally optimal choice at each process with the hope of finding a global optimum.
  • a greedy heuristic may yield locally optimal solutions that approximate a global optimal solution in a reasonable time.
  • a greedy strategy for the traveling salesman problem may be the following heuristic: "At each process visit an unvisited city nearest to the current city.” This heuristic may not find an exhaustive solution, but terminate in a reasonable number of iterations; finding a globally optimal solution typically may include a large number of iterations, and may thus be impractical. In mathematical optimization, greedy algorithms solve combinatorial problems having the properties of matroids.
  • Greedy algorithms may have five components: a candidate set, from which a solution is created; a selection function, which chooses the candidate to be added to the solution; a feasibility function, that is used to determine if a candidate can be used to contribute to a solution; an objective function, which assigns a value to a solution, or a partial solution, and a solution function, which will indicate when a complete solution is discovered.
  • MTMAP, GMP, E-GMP, and FD-MFMAP may utilize the measurements or estimates of a measured seismic signal and its spatial gradient.
  • any combination of multichannel data that can support an effective spectral estimate from aliased input, and that may return in output even a processed version of the input. For example, by applying GMP using measured Y and Z gradients in addition to the pressure as the input, a deghosted version of the input pressure maybe returned in output, resulting in the noise attenuation effect described above.
  • the spectral estimate and/or reconstruction of seismic wavefield are possible also if gradient measurements are not available.
  • the noise attenuation may be achieved by the separation of noisy components from signal components in the estimation, e.g., rather by the technique of estimating the spectrum.
  • spectral estimate and/or reconstruction may be achieved without processing gradient measurements.
  • a signal of interest and a noise signal may be identified in the reconstructed wavefield, based at least partially on different spectral characteristics of the signal of interest and the noise signal.
  • the signal of interest may be removed from the reconstructed wavefield.
  • the multi-channel interpolation and/or spatial analysis techniques can provide an identification of the one or more of the wavefield components associated with a signal of interest in the seismic data, which can be used to identify and remove the signal of interest. Once the signal of interest is identified and removed, at the measurement positions, what remains is a representation of the measured near surface noise which forms the noise model.
  • adaptive filters may be utilized to subtract the noise model from the seismic data.
  • the spatial gradient may be used to allow multi-channel interpolation and/or spatial analysis techniques to distinguish between what is signal and what is noise on the input measurements.
  • a noise model may be generated from the energy that has been classified as noise in this process. This may not call for an adaptive filter to be used to match the noise model to the noise, and subtract it, as the noise model is already an accurate estimate of the noise.
  • one or more filters may be applied during the multi-channel interpolation and/or spatial analysis techniques to identify and remove the signal of interest from the seismic data.
  • near surface noise may be more challenging to interpolate: its reconstruction in between input position may still be subject to residual aliasing. In this example, it may easier to isolate the signal of interest than the near surface noise when filtering the output.
  • the multi-channel interpolation and/or spatial analysis techniques may be based on filtering the reflected signal out of the reconstructed wavefield.
  • a reconstruction technique to interpolate the data may be performed from a coarse grid to a fine grid.
  • a reconstruction technique to interpolate the data may be performed from a coarse grid to a course grid.
  • the multi-channel interpolation and/or spatial analysis techniques may selecting as output the basis functions that are associated to the noise model to be matched, ignoring the components that are likely describing the signals of interest. This example may produce an output of the model of noise, which can be returned at the input positions, and not call for a fine grid of interpolated traces and a directional filter post- reconstruction.
  • a noise attenuated signal for the seismic data may be determined based at least partially on the signal of interest from the reconstructed wavefield.
  • the noise model may be directly subtracted from the collected seismic data, at each measurement position, thereby removing the noise and leaving the noise attenuated signal representing the subsurface volume.
  • the subtraction may be optional, for example, the signal of interest in the reconstructed wavefield may be accurately determined by the multi-channel interpolation and/or spatial analysis techniques and used as the output.
  • the signal of interest from the reconstructed wavefield may be selected as the noise attenuated signal representing the subsurface volume.
  • Figures 5C and 5D illustrates a plot showing the seismic data shown in Figures 5A and 5B, after the noise model was subtracted from the measurement.
  • the vertical axis is time, measured in seconds, and the horizontal is space, in meters.
  • the vertical axis represents frequency, in Hz, and the horizontal represents the wavenumber, in 1/m units.
  • the method 400 can end, repeat, or return to any process.
  • Figure 6 illustrates another example of a method 600 of removing noise from seismic data, according to various embodiments.
  • the illustrated processes of the method 600 are examples and that any of the illustrated processes can be removed, additional processes may be added, and the order of the illustrated processes may be changed.
  • any of the illustrated processes of the method 600 may be performed by one or more computer systems, for example, the computer systems illustrated in Figures 1 A-1D, Figure 2, and Figure 3, described above, and/or the computer systems described below in Figures 15A and 15B.
  • seismic data including measurements of a wavefield sampled on a course grid and an estimate of a spatial gradient of the wavefield by combining a group of measurements may be received.
  • the seismic data may be received by one or more computer systems, for examples, any of the computer systems described herein.
  • the seismic data may include measurements measured directly by seismic survey tools and measurements estimated.
  • the measurements measured directly by seismic survey tools may include direct sensor measurements of acoustic return signal, and direct measurement of sensor locations.
  • measurements that are estimated may include: input wavefield estimated as the average of two measurements form a couple of two receivers (i.e. pair of receiver); and input gradient estimated by finite difference operator applied to measurements of a couple of two receivers.
  • Figures 5 A and 5B illustrated an example of the seismic data.
  • the seismic data may be generated and collected by using the measurements of couples of receivers in a land seismic acquisition.
  • the couples of receivers may be positioned in a layout of a two-dimensional (2D) line.
  • the couples of receivers may be positioned at a short distance from each other, e.g., approximately 5m or less.
  • Each couple of receivers may be positioned at large distance from the closest ones, e.g., approximately 50m.
  • the average of each couple of traces collected may associated to the mid-position between the couples of receivers, and considered as an estimate of the wavefield there.
  • the normalized difference of each couple of traces was used as an estimate of the gradient of the wavefield at the midpoint between the couples of two receivers.
  • MFMAP/GMP/ FD-MTMAP may be performed to interpolate the seismic data.
  • MTMAP may be used to interpolate the seismic data.
  • GMP may be utilized to interpolate the data.
  • FD-MTMAP may be utilized to interpolate the data.
  • both estimated wavefield and gradients may be input to MFMAP to reconstruct the wavefield on a densely spaced grid of virtual receivers.
  • the nominal position of the input traces may be part of the output grid.
  • the nominal position of the input traces may not be used, as long as the information related to these positions is recoverable by processing the data.
  • the seismic wavefield of interest is characterized by seismic reflections perceived with a relatively broad temporal bandwidth and a relatively fast apparent velocity.
  • the available measurements are affected by strong near surface effects, particularly energetic at low frequency, with a slow apparent velocity, hence resulting severely aliased when measured at 50m interval.
  • MFMAP may reconstruct the aliased near- surface energy.
  • both, seismic signal of interest and near surface noise are finely sampled and can be processed with linear digital filters operating in time and space.
  • linear digital filters operating in time and space.
  • the property of MTMAP may estimate the reconstructed wavefield by matching iteratively the input measurements until these are fully explained by a set of basis functions. These basis functions are used to reconstruct the wavefield in between the input positions. The match with the measured traces at input positions may provide an accurate reconstruction.
  • the gradient may be used to allow MTMAP to distinguish between what is signal and what is noise on the input measurements, and to generate a noise model from the energy that has been classified as noise in this process. This may not call for an adaptive filter to be used to match the noise model to the noise, and subtract it, as the noise model is already an accurate estimate of the noise.
  • GMP or FD-MFMAP may be utilized.
  • the relation between the data may be taken into account and finite difference approximations of the signal and gradients may be obtained by processing the measurements.
  • the interpolated data may be filtered to separate signal of interest and noise with different spectral characteristics.
  • a noise model may be formed based on the filtered noise data.
  • the filtered noise data may be decimated at the input locations.
  • the noise model may be formed by removing the filtered noise at any positions except the input sample positions, e.g., location of the seismic sensors that collected the seismic data.
  • the noise model may be subtracted from the measurements.
  • the noise model may be directly subtracted from the collected seismic data, at each measurement position, thereby removing the noise and leaving the signal of interest representing the subsurface volume. After 610, the method 600 can end, repeat, or return to any process.
  • Figures 5C and 5D illustrates a plot showing the seismic data shown in Figures 5A and 5B, after the noise model was subtracted from the measurement.
  • the vertical axis is time, measured in seconds, and the horizontal is space, in meters.
  • the vertical axis represents frequency, in Hz, and the horizontal represents the wavenumber, in 1/m units.
  • Figure 7 illustrates another example of a method 700 of removing noise from seismic data, according to various embodiments.
  • the illustrated processes of the method 700 are examples and that any of the illustrated processes can be removed, additional processes may be added, and the order of the illustrated processes may be changed.
  • any of the illustrated processes of the method 700 may be performed by one or more computer systems, for example, the computer systems illustrated in Figures 1 A-1D, Figure 2, and Figure 3, described above, and/or the computer systems described below in Figures 15A and 15B.
  • seismic data including measurements of a wavefield and a spatial gradient of the wavefield sampled on a course grid may be received.
  • the seismic data may be received by one or more computer systems, for example, any of the computer systems described herein.
  • the seismic data may include measurements measured directly by seismic survey tools and measurements estimated.
  • the measurements measured directly by seismic survey tools may include direct sensor measurements of acoustic return signal, the direct sensor measurements of the spatial gradient, and direct measurement of sensor locations.
  • measurements that are estimated may include: input wavefield estimated as the average of two measurements form a couple of two receivers (i.e. pair of receivers); and input gradient estimated by finite difference operator applied to measurements of a couple of two receivers.
  • the seismic data may be interpolated by processing wavefield and spatial gradient.
  • MTMAP may be used to interpolate the seismic data.
  • GMP may be utilized to interpolate the data.
  • FD-MTMAP may be utilized to interpolate the data.
  • MEVIAP, GMP, and FD-MEVIAP may utilize both the directly measured seismic signal and gradient, and the estimated gradient. The property of MEVIAP may estimate the reconstructed wavefield by matching iteratively the input measurements until these are fully explained by a set of basis functions. These basis functions are used to reconstruct the wavefield in between the input positions. The match with the measured traces at input positions may provide an accurate reconstruction.
  • the gradient may be used to allow MEVIAP to distinguish between what is signal and what is noise on the input measurements, and to generate a noise model from the energy that has been classified as noise in this process. This may not call for an adaptive filter to be used to match the noise model to the noise, and subtract it, as the noise model is already an accurate estimate of the noise.
  • GMP or FD-MEVIAP may be utilized.
  • the relation between the data may be taken into account and finite difference approximations of the signal and gradients may be obtained by processing the measurements.
  • the interpolated data may be filtered to separate signal of interest and noise with different spectral characteristics.
  • a noise model may be formed based on the filtered noise.
  • the filtered noise may be decimated, at the input sample positions, to form the filtered noise.
  • the noise model may be formed by removing the filtered noise at any positions expect the input sample positions, e.g., location of the seismic sensors that collected the seismic data.
  • the noise model may be subtracted from the measurements.
  • the noise model may be directly subtracted from the collected seismic data, at each measurement position, thereby removing the noise and leaving the signal of interest representing the subsurface volume.
  • the method 700 can end, repeat, or return to any process.
  • Figure 8 illustrates another example of a method 800 of removing noise from seismic data, according to various embodiments.
  • the illustrated processes of the method 800 are examples and that any of the illustrated processes can be removed, additional processes may be added, and the order of the illustrated processes may be changed.
  • any of the illustrated processes of the method 800 may be performed by one or more computer systems, for example, the computer systems illustrated in Figures 1 A-1D, Figure 2, and Figure 3, described above, and/or the computer systems described below in Figures 15A and 15B.
  • seismic data including measurements of a wavefield and a spatial gradient of the wavefield sampled on a course grid may be received.
  • the seismic data may be received by one or more computer systems, for examples, any of the computer systems described herein.
  • the seismic data may include measurements measured directly by seismic survey tools and measurements estimated.
  • the measurements measured directly by seismic survey tools may include direct sensor measurements of acoustic return signal, the direct sensor measurements of the spatial gradient, and direct measurement of sensor locations.
  • MEVIAP/ GMP/FD-MFM AP may be performed to interpolate the seismic data.
  • MTMAP may be used to interpolate the seismic data.
  • GMP may be utilized to interpolate the data.
  • FD-MTMAP may be utilized to interpolate the data.
  • MTMAP, GMP, and FD-MFMAP may utilize both the directly measured seismic signal and gradient.
  • the property of MIMAP may estimate the reconstructed wavefield by matching iteratively the input measurements until these are fully explained by a set of basis functions. These basis functions are used to reconstruct the wavefield in between the input positions. The match with the measured traces at input positions may provide an accurate reconstruction.
  • the gradient may be used to allow MFMAP to distinguish between what is signal and what is noise on the input measurements, and to generate a noise model from the energy that has been classified as noise in this process. This may not call for an adaptive filter to be used to match the noise model to the noise, and subtract it, as the noise model is already an accurate estimate of the noise.
  • the interpolated data may be filtered to separate signal of interest and noise with different spectral characteristics.
  • the filtered noise may be decimated, at the input sample positions, to form a noise model.
  • the noise model may be formed by removing the filtered noise at any positions expect the input sample positions, e.g., location of the seismic sensors that collected the seismic data.
  • the noise model may be subtracted from the measurements.
  • the noise model may be directly subtracted from the collected seismic data, at each measurement position, thereby removing the noise and leaving the signal of interest representing the subsurface volume.
  • the method 800 can end, repeat, or return to any process.
  • Figure 9 illustrates another example of a method 900 of removing noise from seismic data, according to various embodiments.
  • the illustrated processes of the method 900 are examples and that any of the illustrated processes can be removed, additional processes may be added, and the order of the illustrated processes may be changed.
  • any of the illustrated processes of the method 900 may be performed by one or more computer systems, for example, the computer systems illustrated in Figures 1 A-1D, Figure 2, and Figure 3, described above, and/or the computer systems described below in Figures 15A and 15B.
  • seismic data including measurements of a wavefield sampled on a course grid and an estimate of a spatial gradient of the wavefield obtained by combining a group of measurement may be received.
  • the seismic data may be received by one or more computer systems, for examples, any of the computer systems described herein.
  • the seismic data may include measurements measured directly by seismic survey tools and measurements estimated.
  • the measurements measured directly by seismic survey tools may include direct sensor measurements of acoustic return signal, and direct measurement of sensor locations.
  • measurements that are estimated may include: input wavefield estimated as the average of two measurements form a couple of two receivers (i.e. pair of receivers); and input gradient estimated by finite difference operator applied to measurements of a couple of two receivers.
  • Figures 5A and 5B illustrated an example of the seismic data.
  • the seismic data may be generated and collected by using the measurements of couples of receivers in a land seismic acquisition.
  • the couples of receivers may be positioned in a layout of a two-dimensional (2D) line.
  • the couples of receivers may be positioned at a short distance from each other, e.g., approximately 5m or less.
  • Each couple of receivers may be positioned at large distance from the closest ones, e.g., approximately 50m.
  • the average of each couple of traces collected may associated to the mid-position between the couples of receivers, and considered as an estimate of the wavefield there.
  • the normalized difference of each couple of traces was used as an estimate of the gradient of the wavefield at the midpoint between the couples of two receivers.
  • E-GMP may be performed to interpolate the seismic data.
  • E-GMP may estimate a reconstructed wavefield by matching iteratively the input measurements until these are fully explained by a set of basis functions. For example, the E-GMP may match iteratively the input measurements by trying to estimate their spatial spectrum and determine iteratively which basis functions are likely components of the wavefield.
  • E-GMP may achieve reconstruction by estimating the spectrum of the data iteration after iteration.
  • E-GMP may process a more generic set of data, provided that the relation between the input measurements and the desired output is known. For instance, E-GMP may determine the components of the wavefield where spatial gradient is not directly measured, but an estimated spatial gradient is available or determined.
  • an estimated spatial gradient may be generated by a finite difference operator on the measurements of a group of receivers.
  • E-GMP may consider the actual model used to generate the estimated gradient, and may matches the data respecting such model. E-GMP may use this estimated gradient for interpolation purposes without suffering of the approximation error in estimating the gradient.
  • E-GMP may use the finite different operator to model the relation between a signal of interest and its estimated spatial gradient. In E-GMP, the use of an actual operator used to estimate the spatial gradient may provide an improved match than the use of the ideal gradient operator that would be used within GMP.
  • the interpolated data may be filtered to separate signal of interest and noise with different spectral characteristics.
  • both, seismic signal of interest and near surface noise are finely sampled and can hence be processed with linear digital filters operating in time and space.
  • linear digital filters operating in time and space.
  • a noise model may be formed based on the filtered noise.
  • the filtered noise may be decimated, at the input sample positions, to form the noise model.
  • the noise model may be formed by removing the filtered noise at any positions expect the input sample positions, e.g., location of the seismic sensors that collected the seismic data.
  • the noise model may be subtracted from the measurements.
  • the noise model may be directly subtracted from the collected seismic data, at each measurement position, thereby removing the noise and leaving the signal of interest representing the subsurface volume.
  • the method 900 can end, repeat, or return to any process.
  • Figure 10 illustrates another example of a method 1000 of removing noise from seismic data, according to various embodiments.
  • the illustrated processes of the method 1000 are examples and that any of the illustrated processes can be removed, additional processes may be added, and the order of the illustrated processes may be changed.
  • any of the illustrated processes of the method 1000 may be performed by one or more computer systems, for example, the computer systems illustrated in Figures 1 A- ID, Figure 2, and Figure 3, described above, and/or the computer systems described below in Figures 15A and 15B.
  • seismic data including measurements of a wavefield and a spatial gradient of the wavefield sampled on a course grid may be received.
  • the seismic data may be received by one or more computer systems, for examples, any of the computer systems described herein.
  • the seismic data may include measurements measured directly by seismic survey tools and measurements estimated.
  • the measurements measured directly by seismic survey tools may include direct sensor measurements of acoustic return signal, the direct sensor measurements of the spatial gradient, and direct measurement of sensor locations.
  • measurements that are estimated may include: input wavefield estimated as the average of two measurements form a couple of two receivers (i.e. pair of receiver); and input gradient estimated by finite difference operator applied to measurements of a couple of two receivers.
  • a spatial spectrum of the seismic data may be estimated based at least partially on the wavefield and spatial gradient.
  • spectral components likely containing signal of interest may be separated from spectral components likely containing noise.
  • the spectral components of the noise at measurement positions may be evaluated to form a noise model.
  • MTMAP may be used to interpolate the seismic data.
  • GMP may be utilized to interpolate the data.
  • FD-MTMAP may be utilized to interpolate the data.
  • E-GMP may be utilized to interpolate the data.
  • MEVI AP/ GMP/E-GMP/FD-MEVI AP may be performed by selecting as output the basis functions that are associated to the noise model to be matched, ignoring the components that are likely describing reflected signals.
  • MEVIAP/ GMP/E-GMP/FD-MEVI AP may produce as output the noise model. This can be returned at the input positions, thus avoiding a fine grid of interpolated traces and a directional filter post- reconstruction.
  • the noise model may be subtracted from the measurements.
  • the noise model may be directly subtracted from the collected seismic data, at each measurement position, thereby removing the noise and leaving the signal of interest representing the subsurface volume.
  • the method 1000 can end, repeat, or return to any process.
  • Figure 1 1 illustrates another example of a method 1100 of removing noise from seismic data, according to various embodiments.
  • the illustrated processes of the method 1 100 are examples and that any of the illustrated processes can be removed, additional processes may be added, and the order of the illustrated processes may be changed.
  • any of the illustrated processes of the method 1 100 may be performed by one or more computer systems, for example, the computer systems illustrated in Figures 1 A- ID, Figure 2, and Figure 3, described above, and/or the computer systems described below in Figures 15A and 15B.
  • seismic data including measurements of a wavefield and a spatial gradient of the wavefield sampled on a course grid may be received.
  • the seismic data may be received by one or more computer systems, for examples, any of the computer systems described herein.
  • the seismic data may include measurements measured directly by seismic survey tools and measurements estimated.
  • the measurements measured directly by seismic survey tools may include direct sensor measurements of acoustic return signal, the direct sensor measurements of the spatial gradient, and direct measurement of sensor locations.
  • measurements that are estimated may include: input wavefield estimated as the average of two measurements form a couple of two receivers (i.e. pair of receiver); and input gradient estimated by finite difference operator applied to measurements of a couple of two receivers.
  • a spatial spectrum of the seismic data may be estimated using MFMAP/GMP/FD- MFMAP.
  • spectral components likely containing signal of interest may be separated within MFMAP/GMP/FD- MFMAP from spectral components likely containing noise.
  • a noise model at the measurement position may be created while running MFMAP/ GMP/FD-MFM AP .
  • MFMAP may be used to interpolate the seismic data.
  • GMP may be utilized to interpolate the data.
  • FD-MFMAP may be utilized to interpolate the data.
  • MIMAP/GMP/FD-MIMAP may be performed by selecting as output the basis functions that are associated to the noise model to be matched, ignoring the components that are likely describing reflected signals.
  • MIMAP/GMP/FD- MIMAP may produce as output the noise model: this can be returned at the input positions, hence not requiring a fine grid of interpolated traces and a directional filter post-reconstruction.
  • the noise model may be subtracted from the measurements.
  • the noise model may be directly subtracted from the collected seismic data, at each measurement position, thereby removing the noise and leaving the signal of interest representing the subsurface volume.
  • the method 1100 can end, repeat, or return to any process.
  • Figure 12 illustrates another example of a method 1200 of removing noise from seismic data, according to various embodiments.
  • the illustrated processes of the method 1200 are examples and that any of the illustrated processes can be removed, additional processes may be added, and the order of the illustrated processes may be changed.
  • any of the illustrated processes of the method 1200 may be performed by one or more computer systems, for example, the computer systems illustrated in Figures 1 A- ID, Figure 2, and Figure 3, described above, and/or the computer systems described below in Figures 15A and 15B.
  • seismic data including measurements of a wavefield sampled on a course grid and an estimate of a spatial gradient of the wavefield obtained by combining a group of measurements may be received.
  • the seismic data may be received by one or more computer systems, for examples, any of the computer systems described herein.
  • the seismic data may include measurements measured directly by seismic survey tools and measurements estimated.
  • the measurements measured directly by seismic survey tools may include direct sensor measurements of acoustic return signal, the direct sensor measurements of the spatial gradient, and direct measurement of sensor locations.
  • measurements that are estimated may include: input wavefield estimated as the average of two measurements form a couple of two receivers (i.e. pair of receiver); and input gradient estimated by finite difference operator applied to measurements of a couple of two receivers.
  • spatial spectrum of the seismic data may be estimated using E-GMP.
  • spectral components likely containing signal of interest may be separated within E-GMP from spectral components likely containing noise.
  • a noise model at the measurement position may be created while running E-GMP.
  • E-GMP may be performed by selecting as output the basis functions that are associated to the noise model to be matched, ignoring the components that are likely describing reflected signals.
  • E-GMP may produce as output the noise model: this can be returned at the input positions, hence not requiring a fine grid of interpolated traces and a directional filter post-reconstruction.
  • the noise model may be subtracted from the measurements.
  • the noise model may be directly subtracted from the collected seismic data, at each measurement position, thereby removing the noise and leaving the signal of interest representing the subsurface volume.
  • the method 1200 can end, repeat, or return to any process.
  • Figure 13 illustrates another example of a method 1300 of removing noise from seismic data, according to various embodiments.
  • the illustrated processes of the method 1300 are examples and that any of the illustrated processes can be removed, additional processes may be added, and the order of the illustrated processes may be changed.
  • any of the illustrated processes of the method 1300 may be performed by one or more computer systems, for example, the computer systems illustrated in Figures 1 A- ID, Figure 2, and Figure 3, described above, and/or the computer systems described below in Figures 15A and 15B.
  • seismic data including measurements of a wavefield and a spatial gradient of the wavefield sampled on a course grid may be received.
  • the seismic data may be received by one or more computer systems, for examples, any of the computer systems described herein.
  • the seismic data may include measurements measured directly by seismic survey tools and measurements estimated.
  • the measurements measured directly by seismic survey tools may include direct sensor measurements of acoustic return signal, the direct sensor measurements of the spatial gradient, and direct measurement of sensor locations.
  • measurements that are estimated may include: input wavefield estimated as the average of two measurements form a couple of two receivers (i.e. pair of receiver); and input gradient estimated by finite difference operator applied to measurements of a couple of two receivers.
  • the seismic data may be interpolated by processing wavefield and spatial gradient.
  • the interpolated data may be filtered to separate signal of interest and noise with different spectral characteristics.
  • the one or more of the wavefield components associated with a signal of interest in the seismic data may be accurately determined by the multi-channel interpolation and/or spatial analysis techniques and used as the output.
  • the multichannel interpolation and/or spatial analysis techniques may reconstruct the wavefield for the signal of interest and the noise with different spectral characteristics. Then, the multi-channel interpolation and/or spatial analysis techniques may filter out the noise with different spectral characteristics.
  • MTMAP may be used to interpolate the seismic data.
  • GMP may be utilized to interpolate the data.
  • FD-MTMAP may be utilized to interpolate the data.
  • E-GMP may be utilized to interpolate the data.
  • an estimated signal of interest may be decimated at the input sample positions.
  • the estimated signal of interest may be formed by removing the signals of interest at any positions expect the input sample positions, e.g., location of the seismic sensors that collected the seismic data.
  • the method 1300 can end, repeat, or return to any process.
  • Figure 14 illustrates another example of a method 1400 of removing noise from seismic data, according to various embodiments.
  • the illustrated processes of the method 1400 are examples and that any of the illustrated processes can be removed, additional processes may be added, and the order of the illustrated processes may be changed.
  • any of the illustrated processes of the method 1400 may be performed by one or more computer systems, for example, the computer systems illustrated in Figures 1 A- ID, Figure 2, and Figure 3, described above, and/or the computer systems described below in Figures 15A and 15B.
  • seismic data including measurements of a wavefield and a spatial gradient of the wavefield sampled on a course grid may be received.
  • the seismic data may be received by one or more computer systems, for examples, any of the computer systems described herein.
  • the seismic data may include measurements measured directly by seismic survey tools and measurements estimated.
  • the measurements measured directly by seismic survey tools may include direct sensor measurements of acoustic return signal, the direct sensor measurements of the spatial gradient, and direct measurement of sensor locations.
  • measurements that are estimated may include: input wavefield estimated as the average of two measurements form a couple of two receivers (i.e. pair of receiver); and input gradient estimated by finite difference operator applied to measurements of a couple of two receivers.
  • a spatial spectrum of the seismic data may be estimated based at least partially on the wavefield and spatial gradient.
  • spectral components likely containing signal of interest may be separated from spectral components likely containing noise.
  • the spectral components of the signal of interest may be evaluated at measurement positions.
  • the one or more of the wavefield components associated with a signal of interest in the seismic data may be accurately determined by the multi-channel interpolation and/or spatial analysis techniques and used as the output.
  • the multichannel interpolation and/or spatial analysis techniques may reconstruct spectral components likely containing signal of interest may be separated from spectral components likely containing noise. Then, the multi-channel interpolation and/or spatial analysis techniques may remove spectral components likely containing noise.
  • MTMAP may be used to interpolate the seismic data.
  • GMP may be utilized to interpolate the data.
  • FD-MTMAP may be utilized to interpolate the data.
  • E-GMP may be utilized to interpolate the data.
  • FIGs 15A and 15B illustrate another flowchart of an example of a method 1500 for noise attenuation, according to an embodiment.
  • the method 1500 may include obtaining seismic data including measurements of a seismic wavefield for at least one subsurface volume, as at 1502.
  • obtaining the seismic data may include rreceiving acoustic signals collected for the at least one subsurface volume, as at 1504.
  • the method 1500 may include obtaining at least one component of a spatial gradient of the measurements of the seismic wavefield, as at 1506.
  • obtaining the at least one component of the spatial gradient may include determining the at least one component of the spatial gradient from the acoustic signals, as at 1508.
  • obtaining the at least one component of the spatial gradient may include receiving the at least one component of the spatial gradient collected for the at least one subsurface volume, as at 1510.
  • the method 1500 may include determining a representation of the seismic data based at partially on the measurements of the seismic wavefield and the at least one component of the spatial gradient, as at 1512.
  • determining the representation of the seismic data may include processing the measurements and the at least one component of the spatial gradient using an analysis process, as at 1514.
  • the analysis process may include at least one of multichannel interpolation by matching pursuit, generalized matching pursuit, an extended generalized matching pursuit, a finite difference multichannel interpolation by matching pursuit, or a greedy algorithm, as at 1516.
  • the method 1500 may include identifying a signal of interest and noise in the representation of the seismic data based at least partially on different characteristics of the signal of interest and the noise in a domain of the representation, as at 1518.
  • the method 1500 may include calculating at least one of a signal model or a noise model from the signal of interest and noise identified in the representation of the seismic data, as at 1520.
  • calculating at least one of the signal model or the noise model may include removing the signal of interest from the representation of the seismic data to form the noise model, as at 1522.
  • the method 1500 may include determining a noise attenuated signal for the seismic data based at least partially on the at least one of the signal model or the noise model, as at 1524.
  • determining the noise attenuated signal for the seismic data may include subtracting the noise model from the measurements of the seismic wavefield, as at 1526.
  • determining the noise attenuated signal for the seismic data may include selecting the signal of interest from the representation of the seismic data as the noise attenuated signal, as at 1528.
  • the methods of the present disclosure may be executed by a computing system.
  • Figure 16 illustrates an example of such a computing system 1600, in accordance with some embodiments.
  • the computing system 1600 may include a computer or computer system 1601 A, which may be an individual computer system 1601 A or an arrangement of distributed computer systems.
  • the computer system 1601 A includes one or more analysis modules 1602 that are configured to perform various tasks according to some embodiments, such as one or more methods disclosed herein. To perform these various tasks, the analysis module 1602 executes independently, or in coordination with, one or more processors 1604, which is (or are) connected to one or more storage media 1606.
  • the processor(s) 1604 is (or are) also connected to a network interface 1607 to allow the computer system 1601 A to communicate over a data network 1609 with one or more additional computer systems and/or computing systems, such as 160 IB, 1601C, and/or 160 ID (note that computer systems 160 IB, 1601C and/or 160 ID may or may not share the same architecture as computer system 1601A, and may be located in different physical locations, e.g., computer systems 1601 A and 160 IB may be located in a processing facility, while in communication with one or more computer systems such as 1601C and/or 160 ID that are located in one or more data centers, and/or located in varying countries on different continents).
  • additional computer systems and/or computing systems such as 160 IB, 1601C, and/or 160 ID
  • computer systems 1601 A and 160 IB may be located in a processing facility, while in communication with one or more computer systems such as 1601C and/or 160 ID that are located in one or more data centers, and/or located in varying
  • a processor can include a microprocessor, microcontroller, processor module or subsystem, programmable integrated circuit, programmable gate array, or another control or computing device.
  • the storage media 1606 can be implemented as one or more computer-readable or machine-readable storage media. Note that while in some example embodiments of Figure 16 storage media 1606 is depicted as within computer system 1601A, in some embodiments, storage media 1606 may be distributed within and/or across multiple internal and/or external enclosures of computing system 1601A and/or additional computing systems.
  • Storage media 1606 may include one or more different forms of memory including semiconductor memory devices such as dynamic or static random access memories (DRAMs or SRAMs), erasable and programmable read-only memories (EPROMs), electrically erasable and programmable read-only memories (EEPROMs) and flash memories, magnetic disks such as fixed, floppy and removable disks, other magnetic media including tape, optical media such as compact disks (CDs) or digital video disks (DVDs), BLUERAY ® disks, or other types of optical storage, or other types of storage devices.
  • semiconductor memory devices such as dynamic or static random access memories (DRAMs or SRAMs), erasable and programmable read-only memories (EPROMs), electrically erasable and programmable read-only memories (EEPROMs) and flash memories
  • magnetic disks such as fixed, floppy and removable disks, other magnetic media including tape
  • optical media such as compact disks (CDs) or digital video disks (DVDs), BLUERAY ® disks,
  • Such computer-readable or machine-readable storage medium or media is (are) considered to be part of an article (or article of manufacture).
  • An article or article of manufacture can refer to any manufactured single component or multiple components.
  • the storage medium or media can be located either in the machine running the machine-readable instructions, or located at a remote site from which machine-readable instructions can be downloaded over a network for execution.
  • computing system 1600 contains one or more noise mitigation module(s) 1608.
  • computer system 1601 A includes the noise mitigation module 1608.
  • a single noise mitigation module may be used to perform at least some aspects of one or more embodiments of the methods disclosed herein.
  • a plurality of noise mitigation modules may be used to perform at least some aspects of methods disclosed herein.
  • computing system 1600 is but one example of a computing system, and that computing system 1600 may have more or fewer components than shown, may combine additional components not depicted in the example embodiment of Figure 16, and/or computing system 1600 may have a different configuration or arrangement of the components depicted in Figure 16.
  • the various components shown in Figure 16 may be implemented in hardware, software, or a combination of both hardware and software, including one or more signal processing and/or application specific integrated circuits.
  • the steps in the processing methods described herein may be implemented by running one or more functional modules in information processing apparatus such as general purpose processors or application specific chips, such as ASICs, FPGAs, PLDs, or other appropriate devices.
  • ASICs general purpose processors or application specific chips, such as ASICs, FPGAs, PLDs, or other appropriate devices.
  • Geologic interpretations, models and/or other interpretation aids may be refined in an iterative fashion; this concept is applicable to methods as discussed herein.
  • This can include use of feedback loops executed on an algorithmic basis, such as at a computing device (e.g., computing system 1600, Figure 16), and/or through manual control by a user who may make determinations regarding whether a given step, action, template, model, or set of curves has become sufficiently accurate for the evaluation of the subsurface three-dimensional geologic formation under consideration.
  • embodiments may be utilized in conjunction with a handheld system (i.e., a phone, wrist or forearm mounted computer, tablet, or other handheld device), portable system (i.e., a laptop or portable computing system), a fixed computing system (i.e., a desktop, server, cluster, or high performance computing system), or across a network (i.e., a cloud- based system).
  • a handheld system i.e., a phone, wrist or forearm mounted computer, tablet, or other handheld device
  • portable system i.e., a laptop or portable computing system
  • a fixed computing system i.e., a desktop, server, cluster, or high performance computing system
  • a network i.e., a cloud- based system

Landscapes

  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Acoustics & Sound (AREA)
  • Environmental & Geological Engineering (AREA)
  • Geology (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

La présente invention concerne un procédé consistant à obtenir des données sismiques comprenant des mesures d'un champ d'ondes sismiques et à obtenir au moins une composante d'un gradient spatial des mesures du champ d'ondes sismiques. Le procédé comprend la détermination d'une représentation des données sismiques au moins partiellement sur la base des mesures et de la ou des composantes du gradient spatial et l'identification d'un signal d'intérêt et d'un bruit sur la représentation des données sismiques au moins partiellement sur la base des différentes caractéristiques du signal d'intérêt et du bruit dans un domaine de la représentation. Le procédé comprend le calcul d'un modèle de signal et/ou d'un modèle de bruit à partir du signal d'intérêt et du bruit identifiés sur la représentation des données sismiques et la détermination d'un signal à bruit atténué pour les données sismiques au moins partiellement sur la base du modèle de signal et/ou du modèle de bruit.
PCT/US2016/027261 2015-04-14 2016-04-13 Génération d'un modèle précis de bruit et sa soustraction de données sismiques WO2016168280A1 (fr)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
US201562147087P 2015-04-14 2015-04-14
US62/147,087 2015-04-14

Publications (1)

Publication Number Publication Date
WO2016168280A1 true WO2016168280A1 (fr) 2016-10-20

Family

ID=57126240

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/US2016/027261 WO2016168280A1 (fr) 2015-04-14 2016-04-13 Génération d'un modèle précis de bruit et sa soustraction de données sismiques

Country Status (1)

Country Link
WO (1) WO2016168280A1 (fr)

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20190293813A1 (en) * 2017-11-20 2019-09-26 Conocophillips Company Offshore application of non-uniform optimal sampling survey design
US10539695B2 (en) 2014-04-28 2020-01-21 Westerngeco L.L.C. Wavefield reconstruction
US10775522B2 (en) 2016-06-15 2020-09-15 Schlumberger Technology Corporation Systems and methods for attenuating noise in seismic data and reconstructing wavefields based on the seismic data
US10928535B2 (en) 2015-05-01 2021-02-23 Reflection Marine Norge As Marine vibrator directive source survey
US10948615B2 (en) 2015-12-02 2021-03-16 Westerngeco L.L.C. Land seismic sensor spread with adjacent multicomponent seismic sensor pairs on average at least twenty meters apart
US10996359B2 (en) 2015-05-05 2021-05-04 Schlumberger Technology Corporation Removal of acquisition effects from marine seismic data
US20220050222A1 (en) * 2019-01-15 2022-02-17 Schlumberger Technology Corporation Residual signal detection for noise attenuation

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20050075791A1 (en) * 2002-08-05 2005-04-07 Patrick Lailly Method of forming a model representative of the distribution of a physical quantity in an underground zone, free of the effect of correlated noises contained in exploration data
US20130182533A1 (en) * 2012-01-12 2013-07-18 Westerngeco L.L.C. Attentuating noise acquired in an energy measurement
US20140200820A1 (en) * 2013-01-11 2014-07-17 Westerngeco L.L.C. Wavefield extrapolation and imaging using single- or multi-component seismic measurements
US20140200854A1 (en) * 2013-01-15 2014-07-17 Westerngeco L.L.C. Methods and computing systems for survey design and modeling workflow for towed multimeasurement seismic streamer data
US20140286127A1 (en) * 2013-03-21 2014-09-25 Westerngeco L.L.C. Processing seismic data to attenuate vibration noise

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20050075791A1 (en) * 2002-08-05 2005-04-07 Patrick Lailly Method of forming a model representative of the distribution of a physical quantity in an underground zone, free of the effect of correlated noises contained in exploration data
US20130182533A1 (en) * 2012-01-12 2013-07-18 Westerngeco L.L.C. Attentuating noise acquired in an energy measurement
US20140200820A1 (en) * 2013-01-11 2014-07-17 Westerngeco L.L.C. Wavefield extrapolation and imaging using single- or multi-component seismic measurements
US20140200854A1 (en) * 2013-01-15 2014-07-17 Westerngeco L.L.C. Methods and computing systems for survey design and modeling workflow for towed multimeasurement seismic streamer data
US20140286127A1 (en) * 2013-03-21 2014-09-25 Westerngeco L.L.C. Processing seismic data to attenuate vibration noise

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10539695B2 (en) 2014-04-28 2020-01-21 Westerngeco L.L.C. Wavefield reconstruction
US10928535B2 (en) 2015-05-01 2021-02-23 Reflection Marine Norge As Marine vibrator directive source survey
US10996359B2 (en) 2015-05-05 2021-05-04 Schlumberger Technology Corporation Removal of acquisition effects from marine seismic data
US10948615B2 (en) 2015-12-02 2021-03-16 Westerngeco L.L.C. Land seismic sensor spread with adjacent multicomponent seismic sensor pairs on average at least twenty meters apart
US10775522B2 (en) 2016-06-15 2020-09-15 Schlumberger Technology Corporation Systems and methods for attenuating noise in seismic data and reconstructing wavefields based on the seismic data
US20190293813A1 (en) * 2017-11-20 2019-09-26 Conocophillips Company Offshore application of non-uniform optimal sampling survey design
US12259511B2 (en) * 2017-11-20 2025-03-25 Shearwater Geoservices Software Inc. Offshore application of non-uniform optimal sampling survey design
US20220050222A1 (en) * 2019-01-15 2022-02-17 Schlumberger Technology Corporation Residual signal detection for noise attenuation
US11988790B2 (en) * 2019-01-15 2024-05-21 Schlumberger Technology Corporation Residual signal detection for noise attenuation

Similar Documents

Publication Publication Date Title
EP3209859B1 (fr) Fonction objective de temps de déplacement pour inversion de forme d'ondes complètes
WO2016168280A1 (fr) Génération d'un modèle précis de bruit et sa soustraction de données sismiques
EP3114312B1 (fr) Modélisation et filtrage de bruit cohérent dans des relevés sismiques
EP3058169B1 (fr) Sommation itérative de partitions d'images sismiques
WO2016161096A1 (fr) Traitement de données sismiques
EP3908859A1 (fr) Combinaison d'atténuation de bruit et de reconstruction de champ d'ondes dans un traitement sismique
US11531128B2 (en) Reconstruction of multi-shot, multi-channel seismic wavefields
EP3423870B1 (fr) Amélioration d'image au moyen d'images de partition sismique
US20240241279A1 (en) Data-driven separation of upgoing free-surface multiples for seismic imaging
US20170235002A1 (en) Joint Interpolation and Deghosting of Seismic Data
US11774620B2 (en) Interpolation of seismic data with time variant locations
EP3271549B1 (fr) Déparasitage de flûte sismique unique à espace de modèle étendu
WO2024145017A1 (fr) Séparation de source itérative à étages multiples avec acquisition antérieure à intervalles de temps
EP4363902A1 (fr) Remplacement direct d'arrivée pour imagerie sismique
WO2024049404A1 (fr) Déconvolution de champs d'ondes sismiques
WO2024049406A1 (fr) Nouveau référencement interférométrique, interpolation, et élimination de surface libre pour des données sismiques de fond océanique

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: 16780611

Country of ref document: EP

Kind code of ref document: A1

NENP Non-entry into the national phase

Ref country code: DE

122 Ep: pct application non-entry in european phase

Ref document number: 16780611

Country of ref document: EP

Kind code of ref document: A1