Estimación de parámetros de acuíferos
utilizando la recesión del flujo de base


Victor M. Ponce, Rosa D. Aguilar, y Sudhir Kumar


14 Junio 2013



1.  INTRODUCCIÓN

En la hidrología de agua subterránea hay necesidad de estimar los parámetros del acuífero, tales como el tiempo de almacenamiento, la constante de la cuenca, y la difusividad hidráulica.

Es posible estimar estos parámetros utilizando la recesión del flujo de base, mediante el uso de la fórmula de Rorabaugh y Simons (1966).

Rorabaugh (1963) convirtió la ecuación de difusión de calor a unidades de flujo de agua subterránea, haciendo posible de esta manera relacionar el flujo de base con los parámetros del acuífero.

En este trabajo hemos aplicado el modelo de Rorabaugh a la recesión del hidrograma en la cuenca del río Papaloapan, en el sur de México.

El objetivo es mostrar que se pueden utilizar datos de aforos de alta calidad para estimar parámetros del acuífero, ligando de esta manera los procesos de agua superficial y agua subterránea.

Estas estimaciones complementan aquéllas obtenidas mediante pruebas de bombeo, pudiendo eventualmente ser usadas en ausencia de éstas.


2.  EL MODELO DE RORABAUGH

Rorabagh (1963) desarrolló un modelo teórico del flujo de agua subterránea a una corriente, asumiendo una cuenca homogénea, isotrópica, y uniforme, es decir, con conductividad hidráulica, coeficiente de almacenamiento, y espesor del acuífero constantes.

La ecuación es:

q  =  2T (ho /a) e - (π 2 T t ) / (4 a 2 S )
(1)

en la cual

q = descarga de agua subterránea, por unidad de longitud de la corriente, de un lado solamente;

T = transmisibilidad, definida como el producto de la conductividad hidráulica y el espesor del acuífero;

ho = carga hidráulica instantánea entre la superfice del agua en el acuífero y la corriente, medida desde la divisoria de agua en la cuenca;

a = distancia de la corriente a la divisoria de la cuenca (la mitad del ancho del acuífero);

S = coeficiente de almacenamiento, y t = tiempo.

Esta ecuación es válida después que se ha establecido la recesión del hidrograma, es decir, después de un tiempo crítico tc de tal manera que (T tc ) / (a 2 S ) = 0.2.

Utilizando la ecuación de Darcy, se puede demostrar que la cantidad 2T(ho /a) es esencialmente una descarga por unidad de ancho qo.

Como la Ecuación 1 modela la atenuación exponencial de la descarga de agua subterránea, Rorabaugh y Simons extendieron esta ecuación para obtener su modelo de flujo de base en una estación:

Q  =  Qo e - (π 2 T t ) / (4 a 2 S )
(2)

en la cual

Q = flujo de base en el tiempo t; y

Qo = flujo de base en el tiempo t = 0.

En este estudio aplicamos el modelo de Rorabaugh para estimar las características del acuífero en una cuenca tropical grande, caracterizada por una gran diversidad de escenarios hidrogeológicos y climáticos.


3.  CARACTERÍSTICAS DEL ACUÍFERO

Boussinesq (1877) linearizó la ecuación de flujo de agua subterránea, expresándola en forma de una ecuación de difusión, la cual puede ser resuelta en forma más conveniente.

Con la analogía de difusión, una cuenca de agua subterránea (manto acuífero) puede ser caracterizada en términos de los siguientes parámetros:

1. Tiempo de almacenamiento ts, una constante de recesión igual a:

           4 a 2 S
ts  =  ________

            π 2 T
(3)

El tiempo de almacenamiento es tal que cuando t = ts, la descarga se ha reducido al 37% del valor correspondiente a t = 0, es decir, es una medida de la velocidad relativa de la recesión del hidrograma.

2. Constante de cuenca Kb:

                T
Kb  =  _______

            a 2 S
(4)

La constante de cuenca Kb (en unidades de T -1) combina la propiedades geométricas e hidrogeológicas del acuífero.

3. Difusividad hidráulica D:

            T
D  =   ____

           S
(5)

La difusividad hidráulica (en unidades L2 T -1) combina las propiedades hidrogeológicas del acuífero en un parámetro que caracteriza el proceso de difusión.


4.  APLICACIÓN A LA CUENCA DEL RÍO PAPALOAPAN

El modelo de Rorabaugh fue aplicado a la recesión del hidrograma en la cuenca del río Papaloapan, en los estados de Veracruz y Oaxaca, México (Fig. 1).

Esta cuenca es extensa, de 46,517 km2, con una gran variedad de escenarios hidrogeológicos y climáticos.

La cuenca del río Papaloapan, en los estados de Veracruz y Oaxaca, México

Fig. 1  La cuenca del río Papaloapan, en los estados de Veracruz y Oaxaca, México.

El Cuadro 1 muestra diez (10) estaciones de aforo seleccionadas, nombre de la corriente o río, localización geográfica, precipitación media anual y propiedades geométricas de la cuenca.

Cuadro 1.  Estaciones de aforo y propiedades geométricas en la cuenca del río Papaloapan.
Estación de aforo
Corriente
o río
Latitud
Longitud
Precipitación
media anual 1
(mm)
Area de
drenaje

(km2)
Longitud
hidráulica
(km)
Ancho del
acuífero
(km)
Achotal La Trinidad 17° 46' 95° 09' 1620 2333 124.4 18.7
Angel R. Cabadas Tecolapa 18° 35' 95° 26' 2255 125 24.4 5.1
Azueta Tesechoacan 18° 05' 95° 43' 1533 1656 83.7 19.8
Bellaco Lalana 17° 46' 95° 11' 1587 2917 115.5 25.3
Cuatotolapan San Juan 18° 09' 95° 18' 1305 7090 167.4 42.3
Jacatepec Valle Nacional 17° 52' 96° 12' 3906 1117 58.3 19.2
La Estrella Usila 17° 55' 96° 26' 4805 774 34.1 22.7
Monterrosa Cajones 17° 48' 95° 56' 2288 2870 111.6 25.7
Qiotepec Grande 17° 54' 96° 59' 636 4832 76.3 63.3
Xiquila Xiquila 18° 02' 97° 09' 354 1073 55.1 19.5
1 Isohyetas anuales de México, 1931-90, Regiones 28 (Papaloapan) y 29 (Coatzacoalcos).

La longitud hidráulica, es decir, la longitud medida a lo largo de la corriente principal, se obtuvo de mapas topográficos.

El ancho del acuífero se estimó como la relación de área de drenaje dividida por la longitud hidráulica.

Las condiciones climáticas varían de hiperhúmedo a árido; por ejemplo, río Usila en La Estrella, con 4,805 mm de precipitación media anual, al río Xiquila en Xiquila, con 354 mm.

La Figura 2 muestra un hidrograma típico de la cuenca del río Papaloapan, el del río Tesechoacán en Azueta (Comisión del Papaloapan 1972).

Río Tesechoacán en Azueta

Fig. 2  Río Tesechoacán en Azueta (Comisión del Papaloapan 1972).

Para las diez estaciones seleccionadas, se coleccionaron los datos de recesión del hidrograma para el período de dos años 1971-72.

Para cada estación, se seleccionaron varios períodos de recesión del hidrograma, con duraciones de 5 a 18 días.

Para cada período, se calcularon valores medios diarios de:

(a) tiempo de almacenamiento,

(b) constante de la cuenca, y

(c) difusividad hidráulica.

Para cada estación, los valores medios de cada período se utilizaron para calcular un valor promedio correspondiente a la estación.

El Cuadro 2 muestra los parámeteros de acuífero calculados para las diez (10) estaciones seleccionadas, las cuales comprenden una gran diversidad de escenarios hidrogeológicos y climáticos.

El tiempo de almacenamiento varía entre 23.49 and 116.64 días; la constante de la cuenca varía entre 0.00347 and 0.01725 días -1; la difusividad hidráulica varía entre 0.058 and 15.619 km2 día -1.

Cuadro 2.  Características de los acuíferos en la cuenca del río Papaloapan.
Estación de aforo
Corriente
o río
Tiempo de
almacenamiento1
Constante
de la cuenca

Difusividad
hidráulica
__
x
(días)
s
(días)
Cv
días-1
km2 día-1
m2 seg-1
Achotal La Trinidad 51.28 8.51 0.166 0.00790 0.694 8.037
Angel R. Cabadas Tecolapa 46.03 9.30 0.202 0.00880 0.058 0.666
Azueta Tesechoacan 57.42 8.35 0.145 0.00706 0.690 7.983
Bellaco Lalana 53.97 14.93 0.277 0.00751 1.196 13.845
Cuatotolapan San Juan 69.12 16.20 0.234 0.00586 2.627 30.403
Jacatepec Valle Nacional 35.20 9.54 0.271 0.01151 1.056 12.228
La Estrella Usila 25.05 5.21 0.208 0.01618 2.079 24.061
Monterrosa Cajones 23.49 6.04 0.257 0.01725 2.850 32.988
Qiotepec Grande 26.00 6.94 0.267 0.01559 15.619 180.778
Xiquila Xiquila 116.64 18.28 0.157 0.00347 0.329 3.811
1 __
    
x = media; s = desviación estándar; Cv = coeficiente de variación.

El Cuadro 3 muestra las rocas predominantes, las cuales varían desde rocas ígneas (basaltos) a sedimentarias (areniscas, calizas) a metamórficas (esquisto).

Cuadro 3.  Rocas predominantes en las subcuencas del río Papaloapan. 1
Estación de aforo
Corriente
o río
Tipos de roca predominante
Achotal La Trinidad Arenisca, esquisto, arenisca calcárea, caliza
Angel R. Cabadas Tecolapa Bssalto, tufa basáltica
Azueta Tesechoacan Arenisca, arenisca calcárea, caliza, conglomerado
Bellaco Lalana Arenisca, esquisto, arenisca calcárea, conglomerado
Cuatotolapan San Juan Arenisca, caliza, esquisto, arenisca calcárea, conglomerado
Jacatepec Valle Nacional Esquisto, arenisca calcárea, caliza
La Estrella Usila Esquisto, arenisca calcárea, caliza
Monterrosa Cajones Esquisto, andesita, arenisca calcárea, caliza, monzonita
Qiotepec Grande Granito metamorfoseado, esquisto, caliza,
lutita, arenisca calcárea, arenisca, conglomerado
Xiquila Xiquila Caliza, arenisca, conglomerado, andesita
1 Mapas Orizaba (E14-6), Oaxaca (E14-9), Coatzacoalcos (E15-4), y Minatitlán (E15-7).

El Cuadro 4 muestra los tipos de roca predominantes agrupados en base al tiempo de almacenamiento.

Cuadro 4.  Tipos de roca predominante agrupados en base al tiempo de almacenamiento.
Grupo Estación de aforo Tiempo de
almacenamiento
(días)
 Tipos de roca predominante
I Jacatepec, La Estrella,
Monterrosa, Quiotepec
23.5 - 35.2 Esquisto, granito metamorfoseado,
arenisca calcárea, caliza, arenisca, lutita
II Angel R. Cabadas 46.0 Basalto, tufa basáltica
III Achotal, Azueta,
Bellaco, Cuatotolapan
51.3 - 69.1 Arenisca, arenisca calcárea, caliza,
conglomerado, esquisto
IV Xiquila 116.6 Caliza, arenisca, conglomerado,
andesita

El Grupo I Grupo I (Jacatepec, La Estrella, Monterrosa, y Quiotepec) consiste de subcuencas con un tiempo de almacenamiento relativamente corto (23 a 35 días).

Estas subcuencas están localizadas lejos de la costa, en dirección sur a sureste (Fig. 1), a través de las montañas, con climas que varían de hiperhúmedo a árido, y consistiendo primeramente de rocas metamórficas, con la presencia de algunas rocas sedimentarias.

Estos acuíferos drenan en forma relativamente rápida.

El Grupo II (Angel R. Cabadas) consiste de una cuenca con tiempo de almacenamiento intermedio (46 días). Esta subcuenca está localizada en la costa, hacia el noreste, presentando un clima húmedo y rocas basálticas.

El Grupo III (Achotal, Azueta, Bellaco, y Cuatotolapan) consiste de subcuencas con un tiempo de almacenamiento relativamente largo (51 a 69 días).

Estas subcuencas están localizadas en las llanuras del este, con clima húmedo y mayormente rocas sedimentarias y algunas metamórficas.

El Grupo IV (Xiquila) consiste de una cuenca con un tiempo de almacenamiento muy largo (116 días).

Esta subcuenca está localizada lejos de la costa, en dirección suroeste, drenando la Sierra Madre, con un clima árido, y mayormente rocas sedimentarias y algunas rocas volcánicas.

Los acuíferos del Grupo I, donde predomina el esquisto, son incapaces de mantener flujos de base por períodos largos.

Contrariamente a esto, en el Grupo IV, donde predomina la caliza, el acuífero es capaz de mantener flujos de base por períodos relativamente largos.


5.  VERIFICACIÓN DE CAMPO

El Cuadro 5 muestra datos geográficos e hidrogeológicos en siete (7) pozos localizados en la cuenca del río Papaloapan.

Los valores de difusividad hidráulica, aunque restringidos a la parte norte-centro de la cuenca, se comparan favorablemente en forma preliminar con aquéllos mostrados en el Cuadro 2.

La ausencia de datos en otras partes de la cuenca no permite una comparación más exhaustiva.

 Cuadro 5.  Estaciones de bombeo, localización geográfica y características hidrogeológicas
en la cuenca del río Papaloapan.1
Estación
de bombeo
Número
del pozo
Latitud Longitud Transmisibilidad T
(m2 seg-1)
Coeficiente de
almacenamiento S
Difusividad
hidráulica D
(m2 seg-1)
San José Independencia 03 18° 23' 96° 03' 0.0550 0.00250 22.000
Cosamaloapan 05 18° 22' 95° 48' 0.0021 0.01000 0.210
Paso Carretas 09 18° 41' 96° 08' 0.0470 0.06800 0.690
Río Moreno 11 18° 38' 96° 14' 0.0260 0.00790 3.290
Cuyucuenda 13 18° 47' 96° 16' 0.0082 0.00043 19.070
Piedras Negras 14 18° 46' 96° 10' 0.0460 0.09400 0.489
Ignacio de la Llave 15 18° 43' 96° 59' 0.0043 0.09400 0.046
1 Fuente:  Comisión Nacional del Agua, Jalapa, México.


6.  CONCLUSIONES

El modelo de Rorabaugh se ha utilizado para estimar parámetros regionales de acuíferos en la cuenca del río Papaloapan, localizada en los estados de Veracruz y Oaxaca, en el sur de México.

Se han calculado el tiempo de almacenamiento, la constante de la cuenca, y la difusividad hidráulica para diez (10) subcuencas.

Los resultados muestran un acuerdo razonable con la geología local (extraída de mapas geológicos) y las pruebas de bombeo disponibles.


130521 05:00