CN105260798A - 一种多维负荷特性分析的大数据挖掘器 - Google Patents
一种多维负荷特性分析的大数据挖掘器 Download PDFInfo
- Publication number
- CN105260798A CN105260798A CN201510685672.5A CN201510685672A CN105260798A CN 105260798 A CN105260798 A CN 105260798A CN 201510685672 A CN201510685672 A CN 201510685672A CN 105260798 A CN105260798 A CN 105260798A
- Authority
- CN
- China
- Prior art keywords
- formula
- day
- load
- power consumer
- demand response
- 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.)
- Pending
Links
- 238000004458 analytical method Methods 0.000 title claims abstract description 53
- 230000004044 response Effects 0.000 claims abstract description 125
- 230000005856 abnormality Effects 0.000 claims abstract description 17
- 238000003745 diagnosis Methods 0.000 claims abstract description 17
- 230000007246 mechanism Effects 0.000 claims abstract description 17
- 238000011156 evaluation Methods 0.000 claims abstract description 13
- 230000002354 daily effect Effects 0.000 claims description 50
- 239000011159 matrix material Substances 0.000 claims description 39
- 230000005611 electricity Effects 0.000 claims description 28
- 238000000034 method Methods 0.000 claims description 15
- 238000012549 training Methods 0.000 claims description 15
- 239000010755 BS 2869 Class G Substances 0.000 claims description 12
- 230000007115 recruitment Effects 0.000 claims description 12
- 230000009466 transformation Effects 0.000 claims description 12
- 238000005259 measurement Methods 0.000 claims description 6
- 230000009467 reduction Effects 0.000 claims description 6
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 claims description 3
- 230000003203 everyday effect Effects 0.000 claims description 3
- 238000005516 engineering process Methods 0.000 abstract description 7
- 230000000694 effects Effects 0.000 abstract 1
- 230000002452 interceptive effect Effects 0.000 abstract 1
- 238000005065 mining Methods 0.000 abstract 1
- 238000011161 development Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000009412 basement excavation Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000012512 characterization method Methods 0.000 description 1
- 238000010224 classification analysis Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 239000013589 supplement Substances 0.000 description 1
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明涉及一种多维负荷特性分析的大数据挖掘器,包括:所述电力用户用电行为模式分析模块分别连接所述负荷特性指标分析模块、用电行为模式的影响因素和影响机理分析模块和负荷曲线异常诊断模块,所述用电行为模式的影响因素和影响机理分析模块分别连接所述电力用户需求响应基线计算模块和需求响应潜力分析模块,所述电力用户需求响应基线计算模块连接所述电力用户需求响应效果评估模块,所述负荷曲线异常诊断模块连接所述突发事件负荷预测模块,本发明提供一种多维负荷特性分析的大数据挖掘器,对大数据环境下电力用户的负荷特性进行了深入挖掘,为双向互动需求响应技术提供了科学理论支持,充分发挥了其运行效率。
Description
技术领域
本发明涉及电力系统及其自动化技术领域,具体涉及一种多维负荷特性分析的大数据挖掘器。
背景技术
负荷特性分析是智能电网研究的基础,传统的电力系统用电特性分析方法主要有两类:一是按影响因素分析,即在其余变量保持不变的前提下逐个提取因变量,定性或者定量描述其对自变量的影响程度,但对象都是区域电网,远远不及电力供应侧的研究精度和深度,无法适应推行智能需求侧管理的需要;二是按行业分类分析,即先对各行业或者细化到某类型用户的用电特性进行分析,定性或者定量得到各类型用户对区域电网用电特性的影响。但是其精确度不高,并且很难根据分析结果进行负荷预测。
总体而言现有负荷特性分析技术数据利用率低,分析结果粗糙,精度不高,且无法深入到用户层面,不利于需求响应策略制定的精细化。
大数据技术的迅速发展给也负荷分析工作带来了更多的机遇和挑战。研究大数据环境下,负荷特性分析的新思路、新方法具有重要意义。
发明内容
针对现有技术的不足,本发明提供一种多维负荷特性分析的大数据挖掘器,对大数据环境下电力用户的负荷特性进行了深入挖掘,为双向互动需求响应技术提供了科学理论支持,充分发挥了其运行效率。
本发明的目的是采用下述技术方案实现的:
一种多维负荷特性分析的大数据挖掘器,其改进之处在于,所述大数据挖掘器包括:
电力用户用电行为模式分析模块、负荷特性指标分析模块、用电行为模式的影响因素和影响机理分析模块、负荷曲线异常诊断模块、电力用户需求响应基线计算模块、电力用户需求响应效果评估模块、需求响应潜力分析模块、突发事件负荷预测模块;
所述电力用户用电行为模式分析模块分别连接所述负荷特性指标分析模块、用电行为模式的影响因素和影响机理分析模块和负荷曲线异常诊断模块,所述用电行为模式的影响因素和影响机理分析模块分别连接所述电力用户需求响应基线计算模块和需求响应潜力分析模块,所述电力用户需求响应基线计算模块连接所述电力用户需求响应效果评估模块,所述负荷曲线异常诊断模块连接所述突发事件负荷预测模块。
优选的,所述电力用户用电行为模式分析模块,用于采用聚类分析方法对电力用户的日负荷曲线进行分类,每一类对应一种用电模式,并使用类的聚类中心表征该类对应用电模式的负荷;所述负荷特性指标分析模块,用于根据所述用电模式的负荷确定所述用电模式的负荷特性指标;所述用电行为模式的影响因素和影响机理分析模块,用于根据所述用电模式的负荷提取所述用电模式的影响因素及影响因素对应的特征量;所述负荷曲线异常诊断模块,用于根据所述用电模式的负荷获取所述电力用户的日负荷曲线的奇异点及产生所述奇异点的原因;所述电力用户需求响应基线计算模块,用于预测电力用户需求响应基线值;所述电力用户需求响应效果评估模块,用于确定所述电力用户参与需求响应后实际负荷的削减量;所述需求响应潜力分析模块,用于利用贝叶斯分类法预测指定预测日用户行为模式的聚类中心Cs;所述突发事件负荷预测模块,用于根据所述奇异点的原因预测奇异点对应的社会事件下的负荷。
进一步的,所述电力用户用电行为模式分析模块采用最短距离法对日负荷数据进行系统聚类,获取初始聚类中心c,所述最短距离法的类间距离公式为:
DKL=min{dij:pi∈Gk,pj∈GL}(1)
式(1)中,pi为类Gk中元素,对应第i日用户的负荷数据序列,pj为类GL中元素,对应第j日用户的负荷数据序列,DKL为类Gk与类GL之间的距离,dij为元素pi与元素pj之间的距离,公式为:
式(2)中,pi,t为用户第i日t时刻的负荷值,pj,t为用户第j日t时刻的负荷值;
采用模糊C均值聚类法进行第二次聚类,包括:
计算元素pj关于第i类初始聚类中心ci的隶属度uij,公式为:
式(3)中,uij介于0至1之间,m为初始聚类中心个数,q为加权指数,q∈[1,∞),dij′为pj与第i类初始聚类中心ci之间的距离,公式为:
式(4)中,pj,t为用户第j日t时刻的负荷值,cit为第i类初始聚类中心ci中第t个数据点;
计算价值函数J(U,c1,c2,……,cm),公式为:
式(5)中,U为隶属矩阵,n为聚类中心中元素个数;价值函数值若小于阈值A或相对于上一次迭代中的价值函数值的变量小于阈值B,则停止迭代;
更新所述初始聚类中心c,其中,新的聚类中心Ci公式为:
上述公式中,m为初始聚类中心个数,n为聚类中心中元素个数。
进一步的,所述用电行为模式的影响因素和影响机理分析模块构建观测样本矩阵X:
式(7)中,m为观测日总数,n为待分析因素总数,xij为第i个观测日在第j个待分析影响因素上的观测值;
对所述观测样本矩阵中元素xij进行标准化变换,公式为:
式(8)中,x′ij为第i个观测日在第j个待分析影响因素上经标准化变换后的观测值,为观测日在第j个待分析影响因素上的平均值,公式为:
构建观测样本矩阵X的相关矩阵R,其中,所述矩阵R中元素rjg的公式为:
式(10)中,x′ij为第i个观测日在第j个待分析影响因素上经标准化变换后的观测值,x′ig为第i个观测日在第g个待分析影响因素上经标准化变换后的观测值;
根据所述相关矩阵R的特征方程|λI-R|=0,获取所述相关矩阵R的第j个待分析因素的特征量λj,其中,I为所述相关矩阵R的单位矩阵,0<j<n,n为待分析因素总数;
将所述相关矩阵R的待分析因素的特征量从大到小进行排序,提取前p个且满足公式(11)的所述相关矩阵R的待分析因素的特征量,公式(11)为:
式(11)中,p为提取所述相关矩阵R的待分析因素的特征量的总数;
上述公式中,m为观测日总数。
进一步的,所述负荷曲线异常诊断模块获取用户t时刻的平均负荷值公式为:
式(12)中,n为观测日总数,t为每日观测时间点数,pi,t为电力用户第i日t时刻的负荷;
获取电力用户t时刻的负荷的方差公式为:
若电力用户第i日t时刻的负荷pi,t满足则所述电力用户第i日t时刻的负荷pi,t为电力用户的日负荷曲线的奇异点,其中,ε为阈值,取1-1.5,为用户t时刻的平均负荷值;
获取电力用户的日负荷曲线的奇异点对应的社会事件,确定电力用户的日负荷曲线的奇异点对应的社会事件的奇异程度β,公式为:
式(14)中,βk为第k类社会事件的负荷奇异程度,Nk为第k类社会事件的奇异点的日负荷曲线样本总数,Δsumpk,h为第k类社会事件的日负荷曲线h中t时刻的奇异点与第k类社会事件的日负荷曲线h对应的相似日的日负荷曲线中t时刻的正常值距离之和,公式为:
式(15)中,ph,t第k类社会事件的日负荷曲线h中t时刻的奇异点负荷数据,为第k类社会事件的日负荷曲线h对应的相似日日负荷曲线中t时刻的正常值;
上述公式中,n为观测日总数。
进一步的,所述电力用户需求响应基线计算模块获取w个电力用户非需求响应日的特征量序列xi={xi(1),xi(2),xi(3),……,xi(p)},i=1,2,……,w,xi(p)为电力用户非需求响应日i的第p个特征量;
获取电力用户需求响应日的特征量序列y={y(1),y(2),y(3),……,y(p)},y(p)为电力用户需求响应日第p个特征量;
确定电力用户非需求响应日的特征量序列与电力用户需求响应日的特征量序列的关联系数,公式为:
式(16)中,ξi(h)为电力用户需求响应日第h个特征量与电力用户非需求响应日i的第h个特征量的关联系数,h=1,2,……,p;
确定所述电力用户需求响应日的特征量序列与所述电力用户非需求响应日的特征量序列的关联度,公式为:
式(17)中,γi为电力用户需求响应日的特征量序列与电力用户非需求响应日的特征量序列xi的关联度,即电力用户需求响应日的特征量序列与电力用户非需求响应日的特征量序列xi的相似度,选择γi>0.9对应的电力用户非需求响应日作为相似日;
根据相似日负荷预测需求响应日的基线负荷Pbasic(t),公式为:
式(18)中,Psimi,k(t)为用户第k个相似日的96点实际负荷,t=1,2,……96,N为提取的相似日总数,k=1,2,……N,ωk为第k个相似日的加权系数,公式为:
式(19)中,γk为第k个相似日的特征量序列与所述电力用户需求响应日的特征量序列的关联度。
进一步的,所述电力用户需求响应效果评估模块确定所述电力用户参与需求响应后实际负荷的削减量PDR(t),公式为:
PDR(t)=Pbasic(t)-Preal(t)(20)
式(20)中,Pbasic(t)为电力用户需求响应日的基线负荷,Preal(t)为电力用户实施需求响应后的96点实际负荷,t=1,2,……96。
进一步的,所述需求响应潜力分析模块根据预测日的用户特征量测量值形成预测元祖X={x1,x2,……,xl},根据观测日的用户特征量测量值形成|D|个训练元祖,训练元组和相关联的类标号形成训练集D;
计算第j类聚类中心Cj先验概率P(Cj),j=1,2,……m,若类的先验概率未知,则假定P(C1)=P(C2)=……=P(Cm),否则:
P(Cj)=|Cj,D|/|D|(21)
式(21)中,|Cj,D|为D中Cj类的训练元组数,|D|为D中训练元祖数;
假定属性值有条件地相互独立,即在属性间不存在依赖关系,计算发生Cj类情况时发生X的概率P(X|Cj),公式为:
式(22)中,xk为元组X第k个属性的值,若xk为分类属性,则P(xk|Cj)为D中第k个属性的值为xk的Cj类的元组数除以D中Cj类的元组数;若xk为连续值属性,则P(xk|Cj)为:
式(23)中,和分别为Cj类训练元组第k个属性的值的均值和标准差;
根据贝叶斯定理计算发生X情况时发生Cj的概率P(Cj|X),公式为:
式(24)中,P(X)为发生X的概率,若类的先验概率未知,则使P(X|Cj)最大的类为所述预测日用户行为模式的聚类中心Cs,否则使P(X|Cj)P(Cj)最大的类为所述预测日用户行为模式的聚类中心Cs;
计算用户预测日的需求响应潜力,包括:
获取最小负荷用电模式的用户负荷Pmin(t),公式为:
式(25)中,Cg为第g类聚类中心,Cg,t为第g类聚类中心第t个数据点;
用户预测日的需求响应力R(t)为:
R(t)=Cs(t)-Pmin(t)(26)
式(26)中,Cs(t)表示t时段预测日用户行为模式的聚类中心。
进一步的,所述突发事件负荷预测模块根据所述奇异点的原因预测奇异点对应的社会事件下的负荷,公式为:
式(27)中,βk为第k类社会事件的负荷奇异程度,为第k类社会事件的日负荷曲线h对应的相似日日负荷曲线中t时刻的正常值,pk,h'为第k类社会事件的日负荷曲线h中t时刻的预测值。
本发明的有益效果:
与现有技术相比,本发明提供一种多维负荷特性分析的大数据挖掘器,能够对电力用户的负荷特性进行全面深入的分析,解决了传统负荷特性分析数据利用率低,分析结果粗糙,精度不高,且无法深入到用户层面,不利于需求响应策略制定的精细化的缺点,为需求响应运行提供科学理论支持。
附图说明
图1是本发明一种多维负荷特性分析的大数据挖掘器的结构示意图。
具体实施方式
下面结合附图对本发明的具体实施方式作详细说明。
为使本发明实施例的目的、技术方案和优点更加清楚,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其它实施例,都属于本发明保护的范围。
本发明提供的一种多维负荷特性分析的大数据挖掘器,如图1所示,所述大数据挖掘器包括:
电力用户用电行为模式分析模块、负荷特性指标分析模块、用电行为模式的影响因素和影响机理分析模块、负荷曲线异常诊断模块、电力用户需求响应基线计算模块、电力用户需求响应效果评估模块、需求响应潜力分析模块、突发事件负荷预测模块;
其中,电力用户用电行为模式分析模块、负荷特性指标分析模块、用电行为模式的影响因素和影响机理分析模块和负荷曲线异常诊断模块为基础分析模块,电力用户需求响应基线计算模块、电力用户需求响应效果评估模块、需求响应潜力分析模块和突发事件负荷预测模块为高级应用分析模块;
所述电力用户用电行为模式分析模块分别连接所述负荷特性指标分析模块、用电行为模式的影响因素和影响机理分析模块和负荷曲线异常诊断模块,所述用电行为模式的影响因素和影响机理分析模块分别连接所述电力用户需求响应基线计算模块和需求响应潜力分析模块,所述电力用户需求响应基线计算模块连接所述电力用户需求响应效果评估模块,所述负荷曲线异常诊断模块连接所述突发事件负荷预测模块。
具体的,所述电力用户用电行为模式分析模块,用于采用聚类分析方法对电力用户的日负荷曲线进行分类,每一类对应一种用电模式,并使用类的聚类中心表征该类对应用电模式的负荷包括:
采用最短距离法对日负荷数据进行系统聚类,获取初始聚类中心c,所述最短距离法的类间距离公式为:
DKL=min{dij:pi∈Gk,pj∈GL}(1)
式(1)中,pi为类Gk中元素,对应第i日用户的负荷数据序列,pj为类GL中元素,对应第j日用户的负荷数据序列,DKL为类Gk与类GL之间的距离,dij为元素pi与元素pj之间的距离,公式为:
式(2)中,pi,t为用户第i日t时刻的负荷值,pj,t为用户第j日t时刻的负荷值;
采用模糊C均值聚类法进行第二次聚类,包括:
计算元素pj关于第i类初始聚类中心ci的隶属度uij,公式为:
式(3)中,uij介于0至1之间,m为初始聚类中心个数,q为加权指数,q∈[1,∞),dij′为pj与第i类初始聚类中心ci之间的距离,公式为:
式(4)中,pj,t为用户第j日t时刻的负荷值,cit为第i类初始聚类中心ci中第t个数据点;
计算价值函数J(U,c1,c2,……,cm),公式为:
式(5)中,U为隶属矩阵,n为聚类中心中元素个数;价值函数值若小于阈值A或相对于上一次迭代中的价值函数值的变量小于阈值B,则停止迭代;
更新所述初始聚类中心c,其中,新的聚类中心Ci公式为:
上述公式中,m为初始聚类中心个数,n为聚类中心中元素个数。
所述负荷特性指标分析模块,用于根据所述用电模式的负荷确定所述用电模式的负荷特性指标;
其中,所述用电模式的负荷特性指标包括:传统负荷特性指标和新型负荷特性指标;
传统负荷特性指标包括:
a.描述类(绝对量):日最大/小负荷、日平均负荷、日峰谷差;
b.比较类(相对量):日峰谷差率、日负荷率、日最小负荷率;
c.曲线类:日负荷曲线。
新型负荷特性指标包括:
a.负荷重要性类别:对工商业用户内部设备按安全保障负荷、主要生产性负荷、辅助生产负荷和非生产性负荷分成A、B、C、D四个等级;
b.负荷可调能力:包括技术可调性和经济可调性,技术可调性是指负荷的直接或间接控制是否可以通过智能电网技术实现,经济可调性是指在技术可调且补偿有限的前提下,考虑负荷调整所带来的收益是否可以弥补损失;
c.响应时间:终端设备对电网调度的响应时间,主要为控制信号的传输时间和智能负荷控制程序的启动时间;
d.可削减/可中断的持续时间:终端设备提供需求响应的过程中,在满足如用户舒适度等人为约束的基础上,可以持续提供需求响应资源的时间,如空调负荷参加直接负荷控制项目,在满足室温要求基础上可以持续关闭的时间;
e.负荷可调整时段:终端设备开启并正常运行,且用户有意愿参与需求响应的时段,如夏季高温时段用户空调负荷正常运行且用户签订了合同参与直接负荷控制项目;
f.爬坡速率:单位时间终端用电设备能够增加或削减的功率。
所述用电行为模式的影响因素和影响机理分析模块,用于根据所述用电模式的负荷提取所述用电模式的影响因素及影响因素对应的特征量包括:
构建观测样本矩阵X:
式(7)中,m为观测日总数,n为待分析因素总数,xij为第i个观测日在第j个待分析影响因素上的观测值;
对所述观测样本矩阵中元素xij进行标准化变换,公式为:
式(8)中,x′ij为第i个观测日在第j个待分析影响因素上经标准化变换后的观测值,xj为观测日在第j个待分析影响因素上的平均值,公式为:
构建观测样本矩阵X的相关矩阵R,其中,所述矩阵R中元素rjg的公式为:
式(10)中,x′ij为第i个观测日在第j个待分析影响因素上经标准化变换后的观测值,x′ig为第i个观测日在第g个待分析影响因素上经标准化变换后的观测值;
根据所述相关矩阵R的特征方程|λI-R|=0,获取所述相关矩阵R的第j个待分析因素的特征量λj,其中,I为所述相关矩阵R的单位矩阵,0<j<n,n为待分析因素总数;
将所述相关矩阵R的待分析因素的特征量从大到小进行排序,提取前p个且满足公式(11)的所述相关矩阵R的待分析因素的特征量,公式(11)为:
式(11)中,p为提取所述相关矩阵R的待分析因素的特征量的总数;
上述公式中,m为观测日总数。
所述负荷曲线异常诊断模块,用于根据所述用电模式的负荷获取所述电力用户的日负荷曲线的奇异点及产生所述奇异点的原因包括:
获取用户t时刻的平均负荷值公式为:
式(12)中,n为观测日总数,t为每日观测时间点数,pi,t为电力用户第i日t时刻的负荷;
获取电力用户t时刻的负荷的方差公式为:
若电力用户第i日t时刻的负荷pi,t满足则所述电力用户第i日t时刻的负荷pi,t为电力用户的日负荷曲线的奇异点,其中,ε为阈值,取1-1.5,为用户t时刻的平均负荷值;
获取电力用户的日负荷曲线的奇异点对应的社会事件,确定电力用户的日负荷曲线的奇异点对应的社会事件的奇异程度β,公式为:
式(14)中,βk为第k类社会事件的负荷奇异程度,Nk为第k类社会事件的奇异点的日负荷曲线样本总数,Δsumpk,h为第k类社会事件的日负荷曲线h中t时刻的奇异点与第k类社会事件的日负荷曲线h对应的相似日的日负荷曲线中t时刻的正常值距离之和,公式为:
式(15)中,ph,t第k类社会事件的日负荷曲线h中t时刻的奇异点负荷数据,为第k类社会事件的日负荷曲线h对应的相似日日负荷曲线中t时刻的正常值;
上述公式中,n为观测日总数。
所述电力用户需求响应基线计算模块,用于预测电力用户需求响应基线值包括:
获取w个电力用户非需求响应日的特征量序列xi={xi(1),xi(2),xi(3),……,xi(p)},i=1,2,……,w,xi(p)为电力用户非需求响应日i的第p个特征量;
获取电力用户需求响应日的特征量序列y={y(1),y(2),y(3),……,y(p)},y(p)为电力用户需求响应日第p个特征量;
确定电力用户非需求响应日的特征量序列与电力用户需求响应日的特征量序列的关联系数,公式为:
式(16)中,ξi(h)为电力用户需求响应日第h个特征量与电力用户非需求响应日i的第h个特征量的关联系数,h=1,2,……,p;
确定所述电力用户需求响应日的特征量序列与所述电力用户非需求响应日的特征量序列的关联度,公式为:
式(17)中,γi为电力用户需求响应日的特征量序列与电力用户非需求响应日的特征量序列xi的关联度,即电力用户需求响应日的特征量序列与电力用户非需求响应日的特征量序列xi的相似度,选择γi>0.9对应的电力用户非需求响应日作为相似日;
根据相似日负荷预测需求响应日的基线负荷Pbasic(t),公式为:
式(18)中,Psimi,k(t)为用户第k个相似日的96点实际负荷,t=1,2,……96,N为提取的相似日总数,k=1,2,……N,ωk为第k个相似日的加权系数,公式为:
式(19)中,γk为第k个相似日的特征量序列与所述电力用户需求响应日的特征量序列的关联度。
所述电力用户需求响应效果评估模块,用于确定所述电力用户参与需求响应后实际负荷的削减量包括:
确定所述电力用户参与需求响应后实际负荷的削减量PDR(t),公式为:
PDR(t)=Pbasic(t)-Preal(t)(20)
式(20)中,Pbasic(t)为电力用户需求响应日的基线负荷,Preal(t)为电力用户实施需求响应后的96点实际负荷,t=1,2,……96。
所述需求响应潜力分析模块,用于利用贝叶斯分类法预测指定预测日用户行为模式的聚类中心Cs包括:
根据预测日的用户特征量测量值形成预测元祖X={x1,x2,……,xl},根据观测日的用户特征量测量值形成|D|个训练元祖,训练元组和相关联的类标号形成训练集D;
计算第j类聚类中心Cj先验概率P(Cj),j=1,2,……m,若类的先验概率未知,则假定P(C1)=P(C2)=……=P(Cm),否则:
P(Cj)=|Cj,D|/|D|(21)
式(21)中,|Cj,D|为D中Cj类的训练元组数,|D|为D中训练元祖数;
假定属性值有条件地相互独立,即在属性间不存在依赖关系,计算发生Cj类情况时发生X的概率P(X|Cj),公式为:
式(22)中,xk为元组X第k个属性的值,若xk为分类属性,则P(xk|Cj)为D中第k个属性的值为xk的Cj类的元组数除以D中Cj类的元组数;若xk为连续值属性,则P(xk|Cj)为:
式(23)中,和分别为Cj类训练元组第k个属性的值的均值和标准差;
根据贝叶斯定理计算发生X情况时发生Cj的概率P(Cj|X),公式为:
式(24)中,P(X)为发生X的概率,若类的先验概率未知,则使P(X|Cj)最大的类为所述预测日用户行为模式的聚类中心Cs,否则使P(X|Cj)P(Cj)最大的类为所述预测日用户行为模式的聚类中心Cs;
计算用户预测日的需求响应潜力,包括:
获取最小负荷用电模式的用户负荷Pmin(t),公式为:
式(25)中,Cg为第g类聚类中心,Cg,t为第g类聚类中心第t个数据点;
用户预测日的需求响应力R(t)为:
R(t)=Cs(t)-Pmin(t)(26)
式(26)中,Cs(t)表示t时段预测日用户行为模式的聚类中心。
所述突发事件负荷预测模块,用于根据所述奇异点的原因预测奇异点对应的社会事件下的负荷包括:
根据所述奇异点的原因预测奇异点对应的社会事件下的负荷,公式为:
式(27)中,βk为第k类社会事件的负荷奇异程度,为第k类社会事件的日负荷曲线h对应的相似日日负荷曲线中t时刻的正常值,pk,h'为第k类社会事件的日负荷曲线h中t时刻的预测值。
最后应当说明的是:以上实施例仅用以说明本发明的技术方案而非对其限制,尽管参照上述实施例对本发明进行了详细的说明,所属领域的普通技术人员应当理解:依然可以对本发明的具体实施方式进行修改或者等同替换,而未脱离本发明精神和范围的任何修改或者等同替换,其均应涵盖在本发明的权利要求保护范围之内。
Claims (9)
1.一种多维负荷特性分析的大数据挖掘器,其特征在于,所述大数据挖掘器包括:
电力用户用电行为模式分析模块、负荷特性指标分析模块、用电行为模式的影响因素和影响机理分析模块、负荷曲线异常诊断模块、电力用户需求响应基线计算模块、电力用户需求响应效果评估模块、需求响应潜力分析模块、突发事件负荷预测模块;
所述电力用户用电行为模式分析模块分别连接所述负荷特性指标分析模块、用电行为模式的影响因素和影响机理分析模块和负荷曲线异常诊断模块,所述用电行为模式的影响因素和影响机理分析模块分别连接所述电力用户需求响应基线计算模块和需求响应潜力分析模块,所述电力用户需求响应基线计算模块连接所述电力用户需求响应效果评估模块,所述负荷曲线异常诊断模块连接所述突发事件负荷预测模块。
2.如权利要求1所述的大数据挖掘器,其特征在于,所述电力用户用电行为模式分析模块,用于采用聚类分析方法对电力用户的日负荷曲线进行分类,每一类对应一种用电模式,并使用类的聚类中心表征该类对应用电模式的负荷;所述负荷特性指标分析模块,用于根据所述用电模式的负荷确定所述用电模式的负荷特性指标;所述用电行为模式的影响因素和影响机理分析模块,用于根据所述用电模式的负荷提取所述用电模式的影响因素及影响因素对应的特征量;所述负荷曲线异常诊断模块,用于根据所述用电模式的负荷获取所述电力用户的日负荷曲线的奇异点及产生所述奇异点的原因;所述电力用户需求响应基线计算模块,用于预测电力用户需求响应基线值;所述电力用户需求响应效果评估模块,用于确定所述电力用户参与需求响应后实际负荷的削减量;所述需求响应潜力分析模块,用于利用贝叶斯分类法预测指定预测日用户行为模式的聚类中心Cs;所述突发事件负荷预测模块,用于根据所述奇异点的原因预测奇异点对应的社会事件下的负荷。
3.如权利要求2所述的大数据挖掘器,其特征在于,所述电力用户用电行为模式分析模块采用最短距离法对日负荷数据进行系统聚类,获取初始聚类中心c,所述最短距离法的类间距离公式为:
DKL=min{dij:pi∈Gk,pj∈GL}(1)
式(1)中,pi为类Gk中元素,对应第i日用户的负荷数据序列,pj为类GL中元素,对应第j日用户的负荷数据序列,DKL为类Gk与类GL之间的距离,dij为元素pi与元素pj之间的距离,公式为:
式(2)中,pi,t为用户第i日t时刻的负荷值,pj,t为用户第j日t时刻的负荷值;
采用模糊C均值聚类法进行第二次聚类,包括:
计算元素pj关于第i类初始聚类中心ci的隶属度uij,公式为:
式(3)中,uij介于0至1之间,m为初始聚类中心个数,q为加权指数,q∈[1,∞),dij′为pj与第i类初始聚类中心ci之间的距离,公式为:
式(4)中,pj,t为用户第j日t时刻的负荷值,cit为第i类初始聚类中心ci中第t个数据点;
计算价值函数J(U,c1,c2,……,cm),公式为:
式(5)中,U为隶属矩阵,n为聚类中心中元素个数;价值函数值若小于阈值A或相对于上一次迭代中的价值函数值的变量小于阈值B,则停止迭代;
更新所述初始聚类中心c,其中,新的聚类中心Ci公式为:
上述公式中,m为初始聚类中心个数,n为聚类中心中元素个数。
4.如权利要求2所述的大数据挖掘器,其特征在于,所述用电行为模式的影响因素和影响机理分析模块构建观测样本矩阵X:
式(7)中,m为观测日总数,n为待分析因素总数,xij为第i个观测日在第j个待分析影响因素上的观测值;
对所述观测样本矩阵中元素xij进行标准化变换,公式为:
式(8)中,x′ij为第i个观测日在第j个待分析影响因素上经标准化变换后的观测值,为观测日在第j个待分析影响因素上的平均值,公式为:
构建观测样本矩阵X的相关矩阵R,其中,所述矩阵R中元素rjg的公式为:
式(10)中,x′ij为第i个观测日在第j个待分析影响因素上经标准化变换后的观测值,x′ig为第i个观测日在第g个待分析影响因素上经标准化变换后的观测值;
根据所述相关矩阵R的特征方程|λI-R|=0,获取所述相关矩阵R的第j个待分析因素的特征量λj,其中,I为所述相关矩阵R的单位矩阵,0<j<n,n为待分析因素总数;
将所述相关矩阵R的待分析因素的特征量从大到小进行排序,提取前p个且满足公式(11)的所述相关矩阵R的待分析因素的特征量,公式(11)为:
式(11)中,p为提取所述相关矩阵R的待分析因素的特征量的总数;
上述公式中,m为观测日总数。
5.如权利要求2所述的大数据挖掘器,其特征在于,所述负荷曲线异常诊断模块获取用户t时刻的平均负荷值公式为:
式(12)中,n为观测日总数,t为每日观测时间点数,pi,t为电力用户第i日t时刻的负荷;
获取电力用户t时刻的负荷的方差公式为:
若电力用户第i日t时刻的负荷pi,t满足则所述电力用户第i日t时刻的负荷pi,t为电力用户的日负荷曲线的奇异点,其中,ε为阈值,取1-1.5,为用户t时刻的平均负荷值;
获取电力用户的日负荷曲线的奇异点对应的社会事件,确定电力用户的日负荷曲线的奇异点对应的社会事件的奇异程度β,公式为:
式(14)中,βk为第k类社会事件的负荷奇异程度,Nk为第k类社会事件的奇异点的日负荷曲线样本总数,Δsumpk,h为第k类社会事件的日负荷曲线h中t时刻的奇异点与第k类社会事件的日负荷曲线h对应的相似日的日负荷曲线中t时刻的正常值距离之和,公式为:
式(15)中,ph,t第k类社会事件的日负荷曲线h中t时刻的奇异点负荷数据,为第k类社会事件的日负荷曲线h对应的相似日日负荷曲线中t时刻的正常值;
上述公式中,n为观测日总数。
6.如权利要求2所述的大数据挖掘器,其特征在于,所述电力用户需求响应基线计算模块获取w个电力用户非需求响应日的特征量序列xi={xi(1),xi(2),xi(3),……,xi(p)},i=1,2,……,w,xi(p)为电力用户非需求响应日i的第p个特征量;
获取电力用户需求响应日的特征量序列y={y(1),y(2),y(3),……,y(p)},y(p)为电力用户需求响应日第p个特征量;
确定电力用户非需求响应日的特征量序列与电力用户需求响应日的特征量序列的关联系数,公式为:
式(16)中,ξi(h)为电力用户需求响应日第h个特征量与电力用户非需求响应日i的第h个特征量的关联系数,h=1,2,……,p;
确定所述电力用户需求响应日的特征量序列与所述电力用户非需求响应日的特征量序列的关联度,公式为:
式(17)中,γi为电力用户需求响应日的特征量序列与电力用户非需求响应日的特征量序列xi的关联度,即电力用户需求响应日的特征量序列与电力用户非需求响应日的特征量序列xi的相似度,选择γi>0.9对应的电力用户非需求响应日作为相似日;
根据相似日负荷预测需求响应日的基线负荷Pbasic(t),公式为:
式(18)中,Psimi,k(t)为用户第k个相似日的96点实际负荷,t=1,2,……96,N为提取的相似日总数,k=1,2,……N,ωk为第k个相似日的加权系数,公式为:
式(19)中,γk为第k个相似日的特征量序列与所述电力用户需求响应日的特征量序列的关联度。
7.如权利要求2所述的大数据挖掘器,其特征在于,所述电力用户需求响应效果评估模块确定所述电力用户参与需求响应后实际负荷的削减量PDR(t),公式为:
PDR(t)=Pbasic(t)-Preal(t)(20)
式(20)中,Pbasic(t)为电力用户需求响应日的基线负荷,Preal(t)为电力用户实施需求响应后的96点实际负荷,t=1,2,……96。
8.如权利要求2所述的大数据挖掘器,其特征在于,所述需求响应潜力分析模块根据预测日的用户特征量测量值形成预测元祖X={x1,x2,……,xl},根据观测日的用户特征量测量值形成|D|个训练元祖,训练元组和相关联的类标号形成训练集D;
计算第j类聚类中心Cj先验概率P(Cj),j=1,2,……m,若类的先验概率未知,则假定P(C1)=P(C2)=……=P(Cm),否则:
P(Cj)=|Cj,D|/|D|(21)
式(21)中,|Cj,D|为D中Cj类的训练元组数,|D|为D中训练元祖数;
假定属性值有条件地相互独立,即在属性间不存在依赖关系,计算发生Cj类情况时发生X的概率P(X|Cj),公式为:
式(22)中,xk为元组X第k个属性的值,若xk为分类属性,则P(xk|Cj)为D中第k个属性的值为xk的Cj类的元组数除以D中Cj类的元组数;若xk为连续值属性,则P(xk|Cj)为:
式(23)中,和分别为Cj类训练元组第k个属性的值的均值和标准差;
根据贝叶斯定理计算发生X情况时发生Cj的概率P(Cj|X),公式为:
式(24)中,P(X)为发生X的概率,若类的先验概率未知,则使P(X|Cj)最大的类为所述预测日用户行为模式的聚类中心Cs,否则使P(X|Cj)P(Cj)最大的类为所述预测日用户行为模式的聚类中心Cs;
计算用户预测日的需求响应潜力,包括:
获取最小负荷用电模式的用户负荷Pmin(t),公式为:
式(25)中,Cg为第g类聚类中心,Cg,t为第g类聚类中心第t个数据点;
用户预测日的需求响应力R(t)为:
R(t)=Cs(t)-Pmin(t)(26)
式(26)中,Cs(t)表示t时段预测日用户行为模式的聚类中心。
9.如权利要求2所述的大数据挖掘器,其特征在于,所述突发事件负荷预测模块根据所述奇异点的原因预测奇异点对应的社会事件下的负荷,公式为:
式(27)中,βk为第k类社会事件的负荷奇异程度,为第k类社会事件的日负荷曲线h对应的相似日日负荷曲线中t时刻的正常值,pk,h'为第k类社会事件的日负荷曲线h中t时刻的预测值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510685672.5A CN105260798A (zh) | 2015-10-21 | 2015-10-21 | 一种多维负荷特性分析的大数据挖掘器 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510685672.5A CN105260798A (zh) | 2015-10-21 | 2015-10-21 | 一种多维负荷特性分析的大数据挖掘器 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN105260798A true CN105260798A (zh) | 2016-01-20 |
Family
ID=55100477
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510685672.5A Pending CN105260798A (zh) | 2015-10-21 | 2015-10-21 | 一种多维负荷特性分析的大数据挖掘器 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105260798A (zh) |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106056271A (zh) * | 2016-05-17 | 2016-10-26 | 珠海许继芝电网自动化有限公司 | 一种用户群用电负荷响应的智能控制方法 |
CN107404120A (zh) * | 2017-07-25 | 2017-11-28 | 南京工程学院 | 一种无功优化在线控制中的设备动作次数挖掘方法 |
CN108764353A (zh) * | 2018-05-25 | 2018-11-06 | 中国电力科学研究院有限公司 | 基于匹配度的大规模电力系统功率曲线分类方法 |
CN110570090A (zh) * | 2019-08-12 | 2019-12-13 | 国网上海市电力公司 | 一种基于贝叶斯分类法的基线负荷模型的负荷辨识方法 |
CN111126498A (zh) * | 2019-12-25 | 2020-05-08 | 国网河北省电力有限公司 | 一种基于分类分析的客户用电行为特征分析方法 |
CN111967723A (zh) * | 2020-07-24 | 2020-11-20 | 南昌大学 | 一种基于数据挖掘的用户调峰潜力分析方法 |
CN114662761A (zh) * | 2022-03-24 | 2022-06-24 | 国网江苏省电力有限公司南通供电分公司 | 一种基于负荷预测的可参与需求响应负荷辨识方法及系统 |
CN116485071A (zh) * | 2023-04-12 | 2023-07-25 | 山东科技大学 | 一种基于概率基线负荷的电力用户需求响应潜力评估方法 |
CN118014166A (zh) * | 2024-04-09 | 2024-05-10 | 陕西德联新能源有限公司 | 基于大数据的供热系统的负荷预测方法及系统 |
-
2015
- 2015-10-21 CN CN201510685672.5A patent/CN105260798A/zh active Pending
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106056271A (zh) * | 2016-05-17 | 2016-10-26 | 珠海许继芝电网自动化有限公司 | 一种用户群用电负荷响应的智能控制方法 |
CN107404120A (zh) * | 2017-07-25 | 2017-11-28 | 南京工程学院 | 一种无功优化在线控制中的设备动作次数挖掘方法 |
CN107404120B (zh) * | 2017-07-25 | 2020-04-24 | 南京工程学院 | 一种无功优化在线控制中的设备动作次数挖掘方法 |
CN108764353B (zh) * | 2018-05-25 | 2021-11-30 | 中国电力科学研究院有限公司 | 基于匹配度的大规模电力系统功率曲线分类方法 |
CN108764353A (zh) * | 2018-05-25 | 2018-11-06 | 中国电力科学研究院有限公司 | 基于匹配度的大规模电力系统功率曲线分类方法 |
CN110570090A (zh) * | 2019-08-12 | 2019-12-13 | 国网上海市电力公司 | 一种基于贝叶斯分类法的基线负荷模型的负荷辨识方法 |
CN111126498A (zh) * | 2019-12-25 | 2020-05-08 | 国网河北省电力有限公司 | 一种基于分类分析的客户用电行为特征分析方法 |
CN111967723A (zh) * | 2020-07-24 | 2020-11-20 | 南昌大学 | 一种基于数据挖掘的用户调峰潜力分析方法 |
CN111967723B (zh) * | 2020-07-24 | 2022-05-20 | 南昌大学 | 一种基于数据挖掘的用户调峰潜力分析方法 |
CN114662761A (zh) * | 2022-03-24 | 2022-06-24 | 国网江苏省电力有限公司南通供电分公司 | 一种基于负荷预测的可参与需求响应负荷辨识方法及系统 |
CN116485071A (zh) * | 2023-04-12 | 2023-07-25 | 山东科技大学 | 一种基于概率基线负荷的电力用户需求响应潜力评估方法 |
CN118014166A (zh) * | 2024-04-09 | 2024-05-10 | 陕西德联新能源有限公司 | 基于大数据的供热系统的负荷预测方法及系统 |
CN118014166B (zh) * | 2024-04-09 | 2024-06-21 | 陕西德联新能源有限公司 | 基于大数据的供热系统的负荷预测方法及系统 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105260798A (zh) | 一种多维负荷特性分析的大数据挖掘器 | |
Liu et al. | Energy consumption prediction and diagnosis of public buildings based on support vector machine learning: A case study in China | |
CN110097297A (zh) | 一种多维度窃电态势智能感知方法、系统、设备及介质 | |
Bian et al. | Study on power consumption load forecast based on K-means clustering and FCM–BP model | |
CN106096830B (zh) | 基于广义证据理论的继电保护状态评价方法及系统 | |
He et al. | Comprehensive evaluation of regional clean energy development levels based on principal component analysis and rough set theory | |
CN104881609A (zh) | 一种复杂软件系统软件单元可信性评估方法 | |
CN103226736B (zh) | 基于聚类分析与灰靶理论的中长期电力负荷预测方法 | |
CN106909933A (zh) | 一种三阶段多视角特征融合的窃电分类预测方法 | |
CN103632203A (zh) | 一种基于综合评价的配电网供电区域划分方法 | |
CN104318482A (zh) | 一套智能配电网综合评估体系和方法 | |
CN110782153A (zh) | 一种企业园区综合能效评估体系建模方法及系统 | |
CN103093285A (zh) | 基于人工神经网络的短期负荷预测方法 | |
CN112149890A (zh) | 基于用户用能标签的综合能源负荷预测方法及系统 | |
CN110109971A (zh) | 一种低压台区用户用电负荷特性分析方法 | |
CN106410781A (zh) | 一种电力用户需求响应潜力确定方法 | |
CN104751254A (zh) | 基于非等距加权灰色模型和模糊聚类排序线损率预测方法 | |
Guan et al. | Customer load forecasting method based on the industry electricity consumption behavior portrait | |
CN104143117A (zh) | 一种电网特殊负荷与日负荷之间的相关系数提取方法 | |
CN117117857A (zh) | 一种考虑用户用电和季节特性的区域短期电力负荷预测方法 | |
Li et al. | Multi‐level refined power system operation mode analysis: A data‐driven approach | |
CN104463683A (zh) | 一种含多源电网中长期负荷预测装置及方法 | |
CN114048200A (zh) | 一种计及缺失数据补齐的用户用电行为分析方法 | |
CN115049072A (zh) | 一种在二分类中同时实现差分隐私和机器学习公平的方法 | |
Liao et al. | Building energy efficiency assessment base on predict-center criterion under diversified conditions |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20160120 |