<?xml version="1.0" encoding="utf-8"?>
<!DOCTYPE article
  PUBLIC "-//NLM//DTD JATS (Z39.96) Journal Publishing DTD v1.0 20120330//EN" "http://jats.nlm.nih.gov/publishing/1.0/JATS-journalpublishing1.dtd">
<article article-type="research-article" dtd-version="1.0" specific-use="sps-1.8" xml:lang="es"
    xmlns:mml="http://www.w3.org/1998/Math/MathML"
    xmlns:xlink="http://www.w3.org/1999/xlink">
    <front>
        <journal-meta>
            <journal-id journal-id-type="publisher-id">cein</journal-id>
            <journal-title-group>
                <journal-title>Ciencia e Ingeniería Neogranadina</journal-title>
                <abbrev-journal-title abbrev-type="publisher">Cienc. Ing. Neogranad.</abbrev-journal-title>
            </journal-title-group>
            <issn pub-type="epub">1909-7735</issn>
            <issn pub-type="ppub">0124-8170</issn>
            <publisher>
                <publisher-name>Universidad Militar Nueva Granada</publisher-name>
            </publisher>
        </journal-meta>
        <article-meta>
            <article-id pub-id-type="publisher-id">3635</article-id>
            <article-id pub-id-type="doi">10.18359/rcin.3635</article-id>
            <article-categories>
                <subj-group subj-group-type="heading">
                    <subject>Research Article</subject>
                </subj-group>
            </article-categories>
            <title-group>
                <article-title>Use of energy storage systems in the optimal operation of distribution networks</article-title>
            </title-group>
            <contrib-group>
                <contrib contrib-type="author">
                    <name>
                        <surname>Valencia D&#x00ED;aza</surname>
                        <given-names>Alejandro</given-names>
                    </name>
                    <xref ref-type="aff" rid="aff1"><sup>a</sup></xref>
                </contrib>
                    <contrib contrib-type="author">
                        <name>
                            <surname>Gallego Rend&#x00F3;n</surname>
                            <given-names>Ram&#x00F3;n Alfonso</given-names>
                        </name>
                    <xref ref-type="aff" rid="aff2"><sup>b</sup></xref>
                    </contrib>
                    <contrib contrib-type="author">
                        <name>
                            <surname>Hincapi&#x00E9; Isaza</surname>
                            <given-names>Ricardo Alberto</given-names>
                        </name>
                    <xref ref-type="aff" rid="aff3"><sup>c</sup></xref>
                    </contrib>
                </contrib-group>
                <aff id="aff1">
                    <label>a</label>
                    <institution content-type="original">Universidad Tecnológica de Pereira, Colombia. Correo electrónico: alejovd4512@utp.edu.co.</institution>
                    <institution content-type="normalized">universidad Tecnológica de Pereira</institution>
                    <institution content-type="orgname">universidad Tecnológica de Pereira</institution>
                    <country country="CO">Colombia</country>
                    <email>alejovd4512@utp.edu.co</email>
                </aff>
                <aff id="aff2">
                    <label>b</label>
                    <institution content-type="original">Universidad Tecnológica de Pereira, Colombia. Correo electrónico: ragr@utp.edu.co</institution>
                    <institution content-type="normalized">universidad Tecnológica de Pereira</institution>
                    <institution content-type="orgname">universidad Tecnológica de Pereira</institution>
                    <country country="CO">Colombia</country>
                    <email>ragr@utp.edu.co</email>
                </aff>
                <aff id="aff3">
                    <label>c</label>
                    <institution content-type="original">Universidad Tecnológica de Pereira, Colombia. Correo electrónico: ricardohincapie@utp.edu.co.</institution>
                    <institution content-type="normalized">universidad Tecnológica de Pereira</institution>
                    <institution content-type="orgname">universidad Tecnológica de Pereira</institution>
                    <country country="CO">Colombia</country>
                    <email>ricardohincapie@utp.edu.co</email>
                </aff>
                <!--<pub-date pub-type="epub">
                    <day>20</day>
                    <month>06</month>
                    <year>2019</year>
                </pub-date>
                <pub-date pub-type="collection">-->
            <pub-date pub-type="epub-ppub">
                <season>Jul-Dec</season>
                    <year>2019</year>
                </pub-date>
                <volume>29</volume>
                <issue>2</issue>
                <fpage>115</fpage>
                <lpage>130</lpage>
                <history>
                    <date date-type="received">
                        <day>08</day>
                        <month>09</month>
                        <year>2018</year>
                    </date>
                    <date date-type="accepted">
                        <day>09</day>
                        <month>03</month>
                        <year>2019</year>
                    </date>
                </history>
                <permissions>
                    <license license-type="open-access" xlink:href="https://creativecommons.org/licenses/by-nc-nd/4.0/" xml:lang="en">
                        <license-p>Este obra est&#x00E1; bajo una licencia de Creative Commons Reconocimiento-NoComercial-SinObrasDerivadas 4.0.</license-p>
                    </license>
                </permissions>
                <abstract>
                    <title>
                        <bold>Resumen:</bold>
                    </title>
                    <p>en este art&#x00ED;culo se propone una metodolog&#x00ED;a para la compra y venta de energ&#x00ED;a que considera elementos almacenadores de energ&#x00ED;a en la operaci&#x00F3;n de sistemas de distribuci&#x00F3;n de energ&#x00ED;a el&#x00E9;ctrica, y su objetivo es mejorar las utilidades del operador de red con un adecuado manejo de la curva de carga. En la formulaci&#x00F3;n del problema son propuestos dos modelos matem&#x00E1;ticos cuya funci&#x00F3;n es minimizar el costo de compra de la energ&#x00ED;a el&#x00E9;ctrica. El primer modelo (lineal) incluye la compra de energ&#x00ED;a y el segundo (no lineal) considera, adem&#x00E1;s de lo anterior, las p&#x00E9;rdidas t&#x00E9;cnicas de energ&#x00ED;a del sistema de distribuci&#x00F3;n. Ambos modelos matem&#x00E1;ticos consideran restricciones operativas de los almacenadores de energ&#x00ED;a y de la red de energ&#x00ED;a el&#x00E9;ctrica. La metodolog&#x00ED;a propuesta puede ser usada por los operadores de red como una estrategia para la compra y venta de energ&#x00ED;a; adem&#x00E1;s, proporciona una herramienta matem&#x00E1;tica que puede ser aplicada en la operaci&#x00F3;n y el planeamiento de sistemas de distribuci&#x00F3;n de energ&#x00ED;a, considerando elementos almacenadores de energ&#x00ED;a. La metodolog&#x00ED;a propuesta es verificada empleando dos sistemas de prueba de diferente tama&#x00F1;o, con lo cual se obtienen resultados que corroboran un beneficio econ&#x00F3;mico cuando se usan almacenadores de energ&#x00ED;a en la operaci&#x00F3;n de los sistemas de distribuci&#x00F3;n.</p>
                </abstract>
                <trans-abstract xml:lang="en">
                    <title>
                        <bold>Abstract:</bold>
                    </title>
                    <p>This paper proposes a methodology for the purchase and sale of energy which consider energy storage elements in the operation of distribution systems and whose objective is to improve the utilities of the network operator by the load curve management. In the formulation of the problem are proposed two mathematical models whose objective function is the minimization of the cost of the purchase of electric energy. The first (linear) model includes the purchase of energy and the second (non-linear) considers, in addition to the above, the energy losses of the distribution system. Both mathematical models consider operational restrictions on energy storage and electrical power grid. The proposed methodology can be used by network operators as a strategy for the purchase and sale of energy, additionally is proposed a mathematical tool that can be applied in the operation and planning of distribution systems with energy storage elements. The proposed methodology is verified using test systems of different sizes, obtaining results that corroborate an economic benefit when energy storage is used in the operation of the distribution systems.</p>
                </trans-abstract>
                <kwd-group xml:lang="es">
                    <title>
                        <bold>Palabras clave:</bold>
                    </title>
                    <kwd>almacenadores de energ&#x00ED;a</kwd>
                    <kwd>manejo de la curva de carga</kwd>
                    <kwd>operaci&#x00F3;n de sistemas de distribuci&#x00F3;n;</kwd>
                    <kwd>p&#x00E9;rdidas t&#x00E9;cnicas</kwd>
                    <kwd>sistema de distribuci&#x00F3;n</kwd>
                </kwd-group>
                <kwd-group xml:lang="en">
                    <title>
                        <bold>Keywords:</bold>
                    </title>
                    <kwd>energy storage systems</kwd>
                    <kwd>load curve management</kwd>
                    <kwd>distribution system operation</kwd>
                    <kwd>technical losses</kwd>
                    <kwd>distribution system</kwd>
                </kwd-group>
            </article-meta>
        </front>
<body>
<sec sec-type="sec-1-3635">
<title><bold>Introducci&#x00F3;n</bold></title>
<p>Los sistemas de distribuci&#x00F3;n tienen como objetivo principal entregar la energ&#x00ED;a que proviene de los sistemas de transmisi&#x00F3;n o subtransmisi&#x00F3;n a los usuarios finales, y su operaci&#x00F3;n se ha enfocado en garantizar criterios de calidad, confiabilidad y continuidad del servicio, al menor costo posible. Sin embargo, en los &#x00FA;ltimos a&#x00F1;os han aparecido diversos elementos que afectan los aspectos operativos tradicionales. Uno de estos son los almacenadores de energ&#x00ED;a, los cuales, a pesar de afectar dicha operaci&#x00F3;n, proporcionan una gran variedad de beneficios para los operadores de red.</p>
<p>El uso de almacenadores genera mejoras significativas econ&#x00F3;micas y operativas, entre las que se destacan menores potencias transportadas por las l&#x00ED;neas de transmisi&#x00F3;n, menores p&#x00E9;rdidas de energ&#x00ED;a, mejores niveles de tensi&#x00F3;n y aplazamiento de los proyectos de expansi&#x00F3;n en el tiempo. Adicionalmente, afectan las curvas de carga del sistema, ya que permiten almacenar energ&#x00ED;a en periodos de baja demanda e inyectarla en periodos de alta carga. Este comportamiento ocasiona una correcci&#x00F3;n de la curva de carga por la disminuci&#x00F3;n de los picos de demanda del sistema, lo que se refleja en un alivio en la cargabilidad de los elementos de la red.</p>
<p>El uso de estos elementos resulta de gran inter&#x00E9;s para los operadores de red debido a la disminuci&#x00F3;n de su costo, el avance de nuevas tecnolog&#x00ED;as y los beneficios econ&#x00F3;micos y operativos que generan, por lo que se han convertido en una nueva y atractiva alternativa de inversi&#x00F3;n. Gracias a las ventajas que presentan, en los &#x00FA;ltimos a&#x00F1;os se han realizado diversos trabajos, entre los que se destacan las siguientes referencias.</p>
<p>En [<xref ref-type="bibr" rid="ref-1-3635">1</xref>] se plantea una operaci&#x00F3;n &#x00F3;ptima de los almacenadores de energ&#x00ED;a en sistemas de distribuci&#x00F3;n con generaci&#x00F3;n e&#x00F3;lica utilizando una funci&#x00F3;n multiobjetivo, en la cual uno de los objetivos es minimizar el costo de compra de energ&#x00ED;a el&#x00E9;ctrica. Esto se soluciona usando un algoritmo gen&#x00E9;tico. En [<xref ref-type="bibr" rid="ref-2-3635">2</xref>] se propone un modelo matem&#x00E1;tico que plantea el aplanamiento de la curva de carga, usando almacenadores de energ&#x00ED;a, y cuyo objetivo es disminuir las p&#x00E9;rdidas de energ&#x00ED;a el&#x00E9;ctrica. En [<xref ref-type="bibr" rid="ref-3-3635">3</xref>] se plantea un modelo con varias funciones objetivo que permita satisfacer las necesidades de un operador de red; en este modelo utilizan un flujo de carga trif&#x00E1;sico desbalanceado con almacenadores de energ&#x00ED;a y generaci&#x00F3;n solar para operar de forma &#x00F3;ptima los sistemas de distribuci&#x00F3;n.</p>
<p>En [<xref ref-type="bibr" rid="ref-4-3635">4</xref>] se presentan diversas estrategias para utilizar el despacho econ&#x00F3;mico con almacenadores de energ&#x00ED;a en sistemas de distribuci&#x00F3;n. Para su soluci&#x00F3;n los autores proponen un modelo de optimizaci&#x00F3;n basado en la operaci&#x00F3;n de almacenadores de energ&#x00ED;a en coordinaci&#x00F3;n con energ&#x00ED;as renovables y variaciones de carga. En [<xref ref-type="bibr" rid="ref-5-3635">5</xref>] se realiza un despacho econ&#x00F3;mico de un sistema con energ&#x00ED;as renovables y almacenadores de energ&#x00ED;a. Emplean un modelo matem&#x00E1;tico lineal que es resuelto utilizando un software de optimizaci&#x00F3;n comercial. En [<xref ref-type="bibr" rid="ref-6-3635">6</xref>] se propone un modelo no lineal entero mixto para encontrar la operaci&#x00F3;n &#x00F3;ptima de los almacenadores de energ&#x00ED;a en un sistema de distribuci&#x00F3;n con generaci&#x00F3;n solar. Este problema tiene dos funciones objetivo: el costo de compra de energ&#x00ED;a y costos ambientales por emisiones de gases. Este problema es resuelto utilizando el software de optimizaci&#x00F3;n comercial gams.</p>
<p>En [<xref ref-type="bibr" rid="ref-7-3635">7</xref>] se propone una metodolog&#x00ED;a para minimizar el costo de p&#x00E9;rdidas de energ&#x00ED;a el&#x00E9;ctrica en sistemas de distribuci&#x00F3;n, se utiliza un modelo no lineal que se resuelve con el uso del software de optimizaci&#x00F3;n comercial gams. En [<xref ref-type="bibr" rid="ref-8-3635">8</xref>] se plantea una metodolog&#x00ED;a para operar de forma coordinada los almacenadores de energ&#x00ED;a en sistemas de distribuci&#x00F3;n. En esta metodolog&#x00ED;a los almacenadores mejoran la flexibilidad y eficiencia de la operaci&#x00F3;n de los sistemas de distribuci&#x00F3;n. En [<xref ref-type="bibr" rid="ref-9-3635">9</xref>] proponen una metodolog&#x00ED;a que analiza los impactos econ&#x00F3;micos de los almacenadores de energ&#x00ED;a en los sistemas de distribuci&#x00F3;n para poder generar compensaciones por parte del Gobierno y de los operadores de red para los inversionistas privados.</p>
<p>En [<xref ref-type="bibr" rid="ref-10-3635">10</xref>] se propone una estrategia de aplanamiento de la curva de carga mediante el uso de almacenadores de energ&#x00ED;a y cuyo objetivo es mejorar la confiabilidad y capacidad del sistema. En [<xref ref-type="bibr" rid="ref-11-3635">11</xref>] se usa un modelo probabil&#x00ED;stico de una red neuronal para predecir la curva de demanda y generaci&#x00F3;n fotovoltaica para carga y descarga de los almacenadores de energ&#x00ED;a, con lo cual se obtiene una curva de carga aplanada. En [<xref ref-type="bibr" rid="ref-12-3635">12</xref>] se desarrolla una metodolog&#x00ED;a que permita a los operadores de red manejar los picos de demanda mediante la operaci&#x00F3;n coordinada de almacenadores de energ&#x00ED;a.</p>
<p>En [<xref ref-type="bibr" rid="ref-13-3635">13</xref>] se presenta una estrategia de control de carga y descarga de los almacenadores de energ&#x00ED;a en sistemas de distribuci&#x00F3;n, con el fin de mejorar el valor de compra de la energ&#x00ED;a en el mercado el&#x00E9;ctrico. Por su parte, [<xref ref-type="bibr" rid="ref-14-3635">14</xref>] proponen una coordinaci&#x00F3;n &#x00F3;ptima de almacenadores de energ&#x00ED;a en sistemas de distribuci&#x00F3;n con generaci&#x00F3;n distribuida para tener una operaci&#x00F3;n segura y flexible de la red de distribuci&#x00F3;n.</p>
<p>A diferencia de los trabajos anteriores, este art&#x00ED;culo propone una metodolog&#x00ED;a para la compra y venta de energ&#x00ED;a con un adecuado manejo de la curva de carga haciendo uso de almacenadores. As&#x00ED;, pues, son propuestos dos modelos matem&#x00E1;ticos (lineal y no lineal), cuya funci&#x00F3;n objetivo es la minimizaci&#x00F3;n del costo de compra de la energ&#x00ED;a el&#x00E9;ctrica. La metodolog&#x00ED;a tambi&#x00E9;n considera el impacto en las p&#x00E9;rdidas de energ&#x00ED;a en la red que abastece el sistema de distribuci&#x00F3;n y proporciona una herramienta que puede ser aplicada en la operaci&#x00F3;n y el planeamiento de sistemas de distribuci&#x00F3;n con almacenadores de energ&#x00ED;a.</p>
<p>El primer modelo es lineal y considera la minimizaci&#x00F3;n de compra de energ&#x00ED;a, y est&#x00E1; sujeto a un conjunto de restricciones t&#x00E9;cnicas y operativas. El segundo modelo es no lineal e incluye, adem&#x00E1;s de lo anterior, la reducci&#x00F3;n de las p&#x00E9;rdidas t&#x00E9;cnicas de la red; en este, el conjunto de restricciones es similar al del primer modelo. Para la soluci&#x00F3;n de ambos modelos se emplea el software de optimizaci&#x00F3;n comercial gams.</p>
<p>Con esta metodolog&#x00ED;a se proporciona una herramienta de toma de decisiones para la operaci&#x00F3;n de sistemas de distribuci&#x00F3;n que consideran elementos almacenadores de energ&#x00ED;a. Adem&#x00E1;s de la disminuci&#x00F3;n del valor de compra de energ&#x00ED;a y las p&#x00E9;rdidas t&#x00E9;cnicas de la red, el uso de esta metodolog&#x00ED;a permite mejorar la cargabilidad de las l&#x00ED;neas de transmisi&#x00F3;n, as&#x00ED; como diferir en el tiempo los planes de expansi&#x00F3;n del operador de red.</p>
</sec>
<sec sec-type="sec-2-3635">
<title><bold>Metodolog&#x00ED;a</bold></title>
<sec sec-type="sec-3-3635">
<title><bold>Formulaci&#x00F3;n del problema</bold></title>
<p>El problema de la operaci&#x00F3;n de sistemas de distribuci&#x00F3;n con almacenadores de energ&#x00ED;a se presenta mediante dos modelos matem&#x00E1;ticos (uno lineal y otro no lineal), los cuales se explican detalladamente a continuaci&#x00F3;n.</p>
</sec>
<sec sec-type="sec-4-3635">
<title><bold>Modelo lineal</bold></title>
<p>El modelo matem&#x00E1;tico se presenta en las ecuaciones (1)-(15); en estas, la funci&#x00F3;n objetivo por minimizar considera el valor de costo de compra de energ&#x00ED;a el&#x00E9;ctrica (1). El conjunto de restricciones presentado en las ecuaciones (2)-(15) tiene en cuenta las caracter&#x00ED;sticas t&#x00E9;cnicas y operativas de la red, y de los almacenadores de energ&#x00ED;a.</p>
<p>
<disp-formula id="fig-01-3635">
    <alternatives>
<graphic xlink:href="1909-7735-cein-29-02-115-130-e1.png"/>
</alternatives>
</disp-formula>
</p>
<p>
<disp-formula id="fig-02-3635">
    <alternatives>
<graphic xlink:href="1909-7735-cein-29-02-115-130-e2.png"/>
</alternatives>
</disp-formula>
</p>
<p>
<disp-formula id="fig-03-3635">
    <alternatives>
<graphic xlink:href="1909-7735-cein-29-02-115-130-e3.png"/>
</alternatives>
</disp-formula>
</p>
<p>
<disp-formula id="fig-04-3635">
    <alternatives>
<graphic xlink:href="1909-7735-cein-29-02-115-130-e4.png"/>
</alternatives>
</disp-formula>
</p>
<p>
<disp-formula id="fig-05-3635">
    <alternatives>
<graphic xlink:href="1909-7735-cein-29-02-115-130-e5.png"/>
</alternatives>
</disp-formula>
</p>
<p>
<disp-formula id="fig-06-3635">
    <alternatives>
<graphic xlink:href="1909-7735-cein-29-02-115-130-e6.png"/>
</alternatives>
</disp-formula>
</p>
<p>
<disp-formula id="fig-07-3635">
    <alternatives>
<graphic xlink:href="1909-7735-cein-29-02-115-130-e7.png"/>
</alternatives>
</disp-formula>
</p>
<p>
<disp-formula id="fig-08-3635">
    <alternatives>
<graphic xlink:href="1909-7735-cein-29-02-115-130-e8.png"/>
</alternatives>
</disp-formula>
</p>
<p>
<disp-formula id="fig-09-3635">
    <alternatives>
<graphic xlink:href="1909-7735-cein-29-02-115-130-e9.png"/>
</alternatives>
</disp-formula>
</p>
<p>
<disp-formula id="fig-010-3635">
    <alternatives>
<graphic xlink:href="1909-7735-cein-29-02-115-130-e10.png"/>
</alternatives>
</disp-formula>
</p>
<p>
<disp-formula id="fig-011-3635">
    <alternatives>
<graphic xlink:href="1909-7735-cein-29-02-115-130-e11.png"/>
</alternatives>
</disp-formula>
</p>
<p>
<disp-formula id="fig-012-3635">
    <alternatives>
<graphic xlink:href="1909-7735-cein-29-02-115-130-e12.png"/>
</alternatives>
</disp-formula>
</p>
<p>
<disp-formula id="fig-013-3635">
    <alternatives>
<graphic xlink:href="1909-7735-cein-29-02-115-130-e13.png"/>
</alternatives>
</disp-formula>
</p>
<p>
<disp-formula id="fig-014-3635">
    <alternatives>
<graphic xlink:href="1909-7735-cein-29-02-115-130-e14.png"/>
</alternatives>
</disp-formula>
</p>
<p>
<disp-formula id="fig-015-3635">
    <alternatives>
<graphic xlink:href="1909-7735-cein-29-02-115-130-e15.png"/>
</alternatives>
</disp-formula>
</p>
<p>Las ecuaciones (2) y (3) corresponden al balance de potencia activa y reactiva, y las (4)-(6) describen la existencia de las variables en los nodos. La ecuaci&#x00F3;n (7) es el estado de carga del almacenador de energ&#x00ED;a b conectado al nodo i en el tiempo t. Las ecuaciones (8) y (9) representan el estado inicial y final de carga del almacenador b conectado al nodo i. Las ecuaciones (10) y (11) representan los l&#x00ED;mites de potencia activa y reactiva del generador g en el nodo i. La ecuaci&#x00F3;n (12) es el l&#x00ED;mite de potencia activa de carga y descarga de los almacenadores de energ&#x00ED;a. Las ecuaciones (13) y (14) son los l&#x00ED;mites de capacidad de las l&#x00ED;neas y la (15) son los l&#x00ED;mites de carga de los almacenadores de energ&#x00ED;a.</p>
</sec>
<sec sec-type="sec-5-3635">
<title><bold>Modelo no lineal</bold></title>
<p>Con el fin de involucrar las p&#x00E9;rdidas t&#x00E9;cnicas de energ&#x00ED;a en el modelo presentado en las ecuaciones (1)-(15), se modifican (2), (3), (13) y (14). En (16) y (17) se plantean las ecuaciones de balance nodal de potencia activa y reactiva, y en (18) y (19) los flujos de potencia activa y reactiva por los tramos de red del sistema [<xref ref-type="bibr" rid="ref-16-3635">16</xref>].</p>
<p>Las ecuaciones (2) y (3) son sustituidas por (16) y (17), respectivamente. En estas ecuaciones, dado que la diferencia angular <italic>&#x03B8;<sub>ij</sub></italic> entre dos nodos es peque&#x00F1;a (en radianes), se asume que cos <italic>&#x03B8;<sub>ij</sub> &#x2248; 1</italic> y <italic>sen</italic> <italic>&#x03B8;<sub>ij</sub> &#x2248; &#x03B8;<sub>ij</sub></italic>, por lo que pueden ser escritas como sigue:</p>
<p>
<disp-formula id="fig-016-3635">
    <alternatives>
<graphic xlink:href="1909-7735-cein-29-02-115-130-e16.png"/>
</alternatives>
</disp-formula>
</p>
<p>
<disp-formula id="fig-017-3635">
    <alternatives>
<graphic xlink:href="1909-7735-cein-29-02-115-130-e17.png"/>
</alternatives>
</disp-formula>
</p>
<p>Las ecuaciones (13) y (14) son sustituidas por (18) y (19). Como en el caso anterior, se asume que la diferencia angular <italic>&#x03B8;<sub>ij</sub></italic> es peque&#x00F1;a, y las ecuaciones (18) y (19) son expresadas en funci&#x00F3;n de voltajes nodales y par&#x00E1;metros de la red. Por lo tanto:</p>
<p>
<disp-formula id="fig-018-3635">
    <alternatives>
<graphic xlink:href="1909-7735-cein-29-02-115-130-e18.png"/>
</alternatives>
</disp-formula>
</p>
<p>
<disp-formula id="fig-019-3635">
    <alternatives>
<graphic xlink:href="1909-7735-cein-29-02-115-130-e19.png"/>
</alternatives>
</disp-formula>
</p>
<p>Asimismo, es necesario adicionar la ecuaci&#x00F3;n (20), con el fin de incluir los l&#x00ED;mites de tensi&#x00F3;n de la red.</p>
<p>
<disp-formula id="fig-020-3635">
    <alternatives>
<graphic xlink:href="1909-7735-cein-29-02-115-130-e20.png"/>
</alternatives>
</disp-formula>
</p>
<p>Considerando las ecuaciones anteriores, el modelo no lineal completo se presenta de la siguiente forma:</p>
<p>
<disp-formula id="fig-021-3635">
    <alternatives>
<graphic xlink:href="1909-7735-cein-29-02-115-130-e21.png"/>
</alternatives>
</disp-formula>
</p>
</sec>
<sec sec-type="sec-6-3635">
<title><bold>Metodolog&#x00ED;a propuesta</bold></title>
<p>La metodolog&#x00ED;a propuesta puede ser observada en el diagrama de flujo de la <xref ref-type="fig" rid="fig-1-3635">Figura 1</xref>. Inicialmente se define el sistema de prueba considerando todos los datos que lo describen por completo. Luego se selecciona el tipo de modelo para formular matem&#x00E1;ticamente el problema (lineal o no lineal). M&#x00E1;s adelante, se modela el problema en el software gams. Finalmente, se ejecuta el software y se analizan los resultados obtenidos.</p>
</sec>
</sec>
<sec sec-type="sec-7-3635">
<title><bold>Aplicacion y resultados</bold></title>
<p>Para verificar los modelos matem&#x00E1;ticos propuestos, usan dos sistemas de prueba de 7 y 34 nodos, respectivamente [<xref ref-type="bibr" rid="ref-15-3635">15</xref>]. Con el objetivo de realizar un an&#x00E1;lisis comparativo, en ambos modelos matem&#x00E1;ticos se consideran dos casos de estudio: 1) sin almacenadores de energ&#x00ED;a y 2) con almacenadores de energ&#x00ED;a. Se considera una vida &#x00FA;til de los almacenadores de energ&#x00ED;a de 10 a&#x00F1;os y una tasa de inter&#x00E9;s del 10 %. Todos los casos de estudio fueron implementados en <sc>gams</sc>.</p>
<p>
<fig id="fig-1-3635">
<label><bold>Fig. 1.</bold></label>
<caption><title>Diagrama de flujo que representa la metodolog&#x00ED;a</title></caption>
<graphic xlink:href="1909-7735-cein-29-02-115-130-gf1.png"/>
<attrib>Fuente: elaboraci&#x00F3;n propia.</attrib>
</fig>
</p>
</sec>
<sec sec-type="sec-8-3635">
<title><bold>Sistema de prueba 1</bold></title>
<p>Este sistema tiene 7 nodos y se ilustra en la <xref ref-type="fig" rid="fig-2-3635">Figura 2</xref>, en la cual se muestra que la fuente est&#x00E1; conectada en el nodo 1 y el sistema de almacenadores de energ&#x00ED;a est&#x00E1; conectado en el nodo 2 (subestaci&#x00F3;n). Los nodos 3-7 tienen demandas asociadas y la l&#x00ED;nea entre los nodos 1-2 representa el equivalente de Th&#x00E9;venin de la red conectada aguas arriba del nodo 2 (subestaci&#x00F3;n). En este sistema el costo de los almacenadores de energ&#x00ED;a es de USD 250.000 y su costo anual de operaci&#x00F3;n y mantenimiento, de USD 10.500.</p>
<p>
<fig id="fig-2-3635">
<label><bold>Fig. 2.</bold></label>
<caption><title>Sistema de prueba 1</title></caption>
<graphic xlink:href="1909-7735-cein-29-02-115-130-gf2.png"/>
<attrib>Fuente: elaboraci&#x00F3;n propia.</attrib>
</fig>
</p>
<sec sec-type="sec-9-3635">
<title><italic>Modelo lineal</italic></title>
<p>Al solucionar este modelo se obtienen los resultados presentados en la <xref ref-type="table" rid="tabw-1-3635">Tabla 1</xref>; en esta se observa que la situaci&#x00F3;n m&#x00E1;s econ&#x00F3;mica es aquella en la cual se usan almacenadores de energ&#x00ED;a (caso 2), con una utilidad diaria de USD 238,75 (USD 87.143,75 anuales). Las curvas de carga de ambos casos se presentan en la <xref ref-type="fig" rid="fig-3-3635">Figura 3</xref>, donde las l&#x00ED;neas con c&#x00ED;rculos peque&#x00F1;os y l&#x00ED;neas con x corresponden a los casos 1 y 2, respectivamente.</p>
<p>
<table-wrap id="tabw-1-3635">
<label><bold>Tabla 1.</bold></label>
<caption><title>Resultados del sistema de prueba 1, modelo lineal (USD)</title></caption>
<alternatives>
    <graphic xlink:href="t1.jpg"/>
<table id="tab-1-3635" frame="hsides" border="1" rules="all">
<colgroup width="100%">
<col width="50%"/>
<col width="50%"/>
</colgroup>
<thead>
<tr>
<th valign="middle" align="center" style="background: black"><bold>Caso</bold></th>
<th valign="middle" align="center" style="background: black"><bold>Funci&#x00F3;n objetivo</bold></th>
</tr>
</thead>
<tbody>
<tr>
<td valign="middle" align="center">1</td>
<td valign="middle" align="center">13.503,63</td>
</tr>
<tr>
<td valign="middle" align="center">2</td>
<td valign="middle" align="center">13.264,88</td>
</tr>
</tbody>
</table>
</alternatives>
<attrib>Fuente: elaboraci&#x00F3;n propia.</attrib>
</table-wrap>
</p>
<p>
<fig id="fig-3-3635">
<label><bold>Fig. 3.</bold></label>
<caption><title>Curvas de carga para el sistema de prueba 1, modelo lineal</title></caption>
<graphic xlink:href="1909-7735-cein-29-02-115-130-gf3.png"/>
<attrib>Fuente: elaboraci&#x00F3;n propia.</attrib>
</fig>
</p>
<p>La curva del caso 1 (l&#x00ED;neas con o) es vista aguas abajo de la subestaci&#x00F3;n, y la curva del caso 2 (l&#x00ED;neas con x) es vista por la red aguas arriba de la subestaci&#x00F3;n. La diferencia en energ&#x00ED;a entre las dos curvas de carga corresponde al ciclo de carga y descarga del almacenador de energ&#x00ED;a. De esta figura se obtiene un factor de carga de 0,64 para el caso 1 y de 0,76 para el caso 2, lo cual refleja una mejor distribuci&#x00F3;n de la curva de carga del sistema; de esta manera, se obtiene una curva de carga m&#x00E1;s aplanada.</p>
<p>Respecto a la demanda m&#x00E1;xima presentada en la hora 20, se aprecia que pasa de 7,81 mw en el caso 1 a 6,56 mw en el caso 2, y as&#x00ED; se consigue una disminuci&#x00F3;n de 1,25 mw. Esta reducci&#x00F3;n representa una mejora en la cargabilidad de los elementos del sistema y, por ende, un mejoramiento de las p&#x00E9;rdidas de energ&#x00ED;a el&#x00E9;ctrica de la red.</p>
<p>En la <xref ref-type="fig" rid="fig-4-3635">Figura 4</xref> se presenta un diagrama de flujo con los costos anualizados, que incluye el costo del equipo, su mantenimiento y operaci&#x00F3;n, as&#x00ED; como la utilidad por compra y venta de energ&#x00ED;a. Con base en lo anterior, el costo anualizado del almacenador es de USD 40.686,3. En la figura se evidencia que la ganancia anual (diferencia entre utilidad y costo de inversi&#x00F3;n, operaci&#x00F3;n y mantenimiento) es de USD 35.957,5.</p>
<p>
<fig id="fig-4-3635">
<label><bold>Fig. 4.</bold></label>
<caption><title>Diagrama de flujo del caso 2, modelo lineal (USD)</title></caption>
<graphic xlink:href="1909-7735-cein-29-02-115-130-gf4.png"/>
<attrib>Fuente: elaboraci&#x00F3;n propia.</attrib>
</fig>
</p>
</sec>
<sec sec-type="sec-10-3635">
<title><italic>Modelo no lineal</italic></title>
<p>Con el prop&#x00F3;sito de incluir en el modelo matem&#x00E1;tico el efecto de las p&#x00E9;rdidas de energ&#x00ED;a aguas arriba de la subestaci&#x00F3;n, se usa el equivalente de Th&#x00E9;venin. De esta forma, se incluye en el an&#x00E1;lisis de forma aproximada el efecto de las p&#x00E9;rdidas de energ&#x00ED;a durante las horas del d&#x00ED;a, as&#x00ED; como tambi&#x00E9;n el valor de las p&#x00E9;rdidas en la hora pico. En la <xref ref-type="table" rid="tabw-2-3635">Tabla 2</xref> se presentan los resultados de la funci&#x00F3;n objetivo del modelo no lineal, para ambos casos de estudio.</p>
<p>
<table-wrap id="tabw-2-3635">
<label><bold>Tabla 2.</bold></label>
<caption><title>Resultados del sistema de prueba 1, modelo no lineal (USD)</title></caption>
<alternatives>
    <graphic xlink:href="t2.jpg"/>
<table id="tab-2-3635" frame="hsides" border="1" rules="all">
<colgroup width="100%">
<col width="50%"/>
<col width="50%"/>
</colgroup>
<thead>
<tr>
<th valign="middle" align="center" style="background: black"><bold>Caso</bold></th>
<th valign="middle" align="center" style="background: black"><bold>Funci&#x00F3;n objetivo</bold></th>
</tr>
</thead>
<tbody>
<tr>
<td valign="middle" align="center">1</td>
<td valign="middle" align="center">13.619,98</td>
</tr>
<tr>
<td valign="middle" align="center">2</td>
<td valign="middle" align="center">13.378,39</td>
</tr>
</tbody>
</table>
</alternatives>
<attrib>Fuente: elaboraci&#x00F3;n propia.</attrib>
</table-wrap>
</p>
<p>En la <xref ref-type="table" rid="tabw-2-3635">Tabla 2</xref> se observa que los valores obtenidos con el modelo no lineal son mayores que los obtenidos con el modelo lineal, por causa del efecto de las p&#x00E9;rdidas; asimismo, que el caso m&#x00E1;s econ&#x00F3;mico es el que usa almacenadores de energ&#x00ED;a (caso 2), con una utilidad diaria de USD 241,59 (USD 88.180,4 anuales). Este valor es mayor al del modelo lineal debido al efecto de las p&#x00E9;rdidas aguas arriba de la subestaci&#x00F3;n.</p>
<p>En la <xref ref-type="fig" rid="fig-5-3635">Figura 5</xref> se muestran las curvas de carga con la energ&#x00ED;a que entrega la red, para ambos casos de estudio. El factor de carga para los casos 1 y 2 es 0,63 y 0,75 respectivamente, lo cual indica un aplanamiento de la curva de carga por a la presencia de los almacenadores de energ&#x00ED;a. En esta misma figura se muestra que la potencia m&#x00E1;xima en la hora 20 pasa de 7,91 mw (caso 1) a 6,65 mw (caso 2), con una disminuci&#x00F3;n de 1,26 mw. El incremento de este valor respecto al modelo anterior es 8 kW, lo cual se debe al valor considerado de las p&#x00E9;rdidas de la red.</p>
<p>
<fig id="fig-5-3635">
<label><bold>Fig. 5.</bold></label>
<caption><title>Curvas de carga para el sistema de prueba 1, modelo no lineal</title></caption>
<graphic xlink:href="1909-7735-cein-29-02-115-130-gf5.png"/>
<attrib>Fuente: elaboraci&#x00F3;n propia.</attrib>
</fig>
</p>
<p>En la <xref ref-type="fig" rid="fig-6-3635">Figura 6</xref> se describen los costos anualizados de la inversi&#x00F3;n inicial, as&#x00ED; como el costo de mantenimiento y operaci&#x00F3;n. Tambi&#x00E9;n se presenta la ganancia por compra y venta de energ&#x00ED;a, con lo que se obtiene un valor anual de USD 36.994,1.</p>
<p>
<fig id="fig-6-3635">
<label><bold>Fig. 6.</bold></label>
<caption><title>Diagrama de flujo del caso 2, modelo no lineal (USD)</title></caption>
<graphic xlink:href="1909-7735-cein-29-02-115-130-gf6.png"/>
<attrib>Fuente: elaboraci&#x00F3;n propia.</attrib>
</fig>
</p>
</sec>
</sec>
<sec sec-type="sec-11-3635">
<title><bold>Sistema de prueba 2</bold></title>
<p>Este sistema tiene 34 nodos y se ilustra en la <xref ref-type="fig" rid="fig-7-3635">Figura 7</xref>, la cual evidencia que la fuente est&#x00E1; conectada en el nodo 1 y el sistema de almacenadores de energ&#x00ED;a, en el nodo 2 (subestaci&#x00F3;n). Los nodos 3-34 tienen demandas asociadas y la l&#x00ED;nea entre los nodos 1-2 representa el equivalente de Th&#x00E9;venin de la red conectada aguas arriba del nodo 2 (subestaci&#x00F3;n). En este sistema el costo de los almacenadores de energ&#x00ED;a es de USD 150.000 y su costo anual de operaci&#x00F3;n y mantenimiento, de USD 6.250.</p>
<p>
<fig id="fig-7-3635">
<label><bold>Fig. 7.</bold></label>
<caption><title>Sistema de prueba 2</title></caption>
<graphic xlink:href="1909-7735-cein-29-02-115-130-gf7.png"/>
<attrib>Fuente: elaboraci&#x00F3;n propia.</attrib>
</fig>
</p>
<sec sec-type="sec-12-3635">
<title><italic>Modelo lineal</italic></title>
<p>Al solucionar este modelo se obtienen los resultados de la <xref ref-type="table" rid="tabw-1-3635">Tabla 1</xref>; en esta se observa que la situaci&#x00F3;n m&#x00E1;s econ&#x00F3;mica es aquella en la cual se usan almacenadores de energ&#x00ED;a (caso 2), con una utilidad diaria de USD 143,25 (USD 52.286,25 anuales). Las curvas de carga de ambos casos se ilustran en la <xref ref-type="fig" rid="fig-8-3635">Figura 8</xref>, donde las l&#x00ED;neas con c&#x00ED;rculos peque&#x00F1;os y l&#x00ED;neas con x corresponden a los casos 1 y 2, respectivamente.</p>
<p>
<table-wrap id="tabw-3-3635">
<label><bold>Tabla 3.</bold></label>
<caption><title>Resultados del sistema de prueba 2, modelo lineal (USD)</title></caption>
<alternatives>
    <graphic xlink:href="t3.jpg"/>
<table id="tab-3-3635" frame="hsides" border="1" rules="all">
<colgroup width="100%">
<col width="50%"/>
<col width="50%"/>
</colgroup>
<thead>
<tr>
<th valign="middle" align="center" style="background: black"><bold>Caso</bold></th>
<th valign="middle" align="center" style="background: black"><bold>Funci&#x00F3;n objetivo</bold></th>
</tr>
</thead>
<tbody>
<tr>
<td valign="middle" align="center">1</td>
<td valign="middle" align="center">5.849,70</td>
</tr>
<tr>
<td valign="middle" align="center">2</td>
<td valign="middle" align="center">5.706,45</td>
</tr>
</tbody>
</table>
</alternatives>
<attrib>Fuente: elaboraci&#x00F3;n propia.</attrib>
</table-wrap>
</p>
<p>
<fig id="fig-8-3635">
<label><bold>Fig. 8.</bold></label>
<caption><title>Curvas de carga para el sistema de prueba 2 , modelo lineal</title></caption>
<graphic xlink:href="1909-7735-cein-29-02-115-130-gf8.png"/>
<attrib>Fuente: elaboraci&#x00F3;n propia.</attrib>
</fig>
</p>
<p>De la <xref ref-type="fig" rid="fig-8-3635">Figura 8</xref> se obtiene un factor de carga de 0,64 para el caso 1 y 0,74 para el caso 2; esto refleja una mejor distribuci&#x00F3;n de la curva de carga del sistema, y se obtiene as&#x00ED; una curva de carga m&#x00E1;s aplanada.</p>
<p>Respecto a la demanda m&#x00E1;xima presentada en la hora 20, se aprecia que pasa de 3,71 en el caso 1 a 2,96 mw en el caso 2; as&#x00ED;, se consigue una disminuci&#x00F3;n de 0,75 mw. Esta reducci&#x00F3;n representa una mejora en la cargabilidad de los elementos del sistema y, por ende, un mejoramiento de las p&#x00E9;rdidas de energ&#x00ED;a el&#x00E9;ctrica de la red.</p>
<p>En la <xref ref-type="fig" rid="fig-9-3635">Figura 9</xref> se muestra un diagrama de flujo con los costos anualizados, que incluye el costo del equipo, su mantenimiento y operaci&#x00F3;n, as&#x00ED; como la ganancia por compra y venta de energ&#x00ED;a. Con base en lo anterior, el costo anualizado del almacenador es de USD 24.411,8. En esta figura se aprecia que la ganancia anual (diferencia entre utilidad y costo de inversi&#x00F3;n, operaci&#x00F3;n y mantenimiento) es de USD 21.624,5.</p>
<p>
<fig id="fig-9-3635">
<label><bold>Fig. 9.</bold></label>
<caption><title>Diagrama de flujo del caso 2, modelo lineal (USD)</title></caption>
<graphic xlink:href="1909-7735-cein-29-02-115-130-gf9.png"/>
<attrib>Fuente: elaboraci&#x00F3;n propia.</attrib>
</fig>
</p>
</sec>
<sec sec-type="sec-13-3635">
<title><italic>Modelo no lineal</italic></title>
<p>Para incluir en el modelo matem&#x00E1;tico el efecto de las p&#x00E9;rdidas de energ&#x00ED;a aguas arriba de la subestaci&#x00F3;n, se usa el equivalente de Th&#x00E9;venin. As&#x00ED;, es comprendido en el an&#x00E1;lisis de forma aproximada el efecto de las p&#x00E9;rdidas de energ&#x00ED;a durante las horas del d&#x00ED;a, as&#x00ED; como el valor de las p&#x00E9;rdidas en la hora pico. En la <xref ref-type="bibr" rid="ref-4-3635">Tabla 4</xref> se presentan los resultados de la funci&#x00F3;n objetivo del modelo no lineal, para ambos casos de estudio.</p>
<p>
<table-wrap id="tabw-4-3635">
<label><bold>Tabla 4.</bold></label>
<caption><title>Resultados del sistema de prueba 2, modelo no lineal (USD)</title></caption>
<alternatives>
    <graphic xlink:href="t4.jpg"/>
<table id="tab-4-3635" frame="hsides" border="1" rules="all">
<colgroup width="100%">
<col width="50%"/>
<col width="50%"/>
</colgroup>
<thead>
<tr>
<th valign="middle" align="center" style="background: black"><bold>Caso</bold></th>
<th valign="middle" align="center" style="background: black"><bold>Funci&#x00F3;n objetivo</bold></th>
</tr>
</thead>
<tbody>
<tr>
<td valign="middle" align="center">1</td>
<td valign="middle" align="center">6.184,39</td>
</tr>
<tr>
<td valign="middle" align="center">2</td>
<td valign="middle" align="center">6.040,01</td>
</tr>
</tbody>
</table>
</alternatives>
<attrib>Fuente: elaboraci&#x00F3;n propia.</attrib>
</table-wrap>
</p>
<p>De la <xref ref-type="bibr" rid="ref-4-3635">Tabla 4</xref> se puede concluir que los valores obtenidos con el modelo no lineal son mayores que los obtenidos con el modelo lineal, debido al efecto de las p&#x00E9;rdidas. Tambi&#x00E9;n se observa que el caso m&#x00E1;s econ&#x00F3;mico es el que usa almacenadores de energ&#x00ED;a (caso 2), con una utilidad diaria de USD 144,38 (USD 52.698,7 anuales). Este valor es mayor al del modelo lineal por el efecto de las p&#x00E9;rdidas aguas arriba de la subestaci&#x00F3;n.</p>
<p>En la <xref ref-type="fig" rid="fig-10-3635">Figura 10</xref> se muestran las curvas de carga con la energ&#x00ED;a que entrega la red, para ambos casos de estudio. El factor de carga para los casos 1 y 2 es 0,64 y 0,77, respectivamente, lo cual indica un aplanamiento de la curva de carga por causa de la presencia de los almacenadores de energ&#x00ED;a. En esta misma figura se observa que la potencia m&#x00E1;xima en la hora 20 pasa de 3,89 mw (caso 1) a 3,14 mw (caso 2), con una disminuci&#x00F3;n de 0,75 mw. El incremento de este valor respecto al modelo lineal es 2 kW, lo cual se debe al valor considerado de las p&#x00E9;rdidas de la red.</p>
<p>
<fig id="fig-10-3635">
<label><bold>Fig. 10.</bold></label>
<caption><title>Curvas de carga para el sistema de prueba 2, modelo no lineal</title></caption>
<graphic xlink:href="1909-7735-cein-29-02-115-130-gf10.png"/>
<attrib>Fuente: elaboraci&#x00F3;n propia.</attrib>
</fig>
</p>
<p>En la <xref ref-type="fig" rid="fig-11-3635">Figura 11</xref> se presenta un diagrama de flujo con los costos anualizados, que incluye el costo del equipo, su mantenimiento y operaci&#x00F3;n, as&#x00ED; como la ganancia por compra y venta de energ&#x00ED;a. En esta figura se aprecia que la ganancia anual (diferencia entre utilidad y costo de inversi&#x00F3;n, operaci&#x00F3;n y mantenimiento) es de USD 22.036,9.</p>
<p>
<fig id="fig-11-3635">
<label><bold>Fig. 11.</bold></label>
<caption><title>Diagrama de flujo del caso 2, modelo no lineal (USD)</title></caption>
<graphic xlink:href="1909-7735-cein-29-02-115-130-gf11.png"/>
<attrib>Fuente: elaboraci&#x00F3;n propia.</attrib>
</fig>
</p>
</sec>
</sec>
<sec sec-type="sec-14-3635">
<title><bold>Comentarios finales</bold></title>
<p>Para trabajos futuros se puede considerar la localizaci&#x00F3;n &#x00F3;ptima de sistemas de almacenamiento de energ&#x00ED;a en sistemas de distribuci&#x00F3;n utilizando el modelo no lineal, que permite representar las p&#x00E9;rdidas de energ&#x00ED;a. Este nuevo modelo ser&#x00ED;a un problema de programaci&#x00F3;n no lineal entera mixta (PNLEM), que podr&#x00ED;a ser resuelto mediante el uso de una metaheur&#x00ED;stica, utilizando como flujo de carga el modelo no lineal planteado en este art&#x00ED;culo. Tambi&#x00E9;n es posible incluir generaci&#x00F3;n distribuida en este problema, agregando las restricciones generadas por cada generador, dependiendo de su tecnolog&#x00ED;a; de esta forma, surge un problema m&#x00E1;s general que tiene en cuenta almacenadores de energ&#x00ED;a y energ&#x00ED;as renovables. Adem&#x00E1;s, se puede adicionar esta metodolog&#x00ED;a con la propuesta de [<xref ref-type="bibr" rid="ref-17-3635">17</xref>], para considerar el planeamiento integrado de sistemas dedistribuci&#x00F3;n con almacenadores de energ&#x00ED;a, el cual podr&#x00ED;a ser resuelto mediante el uso de una metaheur&#x00ED;stica, utilizando como flujo de carga el modelo no lineal o lineal planteado en este art&#x00ED;culo.</p>
<p>El an&#x00E1;lisis de la ubicaci&#x00F3;n de almacenadores de energ&#x00ED;a en transmisi&#x00F3;n debido al impacto de los almacenadores ubicados en distribuci&#x00F3;n tambi&#x00E9;n puede ser estudiado a partir de un equivalente de Th&#x00E9;venin para transmisi&#x00F3;n y generando una metodolog&#x00ED;a bi-nivel. Este ser&#x00ED;a un problema interesante dado que en transmisi&#x00F3;n se busca mejorar solo los picos de alta demanda, mientras que en distribuci&#x00F3;n se ve desde la perspectiva de compra y venta de energ&#x00ED;a el&#x00E9;ctrica y puede que estos dos objetivos se encuentren en conflicto.</p>
</sec>
<sec sec-type="sec-15-3635">
<title><bold>Conclusiones</bold></title>
<p>Se propone una metodolog&#x00ED;a para resolver el problema de la operaci&#x00F3;n en sistemas de distribuci&#x00F3;n con almacenadores de energ&#x00ED;a. El modelo matem&#x00E1;tico planteado considera como funci&#x00F3;n objetivo la compra de energ&#x00ED;a y las p&#x00E9;rdidas de energ&#x00ED;a en la red, y est&#x00E1; sujeto a restricciones t&#x00E9;cnicas y operativas de la red y los almacenadores. Son propuestos, entonces, un modelo lineal y uno no lineal; el primeroconsidera la compra de energ&#x00ED;a, y el segundo, adem&#x00E1;s de lo anterior, las p&#x00E9;rdidas t&#x00E9;cnicas del sistema. Los modelos son analizados con dos casos de prueba y resueltos con el software de optimizaci&#x00F3;n comercial gams.</p>
<p>A pesar de que con los dos modelos se obtienen soluciones relativamente semejantes, el modelo no lineal presenta resultados m&#x00E1;s pr&#x00F3;ximos a los de la operaci&#x00F3;n real de las redes el&#x00E9;ctricas, ya que tiene encuenta para su an&#x00E1;lisis el efecto de las p&#x00E9;rdidas. Las p&#x00E9;rdidas de la red de transmisi&#x00F3;n son representadas a trav&#x00E9;s de un equivalente de Th&#x00E9;venin localizado en la subestaci&#x00F3;n. El modelo lineal es de f&#x00E1;cil implementaci&#x00F3;n y bajo requerimiento de c&#x00F3;mputo, por lo que tambi&#x00E9;n resulta atractiva su implementaci&#x00F3;n.</p>
<p>La instalaci&#x00F3;n de almacenadores de energ&#x00ED;a es un negocio viable para los operadores de red y representa una alta rentabilidad. El tiempo de repago del equipo almacenador es bajo. Anualmente son obtenidas altas ganancias, como consecuencia de la compra y venta de energ&#x00ED;a; adem&#x00E1;s, se logra una mejora en el nivel de p&#x00E9;rdidas de energ&#x00ED;a el&#x00E9;ctrica y de confiabilidad en la red de transmisi&#x00F3;n.</p>
</sec>
</body>
<back>
<ref-list>
<ref id="ref-1-3635">
<mixed-citation>M. D. Zhang, J. Xie y D. Q. Gan, “Optimal ESS operation
strategy in distribution system with wind
power,” en International Conference on Sustainable
Power Generation and Supply (SUPERGEN 2012),
Hangzhou, China, sep. 2012. DOI: https://doi.
org/10.1049/cp.2012.1791</mixed-citation>
<element-citation publication-type="book">
<person-group person-group-type="author">
<name>
     <surname>Zhang</surname>
    <given-names>M. D.</given-names>
    </name>
<name>
     <surname>Xie</surname>
    <given-names>J.</given-names>
    </name>
<name>
     <surname>Gan</surname>
    <given-names>D. Q.</given-names>
    </name>
    </person-group>
<article-title>Optimal ESS operation strategy in distribution system with wind power,</article-title>
<source><italic>International Conference on Sustainable Power Generation and Supply (SUPERGEN 2012)</italic></source>
<publisher-loc>Hangzhou, China</publisher-loc>
<month>09</month>
<year>2012</year>
<pub-id pub-id-type="doi">10.1049/cp.2012.1791</pub-id>
</element-citation>
</ref>
<ref id="ref-2-3635">
<mixed-citation>H. Saboori y H. Abdi, “Application of a grid scale energy
storage system to reduce distribution network losses,”
en 18th Electric Power Distribution Conference,
Kermanshah, Iran, abr.-may. 2013. DOI: https://doi.
org/ 10.1109/EPDC.2013.6565963</mixed-citation>
<element-citation publication-type="book">
<person-group person-group-type="author">
<name>
     <surname>Saboori</surname>
    <given-names>H.</given-names>
    </name>
<name>
     <surname>Abdi</surname>
    <given-names>H.</given-names>
 </name>
</person-group>
<article-title>Application of a grid scale energy storage system to reduce distribution network losses,</article-title>
<source><italic>18th Electric Power Distribution Conference</italic></source>
<publisher-loc>Kermanshah, Iran</publisher-loc>
<comment>abr.</comment>
<month>05</month>
<year>2013</year>
<pub-id pub-id-type="doi">10.1109/EPDC.2013.6565963</pub-id>
</element-citation>
</ref>
<ref id="ref-3-3635">
<mixed-citation>M. C. Bozchalui y R. Sharma, “Optimal operation of
energy storage in distribution systems with renewable
energy resources,” en Clemson University Power Systems
Conference, Clemson, Estados Unidos, mar. 2014.
DOI: https://doi.org/10.1109/PSC.2014.6808125</mixed-citation>
<element-citation publication-type="book">
<person-group person-group-type="author">
<name>
     <surname>Bozchalui</surname>
    <given-names>M. C.</given-names>
    </name>
<name>
     <surname>Sharma</surname>
    <given-names>R.</given-names>
 </name>
</person-group>
<article-title>Optimal operation of energy storage in distribution systems with renewable energy resources,</article-title>
<source><italic>Clemson University Power Systems Conference</italic></source>
<publisher-loc>Clemson, Estados Unidos</publisher-loc>
<month>03</month>
<year>2014</year>
<pub-id pub-id-type="doi">10.1109/PSC.2014.6808125</pub-id>
</element-citation>
</ref>
<ref id="ref-4-3635">
<mixed-citation>M. C. Bozchalui y R. Sharma, “Operation strategies for
energy storage systems in distribution networks,” en
IEEE Power and Energy Society General Meeting, National
Harbor, Estados Unidos, jul. 2014. DOI: https://
doi.org/10.1109/PESGM.2014.6939483</mixed-citation>
<element-citation publication-type="book">
<person-group person-group-type="author">
<name>
     <surname>Bozchalui</surname>
    <given-names>M. C.</given-names>
    </name>
<name>
     <surname>Sharma</surname>
    <given-names>R.</given-names>
 </name>
</person-group>
<article-title>Operation strategies for energy storage systems in distribution networks,</article-title>
<source><italic>IEEE Power and Energy Society General Meeting</italic></source>
<publisher-loc>National Harbor, Estados Unidos</publisher-loc>
<month>07</month>
<year>2014</year>
<pub-id pub-id-type="doi">10.1109/PESGM.2014.6939483</pub-id>
</element-citation>
</ref>
<ref id="ref-5-3635">
<mixed-citation>A. C. Luna et al., “Economic power dispatch of distributed
generators in a grid-connected microgrid,” en
9th International Conference on Power Electronics and
ECCE Asia (ICPE-ECCE Asia), Seúl, Corea del Sur, jul.
2015. DOI: https://doi.org/10.1109/ICPE.2015.7167927</mixed-citation>
<element-citation publication-type="book">
<person-group person-group-type="author">
<name>
     <surname>Luna</surname>
    <given-names>A. C.</given-names>
    </name>
</person-group>
<article-title>Economic power dispatch of distributed generators in a grid-connected microgrid,</article-title>
<source><italic>9th International Conference on Power Electronics and ECCE Asia (ICPE-ECCE Asia)</italic></source>
<publisher-loc>Se&#x00FA;l, Corea del Sur</publisher-loc>
<month>07</month>
<year>2015</year>
<pub-id pub-id-type="doi">10.1109/ICPE.2015.7167927</pub-id>
</element-citation>
</ref>
<ref id="ref-6-3635">
<mixed-citation>R. Bao y W. Zhang, “Optimal scheduling of active
distribution networks with energy storage systems
considering environmental cost,” en 5th International
Conference on Electric Utility Deregulation
and Restructuring and Power Technologies (DRPT),
Changsha, China, nov. 2015. Doi: https://doi.
org/10.1109/DRPT.2015.7432321</mixed-citation>
<element-citation publication-type="book">
<person-group person-group-type="author">
<name>
     <surname>Bao</surname>
    <given-names>R.</given-names>
    </name>
<name>
     <surname>Zhang</surname>
    <given-names>W.</given-names>
 </name>
</person-group>
<article-title>Optimal scheduling of active distribution networks with energy storage systems considering environmental cost,</article-title>
<source><italic>5th International Conference on Electric Utility Deregulation and Restructuring and Power Technologies (DRPT)</italic></source>
<publisher-loc>Changsha, China</publisher-loc>
<month>11</month>
<year>2015</year>
<pub-id pub-id-type="doi">10.1109/DRPT.2015.7432321</pub-id>
</element-citation>
</ref>
<ref id="ref-7-3635">
<mixed-citation>A. Soroudi, P. Siano y A. Keane, “Optimal DR and ESS
scheduling for distribution losses payments minimization
under electricity price uncertainty,” IEEE Transactions
on Smart Grid, vol. 7, no. 1, pp. 261-272, 2015. DOI:
https://doi.org/ 10.1109/TSG.2015.2453017</mixed-citation>
<element-citation publication-type="journal">
<person-group person-group-type="author">
<name>
     <surname>Soroudi</surname>
    <given-names>A.</given-names>
    </name>
<name>
     <surname>Siano</surname>
    <given-names>P.</given-names>
    </name>
<name>
     <surname>Keane</surname>
    <given-names>A.</given-names>
 </name>
</person-group>
<article-title>Optimal DR and ESS scheduling for distribution losses payments minimization under electricity price uncertainty,</article-title>
<source><italic>IEEE Transactions on Smart Grid</italic></source>
<volume>7</volume><issue>1</issue>
<fpage>261</fpage>
<lpage>272</lpage>
<year>2015</year>
<pub-id pub-id-type="doi">10.1109/TSG.2015.2453017</pub-id>
</element-citation>
</ref>
<ref id="ref-8-3635">
<mixed-citation>I. Miranda, H. Leite y N. Silva, “Coordination of multifunctional
distributed energy storage systems in
distribution networks,” IET Generation, Transmission
&amp; Distribution, vol. 10, no. 3, pp. 726-735, 2016. DOI:
https://doi.org/10.1049/iet-gtd.2015.0398</mixed-citation>
<element-citation publication-type="journal">
<person-group person-group-type="author">
<name>
     <surname>Miranda</surname>
    <given-names>I.</given-names>
    </name>
<name>
     <surname>Leite</surname>
    <given-names>H.</given-names>
    </name>
<name>
     <surname>Silva</surname>
    <given-names>N.</given-names>
 </name>
</person-group>
<article-title>Coordination of multifunctional distributed energy storage systems in distribution networks,</article-title>
<source><italic>IET Generation, Transmission &#x0026; Distribution</italic></source>
<volume>10</volume><issue>3</issue>
<fpage>726</fpage>
<lpage>735</lpage> 
<year>2016</year>
<pub-id pub-id-type="doi">10.1049/iet-gtd.2015.0398</pub-id>
</element-citation>
</ref>
<ref id="ref-9-3635">
<mixed-citation>J. Ye, J. Xue, B. Sang, D. Lu y H. Liu, “Economic value
and government compensation calculative method of
energy storage system,” en IEEE 8th International Power
Electronics and Motion Control Conference (IPEMC-ECCE
Asia), Beijing, China, may. 2016. DOI: https://doi.
org/10.1109/IPEMC.2016.7512343</mixed-citation>
<element-citation publication-type="book">
<person-group person-group-type="author">
<name>
     <surname>Ye</surname>
    <given-names>J.</given-names>
    </name>
<name>
     <surname>Xue</surname>
    <given-names>J.</given-names>
    </name>
<name>
     <surname>Sang</surname>
    <given-names>B.</given-names>
    </name>
<name>
     <surname>Lu</surname>
    <given-names>D.</given-names>
    </name>
<name>
     <surname>Liu</surname>
    <given-names>H.</given-names>
 </name>
</person-group>
<article-title>Economic value and government compensation calculative method of energy storage system,</article-title>
<source><italic>IEEE 8th International Power Electronics and Motion Control Conference (IPEMC-ECCE Asia)</italic></source>
<publisher-loc>Beijing, China</publisher-loc>
<month>05</month>
<year>2016</year>
<pub-id pub-id-type="doi">10.1109/IPEMC.2016.7512343</pub-id>
</element-citation>
</ref>
<ref id="ref-10-3635">
<mixed-citation>R. W. Yan y H.-M. Xiao, “Energy storage system optimization
strategy in the distribution network based
on active set method,” en IEEE International Conferences
on Big Data and Cloud Computing (BDCloud),
Social Computing and Networking (SocialCom), Sustainable
Computing and Communications (SustainCom)
(BDCloud-SocialCom-SustainCom), Atlanta, Estados
Unidos, oct. 2016. DOI: https://doi.org/10.1109/BDCloud-
SocialCom-SustainCom.2016.45</mixed-citation>
<element-citation publication-type="book">
<person-group person-group-type="author">
<name>
     <surname>Yan</surname>
    <given-names>R. W.</given-names>
    </name>
<name>
     <surname>Xiao</surname>
    <given-names>H.-M.</given-names>
 </name>
</person-group>
<article-title>Energy storage system optimization strategy in the distribution network based on active set method,</article-title>
<source><italic>IEEE International Conferences on Big Data and Cloud Computing (BDCloud), Social Computing and Networking (SocialCom), Sustainable Computing and Communications (SustainCom) (BDCloud-SocialCom-SustainCom)</italic></source>
<publisher-loc>Atlanta, Estados Unidos</publisher-loc>
<month>10</month>
<year>2016</year>
<pub-id pub-id-type="doi">10.1109/BDCloud-SocialCom-SustainCom.2016.45</pub-id>
</element-citation>
</ref>
<ref id="ref-11-3635">
<mixed-citation>B. R. Ke, T. T. Ku, Y. L. Ke, C. Y. Chung y H. Z. Chen,
“Sizing the battery energy storage system on a university
campus with prediction of load and photovoltaic
generation,” IEEE Transactions on Industry Applications,
vol. 52, no. 2, pp. 1136-1147, 2016. DOI: https://
doi.org/10.1109/TIA.2015.2483583</mixed-citation>
<element-citation publication-type="journal">
<person-group person-group-type="author">
<name>
     <surname>Ke</surname>
    <given-names>B. R.</given-names>
    </name>
<name>
     <surname>Ku</surname>
    <given-names>T. T.</given-names>
    </name>
<name>
     <surname>Ke</surname>
    <given-names>Y. L.</given-names>
    </name>
<name>
     <surname>Chung</surname>
    <given-names>C. Y.</given-names>
    </name>
<name>
     <surname>Chen</surname>
    <given-names>H. Z.</given-names>
    </name>
</person-group>
<article-title>Sizing the battery energy storage system on a university campus with prediction of load and photovoltaic generation,</article-title>
<source><italic>IEEE Transactions on Industry Applications</italic></source>
<volume>52</volume><issue>2</issue>
<fpage>1136</fpage>
<lpage>1147</lpage>
<year>2016</year>
<pub-id pub-id-type="doi">10.1109/TIA.2015.2483583</pub-id>
</element-citation>
</ref>
<ref id="ref-12-3635">
<mixed-citation>M. Cheng, R. Moreira y Q. Yang, “Benefit analysis of
coordinated operation of distributed energy storage:
The distributed network operator’s perspective,” en
IEEE Conference on Energy Internet and Energy System
Integration (EI2), Hefei, China, nov. 2017. DOI: https://
doi.org/10.1109/EI2.2017.8245539</mixed-citation>
<element-citation publication-type="book">
<person-group person-group-type="author">
<name>
     <surname>Cheng</surname>
    <given-names>M.</given-names>
    </name>
<name>
     <surname>Moreira</surname>
    <given-names>R.</given-names>
    </name>
<name>
     <surname>Yang</surname>
    <given-names>Q.</given-names>
    </name>
</person-group><article-title>Benefit analysis of coordinated operation of distributed energy storage: The distributed network operator&#x2019;s perspective,</article-title>
<source><italic>IEEE Conference on Energy Internet and Energy System Integration (EI2)</italic></source>
<publisher-loc>Hefei, China</publisher-loc>
<month>11</month>
<year>2017</year>
<pub-id pub-id-type="doi">10.1109/EI2.2017.8245539</pub-id>
</element-citation>
</ref>
<ref id="ref-13-3635">
<mixed-citation>Y. Zheng et al. “Optimal operation of battery energy
storage system considering distribution system uncertainty,”
IEEE Transactions on Sustainable Energy,
vol. 9, no. 3, pp. 1051-1060, 2018. DOI: https://doi.
org/10.1109/TSTE.2017.2762364</mixed-citation>
<element-citation publication-type="journal">
<person-group person-group-type="author">
<name>
     <surname>Zheng</surname>
    <given-names>Y.</given-names>
    </name>
</person-group>
<article-title>Optimal operation of battery energy storage system considering distribution system uncertainty,</article-title>
<source><italic>IEEE Transactions on Sustainable Energy</italic></source>
<volume>9</volume><issue>3</issue>
<fpage>1051</fpage>
<lpage>1060</lpage>
<year>2018</year>
<pub-id pub-id-type="doi">10.1109/TSTE.2017.2762364</pub-id>
</element-citation>
</ref>
<ref id="ref-14-3635">
<mixed-citation>T. Wang, “An analytical model of distributed energy
storage systems in power distribution networks,” en
17th International Symposium on Distributed Computing
and Applications for Business Engineering and Science
(DCABES), Wuxi, China, oct. 2018. DOI: https://doi.
org/10.1109/DCABES.2018.00011</mixed-citation>
<element-citation publication-type="book">
<person-group person-group-type="author">
<name>
     <surname>Wang</surname>
    <given-names>T.</given-names>
 </name>
</person-group>
<article-title>An analytical model of distributed energy storage systems in power distribution networks,</article-title>
<source><italic>17th International Symposium on Distributed Computing and Applications for Business Engineering and Science (DCABES)</italic></source>
<publisher-loc>Wuxi, China</publisher-loc>
<month>10</month>
<year>2018</year>
<pub-id pub-id-type="doi">10.1109/DCABES.2018.00011</pub-id>
</element-citation>
</ref>
<ref id="ref-15-3635">
<mixed-citation>A. Valencia, “Impacto de los almacenadores de energía
en la operación de sistemas de distribución,” tesis de
pregrado, Facultad de Ingenieria, Universidad Tecnológica
de Pereira, Colombia, 2017.</mixed-citation>
<element-citation publication-type="book">
<person-group person-group-type="author">
<name>
     <surname>Valencia</surname>
    <given-names>A.</given-names>
 </name>
</person-group>
<article-title>Impacto de los almacenadores de energ&#x00ED;a en la operaci&#x00F3;n de sistemas de distribuci&#x00F3;n,</article-title>
<comment>tesis de pregrado, Facultad de Ingenieria, Universidad Tecnol&#x00F3;gica de Pereira</comment>
<publisher-loc>Colombia</publisher-loc>
<year>2017</year>
</element-citation>
</ref>
<ref id="ref-16-3635">
<mixed-citation>R. A. Gallego, A. H. Escobar y M. Granada. Flujo de
carga en sistemas de transmision. Pereira, Colombia:
Editorial UTP, 2016.</mixed-citation>
<element-citation publication-type="book">
<person-group person-group-type="author">
<name>
     <surname>Gallego</surname>
    <given-names>R. A.</given-names>
    </name>
<name>
     <surname>Escobar</surname>
    <given-names>A. H.</given-names>
    </name>
<name>
     <surname>Granada</surname>
    <given-names>M.</given-names>
    </name>
</person-group>
<source xml:lang="es"><italic>Flujo de carga en sistemas de transmision</italic></source>
<publisher-loc>Pereira, Colombia</publisher-loc>
<publisher-name>Editorial UTP</publisher-name>
<year>2016</year>
</element-citation>
</ref>
<ref id="ref-17-3635">
<mixed-citation>L. M. Londoño, R. A. Hincapié R. A. y Gallego,
“Planeamiento de redes de baja tensión, utilizando
un modelo trifásico,” Ciencia e Ingenieria Neogranadina,
vol. 21, no.2, pp. 41-56, 2011. DOI: https://doi.
org/10.18359/rcin.259</mixed-citation>
<element-citation publication-type="journal">
<person-group person-group-type="author">
<name>
     <surname>Londo&#x00F1;o</surname>
    <given-names>L. M.</given-names>
    </name>
<name>
     <surname>Hincapi&#x00E9;</surname>
    <given-names>R. A.</given-names>
    </name>
<name>
     <surname>Gallego</surname>
    <given-names>R. A.</given-names>
    </name>
    </person-group>
<article-title>Planeamiento de redes de baja tensi&#x00F3;n, utilizando un modelo trif&#x00E1;sico,</article-title>
<source xml:lang="es"><italic>Ciencia e Ingenieria Neogranadina</italic></source>
<volume>21</volume><issue>2</issue>
<fpage>41</fpage>
<lpage>56</lpage>
<year>2011</year>
<pub-id pub-id-type="doi">10.18359/rcin.259</pub-id>
</element-citation>
</ref>
</ref-list>
<sec sec-type="sec-16-3635">
<title><bold>Nomenclatura</bold></title>
<sec sec-type="sec-17-3635">
<title><bold>Conjuntos</bold></title>
<p><graphic xlink:href="1909-7735-cein-29-02-115-130-en1.png"/>: Conjunto de todos los tramos de red del sistema.</p>
<p><graphic xlink:href="1909-7735-cein-29-02-115-130-en2.png"/>: Conjunto de todos los nodos del sistema.</p>
<p><graphic xlink:href="1909-7735-cein-29-02-115-130-en3.png"/>: Conjunto de los nodos con almacenadores de energ&#x00ED;a.</p>
<p><graphic xlink:href="1909-7735-cein-29-02-115-130-en4.png"/>: Conjunto de los nodos con generaci&#x00F3;n.</p>
<p><graphic xlink:href="1909-7735-cein-29-02-115-130-en5.png"/>: Conjunto de los intervalos considerados en el periodo de estudio.</p>
</sec>
<sec sec-type="sec-18-3635">
<title><bold>&#x00CD;ndices y par&#x00E1;metros</bold></title>
<p><italic>B</italic>(<italic>i</italic>, <italic>j</italic>) : Susceptancia nodal entre los nodos <italic>i-j</italic>.</p>
<p><italic>b</italic>(<italic>i</italic>, <italic>j</italic>) : Susceptancia primitiva entre los nodos <italic>i-j</italic>.</p>
<p><italic>C<sub>x</sub></italic> (<italic>i</italic>, <italic>j</italic>): Costo de compra de la energ&#x00ED;a al generador <italic>g</italic> del nodo <italic>i</italic> en el periodo <italic>t</italic>.</p>
<p><graphic xlink:href="1909-7735-cein-29-02-115-130-en6.png"/>: Flujo de potencia activa m&#x00E1;xima permitida del tramo de red <italic>i-j</italic>.</p>
<p><graphic xlink:href="1909-7735-cein-29-02-115-130-en7.png"/>: Flujo de potencia reactiva m&#x00E1;xima permitida del tramo de red <italic>i-j</italic>.</p>
<p><italic>G</italic>(<italic>i</italic>, <italic>j</italic>) : Conductancia nodal entre los nodos <italic>i-j</italic>.</p>
<p><italic>g</italic>(<italic>i</italic>, <italic>j</italic>) : Conductancia primitiva entre los nodos <italic>i-j</italic>.</p>
<p><italic>i</italic>: &#x00CD;ndice que recorre los nodos.</p>
<p><italic>ij</italic>: &#x00CD;ndice que recorre los tramos.</p>
<p><italic>P<sub>carg ab</sub></italic>: Potencia activa de carga m&#x00E1;xima para el almacenador de energ&#x00ED;a <italic>b</italic> del nodo <italic>i</italic>.</p>
<p><graphic xlink:href="1909-7735-cein-29-02-115-130-en8.png"/>: Demanda de potencia activa del nodo <italic>i</italic> en el periodo <italic>t</italic>.</p>
<p><italic>P<sub>descarg ab</sub></italic>: Potencia activa de descarga m&#x00E1;xima para el almacenador de energ&#x00ED;a <italic>b</italic> del nodo <italic>i</italic>.</p>
<p><graphic xlink:href="1909-7735-cein-29-02-115-130-en9.png"/>: M&#x00E1;xima generaci&#x00F3;n de potencia activa permitida para el generador <italic>g</italic> del nodo <italic>i</italic>.</p>
<p><graphic xlink:href="1909-7735-cein-29-02-115-130-en10.png"/>: M&#x00ED;nima generaci&#x00F3;n de potencia activa permitida para el generador <italic>g</italic> del nodo <italic>i</italic>.</p>
<p><italic>Q<sub>d</sub> (t, j)</italic>: Demanda de potencia reactiva en el nodo <italic>i</italic> en el tiempo <italic>t</italic>.</p>
<p><graphic xlink:href="1909-7735-cein-29-02-115-130-en11.png"/>: M&#x00E1;xima generaci&#x00F3;n de potencia reactiva permitida para el generador <italic>g</italic> del nodo <italic>i</italic>.</p>
<p><graphic xlink:href="1909-7735-cein-29-02-115-130-en12.png"/>: M&#x00ED;nima generaci&#x00F3;n de potencia reactiva permitida para el generador <italic>g</italic> del nodo <italic>i</italic>.</p>
<p><graphic xlink:href="1909-7735-cein-29-02-115-130-en13.png"/>: Estado de carga m&#x00E1;xima permitida para el almacenador de energ&#x00ED;a <italic>b</italic> del nodo <italic>i</italic>.</p>
<p><graphic xlink:href="1909-7735-cein-29-02-115-130-en14.png"/>: Estado de carga m&#x00ED;nima permitida para el almacenador de energ&#x00ED;a <italic>b</italic> del nodo <italic>i</italic>.</p>
<p><graphic xlink:href="1909-7735-cein-29-02-115-130-en15.png"/>: Estado de carga inicial del almacenador de energ&#x00ED;a <italic>b</italic> del nodo <italic>i</italic>.</p>
<p><graphic xlink:href="1909-7735-cein-29-02-115-130-en16.png"/>: Estado de carga final del almacenador de energ&#x00ED;a <italic>b</italic> del nodo <italic>i</italic>.</p>
<p><italic>T</italic>: Periodo de estudio.</p>
<p><italic>t</italic>: &#x00CD;ndice que recorre cada intervalo.</p>
<p><italic>V&#x00AF;</italic>: Tensi&#x00F3;n m&#x00E1;xima del sistema.</p>
<p><italic>V&#x005F;</italic>: Tensi&#x00F3;n m&#x00ED;nima del sistema.</p>
<p><graphic xlink:href="1909-7735-cein-29-02-115-130-en17.png"/>: Intervalo entre cada periodo <italic>t</italic>.</p>
<p>f<sub>b</sub> : Pendiente de carga/descarga del almacenador de energ&#x00ED;a <italic>b</italic> del nodo <italic>i</italic>.</p>
</sec>
<sec sec-type="sec-19-3635">
<title><bold>Variables</bold></title>
<p><italic>F<sub>PL</sub> (t, ij)</italic>: Flujo de potencia activa en el tramo de red <italic>i-j</italic> en el periodo <italic>t</italic>.</p>
<p><italic>F<sub>Qt</sub> (t, ij)</italic>: Flujo de potencia reactiva en el tramo de red <italic>i-j</italic> en el periodo <italic>t</italic>.</p>
<p><italic>P<sub>b</sub> (t, i)</italic>: Potencia activa entregada/consumida por el almacenador de energ&#x00ED;a <italic>b</italic> del nodo <italic>i</italic> en el tiempo <italic>t</italic>.</p>
<p><italic>P<sub>b</sub> (t, i)</italic>: Potencia activa generada por el generador <italic>g</italic> del nodo <italic>i</italic> en el tiempo <italic>t</italic>.</p>
<p><italic>Q<sub>g</sub> (t, i)</italic>: Potencia reactiva generada por el generador <italic>g</italic> del nodo <italic>i</italic> en el tiempo <italic>t</italic>.</p>
<p><italic>SoC<sub>b</sub> (t, i)</italic>: Estado de carga del almacenador de energ&#x00ED;a <italic>b</italic> del nodo <italic>i</italic> en el tiempo <italic>t</italic>.</p>
<p><italic>V</italic>(<italic>t</italic>,<italic>i</italic>) : Magnitud de tensi&#x00F3;n del nodo <italic>i</italic> en el tiempo <italic>t</italic>.</p>
<p><italic>&#x03B8;<sub>i</sub></italic>: &#x00C1;ngulo del nodo <italic>i</italic>.</p>
<p><italic>&#x03B8;<sub>j</sub></italic>: &#x00C1;ngulo del nodo <italic>j</italic>.</p>
</sec>
</sec>
</back>
</article>
