INVERSIÓN 1D DE DATOS MAGNETOTELÚRICOS DE RESISTIVIDAD APARENTE, UTILIZANDO DOS PARAMETRIZACIONES DIFERENTES: A) RESISTIVIDAD Y CIMA DE CADA CAPA, B) RESISTIVIDAD Y ESPESOR DE CADA CAPA
En la interpretación de datos geofísicos, generalmente se aplica la inversión (modelado inverso) de los datos de campo resolviendo un sistema y=Ax, donde la matriz A relaciona los vectores x e y, que son los parámetros a determinar y los datos de campo, respectivamente. La matriz A (generalmente es rectangular al tener más datos que parámetros a determinar) se construye a partir de la solución al problema directo en el que se “matematiza” el problema físico para determinar la anomalía que se obtiene a partir de ciertos parámetros del modelo. El sistema se resuelve por el método de mínimos cuadrados Se dice que el éxito de la implementación de la solución al problema inverso, depende de una buena elección de los parámetros a determinar, entre otras cosas.
Se presentan dos algoritmos aplicados para la solución del problema inverso (caso 1D) de datos de campo provenientes del método magnetotelúrico; los datos de campo son resistividades aparentes asociadas a ciertos periodos de las señales medidas. En este caso se trata de determinar la estructura de resistividades que puede ser:
a) resistividad y profundidad de la cima de n capas, los parámetros a determinar son 2n; es decir, n resistividades y n profundidades a la cima del modelo de capas a obtener;
b) resistividad y espesor de las capas, los parámetros a determinar son 2n-1; es decir n resistividades y n-1 espesores (la última capa tiene un espesor infinito) del modelo de capas a obtener.
Se presentan varios resultados de los dos algoritmos de inversión, en los cuales se destaca que el algoritmo que funciona mejor es el que está relacionado con la parametrización de resistividades y espesores, que es como se trabaja.