130 likes | 272 Views
Gradient based technique for electromagnetic layered earth model data inversion. Claudio Patriarca, Andrea Di Matteo and Evert Slob. 2011 IEEE Internation Geoscience and Remote Sensing Symposium, Vancouver BC, July 2011. Local or Global inversion. Inversion Schemes. Gradient descent.
E N D
Gradient based technique for electromagnetic layered earth model data inversion Claudio Patriarca, Andrea Di Matteo and Evert Slob 2011 IEEE Internation Geoscience and Remote Sensing Symposium, Vancouver BC, July 2011
Local or Global inversion Inversion Schemes Gradient descent Global Multilevel Coordinate search
Full-waveform inversion EM Inversion Schemes – non linear • Size of parameter space – nr. layers • Number of local minima - nonlinearity • Parameter correlations – limit sensitivity • UWB – many frequency points
Towards local inversion Jacobian calculation Objective Function † Jacobian
Jacobian calculation Numerical vs Explicit gradient • Numerically evaluated • Explicitly given
Local Full-waveform inversion Numerical test case #1: smooth model
Local Full-waveform inversion Numerical test case #1 result
Full-waveform inversion Real data test case • Large datasets • UWB data (SFCW) Marble floor element in the Ancient theatre of Megalopolis, Greece
Local Full-waveform inversion Real data test case
Local Full-waveform inversion Real data test case results
Local Full-waveform inversion Numerical test case #2: piecewise layered half-space
Objective Function topography Numerical test case #2: piecewise layered half-space Global Minimum ε2 h2 ε1
Conclusion • Local inversion advantages in moving downhill; problems with local minima • Comparison Global and Local search numerical and real data: success of local methods in specialized applications • Use explicit gradient • Deep interfaces? • Convenient implementing Hessian?