Método interactivo para la extrapolación de ondas sísmicas. Aplicación a la prospección geofísica (Parte 1)

En los análisis sísmicos por reflexión existe una creciente demanda de estimaciones geológicas y geométricas más precisas y sofisticadas. En último término se pretende extraer estimaciones fiables de las secciones de reflexión, estimaciones que afectan a parámetros tales como la velocidad de propagación, densidad y, por supuesto, contenido de hidrocarburos. Ello exige gran precisión en las técnicas de procesado y representación gráfica, pudiendo decirse que los métodos gráficos interactivos son clave en las tendencias futuras de prospección. Un área de investigación actual es el procesado y representación de datos representados por "horizontes". Un horizonte sísmico se define como la superficie que separa dos capas diferentes de roca. Estas interfases están asociadas a reflexiones que se pueden identificar claramente en las secciones sísmicas. Usando una pantalla gráfica adecuada estos horizontes aparecen coloreados en la sección,pudiendo tornarse puntos selectos para definirlas. Si se utilizan secciones concentradas ("stacked") es preciso determinar la posición correcta de las interfases. Es el proceso llamado "emigración". Tradicionalmente se han aplicado fórmulas de representación o transformaciones integrales para realizar este proceso. Sin embargo es suficiente hacer emigrar las curvas seleccionadas (superficies en el caso 3-D) que representan menos del 1% del total de los datos, con lo que el tiempo de CPU se reduce considerablemente. Si se desprecia la información sobre amplitudes, fases y frecuencias la teoría de rayos es el método más rápido y simple para la conversión en profundidad, poniéndose todo el énfasis en el proceso interactivo y en la representación gráfica. Puesto que en estructuras complejas, como las del Mar del Norte o las del Norte de Alemania, las funciones de velocidad nunca se conocen con precisión, es preciso repetir la emigración mediante diferentes velocidades especificadas por el usuario en forma interactiva. En este artículo se presenta la aplicación del procedimiento para horizontes de separación nula (zero-offset sections) y al final se indica la posibilidad de extensión del estudio a casos con separación constante y al análisis de velocidad.

Método interactivo para la extrapolación de ondas sísmicas. Aplicación a la prospección geofísica (Parte 2)

Este artículo es continuación del publicado en el número 2 de la revista "Anales de Ingeniería Mecánica"(diciembre 1982), sobre los métodos de análisis sísmico por reflexión. Como allí se indicó, con estos métodos se pretende determinar la geometría de los horizontes sísmicos en una determinada sección así como las características físicas de las interfases que separa. En el presente artículo se introducen mejoras en lo que se refiere a tiempo compucional y tipos de campos de velocidades. En este nuevo método propuesto para la resolución del proceso de emigración, mediante un método gráfico interactivo, pueden señalarse como ventjas fundamentales, las siguientes: a) Tiempo de resolución acorde con la complejidad del campo de velocidades, produciéndose economía de dicho tiempo, sobre todo en los casos sencillos; b) El número de elementos de la malla y las dimensiones de los mismos son función de la complejidad del campo de velocidades; c) El estudio en cada elemento se realiza de forma general e independiente. General, en cuanto a que los resultados que se obtienen en cada elemento son valores globales, referidos a la malla y al problema en conjunto, e independiente, puesto que los datos necesarios en cada elemento se obtienen a partir de valores nodales; d) Posibilidad de resolución del problema para un campo de velocidades cualesquiera, y por lo tanto, incluyendo los casos de discontinuidades de velocidad, que se evalúan mediante la ley de Snell; e) Posibilidad de extender el método al caso de tres dimensiones, sin más que elegir como elementos de la malla, tetraedros en vez de triángulos. En conjunto, resulta un método flexible, de fácil utilización e interpretación de resultados.

Análisis de parámetros dinámicos y estáticos en la prospección de áridos. Caso práctico: cantera de arcillas en Madrid (España)

En el presente trabajo se hace una valoración de los recursos existentes en una explotación de áridos estudiando parámetros dinámicos mediante técnicas geofísicas y parámetros estáticos mediante ensayos de caracterización de arcillas en laboratorio. Esta doble actuación obedece a la nueva normativa aplicada a las explotaciones mineras en las que se indica la necesidad de la realización de un informe geotécnico incluido dentro de la valoración arqueológica de la explotación. Para ilustrar este “modus operandi” se ha elegido un caso práctico real: una explotación sita en el municipio de Torres de la Alameda, Madrid (España). En cuanto al análisis de los parámetros dinámicos se ha realizado una campaña de prospección geofísica utilizando Corriente Continua mediante la técnica de los Sondeos Eléctricos Verticales (S.E.V.). La selección de los lugares de prospección ha obedecido a criterios estadísticos apoyados en las observaciones de la geología de superficie y las características morfológicas de la explotación. La aplicación de los S.E.V. ha permitido determinar la distribución en profundidad de los diferentes niveles geoeléctricos del área estudiada, además de la distribución de resistividades en cada uno de los niveles. Para el estudio de parámetros estáticos se realizó una investigación de las características geotécnicas del terreno mediante la realización de calicatas geomecánicas y ensayos de penetración dinámica además de ensayos de caracterización de muestras. Dicha información, conjuntamente con las descripciones sobre las muestras obtenidas, han permitido realizar una caracterización fiable y barata de los recursos presentes en la zona de estudio.

Investigación sobre recursos áridos mediante tomografía geoeléctrica en Malpartida de Plasencia (Cáceres)

El objetivo de este proyecto es la caracterización de los materiales presentes en una finca situada de cara a la valoración de la misma. Con este fin se realiza, por una parte, una campaña de prospección geofísica consistente en perfiles de Tomografía Geoeléctrica cuya inversión proporciona las secciones de resistividad del subsuelo. Su análisis permite establecer los distintos horizontes y su posición relativa. Se realiza una revisión detallada del método de tomografía geoeléctrica, desde la adquisición de datos hasta su interpretación pasando por su procesado. Por otra parte se llevan a cabo una serie de calicatas para el reconocimiento de los materiales y la toma de muestras sobre las que realizar ensayos de caracterización. Los resultados permiten asignar una litología a los valores espaciales de resistividad que corresponden distintas granulometrías y grados de compactación. Por último, se realiza una estimación de las reservas de la finca de cara a su valoración económica.

Inspección técnica de la presa Santa Marta, mediante medidas subacuáticas de potencial espontáneo

Las filtraciones de agua, con la consecuente erosión interna en presas de materiales sueltos, es una de las causas principales de fallos y accidentes. Las consecuencias del fallo de estas estructuras, pueden ser, pérdidas tanto económicas como de vidas humanas. Por lo cual en este proyecto se describe la aplicación de un método de prospección geofísica no invasiva, medidas de potencial espontáneo, para detectar posibles filtraciones de agua en el cuerpo de la presa. El flujo de agua a través de un material poroso y permeable crea un campo de potencial eléctrico de una magnitud de decenas o centenas de milivoltios, el cual puede ser medido y así detectar infiltraciones de agua en presas de materiales sueltos. Se ha aplicado esta técnica en la Presa Santa Marta, y mediante una interpretación cualitativa de los datos medidos, tomados en la cara aguas arriba de la presa (medidas subacuáticas), se logró identificar un flujo de agua vertical y otro subhorizontal, que estaban ingresando en el cuerpo de la presa, los cuales estaban causando erosión interna y la formación de una tubificación. ABSTRACT Water leakages and internal erosion in embankment dams is one of the main causes of failures and accidents. The consequences of the failure of these structures may cause losses both, economical and of human lives. Therefore, this project describes the application of a noninvasive geophysical prospecting method, self potential measurements, to detect water leakages in the body of the dam. Water flow through a porous and pervious medium creates an electric potential field with a magnitude of tens or hundreds of milivolts, which can be measured and thus detect water leakage in embankment dams. This technique has been applied to the Santa Marta dam, and through a qualitative self potential data interpretation, of the measurements obtained in an upstream direction (underwater measurements), a vertical and sub horizontal water flows entering in the body dam were identified, which were causing internal erosion and developing a piping

Estimación del riesgo sísmico en Puerto Príncipe (Haití) como base para la elaboración de planes de emergencia

En este trabajo se realiza una propuesta de cálculo de riesgo sísmico en términos de daño físico para un escenario realista en Puerto Prínicpe (Haití). Los resultados del modelo propuesto se calibran con datos del sismo de enero de 2010.

Desarrollo de la aplicación web SISMO-HAITÍ para difusión de resultados de peligrosidad sísmica en la isla de La Española

Tras el catastrófico terremoto ocurrido en Haití el 12 de enero de 2010, de magnitud Mw 7 y profundidad de 10 km, (fuente: USGS) con un epicentro próximo a la capital, Puerto Príncipe (15 km), el país quedo en una situación catastrófica y de extrema pobreza, con necesidades básicas en salud, nutrición, educación y habitabilidad. Pocos meses después se inició el proyecto de cooperación SISMO-HAITI, financiado y coordinado por el Grupo de Investigación en Ingeniería Sísmica (GIIS) de la Universidad Politécnica de Madrid (UPM), con participación de otras universidades españolas y del CSIC y siendo la contraparte Haitiana el Observatorio de Vulnerabilidad y Medio Ambiente (ONEV). Uno de los objetivos del proyecto es el cálculo de peligrosidad sísmica en la Isla de La Española que constituya la base para la elaboración del primer código sísmico del país. El trabajo que aquí se presenta es una aplicación web desarrollada con el Sistema de Información Geográfica (SIG) del proyecto SISMO-HAITI. En esta aplicación se integran los diferentes mapas generados para el cálculo de la peligrosidad sísmica, así como los mapas resultantes, que pueden ser analizados e interpretados con mayor facilidad gracias a la aplicación. Para analizar la influencia de los diferentes inputs de cálculo se ha introducido el catálogo sísmico, las diferentes zonificaciones sismo genéticas y las principales fallas tectónicas. Toda esta información se puede superponer geográficamente con posibilidad de realizar consultas cruzadas en las correspondientes bases de datos, permitiendo el análisis de sensibilidad de éstos en los resultados. El desarrollo de esta aplicación web se ha creado a través de ArcGis Server 10

Mapa preliminar de períodos predominantes del suelo en Puerto Príncipe (Haití) a partir de medidas de ruido ambiental.

The 2010 Haiti earthquake, occurred on January 12th at 16:53:09 local time (21:53:09 UTC) with epicentral distance of 15 km from the capital Port au Prince, MW 7.0 and 13 km hypocenter deep, was the strongest event in the area since happened in 1770. The maximum macroseismic intensity was estimated as X (MMI scale). The aim of this research is to obtain a preliminary zonation of Port-au-Prince in terms of predominant resonance periods of ground. A total of 36 short-period ambient noise records have been carried out on a grid of about 500x500m. H/V spectral ratio method (HVSR) has been applied to determine the predominant period at each point. The lowest values (<0.2s) predominate in the southern part of the city, composed by Miocene conglomerates, while highest values (> 0.45s) correspond to the center and western parts, composed of Pleistocene-Holocene alluvial deposits and anthropogenic land reclaimed from the sea. We have determined the ground VS30 structure inside National Palace garden, using simultaneous ambient noise measurements. An array made up of 6 sensors were used, with 5 of them uniformly distributed along a circumference and a sixth one placed in its centre. The records were analyzed by using the spatial autocorrelation method (SPAC). The VS 30 value obtained was 331m/sec, in good agreement with the average values obtained for this area by other authors, using prospecting techniques.

Efectividad de trincheras y pantallas en el aislamiento de las vibraciones generadas por el tráfico ferroviario subterráneo

Se aborda un caso práctico de ferrocarril en un túnel soterrado en zona urbana con edificaciones tipo vivienda en las inmediaciones. Como procedimiento numérico para el desarrollo del estudio se ha adoptado el Método de los Elementos de Contorno, debido a las ventajas que presenta, frente a otros más ampliamente difundidos (Elementos Finitos, Diferencias Finitas, etc.), en el modelado de medios infinitos o semi-infinitos. Se ha definido un modelo bidimensional del terreno, en el que se incluye el túnel, con la fuente de excitación, y un edificio tipo. Las características del terreno se han obtenido mediante prospección geofísica. En base a un estudio de sensibilidad se ha simplificado el modelo, verificando su validez mediante comparación de resultados con registros de vibraciones recogidos in situ. Se ha llevado a cabo un estudio paramétrico tomando como referencia el modelo simplificado, en el que se incorporan trincheras así como pantallas de distintos materiales con diversas profundidades y en dos emplazamientos diferentes, obteniendo la reducción del nivel de vibración conseguido en los puntos más significativos del modelo para las distintas soluciones.

Application of Wignner-Ville distribution to identify anomalies in GPR profiles

Con base en la Distribución de Wigner-Ville(WVO) se realizó un análisis en tiempo y frecuencia de datos obtenidos con el Radar de Penetración Terrestre (GPR), basado en el estudio de la descomposición de la señal espectral. Se calcula una correlación entre la señal original y las componentes de tiempo-frecuencia para obtener anomalías estructurales de la información contenida en el radargrama relacionándola con la geología disponible. En primer lugar se describe la aplicación de un ejemplo teórico constituido por lo que representaría un túnel (tubería). Se obtuvieron las firmas correspondientes en el dominio del tiempo y en el dominio de la frecuencia. Finalmente se analiza esta metodología en un sido de prueba en la detección de un tambo enterrado donde son conocidas la geometría y su profundidad. Este especial sitio fue facilitado por la Universidad Nacional Autónoma de México, en los terrenos del Observatorio Magnético de Teoloyucan, Estado de México. Los resultados obtenidos son bastante alentadores, ya que la WVD es capaz de definir los rasgos morfofógicos relacionados con el tambo y abre la posibilidad de localizar este tipo de estructuras.

Nociones de reología

Esta publicación está destinada a facilitar la labor de los alumnos de la asignatura "Mecánica de los medios continuos" que el Profesor D. Luis Lozano Calvo explica en 4 º curso de la especialidad GEOFISICA, Sección Ciencias Físicas,de la Universidad Complutense. En la especialidad "Geofísica" se destina un cuatrimestre al estudio de los problemas de la mécánica de los sólidos. Tras un primer análisis de la elasticidad tradicional, se pasa revista a los fenómenos en que intervienen las características plásticas y el tiempo. Este libro tiene por finalidad facilitar la labor de los alumnos que deben estudiar estas materias, aliviándoles de la obligada consulta de textos foráneos. El lector experto no encontrará, pues, ninguna aportación original, y sí, en cambio, numerosas omisiones obligadas por el destino de lo escrito y la brevedad del periodo docente. En particular, ha sido especialmente dolorosa la supresión de los estupendos teoremas de Colonetti, de las doctrinas de Volterra, de los fenómenos de envejecimiento, de los problemas dinámicas, etc. etc. Realmente nos limitamos a apuntar los temas y a comunicar con el lenguaje necesario para que el lector pueda continuar en textos mejores el análisis del apasionante campo, aquí apenas desvelado. Respecto a la bibliografía, se ha procurado recoger las pocas publicaciones en castellano que existen, así como las que han servido para redactar los apuntes y aquellas que se estiman imprescindibles para estudios posteriores. Estas últimas han sido señaladas con un asterisco.

Hidrogeología e hidrodinámica del acuífero de los manantiales de Gormaz y su importancia en el caudal base del río Duero, España = Hydrogeology and hydrodynamic of the Gormaz Springs Aquifer and its importance to the base flow of the Duero River, Spain

Se cuantifican las descargas subterráneas de un acuífero a un río que lo atraviesa utilizando correlaciones estadísticas. El río Duero, España, incrementa su caudal base en varios m3/s, al atravesar unos afloramientos carbonatados mesozoicos en un pequeño tramo de su cabecera; esto es de especial importancia en época de estiaje, cuando la mayor parte del caudal base del río procede de manantiales que allí se sitúan. Dichos afloramientos corresponden a uno de los dos acuíferos calcáreos confinados, que se desarrollan en paralelo y están hidráulicamente desconectados por una capa impermeable, que forman el sistema acuífero de los manantiales de Gormaz. Este sistema se encuentra en estado de régimen natural y está apenas explotado. Se define el modelo conceptual de funcionamiento hidrogeológico, considerando el papel hidrogeológico de la falla de Gormaz, situada en la zona de descarga del sistema. Analizando información geológica antecedente y la geofísica exploratoria realizada, se obtuvo un mejor conocimiento de la geometría y los límites de los acuíferos, definiéndose un sistema acuífero con una zona de recarga en el sur, correspondiente a los afloramientos calcáreos, los cuales se confinan hacia el norte bajo el Terciario, hasta intersecar con la falla normal de Gormaz. El salto de falla genera una barrera para las formaciones permeables situadas al extremo norte (margen derecha del río Duero); a su vez, el plano de falla facilita el ascenso del agua subterránea del sistema acuífero en estudio y pone en conexión hidráulica los dos acuíferos. Se estimaron, además, los parámetros hidráulicos de los acuíferos en los alrededores de la falla. La buena correlación entre los niveles piezométricos y las descargas subterráneas al río Duero han permitido la reconstrucción del hidrograma de los manantiales de Gormaz en el periodo 1992-2006. Se calcula así que la contribución subterránea al río Duero es de 135.9 hm3/año, que supone el 18.9% de la aportación total del río. In a short stretch of its headwaters, the base flow of the River Duero increases by several m3/s as it traverses some Mesozoic carbonate outcrops. This is of special importance during the dry season, when the majority of the base flow of the river proceeds from springs in this reach. The outcrops correspond to one of two confined calcareous aquifers that developed in parallel but which are not hydraulically connected because of an impermeable layer. Together, they constitute the aquifer system of the Gormaz Springs. The system is still in its natural regime and is hardly exploited. This study defines the conceptual model of hydrogeological functioning, taking into consideration the role of the Gormaz Fault, which is situated in the discharge zone of the system. Analysis of both antecedent geological information and geophysical explorations has led to a better understanding of the geometry and boundaries of the aquifers, defining an aquifer system with a recharge zone in the south corresponding to in the calcareous outcrops. These calcareous outcrops are confined to the north below Tertiary formations, as far as their intersection with the normal fault of Gormaz. The throw of the fault forms the barrier of the permeable formations situated in the extreme north (right bank of the River Duero). In turn, the fault plane facilitates the upflow of groundwater from the aquifer system and creates hydraulic connection between the two aquifers. In addition, the study estimated the hydraulic parameters of the aquifer around the fault. The close correlation between piezometric levels and the groundwater discharges to the River Duero has enabled the reconstruction of the hydrogram of Gormaz springs over the period 1992-2006. By this means, it is calculated that the groundwater contribution to the River Duero is 135.9 hm3/year, or 18.9% of the total river inflow.

Comparación de técnicas geofísicas para determinación de contaminación de suelos agrícolas

El objetivo del presente trabajo es evaluar los métodos de prospección geofísica que, además de servir para la caracterización hidrológica superficial en una finca de investigación vitivinícola, alcance la mayor resolutividad para la localización de contaminantes agrícolas. Con este fin se ha realizado, en primer lugar una revisión teórica de los métodos de prospección geofísica que, a priori, pueden considerarse más adecuados por cuanto reflejan en mayor medida las variaciones físicas derivadas de los efluentes químicos agrícolas de baja salinidad. Para alcanzar el objetivo del proyecto se ha realizado una serie de perfiles eléctricos de Tomografía Eléctrica y Polarización Inducida, además de un Sondeo de Resonancia Magnética Nuclear Se han obtenido mediante las técnicas tomográficas una caracterización del subsuelo para la determinación de la litología en profundidad con variaciones laterales bruscas, así como la porosidad de estas formaciones. La sección de cargabilidad (PI) ha permitido la detección de iones disueltos en el agua subterránea. Combinando estos resultados, con los obtenidos por Resonancia Magnética en la localización de las masas de agua y la permeabilidad de las rocas circundantes, somos capaces de determinar los agentes contaminantes así como su movilidad y distribución. La conclusión es que el método que resulta más efectivo para detección de contaminantes sería la combinación de todos ellos, puesto que son complementarios entre sí. ABSTRACT The main goal of this project is to evaluate the methods of geophysical prospection, also serving for surface hydrologic characterization in a research farm wine, reaching the most resoluteness for locating agricultural pollutants. To this end was made, first a theoretical review of geophysical prospecting methods that may be considered more suitable, a priori, as they reflect more closely the physical changes resulting from agricultural chemical measure low salinity. Reaching the objective of the project has made a series of Electrical Tomography profiles (ET) and Induced Polarization profiles (IP) also has conducted a drilling of Nuclear Magnetic Resonance (RMN). It was obtained by tomography characterization techniques for determining subsurface litho logy in depth and abrupt lateral variations, as well as the porosity of these formations. The chargeability Section (IP) has enabled detection of ions dissolved in groundwater. Combining these results with those that was obtained by Nuclear Magnetic Resonance with the location of groundwater body and the surrounding rocks permeability, we are able to determine the pollutants and the mobility and distribution of contaminants. The bottom line is the method more effective for detection of contaminants would be the combination of all of them because they are complementary in the environmental study.

El análogo natural de almacenamiento y escape de Co2 de la cuenca de Gañuelas-Mazarrón : implicaciones para el comportamiento y la seguridad de un almacenamiento de Co2 en estaso supercrítico

El aumento de la temperatura media de la Tierra durante el pasado siglo en casi 1 ºC; la subida del nivel medio del mar; la disminución del volumen de hielo y nieve terrestres; la fuerte variabilidad del clima y los episodios climáticos extremos que se vienen sucediendo durante las ultimas décadas; y el aumento de las epidemias y enfermedades infecciosas son solo algunas de las evidencias del cambio climático actual, causado, principalmente, por la acumulación de gases de efecto invernadero en la atmósfera por actividades antropogénicas. La problemática y preocupación creciente surgida a raíz de estos fenómenos, motivo que, en 1997, se adoptara el denominado “Protocolo de Kyoto” (Japón), por el que los países firmantes adoptaron diferentes medidas destinadas a controlar y reducir las emisiones de los citados gases. Entre estas medidas cabe destacar las tecnologías CAC, enfocadas a la captura, transporte y almacenamiento de CO2. En este contexto se aprobó, en octubre de 2008, el Proyecto Singular Estratégico “Tecnologías avanzadas de generación, captura y almacenamiento de CO2” (PSE-120000-2008-6), cofinanciado por el Ministerio de Ciencia e Innovación y el FEDER, el cual abordaba, en su Subproyecto “Almacenamiento Geológico de CO2” (PSS-120000-2008-31), el estudio detallado, entre otros, del Análogo Natural de Almacenamiento y Escape de CO2 de la cuenca de Ganuelas-Mazarrón (Murcia). Es precisamente en el marco de dicho Proyecto en el que se ha realizado este trabajo, cuyo objetivo final ha sido el de predecir el comportamiento y evaluar la seguridad, a corto, medio y largo plazo, de un Almacenamiento Geológico Profundo de CO2 (AGP-CO2), mediante el estudio integral del citado análogo natural. Este estudio ha comprendido: i) la contextualización geológica e hidrogeológica de la cuenca, así como la investigación geofísica de la misma; ii) la toma de muestras de aguas de algunos acuíferos seleccionados con el fin de realizar su estudio hidrogeoquímico e isotópico; iii) la caracterización mineralógica, petrográfica, geoquímica e isotópica de los travertinos precipitados a partir de las aguas de algunos de los sondeos de la cuenca; y iv) la medida y caracterización química e isotópica de los gases libres y disueltos detectados en la cuenca, con especial atención al CO2 y 222Rn. Esta información, desarrollada en capítulos independientes, ha permitido realizar un modelo conceptual de funcionamiento del sistema natural que constituye la cuenca de Ganuelas-Mazarrón, así como establecer las analogías entre este y un AGP-CO2, con posibles escapes naturales y/o antropogénicos. La aplicación de toda esta información ha servido, por un lado, para predecir el comportamiento y evaluar la seguridad, a corto, medio y largo plazo, de un AGP-CO2 y, por otro, proponer una metodología general aplicable al estudio de posibles emplazamientos de AGP-CO2 desde la perspectiva de los reservorios naturales de CO2. Los resultados más importantes indican que la cuenca de Ganuelas-Mazarrón se trata de una cubeta o fosa tectónica delimitada por fallas normales, con importantes saltos verticales, que hunden al substrato rocoso (Complejo Nevado-Filabride), y rellenas, generalmente, por materiales volcánicos-subvolcánicos ácidos. Además, esta cuenca se encuentra rellena por formaciones menos resistivas que son, de muro a techo, las margas miocenas, predominantes y casi exclusivas de la cuenca, y los conglomerados y gravas pliocuaternarias. El acuífero salino profundo y enriquecido en CO2, puesto de manifiesto por la xx exploración geotérmica realizada en dicha cuenca durante la década de los 80 y objeto principal de este estudio, se encuentra a techo de los materiales del Complejo Nevado-Filabride, a una profundidad que podría superar los 800 m, según los datos de la investigación mediante sondeos y geofísica. Por ello, no se descarta la posibilidad de que el CO2 se encuentre en estado supe critico, por lo que la citada cuenca reuniría las características principales de un almacenamiento geológico natural y profundo de CO2, o análogo natural de un AGP-CO2 en un acuífero salino profundo. La sobreexplotación de los acuíferos mas someros de la cuenca, con fines agrícolas, origino, por el descenso de sus niveles piezométricos y de la presión hidrostática, el ascenso de las aguas profundas, salinas y enriquecidas en CO2, las cuales son las responsables de la contaminación de dichos acuíferos. El estudio hidrogeoquímico de las aguas de los acuíferos investigados muestra una gran variedad de hidrofacies, incluso en aquellos de litología similar. La alta salinidad de estas aguas las hace inservibles tanto para el consumo humano como para fines agrícolas. Además, el carácter ligeramente ácido de la mayoría de estas aguas determina que tengan gran capacidad para disolver y transportar, hacia la superficie, elementos pesados y/o tóxicos, entre los que destaca el U, elemento abundante en las rocas volcánicas ácidas de la cuenca, con contenidos de hasta 14 ppm, y en forma de uraninita submicroscópica. El estudio isotópico ha permitido discernir el origen, entre otros, del C del DIC de las aguas (δ13C-DIC), explicándose como una mezcla de dos componentes principales: uno, procedente de la descomposición térmica de las calizas y mármoles del substrato y, otro, de origen edáfico, sin descartar una aportación menor de C de origen mantélico. El estudio de los travertinos que se están formando a la salida de las aguas de algunos sondeos, por la desgasificación rápida de CO2 y el consiguiente aumento de pH, ha permitido destacar este fenómeno, por analogía, como alerta de escapes de CO2 desde un AGP-CO2. El análisis de los gases disueltos y libres, con especial atención al CO2 y al 222Rn asociado, indican que el C del CO2, tanto disuelto como en fase libre, tiene un origen similar al del DIC, confirmándose la menor contribución de CO2 de origen mantélico, dada la relación R/Ra del He existente en estos gases. El 222Rn sería el generado por el decaimiento radiactivo del U, particularmente abundante en las rocas volcánicas de la cuenca, y/o por el 226Ra procedente del U o del existente en los yesos mesinienses de la cuenca. Además, el CO2 actúa como carrier del 222Rn, hecho evidenciado en las anomalías positivas de ambos gases a ~ 1 m de profundidad y relacionadas principalmente con perturbaciones naturales (fallas y contactos) y antropogénicas (sondeos). La signatura isotópica del C a partir del DIC, de los carbonatos (travertinos), y del CO2 disuelto y libre, sugiere que esta señal puede usarse como un excelente trazador de los escapes de CO2 desde un AGPCO2, en el cual se inyectara un CO2 procedente, generalmente, de la combustión de combustibles fósiles, con un δ13C(V-PDB) de ~ -30 ‰. Estos resultados han permitido construir un modelo conceptual de funcionamiento del sistema natural de la cuenca de Ganuelas-Mazarrón como análogo natural de un AGP-CO2, y establecer las relaciones entre ambos. Así, las analogías mas importantes, en cuanto a los elementos del sistema, serian la existencia de: i) un acuífero salino profundo enriquecido en CO2, que seria análoga a la formación almacén de un AGPxxi CO2; ii) una formación sedimentaria margosa que, con una potencia superior a 500 m, se correspondería con la formación sello de un AGP-CO2; y iii) acuíferos mas someros con aguas dulces y aptas para el consumo humano, rocas volcánicas ricas en U y fallas que se encuentran selladas por yesos y/o margas; elementos que también podrían concurrir en un emplazamiento de un AGP-CO2. Por otro lado, los procesos análogos mas importantes identificados serian: i) la inyección ascendente del CO2, que seria análoga a la inyección de CO2 de origen antropogénico, pero este con una signatura isotópica δ13C(V-PDB) de ~ -30 ‰; ii) la disolución de CO2 y 222Rn en las aguas del acuífero profundo, lo que seria análogo a la disolución de dichos gases en la formación almacén de un AGP-CO2; iii) la contaminación de los acuíferos mas someros por el ascenso de las aguas sobresaturadas en CO2, proceso que seria análogo a la contaminación que se produciría en los acuíferos existentes por encima de un AGP-CO2, siempre que este se perturbara natural (reactivación de fallas) o artificialmente (sondeos); iv) la desgasificación (CO2 y gases asociados, entre los que destaca el 222Rn) del acuífero salino profundo a través de sondeos, proceso análogo al que pudiera ocurrir en un AGP-CO2 perturbado; y v) la formación rápida de travertinos, proceso análogo indicativo de que el AGP-CO2 ha perdido su estanqueidad. La identificación de las analogías más importantes ha permitido, además, analizar y evaluar, de manera aproximada, el comportamiento y la seguridad, a corto, medio y largo plazo, de un AGP-CO2 emplazado en un contexto geológico similar al sistema natural estudiado. Para ello se ha seguido la metodología basada en el análisis e identificación de los FEPs (Features, Events and Processes), los cuales se han combinado entre sí para generar y analizar diferentes escenarios de evolución del sistema (scenario analysis). Estos escenarios de evolución identificados en el sistema natural perturbado, relacionados con la perforación de sondeos, sobreexplotación de acuíferos, precipitación rápida de travertinos, etc., serian análogos a los que podrían ocurrir en un AGP-CO2 que también fuera perturbado antropogénicamente, por lo que resulta totalmente necesario evitar la perturbación artificial de la formación sello del AGPCO2. Por último, con toda la información obtenida se ha propuesto una metodología de estudio que pueda aplicarse al estudio de posibles emplazamientos de un AGP-CO2 desde la perspectiva de los reservorios naturales de CO2, sean estancos o no. Esta metodología comprende varias fases de estudio, que comprendería la caracterización geológico-estructural del sitio y de sus componentes (agua, roca y gases), la identificación de las analogías entre un sistema natural de almacenamiento de CO2 y un modelo conceptual de un AGP-CO2, y el establecimiento de las implicaciones para el comportamiento y la seguridad de un AGP-CO2. ABSTRACT The accumulation of the anthropogenic greenhouse gases in the atmosphere is the main responsible for: i) the increase in the average temperature of the Earth over the past century by almost 1 °C; ii) the rise in the mean sea level; iii) the drop of the ice volume and terrestrial snow; iv) the strong climate variability and extreme weather events that have been happening over the last decades; and v) the spread of epidemics and infectious diseases. All of these events are just some of the evidence of current climate change. The problems and growing concern related to these phenomena, prompted the adoption of the so-called "Kyoto Protocol" (Japan) in 1997, in which the signatory countries established different measurements to control and reduce the emissions of the greenhouse gases. These measurements include the CCS technologies, focused on the capture, transport and storage of CO2. Within this context, it was approved, in October 2008, the Strategic Singular Project "Tecnologías avanzadas de generación, captura y almacenamiento de CO2" (PSE-120000-2008-6), supported by the Ministry of Science and Innovation and the FEDER funds. This Project, by means of the Subproject "Geological Storage of CO2" (PSS- 120000-2008-31), was focused on the detailed study of the Natural Analogue of CO2 Storage and Leakage located in the Ganuelas-Mazarron Tertiary basin (Murcia), among other Spanish Natural Analogues. This research work has been performed in the framework of this Subproject, being its final objective to predict the behaviour and evaluate the safety, at short, medium and long-term, of a CO2 Deep Geological Storage (CO2-DGS) by means of a comprehensive study of the abovementioned Natural Analogue. This study comprises: i) the geological and hydrogeological context of the basin and its geophysical research; ii) the water sampling of the selected aquifers to establish their hydrogeochemical and isotopic features; iii) the mineralogical, petrographic, geochemical and isotopic characterisation of the travertines formed from upwelling groundwater of several hydrogeological and geothermal wells; and iv) the measurement of the free and dissolved gases detected in the basin, as well as their chemical and isotopic characterisation, mainly regarding CO2 and 222Rn. This information, summarised in separate chapters in the text, has enabled to build a conceptual model of the studied natural system and to establish the analogies between both the studied natural system and a CO2-DGS, with possible natural and/or anthropogenic escapes. All this information has served, firstly, to predict the behaviour and to evaluate the safety, at short, medium and long-term, of a CO2-DGS and, secondly, to propose a general methodology to study suitable sites for a CO2-DGS, taking into account the lessons learned from this CO2 natural reservoir. The main results indicate that the Ganuelas-Mazarron basin is a graben bounded by normal faults with significant vertical movements, which move down the metamorphic substrate (Nevado-Filabride Complex), and filled with acid volcanic-subvolcanic rocks. Furthermore, this basin is filled with two sedimentary formations: i) the Miocene marls, which are predominant and almost exclusive in the basin; xxiv and ii) the Plio-Quaternary conglomerates and gravels. A deep saline CO2-rich aquifer was evidenced in this basin as a result of the geothermal exploration wells performed during the 80s, located just at the top of the Nevado-Filabride Complex and at a depth that could exceed 800 m, according to the geophysical exploration performed. This saline CO2-rich aquifer is precisely the main object of this study. Therefore, it is not discarded the possibility that the CO2 in this aquifer be in supercritical state. Consequently, the aforementioned basin gathers the main characteristics of a natural and deep CO2 geological storage, or natural analogue of a CO2-DGS in a deep saline aquifer. The overexploitation of the shallow aquifers in this basin for agriculture purposes caused the drop of the groundwater levels and hydrostatic pressures, and, as a result, the ascent of the deep saline and CO2-rich groundwater, which is the responsible for the contamination of the shallow and fresh aquifers. The hydrogeochemical features of groundwater from the investigated aquifers show the presence of very different hydrofacies, even in those with similar lithology. The high salinity of this groundwater prevents the human and agricultural uses. In addition, the slightly acidic character of most of these waters determines their capacity to dissolve and transport towards the surface heavy and/or toxic elements, among which U is highlighted. This element is abundant in the acidic volcanic rocks of the basin, with concentrations up to 14 ppm, mainly as sub-microscopic uraninite crystals. The isotopic study of this groundwater, particularly the isotopic signature of C from DIC (δ13C-DIC), suggests that dissolved C can be explained considering a mixture of C from two main different sources: i) from the thermal decomposition of limestones and marbles forming the substrate; and ii) from edaphic origin. However, a minor contribution of C from mantle degassing cannot be discarded. The study of travertines being formed from upwelling groundwater of several hydrogeological and geothermal wells, as a result of the fast CO2 degassing and the pH increase, has allowed highlighting this phenomenon, by analogy, as an alert for the CO2 leakages from a CO2-DGS. The analysis of the dissolved and free gases, with special attention to CO2 and 222Rn, indicates that the C from the dissolved and free CO2 has a similar origin to that of the DIC. The R/Ra ratio of He corroborates the minor contribution of CO2 from the mantle degassing. Furthermore, 222Rn is generated by the radioactive decay of U, particularly abundant in the volcanic rocks of the basin, and/or by 226Ra from the U or from the Messinian gypsum in the basin. Moreover, CO2 acts as a carrier of the 222Rn, a fact evidenced by the positive anomalies of both gases at ~ 1 m depth and mainly related to natural (faults and contacts) and anthropogenic (wells) perturbations. The isotopic signature of C from DIC, carbonates (travertines), and dissolved and free CO2, suggests that this parameter can be used as an excellent tracer of CO2 escapes from a CO2-DGS, in which CO2 usually from the combustion of fossil fuels, with δ13C(V-PDB) of ~ -30 ‰, will be injected. All of these results have allowed to build a conceptual model of the behaviour of the natural system studied as a natural analogue of a CO2-DGS, as well as to establish the relationships between both natural xxv and artificial systems. Thus, the most important analogies, regarding the elements of the system, would be the presence of: i) a deep saline CO2-rich aquifer, which would be analogous to the storage formation of a CO2-DGS; ii) a marly sedimentary formation with a thickness greater than 500 m, which would correspond to the sealing formation of a CO2-DGS; and iii) shallow aquifers with fresh waters suitable for human consumption, U-rich volcanic rocks, and faults that are sealed by gypsums and/or marls; geological elements that could also be present in a CO2-DGS. On the other hand, the most important analogous processes identified are: i) the upward injection of CO2, which would be analogous to the downward injection of the anthropogenic CO2, this last with a δ13C(V-PDB) of ~ -30 ‰; ii) the dissolution of CO2 and 222Rn in groundwater of the deep aquifer, which would be analogous to the dissolution of these gases in the storage formation of a CO2-DGS; iii) the contamination of the shallow aquifers by the uprising of CO2-oversaturated groundwater, an analogous process to the contamination that would occur in shallow aquifers located above a CO2-DGS, whenever it was naturally (reactivation of faults) or artificially (wells) perturbed; iv) the degassing (CO2 and associated gases, among which 222Rn is remarkable) of the deep saline aquifer through wells, process which could be similar in a perturbed CO2- DGS; v) the rapid formation of travertines, indicating that the CO2-DGS has lost its seal capacity. The identification of the most important analogies has also allowed analysing and evaluating, approximately, the behaviour and safety in the short, medium and long term, of a CO2-DGS hosted in a similar geological context of the natural system studied. For that, it has been followed the methodology based on the analysis and identification of FEPs (Features, Events and Processes) that have been combined together in order to generate and analyse different scenarios of the system evolution (scenario analysis). These identified scenarios in the perturbed natural system, related to boreholes, overexploitation of aquifers, rapid precipitation of travertines, etc., would be similar to those that might occur in a CO2-DGS anthropogenically perturbed, so that it is absolutely necessary to avoid the artificial perturbation of the seal formation of a CO2-DGS. Finally, a useful methodology for the study of possible sites for a CO2-DGS is suggested based on the information obtained from this investigation, taking into account the lessons learned from this CO2 natural reservoir. This methodology comprises several phases of study, including the geological and structural characterisation of the site and its components (water, rock and gases), the identification of the analogies between a CO2 storage natural system and a conceptual model of a CO2-DGS, and the implications regarding the behaviour and safety of a CO2-DGS.

Evaluación de la protección conseguida en la cuenca de un torrente de montaña ante los daños causados por eventos torrenciales mediante su restauración hidrológico-forestal. Cuenca del torrente de Arratiecho en el Pirineo Aragonés, España.

La demanda continua e insistente de innovación en el ámbito cotidiano, nos hace perder la perspectiva de análisis cuando los procesos requieren un periodo prolongado para resolverse. Esto ocurre con los bosques ubicados en cuencas vertientes a torrentes de montaña, que tienen su origen en la restauración hidrológico-forestal llevada a cabo en las mismas, como complemento de los trabajos de corrección hidráulica efectuados en los torrentes, para controlar los procesos de erosión, inundaciones y aterramientos, desencadenados en diferentes zonas de la cuenca, cuando se dinamizan en ella los ciclos del agua y de los sedimentos a causa de los aguaceros o de la repentina fusión del manto de nieve. La cuenca vertiente al torrente de Arratiecho, situada en el municipio de Biescas, Pirineo Aragonés, fue restaurada y el barranco corregido entre 1903-05; para evitar su desertificación ante la erosión generalizada que sufría la misma, frenar el deslizamiento de una ladera morrénica y reducir los efectos de las inundaciones y aterramientos en su cono de sedimentación, donde se ubican los terrenos más fértiles y le atraviesa la carretera nacional N-260. Se comentan los trabajos realizados y los resultados obtenidos un siglo después, atendiendo especialmente a la evolución de las reforestaciones, siguiendo los criterios utilizados en el análisis de los bosques protectores.