US20100135115A1 - Multiple anisotropic parameter inversion for a tti earth model - Google Patents
Multiple anisotropic parameter inversion for a tti earth model Download PDFInfo
- Publication number
- US20100135115A1 US20100135115A1 US12/327,671 US32767108A US2010135115A1 US 20100135115 A1 US20100135115 A1 US 20100135115A1 US 32767108 A US32767108 A US 32767108A US 2010135115 A1 US2010135115 A1 US 2010135115A1
- Authority
- US
- United States
- Prior art keywords
- model
- tti
- data
- parameters
- values
- 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.)
- Abandoned
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. analysis, for interpretation, for correction
- G01V1/30—Analysis
- G01V1/303—Analysis for determining velocity profiles or travel times
Definitions
- the present invention relates to a method for determining values of anisotropic model parameters of a Tilted Transversely Isotropic (TTI) earth model.
- TTI Transversely Isotropic
- exploration geophysicists make many assumptions. One of them is that the subterranean formation is isotropic while in fact it is fundamentally anisotropic. This faulty assumption may result in erroneous imaging and interpretation of the geological structure. To extend the seismic processing techniques to anisotropic media, it is desirable to obtain a measure of the anisotropy of the geological structure.
- Seismic anisotropy can be defined as the dependence of seismic velocity on the direction of wave propagation. It is known that a transverse isotropy with tilted axis earth model or TTI earth model can be used to model the propagation of waves and obtain an image of the subterranean formation in anisotropic media.
- the physical parameters to describe a TTI earth model include (1) the symmetry axis, (2) P-wave (compressional) velocity along symmetry axis ⁇ Vp 0 , (3) a parameter that specifies how the velocities vary for small angles from the symmetry axis ⁇ , and (4) a parameter that determines the velocity at large angles from the axis of symmetry ⁇ (See Thomsen, “Weak Elastic Anisotropy”, Geophysics, vol. 51, no. 10, October 1986 and Alkhalifah and Tsvankin, “Velocity analysis for transversely isotropic Media”, Geophysics, vol. 60, 1550-1566, 1995).
- TTI earth models also use anisotropic parameter ⁇ to describe the propagation of waves in an anisotropic medium.
- S (shear) wave velocity is required to completely describe a TTI earth model, but in P-wave processing, S-wave velocity is usually obtained using an empirical relationship with P-wave velocity.
- TTI earth model is a three directional model. Each point in the model is described by its coordinates and the values of anisotropic parameters. In certain situations, only a few quantities of anisotropic parameters may be needed to fully define a model if the properties of the anisotropic medium do not change from point to point. However, in most situations, the TTI earth model requires a large number of spatially varying values of anisotropic parameters to accurately define the model.
- the anisotropic parameters of a TTI earth model may be directly measured from core data.
- drilling a well and coring are very expensive processes and direct measurements are only possible at very few well locations.
- 3D imaging it is desirable to determine the anisotropic parameters of the TTI earth model using also laterally extended data.
- the method including: acquiring input data for a geological volume of interest; determining a theoretical relationship between the input data and the anisotropic model parameters; and calculating the values of the anisotropic model parameters at each of a plurality of subsurface locations in the geological volume of interest based on the theoretical relationships and the input data using workflows involving iterative or sequential combinations of processes including input data preprocessing, conventional tomographic inversion, three dimensional tomographic inversion based on a tilted trans
- a computer implemented method for determining values of anisotropic parameters of a three dimensional TTI earth model, the anisotropic parameters including P-wave velocity (Vp 0 ) along a tilted symmetry axis, the Thomsen anisotropy parameters ⁇ , ⁇ (or ⁇ ( ⁇ )/(1+2 ⁇ )) representative of variation of wave velocities as a function of wave propagation angle from the symmetry axis, the method including (a) obtaining an initial migration velocity model that substantially flattens common-imaging-point gathers and substantially ties seismic data to well data, the initial migration velocity model including initial values Vp 0 , ⁇ , ⁇ (or ⁇ ) at each of a plurality of subsurface locations in a geological volume of interest; (b) inputting checkshot data, or VSP data, or both checkshot and VSP data at well locations into a three dimensional tomographic inversion to determine updated values of Vp 0 near the well locations, the values of Vp 0 being updated by a
- FIG. 1 shows a method for determining values of anisotropic model parameters of a Tilted Transversely Isotropic (TTI) Earth model in accordance with an embodiment of the invention
- FIG. 2 shows a workflow for determining values of anisotropic parameters in accordance with an embodiment of the invention.
- FIG. 3 shows various scenarios that can be used in the workflow of FIG. 2 , in accordance with an embodiment of the invention.
- Embodiments of the invention provide a method for accurately determining values of anisotropic parameters of a Titled Transversely Isotropic (TTI) earth model at each of a plurality of locations in a geological volume of interest.
- the estimated anisotropic parameters are calibrated with direct measurements at well locations.
- Embodiments of the invention also relate to a computer product having machine executable instructions, the instructions being executable by a machine to perform a method for determining values of anisotropic parameters of a TTI earth model.
- the method for accurately determining values of anisotropic parameters of a TTI earth model include three main aspects which are (1) the use of input data to determine the anisotropic parameters, (2) the relationship between the input data and the anisotropic parameters, and (3) the process employed to convert the input data to values of the anisotropic parameters.
- the process that is used to convert input data to anisotropic parameters may be referred to as an inversion process.
- FIG. 1 this figure shows a method 10 for determining values of anisotropic parameters of a TTI earth model in accordance with an embodiment of the invention.
- the values of the anisotropic parameters (Vp 0 , ⁇ , ⁇ and/or ⁇ ) are determined for each coordinate (x, y, z) in the geological volume of interest covered by the TTI earth model.
- the method starts at procedure 15 where input data are acquired for a geological volume of interest.
- Input data may include, for example, surface seismic data, vertical seismic profile data, VSP (vertical seismic profile) data, check shot data, well log data, interpretational data, regional trend, a priori data, or any combination of the foregoing.
- the interpretational data includes picks of horizon and shape.
- VSP and check shot data which record the direct arrival of waves from different directions, are also sensitive to anisotropic parameters. The measurement is travel-time at different borehole locations.
- VSP data are acquired by positioning the wave receiver down the borehole and the wave emitter near the surface.
- Check shot data are acquired by positioning both the wave emitter and receiver down the borehole.
- the method 10 proceeds to procedure 20 where a theoretical relationship is determined between the input data and the anisotropic parameters.
- the ray theory can be used to derive the theoretical relationship between the recorded travel-time (input data) and anisotropic parameters of the TTI earth model.
- the theoretical relationship between the input data and the anisotropic parameters can be obtained by first implementing the TTI specification techniques described in Thomsen, “ Weak elastic anisotropy” Geophysics, vol. 51, pgs. 1954-66 (1986); and/or Alkahalifah, et al., “ Velocity analysis for transversely isotropic media” Geophysics, vol. 60, pgs.
- anisotropic Earth model velocities 1550-1566 (1995) to specify anisotropic Earth model velocities.
- the anisotropic ray tracing is then performed by one or more techniques described Cerveny, “ Seismic rays and ray intensities in inhomogeneous anisotropic media,” Geophysical Journal, vol. 29, pgs. 1-13 (1972) and/or Gajewski et al., “ Vertical seismic profile synthetics by dynamic ray tracing in laterally varying layered anisotropic structures” Journal of Geophysics Res., vol. 95, pgs. 11301-11315, (1990). Additional information regarding the use of ray theory to derive the theoretical relationship between the input data and the anisotropic parameters can be gleaned from U.S.
- the method proceeds to procedure 25 where the values of the anisotropic model parameters at each of a plurality of subsurface locations in the geological volume of interest are calculated based on the theoretical relationships and the input data.
- the calculation of the anisotropic parameters uses workflows involving iterative or sequential combinations of processes including input data preprocessing, conventional tomographic inversion or three dimensional tomographic inversion based on a tilted transversely isotropic model, and three dimensional pre-stack depth migration using a tilted transversely isotropic model.
- the iterative or sequential combinations of processes may include interpretive picking.
- the estimation of the anisotropic parameters of procedure 25 is a computational process which uses the recorded data and a theoretical relationship to solve the model parameters. This process is called tomographic inversion.
- the inversion can be performed differently using various workflows.
- the anisotropic parameters may be estimated sequentially (i.e. one at a time) or, preferably, various parameters may be estimated simultaneously. It is also possible to use one type of data to estimate one or multiple parameters. In another embodiment, all available input data are used to estimate all anisotropic parameters simultaneously.
- the values of the anisotropic model parameters can be determined using various combinations of inversions with full sets or subsets of the input data and full sets or subsets of the output data.
- the tomographic inversion of procedure 25 is recursive and iterative, and may use three dimensional tilted transversely isotropic ray tracing to model three dimensional wave propagation.
- Input data preprocessing may include normalizing the input before performing the three dimensional TTI tomographic inversion.
- Input data preprocessing may also include, prior to performing the three dimensional TTI tomographic inversion, migrating the seismic data using 3D TTI prestack depth migration algorithms, sorting the migrated seismic data according to their subsurface location and their migration offset/angle into common-image-point (CIP) gathers, and selecting and quantifying residual moveouts in the common-image-point gathers domain.
- TTI prestack depth migration algorithms are known in the art.
- the values of the anisotropic parameters are calculated by iteratively selecting values of Vp 0 , ⁇ and ⁇ . Specifically, the values of the anisotropic parameters are calculated until (a) seismic image positions of subsurface structures in the geological volume of interest tie their spatial positions recognized in drilled wells and (b) residual moveouts in common-imaging-point gathers are minimum at every point in the model.
- the subsurface structures include rock boundaries.
- this figure shows a workflow 100 for determining values of anisotropic parameters of a Titled Transversely Isotropic (TTI) earth model in accordance with an embodiment of the invention.
- the iterative workflow can begin at procedure 105 where three dimensional seismic reflection data are measured.
- the workflow then proceeds to procedure 110 where the seismic reflection data are used to determine an initial TTI earth model for the geological volume of interest.
- the initial TTI earth model may be determined using various means.
- the initial TTI earth model can be determined from velocity information and various processing performed on measured data as well as knowledge, either regionally or globally, about the area in which the geological volume of interest is located.
- the initial TTI model building also includes estimating subsurface reflector structural dips and obtaining TTI symmetry axes in the geological volume of interest.
- a value of each parameter Vp 0 , ⁇ and ⁇ is provided at each of a plurality of locations (x, y, z) in the geological volume of interest.
- the workflow proceeds to block 115 and procedure 120 where well data (block 115 ) are used to perform a three dimensional (3D) TTI tomography (procedure 120 ) based on ray tracing to update the value of Vp 0 of the TTI earth model near the well(s).
- Well data may include VSP data or checkshot data, which represent information about wave velocities in the well(s).
- the updated values of Vp 0 in the geological volume of interest near the well(s) corresponds to Vp 0 (initial)+ ⁇ Vp 0 .
- the values of Vp 0 near the well(s) are updated in the initial TTI earth model.
- the VSP data and checkshot data are used in the embodiment of FIG. 2 to tie the seismic data to the well data.
- the transformations of procedures 110 - 125 enable one to refine the initial 3D TTI earth model by providing a more accurate TTI model near the well(s).
- the workflow 100 proceeds to procedure 130 where the remaining points in the model, i.e. the points outside the location(s) near the well(s), are populated.
- the remaining points in the model are populated using a three dimensional extrapolation process and the updated values of Vp 0 and ⁇ near the well(s).
- the result of the extrapolation process provides a three dimensional TTI earth model with well tie having updated values of Vp 0 and ⁇ at each point of the model.
- the three dimensional TTI earth model of block 135 includes the following parameters at each point (x, y, z) of the model: Vp 0 (initial)+ ⁇ Vp 0 , ⁇ (initial)+ ⁇ and ⁇ (initial).
- the model determined at block 135 is used to perform a three dimensional (3D) pre-stack TTI depth imaging of seismic reflection data at procedure 140 .
- seismic reflection data that have been recorded for the geological volume of interest are processed/migrated by the model to provide a new image of the sub-surface of the geological volume of interest.
- the migration of seismic data will typically enable an image to be formed of the geological volume of interest from the migrated seismic data that is a more accurate depiction of the geological features present in the seismic volume of interest than an image formed from migrated data with an isotropic algorithm.
- the three dimensional pre-stack TTI depth imaging of seismic reflection data of procedure 140 provides common image point gathers (block 145 ). As known in the art, common image point gathers correspond to the migrated seismic data that end up at a same image position in the geological volume of interest.
- residual moveout analysis is a basic step in velocity model refinement.
- the analysis is generally carried out using semblances for a grid of inline and xline locations.
- One semblance panel is computed for each common-image-point gather as a function of depth and offset.
- the residual moveout analysis on the common image point gathers is carried out for the near to mid offsets at procedure 150 while the residual moveout analysis on the common image point gathers is carried out for the near to far mid offsets at procedure 155 .
- Near and far offsets correspond, respectively, to a small and a large distance between the source and the receiver.
- the procedures 150 , 155 of FIG. 2 are intended to cover multiple moveout analyses. That is, the residual moveout analysis of FIG. 2 can be performed only for the near to mid offsets (procedure 150 ) or only for the near to far offsets (procedure 155 ) in embodiments of the invention. Alternatively, in an embodiment of the invention, and as shown in FIG. 2 , the residual moveout analysis can be performed for both the near to mid offsets and the near to far offsets simultaneously.
- the residual moveout analysis provides updated wave velocities Vp 0 for each point of the 3D TTI model defined at procedure 135 .
- a determination is made at blocks 165 a and/or 165 b as to whether convergence is obtained, i.e. whether the results of the residual moveout analyses of procedures 150 and/or procedure 155 is/are within a predetermined tolerance. If convergence is obtained, the workflow 100 ends at block 166 .
- the workflow 100 proceeds to block 170 and/or block 180 wherein it is determined whether the values of ⁇ and/or ⁇ should be updated for each point of the model. If the result of the inquiry is positive at block 170 , the workflow proceeds to procedure 185 where a three dimensional (3D) TTI tomography is performed to update the value of ⁇ for each point of the model using the results of the residual moveout analysis for the near to mid offsets. It is greatly beneficial to update the values of ⁇ only because such an update does not modify the values of wave velocities Vp 0 in the model. As a result, the tie between well and seismic data that was previously obtained with the model constructed at block 135 is not qualitatively destroyed in a low dip setting. In that way, it is possible to obtain a faster and more robust conversion for the values of the TTI model.
- 3D three dimensional
- ⁇ it is determined whether the values of ⁇ should be updated for each point of the model. If the result of the inquiry is positive at block 180 , the workflow proceeds to procedure 195 where a three dimensional (3D) TTI tomography is performed to update the value of ⁇ for each point of the model using the results of the residual moveout analysis for the far to mid offsets.
- 3D three dimensional
- the results of the 3D TTI tomography analyses of procedures 185 , 190 and 195 provide a new three dimensional (3D) TTI model at block 197 with updated values of ⁇ , Vp 0 and ⁇ at each point (x, y, z) of the model. It is then determined whether convergence is obtained for the updated values of ⁇ , Vp 0 and ⁇ in the new model of block 197 (block 198 ). Various tests may be used to determine whether convergence is obtained. For example, it is determined whether the obtained values are below a predetermined threshold. If the result of the inquiry is positive, the workflow of FIG. 2 ends at procedure 199 .
- the workflow 100 proceeds back to block 115 where the well data are used to perform a three dimensional (3D) TTI tomography (procedure 120 ) based on ray tracing to update the value of Vp 0 of the TTI earth model obtained at block 197 near the well(s).
- the TTI earth model obtained at block 197 is used as a new initial model for the subsequent iteration.
- this new initial TTI earth model is more refined than the initial model used at the first iteration in that it substantially flattens common-imaging-point gathers and substantially ties seismic data to well data.
- the workflow 100 After determining a modified model with updated values of Vp 0 near the well(s), the workflow 100 proceeds to procedures 125 - 198 where updated values of ⁇ and/or Vp 0 and/or ⁇ are determined. The workflow is then iterated until convergence is obtained for the values of ⁇ and/or Vp 0 and/or ⁇ (i.e. the values of ⁇ and/or Vp 0 and/or ⁇ substantially do not change between two subsequent iterations).
- FIG. 2 shows the various scenarios that can be applied at each iteration.
- ⁇ is the only parameter that is updated at each location in the model.
- the first scenario is preferred because the update of ⁇ does not change the values of Vp 0 in a low dip setting and, as a result, the quantitative tie between the well data and the seismic data is not destroyed during the optimization of ⁇ .
- parameters ⁇ and Vp 0 are updated.
- parameters ⁇ and ⁇ are updated.
- the velocities Vp 0 are updated.
- Vp 0 and ⁇ are updated.
- ⁇ is updated and, in the seventh scenario, ⁇ , ⁇ , and Vp 0 are updated.
- the one or more parameters updated at a given iteration may not be the same as the one or more parameters updated at a subsequent iteration.
- the anisotropic parameters ( ⁇ , ⁇ and Vp 0 ) optimized at each iteration may be different.
- various scenarios could be used to optimize the three dimensional TTI model.
- the workflow of FIG. 2 may be pursued to optimize a second parameter, e.g. ⁇ , and then a third parameter, e.g. Vp 0 . Therefore, the parameters ⁇ , ⁇ and Vp 0 can be optimized simultaneously or sequentially.
- the selection of the parameters updated at each iteration is highly dependent on the maturity of the overall model building process that includes both depth imaging and tomography and the type and the quality of well and seismic data that are used to construct the 3D TTI model.
- Tilted Transversely Isotropic (TTI) Earth model may be executed according to machine executable instructions or codes. These machine executable instructions may be embedded in a data storage medium. A processor may be configured to execute the instructions.
- TTI Transversely Isotropic
- Software functionalities of a computer system involving programming, including executable codes, may be used to implement the above described model.
- the software code may be executable by a general-purpose computer.
- the code and possibly the associated data records may be stored within a general-purpose computer platform.
- the software may be stored at other locations and/or transported for loading into an appropriate general-purpose computer-system.
- the embodiments discussed above involve one or more software or computer products in the form of one or more modules of code carried by at least one machine-readable medium. Execution of such codes by a processor of the computer system enables the platform to implement the functions in essentially the manner performed in the embodiments discussed and illustrated herein.
- Non-volatile media include, for example, optical or magnetic disks, such as any of the storage devices in any computers) operating as discussed above.
- Volatile media include dynamic memory, such as the main memory of a computer system.
- Physical transmission media include coaxial cables, copper wires and fiber optics, including the wires that comprise a bus within a computer system.
- Carrier-wave transmission media can take the form of electric or electromagnetic signals, or acoustic or light waves such as those generated during radio frequency (RF) and infrared (IR) data communications.
- Common forms of computer-readable media therefore include, for example: a floppy disk, a flexible disk, hard disk, magnetic tape, any other magnetic medium, a CD-ROM, DVD, any other optical medium, less commonly used media such as punch cards, paper tape, any other physical medium with patterns of holes, a RAM, a PROM, and EPROM, a FLASH-EPROM, any other memory chip or cartridge, a carrier wave transporting data or instructions, cables or links transporting such a carrier wave, or any other medium from which a computer can read or send programming codes and/or data.
- Many of these forms of computer readable media maybe involved in carrying one or more sequences of one or more instructions to a processor for execution.
Abstract
A method for determining values of anisotropic model parameters of a Tilted Transversely Isotropic (TTI) Earth model, the anisotropic parameters including P-wave velocity (Vp0) along a tilted symmetry axis, the Thomsen anisotropy parameters δ and ε (or η=(ε−δ)/(1+2δ)) representative of variations of wave velocities as a function of wave propagation angle from the symmetry axis, the method including acquiring input data for a geological volume of interest; determining a theoretical relationship between the input data and the anisotropic model parameters; and calculating the values of the anisotropic model parameters at each of a plurality of subsurface locations in the geological volume of interest based on the theoretical relationships and the input data using workflows involving iterative or sequential combinations of processes including input data preprocessing, conventional tomographic inversion, three dimensional tomographic inversion based on a tilted transversely isotropic model, and three dimensional pre-stack depth migration using a tilted transversely isotropic model.
Description
- The present invention relates to a method for determining values of anisotropic model parameters of a Tilted Transversely Isotropic (TTI) earth model.
- In order to analyze a geological structure of a subterranean formation, exploration geophysicists make many assumptions. One of them is that the subterranean formation is isotropic while in fact it is fundamentally anisotropic. This faulty assumption may result in erroneous imaging and interpretation of the geological structure. To extend the seismic processing techniques to anisotropic media, it is desirable to obtain a measure of the anisotropy of the geological structure.
- Seismic anisotropy can be defined as the dependence of seismic velocity on the direction of wave propagation. It is known that a transverse isotropy with tilted axis earth model or TTI earth model can be used to model the propagation of waves and obtain an image of the subterranean formation in anisotropic media. The physical parameters to describe a TTI earth model include (1) the symmetry axis, (2) P-wave (compressional) velocity along symmetry axis −Vp0, (3) a parameter that specifies how the velocities vary for small angles from the symmetry axis −δ, and (4) a parameter that determines the velocity at large angles from the axis of symmetry −η (See Thomsen, “Weak Elastic Anisotropy”, Geophysics, vol. 51, no. 10, October 1986 and Alkhalifah and Tsvankin, “Velocity analysis for transversely isotropic Media”, Geophysics, vol. 60, 1550-1566, 1995).
- Some TTI earth models also use anisotropic parameter ε to describe the propagation of waves in an anisotropic medium. Parameter ε satisfies the following relationship η=(ε−δ)/(1+2δ). S (shear) wave velocity is required to completely describe a TTI earth model, but in P-wave processing, S-wave velocity is usually obtained using an empirical relationship with P-wave velocity.
- Usually a TTI earth model is a three directional model. Each point in the model is described by its coordinates and the values of anisotropic parameters. In certain situations, only a few quantities of anisotropic parameters may be needed to fully define a model if the properties of the anisotropic medium do not change from point to point. However, in most situations, the TTI earth model requires a large number of spatially varying values of anisotropic parameters to accurately define the model.
- The anisotropic parameters of a TTI earth model may be directly measured from core data. However, drilling a well and coring are very expensive processes and direct measurements are only possible at very few well locations. For 3D imaging, it is desirable to determine the anisotropic parameters of the TTI earth model using also laterally extended data.
- In an aspect of the invention, there is provided a computer implemented method for determining values of anisotropic model parameters of a Tilted Transversely Isotropic (TTI) Earth model, the anisotropic parameters including P-wave velocity (Vp0) along a tilted symmetry axis, the Thomsen anisotropy parameters δ and ε (or η=(ε−δ)/(1+2δ)) representative of variations of wave velocities as a function of wave propagation angle from the symmetry axis, the method including: acquiring input data for a geological volume of interest; determining a theoretical relationship between the input data and the anisotropic model parameters; and calculating the values of the anisotropic model parameters at each of a plurality of subsurface locations in the geological volume of interest based on the theoretical relationships and the input data using workflows involving iterative or sequential combinations of processes including input data preprocessing, conventional tomographic inversion, three dimensional tomographic inversion based on a tilted transversely isotropic model, and three dimensional pre-stack depth migration using a tilted transversely isotropic model.
- In another aspect of the invention, there is provided a computer product having machine executable instructions, the instructions being executable by a machine to perform a tomographic inversion method for determining values of anisotropic parameters of a TTI earth model, the anisotropic parameters including P-wave velocity (Vp0) along a tilted symmetry axis, the Thomsen anisotropy parameters δ and ε (or η=(ε−δ)/(1+2δ)) representative of a variation of wave velocities as a function of wave propagation angle from the symmetry axis, the method including determining a relationship between input data and the anisotropic parameters, the input data being acquired for a geological volume of interest; and calculating the values of the anisotropic parameters at each of a plurality of subsurface locations in the geological volume of interest based on the relationship and the input data using workflows, involving iterative or sequential combinations of processes including input data preprocessing, three dimensional tomographic inversion, and three dimensional TTI pre-stack depth migration.
- A computer implemented method for determining values of anisotropic parameters of a three dimensional TTI earth model, the anisotropic parameters including P-wave velocity (Vp0) along a tilted symmetry axis, the Thomsen anisotropy parameters δ, ε (or η=(ε−δ)/(1+2δ)) representative of variation of wave velocities as a function of wave propagation angle from the symmetry axis, the method including (a) obtaining an initial migration velocity model that substantially flattens common-imaging-point gathers and substantially ties seismic data to well data, the initial migration velocity model including initial values Vp0, δ, ε (or η) at each of a plurality of subsurface locations in a geological volume of interest; (b) inputting checkshot data, or VSP data, or both checkshot and VSP data at well locations into a three dimensional tomographic inversion to determine updated values of Vp0 near the well locations, the values of Vp0 being updated by a correction ΔVp0, wherein Vp0=Vp0+ΔVp0; (c) determining an incremental improvement Δδ to δ(initial) using the relative change Δδ=(ΔVp0)/Vp0; (d) extrapolating the relative change Δδ from near-well locations to the entire three dimensional TTI earth model at each of the plurality of subsurface locations, taking into account geological consistency and regularization, to determine updated values of δ, wherein δ=δ+Δδ; (e) determining updated values of Vp0=Vp0 (1−Δδ) using the extrapolated three dimensional Δδ at each of the plurality of subsurface locations and obtaining the three dimensional extended incremental update ΔVp0=−Δδ Vp0; (f) inputting near-to-mid-offset/angle residual moveout information in the common-imaging-point gathers obtained with an improved migration velocity model that is defined with the updated values of Vp0 and δ into a TTI tomographic inversion process to further provide updated values of δ at each of the plurality of subsurface locations; and (g) inputting near-to-far-offset/angle residual moveout information in the common-imaging-point gathers obtained with the improved model the TTI tomographic inversion process to provide updated values of η at each of the plurality of subsurface locations.
- These and other objects, features, and characteristics of the present invention, as well as the methods of operation and functions of the related elements of the structure and the combination of parts and economies of manufacture, will become more apparent upon consideration of the following description and the appended claims with reference to the accompanying drawings, all of which form a part of this specification, wherein like reference numerals designate corresponding parts in the various figures. It is to be expressly understood, however, that the drawings are for the purpose of illustration and description only and are not intended as a definition of the limits of the invention. As used in the specification and in the claims, the singular form of “a”, “an”, and “the” include plural referents unless the context clearly dictates otherwise.
-
FIG. 1 shows a method for determining values of anisotropic model parameters of a Tilted Transversely Isotropic (TTI) Earth model in accordance with an embodiment of the invention; -
FIG. 2 shows a workflow for determining values of anisotropic parameters in accordance with an embodiment of the invention; and -
FIG. 3 shows various scenarios that can be used in the workflow ofFIG. 2 , in accordance with an embodiment of the invention. - Embodiments of the invention provide a method for accurately determining values of anisotropic parameters of a Titled Transversely Isotropic (TTI) earth model at each of a plurality of locations in a geological volume of interest. In an embodiment, the estimated anisotropic parameters are calibrated with direct measurements at well locations. Embodiments of the invention also relate to a computer product having machine executable instructions, the instructions being executable by a machine to perform a method for determining values of anisotropic parameters of a TTI earth model.
- In an embodiment, the method for accurately determining values of anisotropic parameters of a TTI earth model include three main aspects which are (1) the use of input data to determine the anisotropic parameters, (2) the relationship between the input data and the anisotropic parameters, and (3) the process employed to convert the input data to values of the anisotropic parameters. The process that is used to convert input data to anisotropic parameters may be referred to as an inversion process.
- Referring now to
FIG. 1 , this figure shows amethod 10 for determining values of anisotropic parameters of a TTI earth model in accordance with an embodiment of the invention. With the method ofFIG. 1 , the values of the anisotropic parameters (Vp0, δ, ηand/or ε) are determined for each coordinate (x, y, z) in the geological volume of interest covered by the TTI earth model. The method starts atprocedure 15 where input data are acquired for a geological volume of interest. Input data may include, for example, surface seismic data, vertical seismic profile data, VSP (vertical seismic profile) data, check shot data, well log data, interpretational data, regional trend, a priori data, or any combination of the foregoing. The interpretational data includes picks of horizon and shape. - Surface seismic data, which record reflections from the earth, have great coverage and are readily available in most areas. The wave reflected at different angles is sensitive to anisotropic parameters. Therefore, surface seismic data can be used as input data for the anisotropic parameter estimation. The measurable data include travel-times at different offsets and depth error of migrated gathers at different offsets/angles. VSP and check shot data, which record the direct arrival of waves from different directions, are also sensitive to anisotropic parameters. The measurement is travel-time at different borehole locations. As known in the art, VSP data are acquired by positioning the wave receiver down the borehole and the wave emitter near the surface. Check shot data are acquired by positioning both the wave emitter and receiver down the borehole.
- After acquiring input data, the
method 10 proceeds toprocedure 20 where a theoretical relationship is determined between the input data and the anisotropic parameters. The ray theory can be used to derive the theoretical relationship between the recorded travel-time (input data) and anisotropic parameters of the TTI earth model. Specifically, the theoretical relationship between the input data and the anisotropic parameters can be obtained by first implementing the TTI specification techniques described in Thomsen, “Weak elastic anisotropy” Geophysics, vol. 51, pgs. 1954-66 (1986); and/or Alkahalifah, et al., “Velocity analysis for transversely isotropic media” Geophysics, vol. 60, pgs. 1550-1566 (1995) to specify anisotropic Earth model velocities. The anisotropic ray tracing is then performed by one or more techniques described Cerveny, “Seismic rays and ray intensities in inhomogeneous anisotropic media,” Geophysical Journal, vol. 29, pgs. 1-13 (1972) and/or Gajewski et al., “Vertical seismic profile synthetics by dynamic ray tracing in laterally varying layered anisotropic structures” Journal of Geophysics Res., vol. 95, pgs. 11301-11315, (1990). Additional information regarding the use of ray theory to derive the theoretical relationship between the input data and the anisotropic parameters can be gleaned from U.S. patent application Ser. No. 12/079,170, System and Method for Migrating Seismic Data, filed on Mar. 24, 2008. - After determining the theoretical relationship between the input data and the anisotropic parameters, the method proceeds to
procedure 25 where the values of the anisotropic model parameters at each of a plurality of subsurface locations in the geological volume of interest are calculated based on the theoretical relationships and the input data. In an embodiment, the calculation of the anisotropic parameters uses workflows involving iterative or sequential combinations of processes including input data preprocessing, conventional tomographic inversion or three dimensional tomographic inversion based on a tilted transversely isotropic model, and three dimensional pre-stack depth migration using a tilted transversely isotropic model. Further, the iterative or sequential combinations of processes may include interpretive picking. - The estimation of the anisotropic parameters of
procedure 25 is a computational process which uses the recorded data and a theoretical relationship to solve the model parameters. This process is called tomographic inversion. The inversion can be performed differently using various workflows. For example, the anisotropic parameters may be estimated sequentially (i.e. one at a time) or, preferably, various parameters may be estimated simultaneously. It is also possible to use one type of data to estimate one or multiple parameters. In another embodiment, all available input data are used to estimate all anisotropic parameters simultaneously. Further, the values of the anisotropic model parameters can be determined using various combinations of inversions with full sets or subsets of the input data and full sets or subsets of the output data. - The tomographic inversion of
procedure 25 is recursive and iterative, and may use three dimensional tilted transversely isotropic ray tracing to model three dimensional wave propagation. Input data preprocessing may include normalizing the input before performing the three dimensional TTI tomographic inversion. Input data preprocessing may also include, prior to performing the three dimensional TTI tomographic inversion, migrating the seismic data using 3D TTI prestack depth migration algorithms, sorting the migrated seismic data according to their subsurface location and their migration offset/angle into common-image-point (CIP) gathers, and selecting and quantifying residual moveouts in the common-image-point gathers domain. TTI prestack depth migration algorithms are known in the art. - The values of the anisotropic parameters are calculated by iteratively selecting values of Vp0, δ and η. Specifically, the values of the anisotropic parameters are calculated until (a) seismic image positions of subsurface structures in the geological volume of interest tie their spatial positions recognized in drilled wells and (b) residual moveouts in common-imaging-point gathers are minimum at every point in the model. The subsurface structures include rock boundaries.
- Referring now to
FIG. 2 , this figure shows aworkflow 100 for determining values of anisotropic parameters of a Titled Transversely Isotropic (TTI) earth model in accordance with an embodiment of the invention. The iterative workflow can begin atprocedure 105 where three dimensional seismic reflection data are measured. The workflow then proceeds toprocedure 110 where the seismic reflection data are used to determine an initial TTI earth model for the geological volume of interest. The initial TTI earth model may be determined using various means. For example, the initial TTI earth model can be determined from velocity information and various processing performed on measured data as well as knowledge, either regionally or globally, about the area in which the geological volume of interest is located. Furthermore, the initial TTI model building also includes estimating subsurface reflector structural dips and obtaining TTI symmetry axes in the geological volume of interest. In the initial TTI earth model, a value of each parameter Vp0, δ and η is provided at each of a plurality of locations (x, y, z) in the geological volume of interest. - After determining an initial TTI earth model, the workflow proceeds to block 115 and
procedure 120 where well data (block 115) are used to perform a three dimensional (3D) TTI tomography (procedure 120) based on ray tracing to update the value of Vp0 of the TTI earth model near the well(s). Well data may include VSP data or checkshot data, which represent information about wave velocities in the well(s). The updated values of Vp0 in the geological volume of interest near the well(s) corresponds to Vp0(initial)+ΔVp0. With this operation, the values of Vp0 near the well(s) are updated in the initial TTI earth model. Specifically, the VSP data and checkshot data are used in the embodiment ofFIG. 2 to tie the seismic data to the well data. - After updating the values of Vp0 near the well(s), the
workflow 100 proceeds toprocedure 125 where the initial values of δ near the well(s) are updated using the updated values of Vp0 according to the following transformation Δδ=ΔVp0/Vp0 with δ(updated)=δ(initial)+Δδ. The transformations of procedures 110-125 enable one to refine the initial 3D TTI earth model by providing a more accurate TTI model near the well(s). - Once the initial TTI earth model has been updated near the well(s) with the updated values of Vp0 and δ, the
workflow 100 proceeds toprocedure 130 where the remaining points in the model, i.e. the points outside the location(s) near the well(s), are populated. In the embodiment ofFIG. 2 , the remaining points in the model are populated using a three dimensional extrapolation process and the updated values of Vp0 and δ near the well(s). The three dimensional extrapolation process is performed by first extrapolating δ or Δδ in the remaining locations of the model and then obtaining a three dimensional update of Vp0 using the transformation ΔVp0=−Δδ*Vp0. It will be appreciated that various types of algorithms can be used to extrapolate the values of δ or Δδ outside the well locations and to obtain extended updates of Vp0 and δ. The result of the extrapolation process (see block 135) provides a three dimensional TTI earth model with well tie having updated values of Vp0 and δ at each point of the model. The three dimensional TTI earth model ofblock 135 includes the following parameters at each point (x, y, z) of the model: Vp0(initial)+ΔVp0, δ(initial)+Δδ and η(initial). - The model determined at
block 135 is used to perform a three dimensional (3D) pre-stack TTI depth imaging of seismic reflection data atprocedure 140. With this procedure, seismic reflection data that have been recorded for the geological volume of interest are processed/migrated by the model to provide a new image of the sub-surface of the geological volume of interest. As will be appreciated, the migration of seismic data will typically enable an image to be formed of the geological volume of interest from the migrated seismic data that is a more accurate depiction of the geological features present in the seismic volume of interest than an image formed from migrated data with an isotropic algorithm. The three dimensional pre-stack TTI depth imaging of seismic reflection data ofprocedure 140 provides common image point gathers (block 145). As known in the art, common image point gathers correspond to the migrated seismic data that end up at a same image position in the geological volume of interest. - Once the common image point gathers have been identified by the 3D pre-stack TTI depth imaging analysis, the
workflow 100 then proceeds toprocedures 150 and/or 155 wherein a residual moveout analysis is performed. As known in the art, residual moveout analysis is a basic step in velocity model refinement. The analysis is generally carried out using semblances for a grid of inline and xline locations. One semblance panel is computed for each common-image-point gather as a function of depth and offset. The residual moveout analysis on the common image point gathers is carried out for the near to mid offsets atprocedure 150 while the residual moveout analysis on the common image point gathers is carried out for the near to far mid offsets atprocedure 155. Near and far offsets correspond, respectively, to a small and a large distance between the source and the receiver. - It will be appreciated that the
procedures FIG. 2 are intended to cover multiple moveout analyses. That is, the residual moveout analysis ofFIG. 2 can be performed only for the near to mid offsets (procedure 150) or only for the near to far offsets (procedure 155) in embodiments of the invention. Alternatively, in an embodiment of the invention, and as shown inFIG. 2 , the residual moveout analysis can be performed for both the near to mid offsets and the near to far offsets simultaneously. - The residual moveout analysis provides updated wave velocities Vp0 for each point of the 3D TTI model defined at
procedure 135. After determining the updated velocities Vp0 for each point of the model, a determination is made atblocks 165 a and/or 165 b as to whether convergence is obtained, i.e. whether the results of the residual moveout analyses ofprocedures 150 and/orprocedure 155 is/are within a predetermined tolerance. If convergence is obtained, theworkflow 100 ends atblock 166. - If convergence is not obtained, the
workflow 100 proceeds to block 170 and/or block 180 wherein it is determined whether the values of δ and/or η should be updated for each point of the model. If the result of the inquiry is positive atblock 170, the workflow proceeds toprocedure 185 where a three dimensional (3D) TTI tomography is performed to update the value of δ for each point of the model using the results of the residual moveout analysis for the near to mid offsets. It is greatly beneficial to update the values of δ only because such an update does not modify the values of wave velocities Vp0 in the model. As a result, the tie between well and seismic data that was previously obtained with the model constructed atblock 135 is not qualitatively destroyed in a low dip setting. In that way, it is possible to obtain a faster and more robust conversion for the values of the TTI model. - Alternatively or additionally, it is determined whether the values of η should be updated for each point of the model. If the result of the inquiry is positive at
block 180, the workflow proceeds toprocedure 195 where a three dimensional (3D) TTI tomography is performed to update the value of η for each point of the model using the results of the residual moveout analysis for the far to mid offsets. - If the values of δ and/or η are not updated, it is determined whether the values of wave velocities Vp0 should be updated (block 175). If the result of the inquiry is negative, the workflow ends at
block 166. If the result of the inquiry is positive, the workflow proceeds toprocedure 190 where a 3D TTI tomography is performed to update the values of velocities Vp0 in the geological volume of interest. - The results of the 3D TTI tomography analyses of
procedures block 197 with updated values of δ, Vp0 and η at each point (x, y, z) of the model. It is then determined whether convergence is obtained for the updated values of δ, Vp0 and η in the new model of block 197 (block 198). Various tests may be used to determine whether convergence is obtained. For example, it is determined whether the obtained values are below a predetermined threshold. If the result of the inquiry is positive, the workflow ofFIG. 2 ends atprocedure 199. If the result of the inquiry is negative, theworkflow 100 proceeds back to block 115 where the well data are used to perform a three dimensional (3D) TTI tomography (procedure 120) based on ray tracing to update the value of Vp0 of the TTI earth model obtained atblock 197 near the well(s). Thus, the TTI earth model obtained atblock 197 is used as a new initial model for the subsequent iteration. In an embodiment, this new initial TTI earth model is more refined than the initial model used at the first iteration in that it substantially flattens common-imaging-point gathers and substantially ties seismic data to well data. After determining a modified model with updated values of Vp0 near the well(s), theworkflow 100 proceeds to procedures 125-198 where updated values of δ and/or Vp0 and/or η are determined. The workflow is then iterated until convergence is obtained for the values of δ and/or Vp0 and/or η (i.e. the values of δ and/or Vp0 and/or η substantially do not change between two subsequent iterations). - It will be appreciated that the workflow of
FIG. 2 is intended to encompass several scenarios for optimizing the three dimensional (3D) TTI model. Referring now toFIG. 3 , this figure shows the various scenarios that can be applied at each iteration. In the first scenario, δ is the only parameter that is updated at each location in the model. The first scenario is preferred because the update of δ does not change the values of Vp0 in a low dip setting and, as a result, the quantitative tie between the well data and the seismic data is not destroyed during the optimization of δ. In the second scenario, parameters δ and Vp0 are updated. In the third scenario, parameters δ and η are updated. In the fourth scenario, the velocities Vp0 are updated. In the fifth scenario, Vp0 and η are updated. In the sixth scenario, η is updated and, in the seventh scenario, δ, η, and Vp0 are updated. - It will be appreciated that the one or more parameters updated at a given iteration may not be the same as the one or more parameters updated at a subsequent iteration. In other words, the anisotropic parameters (δ, η and Vp0) optimized at each iteration may be different. Thus, in an embodiment, it is envisioned that various scenarios could be used to optimize the three dimensional TTI model. Further it is envisioned that after optimizing a first parameter, e.g. δ according to the first scenario, the workflow of
FIG. 2 may be pursued to optimize a second parameter, e.g. η, and then a third parameter, e.g. Vp0. Therefore, the parameters δ, η and Vp0 can be optimized simultaneously or sequentially. Further, it will be appreciated that the selection of the parameters updated at each iteration is highly dependent on the maturity of the overall model building process that includes both depth imaging and tomography and the type and the quality of well and seismic data that are used to construct the 3D TTI model. - Although the invention has been described in detail for the purpose of illustration based on what is currently considered to be the most practical and preferred embodiments, it is to be understood that such detail is solely for that purpose and that the invention is not limited to the disclosed embodiments, but, on the contrary, is intended to cover modifications and equivalent arrangements that are within the spirit and scope of the appended claims.
- It will be appreciated that the different acts involved in determining values of anisotropic model parameters of a Tilted Transversely Isotropic (TTI) Earth model may be executed according to machine executable instructions or codes. These machine executable instructions may be embedded in a data storage medium. A processor may be configured to execute the instructions.
- Software functionalities of a computer system involving programming, including executable codes, may be used to implement the above described model. The software code may be executable by a general-purpose computer. In operation, the code and possibly the associated data records may be stored within a general-purpose computer platform. At other times, however, the software may be stored at other locations and/or transported for loading into an appropriate general-purpose computer-system. Hence, the embodiments discussed above involve one or more software or computer products in the form of one or more modules of code carried by at least one machine-readable medium. Execution of such codes by a processor of the computer system enables the platform to implement the functions in essentially the manner performed in the embodiments discussed and illustrated herein.
- As used herein, terms such as computer or machine “Readable medium” refer to any medium that participates in providing instructions to a processor for execution. A computer or a machine “readable medium” may be broadly termed a “computer product.” Such a medium may take many forms, including but not limited to, non-volatile media, volatile media, and transmission media. Non-volatile media include, for example, optical or magnetic disks, such as any of the storage devices in any computers) operating as discussed above. Volatile media include dynamic memory, such as the main memory of a computer system. Physical transmission media include coaxial cables, copper wires and fiber optics, including the wires that comprise a bus within a computer system. Carrier-wave transmission media can take the form of electric or electromagnetic signals, or acoustic or light waves such as those generated during radio frequency (RF) and infrared (IR) data communications. Common forms of computer-readable media therefore include, for example: a floppy disk, a flexible disk, hard disk, magnetic tape, any other magnetic medium, a CD-ROM, DVD, any other optical medium, less commonly used media such as punch cards, paper tape, any other physical medium with patterns of holes, a RAM, a PROM, and EPROM, a FLASH-EPROM, any other memory chip or cartridge, a carrier wave transporting data or instructions, cables or links transporting such a carrier wave, or any other medium from which a computer can read or send programming codes and/or data. Many of these forms of computer readable media maybe involved in carrying one or more sequences of one or more instructions to a processor for execution.
- It is to be understood that the present invention contemplates that, to the extent possible, one or more features of any embodiment can be combined with one or more features of any other embodiment.
Claims (16)
1. A computer implemented method for determining values of anisotropic model parameters of a Tilted Transversely Isotropic (TTI) Earth model, the anisotropic parameters including P-wave velocity (Vp0) along a tilted symmetry axis, the Thomsen anisotropy parameters δ and δ (or η=(ε−δ)/(1+2δ)) representative of variations of wave velocities as a function of wave propagation angle from the symmetry axis, the method comprising:
acquiring input data for a geological volume of interest;
determining a theoretical relationship between the input data and the anisotropic model parameters; and
calculating the values of the anisotropic model parameters at each of a plurality of subsurface locations in the geological volume of interest based on the theoretical relationships and the input data using workflows involving iterative or sequential combinations of processes including input data preprocessing, conventional tomographic inversion, three dimensional tomographic inversion based on a tilted transversely isotropic model, and three dimensional pre-stack depth migration using a tilted transversely isotropic model, thereby providing an accurate model of geological features present in the geological volume of interest.
2. The method of claim 1 , wherein the iterative or sequential combinations of processes includes interpretive picking.
3. The method of claim 1 , wherein the tomographic inversion is recursive and iterative, and uses three dimensional tilted transversely isotropic ray tracing to model three dimensional wave propagation.
4. The method of claim 1 , wherein the input data include surface seismic data, vertical seismic profile data, checkshot data, VSP data, well log data, interpretational data, regional trend, a priori data, or any combination of the foregoing.
5. The method of claim 4 , wherein the interpretational data include picks of horizon positions and shapes.
6. The method of claim 1 , further comprising creating an initial TTI model using the input data, the initial TTI model being created by estimating subsurface reflector structural dips and obtaining TTI symmetry axes in the geological volume of interest.
7. The method of claim 1 , further comprising normalizing the input data prior to performing the inversion.
8. The method of claim 1 , wherein the input data includes seismic data, the method further comprising, prior to performing the three dimensional TTI tomographic inversion, migrating the seismic data using 3D TTI prestack depth migration algorithm, sorting the migrated seismic data according to their subsurface location and their migration offset/angle into common-image-point (CIP) gathers, and selecting and quantifying residual moveouts in the common-image-point gathers domain.
9. The method of claim 1 , wherein the values of the anisotropic model parameters at each of the plurality of subsurface locations are calculated using three dimensional TTI tomographic inversions, simultaneously, or separately in a sequential manner, using various combinations of inversions with full sets or subsets of the input data and full sets or subsets of output data.
10. The method of claim 9 , wherein the output data includes Vp0, δ and η and/or ε.
11. The method of claim 1 further comprising performing a TTI pre-stack depth migration process using a TTI model with Vp0, δ and η.
12. The method of claim 1 , wherein the calculation includes calculating values of the anisotropic parameters by iteratively selecting values of Vp0, δ and η until (a) seismic image positions of subsurface structures in the geological volume of interest tie their spatial positions recognized in drilled wells and (b) residual moveouts in common-imaging-point gathers are minimum at every point in the model.
13. The method of claim 12 , wherein the subsurface structures include rock boundaries.
14. A computer implemented method for determining values of anisotropic model parameters of a three dimensional TTI earth model, the anisotropic parameters including P-wave velocity (Vp0) along a tilted symmetry axis, the Thomsen anisotropy parameters δ, ε (or η=(ε−δ)/(1+2δ)) representative of variation of wave velocities as a function of wave propagation angle from the symmetry axis, the method comprising:
a. obtaining an initial TTI earth model that substantially flattens common-imaging-point gathers and substantially ties seismic data to well data, the initial migration velocity model including initial values Vp0, δ, ε (or η) at each of a plurality of subsurface locations in a geological volume of interest;
b. inputting checkshot data, or VSP data, or both checkshot and VSP data at well locations into a three dimensional tomographic inversion to determine updated values of Vp0 near the well locations, the values of Vp0 being updated by a correction ΔVp0, wherein Vp0=Vp0+ΔVp0;
c. determining an incremental improvement Δδ to δ(initial) using the relative change Δδ=(ΔVp0)/Vp0;
d. extrapolating the relative change Δδ from near-well locations to the entire three dimensional TTI earth model at each of the plurality of subsurface locations, taking into account geological consistency and regularization, to determine updated values of δ, wherein δ=δ+Δδ;
e. determining updated values of Vp0=Vp0(1−Δδ) using the extrapolated three dimensional Δδ at each of the plurality of subsurface locations and obtaining the three dimensional extended incremental update ΔVp0=−Δδ Vp0;
f. inputting near-to-mid-offset/angle residual moveout information in the common-imaging-point gathers obtained with an improved migration velocity model that is defined with the updated values of Vp0 and δ into a TTI tomographic inversion process to further provide updated values of δ at each of the plurality of subsurface locations; and
g. inputting near-to-far-offset/angle residual moveout information in the common-imaging-point gathers obtained with the improved model the TTI tomographic inversion process to provide updated values of η at each of the plurality of subsurface locations.
15. The method of claim 14 , further comprising iteratively repeating (a), (b), (c), (d), (e), (f) and (g), wherein the improved parameters of the TTI model determined at the end of an iteration corresponds to the initial parameters of the TTI model at the next iteration.
16. A computer product having machine executable instructions, the instructions being executable by a machine to perform a tomographic inversion method for determining values of anisotropic parameters of a TTI earth model, the anisotropic parameters including P-wave velocity (Vp0) along a tilted symmetry axis, the Thomsen anisotropy parameters δ and ε (or η=(ε−δ)/(1+2δ)) representative of variation of wave velocities as a function of wave propagation angle from the symmetry axis, the method comprising:
determining a relationship between input data and the anisotropic parameters, the input data being acquired for a geological volume of interest; and
calculating the values of the anisotropic parameters at each of a plurality of subsurface locations in the geological volume of interest based on the relationship and the input data using workflows involving iterative or sequential combinations of processes including input data preprocessing, three dimensional tomographic inversion, and three dimensional TTI pre-stack depth migration.
Priority Applications (10)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US12/327,671 US20100135115A1 (en) | 2008-12-03 | 2008-12-03 | Multiple anisotropic parameter inversion for a tti earth model |
SG2011040144A SG171913A1 (en) | 2008-12-03 | 2009-10-22 | Multiple anisotropic parameter inversion for a tti earth model |
CA2745076A CA2745076A1 (en) | 2008-12-03 | 2009-10-22 | Multiple anisotropic parameter inversion for a tti earth model |
EA201170747A EA021972B1 (en) | 2008-12-03 | 2009-10-22 | Multiple anisotropic parameter inversion for a tilted transversely isotropic (tti) earth model |
BRPI0922702A BRPI0922702A2 (en) | 2008-12-03 | 2009-10-22 | computer-implemented method for determining anisotropic model parameter values of a soil model, and, computer product |
AU2009322868A AU2009322868A1 (en) | 2008-12-03 | 2009-10-22 | Multiple anisotropic parameter inversion for a TTI earth model |
CN200980148505XA CN102239429B (en) | 2008-12-03 | 2009-10-22 | Multiple anisotropic parameter inversion for a tti earth model |
EP09830791A EP2370924A2 (en) | 2008-12-03 | 2009-10-22 | Multiple anisotropic parameter inversion for a tti earth model |
PCT/US2009/061612 WO2010065209A2 (en) | 2008-12-03 | 2009-10-22 | Multiple anisotropic parameter inversion for a tti earth model |
US13/368,134 US8666668B2 (en) | 2008-12-03 | 2012-02-07 | Multiple anisotropic parameter inversion for a TTI earth model using well data |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US12/327,671 US20100135115A1 (en) | 2008-12-03 | 2008-12-03 | Multiple anisotropic parameter inversion for a tti earth model |
Related Child Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US13/368,134 Division US8666668B2 (en) | 2008-12-03 | 2012-02-07 | Multiple anisotropic parameter inversion for a TTI earth model using well data |
Publications (1)
Publication Number | Publication Date |
---|---|
US20100135115A1 true US20100135115A1 (en) | 2010-06-03 |
Family
ID=42222700
Family Applications (2)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US12/327,671 Abandoned US20100135115A1 (en) | 2008-12-03 | 2008-12-03 | Multiple anisotropic parameter inversion for a tti earth model |
US13/368,134 Expired - Fee Related US8666668B2 (en) | 2008-12-03 | 2012-02-07 | Multiple anisotropic parameter inversion for a TTI earth model using well data |
Family Applications After (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US13/368,134 Expired - Fee Related US8666668B2 (en) | 2008-12-03 | 2012-02-07 | Multiple anisotropic parameter inversion for a TTI earth model using well data |
Country Status (9)
Country | Link |
---|---|
US (2) | US20100135115A1 (en) |
EP (1) | EP2370924A2 (en) |
CN (1) | CN102239429B (en) |
AU (1) | AU2009322868A1 (en) |
BR (1) | BRPI0922702A2 (en) |
CA (1) | CA2745076A1 (en) |
EA (1) | EA021972B1 (en) |
SG (1) | SG171913A1 (en) |
WO (1) | WO2010065209A2 (en) |
Cited By (35)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102141633A (en) * | 2010-12-10 | 2011-08-03 | 中国科学院地质与地球物理研究所 | Anisotropic three-dimensional prestack time migration method |
US20120059633A1 (en) * | 2010-07-08 | 2012-03-08 | Nader Dutta | Methods and Devices for Transformation of Collected Data for Improved Visualization Capability |
US20120099396A1 (en) * | 2010-10-22 | 2012-04-26 | Chevron U.S.A. Inc. | System and method for characterization with non-unique solutions of anisotropic velocities |
US20120195165A1 (en) * | 2011-01-31 | 2012-08-02 | Chevron U.S.A. Inc. | Exploitation of self-consistency and differences between volume images and interpreted spatial/volumetric context |
CN102721977A (en) * | 2012-05-31 | 2012-10-10 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | High-precision converted wave anisotropy stacking velocity analysis method |
FR2976088A1 (en) * | 2011-05-30 | 2012-12-07 | Cggveritas Services Sa | NON-LINEAR TOMOGRAPHY METHOD FOR A MAIN SYMMETRY AXIS OF AN ANISOTROPIC SPEED MODEL AND DEVICE |
CN102830433A (en) * | 2012-08-16 | 2012-12-19 | 中国石油化工股份有限公司 | Method for calculating anisotropy azimuth angle by using dipole transverse wave logging information based on frequency domain |
US20130060544A1 (en) * | 2010-05-12 | 2013-03-07 | Petrus Maria Bakker | Seismic p-wave modelling in an inhomogeneous transversely isotropic medium with a tilted symmetry axis |
US20130235693A1 (en) * | 2012-03-09 | 2013-09-12 | Chevron U.S.A. Inc | Correction of shear log for elastic anisotropy |
CN103352691A (en) * | 2013-07-05 | 2013-10-16 | 天津大学 | Orthogonal dipole acoustic logging sound system receiving device |
WO2014204440A1 (en) * | 2013-06-18 | 2014-12-24 | Halliburton Energy Services, Inc. | Methods and systems for seismic data analysis using a tilted transversely isotropic (tti) model |
WO2015164917A1 (en) * | 2014-05-01 | 2015-11-05 | Curtin University Of Technology | A method of estimating a volume fraction of coal in a layered geological structure |
US20160187512A1 (en) * | 2014-01-14 | 2016-06-30 | Cgg Services Sa | Detecting and estimating anisotropy errors using full waveform inversion and ray based tomography |
RU2596182C1 (en) * | 2012-08-17 | 2016-08-27 | Лэндмарк Графикс Корпорейшн | System and method for displaying seismic data |
US9442207B1 (en) * | 2015-06-12 | 2016-09-13 | Chevron U.S.A. Inc. | System and method for computing residual moveout from seismic images |
US9562983B2 (en) * | 2014-04-17 | 2017-02-07 | Saudi Arabian Oil Company | Generating subterranean imaging data based on vertical seismic profile data |
CN106896408A (en) * | 2017-03-23 | 2017-06-27 | 中国石油天然气股份有限公司 | A kind of angle domain prestack time migration method |
US20170350245A1 (en) * | 2015-05-01 | 2017-12-07 | Halliburton Energy Services, Inc. | Anisotropic parameter estimation from walkaway vsp data using differential evolution |
CN108710148A (en) * | 2018-05-29 | 2018-10-26 | 中国科学院地质与地球物理研究所 | The steady phase prestack depth migration method in three-dimensional dip domain and device |
US10267937B2 (en) | 2014-04-17 | 2019-04-23 | Saudi Arabian Oil Company | Generating subterranean imaging data based on vertical seismic profile data and ocean bottom sensor data |
US20190235113A1 (en) * | 2016-10-25 | 2019-08-01 | Landmark Graphics Corporation | Estimating interval anisotropy parameter for pre-stack depth migration using a least-sqares method |
CN111596346A (en) * | 2019-02-20 | 2020-08-28 | 中国石油天然气集团有限公司 | Elastic wave velocity inversion method and device |
CN111624647A (en) * | 2020-06-05 | 2020-09-04 | 中油奥博(成都)科技有限公司 | Integrated prestack time migration method and device for variable offset VSP ray tracing |
CN112130202A (en) * | 2020-08-14 | 2020-12-25 | 中国石油天然气集团有限公司 | Method and system for inverting orthotropic velocity |
US10901103B2 (en) | 2018-03-20 | 2021-01-26 | Chevron U.S.A. Inc. | Determining anisotropy for a build section of a wellbore |
CN112305589A (en) * | 2020-09-22 | 2021-02-02 | 中国石油天然气集团有限公司 | Method and device for imaging depth domain of anisotropic medium |
CN112444861A (en) * | 2019-08-27 | 2021-03-05 | 中国石油化工股份有限公司 | Speed model updating method, computer storage medium and computer system |
CN112540407A (en) * | 2020-12-01 | 2021-03-23 | 中国煤炭地质总局地球物理勘探研究院 | Method for establishing prestack depth migration anisotropic field |
CN113009580A (en) * | 2019-12-20 | 2021-06-22 | 中国石油天然气集团有限公司 | Variable offset VSP first-break inversion VTI anisotropic parameter method and device |
CN113221228A (en) * | 2021-06-04 | 2021-08-06 | 中国电建集团成都勘测设计研究院有限公司 | Hydropower station underground cave group surrounding rock mechanical parameter inversion method |
CN113296153A (en) * | 2020-02-24 | 2021-08-24 | 中国石油天然气集团有限公司 | Method and device for determining anisotropy parameters of axisymmetric medium |
CN113568045A (en) * | 2021-07-26 | 2021-10-29 | 中南大学 | Inversion method of TTI medium crack inclination angle |
CN114594515A (en) * | 2020-12-07 | 2022-06-07 | 中国石油化工股份有限公司 | Well control speed inversion method based on slowly varying anisotropy |
CN115951405A (en) * | 2023-03-08 | 2023-04-11 | 中国地震局地球物理研究所 | Speed change imaging method based on sparse convolution and adaptive moment estimation |
RU2806537C1 (en) * | 2023-04-04 | 2023-11-01 | Акционерное общество "Сибирский научно-исследовательский институт геологии, геофизики и минерального сырья" | Seismic method for identifying vertical and steeply inclined layers by reverse scanning on longitudinal refracted waves |
Families Citing this family (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CA2906362A1 (en) * | 2013-04-02 | 2014-10-09 | Halliburton Energy Services, Inc. | Anisotropy analysis using direct and reflected arrivals in seismic survey data |
US20140301165A1 (en) * | 2013-04-03 | 2014-10-09 | Westerngeco L.L.C. | Seismic data processing using joint tomography |
CN104216011B (en) * | 2013-06-05 | 2017-08-04 | 上海青凤致远地球物理地质勘探科技有限公司 | A kind of stable qP ripple reverse-time migration methods of TTI media |
DE112013007359T5 (en) * | 2013-08-23 | 2016-05-25 | Landmark Graphics Corporation | Local update of a geocellular 3D model |
CN105493100A (en) * | 2013-08-29 | 2016-04-13 | 兰德马克绘图国际公司 | Static earth model calibration methods and systems |
CN104570119B (en) * | 2013-10-29 | 2017-11-17 | 中国石油化工股份有限公司 | A kind of three-dimensional perpendicular seismic profile back wave stretches bearing calibration |
CN104076391B (en) * | 2014-04-16 | 2015-12-02 | 孙学凯 | Based on the Local angle domain anisotropy offset method of TTI medium quadravalence travel-time equation |
CN107966729B (en) * | 2016-10-19 | 2019-07-19 | 中国石油化工股份有限公司 | A kind of three-dimensional TTI medium ray-tracing procedure and system |
CN106709031A (en) * | 2016-12-29 | 2017-05-24 | 天津瀚海星云数字科技有限公司 | High-precision 3D GIS Earch implementing method and system |
CN108732623B (en) * | 2017-04-17 | 2020-05-19 | 中国石油化工股份有限公司 | TTI medium tomography inversion imaging method and computer readable storage medium |
CN109581521B (en) * | 2017-09-28 | 2020-12-01 | 中国石油化工股份有限公司 | Local chromatographic method and system for TTI anisotropy |
CN109655889B (en) * | 2017-10-11 | 2021-04-02 | 中国石油化工股份有限公司 | Anisotropy parameter joint inversion method and system |
CN109581499A (en) * | 2018-11-08 | 2019-04-05 | 成都捷科思石油天然气技术发展有限公司 | A method of structural map is generated using anisotropy pre-stack depth migration |
CN110161562B (en) * | 2019-04-29 | 2020-09-18 | 中南大学 | Crack weakness inversion method and system in inclined transverse isotropic medium |
EP3859405A1 (en) * | 2020-01-28 | 2021-08-04 | Repsol Exploración, S.A. | Computer-implemented method for determining a velocity image of a domain of the subsurface structural geology in an oil and gas reservoir |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6714480B2 (en) * | 2002-03-06 | 2004-03-30 | Schlumberger Technology Corporation | Determination of anisotropic moduli of earth formations |
US6785612B1 (en) * | 2003-05-29 | 2004-08-31 | Pgs Americas, Inc. | Seismic velocity update for anisotropic depth migration |
US6920084B2 (en) * | 2003-01-15 | 2005-07-19 | Western Geco, L.L.C. | Offset transformation to zero dip that preserves angle of incidence |
US7072767B2 (en) * | 2003-04-01 | 2006-07-04 | Conocophillips Company | Simultaneous inversion for source wavelet and AVO parameters from prestack seismic data |
US20080273421A1 (en) * | 2007-03-05 | 2008-11-06 | Zvi Koren | Model-based time-preserving tomography |
US20100128562A1 (en) * | 2008-11-25 | 2010-05-27 | Baker Hughes Incorporated | Anisotropic Parameter Determination |
US20100312534A1 (en) * | 2008-02-28 | 2010-12-09 | Shiyu Xu | Rock Physics Model For Simulating Seismic Response In Layered Fractured Rocks |
US7924652B2 (en) * | 2007-06-01 | 2011-04-12 | Baker Hughes Incorporated | Method for determining seismic anisotropy |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8798967B2 (en) * | 2011-03-30 | 2014-08-05 | Chevron U.S.A. Inc. | System and method for computations utilizing optimized earth model representations |
-
2008
- 2008-12-03 US US12/327,671 patent/US20100135115A1/en not_active Abandoned
-
2009
- 2009-10-22 SG SG2011040144A patent/SG171913A1/en unknown
- 2009-10-22 CN CN200980148505XA patent/CN102239429B/en not_active Expired - Fee Related
- 2009-10-22 WO PCT/US2009/061612 patent/WO2010065209A2/en active Application Filing
- 2009-10-22 EP EP09830791A patent/EP2370924A2/en not_active Withdrawn
- 2009-10-22 AU AU2009322868A patent/AU2009322868A1/en not_active Abandoned
- 2009-10-22 CA CA2745076A patent/CA2745076A1/en not_active Abandoned
- 2009-10-22 BR BRPI0922702A patent/BRPI0922702A2/en not_active IP Right Cessation
- 2009-10-22 EA EA201170747A patent/EA021972B1/en not_active IP Right Cessation
-
2012
- 2012-02-07 US US13/368,134 patent/US8666668B2/en not_active Expired - Fee Related
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6714480B2 (en) * | 2002-03-06 | 2004-03-30 | Schlumberger Technology Corporation | Determination of anisotropic moduli of earth formations |
US6920084B2 (en) * | 2003-01-15 | 2005-07-19 | Western Geco, L.L.C. | Offset transformation to zero dip that preserves angle of incidence |
US7072767B2 (en) * | 2003-04-01 | 2006-07-04 | Conocophillips Company | Simultaneous inversion for source wavelet and AVO parameters from prestack seismic data |
US6785612B1 (en) * | 2003-05-29 | 2004-08-31 | Pgs Americas, Inc. | Seismic velocity update for anisotropic depth migration |
US20080273421A1 (en) * | 2007-03-05 | 2008-11-06 | Zvi Koren | Model-based time-preserving tomography |
US7924652B2 (en) * | 2007-06-01 | 2011-04-12 | Baker Hughes Incorporated | Method for determining seismic anisotropy |
US20100312534A1 (en) * | 2008-02-28 | 2010-12-09 | Shiyu Xu | Rock Physics Model For Simulating Seismic Response In Layered Fractured Rocks |
US20100128562A1 (en) * | 2008-11-25 | 2010-05-27 | Baker Hughes Incorporated | Anisotropic Parameter Determination |
Cited By (46)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9285491B2 (en) * | 2010-05-12 | 2016-03-15 | Shell Oil Company | Seismic P-wave modelling in an inhomogeneous transversely isotropic medium with a tilted symmetry axis |
US20130060544A1 (en) * | 2010-05-12 | 2013-03-07 | Petrus Maria Bakker | Seismic p-wave modelling in an inhomogeneous transversely isotropic medium with a tilted symmetry axis |
US20120059633A1 (en) * | 2010-07-08 | 2012-03-08 | Nader Dutta | Methods and Devices for Transformation of Collected Data for Improved Visualization Capability |
US8914269B2 (en) * | 2010-07-08 | 2014-12-16 | Westerngeco L.L.C. | Methods and devices for transformation of collected data for improved visualization capability |
US20120099396A1 (en) * | 2010-10-22 | 2012-04-26 | Chevron U.S.A. Inc. | System and method for characterization with non-unique solutions of anisotropic velocities |
CN102141633A (en) * | 2010-12-10 | 2011-08-03 | 中国科学院地质与地球物理研究所 | Anisotropic three-dimensional prestack time migration method |
US8861309B2 (en) * | 2011-01-31 | 2014-10-14 | Chevron U.S.A. Inc. | Exploitation of self-consistency and differences between volume images and interpreted spatial/volumetric context |
US20120195165A1 (en) * | 2011-01-31 | 2012-08-02 | Chevron U.S.A. Inc. | Exploitation of self-consistency and differences between volume images and interpreted spatial/volumetric context |
US9188690B2 (en) | 2011-05-30 | 2015-11-17 | Cggveritas Services Sa | Non-linear tomography method for main axis of symmetry of anisotropic velocity model and apparatus |
FR2981170A1 (en) * | 2011-05-30 | 2013-04-12 | Cggveritas Services Sa | NONLINEAR TOMOGRAPHY METHOD FOR A MAIN SYMMETRY AXIS OF AN ANISOTROPIC SPEED MODEL AND DEVICE. |
FR2976088A1 (en) * | 2011-05-30 | 2012-12-07 | Cggveritas Services Sa | NON-LINEAR TOMOGRAPHY METHOD FOR A MAIN SYMMETRY AXIS OF AN ANISOTROPIC SPEED MODEL AND DEVICE |
US9207342B2 (en) * | 2012-03-09 | 2015-12-08 | Chevron U.S.A. Inc. | Correction of shear log for elastic anisotropy |
US20130235693A1 (en) * | 2012-03-09 | 2013-09-12 | Chevron U.S.A. Inc | Correction of shear log for elastic anisotropy |
CN102721977A (en) * | 2012-05-31 | 2012-10-10 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | High-precision converted wave anisotropy stacking velocity analysis method |
CN102830433A (en) * | 2012-08-16 | 2012-12-19 | 中国石油化工股份有限公司 | Method for calculating anisotropy azimuth angle by using dipole transverse wave logging information based on frequency domain |
RU2596182C1 (en) * | 2012-08-17 | 2016-08-27 | Лэндмарк Графикс Корпорейшн | System and method for displaying seismic data |
EP2997401A4 (en) * | 2013-06-18 | 2017-02-22 | Halliburton Energy Services, Inc. | Methods and systems for seismic data analysis using a tilted transversely isotropic (tti) model |
US10768323B2 (en) * | 2013-06-18 | 2020-09-08 | Halliburton Energy Services, Inc. | Methods and systems for seismic data analysis using a tilted transversely isotropic (TTI) model |
US20160091622A1 (en) * | 2013-06-18 | 2016-03-31 | Halliburton Energy Services, Inc. | Methods and systems for seismic data analysis using a tilted transversely isotropic (tti) model |
WO2014204440A1 (en) * | 2013-06-18 | 2014-12-24 | Halliburton Energy Services, Inc. | Methods and systems for seismic data analysis using a tilted transversely isotropic (tti) model |
AU2013392659B2 (en) * | 2013-06-18 | 2017-06-15 | Halliburton Energy Services, Inc. | Methods and systems for seismic data analysis using a tilted transversely isotropic (TTI) model |
CN103352691A (en) * | 2013-07-05 | 2013-10-16 | 天津大学 | Orthogonal dipole acoustic logging sound system receiving device |
US20160187512A1 (en) * | 2014-01-14 | 2016-06-30 | Cgg Services Sa | Detecting and estimating anisotropy errors using full waveform inversion and ray based tomography |
US9562983B2 (en) * | 2014-04-17 | 2017-02-07 | Saudi Arabian Oil Company | Generating subterranean imaging data based on vertical seismic profile data |
US10267937B2 (en) | 2014-04-17 | 2019-04-23 | Saudi Arabian Oil Company | Generating subterranean imaging data based on vertical seismic profile data and ocean bottom sensor data |
WO2015164917A1 (en) * | 2014-05-01 | 2015-11-05 | Curtin University Of Technology | A method of estimating a volume fraction of coal in a layered geological structure |
US20170350245A1 (en) * | 2015-05-01 | 2017-12-07 | Halliburton Energy Services, Inc. | Anisotropic parameter estimation from walkaway vsp data using differential evolution |
US9442207B1 (en) * | 2015-06-12 | 2016-09-13 | Chevron U.S.A. Inc. | System and method for computing residual moveout from seismic images |
US10908309B2 (en) * | 2016-10-25 | 2021-02-02 | Landmark Graphics Corporation | Estimating interval anisotropy parameter for pre-stack depth migration using a least-squares method |
US20190235113A1 (en) * | 2016-10-25 | 2019-08-01 | Landmark Graphics Corporation | Estimating interval anisotropy parameter for pre-stack depth migration using a least-sqares method |
CN106896408A (en) * | 2017-03-23 | 2017-06-27 | 中国石油天然气股份有限公司 | A kind of angle domain prestack time migration method |
US10901103B2 (en) | 2018-03-20 | 2021-01-26 | Chevron U.S.A. Inc. | Determining anisotropy for a build section of a wellbore |
CN108710148A (en) * | 2018-05-29 | 2018-10-26 | 中国科学院地质与地球物理研究所 | The steady phase prestack depth migration method in three-dimensional dip domain and device |
CN111596346A (en) * | 2019-02-20 | 2020-08-28 | 中国石油天然气集团有限公司 | Elastic wave velocity inversion method and device |
CN112444861A (en) * | 2019-08-27 | 2021-03-05 | 中国石油化工股份有限公司 | Speed model updating method, computer storage medium and computer system |
CN113009580A (en) * | 2019-12-20 | 2021-06-22 | 中国石油天然气集团有限公司 | Variable offset VSP first-break inversion VTI anisotropic parameter method and device |
CN113296153A (en) * | 2020-02-24 | 2021-08-24 | 中国石油天然气集团有限公司 | Method and device for determining anisotropy parameters of axisymmetric medium |
CN111624647A (en) * | 2020-06-05 | 2020-09-04 | 中油奥博(成都)科技有限公司 | Integrated prestack time migration method and device for variable offset VSP ray tracing |
CN112130202A (en) * | 2020-08-14 | 2020-12-25 | 中国石油天然气集团有限公司 | Method and system for inverting orthotropic velocity |
CN112305589A (en) * | 2020-09-22 | 2021-02-02 | 中国石油天然气集团有限公司 | Method and device for imaging depth domain of anisotropic medium |
CN112540407A (en) * | 2020-12-01 | 2021-03-23 | 中国煤炭地质总局地球物理勘探研究院 | Method for establishing prestack depth migration anisotropic field |
CN114594515A (en) * | 2020-12-07 | 2022-06-07 | 中国石油化工股份有限公司 | Well control speed inversion method based on slowly varying anisotropy |
CN113221228A (en) * | 2021-06-04 | 2021-08-06 | 中国电建集团成都勘测设计研究院有限公司 | Hydropower station underground cave group surrounding rock mechanical parameter inversion method |
CN113568045A (en) * | 2021-07-26 | 2021-10-29 | 中南大学 | Inversion method of TTI medium crack inclination angle |
CN115951405A (en) * | 2023-03-08 | 2023-04-11 | 中国地震局地球物理研究所 | Speed change imaging method based on sparse convolution and adaptive moment estimation |
RU2806537C1 (en) * | 2023-04-04 | 2023-11-01 | Акционерное общество "Сибирский научно-исследовательский институт геологии, геофизики и минерального сырья" | Seismic method for identifying vertical and steeply inclined layers by reverse scanning on longitudinal refracted waves |
Also Published As
Publication number | Publication date |
---|---|
AU2009322868A1 (en) | 2010-06-10 |
US8666668B2 (en) | 2014-03-04 |
CA2745076A1 (en) | 2010-06-10 |
EA201170747A1 (en) | 2011-12-30 |
WO2010065209A3 (en) | 2010-08-12 |
EA021972B1 (en) | 2015-10-30 |
SG171913A1 (en) | 2011-07-28 |
CN102239429B (en) | 2013-11-20 |
BRPI0922702A2 (en) | 2016-01-05 |
WO2010065209A2 (en) | 2010-06-10 |
EP2370924A2 (en) | 2011-10-05 |
US20120185171A1 (en) | 2012-07-19 |
CN102239429A (en) | 2011-11-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US8666668B2 (en) | Multiple anisotropic parameter inversion for a TTI earth model using well data | |
US8880384B2 (en) | Artifact reduction in iterative inversion of geophysical data | |
CN108072892B (en) | Automatic geological structure constraint chromatography inversion method | |
US9022129B2 (en) | Tracking geologic object and detecting geologic anomalies in exploration seismic data volume | |
WO2017035104A1 (en) | Velocity model seismic static correction | |
Waheed et al. | First-arrival traveltime tomography for anisotropic media using the adjoint-state method | |
US10838092B2 (en) | Estimating multiple subsurface parameters by cascaded inversion of wavefield components | |
US10310117B2 (en) | Efficient seismic attribute gather generation with data synthesis and expectation method | |
RU2631407C1 (en) | Method and device for seismic signals processing | |
US20110098996A1 (en) | Sifting Models of a Subsurface Structure | |
US10739481B2 (en) | 2D multiline seismic reflection tomography with seismic-tie constraint | |
CN115877449B (en) | Computer-implemented method for obtaining subsurface superimposed images within a survey area | |
CN111123359B (en) | Logging while drilling and stratum grid constrained well periphery seismic imaging detection method and device | |
US20200271810A1 (en) | Hybrid Residual Moveout Error Estimation | |
US20150378039A1 (en) | System and method for seismic imaging of a complex subsurface | |
US20210018642A1 (en) | System and method for seismic data calibration | |
WO2021191722A1 (en) | System and method for stochastic full waveform inversion | |
CN114861515A (en) | Method, device, equipment and medium for calculating layer speed data volume | |
US10620341B2 (en) | System and method for modifying an earth model | |
US11867856B2 (en) | Method and system for reflection-based travel time inversion using segment dynamic image warping | |
WO2024067458A1 (en) | While-drilling vsp well-driven seismic imaging method and apparatus | |
Mildner | A Self-Consistent Marchenko Method for Enhanced Imaging with Redatumed Data | |
CN113267810A (en) | Seismic exploration full-depth velocity modeling method and device |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
AS | Assignment |
Owner name: CHEVRON U.S.A. INC,CALIFORNIA Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNORS:XU, TONG;SUN, YONGHE;ZHANG, LEONARD LIN;AND OTHERS;SIGNING DATES FROM 20090116 TO 20090120;REEL/FRAME:022172/0946 |
|
STCB | Information on status: application discontinuation |
Free format text: ABANDONED -- FAILURE TO RESPOND TO AN OFFICE ACTION |