Dear Professor Zimmerman,
I think that there is a probable mistake in the auction method
(smartmarket), specifically in DC OPF. I run a smartmarket and OPF (DC and
AC) without fixed load, i.e., the two nodes system only has generators
offers and demand bids (inelastic).
When, I run AC OPF the results are correct, but with DC OPF the results may
be wrong, because the inyections and extractions electricity are very
different to the AC responses:
Node DC OPF DC OPF
Pot Lambda Pot Lambda
MW $/MW MW $/MW
1 70.8298 15 60.0394 15.00
2 -70.8298 15 -60.0067 15.016
I think the correct responses for DC OPF is 60 MW in the two nodes. Why are
the AC and DC results so different?
When, the capacity line decreases to 50 MW, the inyections and extractions
electricity are correct, but the prices changes:
Node DC OPF DC OPF
Pot Lambda Pot Lambda
MW $/MW MW $/MW
1 50 11.32 49.99 10
2 -50 18.00 -49.97 18
I think the correct responses for DC OPF is 10 and 18 ($/MW). Why are the
prices (AC and DC) in the node 1 so different?
Please, explain me What is my mistake? or Where can I modify the mistake?.
Because I need the model for my investigation.
Thank
function mpc = subasta
%% MATPOWER Case Format : Version 2
mpc.version = '2';
%% system MVA base
mpc.baseMVA = 100;
%% bus data
% bus_i type Pd Qd Gs Bs area Vm Va
baseKV zone Vmax Vmin
mpc.bus = [
1 3 0 0 0 0 1 1 0
135 1 1.05 0.95;
2 2 0 0 0 0 1 1 0
135 1 1.05 0.95;
];
%% generator data
% bus Pg Qg Qmax Qmin Vg mBase status Pmax
Pmin
mpc.gen = [
1 30 0 60 -15 1 100 1 100
0;
2 -30 0 0 0 1 100 1 0
-90;
];
%% branch data
% fbus tbus r x b rateA rateB rateC ratio
angle status
mpc.branch = [
1 2 0.001 0.01 0.0 150 150 150 0
0 1;
];
%%----- OPF Data -----%%
%% generator cost data
% 1 startup shutdown n x1 y1 ... xn
yn
% 2 startup shutdown n c(n-1) ... c0
mpc.gencost = [
1 0 0 4 0 0 25 125 50
375 100 1125;
1 0 0 4 -90 -1590 -60 -1140 -30
-600 0 0;
];
%offers.P.qty = [ 25 25 50];
%offers.P.prc = [ 5 10 15];
%bids.P.qty = [ 30 30 30 ];
%bids.P.prc = [ 20 18 15];
