Artculo: COMEII-16015

II CONGRESO NACIONAL
DE RIEGO Y DRENAJE COMEII 2016
Chapingo, Edo. de Mxico, del 08 al 10 de septiembre

IMPLEMENTACIN DE ALGORITMO PARA DISEO DE PEQUEAS PRESAS

DE GAVIONES

Eduardo Jimnez Hernndez1

1Soluciones en Ingeniera y Tecnologas del Agua S.A. de C.V., Texcoco, Estado de Mxico.
eddojh@gmail.com. (*Autor de correspondencia).

Resumen

En este trabajo se presenta el desarrollo de una librera de software libre para el

diseo de pequeas presas de gaviones con propsitos de conservacin de
suelos. Misma que se implementa en un programa para facilitar la tarea de diseo.
La librera y el programa se desarrollaron en el lenguaje de programacin Python,
que es un lenguaje de alto nivel, libre, interpretado y multiplataforma. Se
implementaron algoritmos para el clculo del escurrimiento mximo del cauce,
(usando el mtodo racional o con el mtodo simplificado de la huella mxima),
para el diseo del vertedor rectangular de pared gruesa y para el anlisis
estabilidad en condiciones de deslizamiento y volcamiento. Se desarroll un
algoritmo propio para el dimensionamiento automtico de la presa. Los resultados
de los clculos obtenidos con el programa son iguales a los que se obtienen con el
uso de programas de hoja de clculo. El uso del programa demostr tener
ventajas en reducir el tiempo requerido para el diseo y los errores debidos a
descuidos. Adems, el programa proporciona al usuario un ambiente de trabajo
agradable y especializado para este propsito.

Palabras clave adicionales: Software libre, Python, conservacin de suelos,

erosin, presas de gaviones.
 II Congreso Nacional de Riego y Drenaje 2016, Chapingo, Edo. de Mxico, 08-10 de septiembre

Introduccin

La formacin de crcavas en el terreno es un indicador del grado avanzado de

erosin hdrica, su restauracin se hace por medio de diversas prcticas de control
de azolves, entre ellas se encuentran las presas de gaviones. El procedimiento de
diseo de estas obras, requiere de una previa planeacin de su ubicacin en el
cauce a ser restaurado, posteriormente se procede a realizar los clculos
necesarios para determinar sus dimensiones y estabilidad.

El dimensionamiento de presas de gaviones se realiza a partir de datos de la

seccin transversal del cauce en el sitio donde se desea construir la obra. Una vez
dimensionadas, el anlisis de estabilidad consiste en la aplicacin de ecuaciones
previamente establecidas. Por lo tanto, el procedimiento resulta sencillo de
codificar en un programa de cmputo.

Al no existir un programa para realizar el diseo de presas de gaviones, este

procedimiento se realiza con el uso de programas de cmputo como hojas de
clculo y de dibujo asistido por computadora. Sin embargo, resulta ser un
procedimiento largo y propenso a errores de clculo.

En el contexto anterior, el objetivo del presente trabajo fue implementar un

algoritmo para automatizar la realizacin de los clculos requeridos para el diseo
de presas de gaviones y generar los reportes y planos correspondientes. El
programa y su correspondiente librera se desarrollaron en el lenguaje de
programacin Python, que es un lenguaje de alto nivel, libre, interpretado y
multiplataforma.

Materiales y mtodos

La metodologa que se utiliz para el desarrollo de la librera para el diseo de

presas de gaviones es la presentada por Oropeza-Mota, (2012) en el curso de
Conservacin del suelo y del agua en el Departamento de Irrigacin de la
Universidad Autnoma Chapingo. Para propsitos de este trabajo, la metodologa
se dividi en diferentes partes, mismas que a su vez se codificaron en mdulos
que realizan una tarea especfica.

Escurrimiento mximo del cauce

Los escurrimientos mximos en pequeas cuencas o reas de drenaje pueden

estimarse mediante dos procedimientos (Colegio de Postgraduados, 1991):
1. Mtodo racional.
2. Mtodo simplificado de las huellas mximas.

Mtodo racional

Es un mtodo muy utilizado particularmente en el diseo de drenajes urbanos,

toma en cuenta el rea de la cuenca, la altura o intensidad de la precipitacin.

2
 II Congreso Nacional de Riego y Drenaje 2016, Chapingo, Edo. de Mxico, 08-10 de septiembre

El tiempo que transcurre entre el inicio de la lluvia y el establecimiento del gasto

de equilibrio se denomina tiempo de concentracin, y equivale al tiempo que tarda
el agua en pasar del punto ms alejado hasta la salida de la cuenca (Aparicio,
1992). Una forma de estimar el tiempo de concentracin (tc) es mediante la
frmula de Kirpich (ecuacin 1):

(1)

Donde tc es el tiempo de concentracin (en horas); Lc es la longitud del cauce

principal de la cuenca (en m); v es la velocidad media del agua en el cauce
principal (m s-1).

La ecuacin de Kirpich es vlida para cuencas cuya rea es mayor a 200 ha, para
cuencas de menor superficie se recomienda usar la ecuacin de Giandotti (1934)
(citado por Vlez & Botero, 2011), que se muestra en la ecuacin (2):

(2)

Donde tc es el tiempo de concentracin (h); A es el rea de la cuenca (km 2); Lc es

la longitud del cauce ms largo (km); S es la pendiente promedio del cauce
principal (m m-1).
La frmula racional se expresa como se muestra en la ecuacin (3):

(3)

Donde Ce es el coeficiente de escurrimiento, que representa la fraccin de la lluvia

que escurre de forma directa; i es la intensidad de la lluvia de una frecuencia
elegida para una duracin igual al tiempo de concentracin tc (Bedient, Huber, &
Vieux, 2008); Q es el gasto mximo posible que puede producirse con una lluvia
de intensidad i en una cuenca de rea A y coeficiente de escurrimiento Ce.

El coeficiente de escurrimiento (Ce) se puede determinar por tres mtodos, de los

que se seleccion el que est en funcin del tipo y uso de suelo y del volumen de
precipitacin anual de la cuenca de estudio, por su facilidad de programacin y
disponibilidad de datos (Poder Ejecutivo-SEMARNAT, 2002). Se clasifican los
suelos de la cuenca de estudio en tres diferentes tipos: A (suelos permeables), B
(suelos medianamente permeables) y C (suelos casi impermeables) que se
especifican en el Cuadro 1 y al tomar en cuenta el uso actual del suelo, se obtiene
el valor del parmetro K del Cuadro 2.

Si en la cuenca de estudio existen diferentes tipos y usos de suelo, el valor de K

se calcula como la resultante de subdividir la cuenca en zonas homogneas y

3
 II Congreso Nacional de Riego y Drenaje 2016, Chapingo, Edo. de Mxico, 08-10 de septiembre

obtener el promedio ponderado de todas ellas. Una vez obtenido el valor de K, el

coeficiente de escurrimiento anual Ce, se calcula mediante la ecuacin (4) si K
resulta menor o igual que 0.15 y con la ecuacin (5) si K es mayor que 0.15.

Cuadro 1. Tipos de suelo para determinar el valor de K.

Tipo de
Caractersticas
suelo
A Suelos permeables, tales como arenas profundas y loess poco compactos.
Suelos medianamente permeables, tales como arenas de mediana
B profundidad: loess algo ms compactos que los correspondientes a los
suelos A; terrenos migajosos.
Suelos casi impermeables, tales como arenas o loess muy delgados sobre
C
una capa impermeable, o bien arcillas.
Fuente: (Poder Ejecutivo-SEMARNAT, 2002).

(4)

(5)

Donde Ce es el coeficiente de escurrimiento, que representa la fraccin de la lluvia

que escurre de forma directa; K es un parmetro que se obtiene del Cuadro 2; P
es la precipitacin media anual (en mm).

Rango de validez: Las frmulas se considerarn vlidas para valores de

precipitacin anual (P) entre 350 y 2,150 mm.

Cuadro 2. Valores de K en funcin del tipo y uso del suelo.

Tipo de suelo
Uso de suelo
A B C
Barbecho, reas incultas y desnudas 0.26 0.26 0.26
Cultivos:
En hilera 0.24 0.24 0.24
Legumbres o rotacin de pradera 0.27 0.27 0.27
Granos pequeos 0.30 0.30 0.30
Pastizal (Porcentaje del suelo cubierto o pastoreo):
Ms del 75% -Poco- 0.14 0.14 0.14
Del 50 al 75 % -Regular- 0.20 0.20 0.20

4
 II Congreso Nacional de Riego y Drenaje 2016, Chapingo, Edo. de Mxico, 08-10 de septiembre

Tipo de suelo
Uso de suelo
A B C
Menos del 50% -Excesivo- 0.28 0.28 0.28
Bosque:
Cubierto ms del 75% 0.07 0.07 0.07
Cubierto del 50 al 75 % 0.16 0.16 0.16
Cubierto del 25 al 50 % 0.24 0.24 0.24
Cubierto menos del 25% 0.12 0.12 0.12
Zonas urbanas 0.22 0.22 0.22
Caminos 0.26 0.26 0.26
Pradera permanente 0.17 0.17 0.17
Fuente: (Poder Ejecutivo-SEMARNAT, 2002)

Mtodo simplificado de la huella mxima

El mtodo simplificado de la huella mxima tiene la ventaja de no requerir datos de

precipitacin y se recomienda en crcavas donde se observe claramente las
huellas dejadas por los escurrimientos mximos que se hayan presentado
(Colegio de Postgraduados, 1991).
El mtodo parte de la ecuacin de continuidad (ecuacin 6):

(6)

Donde Qmax es el gasto mximo de diseo (m3 s-1); Ah es el rea hidrulica del
cauce (m2); v es la velocidad del escurrimiento (m s-1).
La seccin transversal se obtiene por cualquier procedimiento topogrfico. El rea
hidrulica se obtiene por medio de figuras geomtricas conocidas. La velocidad
del escurrimiento se determina en funcin de la ecuacin de Manning (7):

(7)

Donde v es la velocidad del escurrimiento (m s-1); S es la pendiente del cauce (m

m-1); R es el radio hidrulico (m); n es el coeficiente de rugosidad de Manning.
Y el radio hidrulico se calcula mediante la relacin (8):

Ah
R (8)
Pm

Donde Ah es el rea hidrulica (m2) y Pm es el permetro mojado (m).

5
 II Congreso Nacional de Riego y Drenaje 2016, Chapingo, Edo. de Mxico, 08-10 de septiembre

Diseo del vertedor

En funcin del escurrimiento superficial calculado se dimensiona el vertedor. En

estructuras de gaviones, muros secos y mampostera hidrulica se utilizan
vertedores rectangulares de cresta gruesa, cuyas dimensiones se encuentran
mediante la ecuacin (9):

(9)

Donde, Hv es la carga sobre el vertedor (m); Q es el escurrimiento mximo de

diseo (m3 s-1); C es el coeficiente del vertedor (adimensional = 1.45) y L es la
longitud del vertedor (m).
Como el valor del escurrimiento mximo Q es conocido, dado que se determina
con alguno de los mtodos presentados anteriormente, se proponen valores de
longitud L y se obtiene como resultado la carga Hv.

Diseo de la colocacin de gaviones

Finalmente, el mtodo indica que se debe determinar la colocacin de gaviones en

funcin de la seccin transversal del cauce, as como el nmero y tamao de los
mismos. El programa realizado no incluye ninguna funcin de colocacin de
gaviones.

Diseo estructural

El clculo estructural de una presa de gaviones se realiza a partir de la seccin

crtica unitaria, es decir, que las dimensiones estn referidas a la unidad de ancho
del muro considerando el perfil de la Figura 1. La ubicacin de la seccin crtica
unitaria se realiza mediante el clculo del centroide de la presa de gaviones.

6
 II Congreso Nacional de Riego y Drenaje 2016, Chapingo, Edo. de Mxico, 08-10 de septiembre

Figura 1. Elementos de una seccin crtica unitaria tpica de una presa de gaviones

De acuerdo con la figura, se define la siguiente nomenclatura para las ecuaciones

que se presentan ms adelante:
H = altura de la presa (m).
h = altura de la carga sobre la cresta vertedora (m).
hi = altura del tendido i de la presa (m).
B = longitud de la base de la presa (m).
b = longitud de la corona de la presa (m).
K = longitud constante del escaln (m).
k' = longitud del colchn amortiguador (m).
bi = largo del tendido i (m).
q = peso de la lmina vertiente (t).
P = peso total de la seccin crtica unitaria (t).
Zp = brazo de palanca del peso de la seccin crtica unitaria (m).
E = empuje hidrosttico del agua (t).
= coeficiente de friccin (adimensional = 0.75 para piedra).
N, M = puntos de referencia aguas arriba y aguas abajo de la presa.
Un elemento que no aparece en la figura pero que es importante es el valor del
ancho de la seccin a = 1.

Peso de la lmina vertiente (q)

Es el peso del agua que pasa por el vertedor, de manera terica, es la lmina
formada por la longitud del vertedor, el ancho de la corona y la carga sobre la
cresta vertedora. Esto se expresa con la ecuacin (10):

(10)

Donde todas las variables se definieron en la nomenclatura de la Figura 1.

7
 II Congreso Nacional de Riego y Drenaje 2016, Chapingo, Edo. de Mxico, 08-10 de septiembre

Su brazo de palanca (Xq) con respecto al punto N de la Figura 1, se expresa con la

ecuacin (11):

(11)
Donde todas las variables se definieron en la nomenclatura de la Figura 1.

Peso de la seccin crtica unitaria (P)

Se hace el clculo del volumen unitario multiplicando la altura de cada tendido por
su ancho y sumando los resultados, esto se expresa con la ecuacin ((12). Ntese
que como el ancho del muro es unitario, este se omite en la ecuacin.

(12)

Donde V es el volumen de la seccin crtica unitaria, n es el nmero de tendidos

de la presa y las dems variables se definieron en la nomenclatura de la Figura 1.
El clculo del peso especfico aparente se hace mediante la ecuacin (13):

(13)

Donde a es el peso especfico aparente (t m-3); es el peso especfico de la

piedra (t m-3); es el peso especfico del agua con sedimentos (t m-3).
Finalmente, una vez que se ha obtenido el volumen unitario y peso especfico
aparente, se puede obtener el peso de la seccin crtica unitaria:

(14)

Donde P es el peso total de la seccin crtica unitaria (t) y V es volumen de la

seccin crtica unitaria (m3).
Se obtiene su brazo de palanca (XP) con respecto al punto N de la Figura 1
usando la ecuacin (15):

(15)
Donde Zp es el brazo de palanca del peso de la seccin crtica unitaria; bi es la
longitud de cada uno de los tendidos (m); hi es la altura de cada uno de los
tendidos (m); n = nmero de tendidos de la presa; V es el volumen de la seccin
crtica unitaria (m3).
Y Zp se obtiene con la ecuacin (16):

(16)

Donde todas las variables se definieron en la nomenclatura de la Figura 1.

8
 II Congreso Nacional de Riego y Drenaje 2016, Chapingo, Edo. de Mxico, 08-10 de septiembre

Empuje hidrosttico (E)

El empuje hidrosttico que ejerce el agua sobre la estructura se calcula con la
ecuacin (17):

(17)

Donde Sm es la superficie de mojado de la seccin crtica unitaria (m 2); h es el

centro de gravedad de la superficie de mojado (m). La superficie de mojado de la
seccin crtica unitaria se obtiene con la ecuacin ((18):

(18)

Y su centro de gravedad con la ecuacin ((19):

(19)

Su brazo de palanca (XE) con respecto al punto N de la Figura 1, se obtiene con la

ecuacin (20):

(20)

Condicin de deslizamiento
La condicin de deslizamiento establece que el peso del cuerpo de la presa y la
lmina de agua vertiente (afectados por el coeficiente de friccin) deben ser
mayores que el empuje hidrosttico.
(21)

Donde todas las variables se definieron en la nomenclatura de la Figura 1.

Condicin de volcamiento
La condicin de volcamiento a la que se debe someter la presa es que la relacin
entre el peso unitario de la misma y el empuje hidrosttico sea mayor a la unidad,
como se expresa en la ecuacin (22):

(22)

Donde XP es el brazo de palanca del peso de la presa; XE es el brazo de palanca

del empuje hidrosttico y las otras variables se definieron en la nomenclatura de la
Figura 1.

9
 II Congreso Nacional de Riego y Drenaje 2016, Chapingo, Edo. de Mxico, 08-10 de septiembre

Condicin de ncleo central

Establece que la resultante de todas las fuerzas debe estar en el tercio medio de
la base de la presa para que exista estabilidad. Esto se expresa en la ecuacin
(23):
(23)

Donde todas las variables se definieron en la nomenclatura de la Figura 1.

Figura 2. Diagrama de flujo para el desarrollo del programa de diseo de presas de

gaviones.

Con base en la metodologa expuesta, se procedi a desarrollar los algoritmos

necesarios en el lenguaje de programacin Python, llegando a obtener una librera
con funciones especializadas para el diseo de presas de gaviones denominada
gdamslib, adems de otros archivos requeridos para la interfaz grfica y el
programa principal, los cuales se agruparon en un sistema computacional
denominado Gabbioni.

10
 II Congreso Nacional de Riego y Drenaje 2016, Chapingo, Edo. de Mxico, 08-10 de septiembre

Anlisis y discusin de resultados

Como resultado principal de la implementacin del algoritmo de Oropeza-Mota,

(2012) para diseo de pequeas presas de gaviones se obtuvo un programa de
cmputo. Mismo que se integra de los siguientes mdulos: mdulo de entrada de
datos, dibujo de la seccin transversal del cauce, mdulo de clculo del
escurrimiento mximo usando a) el mtodo racional o b) el mtodo simplificado de
la huella mxima; mdulo de diseo del vertedor, mdulo de dimensionamiento de
la presa, dibujo de la presa, mdulo de anlisis de estabilidad y reporte de
resultados. Adems se incluye un mdulo bsico de dibujo de planos.

Una vez codificado el programa se procedi a comparar los resultados obtenidos

con un ejemplo realizado por Oropeza-Mota, (2012), mostrado en la Figura 3:

Figura 3. Vista de frente y de perfil de la seccin crtica unitaria.

La comparacin entre los resultados del ejemplo y los del programa se muestran a
continuacin:

Cuadro 3. Comparacin resultados entre un ejemplo y el programa.

Valor de comparacin Resultados ejemplo Resultados programa
Coordenada centroidal X 5.21 5.217
Coordenada centroidal Y 2.06 2.439
Coordenada centroidal Z 1.05 1.053
Peso de la lmina vertiente (t) 1.2 1.2
Peso de la seccin crtica unitaria (t) 9.6 9.6
Empuje hidrosttico (t) 5.4 5.4
Condicin de deslizamiento 1.5 1 1.5 1
Condicin de volcamiento 4.83 1 4.8333 1
Condicin de ncleo central 1.33 < 1.84 < 2.66 1.333 1.843 2.667

Adems de estos valores, los resultados de las condiciones contra volcamiento y

deslizamiento son positivos, es decir, la presa es estable ante estas condiciones.
Un resumen de los resultados del diseo de presas de gaviones que presenta el
programa se muestra en la Figura 4.

11
 II Congreso Nacional de Riego y Drenaje 2016, Chapingo, Edo. de Mxico, 08-10 de septiembre

Figura 4. Resumen de resultados del diseo de presas de gaviones.

Conclusiones

Con el software desarrollado se proporciona al usuario, una herramienta que

realiza de manera automtica los clculos del diseo de una presa de gaviones;
as como los dibujos o grficos necesarios. Los resultados del diseo realizado
con el programa son presentados en formato de texto e imagen, lo que facilita la
creacin de memorias de clculo y reportes.
Con el uso de esta herramienta el procedimiento de diseo de una presa de
gaviones se realiza en menos tiempo, comparado con el uso de programas de
hoja de clculo, dibujo asistido por computadora y procesador de textos por
separado. Adems, se reduce considerablemente la ocurrencia de errores
provocados por descuidos.
Despus de haber comparado los resultados obtenidos con el programa Gabbioni
y los clculos efectuados en el diseo con hojas de clculo se observa que
efectivamente el resultado es prcticamente igual. Las pequeas variaciones que
se observan son causadas por errores de redondeo al efectuar los clculos sin el
programa. Por lo anterior se concluye que el programa Gabbioni es una
herramienta confiable, precisa y funcional para el diseo de presas de gaviones

Referencias bibliogrficas

Aparicio, M. F. (1992). Fundamentos de hidrologa de superficie (primera ed.).

Mxico D.F: LIMUSA.

Bedient, P. B., Huber, W. C., & Vieux, B. E. (2008). Hydrology and foodplain
analysis. Prentice Hall.

Colegio de Postgraduados. (1991). Manual de conservacin del suelo y del agua

(tercera ed.). Chapingo, Mxico: Colegio de postgraduados.

12
 II Congreso Nacional de Riego y Drenaje 2016, Chapingo, Edo. de Mxico, 08-10 de septiembre

Oropeza-Mota, J. (2012). Notas del curso de Conservacin del suelo y del agua.
Mxico: Departamento de Irrigacin, Universidad Autnoma Chapingo.

Poder Ejecutivo-SEMARNAT. (2002). Conservacin del recurso agua. Norma

Oficial Mexicana NOM-011-CNA-2000. Diario Oficial de la Federacin.

Vlez, U. J., & Botero, G. A. (2011). Estimacin del tiempo de concentracin y

tiempo de rezago en la cuenca experimental urbana de La Quebrada San
Luis, Manizales. Dyna(165), 59-71.

13