WO2008147657A2 - Method for identifying and removing multiples for imaging with beams - Google Patents

Method for identifying and removing multiples for imaging with beams Download PDF

Info

Publication number
WO2008147657A2
WO2008147657A2 PCT/US2008/062877 US2008062877W WO2008147657A2 WO 2008147657 A2 WO2008147657 A2 WO 2008147657A2 US 2008062877 W US2008062877 W US 2008062877W WO 2008147657 A2 WO2008147657 A2 WO 2008147657A2
Authority
WO
WIPO (PCT)
Prior art keywords
pegleg
multiples
selecting
beams
convolved
Prior art date
Application number
PCT/US2008/062877
Other languages
French (fr)
Other versions
WO2008147657A3 (en
Inventor
Tamas Nemeth
Yue Wang
Fuhao Qin
Ross N. Hill
Original Assignee
Chevron U.S.A. Inc.
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 Chevron U.S.A. Inc. filed Critical Chevron U.S.A. Inc.
Priority to BRPI0811841-8A2A priority Critical patent/BRPI0811841A2/en
Priority to CA2687095A priority patent/CA2687095C/en
Priority to CN200880021551.9A priority patent/CN101688922B/en
Priority to EA200971083A priority patent/EA017379B1/en
Priority to AU2008257075A priority patent/AU2008257075B2/en
Publication of WO2008147657A2 publication Critical patent/WO2008147657A2/en
Publication of WO2008147657A3 publication Critical patent/WO2008147657A3/en
Priority to NO20093563A priority patent/NO20093563L/en

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28Processing seismic data, e.g. analysis, for interpretation, for correction
    • G01V1/36Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
    • G01V1/364Seismic filtering
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00Details of seismic processing or analysis
    • G01V2210/50Corrections or adjustments related to wave propagation
    • G01V2210/56De-ghosting; Reverberation compensation

Definitions

  • This invention relates to seismic exploration and processing, and more specifically to imaging with beams and a method to predict multiples based on primaries by combining both model-driven and data-driven methodologies.
  • seismic prospecting techniques are commonly used to aid in the search for and the evaluation of subterranean hydrocarbon deposits.
  • one or more sources of seismic energy emit waves into a subsurface region of interest, such as a geologic formation. These waves enter the formation and may be scattered, e.g., by reflection or refraction, by subsurface seismic reflectors (i.e., interfaces between underground formations having different elastic properties).
  • the reflected signals are sampled or measured by one or more receivers, and the resultant data are recorded.
  • the recorded samples may be referred to as seismic data or a set of "seismic traces".
  • the seismic data may be analyzed to extract details of the structure and properties of the region of the earth being explored.
  • Seismic prospecting consists of three separate stages: data acquisition, data processing and data interpretation. The success of a seismic prospecting operation depends on satisfactory completion of all three stages.
  • the purpose of seismic exploration is to map or image a portion of the subsurface of the earth (a formation) by transmitting energy down into the ground and recording the "reflections” or “echoes” that return from the rock layers below.
  • the energy transmitted into the formation is typically sound and shear wave energy.
  • the downward-propagating sound energy may originate from various sources, such as explosions or seismic vibrators on land or air guns in marine environments.
  • Seismic exploration typically uses one or more energy sources and typically a large number of sensors or detectors.
  • the sensors that may be used to detect the returning seismic energy are usually geophones (land surveys) or hydrophones (marine surveys).
  • the energy source may be positioned at one or more locations near the surface of the earth above a geologic structure or formation of interest, referred to as shotpoints.
  • shotpoints Each time the source is activated, the source generates a seismic signal that travels downward through the earth and is at least partially reflected from discontinuities of various types in the subsurface, including reflections from "rock layer" boundaries.
  • a partial reflection of seismic signals may occur each time there is a change in the elastic properties of the subsurface materials. Reflected seismic signals are transmitted back to the surface of the earth, where they are recorded as a function of traveltime at a number of locations. The returning signals are digitized and recorded as a function of time (amplitude vs.. time).
  • a primary reflection is a sound wave that passes from the source to a receiver with a single reflection from a subsurface seismic reflector.
  • a multiple reflection is a wave that has reflected at least three times (up, down and back up again) before being detected by a receiver.
  • multiples are commonly characterized as short-path, implying that they interfere with their own primary reflection, or long-path, where they appear as separate events.
  • Model-driven methodologies generally use an earth model and the recorded data to predict or simulate multiples utilizing the estimated sea-bottom reflectivity function and calculated amplitude functions to model water-layer multiple reflections, those predicted multiples are then subtracted from the original data.
  • Other model-driven technologies utilize an earth model or reflectivity model to predict the stationary multiples.
  • the data-driven methodologies exploit the fact that primaries and multiples are physically related through a convolutional relationship and predict multiples by crossconvolving the relevant primaries thought to contain the stationary contributions for multiples.
  • Data- driven methodologies can generally handle complex geometries and need little or no information about the properties of the subsurface.
  • the model-based technologies are typically cost-effective compared to data-driven technologies, while the latter are typically more flexible.
  • Some model-driven methodologies require structural information, i.e., information about the subsurface structure, the determination of which is the reason for doing seismic exploration in the first place.
  • Other model-driven methodologies require the shape of the source wavelet that will not be a pure delta function because of the reverberations and frequency bandwidth limitation.
  • Some model-driven methodologies require both structural and source wavelet information while others use a matching filter to account for a distorted source wavelet.
  • SRME surface-related multiple elimination
  • SRME techniques are attractive solutions for predicting multiples in complex geologic settings, they do no require any a-priori knowledge of the subsurface (reflectivity, structures and velocities). However, these methods do require one shot location for each receiver position, and this is not the case for most three dimensional ("3D") acquisition geometries.
  • SRME methodologies are generally challenged by complex 3D multiples because of large shot spacing, narrow spread length and/or wide cable spacing. The missing data can be interpolated or extrapolated from the existing data, but interpolation or extrapolation has trouble with aliased seismic data caused by the large shot and/or receiver spacing. Advanced interpolation or extrapolation methods can also be difficult to implement and expensive.
  • a common cause of these complex 3D data that challenge 3D SRME methods is rugosity on the top of salt. But, any type of complex overburden can cause complex 3D seismic data that is hard to interpolate.
  • Another data-driven methodology utilizes predictive deconvolution which is a filtering method that assumes that multiples are periodic while primaries are not. This assumption is usually met for data from water depths less than 500 msec (approximately 1 ,200 feet) and approximately layered subsurface geology.
  • velocity-filtering methods in areas of water depths greater than 500 msec where the velocity difference between primaries and multiples are significant, velocity-filtering methods (as opposed to predictive methods) such as tctu-p and f-k filtering can be used, where the variable / represents frequency, k represents the wavenumber, p represents the ray parameter, and tau represents the zero offset intercept time.
  • filtering methods generally require determination, or at least an educated guess, of the apparent wave propagation velocities in the subsurface media through which the reflected seismic waves pass in their journey from the seismic source to a receiver. These velocities can differ significantly due to the combination of the variations of the subsurface structure and rock properties.
  • predictive deconvolution often leads to inadvertent damage to the primaries due to the difficulty in separating the multiples and primaries.
  • predictive deconvolution often fails to take into account the nonlinear factor in the reflectivity, which are generally caused by peg-leg multiples.
  • One prior art method which has extended predictive deconvolution for applications in deep water has utilized beam techniques. That method applies local slant stacking (or other dip-discriminating methods) to the data to decompose the recorded wavefields into beam components. These components travel approximately along raypaths. Simple raytracing within the water layer describes the long-period reverberations and relate primary and multiple events occurring in the beam components of the wavefield. Based on the information from the raytracing the time series of the beam-component of the primary can be shifted according to raytraced traveltimes and then analyzed with a multi-channel prediction filter. The predicted time series is considered as multiple energies and is removed from the beam components of the original data after a multi-channel matched filtering. Fig.
  • the prior art method includes initializing an earth model 4 which relates to a geological region of interest, and initiating a beam dataset 6 that has been determined from seismic data of the geological region of interest.
  • the prior art method further includes a series of loops wherein an input beam 8, a multiple-generating surface 10 and a time gate 12 are selected.
  • One or more time gates (or windows) are selected to ensure that the signal within each gate is stationary.
  • Trial rays are then "sprayed" from a detector location 14 and a stationary pegleg is determined 16.
  • the stationary pegleg is the pegleg that satisfies the Snell's law for reflection at the multiple-generating surface.
  • a primary beam corresponding to the stationary pegleg is obtained 18, and the primary beam is transformed into a predicted multiples beam 20.
  • the predicted multiples beam is then deconvolved to remove multiples which are present in the input beam 22.
  • the prior beam techniques assume that there is a single dominant multiple-generating surface 10 and the predicted multiples beams are related only to this multiple-generating surface 10 and do not contain predicted multiples from other multiple-generating surfaces.
  • the current invention improves prior art beam techniques to incorporate predicted multiples beams from multiple-generating surfaces that were not explicitly utilized to determine to stationary peglegs.
  • the present invention overcomes the above-described and other shortcomings of the prior art by providing a novel and improved method of predicting multiples based upon primaries which combines features from both model-driven and data-driven methodologies. It is accomplished by determining a model-driven stationary prediction based on the earth model and augmenting that prediction by a data-driven prediction around the stationary prediction. It should be appreciated that the model-driven stationary prediction can be replaced by an a-priori determination of stationary predictions, such as assuming a layered model.
  • One embodiment of the present invention includes a method for generating a prediction model of multiply-reflected, surface-related seismic waves which includes initializing an earth model related to a geological volume and selecting a beam dataset derived from seismic data related to the geological volume. The method also includes selecting an input beam from the beam dataset; a multiple-generating surface from the earth model, and a time gate. A stationary pegleg is determined utilizing the input beam, the multiple generating surface and the time gate. A primary beam which corresponds to the stationary pegleg is then obtained, and a modeled pegleg beam related to the primary beam is determined. The modeled pegleg beam is convolved with the primary beam to generate a convolved multiples beam. The convolved multiples beam is compared with the input beam to remove the multiples in the input beam by matched filtering.
  • the convolved multiples beam is utilized to provide a more accurate method of predicting and removing multiples than prior art methods.
  • the step of convolving the modeled pegleg beam with the primary beam is not included in the prior art methods, and enables the present invention to more accurately predict multiples.
  • convolving the primary beam and the modeled pegleg beam to obtain the convolved multiples beam includes transforming the primary' beam into a predicted multiples beam by traveltime shifting, and convolving the predicted multiples beam with a segment of the modeled pegleg beam to obtain a convolved multiples beam, the segment of the modeled pegleg beam starting at the multiple-generating surface and ending at a detector location.
  • Convolution is known in the art. In general, it is a mathematical operation on two functions that represents the process of linear filtering. Convolution can be applied to any two functions of time or space (or other variables) to yield a third function, the output of the convolution. Although the mathematical definition is symmetric with respect to the two input functions, it is common in signal processing to say that one of the functions is a filter acting on the other function. The response of many physical systems can be represented mathematically by a convolution. For example, a convolution is commonly used to model the filtering of seismic energy by the various rock layers in the Earth.
  • One embodiment of the present invention determines the stationary pegleg utilizing rays which are sprayed from a detector location, the detector location being based on the input beam, and the stationary pegleg is related to one of the rays sprayed from the detector location.
  • time gate is used herein to describe the entire beam or one of a plurality of segments included in the beam.
  • Beams in general are defined as energy components that are partially localized in space and dip.
  • Some examples of beams are Gaussian beams and other non-Gaussian beams, such as beams with complex rays, controlled beams and beams as a finite-difference solution of some version of the wave equation.
  • a segment of a modeled pegleg beam is convolved with the predicted multiples beam to generate a convolved multiples beam.
  • the segment of the modeled pegleg beam starts at the multiple-generating surface and involves a time interval that is either explicitly specified or determined by another horizon in the model.
  • the convolved multiples beam can be obtained by directly convolving the time sequences of the predicted multiples and the predicted modeled pegleg beams.
  • the present invention is intended to be used with a system which includes, in general, a computer configuration including at least one processor, at least one memory device for storing program code or other data, a video monitor or other display device (i.e., a liquid crystal display) and at least one input device.
  • the processor is preferably a microprocessor or microcontroller-based platform which is capable of displaying images and processing complex mathematical algorithms.
  • the memory device can include random access memory (RAM) for storing event or other data generated or used during a particular process associated with the present invention.
  • the memory device can also include read only memory (ROM) for storing the program code for the controls and processes of the present invention.
  • Fig. 1 illustrates a flow chart of a prior art method of attenuating multiples in seismic data
  • Fig. 2 illustrates a flow chart of one embodiment of the present invention for attenuating multiples in seismic data
  • Fig. 3 illustrates one embodiment of the present invention which includes transforming a primary beam into a predicted multiples beam
  • Fig. 4 illustrates one embodiment of the present invention wherein rays are sprayed from a detector location and are reflected off the water bottom surface of a marine environment;
  • Fig. 5 illustrates one embodiment of the present invention wherein a primary beam corresponding to a given pegleg is determined through raytracing
  • Fig. 6 illustrates a method utilized by one embodiment of the present invention to determine the stationary pegleg
  • Fig. 7 illustrates a schematic drawing of the above method of determining the stationary pegleg used by one embodiment of the present invention
  • Fig. 8 illustrates one embodiment of the present invention wherein the primary beam is being convolved with the pegleg beam
  • Fig. 9 illustrates a schematic drawing of a single-beam deconvolution that is utilized by one embodiment of the present invention
  • Fig. 10 illustrates a flow chart of another embodiment of the present invention for attenuating multiples in seismic data
  • Fig. 1 1 illustrates a flow chart of a further embodiment of the present invention for attenuating multiples in seismic data
  • Fig. 12 illustrates a flow chart of one embodiment of the present invention for attenuating multiples in seismic data
  • Fig. 13 illustrates a flow chart of another embodiment of the present invention for attenuating multiples in seismic data
  • Fig. 14 illustrates one embodiment of the invention wherein a primary beam corresponding to a pegleg is determined without raytracing.
  • the embodiment includes initializing an earth model that corresponds to a geological area of interest 32 and selecting a beam dataset derived from seismic data collected in the geological area of interest 34.
  • An input beam from the beam dataset 36, a multiple- generating surface from the earth model 38, and a time gate 40 are selected.
  • a stationary pegleg is determined utilizing the input beam, the multiple generating surface and the time gate 42.
  • a primary beam corresponding to the stationary pegleg is then obtained 44.
  • a modeled pegleg beam related to the primary beam is obtained 46.
  • the modeled pegleg beam is convolved with the primary beam to determine a convolved multiples beam 48.
  • the convolved multiples beam is deconvolved with the input beam to remove the multiples in the input beam 50.
  • the predicted primary beam can be directly convolved with the modeled pegleg beam.
  • rays are sprayed from a detector location and the stationary pegleg is determined based upon one of the rays.
  • the detector location in that embodiment is based upon the input beam.
  • a segment of the modeled pegleg beam is convolved with the primary beam to obtain a convolved multiples beam.
  • the segment of the modeled pegleg beam starts at the multiple-generating surface and involves a time interval that is either explicitly specified or is determined by another horizon in the model.
  • the convolved multiples beam is then deconvolved with the input beam to remove the multiples in the input beam.
  • prior art methods have used local slant stacking or other dip-discriminating methods for seismic traces to separate a recorded wavefield into beam components and those methods are known in the art.
  • the present invention utilizes local slant stacking to separate the recorded wavefield into components that are localized in both position and dip.
  • the local slant stacking described-above is used for 3D acquisition, the recorded energy cannot be completely steered into beams because the wavefield is not densely sampled along all recording directions, thus there is an issue of missing data or severe aliasing that is needed to accurately determine the raytrace traveltime T AB -
  • the local slant stacking is accomplished in the common-offset domain but not in the common-midpoint domain.
  • an assumption is made that stacking velocities describe the dip of primary events in common-midpoint (“CMP") gathers.
  • CMP common-midpoint
  • stacking velocities are a reasonable description of primaries from geological structures above the subsurface salt formations, which can reverberate in marine environments to become the strongest multiples.
  • Figs. 4, 5 and 6 illustrate the manner in which this embodiment of the present invention calculates the raytrace traveltime for the determination of the stationary pegleg.
  • An input beam is selected with a source location "S" 74 and a detector location “D” 76 as illustrated in Fig. 4.
  • a multiple generating surface is determined, in this embodiment it is the water bottom ("wb") 78.
  • Rays 80 are sprayed from the detector D 76 which are reflected back from the multiple-generating surface wb 78 and reach the free surface 82.
  • the angular interval for the rays 80 is predetermined to be dp x in the x (vertically downward from the free surface 82) direction and dp y in the y direction (perpendicular to the sheet).
  • An individual ray 84 is selected with a given ray parameter or detector raypath dip p d and an arrival location Q 86 of the ray 84 is determined at the free surface 82 as illustrated in Fig. 5.
  • the arrival direction (ray parameter p q' 88) of the ray 84 is determined at the arrival location Q 86.
  • the reflection at location Q 86 is determined and an outgoing ray parameter p q 90 is calculated.
  • the locations S 74 and Q 86 and their ray parameters p s 92 and p q 90 determine the primary corresponding to the pegleg with locations D 76 and Q 86 and ray parameters p d 84 and p q' 88.
  • Finding the stationary raypaths that describe the reflection at location Q 86 requires a raypath search.
  • the search is performed to find the separate reflections that occur at the various locations S 74, Q 86 and D 76.
  • the search compares the raytraced p h with the calculated p h ( p h being the offset dip).
  • the raytraced p h is calculated as;
  • t time
  • h the half offset
  • V the NMO velocities
  • T o zero-offset traveltime
  • V and T o are derived from the stacking velocities obtained from the earth model.
  • Fig. 6 illustrates the determination of p h where p h is the local slope 96 of the curve 98 in the h-t plot 94. The beam with the closest match between the raytraced p h and the calculated p h is selected as the stationary pegleg.
  • Fig. 7 illustrates a schematic of the method in which the stationary pegleg is determined by this embodiment of the present invention.
  • a midpoint dip p m is selected 102, and a trial direction p d is selected from the detector location D that corresponds to the p m is chosen 104.
  • the parameter p d is used to perform a raytrace from the detector location D to a point Q selected between the source location S and the detector location D 106.
  • the raytrace results in a raytraced offset dip p h which is compared to a calculated offset dip p h 108, if the raytraced p h closely matches the calculated p h then this process is complete and the beam corresponding to p m is selected as the stationary pegleg 1 10. If the raytraced p h does not match the calculated p h then another trial direction p d is chosen and the process is performed again 1 12 until a satisfactory match between the raytraced p h and the calculated p h is obtained.
  • This embodiment of the present invention also includes convolving the predicted multiples beam with a segment of the modeled beam to obtain the convolved multiples beam.
  • the predicted multiples beam B 1 14 is convolved with the a segment of the modeled pegleg beam C 1 16 starting at the multiple-generating surface 1 18.
  • the result of the convolution is beam E 120 which is the convolved multiples beam that is deconvolved with input beam A 122 to remove multiples that are present in the input beam A 122.
  • FIG. 9 A schematic illustration of this embodiment is provided in Fig. 9, wherein beam B 1 14 is utilized as the source side prediction 124, and beam C 1 16 is utilized for the detector side prediction 126.
  • This embodiment utilizes a Wiener Filter 128 and inputs from beams A 122, B 1 14 and C 116 to generate an estimation of the multiples present in beam A 122.
  • a Wiener filter is known in the art. In general, it is a causal filter which will transform an input into a desired output as closely as possible, subject to certain constraints. As one skilled in the art will appreciate there are other filters or means that can perform this particular function and they are intended to be within the scope of the present invention.
  • a 122 have been determined, the multiples are then removed 130 from Beam A 122.
  • One embodiment of the present invention is illustrated in Fig. 10 wherein an earth model is initialized 134 which correlates to a specific geological region of interest.
  • a beam dataset 136 that has been determined from seismic data of the geological region of interest is also initiated.
  • This embodiment of the present invention includes a series of loops wherein an input beam 138, a multiple-generating surface 140 and a time gate 142 are selected. Rays are sprayed from a detector location that is based on the input beam 144 and a stationary pegleg is selected from one of the rays 146.
  • a primary beam corresponding to the pegleg is obtained 148, and the primary beam is transformed into a predicted multiples beam 150 by a shift corresponding to the traveltime of the ray corresponding to the stationary pegleg.
  • a modeled pegleg beam related to the predicted multiples beam is then generated 152.
  • a segment of the modeled pegleg beam which starts at the multiple-generating surface is convolved with the predicted multiples beam to obtain a convolved multiples beam 154.
  • the convolved multiples beam is then either accumulated or deconvolved with the input beam to remove the multiples in the input beam 162.
  • the accumulated beams can be used to reconstruct the multiple prediction as a seismic trace or be used to deconvolve with the input beam at a later time.
  • This embodiment of the present invention allows a number of different points after convolution to either accumulate the convolved multiples beam or deconvolve the convolved multiples beam to remove the multiples in the input beam 162. Depending on the data being processed, those steps can occur before the end of the For loop for selecting the time gate 142-156 or immediately after that loop 142-156. Those steps can also occur after For loop for selecting the multiple-generating surface 140-158 or after the For loop for selecting the input beam 138- 160.
  • FIG. 1 Another embodiment of the present invention is illustrated in Fig. 1 1, that embodiment includes choosing a pegleg which is in a narrow range around a selected stationary pegleg 176.
  • a primary beam which corresponds to the chosen pegleg is obtained 178, and that primary beam is transformed into a predicted multiples beam 180.
  • the modeled pegleg beam which relates to the predicted multiples beam is obtained 182.
  • the predicted multiples beam is then convolved with a segment of the modeled pegleg beam to obtain a convolved multiples beam.
  • the segment of the modeled pegleg beam utilized by this step starts at the multiple-generating surface and ends at the detector location.
  • the predictions or convolved multiples beam for the narrow range around the stationary are stacked.
  • the steps of accumulating the convolved multiples beam or deconvolving the convolved multiples beam to remove the multiples in the input beam 194 can occur in this embodiment after the For loop for selecting the pegleg 176-186, the For loop for selecting the time gate 172-188, the For loop for selecting the multiple-generating surface 170-190, or the For loop for selecting the input beam 168-192.
  • FIG. 12 Another embodiment of the present invention is illustrated in Fig. 12, wherein the pegleg that is selected 208 is not tied to a stationary pegleg. That pegleg is used to obtain a corresponding primary beam 210.
  • the primary beam is transformed into a predicted multiples beam 212 and a modeled pegleg beam is obtained 214.
  • a segment of the modeled pegleg beam that starts at the multiple-generating surface is convolved with the predicted multiples beam to obtain a convolved multiples beam 216.
  • the convolved multiples beams that are generated within the For loop pegleg 208-218 are accumulated.
  • the steps of accumulating the convolved multiples beam or deconvolving the convolved multiples beam to remove the multiples in the input beam 226 can occur in this embodiment after the For loop for selecting the pegleg 208-218, the For loop for selecting the time gate 204-220, the For loop for selecting the multiple-generating surface 202-222, or the For loop for selecting the input beam 200-224.
  • Both the stationary pegleg and the pegleg can be determined for a variety of multiples.
  • One embodiment of the present invention determines the stationary pegleg or pegleg for a source-side multiple.
  • Another embodiment determines the stationary pegleg or pegleg for a detector-side multiple.
  • a further embodiment determines the stationary pegleg or pegleg for both source-side and the detector-side multiples.
  • an input beam is selected 230 from an opened beam dataset 228 as illustrated in Fig. 13.
  • a time gate is selected 232 and the pegleg selection 236 is aided by an a-priori determination of stationary predictions 234, such as assuming a layered earth model or determining the areal coverage of the source and detector locations.
  • a pegleg is selected 236 and a primary beam corresponding to the pegleg is obtained 238.
  • a modeled pegleg beam corresponding to the primary beam is calculated 240, and the primary beam is convolved with the modeled pegleg beam 242.
  • the above-described steps can either be repeated, 230- 248, 232-246, 336-244 with the convolved multiples beam being accumulated or the convolved multiples beam is deconvolved with the input beam to remove the multiples from the input beam 250.
  • the above-described embodiment convolves a plurality of primary and pegleg beams for a range of locations Q 252 and a range of ray parameter values P q 254 and P s 256 as illustrated in Fig. 14.
  • the range of Q 252 and the ranges P q 254 and P s 256 (or P d 258) is predetermined from an analysis of the input data.
  • the ranges are determined so that enough beams are included to contain the stationary contribution for the modeled multiples that is amplified after summation of the stationary contribution with other non-stationary contributions.
  • the primary beam can then be convolved with the pegleg beam.
  • the primary beam and the pegleg beam are convolved for a range of ray parameter P q 254 values.
  • Figs. 10-13 incorporate loops which illustrate that certain steps of those embodiments can be repeated depending on the data which is being processed.
  • Certain embodiments of the present invention described-above include spraying rays from the detector location to select the stationary pegleg or pegleg, it should be understood that there are another means of selecting the stationary pegleg or pegleg, for example, random selection, and those means are considered to within the scope of the present invention.

Abstract

The present invention incorporates the use of model-driven and data-driven methodologies to attenuate multiples in seismic data utilizing a prediction model which includes multiply-reflected, surface-related seismic waves. The present invention includes beam techniques and convolving a predicted multiples beam with a segment of a modeled pegleg beam to obtain a convolved multiples beam. The convolved multiples beam can then be deconvolved to attenuate the multiples that are present in the original input beam.

Description

METHOD FOR IDENTIFYING AND REMOVING MULTIPLES FOR
IMAGING WITH BEAMS
BACKGROUND OF THE INVENTION
This invention relates to seismic exploration and processing, and more specifically to imaging with beams and a method to predict multiples based on primaries by combining both model-driven and data-driven methodologies.
In the petroleum industry, seismic prospecting techniques are commonly used to aid in the search for and the evaluation of subterranean hydrocarbon deposits. In seismic prospecting, one or more sources of seismic energy emit waves into a subsurface region of interest, such as a geologic formation. These waves enter the formation and may be scattered, e.g., by reflection or refraction, by subsurface seismic reflectors (i.e., interfaces between underground formations having different elastic properties). The reflected signals are sampled or measured by one or more receivers, and the resultant data are recorded. The recorded samples may be referred to as seismic data or a set of "seismic traces". The seismic data may be analyzed to extract details of the structure and properties of the region of the earth being explored.
Seismic prospecting consists of three separate stages: data acquisition, data processing and data interpretation. The success of a seismic prospecting operation depends on satisfactory completion of all three stages.
In general, the purpose of seismic exploration is to map or image a portion of the subsurface of the earth (a formation) by transmitting energy down into the ground and recording the "reflections" or "echoes" that return from the rock layers below. The energy transmitted into the formation is typically sound and shear wave energy. The downward-propagating sound energy may originate from various sources, such as explosions or seismic vibrators on land or air guns in marine environments. Seismic exploration typically uses one or more energy sources and typically a large number of sensors or detectors. The sensors that may be used to detect the returning seismic energy are usually geophones (land surveys) or hydrophones (marine surveys). During a surface seismic survey, the energy source may be positioned at one or more locations near the surface of the earth above a geologic structure or formation of interest, referred to as shotpoints. Each time the source is activated, the source generates a seismic signal that travels downward through the earth and is at least partially reflected from discontinuities of various types in the subsurface, including reflections from "rock layer" boundaries. In general, a partial reflection of seismic signals may occur each time there is a change in the elastic properties of the subsurface materials. Reflected seismic signals are transmitted back to the surface of the earth, where they are recorded as a function of traveltime at a number of locations. The returning signals are digitized and recorded as a function of time (amplitude vs.. time).
One prevalent issue with the seismic energy recorded by the receivers during the data acquisition stage is that the seismic traces often contain both the desired seismic reflections (the "primary" reflections) and unwanted multiple reflections which can obscure or overwhelm the primary seismic reflections. A primary reflection is a sound wave that passes from the source to a receiver with a single reflection from a subsurface seismic reflector. A multiple reflection is a wave that has reflected at least three times (up, down and back up again) before being detected by a receiver. Depending on their time delay from primary events with which they are associated, multiples are commonly characterized as short-path, implying that they interfere with their own primary reflection, or long-path, where they appear as separate events.
There are also a variety of multiple events which are well known in the art. There are signals which are "trapped" in the water layer between two strong reflectors, the free surface and the bottom of the water layer. There are "peg-leg" multiple events, which are reflections that are characterized by an additional roundtrip through the water layer just after emission or just before detection. There are "remaining" surface-related multiple events, where the first and last upward reflections are below the first (water) layer, and there is at least one reflection at the free surface in between. There are also "interbed" multiples which has a downward reflection occurring from a subsurface reflector. In most cases, multiples do not contain any useful information that is not more easily extracted from primaries. Moreover, water-bottom multiples have been recognized as the most serious noise problem in seismic data processing in many offshore areas. Multiples can severely mask primary reflection events for structural imaging and contaminate Amplitude vs. Offset ("AVO") information. For those reasons, removal of multiples, or at least attenuation of multiples is a necessary part of the seismic data processing stage in many environments, particularly in marine settings where multiples are especially strong relative to the primaries.
In the case of deep-water data, suppression of first-order and the next few orders of sea-bottom multiple and peg-leg reflections are of great importance. These rather strong multiples may have the same traveltime as the primary reflections of target reflectors.
There are several prior art methods to attenuate multiples depending on the attributes of the multiples utilized. One class of multiple attenuation methods is the predictive methods where the multiples are predicted from their respective primaries. Prior art predictive multiple attenuation techniques can be generally divided into two categories; model-driven methodologies and data-driven methodologies. Model-driven methodologies generally use an earth model and the recorded data to predict or simulate multiples utilizing the estimated sea-bottom reflectivity function and calculated amplitude functions to model water-layer multiple reflections, those predicted multiples are then subtracted from the original data. Other model-driven technologies utilize an earth model or reflectivity model to predict the stationary multiples. The data-driven methodologies exploit the fact that primaries and multiples are physically related through a convolutional relationship and predict multiples by crossconvolving the relevant primaries thought to contain the stationary contributions for multiples. Data- driven methodologies can generally handle complex geometries and need little or no information about the properties of the subsurface. The model-based technologies are typically cost-effective compared to data-driven technologies, while the latter are typically more flexible. Some model-driven methodologies require structural information, i.e., information about the subsurface structure, the determination of which is the reason for doing seismic exploration in the first place. Other model-driven methodologies require the shape of the source wavelet that will not be a pure delta function because of the reverberations and frequency bandwidth limitation. Some model-driven methodologies require both structural and source wavelet information while others use a matching filter to account for a distorted source wavelet.
Data-driven methodologies rely on the predictability of multiples from primary components. In effect, that methodology utilizes existing seismic data to generate multiples and those generated multiples are then subtracted from the existing data. One such prior art methodology that is data-driven is known as "surface-related multiple elimination" or "SRME". In brief, this method operates by utilizing the existing data to create a dataset that contains only predictions of the multiples that are present in the data. Specifically, the method seeks to predict the seismic expression of multiples, and after adaptation to the existing multiples in the data the predicted multiples are subtracted from the original data leaving behind (at least theoretically) only the primary energy.
Data-driven SRME techniques are attractive solutions for predicting multiples in complex geologic settings, they do no require any a-priori knowledge of the subsurface (reflectivity, structures and velocities). However, these methods do require one shot location for each receiver position, and this is not the case for most three dimensional ("3D") acquisition geometries. SRME methodologies are generally challenged by complex 3D multiples because of large shot spacing, narrow spread length and/or wide cable spacing. The missing data can be interpolated or extrapolated from the existing data, but interpolation or extrapolation has trouble with aliased seismic data caused by the large shot and/or receiver spacing. Advanced interpolation or extrapolation methods can also be difficult to implement and expensive. A common cause of these complex 3D data that challenge 3D SRME methods is rugosity on the top of salt. But, any type of complex overburden can cause complex 3D seismic data that is hard to interpolate. Another data-driven methodology utilizes predictive deconvolution which is a filtering method that assumes that multiples are periodic while primaries are not. This assumption is usually met for data from water depths less than 500 msec (approximately 1 ,200 feet) and approximately layered subsurface geology. In areas of water depths greater than 500 msec where the velocity difference between primaries and multiples are significant, velocity-filtering methods (as opposed to predictive methods) such as tctu-p and f-k filtering can be used, where the variable / represents frequency, k represents the wavenumber, p represents the ray parameter, and tau represents the zero offset intercept time.
However, filtering methods generally require determination, or at least an educated guess, of the apparent wave propagation velocities in the subsurface media through which the reflected seismic waves pass in their journey from the seismic source to a receiver. These velocities can differ significantly due to the combination of the variations of the subsurface structure and rock properties. In addition, predictive deconvolution often leads to inadvertent damage to the primaries due to the difficulty in separating the multiples and primaries. Moreover, predictive deconvolution often fails to take into account the nonlinear factor in the reflectivity, which are generally caused by peg-leg multiples.
One prior art method which has extended predictive deconvolution for applications in deep water has utilized beam techniques. That method applies local slant stacking (or other dip-discriminating methods) to the data to decompose the recorded wavefields into beam components. These components travel approximately along raypaths. Simple raytracing within the water layer describes the long-period reverberations and relate primary and multiple events occurring in the beam components of the wavefield. Based on the information from the raytracing the time series of the beam-component of the primary can be shifted according to raytraced traveltimes and then analyzed with a multi-channel prediction filter. The predicted time series is considered as multiple energies and is removed from the beam components of the original data after a multi-channel matched filtering. Fig. 1 illustrates a flowchart for one example of a prior art method wherein deconvolution is utilized with a beam technique for attenuating multiples. The prior art method includes initializing an earth model 4 which relates to a geological region of interest, and initiating a beam dataset 6 that has been determined from seismic data of the geological region of interest. The prior art method further includes a series of loops wherein an input beam 8, a multiple-generating surface 10 and a time gate 12 are selected. One or more time gates (or windows) are selected to ensure that the signal within each gate is stationary. Trial rays are then "sprayed" from a detector location 14 and a stationary pegleg is determined 16. The stationary pegleg is the pegleg that satisfies the Snell's law for reflection at the multiple-generating surface. A primary beam corresponding to the stationary pegleg is obtained 18, and the primary beam is transformed into a predicted multiples beam 20. The predicted multiples beam is then deconvolved to remove multiples which are present in the input beam 22.
While the beam techniques have improved prior art multiple attenuation techniques, there is still a need for an improved method which provides a more accurate prediction of multiples and therefore allows for more accurate subtraction of those multiples from the data. The prior beam techniques assume that there is a single dominant multiple-generating surface 10 and the predicted multiples beams are related only to this multiple-generating surface 10 and do not contain predicted multiples from other multiple-generating surfaces. The current invention improves prior art beam techniques to incorporate predicted multiples beams from multiple-generating surfaces that were not explicitly utilized to determine to stationary peglegs.
SUMMARY OF THE INVENTION
The present invention overcomes the above-described and other shortcomings of the prior art by providing a novel and improved method of predicting multiples based upon primaries which combines features from both model-driven and data-driven methodologies. It is accomplished by determining a model-driven stationary prediction based on the earth model and augmenting that prediction by a data-driven prediction around the stationary prediction. It should be appreciated that the model-driven stationary prediction can be replaced by an a-priori determination of stationary predictions, such as assuming a layered model.
One embodiment of the present invention includes a method for generating a prediction model of multiply-reflected, surface-related seismic waves which includes initializing an earth model related to a geological volume and selecting a beam dataset derived from seismic data related to the geological volume. The method also includes selecting an input beam from the beam dataset; a multiple-generating surface from the earth model, and a time gate. A stationary pegleg is determined utilizing the input beam, the multiple generating surface and the time gate. A primary beam which corresponds to the stationary pegleg is then obtained, and a modeled pegleg beam related to the primary beam is determined. The modeled pegleg beam is convolved with the primary beam to generate a convolved multiples beam. The convolved multiples beam is compared with the input beam to remove the multiples in the input beam by matched filtering.
The convolved multiples beam is utilized to provide a more accurate method of predicting and removing multiples than prior art methods. The step of convolving the modeled pegleg beam with the primary beam is not included in the prior art methods, and enables the present invention to more accurately predict multiples. For example, in one embodiment of the present invention convolving the primary beam and the modeled pegleg beam to obtain the convolved multiples beam includes transforming the primary' beam into a predicted multiples beam by traveltime shifting, and convolving the predicted multiples beam with a segment of the modeled pegleg beam to obtain a convolved multiples beam, the segment of the modeled pegleg beam starting at the multiple-generating surface and ending at a detector location.
"Convolution" is known in the art. In general, it is a mathematical operation on two functions that represents the process of linear filtering. Convolution can be applied to any two functions of time or space (or other variables) to yield a third function, the output of the convolution. Although the mathematical definition is symmetric with respect to the two input functions, it is common in signal processing to say that one of the functions is a filter acting on the other function. The response of many physical systems can be represented mathematically by a convolution. For example, a convolution is commonly used to model the filtering of seismic energy by the various rock layers in the Earth.
One embodiment of the present invention determines the stationary pegleg utilizing rays which are sprayed from a detector location, the detector location being based on the input beam, and the stationary pegleg is related to one of the rays sprayed from the detector location.
As one skilled in the art will appreciate, the phrase "time gate" is used herein to describe the entire beam or one of a plurality of segments included in the beam.
It should also be appreciated that utilizing beam techniques to process seismic data are well known in the art, and those techniques are within the scope of the present invention. Beams in general are defined as energy components that are partially localized in space and dip. Some examples of beams are Gaussian beams and other non-Gaussian beams, such as beams with complex rays, controlled beams and beams as a finite-difference solution of some version of the wave equation.
In one embodiment of the present invention, a segment of a modeled pegleg beam is convolved with the predicted multiples beam to generate a convolved multiples beam. The segment of the modeled pegleg beam starts at the multiple-generating surface and involves a time interval that is either explicitly specified or determined by another horizon in the model.
In another embodiment of the present invention, the convolved multiples beam can be obtained by directly convolving the time sequences of the predicted multiples and the predicted modeled pegleg beams.
It should also be appreciated that the present invention is intended to be used with a system which includes, in general, a computer configuration including at least one processor, at least one memory device for storing program code or other data, a video monitor or other display device (i.e., a liquid crystal display) and at least one input device. The processor is preferably a microprocessor or microcontroller-based platform which is capable of displaying images and processing complex mathematical algorithms. The memory device can include random access memory (RAM) for storing event or other data generated or used during a particular process associated with the present invention. The memory device can also include read only memory (ROM) for storing the program code for the controls and processes of the present invention.
Additional features and advantages of the present invention are described in, and will be apparent from, the following Detailed Description of the Invention and the Figures.
BRIEF DESCRIPTION OF THE DRAWINGS
These and other objects, features and advantages of the present invention will become better understood with regard to the following description, pending claims and accompanying drawings where:
Fig. 1 illustrates a flow chart of a prior art method of attenuating multiples in seismic data;
Fig. 2 illustrates a flow chart of one embodiment of the present invention for attenuating multiples in seismic data;
Fig. 3 illustrates one embodiment of the present invention which includes transforming a primary beam into a predicted multiples beam;
Fig. 4 illustrates one embodiment of the present invention wherein rays are sprayed from a detector location and are reflected off the water bottom surface of a marine environment;
Fig. 5 illustrates one embodiment of the present invention wherein a primary beam corresponding to a given pegleg is determined through raytracing;
Fig. 6 illustrates a method utilized by one embodiment of the present invention to determine the stationary pegleg;
Fig. 7 illustrates a schematic drawing of the above method of determining the stationary pegleg used by one embodiment of the present invention;
Fig. 8 illustrates one embodiment of the present invention wherein the primary beam is being convolved with the pegleg beam;
Fig. 9 illustrates a schematic drawing of a single-beam deconvolution that is utilized by one embodiment of the present invention; Fig. 10 illustrates a flow chart of another embodiment of the present invention for attenuating multiples in seismic data;
Fig. 1 1 illustrates a flow chart of a further embodiment of the present invention for attenuating multiples in seismic data;
Fig. 12 illustrates a flow chart of one embodiment of the present invention for attenuating multiples in seismic data;
Fig. 13 illustrates a flow chart of another embodiment of the present invention for attenuating multiples in seismic data; and
Fig. 14 illustrates one embodiment of the invention wherein a primary beam corresponding to a pegleg is determined without raytracing.
DETAILED DESCRIPTION OF THE INVENTION
While this invention is susceptible to embodiments in many different forms, there are shown in the drawings, and will herein be described in detail, preferred embodiments of the invention with the understanding that the present disclosure is to be considered as an exemplification of the principles of the invention and is not intended to limit the broad aspect of the invention to the embodiments illustrated.
One embodiment of the present invention 30 is illustrated in Fig. 2. The embodiment includes initializing an earth model that corresponds to a geological area of interest 32 and selecting a beam dataset derived from seismic data collected in the geological area of interest 34. An input beam from the beam dataset 36, a multiple- generating surface from the earth model 38, and a time gate 40 are selected. A stationary pegleg is determined utilizing the input beam, the multiple generating surface and the time gate 42. A primary beam corresponding to the stationary pegleg is then obtained 44. A modeled pegleg beam related to the primary beam is obtained 46. The modeled pegleg beam is convolved with the primary beam to determine a convolved multiples beam 48. The convolved multiples beam is deconvolved with the input beam to remove the multiples in the input beam 50.
In another embodiment of the present invention, the predicted primary beam can be directly convolved with the modeled pegleg beam.
In yet another embodiment of the present invention, rays are sprayed from a detector location and the stationary pegleg is determined based upon one of the rays. The detector location in that embodiment is based upon the input beam.
In a further embodiment of the present invention, a segment of the modeled pegleg beam is convolved with the primary beam to obtain a convolved multiples beam. The segment of the modeled pegleg beam starts at the multiple-generating surface and involves a time interval that is either explicitly specified or is determined by another horizon in the model. The convolved multiples beam is then deconvolved with the input beam to remove the multiples in the input beam. As described-above, prior art methods have used local slant stacking or other dip-discriminating methods for seismic traces to separate a recorded wavefield into beam components and those methods are known in the art. The present invention utilizes local slant stacking to separate the recorded wavefield into components that are localized in both position and dip. These components are what would be recorded at the center of a beam arrival 54, 56 at locations A 58 and B 60 as illustrated in Fig. 3. The beam energy 54 that arrives at location B 60 reverberates in the water layer 64 and is assumed to arrive as multiples 62 within a beam 56 recorded at location A 58. Shifting beam B 54 by the raytrace traveltime TAB 66 from location B 60 to location A 58 will line up the events in beam B 54 with the multiples 68, 70 in beam A 56. Once the events have been aligned and the multiples 68, 70 are identified, the multiples 68, 70 can then be removed.
When the local slant stacking described-above is used for 3D acquisition, the recorded energy cannot be completely steered into beams because the wavefield is not densely sampled along all recording directions, thus there is an issue of missing data or severe aliasing that is needed to accurately determine the raytrace traveltime TAB- For example, as one skilled in the art will appreciate, the local slant stacking is accomplished in the common-offset domain but not in the common-midpoint domain. In one embodiment of the present invention, an assumption is made that stacking velocities describe the dip of primary events in common-midpoint ("CMP") gathers. In general, stacking velocities are a reasonable description of primaries from geological structures above the subsurface salt formations, which can reverberate in marine environments to become the strongest multiples.
Figs. 4, 5 and 6 illustrate the manner in which this embodiment of the present invention calculates the raytrace traveltime for the determination of the stationary pegleg. An input beam is selected with a source location "S" 74 and a detector location "D" 76 as illustrated in Fig. 4. A multiple generating surface is determined, in this embodiment it is the water bottom ("wb") 78. Rays 80 are sprayed from the detector D 76 which are reflected back from the multiple-generating surface wb 78 and reach the free surface 82. The angular interval for the rays 80 is predetermined to be dpx in the x (vertically downward from the free surface 82) direction and dpy in the y direction (perpendicular to the sheet). An individual ray 84 is selected with a given ray parameter or detector raypath dip pd and an arrival location Q 86 of the ray 84 is determined at the free surface 82 as illustrated in Fig. 5. The arrival direction (ray parameter pq' 88) of the ray 84 is determined at the arrival location Q 86. The reflection at location Q 86 is determined and an outgoing ray parameter pq 90 is calculated. The locations S 74 and Q 86 and their ray parameters ps 92 and pq 90 determine the primary corresponding to the pegleg with locations D 76 and Q 86 and ray parameters pd 84 and pq' 88.
Finding the stationary raypaths that describe the reflection at location Q 86 requires a raypath search. The search is performed to find the separate reflections that occur at the various locations S 74, Q 86 and D 76. In this particular embodiment, the search compares the raytraced ph with the calculated ph ( ph being the offset dip). The raytraced ph is calculated as;
Figure imgf000016_0001
where pm (the midpoint dip) at a given location corresponds to a particular beam. The calculated/?/, is obtained from the normal moveout equation ("NMO"):
Figure imgf000016_0002
where: t is time; h is the half offset; V is the NMO velocities; and To is zero-offset traveltime.
V and To are derived from the stacking velocities obtained from the earth model. Fig. 6 illustrates the determination of ph where ph is the local slope 96 of the curve 98 in the h-t plot 94. The beam with the closest match between the raytraced ph and the calculated ph is selected as the stationary pegleg.
Fig. 7 illustrates a schematic of the method in which the stationary pegleg is determined by this embodiment of the present invention. A midpoint dip pm is selected 102, and a trial direction pd is selected from the detector location D that corresponds to the pm is chosen 104. The parameter pd is used to perform a raytrace from the detector location D to a point Q selected between the source location S and the detector location D 106. The raytrace results in a raytraced offset dip ph which is compared to a calculated offset dip ph 108, if the raytraced ph closely matches the calculated ph then this process is complete and the beam corresponding to pm is selected as the stationary pegleg 1 10. If the raytraced ph does not match the calculated ph then another trial direction pd is chosen and the process is performed again 1 12 until a satisfactory match between the raytraced ph and the calculated ph is obtained.
This embodiment of the present invention also includes convolving the predicted multiples beam with a segment of the modeled beam to obtain the convolved multiples beam. As illustrated in Fig. 8, the predicted multiples beam B 1 14 is convolved with the a segment of the modeled pegleg beam C 1 16 starting at the multiple-generating surface 1 18. The result of the convolution is beam E 120 which is the convolved multiples beam that is deconvolved with input beam A 122 to remove multiples that are present in the input beam A 122.
A schematic illustration of this embodiment is provided in Fig. 9, wherein beam B 1 14 is utilized as the source side prediction 124, and beam C 1 16 is utilized for the detector side prediction 126. This embodiment utilizes a Wiener Filter 128 and inputs from beams A 122, B 1 14 and C 116 to generate an estimation of the multiples present in beam A 122. A Wiener filter is known in the art. In general, it is a causal filter which will transform an input into a desired output as closely as possible, subject to certain constraints. As one skilled in the art will appreciate there are other filters or means that can perform this particular function and they are intended to be within the scope of the present invention. Once the multiples in Beam
A 122 have been determined, the multiples are then removed 130 from Beam A 122. One embodiment of the present invention is illustrated in Fig. 10 wherein an earth model is initialized 134 which correlates to a specific geological region of interest. A beam dataset 136 that has been determined from seismic data of the geological region of interest is also initiated. This embodiment of the present invention includes a series of loops wherein an input beam 138, a multiple-generating surface 140 and a time gate 142 are selected. Rays are sprayed from a detector location that is based on the input beam 144 and a stationary pegleg is selected from one of the rays 146. A primary beam corresponding to the pegleg is obtained 148, and the primary beam is transformed into a predicted multiples beam 150 by a shift corresponding to the traveltime of the ray corresponding to the stationary pegleg. A modeled pegleg beam related to the predicted multiples beam is then generated 152. A segment of the modeled pegleg beam which starts at the multiple-generating surface is convolved with the predicted multiples beam to obtain a convolved multiples beam 154. The convolved multiples beam is then either accumulated or deconvolved with the input beam to remove the multiples in the input beam 162. The accumulated beams can be used to reconstruct the multiple prediction as a seismic trace or be used to deconvolve with the input beam at a later time.
This embodiment of the present invention allows a number of different points after convolution to either accumulate the convolved multiples beam or deconvolve the convolved multiples beam to remove the multiples in the input beam 162. Depending on the data being processed, those steps can occur before the end of the For loop for selecting the time gate 142-156 or immediately after that loop 142-156. Those steps can also occur after For loop for selecting the multiple-generating surface 140-158 or after the For loop for selecting the input beam 138- 160.
Another embodiment of the present invention is illustrated in Fig. 1 1, that embodiment includes choosing a pegleg which is in a narrow range around a selected stationary pegleg 176. A primary beam which corresponds to the chosen pegleg is obtained 178, and that primary beam is transformed into a predicted multiples beam 180. The modeled pegleg beam which relates to the predicted multiples beam is obtained 182. The predicted multiples beam is then convolved with a segment of the modeled pegleg beam to obtain a convolved multiples beam. The segment of the modeled pegleg beam utilized by this step starts at the multiple-generating surface and ends at the detector location. Within the For loop for selecting the pegleg 176-186, the predictions or convolved multiples beam for the narrow range around the stationary are stacked. The steps of accumulating the convolved multiples beam or deconvolving the convolved multiples beam to remove the multiples in the input beam 194 can occur in this embodiment after the For loop for selecting the pegleg 176-186, the For loop for selecting the time gate 172-188, the For loop for selecting the multiple-generating surface 170-190, or the For loop for selecting the input beam 168-192.
Another embodiment of the present invention is illustrated in Fig. 12, wherein the pegleg that is selected 208 is not tied to a stationary pegleg. That pegleg is used to obtain a corresponding primary beam 210. The primary beam is transformed into a predicted multiples beam 212 and a modeled pegleg beam is obtained 214. A segment of the modeled pegleg beam that starts at the multiple-generating surface is convolved with the predicted multiples beam to obtain a convolved multiples beam 216. The convolved multiples beams that are generated within the For loop pegleg 208-218 are accumulated. In this embodiment, the steps of accumulating the convolved multiples beam or deconvolving the convolved multiples beam to remove the multiples in the input beam 226 can occur in this embodiment after the For loop for selecting the pegleg 208-218, the For loop for selecting the time gate 204-220, the For loop for selecting the multiple-generating surface 202-222, or the For loop for selecting the input beam 200-224.
Both the stationary pegleg and the pegleg can be determined for a variety of multiples. One embodiment of the present invention determines the stationary pegleg or pegleg for a source-side multiple. Another embodiment determines the stationary pegleg or pegleg for a detector-side multiple. A further embodiment determines the stationary pegleg or pegleg for both source-side and the detector-side multiples.
As one skilled in the art will appreciate, there may be situations in which an earth model is not readily available, in such instances, the present invention is still able to predict and attenuate multiples. One embodiment of the present invention does not include the use of an earth model. In that embodiment, an input beam is selected 230 from an opened beam dataset 228 as illustrated in Fig. 13. A time gate is selected 232 and the pegleg selection 236 is aided by an a-priori determination of stationary predictions 234, such as assuming a layered earth model or determining the areal coverage of the source and detector locations. A pegleg is selected 236 and a primary beam corresponding to the pegleg is obtained 238. A modeled pegleg beam corresponding to the primary beam is calculated 240, and the primary beam is convolved with the modeled pegleg beam 242. The above-described steps can either be repeated, 230- 248, 232-246, 336-244 with the convolved multiples beam being accumulated or the convolved multiples beam is deconvolved with the input beam to remove the multiples from the input beam 250.
The above-described embodiment convolves a plurality of primary and pegleg beams for a range of locations Q 252 and a range of ray parameter values Pq 254 and Ps 256 as illustrated in Fig. 14. The range of Q 252 and the ranges Pq 254 and Ps 256 (or Pd 258) is predetermined from an analysis of the input data. The ranges are determined so that enough beams are included to contain the stationary contribution for the modeled multiples that is amplified after summation of the stationary contribution with other non-stationary contributions. In Fig. 14 a beam with a given ray parameter Pm = Ps 280 + Pd 258, a location Q 252, and ray parameters Pq 254 and Ps 256 (or Pd 258) are selected. The primary beam, Pn, i = Ps 256 + Pq 254, and the pegleg beam, Pm2 = Pd 258 + Pq 254, are determined. In this embodiment, there is an assumption made that the surface reflection at Q 252 is from Pq 254 to -Pq 260. The primary beam can then be convolved with the pegleg beam. In this embodiment, the primary beam and the pegleg beam are convolved for a range of location Q 252 values beginning at Q 252 = (S 262 + 264)/2 or some other predetermined locations. In addition, the primary beam and the pegleg beam are convolved for a range of ray parameter Pq 254 values.
As described-above, the embodiments of the present invention depicted in Figs. 10-13 incorporate loops which illustrate that certain steps of those embodiments can be repeated depending on the data which is being processed. Certain embodiments of the present invention described-above include spraying rays from the detector location to select the stationary pegleg or pegleg, it should be understood that there are another means of selecting the stationary pegleg or pegleg, for example, random selection, and those means are considered to within the scope of the present invention.
While in the foregoing specification this invention has been described in relation to certain preferred embodiments thereof, and many details have been set forth for purpose of illustration, it will be apparent to those skilled in the art that the invention is susceptible to alteration and that certain other details described herein can vary considerably without departing from the basic principles of the invention.

Claims

WHAT IS CLAIMED IS:
1. A method of attenuating multiples in multiply-reflected, surface-related seismic waves, the method comprising;
(a) initializing an earth model related to a geological volume;
(b) selecting a beam dataset derived from seismic data related to the geological volume;
(c) selecting an input beam from the beam dataset;
(d) selecting a multiple-generating surface from the earth model;
(e) selecting a time gate;
(f) determining a stationary pegleg utilizing the input beam, the multiple- generating surface and the time gate;
(g) obtaining a primary beam which corresponds to the stationary pegleg;
(h) obtaining a modeled pegleg beam that is related to the primary beam;
(i) convolving the primary beam and the modeled pegleg beam to obtain a convolved multiples beam; and
(j) deconvolving the convolved multiples beam with the input beam to remove the multiples in the input beam.
2. The method of claim 1 wherein the steps (c) to (j) are repeated for a plurality of beams, steps (d) to (j) are repeated for a plurality of multiple-generating surfaces, and steps (e) to (j) are repeated for a plurality of time gates.
3. The method of claim 2 which includes accumulating the convolved multiples beams before deconvolving the input beam with the convolved multiples beam to remove the multiples in the input beam.
4. The method of claim 1 wherein the step of selecting the time gate includes selecting the entire input beam.
5. The method of claim 1 wherein the step of selecting the time gate includes selecting a portion of the input beam.
6. The method of claim 1 wherein the stationary pegleg is determined for a source- side multiple.
7. The method of claim 1 wherein the stationary pegleg is determined for a detector-side multiple.
8. The method of claim 1 wherein the stationary pegleg is determined for both a source-side multiple and a detector-side multiple.
9. The method of claim 1 wherein determining the stationary pegleg includes spraying rays from a detector location, the detector location being based on the input beam, and determining a stationary pegleg that is related to one of the rays sprayed from the detector location.
10. The method of claim 1 wherein convolving the primary beam and the modeled pegleg beam to obtain the convolved multiples beam includes transforming the primary beam into a predicted multiples beam, and convolving the predicted multiples beam with a segment of the modeled pegleg beam to obtain a convolved multiples beam, the segment of the modeled pegleg beam starting at the multiple-generating surface and ending at a deeper surface or at the end of a pre-determined time interval.
1 1. The method of claim 1 wherein the convolved multiples beam is obtained as a direct convolution of the primary beam with the modeled pegleg beam.
12. The method of claim 1 wherein the beams are defined as energy components that are partially localized in space and dip.
13 The method of claim 1 wherein the beam dataset includes Gaussian beams.
14. The method of claim 1 wherein the beam dataset includes beams with complex rays.
15. The method of claim 1 wherein the beam dataset includes controlled beams.
16. The method of claim 1 wherein the beam dataset includes finite-difference beams.
17. A method of attenuating multiples in multiply-reflected, surface-related seismic waves, the method comprising;
(a) initializing an earth model related to a geological volume;
(b) selecting a beam dataset derived from seismic data related to the geological volume;
(c) selecting an input beam from the beam dataset;
(d) selecting a multiple-generating surface from the earth model;
(e) selecting a time gate;
(f) selecting a pegleg;
(g) obtaining a primary beam which corresponds to the pegleg;
(h) obtaining a modeled pegleg beam that is related to the primary beam; (i) convolving the primary beam and the pegleg beam to obtain a convolved multiples beam; and
(j) deconvolving the convolved multiples beam with the input beam to remove the multiples in the input beam.
18. The method of claim 17 wherein steps (c) to (j) are repeated for a plurality of beams, steps (d) to (j) are repeated for a plurality of multiple-generating surfaces; and steps (e) to (j) are repeated for a plurality of time gates.
19. The method of claim 18 which includes accumulating the convolved multiples beams before deconvolving the input beam with the convolved multiples beam to remove the multiples in the input beam.
20. The method of claim 17 wherein the step of selecting the time gate includes selecting the entire input beam.
21. The method of claim 17 wherein the step of selecting the time gate includes selecting a portion of the input beam.
22. The method of claim 17 wherein the pegleg is selected for a source-side multiple.
23. The method of claim 17 wherein the pegleg is selected for a detector-side multiple.
24. The method of claim 17 wherein the pegleg is selected for both a source-side multiple and a detector-side multiple.
25. The method of claim 17 wherein a stationary pegleg utilizing the input beam, the multiple-generating surface and the time gate is determined, and the selected pegleg is related to the stationary pegleg.
26. The method of claim 25 wherein the stationary pegleg is determined for a source-side multiple.
27. The method of claim 25 wherein the stationary pegleg is determined for a detector-side multiple.
28. The method of claim 25 wherein the stationary pegleg is determined for both a source-side multiple and a detector-side multiple.
29. The method of claim 25 wherein the pegleg is chosen by selecting a range of ray parameter values around a ray parameter value of the stationary pegleg.
30. The method of claim 25 wherein the pegleg is chosen by selecting a range of locations around an arrival location of the stationary pegleg.
31. The method of claim 17 wherein the beams are defined as energy components that are partially localized in space and dip.
32. The method of claim 17 wherein the beam dataset includes Gaussian beams.
33. The method of claim 17 wherein the beam dataset includes complex beams.
34. The method of claim 17 wherein the beam dataset includes controlled beams.
35. The method of claim 17 wherein the beam dataset includes finite-difference beams.
36. A method of attenuating multiples in multiply-reflected, surface-related seismic waves, the method comprising;
(a) selecting an input beam from a beam dataset; (b) selecting a time gate;
(c) selecting a pegleg;
(d) obtaining a primary beam which corresponds to the pegleg;
(e) obtaining a pegleg beam that is related to the primary beam;
(f) convolving the primary beam and the pegleg beam to obtain a convolved multiples beam; and
(g) deconvolving the convolved multiples beam with the input beam to remove the multiples in the input beam.
37. The method of claim 36 wherein steps (c) to (g) are repeated for a plurality of beams, steps (d) to (g) are repeated for a plurality of multiple-generating surfaces; and steps (e) to (g) are repeated for a plurality of time gates.
38. The method of claim 37 which includes accumulating the convolved multiples beams before deconvolving the input beam with the convolved multiples beam to remove the multiples in the input beam.
39. The method of claim 36 wherein the step of selecting the time gate includes selecting the entire input beam.
40. The method of claim 36 wherein the step of selecting the time gate includes selecting a portion of the input beam.
41. The method of claim 36 wherein the pegleg is selected for a source-side multiple.
42. The method of claim 36 wherein the pegleg is selected for a detector-side multiple.
43. The method of claim 36 wherein the pegleg is selected for both a source-side multiple and a detector-side multiple
44. The method of claim 36 wherein the beams are defined as energy components that are partially localized in space and dip.
45. The method of claim 36 wherein the beam dataset includes Gaussian beams.
46. The method of claim 36 wherein the beam dataset includes beams with complex rays.
47. The method of claim 36 wherein the beam dataset includes controlled beams.
48. The method of claim 36 wherein the beam dataset includes finite-difference beams.
PCT/US2008/062877 2007-05-22 2008-05-07 Method for identifying and removing multiples for imaging with beams WO2008147657A2 (en)

Priority Applications (6)

Application Number Priority Date Filing Date Title
BRPI0811841-8A2A BRPI0811841A2 (en) 2007-05-22 2008-05-07 METHOD FOR ATTENUATING MULTIPLES IN SURFACE-RELATED AND MULTIPLE REFLECTED Seismic Waves
CA2687095A CA2687095C (en) 2007-05-22 2008-05-07 Method for identifying and removing multiples for imaging with beams
CN200880021551.9A CN101688922B (en) 2007-05-22 2008-05-07 Method for identifying and removing multiples for imaging with beams
EA200971083A EA017379B1 (en) 2007-05-22 2008-05-07 Method for identifying and removing multiples for imaging with beams
AU2008257075A AU2008257075B2 (en) 2007-05-22 2008-05-07 Method for identifying and removing multiples for imaging with beams
NO20093563A NO20093563L (en) 2007-05-22 2009-12-21 Method for identifying and removing multiples for imaging with rays

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
US11/751,975 US7715986B2 (en) 2007-05-22 2007-05-22 Method for identifying and removing multiples for imaging with beams
US11/751,975 2007-05-22

Publications (2)

Publication Number Publication Date
WO2008147657A2 true WO2008147657A2 (en) 2008-12-04
WO2008147657A3 WO2008147657A3 (en) 2009-03-26

Family

ID=39730746

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/US2008/062877 WO2008147657A2 (en) 2007-05-22 2008-05-07 Method for identifying and removing multiples for imaging with beams

Country Status (9)

Country Link
US (1) US7715986B2 (en)
EP (1) EP1998189A3 (en)
CN (1) CN101688922B (en)
AU (1) AU2008257075B2 (en)
BR (1) BRPI0811841A2 (en)
CA (1) CA2687095C (en)
EA (1) EA017379B1 (en)
NO (1) NO20093563L (en)
WO (1) WO2008147657A2 (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102478664A (en) * 2010-11-23 2012-05-30 中国石油天然气集团公司 Spatial sampling interval determining method without polluting effective signals

Families Citing this family (38)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8326543B2 (en) * 2009-11-20 2012-12-04 Chevron U.S.A. Inc. System and method for attenuating multiples in seismic data that accounts for beam azimuth
US8660798B2 (en) * 2010-02-25 2014-02-25 Chevron U.S.A. Inc. System and method for attenuating aliasing in seismic data caused by acquisition geometry
US8694299B2 (en) 2010-05-07 2014-04-08 Exxonmobil Upstream Research Company Artifact reduction in iterative inversion of geophysical data
SG188191A1 (en) 2010-09-27 2013-04-30 Exxonmobil Upstream Res Co Simultaneous source encoding and source separation as a practical solution for full wavefield inversion
US8437998B2 (en) 2010-09-27 2013-05-07 Exxonmobil Upstream Research Company Hybrid method for full waveform inversion using simultaneous and sequential source method
CN102169189B (en) * 2011-01-19 2015-01-14 中国海洋石油总公司 Method for eliminating deep water interlamination multiple
BR112013018994A2 (en) 2011-03-30 2017-02-21 Exxonmobil Upstream Res Co full wave field inversion convergence rate employing spectral conformation
WO2012134609A1 (en) 2011-03-31 2012-10-04 Exxonmobil Upstream Research Company Method of wavelet estimation and multiple prediction in full wavefield inversion
US9176930B2 (en) 2011-11-29 2015-11-03 Exxonmobil Upstream Research Company Methods for approximating hessian times vector operation in full wavefield inversion
CN104204857B (en) 2012-03-08 2017-12-08 埃克森美孚上游研究公司 Orthogonal source and receiver coding
CA2892041C (en) 2012-11-28 2018-02-27 Exxonmobil Upstream Research Company Reflection seismic data q tomography
WO2014177522A2 (en) * 2013-04-29 2014-11-06 Cgg Services Sa Device and method for wave-field reconstruction
MY169125A (en) 2013-05-24 2019-02-18 Exxonmobil Upstream Res Co Multi-parameter inversion through offset dependent elastic fwi
US10459117B2 (en) 2013-06-03 2019-10-29 Exxonmobil Upstream Research Company Extended subspace method for cross-talk mitigation in multi-parameter inversion
US10054704B2 (en) 2013-07-01 2018-08-21 Westerngeco L.L.C. Predicting multiples in survey data
US9702998B2 (en) 2013-07-08 2017-07-11 Exxonmobil Upstream Research Company Full-wavefield inversion of primaries and multiples in marine environment
US9588244B2 (en) * 2013-07-25 2017-03-07 Chevron U.S.A. Inc. Predicting interbed multiples in seismic data using beam decomposition
ES2671789T3 (en) 2013-08-23 2018-06-08 Exxonmobil Upstream Research Company Simultaneous search of sources during seismic acquisition and seismic inversion
US10036818B2 (en) 2013-09-06 2018-07-31 Exxonmobil Upstream Research Company Accelerating full wavefield inversion with nonstationary point-spread functions
US9817143B2 (en) * 2013-10-30 2017-11-14 Pgs Geophysical As Methods and systems for constraining multiples attenuation in seismic data
US9910189B2 (en) 2014-04-09 2018-03-06 Exxonmobil Upstream Research Company Method for fast line search in frequency domain FWI
EP3140675A1 (en) 2014-05-09 2017-03-15 Exxonmobil Upstream Research Company Efficient line search methods for multi-parameter full wavefield inversion
US10185046B2 (en) 2014-06-09 2019-01-22 Exxonmobil Upstream Research Company Method for temporal dispersion correction for seismic simulation, RTM and FWI
MX362753B (en) 2014-06-17 2019-02-07 Exxonmobil Upstream Res Co Fast viscoacoustic and viscoelastic full-wavefield inversion.
US10838092B2 (en) 2014-07-24 2020-11-17 Exxonmobil Upstream Research Company Estimating multiple subsurface parameters by cascaded inversion of wavefield components
US10422899B2 (en) 2014-07-30 2019-09-24 Exxonmobil Upstream Research Company Harmonic encoding for FWI
US10386511B2 (en) 2014-10-03 2019-08-20 Exxonmobil Upstream Research Company Seismic survey design using full wavefield inversion
EP3210050A1 (en) 2014-10-20 2017-08-30 Exxonmobil Upstream Research Company Velocity tomography using property scans
AU2015363241A1 (en) 2014-12-18 2017-06-29 Exxonmobil Upstream Research Company Scalable scheduling of parallel iterative seismic jobs
US10520618B2 (en) 2015-02-04 2019-12-31 ExxohnMobil Upstream Research Company Poynting vector minimal reflection boundary conditions
WO2016130208A1 (en) 2015-02-13 2016-08-18 Exxonmobil Upstream Research Company Efficient and stable absorbing boundary condition in finite-difference calculations
US10670750B2 (en) 2015-02-17 2020-06-02 Exxonmobil Upstream Research Company Multistage full wavefield inversion process that generates a multiple free data set
US10416327B2 (en) 2015-06-04 2019-09-17 Exxonmobil Upstream Research Company Method for generating multiple free seismic images
US10838093B2 (en) 2015-07-02 2020-11-17 Exxonmobil Upstream Research Company Krylov-space-based quasi-newton preconditioner for full-wavefield inversion
MX2018002797A (en) 2015-10-02 2018-04-13 Exxonmobil Upstream Res Co Q-compensated full wavefield inversion.
CN108139498B (en) 2015-10-15 2019-12-03 埃克森美孚上游研究公司 FWI model domain angular stack with amplitude preservation
US10768324B2 (en) 2016-05-19 2020-09-08 Exxonmobil Upstream Research Company Method to predict pore pressure and seal integrity using full wavefield inversion
US20220137248A1 (en) * 2020-11-02 2022-05-05 China Petroleum & Chemical Corporation Computing program product and method for prospecting and eliminating surface-related multiples in the beam domain with deghost operator

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20070032954A1 (en) * 2003-05-02 2007-02-08 Westerngeco, L.L.C. Method for computing and using timing errors that occur in multiples predicted by multiple prediction algorithm
GB2429528A (en) * 2005-08-26 2007-02-28 Westerngeco Seismic Holdings Handling of static corrections in multiple prediction

Family Cites Families (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4887243A (en) * 1982-12-27 1989-12-12 Mobil Oil Corporation Removal of surface multiples
US5274605A (en) * 1992-06-26 1993-12-28 Chevron Research And Technology Company Depth migration method using Gaussian beams
US5986973A (en) * 1997-09-04 1999-11-16 Western Atlas International, Inc. Energy minimization in surface multiple attenuation
GB2348003B (en) * 1999-03-19 2001-02-07 Geco Prakla Seismic data processing method for data acquired using overlapping vibratory sweeps
US6169959B1 (en) * 1999-05-27 2001-01-02 Baker Hughes Incorporated Method of predicting kinematics for surface multiples
US6507787B1 (en) 1999-10-18 2003-01-14 Petroleo Brasileiro S.A.-Petrobras Method for the suppression of multiple reflections from marine seismic data
GB2368017B (en) * 2000-06-20 2004-05-12 Bournemouth University Higher Apparatus for electrical stimulation of the leg
JP2002357905A (en) * 2001-03-28 2002-12-13 Sumitomo Chem Co Ltd Resist composition
US6735527B1 (en) * 2003-02-26 2004-05-11 Landmark Graphics Corporation 3-D prestack/poststack multiple prediction
US6832161B1 (en) * 2003-08-15 2004-12-14 Westerngeco, L.L.C. Method for attenuating water layer multiples
US7149630B2 (en) * 2005-01-13 2006-12-12 Bp Corporation North America Inc. Method of DMO calculation for use in seismic exploration

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20070032954A1 (en) * 2003-05-02 2007-02-08 Westerngeco, L.L.C. Method for computing and using timing errors that occur in multiples predicted by multiple prediction algorithm
GB2429528A (en) * 2005-08-26 2007-02-28 Westerngeco Seismic Holdings Handling of static corrections in multiple prediction

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
BERKHOUT A J ET AL: "IMAGING OF MULTIPLE REFLECTIONS" 1 July 2006 (2006-07-01), GEOPHYSICS, SOCIETY OF EXPLORATION GEOPHYSICISTS, TULSA, OK, US, PAGE(S) SI209 - SI220 , XP001245339 ISSN: 0016-8033 the whole document *
GUITTON ET AL: "Multiple attenuation in complex geology with a pattern-based approach" GEOPHYSICS,, vol. 70, no. 4, 7 July 2005 (2005-07-07), pages V97-V107, XP002503073 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102478664A (en) * 2010-11-23 2012-05-30 中国石油天然气集团公司 Spatial sampling interval determining method without polluting effective signals

Also Published As

Publication number Publication date
CN101688922A (en) 2010-03-31
AU2008257075A1 (en) 2008-12-04
CN101688922B (en) 2012-05-30
WO2008147657A3 (en) 2009-03-26
EP1998189A3 (en) 2009-02-18
EP1998189A2 (en) 2008-12-03
AU2008257075B2 (en) 2013-08-08
US20080294346A1 (en) 2008-11-27
BRPI0811841A2 (en) 2014-11-18
EA017379B1 (en) 2012-12-28
US7715986B2 (en) 2010-05-11
EA200971083A1 (en) 2010-06-30
CA2687095C (en) 2016-01-05
CA2687095A1 (en) 2008-12-04
NO20093563L (en) 2010-02-19

Similar Documents

Publication Publication Date Title
US7715986B2 (en) Method for identifying and removing multiples for imaging with beams
US8326543B2 (en) System and method for attenuating multiples in seismic data that accounts for beam azimuth
US10371844B2 (en) Processing seismic data acquired using moving non-impulsive sources
US7477992B2 (en) Method for combining seismic data sets
US7953556B2 (en) Geophone noise attenuation and wavefield separation using a multi-dimensional decomposition technique
EP2787370B1 (en) Device and method for de-blending simultaneous shot data
US9817143B2 (en) Methods and systems for constraining multiples attenuation in seismic data
EP2626727A2 (en) Method and System for Determining Source Signatures After Source Ghost Removal
US20030048696A1 (en) Trapped water bottom multiple and peg-leg multiple suppression for ocean bottom seismic data
EP2730949A2 (en) Interference noise attenuation method and apparatus
MX2010005019A (en) Method for calculation of seismic attributes from seismic signals.
EP3025169B1 (en) Predicting interbed multiples in seismic data using beam decomposition
EP3446156B1 (en) Estimating an earth seismic response
EP1532472A2 (en) Removal of noise from seismic data using improved radon transformations
Tylor-Jones et al. Processing Essentials
Edme et al. Near-surface imaging from point-receiver ambient noise
Edme et al. Ambient noise for near-surface imaging
AU2010202339A1 (en) Method for dynamic aperture determination for three-dimensional surface-related multiple elimination

Legal Events

Date Code Title Description
WWE Wipo information: entry into national phase

Ref document number: 200880021551.9

Country of ref document: CN

121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 08747773

Country of ref document: EP

Kind code of ref document: A2

WWE Wipo information: entry into national phase

Ref document number: 2687095

Country of ref document: CA

NENP Non-entry into the national phase

Ref country code: DE

WWE Wipo information: entry into national phase

Ref document number: 2008257075

Country of ref document: AU

WWE Wipo information: entry into national phase

Ref document number: 2261/MUMNP/2009

Country of ref document: IN

ENP Entry into the national phase

Ref document number: 2008257075

Country of ref document: AU

Date of ref document: 20080507

Kind code of ref document: A

WWE Wipo information: entry into national phase

Ref document number: 200971083

Country of ref document: EA

122 Ep: pct application non-entry in european phase

Ref document number: 08747773

Country of ref document: EP

Kind code of ref document: A2

ENP Entry into the national phase

Ref document number: PI0811841

Country of ref document: BR

Kind code of ref document: A2

Effective date: 20091119