Conductividad Hidráulica en Un Suelo No Saturado

UNIVERSIDAD AUTÓNOMA DE QUERÉTARO
FACULTAD DE INGENIERÍA
POSGRADO
Unsaturated Soil Mechanics
Dr. Eduardo Rojas González
Problema
Determinar la curva de permeabilidad a partir de la siguiente curva característica:
En base a la siguiente ecuación:
𝑚
𝑘𝑠
𝑘𝑤 (𝜃𝑤 ) = 𝐴 ∑{(2𝑗 + 1 − 2𝑖)(𝑢𝑎 − 𝑢𝑤 )−2
𝑗 } , 𝑖 = 1, 2, … , 𝑚
𝑘𝑤 𝑖
𝑖=1
Dónde:
𝑘𝑤 (𝜃𝑤 ) Coeficiente de permeabilidad para un contenido volumétrico de agua dado, (𝜃𝑤 )𝑖 ,
correspondiente al i-ésimo intervalo en m/s.
𝑖 Número de intervalos correspondientes a un decremento del contenido de agua. Por
ejemplo, i = 1 es el primer intervalo para el cual está dado un contenido volumétrico
de agua 𝜃𝑠 ; 𝑖 = 𝑚 es el último intervalo correspondiente con el menor contenido
volumétrico de agua en la curva característica experimental, 𝜃𝐿 .
𝑗 Es un contador desde “i” hasta “m”
𝑚 Número total de intervalos entre el contenido volumétrico de agua para el caso
saturado, 𝜃𝑠 , y el menor contenido volumétrico de agua obtenido experimentalmente
de la curva característica, 𝜃𝐿 .
Ing. Rafael Ortiz Hernández 1

POSGRADO
Solución
Los valores de esta curva SWRC se determinaron al colocar la gráfica en EXCEL y aproximarlos a
la unidad más cercana en el eje horizontal (succión) y en el eje vertical (contenido de agua
volumétrica) al decimal más cercano (0.1%). Los resultados se presentan en la Tabla 1 y su grafica
se muestra en la Figura 2.
Tabla 1. Datos de la curva SWRC obtenidos de análisis gráfico de EXCEL.
Dato 𝑈𝑎 − 𝑈𝑤 (kPa) 𝜃 (%)

1 10 38.7
2 76 37.7
3 130 33.6
4 180 31.0
5 243 28.5
6 300 26.0
7 381 23.3
8 485 20.6
9 660 18.2
10 960 15.6
11 1670 12.9
12 2830 10.5
13 4730 7.8
14 16500 5.2
Considerando que se tienen 14 datos, el valor “m” = 14.
Figura 1. Datos de la Tabla 1 (en verde) sobrepuestos sobre la gráfica original.

POSGRADO
Se proceden a obtener los valores de permeabilidad saturada teóricas 𝐾𝑠𝑎𝑡(𝑡) de acuerdo a la formula:
𝑚
(2𝑗 + 1 − 2𝑖)
𝐾𝑠𝑎𝑡(𝑡) = 𝐴𝑑 ∑
𝑆𝑗2
𝑗=1
𝑇𝑠2 𝛾𝑤
𝐴𝑑 = Δ𝜃
2𝜈
Donde:
𝑇𝑠 Fuerza tensional interface aire-agua, constante igual a 0.0728 N/m (Wolfram|Alpha, 2022),
𝛾𝑤 Peso volumétrico de agua, a 20° C se tiene 9.789 kN/m3 (Wolfram|Alpha, 2022),,
𝜈 Viscosidad absoluta del agua, a 20° C se tiene como constante de 0.001 N∙s/m2,
𝜃
Δ𝜃 Intervalo de contenido volumétrico de agua, considerado como 𝑛𝑠 , siendo el 𝜃𝑠 el valor
máximo de contenido volumétrico de agua y 𝑛 el número de datos disponibles:
38.7
Δ𝜃 = = 2.764
14
𝑆𝑗 Succión obtenida en el intervalo j, convertido a N/m2 (1 kPa = 1000 N/m2).
La tabla 2 muestra el cálculo del factor (2j+1-2i) para cada intervalo de j en los 14 puntos de medición.
La tabla 3 muestra los valores de succión a usar en cada
Tabla 2. Calculo del factor 2𝑗 + 1 − 2𝑖
i/j 1 2 3 4 5 6 7 8 9 10 11 12 13 14
1 1
2 3 1
3 5 3 1
4 7 5 3 1
5 9 7 5 3 1
6 11 9 7 5 3 1
7 13 11 9 7 5 3 1
8 15 13 11 9 7 5 3 1
9 17 15 13 11 9 7 5 3 1
10 19 17 15 13 11 9 7 5 3 1
11 21 19 17 15 13 11 9 7 5 3 1
12 23 21 19 17 15 13 11 9 7 5 3 1
13 25 23 21 19 17 15 13 11 9 7 5 3 1
14 27 25 23 21 19 17 15 13 11 9 7 5 3 1

POSGRADO
Tabla 3. Definición de cada valor de succión (kPa) por cada intervalo de la curva.
i/j 1 2 3 4 5 6 7 8 9 10 11 12 13 14
1 10000
2 76000 76000
3 130000 130000 130000
4 180000 180000 180000 180000
5 243000 243000 243000 243000 243000
6 300000 300000 300000 300000 300000 300000
7 381000 381000 381000 381000 381000 381000 381000
8 485000 485000 485000 485000 485000 485000 485000 485000
9 660000 660000 660000 660000 660000 660000 660000 660000 660000
10 960000 960000 960000 960000 960000 960000 960000 960000 960000 960000
11 1670000 1670000 1670000 1670000 1670000 1670000 1670000 1670000 1670000 1670000 1670000
12 2830000 2830000 2830000 2830000 2830000 2830000 2830000 2830000 2830000 2830000 2830000 2830000
13 4730000 4730000 4730000 4730000 4730000 4730000 4730000 4730000 4730000 4730000 4730000 4730000 4730000
14 16500000 16500000 16500000 16500000 16500000 16500000 16500000 16500000 16500000 16500000 16500000 16500000 16500000 16500000
(2𝑗+1−2𝑖)
Tabla 4. Sumatoria de término ∑𝑚
𝑗=1 𝑆𝑗2
i/j 1 2 3 4 5 6 7 8 9 10 11 12 13 14
1 1.000E-08
2 5.194E-10 1.731E-10
3 2.959E-10 1.775E-10 5.917E-11
4 2.160E-10 1.543E-10 9.259E-11 3.086E-11
5 1.524E-10 1.185E-10 8.468E-11 5.081E-11 1.694E-11
6 1.222E-10 1.000E-10 7.778E-11 5.556E-11 3.333E-11 1.111E-11
7 8.956E-11 7.578E-11 6.200E-11 4.822E-11 3.444E-11 2.067E-11 6.889E-12
8 6.377E-11 5.527E-11 4.676E-11 3.826E-11 2.976E-11 2.126E-11 1.275E-11 4.251E-12
9 3.903E-11 3.444E-11 2.984E-11 2.525E-11 2.066E-11 1.607E-11 1.148E-11 6.887E-12 2.296E-12
10 2.062E-11 1.845E-11 1.628E-11 1.411E-11 1.194E-11 9.766E-12 7.595E-12 5.425E-12 3.255E-12 1.085E-12
11 7.530E-12 6.813E-12 6.096E-12 5.378E-12 4.661E-12 3.944E-12 3.227E-12 2.510E-12 1.793E-12 1.076E-12 3.586E-13
12 2.872E-12 2.622E-12 2.372E-12 2.123E-12 1.873E-12 1.623E-12 1.373E-12 1.124E-12 8.740E-13 6.243E-13 3.746E-13 1.249E-13
13 1.117E-12 1.028E-12 9.386E-13 8.492E-13 7.598E-13 6.705E-13 5.811E-13 4.917E-13 4.023E-13 3.129E-13 2.235E-13 1.341E-13 4.470E-14
14 9.917E-14 9.183E-14 8.448E-14 7.713E-14 6.979E-14 6.244E-14 5.510E-14 4.775E-14 4.040E-14 3.306E-14 2.571E-14 1.837E-14 1.102E-14 3.673E-15
Σ 1.153E-08 9.180E-10 4.786E-10 2.715E-10 1.544E-10 8.517E-11 4.395E-11 2.074E-11 8.660E-12 3.131E-12 9.823E-13 2.773E-13 5.572E-14 3.673E-15
Se procede a calcular la constante 𝐴𝑑 :
𝑇𝑠2 𝛾𝑤
𝐴𝑑 = Δ𝜃
2𝜈
N 2 𝑁 𝑁3
(0.0728 ) (9.789 × 103 3 ) 1.43397 5 2
𝐴𝑑 =
m 𝑚 (2.764 %) = 𝑚 = 716.98 𝑁
𝑁∙𝑠 𝑁∙𝑠 𝑚3 ∙ 𝑠
2 (0.001 2 ) 0.002 2
𝑚 𝑚

POSGRADO
Con esta constante se procede a calcular la permeabilidad de cada intervalo, y se considera la

permeabilidad en el primer punto como la permeabilidad experimental de saturación, 𝐾𝑠𝑎𝑡(𝑒) .
13
(2𝑗 + 1 − 2𝑖) 𝑁2 −8
𝑚4 𝑚
𝑘𝑠𝑎𝑡(𝑒) = 𝐴𝑑 ∑ 2 = 716.98 3 ∙ 1.153 × 10 2
= 8.267 × 10−6
𝑆𝑗 𝑚 ∙ 𝑠 𝑁 𝑠
𝑗=1
𝑐𝑚
𝑘𝑠𝑎𝑡(𝑒) = 8.267 × 10−4
𝑠
Se propone un 𝑘𝑠 = 10-7 cm/s (10-9 m/s) debido a que el rango de succiones y el contenido de agua
volumétrico es similar al de un suelo arcilloso. La permeabilidad de una arcilla en estado saturado es
de 10-8 a 10-11 m/s (Briaud, 2013).
Figura 2. Curvas SWRC para distintos materiales (Briaud, 2013) y la curva SWRC del problema.
La relación entre permeabilidad experimental y teórica se tiene como:

𝑐𝑚
𝑘𝑠𝑎𝑡(𝑒) 1 × 10−7
= 𝑠 = 1.21 × 10−4
𝑘𝑠𝑎𝑡(𝑡) 8.267 × 10−4 𝑐𝑚
𝑠
Se procede a realizar el cálculo para cada grado de succión:
Tabla 5. Resultados del cálculo de permeabilidad de acuerdo al grado de saturación
θ(%) k(θ) (cm/s) θ(%) k(θ) (cm/s)

38.7 1.000E-09 20.6 1.798E-12
37.7 7.961E-11 18.2 7.511E-13
33.6 4.151E-11 15.6 2.715E-13
31.0 2.355E-11 12.9 8.520E-14
28.5 1.339E-11 10.5 2.405E-14
26.0 7.386E-12 7.8 4.832E-15
23.3 3.812E-12 5.2 3.186E-16

POSGRADO
Permeabilidad de acuerdo a Fredlund & Xing

Log (Succión (kPa))
10 100 1,000 10,000 100,000
1.E+00
1.E-01
1.E-02
1.E-03
1.E-04
1.E-05
1.E-06
Permeabilidad (cm/s)
1.E-07
1.E-08
1.E-09
1.E-10
1.E-11
1.E-12
1.E-13
1.E-14
1.E-15
1.E-16
Figura 3. Gráfica de Permeabilidad de acuerdo a la succión con la fórmula de Fredlund & Xing (1994)

POSGRADO
Modelo de SWRC de Van Genuchten
La ecuación de Van Genuchten es una relación del contenido de agua volumétrico normalizada, Θ,
y la succión, 𝜓:
𝑚
𝜃 − 𝜃𝑟 1
Θ=( )=[ ]
𝜃𝑠 − 𝜃𝑟 1 + (𝑝𝜓)𝑛
Donde:
𝜃 Contenido de agua volumétrico en un punto de la curva,
𝜃𝑠 Contenido de agua volumétrico al punto de saturación, 𝜃𝑠 = 38.7%,
𝜃𝑟 Contenido de agua volumétrico residual, 𝜃𝑟 = 5.2%,

1
𝑝, 𝑛, 𝑚 Parámetros de ajuste para la curva, 𝑚 = 1 − 𝑛.n
Utilizando EXCEL se selecciona por medio gráfico los parámetros p, n y m que generen una curva
similar a la original, considerando reducir el error relativo lo más posible.
Comparativa de SWRC experimental vs. SWRC de Van Genuchten

1.00
0.90
0.80
0.70
0.60
0.50
0.40
0.30 Experimental
Van Genuchten
0.20
p = 0.00436
0.10 n = 1.9
m = 0.47368
0.00
1 10 100 1,000 10,000 100,000
Figura 4. Comparativa de SWRC experimental con la SWRC de Van Genuchten

POSGRADO
El error relativo se considera como:
|𝜃𝑒𝑥𝑝 − 𝜃𝑉𝑎𝑛 𝐺𝑒𝑛𝑢𝑐ℎ𝑡𝑒𝑛 |

𝐸𝑟𝑟𝑜𝑟 𝑟𝑒𝑙𝑎𝑡𝑖𝑣𝑜 (%) = × 100
𝜃𝑒𝑥𝑝
En este caso para la curva obtenida tenemos:
Tabla 6. Resultados del cálculo de la SWRC con el modelo de Van Genuchten y su error relativo.
θexp ψexp θVan Genuchten Error relativo

Θ ΘVan Genuchten
(%) (kPa) (%) (%)
38.7 1.00 10 1.00 38.7 0.11%
37.7 0.97 76 0.93 36.9 2.09%
33.6 0.85 130 0.85 34.4 2.27%
31.0 0.77 180 0.78 31.8 2.49%
28.5 0.70 243 0.69 28.7 0.66%
26.0 0.62 300 0.62 26.3 0.98%
23.3 0.54 381 0.55 23.4 0.46%
20.6 0.46 485 0.48 20.6 0.07%
18.2 0.39 660 0.39 17.4 4.45%
15.6 0.31 960 0.30 14.2 9.23%
12.9 0.23 1670 0.20 10.8 16.66%
10.5 0.16 2830 0.14 8.7 17.36%
7.8 0.08 4730 0.09 7.4 5.19%
5.2 0.00 16500 0.04 5.9 13.73%
Se observa que el modelo de Van Genuchten tiene resultados muy cercanos a los experimentales en
la rama descendente de la curva experimental, pero presenta errores en el valor de contenido de agua
volumétrico a valores de succión de 12.9 kPa, 10.5 kPa y el valor residual a 5.2 kPa.

POSGRADO
Ahora se realiza el cálculo de la permeabilidad utilizando la fórmula de Van Genuchten, que

determina la permeabilidad con un factor utilizando los parámetros de ajuste de la curva SWRC:
𝑘 = 𝑘𝑠 ∙ 𝑘𝑟
[1 − (𝑝𝜓)𝑛−1 (1 + (𝑝𝜓)𝑛 )−𝑚 ]2

𝑘𝑟 = 𝑚
[1 + (𝑝𝜓)𝑛 ] 2
Recordando 𝑘𝑠 = 10-9 cm/s:
Tabla 7. Resultados de cálculo de permeabilidad utilizando la función de Van Genuchten.
ψexp θVan Genuchten 𝒌

ΘVan Genuchten 𝒌𝒓
(kPa) (%) (cm/s)
10 1.00 38.6 0.883870972 8.839E-10
76 0.93 36.5 0.410667673 4.107E-10
130 0.85 33.8 0.212962616 2.130E-10
180 0.78 31.2 0.116874017 1.169E-10
243 0.69 28.3 0.057230751 5.723E-11
300 0.62 26.1 0.031612583 3.161E-11
381 0.55 23.6 0.014837486 1.484E-11
485 0.48 21.1 0.006412599 6.413E-12
660 0.39 18.2 0.002023425 2.023E-12
960 0.30 15.3 0.000458150 4.581E-13
1670 0.20 11.9 0.000046742 4.674E-14
2830 0.14 9.8 0.000005091 5.091E-15
4730 0.09 8.3 0.000000579 5.789E-16
16500 0.04 6.4 0.000000003 2.875E-18

POSGRADO
Ahora se comparan los valores de permeabilidad obtenidos de ambos métodos, en forma de tabla y
de forma gráfica:
Tabla 8. Tabla comparativa de valores de permeabilidad entre los métodos de Fredlung & Xing y Van Genuchten.
ψexp 𝒌𝑭𝒓𝒆𝒅𝒍𝒖𝒏𝒅 & 𝑿𝒊𝒏𝒈 𝒌𝑽𝒂𝒏 𝑮𝒆𝒏𝒖𝒄𝒉𝒕𝒆𝒏

(kPa) (cm/s) (cm/s)
10 1.000E-09 8.839E-10
76 7.961E-11 4.107E-10
130 4.151E-11 2.130E-10
180 2.355E-11 1.169E-10
243 1.339E-11 5.723E-11
300 7.386E-12 3.161E-11
381 3.812E-12 1.484E-11
485 1.798E-12 6.413E-12
660 7.511E-13 2.023E-12
960 2.715E-13 4.581E-13
1670 8.520E-14 4.674E-14
2830 2.405E-14 5.091E-15
4730 4.832E-15 5.789E-16
16500 3.186E-16 2.875E-18
Comparativa de valores de Permeabilidad obtenidos con diferentes métodos

Log (Succión (kPa))
10 100 1,000 10,000 100,000
1.E+00
1.E-01
1.E-02
1.E-03
1.E-04
1.E-05
Permeabilidad (cm/s)
1.E-06
1.E-07
1.E-08
1.E-09
1.E-10
1.E-11
1.E-12
1.E-13
1.E-14 Fredlund & Xing
1.E-15
Van Genuchten
1.E-16
1.E-17
1.E-18
Figura 5. Comparativa de SWRC experimental con la SWRC de Van Genuchten

POSGRADO
Conclusión
Las fórmulas de Fredlund & Xing y la de Van Genuchten para determinar el valor de permeabilidad
tienen una tendencia de curva similar pero están separados casi por 1 magnitud de escala (1 a 4 veces),
por lo que se debe de tener cuidado al hacer comparativas de resultados en escalas logarítmicas.
Ambos entran en concordancia con un material impermeable, como lo es una arcilla de grano muy
fino (bentonita, alguna arcilla del grupo esmectita) (Shirazi, 2010).
La explicación de estas diferencias de valores es de como los autores de las ecuaciones formularon
las variables de ajuste de curva, mientras Van Genuchten propone un modelo similar al
comportamiento de una hipérbola, Fredlund & Xing usan un modelo que usa logaritmos, exponentes
y potencias para ajustar su curva.
Referencias
Briaud, J. L. (2013). Geotechnical engineering: unsaturated and saturated soils. John Wiley & Sons.
Fredlund, D. G., Xing, A., & Huang, S. (1994). Predicting the permeability function for unsaturated
soils using the soil-water characteristic curve. Canadian Geotechnical Journal, 31(4), 533-546.
Swindells, J. F., Coe Jr, J. R., & Godfrey, T. B. (1952). Absolute viscosity of water at 20 C. Journal
of Research of the National Bureau of Standard, 48(1).
https://www.wolframalpha.com/input?i=density+of+water+at+20%C2%B0+C
https://www.wolframalpha.com/input?i=surface+tension+of+water
https://www.wolframalpha.com/input?i=absolute+viscosity+of+water+at+20%C2%B0C
Shirazi, S. M., Kazama, H., Salman, F. A., Othman, F., & Akib, S. (2010). Permeability and swelling
characteristics of bentonite. International Journal of Physical Sciences, 5(11), 1647-1659.

Conductividad Hidráulica en Un Suelo No Saturado

Cargado por

Copyright:

Formatos disponibles

Conductividad Hidráulica en Un Suelo No Saturado

Cargado por

Información del documento

Título original

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Conductividad Hidráulica en Un Suelo No Saturado

Cargado por

Copyright:

Formatos disponibles

UNIVERSIDAD AUTÓNOMA DE QUERÉTARO

Determinar la curva de permeabilidad a partir de la siguiente curva característica:

En base a la siguiente ecuación:

Ing. Rafael Ortiz Hernández 1

Tabla 1. Datos de la curva SWRC obtenidos de análisis gráfico de EXCEL.

Dato 𝑈𝑎 − 𝑈𝑤 (kPa) 𝜃 (%)

Considerando que se tienen 14 datos, el valor “m” = 14.

Figura 1. Datos de la Tabla 1 (en verde) sobrepuestos sobre la gráfica original.

Ing. Rafael Ortiz Hernández 2

𝛾𝑤 Peso volumétrico de agua, a 20° C se tiene 9.789 kN/m3 (Wolfram|Alpha, 2022),,

𝑆𝑗 Succión obtenida en el intervalo j, convertido a N/m2 (1 kPa = 1000 N/m2).

Tabla 2. Calculo del factor 2𝑗 + 1 − 2𝑖

Ing. Rafael Ortiz Hernández 3

Se procede a calcular la constante 𝐴𝑑 :

Ing. Rafael Ortiz Hernández 4

Con esta constante se procede a calcular la permeabilidad de cada intervalo, y se considera la

La relación entre permeabilidad experimental y teórica se tiene como:

Se procede a realizar el cálculo para cada grado de succión:

Tabla 5. Resultados del cálculo de permeabilidad de acuerdo al grado de saturación

θ(%) k(θ) (cm/s) θ(%) k(θ) (cm/s)

Ing. Rafael Ortiz Hernández 5

Permeabilidad de acuerdo a Fredlund & Xing

Ing. Rafael Ortiz Hernández 6

Modelo de SWRC de Van Genuchten

𝜃 Contenido de agua volumétrico en un punto de la curva,

𝜃𝑠 Contenido de agua volumétrico al punto de saturación, 𝜃𝑠 = 38.7%,

𝜃𝑟 Contenido de agua volumétrico residual, 𝜃𝑟 = 5.2%,

Comparativa de SWRC experimental vs. SWRC de Van Genuchten

Ing. Rafael Ortiz Hernández 7

El error relativo se considera como:

|𝜃𝑒𝑥𝑝 − 𝜃𝑉𝑎𝑛 𝐺𝑒𝑛𝑢𝑐ℎ𝑡𝑒𝑛 |

En este caso para la curva obtenida tenemos:

θexp ψexp θVan Genuchten Error relativo

Ing. Rafael Ortiz Hernández 8

Ahora se realiza el cálculo de la permeabilidad utilizando la fórmula de Van Genuchten, que

[1 − (𝑝𝜓)𝑛−1 (1 + (𝑝𝜓)𝑛 )−𝑚 ]2

Recordando 𝑘𝑠 = 10-9 cm/s:

Tabla 7. Resultados de cálculo de permeabilidad utilizando la función de Van Genuchten.

ψexp θVan Genuchten 𝒌

Ing. Rafael Ortiz Hernández 9

ψexp 𝒌𝑭𝒓𝒆𝒅𝒍𝒖𝒏𝒅 & 𝑿𝒊𝒏𝒈 𝒌𝑽𝒂𝒏 𝑮𝒆𝒏𝒖𝒄𝒉𝒕𝒆𝒏

Comparativa de valores de Permeabilidad obtenidos con diferentes métodos

Ing. Rafael Ortiz Hernández 10

Ing. Rafael Ortiz Hernández 11

También podría gustarte