Introducción
El riego por goteo es el sistema de riego localizado a presión más eficiente en el suministro de agua ―a veces con fertilizante para el suministro a los cultivos― cuando se aplica en cantidad y tiempo adecuados. El suministro de agua y nutrientes se hace por medio de emisores insertados en una red extensa de tuberías, la cual se debe diseñar para su correcto funcionamiento. Los laterales de riego son parte de esta red de tuberías y para su diseño se requieren estimar con precisión las pérdidas por fricción y localizadas.
Las ecuaciones para determinar las pérdidas por fricción se han desarrollado con base en numerosas pruebas de laboratorio, en tuberías de distinto diámetro y tipo, las cuales han sido propuestas para expresar la relación entre la velocidad del agua, el diámetro y la rugosidad interna del tubo (Christiansen, 1942), donde la ecuación que más se utiliza para determinar las pérdidas por fricción en riego por goteo es la de Darcy Weisbach (ecuación (1)) debido a su exactitud y completa gama de aplicaciones:
El único inconveniente para utilizar la ecuación de Darcy Weisbach radica en el cálculo del factor de fricción (f), el cual se solucionó con la ecuación empírica de Blasius (ecuación (2)), para tuberías de diámetro inferiores a 25 mm, con número de Reynolds entre 3 000 < Re < 100 000 (Howell & Barinas, 1980; Hathoot, Al-Amoud, & Mohammad, 1993; Brown, 2003; Juana, Rodriguez-Sinobas, & Losada, 2002):
El número de Reynolds (Re) puede ser calculado con las siguientes expresiones:
O equivalente a:
La determinación de pérdidas por fricción en tuberías con diámetro interno inferior a 25 mm para laterales de riego se estima con la ecuación (5), la cual es el resultado de la combinación algebraica entre las ecuaciones (1) y (4), con temperatura del agua a 18 °C:
El lateral de riego es una tubería con salidas múltiples y la determinación exacta de las pérdidas por fricción en este tipo de tuberías requiere de un análisis tramo por tramo a partir de la última salida y teniendo en cuenta el caudal al inicio del lateral (Figura 1).
Para facilitar el análisis, algunos investigadores han propuesto expresiones que agilizan dicho proceso, entre ellos se destacan Christiansen (1942); DeTar (1982); Scaloppi (1988); Cuenca (1989); Anwar, (1999a); Anwar, (1999b); Chinea y Domínguez (2006), y los más recientes, Ángeles, Carrillo, Ibáñez, Arteaga y Vázquez, (2009); sin embargo, la ecuación que más se utiliza es la de Christiansen:
Las pérdidas por fricción en tuberías con salidas múltiples se determinan con la ayuda de la siguiente ecuación:
a partir de la hipótesis de que las pérdidas menores que se generan por la deformación del tubo y la barra sólida del emisor insertado se pueden despreciar. Algunos diseñadores de sistemas de riego han ocasionado deficiencias en dichos sistemas, pues un número considerable de emisores puede convertirse en un porcentaje importante en las pérdidas de carga total (Yildirim, 2010). Por lo anterior, algunos investigadores han propuesto expresiones para determinar dichas pérdidas, haciendo relaciones entre el área del tubo y la que ocupa el gotero dentro del tubo (Bagarello et al., 1997; Provenzano & Pumo, 2004; Demir, Yurdem, & Degirmencioglu, 2007).
Las pérdidas menores (h m ) se expresan en la ecuación clásica de la carga cinética multiplicada por el coeficiente de pérdida locales (∝).
Bagarello et al. (1997), en su investigación experimental, deducen un procedimiento para evaluar las pérdidas locales debido a la protuberancia de la parte sólida del emisor en la línea del flujo de riego por goteo. El coeficiente de las pérdidas locales (ecuación (10)) corresponde a diferentes sistemas de emisores sobre línea, con diferentes valores de número de Reynolds, y cada sistema de emisor se caracterizó por un índice de obstrucción (OI) (ecuación (9)):
El coeficiente de pérdidas locales ∝ (ecuación (11)) incluye ambas características: la de conexión del emisor y la deformación de la tubería que rodea el vástago (parte insertada del emisor), definida en los rangos 1.0 < A p /A g < 1.44 (Demir et al., 2007):
Howell y Barinas (1980) trabajaron con seis diferentes tipos de emisores sobre línea, donde establecieron la siguiente expresión empírica para determinar las pérdidas localizadas como longitud equivalente:
Los parámetros (C e = 0.0000066 y M = 1.32) dependen del tipo de emisor sobre línea.
Keller y Bliesner (1990) proponen la ecuación (13) para expresar las pérdidas de carga por conexión en longitud equivalente del lateral:
Esta última ecuación tiene el inconveniente de calcular la longitud equivalente de pérdidas de carga por conexión por medio de un gráfico que se debe reajustar (f e = 0.18).
En este trabajo se muestran los resultados de la investigación experimental sobre las pérdidas localizadas que se generan en laterales de riego con emisores sobre línea, después de verificar la aplicabilidad de las diferentes ecuaciones para determinar las pérdidas por fricción y localizadas descritas anteriormente.
Materiales y métodos
En el experimento se utilizaron dos laterales de polietileno (PE) de 50 metros de longitud, con diámetro interior de 13.5 mm; en el primer lateral se insertaron goteros sobre línea espaciados 0.40 metros; y en el segundo se insertaron a cada 0.40 metros en el primer tramo (10 metros), y 0.35 metros los restantes 40 metros del lateral (Figura 1).
El área de la tubería (A p ) y el área que ocupa el gotero dentro de la tubería (A g ), que se utilizan en la estimación de las pérdidas localizadas, se evaluaron con la digitalización de cinco fotografías tomadas al sistema tubería-gotero, como muestra la Figura 2, con la ayuda del programa AutoCAD 2007.
El diámetro interno de la tubería de PE se midió en cinco cortes de muestras de 0.10 metros por medio del calibrador vernier de pie de rey digital (con 0.02 mm/0.001” de exactitud). Esta medida se realizó dos veces por cada muestra, posteriormente se establece el promedio.
En la Figura 3, a cada 10 metros se instaló una toma de presión para hacer tres lecturas y estimar su promedio, y se compararon dichas pérdidas (por fricción y localizadas) obtenidas experimentalmente y las calculadas con las ecuaciones anteriormente descritas (Cuadro 1).
Metodología | Pérdidas en tubería con salidas múltiples | Pérdidas menores | Pérdidas totales en lateral |
---|---|---|---|
Bagarello et al. (1997) | H f = h f F |
|
H f = h f * F + h m |
Howell y Barinas (1980) | H f = h f F | l e = C e * Q M | H f = h f * F + l e |
Keller y Bliesner (1990) | H f = h f F | (f
e
= 0.18) y |
H f = h f * F + J |
Keller y Bliesner (1990) modificada | H f = h f F |
|
H f = h f * F + J * |
Las pérdidas localizadas experimentalmente se calcularon como la diferencia entre la presión arrojada cada 10 metros y las pérdidas por fricción que se generaron en cada tramo. Los datos se analizaron con la ayuda del programa Excel 2010 de Microsoft y la herramienta Cftools del programa Matlab 2010 (curve fitting tool).
La presión de operación de la entrada del lateral fue de 1.83 kg cm-2, que se controló con una válvula de compuerta y se registró en un manómetro de glicerina (Figura 4). Los goteros empleados fueron autocompensados de 8 l/h, marca Naandanjain, cuya presión de apertura es de 8 m.c.a.
Se determinó la ecuación del emisor o relación carga-gasto con presiones superiores e inferiores a su presión nominal de trabajo, sometiendo a los goteros a 2.11, 1.76, 1.41, 1.05 y 0.7 kg cm-2, con intervalos que no rebasaron los cinco metros.
Resultados
Perdidas por fricción
Registradas las presiones en las cinco tomas, ubicadas cada 10 metros (m) a lo largo de la tubería de polietileno (PE) en los dos laterales evaluados, se establece la distribución de las mismas (Cuadros 2 y 3).
Promedios de presión (kg cm-2) | % longitud | Longitud acumulada (m) | Pérdidas de carga cada 10 m1 | Pérdidas acumuladas |
---|---|---|---|---|
1.828 | 0.00 | |||
1.195 | 20.00 | 10 | 48.65% | 48.65% |
0.861 | 40.00 | 20 | 25.68% | 74.32% |
0.668 | 60.00 | 30 | 14.86% | 89.19% |
0.562 | 80.00 | 40 | 8.11% | 97.30% |
0.527 | 100.00 | 50 | 2.70% | 100.00% |
Pérdidas porcentuales de presión cada 10 metros respecto al rango de presión (1.828-0.527 kg cm-2) primer tramo;
[(1.828-1.195) /1.301] * 100 = 48.65%.
Promedio de presión (kg cm-2) | % Longitud | Longitud acumulada (m) | Pérdidas de carga cada 10 m1 | Pérdidas acumuladas |
---|---|---|---|---|
1.828 | 0.00 | |||
1.160 | 20.00 | 10 | 50.00% | 50.00% |
0.773 | 40.00 | 20 | 28.95% | 78.95% |
0.536 | 60.00 | 30 | 17.76% | 96.71% |
0.501 | 80.00 | 40 | 2.63% | 99.34% |
0.492 | 100.00 | 50 | 0.66% | 100.00% |
Pérdidas porcentuales de presión cada 10 metros respecto al rango de presión (1.828-0.492 kg cm-2) primer tramo;
[(1.828-1.160) /1.336] * 100 = 50.00%
Estos resultados de la columna 4 de los Cuadros 2 y 3 son similares a los publicados por Boswell (1990), citado por Andrade (2009), el cual dice que: “la caída media de la presión total se produce en los primeros 22% de la longitud y la presión media está a 38% de la longitud de entrada, cuando el lateral está a nivel”.
De igual forma, en los últimos 10 metros del lateral (80 y 100% l) ocurre entre el 3.36 y 0.66% de las pérdidas totales, resultados equivalentes a los publicados por Valiantzas (2005), en el cual manifiesta que: “el último 20% de la longitud de la manguera de goteros contribuye aproximadamente con el 1% del total de las pérdidas por fricción a lo largo del lateral”.
Pérdidas menores
Las pérdidas menores se establecieron utilizando las ecuaciones de la columna 3 del Cuadro 1. En esta columna se propone la ecuación (14) para calcular la longitud equivalente de pérdidas de carga por conexión después de analizar los datos experimentales. La herramienta Cftools de Matlab 2010 mostró que el mejor método experimental fue el polinomio cuadrático, con un coeficiente de determinación (R
2 = 0.9997), en comparación con el Gaussiano y el potencial. Es importante mencionar que después de realizar la digitalización de las fotografías, se estimó en 1.375 la relación de áreas
En el Cuadro 4 se comparan las pérdidas que se generaron de forma experimental o la caída de la presión a lo largo del lateral (columna 1) con los resultados, empleando diferentes metodologías (columnas 6, 7, 8 y 9). Se afirma que las metodologías de Bagarello et al. (1997), y la propuesta de Keller y Bliesner (1990), modificada en este trabajo (columnas 6 y 9), son las que mejor estiman las pérdidas en el lateral tradicional, con errores inferiores al 11.40% en los primeros 20% de la longitud de la tubería, y errores inferiores al 6.0% en los tramos comprendidos entre 80 y 100% de la longitud de la tubería. Es importante recordar que en este tramo ya se ha producido más de 97% de las pérdidas de carga total (columna 4, Cuadro 2).
Presión1 | Pérdidas | Diferencias | ||||||
---|---|---|---|---|---|---|---|---|
(m.c.a.) | (HT) 2 (m) | (H T ) 3 (m) | (H T ) 4 (m) | (H T ) 5 (m) | (D T ) 6 (m) | (D T ) 7 (m) | (D T ) 8 (m) | (D T ) 9 (m) |
18.28 | ||||||||
11.95 | 5.09 | 4.49 | 4.54 | 4.98 | 13.19 | 13.79 | 13.74 | 13.30 |
8.44 | 3.23 | 3.01 | 2.95 | 3.24 | 8.73 | 8.94 | 9.00 | 8.71 |
6.42 | 1.75 | 1.76 | 1.65 | 1.81 | 6.69 | 6.67 | 6.79 | 6.63 |
5.62 | 0.70 | 0.80 | 0.69 | 0.75 | 5.72 | 5.62 | 5.73 | 5.66 |
5.19 | 0.12 | 0.17 | 0.13 | 0.14 | 5.50 | 5.45 | 5.50 | 5.49 |
Suma | 10.88 | 10.23 | 9.94 | 10.92 |
Presión a lo largo del lateral.
Pérdidas totales metodología Bagarello et al. (1997).
Pérdidas totales metodología Howell y Barinas (1980).
Pérdidas totales metodología Keller y Bliesner (1990).
Pérdidas totales metodología Keller y Bliesner (1990) y modificada.
Diferencia entre la presión manométrica (columna 1) y las pérdidas totales de la columna 2.
Diferencia entre la presión manométrica (columna 1) y las pérdidas totales de la columna 3.
Diferencia entre la presión manométrica (columna 1) y las pérdidas totales de la columna 4.
Diferencia entre la presión manométrica (columna 1) y las pérdidas totales de la columna 5.
Sin embargo, de estas dos metodologías, la que mejor estimó las pérdidas por conexión en el lateral tradicional es la propuesta por Keller y Bliesner (1990) modificada, pues alcanzó un error promedio en todo el lateral de 4.85%.
En el Cuadro 5 se comparan las pérdidas generadas de forma experimental (columna 1) con los resultados, empleando cuatro metodologías (columnas 6, 7, 8 y 9). Se confirma que las metodologías propuestas por Bagarello et al. (1997), y por Keller y Bliesner (1990) (columnas 6 y 8) son las que mejor estiman las pérdidas en el lateral propuesto, arrojando errores inferiores a 6.0% en los primeros 20% de la longitud del lateral, y errores inferiores al 13% en los tramos comprendidos entre 60 y 80% de la longitud de la tubería. En dicho tramo ya se ha producido más de 96% de las pérdidas de carga total (columna 5, Cuadro 3).
Presión1 | Pérdidas | Diferencias | ||||||
---|---|---|---|---|---|---|---|---|
(m.c.a) | (HT) 2 (m) | (H T ) 3 (m) | (H T ) 4 (m) | (H T ) 5 (m) | (D T ) 6 (m) | (D T ) 7 (m) | (D T ) 8 (m) | (D T ) 9 (m) |
18.28 | ||||||||
11.60 | 6.56 | 5.61 | 6.02 | 6.66 | 11.72 | 12.67 | 12.26 | 11.62 |
7.73 | 4.40 | 3.93 | 3.93 | 4.35 | 7.20 | 7.67 | 7.68 | 7.26 |
5.36 | 2.45 | 2.38 | 2.26 | 2.50 | 5.28 | 5.35 | 5.47 | 5.23 |
5.01 | 0.94 | 1.04 | 0.9 | 1.00 | 4.42 | 4.32 | 4.46 | 4.36 |
4.92 | 0.17 | 0.23 | 0.17 | 0.19 | 4.84 | 4.78 | 4.84 | 4.82 |
Suma | 14.52 | 13.19 | 13.28 | 14.70 |
Presión a lo largo del lateral.
Pérdidas totales, metodología Bagarello et al. (1997).
Pérdidas totales, metodología Howell y Barinas (1980).
Pérdidas totales, metodología Keller y Bliesner (1990).
Pérdidas totales, metodología Keller y Bliesner (1990) y modificada.
Diferencia entre la presión manométrica (columna 1) y las pérdidas totales de la columna 2.
Diferencia entre la presión manométrica (columna 1) y las pérdidas totales de la columna 3.
Diferencia entre la presión manométrica (columna 1) y las pérdidas totales de la columna 4.
Diferencia entre la presión manométrica (columna 1) y las pérdidas totales de la columna 5.
No obstante, de estas dos metodologías, la que mejor predijo las pérdidas por conexión en el lateral propuesto fue la de Keller y Bliesner (1990), logrando un error promedio en todo el lateral de 4.20%.
Conclusiones
La mejor ecuación para predecir las pérdidas por conexión en los laterales evaluados fueron la propuesta por Keller y Bliesner (1990), modificada en este trabajo, pues alcanza errores promedios en todo el lateral inferiores a 5.0%. De igual forma, esta ecuación tiene la versatilidad de emplearse para calcular las pérdidas menores tramo por tramo (procedimiento laborioso, que se efectúa con hojas de cálculo), y también se puede utilizar como factor multiplicador con las pérdidas por fricción en tuberías con salidas múltiples para estimar las pérdidas totales, procedimiento que no se puede llevar a cabo con las otras ecuaciones utilizadas, al ser necesario calcular las pérdidas por conexión tramo por tramo.