Servicios personalizados
Articulo en XML
Referencias del artículo
Enviar articulo por email
Citado por SciELO 
Similares en
SciELO 
Permalink
Introducción
Por lo general, las actividades de desarrollo en las estructuras de captación siempre han tenido orientación informativa, tanto técnica como experimental, en relación con los procesos comúnmente utilizados (mecánicos y químicos) en dichas actividades que, en la realidad, resultan bastante abstractos dentro de ese punto de vista (Johnson 1966). No obstante, aunque no se realizó una investigación teórica exhaustiva a los efectos producidos por las actividades de desarrollo, si se contempló que dicha actividad generaba diversos tipos de beneficios para obtener un rendimiento hidráulico eficiente del pozo (Vélez 2004).
Por lo tanto, en el congreso de Cambridge de 1985 (Iglesias 1989), la comunidad científica reveló su preocupación e interés por el estudio sobre las pérdidas de carga que están presentes en los pozos, ya que los estudios teóricos se estaban completando a través del tiempo.
Por lo anterior, uno de los avances relevantes en relación con la cuantificación teórica de las actividades de desarrollo en los pozos, surgió con los trabajos de Iglesias (1989), en los cuales, se dedujeron una serie de expresiones que describen el abatimiento, teniendo en cuenta que; una vez el pozo sea desarrollado, se formará una zona de desarenado más permeable que alcanzará una distancia con respecto al radio del pozo r p , n veces mayor.
Dentro de ese orden de ideas, Pérez Franco (1995) describe la imagen general de un pozo para un acuífero confinado de la siguiente manera:
Zona del empaque de gravas (entre r W y r p ).
Zona de desarenado (entre r p y r da ).
Zona del acuífero sin alteraciones (entre r da y r o ).
donde: r W , distancia desde el eje del pozo a la tubería de revestimiento o rejilla.
r p , distancia desde el eje del pozo hasta el límite exterior del empaque de grava.
r da , distancia desde el eje del pozo hasta el comienzo de la formación acuífera no alterada.
r o , radio de influencia.
Ahora bien, considerando que el flujo es no lineal, sin llegar a ser cuadrático o totalmente desarrollado, se presenta entre el radio turbulento r T y el radio de Darcy r D , donde r D indica el comienzo de la zona donde el régimen es lineal (Pérez Franco 1995) y resultaría probable que la zona de desarenado se encuentre entre el radio turbulento r T y el radio de Darcy r D .
No obstante, de la misma manera, es probable también que dicha zona se encuentre entre el radio del pozo r p y el radio turbulento r T , es decir, donde el flujo se encuentra totalmente desarrollado (turbulento puro). Por otra parte, lógicamente no es posible que dicha zona se desarrolle donde el flujo se puede considerar lineal, o sea, más allá de r D , cuando existe la componente no lineal en el abatimiento, a menos que el radio de Darcy sea menor que el radio del pozo (r D < r p ) (Pérez Franco 1995).
Consideraciones generales al analizar el flujo hacia el pozo
Para lograr un modelo simple que permita analizar las situaciones más frecuentes que se presentan en la realidad, se han supuesto a través del tiempo una serie de condiciones ideales que son las siguientes:
El acuífero tiene un espesor constante.
El pozo atraviesa todo el espesor del acuífero.
El acuífero es homogéneo, isotrópico y de extensión infinita.
El pozo solo recibe agua de todo el espesor del acuífero.
El pozo de extracción no tiene ningún tipo de estructura de revestimiento (filtros, rejillas o camisas), en la zona comprendida del espesor saturado del acuífero.
La superficie piezométrica del acuífero es horizontal, lo que quiere decir que no se tiene en cuenta el gradiente natural del flujo, o sea, las componentes verticales de la velocidad que, en los acuíferos freáticos, en ocasiones pueden ser considerables.
Las líneas de corriente son radiales y están sobre un plano unidireccional horizontal.
Por otra parte, para efectos prácticos se puede eliminar de la imagen general del pozo el empaque de gravas, ahora bien, dicha simplificación de la misma manera es asumida por diversos investigadores e, incluso, en el aporte realizado por Iglesias (1989), en relación con la cuantificación de las actividades de desarrollo en un pozo y su influencia sobre el flujo radial hacia este.
Función de desarenado y permeabilidad
El concepto de factor de mejora de la permeabilidad 
y, por consiguiente, de la función de permeabilidad 
, fue adoptado y desarrollado por Iglesias (1989), teniendo en cuenta el intervalo 
. Dichos desarrollos se pueden expresar en la notación del autor y añadiendo las equivalencias de la componente no lineal (Pérez Franco 1995), de la siguiente manera:
(1)
(2)
donde
μ |
coeficiente de uniformidad |
F |
fracción desarenada en el borde de la rejilla |
n |
nro. de veces donde se prolonga la zona de desarenado partiendo desde el radio del pozo |
rp |
radio del pozo |
KD, KDa |
conductividad hidráulica lineal y mejorada respectivamente |
KT, KTa |
conductividad hidráulica turbulenta y mejorada respectivamente |
TD, TDa |
transmisibilidad hidráulica lineal y mejorada respectivamente |
TT, TTa |
transmisibilidad hidráulica turbulenta y mejorada respectivamente |
Flujo no lineal permanente en un acuífero confinado y desarrollado
Antes que nada, considerando la generalización anterior sobre la imagen general de un pozo (figura 1) y denominando las variables asociadas a la zona desarenada como lo realizó Pérez Franco (1995), se tiene que:
Fig. 1 Imagen general de un pozo sin camisa y filtro artificial de grava en un acuífero confinado, considerando las actividades de desarrollo.
El abatimiento en cualquier punto del radio de influencia, utilizando la ley binómica de flujo, vendrá dado por una componente lineal y cuadrática (Pérez Franco 1995), por lo cual:
(3)
Ahora bien, como se sabe que el desarrollo en el pozo alcanza una distancia 
, si se integra a ambos lados entre los límites 
y se introduce el nuevo valor de la transmisibilidad mejorada 
, resulta que:
De hecho, si se introduce la expresión de la función de transmisibilidad 
y 
(2), se tiene que:
Es importante mencionar que Iglesias (1989), resolvió la expresión de la componente lineal en presencia de la mejora de la permeabilidad. Dicha expresión en la notación del autor vendrá de la siguiente manera:
donde, la función lineal del pozo desarrollado es:
(4)
Por otro lado, si se integra la componente cuadrática faltante, se obtiene lo siguiente:
donde, la función no lineal del pozo desarrollado es:
(5)
Por consiguiente, el abatimiento dentro del intervalo r p - r da , vendrá dado por:
(6)
Por otra parte, el abatimiento entre r da y r o , según Pérez Franco (1995), viene dado por:
(7)
Ahora, si se suman las ecuaciones (6) y (7), correspondientes al abatimiento producido en la zona desarenada y en la zona no alterada, se tiene que:
(8)
donde, según Iglesias (1989), la función de entorno E, será:
Flujo no lineal permanente en un acuífero libre y desarrollado
Considerando de la misma manera la imagen general de un pozo adoptada anteriormente y las condiciones adoptadas por Dupuit (1863), resulta que (figura 2):
Fig.2 Imagen general de un pozo sin camisa y filtro artificial de grava en un acuífero libre, considerando las actividades de desarrollo.
Según los trabajos de Huyakorn y Dudgeon (1976), el abatimiento S r , para diferentes condiciones de confinamiento, se compone de dos partes: una, denominada S D , equivalente a la expresión lineal del abatimiento y, la otra, denominada S ND , equivalente a la expresión no lineal del abatimiento. Expresando lo anterior en la notación deducida por Pérez Franco (1995), se tiene que:
(9)
Ahora bien, Iglesias (1989), encontró una expresión equivalente en régimen lineal en presencia de desarrollo, por lo que:
(10)
No obstante, según Pérez Franco (1995), para formular la ecuación (10) en función del abatimiento y el espesor saturado, se puede deducir con facilidad que:
(11)
Por lo tanto, si se introduce las relaciones anteriores (11) en el primer miembro de la ecuación (10) y, si se despeja S Dp , se tiene que:
(12)
donde, Pérez Franco (1995), corroboró que la relación 
, no era más que el espesor promedio del acuífero saturado, 
. Ecuación que tiene similitud con la componente lineal de la ecuación (8) para acuíferos confinados.
Por otro lado, si se generaliza el segundo miembro de la segunda adición de la ecuación (9), en relación con los abatimientos producidos en la zona de desarenado y en la zona no alterada dentro del radio de influencia en régimen turbulento puro, se tiene que:
O sea que, para el segundo miembro del primer sumando 
, si se introduce la expresión de la función de permeabilidad
(2), en régimen turbulento puro (Pérez Franco 1995), resulta que:
Si se integra la ecuación anterior entre los límites (hrda - hp y rda - rp), se tiene que:
(13)
Ahora bien, para el segundo miembro del segundo sumando ST, Pérez Franco (1995) obtuvo la siguiente expresión, ajustada por el autor para la zona comprendida entre rda y ro.
(14)
En este orden de ideas, si se suman las ecuaciones (13) y (14), resulta que:
(15)
Pero, si se ajusta la ecuación (11) en relación con el abatimiento , se obtiene lo siguiente:
(16)
Por lo tanto, si se introduce de nuevo las relaciones anteriores en el primer miembro de la ecuación (15) y, si se despeja , se tiene que:
(17)
Ahora, si se examina la ecuación (9) y se introducen los valores de 
(12) y
(17), se puede concluir lo siguiente:
(18)
La ecuación anterior no es producto de aproximaciones como se acostumbra a realizar (h o + h ≈ 2h o ), cuando se considera que el espesor saturado del acuífero libre es lo suficientemente grande en relación con los abatimientos observados.
Análisis de las expresiones equivalentes para los acuíferos confinados y libres
Como se sabe, las expresiones del abatimiento registrado en la cara del pozo, utilizando la ley binómica de flujo para acuíferos confinados y libres respectivamente, son:
(18)
Si se realiza un análisis somero a la ecuación (18), en relación con la componente lineal del acuífero libre, se puede observar, tal y como lo corroboró Pérez Franco (1995), que el término (h o - S Dp /2) resulta ser el espesor promedio del acuífero saturado h prD = (h o + h p ). En este orden de ideas, se puede corroborar lógicamente que el espesor promedio del acuífero saturado h prD para la componente lineal, difiere con el denominado como espesor aparente h prT , para la componente turbulenta, sin embargo, las expresiones obtenidas para los acuíferos confinados y libres resultan bastante similares.
Por lo anterior, se puede expresar ambas componentes de la siguiente manera:
Ahora, si se establece una relación entre las ecuaciones anteriores, resulta que (19):
(19)
Conclusiones
Antes que todo, es importante mencionar que las expresiones obtenidas en régimen permanente, utilizando la ley binómica de flujo para acuíferos confinados y libres, resultan bastante similares; tanto para la componente lineal como para la turbulenta. No obstante, cabe resaltar la mayor similitud en relación con la componente lineal para ambos casos, por lo tanto, es posible utilizar la expresión lineal de un acuífero confinado para un acuífero libre, sustituyendo el valor de m por h prD , tal y como lo sugirió Pérez Franco (1995) en su momento.
Por el contrario, para la componente cuadrática se observa una diferencia particular con respecto a la formulación obtenida para los acuíferos confinados, debido a la variación del espesor saturado del acuífero libre en relación con la distancia, teniendo en cuenta las suposiciones de Dupuit (1863).
Aunado a ello se puede enfatizar que las expresiones obtenidas para un acuífero libre no son producto de aproximaciones, como cuando se considera que el espesor promedio del acuífero saturado es lo suficientemente grande en relación con los abatimientos. Por lo anterior, se puede concluir que el espesor promedio del acuífero saturado h prD , cuando el régimen es lineal, difiere del denominado como espesor aparente h prT , cuando el régimen es turbulento puro en un acuífero libre.
Para finalizar, se puede mencionar que, con suma importancia, la función del pozo desarrollado ∅ D y ∅ T (lineal y no lineal) resulta ser diferente para ambos regímenes, además que, dicha función aparece de la misma manera para las formulaciones propuestas en los acuíferos confinados y libres.
