A seismic modeling system is provided which may include a seismic model data storage device, and a processor cooperating with the seismic model data storage device to identify an obstruction within a seismic model spatial domain data set having a lateral boundary, and define a suspect area including spatial domain data within the lateral boundary and directly below the obstruction. The processor may further inpaint the suspect area in the seismic model spatial domain data set based upon an exemplar inpainting algorithm.

An example method for displaying sonic logging data associated with a formation surrounding a borehole can include acquiring sonic data at a plurality of depths using an acoustic array located in the borehole and transforming the acquired sonic data from a time-space domain to a frequency-wave number domain at a limited number of discrete frequencies. The method can also include estimating slowness values at the limited number of discrete frequencies from the transformed sonic data, interpolating the estimated slowness values to obtain a projection of one or more slowness-frequency dispersions of the acquired sonic data and displaying the projection of the slowness-frequency dispersions. The projection of the slowness-frequency dispersions can include a plurality of color bands corresponding to each of the limited number of discrete frequencies.

This disclosure describes processes and systems for generating a high-resolution velocity model of a subterranean formation from recorded seismic data gathers obtained in a marine seismic survey of the subterranean formation. A velocity model is computed by iterative FWI using reflections, resolving the velocity field of deep subterranean targets without requiring ultralong offsets. The processes and systems use of an impedance sensitivity kernel to characterize reflections in a modeled wavefield, and then use the reflections to compute a velocity sensitivity kernel that is used to produce low-wavenumber updates to the velocity model. The iterative process is applied in a cascade such that position of reflectors and background velocity are simultaneously updated. Once the low-wavenumber components of the velocity model are updated, the velocity model is used as an input of conventional FWI to introduce missing velocity components (i.e., high-wavenumber) to increase the resolution of the velocity model.

A data visualization method for visualizing data acquired along a non-linear acquisition path or sail line. The data consists of CMP lines that follow the non-linear acquisition path. The data is arranged such that the in-lines in the binning grid follow the acquisition path and the cross-lines are perpendicular, or near perpendicular, to the in-lines, the method comprising the steps of: creating a binning grid covering the CMP lines of the acquired data, the binning grid comprising a straight portion and a curved portion; calculating bins for each portion; loading the seismic data into the a visualization software; and creating a set of linked windows, wherein a field of view of the different set of linked windows is synchronized, and wherein a marker is provided to visualize the field of view of data in at least two of the linked windows.

Seismic image filtering machines, systems, program products, and computer implemented methods are provided to generate a filtered seismic image responsive to filtered seismic image data generated by attenuating coherent seismic noise from surface waves of an unfiltered wavefield constructed from unfiltered seismic image data through a single downward extrapolation of the unfiltered wavefield using a plurality of nonstationary convolution operators to perform localized filtering at each of a plurality of spatial locations of the unfiltered wavefield. Various embodiments, for example, can beneficially handle strong lateral velocity variations thus making various embodiments effective tools to remove complicated coherent seismic noise which is typically in the form of exponentially decaying evanescent waves. Embodiments of the present invention, for example, can use, as a part of the filtering mechanism, specially designed nonstationary convolution operators that are implemented in the space-frequency domain as nonstationary filters.

A method for displaying slowness data after defects of anisotropy and heterogeneity have been removed includes visually depicting a slowness image on a computer monitor in which the slowness image includes a log of slowness values at a plurality of depths of a wellbore and a plurality of azimuth angles about a periphery of the wellbore. The image is processed via filtering such as with a bandpass filter. The filter may be obtained via fitting a periodic function to the slowness values at each depth or via a selected periodicity value and center angle. A resultant image is displayed.

Systems and methods for geochemical sampling grid locations on a seafloor. At least one of the methods includes generating, using received seismic data, an image representing an interpretation of a seafloor horizon surface; extracting, from the image and based on the seismic data, one or more discontinuity attributes of the seafloor horizon surface; extracting, from the image and based on the seismic data, one or more amplitude attributes of a window extending below the seafloor horizon surface; combining the one or more discontinuity attributes and the one or more amplitude attributes; and selecting, using the image and based at least partly on the combining, one or more locations of the seafloor horizon surface for sampling.

Disclosed are systems and methods for high precision acoustic logging processing for compressional and shear slowness. The method comprises measuring, by a sonic logging tool, sonic data associated with a formation within a borehole, attempting a detection of a first arrival within the sonic data determining whether the attempted detection of the first arrival is accurate, and in response to an accurate detection of the first arrival determining a travel time of the first arrival, generating a coherence map including the first arrival, and determining, based on the coherence map, a characteristic of the formation.

A method, computer program product, and computing system for generating high resolution slowness logs. The method computer program product, and computing system includes receiving a plurality of sonic logs from at least one sensor array and generating at least one high-resolution slowness log from the plurality of sonic logs based upon, at least in part, monopole and dipole data from the plurality of sonic logs.

According certain aspects, embodiments of the invention consider the problem of microseismic event localization from a parameter estimation perspective, and include a method and system for computing and displaying characteristics of event localization errors. According to certain other aspects, embodiments of the invention include techniques for deriving aggregate statistics from a set of event location estimates, including methods for computing and displaying the probability that an event occurred in any given volume, and methods for describing and displaying the smallest volume that contains a specified percentage of the event probability or expected to contain the specified percentage of the events.