Methods and systems for determining an interbed multiple attenuated pre-stack seismic dataset are disclosed. The methods include forming a post-stack seismic image composed of post-stack traces from the pre-stack seismic dataset and identifying a first, second, and third post-stack horizon on each of the post-stack traces. The methods further include for each pre-stack trace, generating a first, second, and third multiple-generator trace based on the first, second and third post-stack horizon and determining a correlation trace based, at least in part, on a correlation between the first multiple-generator trace and the second multiple-generator trace. The methods still further include predicting an interbed multiple trace by convolving the correlation trace and the third multiple-generator trace, determining an interbed multiple attenuated trace by subtracting the interbed multiple trace from a corresponding pre-stack seismic trace, and determining the interbed multiple attenuated pre-stack seismic dataset by combining the interbed multiple attenuated traces.

Systems and a method are disclosed. The method includes obtaining a plurality of raw seismic datasets for a subterranean region of interest, wherein each raw seismic dataset is generated by a high-speed train traversing a train track at a unique speed. The method further includes determining a plurality of processed seismic datasets by processing each of the plurality of raw seismic datasets and determining a final seismic dataset by combining the plurality of processed seismic datasets. The method still further includes identifying subterranean features within the subterranean region of interest using the final seismic dataset.

There is disclosed in the present disclosure a seismic full horizon tracking method, a computer device and a computer-readable storage medium. The method includes: acquiring three-dimensional seismic data; extracting horizon extreme points from the three-dimensional seismic data to construct a sample space; equally dividing the sample space into a plurality of sub-spaces with overlapping portions, and performing a clustering process on the horizon extreme points in each sub-space to obtain horizon fragments corresponding to each horizon of the three-dimensional seismic data; establishing a topological consistency between the horizon fragments; and fusing the horizon fragments corresponding to each horizon of the three-dimensional seismic data based on the topological consistency, to obtain a full horizon tracking result of the three-dimensional seismic data. In the disclosure, a layer crossing phenomenon occurring in seismic full horizon tracking can be avoided, and a better full horizon tracking effect can be achieved.

A method for modeling a subsurface volume includes receiving a plurality of ordered seismic images including representations of objects in the subsurface volume, generating flow fields based on a difference between individual images of the plurality of ordered seismic images, and identifying the objects in the seismic images based on the flow fields and the plurality of ordered seismic images.

In seismic imaging, accurate velocity functions (velocity model) defining seismic velocity as a function of depth in the earth are required. The velocity model is obtained as a result of seismic surveying. Delayed travel times in near surface refraction seismic surveys, an effect known as shingling, can result from an anomalous condition, seismic velocity decreasing with depth. Inclusion of such delayed travel times in a tomographic process for seismic imaging would otherwise cause large errors in determination of a seismic velocity model for seismic imaging of subsurface features. At locations (source-receiver offset) in the survey where the shingling occurs, the velocity inversions are identified. The undesirable effects the delayed travel times caused by the velocity inversions are removed from the survey dataset.

Techniques to reduce noise in seismic data by receiving a set of seismic data comprising a plurality of input volumes each inclusive of positional data and at least one additional attribute related to the seismic data, selecting a first input volume of the plurality of input volumes having a first additional attribute related to the seismic data, and generating a pilot volume by selecting a range of input volumes of the plurality of input volumes and stacking input volumes of the range of input volumes with the first input volume. Additionally, generating a trained dictionary based upon transformation of the pilot volume, transforming the first input volume into transformed data, imposing a sparse condition on the transformed data utilizing the trained dictionary to generate sparsified data, and inverse transforming the sparsified data to generate an output data volume as a portion of a set of modified seismic data.

In contrast to existing methods wherein derived horizons are interpreted in isolation, the disclosure provides a process that does not interpret patches themselves but determines the relationships between patches, in order to associate and link patches to derive a holistic geological interpretation. Predefined patches, such as from a pre-interpreted suite, are received as inputs to determine the relationships and derive an interpretation for a complete volume. In one aspect the disclosure provides an automated method of generating a geological age model for a subterranean area. In one example, the automated method includes: (1) abstracting seismic data of a subsurface into a limited number of patches, (2) abstracting the patches by defining patch-links between the patches, and (3) generating a geological age model of the subsurface by solving for the relative geological age of each of the patches using the patch-links.

A system, computer-readable medium, and method for determining a potential drilling location, of which the method includes obtaining data representing a subterranean domain. The data includes at least seismic data. The method also includes inverting the seismic data, creating a petroleum systems model of the subterranean domain based at least in part on a result of inverting the seismic data, simulating a dynamic reservoir model of the subterranean domain based at least in part on the petroleum systems model, and identifying the potential drilling location based on a combination of the inverting of the seismic data, creating the petroleum systems model, and simulating the dynamic reservoir model.

A system, computer-readable medium, and method for determining a potential drilling location, of which the method includes obtaining data representing a subterranean domain. The data includes at least seismic data. The method also includes inverting the seismic data, creating a petroleum systems model of the subterranean domain based at least in part on a result of inverting the seismic data, simulating a dynamic reservoir model of the subterranean domain based at least in part on the petroleum systems model, and identifying the potential drilling location based on a combination of the inverting of the seismic data, creating the petroleum systems model, and simulating the dynamic reservoir model.

Methods of and systems for forming an image of a subterranean region of interest are disclosed. The method includes obtaining an observed seismic dataset and a seismic velocity model for the subterranean region of interest and generating a simulated seismic dataset based on the seismic velocity model and the source and receiver geometry of the observed seismic dataset. The method also includes forming a plurality of time-windowed trace pairs from the simulated and the observed seismic datasets, and forming an objective function based on a penalty function and a cross-correlation between the members of each pair. The method further includes determining a seismic velocity increment based on the extremum of the objective function and forming an updated seismic velocity model by combining the seismic velocity increment and the seismic velocity model, and forming the image of the subterranean region of interest based on the updated seismic velocity model.