Métodos para la estimación de la temperatura estática de la formación en yacimientos de petróleo y geotérmicos

Las empresas explotadoras de yacimientos de petróleo y geotérmicos diseñan y construyen sus instalaciones optimizándolas según dos características especialmente representativas: su temperatura y su presión. Por ello, tener información precisa sobre ambos aspectos lo antes posible puede permitir ahorros importantes en las inversiones en infraestructuras y también acortar los plazos de comienzo de explotación; ambos, aspectos esenciales en la rentabilidad del yacimiento. La temperatura estática de una formación es la temperatura existente antes del inicio de las tareas de perforación del yacimiento. Las operaciones de perforación tienen como objetivo perforar, evaluar y terminar un pozo para que produzca de una forma rentable. Durante estas tareas se perturba térmicamente la formación debido a la duración de la perforación, la circulación de fluidos de refrigeración, la diferencia de temperaturas entre la formación y el fluido que perfora, el radio, la difusividad térmica de las formaciones y la tecnología de perforación [7]. La principal herramienta para determinar la temperatura estática de formación son las medidas de temperaturas en los pozos. Estas medidas de temperatura se realizan una vez cesados los trabajos de perforación y circulación de fluidos. El conjunto de medidas de temperatura obtenidas es una serie de valores ascendentes, pero no llega a alcanzar el valor de la temperatura estática de la formación. Para estimar esta temperatura se plantean las ecuaciones diferenciales que describen los fenómenos de transmisión de calor en el yacimiento [14]. Estos métodos se emplean tanto en yacimientos geotérmicos como en yacimientos de petróleo indistintamente [16]. Cada uno de ellos modela el problema de transmisión de calor asociado de una forma distinta, con hipótesis simplificadoras distintas. Se ha comprobado que la aplicación de los distintos métodos a un yacimiento en concreto presenta discrepancias en cuanto a los resultados obtenidos [2,16]. En muchos de los yacimientos no se dispone de información sobre los datos necesarios para aplicar el método de estimación. Esto obliga a adoptar una serie de hipótesis que comprometen la precisión final del método utilizado. Además, puede ocurrir que el método elegido no sea el que mejor reproduce el comportamiento del yacimiento. En el presente trabajo se han analizado los distintos métodos de cálculo. De la gran variedad de métodos existentes [9] se han seleccionado los cuatro más representativos [2,7]: Horner (HM), flujo radial y esférico (SRM), de las dos medidas (TLM) y de fuente de calor cilíndrica (CSM). Estos métodos se han aplicado a una serie de yacimientos de los cuales se conoce la temperatura estática de la formación. De cada yacimiento se disponía de datos tanto de medidas como de las propiedades termo-físicas. Estos datos, en ocasiones, eran insuficientes para la aplicación de todos los métodos, lo cual obligo a adoptar hipótesis sobre los mismos. Esto ha permitido evaluar la precisión de cada método en cada caso. A partir de los resultados obtenidos se han formulado una colección de criterios que permiten seleccionar qué método se adapta mejor para el caso de un yacimiento concreto, de tal manera que el resultado final sea lo más preciso posible. Estos criterios se han fijado en función de las propiedades termo-físicas del yacimiento, que determinan el tipo de yacimiento, y del grado de conocimiento que se dispone del yacimiento objeto de estudio.

Estimación del riesgo sísmico en puerto príncipe (Haití)

Tras el devastador terremoto del 12 de enero de 2010 en Puerto Príncipe, Haití, las autoridades locales, numerosas ONGs y organismos nacionales e internacionales están trabajando en el desarrollo de estrategias para minimizar el elevado riesgo sísmico existente en el país. Para ello es necesario, en primer lugar, estimar dicho riesgo asociado a eventuales terremotos futuros que puedan producirse, evaluando el grado de pérdidas que podrían generar, para dimensionar la catástrofe y actuar en consecuencia, tanto en lo referente a medidas preventivas como a adopción de planes de emergencia. En ese sentido, este Trabajo Fin de Master aporta un análisis detallado del riesgo sísmico asociado a un futuro terremoto que podría producirse con probabilidad razonable, causando importantes daños en Puerto Príncipe. Se propone para ello una metodología de cálculo del riesgo adaptada a los condicionantes de la zona, con modelos calibrados empleando datos del sismo de 2010. Se ha desarrollado en el marco del proyecto de cooperación Sismo-Haití, financiado por la Universidad Politécnica de Madrid, que comenzó diez meses después del terremoto de 2010 como respuesta a una petición de ayuda del gobierno haitiano. El cálculo del riesgo requiere la consideración de dos inputs: la amenaza sísmica o movimiento esperado por el escenario definido (sismo de cierta magnitud y localización) y los elementos expuestos a esta amenaza (una clasificación del parque inmobiliario en diferentes tipologías constructivas, así como su vulnerabilidad). La vulnerabilidad de estas tipologías se describe por medio de funciones de daño: espectros de capacidad, que representan su comportamiento ante las fuerzas horizontales motivadas por los sismos, y curvas de fragilidad, que representan la probabilidad de que las estructuras sufran daños al alcanzar el máximo desplazamiento horizontal entre plantas debido a la mencionada fuerza horizontal. La metodología que se propone especifica determinadas pautas y criterios para estimar el movimiento, asignar la vulnerabilidad y evaluar el daño, cubriendo los tres estados del proceso. Por una parte, se consideran diferentes modelos de movimiento fuerte incluyendo el efecto local, y se identifican los que mejor ajustan a las observaciones de 2010. Por otra se clasifica el parque inmobiliario en diferentes tipologías constructivas, en base a la información extraída en una campaña de campo y utilizando además una base de datos aportada por el Ministerio de Obras Públicas de Haití. Ésta contiene información relevante de todos los edificios de la ciudad, resultando un total de 6 tipologías. Finalmente, para la estimación del daño se aplica el método capacidad-demanda implementado en el programa SELENA (Molina et al., 2010). En primer lugar, utilizado los datos de daño del terremoto de 2010, se ha calibrado el modelo propuesto de cálculo de riesgo sísmico: cuatro modelos de movimiento fuerte, tres modelos de tipo de suelo y un conjunto de funciones de daño. Finalmente, con el modelo calibrado, se ha simulado un escenario sísmico determinista correspondiente a un posible terremoto con epicentro próximo a Puerto Príncipe. Los resultados muestran que los daños estructurales serán considerables y podrán llevar a pérdidas económicas y humanas que causen un gran impacto en el país, lo que pone de manifiesto la alta vulnerabilidad estructural existente. Este resultado será facilitado a las autoridades locales, constituyendo una base sólida para toma de decisiones y adopción de políticas de prevención y mitigación del riesgo. Se recomienda dirigir esfuerzos hacia la reducción de la vulnerabilidad estructural - mediante refuerzo de edificios vulnerables y adopción de una normativa sismorresistente- y hacia el desarrollo de planes de emergencia. Abstract After the devastating 12 January 2010 earthquake that hit the city of Port-au-Prince, Haiti, strategies to minimize the high seismic risk are being developed by local authorities, NGOs, and national and international institutions. Two important tasks to reach this objective are, on the one hand, the evaluation of the seismic risk associated to possible future earthquakes in order to know the dimensions of the catastrophe; on the other hand, the design of preventive measures and emergency plans to minimize the consequences of such events. In this sense, this Master Thesis provides a detailed estimation of the damage that a possible future earthquake will cause in Port-au-Prince. A methodology to calculate the seismic risk is proposed, adapted to the study area conditions. This methodology has been calibrated using data from the 2010 earthquake. It has been conducted in the frame of the Sismo-Haiti cooperative project, supported by the Technical University of Madrid, which started ten months after the 2010 earthquake as an answer to an aid call of the Haitian government. The seismic risk calculation requires two inputs: the seismic hazard (expected ground motion due to a scenario earthquake given by magnitude and location) and the elements exposed to the hazard (classification of the building stock into building typologies, as well as their vulnerability). This vulnerability is described through the damage functions: capacity curves, which represent the structure performance against the horizontal forces caused by the seisms; and fragility curves, which represent the probability of damage as the structure reaches the maximum spectral displacement due to the horizontal force. The proposed methodology specifies certain guidelines and criteria to estimate the ground motion, assign the vulnerability, and evaluate the damage, covering the whole process. Firstly, different ground motion prediction equations including the local effect are considered, and the ones that have the best correlation with the observations of the 2010 earthquake, are identified. Secondly, the classification of building typologies is made by using the information collected during a field campaign, as well as a data base provided by the Ministry of Public Works of Haiti. This data base contains relevant information about all the buildings in the city, leading to a total of 6 different typologies. Finally, the damage is estimated using the capacity-spectrum method as implemented in the software SELENA (Molina et al., 2010). Data about the damage caused by the 2010 earthquake have been used to calibrate the proposed calculation model: different choices of ground motion relationships, soil models, and damage functions. Then, with the calibrated model, a deterministic scenario corresponding to an epicenter close to Port-au-Prince has been simulated. The results show high structural damage, and therefore, they point out the high structural vulnerability in the city. Besides, the economic and human losses associated to the damage would cause a great impact in the country. This result will be provided to the Haitian Government, constituting a scientific base for decision making and for the adoption of measures to prevent and mitigate the seismic risk. It is highly recommended to drive efforts towards the quality control of the new buildings -through reinforcement and construction according to a seismic code- and the development of emergency planning.

Estimación de la frecuencia y duración del riego mediante el uso de sensores de humedad

Este trabajo presenta los resultados de la campaña de riego 2010 realizada en la finca Villavieja, ubicada en Tordesillas (Valladolid) y perteneciente a la asociación para la investigación y la mejora del cultivo de remolacha azucarera (AIMCRA). El objetivo general es determinar el tiempo y frecuencia de riego óptimos, para mejorar así la eficiencia del riego y reducir el coste de los agricultores. Para ello se considera el manejo del riego tradicional que aplica el regante de la zona al cultivo de remolacha azucarera y se mide el contenido de humedad en puntos determinados de la finca y la distribución del agua en el suelo, a partir de esas medidas. Además, se realizaron varias evaluaciones de las unidades de riego a lo largo de la campaña. Un sensor de frecuencia FDR (EnviroScan) se utilizó para medir el contenido de humedad en un perfil de suelo de la finca. Previamente se determinó la textura del suelo. Los regantes de la zona conocían, con cierta aproximación, las necesidades hídricas semanales del cultivo, a partir de la información de la evapotranspiración semanal de la remolacha, recibida vía mensajes cortos de texto de telefonía móvil (SMS). Ésta se determinada a partir de las estación meteorológicas más cercanas (pertenecientes a la red de Inforiego) y de la red de tanques de evaporación repartidos por las diferentes zonas. Sin embargo, y debido a la necesidad de reducir costes, muchos regantes no se ajustan a dichas necesidades y aplican un calendario de riego ajustándose a la tarifa eléctrica mínima . Esta situación puede provocar estrés hídrico, durante algunos días de la campaña y, como consecuencia, riegos excesivamente largos que intentan compensar el estrés acumulado. Los suelos, principalmente arenosos mostraron valores altos de la infiltración; el tiempo medio de riego, de unas 6 horas, llegó incluso a duplicar el necesario para la saturación de los primeros 60 cm de suelo. Además, la frecuencia de riego no fue siempre la adecuada y, en varias ocasiones, se observaron visualmente síntomas de estrés hídrico en el cultivo.

Estimación de la infiltración en suelos mediante la utilización del"AHFO method"

Se ha caracterizado la infiltración de un suelo colocado en una columna de suelo de metacrilato, de base hexagonal de diagonal 1 m y 0,6 m de alto, con una densidad aparente de 1,5 g/cm3. El procedimiento utilizado ha sido la fibra óptica y el método denominado “Active Heating pulse method with Fiber Optic temperature sensing” AHFO method, que consiste en emitir un pulso óptico con láser y medir en el tiempo la señal reflejada, de baja intensidad, en diferentes puntos de la fibra óptica. Del espectro de luz reflejada solamente un rango de frecuencias específico, determinado por análisis de frecuencia, se correlaciona con la temperatura. La precisión en la medida es de ± 0,1ºC en una distancia de ± 12,5 cm. En el interior de la columna se colocó el cable de fibra óptica formando tres hélices concéntricas separadas 20 cm , 40 cm y 60 cm del centro. Asimismo, se cubrió la superficie del suelo con una altura media de agua que osciló entre 1,5 a 2,5 cm a lo largo de los 140 min que duró el proceso de calentamiento del cable. El incremento de temperatura antes y después del calentamiento se utilizó para determinar la infiltración instantánea a partir de la expresión de Perzlmaeir et al (2004) y de los números adimensional de Nusselt y Prandtl. Considerando los errores inherentes al procedimiento de cálculo, los resultados muestran que el AHFO method es una herramienta útil en el estudio de la variabilidad espacial de la infiltración en el suelo que permite, además, determinar su valor. Asimismo, muestra su potencial para incluir dichas estimaciones en la calibración de modelos relacionados con la gestión de recursos hídricos.

Modelo de estimación de costes para tolerancias de diseño en la ingeniería de los procesos de producción

La disponibilidad los valores coste-tolerancia de un proceso de producción es un aspecto de interés desde el punto de vista de la Ingeniería de Proyectos en general y de la Ingeniería de Procesos de Producción en particular. Tanto en la fase de planificación de un proyecto como en las etapas de diseño de cada uno de los procesos de producción que posibilitarán su ejecución, las especificaciones del producto deben estar contenidas en un intervalo de tolerancia cuyo coste de obtención es tanto mayor cuanto menor es el margen de variación admisible. Por este motivo, el conocimiento del coste asociado a cada tolerancia es una herramienta de gran utilidad para el ingeniero de proyectos. A pesar de su interés, las curvas de coste-tolerancia se manejan con menor profusión de la que sería deseable, dada la gran interdependencia existente entre ellas y la localización espacio-temporal de los recursos a las que se asignan. El presente trabajo propone un método que permite establecer una relación coste-tolerancia a los procesos de producción de la Ingeniería de Proyectos, basada de la variabilidad derivada de los recursos de cada proceso.

Comparación de dos métodos para la estimación de los daños por sequía en pastos de dehesa / Comparative methods for assessment of dehesa pasture drought insurance in Spain

El Sistema de Seguros Agrarios con el Seguro de cobertura de los daños por sequía en los pastos aprovechados por el ganado en régimen extensivo (línea de seguro 133) aplica la teledetección mediante un índice de vegetación (NDVI), con el fin de solucionar los problemas de peritación que surgen cuando se tiene que determinar la cantidad y calidad del pasto afectado por la sequía. Por ello el seguro de cobertura de los daños por sequía en pastos es el principal instrumento para hacer frente al gasto que supone la necesidad de suplemento de alimentación del ganado reproductor debido a la sequía. En las comarcas de Vitigudino, Trujillo y Valle de los Pedroches (España) se comparó la evolución del seguro de sequía en pastos desde 2006 a 2010 con un modelo matemático de crecimiento del pasto en función de las variables ecofisiológicas y ambientales. Sumadas las decenas de sequía extrema y sequía leve, el modelo matemático contabilizó un número mayor de decenas que las proporcionadas por Agroseguro. La recomendación es comparar las curvas de crecimiento del pasto con las curvas de evolución del NDVI, para ajustar ambos modelos

Estimación de la Reflectividad del Multitrayecto en Entornos Urbanos en Radares de Apertura Sintética

Synthetic Aperture Radar (SAR) images a target region reflectivity function in the multi-dimensional spatial domain of range and cross-range. SAR synthesizes a large aperture radar in order to achieve finer azimuth resolution than the one provided by any on-board real antenna. Conventional SAR techniques assume a single reflection of transmitted waveforms from targets. Nevertheless, today¿s new scenes force SAR systems to work in urban environments. Consequently, multiple-bounce returns are added to direct-scatter echoes. We refer to these as ghost images, since they obscure true target image and lead to poor resolution. By analyzing the quadratic phase error (QPE), this paper demonstrates that Earth¿s curvature influences the defocusing degree of multipath returns. In addition to the QPE, other parameters such as integrated sidelobe ratio (ISLR), peak sidelobe ratio (PSLR), contrast and entropy provide us with the tools to identify direct-scatter echoes in images containing undesired returns coming from multipath.

Estimación del error de discretización con principios variacionales multicampo: I - Elasticidad

Desarrollos recientes para encajar dentro de un marco variacional la llamada Formulación Libre sugieren la posibilidad de introducir un nuevo tipo de estimador de error para cálculos por elementos finitos. Este estimador se basa en una diferencia entre ciertos funcionales multicampo, que toman el mismo valor para la solución exacta del problema. En el presente articulo, dividido en dos partes, se introduce la formulación del estimador para problemas de elasticidad y de flexión de placas según las hipótesis clásicas de Kirchhoff. Se presentan también algunos ejeinplos para dar idea de los comportamientos numéricos observados.

Estimación de ahorro por inspección del regulador del mercado energético a las instalaciones de energías renovables

La Administración General del Estado, a través de la Comisión Nacional de la Energía, y en colaboración con los órganos competentes de las Comunidades Autónomas correspondientes, realiza inspecciones periódicas y aleatorias cada año, sobre las instalaciones de generación eléctrica, siguiendo los criterios de elección e indicaciones que la Secretaria General de la Energía del Ministerio de Industria, Turismo y Comercio imponga en cada caso. En este proyecto se analizan las distintas incidencias de las distintas tecnologías estudiadas, así como las consecuencias jurídicas y económicas. También se ofrecen las distintas propuestas de sanción como consecuencia de estas incidencias. El objeto del presente proyecto es el estudio de la reducción de costes al estado si el organismo competente (Comisión Nacional de Energía) liquidase a las diferentes instalaciones en función de las irregularidades detectadas como consecuencia de las inspecciones realizadas durante el plan 2011/2012, así como una estimación de dicho ahorro al extrapolar dichos resultados al total de instalaciones ABSTRACT The Central Government, through the National Energy Commission, in cooperation with the competent bodies of the Autonomous Communities concerned, realizes periodic and random inspections every year on electricity generation facilities, following the criteria of choice and indications that the General Secretary of Energy from the Ministry of Industry, Tourism and Commerce impose in each case. This project analyses the different incidences of the different technologies studied as well as the legal and economic consequences. It is also offered different proposals of sanction as a result of these incidents. The purpose of this project is to study the reduction of costs to the state if the competent agency (CNE) decided to close the different facilities depending on the irregularities as a result of inspections conducted during the plan 2011/2012 and as an estimate of the savings to extrapolate these results to all the facilities.

Estimación del error de discretización con principios variacionales multicampo: II. Flexión de placas

Desarrollos recientes para encajar dentro de un marco variacional la llamada Formulación Libre sugieren la posibilidad de introducir un nuevo tipo de estimador de error para cálculos por elementos finitos. Este estimador se basa en una diferencia entre ciertos funcionales multicampo, que toman el mismo valor para la solución exacta del problema. En el presente artículo, dividido en dos partes, se introduce la formulación del estimador para problemas de elasticidad y de flexión de placas según las hipótesis clásicas de Kirchhoff. Se presentan también algunos ejemplos para dar idea de los comportamientos numéricos observados.

Contraste de las formulaciones existentes para la estimación de la tasa de rebase en diques de abrigo portuario

La investigación que se resume en este artículo, resultante de diversos trabajos realizados en el CEDEX [10, 12], tiene por objeto contrastar las formulaciones desarrolladas por diferentes autores en materia de rebases en distintas tipologías de diques de abrigo. Para ello, se realizó 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 dos ensayos realizados en el Centro de Estudios de Puertos y Costas del CEDEX, en Cabo Prioriño, Ferrol y en la alineación tercera de Melilla. The research that is summarized in this article, resultant of several studies carried out in the CEDEX [10, 12], is focused in the existing methods to estimate overtopping rates developed by different authors. For that, the summary was carry out firstly we collected and analyzed the existing formulae, particularly for rubble mound and vertical breakwaters. These formulations was compared with the results obtained in two hydraulic model tests of the Hydraulic Research Laboratory (the Center of Studies of Ports and Coasts of the CEDEX, Madrid, Spain). Keywords: overtopping rates, rubble mound breakwater, vertical breakwater, crown wall, run-up.

La estimación de la biodiversidad forestal en el Inventario Forestal Nacional. Aplicación en el IFN-4 en Galicia

Teniendo como punto de partida la Convención de la Diversidad biológica de 1992, desde la segunda mitad del siglo XX las perspectivas en política y gestión medioambiental han dado un giro sin precedentes. Variables intangibles, difíciles de medir, como es la biodiversidad, son cada vez más importantes para la sociedad. En consecuencia, enmarcado en políticas y requerimientos internacionales, desde el Tercer Inventario Forestal Nacional (IFN-3) se está aplicando una nueva metodología en continuo desarrollo que trata de estimar nuevas variables relacionadas con la biodiversidad. El objetivo más importante es poder analizar mediante estos nuevos indicadores la evolución del estado de la biodiversidad forestal de nuestros bosques a lo largo de los sucesivos IFN. Para llevar a cabo este reto, la toma de datos de campo en el inventario se centra principalmente en indicadores relativos a la composición florística y a la estructura de las masas forestales. En este artículo se describe la evolución y el desarrollo de la metodología llevada a cabo hasta el momento para la estimación de la biodiversidad forestal en el IFN. Además, se detallan los resultados derivados de su aplicación en las formaciones forestales más representativas de Galicia. Finalmente, se exponen los nuevos indicadores relacionados con la biodiversidad forestal en los que se está investigando, así como nuevas perspectivas de análisis.

Estimación de la resistencia in situ de un HAC mediante la resistencia de probetas testigo

El presente trabajo tiene como objetivo establecer algunos factores de corrección que afectan el resultado de la resistencia de las probetas testigo a fin de pronosticar la resistencia a compresión in situ. Para el cumplimiento de dicho objetivo se han fabricado 30 bloques de 500x500x1000 mm en Hormigón Autocompactante (HAC) cubriendo un amplio rango de resistencias de 20 MPa a 80 MPa. De estos bloques se han extraído 360 probetas testigo con variaciones en el diámetro (50 mm, 75 mm y 100 mm), en la dirección de extracción (horizontal y vertical) y en la condición de humedad secas y húmedas). Para el pronóstico de la resistencia in situ se ha considerado el principio de superposición de efectos establecido en normas como la ACI-214.4R-10. En este sentido, para las condiciones del ensayo, se han cuantificado mediante factores de corrección los efectos generados por el diámetro y la condición de humedad del testigo en el momento del ensayo y por el nivel resistente del hormigón. Además, se ha valorado de manera conjunta la influencia de las variables intrínsecas de las probetas testigo, tales como, microfisuración, efecto pared, condiciones de contorno, etc. Los resultados indicaron que algunos factores de corrección utilizados para pronosticar la resistencia in situ, inciden de forma muy diferente en los HAC con respecto a los hormigones vibrados tradicionalmente (HVT).

Aplicación de la metodología bootstrap para la estimación de la movilidad de furgonetas

En este trabajo se muestran los resultados de la aplicación de la metodología bootstrap a datos de 3369 encuestas realizadas en 2009 a nivel nacional entre conductores de furgonetas, para obtener datos de movilidad interurbana y total, según edad de los vehículos, uso, conductores y otras características de este tipo de vehículo. Se obtienen estimaciones puntuales e intervalos de confianza para la movilidad total de furgonetas, así como para los cuatro tipos de furgonetas según la clasificación realizada en el proyecto de referencia. Se comparan los resultados obtenidos con estimaciones alternativas realizadas con otras fuentes de datos para el mismo colectivo (encuestas realizadas en inspecciones en carretera realizadas por la ATGC de la DGT e inspecciones en ITV) y datos publicados por fuentes oficiales. Estos resultados de movilidad (en término de vehículo-kilómetro) se usarán para la estimación de ratios de accidentalidad en un estudio comparado con otros colectivos de vehículos.

Estimación de modelos borrosos y su aplicación al control óptimo

El objetivo de la presentación es dar a conocer los últimos trabajos del Grupo de Investigación sobre últimos trabajos del Grupo de Investigación sobre Control Borroso. • Obtención de modelos precisos de sistemas no lineales basados en sistemas borrosos – Mamdani – Takagi-Sugeno – Linealización • Generalización del método propuesto por T-S • Identificación iterativa basada en el Filtro de Kalman • Sistemas de control basados en el modelo TS obtenido – LQR