CN104112168B - 一种基于多代理***的智能家居优化方法 - Google Patents
一种基于多代理***的智能家居优化方法 Download PDFInfo
- Publication number
- CN104112168B CN104112168B CN201410267513.9A CN201410267513A CN104112168B CN 104112168 B CN104112168 B CN 104112168B CN 201410267513 A CN201410267513 A CN 201410267513A CN 104112168 B CN104112168 B CN 104112168B
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- mfrac
- msup
- mtd
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 27
- 238000005457 optimization Methods 0.000 title claims abstract description 19
- 230000005611 electricity Effects 0.000 claims abstract description 70
- 238000010248 power generation Methods 0.000 claims description 47
- 238000004146 energy storage Methods 0.000 claims description 24
- 238000004364 calculation method Methods 0.000 claims description 19
- 238000009826 distribution Methods 0.000 claims description 19
- 230000006870 function Effects 0.000 claims description 13
- 238000004422 calculation algorithm Methods 0.000 claims description 10
- 238000011156 evaluation Methods 0.000 claims description 10
- 230000002068 genetic effect Effects 0.000 claims description 9
- 238000007599 discharging Methods 0.000 claims description 8
- 230000008569 process Effects 0.000 claims description 8
- 238000005286 illumination Methods 0.000 claims description 6
- 238000005520 cutting process Methods 0.000 claims description 5
- 238000010276 construction Methods 0.000 claims description 4
- 230000004044 response Effects 0.000 claims description 4
- 238000006243 chemical reaction Methods 0.000 claims description 3
- 238000004519 manufacturing process Methods 0.000 abstract description 3
- 230000007613 environmental effect Effects 0.000 abstract 1
- 238000003860 storage Methods 0.000 abstract 1
- 239000003795 chemical substances by application Substances 0.000 description 24
- 230000006399 behavior Effects 0.000 description 3
- 238000009499 grossing Methods 0.000 description 3
- 230000003993 interaction Effects 0.000 description 3
- 230000003466 anti-cipated effect Effects 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 239000000446 fuel Substances 0.000 description 2
- 230000000694 effects Effects 0.000 description 1
- 238000005265 energy consumption Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 230000008450 motivation Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Secondary Cells (AREA)
Abstract
本发明公开了一种基于多代理***的智能家居优化方法,该方法考虑了动态的环境因素,把每个用户考虑为一个自动代理,以所有代理组成一个***,基于目前和预测的用电负荷总量、生产电能总量、所储存的电能总量、以及相关电价和用户满意度函数进行合理的出售、购买或者储存电能,得到***经济性最优的方案。本发明能够有效地最小化居民用户的电费,同时能够使总负荷需求曲线更为平滑。
Description
技术领域
本发明涉及一种基于多代理***的智能家居优化方法,属于电力***技术领域。
背景技术
智能电网的一个内涵就是要求把用户的积极性调动起来,实现用户与电网的双向互动,最终达到削峰填谷、提高能源利用率的目的。
电力***作为主要的电能提供和结构,必须不停的适应新科技和新概念来改进电能的适用效率。
电力提供方必须不停的研究电能使用方的需求和行为,然后开发可靠,合理,效率的***。由于大量的电力用户分散在***中,很难总体研究消费用户的用电行为,所以采用模拟的方法来研究和预测用电负荷大小。
随着智能电网和智能分布式***的产生,许多目前的研究主要集中于模拟智能家居和电网的互动,以及基于需求响应的负荷管理。
分布式电源和储能***能够使得居家和小企业通过向电网售卖多余的电能获得利润。通过智能自动管理设备,用户能够基于预测的负荷需求总量、产能总量、和电价等做出合理的决定,购电、售电、或储能。
传统的家居设备不具有小型生产电能设备(风能、太阳能***),和储能设备,也不能向电网售电,但是随着智能电网的发展,智能家居的概念随之产生,凭借智能设备、监控设备和高级量测***的开发,能够较好的与电网进行互动。
智能家居,作为与智能电网互动的主动参与者,不再被简单的视为负荷曲线,由于目前他们能够很好的和电网进行交互,能够影响电力需求,电能生产,和电力价格。
智能家居不仅能消费电能,而且能够产生电能,储存起来或者自己使用。这些智能家居用户能够智能合理的做出自动决策来管理自己的负荷,生产的电能和电力储能,同时和电网交互,降低自己的消费,同时能达到消峰填谷的作用,平滑负荷曲线。
发明内容
本发明提供一种基于多代理***的智能家居优化方法,以每个居家用户为一个代理,以所有代理组成的多代理作为一个***,利用目前和预测的家电用电总量和发电总量以及储能总量进行自动智能处理,最小化家居用户分电费费用。
本发明具体采用以下技术方案实现:
一种基于多代理***的智能家居优化方法,包括以下步骤:
1)建立具有随机性出力的风力发电***出力模型,太阳能光伏发电***出力模型和电动汽车充放电***出力模型;
2)计算居家用户满意度,具体包括以下步骤:
2-1)计算居家用电需求反馈前的电价P(h-)和需求响应后的实际电价P(h):
其中:kτ是电价的权重系数,代表了目前和τ天前电价之间的相关程度,
其中:∝1·、∝2·、∝3·、∝4·是相关系数,代表了平均每个居家用户在h小时内的实际的负荷需求;
2-2)计算居家用户的售电倾向Sell(h),
其中,gi(h)为风力发电***生产的电能总量,si(h)为太阳能发电***生产的电能总量,li(h)为电动汽车充电负荷,P(h-)表示用电需求反馈前的电价,max(P(h-))表示用电需求反馈前的最大电价,li(h-)为储能***的电能总量;
2-3)计算居家用户的储能倾向Store(h),表达式如下:
2-4)计算居家用户满意度计算公式如下:
其中,表示在h时段内设备a消耗的总电能,Eh为用户在h时段内储存的电能,Stores(h)是储能意向判断边界条件数值,Sells(h)是售电意向判断边界条件数值,n为所有的居家用户;
3)以每个居家用户为1个代理,所有用户作为1个多代理***,以***最小化电费为目标,以居家用户满意度最大为目标函数,采用遗传算法进行整个***的经济性优化,根据评判标准值,做出相应的操作。
前述的步骤1)的随机性出力模型,是指根据所选取的时间段内的历史数据,生成随机模型的均值和方差。
前述的步骤1)风力发电***出力模型的构建过程为:
1-11)根据历史风速数据,计算风力发电***的风速;
1-12)计算风力发电***风速的概率密度f(v),计算公式如下:
其中:v为风速,k和c为威布尔分布的两个参数;
1-13)计算风力发电***的风机出力Pw,计算公式如下:
其中:k2=-k1·vci,Pr为风机额定功率,vci为切入风速,vr为额定风速,vco为切出风速;
1-14)计算风力发电***的风机出力的概率密度f(Pw),计算公式如下:
前述的步骤1)太阳能光伏发电***出力模型的构建方法为:
1-21)根据历史太阳能光照强度数据,计算太阳能光伏发电***的日照强度;
1-22)计算太阳能光伏发电***日照强度的概率密度f(r),计算公式如下:
其中:r和rmax分别是计算时间段内的实际光强和最大光强,α和β为Beta分布的形状参数,Γ为Gamma函数;
1-23)计算太阳能光伏发电***的输出功率PM,计算公式如下:
PM=r·A·η
其中:A和η分别为太阳能光伏发电***的面积和光电转换效率;
1-24)计算太阳能光伏发电***输出功率的概率密度f(PM),计算公式如下:
其中:RM表示太阳能光伏发电***输出的最大输出功率。
前述的步骤1)电动汽车充放电***出力模型的构建过程为:
1-31)根据历史电动汽车充电负荷数据,计算电动汽车充电功率的均值和方差;
1-32)电动汽车类型kev,电池的剩余电量Cbat服从二维正太分布:
其中,X表示变量,μ表示电动汽车充电功率的均值,C电池的剩余电量参数;
电动汽车充电功率PCI服从正太分布:
其中,d表示充电效率,σpci表示充电功率的方差;
电动汽车充电时刻T服从正太分布:
其中,μt为电动汽车充电时刻均值;
电动汽车日行驶路程M服从正太分布:
其中,dd表示日行驶里程数损耗;σM表示日行驶路程数方差;
前述的步骤3)中的操作包括:Ⅰ)从电网购电;Ⅱ)对储能电池进行充电或者放电;Ⅲ)向电网售电。
前述的做出相应的操作的评判标准为:
Ⅰ)从电网购电的评断标准为:
且Store(h)-Stores(h)>0或者
Ⅱ-1)储能电池进行充电的评断标准为:
且Store(h)-Stores(h)>0
Ⅱ-2)储能电池进行放电的评断标准为:
且Store(h)-Stores(h)>0或者
Ⅲ)向电网售电的评断标准为:
且Store(h)-Stores(h)>0。
前述的步骤3)中采用遗传算法进行整个***的经济性优化的具体过程如下:
(a)对参数Store(h)、Sell(h)、Eh进行采用遗传算法的二进制编码;
(b)初始化种群迭代次数t1=0,个体迭代次数t2=0,计数参数n=0,种群数组Node,设置种群大小为Nall,种群和个体最大迭代次数T1,T2;
(c)生成初始种群;
(d)计算种群中每个个体的适应度F和聚集距离P;
(e)采用选择操作从种群中选取n1个个体,并对聚集距离从大到小进行排序,选取前n2个个体,其中,n1和n2满足n1+n2=Nall;
(f)对种群进行遗传和变异操作;
(g)重复步骤(d)~(f)直至达到最大迭代次数;
(h)输出最优解并退出。
本发明利用基于JAVA语言的Repast Simphony软件行建模,建立了智能家居代理***,定义了若干智能家居代理***进行自动选择方案的公式及概念,即售电意愿,购电意愿,储能意愿,通过目前的生产电能总量,电价和预测的生产电能总量,负荷需求总量,预测电价进行计算,根据判据进行方案选择,得到售电,买电或储能的方案,能够最大化的减少居家用户的用电费用,同时能够消峰填谷,平滑负荷曲线。
附图说明
图1为本发明的智能家居优化方法的流程图;
图2为采用遗传算法进行整个***的经济性优化的流程图。
具体实施方式
下面对本发明的技术方案进行详细说明:
本发明把智能电网中的每一个用电用户看成一个代理,利用目前和预测的家电用电总量、发电总量、储能总量以及用户对用电满意度情况进行自动智能处理,最小化家居用户分电费费用,消峰填谷,平滑负荷曲线。
每一个代理具有各种用电设备,包括产生电能的分布式电源(风力发电***、太阳能发电***),使用电能的设备(空调、电视等)以及储存电能的设备(储能电池等)。
本发明采用基于JAVA语言的Repast Simphony软件进行建模,由于输入模型的数据(风力发电***、负荷预测值,预测电价值)具有较大的随机性,智能家居智能代理***自动控制的方案受这些输入变量的影响较大,为了使得操作方案更为精确,更好的减少居家用户的电费,为此,本发明建立了描述随机性的模型,即根据所选取的时间段内的历史数据,生成随机模型的基本参数(均值和方差)以及确定其概率分布。
根据居家用户发电总量和用电总量比较得到储能意向的大小,该值越大,表明居家用户越有意向储存电能。智能家居代理想要储存电能,在未来产生足够的电能以一个较高的价格卖给电网。如果预期生产的电能小于预期的负荷大小,电能不足或电价较高的时,智能用户代理将会更愿意储存电能。
如图1所示,本发明包括以下步骤:
一、首先建立具有随机性出力的风力发电***出力模型,太阳能光伏发电***出力模型和电动汽车充放电***出力模型;
1)风力发电***
根据历史风速数据,计算风力发电***的风速,得到风速的概率密度函数如下:
其中,历史风速数据可由当地的风力发电***管理单位获得,
式中:v为风速,k和c为威布尔分布的两个参数,可以由平均风速μ和标准差σ近似算出,
风机的出力大小Pw可由下式得到:
其中:k2=-k1·vci,Pr为风机额定功率,vci为切入风速,vr为额定风速,vco为切出风速。
进而可以得到风力发电***的风机出力的概率密度f(Pw):
根据风机出力的概率密度f(Pw)可以计算出风力发电***生产的电能总量。
2)太阳能发电***
根据历史太阳能光照强度数据,计算太阳能光伏发电***的日照强度,然后得到太阳能光照强度的概率密度函数f(r)如下:
其中,历史太阳能光照强度数据,可由当地的太阳能发电***管理单位获得,
式中:r和rmax(W/m2)分别是计算时间段内的实际光强和最大光强,α和β都是Beta分布的形状参数;Γ为Gamma函数。
Beta分布的参数可由该段时间内的光照强度平均值μ和方差σ得到,关系如下:
太阳能光伏发电***的输出功率PM为:
PM=r·A·η
式中:A和η分别为太阳能光伏发电***的面积和光电转换效率。
进而可以得到太阳能光伏发电***输出功率的概率密度f(PM):
其中,RM表示太阳能光伏发电***输出的最大输出功率,
根据f(PM)可以计算出太阳能发电***生产的电能总量。
3)电动汽车充放电***出力模型
电动汽车的随机模型主要受电动汽车类型、充电方式、用户行为因素。
汽车的类型kev和电池的剩余电量Cbat存在关系:
其中,CFuel为燃料剩余量;
电动汽车类型kev,电池的剩余电量Cbat服从二维正太分布:
其中,X表示kev,Cbat这两个变量,μ表示电动汽车充电功率的均值,C电池的剩余电量参数。
电动汽车充电功率PCI服从正太分布:
其中,d表示充电效率,σpci表示充电功率的方差。
电动汽车充电时刻T服从正太分布:
其中,μt为电动汽车充电时刻均值;
电动汽车日行驶路程M服从正太分布:
其中,dd表示日行驶里程数损耗;σM表示日行驶路程数方差。
根据f(PCI)可以计算得到电动汽车充电负荷。
二、目标函数的计算,目标函数指的是居家用户满意度。
本发明基于实时电价、设备的运行状态(储存的电能大小)和未来的电价建立用户的满意度评价函数:
1)计算居家用电需求反馈前的电价P(h-)和需求响应后的实际电价P(h)
其中:kτ是电价的权重系数,代表了目前和τ天前电价之间的相关程度,
其中:∝1·、∝2·、∝3·、∝4·是相关系数,代表了平均每个居家用户在h小时内的实际的负荷需求。
2)售电倾向Sell(h)的计算
售电倾向代表了居家代理出售自己生产的额外电能的意愿大小,其表达式如下所示:
式中:gi(h)为风力发电***生产的电能总量,si(h)为太阳能发电***生产的电能总量,li(h)为电动汽车充电负荷,P(h-)表示用电需求反馈前的电价,max(P(h-))表示用电需求反馈前的最大电价,li(h-)为储能***的电能总量;
如果目前时间内居家的生产的电能大小大于负荷需求,那么,对于居家用户来说,就有更高的意愿来出售电力给电网,应为目前居家有较多的电能剩余,或者出售的电价比预测的高。
为了计算Sell(h)智能家居必须获得目前时刻的售电电价P(h-),并且使用概率模型来预测需要的变量,越高的Sell(h)表明,在该时刻售电给电网,居家用户能够得到更多的好处。
3)计算居家用户的储能倾向Store(h),
储能倾向表达了居家用户代理储存额外电能的意愿大小,其表达式如下所示:
智能家居代理想要储存电能,然后在未来产生足够的电能以一个较高的价格出售给电网。如果预期的产能小于预期的负荷大小,代理将会更愿意储存电能。如果他们预测未来有一个较大的电能不足时刻或较高的电价,同样的,代理将会更愿意储存电能。
4)用户满意度函数的计算
用户根据负荷需求预测值和电价预测值在h时段储存的电能为Eh,这部分的电能会产生预先储存电能费用P(h-)*Eh,在实际调度阶段,该部分的电能会以实时电价提供给使用者。当储存的电能Eh不能满足用电需求时,就需要从电力市场购买额外的电能来满足需求,这将会产生额外的用以平衡负荷的费用
满意度函数表示用户对于该时段调度的满意程度,
其中,表示在h时段内设备a消耗的总电能,Eh为用户在h时段内储存的电能,Stores(h)是储能意向判断边界条件数值,Sells(h)是售电意向判断边界条件数值,n为所有的居家用户;Eh通过计算机对设备能耗进行记录。
三、以每个居家用户为1个代理,所有用户作为1个多代理***,以***最小化电费为目标,以居家用户满意度最大为目标函数,采用遗传算法进行整个***的经济性优化,根据评判标准值,做出相应的操作。
其中,所述操作包括:Ⅰ)从电网购电;Ⅱ)对储能电池进行充电或者放电;Ⅲ)向电网售电;
上述操作的是在目标函数最优化的情况下得到的参变量(Eh、Sell(h)、Store(h),然后根据下述评判标准进行的:
表1 操作评判标准
如图2所示,采用遗传算法进行整个***的经济性优化,具体过程如下:
(a)对参数Store(h)、Sell(h)、Eh进行采用遗传算法的二进制编码;
(b)初始化种群迭代次数t1=0,个体迭代次数t2=0,计数参数n=0,种群数组Node,设置种群大小为Nall,种群和个体最大迭代次数T1,T2;
(c)生成初始种群;
(d)计算种群中每个个体的适应度F和聚集距离P;
(e)采用选择操作从种群中选取n1个个体,并对聚集距离从大到小进行排序,选取前n2个个体,其中,n1和n2满足n1+n2=Nall;
(f)对种群进行遗传和变异操作;
(g)重复步骤(d)~(f)直至达到最大迭代次数;
(h)输出最优解并退出。
Claims (7)
1.一种基于多代理***的智能家居优化方法,其特征在于,包括以下步骤:
1)建立具有随机性出力的风力发电***出力模型,太阳能光伏发电***出力模型和电动汽车充放电***出力模型;
2)计算居家用户满意度,具体包括以下步骤:
2-1)计算居家用电需求反馈前的电价P(h-)和需求响应后的实际电价P(h):
<mrow>
<mi>P</mi>
<mrow>
<mo>(</mo>
<mi>h</mi>
<mo>-</mo>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msubsup>
<mi>&Sigma;</mi>
<mrow>
<mi>&tau;</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<mi>m</mi>
</msubsup>
<msup>
<mi>k</mi>
<mi>&tau;</mi>
</msup>
<mo>&CenterDot;</mo>
<mi>P</mi>
<mrow>
<mo>(</mo>
<mi>h</mi>
<mo>-</mo>
<mi>&tau;</mi>
<mo>)</mo>
</mrow>
<mo>,</mo>
<msubsup>
<mi>&Sigma;</mi>
<mrow>
<mi>&tau;</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<mi>m</mi>
</msubsup>
<msup>
<mi>k</mi>
<mi>&tau;</mi>
</msup>
<mo>=</mo>
<mn>1</mn>
</mrow>
其中:kτ是电价的权重系数,代表了目前和τ天前电价之间的相关程度,
<mrow>
<mi>P</mi>
<mrow>
<mo>(</mo>
<mi>h</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mo>&Proportional;</mo>
<mrow>
<mn>1</mn>
<mo>&CenterDot;</mo>
</mrow>
</msub>
<msup>
<mi>e</mi>
<mrow>
<msub>
<mo>&Proportional;</mo>
<mrow>
<mn>2</mn>
<mo>&CenterDot;</mo>
</mrow>
</msub>
<mover>
<mi>I</mi>
<mo>&OverBar;</mo>
</mover>
<mrow>
<mo>(</mo>
<mi>h</mi>
<mo>)</mo>
</mrow>
</mrow>
</msup>
<mo>+</mo>
<msub>
<mo>&Proportional;</mo>
<mrow>
<mn>3</mn>
<mo>&CenterDot;</mo>
</mrow>
</msub>
<msup>
<mi>e</mi>
<mrow>
<msub>
<mo>&Proportional;</mo>
<mrow>
<mn>4</mn>
<mo>&CenterDot;</mo>
</mrow>
</msub>
<mover>
<mi>I</mi>
<mo>&OverBar;</mo>
</mover>
<mrow>
<mo>(</mo>
<mi>h</mi>
<mo>)</mo>
</mrow>
</mrow>
</msup>
</mrow>
其中:∝1·、∝2·、∝3.、∝4.是相关系数,代表了平均每个居家用户在h小时内的实际的负荷需求;
2-2)计算居家用户的售电倾向Sell(h),
<mrow>
<mi>S</mi>
<mi>e</mi>
<mi>l</mi>
<mi>l</mi>
<mrow>
<mo>(</mo>
<mi>h</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mi>P</mi>
<mrow>
<mo>(</mo>
<mi>h</mi>
<mo>-</mo>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
<mrow>
<mo>(</mo>
<mi>P</mi>
<mo>(</mo>
<mrow>
<mi>h</mi>
<mo>-</mo>
</mrow>
<mo>)</mo>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>&CenterDot;</mo>
<mfrac>
<mrow>
<mo>(</mo>
<mo>(</mo>
<mrow>
<msub>
<mi>g</mi>
<mi>i</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>h</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>s</mi>
<mi>i</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>h</mi>
<mo>)</mo>
</mrow>
</mrow>
<mo>)</mo>
<mo>-</mo>
<msub>
<mi>l</mi>
<mi>i</mi>
</msub>
<mo>(</mo>
<mrow>
<mi>h</mi>
<mo>-</mo>
</mrow>
<mo>)</mo>
<mo>)</mo>
<mo>&CenterDot;</mo>
<mi>P</mi>
<mo>(</mo>
<mi>h</mi>
<mo>-</mo>
<mo>)</mo>
</mrow>
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
<mrow>
<mo>(</mo>
<mo>(</mo>
<mrow>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>g</mi>
<mi>i</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>h</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>s</mi>
<mi>i</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>h</mi>
<mo>)</mo>
</mrow>
</mrow>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mi>l</mi>
<mi>i</mi>
</msub>
<mrow>
<mo>(</mo>
<mrow>
<mi>h</mi>
<mo>-</mo>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
<mo>)</mo>
<mo>&CenterDot;</mo>
<mi>P</mi>
<mo>(</mo>
<mrow>
<mi>h</mi>
<mo>-</mo>
</mrow>
<mo>)</mo>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
</msup>
<mn>...............</mn>
<mrow>
<mo>(</mo>
<msub>
<mi>g</mi>
<mi>i</mi>
</msub>
<mo>(</mo>
<mi>h</mi>
<mo>)</mo>
<mo>+</mo>
<msub>
<mi>s</mi>
<mi>i</mi>
</msub>
<mo>(</mo>
<mi>h</mi>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mo>></mo>
<msub>
<mi>l</mi>
<mi>i</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>h</mi>
<mo>-</mo>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mfrac>
<mrow>
<mi>P</mi>
<mrow>
<mo>(</mo>
<mi>h</mi>
<mo>-</mo>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
<mrow>
<mo>(</mo>
<mi>P</mi>
<mo>(</mo>
<mrow>
<mi>h</mi>
<mo>-</mo>
</mrow>
<mo>)</mo>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mn>.................................................................................</mn>
<mrow>
<mo>(</mo>
<msub>
<mi>g</mi>
<mi>i</mi>
</msub>
<mo>(</mo>
<mi>h</mi>
<mo>)</mo>
<mo>+</mo>
<msub>
<mi>s</mi>
<mi>i</mi>
</msub>
<mo>(</mo>
<mi>h</mi>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mo><</mo>
<msub>
<mi>l</mi>
<mi>i</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>h</mi>
<mo>-</mo>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
其中,gi(h)为风力发电***生产的电能总量,si(h)为太阳能发电***生产的电能总量,li(h)为电动汽车充电负荷,P(h-)表示用电需求反馈前的电价,max(P(h-))表示用电需求反馈前的最大电价,li(h-)为储能***的电能总量;
2-3)计算居家用户的储能倾向Store(h),表达式如下:
<mrow>
<mi>S</mi>
<mi>t</mi>
<mi>o</mi>
<mi>r</mi>
<mi>e</mi>
<mrow>
<mo>(</mo>
<mi>h</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mfrac>
<mrow>
<mi>A</mi>
<mi>v</mi>
<mi>e</mi>
<mi>r</mi>
<mi>a</mi>
<mi>g</mi>
<mi>e</mi>
<mrow>
<mo>(</mo>
<mi>P</mi>
<mo>(</mo>
<mi>h</mi>
<mo>)</mo>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>max</mi>
<mrow>
<mo>(</mo>
<mi>P</mi>
<mo>(</mo>
<mi>h</mi>
<mo>)</mo>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mn>..................................................................................</mn>
<mrow>
<mo>(</mo>
<msub>
<mi>g</mi>
<mi>i</mi>
</msub>
<mo>(</mo>
<mi>h</mi>
<mo>)</mo>
<mo>+</mo>
<msub>
<mi>s</mi>
<mi>i</mi>
</msub>
<mo>(</mo>
<mi>h</mi>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mo>></mo>
<msub>
<mi>l</mi>
<mi>i</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>h</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mi>A</mi>
<mi>v</mi>
<mi>e</mi>
<mi>r</mi>
<mi>a</mi>
<mi>g</mi>
<mi>e</mi>
<mrow>
<mo>(</mo>
<mi>P</mi>
<mo>(</mo>
<mi>h</mi>
<mo>)</mo>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>max</mi>
<mrow>
<mo>(</mo>
<mi>P</mi>
<mo>(</mo>
<mi>h</mi>
<mo>)</mo>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>&CenterDot;</mo>
<mfrac>
<mrow>
<mi>A</mi>
<mi>v</mi>
<mi>e</mi>
<mi>r</mi>
<mi>a</mi>
<mi>g</mi>
<mi>e</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>l</mi>
<mi>i</mi>
</msub>
<mo>(</mo>
<mi>h</mi>
<mo>)</mo>
<mo>-</mo>
<mo>(</mo>
<mrow>
<msub>
<mi>g</mi>
<mi>i</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>h</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>s</mi>
<mi>i</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>h</mi>
<mo>)</mo>
</mrow>
</mrow>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mo>&CenterDot;</mo>
<mrow>
<mo>(</mo>
<mi>P</mi>
<mo>(</mo>
<mi>h</mi>
<mo>)</mo>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>max</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>l</mi>
<mi>i</mi>
</msub>
<mo>(</mo>
<mi>h</mi>
<mo>)</mo>
<mo>-</mo>
<mo>(</mo>
<mrow>
<msub>
<mi>g</mi>
<mi>i</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>h</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>s</mi>
<mi>i</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>h</mi>
<mo>)</mo>
</mrow>
</mrow>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mo>&CenterDot;</mo>
<mrow>
<mo>(</mo>
<mi>P</mi>
<mo>(</mo>
<mi>h</mi>
<mo>)</mo>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
</msup>
<mn>...............</mn>
<mrow>
<mo>(</mo>
<msub>
<mi>g</mi>
<mi>i</mi>
</msub>
<mo>(</mo>
<mi>h</mi>
<mo>)</mo>
<mo>+</mo>
<msub>
<mi>s</mi>
<mi>i</mi>
</msub>
<mo>(</mo>
<mi>h</mi>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mo><</mo>
<msub>
<mi>l</mi>
<mi>i</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>h</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
2-4)计算居家用户满意度计算公式如下:
<mrow>
<msup>
<mi>W</mi>
<mi>h</mi>
</msup>
<mrow>
<mo>(</mo>
<msubsup>
<mi>X</mi>
<mi>a</mi>
<mi>h</mi>
</msubsup>
<mo>,</mo>
<msub>
<mi>E</mi>
<mi>h</mi>
</msub>
<mo>,</mo>
<mi>S</mi>
<mi>e</mi>
<mi>l</mi>
<mi>l</mi>
<mo>(</mo>
<mi>h</mi>
<mo>)</mo>
<mo>,</mo>
<mi>S</mi>
<mi>t</mi>
<mi>o</mi>
<mi>r</mi>
<mi>e</mi>
<mo>(</mo>
<mi>h</mi>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msubsup>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</msubsup>
<mrow>
<mo>(</mo>
<mo>|</mo>
<mrow>
<mi>S</mi>
<mi>e</mi>
<mi>l</mi>
<mi>l</mi>
<mrow>
<mo>(</mo>
<mi>h</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mi>Sell</mi>
<mi>s</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>h</mi>
<mo>)</mo>
</mrow>
</mrow>
<mo>|</mo>
<mo>+</mo>
<mo>|</mo>
<mrow>
<mi>S</mi>
<mi>t</mi>
<mi>o</mi>
<mi>r</mi>
<mi>e</mi>
<mrow>
<mo>(</mo>
<mi>h</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mi>Store</mi>
<mi>s</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>h</mi>
<mo>)</mo>
</mrow>
</mrow>
<mo>|</mo>
<mo>-</mo>
<mi>P</mi>
<mo>(</mo>
<mi>h</mi>
<mo>)</mo>
<mo>*</mo>
<mo>&lsqb;</mo>
<msub>
<mi>&Sigma;</mi>
<mrow>
<mi>a</mi>
<mo>&Element;</mo>
<mi>Z</mi>
</mrow>
</msub>
<msubsup>
<mi>X</mi>
<mi>a</mi>
<mi>h</mi>
</msubsup>
<mo>-</mo>
<msub>
<mi>E</mi>
<mi>h</mi>
</msub>
<mo>&rsqb;</mo>
<mo>-</mo>
<mi>P</mi>
<mo>(</mo>
<mrow>
<mi>h</mi>
<mo>-</mo>
</mrow>
<mo>)</mo>
<mo>*</mo>
<msub>
<mi>E</mi>
<mi>h</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
其中,表示在h时段内设备a消耗的总电能,Eh为用户在h时段内储存的电能,Stores(h)是储能意向判断边界条件数值,Sells(h)是售电意向判断边界条件数值,n为所有的居家用户;
3)以每个居家用户为1个代理,所有用户作为1个多代理***,以***最小化电费为目标,以居家用户满意度最大为目标函数,采用遗传算法进行整个***的经济性优化,根据评判标准值,做出相应的操作;
所述做出相应的操作的评判标准为:
Ⅰ)从电网购电的评断标准为:
且Store(h)-Stores(h)>0或者
Ⅱ-1)储能电池进行充电的评断标准为:
且Store(h)-Stores(h)>0
Ⅱ-2)储能电池进行放电的评断标准为:
且Store(h)-Stores(h)>0或者
Ⅲ)向电网售电的评断标准为:
且Store(h)-Stores(h)>0。
2.根据权利要求1所述的一种基于多代理***的智能家居优化方法,其特征在于,所述步骤1)的三种随机性出力模型,是指根据所选取的时间段内的历史数据,生成随机模型的均值和方差。
3.根据权利要求1所述的一种基于多代理***的智能家居优化方法,其特征在于,所述步骤1)风力发电***出力模型的构建过程为:
1-11)根据历史风速数据,计算风力发电***的风速;
1-12)计算风力发电***风速的概率密度f(v),计算公式如下:
<mrow>
<mi>f</mi>
<mrow>
<mo>(</mo>
<mi>v</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mi>k</mi>
<mi>c</mi>
</mfrac>
<mo>&CenterDot;</mo>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<mi>v</mi>
<mi>c</mi>
</mfrac>
<mo>)</mo>
</mrow>
<mrow>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mo>&CenterDot;</mo>
<mi>exp</mi>
<mo>&lsqb;</mo>
<mo>-</mo>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<mi>v</mi>
<mi>c</mi>
</mfrac>
<mo>)</mo>
</mrow>
<mi>k</mi>
</msup>
<mo>&rsqb;</mo>
</mrow>
其中:v为风速,k和c为威布尔分布的两个参数;
1-13)计算风力发电***的风机出力Pw,计算公式如下:
<mrow>
<msub>
<mi>P</mi>
<mi>w</mi>
</msub>
<mo>=</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mrow>
<mi>v</mi>
<mo>&le;</mo>
<msub>
<mi>v</mi>
<mrow>
<mi>c</mi>
<mi>i</mi>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>k</mi>
<mn>1</mn>
</msub>
<mo>&CenterDot;</mo>
<mi>v</mi>
<mo>+</mo>
<msub>
<mi>k</mi>
<mn>2</mn>
</msub>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>v</mi>
<mrow>
<mi>c</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>&le;</mo>
<mi>v</mi>
<mo>&le;</mo>
<msub>
<mi>v</mi>
<mi>r</mi>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>P</mi>
<mi>r</mi>
</msub>
</mtd>
<mtd>
<mrow>
<msub>
<mi>v</mi>
<mi>r</mi>
</msub>
<mo>&le;</mo>
<mi>v</mi>
<mo>&le;</mo>
<msub>
<mi>v</mi>
<mrow>
<mi>c</mi>
<mi>o</mi>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mrow>
<msub>
<mi>v</mi>
<mrow>
<mi>c</mi>
<mi>o</mi>
</mrow>
</msub>
<mo>&le;</mo>
<mi>v</mi>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
其中:k2=-k1·vci,Pr为风机额定功率,vci为切入风速,vr为额定风速,vco为切出风速;
1-14)计算风力发电***的风机出力的概率密度f(Pw),计算公式如下:
<mrow>
<mi>f</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>P</mi>
<mi>w</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mi>exp</mi>
<mo>&lsqb;</mo>
<mo>-</mo>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<msub>
<mi>P</mi>
<mi>w</mi>
</msub>
<mo>-</mo>
<msub>
<mi>k</mi>
<mn>2</mn>
</msub>
</mrow>
<mrow>
<msub>
<mi>k</mi>
<mn>1</mn>
</msub>
<mi>c</mi>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mi>k</mi>
</msup>
<mo>&rsqb;</mo>
<mo>*</mo>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<msub>
<mi>P</mi>
<mi>w</mi>
</msub>
<mo>-</mo>
<msub>
<mi>k</mi>
<mn>2</mn>
</msub>
</mrow>
<mrow>
<msub>
<mi>k</mi>
<mn>1</mn>
</msub>
<mi>c</mi>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mrow>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mo>*</mo>
<mfrac>
<mi>k</mi>
<mrow>
<msub>
<mi>k</mi>
<mn>1</mn>
</msub>
<mi>c</mi>
</mrow>
</mfrac>
<mo>.</mo>
</mrow>
4.根据权利要求1所述的一种基于多代理***的智能家居优化方法,其特征在于,所述步骤1)太阳能光伏发电***出力模型的构建方法为:
1-21)根据历史太阳能光照强度数据,计算太阳能光伏发电***的日照强度;
1-22)计算太阳能光伏发电***日照强度的概率密度f(r),计算公式如下:
<mrow>
<mi>f</mi>
<mrow>
<mo>(</mo>
<mi>r</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<mi>&Gamma;</mi>
<mrow>
<mo>(</mo>
<mi>&alpha;</mi>
<mo>+</mo>
<mi>&beta;</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>&Gamma;</mi>
<mrow>
<mo>(</mo>
<mi>&alpha;</mi>
<mo>)</mo>
</mrow>
<mo>&CenterDot;</mo>
<mi>&Gamma;</mi>
<mrow>
<mo>(</mo>
<mi>&beta;</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>&CenterDot;</mo>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<mi>r</mi>
<msub>
<mi>r</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mrow>
<mi>&alpha;</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mo>&CenterDot;</mo>
<msup>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<mfrac>
<mi>r</mi>
<msub>
<mi>r</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mrow>
<mi>&beta;</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
</mrow>
其中:r和rmax分别是计算时间段内的实际光强和最大光强,α和β为Beta分布的形状参数,Γ为Gamma函数;
1-23)计算太阳能光伏发电***的输出功率PM,计算公式如下:
PM=r·A·η
其中:A和η分别为太阳能光伏发电***的面积和光电转换效率;
1-24)计算太阳能光伏发电***输出功率的概率密度f(PM),计算公式如下:
<mrow>
<mi>f</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>P</mi>
<mi>M</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msup>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<mfrac>
<msub>
<mi>P</mi>
<mi>M</mi>
</msub>
<msub>
<mi>R</mi>
<mi>M</mi>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mrow>
<mi>&beta;</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mo>*</mo>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<msub>
<mi>P</mi>
<mi>M</mi>
</msub>
<msub>
<mi>R</mi>
<mi>M</mi>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mrow>
<mi>&alpha;</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mo>*</mo>
<mfrac>
<mrow>
<mi>&Gamma;</mi>
<mrow>
<mo>(</mo>
<mi>&alpha;</mi>
<mo>+</mo>
<mi>&beta;</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>&Gamma;</mi>
<mrow>
<mo>(</mo>
<mi>&alpha;</mi>
<mo>)</mo>
</mrow>
<mo>&CenterDot;</mo>
<mi>&Gamma;</mi>
<mrow>
<mo>(</mo>
<mi>&beta;</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
</mrow>
其中:RM表示太阳能光伏发电***输出的最大输出功率。
5.根据权利要求1所述的一种基于多代理***的智能家居优化方法,其特征在于,所述步骤1)电动汽车充放电***出力模型的构建过程为:
1-31)根据历史电动汽车充电负荷数据,计算电动汽车充电功率的均值和方差;
1-32)电动汽车类型kev,电池的剩余电量Cbat服从二维正太分布:
<mrow>
<mi>f</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>k</mi>
<mrow>
<mi>e</mi>
<mi>v</mi>
</mrow>
</msub>
<mo>,</mo>
<msub>
<mi>C</mi>
<mrow>
<mi>b</mi>
<mi>a</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mrow>
<mn>2</mn>
<mi>&pi;</mi>
<msqrt>
<mrow>
<mi>det</mi>
<mi> </mi>
<mi>C</mi>
</mrow>
</msqrt>
</mrow>
</mfrac>
<mi>exp</mi>
<mo>{</mo>
<mo>-</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<msup>
<mrow>
<mo>(</mo>
<mi>X</mi>
<mo>-</mo>
<mi>&mu;</mi>
<mo>)</mo>
</mrow>
<mi>T</mi>
</msup>
<msup>
<mi>C</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mrow>
<mo>(</mo>
<mi>X</mi>
<mo>-</mo>
<mi>&mu;</mi>
<mo>)</mo>
</mrow>
<mo>}</mo>
</mrow>
其中,X表示变量,μ表示电动汽车充电功率的均值,C电池的剩余电量参数;
电动汽车充电功率PCI服从正太分布:
<mrow>
<mi>f</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>C</mi>
<mi>I</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mrow>
<msqrt>
<mrow>
<mn>2</mn>
<mi>&pi;</mi>
<mi>d</mi>
</mrow>
</msqrt>
<msub>
<mi>&sigma;</mi>
<mrow>
<mi>p</mi>
<mi>c</mi>
<mi>i</mi>
</mrow>
</msub>
</mrow>
</mfrac>
<mi>exp</mi>
<mo>{</mo>
<mo>-</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>C</mi>
<mi>I</mi>
</mrow>
</msub>
<mo>-</mo>
<mi>&mu;</mi>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<msup>
<msub>
<mi>&sigma;</mi>
<mrow>
<mi>p</mi>
<mi>c</mi>
<mi>i</mi>
</mrow>
</msub>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mo>}</mo>
</mrow>
其中,d表示充电效率,σpci表示充电功率的方差;
电动汽车充电时刻T服从正太分布:
<mrow>
<mi>f</mi>
<mrow>
<mo>(</mo>
<mi>T</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mfrac>
<mn>1</mn>
<mrow>
<msqrt>
<mrow>
<mn>2</mn>
<mi>&pi;</mi>
<mi>d</mi>
</mrow>
</msqrt>
<msub>
<mi>&sigma;</mi>
<mrow>
<mi>p</mi>
<mi>c</mi>
<mi>i</mi>
</mrow>
</msub>
</mrow>
</mfrac>
<mi>exp</mi>
<mo>{</mo>
<mo>-</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<msup>
<mrow>
<mo>(</mo>
<mi>T</mi>
<mo>-</mo>
<msub>
<mi>&mu;</mi>
<mi>t</mi>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<msup>
<msub>
<mi>&sigma;</mi>
<mrow>
<mi>p</mi>
<mi>c</mi>
<mi>i</mi>
</mrow>
</msub>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mo>}</mo>
<mn>..................</mn>
<mrow>
<mo>(</mo>
<msub>
<mi>&mu;</mi>
<mi>t</mi>
</msub>
<mo>-</mo>
<mn>12</mn>
<mo>)</mo>
</mrow>
<mo><</mo>
<mi>T</mi>
<mo><</mo>
<mn>24</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mfrac>
<mn>1</mn>
<mrow>
<msqrt>
<mrow>
<mn>2</mn>
<mi>&pi;</mi>
<mi>d</mi>
</mrow>
</msqrt>
<msub>
<mi>&sigma;</mi>
<mrow>
<mi>p</mi>
<mi>c</mi>
<mi>i</mi>
</mrow>
</msub>
</mrow>
</mfrac>
<mi>exp</mi>
<mo>{</mo>
<mo>-</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<msup>
<mrow>
<mo>(</mo>
<mi>T</mi>
<mo>+</mo>
<mn>24</mn>
<mo>-</mo>
<msub>
<mi>&mu;</mi>
<mi>t</mi>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<msup>
<msub>
<mi>&sigma;</mi>
<mrow>
<mi>p</mi>
<mi>c</mi>
<mi>i</mi>
</mrow>
</msub>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mo>}</mo>
<mn>............0</mn>
<mo><</mo>
<mi>T</mi>
<mo><</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>&mu;</mi>
<mi>t</mi>
</msub>
<mo>-</mo>
<mn>12</mn>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
其中,μt为电动汽车充电时刻均值;
电动汽车日行驶路程M服从正太分布:
<mrow>
<mi>f</mi>
<mrow>
<mo>(</mo>
<mi>M</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mrow>
<msqrt>
<mrow>
<mn>2</mn>
<msub>
<mi>&pi;d</mi>
<mi>d</mi>
</msub>
</mrow>
</msqrt>
<msub>
<mi>&sigma;</mi>
<mi>M</mi>
</msub>
</mrow>
</mfrac>
<mi>exp</mi>
<mo>{</mo>
<mo>-</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<msup>
<mrow>
<mo>(</mo>
<mi>I</mi>
<mi>n</mi>
<mi>M</mi>
<mo>-</mo>
<mi>&mu;</mi>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<msup>
<msub>
<mi>&sigma;</mi>
<mi>M</mi>
</msub>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mo>}</mo>
</mrow>
其中,dd表示日行驶里程数损耗;σM表示日行驶路程数方差。
6.根据权利要求1所述的一种基于多代理***的智能家居优化方法,其特征在于,所述步骤3)中的操作包括:Ⅰ)从电网购电;Ⅱ)对储能电池进行充电或者放电;Ⅲ)向电网售电。
7.根据权利要求1所述的一种基于多代理***的智能家居优化方法,其特征在于,所述步骤3)中采用遗传算法进行整个***的经济性优化的具体过程如下:
(a)对参数Store(h)、Sell(h)、Eh进行采用遗传算法的二进制编码;
(b)初始化种群迭代次数t1=0,个体迭代次数t2=0,计数参数n=0,种群数组Node,设置种群大小为Nall,种群和个体最大迭代次数T1,T2;
(c)生成初始种群;
(d)计算种群中每个个体的适应度F和聚集距离P;
(e)采用选择操作从种群中选取n1个个体,并对聚集距离从大到小进行排序,选取前n2个个体,其中,n1和n2满足n1+n2=Nall;
(f)对种群进行遗传和变异操作;
(g)重复步骤(d)~(f)直至达到最大迭代次数;
(h)输出最优解并退出。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410267513.9A CN104112168B (zh) | 2014-06-16 | 2014-06-16 | 一种基于多代理***的智能家居优化方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410267513.9A CN104112168B (zh) | 2014-06-16 | 2014-06-16 | 一种基于多代理***的智能家居优化方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104112168A CN104112168A (zh) | 2014-10-22 |
CN104112168B true CN104112168B (zh) | 2017-08-29 |
Family
ID=51708952
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410267513.9A Active CN104112168B (zh) | 2014-06-16 | 2014-06-16 | 一种基于多代理***的智能家居优化方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104112168B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104376364B (zh) * | 2014-11-21 | 2017-10-24 | 国家电网公司 | 基于遗传算法的智能家居负荷管理优化方法 |
CN110246041B (zh) * | 2019-06-21 | 2020-08-25 | 贵州电网有限责任公司 | 一种基于区块链的p2p能源交易平台的交易方法 |
CN111291920B (zh) * | 2020-01-10 | 2022-04-08 | 国网重庆市电力公司电力科学研究院 | 一种考虑多因素多场景的电动汽车充电站推荐方法 |
CN115268269B (zh) * | 2022-07-29 | 2023-06-02 | 无锡市低碳研究院有限公司 | 一种基于新能源低碳的家居耗能优化***及方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102130507A (zh) * | 2011-03-26 | 2011-07-20 | 河南省电力公司 | 接入智能电网与新能源的智能家居***及能效管理方法 |
CN102361328A (zh) * | 2011-10-25 | 2012-02-22 | 中国科学技术大学 | 一种利用风能、光能互补并与市电综合利用的分布式微网*** |
CN102946100A (zh) * | 2012-10-25 | 2013-02-27 | 国网能源研究院 | 一种电力用户响应***的智能体 |
CN103617566A (zh) * | 2013-12-05 | 2014-03-05 | 中国能源建设集团广东省电力设计研究院 | 一种基于实时电价的智能用电*** |
CN103679299A (zh) * | 2013-12-30 | 2014-03-26 | 华北电力大学(保定) | 兼顾车主满意度的电动汽车最优峰谷分时电价定价方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2009046132A1 (en) * | 2007-10-01 | 2009-04-09 | Gridpoint, Inc. | Modular electrical grid interface device |
US20110047203A1 (en) * | 2009-08-20 | 2011-02-24 | Parmet Marc J | Content Aware Augmented Transform and Transfer Organizer (CAATTO) |
-
2014
- 2014-06-16 CN CN201410267513.9A patent/CN104112168B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102130507A (zh) * | 2011-03-26 | 2011-07-20 | 河南省电力公司 | 接入智能电网与新能源的智能家居***及能效管理方法 |
CN102361328A (zh) * | 2011-10-25 | 2012-02-22 | 中国科学技术大学 | 一种利用风能、光能互补并与市电综合利用的分布式微网*** |
CN102946100A (zh) * | 2012-10-25 | 2013-02-27 | 国网能源研究院 | 一种电力用户响应***的智能体 |
CN103617566A (zh) * | 2013-12-05 | 2014-03-05 | 中国能源建设集团广东省电力设计研究院 | 一种基于实时电价的智能用电*** |
CN103679299A (zh) * | 2013-12-30 | 2014-03-26 | 华北电力大学(保定) | 兼顾车主满意度的电动汽车最优峰谷分时电价定价方法 |
Non-Patent Citations (1)
Title |
---|
基于概率潮流的间歇性电源优化配置研究;邓威;《中国博士学位论文全文数据库工程科技Ⅱ辑》;20130415(第4期);论文第3.4.1.2节,第5.2.2节 * |
Also Published As
Publication number | Publication date |
---|---|
CN104112168A (zh) | 2014-10-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Zhang et al. | Multi-objective load dispatch for microgrid with electric vehicles using modified gravitational search and particle swarm optimization algorithm | |
Xu et al. | Data-driven configuration optimization of an off-grid wind/PV/hydrogen system based on modified NSGA-II and CRITIC-TOPSIS | |
Azuatalam et al. | Energy management of small-scale PV-battery systems: A systematic review considering practical implementation, computational requirements, quality of input data and battery degradation | |
Jannati et al. | Optimal performance of electric vehicles parking lot considering environmental issue | |
Mei et al. | Game approaches for hybrid power system planning | |
Liu et al. | An optimization strategy of controlled electric vehicle charging considering demand side response and regional wind and photovoltaic | |
CN105071389B (zh) | 计及源网荷互动的交直流混合微电网优化运行方法及装置 | |
Petrucci et al. | Modelling of a multi-stage energy management control routine for energy demand forecasting, flexibility, and optimization of smart communities using a Recurrent Neural Network | |
Barelli et al. | Residential micro-grid load management through artificial neural networks | |
Yang et al. | Resilient residential energy management with vehicle-to-home and photovoltaic uncertainty | |
CN105162113B (zh) | 一种基于灵敏度分析的微电网与配电网互动成本计算方法 | |
Zhang et al. | Grid-connected photovoltaic battery systems: A comprehensive review and perspectives | |
CN104112168B (zh) | 一种基于多代理***的智能家居优化方法 | |
CN112101632A (zh) | 基于分布式估计算法楼宇负荷优化调度方法 | |
CN112836849A (zh) | 一种考虑风电不确定性的虚拟电厂调度方法 | |
Banu et al. | Artificial intelligence with attention based BiLSTM for energy storage system in hybrid renewable energy sources | |
Ceccon et al. | Intelligent electric power management system for economic maximization in a residential prosumer unit | |
Barone et al. | The role of energy communities in electricity grid balancing: A flexible tool for smart grid power distribution optimization | |
CN108921368B (zh) | 基于虚拟发电厂的均衡合作博弈控制器 | |
CN113988471A (zh) | 一种微电网运行多目标优化方法 | |
Zhang et al. | A review on capacity sizing and operation strategy of grid-connected photovoltaic battery systems | |
CN112510690B (zh) | 考虑风火储联合和需求响应奖惩的优化调度方法及*** | |
CN117543581A (zh) | 考虑电动汽车需求响应的虚拟电厂优化调度方法及其应用 | |
Kang et al. | Multi-objective sizing and real-time scheduling of battery energy storage in energy-sharing community based on reinforcement learning | |
Wang et al. | Learning-based energy management policy with battery depth-of-discharge considerations |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |