A method can include receiving spatially located geophysical data of a geologic region as acquired by one or more sensors; solving a system of equations for multi-dimensional implicit function values within a multi-dimensional space that represents the geologic region where the system of equations are subject to a smoothness constraint and subject to a weighted curvature minimization criterion at a plurality of spatially located points based on the spatially located geophysical data; and rendering to a display, a structural model of the geologic region based at least in part on the multi-dimensional implicit function values where the structural model characterizes stratigraphy of the geologic region.

A system includes a processor and a memory. The memory includes instructions that are executable by the processor to access a plurality of seismic images of a subterranean formation in a first geological area. The instructions are also executable to generate a plurality of fault estimates from each of the plurality of seismic images. Further, the instructions are executable to generate a processed seismic image of the first geological area by normalizing and merging the plurality of seismic images and the plurality of fault estimates. Additionally, the instructions are executable to generate a statistical fault uncertainty volume of the first geological area using the processed seismic image. Furthermore, the instructions are executable to control a drilling operation in the first geological area using the statistical fault uncertainty volume of the first geological area.

Prestack images from the object store are hierarchically combined to generate a hierarchically stacked image. The hierarchically stacked image is generated by combining stacked images that includes a stacked image. The stacked image is generated by combining at least the prestack images. Based at least on the hierarchically stacked image, a quality measure of a prestack image is generated. Prior to deleting at least a subset of the prestack images from the object store and based at least on the quality measure, the prestack images are further combined to generate an enhanced stacked image. The stacked image is substituted using the enhanced stacked image. Subsequent to the substituting and prior to deleting at least the subset of the stacked images from the object store, the stacked images are combined to generate an enhanced hierarchically stacked image. The enhanced stacked image and the enhanced hierarchically stacked image are generated using failure recovery metadata. The enhanced hierarchically stacked image is presented.

A method for imaging a downhole formation. The method includes combining the captured images to generate a partial image of the formation, wherein the partial image includes captured images separated by gaps representing portions of the formation not captured with sensors what were disposed downhole. The method includes locating dips in the formation within the partial image and interpolating the partial image using the located dips within the partial image.

A method of compressing data from seismic waves using Gabor frames utilizes a plurality of geophones positioned within a region of interest. Each of the plurality of geophones is communicably coupled with at least one remote server. Thus, a plurality of reflected-seismic signals received through the plurality of geophones can be transmitted to the at least one remote server for analyzing and calculations. The plurality of reflected-seismic signals is converted into a set of Gabor frames, wherein the Gabor frames is generated via a plurality of prolate spheroidal wave functions (PSWF). A Gabor frame-generating calculation module utilizes the plurality of PSWF to generate the set of Gabor frames. A dual frame for each of the set of Gabor frames is derived and used for quantization purposes. Preferably, a tree structured vector quantization process is followed.

A method of compressing data from seismic waves using Gabor frames utilizes a plurality of geophones positioned within a region of interest. Each of the plurality of geophones is communicably coupled with at least one remote server. Thus, a plurality of reflected-seismic signals received through the plurality of geophones can be transmitted to the at least one remote server for analyzing and calculations. The plurality of reflected-seismic signals is converted into a set of Gabor frames, wherein the Gabor frames is generated via a plurality of prolate spheroidal wave functions (PSWF). A Gabor frame-generating calculation module utilizes the plurality of PSWF to generate the set of Gabor frames. A dual frame for each of the set of Gabor frames is derived and used for quantization purposes. Preferably, a tree structured vector quantization process is followed.

Computing device, computer instructions and method for identifying seismic traces prone to cycle-skipping in a full waveform inversion method. The method includes receiving recorded seismic data recorded with seismic sensors over a subsurface of interest; selecting a model that describes the subsurface; calculating, based on the model and the recorded seismic data, estimated seismic data; and choosing a probabilistic measure that characterizes a relationship between the recorded seismic data and the estimated seismic data. The probabilistic measure includes at least one statistical function.

A method, apparatus, and program product render a seismic slice in real time and in a computationally-efficient manner using a displacement mapping technique implemented in one or more programmable shaders of a graphics processing unit (GPU), e.g., using programmable shaders in a GPU to perform both tessellation and displacement of primitives in connection with rendering a displacement-mapped visualization of the seismic slice for display in an interactive 3D visualization environment.

A method for deblending seismic signals includes entering as input to a computer recorded signals comprising seismic energy from a plurality of actuations of one or more seismic energy sources. A model of deblended seismic data and a blending matrix are initialized. A blending matrix inversion is performed using the initialized model. The inversion includes using a scaled objective function. The inversion is constrained by a thresholding operator. The thresholding operator is arranged to recover coefficients of the model of the deblended seismic data that are substantially nonzero, against a Gaussian white noise background. The thresholded model is projected into data space. Performing the blending matrix inversion is repeated if a data residual exceeds a selected threshold and the inversion is terminated if the data residual is below the selected threshold. At least one of storing and displaying an output of the blending matrix inversion is performed when the blending matrix inversion is terminated.

A method may include obtaining an initial velocity model regarding a subterranean formation of interest. The method may further include generating various seismic migration gathers with different cross-correlation lag values based on a migration-velocity analysis and the initial velocity model. The method may further include selecting a predetermined cross-correlation lag value automatically using the seismic migration gathers and based on a predetermined criterion. The method may further include determining various velocity boundaries within the initial velocity model using a trained model, wherein the trained model is trained by human-picked boundary data and augmented boundary data. The method may further include updating, by the computer processor, the initial velocity model using the velocity boundaries, the automatically-selected cross-correlation lag value, and the migration-velocity analysis to produce an updated velocity model. The method may further include generating an image of the subterranean formation of interest using the updated velocity model.