Abstract: Method for identifying one or more fracture clusters in a formation surrounding a reservoir. In one implementation, the method may include generating a P to S image, comparing the P to S image to one or more images from a borehole, and identifying one or more fracture clusters using the P to S image and the borehole images.

Abstract: Cross-dipole measurements are obtained in a borehole. By estimating a direction of polarization of the fast shear mode at low and high frequencies and comparing the estimated distances, a cause of anisotropy is established. Formation stresses and directions may be estimated.

Abstract: A method for estimating principal stresses of a subterranean formation from seismic data. In one embodiment, rock strength parameters from seismic data of the formation is first determined to calculate the anisotropic elastic properties of the formation. The three principal stresses of the formation: vertical stress, minimum horizontal stress, and maximum horizontal stress, is determined using at least the calculated anisotropic elastic properties and the rock strength parameters of the formation. From the estimated principal stresses, the differential ratio of the maximum and minimum horizontal stresses can be determined to indicate optimal zones for hydraulic fracturing. In another embodiment, a tectonic strain term is introduced to calibrate the estimated principal estimated stress to a known reference point. In yet another embodiment, hoop stress is incorporated to estimate the fracture initiation pressures.

Abstract: A technique includes a method and apparatus for simultaneous joint inversion of surface wave and refraction data to identify near surface geophysical and geological properties.

Abstract: There is described a method of moveout or velocity analysis of seismic signals using the steps of obtaining such signals 5 from a plurality of receivers, identifying receiver functions within the acoustic signals, analyzing said receiver functions for velocity or moveout characteristics, using the result of said analyzing step to determine, properties of multiple layers of earth located below said 10 receivers. The analyses can involve the use of representation of the traveltime differences as approximated power series of slowness or horizontal distances. The method is the first to comprehensively deal with a multi-layered earth or velocity model.

Abstract: A microseismic method of monitoring fracturing operation or other passive seismic events in hydrocarbon wells is described using the steps of obtaining multi-component s-wave signals of the event; and using a linear derivative of S-wave arrival times of the signals in a first direction, an S-wave velocity and an s-wave polarization to determine at least two components of the S-wave slowness vector.

Abstract: A measurement method for a granular compaction pile 10 using a crosshole seismic test. The measurement method includes forming first to fifth measurement holes around the granular compaction pile, injecting grouting material into the first to fifth measurement holes and inserting casings into the first to fifth measurement holes, respectively, and installing an oscillator and a detector of a crosshole seismic tester in the first to fifth holes and detecting a velocity of a shearing wave, thereby measuring a diameter of the granular compaction pile according to a construction depth thereof. The shape and the diameter of the granular compaction pile are easily checked without causing damage to the granular compaction pile.

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: Method for identifying one or more fracture clusters in a formation surrounding a reservoir. In one implementation, the method may include generating a P to S image, comparing the P to S image to one or more images from a borehole, and identifying one or more fracture clusters using the P to S image and the borehole images.

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 p

Abstract: A method to generate images of acoustic contrasts for structures located between at least one acoustic source and at least one receiver, said structures converting a part of the compressional energy to shear.

Abstract: A method for generating an average gamma factor for recorded PS data and recorded PP data. The method includes mapping one or more arrival times of the recorded PS data onto one or more expected arrival times of the recorded PP data using an estimated average gamma factor, identifying one or more corresponding pairs of events in the recorded PS data and the recorded PP data using forward modeling, applying geophysical reasoning on the corresponding pairs of events and generating the average gamma factor.

Abstract: Maximum and minimum horizontal stresses, and horizontal to overburden stress ratio, are estimated using radial profiles of shear moduli. Inversion enables estimation of maximum and minimum horizontal stresses using radial profiles of three shear moduli associated with an orthogonal set of axes defined by the three principal stress directions. Differences in the far-field shear moduli are inverted together with two difference equations obtained from the radial profiles of the dipole shear moduli C44 and C55, and borehole stresses in the near-wellbore region. The horizontal to overburden stress ratio is estimated using differences in the compressional, dipole shear, and Stoneley shear slownesses at two depths in the same lithology interval where the formation exhibits azimuthal isotropy in cross-dipole dispersions, implying that horizontal stresses are nearly the same at all azimuths.

Abstract: A method of determining the nature of a submarine or subterranean reservoir having an electromagnetic characteristic is described. The method comprises analyzing data associated with the application of a time varying electromagnetic field from above the reservoir using a data analysis methodology that does not utilize an estimate of the reservoir's geometry in determining one or more characteristics of the reservoir.

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: Methods and apparatus facilitating measurement of anisotropy are disclosed. According to some aspects, anisotropy can be determined while drilling. Monopole/quadrupole interactions, as well as monitoring dipole excitations and other methods, may be used to find the principal shear directions of a formation while drilling or during wireline or other operations.

Abstract: A wave generating device and method of using the same wherein the device is for producing a shear wave in a ground layer wherein the generated wave emanates from a bore hole in the ground layer and passes through the ground layer toward a receiver spaced from the bore hole. This wave having a given waveform and the bore hole being defined by at least one bore wall that extends from a surface of the ground layer to a bottom extent of the bore hole. The device including a driver having an outer layer extending in a longitudinal direction along a longitudinal axis and having a radially outwardly facing surface generally parallel to the longitudinal axis which is configured to engage an associated bore wall of an associated bore hole.

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

Abstract: A method for identifying a local stress anomaly in a subsurface region, which method comprises the steps of obtaining a model of the subsurface region, which model includes a salt layer in between adjacent layers; identifying a salt weld in the model; and attributing a stress anomaly to an area surrounding the salt weld.

Abstract: Methods and apparatus facilitating radial profiling of shear slowness are disclosed. According to some aspects of the invention, acoustic tool bias is accounted for in the calculation of radial profiles. According so some aspects, acoustic tool bias is accounted for by replacing acoustic tool structure with a resonance-impedance model. The resonance-impedance modeling according to principles of the present invention is applicable to vertical, deviated, and horizontal boreholes.

Abstract: A method of monitoring heavy oil recovery in a hydrocarbon reservoir including the steps of obtaining a set of acoustic signals; deriving from the obtained signals a velocity model including shear and compressional velocities for at least a section of the reservoir; defining a relation relating shear and compressional velocities to at least one material parameter of the heavy oil; performing a reservoir treatment process; obtaining an in-situ or post-treatment set of acoustic signals; and deriving from the in-situ or post-treatment set of acoustic signals shear and compressional velocities and the at least one material parameter using the relation.

Abstract: A P-wave slowness-polarization data set resolving P-wave slowness against P-wave polarization angle is extracted from a multi-offset vertical seismic profiling data set. A function relating P-wave polarization angle to P-wave slowness, and being dependent on at least three free parameters, wherein in at least two of these free parameters information of shear wave velocity is absorbed together with one or more anisotropy parameters that characterize formation anisotropy is selected, and regression values found for the at least three free parameters by fitting the selected function to the P-wave slowness-polarization data set. The regression values may be utilized in a method of producing a mineral hydrocarbon fluid from an earth formation.

Abstract: A technique includes determining at least one attribute of a slowness vector associated with a seismic gather based on pressure data and an indication of particle motion that is measured by at least one seismic sensor while in tow.

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: A method of seismic processing for obtaining information on the geophysics of the subsoil includes acquiring seismic traces at at least one point on the surface of the subsoil or in the subsoil, the seismic traces corresponding on each occasion to two perpendicular components of a shear wave emitted into the subsoil and reflected by different interfaces therein; applying a succession of transformations (?) at least to a temporal portion of the traces for each of these assumptions, determining the value of a parameter representative of the coherence/similarity between the result traces obtained in this way; and election as a function of the values obtained in this way that one of the hypotheses which is considered as being most representative of the subsoil, the two result traces obtained for said hypothesis being compensated traces of the subsoil birefringency.

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: The present invention provides methods and systems for microseismic hydraulic fracture monitoring in real-time. The methods and systems of the present invention may include continuous map migration of recorded microseismic signals. The methods and systems provide robust automated simultaneous detection and location of microseismic events.

Abstract: A downhole tool for gathering formation data from inside a borehole includes a substantially tubular housing adapted for axial connection to a drill string and multiple sensors coupled to the tubular housing. The sensors include a first pair of sensors aligned along a first axis and adapted to measure a spatial derivative along the first axis, a second pair of sensors aligned along a second axis and adapted to measure a spatial derivative along the second axis, and a third pair of sensors aligned along a third axis and adapted to measure a spatial derivative along the third axis. In selected embodiments, the spatial derivatives are used to differentiate seismic or sonic compression and shear waves measured in a downhole environment.

Abstract: A method of determining seismic properties of a layer of the seabed, in particular a surface or near-surface layer (5), comprises directing seismic energy propagating in a first mode at a boundary face of the layer so as to cause partial mode conversion of the seismic energy at the boundary face. For example partial mode conversion may occur when seismic energy propagates upwards through the interface between a surface or near-surface layer (5) and the basement (6), owing to the difference in seismic properties between the surface or near-surface layer (5) and the basement (6). In the invention, the two modes of seismic energy—that is the initial mode and the mode generated by mode conversions at the interface—are received at a receiver. The difference in travel time of the two modes between the interface and the receiver is determined from the seismic data acquired by the receiver.

Abstract: A method is disclosed for seismic imaging of subsurface diffractors. The method includes performing migration velocity analysis on a seismic time record section and depth migrating the time section for offsets exceeding one-half a distance between a seismic energy source and a seismic receiver most distant from the source during acquisition of seismic data used to generate the time record section.

Abstract: A method of seismic data processing is described in which a particular wavelet is selected from a plurality of wavelets as being most characteristic of a received seismic signal. A subtracted signal can be determined by subtracting a weighted signal of the particular wavelet from the received signal. From the subtracted signal, an additional particular wavelet can be chosen. The process of subtracting a signal and determining an additional particular one of the plurality of wavelets can be repeated until a criterion is met. The method can be repeated at several depths. The resultant spectral analysis can be used to determine, for example, an absorption coefficient.

Abstract: The logging tool of this invention includes a transmitter conveyed on a drilling collar for exciting a quadrupole signal in a borehole being drilled by a drill bit and a receiver for receiving the signal. The transmitter is operated at a frequency below the cut-off frequency of the quadrupole collar mode. The received signal consists primarily of the formation quadrupole mode which, at low frequencies, has a velocity that approaches the formation shear velocity. The transmitter, in one embodiment, consists of eight equal sectors of a piezoelectric cylinder mounted on the rim of the drilling collar. The value of the cut-off frequency is primarily dependent on the thickness of the drilling collar. Alternatively, the transmitter may be operated to produce both the collar mode and the formation mode and a processor may be used to filter out the collar mode.

Abstract: A method of determining a fluid property in a subsurface region of interest of an earth formation uses measurements of seismic attributes on seismic data. For a test region, a plurality of realizations of rock properties are specified, and for each of the realizations and a selected value of a fluid property, the seismic attribute is modeled. This defines a probability density function (PDF). Comparison of the PDF of the model output with the PDF on the measured seismic data is used to determine the likelihood of the selected fluid property.

Abstract: A method of diagnostics of oil-and-gas deposit for seismic survey on the basis of the complex of waves of the different types. It is performed by means of the insonification of the vicinity of the borehole by the seismic waves, the recording of the displacement vectors of the direct longtidunial and shear waves using the three-component seismic profiling in the investigated borehole, the data processing of the detailed profiling for the obtaining of the dynamical, rheological and absorbing parameters of the rock, and the geological interpretation of the above-mentioned parameters for the obtaining of the quantitative estimates of the reservoir and fluid-saturated properties of the rocks for each observation point lengthways depth of the borehole.

Abstract: Array quadrupole data obtained in a LWD environment are processed to determine the shear velocity of an earth formation taking into account dispersion effects caused by the logging tool. When this processing is done in an azimuthally anisotropic formation, the shear velocity that is obtained is the slow shear velocity. When 4C (cross-dipole) data are also obtained, then both the fast and slow shear velocities of an azimuthally anisotropic formation can be obtained.

Abstract: An acoustic logging apparatus comprises a drill collar conveyed on a drilling tubular in a borehole within a formation. At least one transmitter is disposed in the drill collar. The transmitter includes at least one magnetostrictive actuator cooperatively coupled by a flexure ring to a piston for converting a magnetostrictive actuator displacement into a related piston displacement for transmitting an acoustic signal in the formation.

Abstract: The logging tool of this invention includes a transmitter conveyed on a drilling collar for exciting a quadrupole signal in a borehole being drilled by a drill bit and a receiver for receiving the signal. The transmitter is operated at a frequency below the cut-off frequency of the quadrupole collar mode. The received signal consists primarily of the formation quadrupole mode which, at low frequencies, has a velocity that approaches the formation shear velocity. The transmitter, in one embodiment, consists of eight equal sectors of a piezoelectric cylinder mounted on the rim of the drilling collar. The value of the cut-off frequency is primarily dependent on the thickness of the drilling collar. Alternatively, the transmitter may be operated to produce both the collar mode and the formation mode and a processor may be used to filter out the collar mode.

Abstract: A method of determining the polarisation directions of the fast and slow shear waves arising from shear wave splitting due to anisotropy, said directions defining a natural coordinate system, the method comprising the steps of: a) recording at least two components of each shear wave, in a recording coordinate system, b) calculating the value of &thgr;, being the angle of rotation between the natural coordinate system and the recording coordinate system, for which the Lp norm of the rotated traces is minimised if p is less than 2, or maximised if p is greater than 2.

Abstract: The patent discloses a signal processing technique for determining the fast and slow shear wave polarizations, and their orientation, for acoustic waves in an anisotropic earth formation. The signal processing method decomposes composite received waveforms a number of times using a number of different strike angles. The decomposed signals are used to create estimated source signals. The estimated source signals are compared in some way to obtain an objective function. Locations in a plot where the objective function reaches minimum values are indicative of the acoustic velocity of the fast and slow polarizations within the formation.

Abstract: A seismic prospection method in which a compression seismic wave is emitted into the subsoil and sensors are used to collect seismic data having at least a shear component, and in which the data corresponding to said shear component is processed to deduce information about the geology of the subsoil, the method being characterized in that an estimate of the ratio: ∫ Z 0 Z ⁢ v p ·   ⁢ ⅆ l / ∫ Z 0 Z ⁢ v s ·   ⁢ ⅆ l is determined where vp and vs are values for real loca

Abstract: A method for determining formation shear slowness of formations surrounding an earth borehole that includes: transmitting sonic energy into the formations from a dipole source in the borehole; measuring signals from orthogonal wave components that have traveled through the formations at each of a plurality of receiver locations in the borehole, spaced at a respective plurality of distances from the transmitter location; transforming the signals to the frequency domain and separating dispersive and non-dispersive portions of the transformed signals; and determining fast and slow shear slowness of the formation from the low frequency asymptotes of the dispersive portions of the transformed signals. Related methods for determining shear directionality and the magnitude of anisotropy of anisotropic formations surrounding an earth borehole and associated apparatus are also described.

Abstract: A method of processing seismic traces corresponding to two modes of propagation in the subsoil, the method being characterized by the following steps: determining at least one pair of correlated events from the seismic traces corresponding to the two modes respectively; determining the ratio of the differences between the vertical travel times Ts and Tp in the S domain and in the P domain, corresponding to said correlated pairs of events; and determining the values of Ts and Tp over the time interval defined by the two pairs of correlated events so that throughout said time interval, the ratio rt=Ts/Tp is equal to the ratio Vp/Vs of the speeds in the P domain and in the S domain, and so that the mean value of said ratio r over said time interval corresponds to the ratio of the differences between the vertical travel times corresponding to the pairs of correlated events.

Abstract: The present invention provides a system for, and method of determining anisotropy in a stratum using scattered horizontal shear and vertical shear modes. The method includes transforming seismic energy wave generated by a seismic energy source and received by seismic receivers into radial/transverse coordinate space, determining a seismic energy wave corridor along a radial path between the seismic energy source and the seismic receivers, and gathering seismic data received by the seismic receivers within the corridor wherein the data includes horizontal and vertical shear components.

Abstract: Controlled-polarization marine seismic shear-wave source has a base suitable for positioning on a seafloor, beneath a water body. Twin water-blast tubes (seismic-impulse-generation thrusters) are mounted on the base, aimed upwardly, outwardly in opposite directions at equal acute angles &bgr; relative to the base. These blast tubes are positioned symmetrically on opposite sides of vertical centerline CL of the source. When submerged, water enters open mouths at upper ends of the tubes. The base has multiple downwardly projecting blade-like teeth for penetrating into the seafloor. Firing a compressed-air discharger mounted in a lower end of a water-blast tube hurls a powerful slug of water upwardly along the tube, producing a powerful reaction impulse acting downwardly generally parallel with the tube axis.

Abstract: A seismic shear wave generator comprising a seismic wave source (3) arranged on a base plate (11) arranged for being placed at the seabed, and a coupling device (1) arranged for coupling the wave source (3) mechanically to the seabed in order for shear waves to be transmitted to the seabed. The novelty of the invention is as follows: a cylindrical skirt (1) closely arranged along the base plate's (11) circumference and extending down from the underside of the base plate (11), and open downwards, and a pump device (2) arranged for pumping water respectably out of, and into, the skirt (1) and form vacuum or overpressure, respectively, in a closed cavity formed by the skirt (1) and the base plate (11) and the seabed (1) is lowered onto the seabed. By means of the device mentioned above the shear wave generator may be suctioned fixed to the sea bed and thus transfer shear waves via the skirt (1) to the seabed.

Abstract: A system for, and method of, orienting seismic energy sources and seismic energy receivers to yield discriminated horizontal and vertical shear modes. In one embodiment, the system includes: (1) an initial matrix representing wave energy produced by a selected seismic energy source and received by a selected seismic energy receiver, (2) a first transforming matrix that cooperates with the initial matrix to orient the selected seismic energy source with respect to the selected seismic energy receiver and (3) a second transforming matrix that cooperates with the initial matrix to orient the selected seismic energy receiver with respect to the selected seismic energy source, the initial matrix thereby transformed into a successive matrix containing the discriminated horizontal and vertical shear modes.

Abstract: Streamlined, towable, marine seismic energy vibrator for creating intense swept-frequency and pulse-coded seismic signals in a body of water has a sleek, fish-like configuration designed for towing with minimum drag. The vibrator has a streamlined hollow towing head and a streamlined hollow tail head mounted onto front and rear of a long cylindrical tubular wall which is modular, comprising cylinder sections joined in end-to-end axial alignment. Within this long tubular cylinder wall is an axially vibratable multi-piston assembly having a plurality of pistons on a long piston rod. One piston is positioned in each of the cylinder chambers. These chambers hive multiple ports opening out through the long cylindrical wall. An elongated circular cylindrical elastomeric bladder forms a water-filled bladder chamber encircling the wall.

Abstract: The present invention includes a method for overpressure detection and pore pressure change monitoring in subsurface gas, liquid hydrocarbon, or water reservoirs from compressional- and shear-wave measurement data. As part of this method, one or more Poisson's ratios are determined from field-based measurement data and are then compared against known Poisson's ratio values representative of the particular subsurface formation type. By applying a Poisson's ratio—pore pressure criterion that is appropriate for that type of formation, an overpressure in the formation is identified.

Abstract: A method for the study of elastic wave phenomena in the ground, where the detection of said phenomena is carried out with the use of at least a geophone comprising at least a couple of transducers (20, 21; 37, 38), each transducer (20, 21; 37, 38) being apt to detect horizontally polarized shear waves, represented by SH signals, and compression waves, represented by P signals.

Abstract: The method of the invention is employed for determining the state of saturation in a subterranean formation using only seismic velocity measurements (e.g., shear and compressional wave velocity data). Seismic velocity data collected from a region of the formation of like solid material properties can provide relatively accurate partial saturation data derived from a well-defined triangle plotted in a (&rgr;/&mgr;, &lgr;/&mgr;)-plane. When the seismic velocity data are collected over a large region of a formation having both like and unlike materials, the method first distinguishes the like materials by initially plotting the seismic velocity data in a (&rgr;/&lgr;, &mgr;/&lgr;)-plane to determine regions of the formation having like solid material properties and porosity.