Sistema de Modelamiento de Cuencas
Filosofa del Modelo
El Modelo W.M.S. Integra varios esquemas de generacin de hidrogramas mediante el HEC-1 el cual esta diseado para simular la escorrenta de una cuenca mediana. Esta cuenca responde a la precipitacin , que responde a la precipitacin como un sistema de interconeccin de componentes hidrolgicas e hidrulicas. Cada componente del modelo cumple un rol en el proceso de precipitacinescorrenta dentro de una porcin la cuenca, comunmente referenciada como Sub-cuenca. Cada componente representa una unidad de superficie de escorrenta con un cauce un punto de control. La representacin de una componente requiere, un conjunto de parmetros los cuales dan particulares caractersticas a los componentes y relaciones matemticas, que describen el proceso fsico. El resultado del proceso de modelamiento es el hidrogramas de escorrenta de flujo ubicado en un punto especifico del cauce principal de la cuenca.
Introduccin:
Una vez que la cuenca ha sido apropiadamente sudividida en sub cuencas componentes se utilizar el W.M.S para definir e identificarlos parmetros que permitan un completo archivo de entrada para la simulacn mediante el HEC-1. Los parametros geomtricos de cada componente como son reas, longitudes y gradientes son obtenidos automaticamente mediante la sub rutina TIN Geometry. Mientras que los parmetros tales como: Tasa de perdidas, Caudal base, Nmero curva,, y tambin es posible la evaluacin del IU y el Ruteo, introducidos atravez de Caja de dilogo interactivo.
Objetivo:
Generar Hidrogramas Sintticos, de precipitacin y caudales obtenidos en campo en este caso en a la cuenca del ro Torata. Estos resultados sinteticos se analizan y calibran con la informacin de caudales registrados en el punto control establecido dentro del cauce principal de la cuenca denominado Arondaya.
Meta:
Generar Hidrogramas de escorrenta directa en tiempo real, con el objeto de determinar el Caudal pico y el Volumen asociado a este en el punto de control para obtener el incremento de nivel volumen en el embalse y desarrollar la operacin de descarga..
Pasos Bsicos para Construir el Modelo
Definir Objetivos Disear la Malla del Modelo
Recopilar Datos Asignar/Modificar Parmetro, Variables y condiciones de Borde, a cada componente Identificar Bordes del Modelo Calibrar y Validar el Modelo por componente Construir un Modelo Conceptual
NO
Es correcta la integracin SI
NO
Suficientes Datos?
SI
Anlisis de sensibilidad
Resultados Pronosticados
Fase I
Fase II
PRESENTACION DEL SOFTWARE: W.M.S Comandos Usados para delimitacin y edicin de la cuenca.
Menus Modules
Tools
Macros
MODULOS Tin Module.
Dem Module Map Module. Tree Module. 2D Grid Module. 3D Scatter Module. HERRAMIENTAS
Definicin de Parmetros y variables usados en el Modelo
Los parmetros y variables utilizados en el modelamiento de la cuenca son los siguientes:
Parmetros:
Area de la cuenca (A). Pendiente del cauce principal (S). Longitud del cauce principal (L). Longitud del centroide a la salida de la cuenca (Lca). Deficin de escala real en base a tres puntos de coordenadas para definir los lmites de la cuenca imagen TIF. Coeficientes de rugosidad de Manning en tramos a lo largo del cauce de la componente (n). Nmero de Curva, asignado en cada componente (CN). Tiempo de Retardo (Tlag) y tiempo de Concentracin a cada componente. (Tc). Coeficiente de almacenamiento del cauce: k, x. Ancho del cauce principal en cada Sub-cuenca (B). Precipitacin (P), Evaporacin (ETP), Aporte de nieve (N). Descarga (Q)
Variables:
Area de la cuenca.
El rea de la cuenca puede ser generada en forma automtica por el W.M.S. mediante el uso AutoCad. Se han establecido 5 Sub-cuencas componentes cuyos puntos de control los definen secciones de registro ya establecidos.
Area (Km 2) 9.0 16.9 22.1 107.2 125.8 281.0
S4
EST. ARONDAYA
Cuenca S1 S2 S3 S4 S5 TOTAL
EST. TITIJONES
S5
TUNEL DE DERIVACION
S3
S2 S1
EST. ICHUPAMPA RESERVORIO TORATA
EST. EL CAON
Pendientes del cauce principal.
Las pendientes han sido calculadas con los datos de coordenadas a partir de un plano de AutoCad:
Sub-Cuenca
Titijones
S1 S2 S3 S4 S5
Arondaya R. Torata
Pendiente (m/m) 0.036 0.039 0.053 0.057 0.026
Ichupampa
El Can
Longitud del cauce principal.
La longitud del cauce principal es obtenida mediante el uso del Autocad.
Longitud Longitud Sub-cuenca Acumulada (Km) (Km) S5 10.82 10.82 S4 12.07 22.89 S3 5.20 28.09 S2 7.72 35.81 S1 4.45 40.26
S4 L5
S5
S3 S2 S1 L2 L1 L3
L4
Longitud del centroide a la salida de la cuenca.
Los Lca son obtenidos, definiendo primero los centros de gravedad de cada cuenca.
Longitud Sub-Cuenca del cauce (Km) S1 4.45 S2 7.72 S3 5.20 S4 12.07 S5 10.82 TOTAL 40.260 Longitud (Lca) (Km) 2.660 4.210 3.160 6.940 7.360 24.330
S5
C.G
S4 Lca5
C.G
S3
C.G
Lca4
S2
S1
C.G C.G
Lca3 Lca2
Lca1
Coeficiente de rugosidad de Manning.
Los coeficientes de rugosidad de manning son definidos segn las caractersticas del cauce del ro, como son: material de fondo en relacin al tirante, la cobertura vegetal, sinuosidad, etc. todos estas caractersticas fueron observadas durante el trabajo monitoreo de flujo en campo. El coeficiente puede modificarse en funcin al tirante.
Tramo S1 S2 S3 S4 S5 Coeficiente de Manning 0.050 0.055 0.070 0.075 0.045
Seccin Transversal del cauce
Nota:
Los coeficientes medidos son valores promedios puesto que varian a partir del tirante.
Nmero de Curva.
Al representar en grficas la informacin de Pe (Escorrenta efectiva) y P (precipitacin total) para muchas cuencas, el SCS encontr curvas tipos como los que se muestran en la figura. Para estandarizar estas curvas, se define un nmero adimensional de curva CN tal que 0 CN 100. Para superficies impermeables y superficies de agua CN=100, para superficies naturales CN 100.
Curva Numero (CN) 65 70 85 85 85
Sub-Cuenca S1 S2 S3 S4 S5
La sub-cuenca S5 tiene una variacin de CN entre 65 y 85, observado en base a registros y simulacin.
Tiempo de Retardo (Tlag) y Tiempo de Concentracin (Tc).
Con los parmetros de longitud del cauce principal y longitud del centroide a la salida de cada Sud-cuenca se obtiene: Mtodo Snyder:
Tlag 1.35 (L Lca)0.3
Mtodos
Sub-Area S1 S2 S3 S4 S5
L: Km CN: Adimencional S: Km/Km T: Horas
Snyder
SCS
Taylor
Mtodo de SCS:
1000 ( 9)0.7 Tlag L0.8 CN 1900 S
Tc 0.6* Tlag
Riverside County 0.8 1.1 0.8 1.5 1.7
2.8 3.8 3.1 5.1 5.0
Tlag (Hora) 1.4 5.0 2.0 6.5 1.1 4.5 2.1 7.1 2.9 10.4
Precipitacin:
Los datos de precipitacin son obtenidos de todas las estaciones pluviomtricas y pluviogrficas instaladas en la cuenca, con estos datos se determinar la escorrenta superficial, en cada Sub-cuenca componente.
FEBRERO: 2001 Estaciones Meteorologicas FECHA Villa Cuajone (mm) 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 TOTAL 2.3 8.0 6.9 3.5 3.0 1.2 0.7 0.7 1.8 0.0 0.0 0.0 0.4 0.0 2.0 11.0 0.7 0.6 0.0 2.8 2.3 10.4 4.9 8.4 7.8 0.0 32.0 6.7 118.10 Mina Cuajone (mm) 6.3 1.2 10.6 5.8 3.7 22.8 1.9 2.3 6.0 3.9 2.0 1.6 4.7 8.5 14.7 12.3 1.8 6.3 1.4 2.3 7.8 4.0 2.8 6.0 6.1 3.6 22.7 5.6 178.70 Via Blanca (mm) 6.8 1.8 15.8 6.3 8.5 4.6 2.9 4.4 6.3 6.9 2.2 9.7 20.0 6.3 14.0 17.5 3.1 1.9 2.3 3.3 9.8 5.8 3.2 21.8 7.4 3.4 19.2 8.5 223.70 TW4 (mm) 4.1 0.2 22.1 3.8 3.4 0.4 2.3 3.9 18.9 4.3 5.3 7.4 2.2 2.3 19.1 2.6 3.5 2.0 5.7 6.5 0.6 0.2 1.0 17.2 4.6 1.6 20.0 6.9 172.10 Titijones Apacheta (mm) 7.4 0.4 15.4 4.4 1.6 0.5 5.0 1.6 10.6 3.8 9.8 5.1 3.8 6.8 22.4 12.0 0.1 7.2 9.3 2.0 0.7 13.2 2.2 27.7 11.5 0.2 16.1 11.0 211.80 Suches lado cuajone (mm) 3.0 0.8 14.7 1.4 1.6 0.1 2.3 0.0 13.0 2.0 4.5 5.6 3.7 9.6 24.3 10.4 0.4 2.3 5.4 0.6 0.7 4.9 4.8 18.2 10.5 0.5 11.8 5.4 162.50
PRECIPITACION - ALTITUD
No. 1 2 3 4 5 6 ESTACION Villa Cuajone Mina Cuajone Via Blanca Titijones TW-4 Titijones Apacheta Suches - Cuajone PRECIP. ALTITUD MENSUAL msnm mm 2,795 118 3,580 179 3,945 224 4,553 172 4,652 212 4,452 163
RELACION PRECIPITACION - ALTITUD
ESTACIONES METEOROLOGICAS
1,000
PRECIPITACION (mm)
1 100
3 6 4
y = 5E-08x3 - 0.0006x2 + 2.5628x - 3203.5 R = 0.81
10 2,500
2,700
2,900
3,100
3,300
3,500
3,700
3,900
4,100
4,300
4,500
4,700
4,900
ALTITUD (msnm)
Aplicacin del Modelo
FASE I:
1.- Recopilacin de informacin: Estacin de control, Precipitacin, Caudales, Cartografa de la cuenca. 2.- Identificacin de la cuenca y sub-cuencas. 3.- Definir el modelo conceptual: Modelo serie.
Aplicacin del Modelo
FASE II:
1.- Preparar TIN, MAP para delimitar las componentes de la cuenca y la base de datos para el desarrollo de los hidrogramas. 2.- Calcular los atributos geomtricos de cada Sub-Cuenca. 3.- Definir las caracteristicas de la variable de precipitacin . a las estaciones y registros para la calibracin. 4.- Calcular la precipitacin efectiva en cada sub cuenca utilizando el procedimiento del SCS en base al Nmero de Curva (CN) otro alternativo: . 5.- Definir el mtodo del clculo del Hidrograma Unitario ms adecuado a cada Sub-cuencas. 7.- Segun el archivo de entrada para la integracin de los flujos de cada Subcuenca mediante el HEC-1. 8.- Generar el hidrograma resultante y mostrarlo grficamente. 10.- Anlizar en forma comparativa el hidrograma simulado y el observado: Proceso calibracin. Hacer la modificacin necesaria por Sub-cuenca.
Escalamiento de la Cuenca Registro de Imagen:
Seleccionando el icono Map Module y luego selecccionando el Import del comando del menu Image, se registra la imagen a escala real que podra ser usada en la definicin de los lmites de la cuenca y puntos de control.
Creando objetos mediante el TIFF Images:
Registrado la imagen en el Terrain Model, seleccionamos el Map Module , usando la herramientas en este modulo, limitamos la cuenca, definimos nuestros puntos de control y el curso principal de la cuenca, mediante Streams y Generic, que vienen hacer los tipos de curvas que generan el lmite y el curso del flujo de la cuenca.
Oulet Point
Basin Boundary
Basin Polygons Stream
Identificacin del Job de Control:
Job Control. El dilogo Job Control se usa para definir informacin general acerca del HEC-1. Se ingresa a este cuadro de dilogo seleccionando el comando Job control del menu del HEC-1, se usa para definir en el box de dialogo el hidrograma, utiliza los factores de tormenta/Area y cambiar los hidrogramas de ingreso
Seleccin de Unidades:
este cuadro permite seleccionar las unidades con las que se desea trabajar, area, distancias, y si se desea clcular otros parmetros de la cuenca, solo se selecciona la opcin Display Opts y seleccione el parmetro deseado.
Introduccin de los datos:
Los atributos o parmetros para las Sub-cuencas son definidos y/o editados en este cuadro de dilogo, luego de haber definidos los lmites y puntos de control de la cuenca. Para ingresar a este cuadro de dilogo, seleccionamos en el menu principal el Tree y luego el comando HEC-1.
Tipos de datos
Otra forma, es ejerciendo doble click sobre la figura de la Subcuenca definida en arbol logico del Hydrologic Modeling en pantalla.
Los parmetros de ruteo apareceran al hacer click sobre el punto de control.
INGRESO DE DATOS:
1.-
Datos de la Sub-cuenca: (Basin Data) la informacin para cada Sub-cuenca es introducida por cada una en base al dialogo de datos de lectura. Los datos de rea, caudal base etc. pueden ser ingresados o generados automticamente al tener un Triangulated Irregular Network ( TIN).
STRTQ
Edicin del Hidrograma registrado:
El editor X-Y se usa para introducir la curva definida Q-T. La curva de hidrograma puede ser importada o exportada de xy hacia un archivo de texto (txt).
2.- Dilogo de Salida:(Ouput Control)
Este control permite mostrar la informacin sobre las estaciones hidrogrficas, y esta escrita en archivo ASCII procedente del HEC-1. Por defecto la informacin es registrada en el IO record. Por defecto la opcin para mostrar el hidrograma de salida usa el campo TAPE22, del dialogo mostrado.
Ubicacin de estaciones:
El Gage dialog permite definir de una lista de estaciones segn la necesidades del trabajo. Las estaciones pluviomtricas con o sin registrador son definidas en este dilogo: ubicacin, precipitacin total y la distribucin de la precipitacin.
S5
LCUAJ P. TW 4
S4
HIDRO 1
P. TA
S3
ICHU CHUNT
S2
S1
MINA CUAJONE VC CM VB
Distribucin de la precipitacin:
Poligonos de thiessen: Los poligonos de thiessen se generan automaticamente al ingresar todas las estaciones con sus respectivos coordenadas. No existe opcin de correccin.
3.- Precipitacin: (Precipitation)
El patrn de precipitacion es asignado a cada estacin seleccionando el comando Precipitation del Edit HEC-1 parameter.
las estaciones Pluviogrficas y no Pluviometricas que intervienen en cada Sub-cuenca provenientes del dialogo HEC-1, las estaciones aparecen y se identifican con sus pesos ponderados.
Nombre de c/u de las esatciones y el tipo.
La identificacin de las estaciones aparesen en el dialogo del HEC-1, Gage.
4.- Seleccin del mtodo para el
clculo de prdidas: (Loos Method) Se selecciona el mtodos ms adecuado de prdidas para generar los Hidrograma sinttico. Un mtodo de prdida es asignado de acuerdo a cada Sub-cuenca y por seleccin del botn de Loss Method.
En este caso se opt por el mtodo del Soil Conservation Service (SCS), el cual utiliza el numero de curva (CN).
5.- Hidrograma Unitario:
(UH Method) Seleccin del mtodo de clculo Se puede definir el mtodo de clculo del (HU) a partir de este box de dilogo. El mtodo puede ser asignado a una Sub-cuenca escogiendo el Unit Hydrograph Methods.
Si se dispone de archivos Gis y segun el archivo TIN, es posible el clculo directo con el boton Calculate Unit Hydrograph Hydrograph Parameters
Generacin de parmetros: Tlag y Tc:
Si contamos con los datos TIN, este dilogo calcular los parmetros directamente al utilizar el box anterior en los mtodos indicados.
6.- Aporte de Nieve: (Snow Melt Data)
Cuando en el anlisis de escorrenta es necesario considerar, los datos de derretimiento de nieve para una Subcuenca se utiliza el este dilogo. El HEC-1 tiene dos diferentes mtodos para calcular snowfall/melt: Degree-Day Energy-Budget.
Se puede alegir cualquiera de los mtodos e introducirse los parmetros para el clculo.
7.- Datos para el routing:
(Routing Data) Routeo de integrado. hidrograma
Este cuadro permite definir los puntos donde los hidrgramas van a ser combinados, adems de las caractersticas geomtricas del cauce y de prdida a lo largo del trayecto del flujo, entre puntos de control definir las caractersticas de algunas estructuras de regulacin y derivacin que se encuentran en el trayecto del flujo. Los datos de ruteo son introducidos a travez de la seleccin de un punto de control y luego con el Routing Data.
Incorporacin de las derivaciones: Las derivaciones son creadas usando una combinacin del comando Add y Retrieve Diversion definidas en el Tree menu. Una vez creada la seccin de derivacin, los datos del detalle pueden ser definidos y/o editados por seleccin del botn Diversion Data.
Chequeo control de datos para procesar el Modelo:
El comando Model Check puede ser accionado una vez que se tengan todos los datos necesarios para el hidrgrama y el HEC-1. Este comando puede reportar cualquier posible error/inconsistencia en el modelo, la correccin de cualquier error/inconsistencia se realiza haciendo un click donde aparezca un error en la ventana que se muestra:
Corrida del Modelo:
Esta versin puede correrse directamente del W.M.S. usando el comando Run HEC-1 en el HEC-1 menu.
SALIDA GRAFICA DE ARCHIVOS
Exportar a Txt o Dat
Hidrograma Observado Hidrograma Simulados
Tlag = 1.5 horas N=94 f=0 Tlag = 2.8 horas N=85 f = 2.57
Tlag = 4.9 horas N=83 f = 1.49
Hidrograma Simulados Tlag = 4.9 horas N=85 f = 1.32 Hidrograma Observado
Tlag = 4.8 horas N=89 f=0
Concluciones:
La distribucin de la precipitacin al paracer no es suficiente para generar el hidrograma que se aproxime al hidrograma registrado en las estaciones
hidrometricas.
Los Numero de Curva, resultado del modelamiento son muy altos, lo cual indica una superficie impermeable, lo cual no es real. Las estaciones pluviometricas, no estan ubicadas en forma homogenea en la
cuenca, lo cual provoca que la distribucin sea ineficiente.
Bibliografa
Tutorial del W.M.S.; Engineering Computer Graphics Laboratory; Brigham young University, 1997. Manual del W.M.S.; Engineering Computer Graphics Laboratory; Brigham young University, 1997. Hidrologa Aplicada; Vente chow, Daniel R. Maindmont, Larry W. Mays, 2000. Manual Silvo Agropecuario; Servicio silvo agropecuario (SESA); Universidad Nacional Cajamarca (UNC) . Manual de HEC-1; U.S-Army Corp of Engineering; Hidrologic Engineering Center, 1998. Informe de afluentes y afloramientos del ro Torata; Recursos HdricosCuajone; SPCC, 2001