Universidad de Los Andes Facultad de Ingeniera Departamento de Vas Geotecnia para Hidrulica

Flujo hacia zanjas - pozos

Prof. Silvio Rojas Septiembre, 2009

A la memoria del Profesor Juan Francisco Lupini Bianchi

Prof. Silvio Rojas

Universidad de Los Andes Facultad de Ingeniera Departamento de Vas Geotecnia para hidrulica I.- FLUJO PROVENIENTE DE UNA LINEA FUENTE DE FILTRACIN I.1-FLUJO HACIA ZANJAS DE PENETRACIN TOTAL I.1.1- Acuifero confinado - Flujo artesiano

Lnea fuente

Prof. Silvio Rojas

Universidad de Los Andes Facultad de Ingeniera Departamento de Vas Geotecnia para hidrulica Para flujo estacionario (flujo establecido de bombeo), se plantea:

Lnea fuente

Donde: Q: Caudal de bombeo proveniente del estrato confinado k: Permeabilidad del estrato arenoso i:Gradiente de la superficie de abatimiento debido al bombeo (extraccin del caudal Q) dh: Variacin de la carga en un diferencial de longitud dy D: Espesor del estrato confinado x: Longitud de zanja considerada D.x: Area de filtracin considerada hacia la zanja, en el acuifero confinado.
Prof. Silvio Rojas

Universidad de Los Andes Facultad de Ingeniera Departamento de Vas Geotecnia para hidrulica

H:Carga total disponible en la lnea fuente, antes del abatimiento (carga correspondiente al acuifero confinado). L: Longitud de influencia en el abatimiento del nivel piezomtrico inicial, por el bombeo establecido en la zanja. h: Carga total en el estrato de arena, a una distancia "y" de la zanja y dentro de la zona de influencia del bombeo. hw: Altura de agua establecida dentro de la zanja para flujo estacionario.
Prof. Silvio Rojas

Integrando entre las condiciones de borde, establecidas en la figura 1, resulta:

Resolviendo la ec. 2:

Veamos la expresin de la carga a la distancia "y" a partir de la zanja:

Resolviendo la ec. 4:

Sustituyendo la ec. 3 en la ec. 5, se tiene: despejando la carga "h", se escribe: Esta ec. 6, define la geometra de la lnea piezomtrica abatida por el establecimiento del bombeo. Prof. Silvio Rojas

Universidad de Los Andes Facultad de Ingeniera Departamento de Vas Geotecnia para hidrulica Para obtener la expresin de "h", tambin se puede integrar entre:

Sustituyendo la ec. 3 en la ec. 8: Despejando "h", se obtiene: Esta ec. 9 es equivalente a la ec. 6, y tambin define la geometra de la lnea piezomtrica abatida por el establecimiento del bombeo.

Prof. Silvio Rojas

Universidad de Los Andes Facultad de Ingeniera Departamento de Vas Geotecnia para hidrulica Qu expresiones se obtienen, si se considera que el flujo es turbulento a travs del estrato de arena confinado? M. Anandakrishnan and Varadarajalu (1963) Journal of the Soil Mechanics and Foundations Division Proceedings of the American Society of Civil Engineers SM 5 "Laminar and Turbulent Flow of Water Through Sand"

donde: v: Velocidad del flujo a travs del suelo k: Coeficiente de flujo turbulento n: Exponente turbulento i: Gradiente hidrulico. De la ec. 10, se obtiene la velocidad "v":
Prof. Silvio Rojas

Universidad de Los Andes Facultad de Ingeniera Departamento de Vas Geotecnia para hidrulica Sabiendo que la expresin del caudal de infiltracin es el producto de la velocidad por el rea, se tiene entonces que para es te caso de flujo turbulento artesiano, la expresin del caudal ser ser:

Sustituyendo la expresin del gradiente, la ec. 12, queda:

Ordenando trminos e integrando, resulta:

Expresin del caudal para flujo turbulento

Universidad de Los Andes Facultad de Ingeniera Departamento de Vas Geotecnia para hidrulica Para la expresin de la lnea de abatimiento cuando se considera, flujo turbulento se plantea:

Igulando la ec. 16 y la ec. 18, se escribe:

La ec. 20 es la misma ec. 6 correspondiente a flujo laminar. Significa que la turbulencia no tiene influencia en el abatimiento de la lnea piezomtrica. Sin embargo la turbulencia debe ser determinante en la longitud "L" y la altura de agua dentro de la zanja.

Prof. Silvio Rojas

Universidad de Los Andes Facultad de Ingeniera Departamento de Vas Geotecnia para hidrulica I.1.2- Acuifero libre - Flujo gravitacional

Lnea fuente

Observaciones: 1.- Las equipotenciales no son verticales 2.- Las lneas de corriente no son horizontales 3.- El gradiente no es constante con la profundidad. Sin embargo se aplica que el gradiente se expresa como:
Prof. Silvio Rojas

4.- La altura de agua en la zanja (hw), no coincide con la altura de la superficie de abatimiento del N.F (superficie piezomtrica). Entonces se observa que existe una supercie libre de filtracin. Consideremos que "hw" coincide con la superficie de abatimiento una vez establecido un flujo estacionario por bombeo a travs de la zanja:

Donde: Q: Caudal de bombeo proveniente del acuifero libre k: Permeabilidad del estrato arenoso i:Gradiente de la superficie de abatimiento debido al bombeo (extraccin del cuadal Q) dh: Variacin de la carga en un diferencial de longitud dy h: Carga total en determinada seccin y que determina el rea de filtracin hacia la zanja x: Longitud de zanaja considerada h.x: Area de filtracin considerada hacia la zanja, en el acuifero libre. Prof. Silvio Rojas

H:Carga total disponible en la lnea fuente, antes del abatimiento (carga correspondiente al acuifero libre). L: Longitud de influencia en la rebaja del N.F, por el bombeo establecido en la zanja. h: Carga total en el estrato de arena a una distancia "y" de la zanja y dentro de la zona de influencia del bombeo. hw: Altura de agua establecida dentro de la zanja para flujo estacionario. Integrando la ec. 21, se tiene:

Carga a la distancia "y" a partir de la zanja: Resolviendo la ec. 4: Igualando la ec. 23 en la ec. 25:

Prof. Silvio Rojas

Universidad de Los Andes Facultad de Ingeniera Departamento de Vas Geotecnia para hidrulica Esta define la geometra de la lnea piezomtrica abatida por el establecimiento del bombeo. Sin embargo ella debe ser corregida por la altura (hs) de la superficie de infiltracin libre indica en la fig. 2. La expresin definitiva ser: Correccin de h por la altura (hs) de la superficie de infiltracin Para obtener la expresin de "h", tambin se puede integrar entre:

Sustituyendo la ec.23 en la ec. 28:

Prof. Silvio Rojas

Universidad de Los Andes Facultad de Ingeniera Departamento de Vas Geotecnia para hidrulica Despejando "h", se obtiene: Esta ecuacin, tambin debe ser corregida por la altura (hs) de la superficie de filtracin libre. Esta ec. 29 es equivalente a la ec. 26, para determinar la lnea de abatimiento del N.F en el caso de acuiferos libres. Qu expresiones se obtienen, si se considera que el flujo es turbulento a travs del estrato de arena confinado? Ecuacin que representa la turbulento en la masa de suelo. velocidad del flujo

Para el caso de flujo gravitacional, la ec. 13 se modifica encuanto al area de filtracin en el estrato, resultando: Ordenando trminos e integrando, resulta:
Prof. Silvio Rojas

Universidad de Los Andes Facultad de Ingeniera Departamento de Vas Geotecnia para hidrulica

La expresin del caudal ser: Caudal cuando se considera, flujo turbulento:

Igualando la ec. 33 y la ec. 35, se escribe:

Prof. Silvio Rojas

Lnea de abatimiento cuando se considera flujo turbulento Corriendo por la altura de filtracin en la zanja, se tiene:

La ec. 38 no es la misma ec. 29 o 26 correspondiente a flujo laminar. Significa que la turbulencia tiene influencia en el abatimiento de la lnea piezomtrica.

hs / H

hw / H

Grfica para flujo gravitacional en zanjas de penetracin total, para determinar la altura de la superficie de filtracin.

L/H

Prof. Silvio Rojas

Universidad de Los Andes Facultad de Ingeniera Departamento de Vas Geotecnia para hidrulica I.1.3- Flujo artesiano y Flujo gravitacional

Lnea fuente

Se estudia el caso de un acuifero confinado, pero donde la altura de agua (hw) en la zanja que se establece, debido al bombeo, alcanza una altura menor al espesor (D) del estrato confinado de arena.
Prof. Silvio Rojas

Lnea fuente

La red de flujo en la figura 3, indica que cerca de la zanja se presenta flujo gravitacional y a cierta distancia de la zanja ya el flujo se hace artesiano. Se considera que el punto donde la curva de abatimiento del nivel piezomtrico intersecta la lnea de estratificacin, que separa el estrato permeable del impermeable, es donde existe el cambio de artesiano a gravitacional. Este punto est ubicado a partir del borde de la zanaja con la distancia Li. Caso artesiano En este caso la carga vara entre h=D y h = H, entre respectivamente. las distancias y= Li, y=L,

Prof. Silvio Rojas

Universidad de Los Andes Facultad de Ingeniera Departamento de Vas Geotecnia para hidrulica Caso gravitacional En esta caso la carga vara entre h=hw y h=D, entre las distancias y=0 y y=Li, respectivamente.

Distancia Li De la ec. 40 y 42, se determina la longitud "Li":

Expresin del caudal La sustitucin de la ec. 44 en 42, resulta:

Prof. Silvio Rojas

Universidad de Los Andes Facultad de Ingeniera Departamento de Vas Geotecnia para hidrulica

Altura piezomtrica para el caso artesiano ( distancia "y" entre y=Li, y= L)

Igulando la ec. 40 y la ec. 46, resulta:

Prof. Silvio Rojas

Universidad de Los Andes Facultad de Ingeniera Departamento de Vas Geotecnia para hidrulica Altura piezomtrica para el caso gravitacional (distancia "y" entre y = 0, y = Li)

Igualando la ec. 42 y la ec. 49, resulta:

Corrigiendo por superficie libre de filtracin

Prof. Silvio Rojas

Universidad de Los Andes Facultad de Ingeniera Departamento de Vas Geotecnia para hidrulica I.2-FLUJO HACIA ZANJAS DE PENETRACIN PARCIAL I.2.1- Acuifero confinado - Flujo artesiano Generalmente este caso se presenta, debido a que el espesor del estrato permeable es considerable, y mecnicamente es costoso la penetracin completa del estrato.

Lnea fuente

Cuando la zanja es de penetracin parcial, la carga hidrulica residual "hd" aguas abajo de la zanja, a una distancia EA de la zanja, es mayor que la altura "hw" dentro de la zanja.
Prof. Silvio Rojas

Caudal hacia la zanja se determina a travs de:

Comentario (s.r) Es un planteamiento para llegar a la ec. del autor

La carga "hd" aguas abajo viene dada por se obtiene a travs de:

Comentario (sr)

Sustituyendo ec. 52 en 54:

Prof. Silvio Rojas

La distancia EA se estima de la figura:

Poca la diferencia

W: Penetracin del pozo en el acuifero confinado.

b: Ancho zanja.

Fig. Flujo artesiano para una zanja de penetracin parcial desde una lnea fuente.

Comentario (sr) El bombeo que extrae el caudal "Q a travs de la zanja, produce tanto la curva de abatimiento a la derecha de la zanja, como la curva de abatimiento a la izquierda. Por tanto se puede plantear:
Prof. Silvio Rojas

Es un planteamie nto de este desarrollo! Sustituyendo 52, 55 y 55, en 57, resulta:

Es un planteamiento de este desarrollo!

Prof. Silvio Rojas

I.2.2- Acuifero libre - Flujo gravitacional El flujo y descenso del nivel del agua, causado por el bombeo en una zanja de penetracin parcial en un estrato permeable, el cual es alimentado por una lnea fuente, es mostrado en la figura 5. Lo indicado son los resultados del modelo estudiado por Chapman ( ):

Lnea fuente

Lo indicado son los resultados del modelo estudiado por Chapman, vlido para:

Lo cual comprende el rango de la mayora de los problemas

El caudal hacia la zanaja puede ser determinado por:

La mxima carga aguas abajo de la zanja, puede ser obtenida por:

Comentario (sr) Estimando los flujos por la derecha y por la izquierda, se tiene: Es un planteamiento de este desarrollo!

Es un planteamient o de este desarrollo!

Prof. Silvio Rojas

II.- FLUJO PROVENIENTE DE DOS LINEAS FUENTES DE FILTRACIN II.1-FLUJO HACIA ZANJAS DE PENETRACIN TOTAL Generalmente el flujo hacia zanjas de longitud infinita se origina de dos lneas fuentes. Si se considera que la zanja est ubicada simetricamente y paralela respecto a las dos lneas fuentes, el flujo hacia la zanja se estima como dos veces el correspondiente a una sola lnea fuente. II.1.1- Acuifero confinado - Flujo artesiano Para una lnea expresiones fueron: fuente las

Para dos lneas fuentes, ser:

Prof. Silvio Rojas

II.1.2- Acuifero libre - Flujo gravitacional

Para una lnea fuente:

Para dos lneas fuentes:

Prof. Silvio Rojas

II.1.3- Flujo artesiano y Flujo gravitacional

Comportamiento del caso artesiano Para una lnea fuente: Comportamiento del Caso gravitacional Para una lnea fuente:
Prof. Silvio Rojas

Distancia Li De la ec. 40 y 42, se determina la longitud "Li":

Expresin del caudal La sustitucin de la ec. 44 en 42, resulta: Para una lnea fuente: Para dos lneas fuentes:

II.2-FLUJO HACIA ZANJAS DE PENETRACIN PARCIAL El bombeo en zanjas de penetracin parcial de gran longitud, en estratos permeables, ser simtrico con respecto a la zanja y puede ser considerado como el originado por dos lneas fuentes equidistantes y paralelas a la zanja.

Prof. Silvio Rojas

II.2.1- Acuifero confinado - Flujo artesiano

Caudal hacia la zanja se determina a travs de: :Factor el cual depende de la relacin entre W/D.

Prof. Silvio Rojas

Para zanjas de penetracin parcial, flujo artesiano.

Se considera que para distancias mayores de 1.3D, la carga "h" incrementa linealmente, y se calcula a partir de:

Para y>1.3D
Para y<1.3D, existe convergencia de flujo hacia la zanja, y el problema se resuelve a travs de la red de flujo. Aqu ya la carga "h" no vara linealmente con "y".

Prof. Silvio Rojas

Para y<1.3D

En este caso el caudal ser:

donde: h representa prdida de carga en la zona.

Prof. Silvio Rojas

II.2.2- Acuifero libre - Flujo gravitacional El flujo y descenso del nivel del agua, causado por el bombeo en una zanja de penetracin parcial en un estrato permeable, el cual es alimentado por dos lneas fuente, es mostrado en la figura 10. Lo indicado son los resultados del modelo estudiado por Chapman ( ):

Lo indicado son los resultados del modelo estudiado por Chapman, vlido para: Lo cual comprende el rango de la mayora de los problemas
Prof. Silvio Rojas

El caudal hacia la zanja proveniente de una lnea fuente paralela se determin a travs:

Para dos fuentes lineales:

La ec. (60) de Chapman, para flujo gravitacional proveniente de una lnea fuente en zanja de penetracin parcial, el flujo pas por debajo de la zanja y penetr aguas abajo de la zanja. Por tanto la ec. 70 para flujo simtrico respecto a la zanja no es exacto.

Prof. Silvio Rojas

II.3-FLUJO HACIA DOS ZANJAS DE PENETRACIN PARCIAL UBICADAS EQUIDISTANTES DEL CENTRO Y PARALELAS A DOS LINEAS FUENTES II.3.1- Acuifero confinado - Flujo artesiano Se rebaja el nivel piezomtrico del acuifero artesiano en una cierta longitud de excavacin (2l), lo cual se realiza con dos lneas de tubos filtrantes de penetracin parcial, espaciados cerradamente y simulados por dos zanjas como se muestra en la fig. 11. Se considera que el flujo es simtrico hacia las zanjas, y por tanto se estima que es originado por dos fuentes lineales ubicadas a la misma distancia de las zanjas.

Prof. Silvio Rojas

El caudal hacia las zanjas se estima por: Para una fuente lineal:

Para dos fuentes lineales: La carga "hd" en el centro de la excavacin, puede ser estima a travs de:

AL estimar la carga en este caso a travs de la ec. 53, se est considerando que las zanjas estn lo suficientemente separadas, de manera que no existe interferencia en la distribucin de las presiones de una zanja respecto a la otra. Si las zanjas estn prximas el estimado de "hd" es conservador.
Prof. Silvio Rojas

II.3.2- Acuifero libre - Flujo gravitacional Para hacer las estimaciones del caudal y carga, se aplicar el modelo de Chapman, vlido para: La rebaja del nivel fretico est hecho por las dos lneas de tubos filtrantes de penetracin parcial, cerradamente espaciados, equidistantes del centro de la excavacin, simuladas por dos zanjas y donde se considera adems que ambas lneas son alimentadas por dos fuentes lineales paralelas a las zanjas.

Prof. Silvio Rojas

La carga "hd" en el centro de la excavacin se estima a travs de:

donde: C1, C2: constantes obtenidas de las figuras indicadas abajo. El caudal hacia las zanjas se estima por: Para una lnea fuente. Para dos lneas fuentes

Constantes para corregir el flujo hacia dos zanjas de penetracin parcial en acuiferos gravitacional.

Prof. Silvio Rojas

IV.- CORRECCIN A UN SISTEMA DE POZOS REMPLAZANDO ZANJAS Engelud, demostr que a la reduccin de la carga en zanja, simulando un sistema de pozos, debera corregirse cuando se utiliza un grupo finito de pozos, remplazando la zanja. El procedimiento de Engelud, fue desarrollado para sistema de pozos de penetracin total; ellos pueden ser extendidos para aplicar a pozos de penetracin parcial. Fig. 13.- Lnea IV.1-FLUJO ARTESIANO de pozos de penetracin total representando a una zanja, en el caso de flujo artesiano. La fig 13a muestra la planta, fig. 13b se ilustra el perfil a lo largo de la hilera de pozos y la fig. 13c, se indica el corte desde la fuente al sistema de pozos. Prof. Silvio Rojas

Algunas observaciones: distancia "L" de la fuente a la hilera de pozos

1.- Considere una lnea infinita de pozos de penetracin total separados a una distancia "a"

la carga media "hm" el punto medio entre pozos carga "hd" aguas abajo de cada pozo. carga "hw" dentro de cada pozo

2.- La lnea de pozos estn alimentados por una lnea fuente de longitud infinita.

Prof. Silvio Rojas

Si los pozos son remplazados por una zanja de penetracin total de long. infinita, entonces la reduccin de la lnea piezomtrica del acuifero artesiano en una longitud de zanja "a", se estima a travs de la ec. 3, resultando:

donde: H: Carga total en la lnea fuente hw: Carga en el centro de la zanja. Qw: Caudal hacia la zanja de longitud "a". Todos los dems se conocen. Sin embargo existe una prdida de carga adicional en el caso de los pozos debido a la convergencia del flujo hacia los pozos. Esta prdida de carga adicional se estima, a travs de las siguientes expresiones: IV.1.1.- Correcciones para transformar la estimacin referidas a zanjas en estimaciones referidas a pozos: Para el centro del pozo de penetracin total:

Prof. Silvio Rojas

Para el centro del pozo de penetracin parcial:

donde: rw: Radio efectivo del pozo. a: Espaciamiento de los pozos. qa: Factor de correccin (fig. ) En el medio de entre pozos de penetracin total:

En el medio de entre pozos de penetracin parcial:

donde: qm: Factor de correccin (fig. ) Aguas abajo del pozo de penetracin total:

Prof. Silvio Rojas

Aguas abajo del pozo de penetracin parcial:

IV.1.2.- Abatimiento en pozos y cargas en pozos Ahora se aplicar las correcciones a los abatimientos correspondientes a las zanjas para estimar los abatimientos t cargas en pozos. Nota: A partir de este punto el caudal Qw, ya est referido al caudal del pozo de bombeo. Para el centro del pozo de penetracin total: Aplicando la ec. 75 a la ec. 74, se tiene:

Abatimiento correspondiente a un pozo de bombeo

Carga en el pozo de bombeo

Prof. Silvio Rojas

Para el centro del pozo de penetracin parcial: Aplicando la ec. 76 a la ec. 74, resulta: comentario(sr): comentario(sr): Se le debe sumar EA a la long. L? Para pozos de penetracin parcial en flujo artesiano, el caudal se expresaba, como:

Abatimiento correspondiente al centro del pozo de bombeo

Carga en el pozo de bombeo

En el medio de entre pozos de penetracin total: Carga "hm" La encontrar "hm", a la altura "hw" de la ec. 82, se le suma la correccin Dhm de la ec. 77

Prof. Silvio Rojas

Sustituyendo las expresiones mencionadas

Desarrollando:

hm

........... hm hw 85 ( )

En el medio de entre pozos de penetracin parcial: La encontrar "hm", a la altura "hw" de la ec. 84, se le suma la correccin Dhm de la ec. 78

Sustituyendo las ecuaciones mencionadas en 85, se tiene:

Prof. Silvio Rojas

Aguas abajo del pozo de penetracin total: Carga "hd" La encontrar "hd", a la altura "hw" de la ec. 82, se le suma la correccin Dhd de la ec. 80

Sustituyendo:

Se observa que la carga aguas abajo de la zanja es la misma hw de la ec 74, correspondiente a zanjas. Al comparar las cargas hw, hm y hd, determinadas con las ecuaciones 82, 86 y 90, para pozos de pemetracin total, se observa que: hw < hd < hm La ec. 88 y 89, son equivalentes, por tanto:

Prof. Silvio Rojas

De la ec. 90, se tiene:

Combinando la ec. 91 y 92:

Aguas abajo del pozo de penetracin parcial: Carga "hd" La encontrar "hd", a la altura "hw" de la ec. 84, se le suma la correccin de la ec. 80

Sustituyendo en la ec. 88, resulta:

La misma ec. 90, y valen los comentarios hechos anteriormente. De la ec. 94, se tiene:
Prof. Silvio Rojas

Relacionando 96 y 97, se tiene:

Fig. 14.- para obtener los factores de correccin