Aplicación de método de los elementos finitos a los problemas de amplificación dinámica de suelos

En el análisis de estructuras situadas en un emplazamiento de riesgo sísmico, es necesario conocer las característiias de un posible movimiento sísmico actuante, en la cota de cimentación de la estructura. La importancia del papel que juega el suelo circundante a la estructura en la modificación de las características del movimientos sísmico (principalmente en cuanto a amplitudes y contenido frecuencial), es bien conocida, denominándose amplificación sísmica a este fenómeno. En el caso de estructuras importantes, como Centrales Nucleares, presas, etc., o terrenos muy blandos, la presencia de la propia estructura también afecta de forma significativa al movimiento, dando lu gar al fenómeno de la interacción suelo-estructura. En esta comunicación se tratará únicamente del problema de la amplificación de las ondas sísmicas causantes del movimiento del suelo, y del análisis de la respuesta del mismo en ausencia de la estructura. Este movimiento del suelo denominado "movimiento del campo libre"(free field motion) puede, bien ser utilizado directamente en el análisis de la estructura, despreciando los efectos de la interacción, o emplearse como "input" en ciertas técnicas de análisis de interacción suelo-estructura.

Evaluación de los modificadores por comportamiento urbanos que afectan al daño en los terremotos. Aplicación al terremoto de Lorca.

Habitualmente, la estimación de la vulnerabilidad sísmica se centra en el comportamiento estructural de los edificios. Sólo algunas metodologías, como el proyecto Risk‐UE,consideran la influencia de otros factores no estructurales o urbanísticos, tales como el piso blando, la irregularidad en alzado, la irregularidad en planta, etc. Estos factores,denominados también modificadores por comportamiento, pueden tener una incidencia en el daño observado, y la confluencia de varios de ellos puede variar sustancialmente la vulnerabilidad. Los modificadores por comportamiento se han identificado de forma empírica, a través de la observación de patrones de daño típicos en terremotos, teniendo en cuenta las inspecciones visuales (ATC 21 1988, Benedetti y Petrini 1984, UNDP/UNIDO 1985) y otras propuestas (Coburn y Spence 1992). La puntuación del modificador ha sido dada por el conocimiento de expertos en terremotos tras analizar evaluaciones de vulnerabilidad anteriores y bases de datos del daño producido en edificios. En esta comunicación estudiaremos los modificadores que derivan de características urbanísticas. Esta línea de investigación considera que un parámetro modificador deriva de características urbanísticas si puede ser regulado en la Normativa Urbanística de un Plan General de Ordenación Urbana. Se realiza una descripción de cada modificador según cada metodología o investigador (Risk‐UE, Giovinazzi, Lantada y Feriche) y una comparativa entre las distintas ponderaciones de los modificadores. Este análisis nos permite poder tener una primera visión de la posible cuantificación de cada modificador y la tendencia que ha tenido la calibración desde el año 2003 con el proyecto Risk‐UE hasta el año 2012 con la tesis de Feriche. Finalmente se presentan los resultados del estudio exploratorio de los parámetros urbanísticos de tres zonas seleccionados de la ciudad de Lorca según el tipo de suelo en el que se encuentren y se indican aquellos parámetros que han podido influir en el daño provocado por el terremoto de mayo de 2011.

Edificio judicial en El Ejido

La sede judicial en El Ejido, edificio de planta trapezoidal de 47 m x 55 / 26 m, presenta dos volúmenes claramente diferenciados, unidos a nivel de sótano y conectados mediante pasarelas en niveles superiores. Una tercera pieza, de planta trapezoidal y fachada frontal de muro cortina de vidrio, vuela 8 m del plano vertical del edificio más alto. La fachada principal desempeña una función estructural de pantalla y está constituida por 9 filas de pilares prefabricados de hormigón armado entre las que se intercalan unos cordones horizontales de hormigón autocompactante, ejecutados in situ. Debido a la ubicación del edificio en una zona sísmica y el reducido número de apoyos horizontales, la fachada resulta potencialmente vulnerable. Para mitigar los elevados riesgos, particularmente durante las fases constructivas, se hacía necesaria una planificación rigurosa de todo el proceso sobre la base de un análisis específico de la interacción entre la estructura evolutiva de la fachada y la cimbra empleada en su construcción.

Evaluación del riesgo sísmico en España a escala municipal y su evolución temporal.

La acción producida por fenómenos naturales catastróficos en entornos urbanos supone un riesgo importante, pues conlleva innumerables pérdidas, tanto económicas como humanas. Los terremotos son uno de los fenómenos que mayor cantidad de pérdidas ha producido en todo el mundo, debido a su aleatoriedad y su complicada predicción exacta. Por este motivo, el conocimiento del riesgo sísmico de una zona es fundamental para la adopción de medidas de prevención conducentes a la mitigación del riesgo. En España residen cerca de 47.000.000 de habitantes, en un área de, aproximadamente, 500.000 km2. Casi 3.000.000 de edificios residenciales fueron construidos (suponiendo el 33% del total) cuando no existían los códigos de diseño antisísmico. En la actualidad, estos edificios antiguos coexisten con nuevas construcciones que han aplicado el diseño sismorresistente, especialmente en ciudades situadas en zonas sísmicamente activas, que pueden presentar un alto riesgo incluso bajo los efectos de un terremoto moderado. La evaluación del riesgo sísmico es un método de valorar los posibles daños que puede provocar una acción sísmica. Para su estimación, se precisaevaluar la peligrosidad sísmica de la zona,y la vulnerabilidad de los elementos expuestos. Si bien la peligrosidad responde a un proceso natural que no se puede controlar, la vulnerabilidad sí se puede reducir (por ejemplo, ejecutando medidas de construcción sismorresistente). En definitiva, el conocimiento del riesgo sísmico de una zona puede servir para aumentar la conciencia sísmica en las ciudades e impulsar el desarrollo de planes de acción necesarios para la reducción del riesgo sísmico. La presente memoria expone de forma detallada todos los aspectos técnicos y metodológicos en los que nos hemos servido para la realización de este Trabajo Final de Grado.

Diseño sísmico de instalaciones nucleares y obras públicas

Aquí se comentarán algunos criterios que se deben tener en cuenta en el diseño de obras singulares. En particular, se pondrá un especial enfásis en la elaboración de los datos precisos para un análisis sísmico, aspecto de excepcional importancia y que contiene numerosas cuestiones abiertas, en espera de acuerdo entre los expertos. Se finalizará esta exposición con un breve resumen de los métodos y modelos de cálculo más utilizados en el análisis y su rango de aplicación

Un método de obtención del terremoto de diseño en un emplazamiento

Se describe una metodología orientada a la determinación del período de retorno de cada intensidad sísmica sentida en un emplazamiento dado. El procedimiento considera la regionalización sismotectónica del área total susceptible de producir terremotos con efectos sentidos apreciables en el emplazamiento. El modelo de distribución espacial de aparición de terremotos se simula mediante una adecuada definición de las distintas unidades sismotectónicas, en las que se divide el área total y para la que esta distribución se considera uniforme. La aparición en el tiempo de los seismos se supone sigue un modelo de Poisson con un parámetro variable aleatoria, descrito mediante una distribución conjugada garnma-1, cuyas estadísticos se ajustan de acuerdo con el teorema de Bayes a partir de la información histórica existente. La integración de los dos modelos anteriores permlte deducir la distribución de probabilidad de aparición de un movimiento sísmico de intensidad determinada en el emplazamiento que se estudia

Movimientos fuertes y espectros para el Centro Norte de Quito asociados a las fallas ciegas inversas.

En 1587 un sismo de magnitud 6.4 asociado a las fallas ciegas de Quito, causó gran daño en la ciudad; esto es muy preocupante toda vez que existen estudios que indican que el período de recurrencia es de aproximadamente 200 años. Lo que implica que al 2014 debieron haberse producido dos sismos de igual magnitud, que gracias a Dios no ha sucedido, ya que las consecuencias habrían sido desastrosas porque buena parte de la ciudad se encuentra sobre estas fallas. Al no haberse registrado sismos de la magnitud indicada en los cuatro últimos siglos se tiene una gran acumulación de energía toda vez que la tasa de movimiento de estas fallas se encuentra entre los 3 y 4 mm al año. De tal manera que existe una gran probabilidad de tener un sismo muy fuerte. En estas condiciones es muy importante obtener espectros de respuesta elástica que se generarían por un sismo asociado a alguno de los segmentos de las fallas ciegas inversas de Quito y que mejor hacerlo con modelos de movimientos fuertes que han sido publicados en los últimos años y reconocidos a nivel mundial por la gran base de datos con la que fueron obtenidos. En este artículo se van a presentar tres modelos de movimientos fuertes y son los desarrollados por: Abrahamson, Silva y Kamai (2013); Campbell y Borzognia (2013); y el modelo de Zhao et al. (2006). Para cada uno de estos modelos se indican la base de datos y las ecuaciones que definen las ramas del espectro. Se destaca que estos modelos sirven también para encontrar leyes de atenuación del movimiento del suelo. Posteriormente como una aplicación se obtendrá espectros para cinco Parroquias del Centro Norte de Quito, a saber: Rumipamba, Belisario Quevedo, Mariscal Sucre, Iñaquito y Jipijapa; asociados al sismo máximo probable generados en el segmento de falla que está más cerca y es el denominado Ilumbisí-La Bota, con cada uno de los modelos indicados en el párrafo anterior y se encontrarán espectros ponderados para un nivel de confiabilidad del 84%. Las ordenadas espectrales que se hallan son mayores a las ordenadas espectrales que reporta la Norma Ecuatoriana de la Construcción de 2011 y los que se obtienen con los factores de sitio hallados en el estudio de Microzonificación Sísmica de Quito por ERN en el 2012.

Real-Time Dynamic PGD Calculation of Non-Linear Soil Behavior

El estudio sísmico en los últimos 50 años y el análisis del comportamiento dinámico del suelo revelan que el comportamiento del suelo es altamente no lineal e histéretico incluso para pequeñas deformaciones. El comportamiento no lineal del suelo durante un evento sísmico tiene un papel predominante en el análisis de la respuesta de sitio. Los análisis unidimensionales de la respuesta sísmica del suelo son a menudo realizados utilizando procedimientos lineales equivalentes, que requieren generalmente pocos parámetros conocidos. Los análisis de respuesta de sitio no lineal tienen el potencial para simular con mayor precisión el comportamiento del suelo, pero su aplicación en la práctica se ha visto limitada debido a la selección de parámetros poco documentadas y poco claras, así como una inadecuada documentación de los beneficios del modelado no lineal en relación al modelado lineal equivalente. En el análisis del suelo, el comportamiento del suelo es aproximado como un sólido Kelvin-Voigt con un módulo de corte elástico y amortiguamiento viscoso. En el análisis lineal y no lineal del suelo se están considerando geometrías y modelos reológicos más complejos. El primero está siendo dirigido por considerar parametrizaciones más ricas del comportamiento linealizado y el segundo mediante el uso de multi-modo de los elementos de resorte-amortiguador con un eventual amortiguador fraccional. El uso del cálculo fraccional está motivado en gran parte por el hecho de que se requieren menos parámetros para lograr la aproximación exacta a los datos experimentales. Basándose en el modelo de Kelvin-Voigt, la viscoelasticidad es revisada desde su formulación más estándar a algunas descripciones más avanzada que implica la amortiguación dependiente de la frecuencia (o viscosidad), analizando los efectos de considerar derivados fraccionarios para representar esas contribuciones viscosas. Vamos a demostrar que tal elección se traduce en modelos más ricos que pueden adaptarse a diferentes limitaciones relacionadas con la potencia disipada, amplitud de la respuesta y el ángulo de fase. Por otra parte, el uso de derivados fraccionarios permite acomodar en paralelo, dentro de un análogo de Kelvin-Voigt generalizado, muchos amortiguadores que contribuyen a aumentar la flexibilidad del modelado para la descripción de los resultados experimentales. Obviamente estos modelos ricos implican muchos parámetros, los asociados con el comportamiento y los relacionados con los derivados fraccionarios. El análisis paramétrico de estos modelos requiere técnicas numéricas eficientemente capaces de simular comportamientos complejos. El método de la Descomposición Propia Generalizada (PGD) es el candidato perfecto para la construcción de este tipo de soluciones paramétricas. Podemos calcular off-line la solución paramétrica para el depósito de suelo, para todos los parámetros del modelo, tan pronto como tales soluciones paramétricas están disponibles, el problema puede ser resuelto en tiempo real, porque no se necesita ningún nuevo cálculo, el solucionador sólo necesita particularizar on-line la solución paramétrica calculada off-line, que aliviará significativamente el procedimiento de solución. En el marco de la PGD, parámetros de los materiales y los diferentes poderes de derivación podrían introducirse como extra-coordenadas en el procedimiento de solución. El cálculo fraccional y el nuevo método de reducción modelo llamado Descomposición Propia Generalizada han sido aplicado en esta tesis tanto al análisis lineal como al análisis no lineal de la respuesta del suelo utilizando un método lineal equivalente. ABSTRACT Studies of earthquakes over the last 50 years and the examination of dynamic soil behavior reveal that soil behavior is highly nonlinear and hysteretic even at small strains. Nonlinear behavior of soils during a seismic event has a predominant role in current site response analysis. One-dimensional seismic ground response analysis are often performed using equivalent-linear procedures, which require few, generally well-known parameters. Nonlinear analyses have the potential to more accurately simulate soil behavior, but their implementation in practice has been limited because of poorly documented and unclear parameter selection, as well as inadequate documentation of the benefits of nonlinear modeling relative to equivalent linear modeling. In soil analysis, soil behaviour is approximated as a Kelvin-Voigt solid with a elastic shear modulus and viscous damping. In linear and nonlinear analysis more complex geometries and more complex rheological models are being considered. The first is being addressed by considering richer parametrizations of the linearized behavior and the second by using multi-mode spring-dashpot elements with eventual fractional damping. The use of fractional calculus is motivated in large part by the fact that fewer parameters are required to achieve accurate approximation of experimental data. Based in Kelvin-Voigt model the viscoelastodynamics is revisited from its most standard formulation to some more advanced description involving frequency-dependent damping (or viscosity), analyzing the effects of considering fractional derivatives for representing such viscous contributions. We will prove that such a choice results in richer models that can accommodate different constraints related to the dissipated power, response amplitude and phase angle. Moreover, the use of fractional derivatives allows to accommodate in parallel, within a generalized Kelvin-Voigt analog, many dashpots that contribute to increase the modeling flexibility for describing experimental findings. Obviously these rich models involve many parameters, the ones associated with the behavior and the ones related to the fractional derivatives. The parametric analysis of all these models require efficient numerical techniques able to simulate complex behaviors. The Proper Generalized Decomposition (PGD) is the perfect candidate for producing such kind of parametric solutions. We can compute off-line the parametric solution for the soil deposit, for all parameter of the model, as soon as such parametric solutions are available, the problem can be solved in real time because no new calculation is needed, the solver only needs particularize on-line the parametric solution calculated off-line, which will alleviate significantly the solution procedure. Within the PGD framework material parameters and the different derivation powers could be introduced as extra-coordinates in the solution procedure. Fractional calculus and the new model reduction method called Proper Generalized Decomposition has been applied in this thesis to the linear analysis and nonlinear soil response analysis using a equivalent linear method.

Diseño de un sistema regional de alerta de tsunamis como parte integral de la operación portuaria

En este trabajo se presenta el desarrollo de una metodología para obtener un universo de funciones de Green y el algoritmo correspondiente, para estimar la altura de tsunamis a lo largo de la costa occidental de México en función del momento sísmico y de la extensión del área de ruptura de sismos interplaca localizados entre la costa y la Trinchera Mesoamericana. Tomando como caso de estudio el sismo ocurrido el 9 de octubre de 1995 en la costa de Jalisco-Colima, se estudiaron los efectos del tsunami originados en la hidrodinámica del Puerto de Manzanillo, México, con una propuesta metodológica que contempló lo siguiente: El primer paso de la metodología contempló la aplicación del método inverso de tsunamis para acotar los parámetros de la fuente sísmica mediante la confección de un universo de funciones de Green para la costa occidental de México. Tanto el momento sísmico como la localización y extensión del área de ruptura de sismos se prescribe en segmentos de planos de falla de 30 X 30 km. A cada uno de estos segmentos del plano de falla corresponde un conjunto de funciones de Green ubicadas en la isobata de 100 m, para 172 localidades a lo largo de la costa, separadas en promedio 12 km entre una y otra. El segundo paso de la metodología contempló el estudio de la hidrodinámica (velocidades de las corrientes y niveles del mar en el interior del puerto y el estudio del runup en la playa) originada por el tsunami, la cual se estudió en un modelo hidráulico de fondo fijo y en un modelo numérico, representando un tsunami sintético en la profundidad de 34 m como condición inicial, el cual se propagó a la costa con una señal de onda solitaria. Como resultado de la hidrodinámica del puerto de Manzanillo, se realizó un análisis de riesgo para la definición de las condiciones operativas del puerto en términos de las velocidades en el interior del mismo, y partiendo de las condiciones iniciales del terremoto de 1995, se definieron las condiciones límites de operación de los barcos en el interior y exterior del puerto. In this work is presented the development of a methodology in order to obtain a universe of Green's functions and the corresponding algorithm in order to estimate the tsunami wave height along the west coast of Mexico, in terms of seismic moment and the extent of the area of the rupture, in the interplate earthquakes located between the coast and the Middle America Trench. Taking as a case of study the earthquake occurred on October 9, 1995 on the coast of Jalisco-Colima, were studied the hydrodynamics effects of the tsunami caused in the Port of Manzanillo, Mexico, with a methodology that contemplated the following The first step of the methodology contemplated the implementation of the tsunami inverse method to narrow the parameters of the seismic source through the creation of a universe of Green's functions for the west coast of Mexico. Both the seismic moment as the location and extent of earthquake rupture area prescribed in segments fault planes of 30 X 30 km. Each of these segments of the fault plane corresponds a set of Green's functions located in the 100 m isobath, to 172 locations along the coast, separated on average 12 km from each other. The second step of the methodology contemplated the study of the hydrodynamics (speed and directions of currents and sea levels within the port and the study of the runup on the beach Las Brisas) caused by the tsunami, which was studied in a hydraulic model of fix bed and in a numerical model, representing a synthetic tsunami in the depth of 34 m as an initial condition which spread to the coast with a solitary wave signal. As a result of the hydrodynamics of the port of Manzanillo, a risk analysis to define the operating conditions of the port in terms of the velocities in the inner and outside of the port was made, taken in account the initial conditions of the earthquake and tsunami ocurred in Manzanillo port in 1995, were defined the limits conditions of operation of the ships inside and outside the port.

Diseño y evaluación de un procedimiento on-line para la biorremediación de aguas radiactivas en centrales nucleares

Recientemente se ha demostrado la existencia de microorganismos en las piscinas de almacenamiento de combustible nuclear gastado en las centrales nucleares utilizando técnicas convencionales de cultivo en el laboratorio. Estudios posteriores han puesto de manifiesto que los microorganismos presentes eran capaces de colonizar las paredes de acero inoxidable de las piscinas formando biopelículas. Adicionalmente se ha observado la capacidad de estas biopelículas de retener radionúclidos, lo que hace pensar en la posibilidad de utilizarlas en la descontaminación de las aguas radiactivas de las piscinas. En la presente tesis se plantea conocer más profundamente la biodiversidad microbiana de las biopelículas utilizando técnicas de biología molecular como la clonación, además de desarrollar un sistema de descontaminación a escala piloto con el objetivo de valorar si el proceso podría resultar escalable a nivel industrial. Para ello se diseñaron y fabricaron dos biorreactores en acero inoxidable compatibles con las condiciones específicas de seguridad sísmica y protección frente a la radiación en la zona controlada de una central nuclear. Los biorreactores se instalaron en la Central Nuclear de Cofrentes (Valencia) en las proximidades de las piscinas de almacenamiento de combustible nuclear gastado y precediendo a las resinas de intercambio iónico, de forma que reciben el agua de las piscinas permitiendo el análisis in situ de la radiación eliminada del agua de las mismas. Se conectó una lámpara de luz ultravioleta a uno de los biorreactores para poder comparar el desarrollo de bipelículas y la retención de radiactividad en ambas condiciones. En estos biorreactores se introdujeron ovillos de acero inoxidable y de titanio que se extrajeron a diversos tiempos, hasta 635 días para los ovillos de acero inoxidable y hasta 309 días para los ovillos de titanio. Se analizaron las biopelículas desarrolladas sobre los ovillos por microscopía electrónica de barrido y por microscopía de epifluorescencia. Se extrajo el ADN de las biopelículas y, tras su clonación, se identificaron los microorganismos por técnicas independientes de cultivo. Asimismo se determinó por espectrometría gamma la capacidad de las biopelículas para retener radionúclidos. Los microorganismos radiorresistentes identificados pertenecen a los grupos filogenéticos Alpha-proteobacteria, Gamma-proteobacteria, Actinobacteria, Deinococcus-Thermus y Bacteroidetes. Las secuencias de estos microorganismos se han depositado en el GenBank con los números de acceso KR817260-KR817405. Se ha observado una distribución porcentual ligeramente diferente en relación con el tipo de biorreactor. Las biopelículas han retenido fundamentalmente radionúclidos de activación. La suma de Co-60 y Mn-54 ha llegado en ocasiones al 97%. Otros radionúclidos retenidos han sido Cr-51, Co-58, Fe-59, Zn-65 y Zr-95. Se sugiere un mecanismo del proceso de retención de radionúclidos relacionado con el tiempo de formación y desaparición de las biopelículas. Se ha valorado que el proceso escalable puede ser económicamente rentable. ABSTRACT The existence of microorganisms in spent nuclear fuel pools has been demonstrated recently in nuclear power plants by using conventional microbial techniques. Subsequent studies have revealed that those microorganisms were able to colonize the stainless steel pool walls forming biofilms. Additionally, it has been observed the ability of these biofilms to retain radionuclides, which suggests the possibility of using them for radioactive water decontamination purposes. This thesis presents deeper knowledge of microbial biofilms biodiversity by using molecular biology techniques such as cloning, and develops a decontamination system on a pilot scale, in order to assess whether the process could be scalable to an industrial level. Aiming to demonstrate this was feasible, two stainless steel bioreactors were designed and manufactured, both were compatible with seismic and radiation protection standards in the controlled zone of a nuclear plant. These bioreactors were installed in the Cofrentes Nuclear Power Plant (Valencia) next to the spent nuclear fuel pools and preceding (upstream) ion exchange resins. This configuration allowed the bioreactors to receive water directly from the pools allowing in situ analysis of radiation removal. One ultraviolet lamp was connected to one of the bioreactors to compare biofilms development and radioactivity retention in both conditions. Stainless steel and titanium balls were introduced into these bioreactors and were removed after different time periods, up to 635 days for stainless steel balls and up to 309 days for titanium. Biofilms developed on the balls were analyzed by scanning electron microscopy and epifluorescence microscopy. DNA was extracted from the biofilms, was cloned and then the microorganisms were identified by independent culture techniques. Biofilms ability to retain radionuclides was also determined by gamma spectrometry. The identified radioresistant organisms belong to the phylogenetic groups Alphaproteobacteria, Gamma-proteobacteria, Actinobacteria, Deinococcus-Thermus and Bacteroidetes. The sequences of these microorganisms have been deposited in GenBank (access numbers KR817260-KR817405). A different distribution of microorganisms was observed in relation to the type of bioreactor. Biofilms have essentially retained activation radionuclides. Sometimes the sum of Co-60 and Mn-54 reached 97%. Cr-51, Co-58, Fe-59, Zn-65 and Zr-95 have also been retained. A radionuclide retention process mechanism related to biofilms formation and disappearance time is suggested. It has been assessed that the scalable process can be economically profitable.

Análisis de Riesgos y recomendaciones para la aplicación de la técnica de Fracturación Hidráulica.

Como es bien sabido, actualmente la Unión Europea pasa por un momento de crisis energética y en concreto España tiene el problema de su gran dependencia energética de otros países, al mismo tiempo que el consumo de gas natural ha aumentado. El hecho de que Estados Unidos se haya autoabastecido energéticamente gracias al gas de esquisto extraído por fracturación hidráulica, hace que en Europa y consecuentemente en España, se vea esta técnica de forma favorable dado que reduciría esta dependencia energética. Pero en contra, la fracturación hidráulica ha dado lugar en Estados Unidos a una serie de accidentes cuyas consecuencias han tenido impacto en el medio ambiente y esto ha provocado un gran debate social a cerca de las consecuencias medioambientales de la extracción de gas de esquisto. Es inevitable, que al igual que ocurre con la extracción de hidrocarburos convencionales, también durante la extracción de hidrocarburos no convencionales existan unos riesgos asociados. Este Trabajo Fin de Máster se engloba dentro de un proyecto titulado “Seguridad en un proyecto de gas de esquisto” concretamente en la convocatoria de 2014 de “Proyectos I+D+I, del Programa Estatal de Investigación y Desarrollo e Innovación orientado a los retos de la sociedad”, por el Ministerio de Economía y Competitividad. Se ha realizado un estudio de incidentes ocurridos en Estados Unidos que es donde más experiencia hay en la práctica de la técnica de la fracturación hidráulica y una vez conocidas sus causas y el tipo de impacto se ha asignado las probabilidades de que ocurran y en qué fase del proyecto de extracción de gas de esquisto es más probable que tengan lugar. Las principales preocupaciones en relación con el medio ambiente que se plantean respecto a la explotación del gas no convencional mediante la técnica de la fracturación hidráulica son: contaminación de los acuíferos (bien sea por el fluido de fracturación empleado o por el metano); consumo de agua necesaria durante la fracturación hidráulica; el tratamiento, control y posible radiactividad de las aguas de retorno; sismicidad inducida por la fracturación hidráulica y contaminación atmosférica por emisiones incontroladas, especialmente metano. Con el objetivo de controlar, prevenir o mitigar los riesgos que conlleva el desarrollo de estas prácticas, la Unión Europea estableció la Recomendación (2014/70/UE), de 22 de enero de 2014. Del mismo modo, debido al incremento del número de permisos de investigación de esta técnica que se solicitaron en España, las actividades de fracturación hidráulica se regularon en el artículo 9.5 de la Ley 17/2013 del Sector de Hidrocarburos y por la Ley 21/2013 de Evaluación de Impacto Ambiental. Consecuentemente se han propuesto en este trabajo una serie de recomendaciones para minimizar los incidentes que pudiesen ocurrir en España provocados por esta técnica. Las más destacadas son: correcto diseño del pozo, realización de estudios hidrológicos y redes de control, reducción del consumo de agua mediante reutilización de las aguas de retorno y de producción después de que se traten en plantas de tratamiento de aguas residuales, utilización de equipos de medida de sustancias químicas en la planta, monitorización sísmica en tiempo real 3D del pozo, así como planes de vigilancia de tipo semáforo para reducir la sismicidad inducida por la aplicación de la técnica.

Obtención de la función receptora en las estaciones sísmicas de banda ancha del Instituto Geográfico Nacional.

La estructura de velocidad en los primeros metros debajo de una estación sísmica produce un cierto efecto sobre el sismograma que ésta registra, dicho efecto se denomina función receptora. Si se conoce esta función, es posible, mediante modelización, obtener esa estructura de velocidad. La onda P que genera un telesismo, queda “atrapada” en los últimos kilómetros de su trayectoria en la estructura próxima a la estación que la registra, reverberando entre las discontinuidades de la misma y convirtiéndose, en algunas partes, en onda S. La diferencia de registro existente entre el sismograma del rayo no atrapado (componente vertical del registro) y el producido por este fenómeno (componente radial), sería la función receptora. En este proyecto, se pretende obtener mediante la técnica de Langston (1979),las funciones receptoras de las 55 estaciones que forman la Red Sísmica de banda ancha del Instituto Geográfico Nacional, con el fin de que en trabajos venideros, sea posible conocer de una manera más realista, la estructura de la corteza terrestre bajo el territorio español. Para la obtención de estas funciones, será necesario trabajar con un conjunto de sismogramas registrados por estaciones de banda ancha a partir de telesismos cuyos parámetros focales reúnan una serie de condiciones, definidas y detalladas ampliamente en esta memoria. Tras un proceso de selección, que será constante a lo largo de todo este trabajo,se van a emplear para el cálculo de las funciones receptoras, un total de 5231 sismogramas, generados por 186 telesismos. Para terminar, a modo de ejemplo y con el objetivo de mostrar la finalidad para la que se ha llevado a cabo este proyecto, se procederá a realizar la inversión de la función receptora obtenida para una de las estaciones sísmicas de la red, obteniendo así, un modelo de velocidades aproximado de la estructura de la litosfera existente bajo dicha estación. Abstract: The speed structure in the first layers under a seismic station produces an effect on the seismogram that is registered. This effect is known as receiver function. It is possible to obtain the speed structure if this function is detected through modelling. The P wave, which generates a teleseism, gets trapped in the final kilometres of its trajectory in the nearby station structure, reverberating among its discontinuities and turning into a S wave. The registered difference between the not caught ray seismogram (vertical component of the register) and the one that is produced by this phenomenon (radial component) would be known as receiver function. This project aims to obtain, through the Langston technique (1979), the receiver function of the 55 stations that make up the Broadband Seismic Network of the National Geography Institute to get to know, in a more realistic way, the structure of the Earth’s crust under the Spanish territory. In order to obtain these functions, it will be necessary to work with a number of seismograms registered by the broadband stations whose focal parameters meet the requirements that are well defined and detailed in the development of this memory. After selecting the appropriate data, that would be constant throughout the development of this project, we are going to apply a total of 5231 seismograms obtained by 186 teleseism to calculate the receivers functions. To conclude, and with the aim of demonstrating the purposes and uses that lie behind the development of this project, we will reverse the receiver function for one of the seismic stations, obtaining a approximate speeds’ model of the lithosphere structure that can be found under the previously mentioned station.