Abstract: Systems and methods for estimating reservoir productivity as a function of position in a subsurface volume of interest are disclosed. The systems and methods may: obtain subsurface data and well data corresponding to a subsurface volume of interest; obtain a parameter model; use the subsurface data and the well data to generate multiple production parameter maps; apply the parameter model to the multiple production parameter maps to generate refined production parameter values; generate multiple refined production parameter graphs; display the multiple refined production parameter graphs; generate one or more user input options; receive a defined well design and the one or more user input options selected by a user to generate limited production parameter values; generate a representation of estimated reservoir productivity as a function of position in the subsurface volume of interest using the defined well design and visual effects; and display the representation.

Abstract: Systems, methods, and storage media for detecting a given subsurface event in a subsurface volume of interest are disclosed. Exemplary implementations may: receive a subsurface fiber optic data set; receive a sensor data set using one or more sensors; constrain the subsurface fiber optic data set based on a given parameter value of a given parameter within a certain range to generate a constrained subsurface fiber optic data set; use sets of models to refine the constrained subsurface fiber optic data set to generate a refined subsurface fiber optic data set; estimate an event location of the given subsurface event based on the refined subsurface fiber optic data set; and estimate an origin time based on the event location.

Abstract: A method is described for processing residual moveout in seismic image data gathers representing critical reflections. The method includes receiving seismic image data arranged as a function of an angle or offset parameter including a high-velocity-contrast event with post-critical. The method also includes applying a wavelet de-stretch filter to the seismic data to correct wavelet stretching. The method also includes applying a fan-filter to remove coherent noise in the one or more post-critical traces; picking residual moveout of the high-velocity-contrast event; adjacent-trace differencing to detect the impact of phase change at critical reflections in residual moveouts, and applying a median-filter to the residual moveout to reduce the impact of phase change of the high-velocity-contrast event in the one or more post-critical traces. The median-filtered and reconstructed residual moveout is used for improving a velocity model used for generating the seismic image gathers.

Abstract: A method is described for seismic imaging of complex subsurface volumes including an optimized partial stack based on ray coverage determined for vector image partition tiles. The method may be executed by a computer system.

Abstract: A method is described for residual moveout analysis using dynamic warping. The residual moveout curves or surfaces may be used to flatten common image gathers or derive velocity models. The method may be executed by a computer system.

Abstract: A method is described for processing residual moveout in seismic image data gathers representing critical reflections. The method includes receiving seismic image data arranged as a function of an angle or offset parameter including a high-velocity-contrast event with post-critical. The method also includes applying a wavelet de-stretch filter to the seismic data to correct wavelet stretching. The method also includes applying a fan-filter to remove coherent noise in the one or more post-critical traces; picking residual moveout of the high-velocity-contrast event; adjacent-trace differencing to detect the impact of phase change at critical reflections in residual moveouts, and applying a median-filter to the residual moveout to reduce the impact of phase change of the high-velocity-contrast event in the one or more post-critical traces. The median-filtered and reconstructed residual moveout is used for improving a velocity model used for generating the seismic image gathers.

Abstract: A method for determining values of anisotropic model parameters of a Tilted Transversely Isotropic (TTI) Earth model, the method including obtaining an initial TTI earth model that substantially flattens common-imaging-point gathers and substantially ties seismic data to well data; inputting checkshot data and/or VSP data to determine updated values of Vp0 near the well locations; determining an incremental improvement ??; extrapolating the relative change ?? from near-well locations to the entire three dimensional TTI earth model; determining updated values of Vp0=Vp0 (1???); inputting near-to-mid-offset/angle and mid-to-far-offset/angle residual moveout information; and providing updated values of ? and ?.

Abstract: A method for determining values of anisotropic model parameters of a Tilted Transversely Isotropic (TTI) Earth model, the method including obtaining an initial TTI earth model that substantially flattens common-imaging-point gathers and substantially ties seismic data to well data; inputting checkshot data and/or VSP data to determine updated values of Vp0 near the well locations; determining an incremental improvement ??; extrapolating the relative change ?? from near-well locations to the entire three dimensional TTI earth model; determining updated values of Vp0=Vp0 (1???); inputting near-to-mid-offset/angle and mid-to-far-offset/angle residual moveout information; and providing updated values of ? and ?.

Abstract: Seismic data are processed to reduce or eliminate aliasing due, for example to sparse or irregular sampling. An iterative method includes an inhibiting function used in conjunction with a function evaluating a magnitude of Fourier coefficients that together act to reduce the effects of aliased energies and preferentially select true energies. Computational steps are conducted primarily in k-space, without returning to x-space, thereby reducing computational costs.

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: Seismic data are processed to reduce or eliminate aliasing due, for example to sparse or irregular sampling. An iterative method includes an inhibiting function used in conjunction with a function evaluating a magnitude of Fourier coefficients that together act to reduce the effects of aliased energies and preferentially select true energies. Computational steps are conducted primarily in k-space, without returning to x-space, thereby reducing computational costs.