发明名称 面向含有界不确定参数的家庭负荷调度方法
摘要 本发明涉及智能用电,为通过对时变电价机制的潜力挖掘以及对智能电器工作方式的合理调度,家庭能量管理不仅能够在满足用户舒适性的前提下,有效减少其用电费用,还能够通过对电价机制的响应,在整个系统起层面上起到削峰填谷的作用。本发明面向含有界不确定参数的家庭负荷调度方法,步骤是:101:构建含不确定性的家庭负荷模型;<maths num="0001"><math><![CDATA[ <mrow> <msubsup> <mi>&theta;</mi> <mrow> <mi>n</mi> <mo>+</mo> <mn>1</mn> </mrow> <mi>I</mi> </msubsup> <mo>=</mo> <mi>f</mi> <mrow> <mo>(</mo> <msubsup> <mi>&theta;</mi> <mi>n</mi> <mi>I</mi> </msubsup> <mo>,</mo> <msub> <mi>t</mi> <mi>n</mi> </msub> <mo>,</mo> <mi>Q</mi> <mo>,</mo> <mi>C</mi> <mo>,</mo> <mi>R</mi> <mo>,</mo> <msubsup> <mi>d</mi> <mi>n</mi> <mi>I</mi> </msubsup> <mo>,</mo> <msub> <mi>X</mi> <mrow> <mi>W</mi> <mi>H</mi> <mo>,</mo> <mi>n</mi> </mrow> </msub> <mo>,</mo> <msubsup> <mi>&theta;</mi> <mrow> <mi>e</mi> <mi>n</mi> </mrow> <mi>I</mi> </msubsup> <mo>)</mo> </mrow> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>3</mn> <mo>)</mo> </mrow> </mrow>]]></math><img file="DDA0001118767900000011.GIF" wi="941" he="70" /></maths>102:利用区间方程构建家庭负荷联合调度模型;103:对连续任务型负荷模型和能量控制优化决策模型进行求解。本发明主要应用于家庭等场合的智能用电。
申请公布号 CN106300368A 申请公布日期 2017.01.04
申请号 CN201610841369.4 申请日期 2016.09.22
申请人 天津大学 发明人 王继东;李银奇
分类号 H02J3/14(2006.01)I 主分类号 H02J3/14(2006.01)I
代理机构 天津市北洋有限责任专利代理事务所 12201 代理人 刘国威
主权项 一种面向含有界不确定参数的家庭负荷调度方法,其特征是,步骤如下:101:构建含不确定性的家庭负荷模型:负荷调度过程中所涉及到的设备包括电热水器,洗衣机,洗碗机,水泵,电动汽车,而其他部分设备被集中处理为固定负荷;当仅考虑热交换时所带来的热量损失时,电热水器的区间数热力学模型如式(1)所示:<maths num="0001"><math><![CDATA[<mrow><msubsup><mi>&theta;</mi><mrow><mi>n</mi><mo>+</mo><mn>1</mn></mrow><mi>I</mi></msubsup><mo>=</mo><msubsup><mi>&theta;</mi><mrow><mi>e</mi><mi>n</mi></mrow><mi>I</mi></msubsup><mo>+</mo><msub><mi>X</mi><mrow><mi>W</mi><mi>H</mi><mo>,</mo><mi>n</mi></mrow></msub><mi>Q</mi><mi>R</mi><mo>-</mo><mrow><mo>(</mo><msubsup><mi>&theta;</mi><mrow><mi>e</mi><mi>n</mi></mrow><mi>I</mi></msubsup><mo>+</mo><msub><mi>X</mi><mrow><mi>W</mi><mi>H</mi><mo>,</mo><mi>n</mi></mrow></msub><mi>Q</mi><mi>R</mi><mo>-</mo><msubsup><mi>&theta;</mi><mi>n</mi><mi>I</mi></msubsup><mo>)</mo></mrow><mi>exp</mi><mo>&lsqb;</mo><mo>-</mo><mrow><mo>(</mo><msub><mi>t</mi><mrow><mi>n</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>-</mo><msub><mi>t</mi><mi>n</mi></msub><mo>)</mo></mrow><mo>/</mo><mrow><mo>(</mo><mi>R</mi><mi>C</mi><mo>)</mo></mrow><mo>&rsqb;</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001118767870000011.GIF" wi="1454" he="70" /></maths>时间步长Δt=t<sub>n+1</sub>‑t<sub>n</sub>,<img file="FDA0001118767870000012.GIF" wi="377" he="71" />表示在n时刻热水的温度区间,X<sub>WH,n</sub>为电热水器在n时刻的运行状态,<img file="FDA0001118767870000013.GIF" wi="70" he="59" />表示环境温度,Q为电热水器的等效功率,R为其等效热阻,C为其等效热容。当仅考虑用户的用水行为时,电热水器的区间数热力学模型如式(2)所示<maths num="0002"><math><![CDATA[<mrow><msubsup><mi>&theta;</mi><mi>n</mi><mi>I</mi></msubsup><mo>=</mo><mo>&lsqb;</mo><msubsup><mi>&theta;</mi><mrow><mi>c</mi><mi>u</mi><mi>r</mi></mrow><mi>I</mi></msubsup><mrow><mo>(</mo><mi>M</mi><mo>-</mo><msubsup><mi>d</mi><mi>n</mi><mi>I</mi></msubsup><mo>)</mo></mrow><mo>+</mo><msubsup><mi>&theta;</mi><mrow><mi>e</mi><mi>n</mi></mrow><mi>I</mi></msubsup><msubsup><mi>d</mi><mi>n</mi><mi>I</mi></msubsup><mo>&rsqb;</mo><mo>/</mo><mi>M</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001118767870000014.GIF" wi="1430" he="66" /></maths>其中,<img file="FDA0001118767870000015.GIF" wi="444" he="68" />表示用水时刻热水的水温区间,<img file="FDA0001118767870000016.GIF" wi="379" he="70" />为用水量区间,M为水箱容量。综合(1)和(2)可得式(3):<maths num="0003"><math><![CDATA[<mrow><msubsup><mi>&theta;</mi><mrow><mi>n</mi><mo>+</mo><mn>1</mn></mrow><mi>I</mi></msubsup><mo>=</mo><mi>f</mi><mrow><mo>(</mo><msubsup><mi>&theta;</mi><mi>n</mi><mi>I</mi></msubsup><mo>,</mo><msub><mi>t</mi><mi>n</mi></msub><mo>,</mo><mi>Q</mi><mo>,</mo><mi>C</mi><mo>,</mo><mi>R</mi><mo>,</mo><msubsup><mi>d</mi><mi>n</mi><mi>I</mi></msubsup><mo>,</mo><msub><mi>X</mi><mrow><mi>W</mi><mi>H</mi><mo>,</mo><mi>n</mi></mrow></msub><mo>,</mo><msubsup><mi>&theta;</mi><mrow><mi>e</mi><mi>n</mi></mrow><mi>I</mi></msubsup><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001118767870000017.GIF" wi="1430" he="71" /></maths>102:利用区间方程构建家庭负荷联合调度模型:洗衣机,洗碗机,水泵,电动汽车四者所需要满足的约束如式(4)所示:<maths num="0004"><math><![CDATA[<mrow><munderover><mo>&Sigma;</mo><mrow><mi>n</mi><mo>=</mo><msub><mi>t</mi><mrow><mi>h</mi><mo>,</mo><mi>b</mi></mrow></msub></mrow><msub><mi>t</mi><mrow><mi>h</mi><mo>,</mo><mi>e</mi></mrow></msub></munderover><msub><mi>X</mi><mrow><mi>h</mi><mo>,</mo><mi>n</mi></mrow></msub><mo>=</mo><msub><mi>L</mi><mi>h</mi></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001118767870000018.GIF" wi="1433" he="146" /></maths>其中,t<sub>h,b</sub>和t<sub>h,e</sub>分别表示可中断负荷运行区间的起始时间和终止时间,x<sub>h,n</sub>表示设备在时间区间n时的开关状态,L<sub>h</sub>表示设备需要运行的时间区间总数;建立家庭能量管理模型,如式(5)‑式(9)所示,式(5)表示以总电费最小为目标的整体目标函数,式(6)和(7)表示电热水器的约束,为涉及区间数的不确定性约束;式(8)表示可中断负荷的运行约束,为确定性约束;式(9)为家庭需求限制所带来的约束;<maths num="0005"><math><![CDATA[<mrow><mtable><mtr><mtd><mrow><mi>M</mi><mi>i</mi><mi>n</mi></mrow></mtd><mtd><mrow><munderover><mo>&Sigma;</mo><mrow><mi>n</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><mrow><mo>(</mo><munderover><mo>&Sigma;</mo><mi>h</mi><mi>H</mi></munderover><msub><mi>p</mi><mi>n</mi></msub><mo>&CenterDot;</mo><msub><mi>X</mi><mrow><mi>h</mi><mo>,</mo><mi>n</mi></mrow></msub><mo>&CenterDot;</mo><msub><mi>P</mi><mrow><mi>s</mi><mi>t</mi><mi>e</mi><mi>p</mi><mo>,</mo><mi>h</mi></mrow></msub><mo>)</mo></mrow></mrow></mtd></mtr></mtable><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001118767870000019.GIF" wi="1507" he="131" /></maths><maths num="0006"><math><![CDATA[<mrow><mtable><mtr><mtd><mrow><mi>S</mi><mo>.</mo><mi>t</mi></mrow></mtd><mtd><mrow><msubsup><mi>&theta;</mi><mrow><mi>n</mi><mo>+</mo><mn>1</mn></mrow><mi>I</mi></msubsup><mo>=</mo><mi>f</mi><mrow><mo>(</mo><msubsup><mi>&theta;</mi><mi>n</mi><mi>I</mi></msubsup><mo>,</mo><msub><mi>t</mi><mi>n</mi></msub><mo>,</mo><mi>Q</mi><mo>,</mo><mi>C</mi><mo>,</mo><mi>R</mi><mo>,</mo><msub><mi>d</mi><mi>n</mi></msub><mo>,</mo><msub><mi>X</mi><mrow><mi>W</mi><mi>H</mi><mo>,</mo><mi>n</mi></mrow></msub><mo>,</mo><msub><mi>&theta;</mi><mrow><mi>e</mi><mi>n</mi></mrow></msub><mo>)</mo></mrow></mrow></mtd></mtr></mtable><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA00011187678700000110.GIF" wi="1510" he="69" /></maths><maths num="0007"><math><![CDATA[<mrow><msubsup><mi>&theta;</mi><mi>n</mi><mrow><mi>d</mi><mi>o</mi><mi>w</mi><mi>n</mi></mrow></msubsup><mo>&le;</mo><msubsup><mi>&theta;</mi><mi>n</mi><mi>I</mi></msubsup><mo>&le;</mo><msubsup><mi>&theta;</mi><mi>n</mi><mrow><mi>u</mi><mi>p</mi></mrow></msubsup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>7</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA00011187678700000111.GIF" wi="1326" he="63" /></maths><maths num="0008"><math><![CDATA[<mrow><munderover><mo>&Sigma;</mo><mrow><mi>n</mi><mo>=</mo><msub><mi>t</mi><mrow><mi>h</mi><mo>,</mo><mi>b</mi></mrow></msub></mrow><msub><mi>t</mi><mrow><mi>h</mi><mo>,</mo><mi>e</mi></mrow></msub></munderover><msub><mi>x</mi><mrow><mi>h</mi><mo>,</mo><mi>n</mi></mrow></msub><mo>=</mo><msub><mi>L</mi><mi>h</mi></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>8</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA00011187678700000112.GIF" wi="1341" he="144" /></maths><maths num="0009"><math><![CDATA[<mrow><munderover><mo>&Sigma;</mo><mi>h</mi><mi>H</mi></munderover><msub><mi>P</mi><mi>h</mi></msub><mo>&CenterDot;</mo><msub><mi>X</mi><mrow><mi>h</mi><mo>,</mo><mi>n</mi></mrow></msub><mo>+</mo><msub><mi>P</mi><mrow><mi>m</mi><mi>u</mi><mi>s</mi><mi>t</mi><mi>r</mi><mi>u</mi><mi>n</mi><mo>,</mo><mi>n</mi></mrow></msub><mo>&le;</mo><msub><mi>P</mi><mrow><mi>d</mi><mi>e</mi><mi>m</mi><mi>a</mi><mi>n</mi><mi>d</mi></mrow></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>9</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA00011187678700000113.GIF" wi="1350" he="127" /></maths>其中,p<sub>n</sub>为实时电价,P<sub>h</sub>为各个负荷的功率,P<sub>step,h</sub>为各个负荷的在一个步长之内消耗的能量,H为负荷集合,P<sub>mustrun,n</sub>为固定负荷在第n个时间段内的能耗;103:对连续任务型负荷模型和能量控制优化决策模型进行求解定义区间数之间定量比较的新的数学工具:区间可能度P(a≤B<sup>I</sup>)如式(10)所示:<maths num="0010"><math><![CDATA[<mrow><mi>P</mi><mrow><mo>(</mo><mi>a</mi><mo>&le;</mo><msup><mi>B</mi><mi>I</mi></msup><mo>)</mo></mrow><mo>=</mo><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><mn>1</mn><mo>,</mo><mi>a</mi><mo>&le;</mo><msup><mi>B</mi><mrow><mi>d</mi><mi>o</mi><mi>w</mi><mi>n</mi></mrow></msup></mrow></mtd></mtr><mtr><mtd><mrow><mfrac><mrow><msup><mi>B</mi><mrow><mi>u</mi><mi>p</mi></mrow></msup><mo>-</mo><mi>a</mi></mrow><mrow><msup><mi>B</mi><mrow><mi>u</mi><mi>p</mi></mrow></msup><mo>-</mo><msup><mi>B</mi><mrow><mi>d</mi><mi>o</mi><mi>w</mi><mi>n</mi></mrow></msup></mrow></mfrac><mo>,</mo><msup><mi>B</mi><mrow><mi>d</mi><mi>o</mi><mi>w</mi><mi>n</mi></mrow></msup><mo>&lt;</mo><mi>a</mi><mo>&le;</mo><msup><mi>B</mi><mrow><mi>u</mi><mi>p</mi></mrow></msup></mrow></mtd></mtr><mtr><mtd><mrow><mn>0</mn><mo>,</mo><mi>a</mi><mo>&gt;</mo><msup><mi>B</mi><mrow><mi>u</mi><mi>p</mi></mrow></msup></mrow></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>10</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001118767870000021.GIF" wi="1348" he="310" /></maths>其中a表示一个实数,B<sup>I</sup>为一个区间数,将式(7)中的一个约束转换为两个等价的不等式约束,如式(11)所示:<maths num="0011"><math><![CDATA[<mrow><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><mi>P</mi><mrow><mo>(</mo><msubsup><mi>&theta;</mi><mi>n</mi><mi>I</mi></msubsup><mo>&GreaterEqual;</mo><msubsup><mi>&theta;</mi><mi>n</mi><mrow><mi>d</mi><mi>o</mi><mi>w</mi><mi>n</mi></mrow></msubsup><mo>)</mo></mrow><mo>&GreaterEqual;</mo><mn>1</mn><mo>-</mo><msub><mi>&sigma;</mi><mn>1</mn></msub></mrow></mtd></mtr><mtr><mtd><mrow><mi>P</mi><mrow><mo>(</mo><msubsup><mi>&theta;</mi><mi>n</mi><mi>I</mi></msubsup><mo>&GreaterEqual;</mo><msubsup><mi>&theta;</mi><mi>n</mi><mrow><mi>u</mi><mi>p</mi></mrow></msubsup><mo>)</mo></mrow><mo>&GreaterEqual;</mo><msub><mi>&sigma;</mi><mn>2</mn></msub></mrow></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>11</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001118767870000022.GIF" wi="1317" he="157" /></maths>其中σ表示用户的约束违背容忍度,σ<sub>1</sub>和σ<sub>2</sub>分别表示用户对于水温约束下限违背和上限违背的容忍度,两者的值越大,代表用户能够容忍越大的约束触犯,反之亦然;采用二进制粒子群算法和整数规划的联合求解方案。
地址 300072 天津市南开区卫津路92号