Análisis Deformación Frágil PDF
Análisis Deformación Frágil PDF
Análisis Deformación Frágil PDF
El desarrollo de este documento se centra en la identificación de todos los rasgos sobre el modelo digital de
sombras que permitan establecer que existe alguna condición que se manifiesta a través de la deformación
frágil del terreno. Para ello, se sigue un procedimiento metodológico que se centra en el procesamiento de
los segmentos rectos que pueden considerarse deformaciones frágiles. Estos segmentos se marcan sobre el
modelo digital de sombras y se busca establecer la distribución preferente de dirección y concentración.
Con el propósito de evaluar los resultados obtenidos aplicando esta metodología, se realiza el mismo
procesamiento de información sobre registros oficiales como el mapa geológico de Colombia y la plancha
geológica 427, ambos productos han sido emitidos por el Servicio Geológico Colombiano.
Metodología.
El procedimiento para obtener las zonas preferentes de ubicación de las fracturas frágiles y su distribución
de direcciones consta de una primera etapa común que corresponde al tratamiento del modelo digital de
elevaciones para obtener el mapa de sombras. Esta transformación se logra empleando un software de
sistemas de información geográfica que permita cambiar el estilo de representación del ráster que contiene
las celdas con la información de elevación en píxeles de 12,5m por 12,5m (ASF DAAC, 2015) Se decide
acompañar la representación del mapa de sombras con una rampa de colores continua que permite hacerse
una idea de los valores de elevaciones presentes en el área que se va a estudiar y su distribución espacial.
La rampa de colores seleccionada va desde colores verdes hasta un color café oscuro. El producto obtenido
de realizar esta interpretación corresponde a la Figura 2 donde se especifican algunos municipios aledaños
al CVG y se etiqueta la ubicación del Volcán Galeras. Luego de contar con el modelo de sombras, para el
análisis propio se procede a la identificación y señalización de todas las agrupaciones geométricas o
conjunto de características rectilíneas que probablemente correspondan a efectos del tectonismo regional
que produce fracturas, estas son señaladas con líneas discontinuas, el producto de realizar este estudio
corresponde a la Figura 3 donde se pueden apreciar todos los trazos correspondientes a las fracturas
identificadas. Para facilitar la visualización de los segmentos, se decide representar el mapa de sombras en
escalas de grises.
Para el caso donde se trabaja con información de referencia, se extrae información contenida en la plancha
429 del Servicio Geológico Colombiano que se plasma en la Figura 4 ubicando esta información sobre el
modelo de sombras en escala de grises.
El siguiente paso corresponde a la definición de una cuadrícula sobre toda la extensión del área que se
estudia. La creación de esta cuadrícula permite discretizar el espacio de trabajo y asignar un valor
representativo al área que encierra cada celda, el tamaño de celda escogido es de 2,5km por 2,5km. Para
este caso, se define el centroide de cada una de las celdas de la cuadrícula como se muestra en la Figura 1
con un punto negro donde se asigna la razón entre la suma de las longitudes de alineamientos de fracturas
contenido en su celda y el área de la celda correspondiente, lo que al final arroja valores de longitudes de
fracturas por área. Los centroides se definen con el propósito de evitar errores topológicos al momento de
interpolar sobre la superficie para definir un resultado continuo sobre la extensión.
Con la cuadrícula y los centroides definidos, se procede a realizar la intersección entre los segmentos o
alineamientos de fracturas identificados y la cuadrícula, se realiza la suma de todos los segmentos que
quedan contenidos en cada celda y se divide entre el área de la celda para asignar este valor a cada uno de
los centroides. Una vez se cuenta con el valor que se requiere en cada uno de los puntos, se procede a usar
el complemento de mapa de calor que trae el software Qgis que trabaja usando estimaciones de densidad de
Kernel, es decir que se debe suministrar una forma o matriz preferente como Kernel o filtro para que el
software sea capaz de crear una matriz en formato ráster de densidad a partir de la capa de puntos. El cálculo
de la densidad se realiza con la concentración de puntos agrupados, como se cuenta con una distribución
uniforme se debe asignar el peso de la densidad por el valor de la relación descrita anteriormente y debido
a la forma rectangular de la cuadrícula, lo más apropiado es definir un kernel con esta misma geometría. El
resultado de este procesamiento de la información corresponde a la Figura 5 y la Figura 6 que emplea la
rampa de colores característica en mapas de calor que va de colores fríos para los valores de densidad más
bajos y colores calientes para los valores de densidad más altos.
Finalmente, para el análisis de la tendencia de direcciones de los segmentos identificados y de las fallas
reportadas en la información regional de referencia, se realiza un script en Python que toma cada uno de los
segmentos de recta e identifica las coordenadas de sus extremos y con estos valores calcula su azimut.
Luego, recopila todos los valores de azimut calculados y gráfica en coordenadas polares la distribución de
frecuencias relativas dichos valores. Los resultados de este procesamiento se encuentran en la Figura 7 para
las fracturas identificadas en el modelo digital de terreno, la Figura 8 para las fallas recopiladas de la
información de referencia y en la Figura 9 donde se realiza una división por cuadrantes de los resultados de
las fractura identificadas sobre el MDT.
Figura 2. Mapa de elevaciones del área de estudio. Creación propia a partir de: ASF DAAC 2015. SISTEMA DE REFERENCIA:
MAGNA-SIRGAS 3115, Proyección conforme de Gauss - Origen Bogotá, Dátum MAGNA-SIRGAS, Coordenadas Planas
1'000.000 metros Norte - 1'000.000 metros Este, Coordenadas Geográficas 4º35'46''.3215 Lat. N.-77º04'39''.0285Long. W.
De acuerdo con lo que se puede observar en la figura de arriba, se puede decir que existe un levantamiento
con dirección noreste – suroeste que incluye al CVG, este levantamiento corresponde a los Andes
colombianos y se observa cómo este tiende a ser más pronunciado en la parte sureste, esto se explica por la
continuidad que tienen las cordilleras colombianas al norte de la imagen que deriva en una región con
elevaciones más uniformes. La prevalencia de las tonalidades del color café en casi el 50% de la imagen
evidencia que el área de estudio corresponde a una zona elevada a más de 2000 msnm. Se pueden identificar
unas intrusiones muy pronunciadas en la parte norte y al sureste que muestran grandes pendientes con zonas
muy encañonadas que se profundizan bastante con respecto a la zona elevada, uno de estos rasgos se
encuentra muy cerca del volcán Galeras. Finalmente, se puede observar que el CVG hace parte de las zonas
más elevadas del área de estudio.
Figura 3. Interpretación de fracturas frágiles. Creación propia a partir de: ASF DAAC 2015. SISTEMA DE REFERENCIA :
MAGNA-SIRGAS 3115, Proyección conforme de Gauss - Origen Bogotá, Dátum MAGNA-SIRGAS, Coordenadas Planas
1'000.000 metros Norte - 1'000.000 metros, Coordenadas Geográficas 4º35'46''.3215 Lat. N.-77º04'39''.0285Long. W.
Figura 4. Fallas geológicas identificadas en la información de referencia Plancha 429 S.G.C. Cepeda & Murcia, 2008. Creación
propia a partir de: (Cepeda & Murcia, 2008), ASF DAAC 2015. SISTEMA DE REFERENCIA : MAGNA-SIRGAS 3115,
Proyección conforme de Gauss - Origen Bogotá, Dátum MAGNA-SIRGAS, Coordenadas Planas 1'000.000 metros Norte -
1'000.000 metros, Coordenadas Geográficas 4º35'46''.3215 Lat. N.-77º04'39''.0285Long. W.
En las fallas extraídas de la información de referencia se puede observar que corresponden a alineamientos
mucho más extensos que los identificados en el estudio del modelo digital de elevaciones. Se puede ver que
existe una fuerte correspondencia entre las fallas identificadas y lo que se puede observar en el mapa de
sombras. Sin embargo, se ve que existen zonas grandes donde hay rasgos muy fuertes que evidencian
fracturamiento frágil del terreno y que no se encuentran fallas identificadas. Puede ser que el carácter
regional de la información de referencia se centre en la identificación de las extensas zonas de falla con
dirección norte – este que se puede identificar a simple vista sobre la imagen. Además, se puede ver que
hay una falla identificada que atraviesa el CVG y otra que al parecer finaliza en este lugar.
Figura 5. Densidad de fracturamiento inferida por el DEM. Creación propia a partir de: ASF DAAC 2015. SISTEMA DE
REFERENCIA : MAGNA-SIRGAS 3115, Proyección conforme de Gauss - Origen Bogotá, Dátum MAGNA-SIRGAS,
Coordenadas Planas 1'000.000 metros Norte - 1'000.000 metros, Coordenadas Geográficas 4º35'46''.3215 Lat. N.-
77º04'39''.0285Long. W.
El análisis de densidad de fracturamiento a partir de las fracturas frágiles identificadas del modelo digital
de elevaciones muestra que, como ya se había mencionado anteriormente, la mayor concentración de
longitud de fracturas por área analizada se encuentra en la parte sur de la imagen donde se había identificado
el levantamiento de los Andes colombianos. Además, se observa cómo la zona más elevada presenta la
menor cantidad de fracturas agrupadas, cerca al volcán Galeras se observa que la tendencia de densidad
fracturamiento presenta colores cálidos que presentan alrededor de 1km de fracturas por cada km 2 de área.
Un sector de interés corresponde a inmediaciones del municipio de Linares donde se presenta una de las
máximas concentraciones de fracturas, una situación similar sucede en el municipio de Nariño al norte de
la ciudad de Pasto, esta última presenta una baja densidad de fracturas encontrándose en el rango inferior.
La distribución en la paleta de colores se realizó por cuantiles de datos para obtener una mejor representación
de los datos.
Figura 6. Densidad de fracturamiento inferida de la información de referencia. Creación propia a partir de: (Cepeda & Murcia,
2008), ASF DAAC 2015. SISTEMA DE REFERENCIA : MAGNA-SIRGAS 3115, Proyección conforme de Gauss - Origen Bogotá,
Dátum MAGNA-SIRGAS, Coordenadas Planas 1'000.000 metros Norte - 1'000.000 metros, Coordenadas Geográficas
4º35'46''.3215 Lat. N.-77º04'39''.0285Long. W.
Figura 7. Rosa de vientos de la distribución de frecuencias Figura 8. Rosa de vientos de la distribución de frecuencias
relativas de la dirección de las fracturas identificadas. relativas de la dirección de las fracturas encontradas en la
Fuente: Creación propia. información de referencia. Fuente: Creación propia.
Por último, debido a que con la observación de la Figura 3 se pudo identificar algunas zonas donde existe
un cambio en la orientaciones de las fracturas en zonas delimitadas por algunas fracturas mayores, se decide
realizar una sectorización en cuatro cuadrantes de esta figura y definir para cada uno de estos una rosa de
vientos que muestre la tendencia de dirección desagregada por cuadrantes, el resultado de este
procedimiento se muestra en la Figura 9, allí es posible observar que en los cuadrantes NW y SE la tendencia
de las direcciones es muy similar a la obtenida para toda la región. Sin embargo, en los cuadrantes NE y
SW la tendencia de las direcciones presenta un aumento significativo de las fallas con dirección más cercana
a N-S. De esta manera se observa la heterogeneidad en la distribución de direcciones de las deformaciones
frágiles que puede estar asociada con las fallas mayores identificadas anteriormente, que están moldeando
la configuración del fracturamiento. Con esta desagregación con cuadrantes se puede observar la
contribución de cada área de estas a la conformación de la rosa de vientos general.
Figura 9. Sectorización por cuadrantes de la distribución de frecuencias relativas del azimut de los alineamientos de fracturas
inferidas representados en rosas de vientos. Fuente: Creación propia.
Conclusiones
Gómez, J., Montes, N.E., Nivia, Á., Tejada, M.L., Jiménez, D.M., Sepúlveda, M.J., Osorio, J.A.,
Gaona, T., Diederix, H., Mora, M. & Uribe, H., compilers, 2007. Geological Map of Colombia. Scale
1:1,000,000. INGEOMINAS. Bogotá.
Cepeda, H., Murcia, A. 2008. Geología de la Plancha 429 Pasto. Escala: 1:100.000. INGEOMINAS.
Bogotá