Constrained optimization framework for joint inversion of geophysical data sets

Many experimental techniques in geophysics advance the understanding of Earth processes by estimating and interpreting Earth structure (e.g. velocity and/or density structure). Different types of geophysical data can be collected and analysed separately, sometimes resulting in inconsistent models of...

Full description

Autores:
Sosa Aguirre, Uram Anibal
Tipo de recurso:
Article of investigation
Fecha de publicación:
2013
Institución:
Universidad ICESI
Repositorio:
Repositorio ICESI
Idioma:
eng
OAI Identifier:
oai:repository.icesi.edu.co:10906/78332
Acceso en línea:
http://www.scopus.com/inward/record.url?eid=2-s2.0-84887566536&partnerID=tZOtx3y1
http://hdl.handle.net/10906/78332
Palabra clave:
Soluciones numéricas
Teoría inversa
Sismología computacional
Numerical solutions
Inverse theory
Computational seismology
Automatización y sistemas de control
Automation
Control system
Rights
openAccess
License
https://creativecommons.org/licenses/by-nc-nd/4.0/
Description
Summary:Many experimental techniques in geophysics advance the understanding of Earth processes by estimating and interpreting Earth structure (e.g. velocity and/or density structure). Different types of geophysical data can be collected and analysed separately, sometimes resulting in inconsistent models of the Earth depending on the data used. We present a constrained optimization approach for a joint inversion least-squares (LSQ) algorithm to characterize 1-D Earth's structure. We use two geophysical data sets sensitive to shear velocities: receiver function and surface wave dispersion velocity observations. We study the use of bound constraints on the regularized inverse problem, which are more physical than the regularization parameters required by conventional unconstrained formulations. Specifically, we develop a constrained optimization formulation that is solved with a primal-dual interior-point (PDIP) method, and validate our results with a traditional, unconstrained formulation that is solved with a truncated singular value decomposition (TSVD) for a set of numerical experiments with synthetic crustal velocity models. We conclude that the PDIP results are as accurate as those from the regularized TSVD approach, are less affected by noise, and honour the geophysical constraints. © The Authors 2013. Published by Oxford University Press on behalf of The Royal Astronomical Society.