Cuenca Hidrografica
Cuenca Hidrografica
Cuenca Hidrografica
FACULTAD DE INGENIERÍA
TESIS
PRESENTA
JOSÉ LUIS ARAGÓN HERNÁNDEZ
DIRECTOR DE TESIS
DR. ÓSCAR ARTURO FUENTES MARILES
JURADO ASIGNADO
DIRECTOR DE TESIS
FIRMA
Agradecimientos
AGRADECIMIENTOS
Con todo cariño y respeto para mis padres Gonzalo (†) y Francisca, quienes me brindaron su
apoyo, me enseñaron a ser perseverante en la vida y luchar por mis ideales. Para ti Papá
porque en tu vida diste un ejemplo de integridad, de cariño con tu familia y sobre todo de valor
ante la vida.
A mis hermanos Epifanio, Isabel, Bricia, Alejandra, Felicítos, Federico, Francisco, Vicente,
Justino, Félix, Cristina, Edilberto, cuñados y sobrinos, que en algún momento de mi vida me
han ayudado a ser mejor, y quisiera compartir con ustedes este logro.
A mi futura esposa Flor por permitirme disfrutar de momentos muy felices a su lado y
comprenderme en momentos difíciles.
A todos mis amigos y compañeros del Instituto de Ingeniería, UNAM, Faustino de Luna y Juan
Ansberto Cruz, por su colaboración en la revisión y comentarios para la mejora de este trabajo,
Juan Javier Carrillo Enedino Mendoza, Martín Jiménez, Adriana Rodríguez, Edith Vega,
Abraham Ávila, Alicia Ravelo, Juan Gabriel López, José Luis Herrera, por su amistad y apoyo
desinteresado.
A mis amigos de Zoquitlán Octavio, Ramiro, Justino, Donaciano, Naum, Hugo, Daniel, Rosalba,
Simón, de Río Seco Isaías, Javier, Régulo, Naum, Gil, por compartir con ellos momentos
inolvidables.
A mis profesores y sinodales Dr. Óscar Fuentes, M. I. Víctor Franco, Dr. Ramón Domínguez, Dr.
Jesús Gracia, Dr. Carlos Escalante, por su paciencia en la formación de nuevos estudiosos de
este bello campo del conocimiento, la Hidráulica.
A todas aquellas personas que no menciono pero que siempre han estado conmigo y que
saben que siempre las llevo en mi corazón.
RESUMEN
Para que se presenten los escurrimientos máximos en un sitio, en respuesta a una
precipitación, deben transcurrir varios minutos u horas después de iniciada la tormenta (tiempo
de retraso), por lo que es posible anticipar la ocurrencia de los mayores flujos de agua para
advertir el peligro por desbordamientos de los ríos o por la intensidad de sus corrientes. Esta
idea es la base teórica de un Sistema de Alerta Hidrometeorológica.
i
Resumen
ii
Índice
ÍNDICE
1 INTRODUCCIÓN…………………………………………………………………………………...…..1
1.1 Justificación………………………………………………………………………………………… 3
1.2 Objetivos………………………………………………………………………………………….… 4
1.3 Contenido…………………………………………………………………………………….…….. 4
4 PROCESAMIENTO HIDRÁULICO……………………………………………………………….... 33
6 CASO DE APLICACIÓN………………………...………………………………….…......….…….. 53
7 CONCLUSIONES…………………………………………………………………………….……..... 93
8 BIBLIOGRAFÍA Y REFERENCIAS………………………………………………………….….….. 97
ANEXO A………………………………………………………………………………………………… 99
iv
Índice de figuras
ÍNDICE DE FIGURAS
1 INTRODUCCIÓN
4 PROCESAMIENTO HIDRÁULICO
v
Índice de figuras
6 CASO DE APLICACIÓN
6.1. Inundación del 04 al 06 de octubre de 2003. El río Sabinal se desbordó cerca del
estrechamiento del puente ubicado en las calles 3ª Ote. Y 7ª. Nte.……….……………………… 53
6.2 Ubicación geográfica de la cuenca del río Sabinal………………….………….……..……… 54
6.3 Manchas urbanas más importantes………………………………………...……………….…. 55
6.4 Localidades ubicadas dentro de la cuenca en estudio……………..……….……………..… 55
6.5 Precipitación media mensual…………………………...……………………………………….. 56
6.6 Delimitación de subcuencas y red de drenaje……………….…………………………………58
6.7 Determinación del número, longitud y área de corrientes para orden 1, 2, 3 y 4 de la
subcuenca San Francisco……………………………………………...……………………………….58
6.8 Delimitación de microcuencas de orden 1, 2, 3 y 4 de la subcuenca San Francisco…...... 59
6.9 Hidrogramas unitarios por el método del HUIG……………..…………………………...…….60
6.10 Hidrogramas unitarios por el método del HUA………………….……………………...……. 61
6.11 Estaciones climatológicas y sus área de influencia…………………………………………. 63
6.12 Ajuste de curvas tipo a los FRD……………………………………………………….………. 64
6.13 Hietogramas base para duraciones de 1, 2 y 3 horas…………………………………...…..65
6.14 Uso de suelo y vegetación………………………………………………….………………….. 65
6.15 Coeficientes de escurrimiento…………………………………………………….…………… 67
6.16 Ubicación espacial de las estaciones y PCR en la cuenca del río Sabinal….…………….67
6.17 Líneas de vista entre las estaciones y el PCR-02 con el PCR-01………………….……… 69
6.18 Instrumentación de las estaciones pluviométricas……………………………………….….. 70
6.19 Instrumentación de las estaciones hidrométricas.……………………………………….….. 71
6.20 Instrumentación de los PCR………………………..……...……………………………….…..72
6.21 Áreas de influencia para la cuenca total y subcuencas………………..…………………… 72
6.22 Áreas de influencia suponiendo que no funcionan 1, 2, 3, 4 y 5 estaciones al mismo
tiempo……………………………………………………………………………….……………………. 74
6.23 Malla de celdas de 500 m por lado ajustada al perímetro de las subcuencas……..…….. 75
6.24 Ubicación de las secciones a lo largo del río Sabinal………………………………….……. 77
6.25 Perfil del thalweg, bordo derecho e izquierdo…………………………………….………….. 78
6.26 Características y condiciones actuales del río Sabinal en diferentes secciones………… 79
6.27 Hidrogramas de escurrimiento directo…………………………………………………………81
6.28 Esquema idealizado del funcionamiento hidráulico del río…………………………………. 81
6.29 Simulación en flujo permanente para las secciones con desbordamiento…….....………. 83
6.30 Programa de simulación en flujo no permanente del II, UNAM…………….……………… 83
6.31 Simulación para un tr = 5 años……………………...………………………………………….84
6.32 Descripción de la pantalla del programa de cómputo para la recopilación,
almacenamiento y envío de datos……………………………………………………………..…..…..85
6.33 Descripción de la pantalla del programa de cómputo para el procesamiento y análisis de
datos…………………………………………...…………………………………………………………. 86
6.34 Hietogramas de precipitación de las estaciones pluviométricas………….................……. 87
6.35 Curvas masa de las estaciones pluviométricas……………………………………………… 88
6.36 Hietogramas de precipitación media de las subcuencas…………………………………… 89
6.37 Curvas masa de las subcuencas……..………………………..……………………………… 89
6.38 Hidrogramas de escurrimiento directo……………………………………………………..…. 91
6.39 Resultados de la simulación en flujo no permanente…………………………………….…. 92
6.40 Hidrogramas de entrada y salida del cauce……………………………………………..…… 92
7 CONCLUSIONES
8 BIBLIOGRAFÍA Y REFERENCIAS
ANEXO A
vi
Índice de tablas
ÍNDICE DE TABLAS
1 INTRODUCCIÓN
4 PROCESAMIENTO HIDRÁULICO
6 CASO DE APLICACIÓN
7 CONCLUSIONES
8 BIBLIOGRAFÍA Y REFERENCIAS
ANEXO A
viii
Introducción
1 INTRODUCCIÓN
El agua es el recurso natural más importante de cualquier país debido a los beneficios sociales
y económicos que se derivan de su conciente explotación; sin embargo, junto con sus ventajas
también existen desventajas, como cuando su exceso origina inundaciones.
Una inundación se produce cuando una zona terrestre queda momentáneamente cubierta por el
agua. Ello se debe al almacenamiento de agua de lluvia en zonas con escaso drenaje
(inundaciones pluviales), al desbordamiento de arroyos o ríos (inundaciones fluviales), a mareas
de tormenta (inundaciones costeras), a la descarga de obras de excedencias de presas o bien a
la falla de una estructura hidráulica.
La República Mexicana se encuentra sujeta a gran variedad de fenómenos que pueden causar
desastres, como las precipitaciones convectivas, frentes fríos o ciclones tropicales que se
generan tanto en el Océano Atlántico como en el Pacífico, debido a su ubicación geográfica en
una región intertropical y características climáticas.
Los desastres por estos fenómenos han causado un importante número de victimas fatales y
cuantiosas pérdidas materiales no solo a nivel mundial sino también en nuestro país, que se
caracterizan por su frecuencia y carácter recurrente cada vez incrementándose más, los daños
que ocasionan.
Tabla 1.1 Desastres hidrometeorológicos más importantes en México, entre 1943 y 2004
Evento Año No. de Muertos Daños (MDD)
Huracanes varios 1954 110 7.50
Ciclón de Manzanillo 1959 1,500 No cuantificado
Huracán Liza 1976 600 3.10
Huracán Diana 1990 139 90.70
Inundaciones en Chihuahua 1990 200 2.50
Tormentas de invierno en Nayarit 1992 64 78.00
Huracán Ismael 1995 56 26.00
Huracán Paulina 1997 228 447.80
Lluvias torrenciales en Chiapas 1998 229 602.70
Lluvias torrenciales en Tijuana 1998 92 65.60
Inundaciones en Veracruz 1999 124 293.40
Inundaciones en Puebla 1999 263 244.80
Otros 3,778 4,134.40
Total 5,173 5,985.90
MDD. Millones de dólares.
(Bitrán, 2000 y Salas, 2005). En la tabla 1.1 se enlistan los eventos más importantes en este
periodo y en la figura 1.1 se representa el número de personas fallecidas en los años en que
ocurrieron dichos eventos.
1,400
1,000
800
600
400
200
0
1943
1949
1955
1959
1960
1967
1968
1976
1980
1982
1985
1987
1988
1990
1991
1992
1993
1995
1996
1997
1998
1999
2000
2001
2002
2003
2004
Tiempo (Año)
Históricamente, los estados que han sido afectado por estos fenómenos que se caracterizan
por afectar un área geográfica extensa son: los once estados de la costa del océano pacifico,
los seis de atlántico y cinco del centro del país (Coahuila, Hidalgo, Nuevo León, Puebla y San
Luis Potosí), que en área abarcan más del 75 % del territorio nacional, con una población de
más de 64 millones de personas de las cuales alrededor del 30 % están en riesgo por estos
fenómenos (figuras 1.2 y 1.3).
Los devastadores efectos de estos eventos, se magnifican por una elevada vulnerabilidad que
afecta a la sociedad que se deriva entre otros, de la insuficiencia de medidas preventivas y de
2
Introducción
Figura 1.3 Localidades con más de 10,000 habitantes ubicadas en cuencas con peligro de
inundación
(fuente: CENAPRED)
Entre los factores que condicionan las consecuencias de las inundaciones están la distribución
espacial y temporal de la lluvia, la topografía, pendiente del terreno, el tipo y uso de suelo,
cobertura vegetal, las características físicas de los cauces de los arroyos y ríos y la ubicación
de las estructuras hidráulicas.
El tiempo de respuesta de una cuenca depende entre otros factores de las características
fisiográficas y geomorfológicas de la misma; por ello, básicamente las inundaciones se
clasifican en lentas y rápidas. Las primeras significan que se generan avenidas en un tiempo
relativamente largo (del orden de varias horas o días); en ellas ocurren principalmente daños
materiales, mientras que, cuando la inundación es rápida o súbita, esta se forma en poco
tiempo (desde unos minutos hasta unas cuantas horas), causando la pérdida principalmente de
vidas humanas y desde luego daños materiales.
1.1 Justificación
Para prevenir y mitigar futuros desastres en estos sitios y en muchos otros, es indispensable
emprender acciones de protección. Estas pueden ser de dos tipos: estructurales (construcción
de obras) y no estructurales (indirectas o institucionales). El objetivo de las obras estructurales
3
Introducción
es reducir los daños provocados por las inundaciones mediante la construcción de obras; por su
parte, la finalidad de las medidas no estructurales es informar oportunamente a las poblaciones
ribereñas de la ocurrencia de una avenida, para que no haya decesos de personas y se
minimicen los daños materiales. Dentro de esta última acción, se encuentra el Sistema de Alerta
Hidrometeorológica, propósito de este trabajo.
1.2 Objetivos
En primer lugar, describir los aspectos técnicos más importantes que se deben de considerar en
la implantación de un Sistema de Alerta Hidrometeorológica Temprana Municipal.
En segundo, su aplicación a la cuenca del río Sabinal en Tuxtla Gutiérrez, Chiapas, localidad
que ha sufrido inundaciones súbitas en los últimos años, con grandes daños materiales y
contribuir para que en un futuro no existan pérdidas de vidas humanas y disminuir las
materiales.
1.3 Contenido
El segundo capítulo trata de los aspectos hidrológicos como las características fisiográficas, la
precipitación y pérdidas de la misma, modelos para transformar la lluvia en escurrimiento, con lo
cual es posible determinar el escurrimiento en respuesta a una tormenta de una zona en
estudio.
En muchos casos, los asentamientos humanos están ubicados en el cauce de ríos que reciben
el agua de varias subcuencas, es por ello necesario contar con modelos matemáticos para
conocer el comportamiento del agua dentro de dicho cauce, como se describe en el capítulo
cinco.
Con la conjunción de los elementos de los capítulos anteriores se hace una aplicación a la
cuenca del río Sabinal en Tuxtla Gutiérrez, Chis.
4
Estimación de la magnitud de los escurrimientos en un río
En este capítulo se hace una descripción de los conceptos necesarios para estimar la magnitud
de los escurrimientos en un río. Para disponer de mayor tiempo para avisar de una condición de
peligro se utilizan modelos que relacionan la lluvia con el escurrimiento.
A Área de la cuenca
B Corriente principal
La corriente principal de una cuenca es la que pasa por su salida y para determinarla se parte
del punto de salida de la cuenca hacia aguas arriba, siguiendo a la corriente de más alto orden,
hasta alcanzar una bifurcación de dos corrientes de igual orden y se continúa por el cauce que
tenga una mayor área de aportación; a partir de tal punto el proceso se repite hasta terminar en
un tributario de orden uno.
La longitud del cauce principal es la distancia horizontal que recorre la corriente principal.
La pendiente del cauce principal es el desnivel que existe entre los puntos inicial y final del
mismo, divido por su longitud horizontal. Dado que la pendiente varía a lo largo del cauce es
necesario definir una pendiente media y para ello se utiliza el criterio de Taylor y Schwarz
definido por la ecuación siguiente:
5
Estimación de la magnitud de los escurrimientos en un río
2
⎡ ⎤
⎢ ⎥
=⎢ ⎥
L
Sm
⎢ l1 l2 l ⎥
⎢ + + ... + i ⎥
⎣⎢ S1 S2 Si ⎦⎥
donde:
Sm , pendiente media del cauce principal, adimensional.
L , longitud total del cauce, en m.
h
S i = i , pendiente del tramo i , adimensional.
li
hi , diferencia de elevaciones del tramo i , en m.
l i , longitud del tramo i , en m.
i , número de tramos que hay en la longitud L .
La red de drenaje de una cuenca esta integrada por una corriente principal y una serie de
tributarios cuyas ramificaciones se extienden hacia las partes más altas de la cuenca. Las
corrientes tributarias son las corrientes secundarias que llegan a la corriente principal.
3 2
3 3 1
2
4
4
1
1
F Leyes de Horton
Horton (1945) fue el primero en sugerir algunas leyes empíricas que se consideran constantes
de un orden a otro, para que posteriormente fueran confirmadas usando el sistema de
ordenación de Strahler.
Ni
Ley del número de corrientes RB =
N i +1
6
Estimación de la magnitud de los escurrimientos en un río
Li +1
Ley de las longitudes de las corrientes RL =
Li
A
Ley de áreas de las corrientes R A = i +1
Ai
donde:
n
∑ Li , longitud promedio de una corriente de orden i , en km.
1
Li =
n i =1
n
∑A
1
Ai = , área promedio que contribuye al escurrimiento de una corriente de orden
n i =1 i
i , en km2.
Li , longitud de una corriente de orden i , en km.
Ai , área que contribuye al escurrimiento de una corriente de orden i , en km2.
N i , número de corrientes de orden i .
En realidad, las leyes de Horton indican una progresión geométrica de número, longitud y área
de las corrientes de una cuenca y por lo tanto, gráficamente corresponden a las relaciones
lineales entre el número de orden y los logaritmos del número, longitud y área de corrientes; de
donde las relaciones R B , RL y R A se obtienen como el antilogaritmo de la pendiente de dichas
rectas. A través de los años se ha observado que para cuencas naturales tienen cierto rango
de variación, para R B de 3 a 5, para RL de 1.5 a 3.5, y para R A de 3 a 6.
2.2 Precipitación
En la actualidad existen una gran variedad de instrumentos y técnicas para obtener información
de las diferentes fases de la precipitación. Sin embargo, los aparatos que miden la cantidad y la
intensidad de la lluvia son los más importantes en la mayoría de los casos. Todas las formas de
precipitación se miden sobre la base de una columna vertical que se acumularía sobre una
superficie a nivel si ésta permaneciese en el lugar donde cae.
Los aparatos más usuales en nuestro país para medir la precipitación en forma de lluvia son
los pluviómetros y los pluviógrafos, aunque últimamente se están usando los radares
meteorológicos e imágenes de satélite.
7
Estimación de la magnitud de los escurrimientos en un río
colocado dentro del recipiente medidor; el ascenso del flotador se produce por el aumento de
lluvia captada, hasta que el recipiente se llena y entonces se vacía automáticamente por medio
de un sifón autocebante. El flotador se une mediante un juego de varillas a una plumilla, la cual
registra las marcas de altura de precipitación sobre un papel graduado, este último colocado
sobre un tambor que gira a velocidad constante (figura 2.2b). Estos instrumentos miden la
precipitación a lo largo del tiempo. El registro obtenido se llama pluviograma.
b)
a)
Acotaciones en cm.
Figura 2.2 Instrumentos más utilizados en México para medir la precipitación
En los últimos años en México se ha iniciado con el uso de pluviográfos de tipo balancín, que
son instrumentos automáticos y realizan la medición de la precipitación a cada de 10 minutos.
A Curva masa
Es la curva que relaciona la altura de precipitación acumulada a través del tiempo de una
tormenta dada; se obtiene del registro de un pluviógrafo (figura 2.3).
B Altura de precipitación
La altura de precipitación o lámina de lluvia es la cantidad de agua que se precipita durante una
tormenta en un determinado tiempo. Sus unidades son de longitud, generalmente se expresa
en mm.
8
Estimación de la magnitud de los escurrimientos en un río
C Duración
Existen una gran cantidad de ecuaciones para estimarlo en función de la longitud y la pendiente
del cauce principal como la de Kirpich, que se expresa de la forma siguiente:
0.77
⎡ L ⎤
t c = 0.000325 ⎢ ⎥
⎣ S⎦
donde:
t c , tiempo de concentración, en h.
L , longitud del cauce principal, en m.
S , pendiente del cauce principal, adimensional.
D Intensidad
40
30
20
Δh p
10 i =
Δt
0 2 4 6 8 10 12
t (h)
d (h)
Figura 2.3 Curva masa de una tormenta
E Hietograma
F Precipitación media
posteriormente obtener un promedio pesado con respecto al área total, intentando con ello
tomar en cuenta la distribución espacial de la lluvia (figura 2.5).
i (mm/h)
hp (mm)
Δt Δt
t (h) t (h)
a) Hietograma, expresado en b) Hietograma, expresado en
altura de precipitación intensidad de lluvia
Figura 2.4 Hietograma de una tormenta
n
∑ Ai hP
1
hP = i
AT i =1
donde:
hP , altura de precipitación media, en mm.
hPi , altura de precipitación registrada en la estación i, en mm.
AT , área total de la cuenca, en km2.
Ai , área de influencia de la estación, en km2.
n , número de estaciones bajo análisis.
Estaciones
Ai
En general, las pérdidas están constituidas por la retención de la precipitación por la vegetación
y estructuras (intercepción), una parte llega al suelo y otra se evapora; otra parte de ella se
infiltra (infiltración y transpiración), otra parte llena las pequeñas depresiones naturales del
terreno (almacenamiento en depresiones) y luego se evapora (evaporación). De esta forma, las
pérdidas de la precipitación en una cuenca se pueden cuantificar como la diferencia entre el
volumen de precipitación y el que escurre por su salida.
Para la estimación de la avenida que se puede presentar en un punto determinado de un río por
consecuencia de una lluvia, se realiza un estudio de la precipitación y el escurrimiento con base
en los hietogramas e hidrogramas existentes, seguidos de un conocimiento minucioso de la
capacidad de infiltración de la cuenca. La evaporación y transpiración se desprecian debido al
corto tiempo en que se desarrolla el proceso, la intercepción y el almacenamiento en
depresiones se ignora debido a su reducida magnitud con respecto al que se analiza.
10
Estimación de la magnitud de los escurrimientos en un río
2.3.1 Infiltración
La infiltración es el proceso mediante el cual el agua penetra desde la superficie del terreno
hacia el suelo, por la acción de las fuerzas gravitacionales y capilares. El agua primero satisface
la deficiencia de humedad del suelo y queda retenida en él, después cualquier exceso alcanza
un manto acuífero. Por lo general, la infiltración presenta un valor inicial alto el cual disminuye
durante las lluvias continuas, hacia una magnitud reducida y prácticamente constante (figura
2.6).
La cantidad máxima de agua que puede absorber el suelo durante una tormenta en
determinadas condiciones se llama capacidad de infiltración y se expresa normalmente en
mm/h. Durante una tormenta sólo se satisface la capacidad de infiltración mientras ocurre la
lluvia en exceso. El volumen infiltrado es igual al área bajo la curva de capacidad de infiltración,
para cualquier intervalo y representa la lámina de agua infiltrada durante dicho periodo (figura
2.6) (Campos, 1988).
i ≥f
Capacidad de infiltración (f), en mm/h
3
fc
Tiempo, en min
Figura 2.6 Proceso de infiltración
VP = VPT − VED
donde:
VP , volumen de pérdidas, en m3.
VPT , volumen de precipitación total, en m3.
VED , volumen de escurrimiento directo, en m3.
VED
Ce =
VPT
donde:
Ce , coeficiente de escurrimiento, adimensional.
VED , volumen de escurrimiento directo, en m3.
VPT , volumen de precipitación total, en m3.
Donde evidentemente, este cociente tendrá valores entre 0 y 1. Será 0 si toda la precipitación
se convierte en pérdidas y será 1 si la cuenca es totalmente impermeable, además es variable
en el tiempo, de acuerdo a la condición de humedad antecedente.
El Soil Conservation Service (1972) desarrolló un método para calcular las pérdidas de
precipitación de una tormenta, en donde el número de escurrimiento depende del tipo, uso, y
tratamiento del suelo, pendiente y estado de humedad del terreno. La precipitación total se
relaciona con la precipitación en exceso mediante las curvas de la figura 2.7 y se pueden
expresar algebraicamente como:
12
Estimación de la magnitud de los escurrimientos en un río
2
⎛ 5,080 ⎞
⎜P − + 50.8 ⎟
= Pe = ⎝ ⎠
N
Qed
20,320
P+ − 203.2
N
donde:
Qed , escurrimiento directo, en mm.
Pe , precipitación efectiva, en mm.
P , precipitación total, en mm.
N , número de escurrimiento, adim.
El número de escurrimiento N, se determina con la estimación del tipo de suelo que se muestra
en la tabla 2.2 y en función del uso, tratamiento del suelo y pendiente del terreno con ayuda de
la tabla 2.3.
Para tomar en cuenta la condición de humedad antecedente, se hace una corrección al número
de escurrimiento, según la altura de precipitación acumulada cinco días antes de la fecha en
cuestión, utilizando la tabla 2.4.
En ambos métodos, cuando el área de la cuenca está constituída por superficies diferentes, el
coeficiente de escurrimiento puede obtenerse en función de las características de cada porción
del área como un promedio pesado, es decir:
Ce1A1 + Ce 2 A2 + ... + Ce i Ai
Ce =
AT
Donde:
Ce , Coeficiente de escurrimiento ponderado, adim.
Cei , coeficiente de escurrimiento correspondiente al área Ai , adim.
Ai , área de superficie i que tiene un cierto tipo de superficie, en km2.
AT , área de la cuenca, en km2.
2.4 Escurrimiento
13
Estimación de la magnitud de los escurrimientos en un río
es el que más rápido llega hasta la salida de la cuenca, por ello está relacionado directamente
con una tormenta en particular y entonces se dice que proviene de la precipitación efectiva y
por lo tanto constituye el escurrimiento directo.
Actualmente existe un gran número de técnicas para estimar el escurrimiento o gasto de un río
en una sección e instante determinados. Por lo general, estas mediciones se llevan a cabo en
una estación hidrométrica. Los métodos disponibles pueden ser incluidos en alguno de los
grupos siguientes:
a) Sección de control.
b) Relación sección-pendiente.
c) Relación sección-velocidad.
El método más adecuado, depende principalmente del tamaño del río, la precisión requerida,
del tiempo y equipo disponible.
Una vez definido el método y sección adecuados, se realizan varios aforos y se construye una
curva elevaciones-gastos, ajustando los puntos medidos a una curva media que tiene una
ecuación del tipo:
Q = C (E − E 0 )n
donde:
Q , gasto que pasa por la sección, en m3/s.
C y n , constates del ajuste.
E , elevación, en msnm.
E 0 , elevación para la que el gasto es nulo, en msnm.
Limnímetro. Es una regla graduada que se coloca en una de las márgenes del cauce, en la que
normalmente se lee la elevación de la superficie libre del agua en intervalos de tiempo iguales.
Limnígrafo. Es un instrumento que consta de un flotador unido a una plumilla que marca el nivel
de la superficie libre del agua en el río en un papel graduado, fijado a un tambor que gira a
15
Estimación de la magnitud de los escurrimientos en un río
El hidrograma es una gráfica o tabla que representa los gastos con respecto al tiempo, es decir,
que relaciona el gasto que se mide al pasar de manera continua durante un periodo de tiempo
por una sección transversal de un río. Dos tipos de hidrogramas son particularmente
importantes: el hidrograma anual y el hidrograma de una tormenta.
El hidrograma asociado a una tormenta aislada se muestra en la figura 2.8, los cuales varían de
tormenta a tormenta y de una cuenca a otra, pero en general es posible distinguir las partes
siguientes:
A , inicio del escurrimiento directo.
B , pico del hidrograma.
C , final del escurrimiento directo e inicio del escurrimiento base.
Q p , Gasto al pico, en m3/s.
Qb , Gasto base, en m3/s.
t p , tiempo al pico del hidrograma, en h.
t b , tiempo base, en h.
Q(m³/s)
Curva de
Rama
ascendente recesión
Vólumen de escurrimiento
directo (Ve)
Qp
C
A
Qb
t(h)
tp
tb
En la actualidad existen una gran cantidad de métodos que se pueden clasificar de acuerdo a la
información que requieren para su aplicación, que van desde los que necesitan una sola
característica fisiográfica hasta los que utilizan, características fisiográficas, datos de
precipitación y de escurrimiento.
Para su análisis, el viaje de una gota está gobernado por las reglas siguientes:
1. Para una gota que cae en una ladera, su estado correspondiente es e i , donde i es el
orden de la corriente asociada.
2. Del estado e i necesariamente se pasa al estado r j correspondiente.
3. De un estado r j , se puede pasar a cualquier estado rk si k > j .
4. Necesariamente se pasa por rn y de ahí con probabilidad de 1 al estado n + 1 que es la
salida.
Las reglas anteriores definen un conjunto finito de trayectorias que puede seguir una gota hasta
alcanzar la salida de la cuenca. Así, para una cuenca de orden 3, las trayectorias posibles son:
S1 : a1 r1 r2 r3 salida.
S2 : a1 r1 r3 salida.
S3 : a2 r2 r3 salida.
S4 : a3 r3 salida.
Con estas condiciones, Gupta (1980) demostró que el HUIG es igual a la función de densidad
de probabilidad del tiempo de viaje TB a la salida de la cuenca de una gota de agua que cae en
ésta aleatoriamente con distribución espacial uniforme y está dada por:
P (TB ≤ t ) = ∑ P (TS i
)
≤ t P (S i )
Si ∈S
donde:
17
Estimación de la magnitud de los escurrimientos en un río
El tiempo de viaje TSi de una trayectoria en particular, es igual a la suma de los tiempos de
viaje de los elementos de esa trayectoria. Así, para las trayectorias anteriores sus
correspondientes tiempos de viaje son:
Por lo tanto, la función de densidad de probabilidad del tiempo de viaje total de una trayectoria
TS i , está dada por la convolución de las funciones de densidad de cada tramo de la
trayectoria, esto es:
fTS (t ) = fTr (t ) * fTr (t ) * ... * fTr (t ) * fTr (t )
i 1 2 n n +1
Así, para una cuenca de orden tres, la función de distribución de probabilidad de los tiempos de
viaje es:
( ) ( ) ( )
P (TB ≤ t ) = P TS1 ≤ t P (S1 ) + P TS2 ≤ t P (S 2 ) + P TS3 ≤ t P (S 3 )
Donde la probabilidad de los tiempos de viaje de cada una de las trayectorias se obtiene
realizando sus respectivas convoluciones.
Rodríguez-Iturbe y Valdés (1979) consideran que los tiempos de viaje de una gota en una
corriente de orden ω obedecen a una función de densidad de tipo exponencial, dada por:
fTr (t ) = K ω e (−Kω t )
ω
donde:
1
, tiempo promedio de escurrimiento, en h.
Kω
v v K K
Kω = , K1 = , K 2 = 1 , K 3 = 12
Lω L1 RL RL
V , velocidad de pico para cualquier evento dentro de la cuenca, en m/s.
Este método se considera adecuado para cuencas cuyo orden de red de drenaje es
representativo de la misma, esto se logra para un orden igual o mayor a tres.
Cuando para una cuenca el orden de la cuenca es menor a tres se recomienda el uso del
hidrograma unitario sintético adimensional (HUA). El cual fue propuesto por el Soil Conservation
Service (1975), a partir de varios hidrogramas registrados en una gran variedad de cuencas, en
el cual los gastos se expresan en función del gasto de pico, y los tiempos en función del tiempo
al pico (figura 2.9).
18
Estimación de la magnitud de los escurrimientos en un río
Para obtener un hidrograma unitario a partir del HUA para una cuenca, se multiplican las
abcisas y ordenadas del hidrograma unitario sintético adimensional por el gasto de pico y
tiempo al pico respectivamente, en donde estos últimos se calculan con las expresiones
siguientes:
0.208 A d
qP = t P = e + 0.6t c
tP 2
donde:
q P , gasto de pico en m3/s/mm.
A , área de la cuenca, en km2.
t P , tiempo al pico, en h.
t c , tiempo de concentración, en h.
d e , duración en exceso, en h.
1 .0 0 t/tp q/qp
0.0 0.00
0.1 0.03
0 .9 0
0.2 0.10
0.3 0.19
0 .8 0 0.4 0.31
0.6 0.66
0.8 0.93
0 .7 0 1.0 1.00
1.2 0.93
0 .6 0 1.4 0.78
1.6 0.56
q / qp
1.8 0.39
0 .5 0 2.0 0.28
2.2 0.207
2.4 0.147
0 .4 0
2.6 0.107
2.8 0.077
0 .3 0 3.0 0.055
3.2 0.040
3.4 0.029
0 .2 0
3.6 0.021
3.8 0.015
0 .1 0 4.0 0.011
4.2 0.010
4.4 0.007
0 .0 0 4.6 0.003
0 .0 0 .5 1 .0 1 .5 2 .0 2 .5 3 .0 3 .5 4 .0 4 .5 5 .0 4.8 0.0015
t / t p 5.0 0.00
Figura 2.9 Hidrograma unitario sintético adimensional
Cuando la curva de recesión del hidrograma no se considera importante, el HUA puede ser
modificado con un tiempo base igual a cuatro veces el tiempo al pico (figura 2.10).
t/tp q/qp
y = x (0.6704 + x(2.6905 - 2.037x)) - 0.3249 0.00 0.00
0.17 0.05
1.0 0.33 0.20
0.50 0.43
0.9 0.67 0.71
y = x(7.948 + x(-6.0952 + x(1.3889))) - 2.24 0.83 0.92
0.8 1.00 1.00
1.17 0.94
0.7
1.33 0.81
1.50 0.66
0.6
1.67 0.51
q/qp 1.83 0.40
0.5
2.00 0.32
0.4 2.17 0.25
y = x(-4.2829 + x(1.8039 +x(-0.3532 + 0.0266x))) + 4.0679 2.33 0.20
0.3 2.50 0.16
2.67 0.12
0.2 2.83 0.09
y = x(-0.0236 + x (2.0417 - x 0.5556)) 3.00 0.07
0.1 3.17 0.05
3.33 0.04
3.50 0.02
0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 3.67 0.13
t/tp 3.83 0.01
0.67 0.00
Figura 2.10 Hidrograma unitario sintético adimensional modificado
19
Estimación de la magnitud de los escurrimientos en un río
Una vez conocidos el hietograma de una tormenta y el hidrograma unitario obtenido por
cualquiera de los métodos anteriores, el siguiente paso consiste en calcular el hidrograma de
escurrimiento directo. Utilizando el concepto del hidrograma unitario y la ecuación de
convolución que en su forma discreta se expresa como:
n ≤m
Qn = ∑ PmU n−m+1
m =1
de de de de de de
1
t (h)
q(m³/s/mm)
U2
a)
U1 U3
t (h)
hp(mm)
P2
b)
P1
P3
t (h)
Q(m³/s)
P 1U 2
c)
P 1U 1 P 1U 3
+
P 2U 2 t (h)
Q(m³/s)
P 2U 1 P2U3 d)
+
t (h)
Q(m³/s)
P3U2
e)
P 3U 1
P 3U 3
= Q4
Q3 t (h)
Q(m³/s)
Q2 f)
Q5
Q1
t (h)
Figura 2.11 Aplicación de la ecuación de convolución discreta
20
Estimación de la magnitud de los escurrimientos en un río
Q1 = P1U1
Q2 = P1U 2 + P2U1
Q3 = P1U 3 + P2U 2 + P3U1
Q4 = 0 + P2U 3 + P3U 2
Q5 = 0 + 0 + P3U 3
De donde, el número de ordenadas del hidrograma final NQ está ligado con el número de
barras del hietograma de precipitación efectiva N P , y al número de ordenadas del hidrograma
unitario NU por medio de la ecuación:
NQ = N p + NU − 1
donde:
Qi , ordenada del hidrograma de escurrimiento directo para el i-ésimo intervalo de
tiempo.
Pi , i-ésima altura de precipitación efectiva.
U i , i-ésima ordenada del hidrograma unitario de duración en exceso d e .
NQ , número total de ordenadas del hidrograma de escurrimiento directo.
NU , número total de ordenadas del hidrograma unitario de duración en exceso d e .
N P , número total de barras del hietograma de lluvia efectiva, con duración en exceso
d e cada una de ellas.
El periodo de retorno se define como el número promedio de años que transcurren para que un
evento de cierta magnitud sea igualado o excedido por lo menos una vez (si los eventos son
anuales), esto es:
1 1
Tr = =
P ( X > x ) [1 − P ( X ≤ x )]
Por otra parte, la probabilidad de que un evento con un periodo de retorno sea igualado o
excedido en cualquier año es:
1
P=
Tr
En un conjunto de eventos máximos anuales, el periodo de retorno que se asocia a cada uno
de ellos puede ser estimado con la formula de Weibull.
n +1
Tr =
m
Donde:
n , número total de datos o años de registro.
m , número de orden. El orden es el lugar asignado en una lista de valores de la
muestra, ordenados de mayor a menor.
2.9 Tormentas de diseño
21
Estimación de la magnitud de los escurrimientos en un río
22
Instrumentación de la zona de estudio
3 INSTRUMENTACIÓN DE LA
ZONA DE ESTUDIO
Para el estudio de un fenómeno hidrometeorológico es necesario que se disponga con
información de los mismos; para obtener dicha información es necesario hacer mediciones de
algunos aspectos del fenómeno con las técnicas e instrumentos adecuados, para que
posteriormente sirvan de datos en modelos de simulación hidrológica e hidráulica. Los avances
de la electrónica permiten que la información se mida, transmita, almacene y se analice casi al
mismo tiempo en que ocurre el evento con diferentes propósitos, como por ejemplo la
predicción de avenidas.
23
Instrumentación de la zona de estudio
La selección de los lugares en los que se deben instalar las estaciones de medición de
precipitación y del nivel de la superficie libre del agua en un río se fundamenta tanto en criterios
hidrológicos como topográficos, los cuales se comentan a continuación:
Igualdad. Las áreas de influencia correspondiente a cada estación deben tener un peso similar
con respecto al área total de la cuenca.
Visibilidad. Para tener enlace de telemetría entre las estaciones pluviométricas, repetidores de
señales de radio y el PCR, o en su caso entre dos estaciones de las cuales una de ellas debe
de contar con una antena repetidora para la retransmisión de datos hacia el PCR, es importante
que exista línea de vista. Una línea de vista, es un perfil topográfico entre dos puntos que
permite principalmente tener visibilidad entre ellos y por consecuencia comunicación vía radio.
Accesibilidad. La estación debe ser accesible en la mayoría de las veces, por ello se debe de
procurar ubicarlas en áreas urbanas o en caso dado cerca de alguna vía de comunicación.
Importancia. Debe cubrir en orden de importancia las partes altas, medias y bajas de la cuenca
para tener un periodo de tiempo lo más amplio posible para dar aviso a las autoridades
correspondientes.
Permanencia. La estación debe estar situada de tal manera que nunca sea destruida por el
vandalismo o por un fenómeno hidrometeorológico, por ello es importante aprovechar sitios
seguros de la infraestructura existente como escuelas, tanques elevados, dependencias de
gobierno, etc.
3.1.2 Sitios de medición del nivel de la superficie libre del agua en un río
Permanencia. La estación debe estar situada de tal manera que nunca sea destruida por una
avenida. Una de las características más deseables de un registro es que sea continuo y que
esté formado en un mismo sitio. Además no debe estar afectado por tomas o desvíos, por lo
que la estación debe de situarse, en lo posible, aguas arriba de ellos y con fines de alerta es
necesario colocarlas aguas abajo de confluencias de una corriente y en sitios convenientes
para mejorar la calibración del modelo lluvia-escurrimiento.
Accesibilidad. La estación debe ser accesible en cualquier tiempo y bajo cualquier condición,
especialmente durante avenidas.
Suficiencia. La estación debe ser capaz de cubrir todo el rango de gastos que puedan ocurrir. El
nivel mínimo de la regla en los limnímetros, debe estar por debajo de la elevación
correspondiente al gasto nulo y la posición máxima del flotador o de la regla debe quedar arriba
de la elevación correspondiente al gasto máximo o capacidad máxima del cauce.
24
Instrumentación de la zona de estudio
Estabilidad. La sección transversal donde se instale la estación debe estar en un tramo recto
del río, lo más estable posible, de manera que las variaciones que tenga la curva elevaciones-
gastos sean razonablemente pequeñas.
Dada la importancia de operar las estaciones en forma confiable y continua, especialmente bajo
condiciones ambientales muy adversas, se diseñaron en forma robusta y con protecciones
especiales.
Los equipos se diseñan y dimensionan de manera que, aún bajo condiciones prolongadas de
intensos nublados en los cuales la celda es poco eficiente en su proceso de conversión de luz
solar en energía eléctrica, la estación pudiese operar sin problemas por lapsos de varias
semanas.
Figura. 3.4 Diagrama esquemático del interior del gabinete metálico y su interconexión con los equipos externos
(fuente: CENAPRED)
El pluviógrafo empleado es del tipo balancín, similar al mostrado en las figuras 3.5a y 3.5b.
Consiste de un recipiente cilíndrico en cuya parte superior se encuentra expuesto un embudo a
través del cual se capta el agua. Esta se conduce hacia una balanza o balancín con dos
recipientes (figura No. 3.5b). Cuando el agua acumulada en uno de los recipientes llega a una
cierta cantidad, el peso de la misma hace que el balancín salga de equilibrio, con lo cual se
vierte el agua y queda en posición para captar agua el otro recipiente. Este movimiento del
balancín activa momentáneamente un interruptor y produce una señal que es registrada como
un pulso o cuenta por la tarjeta de control y registro. Cada movimiento del balancín equivale por
tanto a una determinada cantidad de agua de lluvia (por lo general 0.25 mm.).
b)
Embudo
a) Balancín
Interruptor
Vertedores
Con el fin de conocer los gastos que fluyen en los ríos y prevenirse contra inundaciones
producidas por fenómenos hidrometeorológicos intensos, se emplea un aparato denominado
limnímetro para medir el nivel de la superficie libre del agua de los flujos que corren en los
cauces. A partir de este nivel, se inferirá el gasto del escurrimiento de las corrientes que
atraviesan aquellas secciones de interés.
El sensor de nivel consta de un tubo de PVC en cuyo interior se tienen puntos de prueba
equidistantes entre si, los cuales se fabrican a la distancia requerida. Se basan en la
conductividad del agua para operar una serie de circuitos con los que se conoce el nivel del
27
Instrumentación de la zona de estudio
En la figura 3.6 se muestra la instalación típica de una estación hidrométrica. Consta de una
caseta de concreto que aloja y protege al equipo electrónico de registro y comunicaciones. En
el techo de la caseta se coloca la celda solar que proporciona energía a los componentes de la
electrónica y un mástil para la antena de radio. Los principios de construcción son los mismos
que los utilizados para las estaciones pluviométricas
Para los enlaces de telemetría se requiere que exista línea de vista entre los sitios donde están
las estaciones pluviométricas e hidrométricas con el PCR, en algunos casos, debido a las
limitaciones físicas y topográficas, no es posible encontrar dicha línea, por lo que es necesario
hacerlo a través de un repetidor de señales de radio.
Un repetidor de señales de radio sirve como enlace entre las estaciones y el PCR, en ocasiones
cuando la topografía lo permite, esta función la puede realizar una estación pluviométrica.
Como se puede ver en el diagrama de la figura 3.7, el PCR primario está compuesto
principalmente por dos computadoras denominadas “Telemetría” e “Hidráulica”, sus periféricos
asociados, un panel de comunicaciones y una antena.
28
Instrumentación de la zona de estudio
Hidráulica
PC
Impresora
120 VCA
Conexión
de red
Antena omnidireccional
energía
Telemetría
Tierra Física
PC
Radio-módem
Interruptor general
(fuente: CENAPRED)
Conexión al sistema de
respaldo de energía
descargas
Protección
eléctricas
Gabinete
contra
+12 Vdc
Control de
alarma
Fuente de voltaje
120 Vac / 12 Vdc
+ 12 Vdc
Botón para
prueba de
Ventilador
alarma
+
+12 Vdc
deshabilitar alarma
Interruptor para
Alarma
29
Instrumentación de la zona de estudio
A fin de garantizar la operación continua del puesto central de registro, en caso de que llegue a
fallar el suministro de energía eléctrica, se cuenta con un sistema de respaldo que puede
mantener operando los equipos por espacio de 4 horas, aproximadamente. Asimismo se cuenta
con un regulador de voltaje, el cual mantiene estable la alimentación ante variaciones de
voltaje.
En caso de que existiera uno o más PCR secundarios se realizan las mismas actividades
excepto que no se interroga a las estaciones. En vez de ello espera recibir cada 10 minutos el
archivo con los datos de precipitación proveniente del PCR primario, una vez terminado un ciclo
de interrogación. Con estos datos procesa, analiza y presenta la información de la misma
manera a como lo hace la computadora “Hidráulica” en el PCR.
Las estaciones remotas fueron diseñadas como estaciones "inteligentes". Cada una tiene como
parte central un procesador que controla en forma automática y autónoma la medición y el
registro de la precipitación y almacena los datos temporalmente en su memoria hasta ser
solicitados.
La comunicación entre el puesto central de registro y las estaciones remotas es en forma digital,
es decir, a través del intercambio de mensajes vía radio-módem entre la computadora
denominada “Telemetría” del PCR primario y el procesador de cada estación. Para interrogar a
todas las estaciones de la red, la computadora encargada de la telemetría en el PCR primario
está programada con un algoritmo especial que se comunica a intervalos de diez minutos
(tiempo controlado por su reloj interno) con cada estación, para solicitarle su envío de datos.
Inicia con la estación No. 1, después de haber recibido correctamente los datos, se comunica
con la estación No. 2 y así sucesivamente con las demás estaciones de la red (figura 3.8). En
caso de que al primer intento de interrogación con una estación el PCR primario no logre recibir
sus datos, esperará unos segundos e intentará nuevamente establecer la comunicación, hasta
tres veces. Si no tiene éxito dejará pendiente esa estación y continuará con la siguiente
estación en turno. Como se verá más adelante, en los próximos ciclos de interrogación podrá
recuperarse esta información.
30
Instrumentación de la zona de estudio
sucesivamente; el intervalo No. 144 recoge los datos registrados entre las 7:50 y las 8:00 del
próximo día.
Todas las estaciones reciben este mensaje, pero sólo la estación que reconoce su identificación
responde trasmitiendo un paquete codificado de datos. En el caso de las estaciones de registro
de precipitación, este paquete consta de la siguiente información:
Este último es de gran importancia como signo vital de la estación para verificar su
funcionamiento.
Para el caso de las estaciones hidrométricas, el paquete de información es muy similar, solo
que en lugar de precipitación, contiene datos del nivel de la superficie libre del agua del río.
Una vez que los datos de un intervalo han sido recolectados, es decir, recibidos, decodificados
y desplegados en forma tabular, se almacenan en un archivo para ser procesados y analizados
por las computadoras denominadas “Hidráulica”. De esta manera se abren 144 archivos diarios.
proceso de interrogación, sólo recibe el archivo con todos los datos de precipitación ya
recolectados. Una vez recibidos, los analiza y los presenta en pantalla en la misma forma que
en el PCR primario.
32
Procesamiento hidráulico
4 PROCESAMIENTO
HIDRÁULICO
Con los métodos lluvia-escurrimiento es posible determinar el hidrograma en respuesta a una
precipitación en una sección del cauce de un río de una cuenca en estudio. En ocasiones es
necesario conocer el cambio que sufre dicho hidrograma al trasladarse en un tramo del río
dentro del cauce, hasta llegar a una sección aguas abajo; para lo cual se emplea el concepto de
tránsito de avenidas en cauces.
El tránsito de avenidas en cauces, se define como el procedimiento a través del cual se puede
determinar la evolución de un hidrograma en una sección dada y a lo largo del cauce de un río
o canal. De este modo se conoce la variación detallada del gasto, tirante, velocidad a lo largo
del cauce en diferentes tiempos.
Retraso
Q(m³/s)
Atenuación
Gasto de salida
Gasto de entrada
0
t (h)
Figura 4.1 Relación de gasto de entrada y salida en un tramo de cauce al presentarse una avenida
Las ecuaciones de la Mecánica de Fluidos que describen el movimiento de los líquidos o gases
son las ecuaciones de Navier-Stokes, pero debido a la dificultad para resolverlas se hacen
simplificaciones a las mismas.
Para el flujo el flujo en ríos o canales, como una aproximación aceptable desde el punto de vista
ingenieril, se usan las hipótesis de Saint-Venant, que se describen a continuación:
d) La pendiente del fondo del cauce o canal es pequeña y el lecho es fijo, esto es, no hay
erosión ni sedimentación del mismo.
e) El fluido es incompresible y de densidad constante a lo largo del flujo.
Para este caso, las ecuaciones fundamentales son las de continuidad, energía y cantidad de
movimiento pero sólo son necesarias dos (continuidad y cantidad de movimiento), con las
cuales se obtienen resultados mas generales. En su forma diferencial conservativa las
ecuaciones de Saint-Venant son:
Ecuación de continuidad:
∂A ∂Q
+ =0 (4.1)
∂t ∂x
∂Q ∂ ⎛⎜ Q 2 ⎞⎟ ∂y
+ + gA + gA(Sf − S0 ) = 0 (4.2)
⎜
∂t ∂x ⎝ A ⎠ ⎟ ∂x
donde:
y , tirante, en m.
A , área hidráulica, en m2.
Q , gasto, en m3/s.
g , aceleración de la gravedad, en m/s2.
S0 , pendiente del fondo, adim.
S f , pendiente de la línea de energía, adim.
x , distancia longitudinal, en m.
t , tiempo, en s.
El método de Muskingum, es un método hidrológico que fue presentado por McCarthy (1938) y
maneja relaciones gastos-almacenamientos variables. Este método usa una relación algebraica
lineal entre el almacenamiento, las entradas y las salidas mediante la combinación de
almacenamiento en forma de prisma y cuña (figura 4.2).
34
Procesamiento hidráulico
S = SP + SC = K 0 + KX (I − O ) = K [XI + (1 − X )O ] (4.3)
I Cuña de almacenamiento
(KX(I-O))
Prisma de almacenamiento O
(KO)
donde:
K , coeficiente de proporcionalidad y se considera igual al tiempo de tránsito de la onda
de la avenida a través del tramo.
X , factor de ponderación dentro del rango 0 ≤ X ≤ 0.5 .
S i +1 − S i = K [X (I i +1 − I i ) + (1 − X )(Oi +1 − Oi )] (4.4)
Si +1 − Si =
(I i + I i +1 ) (O + Oi +1 ) Δt
Δt − i (4.6)
2 2
Δt Δt Δt
KX + − KX + K (1 − X ) −
Δt
C1 = 2 C2 = 2 C3 = 2 α = K (1 − X ) +
α α α 2
35
Procesamiento hidráulico
A Información escasa
Cuando no existe información se puede decir en términos muy generales que X se aproxima a
0 en cauces muy caudalosos y de pendiente pequeña y a 0.5 en caso contrario. A falta de otros
datos es recomendable tomar 0.25 como un valor medio y K se hace igual al tiempo entre los
gastos de pico de los hidrogramas de entrada y de salida del tramo del río o bien se estima
(Linsley et al, 1975), como:
L1.6
K = 0.014 (4.8)
H 0 .6
donde:
H , desnivel entre las dos secciones de los extremos del tramo del río, en km.
L , longitud del tramo del río, en km.
B Información de campo
Cuando se dispone de los hidrogramas de entrada y salida para un tramo de un río los valores
de K y X pueden determinarse utilizando el método tradicional, bajo el procedimiento
siguiente:
Y = X (I i +1 − I i ) + (1 − X )(Oi +1 − Oi )
4. Se grafican los valores calculados en el paso 2 en el eje de las ordenadas y los del
paso 3 en el eje de las abcisas. El valor de X que produzca un ciclo lo más parecido a
una recta se utiliza para calcular el valor de K , que es la pendiente de dicha recta.
El programa de cómputo desarrollado, considera las variables necesarias para el cálculo del
flujo permanente gradualmente variado en cauces naturales. Este programa está integrado en
tres partes. La primera (preproceso), lee las condiciones iniciales o de frontera y los datos de
las secciones transversales; calcula las variables hidráulicas tales como son el perímetro
mojado, área, radio hidráulico, tirante crítico y factor de fricción para diferentes tirantes en cada
una de las secciones, y a partir de éstas forma una base de datos; además permite realizar una
interpolación entre secciones, cuando la distancia de separación entre ellas es razonablemente
grande. La segunda (proceso) realiza los cálculos necesarios para determinar el perfil del flujo.
36
Procesamiento hidráulico
La última parte (postproceso) consiste en mostrar en forma gráfica los perfiles del fondo, bordos
del cauce y del flujo, las secciones transversales del cauce mostrando la superficie libre del
agua, en planta las secciones que tienen derrames, y en forma tabular, para cada sección, la
distancia entre secciones, elevación del fondo del cauce, tirante, área, elevación de la superficie
libre del agua, velocidad, tirante crítico, número de Froude, bordo de la sección que desborda,
pérdidas de energía y la elevación de la línea de energía; esta información es guardada junto
con otros archivos con formato de texto que son los resultados de la simulación y que
posteriormente serán utilizados en el cálculo del flujo no permanente. A continuación se
describen los elementos teóricos que conforman el modelo matemático desarrollado.
El modelo matemático requiere información de las secciones transversales del tramo del río a
estudiar para realizar el cálculo hidráulico del flujo permanente gradualmente variado, como la
mostrada en la figura 4.3; se debe notar que se identifican los anchos Xd (i ) a la derecha de un
eje vertical que pasa por el punto más bajo de la sección, y por Xi (i ) a la izquierda del mismo
eje; estos anchos se han medido a equidistancias Δy del fondo.
Xi(n) Xd(n)
... Δy ...
Xi(3) Xd(3)
Xi(2) Xd(2)
Xi(1) Xd(1)
De acuerdo con esta nomenclatura, las expresiones para determinar las variables hidráulicas de
cada desnivel Δy en la sección j , se presentan en la tabla 4.1.
hf = Sf Δ l =
1
(
Sf + Sf2 Δ l
2 1
) (4.13)
donde:
4
V 2n 2 2 2⎛
p⎞ 3
Sf = =V n ⎜ ⎟ (4.14)
r
4
3 ⎝ A⎠
37
Procesamiento hidráulico
4
1 ⎛p⎞ 3 2 2
hf = Δ l n 2 ⎜ ⎟ V = FFV (4.15)
2 ⎝ A⎠
⎛ p (i ) ⎞
4
n2 Δ l 3
FF (i ) = ⎜⎜ ⎟⎟ (4.16)
2 ⎝ A (i ) ⎠
Para todos los valores del tirante y (i ) con sus respectivos valores de B (i ) , A(i ) , p(i ) y FF (i ) se
forma una base de datos que será utilizada posteriormente en el proceso de cálculo.
Generalmente, en los ríos de planicie (pendiente pequeña, Cosθ ≅ 1 ) donde ocurren problemas
de derrames se presenta régimen subcrítico, por lo que los tirantes en todas las secciones
transversales son mayores o iguales al tirante crítico y c . Por lo anterior, en los cálculos que se
exponen más adelante, tendrán como límite el tirante crítico y c , por lo cual es necesario
calcularlo previamente para cada sección.
La energía específica en una sección de un canal natural es la que corresponde al flujo por
unidad de peso a través de ella, y se mide con respecto al fondo de esta, por ello según la
figura 4.4, tenemos:
Q2
E = y +α (4.17)
2 g A2
La condición límite entre los dos regímenes alternos se llama sección crítica en un
escurrimiento a superficie libre y es aquella en la que la energía específica es la mínima posible
para que pueda fluir el gasto de dicho escurrimiento. Pueden encontrarse sus características
calculando la derivada respecto a y , de la expresión 4.17 y considerando que α también es
constante, obtenemos que:
dE α Q 2 dA
= 1− (4.18)
dy g A 3 dy
Podemos decir que para cualquier sección, el ancho de la superficie libre es B = dA dy , la cual
sustituyendo en la expresión 4.18, y aplicando el criterio de la primera derivada (dE dy ) = 0 , se
obtiene que en la sección donde la energía específica es mínima, se cumple que:
Q2 A3
α = (4.19)
g B
38
Procesamiento hidráulico
La igualdad anterior nos permite calcular el gasto con tirante crítico yc y α = 1, como:
g A (i )3
Q= (4.20)
B (i )
Por último, se puede resumir que en estas condiciones, para cada elevación i , existirán
relaciones de la forma:
A(i ) , es función del tirante y (i ) .
FF (i ) , es función del tirante y (i ) .
y c (i ) , es función del gasto Q .
Esta última relación existe, aunque de la forma Q = f (y ) , pero se puede hacer una relación
inversa, es decir, a partir del gasto que sea el caso de estudio y con las relaciones anteriores,
se puede encontrar el y c correspondiente a dicho gasto.
La energía total por unidad de peso que posee el agua en su movimiento medida entre las
secciones j y j + 1 de la figura 4.4, se puede determinar aplicando la ecuación de la energía,
de esta forma se obtiene:
v 2j +1 v 2j
y j +1 + + Δzj = y j + + hf (4.21)
2g 2g
yj+1 yj
2
Línea de hf
Vj+1 energía
2g
2
Superf Vj
icie de
l ag ua 2g
y j+1 Vj+1
Vj
Fondo yj
del cau
Δz c e
Para el término hf , podemos sustituir la combinación de las ecuaciones 4.13 y 4.14, de lo cual
resulta:
v 2j +1 v 2j j +1 v 2n 2
y j +1 +
2g
+ Δzj = y j +
2g
+ ∫j 4
3
dl (4.22)
r
⎡ 2 2⎤
⎛ ⎞ ⎛ ⎞
dl ⎢ 2 ⎥
y j +1 +
v 2j +1
+ Δzj = y j + + ⎢ v j +1
v 2j
( ) ⎜ n ⎟ ⎜ n
⎜ 2 ⎟ + vj ⎜ 2
2
( ) ⎟
⎟ ⎥ (4.23)
2g 2g 2 ⎢ ⎜r 3 ⎟ ⎜r 3 ⎟ ⎥
⎢⎣ ⎝ j +1 ⎠ ⎝ j ⎠ ⎥⎦
39
Procesamiento hidráulico
2 2
⎛ ⎞ ⎛ ⎞
y j +1 +
v 2j +1
( )
− v 2j + 1
⎜ n ⎟ dl
⎜ 2 ⎟ =yj +
v 2j
+ v 2j( ) ⎜ n
⎜ 2
⎟ dl
⎟ −Δzj (4.24)
2g ⎜r 3 ⎟ 2 2g ⎜r 3 ⎟ 2
⎝ j +1 ⎠ ⎝ j ⎠
Q2
De la ecuación de continuidad, v 2 = , por lo que sustituyendo en la ecuación 4.24, resulta:
A2
⎡ ⎛ ⎞
2 ⎤ ⎡ ⎛ ⎞
2 ⎤
⎛ Q 2 ⎞ ⎢ 1 ⎜ n ⎟ dl ⎥ ⎛ Q2 ⎞ ⎢ 1 ⎜ n ⎟ dl ⎥
y j +1 + ⎜ ⎟ ⎢ − ⎜ ⎟ ⎥ = y + ⎜ ⎟ ⎢ + ⎜ 2 ⎟ ⎥ −Δzj (4.25)
⎜ A 2 ⎟ ⎢ 2g ⎜ 2 3 ⎟ 2 ⎥ j ⎜ A 2 ⎟ ⎢ 2g ⎜r 3 ⎟ 2⎥
⎝ j +1 ⎠ r ⎝ j ⎠
⎢⎣ ⎝ j +1 ⎠ ⎥⎦ ⎢⎣ ⎝ j ⎠ ⎥⎦
2/3
1 ⎛P ⎞
Además, =⎜ ⎟ , sustituyendo en la ecuación anterior, y reordenando, se tiene
r 2/3 ⎝ A⎠
⎛ Q2 ⎞ ⎡ 4
3 2 ⎤ ⎛ Q2 ⎞ ⎡ 1
4 ⎤
⎛ p ⎞ ⎛ pj ⎞ 3
n 2 dl ⎥
= y j + ⎜ 2 ⎟ ⎢⎢
⎢ n dl ⎥
+⎜ 2 ⎟
1 j +1 ⎟
y j +1 ⎜
⎢ 2g − ⎜ A ⎟ + ⎜ ⎟ −Δzj (4.26)
⎜A ⎟
⎝ j +1 ⎠ 2 ⎥ ⎜ A ⎟ 2g ⎜ Aj
⎝
⎟
⎠ 2 ⎥
⎝ j +1 ⎠ ⎢⎣ ⎥⎦ ⎝ j ⎠ ⎢⎣ ⎥⎦
⎛ Q2 ⎞ ⎛ 1 ⎞ ⎛ Q2 ⎞ ⎛ 1 ⎞
y j +1 + ⎜ 2 ⎟ ⎜⎜ − FF j +1 ⎟⎟ = y j + ⎜ 2 ⎟ ⎜⎜ + FF j ⎟⎟ − Δ z j (4.27)
⎜ A ⎟ ⎝ 2g ⎠ ⎜ A ⎟ ⎝ 2g ⎠
⎝ j +1 ⎠ ⎝ j ⎠
K j +1 Kj
De igual forma, que para el flujo permanente, el programa de cómputo para el cálculo del flujo
no permanente se compone de tres partes esenciales: la primera (preproceso), lee las
condiciones iniciales o de frontera, las variables hidráulicas y realiza un cálculo de
calentamiento con un gasto base ambos calculados en el flujo permanente gradualmente
variado. La segunda (proceso) consiste en calcular las variables hidráulicas dentro del cauce
como son los perfiles del flujo, áreas, gasto, velocidades con respecto al tiempo. La tercera
(posproceso) permite ver en forma gráfica los hidrogramas de entrada al sistema de ríos,
mostrando el avance de cálculo de los mismos, los perfiles del fondo, bordos del cauce y del
flujo, las secciones transversales con su correspondiente nivel de superficie libre del agua, en
planta las secciones que tienen derrames, y en forma tabular, para cada sección el
cadenamiento, elevación del fondo del cauce, bordos, superficie libre del agua, tirantes, áreas,
velocidades, gastos de entrada al cauce, gastos de derrame derecho e izquierdo. La teoría que
justifica el modelo matemático se describe en las líneas siguientes.
Este método considera como flujo unidimensional al movimiento del agua en el cauce natural;
se limita al flujo en cauces y a volúmenes de control constantes en el espacio en los que sólo
se modifica el tirante hidráulico. El modelo matemático resuelve las ecuaciones de continuidad y
conservación de la cantidad de movimiento utilizando diferencias finitas en su forma implícita.
∂ (AV ) ∂Y
=B (4.28)
∂x ∂t
∂Y ⎡V ∂V 1 ∂V ⎤
= −⎢ + + Sf ⎥ (4.29)
∂x ⎣ g ∂x g ∂t ⎦
donde:
Y , elevación de la superficie del agua respecto a un plano horizontal de referencia
(suma del tirante más la cota de plantilla), en m.
A , área hidráulica, en m2.
V , velocidad media, en m/s.
B , ancho de la superficie libre, en m.
g , aceleración de la gravedad, en m/s2.
Sf , pendiente de la línea de energía, adim.
x , distancia longitudinal, en m.
t , tiempo, en s.
Para representar en diferencias finitas a las ecuaciones 4.28 y 4.29 se toman en cuenta los
volúmenes de control, según la figura 4.5.
Debido a que las elevaciones y velocidades varían con el tiempo y la posición, los distintos
términos de las ecuaciones 4.28 y 4.29 se escriben en diferencias finitas de la manera
siguiente:
41
Procesamiento hidráulico
(i-1)
(i)
Gi,j
(i+1)
Yi-1,j
Yi, j
Ui-1,j
Vi, j Yi+1,j
Ui, j
Vi+1,,j
PHC Δx Δx
Figura 4.5 Secciones a lo largo del cauce en diferencias finitas
∂Y
∂x
≅
θ
Δx i
(
Yi +1, j +1 − Yi , j +1 +
1− θ
Y
Δx i i +1, j
)
− Yi , j ( ) (4.30)
V ∂V Vi + 1, j + 1 + U i , j + 1 ⎛
≅ ⎜Vi + 1, j − U i , j ⎞⎟ (4.31)
g ∂t 2gΔx i ⎝ ⎠
1 ∂V
g ∂t
≅
1
2gΔt
[
(U i , j +1 + Vi +1, j +1) − (U i , j + Vi +1, j ) ] (4.32)
1 n2
Sf ≅ U i . j + Vi +1, j (U i , j + 1 + Vi + 1, j +1 ) (4.33)
4 ri4,j/ 3
∂( AV ) Ai + 1, j Vi + 1, j + 1 − Ai , j U i , j + 1
≅ (4.34)
∂x Δx i
∂Y Bi + 1, j + Bi , j
B
∂t
≅
4 Δx i
[
(Yi + 1, j + 1 + Yi , j + 1 ) − (Yi +1, j + Yi , j ) ] (4.35)
donde:
θ , es un factor de peso que sirve para calcular promedios ponderados en el tiempo con
la finalidad de mejorar la aproximación de las derivadas temporales.
Sustituyendo las ecuaciones 4.34 y 4.35 en la ecuación 4.28 y las ecuaciones 4.30 a la 4.33 en
la 4.29 se obtiene:
(
Ai , j U i , j + 1 − Ai + 1, j Vi +1, j + 1 = Fi , j Yi + 1, j + 1 + Yi , j + 1 − Yi + 1, j − Yi , j ) (4.36)
C i , j U i , j + 1 − C i , jVi + 1, j + 1 = Yi + 1, j + 1 − Yi , j + 1 + Di , j (4.37)
donde:
Fi , j =
Δx i
4Δt
(
Bi +1, j + Bi , j ) (4.38)
U i , j − Vi +1, j Δx i Δx n 2
Ci , j = − − i 4 / 3 U i , j + Vi +1, j (4.39)
2gθ 2gθΔt 4θ ri , j
42
Procesamiento hidráulico
1− θ Δx i
Di , j = (Yi + 1, j − Yi , j ) − (U i , j + Vi , j ) (4.40)
θ 2gθΔt
Las ecuaciones 4.36 y 4.37 forman un sistema lineal no homogéneo en U i , j +1 y Vi +1, j +1 cuya
solución se escribe como:
U i , j + 1 = Pi , j Yi +1, j + 1 + Qi , j Yi , j + 1 + R i , j (4.41)
Vi , j +1 = S i , j Yi +1, j +1 + Ti , j Yi , j +1 + W i , j (4.42)
En estas expresiones:
Fi , j Ai + 1
Pi , j = + (4.43)
Ai , j + Ai + 1, j C i , j ( Ai , j + Ai + 1, j )
Fi , j Ai + 1, J
Qi , j = + (4.44)
Ai , j + Ai + 1, j C i , j ( Ai , j + Ai + 1, j )
Di , j Ai + 1, j Fi , j (Yi + 1, j + Yi , j )
Ri , j = − (4.45)
C i , j + ( Ai , j + Ai + 1, j ) ( Ai , j + Ai + 1, j )
Ai , j Fi , J
Si , j = − (4.46)
C i , j ( Ai , j + Ai +1, j ) ( Ai , j + Ai +1, j )
Ai , j Fi , J
Ti , j = + (4.47)
C i , j Ai , j + Ai + 1, j Ai , j + Ai + 1, j
Di , j Ai + 1, j Fi , j (Yi + 1, j + Yi , j )
Wi , j = + (4.48)
C i , j + ( Ai , j + Ai +1, j ) ( Ai , j + Ai + 1, j )
1
Ti − j Yi , j −1 + (Qi , j − Si −1, j )Yi , j +1 + Pi . j Yi +1, j +1] = Wi −1, j − R i , j + (θGi , j +1 + (1 − θ )Gi , j (4.50)
Ai , j
La ecuación 4.49 junto con las de las secciones correspondientes a sus extremos, forman un
sistema de ecuaciones lineales tridiagonal, cuyas incógnitas son las elevaciones Yi en la etapa
de cálculo j + 1 . Una vez obtenidas las elevaciones, las velocidades de llegada Vi , j + 1 y de
partida U i , j + 1 , se calculan con las ecuaciones 4.41 y 4.42.
43
Procesamiento hidráulico
En las secciones inicial y final deben fijarse las condiciones de frontera, las cuales se
determinan de acuerdo con el problema en estudio. Por ejemplo, para el caso particular de flujo
subcrítico la condición de frontera aguas arriba se obtiene de la manera siguiente:
Como condición de frontera aguas abajo, se considera una elevación de la superficie libre del
agua conocida, sea Yf 1 , las velocidades del último tramo son las siguientes:
3TM−2, j YM−2, j +1 + (QM −1, j − SMi −2, j )YMi −1, j = WM −2, j − RM −1, j −
1
AM −1
[θGM−1, j +1 + (1− θ )GM−1, j ] − PM−1, jYF (4.54)
Para asegurar que la descarga de esta sección sea con una elevación del agua mayor o igual a
la mínima, la asociada al tirante crítico y cj para el instante j se emplea la ecuación siguiente:
Acj VMj
= (4.55)
Bcj g
donde:
Acj , área correspondiente al tirante crítico, en m2.
Bcj , ancho de la superficie libre correspondiente al tirante crítico, en m.
Vmj , velocidad de entrada en la sección M para el instante j , en m/s.
Para obtener las elevaciones del agua en la etapa de cálculo se forma un sistema de
ecuaciones lineales, en la cual se deben considerar los siguientes puntos:
1. Sección aguas arriba del tramo en estudio i = 1 , se debe utilizar la ecuación de condición
de frontera aguas arriba, ecuación 4.51.
2. Para las secciones intermedias i = 1,2,..., L − 1 , se utiliza la ecuación 4.50.
44
Procesamiento hidráulico
45
Procesamiento hidráulico
46
Sistema de Alerta Hidrometeorológica
5 SISTEMA DE ALERTA
HIDROMETEOROLÓGICA
Para que se presenten los escurrimientos máximos en un sitio, en respuesta a una
precipitación, deben transcurrir varios minutos u horas después de iniciada la tormenta (tiempo
de retraso), por lo que es posible anticipar la ocurrencia de los mayores flujos de agua para
advertir el peligro por desbordamientos de los ríos o por la intensidad de sus corrientes.
El tiempo de retraso es el tiempo transcurrido entre le centro geométrico del hietograma efectivo
y el instante en que se produce el gasto de pico del hidrograma de escurrimiento directo (figura
5.1) y se evalúa en función del tiempo de concentración (t r = 0.6t c ) . Esta idea es la base
teórica de un Sistema de Alerta Hidrometeorológica.
hp (mm)
tr
t (h)
Q(m³/s)
t(h)
Figura 5.1 Ilustración del tiempo de retraso
Para definir áreas relativamente pequeñas (menores a 50 km2) expuestas a sufrir daños por
inundaciones, como es el área de un municipio o una ciudad, se puede tomar en cuenta la
medición de la precipitación y el nivel de la superficie libre del agua en los ríos que están
ocurriendo en la zona de estudio, debido a ello, es mayor la confiabilidad en los avisos de
alertamiento. A un sistema de este tipo se le llamará sistema de alerta hidrometeorológica
temprana municipal (Fuentes et al, 2002).
47
Sistema de Alerta Hidrometeorológica
Un SAHTM se justifica sobre todo en las regiones donde ocurren crecidas rápidas y severas
cercanas a las ciudades, donde los flujos de agua sobre el terreno tienen una gran masa y
velocidad, los cuales pueden provocar daños materiales y en ocasiones, la pérdida de vidas
humanas.
El tiempo de anticipación con que se cuenta para estimar las intensidades de precipitación y
escurrimientos es de entre 20 a 180 minutos, mediante un sistema de medición remota
instalado en las cuencas de la zona de estudio.
5.2.1 Componentes
Un SAHTM está formado por varias partes que se complementan entre sí (figura 5.2).
SUBSISTEMA DE
INFORMACIÓN A LA
POBLACIÓN
SUBSISTEMA DE MEDICIÓN
Y PROCESAMIENTO HIDRÁULICO
El subsistema de información incluye los procedimientos para enterar a la población sobre las
posibles zonas de peligro, la ubicación de los albergues en los lugares seguros, las rutas hacia
estos sitios desde sus viviendas, así como la puesta en marcha del plan operativo contra
inundaciones y el momento oportuno de iniciar una evacuación. Las zonas o áreas de peligro
por inundación se determinan realizando simulaciones de flujo bidimensional en planicies de
inundación considerando diferentes escenarios.
48
Sistema de Alerta Hidrometeorológica
Es necesario destacar que los componentes que integran el SAHTM deben estar elaborados
antes de la ocurrencia de un fenómeno hidrometeorológico importante, tomando en cuenta los
tipos de peligro en las zonas de la región donde pueden ocurrir daños (figura 5.3).
Vigilancia del
Reconocimiento del Ocurrencia del
fenómeno. Medición
fenómeno y de zonas fenómeno. Avenidas e
de lluvia o niveles de
de peligro inundaciones
agua en los ríos
Detección de
Estudio hidrológico precipitaciones y Daños
niveles críticos
Decisión de poner en
Elaboración de un
marcha el plan Aplicación del plan de
plan operativo contra
operativo contra reconstrucción
inundaciones
inundaciones
El subsistema de medición y procesamiento hidráulico estima los efectos que producirá la lluvia
en una región, en los minutos u horas posteriores a su ocurrencia, para advertir del peligro que
podría generarse en algunas zonas de una ciudad cercanas a los cauces de los ríos. Su
objetivo es predecir, con base en mediciones de precipitación o nivel de la superficie libre del
agua en los ríos, la ocurrencia de los eventos que causan daños a la población, para que se
actúe anticipadamente con la intención de disminuir sus consecuencias.
Para calcular los escurrimientos de los ríos en cada una de las cuencas de la zona de estudio
durante una tormenta se hace bajo el procedimiento siguiente:
1. Los hidrogramas unitarios se pueden obtener por diversos métodos, en este caso se
consideraron los métodos del HUIG y HUA, los cuales se obtienen a partir de las
características fisiográficas y geomorfológicas de las cuencas de estudio.
2. Se estiman las pérdidas de precipitación utilizando el concepto del coeficiente de
escurrimiento; para determinar el valor de este coeficiente existen valores teóricos de
49
Sistema de Alerta Hidrometeorológica
acuerdo con características del tipo y uso de suelo y vegetación; es importante mencionar
que este coeficiente varía de acuerdo con la condición de humedad antecedente a la
tormenta que se esté analizando.
3. La precipitación media se obtiene con el método de los polígonos de Thiessen; se
determinan factores de peso definidos como el cociente del área de influencia de cada
estación, entre el área total de la cuenca. Para este caso es necesario determinar varias
configuraciones para prever que una o más estaciones de medición de la precipitación no
puedan transmitir las alturas de precipitación medidas.
4. Se calcula la precipitación media efectiva multiplicando la precipitación media por el
coeficiente de escurrimiento; dado que la información de medición de la precipitación es
recibida cada 10 minutos; con el paso del tiempo se va formando el hietograma de
precipitación efectiva.
5. Se obtiene el hidrograma de escurrimiento directo con la convolución del hidrograma
unitario y el hietograma de precipitación efectiva.
6. En caso de ser necesario, se transita el hidrograma de escurrimiento directo desde el
punto donde se determinó hasta el sitio de interés.
7. Finalmente, los escurrimientos determinados se comparan con umbrales preestablecidos,
para que, en caso de ser necesario, se activen las señales de alerta.
Los umbrales son valores críticos ya sea de precipitación, gasto o nivel de la superficie libre del
agua en un río, que nos indican a partir de que magnitudes se inicia el probable
desbordamiento de un río. Para obtenerlos es necesario realizar el estudio hidrológico de la
zona de estudio y un estudio hidráulico del flujo dentro de los cauces, siempre y cuando se
cuente con información topobatimétrica de los mismos.
Determinar un umbral de precipitación único en una cuenca en particular para advertir de los
efectos que puede desencadenar un fenómeno hidrometeorológico, es difícil, dadas las
múltiples circunstancias que se pueden presentar en el lugar afectado (asentamientos
irregulares, escasa infraestructura hidráulica, etc.).
Por otra parte, la respuesta hidrológica de una cuenca a la precipitación, depende, entre otros
factores, del estado de saturación inicial del suelo, la pendiente del terreno, estado de los
cauces, adecuación de la infraestructura, así como también del uso y tipo de suelo; debido a
esto, las características de la precipitación que se deben tener en cuenta para definir un umbral
de precipitación son la intensidad, duración y extensión de la misma. Por lo tanto, para
determinar un umbral adecuado a las condiciones del sitio, es necesario tomar en cuenta las
consideraciones mencionadas anteriormente, así como también de un análisis estadístico de
los registros históricos de precipitación y su origen para proponer una lamina de precipitación
que se presenta con periodo de retorno entre dos y cinco años.
Para estimar el gasto o nivel de la superficie libre del agua máximos que puede conducir un río
antes de desbordarse, existe más de un procedimiento; en una primera aproximación se
pueden utilizar si es posible, alguno de los métodos mencionados en el inciso 2.4.2.
Un segundo procedimiento con buena aproximación consiste en utilizar los hidrogramas que se
obtienen para cada subcuenca al utilizar un modelo lluvia-escurrimiento y la lamina de
precipitación propuesta como umbral o los correspondientes a periodos de retorno mayores; a
partir del punto donde se determinó el hidrograma se transita por el cauce con un método
hidrológico hasta el sitio de interés, y de acuerdo a la geometría de las secciones del sitio se
puede inferir la capacidad hidráulica máxima del cauce. Es importante destacar que cuando la
50
Sistema de Alerta Hidrometeorológica
zona de estudio esta integrada por varias subcuencas se deben considerar las entradas de los
hidrogramas en los puntos correspondientes. Los procedimientos anteriores son adecuados
cuando no existe información detallada de los cauces.
Para reducir o restaurar los daños en una región afectada por un fenómeno natural
hidrometeorológico, se pueden poner en práctica medidas estructurales y no estructurales
(institucionales).
Las medidas estructurales se refieren a obras planeadas y diseñadas cuidadosamente para ser
construidas en los cauces o en la cuenca para derivar, confinar, retener o almacenar los
escurrimientos para disminuir los daños materiales y resguardar a los pobladores. Usualmente
son construidas por dependencias gubernamentales, ya que requieren de fuertes inversiones
para su realización y conservación. Entre estas acciones se pueden citar:
51
Sistema de Alerta Hidrometeorológica
Las medidas no estructurales, son todas aquellas que no están relacionadas con la
construcción directa de obras, y que permiten avisar a tiempo a la población que puede sufrir
una inundación, así como el control y manejo de los escurrimientos, cuando ello es factible,
para minimizar los daños, y consisten principalmente de:
Desde el punto de vista económico, tanto las medidas estructurales como las no estructurales
tienen aplicación en las zonas desarrolladas y densamente pobladas, mientras que en áreas
poco desarrolladas solamente las segundas se pueden justificar.
En resumen, se puede mencionar que las acciones estructurales tienden a reducir o eliminar los
daños por inundaciones en áreas más o menos extensas, mientras que las no estructurales
permiten disminuir o evitar la pérdida de vidas humanas y disminuir los daños materiales.
Todas las acciones que sirven para reducir o evitar los daños causado por fenómenos naturales
que pueden ocasionar un desastre, así como las acciones que permiten auxiliar a la población y
proporcionarle ayuda oportuna durante la presencia de esos fenómenos reciben el nombre de
“medidas para la prevención de desastres”.
52
Caso de aplicación
6 CASO DE APLICACIÓN
El crecimiento irregular de las zonas urbanas ha tenido una notable expansión en los últimos 20
años en todo el país con los cambios que ello implica; Tuxtla Gutiérrez, Chis., no es la
excepción y en los últimos 10 años ha tenido un crecimiento demográfico medio anual del 4 %.
Por la necesidad de ampliar la superficie de suelo para vivienda y la construcción de servicios
básicos, la población se ha establecido en las márgenes del cauce del río Sabinal y sus
afluentes (zonas federales) que son zonas de peligro por inundación, lo que además genera
una degradación de los ecosistemas por la acelerada deforestación de la cuenca y la
impermeabilización de la superficie y las consecuencias que ello conlleva. Aunado a ello, se
combinan fenómenos hidrometeorológicos tales como lluvias intensas y ciclones tropicales
provocando que el comportamiento del río Sabinal y algunos de sus afluentes durante la época
de lluvia se torne peligroso habiéndose generado en los últimos 10 años inundaciones
recurrentes que han traído consigo la pérdida de de vidas humanas y cuantiosos daños
materiales.
Figura 6.1 Inundación del 04 al 06 de octubre de 2003. El río Sabinal se desbordó cerca del
estrechamiento del puente ubicado en las calles 3ª Ote. y 7ª Nte.
53
Caso de aplicación
Entre los días 04 y 06 de octubre de 2003, la tormenta tropical Larry causó lluvias intensas en
el centro y sur de la República Mexicana, en especial en la parte alta de la cuenca del río
Sabinal (225 mm), generando grandes gastos durante su recorrido por los municipios de
Berriozábal y Tuxtla Gutiérrez en donde la capacidad del cauce fue insuficiente a la avenida
estimada en 300 m3/s, afectando una superficie aproximada de 290 hectáreas, siendo las
colonias FOVISSSTE, Los Laureles, Jardines de Tuxtla y San Francisco Sabinal las más
afectadas.
Como puede observarse, la cuenca del río Sabinal se ve frecuentemente afectada por
inundaciones, además de que por lo menos 2,000 viviendas se encuentran construidas en las
márgenes de los cauces (GRFS CNA, 2003).
6.1.1 Localización
ZONA EN
ESTUDIO
La denominada cuenca del río Sabinal abarca porciones de cinco municipios, entre ellos la
capital del estado de Chiapas, Tuxtla Gutiérrez, que aporta casi el 40 % del área (141.61 km2)
de 357.76 km2, que es el área de aportación de dicha cuenca.
54
Caso de aplicación
6.1.2 Población
La cuenca cuenta con tres manchas urbanas importantes: Tuxtla Gutiérrez con un área de
78.08 km2, Berriozábal con 2.48 km2 y San Fernando con 1.98 km2, las cuales generan un área
total urbanizada de 82.54 km2 que representa el 23.07 % del área total de la misma, la mayor
parte de ésta corresponde a Tuxtla Gutiérrez (94.60 %) ubicada en la parte baja (figura 6.3).
Zona
urbana
Los municipios que abarca la cuenca en estudio tienen una población total de 712, 683
habitantes de los cuales 539, 380 (tabla 6.2) habitan dentro de la misma, según el II Conteo de
población y vivienda 2005 (INEGI, 2006), la mayor parte de ésta (91.38 %) se encuentra
asentada en su parte baja (figuras 6.3 y 6.4).
N #
#
#
# #
# #
## #
# #
# #
W E # #
# #
#
#
#
# #
S # #
#
#
# #
#
## # # #
# #
SAN
# FERNANDO
# #
#
# #
#
#
# #
# #
## #
## #
#
# #
# #
## # #
# ## #
## # #
##
#
# ## # # #
BERRIOZABAL
# # #
# ######
#
# # #
## #
##
#
## #
# # ##
#
## # # ### #
## #
## # #
# # #
#
# #
# ## #
#
# ####
# #
# ## ### #
## CHIAPA DE CORZO
# ## # # ## # #
# # # #
## # # # # #
# #
# #
# # #
#
## ## #
## ### # #
# # #
# ## # #
#
# # ## # #
#
# # #
# TUXTLA GUTIÉRREZ
#
#
#
# # #
#
#
OCOZOCOAUTLA
# DE ESPINOSA
#
55
Caso de aplicación
6.1.3 Precipitación
250.0
200.0
150.0
hp (mm)
100.0
50.0
0.0
Ene Feb Mar Abr May Jun Jul Ago Sep Oct Nov Dic
6.1.4 Hidrografía
El río Sabinal, llamado así porque a finales del siglo pasado y principios del actual existían en
sus márgenes árboles de Sabino, es un afluente del río Grijalva. Nace en la loma El Chupadero,
unos 5 km al noroeste de la población de Berriozábal, Chiapas, a una altitud del orden de 1,100
msnm. Corre en dirección sureste atravesando dicha población; adelante recibe por la margen
izquierda a la altura del poblado río Dulce una corriente que drena la sierra La Tienda. Aguas
abajo cambia su curso al oriente y posteriormente al sur para recibir la aportación del arroyo
San Francisco y sigue su curso al oriente para cruzar por la población de Terán. Un poco antes
de dejar ésta localidad (cruce con el Boulevard Belisario Domínguez) se incorpora al río Sabinal
por la margen derecha, el caudal de un pequeño escurrimiento denominado arroyo La Laguna,
que proviene de las faldas del cerro Mactumatzá.
Aguas abajo después de cruzar el boulevard Belisario Domínguez, recibe por la margen
izquierda la corriente conocida como arroyo San Agustín, el que a su vez recibe por su margen
izquierda dos subafluentes que se originan en las inmediaciones del poblado San Fernando, a
saber, los arroyos El Ojo y La Chacona, los cuales se originan como manantiales; estas dos
corrientes drenan una cuenca que se caracteriza por sus condiciones geológicas y topográficas
que provocan que algunos escurrimientos superficiales tiendan a fluir en “resumideros”.
56
Caso de aplicación
El río Sabinal sigue su curso al oriente atravesando el valle de Tuxtla Gutiérrez, donde a la
altura de la 23ª calle poniente norte recibe por la margen derecha al arroyo El Cocal y en la 16ª
calle poniente norte capta una pequeña corriente que actualmente se encuentra embovedada;
por la margen izquierda pero ahora entre las calles 13ª y 12ª poniente norte se le incorpora el
arroyo Lupanaspak o Pomarrosa o Mirador; posteriormente sobre su margen derecha, cerca de
la 4ª calle poniente norte se encuentra la confluencia del arroyo San Francisco, también
embovedado; aguas debajo de esa confluencia, aproximadamente en la 1ª calle oriente norte
sobre la margen izquierda toma la aportación del arroyo Potinaspak.
Recoge aguas sobre el arroyo San Roque o El Zope por la margen derecha, entubada en su
parte baja, al igual se encuentra otra corriente que se le suma sobre la calle central norte sobre
la misma margen; en la 15ª calle oriente norte, recibe por la margen izquierda las aportaciones
del arroyo Potinaspak o el Poti; que se encuentra entubada en su parte baja; continuando hacia
aguas abajo aproximadamente 500 m antes de llegar al cruce con el libramiento Norte de la
ciudad, recibe por su margen derecha un afluente que baja de cerro Hueco, al cual se le
agrega otra corriente llamada el Zapotal o Santa Ana, posteriormente sale de la mancha urbana
y recibe por ambas márgenes una serie de arroyos que descienden de los cerros Nido de
Águilas y la Ceiba, para descargar finalmente, en el río Grijalva (Bañuelos, 1997).
Para obtener las características fisiográficas de la cuenca del río Sabinal se utilizó la
información topográfica del Instituto Nacional de Estadística Geografía e Informática, que en
este caso fueron las cartas topográficas impresas escala 1:50,000 y los Modelos Digitales de
Elevación con resolución de 50 m, Chicoasén (E15C59) y Tuxtla Gutiérrez (E15C69). Se trabajó
con ambas fuentes de información y debido a su origen, con las cartas topográficas se utilizó la
metodología tradicional y con los Modelos Digitales de Elevación se utilizaron metodologías
utilizando Sistemas de Información Geográfica, obteniéndose dos grupos de información las
cuales fueron comparadas y validadas (Aragón, et al, 2004). En el presente trabajo se utilizará
una combinación de ambos resultados; es importante destacar que el uso de las metodologías
con Sistemas de Información Geográfica permite un ahorro significante de tiempo en horas-
hombre.
57
Caso de aplicación
SAN FERNANDO
(CHACONA)
EL SABINAL POTI
EL POTI
SAN AGUSTÍN
TOTOPOSTE
BERRIZOABAL
RÍO GRIJALVA
SA
NA 24 DE JUNIO
GU
EL SABINAL ST SABINAL
ÍN
CENTRO
EL SABINAL PATRIA
SAN FRANCISCO SAN FRANCISCO EL SABINAL
NUEVA
jito
L
SAN JOSÉ ZAPOTA A. GRANDE
gua
(C. HUECO)
SABIN
A. A
A. Grande
ZAPOTAL
(POC POC)
SABINO
ESCALA GRÁFICA (SAN ROQUE)
0 1.5 5.0 km
0.5 2.5
Para cada una de las subcuencas de orden tres y cuatro (San Francisco, Berriozábal, San
Agustín, San Fernando, El Poti y Patria Nueva) se realizó una discretización en microcuencas
de orden 1 hasta el orden de la mismas, determinando a partir de estas, el número, longitud y
área de corrientes; los resultados obtenidos se muestran en la tabla 6.4; con las figuras 6.7 se
ilustra el procedimiento que se realizó y en la figura 6.8 se muestran la integración de la
delimitación de las microcuencas para la subcuenca San Francisco. Con estos datos se
calcularon los parámetros geomorfológicos de Horton ( RB , RL y R A ), realizando un ajuste por
mínimos cuadrados entre el orden de las corrientes y el logaritmo del número, promedio de las
longitudes y áreas de las mismas, los cuales serán utilizados posteriormente en la aplicación del
HUIG; los resultados se presentan el la tabla 6.5.
58
Caso de aplicación
2 2 2
ORDEN No. L (km) A (km ) ORDEN No. L (Km) A (km ) ORDEN No. L (km) A (km )
SAN FRANCISCO BERRIOZABAL SAN AGUSTÍN
1 24 42.20 31.47 1 36 50.68 34.79 1 17 32.78 18.13
2 7 27.82 35.99 2 9 33.46 43.18 2 5 12.11 10.32
3 2 12.80 29.91 3 3 24.87 49.09 3 1 14.68 30.13
4 1 15.76 57.00 4 1 20.75 64.38
SAN FERNANDO EL POTI PATRIA NUEVA
1 39 58.30 43.13 1 20 26.66 15.67 1 12 19.40 9.36
2 6 21.62 23.78 2 4 10.75 9.92 2 3 12.56 12.80
3 1 21.00 75.32 3 1 13.84 27.07 3 1 9.30 15.44
INSTITUTO DE INGENIERÍA, UNAM INSTITUTO DE INGENIERÍA, UNAM INSTITUTO DE INGENIERÍA, UNAM INSTITUTO DE INGENIERÍA, UNAM
a) SISTEMA DE ALERTA HIDROMETEOROLÓGICA
TEMPRANA MUNICIPAL
b) SISTEMA DE ALERTA HIDROMETEOROLÓGICA
TEMPRANA MUNICIPAL
c) SISTEMA DE ALERTA HIDROMETEOROLÓGICA
TEMPRANA MUNICIPAL d) SISTEMA DE ALERTA HIDROMETEOROLÓGICA
TEMPRANA MUNICIPAL
01 SAN FRANCISCO
01 SAN FRANCISCO 01 SAN FRANCISCO 01 SAN FRANCISCO
Figura 6.7 Determinación del número, longitud y área de corrientes para orden 1, 2, 3 y 4 de la
subcuenca San Francisco
Para todas las subcuencas y con los datos de la tabla 6.3, se determinó el tiempo de
concentración con la formula de Kirpich y considerando una duración en exceso de diez
minutos se calculó el tiempo al pico, el gasto de pico y velocidad promedio a la salida de la
cuenca (tabla 6.6), datos necesarios para la aplicación del HUA en las subcuencas de orden 2
(San José, Sabinal Centro, El Sabino, Totoposte, 24 de Junio y Arroyo Grande), y el HUIG para
las subcuencas de orden 3 y 4 (San Francisco, Berriozábal, San Agustín, San Fernando, El Poti
y Patria Nueva), además del tiempo de retraso para tener en cuenta el tiempo de respuesta de
cada una de las subcuencas.
Tabla 6.6. Tiempo de concentración y de retraso, tiempo y gasto de pico y velocidad promedio
59
Caso de aplicación
Los hidrogramas unitarios se determinaron para un Δt = 10 minutos por el método del HUIG
para las subcuencas de orden tres y cuatro (San Francisco, Berriozábal, San Agustín, San
Fernando, El Poti y Patria Nueva), con los datos de las tablas 6.3 y 6.4; para las cuencas de
orden dos (San José, Sabinal Centro, El Sabino, Totoposte, El Zapotal, 24 de junio y Arroyo
Grande), se utilizó el método de HUA con la información de la tabla 6.4; los resultados
obtenidos se pueden apreciar en las figuras 6.9 y 6.10, respectivamente.
a) HIDROGRAMA UNITARIO INSTANTÁNEO GEOMORFOLÓGICO b) HIDROGRAMA UNITARIO INSTANTÁNEO GEOMORFOLÓGICO
SAN FRANCISCO BERRIOZABAL
4.0
4.0
3.5
3.5
3.0
3.0
2.5
q (m 3/s/mm)
q (m 3/s/mm)
2.5
2.0 2.0
1.5 1.5
1.0 1.0
0.5 0.5
0.0
0.0
0.0 5.0 10.0 15.0 20.0 25.0
0.0 5.0 10.0 15.0 20.0 25.0
Tiem po (h)
Tiem po (h)
4.0
4.0
3.5
3.5
3.0 3.0
2.5
q (m 3/s/mm)
q (m 3/s/mm)
2.5
2.0 2.0
1.5 1.5
1.0
1.0
0.5
0.5
0.0
0.0 0.0 5.0 10.0 15.0 20.0 25.0
0.0 5.0 10.0 15.0 20.0 25.0 Tiem po (h)
Tiem po (h)
60
Caso de aplicación
4.0 4.0
3.5 3.5
3.0 3.0
2.5 2.5
q (m 3/s/mm)
q (m 3/s/mm)
2.0 2.0
1.5 1.5
1.0 1.0
0.5 0.5
0.0 0.0
0.0 5.0 10.0 15.0 20.0 25.0 0.0 5.0 10.0 15.0 20.0 25.0
Tiem po (h) Tiem po (h)
4.0
3.5
3.0
q (m 3/s/mm)
2.5
2.0
1.5
1.0
0.5
0.0
0.0 5.0 10.0 15.0 20.0 25.0
Tiem po (h)
4.0 4.0
3.5 3.5
3.0 3.0
q (m 3/s/mm)
2.5
q (m 3/s/mm)
2.5
2.0 2.0
1.5 1.5
1.0 1.0
0.5 0.5
0.0 0.0
0.0 5.0 10.0 15.0 20.0 25.0 0.0 5.0 10.0 15.0 20.0 25.0
Tiem po (h) TIEMPO (h)
4.0 4.0
3.5 3.5
3.0 3.0
q (m 3/s/mm)
2.5
q (m 3/s/mm)
2.5
2.0 2.0
1.5 1.5
1.0 1.0
0.5 0.5
0.0 0.0
0.0 5.0 10.0 15.0 20.0 25.0 0.0 5.0 10.0 15.0 20.0 25.0
Tiem po (h) Tiem po (h)
61
Caso de aplicación
4.0 4.0
3.5 3.5
3.0 3.0
2.5
q (m 3/s/mm)
q (m 3/s/mm)
2.5
2.0 2.0
1.5 1.5
1.0 1.0
0.5 0.5
0.0 0.0
0.0 5.0 10.0 15.0 20.0 25.0 0.0 5.0 10.0 15.0 20.0 25.0
Tiem po (h) Tiem po (h)
De la revisión de las estaciones climatológicas existentes se determinó que sólo 5 (tabla 6.7 y
figura 6.11) tienen influencia en la cuenca en estudio, y cuentan con precipitación acumulada en
24 horas. Para el procesamiento de la información se obtuvo la precipitación máxima anual,
cuidando la calidad y cantidad de la misma.
62
Caso de aplicación
m OCOZOCUAUTLA TUXTLA GTZ. (CFE) TUXTLA GTZ. (DGE) SAN FERNANDO BERRIOZABAL
1 153.4 55.00 111.8 35.00 148.1 54.00 131 26.00 246 17.00
2 120 27.50 100.2 17.50 106.4 27.00 110 13.00 225.5 8.50
3 120 18.33 87.5 11.67 104.6 18.00 96.1 8.67 170 5.67
4 113 13.75 82 8.75 103.8 13.50 95 6.50 164 4.25
5 106.9 11.00 79 7.00 97 10.80 93.5 5.20 152 3.40
6 98 9.17 78 5.83 96 9.00 90.2 4.33 136 2.83
7 97 7.86 77.4 5.00 94.5 7.71 88 3.71 135 2.43
8 96 6.88 76.4 4.38 87.7 6.75 88 3.25 130.23 2.13
9 93.5 6.11 75.3 3.89 87.2 6.00 80 2.89 120.23 1.89
10 93.2 5.50 72 3.50 86.5 5.40 78.68 2.60 116.8 1.70
11 91.5 5.00 70.8 3.18 85.8 4.91 76 2.36 96.9 1.55
12 89.69 4.58 65.95 2.92 82 4.50 75.2 2.17 82 1.42
13 88 4.23 65.8 2.69 81.5 4.15 73.68 2.00 79 1.31
14 86.2 3.93 64.41 2.50 81 3.86 72 1.86 69.5 1.21
15 83 3.67 64.17 2.33 79.9 3.60 68.68 1.73 46.5 1.13
16 82.3 3.44 63.69 2.19 79.5 3.38 65 1.63 34 1.06
17 81 3.24 63.5 2.06 78.8 3.18 60.4 1.53
18 81 3.06 61.95 1.94 77.7 3.00 60 1.44
19 80 2.89 60.17 1.84 77.4 2.84 58 1.37
20 80 2.75 60 1.75 77 2.70 55 1.30
21 80 2.62 60 1.67 75.9 2.57 51.5 1.24
22 80 2.50 59 1.59 75.5 2.45 50 1.18
23 79.6 2.39 57.95 1.52 74.3 2.35 46 1.13
24 78.5 2.29 56 1.46 73.1 2.25 40 1.08
25 78.41 2.20 55 1.40 72.5 2.16 40 1.04
26 76 2.12 55 1.35 72 2.08
27 72.5 2.04 54.8 1.30 71.5 2.00
28 71.6 1.96 53.5 1.25 68.6 1.93
29 70.41 1.90 53.3 1.21 66.3 1.86
30 68 1.83 51.6 1.17 65.5 1.80
31 68 1.77 47.3 1.13 64 1.74
32 68 1.72 42.5 1.09 63.7 1.69
33 68 1.67 36.5 1.06 63.5 1.64
34 67 1.62 27.3 1.03 63.2 1.59
35 63 1.57 62.2 1.54
36 63 1.53 62 1.50
37 62 1.49 61.3 1.46
38 61.7 1.45 60.4 1.42
39 60.5 1.41 60.3 1.38
40 58 1.38 60 1.35
41 58 1.34 59 1.32
42 57.5 1.31 59 1.29
43 54.5 1.28 58.7 1.26
44 53.5 1.25 56.1 1.23
45 53.5 1.22 54.1 1.20
46 53 1.20 52.5 1.17
47 52.7 1.17 52.1 1.15
48 52.7 1.15 50.5 1.13
49 51.5 1.12 49.6 1.10
50 49 1.10 48.7 1.08
51 46 1.08 46.5 1.06
52 45 1.06 46.3 1.04
53 42.4 1.04 39.8 1.02
54 20.9 1.02
63
Caso de aplicación
W E
#
S #
SAN FERNANDO
#
BERRIOZABAL
Con los datos de la tabla 6.9 se calculó la precipitación media en cada subcuenca a través del
método de los polígonos de Thiessen (figura 6.11), y los resultados obtenidos se presentan en
la tabla 6.10.
64
Caso de aplicación
2
No. SUBCUENCA Pm2 (mm) Pm5 (mm) A (km ) FRA Pm2 (mm) Pm5 (mm) FRD Pm2 (mm) Pm5 (mm)
01 SAN FRANCISCO 125.45 167.87 57.00 0.79 99.33 132.92 0.73 72.14 96.54
02 BERRIOZABAL 134.61 181.13 64.38 0.78 105.61 142.11 0.77 81.42 109.56
03 SAN AGUSTÍN 134.24 180.38 30.13 0.83 111.52 149.86 0.73 81.00 108.85
04 SAN FERNANDO 87.07 109.98 75.32 0.78 67.51 85.27 0.77 52.05 65.74
05 SAN JOSÉ 112.15 145.80 17.75 0.86 96.96 126.06 0.73 70.43 91.56
06 SABINAL CENTRO 80.80 97.10 27.70 0.84 67.55 81.18 0.73 49.06 58.96
07 EL POTI 81.31 98.71 27.07 0.84 68.09 82.67 0.73 49.46 60.04
08 EL SABINO 80.80 97.10 7.18 0.93 74.79 89.88 0.65 48.61 58.43
09 TOTOPOSTE 74.99 91.47 11.54 0.89 66.98 81.70 0.65 43.54 53.10
10 EL ZAPOTAL 72.24 88.80 10.97 0.90 64.77 79.62 0.65 42.10 51.75
11 24 DE JUNIO 70.85 87.46 6.67 0.93 65.95 81.41 0.65 42.87 52.92
12 ARROYO GRANDE 70.85 87.46 6.63 0.93 65.98 81.45 0.65 42.89 52.94
13 PATRIA NUEVA 70.85 87.46 15.44 0.87 61.91 76.42 0.73 44.96 55.50
1
y = 0.288Ln(x) + 0.9889
R2 = 0.9953
0.8
y = 0.65x0.3921
0.6 R2 = 1
FRD
0.4
y = 0.65x0.6215
0.2 R2 = 1
0
0 1 1 2 2 3 3
Duración (h)
d =1 h d=2h d =3 h
Logarítmica (d =1 h) Potencial (d = 2 h) Potencial (d =3 h)
65
Caso de aplicación
Con las diferencias de los porcentajes de precipitación obtenidos en la tabla 6.12 y utilizando la
metodología del bloque alterno, se formaron los hietogramas base para 1, 2 y 3 horas de
duración (figura 6.13).
Por último es necesario mencionar que para determinar el porcentaje de la precipitación que
corresponde a las duraciones de 1, 2 y 3 horas del total de de la de 24 horas, se ajustó una
curva tipo logarítmica a los datos, donde se consideró un FRD = 0.65 , para una hora y un
FRD = 1 para 24 horas, obteniéndose un FRD = 0.73 para una duración de dos horas y un
FRD = 0.77 para 3 horas.
%
%
%
0.20
0.20 0.20
0.10 0.10
0.10
0.00 0.00
0.00
10
20
30
40
50
60
70
80
90
100
110
120
130
140
150
160
170
180
10
20
30
40
50
60
70
80
90
100
110
120
10
20
30
40
50
60
W E
66
Caso de aplicación
Para determinar el número de escurrimiento por subcuenca, con las de áreas de uso de suelo y
vegetación para cada una de ellas se realizó un promedio pesado; los resultados obtenidos se
muestran en la tabla 6.14.
Como podemos observar, los coeficientes de escurrimiento (figura 6.15) que se obtuvieron con
el método de los números de escurrimiento son menores comparados con los de la tabla 6.14,
pero hay que recordar que estos últimos son un poco conservadores por lo que se consideran
como adecuados.
Para considerar las condiciones iniciales de humedad del suelo, se hace una corrección al
coeficiente de escurrimiento, esto es, si existe una precipitación acumulada en los últimos 5
días antes del evento menor a 25 mm el coeficiente de escurrimiento se debe de incrementar
en un 5 %, si ésta es mayor a 25 mm se debe de incrementar en un 10 %; éstos valores son
aproximados y no deben de salirse del rango de los valores máximos según la tabla 2.1.
67
Caso de aplicación
W E
S #
#
0.10
0.27
0.08 0.17
# # 0.15
0.25
0.17
# 0.10
%$
a $
# 0.20 0.34 #
a
%
0.16
0.10 0.21
Bajo los criterios establecidos para la selección de sitios de medición, en gabinete se realizó
una configuración preliminar de la ubicación de las estaciones pluviométricas e hidrométricas y
los puestos centrales de registro, tomando en cuenta que de éstos últimos, la ubicación se fijó
de acuerdo a las necesidades de las autoridades que operarán el Sistema de Alerta
Hidrometeorológica Temprana Municipal, Protección Civil Municipal para el caso del PCR-01 y
la Comisión Nacional del Agua, Gerencia Regional Frontera Sur en el PCR-02.
Figura 6.16 Ubicación espacial de las estaciones y PCR en la cuenca del río Sabinal
68
Caso de aplicación
Tabla 6.16 Datos generales de las estaciones y PCR en la cuenca del río Sabinal
COORD. GEOGRÁFICAS COORD. UTM, NAD 27, ZONA 15
No. NOMBRE TIPO
LAT. LONG. x y
AB-01 SAN ANTONIO BOMBANÓ LLUVIA 16° 45’ 12.5” 93° 18’ 31.4” 467,095.354 1,852,196.211
SO-02 SOLIDARIDAD LLUVIA 16° 46’ 07.9” 93° 15’ 35.9” 472,293.552 1,853,891.056
BE-03 BERRIOZABAL LLUVIA 16° 48’ 01.1” 93° 16’ 24.1” 470,871.419 1,857,371.284
CA-04 LA CARIDAD LLUVIA 16° 48’ 05.5” 93° 14’ 14.1” 474,719.498 1,857,501.527
VH-05 VISTA HERMOSA LLUVIA 16° 50’ 44.7” 93° 18’ 09.2” 467,768.202 1,862,402.792
VC-06 VIVA CÁRDENAS LLUVIA 16° 50’ 02.7” 93° 11’ 40.1” 479,281.296 1,861,097.778
MI-07 EL MIRADOR LLUVIA 16° 53’ 07.5” 93° 14’ 23.1” 474,464.357 1,866,781.475
OM-08 OBSEVATORIO METEOROLÓGICO LLUVIA 16° 45’ 47.0” 93° 08’ 50.8” 484,285.714 1,853,236.563
CM-09 CAMPESTRE NIVEL 16° 45’ 14.4” 93° 11’ 12.5” 480,089.746 1,852,238.395
JM-10 JOYYO MAYU NIVEL 16° 45’ 34.0” 93° 08’ 44.8” 484,463.051 1,852,836.983
PO-11 PARQUE ORIENTE NIVEL 16° 45’ 37.1” 93° 05’ 16.6” 490,626.957 1,852,928.610
PCR-01 PROTECCIÓN CIVIL PCR-01 16° 44’ 26.8” 93° 05’ 14.4” 490,691.139 1,850,768.490
PCR-02 CNA PCR-02 16° 45’ 36.0” 93° 09’ 00.0” 484,013.093 1,852,898.771
500.00 500.00
0.00
2,892.53
5,785.06
8,677.59
11,570.11
14,462.64
17,355.17
20,247.70
23,140.23
0.00
3,093.49
6,186.97
9,280.45
12,373.94
15,467.42
18,560.91
Distancia (m ) Distancia (m )
LÍNEA DE VISTA PERFIL DEL TERRENO LÍNEA DE VISTA PERFIL DEL TERRENO
1,000.00 900.00
900.00 800.00
800.00
700.00
700.00
600.00 600.00
500.00 500.00
0.00
1,899.15
3,798.30
5,697.44
7,596.59
9,495.74
11,394.89
13,294.03
15,193.18
17,092.33
18,991.48
20,890.63
0.00
1,693.44
3,386.88
5,080.31
6,773.75
8,467.19
10,160.63
11,854.06
13,547.50
15,240.94
16,934.38
Distancia (m ) Distancia (m )
LÍNEA DE VISTA PERFIL DEL TERRENO LÍNEA DE VISTA PERFIL DEL TERRENO
1,200.00 1,000.00
1,100.00
900.00
1,000.00
900.00 800.00
800.00 700.00
700.00
600.00 600.00
500.00 500.00
0.00
3,200.80
6,401.59
9,602.38
12,803.18
16,003.97
19,204.77
22,405.56
25,606.36
0.00
2,198.70
4,397.39
6,596.09
8,794.78
10,993.48
13,192.17
15,390.87
Distancia (m ) Distancia (m )
LÍNEA DE VISTA PERFIL DEL TERRENO LÍNEA DE VISTA PERFIL DEL TERRENO
69
Caso de aplicación
3,799.58
7,599.15
11,398.73
15,198.31
18,997.88
22,797.46
0.00
2,288.16
4,576.31
6,864.46
Distancia (m ) Distancia (m )
LÍNEA DE VISTA PERFIL DEL TERRENO LÍNEA DE VISTA PERFIL DEL TERRENO
1,500.39
3,000.79
4,501.18
6,001.58
7,501.97
9,002.36
10,502.76
0.00
1,093.77
2,187.53
3,281.30
4,375.07
5,468.84
6,562.60
Distancia (m ) Distancia (m )
LÍNEA DE VISTA PERFIL DEL TERRENO LÍNEA DE VISTA PERFIL DEL TERRENO
575.00 575.00
550.00 550.00
525.00 525.00
500.00 500.00
0.00
1,080.54
2,161.07
0.00
1,001.37
2,002.74
3,004.11
4,005.48
5,006.85
6,008.22
LÍNEA DE VISTA PERFIL DEL TERRENO LÍNEA DE VISTA PERFIL DEL TERRENO
Figura 6.17 Líneas de vista entre las estaciones y el PCR-02 con el PCR-01
Una vez definida la ubicación final, se construyeron las casetas para la protección de los
equipos; posteriormente personal del Instituto de Ingeniería, UNAM y del Centro Nacional de
Prevención de Desastres (CENAPRED) realizaron la instalación e instrumentación de las
estaciones pluviométricas, hidrométricas y puestos centrales de registro.
70
Caso de aplicación
b)
a)
c)
d)
f) g)
e)
71
Caso de aplicación
a) b) c)
d)
e)
a)
b)
72
Caso de aplicación
c) d)
Para determinar los factores de peso se tiene que calcular el área de influencia de cada
estación pluviométrica (tabla 6.16 y figura 6.16), para ello se utiliza el método de los polígonos
de Thiessen (figura 6.21a); los resultados aparecen en la tercera fila de la tabla 6.17.
N N
a) b)
W E W E
# #
# #
EL MIRADOR EL MIRADOR
S S
# #
VISTA HERMOSA VISTA HERMOSA
# #
VIVA CÁRDENAS VIVA CÁRDENAS
# # # #
BERRIOZABAL LA CARIDAD BERRIOZABAL LA CARIDAD
# #
SOLIDARIDAD # SOLIDARIDAD #
OBSEVATORIO METEOROLÓGICO OBSEVATORIO METEOROLÓGICO
# #
SAN ANTONIO BOMBANÓ SAN ANTONIO BOMBANÓ
73
Caso de aplicación
Se parte de suponer que no funciona una estación a la vez, tal como se muestra en la figura
6.22a, donde se tiene una distribución de áreas para cuando no funciona la estación San
Antonio Bombanó; de esta forma tendríamos 8 nuevas distribuciones de áreas de influencia.
De forma similar, ahora se supone que no funcionan dos estaciones a la vez, en este caso son
las estaciones San Antonio Bombanó y Solidaridad; en la figura 6.22b se muestra la distribución
de las áreas de influencia; con esta suposición se tendrían 28 distribuciones de áreas de
influencia. Continuando de esta forma, en las figuras 6.22c, 6.22d y 6.22e se presenta la
configuración de áreas de influencia, para cuando no funcionan 3, 4 y 5 estaciones
respectivamente, siendo estas, San Antonio Bombanó, Solidaridad, Vista Hermosa, Berriozábal
y La Caridad.
74
Caso de aplicación
N
a)
W E
#
EL MIRADOR
S #
#
VISTA HERMOSA
#
VIVA CÁRDENAS
# #
BERRIOZABAL LA CARIDAD
#
SOLIDARIDAD #
OBSEVATORIO METEOROLÓGICO
N N
b) c)
W E W E
# #
EL MIRADOR EL MIRADOR
S #
S #
#
VISTA HERMOSA
# #
VIVA CÁRDENAS VIVA CÁRDENAS
# # # #
# #
OBSEVATORIO METEOROLÓGICO OBSEVATORIO METEOROLÓGICO
N N
d) e)
W E W E
# #
EL MIRADOR EL MIRADOR
S #
S #
# #
VIVA CÁRDENAS VIVA CÁRDENAS
#
LA CARIDAD
# #
OBSEVATORIO METEOROLÓGICO OBSEVATORIO METEOROLÓGICO
75
Caso de aplicación
Figura 6.23 Malla de celdas de 500 m por lado ajustada al perímetro de las subcuencas
76
Caso de aplicación
A Umbrales de precipitación
Según las condiciones físicas y eventos suscitados en la cuenca, se consideró que un periodo
de retorno de 5 años es el más adecuado para determinar los umbrales de precipitación. Para
ello, con los datos de precipitación máxima anual acumulada en 24 horas asociada a un periodo
de retorno de 5 años (tabla 6.11) y los hietogramas base (figura 6.13), se determinaron los
hietogramas de precipitación para cada subcuenca (tabla 6.21) y con los coeficientes de
escurrimiento propuestos (tabla 6.14), los hietogramas de precipitación efectiva (tabla 6.22),
que corresponden a los umbrales de precipitación en ambos casos.
Al momento en que sean sobrepasados los valores de los umbrales asociados a un periodo de
retorno de cinco años para cada subcuenca, se activará el color amarillo del semáforo de
alarmas por precipitación, con la finalidad de que las autoridades correspondientes tomen las
medidas adecuadas
La información para determinar la capacidad hidráulica del río Sabinal, es la geometría de las
de secciones, el valor de la n de Manning y los hidrogramas de escurrimiento directo.
77
Caso de aplicación
Se contó con la geometría de secciones a lo largo de más de 16 km del río Sabinal a una
distancia de 100 m en promedio (anexo A); la ubicación de las mismas como lo muestra la
figura 6.24 se encuentran dentro de la zona urbana, y en la figura 6.25 se presenta el perfil del
thalweg, bordo derecho y bordo izquierdo.
B.2 n de Manning
En el cálculo las pérdidas por fricción es indispensable contar con el valor del coeficiente de
Manning, el cual depende de la rugosidad en las paredes y fondo del cauce, geometría de las
78
Caso de aplicación
secciones y de las características del flujo, para evaluar dichas variables se realizó un recorrido
del mismo (figuras 6.26a a 6.26f) y con base en ello se propuso igual a 0.04.
610.00
600.00
590.00
580.00
ELEVACIÓN (m s n.m)
570.00
560.00
THALWEG
550.00
BOR. DER.
540.00
BOR. IZQ.
530.00
520.00
510.00
S = 0.00683
500.00
490.00
480.00
0+000 1+000 2+000 3+000 4+000 5+000 6+000 7+000 8+000 9+000 10+000 11+000 12+000 13+000 14+000 15+000 16+000
CADENAMIENTO (Km)
79
Caso de aplicación
Figura 6.26. Características y condiciones actuales del río Sabinal en diferentes secciones
Para obtener los hidrogramas de escurrimiento directo de cada subcuenca (figuras 6.27) se
realizó la convolución entre de los hietogramas asociados a un periodo de retorno de cinco
años (tabla 6.22) y los hidrogramas unitarios correspondientes (figuras 6.9 y 6.10).
80.0
70.0
60.0
50.0
Q (m3/s)
40.0
30.0
20.0
10.0
0.0
0.0 5.0 10.0 15.0 20.0 25.0 30.0
Tiem po (h)
80.0 80.0
70.0 70.0
60.0 60.0
50.0 50.0
Q (m3/s)
Q (m3/s)
40.0 40.0
30.0 30.0
20.0 20.0
10.0 10.0
0.0 0.0
0.0 5.0 10.0 15.0 20.0 25.0 30.0 0.0 5.0 10.0 15.0 20.0 25.0 30.0
Tiem po (h) Tiem po (h)
80
Caso de aplicación
80.0 80.0
70.0 70.0
60.0 60.0
50.0 50.0
Q (m3/s)
Q (m3/s)
40.0 40.0
30.0 30.0
20.0 20.0
10.0 10.0
0.0 0.0
0.0 5.0 10.0 15.0 20.0 25.0 30.0 0.0 5.0 10.0 15.0 20.0 25.0 30.0
Tiem po (h) Tiem po (h)
80.0 80.0
70.0 70.0
60.0 60.0
50.0 50.0
Q (m3/s)
Q (m3/s)
40.0 40.0
30.0 30.0
20.0 20.0
10.0 10.0
0.0 0.0
0.0 5.0 10.0 15.0 20.0 25.0 30.0 0.0 5.0 10.0 15.0 20.0 25.0 30.0
Tiem po (h) Tiem po (h)
80.0 80.0
70.0 70.0
60.0 60.0
50.0 50.0
Q (m3/s)
Q (m3/s)
40.0 40.0
30.0 30.0
20.0 20.0
10.0 10.0
0.0 0.0
0.0 5.0 10.0 15.0 20.0 25.0 30.0 0.0 5.0 10.0 15.0 20.0 25.0 30.0
Tiem po (h) Tiem po (h)
80.0 80.0
70.0 70.0
60.0 60.0
50.0 50.0
Q (m3/s)
Q (m3/s)
40.0 40.0
30.0 30.0
20.0 20.0
10.0 10.0
0.0 0.0
0.0 5.0 10.0 15.0 20.0 25.0 30.0 0.0 5.0 10.0 15.0 20.0 25.0 30.0
Tiem po (h) Tiem po (h)
81
Caso de aplicación
80.0 80.0
70.0 70.0
60.0 60.0
50.0 50.0
Q (m3/s)
Q (m3/s)
40.0 40.0
30.0 30.0
20.0 20.0
10.0 10.0
0.0 0.0
0.0 5.0 10.0 15.0 20.0 25.0 30.0 0.0 5.0 10.0 15.0 20.0 25.0 30.0
Tiem po (h) Tiem po (h)
De esta forma, para llevar a cabo la simulación del funcionamiento hidráulico se consideró que
por cada subcuenca (figura 6.6), ingresa un hidrograma al cauce del río según el esquema
simplificado de la figura 6.28; de acuerdo a ello el cauce del río quedó dividido en 10 tramos.
T-01 T-02 T-03 T-04 T-05 T-06 T-07 T-08 T-09 T-10
H-01 H-02 H-03 H-04 H-05 H-06 H-07 H-08 H-09 H-10
En la tabla 6.23 podemos ver la suma de hidrogramas que representan el comportamiento del
escurrimiento dentro del cauce, donde: h − i es el hidrograma correspondiente a la subcuenca
i según la figura 6.6 y tabla 6.3, y H − i , corresponde al hidrograma que se transita en el tramo
del cauce T − i .
Tabla 6.23 Comportamiento del escurrimiento dentro del cauce del río
H − 01 = h − 01 + h − 02
H − 02 = H − 01 + h − 05
H − 03 = H − 02 + h − 03 + h − 04 + h − 06
H − 04 = H − 03 + h − 07
H − 05 = H − 04 + h − 08
H − 06 = H − 05 + h − 09
H − 07 = H − 06 + h − 10
H − 08 = H − 07 + h − 11
H − 09 = H − 08 + h − 12
H − 13 = H − 09 + h − 13
B.4 Flujo permanente
Por lo tanto, para determinar la capacidad hidráulica del cauce en diferentes secciones (tramos)
en una primera estimación se utilizó el programa de simulación en flujo permanente del Instituto
de Ingeniería, UNAM, el cual utiliza de entrada la información descrita anteriormente. Se inició
proponiendo un gasto de 40 m3/s, con el cual no se desborda ninguna sección, a partir de éste
se incrementaron hasta que se tuvieran desbordamiento en prácticamente todas las secciones;
de esta forma se determinó la capacidad hidráulica de cada uno de los tramos del río Sabinal,
presentándose en la figura 6.29 algunos resultados de las simulaciones.
82
Caso de aplicación
a) b)
c) d)
e) f)
g) h)
83
Caso de aplicación
i) j)
Figura 6.29 Simulación en flujo permanente para las secciones con desbordamiento
Dado que en un tramo del cauce hay más de una sección, se consideró que la capacidad del
tramo es la correspondiente a la primer sección que comienza a desbordarse, posteriormente
de acuerdo con la geometría de la sección en cuestión se determinó el umbral del nivel de la
superficie libre del agua, para ello en la tabla 6.24 se presentan el tramo del río y la sección
crítica con su respectivos umbrales.
Tabla 6.24 Umbrales para cada tramo del río en flujo permanente
CADENAMIENTO UMBRALES
TRAMO
INI. FIN SECCIÓN Q (m3/s) Elev. Fondo SLA (msnm) y (m)
01 0+000 2+000 1+800 44.00 586.82 588.64 1.82
02 2+000 5+600 3+100 43.00 574.88 576.89 2.01
03 5+600 10+520 10+520 111.50 521.32 524.38 3.06
04 10+520 11+900 11+480 87.50 517.26 520.42 3.16
05 11+900 13+200 12+100 103.00 513.73 516.55 2.82
06 13+200 14+400 13+800 196.50 503.54 507.48 3.94
07 14+400 14+600 14+500 136.50 499.90 503.50 3.60
08 14+600 15+100 14+820 147.50 498.77 503.23 4.46
09 15+100 15+800 15+800 213.00 492.37 495.11 2.74
10 15+800 16+300 16+000 242.50 489.96 494.03 4.07
La variable tiempo, en el
tránsito de avenidas en cauces
es de gran importancia para
obtener resultados más
realistas, es por ello que, para
afinar los resultados
anteriores, se utilizó el
programa de simulación en
flujo no permanente del
Instituto de Ingeniería, UNAM
para transitar los hidrogramas
dentro del cauce (figura 6.30);
en la figura 6.31 se muestra la
superficie libre del agua, así
como el respectivo gasto,
ambos para el tiempo de pico
y en tabla 6.25 la capacidad
hidráulica de cada tramo.
Figura 6.30 Programa de simulación en flujo no permanente del II, UNAM
84
Caso de aplicación
Es necesario mencionar que en algunos tramos del río no se pudo determinar la capacidad
hidráulica con los hidrogramas asociados a un periodo de retorno de 5 años (figura 6.27), por lo
que fue necesario utilizar periodos de retorno mayores.
125
605
595 115
585 105
575
ELEVACIÓN (msnm)
95 THALWEG
565 BOR. DER.
GASTO (m3/s)
85
555 BOR. IZQ.
545 75 SLA
535 65 GASTO
525 55
515
45
505
495 35
485 25
0+000
1+000
2+000
3+000
4+000
5+000
6+000
7+000
8+000
9+000
10+000
11+000
12+000
13+000
14+000
15+000
16+000
CADENAMIENTO (Km)
Tabla 6.25 Umbrales para cada tramo del río en flujo no permanente
CADENAMIENTO UMBRALES
TRAMO 3
INI. FIN SECCIÓN Q (m /s) Elev. Fondo SLA (msnm) y (m)
01 0+000 2+000 2+000 40.35 583.93 586.23 2.30
02 2+000 5+600 3+100 62.63 574.88 576.89 2.01
03 5+600 10+520 7+700 84.67 541.06 545.26 4.20
04 10+520 11+900 11+480 77.51 517.26 520.42 3.16
05 11+900 13+200 12+100 107.81 513.73 516.55 2.82
06 13+200 14+400 14+400 199.02 500.60 505.22 4.62
07 14+400 14+600 14+500 64.28 499.90 503.50 3.60
08 14+600 15+100 14+820 57.48 498.77 503.23 4.46
09 15+100 15+800 15+800 190.26 492.37 495.11 2.74
10 15+800 16+300 16+000 139.12 489.96 494.03 4.07
85
Caso de aplicación
Tabla de precipitaciones
Tabla de niveles de la
superficie libre del
agua en el río
PA24, precipitación acumulada en las últimas 24 horas (a partir de las 08:00), en mm.
PA10, precipitación acumulada en los últimos 10 minutos, en mm.
PA60, precipitación acumulada en los últimos 60 minutos, en mm.
PAA1, precipitación acumulada en el intervalo t-10 minutos (en los últimos 20 minutos), en mm.
PAA2, precipitación acumulada en el intervalo t-20 minutos (en los últimos 30 minutos), en mm.
PAA3, precipitación acumulada en el intervalo t-30 minutos (en los últimos 40 minutos), en mm.
VB, voltaje de la batería, en volts.
STAT, código de estado.
ERAT, errores acumulados en la telemetría.
UITE, último intervalo con transmisión exitosa.
ESTADO, estado general de la telemetría.
Se muestran los niveles del agua registrados en cada una de las estaciones hidrométricas.
N10, nivel de la superficie libre del agua en los últimos 10 minutos, en msnm.
NA1, nivel en el intervalo t-10 minutos (en los últimos 20 minutos), en msnm.
NA2, nivel en el intervalo t-20 minutos (en los últimos 30 minutos), en msnm.
NA3, nivel en el intervalo t-30 minutos (en los últimos 40 minutos), en msnm.
VB, voltaje de la batería, en volts.
STAT, código de estado.
ERAT, errores acumulados en la telemetría.
UITE, último intervalo con transmisión exitosa.
ESTADO, estado general de la telemetría.
Una vez que se almacenan los datos en un archivo con formato de texto en la computadora
“Telemetría” son enviados a la computadora “Hidráulica” para que el programa de cómputo para
86
Caso de aplicación
Tabla de precipitaciones
Alarmas por
precipitación y
umbrales y gastos de
los tramos del cauce
Alarmas por
escurrimiento y
niveles de la
superficie libre del
agua en el río
Figura 6.33. Descripción de la pantalla del programa de cómputo para el procesamiento y análisis
de datos
HpA24, precipitación acumulada en las últimas 24 horas (a partir de las 08:00), en mm.
Hp10, precipitación acumulada en los últimos 10 minutos, en mm.
HpM60, precipitación acumulada en los últimos 60 minutos, en mm.
Hpi10, precipitación acumulada en el intervalo t-10 minutos (en los últimos 20 minutos), en mm.
Hpi20, precipitación acumulada en el intervalo t-20 minutos (en los últimos 30 minutos), en mm.
Hpi30, precipitación acumulada en el intervalo t-30 minutos (en los últimos 40 minutos), en mm.
Se muestran los gastos calculados para cada tramo del cauce del río Sabinal utilizando las
precipitaciones en intervalos de 10 minutos de la tabla de precipitaciones. Además, se muestran
los umbrales de gasto para cada tramo del río, valores determinados en la sección 6.2.7.
En el cuadro “ACTIVADA” está en verde cuando los gastos calculados están muy por debajo de
los umbrales de gasto, en amarillo cuando se encuentran muy cerca del umbral
(0.75Umbral ≤ Q < Umbral ) y en rojo una vez que estos son igualados o excedidos.
87
Caso de aplicación
a) HIETOGRAMAS DE PRECIPITACIÓN
14.0
12.0
Precipitación (mm)
10.0
AB-01
8.0
SO-02
6.0 BE-03
CA-04
4.0
2.0
0.0
13:40
13:50
14:00
14:10
14:20
14:30
14:40
14:50
15:00
15:10
15:20
15:30
15:40
15:50
16:00
16:10
16:20
16:30
16:40
16:50
17:00
17:10
17:20
17:30
17:40
17:50
18:00
18:10
18:20
18:30
18:40
18:50
19:00
19:10
19:20
19:30
19:40
19:50
Tiempo (h)
b) HIETOGRAMAS DE PRECIPITACIÓN
14.0
12.0
Precipitación (mm)
10.0
VH-05
8.0 VC-06
6.0 MI-07
OM-08
4.0
2.0
0.0
13:40
13:50
14:00
14:10
14:20
14:30
14:40
14:50
15:00
15:10
15:20
15:30
15:40
15:50
16:00
16:10
16:20
16:30
16:40
16:50
17:00
17:10
17:20
17:30
17:40
17:50
18:00
18:10
18:20
18:30
18:40
18:50
19:00
19:10
19:20
19:30
19:40
19:50
Tiempo (h)
88
Caso de aplicación
80.0
70.0 AB-01
Precipitación (mm)
60.0 SO-02
BE-03
50.0
CA-04
40.0
VH-05
30.0
VC-06
20.0 MI-07
10.0 OM-08
0.0
13:30
13:40
13:50
14:00
14:10
14:20
14:30
14:40
14:50
15:00
15:10
15:20
15:30
15:40
15:50
16:00
16:10
16:20
16:30
16:40
16:50
17:00
17:10
17:20
17:30
17:40
17:50
18:00
18:10
18:20
18:30
18:40
18:50
19:00
19:10
19:20
19:30
19:40
19:50
Tiempo (h)
Para cada intervalo de tiempo (10 minutos), una vez que la computadora “Hidráulica” recibe el
paquete de datos enviados por la computadora “Telemetría”, el programa de procesamiento y
análisis de datos, se encarga de calcular la precipitación media por subcuenca (figura 6.36),
además en la figura 6.37 se muestran las curvas masas de precipitación para cada subcuenca.
a) HIETOGRAMAS DE PRECIPITACIÓN
10.0
9.0
Precipitación (mm)
8.0
7.0
6.0 SAN FRANCISCO
5.0 BERRIOZABAL
4.0 SAN AGUSTÍN
3.0
2.0
1.0
0.0
13:40
14:00
14:20
14:40
15:00
15:20
15:40
16:00
16:20
16:40
17:00
17:20
17:40
18:00
18:20
18:40
19:00
19:20
19:40
Tiempo (h)
89
Caso de aplicación
b) HIETOGRAMAS DE PRECIPITACIÓN
8.0
7.0
Precipitación (mm)
6.0
5.0 SAN FERNANDO
4.0 SAN JOSÉ
3.0 EL POTI
2.0
1.0
0.0
13:40
13:50
14:00
14:10
14:20
14:30
14:40
14:50
15:00
15:10
15:20
15:30
15:40
15:50
16:00
16:10
16:20
16:30
16:40
16:50
17:00
17:10
17:20
17:30
17:40
17:50
18:00
18:10
18:20
18:30
18:40
18:50
19:00
19:10
19:20
19:30
19:40
19:50
Tiempo (h)
c) HIETOGRAMAS DE PRECIPITACIÓN
8.0
7.0 SABINAL CENTRO
Precipitación (mm)
6.0 EL SABINO
5.0 TOTOPOSTE
4.0 EL ZAPOTAL
3.0 24 DE JUNIO
2.0 ARROYO GRANDE
1.0 PATRIA NUEVA
0.0
13:40
13:50
14:00
14:10
14:20
14:30
14:40
14:50
15:00
15:10
15:20
15:30
15:40
15:50
16:00
16:10
16:20
16:30
16:40
16:50
17:00
17:10
17:20
17:30
17:40
17:50
18:00
18:10
18:20
18:30
18:40
18:50
19:00
19:10
19:20
19:30
19:40
19:50
Tiempo (h)
35.0
SAN FERNANDO
30.0
25.0 SAN JOSÉ
20.0 SABINAL CENTRO
15.0 EL POTI
10.0 EL SABINO
5.0
TOTOPOSTE
0.0
EL ZAPOTAL
13:30
13:40
13:50
14:00
14:10
14:20
14:30
14:40
14:50
15:00
15:10
15:20
15:30
15:40
15:50
16:00
16:10
16:20
16:30
16:40
16:50
17:00
17:10
17:20
17:30
17:40
17:50
18:00
18:10
18:20
18:30
18:40
18:50
19:00
19:10
19:20
19:30
19:40
19:50
24 DE JUNIO
ARROYO GRANDE
Tiempo (h)
PATRIA NUEVA
90
Caso de aplicación
ESTACIÓN AB-01 SO-02 BE-03 CA-04 VH-05 VC-06 MI-07 OM-08 Precipitación media (mm)
Precipitación (mm) P PTR=5 Pe
No. SUBCUENCA Ce
4.60 21.60 70.25 30.00 50.25 10.00 7.25 24.25 (mm) (mm) (mm)
01 SAN FRANCISCO 2.77 7.31 0.00 1.78 0 0 0 0 11.86 96.54 0.20 2.37
02 BERRIOZABAL 0 4.11 17.35 11.44 8.90 0.04 0 0 41.84 109.56 0.25 10.51
03 SAN AGUSTÍN 0 0 0 0 39.26 0 0.67 3.06 42.99 108.85 0.27 11.41
04 SAN FERNANDO 0 0 0 0 0.20 3.30 3.89 3.14 10.53 65.74 0.10 1.10
05 SAN JOSÉ 0 2.47 0 0 4.76 0 0 19.19 26.41 91.56 0.34 9.10
06 SABINAL CENTRO 0 0 0 0 0 0 0 24.25 24.25 58.96 0.17 4.01
07 EL POTI 0 0 0 0 0 0 2.71 15.18 17.89 60.04 0.08 1.44
08 EL SABINO 0 0 0 0 0 0 0 24.25 24.25 58.43 0.10 2.33
09 TOTOPOSTE 0 0 0 0 0 0 0 24.25 24.25 53.10 0.17 4.08
10 EL ZAPOTAL 0 0 0 0 0 0 0 24.25 24.25 51.75 0.16 3.96
11 24 DE JUNIO 0 0 0 0 0 0 0 24.25 24.25 52.92 0.15 3.52
12 ARROYO GRANDE 0 0 0 0 0 0 0 24.25 24.25 52.94 0.21 5.00
13 PATRIA NUEVA 0 0 0 0 0 0 0 24.25 24.25 55.50 0.10 2.54
Con los hietogramas de precipitación efectiva y los hidrogramas unitarios (figura 6.9 y 6.10
efectúa las convoluciones para obtener los hidrogramas de escurrimiento directo (figura 6.38).
30.0
20.0
Q (m3/s)
10.0
0.0
0.0 5.0 10.0 15.0 20.0 25.0 30.0
Tiem po (h)
30.0 30.0
20.0 20.0
Q (m3/s)
Q (m3/s)
10.0 10.0
0.0 0.0
0.0 5.0 10.0 15.0 20.0 25.0 30.0 0.0 5.0 10.0 15.0 20.0 25.0 30.0
Tiem po (h) Tiem po (h)
91
Caso de aplicación
30.0 30.0
20.0 20.0
Q (m3/s)
Q (m3/s)
10.0 10.0
0.0 0.0
0.0 5.0 10.0 15.0 20.0 25.0 30.0 0.0 5.0 10.0 15.0 20.0 25.0 30.0
Tiem po (h) Tiem po (h)
30.0 30.0
20.0 20.0
Q (m3/s)
0.0 0.0
0.0 5.0 10.0 15.0 20.0 25.0 30.0 0.0 5.0 10.0 15.0 20.0 25.0 30.0
Tiem po (h) Tiem po (h)
30.0 30.0
20.0 20.0
Q (m3/s)
Q (m3/s)
10.0 10.0
0.0 0.0
0.0 5.0 10.0 15.0 20.0 25.0 30.0 0.0 5.0 10.0 15.0 20.0 25.0 30.0
Tiem po (h) Tiem po (h)
30.0 30.0
20.0 20.0
Q (m3/s)
Q (m3/s)
10.0 10.0
0.0 0.0
0.0 5.0 10.0 15.0 20.0 25.0 30.0 0.0 5.0 10.0 15.0 20.0 25.0 30.0
Tiem po (h) Tiem po (h)
92
Caso de aplicación
30.0 30.0
20.0 20.0
Q (m3/s)
Q (m3/s)
10.0 10.0
0.0 0.0
0.0 5.0 10.0 15.0 20.0 25.0 30.0 0.0 5.0 10.0 15.0 20.0 25.0 30.0
Tiem po (h) Tiem po (h)
A continuación realiza el tránsito de los hidrogramas dentro del cauce considerando las
entradas en las diferentes secciones y realizando la suma de los hidrogramas según el
esquema idealizado de la figura 6.28.
Por último, si la capacidad en algunos de los tramos del cauce no es suficiente para conducir el
flujo producto de la precipitación, se activan las alarmas por escurrimiento.
Con los hidrogramas de escurrimiento directo producto de esta tormenta (figura 6.38), se llevó a
cabo el transito de los mismos dentro del cauce del río Sabinal con el programa de simulación
de flujo no permanente del Instituto de Ingeniería, UNAM, los resultados obtenidos nos indican
que no fue necesario que se activaran las alarmas por escurrimiento, dado que ningún umbral
del nivel de la superficie libre del agua fue superado según lo muestra la figura 6.39.
565 THALWEG
GASTO (m3/s) BOR. DER.
555 45
BOR. IZQ.
545
40 SLA
535
GASTO
525
35
515
505 30
495
485 25
0+000
1+000
2+000
3+000
4+000
5+000
6+000
7+000
8+000
9+000
10+000
11+000
12+000
13+000
14+000
15+000
16+000
CADENAMIENTO (Km)
93
Caso de aplicación
70.0
60.0
94
Conclusiones
7 CONCLUSIONES
Dentro de las acciones no estructurales o institucionales que se han implantado en México para
reducir los daños por inundación están los Sistemas de Alerta Hidrometeorológica Temprana
Municipal, en ciudades donde han ocurrido desastres de origen hidrometeorológico. Se han
instalado en Acapulco, Guerrero, Tijuana, Baja California, Motozintla y Tapachula, Chiapas, Villa
Hermosa, Tabasco, Álamo, Gutiérrez Zamora, Martínez de la Torre, Misantla, Poza Rica y
Tecolutla, Veracruz, Monterrey, Nuevo León, y Chalco, Iztapaluca y Los Reyes La paz, México.
Otro aspecto importante que se debe tener en cuenta, es que casi todas las obras de
ingeniería necesitan mantenimiento y tienen un periodo de vida útil, los SAHTM no son la
excepción y por lo tanto estos requieren de una supervisión y mantenimiento preventivo de por
lo menos una vez cada año y de preferencia antes de la temporada de lluvias para revisar el
funcionamiento adecuado del equipo electrónico, reparar fallas de comunicación por
interferencias o por descargas eléctricas, revisar el funcionamiento mecánico de los balancines
de los pluviómetros que son fallas comunes al atascarse o taparse por la acumulación de
basura o simplemente para corroborar el funcionamiento adecuado y la existencia y condiciones
del equipo, sin lo cual el SAHTM funcionaría de forma incorrecta, llegando incluso a dejar de
operar.
Por otro lado, no solamente desde el punto de vista de Protección Civil es importante el
SAHTM, sino que además, los datos almacenados de la medición tanto de precipitación como
del escurrimiento a intervalos de 10 minutos se han utilizado en estudios hidrológicos e
hidráulicos con la finalidad de determinar avenidas de diseño para la construcción de obras
hidráulicas con mayor certidumbre.
El subsistema de medición y procesamiento hidráulico del SAHTM emite los pronósticos con
base en un estudio hidrológico minucioso de la zona en estudio, como el presentado en este
trabajo en donde a partir del uso de suelo y la fisiografía de la cuenca se determinaron los
93
Conclusiones
También fue necesario realizar un estudio hidráulico de los cauces del río Sabinal para conocer
el comportamiento de la avenida al trasladarse de un sitio a otro aguas abajo por medio del
tránsito de avenidas; en este trabajo para determinar los umbrales del nivel de la superficie libre
del agua los cálculos se realizaron con un método de tipo hidráulico y el subsistema de
medición y procesamiento hidráulico los realiza con uno de tipo hidrológico, que aunque es
menos preciso en los resultados por tomar en cuenta menos variables, en estos casos son más
adecuados porque el tiempo para realizar los cálculos es mucho menor y esto favorece en la
toma de decisiones que tienen que hacer las autoridades, ya que en la cuenca del río Sabinal
los tiempos para que se presenten los escurrimientos máximos van desde los 30 a los 140
minutos.
Las estaciones hidrométricas consisten en un limnímetro de dos metros de longitud el cual toma
lecturas de la superficie libre del agua a distancias de 10 cm., cada 10 minutos. Aunque la idea
de instalar limnímetros para calibrar los modelos lluvia-escurrimiento y con esto aumentar la
confiabilidad de los pronósticos es valida, cuando se llegue a presentar una avenida, con las
lecturas de dichos instrumentos no se podrá reproducirla totalmente porque no cubren el rango
total de elevaciones de la superficie libre del agua, dado que la altura de los bordos del cauce
del río Sabinal en las secciones donde fueron instalados llegan a tener hasta 5 metros a partir
del thalweg, con lo cual se desaprovecharía una oportunidad de conocer y comprender la
hidrología de las subcuencas y comportamiento de los cauces de forma adecuada, ya que con
las mediciones tanto de precipitación como de escurrimiento se podrán estimar las pérdidas de
precipitación con mayor precisión, además de ajustar algunos parámetros de los hidrogramas
unitarios.
Un SAHTM se integra de otras partes tan importantes como son: el plan operativo contra
inundaciones, el plan para restaurar las condiciones previas al desastres y el subsistema de
información a la población, los cuales deben implantarse porque sin ellos, este tipo de
proyectos quedaría incompleto y no se podrían aprovechar todas su bondades, dentro de las
que se puede destacar están el conocer las zonas de peligro por inundación, a partir de las
cuales en una primera instancia serviría para regular el uso de suelo y evitar los asentamientos
humanos, y si éstos ya existiesen, se podrían realizar algunas acciones estructurales.
Como una forma de ejemplificar lo referente a las zonas con peligro por inundación
mencionadas en el párrafo anterior, Fuentes et al, 2006, realizó un estudio de flujo
bidimensional para la cuenca del río Sabinal en donde se concluye que para una tormenta con
un periodo de retorno de 10 años en condiciones actuales, el río Sabinal desborda 2.884 hm3,
inunda un área 2.653 km2 y afecta aproximadamente a 16,744 personas con sus respectivas
consecuencias. El 24.13 % del área inundada (0.64 km2), tiene una profundidad mayor a 0.60 m
con la cual ocurre la mayoría de las pérdidas materiales en viviendas tipo I, II y III según la
clasificación realizada por el CENAPRED; este tipo de viviendas son las predominantes en las
zonas con peligro por inundación ubicadas en las márgenes del río Sabinal, inclusive llegando a
alcanzar profundidades mayores de 5 m. en mucho menor proporción y que combinándose con
las altas velocidades que en algunos puntos alcanza más de 2.0 m/s, como lo sucedido en años
recientes, se podrían tener pérdida de vidas humanas; con todo esto, la implementación del
Sistema de Alerta Hidrometeorológica Temprana Municipal en Tuxtla Gutiérrez, se vuelve
fundamental para evitar o minimizar éstos daños.
94
Conclusiones
95
Conclusiones
96
Bibliografía
8 BIBLIOGRAFÍA Y REFERENCIAS
1. Aparicio, M. F. J. (2001). Fundamentos de hidrología de superficie. Editorial LIMUSA,
México, 300 pp.
10. Comisión Nacional del Agua (2005). Datos de las estaciones climatológicas del sistema
CLICOM.
11. Domínguez, M. R., Fuentes, M. O. A., De Luna, C. F., Cruz, G. J. A., López, E. J. G. y
Aragón, H. J. L. (2005). Estudio hidrológico, hidráulico y de control de inundaciones del río
Sabinal (Aspectos Hidrológicos). Instituto de Ingeniaría, UNAM, México, 51 pp.
13. Fuentes, M. O. A., Aragón, H. J. L., Cruz, G. J. J., De Luna, C. F. (2006). “Modelo numérico
de flujo bidimensional para la planeación de acciones en llanuras de inundación de zonas
urbanas”. XIX Congreso Nacional de Hidráulica, Cuernavaca, Morelos, México.
97
Bibliografía
14. Fuentes, M. O. A., Domínguez, M. R., De Luna, C. F., Cruz, G. J. A., Aragón, H. J. L. y
López, E. J. G. (2005). Estudio hidrológico, hidráulico y de control de inundaciones del río
Sabinal (Aspectos Hidráulicos). Instituto de Ingeniaría, UNAM, México.
15. Fuentes, M. O. A., Franco, S. M. A., González, P. J., Jiménez, E. M., Quass, W. R., Cruz,
G. J. A, Aragón, H. J. L., López, E. J. G., Ávila, L. A. A. y De Luna, C. F. (2006). Sistema de
Alerta Hidrometeorológica para la cuenca del río Sabinal en la ciudad de Tuxtla Gutiérrez,
Chiapas. Instituto de Ingeniería, UNAM, México.
17. Fuentes, M. O. A., Quaas, W. R., Jiménez, E. M., Franco, S. M. A., Eslava, M. H. y
González, P. J. (2002). Sistemas de Alerta Hidrometeorológica en Acapulco, Tijuana,
Motozintla, Tapachula y Monterrey. CENAPRED, México, 97 pp.
28. UNACH (2004). Topografía del cauce del río Sabinal en la zona urbana.
98
Anexo A
ANEXO A
Graficas de las secciones correspondientes al cauce del río Sabinal a lo largo de más de 16
kilómetros dentro de la zona urbana.
99
Anexo A
612.50 612.00
608.50 608.00
606.50 606.00
604.50 604.00
602.50 602.00
0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00 0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00
610.50 608.50
Elevación (m. s. n. m.)
606.50 604.50
604.50 602.50
602.50 600.50
600.50 598.50
0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00 0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00
607.50 606.50
Elevación (m. s. n. m.)
Elevación (m. s. n. m.)
605.50 604.50
603.50 602.50
601.50 600.50
599.50 598.50
597.50 596.50
0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00 0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00
606.00 605.50
Elevación (m. s. n. m.)
604.00 603.50
602.00 601.50
600.00 599.50
598.00 597.50
596.00 595.50
0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00 0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00
Distancia (m) Distancia (m)
604.50 603.00
Elevación (m. s. n. m.)
602.50 601.00
600.50 599.00
598.50 597.00
596.50 595.00
594.50 593.00
0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00 0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00
100
Anexo A
602.50 602.50
598.50 598.50
596.50 596.50
594.50 594.50
592.50 592.50
0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00 0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00
601.50 601.00
599.50 599.00
597.50 597.00
595.50 595.00
593.50 593.00
591.50 591.00
0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00 0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00
600.00 598.50
Elevación (m. s. n. m.)
598.00 596.50
596.00 594.50
594.00 592.50
592.00 590.50
590.00 588.50
0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00 0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00
598.00 597.50
Elevación (m. s. n. m.)
Elevación (m. s. n. m.)
596.00 595.50
594.00 593.50
592.00 591.50
590.00 589.50
588.00 587.50
0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00 0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00
596.50 595.50
Elevación (m. s. n. m.)
594.50 593.50
592.50 591.50
590.50 589.50
588.50 587.50
586.50 585.50
0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00 0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00
101
Anexo A
593.50 592.00
589.50 588.00
587.50 586.00
585.50 584.00
583.50 582.00
0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00 0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00
591.00 590.00
Elevación (m. s. n. m.)
587.00 586.00
585.00 584.00
583.00 582.00
581.00 580.00
0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00 0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00
589.50 589.00
Elevación (m. s. n. m.)
Elevación (m. s. n. m.)
587.50 587.00
585.50 585.00
583.50 583.00
581.50 581.00
579.50 579.00
0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00 0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00
588.00 588.00
Elevación (m. s. n. m.)
586.00 586.00
584.00 584.00
582.00 582.00
580.00 580.00
578.00 578.00
0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00 0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00
587.00 586.00
Elevación (m. s. n. m.)
585.00 584.00
583.00 582.00
581.00 580.00
579.00 578.00
577.00 576.00
0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00 0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00
102
Anexo A
585.00 584.50
581.00 580.50
579.00 578.50
577.00 576.50
575.00 574.50
0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00 0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00
583.50 583.50
Elevación (m. s. n. m.)
579.50 579.50
577.50 577.50
575.50 575.50
573.50 573.50
0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00 0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00
582.00 581.50
Elevación (m. s. n. m.)
Elevación (m. s. n. m.)
580.00 579.50
578.00 577.50
576.00 575.50
574.00 573.50
572.00 571.50
0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00 0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00
581.00 580.50
Elevación (m. s. n. m.)
579.00 578.50
577.00 576.50
575.00 574.50
573.00 572.50
571.00 570.50
0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00 0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00
579.50 579.00
Elevación (m. s. n. m.)
577.50 577.00
575.50 575.00
573.50 573.00
571.50 571.00
569.50 569.00
0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00 0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00
103
Anexo A
577.50 577.00
573.50 573.00
571.50 571.00
569.50 569.00
567.50 567.00
0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00 0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00
577.00 576.50
Elevación (m. s. n. m.)
573.00 572.50
571.00 570.50
569.00 568.50
567.00 566.50
0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00 0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00
576.00 575.50
Elevación (m. s. n. m.)
574.00 573.50
572.00 571.50
570.00 569.50
568.00 567.50
566.00 565.50
0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00 0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00
574.50 573.50
Elevación (m. s. n. m.)
572.50 571.50
570.50 569.50
568.50 567.50
566.50 565.50
564.50 563.50
0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00 0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00
572.50 572.50
Elevación (m. s. n. m.)
Elevación (m. s. n. m.)
570.50 570.50
568.50 568.50
566.50 566.50
564.50 564.50
562.50 562.50
0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00 0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00
Distancia (m) Distancia (m)
104
Anexo A
572.00 571.00
568.00 567.00
566.00 565.00
564.00 563.00
562.00 561.00
0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00 0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00
569.50 569.00
567.50 567.00
565.50 565.00
563.50 563.00
561.50 561.00
559.50 559.00
0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00 0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00
568.50 568.00
Elevación (m. s. n. m.)
Elevación (m. s. n. m.)
566.50 566.00
564.50 564.00
562.50 562.00
560.50 560.00
558.50 558.00
0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00 0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00
567.50 567.50
Elevación (m. s. n. m.)
565.50 565.50
563.50 563.50
561.50 561.50
559.50 559.50
557.50 557.50
0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00 0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00
567.00 565.50
Elevación (m. s. n. m.)
Elevación (m. s. n. m.)
565.00 563.50
563.00 561.50
561.00 559.50
559.00 557.50
557.00 555.50
0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00 0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00
105
Anexo A
565.50 565.50
561.50 561.50
559.50 559.50
557.50 557.50
555.50 555.50
0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00 0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00
564.50 563.00
Elevación (m. s. n. m.)
560.50 559.00
558.50 557.00
556.50 555.00
554.50 553.00
0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00 0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00
562.50 562.00
Elevación (m. s. n. m.)
560.50 560.00
558.50 558.00
556.50 556.00
554.50 554.00
552.50 552.00
0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00 0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00
Distancia (m) Distancia (m)
562.00 560.50
Elevación (m. s. n. m.)
560.00 558.50
558.00 556.50
556.00 554.50
554.00 552.50
552.00 550.50
0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00 0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00
558.50 558.50
Elevación (m. s. n. m.)
556.50 556.50
554.50 554.50
552.50 552.50
550.50 550.50
548.50 548.50
0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00 0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00
106
Anexo A
558.00 557.00
554.00 553.00
552.00 551.00
550.00 549.00
548.00 547.00
0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00 0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00
555.50 555.00
553.50 553.00
551.50 551.00
549.50 549.00
547.50 547.00
545.50 545.00
0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00 0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00
554.00 554.50
Elevación (m. s. n. m.)
552.00 552.50
550.00 550.50
548.00 548.50
546.00 546.50
544.00 544.50
0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00 0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00
552.50 551.00
Elevación (m. s. n. m.)
550.50 549.00
548.50 547.00
546.50 545.00
544.50 543.00
542.50 541.00
0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00 0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00
550.50 550.50
Elevación (m. s. n. m.)
548.50 548.50
546.50 546.50
544.50 544.50
542.50 542.50
540.50 540.50
0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00 0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00
107
Anexo A
548.50 544.00
544.50 540.00
542.50 538.00
540.50 536.00
538.50 534.00
0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00 0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00
543.50 543.50
541.50
541.50
539.50
539.50
537.50
537.50
535.50
535.50 533.50
533.50 531.50
0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00 0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00
540.50 540.00
Elevación (m. s. n. m.)
538.50 538.00
536.50 536.00
534.50 534.00
532.50 532.00
530.50 530.00
0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00 0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00
539.50 539.00
Elevación (m. s. n. m.)
537.50 537.00
535.50 535.00
533.50 533.00
531.50 531.00
529.50 529.00
0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00 0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00
538.50 538.00
Elevación (m. s. n. m.)
536.50 536.00
534.50 534.00
532.50 532.00
530.50 530.00
528.50 528.00
0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00 0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00
108
Anexo A
538.00 537.00
534.00 533.00
532.00 531.00
530.00 529.00
528.00 527.00
0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00 0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00
537.00 536.00
Elevación (m. s. n. m.)
533.00 532.00
531.00 530.00
529.00 528.00
527.00 526.00
0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00 0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00
535.00 534.00
Elevación (m. s. n. m.)
533.00 532.00
531.00 530.00
529.00 528.00
527.00 526.00
525.00 524.00
0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00 0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00
532.50 531.50
Elevación (m. s. n. m.)
530.50 529.50
528.50 527.50
526.50 525.50
524.50 523.50
522.50 521.50
0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00 0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00
533.00 533.00
Elevación (m. s. n. m.)
531.00 531.00
529.00 529.00
527.00 527.00
525.00 525.00
523.00 523.00
0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00 0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00
109
Anexo A
533.00 532.00
529.00 528.00
527.00 526.00
525.00 524.00
523.00 522.00
0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00 0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00
531.50 532.00
Elevación (m. s. n. m.)
527.50 528.00
525.50 526.00
523.50 524.00
521.50 522.00
0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00 0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00
531.00 531.00
Elevación (m. s. n. m.)
529.00 529.00
527.00 527.00
525.00 525.00
523.00 523.00
521.00 521.00
0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00 0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00
531.00 530.50
Elevación (m. s. n. m.)
Elevación (m. s. n. m.)
529.00 528.50
527.00 526.50
525.00 524.50
523.00 522.50
521.00 520.50
0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00 0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00
530.50 530.50
Elevación (m. s. n. m.)
Elevación (m. s. n. m.)
528.50 528.50
526.50 526.50
524.50 524.50
522.50 522.50
520.50 520.50
0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00 0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00
110
Anexo A
529.00 528.50
525.00 524.50
523.00 522.50
521.00 520.50
519.00 518.50
0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00 0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00
528.50 528.00
Elevación (m. s. n. m.)
524.50 524.00
522.50 522.00
520.50 520.00
518.50 518.00
0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00 0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00
527.50 527.00
Elevación (m. s. n. m.)
Elevación (m. s. n. m.)
525.50 525.00
523.50 523.00
521.50 521.00
519.50 519.00
517.50 517.00
0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00 0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00
526.50 526.00
Elevación (m. s. n. m.)
Elevación (m. s. n. m.)
524.50 524.00
522.50 522.00
520.50 520.00
518.50 518.00
516.50 516.00
0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00 0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00
524.50 524.00
Elevación (m. s. n. m.)
Elevación (m. s. n. m.)
522.50 522.00
520.50 520.00
518.50 518.00
516.50 516.00
514.50 514.00
0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00 0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00
111
Anexo A
523.50 523.00
519.50 519.00
517.50 517.00
515.50 515.00
513.50 513.00
0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00 0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00
522.00 521.50
Elevación (m. s. n. m.)
518.00 517.50
516.00 515.50
514.00 513.50
512.00 511.50
0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00 0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00
521.00 521.00
Elevación (m. s. n. m.)
519.00 519.00
517.00 517.00
515.00 515.00
513.00 513.00
511.00 511.00
0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00 0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00
520.00 519.00
Elevación (m. s. n. m.)
518.00 517.00
516.00 515.00
514.00 513.00
512.00 511.00
510.00 509.00
0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00 0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00
519.00 518.50
Elevación (m. s. n. m.)
517.00 516.50
515.00 514.50
513.00 512.50
511.00 510.50
509.00 508.50
0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00 0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00
112
Anexo A
517.50 517.50
513.50 513.50
511.50 511.50
509.50 509.50
507.50 507.50
0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00 0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00
516.50 516.50
Elevación (m. s. n. m.)
512.50 512.50
510.50 510.50
508.50 508.50
506.50 506.50
0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00 0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00
515.50 515.00
Elevación (m. s. n. m.)
513.50 513.00
511.50 511.00
509.50 509.00
507.50 507.00
505.50 505.00
0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00 0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00
514.00 513.00
Elevación (m. s. n. m.)
512.00 511.00
510.00 509.00
508.00 507.00
506.00 505.00
504.00 503.00
0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00 0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00
512.50 512.50
Elevación (m. s. n. m.)
510.50 510.50
508.50 508.50
506.50 506.50
504.50 504.50
502.50 502.50
0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00 0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00
113
Anexo A
512.00 511.50
508.00 507.50
506.00 505.50
504.00 503.50
502.00 501.50
0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00 0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00
511.00 510.00
Elevación (m. s. n. m.)
507.00 506.00
505.00 504.00
503.00 502.00
501.00 500.00
0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00 0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00
509.50 509.50
Elevación (m. s. n. m.)
507.50 507.50
505.50 505.50
503.50 503.50
501.50 501.50
499.50 499.50
0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00 0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00
508.50 508.50
Elevación (m. s. n. m.)
Elevación (m. s. n. m.)
506.50 506.50
504.50 504.50
502.50 502.50
500.50 500.50
498.50 498.50
0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00 0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00
508.00 507.00
Elevación (m. s. n. m.)
Elevación (m. s. n. m.)
506.00 505.00
504.00 503.00
502.00 501.00
500.00 499.00
498.00 497.00
0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00 0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00
114
Anexo A
505.00 506.00
501.00 502.00
499.00 500.00
497.00 498.00
495.00 496.00
0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00 0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00
505.50 505.00
Elevación (m. s. n. m.)
501.50 501.00
499.50 499.00
497.50 497.00
495.50 495.00
0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00 0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00
504.00 503.50
Elevación (m. s. n. m.)
502.00 501.50
500.00 499.50
498.00 497.50
496.00 495.50
494.00 493.50
0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00 0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00
503.00 502.00
Elevación (m. s. n. m.)
Elevación (m. s. n. m.)
501.00 500.00
499.00 498.00
497.00 496.00
495.00 494.00
493.00 492.00
0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00 0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00
501.50 499.50
Elevación (m. s. n. m.)
Elevación (m. s. n. m.)
499.50 497.50
497.50 495.50
495.50 493.50
493.50 491.50
491.50 489.50
0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00 0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00
115
Anexo A
500.50 499.50
496.50 495.50
494.50 493.50
492.50 491.50
490.50 489.50
488.50 487.50
0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00 0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00
CADENAMIENTO 16+300
497.50
Elevación (m. s. n. m.)
495.50
493.50
491.50
489.50
487.50
0.00 10.00 20.00 30.00 40.00 50.00 60.00 70.00 80.00
Distancia (m)
116