Equivalente de Red Thevenin para Estudios EMT: De Potencias de Cortocircuito a ATPDraw y PSCAD

Publicada el por fmoreno

Cuando realizamos estudios de transitorios electromagnéticos (EMT) en programas especializados como ATPDraw / EMTP o PSCAD / EMTDC (por ejemplo, para dimensionar pararrayos, evaluar sobretensiones de maniobra, energización de transformadores o coordinar aislamiento), no es viable modelar la totalidad del sistema eléctrico interconectado del país.

La práctica habitual de ingeniería consiste en acotar la frontera del estudio y representar la red aguas arriba mediante un Equivalente de Red Thevenin acoplado en secuencia positiva y secuencia cero: una fuente ideal de tensión trifásica en serie con una rama inductivo-resistiva (R_1, L_1 y R_0, L_0).

Sin embargo, el operador de la red (TSO o distribuidor) rara vez suministra los valores de inductancia en milihenrios (\text{mH}). Lo habitual es recibir una hoja técnica con:

  • Tensión nominal de línea: U_n [kV]
  • Frecuencia del sistema: f [Hz] (50\text{ Hz} o 60\text{ Hz})
  • Potencia de cortocircuito trifásica: S_{k3}^{\prime\prime} [MVA]
  • Potencia de cortocircuito monofásica a tierra: S_{k1}^{\prime\prime} [MVA]
  • Relaciones X/R de secuencia positiva (X_1/R_1) y secuencia cero (X_0/R_0)

A continuación explicamos la deducción analítica para desacoplar estos datos y cómo resolverlo de forma instantánea en CALCELEC.

1. Deducción de Impedancias de Secuencia

A. Secuencia Positiva (Z_1, R_1, X_1, L_1)

En un cortocircuito trifásico simétrico franco en barras, la corriente de falla depende exclusivamente de la impedancia de secuencia positiva del sistema:

I_{cc3}^{\prime\prime} = \frac{U_n}{\sqrt{3} \cdot Z_1}

Dado que la potencia aparente de cortocircuito trifásica se define como S_{k3}^{\prime\prime} = \sqrt{3} \cdot U_n \cdot I_{cc3}^{\prime\prime}, sustituyendo obtenemos directamente el módulo de la impedancia de secuencia positiva:

Z_1 = \frac{U_n^2}{S_{k3}^{\prime\prime}}

A partir de la relación (X_1 / R_1):

\theta_1 = \arctan\left(\frac{X_1}{R_1}\right) R_1 = Z_1 \cdot \cos(\theta_1) X_1 = Z_1 \cdot \sin(\theta_1) L_1 = \frac{X_1}{2\pi f} \times 1000 \quad [\text{mH}]

B. Secuencia Cero (Z_0, R_0, X_0, L_0)

El cortocircuito monofásico a tierra involucra la trayectoria de retorno por tierra y neutro. De acuerdo con la teoría de componentes simétricas:

I_{cc1}^{\prime\prime} = \frac{3 \cdot \frac{U_n}{\sqrt{3}}}{2 Z_1 + Z_0} = \frac{\sqrt{3} \cdot U_n}{2 Z_1 + Z_0}

La potencia de cortocircuito monofásica reportada por las empresas eléctricas se formaliza como:

S_{k1}^{\prime\prime} = \sqrt{3} \cdot U_n \cdot I_{cc1}^{\prime\prime} = \frac{3 \cdot U_n^2}{2 Z_1 + Z_0}

Despejando la impedancia de secuencia cero (Z_0):

Z_0 = \frac{3 \cdot U_n^2}{S_{k1}^{\prime\prime}} - 2 Z_1

Condición de consistencia física: Para que la red tenga significado físico real, Z_0 debe ser mayor que cero (Z_0 > 0). Si S_{k1}^{\prime\prime} > 1.5 \cdot S_{k3}^{\prime\prime}, la impedancia resultante resultaría matemáticamente negativa, lo cual indica que los datos de entrada suministrados por la empresa eléctrica presentan inconsistencias.

Una vez obtenido Z_0 y conociendo la relación (X_0 / R_0):

\theta_0 = \arctan\left(\frac{X_0}{R_0}\right) R_0 = Z_0 \cdot \cos(\theta_0) X_0 = Z_0 \cdot \sin(\theta_0) L_0 = \frac{X_0}{2\pi f} \times 1000 \quad [\text{mH}]

2. Cómo Implementar el Equivalente en ATPDraw y PSCAD

Una vez calculados R_1, L_1, R_0, L_0, ¿cómo se configuran en los softwares de simulación transitoria?

En ATPDraw:

  • Puedes utilizar el bloque estándar de fuente trifásica con impedancia interna acoplada (Type 14 en conjunto con una rama acoplada RL Coupled / R-L Matrix).


  • La conversión de componentes simétricas a componentes de fase (R_s, R_m, L_s, L_m) se realiza mediante:


R_{\text{propia}} (R_s) = \frac{2 R_1 + R_0}{3}, \quad R_{\text{mutua}} (R_m) = \frac{R_0 - R_1}{3} L_{\text{propia}} (L_s) = \frac{2 L_1 + L_0}{3}, \quad L_{\text{mutua}} (L_m) = \frac{L_0 - L_1}{3}

En PSCAD:

  • En el componente estándar de fuente trifásica (Three-Phase Voltage Source Model 1 / Model 2), selecciona el modo de especificación por impedancias de secuencia (Zero sequence impedance entered) e ingresa directamente R_1, L_1, R_0, L_0 sin necesidad de realizar conversiones manuales adicionales.

3. Uso del Módulo en CALCELEC

El módulo Equivalente de Red (EMT) de CALCELEC (ubicado en el menú QuickCalc) te evita tener que programar scripts o recurrir a hojas de cálculo con riesgo de errores de trigonometría:

  1. Introduce:
    • Tensión nominal de barra U_n (kV) y frecuencia de red (50 o 60\text{ Hz}).
    • S_{k3}^{\prime\prime} y S_{k1}^{\prime\prime} en MVA.
    • Factores X_1/R_1 y X_0/R_0.
  2. Validación automática: El motor comprueba que Z_0 > 0 antes de procesar los resultados.
  3. Salida lista para modelado: Visualizas de inmediato Z_1, R_1, X_1, L_1 y Z_0, R_0, X_0, L_0 con cuatro decimales de precisión.

4. Ejemplo Numérico Real

Consideremos un punto de interconexión en una subestación de 220\text{ kV}, 50\text{ Hz}:

  • S_{k3}^{\prime\prime} = 5000\text{ MVA}, con X_1/R_1 = 10
  • S_{k1}^{\prime\prime} = 4200\text{ MVA}, con X_0/R_0 = 3

Al procesar los datos en CALCELEC obtenemos:

  • Secuencia Positiva:
    Z_1 = 9.6800\ \Omega \rightarrow R_1 = 0.9632\ \Omega, X_1 = 9.6319\ \Omega, L_1 = 30.6593\text{ mH}
  • Secuencia Cero:
    Z_0 = 15.2114\ \Omega \rightarrow R_0 = 4.8103\ \Omega, X_0 = 14.4308\ \Omega, L_0 = 45.9348\text{ mH}

5. Conclusión

Garantizar la exactitud del equivalente Thevenin de la red es el cimiento de cualquier estudio de transitorios electromagnéticos de alta calidad.

¿Preparando tu modelo en ATPDraw o PSCAD?
Obtén tus parámetros equivalentes de red al instante con CALCELEC.com.

Deja una respuesta

Entradas relacionadas