Modelado Hidráulico Amazónico HEC-RAS
Modelado Hidráulico Amazónico HEC-RAS
I. INTRODUCCIÓN
Los ríos navegables de la cuenca amazónica del Perú son el medio de comunicación más
importante después de la aérea, desde el punto de vista comercial y traslado de grandes
volúmenes, representa el principal medio de transporte, por tanto tiene influencia directa en
el desarrollo y bienestar de las ciudades y centros poblados ubicados en sus riberas. Es
necesario conocer las características hidráulicas, hidrológicas, topográficas, batimétricas,
ambientales del río y su entorno; para dar seguridad y agilizar la navegación, alertar de
eventos extraordinarios, como el periodo de vaciante extrema y las épocas de avenidas
máximas.
El programa HEC-RAS será nuestra herramienta para realizar las simulaciones del flujo en
estado permanente, se hace una breve descripción de otras simulaciones que es factible
realizar, el objetivo es, que los participantes del curso aprendan a usar el software HEC-
RAS y los conceptos básicos de la ingeniería sobre el cual está fundamentado el
programa.
Los canales de flujo de aguas de superficie libre o lámina libre, son los que tienen el
espejo de agua en contacto con la atmosfera, la principal fuerza que produce el
movimiento de las aguas es la componente horizontal de la fuerza de gravedad. Los
canales se dividen en dos grandes grupos: artificiales y naturales.
………………………………………….(1)
Rh = Radio hidráulico
……………………………………………………………………(2)
Los canales artificiales son construidos por el hombre para diferentes aprovechamientos
hidráulicos, como la agricultura, generación hidroeléctrica, entre otros. Se les conoce
como canales prismáticos, porque la sección no varía longitudinalmente. Las secciones
transversales más comunes son: trapezoidal, rectangular, triangular, semicircular.
Tenemos canales de tierra, mampostería de piedra y revestidos de concreto de fc= 175
kg/cm2. En el diseño de canales artificiales se toma en cuenta diversos factores
hidráulicos y topográficos. A continuación una vista de un canal revestido.
Figura 2: Vista del espejo de agua del río Amazonas cerca a Iquitos
III. CLASIFICACIÓN DEL FLUJO DE AGUAS SUPERFICIALES SEGÚN LA VARIABLE
TIEMPO
1 2
V2
Ht
2g Sw Línea de energía
Z So
x
Zb
Nivel de Referencia
…………………………………………………………. (3)
Ecuación de continuidad:
……………………………………………(5)
Q Z
T 0 ………………………………………………………………………………….(6)
x t
Q ( Q 2 / A) Z
gA gAS f 0 ……………………………………………………….(7)
x x x
La Hidrología es una ciencia que estudia o versa sobre el agua que existe en sus
diferentes estados sobre la superficie terrestre y capa atmosférica adyacente a la tierra,
trata sobre su distribución, propiedades físicas y químicas, su relación con el medio
ambiente, que incluye su relación con los seres vivos, su circulación sobre y debajo de la
superficie terrestre.
La Hidrología se divide en dos grandes grupos: Hidrología Superficial, Hidrología
Subterránea
permiten saber cuánto asciende o desciende el nivel del rio, cuando correlacionamos
el nivel y caudal podemos estimar un comportamiento futuro, también lo usamos
como datos de entrada en simulaciones de flujo de corrientes naturales. A
continuación se muestra una regla limnimétrica.
Una cuenca hidrográfica es un área de terreno que drena agua en un punto común,
como un riachuelo, arroyo, río o lago cercano. Cada cuenca pequeña drena agua en
una cuenca mayor que luego desemboca en el océano.
La morfología trata de las formas en planta que adopta un río, los ríos de la cuenca
Amazónica correspondiente a la selva baja, son bastante meandricos debido a la baja
pendiente y al gran transporte de sedimentos de sus aguas, es raro encontrar tramos
rectos, la mayoría son trenzados o sinuosos, el cauce sufre modificaciones a causa de la
sedimentación y erosión.
Se muestra algunos tipos de cauce en vista en planta para ríos de baja pendiente y
lechos inestables.
Tenemos muchos lenguajes de programación para ingeniería y ciencia, que nos permite
automatizar los procesos y ahorrar tiempo en los cálculos, permite hacer simulaciones
con salidas graficas y tabulares, constituye una herramienta importante para las
soluciones de problemas ingenieriles y científicos, podemos citar algunos, como Fortran
90, C++, Java, Visual Basic, Matlab, etc.
Para tener una idea del tramo que vamos a modelar, es necesario tener una vista en
planta de la rivera, islas del río, es recomendable que sea en Autocad, ya que nos
permite hacer mediciones entre secciones, ver la sinuosidad del rio, plotear los puntos del
eje, de aquí obtenemos algunas mediciones que serán datos de entrada para el
programa Hec-Ras. Los datos para referenciar en planta se seleccionan de la vista en
planta y del procesamiento en el programa Hypack. Nos permite distinguir la posición
correcta de las márgenes derecha e izquierda de cada sección transversal, cuando solo
se tiene datos tabulares y no tenemos la vista en planta, es difícil saber cómo se
distribuye los datos en la vista en planta. En seguida se muestra la vista en planta del
tramo que vamos a usar para nuestro modelamiento hidráulico.
270859.82 271047.43
9476852.40 9476829.73
AFORO PUERTO AMERICA
La obtención de los puntos de las secciones transversales se realiza a partir de los datos
de control vertical. Son datos obtenidos en el levantamiento hidrográfico, tenemos a los
sondajes o profundidades de cada sección transversal del río. Las lecturas de las reglas
limnimetricas en por lo menos dos estaciones, de preferencia al inicio y final del tramo.
Convertimos los datos de niveles del río en cotas geoidales o referidos al nivel medio del
mar, en seguida establecemos una pendiente hidráulica promedio del tramo para la fecha
de levantamiento hidrográfico, luego usamos el concepto de pérdida de carga por fricción
longitudinal, esto nos permite calcular las cotas referenciadas del nivel de espejo de agua
en cada sección transversal para la fecha del levantamiento hidrográfico, a esta altura
geoidal se resta las profundidades de su respectiva sección, así obtenemos el eje de las
ordenadas, que son las cotas de fondo de la sección transversal.
Para las abscisas calculamos la distancia entre los puntos (X, Y) de referencia en planta
de cada sección transversal y con estas distancias acumuladas se forma el eje horizontal
de la sección. Se ilustra una sección transversal en su respectiva cota o altitud.
Al realizar el paso anterior hemos tabulado estos datos, simplemente son los puntos (X,
Y) donde se hizo los sondajes en cada sección transversal, estos permite referenciar en
planta las secciones, el programa nos solicita como mínimo dos puntos que deben ser los
extremos de cada sección. Se ilustra estos datos con una grafica.
Figura 10: Se muestra los coordenadas (X, Y); para referenciar en planta
El programa nos pide calcular las distancias entre el extremo izquierdo, central y derecho
de las secciones adyacentes de todo el tramo de canal elegido, esta longitud se traza con
el comando línea y luego haciendo clic derecho elegimos propiedades, esto nos permite
visualizar la longitud de esta línea. Para nuestro caso empezamos a medir desde aguas
abajo hacia arriba. A continuación se ilustra con un esquema.
Figura 11: La línea celeste representa las distancias izquierda, central y derecha.
6.5 Elegir la Ubicación de Bancos Izquierdos y derechos de las Secciones del Tramo de Río
Estos puntos establecen el ancho de las riberas y del canal principal, en nuestro caso
tenemos datos del canal principal, por lo tanto elegimos los extremos del ancho de la
sección transversal del flujo
Figura 12: Los puntos rojos indica la ubicación del banco izquierdo y derecho.
Cuando creamos un nuevo proyecto en HEC-RAS, tenemos que seguir una secuencia de
ingreso de datos, en seguida se explica los pasos más importantes del ingreso de datos de
referencia vertical y horizontal del tramo de río.
Debemos tener los puntos (X, Y) en coordenadas UTM del eje del tramo que hemos
elegido, estos puntos deben estar ordenados desde aguas arriba hacia aguas abajo.
Vamos a Gis Tools, aquí cliqueamos en Reach Invert Lines Table y luego hacemos un
copia y pega de estos puntos, que previamente lo tenemos tabulado en Excel. Luego
tenemos que referenciar el eje del canal con estos puntos. La figura 17 muestra el tramo
referenciado.
Hacemos clic en Cross Sección Data, opción Add New Section, escribimos el número de
estación, luego copiamos los datos de Excel a la cuadricula de Hec Ras, donde
previamente lo hemos puesto en cotas y ancho de la sección. Si nuestro tramo de rio
presenta un flujo lento o subcritico, entonces empezamos a ingresar las secciones desde
aguas abajo hacia arriba.
7.6 Ingresamos las Distancias entre Secciones, entre Márgenes Izquierdas, Derechas y
Central
Previo cálculo de las distancias de las márgenes izquierdas, centrales y derechas entre
las secciones contiguas, estos son distancias respecto a una sección que esta aguas
abajo. Teniendo la tabla de distancias en tres columnas y correspondiente a cada
estación, vamos al menú Tables>> Reach Lengths, como se muestra en la figura 20.
Para el cálculo de las distancias se usa las herramientas de dibujo de Autocad, otra
manera de calcular estas distancias es ordenando los puntos extremos izquierdos,
centros y derechos respectivamente, luego aplicamos la formula de distancia de un punto
a otro.
7.7 Ingresamos Datos de los Bancos Izquierdos y Derechos de cada Sección Transversal
Los bancos forman parte de la ribera o planicie de inundación, los puntos izquierdo y derecho
delimitan el canal de los bancos. Los puntos deben ser elegidos teniendo en cuenta la información
topográfica y batimétrica de cada sección disponible. Después de corrido el programa estos
puntos pueden editarse usando la herramienta correspondiente, se debe tener en cuenta que la
ubicación de estos puntos es de izquierda a derecha, y el banco izquierdo es menor que el
derecho.
7.8 Ingresamos los Valores del Coeficiente de Manning para cada Sección Transversal
El valor del coeficiente de manning, mayormente lo obtenemos de tablas que han sido elaborados
en base a experimentos para diversos tipos de cauces. Si tenemos datos de aforos, geometría y
parámetros hidráulicos de secciones de medición de caudal, podemos determinar el valor del
coeficiente, usando la ecuación de manning. Esta ecuación calcula la velocidad y caudal para flujo
permanente. La ecuación de manning se escribe así:
…………………………………………………………….......................(8)
Estos datos son los que nos mostrará las características hidráulicas del flujo, establecer una
situación del nivel del rio, en síntesis nos permite modelar el comportamiento hidrodinámico del
río para situaciones extremas o eventos extraordinarios que suceden en un periodo de tiempo.
Podemos ingresar caudales cada cierta distancia en las secciones ingresadas, pero no en secciones
interpoladas, en cualquier rio el caudal aumenta hacia aguas abajo, esto se debe al aporte de sus
tributarios o afluentes.
El programa nos pide establecer las condiciones de contorno o limite, para nuestro caso, nos pide
la pendiente hidráulica aguas abajo de nuestro tramo, está pendiente lo calculamos haciendo el
cociente del desnivel de espejo de los extremos y la longitud del tramo.
Figura 27: Ingreso de datos para condiciones de contorno o límite aguas abajo.
Servan Vargas, Mario 2009 SEHINAV - IIAP- UNAP Página 27
Modelamiento Hidráulico en Ríos de la Amazonía con HEC- RAS
Después de haber ingresado los datos geométricos de referencia horizontal y vertical, de flujo,
pendiente hidráulica, las condiciones limite, nivel del espejo aguas abajo y demás datos que
nos solicito el programa, procedemos a realizar el análisis del flujo en estado permanente, para
ello creamos un plan de simulación para los datos ingresados, en seguida corremos el
programa, y ya estaremos listo para visualizar los resultados usando los respectivos iconos.
Visualizamos los cortes transversales del flujo, perfiles, vista en planta, variación de la velocidad
en el tramo, espejo de agua versus caudal para cada estación, tabla de resultados de las
variables y parámetros calculados por el programa.
Podemos variar los datos de flujo y tendremos un nuevo comportamiento del flujo, si tenemos
datos para condiciones de vaciante o avenidas, el programa simulara un perfil para cada
situación o estado del rio, esto permite tomar decisiones en determinados problemas de
eventos extraordinarios en ríos.
128 Ground
126 Bank Sta
124
122
120
118
0 20 40 60 80 100 120 140
Station (m)
Figura 31: Vista de la sección transversal del flujo.
128 EG PF 1
WS PF 1
126
0.0 m/s
Elevation (m)
0.6 m/s
120
0.8 m/s
118 1.0 m/s
Ground
116
Bank Sta
114
0 20 40 60 80 100 120 140
Station (m)
Figura 32: Sección transversal mostrando distribución de velocidades
EG PF 1
130 WS PF 1
Crit PF 1
Elevation (m)
Ground
125
120
115
110
0 10000 20000 30000 40000 50000
Main Channel Distance (m)
Figura 33: Perfil longitudinal del espejo de agua y fondo del tramo
Legend
1.2
1.0
0.8
0.6
0.4
0.2
0.0
0 10000 20000 30000 40000 50000
Main Channel Distance (m)
Figura 34: Variación de las velocidades en el tramo
WS PF 1
700
695.384* Ground
690.769*
Bank Sta
686.153*
681.538*
Ground
674.666*
668.*
660
654.666*
649.333*
644.*
639.047*
635.238*
632.381*
628.571*
625.714*
621.904*
616.*
608.*
600
593.75*
587.5*
581.25*
568.888*
558.461*
550.769*
541.538*
535.*
528.75*
522.5*
513.333*
503.333*
494.545*
487.272*
480
466.666*
Para la simulación anterior hemos considerado un flujo subcritico, es decir que el flujo de las
corrientes naturales es lento, y el número de froude es menor a 1, además se estableció las
condiciones límite aguas abajo. Si el flujo es supercritico (Fr>1), entonces el flujo es rápido, esto lo
encontramos en ríos de la selva alta o ceja de selva, los ríos de la sierra tienen mayor pendiente
hidráulica, en consecuencia son bastante rápidos.
En la naturaleza encontramos ríos que tienen tramos de flujo lento y otros de flujo rápido, para
este tipo de flujo nuestro análisis tiene que ser mixto, y requiere las condiciones límite aguas
arriba y aguas abajo. Si nuestro flujo a modelar es supercritico entonces empezamos a numerar
las estaciones desde aguas arriba.
Si Fr=1 el flujo adopta un valor critico que nos sirve como referencia para el análisis de los demás
tipos de flujo, el programa calcula las variables para un estado crítico.
La modelación de la calidad de agua que realiza el Hec- Ras está fundamentado en un esquema
numérico explicito. Resuelve la ecuación de advección y dispersión en una dimensión espacial,
significa que el análisis de la calidad de agua es desde el punto de vista físico y numérico.
Con esta herramienta podemos hacer un análisis de la concentración y dispersión de ciertos
elementos químicos y nutrientes que forman parte de la masa de agua que trasladan los ríos o
corrientes naturales.
Los ríos de la cuenca amazónica, específicamente de la región Loreto, presentan gran transporte
de sedimentos, esto genera cambio continuo del cauce y una nueva configuración en planta del
río. En consecuencia es un factor importante que debemos tomar en cuenta para un análisis más
real del flujo. Tenemos transporte de fondo y transporte en suspensión, todos estos procesos
naturales hace que un cauce de rio sea inestable en el tiempo. Modelar el transporte de
sedimentos es un trabajo que tiene muchas complejidades. Debido al transporte de sedimentos
tenemos zonas de erosión y deposición.
Figura 39: Variación de la sección debido al transporte de sedimentos
facilitará el manejo con criterio del software y a la vez permitirá interpretar los
B) En la obtención de los datos para ingreso al Hec Ras se debe evitar cometer errores de
orden o mala ubicación de puntos, si sucede esto el programa nos notificará del error.
C) Se realiza simulaciones para flujo permanente, es decir las variables cambian a lo largo
D) Las simulaciones son para un flujo subcritico o de corrientes naturales con poca
E) Las simulaciones son para un flujo en época de vaciante, es decir no tenemos las
planicies de inundación.
F) Consideramos un tramo de río que no tiene tributario, esto para hacer más fácil el
ejemplo de simulación.
G) Se debe tener en cuenta que la modelación del flujo permanente con Hec Ras está
gobernado por la ecuación de energía, y por tanto no se considera todas las variables
H) Al margen de simplificaciones del flujo real, el programa Hec Ras nos da una buena
estimación del flujo en un canal natural o artificial, esto nos permite tomar decisiones
XI. RECOMENDACIONES
B) La modelación del río con estructuras hidráulicas queda para un curso avanzado o
permanente, además se debe conocer la teoría que usa la modelación del flujo
impermanente.
XII. ANEXOS
[5] http://www.hec.usace.army.mil/
[6] http://www.construaprende.com/tesis03/2006/09/definicion-y-objetivo-de-la-hi.html
[7] http://www.dhi.es/Software/RecursosH%C3%ADdricos.aspx