Electa-Git / PowerModelsDistributionStateEstimation.jl

A Julia Package for Power System State Estimation.
BSD 3-Clause "New" or "Revised" License
32 stars 12 forks source link

polar and ivr wls on case3_unbalance "locally_infeasible" but actually not #2

Closed MartaVanin closed 4 years ago

MartaVanin commented 4 years ago

Look into ipopt settings: ipopt returns a "locally_infeasible" exit status but also a correct solution.

timmyfaraday commented 4 years ago

Is this always the case? I only got such behavior when really pushing down the variance of the uncertainty. The documentation suggests looking into a different starting point. The exploitation of the starting point given by a power flow is also something we need to look into, in my opinion.

MartaVanin commented 4 years ago

it does not happen with higher variances, but at that point the errors I generate on the fake measurement are so high that the estimation error is not negligible anymore. One thing we could do is replace the weight in the residual definition, using instead of the (square of the) variance ,100 times its amount. I think it is acceptable to do so and it won't change the result, as the ratios between the weights of more/less accurate measurements won't change, but I am not sure the literature supports this approach, that needs to be double checked. The starting point is absolutely something to discuss, we should choose an option between:

timmyfaraday commented 4 years ago

I'm not sure changing the expression of the residual is such a good idea. Additionally, a possible other startpoint could be the result of the preceeding state estimation given a number of them are performed in sequence.

MartaVanin commented 4 years ago

I added a rescaler possibility to the weights. If not explicitly declared otherwise, it defaults to 1, so the expression of the residuals remains unchanged. We can discuss this on Wednesday. Indeed, if you do SE for two succeeding time steps I think what you suggested is possibly the best approach. I reckon a better starting point than that could only (maybe) be achieved via statistical/machine learning techniques.

MartaVanin commented 4 years ago

old issue, these problems seem to be gone with the enwl dataset and the recent changes. I'll just close it