LAS TRINCHERAS : Son unas de las obras de captación horizontal de agua subterránea.

Se emplean muy
eficazmente en zonas costeras debido a que el agua dulce flota por encima del agua salada y se puede extraer
con relativa facilidad. Hay peligro de contaminación si se deja expuesta. También pueden emplearse drenes
horizonales. 2
ln⁡(𝑥/𝑟) (𝑙𝑛(𝑥/𝑟)) 𝑥−𝑟 1
−( ) + √[ + 4( ) ( 2 2 ) 𝑦]
𝐾𝐷 𝜋 𝐾𝐷2 𝜋 2 𝑥𝑟 𝜋 𝐾𝑇
𝑞1 =
𝑥−𝑟 1
2 ( 𝑥𝑟 ) ( 2 2 )
𝜋 𝐾𝑇

Recuérdese que el 𝑞𝑡 es por metro de trinchera. Esta

ecuación es un caso particular de la general, cuando uno
de los puntos coincide con la cara de la trinchera.

a) Para régimen permanente, lineal y flujo unidireccional

en acuífero libre.

Cuando 𝐼 < 𝐼𝐶𝑅𝐷 el régimen de circulación es lineal

(darciano) y se modela con la parábola de Dupuit luego:
𝐾𝐷 2
𝑞= (ℎ − ℎ𝑔2 )
2𝑥𝑜 𝑜
Y se sabe que 𝑞𝑡 = 2𝑞 es el caudal total por metro de
trinchera.

b) Para régimen permanente, no lineal y flujo

unidireccional en acuífero libre.

a) Para régimen permanente, no lineal y flujo

unidireccional en acuífero confinado.
1 1 4𝑆𝑔
− + +
𝐾𝐷 𝑚 √𝐾𝐷2 𝑚2 𝐾𝑇2 𝑚2 𝑥𝑜
𝑞=
2/𝐾𝑇2 𝑚2
b) Para régimen permanente, lineal y flujo
unidireccional en acuífero confinado.
𝑆𝑔 𝐾𝐷 𝑚
𝑞=
𝑥𝑜
“La hidrología del agua subterránea puede ser definida como la ciencia de la ocurrencia, distribución y
movimiento del agua por debajo de la superficie de la tierra.” Todd & Mays (2005)
Prop hidrogeológicas relacionadas con la conducción del agua.
Hidráulica de las captaciones del agua subterránea.
Pérez F (2001). La explotación del agua subterránea, un nuevo enfoque.
Martínez-Santos, P. (2018). Hidrogeología: principios y aplicaciones. McGraw Hill.
Fitts, C. R. (2012). Groundwater science, Second Edition. Elsevier.
Bear, J. (2012). Hydraulics of groundwater, Dover Edition. McGraw Hill.
Homogeneidad

Cuando las propiedades del acuífero son iguales en cualquier punto o

lugar de este, es homogéneo. Cuando cambian en el espacio se dice que
es heterogéneo.
Isotropía
Cuando las propiedades físicas en un punto determinado en un medio se
mantienen iguales en todas las direcciones se dice que el medio es isótropo
o isotrópico.
Conclusión importante
De manera estricta, todo acuífero es heterogéneo y anisotrópico.
En algunos casos es posible transformar:
a) Un medio heterogéneo en otro homogéneo equivalente.
b) Un medio anisótropo en otro equivalente isótropo.
Flujo paralelo a los estratos.
La conductividad 𝑲𝑫𝒆𝒒 representa la conductividad
de un estrato homogéneo de espesor saturado ∑ 𝒎𝒊 .
Las sumatorias se aplican desde 𝒊 = 𝟏 hasta 𝒏
estratos.
Flujo perpendicular a los estratos

La caída total de la carga hidráulica está dada por:

𝚫𝒉 = 𝒉𝑨 − 𝒉𝑩 = 𝚫𝒉𝟏 + 𝚫𝒉𝟐 + 𝚫𝒉𝟑 𝑊 ∑𝑛𝑖=1 𝐾𝐷𝑖 𝑚𝑖

𝐾𝐷𝑒𝑞 = 𝑚𝑖 𝐾𝐷𝑒𝑞 =
∑𝑛𝑖=1
𝐾𝐷𝑖
∑𝑛𝑖=1 𝑚𝑖
Como: 𝑰𝒏 = 0.002 < 𝐼𝐶𝑅𝐷 = 0.032 el régimen de flujo en el estrato 1 es lineal, luego se cumple la Ley de Darcy.
∆ℎ
𝐼𝑛 =
𝐿 𝑈
Flujo natural permanente en acuífero confinado. 𝑄 = 𝑈𝐴 ℎ=−
𝐾𝐷
𝑥 + ℎ𝑜
a) Para régimen permanente, lineal y flujo unidireccional en acuífero confinado.
a) Para régimen permanente, no lineal y flujo unidireccional en acuífero confinado. ℎ = − ( 1 𝑈 + 1 𝑈 2 ) 𝑥 + ℎ
2 𝑜
𝐾𝐷 𝐾𝑇
Flujo natural permanente en acuífero libre.
Ley de Darcy, el caudal por unidad de ancho de una sección vertical
saturada de altura ℎ como se ve en la siguiente imagen:
a) Para régimen permanente, lineal y flujo unidireccional en acuífero libre.
a) Que, para pequeñas inclinaciones de la superficie libre, las líneas de
corriente eran horizontales (las líneas equipotenciales serán entonces
prácticamente verticales). 2𝑞𝑥
ℎ2 = ℎ𝑜2 −
𝐾𝐷
b) Que el gradiente hidráulico era igual a la pendiente de la superficie
freática (tangente en lugar del seno del ángulo de la línea de corriente) y
que no variaba con la profundidad.

Red de flujo. y mapas de hidroisohipsas. 𝐾𝐷 2

𝑞= (ℎ − ℎ 2 )
2𝑥 𝑜
La existencia en el régimen lineal de un Potencial Velocidad (𝑥, 𝑦, 𝑧, 𝑡)
implica que el flujo es irrotacional y por tanto existe la Función
Corriente  (𝑥, 𝑦, 𝑧, 𝑡) que cumple con la condición de ser ortogonal
con la anterior, dando lugar a la existencia de una red de flujo que en
el caso del flujo bidimensional resultaría:
Los mapas de la superficie freática/libre (acuíferos libres) o de la
superficie potenciométrica/piezométrica (acuíferos confinados) son
una excelente herramienta para comprender y estimar los patrones del
flujo del agua subterránea en acuíferos.
La superficie del agua subterránea (freática o potenciométrica) es el
lugar geométrico de los puntos que corresponden a la altura del agua
que representa la carga hidráulica en cada punto.
Estas superficies se representan en un plano horizontal por curvas que
unen puntos de igual cota del agua subterránea.
Términos equivalentes: isopiésticas,
isopiezas, hidroisohipsas, líneas
equipotenciales, mapa de contorno,
mapa con curvas de nivel del agua
subterránea.
Estructura y elementos de los pozos de bombeo.

Modelo de flujo hacia un pozo de bombeo. Zonas y límites.

Si se supone que la superficie piezométrica del acuífero es

horizontal, ocurrirá que, en cada momento, el cono de
abatimiento abarcará el área del acuífero que haga falta para
suministrar el caudal extraído. Esta área se conoce como área
de influencia y la distancia desde el centro del pozo de bombeo
hasta sus límites externos, como radio de influencia del pozo

Se ha demostrado que en el flujo hacia un pozo se producen

desviaciones importantes de la ley de Darcy y el flujo no lineal
se presenta en zonas más o menos alejadas del pozo de bombeo
y en acuíferos de alta conductividad hidráulica (granulometría
gruesa y rocas fisuradas, fracturadas o con canales de
disolución). El flujo no lineal también se presenta en acuíferos
de baja conductividad hidráulica, en zonas alejadas al pozo de bombeo por lo
que el análisis del flujo hacia los pozos se debe hacer siempre partiendo del
enfoque no lineal.

De este modo pueden aparecer alrededor del pozo de bombeo los distintos
regímenes de circulación del agua subterránea (desde el darciano hasta el
turbulento puro)

¿Cómo determinar las zonas en las que ocurren los diferentes tipos de flujo?
𝑄 𝑇𝐷 2 𝑄 𝑇𝐷 2
𝑟𝑡 = ( ) 𝑟𝐷 = ( )
38𝜋𝑇𝐷 𝑇𝑇 0.1𝜋𝑇𝐷 𝑇𝑇
Radio turbulento Radio de Darcy
Zona de Régimen Turbulento puro: Ubicada entre 𝑟𝑝 y 𝑟𝑡 .
Zona de Régimen no lineal: Ubicada entre 𝑟𝑝 y 𝑟𝐷 .
Zona de Régimen lineal o darciano: Ubicada entre 𝑟𝐷 y 𝑟𝑜 .

Si 𝑟𝐷 < 𝑟𝑝 el flujo es darciano en todo el campo y no existen las otras dos zonas.

Si 𝑟𝐷 > 𝑟𝑝 y 𝑟𝑡 < 𝑟𝑝 existen solamente dos zonas: no lineal alrededor del pozo,

lineal más allá del radio 𝑟𝐷 .

Si 𝑟𝑡 > 𝑟𝑝 teóricamente pueden aparecer las tres zonas, pero en la práctica, la

zona de flujo darciano es posible que no exista en realidad.
Flujo lineal permanente hacia pozos en acuíferos confinados.

𝑸 𝒓𝒐
𝒉𝟐𝒐 − 𝒉𝟐 = 𝒍𝒏 ( )
𝝅𝑲𝑫 𝒓
𝑄 𝑟𝑜
𝑆= 𝑙𝑛 ( )
2𝜋𝐾𝐷 (ℎ𝑜 − 𝑆/2) 𝑟
𝑄 𝑟𝑜
𝑆 = 0.366 𝑙𝑜𝑔 ( )
𝐾𝐷 (ℎ𝑜 − 𝑆/2) 𝑟
𝑄 𝑟𝑜
𝑆 = 0.366 𝑙𝑛 ( )
ℎ𝑜 𝐾𝐷 𝑟

𝑄 𝑄
∆𝑆10 = 0.366 𝑇𝐷 = 0.366
𝑇𝐷 ∆𝑆10
 Flujo no lineal permanente hacia pozos en acuíferos confinados

Se puede demostrar que: Recordar zonas de flujo

𝑄 𝑟 > 𝑟𝐷 Régimen lineal o darciano
∆𝑆10 = 0.733 𝑟𝑡 < 𝑟 < 𝑟𝐷 Régimen no lineal
𝐾𝐷
𝑟 < 𝑟𝑡 Régimen turbulento puro
𝑄
𝑇𝐷 = 0.733
∆𝑆10

𝑈 = 𝑄/2𝜋𝑟𝑚

Para distancias 𝑟 muy próximas al pozo de bombeo 𝑟𝑜 − 𝑟 1

( )=
𝑟𝑟𝑜 𝑟
𝑄 𝑟𝑜 𝑄 2 1
𝑆 = 0.366 𝑙𝑜𝑔 ( ) + ( ) ( )
𝑇𝐷 𝑟 2𝜋𝑇𝑇 𝑟

𝐴
𝑄𝐴 𝑟2 𝑄𝐴 2 𝑟2 − 𝑟1
∆𝑆1−2 = 0.366 𝑙𝑜𝑔 ( ) + ( ) ( )
𝑇𝐷 𝑟1 2𝜋𝑇𝑇 𝑟2 𝑟1

𝐵
𝑄𝐵 𝑟2 𝑄𝐵 2 𝑟2 − 𝑟1
∆𝑆1−2 = 0.366 𝑙𝑜𝑔 ( ) + ( ) ( )
𝑇𝐷 𝑟1 2𝜋𝑇𝑇 𝑟2 𝑟1
Determinación de propiedades del acuífero con pruebas de pozo en RPNL

1/2
𝑟 − 𝑟1 1
(2 ) [ 𝑄 − 𝑄 ]
𝑟𝑟 (2𝜋 )2 𝐵 𝐴
𝑇𝑇 = [ 2 1 𝐵 𝐴 ]
∆𝑆1−2 ∆𝑆1−2
𝑆1 𝑇𝐷 ( )
(
0.366𝑄
) 𝑄𝐵 − 𝑄𝐴
𝑟𝑜 = 𝑟1 10
Flujo lineal impermanente hacia un pozo en acuífero confinado
𝑄 𝑄 𝑡2
∆𝑆𝑡10 = 0.183 ∆𝑆𝑡2−𝑡1 = 0.183 𝑙𝑜𝑔 ( )
𝑇𝐷 𝑇𝐷 𝑡1

𝑄 2.246𝑇𝐷 𝑡
𝑇𝐷 = 0.183 𝐾𝐷 = 𝑇𝐷 /𝑚 𝐸=
∆𝑆𝑡10 𝑟2
𝑄 2.246𝑇𝐷 𝑡
𝑆𝐷 = 0.183 𝑙𝑜𝑔 ( )
𝑇𝐷 𝑟2𝐸
 2.246𝑇𝐷 𝑡 1/2
𝑟𝑜 = ( )
𝐸
𝑄 2.246𝑇𝐷 𝑡
𝑆𝐷 = 0.183 𝑙𝑜𝑔 ( )
𝑇𝐷 𝑟 2𝐸
𝑄 𝑟𝑜
𝑆𝑐 = 0.366 𝑙𝑜𝑔 ( )
ℎ𝑜 𝐾𝐷 𝑟

Flujo no lineal impermanente hacia un pozo en acuífero confinado.

- El principal hallazgo es que para flujo impermanente 𝑆𝐷 es proporcional a 𝑄 y varía con

el tiempo, mientras que 𝑆𝑁𝐷 es proporcional a 𝑄2 y no depende del tiempo, luego:
𝑆𝑟 = 𝐵(𝑡 )𝑄 + 𝐷𝑄2
𝑄 2.246𝑇𝐷 𝑡 𝑄 2 𝑟𝑜 − 𝑟
𝑆𝑟 = 0.183 𝑙𝑜𝑔 ( )+( ) ( )
𝑇𝐷 𝑟2𝐸 2𝜋𝑇𝑇 𝑟𝑜 𝑟
𝑄 2.246𝑇𝐷 𝑡 𝑄 2 1
𝑆𝑟 = 0.183 𝑙𝑜𝑔 ( )+( ) ( )
𝑇𝐷 𝑟2𝐸 2𝜋𝑇𝑇 𝑟
Determinación de prop hidrogeológicas en régimen impermanente

2.246𝑇𝐷 𝑡 1/2
𝑟𝑜 = ( )
𝐸
𝑄 𝑄 𝑇𝐷 2
∆𝑆𝑡10 = 0.183
𝑇𝐷 𝑟𝐷 = 0.1𝜋𝑇𝐷 (𝑇𝑇 )

1/2
𝑄 2 𝑟2 − 𝑟1
( ) ( )
𝑇𝑇 =
2𝜋 𝑟1 𝑟2 2.246𝑇𝐷 𝑡 1/2
𝑄 𝑟 𝑟𝑜 = ( )
(𝑆𝑃𝑂1𝐴 − 𝑆𝑃𝑂2𝐴 ) − 0.366 𝑙𝑜𝑔 ( 2 ) 𝐸
[ 𝑇𝐷 𝑟1 ]

𝑄 2.246𝑇𝐷 𝑡𝐴 𝑄 2 𝑟𝑜 − 𝑟1
𝑆𝑃𝑂1𝐴 = 0.183 𝑙𝑜𝑔 ( )+( ) ( )
𝑇𝐷 𝑟12 𝐸 2𝜋𝑇𝑇 𝑟𝑜 𝑟1

𝑄 𝑇𝐷 2
𝑟𝐷 = ( ) 2.246𝑇𝐷 𝑡𝑜
0.1𝜋𝑇𝐷 𝑇𝑇 𝐸=
𝑟12
 2.246𝑇𝐷 𝑡𝐴
𝐸=
𝑄 2 1
𝑆𝑃𝑂1𝐴 −[(2𝜋𝑇 ) (𝑟 )]
𝑇 1
∆𝑆𝑡10 2.246𝑇𝐷 𝑡𝐴
𝐸=
𝑟12 ∙ 10( ) 𝑆𝑃𝑂 −𝑆𝑇
( ∆𝑆1 )
𝑟12 ∙ 10 𝑡10

Adaptación de la HP a las condiciones de campo

Pozos múltiples o interferencia de pozos.

𝐴
𝑄1 𝑟𝑜 𝑄2 𝑟𝑜 𝑄3 𝑟𝑜
𝑆𝑇𝑜𝑡𝑎𝑙 = 0.366 𝑙𝑜𝑔 ( 1 ) + 0.366 𝑙𝑜𝑔 ( 2 ) + 0.366 𝑙𝑜𝑔 ( 3 )
𝑇𝐷 𝑟1 𝑇𝐷 𝑟2 𝑇𝐷 𝑟3
-Si fuera libre el acuífero, estos abatimientos serían los corregidos y hay que aplicarles la linealidad.
-En los pozos de bombeo hay que adicionar la parte componente del propio pozo.

𝑃𝐵1
𝑆𝑇𝑜𝑡𝑎𝑙 = 𝑆𝐷𝑃𝐵1 + 𝑆𝑇𝑃𝐵1 + 𝑆𝐷𝑃𝐵3−1
Presencia de fronteras hidrogeológicas.

FLP hacia un pozo en acuífero con gradiente natural.