papaloapan river basin map
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


Junio 2013


papaloapan river basin hydrograph

RESUMEN

El modelo de flujo de base de Rorabaugh (1963) se utiliza para estimar parámetros regionales de acuíferos en la cuenca del río Papaloapan, ubicada en el sur de México. Ésta es una cuenca tropical extensa, con un área de 46,527 km2, una gran diversidad de condiciones climáticas, y una importante interacción entre la hidrología superficial y la hidrología subterránea. Se han utilizado curvas de recesión de caudales para calcular el tiempo de almacenamiento, la constante de la cuenca, y la difusividad hidráulica del acuífero, en diez (10) puntos localizados a lo largo y ancho de la cuenca. Los resultados preliminares muestran un acuerdo razonable entre los parámetros calculados y los correspondientes parámetros hidrogeológicos obtenidos del bombeo y la cartografía.


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 (Hall 1968; Ponce 1989), la constante de la cuenca (Rorabaugh y Simons 1966), y la difusividad hidráulica (Freeze y Cherry 1979). 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 (Moore 1992).


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 (Freeze y Cherry 1979); 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 (Rorabaugh 1963; Rutledge y Daniel 1994; Mau y Winter 1997).

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. Trainer y Watkins (1974) han usado el modelo de Rorabaugh para estimar transmisibilidades en la cuenca del río Alto Potomac. Estudios más recientes han aplicado el modelo de Rorabaugh para estimar la recarga de agua subterránea en varios escenarios hidrogeólogicos (Rutledge y Daniel 1994; Mau y Winter 1997; Sanz 1997). 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 (Hall 1968). 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 (Hall 1968; Ponce 1989; Tallaksen 1995), 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 (Rorabaugh y Simons 1966):

                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 (Freeze y Cherry 1979):

            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. 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.

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 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) (Mapa Geológico de México, escala 1:250,000, publicado por INEGI).

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. Esta experiencia subraya la promesa de este método para estimar parámetros regionales de acuíferos utilizando los datos de recesión del hidrograma.


BIBLIOGRAFÍA

Boussinesq, J. 1877. Essai sur la theories des eaux courantes. Memoires presentes par divers savants a l'Academic des Sciences de l'Institut National de France, Tome XXIII,No. 1.

Comisión del Papaloapan. 1972. Boletín Hidrométrico No. 19: 1971-1972. Secretaría de Recursos Hidráulicos, Mexico City, Mexico.

Freeze, R. A., and Cherry. J. A. 1979. Groundwater. Prentice Hall, Englewood Clis, New Jersey.

Hall, F. R. 1968. Base-flow recessions - A review. Water Resources Research, 4(5), 973-983.

Mau, D. P., and Winter, T. C. 1997. Estimating groundwater recharge from streamflow hydrographs for a small mountain watershed in a temperate humid climate, New Hampshire, USA. Ground Water, 35(2), March-April, 291-304.

Moore, G. K. 1992. Hydrograph analysis in a fractured rock terrain. Ground Water, 30(3), May-June, 390-395.

Ponce, V. M. 1989. Engineering Hydrology: Principles and Practices. Prentice Hall, Englewood Cliffs, New Jersey.

Rorabaugh, M. I. 1963. Estimating changes in bank storage in groundwater contributions to streamflow. International Association of Scientic Hydrology, Publication No. 63, 432-441.

Rorabaugh, M. I., and Simons, W. D. 1966. Exploration of methods of relating groundwater to surface water, Columbia river basin - second phase. U.S. Geological Survey Open-file Report, March.

Rutledge, A. T., and Daniel, III, C. C. 1994 Testing an automated method to estimate groundwater recharge from streamflow records. Ground Water, 32(2), March-April, 180-189.

Sanz Pérez, E. 1997. Estimation of basinwide recharge rates using springflow, precipitation,and temperature data. Ground Water, 35(6), November-December, 1958-1965.

Tallaksen, L. M. 1995. A review of baseflow recession analysis. Journal of Hydrology, 165,349-370.

Trainer, F. W., and Watkins Jr., F. A. 1974. Use of base-runoff recession curves to determine real transmissivities in the Upper Potomac River Basin. U.S. Geological Survey Journal of Research, 2(1), 125-131.


1  Profesor de Ingeniería Civil y Ambiental, Universidad Estatal de San Diego, California.
2  Estudiante de postgrado en ingeniería civil, Universidad Estatal de San Diego, California.
3  Científico F, Instituto Nacional de Hidrología, Roorke, Uttaranchal, India.


130515 18:15