Abstract: A method, having application to petroleum exploration or production, for picking the arrival time of seismic waves and use thereof for orienting the components of a multi-component sensor. After acquisition of seismic data using a VSP type method, with a multi-component sensor, a module signal is constructed by calculating the square root of the sum of the squares of at least two orthogonal seismic components. Arrival times of a direct seismic wave are then picked on an amplitude extremum of this module signal. Based on this picking, the seismic components can then be oriented in a unique reference frame whatever the depth of the sensor. A time window is defined on either side of the picked arrival times and the azimuthal direction is determined by maximizing the energy of the horizontal components within this time window. Finally, the three components are oriented in a reference frame defined with respect to the azimuthal direction, which is identical for each depth.

Abstract: This invention relates to a method for processing seismic data comprising a collection of seismic traces with different offsets, comprising the steps of: a) breaking up one or several first trace(s) of the collection of seismic traces into a series of trace segment's according to a predetermined segmentation interval; b) defining a series of expansion coefficients, each expansion coefficient being associated with a segment of the first trace or traces; c) applying the associated expansion coefficient to each segment of the first trace or traces; d) comparing the first trace or traces thus expanded with a second trace from the collection of seismic traces to evaluate their similarity; e) repeating steps b), and d) with a new series of expansion coefficients, f) determining an optimum series of expansion coefficients that maximizes the similarity between the first expanded trace and the second trace in order to obtain one or several corrected first trace(s).

Abstract: A stable method for using fat-ray tomography to determine a high-resolution velocity model of the subsurface from seismic data (71). The velocity model (72) may be used in migrating the seismic data (76) to image the subsurface. Rays are traced from a subsurface reflection point to surface source and receiver locations (73), using Fresnel zone construction methods (74) that honor correct initial conditions, with the Fresnel radius being a function of velocity.

Abstract: Seismic data are assembled and stored for a set of cross-correlation lag times to form an array of common image gathers over depth levels of interest. The two dimensional gathers assembled over different lag times form a three-dimensional cube of common image data. The data are analyzed to determine the travel time shift required to equalize upgoing and downgoing wavefields. Events in the common image gathers are then modeled using Green's functions to generate a data set representing the data resulting from processing had a precise velocity model been obtainable from the seismic data. The generated data are then processed with inversion techniques to form a velocity model for seismic data analysis.

Abstract: The invention relates to methods and computer-readable medium for computing the propagation velocity of seismic waves in the earth. The invention computes the true propagation velocity of seismic waves in the earth, which is a condition of obtaining an accurate image of subsurface geology that can be used to prospect for oil and gas deposits. In an embodiment, the method of computing the propagation velocity of seismic waves in earth, includes providing an estimate of the propagation velocity, generating a time shift gather using a depth migration at a plurality of locations of the earth, converting each of the time shift gathers to a semblance gather, transforming each semblance gather into a velocity gather whose energy peaks represent a root-mean-square average of the propagation velocity along the forward and backward path between earth's surface and a point of the subsurface geology, and converting the energy peaks to the propagation velocity.

Abstract: A method for processing seismic data, characterized in that with a view to analyzing them as a function of at least four dimensions. Picking and/or picking propagation is carried out on at least one subset of data which corresponds to a collection of data sampled in three or two dimensions. The quantity or quantities which correspond to the other dimension or dimensions are fixed values. In applying a given projection law to at least one of the points highlighted by the picking and/or picking propagation to determine a point that corresponds to it in another subset corresponding to a collection of data sampled in the same dimensions for a different value of at least one quantity that corresponds to at least one other dimension.

Abstract: A disclosed direct velocity seismic sensor includes a housing, a proof mass suspended in the housing by a resilient component, and a motion dampener that damps oscillation of the proof mass to a degree that displacement of the proof mass relative to the housing is substantially linearly proportional to a rate of change of seismic displacements of the housing over a frequency range of interest. A described method for constructing a seismic sensor includes using a calculated resonant frequency to determine a damping factor that causes the displacement of the proof mass to be substantially proportional to the rate of change of seismic displacement of the housing. One illustrative disclosed system includes an optical velocity sensor and a detector where a light beam produced by the velocity sensor and a reference beam interfere at the detector, and the detector produces a signal indicative of a velocity experienced by the velocity sensor.

Abstract: A calculated vertical velocity sensor signal is determined from a recorded pressure sensor signal. A constructed vertical velocity sensor signal is determined as a linear combination of the calculated vertical velocity sensor signal and a recorded vertical velocity sensor signal in dual-sensor seismic streamer data, using a mixture coefficient as a proportionality constant. An upgoing pressure wavefield component is determined as one half of a difference of the recorded pressure sensor signal and the constructed vertical velocity sensor signal, as a function of the mixture coefficient. An error in the upgoing pressure wavefield component is determined by propagating errors in the recorded pressure sensor signal and constructed vertical velocity sensor signal terms. A value of the mixture coefficient is determined that minimizes the error in the upgoing pressure wavefield component.

Abstract: A method for evaluating fluid pressures and detecting overpressures in an underground medium is disclosed having application to petroleum exploration for detection of overpressure zones while drilling for example. A seismic P wave velocity cube and a seismic S wave velocity cube are constructed by a stratigraphic inversion of seismic data, and a lithology cube identifying argillaceous lithologies and non-argillaceous lithologies is deduced therefrom. A relationship for estimating the fluid pressure from seismic P wave velocities is then determined from well data and for each one of the two lithologies. Finally, the fluid pressures in the underground medium are assessed by constructing a fluid pressure cube by applying the relationship to the seismic P wave velocity cube as a function of the lithology cube.

Abstract: A method for transferring data from an unmanned offshore energy platform to a manned offshore energy platform. The method includes obtaining the data on the unmanned offshore energy platform from a monitoring device, wherein the monitoring device is configured to monitor at least one environmental condition surrounding the unmanned offshore energy platform. The method further includes formatting the data for transmission. In addition, the method includes wirelessly transmitting the data from the unmanned offshore energy platform to the manned offshore energy platform.

Abstract: The present invention relates to a method and system for processing multi-dimensional signal data to determine frequency dependent features therefrom. The multi-dimensional signal data are transformed into space-frequency or time-space-frequency domain, providing second signal data. At predetermined locations of at least a portion of the one of space and time-space of the second signal data a dominant feature corresponding to a largest value of the second signal data is determined. This is followed by the determination of a wave-vector corresponding to the dominant feature at each of the predetermined locations. Finally, a dip map, a frequency map, and an amplitude map are generated using the wave-vectors.

Abstract: A method of predicting subsurface properties of a geologic formation includes acquiring seismic data for a subsurface region including the geologic formation of interest, computing seismic attributes of the measured seismic data over at least part of this geologic formation, determining internally consistent rock properties representative of the geologic formation, generating models of the same part of the geologic formation with these rock properties, computing synthetic seismic data from the models, computing the same attributes from these synthetic seismic data, and using Bayesian analysis to predict, from the probability of modeled attributes given the models, the probability of the actual subsurface properties given the measured attributes.

Abstract: A method of forming a geological model of a region of the earth includes obtaining seismic data relating to the region, the seismic data including seismic traveltime uncertainty. A seismic velocity model of the region may also be provided and includes velocity uncertainty. Ray tracing may be performed on the seismic data using the velocity model to determine the three dimensional positions of a plurality of points of the region. The three dimensional positional uncertainties (7-12) of at least some of the points can be calculated from the traveltime uncertainty, the velocity uncertainty, and uncertainty in the ray propagation direction. This can be combined with the positions determined by the ray tracing to form a geological model.

Abstract: A method and apparatus for estimating velocity in a subsurface region. Seismic data for a subsurface region may be received. One or more attributes for the seismic data may be calculated. A posterior distribution may be generated. The posterior distribution may represent one or more probabilities of one or more velocities for the attributes. A velocity with uncertainty may be determined for the subsurface region based on the posterior distribution. A pore pressure with uncertainty may be determined based on the velocity with uncertainty.

Abstract: Systems and methods create velocity models for a single well or for a set of wells. In one implementation, a system optimizes a time-depth relationship applied to data points from a single well to estimate coefficients for a linear-velocity-in-time function that models the data points. The system optimizes by reducing the influence of outlier data points, for example, by weighting each data point to decrease the influence of those far from the velocity function. The system also reduces the influence of top and bottom horizons of geological layers by applying data driven techniques that estimate the velocity function in a way that reduces dependence on the boundary conditions. The systems and methods can also create velocity models based on data from a set of wells, applying a well weights method to reduce the influence of outlier wells and thereby prevent wells with aberrant data from degrading a correct velocity model.

Abstract: Embodiments are directed to systems and methods (200, 300) that enable spatial variability of surface waves to be accounted for in dispersion correction in seismic data processing. This yields superior surface wave noise mitigation, with reduced likelihood of attenuating signal. Embodiments are operative with spatially inhomogeneous media.

Abstract: Embodiments use seismic processing methods that account for the spatial variability of surface wave velocities. Embodiments analyze surface wave properties by rapidly characterizing spatial variability of the surface waves in the seismic survey data (302). Filtering criteria are formed using the spatial variability of the surface waves (204). The filtering criteria can then be used to remove at least a portion of the surface waves from the seismic data (206, 319). The rapid characterization involves estimating a local group velocity of the surface waves by cross-correlation of the analytic signals (302).

Abstract: Method for processing seismic data. In one implementation, the method includes converting a common midpoint (CMP) gather of seismograms into one or more interferogram common midpoint (ICMP) gathers, generating a semblance spectrum for each ICMP gather, stacking the semblance spectrum from each ICMP gather to generate a combined semblance spectrum and deriving a normal moveout (NMO) velocity profile from the combined semblance spectrum.

Abstract: Methods and apparatus are provided for generating velocity models for pre-stack depth migration (“PSDM”). Seismic, gravity and electromagnetic joint inversion input data are generated based on observed seismic, gravity and electromagnetic data (e.g., magnetotelluric and/or controlled source electromagnetic), and velocity, density and resistivity models. A joint inversion is performed to produce a multiparametric model that is a function of velocity, density and resistivity parameter distributions. The separate parameter distributions are extracted from the multiparametric model, and the extracted velocity model is used to perform a PSDM. A migration velocity analysis is performed on PSDM output to generate seismic image residuals. If the seismic image residuals meet predetermined quality objectives, the extracted velocity model is output as the final velocity model for PSDM.

Abstract: Systems and methods in which data compression techniques are utilized to fill a predetermined channel capacity are shown. According to one configuration, event data points within test data are selected for communication via the data communication channel and a data decimator is utilized to identify other data points within the test data to fill or substantially fill the predetermined channel capacity. The foregoing data decimator may employ one or more variables for selecting data for communication, wherein one or more of the variables are preferably adjusted in decimator iterations to select an optimum or otherwise desirable subset of data for communication. Data decimators may additionally or alternatively implement a suitable “growth” function to select the particular data for communication and/or the amount of data communicated.

Abstract: A method for adjusting an isotropic depth image based on a mis-tie volume is provided. The method generally includes obtaining an isotropic velocity volume for a geophysical volume, obtaining an isotropic depth image of the geophysical volume, obtaining time-depth pairs at downhole locations in the geophysical volume, generating mis-tie values based on the time-depth pairs and the isotropic velocity volume, assigning uncertainties to the mis-tie values, generating a smoothest mis-tie volume that satisfies a target goodness of fit with the mis-tie values. Adjustment of the isotropic depth image may be achieved based on the mis-tie volume or a calibration velocity obtained from the mis-tie volume.

Abstract: A method of TI formation evaluation is disclosed. The method comprises receiving a plurality of borehole measurements; deriving a correlation between a first TI stiffness parameter and other TI stiffness parameters where the first and other TI stiffness parameters representing mechanical behavior of the TI formation; and computing the first and other TI stiffness parameters based on the borehole measurements and the derived correlation. The method further comprises evaluating TI formation elastic properties based on the computed first and other TI stiffness parameters. The method further comprises assuming that the shear modulus parallel to TI symmetric axis can be approximated from other moduli.

Abstract: A method of subsurface imaging is disclosed, the method involving using input data relating to the geometry of a fracture in a subsurface medium, mechanical properties of the medium and in-situ stress states in the medium, to model the behaviour of the fracture and obtain predicted information relating to deformation of the fracture. The predicted information is then used to constrain the interpretation of measured microseismic data caused by movement of the fracture to determine characteristics of the movement and/or failure of the fracture.

Abstract: The invention relates to methods and computer-readable medium to implement computing the propagation velocity of seismic waves in the earth. The invention computes the true propagation velocity of seismic waves in the earth, which is a condition of obtaining an accurate image of subsurface geology that can be used to prospect for oil and gas deposits. In an embodiment, the method of computing the propagation velocity of seismic waves in earth, includes providing an estimate of the propagation velocity, generating a time shift gather using a depth migration at a plurality of locations of the earth, converting each of the time shift gathers to a semblance gather, transforming each semblance gather into a velocity gather whose energy peaks represent a root-mean-square average of the propagation velocity along the forward and backward path between earth's surface and a point of the subsurface geology, and converting the energy peaks to the propagation velocity.

Abstract: A method for predicting a processing result includes a process of performing a principal component analysis on a plurality of detected data obtained during a first standard processing, to construct a principal component analysis model; a process of obtaining residuals of the principal component analysis model as first residuals; a process of performing a second standard processing; and a process of obtaining a plurality of detected data from a plurality of detectors during the second standard processing. And, the detected data obtained during the second standard processing are applied to the principal component analysis model, to obtain second residuals. The method further includes a process of weighting the second residuals based on weighting references, and constructing a new multivariate analysis model with use of the weighted second residuals; and a process of predicting a processing result of the second standard processing with use of the multivariate analysis model.

Abstract: Method for subsurface exploration. In one implementation, the method may include receiving a gamma value at a subsurface interface, wherein the gamma value is a ratio of Vp/Vs at the subsurface interface. Vp is the average P wave velocity and Vs is the average S wave velocity. The method may further include receiving a converted wave reflection coefficient at the subsurface interface. The converted wave reflection coefficient is a function of a density contrast at the subsurface interface and a shear velocity contrast at the subsurface interface. The method may further include determining a P wave incident angle such that the converted wave reflection coefficient is a function of substantially only the density contrast at the P wave incident angle. The P wave incident angle is referred to herein as ?. The method may further include calculating the density contrast based on the converted wave reflection coefficient at ?.

Abstract: A method for the quantitative characterization of a seabed sediment composition and a seabed's layered subbottom structure from at least one normal-incidence, single-channel reflection acoustic amplitude time series seismogram is disclosed. The method detects a plurality of reflections from subbottom interfaces in said seismograms, determines the traveltime, the polarity and the reflectivity of each detected reflection, determines the intrinsic attenuation of the sediment layer between pairs of adjacent reflections, and determines the acoustic properties, layer thicknesses and material properties of the seabed's layered subbottom structure as a function of said traveltimes, polarities and reflectivities of the detected reflections and said intrinsic attenuation of the sediment layer between pairs of adjacent reflections. A forward model describing the physical relationship between the material properties and the acoustic properties of seabed sediments is also disclosed.

Abstract: The invention relates to methods and computer-readable medium to implement computing the propagation velocity of seismic waves in the earth. The invention computes the true propagation velocity of seismic waves in the earth, which is a condition of obtaining an accurate image of subsurface geology that can be used to prospect for oil and gas deposits. In an embodiment, the method of computing the propagation velocity of seismic waves in earth, includes providing an estimate of the propagation velocity, generating a time shift gather using a depth migration at a plurality of locations of the earth, converting each of the time shift gathers to a semblance gather, transforming each semblance gather into a velocity gather whose energy peaks represent a root-mean-square average of the propagation velocity along the forward and backward path between earth's surface and a point of the subsurface geology, and converting the energy peaks to the propagation velocity.

Abstract: A visualization system for wellbore and drillstring data includes a graphics processor for creating a wire mesh model of a well and drillstring based on datasets of depth-varying parameters of the well. A graphics system maps appropriate textures to the wire mesh models, which are then displayed on a graphics display. A user interface facilitates user navigation along the length of the well to any selected location therein, and further permits user adjustment of orientation of the displayed renderings. The data is sufficient to permit calculation of fluid velocity in the wellbore at any selected location. The fluid velocity is presented as a velocity profile in the rendered visualization of the wellbore and drillstring to provide the user with a visual indication of fluid velocity in the wellbore as the user navigates the visualization along the length of the wellbore and drillstring.

Abstract: A method for updating a velocity model of a subsurface of the earth is described herein. A tomographic update to the velocity model of the subsurface may be performed to generate a tomographic velocity model update. A geomechanical velocity model update of the subsurface may be calculated. The geomechanical velocity model update may be combined with the tomographic velocity model update.

Abstract: Method of processing seismic data corresponding to acquisitions carried out in a medium exhibiting azimuthal anisotropy, characterized in that for two sets of data corresponding, for one and the same seismic variable and for one and the same ground zone, to at least two different source/azimuth receiver pair sectors, an average between the data which correspond to the azimuth considered and data which correspond, for the same variable and the same ground zone, to a perpendicular azimuth sector is determined for each of the azimuths considered, and in that an estimate of the part common to the two sets thus obtained of averaged data is determined.

Abstract: The present invention relates to a method and system for processing seismic data. Seismic data indicative of at least a time-dependent seismic trace are transformed into time-frequency domain using the S-transform. The transformed seismic data are processed for determining at predetermined time instances an instantaneous amplitude and a dominant instantaneous frequency. Based on a quotient of the instantaneous amplitude to the dominant instantaneous frequency seismic data indicative of a likelihood of a presence of hydrocarbons are then determined which are then graphically displayed. The method and system for processing seismic data provides a powerful tool for determining the likelihood of a presence of hydrocarbons based on a double peak hydrocarbon signature in the quotient of the instantaneous amplitude to the dominant instantaneous frequency.

Abstract: A visualization system for a wellbore environment includes a graphics processor for creating a wire mesh model of a well and optionally a drillstring based on datasets of depth-varying parameters of the well. A graphics system maps appropriate textures to the wire mesh models, which are then displayed on a graphics display. A user interface facilitates user navigation along the length of the well to any selected location therein, and further permits user adjustment of orientation of the displayed renderings. Simulated, real or a combination of simulated and real wellbore data may be visually depicted in the wellbore at any selected location. This provides the user with a visual indication of the wellbore environment as the user navigates the visualization along the length of the wellbore.

Abstract: A visualization system for a wellbore environment includes a graphics processor for creating a computer rendered visual model of a well, and optionally a drill string, based on data sets of depth-varying and/or time-varying parameters of the well. The model is then displayed on a graphics display. A user interface facilitates user navigation along the length of the well to any selected region therein, and further permits user adjustment of orientation of the displayed renderings as well as a temporal selection of the time-varying data to be displayed. Simulated, real, or a combination of simulated and real wellbore data, which may be steady state, transient, or real-time data, may be visually depicted at any selected region. This provides the user with a visual indication of the wellbore environment as the user navigates the visualization spatially and temporally.

Abstract: In general, in one aspect, the invention relates to a method for determining properties of a subterranean formation.

Abstract: Method for building a subsalt velocity model. In one implementation, the method may include determining a velocity of a sedimentary area surrounding a salt body, determining an initial velocity of a subsalt sedimentary area disposed below the salt body by interpolating the velocity of the surrounding sedimentary area with one or more boundaries defining the salt body, and updating the initial velocity of the subsalt sedimentary area.

Abstract: A method of determining in-situ a relation between the seismic velocity and the state of stress in an underground formation located under a surface subjected to time-changing surface loading conditions. A relation is selected between the seismic velocity and the state of stress containing at least one unknown parameter. A seismic source is arranged at surface or in a borehole penetrating the underground formation, and a seismic receiver is arranged at a distance from the seismic source at surface or in a second borehole. At two different times the seismic velocity of the formation along a path from the seismic source to the seismic receiver is determined. The difference in surface loading conditions at the two times is converted in a difference in states of stress in the underground formation. The unknown parameter(s) are calculated to obtain the relation between the seismic velocity and the state of stress in the underground formation.

Abstract: The invention relates to a method for determining maximum horizontal stress in an earth formation. The method includes obtaining fast shear wave velocities (Vs1) and slow shear wave velocities (Vs2) for various depths in the earth formation, calculating shear wave anisotropy (Adata) using Vs1 and Vs2, obtaining vertical stress (Sv) and minimum horizontal stress (Sh) for the formation, representing maximum horizontal stress (SH) using a parameterized function having at least one parameter and using Sh and Sv as input, determining a value of the at least one parameter by minimizing a cost function that represents a measure of difference between Adata and Apred for the various depths and Apred is predicted shear wave anisotropy determined using Sv, Sh, and SH, calculating SH using the parameterized function and the value of the at least one parameter, and storing SH in relation to the earth formation.

Abstract: A method for modeling geological structures includes identifying fault planes in a geologic structure, and creating a set of substantially vertically oriented pillars along the fault plane. The pillar grids are then combined into a common three-dimensional grid pillar network, and a two-dimensional plane is defined that intersects the middle node of the pillars. Planes are identified that are oriented in the substantially horizontal direction and horizontal node points are created, and such a process is repeated for other node levels. A skeleton grid is then created by drawing a substantially vertically oriented pillar through each corresponding grid intersection, and then primary horizons as defined by interpretation of seismic data are inserted into the skeleton grid. And, a final scale resolution is created by dividing the spaces bounded between actual horizons and fault planes, and eventual actual part of the outer boundary into one or more horizontally oriented fine layers.

Abstract: Method for conducting a seismic survey (23) wherein the acquisition parameters are modified as the survey progresses using information (24) from sensors located to monitor structures and fixtures located in the survey area. The sensor output is compared to acceptable levels, and acquisition parameters such as the seismic source strength are adjusted (26) to prevent exceeding the acceptable levels. An automated arbitration scheme (22) is used to resolve conflicting system goals within survey resource constraints.

Abstract: A method of processing geophysical data including determining the azimuth of a mirror symmetry plane within the earth from the geophysical data having sets of geophysical data acquired with different source-receiver azimuths. The data are processed to determine attributes of the geophysical data that are azimuth-dependent. One azimuth-dependent attribute of the geophysical data is selected, and a value of the source-receiver azimuth about which the selected attribute exhibits mirror symmetry is determined. This locates a mirror symmetry plane within the earth. The value of the source-receiver azimuth may be determined by generating an objective function indicative of the difference between the actual value of an attribute at one azimuth and the value predicted for that attribute at that azimuth using a trial estimate of the azimuth of a mirror symmetry plane, and finding the azimuth of a mirror symmetry plane that minimizes the objective function.

Abstract: A Laplace transform system comprising a processor, a measured time domain wavefield, a velocity model, and Laplace damping constants, wherein the processor is programmed to calculate a substantially about zero frequency component of a Fourier transform of a time domain damped wavefield, wherein the time domain damped wavefield is damped by the Laplace damping constants to obtain long wavelength velocity information for deeper subsurface regions.

Abstract: Data processing means are described for calculating an estimated time of arrival of a seismic or microseismic P or S wave at a sensor station, based on a P to S wave velocity ratio, a calculated estimated time of origin of the event and, where the estimated arrival time of a P wave is to be calculated, a picked arrival time of the S wave at the sensor station or, where the estimated arrival time of a S wave is to be calculated, a picked arrival time of the P wave at the sensor station.

Abstract: Method for processing seismic data. In one implementation, the method includes converting a common midpoint (CMP) gather of seismograms into one or more interferogram common midpoint (ICMP) gathers, generating a semblance spectrum for each ICMP gather, stacking the semblance spectrum from each ICMP gather to generate a combined semblance spectrum and deriving a normal moveout (NMO) velocity profile from the combined semblance spectrum.

Abstract: Method for processing seismic data to correct for errors caused by variation in water conditions.

Abstract: A method of and apparatus for processing seismic data. A method of processing seismic data comprising the steps acquiring seismic data in at least two shots that have different source-receiver azimuths from one another, and selecting events arising from interface seismic waves in each shot. The interface wave events are processed to provide information about azimuthal dependence of the velocity of seismic energy in the near-surface of the earth's interior. The invention may be used to process events arising from Scholte waves in marine seismic data, or to processing events arising from Love waves or Rayleigh waves in land seismic data. The method may be applied to seismic data acquired specifically to provide information about the near-surface or to seismic data acquired to provide information about the underlying geological structure.

Abstract: System and method for performing anisotropy corrections of post-imaging seismic data for a subsurface formation. The method may receive seismic data, preferably pre-stack seismic data comprising a plurality of traces, e.g., collected from a plurality of source and receiver locations. The seismic data may be imaged/migrated to produce imaged seismic data, which is organized into an arrangement that preserves aspects of the relative seismic propagation angle in the subsurface. One or more anisotropic parameters and corresponding corrections may be determined by analyzing the organized imaged seismic data. The determined parameters or corrections may be used to correct at least a subset of the imaged seismic data, thereby producing corrected seismic data which is useable in analyzing the formation. The corrected data may then optionally be stacked to produce a collection of corrected stacked traces, and/or analyzed as desired. The pre-stack data and/or the corrected seismic data may optionally be displayed.

Abstract: Method of determining the specular source/receiver pairs associated with seismic images produced from a transformation of multi-offset seismic wave records as a function of time in the depth domain. When using a prestack depth migration algorithm of Kirchhoff type, one also calculates and saves, in addition to the necessary traveltimes of the rays from the sources and receivers of the acquisition device to the image points of the image to be migrated, the line parameters of these rays and the desired quantities to be reached in fine. Then, once the prestack migrated volume produced, the local slopes of seismic events of interest for the image points considered are extracted. Comparison between line parameters saved during migration and the normal to the local slopes at the image points considered shows which ray is specular among those for which information has been saved, and allows to produce then the required quantity, which is then referred to as specular quantity.

Abstract: Methods and apparatus are provided for generating velocity models for pre-stack depth migration (“PSDM”). Seismic, gravity and electromagnetic joint inversion input data are generated based on observed seismic, gravity and electromagnetic data (e.g., magnetotelluric and/or controlled source electromagnetic), and velocity, density and resistivity models. A joint inversion is performed to produce a multiparametric model that is a function of velocity, density and resistivity parameter distributions. The separate parameter distributions are extracted from the multiparametric model, and the extracted velocity model is used to perform a PSDM. A migration velocity analysis is performed on PSDM output to generate seismic image residuals. If the seismic image residuals meet predetermined quality objectives, the extracted velocity model is output as the final velocity model for PSDM.

Abstract: A technique for calculating traveltime of a seismic wave in three dimensional tilted transversely isotropic (3D TI) media includes determining a wave vector, defining a unit vector, calculating an angle of the wave vector from an axis and performing a slowness determination, wherein a unit vector for a symmetry axis is defined as: (cos ? sin ?, sin ? sin ?, cos ?); where ? represents the azimuth of the symmetry axis measured from the x direction; and, ? represents the dip angle of the symmetry axis measured from the z direction. The technique may be practiced as a computer implemented set of instructions, and may be incorporated into measurement equipment.