Artículos técnicos

 

Relación gasto específico y transmisividad con flujo no lineal y pozo de penetración parcial

 

Relation between Specific Capacity and Transmissivity with Non-linear Flow and Partial Penetration Well

 

Armando O. Hernández-Valdés*
Instituto Superior Politécnico José Antonio Echeverría, Cuba.
*Autor de correspondencia.

 

Dirección institucional del autor

Dr. Armando O. Hernández-Valdés
Facultad de Ingeniería Civil
Centro de Investigaciones Hidráulicas
Instituto Superior Politécnico José Antonio Echeverría
Calle 114, núm. 11901 entre 119 y 127, Marianao
19390 La Habana, CUBA
Teléfonos: +53 (7) 2601 416 y 2603 636
Fax: +53 (7) 2672 013
ahernandez@cih.cujae.edu.cu

 

Recibido: 27/03/12
Aceptado: 22/04/13

 

Resumen

El gasto específico ha sido utilizado por numerosos autores, incluido el del presente trabajo, para estimar la transmisividad darciana del acuífero, dada la carencia de datos de estos sistemas y de los costos que implican la realización de ensayos de bombeo para determinarla. Existen numerosos factores que pueden incidir en errores significativos de dicha estimación. En este trabajo se presentan expresiones analíticas y resultados de ensayos de bombeo para demostrar las magnitudes de los errores, sobre todo cuando hay presencia de flujo no lineal en los alrededores del pozo de bombeo, si es de penetración parcial o hay efectos de recarga o cavernas que provocan en los pozos una rápida estabilización de los niveles durante el bombeo. En los dos primeros casos, las estimaciones subvaloran la transmisividad darciana del acuífero y en el último se sobrevalora este parámetro. Se muestran los resultados de una prueba de bombeo con dos satélites o pozos de observación, situación poco frecuente, pero que permite realizar la interpretación utilizando una de las variantes propuestas anteriormente para calcular las propiedades hidrogeológicas con régimen no lineal de flujo, donde se evidencia el error que se comete en la determinación del coeficiente de almacenamiento y el uso de un inadecuado gasto específico para estimar la transmisividad.

Palabras clave: régimen no lineal, penetración parcial, gasto específico-transmisividad, hidráulica de pozos.

 

Abstract

Specific capacity has been used by several authors, including the author of the present work, to estimate the Darcian transmissivity of an aquifer, given the lack of data related to these systems because of the costs involved in conducting pumping tests to determine transmissivity. Several factors can cause significant errors in this estimate. This work presents the analytical expressions and results from pumping tests in order to demonstrate the error magnitudes, especially when non-linear fluid is present around pumping wells, in the case of a partial penetration well, or there are recharge effects or caves that quickly stabilize the levels of the wells during pumping. In the first two cases the Darcian transmissivity of the aquifer is underestimated, whereas in the latter it is overestimated. Results from a pumping test using two satellites, or observation wells are also shown. Although this is an uncommon situation, interpretation is possible based on one of the proposals above to calculate hydrogeological properties under a non-linear flow regime. Errors in determining the storage coefficient and the use of an inadequate specific capacity parameter to estimate transmissivity were found.

Keywords: Non linear regime, partial penetration, specific capacity-transmissivity, well hydraulics.

 

Introducción

Numerosos autores han utilizado con frecuencia los datos del gasto específico (Qs) (relación entre el caudal de bombeo y el abatimiento o descenso que se produce en éste) para estimar la transmisividad del acuífero, dada la carencia de pruebas de bombeo que permitan caracterizar con mayor detalle a los sistemas acuíferos. La principal causa que motiva lo anterior es el costo de los ensayos de bombeo y la facilidad con que se determina el Qs durante el desarrollo del pozo o cuando se afora.

Aunque el Qs no sea un parámetro constante para cada pozo de bombeo ni pueda utilizarse para estimar la transmisividad del acuífero, según se verá más adelante, en cierta medida permite realizar una evaluación cualitativa de la calidad constructiva del pozo y de las características acuíferas de la formación geológica donde se encuentra ubicado. Por otra parte, si se quieren aplicar técnicas de optimización donde las restricciones de niveles mínimos se impongan en los pozos de bombeo, como es necesario hacer en los acuíferos costeros, entonces se puede aplicar una corrección al abatimiento simulado en puntos correspondientes a los pozos de bombeo, procedimiento utilizado por Cabrera-Estupiñán (2010) en función del gasto específico.

Pérez (2001) señala: "el gasto específico es función del tiempo y del caudal, y que solamente se podría considerar constante en el caso de flujo lineal en todo el acuífero, en un pozo sin estructura, caso que resulta poco probable, ya que en aquellos acuíferos donde el flujo pudiera ser lineal en la cara del pozo, es casi seguro que resulta necesario que el pozo tenga una estructura, al menos, formada por camisa y rejilla. Por otra parte, en los casos en que no se necesite estructura en el pozo, lo más probable es que el flujo sea no lineal en la cara del pozo. Es decir, que la probabilidad de que el caudal específico sea constante es muy pequeña. No obstante, se puede utilizar como índice de fertilidad relativa, aunque no con la absolutez que se ha usado hasta ahora".

Por lo general, se utiliza el gasto específico para caracterizar de manera cuantitativa las características de los pozos y de los acuíferos en diferentes cuencas hidrogeológicas, pero sin asociarlo con las propiedades hidrogeológicas de los acuíferos correspondientes, como reportan Capucci et al. (2001).

En adelante se abordarán los factores que inciden en las magnitudes del Qs y cómo inciden en la estimación de la transmisividad a partir de él, aunque hay trabajos, como el de Chenini, Silvain y Ben-Mammou (2008) que sólo consideran la variación con el tiempo debido al aumento del radio de influencia.

 

El gasto específico y factores que lo afectan

Son muchos los esfuerzos que se han realizado para tratar de buscar correlaciones entre la transmisividad T_D y el gasto específico Q_S, definido por la relación entre el caudal de bombeo y el abatimiento que se produce en dicho pozo.

Según señalan Pérez (2001) y Walton (1970), el abatimiento en los pozos de bombeo y, por ende, el gasto específico (Q_s) depende de varios factores:

• Tiempo de bombeo (A).

• Flujo no lineal en las proximidades del pozo de bombeo (B).

• Pérdidas por la estructura del pozo de bombeo (D).

Además de los anteriores factores se debe incluir:

• Penetración parcial.

• Superficie de filtración en acuíferos libres.

• Efectos de interferencia entre pozos y fronteras impermeables o de recarga.

• Heterogeneidad y anisotropía del medio geológico, etcétera.

Aunque es cierto lo que señalan Aitchison-Earl y Smith (2008) de que el gasto específico se reduce con el incremento del tiempo y caudal de bombeo, debe especificarse que el efecto del tiempo de bombeo sólo se refleja en la componente lineal del abatimiento, siendo más significativa su variación al principio y menor para después de varias horas de bombeo, donde ya sea aplicable la aproximación de Jacob, u ≤ 0.01, donde siendo E el coeficiente de almacenamiento, t el tiempo de bombeo, r el radio del pozo y T_D la transmisividad del acuífero. En relación con los incrementos de caudal, es importante destacar que tienen efectos importantes en el gasto específico, pues influyen tanto en la componente turbulenta del abatimiento como en las pérdidas por estructura y los efectos de la penetración parcial.

En el estudio realizado por Knopman y Hollyday (1993) es posible constatar la gran cantidad de variables que pueden influir en el gasto específico y realizan una valoración del comportamiento estadístico de dichas variables en un número considerable de pozos con datos de gasto específico, pero sin llegar a proponer expresiones matemáticas que relacionen éste con la transmisividad.

Lari, Knochenmus y Bowman (1998) señalan que en la zona de estudio, la gran variabilidad de la relación entre transmisividad y gasto específico puede ser causada por las pérdidas por estructura del pozo, la penetración de éste en el acuífero, la heterogeneidad del acuífero y la influencia de las fracturas en su interconexión.

Por su parte, El-Naga (1994) justifica que una correlación logarítmica entre transmisividad y gasto específico tiene un mejor ajuste que una lineal debido a que ambas son distribuciones lognormales y que la transmisividad de la matriz acuífera está afectada por los fenómenos de fracturación y carstificación. Como se demostrará más adelante, una buena correlación no necesariamente implica una buena estimación de un parámetro.

En el WRI Report 87-4034 (2008) se señalan las grandes limitaciones de estimar la transmisividad a partir del gasto específico, considerando que el radio efectivo del pozo de bombeo puede ser sustancialmente muy superior al nominal, ya que muchos pozos en formaciones cársticas pueden estar perforados en cavernas, lo que resultaría un alto gasto específico y por tanto estimarse altas transmisividades; por otra parte, la presencia del flujo turbulento y las pérdidas por estructura del pozo de bombeo pueden causar mayores abatimientos, al igual que la penetración parcial, reduciendo ambos parámetros. Se señala que los efectos anteriores pueden compensar los errores de estimación, cosa que es cierta para evaluar valores medios de la transmisividad, pero no para estimar la variabilidad espacial. Además se señala que en muchos casos el periodo de bombeo no se conoce y que virtualmente todo el abatimiento ocurre en pocos minutos.

En el presente trabajo se analizarán los efectos del régimen no lineal de flujo en las proximidades de los pozos de bombeo y el efecto de la penetración parcial del pozo en el acúifero, lo que hace coincidir con lo planteado en el reporte anterior y se demuestran teórica y prácticamente varios de dichos señalamientos.

 

Relación entre gasto específico y transmisividad T_D con régimen lineal alrededor del pozo de bombeo

a) Pozo de penetración total

A partir de la fórmula conocida de Dupuit-Thiem:

Se obtiene de despejar T_D y realizando el cambio de unidades:

La ecuación anterior pudiera expresarse en otras unidades, como las utilizadas en muchos de los trabajos a que se hacen referencia, de la siguiente manera:

Para un pozo de bombeo de 50 cm de diámetro y un radio de influencia de 370 m, la ecuación (3) se convertiría en:

Frecuentemente se ha utilizado en Cuba la expresión anterior para estimar la transmisividad hidráulica de Darcy, aunque entre otros factores depende del tiempo de bombeo.

Para los anteriores valores de diámetro y radio de influencia en la ecuación (3) se obtendría:

Para un radio de influencia de r_o = 1 250 m se obtendría de las dos ecuaciones anteriores:

De lo anterior se evidencia que el tiempo de bombeo afecta al gasto específico, disminuyendo a medida que el tiempo aumenta, pero que su variación es más significativa para tiempos pequeños de bombeo.

En las ecuaciones (4a) y (5a), el coeficiente obtenido está en el rango de las fórmulas utilizadas por Abbott y Dehay (2008), quienes emplearon la fórmula de Parkhurst, Christenson y Breit (1996= T = 242 Q/s, y señalan que los valores estimados por la relación anterior son de 1.1 a 4.3 veces mayores que los valores de las transmisividades obtenidas a partir de los ensayos de bombeo.

En el WRI Report 87-4034 (2008) se utiliza la relación T = 270 Q/s, indicando que para esa zona el coeficiente está entre 170 y 370, con un valor promedio de 270. Se señala que este coeficiente es muy variable y que acuíferos con cavernas como el estudiado provocan altos gastos específicos que pueden inducir a estimar transmisividades muy grandes y por lo contrario, el efecto del flujo turbulento estima transmisividades muy bajas; ambos aspectos serán analizados y justificados analíticamente en el presente trabajo.

 

b) Pozo de penetración parcial

De las expresiones que relacionan la transmisividad T_D con el gasto específico Q_S en pozos de penetración parcial con régimen lineal se tienen las tres primeras citadas por Pérez (2001):

En las expresiones anteriores, Fp = hs/m, α = r_o/r_p y x = hs/2r_p , siendo hs la penetración del pozo dentro del acuífero (m); m, el espesor del acuífero; r_p, el radio del pozo de bombeo, y r el radio de influencia.

En el cuadro 1 aparecen los resultados de evaluar las anteriores expresiones para la mismas magnitudes de los parámetros, obteniéndose resultados similares, salvo el caso de Babouchkine.

Tanto el efecto de la penetración parcial de los pozos de bombeo como el flujo no lineal que se puede producir alrededor de ellos provoca descensos superiores a los producidos si el régimen fuese lineal y el pozo de penetración total. En la simulación de los efectos locales en un modelo regional, ambos factores pueden ser simulados mediante un refinamiento de malla y medio heterogéneo (Cabrera-Estupiñán & Hernández-Valdés, 2011).

 

Relación entre gasto específico y transmisividad hidráulica de Darcy con régimen no lineal alrededor del pozo de bombeo

Se presentan de manera frecuente dos situaciones:

a) Pozo de penetración total

Para este análisis se considera importante introducir las expresiones que determinan el abatimiento a cualquier distancia de un pozo que bombea a caudal constante en un medio hiomogéneo e isótropo (Pérez, 2001):

O su equivalente:

Donde y se conoce como radio de influencia.

La distancia medida desde el centro del pozo de bombeo a partir de la cual puede considerarse el régimen lineal de flujo se determina por el radio de Darcy (r_D), definido por Pérez (2001):

La ecuación (6a) para el pozo de bombeo en función del r_D se expresa de la manera siguiente:

En las expresiones (6), (6a) y (8), el primer sumando de la derecha representa a la componente darciana del abatimiento y el segundo a la componente turbulenta, siendo Q el caudal de bombeo; S, el abatimiento; T_D = mK_D, la transmisividad darciana; T_T = mK_T, la transmisividad turbulenta, donde m es el espesor del acuífero; E, el coeficiente de almacenamiento; t, el tiempo de bombeo, y r es la distancia desde el pozo de observación al centro del pozo de bombeo.

De las expresiones anteriores y a partir de las conductividades hidráulicas darciana y turbulenta y por las siguientes ecuaciones: donde v es la viscosidad cinemática del fluido y g la acelaración de la gravedad, se pudieran obtener las propiedades hidrogeológicas características de los acuíferos: permeabilidad intrínseca o geométrica (k), y rugosidad equivalente del medio poroso o fisurado (C).

La representación en un gráfico semilogarítmico S vs log(r) de la ecuación (5a) tendría un comportamiento como el indicado en la figura 1.

Para obtener el mismo abatimiento en el pozo de bombeo con régimen lineal de flujo, la pendiente del gráfico semilogarítmico de la figura 1 sería ΔS'_r10 > ΔS'_r10, lo que implica que la transmisividad correspondiente que simula el efecto del régimen no lineal es menor que la del acuífero T'_D < T_D.

De igualar (8) con (9) se obtiene la relación entre la transmisividad real del acuífero y la estimada a partir del gasto específico si hay componente no lineal de flujo

Si se despeja T_D en la ecuación (8) se obtiene:

Sustituyendo r_p = 0.25 m y para un radio de influencia de r_o = 1 250 m se obtendría:

Observe que para régimen lineal por lo que el Factor_RNL coincide con el del régimen lineal de la ecuación (5).

 

b) Pozo de penetración parcial

El uso de pozos penetración parcial es lo más frecuente que se presenta como situación real en los pozos de bombeo; así, los efectos del régimen no lineal se manifiestan con mayor frecuencia por inducirse mayores gradientes en las proximidades del pozo de bombeo, lo que implica mayores abatimientos en éstos.

En la figura 2 se muestra un esquema representativo de un pozo de penetración parcial utilizado por Hernández-Valdés (1984) para calcular el abatimiento en este tipo de pozo, considerando el régimen no lineal de flujo. Se realizó la hipótesis de suponer que la línea de corriente inferior tiene forma parabólica, con vértice en la pared del pozo e intercepto con el fondo impermeable del acuífero a una distancia de dos veces el espesor de éste, por lo que el espesor del acuífero varía en función de la penetración del pozo de bombeo desde la pared de éste hasta dos veces su espesor.

A partir del esquema de la figura 2, en dicho trabajo se obtiene por integración la ecuación que permite determinar el abatimiento en el pozo de bombeo, considerando un régimen no lineal de flujo, que sigue la ley binómica que relaciona la velocidad con el gradiente hidráulico (Pérez, 2001):

Entre el pozo de bombeo y la distancia de dos veces el espesor del acuífero (2 m), la ecuación de continuidad presenta la siguiente forma:

Donde U es la velocidad media aparente.

Considerando despejando U de la ecuación (14), sustituyendo en la ecuación (13) y mediante un proceso de integración, se puede determinar la diferencia de abatimientos entre el pozo de bombeo y la distancia de dos veces el espesor del acuífero mediante la siguiente expresión:

El espesor es constante a partir de la distancia de dos veces el espesor del acuífero (2 m), hasta el radio de influencia (r_o), por lo que el abatimiento en éste se determinaría por:

La ecuación (17) se obtiene de la integración de la ecuación (15) con la sustitución de la ecuación (16) y mediante el empleo de variables adimensionales.

Despejando T_D de la ecuación obtenida del abatimiento se obtiene:

donde: Fp = hs/m, β = r_D/r_p, α = r_o/r_p y x = hs/2r_p siendo hs la penetración del pozo dentro del acuífero (m); r_D, el radio de Darcy; m, el espesor del acuífero; r_p, el radio del pozo de bombeo, y r_o, el radio de influencia.

La ecuación (17) puede sustituirse por el producto del Q_S por un factor que considera el régimen no lineal y la penetración parcial (Factor_RNLpp).

Una expresión aproximada del Factor_RNLpp que se obtiene por regresión de la ecuación (18) sería la siguiente:

 

Para el caso de penetración total (F_p = 1), el Factor_RNLpp ≈ Factor_RNL en la ecuación (12).

La correlación entre el factor de la fórmula analítica (17) con el ajustado representado por la ecuación (20), para un amplio rango de valores de Fp = hs/m y β = r_D/r , se observa en la figura 3.

Si se divide el coeficiente que afecta al Qs en pozos de penetración parcial o con efectos del régimen no lineal dado por (Factor_RNLpp) entre el (Factor_RL), se obtendrían los resultados del cuadro 2, que representan al coeficiente por el que se multiplicaría la transmisividad obtenida a partir del gasto específico, suponiendo régimen lineal de flujo y pozo de penetración total, en función de la penetración del pozo (Fp) y de la presencia de un régimen no lineal de flujo.

 

Valoración de la aplicación de los resultados anteriores en acuíferos cársicos de la antigua provincia La Habana, Cuba

Los principales acuíferos que abastecen a la ciudad de La Habana están en formaciones cársicas con diversos grados de fisuramiento y cavernosidad. La aplicación de ecuaciones deducidas para medios no consolidados en acuíferos cársicos ha encontrado en la comunidad científica diversos cuestionamientos; no obstante, modelos conceptuales que consideran estas formaciones geológicas como medios porosos equivalentes han demostrado ser aceptables en la modelación matemática a escala regional y en el flujo hacia pozos a escala local. No obstante, en este último caso es importante considerar la componente del abatimiento, que por lo general se manifiesta en los alrededores del pozo de bombeo y se debe a un régimen no lineal que no cumple con la ley de Darcy.

Para evaluar el comportamiento de la relación entre gasto específico y transmisividad, se incluyen los resultados de un conjunto de pruebas de bombeo realizadas en acuíferos de la antigua provincia La Habana, en Cuba, y fundamentalmente del acuífero costero del sur de la Habana, obtenidos por Molina-Rosabal y Ferrás-Martínez (1989), y que se muestran en el cuadro 3.

Las muestras no seleccionadas (148) del citado estudio presentaban como principal característica que en los pozos de bombeo se producía una estabilización muy rápida (85, inmediata estabilización; 63 en menos de 15 minutos), típico de un efecto de recarga, que puede ser causado por el sistema de cavernamiento de esta región cársica; tales pruebas daban un gasto específico muy grande y si se pronosticara la transmisividad T_D del acuífero por la ecuación (5), se estaría sobrevalorando este parámetro.

En la figura 4 se puede apreciar el comportamiento del factor que relaciona la transmisividad con el gasto específico (T_D/Qs) en las 90 muestras seleccionadas. Sólo el 34% de la muestra tiene un factor en el rango corrientemente utilizado entre 100 y 140. El 58% de la muestra tiene un factor en el rango entre 80 y 160, con lo que se demuestra que no es confiable realizar estimaciones de la transmisividad a partir del gasto específico. Las desviaciones evidencian los efectos del régimen no lineal de flujo, las pérdidas por estructura o la penetración parcial de los pozos de bombeo, no pudiéndose determinar la causa específica sin otros detalles adicionales.

En el propio trabajo de referencia (Molina-Rosabal & Ferrás-Martínez, 1989), se incluyeron los datos de pozos que disponían de satélites para tratar de determinar la propiedad del acuífero, indicativo del régimen no lineal, utilizando el procedimiento propuesto por Pérez (2001) y modificado junto con otro procedimiento por Hernández-Valdés (2008). En este último trabajo, el autor propone que el pozo de bombeo no sea utilizado para el cálculo de propiedades, ya que se le pueden achacar al régimen no lineal y a la propiedad que lo caracteriza causas como la penetración parcial, pérdidas por estructuras u otras, que realmente no se deben a él.

De los ensayos que debajo se muestran, en cinco de ellos se utilizó el pozo de bombeo para determinar T_T, por lo que el radio de Darcy obtenido puede estar afectado por las causas antes indicadas; no obstante, se aprecia que la relación equivalente al factor, puede ser explicada en su magnitud por la presencia del régimen no lineal de flujo y/o por la presencia de la penetración parcial del pozo de bombeo, sin poder esclarecer la causa sin otras informaciones.

En el cuadro 4 se muestran los resultados de cinco ensayos de bombeo con pozos satélites tomados del trabajo de Molina-Rosabal y Ferrás-Martínez (1989); se agregó el 966 de la misma cuenca.

 

Comentarios sobre resultados de ensayos de bombeo con pozos satélites

En los casos donde se disponía de dos pozos de observación (pozo 15 y 966), se realizan los siguientes comentarios:

• En el primer caso no se manifestaba el régimen no lineal entre los dos satélites, por lo que se utilizó el pozo de bombeo y el más cercano, situado a 38 m; esto implicó que se consideraran pérdidas de carga como producto del régimen no lineal de flujo, siendo otras las causas y se cometiera un error en el cálculo de T_T, dando un radio de Darcy superior a las distancias a las que se encontraban los pozos de observación, donde no se manifestaba el régimen lineal, por lo que de acuerdo con la relación entre T_D/ Q_S, el valor del radio de Darcy debía ser teóricamente de 53 m, despejando del

• En el segundo caso se aprecia que el gasto específico obtenido de la prueba está afectado por efectos de recarga o cavernamiento, y que el abatimiento que debía obtenerse en el pozo de bombeo resultaría un gasto específico cuatro veces menor y que por ende existiría entonces correspondencia entre la relación y el Factor_RNL. En la figura 5 se ilustra el cálculo de propiedades con la metodología propuesta por Hernández-Valdés (2008) y se muestran graficamente los datos del ensayo de bombeo.

Esta prueba de bombeo se realizó con un caudal de 26 lps en un pozo de 40 cm de diámetro y se midieron los abatimientos en dos satélites situados a 9.7 y 15 m, respectivamente, del pozo de bombeo.

 

Cálculo de la transmisividad darciana T_D

Se sigue el mismo procedimiento enunciado por Cooper y Jacob (1946), citado por Todd (1959), obteniendo de la representación en un gráfico S vs. log(t), la pendiente por ciclo ΔSt ₁₀, y de aquí se determina la T_D por la siguiente expresión:

 

Cálculo de la transmisividad turbulenta T_T

Como las rectas ajustadas en los gráficos semilogarítmicos del inciso anterior deben ser paralelas, se determina la diferencia de abatimientos entre ellas para cualquier tiempo y se despeja de la ecuación (6a) el valor de T_T.

Siendo r₁ = 9.7 m y r₂ = 15 m distancias desde el pozo de bombeo a los pozos de observación; S₁ = 2.96 m y S₂ = 2.16 m, abatimientos correspondientes a las distancias anteriores para un tiempo de 239 minutos.

 

Cálculo del radio de Darcy r_D

Este parámetro ya se definió como y daría r_D = 224 m.

 

Determinación del coeficiente de almacenamiento E

Como los dos pozos de observación están en la zona de régimen no lineal (inferiores al radio de Darcy), entonces se sigue el procedimiento empleado por Jacob (Chenini et al., 2008) para calcular E_j, obteniéndose a partir de este valor el coeficiente de almacenamiento E:

Para r₁ = 9.7 m, t₀ = 3.47 *10^-4 días y E_j = 0.0031, por tanto E = 0.031.

Para r₁ = 15 m, t₀ = 1.8 *10^-3 días y E_j = 0.0069, por tanto E = 0.031.

El factor que afecta a E_j en la ecuación anterior es precisamente el error que se comete al determinar el coeficiente de almacenamiento E, por el procedimiento propuesto por Jacob en un pozo de observación que se encuentra en la zona no lineal de flujo.

 

Utilización del factor Factor_RNLpp para estimar abatimientos en los pozos de bombeo

De igual manera, se puede decir que el abatimiento en el pozo de bombeo puede ser determinado en función del caudal (Q) y la transmisividad (T_D), despejando de la ecuación (18) por:

En la propuesta de Hernández-Valdés y Llanusa-Ruiz (2009), aplicada por Cabrera-Estupiñán (2010), para considerar como restricción el abatimiento en algunos nodos representativos de pozos de bombeo si se conoce la penetración del pozo, se calcularía el Factor_RNLpp para el caudal de bombeo de máxima capacidad instalada, a fin de evitar iterar, dado que el radio de Darcy (7) depende de éste y estar del lado de la seguridad; ello permitiría utilizar el modelo de optimización MADA (Cabrera-Estupiñán & Dilla-Salvador, 2011), con un refinamiento semilocal de la malla en el modelo regional del acuífero.

 

Conclusiones

De este trabajo se pueden obtener las siguientes conclusiones:

• Las estimaciones de la transmisividad (T_D) a partir del gasto específico (Q_S) induce a errores significativos de este parámetro, ya que se puede sobrevalorar en los casos donde se logre una rápida estabilización de los niveles en los pozos de bombeo, o se puede subvalorar cuando efectos de penetración parcial, régimen no lineal, pérdidas por estructura del pozo, etcétera, estén influyendo en la magnitud del factor que los correlaciona.

• El efecto del régimen no lineal de flujo en los pozos de bombeo induce a calcular transmisividades menores que las reales de utilizarse el gasto específico para estimarlas y éstas serían cada vez menores a medida que aumenten los caudales de bombeo.

• La magnitud del coeficiente que relaciona T_D (m²/d) con Q_S (lps/m) puede indicar si hay efectos de recarga o cavernamiento, régimen no lineal de flujo, penetración parcial, pérdidas por estructura del pozo de bombeo, etcétera.

• Para identificar las verdaderas causas que originan los incrementos del abatimiento en los pozos de bombeo es necesario conocer las caracteristicas constructivas de la obra de captación, las propiedades hidrogeológicas del acuífero y las condiciones de explotación a que está siendo sometida dicha obra.

• Para diferentes caudales de bombeo, los abatimientos en los pozos de bombeo pueden ser estimados en pozos de penetración parcial y con régimen no lineal de flujo a partir del factor Factor_RNLpp y de la transmisividad T_D, conocidas las propiedades del acuífero y el factor de penetración.

• Se recomienda realizar ensayos de bombeo que permitan determinar las propiedades del acuífero con régimen no lineal, realizando dos ensayos de bombeo con caudales diferentes y un pozo de observación, o un ensayo si se dispone de dos pozos de observación próximos al de bombeo y al menos uno donde se manifieste el efecto del flujo no lineal. No es aconsejable utilizar los datos del pozo de bombeo, salvo para calcular la transmisividad T_D.

 

Referencias

Abbott, M. M., & Dehay, K. (2008). Aquifer Tests and Characterization of Transmissivity, Ada-Vamoosa Aquifer on the Osage Reservation (10 pp.). Osage County, USA: Geological Survey Scientific Investigations Report 20085118. Recuperado de pubs.usgs.gov/sir/2008/5118/.         [ Links ]

Aitchison-Earl, P., & Smith, M. (2008). Aquifer Test Guidelines (2nd edition). Report No. R08/25. Environment Canterbury Technical Report. Recuperado de www.ecan.govt.nz.         [ Links ]

Cabrera-Estupiñán, E. (2010). Sistema para la administración de la explotación del agua subterránea. Tesis presentada en opción al grado científico de Doctor en Ciencias Técnicas. La Habana: Instituto Superior Politécnico José Antonio Echeverría.         [ Links ]

Cabrera-Estupiñán, E., & Hernández-Valdés, A. O. (2011). Modelación del agua subterránea a escala regional con refinamiento local de la malla. Planteamiento y validación del algoritmo. Tecnología y Ciencias del Agua, antes Ingeniería Hidráulica en México, 2(1), 65-82.         [ Links ]

Cabrera-Estupinán, E., & Dilla-Salvador, F. (octubre-diciembre, 2011). Modelo de administración de acuíferos: MADA. Tecnología y Ciencias del Agua, antes Ingeniería hidráulica en México, 2(4), 5-24.         [ Links ]

Chenini, I., Silvain, R., & Ben-Mammou, A. (2008). A simple method to estimate Transmissibility and Storativity of Aquifer Using Specific Capacity of Wells. J. Applied Sci., 8(14), 2640-2643.         [ Links ]

De Glee, G. J. (1930). Over Groundwaterstominger Bij Water Onnttrekking Door Middle Van Put Ten. Thesis. Delft Holland: J. Waltman.         [ Links ]

El-Naga, A. (1994). Estimation of transmissivity from specific capacity data in fractured carbonate rock aquifer, central Jordan. Environmental Geology, 23(1), 73-80.         [ Links ]

Hernández-Valdés, A. O. (1984). Efectos de la penetración parcial e pozos con régimen no lineal. Ingeniería Hidráulica (ISPJAE), 5(2).         [ Links ]

Hernández-Valdés, A. O., & Llanusa-Ruiz, H. (2009). "La modelación matemática y la explotación de los acuíferos costeros". Congreso Internacional de Hidráulica. Cayo Santa María. Cuba.         [ Links ]

Hernández-Valdés, A. O. (2008). Nuevas formulaciones para determinar las propiedades hidrogeológicas con régimen no lineal. XXIII Congreso Latinoamericano de Hidráulica, Cartagena de Indias, Colombia.         [ Links ]

Knopman, D. S., & Hollyday, E. F. (1993). Variation in Specific Capacity in Fractured Rocks, Pennsylvania. Ground Water, 31(1), 135-145.         [ Links ]

Knochenmus, A., & Bowman, G. (MOE) (1998). Transmissivity and Water Quality of Water-Producing Zones in the Intermediate Aquifer System. Sarasota County, USA: Geological Survey Water-Resources Investigations Report 98-4091.         [ Links ]

Molina-Rosabal, Y., & Ferrás-Martínez, L. (1989). Características de las propiedades hidrogeológicas de la Cuenca Sur: Tramo Güira-Quivicán. Trabajo de Diploma. La Habana: Facultad de Ingeniería Civil, ISPJAE.         [ Links ]

Parkhurst, D. L., Christenson, S., & Breit, G. N. (1996). Ground-Water-Quality Assessment of the Central Oklahoma Aquifer, Oklahoma: Geochemical and Geohydrologic Investigations (101 pp.). Paper 2357-C. US Geological Survey Water-Supply.         [ Links ]

Pérez, D. (2001). La explotación del agua subterránea. Un nuevo enfoque. La Habana: Editorial Félix Varela.         [ Links ]

Schneebeli, G. (1966). Hydraulique Souterraine. Paris: Eyrolles.         [ Links ]

Walton, W. (1970). Groundwater Resource Evaluation (664 pp.). McGraw-Hill Series in Water Resources and Environmental Engineering.         [ Links ]

WRI Report 87-4034 (2008). Hydrogeology, Aquifer Characteristics, and Ground-Water Flow of the Surficial Aquifer System, Broward County. Disponible en: sofia.usgs.gov/publications/wri/87-4034/spcapacity.html.         [ Links ]