Hypersurface curvatures of geological features. Issue 1 (23rd January 2023)

• Record Type: Journal Article
• Authors: Ravve, Igor; Tertois, Anne-Laure; de Ribet, Bruno; Koren, Zvi
• Is Part Of: Geophysical journal international. Volume 234:Issue 1(2023)
• Page Start: 88
• Page End: 123
• Publication Date: 2023-01-23
• Subjects: Image processing -- Numerical modelling -- Spatial analysis -- Fractures, faults, and high strain deformation zones -- Marine geosciences and applied geophysics; Geophysics -- Periodicals; 550
• DOI: 10.1093/gji/ggad021
• Languages: English
• ISSNs: 0956-540X

Abstract

SUMMARY: Reflector-normal angles and reflector-curvature parameters are the principal geometric attributes used in seismic interpretation for characterizing the orientations and shapes, respectively, of geological reflecting surfaces. Commonly, the input data set for their computation consists of fine 3-D grids of scalar fields representing either the seismic-driven reflectivities (e.g. amplitudes of 3-D seismic migrated volumes) or model-driven reflectivities, computed, for example, from the derived elastic impedance parameters. Conventionally, the computation of curvature parameters at each gridpoint is based on analysing the local change in the inline/crossline dips, considering the potential existence of a local quadratic reflecting surface in the vicinity of that point. This assumption breaks down for subsurface points in the vicinity of either complex reflecting surfaces (e.g. brittle/rough/tilted synclines/anticlines, ridges/troughs and saddles) and/or sharp, discontinuous geological features (e.g. fault edges/tips, pinch-outs, fracture systems, channels and small geobodies), where the values of the computed curvature become extremely high. However, while these high values can indicate the existence of non-reflecting objects, they do not deliver their specific geometric characteristics. In this study, we present a novel method that better characterizes the shapes of these complex geological features by extending the assumption of local surfaces (2-D surfaces in 3-D space) into local hypersurfaces (3-D hypersurfaces in 4-D space), with their corresponding (three rather than two) principal (and effective) curvature parameters. We demonstrate the advantages of our method by comparing the conventional dip-based surface curvature parameters with the hypersurface curvature parameters, using a synthetic model/image with different types and shapes of geological features and a seismic image of real data containing a complex fault and hidden buried channels.