XXIICONGRESONACIONALDEHIDRULICA

AMH AMH
ACAPULCO, GUERRERO, MXICO, MARZO2017

MODELACIN NUMRICA DE LOS SEDIMENTOS ENTRE EL P.H.LAS CRUCES Y LA

PRESA DE CAMBIO DE REGIMEN (PCR)
lvarez Celso Irving Juvenal, Reyes Soto Fernando, Snchez Cruz Fidel, Solano Ibarra Cristian
Miguel ngel
Laboratorio de Hidrulica, Comisin Federal de Electricidad, Cananea 101 Col. Lomas de la Selva C.P. 62270, Cuernavaca
Morelos.
irving.alvarez@cfe.gob.mx, fernando.reyes10@cfe.gob.mx, fidel.sanchez01@cfe.gob.mx, cristian.solano@cfe.gob.mx

Introduccin
Los sedimentos son un aspecto muy importante a tener en 29,098 km, que drenan al Ocano Pacfico un volumen medio
cuenta en la hidrulica fluvial, ya que los parmetros anual superior a los 2,700 hm, que se traduce en alrededor de
caractersticos de un rio no son constantes. El clculo vara 87 m/s de gasto medio anual.
segn el caudal, las caractersticas fisiogrficas de la cuenca,
pendiente del ro, rugosidad y granulometra. Todo esto en
conjunto, hace que los clculos de sedimentos no sean del todo
precisos. Sin embargo, se dispone de modelos numricos para
predecir aproximadamente, en lneas generales, el
comportamiento de un ro. En este trabajo se estudia, mediante
modelos numricos, el comportamiento de los sedimentos
alojados en el tramo del Proyecto Hidroelctrico Las Cruces
(PHC) y la Presa de Cambio de Rgimen (PCR), mediante la
aplicacin de tres software, HEC-RAS, IBER y FlOW3D.

Antecedentes
Figura 1. Localizacin del PHC y la PCR.
El PHC tiene como objetivo contribuir a satisfacer la demanda
de energa elctrica pronosticada para el rea Occidente de
Datos del Proyecto
Mxico. La generacin de energa elctrica en el PHC no ser
continua a lo largo del da, variando el nmero de horas de
Informacin topogrfica
operacin dependiendo de la temporada del ao en la que se
Se utiliz informacin Lidar de curvas de nivel a cada metro en
encuentre el embalse. En esta situacin, la descarga de agua de
el tramo entre el PHC y la PCR, esta informacin fue
la central hidroelctrica variar de igual manera, por ello, con
proporcionada por la Gerencia de Estudios de Ingeniera Civil
la finalidad de reducir de manera significativa estas
de la CFE (Figura 2).
fluctuaciones en el caudal del ro, se ha incorporado al diseo
del PHC una PCR ubicada aguas abajo de la cortina, donde se
considera instalar una mini central generadora. Por lo anterior,
la Subgerencia Tcnica de Proyectos Hidroelctricos de la
Comisin Federal de Electricidad (CFE), instruy al
Laboratorio de Hidrulica desarrollar un proyecto en modelo
numrico, para representar la caracterizacin y determinar el
comportamiento de los sedimentos, en el tramo del PHC y la
PCR.

Objetivo
Determinar la caracterizacin y comportamiento de los
sedimentos, en modelo numrico, en el tramo del PHC y de la
PCR. Analizar la erosin local aguas abajo del tanque
amortiguador de la PCR y determinar si la colocacin de la
obra de toma de la PCR es la adecuada en la margen derecha.
Figura 2. Levantamiento Lidar.
Localizacin
Informacin hidrolgica e hidrolgica
La PCR asociada al PHC se ubica aproximadamente 6 km
La informacin hidrulica e hidrolgica fue proporcionada por
hacia aguas abajo de ste, en la Regin Hidrolgica No. 11
la Gerencia de Tcnica de Proyectos Hidroelctricos de la CFE,
Presidio - San Pedro de la cuenca del ro San Pedro-Mezquital
la cual se muestra a continuacin:
en el estado de Nayarit, Mxico (Figura 1). Contiene una
superficie de
 XXIICONGRESONACIONALDEHIDRULICA
AMH AMH
ACAPULCO, GUERRERO, MXICO, MARZO2017

Tabla 1. Datos hidrolgicos del PHC y la PCR Pre-proceso

Proyecto Hidroelctrico Presa de Cambio de Esta etapa comprende la elaboracin de las geometras,
Las Cruces (PHC) Rgimen (PCR) generacin del mallado y la definicin de condiciones de
Gasto de Gasto de frontera e iniciales para los tres modelos (Figura 5, 6 y 7).
Elevacin Elevacin
diseo diseo
[msnm] [msnm]
[m3/s] [m3/s]
238 360 84.28 530
243 10,930 92.60 11,300

Figura 5. Pre-procesamiento del modelo 1D.

Figura 3. Gastos medios mensuales del PHC.

La granulometra fue extrada de un informe realizado por la

empresa consultora Fluvitecno S.A, de C.V. (FLUVITECNO,
2012), que la Figura 4 muestra la Curva granulomtrica
correspondiente.

Figura 6. Pre-procesamiento del modelo 2D.

Figura 4. Curva granulomtrica en la PHC.

Metodologa

Se utilizaron tres modelos numricos para determinar la

caracterizacin y comportamiento de los sedimentos, estos
dependen del dominio del problema y de los fenmenos locales
que se pueden representar. En un primer modelo
unidimensional o 1D (HEC-RAS) se determinar el manejo de
los sedimentos de arrastre de fondo en los tramos entre el PHC
y la PCR, el cual permita conocer la evolucin de tales
sedimentos. Un segundo modelo numrico bidimensional o 2D
(IBER) determinar la distribucin de los sedimentos de
arrastre de fondo y en suspensin en las cercanas de la
captacin, para determinar la concentracin y el tamao de los
Figura 7. Pre-procesamiento del modelo 3D.
sedimentos que llegarn a la toma de generacin de la PCR.
Finalmente, un tercer modelo numrico tridimensional o 3D Procesamiento
(FLOW3D) permitir estudiar el comportamiento y tamaos de El procesamiento son la(s) simulacin(es) propiamente, en esta
los sedimentos mediante la operacin del vertedor y etapa se debe cuidar y vigilar la convergencia y estabilidad de
desarenador de la PCR, as como la zona de erosin aguas los modelos para conocer el desempeo de los clculos
debajo de su tanque amortiguador. realizados por los programas. Para el modelo 1D se realizaron
Cada modelo consta de tres etapas para su correcta ejecucin, a dos escenarios de simulacin, en el primero tom en cuenta los
continuacin se enlistan dichas etapas: gastos medios mensuales y el segundo escenario ocup la
 XXIICONGRESONACIONALDEHIDRULICA
AMH AMH
ACAPULCO, GUERRERO, MXICO, MARZO2017

avenida de diseo de 11,300 m3/s. Con el modelo 2D se

ejecutaron dos simulaciones, la primera para el anlisis
hidrodinmico en el tramo del PHC y la PCR para la avenida
de diseo y la segunda para un anlisis de sedimentos en el
mismo tramo, tambin para la avenida de diseo. Finalmente,
con el modelo 3D se analiz el comportamiento y se
determin los tamaos de sedimento que pasan a travs del
vertedor y desarenador de la PCR ante la avenida de diseo.
Pos Procesamiento
En la etapa de pos-procesamiento se realiz la extraccin de
datos generados en la etapa de procesamiento. Los resultados Figura 9. Calibracin entre los modelos 2D y 3D (lvarez y Reyes,
obtenidos posteriormente se graficaron y analizaron para 2016) de la PCR
conocer el comportamiento de la estructura en cuestin. Sin
embargo, para tener confiabilidad en los resultados obtenidos La calibracin del modelo tridimensional (Tabla 3) se realiz
se realiz una calibracin de la modelacin numrica. conforme a dos aspectos, el hidrodinmico y el de socavacin. El
Calibracin numrica hidrodinmico se calibr con ayuda de las curvas elevacin-gasto
La calibracin del modelo numrico unidimensional (Tabla 2) de la cortina de la PCR calculadas por la Gerencia Tcnica de
se realiz al hacer la comparativa con las curvas elevacin- Proyectos Hidroelctricos de la CFE.
gasto de las dos presas, obtenidas por la Gerencia Tcnica de
Proyectos Hidroelctricos de la CFE. Tabla 3. Calibracin del modelo 3D. Anlisis hidrodinmico
Estructura en
Tabla 2. Calibracin del modelo 1D Modelo 3D Datos de Diseo
funcionamiento
Modelo 1D Datos de Diseo Elevacin Gasto Elevacin Gasto
Elevacin Gasto Elevacin Gasto [msnm] [m3/s] [msnm] [m3/s]
Presa
[msnm] [m3/s] [msnm] [m3/s] 10.8 Vertedor
92.60 92.60 10,830
238.40 360 238 360 45
P.H.C. 11,3 Vertedor +
243.50 11,300 243 10,930 92.60 92.60 11,300
84.60 550 84.28 530 95 desarenador
P.C.R. 92.60 550 92.60 530 Desarenador
93.20 11,300 92.60 11,300

La calibracin del modelo bidimensional se hizo mediante la Para evaluar la calidad de los resultados de la socavacin total
comparacin con dos modelos tridimensionales (Figura 8) la aguas abajo del tanque amortiguador obtenida con el modelo
correlacin entre los modelos comparados ya que la numrico 3D, se compararon dichos valores con los calculados
distribucin de velocidades fue similar y la correlacin media mediante algunos mtodos empricos.
en la magnitud de esta variable fue del 97% respecto a la
velocidad media correspondiente al modelo3D. El primer Tabla 4. Calibracin del modelo 3D. Anlisis de socavacin
modelo fue el correspondiente al estudio en modelo numrico Socavacin
3D del PHC (lvarez, 2014), que fue corroborado con un Mtodo hs
[m]
modelo fsico del mismo proyecto (Montero, 2015). El
Schoklitsch (1932) -2.28
segundo, con los resultados hidrodinmicos del modelo
Lischtvan-Lebediev (1959) 3.79
numrico 3D del PHC (lvarez y Reyes 2016). Ambos estudios
Laursen (1963) 3.19
fueron elaborados por rea de Modelos Numricos del
Maza lvarez y Garca Flores (1978) 3.20
Laboratorio de Hidrulica de la CFE (Figura 9).
Neill (1980) 4.72
Lacey (1930) 3.33
Modelo 3D 3.95

Anlisis de resultados
Para analizar el comportamiento hidrulico y de sedimentos en el
tramo del PHC y la PCR, as como dela obra de excedencias de
esta ltima, se realizaron diferentes condiciones de simulacin, a
saber:
Anlisis Hidrodinmico
Tirantes

Figura 8. Calibracin entre los modelos 2D y 3D (lvarez, 2014) Los resultados del programa HEC-RAS fueron exportados a un
de la PHC GIS, en los cuales se puede observar en Figura 10 los tirantes para
los dos escenarios realizados con este modelo, del lado izquierdo
se visualizan los tirantes asociados a los gastos medios
mensuales, en el cual se aprecia que la PCR no vierte un gasto
significativo; estos gastos corresponden al nivel del NAMO
(84.28 msnm). Del lado derecho se muestran los tirantes para la
avenida de 11,300 m3/s, donde se aprecia que la zona aguas abajo
 XXIICONGRESONACIONALDEHIDRULICA
AMH AMH
ACAPULCO, GUERRERO, MXICO, MARZO2017

de la cortina de la PCR tiene un flujo apreciable, que m/s, entre el PHC y la PCR es 3.3 m/s, mientras que en el vertedor
corresponde al nivel del NAME (92.60 msnm), en de cresta libre de la PCR es de 10.45 m/s.
consecuencia, tambin se presentan las velocidades ms
grandes de este modelo numrico. En la Figura 12 se muestra la distribucin de velocidades en la
presa y vertedor del PHC, donde los valores mximos se
encuentran en dicho vertedor y, sobre todo, en la descarga de este,
siendo estas del orden de 40 m/s. En el campo de velocidades de
la PCR predominan velocidades bajas antes y despus del
vertedor (de 0 a 5 m/s), y la mayores se presentan en la cada de
dicho vertedor, con valores prximos a los 18 m/s.

Figura 10. Tirantes en el modelo 1D

Del modelo 2D se presentan los resultados obtenidos del
anlisis hidrodinmico para el caudal de diseo (Q=10,900
m3/s). En la Figura 11 se puede observar de manera especfica
el embalse y el vertedor del PHC, con la variacin de los niveles
del agua. El nivel medio en embalse se calcul en 243 msnm,
que corresponde al NAME del proyecto mencionado, al igual
se pueden apreciar los niveles de agua contiguos a la PCR. El
nivel en el embalse calculado fue de 92.6 msnm, que tambin
corresponde al de NAME conforme a los datos de proyecto de
dicha estructura.
Figura 12. Velocidades en el modelo 2D

La Figura 13 muestra las velocidades del modelo numrico 3D

para un gasto de 11,300 m3/s, donde se observa una zona muerta
en las mrgenes de la PCR y aguas abajo del tanque
amortiguador, esto ltimo se debe a la disminucin de la
velocidad inducida por el tanque amortiguador. En el tanque
amortiguador se observan velocidades que oscilan en 16.00 m/s
en la zona del cimacio y de 1.00 m/s antes de salir de la estructura.
La velocidad mxima se encuentra antes del choque del flujo en
el tanque amortiguador (pie del cimacio), con un valor de 18.00
m/s; en tanto que a la salida del tanque se encuentran velocidades
en torno a los 8.00 m/s y no se aprecia un recargamiento en las
mrgenes del ro. Se realizaron dos cortes para visualizar en perfil
las velocidades en el cimacio y en el desarenador de la PCR.

Figura 11. Tirantes en el modelo 2D

Velocidades
Las velocidades calculadas con el programa HEC-RAS varan
segn sea el escenario simulado. Para el escenario uno con los
gastos medios mensuales, las velocidades aguas arriba de la
cortina de la PHC son del orden de 0.10 m/s, entre el tramo entre Figura 13. Campo de velocidades en el modelo 3D
el PHC y la PCR son de 0.15 m/s y a la salida de sta ltima son
de 1.2 m/s. Para el escenario dos con la avenida de diseo, la En la Figura 14 se presenta el corte uno, donde se observa el
velocidad media a la entrada del vaso del PCH es de 0.10 m/s, campo de velocidades en la rpida del vertedor de excedencias y
en la descarga del vertedor de excedencias del PCH es de 25.35 la zona inferior del tanque amortiguador, donde se presentan las
velocidades mximas que alcanzan los 18.00 m/s; se aprecia la
 XXIICONGRESONACIONALDEHIDRULICA
AMH AMH
ACAPULCO, GUERRERO, MXICO, MARZO2017

formacin de un salto hidrulico incipiente al pie del cimacio, finos van en suspensin. En la Figura17 se muestra un
as como el cambio de rgimen de supercrtico a sub-crtico, acercamiento (zoom) de la zona de la PCR, donde se aprecia
provocando al final una disminucin en la magnitud de las una socavacin de 2.0m de profundidad aguas abajo del tanque
velocidades de unos 8.00 m/s. amortiguador del vertedor de cresta libre, y una pequea
acumulacin de material de unos 0.07 m a la entrada de la
cortina.

Figura 14. Campo de velocidades en el vertedor de cresta libre de

la PCR

En la Figura 15 se aprecia el corte dos, donde se muestra el

campo de velocidades en la zona de la compuerta del
desarenador, donde las velocidades a la entrada de la compuerta
son de 8.00 m/s, aproximadamente y las mximas del orden de
17.00 m/s a la salida de la contraccin. Figura 16. Distribucin de sedimento en el tramo PHC-PCR

Figura 15. Campo de velocidades en el desarenador de la PCR

Anlisis de sedimentos
Con el modelo unidimensional se calcul, para el escenario uno,
una cantidad de sedimentos en el tramo del PHC PCR de Figura 17. Erosin y sedimentacin en la PCR.
5,684.90 m3/s para un ao; esto se debe a la cortina del PHC, la
cual acumula los sedimentos gruesos en la parte inicial de su A continuacin se presentan los resultados referentes al estudio
almacenamiento. El mtodo de clculo para este escenario fue del anlisis de la evolucin de los sedimentos entre el PHC y la
el de Engelnd Hasen debido a las caractersticas PCR del modelo bidimensional
morfolgicas de la cuenca, la granulometra del lecho y la
pendiente promedio del ro en estudio. En la Tabla 6 se En la Figura 18 se presenta una vista general de la erosin en la
muestran los volmenes de sedimentos estimados con cada una zona de estudio, en la que se observa que el mayor movimiento
de las siete frmulas de transporte que tiene el software de material se presenta en la descarga del vertedor del PHC, con
HEC.RAS, incluyendo la de Engelnd Hasen. erosiones del orden de 60 m que forman un cuenco de
socavacin y la posterior deposicin de sedimento aguas abajo
Se utiliz la frmula de Engelnd Hasen para el clculo de de dicho cuenco con alturas de 20 m; mientras que aguas arriba
arrastre de sedimentos para el gasto de 11,300 m3/s, de la PCR no se presenta gran acumulacin de sedimento, cabe
obteniendo un valor de 73,269.00 m3 en 100 horas de sealar que esto ocurre slo despus del paso de la avenida de
acumulacin de sedimento a lo largo del tramo entre la presa diseo (48 horas). En consecuencia, en esta simulacin no se
del PHC y la PCR; este valor es 13 veces mayor, al estimado alcanza a ver acumulacin de material en alguna de las
con los gastos medios mensuales. En la Figura 16 se muestra mrgenes cercano a la obra de toma, fenmeno que pudiera
la configuracin final de la plantilla antes y despus del paso observarse al cabo de varios aos. Por ltimo, aguas abajo del
de la avenida. La mayor erosin o remocin en este caso, se vertedor de cresta libre de la PCR se observa una erosin sobre
encuentra en la seccin transversal del cadenamiento 0+4800 la margen derecha alrededor de dos metros cercana a la frontera
y es del orden de 15 metros, esto se debe a la descarga del de clculo.
vertedor de demasas del PHC. Para el caso de la
sedimentacin o acumulacin de material, esta se localiza en
la seccin transversal del cadenamiento 0+1350 con un valor
de 2.10 m, debido a la configuracin del ro, en especial a las
curvaturas y la disminucin de la velocidad antes de la llegada
a la cortina de la PCR, ya que los materiales ms grandes solo
son transportados por arrastre de fondo, mientras que los ms
 XXIICONGRESONACIONALDEHIDRULICA
AMH AMH
ACAPULCO, GUERRERO, MXICO, MARZO2017

Figura 21. Lecho del ro antes y despus de la avenida de diseo

Las dimensiones de la erosin en sentido del flujo dominante

(X) es de 175 metros aproximadamente, en Y igual a 200
metros, equivalente al ancho de la cortina de la PCR y
finalmente, la profundidad es variable, siendo mayor la erosin
en la cercana al Tanque amortiguador y disminuye
gradualmente conforme se aleja de esta estructura.
En la Figura 22 se muestra un corte (plano X-Z) en margen
Figura 18. Erosin y sedimentacin en el tramo del PCH y la izquierda, donde aparecen zonas de erosin, en las cuales tienen
PCR. Modelo Bidimensional una profundidad mxima de 3.95 m en la zona inmediatamente
aguas abajo del tanque amortiguador.
El modelo tridimensional tiene la capacidad de representar de
forma individual cada dimetro de sedimento as como la
concentracin que se presente. En la Figura 19 se aprecia que
el material de 0.004 m pasa con dificultad y en pocas
cantidades sobre la cresta del vertedor, es removido
completamente en el tanque amortiguador y distribuido hacia
aguas abajo.

Figura 19 vertedor de excedencias Figura 22. Lecho del ro antes y despus de la avenida de diseo

En la Figura 20 se aprecia que los dimetros mayores de 0.008 Conclusiones

m difcilmente ocurre su desalojo por la compuerta del
desarenador, a pesar de ocurrir un gasto de 11,300 m3/s y el Modelo unidimensional (HEC-RAS)
NAME = 92.60 msnm. As mismo, se aprecian
Desde el punto de vista hidrulico y de comportamiento de
concentraciones pequeas del orden de 0.35 Kg/m3 al inicio
sedimentos, en el presente estudio se plantearon dos escenarios
del desarenador.
de simulacin:
En el escenario 1, anlisis hidrodinmico y de sedimentos para
gastos medios mensuales, se muestra un comportamiento lento
del flujo (sub-crtico) en el tramo del PHC y la PCR, esto se
debe al vaso de sta ltima. El tener un flujo sub-critico, hace
que el transporte de sedimentos gruesos (tamaos de 0.1 a 1
mm) sea casi nulo para los gastos medios mensuales, llegando
as a tener un volumen de 5,685 m3/ao (Mtodo de Engelnd-
Figura 20. Concentraciones de sedimento D50 en la compuerta del Hasen); con este valor anual se estima que en 10 aos el nivel
desarenador de sedimento en el vaso de la PCR tendr una elevacin de
apenas 0.1 metros. En cuanto a la socavacin, no se encontraron
Erosin local erosiones significativas.
En la Figura 21 se muestra la configuracin del lecho del ro En el escenario 2, anlisis hidrodinmico y de sedimentos para
antes y despus de la avenida de 11,300 m3/s, considerando el gasto de 11,300 m3/s, se presenta un comportamiento rpido
rgida toda la cortina y las bermas de la PCR y fondo mvil el del flujo (rgimen sper-crtico) en el tramo de estudio, pues el
lecho del ro. rango de las velocidades calculado entre 3.3 y 26.0 m/s, da
como consecuencia, que el promedio del nmero de Froude en
todas las secciones sea mayor a la unidad, por lo que el arrastre
de fondo es mayor, con un volumen acumulado de 73,269.00
 XXIICONGRESONACIONALDEHIDRULICA
AMH AMH
ACAPULCO, GUERRERO, MXICO, MARZO2017

m3 (Mtodo de Engelnd-Hasen) a lo largo del tramo entre la para trasladar este dimetro hacia aguas abajo, ms an, se
presa del PH Cruces y la PCR, que es aproximadamente 13 encontr que pueden pasar dimetros mayores dependiendo de
veces mayor al sedimento estimado con gastos medios la poltica de operacin dela obra. Para el caso del desarenador,
mensuales. En esta ltima avenida, como se espera para este tiene la capacidad para desalojar dimetros hasta de 8.00
pequeas avenidas que son las ms frecuentes, se tiene una mm, mientras que el vertedor de cresta libre puede desalojar
acumulacin pequea de material a la entrada de la cortina de dimetros hasta de 4.00 mm. Tambin existe una erosin
la PCR del orden de 0.07 metros, aproximadamente. En significativa aguas abajo del tanque amortiguador, del orden de
cuanto a la socavacin, la ms importante se localiza a la 4.0 metros, aproximadamente, por lo que se recomienda
salida del vertedor de excedencias de la PCH con una aumentar la longutud de dicha estructura para evitar este
magnitud de 14.50 metros aproximadamente, para la avenida fenmeno.
de diseo. Para la PCR, la erosin aguas abajo del tanque
amortiguador se estim en 2.20 metros y es congruente con lo Referencias
obtenido empricamente.
Blad, E., Cea, L., Corestein, G., Escolano, E., Puertas, J.,
Modelo bidimensional (IBER 2.0) VzquezCendn, M.E., Dolz, J., Coll, A. (2014). "Iber:
Como resultado del procedimiento de calibracin del modelo herramienta de simulacin numrica del flujo en ros. Revista
2D respecto a los modelos numrico 3D y fsico, se determina Internacional de Mtodos Numricos para Clculo y Diseo en
que el modelo en cuestin representa adecuadamente el Ingeniera", Vol.30 (1) pp.110. Espaa.
comportamiento hidrodinmico de la zona bajo estudio. Al Brethour J. (2003) Modeling Sediment Scour, Technical note
analizar los resultados obtenidos con los modelos fsicos y FSI-03-TN-62, Flow Science.
numrico 3D, en general los resultados del modelo 2D
permiten tener una distribucin espacial aceptable de los Calzas Prez, Mario (2013). Modelizacin del canal de la
diferentes parmetros, tales como los campos de velocidades planta TEL en la U.P.T. de As Pontes (Tesis de
y la evolucin del comportamiento (socavacin y maestra).Universidad de Corua.
sedimentacin) de los sedimentos. Bajo estas afirmaciones se Caracterizacin geomecnica y determinacin de mdulos de
puede dar por cumplido el objetivo de determinar el anlisis deformabilidad por medio de ensaye de gato Goodman en la
de comportamiento de los sedimentos entre las presas P.H. presa cambio de rgimen alternativa II de P.H. Las Cruces,
Las Cruces y la P.C.R, donde se observ que grandes Nayarit. GEIC, CFE. Junio 2015. (GEIC, 2015)
afectaciones y movimientos ocurren durante el paso de
Cea, L., Puertas, J., Vzquez-Cendn, M.E. (2007). Depth
avenidas grande y al pie del vertedor de excedencias con
averaged modelling of turbulent shallow water flow with wet-
grandes erosiones en el lecho del rio. Las acumulaciones de
dry fronts. Archives of Computational Methods in Engineering,
sedimento en el vaso de la P.C.R durante las avenidas son
State of the art reviews, Vol.14 (3) pp.303-341.Espaa.
mnimas y no mayores a 1 m a 150 m aguas arriba dela cortina.
Estas acumulaciones no se ven afectadas donde se encuentra Cea Gmez, Luis (2008). Modelizacin de la turbulencia en
la obra de toma, ya que en un periodo de avenidas los aguas poco profundas. Boletn de la Sociedad Espaola de
sedimentos no alcanzan a llegar a esta estructura. Cabe Matemtica Aplicada, Num.42.
mencionar que en este trabajo solo se realiz para la avenida Cea, L., Vzquez-Cendn, M. E., Puertas, J. (2009). El
de diseo. Por lo anterior se recomienda hacer un anlisis de mtodo de volmenes finitos aplicado a problemas de
una simulacin con una duracin ms larga para apreciar en ingeniera fluvial y costera. RSME, Vol. 12. Espaa.
cuanto tiempo llegan los sedimentos a la P.C.R. As como
probar el impacto que diferentes gastos tienen en las CIMNE-UdC-UPC. Modelo IBER, Tema 7 Turbulencia [en
estructuras consideradas lnea]. [Fecha de consulta: 23 octubre 2014]. Disponible en:
www.iberaula.com
Modelo tridimensional (FLOW3D)
CIMNE-UdC-UPC. Modelo IBER, Tema 5 Esquemas
El resultado hidrodinmico de la modelacin 3D muestra una
numricos y clculo [en lnea]. [Fecha de consulta: 23 octubre
buena correlacin entre los niveles de operacin y los gastos
de las estructuras. Este primer anlisis asegur que el 2014]. Disponible en: www.iberaula.com
funcionamiento de las obras es bueno y sirvi como base para CEDEX, GEAMA, FLUMEN y CIMNE 2010: IBER.
la simulacin del transporte de sedimentos. Modelizacin bidimensional del flujo en lmina libre en aguas
En las simulaciones del trasporte de sedimentos se pudo poco profundas. Manual bsico de usuario.
conocer la distribucin de los mismos para la avenida de Chanson, H. (2004) The hydraulics of open channel flow: an
11,300 m3/s. Se observa que la sedimentacin no se concentra introduction. Elsevier
en alguna de las mrgenes y no es apreciable para el caso de
Cuervo, Astrid E. (2012). Comparacin de modelos
la avenida de 11,300 m3/s, por lo que es factible construir la
hidrulicos unidimensional (HEC-RAS) y bidimensional
obra de toma en cualquiera de las mrgenes, aunque del
(IBER) en el anlisis del rompimiento de presas de materiales
anlisis hidrodinmico se observ que el flujo entrante se
sueltos. Tesis de maestra. Universitat Politcnica de Catalunya.
estanca sobre la margen derecha del ro.
Espaa.
Con la simulacin numrica se pudo calcular la distribucin
de la concentracin de sedimentos para doce dimetros ESRI. ArcGISResources. Help Library [en lnea].
Actualizada: 07 julio 2014. [Fecha de consulta: 04 enero
distintos de partculas. Un requerimiento importante de la
2015]. Disponible en:
PCR es que los dimetros de 1.00 mm pasen por dicha
estructura hacia la zona de Marismas Nacionales. As las Estudio complementario para la estimacin del transporte y
cosas, se observ favorablemente que el paso de la avenida de depsito de sedimentos desde la presa hidroelctrica Las Cruces
11,300 m3/s, y otras avenidas grandes, no tienen problema Nayarit, en la zona aguas abajo en el cauce del ro San Pedro.
 XXIICONGRESONACIONALDEHIDRULICA
AMH AMH
ACAPULCO, GUERRERO, MXICO, MARZO2017

Informe final. Instituto de Ingeniera, UNAM. Diciembre

2015. (Fuentes, 2015)
Estudios complementarios en materia de polticas de
operacin, diagnstico de zonas costeras, transporte de
sedimentos y manejo de cuenca, del proyecto hidroelctrico
Las Cruces, Nayarit (2da etapa). Informe de modelos fsicos.
Instituto de Ingeniera, UNAM. Abril2016. (Gracia, 2016)
FLUMEN, CEDEX, GEAMA, CIMNE (2014). Manual de
referencia hidrulico IBER. Modelizacin bidimensional del
flujo en lmina libre en aguas poco profundas. Espaa.
FLUMEN, CEDEX, GEAMA, CIMNE (2012). Manual
bsico de usuario IBER. Modelizacin bidimensional del
flujo en lmina libre en aguas poco profundas. Espaa.
FLUMEN. Tutorial Iber [en lnea]. [Fecha de consulta: 25
noviembre 2014]. Disponible en:
Flow Sciences Incorporated. (2011) FLOW-3D Users
Manual Version 10.0, Santa Fe, New Mxico.
Flumen, et al. (2014). Manual de referencia hidrulico IBER.
Modelizacin bidimensional del flujo en lmina libre en aguas
poco profundas.
Flumen, et al. (2012). Manual bsico de usuario IBER.
Modelizacin bidimensional del flujo en lmina libre en aguas
poco profundas.
Garca Flores, Manuel, Maza lvarez Jos A (1994).
Captulo 7, Origen y propiedades de los sedimentos. Manual
de Ingeniera de Ros. Instituto de Ingeniera (UNAM);
Comisin Nacional del Agua (CNA). Mxico.
Garca Flores, Manuel, Maza lvarez Jos A (1994).
Captulo 10, Transporte de sedimentos. Manual de Ingeniera
de Ros. Instituto de Ingeniera (UNAM); Comisin Nacional
del Agua (CNA). Mxico.
Garca Flores, M. y Maza A., Jos A. (1996). Morfologa de
ros. Captulo 11 del Manual de Ingeniera de Ros. Instituto
de Ingeniera, UNAM. CONAGUA. Mxico.
Garca Flores, M. y Maza A., Jos A. (1996). Erosin en
ros. Captulo 13 del Manual de Ingeniera de Ros. Instituto
de Ingeniera, UNAM. CONAGUA. Mxico.
Gracia Snchez, Jess (1994). Captulo 18, Sedimentacin
en embalses. Manual de Ingeniera de Ros. Instituto de
Ingeniera (UNAM); Comisin Nacional del Agua (CNA).
Mxico.
Laboratorio de Hidrulica de CFE (1992,1995 y 2004),
Informe de estudios. (CFE, 2004)
May, R. W. P., Ackers, J. C., & Kirby, A. M.
(2002).Manual on scour at bridges and other hydraulic
structures. (Vol. 551). Ciria. Londres, Inglaterra.
Santo, Mara, Lecumberry, Graciela (2005). El proceso de
medicin. Anlisis y comunicacin de datos experimentales.
Universidad Nacional de Ro Cuarto. Facultad de Ciencias
Exactas. Fsico-Qumicas y Naturales. Departamento de
Fsica. Argentina.
Surez, Jaime (2001). Control de erosin en zonas tropicales.
Instituto de erosin. Bucaramanga, Colombia.
U.S. Army Corps of Engineers. Hydrologic engineering
center (2010). HEC-RAS, River Analysis System. Hydraulic
Reference Manual. Version 4.1. Davis, California