Ingeniería Eléctrica · Operación de Sistemas de Potencia

Despacho Económico de Carga: Formulación Matemática, Algoritmos y Optimización

Cómo determinan los operadores del sistema la asignación óptima de potencia en cada generador para minimizar el costo total de combustible garantizando la seguridad de la red.

Análisis y Formulación Técnica · Operación y Control de Sistemas Eléctricos de Potencia (SEP)
Gen G1 C₁(P₁) Gen G2 C₂(P₂) OPERADOR (CENACE/ISO) Min ∑ Cᵢ(Pᵢ) dC₁/dP₁ = dC₂/dP₂ = λ ∑ Pᵢ = P_Demanda + P_L DEMANDA P_D (MW)
Publicidad

En cualquier instante, la energía eléctrica consumida por millones de usuarios en una red nacional debe ser generada de manera idéntica y en tiempo real. Sin embargo, no todas las plantas de generación cuestan lo mismo: quemar gas natural, diésel o turbinar agua tiene costos operativos drásticamente diferentes. El Despacho Económico de Carga (Economic Dispatch) es el problema de optimización matemática central que resuelve qué potencia exacta debe entregar cada central para satisfacer la demanda con el mínimo costo operativo posible, sin vulnerar los límites físicos de la red.

1. Formulación Matemática Fundamental

Consideremos un sistema con $N$ unidades generadoras térmicas conectadas a una red para alimentar una demanda total $P_D$. El objetivo primordial es minimizar la suma del costo total de combustible de todas las plantas operativas en un horizonte horario dado:

Función Objetivo:
Minimizar C_total = ∑ [ Cᵢ(Pᵢ) ] para i = 1, 2, ..., N

Para cada unidad térmica individual $i$, la función de costo de combustible suele aproximarse con alta fidelidad mediante una ecuación cuadrática en función de la potencia activa generada $P_i$ (en MW):

Cᵢ(Pᵢ) = αᵢ + βᵢ · Pᵢ + γᵢ · Pᵢ² [$ / hora]

Donde los coeficientes representan:

Costo Incremental de Combustible (IC)

El costo incremental (Incremental Cost) se obtiene calculando la primera derivada de la función de costo con respecto a la potencia activa generada:

ICᵢ = dCᵢ / dPᵢ = βᵢ + 2 · γᵢ · Pᵢ [$ / MWh]

El costo incremental indica cuántos dólares adicionales costará producir un megavatio-hora (MWh) más en esa central específica.

Publicidad

2. Despacho sin Pérdidas de Transmisión: Criterio de Costos Incrementales Iguales

En una primera aproximación ideal donde las distancias son cortas o las pérdidas en las líneas de transmisión se desprecian ($P_L \approx 0$), la única restricción de igualdad es el balance exacto entre generación y demanda:

Restricción de Igualdad:
∑ Pᵢ - P_D = 0

Y cada unidad tiene límites técnicos de capacidad mínima y máxima por estabilidad de combustión y refrigeración del estator:

Restricciones de Desigualdad:
Pᵢ_min ≤ Pᵢ ≤ Pᵢ_max

Deducción con Multiplicadores de Lagrange

Para encontrar el punto óptimo sujeto a una restricción de igualdad, se formula la función de Lagrange $\mathcal{L}$ introduciendo el multiplicador $\lambda$:

ℒ(P₁, P₂, ..., P_N, λ) = ∑ Cᵢ(Pᵢ) + λ · [ P_D - ∑ Pᵢ ]

Las condiciones necesarias de primer orden (Condiciones de Karush-Kuhn-Tucker / Euler-Lagrange) establecen que las derivadas parciales respecto a cada variable independiente deben anularse:

∂ℒ / ∂Pᵢ = (dCᵢ / dPᵢ) - λ = 0 ===> dCᵢ / dPᵢ = λ

El Teorema Fundamental del Despacho Económico establece que la condición para mínima factura de combustible ocurre cuando todas las unidades generadoras no saturadas operan exactamente al mismo costo incremental de combustible ($\lambda$).

Si una unidad alcanza su límite técnico superior o inferior, se fija en dicho límite ($P_i = P_{i,\max}$ si $dC_i/dP_i < \lambda$, o $P_i = P_{i,\min}$ si $dC_i/dP_i > \lambda$) y la potencia restante se redistribuye entre las unidades libres hasta igualar sus costos incrementales.

3. Despacho considerando Pérdidas de Transmisión (Factores de Penalización)

En redes eléctricas extensas de alta tensión (230 kV, 500 kV), transportar energía desde una central distante hasta los centros de consumo genera pérdidas por efecto Joule ($I^2 R$) que suelen oscilar entre el 2% y el 8% de la demanda total. Por ende, la restricción de balance real es:

∑ Pᵢ - P_D - P_L = 0

Donde $P_L = f(P_1, P_2, ..., P_N)$ es una función no lineal de las inyecciones de los generadores, descrita analíticamente mediante los coeficientes de pérdidas de Kron (Matriz B):

P_L = ∑∑ Pᵢ · Bᵢⱼ · Pⱼ + ∑ B₀ᵢ · Pᵢ + B₀₀

Al formular nuevamente el Lagrangiano considerando $P_L$:

ℒ = ∑ Cᵢ(Pᵢ) + λ · [ P_D + P_L - ∑ Pᵢ ]

∂ℒ / ∂Pᵢ = (dCᵢ / dPᵢ) + λ · (∂P_L / ∂Pᵢ) - λ = 0

(dCᵢ / dPᵢ) · [ 1 / (1 - ∂P_L / ∂Pᵢ) ] = λ

El Factor de Penalización ($L_i$)

Se define el Factor de Penalización de la barra $i$ como:

Lᵢ = 1 / ( 1 - ∂P_L / ∂Pᵢ )

Condición de Optimalidad: Lᵢ · (dCᵢ / dPᵢ) = λ para toda unidad i

Interpretación Física del Factor de Penalización ($L_i$)

El término $\frac{\partial P_L}{\partial P_i}$ es la pérdida incremental de transmisión. Si incrementar la potencia en la central $i$ aumenta fuertemente las pérdidas de la red hacia el centro de carga, $\frac{\partial P_L}{\partial P_i} > 0$, lo que hace que $L_i > 1$. Como resultado, su costo incremental aparente se "penaliza" (se encarece artificialmente), obligando al algoritmo a despachar menos de esa central lejana y más de una central cercana a la carga, incluso si su costo puro en bornes es ligeramente superior.

4. Algoritmo de Resolución: Iteración Lambda ($\lambda$-iteration)

Dado que la condición óptima relaciona todas las potencias a través del multiplicador común $\lambda$, el método numérico estándar más eficiente es la Iteración Lambda:

1
Establecer límites $\lambda_{min}$ y $\lambda_{max}$
2
Calcular $\lambda = \frac{\lambda_{min} + \lambda_{max}}{2}$
3
Calcular $P_i = \frac{\lambda - \beta_i}{2\gamma_i}$
4
Acotar $P_i$ a $[P_{i,min}, P_{i,max}]$
5
Evaluar error de balance: $\Delta P = \sum P_i - P_D - P_L$
6
Ajustar $\lambda$ hasta que $|\Delta P| < \varepsilon$
Potencia P (MW) λ ($/MWh) IC₁ = β₁ + 2γ₁P₁ IC₂ = β₂ + 2γ₂P₂ λ_óptimo P₁* P₂*
Convergencia del punto de despacho: para un mismo valor de λ, se proyectan las potencias óptimas P₁* y P₂* cuya suma cumple con la demanda.

5. Ejemplo Numérico Práctico Resuelto

Analicemos un caso de estudio real con dos plantas térmicas que deben alimentar una demanda consolidada de $P_D = 800\text{ MW}$ despreciando pérdidas:

Unidad Función de Costo Cᵢ(Pᵢ) [$/h] Costo Incremental ICᵢ [$/MWh] Límites Técnicos
Generador 1 (G1) C₁ = 500 + 15·P₁ + 0.005·P₁² IC₁ = 15 + 0.010·P₁ 100 MW ≤ P₁ ≤ 600 MW
Generador 2 (G2) C₂ = 400 + 18·P₂ + 0.008·P₂² IC₂ = 18 + 0.016·P₂ 100 MW ≤ P₂ ≤ 500 MW

Paso 1: Aplicar la condición de igualdad de costos incrementales:

IC₁ = IC₂ = λ
15 + 0.010·P₁ = λ ===> P₁ = (λ - 15) / 0.010
18 + 0.016·P₂ = λ ===> P₂ = (λ - 18) / 0.016

Paso 2: Sustituir en la ecuación de balance de demanda ($P_1 + P_2 = 800$):

(λ - 15) / 0.010 + (λ - 18) / 0.016 = 800
100·(λ - 15) + 62.5·(λ - 18) = 800
100λ - 1500 + 62.5λ - 1125 = 800
162.5·λ - 2625 = 800
162.5·λ = 3425 ===> λ = 21.0769 $ / MWh

Paso 3: Calcular la potencia asignada a cada generador:

P₁* = (21.0769 - 15) / 0.010 = 607.69 MW (¡Supera P₁_max = 600 MW!)

Al violarse el límite de G1, se fija en su cota máxima: $P_1 = 600\text{ MW}$. Por tanto, la unidad G2 debe asumir la diferencia:

P₂* = 800 - 600 = 200 MW
(Cumple: 100 ≤ 200 ≤ 500 MW)

El costo incremental del sistema ahora queda determinado por la unidad en margen (G2):

λ_sistema = IC₂(200) = 18 + 0.016·(200) = 21.20 $ / MWh

6. Despacho Moderno: Hidrotérmico, Renovables y Restricciones de Red (SCED)

En las redes eléctricas contemporáneas, el despacho económico ya no se ejecuta de forma aislada, sino dentro de un conjunto integrado de herramientas computacionales de alta velocidad:

A. Despacho Hidrotérmico y el Costo de Oportunidad del Agua

Las centrales hidroeléctricas no consumen combustible fósil, por lo que su costo marginal directo es prácticamente cero. Sin embargo, el agua embalsada es un recurso finito: usar 1 m³ de agua hoy impide usarlo durante el estiaje del próximo mes. Por ello, los operadores calculan el Costo Futuro del Agua (mediante programación dinámica estocástica como el modelo SDDP), asignándole un costo de oportunidad sombra que compite en el orden de mérito frente a las plantas térmicas.

B. Despacho Económico con Restricciones de Seguridad (SCED)

El Security-Constrained Economic Dispatch (SCED) incorpora las leyes de Kirchhoff mediante un Flujo de Carga Óptimo (OPF), garantizando que las líneas de transmisión no superen sus límites térmicos de transferencia y que el sistema soporte contingencias N-1 (la desconexión imprevista de la línea o generador más grande sin provocar colapso de frecuencia o voltaje).

C. Despacho con Energías Renovables Variables (Solar y Eólica)

Debido a la alta penetración de energía fotovoltaica y eólica de costo marginal cero, los operadores ejecutan despachos económicos en intervalos ultra rápidos (cada 5 a 15 minutos en mercados de tiempo real). Esto mitiga el fenómeno de la Curva de Pato, donde las plantas térmicas flexibles y los sistemas de almacenamiento BESS deben elevar su potencia bruscamente al atardecer.

7. Casos de Operación en el Mundo Real

Operadores del Sistema y Despacho en Latinoamérica

Despacho Económico
Ingeniería Eléctrica
Sistemas de Potencia
Optimización Matemática
Costo Incremental
Iteración Lambda
Matriz B de Kron
CENACE
XM Colombia
SCED
Publicidad