Gecttsic a Colombiana pp. 9-16 diciembre de 2002 Bogota, D.C. ISSN - 0121 - 2974 Amenaza sismica por metoda determinista en el triangulo de fallas Valera-Bocon6-0ca en el occidente de Venezuela RAYMI A. CASTILLA, FRANK A. AUDERMARD M. Y JOSE ANTONIO RODRiGUEZ Departamento de Ciencias de la Tierra, Funvisis, EI Llanito, Caracas, Venezuela, rcastilla@funvisis org.ve RESUMEN EI desarrollo de multiples proyectos de geologia de fallas activas en Venezuela ha permitido que el pais cuente hoy dia con una base de datos completa y de gran cali dad, usada como complemento a la informacion instrumental en los calculos de amenaza sismica del pais con metodos probabilistas, En el presente trabajo se utiliza esta informacion geologica disponible para realizar calculos de amenaza con metodos deterministas para luego comparar los resultados y proponer la irnplantacion de metodolo- gias mixtas que permitan aprovechar mejor los datos, tanto en la ventana instrumental como en la geologica. Sobre la base de disponibilidad de informacion publicada y de complejidad en la geometria de las mayores estructuras, se ha escogido una zona del occidente del pais limitada por las fallas de Valera, Bocono y Oca. Se utilizan sismos maximos creibles de cada falla, distancias rninimas a la superficie de ruptura y las leyes de atenuacion necesarios para crear un marco multiescenario combinando todas las estructuras que puedan, potencial mente, generar amen a- za. La zona de estudio, con un area aproximada de 45.000 km2 que incluye tres capitales de estado y una actividad industrial importante, es cubierta con una red espaciada cada 5 kilometres para el calculo de los val ores de aceleracion. Los resultados se muestran en mapas con areas de isoaceleraciones donde se puede observar un control marcado de la geo- metria de las fallas activas sobre los val ores de aceleracion esperados, por encima de la sismicidad asociada a cada una de elias. A difercncia de los metodos de uso mas extendido, el calculo por metodos deterministas no permite la visualizacion de proba- bilidades de no excedencia, pero permite analizar escenarios de amenaza asociada a fallas que presenten poca sismicidad en la ventana instrumental. PALABRAS CLAVE: AMENAZA SiSMICA, PROBABILISTA, DETERMINISTA, VALOR B, BOCONO, VALERA, OCA, VENEZUELA ABSTRACT The development of several studies regarding active faults geology has allowed us to build a complete database wich is used together with the instrumental catalog, in the probabilistic seismic hazard assesment (PSHA). The geologic data is used here to develop a deterministic seismic hazard assesment (DSHA) and then the results are compared with former PSHA results to pro- pose hybrid approaches that lead us to make good use of all the available information. Maximum credible earthquakes and minimum actual distances to the rupture surfaces are used with the attenuation rela- tionships to create a scenario describing the hazard in terms of the maximum acceleration. These strong ground motion values are the result of the maximum earthquake that could be triggered by the faults in the area. The whole zone, including three main cities and important economic activities, is covered in this study by a grid with points every 5 km. The. results are shown in acceleration maps, where the maximum values follow the active faults trends. When the return period of the maximum credible earthquake is included in the analysis this correlation changes and a more likely scenario is created. The pros and cons of both approaches were analized in terms of the available data for Venezuela. This has shown that more efforts must be done in the treatment and processing of the data and in the design and application of better metho- dologies as well. KEY WORDS: SEISMIC HAZARD, PROBABILISTIC, DETERMINISTIC, B-VALUE, BDCDNO, VALERA, OCA, VENEZUELA METODOLOGiA Se ha realizado un analisis del catalogo sismico nacional en termi- nos del valor b y de la relaci6n Gutenberg-Richter para evaluar las limitaciones que pudieran estar asociadas a la utilizaci6n de dicha base de datos en los calculos de amenaza sismica que se realizan en el pais. Dicho analisis se ha concentrado en las variaciones de estos GEOFislCA CDLOMBIANA, 6, DICIEMBRE DE 2002 9 CASTILLA ET AL. pararnetros a traves del tiempo y en la consideracion que dichas va- riaciones son el producto de cambios en la instrumentacion, y no de cambios naturales en el comportamiento de la sismicidad del pais. Se realizaron tambien dos ensayos de calculo de amenaza sismi- ca por metodos deterministas, los cuales han permitido, de manera similar al caso anterior, exponer las ventajas y desventajas, asi como las limitaciones que el metodo impone. Para este fin se ha to- rnado la informacion de geologia de fallas activas hasta ahora gene- rada en el pais, en terminos de sismos maximos creibles con sus respectivos periodos de retorno, la ubicacion de todas las fallas y su segmentacion (Audemard et al., 2000). Se utilize la metodologia clasica (figura I) para el enfoque de- terminista (Krinitzsky, 199 ; Reiter, 1990 y otros) con la diferencia que, en el presente ensayo, se toman en cuenta multiples puntos para generar una malla y mapas de isoaceleraciones. En ambos casos se ha escogido una reticula de 5 km de espacia- miento, y para todos los puntos se ha calculado la aceleracion que produce cada falla al desencadenar su sismo maximo probable. Usando como parametro de distancia, la minima entre el punto y la falla, dichos val ores de aceleracion son comparados entre si para escoger el maximo y asignarlo al punto. ~nle3 Si~ ~'1 PASO I PASO) • U , '0 y, ~ >- 0 " 0 Y,'s ~ y = .~ ~ 0 YN::E R, PASO 4 PASO 3 Figura 1. Etapas basicas del anal isis de amenaza sismica por medio del enfoque pro- babilista. Modificado de Smit (2002). Para considerar el periodo de retorno de cada uno de los sismos maximos probables, se multiplico la aceleracion por la frecuencia anual de dicho evento. De esta manera se obtiene un valor conside- rado adimensional, utilizado solo para escoger el valor de acelera- cion maxima que muy probablemente afectaria a cada punta. 10 METODO PROBABILISTA Y EL CATALOGO SiSMICO NACIONAL En Venezuela, la sismicidad puede ser considerada moderada cuan- do se Ie compara dentro de un contexto mundial (Romero et al., 1999). EI 72,4% del catalogo instrumental, desde 1910 hasta hoy, corresponde a sismos con magnitudes 36.0) en estos 92 afios, solo constituye un 0,46% del total del catalogo. No obstante, al amp liar la ventana a observaciones historicas (desde 1530 hasta hoy), se obtiene un re- gistro de 16 sismos destructores con perdidas materiales y de vidas humanas. Sismicidad - Venezuela (1910-2002) 1200~ 1000f-- 400-- 4 S 6 Magnitudes Figura 2. Distribuci6nde sismicidad instrumental en Venezuela, segun lamagnitud. La figura 3 muestra la distribucion de los epicentros en el catalo- go de sismicidad, discriminados segun su magnitud y la profundi- dad hipocentral, respectivamente. En las mencionadas figuras se puede observar el marcado control que ejerce la presencia de las fa- lias mas importantes en la concentracion de la mayoria de los epi- centros para eventos registrados en el catalogo. EI catalogo instrumental puede ser dividido en tres periodos que reflejan, en diversos sentidos, las distintas etapas evolutivas de la red sismologica nacional. El primer periodo, comprendido entre 190 I Y 1982, durante el cual son puestos en funcionamiento distin- tos sismografos para monitoreo. Para el momenta de la creacion de Funvisis, luego del sismo de Caracas de 1967, existian 8 estaciones sisrnologicas con operadores en sitio, pertenecientes a distintos or- ganismos (Rendon, 1997). El segundo periodo comienza en 1981 cuando se inicia el Pro- yecto Resvac, que permitiria la incorporacion de las estaciones con- vencionales existentes a una red con 10 estaciones telemetricas de sensores de periodo corto distribuidas en la region central. EI arre- glo resultante de este proyecto, junto con otras redes de caracter re- GEOFislCA COlOMBIANA, 6. OICIEMBRE OE 2002 AMENAZA SISMICA POR METODO DETERMINISTA EN EL TRIANGUlO DE FAllAS VAlERA.BOCONO·OCA EN El OCCIDENTE DE VENEZUELA a) b) Figura 3. a) Distribuci6n espacial del catalogo de sismicidad en Venezuela con epicentros discriminados segun la magnitud. Azul = mb<5.0; Verde = 57. b) Distribucion espacial del catalogo de sisrnicidad en Venezuela con epicentros discriminados segun su profundidad. Amarillo = somera; Naranja = Intermedia; Negro = Pro- gional implantadas por otros organism os, ha permitido el procesamiento de un numero importante de eventos, 10 cual repre- senta un mejoramiento sustancial en el catalogo sismico nacional (Rend6n, op. cit.). EI comienzo del tercer periodo puede ubicarse en el ana 1999, cuando se pone en marcha el proyecto de modemizaci6n de la red sismo16gica naciona!. Dicho proyecto, todavia en fase de ejecucion, contempla la instalaci6n de sensores de tres componentes del tipo banda ancha. Todos estos cambios en la configuraci6n y densidad de la red na- cional han tenido repercusiones importantes en el catalogo sismico, La figura 4 refleja dichos cambios enla frecuencia de registro por ano. Cambios importantes se pueden observar despues de 1980, cuando la Resvac entra en funcionamiento, y poco antes del afio 2000, cuando comienza la modemizaci6n de la red naciona!. Sismicidad - Venezuela (1910 - 2002) 1200 ,--,----,--,-~,__,_--,--,----,_-,-__,__, 1000 800 200 1910 1920 1930 1940 1950 Figura 4.Registro de actividad sismica anual. GEOFlslCA COLOMBIANA, 6, OICIEMBRE OE 2002 Otros picos importantes pueden ser observados en asociaci6n con eventos sismicos. Se puede observar un aumento considerable de la frecuencia en el afio 1967 asociado al sismo de Caracas y, en 1989, asociado a la tormenta sismica que tuvo lugar en Falc6n oriental. Los cambios antes descritos en la red sismol6gica nacional, al afectar la frecuencia y el nivel de detecci6n de los multiples even- tos que suceden cada ano, produjeron tam bien cambios en la infor- maci6n que se extrae del catalogo, para calcular arnenazas. Se puede observar, de distintas maneras, la influencia que ejercen es- tos cambios del catalogo disponible al realizar el procesamiento de los datos. La figura 5 muestra la relacionnes Gutenberg-Richter (b-lines) para Venezuela en las etapas previa a la instalaci6n de Resvac y durante el funcionamiento de la misma, respectivamente. 4 Magnitud Re laci6n de recurrencia 10' ,--~-~-----_---~-_-- '~ LogN-a,344-1,15M o 1983·1998 1910·1979 I ~ 10' Z 0/) 0 oJ 10' 10.2 10'" 0 LogN-4,451.0.84H.1' Figura 5. Relaci6n Gutenberg-Richter para la primera y segunda etapa descritas an- teriormente. 11 CASTILLA ET AL. Ambas relaciones se ajustan bien a los datos del catalogo entre las magnitudes 4 y 6. Por encima de este intervalo, el ajuste es mas bien deficiente y las frecuencias son subestimadas. Adernas de este Iimitado intervalo donde el ajuste puede considerarse confiable, las diferencias que presentan ambas relaciones entre si es otra caracte- ristica poco satisfactoria. Por ejemplo, para un sismo de magnitud 7.5, cada una de las relaciones propone un periodo de retorno de 72 y 191 afios, respecti vamente. La figura 6 permite observar distintos valores b calculados para el catalogo de sismicidad del pais. Cada valor es calculado para to do el periodo comprendido entre 1910 y el afio indicado en el eje de las abscisas. EI comportamiento de la curva permite observar, de manera similar a las anteriores, una diferenciacion entre los dos pe- riodos de registro. En este ejemplo, el cambio en el comportamiento de la grafica se encuentra desplazado con respecto al afio de la in- corporacion de la nueva red debido al caracter acumulativo de la curva, 10 que no permite el retlejo instantaneo de las variaciones. A partir del periodo 1980- I985, la curva tiende a la convergencia alre- dedor del valor 0,9. Aunque dicha convergencia parece indicar una estabilizacion del valor b, 10 que refleja el mejoramiento de la red sismologica nacional, se debe seguir prestando atencion a los pri- meros afios del catalogo, pues estes influyen de manera significati- va en los estudios realizados tomando en cuenta que dicho catalogo es muy completo. Cuando el termino b se calcula de modo similar para period os anuales (figura 7), se observa no solo una variacion importante en el valor del mismo en terminos absolutos, con una variacion de 1,4 unidades entre el minimo y maximo de la curva, sino tambien el re- flejo de los cambios en la configuracion de la red como un cambio en la frecuencia de los picos en la curva. EI aumento en el valor b medio, luego de la instalacion de la red Resvac (de 0,94 a 1,05), re- fleja el cambio de umbral de deteccion, 10 cual permite registrar en el catalogo mayor cantidad de sismos de menor magnitud. Valor b cal cui ado para el catalogo desde 19I0 0.951------::------------------- 0.9 -c 0.85r--(f--+----------:r"-------:-------- Ci 'Oil 0.8r----\--------:;;/"--------------- > 0.75 r----\---rJf'''---------------- 0.7r----t--j------------------j 0.65r----~.-----------------j 0.6';";----;::-:-:-_--:--:-:-:----:-=-:-_=-:_-:-':--_---'-_----L_----' __ '-----.I 19551960196519701975198019851990199520002005 Alios Figura 6. Variacion del valor b en period os cuya duracion se incrementa a razon de cuatro afios. 12 Variaci6n anual del valor b 1970 1975 1980 1985 1990 1995 2000 Afio Figura 7. Variacion del valor b calculado para periodos anuales. Todo 10 antes descrito permite observar un hecho logico espera- do: existe una correlacion directa entre los cam bios en la configura- cion de la red nacional y los cambios en las caracteristicas del catalogo de sismicidad del pais. De esta manera, cualquier estudio que implique el procesamiento y uso del catalogo sismologico na- cional (entre ell os el calculo de amenaza sismica), debe pasar nece- sariamente por una revision minuciosa del mismo, aplicando criterios logicos y coherentes que permitan obtener resultados cada vez mas apegados a la realidad sismica del pais. De igual manera, la metodologia de calculo de amenaza por meto- dos probabilistas, cuando es aplicada en su forma basica, presenta se- rias limitaciones que no deben ser menospreciadas. En Krinitzsky (1993; 1995), entre otros, se han desarrollado extensas discusiones acerca de estas limitaciones. La relacion Gutenberg-Richter, en su concepcion original, describe de manera eficiente la sismicidad de un area extensa. Cuando se realiza un analisis de amenaza, se supone que esta relacion regional describe de manera satisfactoria el com- portamiento de fallas individuales. Esta asuncion elimina la-nocion segun la cual, en un momenta dado, cada falla, 0 incluso cada seg- mento de una misma falla, puede encontrarse en momentos distintos dentro de su propio cicio sismico. La accion de uniformar los val ores b de toda una region impide la posibilidad de observar singularidades inherentes a cada falla, que en definitiva pudieran ser controles importantes para el com- portamiento sismogenico de la misma. Diversos ensayos han sido realizados en la busqueda de modelos y metodologias que reflejen de manera mas efectiva estas propiedades, entre los que podemos nombrar los realizados para Venezuela por Sobiesiak et al. (200 I) y Grosser et al. (2001). GEOFislCA COLOMBIANA, 6, DICIEMBRE DE 2002 AMENAZA SiSMICA POR METODO DETERMINISTA EN El TRIANGULO DE FALLAS VALERA.BOCON6·oCA EN EL OCCIDENTE DE VENEZUElA TABlA I. TABLA RESUMEN CON LASCARACTERisTICAS MAs IMPORTANTES DE LAS 13FALLAS INCLUIDAS EN IA ZONA DE ESTUDIO. MODIFICADA DEAUDEMARD ETAL. (2000). Nombre de la Sentido de Movimiento Edad del ultimo Tasa de movimiento Sismo maximo Periodo de Estructura (mavor/menor) movimiento (mm/afio) creible retorno (afios) Oca Dextral <15 ka 045-20 74 4295±1020 Ancon Dextral <15 ka 16 74 1900 Camare-Paraiso Dextral/lnversa Probablemente <15 ka ?20 7 I 670 --- DextrallN ormal Probablemente <15 ka ?30 67 390 Socremo Inversa/Dextral <16Ma <20 67 390 Urumaco (oeste) Dextral <16 Ma 005 66 14000 Urumaco (este) Dextral <15 Ka 005 64 12000- 15000 Rio Seco Dextral <15 Ka 035 70 1500 Valera Sinestral <15 Ka 10 69 2000 Bocono (Anzoategui- Dextral <15 Ka 5 7,5 400 Barquisimeto) Bocono (Cabudare- Dextralllnversa <15 Ka 1-3 7,5 700-2000 Moron) Flexura Frontal Reverse <1,6 Ma 0,5 7,0 500-750 Surandina La Victoria-Guacamaya Dextral <1,6Ma 0,6 7,0 2000 EI metoda determinista. Ensayo en el occidente del pais EI triangulo de fallas de Valera-Bocono-Oca, que incluye territorio de varios estados, (Merida, Zulia, Trujillo, Lara, Falcon y Portugue- sa), presenta una complejidad estructural y cinematica que merece especial atencion. Se encuentra limitada por tres fall as de movi- miento transcurrente, de velocidades distintas, 10 que sugiere la existencia de complejidades cinematicas intemas y extemas. Dichos Iimites los constituyen las fallas de Oca-Ancon, Bocono en sus seg- mentos de Mucuchies-Anzoategui, Anzoategui-Barquisimeto y Ca- budare-Moron, y Valera en sus dos segmentos (Valera y rio Momboy) (Audemard et al., 2000). Con el fin de lograr una mejor y mas completa cobertura de la mall a dentro de la zona del triangulo, el area de estudio se amplio, incluyendo un total de 13. fallas activas, potenciales generadoras de sismos en la region. La figura 8 muestra la distribucion de la sismici- dad en la zona de estudio, evidenciando la ocurrencia de eventos dispersos, no necesariamente asociados de manera directa a las es- tructuras mayores. La relacion de recurrencia para la region presenta caracteristicas similares a la calculada para la completitud del cata- logo, el intervalo con comportamiento lineal es corto y por encima de magnitud 6.0 el ajuste de los datos con la lineal es deficiente. Se han generado mapas de amenaza sismica usando la informa- cion disponible de geologia de fallas activas. Todos los datos usa- dos, son el resultado de estudios de paleosismologia, y en general de neotectonica en el area escogida. Datos como el sismo maximo crei- ble, su periodo de retorno y la geometria detallada de las trazas de GEOFlslCA COLOMBIANA, 6, OICIEMBRE DE 2002 cada una de las fallas de la region fueron usados para calcular,junto con leyes de atenuacion, las aceleraciones maximas en roca espera- das en cada punto. Claramente se puede observar en los mapas, un control impor- tante que ejerce la geometria de las estructuras consideradas en las tendencias de las zonas con mayores aceleraciones. En eI caso de la figura 9 se tomaron en cuenta solo las magnitudes de los sismos maximos creibles asignados a cada segmento de falla (Audemard, . 2000). La aceleracion asignada a cada punto es la maxima posible, producida por un sismo de magnitud maxima creible para un seg- mento especifico de falla. De esta manera, tal valor solo depende de la magnitud y la distancia minima entre el punto y la falla. El valor maximo de aceleracion en este mapa alcanza los 0,46g, mientras que el minimo es de O,06g. De esta manera es facil entender que los valores maximos de aceleracion se encuen!ren alrededor de las fallas de Oca, Ancon y Bocono asociados a magnitudes 7.4 y 7.5 respectivamente. Mien- tras que los val ores minimos son asignados a los puntos mas distan- tes de las estructuras activas (ej: centro del triangulo), La figura 9 representa entonces, segun los datos geologicos, el peor escenario posible con las caracteristicas que presentan las fa- lias de la zona. Cada valor asignado a cada punta es equivalente a la maxima aceleracion esperada en roca para ese sitio, segun la ley de atenuacion usada. Con el fin de plantear un escenario en algo mas real se incorporo al calculo el periodo de retorno de cada sismo ma- ximo creible, tambien calculado a partir de datos geologicos. 13 CASTILLA ET AL, • ••• •. "', \, $ ,••..... .. \•• ••• •.... .... '. Rclacion de Recurrencia Zona de Estudio 10' ,------.--~-~--~-~--~--,.__, 10' Log N: 4.226 ·O.976M o 10' 0 0 0 0 0 a z ~o -' 10° 10' 10' 10' 10" 0~-~--=-2----:-3----":4----:5'---------:6:-------=----.J Magnitud bl Figura 8, Caracteristicas de la sismicidad en la zona del triangulo, A)Ubicaci6n de epicentros reportados en el catalogo. B) Relaci6n de recurrencia correspondiente, La figura 10 muestra el resultado de la incorporacion de este dato adicional. En este caso la aceleracion es calculada del mismo modo que en el anterior, 10 que varia es el criterio de escogencia del valor. Una vez calculada las distintas aceleraciones que para un punto pudieran producir los distintos sisrnos, se multiplican por la frecuencia anual del sismo correspondiente. Se obtiene asi un valor adimensional cuyo maximo es usado para escoger la aceleracion a asignar. De esta rnanera, se obtiene un mapa que permite observar dentro de la zona, las aceleraciones que de manera mas probable liON Acslsraclonss • 0,4 • 0,3o 0,2 • 0,15 • 0,05 Figura 9. Mapa de aceleraciones calculado con los sismos maxirnos creibles para cada falla. 14 pueden afectar a un punto dado dentro de un escenario de sismos maximos creibles. Los cambios mas importantes se observan en el area alrededor de la falla de Valera, La Victoria y Bocono en su segmento Cabu- dare-Moron con periodos de retorno de 2000 afios para las dos pri- meras y 1000 afios para la ultima. Mientras que las zonas con menos variaciones corresponden a las mas cercanas a fallas con periodos de retorno menores como Bocono en su segmento Anzoate- gui-Barquisimeto (400 afios) y la Flexura Frontal Surandina (625 Aceleraciones • 0,4 • 0,3o 0,2 • 0,15 • 0,05 Figura 10. Mapa de aceleraciones cal cui ado con los sismos rnaximos creibles y el pe- dodo de retorno para cada falla. GEOFlslCA COLOMBIANA, 6, OICIEMBRE OE 2002 AMENAZA SiSMICA POR METODO DETERMINISTA EN EL TRIANGULO DE FALLAS VALERA.BOCONO·OCA EN EL OCClDENTE DE VENEZUELA afios), De igual modo fallas cap aces de producir grandes sismos, que se encuentren lejos de otras estructuras importantes mantienen su influencia en sus zonas aledafias, Del mismo modo que para el enfoque probabilista, los rnetodos deterministas presentan multiples limitaciones para el calculo de amenaza. En este tipo de anal isis, la calidad de la informacion e in- terpretacion geologica surgen como limitantes importantes, de he- cho, en un mismo mapa pudieran coexistir datos de distinta calidad y procedencia, resultando asi en valores menos confiables. La dificultad de contar con sismos de magnitud media, pues su presencia en el registro geologico es usualmente enmascarada por sismos mayores, constituye un hecho que limita esta metodologia al uso de sismos maxim os, resultando en valores de movimiento del te- rreno altos. Dichos valores altos de movimiento del terreno, pueden representar una sobreestimacion de la amenaza para estructuras de corta vida util 0 de poca importancia en el manejo del riesgo. Debido a esto, los mapas aqui generados son de utilidad cuando se calcula la amenaza para obras de ingenieria de gran importancia en don de se deben considerar sismos con largos periodos de retorno. La estrecha asociacion de los val ores de movimiento del terreno con el tamafio del sismo y la fall a que los produjo, permite acotar di- chos valores en un nivel superior. Es decir, mientras el enfoque pro- babilista permite realizar mapas incrementando las magnitudes tanto como el analista 10 desee, la metodologia aqui usada impone un "techo" mediante la magnitud maxima que es capaz de-producir cada accidente tectonico, EI problema de la perdida de nocion temporal dentro del cicio sismico puede ser subsanado cuando se posee un estimado de la fe- cha del ultimo sismo equivalente al maximo creible, de 10 contrario esta limitante sigue estando presente en esta metodologia tal y como 10 esta en el enfoque probabilista. Una ventaja importante que presenta esta metodologia es la posi- bilidad de asociar cad a valor de aceleracion con los parametres fuente y del sismo que 10produce. Esto permitiria en un paso subsi- guiente conocer en mas detalle el contenido de frecuencia responsa- ble de la excitacion del terreno. CONCLUSIONES Los ensayos con enfoques deterministas aqui realizados han demos- trado su utilidad para amp liar la vision de la amenaza y acotar los va- lores segun las limitaciones naturales que imponen las mismas estructuras genreadoras de sismos. EI uso de esta metodoiogia en re- giones con bases de datos geologicas completas y de buena calidad, en el cmapo de las fallas activas, permite generar resultados que son de gran unlidad como complemento en enfoques integrales. No obstante, el igual que la metodologia que ooohace uso de las probabilidades, los enfoques detenninistas han mostrado diversas desventajas. En el ensayo realizado se comprobo que debido al ca- racter netamente geologico de los datos utilizados, la disponibilidad GEOFislCA COLOMBIANA, 6, OICIEMBRE DE 2002 y calidad de los mismos puede cambiar de una falla a otra. La limi- tacion que 10 anterior impone radica en que para la realizacion del anal isis es necesario contar con la uniformidad de todos los para- metros en todas las fallas. Es decir, todas y cad a una de las estructu- ras consideradas como potenciales generadoras de sismos, deben contar con estudios de geologia de fallas activas que hayan deter- minado sismos maxirnos creibles, periodos de retorno y edad del ultimo sismo, todos con niveles de confiabilidad similar. A medida que mas y mejores datos se hacen disponibles, este tipo de mapas pueden ser mejorados de manera significativa, incor- porando criterios adicionales como el estado de esfuerzos alrede- dor de las fall as, que permitan plantear escenarios mas reales, usando igualmente datos con origenes distintos a los estudios geo- logicos. EI analisis del catalogo sismico nacional en terminos del valor b y de la relacion Gutenberg-Richter ha permitido la identificacion de diversas variaciones en el mismo, asociadas a los distintos cam- bios realizados en la instrumentacion de la red nacional y no a un cambio real de la sismicidad del pais. Si dichas variaciones son trasladadas a los calculos de amenaza sismica sin un exam en critico de sus efectos, los resultados obtenidos pudieran representar s610 un reflejo muy aproximado de la realidad. Los mapas asi generados pueden ser de utilidad para la construccion de obras de interes ge- neral 0 para el manejo de emergencias. No obstante, la aplicacion de estos resultados en estructuras vitales 0 que impliquen un bajo nivel de aceptacion al riesgo constituiria un hecho preocupante desde todo pun to de vista. Los analisis aqui realizados permiten visualizar la evolucion de la red sisrnologica nacional y como sus modificaciones han influi- do en el catalogo sismico nacional. En este sentido, los resultados aqui obtenidos son una muestra de los avances que en esa materia se han hecho. Se hace absolutamente necesario adoptar metodologias mixtas que incorporen otros parametres. Esto permitira que los estimados de amenaza sismica arrojen resultados verdaderamente dinarnicos que reflejen de manera mas eficiente las propiedades fisicas que controlan la ocurrencia de sismos y su cambio en el tiempo. REFERENCIAS AUDEMARD,F. A., M. N. MACHETTE,J. W. Cox, R. L.DART & K. M. HALLER, 2000, Map and database of Quaternary faults in Venezuela and its offshore regions, Open-File Report 00-018 USGS. GROSSER, H., A. RIETBROCK, M. BAUMBACH, G. ROMERO, J. ROJAS, 200 I, Seismic Tomography in the aftershock region of the Cariaco earthquake, 1997, Eos. Trans. AGU, 82(47), Fall Meet. Suppl., Abstract S32E-03. KRINITZSKY, E. L., 1993, Earthquake probability in enginee- ring-Part 2: Earthquake recurrence and limitations of Gu- 15 CASTILLA ET AL. temberg-Richter b-values for the engineering of critical structures. The third Richards H. Jahns Distinguished Lectu- re in Engineering Geology, Engineering Geology, v.36, pp. I-52. KR!NITZSKY, E. L., 1995, Deterministic versus probabilistic seis- mic hazard analysis for critical structures, Engineering Geo- logy, vAO, pp. 1-7. REITER, L., 1990, Earthquake hazard analysis. Issues and in- sights, New York, Columbia University Press, 254 pp. RENDON, H., 1997, Red sismica de Venezuela, en AA VV Grases, J., Dise- no Sismorresistente, especificaciones y criterios empleados en Vene- zuela. Caracas, 662 pp. 16 ROMERO, G., J. ROJAS, M. BAUMBACH, H. GROSSER, M. SOBIESIAK y grupo de trabajo RESICA 97, 1999, Distribuci6n en profundidad de la actividad sismica y su relaci6n sismotectonica, VI Congo Venez. Sism. e Ing. Sismica., Merida, Venezuela. SMIT, P, 2002, Seismic Hazard and Risk Assesment, Lecture Notes Certifi- cat d'Etudes des Risques Geologiques 2002. SOBIESIAK, M., C. PALME, J. CHOY, M. T. MORANDI, M. CAMPO, C. GRANADO, 200 I, Seismological parameters in the Northern Andes, Ve- nezuela, Eos. Trans. AGU, 82(47), Fall Meet. Suppi., Abstract SI2B-0599. GEOFislCA COLOMBIANA, 6. OICIEMBRE OE 2002