Modelamineto Numerico de Una Onda de Avenida Generada Por L A Rotura de Un Presa de Tierra y Su Estimacion de Riesgo
Modelamineto Numerico de Una Onda de Avenida Generada Por L A Rotura de Un Presa de Tierra y Su Estimacion de Riesgo
Modelamineto Numerico de Una Onda de Avenida Generada Por L A Rotura de Un Presa de Tierra y Su Estimacion de Riesgo
LIMA PERU
2015
1.5 Limitaciones
rea geogrfica: Para realizar el modelamiento numrico de rotura de presa, se
escogi como elemento de estudio a la Presa Yuracmayo ubicada en el distrito de San
Mateo, provincia de Huarochir, departamento de Lima, y para simular la inundacin
provocada por la propagacin de la onda de rotura, se consider las caractersticas
topogrficas del embalse de Yuracmayo y valle aguas abajo tal como se muestra en
la Figura.
Ubicacin:
El embalse de Yuracmayo se encuentra ubicado en la naciente del rio Blanco, que es
tributario por la margen izquierda del rio Rmac a unos 120 Km. al este de Lima y a
4315 m.s.n.m. La presa de Yuracmayo se encuentra en la latitud -11.84 (374224.83
E) y longitud -76.15 (8691076.09 S), a unos 500 m. aguas arriba de la poblacin de
Yuracmayo, distrito de San Mateo, provincia de Huarochir, departamento de Lima.
Extensin:
El rea de estudio a considerarse en el modelamiento hidrulico est delimitado desde
laguna y presa Yuracmayo hasta la confluencia del rio Rmac y el rio Blanco. Tal
como se muestra en la Figura.
2.- Objetivos
2.1 objetivos generales
Simular el comportamiento de una onda de avenida generada por la rotura de una
presa de tierra para realizar su estimacin de riesgo en la presa Yuracmayo.
2.2 objetivos especficos
Analizar el modelo HEC RAS para realizar la simulacin de una rotura de presa y la
onda de avenida generada.
Modelar distintas opciones de colapso de una presa de tierra y analizar los resultados.
Analizar y comparar distintas simulaciones, generados por la variacin en las
condiciones de contorno y frontera.
Describir y analizar las caractersticas hidrulicas de la inundacin valle aguas abajo
generadas por la onda de avenida.
Determinar un plano de riesgo como resultado del modelamiento de la onda de
avenida y realizar la estimacin del riesgo de inundacin por rotura de la presa
Yuracmayo.
d) Elementos finitos
El mtodo de los elementos finitos tambin se ha utilizado para la
resolucin de las ecuaciones de Saint Venant unidimensionales.
II. Esquemas de resolucin de las ecuaciones completas de SaintVenant simplificadas
a) Mtodos hidrolgicos
El mtodo no toma en cuenta la ecuacin de movimiento, siendo el
mtodo de Muskingum el ms destacado.
b) Mtodo de la onda cinemtica
El mtodo, solo considera la friccin y pendiente de la solera en la
ecuacin del movimiento.
c) Mtodo de la onda difusiva
Adems de la friccin y pendiente de solera, el mtodo considera la
presin en la ecuacin de movimiento.
d) Mtodo de la onda dinmica cuasi-permanente
EL mtodo no considera la aceleracin local. Este mtodo ha ido
perdiendo uso, debido al avance de la informtica que hace que los
procesos de clculo sean ms rpidos.
3.2.3 Modelos Cuasi-Bidimensionales
Los esquemas bidimensionales sirven para describir muchos fenmenos
naturales como pueden ser la inundacin de una gran llanura, la confluencia
de dos cauces, el cruce de dos corrientes de agua, el flujo en un cauce ancho e
irregular, etc.
Modelos que incorporan la aproximacin cuasi-bidimensional para llanuras
de inundacin son MIKE11, HEC-RAS y GISPLANA.
3.2.4 Modelos Bidimensionales
Igual que para el caso unidimensional, para la resolucin de las ecuaciones de
Navier Stokes en los ejes x e y, se ha utilizado los mtodos de diferencias
finitas y los elementos finitos, pero para el caso bidimensional tambin, se ha
utilizado la tcnica de discretizacin de volmenes finitos.
Los modelos comerciales ms populares en este campo son PHOENICS,
FLUENT, FLOW3D y STAR-CD, que utilizan el mtodo de volmenes
finitos.
3.3 Hiptesis
Con los resultados obtenidos de la aplicacin del modelo numrico a un escenario
de rotura de presa, se puede evaluar y disminuir las prdidas econmicas, y sobre
todo, la prdida de vidas humanas.
permanente
Donde:
Z1 y Z2 representan las cotas de las secciones respecto a un plano de
referencia arbitrario. En caso de geometras irregulares, como en general
sucede en un ro, se toma la cota del punto ms bajo de la seccin.
Y1 y Y2 son los calados en cada una de las secciones consideradas. En caso
de geometras irregulares, se considera la profundidad respecto del punto ms
bajo de la seccin.
V1 y V2 son las velocidades medias en cada seccin. En el caso de rgimen
permanente se obtiene como el cociente entre el caudal circulante y el rea de
la seccin.
1 y 2 son los coeficientes de Coriolis estimados en cada seccin, lo cual
permite corregir el hecho de que la distribucin de velocidad en la seccin se
aleja de una distribucin uniforme.
H es el trmino que estima la energa por unidad de peso que se disipa entre
las secciones 1 y 2.
4.1.3 Conservacin de la cantidad de movimiento o momentum
En el caso de que la superficie del agua pase por el rgimen crtico, el flujo
deja de ser gradualmente variado y pasa a ser rpidamente variado. En dicha
situacin la ecuacin de la energa deja de ser aplicable. Dichos casos
Donde:
P1 y P2 son las componentes resultantes de la distribucin de presiones
hidrostticas, en la direccin del movimiento, en las secciones 1 y 2
respectivamente.
W es el peso contenido en un volumen de control en la direccin del
movimiento.
F es la fuerza debida a la friccin sobre el contorno entre las secciones 1 y 2.
Q es el caudal circulante
1Y2 son los coeficientes de momentum o coeficientes de Boussines para
las secciones 1 y 2 respectivamente.
es la densidad del agua.
V1 y V2 son las velocidades medias presentes entre las secciones
1 y 2 respectivamente.
5.- Metodologa
5.1 Mtodo empleado
Para realizar el modelamiento numrico de una rotura de presa se requiere algunos
elementos previos como son las caractersticas geomtricas de la presa, una topografa
a escala adecuada de la zona de estudio, valores de rugosidad del valle aguas abajo,
las condiciones iniciales y de contorno del rio en estudio, adems, se necesitar de
una extensin del software ArcGIS, llamado HEC-GeoRAS que ha sido desarrollado
por el Hidrologic Engineering Center (HEC) con la colaboracin del Instituto que
desarrollo el ArcGIS (ESRI).
5.2 Muestra seleccionada y caractersticas
Modelo digital del terreno en formato TIN.
La creacin del modelo digital del terreno MDT en formato TIN (Triangular Irregular
Network o red de tringulos irregulares) de la zona de estudio se hace con las curvas
nivel adecuadas (ver el tem 4.5) Este MDT ser una representacin digital
aproximada de la zona de estudio, tal como se muestra en la figura
Delimitacin de las zonas lmites de inundacin Flowpaths Los flowpaths son las
lneas que definen el flujo de agua, tanto del margen izquierdo, como del derecho, y
del centro del ro.
Los flowpaths son importantes pues delimitan la zona del cauce principal y la llanura
de inundacin.
Creacin de secciones transversales XS CutLines
Se definen cada una de las secciones transversales, dibujando lneas de izquierda a
derecha a travs del cauce principal, dentro de los lmites de la llanura de inundacin.
Cada seccin transversal deber ser perpendicular al cauce y estas lneas no se
debern cruzar. Es importante que la distancia entre secciones sea muy estrecha para
garantizar un buen detalle del estudio.
contorno, siempre producir resultados. Por este motivo, una vez ejecutada la
primera simulacin, es necesario realizar lo que se llama el ajuste del
modelo; es decir, la verificacin y contraste de los datos obtenidos.
5.4 instrumentos de la investigacin
Modelo HEC RAS
El programa HEC-RAS es un modelo hidrulico unidimensional, desarrollado por el
centro de ingeniera Hidrolgica (Hydrologic Engineering Center (HEC)) del cuerpo
de ingenieros de la armada de los EE.UU (US Army Corps of Engineers). Es una
aplicacin que permite la modelacin hidrulica en rgimen permanente y no
permanente en lmina libre de cauces abiertos, ros y canales artificiales.
Para rgimen no permanente, HEC-RAS resuelve las ecuaciones completas de Saint
Venant mediante el mtodo de diferencias finitas implcitas utilizando el esquema de
Preissmann. Por otro lado HEC-RAS ha desarrollado una interfaz grfica de usuario
que facilita las labores de Pre-proceso y Pos-proceso, gracias a esto la edicin y
visualizacin datos se realiza de manera fcil.
Se muestra el esquema del proceso de clculo que sigue HEC-RAS.
Fuente: E.Blad M. Snchez Juny, Modelacin numrica en ros en rgimen permanente y variable. Pg. 180 .
En la Figura se observa que las viviendas del poblado de Caruya estn cerca al rio a
distancias de 5 a 10 metros, lo cual representa un riesgo alto, debido a ocurrencias de
posibles inundaciones que se podran generar por fuertes lluvias, deslizamientos, o
peor an, la rotura de la presa Yuracmayo.
VII Bibliografa
Agencia Catalana de Agua. (2003). lAigua, Delimitaci de zones inundables per
a la redacci de lINUNCAT. Departamento de medio Ambiente. Documentos
tcnicos.
Bautista, Carlos, F. (2007). Deduccin de la ecuacin dinmica del flujo
gradualmente variado a partir de las ecuaciones de Saint Venant.. Revista digital
de la Universidad Autnoma de Zacatecas, Volumen 3.
Blad, E. Snchez, J. Snchez H.P. Nierola, D. Gmez, M (2009),
Modelacin numrica en ros en rgimen permanente y variable. Una visin a
partir del modelo HEC-RAS. Edicion UPC.
Cea, Luis. Blade, Ernest. Modelizacin matemtica en lecho fijo del flujo en ros
Modelos 1D y 2D en rgimen permanente y variable. Grupo
GIAMA. Grupo Flumen. E.T.S. de Ing. De Caminos.
Fox, Robert. McDonald, T. Alan. (1995). Introduccin a la mecnica de los
fluidos. Cuarta edicin. McGraw-Hill. Mxico. Pgs. 19 43.
HEC-RAS (2010). River analysis system. Gua de Aplicacin. U.S. Army
Corps of Engineer. Version 4.1.
HEC-RAS (2010). River analysis system. Gua de Referencias. U.S. Army
Corps of Engineer. Version 4.1.
HEC-RAS (2010). River analysis system. Notas de Realizacin. U.S. Army
Corps of Engineer. Versin 4.1.
HEC-RAS (2010) River analysis system. Manual de Usuario. U.S. Army
Corps of Engineer. Versin 4.1.
HEC-GeoRAS (2011) GIS Tools for Support of HEC-RAS using ArcGIS
ICOLD. (1998). Dam Break Flood Analysis, Review and recommendations.
Boletn 111. Pars.
Ministerio de Medio Ambiente de Espaa. (1996). Clasificacin de presas en
funcin del riesgo potencial. Direccin General de Obras
Hidrulicas y Calidad de las Aguas. Gua tcnica
VII Anexos
HIDROGRAMAS OBTENIDOS POR EL MODELAMIENTO HIDRULICO
PARA LOS ESCENARIOS DE SOBREVERTIMIENTO Y TURIFICACIN
Hidrogramas
Los hidrogramas siguientes son los obtenidos luego del modelamiento numrico, estos
pertenecen a las secciones transversales ms importantes.
Primer escenario: Por sobrevertimiento
El la Figura A-1 se observa el hidrograma para la seccin 19433 (ubicacin de la presa),
donde se observa la variacin en funcin del tiempo del caudal (Flow), el nivel del agua
(Stage HW) y el nivel del agua en la seccin siguiente aguas abajo (Stage TW).
El la Figura A-2 se observa el hidrograma para la primera seccin ubicada aguas debajo de la
presa (seccin 19432) donde la variacin est en funcin del tiempo, del caudal (Flow) y el
nivel del agua (Stage)
El la Figura A-3 se observa el hidrograma para la ltima seccin ubicada aguas abajo
(seccin 27) donde la variacin est en funcin del tiempo, del caudal (Flow) y el nivel
del agua (Stage)
El la Figura A-2 se observa el hidrograma para la primera seccin ubicada aguas debajo de la
presa (seccin 19432) donde la variacin est en funcin del tiempo, del caudal (Flow) y el
nivel del agua (Stage)
El la Figura A-6 se observa el hidrograma para la ltima seccin ubicada aguas abajo (seccin
27) donde la variacin est en funcin del tiempo, del caudal (Flow) y el nivel del agua (Stage)