<?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">WJET</journal-id><journal-title-group><journal-title>World Journal of Engineering and Technology</journal-title></journal-title-group><issn pub-type="epub">2331-4222</issn><publisher><publisher-name>Scientific Research Publishing</publisher-name></publisher></journal-meta><article-meta><article-id pub-id-type="doi">10.4236/wjet.2016.43D007</article-id><article-id pub-id-type="publisher-id">WJET-71324</article-id><article-categories><subj-group subj-group-type="heading"><subject>Articles</subject></subj-group><subj-group subj-group-type="Discipline-v2"><subject>Chemistry&amp;Materials Science</subject><subject> Engineering</subject></subj-group></article-categories><title-group><article-title>
 
 
  Reliability of Power System Integrated with Wind Generation Considering Carbon Tax
 
</article-title></title-group><contrib-group><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Xinwen</surname><given-names>Ni</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Kwok</surname><given-names>Lun Lo</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib></contrib-group><aff id="aff1"><addr-line>Power System Research Group, Department of Electronic and Electrical Engineering, Glasgow, UK</addr-line></aff><pub-date pub-type="epub"><day>20</day><month>10</month><year>2016</year></pub-date><volume>04</volume><issue>03</issue><fpage>50</fpage><lpage>57</lpage><history><date date-type="received"><day>July</day>	<month>13,</month>	<year>2016</year></date><date date-type="rev-recd"><day>Accepted:</day>	<month>October</month>	<year>13,</year>	</date><date date-type="accepted"><day>October</day>	<month>20,</month>	<year>2016</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>
 
 
   
   This paper proposes a method to evaluate the reliability of power system with different capacities of wind power while considering carbon tax. The proposed method is a hybrid approach which combines Frequency and Duration (F&amp;D) method and Monte Carlo Simulation (MCS) method. MCS method is used to achieve a model to simulate the random status of power system. Also, the proposed method is applied on the IEEE 14-bus test system to investigate the effects of integrating different capacities of wind energy to the reliability of power system with considering carbon tax. 
  
 
</p></abstract><kwd-group><kwd>Wind Power Integration</kwd><kwd> Power System Reliability</kwd><kwd> Carbon Tax</kwd></kwd-group></article-meta></front><body><sec id="s1"><title>1. Introduction</title><p>Renewable energy replacing traditional energy sources is an inevitable fact as more and more people have come to realize that the amount of easily obtainable fossil fuel is limited and keeps decreasing as well as the severe influences of by product of burning fossil fuels. Many countries had announced plans to replace traditional energy sources partially or entirely with renewable energy. According to UK government, 15% of the total energy demand will be generated using renewable sources by 2020 [<xref ref-type="bibr" rid="scirp.71324-ref1">1</xref>]. However, due to the intermittent and inconsistent nature of the wind, implementation of wind power into the current power grid could severely affect its performance. Moreover, integration of wind energy requires a significant amount of capital investment for network reinforcement and ancillary services. Approximately 75% of the total cost of energy for a wind turbine is related to upfront costs such as the expense of the turbine, foundation, electrical equipment, grid-connection and so on [<xref ref-type="bibr" rid="scirp.71324-ref3">3</xref>]. Also, implementation of wind energy helps to reduce the amount of CO<sub>2</sub> emission which leads to less expense of carbon tax.</p><p>The primary function of the power system is to provide electrical power to its customers as economically as possible with an acceptable degree of quality [<xref ref-type="bibr" rid="scirp.71324-ref2">2</xref>]. In other words, providing cheap electricity to consumers while keeping a certain level of reliability.</p><p>Reliability has always been and still is one of the most important factors when operating and maintaining a power system [<xref ref-type="bibr" rid="scirp.71324-ref8">8</xref>]. System reliability evaluation methods can be categorized into two groups, deterministic and probabilistic methods. Probabilistic methods are more widely used as it reflects the operating process of the power system in a better way. There are two aspects of probabilistic methods, analytical methods and Monte Carlo Simulation (MCS) methods. Analytical methods like Markov Chain Method and Frequency and Duration (F &amp; D) Method utilize mathematical models to represent power system and calculate different reliability indices through equations. Analytical methods can provide accurate results of power system reliability directly once given input data. The process can be very quick as long as the test system is simple. However, when it comes to large systems that contain lots of generating units and loads with various working status, the computation speed can slow down significantly, and the memory required to process and store data also increases considerably. MCS methods, on the other hand, represent the entire system with software modules and simulate the operating process of the power system while considering random behaviors of all participants. MCS method is based on repeating random sampling simulations to obtain results that as close to the real one as possible. Thus it can reflect the intermittent nature of wind power. This paper proposes a new method that combines MCS with F&amp;D while taking carbon tax into consideration. And this method is tested using modified IEEE 14-bus system in MATPOWER [<xref ref-type="bibr" rid="scirp.71324-ref4">4</xref>], a package developed within MATLAB.</p></sec><sec id="s2"><title>2. Simulation Process</title><p>Loss of Load Expectation (LOLE) is one of the most widely used probabilistic methods in power system reliability analysis. It is defined as the average hours/days per year that load demand exceeds power supply. F&amp;D method is an extension of LOLE as it also provides information about the expected frequency and duration of inadequacy. However, it cannot reflect the severeness of each encounter by showing the amount of power shortage. Loss of Energy Expectation measures the expected amount of electricity not to be met by generation in one year. Combining these two approaches together will give a more accurate method that provides all of the information. Nevertheless, the computation burden increases significantly when taking wind power into account due to the characteristics of wind energy. On the contrary, MCS method can be applied to simulate the random behaviors of the system without considering the generating unit’s output nature [<xref ref-type="bibr" rid="scirp.71324-ref5">5</xref>].</p><sec id="s2_1"><title>2.1. Test System</title><p>The test system used in this paper is modified based on the IEEE 14-bus test system as shown in <xref ref-type="fig" rid="fig1">Figure 1</xref>. This test system contains five generators and 11 loads. The maximum power output is 772.4 MW and the total amount of electricity demand is 259 MW. And a wind farm is connected to bus 3. The system is tested mainly in two scenarios, with and without considering line constraints.</p></sec><sec id="s2_2"><title>2.2. Wind Power Generation Model</title><p>Wind power employed in this paper is generated using Wei bull distribution function. By setting the scale factor α to 7 and shape factor β to 2, wind speed can be produced as <xref ref-type="fig" rid="fig2">Figure 2</xref> shows. And for the sake of simplicity, only partial wind speed data is presented in the illustration. <xref ref-type="fig" rid="fig3">Figure 3</xref> illustrates the sorted wind speed which fits the Wei</p><fig id="fig1"  position="float"><label><xref ref-type="fig" rid="fig1">Figure 1</xref></label><caption><title> IEEE 14-bus test system [<xref ref-type="bibr" rid="scirp.71324-ref6">6</xref>]</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/71324x2.png"/></fig><fig id="fig2"  position="float"><label><xref ref-type="fig" rid="fig2">Figure 2</xref></label><caption><title> Generated wind speed</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/71324x3.png"/></fig><fig id="fig3"  position="float"><label><xref ref-type="fig" rid="fig3">Figure 3</xref></label><caption><title> Sorted wind speed</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/71324x4.png"/></fig><p>bull distribution roughly. Then wind power output can be calculated using Equation (1). In this paper, cut-in speed is set to 4 m/s, rated speed is 10 m/s, cut-out speed is 25 m/s and rated power is 30 MW. And the price of wind energy is set to zero since it is considered as a must-taken source of energy.</p><disp-formula id="scirp.71324-formula127"><label>(1)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/71324x5.png"  xlink:type="simple"/></disp-formula><p>where:</p><p><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/71324x6.png" xlink:type="simple"/></inline-formula>is the wind speed</p><p><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/71324x7.png" xlink:type="simple"/></inline-formula>is the cut-in speed of wind turbine</p><p><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/71324x8.png" xlink:type="simple"/></inline-formula>is the rated speed of wind turbine</p><p><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/71324x9.png" xlink:type="simple"/></inline-formula>is the cut-out speed of wind turbine</p><p><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/71324x10.png" xlink:type="simple"/></inline-formula>is the rated power output of the wind turbine</p><p>A, B and C are constants calculated using cut-in and rated speed.</p></sec></sec><sec id="s3"><title>3. Case Analysis</title><p>A few cases are carried out to study the effect of wind power integration on power system while considering the cost of carbon tax. Assuming that [<xref ref-type="bibr" rid="scirp.71324-ref7">7</xref>] 6.89551 &#215; 10<sup>−4</sup> metric tons of CO<sub>2</sub> is generated to produce one kWh electricity by burning coal. And the carbon tax is 0.028 dollar per kWh.</p><sec id="s3_1"><title>3.1. Without Considering Line Constraints</title><p>In this section, the system is tested when connected with 5% and 10% wind power without taking transmission line constraints into account. <xref ref-type="fig" rid="fig4">Figure 4</xref> and <xref ref-type="fig" rid="fig5">Figure 5</xref> shows</p><fig id="fig4"  position="float"><label><xref ref-type="fig" rid="fig4">Figure 4</xref></label><caption><title> System operating cost,5% and 10% wind power, no constratints</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/71324x11.png"/></fig><fig id="fig5"  position="float"><label><xref ref-type="fig" rid="fig5">Figure 5</xref></label><caption><title> Carbon cost of the test system, 5% and 10% wind power, no constratints</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/71324x12.png"/></fig><p>the system operating cost and carbon cost respectively. The black line indicates the base case without connecting to wind farms while the red and green broken line represents system connected to 30MW and 60 MW wind farm. <xref ref-type="fig" rid="fig4">Figure 4</xref> and <xref ref-type="fig" rid="fig5">Figure 5</xref> looks very similar to operating cost and carbon cost are both calculated based on the total power generation. <xref ref-type="fig" rid="fig4">Figure 4</xref> shows that after connecting wind farm to bus, the operating cost is reduced significantly. And the system connected with 60 MW wind power has larger reduction than 30 MW. The maxim reduction of operating cost of 30 MW and 60 MW penetration level are 1178$/hr and 2316$/hr respectively. <xref ref-type="fig" rid="fig5">Figure 5</xref> indicates that the system has larger wind power capacity tend to have less expense of carbon cost. However, in <xref ref-type="fig" rid="fig4">Figure 4</xref>, there are some points that the operating costs are higher than before. This happens when the wind output is low and cannot provide enough electric power to consumers and generators with much greater price have to be brought online. <xref ref-type="table" rid="table1">Table 1</xref> showcases one extreme situation where the wind output is zero.</p></sec><sec id="s3_2"><title>3.2. Considering Line Constraints</title><p>In this case, a 30 MW maxim limit is introduced to all of the transmission lines. The blue and red line indicate 5% and 10% wind power penetration level respectively. <xref ref-type="fig" rid="fig6">Figure 6</xref> and <xref ref-type="fig" rid="fig7">Figure 7</xref> have very similar look comparing to <xref ref-type="fig" rid="fig4">Figure 4</xref> and <xref ref-type="fig" rid="fig5">Figure 5</xref>. The operating cost is reduced and the system with higher wind power penetration level tend to have a larger price drop. As <xref ref-type="table" rid="table2">Table 2</xref> shows, when transmission line limit is considered, system operating cost has smaller decrease comparing to it without taking constraint into account. This happens when wind farm has abundant output. After meeting the local demand, the extra electricity is prevented being distributed to other buses by</p><table-wrap id="table1" ><label><xref ref-type="table" rid="table1">Table 1</xref></label><caption><title> Bus data, unconstratined, 30 MW wind</title></caption><table><tbody><thead><tr><th align="center" valign="middle" >Bus No.</th><th align="center" valign="middle" >Generation (MW)</th><th align="center" valign="middle" >Load (MW)</th><th align="center" valign="middle" >LMP ($/MWh)</th></tr></thead><tr><td align="center" valign="middle" >1</td><td align="center" valign="middle" >195.44</td><td align="center" valign="middle" >-</td><td align="center" valign="middle" >36.819</td></tr><tr><td align="center" valign="middle" >2</td><td align="center" valign="middle" >36.37</td><td align="center" valign="middle" >21.70</td><td align="center" valign="middle" >38.382</td></tr><tr><td align="center" valign="middle" >3</td><td align="center" valign="middle" >0</td><td align="center" valign="middle" >24.20</td><td align="center" valign="middle" >39.818</td></tr><tr><td align="center" valign="middle" >4</td><td align="center" valign="middle" >-</td><td align="center" valign="middle" >67.80</td><td align="center" valign="middle" >40.574</td></tr><tr><td align="center" valign="middle" >5</td><td align="center" valign="middle" >-</td><td align="center" valign="middle" >17.60</td><td align="center" valign="middle" >40.085</td></tr><tr><td align="center" valign="middle" >6</td><td align="center" valign="middle" >10.33</td><td align="center" valign="middle" >11.20</td><td align="center" valign="middle" >40.206</td></tr><tr><td align="center" valign="middle" >7</td><td align="center" valign="middle" >-</td><td align="center" valign="middle" >-</td><td align="center" valign="middle" >40.546</td></tr><tr><td align="center" valign="middle" >8</td><td align="center" valign="middle" >26.83</td><td align="center" valign="middle" >-</td><td align="center" valign="middle" >40.537</td></tr><tr><td align="center" valign="middle" >9</td><td align="center" valign="middle" >-</td><td align="center" valign="middle" >29.50</td><td align="center" valign="middle" >40.550</td></tr><tr><td align="center" valign="middle" >10</td><td align="center" valign="middle" >-</td><td align="center" valign="middle" >19.00</td><td align="center" valign="middle" >41.048</td></tr><tr><td align="center" valign="middle" >11</td><td align="center" valign="middle" >-</td><td align="center" valign="middle" >13.50</td><td align="center" valign="middle" >41.106</td></tr><tr><td align="center" valign="middle" >12</td><td align="center" valign="middle" >-</td><td align="center" valign="middle" >16.10</td><td align="center" valign="middle" >41.619</td></tr><tr><td align="center" valign="middle" >13</td><td align="center" valign="middle" >-</td><td align="center" valign="middle" >13.50</td><td align="center" valign="middle" >41.330</td></tr><tr><td align="center" valign="middle" >14</td><td align="center" valign="middle" >-</td><td align="center" valign="middle" >24.90</td><td align="center" valign="middle" >42.392</td></tr><tr><td align="center" valign="middle" >Operating cost</td><td align="center" valign="middle" >8119.82 ($/hr)</td><td align="center" valign="middle" >Carbon tax</td><td align="center" valign="middle" >17.7750$/hr</td></tr><tr><td align="center" valign="middle" >LOLE</td><td align="center" valign="middle" >-</td><td align="center" valign="middle" >LOEE</td><td align="center" valign="middle" >-</td></tr></tbody></table></table-wrap><fig id="fig6"  position="float"><label><xref ref-type="fig" rid="fig6">Figure 6</xref></label><caption><title> System operating cost connecting 5% and 10% wind power, considering constratints</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/71324x13.png"/></fig><fig id="fig7"  position="float"><label><xref ref-type="fig" rid="fig7">Figure 7</xref></label><caption><title> System operating cost connecting 5% and 10% wind power, considering constratints</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/71324x14.png"/></fig><p>the transmission line limit. The LOLE and LOEE when 30 MW of wind power connected is 4.2 hours per week and 15.39 MW respectively. When 60 MW of wind power is integrated into the system, the indices increases to 5.56 and 18.11.</p><table-wrap id="table2" ><label><xref ref-type="table" rid="table2">Table 2</xref></label><caption><title> Camparision of operating cost and carbon cost</title></caption><table><tbody><thead><tr><th align="center" valign="middle" ></th><th align="center" valign="middle" >Max operating cost</th><th align="center" valign="middle" >Min operating Cost</th><th align="center" valign="middle" >Max carbon cost</th><th align="center" valign="middle" >Min carbon cost</th></tr></thead><tr><td align="center" valign="middle" >Base case</td><td align="center" valign="middle" >8108</td><td align="center" valign="middle" >8108</td><td align="center" valign="middle" >17.87</td><td align="center" valign="middle" >17.87</td></tr><tr><td align="center" valign="middle" >Unconstrained 5% wind power</td><td align="center" valign="middle" >8119</td><td align="center" valign="middle" >6941</td><td align="center" valign="middle" >17.77</td><td align="center" valign="middle" >16.83</td></tr><tr><td align="center" valign="middle" >Unconstrained 10% wind power</td><td align="center" valign="middle" >8119</td><td align="center" valign="middle" >5803</td><td align="center" valign="middle" >17.77</td><td align="center" valign="middle" >15.58</td></tr><tr><td align="center" valign="middle" >Constrained 5% wind power</td><td align="center" valign="middle" >8119</td><td align="center" valign="middle" >6975</td><td align="center" valign="middle" >17.77</td><td align="center" valign="middle" >16.28</td></tr><tr><td align="center" valign="middle" >Constrained 10% wind power</td><td align="center" valign="middle" >8119</td><td align="center" valign="middle" >6934</td><td align="center" valign="middle" >17.77</td><td align="center" valign="middle" >15.24</td></tr></tbody></table></table-wrap></sec></sec><sec id="s4"><title>4. Conclusion</title><p>This paper presents a way to evaluate system reliability while considering carbon tax. It shows that wind power integration help to reduce the operating cost and carbon cost by cutting the emission of CO<sub>2</sub> though it requires a large amount of capital investment. Power system reliability, on the other hand, is also affected by integrating wind power.</p></sec><sec id="s5"><title>Cite this paper</title><p>Ni, X.W. and Lo, K.L. (2016) Reliability of Power System In- tegrated with Wind Generation Considering Carbon Tax. World Journal of Engineering and Technology, 4, 50-57. http://dx.doi.org/10.4236/wjet.2016.43D007</p></sec></body><back><ref-list><title>References</title><ref id="scirp.71324-ref1"><label>1</label><mixed-citation publication-type="other" xlink:type="simple">Department of Energy and Climate Change (2011) UK Renewable Energy Road Map.  
https://www.gov.uk/government/uploads/system/uploads/attachment_data/file/48128/2167-uk-renewable-energy-roadmap.pdf</mixed-citation></ref><ref id="scirp.71324-ref2"><label>2</label><mixed-citation publication-type="other" xlink:type="simple">Li, W.Y. and Billinton, R. Reliability Assessment of Electric Power Systems Using Monte Carlo Methods.</mixed-citation></ref><ref id="scirp.71324-ref3"><label>3</label><mixed-citation publication-type="other" xlink:type="simple">The European Wind Energy Associa-tion. The Economics of Wind Energy.  
http://www.ewea.org/fileadmin/files/library/publications/reports/Economics_of_Wind_Energy.pdf</mixed-citation></ref><ref id="scirp.71324-ref4"><label>4</label><mixed-citation publication-type="other" xlink:type="simple">Zimmerman, R.D., Murillo-Sánchez, C.E. and Thomas, R.J. (2011) MATPOWER: Steady- State Operations, Planning and Analysis Tools for Power Systems Research and Education. IEEE Transactions on Power Systems, 26, 12-19.</mixed-citation></ref><ref id="scirp.71324-ref5"><label>5</label><mixed-citation publication-type="other" xlink:type="simple">Shi, S. and Lo, K.L. Reliability Assessment of Power System Considering the Impact of Wind Energy.</mixed-citation></ref><ref id="scirp.71324-ref6"><label>6</label><mixed-citation publication-type="other" xlink:type="simple">https://www.ee.washington.edu/resear/pstca/pf14/pg_tcal14bus.thm</mixed-citation></ref><ref id="scirp.71324-ref7"><label>7</label><mixed-citation publication-type="other" xlink:type="simple">US Department of Energy, Energy Information Administration. Fuel and Energy Source Codes and Emission Coefficient.</mixed-citation></ref><ref id="scirp.71324-ref8"><label>8</label><mixed-citation publication-type="other" xlink:type="simple">Nikzad, M., Shams, S., Farahani, S. and Tabar, M.B. (2012) Calculation of Generation System Reliability Index: Loss of Load Expectation. Life Science Journal, 9.</mixed-citation></ref></ref-list></back></article>