Skip to content

Commit

Permalink
Stepuptransformers reactive bound fix (#299)
Browse files Browse the repository at this point in the history 
* checks the step-up transformer impedance is compatible with the reactive limit bounds

* sets INFINITE_REACTIVE_LIMIT to 99999
  • Loading branch information
@CBiasuzzi @mathbagu
CBiasuzzi authored and mathbagu committed Nov 4, 2017
1 parent e13e9af commit 849b8f7
Showing 1 changed file with 33 additions and 7 deletions.
40 changes: 33 additions & 7 deletions .../src/main/java/eu/itesla_project/iidm/ddb/eurostag/EurostagStepUpTransformerInserter.java
View file
Original file line number Diff line number Diff line change
@@ -1,15 +1,13 @@
/**
* Copyright (c) 2016, All partners of the iTesla project (http://www.itesla-project.eu/consortium)
* Copyright (c) 2017, RTE (http://www.rte-france.com)
* This Source Code Form is subject to the terms of the Mozilla Public
* License, v. 2.0. If a copy of the MPL was not distributed with this
* file, You can obtain one at http://mozilla.org/MPL/2.0/.
*/
package eu.itesla_project.iidm.ddb.eurostag;

import com.powsybl.computation.ComputationManager;
import eu.itesla_project.iidm.ddb.eurostag.model.PowerFlow;
import eu.itesla_project.iidm.ddb.eurostag.model.StateVariable;
import eu.itesla_project.iidm.ddb.eurostag.model.TransformerModel;
import com.powsybl.iidm.network.*;
import com.powsybl.iidm.network.ReactiveCapabilityCurve.Point;
import com.powsybl.iidm.network.util.Identifiables;
Expand All @@ -18,6 +16,9 @@
import com.powsybl.loadflow.LoadFlowFactory;
import com.powsybl.loadflow.LoadFlowParameters;
import com.powsybl.loadflow.LoadFlowResult;
import eu.itesla_project.iidm.ddb.eurostag.model.PowerFlow;
import eu.itesla_project.iidm.ddb.eurostag.model.StateVariable;
import eu.itesla_project.iidm.ddb.eurostag.model.TransformerModel;
import org.slf4j.Logger;
import org.slf4j.LoggerFactory;

Expand All @@ -39,7 +40,7 @@ public class EurostagStepUpTransformerInserter {

private final static float SB = 100;
private final static float ACCEPTABLE_VOLTAGE_DIFF = 0.1f;
private final static float INFINITE_REACTIVE_LIMIT = 9999f;
private final static float INFINITE_REACTIVE_LIMIT = 99999f;

public enum InsertionStatus {
OK("ok"),
Expand Down Expand Up @@ -301,7 +302,7 @@ private static void fillGeneratorState(Generator g, StateVariable sv) {
}

private static Generator moveGenerator(Generator srcGen, StateVariable srcSv, VoltageLevel destVl, Bus destBus, boolean connected,
Function<StateVariable, StateVariable> fct, EurostagStepUpTransformerConfig config) {
Function<StateVariable, StateVariable> fct, EurostagStepUpTransformerConfig config, TwoWindingsTransformer twt) {

fillGeneratorState(srcGen, srcSv);

Expand Down Expand Up @@ -367,6 +368,17 @@ private static Generator moveGenerator(Generator srcGen, StateVariable srcSv, Vo
float newP = (float) -fct.apply(new StateVariable(-point.getP(), 0, srcSv.u, srcSv.theta)).p;
float newMinQ = config.isNoReactiveLimits() ? -INFINITE_REACTIVE_LIMIT : (float) -fct.apply(new StateVariable(srcSv.p, -point.getMinQ(), srcSv.u, srcSv.theta)).q;
float newMaxQ = config.isNoReactiveLimits() ? INFINITE_REACTIVE_LIMIT : (float) -fct.apply(new StateVariable(srcSv.p, -point.getMaxQ(), srcSv.u, srcSv.theta)).q;

//test if the step-up transformer impedance is compatible with the reactive limit bounds: the maximun nominal power with a very high
//short-circuit power should be higher than the reactive bounds. If not, remove the bounds
float sMax = computeMaxRate(twt);
if ((sMax < Math.abs(newMinQ)) || (sMax < Math.abs(newMaxQ))) {
LOGGER.warn("Infinite limit mode activated for '{}': either abs({}) or abs({}) higher than {} (maximum nominal power with a short-circuit power of 40%)",
lvGen.getId(), newMinQ, newMaxQ, sMax);
newMinQ = -INFINITE_REACTIVE_LIMIT;
newMaxQ = INFINITE_REACTIVE_LIMIT;
}

LOGGER.trace("Resizing reactive limits of '{}': [{}, {}] -> [{}, {}]",
lvGen.getId(), point.getMinQ(), point.getMaxQ(), newMinQ, newMaxQ);
rcca.beginPoint()
Expand Down Expand Up @@ -455,13 +467,27 @@ public static InsertionStatus insert(Generator hvGen, TG tg, IdDictionary auxDic
Function<StateVariable, StateVariable> fct = sv -> toLvGenPf(transformerModel, sv, hvAuxPf, lvAuxPf);

StateVariable hvGenSv = new StateVariable();
moveGenerator(hvGen, hvGenSv, lvVl, lvBus, hvGenBus != null, fct, config);
moveGenerator(hvGen, hvGenSv, lvVl, lvBus, hvGenBus != null, fct, config, twt);

stateBefore.injections.put(twt.getId(), new PowerFlow((float) (hvGenSv.p + hvAuxPf.p), (float) (hvGenSv.q + hvAuxPf.q)));

return InsertionStatus.OK;
}

/**
* Compute a reasonable maximum nominal power ('rate') for a TwoWindingsTransformer
* with a very high short-circuit power of 40%
*
* @param twt two winding transformer
* @return maximum nominal power
*/
private static float computeMaxRate(TwoWindingsTransformer twt) {
Objects.requireNonNull(twt);
VoltageLevel u2 = twt.getTerminal2().getVoltageLevel();
float z = SV.getX(twt);
return 0.4f * (float) Math.pow(u2.getNominalV(), 2) / z;
}

public static InsertionStatus insert(Generator g, Path ddbFile, IdDictionary auxDict, EurostagStepUpTransformerConfig config, StateBefore stateBefore) throws IOException {
if (!Files.exists(ddbFile) && !Files.isRegularFile(ddbFile)) {
throw new IllegalArgumentException(ddbFile + " must exist and be a file");
Expand Down Expand Up @@ -612,7 +638,7 @@ public void visitThreeWindingsTransformer(ThreeWindingsTransformer transformer,
boolean connected = twt.getTerminal1().isConnected() && twt.getTerminal2().isConnected() && lvGen.getTerminal().isConnected();

StateVariable lvGenSv = new StateVariable();
moveGenerator(lvGen, lvGenSv, hvVl, hvBus, connected, fct, config);
moveGenerator(lvGen, lvGenSv, hvVl, hvBus, connected, fct, config, twt);
}

for (Load aux : auxLs) {
Expand Down

0 comments on commit 849b8f7

Please sign in to comment.