REDUCING NON-UNIQUENESS IN SATELLITE GRAVITY INVERSION USING SEISMIC TOMOGRAPHY AND 3D OBJECT ORIENTED IMAGE ANALYSIS TECHNIQUES ISLAM FADEL MARK VAN DER MEIJDE NORMAN KERLE
Research Topic: Reducing the non-uniqueness of satellite gravity inversion using constraints based on 3D seismic tomography data applying 3D Object Oriented Image Analysis (OOA).
3D OOA (Bio-medical) Recent Image analysis technique. Ruleset = Set of classification criteria Extract 3D objects from 3D image stacks. Schӧnmeyer et al. (2006)
3D OOA (Geosciences) Extract 3D objects from LiDAR point cloud. Example by ecognition community: http://community.ecognition.com/home/ecognition-labs/3d-lidar-point-cloud-analysis/?searchterm=3d
How did the non-uniqueness of satellite gravity been recently reduced? The non-uniqueness has been reduced using seismic data. 3D Seismic Tomography Cammarano et al., 2011 Conversion functions Velocity Density Mishra et al. (2012) Ebbing et al. (2001) Subsurface Interpretation of layers and objects then forward modeling using estimated density values. Additional uncertainties due to the dependency on other physical parameters: Temperature Pressure Composition Trial and error. Usually 2D over profiles. NEW 3D structures from 3D seismic tomography data constrain the inversion of satellite gravity.
The Methodology The problem of the non-uniqueness has been divided into 2 sub-problems: The non-uniqueness The shape of the subsurface objects The values of the density contrasts of these objects?? Reduced using 3D OOA of 3D seismic tomography model. Reduced using Object-based inversion and Forward Modeling.
The Methodology 3D seismic tomography model (4) Object-based inversion (1) 3D OOA (3) Reconstruct the objects in IGMAS+ (5) Forward modeling (2) Classify the subsurface objects (6) Evaluate
Study Area From 3D tomography model and literature the study area mainly composed of : Tanzania Craton. Eastern and Western rift branches. Low velocity anomaly (start ~ at 250 km depth). (Adams et al., 2012) (Adams et al., 2012)
1-3D OOA results Results 2- Objects-Reconstruction 3- Object-based inversion
3D OOA Results Use of 3D histogram and visual interpretation. Homogeneity and discontinuity characteristics to extract 3D objects. 3D histogram upper mantle part 40-500 km
Object-Based inversion (Upper Mantle) Object-Based Inversion Measured signal EIGEN6C2 Craton Rift Rift 70 % correlation Boundary Calculated signal from the upper mantle part. Low velocity Boundary 10000 9500 9000 km Deep High velocity Construction Process km 0.0 100 200 300 400 0.0 500 km 1000
Objects-Based Inversion (Crustal Part) Object-Based Inversion MOHO Measured signal EIGEN6C2 95 % correlation Crust Calculated signal from the full model. Moho 0 00 10 Upper Mantle 00 95 m k Add layers represent the crustal part from the 3D tomography model. 0.0 km 10 00 100 0 50 m k 200 00 90 300 400 0.0 The Moho (Tugume et al. (in press))
Objects-Based Inversion Calculated signal Measured signal EIGEN6C2 The difference between the measured and the calculated signal from the full model.
Upper Mantle Crust Density Contrast Values Object Density Contrast (t/m3) Unit 1 (U1) -0.215 Unit 2 (U2) -0.243 Unit 3 (U3) -0.277 Unit 4 (U4) -0.312 Unit 5 (U5) -0.374 MOHO MOHO 5.154 Craton (C) 0.008 Shallow High Velocity (SHV) 0.011 Boundary Shallow High Velocity (BSHV) 0.005 Rift (R) -0.002 Boundary Shallow Low Velocity (BSLV) -0.015 Low Velocity Zone (LV) -0.053 Inner Low Velocity (ILV) -0.002 Boundary Deep Low Velocity (BDLV) -0.067 Deep High Velocity (DHV) 0.173 Reference 0 Correlation 0.949 Standard deviation 11
Conclusion and Recommendations A new methodology was developed and tested to reduce the non-uniqueness of the gravity modeling. 3D OOA can extract the 3D objects from 3D geophysical data. The methodology was able to reduce the non-uniqueness of the shape of the subsurface objects. However, the density contrast estimation needs to be constrained. The object based inversion approach is very promising and fast. However, IGMAS+ needs to be developed to allow the user to put boundary conditions for density contrast estimations.
THANK YOU For more details: Fadel, I., Kerle, N., van der Meijde, M., 2014. 3-D object-oriented image analysis of geophysical data. Geophysical Journal International, 10.1093/gji/ggu139. Fadel, I., et al., 3D object-oriented image analysis in 3D geophysical modelling: Analysing the central part of the East African Rift System. Int. J. Appl. Earth Observ. Geoinf. (2013), http://dx.doi.org/10.1016/j.jag.2013.11.004. References Adams, A., Nyblade, A., Weeraratne, D., 2012. Upper mantle shear wave velocity structure beneath the East African plateau: evidence for a deep, plateauwide low velocity anomaly. Geophysical Journal International 189 (1), 519 123 142. Götze, H., Lahmeyer, B., 1988. Application of three-dimensional interactive modeling in gravity and magnetics. Geophysics 53 (8), 1096 1108. Schmidt, S., Götze, H., Fichler, Ch.,A.M., 2010. IGMAS+: a new 3D gravity, FTG and magnetic modelling software. Extended abstract. Geoinformatik, 57 63. Schmidt, S., Plonka, C., Götze, H., Lahmeyer, B., 2011. Hybrid modelling of gravity, gravity gradients and magnetic fields. Geophysical Prospecting 59 (6, SI), 1046 1051. Schӧenmeyer, R., Prvulovic, D., Rotarska-Jagiela, A., Haenschel, C., Linden, D. E. J., 2006. Automated segmentation of lateral ventricles from human and primate magnetic resonance images using cognition network technology. Magnetic Resonance Imaging 24 (10), 1377-1387. Tugume, F., Nyblade, A., Julia, J., van der Meijde, M., 2013. Precambrian crustal structure in Africa and Arabia: evidence lacking for secular variation. Tectonophysics (in press).