Sensitivity analysis

Introduction

The sensitivity analysis module is dedicated to computing the linearized impact of small network variations on the state variables of some equipments.

A sensitivity value is the numerical estimation of the partial derivative of the observed function with respect to the variable of impact. The sensitivity analysis can also be seen as the computation of partial derivatives on the network model. For example, it may be used to know, among a group of selected lines, which are the most impacted by a change in a generator production or a change of tap on a phase tap changer. The user story about RSC capacity calculation provides an example of application of the sensitivity analysis.

Sensitivity analysis inputs

Network

The first input for the sensitivity analysis module is an IIDM network.

Sensitivity factors

Aside from providing an input network, it is necessary to specify which equipments are going to be studied:

  • what impacted equipments are selected to be monitored (lines for example)
  • according to a change on which equipment (a generator’s production or a group of generator’s production, or the tap position of a phase tap changer, etc.)

It is also necessary to specify which quantity is being observed: the active power or the current on the monitored equipments, the voltage at a bus.

This set of information constitutes the sensitivity factors (SensitivityFactor). These factors correspond to the definition of the expected partial derivatives to be extracted from the input network. A standard sensitivity analysis input thus comprises a list of sensitivity factors, each one constituted of:

  • a sensitivity variable (the variable of impact) which type is defined by a SensitivityVariableType.
  • a sensitivity function (the observed function) which type is defined by a SensitivityFunctionType.
  • a contingency context. Usually we compute the impact of an injection increase on a branch flow or current, the impact of a shift of a phase tap changer on a branch flow or current or the impact of a voltage target increase on a bus voltage.

A sensitivity variable represents a change on a equipment or on a group of equipments. The supported variable types are:

  • Use INJECTION_ACTIVE_POWER to model a change on the production of a generator or on a group of generators, on the consumption of a load or on a group of loads or on GLSK (for Generation and Load Shift keys) that describes a linear combination of power injection shifts on generators and loads. The variable increase is in MW.
  • Use TRANSFORMER_PHASE to model the change of the tap position of a phase tap changer of a two windings transformer. The increase is in degree.
  • Use BUS_TARGET_VOLTAGE to model an increase of the voltage target of a generator, a static var compensator, a two or three windings transformer, a shunt compensator or a VSC converter station. The increase is in KV.
  • Use HVDC_LINE_ACTIVE_POWER to model the change of the active power set point of an HVDC line. The increase is in MW.
  • Use TRANSFORMER_PHASE_1, TRANSFORMER_PHASE_2 or TRANSFORMER_PHASE_3 to model the change of the tap position of a phase tap changer of a three windings transformer that contains several phase tap changers.

The supported sensitivity function types, related to the equipment to monitor, are:

  • Use BRANCH_ACTIVE_POWER_1 and BRANCH_ACTIVE_POWER_2 if you want to monitor the active power in MW of a network branch (lines, two windings transformer, dangling lines, etc.). Use 1 for side 1 and 2 for side 2. In case of a three windings transformer, use BRANCH_ACTIVE_POWER_3 to monitor the active power in MW of the leg 3 (network side).
    • Use BRANCH_REACTIVE_POWER_1 and BRANCH_REACTIVE_POWER_2 if you want to monitor the reactive power in MVar of a network branch (lines, two-winding transformer, dangling lines, etc.). Use 1 for side 1 and 2 for side 2. In case of a three-winding transformer, use BRANCH_REACTIVE_POWER_3 to monitor the reactive power in MVar of the leg 3 (network side).
  • Use BRANCH_CURRENT_1 and BRANCH_CURRENT_2 if you want to monitor the current in A of a network branch (lines, two windings transformer, dangling lines, etc.). Use 1 for side 1 and use 2 for side 2. In case of a three windings transformer, use BRANCH_CURRENT_3 to monitor the current in A of the leg 3 (network side).
  • BUS_VOLTAGE if you want to monitor the voltage in KV of a specific network bus.

A sensitivity variable can group some equipments and has to be modeled as a variable set. In a SensitivityVariableSet, we have a list of individual variables, each one with a weight (called WeightedSensitivityVariable). We use variable sets to model what it commonly called GLSK.

How to provide the sensitivity factors input

The sensitivity factors may be created directly through Java code, or be provided to PowSyBl via a JSON file. This file should contain a list of JSON objects, each one representing a sensitivity factor. The example below shows how to write a JSON file to perform a sensitivity analysis on the active power through a line, with respect to a GLSK and to an injection on the network.

[ {
  "functionType" : "BRANCH_ACTIVE_POWER_1",
  "functionId" : "l45",
  "variableType" : "INJECTION_ACTIVE_POWER",
  "variableId" : "glsk",
  "variableSet" : true,
  "contingencyContextType" : "ALL"
}, {
  "functionType" : "BRANCH_ACTIVE_POWER_1",
  "functionId" : "l12",
  "variableType" : "INJECTION_ACTIVE_POWER",
  "variableId" : "g2",
  "variableSet" : false,
  "contingencyContextType" : "ALL"
} ]

Contingencies

The sensitivity analysis may also take, optionally, a list of contingencies as an input. When contingencies are provided, the sensitivity values shall be calculated on the network at state N, but also after the application of each contingency. The contingencies are provided in the same way than for the security analysis. This then constitutes a systematic sensitivity analysis.

{
  "type" : "default",
  "version" : "1.0",
  "name" : "default",
  "contingencies" : [ {
    "id" : "l34",
    "elements" : [ {
      "id" : "l34",
      "type" : "BRANCH"
    } ]
  } ]
}

At the moment the only available sensitivity simulator officially compatible with PowSyBl is the one available through OpenLoadFlow. In this case, the network is provided only once in state N, and then all the calculations are done successively by modifying the Jacobian matrix directly in the solver based on the contingencies input. The network is thus loaded only once, which improves performance.

Sensitivity analysis outputs

Sensitivity values

The outputs of the sensitivity analysis are called sensitivity values. A sensitivity value represents an elementary result given a sensitivity factor and a contingency, and contains:

  • The actual value of the partial derivative
  • The reference value of the function at linearization point in case of contingency context NONE or in its post-contingency state in case of a factor associated to a contingency.

These results may be serialized in JSON format.

Example of interpretation

Let’s imagine that one wants to compute the impact of an increase of active power generation of the generator G on the branch B. The sensitivity analysis input will contain one sensitivity factor, with sensitivity function type BRANCH_ACTIVE_POWER_1 and sensitivity variable type INJECTION_ACTIVE_POWER, and we do not provide any input contingencies.

After the computation, let us consider that the values of the three elements of the sensitivity result are:

  • a value of -0.05 for the partial derivative
  • a variable reference value of 150
  • a function reference value of 265

This can be interpreted in the following way:

  • an increase of 100 MW on generator G may be approximated on branch B as a 5MW decrease of the active flow from side 1 to side 2
  • the initial generation on generator G is 150MW
  • the initial active flow on branch B is 265MW from side 1 to side 2

Implementations

The following sensitivity analysis implementations are supported:

Going further

To go further about the sensitivity analysis, check the following content: