Diseño de Dique y Bocatoma para el municipio de La Cumbre

Presentado a:
Ing. Hector Matamoros Rodriguez

Presentado por:
Breydi Santiago Becerra
Silvia Fernanda Cortes
Cristian David Molina

Escuela Colombiana de ingeniería Julio Garavito

Programa de ingeniería civil
Acueductos y alcantarillados

Bogotá D.C
28 de octubre de 2021
Tabla de contenido
1. Introducción........................................................................................................................3
2. Determinación de caudales.................................................................................................3
2.1 Análisis hidrológico..........................................................................................................3
2.2 Curva de duración de caudales..........................................................................................4
2.3 Pruebas de bondad y ajuste...............................................................................................6
2.3.1 Prueba de bondad y ajuste para caudales máximos.......................................................6
[Link] Método Normal...........................................................................................................6
[Link] Método Log-normal....................................................................................................8
[Link] Método de Gumbel......................................................................................................9
[Link] Método Log-Gumbel.................................................................................................10
2.3.2 Prueba de bondad y ajuste para caudales mínimos......................................................12
[Link] Método Normal.........................................................................................................12
[Link] Método Log-normal..................................................................................................13
[Link] Método de Gumbel....................................................................................................15
[Link] Método de Log-Gumbel............................................................................................16
2.3.3 Prueba de bondad y ajuste para caudales medios........................................................17
[Link] Método Normal.........................................................................................................17
[Link] Método Log-Normal.................................................................................................19
[Link] Método de Gumbel....................................................................................................20
[Link] Método de Log-Gumbel............................................................................................21
2.4 Distribución de probabilidad...........................................................................................23
2.4.1 Caudales máximos.......................................................................................................23
2.4.2 Caudales mínimos........................................................................................................24
2.4.3 Caudales medios..........................................................................................................25
3. Traslado de los caudales al punto de bocatoma................................................................27
4. Modelación en Hec-Ras del tramo del río........................................................................27
5. Diseño del Dique...............................................................................................................27
5.1 Diseño del dique WES....................................................................................................27
5.2 Diseño de Dique Ogee....................................................................................................31
 1. Introducción

El trabajo que se presenta a continuación contiene los diseños de las obras de captación, las
cuales son estructuras hidráulicas que tienen como objetivo extraer una cantidad de agua
determinada de una fuente para satisfacer una necesidad y encausarla a un sistema de
transporte. Se presentarán diferentes modelos y tipos de diques y bocatomas que se
requieren para los caudales de diseño calculados a partir de los datos de consumo obtenidos
del municipio de La Cumbre, ubicado en el departamento del Valle del Cauca.
Este trabajo tiene dos propósitos principales, el primero es establecer los caudales de diseño
de consumo y del río del cual se abastecerá el sistema de acueducto y el segundo es,
caracterizar el comportamiento de flujo a partir de los índices de regulación y caudales
ecológicos, con el fin de aseguración la menor afectación posible al río y a la fauna.
Para establecer los comportamientos reales de las láminas y del flujo se implementó el
software Hec-Ras, en el cual se modelo teniendo en cuenta los criterios necesarios para
localizar el sitio de captación del dique y de la bocatoma. Es importante resaltar que la
modelación está basada en las suposiciones iniciales de las dimensiones de las estructuras.

2. Determinación de caudales

2.1 Análisis hidrológico

Se tienen los registros hidrológicos de los caudales diarios por 25 años del municipio de La
Cumbre que empiezan en 1994 y terminan en 2019. Como se observa en la tabla 1 donde se
encuentran los valores de los caudales máximos, mínimos y medios.
 Caudales diarios anuales
Año QMed (m³/s) Qmin (m³/s) QMax (m³/s)
1994 33,397 0,930 149,396
1995 26,803 0,715 114,920
1996 26,450 0,706 112,268
1997 27,068 0,706 114,036
1998 27,288 0,776 114,920
1999 26,784 0,715 105,196
2000 26,868 0,736 110,500
2001 26,542 0,715 112,268
2002 26,515 0,706 113,152
2003 26,594 0,706 106,080
2004 26,658 0,847 114,920
2005 26,591 0,736 110,500
2006 26,868 0,776 114,920
2007 26,936 0,756 112,268
2008 26,517 0,746 109,616
2009 26,909 0,766 114,920
2010 26,512 0,776 106,080
2011 26,366 0,764 106,080
2012 26,885 0,756 110,500
2013 27,233 0,827 112,268
2014 26,029 0,726 107,848
2015 26,386 0,715 111,384
2016 33,851 0,858 149,396
2017 33,693 1,022 134,722
2018 33,010 0,944 145,948
2019 33,387 0,930 135,782
Tabla 1. Resumen de caudales diarios mensuales

2.2 Curva de duración de caudales

Se realiza la curva de duración de caudales utilizando los datos de caudales del sistema de
información nacional ambiental para poder determinar el índice de retención y regulación
hídrica, como también el caudal ecológico.
No. datos 9503
Qmax 149,396
Qmin 0,706
Rango 148,69
14,12694
No. intervalos
15
Tamaño intervalo 9,96
Tabla 2. Parámetros de los intervalos
Para poder calcular la curva de duración de caudales se utilizan los valores de la Tabla 2.
Rango Superior Rango Inferior Rango medio Frecuencia Frecuencia acum. Frec. Relativa acum (%) Q (m³/s)
1 149,396 139,4 144,416133 6 6 0,06 144,4161333
2 139,4 129,5 134,4564 19 25 0,26 134,4564
3 129,5 119,5 124,496667 31 56 0,59 124,4966667
4 119,5 109,6 114,536933 55 111 1,17 114,5369333
5 109,6 99,6 104,5772 159 270 2,84 104,5772
6 99,6 89,6 94,6174667 202 472 4,97 94,61746667
7 89,6 79,7 84,6577333 340 812 8,54 84,65773333
8 79,7 69,7 74,698 512 1324 13,93 74,698
9 69,7 59,8 64,7382667 512 1836 19,32 64,73826667
10 59,8 49,8 54,7785333 368 2204 23,19 49,8
11 49,8 39,8 44,8188 276 2480 26,10 39,8
12 39,8 29,9 34,8590667 539 3019 31,77 29,9
13 29,9 19,9 24,8993333 1262 4281 45,05 19,9
14 19,9 10,0 14,9396 1529 5810 61,14 10,0
15 10,0 0,0 4,97986667 3693 9503 100,00 0,0
Suma 9503
Tabla 3. Datos de curva de regulación de caudales.

Grafico 1. Curva de duración de caudales diarios

A partir del Grafico 1 se determina el valor del volumen bajo la recta definida por el Q 50 y
el valor del volumen bajo la curva (VT). Con un valor de Q50=16.8534 se realiza el
respectivo calculo del valor de VP.

VP=Q 50∗50

VP=16.8534∗50=842.67 m3
Para calcular el VT se usa el método de los trapecios, en el cual se usa la siguiente ecuación
( B+b )∗h
A=
2
 Calculo VT
Rango medio Frec. Relativa acum (%) Área de trapecios
144,4161333 0,06 27,88
134,4564 0,26 42,24
124,4966667 0,59 69,17
114,5369333 1,17 183,31
104,5772 2,84 211,71
94,61746667 4,97 320,71
84,65773333 8,54 429,29
74,698 13,93 375,63
64,73826667 19,32 231,41
54,77853333 23,19 144,63
44,8188 26,10 225,96
34,85906667 31,77 396,80
24,89933333 45,05 320,50
14,9396 61,14 387,05
4,979866667 100,00 193,52
Área total 3559,80
Tabla 4. Datos para calcular VT

Con esto nos deja un VT con valor de 3559.80

Se procede a calcular el valor de IRH por el siguiente método
VP
IRH=
VT
842.67
IRH= =0.2
3559.8
IRH
VT 3559,80
VP 842,67
IRH 0,2
Tabla 5. Resumen cálculos IRH

Teniendo en cuenta el valor obtenido de IRH se puede decir que la corriente tiene retención
y regulación bajas, y que al ser menor a 0.7 el caudal ecológico requerido va a ser igual al
caudal correspondiente al Q 75 en la curva de duración de caudales.
 Grafico 2. Caudal ecológico

Teniendo en cuenta el Q 75, el valor del caudal ecológico es de 6.407 m3 /s

2.3 Pruebas de bondad y ajuste

Se realizan distintas funciones de distribuciones de probabilidad teórica y de estos se elige

la que mejor se ajusta a los datos que se tienen. Para el proyecto se usarán los métodos
Normal, Log-normal, Gumbel y Log-gumbel.

2.3.1 Prueba de bondad y ajuste para caudales máximos

[Link] Método Normal
 Caudales máximos
105,20
106,08
106,08
106,08
107,85
109,62
110,50
110,50
110,50
111,38
112,27
112,27
112,27
112,27
113,15
114,04
114,92
114,92
114,92
114,92
114,92
134,72
145,95
149,40
149,40
Tabla 6. Caudales máximos para función normal

Valor máximo 149,40

Valor mínimo 105,20
N° Datos 25
Rango 44,20
5,281
N° Intervalos
6
Tamaño intervalo 8,840
Promedio 116,56
Desv. Estandar 13,189
Tabla 7. Parámetros de cálculo de función normal para caudales máximos
 FRECUENCIA O = FRECUENCIA RELATIVA E = PROBABILIDAD POR LA FUNCIÓN ܱെ‫ܧ‬ ଶ
INTERVALOS FRECUENCIA
RELATIVA ACUMULADA NORMAL
ܱ
140,56 149,40 3 0,1200 0,1200 0,0281 0,070456
131,72 140,56 1 0,0400 0,1600 0,0909 0,029873
122,88 131,72 0 0,0000 0,1600 0,1908 0,005933
114,04 122,88 6 0,2400 0,4000 0,2599 0,049080
105,20 114,04 15 0,6000 1,0000 0,2296 0,593463
96,36 105,20 0 0,0000 1,0000 0,1316 0,754081
Suma 25 1,502887
Tabla 8. Prueba de bondad del ajuste normal

[Link] Método Log-normal

Caudales máximos LN
105,20 4,66
106,08 4,66
106,08 4,66
106,08 4,66
107,85 4,68
109,62 4,70
110,50 4,71
110,50 4,71
110,50 4,71
111,38 4,71
112,27 4,72
112,27 4,72
112,27 4,72
112,27 4,72
113,15 4,73
114,04 4,74
114,92 4,74
114,92 4,74
114,92 4,74
114,92 4,74
114,92 4,74
134,72 4,90
145,95 4,98
149,40 5,01
149,40 5,01
Tabla 9. Caudales máximos para función Log-normal
 Valor máximo 5,01
Valor mínimo 4,66
N° Datos 25
Rango 0,35
5,281
N° Intervalos
6
Tamaño intervalo 0,070
Promedio 4,75
Desv. Estandar 0,10
Tabla 10. Parámetros de cálculo de función Log-normal para caudales máximos

FRECUENCIA O = FRECUENCIA RELATIVA E = PROBABILIDAD POR LA FUNCIÓN ଶ

INTERVALOS FRECUENCIA ܱെ‫ܧ‬
RELATIVA ACUMULADA LOG NORMAL ܱ
4,94 5,01 3 0,1200 0,1200 0,0315 0,065323
4,87 4,94 1 0,0400 0,1600 0,0991 0,023211
4,80 4,87 0 0,0000 0,1600 0,2011 0,010553
4,73 4,80 7 0,2800 0,4400 0,2633 0,070926
4,66 4,73 14 0,5600 1,0000 0,2225 0,604474
4,59 4,66 0 0,0000 1,0000 0,1213 0,772103
Suma 25 1,546591
Tabla 11. Prueba de bondad del ajuste Log-normal

[Link] Método de Gumbel

Caudales máximos
105,20
106,08
106,08
106,08
107,85
109,62
110,50
110,50
110,50
111,38
112,27
112,27
112,27
112,27
113,15
114,04
114,92
114,92
114,92
114,92
114,92
134,72
145,95
149,40
149,40
 Tabla 12. Caudales máximos para función Gumbel

Valor máximo 149,40

Valor mínimo 105,20
N° Datos 25
Rango 44,20
5,281
N° Intervalos
6
Tamaño intervalo 8,840
Promedio 116,56
Desv. Estandar 13,189
Tabla 13. Parámetros de cálculo de función Gumbel para caudales máximos

FRECUENCIA O = FRECUENCIA RELATIVA E = PROBABILIDAD POR LA FUNCIÓN ܱ െ‫ܧ‬ ଶ

INTERVALOS FRECUENCIA
RELATIVA ACUMULADA GUMBEL
ܱ
140,56 149,40 3 0,1200 0,1200 0,039467 0,054046
131,72 140,56 1 0,0400 0,1600 0,076924 0,043135
122,88 131,72 0 0,0000 0,1600 0,139969 0,002508
114,04 122,88 6 0,2400 0,4000 0,221180 0,079941
105,20 114,04 15 0,6000 1,0000 0,262713 0,543592
96,36 105,20 0 0,0000 1,0000 0,177988 0,675704
Suma 25 1,398927
Tabla 14. Prueba de bondad del ajuste Gumbel

[Link] Método Log-Gumbel

 Caudales máximos LN
105,20 4,66
106,08 4,66
106,08 4,66
106,08 4,66
107,85 4,68
109,62 4,70
110,50 4,71
110,50 4,71
110,50 4,71
111,38 4,71
112,27 4,72
112,27 4,72
112,27 4,72
112,27 4,72
113,15 4,73
114,04 4,74
114,92 4,74
114,92 4,74
114,92 4,74
114,92 4,74
114,92 4,74
134,72 4,90
145,95 4,98
149,40 5,01
149,40 5,01

Tabla 15. Caudales máximos para función Log-Gumbel

Valor máximo 5,01

Valor mínimo 4,66
N° Datos 25
Rango 0,35
5,281
N° Intervalos
6
Tamaño intervalo 0,070
Promedio 4,75
Desv. Estandar 0,10
Tabla 16. Parámetros de cálculo de función Log-Gumbel para caudales máximos
 FRECUENCIA O = FRECUENCIA RELATIVA E = PROBABILIDAD POR LA FUNCIÓN ܱെ ‫ ܧ‬ଶ
INTERVALOS FRECUENCIA
RELATIVA ACUMULADA LOG GUMBEL
ܱ
4,94 5,01 3 0,1000 0,1000 0,041925 0,033727
4,87 4,94 1 0,0333 0,1333 0,081740 0,019964
4,80 4,87 0 0,0000 0,1333 0,147980 0,001609
4,73 4,80 7 0,2333 0,3667 0,229971 0,050961
4,66 4,73 14 0,4667 0,8333 0,262228 0,391394
4,59 4,66 0 0,0000 0,8333 0,162697 0,539704
Suma 25 1,037359
Tabla 17. Prueba de bondad del ajuste Log-Gumbel

2.3.2 Prueba de bondad y ajuste para caudales mínimos

[Link] Método Normal
Caudales mínimos
0,71
0,71
0,71
0,71
0,72
0,72
0,72
0,72
0,73
0,74
0,74
0,75
0,76
0,76
0,76
0,77
0,78
0,78
0,78
0,83
0,85
0,86
0,93
0,94
1,02
Tabla 18. Caudales mínimos para función normal
 Valor máximo 1,02
Valor mínimo 0,71
N° Datos 25
Rango 0,32
5,281
N° Intervalos
6
Tamaño intervalo 0,063
Promedio 0,78
Desv. Estandar 0,084
Tabla 19. Parámetros de cálculo de función normal para caudales mínimos

FRECUENCIA O = FRECUENCIA RELATIVA E = PROBABILIDAD POR LA FUNCIÓN ܱെ‫ܧ‬ ଶ

INTERVALOS FRECUENCIA
RELATIVA ACUMULADA NORMAL
ܱ
0,96 1,02 1 0,0400 0,0400 0,0134 0,017709
0,90 0,96 2 0,0800 0,1200 0,0637 0,026448
0,83 0,90 2 0,0800 0,2000 0,1759 0,002912
0,77 0,83 4 0,1600 0,3600 0,2826 0,016659
0,71 0,77 16 0,6400 1,0000 0,2642 0,541386
0,64 0,71 0 0,0000 1,0000 0,1438 0,733125
Suma 25 1,338238
Tabla 20. Prueba de bondad del ajuste normal

[Link] Método Log-normal

 Caudales mínimos LN
0,71 -0,35
0,71 -0,35
0,71 -0,35
0,71 -0,35
0,72 -0,34
0,72 -0,34
0,72 -0,34
0,72 -0,34
0,73 -0,32
0,74 -0,31
0,74 -0,31
0,75 -0,29
0,76 -0,28
0,76 -0,28
0,76 -0,27
0,77 -0,27
0,78 -0,25
0,78 -0,25
0,78 -0,25
0,83 -0,19
0,85 -0,17
0,86 -0,15
0,93 -0,07
0,94 -0,06
1,02 0,02
Tabla 21. Caudales mínimos para función Log-normal

Valor máximo 0,02

Valor mínimo -0,35
N° Datos 25
Rango 0,37
5,281
N° Intervalos
6
Tamaño intervalo 0,074
Promedio -0,26
Desv. Estandar 0,101
Tabla 22. Parámetros de cálculo de función Log-normal para caudales mínimos

INTERVALOS FRECUENCIA
FRECUENCIA O = FRECUENCIA RELATIVA E = PROBABILIDAD POR LA FUNCIÓN ܱെ ‫ ܧ‬ଶ
RELATIVA ACUMULADA LOG NORMAL ܱ
-0,05 0,02 1 0,0400 0,0400 0,0182 0,011895
-0,13 -0,05 2 0,0800 0,1200 0,0758 0,016288
-0,20 -0,13 3 0,1200 0,2400 0,1886 0,011001
-0,27 -0,20 5 0,2000 0,4400 0,2806 0,057724
-0,35 -0,27 10 0,4000 0,8400 0,2497 0,414764
-0,42 -0,35 4 0,1600 1,0000 0,1329 0,751825
Suma 25 1,263496
 Tabla 23. Prueba de bondad del ajuste Log-normal

[Link] Método de Gumbel

Caudales mínimos
0,71
0,71
0,71
0,71
0,72
0,72
0,72
0,72
0,73
0,74
0,74
0,75
0,76
0,76
0,76
0,77
0,78
0,78
0,78
0,83
0,85
0,86
0,93
0,94
1,02
Tabla 24. Caudales mínimos para función Gumbel
Valor máximo 1,02
Valor mínimo 0,71
N° Datos 25
Rango 0,32
5,281
N° Intervalos
6
Tamaño intervalo 0,063
Promedio 0,78
Desv. Estandar 0,084
Tabla 25. Parámetros de cálculo de función Gumbel para caudales mínimos
 FRECUENCIA O = FRECUENCIA RELATIVA E = PROBABILIDAD POR LA FUNCIÓN ܱെ‫ܧ‬ ଶ
INTERVALOS FRECUENCIA
RELATIVA ACUMULADA GUMBEL
ܱ
0,96 1,02 1 0,0400 0,0400 0,029772 0,002615
0,90 0,96 2 0,0800 0,1200 0,064408 0,025754
0,83 0,90 2 0,0800 0,2000 0,130698 0,024014
0,77 0,83 4 0,1600 0,3600 0,229723 0,047144
0,71 0,77 16 0,6400 1,0000 0,294208 0,498142
0,64 0,71 0 0,0000 1,0000 0,192887 0,651431
Suma 25 1,249101
Tabla 26. Prueba de bondad del ajuste Gumbel

[Link] Método de Log-Gumbel

Caudales mínimos LN
0,71 -0,35
0,71 -0,35
0,71 -0,35
0,71 -0,35
0,72 -0,34
0,72 -0,34
0,72 -0,34
0,72 -0,34
0,73 -0,32
0,74 -0,31
0,74 -0,31
0,75 -0,29
0,76 -0,28
0,76 -0,28
0,76 -0,27
0,77 -0,27
0,78 -0,25
0,78 -0,25
0,78 -0,25
0,83 -0,19
0,85 -0,17
0,86 -0,15
0,93 -0,07
0,94 -0,06
1,02 0,02
Tabla 27. Caudales mínimos para función Log-Gumbel
 Valor máximo 0,02
Valor mínimo -0,35
N° Datos 25
Rango 0,37
5,281
N° Intervalos
6
Tamaño intervalo 0,074
Promedio -0,26
Desv. Estandar 0,101
Tabla 28. Parámetros de cálculo de función Log-Gumbel para caudales mínimos

FRECUENCIA O = FRECUENCIA RELATIVA E = PROBABILIDAD POR LA FUNCIÓN ܱെ‫ܧ‬ ଶ

INTERVALOS FRECUENCIA
RELATIVA ACUMULADA LOG GUMBEL
ܱ
-0,05 0,02 1 0,0333 0,0333 0,033618 0,000002
-0,13 -0,05 2 0,0667 0,1000 0,070717 0,008575
-0,20 -0,13 3 0,1000 0,2000 0,138864 0,018688
-0,27 -0,20 5 0,1667 0,3667 0,234464 0,047666
-0,35 -0,27 10 0,3333 0,7000 0,285924 0,244941
-0,42 -0,35 4 0,1333 0,8333 0,177732 0,515776
Suma 25 0,835649
Tabla 29. Prueba de bondad del ajuste Log-Gumbel

2.3.3 Prueba de bondad y ajuste para caudales medios

[Link] Método Normal
 Caudales medios
26,03
26,37
26,39
26,45
26,51
26,52
26,52
26,54
26,59
26,59
26,66
26,80
26,80
26,87
26,87
26,88
26,91
26,94
27,07
27,23
27,29
33,01
33,40
33,69
33,85
Tabla 30. Caudales medios para función Normal

Valor máximo 33,85

Valor mínimo 26,03
N° Datos 25
Rango 7,82
5,281
N° Intervalos
6
Tamaño intervalo 1,564
Promedio 27,79
Desv. Estandar 2,56
Tabla 31. Parámetros de cálculo de función normal para caudales medios

FRECUENCIA O = FRECUENCIA RELATIVA E = PROBABILIDAD POR LA ܱെ‫ܧ‬ ଶ

INTERVALOS FRECUENCIA
RELATIVA ACUMULADA FUNCIÓN NORMAL
ܱ
32,29 33,85 4 0,1600 0,1600 0,0304 0,104929
30,72 32,29 0 0,0000 0,1600 0,0864 0,033834
29,16 30,72 0 0,0000 0,1600 0,1707 0,000710
27,59 29,16 0 0,0000 0,1600 0,2344 0,034573
26,03 27,59 20 0,8000 0,9600 0,2239 0,564441
24,46 26,03 1 0,0400 1,0000 0,1488 0,724618
Suma 25 1,463104
Tabla 32. Prueba de bondad del ajuste normal
[Link] Método Log-Normal

Caudales medios LN
26,03 3,26
26,37 3,27
26,39 3,27
26,45 3,28
26,51 3,28
26,52 3,28
26,52 3,28
26,54 3,28
26,59 3,28
26,59 3,28
26,66 3,28
26,80 3,29
26,80 3,29
26,87 3,29
26,87 3,29
26,88 3,29
26,91 3,29
26,94 3,29
27,07 3,30
27,23 3,30
27,29 3,31
33,01 3,50
33,40 3,51
33,69 3,52
33,85 3,52
Tabla 33. Caudales medios para función Log-normal

Valor máximo 3,52

Valor mínimo 3,26
N° Datos 25
Rango 0,26
5,281
N° Intervalos
6
Tamaño intervalo 0,053
Promedio 3,32
Desv. Estandar 0,085
Tabla 34. Parámetros de cálculo de función Log-normal para caudales medios
 FRECUENCIA O = FRECUENCIA RELATIVA E = PROBABILIDAD POR LA ܱെ‫ܧ‬ ଶ
INTERVALOS FRECUENCIA
RELATIVA ACUMULADA FUNCIÓN LOG NORMAL
ܱ
3,47 3,52 4 0,1600 0,1600 0,0318 0,102664
3,42 3,47 0 0,0000 0,1600 0,0897 0,030851
3,36 3,42 0 0,0000 0,1600 0,1752 0,001449
3,31 3,36 0 0,0000 0,1600 0,2370 0,037099
3,26 3,31 20 0,8000 0,9600 0,2222 0,567014
3,21 3,26 1 0,0400 1,0000 0,1443 0,732157
Suma 25 1,471234
Tabla 35. Prueba de bondad del ajuste normal

[Link] Método de Gumbel

Caudales medios
26,03
26,37
26,39
26,45
26,51
26,52
26,52
26,54
26,59
26,59
26,66
26,80
26,80
26,87
26,87
26,88
26,91
26,94
27,07
27,23
27,29
33,01
33,40
33,69
33,85
Tabla 36. Caudales medios para función Gumbel
 Valor máximo 33,85
Valor mínimo 26,03
N° Datos 25
Rango 7,82
5,281
N° Intervalos
6
Tamaño intervalo 1,564
Promedio 27,79
Desv. Estandar 2,56
Tabla 37. Parámetros de cálculo de función Gumbel para caudales medios

FRECUENCIA O = FRECUENCIA RELATIVA E = PROBABILIDAD POR LA ܱെ‫ܧ‬ ଶ

INTERVALOS FRECUENCIA
RELATIVA ACUMULADA FUNCIÓN GUMBEL
ܱ
32,29 33,85 4 0,1600 0,1600 0,039370 0,090948
30,72 32,29 0 0,0000 0,1600 0,072185 0,048197
29,16 30,72 0 0,0000 0,1600 0,124851 0,007722
27,59 29,16 0 0,0000 0,1600 0,192928 0,006777
26,03 27,59 20 0,8000 0,9600 0,240270 0,539596
24,46 26,03 1 0,0400 1,0000 0,199422 0,640925
Suma 25 1,334164
Tabla 38. Prueba de bondad del ajuste Gumbel

[Link] Método de Log-Gumbel

 Caudales medios LN
26,03 3,26
26,37 3,27
26,39 3,27
26,45 3,28
26,51 3,28
26,52 3,28
26,52 3,28
26,54 3,28
26,59 3,28
26,59 3,28
26,66 3,28
26,80 3,29
26,80 3,29
26,87 3,29
26,87 3,29
26,88 3,29
26,91 3,29
26,94 3,29
27,07 3,30
27,23 3,30
27,29 3,31
33,01 3,50
33,40 3,51
33,69 3,52
33,85 3,52
Tabla 39. Caudales medios para función Log-Gumbel

Valor máximo 3,52

Valor mínimo 3,26
N° Datos 25
Rango 0,26
5,281
N° Intervalos
6
Tamaño intervalo 0,053
Promedio 3,32
Desv. Estandar 0,085
Tabla 40. Parámetros de cálculo de función Log-Gumbel para caudales medios

FRECUENCIA O = FRECUENCIA RELATIVA E = PROBABILIDAD POR LA ܱെ‫ܧ‬ ଶ

INTERVALOS FRECUENCIA
RELATIVA ACUMULADA FUNCIÓN LOG GUMBEL
ܱ
3,47 3,52 4 0,1333 0,1333 0,040392 0,064786
3,42 3,47 0 0,0000 0,1333 0,074188 0,026237
3,36 3,42 0 0,0000 0,1333 0,128261 0,000193
3,31 3,36 0 0,0000 0,1333 0,197216 0,030607
3,26 3,31 20 0,6667 0,8000 0,242128 0,389026
3,21 3,26 1 0,0333 0,8333 0,194540 0,489669
Suma 25 1,000518
 Tabla 41. Prueba de bondad del ajuste Log-Gumbel

2.4 Distribución de probabilidad

El método que más se acomodo fue la distribución de Log-Gumbel para todos los caudales.

2.4.1 Caudales máximos

CAUDAL MAXIMO DIARIO ANUAL LN

Año Qmax Log(Qmax)

1995 149,396 5,007
1996 114,92 4,744
1997 112,268 4,721
1998 114,036 4,737
1999 114,92 4,744
2000 105,196 4,656
2001 110,5 4,705
2002 112,268 4,721
2003 113,152 4,729
2004 106,08 4,664
2005 114,92 4,744
2006 110,5 4,705
2007 114,92 4,744
2008 112,268 4,721
2009 109,616 4,697
2010 114,92 4,744
2011 106,08 4,664
2012 106,08 4,664
2013 110,5 4,705
2014 112,268 4,721
2015 107,848 4,681
2016 111,384 4,713
2017 149,396 5,007
2018 134,722 4,903
2019 145,948 4,983
Tabla 42. Caudales máximos diarios anuales

Parámetros del ajuste

µ 4,70
α 10,71
Tabla 43. Parámetros estadísticos función Log-Gumbel
De acuerdo con lo anterior, la función de densidad de probabilidad se define como:
0.475 ( xi −5.259)

F ( x ≤ x i )=e−e

Se definen los periodos de retorno 2, 5, 10, 20, 25, 50 y 100 años. Se calcula la
probabilidad de ocurrencia, de excedencia y se calcula la variable estándar para el cálculo
del caudal

LOG GUMBEL
TR F(x ≥ xi) F(x ≤ xi) Ln(xi) xi
2 0,5 0,5 4,74 114,16
5 0,2 0,8 4,84 126,90
10 0,1 0,9 4,91 136,11
15 0,1 0,9 4,95 141,60
20 0,05 0,95 4,98 145,57
25 0,04 0,96 5,00 148,71
50 0,02 0,98 5,07 158,80
100 0,01 0,99 5,13 169,49
Tabla 44. Aplicación de distribución de Log-Gumbel

2.4.2 Caudales mínimos

 CAUDAL M INIM O DIARIO ANUAL LN

Año Qm in Log(Qm in)

1995 0,93 -0,073
1996 0,715 -0,335
1997 0,706 -0,348
1998 0,706 -0,348
1999 0,776 -0,254
2000 0,715 -0,335
2001 0,736 -0,307
2002 0,715 -0,335
2003 0,706 -0,348
2004 0,706 -0,348
2005 0,847 -0,166
2006 0,736 -0,307
2007 0,776 -0,254
2008 0,756 -0,280
2009 0,746 -0,293
2010 0,766 -0,267
2011 0,776 -0,254
2012 0,764 -0,269
2013 0,756 -0,280
2014 0,827 -0,190
2015 0,726 -0,320
2016 0,715 -0,335
2017 0,858 -0,153
2018 1,022 0,022
2019 0,944 -0,058
Tabla 45. Caudales mínimos diarios anuales

Parámetros del ajuste

µ -0,305
α 11,040
Tabla 46. Parámetros estadísticos función Log-Gumbel

De acuerdo con lo anterior, la función de densidad de probabilidad se define como:

0.475 ( xi −5.259)

F ( x ≤ x i )=e−e

Se definen los periodos de retorno 2, 5, 10, 20, 25, 50 y 100 años. Se calcula la
probabilidad de ocurrencia, de excedencia y se calcula la variable estándar para el cálculo
del caudal
LOG GUMBEL
TR F(x ≥ xi) F(x ≤ xi) Ln(xi) xi
2 0,5 0,5 -0,27 0,76
5 0,2 0,8 -0,17 0,84
10 0,1 0,9 -0,10 0,90
15 0,1 0,9 -0,06 0,94
20 0,05 0,95 -0,04 0,96
25 0,04 0,96 -0,02 0,98
50 0,02 0,98 0,05 1,05
100 0,01 0,99 0,11 1,12
 Tabla 47. Aplicación de distribución de Log-Gumbel

2.4.3 Caudales medios

CAUDAL MEDIO DIARIO ANUAL LN

Año Qmed Log(Qmed)

1995 33,397 3,508
1996 26,803 3,289
1997 26,450 3,275
1998 27,068 3,298
1999 27,288 3,306
2000 26,803 3,289
2001 26,868 3,291
2002 26,542 3,279
2003 26,515 3,278
2004 26,594 3,281
2005 26,658 3,283
2006 26,591 3,281
2007 26,868 3,291
2008 26,936 3,293
2009 26,517 3,278
2010 26,909 3,292
2011 26,512 3,278
2012 26,366 3,272
2013 26,885 3,292
2014 27,233 3,304
2015 26,029 3,259
2016 26,386 3,273
2017 33,851 3,522
2018 33,693 3,517
2019 33,010 3,497
Tabla 48. Caudales medios diarios anuales

Parámetros del ajuste

µ 3,28032
α 13,04408
Tabla 49. Parámetros estadísticos función Log-Gumbel

De acuerdo con lo anterior, la función de densidad de probabilidad se define como:

0.475 ( xi −5.259)

F ( x ≤ x i )=e−e
Se definen los periodos de retorno 2, 5, 10, 20, 25, 50 y 100 años. Se calcula la
probabilidad de ocurrencia, de excedencia y se calcula la variable estándar para el cálculo
del caudal
LOG GUMBEL
TR F(x ≥ xi) F(x ≤ xi) Ln(xi) xi
2 0,5 0,5 3,31 27,34
5 0,2 0,8 3,40 29,82
10 0,1 0,9 3,45 31,59
15 0,1 0,9 3,49 32,63
20 0,05 0,95 3,51 33,38
25 0,04 0,96 3,53 33,97
50 0,02 0,98 3,58 35,85
100 0,01 0,99 3,63 37,83
Tabla 50. Aplicación de distribución de Log-Gumbel

Se decide trabajar un periodo de retorno de 100 años, por lo cual se procede a ajustar los
caudales para ser usados en el diseño

CAUDALES DE DISEÑO PARA 100 AÑOS DE PERIODO DE RETORNO

Qmax [m 3/s] Qmin [m 3/s] Qmed [m 3/s]
16,95 0,11 3,78
Tabla 51. Caudales de diseño

3. Traslado de los caudales al punto de bocatoma

3.1 Zona de captación
A continuación, se presenta una ilustración de la cuenca en dónde va a estar ubicado
nuestro proyecto y la Bocatoma:
 Ilustracion 1. Zona de captacion

Teniendo en cuenta la topografía e hidrología presentes en la cuenca, se optó por la zona de

bocatoma tal cómo se puede apreciar en la ilustración siguiente, y en la cual se presenta un
gran potencial hidrológico, ya que ventaja, está la altitud del lugar, puesto que la Bocatoma
está ubicada en la cota 819.00, la cual nos proporciona unas diferencias de alturas con
respecto a las cotas dónde se encuentra ubicado el municipio de interés y una pendiente de
solera en condiciones convenientes para la instalación del mismo, en esta situación lo que
facilitará en parte la eficiencia del sistema y los costos del transporte del recurso hídrico
hacia los tanques de abastecimiento.
 Ilustracion 2. Zona de bocatoma

3.2 Traslado de caudales

Una vez definido el sitio de captación, es necesario conocer cuáles serán los caudales que
se presentarán en dicho punto de Bocatoma. Para el traslado de los caudales se utilizó el
método de rendimiento de la cuenca y las microcuencas de la sección estudiada y por donde
esta alimenta, este método permite representar la escorrentía por unidad de área. Es
necesario realizar un traslado de caudales teniendo en cuenta los datos conocidos por la
estación hidrológica ubicada aguas abajo del punto de interés sabiendo que los caudales
obtenidos en dicha estación fueron:

CAUDALES DE DISEÑO PARA 100 AÑOS DE PERIODO DE RETORNO

Qmax [m 3/s] Qmin [m 3/s] Qmed [m 3/s]
16,95 0,11 3,78

Tabla 52. Traslados de Caudales

De acuerdo con lo supuesto del método este se da inicio según hipótesis de:

Qmedio m3
Rmedio = (
Acuenca s
/Km 2 )
Para el traslado de los caudales medios se establece la siguiente relación:

Restacion ≅ Rinteres

De la anterior relación obtenemos lo siguiente:

Qestacion∗A interes
Qinteres=
A estacion

Teniendo en cuenta el área de aporte de las subcuentas a el sitio de bocatoma, y el área total
de la cuenca que aporta el agua a la estación hidrológica. Se obtuvieron los siguientes
resultados
Microcuenca Área (m²)
Área cuenca (m²) 941099536.97
1 1154448.00
2 3199536.00
3 797328.00
4 614592.00
5 1254384.00
6 1860336.00
7 4914864.00
8 1258129.00
9 847440.00
10 1075824.00
11 952272.00
12 1462176.00
13 723024.00
14 547200.00
15 1515024.00
16 227664.00
17 1672848.00
18 1935432.00
Subcuenca 26012521.00
 Tabla 53

Dichas áreas se muestran a continuación con el objetivo de brindar mayor claridad a la

descripción dada anteriormente:

Ilustracion 3.

Teniendo claramente definidas las áreas y los caudales de diseño obtenidos de la estación
de medición, se procede a hacer el cálculo para el traslado de los caudales como se mostró
en la anterior ecuación descrita

A subcuencu =26012521.00m2
Acuenca =941099536.97 m2
16.95∗26012521.00
Q max =
941099536.97

m3
Q max =0.47
s
Repitiendo este mismo paso con cada uno de los caudales de diseño para hacer uso en el
diseño de este de la bocatoma se obtuvo:
CAUDALES DE DISEÑO PARA 100 AÑOS DE PERIODO DE RETORNO
3 3 3
Qmax [m /s] Qmin [m /s] Qmed [m /s]
0.47 0.0030892 0.1046

Tabla 54

4. Modelación en Hec-Ras del tramo del río

Mediante la modelación del sistema el en programa HEC-RAS es posibles realizar el

diseño de la bocatoma y escoger la sección más adecuada del tramo del rio optada para
poder ubicar pertinentemente la zona de captación, con el fin de obtener datos del flujo en
la sección escogida para la zona de captación utilizando las curvas de nivel que sigue el
trayecto del rio. Una vez importados los datos a HEC-RAS se procede a realizar la
modelación del flujo en el trayecto del rio seleccionado
 Gráfico 3. Tramo del rio

Para la modelación del ancho de banca se asumen los valores de 0.027, coeficiente de
rugosidad correspondiente a terreno con pastos cortos y con pocas malezas, y para el canal
un coeficiente de rugosidad de 0.017 correspondiente a concreto sin pulir de acuerdo con la
tabla 5.6 del libro Open Chanel Flow Hydraulics de Ven Te Chow. Es importante tener en
cuenta que el ancho de banca se ubicó en un proceso iterativo, esto se realizó corriendo el
programa varias veces y modificándolo con el fin de obtener resultados más precisos.
Una vez se tenga la información de las secciones transversales se introducen los valores
correspondientes a los caudales transpuestos calculados anteriormente como se muestra en
la siguiente ilustración:
Finalmente se obtienen las vistas en dos dimensiones para los caudales mínimo, medio y
máximo como se muestran en las siguientes imágenes:

Gráfico 4. Caudal mínimo

Gráfico 5. Caudal medio

 Gráfico 6. Caudal máximo

A partir de las anteriores ilustraciones y de los datos obtenidos mediante HEC RAS, se
escoge la sección 90 debido al comportamiento del flujo, la pendiente de la solera y las
dimensiones de esta.

Gráfico 7. Sección
 Tabla 55

5. Diseño del Dique

5.1 Diseño del dique WES

Para el diseño de las dimensiones del dique se tomó Lo=7 m y h=2.8 m. Con estos datos ya
selectos se pueden definir la carga del agua bajo condiciones máximas y mínimas en la
cresta del dique, también se calculan la velocidad y el valor de la energía en el punto
elegido para el análisis.
h
Se tomo que >1.33 para que el valor de H e =H d , con un valor de C=2.225 m−1 /s
Hd

Q max 23
Hd= ( C∗Lo ) =1.058 m=H e

h
=2.65>1.33
Hd

Como se cumple el parámetro no se debe tener en cuenta la corrección de H v

2
Q max
H v=
(T∗( h+ H d )∗2∗g )
=0.02m
 Ilustración 4. Factor de corrección

Debido a que no se hace corrección se tomo como un dique vertical, por lo cual es de tipo 1
donde no hay talud de ningún tipo.

K 2,00
n 1,850
r1 0,529 m
r2 0,212 m
M 0,185 m
N 0,298 m
Tabla 56. Parámetros de diseño dique WES

Usando los valores de la tabla 52 se puede graficar el comportamiento del vertedero por
medio de la siguiente ecuación

xn
y=
K H n−1
d
 x(m) y(m)
0 0,000
0,1 0,007
0,2 0,024
0,3 0,051
0,4 0,087
0,5 0,132
0,6 0,185
0,7 0,246
0,8 0,315
0,9 0,392
1 0,477
1,1 0,568
1,2 0,668
1,3 0,774
1,4 0,888
1,5 1,009
1,6 1,137
1,7 1,272
1,8 1,414
1,9 1,563
2 1,718
2,1 1,880
2,2 2,049
2,3 2,225
2,4 2,407
2,5 2,596
Tabla 57. Datos de perfil dique WES

Gráfico 8. Vertedero Dique WES

Para la salida del vertedero se verifica el resalto hidráulico que se presenta en el cuenco
amortiguador sea estable, se asume ∆ B=0.03 m
Se realiza el cálculo de la altura y la velocidad por medio de las siguientes ecuaciones
Z=h+ H o + ∆ B

V 1= √2∗g∗( Z−0.5∗H d )

∆B 0,03 m
Z 3,88802135 m
V1 8,11645966 m/s
Y1 0,29833556 m
FR1 4,74534919 Cumple
Tabla 58. Valores de resalto dique WES

No se consideró hacer iteraciones ya que el valor del número de Froude nos indica que el
resalto hidráulico es estable.
Una vez establecido que el resalto es estable se procede a calcular la altura del obstáculo
para poder determinar la lamina de agua que hay después del resalto, esto se hace por
medio de la siguiente ecuación

(√ 1+ 8 NF 2−1− 2yh)
21.33 NF 2 1
=
√1+8 NF 2−1
( √1+8 NF −1+ 2yh )
2

Con la cual se obtiene un valor de 0.58m

Se procede a hallar el valor de Y 2 con la siguiente ecuación

( )(
3
y y2 h
2
266 NF 1+
y2
i
= −
y1 y1 )
y1

El cual da un valor de 2.25m

Se plantea la ecuación de fuerza específica para hallar el valor de y 3que es el valor que se
encuentra justo sobre encima del vertedero y el y 4 que se calcula planteando la energía
especifica entre estos dos últimos puntos.
Por ultimo se plantea conocer la longitud del cuenco y para esto se hace uso de las
siguientes ecuaciones
x=5 ( H d − y 4 )

x=6 y 2

x=6.9 ( y 2− y1 )

Y se calcula el radio de curvatura con la siguiente expresión

V +6.4 H +16
3.6H +64
R=10

hobs 0,57964722 m
Y2 2,25389016 m
Yc 0,84246231 m
Y3 1,38175142 m
Y4 1,88824991 m
x 12,3394857 m
x 13,523341 m
x 13,5691209 m
X 14 m
R 2,31748 m
Tabla 59. Resumen diseño de resalto dique WES

Para la longitud del cuenco se adopta un valor de 14m ya que con este podemos tener una
certeza que el resalto se producirá en esa zona.

5.2 Diseño de Dique Ogee

Para el diseño del dique tipo Ogee se usan los mismos valores que se tomaron para el
diseño de dique WES, como lo son caudales y alturas de diseño.
El primer paso es calcular el coeficiente de descarga para esto se adopta un valor entre 3,0 y
4,0 para empezar la iteración.
Se toma un valor de C=3.8 y se procede a calcular la carga de diseño sobre el vertedero
con la siguiente ecuación
2
Q
( )
H o = max
C∗L
3

H o =0.741m
Y con este valor se realiza la relación P/ H o y se usa el grafico 4 para hallar el coeficiente
que corresponde con los parámetros de diseño.
P 2.8
= =3.78
H o 0.741

Gráfico 9. Coeficiente de descarga

Como el primer valor de la relación P/ H o ≥3 se toma el valor del coeficiente de descarga

como 3.95 y se hace el cálculo de H o nuevamente al igual que la relación P/ H o

Lo cual deja un valor para H o =0.722 y P/ H o=3.88 , teniendo en cuenta estos parámetros
se tiene un vertedero vertical y un C corregido =3.95
Se realiza el proceso de cálculo del coeficiente de descarga, pero esta vez para el caudal
mínimo ( Q min=0.11 m3 /s )

En este caso se elige un valor de C entre 3,0 y 4,0 para iniciar el proceso, se toma un valor
de C=3,8 y se procede a calcular H e de la siguiente manera

2
Q max
H e= ( )
C∗L
3

H e =0.026

Se realiza la relación H e / H oy su valor se utiliza en el grafico 5 para obtener C /C o

He
=0.0357
Ho

Gráfico 10. Coeficiente de descarga para carga diferentes

Debido a que la relación H e / H o da un valor muy pequeño, se toma C /C o =0.8 y se

procede a calcular el Coeficiente de descarga para el diseño
C
C= ∗C corregido =3.16
Co

Con el valor de C recalculado se procede a hacer el mismo proceso antes mencionado y se

obtiene un valor de coeficiente de descarga para el caudal mínimo de 3.24.

Se realiza el mismo proceso con el caudal medio ( Q med=3.78 m3 / s ) y su valor de coeficiente

de descarga es de 3.55
Una vez obtenido el valor del coeficiente de descarga se busca calcular H a ,que será
necesario para calcular los radios y la altura critica, para esto se usa la siguiente expresión
2
Q max

H a=
( )
L
2
( h+ H o ) 2 g
El valor obtenido se resta al H oy se calcula el valor de h o y se realiza una iteración hasta
que el valor de H ase mantenga constante

El valor de H aes igual a 0.0244, con este valor se realiza la relación H a / H opara usarlo en el
grafico 6 para calcular los valores restantes del diseño
Ha
=0.033
Ho
Gráfico 11. Coeficientes de diseño Dique Ogee
Se toman lo valores que correspondan al diseño de dique vertical
Xc
=0.167
Ho

Lo cual nos deja un valor de X c =0.121 m

Yc
=0.113
Ho

Lo cual nos deja un valor de Y c =0.082 m

R1
=0.515
Ho

Lo cual nos deja un valor de R1=0.372m

R2
=0.215
Ho

Lo cual nos deja un valor de R2=0.155 m

Se usa el mismo valor de H a / H oen el grafico 7 para obtener el valor de K y n

 Gráfico 12. Coeficiente de perfil del vertedero estricto
Del grafico 7 se obtiene que K=0.508 y n=1.855
ha 0,0244
Ha/Ho 0,0339
Xc/Ho 0,167
Xc 0,121
Yc/Ho 0,113
Yc 0,082
R1/Ho 0,515
R1 0,372
R2/Ho 0,215
R2 0,155
Lectura K 0,508
Lectura n 1,855
Tabla 60. Resumen de parámetros de diseño dique Ogee

Se procede a graficar el comportamiento del perfil del dique y se usa la siguiente ecuación

x n
y=−K ( )
Ho
Ho
 x(m) y(m)
0 0
0,1 -0,0094
0,2 -0,0339
0,3 -0,0720
0,4 -0,1227
0,5 -0,1856
0,6 -0,2603
0,7 -0,3465
0,8 -0,4438
0,9 -0,5522
1 -0,6714
1,1 -0,8013
1,2 -0,9416
1,3 -1,0924
1,4 -1,2533
1,5 -1,4245
1,6 -1,6056
1,7 -1,7967
1,8 -1,9977
1,9 -2,2085
2 -2,4289
2,1 -2,6590
2,2 -2,8987
2,3 -3,1478
2,4 -3,4064
2,5 -3,6744
Tabla 61. Datos perfil dique Ogee

Gráfico 13. Perfil dique Ogee

Se realiza el mismo proceso para el cálculo de la longitud del cuenco que se usa en el
diseño de dique tipo WES
∆b 0,03 m
Z 3,552 m
V1 7,912257555 m/s
Y1 0,306035105 m
FR1 4,566466633 Cumple
Tabla 62. Valores de resalto dique Ogee

hobs 0,553604836 m
Y2 2,201821205 m
Yc 0,842347787 m
Y3 1,337504461 m
Y4 1,814758723 m
x 11,8418178 m
x 13,21092723 m
x 13,09969594 m
X 14 m
R 2,326062464 m
Tabla 63. Resumen diseño de resalto dique Ogee

Como se puede observar el resalto es estable y el valor de la longitud del cuenco y los
demás valores calculados son muy similares a los del diseño tipo WES.

6. Diseño de Bocatoma
6.1 Bocatoma de orificio

Se realiza el diseño de la bocatoma de orificio y para esto se hace uso del QMD y la
viscosidad cinemática del agua a una temperatura de 18°C y para el taponamiento por
solidos transportados en la corriente y las perdidas en el sistema se amplia el caudal de
diseño en un porcentaje de 100%.
Se realiza el calculo del caudal derivado por medio de la siguiente expresión

ampliación(%) m3
Qd =QMD 1+ ( 100 )
=0.0928
s

Se procede a seleccionar el valor de la altura del orificio (a) y la carga sobre el centroide del
orificio (H). Los valores adoptados son a=0.15 m y H=0.5 m, se aconseja tomar el valor de
H como mínimo de 2.5a.
Para conocer el tipo de orificio y condiciones de operación de calcula la relación H /a, de la
cual se obtiene un valor de 3.33 lo cual significa que será un orificio de grandes
dimensiones.
Se realiza el calculo del número de Reynolds con el fin de obtener el coeficiente de
descarga del grafico 9, este se calcula por medio de la siguiente ecuación
a
ℜ=√ 2 gH
v
Con la cual se obtiene un valor de ℜ=4.45 x 10 5

Gráfico 14. Coeficientes en termino de Re

Del grafico 9 se obtiene C d=0.6
Se define el ancho y la separación entre las barras que serán usadas para evitar la entrada de
solidos a la cámara de captación, se toma una separación(b) de 1 pulgada y un ancho(s) de
0.5 pulgadas.
Con los valores definidos anteriormente se calcula el ancho de la estructura(B) se la
siguiente manera
Qd
B=
( b+b s )√ 2 gH 1− 96 1a
Cd a
( (H)
2

)
Con el cual se obtiene un valor B de 0.494m
Se realiza el calculo de la cantidad de barras(N) que son necesarias para el orificio con la
siguiente expresión
B−b
N=
b+ s
Se obtiene N=12.29 por lo cual se toman 13 barras y 14 espacios.
Con el número de barras se recalcula el valor del ancho de la estructura(B) con el número
de barras el cual se obtiene que es B=0.521 m y se calcula el ancho neto con el número de
espacios dando un valor de Bneto =0.356 m
Con el valor del ancho neto se procede a verificar los valores del número de Reynolds y
con este el coeficiente de descarga.
Qd a
ℜ=
Bneto ∗a v()
El número de Reynolds verificado tienen un valor de 2.47 x 10 5 y buscando este valor en el
grafico 9 obtenemos un valor para el coeficiente de descarga de 0.6, por lo cual se concluye
que queda verificado y no es necesario recalcular el ancho B.

6.1.1 Cámara de captación

Se realiza el diseño de la cámara de captación teniendo en cuenta que el ancho mínimo

debe ser de 1.2m ya que con esta distancia se tiene espacio suficiente para que una persona
entre y pueda realizar el mantenimiento.
Se realiza el calculo de las condiciones de nivel, carga, caudal captado y la velocidad inicial
con el fin de obtener el trayecto del chorro de agua en dos condiciones de flujo, es decir,
para caudales máximo y mínimos.

Ancho de rio Flujo máx Flujo mín

Nivel 3,8580 2,8368
Carga 1,5212 0,5
Q Cap (m3/s) 0,1707 0,0979
V inicial (m/s) 3,1995 1,8360
Tabla 64. Diseño de cámara de captación

Se realiza el calculo de la trayectoria mediante la siguiente formula

−g x 2
y=
( )
2 V 2o
Cambiando en la velocidad inicial ( V 2o ) por el valor de cada tipo de flujo.
X(m) Ymáx (m) Ymín (m)
0,000 0 0
0,100 -0,005 -0,015
0,200 -0,019 -0,058
0,300 -0,043 -0,131
0,400 -0,077 -0,233
0,500 -0,120 -0,364
0,600 -0,172 -0,524
0,700 -0,235 -0,713
0,800 -0,307 -0,931
0,900 -0,388 -1,178
1,000 -0,479 -1,454
1,100 -0,580 -1,760
1,200 -0,690 -2,094
1,300 -0,809 -2,458
1,327 -0,844 -2,563
1,500 -1,078 -3,273
1,600 -1,226 -3,723
1,700 -1,384 -4,203
1,800 -1,552 -4,713
1,900 -1,729 -5,251
2,000 -1,916 -5,818
2,100 -2,112 -6,414
2,200 -2,318 -7,040
Tabla 65. Trayecto de chorro para flujos máximos y mínimos
 Gráfico 15. Trayectorias de chorro

6.1.2 Vertedero

Como paso final del diseño de la bocatoma se realiza el diseño del vertedero, para esto se
realiza el calculo del caudal de excesos y se debe adoptar el valor del ancho del vertedero.
El caudal de excesos se calcula restando el caudal derivado ( Q d )al caudal captado ( Qc ) el
cual nos da un valor de 0.078 m 3 /s y se adopta un valor de ancho de 0.51m
Se realiza el calculo de H por medio de la ecuación de Francis
2
Qc

( )
3
1.84
H=
L

Del cual se obtiene un valor H=0.190 m

Se hace el cálculo de la aducción para estos de adopta un valor de 1 m/s para la velocidad y
por medio de esto se halla el valor del diámetro para buscar el diámetro comercial que más
se acomode.

Qd

D= 4
√(
π
V
)
De la expresión anterior se obtiene un valor de 0.34m por lo cual se adopta un valor
comercial de 0.35m o mejor conocidas como de 14”. Se hace el cálculo de la sumergencia
(S) tomando que la sumergencia es 2.5 veces el diámetro por lo cual S=0.875 m.
Se hace el calculo de la altura usando la siguiente expresión
D
h= + S+0.25
2
Dando un valor de h=1.3 m
Teniendo el valor de la altura de excesos(h) se realiza la relación H /hde la cual se obtiene
un valor de 0.146, esta relación debe cumplir el parámetro H /h ≪10 por lo cual cumple y
se puede usar la ecuación de Rouse para determinar el coeficiente C.

C=3.27+0.4 ( Hh )
De la cual se obtiene un valor para el coeficiente de 3.33 que se pasa a S.I.
C ( S . I ) =1.84
Se puede considerar que el vertedero tendrá doble contracción por los parámetros
adoptados para el vertedero. A partir de esta contracción se deben calcular los valores de la
longitud efectiva( Le ), la carga(H) y longitud real del sistema( L¿¿ r )¿

Le =( L−contracción ) 0.1∗H
2
Qe

( )
3

C (S . I )
H=
Le

Lr =Le +0.1 Contracciones∗H

Le 0,472 m
H 0,200 m
L 0,512 m
Tabla 66. Resumen longitud real

Se debe determinar como es la descarga hacia la cámara de excesos, para esto se usarán las
ecuaciones de USBR
Qe
V=
Le∗H
 V2
h v=
2g
hv
A=−0.425+0.25
H

hv hv 2 hv
B=0.411−1.603 − 1.568
H H √ ( )
−0.892 + 0.127
H
hv
C=0.150−0.45
H
hv
m= −0.208
H

D=0.57−0.02m2 em
V(m/s) 0,82255224
hv(m) 0,0345
A -0,38197192
B -0,00639313
C 0,07254945
m -0,03588766
D 0,56820087
Tabla 67. Descarga hacia la cámara de excesos

Se realiza la relación h v /H para comprobar que h v /H >2 , el valor obtenido de esta relación
es de 1.72 lo cual nos indica que las ecuaciones de USBR no son las más aptas, pero como
no se tienen otras ecuaciones para representar la trayectoria del vertedero se usan teniendo
algo de incertidumbre.
Se realiza la trayectoria del vertedero usando la siguiente ecuación

x x
(( ) ( )
y= A
H
+B
H
+C + D H
)
 X(m) Y(m)
0 0,128
0,1 0,108
0,2 0,050
0,3 -0,046
0,4 -0,180
0,5 -0,352
0,6 -0,562
0,7 -0,811
0,8 -1,097
0,9 -1,422
1 -1,785
1,1 -2,186
1,2 -2,625
1,3 -3,102
1,4 -3,618
1,5 -4,161
Tabla 68. Datos de la trayectoria del vertedero

Gráfico 16. Trayectoria del vertedero

Se adopta una distancia de sobreancho de 1.7m para evitar erosión como se ve en la línea
negra del grafico 11.

BOCATOMA DE FONDO
Los datos de entrada necesarios para este diseño son: el caudal máximo diario y a partir de
ese calculamos el Qmin, el Qdiseño, y la resta de los dos anteriores es el valor de Q3.
QMD= 16,9500 m³/s
Q1= Qmin= 0,0357 qmin= 0,0128 m³/s*m
Q2= Qdiseño= 33,9 qdiseño= 4,8429 m³/s*m
Q3= Qmin-Qdiseño= -33,8643 qmin-qdiseño= -4,8301 m³/s*m

Tabla 68. Datos iniciales bocatoma

El ancho de dique que tomamos es: 2,8m

Se adopta un ancho de la rejilla sobre el dique:

L: 7m

CÁLCULO DE LAS ALTURAS DEL AGUA AL INICIO Y FINAL DE LA REJILLA

 En la cresta del vertedero se tiene flujo critico:

Yc: 0,0255 m
 La altura al inicio de la rejilla debe ser ligeramente menor a Yc

Y1: 0,0229

 La velocidad inicial es

V1: 0,5556 m/s

 La energía especifica es

E1: 0,0387
y1/E1: 0,5932 OK

 Para la rejilla se asignan anchos de barras y separación entre ellas:

Separación: b: 1”
Ancho: s: 0,5”

Bocatoma Vertedero Lateral

El diseño de la bocatoma por vertedero lateral se realizará con un ángulo de 90° debido a la
facilidad que éste presenta a la hora de la construcción. Adicional, se recomienda una
limpieza de la rejilla constante para evitar obstrucciones en este. Para el diseño de esta
estructura se trabajará nuevamente con el caudal derivado calculado en el diseño de
bocatoma por orificio. En este diseño se deben suponer algunos datos como:

Ángulo de la rejilla 90

Separación de barras b (ft) 1

Ancho de barras s (ft) 0,500

Velocidad neta (m/s) 0,65

Coeficiente de forma (Cf) 1,79

Perdida de rejilla h (m) 0,015

3*Perdida de rejilla h (m) 0,045

QMD (m³/s) 5,85

Ampliación de Caudal 100
Q diseño (m³/s) 11,70

Carga H (m) 0,05

Sumergencia S 0,096
Tabla 69. Cálculos iniciales

Es necesario obtener un coeficiente de forma el cuál se obtiene de la siguiente tabla:

Ilustracion 5. Coeficiente de forma

Se determina el caudal sumergido a partir de la ecuación de Villemonte,

 Q sumergido (m³/s) 11,835
Le (m) 569,10

Tabla 70. Q sumergido

Procedemos, a estimar la contracción doble, calcular el dato de varillas y espacios y la

longitud neta.

Contracción doble

L (m) 569,11

Varillas y espacios

N 22405,00
L (m) 853,66

Longitud neta

Ln 21,71

Tabla 70.

Realizamos la verificación con Engels de la siguiente manera, donde observamos que si

cumple la longitud y la velocidad neta en la rejilla.

Verificando con Engels

L (m) 1606,48698 OK

Caudal Captado Sumergido 0,44163046

Velocidad neta en la rejilla (m/s) 0,01551997 OK

LB 1,2

H (m) 0,51
 Tabla 71. Verificación con Engels

Hacemos las respectivas iteraciones:

 Perdidas iteración 1

Q (m³/s) 17,69 Velocidad (m/s) 1,60

h (m) 0,35
S 0,32
H1 (m) 0,16

Perdidas iteración 2

Q (m³/s) 16,68 Velocidad (m/s) 1,51

h (m) 0,33
S 0,36
H1 (m) 0,18

Perdidas iteración 3

Q (m³/s) 16,43 Velocidad (m/s) 1,48

h (m) 0,32
S 0,37
H1 (m) 0,19

Perdidas iteración 4

Q (m³/s) 16,36 Velocidad (m/s) 1,48

h (m) 0,32
S 0,37
H1 (m) 0,19

Perdidas iteración 5

Q (m³/s) 16,34 Velocidad (m/s) 1,48

h (m) 0,32
S 0,37
H1 (m) 0,19

Perdidas iteración 6

Q (m³/s) 16,34 Velocidad (m/s) 1,47

h (m) 0,32
S 0,37
H1 (m) 0,19
 Tabla 72. Iteraciones

VERTEDERO DE EXCESOS
Q excesos (m³/s) 4,64
Carga (m) 0,19
C 1,84
L (m) 30,49
Tabla 73. Valores iniciales

TANQUE DE EXCESOS
La trayectoria de la napa exterior define el ancho del tanque de excesos, para esto se usa la
trayectoria del chorro de vertedero del U.S. Bureau of Reclamation, entonces se tiene
V (m/s) 1,00
D (m) 3,86
D (m) 4,00
S (m) 10,00
h (m) 12,25

Sumergencia 2,5D
Altura del vertedero 0,02
Ok
C (m1/2/s) 4,27
C (S.I) (m1/2/s) 2,35

Contracciones 2
Le (m) 30,49
L (m) 30,53
Q (m³/S) 5,94

Tabla 74. Tanque de excesos

TRAYECTORIA NAPA EXTERIOR

V (m/S) 0,802
hv (m/s) 0,033
 A -0,38
B -0,01
C 0,07
m -0,04
D 0,57

X (m) Y (m)
0 0,1216
0,2 0,0399
0,4 -0,2027
0,6 -0,6062
0,8 -1,1706
1 -1,8959
1,2 -2,7821
1,4 -3,8292
1,6 -5,0372
1,8 -6,4061
2 -7,9359

Tabla 75. Trayectoria Napa Exterior

La grafica quedaría:

chorro libre vertedero

1
0
0 0.5 1 1.5 2 2.5
-1
-2
-3
-4
-5
-6
-7
-8
-9

Gráfico 17. Chorro libre vertedero

Compuerta de Limpieza
Caudal de diseño
Caudal medio en el rio
3
Qm= 2,977 m /s
3
QMD= 0,014 m /s

Tabla 76. Caudal medio del rio

 Ampliamos el caudal de diseño en un % de 100

3
Qdiseño= 2,949 m /s
1/2
C= 2,225 m /s
Hd= 0,302 m
y1= 1,802 m

Tabla 77. Ajuste del caudal

 Se le asigna el ancho del canal de limpieza

L = 0,6m
 Se le asigna un valor inicial para la apertura de la compuerta
a = 0,6m

Ilustracion [Link] de gasto

 y1/a= 3,003
Cd= 0,5
a= 1,654 m

Tabla 78. Coeficiente de gasto