Abstract: Seismic data recorded by subsurface seismic sensors placed in a borehole, such as an oil or gas well, are transformed via a process of upward wavefield propagation to pseudo-receivers at the surface of the earth. The seismic data thus transformed can be processed as though the data had been recorded by the pseudo-receivers at the surface rather than in the borehole where the data were actually recorded. This method accurately accounts for seismic source statics, anisotropy, and all velocity effects between the real receivers in the borehole and the pseudo-receivers at the surface of the earth.

Abstract: A borehole receiver for seismic testing has a support frame and an inflatable membrane configured to form a single, compact unit. The inflatable membrane has at least one sensor, such as an accelerometer affixed to the membrane. The at least one accelerometer may be adhesively affixed to an interior portion of the membrane with a cyanoacrylate adhesive. A pneumatic inflation system is used to inflate the membrane with pressurized air. In use, the membrane is slipped over the support frame and secured with a pair of O-rings. Once the borehole receiver is properly positioned in the borehole, pressurized air from the pneumatic inflation system inflates the membrane such that the at least one accelerometer achieves intimate coupling with a wall of the borehole to accurately measure dynamic seismic waves.

Abstract: The methods described herein are conceptually similar to classical migration/imaging for surface seismic data. However, instead of using computed (estimated) Green's functions in the imaging process, the instant invention utilizes measured (near-exact) Green's functions from VSP data to image the surface seismic data. Although the instant invention is best utilized where the velocity profile is approximately 1D (i.e., v(z)), the methods disclosed herein can also be extended to instances where there are some lateral velocity variations. Under these conditions, the instant invention allows for imaging surface seismic data and ‘self-imaging’ VSP data without first having to estimate a velocity model. The measurements obtained from the VSP data can also be used as a tool for calibrating computed Green's functions and migration operators.

Abstract: Channel estimation and signal equalization is used in a mud-pulse telemetry system for uplink communication during drilling of wellbores.

Abstract: A walkaway VSP survey is carried out with receivers located in a borehole near the base salt. Reflection tomographic inversion of data from the walkaway VSP is used to derive a velocity model for the subsurface and may be used for imaging of sub-salt reflections.

Abstract: A method of imaging a target area of the subsoil from “walkaway” data having application to development of oil reservoirs or monitoring of geologic storage sites is disclosed. After acquisition of seismic data of walkaway type and estimation of the rate of propagation of the seismic waves in the subsoil, p illumination angles are selected. The seismic measurements are then converted to data Dp by illumination angle. The distribution of acoustic impedances best explaining data Dp is determined within the target by using a non-linear inversion which minimizes a difference between the data Dp obtained from measurements and data by illumination angle resulting from an estimation. This estimation is performed by solving a wave propagation equation from the velocity field, an acoustic impedance distribution and a pressure distribution at the level of the top of the target for each illumination angle.

Abstract: Methods and apparatuses are disclosed for generating three dimensional (3D) vertical seismic profiles (VSPs) in a more time efficient manner. The methods and apparatuses enable faster and more efficient VSP surveys than previous techniques. Real-time updating of the velocity model and real-time reflection point density calculations are carried out, which are used to determine the location of the next seismic source firing. In the event the data of a particular bin has an insufficient fold, common image points (CIPs) or reflection point density, additional source firings may be carried out prior to moving the source. Further, in the event the data of a particular bin is excessive in terms of the fold, CIPs or reflection point density before the planned source firings for the bin are completed, remaining source firings for that bin may be skipped to save time and improve the efficiency of the data collection.

Abstract: A method is disclosed for more accurately determining the relative bearing angle of a directional receiver in a borehole using an existing three-dimensional (3D) geological model, one or more seismic sources and a three component (3C) directional receiver. A disclosed method includes: receiving direct compressional arrivals generated by multiple source events at the directional receiver disposed in the borehole; rotating the seismic data into the true earth frame using an estimated relative bearing angle; measuring a polarization vector of the rotated seismic data; estimating an incident ray vector of the direct compressional arrivals at the directional receiver using ray-tracing through the 3D model; calculating the weighted sum of an angular difference between the polarization vector and the incident ray vector; and adjusting the estimated relative bearing angle until the angular difference between the polarization and incident ray vectors is minimized.

Abstract: An apparatus and a method for migration of three components, 3-Dimensions seismic (3-C, 3-D) data acquired by down-hole receivers and surface seismic sources. This method utilizes full 3 components reflection wave field. It uses a dynamic, vector energy mapping method to image a reflection position and maps each time sample only to its reflected image point. Therefore, this method reduces unwanted data smearing and false mirror images. This method overcomes the weakness of using only a single component trace or pre-rotated three-component traces in the 1-C or 3-C 3-D VSP migration and produces better 3-D image.

Abstract: A walkaway VSP survey is carried out with receivers located in a borehole near the base salt. Reflection tomographic inversion of data from the walkaway VSP is used to derive a velocity model for the subsurface and may be used for imaging of sub-salt reflections.

Abstract: Seismic data recorded by subsurface seismic sensors placed in a borehole, such as an oil or gas well, are transformed via a process of upward wavefield propagation to pseudo-receivers at the surface of the earth. The seismic data thus transformed can be processed as though the data had been recorded by the pseudo-receivers at the surface rather than in the borehole where the data were actually recorded. This method accurately accounts for seismic source statics, anisotropy, and all velocity effects between the real receivers in the borehole and the pseudo-receivers at the surface of the earth.

Abstract: A walkaway VSP survey is carried out using a receiver array. Using a vertical VSP survey and arrival times of surface multiples on the walkaway VSP, vertical interval velocities and the anisotropy parameters ? and ? are estimated. This may then be used to process surface seismic data to do a prestack depth migration of surface seismic data and used for interpretation. For multi-azimuthal walkaway or 3D VSP data, we determine two VTI parameters ? and ? for multi-azimuth vertical planes. Then we determine five anisotropic interval parameters that describe P-wave kinematics for orthorhombic layers. These orthorhombic parameters may then be used to process surface seismic data to give a stacked image in true depth and for the interpretation purposes.

Abstract: A method includes providing a first data set representative of a first 3D VSP of a first region of a subterranean formation, and providing a second data set representative of a second 3D VSP of a second region of the subterranean formation. The first data set and the second data set define a common region of the subterranean formation. The first data set and the second data set are merged within the common region to produce a third data set representative of a third 3D VSP of the first and second regions of subterranean formation. The third data set is then stored on a computer readable medium.

Abstract: A method includes seismic wave field continuation, imaging and data analysis steps that are applied in a near well region.

Abstract: A method includes providing a first velocity model obtained from a VSP survey representative of an upper region of a subterranean formation. Wavefield equations from the first velocity model are datumed to a datum line between the upper region and a target area beneath the upper region to obtain datumed wavefield equations. The method further includes obtaining interferometric common shot data and interferometric common midpoint data from the datumed wavefield equations at the datum line. The first velocity model, the datumed wavefield equations, and the interferometric common midpoint data are then used to generate a second velocity model representative of velocities in the target area.

Abstract: A while-drilling Vertical Seismic Profiling (VSP) data acquisition system utilizing the same seismic shots for three purposes is disclosed. First, the seismic shots provide a means for synchronizing a downhole clock in the VSP receiver to a master clock at the surface, thereby enabling correct determination of seismic travel times. Second, the same seismic shots are also used to communicate commands and other information to the downhole VSP receiver, such commands controlling the actions of the VSP receiver or associated devices. Third, the same seismic shots are utilized for purposes of the VSP survey itself, i.e. determination of seismic travel times, forming of seismic images, and determination of geologic and formation fluid properties using the VSP methods.

Abstract: 1st order free-surface multiples recorded in VSP data or reverse VSP data are processed using a 3-C 3-D vector migration method to produce an image of the subsurface. This image produces a larger coverage than that obtained in 3-C 3-D processing of reflection data acquired in the VSP.

Abstract: Methods and apparatus for creating a velocity profile of a formation surrounding a borehole by checkshot measurements while moving the tool along the borehole. A conveyance and a sensor section are configured to move the sensor section in the borehole. At least one receiver is configured to detect signals generated at or near the surface while the sensor section is moving in the borehole.

Abstract: A method and system for acquiring seismic data while conducting drill string operations in a wellbore. The seismic receiver combination set comprises a combination of orthogonal geophones and accelerometers, and an array of hydrophones.

Abstract: Resonance scattering analysis of at least 3-component (3-C) VSP data detects heterogeneities in the proximity of a borehole. A method for seismic exploration of a pre-determined volume of the earth for assessing features of the volume using at least 3-C VSP data generated for the volume comprises: computing a resonance spectra indicating resonance scattering of the at least 3-C VSP data; and determining a lateral continuity of said features in accordance with the resonance spectra. Zero amplitude in a resonance spectrum indicates definite polarization of the direct pressure wave into the ray direction and very weak lateral heterogeneity along the path of the direct wave. High amplitudes in a resonance spectrum are observed if energy of the direct wave is observed on the horizontal components due to scattering at small-scale lateral heterogeneities near the receiver. Peak frequency provides information on the size and location of the scattering structure.

Abstract: A system and method for microseismic fracture mapping using seismic source timing measurements for velocity calibration is disclosed. The system may include a seismic source coupled to a wireline and a seismic source trigger, a sensor capable of detecting a first signal from the seismic source trigger, a transmitter coupled to the sensor, capable of transmitting a time value associated with the first signal, a receiver capable of detecting an event generated by the seismic source; and an analyzer capable of calculating a microseismic velocity of the event. In one embodiment, a first signal is transmitted through a wireline to trigger a seismic source. The first signal is detected, and a first time value associated with the first signal is transmitted. An event associated with the seismic source is detected and a second signal associated with the event is transmitted to an analyzer.

Abstract: A procedure to determine VSP-WD (Vertical Seismic Profiling—While Drilling) first breaks in absolute time is described. By combining a seismic surface source and firing complex seismic pattern unique in time (with a highly accurate surface—downhole time reference system), an automatic downhole procedure for first break detection on the downhole data, even under harsh conditions (low signal to noise ratio), is established. It is emphasized that this abstract is provided to comply with the rules requiring an abstract which will allow a searcher or other reader to quickly ascertain the subject matter of the technical disclosure. It is submitted with the understanding that it will not be used to interpret or limit the scope or meaning of the claims. 37 CFR 1.72(b).

Abstract: The invention is directed to a system for detecting seismic waves. The system has one or more sensor modules. Each sensor module has a detection unit, a positioning module, a digitizer, a radio transmitter, and a power supply. The system also includes a communications interface including a receiver, a data storage device, and a data relay module, and a data processor. The system may be used to detect seismic events by positioning sensor modules in an area, positioning a communications interface module in an area, establishing communication, polling the sensor modules for data, and relaying the data. The polling and relaying may be repeated at predetermined time intervals. Then, analysis may be performed on the data, and the seismic event may be identified as a precursor.

Abstract: The detailed analysis of cross well seismic data for a gas reservoir in Texas revealed two newly detected seismic wave effects, recorded approximately 2000 feet above the reservoir. A tube-wave (150) is initiated in a source well (110) by a source (111), travels in the source well (110), is coupled to a geological feature (140), propagates (151) through the geological feature (140), is coupled back to a tube-wave (152) at a receiver well (120), and is and received by receiver(s) (121) in either the same (110) or a different receiving well (120). The tube-wave has been shown to be extremely sensitive to changes in reservoir characteristics. Tube-waves appear to couple most effectively to reservoirs where the well casing is perforated, allowing direct fluid contact from the interior of a well case to the reservoir.

Abstract: An apparatus and a method for processing of three components, 3-Dimensions seismic (3-C, 3-D) data acquired by down-hole receivers and surface seismic sources. Automatic velocity analysis is used to identify the velocities of dominant events in a VSP panel. Different wave-types (downgoing P, downgoing PS, upcoming PS and upcoming PP) are identified and sequentially removed.

Abstract: A while-drilling Vertical Seismic Profiling (VSP) data acquisition system utilizing the same seismic shots for three purposes is disclosed. First, the seismic shots provide a means for synchronizing a downhole clock in the VSP receiver to a master clock at the surface, thereby enabling correct determination of seismic travel times. Second, the same seismic shots are also used to communicate commands and other information to the downhole VSP receiver, such commands controlling the actions of the VSP receiver or associated devices. Third, the same seismic shots are utilized for purposes of the VSP survey itself, i.e. determination of seismic travel times, forming of seismic images, and determination of geologic and formation fluid properties using the VSP methods.

Abstract: An apparatus and a method for migration of three components, 3-Dimensions seismic (3-C, 3-D) data acquired by down-hole receivers and surface seismic sources. This method utilizes full 3 components reflection wave field. It uses a dynamic, vector energy mapping method to image a reflection position and maps each time sample only to its reflected image point. Therefore, this method reduces unwanted data smearing and false mirror images. This method overcomes the weakness of using only a single component trace or pre-rotated three-component traces in the 1-C or 3-C 3-D VSP migration and produces better 3-D image.

Abstract: In-sea seismic energy sources, deployed at or near the seabed, generate seismic waves in earth formations beneath the seabed when activated by a triggering mechanism at a precise time synchronized with seismic data acquisition by sensors deployed at predetermined locations around a well drilled for oil and gas exploration. Seismic waves generated at the seabed travel with more direct and predictable ray-paths through the earth to the sensors. Acquired data, including a precise time and location of the seismic events, are used to provide enhanced borehole and sea-bottom seismic surveys.

Abstract: 1st order free-surface multiples recorded in VSP data or reverse VSP data are processed using a 3-C 3-D vector migration method to produce an image of the subsurface. This image produces a larger coverage than that obtained in 3-C 3-D processing of reflection data acquired in the VSP.

Abstract: The detailed analysis of cross well seismic data for a gas reservoir in Texas revealed two newly detected seismic wave effects, recorded approximately 2000 feet above the reservoir. A tube-wave (150) is initiated in a source well (110) by a source (111), travels in the source well (110), is coupled to a geological feature (140), propagates (151) through the geological feature (140), is coupled back to a tube-wave (152) at a receiver well (120), and is and received by receiver(s) (121) in either the same (110) or a different receiving well (120). The tube-wave has been shown to be extremely sensitive to changes in reservoir characteristics. Tube-waves appear to couple most effectively to reservoirs where the well casing is perforated, allowing direct fluid contact from the interior of a well case to the reservoir.

Abstract: A method of acquiring seismic data that includes deploying a first array of seismic receivers and a second array of seismic receivers, and simultaneously receiving drill noise seismic energy produced as a wellbore is drilled relatively near the first and second arrays of seismic receivers using the first and second array of seismic receivers, wherein the first array of seismic receivers is closer to the source of the drill noise seismic energy than the second array of seismic receivers. A related method for processing seismic data and computer useable media are also described. Also a method of acquiring seismic data that includes deploying a first array of seismic receivers within a borehole, receiving seismic energy produced as a wellbore is drilled relatively near first array using said seismic receivers, and recording seismic data associated with the received seismic energy. A related method for acquiring and processing seismic data is also described.

Abstract: 1storder free-surface multiples recorded in VSP data or reverse VSP data are processed using a 3-C 3-D vector migration method to produce an image of the subsurface. This image produces a larger coverage than that obtained in 3-C 3-D processing of reflection data acquired in the VSP.

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 multi-scale finite-volume (MSFV) method to solve elliptic problems with a plurality of spatial scales arising from single or multi-phase flows in porous media is provided. The method efficiently captures the effects of small scales on a coarse grid, is conservative, and treats tensor permeabilities correctly. The underlying idea is to construct transmissibilities that capture the local properties of a differential operator. This leads to a multi-point discretization scheme for a finite-volume solution algorithm. Transmissibilities for the MSFV method are preferably constructed only once as a preprocessing step and can be computed locally.

Abstract: A method and apparatus for enhancing the moveout between a direct wave and a reflected wave. The method involves transmitting imaging signals into a body to be imaged and receiving the resulting signals propagated from the signal source. The step of receiving the propagated signals includes selectively adjusting the distance between the signal source and the signal receivers between successive signals. The method further comprises separating the reflected signals from the total received signals and enhanced stacking of the measured reflected signals.

Abstract: A while-drilling Vertical Seismic Profiling (VSP) data acquisition system utilizing the same seismic shots for three purposes is disclosed. First, the seismic shots provide a means for synchronizing a downhole clock in the VSP receiver to a master clock at the surface, thereby enabling correct determination of seismic travel times. Second, the same seismic shots are also used to communicate commands and other information to the downhole VSP receiver, such commands controlling the actions of the VSP receiver or associated devices. Third, the same seismic shots are utilized for purposes of the VSP survey itself, i.e. determination of seismic travel times, forming of seismic images, and determination of geologic and formation fluid properties using the VSP methods.

Abstract: A seismic surveying arrangement comprises a vertical source array (26) having at least two sources (27,28) of seismic energy. The source array (26) is disposed such that a first source (28) is located within a weathered layer (35) of the earth's interior and such that a second source (27) is located within the weathered layer at a greater depth than the first source or within a high velocity layer (8) underlying the weathered layer (35). The surveying arrangement also includes a first horizontal receiver array (29) located in the weathered layer, and means for processing seismic data acquired by the first receiver array. The processing means are able to process the acquired seismic data so as to attenuate seismic events involving reflection of seismic energy at at least one boundary face of the high velocity layer. The seismic surveying arrangement may also comprise a second receiver array (30) disposed at a lower depth than the first receiver array (29).

Abstract: A method and system for seismic data processing utilizes azimuthal variations in the velocity of seismic signals. The system and method utilizes a plurality of seismic energy sources that are located at known positions at the surface of the earth. The seismic energy sources generate seismic signals that propagate downward into the earth. Some of the seismic signals are reflected and diffracted by various sub-surface layers and are returned to the surface of the earth. The returned seismic signals are received by a plurality of receivers. The method includes the step of determining the distance from an energy source to an image point. A fast travel time of the seismic signal from the energy source to the image point is determined, and a slow travel time of the seismic signal from the energy source to the image point is determined. The azimuth angle between the energy source and the surface location of the image point is calculated.

Abstract: Methods and apparatus for creating a velocity profile of a formation surrounding a borehole by checkshot measurements while moving the tool along the borehole. A conveyance and a sensor section are configured to move the sensor section in the borehole. At least one receiver is configured to detect signals generated at or near the surface while the sensor section is moving in the borehole.

Abstract: A method for improving the seismic resolution in the field of petroleum exploration and development, using the steps of: (1) ensuring the depth and the velocity of the low-velocity layer and the high-velocity layer and the characteristic of the frequency spectrum of the seismic wave at the near-surface using the borehole seismic survey; (2) separating the ghost wave from the seismic record in the high-velocity layer using the processing method of the vertical seismic profiling, and improving the SNR of the seismic record of the direct wave at the high-velocity layer using the uphole stacking of the vertical seismic profiling; (3) evaluating the de-convolution operator of the near-surface by the convolution math-physics equation; (4) evaluating the seismic record with high resolution by the de-convolution math-physics equation; (5) providing the band-pass filtering and geologic interpretation to the seismic record with high resolution.

Abstract: A technique includes providing a source in a first well and a seismic receiver in a second well to detect a seismic event that is caused by the source. The technique includes referencing clocks in the source and receiver to a common reference time frame and determining a time in the reference time frame at which the seismic source generates the seismic event.

Abstract: The methods described herein are conceptually similar to classical migration/imaging for surface seismic data. However, instead of using computed (estimated) Green's functions in the imaging process, the instant invention utilizes measured (near-exact) Green's functions from VSP data to image the surface seismic data. Although the instant invention is best utilized where the velocity profile is approximately 1D (i.e., v(z)), the methods disclosed herein can also be extended to instances where there are some lateral velocity variations. Under these conditions, the instant invention allows for imaging surface seismic data and ‘self-imaging’ VSP data without first having to estimate a velocity model. The measurements obtained from the VSP data can also be used as a tool for calibrating computed Green's functions and migration operators.

Abstract: Implementations of various techniques for processing seismic data acquired by three or more over/under streamers. In one implementation, the seismic data may be processed by grouping the seismic data into one or more sets of seismic data based on all possible pairs of the three or more over/under streamers, wherein each set of seismic data comprises seismic data that had been acquired by a pair of the three or more over/under streamers, applying a shift and subtract algorithm to the one or more sets of seismic data, and applying a dephase and sum algorithm to the one or more sets of seismic data.

Abstract: Seismic data recorded by subsurface seismic sensors placed in a borehole, such as an oil or gas well, are transformed via a process of upward wavefield propagation to pseudo-receivers at the surface of the earth. The seismic data thus transformed can be processed as though the data had been recorded by the pseudo-receivers at the surface rather than in the borehole where the data were actually recorded. This method accurately accounts for seismic source statics, anisotropy, and all velocity effects between the real receivers in the borehole and the pseudo-receivers at the surface of the earth.

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.

Abstract: Methods for determining the existence and characteristics of a gradational pressurized zone within a subterranean formation are disclosed. One embodiment involves employing an attenuation relationship between a seismic response signal and increasing wavelet wavelength, which relationship may be used to detect a gradational pressurized zone and/or determine characteristics thereof. In another embodiment, a method for analyzing data contained within a response signal for signal characteristics that may change in relation to the distance between an input signal source and the gradational pressurized zone is disclosed. In a further embodiment, the relationship between response signal wavelet frequency and comparative amplitude may be used to estimate an optimal wavelet wavelength or range of wavelengths used for data processing or input signal selection. Systems for seismic exploration and data analysis for practicing the above-mentioned method embodiments are also disclosed.

Abstract: Methods for determining the existence and characteristics of a gradational pressurized zone within a subterranean formation are disclosed. One embodiment involves employing an attenuation relationship between a seismic response signal and increasing wavelet wavelength, which relationship may be used to detect a gradational pressurized zone and/or determine characteristics thereof. In another embodiment, a method for analyzing data contained within a response signal for signal characteristics that may change in relation to the distance between an input signal source and the gradational pressurized zone is disclosed. In a further embodiment, the relationship between response signal wavelet frequency and comparative amplitude may be used to estimate an optimal wavelet wavelength or range of wavelengths used for data processing or input signal selection. Systems for seismic exploration and data analysis for practicing the above-mentioned method embodiments are also disclosed.

Abstract: A method of obtaining a spatial model of a property of part of a subsurface formation located between underground seismic receivers in which at least two sets of pairs of seismic receivers are utilized and one pair of receivers is used to record a signal from a seimic source and obtaining a response by solving (s11(?t){circle around (x)}s21(t))=r11,21(t){circle around (x)}(s11(?t){circle around (x)}s21(t)), wherein the symbol {circle around (x)} denotes convolution and wherein s11(?t) is the time-reverse of the signal s11(t). A path-related attribute is selected from transmission response r11,21(t) that corresponds to the property of the subsurface formation and a tomographic reconstruction technique is applied to the path-related attribute to obtain the spatial model of the property of part of the subsurface formation.

Abstract: A method is described for acquiring seismic data while removing, or tripping, a drill string from a wellbore. The method comprises, conveying an autonomous seismic receiver down a fluid passage in the drill string to a predetermined location proximate a lower end of the drill string; generating signals by a seismic source near a surface location; detecting the seismic signals with at least one sensor in the seismic receiver at least one predetermined location of interest in the wellbore as the drill string is removed from the wellbore; and storing the detected seismic signals in the seismic receiver.

Abstract: A method of applying an effective velocity model to vertical seismic profile (VSP) seismic data comprises correcting for offset using a non-hyperbolic effective velocity model so as to take account of the earth's layering and anisotropy. One preferred non-hyperbolic model for the relationship between offset and travel time is: formula (I) where t is the travel time of seismic energy from the source to the receiver, x is the offset between the source and the receiver, and z is the depth of the receiver.