ltbams 0.9.9__py3-none-any.whl → 1.0.2__py3-none-any.whl

ams/routines/rted.py

@@ -0,0 +1,519 @@
+"""
+Real-time economic dispatch.
+"""
+import logging
+from collections import OrderedDict
+import numpy as np
+
+from ams.core.param import RParam
+from ams.core.service import ZonalSum, VarSelect, NumOp, NumOpDual
+from ams.routines.dcopf import DCOPF
+
+from ams.opt import Var, Constraint
+
+logger = logging.getLogger(__name__)
+
+
+class RTEDBase:
+    """
+    Base class for real-time economic dispatch (RTED).
+    """
+
+    def __init__(self):
+        # --- zone ---
+        self.zg = RParam(info='Gen zone',
+                         name='zg', tex_name='z_{one,g}',
+                         model='StaticGen', src='zone',
+                         no_parse=True)
+        self.zd = RParam(info='Load zone',
+                         name='zd', tex_name='z_{one,d}',
+                         model='StaticLoad', src='zone',
+                         no_parse=True)
+        self.gs = ZonalSum(u=self.zg, zone='Zone',
+                           name='gs', tex_name=r'S_{g}',
+                           info='Sum Gen vars vector in shape of zone',
+                           no_parse=True, sparse=True)
+        self.ds = ZonalSum(u=self.zd, zone='Zone',
+                           name='ds', tex_name=r'S_{d}',
+                           info='Sum pd vector in shape of zone',
+                           no_parse=True,)
+        self.pdz = NumOpDual(u=self.ds, u2=self.pd,
+                             fun=np.multiply,
+                             rfun=np.sum, rargs=dict(axis=1),
+                             expand_dims=0,
+                             name='pdz', tex_name=r'p_{d,z}',
+                             unit='p.u.', info='zonal total load',
+                             no_parse=True,)
+        # --- generator ---
+        self.R10 = RParam(info='10-min ramp rate',
+                          name='R10', tex_name=r'R_{10}',
+                          model='StaticGen', src='R10',
+                          unit='p.u./h',)
+
+
+class SFRBase:
+    """
+    Base class for SFR used in DCED.
+    """
+
+    def __init__(self):
+        #  --- SFR cost ---
+        self.cru = RParam(info='RegUp reserve coefficient',
+                          name='cru', tex_name=r'c_{r,u}',
+                          model='SFRCost', src='cru',
+                          indexer='gen', imodel='StaticGen',
+                          unit=r'$/(p.u.)',)
+        self.crd = RParam(info='RegDown reserve coefficient',
+                          name='crd', tex_name=r'c_{r,d}',
+                          model='SFRCost', src='crd',
+                          indexer='gen', imodel='StaticGen',
+                          unit=r'$/(p.u.)',)
+        # --- reserve requirement ---
+        self.du = RParam(info='RegUp reserve requirement in percentage',
+                         name='du', tex_name=r'd_{u}',
+                         model='SFR', src='du',
+                         unit='%', no_parse=True,)
+        self.dd = RParam(info='RegDown reserve requirement in percentage',
+                         name='dd', tex_name=r'd_{d}',
+                         model='SFR', src='dd',
+                         unit='%', no_parse=True,)
+        self.dud = NumOpDual(u=self.pdz, u2=self.du, fun=np.multiply,
+                             rfun=np.reshape, rargs=dict(newshape=(-1,)),
+                             name='dud', tex_name=r'd_{u, d}',
+                             info='zonal RegUp reserve requirement',)
+        self.ddd = NumOpDual(u=self.pdz, u2=self.dd, fun=np.multiply,
+                             rfun=np.reshape, rargs=dict(newshape=(-1,)),
+                             name='ddd', tex_name=r'd_{d, d}',
+                             info='zonal RegDn reserve requirement',)
+        # --- SFR ---
+        self.pru = Var(info='RegUp reserve',
+                       unit='p.u.', name='pru', tex_name=r'p_{r,u}',
+                       model='StaticGen', nonneg=True,)
+        self.prd = Var(info='RegDn reserve',
+                       unit='p.u.', name='prd', tex_name=r'p_{r,d}',
+                       model='StaticGen', nonneg=True,)
+        # NOTE: define e_str in scheduling routine
+        self.rbu = Constraint(name='rbu', is_eq=True,
+                              info='RegUp reserve balance',)
+        self.rbd = Constraint(name='rbd', is_eq=True,
+                              info='RegDn reserve balance',)
+        self.rru = Constraint(name='rru', is_eq=False,
+                              info='RegUp reserve source',)
+        self.rrd = Constraint(name='rrd', is_eq=False,
+                              info='RegDn reserve source',)
+        self.rgu = Constraint(name='rgu', is_eq=False,
+                              info='Gen ramping up',)
+        self.rgd = Constraint(name='rgd', is_eq=False,
+                              info='Gen ramping down',)
+
+
+class RTED(DCOPF, RTEDBase, SFRBase):
+    """
+    DC-based real-time economic dispatch (RTED).
+    RTED extends DCOPF with:
+
+    - Mapping dicts to interface with ANDES
+    - Function ``dc2ac`` to do the AC conversion
+    - Vars for SFR reserve: ``pru`` and ``prd``
+    - Param for linear SFR cost: ``cru`` and ``crd``
+    - Param for SFR requirement: ``du`` and ``dd``
+    - Param for ramping: start point ``pg0`` and ramping limit ``R10``
+    - Param ``pg0``, which can be retrieved from dynamic simulation results.
+
+    The function ``dc2ac`` sets the ``vBus`` value from solved ACOPF.
+    Without this conversion, dynamic simulation might fail due to the gap between
+    DC-based dispatch results and AC-based dynamic initialization.
+
+    Notes
+    -----
+    1. Formulations has been adjusted with interval ``config.t``, 5/60 [Hour] by default.
+
+    2. The tie-line flow has not been implemented in formulations.
+    """
+
+    def __init__(self, system, config):
+        DCOPF.__init__(self, system, config)
+        RTEDBase.__init__(self)
+        SFRBase.__init__(self)
+
+        self.config.add(OrderedDict((('t', 5/60),
+                                     )))
+        self.config.add_extra("_help",
+                              t="time interval in hours",
+                              )
+        self.config.add_extra("_tex",
+                              t='T_{cfg}',
+                              )
+
+        self.info = 'Real-time economic dispatch'
+        self.type = 'DCED'
+
+        # --- Mapping Section ---
+        # Add p -> pg0 in from map
+        self.map1.update({
+            'pg0': ('StaticGen', 'p'),
+        })
+        # nothing to do with to map
+
+        # --- Model Section ---
+        # --- SFR ---
+        # RegUp/Dn reserve balance
+        self.rbu.e_str = 'gs @ mul(ug, pru) - dud'
+        self.rbd.e_str = 'gs @ mul(ug, prd) - ddd'
+        # RegUp/Dn reserve source
+        self.rru.e_str = 'mul(ug, (pg + pru)) - mul(ug, pmaxe)'
+        self.rrd.e_str = 'mul(ug, (-pg + prd)) + mul(ug, pmine)'
+        # Gen ramping up/down
+        self.rgu.e_str = 'mul(ug, (pg-pg0-R10))'
+        self.rgd.e_str = 'mul(ug, (-pg+pg0-R10))'
+
+        # --- objective ---
+        self.obj.info = 'total generation and reserve cost'
+        # NOTE: the product involved t should use ``dot``
+        cost = 't**2 dot sum(mul(c2, pg**2)) + sum(ug * c0)'
+        _to_sum = 'c1 @ pg + cru * pru + crd * prd'
+        cost += f'+ t dot sum({_to_sum})'
+        self.obj.e_str = cost
+
+    def dc2ac(self, kloss=1.0, **kwargs):
+        """
+        Convert the RTED results with ACOPF.
+
+        Parameters
+        ----------
+        kloss : float, optional
+            The loss factor for the conversion. Defaults to 1.2.
+        """
+        exec_time = self.exec_time
+        if self.exec_time == 0 or self.exit_code != 0:
+            logger.warning(f'{self.class_name} is not executed successfully, quit conversion.')
+            return False
+        # set pru and prd into pmin and pmax
+        pr_idx = self.pru.get_all_idxes()
+        pmin0 = self.system.StaticGen.get(src='pmin', attr='v', idx=pr_idx)
+        pmax0 = self.system.StaticGen.get(src='pmax', attr='v', idx=pr_idx)
+        p00 = self.system.StaticGen.get(src='p0', attr='v', idx=pr_idx)
+
+        # --- ACOPF ---
+        # scale up load
+        pq_idx = self.system.StaticLoad.get_all_idxes()
+        pd0 = self.system.StaticLoad.get(src='p0', attr='v', idx=pq_idx).copy()
+        qd0 = self.system.StaticLoad.get(src='q0', attr='v', idx=pq_idx).copy()
+        self.system.StaticLoad.set(src='p0', idx=pq_idx, attr='v', value=pd0 * kloss)
+        self.system.StaticLoad.set(src='q0', idx=pq_idx, attr='v', value=qd0 * kloss)
+        # preserve generator reserve
+        ACOPF = self.system.ACOPF
+        pmin = pmin0 + self.prd.v
+        pmax = pmax0 - self.pru.v
+        self.system.StaticGen.set(src='pmin', idx=pr_idx, attr='v', value=pmin)
+        self.system.StaticGen.set(src='pmax', idx=pr_idx, attr='v', value=pmax)
+        self.system.StaticGen.set(src='p0', idx=pr_idx, attr='v', value=self.pg.v)
+        # run ACOPF
+        ACOPF.run()
+        # scale load back
+        self.system.StaticLoad.set(src='p0', idx=pq_idx, attr='v', value=pd0)
+        self.system.StaticLoad.set(src='q0', idx=pq_idx, attr='v', value=qd0)
+        if not ACOPF.exit_code == 0:
+            logger.warning('<ACOPF> did not converge, conversion failed.')
+            self.vBus.optz.value = np.ones(self.system.Bus.n)
+            self.aBus.optz.value = np.zeros(self.system.Bus.n)
+            return False
+
+        self.pg.optz.value = ACOPF.pg.v
+        self.vBus.optz.value = ACOPF.vBus.v
+        self.aBus.optz.value = ACOPF.aBus.v
+        self.exec_time = exec_time
+
+        # reset pmin, pmax, p0
+        self.system.StaticGen.set(src='pmin', idx=pr_idx, attr='v', value=pmin0)
+        self.system.StaticGen.set(src='pmax', idx=pr_idx, attr='v', value=pmax0)
+        self.system.StaticGen.set(src='p0', idx=pr_idx, attr='v', value=p00)
+
+        # --- set status ---
+        self.system.recent = self
+        self.converted = True
+        logger.warning(f'<{self.class_name}> converted to AC.')
+        return True
+
+
+class DGBase:
+    """
+    Base class for DG used in DCED.
+    """
+
+    def __init__(self):
+        # --- params ---
+        self.gendg = RParam(info='gen of DG',
+                            name='gendg', tex_name=r'g_{DG}',
+                            model='DG', src='gen',
+                            no_parse=True,)
+        info = 'Ratio of DG.pge w.r.t to that of static generator'
+        self.gammapdg = RParam(name='gammapd', tex_name=r'\gamma_{p,DG}',
+                               model='DG', src='gammap',
+                               no_parse=True, info=info)
+
+        # --- vars ---
+        # TODO: maybe there will be constraints on pgd, maybe upper/lower bound?
+        # TODO: this might requre new device like DGSlot
+        self.pgdg = Var(info='DG output power',
+                        unit='p.u.', name='pgdg',
+                        tex_name=r'p_{g,DG}',
+                        model='DG',)
+
+        # --- constraints ---
+        self.cdg = VarSelect(u=self.pg, indexer='gendg',
+                             name='cd', tex_name=r'C_{DG}',
+                             info='Select DG power from pg',
+                             gamma='gammapdg',
+                             no_parse=True, sparse=True,)
+        self.cdgb = Constraint(name='cdgb', is_eq=True,
+                               info='Select DG power from pg',
+                               e_str='cdg @ pg - pgdg',)
+
+
+class RTEDDG(RTED, DGBase):
+    """
+    RTED with distributed generator :ref:`DG`.
+
+    Note that RTEDDG only inlcudes DG output power. If ESD1 is included,
+    RTEDES should be used instead, otherwise there is no SOC.
+    """
+
+    def __init__(self, system, config):
+        RTED.__init__(self, system, config)
+        DGBase.__init__(self)
+        self.info = 'Real-time economic dispatch with DG'
+        self.type = 'DCED'
+
+
+class ESD1Base(DGBase):
+    """
+    Base class for ESD1 used in DCED.
+    """
+
+    def __init__(self):
+        DGBase.__init__(self)
+
+        # --- params ---
+        self.En = RParam(info='Rated energy capacity',
+                         name='En', src='En',
+                         tex_name='E_n', unit='MWh',
+                         model='ESD1', no_parse=True,)
+        self.SOCmax = RParam(info='Maximum allowed value for SOC in limiter',
+                             name='SOCmax', src='SOCmax',
+                             tex_name=r'SOC_{max}', unit='%',
+                             model='ESD1',)
+        self.SOCmin = RParam(info='Minimum required value for SOC in limiter',
+                             name='SOCmin', src='SOCmin',
+                             tex_name=r'SOC_{min}', unit='%',
+                             model='ESD1',)
+        self.SOCinit = RParam(info='Initial SOC',
+                              name='SOCinit', src='SOCinit',
+                              tex_name=r'SOC_{init}', unit='%',
+                              model='ESD1',)
+        self.EtaC = RParam(info='Efficiency during charging',
+                           name='EtaC', src='EtaC',
+                           tex_name=r'\eta_c', unit='%',
+                           model='ESD1', no_parse=True,)
+        self.EtaD = RParam(info='Efficiency during discharging',
+                           name='EtaD', src='EtaD',
+                           tex_name=r'\eta_d', unit='%',
+                           model='ESD1', no_parse=True,)
+
+        # --- service ---
+        self.REtaD = NumOp(name='REtaD', tex_name=r'\frac{1}{\eta_d}',
+                           u=self.EtaD, fun=np.reciprocal,)
+        self.Mb = NumOp(info='10 times of max of pmax as big M',
+                        name='Mb', tex_name=r'M_{big}',
+                        u=self.pmax, fun=np.max,
+                        rfun=np.dot, rargs=dict(b=10),
+                        array_out=False,)
+
+        # --- vars ---
+        self.SOC = Var(info='ESD1 State of Charge', unit='%',
+                       name='SOC', tex_name=r'SOC',
+                       model='ESD1', pos=True,
+                       v0=self.SOCinit,)
+        self.SOClb = Constraint(name='SOClb', is_eq=False,
+                                info='SOC lower bound',
+                                e_str='-SOC + SOCmin',)
+        self.SOCub = Constraint(name='SOCub', is_eq=False,
+                                info='SOC upper bound',
+                                e_str='SOC - SOCmax',)
+        self.pce = Var(info='ESD1 charging power',
+                       unit='p.u.', name='pce',
+                       tex_name=r'p_{c,ESD}',
+                       model='ESD1', nonneg=True,)
+        self.pde = Var(info='ESD1 discharging power',
+                       unit='p.u.', name='pde',
+                       tex_name=r'p_{d,ESD}',
+                       model='ESD1', nonneg=True,)
+        self.uce = Var(info='ESD1 charging decision',
+                       name='uce', tex_name=r'u_{c,ESD}',
+                       model='ESD1', boolean=True,)
+        self.ude = Var(info='ESD1 discharging decision',
+                       name='ude', tex_name=r'u_{d,ESD}',
+                       model='ESD1', boolean=True,)
+        self.zce = Var(name='zce', tex_name=r'z_{c,ESD}',
+                       model='ESD1', nonneg=True,)
+        self.zce.info = 'Aux var for charging, '
+        self.zce.info += ':math:`z_{c,ESD}=u_{c,ESD}*p_{c,ESD}`'
+        self.zde = Var(name='zde', tex_name=r'z_{d,ESD}',
+                       model='ESD1', nonneg=True,)
+        self.zde.info = 'Aux var for discharging, '
+        self.zde.info += ':math:`z_{d,ESD}=u_{d,ESD}*p_{d,ESD}`'
+
+        # NOTE: to ensure consistency with DG based routiens,
+        # here we select ESD1 power from DG rather than StaticGen
+        self.genesd = RParam(info='gen of ESD1',
+                             name='genesd', tex_name=r'g_{ESD}',
+                             model='ESD1', src='idx',
+                             no_parse=True,)
+        self.ces = VarSelect(u=self.pgdg, indexer='genesd',
+                             name='ces', tex_name=r'C_{ESD}',
+                             info='Select ESD power from DG',
+                             no_parse=True)
+        self.cescb = Constraint(name='cescb', is_eq=True,
+                                info='Select pce from DG',
+                                e_str='ces @ pgdg - pce',)
+        self.cesdb = Constraint(name='cesdb', is_eq=True,
+                                info='Select pde from DG',
+                                e_str='ces @ pgdg - pde',)
+
+        # --- constraints ---
+        self.cdb = Constraint(name='cdb', is_eq=True,
+                              info='Charging decision bound',
+                              e_str='uce + ude - 1',)
+
+        self.zce1 = Constraint(name='zce1', is_eq=False, info='zce bound 1',
+                               e_str='-zce + pce',)
+        self.zce2 = Constraint(name='zce2', is_eq=False, info='zce bound 2',
+                               e_str='zce - pce - Mb dot (1-uce)',)
+        self.zce3 = Constraint(name='zce3', is_eq=False, info='zce bound 3',
+                               e_str='zce - Mb dot uce',)
+
+        self.zde1 = Constraint(name='zde1', is_eq=False, info='zde bound 1',
+                               e_str='-zde + pde',)
+        self.zde2 = Constraint(name='zde2', is_eq=False, info='zde bound 2',
+                               e_str='zde - pde - Mb dot (1-ude)',)
+        self.zde3 = Constraint(name='zde3', is_eq=False, info='zde bound 3',
+                               e_str='zde - Mb dot ude',)
+
+        SOCb = 'mul(En, (SOC - SOCinit)) - t dot mul(EtaC, zce)'
+        SOCb += '+ t dot mul(REtaD, zde)'
+        self.SOCb = Constraint(name='SOCb', is_eq=True,
+                               info='ESD1 SOC balance',
+                               e_str=SOCb,)
+
+
+class RTEDES(RTED, ESD1Base):
+    """
+    RTED with energy storage :ref:`ESD1`.
+    The bilinear term in the formulation is linearized with big-M method.
+    """
+
+    def __init__(self, system, config):
+        RTED.__init__(self, system, config)
+        ESD1Base.__init__(self)
+        self.info = 'Real-time economic dispatch with energy storage'
+        self.type = 'DCED'
+
+
+class VISBase:
+    """
+    Base class for virtual inertia scheduling.
+    """
+
+    def __init__(self) -> None:
+        # --- Data Section ---
+        self.cm = RParam(info='Virtual inertia cost',
+                         name='cm', src='cm',
+                         tex_name=r'c_{m}', unit=r'$/s',
+                         model='VSGCost',
+                         indexer='reg', imodel='VSG')
+        self.cd = RParam(info='Virtual damping cost',
+                         name='cd', src='cd',
+                         tex_name=r'c_{d}', unit=r'$/(p.u.)',
+                         model='VSGCost',
+                         indexer='reg', imodel='VSG',)
+        self.zvsg = RParam(info='VSG zone',
+                           name='zvsg', tex_name='z_{one,vsg}',
+                           model='VSG', src='zone',
+                           no_parse=True)
+        self.Mmax = RParam(info='Maximum inertia emulation',
+                           name='Mmax', tex_name='M_{max}',
+                           model='VSG', src='Mmax',
+                           unit='s',)
+        self.Dmax = RParam(info='Maximum damping emulation',
+                           name='Dmax', tex_name='D_{max}',
+                           model='VSG', src='Dmax',
+                           unit='p.u.',)
+        self.dvm = RParam(info='Emulated inertia requirement',
+                          name='dvm', tex_name=r'd_{v,m}',
+                          unit='s',
+                          model='VSGR', src='dvm',)
+        self.dvd = RParam(info='Emulated damping requirement',
+                          name='dvd', tex_name=r'd_{v,d}',
+                          unit='p.u.',
+                          model='VSGR', src='dvd',)
+
+        # --- Model Section ---
+        self.M = Var(info='Emulated startup time constant (M=2H)',
+                     name='M', tex_name=r'M', unit='s',
+                     model='VSG', src='M',
+                     nonneg=True,)
+        self.D = Var(info='Emulated damping coefficient',
+                     name='D', tex_name=r'D', unit='p.u.',
+                     model='VSG', src='D',
+                     nonneg=True,)
+
+        self.gvsg = ZonalSum(u=self.zvsg, zone='Zone',
+                             name='gvsg', tex_name=r'S_{g}',
+                             info='Sum VSG vars vector in shape of zone',
+                             no_parse=True)
+        self.Mub = Constraint(name='Mub', is_eq=False,
+                              info='M upper bound',
+                              e_str='M - Mmax',)
+        self.Dub = Constraint(name='Dub', is_eq=False,
+                              info='D upper bound',
+                              e_str='D - Dmax',)
+        self.Mreq = Constraint(name='Mreq', is_eq=True,
+                               info='Emulated inertia requirement',
+                               e_str='-gvsg@M + dvm',)
+        self.Dreq = Constraint(name='Dreq', is_eq=True,
+                               info='Emulated damping requirement',
+                               e_str='-gvsg@D + dvd',)
+
+        # NOTE: revise the objective function to include virtual inertia cost
+
+
+class RTEDVIS(RTED, VISBase):
+    """
+    RTED with virtual inertia scheduling.
+
+    This class implements real-time economic dispatch with virtual inertia scheduling.
+    Please ensure that the parameters `dvm` and `dvd` are set according to the system base.
+
+    Reference:
+
+    [1] B. She, F. Li, H. Cui, J. Wang, Q. Zhang and R. Bo, "Virtual
+    Inertia Scheduling (VIS) for Real-time Economic Dispatch of
+    IBRs-penetrated Power Systems," in IEEE Transactions on
+    Sustainable Energy, doi: 10.1109/TSTE.2023.3319307.
+    """
+
+    def __init__(self, system, config):
+        RTED.__init__(self, system, config)
+        VISBase.__init__(self)
+        self.info = 'Real-time economic dispatch with virtual inertia scheduling'
+        self.type = 'DCED'
+
+        # --- objective ---
+        self.obj.info = 'total generation and reserve cost'
+        vsgcost = '+ t dot sum(cm * M + cd * D)'
+        self.obj.e_str += vsgcost
+
+        self.map2.update({
+            'M': ('RenGen', 'M'),
+            'D': ('RenGen', 'D'),
+        })

ams/routines/type.py

@@ -0,0 +1,160 @@
+import logging
+import inspect
+from collections import OrderedDict
+
+logger = logging.getLogger(__name__)
+
+
+class TypeBase:
+    """
+    Base class for types.
+    """
+
+    def __init__(self):
+
+        self.common_rparams = []
+        self.common_vars = []
+        self.common_constrs = []
+
+        self.routines = OrderedDict()
+
+    @property
+    def class_name(self):
+        return self.__class__.__name__
+
+    @property
+    def n(self):
+        """
+        Total number of routines.
+        """
+        return len(self.routines)
+
+    def __repr__(self):
+        dev_text = 'routine' if self.n == 1 else 'routines'
+        return f'{self.class_name} ({self.n} {dev_text}) at {hex(id(self))}'
+
+    def doc(self, export='plain'):
+        """
+        Return the documentation of the type in a string.
+        """
+        out = ''
+        if export == 'rest':
+            out += f'.. _{self.class_name}:\n\n'
+            group_header = '=' * 80 + '\n'
+        else:
+            group_header = ''
+
+        if export == 'rest':
+            out += group_header + f'{self.class_name}\n' + group_header
+        else:
+            out += group_header + f'Type <{self.class_name}>\n' + group_header
+
+        if self.__doc__ is not None:
+            out += inspect.cleandoc(self.__doc__) + '\n\n'
+
+        if len(self.common_rparams):
+            out += 'Common Parameters: ' + ', '.join(self.common_rparams)
+            out += '\n\n'
+        if len(self.common_vars):
+            out += 'Common Vars: ' + ', '.join(self.common_vars)
+            out += '\n\n'
+        if len(self.common_constrs):
+            out += 'Common Constraints: ' + ', '.join(self.common_constrs)
+            out += '\n\n'
+        if len(self.routines):
+            out += 'Available routines:\n'
+            rtn_name_list = list(self.routines.keys())
+
+            if export == 'rest':
+                def add_reference(name_list):
+                    return [f'{item}_' for item in name_list]
+
+                rtn_name_list = add_reference(rtn_name_list)
+
+            out += ',\n'.join(rtn_name_list) + '\n'
+
+        return out
+
+    def doc_all(self, export='plain'):
+        """
+        Return documentation of the type and its routines.
+
+        Parameters
+        ----------
+        export : 'plain' or 'rest'
+            Export format, plain-text or RestructuredText
+
+        Returns
+        -------
+        str
+
+        """
+        out = self.doc(export=export)
+        out += '\n'
+        for instance in self.routines.values():
+            out += instance.doc(export=export)
+            out += '\n'
+        return out
+
+
+class UndefinedType(TypeBase):
+    """
+    The undefined type.
+    """
+
+    def __init__(self):
+        TypeBase.__init__(self)
+
+
+class PF(TypeBase):
+    """
+    Type for power flow routines.
+    """
+
+    def __init__(self):
+        TypeBase.__init__(self)
+        self.common_rparams.extend(('pd',))
+        self.common_vars.extend(('pg',))
+
+
+class DCED(TypeBase):
+    """
+    Type for DC-based economic dispatch.
+    """
+
+    def __init__(self):
+        TypeBase.__init__(self)
+        self.common_rparams.extend(('c2', 'c1', 'c0', 'pmax', 'pmin', 'pd', 'ptdf', 'rate_a',))
+        self.common_vars.extend(('pg',))
+        self.common_constrs.extend(('pb', 'lub', 'llb'))
+
+
+class DCUC(TypeBase):
+    """
+    Type for DC-based unit commitment.
+    """
+
+    def __init__(self):
+        TypeBase.__init__(self)
+        # TODO: add common parameters and variables
+
+
+class DED(TypeBase):
+    """
+    Type for Distributional economic dispatch.
+    """
+
+    def __init__(self):
+        TypeBase.__init__(self)
+        # TODO: add common parameters and variables
+
+
+class ACED(DCED):
+    """
+    Type for AC-based economic dispatch.
+    """
+
+    def __init__(self):
+        DCED.__init__(self)
+        self.common_rparams.extend(('qd',))
+        self.common_vars.extend(('aBus', 'vBus', 'qg',))