geoomorfo

CAPITULO 3
PRUEBA DEL METODO EN DATOS SONTETICOS
Para la puesta a prueba del software que aquí se presenta se usó un modelo de prisma usando los métodos de modelado directo que se analizan en el artículo del Journal de Geodesy vol. 74 año 2000 con el título “The gravitational potential and it´s derivatives for the prism” por Nagy, Papp y Benedek.
Se generó una bacteria de archivos .dat queincluyen el campo sintético y su respectivo tensor para cada punto. En forma de matriz n*n, se calculó la anomalía en las tres direcciones (Gz, Gx, Gy). Para la delimitación de estructuras en planta se usa el Gzz, la señal analítica y el radio invariante. Mientras que para estimar la profundidad, usaremos la deconvolución de Euler en 3-D y el programa que provisto para la deconvolución tensorial.Nuestro algoritmo solo ocupa las componentes del tensor y la anomalía gravimétrica (Gz, pero se muestran todos los componentes del vector gravimétrico para ver el comportamiento del campo potencial gravitacional producido por un prisma cúbico en el vacío.
Los parámetros con los que se generó este modelo son los siguientes:
Se  consideran  observaciones  sobre  un  plano  sin topografía  (Z  =  0.0  km).
Unidades:
1. Distancia:  [km]  
2. Densidad:  [kg/m3]  
3. Componente  gravimétrica  vectorial:  [mGal]  
4. Componente  gradiométrica  tensorial:  [mGal/km]  Eötvös  
Coordenada  X  inicial:  0.0  km  
Coordenada  X  final:  1.0  km  
Coordenada  Y  inicial:  0.0  km  
Coordenada  Y  final:  1.0  km  
No. observaciones  en  X:  251  puntos  
No.  observaciones  en  Y:  251  puntos  
Intervalo  de  muestreo  en  X:  0.004  km  
Intervalo  de  muestreo  en  Y:  0.004  km  
Longitud  de  la  malla de  observaciones  en  Dirección X:  1.0  km  
Longitud  de  la  malla de  observaciones  en Dirección Y:  1.0  km  
Prismarectangular causante de perturbación en el campo gravimétrico detectado:
Dimensiones del  prisma:  0.2  km  x  0.2  km  x  2.0  km
El  prisma  se  extiende  desde  la  coordenada X  =  0.4  km,  hasta  X  =  0.6  km  y  desde la  coordenada  Y  =  0.4  km,  hasta  Y  =  0.6  km
Profundidad  a  la  cima  del  prisma:  0.05 km
Profundidad  a  la  base  del  prisma:  2.05  km  
Contraste  de  densidad:  entre  el  prisma  y  el  medio  circundante:  1000.0  kg/m3  












Fig. 3.1  Vector gravimétrico (Modelo  sintético  de  un  cubo)











Fig.  3.2  Tensor  gradiométrica  gravimétrico  completo  (FTG)
En las imágenesanteriores se puede observar las anomalías  en x,  y,  z  además de su respectivo tensor gravimétrico completo. Se ve como las componentes verticales, tanto del vector como del tensor son los que mejor nos definen la fuente causante. Haciendo un análisis interpretativo se pueden inferir algunas fuentes usando solo la anomalía de Bouguer y su derivada  vertical, pero no es suficiente. Acontinuación se presentan un par métodos de ayuda para la identificación de fuentes anómalas en planta.

3.1.- Delimitación de estructuras en planta
Para éste punto como se mencionó con anterioridad, se usara la señal analítica gravimétrica propuesta por Roest Et al. 1992. También se presentara la gráfica obtenida usando los radios  invariantes usados en el cálculo de la solución de profundidad denuestro programa. Éste nos  da los  archivos para graficar éstas imágenes.
 













Fig. 3.3 Señal Analítica y Radio Invariante  


En éstas últimas imágenes se observa la utilidad de éstas herramientas para la identificación de fuentes anómalas en planta. De hecho por el parecido de la señal analítica con los respectivos gradientes se puede dar uno cuenta, que...