FACULTAD DE INGENIERÍA INDUSTRIAL Y MECÁNICA

ESCUELA PROFESIONAL DE INGENIERÍA CIVIL

INFORME DE LABORATORIO: AFORADOR PARSHALL O MOLINETE HIDROMÉTRICO.

ASIGNATURA

HIDRÁULICA DE CANALES

AUTOR

ACUÑA CHICLAYO CRISTIAN

DOCENTE

ING. VICTOR TEPE ATOCHE

CHICLAYO – PERÚ

2019

Ingeniería del Agua. Vol. 2 Num. 3 (septiembre 1995) p. 53

 FÓRMULAS PRÁCTICAS PARA EL DISEÑO DE
CANALES SIN REVESTIR EN TERRENOS
ALUVIALES

Hector Daniel Paria

Instituto de Recursos Hídricos. Facultad de Ciencias Exactas y Tecnologías, Universidad
Nacional de Santiago del Estero (Argentina)

RESUMEN: Se presentan tres conjuntos de fórmulas prácticas, con una base racional, para
el diseño de canales sin revestir excavados en terrenos aluviales finos. Las ecuaciones, del
tipo "régimen", permiten estimar la geometría hidráulica de canales arenosos en equilibrio
dinámico. Los grupos de fórmulas se obtienen mediante la combinación de funciones de
resistencia al flujo, transporte de sedimentos y estabilidad de márgenes. La capacidad
predictiva de las mismas se evalúa empleando una importante base de datos que cubre los
rangos típicos de variación de las variables involucradas en el diseño. Se demuestra que los
errores que se obtienen aplicando estas fórmulas como herramientas de diseño se
encuentran dentro de los márgenes de tolerancia típicas de la dinámica fluvial. Asimismo,
se presenta una serie de ejemplos prácticos que ilustran la aplicabilidad de la metodología
propuesta en el diseño de canales arenosos de envergadura diversa.

INTRODUCCIÓN
para el cálculo de canales arenosos en régimen, así como
La planificación, explotación, administración y control para el diseño de obras de estabilización en cauces
racional de los recursos hídricos de una cuenca es un fluviales.
tópico de gran trascendencia en regiones semi-áridas. En
particular, los temas relativos a la conducción eficiente Las técnicas desarrolladas se verifican rigurosamente a
del agua a través de canales merecen especial través de una contrastación de sus predicciones teóricas
consideración en tierras bajo riego. Este problema con datos de canales estables reales publicados en la
adquiere un nivel importante de complejidad cuando se literatura.
trata el diseño de canales sin revestir excavados en suelos
aluviales. Ello se debe a la gran cantidad de variables a Las metodologías originales desarrolladas, y los
tener en cuenta y a la falta de un conocimiento adecuado promisorios resultados alcanzados a la fecha, marcan una
en lo que respecta a ciertos fenómenos físicos vinculados perspectiva futura muy alentadora en lo que respecta al
con la formación y conservación de los canales a diseñar. mejoramiento de las prácticas actuales de diseño
hidráulico de canales de tierra, y sus múltiples
aplicaciones a la ingeniería de recursos hídricos:
En este trabajo se presenta un conjunto de fórmulas
conducciones para riego y drenaje, estabilización de
prácticas de diseño, basadas en el análisis de los aspectos
cauces fluviales, canalizaciones de ríos, etc.
esenciales de la estabilidad de canales en planicies
aluviales, como lo son: fricción aluvial, transporte de
material sólido y estabilidad de márgenes. El estudio OBJETIVOS DEL ESTUDIO
permite arribar a la obtención de un menú de ecuaciones
de diseño con base racional, que podrían aplicarse tanto
La finalidad del presente estudio es la de obtener un
conjunto de fórmulas prácticas de diseño para canales

Artículo recibido el 17 de julio de 1995 y aceptado para su publicación el 5 de septiembre de 1995. Pueden ser remitidas
discusiones sobre el artículo hasta seis meses después de la publicación del mismo. En el caso de ser aceptadas, las
discusiones serán publicadas conjuntamente con la respuesta de los autores en el primer número de la revista que aparezca
una vez transcurrido el plazo indicado.

Ingeniería del Agua. Vol. 2 Num. 3 (septiembre 1995) p. 54

 DIS EÑ O DE CANALES SIN REVESTIR

arenosos en régimen. Las mismas estarán basadas en un comportamiento se ilustra en la Figura 1.

análisis cuantitativo racional, apoyado en la mecánica de
los procesos físicos dominantes: resistencia al flujo, Considerando lo antedicho, puede decirse que un canal
transporte de sedimentos y estabilidad de las márgenes en aluvial se encuentra en una condición de estabilidad,
un canal aluvial. Para ello se procede a la identificación equilibrio dinámico o régimen, cuando la capacidad de
de las variables que controlan la actividad morfodinámica transportar sedimentos por parte del flujo se balancea con
(erosión y deposición) del flujo sobre los contornos del la tasa de abastecimiento de material sólido al tramo
cauce, se establecen las relaciones funcionales de considerado.
vinculación entre las mismas y, mediante el uso de
relaciones publicadas en la literatura (y de probada Introducido el concepto de régimen, en la fase de proyecto
capacidad predictiva), tanto en sus versiones originales así el problema consiste en dimensional el canal de forma tal
como también modificadas por el autor, se obtienen varios que su configuración sea lo más aproximada posible a la
conjuntos de fórmulas prácticas de diseño. Finalmente, las geometría estable, lo que minimizará los procesos de
relaciones obtenidas se aplican a situaciones típicas de erosión y sedimentación, y consecuentemente los costos
diseño, y las predicciones teóricas obtenidas se comparan de operación y mantenimiento del sistema.
con datos de canales reales que funcionan en condición de
estabilidad dinámica o régimen. Las variables que se consideran relevantes para describir
el flujo en canales aluviales en régimen son: el

CANALES ALUVIALES. VARIABLES Y

RELACIONES FUNDAMENTALES

Si se excava un canal de geometría inicial dada

(por ejemplo, trapecial; ver Figura 1), sin
revestimiento, en una planicie aluvial, en
primera instancia la fracción sólida del
contorno de la sección transversal del canal
estará constituida por un material de naturaleza
semejante al suelo local. Luego de un período
razonablemente prolongado de operación
continua del sistema (generalmente del orden
de magnitud de un año), se genera un proceso
de "maduración" del canal, y el material del
perímetro comienza a presentar características
comparables a las del sedimento transportado
por el flujo, que a su vez es de naturaleza
semejante a la del material que transporta el río Figura 1. Esquema de definición para las distintas situaciones que se
de donde se toma el agua para el canal en pueden presentar en un canal aluvial.
cuestión. Este fenómeno de maduración es una
consecuencia de la actividad morfodinámica

desarrollada por el flujo, a través de los fenómenos de

erosión y sedimentación. Si el canal se proyecta (y se
construye) de manera tal que su geometría tridimensional caudal líquido dominante Q , la carga volumétrica total de
sea capaz de conducir adecuadamente el rango de sedimentos Qs, el tamaño mediano de las partículas de
caudales líquidos de servicio y transportar eficientemente material sólido d , la pendiente de la superficie libre del
las concentraciones de material sólido asociadas a esos flujo S, la densidad del agua p, la viscosidad cinemática
caudales líquidos, el proceso de maduración se del agua v, la densidad del sedimento ps , la aceleración
desarrollará sin producir grandes variaciones en la de la gravedad g , la tensión de corte promediada en la
morfología general del canal (Figura 1.a). sección , la velocidad media del flujo V , características
geométricas de la sección transversal (que pueden ser el
Por su parte, si el diseño inicial no resulta apropiado perímetro mojado P y el radio hidráulico R , o
(Figuras l.b y l.c), lo más probable es que se produzcan alternativamente, el ancho superficial T y la profundidad
erosiones y/o deposiciones de material sedimentario hasta central H ) y, finalmente, dos conjuntos de parámetros
que el flujo logre "esculpir" (modelar) una geometría adimensionales: uno asociado con la forma del canal (i
(transversal y longitudinal) adecuada para conducir los ) y el otro relacionado con la forma de las partículas de
caudales líquidos y sólidos de la manera más eficiente sedimento y su distribución granulométrica (j ). El
posible (Parias, 1990). Este conjunto de parámetros de forma de la sección del canal,

Ingeniería del Agua. Vol. 2 Num. 3 (septiembre 1995) p. 55

 DIS EÑ O DE CANALES SIN REVESTIR

V  a0 d a1 Ra2 S a3 (2)
i, proporciona las relaciones entre las variables geométricas donde a0, a1 , a2 y a3 son constantes numéricas.
seleccionadas, es decir: T = fT(P,R,i) y H = fH(P,R,i) , o

alternativamente P = fp(T,H,i) y R = fR(T,H,i). Esta ecuación, en la instancia de diseño, debe combinarse

con la de continuidad para flujo permanente uniforme:
Como se induce de los párrafos anteriores, puede afirmarse

que, desde un punto de vista macro-morfológico, un canal

QP R V (3)
aluvial es un sistema físico con tres grados de libertad. En

efecto, el mismo puede experimentar ajustes de su geometría En el marco del presente trabajo, todas las ecuaciones de
transversal (ancho y profundidad), así como también de la fricción a emplear fueron reducidas al formato indicado por
longitudinal (pendiente), ante diferentes estímulos (caudal la Ec. 2. Para ello, algunas simplemente se manipularon
líquido y transporte de sedimentos) impuestos natural o algebraicamente, reemplazando valores típicos de algunas
artificialmente. Desde el punto de vista matemático, ello se constantes (por ejemplo, g = 9.806 m/s2 , = 1.65). En otras
traduce en la necesidad de disponer de tres ecuaciones para situaciones, tal el caso de las fórmulas de Brownlie, se
resolver otras tantas incógnitas (ancho, profundidad, adoptó un valor de 1.5 para la desviación estándar
pendiente; o bien, más precisamente: P,R,S: perímetro geométrica de la distribución granulométrica del sedimento
mojado, radio hidráulico, pendiente). Esas ecuaciones fluvial g -. Este es un valor medio para los datos
resultan ser, teniendo en cuenta los procesos físicos disponibles, y, generalmente no induce grandes errores.
intervinientes, las que se detallan en las secciones siguientes.

Con respecto a este ítem, cabe destacar que en un número Ecuación de Fricción de Brownlie:
importante de casos, se usaron versiones simplificadas de las
ecuaciones originales. El propósito fue el de expresar las La fórmula de fricción de Brownlie (1981) para régimen
relaciones pertinentes como funciones monomias (productos fluvial inferior (rizos y dunas), se presenta originalmente de
de potencias) de las variables independientes. la siguiente forma:

Teniendo en cuenta el rango de variación de las variables 

  0.6539
Rs  0.3724  VRS  S 0.09188 0.1050 (4)
más relevantes, los errores inducidos por estas 
 

3
aproximaciones resultaron muy poco significativos, y su d  gd  g

incidencia resulta despreciable en las fórmulas de diseño  

finales. Los detalles atinentes a cada una de las donde g simboliza la desviación estándar geométrica de la
simplificaciones empleadas se presenta en cada una de las distribución granulométrica del sedimento fluvial en
formulaciones correspondientes. cuestión, que se supone representada fielmente por una
distribución log-normal. Introduciendo las simplificaciones
A continuación se presentan las generalidades, así como antes detalladas, y agrupando términos, se obtiene la
también las ecuaciones específicas empleadas, de cada uno siguiente fórmula:
de los tres procesos físicos gobernantes de la geometría
hidráulica de canales aluviales en régimen: Resistencia al V  13.29 d 0.029 R 0.529 S 0.389 (5)
Flujo, Transporte de Sedimentos y Estabilidad de Márgenes.
Ecuación de Fricción de Karim y Kennedy:

Tal como se presenta originalmente (Karim y Kennedy,

Resistencia al Flujo 1990), esta ecuación puede escribirse de la siguiente manera:

Esta relación puede expresarse a través de una relación

funcional para la velocidad media, tal como la siguiente:

q 

0.376
 
 2.822  q  S 0.310 (6)

V  f v (d , R, S; g, ,  s , v, i ,  j ) (1) 
  
gd gd 3 


De manera análoga al caso anterior y, teniendo presente que
El punto y coma ";" indica que todas las variables a su q = V R , luego de manipular algebraicamente la relación
derecha quedan definidas si se especifican las propiedades anterior, se llega a la siguiente forma:
físicas del agua, forma y distribución granulomé-trica de

Ingeniería del Agua. Vol. 2 Num. 3 (septiembre 1995) p. 56

 DIS EÑ O DE CANALES SIN REVESTIR

los sedimentos y el tipo geométrico de sección transversal V  18.19 d 0.103 R 0.603 S 0.497 (7)
adoptada para el canal.
Ecuación de Fricción de Pavlovsky-Lacey:
Haciendo abstracción de las propiedades del fluido,
sedimento, etc. (que para el caso de canales arenosos varían

dentro de un rango muy estrecho), un caso particular La fórmula de resistencia al flujo de Pavlovsky, tal como se
aceptable de la ecuación anterior puede ser una fórmula presenta en los textos clásicos de hidráulica de canales
potencial del tipo Manning: (e.g., Chow, 1959, pp. 100-101), se escribe como sigue:

Ingeniería del Agua. Vol. 2 Num. 3 (septiembre 1995) p. 57

DIS EÑ O DE CANALES SIN REVESTIR

V
1
n

R 0.5  2.5 n  0.13  0.75 R n  0.10  S 0.5
(8)
Ecuación de Transporte de Brownlie (Simplificada):

donde n es el coeficiente de rugosidad de Manning. Esla La ecuación de transporte propuesta por Brownlie (1981) se
ecuación fue aproximada mediante un análisis de regresión basa en técnicas de análisis dimensional, combinadas con un
múltiple a la siguiente: V = 0.8985 n-1.061 R0.736 S0.515 , la cual ajuste experimental empleando análisis de regresión
produce estimaciones dentro de un 5% con relación a la múltiple sobre un importante volumen de datos que incluyen
original de Pavlovsky (Ec. 8) para el rango de datos de la ensayos de laboratorio, así como también mediciones en
base total a emplear en la verificación de las ecuaciones de canales artificiales y ríos naturales. En su versión original, la
diseño. Por su parte, teniendo en cuenta la definición del ecuación de transporte de Brownlie es la siguiente:
factor de sedimento de Lacey, luego de unas pequeñas
manipulaciones, surge la expresión de n en función del

 7115C 
V  R 0.3301
tamaño mediano de los granos, d, tal como sigue: n = [d 1/8 / C  4.5960.5293S0.14050.1606 
1.978
S 0.6601


   
16.7] (Ackers, 1983). Reemplazando esta relación en la SG F
 *c g    
d
 gd 
fórmula aproximada de Pavlovsky, se arriba a la siguiente 
(13)
ecuación:

donde CSG es la concentración de sedimentos en peso

V  17.81 d 0.133 R 0.736 S 0.515 (9) expresada en ppm, CF es una constante que vale 1.268, y T *c
que se designará como fórmula de Pavlovsky-Lacey a los es la tensión crítica de Shields que, según Brownlie, puede
fines del presente estudio. calcularse con la siguiente ecuación:

Como puede apreciarse en los parágrafos precedentes, es  c  0.22Y  0.06 10 7.7Y * (14)
posible expresar prácticamente todas las fórmulas de fricción * *
con – 0.6
en el formato simplificado de la Ec. 2, a la vez que resulta g   1/2 d 3/2
interesante advertir la afinidad de las Ees. 5, 7 y 9 con la Y*  ( v
) (15)
clásica relación de Manning-Stric-kler (Chow, 1959,

p. 206). Luego de aplicar un análisis de regresión múltiple sobre un

banco teórico de datos obtenido a partir de las fórmulas
precedentes, se obtuvo la siguiente ecuación aproximada:
Transporte de Sedimentos

q s  1.535 10 4 d 0.759V 3.774 R 0.632 S 0.758 (16)

Generalmente, las ecuaciones de transporte se presentan en
la forma de una relación funcional para la predicción de la
carga total de material sólido por unidad de ancho del

canal, tal como la siguiente:

Ecuación de Transporte de Karim y Kennedy
(Simplificada):
qs  f qs (d , R, S,V ; g,  ,  s , , i ,  j ) (10)

De manera análoga al caso de fricción, existen Del menú de ecuaciones de transporte presentado por estos
formulaciones simplificadas de qs , que se expresan como autores (Karim y Kennedy, 1990), se consideró la segunda
productos de potencias, o bien como funciones del exceso de ecuación, la que en su versión original se escribe:
algún indicador de la intensidad del flujo sobre un valor
crítico (por ejemplo, tensión de cizalla-miento, velocidad, 3.369
 
 
qs V
 0.00151 
0.840
potencia del flujo, etc.).   1/2 1/ 2
(17)
Asimismo, la ecuación de continuidad para transporte se escribe
como:
En el caso más sencillo (funciones monomias), la ecuación
de transporte de sedimentos puede escribirse como:

q s  b0 d b1V b2 Rb3 S b4 (11)

Ingeniería del Agua. Vol. 2 Num. 3 (septiembre 1995) p. 58
 DIS EÑ O DE CANALES SIN REVESTIR

  transporte de partículas. Nuevamente, y siguiendo un

  * *c
análogogd
  gd 3
procedimiento al del caso anterior, se obtuvo una
donde t* = (R S)/(d ) es la tensión adimensional de función monomia simplificada de la ecuación precedente, la
Shields, y t* es su valor crítico para iniciación del
C
cual se escribe a continuación.
qs  9.74 10 5 d 0.682 V 3.334 R 0.623 S 0.625 (18)
Cs  qs P / Q (12)
Ecuación de Transporte de Levi (Simplificada):

en la que Cs es la concentración volumétrica de sedimentos

y P el perímetro mojado del canal. Las relaciones de La fórmula propuesta por Levi en 1948 (Goncharov, 1962;
transporte sólido consideradas en el presente estudio son las Maza Alvarez, 1987) puede escribirse como sigue:
siguientes.

 d 0.25   
3

qs  0.002 d(V Vc) R   V 

 (19)

   gd 

Ingeniería del Agua. Vol. 2 Num. 3 (septiembre 1995) p. 59

 DIS EÑ O DE CANALES SIN REVESTIR

donde Vc indica el valor de la velocidad media crítica del líquido circulante, p y  son constantes. Por su parte, se ha
flujo que, según Lcvi, se calcula con las siguientes demostrado que, en la condición estable, la inclinación de los
ecuaciones: taludes depende del caudal líquido y del tamaño del
sedimento: z = z(Q,d) (Farias, 1992). A su vez, estudios
Vc   R  de transporte sólido muestran que: Cs = Cs(V), siendo V la
 1.4 11n     si 10  R / D  60 velocidad media del flujo, que, de acuerdo a
  7d 
investigaciones sobre fricción en lechos aluviales se
Vc R expresa como: V = v(d,R,S) (Stevens, 1989). Todo este
 1.4In  7d
  si R /d  60
razonamiento posibilita la obtención de una relación
 
funcional genérica para el fenómeno de estabilidad de
márgenes de un canal aluvial, la cual se puede simbolizar
Teniendo en cuenta estas relaciones, la fórmula aproximada como: P = p(Q,d,R,S). O bien, de una manera más
encontrada para la tasa total de transporte sólido por unidad general:
de ancho del canal es la siguiente:

P   p (Q, d , R, S; g, ,  s , v, i ,  j ) (21)
q s  4.88 10 7 d 0.774V 5.495 R 0.449 (20)

A los efectos prácticos de proyecto de canales en terrenos

aluviales, es necesario asignarle una estructura matemática
Estabilidad de Márgenes
concreta a la relación funcional genérica anteriormente
encontrada. En general, se asume que puede adoptarse un
Las condiciones de estabilidad de las partículas de sedimento
producto de potencias (Chítale, 1966; Ackers, 1983; Farias,
que componen los taludes de un canal aluvial dependen del
1991), del tipo:
valor de la tensión media de corte generada

por el flujo, la cual se expresa como: TO = y R S , en la que P  c o Q c 1 d c2 R c 3 S c4 (22)

 es el peso específico del agua, R es el radio hidráulico del

cauce y S la pendiente longitudinl de la superficie libre del Como puede apreciarse, la clásica ecuación de Lacey para el
flujo. Estudios teóricos y experimentales (Farias, 1990) perímetro mojado, P = 4.84 Q0.5 (unidades SI), es un caso
permiten establecer que la tensión tangencial máxima sobre particular de la ecuación anterior.
los contornos laterales de un canal trapecial está dada por:
Las ecuaciones de estabilidad lateral consideradas en el
tlM = kt  , donde k = 0.75 . El canal sufrirá erosión si ti
0
marco del presente estudio incluyen las siguientes:
supera a la tensión tractiva crítica ( ) que depende a su vez
C

del tamaño mediano de las partículas de sedimento (d) y de

las características físicas del fluido (densidad p y viscosidad Ecuación de Chitale:
). Cuando ocurra erosión de márgenes, el canal aumentará
su ancho de fondo B y se modificará la inclinación (z) de sus Esta relación se obtuvo a través del análisis estadístico de
taludes, aumentando el perímetro mojado P del cauce y unos 350 conjuntos de datos colectados en canales de la
disminuyendo l , hasta que se alcanza la estabilidad. Por India, Pakistán, Egipto y Estados Unidos (presentados por
otra parte, si el canal es muy ancho, la turbulencia del flujo Chítale en 1966), lo cual permitió determinar los siguientes
genera una transferencia lateral de cantidad de movimiento, valores para las constantes que aparecen en la Ec. 22: c = 0

produciendo la deposición de partículas sólidas (las más 6.592 , c , = 0.414 , c2 = 0.1 15 , c , = 0.209 , c4 = -0.097
1 3

finas) sobre el talud hasta modelar la margen estable (unidades SI) (Ackers, 1983). De esta manera, la ecuación
(Pizzuto, 1984). Estos procesos físicos dependen del caudal de Chítale puede escribirse como:
líquido Q

conducido por el canal y de la concentración total de P  6.592 Q 0.414 d 0.115 R 0.209 S 0.097 (23)
sedimentos Cs transportada, la que a su vez es función de la
velocidad media V del escurrimiento (Ste-vens, 1989).
Ecuación de Mao y Flook:
El problema consiste ahora en determinar de que manera se
vinculan entre sí las variables antes mencionadas. Por La expresión matemática de esta relación fue obtenida por el
definición, el canal habrá adquirido una configuración autor mediante un ajuste del gráfico presentado por Mao y
estable cuando su geometría tridimensional (P,R,S) Flook (1971), en el que se indica la relación entre P/R y el
permanezca estadísticamente inalterable en el tiempo bajo la caudal Q, observada para los canales denominados CHOP
acción del flujo de agua y transporte de sedimentos. El radio (Canal and Headworks Observation Project), en actual
hidráulico R, o la profundidad H del flujo, están controlados territorio de Pakistán. La relación obtenida para la curva es:
por la fricción en el contorno; el gradiente S P/R = 7.420 Q0.2287, la que puede escribirse como:
Ingeniería del Agua. Vol. 2 Num. 3 (septiembre 1995) p. 60
 DIS EÑ O DE CANALES SIN REVESTIR

depende de la capacidad de transporte del flujo. Finalmente, P  7.42 Q 0.229 R (24)

P está gobernado por la actividad morfodinámica del flujo
sobre los taludes, la cual, en la

condición de estabilidad, deberá ser estadísticamente nula. De esta manera, se obtiene una expresión particular de la la
Esta condición se expresa por una función de la forma: P = Ec. 22, con c2 = c4 = 0.
Øp(Q,Cs,d,z,p,u). Si se conocen las propiedades físicas del

Ingeniería del Agua. Vol. 2 Num. 3 (septiembre 1995) p. 61

DIS EÑ O DE CANALES SIN REVESTIR

Tabla 1. Fórmulas obtenidas para el diseño de canales aluviales en régimen.

Código Proceso Autor Ecuaciones Usadas Ecuaciones de Régimen

Fricción Brownlie V=13.29d-0029 R0.529 S0.389 P = 2.70 d0.033 Cs-0.080 Q0.514

BBC Transporte Brownlie qs=1.5351000 d-3.774R0.632 S0.758 R = 0.049 d-0.140 Cs-0.102 Q0.375

Estabilidad Lat Chítale P =6.592 d0.115 Q0.414 R0.209 S-0.097 S = 14.645 d0.541 Cs0.605 Q-0.224

Fricción Karim-Kennedy V=18.19 d-0.103 R0.603 S0.497 P = 0.664 d-0.066 Cs-0.120 Q0.562

KKM Transporte Karim-Kennedy Qs =9.74 10-5 d-0.682 V3.334 R0.623 S0.625 R = 0.089 d-0.066 Cs-0.120 Q0.333

Estabilidad Lat. Mao-Flook (Farias) P= 7.42 Q0.229 R S = 16.04 d0.555 Cs0.630 Q-0.192

Fricción Pavlovsky-Lacey V = 17.81 d-0.133 R0.736 S0.515 P = 0.710 d-0.089 Cs-0.115 Q0.519

PLG Transporte Levi (Modif.) qs = 4.8810-7 d-0.774 V5.495 R-0.449 R = 0.112 d-0.063 Cs-0.081 Q-0.363

Estabilidad Lat. Glushkov (Modif.) P = 16.12 R1.428 S = 11.498 d0.643 Cs0.497 Q-0.292

respectivas, las relaciones de fricción, transporte y

Ecuación de Glushkov: estabilidad lateral (reducidas al formato dado por las Ees.
2, 11 y 22, respectivamente), así como también las
Esta fórmula, muy empleada por autores rusos, en su ecuaciones racionales de régimen obtenidas en un formato
versión original (aparentemente publicada en 1925) generalizado, tal como el siguiente:
expresa una relación entre el ancho del cauce (no se

especifica cual) y la profundidad media del flujo, tal

P  x0d x 1
C xS Q x
2 3
(26)
como: Bm = KG H, donde m y KG son constantes y y2 y3
numéricas que dependen del tipo de material de lecho y
R  y0 d 1 CS Q (27)
márgenes presentes en el cauce en cuestión (Maza
Alvarez, 1987). Puesto que, en gran medida, la fórmula S  z0 d C S z1 Q z2
(28)
z3

de Glushkov fue calibrada con datos de ríos naturales, en donde las constantes xj , yj , zj (j = 0,1,2,3) resultan de la
los que se cumple que P = B y R = H , con un error mínimo, combinación algebraica de las constantes que aparecen en
puede obtenerse una ecuación particular de la Ec. 22, las funciones fundamentales (Ees. 2, 11, 22), luego de
usando los valores adecuados de m y KG. Para cauces con aplicar las ecuaciones de continuidad en el flujo líquido y
lecho arenoso y márgenes cohesivas (caso típico de canales transporte sólido. En el Apéndice 1 se expone el conjunto
en régimen) las constantes adoptan los siguientes valores de transformaciones algebraicas necesarias para la
medios: m = 0.7, KG = 7. Así, se obtiene: obtención de las ecuaciones de diseño (Ees. 26, 27 y 28) a
partir de las relaciones fundamentales, es decir, fricción
P = 16.12 R1.428 (25) aluvial, transporte sólido y estabilidad lateral (Ees. 2, 11 y
22).
Nuevamente, esta expresión es un caso particular de la Ec.
22, con c = c2 = c4 = 0.
1
Para el proyecto completo de un canal estas ecuaciones
deben complementarse con las siguientes, obtenidas de
relaciones geométricas elementales para canales de sección
FÓRMULAS DE DISEÑO transversal trapecial:

Uno de los objetivos principales del presente estudio es el A  PR (29)

de presentar ecuaciones de diseño para canales arenosos en k z  2 1 z 2 (30)
régimen que sirvan como herramientas prácticas, pero con
una base conceptual sólida. Por ello, se adoptaron
P  P 2  4(k z  z)A
relaciones de fricción y transporte que a la fecha parecen H  (31)
tener un nivel razonable de aceptación, debido a su 2(k z  z)
calibración con gran cantidad de datos y a su A  zH 2
 (32)
comportamiento aceptable cuando se comparan sus H
predicciones con valores observados en canales existentes

Ingeniería del Agua. Vol. 2 Num. 3 (septiembre 1995) p. 62

 DIS EÑ O DE CANALES SIN REVESTIR

(Chítale, 1966; Smith, 1971; Brownlie, 1981; Bakker y En relación a la pendiente transversal de los taludes (z), la
Vermaas, 1986). En la Tabla 1 se presentan, codificadas de misma puede estimarse teniendo en cuenta las característi-
acuerdo a los autores de las funciones fundamentales

Ingeniería del Agua. Vol. 2 Num. 3 (septiembre 1995) p. 63

DIS EÑ O DE CANALES SIN REVESTIR

cas físicas y mecánicas del suelo en el que se excava el

canal (Chow, 1959). Sin embargo, en un estudio sobre
la forma de la sección transversal en canales arenosos en
régimen (Parias, 1993a,b,c), se ha encontra- do una
ecuación que permite estimar z en función de parámetros
adimensionales característicos del flujo y del sedimento,
a saber:
z = 0.022 d 0.368 Q 0.161 (33)
* *

1/2 1/2 5/2

donde Q =
* Q/(g D d ) y d =* [(g D / n2)1/3 d] .

EVALUACIÓN DE LAS FÓRMULAS DE DISEÑO

Con el objeto de evaluar el comportamiento de las

ecuaciones de régimen obtenidas como herramientas de
diseño de canales aluviales, se efectuó una comparación
de sus predicciones con datos de campo. Para ese fin se
emplearon las compilaciones de datos presentadas por
Simons y Albertson (1963), Chítale (1966), Smith (1971),
Brownlic (1981) y Bakker y Vermaas (1986). En la Tabla
2 se presenta un resumen de las compilaciones de datos
empleadas, reducidas todas al formato [N,Q,d,P,R,S,Cs],
donde N indica la cantidad de datos contenida en cada
conjunto y las restantes variables (cuyos rangos de
variación se indican el la tabla) han sido definidas
previamente. Los datos fueron extraídos de las
publicaciones citadas, y sometidos a un análisis y
depuración preliminar, que consistió básicamente de los
procedimientos que a continuación se describen.

Algunas compilaciones (fundamentalmente las más

antiguas) presentan la información en unidades inglesas.
En estos casos, los datos dimensionales de interés se
redujeron al sistema métrico SI. Luego se llevó a cabo una
selección, ya que muchas compilaciones presentan
registros múltiples de mediciones practicadas sobre el
mismo canal. En estos casos, se adoptó el criterio de
seleccionar el registro de mayor caudal y descartar los
restantes. Asimismo, se practicó un análisis de
consistencia, verificándose el cumplimiento de la
ecuación de continuidad: Q = P R V, y descartando los
registros que no satisfacían esta condición
(probablemente por errores tipográficos o de
transcripción), con una tolerancia del 2%. Algunos
autores no presentan información explícita del perímetro
mojado P y/o el radio hidráulico R, pero sí del ancho de
fondo B (o de superficie T) y de la profundidad central H,
conjuntamente (en ciertos casos) con información relativa
a la inclinación lateral de los taludes (z). En esos casos, P
y R se determinaron mediante relaciones geométricas para
canales trapeciales. En los casos en los que no se disponía
de información atinente a los taludes, y considerando un
canal prismático rectilíneo de forma trapecial, el criterio
que se propone en el presente estudio es el de definir la
pendiente lateral de una sección trapecial hipotética, de Figura 2. Gráfico del Perímetro Mojado observado versas el
manera tal que el calculado, para los tres grupos de fórmulas: (a) BBC
[o]; (b) KKM [o]; (c) PLG [D].

Ingeniería del Agua. Vol. 2 Num. 3 (septiembre 1995) p. 64

 DIS EÑ O DE CANALES SIN REVESTIR

Tabla 2. Resumen de los bancos de datos usados para la verificación de las fórmulas de diserto

Fuente de N Q d P R S Cs
Datos. - M3/s mm m m X104 ppm
PRC-ACOP 292 0.07-138.03 0.064-0.357 1.70-58.40 0.18-2.865 0.60-4.73 16-727
Sind 44 1.29-256.47 0.020-0.196 6.00-87.50 0.594-3.572 0.54-1.446 39-763
Simons & Bender 24 1.22-29.42 0.029-0.805 3.90-23.60 0.497-1.923 0.58-3.88 10-526
CHOP 17 27.61-407.76 0.146-0.261 23.80-111.80 1.759-3.718 0.746-2.0 51-531
Jhelum 17 0.22-39.84 0.171-0.353 2.50-30.80 0.271-1.682 1.30-3.40 36-100
Chenab 23 0.64-254.86 0.217-0.425 3.50-89.30 0.39-2.89 1.60-4.00 46-179
UttarPradesh 24 0.42-282.15 0.075-0.315 2.20-59.60 0.274-3.706 1.02-4.34 83-908
Smith 57 0.16-262.72 0.135-0.425 2.0-89.30 0.241-3.679 1.07-4.13 59-861
Total 496 0.070- 407.76 0.020-0.805 1.70-111.80 0.180-3.718 0.54-4.73 10-908

término generales, mientras que para el caso de la

canal resultante sea de igual área transversal que la
pendiente, variable fuertemente afectada por el transporte
sección natural estable (en la mayoría de los casos es
de sedimentos, la dispersión observada resulta muy
curvilínea). Mediante el análisis de datos de canales
significativa.
arenosos estables (Simons y Albertson, 1963; Garde y
Ranga Raju, 1985), se obtuvo la siguiente ecuación para
el cálculo de la pendiente lateral: z = 0.022 d, 0368 Q,0161, Para efectuar una evaluación cuantitativa de la bondad de
con Q. = [Q/(g d5)1/2] y dj = [(g /v2)1/3 d]. La misma se las ecuaciones de la Tabla 1 como predictores de la
aplicó en los casos en que no se disponía del valor de la geometría de régimen, se emplearon las razones de
inclinación de taludes z. discrepancia P = PC/P0 , R = RC/R0 , S = SC/S0, en las
que el subíndice "C" indica valor calculado y el "O"
observado. Se determinaron los parámetros estadísticos
Cabe destacar que ciertos conjuntos de datos carecen de
de las discrepancias para cada uno de los conjuntos de
información sobre transporte de sedimentos (Canales UP
ecuaciones y para cada banco de datos analizado.
y CHOP presentados por Chítale, datos de Smith, y
algunos de Simons y Bender). En esos casos, Cs se
determinó como una media aritmética simple de las Con relación al comportamiento de las razones de
predicciones obtenidas con seis relaciones de transporte discrepancia, se pudo observar que las mismas se ajustan
de uso extendido y bondad probada, a saber: Engclund- a una distribución log-normal. En efecto, P y R exhiben
Hansen, Ackers-White, Brownlie, van Rijn, Peterson- un ajuste casi perfecto, mientras que S presenta un mayor
Peterson y Karim-Kennedy. En todos esos casos se aplicó grado de dispersión, pero sin desvirtuar la hipótesis de
la ecuación correspondiente empleando el resto de los log-normalidad. Según Brownlie (1981), una variable x
datos hidráulicos del canal en cuestión. Este aspecto, es log-normalmente distribuida puede describirse a través de
decir, la generación artificial de datos de transporte sólido su media geométrica (xg) y su desviación estándar
es un punto discutible. No obstante, se encontró que para geométrica (σgx). Además, la mediana debe coincidir con
el presente estudio el mismo no genera sesgos de la media geométrica, mientras que el percentil 84% (x84)
importancia en la evaluación de las ecuaciones de diseño.
puede determinarse mediante el producto de (xg) y (σgx)
, y_el percentil 16% (xl6) a través del cociente de (xg)
En las Figuras 2, 3 y 4 se presentan gráficamente los dividido por (σgdx). Todos estos patrones de
resultados obtenidos con cada grupo de ecuaciones para comportamiento fueron verificados para los tres
la predicción de P, R y S, respectivamente, en la forma conjuntos de ecuaciones cuyas discrepancias se analizan,
convencional de gráficos de valores observados versus y para todo el banco de datos empleado. Los parámetros
valores calculados. Como puede advertirse, de la simple estadísticos del análisis de las discrepancias se presentan
inspección visual de estos gráficos, se desprende la en la Tabla 3. En ella, para cada variable a estimar (P,R,S)
evidencia que los predictores de la geometría transversal y para cada grupo de fórmulas de diseño (BBC, KKM,
(P y R) presentan un comportamiento muy aceptable en PLG) se presentan la media y desviación estándar de las
discrepancias, así como también los percentiles más
relevantes.

Ingeniería del Agua. Vol. 2 Num. 3 (septiembre 1995) p. 65

DIS EÑ O DE CANALES SIN REVESTIR

Figura 4. Gráfico de la Pendiente Longitudinal

observada versus la calculada, para los tres grupos
Figura 3. Gráfico del Radio Hidráulico observado versus el de fórmulas: (a) BBC [ o ]; (b) KKM [ o ]; (c) PLG
calculado, para los tres grupos de fórmulas: (a) BBC [D ]
[ o ]; (b) KKM [ o ]; (c) PLG [ D ]

Ingeniería del Agua. Vol. 2 Num. 3 (septiembre 1995) p. 66

 DIS EÑ O DE CANALES SIN REVESTIR

Tabla 3. Resultados del análisis estadístico de las discrepancias para los tres grupos de fórmulas de diseño

Variable Grupo de Parámetros Estadísticos

del Canal Fórmulas  σ8 min 16 50 84 max

BBC 0.9436 0.1493 0.5705 0.7986 0.9315 1.0884 1.8583

P KKM 0.9213 0.1706 0.5161 0.7539 0.9079 1.0885 1.6505
PLG 1.0425 0.1653 0.6418 0.8805 1.0278 1.2041 1.9081

BBC 1.0259 0.0991 0.7752 0.9277 1.0191 1.1239 1.4559

R KKM 1.05840 0.1304 0.7447 0.9301 1.0496 1.1864 1.6429
PLG 0.9401 0.0899 0.7141 0.8506 0.9381 1.0291 1.3432

BBC 1.2409 0.4915 0.4103 0.8312 1.1322 1.6377 4.1119

S KKM 1.0132 0.4218 0.3259 0.6605 0.9225 1.3548 3.4365
PLG 0.9801 0.3538 0.3401 0.6751 0.9312 1.2787 3.3591

• En prácticamente todos los casos, las ecuaciones

ANÁLISIS DE RESULTADOS racionales de régimen presentadas en este estudio
producen mejores predicciones que las resultantes de
Se pudo observar que, en general, ningún grupo de los métodos empíricos clásicos. Sin embargo, debe
ecuaciones tiene supremacía apreciable sobre las otras destacarse la bondad de las ecuaciones originales de
para todo el rango de datos. Sin embargo, se advierten Lacey, que, habiendo sido calibradas con muy pocos
algunos resultados interesantes. datos, aún hoy producen resultados aceptables a nivel
de pre-diseño (de la sección transversal, pero no para
Resulta sorprendente el buen comportamiento predicti- la pendiente).
vo evidenciado por el grupo de ecuaciones simbolizado
como PLG, en el que emplean versiones simplificadas de • Para los grupos de fórmulas analizados, el grupo
ecuaciones de fricción y transporte (Pavlovsky y Lcvi) BBC presenta las mejores predicciones para el radio
prácticamente en desuso en la actualidad, mientras que hidráulico R, lo que confirma las bondades del
para la estabilidad lateral se emplea una versión (también predictor de Brownlic para fricción en canales
modificada respecto de la original) de la ecuación de arenosos.
Glushkov. A la luz de estos resultados, y teniendo en
cuenta que el conjunto de ecuaciones PLG representa en • La ecuación de grupo PLG se comporta
cierto modo una suerte de síntesis de la escuela soviética satisfactoriamente para predecir el ancho en canales
de dinámica fluvial, puede afirmarse el alto grado de grandes (como por ejemplo los canales de trasvase
afinidad (hasta ahora un tanto oculto) entre esta y los CHOP, en Pakistán), mientras que para los canales
conceptos clásicos de la teoría del régimen. menores, el comportamiento de los grupos BBC y
KKM es comparable al PLG.
En la Tabla 4 se presenta el comportamiento de cada
grupo de fórmulas expresado como porcentaje de sus • En general, todos los grupos de ecuaciones exhiben
predicciones contenidas entre bandas características de un considerable grado de dispersión en la predicción
errores, las que se escogieron (subjetivamente) como 20% de la pendiente, observándose una cierta ventaja en
para el perímetro mojado, 10% para el radio hidráulico y lo que respecta al comportamiento del grupo PLG,
25% para la pendiente. con relación a los restantes.

Por otra parte, teniendo en cuenta que el radio hidráulico • Este último punto demuestra una vez más que el
está asociado a la ecuación de fricción empleada, y ramayor nivel de incertidumbre de la dinámica
análogamente la pendiente con el predictor de transporte fluvial sigue presentándose en la predicción del
y el perímetro mojado con la relación empleada para transporte de sedimentos.
cuantificar la estabilidad lateral, pueden enunciarse los
siguientes patrones de comportamiento:

Ingeniería del Agua. Vol. 2 Num. 3 (septiembre 1995) p. 67

DIS EÑ O DE CANALES SIN REVESTIR

Tabla 4. Errores de estimación para los tres grupos de fórmulas gráficamente los patrones generales de comportamiento en
de diseño lo que se refiere a la predicción de la geometría transversal.
En los esquemas, sobre un sistema de coordenadas (x,y), se
dibujan las secciones observadas y se superponen las
calculadas, como para visualizar los eventuales procesos
Estimación de Estimación Estimación morfológicos que se desarrollarían como consecuencia de
P de de los diseños teóricos.
R S
Predictores
% de datos % de datos % de datos Ejemplo 1. Diseño de un gran canal de riego en Egipto: En
entre entre entre este ejemplo se trabaja con los datos del Canal Menufia
0.80 y 1.25 0.90 y 1.11 0.75 y 1.33 (Egipto), presentados por Abou-Seida y Saleh (1987, p.
444): Q = 245 m3/s , d = 0.26 mm , Cs = 50 ppm , z = 1.
BBC 80 71 68
Cabe recordar que las ecuaciones de la Tabla 1, d se expresa
KKM 71 59 65 en metros y Cs se expresa en forma volumétrica
PLG 85 62 68 adimensional, por lo que si se tiene el dato en parte por
millón en peso (como se presentan habitualmentc), se debe
EJEMPLOS DE CASOS DE DISEÑO reducir el mismo a volumen. En este caso, y suponiendo un
valor de 1.65 para la densidad relativa del sedimento
Con la finalidad de ilustrar la manera en que se aplican en sumergido, se tiene: Cs = 50/2.65xl06 = 1 .886810-5 . Hecha
la práctica las fórmulas de diseño de canales de tierra (Tabla esta salvedad, se aplican las ecuaciones en forma directa.
1), se presenta a continuación una serie de ejemplos de Por ejemplo, para el grupo BBC:
diseño en los que se trabaja con datos extraídos de la
literatura y también recolectados por el autor en el marco de P = 2.70·(0.00026)0.0333(1.8868·l0-5)-0080·(245)0514 = 82.99 m
la presente investigación. Con el fin de ilustrar los ejemplos,
los resultados se presentan en tablas (en las que se indican R = 0.049·(0.00026)-140 ·(1.8868 ·10-5)0.605 ·(245)0.375 = 3.71
los errores de predicción en cada caso) y en figuras, en las
m
que se muestran
S = 14.645·(0.00026)0.541·(1.8868·10-5)0.605·(245)-0.224 = 6.81·10-5

Tabla 5. Resumen de resultados obtenidos en la predicción de la geometría estable empleando Luego, empleando las
los tres grupos de fórmulas de diseño. Canal Menufia (Egipto) ecuaciones 29 a 33 se completa
el diseño. Así, se obtiene:
Ecuac P p R R S s B H T
ión M % M % s10 4
% m m m B= 60.72 m; T= 81.23 m; H=
4.34 m
Canal Menufia (Egipto)
BBC 82.99 12.78 3.71 -12.29 0.681 -1.16 60.72 4.34 81.23 Las dimensiones reales del
KKM 92.98 25.37 3.53 -16.41 0.603 -12.43 72.39 4.01 91.35 canal (Abou— . Seida y Saleh,
PLG 89.86 22.12 3.35 -20.8 0.513 -25.54 70.41 3.79 88.32 1987) son:

P = 73.58 m, R = 4.23 m, S
= 6.89 10-5, B = 60.0 m, H =
4.80 m.
para el Canal Menufia (Egipto)

Figura 5. Sección transversal observada y secciones calculadas

Ingeniería del Agua. Vol. 2 Num. 3 (septiembre 1995) p. 68
 DIS EÑ O DE CANALES SIN REVESTIR

Con los grupos de fórmulas restantes (KKM y

PLG) se procede de manera análoga. En la Figura 5 se
presentan gráficamente los resultados obtenidos con
cada grupo para la sección transversal, encontrándose
que, en este caso, el grupo BBC produce resultados más
aproximados. Todos los grupos producen, para este
caso, secciones más anchas y menos profundas que la
sección real del canal. En la Tabla 5 se presenta un
resumen del comportamiento de los grupos de fórmulas
para el presente ejemplo. En la misma se incluyen los
errores de estimación de cada una de las variables
geométricas, los que se calcularon con las fórmulas:

p  Pc  Po  / Po ,R  (R c  R o ) / R o , s  (Sc  So ) / So

donde el sub-índice "c" indica valor calculado y "o"
representa al observado.

Ingeniería del Agua. Vol. 2 Num. 3 (septiembre 1995) p. 69

DIS EÑ O DE CANALES SIN REVESTIR

Ejemplo 2. Diseño de un gran canal de tmsvasamiento en

Pakistán: En este caso, se trata del Canal Nara, cuyos
datos, reportados por Rickard (1989) incluyen un caudal
líquido a conducir de 385 m3/s , un tamaño mediano de
partículas de sedimento de 0.15 mm , con una
concentración total estimada de 150 ppm. Según indica
Rickard, el canal fue originalmente construido con una
sección trapecial, con ancho de fondo de 105.5 m,
profundidad de 3.50 m y taludes laterales 1:1, con lo que
se tiene: P = 115.4 m y R = 3.31 m. La pendiente
longitudinal aproximada es 7.6910-5. Rickard no reporta
las dimensiones estables del canal, ya que los tramos
estudiados del mismo exhiben procesos morfológicos que Figura 6. Sección transversal observada y secciones
denuncian que aún no se ha obtenido una condición de calculadas para el Canal Nara (Pakistán)
régimen. El objetivo del estudio de Rickard fue conseguir
un diseño de régimen para la remodelación del canal Ejemplo 3. Diseño de tres canales de riego en la región
actual. En este estudio, se emplearon los datos del Canal Norte de Argentina: Para este ejemplo se emplearon datos
Nara con la finalidad de ilustrar los resultados que recolectados en el área de riego del denominado "Proyecto
producen las fórmulas de diseño cuando se las aplica a Río Dulce" (PRD), en la provincia de Santiago del Estero,
canales de grandes dimensiones. Los resultados del Argentina. La zona está enclavada en la paleo
planicie aluvial del Río Dulce, y
Tabla 6. Resumen de resultados obtenidos en la predicción de la geometría estable por tanto presenta suelos
empleando los tres grupos de fórmulas de diseño. Canal Nara (Pakistán) superficiales de textura arenosa y
limosa, resultante de depósitos
P p R R S s B H T fluviales. En este estudio se
Ecuación
M % M % s104 % m m m seleccionaron tres canales
Canal Nara (Pakistán) secundarios, que presentan
características de comportamiento
BBC 94.16 - 4.25 28.55 0.888 15.47 66.13 5.05 92.28 cn-marcadas claramente dentro del
18.14
concepto de régimen. Los mismos
KKM 108.95 -5.59 3.73 13.01 0.814 -5.91 85.53 4.22 107.58 poseen trazado cuasi-
PLG 105.16 -8.87 3.74 13.13 0.545 -29.13 81.61 4.24 103.58 rectilíneo, secciones .
transversales razonablemente
cálculo se presentan en la Tabla 6 y en la Figura 6. Puede uniformes, pendientes suaves, lechos arenosos, márgenes
observarse que, en este caso, las fórmulas aquí propuestas estables compuestas por materiales más finos que los del
predicen secciones más estrechas y más profundas que la fondo del canal, transporte de sedimento moderado y están
sec- alimentados por caudales aproximadamente constantes
durante el período de operación (temporada de riego).
ción con la que se construyó inicialmentc el Canal Nara.
Esto marca una buena aproximación a la luz de lo relatado Las operaciones llevadas a cabo para la recolección de
por Rickard (1989), quien indica que, en su búsqueda de información de campo incluyeron: relevamientos
una condición de régimen, el canal ha producido una topográficos (perfiles transversales y longitudinales de
degradación importante del lecho en los tramos lecho y pelo de agua), aforos líquidos a nivel del umbral
superiores. de desborde (bankfull) de los canales (para estimar el
caudal dominante),
Tabla 7. Resumen de datos geométricos, hidráulicos y sedimentológicos de los canales Suri Pozo,
Los Romano y Sud (Argentina) muestreo de material de
Q P R S Cs d16 d 35 d50 d65 d84 d90 lecho y sedimento en
3 4
suspensión. La carga total de
Canal M /s m m X10 ppm mm Mm mm mm mm mm sedimentos se estimó
Suri 9.94 14.51 1.25 1.89 73.21 0.14 0.19 0.22 0.3 0.42 0.46 promediando los resultados
Pozo obtenidos " mediante la
Los 6.66 12.65 0.95 2.86 94.28 0.17 0.26 0.31 0.4 0.49 0.65 aplicación de las seis
Romano funcioncs de transporte antes
indicadas. En la Tabla 7 se
Sud 5.01 9.62 0.75 4.72 256.58 0.15 0.25 0.3 0.35 0.43 0.47 presenta un resumen de la
información relativa a los
tres canales estudiados.

Ingeniería del Agua. Vol. 2 Num. 3 (septiembre 1995) p. 70

 DIS EÑ O DE CANALES SIN REVESTIR

Tabla 8. Resumen de resultados obtenidos en la predicción de la geometría estable empleando menores que 100 ppm, este
los tres grupos de fórmulas de diseño. Canales Suri Pozo, Los Romano y Sud criterio resulta admisible.
(Argentina) Se estima que este
procedimiento podría
Ecuación P  R R S s B H T generalizarse para aquellos
m % m % X10 4
% m m m casos en los que no se
disponga de una estimación
Canal Suri Pozo
previa de la concentración
BBC 15.05 3.72 1.08 -1.36 1.62 -14.28 10.51 1.32 14.21 de
KKM 14.45 -0.42 1.16 -7.45 1.3 -30.94 9.41 1.46 13.52 sedimentos, que es la regla
general en este tipo de
PLG 16.15 11.3 1.01 -19.2 1.44 -23.81 12.05 1.19 15.39
problemas.
Canal Los Romano
BBC 12.43 -1.74 0.88 -7.34 2.46 -13.98 8.94 1.06 11.72
CONCLUSIONES
KKM 11.23 -11.22 0.97 2.61 1.98 -30.69 7.09 1.25 10.38
PLG 12.66 0.08 0.85 -10.53 2.25 -21.33 9.33 1.01 11.98
Se han presentado tres
Canal Sud conjuntos de fórmulas
BBC 9.91 3.01 0.72 -3.99 4.73 0.21 7.13 0.86 9.31 prácticas, con una base
racional, para el diseño de
KKM 8.49 -11.59 0.79 5.09 3.86 -18.17 5.18 1.03 7.78
canales sin revestir
PLG 9/76 1.45 0.71 -5.33 3.95 -16.31 7.02 0.85 9.16 excavados’en terrenos
aluviales finos. Los tres
grupos de ecuaciones
obtenidos, permiten predecir la geometría hidráulica
(ancho, profundidad, pendiente) de canales arenosos que
se encuentren en una condición de estabilidad, equilibrio
dinámico o régimen. Las relaciones se basan en una
combinación algebraica de funciones descriptivas de los
procesos físicos gobernantes de la dinámica de canales
aluviales en régimen, es decir, resistencia al flujo,
transporte de sedimentos y estabilidad de márgenes.

Las fórmulas de diseño obtenidas se han sometido a una

Figura 7. Sección transversal observada y secciones verificación exhaustiva de sus capacidades predic-livas,
calculadas para el Canal Suri Pozo (Argentina)
empleando una importante base de dalos que cubre los
rangos típicos de variación de las variables involucradas
Los resultados de la aplicación de las fórmulas de diseño en el diseño. Se ha demostrado, a través de un análisis
se especifican en la Tabla 8 y en las Figuras 7, 8 y 9. Como estadístico, que los errores que se pueden obtener si se
puede advertirse se obtienen predicciones razonables, emplean estas fórmulas como herramientas de diseño se
dentro de los márgenes de error habitual- mente admitidos encuentran holgadamente dentro de los márgenes de
en dinámica lluvial. tolerancia típicas de la dinámica fluvial. Asimismo, se ha
presentado una serie de ejemplos que ilustran la sencillez
del enfoque propuesto en la instancia de diseño tanto de
Con relación a los ejemplos presentados, debe destacarse pequeños canales de riego, como de grandes
que Cs es un dato que no se dispone en la instancia de conducciones de trasvasamiento de flujos.
diseño. No obstante, se determinó que si se practica un
pre-diseño empleando las ecuaciones clásicas de Laccy
(Farias, 1994), y luego se calcula Cs usando los datos así Los resultados obtenidos demuestran la aplicabilidad de
generados, los resultados que se obtienen con la la metodología propuesta a problemas de cálculo
aplicación de las ecuaciones de régimen son hidráulico de canales aluviales estables y que podrían
prácticamente los mismos que los aquí discutidos para las extenderse al diseño de obras de estabilización de cauces
dimensiones transversales (P y R). Para la pendiente se fluviales, que son materias de gran importancia en la
obtienen discrepancias importantes cuando la planificación, gestión y conservación de recursos hídricos
concentración es alta. Para valores moderados, digamos en regiones áridas y semi-áridas.

Ingeniería del Agua. Vol. 2 Num. 3 (septiembre 1995) p. 71

 DIS EÑ O DE CANALES SIN REVESTIR

Figura 8. Sección transversal observada y secciones Figura 9. Sección transversal observada y secciones
calculadas para el Canal Los Romano (Argentina) calculadas para el Canal Sud (Argentina)

i Subíndice genérico
RECONOCIMIENTO j Subíndice genérico
KG Coeficiente en la ecuación estabilidad lateral de
Parte de los estudios relatados en el presente trabajo Glushkov
fueron desarrollados con financiamiento del Consejo de ki Coeficiente de reducción de la tensión de corte
Investigaciones Científicas y Tecnológicas de la sobre los taludes de un canal
Universidad Nacional de Santiago del Estero (CICYT- m Exponente en la ecuación de Glushkov
UN-SE), bajo el Programa de Recursos Hídricos y N Cantidad de datos contenida en una base de
Ciencias de la Tierra. Asimismo, el autor desea expresar información
su agradecimiento al Consejo Nacional de Investigaciones n Coeficiente de rugosidad de Manning
Científicas y Técnicas (CONICET) de Argentina, por el P Perímetro mojado de la sección de un canal
apoyo brindado mediante una Beca Post-Doctoral. Los Q Caudal líquido
asistentes de Investigación María T. Pilan, Carlos A. Q* Caudal líquido adimensional : Q,= Q/(g1/2 1/2
Roldan y Mariano Macchiarola, brindaron invalorable d5/2)
colaboración en el procesamiento electrónico de datos. Qs Caudal sólido (volumétrico)
q Caudal líquido unitario (por unidad de ancho)
qs Caudal sólido unitario (por unidad de ancho)
LISTA DE SÍMBOLOS R Radio hidráulico
S Pendiente longitudinal
A Área de la sección transversal de un canal T Ancho del canal a nivel de la superficie libre
a¡ Constantes en la ecuación genérica de U Velocidad de cizallamicnto
resistencia al flujo (i = 0,1,2,3) V Velocidad media del flujo en la sección
B Ancho de fondo de un canal o cauce transversal de un canal
bj Constantes en la ecuación genérica de Vc Velocidad crítica para el inicio de transporte en
transporte de sedimentos (j = 0,1,2,3,4) la ecuación de Levi
C Coeficiente de Chézy xi Constantes en la fórmula de diseño del
Cs Concentración total de sedimentos perímetro mojado (i =0,1,2,3)
(volumétrica) Y* Parámetro adimensional de partícula de
CsG Concentración total de sedimentos Brownlie
(gravimétrica) yi Constantes en la fórmula de diseño del radio
cF Factor en la fórmula de transporte de Brownlie hidráulico (i = 0,1,2,3)
ck Constantes en la ecuación genérica de z Pendiente transversal de los taludes (z
estabilidad lateral (k = 0,1,2,3,4) horizontal : 1 vertical) en una sección trapecial
d Tamaño mediano de una partícula de sedimento zi Constantes en la fórmula de diseño de la
dp Diámetro de sedimento tal que el p% de la pendiente (i = 0,1,2,3)
muestra sea igual o más lina que él.  Razón de aspecto de la sección transversal de
d* Diámetro sedimentológico adimensional: d* = un canal
[(g A / n2)10 d] i Parámetro adimensional asociado con la forma
f Función del canal
g Aceleración de la gravedad j Parámetro adimensional relacionado con la
H Profundidad del flujo en el centro del canal forma y granulometría del sedimento

Ingeniería del Agua. Vol. 2 Num. 3 (septiembre 1995) p. 67

 DIS EÑ O DE CANALES SIN REVESTIR

 Peso específico del agua S/13-S/22.

s Peso específico del sedimento Farias,H.D. (1992). "Máximum Efficiency Sandy Channels
 Densidad relativa del sedimento sumergido and theirMorphology", ISRS-5: Sediment Management
X Razón de discrepancia de la variable X (Xcalculado/- – Proccedings of V International Symposium on River
Xobservado) Sedimcntation. Editcd by P. Larsen & N. Eisenhauer,
x Error relativo de la variable X [(Xcalculado- Karlsruhe, Germany, April 1992. Volume I, pp. 163-
Xobservado)/Xobservado] 170.
Ф Función
 Función Farias,H.D. (1993a). "Morphology ofRegime A lluvial
 Viscosidad dinámica del fluido Channels: A Review". Advanccs in Hydro-Scicnce and
 Viscosidad cinemática del fluido Engineering. S.Y. Wang (Editor), University of
ρ Densidad del fluido Mississippi, USA, Vol.l, Part B, pp. 1423-1428.
ρ s Densidad del sedimento Farias,H.D. (1993b). "Menú de Ecuaciones de Diseño para
X Desviación estándar de la variable x Canales Arenosos en Régimen", Revista Latino-
g Desviación estándar geométrica de la distribución americana de Hidráulica. No. 6, pp. 1-12.
gnmulométrica
Farias,H.D. (1993c). "Forma de la Sección Transversal en
 Tensión de cizallamiento
Canales Aluviales Estables", Revista Latinoamericana
c Valor crítico de la tensión de corte para inicio de
de Hidráulica. No. 7, pp. 39-50.
arrastre
 Parámetro de Shields Farias,H.D. (1994). "Geometry of A lluvial Channels in
 *c Valor crítico de  Dynamic Equilibrium", Proc. International Symposium
on State of the Art in River Engineering Methods and
REFERENCIAS BIBLIOGRÁFICAS Design Philosophies. St. Petersburg, Russia, Vol. 3 .
Goncharov,V.N. (1962). Dynamics of Channel Flow.
Abou-Seida,M.M. & Saleh,M. (1987). "Design of Stable A
Israel Program for Scientific Translalions, Jerusa-lem,
lluvial Channels", Journal of Hydraulic Research.
Israel.
IAHR, Vol.25, No. 4, pp.433-446.
Karim,M.F. & Kennedy.J.F. (1990). "Menú ofCoupled
Ackers.P. (1983). "Sediment Transpon Problems in
Velocity and Sediment Discharge Relations for Rivers",
Irrigarían Systems Design", Chapter 4 in Develop-
Journal of Hydraulic Engineering. Vol. 16, No.8,
m.ents in Hydraulic Enginecring. P.Novak (Editor),
August, pp. 978-996.
Applied Science Publishers, London, U.K., Vol.l,
pp.151-195. Maza AlvarezJ.A. (1987). Introduclion to River
Engineering. Universitá Italiana per Stranieri, Perugia,
Bakker,B. & Vermaas.H. (1986). "Study to Establish
Italia.
Hydraulic Design Gritería for Pakistán Irrigation
System". Report PRC Engineering / Checchi, USAID Mao,S.W. & Flook.L.R. (1971). "Link Canal Design
Irrigation Systems Management Project, Lahorc, Practice in West Pakistán", Appcndix A in River
Pakistán. Mechanics. H.W.Shen (Ed.), Fort Collins, Colorado,
U.S.A., pp. A.1-A.40.
Brownlic,W.R. (1981). "Prediction of Flow Depth and
Sediment Discharge in Open Channels", Report KH-R- Pizzuto.J.E. (1984). "Equilibrium Bank Geometry and the
43A, California Institute Technology, Pasa-dena, U.S.A. Width of Shallow Sandbed Streams", Earth Surface
Processes and Landforms. Vol. 9, pp. 199-207.
Chítale,S.V. (1966). "Design of A lluvial Channels", 6th.
Congress of thc International Commission on Irrigation Rickard.C.E. (1989). "Regime Canal Design: Theory and
and Drainage (ICID), New Delhi, India, Report 17, Reality", International Confercnce on Irrigation.
Question 20. University of Southampton, 12-15 Sept. 1989, U.K., pp.
395-404.
Chow,V.T. (1959). Open Channel Hydraulics. McGraw-
Hill Book Co., New York, USA. Simons,D.B. & Albertson,M.L. (1963). "Unifoim Water
Conveyance Channels in A lluvial Material",
Farias,H.D. (1990). Modelos Matemáticos Variaciona-l.es
Transactions ASCE. Vol. 128, Part I, pp. 65-167.
para la Determinación de la Forma y Dimesio-ncs de
Canales Aluviales Estables. Tesis M.Sc., CIDIAT, Smith.K.V.H. (1970). "Similarity in Unlined Irrigation
Universidad de Los Andes, Mérida, Venezuela. Canal Systems", Journal of thc Hydraulics División.
ASCE, Vol. 96, No. HY1, pp. 13-28.
Farias.H.D. (1991). "A Unifica Approach to Extremal
Hypotheses for Morphology of Stable A lluvial Stevens,M.A. (1989). "Width of Straight A lluvial
Channels", Proc. XXIV Congress IAHR, Spain, pp. Channels", Journal of Hydraulic Engineering, ASCE,
Vol.l 15, No.3, pp. 309-326.

Ingeniería del Agua. Vol. 2 Num. 3 (septiembre 1995) p. 67

DIS EÑ O DE CANALES SIN REVESTIR

APÉNDICE 1. Deducción de las Ecuaciones de

Régimen a partir de las relaciones de fricción aluvial, m3  b4  c4 (1 b2 ) ( A18)
transporte sólido y estabilidad lateral.
m4  b2  c2 (1 b2 ) 1 ( A19)
Las ecuaciones que se disponen son las siguientes:

Reagrupando y manipulando algebraicamente estas

* Fricción Aluvial: V = a0 dal Ra2 Sa3 (A01) relaciones, se obtienen las siguientes ecuaciones de
régimen:
* Transporte Sólido: qs = b0 db1 Vb2 Rb3 Sb4 (A02)

P  x0 d x1C Sx2Q x3 ( A20)

* Estabilidad Lateral: P = c0 Qc1 dc2 Rc3 Sc4 (A03)

R  y 0 d y1C Sy 2Q y3
l0
( A21)

S
sujetas al cumplimiento de las ecuaciones de S  Z 0 d z1C z 2Q z3 ( A22)
conservación para flujo permanente uniforme:
Con:
* Caudal Líquido: Q =PRV (A04)

* Concentración de Sedimentos: Cs = q, P/ Q (A05) 

Z 0  m0 k 2
 (m3 k2 m2 )
( A23)
1
m 
0

A partir de la Ec. A05, puede escribirse:

(m  k m ) ( A24)
Z 1 1 2
1
C s  b0 d b1S b4 P (1b2 ) R (b3 b2 )Q (b2 1) (m  k m )
( A06) 3 2 2

De igual modo, a partir de la Ec. A04, se tiene: Z 2 1

(m  k m ) ( A25)

3 2 2
Q (1c2 )  a 0 c0 d (a1 c1 ) R (a2 c3 1) S (a3 c4 ) ( A07)

Combinando las ecuaciones A06 y A03, se obtiene: (m4  k3m2 ) ( A26)

Z
3
(m  k m )
3 2 2
b2 c2 (1b2 )1
C s b c0 0(1b2 )d b1 c1 (1b2 )R(b3 b2 b1 c3 )Sb4 c4 (1b2 )Q

( A08) y0  k0 Z k02 ( A27)

A partir de las Ees. A07 y A08, se obtienen las y1  k1  z1k 2 ( A28)

siguientes:

R  k 0 d k 1 S k 2 Qk 3 ( A09) y2  k3  z3k 2 ( A29)

Cs  m0 d m1 Rm3 S m3 Qm4 ( A10)

y 3  k1  z 1 k 2 ( A30)
Con:
x 0  c 0 y c3 Z c4 ( A31)
0
1 0
k 0 a c0 0
 ( A11)
(a2 c3 1)

k   (a1  c1 ) ( A12) x1  c1  c3 y1  c4 z1 ( A32)

1
(a 2  c3 1)
 DIS EÑ O DE CANALES SIN REVESTIR

( A13) x2  c3 y2  c4 z2 ( A33)
k   (a3  c4 )

2
(a 2  c3 1)

k   (1 c2 ) ( A14) x3  c2  c3 y3  c4 z3 ( A34)

3
(a 2  c3 1)

m0  b0c0 (1b2 ) ( A15)

m1  b1  c1 (1 b2 ) ( A16) Obsérvese que las Ees. A20, A21 y A22 tienen el formato
clásico de las ecuaciones de régimen para el diseño de
m2  b3  b2  c3 (1 b2 ) ( A17)
canales arenosos. Todas las Ees. presentadas en la Tabla
1 han sido reducidas a ese formato.

Ingeniería del Agua. Vol. 2 Num. 3 (septiembre 1995) p. 68

Ingeniería del Agua. Vol. 2 Num. 3 (septiembre 1995) p. 69