INFILTRACIÓN POR ANILLOS INFILTRÓMETROS

Y BOTELLA DE MARIOTTE
Aguirre Hernán Alexis; Bolaños Wladimir; Esponda Maria del Mar;
Mosquera Katherine; Pechené Ossa Lisbeth; Montaño Leon Daniel.
Escuela de Ingeniería Agrícola - Universidad del Valle
Hidrología
Fecha: Junio 12 de 2018

RESUMEN
Se determinó la velocidad de infiltración y la conductividad hidráulica del suelo de la Universidad del
Valle ubicado en el Laboratorio de Aguas y Suelos Agrícolas (LASA), utilizado como método los
anillos infiltrómetros, además se realiza por medio de la botella de Mariotte la velocidad de
infiltración en la estación experimental del departamento de Biología. Para el cálculo se utilizó el
método de Kostiakov el cual se compara con el modelo de Horton, Philp y Green-Ampt. Según los
datos obtenidos se encontró que para el suelo del laboratorio del LASA cuenta con una velocidad de
infiltración (I = 111.16 mm/h) y una conductividad hidráulica (Ks = 3.79mm/h) donde el (Ks) es
moderadamente lenta propia de un suelo Arcillo-Limoso y para el suelo del la estación experimental
de biología la velocidad de infiltración (I = 32,66 mm/h) y una conductividad hidráulica (Ks = 0,26
mm/h) donde el (Ks) corresponde a un suelo arcilloso.

Palabras Claves: Conductividad hidráulica, velocidad de infiltración, anillos infiltrómetros, botella

de Mariotte.

1. INTRODUCCIÓN
El frasco o botella de Mariotte es un dispositivo creado para calcular de manera eficiente la velocidad
de efusión constante de un determinado liquido, para realizar dicho experimento se debe abrir un
orificio en la parte lateral de la botella el cual debe ser de igual longitud al que va en la tapa del
dispositivo esto con el fin de regular de manera efectiva el efecto que causa la presión atmosférica
dentro y fuera de la botella y así poder establecer un goteo constante y regulado con el fin de medir el
tiempo promedio en que se desocupa el recipiente y tener los datos de aforo. Históricamente, el
principio del Frasco de Mariotte fue utilizado habitualmente en el siglo XIX en los quinqués y
lámparas del alumbrado doméstico de Argan, en los que, a pesar de tener el depósito más alto que el
nivel por donde se alimenta la llama, se impide el salto del aceite en forma de surtidor como debería
hacerlo si no fuera por esta ingeniosa disposición.

2. OBJETIVOS
2.1. OBJETIVOS GENERALES
● Calcular la velocidad de infiltración por los diferentes métodos.
● Calcular los cambios en los porcentajes de humedad del suelo por medio del sistema de goteo.

2.2. OBJETIVOS ESPECÍFICOS

● Implementar el uso de la botella de Mariotte previamente aforada como sistema de riego por
goteo.
● Determinar la infiltración en el suelo de las instalaciones del laboratorio de hidráulica, por
medio del método de anillos infiltrómetros.
3. METODOLOGÍA
3.1. Botella de Mariotte.
 Para conocer el caudal de la botella de Mariotte de 3 litros se utilizó una base metálica, mecanismo
artesanal de goteo constante (regulador de goteo) y una probeta de 1000 ml, se lleva a cabo el aforo de
la botella por un intervalo de tiempo de 9 horas 24 minutos (9.4 h). Luego de saber el tiempo que
tardó la botella en vaciarse por medio de un goteo constante se lleva a campo y se instala en un lugar
previamente seleccionado que cumpla con algunas especificaciones, como que: no haya presencia de
árboles cerca o que el terreno sea empedrado, esto con el fin de no afectar el proceso de infiltración
del agua en el suelo.

Se escogió el suelo de la estación experimental del departamento de biología de la universidad del

Valle sede Meléndez, se adaptó el lugar y se ubicó la botella aforada entre 5-10 cm sobre el suelo para
luego dar inicio la prueba registrando los tiempos y lámina de agua en la botella de Mariotte cada 20
minutos donde se revisaba o ajustaba el goteo (para así generar un goteo constante a lo largo de toda
la práctica), al finalizar se mide el radio de saturación y humedecimiento.

Imagen 1. Botella de Mariotte.

Además se tomaron 3 muestras de suelo por medio de una pequeña calicata a una profundidad de
10, 20 y 30 cm respectivamente, el proceso tuvo una duración de 8 horas, logrando que se consumiera
un volumen de 2540 ml . De igual forma se tomó muestra del suelo antes del goteo para ver la
humedad inicial del suelo, estas también se hicieron a 10, 20 y 30 cm en un zona cercana al lugar del
montaje. Estas muestras con y sin goteo fueron llevada a el laboratorio de Química Ambiental del
programa académico de Ingeniería Sanitaria y ambiental al horno por un tiempo de 24 horas a una
temperatura de 110 °C.
 Imagen 2. Muestra del suelo seco.

3.2. Anillos infiltrómetros.

Para la práctica de anillos infiltrómetros se escogió un suelo cerca al Laboratorio de Aguas y

Suelos Agrícolas (LASA) el cual cumplió con las mejores condiciones para realizar el ensayo,
primero se quita la capa vegetal presente en el sitio con ayuda de un palín, después con ayuda de la
placa de accionamiento y el martillo de goma se golpea el anillo interior hasta una profundidad de 15
cm, para el cilindro exterior se realiza la misma operación; con los anillos instalados en el sitio se
pone una regla en la pared del anillo interior seguidamente se cubre con una bolsa plástica el anillo
interior para luego llenar con agua el cilindro exterior hasta una altura aproximada entre 5-10 cm,
enseguida se llena con agua el cilindro interior hasta alcanzar una altura entre 5-10 cm retirando la
bolsa plástica se da inicio a la prueba registrando los datos de altura y tiempo de infiltración cada 30
seg, 1 y 5 min teniendo en cuenta de posibles recargas de agua de los anillos interior y exterior (para
procurar que ambos cilindros tuviesen niveles de agua muy parecidos), la prueba tuvo un tiempo de
duración de 1 hora.

Imagen 3. Instalación de anillos infiltrómetros.

4. DATOS Y RESULTADOS
Para determinar la lámina acumulada se usó la ecuación de Kostiakov:
L=K ×t m Ec.1. Ecuación de Kostiakov
Donde:
L = Lámina Acumulada en el tiempo (mm)
K = Constante que depende de la estructura inicial del suelo
m = Constante que depende de la estabilidad de la estructura del suelo frente al agua, 0 > m < 1.

K en suelos arenosos o francos con grietas presentan valores entre 10 a 30. Suelos con estructuras
estables frente al agua muestran valores de m mayores a 0.6 .
Ahora se transforma la ecuación exponencial (2) en una recta por medio de aplicarle log a ambos
lados de la ecuación:
log L=log K + mlog t
Esta ecuación corresponde a la ecuación de la recta y = b + mX, donde log L es la variable
dependiente, log K es la ordenada al origen, m corresponde a la pendiente de la recta y log t hace
referencia a la variable independiente.
Para determinar la pendiente se usará la ecuación de mínimos cuadrados:

Ec. 2. Mínimos cuadrados.

4.1. Botella de Marriott

Lectura de
Tiempo Lectura de Lámina Lámina Lámina Tiempo Lámina Lámina Lámina
Lámina en
Hora acumulad Lámina en la parcial acumulad acumulad Hora acumulado parcial acumulad acumulad
la botella
o [min] botella [cm] [cm] a [cm] a [mm] [min] [cm] a [cm] a [mm]
[cm]
7:34 0 56,0 0 0 0 11:54 260 42,3 1,6 13,7 137
7:54 20 55,3 0,7 0,7 7 12:14 280 41,4 0,9 14,6 146
8:14 40 54,8 0,5 1,2 12 12:34 300 40,6 0,8 15,4 154
8:34 60 53,4 1,4 2,6 26 12:54 320 39,8 0,8 16,2 162
8:54 80 52,6 0,8 3,4 34 1:14 340 38,9 0,9 17,1 171
9:14 100 51,7 0,9 4,3 43 1:34 360 38,2 0,7 17,8 178
9:34 120 50,5 1,2 5,5 55 1:54 380 37,6 0,6 18,4 184
9:54 140 49,5 1 6,5 65 2:14 400 36,9 0,7 19,1 191
10:14 160 46,4 3,1 9,6 96 2:34 420 36,3 0,6 19,7 197
10:34 180 46,0 0,4 10 100 3:54 440 35,5 0,8 20,5 205
10:54 200 45,0 1 11 110 3:14 460 34,8 0,7 21,2 212
11:14 220 44,2 0,8 11,8 118 3:34 480 33,5 1,3 22,5 225
11:34 240 43,9 0,3 12,1 121
Tabla 1. Datos de la botella de Marriott.
 Gráfico 1. Lámina acumulada vs tiempo para la botella de Marriott.

Después de aplicar la Ec.1 obtenemos la siguiente información.

Nº T [min] L [mm] Log t Log L Log t.Log L (Log t)2 Nº T [min] L [mm] Log t Log L log t.Log L (Log t)2

1 20 7 1,30 0,85 1,10 1,69 13 260 137 2,41 2,14 5,16 5,83
2 40 12 1,60 1,08 1,73 2,57 14 280 146 2,45 2,16 5,30 5,99
3 60 26 1,78 1,41 2,52 3,16 15 300 154 2,48 2,19 5,42 6,14
4 80 34 1,90 1,53 2,91 3,62 16 320 162 2,51 2,21 5,54 6,28
5 100 43 2,00 1,63 3,27 4,00 17 340 171 2,53 2,23 5,65 6,41
6 120 55 2,08 1,74 3,62 4,32 18 360 178 2,56 2,25 5,75 6,53
7 140 65 2,15 1,81 3,89 4,61 19 380 184 2,58 2,26 5,84 6,66
8 160 96 2,20 1,98 4,37 4,86 20 400 191 2,60 2,28 5,94 6,77
9 180 100 2,26 2,00 4,51 5,09 21 420 197 2,62 2,29 6,02 6,88
10 200 110 2,30 2,04 4,70 5,29 22 440 205 2,64 2,31 6,11 6,99
11 220 118 2,34 2,07 4,85 5,49 23 460 212 2,66 2,33 6,19 7,09
12 240 121 2,38 2,08 4,96 5,67 24 480 225 2,68 2,35 6,31 7,19
TOTAL 55,02 47,25 111,65 129,11
Tabla 2. Datos procesados de la botella de Mariotte.
 Gráfico 2. Velocidad de infiltración obtenida con la botella de Mariotte.

Reemplazando en la Ec.2 obtenemos, que:

m = 1,12
Recordando que en la ecuación de la recta (anteriormente descrita) tenemos como variable
independiente (x) a Log t y como variable dependiente (y) a Log L.

47,25 55,02
b=Y Promedio −m X Promedio → b = −1,12.
24 24
b=−0.59
K= Antilog b → K= Antilog−0.59
K=0,26
Finalmente la ecuación que describe la infiltración (mm) en función del tiempo (min) es:
L=0,26T 1,12

La velocidad de infiltración es:

I =Km T m−1 → I =0,29 T 0,12 [mm/min]
I =60. Km T m−1 → I =17,26 T 0,12 [mm/h]

Tiempo acumulado [min] I [mm/h] Tiempo acumulado [min] I [mm/h]

20 24,727 260 33,639
40 26,871 280 33,939
60 28,211 300 34,221
80 29,202 320 34,487
100 29,994 340 34,739
120 30,658 360 34,978
140 31,230 380 35,206
160 31,735 400 35,423
180 32,186 420 35,631
200 32,596 440 35,831
220 32,971 460 36,022
240 33,317 480 36,207
Tabla 3. Velocidad de infiltración determinada por la Ec.1..

La tasa de infiltración para este suelo fue de 32,66 mm/h

 Profundidad [cm] 10 20 30
Suelo húmedo ( S❑H ) [g] 18.56 57.68 66.95

Suelo seco ( S❑S ) [g] 14.00 44.20 51.92

%W 32.59 % 30.50 % 28.95 %
Tabla 4. Datos de contenido de humedad inicial a distintas profundidades.

Profundidad [cm] 10 20 30
Suelo húmedo ( S❑H ) [g] 25,98 48,99 69,97

Suelo seco ( S❑S ) [g] 19.48 38.03 54.50

%W 33.37 % 28.82 % 28.39 %
Tabla 5. Datos de contenido de humedad final a distintas profundidades.

Gráfico 3. Contenido de humedad vs Caudal antes y después del goteo (caudal de 0,319 l/h).

Diámetro de humedad [cm] 35,0

Diámetro de encharcamiento [cm] 20,0
Tabla 6. Diámetro de humedad y encharcamiento.

Por otro lado asumiendo que el disco formado en la superficie adquiere una geometría circular tenemos que la
velocidad de infiltración corresponde a la siguiente ecuación:
Qe
f= Ec. 3 Velocidad de infiltración.
π ⋅r 2

Tomando f como la velocidad de infiltración hallada por el método de Kostiakov tenemos que, el radio del área
por donde pasa el determinado caudal es:
r =√ ❑= 5.58 cm

Según este valor se tiene un error relativo porcentual de:

10,0−5.58
E R %= ⋅100
10,0
E R %=79.211%
4.2. Anillos infiltrómetros.

Tiempo Tiempo Láminas Lámina Tiempo Tiempo Láminas Lámina

Lectura Enrase Lectura Enrase
Hora parcial acumula Parciale acumula Hora parcial acumula Parciale acumula
[cm] [cm] [cm] [cm]
[s] do [s] s [cm] da [cm] [s] do [s] s [cm] da [cm]
8:53:00 0 0 9 9,00 0,00 0,00 9:03:00 60 600 6,90 0,20 2,10
8:53:30 30 30 8,70 0,30 0,30 9:04:00 60 660 6,70 0,20 2,30
8:54:00 30 60 8,60 0,10 0,40 9:05:00 60 720 6,50 0,20 2,50
8:54:30 30 90 8,40 0,20 0,60 9:06:00 60 780 6,40 0,10 2,60
8:55:00 30 120 8,30 0,10 0,70 9:07:00 60 840 6,20 0,20 2,80
8:55:30 30 150 8,20 0,10 0,80 9:08:00 60 900 5,80 0,40 3,20
8:56:00 30 180 8,10 0,10 0,90 9:13:00 300 1200 4,90 0,90 4,10
8:56:30 30 210 8,00 0,10 1,00 9:18:00 300 1500 8,40 9,30 0,90 5,00
8:57:00 30 240 7,90 0,10 1,10 9:23:00 300 1800 7,70 0,70 5,70
8:57:30 30 270 7,75 0,15 1,25 9:28:00 300 2100 6,80 0,90 6,60
8:58:00 30 300 7,60 0,15 1,40 9:33:00 300 2400 6,10 0,70 7,30
8:59:00 60 360 7,50 0,10 1,50 9:38:00 300 2700 5,40 0,70 8,00
9:00:00 60 420 7,45 0,05 1,55 9:43:00 300 3000 10,10 10,20 0,10 8,10
9:01:00 60 480 7,30 0,15 1,70 9:48:00 300 3300 9,30 0,80 8,90
9:02:00 60 540 7,10 0,20 1,90 9:53:00 300 3600 8,50 0,80 9,70
Tabla 7. Datos por el método de los anillos infiltrómetros.

Gráfico 4. Lámina acumulada vs tiempo para los anillos infiltrómetros.

Después de aplicar la Ec.1 obtenemos la siguiente información.

Log t . (Log Log t .

No. t[min] L [mm] Log t Log L No. t[min] L [mm] Log t Log L (Log t)2
Log L t)^2 Log L
1 1 4 0,00 0,60 0,00 0,00 15 11 23 1,04 1,36 1,42 1,08
2 1,5 6 0,18 0,78 0,14 0,03 16 12 25 1,08 1,40 1,51 1,16
3 2 7 0,30 0,85 0,25 0,09 17 13 26 1,11 1,41 1,58 1,24
 4 2,5 8 0,40 0,90 0,36 0,16 18 14 28 1,15 1,45 1,66 1,31
5 3 9 0,48 0,95 0,46 0,23 19 15 32 1,18 1,51 1,77 1,38
6 3,5 10 0,54 1,00 0,54 0,30 20 20 41 1,30 1,61 2,10 1,69
7 4 11 0,60 1,04 0,63 0,36 21 25 50 1,40 1,70 2,38 1,95
8 4,5 12,5 0,65 1,10 0,72 0,43 22 30 57 1,48 1,76 2,59 2,18
9 5 14 0,70 1,15 0,80 0,49 23 35 66 1,54 1,82 2,81 2,38
10 6 15 0,78 1,18 0,92 0,61 24 40 73 1,60 1,86 2,99 2,57
11 7 15,5 0,85 1,19 1,01 0,71 25 45 80 1,65 1,90 3,15 2,73
12 8 17 0,90 1,23 1,11 0,82 26 50 81 1,70 1,91 3,24 2,89
13 9 19 0,95 1,28 1,22 0,91 27 55 89 1,74 1,95 3,39 3,03
14 10 21 1,00 1,32 1,32 1,00 28 60 97 1,78 1,99 3,53 3,16
TOTAL 28,08 38,19 43,58 34,90
Tabla 8. Procesamiento de datos.

Gráfico 5. Velocidad de infiltración obtenida para los anillos infiltrómetros.

Resolviendo por mínimos cuadrados, se tiene:

m=0.7827
Recordando que en la ecuación de la recta (anteriormente descrita) tenemos como variable
independiente (x) a Log t y como variable dependiente (y) a Log L.

b=Y Promedio −m X Promedio

38.19 28.08
b= −0.7827 .
28 28
b=0.58
K= Antilog b
K= Antilog 0.58
K=3.79
Finalmente la ecuación que describe la infiltración (mm) en función del tiempo (min) es:
L=3.79 T 0.7827

La velocidad de infiltración es:

I =Km T m−1 → I =2.97 T −0.2173 [mm/min]
I =60. Km T m−1
→ I =178.13 T −0.2173
[mm/h]

Aplicando esta última ecuación a cada uno de los tiempo se obtuvo la columna de I [mm/h], en la
siguiente tabla.
Tiempo acumulado [min] I [mm/h] Tiempo acumulado [min] I [mm/h]*
1 178.13 11 105.78
1.5 163.11 12 103.80
2 153.22 13 102.01
2.5 145.97 14 100.38
3 140.30 15 98.8
3.5 135.67 20 92.89
4 131.79 25 88.49
4.5 128.46 30 85.06
5 125.55 35 82.25
6 120.67 40 79.90
7 116.70 45 77.88
8 113.36 50 76.12
9 110.50 55 74.56
10 107.99 60 73.16
*método de Kostiakov.
Tabla 9. Velocidad de infiltración método de Kostiakov para los anillos infiltrómetros.

Finalmente el promedio de las velocidades fue 111.16 mm/h

A continuación se observan tres tipos de gráficas con diferentes ajustes para la gráfica 4.

Gráfica 6. Infiltración acumulada vs. tiempo con ajuste exponencial.

Gráfica 7. Infiltración acumulada vs. tiempo con ajuste lineal.

Gráfica 8. Infiltración acumulada vs. tiempo con ajuste polinómico.

A continuación se observan tres tipos de gráficas con diferentes ajustes para la gráfica 5.

Gráfica 9 .Velocidad de infiltración instantánea vs. tiempo con ajuste exponencial.

Gráfica 10 .Velocidad de infiltración instantánea vs. tiempo con ajuste lineal.

Gráfica 11 .Velocidad de infiltración instantánea vs. tiempo con ajuste polinómico.

 Según el coeficiente de determinación (r ❑2) el mejor ajuste se obtiene del ajuste polinómico. con una
ecuación de la forma y=−2 E−7 x 2+ 0.0035 x +0.2025(1) , esto para la gráfica de infiltración
acumulada vs tiempo.
También para la gráfica de la tasa de infiltración vs el tiempo encontramos que la función que mejor
se ajusta respecto a otra es la función polinómica y=0,036 x 2−3,3164 x +152,76( 2)

Respecto a la ecuación polinomial de grado dos encontrada para la grafica de la lamina acumulda vs
tiempo podemos decir que esta tiende a ser cóncava hacia abajo debido al coeficiente que acompaña a
la variables x2, esto podría significar que en el tiempo la lámina tiende a disminuir su altura pero
llegado su punto de mínima altura, vuelve aumentar tal vez debida a que el anillo exterior comparte
un conexión hidráulica con el anillo interior, lo cual ayuda este posible aumento después del punto
donde la lámina es mínima, en conclusión se planteamos la ecuación en forma general ax 2 + bx +c = y,
los terminó a y b que son negativos ayudan a comprender que la lámina de agua siempre tiende a
disminuir
Ahora respecto a la ecuación de la velocidad instantanea vs tiempos podemos inferir que el valor en el
coeficiente b que es negativo genera cierta reducción en la velocidad de infiltración aunque esta
función jamás será negativa y posera un valor de velocidad en dos tiempos distintos, esto se evidencia
al momento de graficar esta funcion. Cabe resaltar que esta gráfica no es la que mejor se ajusta debido
a que el valor de R2, no es muy cercano a 1 y por ellos el modelo no se ajusta tan bien al modelo, esto
talvez se debe a que los datos son muy dispersos.
Como dato extra y según estas funciones se puede encontrar el valor de x donde esta el foco de la
función que para las ecuaciones (1) y (2), son respectivamente 46.06 y 8750 segundos .

En la siguiente tabla se compara el modelo de la ecuación de polinómica de la velocidad de

infiltración vs el tiempo, respecto a varios modelo de infiltración.

Ecuación Modelo de Modelo de Modelo de Modelo de Modelo de Holtan

elegida Philip Horton. kostiakov Green-Ampt

Expresión y= A x 2 + Bx+Cf =s t −1/ 2+C f =f c+(f o−f c) ef−kt

=α t β f = A+C / F f =f c❑+ a ¿
general

El mejor ajuste que se le puede realizar a la gráfica es el ajuste polinómico (gráfica 11) ya que la
función que la describe es de la forma F (x)=a x n siendo a y n valores reales, a equivale al valor de k
(3.79) y n al de m (0.78). La función de la gráfica es:
L=3.79 T 0.7827

K= Constante que depende de la estructura inicial del suelo (en seco). Es la lámina que se infiltra en el
primer instante mayor que cero (mm). Esto se explica ya que al emplear log a ambos lados de la
ecuación que describe la infiltración de la lámina tenemos que obedece a la ecuación de la recta y k se
relaciona a b, la cual hace referencia al intercepto en el eje y (ordenadas) el cual en este caso sería la
lámina infiltrada en el minuto 0, lo que se traduce en el momento inicial.
Lo anteriormente dicho se confirma con : “K en suelos arenosos o francos con grietas presentan
valores entre 10 a 30.Suelos con estructuras estables frente al agua muestran valores de m mayores a
0.6”[1].
m= Constante que depende de la estabilidad de la estructura del suelo frente al agua, 0 > m < 1
Entre mayor sea m, mayor capacidad tendrá en suelo de infiltrar agua en el tiempo.

Ahora respecto al modelo que menos ajusta la función obtenida es la expresión podría ser la modelo
de Horton debido a que esta es una expresión que tiene la fórmula ajustada de :
f = A e−kt lo cual tiene un decrecimiento constante, mientras que la ecuación que polinómica de grado
2 posee un decrecimiento y un crecimiento con punto mas bajo en el eje x dado por x = b/2a, mientras
que el punto más b1ajo de la ecuación de horton es cuando el valor de t que es tiempo tiene al infinito.

Como sugerencia para calibrar mejor la función obtenida se recomienda realizar el experimento en un
tiempo mucho más prolongado, el cual sea suficiente para estabilizar los valores de infiltración en el
suelo y de esta forma poder estimar un comportamiento más realista.

Por otro lado suponiendo que el suelo sobre el cual se desarrolló la práctica es de tipo franco-arenoso
y un contenido inicial de humedad de 0.15, estime la infiltración acumulada y la tasa de
infiltración usando el modelo del Green-Ampt.

Esta ecuación es conocida como la ecuación de Green-Ampt para la infiltración acumulada. Una vez
que F ha sido calculada, la tasa de infiltración f puede obtenerse como:

Por otra parte, en el caso de que la profundidad de encharcamiento h0 no sea despreciable, el valor ψ -
h0 se sustituye por ψ en las anteriores ecuaciones. La ecuación de Green-Ampt es no lineal para F y
puede resolverse mediante el método de sustituciones sucesivas.

La aplicación del modelo de Green-Ampt requiere estimaciones de la conductividad hidráulica (K), de

la porosidad (n) y de la carga de succión de suelo en el frente de mojado (ψ). Experiencias de
laboratorio con una amplia variedad de suelos han mostrado que ψ puede expresarse como una ley de
potencias de la saturación efectiva y se define como:

Donde θr es el contenido residual de humedad de suelo después de que éste haya sido drenado
completamente. Mientras que θ-θr es la humedad disponible y n-θr es la porosidad efectiva (θe).
La saturación efectiva está contenida en el rango 0 ≤se≤ 1 siempre que θr ≤ θ ≤ n. Entonces, el
cambio de contenido de humedad ∆θ cuando pasa el frente de mojado es igual a:

La siguiente tabla muestra los rangos típicos para n, θe, ψ y K.

 Tabla 10. Parámetros de infiltración de Green-ampt para varias clases de suelos.

5. ANÁLISIS DE RESULTADOS

5.1 Análisis de la botella de Mariotte

La humedad del suelo en la estación experimental antes y después de llevar a cabo el goteo por medio
de la botella de Mariotte, demostró que el suelo sí se encuentra más húmedo a mayor profundidad
como se puede observar en el gráfico 3. También se observa que no hay mucha variación entre los
porcentajes de humedad del suelo antes y después del experimento debido a que el suelo contiene un
porcentaje de humedad que depende del tiempo, el para esa semana resultó ser muy lluvioso, lo cual
podría generar una alteración en los datos obtenidos.
De acuerdo con la velocidad media de infiltración de infiltración hallada de I = 32,66 mm/h, se logra
obtener la relación que tiene con el radio de encharcamiento (350 mm), debido a que esta velocidad es
relativamente pequeña y dicho radio es grande. Este dato determina también el tipo de suelo en el que
se está trabajando que para este caso es un suelo pesado, lo cual es corroborado al obtener un valor
para la conductividad hidráulica (Ks) menor a 1, dando como resultado un tipo de suelo arcilloso.

En el gráfico 2, se observa que la curva de los datos es ascendente y además es cóncava hacia abajo,
esto es debido a que el agua que está goteando es constante y por lo tanto el suelo siempre va a recibir
la misma cantidad de agua (observar gráfico 1, donde los datos obtenidos forman una línea no muy
dispersa). También se observa que la velocidad de infiltración se hace constante a partir del minuto
400 cuando I = 35,423 mm/h.
Cuando los datos de la tabla 2 se reemplazaron en la Ec. 2 para determinar el valor de la pendiente
(m), se observó que está dio un valor mayor a 1 el cuál presente inconsistencias según la teoría donde
se afirma que: m debe estar entre 0 y 1. Está inconsistencia se debe al tiempo con el que el suelo
absorbe la lámina de agua ya que prácticamente estos dos datos son constantes, como se mencionó
anteriormente el goteo es controlado de forma que, la frecuencia del goteo siempre sea la misma.
El valor obtenido por medio de la ecuación 3 entrega el área teórica por donde fluye el agua, por tanto
se asume como punto de referencia para evaluar el modelo de Kostiakov anteriormente realizado,
dando un error relativo porcentual de 79.211% lo cual claramente significa que no es un método
confiable para predecir el área por donde circulará el agua infiltrada, ya que el radio de
encharcamiento que se obtuvo es de 10 cm, además también se tiene el area de humedad por donde
también ocurre infiltración de agua.
Lo anterior suscita el hecho de ver la infiltración en el suelo como un hecho que depende de otras
variables, las cuales pueden ser el tipo de suelo, el microrelieve que se presenta el la zona, así como la
humedad inicial del suelo.

Estimación de la capacidad de adsorción del suelo estimado.

Haines ha desarrollado un concepto teórico de cohesión en un suelo ideal sobre la base de las fuerzas
de tensión superficial que surgen de las películas de agua entre las partículas.

El suelo ideal es considerado como formado por esferas uniformes que están arregladas en uno u otro
tipo de empaquetamiento (abierto o cerrado). En el abierto hay 6 puntos de contacto y en el cerrado
hay 12 puntos de contacto por partícula (Imagen 4). A bajos contenidos de humedad la mayoría de las
películas de agua se encuentran como anillos alrededor de los puntos de contacto. Cada una de estas
películas tiende a unir las partículas. La fuerza cohesiva total es equivalente a la suma de todas las
fuerzas individuales ejercidas por cada película en cada punto de contacto Analizando los trabajos de
Haines, Fisher ha demostrado que la fuerza total de unión ejercida por cualquier película de agua es
igual a la presión ejercida por la interfase agua-aire. La fuerza total la determina la siguiente fórmula:

2 xaT
F=
1+ tan(0.5)

a = radio de las partículas, T = tensión superficial (fuerza con la cual las moléculas del líquido se
atraen entre sí, ocupando la menor superficie), 0 = ángulo del menisco de agua.

Imagen 4. Cohesión del suelo en relación con las fuerzas de tensión superficial.
Nichols ha demostrado que la cohesión aumenta hasta un máximo y luego decrece rápidamente al ir
aumentando el contenido de humedad; esto se debe a que la cohesión es función del número y espesor
de los “films”.

Nichols trabajando con suelos sintéticos obtuvo los siguientes datos. La cohesión aumenta con el
contenido de arcilla pero decrece con el de humedad. Esto se debe a que los contenidos de humedad
más bajos citados para cada caso (Tabla 11) corresponden a los puntos de máxima cohesión.

Variación de la cohesión con la humedad Atterberg y sus colaboradores fueron los primeros en
estudiar los cambios en las propiedades cohesivas de los suelos como función de la humedad. La
técnica consiste en hacer ladrillos con el suelo, secándose a varios contenidos de humedad diferentes
y determinando la fuerza necesaria para hacer penetrar una cuña de metal en la masa. Los resultados
están ilustrados en las curvas A y B de la gráfica (Imagen 5). La cohesión aumenta con el descenso en
el contenido de humedad; esto se debe al decrecimiento con el espesor del film de agua. La arcilla
tiene una cohesión mucho mayor que el franco arenoso fino, como resultado de mayor cantidad de
“films” y mayor superficie de contacto.

Para el tipo de suelo franco-arenoso y con una humedad de 0.15 se obtuvo una infiltración acumulada
F = 113.36045 cm y una tasa de infiltración de I = 11.96232 cm/h; estos resultados permiten concluir
que una inflamación muy rápida debido a su alta porosidad con presencia de arenas que permiten un
movimiento del agua muy ligero.

SUELO Contenido de humedad Cohesión [gramos/in]

10,90 17,25
Arena ⅔, arc 1/3
12,90 15,00
12,73 26,40
Arena ⅓, arc 2/3
13,10 22,50
13,55 56,00
Arena fina pura
17,50 49,00
Tabla 11.Relación de cohesión con el contenido de humedad del suelo.

Imagen 5 . Relación de la cohesión y de la adhesión con la humedad del suelo.Las curvas A y B han
sido tomadas de Johannsen (1914), quien empleó suelos batidos; las curvas C y D son de los datos de
Nichols (1929, 1931), quien partió de material seco y suelto y lo humedece paulatinamente.
5.2. Análisis Anillos infiltrómetros.
Con respecto a los datos obtenidos y analizando la serie de gráficas que se encontró podemos analizar
que la gráfica de lámina acumulada de agua vs tiempo presenta unos datos buenos aceptables, los
cuales se ajustan a la un función polinómica descrita anteriormente y que describe un valor de R 2 muy
cercano al dos, lo cual implica que la función se ajusta muy bien a los datos obtenidos. igualmente
podemos agregar que coeficiente que acompaña a la variables x 2, podría significar que en el tiempo la
lámina tiende a disminuir su altura pero llegado su punto de mínima altura, vuelve aumentar tal vez
debida a que el anillo exterior comparte un conexión hidráulica con el anillo interior, lo cual ayuda
este posible aumento después del punto donde la lámina es mínima, en conclusión se planteamos la
ecuación en forma general ax 2 + bx +c = y, los terminó a y b que son negativos ayudan a comprender
que la lámina de agua siempre tiende a disminuir .Todo este análisis se realiza de acuerdo con las con
la gráfica:

Gráfica 11. Ecuación que tiene el mejor ajuste de .con los datos de lamina acumula vs tiempo

Ahora respecto a la ecuación de la velocidad instantanea vs tiempos podemos inferir que el valor en el
coeficiente b que es negativo genera cierta reducción en la velocidad de infiltración aunque esta
función jamás será negativa y posera un valor de velocidad en dos tiempos distintos, esto se evidencia
al momento de graficar esta funcion. Cabe resaltar que esta gráfica no es la que mejor se ajusta debido
a que el valor de R2, no es muy cercano a 1 y por ellos el modelo no se ajusta tan bien al modelo, esto
talvez se debe a que los datos son muy dispersos.
la gráfica de la función que tuvo el mejor ajuste y sobre la cual se hizo el análisis fue la siguiente:
 Gráfica 12. Ecuación que tiene el mejor ajuste de .con los datos de lamina acumula vs tiempo

Como dato extra y según estas funciones se puede encontrar el valor de x donde esta el foco de la
función que para las ecuaciones (1) y (2), es decir las gráfica si fueran completamente ajustadas a los
funcion descritas en la gráficas 11 y 12, tendrán como puntos de tiempo de alcance maximo y minimo
los 46.06 s y 8750 segundos respectivamente.
También se puede analizar que que la tasa de infiltración dio con datos tan dispersos tal vez por las
siguientes causas:
● El terreno no era completamente plano ni homogéneo.
● Debido a que días anterior a la práctica se tuvo fuertes lluvias en el lugar de medicion, esto
puede influir a que las velocidad sean tan variantes.
Respecto al modelo de infiltración que mejor se ajustó que fue el de
Holtan, podemos decir que este sea el posible modelo que describa mejor
la infiltración del suelo en la zona, así mismo la expresión general del
n
modelo f = f c + a(So − F) ;donde a y n son constantes que dependen del tipo de suelo,
superficie y el tipo de cultivo, y So es el almacenamiento potencial por encima de la capa de
impedimento y la expresión (So-F) representa el potencial de infiltración, teniendo en cuenta todo lo
anterior se debería realizar un estudio del suelo y del cultivo que se desee producir.

6. CONCLUSIONES
La capacidad de infiltración varía de acuerdo al tipo de suelo, en suelos desagregados y permeables,
ésta es mayor que en suelos arcillosos y con índices de compactación.

Para el suelo de la estación experimental se obtuvo un suelo arcilloso, y al ser un suelo de este tipo
presentara mala aireación, baja infiltración y fácil encharcamiento, lo que dificultará las labores de
labranza. Como recomendación para determinar la humedad correcta del suelo, el experimento
debería de realizarse durante una época en donde las precipitaciones no estén tan presentes o realizar
diferentes pruebas para diferentes épocas del año y así poder establecer un sistema de riego óptimo.
La adhesión de un objeto extraño a un suelo debería darse solamente a contenidos de humedad por
encima del máximo de cohesión. A estos elevados contenidos de humedad, el agua estaría retenida
menos fuertemente por las partículas y sería atraída sobre la superficie del objeto para formar films de
enlace entre el objeto y el suelo. La adhesión del suelo hacia el objeto es por medio de estos films. El
contenido de humedad al cual ocurre la máxima adhesión depende del total de agua requerida para
satisfacer los “films” entre partículas individuales y las fuerzas atractivas en la superficie de los
objetos extraños.

Se puede concluir que los anillos de infiltración son de gran utilidad al momento de conocer,
características propias de lo suelos como son la tasa de infiltración , es importante resaltar que para
hacer una mejor práctica y por lo observado se debería conocer el tipo de suelo, la conductividad
hidráulica, además que para tener datos que sean mejores se debe realizar en temporada seca para que
el suelo tal vez no presente inconvenientes por saturación.
También para un practica futura se puede hacer un análisis se la infiltración pero aplicando el modelo
de Holtan, claro está que se debería consultar en la literatura cuales son las condiciones en las cuales
se podría aplicar este método y cuales son sus beneficios

7. BIBLIOGRAFÍA

● Jaramillo, D. (2002). Introducción a la ciencia del suelo. Medellín. Recuperado de:

http://www.bdigital.unal.edu.co/2242/1/70060838.2002.pdf.
● Instructivo, MEDICIÓN DE LA INFILTRACIÓN EN EL SUELO MEDIANTE
INFILTRÓMETROS DOBLE ANILLO. Recuperado de:
http://www.prosap.gov.ar/Docs/INSTRUCTIVO%20_R014_%20infiltrometro%20doble
%20anillo.pdf
● UNIVERSIDAD NACIONAL. "SANTIAGO ANTÚNEZ DE MAYOLO". Recuperado de:
http://repositorio.unasam.edu.pe/bitstream/handle/UNASAM/1173/T%20190%202014.pdf?
sequence=1&isAllowed=y