2º Congreso Nacional AMICA 2015
METODOLOGÍA PARA EL CÁLCULO DE CUENCA PROPIA (MÉTODO
INDIRECTO) UTILIZANDO LLUVIA DIARIA
Suárez Medina María de los Ángeles, Aguilar Garduño Ernesto, Hernández Solís Carlos Alberto,
Astudillo Enríquez Citlalli, Martínez Ocampo Delker Emmanuel.
Instituto Mexicano de Tecnología del Agua, Comisión Estatal del Agua de Jalisco, Consorcio Empresarial APDER
msuarez@[Link], eaguilar@[Link], cahernandez@[Link], xitlae@[Link],
[Link]@[Link]
Resumen continua debido a la falta de personal adecuado para la toma
de las mediciones, falta de mantenimiento; vandalismo en
Se entiende por Cuenca Propia (Cp) al volumen medio anual aquellas estaciones que están ubicadas en zonas de alto riesgo,
de escurrimiento natural. Esta variable es fundamental en el por mencionar algunas.
cálculo de balance hídrico que a su vez se fundamenta en el
El balance hídrico es el equilibrio entre los volúmenes de agua
principio de continuidad que dice que “las entradas menos las que entran y salen dentro de un sistema hídrico, si es que
salidas son iguales a cero”. puede almacenarse el recurso dentro del sistema (variación de
volumen V), en un intervalo de tiempo determinado. Las
El objetivo de este trabajo es presentar los pasos que se deben variables de entrada son: el volumen medio anual de
tomar en cuenta para el cálculo de la variable Cp en una escurrimiento natural conocido por Cp, volumen medio anual
cuenca hidrológica. de escurrimiento desde la cuenca aguas arriba (Ar), volumen
anual de retornos (R), volumen anual de importaciones (Im).
La variable Cp es el escurrimiento virgen donde el gasto Las variables de salida son: Evaporación en embalses (Ev),
volumen medio anual de escurrimiento de la cuenca hacia
generado en una cuenca no existen aprovechamientos.
aguas abajo (Ab), volumen anual de extracción de agua
Representa la máxima cantidad de agua superficial que una superficial (U) y volumen anual de exportaciones (Ex), que se
cuenca puede producir. Su valor medio anual es considerado, representa mediante la ecuación 1.
dentro de la planeación hidráulica, como parámetro
fundamental del potencial de los recursos hidráulicos ΔV = (Cp + Ar + R + Im) − (Ev + Ab + U + Ex) (1)
superficiales y se puede calcular de dos formas, una utilizando
los aforos de las estaciones hidrométricas o los registros de las En este artículo sólo se estudia la variable Cp para los casos en
presas ubicadas a la salida de una cuenca (método directo) y donde no se cuenta con suficiente información medida para
los escurrimientos de la cuenca hacia aguas abajo (Ab) y se
en caso de no contar con mediciones, la otra forma es
recurre a un método indirecto.
utilizando los registros de lluvia preferentemente diaria, para
aplicar algún método indirecto, para este estudio se tomó el Existen varios métodos indirectos para calcular el volumen de
método del “coeficiente de escurrimiento”. Este método escurrimiento por cuenca propia o escurrimiento virgen, y en
involucra los registros de lluvia de por lo menos 20 años de términos prácticos todos conceptualmente involucran las
acuerdo con la NOM-011-CNA-2015. mismas variables basados en datos de precipitación y
temperatura. Ciertos métodos usan información adicional
Desafortunadamente, en muchas de las estaciones como tipo de suelo, cobertura vegetal, características del cauce
principal como longitud y pendiente, además de las
climatológicas no se tiene la continuidad de los datos de lluvia
características fisiográficas de la cuenca.
por lo que ha sido necesario considerar los registros diarios
para la generación de los mapas de lluvia de la cuenca en Algunos de los métodos son: Temez, Coutagne, Turc,
estudio. El resultado de este estudio es la definición de los Langbein, Smith y Coeficiente de Escurrimiento, que
pasos a seguir para la obtención de Cp utilizando mapas de contemplan expresiones algebraicas derivadas de analizar el
lluvia diaria. Para definir esta metodología se consultaron comportamiento de una diversidad de cuencas.
(Suarez Medina et al, 2015) y la publicación de DOF de la
Metodología
publicación de la disponibilidad de agua superficial de la
RH12 Lerma Santiago Pacífico. La metodología para el cálculo de Cp es la siguiente:
Introducción 1) Cálculo de la precipitación
En México existen estaciones climatológicas con registros La estimación del volumen total de lluvia que ingresa a la
desde 1930, que miden, entre otras variables, la precipitación, cuenca en estudio se determinó con los rasters de lluvia diaria.
la temperatura, la evaporación y la dirección del viento; sin
embargo, la información registrada no se tiene de forma
2º Congreso Nacional AMICA 2015
Cálculo de la precipitación media diaria
Raster de precipitación diaria. Para obtener el raster de
Se seleccionaron las estaciones climatológicas (de la base de precipitación diaria se utilizó la herramienta de
datos CLICOM) que se ubicaban dentro de la zona de estudio Interpolación IDW (dentro de la caja de herramientas
y 30 km alrededor de ella. Interpolation en Spatial Analyst Tools), se ingresó el
dato de precipitación del shape de estaciones
seleccionadas y se obtuvo el raster.
Con el método de la interpolación IDW (Inverse
Distance Weighted) se asegura que el valor de salida
para una celda se limite al rango de valores utilizados
para la interpolación, es decir, no sea mayor que la
entrada máxima o inferior que la entrada mínima y como
resultado se obtiene el promedio de la distancia
ponderada. Para establecer una función de
proporcionalidad entre el peso y la distancia, la fórmula
general queda como sigue:
Figura 1. Estaciones climatológicas dentro de la zona de 𝑧̂𝑗 = ∑𝑛𝑖=𝑖 𝑘𝑖𝑗 ∙ 𝑧𝑖 (2)
estudio
En un proyecto de ArcGIS 10.1 se colocaron los mapas 𝛽
correspondientes a las estaciones climatológicas seleccionadas Donde 1⁄𝑑𝑖𝑗 y es un exponente de ponderación que
y se creó un módulo (ModelBuilder) para la generación de los controla la forma en la que el peso disminuye con la
raster, en el que se enlazaron las siguientes herramientas: distancia (d).
Select, IDW, CLIP y Zonal Statictics, Table to dBase.
Con la herramienta Clip del ArcToolbox se corta el
raster al contorno de la región.
Figura 2. Proyecto base en ArcGIS
La generación de rásters de precipitación se hizo de la
siguiente manera:
Figura 4. Ejemplo de raster de precipitación diaria
Selección de estaciones con registros para una fecha
Precipitación media por cuenca. Con la herramienta
específica. (Ejemplo el 3 de julio de 1980), a través de
Zonal que utiliza el raster de precipitación y la
una sentencia SQL y la Aplicación ArcMap del ArcGIS
configuración de la cuenca, se calculó, para cada
10.1, se seleccionó el shape de todas las estaciones con
una de las cuencas, el promedio de precipitación
valores de precipitación diferentes a <NULL>, dando
utilizando cada uno de los pixeles que caen dentro
como resultado un shape de estaciones.
del límite, almacenando los resultados en una tabla.
Figura 3. Selección de estaciones con datos para una fecha Figura 5. Exportación a tabla de Excel
2º Congreso Nacional AMICA 2015
Comparación de datos con el método de polígonos
de Thiessen. Los resultados obtenidos con el método
del IDW se compararon con los resultados obtenidos
con el método de polígonos de Thiessen, los cuales
se generaron con la herramienta Create Thiessen
Polygons (en la caja de herramientas Proximity de
Análysis Tools) utilizando las estaciones con
información del 3 de julio de 2015, por continuar
con el ejemplo antes mencionado.
Figura 8. Precipitación mensual y anual por cuenca
2) Determinación del área de la cuenca en estudio
Con la información de los modelos digitales del INEGI y la
herramienta del ArcGis, se calculó el área de la cuenca en
estudio.
3) Cálculo del coeficiente de escurrimiento
Figura 6. Polígonos de Thiessen El método del coeficiente de escurrimiento está en función del
coeficiente K que a su vez depende del tipo y uso del suelo de
Con los dos métodos utilizados se observa lo siguiente: la zona de estudio, para ello tomó en cuenta:
La precipitación media diaria obtenida con los si K resulta menor o igual que 0.15 entonces
rasters en general está por abajo con respecto a los
valores de precipitación obtenida por los polígonos Ce = K(P − 250)⁄2000 (3)
de Thiessen.
Pero si K es mayor que 0.15 entonces
El método empleado para obtener los rasters de
precipitación diaria hasta la cuenca no. 20 es muy Ce = K(P − 250)⁄2000 + (k − 0.15)⁄1.5 (4)
parecida por tener valores bajos.
En donde el coeficiente k se obtiene de las cartas de INEGI de
De la cuenca 21 en adelante las variaciones son más edafología y usos del suelo.
evidentes, siendo más altos los resultados de los
Polígonos de Thiessen. 4) Cálculo del volumen de Cp
Finalmente, se toman en cuenta los valores obtenidos de los
pasos anteriores para calcular el volumen de Cp con la
siguiente ecuación:
Cp = P ∙ A ∙ Ce (5)
dónde:
Cp es igual al volumen medio anual de
escurrimiento natural
P es la precipitación media anual
A es el área de la cuenca
Figura 7. Gráfico de comparación, día 03/07/1980 Ce es el coeficiente de escurrimiento
Integración de datos diarios. Una vez obtenidas las
tablas con la precipitación media diaria, se realizó la
integración de la información e un archivo de Excel
donde se agrupó por cuenca.
Esta integración permitió tener los datos a nivel mensual y
anual y a través de tablas dinámicas generar la consulta de los
mismos.
2º Congreso Nacional AMICA 2015
Conclusiones y recomendaciones Referencias
Esta metodología resulta muy valiosa para el cálculo de Suarez Medina, M. d., Patiño Gómez, C., Velázquez Álvarez,
Cuenca propia cuando las estaciones climatológicas ubicadas J., Rivera Benítez, J., Aguilar Garduño, E., Bautista, G.,
en una zona de interés no cuentan con suficiente información Astudillo Enríquez C. (2015). Aplicación del modelo de datos
para determinar la precipitación media anual. ArcHydro en el cálculo de disponibilidad de agua superficial.
Tecnología y Ciencias del Agua, Volumen VI Núm. 1.
Los mapas de precipitación generados por rásters, muestra la
tendencia de la lluvia según la influencia de la precipitación
DOF 2010. Estudio de disponibilidad de agua superficial de la
medida por otras estaciones climatológicas cercanas a la zona
RH12 Lerma Santiago Pacífico.
de estudio, dando como resultado una mejor configuración de
la lluvia.
Por otro lado, los polígonos de Thiessen, asignan una misma
lámina de lluvia para todo el polígono que tiene influencia en
el área de estudio.
Se recomienda hacer uso de los Sistemas de Información
Geográfica para simplificar la obtención de la lluvia
independientemente del método que se seleccione.
