Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Issue 19 #20

Open
wants to merge 13 commits into
base: master
Choose a base branch
from
+1,826 −43
Open

Issue 19 #20

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
13 commits
Select commit Hold shift + click to select a range
302b45e
fix error creating branches of type PowerTransformer
martinmoraga Feb 23, 2024
1a6940a
added new powerflow example using the 9bus system
martinmoraga Feb 23, 2024
db69164
MOD: np.complex is depricated since numpy 1.20.
priyanka1ag Apr 26, 2024
45ab7f1
MOD: Extracting capacitance of line from CIM. Additional helper routi…
priyanka1ag Apr 26, 2024
9a237e9
function for updating measurement value of existing measurement object
priyanka1ag Jul 31, 2024
c21db34
MOD: Included necessary changes to include PMU node current injection…
priyanka1ag Sep 23, 2024
d80913c
MOD: it is charging susceptance, not capacitance from CIM
priyanka1ag Nov 4, 2024
156f2c6
MOD: Added measurements in actuals. Previously measurements were avai…
priyanka1ag Nov 12, 2024
47c1a90
MOD: to keep it consistent to power flow results, the actuals in meas…
priyanka1ag Nov 12, 2024
2501b7b
MOD: the base voltage and power are in three phase, but power flow an…
priyanka1ag Nov 12, 2024
c20ff45
MOD: Corrected the measurement standard deviation calculation.
priyanka1ag Nov 22, 2024
acaa5e2
MOD: Measurement covariance matrix
priyanka1ag Nov 22, 2024
459daae
MOD: function for getting measurements in actuals
priyanka1ag Nov 25, 2024
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
61 changes: 61 additions & 0 deletions examples/quickstart/run_nv_powerflow_WSCC9.py
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,61 @@
import logging
from pathlib import Path
from pyvolt import network
from pyvolt import nv_powerflow
import numpy
import cimpy
import os
from pyvolt.network import BusType


logging.basicConfig(filename='run_nv_powerflow.log', level=logging.INFO, filemode='w')

#python starts as module in subdirectory, 2 folders up to set the new path
this_file_folder = Path(__file__).parents[2]
p = str(this_file_folder)+"/examples/sample_data/WSCC9"
xml_path = Path(p)


xml_files = [os.path.join(xml_path, "WSCC-09_Dyn_Fourth_EQ.xml"),
os.path.join(xml_path, "WSCC-09_Dyn_Fourth_SV.xml"),
os.path.join(xml_path, "WSCC-09_Dyn_Fourth_TP.xml")]

# Read cim files and create new network.System object
res = cimpy.cim_import(xml_files, "cgmes_v2_4_15")
system = network.System()
base_apparent_power = 100 # MW
system.load_cim_data(res['topology'], base_apparent_power)

# Set node 1 as slack
system.nodes[4].type=BusType.SLACK

# Execute power flow analysis
results_pf, num_iter = nv_powerflow.solve(system)

# Print node voltages
print("Powerflow converged in " + str(num_iter) + " iterations.\n")
print("Results: \n")
voltages = []
for node in results_pf.nodes:
print('Magnitude {}={}, Phase {}={}'.format(node.topology_node.uuid, numpy.abs(node.voltage_pu), node.topology_node.uuid, numpy.angle(node.voltage_pu, deg=True)))
voltages.append(node.voltage_pu)

"""
# validation against matpower
matpower_voltages = {
"BUS1": [1.040, 0.000],
"BUS2": [1.025, 9.693],
"BUS3": [1.025, 4.881],
"BUS4": [0.996, -2.306],
"BUS5": [0.966, -3.737],
"BUS6": [1.004, 2.107],
"BUS7": [0.979, 0.836],
"BUS8": [0.997, 3.974],
"BUS9": [0.951, 4.138]
}

# assertion
for idx, node in enumerate(results_pf.nodes):
print(matpower_voltages[node.topology_node.uuid][0] - numpy.round(numpy.abs(node.voltage_pu),3))
print(matpower_voltages[node.topology_node.uuid][1] - numpy.round(numpy.angle(node.voltage_pu, deg=True),3))
"""
1,268 changes: 1,268 additions & 0 deletions examples/sample_data/WSCC9/WSCC-09_Dyn_Fourth_EQ.xml
View file
Open in desktop

Large diffs are not rendered by default.

78 changes: 78 additions & 0 deletions examples/sample_data/WSCC9/WSCC-09_Dyn_Fourth_SV.xml
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,78 @@
<?xml version="1.0" encoding="utf-8"?>
<rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#" xmlns:cim="http://iec.ch/TC57/2012/CIM-schema-cim16#" xmlns:md="http://iec.ch/TC57/61970-552/ModelDescription/1#" xmlns:entsoe="http://entsoe.eu/Secretariat/ProfileExtension/2#">
<cim:SvVoltage rdf:ID="BUS3_5B20">
<cim:SvVoltage.v>14.145000</cim:SvVoltage.v>
<cim:SvVoltage.angle>4.664751</cim:SvVoltage.angle>
<cim:SvVoltage.TopologicalNode rdf:resource="#BUS3"/>
</cim:SvVoltage>
<cim:SvVoltage rdf:ID="BUS5_5B20">
<cim:SvVoltage.v>228.995097</cim:SvVoltage.v>
<cim:SvVoltage.angle>-3.988800</cim:SvVoltage.angle>
<cim:SvVoltage.TopologicalNode rdf:resource="#BUS5"/>
</cim:SvVoltage>
<cim:SvVoltage rdf:ID="BUS6_5B20">
<cim:SvVoltage.v>232.910495</cim:SvVoltage.v>
<cim:SvVoltage.angle>-3.687390</cim:SvVoltage.angle>
<cim:SvVoltage.TopologicalNode rdf:resource="#BUS6"/>
</cim:SvVoltage>
<cim:SvVoltage rdf:ID="BUS4_5B20">
<cim:SvVoltage.v>235.931330</cim:SvVoltage.v>
<cim:SvVoltage.angle>-2.216780</cim:SvVoltage.angle>
<cim:SvVoltage.TopologicalNode rdf:resource="#BUS4"/>
</cim:SvVoltage>
<cim:SvVoltage rdf:ID="BUS1_5B20">
<cim:SvVoltage.v>17.160000</cim:SvVoltage.v>
<cim:SvVoltage.angle>-0e+000</cim:SvVoltage.angle>
<cim:SvVoltage.TopologicalNode rdf:resource="#BUS1"/>
</cim:SvVoltage>
<cim:SvVoltage rdf:ID="BUS2_5B20">
<cim:SvVoltage.v>18.450000</cim:SvVoltage.v>
<cim:SvVoltage.angle>9.280005</cim:SvVoltage.angle>
<cim:SvVoltage.TopologicalNode rdf:resource="#BUS2"/>
</cim:SvVoltage>
<cim:SvVoltage rdf:ID="BUS7_5B20">
<cim:SvVoltage.v>235.926956</cim:SvVoltage.v>
<cim:SvVoltage.angle>3.719701</cim:SvVoltage.angle>
<cim:SvVoltage.TopologicalNode rdf:resource="#BUS7"/>
</cim:SvVoltage>
<cim:SvVoltage rdf:ID="BUS9_5B20">
<cim:SvVoltage.v>237.441178</cim:SvVoltage.v>
<cim:SvVoltage.angle>1.966716</cim:SvVoltage.angle>
<cim:SvVoltage.TopologicalNode rdf:resource="#BUS9"/>
</cim:SvVoltage>
<cim:SvVoltage rdf:ID="BUS8_5B20">
<cim:SvVoltage.v>233.652994</cim:SvVoltage.v>
<cim:SvVoltage.angle>0.727536</cim:SvVoltage.angle>
<cim:SvVoltage.TopologicalNode rdf:resource="#BUS8"/>
</cim:SvVoltage>
<cim:SvPowerFlow rdf:ID="LOAD8-sv">
<cim:SvPowerFlow.p>100.000000</cim:SvPowerFlow.p>
<cim:SvPowerFlow.q>35.000000</cim:SvPowerFlow.q>
<cim:SvPowerFlow.Terminal rdf:resource="#E-1076922569"/>
</cim:SvPowerFlow>
<cim:SvPowerFlow rdf:ID="LOAD5-sv">
<cim:SvPowerFlow.p>125.000000</cim:SvPowerFlow.p>
<cim:SvPowerFlow.q>50.000000</cim:SvPowerFlow.q>
<cim:SvPowerFlow.Terminal rdf:resource="#E-1076922595"/>
</cim:SvPowerFlow>
<cim:SvPowerFlow rdf:ID="LOAD6-sv">
<cim:SvPowerFlow.p>90.000000</cim:SvPowerFlow.p>
<cim:SvPowerFlow.q>30.000000</cim:SvPowerFlow.q>
<cim:SvPowerFlow.Terminal rdf:resource="#E-1076922611"/>
</cim:SvPowerFlow>
<cim:SvPowerFlow rdf:ID="GEN3-sv">
<cim:SvPowerFlow.p>-85.000000</cim:SvPowerFlow.p>
<cim:SvPowerFlow.q>10.859709</cim:SvPowerFlow.q>
<cim:SvPowerFlow.Terminal rdf:resource="#E-1076922537"/>
</cim:SvPowerFlow>
<cim:SvPowerFlow rdf:ID="GEN1-sv">
<cim:SvPowerFlow.p>-71.641021</cim:SvPowerFlow.p>
<cim:SvPowerFlow.q>-27.045923</cim:SvPowerFlow.q>
<cim:SvPowerFlow.Terminal rdf:resource="#E-1076922656"/>
</cim:SvPowerFlow>
<cim:SvPowerFlow rdf:ID="GEN2-sv">
<cim:SvPowerFlow.p>-163.000000</cim:SvPowerFlow.p>
<cim:SvPowerFlow.q>-6.653660</cim:SvPowerFlow.q>
<cim:SvPowerFlow.Terminal rdf:resource="#E-1076922486"/>
</cim:SvPowerFlow>
</rdf:RDF>
147 changes: 147 additions & 0 deletions examples/sample_data/WSCC9/WSCC-09_Dyn_Fourth_TP.xml
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,147 @@
<?xml version="1.0" encoding="utf-8"?>
<rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#" xmlns:cim="http://iec.ch/TC57/2012/CIM-schema-cim16#" xmlns:md="http://iec.ch/TC57/61970-552/ModelDescription/1#" xmlns:entsoe="http://entsoe.eu/Secretariat/ProfileExtension/2#">
<cim:NameType rdf:ID="_1b855250-0d16-41e4-bf2b-725a23b25b04">
<cim:NameType.name>description</cim:NameType.name>
</cim:NameType>
<cim:TopologicalNode rdf:ID="BUS3">
<cim:IdentifiedObject.name>BUS3</cim:IdentifiedObject.name>
<cim:TopologicalNode.BaseVoltage rdf:resource="#_5efd264d-c4b8-41aa-921f-e423e58fee04"/>
<cim:TopologicalNode.ConnectivityNodeContainer rdf:resource="#_bf963bb0-be6f-44a9-bccd-7469ee94a295"/>
</cim:TopologicalNode>
<cim:TopologicalNode rdf:ID="BUS5">
<cim:IdentifiedObject.name>BUS5</cim:IdentifiedObject.name>
<cim:TopologicalNode.BaseVoltage rdf:resource="#_e8c76dfc-f327-4e44-b9ce-01f62639fcef"/>
<cim:TopologicalNode.ConnectivityNodeContainer rdf:resource="#_ad59543c-9983-487a-83a2-0dc24edfecad"/>
</cim:TopologicalNode>
<cim:TopologicalNode rdf:ID="BUS6">
<cim:IdentifiedObject.name>BUS6</cim:IdentifiedObject.name>
<cim:TopologicalNode.BaseVoltage rdf:resource="#_e8c76dfc-f327-4e44-b9ce-01f62639fcef"/>
<cim:TopologicalNode.ConnectivityNodeContainer rdf:resource="#_2d832774-e278-4b41-ba9e-4f9731998b2f"/>
</cim:TopologicalNode>
<cim:TopologicalNode rdf:ID="BUS4">
<cim:IdentifiedObject.name>BUS4</cim:IdentifiedObject.name>
<cim:TopologicalNode.BaseVoltage rdf:resource="#_e8c76dfc-f327-4e44-b9ce-01f62639fcef"/>
<cim:TopologicalNode.ConnectivityNodeContainer rdf:resource="#_68b2a0e0-1a2e-4e8c-af87-86c54dfb0768"/>
</cim:TopologicalNode>
<cim:TopologicalNode rdf:ID="BUS1">
<cim:IdentifiedObject.name>BUS1</cim:IdentifiedObject.name>
<cim:TopologicalNode.BaseVoltage rdf:resource="#_0cae189e-5250-4cd2-9ab7-8bb9f1a58c8b"/>
<cim:TopologicalNode.ConnectivityNodeContainer rdf:resource="#_747d2a10-a3a7-45cf-abfb-05b933fc78f0"/>
</cim:TopologicalNode>
<cim:TopologicalNode rdf:ID="BUS2">
<cim:IdentifiedObject.name>BUS2</cim:IdentifiedObject.name>
<cim:TopologicalNode.BaseVoltage rdf:resource="#_5c7701c6-e7cd-46f0-a27c-489c5e45afb3"/>
<cim:TopologicalNode.ConnectivityNodeContainer rdf:resource="#_c941c004-7df4-4955-a07e-b3ac863e5901"/>
</cim:TopologicalNode>
<cim:TopologicalNode rdf:ID="BUS7">
<cim:IdentifiedObject.name>BUS7</cim:IdentifiedObject.name>
<cim:TopologicalNode.BaseVoltage rdf:resource="#_e8c76dfc-f327-4e44-b9ce-01f62639fcef"/>
<cim:TopologicalNode.ConnectivityNodeContainer rdf:resource="#_b206adb1-acad-49db-bdb9-3d0a7bd0ea68"/>
</cim:TopologicalNode>
<cim:TopologicalNode rdf:ID="BUS9">
<cim:IdentifiedObject.name>BUS9</cim:IdentifiedObject.name>
<cim:TopologicalNode.BaseVoltage rdf:resource="#_e8c76dfc-f327-4e44-b9ce-01f62639fcef"/>
<cim:TopologicalNode.ConnectivityNodeContainer rdf:resource="#_776687ee-84ed-4b36-b3ad-ba62dcad851b"/>
</cim:TopologicalNode>
<cim:TopologicalNode rdf:ID="BUS8">
<cim:IdentifiedObject.name>BUS8</cim:IdentifiedObject.name>
<cim:TopologicalNode.BaseVoltage rdf:resource="#_e8c76dfc-f327-4e44-b9ce-01f62639fcef"/>
<cim:TopologicalNode.ConnectivityNodeContainer rdf:resource="#_ae2ca5db-bcd5-46d0-9280-7234751944da"/>
</cim:TopologicalNode>
<cim:Terminal rdf:about="#E-1076922586">
<cim:Terminal.connected>true</cim:Terminal.connected>
<cim:Terminal.TopologicalNode rdf:resource="#BUS7"/>
</cim:Terminal>
<cim:Terminal rdf:about="#E-1076922587">
<cim:Terminal.connected>true</cim:Terminal.connected>
<cim:Terminal.TopologicalNode rdf:resource="#BUS5"/>
</cim:Terminal>
<cim:Terminal rdf:about="#E-1076922602">
<cim:Terminal.connected>true</cim:Terminal.connected>
<cim:Terminal.TopologicalNode rdf:resource="#BUS9"/>
</cim:Terminal>
<cim:Terminal rdf:about="#E-1076922603">
<cim:Terminal.connected>true</cim:Terminal.connected>
<cim:Terminal.TopologicalNode rdf:resource="#BUS6"/>
</cim:Terminal>
<cim:Terminal rdf:about="#E-1076922623">
<cim:Terminal.connected>true</cim:Terminal.connected>
<cim:Terminal.TopologicalNode rdf:resource="#BUS6"/>
</cim:Terminal>
<cim:Terminal rdf:about="#E-1076922624">
<cim:Terminal.connected>true</cim:Terminal.connected>
<cim:Terminal.TopologicalNode rdf:resource="#BUS4"/>
</cim:Terminal>
<cim:Terminal rdf:about="#E-1076922618">
<cim:Terminal.connected>true</cim:Terminal.connected>
<cim:Terminal.TopologicalNode rdf:resource="#BUS5"/>
</cim:Terminal>
<cim:Terminal rdf:about="#E-1076922619">
<cim:Terminal.connected>true</cim:Terminal.connected>
<cim:Terminal.TopologicalNode rdf:resource="#BUS4"/>
</cim:Terminal>
<cim:Terminal rdf:about="#E-1076922524">
<cim:Terminal.connected>true</cim:Terminal.connected>
<cim:Terminal.TopologicalNode rdf:resource="#BUS9"/>
</cim:Terminal>
<cim:Terminal rdf:about="#E-1076922525">
<cim:Terminal.connected>true</cim:Terminal.connected>
<cim:Terminal.TopologicalNode rdf:resource="#BUS8"/>
</cim:Terminal>
<cim:Terminal rdf:about="#E-1076922516">
<cim:Terminal.connected>true</cim:Terminal.connected>
<cim:Terminal.TopologicalNode rdf:resource="#BUS7"/>
</cim:Terminal>
<cim:Terminal rdf:about="#E-1076922517">
<cim:Terminal.connected>true</cim:Terminal.connected>
<cim:Terminal.TopologicalNode rdf:resource="#BUS8"/>
</cim:Terminal>
<cim:Terminal rdf:about="#E-1076922554">
<cim:Terminal.connected>true</cim:Terminal.connected>
<cim:Terminal.TopologicalNode rdf:resource="#BUS3"/>
</cim:Terminal>
<cim:Terminal rdf:about="#E-1076922557">
<cim:Terminal.connected>true</cim:Terminal.connected>
<cim:Terminal.TopologicalNode rdf:resource="#BUS9"/>
</cim:Terminal>
<cim:Terminal rdf:about="#E-1076922644">
<cim:Terminal.connected>true</cim:Terminal.connected>
<cim:Terminal.TopologicalNode rdf:resource="#BUS1"/>
</cim:Terminal>
<cim:Terminal rdf:about="#E-1076922647">
<cim:Terminal.connected>true</cim:Terminal.connected>
<cim:Terminal.TopologicalNode rdf:resource="#BUS4"/>
</cim:Terminal>
<cim:Terminal rdf:about="#E-1076922506">
<cim:Terminal.connected>true</cim:Terminal.connected>
<cim:Terminal.TopologicalNode rdf:resource="#BUS2"/>
</cim:Terminal>
<cim:Terminal rdf:about="#E-1076922509">
<cim:Terminal.connected>true</cim:Terminal.connected>
<cim:Terminal.TopologicalNode rdf:resource="#BUS7"/>
</cim:Terminal>
<cim:Terminal rdf:about="#E-1076922537">
<cim:Terminal.connected>true</cim:Terminal.connected>
<cim:Terminal.TopologicalNode rdf:resource="#BUS3"/>
</cim:Terminal>
<cim:Terminal rdf:about="#E-1076922656">
<cim:Terminal.connected>true</cim:Terminal.connected>
<cim:Terminal.TopologicalNode rdf:resource="#BUS1"/>
</cim:Terminal>
<cim:Terminal rdf:about="#E-1076922486">
<cim:Terminal.connected>true</cim:Terminal.connected>
<cim:Terminal.TopologicalNode rdf:resource="#BUS2"/>
</cim:Terminal>
<cim:Terminal rdf:about="#E-1076922569">
<cim:Terminal.connected>true</cim:Terminal.connected>
<cim:Terminal.TopologicalNode rdf:resource="#BUS8"/>
</cim:Terminal>
<cim:Terminal rdf:about="#E-1076922595">
<cim:Terminal.connected>true</cim:Terminal.connected>
<cim:Terminal.TopologicalNode rdf:resource="#BUS5"/>
</cim:Terminal>
<cim:Terminal rdf:about="#E-1076922611">
<cim:Terminal.connected>true</cim:Terminal.connected>
<cim:Terminal.TopologicalNode rdf:resource="#BUS6"/>
</cim:Terminal>
</rdf:RDF>
Loading