[go: up one dir, main page]
Scribd
Upload
9 views4 pages

Time Series Ex

Uploaded by

hod
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as DOCX, PDF, TXT or read online on Scribd
9 views4 pages

Time Series Ex

Uploaded by

hod
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as DOCX, PDF, TXT or read online on Scribd
You are on page 1/ 4

import os

import numpy as np

import pandas as pd

import tempfile

import random

import pandapower as pp

from pandapower.timeseries import DFData

from pandapower.timeseries import OutputWriter

from pandapower.timeseries.run_time_series import run_timeseries

from pandapower.control import ConstControl

random.seed(10)

##########################################a simple example pandapower net consisting of five


buses, a transformer, three lines, a load and a sgen.

# Create a simple power grid

def simple_test_net():

"""

# simple net that looks like:

ext_grid b0---b1 trafo(110/20) b2----b3 load

b4 sgen

"""

net = pp.create_empty_network()

pp.set_user_pf_options(net, init_vm_pu = "flat", init_va_degree = "dc",


calculate_voltage_angles=True)
b0 = pp.create_bus(net, 110)

b1 = pp.create_bus(net, 110)

b2 = pp.create_bus(net, 20)

b3 = pp.create_bus(net, 20)

b4 = pp.create_bus(net, 20)

pp.create_ext_grid(net, b0)

pp.create_line(net, b0, b1, 10, "149-AL1/24-ST1A 110.0")

pp.create_transformer(net, b1, b2, "25 MVA 110/20 kV", name='tr1')

pp.create_line(net, b2, b3, 10, "184-AL1/30-ST1A 20.0")

pp.create_line(net, b2, b4, 10, "184-AL1/30-ST1A 20.0")

pp.create_load(net, b3, p_mw=15., q_mvar=10., name='load1')

pp.create_sgen(net, b4, p_mw=20., q_mvar=0.15, name='sgen1')

return net

net = simple_test_net()

#print(net)

#############################

def timeseries_example(output_dir):

# 1. create test net

net = simple_test_net()

# 2. create (random) data source

n_timesteps = 10

profiles, ds = create_data_source(n_timesteps)
# 3. create controllers (to control P values of the load and the sgen)

create_controllers(net, ds)

# time steps to be calculated. Could also be a list with non-consecutive time steps

time_steps = range(0, n_timesteps)

# 4. the output writer with the desired results to be stored to files.

ow = create_output_writer(net, time_steps, output_dir=output_dir)

# 5. the main time series function

run_timeseries(net, time_steps)

print(net.res_line.loading_percent)

#############The DFData(profiles) converts the Dataframe to the required format for the controllers.

def create_data_source(n_timesteps=24):

profiles = pd.DataFrame()

profiles['load1_p'] = np.random.random(n_timesteps) * 15.

profiles['sgen1_p'] = np.random.random(n_timesteps) * 20.

ds = DFData(profiles)

return profiles, ds

#print(profiles.head())

########create the controllers by telling the function which element_index belongs to which profile.

def create_controllers(net, ds):

ConstControl(net, element='load', variable='p_mw', element_index=[0],

data_source=ds, profile_name=["load1_p"])

ConstControl(net, element='sgen', variable='p_mw', element_index=[0],

data_source=ds, profile_name=["sgen1_p"])

####### create the output writer. Instead of saving the whole net (which takes a lot of time), we extract
only pre defined outputs

def create_output_writer(net, time_steps, output_dir):


ow = OutputWriter(net, time_steps, output_path=output_dir, output_file_type=".xlsx",
log_variables=list())

# these variables are saved to the harddisk after / during the time series loop

ow.log_variable('res_load', 'p_mw')

ow.log_variable('res_bus', 'vm_pu')

ow.log_variable('res_line', 'loading_percent')

ow.log_variable('res_line', 'i_ka')

return ow

##### Now lets execute the code.

output_dir = os.path.join(tempfile.gettempdir(), "time_series_example")

print("Results can be found in your local temp folder: {}".format(output_dir))

if not os.path.exists(output_dir):

os.mkdir(output_dir)

timeseries_example(output_dir)

You might also like