Saltar al contenido

Curva de duración de caudales: construcción, Q95 y aplicaciones

Pablo RojasActualizado el 13 min de lectura

En este artículo
  1. Definición y usos
  2. Construcción de la curva
  3. Ecuaciones y variables
  4. Caudales característicos y sus usos
  5. Forma de la curva e índice de flujo base
  6. Transferencia a cuencas sin datos
  7. Relación con la curva de variación estacional
  8. Ejemplo resuelto: curva y central de pasada
  9. Errores frecuentes
  10. Preguntas frecuentes
  11. Fuentes
Río andino de caudal moderado corriendo entre rocas, con estación de aforo en la ribera

La curva de duración de caudales indica el porcentaje del tiempo en que cada caudal de un río fue igualado o superado durante un período de registro. Se construye ordenando caudales medios diarios o mensuales, y de ella se leen caudales característicos como el Q50 o el Q95, usados en derechos de agua, centrales de pasada, caudal ecológico y calidad del agua.

Definición y usos

La curva relaciona el caudal QQ con la duración pp, el porcentaje del tiempo en que QQ fue igualado o superado. Searcy (1959) la define como una curva de frecuencia acumulada que reúne todo el rango de caudales de un río, desde las crecidas hasta el estiaje; en términos estadísticos es el complemento de la función de distribución acumulada de los caudales (Vogel y Fennessey, 1994).

La curva no conserva el orden cronológico: en un registro de 15 años, los caudales menores que el de 96,7 % pueden concentrarse en un solo semestre de sequía (Searcy, 1959). Por eso la duración no es un período de retorno, que se obtiene con el análisis de frecuencia de valores anuales. La curva describe solo su período y anticipa caudales futuros cuando ese período representa el largo plazo.

Vogel y Fennessey (1994) mencionan usos en energía hidroeléctrica, calidad del agua, sedimentación, hábitat y aumento de caudales bajos.

Construcción de la curva

Se usan caudales medios de una estación fluviométrica, obtenidos con aforos y la curva de descarga, con al menos 10 años de registro (Risley et al., 2008). Si hubo embalses o extracciones, el régimen natural se describe solo con el tramo no regulado (Searcy, 1959). El procedimiento es:

  1. Reunir los nn caudales medios del período, sin lagunas.
  2. Ordenarlos de mayor a menor sin considerar año ni mes, con mm = 1 para el mayor; los valores repetidos reciben órdenes distintos.
  3. Calcular la duración con la posición de graficación de Weibull, p=100m/(n+1)p = 100\,m/(n+1).
  4. Graficar QQ contra pp con escala logarítmica de caudales, para leer crecidas y estiaje en el mismo gráfico (Searcy, 1959).
  5. Interpolar entre los dos valores ordenados que encierran la duración buscada.

Los caudales mensuales suavizan crecidas y estiajes: en un estudio de Carolina del Norte citado por Searcy (1959), las curvas mensual y diaria difirieron hasta en 35 %. Con datos mensuales, el Q90 es el caudal medio superado en el 90 % de los meses, aunque en parte de ellos hubo días con caudal menor.

Curva del período total y curva del año mediano

Pocos años extremos desplazan la curva del período total. Vogel y Fennessey (1994) propusieron construir una curva por año y tratar las nn curvas anuales como una serie de valores anuales. La mediana de cada duración forma la curva del año mediano, que representa un año típico; los demás percentiles dan intervalos de confianza y períodos de retorno (Castellarin et al., 2004). La curva del período total sirve para calcular energía. Con datos mensuales, cada curva anual llega solo a 100×12/13100 \times 12/13 = 92,3 %.

Ecuaciones y variables

Posición de graficación de Weibull (Risley et al., 2008), con pp la duración (%), mm el orden de mayor a menor y nn el número de caudales:

p=mn+1×100p = \frac{m}{n+1} \times 100

Para una duración dada se calcula m=p(n+1)/100m^* = p\,(n+1)/100 y se interpola entre los caudales ordenados Q(k)Q_{(k)} y Q(k+1)Q_{(k+1)}, con kk la parte entera de mm^*:

Qp=Q(k)+(mk)[Q(k+1)Q(k)]Q_p = Q_{(k)} + \left(m^* - k\right)\left[Q_{(k+1)} - Q_{(k)}\right]

El área bajo la curva dividida por 100 es el caudal medio del período (Searcy, 1959). En una central de pasada, el caudal turbinado QtQ_t descuenta el caudal ecológico QecoQ_{eco} y se limita al caudal de diseño QdQ_d, todos en m³/s:

Qt=min[max(QQeco, 0), Qd]P=ρgQtHηQ_t = \min\left[\max\left(Q - Q_{eco},\ 0\right),\ Q_d\right] \qquad P = \rho\,g\,Q_t\,H\,\eta

con PP la potencia (W), ρ\rho = 1.000 kg/m³, gg = 9,81 m/s², HH la caída neta (m) y η\eta la eficiencia global, cercana a 90 % en centrales actuales según el USBR. Con Pˉ\bar{P} la potencia media y PinstP_{inst} la potencia con Qt=QdQ_t = Q_d, ambas en kW, la energía anual EE (kWh) y el factor de planta FPFP son:

E=8760PˉFP=PˉPinstE = 8760\,\bar{P} \qquad FP = \frac{\bar{P}}{P_{inst}}

Caudales característicos y sus usos

CaudalDefiniciónUso y referencia
Q10Superado el 10 % del tiempoCaudal alto; frente al Q95 mide la variabilidad del régimen (NRFA)
Q50MedianaCaudal típico (Searcy, 1959)
Q85Superado el 85 % del tiempoMáximo para derechos permanentes en Chile, con caudales mensuales (DGA, 2024)
Q90Superado el 90 % del tiempoCaudal bajo; con el Q50 describe el tramo de estiaje
Q95Superado el 95 % del tiempoCaudal bajo de referencia para calidad del agua (NRFA)
Q347Superado 347 días al año, promedio de diez años, sin alteracionesCaudal residual mínimo en Suiza

Disponibilidad hídrica y derechos de agua

El manual de la DGA (2024) fija como máximo 85 % de probabilidad de excedencia para el caudal de un derecho permanente, y 10 %, o excepcionalmente 5 % con al menos 30 años de estadística, para uno eventual. La oferta se evalúa mes a mes: el caudal disponible es el de régimen natural menos el caudal ecológico mínimo y los derechos permanentes de aguas arriba, dentro del balance hídrico de la cuenca.

Generación hidroeléctrica de pasada

Una central de pasada genera con el caudal que trae el río, sin almacenar agua; hacia 1920 el USGS adoptó la curva de duración para sus estadísticas de energía hidráulica (Searcy, 1959). Un QdQ_d cercano al Q85 da plena carga cerca del 85 % del tiempo si QecoQ_{eco} es pequeño, pero desaprovecha volumen; un QdQ_d alto genera más con menor factor de planta.

Caudal ecológico

La ley suiza de protección de aguas define el Q347 y, en su artículo 31, fija un caudal residual mínimo que crece con él, desde 50 l/s para un Q347 de hasta 60 l/s hasta 10.000 l/s desde 60.000 l/s. Como 347 días son el 95,1 % del año, el Q347 se aproxima al Q95 de caudales diarios. En Chile, el caudal ecológico mínimo se calcula con caudales medios mensuales y no con el Q95.

Dilución y calidad del agua

El National River Flow Archive (NRFA) del Reino Unido destaca el Q95 para fijar condiciones de calidad en autorizaciones de descarga. Con mezcla completa, la concentración aguas abajo es:

C=QCr+QeCeQ+QeC = \frac{Q\,C_r + Q_e\,C_e}{Q + Q_e}

con QQ y QeQ_e los caudales del río y del efluente (m³/s) y CrC_r y CeC_e sus concentraciones (mg/L). Con carga constante, la concentración calculada con el Q95 se supera cerca del 5 % del tiempo.

Forma de la curva e índice de flujo base

Según Searcy (1959), una pendiente fuerte indica un río variable alimentado por escorrentía directa, y una suave, almacenamiento superficial o subterráneo. Un tramo inferior plano revela gran almacenamiento permanente, ligado a la geología; un tramo superior plano es típico de ríos con caudales altos de deshielo.

La razón Q10/Q95 mide la variabilidad del régimen (NRFA), y Q90/Q50, la pendiente del tramo de estiaje: cuanto más cerca de 1, más sostenido es el caudal base. El índice de flujo base (BFI) es la fracción de la escorrentía que proviene de almacenamientos; el NRFA lo calcula con caudales diarios, con valores de 0,15 a 0,35 en ríos sobre arcillas impermeables y sobre 0,9 en la mayoría de los ríos sobre creta. A diferencia de las razones, requiere la secuencia del hidrograma.

Transferencia a cuencas sin datos

Para un sitio sin estación, el método de la razón de áreas supone igual caudal por unidad de área que en una estación cercana:

Qu=AuAgQgQ_u = \frac{A_u}{A_g}\,Q_g

con QuQ_u y QgQ_g los caudales de igual duración en el sitio y en la estación (m³/s) y AuA_u y AgA_g sus áreas de drenaje (km²). Risley et al. (2008) lo consideran confiable solo si ambos sitios están cerca. La DGA (2024) agrega la precipitación Pr\mathrm{Pr} de cada cuenca, K=A1Pr1/(A2Pr2)K = A_1\,\mathrm{Pr}_1/(A_2\,\mathrm{Pr}_2), en cuencas similares. En estiaje, Searcy (1959) advierte que el caudal por unidad de área de cuencas vecinas rara vez sirve si la geología de aguas subterráneas difiere; unos pocos aforos de flujo base en el sitio dan resultados mucho mejores.

Relación con la curva de variación estacional

La curva de variación estacional hace un análisis de frecuencia por mes y grafica, a lo largo del año hidrológico, los caudales medios mensuales de varias probabilidades de excedencia; la DGA (2024) pide las de 10, 50, 85 y 95 %. Su Q85 de octubre se supera en el 85 % de los octubres, mientras que el Q85 de la curva de duración se supera en el 85 % de todos los meses. Los derechos de agua se evalúan con la primera.

Ejemplo resuelto: curva y central de pasada

Datos ilustrativos: la serie no corresponde a ninguna estación real; se generó con python3 para mostrar el procedimiento. Son 10 años de caudales medios mensuales (m³/s), de abril a marzo.

AñoAbrMayJunJulAgoSepOctNovDicEneFebMarMedia
12,142,944,103,822,773,803,694,856,204,571,652,323,57
20,823,134,013,464,302,893,914,596,314,451,701,263,40
32,333,353,5910,123,353,074,288,5612,435,643,222,525,21
41,402,523,255,024,142,003,895,516,663,491,841,393,43
50,991,192,053,401,851,622,263,224,281,991,010,962,07
62,032,064,355,654,634,034,375,175,566,242,242,594,08
71,506,163,758,544,104,395,176,0612,655,622,432,595,25
81,591,724,214,144,073,583,564,356,444,011,591,363,39
91,442,457,274,424,445,184,696,168,254,602,712,094,48
101,231,572,383,283,812,403,794,294,051,731,911,192,64

Paso 1. Curva y caudales característicos. Los nn = 120 caudales se ordenan de mayor a menor, con p=100m/121p = 100\,m/121. Para el Q85, m=0,85×121=102,85m^* = 0{,}85 \times 121 = 102{,}85, y los caudales de orden 102 y 103 son 1,70 y 1,65 m³/s:

Q85=1,70+0,85×(1,651,70)=1,66 m3/sQ_{85} = 1{,}70 + 0{,}85 \times \left(1{,}65 - 1{,}70\right) = 1{,}66\ \text{m}^3\text{/s}
Caudalmm^*QpQ_p (m³/s)
Q1012,106,20
Q5060,503,59
Q85102,851,66
Q90108,901,45
Q95114,951,19

El caudal medio es 3,75 m³/s, y las razones de forma son Q10/Q95 = 5,2 y Q90/Q50 = 0,40.

Paso 2. Curvas anuales. Con p=100m/13p = 100\,m/13 en cada año, la mediana de las diez curvas da Q50 = 3,72 m³/s y Q85 = 1,82 m³/s, sobre el período total, que concentra los meses secos de los años 5 y 10. En el año 5, el Q85 bajó a 0,99 m³/s.

Paso 3. Energía de una central de pasada. Supuestos: HH = 45 m; η\eta = 0,85, bajo el 90 % del USBR por prudencia; QecoQ_{eco} = 0,40 m³/s; QdQ_d = 3,5 m³/s. El promedio de QtQ_t en los 120 meses es 2,646 m³/s:

Pˉ=1000×9,81×2,646×45×0,85=992,8×103 WE=8760×992,8=8,70×106 kWh\bar{P} = 1000 \times 9{,}81 \times 2{,}646 \times 45 \times 0{,}85 = 992{,}8 \times 10^{3}\ \text{W} \qquad E = 8760 \times 992{,}8 = 8{,}70 \times 10^{6}\ \text{kWh}

La central genera en promedio 8,70 GWh al año, con factor de planta de 75,6 %, y opera a plena carga en el 43,3 % de los meses, aquellos con al menos Qd+QecoQ_d + Q_{eco} = 3,9 m³/s. Otros caudales de diseño:

QdQ_d (m³/s)PinstP_{inst} (kW)EE (MWh/año)Factor de planta (%)Meses a plena carga (%)
1,55634.62793,880,8
2,59386.92784,361,7
3,51.3138.69775,643,3
4,51.6899.66765,420,0
6,02.25110.37052,67,5

Al pasar QdQ_d de 1,5 a 6,0 m³/s, la potencia instalada se multiplica por 4 y la energía solo por 2,2, porque los caudales altos duran poco. Con serie mensual la energía queda sobrestimada: en un mes con caudal medio igual a QdQ_d hay días cuyo exceso se vierte y días bajo plena carga.

Errores frecuentes

  • Leer la duración como período de retorno. El Q95 se supera el 95 % del tiempo del registro, y los días bajo él pueden concentrarse en una sola sequía.
  • No declarar el intervalo de tiempo. El Q95 mensual suele ser mayor que el diario del mismo río.
  • Tratar un registro alterado como régimen natural, pese a embalses o extracciones.
  • Usar menos de 10 años o un período sin años secos, que deja mal definido el estiaje.
  • Confundir el Q85 de la curva con el Q85 mensual que usa la DGA para derechos permanentes.
  • Transferir por área sin revisar precipitación ni geología, que controlan el tramo de estiaje.
  • Calcular la energía con el caudal medio. En el ejemplo, sobrestima la energía en 27 %.

Preguntas frecuentes

¿Qué diferencia hay entre la curva de duración y el análisis de frecuencia?

La curva de duración usa todos los caudales y da la fracción del tiempo en que se supera cada valor; el análisis de frecuencia usa un valor por año y da el período de retorno.

¿Cuántos años de registro se necesitan?

Risley et al. (2008) usan como mínimo 10 años, y la DGA (2024) pide series mensuales de 30 años en lo posible. También importa que el período incluya años secos y húmedos con un régimen estable.

¿El Q347 es lo mismo que el Q95?

Casi: 347 de 365 días son el 95,1 %, pero la ley suiza promedia el valor de diez años y exige caudales no alterados por embalses, extracciones o aportes.

¿Qué caudal de diseño conviene para una central de pasada?

No hay un valor único. Se compara la energía adicional de cada aumento de QdQ_d con su costo, después de descontar el caudal ecológico.

Fuentes

  • Searcy, J. K. (1959). Flow-duration curves. Manual of hydrology: Part 2. Low-flow techniques (Water-Supply Paper 1542-A). U.S. Geological Survey. pubs.usgs.gov
  • Vogel, R. M. y Fennessey, N. M. (1994). Flow-duration curves. I: New interpretation and confidence intervals. Journal of Water Resources Planning and Management, 120(4), 485–504. doi.org
  • Castellarin, A., Vogel, R. M. y Brath, A. (2004). A stochastic index flow model of flow duration curves. Water Resources Research, 40(3). doi.org
  • Risley, J., Stonewall, A. y Haluska, T. (2008). Estimating flow-duration and low-flow frequency statistics for unregulated streams in Oregon (Scientific Investigations Report 2008-5126). U.S. Geological Survey. pubs.usgs.gov
  • Dirección General de Aguas (2024). Manual de normas y procedimientos para la gestión y administración de recursos hídricos (SDT N.º 477). Ministerio de Obras Públicas de Chile. dga.mop.gob.cl
  • Asamblea Federal de la Confederación Suiza (1991). Federal Act on the Protection of Waters (Waters Protection Act, WPA), SR 814.20, estado al 1 de enero de 2022. fedlex.admin.ch
  • UK Centre for Ecology & Hydrology (s. f.). Derived flow statistics. National River Flow Archive. nrfa.ceh.ac.uk
  • U.S. Bureau of Reclamation (s. f.). Hydroelectric power. usbr.gov