Water and nitrogen footprint in an irrigated crop under mineral and organic fertilization

In order to establish rational nitrogen (N) application and reduce groundwater contamination, a clearer understanding of the N distribution through the growing season and its balance is crucial. Excessive doses of N and/or water applied to fertigated crops involve a substantial risk of aquifer contamination by nitrate; but knowledge of N cycling and availability within the soil could assist in avoiding this excess. In central Spain, the main horticultural fertigated crop is the melon type ?piel de sapo¿ and it is cultivated in vulnerable zones to nitrate pollution (Directive 91/676/CEE). However, until few years ago there were not antecedents related to the optimization of nitrogen fertilization together with irrigation. Water and N footprint are indicators that allow assessing the impact generated by different agricultural practices, so they can be used to improve the management strategies in fertigated crop systems. The water footprint distinguishes between blue water (sources of water applied to the crop, like irrigation and precipitation), green water (water used by the crop and stored in the soil), and it is furthermore possible to quantify the impact of pollution by calculating the grey water, which is defined as the volume of polluted water created from the growing and production of crops. On the other hand, the N footprint considers green N (nitrogen consumed by the crops and stored in the soil), blue N (N available for crop, like N applied with mineral and/or organic fertilizers, N applied with irrigation water and N mineralized during the crop period), whereas grey N is the amount of N-NO3- washed from the soil to the aquifer. All these components are expressed as the ratio between the components of water or N footprint and the yield (m3 t-1 or kg N t-1 respectively). The objetives of this work were to evaluate the impact derivated from the use of different fertilizer practices in a melon crop using water and N footprint.

Checking existing formulas for estimating overtopping rates in harbour breakwaters

The purpose of the research work resulting from various studies undertaken in the CEDEX, as summarized in this article, is to make a comparative analysis of methods for calculating overtopping rates developed by different authors. To this effect, in the first place, existing formulas for estimating the overtopping rate on rubble mound and vertical breakwaters were summarised and analysed. Later, the above mentioned formulas were compared using the results obtained in a series of hydraulic model tests at the CEDEX. The results obtained in the Ferrol outer harbour breakwater and Melilla harbour breakwater tests are presented here. A calculation method based on the neural network theory, developed in the European CLASH Project, was applied to a series of sloping breakwater tests in order to complete this research and the results obtained in the Ferrol outer harbour breakwater test are presented in this article. A series of additional tests was also carried out in a physical model on the standard cross section of the Bilbao harbour sloping breakwater’s cross section, the results of which are under study using the empirical formulas applicable to the cross section, as well as the NN-OVERTOPPING neural network

Comparison of Existing Formulae for Estimating the Overtopping Rate in Harbour Breakwaters

The purpose of the research work resulting from various studies undertaken in the CEDEX, as summarized in this article, is to make a comparative analysis of methods for calculating overtopping rates developed by different authors. To this effect, in the first place, existing formulae for estimating the overtopping rate on rubble mound and vertical breakwaters were summarised and analysed. Later, the above mentioned formulae were compared using the results obtained in a series of hydraulic model tests at the CEDEX (the Center of Studies of Ports and Coasts of the CEDEX, Madrid, Spain). A calculation method based on the neural network theory, developed in the European CLASH Project, was applied to a series of sloping breakwater tests in order to complete this research. The results obtained in the Ferrol, Ciervana and Alicante breakwaters tests are presented here.

Overtopping of harbour breakwaters: a comparison of semi-empirical equations, neural networks, and physical model tests

This paper reports extensive tests of empirical equations developed by different authors for harbour breakwater overtopping. First, the existing equations are compiled and evaluated as tools for estimating the overtopping rates on sloping and vertical breakwaters. These equations are then tested using the data obtained in a number of laboratory studies performed in the Centre for Harbours and Coastal Studies of the CEDEX, Spain. It was found that the recommended application ranges of the empirical equations typically deviate from those revealed in the experimental tests. In addition, a neural network model developed within the European CLASH Project is tested. The wind effects on overtopping are also assessed using a reduced scale physical model

Caracterización de parámetros petrofísicos relacionados con la capacidad de almacenamiento de CO2 en el acuífero salino profundo de Hontomín (Burgos)

La conciencia de la crisis de la modernidad -que comienza ya a finales del siglo XIX- ha cobrado más experiencia debido al conocimiento de los límites del desarrollo económico, ya que como parecía razonable pensar, también los recursos naturales son finitos. En 1972, el Club de Roma analizó las distintas opciones disponibles para conseguir armonizar el desarrollo sostenible y las limitaciones medioambientales. Fue en 1987 cuando la Comisión Mundial para el Medio Ambiente y el Desarrollo de la ONU definía por primera vez el concepto de desarrollo sostenible. Definición que posteriormente fue incorporada en todos los programas de la ONU y sirvió de eje, por ejemplo, a la Cumbre de la Tierra celebrada en Río de Janeiro en 1992. Parece evidente que satisfacer la demanda energética, fundamentalmente desde la Revolución Industrial en el s XIX, trajo consigo un creciente uso de los combustibles fósiles, con la consiguiente emisión de los gases de efecto invernadero (GEI) y el aumento de la temperatura global media terrestre. Esta temperatura se incrementó en los últimos cien años en una media de 0.74ºC. La mayor parte del incremento observado desde la mitad del siglo XX en esta temperatura media se debe, con una probabilidad de al menos el 90%, al aumento observado en los GEI antropogénicos, siendo uno de ellos el CO2 que proviene de la transformación del carbono de los combustibles fósiles durante su combustión. Ante el creciente uso de los combustibles fósiles, los proyectos CAC, proyectos de captura, transporte y almacenamiento, se presentan como una contribución al desarrollo sostenible ya que se trata de una tecnología que permite mitigar el cambio climático. Para valorar si la tecnología CAC es sostenible, habrá que comprobar si existe o no capacidad para almacenar el CO2 en una cantidad mayor a la de producción y durante el tiempo necesario que impone la evolución de la concentración de CO2 en la atmósfera para mantenerla por debajo de las 450ppmv (concentración de CO2 que propone el Panel Intergubernamental para el Cambio Climático). El desarrollo de los proyectos CAC completos pasa por la necesaria selección de adecuados almacenes de CO2 que sean capaces de soportar los efectos de las presiones de inyección, así como asegurar la capacidad de dichos almacenes y la estanqueidad del CO2 en los mismos. La caracterización geológica de un acuífero susceptible de ser almacén de CO2 debe conducir a determinar las propiedades que dicho almacén posee para asegurar un volumen adecuado de almacenamiento, una inyectabilidad del CO2 en el mismo a un ritmo adecuado y la estanqueidad del CO2 en dicho acuífero a largo plazo. El presente trabajo pretende estudiar los parámetros que tienen influencia en el cálculo de la capacidad del almacén, para lo que en primer lugar se ha desarrollado la tecnología necesaria para llevar a cabo la investigación mediante ensayos de laboratorio. Así, se ha desarrollado una patente, "ATAP, equipo para ensayos petrofísicos (P201231913)", con la que se ha llevado a cabo la parte experimental de este trabajo para la caracterización de los parámetros que tienen influencia en el cálculo de la capacidad del almacén. Una vez desarrollada la tecnología, se aborda el estudio de los distintos parámetros que tienen influencia en la capacidad del almacén realizando ensayos con ATAP. Estos ensayos definen el volumen del almacenamiento, llegándose a la conclusión de que en la determinación de este volumen, juegan un papel importante el alcance de los mecanismos trampa, físicos o químicos, del CO2 en el almacén. Ensayos que definen la capacidad del almacén de "aceptar" o "rechazar" el CO2 inyectado, la inyectabilidad, y por último, ensayos encaminados a determinar posibles fugas que se pueden dar a través de los pozos de inyección, definidos estos como caminos preferenciales de fugas en un almacén subterráneo de CO2. Queda de este modo caracterizada la estanqueidad del CO2 en el acuífero a largo plazo y su influencia obvia en la determinación de la capacidad del almacén. Unido al propósito de la estimación de la capacidad del almacén, se encuentra el propósito de asegurar la estanqueidad de dichos almacenes en el tiempo, y adelantarse a la evolución de la pluma de CO2 en el interior de dichos almacenes. Para cumplir este propósito, se ha desarrollado un modelo dinámico a escala de laboratorio, mediante el programa ECLIPSE 300, con el fin de establecer una metodología para el cálculo de la capacidad estimada del almacén, así como el estudio de la evolución de la pluma de CO2 dentro del acuífero a lo largo del tiempo, partiendo de los resultados obtenidos en los ensayos realizados en ATAP y con la modelización de la probeta de roca almacén empleada en dichos ensayos. Presentamos por tanto un trabajo que establece las bases metodológicas para el estudio de la influencia de distintos parámetros petrofísicos en el cálculo de la capacidad del almacén unidos al desarrollo tecnológico de ATAP y su utilización para la determinación de dichos parámetros aplicables a cada acuífero concreto de estudio. ABSTRACT The crisis of modernity –which begins at the end of 19th Century- has been more important due to the knowledge of the limits of economic development, since it appeared to be thought reasonable, the natural resources are finite. In 1972, The Club of Rome analyzed the different options available in order to harmonize the sustainability and the environment development. It was in 1987 when The Global Commission on The Environment and the Development of UN, defined for the first time the concept of Sustainable Development. This definition that was fully incorporated in all the UN programs and it was useful as an axis, for example, in La Cumbre de la Tierra summit in Río de Janeiro in 1992. It seems obvious to satisfy energetic demand, basically after The Industrial Revolution in 19th Century, which represented an increasing use of fossil fuels, therefore greenhouse gases emission and the increasing of global average temperature. This temperature increased in the last 100 years up to 0.74ºC. The major part of the temperature increase is due to the increase observed in Greenhouse gases with human origin, at least with 90% of probability. The most important gas is the CO2 because of its quantity. In the face of the increasing use of fossil fuels, the CCS projects, Carbon Capture and Storage projects, appear as a contribution of sustainable development since it is a technology for avoiding the climate change. In order to evaluate if CCS technology is sustainable, it will be necessary to prove if the capacity for CO2 storage is available or not in a quantity greater than the production one and during the time necessary to keep the CO2 concentration in the atmosphere lower than 450ppmv (concentration imposed by IPCC). The development of full CCS projects goes through the selection of good CO2 storages that are able to support the effects of pressure injection, and assure the capacity of such storages and the watertightness of CO2. The geological characterization of the aquifer that could be potential CO2 storage should lead to determine the properties that such storage has in order to assure the adequate storage volume, the CO2 injectivity in a good rate, and the watertightness of the CO2 in the long term. The present work aims to study the parameters that have influence on the calculation of storage capacity, and for that purpose the appropriate technology has been developed for carrying out the research by mean of laboratory tests. Thus, a patent has been developed, "ATAP, equipo para ensayos petrofísicos (P201231913)", that has been used for developing the experimental part of this work. Once the technology has been developed, the study of different parameters, that have influence on the capacity of the storage, has been addressed developing different tests in ATAP. These tests define the storage volume which is related to the scope of different CO2 trap mechanisms, physical or chemical, in the storage. Tests that define the capacity of the storage to “accept” or “reject” the injected CO2, the injectivity, and tests led to determine possible leakages through injection wells. In this way we could talk about the watertightness in the aquifer in the long term and its influence on the storage capacity estimation. Together with the purpose of the storage capacity estimation, is the purpose of assuring the watertightness of such storages in the long term and anticipating the evolution of CO2 plume inside such aquifers. In order to fulfill this purpose, a dynamic model has been developed with ECLIPSE 300, for stablishing the methodology for the calculation of storage capacity estimation and the evolution of the CO2 plume, starting out with the tests carried out in ATAP. We present this work that establishes the methodology bases for the study of the influence of different petrophysics parameters in the calculation of the capacity of the storage together with the technological development of ATAP and its utilization for the determination of such parameters applicable to each aquifer.

Análisis comparativo y nuevas estimaciones de las tasas de rebase en obras marítimas en grandes profundidades

El rebase se define como el transporte de una cantidad importante de agua sobre la coronación de una estructura. Por tanto, es el fenómeno que, en general, determina la cota de coronación del dique dependiendo de la cantidad aceptable del mismo, a la vista de condicionantes funcionales y estructurales del dique. En general, la cantidad de rebase que puede tolerar un dique de abrigo desde el punto de vista de su integridad estructural es muy superior a la cantidad permisible desde el punto de vista de su funcionalidad. Por otro lado, el diseño de un dique con una probabilidad de rebase demasiado baja o nula conduciría a diseños incompatibles con consideraciones de otro tipo, como son las estéticas o las económicas. Existen distintas formas de estudiar el rebase producido por el oleaje sobre los espaldones de las obras marítimas. Las más habituales son los ensayos en modelo físico y las formulaciones empíricas o semi-empíricas. Las menos habituales son la instrumentación en prototipo, las redes neuronales y los modelos numéricos. Los ensayos en modelo físico son la herramienta más precisa y fiable para el estudio específico de cada caso, debido a la complejidad del proceso de rebase, con multitud de fenómenos físicos y parámetros involucrados. Los modelos físicos permiten conocer el comportamiento hidráulico y estructural del dique, identificando posibles fallos en el proyecto antes de su ejecución, evaluando diversas alternativas y todo esto con el consiguiente ahorro en costes de construcción mediante la aportación de mejoras al diseño inicial de la estructura. Sin embargo, presentan algunos inconvenientes derivados de los márgenes de error asociados a los ”efectos de escala y de modelo”. Las formulaciones empíricas o semi-empíricas presentan el inconveniente de que su uso está limitado por la aplicabilidad de las fórmulas, ya que éstas sólo son válidas para una casuística de condiciones ambientales y tipologías estructurales limitadas al rango de lo reproducido en los ensayos. El objetivo de la presente Tesis Doctoral es el contrate de las formulaciones desarrolladas por diferentes autores en materia de rebase en distintas tipologías de diques de abrigo. Para ello, se ha realizado en primer lugar la recopilación y el análisis de las formulaciones existentes para estimar la tasa de rebase sobre diques en talud y verticales. Posteriormente, se llevó a cabo el contraste de dichas formulaciones con los resultados obtenidos en una serie de ensayos realizados en el Centro de Estudios de Puertos y Costas. Para finalizar, se aplicó a los ensayos de diques en talud seleccionados la herramienta neuronal NN-OVERTOPPING2, desarrollada en el proyecto europeo de rebases CLASH (“Crest Level Assessment of Coastal Structures by Full Scale Monitoring, Neural Network Prediction and Hazard Analysis on Permissible Wave Overtopping”), contrastando de este modo la tasa de rebase obtenida en los ensayos con este otro método basado en la teoría de las redes neuronales. Posteriormente, se analizó la influencia del viento en el rebase. Para ello se han realizado una serie de ensayos en modelo físico a escala reducida, generando oleaje con y sin viento, sobre la sección vertical del Dique de Levante de Málaga. Finalmente, se presenta el análisis crítico del contraste de cada una de las formulaciones aplicadas a los ensayos seleccionados, que conduce a las conclusiones obtenidas en la presente Tesis Doctoral. Overtopping is defined as the volume of water surpassing the crest of a breakwater and reaching the sheltered area. This phenomenon determines the breakwater’s crest level, depending on the volume of water admissible at the rear because of the sheltered area’s functional and structural conditioning factors. The ways to assess overtopping processes range from those deemed to be most traditional, such as semi-empirical or empirical type equations and physical, reduced scale model tests, to others less usual such as the instrumentation of actual breakwaters (prototypes), artificial neural networks and numerical models. Determining overtopping in reduced scale physical model tests is simple but the values obtained are affected to a greater or lesser degree by the effects of a scale model-prototype such that it can only be considered as an approximation to what actually happens. Nevertheless, physical models are considered to be highly useful for estimating damage that may occur in the area sheltered by the breakwater. Therefore, although physical models present certain problems fundamentally deriving from scale effects, they are still the most accurate, reliable tool for the specific study of each case, especially when large sized models are adopted and wind is generated Empirical expressions obtained from laboratory tests have been developed for calculating the overtopping rate and, therefore, the formulas obtained obviously depend not only on environmental conditions – wave height, wave period and water level – but also on the model’s characteristics and are only applicable in a range of validity of the tests performed in each case. The purpose of this Thesis is to make a comparative analysis of methods for calculating overtopping rates developed by different authors for harbour breakwater overtopping. First, existing equations were compiled and analysed in order to estimate the overtopping rate on sloping and vertical breakwaters. These equations were then compared with the results obtained in a number of tests performed in the Centre for Port and Coastal Studies of the CEDEX. In addition, a neural network model developed in the European CLASH Project (“Crest Level Assessment of Coastal Structures by Full Scale Monitoring, Neural Network Prediction and Hazard Analysis on Permissible Wave Overtopping“) was also tested. Finally, the wind effects on overtopping are evaluated using tests performed with and without wind in the physical model of the Levante Breakwater (Málaga).

Estudio para un Modelo de Gestión de un Acuífero Costero, mediante Metodologías Participativas y Análisis Geoestadístico en el marco del Desarrollo Local. Manglaralto, Ecuador

El presente estudio se fundamenta en la investigación-acción-participativa (IAP), para buscar alternativas que tiendan al desarrollo local de un territorio. Se centra en la cuenca hidrográfica del rio Manglaralto-Santa Elena-Ecuador, aplicando un sistema metodológico participativo que considera las características peculiares del territorio, que se analizan geoespacialmente reconociendo la influencia de la dinámica de sus cambios y observando los móviles que la propiciaban. A través de mecanismos participativos, se conectan los aspectos técnicos para el conocimiento y el aprovechamiento racional del acuífero costero, con los valores de los habitantes del territorio, para mejorar su abastecimiento de agua y crear nuevas condiciones y oportunidades en el camino del desarrollo local, vislumbrando la sostenibilidad. Cabe indicar que el ente administrativo y propulsor es la Junta de Agua Potable Regional Manglaralto (JAPRM). La hipótesis del estudio considera, que los métodos participativos generan en la comunidad una respuesta basada en su identidad y sus deseos de mejorar, que propiciará una gestión del acuífero costero que conlleve al desarrollo local. Otra hipótesis complementaria estipula que las estrategias del gobierno respecto al turismo propicia un crecimiento en la demanda del agua del acuífero. En Manglaralto-Ecuador, una parroquia de 30.000 habitantes aproximadamente, donde la JAPRM, administra y suministra agua a 23.586 habitantes que cuenta en su organización, llevada por 6 representantes de las comunidades rurales que la conforman, empezaron hace 7 años a buscar una forma de lograr un cambio, de tener agua para el desarrollo de la comunidad. Buscaron ayuda por diferentes medios, políticos, económicas, sociales y encontraron como base fundamental a la cooperación con el Organismo Internacional de Energía Atómica (OIEA) y la Escuela Superior Politécnica del Litoral (ESPOL) para entrelazar aspectos técnicos, ambientales, sociales y culturales. La gestión del acuífero costero, desde la perspectiva del IAP repercute en el desarrollo de Manglaralto. También se realiza un análisis geoespacial-geoestadístico, para vislumbrar aspectos de cambios en el territorio ligados al crecimiento turístico, que afectan a la demanda del recurso agua proveniente del acuífero costero bajo la administración de la JAPRM. La tesis presenta el modelo integral y propio de la comunidad de Manglaralto, que refleja una evolución que alcanzó un apogeo en 2011 y parte del 2012, con 9 pozos de agua que daban servicio los 365 días del año, 24 horas al día ininterrumpidamente. Las condiciones externas (promociones turísticas de la ruta del Spondylus) han repercutido en nuevas problemáticas (crecimiento elevado de la demanda del agua). El acuífero costero se convierte en el emblema y móvil de solución, gracias a la gestión integral y a la interacción IAP que se amolda a la evolución de las condiciones, buscando soluciones para la comunidad y su entorno. El modelo integral del territorio con la participación de sus pobladores, considera el aspecto turístico, como un agente que propicia la mayor demanda del agua. Situación a la que hay que dar respuesta mediante la observación-reflexión en el ciclo del IAP para generar nuevas directrices estratégicas y gestionar el desarrollo local. ABSTRACT The present study is based on the participatory action research (PAR) methodology in order to look for alternatives which tend to the local development of a territory. It focuses on the Manglaralto hydrographic river basin located in Santa Elena-Ecuador through the application of the participatory methodology which considers the peculiar characteristics of the territory. These are geospatially analyzed recognizing the influence of its dynamic of changes and observing the causes that originated them. Through the use of participatory mechanisms, technical aspects are connected for stimulating knowledge and rational use of the coastal aquifer with the values of inhabitants of the territory to improve the water supply and create new conditions of sustainability. It is important to point out that the administrative organism and promoter is the Manglaralto Regional Fresh Water Board (JAPRM). In Manglaralto-Ecuador, a parish of approximately 30,000 inhabitants, the MRFWB manages and supplies water to 23.586 inhabitants. This organization is composed by 6 representatives of rural communities. It started 7 years ago looking for a way to achieve a change, from obtaining water to developing the community. They seeked for help in different fields such as: political, economic and social and they found International Atomic Energy Agency (IAEA) and Escuela Superior Politécnica del Litoral (ESPOL) as a fundamental basis for cooperation to bond technical, environmental, social and cultural aspects. Management of coastal aquifer, from the PAR perspective affects the development of Manglaralto. Also, a geospatial and geostatistical analysis is carried out to distinguish change aspects in territories related to touristy growth which affects the demand of water obtained from the coastal aquifer under the management of the MRFWB. The thesis presents a comprehensive model that belongs to the Manglaralto community and reveals an evolution that reached a peak in 2011 and part of 2012, with 9 water wells that operated the 365 days of the year 24 hours a day without interruption. The external conditions (touristic packages of Spondylus route) have created new problems (higher demand of water). The coastal aquifer is a symbol and solution, thanks to the comprehensive management and PAR interaction which fits the evolution of conditions, looking for solutions for the community and its surroundings. The comprehensive model of territory with the participation of inhabitants considers the touristic aspect as an agent which brings about a higher demand of water. This situation requests a response through the observation-reflection in the PAR cycle to generate new strategic guidelines and promote the local development.

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.

Desarrollo de métodos numérico-interpretativos para la realización de ensayos de bombeo

Los ensayos de bombeo son, sin lugar a dudas, una de las pruebas más fiables y de mayor interés que se hacen en el medio físico. No son pruebas estrictamente puntuales, dado que el bombeo atrae flujo desde distancias lejanas al pozo, la prueba tiene una excelente representatividad espacial. Los métodos de interpretación mediante ensayos de bombeo se empezaron a plantear en la primera mitad del pasado siglo. Con los ensayos de bombeo se puede calcular la transmisividad y coeficiente de almacenamiento de las formaciones acuíferas y suministran información sobre el tipo de acuífero, la calidad constructiva del pozo de extracción, la existencia de barreras impermeable o bordes de recarga próximos, e incluso en algunas circunstancias permiten el cálculo del área de embalse subterráneo. Desde mediados del siglo 20 existe una eficaz y abundante gama de métodos analítico-interpretativos de ensayos de bombeo, tanto en régimen permanente como transitorio. Estos métodos son ampliamente conocidos y están muy experimentados a lo largo de muchos países, sin embargo, hoy día, podrían utilizarse modelos de flujo para la interpretación, logrando la misma fiabilidad e incluso mejores posibilidades de análisis. Muchos ensayos que no pueden interpretarse porque las configuraciones del medio son demasiado complejas y no están disponibles, o no es posible, el desarrollo de métodos analíticos, tienen buena adaptación y en ocasiones muy fácil solución haciendo uso de los métodos numéricos de simulación del flujo. En esta tesis se ha buscado una vía de interpretar ensayos de bombeo haciendo uso de modelos de simulación del flujo. Se utiliza el modelo universal MODFLOW del United States Geological Survey, en el cual se configura una celda de simulación y mallado particularmente adecuados para el problema a tratar, se valida con los métodos analíticos existentes. Con la célula convenientemente validada se simulan otros casos en los que no existen métodos analíticos desarrollados dada la complejidad del medio físico a tratar y se sacan las oportunas conclusiones. Por último se desarrolla un modelo específico y la correspondiente aplicación de uso general para la interpretación numérica de ensayos de bombeo tanto con las configuraciones normales como con configuraciones complejas del medio físico. ABSTRACT Pumping tests are, without doubt, one of the most reliable and most interesting tests done in the physical environment. They are not strictly anecdotal evidence, since pumping flow attracts from far distances to the well, the test has excellent spatial representation. Methods of interpretation by pumping tests began to arise in the first half of last century. With pumping tests, can be calculated transmissivity and storage coefficient of the aquifer formations, and provide information on the type of aquifer, the construction quality of the well, the existence of waterproof barriers or borders next recharge, and even in some circumstances allow calculating the area of underground reservoir. Since the mid-20th century there is effective and abundant range of analytical interpretative pumping tests, both in steady state and transient methods. These methods are very widely known and experienced over many countries, however, nowadays, may flow models used for interpretation, obtaining equally reliable or even better possibilities for analysis. Many trials cannot be interpreted as environmental settings are too complex and are not available, or not possible, the development of analytical methods, have good adaptation and sometimes very easily solved using numerical flow simulation methods. This thesis has sought a way to interpret pumping tests using flow simulation models. MODFLOW universal model of United States Geological Survey, in which a simulation cell and meshing particularly suitable for the problem to be treated, is validated with existing analytical methods used is set. With suitably validated cell other cases where there are no analytical methods developed given the complexity of the physical environment to try and draw appropriate conclusions are simulated. Finally, a specific model and the corresponding application commonly used for numerical interpretation of pumping tests both with normal settings as complex configurations of the physical environment is developed.

The problem of the spreading of a liquid film along a solid surface: A new mathematical formulation

A new mathematical model is proposed for the spreading of a liquid film on a solid surface. The model is based on the standard lubrication approximation for gently sloping films (with the no-slip condition for the fluid at the solid surface) in the major part of the film where it is not too thin. In the remaining and relatively small regions near the contact lines it is assumed that the so-called autonomy principle holds—i.e., given the material components, the external conditions, and the velocity of the contact lines along the surface, the behavior of the fluid is identical for all films. The resulting mathematical model is formulated as a free boundary problem for the classical fourth-order equation for the film thickness. A class of self-similar solutions to this free boundary problem is considered.

Análise de conflitos de usos múltiplos da água auxiliada por simulação computacional: bacia do rio Sapucaí - Mirim/Grande

Na mesma proporção que fontes de água potável são reduzidas, a competição por elas tem aumentado, conseqüência imediata da expansão de áreas urbanas e de novas demandas agro-industriais. Os conflitos gerados por esta competição tendem a comprometer o desenvolvimento sustentável e o progresso de uma forma geral, uma vez que resultam em perda de eficácia para vultuosos investimentos públicos e privados e em prejuízos para os usuários menos estruturados. Considerando os fatos acima, o objeto deste trabalho é o estabelecimento de cenários de alocação das disponibilidades hídricas entre usuários de água da bacia do Sapucaí-Mirim/Grande, visando determinar a convivência harmoniosa das várias finalidades de uso provenientes do mesmo recurso hídrico. Tal meta poderá ser alcançada através da análise dos indicadores de desempenho - confiabilidade, vulnerabilidade e resiliência - do sistema os quais serão obtidos com o auxílio do modelo de simulação IRAS (Interactive River Aquifer Simulation).

Investigação de aquífero fraturado para entendimento de fluxo e transporte de contaminantes clorados: estudo de caso em Valinhos, SP

Casos de contaminação de aquíferos fraturados são bastante complexos, tendo em vista a heterogeneidade das redes de fraturas, e no geral, sua investigação demanda a utilização de técnicas pouco usuais, como por exemplo o imageamento acústico e a perfilagem de velocidade de fluxo de água. Na área de estudo, localizada em Valinhos/SP, o uso inadequado de solventes organoclorados no passado ocasionou a contaminação do aquífero raso em duas áreas, e o aparecimento de concentrações no aquifero profundo levaram a condução do atual trabalho, que teve como principal objetivo a elaboração de um modelo conceitual de fluxo de água e transporte de contaminantes no aquífero cristalino. Previamente à investigação do aquífero fraturado, foi realizada uma análise de trabalhos existentes, incluindo a interpretação de lineamentos, levantamentos geológicos além de perfilagens geofísicas de superfície. Em cada área investigada, foi realizada a perfuração de um poço profundo e aplicadas as técnicas de perfilagens de raios gama, cáliper, flowmeter, imageamento acústico, além da filmagem do poço e realização de ensaios hidráulicos nos dois pontos perfurados. Para caracterização química do aquífero fraturado, foram realizadas coletas de água subterrânea em intervalos selecionados com a utilização de obturadores pneumáticos. As cargas hidráulicas medidas durante a amostragem também auxiliaram no entendimento da direção do fluxo de água. O aquífero cristalino é formado por rochas gnáissicas e se encontra bastante fraturado e intemperizado, principalmente na porção superficial da rocha (até aproximadamente 65,0 m) onde as maiores velocidades de fluxo de água também foram observadas. A rocha sã possui uma menor densidade de fraturas e predominância de minerais mais claros. As fraturas de baixo a médio angulo de mergulho (Grupo 1) são as mais frequentes em ambas as perfurações e possuem direção principal N-S a NE-SW. São observadas, no geral, exercendo grande influência sobre o fluxo de água, principalmente na porção alterada do gnaisse. Fraturas com ângulo elevado de mergulho, classificadas como Grupo 2 (paralelas à foliação) e Grupo 3 (direção NW à W), são também observadas ao longo de toda a perfuração estabelecendo a conexão hidráulica entre as fraturas do Grupo 1. Em menor proporção, são ainda verificadas fraturas com ângulos de mergulho >40 ° pertencente aos Grupos 4 (NE-SW), 5 (E-W), 6 (NW-SE) e 7 (E-W). O fluxo de água subterrânea se mostrou descendente na porção superior da rocha alterada e ascendente na porção mais profunda, possivelmente direcionando a água subterrânea para a região de transição da rocha mais alterada para a rocha sã (entre 61 a 65 m de profundidade). Apesar do fluxo ascendente em profundidade, o bombeamento de poços tubulares existentes no entorno ao longo dos anos, favoreceu a migração dos contaminantes para porções mais profundas. Os contaminantes observados no poço tubular P6 possuem maior semelhança com os contaminantes observados na Área 2, e ambos estão localizados entre lineamentos NW-SE, indicando uma possível influência dos lineamentos no controle sobre o fluxo de água. No entanto, para entendimento do transporte dos contaminantes em área, é necessário um adensamento da rede de monitoramento, levando em consideração a heterogeneidade do meio e as incertezas relacionadas à extrapolação dos dados para áreas não investigadas.

Aplicação de métodos geoestatísticos no estudo das distribuições espaciais de condutividade hidráulica em áreas do Sistema Aquífero Guarani (SAG) e do Sistema Aquífero Bauru (SAB) no Estado de São Paulo

A condutividade hidráulica (K) é um dos parâmetros controladores da magnitude da velocidade da água subterrânea, e consequentemente, é um dos mais importantes parâmetros que afetam o fluxo subterrâneo e o transporte de solutos, sendo de suma importância o conhecimento da distribuição de K. Esse trabalho visa estimar valores de condutividade hidráulica em duas áreas distintas, uma no Sistema Aquífero Guarani (SAG) e outra no Sistema Aquífero Bauru (SAB) por meio de três técnicas geoestatísticas: krigagem ordinária, cokrigagem e simulação condicional por bandas rotativas. Para aumentar a base de dados de valores de K, há um tratamento estatístico dos dados conhecidos. O método de interpolação matemática (krigagem ordinária) e o estocástico (simulação condicional por bandas rotativas) são aplicados para estimar os valores de K diretamente, enquanto que os métodos de krigagem ordinária combinada com regressão linear e cokrigagem permitem incorporar valores de capacidade específica (Q/s) como variável secundária. Adicionalmente, a cada método geoestatístico foi aplicada a técnica de desagrupamento por célula para comparar a sua capacidade de melhorar a performance dos métodos, o que pode ser avaliado por meio da validação cruzada. Os resultados dessas abordagens geoestatísticas indicam que os métodos de simulação condicional por bandas rotativas com a técnica de desagrupamento e de krigagem ordinária combinada com regressão linear sem a técnica de desagrupamento são os mais adequados para as áreas do SAG (rho=0.55) e do SAB (rho=0.44), respectivamente. O tratamento estatístico e a técnica de desagrupamento usados nesse trabalho revelaram-se úteis ferramentas auxiliares para os métodos geoestatísticos.

Reconhecimento de monossílabos em idosos: análise do nível de apresentação da fala

A perda auditiva no idoso acarreta em dificuldade na percepção da fala. O teste comumente utilizado na logoaudiometria é a pesquisa do índice de reconhecimento de fala máximo (IR-Max) em uma única intensidade de apresentação da fala. Entretanto, o procedimento mais adequado seria a realização do teste em diversas intensidades, visto que o índice de acerto depende da intensidade da fala no momento do teste e está relacionado com o grau e configuração da perda auditiva. A imprecisão na obtenção do IR-Max poderá gerar uma hipótese diagnóstica errônea e o insucesso no processo de intervenção na perda auditiva. Objetivo: Verificar a interferência do nível de apresentação da fala, no teste de reconhecimento de fala, em idosos com perda auditiva sensorioneural com diferentes configurações audiométricas. Métodos: Participaram 64 idosos, 120 orelhas (61 do gênero feminino e 59 do gênero masculino), idade entre 60 e 88 anos, divididos em grupos: G1- composto por 23 orelhas com configuração horizontal, G2- 55 orelhas com configuração descendente, G3- 42 orelhas com configuração abrupta. Os critérios de inclusão foram: perda auditiva sensorioneural de grau leve a severo, não usuário de aparelho de amplificação sonora individual (AASI), ou com tempo de uso inferior a dois meses, e ausência de alterações cognitivas. Foram realizados os seguintes procedimentos: pesquisas do limiar de reconhecimento de fala (LRF), do índice de reconhecimento de fala (IRF) em diversas intensidades e do nível de máximo conforto (MCL) e desconforto (UCL) para a fala. Para tal, foram utilizadas listas com 11 monossílabos, para diminuir a duração do teste. A análise estatística foi composta pelo teste Análise de Variância (ANOVA) e teste de Tukey. Resultados: A configuração descendente foi a de maior ocorrência. Indivíduos com configuração horizontal apresentaram índice médio de acerto mais elevado de reconhecimento de fala. Ao considerar o total avaliado, 27,27% dos indivíduos com configuração horizontal revelaram o IR-Max no MCL, assim como 38,18% com configuração descendente e 26,19% com configuração abrupta. O IR-Max foi encontrado no UCL, em 40,90% dos indivíduos com configuração horizontal, 45,45% com configuração descendente e 28,20% com configuração abrupta. Respectivamente, o maior e o menor índice médio de acerto foram encontrados em: G1- 30 e 40 dBNS; G2- 50 e 10 dBNS; G3- 45 e 10 dBNS. Não há uma única intensidade de fala a ser utilizada em todos os tipos de configurações audiométricas, entretanto, os níveis de sensação que identificaram os maiores índices médios de acerto foram: G1- 20 a 30 dBNS, G2- 20 a 50 dBNS; G3- 45 dBNS. O MCL e o UCL-5 dB para a fala não foram eficazes para determinar o IR-Max. Conclusões: O nível de apresentação teve influência no desempenho no reconhecimento de fala para monossílabos em idosos com perda auditiva sensorioneural em todas as configurações audiométricas. A perda auditiva de grau moderado e a configuração audiométrica descendente foram mais frequentes nessa população, seguida da abrupta e horizontal.

Migração de solutos em basalto fraturado: quantificação experimental em laboratório e validação matemática

A avaliação do risco a contaminação e a escolha de técnicas de remediação de poluentes em aquíferos fraturados depende da quantificação dos fenômenos envolvidos no transporte de solutos. A geometria da fratura, usualmente caracterizada pela abertura, é o principal parâmetro que indiretamente controla o transporte nos aquíferos fraturados. A simplificação mais comum desse problema é assumir que as fraturas são um par de placas planas e paralelas, isto é, com uma abertura constante. No entanto, por causa do limitado número de trabalhos experimentais, não está esclarecida a adequabilidade do uso de uma abertura constante para simular o transporte conservativo em fraturas do Aquífero Serra Geral (ASG), Brasil. O objetivo deste trabalho é avaliar a influência da abertura de uma fratura natural do Aquífero Serra Geral sob o transporte conservativo de solutos. Uma amostra natural de basalto fraturado foi usada em um experimento hidráulico e de transporte de um traçador conservativo (escala de laboratório). O campo de abertura foi medido usando a técnica avançada, de alta resolução e tridimensional, chamada microtomografia computadorizada de raios-X. A concentração de traçador medida foi utilizada para validar uma solução analítica unidimensional da Equação de Advecção-dispersão (ADE). O desemprenho do ajuste da ADE às curvas de passagem experimentais foi avaliado para quatro diferentes tipos de aberturas constantes. Os resultados mostraram que o escoamento de água e o transporte de contaminantes pode ocorrer através de fraturas micrométricas, ocasionando, eventualmente, a contaminação do ASG. A abertura de balanço de massa é a única que pode ser chamada propriamente de \"abertura equivalente\". O uso de aberturas constantes na ADE não permitiu representar completamente o formato das curvas de passagem porque o campo de velocidade não é uniforme e intrinsicamente bidimensional. Portanto, na simulação do transporte deve-se incorporar a heterogeneidade da abertura da fratura.