<?xml version="1.0" encoding="UTF-8"?><!DOCTYPE article  PUBLIC "-//NLM//DTD Journal Publishing DTD v3.0 20080202//EN" "http://dtd.nlm.nih.gov/publishing/3.0/journalpublishing3.dtd"><article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" dtd-version="3.0" xml:lang="en" article-type="research article"><front><journal-meta><journal-id journal-id-type="publisher-id">ENG</journal-id><journal-title-group><journal-title>Engineering</journal-title></journal-title-group><issn pub-type="epub">1947-3931</issn><publisher><publisher-name>Scientific Research Publishing</publisher-name></publisher></journal-meta><article-meta><article-id pub-id-type="doi">10.4236/eng.2019.114017</article-id><article-id pub-id-type="publisher-id">ENG-91725</article-id><article-categories><subj-group subj-group-type="heading"><subject>Articles</subject></subj-group><subj-group subj-group-type="Discipline-v2"><subject>Engineering</subject></subj-group></article-categories><title-group><article-title>
 
 
  Taking into Account of Functional Constraints in Optimization of Modes of Power Systems by Genetic Algorithms
 
</article-title></title-group><contrib-group><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Gayibov</surname><given-names>Tulkin Shernazarovich</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref><xref ref-type="corresp" rid="cor1"><sup>*</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Latipov</surname><given-names>Sherkhon Shuxratovich</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib></contrib-group><aff id="aff1"><addr-line>Department of Power Plants, Networks and Systems, Tashkent State Technical University, Tashkent, Uzbekistan</addr-line></aff><pub-date pub-type="epub"><day>11</day><month>04</month><year>2019</year></pub-date><volume>11</volume><issue>04</issue><fpage>240</fpage><lpage>246</lpage><history><date date-type="received"><day>5,</day>	<month>March</month>	<year>2019</year></date><date date-type="rev-recd"><day>8,</day>	<month>April</month>	<year>2019</year>	</date><date date-type="accepted"><day>11,</day>	<month>April</month>	<year>2019</year></date></history><permissions><copyright-statement>&#169; Copyright  2014 by authors and Scientific Research Publishing Inc. </copyright-statement><copyright-year>2014</copyright-year><license><license-p>This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/</license-p></license></permissions><abstract><p>
 
 
  The development of the capabilities of computational tools has created up new possibilities for the effective use of a number of classical mathematical methods and algorithms for solving many important problems in the power engineering. In particular, a set of algorithms are developed to optimize the modes of electric power systems based on genetic algorithms. At the same time, the issues of taking into account functional constraints in solving such problems by genetic algorithms need to be improved. In accordance with it in this article the problems of taking into account of different constraints in optimization of modes of power systems using genetic algorithms are considered. The algorithm of optimization by genetic algorithm taking into account of functional constraints in forms of equality and inequality by penalty functions is proposed. The results of research of proposed algorithm’s efficiency in example of optimization of mode of power system with 8 buses, 4 thermal power plants and 3 transmission lines with controlled power flow are presented.
 
</p></abstract><kwd-group><kwd>Optimization</kwd><kwd> Power System</kwd><kwd> Power Plant</kwd><kwd> Objective Function</kwd><kwd> Constraint</kwd></kwd-group></article-meta></front><body><sec id="s1"><title>1. Introduction</title><p>Optimization of modes of modern power systems on active power is a complex task of nonlinear mathematical programming with a set of simple and functional constraints in the forms of equality and inequality. Therefore, the effectiveness of algorithms of its solution is determined, in particular, by the possibility of effectively taking into account such constraints.</p><p>The main constraints that are taken into account in solving the problem of optimal planning of modes of power systems with thermal power plants for a certain period of time are:</p><p>- constraints in the form of inequality on minimum and maximum possible loads of power plants (simple constraints)</p><p>P i min ≤ P i ≤ P i max ,   i = 1 ,   2 ,   ⋯ ,   n ; (1)</p><p>- condition of active power balance in power system (functional constraint in the form of equality)</p><p>∑ i = 1 n P i = P L , (2)</p><p>where n is the number of power plants, which participate in optimization; P<sub>i</sub>―power of ith plant; P<sub>L</sub>―the total load of consumers in power system;</p><p>- constraints on power flows along certain controlled transmission lines (functional constraints in the form of inequality)</p><p>P l min ≤ P l ≤ P l max ,   l = 1 ,   2 ,   ⋯ ,   L , (3)</p><p>where L is the number of lines in which active power flows are limited; P<sub>l</sub>-active power flow along the lth controlled line.</p><p>At present, the methods for taking into account such constraints in optimization of modes of power systems by traditional algorithms are quite developed [<xref ref-type="bibr" rid="scirp.91725-ref1">1</xref>] [<xref ref-type="bibr" rid="scirp.91725-ref2">2</xref>] [<xref ref-type="bibr" rid="scirp.91725-ref3">3</xref>] [<xref ref-type="bibr" rid="scirp.91725-ref4">4</xref>] . However, the issues of taking into account functional constraints in solving of optimization problems using nontraditional algorithms, in particular, genetic algorithms [<xref ref-type="bibr" rid="scirp.91725-ref5">5</xref>] [<xref ref-type="bibr" rid="scirp.91725-ref6">6</xref>] [<xref ref-type="bibr" rid="scirp.91725-ref7">7</xref>] [<xref ref-type="bibr" rid="scirp.91725-ref8">8</xref>] , require additional research.</p><p>Genetic algorithms offer a new and powerful approach to solving optimization problems. Their use was made possible by expanding the capabilities of computational tools at relatively low costs. Recently, these algorithms are used in solving global problems of search optimization, when traditional optimization algorithms cannot be used. They use parallel and global search methods that mimic natural genetic operators. The probability of convergence of a genetic algorithm to the global solution of the problem is the highest, since it simultaneously evaluates a set of points in the parameter space. These algorithms also do not require differentiability and continuity of the search space. At the same time, the issues of taking into account of functional constraints in the form of equality and inequality in optimization of power system modes by these algorithms are researched not sufficiently.</p></sec><sec id="s2"><title>2. The Algorithms of Taking into Account of Functional Constraints</title><p>The issues of using the genetic algorithms for optimal distribution of consumers load in power system between thermal power plants were researched by author of this paper in [<xref ref-type="bibr" rid="scirp.91725-ref9">9</xref>] In present paper, on the basis of researches carried out the algorithm of taking into account of functional constraints in the form of equality and inequality in optimization of modes of power system by genetic algorithms is proposed. According to this algorithm, the problem of conditional minimization of a function, which is a function of total fuel costs or total consumption of reference fuel B in thermal power plants of the power system, is reduced to unconditional minimization of a generalized function consisting of sum of the objective function of the initial problem and the penalty functions that take into account constraints (2) and (3). Simple constraints (1) are taken into account in accordance with the procedure of the genetic algorithm automatically.</p><p>The value of penalty function that take into account of the constraint should be equal to zero when it is satisfied and increase proportionally to the degree of violation when the constraint is not satisfied.</p><p>In accordance with it, in order to take into account the constraint in the form of equality (2), the penalty function is taken in the form of</p><p>S P = α ( ∑ i = 1 n P i − P L ) 2 . (4)</p><p>Particular difficult is the problem of taking into account of functional constraints in the form of inequality (3). In many traditional algorithms of power systems mode optimization such constraints are also taken into account by quadratic penalty functions as</p><p>S l = α ( P l − P l max ) 2 or S l = α ( − P l + P l min ) 2 . (5)</p><p>This form of the penalty function is adapted for taking into account of constraints in the form of equality. For example, to take into account of constraint P l ≤ P l max in the task of finding the minimum of an objective function, the penalty function has a positive value and grows even when the constraint is satisfied (<xref ref-type="fig" rid="fig1">Figure 1</xref>(a)).</p><p>To overcome this difficulty in traditional algorithms of optimization after each iteration of the optimization process, the constraint is checked out. In the case of its satisfying the next step of calculation is carried out without taking into account of this constraint. Otherwise, it is carried out with taking into account of the constraint using the penalty function.</p><p>The results of research show that taking into account of functional constraints in the form of inequality using such penalty functions in optimization by traditional</p><p>algorithms generally cause some difficulties connected with choosing the penalty coefficients and ensuring the convergence of the iterative calculating process, which is often oscillatory in nature [<xref ref-type="bibr" rid="scirp.91725-ref1">1</xref>] .</p><p>There are some difficulties for taking into account the constraints in the form of inequality (3) using penalty functions (5) in optimization by genetic algorithms. Since the use of the objective function and the corresponding fitness function, which do not include such constraints may lead to unacceptable results, where some of these constraints are not satisfied. Therefore, for overcome of such difficulties in this work the use of the exponential form of penalty function is recommended:</p><p>S l = α e β ( P l − P l max ) or S l = α e β ( − P l + P l min ) , (6)</p><p>where α, β are weight (penalty) coefficients.</p><p>This form of penalty function is very closely meets the requirements set for it―equality to zero when the constraint is satisfied and increase proportionally to degree of violation when it is not satisfied (<xref ref-type="fig" rid="fig1">Figure 1</xref>(b)).</p><p>By selecting the values of the penalty coefficients α and β, it is possible to take into account the constraints with sufficient accuracy.</p><p>Thus, the generalized objective function when optimized by the genetic algorithm taking into account of constraints is presented in the following form:</p><p>F = B + S P + ∑ l S l . (7)</p></sec><sec id="s3"><title>3. Example and Evaluation of Results</title><p>The efficiency of the described algorithm is researched in example of optimization of mode of power system with 8 buses and 13 power transmission lines (PTL) (<xref ref-type="fig" rid="fig2">Figure 2</xref>). In 4 buses there are thermal power plants (TPP) with the following fuel consumption characteristics of equivalent fuel, toe/h:</p><p>B 0 = 70 + 0.14 P 0 + 0.0014 P 0 2 , B 1 = 84 + 0.14 P 1 + 0.00175 P 1 2 ,</p><p>B 6 = 42 + 0.105 P 6 + 0.00105 P 6 2 , B 7 = 56 + 0.175 P 7 + 0.0007 P 7 2 .</p><p>The constraints on the minimum and maximum loads of TPPs are the same and P<sup>min</sup> = 200 MW, P<sup>m</sup><sup>ах</sup> = 700 MW. Capacities of consumers are: P<sub>2</sub> = 400 MW, P<sub>3</sub> = 600 MW, P<sub>4</sub> = 200 MW, P<sub>5</sub> = 500 MW.</p><p>Three power transmission lines have limited active power flows:</p><p>P<sub>6-3</sub> ≤ 450 MW, P<sub>6-5</sub> ≤ 65 MW, P<sub>0-3</sub> ≤ 95 MW.</p><p>Power flows for these controlled transmission lines are determined using distribution coefficients C<sub>li</sub>, given in <xref ref-type="table" rid="table1">Table 1</xref>, on following expression:</p><p>P l = ∑ i = 1 n C l i P i + P l 0 .</p><p>To compare the results of optimization on the proposed algorithm, <xref ref-type="table" rid="table2">Table 2</xref> shows the reference result obtained by the classical gradient method.</p><p>The results of optimization by the genetic algorithm, taking into account of</p><table-wrap id="table1" ><label><xref ref-type="table" rid="table1">Table 1</xref></label><caption><title> Distribution coefficients of buses power over controlled transmission lines</title></caption><table><tbody><thead><tr><th align="center" valign="middle"  rowspan="2"  >Power lines</th><th align="center" valign="middle"  colspan="7"  >Buses</th></tr></thead><tr><td align="center" valign="middle" >1</td><td align="center" valign="middle" >2</td><td align="center" valign="middle" >3</td><td align="center" valign="middle" >4</td><td align="center" valign="middle" >5</td><td align="center" valign="middle" >6</td><td align="center" valign="middle" >7</td></tr><tr><td align="center" valign="middle" >6 - 3</td><td align="center" valign="middle" >0.2536</td><td align="center" valign="middle" >−0.0713</td><td align="center" valign="middle" >−0.1720</td><td align="center" valign="middle" >0.4106</td><td align="center" valign="middle" >0.4890</td><td align="center" valign="middle" >0.6343</td><td align="center" valign="middle" >0.4986</td></tr><tr><td align="center" valign="middle" >6 - 5</td><td align="center" valign="middle" >−0.0701</td><td align="center" valign="middle" >0.02025</td><td align="center" valign="middle" >0.0484</td><td align="center" valign="middle" >−0.2810</td><td align="center" valign="middle" >−0.4223</td><td align="center" valign="middle" >0.1029</td><td align="center" valign="middle" >−0.2986</td></tr><tr><td align="center" valign="middle" >0 - 3</td><td align="center" valign="middle" >−0.1400</td><td align="center" valign="middle" >−0.1735</td><td align="center" valign="middle" >−0.2670</td><td align="center" valign="middle" >−0.1744</td><td align="center" valign="middle" >−0.1917</td><td align="center" valign="middle" >−0.2240</td><td align="center" valign="middle" >−0.1939</td></tr></tbody></table></table-wrap><table-wrap id="table2" ><label><xref ref-type="table" rid="table2">Table 2</xref></label><caption><title> Reference result of optimization</title></caption><table><tbody><thead><tr><th align="center" valign="middle" >TPP</th><th align="center" valign="middle" >Power of TPP, MW</th><th align="center" valign="middle" >Fuel consumption, toe/h.</th></tr></thead><tr><td align="center" valign="middle" >TPP at bus 0</td><td align="center" valign="middle" >234.7</td><td align="center" valign="middle" >179.97</td></tr><tr><td align="center" valign="middle" >TPP at bus 1</td><td align="center" valign="middle" >505.6</td><td align="center" valign="middle" >602.14</td></tr><tr><td align="center" valign="middle" >TPP at bus 6</td><td align="center" valign="middle" >281.3</td><td align="center" valign="middle" >154.55</td></tr><tr><td align="center" valign="middle" >TPP at bus 7</td><td align="center" valign="middle" >678.4</td><td align="center" valign="middle" >496.88</td></tr><tr><td align="center" valign="middle" >Total fuel consumption in power system</td><td align="center" valign="middle" >-</td><td align="center" valign="middle" >1433.54</td></tr><tr><td align="center" valign="middle" >Power flow on controlled lines</td><td align="center" valign="middle"  colspan="2"  >Р<sub>6-3</sub> = 450.0 MW, Р<sub>6-5</sub> = 21.2 MW, Р<sub>0-3</sub> = 95.0 MW.</td></tr></tbody></table></table-wrap><p>constraints using proposed algorithm at values of weighting coefficients α = 10, β = 2 are given in <xref ref-type="table" rid="table3">Table 3</xref>.</p><p>Comparing the obtained results of optimization by proposed algorithm with the reference result, we will verify the sufficiently high accuracy of the proposed algorithm.</p><p>Thus, the proposed algorithm of taking into account of functional constraints in the form of equality and inequality in optimization of power systems mode by genetic algorithms is highly efficient. It allows you to reliably determine the optimal mode of power system taking into account simple and functional constraints in the form of equality and inequality with acceptable accuracy for practical purposes.</p></sec><sec id="s4"><title>4. Conclusions</title><p>1) The efficient algorithm for taking into account of functional constraints in</p><table-wrap id="table3" ><label><xref ref-type="table" rid="table3">Table 3</xref></label><caption><title> Optimization results with taking into account of constraints by proposed algorithm</title></caption><table><tbody><thead><tr><th align="center" valign="middle" >TPP</th><th align="center" valign="middle" >Power of TPP, MW</th><th align="center" valign="middle" >Fuel consumption, toe/h.</th></tr></thead><tr><td align="center" valign="middle" >TPP at bus 0</td><td align="center" valign="middle" >236.9</td><td align="center" valign="middle" >181.74</td></tr><tr><td align="center" valign="middle" >TPP at bus 1</td><td align="center" valign="middle" >498.8</td><td align="center" valign="middle" >589.23</td></tr><tr><td align="center" valign="middle" >TPP at bus 6</td><td align="center" valign="middle" >273.7</td><td align="center" valign="middle" >149.39</td></tr><tr><td align="center" valign="middle" >TPP at bus 7</td><td align="center" valign="middle" >690.5</td><td align="center" valign="middle" >510.59</td></tr><tr><td align="center" valign="middle" >Total fuel consumption in power system</td><td align="center" valign="middle" >-</td><td align="center" valign="middle" >1430.96</td></tr><tr><td align="center" valign="middle" >Power flow on controlled lines</td><td align="center" valign="middle"  colspan="2"  >Р<sub>6-3</sub> = 449.5 MW, Р<sub>6-5</sub> = 17.2 MW, Р<sub>0-3</sub> = 95.3 MW.</td></tr></tbody></table></table-wrap><p>the form of equality and inequality in optimization of modes of power systems by genetic algorithm is proposed.</p><p>2) Genetic algorithm for optimization of modes of power systems taking into account the functional constraints in the form of inequalities by exponential form of penalty function has a reliable convergence of an iterative calculation process. It also does not require differentiability and continuity of the search space.</p><p>3) The proposed algorithm of taking into account of functional constraints in optimization by genetic algorithms can be effectively used for optimal planning of short-term modes of power systems.</p></sec><sec id="s5"><title>Conflicts of Interest</title><p>The authors declare no conflicts of interest regarding the publication of this paper.</p></sec><sec id="s6"><title>Cite this paper</title><p>Shernazarovich, G.T. and Shuxratovich, L.S. (2019) Taking into Account of Functional Constraints in Optimization of Modes of Power Systems by Genetic Algorithms. Engineering, 11, 240-246. https://doi.org/10.4236/eng.2019.114017</p></sec></body><back><ref-list><title>References</title><ref id="scirp.91725-ref1"><label>1</label><mixed-citation publication-type="other" xlink:type="simple">Nasirov, T.Kh. and Gayibov, T.Sh. (2014) Teoreticheskie osnovi optimizatsii rejimov energosystem [Theoretical Foundations of Optimization of Modes of Power Systems]. Fan va Texnologiya Publishing, Tashkent. (In Russian)</mixed-citation></ref><ref id="scirp.91725-ref2"><label>2</label><mixed-citation publication-type="other" xlink:type="simple">Rudenko, Yu.N. and Semyonov, V.A. (2000) Avtomatizatsiya dispetcherskogo upravleniya v elektroenergetike [Automation of Dispatching Control in Electric Power Industry]. MEI Publishing, Moscow. (In Russian)</mixed-citation></ref><ref id="scirp.91725-ref3"><label>3</label><mixed-citation publication-type="other" xlink:type="simple">Murty, P.S.R. (2008) Operation and Control in Power Systems. B.S. Publications, Hyderabad.</mixed-citation></ref><ref id="scirp.91725-ref4"><label>4</label><mixed-citation publication-type="other" xlink:type="simple">Tulkin, G. and Kamal, R. (2017) Optimal Planning of Short-Term Modes of Power Systems with Control of Loads of Electric Consumers Taking into Account of Network Factor. European Science Review, 9-10, 86-91.</mixed-citation></ref><ref id="scirp.91725-ref5"><label>5</label><mixed-citation publication-type="other" xlink:type="simple">Yuryevich, J. and Wong, K.P. (1999) Evolutionary Programming Based Optimal Power Flow Algorithm. IEEE Transaction on Power Systems, 14, 1245-1250.  
https://doi.org/10.1109/59.801880</mixed-citation></ref><ref id="scirp.91725-ref6"><label>6</label><mixed-citation publication-type="other" xlink:type="simple">Lai, L.L., Ma, J.T., Yokoma, R. and Zhao, M. (1997) Improved Genetic Algorithms for Optimal Power Flow under Both Normal and Contingent Operation States. Electrical Power &amp; Energy System, 19, 287-292.  
https://doi.org/10.1016/S0142-0615(96)00051-8</mixed-citation></ref><ref id="scirp.91725-ref7"><label>7</label><mixed-citation publication-type="other" xlink:type="simple">Amin, N., Marowan, A.Q. and Ghosh, B.C. (2017) Application of Genetic Algorithm in Power System Optimization with Multi-Type FACTS. International Journal of Research Publications, 7, 748-754.</mixed-citation></ref><ref id="scirp.91725-ref8"><label>8</label><mixed-citation publication-type="other" xlink:type="simple">Bhaskar, K., Vijay, M., William, A. and Bhattacharya, D. (2013) Mathematical Modeling of Optimizing Power Stream Measurement Using Genetic Algorithm. American Journal of Engineering Research (AJER), 2, 71-79.</mixed-citation></ref><ref id="scirp.91725-ref9"><label>9</label><mixed-citation publication-type="other" xlink:type="simple">Gayibov, T.Sh. (2017) Optimizatsiya rejimov energosistem geneticheskimi algiritmami [Optimization of Power Systems Modes by Genetic Algorithms]. Problems of Energy and Resource Saving, 1-2, 43-48.</mixed-citation></ref></ref-list></back></article>