ENSAYO DE INFILTRACION CON DOBLE ANILLO

HIDROLOGIA
Trabajo final

Por:
Andrés Fabian Higuera Bautista
115529
Freddy Alonso Pérez Gonzales
11353
Edna Daniela Carmona Martínez
115513

JORGE JULIAN VELEZ UPEGUI

Docente

UNIVERSIDAD NACIONAL DE COLOMBIA

Sede Manizales
2019

1. INTRODUCCION
El suelo es el principal ambiente sobre el que se desarrolla la vida, con múltiples
usos desde económicos hasta recreativos, todos esenciales para la vida en la tierra,
posee diferentes características con las cuales es posible clasificarlo y así realizar
su correcto uso en determinadas acciones, la estructura y la forma de sus partículas
son dos importantes características físicas, de estos parámetros se puede partir
para comenzar a decir que tipo de suelo se tiene.
La cantidad de agua que entra en un suelo desde la superficie hasta las capas más
inferiores es una variable importante a tener en cuenta para determinar la aptitud de
un suelo para determinados objetivos, a este fenómeno se le llama infiltración;
depende de factores como la porosidad, granulometría y permeabilidad, los cuales
de igual forma controlaran la estructura y forma del suelo, parámetros útiles para
cómo se mencionó anteriormente comenzar a clasificar el suelo. La infiltración
también es importante en consideraciones ambientales como flujos subterráneos,
tiempo de riego, humedad del suelo etc.… y demás fenómenos relacionados como
la escorrentía.
El presente informe recoge la información acerca de la experiencia de un ensayo de
infiltrometro de doble anillo realizado en un terreno ubicado en el laboratorio de
hidráulica de la Universidad Nacional de Colombia sede Manizales, además de una
comparación de esta práctica con la información que puede obtenerse al aplicar
distintos modelos empíricos como Green y Ampt, Kosiakov, Horton etc. para concluir
cual de estos pude acercarse más a la información obtenida en campo.
2. OBJETIVOS
Mediante los conceptos académicos estudiados sobre la infiltración, realizar en
campo un ensayo de infiltrometro de doble anillo y calcular para el terreno la tasa
de infiltración.
• Investigar acerca de la norma ASTM (D3385 – 94) Double ring, entender el
procedimiento y realizarlo en campo para un suelo determinado.
• Realizar un gráfico de infiltración contra tiempo en unidades de milímetros
por segundo, con la información obtenida en campo.
• Aplicar algunos modelos empíricos y graficar junto con la información de
campo para obtener la ecuación que más se ajuste a los datos reales.

3. MARCO TEORICO
La infiltración es el proceso en el cual el agua superficial penetra a las capas internas
del suelo debido a fuerzas gravitatoria o de tipo capilar.
Es determinante en la recarga de acuíferos y afecta directamente la escorrentía,
puede volver a la atmosfera por evaporación o transpiración o quedar retenida en el
suelo.
Parámetros como la cubierta vegetal de la zona, la intensidad y duración de la lluvia
y las características del suelo como el estar compuesto de un material permeable
ayuda a que la infiltración se desarrolle, se requiere cierta granulometría y una
interconexión entre los poros del material.

Ilustración 1 Evolución de la tasa de infiltración

La capacidad de infiltración o velocidad con la que el agua penetra en el suelo, se

expresa en mm/h y depende directamente del grado de humedad de sus capas, si
las superiores están
saturadas, se dará en
función de la
permeabilidad de estratos
inferiores.
La humedad es medible
mediante el punto de
marchitez donde el agua Ilustración 2 Humedad en el suelo

contenida en el suelo no se puede usar por las plantas y con la capacidad de campo
donde el agua ya no circula por gravedad.

Método de doble anillo

Según Selker et al. (1999), el método del doble
anillo es la prueba más común para la
estimación de propiedades hidráulicas del suelo.
Para su determinación los dos anillos (exterior e
interior) se introducen en el suelo a una
profundidad de 10 a 15 cm, hasta lograr que
estén perfectamente enterrados y nivelados. El
espacio interno de los anillos se llena de agua.
Existen dos tipos básicos de pruebas: carga
constante o carga variable. En esta investigación
se utilizó el caso de carga variable, en el que se Ilustración 3 Modelo ensayo de doble
anillo
permite que el nivel del agua descienda a lo largo
de la prueba. El método consiste en tomar lecturas del descenso del nivel de agua
en el anillo interior. Cuando se consigue que la tasa de la infiltración sea constante,
la prueba se termina.

Modelos empíricos
Algunas de las metodologías más conocidas y útiles en la determinación de la
velocidad de infiltración son las siguientes:
Horton (1939), establece que la curva que representa la capacidad de infiltración se
manifiesta, solo si la intensidad de precipitación es mayor que la capacidad de
infiltración del suelo analizado. Este principio o postulado es debidamente respetado
en el ensayo de infiltración de doble anillo, porque siempre hay agua en superficie
que satisface la capacidad de absorción que tiene el suelo en su grado máximo o
potencial, cualquiera sea el tiempo que transcurre desde el inicio del ensayo.
Kostiakov (1932), propuso un modelo exponencial, el cual desarrollo de manera
empírica, como se muestra en la siguiente ecuación: 𝑖(𝑡) = 𝛼𝑡 −𝛽 , de gran utilidad
en las practicas de riego, cuyas constantes empíricas tienen como interpretación
física la variación de la lámina de agua y la tasa de infiltración en el tiempo.
Mezencev (1948), vence las limitaciones de la ecuación de Kostiakov, que no ofrece
buenos resultados en el caso de periodos de mediciones extendidos.
Philip (1957), soluciono numéricamente la ecuación de Richards, estudio la
infiltración vertical comparándola inicialmente con la horizontal, donde no tiene
efecto la gravedad; la solución indica que, para tiempos pequeños, la infiltración es
proporcional a la raíz cuadrada del tiempo (infiltración horizontal).
Green-Ampt (1911), aunque simple en el proceso, es la primera expresión de
infiltración con base física, se publicó en 1911. Sigue siendo muy utilizada, pues
considera los aspectos fundamentales, y ha sido importante para el desarrollo de la
hidrología y la física del suelo. Aunque presenta y adaptabilidad en su aplicación
para diferentes escenarios, queda por determinar con precisión los parámetros que
lo definen.
4. DESCRIPCION DEL METODO
El 13 de febrero del presente año, el
ensayo se llevó a cabo en horas de
la tarde, en un terreno ya preparado
en el laboratorio de hidráulica de la
misma universidad, los anillos del
infiltrometro se encontraban ya
instalados en el terreno; se realizó el
montaje del equipo compuesto por
dos cilindros calibrados o tubos de
Mariotte de diámetros diferentes que
Ilustración 4 Montaje del ensayo
ayudaran a tomar los datos del
procedimiento, estos deben estar nivelados en el terreno. Con mangueras de
plástico se conectan los cilindros a los anillos del infiltrometro donde el cilindro de
menor diámetro ira conectado al anillo de menor diámetro y de la misma forma con
el otro tubo. Se debe de tener cuidado con el aire que queda dentro de las
mangueras expulsando con ayuda del agua todo residuo antes de comenzar con el
ensayo ya que puede afectar las mediciones.
Posterior a esto se inicia el vaciado de agua en los anillos hasta cierto nivel de
referencia marcado, se espera que parte del líquido pase a los tubos y se estabilice
en un nivel constante (cabe resaltar que se alcanza el nivel constante de los tubos
en un mismo valor), cuando esto ocurre se empieza a tomar el tiempo de medición.
Se tomarán datos por 3 horas continuas de la siguiente manera, la primera hora se
tomarán datos de variación del nivel en los cilindros cada 15 minutos; la siguiente
hora cada 10 minutos y finalmente cada 5 minutos. El nivel del agua debe
mantenerse en su estado inicial, en este caso, en el que se registró el nivel del agua
a medida que esta se infiltraba, la idea fue no dejar que se secaran completamente
los tubos adicionándole agua en una ocasión.
A medida que avanzaba el ensayo se observan los
cambios en la tasa de infiltración; el suelo ya había
sido estudiado previamente por otras personas para
el mismo procedimiento porque lo que se parte de la
premisa de que se encontraba saturado. Al finalizar
el procedimiento de medición se procede a
desmantelar el montaje y a realizar los cálculos y la
interpretación de los datos registrados. Ilustración 5 Procedimiento

5. MATERIALES
• Anillos concéntricos o anillos del infiltrometro: Estos deben ser de un material
resistente y en las ocasiones que lo requiera de materiales fuertes que
soporten fluidos corrosivos. Su diámetro y altura deben ser proporcionales a
las medidas de los tubos de Mariotte.

Ilustración 6 Anillos del

infiltrometro

• Cilindros graduados o tubos de Mariotte: Son tubos de

plástico de altura proporcional al contenido que
pueden mantener los anillos, uno junto al otro unidos
por una misma base para realizar fácilmente la
calibración, deben tener cada uno un grifo que
obstruya el paso del agua cuando sea necesario, por
este mismo deben poder introducirse las mangueras
de plástico que conectaran con los anillos del
infiltrometro, también para practicidad en el
procedimiento es útil que tengan sobre su superficie
una cinta métrica o decámetro que permita determinar Ilustración 7 Tubos de
las variaciones en la elevación, no importan las Mariotte
 unidades que maneje ya que estas pueden ser convertidas en el
procedimiento de cálculos.
• Cronometro: Para mantener un control sobre los tiempos entre mediciones,
buscando obtener mejores resultados para las diferentes horas, se
establecieron diferentes intervalos de tiempo para tomar un dato.
• Agua: Es necesario contar con un suministro de agua cercano para mantener
el nivel requerido en ambos anillos.
• Mangueras de platico: De diferentes longitudes para permitir conectar
fácilmente cada anillo con un tubo respectivo.

Ilustración 8 Mangueras

• Recipiente: De ser necesario tener un control sobre la cantidad de agua

que se vierte en los anillos, se debe contar con un recipiente debidamente
graduado para tener más precisión en el manejo del volumen y la lámina
de agua.

6. RESULTADOS
Luego de la toma de datos se recoge la siguiente información, en las columnas 3 y
4 el tubo Nº1 corresponde al anillo exterior (el tubo de mayor diámetro) y el tubo Nº2
al anillo interior (el tubo de menor diámetro), estos representan los niveles de agua
en los cilindros pasados los intervalos de tiempo definidos, con esta información en
las siguientes columnas, se realizo la diferencia entre los datos, con esto ya
podemos encontrar los valores de la infiltración:

Infiltración Tiempo en
Hora Tiempo (min) Nivel del agua Diferencias
(mm)1 Horas
Tubo 1 Tubo 2 Tubo 1 Tubo 2
3:25 14,7 14,7 0
3:40 13,7 13,9 1 0,8 8 0,25
3:55 15 12,7 13,3 1 0,6 6 0,5
4:10 11,8 12,6 0,9 0,7 7 0,75
4:25 10,9 11,9 0,9 0,7 7 1
4:35 10,2 11,5 0,7 0,4 4 1,16
10
4:45 9,7 11,1 0,5 0,4 4 1,33

1
Solo del anillo interior.
 4:55 9,1 10,7 0,6 0,4 4 1,5
5:05 8,6 10,2 0,5 0,5 5 1,66
5:15 8 9,7 0,6 0,5 5 1,83
5:25 7,6 9,2 0,4 0,5 5 2
5:30 7,2 8,7 0,4 0,5 5 2,08
5:35 6,8 8,4 0,4 0,3 3 2
5:40 6,4 8,1 0,4 0,3 3 2,25
5:45 6,1 7,9 0,3 0,2 2 2,33
5:50 5,8 7,6 0,3 0,3 3 2,42
5:55 5,3 7,3 0,5 0,3 3 2,5
5
6:00 5 7 0,3 0,3 3 2,58
6:05 4,8 6,8 0,2 0,2 2 2,67
6:10 4,4 6,6 0,4 0,2 2 2,75
6:15 4,2 6,3 0,2 0,3 3 2,83
6:20 3,8 6 0,4 0,3 3 2,92
6:25 3,5 5,7 0,3 0,3 3 3
Tabla 1 Mediciones de campo

Sin embargo, por recomendación de la norma ASTM D 3385 – 94, se emplea la

información recogida solamente por el tubo N°2 (anillo interior), las diferencias
calculadas se convierten a unidades de milímetros y a su vez el tiempo a horas
(columnas 7 y 8), para posteriormente realizar un gráfico que las relacione,
obteniendo así los siguientes puntos:

Ensayo de Infiltración
9
8
Infiltración (mm/hr)

7
6
5
4
3 Campo
2
1
0
0 0,5 1 1,5 2 2,5 3 3,5
Tiempo (hr)

Gráfico 1 Datos de campo

Ahora a través de los modelos empíricos aprendidos en el curso se busca obtener

una ecuación que se aproxime a los datos obtenidos en campo; a cada uno de estos
métodos se le realizó su respectivo cálculo de error, para esto se usó el método
cuadrático, medio que hace referencia a:
 𝑛

𝐸𝑟𝑟𝑜𝑟 = ∑(𝑖𝑚 − 𝑖𝑐 )2
𝑛=1

Donde
𝑖𝑚 = Infiltración calculada por el método respectivo
𝑖𝑐 = Infiltración medida en el ensayo.
El cálculo del error nos ayuda a obtener los coeficientes más óptimos, permitiendo
su utilización como parámetro en las iteraciones de coeficientes empíricos revisando
siempre que este se mantenga en su menor valor posible.
Modelo de Kostiakov

Este modelo parte de la ecuación: 𝑖(𝑡) = α𝑡 −𝛽 donde i representa la tasa de

infiltración en un tiempo t, α y β son valores empíricos que deben cumplir α(α>0) y
β(0< β<1), los cuales se obtienen de un proceso iterativo, posterior a la aplicación
de la ecuación se realiza el cálculo del error por el método cuadrático, para
relacionar este valor con el obtenido en campo y así obtener resultados más
precisos buscando estar entre unos márgenes aceptables, para esto se modifican
los valores de α,β hasta obtener el menor error posible (sumatoria de la última
columna cuyos valores cambian al variar los parámetros α y β).
La tabla representa los datos calculados con los parámetros a la izquierda,
finalmente se grafica esto (Kostiakov sobre el Grafico 1

Infiltración
Tiempo en
medida en el Kostiakov Error
Horas
ensayo (mm)
α 0,021916
0
β 0,4164431
8 0,25 8,743844 0,553304
Error 20,209279
6 0,5 6,551497 0,304149
7 0,75 5,533610 2,150299
7 1 4,908838 4,372959
4 1,16 4,614616 0,377753
4 1,33 4,359147 0,128987
4 1,5 4,146166 0,021365
5 1,66 3,974808 1,051018
5 1,83 3,816654 1,400308
5 2 3,678043 1,747570
5 2,08 3,618457 1,908661
3 2 3,557919 0,311274
3 2,25 3,501989 0,251993
 2 2,33 3,451405 2,106575
3 2,42 3,397359 0,157894
3 2,5 3,351655 0,123661
3 2,58 3,307977 0,094850
2 2,67 3,261077 1,590315
2 2,75 3,221229 1,491401
3 2,83 3,182991 0,033486
3 2,92 3,141762 0,020096
3 3 3,106597 0,011363
Tabla 2 Cálculos por Kostiakov

Ensayo de Infiltración
10
9
Infiltración (mm/hr)

8
7
6
5
4 Campo
3 Kosiakov
2
1
0
0 0,5 1 1,5 2 2,5 3 3,5
Tiempo (hr)

Gráfico 2 Modelo de Kostiakov con campo

Modelo de Horton

El modelo de Horton nos presenta la siguiente ecuación:

𝑖(𝑡) = 𝑖𝑓 + (𝑖0 − 𝑖𝑓 )𝑒 −𝛾𝑡

En este caso 𝑖𝑜 representa la infiltración inicial, e 𝑖𝑓 la infiltración final, los cuales
son estimados de la información de campo y por tanto de criterio propio para su
selección, la constante empírica 𝛾, al igual como se hizo en el caso anterior se itera,
y posteriormente se analiza por medio del error si es necesario realizar ajustes para
que este disminuya, hasta encontrar la información más precisa.
Es importante aclarar que para ninguno de estos modelos se está empleando la
ecuación para la infiltración acumulada. Análogamente a Kostiakov se representa la
tabla y se grafican los datos de Horton vs t, en el mismo grafico con los datos de
campo:
 Infiltración
Tiempo en
medida en el Horton Error
Horas
ensayo (mm) Io 10
0 If 2
8 0,25 8,554573 0,307551 γ 0,7971144
6 0,5 7,370303 1,877730 Error 16,499311
7 0,75 6,400005 0,359994
7 1 5,605019 1,945971
4 1,16 5,173353 1,376758
4 1,33 4,771197 0,594745
4 1,5 4,420006 0,176405
5 1,66 4,130234 0,756494
5 1,83 3,860271 1,298982
5 2 3,624521 1,891944
5 2,08 3,524160 2,178103
3 2 3,423177 0,179079
3 2,25 3,331005 0,109564
2 2,33 3,248777 1,559445
3 2,42 3,162328 0,026350
3 2,5 3,090521 0,008194
3 2,58 3,023150 0,000536
2 2,67 2,952320 0,906914
2 2,75 2,893487 0,798320
3 2,83 2,838289 0,026150
3 2,92 2,780257 0,048287
3 3 2,732053 0,071795
Tabla 3 Cálculos por Horton

Ensayo de Infiltración
9
8
Infiltración (mm/hr)

7
6
5
4 Campo
3
2 Horton
1
0
0 0,5 1 1,5 2 2,5 3 3,5
Tiempo (hr)

Gráfico 3 Modelo de Horton con campo

Modelo de Mezencev

𝑖(𝑡) = 𝑖𝑓 + α𝑡 −β
Para utilizar la ecuación de Mezencev se deben conocer parámetros como las
constantes empíricas α y β, el tiempo y la tasa de infiltración final en estado
permanente; para este caso tomamos la tasa de infiltración menor del Grafico 1,
correspondiente a 2 mm y se le asignaron valores a las constantes hasta que estas
dieran el menor rango de error; cabe resaltar que las constantes deben cumplir con
(α > 0) y ( 0 < β <1).

Infiltración
Tiempo en
medida en el Mecenzev Error
Horas
ensayo (mm) α 0,2158985
0 β 0,6539856
8 0,25 8,747296 0,558451 If 2
6 0,5 6,288056 0,082976 Error 23,206575
7 0,75 5,289268 2,926602
7 1 4,725155 5,174920
4 1,16 4,473071 0,223796
4 1,33 4,261486 0,068375
4 1,5 4,090403 0,008173
5 1,66 3,956337 1,089233
5 1,83 3,835490 1,356084
5 2 3,731896 1,608087
5 2,08 3,688038 1,721243
3 2 3,643900 0,414607
3 2,25 3,603500 0,364212
2 2,33 3,567277 2,456356
3 2,42 3,528908 0,279744
3 2,5 3,496732 0,246743
3 2,58 3,466216 0,217357
2 2,67 3,433702 2,055502
2 2,75 3,406287 1,977643
3 2,83 3,380160 0,144521
3 2,92 3,352189 0,124037
3 3 3,328498 0,107911
Tabla 4 Cálculos por Mecenzev
 Ensayo de Infiltración
10
9
8
Infiltración (mm/hr)
7
6
5
Campo
4
3 Mecenzev
2
1
0
0 0,5 1 1,5 2 2,5 3 3,5
Tiempo (hr)

Gráfico 4 Modelo de Mecenzev con campo

Método de Philip
1 −1
𝑓(𝑡) = 𝐴𝑡 2 + 𝐵
2
Philip tiene en cuenta que la intensidad de la precipitación es mayor que la
capacidad de infiltración del terreno. Para el uso de esta fórmula se deben conocer
parámetros A y B que dependen del tipo de terreno y el tiempo; como una de las
dificultades es el hallar correctamente estos parámetros se usó el método de
iteración hasta que el error fuese el menor.

Infiltración
Tiempo en
medida en el Philip Error
Horas
ensayo (mm)
A 8,29102785
0 B 0,70515739
8 0,25 8,996185 0,992385 Error 21,0524932
6 0,5 6,567799 0,322396
7 0,75 5,491985 2,274111
7 1 4,850671 4,619614
4 1,16 4,554170 0,307105
4 1,33 4,299774 0,089864
4 1,5 4,089955 0,008092
5 1,66 3,922700 1,160575
5 1,83 3,769609 1,513861
5 2 3,636478 1,859191
5 2,08 3,579554 2,017666
 3 2 3,521913 0,272393
3 2,25 3,468833 0,219805
2 2,33 3,420974 2,019167
3 2,42 3,369995 0,136896
3 2,5 3,327011 0,106936
3 2,58 3,286042 0,081820
2 2,67 3,242171 1,542988
2 2,75 3,204996 1,452016
3 2,83 3,169410 0,028700
3 2,92 3,131136 0,017197
3 3 3,098571 0,009716
Tabla 5 Cálculos por Philip

Ensayo de Infiltración
10
9
8
Infiltración (mm/hr)

7
6
5
Campo
4
3 Philip
2
1
0
0 0,5 1 1,5 2 2,5 3 3,5
Tiempo (hr)

Gráfico 5 Modelo de Phillip con campo

Finalmente, con el siguiente grafico se relacionan todas las ecuaciones en conjunto

con los datos de campo para observar más ampliamente la relación de proximidad
entre todos los métodos aplicados y los puntos de la información del ensayo.
Por tanto, observando los modelos en el grafico puede afirmarse que los métodos
de Kostiakov y Phillip presentan una mejor distribución con relación a los puntos
correspondientes a los valores de campo, además en ambos casos las
formulaciones requieren de un número más reducido de parámetros que en los
demás, lo cual es muy favorable para la precisión en los cálculos, razones por la
que se prefieren estos métodos sobre los demás para la determinación de la tasa
de infiltración.
 Ensayo de Infiltración
10

7
Infiltración (mm/hr)

6
Campo
5 Kosiakov
Horton
4
Mecenzev
3 Philip

0
0 0,5 1 1,5 2 2,5 3 3,5
Tiempo (hr)

Gráfico 6 Comparación entre modelos empíricos con datos de campo en Infiltración vs tiempo

7. CONCLUSIONES

• Es siempre muy favorable antes de realizar el ensayo, precisar que las

condiciones del terreno donde se pretende realizar sean las mejores
posibles, en búsqueda de obtener resultados más precisos; en este caso
en particular, en la zona se había realizado el mismo ensayo en horas
previas por lo que ya se encontraba el suelo con cierto grado de
saturación, lo que significa que los datos medidos pueden no representar
las condiciones verdaderas u optimas del terreno.

• Es importante la comparación entre métodos empíricos con la información

que puede obtenerse en una práctica de laboratorio como esta para
observar cómo puede variar la información entre sí y el alcance de cada
metodología teniendo en cuenta la limitantes, ya que en cada caso existe
alguna; igualmente es conveniente observar el número de parámetros y
su tipo, ya que estos hacen más precisa a una ecuación, por ejemplo en
la ecuación de Horton, maneja más parámetros y hace que su
funcionamiento sea inferior en muchas ocasiones.
 • Los métodos de iteración, pueden parecer que carecen de un verdadero
fundamento matemático en relación a estas aplicaciones, pero son
confiables gracias a los cálculos del error, que en este caso se pudieron
comparar con la información de las mediciones hechas, aunque por las
condiciones del suelo pueden no ser correctos, en un ambiente diferente
se pueden aplicar estas mismas formulaciones, y tener certeza de su
veracidad.
• El ensayo de infiltración, según fuentes diversas, tiene un gran número
de variaciones en su proceso de aplicación, hay una gran variedad de
formas de realizarlo, y no es necesario contar con gran cantidad de
instrumentos complejos; lo primordial es tener un buen control sobre la
información que se esta tomando y cuidar las condiciones del lugar, una
forma de estar seguro de estar realizando correctamente es notar un
patrón que llevan las mediciones, y estar pendiente de cuando este
cambia (disminuye o aumenta).

8. REFERENCIAS

• Alvarado B. Cesar, Barahona P. Marco; Comparación de tres métodos de

infiltración para calcular el balance hídrico del suelo, en la cuenca del rio
Suquiapa, El Salvador. San Salvador, 2016.
• Sayah Boumediene; Soluciones para la infiltración y redistribución en el
riego por superficie: aproximaciones simplificadas y modelización.
Departamento de ingeniería forestal, Escuela técnica superior de Ingeniería
Agronómica, Alimentaria y de Biosistemas; Madrid, 2017.
• Método de prueba estándar para la tasa de infiltración de suelos en campo
con doble anillo infiltrometro (ASTM D 3385 – 94)