发明名称 一种基于多目标优化的电网规划设计方法
摘要 本发明公开了一种基于多目标优化的电网规划设计方法,其综合考虑了经济性适应度指标、可靠性适应度指标和地理适应度指标,将可靠性指标转化成经济形式加入目标函数,使优化方案的综合效益达到最佳,而且可以对规划方案的经济性和可靠性进行灵活地评价,正确地反映投入资金对可靠性指标的增幅之间的确定关系,从而使电网规划的成本计算更为准确。
申请公布号 CN103049799B 申请公布日期 2016.07.06
申请号 CN201210526881.1 申请日期 2012.12.10
申请人 河海大学;南京河海电力软件有限公司 发明人 丁晓群;周春泉;周玲;吴桂军;史静;陈光宇
分类号 G06F17/15(2006.01)I 主分类号 G06F17/15(2006.01)I
代理机构 南京苏高专利商标事务所(普通合伙) 32204 代理人 柏尚春
主权项 一种基于多目标优化的电网规划设计方法,其特征在于,包括以下步骤:1)对规划区电网运行的历史数据,进行技术合理性、运行安全性和电网转供能力分析,作出电网现状评估;2)收集、分析规划区的经济、人口、政府的有关政策和城市总体发展规划的资料,利用空间电力负荷预测方法进行负荷预测;3)根据步骤2)得出的电力负荷总量,对规划区进行电力电量分析,确定规划区的电压等级序列、电网结构形式和设备选择原则,并且要与当地的地方电网规划原则相适应;4)将规划区按地块功能、地理环境和行政划分因素划分为若干个相互独立、多层次的供电区域;5)利用地理信息系统的数据库,在规划区内建立一个坐标系,确定各负荷中心的坐标;6)基于等负荷原则、初投资最小原则、负荷距最小原则和网络运行费用最小原则对变电站进行优化初步选址,并根据地理情况做相应的站址调整,选出候选变电站站址,并计算其相应的经济性适应度;多源连续选址就是在一个规划区中同时确定几个变电站的站址,在建立模型时,基于不同目标所建立的模型不同,这里设所研究的问题有m个变电站,向n个负荷点供电,则基于等负荷原则、初投资最小原则、负荷距最小原则和网络运行费用最小原则的多源连续选址的目标函数分别为:<maths num="0001"><math><![CDATA[<mrow><msub><mi>minf</mi><mn>1</mn></msub><mo>=</mo><munderover><mo>&Sigma;</mo><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msub><mi>&delta;</mi><mrow><mi>j</mi><mi>i</mi></mrow></msub><mo>&lsqb;</mo><mo>|</mo><msub><mi>u</mi><mi>j</mi></msub><mo>-</mo><msub><mi>x</mi><mi>i</mi></msub><mo>|</mo><mo>+</mo><mo>|</mo><msub><mi>v</mi><mi>j</mi></msub><mo>-</mo><msub><mi>y</mi><mi>i</mi></msub><mo>|</mo><mo>&rsqb;</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000960278230000011.GIF" wi="1277" he="151" /></maths><maths num="0002"><math><![CDATA[<mrow><msub><mi>minf</mi><mn>2</mn></msub><mo>=</mo><munderover><mo>&Sigma;</mo><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msub><mi>&delta;</mi><mrow><mi>j</mi><mi>i</mi></mrow></msub><msub><mi>A</mi><mi>i</mi></msub><mo>&lsqb;</mo><mo>|</mo><msub><mi>u</mi><mi>j</mi></msub><mo>-</mo><msub><mi>x</mi><mi>i</mi></msub><mo>|</mo><mo>+</mo><mo>|</mo><msub><mi>v</mi><mi>j</mi></msub><mo>-</mo><msub><mi>y</mi><mi>i</mi></msub><mo>|</mo><mo>&rsqb;</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000960278230000012.GIF" wi="1270" he="151" /></maths><maths num="0003"><math><![CDATA[<mrow><msub><mi>minf</mi><mn>3</mn></msub><mo>=</mo><munderover><mo>&Sigma;</mo><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msub><mi>&delta;</mi><mrow><mi>j</mi><mi>i</mi></mrow></msub><msub><mi>P</mi><mi>i</mi></msub><mo>&lsqb;</mo><mo>|</mo><msub><mi>u</mi><mi>j</mi></msub><mo>-</mo><msub><mi>x</mi><mi>i</mi></msub><mo>|</mo><mo>+</mo><mo>|</mo><msub><mi>v</mi><mi>j</mi></msub><mo>-</mo><msub><mi>y</mi><mi>i</mi></msub><mo>|</mo><mo>&rsqb;</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000960278230000013.GIF" wi="1286" he="150" /></maths><maths num="0004"><math><![CDATA[<mrow><msub><mi>minf</mi><mn>4</mn></msub><mo>=</mo><munderover><mo>&Sigma;</mo><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msub><mi>&delta;</mi><mrow><mi>j</mi><mi>i</mi></mrow></msub><msub><mi>&beta;</mi><mi>i</mi></msub><msub><mi>P</mi><mi>i</mi></msub><mo>&lsqb;</mo><mo>|</mo><msub><mi>u</mi><mi>j</mi></msub><mo>-</mo><msub><mi>x</mi><mi>i</mi></msub><mo>|</mo><mo>+</mo><mo>|</mo><msub><mi>v</mi><mi>j</mi></msub><mo>-</mo><msub><mi>y</mi><mi>i</mi></msub><mo>|</mo><mo>&rsqb;</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000960278230000014.GIF" wi="1286" he="150" /></maths><maths num="0005"><math><![CDATA[<mrow><mi>s</mi><mo>.</mo><mi>t</mi><mo>.</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><msub><mi>&delta;</mi><mrow><mi>j</mi><mi>i</mi></mrow></msub><mo>=</mo><mn>1</mn><mo>,</mo><mrow><mo>(</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn><mo>,</mo><mn>2</mn><mo>,</mo><mn>3</mn><mo>,</mo><mo>...</mo><mo>,</mo><mi>n</mi></mrow><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000960278230000015.GIF" wi="1198" he="149" /></maths>式中:P<sub>i</sub>—各负荷点功率kWA<sub>i</sub>—给负荷点i供电的配电线路的截面积β<sub>i</sub>—单位距离、单位负荷的费用系数<img file="FDA0000960278230000021.GIF" wi="694" he="151" />n—负荷点的总个数m—变电站的总个数(u<sub>j,</sub>v<sub>j</sub>)—第j个变电站的站址坐标(x<sub>i</sub>,y<sub>i</sub>)—第i个负荷点的坐标将四个原则综合考虑,引入经济性适应度评价函数E(f),即<maths num="0006"><math><![CDATA[<mrow><mi>E</mi><mrow><mo>(</mo><mi>f</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>&omega;</mi><mn>1</mn></msub><mfrac><msub><mi>f</mi><mn>1</mn></msub><mrow><mi>min</mi><mi> </mi><msub><mi>f</mi><mn>1</mn></msub></mrow></mfrac><mo>+</mo><msub><mi>&omega;</mi><mn>2</mn></msub><mfrac><msub><mi>f</mi><mn>2</mn></msub><mrow><mi>min</mi><mi> </mi><msub><mi>f</mi><mn>2</mn></msub></mrow></mfrac><mo>+</mo><msub><mi>&omega;</mi><mn>3</mn></msub><mfrac><msub><mi>f</mi><mn>3</mn></msub><mrow><mi>min</mi><mi> </mi><msub><mi>f</mi><mn>3</mn></msub></mrow></mfrac><mo>+</mo><msub><mi>&omega;</mi><mn>4</mn></msub><mfrac><msub><mi>f</mi><mn>4</mn></msub><mrow><mi>min</mi><mi> </mi><msub><mi>f</mi><mn>4</mn></msub></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000960278230000022.GIF" wi="1261" he="151" /></maths>式中有E(f)=E(U,V),ω<sub>i</sub>(i=1,2,3,4)表示评价函数中各目标函数的相应的权系数;为求解以上权系数,引入均差排序法,求解上述各目标函数赋予的权重,步骤如下:先在其可行域上极小化各目标函数f<sub>j</sub>(u,v),得极小点(u<sup>j</sup>,v<sup>j</sup>),<maths num="0007"><math><![CDATA[<mrow><msub><mi>f</mi><mi>i</mi></msub><mrow><mo>(</mo><msup><mi>u</mi><mi>i</mi></msup><mo>,</mo><msup><mi>v</mi><mi>i</mi></msup><mo>)</mo></mrow><mo>=</mo><munder><mrow><mi>m</mi><mi>i</mi><mi>n</mi></mrow><mrow><mo>(</mo><mi>u</mi><mo>,</mo><mi>v</mi><mo>)</mo><mo>&Element;</mo><mo>(</mo><mi>U</mi><mo>,</mo><mi>V</mi><mo>)</mo></mrow></munder><msub><mi>f</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>u</mi><mo>,</mo><mi>v</mi><mo>)</mo></mrow><mo>,</mo><mi>i</mi><mo>=</mo><mn>1</mn><mo>,</mo><mn>2</mn><mo>,</mo><mo>...</mo><mo>,</mo><mi>m</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>7</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000960278230000023.GIF" wi="1070" he="103" /></maths>利用这m个极小点,可以算出第i(i=1,2,…,m)个目标关于各个极小点的离差<maths num="0008"><math><![CDATA[<mrow><msubsup><mi>&delta;</mi><mi>i</mi><mi>j</mi></msubsup><mo>=</mo><msub><mi>f</mi><mi>i</mi></msub><mrow><mo>(</mo><msup><mi>u</mi><mi>j</mi></msup><mo>,</mo><msup><mi>v</mi><mi>j</mi></msup><mo>)</mo></mrow><mo>-</mo><msub><mi>f</mi><mi>i</mi></msub><mrow><mo>(</mo><msup><mi>u</mi><mi>i</mi></msup><mo>,</mo><msup><mi>v</mi><mi>i</mi></msup><mo>)</mo></mrow><mo>,</mo><mi>j</mi><mo>=</mo><mn>1</mn><mo>,</mo><mn>2</mn><mo>,</mo><mo>...</mo><mo>,</mo><mi>m</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>8</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000960278230000024.GIF" wi="1062" he="87" /></maths>当j=i时,有<img file="FDA0000960278230000025.GIF" wi="607" he="86" />除了上式外,第i个目标关于其他各目标的离差共为m‑1个,对这些离差做算术平均可得到第i个目标关于各极小值点得的平均离差:<img file="FDA0000960278230000026.GIF" wi="540" he="207" />并且Δ<sub>i</sub>≥0,为使权系数规范化,可取<img file="FDA0000960278230000031.GIF" wi="520" he="194" />作为目标f<sub>i</sub>对应的权系数,所以评价函数可表示为<maths num="0009"><math><![CDATA[<mrow><mi>E</mi><mrow><mo>(</mo><mi>u</mi><mo>,</mo><mi>v</mi><mo>)</mo></mrow><mo>=</mo><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mn>4</mn></munderover><msub><mi>&omega;</mi><mi>i</mi></msub><msub><mover><mi>f</mi><mo>&OverBar;</mo></mover><mi>i</mi></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>9</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000960278230000032.GIF" wi="1036" he="151" /></maths>式中<img file="FDA0000960278230000033.GIF" wi="47" he="86" />为f<sub>i</sub>的标量化,通过评价函数评分,选出k个候选变电站站址;7)运用多层次分析法分析候选变电站站址的地理适应度,得出各候选变电站站址的地理适应度;使用AHP法进行地理适应度指标分析,这里引入地理适应度指标F(U,V),即<maths num="0010"><math><![CDATA[<mrow><mi>F</mi><mrow><mo>(</mo><mi>U</mi><mo>,</mo><mi>V</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mn>1</mn><mi>m</mi></mfrac><munder><mo>&Sigma;</mo><mrow><mi>j</mi><mo>&Element;</mo><mi>J</mi><mrow><mo>(</mo><mi>U</mi><mo>,</mo><mi>V</mi><mo>)</mo></mrow></mrow></munder><mfrac><msub><mi>c</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub><msub><mi>c</mi><mrow><mi>i</mi><mi>m</mi><mi>a</mi><mi>x</mi></mrow></msub></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>10</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000960278230000034.GIF" wi="990" he="151" /></maths>式中:m—变电站的总个数c<sub>ij</sub>—第i个新建变电站的第j个候选站址评分c<sub>imax</sub>—第i个新建站所有候选站址中评分最高值对用地性质、交通情况、施工条件、防洪排水、对通信的干扰和地形地质这6个方面进行评分,确定初选地址的地理适应度指标值;8)根据经济性指标和地理适应度组合所形成的评价函数,求解各候选变电站站址的综合适应度;引入综合适应度H(U,V),作为经济适应度指标和地理适应度指标的综合评价函数,即<maths num="0011"><math><![CDATA[<mrow><mi>H</mi><mrow><mo>(</mo><mi>U</mi><mo>,</mo><mi>V</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>&beta;</mi><mn>1</mn></msub><mfrac><mn>1</mn><mrow><mi>E</mi><mrow><mo>(</mo><mi>U</mi><mo>,</mo><mi>V</mi><mo>)</mo></mrow></mrow></mfrac><mo>+</mo><msub><mi>&beta;</mi><mn>2</mn></msub><mi>F</mi><mrow><mo>(</mo><mi>U</mi><mo>,</mo><mi>V</mi><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>11</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000960278230000035.GIF" wi="1131" he="149" /></maths>式中,β<sub>1</sub>,β2表示评价函数中各目标函数的相应的权系数;9)将各候选变电站站址,按照综合适应度的高低进行排序,并将最优和次优的站址作为本次规划的优选站址;10)在步骤9)中选出的优选站址和电网现状评估的基础上,对电网进行网架优化;根据电网规划的经济性和可靠性分析,以供应方开发成本的贴现值最小和需求方缺电成本的贴现值最小为多目标电网规划问题的优化目标分别为:<maths num="0012"><math><![CDATA[<mrow><mi>min</mi><mi> </mi><msub><mi>f</mi><mn>1</mn></msub><mo>=</mo><munderover><mo>&Sigma;</mo><mrow><mi>t</mi><mo>=</mo><mn>1</mn></mrow><mi>T</mi></munderover><mfrac><mrow><msub><mi>C</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>u</mi><mo>(</mo><mrow><mi>t</mi><mo>-</mo><mn>1</mn></mrow><mo>)</mo><mo>)</mo></mrow><mo>+</mo><msub><mi>C</mi><mn>2</mn></msub><mrow><mo>(</mo><mi>s</mi><mo>(</mo><mi>t</mi><mo>)</mo><mo>)</mo></mrow></mrow><msup><mrow><mo>(</mo><mn>1</mn><mo>+</mo><mi>r</mi><mo>)</mo></mrow><mrow><msub><mi>y</mi><mi>a</mi></msub><mrow><mo>(</mo><mrow><mi>t</mi><mo>-</mo><mn>1</mn></mrow><mo>)</mo></mrow></mrow></msup></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>12</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000960278230000041.GIF" wi="1051" he="151" /></maths><maths num="0013"><math><![CDATA[<mrow><mi>min</mi><mi> </mi><msub><mi>f</mi><mn>2</mn></msub><mo>=</mo><munderover><mo>&Sigma;</mo><mrow><mi>t</mi><mo>=</mo><mn>1</mn></mrow><mi>T</mi></munderover><mfrac><mrow><msub><mi>C</mi><mn>3</mn></msub><mrow><mo>(</mo><mi>s</mi><mo>(</mo><mi>t</mi><mo>)</mo><mo>)</mo></mrow></mrow><msup><mrow><mo>(</mo><mn>1</mn><mo>+</mo><mi>r</mi><mo>)</mo></mrow><mrow><msub><mi>y</mi><mi>a</mi></msub><mrow><mo>(</mo><mrow><mi>t</mi><mo>-</mo><mn>1</mn></mrow><mo>)</mo></mrow></mrow></msup></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>13</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000960278230000042.GIF" wi="830" he="151" /></maths>s(t)∈S(t)u(t)∈U(t)<maths num="0014"><math><![CDATA[<mrow><msub><mi>P</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>&le;</mo><msub><mover><mi>P</mi><mo>&OverBar;</mo></mover><mrow><mi>i</mi><mi>j</mi></mrow></msub></mrow>]]></math><img file="FDA0000960278230000043.GIF" wi="228" he="87" /></maths><maths num="0015"><math><![CDATA[<mrow><msubsup><mi>P</mi><mrow><mi>i</mi><mi>j</mi></mrow><mo>&prime;</mo></msubsup><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>&le;</mo><msub><mover><mi>P</mi><mo>&OverBar;</mo></mover><mrow><mi>i</mi><mi>j</mi></mrow></msub></mrow>]]></math><img file="FDA0000960278230000044.GIF" wi="228" he="87" /></maths>式中:T—规划期内的总阶段数s(t)—第t阶段的网络状态u(t)—从第t阶段到第t+1阶段的网络扩展方案C<sub>1</sub>(u(t‑1))—第t阶段新增线路的投资费用C<sub>2</sub>(s(t))—按方案u(t‑1)扩展网络到状态s(t)后网络的运行费用(包括网损费用和维护费用)C<sub>3</sub>(s(t))—第t阶段的缺电成本y<sub>a</sub>(t)—规划初期到第t阶段末的总年数引入评价函数Q(U,V),令<maths num="0016"><math><![CDATA[<mrow><mi>Q</mi><mrow><mo>(</mo><mi>f</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>&eta;</mi><mn>1</mn></msub><mfrac><msub><mi>f</mi><mn>1</mn></msub><mrow><mi>min</mi><mi> </mi><msub><mi>f</mi><mn>1</mn></msub></mrow></mfrac><mo>+</mo><msub><mi>&eta;</mi><mn>2</mn></msub><mfrac><msub><mi>f</mi><mn>2</mn></msub><mrow><mi>min</mi><mi> </mi><msub><mi>f</mi><mn>2</mn></msub></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>14</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000960278230000045.GIF" wi="894" he="150" /></maths>式中η<sub>1</sub>,η<sub>2</sub>分别表示上式目标函数的权系数;网架优化规划分两阶段进行,第一阶段利用遗传算法求解单目标函数,第二阶段利用均差排序法求解出评价函数,其中评价函数值越小越好,表示供应方和需求方的损失达到最小,通过评价函数值大小排序来确定最优和次优方案;11)根据步骤10)中得到的网架优化的结果,选取最优和次优的两套方案,作为优选网架方案;12)根据优选变电站站址和优选网架方案的结果,制定规划方案,并进行技术经济评估,得出一个最佳方案。
地址 210098 江苏省南京市鼓楼区西康路1号