Modelación del Comportamiento Hidrodinámico del agua subterránea

de la zona comprendida entre Carbó, Pesqueira y Zamora, Sonora.

8.- BALANCE INTEGRAL DE AGUAS SUBTERRÁNEAS

Un balance de aguas subterráneas consiste en registrar las entradas, salidas y cambio en el

volumen de almacenamiento, que acontecen en un volumen específico del sistema acuífero en un
determinado lapso de tiempo; por tanto, debemos de conocer las componentes que intervienen en
el área de balance, como pueden ser la evapotranspiración, bombeo, caudal base, flujo subterráneo
que entra o sale lateralmente, etc. Para lograr este fin, se realizó el balance para la zona de valle
dentro del área de estudio.

La ecuación general de balance, de acuerdo a la ley de la conservación de la masa es:

Entradas (E) - Salidas (S) = Cambio de Almacenamiento

Aplicando esta ecuación, las entradas quedan representadas por la recarga total, las salidas
por la descarga total y el cambio de masa por el cambio de almacenamiento de una unidad
hidrogeológica.

Recarga total – Descarga total = Cambio de almacenamiento en la unidad hidrogeológica

De esta forma la ecuación de balance queda de la siguiente manera:

Rv + Eh + Ri – B – Sh – ETR = ±ΔV(S)

Donde:
Rv= Recarga vertical por lluvia
Eh= Recarga por flujo subterráneo
Ri= Recarga inducida
B= Bombeo
Sh= Salidas por flujo horizontal
ETR= Evapotranspiración real
AV(S)= Cambio de volumen almacenado

8.1 Parámetros Hidráulicos

El agua subterránea se desplaza a través de una serie de sedimentos y rocas de origen

diverso que conforman estructuras tectónicas de tipo horst y graben, estos últimos rellenos por
materiales sedimentarios de origen continental, los cuales se agrupan en dos medios físicos
diferenciables. Estos medios son: i) granular y ii) fracturado.

Tesis Para Obtener el Grado de Maestro en Ciencias - 40 -

 Modelación del Comportamiento Hidrodinámico del agua subterránea
de la zona comprendida entre Carbó, Pesqueira y Zamora, Sonora.

Las rocas que conforman a los horst y graben, agrupan alternadamente y en superposición
a los medios fracturados (andesitas, ignimbrítas y tobas) y granular (depósitos sedimentarios
continentales) con conductividades hidráulicas obtenidas a partir de la interpretación de pruebas de
bombeo mediante métodos diversos.

• Para las andesitas, ignimbrítas y tobas se desconoce su conductividad hidráulica debido a que
no existe información para su estimación y valoración.
• El medio granular se consideró anisotrópico, debido a la estratificación horizontal de
sedimentos; por tanto, los valores de conductividad hidráulica vertical se podrían estimar en un
orden de magnitud menor a la Kh, como lo asignan en la literatura mundial.

Transmisividad
Los valores de transmisividad para este sistema acuífero varían de 2x10-4 a 2x10-2 m2/s. La
principal aportación de agua subterránea proviene de la tercera capa que se conforman por gravas
y arenas.

Recarga Natural Inducida o Retorno de Riego

Del número total de aprovechamientos, 192 son de uso agrícola con un volumen de
3
extracción general de 13.03 Hm /año, de esto se asumió un 15% del gasto total es decir 1.95
Hm3/año para el retorno de riego o recarga inducida bajo referencia de Fetter, 1994.

Bombeo (Hm3/año)
En base al estudio (Consolidación del Comité Técnico Administrativo (COTAS) del río
Zanjón) elaborado por la empresa IDEAS (año 2007) se calculó una extracción total de 85 Hm3
anuales en los aprovechamientos activos actualmente, y siendo 89 Hm3 anuales el volumen
concesionado inscrito en el Registro Público de Derechos de Agua (REPDA), con una diferencia de
4Hm3/año de lo calculado con lo concesionado. Cabe mencionar que el balance se llevó a cabo
limitándose al área del acuífero administrativo del Río Zanjón que se encuentra dentro del área de
estudio.

Entradas y Salidas Horizontales (Hm3/año)

Se calculó mediante el método de celdas de flujo, basado en la ecuación de Darcy, donde
se tomaron como base las curvas de igual elevación al nivel estático y se obtuvieron los
diferenciales de carga y de distancia entre las curvas (gradiente dh/dl).

Tesis Para Obtener el Grado de Maestro en Ciencias - 41 -

 Modelación del Comportamiento Hidrodinámico del agua subterránea
de la zona comprendida entre Carbó, Pesqueira y Zamora, Sonora.

Con base en esto, se ubicaron las celdas en las los extremos del acuífero donde es visible
según las curvas que existe un flujo subterráneo de entrada o salida (Figura 8.1). La ley de Darcy,
se expresa de la siguiente manera (Fetter, 1994):

Q = - K A (h1-h2) / L
donde:

Q = Gasto, en m3/s
K = coeficiente de permeabilidad (m/s)
A = Área de la sección transversal por donde pasa el flujo, en m2
h1, h2 = Cargas hidráulicas (m)
L = Distancia entre curvas de igual valor (m), o en su caso distancia entre piezómetros

El signo negativo solo indica que la dirección del flujo es hacia donde el gradiente hidráulico
disminuye.

Conductividad Hidráulica (K)

Se obtiene a partir de las pruebas de bombeo, éstas arrojan un valor de transmisividad,
que viene a ser el producto de la permeabilidad por el espesor saturado:
T=Kb

Figura 8.1- Celdas de entrada de flujo en la zona de estudio.

Tesis Para Obtener el Grado de Maestro en Ciencias - 42 -

 Modelación del Comportamiento Hidrodinámico del agua subterránea
de la zona comprendida entre Carbó, Pesqueira y Zamora, Sonora.

8.2 Cambio de almacenamiento (Hm3/año)

Normalmente se consideran como entradas al acuífero: entradas por flujo horizontal de

otros acuíferos (EH), e infiltración vertical directa de lluvia y corrientes superficiales (IV). Pero como
es muy difícil evaluar las entradas al acuífero, se evalúan las salidas, que son: salidas por flujo
subterráneo horizontal (SH), descarga de manantiales (M) (este incluye el flujo base mencionado
anteriormente, siempre y cuando la descarga sea del acuífero que se está evaluando), extracciones
mediante pozos de bombeo (B), y evapotranspiración (ETR). El cambio de almacenamiento (Δ V)
será la suma de estos parámetros.

El cambio de almacenamiento en el intervalo de tiempo considerado en el balance, se

determina mediante la evolución de los niveles de agua en el intervalo de tiempo considerado y
valores representativos del coeficiente de almacenamiento del acuífero.

Hay varias formas para obtener el cambio del almacenamiento, una manera es obtenerla a
partir de pruebas de bombeo o de consideraciones relativas al tipo y litología de la zona de estudio,
otro modo (el que se utilizó en este trabajo) es crear una superficie y calcular el volumen y
multiplicarlo por un coeficiente de almacenamiento, esto se realiza para el inicio y fin del intervalo
de tiempo considerado a partir de datos de evolución piezométrica, como se muestra la Tabla 8.1

Tabla 8.1. Cálculo del cambio de almacenamiento para el período 2006 – 2007.

Evolución Abatimiento 2
∆V(S) (m /año)
3
Area (m ) S
(m) (m)
-5 -5 992623.28 0.08 -397049.312
-5 a -3 -4 4147939.51 0.08 -1327340.643
-3 a -1 -2 10175281.95 0.08 -1628045.112
-1 -1 28140760.00 0.08 -2251260.800
-5,603,695.9

Las superficies que se utilizaron para calcular el cambio de almacenamiento anteriormente

descrito, fueron sobre las zonas de abatimiento en la parte sur del área de estudio a la altura de los
poblados Estación Pesqueira, Zamora y El Zacatón como se aprecia la Figura 8.2, éstas son las
partes más afectadas por el sobrebombeo en el área de estudio.

Tesis Para Obtener el Grado de Maestro en Ciencias - 43 -

 Modelación del Comportamiento Hidrodinámico del agua subterránea
de la zona comprendida entre Carbó, Pesqueira y Zamora, Sonora.

Figura 8.2- Polígonos para el cálculo del cambio del almacenamiento.

8.3 Recarga Vertical (Rv)

La recarga vertical es uno de los términos que mayor incertidumbre implica a su cálculo.
Debido a que el cambio de almacenamiento (ΔV) fue considerado de acuerdo a la Tabla 8.1,
también se tiene información para calcular las entradas y salidas por flujo subterráneo con base en
la Ley de Darcy, el valor de la recarga vertical fue despejado de la ecuación de balance definida por
la siguiente expresión.

Rv + Eh + Ri – B – Sh – ETR = ± ∆V(S)

De esta manera, despejando la recarga vertical:

Rv = B + Sh + ETR +∆V(S) – Eh – Ri

Rv= 85 + 0 + 0 + (-5.6) - 53.90 - 8.08

Rv= 17.42 Hm3/año

Tesis Para Obtener el Grado de Maestro en Ciencias - 44 -

 Modelación del Comportamiento Hidrodinámico del agua subterránea
de la zona comprendida entre Carbó, Pesqueira y Zamora, Sonora.

La recarga total de aguas subterráneas será:

ENTRADAS m³/año
Entradas Horizontales (celdas de flujo) 53,906,104.06
Retorno de riego 8,085,915.61
Infiltración de cuerpos de agua 0.00
Recarga vertical 17,404,284.47

Recarga Total 79,396,304.13

SALIDAS
Bombeo 85,000,000.00
Salidas Horizontales (celdas de flujo) 0.00
Descarga de manantiales (flujo base) 0.00
Evapotranspiración en niveles someros 0.00

Salidas totales 85,000,000.00

CAMBIO DE ALMACENAMIENTO -5,603,696

La figura 8.3 muestra un modelo conceptual de la zona de balance para el área de estudio, con los
valores obtenidos en el balance de aguas subterráneas.

Figura 8.3 Modelo de balance conceptual en el área de estudio.

Tesis Para Obtener el Grado de Maestro en Ciencias - 45 -