Phase shifting transformers (PSTs)

PSTs induce an equivalent voltage angle in the lines they are connected in series to by means of a voltage induction in quadrature of the line voltage, and control power flows based on that. In practise, the angle range of PSTs are given (instead of the voltage range), which is used in the power flow equations. This model is based on such PSTs.

Parameters

Set of parameters used to model phase shifting transformers as defined in the input data

namesymb.unittypedefaultdefinition
index$pst$-Int-unique index of the PST
f_bus$i$-Int-unique index of the bus to which the PST is originating from
t_bus$j$-Int-unique index of the bus to which the PST is terminating at
pst_r$r_{pst}$p.u.Real-resistance of the PST
pst_x$x_{pst}$p.u.Real-inductance of the PST
rate_a$\overline{S_{pst}}$MVAReal-long term rating of the PST
rate_b$\overline{S^{st}_{pst}}$MVAReal-short term term rating of the PST
rate_c$\overline{S^{em}_{pst}}$MVAReal-emergency rating of the PST
angle$\alpha_{pst}$radReal-current angle set point (used in power flow calculation)
pst_status$\delta_{pst}$-Int-status indicator of the PST
angmin$\underline{\alpha_{pst}}$radReal-minimum voltage angle shift
angmax$\overline{\alpha_{pst}}$radReal-maximum voltage angle shift

Variables

Optimisation variables representing PST behaviour

namesymb.unitformulationdefinition
psta$\alpha_{pst}$radACP, DCPEquivalent phase angle shift induced by PST
ppst$p_{pst,i,j}$p.u.ACP, DCPActive power flow through pst
qpst$q_{pst,i,j}$p.u.ACPReactive power flow through pst

Constraints

Flow and angle limits

Active, reactive, and apparent power limits:

\[\begin{align} - \overline{S_{pst}} \leq p_{pst,i,j} \leq \overline{S_{pst}} \\ - \overline{S_{pst}} \leq q_{pst,i,j} \leq \overline{S_{pst}} \\ p_{pst,i,j}^{2} + q_{pst,i,j}^{2} \leq \overline{S_{pst}}^{2} \end{align}\]

Range for equivalent phase angle shift:

\[\begin{align} \underline{\alpha_{pst}} \leq \alpha_{q} \leq \overline{\alpha_{pst}} \end{align}\]

PST admittance

\[\begin{align} g_{pst} = \Re{\frac{1} {r_{pst} + j \cdot x_{pst}}} \\ b_{pst} = \Im{\frac{1} {r_{pst} + j \cdot x_{pst}}} \\ \end{align}\]

Ohm's law

ACP model:

\[\begin{align} p_{pst,i,j} &= g_{pst} \cdot (v_{i})^{2} - g_{pst} \cdot v_{i} \cdot v_{j} \cdot cos(\theta_{i} - \theta_{j} - \alpha_{pst}) - b_{pst} \cdot v_{i} \cdot v_{j} \cdot sin(\theta_{i} - \theta_{j} - \alpha_{pst}) \\ q_{pst,i,j} &= -b_{pst} \cdot (v_{i})^{2} + b_{pst} \cdot v_{i} \cdot v_{j} \cdot cos(\theta_{i} - \theta_{j} - \alpha_{pst}) - g_{pst} \cdot v_{i} \cdot v_{j} \cdot sin(\theta_{i} - \theta_{j} - \alpha_{pst}) \\ \end{align}\]

\[\begin{align} p_{pst,j,i} &= g_{pst} \cdot (v_{j})^{2} - g_{pst} \cdot v_{i} \cdot v_{j} \cdot cos(\theta_{j} - \theta_{i} + \alpha_{pst}) - b_{pst} \cdot v_{i} \cdot v_{j} \cdot sin(\theta_{j} - \theta_{i} + \alpha_{pst}) \\ q_{pst,j,i} &= -b_{pst} \cdot (v_{j})^{2} + b_{pst} \cdot v_{i} \cdot v_{j} \cdot cos(\theta_{j} - \theta_{i} + \alpha_{pst}) - g_{pst} \cdot v_{i} \cdot v_{j} \cdot sin(\theta_{j} - \theta_{i} + \alpha_{pst}) \\ \end{align}\]

ACR model:

To be derived.....

DCP model:

\[\begin{align} p_{pst,i,j} &= - b_{pst} \cdot (\theta_{i} - \theta_{j} - \alpha_{pst}) \\ p_{pst,i,j} + p_{pst,j,i} &= 0 \end{align}\]