ESTUDIO DE HIDROLOGIA E INUNDACIÓN
PROYECTO PLANTA AGUAS SERVIDAS LLAY-LLAY
Consultor Elaborado para
SEASA, Servicios Ambientales y Sanitarios Municipalidad de Llay-Llay
Gabriel Marco Zalaquett
Informe Nº: Rev_02 Proyecto: “Planta de Tratamiento de Aguas Servidas
Porvenir”
Elaboró: Claudio Gaete, Ingeniero en Medio Ambiente Estudio Riesgo de Inundaciones, Febrero 2022
Contenidos
Estudio de Inundación Planta Aguas Servidas Llay-Llay................................................................................ 1
Introducción .................................................................................................................................................. 3
Objetivos ....................................................................................................................................................... 3
Antecedentes ................................................................................................................................................ 4
Caracterización Hidrológica .......................................................................................................................... 4
Cuenca Aportante Y Precipitaciones Medias ................................................................................................ 6
Temperaturas................................................................................................................................................ 9
Geomorfología y Geología .......................................................................................................................... 10
Análisis Hidráulico Canal “Las Masas” ........................................................................................................ 13
Coeficiente de rugosidad de Manning .................................................................................................... 14
Condiciones de Borde ............................................................................................................................. 15
Implementación en Hec-ras 5.05 ............................................................................................................ 16
Resultados del Modelo ........................................................................................................................... 17
Conclusiones ............................................................................................................................................... 20
Recomendaciones ....................................................................................................................................... 20
Carpeta fotográfica ..................................................................................................................................... 20
Introducción
En la comuna de Llay-Llay, región de Valparaíso se llevará a construcción una Planta de Tratamiento de
Aguas Servidas. La tramitación de esta requiere diversos estudios, entre ellos el de riesgo de inundación
por algún cauce cercano. El proyecto se ubica a unos 500 m de la ruta 5 caletera (Figura 1), predominando
el uso agrícola del suelo y la instalación de empresas de dicho rubro.
Un estudio de inundación debe caracterizar todos los cauces cercanos al proyecto, ya sean naturales o
artificiales. El análisis de caudales puede ser por crecidas asociadas a una precipitación máxima, en el caso
de ríos y cuerpos fluviales naturales. Por el lado de los canales artificiales, los principales datos
corresponderán a los aforos realizados en terreno.
Figura 1. Ubicación del Proyecto Planta de Aguas Servidas.
Objetivos
El objetivo general del presente estudio es poder determinar el riesgo de inundación en el sector donde
se emplazará la plata de tratamiento de aguas servidas de sector El Porvenir Llay-Llay, a partir de esto se
definieron los siguientes objetivos específicos:
• Revisar y analizar la información disponible de la zona del proyecto
• Identificar, delimitar la cueca hidrológica en el área del proyecto y la precipitación asociada para
distintos periodos de retorno
• Identificar los cauces naturales y/o canales artificiales cercanos al proyecto que constituyan un
riesgo de inundación para este
• Estudiar el riesgo de inundación por saturación del suelo
Antecedentes
Se han recopilado los antecedentes disponibles de distintas instituciones, principalmente públicas, como
la Dirección General de Aguas, Dirección de Obras Hidráulicas (DOH) y Declaraciones de Impacto
ambiental (DIA) aprobadas por el Servicio de Evaluación Ambiental (SEA).
• DGA-BNA (2020). Estadísticas de precipitaciones, temperaturas y caudales en las estaciones
DGA próximas a la zona de estudio.
• DGA-IPLA (1987). Balance Hídrico de Chile. Ministerio de Obras Públicas, Santiago, Chile.
• DGA, U. Chile & PUC (2017). Actualización del Balance Hídrico Nacional, SIT N° 417. Ministerio
de Obras Públicas, Dirección General de Aguas, División de Estudios y Planificación. Realizado
por la Universidad de Chile & Pontificia Universidad Católica de Chile.
• DGA-MOP (1986). Estudio del Mapa Hidrogeológico Nacional. Escalas 1:1.000.000 y
1:2.500.000. Santiago.
• DGA (2016). Disponibilidad de recursos hídricos para el otorgamiento de derechos de
aprovechamiento de aguas subterráneas en el Valle del Aconcagua Sectores hidrogeológicos
de San Felipe, Putaendo, Panquehue, Catemu y Llay Llay. SDT N°387. Informe técnico DARH
N° 333. Expediente VAR-0503-1. Santiago.
• Ampliación y Mejoramiento Agrícola e Inmobiliaria Vichiculén (INRHED, 2021)
Caracterización Hidrológica
El área de estudio se encuentra en inserto en la subsubcuenca DGA “Estero Los Loros” (338 km2),
perteneciente a la subcuenca DGA “Aconcagua Bajo” (2.626 km2), y esta a su vez perteneciente a la gran
cuenca DGA “Río Aconcagua” (7.333,8 km2), que actualmente presenta un régimen exorreico activo
(Figura 2).
Figura 2. Cuenca de Ubicación.
No existe presencia de cauces naturales cercanos, al este del proyecto, destaca el canal “Las Masas”
(Figura 3), de carácter intermitente, donde su caudal aportante corresponde aguas sobrantes de riego de
los predios vecinos y posteriormente se conecta a una serie de canales de la misma índole. El caudal del
canal “Las Masas” se determinó mediante aforos, generando un caudal medio de 160 l/s.
Figura 3. Red de drenaje de la zona, natural y artificial.
Cuenca Aportante Y Precipitaciones Medias
Al calcular la cuenca aportante (Figura 4) a la zona de proyecto, cuyas características principales se
resumen en la tabla, es fácil notar que la red de drenaje DGA no alcanza la zona de proyecto, siendo el
Estero “Los Loros” el más cercano, a una distancia de unos 4 km, por lo que no representa riesgo de
inundación. En la zona de proyecto, domina fuertemente la canalización antrópica, cuya administración
recae, en la Sociedad de Canalistas de Llay-Llay.
Figura [Link] Aportante al Proyecto.
La siguiente tabla muestra que el área de la cuenca es de unos 42 km2, dicho cauce relacionado a la
cuenca no conecta con el área del proyecto, ni su entorno.
Tabla 1. Parámetros de la Cuenca.
PARÁMETROS MORFOMÉTRICOS CAP_01
GENERALES
Parámetro Símbolo Valor
Perímetro (km) P 32.06
Área (km2) A 42.03
Longitud cauce principal (km) LCP 13.41
Desnivel altitudinal (m) DA 1,260.00
PARÁMETROS MORFOMÉTRICOS
ASOCIADOS A LA FORMA DEL
RELIEVE
Cota máxima (msnm) cmax 1,649.0
Cota mínima (msnm) cmin 389.0
Altitud media (msnm) Am 751.7
Pendiente promedio de la cuenca Smed 17.8
(%)
Altitud máx del cauce Amaxc 1,466.0
Altitud min del cauce Aminc 389.0
Altitud media del cauce Amedc 583.0
Pendiente promedio del cauce Sc 8.0
(%)
PARÁMETROS MORFOMÉTRICOS
ASOCIADOS A LA RED DE
DRENAJES
Tiempo de concentración (horas) Tc 1.22
[California Culvert Practice]
Tiempo de concentración (horas) Tc 1.29
[Kirpich]
Sobre las precipitaciones, la estación meteorológica más cercana es la Estación “Las Chilcas” y la cual
cuenta con datos desde 1990 hasta la actualidad. Las mayores precipitaciones se dan en los meses
invernales, junio y julio. Dicho registro es característico de la Región de Valparaíso y la zona central de
Chile. En el gráfico (Gráfico 1) de precipitación media anual vs año (Inrhed, 2021) se observa que en los
últimos 10 años ha existido una merma en la precipitación, dicho efecto se correlaciona con la sequía que
experimenta la región de Valparaíso.
Tabla 2. Precipitaciones Medias en Estaciones Cercanas al Proyecto
Gráfico 1. Precipitación media anual vs año. Fuente: Inrhed, 2021.
Temperaturas
El registro de información disponible de temperaturas en las zonas es muy escaso, dentro de la cunca “Río
Aconcagua” sólo se cuenta con cuatro estaciones que poseen registros de temperatura, como son Llui-
Llui embalse, Los Armos, Quillota y Vilcuya. La temperatura media anual en la zona de estudio es de 14 °C
de acuerdo con las isotermas del Balance Hídrico de Chile (DGA, 1987). Dichos valores son coherentes con
el estimado en la estación meteorológica DGA Quillota que tiene un promedio anual de 14.9°C, la que se
toma como referencia debido a su cercanía con la zona (Gráfico 2).
Gráfico 2. Temperatura Media Anual.
Geomorfología y Geología
La Región de Valparaíso está compuesta por 7 unidades geomorfológicas. Las unidades de mayor
extensión de oeste a este son la Planicie Marina y/o Fluviomarina, los Cordones Transversales, la Cordillera
de la Costa y la Cordillera Andina de Retención Crionival. Las unidades de menor área, también de oeste
a este, son el Farellón Costero, los Llanos de Sedimentación Fluvial y/o Aluvional y las Cuencas
Transicionales Semiáridas. (Inrhed, 2021)
A escala local, la Pta Tratamiento de A.S se emplazan en la unidad de la Cordillera de la Costa, que está
representada por los cerros circundantes al Proyecto. Esta se caracteriza por la presencia de importantes
alturas geográficas que varían entre los 1.800 m.s.n.m. y los 2.300 m.s.n.m. (El proyecto se encuentra en
torno a los 390 m.s.n.m en una planicie en dicha unidad geomorfológica) En el límite oeste de dicha unidad
comienza el desarrollo de las Cuencas Transicionales Semiáridas.
Figura 5. Mapa Geomorfológico de la zona.
Geológicamente el Proyecto se ubica en el interior del Valle de Llay-Llay, en una reducida depresión que
se abre en la ribera sur del Río Aconcagua totalmente rodeada de cerros, compuestos por rocas volcánicas,
volcanosedimentarias e intrusivas del Cretácico. En los cerros ubicados al suroeste del valle, afloran las
rocas más antiguas de la zona de estudio, las cuales se agrupan bajo el nombre de Formación Veta Negra
(Kiag. Estas rocas están intruidas por granitoides de la Unidad Chalinga (Ki1), los cuales también afloran
en los cerros ubicados al norte del Proyecto y que rodean a la ciudad de Llay-Llay. Por otro lado, hacia el
sureste del Proyecto, los cerros están conformados las rocas de la Formación Las Chilcas (Klc). Finalmente,
el relleno del valle lo conforma la unidad conocida como Sedimentos Aluviales, Coluviales y Depósitos de
Relleno de Valles (Q1), representada por la acumulación de material proveniente de las laderas de los
altos topográficos circundantes. La (Figura 6) muestra la distribución de las unidades mencionadas
(Inrhed,2021)
Figura 6. Mapa Geológico local del Proyecto.
Análisis Hidráulico Canal “Las Masas”
Ecuación de Energía Hec-ras
Utilizando el software libre HEC-RAS 5.0.5, desarrollado por el Cuerpo de Ingenieros del Ejército de
Estados Unidos (US Army Corps of Engineers), que permite la modelación en 2D para representar el
escurrimiento. El procedimiento computacional básico se basa en la solución de la ecuación de energía
(Ecuación 1).
2V2 2 1V12
Y2 + Z 2 + = Y1 + Z1 + + he
2g 2g Ecuación 1
Donde:
Y1, Y2 = Altura de agua en las secciones transversales 1 y 2
Z1, Z2 = Cota de terreno las secciones transversales 1 y 2
V1, V2 = velocidades promedio de la sección
α1, α2 = coeficientes de ajuste de velocidades
g = aceleración de la gravedad
he = pérdida de carga
Figura 7. Esquema de la ecuación de energía
Coeficiente de rugosidad de Manning
Para el análisis hidráulico, es necesario caracterizar la rugosidad de las superficies mediante el coeficiente
de Manning como antecedente para el modelo hidráulico, parámetro que representa la resistencia que
ejerce la superficie al escurrimiento del agua, y que se relaciona con el caudal, la pendiente y geometría
del cauce.
𝑄 1 2 1
𝑣= = × 𝑅 ⁄3 × 𝑖 ⁄2 Ecuación 2
𝐴 𝑛
Dónde:
𝑣: velocidad media de la sección transversal de escurrimiento en el canal (m/s)
𝑄: caudal en el canal (m3/s)
𝐴: área sección transversal de la corriente (m2)
𝑛: coeficiente de rugosidad de Manning
𝑅: radio hidráulico (área/perímetro mojado) (m)
𝑖: pendiente en uno por uno (m/m)
El coeficiente de rugosidad de Manning del cauce natural se determinó utilizando la Tabla 3.705.1.A del
Manual de Carreteras (MOP, 2018), el libro “Hidráulica de Canales Abiertos”, Ven Te Chow (1994); la
fórmula de Cowan y también considerando las características morfológicas del canal las Masas.
Cowan (1956) desarrolló un procedimiento para estimar el valor del coeficiente de Manning, el cual queda
definido por la siguiente ecuación:
𝑛 = 𝑚 ∙ (𝑛0 + 𝑛1 + 𝑛2 + 𝑛3 + 𝑛4 ) Ecuación 3
Dónde:
𝑛: Coeficiente de rugosidad de Manning
𝑛0 : Rugosidad base para un canal recto, uniforme, prismático y con rugosidad homogénea
𝑛1 : Rugosidad adicional debida a irregularidades superficiales del perímetro mojado a lo largo del tramo
en estudio
𝑛2 : Rugosidad adicional equivalente debida a variación de forma y de dimensiones de las secciones a lo
largo del tramo en estudio
𝑛3 : Rugosidad adicional equivalente debida a obstrucciones existentes en el cauce
𝑛4 : Rugosidad adicional equivalente debida a la presencia de vegetación
𝑚: Factor de corrección para incorporar efecto de sinuosidad del cauce o presencia de meandros
En la Tabla se presentan los valores de los parámetros que intervienen en la fórmula de Cowan.
Tabla 3. Valores para el Cálculo del Coeficiente de Manning Mediante Método de Cowan.
Condiciones del Canal Valor
Material del Lecho Tierra 𝑛0 0,020
Roca Cortada 0,025
Grava Fina 0,024
Grava Gruesa 0,028
Grado de Irregularidad Perímetro Despreciable 𝑛1 0,000
Mojado Leve 0,005
Moderado 0,010
Alto 0,020
Variaciones de las Secciones Graduales 𝑛2 0,000
Alternándose Ocasionalmente 0,005
Alternándose Frecuentemente 0,010 - 0,015
Efecto Relativo de las Despreciable 𝑛3 0,000
Obstrucciones Leve 0,010 - 0,015
Apreciable 0,020 - 0,030
Alto 0,040 - 0,060
Densidad de Vegetación Baja 𝑛4 0,005 - 0,010
Media 0,010 - 0,025
Alta 0,025 - 0,050
Muy alta 0,050 - 0,100
Sinuosidad y Frecuencia de Leve 𝑚 1,000
Meandros Apreciable 1,150
Alto 1,300
Fuente: Manual de Carreteras (MOP, 2018).
Tabla 4. Coeficientes de Manning según fórmula de Cowan
Variable Cauce principal Riberas
𝑛0 0,020 0,030
𝑛1 0,005 0,005
𝑛2 0,005 0,005
𝑛3 0,01 0,01
𝑛4 0,005 0,005
𝑚 1,0 1,0
𝒏 0,045 0,055
Realizando el cálculo, se obtiene un coeficiente de rugosidad de Manning para el cauce principal de 0,045
y para las riberas 0,055.
Condiciones de Borde
Dado que se desconocen las condiciones de escurrimiento producto de los aforos realizados, se supone
que el flujo puede presentar régimen supercrítico, crítico o subcrítico, por lo que se realiza el análisis para
flujo mixto, es decir, se impone una condición de borde aguas arriba y aguas abajo. Por lo tanto, como
condición de borde se establece que se produzca altura normal de escurrimiento aguas abajo del tramo
de análisis. Para calcular esta altura normal el software solicita ingresar la pendiente del cauce; siendo la
pendiente adoptada aguas abajo (2,0%) la que se estiman de acuerdo con la topografía realizada. Para
aguas arriba se introduce como condición de borde un Hidrograma, el cual se basa en el caudal aforado
de ~160 l/s y un factor de seguridad de un 20%, generando un caudal máximo de unos 200 l/s. De todas
maneras, el canal está diseñado para caudales mayores, además es importante señalar que dicho canal
depende de la Sociedad de Canalistas de Llay-Llay por lo que su flujo es controlado.
Implementación en Hec-ras 5.05
Se modela el Canal Las Masas, el cual posee dimensiones de 1,5 metros de ancho hasta unos 2,5,
profundidades entre 1 y 1,5 metros y un talud de 1:2. La grilla generada para el modelo posee áreas que
van desde los 1,96 m2 hasta los 2.28 m2. Lo anterior con el objetivo de construir un modelo que se
encuentre en el óptimo de calidad de resolución, representación del problema y eficiencia computacional
de proceso. La discretización temporal del modelo es de una hora (variación de caudales entre 0 y 20 l/s),
una iteración cada 10 min y una presentación de los resultados cada una hora.
Figura 8. Parámetros de la Grilla en Hec-ras.
Figura 9. Topografía del sector y Canal Las Masas.
Resultados del Modelo
Una vez ejecutado el modelo, los resultados se presentan en función del tiempo, para este caso, es de
interés conocer los valores máximos, tanto de la altura de escurrimiento, el área de escurrimiento y la
velocidad de este. Se puede observar en la Figura 10 que el canal es capaz de soportar dicho caudal sin
entrar en inundación. Las alturas máximas de escurrimiento están en torno a los 2 m, mientras que las
velocidades máximas se registran valores cercanos a los 3 m/s (Figura 11).
Figura 10. Altura de escurrimiento máxima para un caudal de 200 l/s.
Figura [Link] de escurrimiento para el caudal de 200 l/s.
Conclusiones
• Se concluye que el área del proyecto Planta de Aguas Servidas, comuna de Llay-Llay , se encuentra
en una zona principalmente agrícola, con empresas dedicadas a dicho rubro.
• No existen cauces naturales cercanos al proyecto que puedan poner en riesgo las instalaciones
por inundación, siendo el más cercano el estero “Los Loros”, a unos 4 km de distancia.
• Los flujos superficiales corresponden a canales colectores de caudales sobrantes de la agricultura,
al este del proyecto, a unos 50 m, se encuentra el canal “Las Masas”, dependiente de la Sociedad
de Canalistas de Llay-Llay. El caudal que circula es intermitente, incluso caudal cero, en los aforos
realizados se determinó que el caudal es en torno a los 160 l/s.
• De los resultados obtenidos del programa Hec-ras 5.05, se observa que el canal es capas de llevar
un caudal de al menos unos 200 l/s. El cauce tiene un ancho de unos 2,5 m y una profundidad en
torno a los 1,5 m. Por lo que permite el flujo de caudales muy superiores a los modelados y/o
aforados. Este además se encuentra manejado con la asociación, por lo que, el flujo o no flujo de
agua en este, depende de dicha institución (principal función del canal es recolectar las aguas
sobrantes de riego).
Recomendaciones
Se recomienda realizar un estudio completo de la zona que incluya los contenidos mínimos para una DIA:
• Balance Hídrico natural y con proyecto.
• Estudio Hidrogeológico
• Riesgo de afloramiento de la napa y posterior inundación, tanto naturalmente como por el
periodo de construcción y fundaciones. Realizar calicatas y ensayos en terreno.
• Análisis de calidad de aguas, toma de muestras en el estero.
• Aplicabilidad de PAS156 y/o PAS 157, generando las carteras técnicas correspondientes en cada
caso, si fuese necesario.
• Modelación hidrogeológica (solo si DGA pide explícitamente), de las implicancias del proyecto en
el sector y el área de influencia.
• Afectación a terceros
Carpeta fotográfica
Figura 12: Canal ubicado en las inmediaciones del terreno
Figura 13: Límites Área de Emplazamiento
Figura 14: Canal Zona de emplazamiento
Figura 15: Visita a terreno, Febrero 2022
Figura 16: Terreno donde se instalará PTAS
