Since this script uses JuMP to implement the model, some changes were made to improve the performance and take advantage of the JuMP.
The use of JuMP's expressions avoids unnecessary variable and reduce the size of the problem, increasing the performance.
-
All present value costs definitions were substituted by expressions, even the objective function variable
$c^{TPV}$ is an expression in this version. -
The variables
$f_{srktb}^{l}$ and$g_{srtb}^{tr}$ were also replaced by expressions since they are defined by the linearization of the losses.
The paper uses
The paper allows the replacement of planned investment, but this only increases the costs. More recent papers avoid that (Eq. 22 - 24).
During the development of this model, we cannot find the Wind generators' data. Thus, the data used is from this reference.
In this model, we divide the generators into two different sets.
The
We also simulate the case where the renewables are not injecting power into the distribution system.
This consideration adds the index
We simplify the radiality constraints. But it is not expected to change the results.
Some other minor changes were made, but it is not expected to cause differences between
this implementation and the model proposed by the paper. See the src/model.jl
file.