Multicomponent data are acquired using a downhole acoustic tool having transmitters and receiver stations distributed azimuthally in a plane perpendicular to the axis of the tool. The receiver stations are located at several receiving stations along the axis of the tool. At each acquisition depth, waveforms are processed through a multi-dimensional fast Fourier transform, extrapolation and inverse multi-dimensional fast Fourier transform. At each receiver station, waveforms are combined to produce the standard monopole waveforms and the inline and crossline dipole waveforms along fixed azimuths. These oriented waveforms produce a finer azimuthal sampling of the surrounding formation, and can then be used for imaging geological features within the surrounding formation.

Disclosed are methods, systems, and computer-readable medium to perform operations including: receiving for a plurality of common midpoint-offset bins each comprising a respective plurality of seismic traces, respective candidate pilot traces representing the plurality of common midpoint-offset bins; generating, based on the respective candidate pilot traces, a respective plurality of corrected seismic traces for each of the plurality of common midpoint-offset bins; grouping the respective pluralities of corrected seismic traces into a plurality of enhanced virtual shot gathers (eVSGs); generating, based on the plurality of common midpoint-offset bins, a common-midpoint (CMP) velocity model; calibrating the CMP velocity model using uphole velocity data to generate a pseudo-3 dimensional (3D) velocity model; performing, based on the plurality of enhanced virtual shot gathers and the pseudo-3D velocity model, a 1.5-dimensional full waveform inversion (FWI); and determining the subsurface velocity model based on the 1.5 dimensional FWI.

A sonic logging method is provided that transmits acoustic signals using a high order acoustic source and processes waveform data to identify a set of arrival events and time picks by automatic and/or manual methods. Ray tracing inversion is carried out for each arrival event over a number of possible raypath types that include at least one polarized shear raypath type to determine two-dimensional reflector positions and predicted inclination angles of the arrival event for the possible raypath types. One or more three-dimensional slowness-time coherence representations are generated for the arrival event and raypath type(s) and evaluated to determine azimuth, orientation and raypath type of a corresponding reflector. The method outputs a three-dimensional position and orientation for at least one reflector. The information derived from the method can be conveyed in various displays and plots and structured formats for reservoir understanding and also output for use in reservoir analysis and other applications.

Methods and systems including computer programs encoded on a computer storage medium, for utilizing equivalent linear velocity for first arrival picking of seismic refraction. In one aspect, a method includes receiving data for the shot gather record, generating a diving wave equation curve for a particular parameter pair of multiple parameter pairs, and integrating the shot gather record data corresponding to the diving wave equation curve over a selected range of offsets of the shot gather to generate an equivalent linear velocity value for the particular parameter pair and the shot gather record data, selecting, from the equivalent linear velocity values for the plurality of parameter pairs, a greatest equivalent linear velocity value of the equivalent linear velocity values, the greatest equivalent linear velocity value corresponding to a first-arrival parameter pair, and determining, using the first-arrival parameter pair, a set of first-arrival onsets for the selected sub-range of offsets.

A computer-implemented method and computing system apparatus programmed to perform operations of the computer-implemented method for obtaining a subsurface stack image, subsurface angle gathers, and a subsurface velocity model over an entire survey region having high velocity contrast geo-bodies. Particularly, user inputs, input velocity models, and surface-seismic data are obtained by fixed source and receiver pairs and then used by the computer program product embedded within the computing system apparatus to minimize the number of iterations, required to obtain a final velocity model, a final stack image, and final angle gathers wherein their flatness deviation is equal to, or less than, a user-defined flatness value. Therefore, the attributes developed by said computer-implemented method and system can help solve the imaging problem of sub high velocity contrast geo-bodies like subsalt, or salt overhung deep mini basins.

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.

Waves from cement bond logging with a sonic logging-while-drilling tool (LWD-CBL) are often contaminated with tool waves and may yield biased CBL amplitudes. The disclosed LWD-CBL wave processing corrects the first echo amplitudes of LWD-CBL before calculating the BI. The LWD-CBL wave processing calculates a tool wave amplitude and a phase angle difference as the difference of the phases between the tool waves and casing waves. The tool waves are then used to correct the LWD-CBL casing wave amplitude and remove errors introduced from tool waves. In conjunction with the sets of operations described, the LWD-CBL wave processing also include array preprocessing operations. Array preprocessing may employ variation of bandpass filtering and frequency-wavenumber (F-K) filtering operations to suppress tool wave.

Downhole measurement systems and methods include deploying a bottomhole assembly having a multipole transmitter into a formation and transmitting acoustic signals into the formation. The multipole transmitter is of order n ≥ 2. Acoustic signals are received at respective receivers that are circumferentially aligned with the multipole transmitter, and are axially offset from the multipole transmitter, and axially offset from each other. The order of the first and second multipole receivers are equal to the order of the multipole transmitter. A controller is used to obtain first and second acoustic multipole data from the first and second multipole receivers at one or more azimuthal angles of a rotation of the bottomhole assembly in a formation during a drilling operation. Acoustic azimuthal anisotropy of the formation is determined from the first acoustic multipole data and the second acoustic multipole data.