hidraulica de pozos - Docentes.unal.edu.co - Universidad Nacional ...

8CAPITULOHIDRAULICA DE POZOSLEONARDO DAVID DONADO GARZONTABLA DE CONTENIDOPág.1 INTRODUCCION 22 CONCEPTOS BASICOS 23 MOVIMIENTO NO PERMANENTE 33.1 POZOS DE PEQUEÑO DIÁMETRO 43.2 POZOS DE GRAN DIÁMETRO 204 MOVIMIENTO PERMANENTE 234.1 ACUÍFEROS CONFINADOS 234.2 ACUÍFEROS SEMICONFINADOS 264.3 ACUÍFEROS LIBRES 305 PRINCIPIO DE SUPERPOSICION 345.1 CASO DE DOS POZOS 345.2 MÉTODO DE LAS IMÁGENES 356 APLICACIONES 386.1 USO DE LA ECUACIÓN DE THEIS 386.2 USO DE LA ECUACIÓN DE JACOB 396.3 USO DE LA ECUACIÓN DE CHEN 396.4 USO DE LA ECUACIÓN DE PAPADOPULOS & COOPER 396.5 USO DE LA ECUACIÓN DE THIEM 406.6 USO DE LA ECUACIÓN DE DE GLEE - JACOB 406.7 USO DE LA ECUACIÓN DE DUPUIT - FORCHHEIMER 417 REFERENCIAS 41

2 CAPÍTULO 8 — HIDRÁULICA DE POZOS1 IINTRODUCCIIONUna vez determinadas las posibilidades de producción de agua subterránea en una determinada zona, elsiguiente proceso es determinar su adecuada explotación.Para una adecuada producción de los pozos de explotación de los acuíferos fuente, es necesario determinar eluso y así caracterizar de manera económica el beneficio de la explotación del recurso.A continuación, se presentan los diferentes métodos de análisis de pozos en los diferentes tipos de acuíferosexistentes. La intención es mostrar el desarrollo matemático de todas las ecuaciones que gobiernan elmovimiento del agua subterránea en explotación, ya sea bombeo o recarga de acuíferos.La principal aplicación planteada en este capítulo es la de determinar los radios de influencia de los pozos paraasí se necesita determinar que interferencia pueden tener entre ellos. Además con los conceptos explicados, setendrá la capacidad de determinar el abatimiento del nivel freático del acuífero en cualquier punto cuando se estaextrayendo agua.2 CONCEPTOS BASIICOSLa Figura 1 ilustra un pozo en una formación acuífera. En ella se detallan cada uno de los conceptos definidos acontinuación:• Nivel EstáticoEs el nivel de agua presente en la formación acuífera antes de comenzar el bombeo. Este nivel se ve afectadopor efectos meteorológicos (precipitación, infiltración) estacionales o por cargas adicionales (edificaciones), o porla descarga producida por pozos cercanos.• Nivel DinámicoTambién llamada nivel de bombeo, por que es producido cuando comienza la descarga de l acuífero por el pozo.Este nivel depende del caudal de bombeo, del tiempo de bombeo y de las características hidrogeológicas delacuífero. También se debe tener en cuenta la técnica desarrollada en el diseño de pozo.• AbatimientoBajo condiciones de extracción o inyección de un pozo, la carga hidráulica inicial en cualquier punto del acuíferocambia. En condiciones de extracción de un pozo, la distancia vertical entre la carga hidráulica inicial en unpunto en el acuífero y la posición baja de la carga hidráulica para el mismo punto es llamado abatimiento. Paraun acuífero libre el nivel del agua en el nivel freático está determinado por la distancia s(x,y,z,t), la cual es elabatimiento. Para el caso del acuífero confinado, el abatimiento es definido con respecto a la superficiepiezométrica. Este descenso de niveles, define la curva de abatimiento, por lo tanto es claro que el abatimientopresente su menor valor en lejanías del pozo y el mayor valor en el pozo. La dimensión del abatimiento es lalongitud [L]. El abatimiento es generalmente expresado en metros de agua

CAPÍTULO 8 — HIDRÁULICA DE POZOS 3Superficie del terrenobAbatimientoAcuífero libreCapa filtrante confinateAcuífero confinadoLecho impermeable2rwzQPozoSuperficie piezométricaantes del bombeoSuperficie piezométricaal tiempo tSuperficie piezométricaal tiempo t +∆tFigura 1 Esquema representativo del bombeo de un pozo.rQ(r)∆rQ(r+∆r)h0Datumh(r,t)• Cono de depresiónAl producirse el descenso del nivel estáticodel pozo, se establece un gradientehidráulico entre cualquier punto de laformación y el pozo, originándose unmovimiento radial desde todas lasdirecciones hacia el pozo en una formasimétrica y de tal manera que el caudal Qque se extrae del pozo es igual al caudalque pasa por cualquier sección del acuífero.A medida que la velocidad aumenta mayorserá el gradiente hidráulico ya que aumentala fricción existente entre el fluido y laspartículas sólidas en contacto; es por esoque lo que se forma alrededor del pozo se leconoce como cono de depresión que sobreun plano vertical presenta una curvaconocida con el nombre de curva deabatimiento. La forma, alcance yprofundidad de este cono de depresióndependerá de las condicioneshidrogeológicas (transmisividad y coeficiente de almacenamiento del acuífero), del caudal y el tiempo de bombeoo inyección. En el acuífero confinado el cono de depresión es la representación de la variación de los nivelespiezométricos en tanto que en el acuífero libre es además la forma real de la superficie piezométrica.• Capacidad EspecíficaEs la relación que existe entre el caudal que se obtiene de un pozo y el abatimiento producido y se expresa enunidades de caudal por longitud, [L 3 /T/L]. Este valor es contante para acuíferos confinados y variables para losacuíferos libres; es un término que representa el grado de eficiencia de un pozo ya que de dos pozos perforadosen una misma formación acuífera, el de menor capacidad específica tendrá menos eficiencia. El grado deeficiencia de un pozo lo determinaremos con base en la transmisividad y el coeficiente de almacenamiento de laformación acuífera, (con la cual podremos calcular un valor de la capacidad específica teórica) el valor de lacapacidad específica real medida en el pozo.3 MOVIIMIIENTO NO PERMANENTEEn 1935 Theis planteó el modelo matemático para describir el movimiento de agua subterránea en acuíferoshomogéneos e isotrópicos. Este modelo describe el flujo transiente en acuíferos bajo condiciones constantes deextracción de un pozo en acuíferos. A pesar de sus limitaciones tiene muchas aplicaciones en la hidráulica depozos. Trata el pozo como una línea origen y no toma en consideración el agua obtenida del almacenamientodentro del pozo. Papadopulos y Cooper generalizaron la ecuación de Theis considerando los efectos dealmacenamiento.

4 CAPÍTULO 8 — HIDRÁULICA DE POZOS3.1 POZOS DE PEQUEÑO DIÁMETRO3.1.1 Acuíferos confinados3.1.1.1 Consideraciones BásicasPara el cumplimiento del Modelo de Theis hay que tener en cuenta las siguientes consideracionesesquematizadas en la Figura 2.• Acuífero homogéneo e isotrópico• Acuífero horizontal y de espesor constante, b• Descarga contante, Q• No hay goteo• Acuífero de extensión infinita• El diámetro del pozo es infinitesimalmente pequeño, es decir que no existe almacenamiento en el pozo• El pozo penetra todo el acuífero• Antes del bombeo la carga piezométrica en el acuífero en la misma en cada punto del acuífero• La descarga del pozo es obtenida exclusivamente del almacenamiento del acuífero• El agua es inmediatamente liberada del almacenamiento del acuífero al declinar la carga hidráulica• El almacenamiento en el acuífero es proporcional a la carga hidráulicaSuperficie del terrenozQPozoSuperficie piezométricaantes del bombeoSuperficie piezométricaal tiempo tSuperficie piezométricaal tiempo t +∆tCapa confinateh0h(r,t)Acuífero confinador∆rb2rwQ(r)Q(r+∆r)Lecho impermeableDatumFigura 2. Flujo inestable en un pozo que penetra totalmente en un acuíferoconfinado. Sección transversal vertical.3.1.1.2 Ecuación de MovimientoUtilizando la Ecuación de Movimiento que gobierna en flujo en acuíferos isotrópicos:2⎛ ∂ h⎜2⎝ ∂x22∂ h ⎞ ∂ h+ + K2y⎟∂ ⎠ ∂zK 2=ST∂h∂t[3.1]Donde T es la transmisividad, S el coeficiente de almacenamiento y K es la conductividad hidráulica.

CAPÍTULO 8 — HIDRÁULICA DE POZOS 5Sabiendo queecuación:T = K b y S = S b , analizando el problema bidimensional, se obtiene la siguientes2 2∂ h ∂ h+ =2 2∂x∂yST∂h∂t[3.2]Utilizando coordenadas polares, dondecaudal es definida por v v ( r, z, t)radial.rr +2 2= x y y considerando la ley de Darcy, que en términos de( t)∂h r, z,=r= −K, donde K es la conductividad hidráulica en dirección∂rLa tasa total del flujo a una distancia r del pozo es:∂hQ() r = Arvr= −( 2 π r b)qr= 2 π r T[3.3]∂rLa carga piezométrica una distancia r es h(r,t); luego de un tiempo ∆t, la carga piezométrica es será h(r,t+∆t) y ladisminución de la carga piezométrica es:( t + ∆t) - h( r, t)Usando la figura 3.1, y aplicando la ecuación de continuidad:∆ h = h r,[3.4][ Q() r Q( r + ∆r)] ∆t= 2 r ∆r∆hS− π [3.5]y como ∆ r → 0 y ∆t→ 0∂Q∂r∂h= 2 π r S∂t[3.6]Reemplazando la ecuación 3.3 en la ecuación 3.6, se obtiene:Si el abatimiento está definido por:1 ∂ ⎛ ∂h⎞⎜r⎟ =r ∂r⎝ ∂r⎠2∂ h 1 ∂h+ =2∂rr ∂rSTSTs = h − 0h2∂ s 1 ∂sS ∂s+ =2∂rr ∂rT ∂t∂h∂t∂h∂t[3.7][3.8]Que es la ecuación de movimiento en flujo transitorio radial.3.1.1.3 Condiciones de FronteraSegún las suposiciones de Theis, las condiciones son las siguientes:

6 CAPÍTULO 8 — HIDRÁULICA DE POZOS• Para el AbatimientoCuando no se está extrayendo agua en cualquier punto del acuífero el abatimiento es nulo; es decir:∀ r, en t=0, s(r,t) = s(r,0) = 0En condiciones de extracción de agua, se supone que en la distancia más lejana del pozo, el abatimiento es nulo;es decir: ∀ t>0, en un radio r = ∞ , s(r,t) = s(r, ∞ ) = 0.• DescargaSi se tiene que en cuenta que sólo se produce abatimiento cuando se extrae agua, se concluye que:Cuando t < 0, Q = 0Cuando t ≥ 0, Q = constanteAhora, como la tasa de bombeo es constante en el pozo, de la ecuación 3.6, se tiene que para t ≥ 0 :⎛ ∂s⎞ Qlim rr 0⎜ ⎟ = −[3.9]→⎝ ∂r⎠ 2 π T3.1.1.4 Solución de la Ecuación de MovimientoPara encontrar la solución se aplica el método de separación de variables (Piskunov, 1977); es decir se busca lasolución particular de la ecuación 3.8 en forma de un producto de dos funciones:s( r, t) f( r) ⋅ g( t)Remplazando está función en la ecuación 3.8 se obtiene:= [3.10]1 Sf′g + f′g = fg′r Tf′′1 f′S g′+ =f r f T g[3.11]Al demostrar que son separables, estás funciones son iguales a una constante, que se llamará λ . Entoncesigualando λ al lado izquierdo de la ecuación 3.11:f′′1 f′+ = λf r f1f ′ + f′= fλr1f ′+ f′− fλ= 0rAl solucionar por operador cuadrático:1D 2 + D − λ = 0r

CAPÍTULO 8 — HIDRÁULICA DE POZOS 7D = −12r2[ 1 ± 1 + 4λr]Pero dependiendo del valor que tome el discriminante, se obtendrán las diferentes raíces, así:12+ , existen 2 raíces reales diferentes: D − [ 1 ± 1 + 4λr]i) Si 1 4λr2 > 0= .2r− 1 ⎡2 ⎤1+ + λ− ⎡21 1 4 r1−1+4λr ⎤1 2e2r ⎢⎣⎥⎦2r ⎢⎣⎥⎦Entonces la solución particular es: f () r = C e+ Cii) Sí 1 + 4λr2 = 0 , existen 2 raíces reales iguales: D = −−f122r2rEntonces la solución particular es: () r = C e + C r e .112+ , existen 2 raíces imaginarias diferentes: D = − 1 ± i − ( 1 + 4λr)−112r[ ]iii) Si 1 4λr2 < 02rEntonces la solución particular es:⎡ ⎛⎞ ⎛⎞⎤() = − 11212f r e ⎢ ⎜−− ( + λ ) ⎟ +2 ⎜−− ( + λ ) 2 rC1 cos 1 4 r C sen 1 4 r ⎟ ⎥⎦ .⎣ ⎝ 2r⎠ ⎝ 2r⎠Igualando ahora al lado izquierdo de la ecuación 3.11 a λ: λS g′= λT g ,S g′∫ = ∫ λIntegrando: T g , y luego despejando g(t) se llega a:Sln g = λt+TgλS() t P eTt( ) M= , donde P = eM = constantePor lo tanto, como de f se obtienen tres soluciones, la solución de la ecuación 3.8 puede tener tres formas:s r, t⎛= ⎜C⎝λ e1 − ⎡1+⎢⎣21+4λr⎤⎥⎦+ Cλ e1 − ⎡1−⎢⎣21+4λr⎤⎥⎦⎞⎟e⎠2r2rSi) ( ) ( ) ( )1−2rii) s( r, t) = C1( λ) e + C2( λ)λTt1( − ) ⎡ ⎛S 2riii) s( r, t)= e C λ cos11 λ−2r S( r e ) e⎢⎣2Tt12 ⎞ ⎛ 12 ⎞⎤( ) − − ( 1 + 4λr) + C ( λ) sen − − ( 1 + 4λ) ⎥⎦⎜⎝⎟⎠1 2r2rPara cada valor de λ , las constantes arbitrarias C 1 , C 2 y P tienen valores determinado; por eso C 1 y C 2 ; sonfunciones de λ y absorben el valor de P. También se aclara que la suma de las tres formas de solución sonsoluciones de la ecuación 3.8, debido a su linealidad..λTt⎜⎝2r⎟⎠

8 CAPÍTULO 8 — HIDRÁULICA DE POZOS3.1.1.5 Solución de TheisPara encontrar la solución y el valor de las constantes, Theis reemplazó las condiciones iniciales y de frontera enlas anteriores combinaciones, y así encontró la función de abatimiento por analogía de transferencia de calor ensólidos:A −us( r,t) = et[3.12]2r SDonde A es una constante y u = . Para t>0, el volumen total V, de agua tomado del acuífero es:4tTV∫ ∞= 2 π r s S dr0[3.13]Reemplazando 3.12 en 3.13:VV==∫∫∞0∞02 π r2 π rAtAtee−ur s−4tTS dr2[3.14]S drAl solucionar esta integral se tiene que:V= 2 πASt∫ ∞0e2− r s4tTr dr= 2 πASt⎡⎢−Tt e⎢ S⎣2 r2s−4tT⎤⎥⎥⎦r = ∞r = 0[3.15]De donde:VA==4 π T AV4 π T[3.16]Reemplazando 3.16 en 3.12 se tiene que:⎛ ⎞⎜r 2 S− ⎟⎝ 4Tt ⎠Vs ( r, t)= e[3.17]4 π T tEl Volumen de agua V, del acuífero es removido durante el período de tiempo dt. Así que V=Q t. y dV=Q dt, yentonces:

CAPÍTULO 8 — HIDRÁULICA DE POZOS 9⎛ ⎞⎜r 2 S− ⎟⎝ 4Tt ⎠dVds ( r, t)= e[3.18]4 π T tSi el agua es bombeada a una tasa de Q por unidad de tiempo de t=0 a t=t en el origen por integración seobtiene:dss( r, t)( r, t)⎛⎜r 2S ⎞−Q dt ⎜ 4Tt ⎟ ⎝ ⎠= e[3.19]4 π T tt⎛ 2r S ⎞⎜ ⎟−Q dt ⎜ 4Tt ⎟⎝ ⎠= e[3.20]4 π T∫t0Reemplazando:2r Su = , entonces:4Tt∞ -uQ es( r, t) = h0 - h( r, t) =du4 π T∫[3.21]uuDonde:∞ − u∫ueudu =−Ei( − u) = W( u)[3.22]La integral exponencial se conoce como la función de pozo de Theis, y su solución está dada por una serie depotencias:2 3 4u u uW() u = −0.5772− ln()u + u − + − + K2.2! 3.3! 4.4!nn uW() u = −0.5772− ln() u − ( − 1)n.n!∑ ∞n=1[3.23]Ahora se puede definir el abatimiento en términos de la curva de Theis:sQ= [3.24]4 π T( r, t) W()uLa Figura 3 muestra la curva típica de Theis, útil para determinar las parámetros hidrogeológicos de acuíferosconfinados usando datos de pruebas de bombeo. También se pueden trazar isolíneas de tiempo graficando elabatimiento en función del radio e isolíneas de radio, graficando el abatimiento en función del tiempo.

10 CAPÍTULO 8 — HIDRÁULICA DE POZOSFigura 3 Curva de Theis. (Batu, 1998)La ecuación es aplicable a acuíferos libres si el abatimiento es pequeño comparado con el espesor b de laformación. (Batu, 1998)2rcLa ecuación se cumple para la siguiente condición: t > 250 , donde r c es el radio del pozo, por no tener enTcuenta el almacenamiento en el pozo. Si los parámetros K, b, S y Q son conocidos, se puede determinar elabatimiento de la carga hidráulica en el acuífero confinado a cualquier distancia r del pozo, en cualquier tiempo.Lo único necesario es determinar el valor del parámetro u y así encontrar el valor de la función del pozo de Theis,W(u).Figura 4 Isolíneas de tiempo y de radio en función del abatimiento. (Batu, 1998)

CAPÍTULO 8 — HIDRÁULICA DE POZOS 113.1.1.6 Ecuación de JacobCooper & Jacob, 1946, tomaron en cuenta que cuando u,u < 0.01, la suma de los términos más allá de ln (u), enla ecuación 3.23, no es significativa. Los valores de u decrecen cuando el tiempo se incrementa y cuando ladistancia radial r decrece. Bajo esas condiciones:sQ≅ [3.25]4πT( r, t) [ − 0.5772 − ln()u ]( 0.5614)QQ ⎡ ⎤s ( r, t) ≅ [ ln( 0.5614) − ln( u)] =ππ ⎢ln4 T4 T⎥[3.26]⎣ u ⎦Reemplazando, u =2r S4Tt( r, t)⎡Q ⎢≅ ⎢ln4πT⎢⎢⎣⎤⎥ Q ⎡ 2.25 Tt ⎤⎥ =π ⎢ln⎥[3.27]2r S ⎥ 4 T ⎣ r S ⎦4Tt ⎥⎦( 0.5614)s2La ecuación 3.27 es conocida como la ecuación de Jacob. Que expresada en términos del logaritmo en base 10es igual a:2.302 Q ⎡ 2.25 Tt ⎤s( r, t) ≅π ⎢log24 T⎥ [3.28]⎣ r S ⎦Como primera aplicación de la ecuación de Jacob se puede usar para obtener el radio se influencia, cuando elabatimiento es nulo. Entonces despejando el Radio se obtieneQ ⎡ 2.25 Tt ⎤0 = ln24πT⎢R S⎥⎣ ⎦2.25 Tt0 = ln2R S12⎛ Tt ⎞R = 1.5⎜⎟ [3.29]⎝ S ⎠La ecuación de Jacob tiene la ventaja, respecto a la ecuación de Theis, de no requerir la consulta o tablas de lafunción de pozo de Theis.3.1.1.7 Capacidad Específica y Estimación de TransmisividadLa capacidad específica, CE de un pozo es definida como la relación de su descarga con su abatimiento total[CE=Q/s]; en otras palabras es el caudal por unidad de abatimiento. Se puede desarrollar una muy simpleecuación para estimar la transmisividad a partir de la capacidad específica, usando la ecuación de Jacob. Estaderivación está basada en un diámetro medio del pozo en un período promedio de bombeo, y valores típicos delcoeficiente de almacenamiento y producción específica.

12 CAPÍTULO 8 — HIDRÁULICA DE POZOSPara acuíferos confinados, Driscoll en 1986 (Batu, 1998) asumió los siguientes valores típicos:Tabla 1 Valores típicos para acuíferos confinados según Driscoll. (Batu, 1998)Parámetro Valor UnidadesTiempo, t 1 DíaRadio del pozo, rw 0.152 mProducción, S 0.001 AdimensionalTransmisividad, T 373 m 2 /díaSustituyendo estos valores en la ecuación de Jacob, se obtiene:Qs32[ m ][ ][ m2T ]s= CE m =día[ m] díaw1.38522T [ m ] = 1.385 CE [ m ]díadía[3.30][3.31]Para un acuífero libre, con producción específica S y =0.075, como valor típico, y el resto de valores mostrados enla tabla 1, se produce la siguiente relación:Qs32[ m ][ ][ m2T ]s= CE m =día[ m] díaw1.04222T [ m ] = 1.042 CE [ m ]díaSi se tienen múltiples pozos, la información obtenida de las anteriores ecuaciones puede usarse para estimar laconductividad hidráulica promedio (K med [m/d]) del acuífero, mediante la siguiente relación:∑∑día[3.32][3.31]KnLnKmed=[3.32]LDonde K es la conductividad de cada pozo, n es el número del pozo y L es la longitud del filtro.3.1.1.8 Ecuación de ChenEn 1984, Chen extendió la ecuación de Theis, para acuíferos de extensión lateral finita, como islas o meandros.Determinó que la distancia en la cual el abatimiento es nulo, en condiciones de bombeo, es conocida, y la llamaR. Es decir: s(R,t) = 0, donde R es la es la distancia radial donde la energía es cero. La solución encontrada seconoce como la Ecuación de Chen (Batú, 1998):nDonde:Qs ( r,t) = [ W( u) − W( U)+ 2I]4 π T[3.33]2R SU =4Tt

CAPÍTULO 8 — HIDRÁULICA DE POZOS 131⎡ 2⎛ u ⎞⎤J ⎢ ⎥0 ⎜ ⎟ χn∞ ⎢⎣⎝ U ⎠ ⎥ 1I =⎦∑ ( ) ∫ en=0 χnJ1χn0χn= R βnDonde:J 0 , J 1 : función de Bessel de orden cero y uno.β n : es la enésima raíz que satisface J 0 (R χ n ) = 0.( 1−x )⎡ 2U χn⎢−−⎢⎣x 4ULa Figura 5 muestra la gráfica u contra4 π T sQ , que es la usada para efectos prácticos. Se nota que cuando2R SU ≥ 4, la solución es igual a la de Theis. En otras palabras, sólo cuando t ≤ , se justifica usar este16Tmodelo.⎤⎥⎥⎦dxxSolución de Theis3.1.2 Acuíferos SemiconfinadosFigura 5 Curva de Chen. (Batu, 1998)Hantush y Jacob en 1955 (Batu, 1998), desarrollaron el modelo aplicable a acuíferos semiconfinados, isotrópicosy homogéneos, ilustrado en la Figura 6. Estos dos investigadores tuvieron en cuenta las siguientes suposiciones:

14 CAPÍTULO 8 — HIDRÁULICA DE POZOSQSuperficie piezométricaantes del bombeoSuperficie del terrenozSuperficie piezométricadurante el bombeo(cono de depresión)Nivel EstáticoswsrQvSuperficie piezométricaal tiempo t +∆tCapa confinantebAcuíferoconfinadohwH2rwLecho impermeableFigura 6. Acuífero semiconfinado• Acuífero homogéneo e isotrópico• Acuífero horizontal y de espesor constante, b, y su capa confinante posee un espesor constante b’ y unaconductividad hidráulica vertical K’.• Descarga contante, Q• Acuífero de extensión infinita• El diámetro del pozo es infinitesimalmente pequeño, es decir que no existe almacenamiento en el pozo• El pozo penetra todo el acuífero• La capa confinante no almacena agua• El flujo en el acuífero es horizontal y el goteo es vertical• Inicialmente, la tabla de agua posee la misma altura de la carga hidráulica del acuífero y es igual a h 0 .La ecuación diferencial parcial para flujo radial, fue obtenido por Jacob en 1946, es la ecuación que gobierna elmovimiento en este tipo de acuíferos. Aplicando el principio de continuidad, par el anillo dado, se tiene:[ Q() r Q( r + ∆r)] ∆t+ Qv ∆t= ( 2πr) ∆r∆hS− [3.34]( 2πr∆rb) v vQ v= A v =[3.35]h0− hUsando la Ley de Darcy: vv= K'[3.36]b'

CAPÍTULO 8 — HIDRÁULICA DE POZOS 15Y combinado las anteriores ecuaciones, se concluye que:h − h− [3.37]b'0[ Q() r Q( r + ∆r)] ∆t+ ( 2πr∆r) K' ∆t= ( 2πr) S∆r∆hComo ∆r y ∆s, tienden a cero y aplicando la definición de la derivada, se llaga a:∂∂rQ 0+h − hb'( 2πr) K' = ( 2πr)∂hS∂t[3.38]Y sabiendo que el abatimiento es el la diferencia de niveles, s=h 0 -h y que el caudal en el acuífero está dado por:∂hQ()r = 2π rT∂r[3.39]La ecuación es igual a:2∂ s2∂r+1 ∂sr ∂r−K'Tb's =ST∂s∂t[3.40]Y sí se reemplaza: B =Tb'K', la ecuación toma la forma:2∂ s 1 ∂ss+ − =22∂rr ∂rBST∂s∂t[3.41]Las condiciones iniciales y de frontera planteadas, para la cabeza y el abatimiento son:h r,0 = h , para todo r• ( )0• s ( r,0) = 0 , para todo rh ∞ , t = h , para todo t• ( )0• s ( ∞ , t) = 0, para todo tLas condiciones de descarga son:• Q = 0, cuando t=0• Q = constante, cuando t ≥ 0⎛ ∂h⎞ Qlim⎜r⎟ =r→0• ⎝ ∂r⎠ 2πT, para t ≥ 0⎛ ∂h⎞ Qlim⎜s⎟ = −r → 0• ⎝ ∂r⎠ 2πT, para t ≥ 0Al igual que Theis, Hantush y Jacob encontraron la solución a la ecuación de movimiento, la cual es:Q ⎛ r ⎞s ( r, t) = W⎜u,⎟ [3.42]4πT⎝ B ⎠

16 CAPÍTULO 8 — HIDRÁULICA DE POZOS⎛ r ⎞W⎜u,⎟Donde ⎝ B ⎠ es la función de pozo para acuíferos semiconfinados de Hantush y Jacob. Está funcióndescribe una serie, cuya expresión es:⎛⎜u,⎝rB⎞⎟ =⎠∞∫u1ue⎡2⎛ r ⎞ ⎤⎢ ⎜ ⎟ ⎥⎢ ⎝ B ⎠ ⎥⎢−u−u ⎥⎢ ⎥⎢⎣⎥⎦du[3.43]2r SAdemás, u = . La Figura 7 tabula los valores de la función de pozo, que también están en tablas en libros4Ttde matemáticas avanzadas e hidráulica de pozos.Figura 7. Curva de Hantush y Jacob. (Batu, 1998)3.1.3 Acuíferos LibresEn 1972, Neuman, aprovechando desarrollos realizados por Boulton (1954), (Batu, 1998) simplifico la ecuaciónde movimiento en acuíferos libre, ilustrados en la Figura 8. Las consideraciones que él tuvo en cuenta son:

CAPÍTULO 8 — HIDRÁULICA DE POZOS 17QSuperficie del terrenozNivel EstáticoSuperficie piezométricaantes del bombeoSuperficie piezométricadurante el bombeo(cono de depresión)swA 1rsSuperficie piezométricaal tiempo t +∆tKzFSKrbAcuíferolibreξH2rwA 2DatumLecho impermeableFigura 3.7 Flujo a un pozo en un acuífero libreinfinito• La tasa de bombeo es contante, Q• El diámetro del pozo es infinitamente pequeño• El pozo penetra completamente en el acuífero• En la zona saturada del acuífero , la ley de Darcy se cumple siempre• El acuífero tiene extensión lateral infinita• El material del acuífero es homogéneo pero anisotrópico, y su principal conductividad hidráulica estáorientada paralela a los ejes coordenados• El agua es bombeada por compactación del acuífero, expansión del aguay drena por gravedad de lasuperficie libre• El pozo puede ser tratados como una línea hundida• El abatimiento de la tabla de agua es pequeño comparado con el espesor de la zona saturada• Los efectos de capilaridad son despreciablesLa ecuación de movimiento fue plantada en el capítulo 7 y es:22∂ s Kr∂s∂ s S ∂sr+ + K = , 0 < z < ξ2z[3.44]∂rr ∂r∂zT ∂tK2La posición de la superficie libre de los acuíferos libres cambia en el espacio bajo condiciones de flujo transiente,por este motivo, la superficie libre es tratada como una frontera en movimiento. Bajo esta concepción, la fronterade la región de flujo, consiste de tres partes complementarias, mostradas en la Figura 8: La frontera de cargaprescrita, A 1 , la frontera de flujo prescrito, A 2 y frontera de la superficie libre, FS. Las otras fronteras tienden alinfinito. La pared del pozo se incluye en A 1 . Las condiciones iniciales para abatimiento, s(r,z,t) y espesorsaturado ξ(r,t), respectivamente son:

18 CAPÍTULO 8 — HIDRÁULICA DE POZOS• s(r,z,0) = 0• ξ(r,0) = bLa condición de frontera del abatimiento en el infinito es s ( ∞, z, t) = 0y en la frontera A 2 esLa condición de tasa de bombeo constante Q en el pozo está dada por la siguiente expresión:( r,0, t)∂s∂z= 0 .∞limr→0∫0∂sQr dz = −∂r2πKr[3.45]Neuman, simplificó la ecuación de movimiento, llegando a la siguiente expresión:2∂ s2∂r+1r∂s∂r+ KD2∂ s2∂z=1αs∂s,∂t0 < z < b[3.46]Donde:K K KK =zrzD= , αs= , αy[3.47]KrSsSy( r,b, t) 1 ∂s( r,b, t)∂ s∂z= −αy∂t[3.48]La solución encontrada por Neuman, para el abatimiento es:∞∞Q1 ⎡⎤2s( r, z, t) = 4xJ0[ x( KD) ]0( x) n( x) dx4 T∫⎢ω+ ∑ ω ⎥[3.49]π⎣n=⎦0 1Donde J 0 es la función de Bessel de primera clase de orden cero yωω002 2( ){ 1 − exp[ − tsKD( x − β0)]} cosh( β0zDbD)x =2⎪⎧2⎡2 222 b ⎤D ⎪⎫x − ( 1 + σ) β − ( x + β ) cosh( β b )⎨⎪⎩0⎢⎣0⎥⎬σ ⎦⎪⎭2 2( ){ 1 − exp[ − tsKD( x − βn)]} cosh( βnzDbD)x =2⎪⎧2⎡2 222 b ⎤D ⎪⎫x − ( 1 + σ) β − ( x + β ) cosh( β b )⎨n ⎢ n ⎥⎬⎪⎩⎣ σ ⎦⎪⎭Tt Tt b zts= , t2 y= , b2 D= , zD= , σ =Sr S r r byLas Figura 9 y 10 muestran la función de pozo de Neuman, en función del abatimiento relativo y el tiemporelativo.SSy0nDD

CAPÍTULO 8 — HIDRÁULICA DE POZOS 19Curva de TheisFigura 9 Función de Pozo de Neuman. (Batu, 1998)Curva deTheisFigura 10 Función de Pozo de Neuman. (Batu, 1998)

20 CAPÍTULO 8 — HIDRÁULICA DE POZOS3.2 POZOS DE GRAN DIÁMETROLos pozos de pequeño diámetro generalmente varían entre 0.05 m y 0.25 m. Como se mostró anteriormente,esos son representados por una línea en los modelos matemáticos. Esta aproximación es válida para los pozosen este rango de diámetros, pero inapropiada para pozos con un diámetro mayor. En particular, los radios debSuperficie del terrenoNivel EstáticoCapa confinanteAcuífero confinadoLecho impermeableFigura 11 Esquema representativa de un pozo de gran diámetro.mencionada del flujo en pozos de gran diámetro.2rwzswQrcSuperficie piezométricaantes del bombeoSuperficie piezométricadurante el bombeo(cono de depresión)Superficie piezométricaal tiempo t +∆tpozos excavados pueden ser de0.5 m a 2 m o más.La teoría de Theis asume que elpozo es una línea en el origen.Esta suposición no tiene encuenta los efectos significativosde almacenamiento. Los efectosde este almacenamiento en elpozo, llegan a ser importantescuando la transmisividad y elcoeficiente de almacenamientodel acuífero son pequeños ocuando diámetro del pozo debombeo es grande.Papadopulos y Cooper (1967)desarrollaron solucionesanalíticas en y alrededor depozos de gran diámetro enacuíferosconfinadoshomogéneos e isotrópicos,tomando en cuenta los efectosdel almacenamiento dentro delpozo. Después, Moensch (1985)presentó modelos matemáticosque combinaron los acuíferossemiconfinados de Hantush(1985) con la teoría antes3.2.1 Consideraciones BásicasLa Figura 11 muestra la sección transversal de un pozo de gran diámetro que penetra totalmente un acuíferoconfinado. Papadopulos y Cooper (1967) desarrollaron una solución analítica bajo condiciones de explotacióncon las siguientes suposiciones:• El acuífero es un homogéneo e isotrópico• El acuífero es horizontal y tiene un espesor constante (b)• La tasa de descarga (Q) del pozo es constante• El acuífero no tiene goteo y es horizontalmente infinito• El pozo penetra totalmente el acuífero• Las pérdidas en el pozo son despreciables• Antes del bombeo, la carga hidráulica en el acuífero es la misma en todos los puntos del acuífero• La descarga de los pozos es derivada exclusivamente del volumen almacenado en el acuífero• El agua es inmediatamente tomada en el bombeo, lo que hace decaer la carga hidráulica• El almacenamiento en el acuífero es proporcional a la carga hidráulicars

CAPÍTULO 8 — HIDRÁULICA DE POZOS 21La Ecuación de Movimiento es la misma ecuación 3.8, con la condición de que el radio,2∂ s 1 ∂sS ∂s+ =2∂rr ∂rT ∂tr ≥ r w.[3.34]Donde s es el abatimiento en el acuífero a una distancia r en un tiempo t; r es la distancia radial desde el centrodel pozo; S es el coeficiente de almacenamiento del acuífero; T es la transmisividad y r w es el radio efectivo dela pared del pozo.Las condiciones iniciales en el acuífero y el pozo, respectivamente, son:• r ≥ r , cuando s(r,0) = 0, s ww (0) = 0Las condiciones de frontera son:• s(r w ,t) = s w (t)• s( ∞ ,t)= 0• Almacenamiento dentro del pozo:( r ,t) ∂s() tsww2πr T∂ w−πr2c∂t∂t= −Qt ≥ 0Donde s w (t) es el abatimiento en el pozo a un tiempo t y r c el radio del pozo en el intervalo sobre el cual el nivel deagua decae. Las condiciones iniciales muestran que en un comienzo el abatimiento en el acuífero y en el pozoes cero. La primera condición de frontera indica que el abatimiento en el acuífero, en una cara del pozo es igualal abatimiento en el pozo. La segunda señala que el abatimiento en el acuífero en el infinito es cero. Finalmente,se expresa el efecto que tiene la tasa de descarga del pozo, que es iguala la suma de la tasa de flujo de aguadel pozo y la tasa de descenso en el volumen de agua dentro del pozo.3.2.2 Ecuación de Papadopulos & CooperEl problema planteado fue resuelto por Papadopulos & Cooper, mediante la transformada de Laplace (Batu,1998).Qs ( r, t) = F( u, α,ρ)[3.35]4πTDonde⎡F( u, α,ρ)8α=π⎤∞∫0DC()β() β β⎛ 2 ⎞⎜ 2 ρ−β ⎟C0⎜ ⎟() ⎢ ⎝ 4u ⎠β = 1 − e ⎥[ J ( βρ)() A β − Y ( βρ)()B β ]⎢⎣dβ2[3.36]0⎥[3.37]⎦A = βYβ − 2αβ[3.38]() β0() Y1()() β = βJ0 () β − 2αJ1() β22() β = [ A()β ] + [ B()β ] −B [3.39]D22r S rwS ru = , α = , ρ =2[3.41]4Tt r rcwJ 0 y Y 0 son las funciones de Bessel de orden cero y primera clase. Y 1 es la función de Bessel de primer orden yde segunda clase.[3.40]

22 CAPÍTULO 8 — HIDRÁULICA DE POZOS• Abatimiento dentro del PozoEl abatimiento dentro del pozo es obtenido cuando r = r w y puede ser expresado como:Donde:Los valores de ( u, α,ρ)Qs( r, t) = F( uw , α)[3.42]4πT( u , ) = F( u, ,1)F wα α[3.43]2r Su = ww [3.44]4TtF son tabulados por integración numérica de la ecuación 3.36. En la Figura 12, lossvalores son representados como una familia de cinco curvas de wcontra 1/u w ; una curva para cada unoQ4πTde los cinco valores del parámetro α. La curva de Theis, es también mostrada en la Figura 12, de la que seF u, α,ρ :obtienen importantes características de ( )El abatimiento predicho por la ecuación de Theis, se aproxima al abatimiento en el pozo de diámetro finito sólopara valores de tiempo relativamente grandes. Papadopulos (1967) comparó su aproximación con la Theis, así:3210 rcαρ 4F( u, α , ρ) ≈ W( u)para t > 2.5 , > 10[3.45]T uF( u , ) ≈ W( )α parawu w10 rTαu2ct > 2.5 , >w103[3.46]Las aproximaciones en las ecuaciones 3.38 y 3.39, son válidas para ambas condiciones: Para pozos que tienenun pequeño diámetro o acuíferos de transmisividad relativamente alta, el período definido en las anterioresecuaciones es muy pequeño. Así pues, para pozos de gran diámetro y acuíferos de baja transmisividad, esteperíodo es considerablemente largo.Sí 1/u w llega a ser suficientemente pequeño, las curvas se aproximan a líneas rectas que satisfacen la ecuación:Qt Volumen de agua descargada Q αsw= ==2[3.47]πr4 T ucÁrea del pozoπwoαF( uw, α ) =u[3.48]w

CAPÍTULO 8 — HIDRÁULICA DE POZOS 23Curva de TheisFigura 12. Curvas de Papadopulos y Cooper (Batu, 1998)En los primeros períodos, las líneas rectas representan las condiciones bajo la cual todo el agua bombeada esobtenida del almacenamiento dentro del pozo. Como resultado, los datos que están dentro del tramo de línearecta, de las curvas tipo, no dan información acerca de las características hidrogeológicas del acuífero.4 MOVIIMIIENTO PERMANENTEDespués de largos períodos de bombeo o recarga de un pozo, el flujo de aguas subterráneas alrededor de unpozo se aproxima al estado estable. Esto significa que la carga hidráulica del pozo en cualquier punto delacuífero no cambia con el tiempo. El período requerido para alcanzar el estado estable depende de lascaracterísticas hidráulicas del acuífero. Para los acuíferos menos permeables el período es más largo que paralos altamente permeables.Las soluciones de estado estable juegan un papel muy importante en el análisis de datos de abatimiento para ladeterminación de las características hidráulicas del acuífero y hacer el avalúo de la zona de influencia de un pozoo una batería de pozos.4.1 ACUÍFEROS CONFINADOS4.1.1 Consideraciones BásicasThiem (1906) fue el primero en derivar una solución para el flujo hacia un pozo en condiciones estables paraacuíferos confinados con base en las siguientes suposiciones:

24 CAPÍTULO 8 — HIDRÁULICA DE POZOS• Acuífero horizontal y con espesor constante• Acuífero homogéneo e isotrópico y de extensión lateral infinita• La carga hidráulica tiene una superficie horizontal antes del bombeo• La ley de Darcy es válida en el acuífero• El agua es instantáneamente removida del almacenamiento proporcionalmente con el decaimiento de lacarga hidráulica• La tasa del bombeo del pozo es contante• El flujo es simétrico con respecto al eje del pozo• La ecuación de movimiento (3.8), en flujo estable se reduce a:Superficie del terrenozQSuperficie piezométricaantes del bombeoSuperficie piezométricadurante el bombeo(cono de depresión)Nivel EstáticoswsrSuperficie piezométricaal tiempo t +∆tCapa confinantebAcuíferoconfinadohwKH2rwLecho impermeableFigura 13 Esquema representativa de un pozo en un acuífero confinado2∂ h 1 ∂h+ =2∂rr ∂rST∂h∂tT = K br[4.1]Para condiciones de flujo estable, el término de la derecha tiende a cero, entonces:2∂ h 1 ∂h+ = 02∂rr ∂r[4.2]

CAPÍTULO 8 — HIDRÁULICA DE POZOS 25Es necesario conocer las condiciones de frontera de Dirichlet (primer tipo), con referencia en la Figura 13:• h = h w Carga piezométrica conocida en la frontera del pozo• r = r w Radio del pozo• h = H Nivel de la carga piezométrica antes del bombeo• r = R Radio de influencia del pozo en el cual el abatimiento es cero4.1.2 Ecuación de ThiemUtilizando la ecuación de continuidad, a cualquier anillo concéntrico al pozo y teniendo en cuenta que se analizael proceso de bombeo, el caudal es negativo (si el pozo fuera de inyección el caudal sería positivo), se tiene que:− Q = AV = 2πrbv[4.3]Donde v r es la velocidad radial dada por la Ley de Darcy:Entonces:Resolviendo por variables separables:Para evaluar C, se aplican las condiciones de frontera:vr( )r( z, t)∂h r,= vr( r, z, t)= -K[4.4]∂r∂hQ = 2πrbK[4.5]∂rQK∂ h = ∂r[4.6]2πrb()Q ln rh () r = + C[4.7]2 π TSi h = h w , entonces r = r w ; por lo tanto:( )Q ln rwhw = + C2 π TQ ln( rw)C = hw−2 π TReemplazando la ecuación 4.5, se obtiene la Ecuación de ThiemQ ln r Q ln rwh()r = + hw−2 π T 2 π TQh r − h = ln r − ln2 π TQ ⎡ ⎛ r ⎞⎤h() r − hw= ⎢ln⎜⎟⎥ 2 π T ⎣ ⎝ rw⎠ ⎦() ( )() [ () ( )]wr w[4.8]Analizando la Ecuación de Thiem se puede concluir que:

26 CAPÍTULO 8 — HIDRÁULICA DE POZOS• La carga piezométrica h se incrementa asintóticamente con el incremento de la distancia radial r• La superficie piezométrica no puede ascender sobre h(r).• Es válida sólo en la proximidad de un pozo donde el flujo estable ha sido definido.Con la Ecuación de Thiem, se puede predecir el Radio de Influencia de un pozo, en términos del abatimiento enel mismo cuando h = H, r = R,Q ⎡ ⎛ R ⎞⎤Q ⎡ ⎛ R ⎞⎤h( R) − hw= ⎢ln⎜⎟⎥ ∴ H − hw= ⎢ln⎥2 π T ⎣ ⎝ r⎜⎟w ⎠ ⎦ 2 π T ⎣ ⎝ rw⎠ ⎦s = h RQ ⎡ ⎛ R ⎞⎤( ) − h( r ) = ln⎜⎟⎥ ⎦w⎢2 π T ⎣Esta forma de la ecuación de Thiem, posee las siguientes características:• La distancia R, para la cual el abatimiento es cero, es el radio de influencia del pozo.• El parámetro R tiene que ser estimado antes de la predicción de los abatimientos.4.2 ACUÍFEROS SEMICONFINADOSLa Figura 14 muestra un pozo que penetra totalmente un acuífero semiconfinado, a través del cual la filtraciónproviene de un acuitardo superior. La solución propuesta independientemente por De Glee & Jacob, se basa enlas siguientes suposiciones:• El acuífero es limitado abajo por un lecho impermeable, y arriba por una capa semiconfinate.• Sobre la capa semiconfinante, existe un acuífero libre que tiene una tabla de aguas horizontal, cuya cargahidráulica es constante (h 0 ). El suministro de agua al acuífero libre es suficiente para mantener h 0 constante.• El flujo en la capa semiconfinante es vertical• Las mismas suposiciones del acuífero confinadoAplicando la ecuación de continuidad a cualquier anillo de radio r, mostrado en la Figura 14 se tiene que:⎜⎝ rw⎟⎠( r ∆ r) − Q( r) + ( 2 π r ∆ r) v 0[4.9]Q + =[4.10]vDonde v v es la velocidad de goteo desde la capa semiconfinate. Si se divide por ∆r y como ∆r tiende a cero, sellega a:⎡Q rlim∆r→0⎢⎣( + ∆r) − Q( r)∆r+( 2πr) v = 0v⎤⎥⎦∂ Q∂r+ 2 π rvv =0[4.11]La Ley de Darcy por el acuífero semiconfinado, conduce a:

CAPÍTULO 8 — HIDRÁULICA DE POZOS 27Superficie del terrenozQSuperficie piezométricaantes del bombeoSuperficie piezométricadurante el bombeo(cono de depresión)Nivel EstáticoswsrK’ vvSuperficie piezométricaal tiempo t +∆tb’CapasemiconfinantebAcuíferosemiconfinadohwKH02rwDatumLecho impermeableFigura 14. Esquema representativa de un pozo en un acuífero semiconfinado.( + ∆r) − Q( r)Q r( 2 r) vv0∆r+ π =∂hQ() r = ( 2πrb)K∂r∂hQ()r = 2π rT∂r[4.12]La Ley de Darcy también controla la velocidad de goteo:vvh0− h= K'[4.13]b'

28 CAPÍTULO 8 — HIDRÁULICA DE POZOSDonde K’ y b’ son la conductividad hidráulica y el espesor de la capa confinante (acuitardo).Reemplazando 4.12 y 4.13 en 4.10, se obtiene:∂∂⎛ ∂h⎞∂⎜2πrT⎟⎝ ∂r⎠ ⎛ h0 − h ⎞+ 2πr⎜K'⎟ = 0∂r⎝ b' ⎠ ∴⎛ ∂h⎞ ⎛ h0− h ⎞ 1 ∂ ⎛ ∂h⎞ h0− h⎜r⎟ + r⎜⎟ = 0 ∴ ⎜r⎟ + = 02⎝ ∂r⎠ ⎝ B ⎠ r ∂r⎝ ∂r⎠ Br2[4.14]b ⋅b'⋅KDonde B 2 = , es llamado factor de filtración. En la misma ecuación b'/K' es conocida como laK'resistencia hidráulica.La ecuación puede ser escrita como una Ecuación de Movimiento ordinaria porque h sólo depende del radio r.Reemplazando s = h - h 0, , en la ecuación:1rddr⎛ dh ⎞ h − h⎜r⎟ +2⎝ dr ⎠ B0=Si r/B=x, entonces r = B x, y dr = B dx,1 d ⎛ ds ⎞ s0 ⎜r⎟ 02∴ r dr ⎝ dr ⎠+ B=21 ds d s s+ − = 02 2∴ r dr dr B ∴2d s ds 2 s+ r − r = 02[4.15]dr dr B2r222d s ds s2 d s ds 2Bx + ∴ x + x − x s = 02 222[4.16]B dx Bdx B dx dx22 2( ) Bx − B x = 04.2.1.1 Ecuación de De Glee - JacobLa ecuación 4.16 es un a ecuación diferencial no lineal no homogénea que se puede solucionar por el método deCauchy – Euler. Este método consiste en aplicar la regla de la cadena luego de hacer el siguiente reemplazo:u = ln( x).Entonces:e u 2u 2= x ∴ e = xAhora se encontrarán las derivadas:12d s=2dxddxdsdx⎛ 1⎜⎝ x2d s 1= −2 2dx xReemplazando en la ecuación 4.16, se llega a:dsdudu =dx xds du 1 ds= =du dx x dudsdu2⎞ 1 ds 1 d s⎟ = − +22⎠ x du x du221 d s 1 ⎛ d s+ =⎜ −2 2 2 2x du x ⎝ dududxdsdu⎞⎟⎠

CAPÍTULO 8 — HIDRÁULICA DE POZOS 292d s2duds dse2 u2− + = s ∴ s′ − x s = 0[4.17]du duAhora se soluciona esta ecuación por medio del operador cuadrático, explicado en el numeral 3.1.1.4,obteniendo las raíces D = ± x .Entonces, la solución es igual a:s = Cx −x1e+ C2err−Bs C1e+ C2erBQue al reemplazar el valor de x = , se obtiene: = . Para encontrar el valor de lasBconstantes, se sabe que cuando el radio tiende al infinito, el abatimiento es nulo, y que cuando el radio es igual al0 = 1, pero esto esindeterminado, lo cual no permite concluir el valor de C 1 . Es decir que la solución planteada no es compatiblecon las condiciones de frontera dadas. Sólo se sabe que las soluciones son linealmente independientes, por loque usando las función de Bessel se pueda encontrar la solución compatible.radio del pozo, el caudal es constante. Entonces: si s(r) = s( ∞ ) = 0, entonces C ( ∞)⎛ r ⎞ ⎛ r ⎞= C1 I0⎜ ⎟ + C2I⎜ − ⎟[4.18]⎝ B ⎠ ⎝ B ⎠s0Donde:⎛ r ⎞I 0 ⎜ ⎟ : es el modificador de la función de Bessel de orden 0, de primera clase.⎝ B ⎠⎛ r ⎞ ⎛ r ⎞I0 ⎜ − ⎟ = K0 ⎜ ⎟ : es el modificador de la función de Bessel de orden 0, de segunda clase, así que la⎝ B ⎠ ⎝ B ⎠solución queda definida como:⎛ r ⎞ ⎛ r ⎞s = C1 I0⎜ ⎟ + C2K0 ⎜ ⎟[4.19]⎝ B ⎠ ⎝ B ⎠Los valores de los modificadores de Bessel están ya tabulados, y se pueden encontrar en libros de CálculoAvanzado.Con las condiciones de frontera antes mencionadas, se concluye que: ( ∞ ) = ∞,K ( ∞) 0igual a cero.I0 0=, así que C 1 esAhora cuando r = r w:∂sQ = −2πrwbKdr⎛ r ⎞ ∂s∂ ⎛ ⎛ r ⎞⎞1 ⎛ r ⎞ ⎛ r ⎞Si s() r = C2 K0 ⎜ ⎟ entonces, = C2 ⎜K0 ⎜ ⎟⎟= C2K1⎜⎟ , donde K 1 ⎜ ⎟ , es el operador de⎝ B ⎠ ∂r∂r⎝ ⎝ B ⎠⎠B ⎝ B ⎠ ⎝ B ⎠la función de Bessel de primer orden de segunda clase.

30 CAPÍTULO 8 — HIDRÁULICA DE POZOSQ = 2πrw⎛ 1T⎜ C⎝ B2⎛ rw⎞⎞K1⎜⎟⎟ ∴ C⎝ B ⎠⎠Y reemplazando el valor de la constante, se tiene que:2=2πrwQ⎛ 1T⎜ K⎝ B1⎛ rw⎞⎞⎜ ⎟⎟⎝ B ⎠⎠s=Q2πT⎛ r ⎞K0 ⎜ ⎟⎝ B ⎠rw⎛ rwK1⎜B ⎝ B⎞⎟⎠[4.20]La ecuación 4.20 representa la Ecuación de DeGlee – Jacob, para acuíferos semiconfinados. Está ecuaciónpuede simplificarse para usos prácticos; si r w 0. 01B< , se puede aproximar el factor: rw⎛ rw⎞KB B11⎜⎟ ≅ , y⎝ ⎠entonces la ecuación 4.18 se puede escribir como:Q ⎛ r ⎞s = K0 ⎜ ⎟[4.21]2πT⎝ B ⎠rSegún Hantush (Batu, 1998), si ≤ 0.05 la ecuación de De Glee - Jacob se puede escribir como:B2.303 Q ⎛1.12 B ⎞s () r ≅ log⎜⎟ [4.22]2 π T ⎝ r ⎠4.3 ACUÍFEROS LIBRESDupuit y Forchheimer derivaron la expresión sin reconocer quien la hizo primero, por esta razón lleva ambosnombres. Esta ecuación es la simple ecuación para acuíferos libres.4.3.1 Consideraciones BásicasLa Figura 15 muestra un pozo que penetra completamente el acuífero libre. Dupuit y Forchheimer encontraronindependientemente la solución para la carga piezométrica con base en las siguientes suposiciones:• El acuífero es homogéneo e isotrópico y de extensión infinita• La tabla de aguas es horizontal antes del bombeo• La ley de Darcy es valida para el flujo en el acuífero• El agua es instantáneamente removida del almacenamiento, como la carga piezométrica decae.• La tasa de bombeo del pozo es constante• Las condiciones de Dupuit son validas.• El flujo es simétrico, respecto al eje del pozo. La filtración de las paredes del pozo es despreciable y elacuífero recibe una tasa constante de recarga.• Se desprecian las pérdidas en el pozo, H 0 =H w

CAPÍTULO 8 — HIDRÁULICA DE POZOS 31Dupuit en 1863 (Batu, 1998) indicó que la pendiente de la tabla de aguas, de un acuífero libre bajo condicionesde no extracción a lo largo de una sección transversal vertical es muy pequeña. El rango de valores típicos va de1/1000 a 1/10000. Alrededor de un pozo de extracción en un acuífero libre la pendiente es muy alta, con eldescenso de la distancia radial del pozo dependiendo de las conductividades hidráulicas verticales y horizontalesdel acuífero. La condición de una pendiente geométrica pequeña significa que el flujo es esencialmentehorizontal y la carga hidráulica (h) es igual a la elevación de la tabla de aguas.Superficie del terrenoQzNivel EstáticoswsSuperficie piezométricaantes del bombeoSuperficie piezométricadurante el bombeo(cono de depresión)rSuperficie piezométricaal tiempo t +∆tAcuíferolibreHoHwH2rwDatumLecho impermeableFigura 15 Flujo a un pozo en un acuífero libre infinito con filtraciónRealizando el análisis de continuidad, en un anillo de radio r, se tiene que:.Q r + ∆r − Q r + 2πr∆rI =[4.23]( ) ( ) ( ) 0Donde I [L/T] representa el volumen de agua entrando en una unidad de área horizontal del acuífero por unidadde tiempo, debido a la recarga por infiltración. Los valores positivos y negativos de I, representan la recarga y laevaporación respectivamente. Al dividir por ∆r, y haciendo tender este a cero:∂Q+ 2π rI = 0∂r[4.24]

32 CAPÍTULO 8 — HIDRÁULICA DE POZOSLa velocidad radial de Darcy está dada por:Entonces:Qv r( r) = −( 2 rh) q = −( 2πrh)r∂h= −K∂r21 ∂π ∴ () ( )( h )Q r 2 r K∂hK∂r= π [4.25]2 ∂rY reemplazando 4.25 en 4.24:1r2( h )∂ ⎡ ∂r∂r⎢⎣ ∂r⎤⎥ +⎦2I= 0K[4.26]Las condiciones de primer tipo o de Dirichlet son:• Cuando h = H w , entonces r = r w . Carga piezométrica en la cara del pozo• Cuando h = H, entonces r = R. Carga piezométrica del acuífero antes del bombeo.• R es el radio de influencia en el cual el abatimiento es cero.2 222 ∂ ( h ) ∂( h ) 2I 2La ecuación 4.26 puede ser escrita de la forma: r + r = − r , y ser solucionada por el2∂r∂rKmétodo de Cauchy – Euler, ya que es una ecuación no homogénea. Realizando la siguiente cambio de variable,se tiene que:Y las derivadas son:u = lneeu2u= x( x)= x22 2 22 2 2( h ) 1 ∂( h ) ∂ ( h ) 1 ⎛ ∂ ( h ) ∂( h ⎞=y =⎜ −)⎟⎟2 2 2∂∂rr∂u∂rremplazando estas en la ecuación 4.26, se llega a que:2r⎝∂u∂u⎠Integrando dos veces:Reemplazando u por ln (r):2∂∂r2( h )2( h ) 2I 2u2= −K∂ I 2u= − e + C1∂rKIh +2K22u() r = − e + C1u C2e2 I 2h () r = − ( r ) + C1ln() r + C2[4.27]2KReemplazando las condiciones de frontera se obtienen C 1 y C 2 .

CAPÍTULO 8 — HIDRÁULICA DE POZOS 332 2 I 2 2 ⎡ 2 2 I 2 2() = H + ( R − r ) − H -H + ( R − r )h r2K⎢⎣w2Kw⎛ r ⎞ln⎜⎟⎤ ⎝ R ⎠⎥⎦ ⎛ rwln⎜⎝ R⎞⎟⎠[4.28]La descarga Q del pozo puede ser determinada por continuidad y la ley de Darcy:2∂hd( h )Q = -2π rwhqr= 2πrwhK= πKrw, r = rw[4.29]∂rdrDerivando 4.28 y reemplazando en 4.29:Q = −I2⎡2I2 2( R − r )⎤πK2π rw−⎢H - Hw+w⎣⎥[4.30]2K ⎦ ⎛ rw⎞ln⎜⎟⎝ R ⎠Teniendo en cuenta que:2− Iπr w: es la recarga en el pozo mismo y es despreciablemente pequeño comparado con los demás términos:⎡Q = −⎢H⎣I+2K2 2( R − r )⎤ K⎥⎦ ⎛ rwln⎜⎝ R2 2π- Hww⎞⎟⎠[4.31]Reemplazando 4.31 en 4.29:h2= H2+I2K2 2 Q ⎛ r ⎞( R − r ) + ln ⎟ ⎠πK⎜⎝ R[4.32]En el caso particular, en el que I = 0:h2= H2Q+πK⎛ln⎜⎝rR⎞⎟⎠[4.33]Las ecuaciones 4.32 y 4.33 representan la distribución de la carga piezométrica con recarga y sin recargarespectivamente. Para la condición descrita en las condiciones de frontera, la tasa de descarga Q, puede serrepresentada como:2 2πK( H − Hw)Q =⎛ R ⎞ln⎜⎟[4.34]⎝ rw⎠La ecuación 4.34 es la llamada ecuación de descarga de Dupuit - Forchheimer. Esta ecuación es obtenida conbase en las condiciones de Dupuit. Estas suposiciones no toman en cuenta la forma curvilínea del flujo en unplano radial. Los componentes del flujo vertical son despreciados. La ecuación da un resultado con razonableaproximación, si la distancia radial r es suficientemente grande y los efectos curvilíneos son despreciables.Luego, la aplicación de métodos numéricos (Boulton, 1951 (Batu, 1998)) e investigaciones experimentales((Babbit y Cantwell, 1948) (Peterson et al, 1952) Batu, 1998) muestran que la ecuación representa la superficielibre para valores de r ≥ 1.5H, siempre y cuando el nivel de agua del pozo (H o ) sea cero en la Figura 4.3.

34 CAPÍTULO 8 — HIDRÁULICA DE POZOSAdicionalmente ocurre que para valores pequeños de r, H 0 se incrementa. Estas investigaciones tambiénmuestran que la superficie libre cercana y a la misma distancia alrededor del pozo no es correctamente modeladapor la ecuación 4.34 y que la superficie libre cruza la pared del pozo a alguna distancia sobre el nivel del agua enel pozo. Hantush (1962) (Batu, 1998) analizó la validez de la ecuación 4.34, tomando en cuenta el nivel de aguaen el pozo (H o ) y el nivel del agua en la pared del pozo (H w ) y obtuvo la misma ecuación 4.32 con la excepción deque H w = H 0 . Esto significa que al tomar en cuenta la naturaleza curvilínea del flujo, a lo largo de la pendiente defiltración, virtualmente se obtiene la misma ecuación.5 PRIINCIIPIIO DE SUPERPOSIICIIONEste principio (Quintero, 1994) se encarga de analizar la interferencia entre una batería de pozos en unaformación acuífera, y el efecto que presenta este en la producción de los mismos. Como en la realidad, seencuentran los acuíferos con limitaciones hidrogeológicas definidas, que restringen la aplicabilidad de losmétodos analíticos, que suponen la extensión infinita de los acuíferos, como lo muestra las Figuras 16, 17 y 18.El método de las imágenes se utiliza par resolver teóricamente estos casos, aproximando una extensión finita delos acuíferos, con un pozo real y otro imagen. Basado en la linealidad de la Ecuación de Laplace (Para acuíferoslibres, se mantiene si sí la variable de estado es h 2 y no h), suponiendo el trabajo de cada pozo y luegosuperponerlos, para así obtener la resultante de todos los pozos trabajando en conjunto.5.1 CASO DE DOS POZOSSuponiendo que en un acuífero confinado se tienen dos pozos, separados a una distancia 2 a, como lo muestrala Figura 19. Los pozos están diseñados en igual forma, y están localizados en forma tal que a una distanciaradial el potencial permanece constante. El caudal que se extrae de ambos, es el mismo, Q.De acuerdo al principio de superposición el abatimiento total producido en un punto P(x,y) será la suma de losabatimientos que produce cada pozo en su operación individual, por lo tanto:Q ⎛ ⎞0R Qs = s1+ s2= 0.366 log⎜⎟ + 0.366T ⎝ r1⎠ T2Q ⎛ ⎞0Rs = 0.366 log⎜⎟T ⎝ r1r2 ⎠0⎛ R ⎞log⎜⎟⎝ r2⎠[5.1]Por lo tanto el abatimiento total en cada pozo será:Q ⎛ ⎞0 ⎜Rs =⎟p0.366 log[5.2]T⎝ rp2a ⎠El caudal que produce cada uno seráQ0=T s⎛ ⎞⎜R0.366log ⎟⎝ rp2a ⎠p[5.3]Como lo muestra la ecuación 5.3, el caudal disminuye a medida que disminuye la distancia 2 a, entre pozos.

CAPÍTULO 8 — HIDRÁULICA DE POZOS 355.2 MÉTODO DE LAS IMÁGENES5.2.1 Pozo cerca de una zona de recargaEste es el comportamiento típico de un pozo situado en cercanías de un río y perforando un acuífero que está encontacto directo con el río el cual se extiende linealmente en una gran distancia. La Figura 19, representa la zonade recarga como una línea que se extiende a lo largo del eje Y y a una distancia a se encuentra un pozo del cualse bombea un caudal determinado, Q.La zona de recarga se puede simular con dos pozos separados a una distancia 2 a, y en forma tal que uno deellos, el pozo imagen, es un pozo de recarga. Estos dos pozos producen a lo largo del eje y, la condición s=0.La solución está dada por la ecuación:Q ⎛ ⎞0r2s = ln⎜⎟[5.4]2πT⎝ r1⎠Si r 2 = 2a y r 1 = r p , se tiene que el caudal Q, es igual a:encuentra el abatimiento.2πTspQ =. Aplicando el teorema del coseno, se⎛ ⎞⎜2aln ⎟⎝ rp⎠22Por lo tanto, reemplazando en 5.4, se obtiene:rr2= r=2121r+ 4a − 4ar1+ 4a − 4arcosβ1cos β⎛ 2 2Q⎞⎜ r1+ 4a − 4ar1cos β0s = ln⎟[5.5]2πT⎜ r ⎟1⎝⎠Para los puntos paralelos a la línea de recarga, es decir cuando β=90º, el cos (90) = 0, y la ecuación 5.5 sesimplifica:⎛ 2 2Q⎞⎜ r1+ 4a0s = ln⎟[5.6]2πT⎜ r ⎟1⎝ ⎠Para los punto situados sobre la línea perpendicular a la zona de recarga, cuando β=0 o 180ª, el cos (β) es iguala ± 1, y la expresión se simplifica:

36 CAPÍTULO 8 — HIDRÁULICA DE POZOSSuperficie del TerrenoNivelFreáticoBarrera ImpermeableAcuíferoMaterialImpermeableFigura 16. Acuífero limitado por una barrera impermeable. Quintero,1994AcuíferoNivel FreáticoMaterialImpermeableFigura 17. Acuífero limitado por dos barreras impermeables. Quintero, 1994Nivel FreáticoCorrienteAcuíferoMaterial ImpermeableFigura 18. Acuífero limitado por una zona de recarga. Quintero, 1994

CAPÍTULO 8 — HIDRÁULICA DE POZOS 37s =s =s =Q⎛ 2 2r 4a 4ar0ln⎜1+ ±2πT⎜ r1⎝Q⎛0ln⎜2πT⎜⎝Q⎛0ln⎜1+2πT⎜ r1⎝2( r − 2a) ⎞Q ⎛ ( r − 2a)1r11⎞⎟⎟⎠⎟ 0= ln2 T⎜⎟ π⎠ ⎝2( r 2a) ⎞Q ⎛ ( r + 2a)⎟ 0= ln2 T⎜⎟ π⎠ ⎝11rr11⎞⎟,β = 0⎠⎞⎟,β = 180º⎠[5.7]De estas ecuaciones se puede concluir que la pendiente de la curva de abatimiento de la parte que queda haciael río es más fuerte que la que va tierra adentro.5.2.2 Pozo construido en un acuífero que está limitado por una barreraimpermeableyP (x,y)r 2r 1Pozo 22axPozo 1Figura 19. Esquema de la ubicación de dos pozos. Quintero, 1994En la Figura 19, se representa un pozo construido en un acuífero que está limitado por una barrera impermeabley la cual no puede contribuir al bombeo, por lo tanto cuando el cono de abatimiento alcanza la barreraimpermeable y ante la imposibilidad de extenderse más allá de este límite se produce una caída más aceleradade la curva de abatimiento. Como se estudió en anterior numeral, el efecto que producen los pozos separadosuna distancia 2ª, sobre la línea que los divide, es que el abatimiento no varía con la distancia, y por lo tanto lalínea divisoria se comporta como impermeable. Así también se puede decir que el sistema analizado, esequivalente a dos pozos de descarga, funcionando en un acuífero infinito. La solución está dada por:Q ⎛ R ⎞2s = ln⎜⎟[5.8]2πT⎝ r1r2⎠En donde r 2 es la distancia desde el pozo imagen al punto, considerando r 1 la distancia del punto considerado alpozo de bombeo5.2.3 Ley de los tiemposCuando se tiene un piezómetro de monitoréo a una distancia r 0 del pozo de bombeo y sobre la líneaperpendicular a la barrera impermeable. Como en los dos casos anteriores, el sistema es equivalente almostrado en la Figura 19.El tiempo a partir del comienzo del bombeo para el cual se siente algún abatimiento en el piezómetro demonitoréo es cuando:

38 CAPÍTULO 8 — HIDRÁULICA DE POZOS2.25Tt20rSt0S= 1.0 ∴ =2[5.9]r 2.25 T0El tiempo que se necesita para que el pozo imagen tenga alguna influencia en el de monitoréo es cuando2.25Tt20rStiS= 1.0 ∴ =2[5.10]r 2.25 TiIgualando las dos expresiones, se tiene que:t020rt= [5.11]ri2iQue se conoce como la ley de los tiempos de Ingersoll (Quintero, 1994) en donde:t 0:r 0:t 1:r 1Tiempo transcurrido desde el comienzo del bombeo hasta que comienza a sentirse el abatimiento en elpozo de monitoréo.Distancia desde el pozo de monitoréo la pozo de bombeo.Tiempo a partir del cual existe una influencia del pozo imagen (o sea de la berrera impermeable)Distancia desde el pozo imagen al pozo de observación.Si el sistema está compuesto de varios pozos situados en una cierta distancia de la zona de recarga, el problemase resuelve como en los casos anteriores utilizando el método de la imágenes y el principio de superposición.6 APLIICACIIONESLa principal aplicación de la hidráulica de pozos está en determinar las características hidrogeológicas delacuífero, mediante el análisis de pruebas de bombeo, tema que discutirá en el próximo capítulo. Se ilustrará aquíla aplicación práctica de ñas diversas ecuaciones desarrolladas mediante los siguientes ejemplos.6.1 USO DE LA ECUACIÓN DE THEISEn una formación acuífera que tiene un espesor promedio de 12 m, una transmisividad de 8.64 m/d y uncoeficiente de almacenamiento de 0.001. El caudal de producción es de 4 L/s. Se necesita conocer elabatimiento a una distancia de 25 m, 8 horas después de comenzar la extracción de agua.El primer paso es determinar el parámetro u:22r S ( 25m) ( 0.001)-3u = == 4.5 x 10 [Adimensional]4Tt 4( 8.64 d)( )( )m 8 h 1 d24 hA continuación se determina la función del pozo de Theis W(u), mediante la curva de Theis (Figura 2) o en tablas,la cual en este punto posee un valor de 4.83.Ahora con estos valores sólo resta aplicar la ecuación de Theis, para encontrar el abatimiento.

CAPÍTULO 8 — HIDRÁULICA DE POZOS 39ssQ4 π T( r, t) = W( u)( 25 m , 8 h)=4 π, T = Kb0.004s( )( )( ) ( ) 8.64 12 m4.83 1.28 mm 1 ddm 3 =86400 s6.2 USO DE LA ECUACIÓN DE JACOBA manera de comparación se puede resolver el mismo ejemplo que fue resuelto con la ecuación de Theis, en elnumeral anterior.⎡⎛ ⎞⎤sQ ⎡ 2.25 Tt ⎤ln =24πT⎢ r S ⎥⎣ ⎦( r, t) ≅ ( 0.26 m)( 4.82) = 1.28m⎢ 2.25 8.64⎢ln⎢⎢⎣(m)( 12 m)4Ldss( r, t)≅24π( 8.64md)( 12 m)El Radio de Influencia es igual a:⎜⎝824( 25 m) ( 0.001)día⎟⎥⎠⎥⎥⎥⎦R12⎛ Tt ⎞= 1.5⎜⎟⎝ S ⎠⎡⎢= 1.5⎢⎢⎢⎣⎛( 8.64md)( 12 m)⎜⎝( 0.001)824⎞⎤día⎟⎥⎠⎥⎥⎥⎦12= 279 m6.3 USO DE LA ECUACIÓN DE CHENEn el acuífero descrito en el ejemplo anterior se tiene una frontera exterior a aproximadamente 80 m de distancia.Determinar el período durante el cual el acuífero puede ser analizado como un acuífero infinito.2R St ≤16T=162( 80 m) ( 0.001)( 8.64m)( 12 m)d= 38.5 x 10−3días = 5.56 minFísicamente, esto significa que el acuífero puede ser analizado como un acuífero infinito duranteaproximadamente cinco minutos y medio, desde que comienza el bombeo.6.4 USO DE LA ECUACIÓN DE PAPADOPULOS & COOPEREn un acuífero se excava un pozo de 0.2 m de radio, en toda su profundidad, (r w =r c ). Se pretende obtener uncaudal de 432 m 3 /día. La formación acuífera tiene una transmisividad de 86.4 m 2 /día y una capacidad dealmacenamiento de 0.01. Determinar el abatimiento después de 1 hora después de comenzar el bombeo, a unadistancia de 20 m.Usando la ecuación de Papadopulos y Cooper, se tiene que:

40 CAPÍTULO 8 — HIDRÁULICA DE POZOSssQ4πT( r, t) = F( u, α,ρ)( 20 m,60 min)432=4π22( 20 m) ( 0.01)rwS4u = == = 0.2724Tt 4( 86.4m)( 1 d)14.4día2422r ( 0.2 m)w 2S0.01 122rc( 0.2 m)10 −α = = × = ×r 20 mρ = = = 100r 0.2 mm3w4321085.7( ) ( ) día−2F 0.27,1 × 10 ,100 = ( 1.2) m = 0.397m × ( 1.2)=86.40.48 m2mdía6.5 USO DE LA ECUACIÓN DE THIEMA manera de ejemplo, se puede suponer un acuífero con un espesor de 6 m y una conductividad hidráulica de 10 -4m/s, encontrar el abatimiento a una distancia de 100 m del pozo, si el radio del pozo es de 0.1 m y la descargaes de 5 L/s.s = h rs =2 π() − h( r )s = 9.16 mw0.005Q ⎡ ⎛ r ⎞⎤= ⎢ln⎥2 T⎜⎣ r⎟π ⎝ w ⎠⎦ms( 6 m)( 0.0001m)6.6 USO DE LA ECUACIÓN DE DE GLEE - JACOB3s⎡ ⎛100 m ⎞⎤⎢ln⎜⎟⎥⎣ ⎝ 0.1m ⎠⎦Tomando el ejemplo de un acuífero confinado con espesor de 10 m, conductividad de 0.0002 m/s, y con unacapa semiconfinante de espesor 4 m y una conductividad hidráulica de 2 x 10 -8 m/s; se pìde determinar elabatimiento en condiciones estables a 400 metros de distancia de un pozo de 0.2 m de diámetro con unadescarga de 5 L/s, y en la pared del mismo.Para usar la ecuación simplificada de De Glee - Jacob, se necesita conocer el factor B:B12 ⎡= =− 8⎡b⋅ b' ⋅K⎤⎢ K' ⎥⎣ ⎦⎢⎣( 10 m)( 4 m)( 0.0002m)2x10Ahora si se reemplaza en la ecuación de De Glee - Jacob, y se llega a:0.005 s ⎛ 400 m ⎞s =K0⎜⎟ = 0.39 m2π( 0.0002ms)( 10 m)⎝ 632.45m ⎠mss12⎤⎥⎦= 632.45 m( )(.7397) 0.29 mm 3 =El valor de la función de Bessel se encuentra en tablas, en libros de Cálculo

CAPÍTULO 8 — HIDRÁULICA DE POZOS 41Y en la pared del pozo es de:0.005 s ⎛ 0.2 m ⎞s =K0⎜ ⎟ =2π( 0.0002ms)( 10 m)⎝ 632.45m ⎠( 0.39m)( 2.303) 3.26 mm 3 =6.7 USO DE LA ECUACIÓN DE DUPUIT - FORCHHEIMERPara un acuífero libre con H = 20m, radio del pozo 0.1 m y K = 0.0003 m/s, determinar la tasa máxima debombeo para crear un abatimiento de 2 m.Aplicando la ecuación de descarga de Dupuit – Fochheimer, se llega a:2( )2( 0.0003m)( 20 m) − ( 18 m)πs0.0716 3Q = = m . Donde R es el radio de influencia del pozo.⎛ R ⎞ln( 10R[ m]) sln⎜⎟⎝ 0.1 m ⎠7 REFERENCIIASBATU, Vedat. AQUIFER HYDRAULICS. John Wiley & Sons, Inc. USA. 1998.PISKUNOV, N. CÁLCULO INTEGRAL Y DIFERENCIAL. Editorial Mir. Moscú, Rusia. 1977.QUINTERO SAGRE, Jorge. HIDRÁULICA DE POZOS. Curso internacional de manejo y protección deacuíferos. Universidad Nacional de Colombia. Santafé de Bogotá. Agosto de1994.