Universidad de Boyacá-Ingeniería Sanitaria

Curso Estaciones Sostenibles de Tratamiento de Aguas Residuales

Jaime Díaz Gómez

Guia_1_Diseño Tratamiento Preliminar-Cribado

1 Diseño de un sistema de cribado

El siguiente ejemplo de diseño de un sistema de cribado está basado en un ejemplo de Qasim

et al (2018)1.

Se deben diseñar dos rejas de barras idénticas que puedan operar en condiciones máximas de
flujo que serán equipadas con un equipo mecánico para limpieza. El ángulo de inclinación
de la reja es =85 °C.

Factores de diseño de rejas manuales y mecánicas para el tratamiento de aguas

residuales

Factor de diseño Manual Mecánica

Velocidad a través de la reja (m/s) 0.3-0.6 0.6-1.0
Tamaño de las barras
Ancho (mm), ab 5-15 5-15
Profundidad (mm) 25-50 50-75
Abertura entre barras (mm), Ab 25-75 6-75
Inclinación (grados) 45-60 75-85
Velocidad de aproximación (m/s)
Máxima 0.3-0.6 0.6-1.0
Mínima 0.3-0.5
Perdida de carga permisible, reja 150 150
obstruida (mm)
Máxima perdida de carga, reja 800 800
obstruida (mm)
Fuente: Qasim et al (2018); Metcalf y Eddy (2014)

El RAS (2017, Articulo 186) recomienda el diseño de rejas mecánicas para caudales de
diseño mayores a 100 L/s. Además, distingue entre rejas finas y rejas gruesas por el espacio
entre las barras (Ab), así:

1
Qasim, S. 1999. Wastewater treatment plants. Planning, design and operation. Technomic Publishing Co.
Lancaster.

1
Durante la limpieza una de las rejas continuara operando para no interrumpir la operación de
la Planta de Tratamiento de Aguas Residuales (PTAR). En la siguiente figura se presenta una
rejilla manual típica de una pequeña PTAR.

Figura 1 Vista en planta y perfil de una rejilla manual

1.1 Criterios de diseño

Espaciamiento entre barras: 2.5 cm

La velocidad a través de las rejas en diferentes flujos debe ser

Velocidad en flujo máximo: 0.9 m/s

Velocidad en flujo medio no menor a 0.6 m/s
Velocidad en flujo mínimo no menor a 0.4 m/s

Los caudales de diseño son:

Caudal máximo: 1.321 m3/s

Caudal medio: 0.440 m3/s
Caudal mínimo: 0.220 m3/s

El caudal mínimo podría también calcularse con

C ∗ P ∗ q ∗ k3
Q=
86400
donde

q: dotación neta [Lt/hab*dia]

C: relación agua residual agua, se puede tomar como 1.0
k3: coeficiente. Se asume un valor de 0.5
P: población

1.2 Calculo del conducto de entrada a la PTAR

El caudal del conducto a tubo lleno se calcula con la siguiente expresión

2
 0.312 8⁄ 1⁄
Q= D 3S 2
n
donde:

Q: caudal [L3/t]
S: pendiente [L/L]
n: coeficiente de rugosidad empleado en la ecuación de Manning

La velocidad cuando el conducto se encuentra a tubo lleno se calcula con

1 2⁄ 1⁄
V= R 3S 2
n
donde:

R: radio hidráulico (A/P), en una tubería circular que fluye a tubo lleno R=D/4,
[L]
A(D2/4): área del tubo fluyendo lleno [L2]
D: diámetro [L]
V: velocidad de la tubería fluyendo a tubo lleno [L/t]

Las características del ducto de entrada a la PTAR son:

Pendiente (S): 0.00047

Rugosidad (n): 0.015 (concreto)
Diámetro: 60” (1.53 m)

Empleando las anteriores expresiones

0.312 8 1
Q= (1.53) ⁄3 (0.00047) ⁄2
0.015
Q=1.40 m3/s

El radio hidráulico (R) es: R=D/4=0.38 m

1 2⁄ 1
V = n (0.38) 3 (0.00047) ⁄2 V=0.76 m/s

El caudal máximo de diseño (q) es 1.321 m3/s. Por lo tanto, la relación q/Q= 1.321/1.40=0.94

3
 Figura 1 Relaciones hidráulicas en conductos circulares de alcantarillado
Fuente: Qasim (1999)

Tomando en cuenta las relaciones hidráulicas de conductos circulares presentada en la figura

1 la relación v/V es de 1,18. Por lo tanto la velocidad en el conducto de entrada en el caudal
máximo de diseño v=0.89 m/s.

Igualmente, la relación d/D=0.77. Por lo tanto, la profundidad del flujo en el flujo máximo
de diseño (d) es

d=D*o.77=1.53*0.77=1.18 m

1.3 Diseño

1.3.1 Paso 1. Condiciones de flujo en el conducto de entrada

Pendiente (S): 0.00047

Rugosidad (n): 0.015 (concreto)
Diámetro: 60” (1.53 m)
Velocidad en el flujo máximo de diseño: 0.897 m/s
Profundidad del flujo en el conducto en el flujo máximo de diseño: 1.18 m

4
1.3.2 Paso 2. Diseño del canal de cribado

Este se diseña para el flujo máximo. Se calculan las velocidades a través de la reja y el canal,
y la profundidad se calculan

1. Calcular los espaciamientos entre barras y las dimensiones del canal de la reja. Este canal
se diseña para flujo máximo. Las velocidades a través de la reja y el canal deben ser
corroboradas para el flujo promedio y mínimo.

Asumir que la profundidad de flujo en la reja es igual a la profundidad en el canal de

aproximación. Para el ejercicio es 1,18 m

a. Área de los espaciamientos en la reja:

3
Qd 1,321 m ⁄s
Ae = = m = 1,47 m2
vr 0,9 ⁄s
donde:

vr: velocidad a través de la reja [L/t]

b. Espaciamiento total entre barras

Area 1,47 m2
= = 1.25 m
profundidad del flujo 1,18 m

c. Se adoptan 50 separaciones cada 25 mm

El espaciamiento total en la reja (et) es

Separacion ∗ espaciamiento 50 ∗ 25 mm
et = = = 1.25 m
1000 1000

d. El número total de barras (Nb) es 49 (espaciamientos -1)

e. El ancho de las barras (ab) es de 10 mm
f. El ancho del canal de la reja es

espaciamiento total + ancho de la barra ∗ # de barras 1,25 m + 10 mm ∗ 49

Ac = =
1000 mm/m 1000 mm/m
= 1.74 m

La vista en planta, la sección longitudinal y la sección transversal del cribado se muestra a

continuación

5
 Figura 2 Vista en planta y corte transversal del cribado

Asumiendo condiciones de flujo uniforme el tiempo hidráulico de retención en el canal de

cribado puede tener un valor típico de 3 segundos. La longitud del canal de cribado se puede
calcular de acuerdo con

Q
L = ( ) ∗ THR
A

En rejillas manuales se recomienda una velocidad de aproximación alrededor de 0.45 m/s en

el flujo promedio para reducir la sedimentación de arena y otros materiales en el canal de
cribado (Metcalf and Eddy, 2014).

La cantidad de residuos sólidos generados en la rejilla (M) se puede calcular asumiendo un

valor de 20 litros de material por 1000 m3 de agua residual.

20
M (m3 ⁄d) = ∗Q
1000

6
 Figura 3 Seccional transversal de la reja

2. Calcular el coeficiente de eficiencia (CE)

espacimiento total 50 ∗ 25 mm
CE = = = 0.72
ancho del canal 1740 mm

3. Calcular la profundidad actual de flujo y la velocidad en el canal en el flujo máximo

a. la profundidad actual en el canal aguas arriba de la reja se calcula empleando la
ecuación de la energía

v1 2 v2 2
Z1 + d1 + = Z2 + d2 + + hL
2g 2g
donde:

Z1 y Z2: altura (m)

v1 y v2: velocidad en la sección 1 y 2
hL: altura total
d1 y d2: profundidad de flujo en la sección 1 y 2.

En este ejemplo para expansión con v1>v2:

v1 2 v2 2
hL = K e ( − )
2g 2g
Donde:

7
Ke: coeficiente de expansión. Su valor de pende de la configuración de la sección de
aproximación al sistema.

Para expansión con v1<v2:

v2 2 v1 2
hL = K e ( − )
2g 2g

b. Usando la ecuación de la energía entre la sección 1 en el conducto de entrada y la

sección 2 en el canal aguas arriba de la reja se determinan la profundidad y la velocidad
en el canal. La ecuación de la energía se escribe teniendo en cuanta las siguientes
condiciones

 El fondo del canal es horizontal

 El punto de referencia está en el fondo del canal (Z2=0)
 La parte más baja del tubo de entrada esta 8 cm arriba del punto de referencia
 Ke: 0,3
(0,87 m/s)2
0,08 + 1,18 +
2 ∗ 9,81 m/s2
2
1,321 m3 /s
( )
1,74m ∗ d2 m
= d2 m +
2 ∗ 9,81 m/s2
2
1,321 m3 /s
(0,87 m/s)2 ( )
1,74m ∗ d2 m
+ 0,3 2
−
2 ∗ 9,81 m/s 2 ∗ 9,81 m/s2
[ ]

c. Simplificando la anterior ecuación:

d2 3 − 1,288 d2 2 + 0,021 = 0

d. Solucionando por ensayo y error

Con Matlab

>>p=[1 -1.288 0.021]

>>roots(p)

1.2715
0.0165
m3
1,321 s m
d2 = 1,28 m y V2 = = 0,59
1,74m ∗ 1,28 m s

8
Según recomendación del RAS (2017) la velocidad de aproximación a la reja debe ser de 1,2
m/s para caudal máximo y 0,3 m/s para caudal mínimo con el fin de evitar sedimentación de
material en el canal de cribado.

4. Calcular la velocidad V a través de las aberturas en la reja

Q
V=
Area de la reja(espaciamiento total ∗ altura lamina de agua en flujo maximo)
m3
1,321 s
= = 0.83 m/s
1,25 m ∗ 1,28 m

5. Calcular la perdida de carga a través de la reja.

V 2 − V2 2 1
hL = ∗
2g 0.7

m 2 m 2
(0,83 s ) − (0,59 s ) 1
hL = ∗
19,62 0.7

hL = 0,025 m

6. Calcule la profundidad del flujo y la velocidad en el canal aguas abajo de la reja.

a. La profundidad y la velocidad en el canal por ensayo y error con la ecuación de la

energía:

v2 2 v3 2
d2 + = d3 + + hL
2g 2g
Donde:

d2 y v2: profundidad del flujo y velocidad en el canal aguas arriba de la reja

d3 y v3: profundidad del flujo y velocidad aguas abajo de la reja
hL: perdida de carga a través de la reja
2
m3
1,321 s
m 2 ( )
1,74 m ∗ d3 m
(0,59 s )
1,28 m + m = d3 m + m + 0,025 m
19,62 2 19,62 2
s s

b. Simplificando la ecuación:

d3 3 − 1, 273 d3 2 + 0,029 = 0

c. Solucionando esta ecuación por ensayo y error:

9
 m
d3 = 1,25 m y v3 = 0,61
s

7. Calcular la perdida de carga a través de la reja en una obstrucción del 50%

a. En un 50 % de obstrucción, el área libre a través de la reja se reduce a la mitad y la

pérdida de carga a través de la reja se obtiene de la ecuación de la energía

v2´2 v3 2
d´2 + = d3 + + h50
2g 2g

d´2 y v2´2 : profundidad del flujo y velocidad en el canal aguas arriba de la reja
h50: perdida de carga en la reja en un 50 % de obstrucción

b La velocidad y profundidad del flujo en el canal están definidos por la condición en

la salida del canal (el canal que transporta el agua residual depurada). Ya que las
condiciones en el canal de salida no cambian se puede asumir que d3 y v3 tienen el
mismo valor que el calculado cuando la reja está limpia.

c.
(Velocidad a traves de la abertura de la reja)2 − v2´2 1
h50 = ∗
2g 0.7

d. La velocidad a través de la abertura de la reja en una obstrucción del 50%

3
1.321 m ⁄s 2.114
= m/s
1.25 ∗ 0.5 ∗ d´2 m d´2

m
1.51
s
e. La velocidad en el canal aguas arriba de la reja
3
1.321 m ⁄s 0.759
v2´2 = = m/s
1.74 m ∗ d´2 m d´2

0.54 m/s
f.

2 2 2
0.759 2.114 0.759
( ´ ) [( ´ ) − ( ´ ) ]
d2 (0.61 m/s)2 d2 d2
d´2 m + = 1.25 m + + m/s
2 ∗ 9.81 m⁄s2 2 ∗ 9.81 m⁄s2 2 ∗ 0.7 ∗ 9.81 m⁄s 2

g Simplificando

10
 d2 ′´3 − 1, 269 d2 ′´2 − 0,254 = 0

h Solucionando por ensayo y error

d2 ′´ = 1.40
3
1.321 m ⁄s
v2 ′´ = = 0.54 m/s
1.74 m ∗ 1.40 m

i La pérdida de carga en el 50 % de obstrucción es

m
(1.51 s )2 − (0.54 m/s)2
h50 = = 0.15 m
2 ∗ 9.81 m⁄s2 ∗ 0.7

La longitud del canal de cribado se puede calcular asumiendo que el tiempo hidráulico de
retención en el canal puede estar entre 3 y 5 segundos. Teniendo en cuenta la velocidad en el
canal la longitud del canal se calcula con la siguiente expresión.

L= v∗t

1.3.3 Paso 3. Estructura efluente

Es importante diseñar una sección de control al final del canal de cribado para mantener una
profundidad normal ya que la salida de la cámara es de caída libre en un pozo húmedo. Esta
sección de control puede ser un vertedero proporcional o una canaleta Parshall. En este caso
emplearemos un vertedero proporcional. Una de las ventajas de este tipo de vertedero es que
mantiene una velocidad uniforme en el canal en condiciones de flujo variable.

El flujo a través de un vertedero proporcional es

Q = 1.57 Cd √2g L H 3⁄2

donde

Q: flujo a través del vertedero proporcional, m3/s

H: carga sobre el vertedero, m
Cd : coeficiente de descarga. El valor esta en un rango de 0.6-0.9. Debido al crecimiento
de bio película y obstrucciones se emplea comúnmente un valor de 0.6.
L: longitud de la abertura del vertedero en una altura H cresta del vertedero, m
g: aceleración de la gravedad, 9.81 m/s2

Sustituyendo los valores de Cd y g

Q = 4.17 ∗ [LH 1⁄2 ] ∗ H

11
La cresta del vertedero con flujo máximo (H) es de 1.25 (d3 en el canal de cribado)

La longitud de la abertura del vertedero en la cabeza máxima (L) es

Q 1.321 m3 ⁄s
L= = = 0.23 m
4.173 ∗ H 3⁄2 4.17 ∗ (1.25)3⁄2

La geometría del perfil del vertedero se calcula manteniendo el producto LH1/2 constante.
Esto permite mantener una velocidad casi uniforme bajo condiciones de flujo variable ya que
la profundidad en el canal debe ser proporcional al flujo en el canal o a la cabeza H sobre el
vertedero.
4.173[LH 1⁄2 ]H LH 1⁄2
L= = 1⁄2
4.173H 3⁄2 H

El valor constante de k=LH1/2 es

k = LH 1⁄2 = 0.23 √1.25 = 0.25

La velocidad en el canal aguas arriba del vertedero Sutro se calcula de acuerdo con:

Q
v
Ac * H
Cuando se tiene una altura d3=1.25 m, la velocidad en el canal de cribado es:

1.321
v  0.61 m/s
1.74*1.25

Un ejemplo de los cálculos para diferentes caudales en el canal de cribado se presenta a

continuación. Para Q=0.92 m3/s (asumido):

La carga sobre la cresta del vertedero (d) se calcula con la siguiente expresión

Q
dy
4.17*k

Para este caudal

0.92
dy  0.856 m
4.17*0.257

La longitud correspondiente a esa lámina de agua es

L  k * y 1 2

12
L  0.257*1.231 2  0.232 m

La velocidad en el canal de cribado es:

Q
v
Ac *d
0.92
v  0.6 m/s
1.74*0.88

A continuación, se presenta la Tabla que describe los cálculos para el diseño del vertedero
Sutro. Se observa una ligera variación en la velocidad para las diferentes condiciones de flujo.

Tabla 1. Diseño del vertedero Sutro

Q [m3/s] H [m] L superior [m] Ac [m] k

1.321 1.25 0.23 1.74 0.25

Q y=Q/(4.17*k) L=k*y-1/2 v=Q/(Ac*d)

1.321 1.25 0.23 0.61
0.92 0.87 0.24 0.61
0.44 0.42 0.35 0.61
0.22 0.21 0.50 0.61
0.15 0.14 0.60 0.61
0.09 0.09 0.78 0.61
0.053 0.05 1.01 0.61

Un esquema del vertedero se presenta en la siguiente figura

13
Figura 1 Detalle diseño vertedero Sutro y perfil hidráulico del canal de cribado en el flujo
máximo con la reja limpia y con 50 % de obstrucción

El perfil final del canal de cribado se presenta en la siguiente figura

Bibliografía

Ministerio de Vivienda Ciudad y Territorio (2017). Reglamento de Agua Potable y

Saneamiento Basico. Resolucion 0330 de 2017.

Qasim, S. 1998. Wastewater Treatment Plants-Planning design and operation. 2th Edition.
Technomic Publication

Qasim, S., Zhu, G. 2018. Wastewater Treatment and Reuse, Theory and Design Examples,
Volume 1: Principles and Basic Treatment. 1st Edition. CRC Press

14