Validation 04 - Two Generator, Two Bus, Two Load, One Branch, LMP Change

In [38]:

%matplotlib inline

In [39]:

import psst

In [40]:

from psst.case import read_matpower
from psst.network import create_network

Validation of case 1

In [41]:

case = read_matpower('./cases/case4.m')

In [42]:

network = create_network(case, prog='neato')
network.draw()
../../_images/notebooks_validation_Validation04-TwoGeneratorTwoBusTwoLoadOneBranchLMPChange_6_0.png 
In [43]:

case

Out[43]:

<psst.case.PSSTCase(name=case5, Generators=2, Buses=2, Branches=1)>

In [63]:

case.bus

Out[63]:
TYPE PD QD GS BS AREA VM VA BASEKV ZONE VMAX VMIN
Bus1 3 100 131.47 0 0 1 1 0 230 1 1.1 0.9
Bus2 2 300 0.00 0 0 1 1 0 230 1 1.1 0.9 
In [62]:

case.branch

Out[62]:
F_BUS T_BUS BR_R BR_X BR_B RATE_A RATE_B RATE_C TAP SHIFT BR_STATUS ANGMIN ANGMAX
0 Bus1 Bus2 0.00281 0.0281 0.00712 800 800 800 0 0 1 -360 360 
In [44]:

case.gen

Out[44]:
GEN_BUS PG QG QMAX QMIN VG MBASE GEN_STATUS PMAX PMIN PC1 PC2 QC1MIN QC1MAX QC2MIN QC2MAX RAMP_AGC RAMP_10 RAMP_30 RAMP_Q APF
GenCo0 Bus1 200 0 30 -30 1 100 1 200 0 0 0 0 0 0 0 0 0 0 0 0
GenCo1 Bus2 500 0 30 -30 1 100 1 500 0 0 0 0 0 0 0 0 0 0 0 0 
In [45]:

case.gencost

Out[45]:
MODEL STARTUP SHUTDOWN NCOST COST_1 COST_0
GenCo0 1 0 0 2 10 0
GenCo1 1 0 0 2 14 0 
In [46]:

case.load

Out[46]:
Bus1 Bus2
0 100.0 300.0 
In [47]:

from psst.model import build_model

In [48]:

model = build_model(case)

In [49]:

model

Out[49]:

<psst.model.PSSTModel(status=None)>

In [50]:

model.solve(solver='cbc', verbose=True)

Welcome to the CBC MILP Solver
Version: 2.9.6
Build Date: May 27 2016

command line - /usr/local/bin/cbc -mipgap 0.01 -printingOptions all -import /var/folders/wk/lcf0vgd90bx0vq1873tn04knk_djr3/T/tmpYO7VdM.pyomo.lp -import -stat=1 -solve -solu /var/folders/wk/lcf0vgd90bx0vq1873tn04knk_djr3/T/tmpYO7VdM.pyomo.soln (default strategy 1)
No match for mipgap - ? for list of commands
No match for 0.01 - ? for list of commands
Option for printingOptions changed from normal to all
Current default (if $ as parameter) for import is /var/folders/wk/lcf0vgd90bx0vq1873tn04knk_djr3/T/tmpYO7VdM.pyomo.lp
Presolve 7 (-41) rows, 13 (-27) columns and 22 (-77) elements
Statistics for presolved model
Original problem has 2 integers (2 of which binary)


Problem has 7 rows, 13 columns (8 with objective) and 22 elements
There are 8 singletons with objective
Column breakdown:
10 of type 0.0->inf, 2 of type 0.0->up, 0 of type lo->inf,
1 of type lo->up, 0 of type free, 0 of type fixed,
0 of type -inf->0.0, 0 of type -inf->up, 0 of type 0.0->1.0
Row breakdown:
0 of type E 0.0, 0 of type E 1.0, 0 of type E -1.0,
3 of type E other, 0 of type G 0.0, 0 of type G 1.0,
0 of type G other, 2 of type L 0.0, 0 of type L 1.0,
2 of type L other, 0 of type Range 0.0->1.0, 0 of type Range other,
0 of type Free
Continuous objective value is 4800 - 0.00 seconds
Cgl0004I processed model has 7 rows, 13 columns (0 integer (0 of which binary)) and 20 elements
Cbc3007W No integer variables - nothing to do
Cuts at root node changed objective from 4800 to -1.79769e+308
Probing was tried 0 times and created 0 cuts of which 0 were active after adding rounds of cuts (0.000 seconds)
Gomory was tried 0 times and created 0 cuts of which 0 were active after adding rounds of cuts (0.000 seconds)
Knapsack was tried 0 times and created 0 cuts of which 0 were active after adding rounds of cuts (0.000 seconds)
Clique was tried 0 times and created 0 cuts of which 0 were active after adding rounds of cuts (0.000 seconds)
MixedIntegerRounding2 was tried 0 times and created 0 cuts of which 0 were active after adding rounds of cuts (0.000 seconds)
FlowCover was tried 0 times and created 0 cuts of which 0 were active after adding rounds of cuts (0.000 seconds)
TwoMirCuts was tried 0 times and created 0 cuts of which 0 were active after adding rounds of cuts (0.000 seconds)

Result - Optimal solution found

Objective value:                4800.00000000
Enumerated nodes:               0
Total iterations:               0
Time (CPU seconds):             0.00
Time (Wallclock seconds):       0.00

Total time (CPU seconds):       0.00   (Wallclock seconds):       0.01

Input data

In [51]:

import pandas as pd

In [52]:

pd.DataFrame(case.gen['PMAX'])

Out[52]:
PMAX
GenCo0 200
GenCo1 500 
In [53]:

case.load

Out[53]:
Bus1 Bus2
0 100.0 300.0

Model Results

In [54]:

model.results.unit_commitment

Out[54]:
GenCo0 GenCo1
0 1 1 
In [55]:

model.results.power_generated

Out[55]:
GenCo0 GenCo1
0 200 200 
In [56]:

model.results.production_cost

Out[56]:

4800

In [57]:

case.gencost.loc['GenCo0', 'COST_1']*200 + case.gencost.loc['GenCo1', 'COST_1']*200

Out[57]:

4800

In [58]:

model.results.line_power

Out[58]:
0
0 100 
In [59]:

model.results.lmp

Out[59]:
Bus1 Bus2
0 14 14 
In [60]:

from psst.plot import line_power

In [61]:

line_power(network, model.results, hour=0)

Out[61]:

<matplotlib.axes._subplots.AxesSubplot at 0x11128aed0>
../../_images/notebooks_validation_Validation04-TwoGeneratorTwoBusTwoLoadOneBranchLMPChange_29_1.png