Project

Back to overview

New Developments of efficient geophysical tomography algortithms for exploring near-surface structures

English title New Developments of efficient geophysical tomography algortithms for exploring near-surface structures
Applicant Maurer Hansruedi
Number 132745
Funding scheme Project funding
Research institution Institut für Geophysik ETH Zürich
Institution of higher education ETH Zurich - ETHZ
Main discipline Geophysics
Start/End 01.05.2011 - 30.04.2014
Approved amount 166'448.00
Show all

All Disciplines (3)

Discipline
Geophysics
Mathematics
Other disciplines of Physics

Keywords (5)

Seismics; adpative methods; Wavelets; Geophysics; Modelling

Lay Summary (English)

Lead
Lay summary
Inverting complete waveforms of seismic data has the potential to significantly increase our knowledge about the interior of the earth, including the uppermost 100 m, which are of vital interest for numerous problems of high societal relevance (e.g., groundwater resources, storage of dangerous waste, natural hazards and major building projects). However, the large computational expenses together with problematic system effects, such as the generally unknown source characteristics and variable receiver coupling, have precluded application of this potentially very powerful technique on a routine basis. Building on the latest developments in the applicant's research group in the areas of numerical modelling, inversion technology and experimental design, novel waveform inversion algorithms and data acquisition strategies will be developed. In particular, an adaptive-wavelet modelling technique, recently implemented in our new 3D geoelectrical code, will be modified such that the frequency-domain acoustic response of a 2D model can be computed. If this task will be successful, further extensions to visco-acoustic, elastic, visco-elastic and ground penetrating radar (GPR) problems will be performed. Particular emphasis will be put on exploiting common features of the different problems. For example, the GPR problem is formally equivalent to the visco-acoustic case. Likewise, the equations of motion that describe the (visco-)elastic problems can be decomposed into Helmholtz-type scalar equations, on which (visco-)acoustic modelling is based. By sharing as many common elements as possible in the modelling codes, it will be relatively straightforward to apply further developments, such as 2.5D modelling, simultaneously to all types of problems. Besides improving numerical forward modelling, critical issues of waveform inversions, such as pronounced non-linearities and systematic effects caused by variable receiver coupling will be studied and appropriate strategies to account for these problems will be devised. All the theoretical developments will be verified with a 3D field experiment, where comprehensive seismic and GPR data sets will be acquired and analysed.
Direct link to Lay Summary Last update: 21.02.2013

Responsible applicant and co-applicants

Employees

Publications

Publication
A hybrid boundary element-finite element approach to modeling plane wave 3D electromagnetic induction responses in the Earth
Ren Zhengyong, Kalscheuer Thomas, Greenhalgh Stewart Alan, Maurer Hansruedi R. (2014), A hybrid boundary element-finite element approach to modeling plane wave 3D electromagnetic induction responses in the Earth, in Journal of Computational Physics, 258, 705-717.
A critical appraisal of asymptotic 3D-to-2D data transformation in full-waveform seismic crosshole tomography
Auer Ludwig, Nuber Andre Marc, Greenhalgh Stewart Alan, Maurer Hansruedi, Marelli Stefano (2013), A critical appraisal of asymptotic 3D-to-2D data transformation in full-waveform seismic crosshole tomography, in GEOPHYSICS, 78(6), 235-247.
A goal-oriented adaptive finite-element approach for plane wave 3-D electromagnetic modelling
Ren Zhengyong, Kalscheuer Thomas, Greenhalgh Stewart Alan, Maurer Hansruedi R. (2013), A goal-oriented adaptive finite-element approach for plane wave 3-D electromagnetic modelling, in Geophysical Journal International, 194(2), 700-718.
Boundary element solutions for broad-band 3-D geo-electromagnetic problems accelerated by an adaptive multilevel fast multipole method
Ren Zhengyong, Kalscheuer Thomas, Greenhalgh Stewart Alan, Maurer Hansruedi R. (2013), Boundary element solutions for broad-band 3-D geo-electromagnetic problems accelerated by an adaptive multilevel fast multipole method, in Geophysical Journal International, 192(2), 473-499.
Boundary element solutions for broad-band 3-D geo-electromagnetic problems accelerated by an adaptive multilevel fast multipole method
Ren Zhengyong, Kalscheuer Thomas, Greenhalgh Stewart, Maurer Hansruedi (2013), Boundary element solutions for broad-band 3-D geo-electromagnetic problems accelerated by an adaptive multilevel fast multipole method, in GEOPHYSICAL JOURNAL INTERNATIONAL, 192(2), 473-499.
Investigating the dynamics of an alpine glacier using probabilistic icequake locations: Triftgletscher, switzerland
Dalban Canassy Pierre, Walter Fabian, Husen Stephan, Maurer H., Faillettaz Jérôme, Farinotti Daniel (2013), Investigating the dynamics of an alpine glacier using probabilistic icequake locations: Triftgletscher, switzerland, in Journal of Geophysical Research F: Earth Surface, 118(4), 2003-2018.
2.5-D frequency-domain seismic wave modeling in heterogeneous, anisotropic media using a Gaussian quadrature grid technique
Zhou B, Greenhalgh S, Maurer H (2012), 2.5-D frequency-domain seismic wave modeling in heterogeneous, anisotropic media using a Gaussian quadrature grid technique, in COMPUTERS & GEOSCIENCES, 39, 18-33.
3-D electrical resistivity tomography using adaptive wavelet parameter grids
Plattner A, Maurer HR, Vorloeper J, Blome M (2012), 3-D electrical resistivity tomography using adaptive wavelet parameter grids, in GEOPHYSICAL JOURNAL INTERNATIONAL, 189(1), 317-330.
Exploitation of data-information content in elastic-waveform inversions
Manukyan E, Latzel S, Maurer H, Marelli S, Greenhalgh SA (2012), Exploitation of data-information content in elastic-waveform inversions, in GEOPHYSICS, 77(2), 105-115.
GPR Full-Waveform Sensitivity and Resolution Analysis Using an FDTD Adjoint Method
Meles GA, Greenhalgh SA, Green AG, Maurer H, Van der Kruk J (2012), GPR Full-Waveform Sensitivity and Resolution Analysis Using an FDTD Adjoint Method, in IEEE TRANSACTIONS ON GEOSCIENCE AND REMOTE SENSING, 50(5), 1881-1896.
Receiver-coupling effects in seismic waveform inversions
Maurer H, Greenhalgh SA, Manukyan E, Marelli S, Green AG (2012), Receiver-coupling effects in seismic waveform inversions, in GEOPHYSICS, 77(1), 57-63.
Seismic monitoring of radioactive waste repositories
Manukyan E, Maurer H, Marelli S, Greenhalgh SA, Green AG (2012), Seismic monitoring of radioactive waste repositories, in Geophysics, 77(6), EN73-EN83.
Validity of the acoustic approximation in full-waveform seismic crosshole tomography
Marelli S., Maurer H.R., Manukyan E. (2012), Validity of the acoustic approximation in full-waveform seismic crosshole tomography, in Geophysics, 77(3), R129-R139.
A new model for estimating subsurface ice content based on combined electrical and seismic data sets
Hauck C, Bottcher M, Maurer H (2011), A new model for estimating subsurface ice content based on combined electrical and seismic data sets, in CRYOSPHERE, 5(2), 453-468.

Scientific events



Self-organised

Title Date Place
Workshop on seismic waveform inversion 15.02.2013 Karlsruhe, Germany

Associated projects

Number Title Start Funding scheme
121707 New Developments of efficient geophysical tomography algorithms for exploring near-surface structures 01.11.2008 Project funding
153089 New Developments of efficient geophysical tomography algortithms for exploring near-surface structures 01.05.2014 Project funding

Abstract

Inverting complete waveforms of seismic data has the potential to significantly increase our knowledge about the interior of the earth, including the uppermost 100 m, which are of vital interest for numerous problems of high societal relevance (e.g., groundwater resources, storage of dangerous waste, natural hazards and major building projects). However, the large computational expenses together with problematic system effects, such as the generally unknown source characteristics and variable receiver coupling, have precluded application of this potentially very powerful technique on a routine basis. Building on the latest developments in the applicant’s research group in the areas of numerical modelling, inversion technology and experimental design, novel waveform inversion algorithms and data acquisition strategies will be developed. The project involves two subprojects: one devoted to numerical forward modelling (project A) and one devoted to inversion technology and experimental design problems (project B). In project A, an adaptive-wavelet modelling technique, which we originally developed for 3D geoelectrical problems, will be significantly modified such that the frequency-domain acoustic response of a 2D model can be computed. The performance of this algorithm will be compared with state-of-the-art adaptive finite-element modelling codes. Based on these comparisons, the most suitable of the two techniques will be further developed and adapted to tackle visco-acoustic, elastic, visco-elastic and ground penetrating radar (GPR) problems. Particular emphasis will be put on exploiting common features of the different problems. For example, the GPR problem is formally equivalent to the visco-acoustic case. Likewise, the equations of motion that describe the (visco)-elastic problems can be decomposed into Helmholtz-type scalar equations, on which (visco-)acoustic modelling is based. By sharing as many common elements as possible in the modelling codes, it should be more straightforward to apply further developments, such as 2.5D modelling, simultaneously to all types of problems. In project B, optimised inversion model parameterisations will be investigated. Again, this task will benefit from previous research on geoelectrical problems conducted in the applicant’s research group. Moreover, critical issues of waveform inversions, such as pronounced non-linearities and systematic effects caused by variable receiver coupling, will be studied and appropriate strategies to account for these problems will be devised. In a second phase of project B, statistical experimental design will be used to determine optimised field layouts for seismic and GPR waveform inversion experiments. In addition, an information content analysis will help identify particularly important (i.e., information-rich) portions of the seismic and GPR waveforms. The new algorithms developed for projects A and B will be tested on comprehensive 3D seismic and GPR data sets to be acquired and analysed as components of this application.
-