ESCUELA COLOMBIANA DE INGENIERÍA JULIO GARAVITO

ACUEDUCTOS Y ALCANTARILLADO

SUBGRUPO A-B

DISEÑO DE UN SISTEMA DE ACUEDUCTO EN EL MUNICIPIO GRAMALOTE, NORTE DE

SANTANDER, COLOMBIA

ENTREGA 3. DISEÑO DE DIQUE VERTEDERO, BOCATOMA Y COMPUERTAS DE LIMPIEZA

PRESENTADO A:
ING. HÉCTOR MATAMOROS

PRESENTADO POR:

ANDRES STEVEN ACHURY CADENA

[Link]@[Link]
JUAN CAMILO MORENO GARCIA
[Link]-ga @[Link]
JORDAN MATEO PEDRAZA CALDERON
[Link]@[Link]

ABRIL 14 DE 2022
TABLA DE CONTENIDO
DETERMINACIÓN HIDROLÓGICA................................................................................4
1. CAUDALES MÁXIMO, MEDIO Y MÍNIMO.....................................................................4
2. CURVA DE DURACIÓN DE CAUDALES.........................................................................11
3. CAUDAL ECOLOGICO.......................................................¡Error! Marcador no definido.
4. TRASLADO DE CAUDALES..........................................................................................15
MODELO HIDRAULICO: HEC-RAS..............................................................................17
1. TOPOBATIMÉTRIA.....................................................................................................17
DISEÑO: DIQUE VERTEDERO.....................................................................................31
1. WES...........................................................................................................................31
2. OGGY........................................................................................................................38
DISEÑO: BOCATOMA................................................................................................44
1. POR ORIFICIO............................................................................................................44
2. POR VERTEDERO LATERAL.........................................................................................47
3. POR REJILLA DE FONDO.............................................................................................53
DISEÑO: COMPUERTAS DE LIMPIEZA........................................................................58
SISTEMA DE BOCATOMA RECOMENDADO......................¡Error! Marcador no definido.
ESQUEMA: BOCATOMA...........................................................................................59
REFERENCIAS BIBLIOGRAFICAS.................................................................................63
DETERMINACIÓN HIDROLÓGICA

1. CAUDALES MÁXIMO, MEDIO Y MÍNIMO

Para la determinacion de los caudales; maximo, medio y minimo, se hace uso de metodos estadisticos
con el fin de verificar la confiabilidad de los datos partiendo de los caudales diarios obtenidos de la
estacion limnigrafica tomados durante 30 años, aquí se hace una tabla resumen con los caudales
maximos, minimos y medios diarios anuales:

CAUDALES
AÑO MAX MINIMO MEDIO
1992 109.148 0.454 28.09376776
1993 143.834 0.686 29.71206027
1994 143.685 0.728 28.80763014
1995 106.25 0.563 25.34480274
1996 106.223 0.532 24.45968852
1997 103.673 0.588 24.5276274
1998 109.616 0.582 24.84592329
1999 104.496 0.557 25.23921096
2000 93.5 0.566 24.07129508
2001 98.355 0.52 24.65498082
2002 102.632 0.494 25.38150411
2003 99.559 0.476 24.26195342
2004 105.264 0.483 24.27650546
2005 106.964 0.503 24.07349315
2006 99.532 0.513 25.30814521
2007 96.968 0.529 24.92812055
2008 91.494 0.599 24.47879235
2009 104.496 0.57 24.23515342
2010 104.312 0.518 24.34996712
2011 100.382 0.578 25.0388
2012 108.732 0.586 24.55072404
2013 110.527 0.524 24.86473699
2014 101.184 0.529 24.74055616
2015 101.102 0 23.80651344
2016 176.517 0.789 39.31581694
2017 162.939 0.879 39.15052329
2018 173.759 0.845 37.49704658
2019 165.829 0.894 38.36051507
2020 144.169 0.448 32.79764754
2021 142.501 0.605 24.65724932
Tabla 1. Caudales diarios anuales 1992-2021.

Ya con estos datos se procede a realizar sus respectivos ajustes con diferentes metodos estadisticos, como lo
son distribucion normal, distribucion log normal y distribucion pearson III.
Como primer paso se aconseja una organización de los datos, tanto para caudales maximos, minimos y
medios, se calcula atraves de la media, desviacion estandar y caudal maximo, obteniendo un periodo de
retorno de 97 años

Para la distribucion normal se tiene en cuenta el promedio de los datos para un z obtenido a traves de la
siguiente ecuacion, para una probabilidad del 99%

x−xmed
z=
desviacion estandar

Para la distribucion normal:

Caudal Caudal
Caudal Medio
Datos Máximo Mínimo
(m3/s)
(m3/s) (m3/s)
Media (m3/s) 117.255 0.571 27.194

Desviación
25.577 0.162 4.946
Estandar (m3/s)

Q Máx (m3/s) 176.517 0.894 39.316

n (años) 25 25 25
Z 2.317 2.317 2.317
P (X < Xo) 99% 99% 99%
Cs 1.285 0.818 1.763
Tabla 2: Distribucion normal

Para la distribucion log normal, se hace uso de logaritmo natural a los datos de caudales diarios anuales;
maximos, minimos y medios:
 Caudal Caudal Log Caudal Log Caudal LogCaudal
Caudal Max
No. Año Minimo Medio Máximo Mínimo Medio
(m3/s)
(m3/s) (m3/s) (m3/s) (m3/s) (m3/s)
1 1992 109.148 0.454 28.0937678 4.693 -0.790 3.3355
2 1993 143.834 0.686 29.7120603 4.969 -0.377 3.3916
3 1994 143.685 0.728 28.8076301 4.968 -0.317 3.3606
4 1995 106.25 0.563 25.3448027 4.666 -0.574 3.2326
5 1996 106.223 0.532 24.4596885 4.666 -0.631 3.1970
6 1997 103.673 0.588 24.5276274 4.641 -0.531 3.1998
7 1998 109.616 0.582 24.8459233 4.697 -0.541 3.2127
8 1999 104.496 0.557 25.239211 4.649 -0.585 3.2284
9 2000 93.5 0.566 24.0712951 4.538 -0.569 3.1810
10 2001 98.355 0.52 24.6549808 4.589 -0.654 3.2050
11 2002 102.632 0.494 25.3815041 4.631 -0.705 3.2340
12 2003 99.559 0.476 24.2619534 4.601 -0.742 3.1889
13 2004 105.264 0.483 24.2765055 4.656 -0.728 3.1895
14 2005 106.964 0.503 24.0734932 4.672 -0.687 3.1811
15 2006 99.532 0.513 25.3081452 4.600 -0.667 3.2311
16 2007 96.968 0.529 24.9281205 4.574 -0.637 3.2160
17 2008 91.494 0.599 24.4787923 4.516 -0.512 3.1978
18 2009 104.496 0.57 24.2351534 4.649 -0.562 3.1878
19 2010 104.312 0.518 24.3499671 4.647 -0.658 3.1925
20 2011 100.382 0.578 25.0388 4.609 -0.548 3.2204
21 2012 108.732 0.586 24.550724 4.689 -0.534 3.2007
22 2013 110.527 0.524 24.864737 4.705 -0.646 3.2135
23 2014 101.184 0.529 24.7405562 4.617 -0.637 3.2084
24 2015 101.102 0 23.8065134 4.616 - 3.1700
25 2016 176.517 0.789 39.3158169 5.173 -0.237 3.6716
26 2017 162.939 0.879 39.1505233 5.093 -0.129 3.6674
27 2018 173.759 0.845 37.4970466 5.158 -0.168 3.6243
28 2019 165.829 0.894 38.3605151 5.111 -0.112 3.6470
29 2020 144.169 0.448 32.7976475 4.971 -0.803 3.4904
30 2021 142.501 0.605 24.6572493 4.959 -0.503 3.2051
Tabla 3:logaritmo de los valores de caudal

Para esta distribucion log normal se realiza el mismo procedimiento anterior, solo que para determinar los
caudales se hace atraves de:

xpromedio+(z × desviacion estandar)

Q=e

Obteniendo los siguientes caudales:

 Caudal Caudal
Caudal Medio
Datos Máximo Mínimo
(m3/s)
(m3/s) (m3/s)
Q Medio
4.744 -0.544 3.289
(m3/s)

Desviación
Estandar 0.198 0.189 0.161
(m3/s)

Q Máx (m3/s) 5.173 -0.803 3.672

n (años) 25 25 25
Z 2.17 -1.37 2.37
P (X < Xo) 99% 99% 99%
Cs 1.128 1.093 1.654
Q 176.517 0.448 39.31581694
Tabla 4:Distribucion log normal

Una vez obtenido los caudales por medio de la distribucion log normal, se hace uso del metodo de
distribucion pearson III:
Para este metodo se realiza a partir de un orden, en este caso se organizan caudales de mayor a menor, con
un respectivo conteo de 1 al ultimo numero de datos, que en este caso seria de 30:

Caudal Caudal
Caudal Max
M Minimo Medio
(m3/s)
(m3/s) (m3/s)
1 176.517 0.789 39.3158169
2 173.759 0.845 37.4970466
3 165.829 0.894 38.3605151
4 162.939 0.879 39.1505233
5 144.169 0.448 32.7976475
6 143.834 0.686 29.7120603
7 143.685 0.728 28.8076301
8 142.501 0.605 24.6572493
9 110.527 0.524 24.864737
10 109.616 0.582 24.8459233
11 109.148 0.454 28.0937678
12 108.732 0.586 24.550724
13 106.964 0.503 24.0734932
14 106.25 0.563 25.3448027
15 106.223 0.532 24.4596885
16 105.264 0.483 24.2765055
17 104.496 0.557 25.239211
18 104.496 0.57 24.2351534
19 104.312 0.518 24.3499671
20 103.673 0.588 24.5276274
21 102.632 0.494 25.3815041
22 101.184 0.529 24.7405562
23 101.102 0 23.8065134
24 100.382 0.578 25.0388
25 99.559 0.476 24.2619534
26 99.532 0.513 25.3081452
27 98.355 0.52 24.6549808
28 96.968 0.529 24.9281205
29 93.5 0.566 24.0712951
30 91.494 0.599 24.4787923

Ya organizado los datos se procede a determinar unos coeficientes k para cada uno de los caudales,
determinados de la siguiente manera:
 Qi
K=
Qpromedio

Tambien se calcula la probabilidad de manera porcentual para cada uno de los datos, donde N es la cantidad
total de datos:

Mi
%Prob= ×100
N +1

Se procede al calculo del coeficiente de variacion y el coeficiente asimetria para caudales maximos, minimos
y medios

Coef De Variacion=
√ ∑ (k−1)2
N −1

Coef De asimetria=
∑ (k −1)3
3
(N −1) × coeficiente variacion

Datos calculados:

Caudal Caudal
Caudal Max
M Minimo Medio k1 k2 k3 k1-1 k2-1 k3-1
(m3/s) 3
(m /s) (m3/s)
1 176.517 0.789 39.3158169 1.50541471 1.38114132 1.44573431 0.50541471 0.38114132 0.44573431
2 173.759 0.845 37.4970466 1.48189327 1.4791691 1.37885388 0.48189327 0.4791691 0.378853881
3 165.829 0.894 38.3605151 1.41426274 1.5649434 1.41060563 0.41426274 0.5649434 0.410605632
4 162.939 0.879 39.1505233 1.38961554 1.53868596 1.43965608 0.38961554 0.53868596 0.439656077
5 144.169 0.448 32.7976475 1.22953672 0.7842222 1.20604601 0.22953672 -0.2157778 0.206046015
6 143.834 0.686 29.7120603 1.22667969 1.20084024 1.09258177 0.22667969 0.20084024 0.092581773
7 143.685 0.728 28.8076301 1.22540895 1.27436107 1.05932377 0.22540895 0.27436107 0.059323766
8 142.501 0.605 24.6572493 1.21531128 1.05905006 0.90670458 0.21531128 0.05905006 -0.093295419
9 110.527 0.524 24.864737 0.94262293 0.91725989 0.91433439 -0.05737707 -0.08274011 -0.085665613
10 109.616 0.582 24.8459233 0.93485352 1.01878866 0.91364256 -0.06514648 0.01878866 -0.086357437
11 109.148 0.454 28.0937678 0.93086221 0.79472517 1.03307338 -0.06913779 -0.20527483 0.033073381
12 108.732 0.586 24.550724 0.92731438 1.02579064 0.9027874 -0.07268562 0.02579064 -0.097212602
13 106.964 0.503 24.0734932 0.91223609 0.88049947 0.88523851 -0.08776391 -0.11950053 -0.114761495
14 106.25 0.563 25.3448027 0.90614679 0.98552923 0.93198752 -0.09385321 -0.01447077 -0.068012475
15 106.223 0.532 24.4596885 0.90591652 0.93126386 0.89943981 -0.09408348 -0.06873614 -0.10056019
16 105.264 0.483 24.2765055 0.89773775 0.84548956 0.89270374 -0.10226225 -0.15451044 -0.107296259
17 104.496 0.557 25.239211 0.8911879 0.97502626 0.92810467 -0.1088121 -0.02497374 -0.071895332
18 104.496 0.57 24.2351534 0.8911879 0.99778271 0.89118313 -0.1088121 -0.00221729 -0.10881687
19 104.312 0.518 24.3499671 0.88961867 0.90675691 0.8954051 -0.11038133 -0.09324309 -0.104594902
20 103.673 0.588 24.5276274 0.884169 1.02929163 0.90193808 -0.115831 0.02929163 -0.098061919
21 102.632 0.494 25.3815041 0.87529089 0.86474501 0.93333712 -0.12470911 -0.13525499 -0.06666288
22 101.184 0.529 24.7405562 0.86294171 0.92601237 0.90976797 -0.13705829 -0.07398763 -0.090232032
23 101.102 0 23.8065134 0.86224238 0 0.87542104 -0.13775762 -1 -0.124578962
24 100.382 0.578 25.0388 0.8561019 1.01178667 0.92073509 -0.1438981 0.01178667 -0.07926491
25 99.559 0.476 24.2619534 0.84908299 0.83323608 0.89216863 -0.15091701 -0.16676392 -0.107831371
26 99.532 0.513 25.3081452 0.84885273 0.89800443 0.93063954 -0.15114727 -0.10199557 -0.069360459
27 98.355 0.52 24.6549808 0.83881475 0.91025791 0.90662116 -0.16118525 -0.08974209 -0.093378837
28 96.968 0.529 24.9281205 0.8269858 0.92601237 0.91666515 -0.1730142 -0.07398763 -0.083334852
29 93.5 0.566 24.0712951 0.79740917 0.99078072 0.88515768 -0.20259083 -0.00921928 -0.114842323
30 91.494 0.599 24.4787923 0.78030112 1.04854709 0.9001423 -0.21969888 0.04854709 -0.099857698

Tabla 5: datos para determinar los coeficientes

 Caudal Caudal
Caudal Max
M Minimo Medio k1 k2 k3 k1-1 k2-1 k3-1 (k1-1)^2 (k2-1)^2 (k3-1)^2 (k1-1)^3 (k2-1)^3 (k3-1)^3
(m3/s)
(m3/s) (m3/s)
%probab
1 176.517 0.789 39.3158 1.5054 1.3811 1.4457 0.5054 0.3811 0.4457 0.2554 0.1453 0.1987 0.1291 0.0554 0.0886 3.2
2 173.759 0.845 37.4970 1.4819 1.4792 1.3789 0.4819 0.4792 0.3789 0.2322 0.2296 0.1435 0.1119 0.1100 0.0544 6.5
3 165.829 0.894 38.3605 1.4143 1.5649 1.4106 0.4143 0.5649 0.4106 0.1716 0.3192 0.1686 0.0711 0.1803 0.0692 9.7
4 162.939 0.879 39.1505 1.3896 1.5387 1.4397 0.3896 0.5387 0.4397 0.1518 0.2902 0.1933 0.0591 0.1563 0.0850 12.9
5 144.169 0.448 32.7976 1.2295 0.7842 1.2060 0.2295 -0.2158 0.2060 0.0527 0.0466 0.0425 0.0121 -0.0100 0.0087 16.1
6 143.834 0.686 29.7121 1.2267 1.2008 1.0926 0.2267 0.2008 0.0926 0.0514 0.0403 0.0086 0.0116 0.0081 0.0008 19.4
7 143.685 0.728 28.8076 1.2254 1.2744 1.0593 0.2254 0.2744 0.0593 0.0508 0.0753 0.0035 0.0115 0.0207 0.0002 22.6
8 142.501 0.605 24.6572 1.2153 1.0591 0.9067 0.2153 0.0591 -0.0933 0.0464 0.0035 0.0087 0.0100 0.0002 -0.0008 25.8
9 110.527 0.524 24.8647 0.9426 0.9173 0.9143 -0.0574 -0.0827 -0.0857 0.0033 0.0068 0.0073 -0.0002 -0.0006 -0.0006 29.0
10 109.616 0.582 24.8459 0.9349 1.0188 0.9136 -0.0651 0.0188 -0.0864 0.0042 0.0004 0.0075 -0.0003 0.0000 -0.0006 32.3
11 109.148 0.454 28.0938 0.9309 0.7947 1.0331 -0.0691 -0.2053 0.0331 0.0048 0.0421 0.0011 -0.0003 -0.0086 0.0000 35.5
12 108.732 0.586 24.5507 0.9273 1.0258 0.9028 -0.0727 0.0258 -0.0972 0.0053 0.0007 0.0095 -0.0004 0.0000 -0.0009 38.7
13 106.964 0.503 24.0735 0.9122 0.8805 0.8852 -0.0878 -0.1195 -0.1148 0.0077 0.0143 0.0132 -0.0007 -0.0017 -0.0015 41.9
14 106.25 0.563 25.3448 0.9061 0.9855 0.9320 -0.0939 -0.0145 -0.0680 0.0088 0.0002 0.0046 -0.0008 0.0000 -0.0003 45.2
15 106.223 0.532 24.4597 0.9059 0.9313 0.8994 -0.0941 -0.0687 -0.1006 0.0089 0.0047 0.0101 -0.0008 -0.0003 -0.0010 48.4
16 105.264 0.483 24.2765 0.8977 0.8455 0.8927 -0.1023 -0.1545 -0.1073 0.0105 0.0239 0.0115 -0.0011 -0.0037 -0.0012 51.6
17 104.496 0.557 25.2392 0.8912 0.9750 0.9281 -0.1088 -0.0250 -0.0719 0.0118 0.0006 0.0052 -0.0013 0.0000 -0.0004 54.8
18 104.496 0.57 24.2352 0.8912 0.9978 0.8912 -0.1088 -0.0022 -0.1088 0.0118 0.0000 0.0118 -0.0013 0.0000 -0.0013 58.1
19 104.312 0.518 24.3500 0.8896 0.9068 0.8954 -0.1104 -0.0932 -0.1046 0.0122 0.0087 0.0109 -0.0013 -0.0008 -0.0011 61.3
20 103.673 0.588 24.5276 0.8842 1.0293 0.9019 -0.1158 0.0293 -0.0981 0.0134 0.0009 0.0096 -0.0016 0.0000 -0.0009 64.5
21 102.632 0.494 25.3815 0.8753 0.8647 0.9333 -0.1247 -0.1353 -0.0667 0.0156 0.0183 0.0044 -0.0019 -0.0025 -0.0003 67.7
22 101.184 0.529 24.7406 0.8629 0.9260 0.9098 -0.1371 -0.0740 -0.0902 0.0188 0.0055 0.0081 -0.0026 -0.0004 -0.0007 71.0
23 101.102 0 23.8065 0.8622 0.0000 0.8754 -0.1378 -1.0000 -0.1246 0.0190 1.0000 0.0155 -0.0026 -1.0000 -0.0019 74.2
24 100.382 0.578 25.0388 0.8561 1.0118 0.9207 -0.1439 0.0118 -0.0793 0.0207 0.0001 0.0063 -0.0030 0.0000 -0.0005 77.4
25 99.559 0.476 24.2620 0.8491 0.8332 0.8922 -0.1509 -0.1668 -0.1078 0.0228 0.0278 0.0116 -0.0034 -0.0046 -0.0013 80.6
26 99.532 0.513 25.3081 0.8489 0.8980 0.9306 -0.1511 -0.1020 -0.0694 0.0228 0.0104 0.0048 -0.0035 -0.0011 -0.0003 83.9
27 98.355 0.52 24.6550 0.8388 0.9103 0.9066 -0.1612 -0.0897 -0.0934 0.0260 0.0081 0.0087 -0.0042 -0.0007 -0.0008 87.1
28 96.968 0.529 24.9281 0.8270 0.9260 0.9167 -0.1730 -0.0740 -0.0833 0.0299 0.0055 0.0069 -0.0052 -0.0004 -0.0006 90.3
29 93.5 0.566 24.0713 0.7974 0.9908 0.8852 -0.2026 -0.0092 -0.1148 0.0410 0.0001 0.0132 -0.0083 0.0000 -0.0015 93.5
30 91.494 0.599 24.4788 0.7803 1.0485 0.9001 -0.2197 0.0485 -0.0999 0.0483 0.0024 0.0100 -0.0106 0.0001 -0.0010 96.8
prom 117.254733 0.57126667 27.1943584 total 1.37988643 2.331233687 0.95933254 0.36107876 -0.504382709 0.287150345 50
PROM
coef variac 0.21813365 0.283526663 0.18188027
coef de asime 1.19959852 -0.76309828 1.6457148

Con estos datos se procede a determinar el coeficiente de Foster (Tp) el cual se puede obtener a partir
de tablas en funcion del coeficiente de asimetria, teniendo en cuenta el % de probabilidad calculado
anteriormente pero promediandolo para ajustar los tres caudales; max, minimo y medio:
 max min medio

TP -0.19 -0.13 -0.26

Q 112.395065 0.550210654 25.9083679

Para obtener los caudales medios para gumbel se procede a realizar la prueba de bondad

Caudal Medio Diario Anual f(Q˃Qi) =

Año ordenado i Qi Zi
(m³/s) (i/(n+1)) Kimbal

1992 28.093768 39.31581694 1 0.032258065 18.0509844 -1.8485963

1993 29.712060 39.15052329 2 0.064516129 19.6865028 -1.5179292
1994 28.807630 38.36051507 3 0.096774194 20.7636471 -1.3001534
1995 25.344803 37.49704658 4 0.129032258 21.6004105 -1.1309776
1996 24.459689 32.79764754 5 0.161290323 22.3018144 -0.9891686
1997 24.527627 29.71206027 6 0.193548387 22.9164895 -0.8648944
1998 24.845923 28.80763014 7 0.225806452 23.4712735 -0.7527288
1999 25.239211 28.09376776 8 0.258064516 23.9827262 -0.6493239
2000 24.071295 25.38150411 9 0.290322581 24.4619126 -0.5524426
2001 24.654981 25.34480274 10 0.322580645 24.9166984 -0.4604945
2002 25.381504 25.30814521 11 0.35483871 25.3529715 -0.3722894
2003 24.261953 25.23921096 12 0.387096774 25.7753474 -0.2868939
2004 24.276505 25.0388 13 0.419354839 26.1876047 -0.2035442
2005 24.073493 24.92812055 14 0.451612903 26.5929729 -0.1215874
2006 25.308145 24.86473699 15 0.483870968 26.9943349 -0.0404405
2007 24.928121 24.84592329 16 0.516129032 27.3943819 0.04044051
2008 24.478792 24.74055616 17 0.548387097 27.7957438 0.12158738
2009 24.235153 24.65724932 18 0.580645161 28.201112 0.20354423
2010 24.349967 24.65498082 19 0.612903226 28.6133693 0.28689392
2011 25.038800 24.55072404 20 0.64516129 29.0357452 0.37228936
2012 24.550724 24.5276274 21 0.677419355 29.4720184 0.46049454
2013 24.864737 24.47879235 22 0.709677419 29.9268042 0.55244258
2014 24.740556 24.45968852 23 0.741935484 30.4059906 0.64932391
2015 23.806513 24.34996712 24 0.774193548 30.9174432 0.75272879
2016 39.315817 24.27650546 25 0.806451613 31.4722273 0.86489436
2017 39.150523 24.26195342 26 0.838709677 32.0869024 0.98916863
2018 37.497047 24.23515342 27 0.870967742 32.7883062 1.13097761
2019 38.360515 24.07349315 28 0.903225806 33.6250697 1.30015343
2020 32.797648 24.07129508 29 0.935483871 34.702214 1.51792916
2021 24.657249 23.80651344 30 0.967741935 36.3377323 1.84859629

Intervalo Xi Fi Xi * Fi (xi-X)2 fi*(xi-X)2

23.80651344 26.39139736 25.0989554 22 552.1770188 5.41211614 119.066555
26.39139736 28.97628127 27.68383932 2 55.36767863 0.06681625 0.1336325
28.97628127 31.56116519 30.26872323 1 30.26872323 8.08476608 8.08476608
31.56116519 34.14604911 32.85360715 1 32.85360715 29.4659656 29.4659656
34.14604911 36.73093302 35.43849107 0 0 64.2104149 0
36.73093302 39.31581694 38.02337498 4 152.0934999 112.318114 449.272456
30 606.023375

X 27.42535092
s 4.57136279
 Valor
Intervalo ni fi Fi acumulada Zi F(xi) p(xi) X^2
esperado
23.80651344 26.39139736 22 0.73333333 0.733333333 -0.2261806 0.41053048 0.41053048 9.03167057 7.624547456
26.39139736 28.97628127 2 0.06666667 0.8 0.3392709 0.63279717 0.22226669 0.44453337 0.32119804
28.97628127 31.56116519 1 0.03333333 0.833333333 0.90472239 0.81719376 0.59492707 0.59492707 0.530128051
31.56116519 34.14604911 1 0.03333333 0.866666667 1.47017388 0.92924267 0.33431559 0.33431559 0.270972465
34.14604911 36.73093302 0 0 0.866666667 2.03562538 0.97910601 0.64479041 0 #¡DIV/0!
36.73093302 39.31581694 4 0.13333333 1 2.60107687 0.99535342 0.350563 1.40225202 1.148263371
30 2.55739325

Tr p(x<=xi) Qi
96 0.989583333 37.98973074
25 0.96 35.42837209
40 0.975 36.38505735
50 0.98 36.81378227
100 0.99 38.05993103

Tambien se procedio a realizar la distribucion para el caudal maximo y minimo por gumbel, obteniendo los
siguientes caudales para cada metodo de distribucion; normal, log normal, pearson tipo III y gumbel:

Caudal Caudal Caudal

Distribución Máximo Mínimo Medio
(m3/s) (m3/s) (m3/s)
Normal 176.517 0.894 39.316

Log - Normal 176.517 0.448 39.3158169

Pearson tipo III 112.395065 0.55021065 25.9083679

Gumbel 105.3542 31.8050593 37.9897307

Se pudo observar que los datos que mejor se acoplan a la situacion hidrologia son los de la distribucion
Pearson III debido a la asociacion lineal y a la distribucion de las variables de forma aproximada

Por lo tanto los datos de caudales a selección son los siguientes:

Caudal m3/s L/s

Maximo 112.40 112395.06
Minimo 0.55 550.21
Medio 25.91 25908.37

2. CURVA DE DURACIÓN DE CAUDALES

Para la curva de duracion de caudales se hace teniendo en cuenta todos los caudales obtenidos de la
estacion limnigrafica, donde se organizan del mayor al menor caudal y una vez organizado se rocede
acalcular la probabilidad de excedencia para cada uno de los datos, con esto se procede a realizar la
curva de duracion de caudales:
 Grafica de duracion de caudales
200
180
160
140
120
Q (m3/s)

100
80
60
40
20
0
0.0% 10.0% 20.0% 30.0% 40.0% 50.0% 60.0% 70.0% 80.0% 90.0% 100.0%
P >=

Gráfica 1. Curva duracion de caudales. Fuente: propia

una ves graficada la curva de duracion de caudales, procedemos a calcular el indice de regulacion y
retencion hidrica como:

IRH=v p ∕ vT

Ecuación 1. Población proyectada a partir del método de crecimiento lineal

V p :volumen bajo larecta definida por el q 50 en lacurva de duracion

V t :es el volumen bajola curva de duracionde caudales

Ya con esto definido se procede a calcular cada uno de los datos a traves de areas de trapecios para
calcular las areas bajo la curva de duracion de caudales para determinar su respectivo dato Vp y Vt.

T util (Años) 25
Riesgo asumido 0.23
Periodo de retorno (años) 96

Q50 (m3/s) 14.78

Q75 (m3/s) 3.405
Q95 (m3/s) 0.972
Q95 =CAUDAL ECOLOGICO

VT 27.201901
VP 9.97698439
IRH 0.36677526 MUY BAJA

Estos valores permiten obtner el indice de regulacion, por lo tanto se compara el valor obtenido con
tablas para ver su descripcion:
 Tabla. IRH descripcion a partir de rangos de valores

Como se puede observar el valor obtenido de indice de regulacion y retencion hidrica es de 0.36 por lo tanto
debido a que este es menor a 0.50 se le considera como descripcion muy baja.

Como se puede observar tambien se procede a asignar el valor del caudal ecologico un valor de 0.97 M3/S,
debido a que el caudal ecologico corresponde al valor de Q95
3. TRASLADO DE CAUDALES

Con base en la cuenca donde está ubicada la estación limnimetrica de la cual fue extraída la
información histórica de caudales, se procede a ubicar el punto de captación en uno de los cuerpos de
agua que alimentan a esta, a una cota mínima de 625m, donde a partir de este punto si designara la
subcuenca que se empleara para trasladar los caudales mínimo, máximo y medio que serán
empleados dentro del modelo de flujo HEC-RAS y el diseño de las obras de captación. A
continuación, se presenta el plano en planta de la cuenca hidrográfica:

Ilustración 1. Cuenta hidrográfica

El punto verde que se visualiza en la ilustración corresponde a la estación limnimetrica. A partir de

esta se procede a ubicar el punto donde estará la obra de captación y a su vez se delimitará la
subcuenca. Al haber determinado la cuenca hidrográfica y la subcuenca, con la ayuda del software
Civil 3D se obtienen los valores correspondientes a las áreas de estas, con las cuales se podrán
calcular los caudales trasladados. La subcuenca y el punto de captación se muestran a continuación:
 Subcuenca

Bocatoma

Cuenca

Ilustración 2. Localización de bocatoma y subcuenca definida

La figura delineada de color amarillo corresponde a la subcuenca definida, donde dentro de ella
estará ubicada a captación del proyecto, a una cota de 825m, la cual cumple con el limite definido.

Se empleó el método de áreas para el traslado de caudales, para este se apoya de la siguiente
ecuación:

Donde:

QB :Caudal trasladado a punto de captación

Q A :Caudal estimado hidrologico
Área cuenca : áreatotal de la cuenca incluyendo la subcuenca
Área subcuenca :área delimitada de subcuenca donde se ubica la captación

Los datos correspondientes al área se obtienen mediante el uso de la herramienta CAD la cual calcula
el área de la cuenca y subcuenca, mientras que para el Qa se determina del análisis hidrológico
presentado anteriormente. Los datos fueron los siguientes:
 AREA CUENCA TOTAL M2 6274060,83
AREA SUB CUENCA M2 434603,50
Tabla 6. Área de cuencas

Al aplicar el método de áreas se obtienen los siguientes caudales trasladados:

Correccion de caudales (m3/s)

Caudal min Caudal med Caudal max
0.038 1.795 7.786
Tabla 7. Caudales trasladados

MODELO HIDRAULICO: HEC-RAS

1. TOPOBATIMÉTRIA
A partir del levantamiento topo batimétrico del trazado del eje de rio con sus respectivas curvas de nivel, se
procede a ubicar una zona donde estará ubicada la obra de captación, esta deberá cumplir con los siguientes
aspectos: secciones relativamente profundas, eje del rio plano, procurar que no quede en pendientes
pronunciadas, evitar cambios bruscos en la dirección (curvas) y relativamente ancho.

A continuación se señala en el plano la ubicación seleccionada:

Con
base en
la

Bocatoma

Ilustración 3. Ubicación bocatoma

ilustración anterior, se realiza un zoom de la zona seleccionada para la ubicación de bocatoma, la cual
corresponde a una zona por la cual el trazado va entre curvas de nivel, lo que hace referencia a un trayecto
plano y no presenta variaciones bruscas en su dirección aguas abajo.
 Ilustración 4. Zoom Ubicación bocatoma
Mediante el uso
del software Civil 3D se definieron las secciones repartidas sobre el tramo de longitud 500m, a cada 30m,
con un ancho de 50m, asegurando que ninguna de las secciones vistas en planta se interseque lo cual
afectaría el modelo de flujo. En total fueron 16 secciones obtenidas, siendo la encerrada en rojo la de interés,
ya que cumple con las mínimas recomendaciones expresadas anteriormente.

El perfil obtenido para dicha sección generado por Civil 3D se presenta a continuación: Nota: las unidades
se encuentran en pies debido al origen del archivo CAD, pero a realizar su exportación a HEC-RAS se hizo
en metros.

Ilustración 5. Sección de interés. Unidades: ft

A partir de las secciones generadas en Civil 3D se procedió a exportar la data de secciones a un modelo de
HEC-RAS encargado de reproducir el comportamiento del rio que como producto final se obtendrá el ancho
superficial que será necesario para el diseño del dique y demás obras hidráulicas.

Los datos de Right Over Back (ROB), Left Over Back (LOB) y Channel serán parámetros que el modelo de
HEC identificara automáticamente de los datos que provienen del Civil. Por otro lado, el modelo requerirá
ser alimentado por los caudales máximo, medio y mínimo, los cuales se obtuvieron en el traslado de
caudales realizado anteriormente. El modelo también requiere de un coeficiente de Manning de entrada, el
cual se vera modificado luego de correr el programa. Para este caso se utilizo el valor de 0.03
correspondiente a canales dragados.

Una vez corrido el modelo se obtuvieron los siguientes resultados:

Secciones:

Se obtuvieron 16 secciones, correspondientes a las exportadas desde Civil 3D.

Ilustración 6. Tramos de secciones modelado en HEC-RAS

Las secciones llevan como nombre: 3397.95, 3420, 3450, 3480, 3510, 3540, 3570, 3600, 3630, 3660,
3690, 3720, 3750, 3780, 3810, 3840. A continuación, se muestran las secciones obtenidas en el
software:

Ilustración 7. Estación 3420

Ilustración 8. Estación 3450

Ilustración 9. Estación 3480

Ilustración 10. Estación 3510

Ilustración 11. Estación 3540

Ilustración 12. Estación 3570

Ilustración 13. Estación 3600

Ilustración 14. Estación 3630

Ilustración 15. Estación 3660

Ilustración [Link]ón 3690

Ilustración [Link]ón 3720

Ilustración 18. Estación 3750

Ilustración 19. Estación 3780

Ilustración 20. Estación 3810

Ilustración 21. Estación 3840

Como se indicó anteriormente, en la ilustración 5 se designó la sección de interés, que para el caso
del modelamiento en HEC-RAS corresponde a la estación 3450. Las leyendas de las secciones
corresponden a lo siguiente:

PF1: Caudal Máximo

PF2: Caudal Medio
PF3: Caudal Mínimo

Ilustración 22. Estación 3450. Ubicación de bocatoma

Ancho superficial

Con base en las estaciones modeladas, se extrajeron las tablas de salida del modelo, en la cual se
visualizan datos correspondientes a cada sección a partir de los caudales ingresados, y como
resultado nos da el ancho y área superficial. A continuación, se presenta la tabla de salidas:
 Ilustración 23. Tabla de salidas. Sección 3450

El ancho superficial de acuerdo con la tabla corresponde a un valor de 14,45m para el caudal
máximo, 12,46m para el caudal medio y 10,60 para el caudal mínimo.

El criterio de selección de la estación 3450 corresponde a las mínimas exigencias de anchura (medida
moderada), sin cambios bruscos, en una zona preferiblemente plana y con una profundidad
considerable.

Perfil de flujo

A partir del modelamiento, se obtuvo el perfil de flujo que se presenta a continuación:

 Ilustración 24. Perfil de flujo

Ilustración 25. Perfil 2

Omitiendo el primer punto del perfil, el cual por problemas en el CAD corresponden a una cota de valor
cero, de ese punto en adelante se evidencia un flujo subcrítico, semipermanente. Bajo esas condiciones se
corre el modelo, definiendo dos zonas de control, una aguas abajo y otra agua arriba. Estas condiciones se
explican debido al punto donde fue ubicada la bocatoma, siendo una zona plana sin cambios bruscos.

Modelo Río

La siguiente ilustración corresponde al comportamiento del río con base en los caudales trasladados,
lugar donde estará ubicada la obra de captación.
 Ilustración 26. Modelo final del río

Con base en los resultados obtenidos y los criterios de selección optados el modelo de rio fue exitoso
y cumple con lo necesario para que en este tramo sea ubicada la obra de captación.

DISEÑO: DIQUE VERTEDERO

1. WES

En el caso del dique tipo WES, su geométrica se modela a través de la siguientes ecuación
n n−1
X =K∗H ∗Y

Donde:
  X y Y: son las coordenadas del perfil de la cresta con origen en el punto más alto del vertedero
 H: es la altura de carga sobre el vertedero
 K, n: son parámetros que dependen de la pendiente de la cara de aguas arriba de la
estructura

Para este diseño del dique vertedero WES se tiene un radio de curvatura R2 constante un radio de
curvatura R1 constante hasta la cresta del vertedero.

Adicionalmente para completar el diseño de la bocatoma es necesario para no afectar las

condiciones del río, disipar la energía de exceso almacenadas mediante un resalto hidráulico en un
cuenco disipador de energía controlado por una estructura de control de aguas abajo, como lo es un
vertedero de pared gruesa. Para esto en primer lugar se realiza una curva de radio constante, que
entregue al cuenco, este a la vez posea un vertedero de pared gruesa que genere un conflicto de
controles entre flujo en supercrítico aguas arriba a un flujo subcrítico aguas abajo que va a resolver
mediante un resalto hidráulico.
ESCUELA COLOMBIANA DE INGENIERÍA JULIO GARAVITO PROGRAMA DE
INGENIERÍA CIVIL
PROFESOR: Ing. Héctor Matamoros
ASIGNATURA: ACUEDUCTOS Y ALCANTARILLADO

Para el diseño del dique vertedero tipo WES se parte del estudio hidrológico
realizado anteriormente donde se determinaron los valores (Qmax, Qm y Qmin)
en el punto de la bocatoma, y también de la sección del río luego del
levantamiento en Hec-RAS seleccionada como apta para establecer el punto
de la bocatoma, donde se obtienen las láminas mínima, láminas medias y
láminas máximas de agua del río, las cuales se van transformar físicamente
mediante canalizaciones y se va a tener una estructura de toma. Se busca que
la sección transformada reproduzca en el río una condicione de flujo
semejante a las de condición natural.

Q diseño: 7.79 m3/s

L(ancho presa) 15 m
h adoptada: 1.5 m

Se define el ancho de la presa que será una constante en el diseño y se plantea

inicialmente la altura de esta como h0 esta se puede variar para lograr obtener
estabilidad en el diseño y lograr volver semejante nuestro sistema al comportamiento
natural. Inicialmente se supone que h/hd>1.33, razón por la cual la altura de
velocidad aguas arriba del dique va a ser despreciable lo que implica un valor de Cd de
2.225 en el sistema internacional.

Se determino la carga sobre el vertedero a partir de.

2
Qdiseñ o 3
He=
CL

Luego se calculó la columna de velocidad.

2
3
Qdiseñ o
L
Hv=
ho+ He
2∗9.806

Con dicho valor se recalcula He/Hd y se busca el valor en la grafica 14.4 del libro de
Ven Te Chow, para así corregir C.
REVISION
h/Hd 3.958182755 cumple Hv 0.00389124 m
He/Hd 1.010268171 C/Cd 1
C 2.225 Hv 0.00387518 m
Hd 0.378961784 m He 0.38285303 m

Qmin 0.038 m^3/s

Hdmin 0.010903725 m

Qmedio 1.884 m^3/s

Hmedio 0.14715434 m
33
ESCUELA COLOMBIANA DE INGENIERÍA JULIO GARAVITO PROGRAMA DE
INGENIERÍA CIVIL
PROFESOR: Ing. Héctor Matamoros
ASIGNATURA: ACUEDUCTOS Y ALCANTARILLADO

Apartir del Hd Max se calculan los factores de K, n, r1, r2, M, y N para las condiciones
ya establecidas de ancho, y alto del Dique y caudal de diseño cuyos resultados se
muestran en la Tabla 19 Resultado valores de diseño vertedero [Link] 19
Resultado valores de diseño vertedero W.

K n r1 r2 M N
Vertical 1 2 1.85 0.189 0.076 0.066 0.107
3a1 2 1.936 1.836 0.258 0.080 0.053 0.090
3a2 3 1.939 1.81 0.182 0.083 0.044 0.081
3a3 4 1.873 1.776 0.171 0.000 0.045 0.000

Ya que se escogió un vertedero vertical, se tendrán en cuenta los factores de dicha

condición,

K 2
n 1.85
r1 0.189 m
r2 0.076 m
M 0.066 m
N 0.107 m

34
ESCUELA COLOMBIANA DE INGENIERÍA JULIO GARAVITO PROGRAMA DE
INGENIERÍA CIVIL
PROFESOR: Ing. Héctor Matamoros
ASIGNATURA: ACUEDUCTOS Y ALCANTARILLADO

Finalmente, para el cálculo del vertedero WES se procede a aplicar la ecuación de su

funcionamiento
𝑋𝑛 = (𝐾𝐻𝑛−1𝑌) con los valores de K, n, y H calculados. Cuyo resultado se
muestra en la Tabla 20 Vertedero WES. Por efectos prácticos se da el valor
opuesto de Y en la tabla para poder graficar adecuadamente el vertedero
como se muestra en Ilustración 10 Vertedero WES.

X(m) Y(m)
0 0
0.1 -0.0161125
0.2 -0.0580857
0.3 -0.1229811
0.4 -0.2093991
0.5 -0.316416
0.6 -0.4433469
0.7 -0.5896513
0.8 -0.7548843
0.9 -0.9386691
1 -1.1406802
1.1 -1.360631
1.2 -1.5982665
1.3 -1.8533571
1.4 -2.1256941
1.5 -2.4150867
1.6 -2.7213592
1.7 -3.0443488
1.8 -3.383904

35
ESCUELA COLOMBIANA DE INGENIERÍA JULIO GARAVITO PROGRAMA DE
INGENIERÍA CIVIL
PROFESOR: Ing. Héctor Matamoros
ASIGNATURA: ACUEDUCTOS Y ALCANTARILLADO

Adicionalmente para completar el diseño de la bocatoma como se mencionó

anteriormente es necesario para no afectar las condiciones del río, por esta razón para
disipar la energía de exceso almacenada se diseña un cuenco disipador de energía y
una estructura de control de aguas abajo. Dicho proceso se inicia planteando el valor
de ΔB, el cual es la diferencia entre la cota de la sección del río para su punto más
bajo donde se definió el inicio del dique y la cota de la solera que va a corresponder
por ser un sistema horizontal con una altura de vertedero con este punto de entrega,
relacionado con la inclinación del río en ese sector.

Delta B: definido por el diseñador

Z: La diferencia de energía entre el cuenco y aguas arriba del dique
V1: velocidad aguas arriba(m/s)
Y1:lamina de aguas arriba(m)
NF: Numero de froude de Y1

A partir de esto se realiza un proceso iterativo hasta obtener un y1 con tendencia a cierto
valor.
DELTA B 0.050 m
Z 1.929 m
V1 5.841 m/s
Y1 0.089 m
F1 6.257

V1 6.014 m/s
Y1 0.086 m
F1 6.537

V1 6.018 m/s
Y1 0.086 m
F1 6.543

V1 6.018 m/s
Y1 0.086 m
F1 6.544

Posteriormente se obtiene el valor de la altura del obstáculo a partir de la ecuación.

Luego la lámina de Y2.

36
ESCUELA COLOMBIANA DE INGENIERÍA JULIO GARAVITO PROGRAMA DE
INGENIERÍA CIVIL
PROFESOR: Ing. Héctor Matamoros
ASIGNATURA: ACUEDUCTOS Y ALCANTARILLADO

la profundidad critica mediante:

Posteriormente se plantea ecuación de fuerza específica para hallar la lámina y3, que se
encuentra por encima del obstáculo.

Por ultimo se plantea ecuación de energía para hallar la lámina y4, después del
obstáculo.

Ya habiendo calculo los anteriores valores se halla la longitud del cuenco y su radio de
curvatura, a partir de las siguientes ecuaciones.

37
ESCUELA COLOMBIANA DE INGENIERÍA JULIO GARAVITO PROGRAMA DE
INGENIERÍA CIVIL
PROFESOR: Ing. Héctor Matamoros
ASIGNATURA: ACUEDUCTOS Y ALCANTARILLADO

hobsta 0.290 m
Y2 0.974 m
Yc 0.302 m
Y3 0.667 m
Y4 0.668 m
X 4.788 m
X 5.842 m
X 6.124 m
X 6.5 m
R 2.367 m

Apartir de las longitudes obtenidas, se opto por una de 6.5m ya que es mayor a las
halladas. Así mismo se hallaron las velocidades de cada lamina, para conocer con que
velocidad se entrega el flujo al rio.

V2 0.533 m/s
Vc 1.720 m/s
V3 0.778 m/s
V4 0.778 m/s

2. OGEE

El vertedero tipo OGEE se diferencia en sus radios de curvatura R2 y R1, y en la forma

de la curva de descarga. Este también funciona para paramento vertical, para 3 a 1 ,3 a
2, y 3 a 3. En su diseño los coeficientes de descarga de este vertedero son mayores que
del tipo WES. Se procede a diseñar la geometría definida se un dique de forma vertical,
de acuerdo con las tablas de diseño se obtendrán las distancias X, Y, R1, y R2.

Al igual que con el diseño del vertedero tipo WES se diseñará un empate con el cuenco
disipador de energía mediante una curva de radio constante y una estructura de control
de aguas abajo (el vertedero de pared gruesa) que permitirá el desarrollo del resalto
hidráulico para poder entregar las aguas al río de una manera semejante a la condición
natural. Nuevamente se parte del estudio hidrológico realizado donde se determinaron
los valores (Qmax, Qm y Qmin) en el punto de la bocatoma, y también de la sección del
río luego del levantamiento en HecRAS seleccionada como apta para establecer el punto
de la bocatoma.

CAUDAL MAXIMO
Q diseño: 7.79 m3/s
L(ancho presa) 15 m
h adoptada: 1.5 m

38
ESCUELA COLOMBIANA DE INGENIERÍA JULIO GARAVITO PROGRAMA DE
INGENIERÍA CIVIL
PROFESOR: Ing. Héctor Matamoros
ASIGNATURA: ACUEDUCTOS Y ALCANTARILLADO

Co 2.1
Ho (m) 0.394
P/Ho 3.809

Co 2.18
Ho (m) 0.384
P/Ho 3.905

Co 2.18
Ho (m) 0.384
P/Ho 3.905

Factor de correcion 1
Correguido 2.18

39
ESCUELA COLOMBIANA DE INGENIERÍA JULIO GARAVITO PROGRAMA DE
INGENIERÍA CIVIL
PROFESOR: Ing. Héctor Matamoros
ASIGNATURA: ACUEDUCTOS Y ALCANTARILLADO

CAUDAL MINIMO
Qmin (m3/s) 0.04
Lo (m) 15
Altura paramento p (m) 1.5

Co 2.18
He (m) 0.011
He/Ho 0.029
C/Co 0.925
Valor c 2.0165

He 0.012
He/Ho 0.030
C/Co 0.93
Valor c 2.0274

He 0.012
He/Ho 0.030
C/Co 0.93
Valor c 2.0274

Hemin 0.0116012

40
ESCUELA COLOMBIANA DE INGENIERÍA JULIO GARAVITO PROGRAMA DE
INGENIERÍA CIVIL
PROFESOR: Ing. Héctor Matamoros
ASIGNATURA: ACUEDUCTOS Y ALCANTARILLADO

CAUDAL MEDIO
Qmin (m3/s) 1.795
Lo (m) 15
Altura paramento p (m) 1.5
Co 2.18
He 0.144
He/Ho 0.376
C/Co 0.9
Valor c 1.962

He 0.155
He/Ho 0.403
C/Co 0.91
Valor c 1.9838

He 0.154
He/Ho 0.400
C/Co 0.91
Valor c 1.9838

Hemin 0.1537908

DISEÑO DE LA GEOMETRIA

41
ESCUELA COLOMBIANA DE INGENIERÍA JULIO GARAVITO PROGRAMA DE
INGENIERÍA CIVIL
PROFESOR: Ing. Héctor Matamoros
ASIGNATURA: ACUEDUCTOS Y ALCANTARILLADO

Va 0.275 m/s
ha 0.004 m
ho 0.380 m

Va 0.276 m/s
ha 0.004 m
ho 0.380 m

Va 0.276 m/s
ha 0.004 m
ho 0.380 m

ha/Ho 0.01

LECTURAS
Xo/Ho 0.28
Xo (m) 0.108
yo/Ho 0.120
Yo (m) 0.046
R1/Ho 0.525
R1 (m) 0.200
R2/Ho 0.225
R2 (m) 0.086
k 0.500
n 1.865

42
ESCUELA COLOMBIANA DE INGENIERÍA JULIO GARAVITO PROGRAMA DE
INGENIERÍA CIVIL
PROFESOR: Ing. Héctor Matamoros
ASIGNATURA: ACUEDUCTOS Y ALCANTARILLADO

X(m) y(m)
0 0
0.1 -0.015747
0.2 -0.05736
0.3 -0.122185
0.4 -0.208944
0.5 -0.316786
0.6 -0.445081
0.7 -0.593328
0.8 -0.761114
0.9 -0.94809
1 -1.15395
1.1 -1.378429
1.2 -1.621288
1.3 -1.882312
1.4 -2.161304
1.5 -2.458088
1.6 -2.772496
1.7 -3.104377
1.8 -3.453587
1.9 -3.819993
2 -4.203471
2.1 -4.603903
2.2 -5.021176
2.3 -5.455186

GOLA
0
0 0.5 1 1.5 2 2.5
-1

-2
y(m)

-3

-4

-5

-6
x(m)

43
ESCUELA COLOMBIANA DE INGENIERÍA JULIO GARAVITO PROGRAMA DE
INGENIERÍA CIVIL
PROFESOR: Ing. Héctor Matamoros
ASIGNATURA: ACUEDUCTOS Y ALCANTARILLADO

hobsta 0.290 m
Y2 0.974 m
Yc 0.302 m
Y3 0.646 m
Y4 0.951 m
X 6.206 m
X 5.844 m
X 6.721 m
X 6.000 m
R 1.809 m
V2 0.533 m/s
Vc 1.720 m/s
V3 0.803 m/s
V4 0.546 m/s

DISEÑO: BOCATOMA

1. POR ORIFICIO

Q diseño 0.005 m3/s

Viscosidad
cinematica 1.00E-06 m²/s

H 0.5
H/a 2.5
Re 6.24E+05
Cd 0.6

separacion(b) 1"
ancho(s) 0.75 "
# de barras 0.474533343
# de barras 1
L total(m) 0.070
Vn(m/s) 0.984251969
Re 1.96E+05
Cd 0.6
44
ESCUELA COLOMBIANA DE INGENIERÍA JULIO GARAVITO PROGRAMA DE
INGENIERÍA CIVIL
PROFESOR: Ing. Héctor Matamoros
ASIGNATURA: ACUEDUCTOS Y ALCANTARILLADO

Camara de capacitación:

Flujo Máximo Flujo minimo

Energia H 1.88 1.511 m
Carga a 0.868 0.5 m
Q 0.020 0.015 m3/s
V0 1.946 1.480 m2/s

X(m) Y(m) Y(m) Vy(m/s) V(m/s) hc(m/s) H/Lb a b Y/a Cp Ypsicina Y/a
0 0 0 0 1.9462172 0.19
0.1 -0.01 -0.02 2.01 2.80 0.40 38.77 0.5 0.4 5.8 0.049 -0.020 0.098
0.2 -0.05 -0.09 2.19 2.93 0.44 9.69 0.5 0.4 5.8 0.049 -0.022 0.108
0.3 -0.12 -0.20 2.46 3.14 0.50 4.31 0.5 0.4 5.8 0.049 -0.025 0.123
0.4 -0.21 -0.36 2.80 3.41 0.59 2.42 0.5 0.4 5.8 0.049 -0.029 0.146
0.5 -0.32 -0.56 3.18 3.73 0.71 1.55 0.55 0.225 5.8 0.149 -0.106 0.529
0.6 -0.47 -0.81 3.60 4.09 0.85 1.08 0.65 0.174 5.8 0.237 -0.202 1.009
0.7 -0.63 -1.10 4.03 4.47 1.02 0.79 0.65 0.067 5.8 0.441 -0.450 2.248
0.8 -0.83 -1.43 4.48 4.88 1.21 0.61 0.65 0.067 5.8 0.441 -0.535 2.676
0.90 -1.05 -1.81 4.93 5.30 1.43 0.48 0.98 0.07 5.8 0.653 -0.936 4.682
0.99 -1.27 -2.20 5.36 5.70 1.66 0.40 0.98 0.07 5.8 0.653 -1.081 5.405
1.1 -1.57 -2.71 5.87 6.19 1.95 0.32 0.98 0.07 5.8 0.653 -1.274 6.372
1.2 -1.86 -3.22 6.35 6.64 2.25 0.27 0.98 0.07 5.8 0.653 -1.469 7.344
1.3 -2.19 -3.78 6.83 7.10 2.57 0.23 0.98 0.07 5.8 0.653 -1.680 8.400
1.4 -2.54 -4.39 7.32 7.57 2.92 0.20 0.98 0.07 5.8 0.653 -1.908 9.541
1.5 -2.91 -5.04 7.80 8.04 3.30 0.17 0.98 0.07 5.8 0.653 -2.153 10.766
1.6 -3.31 -5.73 8.29 8.52 3.70 0.15 0.98 0.07 5.8 0.653 -2.415 12.075
1.7 -3.74 -6.47 8.78 9.00 4.13 0.13 0.98 0.07 5.8 0.653 -2.694 13.469
1.8 -4.19 -7.25 9.28 9.48 4.58 0.12 0.98 0.07 5.8 0.653 -2.990 14.948
1.9 -4.67 -8.08 9.77 9.96 5.06 0.11 0.98 0.07 5.8 0.653 -3.302 16.511
2 -5.18 -8.96 10.26 10.45 5.56 0.10 0.98 0.07 5.8 0.653 -3.632 18.159
2.1 -5.71 -9.87 10.76 10.93 6.09 0.09 0.98 0.07 5.8 0.653 -3.978 19.891
2.2 -6.27 -10.84 11.25 11.42 6.65 0.08 0.98 0.07 5.8 0.653 -4.341 21.707

Chorro libre
0
0 0.5 1 1.5 2 2.5

-2

-4

-6

-8

-10

-12

Q excesos 0.010 m3/s

L 0.4 m
H 0.056 m

formula francis
H= 0.056076838

45
ESCUELA COLOMBIANA DE INGENIERÍA JULIO GARAVITO PROGRAMA DE
INGENIERÍA CIVIL
PROFESOR: Ing. Héctor Matamoros
ASIGNATURA: ACUEDUCTOS Y ALCANTARILLADO

V= 1 m/s
D= 0.112837917 m
D= 0.2 m
S= 0.5 m
h= 0.85 m

SUMERGENCIA DE 2.5D:
Le= 0.388784632 m
H= 0.057570327 m
L= 0.400298698 m

C 3.2963891
Le 0.388784632 m
H 0.057570327 m
L 0.400298698 m

V= 0.437 m/s
hv= 0.009718329 m
A= -0.38
B= -0.01
C= 0.07
m= -0.04
D= 0.57

X(m) Y(m) X/H

0 0.037 0 mal
0.1 -0.030 1.737005933 ok
0.2 -0.230 3.474011866 ok
0.3 -0.563 5.211017798 ok
0.4 -1.029 6.948023731 ok
0.5 -1.628 8.685029664 ok
0.6 -2.360 10.4220356 ok
0.7 -3.225 12.15904153 ok

46
ESCUELA COLOMBIANA DE INGENIERÍA JULIO GARAVITO PROGRAMA DE
INGENIERÍA CIVIL
PROFESOR: Ing. Héctor Matamoros
ASIGNATURA: ACUEDUCTOS Y ALCANTARILLADO

QMD 0.005 m3/s

Temperatura 20 °C
Viscosidad 1.00E-06 m²/s
% Ampliación 100 %
Cauldal derivado 0.01 m3/s
Altura orificio (a) 0.2 m
Caarga sobre orificio (H) 0.5 m

Re 624417.42
Carga adopatada Cd 0.6

Separación 1
Ancho 0.75

Numero de barras 1
B total 0.070
B neto 0.0508
Velocidad neta 0.984251969
Re 196261.6089
Cd 0.6

2. POR VERTEDERO LATERAL

47
ESCUELA COLOMBIANA DE INGENIERÍA JULIO GARAVITO PROGRAMA DE
INGENIERÍA CIVIL
PROFESOR: Ing. Héctor Matamoros
ASIGNATURA: ACUEDUCTOS Y ALCANTARILLADO

Angulo de rejilla 90 °
separacion(b) 1"
ancho(s) 0.75 "

Velocidad sobre
rejilla 0.7 m/s
Cf 1.79
h(perdida de
rejilla) 0.030474897 m

h(perdida de
rejilla+taponami
ento) 0.091424691 m

Ampliacion de
caudal 100 %
QMD 0.05 m3/s
Qdiseño 0.1 m3/s

H(supuesta) 0.15 m

Sumergencia 0.390502063
Qsumergido 0.111371668 m3/s

formula francis
Le 1.04188554 m
L 1.07188554 m

#varillas 42
Lvertedero 1.8923
L neta 1.092 m

Verificacion de Engles
L 1.277 m

Qcaptado
sumergido 0.101949778 m3/s
Vneta rejilla 0.622290047 m/s
H1 0.058575309

48
ESCUELA COLOMBIANA DE INGENIERÍA JULIO GARAVITO PROGRAMA DE
INGENIERÍA CIVIL
PROFESOR: Ing. Héctor Matamoros
ASIGNATURA: ACUEDUCTOS Y ALCANTARILLADO

CAMARA DE CAPTACION:

LB 1.5 fijo
H(aguas
maxima) 0.51805806 m
Qmax captado 0.638249796 m3/s

hx3 (por
obstrucciones) 0.23740375 m
Sumergencia 0.541742965
H1 0.280654309 m

Qmax captado 0.586672089 m3/s

hx3 (por obstrucciones)0.200584356 m
Sumergencia 0.61281491
H1 0.317473703 m

49
ESCUELA COLOMBIANA DE INGENIERÍA JULIO GARAVITO PROGRAMA DE
INGENIERÍA CIVIL
PROFESOR: Ing. Héctor Matamoros
ASIGNATURA: ACUEDUCTOS Y ALCANTARILLADO

Qmax captado 0.555839795 m3/s

hx3 (por obstrucciones)0.180055121 m
Sumergencia 0.652442197
H1 0.338002938 m

Qmax captado 0.536318609 m3/s

hx3 (por obstrucciones)0.167630074 m
Sumergencia 0.676426085
H1 0.350427985 m

Qmax captado 0.523552605 m3/s

hx3 (por obstrucciones)0.159744846 m
Sumergencia 0.691646828
H1 0.358313214 m

Qmax captado 0.515036747 m3/s

hx3 (por obstrucciones)0.154590442 m
Sumergencia 0.7015963
H1 0.363467618 m

Qmax captado 0.509283083 m3/s

hx3 (por obstrucciones)0.151155762 m
Sumergencia 0.708226214
H1 0.366902298 m

Qmax captado 0.50536278 m3/s

hx3 (por obstrucciones)0.148837618 m
Sumergencia 0.712700893
H1 0.369220442 m

Qmax captado 0.502676503 m3/s

hx3 (por obstrucciones)0.147259519 m
Sumergencia 0.715747075
H1 0.370798541 m

Qmax captado 0.500828733 m3/s

hx3 (por obstrucciones)0.146178896 m
Sumergencia 0.717832985
H1 0.371879163 m

Qmax captado 0.499554397 m3/s

hx3 (por obstrucciones)0.145435952 m
Sumergencia 0.71926708
H1 0.372622108 m

Qmax captado 0.498673953 m3/s

hx3 (por obstrucciones)0.144923754 m
Sumergencia 0.720255769
H1 0.373134306 m

Qmax captado 0.498064896 m3/s

hx3 (por obstrucciones)0.144569964 m
Sumergencia 0.720938684
H1 0.373488096 m

Qmax captado 0.497643213 m3/s

hx3 (por obstrucciones) 0.14432527 m
Sumergencia 0.721411014
H1 0.37373279 m

Qmax captado 0.497351087 m3/s

hx3 (por obstrucciones)0.144155876 m
Sumergencia 0.721737992
H1 0.373902184 m

50
ESCUELA COLOMBIANA DE INGENIERÍA JULIO GARAVITO PROGRAMA DE
INGENIERÍA CIVIL
PROFESOR: Ing. Héctor Matamoros
ASIGNATURA: ACUEDUCTOS Y ALCANTARILLADO

VERTEDERO DE EXCESOS:

Qexcesos 0.397351087 m3/s

H sobre vertedero 0.315326874 m3/s
C 1.84
L 1.22 m

ADUCCION(PRESION)
V(m/s) 1.5
D(m) 0.291346248
D(m) 0.3
S(m) 0.75
(m) 1.15
Hsobre vertedero/h 0.274197282

VERTEDERO CON DOS CONTRACCIONES:

Le 1.22 m
L 1.282658883 m
Q 0.440438449 m3/s

V 1.0332 m/s

Hv 0.0544 m

A -0.3819
B -0.0066
C 0.0723
m -0.0354
D 0.5682
Hv/H 0.1726

51
ESCUELA COLOMBIANA DE INGENIERÍA JULIO GARAVITO PROGRAMA DE
INGENIERÍA CIVIL
PROFESOR: Ing. Héctor Matamoros
ASIGNATURA: ACUEDUCTOS Y ALCANTARILLADO

X(m) Y(m) X/Y

0 0.20199402 0 no cumple
0.04 0.19979432 0.20020589 no cumple
0.08 0.19371942 0.4129684 no cumple
0.12 0.18376933 0.65299252 cumple
0.16 0.16994405 0.94148634 cumple
0.2 0.15224357 1.31368439 cumple
0.24 0.1306679 1.83671743 cumple
0.28 0.10521703 2.66116627 cumple
0.32 0.07589096 4.21657591 cumple
0.36 0.0426897 8.43294709 cumple
0.4 0.00561325 71.2599903 cumple
0.44 -0.0353384 -12.451044 cumple
0.48 -0.0801652 -5.9876322 cumple
0.52 -0.1288673 -4.0351591 cumple
0.56 -0.1814445 -3.0863429 cumple
0.6 -0.2378969 -2.5221005 cumple
0.64 -0.2982246 -2.1460338 cumple
0.68 -0.3624274 -1.8762379 cumple
0.72 -0.4305054 -1.6724529 cumple
0.76 -0.5024586 -1.5125624 cumple
0.8 -0.578287 -1.3833961 cumple

52
ESCUELA COLOMBIANA DE INGENIERÍA JULIO GARAVITO PROGRAMA DE
INGENIERÍA CIVIL
PROFESOR: Ing. Héctor Matamoros
ASIGNATURA: ACUEDUCTOS Y ALCANTARILLADO

FONDO DE
CAPTACION 0.51805806 m
ANCHO DE
CAMARA 0.9 m

3. POR REJILLA DE FONDO

BOCATOMA SUMERGIDA:

QMD 0.010 m^3/s

Q1 Qmin 0.0380 m^3/s qmin 0.00253 m^3/s
Q2 Qdiseño 0.0163 m^3/s qmin 0.02331 m^3/s
Q3 qmin-qdiseño 0.0217 m^3/s qmin-qdiseño -0.02077 m^3/s

53
ESCUELA COLOMBIANA DE INGENIERÍA JULIO GARAVITO PROGRAMA DE
INGENIERÍA CIVIL
PROFESOR: Ing. Héctor Matamoros
ASIGNATURA: ACUEDUCTOS Y ALCANTARILLADO

En la cresta del vertedero tiene flujo critico

ϒc 0.00868222 m
la altura al inicio de la rejilla debe ser ligeramente menor a Yc
Y1 0.007814 m
la velocidad inicial es
V1 0.3242044 m/s
la energia especifica es
E1 0.0132 m
Y1/E1 0.5932
para la rejilla se asignan anchos de barras y separación entre ellos.
separacion(b) 1"
ancho(s) 0.75 "
e 0.571
Según la curvatura del dique la pendiente es
1 a 0.082197598
Cd 0.496

Y2= 0.03530586

La rejilla de la captacion tota es

Lr= 0.0176

Q max
qmax= 2.80298756
En la cresta del vertedero se tiene flujo critico Qmax= 42.04
Yc= 0.92878847
La altura al inicio de la rejilla debe ser ligeramente menor a Yc
Y1= 0.83590963
La velocidad inicial es
V1= 3.353218423
La energia especifica es
E1= 1.40923584
Y1/E 0.59316517 ok

Para el ancho diseñado la altura al final de la reja es=

Y2= 0.5361 #¡REF!
q= 1.2291

Q captado: 1.5739

CANAL DE CAPTACION:

54
ESCUELA COLOMBIANA DE INGENIERÍA JULIO GARAVITO PROGRAMA DE
INGENIERÍA CIVIL
PROFESOR: Ing. Héctor Matamoros
ASIGNATURA: ACUEDUCTOS Y ALCANTARILLADO

Q= 1.1017163
q= 0.85084473

b= 0.3
Yc= 1.11207749
H2= 1.22328524

S= 0.001
H1= 1.93406472
B libre= 0.64468824
H total= 2.57875295

Con V =3.002m/s
Lb= 1.42328524
H/Lb= 0.8595
Cp= 0.4214
Deair= 0.5786

X(m) Y(m) Ye(m) Vy(m/s) V(m/s) hc(m/s)

0.000 0.00000 1.22329 0 0.000 0.0000
0.100 -0.00544 1.21784 3.01978767 4.258 0.9241
0.200 -0.02176 1.20152 3.00206979 4.246 0.9187
0.300 -0.04896 1.17432 3.05491517 4.283 0.9350
0.400 -0.08704 1.13624 3.14101566 4.345 0.9622
0.500 -0.13601 1.08728 3.25773559 4.430 1.0003
0.600 -0.19585 1.02744 3.40192476 4.537 1.0492
0.623 -0.21121 1.01208 3.57025658 4.665 1.1090
0.800 -0.34818 0.87511 3.61219898 4.697 1.1244
0.900 -0.44066 0.78262 3.96663324 4.975 1.2613
1.000 -0.54403 0.67926 4.18903169 5.154 1.3537
1.100 -0.65827 0.56501 4.42438502 5.347 1.4571
1.200 -0.78340 0.43989 4.6707353 5.552 1.5713
1.300 -0.91941 0.30388 4.92643306 5.769 1.6963
1.400 -1.06629 0.15699 5.19009694 5.996 1.8323
1.500 -1.22406 -0.00078 5.46057313 6.231 1.9791
1.600 -1.39271 -0.16942 5.73689818 6.475 2.1368
1.700 -1.57224 -0.34895 6.01826649 6.725 2.3054
1.790 -1.74222 -0.51894 6.30400281 6.982 2.4849
1.900 -1.96394 -0.74065 6.5630906 7.217 2.6548

55
ESCUELA COLOMBIANA DE INGENIERÍA JULIO GARAVITO PROGRAMA DE
INGENIERÍA CIVIL
PROFESOR: Ing. Héctor Matamoros
ASIGNATURA: ACUEDUCTOS Y ALCANTARILLADO

1.50000

1.00000

0.50000

0.00000
0.000 0.200 0.400 0.600 0.800 1.000 1.200 1.400 1.600 1.800 2.000

-0.50000

-1.00000

-1.50000

-2.00000

-2.50000

B= 1.166 m
Y/B= 1 m
Y= 1.2 m
Vf^2/2g= 1.797 m
Hm= 1.924 m

Vertedero de excesos
El caudal de excesos es 1.558 m^3/s
L 1.80 m

v m/s 1.431 m/s

hv m 0.104 m
A -0.411
B 0.034
C 0.124
m -0.150
D 0.560
hv/H 0.058

56
ESCUELA COLOMBIANA DE INGENIERÍA JULIO GARAVITO PROGRAMA DE
INGENIERÍA CIVIL
PROFESOR: Ing. Héctor Matamoros
ASIGNATURA: ACUEDUCTOS Y ALCANTARILLADO

X(m) Y(m) X/H Revisión

0 0.414
0.1 0.410 0.165 no cumple
0.2 0.393 0.331 no cumple
0.3 0.363 0.496 no cumple
0.4 0.319 0.661 cumple
0.5 0.261 0.827 cumple
0.6 0.190 0.992 cumple
0.7 0.105 1.158 cumple
0.8 0.006 1.323 cumple
0.9 -0.106 1.488 cumple
1 -0.231 1.654 cumple
1.1 -0.370 1.819 cumple
1.2 -0.523 1.984 cumple
1.3 -0.689 2.150 cumple
1.4 -0.869 2.315 cumple
1.5 -1.063 2.480 cumple
1.6 -1.270 2.646 cumple
1.7 -1.490 2.811 cumple
1.8 -1.724 2.976 cumple
1.9 -1.972 3.142 cumple
2 -2.233 3.307 cumple

57
ESCUELA COLOMBIANA DE INGENIERÍA JULIO GARAVITO PROGRAMA DE
INGENIERÍA CIVIL
PROFESOR: Ing. Héctor Matamoros
ASIGNATURA: ACUEDUCTOS Y ALCANTARILLADO

DISEÑO: COMPUERTAS DE LIMPIEZA

CAUDAL DE DISEÑO

Qmedio 1.884 m3/s h

QMD 0.05 m3/s Hd
Ampliacion
de caudal 100 %
Qdiseño 1.784
Hd 0.14715434
y1 1.64715434

L(ancho de limpieza) 0.5 m

a(apertura de compuerta)
0.3 m

y1/a 5.49051448
Cd 0.57
a 1.10111668

y1/a 1.49589446
Cd 0.55
a 1.14115728

y1/a 1.44340694
Cd 0.55
a 1.14115728

H (altura
compuerta) 2.17896178 m
Fuerza
hidrostatica 8655.01434 N
Vflujo 3.12665051 m/s

58
ESCUELA COLOMBIANA DE INGENIERÍA JULIO GARAVITO PROGRAMA DE
INGENIERÍA CIVIL
PROFESOR: Ing. Héctor Matamoros
ASIGNATURA: ACUEDUCTOS Y ALCANTARILLADO

SELECCIÓN DE BOCATOMA

Se escoge teoricamente a traves de varios diseños realizados para cada componente,

teniendo como resultado que la bocatoma de vertedero lateral en cuanto a costos resulta
ser la mas economica en comparacion a las demas, ademas de ofrecer datos racionales y
concisos.

ESQUEMA: BOCATOMA

59
ESCUELA COLOMBIANA DE INGENIERÍA JULIO GARAVITO PROGRAMA DE
INGENIERÍA CIVIL
PROFESOR: Ing. Héctor Matamoros
ASIGNATURA: ACUEDUCTOS Y ALCANTARILLADO

CÁMARA
CAPTACIÓN

REJILLA

DIQUE

T. excesos

Ilustración 27. Esquema Bocatoma

60
ESCUELA COLOMBIANA DE INGENIERÍA JULIO GARAVITO PROGRAMA DE
INGENIERÍA CIVIL
PROFESOR: Ing. Héctor Matamoros
ASIGNATURA: ACUEDUCTOS Y ALCANTARILLADO

Ilustración 28. Vista 3D Boccatoma

Ilustración 29. Vista 3D 2

61
ESCUELA COLOMBIANA DE INGENIERÍA JULIO GARAVITO PROGRAMA DE
INGENIERÍA CIVIL
PROFESOR: Ing. Héctor Matamoros
ASIGNATURA: ACUEDUCTOS Y ALCANTARILLADO

Ilustración 30. Vista 3D 3

SELECCIÓN DE BOCATOMA

Se escoge teoricamente a traves de varios diseños realizados para cada componente,

teniendo como resultado que la bocatoma de vertedero lateral en cuanto a costos resulta
ser la mas economica en comparacion a las demas, ademas de ofrecer datos racionales y
concisos.

62
ESCUELA COLOMBIANA DE INGENIERÍA JULIO GARAVITO PROGRAMA DE
INGENIERÍA CIVIL
PROFESOR: Ing. Héctor Matamoros
ASIGNATURA: ACUEDUCTOS Y ALCANTARILLADO

REFERENCIAS BIBLIOGRAFICAS

- Centro Administrativo Municipal (2018, marzo 13), Acuerdo N°4: Por el cual
se adopta la revisión ordinaria del Esquema de Ordenamiento Territorial del
Municipio de Gramalote, Gramalote, Norte de Santander, Colombia.

- Consejo Municipal de Gramalote (2020), Acuerdo No. Por el cual se adopta el

plan de desarrollo “Entre todos podemos fortalecer nuestro municipio” 2020-
2023, Gramalote, Norte de Santander, Colombia.

- Vargas, M. (2022), Informe de Gestión del Año 2022, Gramalote, Norte de

Santander, Colombia.

- Vargas, M. (2021), Informe de Gestión del Año 2021, Gramalote, Norte de

Santander, Colombia.

- Impactos y Valoraciones Consultores (2020, diciembre), Formulación del plan

de gestión integral de residuos sólidos-PGIRS-del municipio de Gramalote (NS),
Contrato de Consultoría No. 134, Gramalote, Norte de Santander, Colombia.

- Viceministro de Agua y Saneamiento Básico, Ministerio de Vivienda, Ciudad y

Territorio, (2010), Reglamento Técnico del Sector de Agua Potable y
Saneamiento Básico – RAS – Titulo B. Sistemas de Acueducto, Recuperado de:
[Link]
[Link]
- COLOMBIA. Censo Municipal: Gramalote 2018, [base de datos en línea].
[Censo 2018]. Disponible en: [Link]
por-tema/demografia-y-poblacion/censo-nacional-de-poblacion-y-vivenda-2018

63