Tesis
Inversión conjunta 3D de múltiples datos geofísicos usando restricciones de gradientes cruzados
Joint 3D inversion of multiple geophysical data using crossgradient
Autor
Emilia Fregoso Becerra
Institución
Resumen
En este trabajo generalizamos la metodología de inversión conjunta mediante gradientes
cruzados a medios tridimensionales e introducimos un procedimiento de solución basado en una
formulación estadística y restricciones de igualdad que resaltan la similitud geométrica entre dos
parámetros físicos del subsuelo. Desarrollamos la formulación y algoritmo matemáticos usando
la metodología de inversión generalizada propuesta por Tarantola y Valette (1982), obteniendo
una solución analítica que incorpora términos de solución del problema de mínimos cuadrados
para inversión separada de ambos parámetros del modelo, así como términos que involucran las
restricciones de igualdad dadas por los gradientes cruzados. Utilizamos descomposición en
valores singulares (DVS) para calcular la pseudoinversa de la matriz no invertible que aparece en
la solución analítica. Aplicamos la solución propuesta a inversión conjunta 3D de combinaciones
de datos gravimétricos, magnéticos y de sísmica de refracción y evaluamos las ventajas de esta
formulación en experimentos sintéticos y de campo. In this work we generalized the cross-gradients joint inversion methodology to tridimensional
media and introduced a solution procedure based on a statistical formulation and equality
constraints that enhances the structural similarity between the distributions of two physical
parameters in the subsurface. We developed the mathematical formulation and computational
algorithm using the generalized inversion methodology proposed by Tarantola and Valette
(1982), obtaining an analytical solution that incorporates least squares terms for separate
inversion of both model parameters, as well as terms that account for cross-gradient equality
constraints. We used singular value decomposition (SVD) to calculate a pseudo-inverse of the
non-invertible matrix that is found in the analytical formulation. We then applied the 3D joint
inversion formulation to paired combinations of gravity, magnetic and seismic refraction data
and tested its advantages on synthetic and field examples.