Systems and methods of placing wells in a hydrocarbon field based on seismic attributes and quality indicators associated with a subterranean formation of the hydrocarbon field can include receiving seismic attributes representing the subterranean formation and seismic data quality indicators. A cutoff is generated for each seismic attribute and seismic data quality indicator. A weight is assigned to each seismic attribute and seismic data quality indicator. The weighted seismic attributes and data quality indicators are aggregated for each location in the hydrocarbon field. A risk ranking is assigned based on the weighted seismic attributes and data quality indicators associated with each location in the hydrocarbon field based on the cutoffs. A map is generated with each location on the surface of the subterranean formation color-coded based on its assigned risk ranking.

A method for generating a seismic image representing a subsurface includes receiving seismic data for the subsurface formation, including receiver wavelet data and source wavelet data. Source wavefield data are generated based on a forward modeling of the source wavelet data. Receiver wavefield data are generated that compensate for distortions in the seismic data by: applying a dispersion-only model to the receiver wavelet data to generate a first reconstructed back-propagated receiver wavefield portion, applying a dissipation-only model to the receiver wavelet data to generate a second reconstructed back-propagated receiver wavefield portion, and combining the first back-propagated receiver wavefield portion and the second back-propagated receiver wavefield portion into the receiver wavefield data. The method includes applying an imaging condition to the receiver wavefield data and the source wavefield data and generating, based on applying the imaging condition, visco-acoustic reverse time migration (VARTM) result data.

A method for training a backpropagation-enabled segmentation process is used for identifying an occurrence of a sub-surface feature. A multi-dimensional seismic data set with an input dimension of at least two is inputted into a backpropagation-enabled process. A prediction of the occurrence of the subsurface feature has a prediction dimension of at least 1 and is at least 1 dimension less than the input dimension.

The disclosure presents processes to automatically generate one or more set of fault segments from a fault plane pointset. The processes can identify a predominant direction and derive a set of fault segments from the fault plane pointset, where the fault segments are generated by using slices of data from the fault plane pointset that are perpendicular to the predominant direction. For each slice of data, the fault segments can be analyzed with neighboring fault segments to determine if they are overlapping. Fault segments that block or overlap other fault segments can be assigned to a different subset of fault segments from the underlying fault segments. Gaps in the fault plane pointset, and the resulting set of fault segments, can be filled in by merging neighboring fault segments above and below the gap if the neighboring fault segments satisfy a criteria for filling the gap.

Disclosed is a method for analyzing and predicting a main fracture orientation of a mining face based on microseismic monitoring, including: collecting microseismic data generated by a coal rock burst; carrying out a hierarchical clustering on the microseismic data to obtain target hypocenter groups, of which the target hypocenter groups comprise several types of hypocenters; acquiring focal mechanism solutions of all the target hypocenter groups in the target hypocenter group, and acquiring a hypocenter azimuth and a hypocenter dip based on the focal mechanism solutions; and carrying out the hierarchical clustering on a hypocenter location, the hypocenter azimuth and the hypocenter dip, and predicting the main fracture orientation of the mining face.

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 three-dimensional velocity geological modeling with structures and velocities randomly arranged, including determining base points in three-dimensional space, building equation according to the base points to determine planar layered model, complicating a tilt layer of planar layered model, and building a fold layer model of a surface in three-dimensional space; building three-dimensional fault folded model based on the three-dimensional surface fold layer model combined with a fault plane of a random reference point and displacement of each point in a global coordinate system; building a velocity model containing a salt body based on the three-dimensional fault folded model, and simulating salt body intrusion in a geological body of a certain depth; and performing a random velocity amplitude to realize three-dimensional velocity modeling according to the layered type which has been set and according to the set velocity range and the velocity difference range between each layer of geology.

A full probability-based seismic risk analysis method for a tunnel under fault dislocation comprises: evaluating a magnitude-frequency relationship of a fault; obtaining a probabilistic seismic risk curve of a fault dislocation; calculating a series of bending moments of a tunnel lining under different fault dislocations; obtaining a series of damage index values R.sub.M of the tunnel; obtaining a vulnerability model of the tunnel damaged by fault dislocation; calculating a probabilistic risk that the tunnel crossing the fault is damaged due to the dislocation of the active fault; obtaining a probability P that the damage state is equal to or higher than a certain damage state within a specified period; and using the results to guide the assessment of the seismic risk of the tunnel crossing the fault. Modeling and analysis can be performed according to the actual situation of the tunnel crossing the fault with different factors.

Systems and methods for interpreting and visualizing multi-Z horizons from seismic data are disclosed. A two-dimensional (2D) representation of seismic data is displayed via a graphical user interface (GUI). User input is received via the GUI for interpreting a multi-Z horizon within a portion of the displayed 2D representation. The user's input is tracked relative to displayed 2D representation within the GUI. Based on the tracking, each of a plurality of surfaces for the multi-Z horizon is determined. At least one intersection point between the multi-Z horizon surfaces is identified. A depth position for each surface relative to other surfaces is determined. The 2D representation of the seismic data is dynamically updated to include visual indications for the plurality of surfaces and the intersection point(s), based on the depth position of each surface, where the visual indications use different visualization styles to represent the surfaces and intersection point(s).

A fault indicator calculator, a method for determining a fault indicator, and a fault indicator calculating system are disclosed herein. One embodiment of a fault indicator calculator includes: 1) an interface configured to receive seismic data, and 2) a processor configured to scan a manifold-shaped operator through said seismic data at a range of dips and azimuths and calculate fault throws at various orientations of said dips and azimuths independent of determining other fault indicators.