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.

Near real-time methodologies for maximizing return-on-fracturing-investment for shale fracturing. An example system can calculate, based on sonic data and density data, mechanical properties and closure stress of a portion of shale rocks for fracture modeling. The system can generate one or more rock mechanical models based on the mechanical properties and closure stress of the portion of shale rocks, and perform one or more fracture modeling simulations based on one or more treatment parameter values. Based on the one or more fracture modeling simulations, the system can generate a neural network model which predicts a fracture productivity indicator of an effective propped area (EPA) and/or an effective propped length (EPL), and calculate a return-on-fracturing-investment (ROFI) based on the EPA or EPL predicted by the neural network model.

A portion of an anisotropy formation through which a wellbore is formed can be identified. An estimate of an elastic anisotropy parameter for the portion can be adjusted based on a first quality control analysis using the elastic anisotropy parameter for the portion. The first signal representing the elastic anisotropy parameter for the portion. The estimate of the elastic anisotropy parameter for the portion can be adjusted based on a second quality control analysis using estimates for the elastic anisotropy parameters for two or more portions of the anisotropy formation.

Systems and methods for estimating orthorhombic anisotropy parameters of subsurface rock layers are provided. An initial three-dimensional (3D) model of layers in a subsurface formation is generated. Particular combinations of seismic source and receiver locations associated with a vertical seismic profile (VSP) survey of the formation are selected based on a simulation of seismic wave propagation through each layer of the initial 3D model. A global inversion is performed using data points selected from travel time data associated with the VSP survey in order to estimate different sets of anisotropy parameters for the layers of the initial 3D model, where the selected data points correspond to the selected source-receiver combinations. The initial 3D model is refined with an optimal set of anisotropy parameters selected from the estimated parameter sets. The refined 3D model is provided for seismic analysis and well planning operations to be performed for the subsurface formation.

Fracture orientation ambiguity in the results of amplitude variation with offset and azimuth inversion is overcome by using additional geological information. Methods, apparatus and executable codes stored on non-transitory media cause, for each interface encountered by traces in a CMP gather, calculating angle-dependent azimuthal Fourier coefficients, performing a nonlinear inversion of amplitude versus offset and azimuth, AVOAz, equations built using the angle-dependent azimuthal Fourier coefficients to determine possible fracture orientations, and selecting one of the possible fracture orientations using constraints based on the additional geological information.

Various examples are provided for wavefield extrapolation in anisotropic media. In one example, among others, a method includes determining an effective isotropic velocity model and extrapolating an equivalent propagation of an anisotropic, poroelastic or viscoelastic wavefield. The effective isotropic velocity model can be based upon a kinematic geometrical representation of an anisotropic, poroelastic or viscoelastic wavefield. Extrapolating the equivalent propagation can use isotopic, acoustic or elastic operators based upon the determined effective isotropic velocity model. In another example, non-transitory computer readable medium stores an application that, when executed by processing circuitry, causes the processing circuitry to determine the effective isotropic velocity model and extrapolate the equivalent propagation of an anisotropic, poroelastic or viscoelastic wavefield. In another example, a system includes processing circuitry and an application configured to cause the system to determine the effective isotropic velocity model and extrapolate the equivalent propagation of an anisotropic, poroelastic or viscoelastic wavefield.

Example computer-implemented method, computer-readable media, and computer system are described for generating subterranean imaging data based on initial isotropic and/or anisotropic velocity models for the vertical seismic profile (VSP) data and stored ocean bottom sensor (OBS) data. In some aspects, VSP data and OBS data of a subterranean region are received. Angle attributes for each image point are computed to image primary reflection and free surface multiples of the received VSP data and OBS data, respectively. Angle-domain common-image gathers (ADCIG) are generated according to a ray-equation method based on the angle attributes computed based on the received VSP data and OBS data, respectively. The ADCIG are further post-processed.

Example computer-implemented method, computer-readable media, and computer system are described for generating subterranean imaging data based on vertical seismic profile (VSP) data. In some aspects, VSP data of a subterranean region can be received. Four angle attributes for each image point can be computed based on the received VSP data. Five-dimensional (5D) angle-domain common-image gathers (ADCIG) can be generated according to a ray-equation method based on the four angle attributes.

Systems and methods for imaging seismic data using hybrid one-way wave-equation-migration in tilted transverse isotropic media and/or hybrid two-way reverse-time-migration in tilted transverse isotropic media.

The present disclosure provides a method for estimating the elastic constants of an anisotropic material such as anisotropic rocks such as gneiss and shale even when using a general loader, the method enabling anisotropic elastic constants to be estimated even with only a single core sample, thereby significantly reducing time and costs.