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...
- 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/
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. |
---|