发明名称 一种基于需求响应的负荷建模和优化控制方法
摘要 本发明公开了一种基于需求响应的负荷建模和优化控制方法,将负荷进行刚性负荷和柔性负荷的解耦分类,建立负荷模型,进一步地,将柔性负荷根据参与的需求响应项目分为三类:电价项目、直接负荷控制项目、可中断负荷项目,分类建立电价、直接负荷和可中断负荷模型,模型的建立为后续的负荷优化提供有力的数据资源;根据可控变量所能提供的最大消减负荷,直接负荷控制项目最大削减负荷、可中断负荷项目最大削减负荷与供电缺口最大值之间的关系建立三种情况的优化控制模型,在通过遗传算法进行求解,求得优化执行结果,进而为实际电力系统运行提供更多资源和策略。
申请公布号 CN103972896B 申请公布日期 2015.09.09
申请号 CN201410200869.0 申请日期 2014.05.13
申请人 国家电网公司;国网河南省电力公司电力科学研究院;国网河南省电力公司;中国电力科学研究院 发明人 张景超;安庆;李强;孙芊;王倩;牛雨;李帅;屈博;陈宋宋
分类号 H02J3/06(2006.01)I 主分类号 H02J3/06(2006.01)I
代理机构 郑州联科专利事务所(普通合伙) 41104 代理人 刘建芳
主权项 一种基于需求响应的负荷建模和优化控制方法,其特征在于:包括以下步骤:a:按照负荷的可控性对用电负荷进行刚性负荷和柔性负荷的解耦分类,柔性负荷是指通过技术手段可以转移或削减的负荷,且该过程具有成本效益,时间跨度符合要求;刚性负荷是指通过技术手段不可以转移或削减的负荷;b:建立负荷模型;负荷模型由刚性负荷和柔性负荷两部分组成,用于表示负荷随时间变化的曲线,以及负荷与需求响应指令之间的函数关系;负荷模型采用公式①,为:L<sub>total</sub>(t,P,DLC,IL)=L<sub>s</sub>(t)+L<sub>f</sub>(t,P,DLC,IL)   ①式中,L<sub>total</sub>表示总负荷;L<sub>s</sub>表示刚性负荷,是时间的函数;L<sub>f</sub>表示柔性负荷,是时间和需求响应指令的函数;t表示时间;需求响应指令:包括电价P,直接负荷控制变量DLC,可中断负荷变量IL;需求响应指令本身也是时间的函数,即不同时刻点有不同的需求响应指令;c:将柔性负荷根据参与的需求响应项目分为三类:电价项目、直接负荷控制项目、可中断负荷项目,并按照三项需求响应项目,根据公式②,L<sub>f</sub>(t,P,DLC,IL)=L<sub>1</sub>(t,P)+L<sub>2</sub>(t,DLC)+L<sub>3</sub>(t,IL)②分类建立模型;(1):与电价项目相关的柔性负荷模型:用户的电力需求随电价的变化而变化,但各类用户的变化趋势和幅度各不相同,整体变化规律可用抛物线模型或对数模型进行表示,表示采用公式③,公式③为:L<sub>1</sub>(t,P)=a+bP(t)+cP<sup>2</sup>(t)或L<sub>1</sub>(t,P)=a+bln[P(t)]③;(2):与直接负荷控制项目相关的柔性负荷模型:直接负荷控制项目由实施机构通过远程关闭或循环控制用户的用电设备,控制过程中需要满足用户满意度的约束条件,采用公式④,公式④为:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>L</mi><mn>2</mn></msub><mrow><mo>(</mo><mi>t</mi><mo>,</mo><mi>D</mi><mi>L</mi><mi>C</mi><mo>)</mo></mrow><mo>=</mo><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><mo>{</mo><mo>&lsqb;</mo><mn>1</mn><mo>-</mo><msub><mi>&alpha;</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>&rsqb;</mo><mo>&CenterDot;</mo><msub><mi>L</mi><mrow><mi>o</mi><mi>r</mi><mi>i</mi><mi>g</mi><mi>i</mi><mi>n</mi><mo>,</mo><mi>i</mi></mrow></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>}</mo></mrow>]]></math><img file="FDA0000743035740000021.GIF" wi="869" he="155" /></maths>④;约束条件为:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mi>s</mi><mo>.</mo><mi>t</mi><mo>.</mo><munderover><mo>&Sigma;</mo><mrow><mi>t</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msub><mi>&alpha;</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>&le;</mo><msub><mi>&alpha;</mi><mrow><mi>i</mi><mo>_</mo><mi>m</mi><mi>a</mi><mi>x</mi></mrow></msub><mo>;</mo></mrow>]]></math><img file="FDA0000743035740000022.GIF" wi="588" he="158" /></maths>式中,i表示参与直接负荷控制项目的用户编号,一共有N组用户;α<sub>i</sub>(t)表示第i组用户在t时段的控制变量,为0‑1决策变量;L<sub>origin,i</sub>表示第i组用户实施直接负荷控制之前的负荷大小;α<sub>i_max</sub>为第i组用户合同中规定的每天调用次数上限,n为一天的时段划分总数;(3):与可中断负荷控制项目相关的柔性负荷模型:可中断负荷是基于合同的需求响应项目,负荷的中断由用户根据请求信号实行,所以存在用户违约率的问题,可中断负荷的实行以合同的形式,规定了最小切除容量、中断持续时间等约束条件,采用公式⑤,公式⑤为:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mi>L</mi><mn>3</mn></msub><mrow><mo>(</mo><mi>t</mi><mo>,</mo><mi>I</mi><mi>L</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>L</mi><mrow><mi>o</mi><mi>r</mi><mi>i</mi><mi>g</mi><mi>i</mi><mi>n</mi></mrow></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>-</mo><munderover><mo>&Sigma;</mo><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><msub><mi>&beta;</mi><mi>j</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>&CenterDot;</mo><msub><mi>L</mi><mrow><mi>I</mi><mi>L</mi><mo>,</mo><mi>j</mi></mrow></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>&CenterDot;</mo><mo>&lsqb;</mo><mn>1</mn><mo>-</mo><msub><mi>&delta;</mi><mi>j</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>&rsqb;</mo></mrow>]]></math><img file="FDA0000743035740000023.GIF" wi="1096" he="163" /></maths>⑤;约束条件为:s.t. L<sub>IL,j</sub>(t)≥L<sub>IL,j_min</sub><maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msub><mi>D</mi><mrow><mi>j</mi><mo>_</mo><mi>min</mi></mrow></msub><mo>&le;</mo><munderover><mo>&Sigma;</mo><mrow><mi>t</mi><mo>=</mo><mn>1</mn></mrow><mn>24</mn></munderover><msub><mi>&beta;</mi><mi>j</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>&lt;</mo><msub><mi>D</mi><mrow><mi>j</mi><mo>_</mo><mi>m</mi><mi>a</mi><mi>x</mi></mrow></msub><mo>;</mo></mrow>]]></math><img file="FDA0000743035740000024.GIF" wi="649" he="160" /></maths>式中,j表示可中断合同编号,合同总数为M;β<sub>j</sub>(t)表示第j个中断合同在t时段的状态,为0‑1决策变量;L<sub>IL,j</sub>(t)表示第j个中断合同在t时段的中断容量;δ<sub>j</sub>(t)为第j个中断合同在t时段的违约率,δ<sub>j</sub>(t)∈[0%,100%];L<sub>IL,j_max</sub>表示第j个中断合同的最小切除容量;D<sub>j_min</sub>和D<sub>j_max</sub>分别表示第j个中断合同的中断持续时间下限和上限;d:负荷模型中柔性负荷的电价项目为不可控变量,直接负荷控制变量α<sub>i</sub>(t)和可中断负荷变量β<sub>j</sub>(t)为可控变量,根据可控变量所能提供的最大消减负荷与供电缺口最大值之间的关系建立三种情况的优化控制模型,即直接负荷控制项目最大削减负荷DLC<sub>max</sub>、可中断负荷项目最大削减负荷IL<sub>max</sub>与供电缺口最大值P<sub>max</sub>之间的关系建立三种情况的优化控制模型;情况一:在P<sub>max</sub>≥DLC<sub>max</sub>+IL<sub>max</sub>时,建立模型一:<maths num="0005" id="cmaths0005"><math><![CDATA[<mfenced open = '' close = ''><mtable><mtr><mtd><mrow><mi>min</mi><mo>|</mo><mi>P</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>-</mo><msub><mi>IL</mi><mi>max</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>-</mo><mi>D</mi><mi>L</mi><mi>C</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>|</mo></mrow></mtd></mtr><mtr><mtd><mrow><mo>=</mo><mo>|</mo><mi>P</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>-</mo><munderover><mo>&Sigma;</mo><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><msub><mi>&beta;</mi><mi>j</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>&CenterDot;</mo><msub><mi>L</mi><mrow><mi>I</mi><mi>L</mi><mo>,</mo><mi>j</mi></mrow></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>&CenterDot;</mo><mo>&lsqb;</mo><mn>1</mn><mo>-</mo><msub><mi>&delta;</mi><mi>j</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>&rsqb;</mo><mo>-</mo><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><mo>&lsqb;</mo><msub><mi>&alpha;</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>&CenterDot;</mo><msub><mi>L</mi><mrow><mi>o</mi><mi>r</mi><mi>i</mi><mi>g</mi><mi>i</mi><mi>n</mi><mo>,</mo><mi>i</mi></mrow></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>&rsqb;</mo><mo>|</mo></mrow></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000743035740000031.GIF" wi="1425" he="288" /></maths><maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><mi>s</mi><mo>.</mo><mi>t</mi><mo>.</mo><munderover><mo>&Sigma;</mo><mrow><mi>t</mi><mo>=</mo><mn>1</mn></mrow><mn>24</mn></munderover><msub><mi>&beta;</mi><mi>j</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mi>m</mi><mi>i</mi><mi>n</mi><mo>{</mo><msub><mi>D</mi><mrow><mi>j</mi><mo>_</mo><mi>m</mi><mi>a</mi><mi>x</mi></mrow></msub><mo>,</mo><msub><mi>T</mi><mrow><mi>m</mi><mi>a</mi><mi>x</mi></mrow></msub><mo>}</mo></mrow>]]></math><img file="FDA0000743035740000032.GIF" wi="757" he="159" /></maths>L<sub>IL,j</sub>(t)≥L<sub>IL,j_min</sub><maths num="0007" id="cmaths0007"><math><![CDATA[<mrow><munderover><mo>&Sigma;</mo><mrow><mi>t</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msub><mi>&alpha;</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>&le;</mo><msub><mi>&alpha;</mi><mrow><mi>i</mi><mo>_</mo><mi>m</mi><mi>a</mi><mi>x</mi></mrow></msub></mrow>]]></math><img file="FDA0000743035740000033.GIF" wi="387" he="152" /></maths>   ⑥‑‑‑‑式中,T<sub>max</sub>为供电缺口最大持续时间;情况二:在DLC<sub>max</sub>≤P<sub>max</sub><DLC<sub>max</sub>+IL<sub>max</sub>时,建立模型二:<maths num="0008" id="cmaths0008"><math><![CDATA[<mfenced open='' close=''><mtable><mtr><mtd><mrow><mi>min</mi><mo>|</mo><mi>P</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>-</mo><mi>IL</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>-</mo><mi>D</mi><mi>L</mi><mi>C</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>|</mo></mrow></mtd></mtr><mtr><mtd><mrow><mo>=</mo><mo>|</mo><mi>P</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>-</mo><munderover><mo>&Sigma;</mo><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>x</mi></munderover><msub><mi>&beta;</mi><mi>j</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>&CenterDot;</mo><msub><mi>L</mi><mrow><mi>I</mi><mi>L</mi><mo>,</mo><mi>j</mi></mrow></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>&CenterDot;</mo><mo>&lsqb;</mo><mn>1</mn><mo>-</mo><msub><mi>&delta;</mi><mi>j</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>&rsqb;</mo><mo>-</mo><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><mo>&lsqb;</mo><msub><mi>&alpha;</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>&CenterDot;</mo><msub><mi>L</mi><mrow><mi>o</mi><mi>r</mi><mi>i</mi><mi>g</mi><mi>i</mi><mi>n</mi><mo>,</mo><mi>i</mi></mrow></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>&rsqb;</mo><mo>|</mo></mrow></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000743035740000041.GIF" wi="1402" he="281" /></maths><maths num="0009" id="cmaths0009"><math><![CDATA[<mrow><mi>s</mi><mo>.</mo><mi>t</mi><mo>.</mo><munderover><mo>&Sigma;</mo><mrow><mi>t</mi><mo>=</mo><mn>1</mn></mrow><mn>24</mn></munderover><msub><mi>&beta;</mi><mi>j</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mi>m</mi><mi>i</mi><mi>n</mi><mo>{</mo><msub><mi>D</mi><mrow><mi>j</mi><mo>_</mo><mi>m</mi><mi>a</mi><mi>x</mi></mrow></msub><mo>,</mo><msub><mi>T</mi><mrow><mi>m</mi><mi>a</mi><mi>x</mi></mrow></msub><mo>}</mo></mrow>]]></math><img file="FDA0000743035740000042.GIF" wi="746" he="156" /></maths><maths num="0010" id="cmaths0010"><math><![CDATA[<mrow><munderover><mo>&Sigma;</mo><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>x</mi></munderover><msub><mi>L</mi><mrow><mi>I</mi><mi>L</mi><mo>,</mo><mi>j</mi></mrow></msub><mrow><mo>(</mo><msub><mi>t</mi><mrow><mi>m</mi><mi>a</mi><mi>x</mi></mrow></msub><mo>)</mo></mrow><mo>=</mo><msub><mi>P</mi><mrow><mi>m</mi><mi>a</mi><mi>x</mi></mrow></msub><mo>-</mo><msub><mi>DLC</mi><mrow><mi>m</mi><mi>a</mi><mi>x</mi></mrow></msub></mrow>]]></math><img file="FDA0000743035740000043.GIF" wi="658" he="160" /></maths>L<sub>IL,j</sub>(t)≥L<sub>IL,j_min</sub><maths num="0011" id="cmaths0011"><math><![CDATA[<mrow><munderover><mo>&Sigma;</mo><mrow><mi>t</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msub><mi>&alpha;</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>&le;</mo><msub><mi>&alpha;</mi><mrow><mi>i</mi><mo>_</mo><mi>m</mi><mi>a</mi><mi>x</mi></mrow></msub></mrow>]]></math><img file="FDA0000743035740000044.GIF" wi="389" he="150" /></maths>   ⑦式中,t<sub>max</sub>为最大供电缺口发生时刻;情况三:在P<sub>max</sub><DLC<sub>max</sub>时,建立模型三:<maths num="0012" id="cmaths0012"><math><![CDATA[<mrow><mi>m</mi><mi>i</mi><mi>n</mi><mo>|</mo><mi>P</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>-</mo><mi>D</mi><mi>L</mi><mi>C</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>|</mo><mo>=</mo><mo>|</mo><mi>P</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>-</mo><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><mo>&lsqb;</mo><msub><mi>&alpha;</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>&CenterDot;</mo><msub><mi>L</mi><mrow><mi>o</mi><mi>r</mi><mi>i</mi><mi>g</mi><mi>i</mi><mi>n</mi><mo>,</mo><mi>i</mi></mrow></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>&rsqb;</mo><mo>|</mo></mrow>]]></math><img file="FDA0000743035740000045.GIF" wi="1174" he="178" /></maths><maths num="0013" id="cmaths0013"><math><![CDATA[<mrow><mi>s</mi><mo>.</mo><mi>t</mi><mo>.</mo><munderover><mo>&Sigma;</mo><mrow><mi>t</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msub><mi>&alpha;</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>&le;</mo><msub><mi>&alpha;</mi><mrow><mi>i</mi><mo>_</mo><mi>m</mi><mi>a</mi><mi>x</mi></mrow></msub></mrow>]]></math><img file="FDA0000743035740000046.GIF" wi="499" he="156" /></maths>   ⑧e:将上述的三个优化控制模型采用优化算法进行求解,最终获得对负荷的最优控制方法。
地址 100031 北京市西城区西长安街86号