Patterns in Geophysical Data and Models Jens Tronicke Angewandte Geophysik Institut für Geowissenschaften Universität Potsdam jens@geo.uni-potsdam.de
Near-surface geophysics Using geophysical tools to explore the shallow subsurface Scale of interest: Tens of centimeters up to tens of meters Common techniques: -Magnetics - Seismics (reflection and refraction) - Electrical resistivity and EM induction methods - Georadar (GPR) - Borehole and logging techniques Typical applications are from the fields of - Geology - Engineering geology and geotechnics - Sedimentology - Hydrogeology - Soil sciences - Environmental sciences - Archaeology - Agriculture - Civil engineering - Research topics are directed towards developing and improving geophysical techniques for specific applications
Patterns in Geophysical Data & Models What is a pattern? - Opposite of chaos (Watanabe, 1985) - A regular repetition of values in an image (Sheriff, 2002) General flow of a geophysical experiment Data gathering (employing magnetics, seismics, electrics, georadar, borehole techniques, ) Processing (editing, analyzing, filtering, inversion, In part specially developed techniques) Analyzing & enhancing specific features in the data space Examples from seismic & georadar: Residual statics (cross-correlation) Velocity analysis (coherency analysis of CMP gathers) NMO stacking Coherency & dip filtering
Patterns in Geophysical Data & Models What is a pattern? - Opposite of chaos (Watanabe, 1985) - A regular repetition of values in an image (Sheriff, 2002) General flow of a geophysical experiment Data gathering (employing magnetics, seismics, electrics, georadar, borehole techniques, ) Processing (editing, analyzing, filtering, inversion, In part specially developed techniques) Interpretation (generating earth models based on one or several geophysical techniques) Analyzing, extracting patterns & features in the model space, e.g.: Visualizing specific structures (faults, sedimentological units, archaeological features, ) Integrating all available models & data (including non-geophysical information, e.g., core data, geological background, geotechnical or hydrological data, )
Overview Selected examples from geophysical data interpretation 3-D georadar data: Attribute analysis to extract fault related features Multi-technique data sets: Integrated earth models based on cluster analysis
Overview Selected examples from geophysical data interpretation 3-D georadar data: Attribute analysis to extract fault related features Multi-technique data sets: Integrated earth models based on cluster analysis
Maleme Fault Zone, New Zealand Goal of the study: Mapping shallow geometry of the fault zone and stratigraphic details within near surface sediments using 2-D and 3-D georadar 3-D data set (100 MHz antennas) covering an area of ~ 20 x 70 m Central cut through migrated 3-D volume (x = 14 m) 0 0 (Coates, 2002) 100 3 Time [ns] 200 6 Depth [m] 300 9 NW SE -30-20 -10 0 10 20 30 y [m]
Attribute analysis (3-D data) Post-processing sequence to visualize near vertical features Migrated data cube Instantaneous phase: cos ( φ ) φ Emphasizing continuity/discontinuity Local gradients: d x and d y Slope: S = + 2 2 d x d y Emphasizing local changes (including information on directions)
Maleme Fault Zone: 3-D data Central cut through 3-D volume (x = 14 m) 0 0 100 3 Time [ns] 200 6 Depth [m] 300 9 NW SE -30-20 -10 0 10 20 30 y [m]
Maleme Fault Zone: 3-D data Central cut through 3-D volume (x = 14 m) 0 0 100 3 195 ns Time [ns] 200 6 Depth [m] 300 9 NW SE -30-20 -10 0 10 20 30 y [m] Reflection A Reflection B
Maleme Fault Zone: 3-D Daten 195 ns (~5.85 m) NW SE
Maleme Fault Zone: 3-D Daten 195 ns (~5.85 m) NW SE
Maleme Fault Zone: 3-D Daten 195 ns (~5.85 m) NW SE
Maleme Fault Zone: 3-D Daten 195 ns (~5.85 m) NW SE
Maleme Fault Zone: 3-D Daten 195 ns (~5.85 m) NW SE
Maleme Fault Zone: 3-D data Main discontinuities picked in 3-D slope cube Depth [m] 4.5 5.0 5.5 6.0 6.5 7.0 7.5 8.0
Overview Selected examples from geophysical data interpretation 3-D georadar data: Attribute analysis to extract fault related features Multi-technique data sets: Integrated earth models based on cluster analysis
Geophysical tools for aquifer characterization For example: contaminant site characterization Groundwater table Unsaturated sediments Depth to bedrock Saturated sediments Bedrock Borders of the aquifer Internal structures Hydraulic relevant parameters (e.g., porosity, permeability)
Crosshole georadar and seismic tomography Exploring the inter-borehole plane - Traveltime of the signals velocity v - Amplitudes of the signals attenuation α Scale of such experiments: several meters to tens of meters Resolution: ~ 1 m Sensitive to hydrological relevant parameters: water content, porosity, Often ideal complement to the more conventional hydrological field techniques n sources m receivers
Synthetic example How to generate an integrated model from different geophysical models?
Fuzzy c-means (FCM) cluster analysis n-dimensional model space FCM cluster analysis m ij : membership of data point x j to cluster i Parameter 2 Parameter 1
Fuzzy c-means (FCM) cluster analysis n-dimensional model space FCM cluster analysis m ij : membership of data point x j to cluster i Parameter 2 Parameter 1 Defuzzification hard cluster separation (Colour saturation proportional to membership) z [m] Zoned multi-parameter model including information regarding its reliability x [m]
Fuzzy c-means (FCM) cluster analysis n-dimensional model space FCM cluster analysis m ij : membership of data point x j to cluster i z [m] x [m] Defuzzification hard cluster separation (Colour saturation proportional to membership) Multi-prameter model: p i p j = c i= 1 m ij p : mean of p for cluster i i Target parameter p (e.g. porosity-logs) Zoned multi-parameter model including information regarding its reliability Spatial model of the petrophysical target parameter
Synthetic example
Results FCM cluster analysis Colour saturation reflects reliability of the zonation
Reconstructed porosity distributions Resampled input model Three-cluster solution
Results cluster analysis Real data example: Boise Hydrogeophysical Research Site (BHRS) Velocity Attenuation Clustered model
Analyzing reliability of the zonation Fuzzy c-means (after defuzzification) Membership function for cluster 1 ( low porosity cluster )
Petrophysical parameter reconstruction Fuzzy c-means (after defuzzification) Reconstructed porosity distribution
Conclusions Pattern & feature recognition and related methods are fundamental steps in the geophysical workflow Data processing Interpretation Examples illustrated the potential to - extract specific features in a data set - integrate different data sets in a quantitative fashion
Acknowledgements Deutsche Forschungsgemeinschaft (DFG) Schweizerischer Nationalfonds (SNF) P. Villamor (GNS, New Zealand) W. Barrash & M. Knoll (Boise State University) K. Holliger & H. Paasche (ETH Zürich)