Cálculo de la Infiltración según el Método de Kostiakov a partir del Ensayo de Doble Anillo

La determinación precisa de la tasa de infiltración del suelo es fundamental para el diseño de sistemas de riego, la evaluación de recarga de acuíferos y el modelado hidrológico de cuencas. Entre los métodos empíricos más utilizados destaca la ecuación de Kostiakov, desarrollada originalmente en 1932, que permite describir el comportamiento de la infiltración acumulada en función del tiempo mediante una relación potencial

En este artículo, explicaremos cómo aplicar el método de Kostiakov a partir de datos obtenidos con el infiltrómetro de doble anillo, siguiendo el procedimiento estandarizado ASTM D3385, y cómo interpretar los parámetros resultantes para la caracterización hidráulica del suelo.

El Ensayo de Doble Anillo: Fundamentos y Procedimiento

El infiltrómetro de doble anillo consiste en dos cilindros metálicos concéntricos que se hincan en el suelo. El anillo exterior actúa como barrera hidráulica para minimizar el flujo lateral, mientras que en el anillo interior se mide la disminución del nivel de agua para calcular la infiltración vertical

Especificaciones típicas del equipo:

• Anillo interior: diámetro 0,28–0,32 m, altura ~0,40 m
• Anillo exterior: diámetro 0,53–0,57 m, altura ~0,30 m

Procedimiento resumido (ASTM D3385):

1. Instalación de los anillos con mínima alteración del suelo.
2. Saturación previa del perfil (opcional, según objetivo).
3. Llenado simultáneo de ambos anillos y registro del nivel de agua en el interior a intervalos de tiempo definidos.
4. Cálculo de la lámina infiltrada y la velocidad instantánea.

💡 Nota técnica: La tasa de infiltración medida en condiciones de saturación tiende asintóticamente a la conductividad hidráulica saturada (Ks) del suelo

La Ecuación de Kostiakov: Formulación Matemática

El modelo empírico de Kostiakov expresa la velocidad de infiltración instantánea (I) como una función potencial del tiempo (t):

$$I(t)=K\cdot t^n$$I(t)=K⋅tn

Donde:

La infiltración acumulada (Z) se obtiene integrando la ecuación anterior:

$$Z(t)=\frac{K}{n+1}\cdot t^{n+1}\text{(para }n\mathrm{\ne }-1\text{)}$$Z(t)=n+1K​⋅tn+1(para n=−1)

📊 Estimación de parámetros por regresión log-log

Para ajustar los parámetros K y n a partir de datos de campo, se linealiza la ecuación aplicando logaritmos:

$$\log (I)=\log (K)+n\cdot \log (t)$$log(I)=log(K)+n⋅log(t)

Esto transforma el problema en una regresión lineal simple donde:

• Pendiente = n
• Ordenada al origen = $\log (K)$→ $K={10}^{\text{intercepto}}$

🎯 Clasificación de la tasa básica de infiltración

Una vez estimada la velocidad de infiltración en un tiempo de referencia (t_base), se clasifica el suelo según criterios hidrológicos ampliamente aceptados:

Esta clasificación orienta decisiones de diseño: desde la selección de emisores en riego por goteo hasta la evaluación de riesgo de escorrentía superficial.

🔹 Bibliografía Recomendada

1. Kostiakov, A.N. (1932). On the dynamics of the coefficient of water percolation in soils. 10th International Congress of Soil Science.
2. ASTM D3385-18. Standard Test Method for Infiltration Rate of Soils in Field Using Double-Ring Infiltrometer. ASTM International 分析测试百科网.
3. Ahuja, L.R. et al. (2007). Scaling Parameters of the Lewis-Kostiakov Water Infiltration Equation. ASABE. www.ars.usda.gov
4. Smith, R.E. (1972). The infiltration envelope: results from a theoretical infiltrometer. Journal of Hydrology.
5. US Salinity Laboratory Staff. (1954). Diagnosis and Improvement of Saline and Alkali Soils. USDA Handbook 60.