Origin of the gases released from the Acqua Passante and Ermeta wells (Mt. Amiata, central Italy) and possible environmental implications for their closure

The Mt. Amiata volcano (Tuscany, central Italy) hosts the second largest geothermal field of Italy. Its SW and NE sectors are characterized by the presence of several CO2-rich (mayor que95% by vol.) gas discharges. An intense Hg mining activity had taken place from the 19th century up to the end of the ?70s, particularly close to Abbadia San Salvatore, during which two drillings (Acqua Passante and Ermeta) intercepted a CO2-rich gas fertile horizon. The related gases are emitted in the atmosphere since 1938 and 1959, respectively, causing severe concerns for the local air quality. In this work the results of a geochemical and isotopic survey carried out on these gas emissions from March 2009 to January 2014 are presented. CO2 fluxes from both the two wells and soil from an area of about 653,500 m2 located between them were measured. The two wells are emitting up to 15,000, 92 and 8 tons y-1 of CO2, CH4 and H2S, respectively, while the computed soil CO2 output was estimated at 4,311 ton y-1. The spatial distribution of the CO2 soil flux suggests the presence of preferential patterns, indicating sites of higher permeability. Since the local municipality is evaluating the possibility to plug the Ermeta vent, a temporarily closure should first be carried out to test the possible influence of this operation on the diffuse soil degassing of deep-originated CO2 in the surrounding area. This implies that diffuse soil gases should carefully be monitored before proceeding with its definitive closure.

Hydrogeochemistry of surface and spring waters in the surroundings of the CO2 injection site at Hontomín–Huermeces (Burgos, Spain)

In this paper the very first geochemical and isotopic data related to surface and spring waters and dissolved gases in the area of Hontomín–Huermeces (Burgos, Spain) are presented and discussed. Hontomín–Huermeces has been selected as a pilot site for the injection of pure (>99%) CO2. Injection and monitoring wells are planned to be drilled close to 6 oil wells completed in the 1980s for which detailed stratigraphical logs are available, indicating the presence of a confined saline aquifer at the depth of about 1500 m into which less than 100,000 tons of iquid CO2 will be injected, possibly starting in 2013. The chemical and features of the spring waters suggest that they are related to a shallow hydrogeological system as the concentration of the Total Dissolved Solids approaches 800 mg/L with a Ca2+(Mg2+)-HCO3− composition, similar to that of the surface waters. This is also supported by the oxygen and hydrogen isotopic ratios that have values lying between those of the Global and the Mediterranean Meteoric Water Lines. Some spring waters close to the oil wells are haracterized by relatively high concentrations of NO3− (up to 123 mg/L), unequivocally suggesting an anthropogenic source that adds to the main water–rock interaction processes. The latter can be referred to Ca-Mg-carbonate and, at a minor extent, Al-silicate dissolution, being the outcropping sedimentary rocks characterized by Palaeozoic to Quaternary rocks. Anomalous concentrations of Cl−, SO42−, As, B and Ba were measured in two springs discharging a few hundred meters from the oil wells and in the Rio Ubierna. These contents are significantly higher than those of the whole set of the studied waters and are possibly indicative of mixing processes, although at very low extent, between deep and shallow aquifers. No evidence of deep-seated gases interacting with the Hontomín–Huermeces waters was recognized in the chemistry of the disolved gases. This is likely due to the fact that they are mainly characterized by an atmospheric source as highlighted by the high contents of N2, O2 and Ar and by N2/Ar ratios that approach that of ASW (Air Saturated Water) and possibly masking any contribution related to a deep source. Nevertheless, significant concentrations (up to 63% by vol.) of isotopically negative CO2 (<−17.7‰ V-PDB) were found in some water samples, likely related to a biogenic source. The geochemical and isotopic data of this work are of particular importance when a monitoring program will be established to verify whether CO2 leakages, induced by the injection of this greenhouse gas, may be affecting the quality of the waters in the shallow hydrological circuits at Hontomín–Huermeces. In this respect, carbonate chemistry, the isotopic carbon of dissolved CO2 and TDIC (Total Dissolved Inorganic Carbon) and selected trace elements can be considered as useful parameters to trace the migration of the injected CO2 into near-surface environments.

Water chemistry at Hontomín-Huermeces (Burgos, Spain): insights for a pre-, intra- and post-CO2 injection geochemical monitoring.

In this study, the very first geochemical and isotopic data related to surface and spring waters and dissolved gases in the area of Hontomín-Huermeces (Burgos, Spain) are presented and discussed. Hontomín-Huermeces was selected as a pilot site for the injection of pure (>99 %) CO2. Injection and monitoring wells are planned to be drilled close to 6 oil wells completed in the 1980’s. Stratigraphical logs indicate the presence of a confined saline aquifer at the depth of about 1,500 m into which less than 100,000 tons of liquid CO2 will be injected, possibly starting in 2013. The chemical and isotopic features of the spring waters suggest the occurrence of a shallow aquifer having a Ca2+(Mg2+)-HCO3- composition, relatively low salinity (Total Dissolved Solids _800 mg/L) and a meteoric isotopic signature. Some spring waters close to the oil wells are characterized by relatively high concentrations of NO3- (up to 123 mg/L), unequivocally indicating anthropogenic contamination that adds to the main water-rock interaction processes. The latter can be referred to Ca-Mg-carbonate and, at a minor extent, Al-silicate dissolution, being the outcropping sedimentary rocks characterized by Palaeozoic to Quaternary rocks. Anomalous concentrations of Cl-, SO42-, As, B and Ba were measured in two springs discharging a few hundreds meters from the oil wells and in the Rio Ubierna, possibly indicative of mixing processes, although at very low extent, between deep and shallow aquifers. Gases dissolved in spring waters show relatively high concentrations of atmospheric species, such as N2, O2 and Ar, and isotopically negative CO2 (<-17.7 h V-PDB), likely related to a biogenic source, possibly masking any contribution related to a deep source. The geochemical and isotopic data of this study are of particular importance when a monitoring program will be established to verify whether CO2 leakages, induced by the injection of this greenhouse gas, may affect the quality of the waters of the shallow Hontomín-Huermeces hydrological circuit. In this respect, carbonate chemistry, the isotopic carbon of dissolved CO2 and TDIC (Total Dissolved Inorganic Carbon) and selected trace elements can be considered as useful parameters to trace the migration of the injected CO2 into near-surface environments.

Estrategias de monitorización de CO2 y otros gases en los estudios de análogos naturales

La Universidad Politécnica de Madrid (UPM) y la Università degli Studi di Firenze (UniFi), bajo la coordinación técnica de AMPHOS21, participan desde 2009 en el proyecto de investigación “Estrategias de Monitorización de CO2 y otros gases en el estudio de Análogos Naturales”, financiado por la Fundación Ciudad de la Energía (CIUDEN) en el marco del Proyecto Compostilla OXYCFB300 (http://www.compostillaproject.eu), del Programa “European Energy Program for Recovery - EEPR”. El objetivo principal del proyecto fue el desarrollo y puesta a punto de metodologías de monitorización superficiales para su aplicación en el seguimiento y control de los emplazamientos donde se realice el almacenamiento geológico de CO2, analizando técnicas que permitan detectar y cuantificar las posibles fugas de CO2 a la atmósfera. Los trabajos se realizaron tanto en análogos naturales (españoles e italianos) como en la Planta de Desarrollo Tecnológico de Almacenamiento de CO2 de Hontomín. Las técnicas analizadas se centran en la medición de gases y aguas superficiales (de escorrentía y manantiales). En cuanto a la medición de gases se analizó el flujo de CO2 que emana desde el suelo a la atmósfera y la aplicabilidad de trazadores naturales (como el radón) para la detección e identificación de las fugas de CO2. En cuanto al análisis químico de las aguas se analizaron los datos geoquímicos e isotópicos y los gases disueltos en las aguas de los alrededores de la PDT de Hontomín, con objeto de determinar qué parámetros son los más apropiados para la detección de una posible migración del CO2 inyectado, o de la salmuera, a los ambientes superficiales. Las medidas de flujo de CO2 se realizaron con la técnica de la cámara de acúmulo. A pesar de ser una técnica desarrollada y aplicada en diferentes ámbitos científicos se estimó necesario adaptar un protocolo de medida y de análisis de datos a las características específicas de los proyectos de captura y almacenamiento de CO2 (CAC). Donde los flujos de CO2 esperados son bajos y en caso de producirse una fuga habrá que detectar pequeñas variaciones en los valores flujo con un “ruido” en la señal alto, debido a actividad biológica en el suelo. La medida de flujo de CO2 mediante la técnica de la cámara de acúmulo se puede realizar sin limpiar la superficie donde se coloca la cámara o limpiando y esperando al reequilibrio del flujo después de la distorsión al sistema. Sin embargo, los resultados obtenidos después de limpiar y esperar muestran menor dispersión, lo que nos indica que este procedimiento es el mejor para la monitorización de los complejos de almacenamiento geológico de CO2. El protocolo de medida resultante, utilizado para la obtención de la línea base de flujo de CO2 en Hontomín, sigue los siguiente pasos: a) con una espátula se prepara el punto de medición limpiando y retirando el recubrimiento vegetal o la primera capa compacta de suelo, b) se espera un tiempo para la realización de la medida de flujo, facilitando el reequilibrio del flujo del gas tras la alteración provocada en el suelo y c) se realiza la medida de flujo de CO2. Una vez realizada la medición de flujo de CO2, y detectada si existen zonas de anomalías, se debe estimar la cantidad de CO2 que se está escapando a la atmósfera (emanación total), con el objetivo de cuantificar la posible fuga. Existen un amplio rango de metodologías para realizar dicha estimación, siendo necesario entender cuáles son las más apropiadas para obtener el valor más representativo del sistema. En esta tesis se comparan seis técnicas estadísticas: media aritmética, estimador insegado de la media (aplicando la función de Sichel), remuestreo con reemplazamiento (bootstrap), separación en diferentes poblaciones mediante métodos gráficos y métodos basados en criterios de máxima verosimilitud, y la simulación Gaussiana secuencial. Para este análisis se realizaron ocho campañas de muestreo, tanto en la Planta de Desarrollo Tecnológico de Hontomón como en análogos naturales (italianos y españoles). Los resultados muestran que la simulación Gaussiana secuencial suele ser el método más preciso para realizar el cálculo, sin embargo, existen ocasiones donde otros métodos son más apropiados. Como consecuencia, se desarrolla un procedimiento de actuación para seleccionar el método que proporcione el mejor estimador. Este procedimiento consiste, en primer lugar, en realizar un análisis variográfico. Si existe una autocorrelación entre los datos, modelizada mediante el variograma, la mejor técnica para calcular la emanación total y su intervalo de confianza es la simulación Gaussiana secuencial (sGs). Si los datos son independientes se debe comprobar la distribución muestral, aplicando la media aritmética o el estimador insesgado de la media (Sichel) para datos normales o lognormales respectivamente. Cuando los datos no son normales o corresponden a una mezcla de poblaciones la mejor técnica de estimación es la de remuestreo con reemplazamiento (bootstrap). Siguiendo este procedimiento el máximo valor del intervalo de confianza estuvo en el orden del ±20/25%, con la mayoría de valores comprendidos entre ±3,5% y ±8%. La identificación de las diferentes poblaciones muestrales en los datos de flujo de CO2 puede ayudar a interpretar los resultados obtenidos, toda vez que esta distribución se ve afectada por la presencia de varios procesos geoquímicos como, por ejemplo, una fuente geológica o biológica del CO2. Así pues, este análisis puede ser una herramienta útil en el programa de monitorización, donde el principal objetivo es demostrar que no hay fugas desde el reservorio a la atmósfera y, si ocurren, detectarlas y cuantificarlas. Los resultados obtenidos muestran que el mejor proceso para realizar la separación de poblaciones está basado en criterios de máxima verosimilitud. Los procedimientos gráficos, aunque existen pautas para realizarlos, tienen un cierto grado de subjetividad en la interpretación de manera que los resultados son menos reproducibles. Durante el desarrollo de la tesis se analizó, en análogos naturales, la relación existente entre el CO2 y los isótopos del radón (222Rn y 220Rn), detectándose en todas las zonas de emisión de CO2 una relación positiva entre los valores de concentración de 222Rn en aire del suelo y el flujo de CO2. Comparando la concentración de 220Rn con el flujo de CO2 la relación no es tan clara, mientras que en algunos casos aumenta en otros se detecta una disminución, hecho que parece estar relacionado con la profundidad de origen del radón. Estos resultados confirmarían la posible aplicación de los isótopos del radón como trazadores del origen de los gases y su aplicación en la detección de fugas. Con respecto a la determinación de la línea base de flujo CO2 en la PDT de Hontomín, se realizaron mediciones con la cámara de acúmulo en las proximidades de los sondeos petrolíferos, perforados en los ochenta y denominados H-1, H-2, H-3 y H-4, en la zona donde se instalarán el sondeo de inyección (H-I) y el de monitorización (H-A) y en las proximidades de la falla sur. Desde noviembre de 2009 a abril de 2011 se realizaron siete campañas de muestreo, adquiriéndose más de 4.000 registros de flujo de CO2 con los que se determinó la línea base y su variación estacional. Los valores obtenidos fueron bajos (valores medios entre 5 y 13 g•m-2•d-1), detectándose pocos valores anómalos, principalmente en las proximidades del sondeo H-2. Sin embargo, estos valores no se pudieron asociar a una fuente profunda del CO2 y seguramente estuvieran más relacionados con procesos biológicos, como la respiración del suelo. No se detectaron valores anómalos cerca del sistema de fracturación (falla Ubierna), toda vez que en esta zona los valores de flujo son tan bajos como en el resto de puntos de muestreo. En este sentido, los valores de flujo de CO2 aparentemente están controlados por la actividad biológica, corroborado al obtenerse los menores valores durante los meses de otoño-invierno e ir aumentando en los periodos cálidos. Se calcularon dos grupos de valores de referencia, el primer grupo (UCL50) es 5 g•m-2•d-1 en las zonas no aradas en los meses de otoño-invierno y 3,5 y 12 g•m-2•d-1 en primavera-verano para zonas aradas y no aradas, respectivamente. El segundo grupo (UCL99) corresponde a 26 g•m-2•d- 1 durante los meses de otoño-invierno en las zonas no aradas y 34 y 42 g•m-2•d-1 para los meses de primavera-verano en zonas aradas y no aradas, respectivamente. Flujos mayores a estos valores de referencia podrían ser indicativos de una posible fuga durante la inyección y posterior a la misma. Los primeros datos geoquímicos e isotópicos de las aguas superficiales (de escorrentía y de manantiales) en el área de Hontomín–Huermeces fueron analizados. Los datos sugieren que las aguas estudiadas están relacionadas con aguas meteóricas con un circuito hidrogeológico superficial, caracterizadas por valores de TDS relativamente bajos (menor a 800 mg/L) y una fácie hidrogeoquímica de Ca2+(Mg2+)-HCO3 −. Algunas aguas de manantiales se caracterizan por concentraciones elevadas de NO3 − (concentraciones de hasta 123 mg/l), lo que sugiere una contaminación antropogénica. Se obtuvieron concentraciones anómalas de of Cl−, SO4 2−, As, B y Ba en dos manantiales cercanos a los sondeos petrolíferos y en el rio Ubierna, estos componentes son probablemente indicadores de una posible mezcla entre los acuíferos profundos y superficiales. El estudio de los gases disueltos en las aguas también evidencia el circuito superficial de las aguas. Estando, por lo general, dominado por la componente atmosférica (N2, O2 y Ar). Sin embargo, en algunos casos el gas predominante fue el CO2 (con concentraciones que llegan al 63% v/v), aunque los valores isotópicos del carbono (<-17,7 ‰) muestran que lo más probable es que esté relacionado con un origen biológico. Los datos geoquímicos e isotópicos de las aguas superficiales obtenidos en la zona de Hontomín se pueden considerar como el valor de fondo con el que comparar durante la fase operacional, la clausura y posterior a la clausura. En este sentido, la composición de los elementos mayoritarios y traza, la composición isotópica del carbono del CO2 disuelto y del TDIC (Carbono inorgánico disuelto) y algunos elementos traza se pueden considerar como parámetros adecuados para detectar la migración del CO2 a los ambientes superficiales. ABSTRACT Since 2009, a group made up of Universidad Politécnica de Madrid (UPM; Spain) and Università degli Studi Firenze (UniFi; Italy) has been taking part in a joint project called “Strategies for Monitoring CO2 and other Gases in Natural analogues”. The group was coordinated by AMPHOS XXI, a private company established in Barcelona. The Project was financially supported by Fundación Ciudad de la Energía (CIUDEN; Spain) as a part of the EC-funded OXYCFB300 project (European Energy Program for Recovery -EEPR-; www.compostillaproject.eu). The main objectives of the project were aimed to develop and optimize analytical methodologies to be applied at the surface to Monitor and Verify the feasibility of geologically stored carbon dioxide. These techniques were oriented to detect and quantify possible CO2 leakages to the atmosphere. Several investigations were made in natural analogues from Spain and Italy and in the Tecnchnological Development Plant for CO2 injection al Hontomín (Burgos, Spain). The studying techniques were mainly focused on the measurements of diffuse soil gases and surface and shallow waters. The soil-gas measurements included the determination of CO2 flux and the application to natural trace gases (e.g. radon) that may help to detect any CO2 leakage. As far as the water chemistry is concerned, geochemical and isotopic data related to surface and spring waters and dissolved gases in the area of the PDT of Hontomín were analyzed to determine the most suitable parameters to trace the migration of the injected CO2 into the near-surface environments. The accumulation chamber method was used to measure the diffuse emission of CO2 at the soil-atmosphere interface. Although this technique has widely been applied in different scientific areas, it was considered of the utmost importance to adapt the optimum methodology for measuring the CO2 soil flux and estimating the total CO2 output to the specific features of the site where CO2 is to be stored shortly. During the pre-injection phase CO2 fluxes are expected to be relatively low where in the intra- and post-injection phases, if leakages are to be occurring, small variation in CO2 flux might be detected when the CO2 “noise” is overcoming the biological activity of the soil (soil respiration). CO2 flux measurements by the accumulation chamber method could be performed without vegetation clearance or after vegetation clearance. However, the results obtained after clearance show less dispersion and this suggests that this procedure appears to be more suitable for monitoring CO2 Storage sites. The measurement protocol, applied for the determination of the CO2 flux baseline at Hontomín, has included the following steps: a) cleaning and removal of both the vegetal cover and top 2 cm of soil, b) waiting to reduce flux perturbation due to the soil removal and c) measuring the CO2 flux. Once completing the CO2 flux measurements and detected whether there were anomalies zones, the total CO2 output was estimated to quantify the amount of CO2 released to the atmosphere in each of the studied areas. There is a wide range of methodologies for the estimation of the CO2 output, which were applied to understand which one was the most representative. In this study six statistical methods are presented: arithmetic mean, minimum variances unbiased estimator, bootstrap resample, partitioning of data into different populations with a graphical and a maximum likelihood procedures, and sequential Gaussian simulation. Eight campaigns were carried out in the Hontomín CO2 Storage Technology Development Plant and in natural CO2 analogues. The results show that sequential Gaussian simulation is the most accurate method to estimate the total CO2 output and the confidential interval. Nevertheless, a variety of statistic methods were also used. As a consequence, an application procedure for selecting the most realistic method was developed. The first step to estimate the total emanation rate was the variogram analysis. If the relation among the data can be explained with the variogram, the best technique to calculate the total CO2 output and its confidence interval is the sequential Gaussian simulation method (sGs). If the data are independent, their distribution is to be analyzed. For normal and log-normal distribution the proper methods are the arithmetic mean and minimum variances unbiased estimator, respectively. If the data are not normal (log-normal) or are a mixture of different populations the best approach is the bootstrap resampling. According to these steps, the maximum confidence interval was about ±20/25%, with most of values between ±3.5% and ±8%. Partitioning of CO2 flux data into different populations may help to interpret the data as their distribution can be affected by different geochemical processes, e.g. geological or biological sources of CO2. Consequently, it may be an important tool in a monitoring CCS program, where the main goal is to demonstrate that there are not leakages from the reservoir to the atmosphere and, if occurring, to be able to detect and quantify it. Results show that the partitioning of populations is better performed by maximum likelihood criteria, since graphical procedures have a degree of subjectivity in the interpretation and results may not be reproducible. The relationship between CO2 flux and radon isotopes (222Rn and 220Rn) was studied in natural analogues. In all emissions zones, a positive relation between 222Rn and CO2 was observed. However, the relationship between activity of 220Rn and CO2 flux is not clear. In some cases the 220Rn activity indeed increased with the CO2 flux in other measurements a decrease was recognized. We can speculate that this effect was possibly related to the route (deep or shallow) of the radon source. These results may confirm the possible use of the radon isotopes as tracers for the gas origin and their application in the detection of leakages. With respect to the CO2 flux baseline at the TDP of Hontomín, soil flux measurements in the vicinity of oil boreholes, drilled in the eighties and named H-1 to H-4, and injection and monitoring wells were performed using an accumulation chamber. Seven surveys were carried out from November 2009 to summer 2011. More than 4,000 measurements were used to determine the baseline flux of CO2 and its seasonal variations. The measured values were relatively low (from 5 to 13 g•m-2•day-1) and few outliers were identified, mainly located close to the H-2 oil well. Nevertheless, these values cannot be associated to a deep source of CO2, being more likely related to biological processes, i.e. soil respiration. No anomalies were recognized close to the deep fault system (Ubierna Fault) detected by geophysical investigations. There, the CO2 flux is indeed as low as other measurement stations. CO2 fluxes appear to be controlled by the biological activity since the lowest values were recorded during autumn-winter seasons and they tend to increase in warm periods. Two reference CO2 flux values (UCL50 of 5 g•m-2•d-1 for non-ploughed areas in autumn-winter seasons and 3.5 and 12 g•m-2•d-1 for in ploughed and non-ploughed areas, respectively, in spring-summer time, and UCL99 of 26 g•m-2•d-1 for autumn-winter in not-ploughed areas and 34 and 42 g•m-2•d-1 for spring-summer in ploughed and not-ploughed areas, respectively, were calculated. Fluxes higher than these reference values could be indicative of possible leakage during the operational and post-closure stages of the storage project. The first geochemical and isotopic data related to surface and spring waters and dissolved gases in the area of Hontomín–Huermeces (Burgos, Spain) are presented and discussed. The chemical and features of the spring waters suggest that they are related to a shallow hydrogeological system as the concentration of the Total Dissolved Solids approaches 800 mg/L with a Ca2+(Mg2+)-HCO3 − composition, similar to that of the surface waters. Some spring waters are characterized by relatively high concentrations of NO3 − (up to 123 mg/L), unequivocally suggesting an anthropogenic source. Anomalous concentrations of Cl−, SO4 2−, As, B and Ba were measured in two springs, discharging a few hundred meters from the oil wells, and in the Rio Ubierna. These contents are possibly indicative of mixing processes between deep and shallow aquifers. The chemistry of the dissolved gases also evidences the shallow circuits of the Hontomín– Huermeces, mainly characterized by an atmospheric source as highlighted by the contents of N2, O2, Ar and their relative ratios. Nevertheless, significant concentrations (up to 63% by vol.) of isotopically negative CO2 (<−17.7‰ V-PDB) were found in some water samples, likely related to a biogenic source. The geochemical and isotopic data of the surface and spring waters in the surroundings of Hontomín can be considered as background values when intra- and post-injection monitoring programs will be carried out. In this respect, main and minor solutes, the isotopic carbon of dissolved CO2 and TDIC (Total Dissolved Inorganic Carbon) and selected trace elements can be considered as useful parameters to trace the migration of the injected CO2 into near-surface environments.

Virtual-bound, filamentary and layered states in a box-shaped quantum dot of square potential form the exact numerical solution of the effective mass Schrodinger equation

The effective mass Schrodinger equation of a QD of parallelepipedic shape with a square potential well is solved by diagonalizing the exact Hamiltonian matrix developed in a basis of separation-of-variables wavefunctions. The expected below bandgap bound states are found not to differ very much from the former approximate calculations. In addition, the presence of bound states within the conduction band is confirmed. Furthermore, filamentary states bounded in two dimensions and extended in one dimension and layered states with only one dimension bounded, all within the conduction band which are similar to those originated in quantum wires and quantum wells coexist with the ordinary continuum spectrum of plane waves. All these subtleties are absent in spherically shaped quantum dots, often used for modeling.

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.

Establishing Quercus ilex under Mediterranean dry conditions: sowing recalcitrant acorns versus planting seedlings at different depths and tube shelter light transmissions

Article New Forests November 2015, Volume 46, Issue 5, pp 869-883 First online: 17 June 2015 Establishing Quercus ilex under Mediterranean dry conditions: sowing recalcitrant acorns versus planting seedlings at different depths and tube shelter light transmissionsJuan A. OlietAffiliated withDepartamento de Sistemas y Recursos Naturales, E.T.S. Ingenieros de Montes, Universidad Politécnica de Madrid Email author View author's OrcID profile , Alberto Vázquez de CastroAffiliated withDepartamento de Sistemas y Recursos Naturales, E.T.S. Ingenieros de Montes, Universidad Politécnica de Madrid, Jaime PuértolasAffiliated withLancaster Environment Centre, Lancaster University $39.95 / €34.95 / £29.95 * Rent the article at a discount Rent now * Final gross prices may vary according to local VAT. Get Access AbstractSuccess of Mediterranean dry areas restoration with oaks is a challenging goal. Testing eco-techniques that mimic beneficial effects of natural structures and ameliorate stress contributes to positive solutions to overcoming establishment barriers. We ran a factorial experiment in a dry area, testing two levels of solid wall transmission of tube shelters (60 and 80 %) plus a control mesh, and two depths (shallow and 15 cm depth) of placing either planted seedlings or acorns of Quercus ilex. Microclimate of the planting or sowing spots was characterized by measuring photosynthetically active radiation, temperature and relative humidity. Plant response was evaluated in terms of survival, phenology, acorn emergence and photochemical efficiency (measured through chlorophyll fluorescence). We hypothesize that tube shelters and deep planting improve Q. ilex post-planting and sowing performance because of the combined effects of reducing excessive radiation and improving access to moist soil horizons. Results show that temperature and PAR was reduced, and relative humidity increased, in deep spots. Midsummer photochemical efficiency indicates highest level of stress for oaks in 80 % light transmission shelter. Optimum acorn emergence in spring was registered within solid wall tree shelters, and maximum summer survival of germinants and of planted seedlings occurred when acorns or seedlings were placed at 15 cm depth irrespectively of light transmission of shelter. Survival of germinants was similar to that of planted seedlings. The importance of techniques to keep high levels of viability after sowing recalcitrant seeds in the field is emphasized in the study

CO2 soil flux baseline at the technological development plant for CO2 injection at Hontomin (Burgos, Spain)

From the end of 2013 and during the following two years, 20 kt of CO2sc are planned to be injected in a saline reservoir (1500 m depth) at the Hontomín site (NE Spain). The target aquifers are Lower Jurassic limestone formations which are sealed by Lower Cretaceous clay units at the Hontomín site (NE Spain). The injection of CO2 is part of the activities committed in the Technology Development phase of the EC-funded OXYCFB300 project (European Energy Program for Recovery – EEPR, http://www.compostillaproject.eu), which include CO2 injection strategies, risk assessment, and testing and validating monitoring methodologies and techniques. Among the monitoring works, the project is intended to prove that present-day technology is able to monitor the evolution of injected CO2 in the reservoir and to detect potential leakage. One of the techniques is the measurement of CO2 flux at the soil–atmosphere interface, which includes campaigns before, during and after the injection operations. In this work soil CO2 flux measurements in the vicinity of oil borehole, drilled in the eighties and named H-1 to H-4, and injection and monitoring wells were performed using an accumulation chamber equipped with an IR sensor. Seven surveys were carried out from November 2009 to summer 2011. More than 4000 measurements were used to determine the baseline flux of CO2 and its seasonal variations. The measured values were low (from 5 to 13 g m−2 day−1) and few outliers were identified, mainly located close to the H-2 oil well. Nevertheless, these values cannot be associated to a deep source of CO2, being more likely related to biological processes, i.e. soil respiration. No anomalies were recognized close to the deep fault system (Ubierna Fault) detected by geophysical investigations. There, the CO2 flux is indeed as low as other measurement stations. CO2 fluxes appear to be controlled by the biological activity since the lowest values were recorded during autumn-winter seasons and they tend to increase in warm periods. Two reference CO2 flux values (UCL50 of 5 g m−2 d−1 for non-ploughed areas in autumn–winter seasons and 3.5 and 12 g m−2 d−1 for in ploughed and non-ploughed areas, respectively, in spring–summer time, and UCL99 of 26 g m−2 d−1 for autumn–winter in not-ploughed areas and 34 and 42 g m−2 d−1 for spring–summer in ploughed and not-ploughed areas, respectively) were calculated. Fluxes higher than these reference values could be indicative of possible leakage during the operational and post-closure stages of the storage project.

CO2and Rn degassing from the natural analog of Campo de Calatrava(Spain): Implications for monitoring of CO2storage sites

Natural analogs offer a valuable opportunity to investigate the long-term impacts associated with thepotential leakage in geological storage of CO2.Degassing of CO2and radon isotopes (222Rn?220Rn) from soil, gas vents and thermal water dischargeswas investigated in the natural analog of Campo de Calatrava Volcanic Field (CCVF; Central Spain) todetermine the CO2?Rn relationships and to assess the role of CO2as carrier gas for radon. Furthermore,radon measurements to discriminate between shallow and deep gas sources were evaluated under theperspective of their applicability in monitoring programs of carbon storage projects.CO2flux as high as 5000 g m?2d?1and222Rn activities up to 430 kBq m?3were measured;220Rn activi-ties were one order of magnitude lower than those of222Rn. The222Rn/220Rn ratios were used to constrainthe source of the Campo de Calatrava soil gases since a positive correlation between radon isotopic ratiosand CO2fluxes was observed. Thus, in agreement with previous studies, our results indicate a deepmantle-related origin of CO2for both free and soil gases, suggesting that carbon dioxide is an efficientcarrier for Rn. Furthermore, it was ascertained that the increase of222Rn in the soil gases was likely pro-duced by two main processes: (i) direct transport by a carrier gas, i.e., CO2and (ii) generation at shallowlevel due to the presence of relatively high concentrations of dissolved U and Ra in the thermal aquiferof Campo de Calatrava.The diffuse CO2soil flux and radon isotopic surveys carried out in the Campo de Calatrava VolcanicFields can also be applicable to geochemical monitoring programs in CCS (Carbon Capture and Storage)areas as these parameters are useful to: (i) constrain CO2leakages once detected and (ii) monitor both theevolution of the leakages and the effectiveness of subsequent remediation activities. These measurementscan also conveniently be used to detect diffuse leakages.

Sea Bottom Classification by Means of Bathymetric LIDAR Data

Light Detection and Ranging (LIDAR) provides high horizontal and vertical resolution of spatial data located in point cloud images, and is increasingly being used in a number of applications and disciplines, which have concentrated on the exploit and manipulation of the data using mainly its three dimensional nature. Bathymetric LIDAR systems and data are mainly focused to map depths in shallow and clear waters with a high degree of accuracy. Additionally, the backscattering produced by the different materials distributed over the bottom surface causes that the returned intensity signal contains important information about the reflection properties of these materials. Processing conveniently these values using a Simplified Radiative Transfer Model, allows the identification of different sea bottom types. This paper presents an original method for the classification of sea bottom by means of information processing extracted from the images generated through LIDAR data. The results are validated using a vector database containing benthic information derived by marine surveys.

Método para el análisis independiente de problemas

habilidades de comprensión y resolución de problemas. Tanto es así que se puede afirmar con rotundidad que no existe el método perfecto para cada una de las etapas de desarrollo y tampoco existe el modelo de ciclo de vida perfecto: cada nuevo problema que se plantea es diferente a los anteriores en algún aspecto y esto hace que técnicas que funcionaron en proyectos anteriores fracasen en los proyectos nuevos. Por ello actualmente se realiza un planteamiento integrador que pretende utilizar en cada caso las técnicas, métodos y herramientas más acordes con las características del problema planteado al ingeniero. Bajo este punto de vista se plantean nuevos problemas. En primer lugar está la selección de enfoques de desarrollo. Si no existe el mejor enfoque, ¿cómo se hace para elegir el más adecuado de entre el conjunto de los existentes? Un segundo problema estriba en la relación entre las etapas de análisis y diseño. En este sentido existen dos grandes riesgos. Por un lado, se puede hacer un análisis del problema demasiado superficial, con lo que se produce una excesiva distancia entre el análisis y el diseño que muchas veces imposibilita el paso de uno a otro. Por otro lado, se puede optar por un análisis en términos del diseño que provoca que no cumpla su objetivo de centrarse en el problema, sino que se convierte en una primera versión de la solución, lo que se conoce como diseño preliminar. Como consecuencia de lo anterior surge el dilema del análisis, que puede plantearse como sigue: para cada problema planteado hay que elegir las técnicas más adecuadas, lo que requiere que se conozcan las características del problema. Para ello, a su vez, se debe analizar el problema, eligiendo una técnica antes de conocerlo. Si la técnica utiliza términos de diseño entonces se ha precondicionado el paradigma de solución y es posible que no sea el más adecuado para resolver el problema. En último lugar están las barreras pragmáticas que frenan la expansión del uso de métodos con base formal, dificultando su aplicación en la práctica cotidiana. Teniendo en cuenta todos los problemas planteados, se requieren métodos de análisis del problema que cumplan una serie de objetivos, el primero de los cuales es la necesidad de una base formal, con el fin de evitar la ambigüedad y permitir verificar la corrección de los modelos generados. Un segundo objetivo es la independencia de diseño: se deben utilizar términos que no tengan reflejo directo en el diseño, para que permitan centrarse en las características del problema. Además los métodos deben permitir analizar problemas de cualquier tipo: algorítmicos, de soporte a la decisión o basados en el conocimiento, entre otros. En siguiente lugar están los objetivos relacionados con aspectos pragmáticos. Por un lado deben incorporar una notación textual formal pero no matemática, de forma que se facilite su validación y comprensión por personas sin conocimientos matemáticos profundos pero al mismo tiempo sea lo suficientemente rigurosa para facilitar su verificación. Por otro lado, se requiere una notación gráfica complementaria para representar los modelos, de forma que puedan ser comprendidos y validados cómodamente por parte de los clientes y usuarios. Esta tesis doctoral presenta SETCM, un método de análisis que cumple estos objetivos. Para ello se han definido todos los elementos que forman los modelos de análisis usando una terminología independiente de paradigmas de diseño y se han formalizado dichas definiciones usando los elementos fundamentales de la teoría de conjuntos: elementos, conjuntos y relaciones entre conjuntos. Por otro lado se ha definido un lenguaje formal para representar los elementos de los modelos de análisis – evitando en lo posible el uso de notaciones matemáticas – complementado con una notación gráfica que permite representar de forma visual las partes más relevantes de los modelos. El método propuesto ha sido sometido a una intensa fase de experimentación, durante la que fue aplicado a 13 casos de estudio, todos ellos proyectos reales que han concluido en productos transferidos a entidades públicas o privadas. Durante la experimentación se ha evaluado la adecuación de SETCM para el análisis de problemas de distinto tamaño y en sistemas cuyo diseño final usaba paradigmas diferentes e incluso paradigmas mixtos. También se ha evaluado su uso por analistas con distinto nivel de experiencia – noveles, intermedios o expertos – analizando en todos los casos la curva de aprendizaje, con el fin de averiguar si es fácil de aprender su uso, independientemente de si se conoce o no alguna otra técnica de análisis. Por otro lado se ha estudiado la capacidad de ampliación de modelos generados con SETCM, para comprobar si permite abordar proyectos realizados en varias fases, en los que el análisis de una fase consista en ampliar el análisis de la fase anterior. En resumidas cuentas, se ha tratado de evaluar la capacidad de integración de SETCM en una organización como la técnica de análisis preferida para el desarrollo de software. Los resultados obtenidos tras esta experimentación han sido muy positivos, habiéndose alcanzado un alto grado de cumplimiento de todos los objetivos planteados al definir el método.---ABSTRACT---Software development is an inherently complex activity, which requires specific abilities of problem comprehension and solving. It is so difficult that it can even be said that there is no perfect method for each of the development stages and that there is no perfect life cycle model: each new problem is different to the precedent ones in some respect and the techniques that worked in other problems can fail in the new ones. Given that situation, the current trend is to integrate different methods, tools and techniques, using the best suited for each situation. This trend, however, raises some new problems. The first one is the selection of development approaches. If there is no a manifestly single best approach, how does one go about choosing an approach from the array of available options? The second problem has to do with the relationship between the analysis and design phases. This relation can lead to two major risks. On one hand, the analysis could be too shallow and far away from the design, making it very difficult to perform the transition between them. On the other hand, the analysis could be expressed using design terminology, thus becoming more a kind of preliminary design than a model of the problem to be solved. In third place there is the analysis dilemma, which can be expressed as follows. The developer has to choose the most adequate techniques for each problem, and to make this decision it is necessary to know the most relevant properties of the problem. This implies that the developer has to analyse the problem, choosing an analysis method before really knowing the problem. If the chosen technique uses design terminology then the solution paradigm has been preconditioned and it is possible that, once the problem is well known, that paradigm wouldn’t be the chosen one. The last problem consists of some pragmatic barriers that limit the applicability of formal based methods, making it difficult to use them in current practice. In order to solve these problems there is a need for analysis methods that fulfil several goals. The first one is the need of a formal base, which prevents ambiguity and allows the verification of the analysis models. The second goal is design-independence: the analysis should use a terminology different from the design, to facilitate a real comprehension of the problem under study. In third place the analysis method should allow the developer to study different kinds of problems: algorithmic, decision-support, knowledge based, etc. Next there are two goals related to pragmatic aspects. Firstly, the methods should have a non mathematical but formal textual notation. This notation will allow people without deep mathematical knowledge to understand and validate the resulting models, without losing the needed rigour for verification. Secondly, the methods should have a complementary graphical notation to make more natural the understanding and validation of the relevant parts of the analysis. This Thesis proposes such a method, called SETCM. The elements conforming the analysis models have been defined using a terminology that is independent from design paradigms. Those terms have been then formalised using the main concepts of the set theory: elements, sets and correspondences between sets. In addition, a formal language has been created, which avoids the use of mathematical notations. Finally, a graphical notation has been defined, which can visually represent the most relevant elements of the models. The proposed method has been thoroughly tested during the experimentation phase. It has been used to perform the analysis of 13 actual projects, all of them resulting in transferred products. This experimentation allowed evaluating the adequacy of SETCM for the analysis of problems of varying size, whose final design used different paradigms and even mixed ones. The use of the method by people with different levels of expertise was also evaluated, along with the corresponding learning curve, in order to assess if the method is easy to learn, independently of previous knowledge on other analysis techniques. In addition, the expandability of the analysis models was evaluated, assessing if the technique was adequate for projects organised in incremental steps, in which the analysis of one step grows from the precedent models. The final goal was to assess if SETCM can be used inside an organisation as the preferred analysis method for software development. The obtained results have been very positive, as SETCM has obtained a high degree of fulfilment of the goals stated for the method.

Dynamic analysis of translational shells

In this paper a dynamic analysis of transnational shells is presented. The general linear shell theory is used in conjunction with additional shallow and curved plate approximations. In order to apply some type of extended Levy solution, the shell is assumed to be limited by a rectangular plan form, with two opposite edges simply supported (gable boundary conditions). First, the shells free vibrations are studied in the usual way, obtaining for each Fourier term the natural frequencies as solutions of a transcendental equation. However, solving these equations arises enormous computational difficulties. This paper deals specifically with this problem, trying to reduce its dimension by a discretization procedure. In the shell dynamic characteristics, namely the mass. The shell mass is lumped along a family of coordinate lines. Therefore, the natural frequencies for each harmonic term can be found from the solution of a typical matrix eigenvalues problem and standard numerical techniques can be applied. The shell response to forced vibrations, particularly to earthquake excitation, can be determined by using conventional procedure either in the time or in the frequency domain. Finally, extending the above procedure, any system of translational shells under dynamic loading can be studied. Then, by using matrix methods, a general computer program is written and applied to some illustrative examples. Numerical results has been obtained in two cases: circular cylindrical shell and box girder bridge.

Propiedades ópticas y estructurales de nanoestructuras de Zn(Cd)(Mg)O

En la presente tesis doctoral se ha realizado un estudio utilizando diferentes técnicas de crecimiento (RPE-MOCVD y spray pyrolysis) y estructuras (nanohilos, pozos y puntos cuánticos y capas) con el objetivo de desarrollar dispositivos que cubran desde el rango visible hasta el ultravioleta. Es por esta razón por la que se han elegido materiales basados en ZnO, debido a la posibilidades que estos ofrecen para variar su bandgap en un amplio rango de energías. Prueba de ello es que en este estudio se ha conseguido cubrir un rango espectral desde 1.86 hasta 4.11 eV, estudiandose además fenómenos físicos como son la difusión e incorporaci ón de la aleación o la adsorción de gases en la super_cie, lo que ha permitido la fabricación de diferentes fotodetectores de gran sensibilidad. Por todo ello, los resultados obtenidos en esta tesis suponen una gran contribución al conocimiento de las propiedades físicas de las aleaciones de Zn(Cd)O y Zn(Mg)O para potenciales aplicaciónes en dispositivos que operen en el rango visible y ultravioleta del espectro, respectivamente. En esta memoria se da en primer lugar una visión de las propiedades de materiales basados en ZnO, entrando en detalle en una de las ventajas que este presenta, la facilidad que tiene este material para formar nanoestructuras. En el capítulo 3 se dan los conceptos teóricos necesarios para comprender las propiedades ópticas de este tipo de materiales, mostrando también los resultados más reseñables obtenidos en ZnO. En los capítulos referentes a los resultados se pueden diferenciar dos grandes bloques. En el primer bloque de resultados se han analizado nanohilos y pozos cuánticos de Zn(Cd)O crecidos por la técnica de RPE-MOCVD (Capítulos 4 y 5). En el segundo se expondrá el estudio realizado sobre capas y puntos cuánticos de Zn(Mg)O crecidos por la técnica spray pyrolysis como se describe en mayor detalle a continuación. Nanohilos y pozos cuánticos de Zn(Cd)O crecidos por RPE-MOCVD Teóricamente aleando el ZnO con CdO es posible disminuir el valor del band- gap desde 3.37 eV hasta 0.95 eV, cubriendo por completo el espectro visible. El desarrollo del ternario Zn(Cd)O permitiría la fabricación de heteroestructuras y pozos cuánticos, muy importantes en el desarrollo de dispositivos optoelectrónicos que cubran la parte visible del espectro. Sin embargo, la diferencia de estructura cristalina entre estos dos materiales junto a la baja solubilidad del Cd y su alta presión de vapor, di_culta la obtención de material de alta calidad cristalina con alto contenido en Cd. En esta tesis doctoral se ha realizado una completa caracterización óptica y estructural de nanohilos de Zn(Cd)O credidos por la técnica de RPE-MOCVD. Estos nanohilos tinene unas longitudes comprendidas entre 1 y 3 _m y diámetros entre 100 y 200 nm. La concentración máxima introducida de Cd en estas estructuras ha sido de hasta un 54% manteniendo la estructura wurtzita del ZnO, siendo este el mayor contenido de Cd introducido hasta la fecha en nanostructuras basada en ZnO. Este hecho se traduce en una variación de la energía de emisión entre 3.31 y 1.86 eV con el aumento en Cd. El uso de diferentes técnicas de alta resoluci ón de caracterización estructural ha permitido demostrar la presencia de una sola fase estructural wurtzita sin observarse ningún indicio de separación de fases ni acumulación de Cd a lo largo del nanohilo para todos los contenidos de Cd. Con el propósito de fabricar dispositivos en nanohilos individuales, parte de esta tesis doctoral ha estado dedicada a estudiar el impacto que el recocido térmico tiene en las propiedades ópticas y eléctricas de nanohilos de Zn(Cd)O. El recocido térmico es un proceso clave en la optimización de dispositivos, ya sea para la obtenci ón de contactos óhmicos, reducción de defectos o difusión de dopantes por ejemplo. En este estudio se ha observado una mejora muy signi_cativa de las propiedades de emisión de los nanohilos cuando estos eran recocidos a temperaturas mayores que la de crecimiento (300 oC). En las muestras con Cd se ha observado además que el recocido también produce un desplazamiento de la emisión hacia mayores energías debido a una reducción homogénea del contenido de Cd. Medidas de fotoluminiscencia con resolución temporal muestran el impacto que tiene la localización del excitón en las _uctuaciones de potencial, debidas a una distribución estadística del Cd, en la dinámica de los portadores. Comparando el tiempo de vida de los portadores entre los nanohilos recocidos y sin recocer se ha observado un aumento de este parámetro en las estructuras recocidas. Este aumento es fundamentalmente debido a una reducción de centros de recombinación no radiativa asociados a defectos presentes a lo largo del nanohilo. Además, se ha estudiado la evolución de los tiempos de vida de los portadores en función de la temperatura, registrándose una menor estabilidad con la temperatura de los tiempos de vida en las muestras recocidas. Este resultado sugiere que el recocido térmico consigue reducir parte del desorden de la aleación en la estructura. Tras haber caracterizados los nanohilos se desarrollaron una serie de procesa dos para la fabricación de dispositivos basados en nanohilos individuales. Se fabricaron en concreto fotodetectores sensibles al UV, en los que se observó también la alta sensibilidad que muestran a la adsorción de gases en la super_cie, incrementada por la gran relación super_cie/volúmen característica de las nanoestructuras. Estos procesos de adsorción observados tienen un impacto directo sobre las propiedades ópticas y electricas de los dispositivos como se ha demostrado. Por ello que en esta tesis se hayan estudiado en detalle este tipo de procesos, ideando maneras para tener un mayor control sobre ellos. Finalmente se crecieron estructuras de pozos cuántico de ZnCdO/ZnO en nanohilos con contenidos de Cd nominales de 54 %. Las medidas ópticas realizadas mostraron como al aumentar la anchura del pozo de 0.7 a 10 nm, la emisión relacionada con el pozo se desplazaba entre 3.30 y 1.97 eV. Este gran desplazamiento representa el mayor obtenido hasta la fecha en pozos cuánticos de ZnCdO/ZnO. Sin embargo, al caracterizar estructuralmente estas muestras se observó la presencia de procesos de difusión de Cd entre el pozo y la barrera. Como se ha podido medir, este tipo de procesos reducen sustancialmente la concentración de Cd en el pozo al difundirse parte a la barrera. cambiando completamente la estructura de bandas nominal de estas estructuras. Este estudio demuestra la importancia del impacto de los procesos de difusión en la interpretación de los efectos de con_namiento cuántico para este tipo de estructuras. Capas y puntos cuánticos de Zn(Mg)O crecidos por spray pyrolysis La técnica de spray pyrolysis, debido a su simplicidad, bajo coste y capacidad de crecer sobre grandes áreas conservando una alta calidad cristalina presenta un gran interés en la comunidad cientí_ca para el potencial desarrollo de dispositivos comerciales. En esta tesis se ha estudiado las propiedades ópticas y eléctricas de capas y puntos cuánticos de Zn(Mg)O crecidos por esta técnica. Al contrario que pasa con el Cd, al introducir Mg en la estructura wurtzita de ZnO se consigue aumentar el bandgap del semiconductor. Sin embargo, al igual que pasa con el CdO, la diferencia de estructura cristalina entre el ZnO y el MgO limita la cantidad de Mg que se puede incorporar, haciendo que para una cierta concentración de Mg aparezcan el fenómeno de separación de fases. En esta tesis se ha conseguido incorporar hasta un contenido de Mg del 35% en la estructura wurtzita del ZnO utilizando la técnica de spray pyrolysis, resultado que representa la mayor concentración de Mg publicada hasta la fecha. Este hecho ha posibilitado variar la energía del borde de absorción desde 3.30 a 4.11 eV. En estas capas se ha realizado una completa caracterización óptica observándose una diferencia entre las energías del borde de absorción y del máximo de emisión creciente con el contenido en Mg. Esta diferencia, conocida como desplazamiento de Stokes, es debida en parte a la presencia de _uctuaciones de potencial producidas por un desorden estadístico de la aleación. Se han fabricado fotodetectores MSM de alta calidad utilizando las capas de Zn(Mg)O previamente caracterizadas, observándose un desplazamiento del borde de absorción con el aumento en Mg desde 3.32 a 4.02 eV. Estos dispositivos muestran altos valores de responsividad (10-103 A/W) y altos contrastes entre la responsividad bajo iluminación y oscuridad (10-107). Estos resultados son en parte debidos a la presencia de mecanismos de ganancia y una reducción de la corriente de oscuridad en las muestras con alto contenido de Mg. Utilizando esta misma técnica de crecimiento se han crecido puntos de Zn(Mg)O con concentraciones nominales de Mg entre 0 y 100 %, con dimensiones medias entre 4 y 6 nm. Las medidas estructurales realizadas muestran que hasta un valor de Mg de 45 %, los puntos están compuestos por una única fase estructural, wurtzita. A partir de esa concentración de Mg aparece una fase cúbica en los puntos, coexistiendo con la fase hexagonal hasta una concentración nominales del 85 %. Para concentraciones mayores de Mg, los puntos muestran una única fase estructural cúbica. Medidas de absorción realizadas en estos puntos de Zn(Mg)O muestran un desplazamiento del borde de absorción entre 3.33 y 3.55 eV cuando la concentraci ón de Mg en los puntos aumenta hasta el 40 %. Este desplazamiento observado es debido solamente a la fase wurtzita del Zn(Mg)O donde se incorpora el Mg. ABSTRACT This PhD theis presents a study using di_erent growth techniques (RPEMOCVD and spray pyrolysis) and structures (nanowires, quantum dots and wells and layers) in order to develop devices that extend from the visible to the ultraviolet range. For this reason ZnO based materials have been choosen, because they o_er the possibility to tunne the bandgap in this energy range. Proof of this is that this study has managed to cover a spectral range from 1.86 to 4.11 eV, also being studied physical phenomena such as di_usion and incorporation of alloy or adsorption of gases on the surface, allowing the develop di_erent highly sensitive photodetectors. Therefore, the results obtained in this thesis are a great contribution two large blockso the knowledge of the physical properties of alloys Zn(Cd)O and Zn(Mg)O for potential applications in devices that operate in the visible and ultraviolet range, respectively. In the _rst chapter, the general properties of ZnO-based materials are presented, showing the facilities that these kind of materials o_er to obtain di_erent nanoestructures. In Chapter 3, optical theoretical concepts are given to understand the optical properties of these materials, also showing the most signi_cant results of ZnO. In the chapters related with the results, two blocks could be distinguish. In the _rst one, Zn(Cd)O nanowires and quantum wells grown by RPE-MOCVD have been analyzed (Chapters 4 and 5). The second block of results shows the study performed in Zn(Mg)O _lms and quantum dots grown by spray pyrolysis. Zn(Cd)O nanowires and quantum wells grown by RPE-MOCVD In summary, the results of the PhD thesis are a great contribution to the knowledge of the physical properties of Zn(Cd)O and Zn(Mg)O alloys and their application for high performance devices operating in the visible and UV ranges, respectively. The performance of the device is still limited due to alloy solubility and p-doping stability, which opens a door for future research in this _eld. Theoretically, annealing ZnO with CdO allows to reduce the bandgap from 3.37 to 0.95 eV, covering the whole visible spectrum. The development of ZnCdO alloys allows the fabrication of heterostructures and quantum wells, necessary for the development of high performance optoelectronic devices. However, the di_erent crystal structures between CdO and ZnO and the low solubility of Cd and its high vapor pressure, hinders the growth of ZnCdO alloys with high Cd contents. In this PhD thesis Zn(Cd)O nanowires have been optically and structurally characterized, obtaining a maximum Cd content of 54% while maintaining their wurtzite structure. This Cd content, which allows lowering the bandgap down to 1.86 eV, is the highest concentration ever reported in nanostructures based on ZnO. The combination of optical and structural characterization techniques used during this thesis has allowed the demonstration of the presence of a single wurtzite structure, without observing any indication of phase separation or Cd accumulation along the nanowire. Annealing processes are essential in the fabrication of optoelectronic devices. For this reason, a complete study of the annealing e_ects in the optical and electrical properties of Zn(Cd)O nanowires has been performed. In the _rst place, annealing nanowires at higher temperatures than their growth temperature (300 oC) allows a signi_cant improvement of their emission properties. However, in the samples that contain Cd a shift in the emission towards higher energies has been observed due to a homogeneous reduction of the Cd content in the nanowires. Time resolved photoluminescence measurements show the impact of the exciton localization in the potential _uctuations due to a statistical alloy disorder. An increase in the carrier lifetime has been obtained for the annealed nanowires. This increase is mainly due to the reduction of non-radiative recombination centers associated with the defects present in the material. Furthermore, temperature dependent time resolved photoluminescence measurements suggest a reduction of the alloy disorder in the annealed samples. In this thesis, single nanowire photodetectors with a high responsivity in the UV range have been demonstrated. Due to the high surface/volume ratio, these structures are very sensitive to gas adsorption at the surface, which largely de_nes the optical and electrical properties of the material and, therefore, of the device. With the aim of obtaining time stable devices, the dynamic adsorption-desorption processes have been studied, developing di_erent approaches that allow a higher control over them. Finally, ZnCdO/ZnO quantum wells have been grown with a nominal Cd concentration of 54% inside the well. The performed optical measurements show that increasing the well width from 0.7 to 10 nm, shifts the emission related with the well from 3.30 to 1.97 eV. This result represents the highest shift reported in the literature. However, a detailed structural characterization shows the presence of di_usion phenomena which substantially reduce the concentration of Cd in the well, while increasing it in the barrier. This type of phenomena should be considered when ac curately interpretating the quantum con_nement e_ects in Zn(Cd)O/ZnO quantum wells. Theoretically, annealing ZnO with CdO allows to decrease the bandgap from 3.37 to 0.95 eV, covering the whole visible spectrum. Zn(Mg)O _lms and quantum dots grown by spray pyrolysis Due to its simplicity, low-cost and capacity to grow over large areas conserving a high crystal quality, spray pyrolysis technique presents a great interest in the scienti_c community for developing comercial devices. In this thesis, a complete study of the optical and structural properties of Zn(Mg)O _lms and quantum dots grown by spray pyrolysis has been performed. Contrary to Zn(Cd)O alloys, when introducing Mg in the ZnO wurtzite structure an increase in the bandgap in obtained. Once again, the di_erence in the crystal structure of ZnO and MgO limits the amount of Mg that can be introduced before phase separation appears. In this PhD thesis, a maximum Mg content of 35% has been incorporated in the wurtzite structure using spray pyrolysis. This variation in the Mg content translates into an increase of the absorption edge from 3.30 to 4.11 eV. Up to this date, this result represents the highest Mg content introduced by spray pyrolysis in a ZnO wurzite structure reported in the literature. The comparison of the emission and absorption spectra shows the presence of an increasing Stokes shift with Mg content. This phenomenon is partialy related with the presence of potential _uctuations due to an statistic alloy disorder. MSM photodetectors have been processed on previously characterized Zn(Mg)O _lms. These devices have shown a shift in the absorption edge from 3.32 to 4.02 eV with the increase in Mg content, high responsivity values (10-103 A/W) and high contrast ratios between illuminated and dark responsivities (10-107). These values are explained by the presence of a gain mechanism and a reduction of dark current in the ZnMgO samples. Zn(Mg)O quantum dots have also been grown using spray pyrolysis with Mg concentrations between 0 and 100% and with average widths ranging 4 to 6 nm. Structural measurements show that at a Mg concentration of 45% the cubic phase appears, coexisting with the hexagonal phase up to an 85% concentration of Mg content. From 85% onwards the quantum dots show only the cubic phase. Absorption measurements performed in these structures reveal a shift in the absorption edge from 3.33 to 3.55 eV when the Mg content increases up to 40 %.

Investigación de métodos de recuperación secundaria en yacimientos kársticos complejos en presencia de acuífero activo infinito

El yacimiento Casablanca es un campo petrolífero maduro en etapa de agotamiento considerado como el más grande del mar Mediterráneo. Lleva en explotación desde 1977 y tiene una producción acumulada de 22.6 MMm3 de petróleo. La formación productiva consiste en carbonatos karstificados del Grupo Basal Terciario y del Mesozoico. El mecanismo de drenaje identificado es por empuje de agua de un gran acuífero activo considerado como infinito ya que ha mantenido la presión al 95% de la original después de casi 40 años de producción. En el año 1999, los pozos asociados al campo Casablanca producían unos 500 m3/d de agua que era tratada y vertida al mar. Para cumplir con las leyes medio ambientales de la época, se convirtió el pozo Casablanca-9 en pozo sumidero con el objetivo de devolver a la formación toda el agua producida de una manera segura, limpia y totalmente respetuosa con el medio ambiente. Años después se observó que ésta inyección no era inocua, sino que tenía un impacto en la producción de petróleo. En la presente tesis se ha definido una metodología que, mediante la experimentación en campo con trazadores, pruebe la existencia de comunicación entre pozos productores y pozos sumidero, rompiendo así el paradigma instaurado en el campo que reza que no es posible la recuperación mejorada mediante inyección de agua en Casablanca debido al gran acuífero existente. Los resultados obtenidos serán el punto de partida para la construcción de un modelo de simulación que permita verificar que es posible la aplicación de técnicas IOR/EOR, y más concretamente la recuperación mejorada mediante inyección de agua en presencia de un acuífero activo infinito. ABSTRACT Casablanca is a brown field in the decline stage and is considered as the largest field in the Mediterranean Sea. It has been on production since 1977 and the cumulative production is 22.6 MMm3 of oil. The productive reservoir formation consists on complex karstified carbonates from Basal Tertiary Group and Mesozoic. The drive mechanism identified is water drive by a large aquifer considered as infinite acting due to the pressure maintenance at 95% of the original after near 40 years of production. In 1999, the wells associated to Casablanca field produced about 500 m3/d of water that was treated and disposed to the sea. In order to comply with the environmental laws at that time, Casablanca-9 was converted from producer to water disposal well with the objective to dispose all the water back to the formation in a safe, clean and environmental fully respectful way. Years later, it was observed that injection was not innocuous, but had an impact on oil production. The methodology defined in this thesis will demonstrate the existence of communication between producers and disposal wells through field experiments with tracers, breaking the paradigm established in the field that says it is not possible to apply waterflooding techniques in Casablanca due to the existence of a strong infinite acting aquifer. The results obtained will be the starting point in order to build a simulation model able to demonstrate that the application of IOR / EOR techniques are suitable, more specifically water flooding techniques in presence of an infinite active aquifer.

Development of a novel autonomous robot for navigation and inspect in oil wells

This paper proposes a novel robotic system that is able to move along the outside of the oil pipelines used in Electric Submersible Pumps (ESP) and Progressive Cavity Pumps (PCP) applications. This novel design, called RETOV, proposes a light weight structure robot that can be equipped with sensors to measure environmental variables avoiding damage in pumps and wells. In this paper, the main considerations and methodology of design and implementation are discussed. Finally, the first experimental results that show RETOV moving in vertical pipelines are analyzed.