Correct inconsistency between disconnected bus and branch status

Author: iagochavarry

src/pwf2pm/branch.jl

@@ -25,6 +25,21 @@ function _create_dict_dshl(data::Dict)
     return dshl_dict
 end
 
+function _handle_br_status!(pm_data::Dict, sub_data::Dict, branch::Dict)
+    f_bus = sub_data["f_bus"]
+    t_bus = sub_data["t_bus"]
+
+    if pm_data["bus"]["$f_bus"]["bus_type"] != 4 && pm_data["bus"]["$t_bus"]["bus_type"] != 4
+        if branch["STATUS"] == branch["OPENING FROM BUS"] == branch["OPENING TO BUS"] == 'L'
+            sub_data["br_status"] = 1
+        end
+    else
+        sub_data["br_status"] = 0
+    end
+    
+    return sub_data
+end
+
 function _pwf2pm_branch!(pm_data::Dict, pwf_data::Dict, branch::Dict; add_control_data::Bool=false)
     sub_data = Dict{String,Any}()
 
@@ -46,11 +61,7 @@ function _pwf2pm_branch!(pm_data::Dict, pwf_data::Dict, branch::Dict; add_contro
     sub_data["angmax"] = 60.0 # PowerModels.jl standard
     sub_data["transformer"] = false
 
-    if branch["STATUS"] == branch["OPENING FROM BUS"] == branch["OPENING TO BUS"] == 'L'
-        sub_data["br_status"] = 1
-    else
-        sub_data["br_status"] = 0
-    end
+    _handle_br_status!(pm_data, sub_data, branch)
 
     if add_control_data
         sub_data["control_data"] = Dict{String,Any}()

src/pwf2pm/transformer.jl

@@ -54,11 +54,7 @@ function _pwf2pm_transformer!(pm_data::Dict, pwf_data::Dict, branch::Dict; add_c
     sub_data["angmax"] = 60.0 # PowerModels.jl standard
     sub_data["transformer"] = true
 
-    if branch["STATUS"] == branch["OPENING FROM BUS"] == branch["OPENING TO BUS"] == 'L'
-        sub_data["br_status"] = 1
-    else
-        sub_data["br_status"] = 0
-    end
+    _handle_br_status!(pm_data, sub_data, branch)
 
     n = 0 # count(x -> x["f_bus"] == sub_data["f_bus"] && x["t_bus"] == sub_data["t_bus"], values(pm_data["branch"])) 
     sub_data["source_id"] = ["transformer", sub_data["f_bus"], sub_data["t_bus"], 0, "0$(n + 1)", 0]

test/data/pwf/test_system.pwf

@@ -22,7 +22,7 @@ LFPO     .1
 99999
 DBAR
 (Num)OETGb(   nome   )Gl( V)( A)( Pg)( Qg)( Qn)( Qm)(Bc  )( Pl)( Ql)( Sh)Are(Vf)M(1)(2)(3)(4)(5)(6)(7)(8)(9)(10
-   21 L1 ABus 1        A10408.88 19.6 -3.2 -130130.4                       11000                                   
+   21 D1 ABus 1        A10408.88 19.6 -3.2 -130130.4                       11000                                   
    22 L2 BBus 2        A1025 9.339.2 -4.6 -99999999.                       11000                                   
    23 L2 BBus 3        A10257.5325.5 -11.4-99999999.                       11000                                   
    24 L  CBus 4        A1060 8.4                           34.6 6.9        11000                                   

test/test_pwf_to_powermodels.jl

@@ -147,14 +147,21 @@
                     @test solution["multiinfrastructure"] == false
                     @test solution["multinetwork"]        == false
                     @test solution["per_unit"]            == true
-                    @test length(solution["bus"])         == 9
+                    @test length(solution["bus"])         == 8
                 end
 
             end
 
         end
 
+        @testset "Bus Disconnected - Branch Status" begin
+            # If a branch whose extremity is in a disconnected bus (bus_type = 4)
+            # the branch status must be off
 
+            file = open(joinpath(@__DIR__,"data/pwf/test_system.pwf"))
+            pm_data = PWF.parse_pwf_to_powermodels(file)
+            @test pm_data["branch"]["5"]["br_status"] == 0
+        end
 
         @testset "Power Flow results" begin
             filenames = ["3bus", "9bus"]