A method for seismic processing includes steps of seismic signal forward propagation and seismic data back propagation. The subsurface medium image is created after correlating and summarizing forward and backward propagation results. To address migration footprint and noise due to the incomplete data acquisition aperture and migration approximation in the migration operator, the iteration inversion strategy incorporates tensor flow calculated from seismic image. A regularization operator based on structure tensor of image is applied to seismic image inversion.

Systems and methods for simulating subterranean regions having multi-scale fracture geometries. Non-intrusive embedded discrete fracture modeling formulations are applied in conjunction with commercial simulators to efficiently and accurately model subsurface transient flow characteristics in regions having complex hydraulic fractures, complex natural fractures, or a combination of both, and geometries including corner point grids.

The present disclosure describes methods and systems for interpreting geological features in a seismic volume based on mono-frequency filtering of the seismic volume. One computer-implemented method includes receiving a seismic data volume, decomposing the seismic data volume into multiple sub-volumes, generating one or more seismic horizons on each sub-volume, analyzing the generated seismic horizons for the multiple sub-volumes including determining a first sub-volume and a second sub-volume from the multiple sub-volumes, and subtracting the generated one or more seismic horizons for the first sub-volume from the generated one or more seismic horizons for the second sub-volume.

Systems, methods, and computer-readable medium for generating a three-dimensional fracture network model are provided. The method can include receiving reflected acoustic signal measurements acquired in response to emission of acoustic waves by one or more sensors disposed in a wellbore formed within a target region. Each reflected acoustic signal measurement represents a strength of a reflected acoustic wave as a function of time measured in at least one predetermined direction oriented with respect to an axis of the wellbore. A fracture extension estimate is generated for each of the reflected acoustic signal measurements. A three-dimensional fracture network model is generated corresponding to the fracture extension estimates generated for each of the plurality of reflected acoustic measurements. The generated fracture network model is output for display or use in modeling environments.

Methods, systems, and apparatus, including computer programs encoded on a computer storage medium, to identify geobodies for exploration and production activities. In one aspect, a method includes receiving seismic volume data regarding a geobody, generating an interpretation object by applying a multi-Z interpretation at representative intervals of edge-detected attributes through the seismic volume data, converting Z-value surfaces to points at each of the representative intervals, the Z-value surfaces generated by gridding polygons having been infilled and converted from the interpretation object, applying a mathematical addition or subtraction to generate a body of points representing the geobody by replicating the points at successive intervals, forming a blended 3D property model of the geobody by combining the seismic volume data resampled into a three-dimensional (3D) cellular grid, and enabling exploration activities by employing the seismic attribute analysis of the geobody identified based on the blended 3D property model.

A method is described for seismic modeling implemented on a computer system including multiple computer nodes, which includes receiving, at a computer processor, a 4D speed of sound in water library and a marine seismic survey geometry; creating, via the computer processor, a full survey 3D Earth model for each seismic shot time based on the 4D speed of sound in water library to generate a set of dynamic 3D Earth models; simulating 4D seismic shots for the set of dynamic 3D Earth models by simulating multiple shots on each node to generate a set of 4D seismic shot gathers; and outputting the 4D seismic shot gathers.

A method for identifying a location of a distributed acoustic system channel in a distributed acoustic system. The method may comprise generating a two or three dimensional layer model interface with an information handling system, preparing a P-wave first arrival pick time table, estimating an initial model layer properties, estimating a location of the distributed acoustic system channels, preparing an overburden file of layer properties, running an anisotropic ray tracing, defining an upper and a lower limits for model parameters, specifying parameters for the inversion, running an inversion, selecting a solution based at least in part on stored error predictions, and calculating a mean and a standard deviation of an inverted model parameter.

Computer vision systems and methods for detecting anomalous building models are provided. The systems and methods can detect anomalies in building models using one or more of an independent univariate Gaussian algorithm, a multivariate Gaussian algorithm, a combination of a multivariate Gaussian algorithm for continuous features and a frequency histogram algorithm for discrete features, and/or a bin frequency model. The system automatically processes computerized models to determine anomalies, and indicates whether the models are accurate and whether correction is required.

A method to select a representative subset of a plurality of horizon surfaces or surface patches from geophysical subsurface imaging data, including: defining a score function on one or more horizon surfaces or surface patches; calculating, by a computer, the score for each of the plurality of horizon surfaces or surface patches with regard to other horizon surfaces or surface patches and whether the other horizon surfaces or surface patches have been selected or not for inclusion or exclusion in the subset of the plurality of horizon surfaces; selecting, by a computer, one or more of the plurality of horizon surfaces or surface patches to be included in the subset of the plurality of horizon surfaces or surface patches or excluded from the subset of the plurality of horizon surfaces or surface patches based on their respective scores; iteratively repeating the selecting and calculating steps until a stopping condition is reached and the subset of the plurality of horizon surfaces or surface patches is determined; and performing interpretation on the subset of the plurality of horizon surfaces or surfaces patches.

The disclosure provides methods and systems for characterizing depositional features. The methods and systems include accessing data encoding seismic waves as seismic traces reflected from cells at various locations within a particular stratum in response to a seismic source. The cells are classified into multiple non-overlapping groups according to the amplitude values or other seismic attributes of the seismic waves reflected from the various locations within the particular stratum. One or more subgroups of adjacent cells are identified. A subgroup area metric is calculated for each subgroup of cells by combining individual area metrics from adjacent cells in a given subgroup and subsequently assigning the calculated subgroup area metric to each cell of the given subgroup. One or more depositional features within the stratum are characterized based at least in part on the variation map based on the subgroup area metric of each cell.