发明名称 受端电网安全域最优潮流模型的建模方法
摘要 本发明公开了一种受端电网安全域最优潮流模型的建模方法,具体步骤为:第一步:通过对潮流方程的稳定性分析即对系统静态电压、相角、振荡频率的稳定性和N-1运行准则来构建系统安全域,系统的N个负荷可由给定发电调度准则变换为负荷方向组成的M个不同集合,生成临界负荷矩阵来对安全域进行近似处理:第二步:构建动态安全约束最优潮流模型;第三步:建立自适应神经模糊推理系统;第四步:训练自适应模糊推理系统及建立最优潮流模型。本发明所述的最优潮流模型很好地阐释了当今电力系统调度的运行准则,基于自适应模糊推理系统的最优潮流安全域逼近技术可运用于分层分区后受端电网区域间的调度优化。
申请公布号 CN102841965B 申请公布日期 2015.10.28
申请号 CN201210301629.0 申请日期 2012.08.23
申请人 山东电力集团公司电力科学研究院;国家电网公司 发明人 麻常辉;张磊;蒋哲;武乃虎
分类号 G06F17/50(2006.01)I 主分类号 G06F17/50(2006.01)I
代理机构 济南圣达知识产权代理有限公司 37221 代理人 张勇
主权项 一种受端电网安全域最优潮流模型的建模方法,其特征是,具体步骤为:第一步:通过对潮流方程的稳定性分析即对系统静态电压、相角、振荡频率的稳定性和N‑1准则来构建系统安全域,为获得离散的安全域表达式,系统的N个负荷由给定发电调度准则变换为负荷方向组成的M个不同集合,N、M为大于等于1的自然数,N为系统中负荷数量,给定的系统具有固定的负荷数量,生成临界负荷矩阵来对安全域进行近似处理;第二步:构建动态安全约束最优潮流模型;第三步:建立自适应神经模糊推理系统;第四步:训练自适应模糊推理系统及建立最优潮流模型;所述第一步的具体步骤为:(1)建立电力系统可微代数方程为:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mfenced open='' close=''><mtable><mtr><mtd><mover><mi>x</mi><mo>&CenterDot;</mo></mover><mo>=</mo><mi>h</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>,</mo><mi>&rho;</mi><mo>,</mo><mi>&lambda;</mi><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mn>0</mn><mo>=</mo><mi>g</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>,</mo><mi>&rho;</mi><mo>,</mo><mi>&lambda;</mi><mo>)</mo></mrow></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>22</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000698761190000011.GIF" wi="461" he="162" /></maths>其中,x是系统状态变量,y代表代数变量,ρ表示系统可控变量;λ是一组不可控参量;(2)确定负荷方向d<sub>i</sub>=[d<sub>i1</sub>d<sub>i2</sub>...d<sub>iN</sub>]<sup>T</sup>,当负荷沿某一特定方向增加时,电力系统将达到运行极限,通过对潮流方程的稳定性分析来构建系统安全域,即根据对系统静态电压、相角、振荡频率进行稳定性和N‑1运行准则来构建系统安全域,<img file="FDA0000698761190000014.GIF" wi="48" he="63" />是功角,有功功率、无功功率和视在功率组成直角三角形,视在功率为斜边,有功与视在功率之间的夹角叫功角,一般用其余弦值表示功率因数,G是一个虚数,m是1~G中的一个数;设λ<sub>i</sub>=[λ<sub>i1</sub>λ<sub>i2</sub>...λ<sub>iN</sub>]<sup>T</sup>为N个负荷中第i个负荷的增长率,i为大于等于1的自然数,λ表示为:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mfenced open='' close=''><mtable><mtr><mtd><msub><mi>&lambda;</mi><mrow><mi>i</mi><mn>1</mn></mrow></msub><mo>=</mo><mi>&alpha;</mi><msub><mi>d</mi><mrow><mi>i</mi><mn>1</mn></mrow></msub></mtd></mtr><mtr><mtd><msub><mi>&lambda;</mi><mrow><mi>i</mi><mn>2</mn></mrow></msub><mo>=</mo><mi>&alpha;</mi><msub><mi>d</mi><mrow><mi>i</mi><mn>2</mn></mrow></msub></mtd></mtr><mtr><mtd><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo></mtd></mtr><mtr><mtd><msub><mi>&lambda;</mi><mi>iN</mi></msub><mo>=</mo><mi>&alpha;</mi><msub><mi>d</mi><mi>iN</mi></msub></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>23</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000698761190000012.GIF" wi="398" he="301" /></maths>其中,标量α≥0代表负荷系数,d<sub>ij</sub>表示负荷j在第i个负荷增长率下的负荷增长方向,i和j均为大于等于1的自然数,负荷方向满足如下条件:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mfenced open='' close=''><mtable><mtr><mtd><mn>0</mn><mo>&le;</mo><msub><mi>d</mi><mi>ij</mi></msub><mo>&le;</mo><mn>1</mn></mtd><mtd><mo>&ForAll;</mo><mi>j</mi></mtd></mtr><mtr><mtd><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><msub><mi>d</mi><mi>ij</mi></msub><mo>=</mo><mn>1</mn></mtd><mtd></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>24</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000698761190000013.GIF" wi="492" he="231" /></maths>(3)通过增加负荷系数α能使系统逐渐达到安全稳定边界,进而确定稳定边界极限值<img file="FDA0000698761190000027.GIF" wi="79" he="91" />系统的N个负荷可由给定发电调度准则变换为负荷方向组成的M个不同集合,生成如下临界负荷矩阵来对安全域进行近似:<img file="FDA0000698761190000021.GIF" wi="1237" he="485" />所述第二步的最优潮流模型如下:目标函数:<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msub><mi>S</mi><mi>b</mi></msub><mo>=</mo><mo>-</mo><mrow><mo>(</mo><msubsup><mi>C</mi><mi>d</mi><mi>T</mi></msubsup><msub><mi>P</mi><mi>d</mi></msub><mo>-</mo><msubsup><mi>C</mi><mi>s</mi><mi>T</mi></msubsup><msub><mi>P</mi><mi>s</mi></msub><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000698761190000026.GIF" wi="595" he="91" /></maths>约束条件:F<sub>PF</sub>(δ,V,Q<sub>g</sub>,P<sub>s</sub>,P<sub>d</sub>)=0   (3)<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><mn>0</mn><mo>&le;</mo><msub><mi>P</mi><mi>s</mi></msub><mo>&le;</mo><msub><mi>P</mi><msub><mi>s</mi><mi>max</mi></msub></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000698761190000025.GIF" wi="369" he="90" /></maths><maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><mn>0</mn><mo>&le;</mo><msub><mi>P</mi><mi>d</mi></msub><mo>&le;</mo><msub><mi>P</mi><msub><mi>d</mi><mi>max</mi></msub></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000698761190000024.GIF" wi="457" he="90" /></maths><maths num="0007" id="cmaths0007"><math><![CDATA[<mrow><msub><mi>I</mi><mi>ij</mi></msub><mrow><mo>(</mo><mi>&delta;</mi><mo>,</mo><mi>V</mi><mo>)</mo></mrow><mo>&le;</mo><msub><mi>I</mi><msub><mi>ij</mi><mi>max</mi></msub></msub><mo>,</mo><mo>&ForAll;</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>,</mo><mi>i</mi><mo>&NotEqual;</mo><mi>j</mi></mrow>]]></math><img file="FDA0000698761190000023.GIF" wi="615" he="99" /></maths>i、j为大于等于1的自然数(6)Q<sub>gmin</sub>≤Q<sub>g</sub>≤Q<sub>gmax</sub>   (7)V<sub>min</sub>≤V≤V<sub>max</sub>   (8)f<sub>NR0</sub>‑f<sub>NR</sub>(V,δ)≤0   (9)其中,C<sub>s</sub>和C<sub>d</sub>分别是电力供给与需求的出价,单位为$/MWh;系统供给与需求功率分别为P<sub>s</sub>和P<sub>d</sub>,单位是MW;F<sub>PF</sub>(·)为系统潮流方程;V和δ分别是节点电压和相交;I<sub>ij</sub>是通过输电线路ij的电流,该约束定义了系统的热稳定极限;Q<sub>g</sub>为发电机无功功率;f<sub>NR</sub>(·)用于表示系统安全域,<img file="FDA0000698761190000022.GIF" wi="81" he="75" />为其合适的临界值;所述第三步中,自适应神经模糊推理系统分为六层:X<sub>1</sub>,X<sub>2</sub>是系统的输入,y是推理系统输出;网络同一层的每个节点具有相似的功能,用O<sub>1i</sub>表示第i个节点输出,i为大于等于1的自然数;第一层:将输入数据进行模糊处理:O<sub>1i</sub>=μ<sub>Ai</sub>(x<sub>1</sub>),O<sub>2i</sub>=μ<sub>Bi</sub>(x<sub>2</sub>),i=1,2   (25)其中,A<sub>i</sub>或B<sub>i</sub>是模糊集;μ<sub>Ai</sub>(x<sub>1</sub>)是模糊集的隶属度函数;第二层:将各输入数据隶属函数相乘,作为本层规则的适用度w<sub>i</sub>:w<sub>i</sub>=μ<sub>Ai</sub>(x<sub>1</sub>)μ<sub>Bi</sub>(x<sub>2</sub>),i=1,2   (26)第三层:计算第i条规则的w<sub>i</sub>及所有适用度之和w<sub>1</sub>+w<sub>2</sub>,并通过两者比值完成各条规则适用度的归一化:<maths num="0008" id="cmaths0008"><math><![CDATA[<mrow><msubsup><mi>w</mi><mi>i</mi><mo>&prime;</mo></msubsup><mo>=</mo><mfrac><msub><mi>w</mi><mi>i</mi></msub><mrow><msub><mi>w</mi><mn>1</mn></msub><mo>+</mo><msub><mi>w</mi><mn>2</mn></msub></mrow></mfrac><mo>,</mo><mi>i</mi><mo>=</mo><mn>1,2</mn><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>27</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000698761190000031.GIF" wi="661" he="132" /></maths>第四层:用于计算各条规则的输出:<maths num="0009" id="cmaths0009"><math><![CDATA[<mrow><msubsup><mi>O</mi><mrow><mn>4</mn><mi>i</mi></mrow><mo>&prime;</mo></msubsup><mo>=</mo><msubsup><mi>w</mi><mi>i</mi><mo>&prime;</mo></msubsup><msub><mi>f</mi><mi>i</mi></msub><mo>=</mo><msubsup><mi>w</mi><mi>i</mi><mo>&prime;</mo></msubsup><mrow><mo>(</mo><msub><mi>p</mi><mi>i</mi></msub><msub><mi>x</mi><mn>1</mn></msub><mo>+</mo><msub><mi>q</mi><mi>i</mi></msub><msub><mi>x</mi><mn>2</mn></msub><mo>+</mo><msub><mi>r</mi><mi>i</mi></msub><mo>)</mo></mrow><mo>,</mo><mi>i</mi><mo>=</mo><mn>1,2</mn><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>28</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000698761190000036.GIF" wi="1008" he="103" /></maths>其中,f<sub>i</sub>为模糊系统的后项结论输出函数,p<sub>i</sub>、q<sub>i</sub>为第i条规则下的系统加权系数,r<sub>i</sub>为第i条规下的常量,当该输出函数为线性函数时,称为“一阶系统”;若为常量,称为“零阶系统”;第五层:用于计算系统的总输出:<maths num="0010" id="cmaths0010"><math><![CDATA[<mrow><mi>y</mi><mo>=</mo><munder><mi>&Sigma;</mi><mi>i</mi></munder><msubsup><mi>w</mi><mi>i</mi><mo>&prime;</mo></msubsup><msub><mi>f</mi><mi>i</mi></msub><mo>=</mo><mfrac><mrow><munder><mi>&Sigma;</mi><mi>i</mi></munder><msub><mi>w</mi><mi>i</mi></msub><msub><mi>f</mi><mi>i</mi></msub></mrow><mrow><munder><mi>&Sigma;</mi><mi>i</mi></munder><msub><mi>w</mi><mi>i</mi></msub></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>29</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000698761190000032.GIF" wi="662" he="215" /></maths>第六层,通过加权平均法将该输出结果进行解模糊化处理,并通过反向传播法和最小二乘法使输入与输出之间的误差最小;所述第四步中的具体步骤为:(1)训练自适应模糊推理系统;使N‑1个负荷的M个极限值<img file="FDA0000698761190000035.GIF" wi="66" he="64" />所组成的边界作为自适应模糊推理系统的输入,N为系统中负荷数量,给定的系统具有固定的负荷数量,M为自然数,定义第i个负荷节点的安全边界值<img file="FDA0000698761190000034.GIF" wi="54" he="68" />表示为:<maths num="0011" id="cmaths0011"><math><![CDATA[<mrow><msubsup><mi>&lambda;</mi><mi>il</mi><mi>c</mi></msubsup><mo>&ap;</mo><mi>f</mi><mrow><mo>(</mo><msub><mi>&lambda;</mi><mrow><mi>i</mi><mn>1</mn></mrow></msub><mo>,</mo><msub><mi>&lambda;</mi><mrow><mi>i</mi><mn>2</mn></mrow></msub><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><msub><mi>&lambda;</mi><mrow><mi>il</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>,</mo><msub><mi>&lambda;</mi><mrow><mi>il</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><msub><mi>&lambda;</mi><mi>iN</mi></msub><mo>)</mo></mrow><mo>&ap;</mo><mi>f</mi><mrow><mo>(</mo><msub><mover><mi>&lambda;</mi><mo>^</mo></mover><mi>i</mi></msub><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000698761190000033.GIF" wi="1188" he="136" /></maths>i为大于等于1的自然数:(30)由公式(9)和公式(25)得负荷增长率的映射函数为:<maths num="0012" id="cmaths0012"><math><![CDATA[<mrow><msubsup><mi>&lambda;</mi><mi>l</mi><mi>c</mi></msubsup><mo>=</mo><mfrac><mrow><munder><mi>&Sigma;</mi><mi>i</mi></munder><msub><mi>w</mi><mi>i</mi></msub><msub><mi>f</mi><mi>i</mi></msub></mrow><mrow><munder><mi>&Sigma;</mi><mi>i</mi></munder><msub><mi>w</mi><mi>i</mi></msub></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>31</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000698761190000041.GIF" wi="571" he="282" /></maths>(2)将公式(26)用于最优潮流方程的安全域约束中,形成安全域约束最优潮流模型为:目标函数:<maths num="0013" id="cmaths0013"><math><![CDATA[<mrow><msub><mi>S</mi><mi>b</mi></msub><mo>=</mo><msubsup><mi>C</mi><mi>s</mi><mi>T</mi></msubsup><msub><mi>P</mi><mi>s</mi></msub><mo>-</mo><msubsup><mi>C</mi><mi>d</mi><mi>T</mi></msubsup><mi>&Delta;</mi><msub><mi>P</mi><mi>d</mi></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>10</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000698761190000048.GIF" wi="543" he="84" /></maths>约束条件:F<sub>PF</sub>(δ,V,Q<sub>g</sub>,P<sub>s</sub>,P<sub>d</sub>,Q<sub>d</sub>)=0   (11)0≤P<sub>s</sub>≤P<sub>smax</sub>   (12)Q<sub>smin</sub>≤Q<sub>s</sub>≤Q<sub>smax</sub>   (13)V<sub>min</sub>≤V≤V<sub>max</sub>   (14)<maths num="0014" id="cmaths0014"><math><![CDATA[<mrow><msub><mi>&lambda;</mi><mi>l</mi></msub><mo>-</mo><mfrac><mrow><munder><mrow><mi>&Sigma;</mi><msub><mi>w</mi><msub><mi>i</mi><mi>m</mi></msub></msub></mrow><mi>i</mi></munder><msub><mi>f</mi><msub><mi>i</mi><mi>m</mi></msub></msub></mrow><mrow><munder><mi>&Sigma;</mi><mi>i</mi></munder><msub><mi>w</mi><msub><mi>i</mi><mi>m</mi></msub></msub></mrow></mfrac><mo>&le;</mo><mn>0</mn><mo>,</mo><mo>&ForAll;</mo><mi>m</mi><mo>=</mo><mn>1</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>G</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>15</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000698761190000042.GIF" wi="845" he="216" /></maths><maths num="0015" id="cmaths0015"><math><![CDATA[<mrow><msub><mi>&Delta;P</mi><mi>dj</mi></msub><mo>&le;</mo><mn>0</mn><mo>,</mo><mo>&ForAll;</mo><mi>j</mi><mo>=</mo><mn>1</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>N</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>16</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000698761190000047.GIF" wi="618" he="89" /></maths><maths num="0016" id="cmaths0016"><math><![CDATA[<mrow><msub><mi>&Delta;P</mi><msub><mi>d</mi><mi>j</mi></msub></msub><mo>=</mo><mrow><mo>(</mo><mi>&alpha;</mi><msub><mi>d</mi><mi>j</mi></msub><mo>-</mo><msub><mi>&alpha;</mi><mn>0</mn></msub><msub><mi>d</mi><mrow><mi>j</mi><mn>0</mn></mrow></msub><mo>)</mo></mrow><msub><mi>P</mi><mrow><mi>dj</mi><mn>0</mn></mrow></msub><mo>,</mo><mo>&ForAll;</mo><mi>j</mi><mo>=</mo><mn>1</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>N</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>17</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000698761190000046.GIF" wi="958" he="88" /></maths><img file="FDA0000698761190000043.GIF" wi="797" he="94" /><maths num="0017" id="cmaths0017"><math><![CDATA[<mrow><mn>0</mn><mo>&le;</mo><msub><mi>d</mi><mi>j</mi></msub><mo>&le;</mo><mn>1</mn><mo>,</mo><mo>&ForAll;</mo><mi>j</mi><mo>=</mo><mn>1</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>N</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>19</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000698761190000045.GIF" wi="618" he="84" /></maths><maths num="0018" id="cmaths0018"><math><![CDATA[<mrow><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><msub><mi>d</mi><mi>j</mi></msub><mo>=</mo><mn>1</mn><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>20</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000698761190000044.GIF" wi="358" he="147" /></maths>α≥0   (21);其中,P<sub>dj0</sub>为第j个节点的初始负荷,P<sub>dj</sub>为第j个节点的负荷,ΔP<sub>d</sub>为负荷改变量;m为所有G个调度方案的第m个系统安全域,约束条件(11)强行使ΔP<sub>d</sub>为0或负数;若ΔP<sub>d</sub>为0,则最优潮流模型有解;反之,若ΔP<sub>d</sub>为负数,则表示最优潮流模型无解,因此,该最优潮流模型很好地阐释了当今电力系统调度的运行准则,<img file="FDA0000698761190000051.GIF" wi="48" he="63" />是功角,有功功率、无功功率和视在功率组成直角三角形,视在功率为斜边,有功与视在功率之间的夹角叫功角,一般用其余弦值表示功率因数,G是一个虚数,m是1~G中的一个数。
地址 250002 山东省济南市市中区二环南路500号