Abstract: Skeletonization and Empirical-Mode Decomposition of 2-D potential-field images

The volumes of modern gridded data collected by the geophysical exploration industry are often large and can benefit from methods for image decomposition, pattern analysis, and interpretation. Automated pattern-recognition methods can be useful for both seismic and potential-field images. In both of these cases, it is important to extract quantitative attributes useful for further data analysis and inversion for the subsurface structure. In this paper, we investigate one of such pattern-analysis approaches called “skeletonization”.

Although using different physical fields and models, many types of 2-D geophysical images, such as seismic sections and slices, gravity and magnetic maps, possess a number of similar geometrical features. These features can be expressed by linear continuity, branching, amplitudes, widths, polarities, curvatures, orientations and/or other attributes and can be subdivided into “background trends” on top of which some kinds of “anomalies” or “wavelets” can be recognized. The types of spatial dimensions of the images may also vary, ranging from the usual distances, elevations and depths to travel times and travel-time lags. Automatic identification of such spatially-connected wavelets and measurement of their parameters is the general objective of skeletonization.

This paper describes a new, azimuthally-uniform skeletonization approach to 2-D potential field data. The majority of the material in the following sections was originally published as Paper A-3 in Summary of Investigations 2012, Volume 1, Saskatchewan Geological Survey Miscellaneous Report 2012-4.1 (Gao and Morozov, 2012). This material is reproduced here, with modifications, with permission from the Saskatchewan Geological Survey.


Le Gao received her B.S in Geology and M.S in Applied Geophysics from Northwest University in China, and a PhD in Geophysics (2017) from the University of Saskatchewan, Canada. Her research focused on 3D surface land seismic data analysis, including time-lapse monitoring, AVO, seismic inversion, receiver functions, and pattern recognition in seismic and potential-field images. Le is currently at CGG, Calgary, AB and has experience in velocity model building, pre-stack time and depth migration. 

Igor B. Morozov received his Ph.D. (1985) in theoretical physics from Moscow State University (Russia). He studies seismology at scales ranging from free oscillations, earthquake codas and nuclear explosions to deep crustal profiling, exploration, and highresolution. Recent research interests include seismic imaging and attributes, potential-field imaging, time-lapse seismic monitoring, seismic attenuation, theoretical rock physics, and integrated geophysical software. 


Share your thoughts with us! Comments are moderated. Spam, inappropriate or abusive comments will be removed.