The present invention proposes a method to represent seismic fault damage zones and fracture density in the geological models of reservoirs in a simple, agile and automated way, so that it can be easily replicated by geologists in any production design. It was developed as a group of workflows, inserted in the commercial software Petrel, widely used in the company for the 3D numerical modeling of reservoirs.

A method can include receiving data associated with a geologic environment; based on at least a portion of the data, estimating relationships for multiple properties of the geologic environment; and based at least in part on the relationships, performing simultaneous joint inversion for at least one property of the geologic environment.

A method of predicting three-dimensional fracture geometry in a subterranean region of interest is disclosed. The method includes obtaining a strain tensor for the subterranean region of interest, calculating a set of principal strain components from the strain tensor, and determining a strain cyclide from the set of principal strain components. The method further includes calculating a set of quadrimodal fault normal vectors from the strain cyclide and determining an in-plane shear strain magnitude and a shear strain orientation from the set of quadrimodal fault normal vectors.

The approaches presently disclosed provide for fault-interpretation in a seismic volume with computer assistance, allowing automatic or semi-automatic determination of a fault surface and associated displacement across the fault. The present fault interpretation approach uses pattern matching algorithms and does not require prior interpretation of the stratigraphic horizons. In certain implementations the fault interpretation approach estimates the 3D fault surface as part of a joint fault surface location and displacement optimization process.

Methods and systems, including computer programs encoded on a computer storage medium can be used to preserve edges while performing Gaussian grid smoothing of noise components in subsurface grids to generate geological maps. A subsurface grid is generated from data indicating properties of subsurface formations. A weighting grid is generated by: i) receiving seismic data representing the subsurface formations; ii) generating seismic attributes associated with discontinuities in the subsurface formations; and iii) assigning a particular weight value to weighting grid points that the seismic attributes associated with discontinuities in the subsurface formations indicate the presence of a discontinuity. The subsurface grid is processed by iteratively computing local averages of grid points in the subsurface grid using a compact Gaussian filter weighted by values in the weighting grid. A geological map of subsurface formations is generated based on the filtered subsurface grid.

A system can receive fault likelihood data about a subterranean environment and apply a binary mask filter using a tuning parameter to convert the fault likelihood data to binary distribution data having a plurality of pixels arranged in a plurality of profiles in at least two directions. The system can perform, for each profile of the plurality of profiles, fault skeletonization on the binary distribution data to form fault skeletonization data with pixels connected that represent part of a fracture. The system can convert the fault skeletonization data to seismic volume data and combine and filter the seismic volume data in the at least two directions to form combined seismic volume data. The system can output the combined seismic volume data as an image for use in detecting objects to plan a wellbore operation.

A method for determining an intersection between a polygon representing a boundary of a surface in an earth formation and a plane includes: receiving a polygon representing a boundary of a surface in an earth formation, the polygon having a series of straight segments with a point at each end of each of the segments; overlaying a cutting grid having grid planes over the polygon; identifying a specific pattern of two adjacent segments in the polygon by proceeding from a first segment to a second segment in a selected rotational direction; matching the specific pattern to a reference pattern; categorizing the point between the two adjacent segments as an intersection point or as a non-intersection point based on the reference pattern; the iterating the identifying, matching, and categorizing for each pair of adjacent segments in the polygon such that each point between adjacent segments in the polygon is categorized.

The present disclosure describes methods and systems, including computer-implemented methods, computer program products, and computer systems, for smoothing seismic data. One computer-implemented method includes obtaining, by a hardware data processing apparatus, a plurality of seismic data samples; forming, by the hardware data processing apparatus, guiding vectors using the plurality of seismic data samples and a plurality of guiding structure attributes; generating, by the hardware data processing apparatus, a structure guided directional weighted vector filter using the guiding vectors and a plurality of weighting factors; filtering, by the hardware data processing apparatus, the seismic data samples using the structure guided directional weighted vector filter to generate smoothed seismic data; and initiating output of the smoothed seismic data.

The present disclosure describes computer-implemented methods, computer-program products, and computer systems, for providing parameters for successful automated fault patch extraction. Attributes are selected for annotating images generated from, and for interpretation of, seismic amplitude volume. Images are generated from layers of a seismic cube, each generated using a different attribute of the plural attributes. The plural images are blended using customized palettes and initial parameters to create a blended image illuminating discontinuities in the layers. Optimal parameters are iteratively determined for automatic derivation of fault discontinuities on an interpreter-selected edge-enhanced sub-volume. The iterations are controlled and terminated based on interpreter inputs. The optimal extraction parameters are applied to an entire edge-enhanced volume. Important extracted fault discontinuities are isolated using commercial filtering tools. Extracted fault patches are refined based on received manual interpretation. Patch results are converted to traditional fault objects for further interpretation and refinement.

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.