US8542554B2 - Method for velocity model building and analysis - Google Patents

Method for velocity model building and analysis Download PDF

Info

Publication number
US8542554B2
US8542554B2 US12/709,868 US70986810A US8542554B2 US 8542554 B2 US8542554 B2 US 8542554B2 US 70986810 A US70986810 A US 70986810A US 8542554 B2 US8542554 B2 US 8542554B2
Authority
US
United States
Prior art keywords
wavefield
coordinates
velocity model
reflection operator
error
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Expired - Fee Related, expires
Application number
US12/709,868
Other languages
English (en)
Other versions
US20110110190A1 (en
Inventor
Colin J. Thomson
Philip Kitchenside
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Westerngeco LLC
Original Assignee
Westerngeco LLC
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 Westerngeco LLC filed Critical Westerngeco LLC
Priority to US12/709,868 priority Critical patent/US8542554B2/en
Assigned to WESTERNGECO L. L. C. reassignment WESTERNGECO L. L. C. ASSIGNMENT OF ASSIGNORS INTEREST (SEE DOCUMENT FOR DETAILS). Assignors: THOMSON, COLIN J., KITCHENSIDE, PHILIP
Priority to PCT/US2010/056066 priority patent/WO2011059984A2/fr
Priority to EP10830600.2A priority patent/EP2499517A4/fr
Publication of US20110110190A1 publication Critical patent/US20110110190A1/en
Application granted granted Critical
Publication of US8542554B2 publication Critical patent/US8542554B2/en
Expired - Fee Related legal-status Critical Current
Adjusted expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/30—Analysis
    • G01V1/303—Analysis for determining velocity profiles or travel times
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00—Details of seismic processing or analysis
    • G01V2210/60—Analysis
    • G01V2210/62—Physical property of subsurface
    • G01V2210/622—Velocity, density or impedance
    • G01V2210/6222—Velocity; travel time

Definitions

  • This invention relates data processing for geophysical exploration, more specifically related to building velocity models for interpreting and refining a subsurface model in geological exploration.
  • Wave-equation migration uses a “one-way” or “two-way” wave equation to propagate the recorded data (sometimes referred to as the “receiver” wavefield) back into the subsurface. Simultaneously, the propagation of the wave field due to the seismic source is modeled throughout the subsurface. This wavefield is commonly referred to as the source wavefield.
  • Migration with the “one-way” equation is commonly referred to as wavefield-extrapolation migration (WEM), whereas migration with the “two-way” equation is referred to as reverse-time migration (RTM).
  • WEM wavefield-extrapolation migration
  • RTM reverse-time migration
  • CIGs common-image gathers
  • CIG common-image gathers
  • one type of CIG could be constructed from the energy at zero-time on the reflection operators, if the operators are computed for a sampled range of depths.
  • Such a CIG would then be a panel of migrated energy displayed as a function of the depth z and x, where x is the lateral separation of the down-going and up-going wavefield components. The latter are conventionally sampled at constant midpoint.
  • This invention pertains to seismic imaging where images of the subsurface geology are constructed using wave-equation migration.
  • the invention uses the propagated source and receiver wavefields modeled by reverse-time migration or wavefield-extrapolation migration to construct “wavefield-coordinates” reflection operators at locations throughout the subsurface.
  • the synthesized reflection operators are used to interpret and refine the subsurface model. This may be done locally by matching parametric approximations of reflection operators obtained from wave theory with the synthesized reflection operators obtained from the data.
  • the same concepts may also be used in a more global tomographic inversion approach.
  • FIG. 1 depicts a basic reflection-operator geometry.
  • FIG. 2 a depicts the asymptotic (stationary phase and branch point) behavior of a wavefield-coordinates reflection operator above, at and below a leftward-dipping plane interface, showing the kinematics of various signals in the (x′,t) domain when the migration velocity equals the true velocity above the reflector.
  • Solid curves indicate stationary-phase signals and dashed lines indicate branch-point signals.
  • FIG. 2 b depicts the asymptotic (stationary phase and branch point) behavior of a wavefield-coordinates reflection operator above, at and below a leftward-dipping plane interface, showing the kinematics of various signals in the (x′,t) domain when the migration velocity is somewhat less than the true velocity above the reflector.
  • FIG. 3 depicts six reflection operators computed during wave-equation migration of synthetic data at various x-locations within the water layer.
  • FIG. 4 depicts the amplitude behavior of a reflection operator located above a flat reflector associated with a velocity contrast.
  • the plot on the left of the Figure shows the reflection operator amplitude as a function of x and t and is for the case where the velocity above the reflector is greater than the velocity below.
  • the plot on the right is for when the velocity above the reflector is smaller.
  • FIG. 5 depicts six reflection operators computed during wave-equation migration of the synthetic data, at a depth below the water bottom.
  • FIG. 6 depicts a fitting of two wave-equation migration reflection operators computed above (left panel) and below (right panel) a planar dipping reflector.
  • FIG. 7 depicts a different fitting of two wave-equation migration reflection operators computed above (left panel) and below (right panel) a planar dipping reflector.
  • FIG. 8 depicts a third fitting of two wave-equation migration reflection operators computed above (left panel) and below (right panel) a planar dipping reflector.
  • FIG. 9 depicts six reflection operators computed at a depth close to structural complexity.
  • FIG. 10 depicts a computer system implementing a method.
  • Wave-equation migration also referred to as wavefield-continuation migration or WEM migration.
  • Wave-equation migration uses a “one-way” or “two-way” wave equation to propagate the recorded data (sometimes referred to as the “receiver” wavefield) back into the subsurface.
  • the propagation of the wave field due to the seismic source is modeled throughout the subsurface. This wavefield is commonly referred to as the source wavefield.
  • Migration with the “one-way” equation is commonly referred to as wavefield-extrapolation migration (WEM), whereas migration with the “two-way” equation is referred to as reverse-time migration (RTM).
  • WEM wavefield-extrapolation migration
  • RTM reverse-time migration
  • the propagated source and receiver wavefields may be combined in a way that provides an image of the subsurface geology.
  • the propagation of the wavefields requires an accurate estimate of the velocity model of the subsurface to be specified.
  • the methods disclosed here use the propagated source and receiver wavefields to construct so-called “wavefield-coordinates” reflection operators at locations throughout the subsurface. A qualitative and quantitative analysis of the reflection operators is then used to improve the model.
  • a particular reflection operator pertains to an (x,y,z) location within the subsurface and comprises a number of modeled seismic traces.
  • the z coordinate usually signifies the depth at which the reflection operator is constructed and x and y denote its lateral position. More generally, the coordinate frame could be tilted so that the surface of constant z is a plane that is not horizontal or flat-lying. Such a plane might arise in a targeted imaging method.
  • the important point and a key assumption are that a set of source and receiver wavefields are known across a sufficiently wide aperture on some plane of constant z (“depth”). In preferred embodiments of the ideas, there will be many such planes throughout the subsurface, but it suffices to consider one example in order to understand basic principles.
  • Each trace of the reflection operator represents the wavefield which would be recorded by a receiver at (x,y,z) due to a localized down-going source wavefield at another x′,y′ (at the same depth).
  • the reflection operator therefore characterizes the response of that part of the medium deeper than z.
  • FIG. 1 shows a flat datum 140 at which the reflection operator is defined or computed.
  • the reflection properties of the dipping interface can be characterized in the neighbourhood shown by the row of small hexagons 154 . Migration or further depth extrapolation with trial velocities is one way of doing this and the apparent reflection operator that is imaged in that way can be used to improve a preliminary velocity model.
  • a plane dipping reflector 130 which serves as the interface between two constant-velocity media 132 and 134 having velocities v 1 and v 2 respectively.
  • the reflection operator for this structure 130 will act upon an arbitrary downgoing or incident wave 110 at the datum level 140 so as to yield the upgoing or reflected wave 120 at the same datum level 140 .
  • the reflection operator is a function of x, x′ and t. If the incident wave 110 is localized at x′ (i.e. a spatial delta-function wavefield), then the response at x is a time series that looks like a quasi-point-source response of the dipping plane interface.
  • the reflection operator concept is independent of how the incident and reflected wavefields are obtained near a boundary, and so similar principles apply to WEM and RTM.
  • wavefields sampled over subsurface planes are used to estimate properties, rather than ray-based point data, as with Kirchhoff migration gathers.
  • the location x′ is used to label a trace within the reflection operator defined at a fixed x and x′-x is therefore the source—receiver subsurface wavefield offset.
  • the x′ are regularly sampled. This representation of the reflection response is referred to as a “wavefield-coordinates” representation because one of either x or x′ is held constant not their midpoint.
  • a subsurface extended image gather or CIG is in essence an estimate of the reflection operator at a point in some neighborhood of the reflector.
  • the CIG may also contain signals from many distant reflectors, but for simplicity, the analysis is done for the nearest or primary reflections.
  • wavefield-coordinates reflection operator or extended gather
  • the wavefield-coordinates reflection operator will appear kinematically similar to the wavefield that would be obtained if a source were to be placed just above depth z at lateral location x, such a wavefield being observed as a function of x′ and t.
  • a source at x′ and observe the operator or extended gather as a function of x and t.
  • wavefield coordinates (x,x′) the effects of dip on the reflection-operator gathers is revealed in a physically interpretable way before a final fully stacked image is formed.
  • a reflection operator is an Earth property and estimates of local reflection operators obtained by wavefield extrapolation/migration methods are potentially dynamically as well as kinematically interpretable, for example including total reflections and headwaves. They therefore provide a connection between amplitude variation with offset (AVO) analysis and imaging.
  • a key result leading to one embodiment is the wave-theory prediction of the dynamic form of the reflection operator (or gather) for a dipping plane interface. This will be discussed later. To put this embodiment in context, it is appropriate to briefly describe some algorithmic features of wave-equation migration.
  • the input data consist of shot records, where a shot record is a collection of seismic traces which are records of acoustic energy reflected from geological interfaces, the energy being due to a single seismic source.
  • WEM computations are usually carried out in the frequency-space ( ⁇ -x) domain, whereas RTM operates in the space-time domain.
  • RTM operates in the space-time domain.
  • an image is constructed by suitably combining two wavefields.
  • One wavefield is the extrapolated hypothetical source wavefield d which is typically an impulse function.
  • the other is the extrapolated recorded wavefield u.
  • Two well-known formulae for constructing a WEM image using frequency-space domain computations are:
  • the form chosen is unimportant except when the amplitude versus offset behaviour of a reflection operator is used to characterize the velocity contrast at the reflector. In this case, the second form should be used. Note that when computing R, the lateral spatial locations of d and u differ.
  • R(x,z,x′,t) is an estimate of the reflection operator that relates the down-going wavefield at location (x′,z) to the up-going wavefield at (x,z). Therefore it describes the reflection response of the Earth for locations deeper than z, if we ignore or remove possible contributions due to reflectors above.
  • R(x′,t) may be referred to as a reflection operator at x. It is a panel of energy which may be visualized as a function of x′ and time t, for each location x of interest at the current depth z.
  • x and x′ are wavefield coordinates, as opposed to the midpoint-offset coordinates that are generally used in imaging work. Significantly different results will be obtained when the wavefields d and u are sampled at constant midpoint, and in particular the effects of reflector dip may be obscured by using such coordinates.
  • the model update part of the method utilizes formulae derived from a wave-theoretical investigation of R(x,z,x′,t). Wave theory enables the reflection operators for simple earth models to be formulated as integrals over a slowness or ray parameter p.
  • the frequency domain form of R(x,z,x′,t) for a dipping reflector lying somewhat below the level z is given by
  • the above formula may be used to predict how the reflection operator gathers will evolve if the data are further extrapolated from the level z downwards towards the reflector and even to levels at which the extrapolated incident and reflected wavefields for each shot are taken beyond the physical limits of the layer above the interface (i.e. analytically continued).
  • FIG. 2 a shows the predicted kinematic behaviour of a reflection operator in the vicinity of a planar leftward-dipping reflector when the migration velocity equals the true local velocity above the interface.
  • the predictions were derived by asymptotic analysis of the full dynamical expression for the reflection operator in such a model.
  • the horizontal axis represents the spatial coordinate x′ and vertical axis represents the time t.
  • the form of the reflection operator is determined by many factors. In particular, the location x of the operator relative to the reflector, the acquisition geometry, the geometry of the reflector, the velocity used for the migration (the migration velocity) and the actual medium velocity. Note that when the reflection operator is computed at a depth greater than that of the reflector, or equivalently at a point x>0 to the right of the reflector, rightmost panel, then the asymmetric hyperbolic curve is inverted.
  • the method combines the described theoretical properties of reflection operators with numerical reflection operators computed during wave-equation migration.
  • FIG. 2 a There are signals in FIG. 2 a that correspond to physical arrivals, such as the dipping-reflector shifted hyperbola (solid) and tilted linear head-wave signals (dashed) in the left frame, which corresponds to points (x,x′) that lie above and to the left of the dipping interface. Other signals are not physical, such as the inverted hyperbola (solid) in the right frame that corresponds to data that have been extrapolated below or to the right of the reflector (over-migrated).
  • FIG. 2 b shows the corresponding predicted kinematic behavior of a reflection operator imaged in the vicinity of the planar leftward-dipping reflector when the migration velocity is somewhat less than the true local velocity above the interface.
  • the predictions were again derived by asymptotic analysis of the full dynamical expression for the wavefield-coordinates reflection operator in such a situation.
  • a primary claim in this method is that geometrical observations such as this can be interpreted in terms of the reflector properties (distance of the operator/gather location x from the interface, true medium velocity, reflector dip, etc.).
  • FIG. 3 shows six reflection operators which have been constructed during the migration of a synthetic dataset.
  • the six reflection operator panels have been displayed side by side as a single panel of data.
  • the operators have all been computed at the same depth, but at different lateral (x) locations within the water layer.
  • the dominant seismic event on each panel corresponds to the water bottom reflection.
  • the hyperbolic nature of this event is characteristic of a reflection operator located above a reflector, where the velocity model above the reflector is mildly laterally varying.
  • the amplitudes of refection operator events also provide information about the velocity contrast at a reflector.
  • FIG. 4 shows how, for the simple example of a flat boundary, the amplitude of a reflection operator varies with source-receiver offset as a function of the velocity contrast. Note that the amplitude of the reflection operator decreases with increasing source-receiver offset when the velocity above the reflector is the largest. If velocity increases, the converse is true for small source-receiver onsets corresponding to near normal incidence partial reflections.
  • FIG. 5 shows six reflection operators computed during wave-equation migration of the synthetic data at a depth below the water bottom.
  • upturned events are evident. From the theory (see FIG. 2 a or 2 b ), these arise when the reflection operator is computed at a depth greater than the true reflector depth.
  • events which have been focussed to an impulse with linear arms are also shown in FIG. 5 . These occur when the depth of the reflection operator coincides with the depth of the corresponding reflector (see centre panel of FIG. 2 a ). The slopes of the linear arms are determined by the velocities above and below the reflector in question.
  • the reflection operator leading-order kinematics can be parameterized in terms of the structure and the velocity model. Furthermore, the effects of velocity error can be obtained by perturbation analysis. For example, assume that the data at the primary datum level z referred to above are error free and consider a further extrapolation of the wavefields towards the reflector using an erroneous migration velocity v m . This erroneous local one-way extrapolation of the wavefields leads to new reflection operators or gathers at the new depth levels that will differ from what would have been obtained using the true velocity above the reflector. More generally the data at the primary datum level z will also have uncertainties, due to uncertainties in an overburden model. In that case, variations in the overburden propagator must be included, but for present purposes it suffices to assume that the primary datum values are accurate and to focus on locally induced distortions of the reflection operator.
  • kinematic approximation may be obtained by asymptotic (stationary phase) analysis of the analytically known reflection operator (for a dipping plane reflector between homogeneous media), and the effects of an erroneous migration velocity can be included.
  • asymptotic stationary phase
  • reflection operator for a dipping plane reflector between homogeneous media
  • the velocity error may be estimated by scanning over ⁇ v and observing the value giving the best agreement between the kinematic approximation and the wave-equation migration reflection operator.
  • a threshold of velocity errors is given as a measure of the quality the corresponding velocity model.
  • the velocity error determination may be done either quantitatively by performing a least-squares fit of x-t values picked from an event on the reflection operator or qualitatively by visual inspection.
  • FIGS. 6 to 8 illustrate examples using the qualitative approach for data from a model consisting of a single planar dipping reflector.
  • the WEM operators have been superimposed on the two aforementioned analytic approximations.
  • the WEM operators were computed using a migration velocity which was 10% smaller than the medium velocity.
  • the figures show results of varying the velocity error parameter in the analytic approximation.
  • FIG. 6 shows the fitting of two WEM reflection operators computed above (left panel) and below (right panel) a planar dipping reflector. It shows the WEM reflection operators 602 and 604 , small-dip small-offset approximation 612 and 614 , and small-velocity error approximation 622 and 624 .
  • the panel on the right of FIG. 6 shows that the lateral extent of the WEM reflection operator is reduced by the synthetic-data aperture limitations. In other words, there are insufficient shot records to the left and to the right sides of the x-location of the reflection operator.
  • FIG. 7 shows the fitting of two WEM reflection operators computed above (left panel) and below (right panel) a planar dipping reflector. It shows the WEM reflection operators 702 and 704 , small-dip small-offset approximation 712 and 714 , and small-velocity error approximation 722 and 724 .
  • FIG. 7 shows the result of using the correct value of the velocity error parameter. Note that the (stationary-phase) approximation breaks down in certain regions of the x-t plane. However, within the region of applicability there is good agreement.
  • FIG. 8 shows the fitting of two WEM reflection operators computed above (left panel) and below (right panel) a planar dipping reflector. It shows the WEM reflection operators 802 and 804 , small-dip small-offset approximation 812 and 814 , and small-velocity error approximation 822 and 824 .
  • FIG. 8 shows the result of using a velocity which is 2.5% too high.
  • FIG. 5 suggests that events on a reflection operator may be matched with the corresponding geological layers above and below the depth at which the reflection operator is calculated. This property together with the assumptions used for deriving the kinematic approximation suggests performing the fitting on a layer by layer basis. This type of fitting procedure is therefore integrated within a WEM based “layer-stripping” method for refining the velocity model.
  • the procedure for each layer is as follows:
  • the fitting procedure (Step 4) provides estimates of the velocity error, dip angle of the reflector (i.e. essentially the dip of the deeper lower boundary) and the vertical distance to the reflector.
  • dip is obtained from the local reflection-operator gather, rather than from a final fully-stacked image, and that this dip parameter is made more apparent by the use of wavefield coordinates.
  • FIG. 9 shows reflection operators computed at a depth proximate to structural complexity.
  • the operators are still interpretable in terms of the structure, but they are too complex to permit fitting with the kinematic approximation used above.
  • the model parameters must be derived using tomographic inversion.
  • the data comprising wavefield-coordinates reflection operators may be interpreted as being due to a seismic quasi-point-source within the subsurface (or, better, a dipole source). Wavefield-coordinates reflection operators are therefore conceptually closer to an actual subsurface or local seismic experiment. Also, dipping structure causes wavefield-coordinates reflection operators to be spatially asymmetric and hence the presence of dip is more easily interpretable. Lastly, the wavefield nature of the reflection operators permits both asymptotic analysis (e.g. involving Snell's Law) and the derivation of dynamic formulae for the purpose of refining the structural model. The latter pertains to AVO, total reflection and head-wave effects.
  • the wavefield-coordinates reflection operator may show the effects of AVO, including in principle total reflection and head-wave signals. This opens up the incident angle range from 0 degree (normal to the interface) to the angle of total reflection. With that information, the above method will lead to estimates of velocity contrasts which can be used for validation or further refinement of the model.
  • the methods described above are typically implemented in a computer system 1900 , one of which is shown in FIG. 10 .
  • the system computer 1930 may be in communication with disk storage devices 1929 , 1931 , 1933 and 1935 , which may be external hard disk storage devices. It is contemplated that disk storage devices 1929 , 1931 , 1933 and 1935 are conventional hard disk drives, and as such, will be implemented by way of a local area network or by remote access. Of course, while disk storage devices are illustrated as separate devices, a single disk storage device may be used to store any and all of the program instructions, measurement data, and results as desired.
  • seismic data from the seismic receivers may be stored in disk storage device 1931 .
  • Various non-seismic data from different sources may be stored in disk storage device 1933 .
  • the system computer 1930 may retrieve the appropriate data from the disk storage devices 1931 or 1933 to process data according to program instructions that correspond to implementations of various techniques described herein.
  • the program instructions may be written in a computer programming language, such as C++, Java and the like.
  • the program instructions may be stored in a computer-readable medium, such as program disk storage device 1935 .
  • Such computer-readable media may include computer storage media.
  • Computer storage media may include volatile and non-volatile, and removable and non-removable media implemented in any method or technology for storage of information, such as computer-readable instructions, data structures, program modules or other data.
  • Computer storage media may further include RAM, ROM, erasable programmable read-only memory (EPROM), electrically erasable programmable read-only memory (EEPROM), flash memory or other solid state memory technology, CD-ROM, digital versatile disks (DVD), or other optical storage, magnetic cassettes, magnetic tape, magnetic disk storage or other magnetic storage devices, or any other medium which can be used to store the desired information and which can be accessed by the system computer 1930 . Combinations of any of the above may also be included within the scope of computer readable media.
  • system computer 1930 may present output primarily onto graphics display 1927 , or alternatively via printer 1928 (not shown).
  • the system computer 1930 may store the results of the methods described above on disk storage 1929 , for later use and further analysis.
  • the keyboard 1926 and the pointing device (e.g., a mouse, trackball, or the like) 1925 may be provided with the system computer 1930 to enable interactive operation.
  • the system computer 1930 may be located at a data center remote from an exploration field.
  • the system computer 1930 may be in communication with equipment on site to receive data of various measurements.
  • the system computer 1930 may also be located on site in a field to provide faster feedback and guidance for the field operation.
  • Such data after conventional formatting and other initial processing, may be stored by the system computer 1930 as digital data in the disk storage 1931 or 1933 for subsequent retrieval and processing in the manner described above.
  • FIG. 19 illustrates the disk storage, e.g. 1931 as directly connected to the system computer 1930 , it is also contemplated that the disk storage device may be accessible through a local area network or by remote access.
  • disk storage devices 1929 , 1931 are illustrated as separate devices for storing input seismic data and analysis results, the disk storage devices 1929 , 1931 may be implemented within a single disk drive (either together with or separately from program disk storage device 1933 ), or in any other conventional manner as will be fully understood by one of skill in the art having reference to this specification.

Landscapes

  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Acoustics & Sound (AREA)
  • Environmental & Geological Engineering (AREA)
  • Geology (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics (AREA)
  • Geophysics And Detection Of Objects (AREA)
  • Radar Systems Or Details Thereof (AREA)
  • Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
US12/709,868 2009-11-10 2010-02-22 Method for velocity model building and analysis Expired - Fee Related US8542554B2 (en)

Priority Applications (3)

Application Number Priority Date Filing Date Title
US12/709,868 US8542554B2 (en) 2009-11-10 2010-02-22 Method for velocity model building and analysis
PCT/US2010/056066 WO2011059984A2 (fr) 2009-11-10 2010-11-09 Méthode de construction et d'analyse de modèles de vitesse
EP10830600.2A EP2499517A4 (fr) 2009-11-10 2010-11-09 Méthode de construction et d'analyse de modèles de vitesse

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
US25994709P 2009-11-10 2009-11-10
US12/709,868 US8542554B2 (en) 2009-11-10 2010-02-22 Method for velocity model building and analysis

Publications (2)

Publication Number Publication Date
US20110110190A1 US20110110190A1 (en) 2011-05-12
US8542554B2 true US8542554B2 (en) 2013-09-24

Family

ID=43974082

Family Applications (1)

Application Number Title Priority Date Filing Date
US12/709,868 Expired - Fee Related US8542554B2 (en) 2009-11-10 2010-02-22 Method for velocity model building and analysis

Country Status (3)

Country Link
US (1) US8542554B2 (fr)
EP (1) EP2499517A4 (fr)
WO (1) WO2011059984A2 (fr)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9470811B2 (en) 2014-11-12 2016-10-18 Chevron U.S.A. Inc. Creating a high resolution velocity model using seismic tomography and impedance inversion
US9513390B2 (en) 2014-11-03 2016-12-06 Westerngeco L.L.C. Compensating for space and slowness/angle blurring of reflectivity
US11320557B2 (en) 2020-03-30 2022-05-03 Saudi Arabian Oil Company Post-stack time domain image with broadened spectrum
US11340368B2 (en) 2019-03-21 2022-05-24 Saudi Arabian Oil Company Generating a velocity model and a density model of a subsurface structure of a reservoir

Families Citing this family (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB2443436B (en) * 2006-11-04 2010-07-21 Westerngeco Seismic Holdings Method of determing properties of the earth using turning wave data
SG185364A1 (en) * 2010-06-02 2012-12-28 Exxonmobil Upstream Res Co Efficient computation of wave equation migration angle gathers
EP2715403A4 (fr) * 2011-05-27 2015-01-28 Conocophillips Co Sélection de données ciblées d'une équation d'onde à deux sens destinée à l'imagerie améliorée de prospections parmi des structures géologiques complexes
US9025414B2 (en) 2011-05-27 2015-05-05 Conocophillips Company Reciprocal method two-way wave equation targeted data selection for seismic acquisition of complex geologic structures
CN102590859B (zh) * 2011-12-31 2014-01-22 中国石油集团西北地质研究所 垂向各向异性介质准p波方程逆时偏移方法
US9665604B2 (en) * 2012-07-31 2017-05-30 Schlumberger Technology Corporation Modeling and manipulation of seismic reference datum (SRD) in a collaborative petro-technical application environment
EP2755055B1 (fr) * 2013-01-11 2022-03-02 CGG Services SAS Tomographie par immersion pour modélisation de vitesse en profondeur en inversant des informations pré-empilées présentes dans des données sismiques migrées
US10379245B2 (en) * 2013-07-03 2019-08-13 Pgs Geophysical As Method and system for efficient extrapolation of a combined source-and-receiver wavefield
CN109188522B (zh) * 2018-10-09 2020-02-14 中国石油天然气股份有限公司 速度场构建方法及装置

Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5103429A (en) 1989-10-27 1992-04-07 Boris Gelchinsky Homeomorphical imaging method of analyzing the structure of a medium
US5260911A (en) * 1990-05-25 1993-11-09 Mason Iain M Seismic surveying
US5285422A (en) * 1992-10-23 1994-02-08 Western Atlas International, Inc. Method for compensating 3D DMO for the effects of an inhomogeneous earth
US5508973A (en) * 1994-06-06 1996-04-16 Western Atlas International, Inc. Method for determining the principal axes of azimuthal anisotropy from seismic P-wave data
US20030167126A1 (en) * 2002-01-15 2003-09-04 Westerngeco L.L.C. Layer stripping converted reflected waveforms for dipping fractures
US6839658B2 (en) 2000-07-17 2005-01-04 Sintef Petroleumsforskning Seismic processing with general non-hyperbolic travel-time corrections
US7039525B2 (en) 2003-09-23 2006-05-02 Pgs Geophysical As Method for seismic migration using explicit depth extrapolation operators with dynamically variable operator length
US20090257308A1 (en) 2008-04-11 2009-10-15 Dimitri Bevc Migration velocity analysis methods
US20100110830A1 (en) * 2006-11-04 2010-05-06 Colin Thomson Method of determining properties of the earth
US20100118652A1 (en) * 2008-11-13 2010-05-13 Jorg Friedrich Schneider Determination of depth moveout and of residual radii of curvature in the common angle domain

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6546339B2 (en) * 2000-08-07 2003-04-08 3D Geo Development, Inc. Velocity analysis using angle-domain common image gathers
WO2007145694A2 (fr) * 2006-06-12 2007-12-21 Exxonmobil Upstream Research Company Mise à jour de modèles de vitesse à l'aide de balayages de vitesse de migration
US7797110B2 (en) * 2007-06-26 2010-09-14 Shin's Geophysics Method for velocity analysis using waveform inversion in Laplace domain for geophysical imaging

Patent Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5103429A (en) 1989-10-27 1992-04-07 Boris Gelchinsky Homeomorphical imaging method of analyzing the structure of a medium
US5260911A (en) * 1990-05-25 1993-11-09 Mason Iain M Seismic surveying
US5285422A (en) * 1992-10-23 1994-02-08 Western Atlas International, Inc. Method for compensating 3D DMO for the effects of an inhomogeneous earth
US5508973A (en) * 1994-06-06 1996-04-16 Western Atlas International, Inc. Method for determining the principal axes of azimuthal anisotropy from seismic P-wave data
US6839658B2 (en) 2000-07-17 2005-01-04 Sintef Petroleumsforskning Seismic processing with general non-hyperbolic travel-time corrections
US20030167126A1 (en) * 2002-01-15 2003-09-04 Westerngeco L.L.C. Layer stripping converted reflected waveforms for dipping fractures
US6862531B2 (en) * 2002-01-15 2005-03-01 Westerngeco, L.L.C. Layer stripping converted reflected waveforms for dipping fractures
US7039525B2 (en) 2003-09-23 2006-05-02 Pgs Geophysical As Method for seismic migration using explicit depth extrapolation operators with dynamically variable operator length
US20100110830A1 (en) * 2006-11-04 2010-05-06 Colin Thomson Method of determining properties of the earth
US20090257308A1 (en) 2008-04-11 2009-10-15 Dimitri Bevc Migration velocity analysis methods
US20100118652A1 (en) * 2008-11-13 2010-05-13 Jorg Friedrich Schneider Determination of depth moveout and of residual radii of curvature in the common angle domain

Non-Patent Citations (8)

* Cited by examiner, † Cited by third party
Title
Biondi, et al, Angle-Domain Common-Image Gathers for Migration Velocity Analysis by Wavefield-Continuation Imaging, Geophysics, Sep.-Oct. 2004, pp. 1283-1298, vol. 69, No. 5.
International Search Report of PCT Application No. PCT/US2010/056066 dated Jun. 30, 2011: pp. 1-4.
N. Ross Hill, 'Prestack Gaussian-beam depth migration', Geophysics, vol. 66 No. 4 p. 1240-1250, Jul.-Aug. 2001. *
Rickett, et al., Offset and Angle-Domain Common-Angle Point Gathers for Shot-Profile Migration, Geophysics, May-Jun. 2002, pp. 883-889, vol. 67, No. 3.
Sava, Efficient Computation of Extended Images by Wavefield-Based Migration, SEG 79th Annual International Meeting, 2009.
Sava, et al., Angle-Domain Common-Image Gathers by Wavefield Continuation Methods, Geophysics, May-Jun. 2003, pp. 1065-1074, vol. 68, No. 3.
Sava, et al., Wave-Equation Migration Velocity Analysis I, Theory, Geophysical Prospecting, 2004, pp. 593-606, vol. 52.
Vasconcelos, et al., Wave-Equation Extended Images via Image-Domain Interferometry, SEG 79th Annual International Meeting, 2009.

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9513390B2 (en) 2014-11-03 2016-12-06 Westerngeco L.L.C. Compensating for space and slowness/angle blurring of reflectivity
US9470811B2 (en) 2014-11-12 2016-10-18 Chevron U.S.A. Inc. Creating a high resolution velocity model using seismic tomography and impedance inversion
US11340368B2 (en) 2019-03-21 2022-05-24 Saudi Arabian Oil Company Generating a velocity model and a density model of a subsurface structure of a reservoir
US11320557B2 (en) 2020-03-30 2022-05-03 Saudi Arabian Oil Company Post-stack time domain image with broadened spectrum

Also Published As

Publication number Publication date
EP2499517A4 (fr) 2014-07-30
EP2499517A2 (fr) 2012-09-19
WO2011059984A2 (fr) 2011-05-19
WO2011059984A3 (fr) 2011-10-27
US20110110190A1 (en) 2011-05-12

Similar Documents

Publication Publication Date Title
US8542554B2 (en) Method for velocity model building and analysis
US8902709B2 (en) Methods for concurrent generation of velocity models and depth images from seismic data
US6839658B2 (en) Seismic processing with general non-hyperbolic travel-time corrections
US8352190B2 (en) Method for analyzing multiple geophysical data sets
US5596548A (en) Seismic imaging using wave equation extrapolation
US7911878B2 (en) Sub-salt reflection tomography and imaging by walkaway VSP survey
Samec et al. Viscoelastic attenuation, anisotropy, and AVO
US5640368A (en) Migration velocity analysis using limited-aperture and monte carlo migration
US8395965B2 (en) Method of determining properties of the earth
US11614554B2 (en) Velocity model building for seismic data processing using PP-PS tomography with co-depthing constraint
US20140200813A1 (en) Systems and methods for seismic data processing using kinematic analysis of source-receive migration adcigs
US20160061975A1 (en) Vertical Seismic Profiling Migration Method
US10838092B2 (en) Estimating multiple subsurface parameters by cascaded inversion of wavefield components
US20160320505A1 (en) Device and method for mitigating cycle-skipping in full waveform inversion
EP3067718A1 (fr) Procédé et dispositif de tomographie de couche limite
US6868039B2 (en) Method for determining a seismic wave velocity model in a heterogeneous underground
US11635540B2 (en) Methods and devices performing adaptive quadratic Wasserstein full-waveform inversion
US8417458B2 (en) Removing ground-roll from geophysical data
Zhu et al. Amplitude and phase versus angle for elastic wide-angle reflections in the τ‐p domain
Li et al. Gaussian beam imaging of fractures near the wellbore using sonic logging tools after removing dispersive borehole waves
Gibson Jr et al. Modeling and velocity analysis with a wavefront-construction algorithm for anisotropic media
Bai Three dimensional seismic kinematic inversion with application to reconstruction of the velocity structure of Rabaul volcano
US20240184008A1 (en) System and method for multiple prediction with angular dependent reflectivity
US20230350089A1 (en) Full-waveform inversion with elastic mitigation using acoustic anisotropy
Wo Tomographic velocity model building in high-fidelity seismic depth imaging of near surface and VSP settings

Legal Events

Date Code Title Description
AS Assignment

Owner name: WESTERNGECO L. L. C., TEXAS

Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNORS:THOMSON, COLIN J.;KITCHENSIDE, PHILIP;SIGNING DATES FROM 20090814 TO 20090817;REEL/FRAME:023970/0333

STCF Information on status: patent grant

Free format text: PATENTED CASE

FPAY Fee payment

Year of fee payment: 4

FEPP Fee payment procedure

Free format text: MAINTENANCE FEE REMINDER MAILED (ORIGINAL EVENT CODE: REM.); ENTITY STATUS OF PATENT OWNER: LARGE ENTITY

LAPS Lapse for failure to pay maintenance fees

Free format text: PATENT EXPIRED FOR FAILURE TO PAY MAINTENANCE FEES (ORIGINAL EVENT CODE: EXP.); ENTITY STATUS OF PATENT OWNER: LARGE ENTITY

STCH Information on status: patent discontinuation

Free format text: PATENT EXPIRED DUE TO NONPAYMENT OF MAINTENANCE FEES UNDER 37 CFR 1.362

FP Lapsed due to failure to pay maintenance fee

Effective date: 20210924