Methods for measuring seismic velocities and for monitoring local changes in inter-well seismic velocities in real time are described. Two or more spaced-apart observation wells are provided. Seismic receiver arrays are placed in the observation wells, and a seismic source array is provided at surface locations away from the well bores and producing areas. Compression (P), vertical shear (Sv) and/or horizontal shear (Sh) seismic wave signals are generated from each element of the seismic source array, and the seismic signals arriving at the receivers in the observation wells are recorded. The virtual source method is then applied to the recorded data to compute emulated cross-well seismic signals of the virtual sources at receiver locations in one observation well propagating toward the receivers at other observation wells. Analysis of direct arrivals of emulated cross-well seismic signals can be completed to extract travel times, inter-well seismic velocities, and rock properties.

An automatic microseismic monitoring-intelligent rockburst early warning integrated method is further provided.

A distributed acoustic sensor is positioned within a wellbore of a geologic formation. Seismic waves are detected using the distributed acoustic sensor. A raw seismic profile is generated based on the detected seismic waves. Resonant noise is identified and reduced in seismic data associated with the raw seismic profile.

A method for determining microseismic events. The method may include measuring a seismic travel time of a microseismic event with a fiber optic line disposed in a first wellbore, forming a probability density function for the microseismic event based at least in part on the seismic travel time measurement, modifying the probability density function by applying one or more constraints to form a modified probability density function, identifying one or more most probable source locations from the modified probability density function, and forming a microseismic event cloud from the one or more most probable source locations.

A data seismic sensing system and method for obtaining seismic refraction data and tomography data. The system may comprise a subsurface sensor array, wherein the subsurface sensor array is a fiber optic cable disposed near a wellbore, a seismic source, wherein the seismic source is a truck-mounted seismic vibrator comprising a base plate, and a surface sensor array, wherein the surface sensor array is coupled to the seismic source. The method may comprise disposing a surface sensor array on a surface, disposing a subsurface sensor array into a wellbore, activating a seismic source, wherein the seismic source is configured to create a seismic wave, recording a reflected seismic wave with the surface sensor array and the subsurface sensor array, and creating a seismic refraction data and a seismic tomography data from the reflected seismic wave.

Methods to estimate formation shear wave slowness from multi-firings of different types of acoustic sources and multi-mode dispersion estimation systems are presented. The method includes obtaining waveform data of waves traversing through a downhole formation, where the waves are generated from multi-firings of different types of acoustic sources. The method also includes performing a multimode dispersion analysis of the waveform data for each firing of the multi-firings, and removing one or more tool waves generated from the multi-firings. The method further includes determining a formation type of the formation the waves traverse based properties of the waves and determining an initial shear wave slowness estimate of the waves. The method further includes generating a modeling of the waves, and reducing a mismatch between the modeling of the waves and a slowness dispersion of the waves to improve the modeling of the waves.

A method of determining a shear-wave attenuated vertical component vertical seismic profile (VSP) dataset is disclosed. The method includes, obtaining a multi-component VSP dataset, including a vertical and a horizontal component, transforming the vertical component into a vertical spectrum and the horizontal component into a horizontal spectrum, and designing a band-pass filter based, at least in part, on an energetic signal of the horizontal spectrum. The method further includes determining a muted vertical amplitude spectrum by applying the pass-band filter to an amplitude spectrum of the vertical spectrum, determining an estimated noise model based on the muted vertical amplitude spectrum and the vertical spectrum; and determining the shear-wave attenuated vertical component VSP dataset by adaptively subtracting the estimated noise model from the vertical component of the multi-component VSP dataset. A system including a seismic source, a plurality of seismic receivers, and a seismic processor for executing the method is disclosed.

A method for performing a formation-related operation based on corrected vertical seismic profile (VSP) data of an earth formation includes performing a VSP survey and applying a normal moveout (NMO) correction equation to the survey data that is a function of source offset to wellhead. The method also includes solving the NMO correction equation using a simulated annealing algorithm having an object function that is a coherence coefficient of semblance analysis of an NMO corrected reflection event within a time window to provide NMO corrected data. The method further includes performing the formation-related operation at at least one of a location, a depth and a depth interval based on the VSP NMO corrected data.

A method for performing a formation-related operation based on corrected vertical seismic profile (VSP) data of an earth formation includes performing a VSP survey and applying a normal moveout (NMO) correction equation to the survey data that is a function of source offset to wellhead. The method also includes solving the NMO correction equation using a simulated annealing algorithm having an object function that is a coherence coefficient of semblance analysis of an NMO corrected reflection event within a time window to provide NMO corrected data. The method further includes performing the formation-related operation at at least one of a location, a depth and a depth interval based on the VSP NMO corrected data.

A seismic attenuation quality factor Q is determined for seismic signals at intervals of subsurface formations between a seismic source at a marine level surface and one or more receivers of a well. Hydrophone and geophone data are obtained. A reference trace is generated from the hydrophone and geophone data. Vertical seismic profile (VSP) traces are received. First break picking of the VSP traces is performed. VSP data representing particle motion measured by a receiver of the well are generated. The reference trace is injected into the VSP data. A ratio of spectral amplitudes of a direct arrival event of the VSP data and the reference trace is determined. From the ratio, a quality factor Q is generated representing a time and depth compensated attenuation value of seismic signals between the seismic source at the marine level surface and the first receiver.