CN107505519A - 一种分布式电源接入电网电能质量分析方法及装置 - Google Patents
一种分布式电源接入电网电能质量分析方法及装置 Download PDFInfo
- Publication number
- CN107505519A CN107505519A CN201710619336.XA CN201710619336A CN107505519A CN 107505519 A CN107505519 A CN 107505519A CN 201710619336 A CN201710619336 A CN 201710619336A CN 107505519 A CN107505519 A CN 107505519A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- mfrac
- msup
- munderover
- 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.)
- Granted
Links
- 238000004458 analytical method Methods 0.000 title claims abstract description 10
- 238000000034 method Methods 0.000 claims abstract description 48
- 238000005070 sampling Methods 0.000 claims description 28
- 238000004422 calculation algorithm Methods 0.000 claims description 10
- 238000005457 optimization Methods 0.000 claims description 7
- 230000008859 change Effects 0.000 claims description 6
- 238000012797 qualification Methods 0.000 claims description 4
- 238000007619 statistical method Methods 0.000 abstract description 2
- 230000006870 function Effects 0.000 description 92
- 238000009826 distribution Methods 0.000 description 13
- 238000012544 monitoring process Methods 0.000 description 10
- 238000010586 diagram Methods 0.000 description 9
- 238000004590 computer program Methods 0.000 description 7
- 238000012545 processing Methods 0.000 description 5
- 238000005315 distribution function Methods 0.000 description 3
- 238000003860 storage Methods 0.000 description 3
- 230000006872 improvement Effects 0.000 description 2
- 238000010248 power generation Methods 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 238000013459 approach Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000002790 cross-validation Methods 0.000 description 1
- 230000001186 cumulative effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000003012 network analysis Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R31/00—Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
- G01R31/001—Measuring interference from external sources to, or emission from, the device under test, e.g. EMC, EMI, EMP or ESD testing
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R31/00—Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
- G01R31/003—Environmental or reliability tests
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Environmental & Geological Engineering (AREA)
- Electromagnetism (AREA)
- Other Investigation Or Analysis Of Materials By Electrical Means (AREA)
- Pharmaceuticals Containing Other Organic And Inorganic Compounds (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明涉及一种分布式电源接入电网电能质量分析方法及装置,所述方法包括:获取电能质量指标的样本集合;优化基于核密度函数的电能质量指标的概率密度函数模型中窗口宽度,其中,所述基于核密度函数的电能质量指标的概率密度函数模型根据所述样本集合确定;利用所述电能质量指标的概率密度函数模型确定电能质量指标的统计特征值,本发明提供的技术方案,采用核密度估计法对电能质量指标进行分析,基于核密度估计方法求出上述电能指标的概率密度分布函数,然后再根据概率密度函数对电能指标进行统计分析。
Description
技术领域
本发明涉及配电网分析领域,具体涉及一种分布式电源接入电网电能质量分析方法及装置。
背景技术
国民经济的迅速发展对电网的电能质量提出了越来越高的要求,为了保证电能质量,国家颁布了有关标准,包括:供电电压允许偏差、电压波动与闪变、公用电网谐波、三相电压允许不平衡度、电力系统频率允许偏差等,并要求电力部门严格执行这些质量指标。
随着大规模光伏发电、风力发电等分布式电源的接入,分布式电源输出功率的随机变化给电网的电能质量带来严重影响,还引起电压偏差、电压波动与闪变、谐波畸变、三相不平衡和频率波动等问题,给电网带来了具有间歇性和随机波动性的冲击。电能质量是影响电力系统安全稳定运行的重要因素,同时也是衡量新能源电源可否接入的重要指标之一。为此,对分布式电源接入电网情况下的电能质量进行全面地分析,评估其接入对电网的电能质量所造成的影响,对最大化地接受新能源电力的有重要的指导意义。
目前,国内的电能质量监测和分析方法是对监测点的电压、电流和频率等进行监测,对采集的数据进行分析,计算出监测点的电压偏差、电压波动与闪变、谐波畸变、三相不平衡和频率波动等电能质量指标,并以日、月、年为统计周期对上述电能指标进行分析,计算出最大值、最小值、平均值和95%概率大值。但这种方法只给出了单项电能质量平均值、最大值、最小值和95%概率大值,没有给出电网电能质量指标的分布类型及统计意义下的描述,不能评估出电能质量指标超出国家标准规定限值的概率。
发明内容
本发明提供一种分布式电源接入电网电能质量分析方法及装置,其目的是采用核密度估计法对电能质量指标进行分析,基于核密度估计方法求出上述电能指标的概率密度分布函数,然后再根据概率密度函数对电能指标进行统计分析。
本发明的目的是采用下述技术方案实现的:
一种分布式电源接入电网电能质量分析方法,其改进之处在于,包括:
获取电能质量指标的样本集合;
优化基于核密度函数的电能质量指标的概率密度函数模型中窗口宽度,其中,所述基于核密度函数的电能质量指标的概率密度函数模型根据所述样本集合确定;
利用所述电能质量指标的概率密度函数模型确定电能质量指标的统计特征值。
优选的,所述电能质量指标包括:电压偏差、电压波动、电压闪变、三相电压不平衡、公用电网谐波和电力系统频率偏差。
优选的,所述基于核密度函数的电能质量指标的概率密度函数模型根据所述样本集合确定,包括:
令电能质量指标为x,则其对应的样本集合为{x1,x2,…,xn},其中,n为采样点总数;
按下式确定基于核密度函数的电能质量指标的概率密度函数模型
上式中,h为窗口宽度,xi为所述样本集合中第i个采样点对应的电能质量指标值,K(u)为核函数,其中,K(u)中的u为自定义变量。
进一步的,采用Gaussian核函数,则K(u)的公式为:
则基于Gaussian核函数的电能质量指标的概率密度函数模型的公式为:
优选的,所述优化所述电能质量指标的概率密度函数模型中窗口宽度,包括:
利用积分方差方式确定使所述电能质量指标的概率密度函数模型与实际电能质量指标的概率密度函数模型的误差为最小的窗口宽度的关于窗口宽度的积分方差表达式;
利用所述积分方差表达式,采用牛顿迭代算法获取最优窗口宽度。
进一步的,所述利用积分方差方式确定使所述电能质量指标的概率密度函数模型与实际电能质量指标的概率密度函数模型的误差为最小的窗口宽度的关于窗口宽度的积分方差表达式,包括:
按下式确定关于窗口宽度的积分方差表达式ISE*(h):
上式中,h为窗口宽度,xi为所述样本集合中第i个采样点对应的电能质量指标值,xj为所述样本集合中第i个采样点对应的电能质量指标值,n为采样点总数;
其中,使ISE*(h)最小的窗口宽度即为使所述电能质量指标的概率密度函数模型与实际电能质量指标的概率密度函数模型的误差为最小的窗口宽度。
进一步的,所述利用所述积分方差表达式,采用牛顿迭代算法获取最优窗口宽度,包括:
a.设迭代次数t=0,设定窗口宽度初始值h0及收敛条件;
b.确定所述积分方差表达式对ht的一阶导数解析式和二阶导数解析式;
c.利用所述积分方差表达式对ht的一阶导数解析式和二阶导数解析式确定迭代窗口宽度ht+1;
d.若ht+1满足所述收敛条件,则ht+1为最优窗口宽度,若ht+1不满足所述收敛条件,则令t=t+1,ht=ht+1,并返回步骤b,其中,所述收敛条件包括:|ht+1-ht|≤ε和迭代次数到达T,t∈T。
进一步的,按下式确定所述积分方差表达式对ht的一阶导数解析式:
按下式确定所述积分方差表达式对ht的二阶导数解析式:
上式中,xi为所述样本集合中第i个采样点对应的电能质量指标值,xj为所述样本集合中第i个采样点对应的电能质量指标值,n为采样点总数;
按下式利用所述积分方差表达式对ht的一阶导数解析式和二阶导数解析式确定迭代窗口宽度ht+1:
优选的,所述利用所述电能质量指标的概率密度函数模型确定电能质量指标的统计特征值,包括:
按下式确定电能质量指标的均值E(X):
按下式确定电能质量指标的方差D(X):
按下式确定电能质量指标的合格概率Pok:
按下式确定电能质量指标的超出下限值的概率PLL:
按下式确定电能质量指标的超出上限值的概率PUL:
上式中,(a,b)为电能质量指标的变化区域,Xmin为电能质量指标的下限值,Xmax为电能质量指标的上限值,为所述电能质量指标的概率密度函数模型;
将所述电能质量指标的概率密度函数模型的变化特性曲线中峰值所在的点作为电能质量指标的众数。
一种分布式电源接入电网电能质量分析装置,其改进之处在于,所述装置包括:
获取单元,用于获取电能质量指标的样本集合;
优化单元,用于优化基于核密度函数的电能质量指标的概率密度函数模型中窗口宽度,其中,所述基于核密度函数的电能质量指标的概率密度函数模型根据所述样本集合确定;
确定单元,用于利用所述电能质量指标的概率密度函数模型确定电能质量指标的统计特征值。
本发明的有益效果:
本发明提供的技术方案,采用核密度估计法对电能质量指标进行分析,基于核密度估计方法估计电能质量指标的概率密度分布函数,利用电能质量指标的概率密度分布函数,能够得出的电能质量指标统计意义上的分布规律,可以更全面的评估分布式电源接入配电网对电能质量产生的影响,让电力系统更加安全稳定的运行,进一步的,根据电能质量指标的概率密度分布函数对电能质量指标进行统计分析,不仅可以求出电能质量指标的均值和方差,还可以求出其合格和超下限及超上限的概率,以及其输出值较多的众数,从而得到一个渐进无偏的密度估计,有着良好的统计学性质,基于本发明提供的技术方案,能够大大减少电能质量指标超出国家标准所规定的状况,对最大化地接受新能源电力的有重要的指导意义。
附图说明
图1是本发明一种分布式电源接入电网电能质量分析方法的流程图;
图2是本发明实施例中基于牛顿法优化窗口宽度的流程图;
图3是本发明一种分布式电源接入电网电能质量分析装置的结构示意图。
具体实施方式
下面结合附图对本发明的具体实施方式作详细说明。
为使本发明实施例的目的、技术方案和优点更加清楚,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其它实施例,都属于本发明保护的范围。
在实际电网中,受电源输出功率变化及负荷波动等的影响,电压偏差、电压波动与闪变、谐波畸变、三相不平衡和频率波动等反映电能质量指标的电气参数是随机变化的,现有的电能质量分析方法主要是根据监测数据计算电能质量指标的平均值、最大值、最小值和95%概率大值,未能给出电网电能质量指标的分布类型及统计意义下的描述。
为了克服上述现有技术的不足,本发明提供的一种分布式电源接入电网电能质量分析方法,根据分布式电源并网点的电压偏差、电压波动与闪变、谐波畸变、三相不平衡和频率波动等监测数据,采用核密度估计方法对电压偏差、电压波动与闪变、谐波畸变、三相不平衡和频率波动等电能质量指标进行分析,求电压偏差、电压波动与闪变、谐波畸变、三相不平衡和频率波动等的概率密度函数估计,以给出电能质量指标统计意义上的分布规律。
概率密度函数估计是指利用统计学习方法从服从未知分布的观察样本中估计出随机变量的概率密度函数,常用的概率密度函数估计方法分为有参数和无参数两种。其中,有参数估计是指在概率密度函数已知的情况下,估计概率密度函数中的未知参数;无参数估计是指直接从观察样本中推断形式未知的概率密度函数。在实际电网中,电能质量参数的概率密度函数一般是未知的,有参数估计不适于其问题求解,采用无参数概率密度估计方法更为适宜。核密度估计方法是一种既直观形象又具有完备数学理论支持的无参数概率密度估计方法,可以真正地从观察样本出发研究随机变量的概率分布,能够得到一个渐进无偏的密度估计,有着良好的统计学性质。
本发明提供的一种分布式电源接入电网电能质量分析方法,可以对电压偏差、电压波动与闪变、谐波畸变、三相不平衡和频率波动等电能质量指标进行分析,估计出上述电能质量指标的核密度函数,进而从统计意义上对电能质量指标进行全面分析,如图1所示,包括:
101.获取电能质量指标的样本集合;
102.优化基于核密度函数的电能质量指标的概率密度函数模型中窗口宽度,其中,所述基于核密度函数的电能质量指标的概率密度函数模型根据所述样本集合确定;
103.利用所述电能质量指标的概率密度函数模型确定电能质量指标的统计特征值。
其中,所述电能质量指标包括:电压偏差、电压波动、电压闪变、三相电压不平衡、公用电网谐波和电力系统频率偏差。
具体的,供电电压偏差用dU表示:
电压波动用d表示:
闪变是电压波动在一段时间内的累积效果,用短时闪变Pst和长时闪变Plt来衡量。短时闪变Pst根据中华人民共和国国家标准GB/T 12326-2008《电能质量电压波动和闪变》附录A所给出的方法进行计算,长时闪变Plt由短时闪变计算而得:
三相电压不平衡用εU2和εU0表示
电压总谐波畸变率用THDU表示:
其中,
电流总谐波畸变率用THDI表示:
其中,
频率偏差用df表示:
df=f-fN
上式中,U为监测点的电压,UN为监测点的额定电压,ΔU为监测点的电压均方根曲线上相临两个极值电压之差,UN为监测点的额定电压,Pstj为2h内第j个闪变值,U1为三相电压正序分量方均根值,U2为三相电压负序分量方均根值,U0为三相电压零序分量方均根值,U1为基波电压方均根值,UH为谐波电压含量,Uh为第h次谐波电压方均根值,I1为基波电流方均根值,IH为谐波电流含量,Ih为第h次谐波电流方均根值,f为测量的配电网实际频率,fN为配电网的标称频率;
上述各项电能质量指标均为随机变量,可采用无参数的核密度估计方法求出其概率密度函数。因具体参数不同,需分别根据实际监测数据构建观察样本,再根据观察样本进行求解计算,但求解方法和过程一致。为不失一般性,用随机变量X表示某项电能质量指标。如果分析电压偏差,X代表dU;如果分析电压波动,X代表d;如果分析电压短时闪变,X代表Pst;如果分析电压长时闪变,X代表Plt;如果分析三相电压不平衡,X代表εU2或εU0;如果分析电压总谐波畸变率,X代表THDU;如果分析电流总谐波畸变率,X代表THDI;如果分析频率差,X代表df。
进一步的,采集电能质量指标之后,所述基于核密度函数的电能质量指标的概率密度函数模型根据所述样本集合确定,包括:
令电能质量指标为x,则其对应的样本集合为{x1,x2,…,xn},其中,n为采样点总数;
采用核密度函数估计随机变量x的概率密度函数,按下式确定基于核密度函数的电能质量指标的概率密度函数模型
上式中,h为窗口宽度,xi为所述样本集合中第i个采样点对应的电能质量指标值,K(u)为核函数,其中,K(u)中的u为自定义变量。
其中,可以采用Gaussian核函数,则K(u)的公式为:
则基于Gaussian核函数的电能质量指标的概率密度函数模型的公式为:
电能质量指标的概率密度函数模型中的窗口宽度h的选取决定核密度估计的表现,表现最优的窗口宽度是使所述电能质量指标的概率密度函数模型与实际电能质量指标的概率密度函数模型的误差为最小的窗口宽度,因此,所述步骤102,包括:
利用积分方差方式确定使所述电能质量指标的概率密度函数模型与实际电能质量指标的概率密度函数模型的误差为最小的窗口宽度的关于窗口宽度的积分方差表达式;
利用所述积分方差表达式,采用牛顿迭代算法获取最优窗口宽度。
其中,所述利用积分方差方式确定使所述电能质量指标的概率密度函数模型与实际电能质量指标的概率密度函数模型的误差为最小的窗口宽度的关于窗口宽度的积分方差表达式,包括:
按下式确定关于窗口宽度的积分方差表达式ISE*(h):
上式中,h为窗口宽度,xi为所述样本集合中第i个采样点对应的电能质量指标值,xj为所述样本集合中第i个采样点对应的电能质量指标值,n为采样点总数;
其中,使ISE*(h)最小的窗口宽度即为使所述电能质量指标的概率密度函数模型与实际电能质量指标的概率密度函数模型的误差为最小的窗口宽度。
ISE*(h)的推导过程如下:
采用积分方差(Integrated Square Error,ISE)确定使所述电能质量指标的概率密度函数模型与实际电能质量指标的概率密度函数模型的误差为最小的窗口宽度,ISE的数学表达式为:
上式中第三项与未知的h无关,因此使ISE(h)为最小的h也必定使下式取得最小值:
采用最小二乘交叉验证处理所述电能质量指标的概率密度函数模型,即用除xi之外的n-1个观察值得到的概率密度函数:
因而可得:
进而可得ISE*(h),为:
由于无法从上式求出使ISE*(h)为最小的最优窗口宽度的解析解,需采用优化算法求解最优窗口宽度,鉴于所求解的问题属于单变量连续非线性问题,经典优化算法即能对问题进行求解,求解思路是将求ISE*(h)的最小值问题转化为求ISE*(h)的导数为0的问题,把ISE*(h)按泰勒展开,忽略3次及3次以上的高次项,得到其2阶形式:
当且仅当Δh无限趋近于0时,ISE*(h+Δh)=ISE*(h),约去这两项,余项对Δh求导,上式等价为:
求解得:
进而得到迭代公式:
因此,如图2所示,所述利用所述积分方差表达式,采用牛顿迭代算法获取最优窗口宽度,包括:
a.设迭代次数t=0,设定窗口宽度初始值h0及收敛条件;
b.确定所述积分方差表达式对ht的一阶导数解析式和二阶导数解析式;
c.利用所述积分方差表达式对ht的一阶导数解析式和二阶导数解析式确定迭代窗口宽度ht+1;
d.若ht+1满足所述收敛条件,则ht+1为最优窗口宽度,若ht+1不满足所述收敛条件,则令t=t+1,ht=ht+1,并返回步骤b,其中,所述收敛条件包括:|ht+1-ht|≤ε和迭代次数到达T,t∈T。
其中,按下式确定所述积分方差表达式对ht的一阶导数解析式:
按下式确定所述积分方差表达式对ht的二阶导数解析式:
上式中,xi为所述样本集合中第i个采样点对应的电能质量指标值,xj为所述样本集合中第i个采样点对应的电能质量指标值,n为采样点总数;
按下式利用所述积分方差表达式对ht的一阶导数解析式和二阶导数解析式确定迭代窗口宽度ht+1:
获取电能质量指标的概率密度函数模型中最优窗口宽度后更新电能质量指标的概率密度函数模型,进一步利用更新后电能质量指标的概率密度函数模型电能质量指标的特征进行分析,所述步骤103,包括:
按下式确定电能质量指标的均值E(X):
按下式确定电能质量指标的方差D(X):
按下式确定电能质量指标的合格概率Pok:
按下式确定电能质量指标的超出下限值的概率PLL:
按下式确定电能质量指标的超出上限值的概率PUL:
上式中,(a,b)为电能质量指标的变化区域,Xmin为电能质量指标的下限值,Xmax为电能质量指标的上限值,为所述电能质量指标的概率密度函数模型;
将所述电能质量指标的概率密度函数模型的变化特性曲线中峰值所在的点作为电能质量指标的众数。
虽然电能质量指标的均值可以反映其数据集中的趋势,但实际情况下,电能质量指标很少服从正态分布,在统计分布上具有明显集中趋势点的数值并不在均值附近。根据核密度估计的极值点,可以求出其众数,即电能质量指标出现次数最多的数值,即使有多个众数,也能够求出。因此,基于核密度的电能质量分析方法能更好地反映其分布规律。
本发明还提供一种分布式电源接入电网电能质量分析装置,如图3所示,所述装置包括:
获取单元,用于获取电能质量指标的样本集合;
优化单元,用于优化基于核密度函数的电能质量指标的概率密度函数模型中窗口宽度,其中,所述基于核密度函数的电能质量指标的概率密度函数模型根据所述样本集合确定;
确定单元,用于利用所述电能质量指标的概率密度函数模型确定电能质量指标的统计特征值。
其中,所述电能质量指标包括:电压偏差、电压波动、电压闪变、三相电压不平衡、公用电网谐波和电力系统频率偏差。
所述确定单元,用于:
令电能质量指标为x,则其对应的样本集合为{x1,x2,…,xn},其中,n为采样点总数;
按下式确定基于核密度函数的电能质量指标的概率密度函数模型
上式中,h为窗口宽度,xi为所述样本集合中第i个采样点对应的电能质量指标值,K(u)为核函数,其中,K(u)中的u为自定义变量。
其中,采用Gaussian核函数,则K(u)的公式为:
则基于Gaussian核函数的电能质量指标的概率密度函数模型的公式为:
所述优化单元,用于:
利用积分方差方式确定使所述电能质量指标的概率密度函数模型与实际电能质量指标的概率密度函数模型的误差为最小的窗口宽度的关于窗口宽度的积分方差表达式;
利用所述积分方差表达式,采用牛顿迭代算法获取最优窗口宽度。
其中,所述利用积分方差方式确定使所述电能质量指标的概率密度函数模型与实际电能质量指标的概率密度函数模型的误差为最小的窗口宽度的关于窗口宽度的积分方差表达式,包括:
按下式确定关于窗口宽度的积分方差表达式ISE*(h):
上式中,h为窗口宽度,xi为所述样本集合中第i个采样点对应的电能质量指标值,xj为所述样本集合中第i个采样点对应的电能质量指标值,n为采样点总数;
其中,使ISE*(h)最小的窗口宽度即为使所述电能质量指标的概率密度函数模型与实际电能质量指标的概率密度函数模型的误差为最小的窗口宽度。
所述利用所述积分方差表达式,采用牛顿迭代算法获取最优窗口宽度,包括:
a.设迭代次数t=0,设定窗口宽度初始值h0及收敛条件;
b.确定所述积分方差表达式对ht的一阶导数解析式和二阶导数解析式;
c.利用所述积分方差表达式对ht的一阶导数解析式和二阶导数解析式确定迭代窗口宽度ht+1;
d.若ht+1满足所述收敛条件,则ht+1为最优窗口宽度,若ht+1不满足所述收敛条件,则令t=t+1,ht=ht+1,并返回步骤b,其中,所述收敛条件包括:|ht+1-ht|≤ε和迭代次数到达T,t∈T。
按下式确定所述积分方差表达式对ht的一阶导数解析式:
按下式确定所述积分方差表达式对ht的二阶导数解析式:
上式中,xi为所述样本集合中第i个采样点对应的电能质量指标值,xj为所述样本集合中第i个采样点对应的电能质量指标值,n为采样点总数;
按下式利用所述积分方差表达式对ht的一阶导数解析式和二阶导数解析式确定迭代窗口宽度ht+1:
所述确定单元,用于:
按下式确定电能质量指标的均值E(X):
按下式确定电能质量指标的方差D(X):
按下式确定电能质量指标的合格概率Pok:
按下式确定电能质量指标的超出下限值的概率PLL:
按下式确定电能质量指标的超出上限值的概率PUL:
上式中,(a,b)为电能质量指标的变化区域,Xmin为电能质量指标的下限值,Xmax为电能质量指标的上限值,为所述电能质量指标的概率密度函数模型;
将所述电能质量指标的概率密度函数模型的变化特性曲线中峰值所在的点作为电能质量指标的众数。
本领域内的技术人员应明白,本申请的实施例可提供为方法、系统、或计算机程序产品。因此,本申请可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且,本申请可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本申请是参照根据本申请实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器,使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品,该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上,使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理,从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
最后应当说明的是:以上实施例仅用以说明本发明的技术方案而非对其限制,尽管参照上述实施例对本发明进行了详细的说明,所属领域的普通技术人员应当理解:依然可以对本发明的具体实施方式进行修改或者等同替换,而未脱离本发明精神和范围的任何修改或者等同替换,其均应涵盖在本发明的权利要求保护范围之内。
Claims (10)
1.一种分布式电源接入电网电能质量分析方法,其特征在于,所述方法包括:
获取电能质量指标的样本集合;
优化基于核密度函数的电能质量指标的概率密度函数模型中窗口宽度,其中,所述基于核密度函数的电能质量指标的概率密度函数模型根据所述样本集合确定;
利用所述电能质量指标的概率密度函数模型确定电能质量指标的统计特征值。
2.如权利要求1所述的方法,其特征在于,所述电能质量指标包括:电压偏差、电压波动、电压闪变、三相电压不平衡、公用电网谐波和电力系统频率偏差。
3.如权利要求1所述的方法,其特征在于,所述基于核密度函数的电能质量指标的概率密度函数模型根据所述样本集合确定,包括:
令电能质量指标为x,则其对应的样本集合为{x1,x2,…,xn},其中,n为采样点总数;
按下式确定基于核密度函数的电能质量指标的概率密度函数模型
上式中,h为窗口宽度,xi为所述样本集合中第i个采样点对应的电能质量指标值,K(u)为核函数,其中,K(u)中的u为自定义变量。
4.如权利要求3所述的方法,其特征在于,采用Gaussian核函数,则K(u)的公式为:
<mrow>
<mi>K</mi>
<mrow>
<mo>(</mo>
<mi>u</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<msqrt>
<mrow>
<mn>2</mn>
<mi>&pi;</mi>
</mrow>
</msqrt>
</mfrac>
<mi>exp</mi>
<mrow>
<mo>(</mo>
<mo>-</mo>
<mfrac>
<msup>
<mi>u</mi>
<mn>2</mn>
</msup>
<mn>2</mn>
</mfrac>
<mo>)</mo>
</mrow>
<mo>,</mo>
<mi>u</mi>
<mo>&Element;</mo>
<mrow>
<mo>(</mo>
<mo>-</mo>
<mi>&infin;</mi>
<mo>,</mo>
<mo>+</mo>
<mi>&infin;</mi>
<mo>)</mo>
</mrow>
</mrow>
则基于Gaussian核函数的电能质量指标的概率密度函数模型的公式为:
5.如权利要求1所述的方法,其特征在于,所述优化所述电能质量指标的概率密度函数模型中窗口宽度,包括:
利用积分方差方式确定使所述电能质量指标的概率密度函数模型与实际电能质量指标的概率密度函数模型的误差为最小的窗口宽度的关于窗口宽度的积分方差表达式;
利用所述积分方差表达式,采用牛顿迭代算法获取最优窗口宽度。
6.如权利要求5所述的方法,其特征在于,所述利用积分方差方式确定使所述电能质量指标的概率密度函数模型与实际电能质量指标的概率密度函数模型的误差为最小的窗口宽度的关于窗口宽度的积分方差表达式,包括:
按下式确定关于窗口宽度的积分方差表达式ISE*(h):
<mrow>
<msup>
<mi>ISE</mi>
<mo>*</mo>
</msup>
<mrow>
<mo>(</mo>
<mi>h</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mrow>
<mn>2</mn>
<msqrt>
<mi>&pi;</mi>
</msqrt>
<mi>n</mi>
<mi>h</mi>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<mn>1</mn>
<mrow>
<mn>2</mn>
<msqrt>
<mi>&pi;</mi>
</msqrt>
<msup>
<mi>n</mi>
<mn>2</mn>
</msup>
<mi>h</mi>
</mrow>
</mfrac>
<mo>&lsqb;</mo>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<munderover>
<munder>
<mi>&Sigma;</mi>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
</munder>
<mrow>
<mi>j</mi>
<mo>&NotEqual;</mo>
<mi>i</mi>
</mrow>
<mi>n</mi>
</munderover>
<mo>&lsqb;</mo>
<mi>exp</mi>
<mo>&lsqb;</mo>
<mo>-</mo>
<mfrac>
<mn>1</mn>
<mn>4</mn>
</mfrac>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<msub>
<mi>x</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<msub>
<mi>x</mi>
<mi>j</mi>
</msub>
</mrow>
<mi>h</mi>
</mfrac>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>&rsqb;</mo>
<mo>&rsqb;</mo>
<mo>-</mo>
<mn>2</mn>
<msqrt>
<mn>2</mn>
</msqrt>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<munderover>
<munder>
<mi>&Sigma;</mi>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
</munder>
<mrow>
<mi>j</mi>
<mo>&NotEqual;</mo>
<mi>i</mi>
</mrow>
<mi>n</mi>
</munderover>
<mo>&lsqb;</mo>
<mi>exp</mi>
<mo>&lsqb;</mo>
<mo>-</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<msub>
<mi>x</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<msub>
<mi>x</mi>
<mi>j</mi>
</msub>
</mrow>
<mi>h</mi>
</mfrac>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>&rsqb;</mo>
<mo>&rsqb;</mo>
<mo>&rsqb;</mo>
</mrow>
上式中,h为窗口宽度,xi为所述样本集合中第i个采样点对应的电能质量指标值,xj为所述样本集合中第i个采样点对应的电能质量指标值,n为采样点总数;
其中,使ISE*(h)最小的窗口宽度即为使所述电能质量指标的概率密度函数模型与实际电能质量指标的概率密度函数模型的误差为最小的窗口宽度。
7.如权利要求5所述的方法,其特征在于,所述利用所述积分方差表达式,采用牛顿迭代算法获取最优窗口宽度,包括:
a.设迭代次数t=0,设定窗口宽度初始值h0及收敛条件;
b.确定所述积分方差表达式对ht的一阶导数解析式和二阶导数解析式;
c.利用所述积分方差表达式对ht的一阶导数解析式和二阶导数解析式确定迭代窗口宽度ht+1;
d.若ht+1满足所述收敛条件,则ht+1为最优窗口宽度,若ht+1不满足所述收敛条件,则令t=t+1,ht=ht+1,并返回步骤b,其中,所述收敛条件包括:|ht+1-ht|≤ε和迭代次数到达T,t∈T。
8.如权利要求7所述的方法,其特征在于,按下式确定所述积分方差表达式对ht的一阶导数解析式:
<mfenced open = "" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mfrac>
<mrow>
<mi>d</mi>
<mi>I</mi>
<mi>S</mi>
<mi>E</mi>
<mo>*</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>h</mi>
<mi>t</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>dh</mi>
<mi>t</mi>
</msub>
</mrow>
</mfrac>
<mo>=</mo>
<mo>-</mo>
<mfrac>
<mn>1</mn>
<mrow>
<mn>2</mn>
<msqrt>
<mi>&pi;</mi>
</msqrt>
<msup>
<msub>
<mi>nh</mi>
<mi>t</mi>
</msub>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<mo>-</mo>
<mfrac>
<mn>1</mn>
<mrow>
<mn>2</mn>
<msqrt>
<mi>&pi;</mi>
</msqrt>
<msup>
<mi>n</mi>
<mn>2</mn>
</msup>
<msup>
<msub>
<mi>h</mi>
<mi>t</mi>
</msub>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<mo>&lsqb;</mo>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<munderover>
<munder>
<mi>&Sigma;</mi>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
</munder>
<mrow>
<mi>j</mi>
<mo>&NotEqual;</mo>
<mi>i</mi>
</mrow>
<mi>n</mi>
</munderover>
<mo>&lsqb;</mo>
<mi>exp</mi>
<mo>&lsqb;</mo>
<mo>-</mo>
<mfrac>
<mn>1</mn>
<mn>4</mn>
</mfrac>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<msub>
<mi>x</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<msub>
<mi>x</mi>
<mi>j</mi>
</msub>
</mrow>
<msub>
<mi>h</mi>
<mi>t</mi>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>&rsqb;</mo>
<mo>&rsqb;</mo>
<mo>-</mo>
<mn>2</mn>
<msqrt>
<mn>2</mn>
</msqrt>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<munderover>
<munder>
<mi>&Sigma;</mi>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
</munder>
<mrow>
<mi>j</mi>
<mo>&NotEqual;</mo>
<mi>i</mi>
</mrow>
<mi>n</mi>
</munderover>
<mo>&lsqb;</mo>
<mi>exp</mi>
<mo>&lsqb;</mo>
<mo>-</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<msub>
<mi>x</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<msub>
<mi>x</mi>
<mi>j</mi>
</msub>
</mrow>
<msub>
<mi>h</mi>
<mi>t</mi>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>&rsqb;</mo>
<mo>&rsqb;</mo>
<mo>&rsqb;</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>+</mo>
<mfrac>
<mn>1</mn>
<mrow>
<mn>4</mn>
<msqrt>
<mi>&pi;</mi>
</msqrt>
<msup>
<mi>n</mi>
<mn>2</mn>
</msup>
<msup>
<msub>
<mi>h</mi>
<mi>t</mi>
</msub>
<mn>4</mn>
</msup>
</mrow>
</mfrac>
<mo>&lsqb;</mo>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<munderover>
<munder>
<mi>&Sigma;</mi>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
</munder>
<mrow>
<mi>j</mi>
<mo>&NotEqual;</mo>
<mi>i</mi>
</mrow>
<mi>n</mi>
</munderover>
<mo>&lsqb;</mo>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<msub>
<mi>x</mi>
<mi>j</mi>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mi>exp</mi>
<mo>&lsqb;</mo>
<mo>-</mo>
<mfrac>
<mn>1</mn>
<mn>4</mn>
</mfrac>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<msub>
<mi>x</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<msub>
<mi>x</mi>
<mi>j</mi>
</msub>
</mrow>
<msub>
<mi>h</mi>
<mi>t</mi>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>&rsqb;</mo>
<mo>&rsqb;</mo>
<mo>-</mo>
<mn>4</mn>
<msqrt>
<mn>2</mn>
</msqrt>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<munderover>
<munder>
<mi>&Sigma;</mi>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
</munder>
<mrow>
<mi>j</mi>
<mo>&NotEqual;</mo>
<mi>i</mi>
</mrow>
<mi>n</mi>
</munderover>
<mo>&lsqb;</mo>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<msub>
<mi>x</mi>
<mi>j</mi>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mi>exp</mi>
<mo>&lsqb;</mo>
<mo>-</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<msub>
<mi>x</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<msub>
<mi>x</mi>
<mi>j</mi>
</msub>
</mrow>
<msub>
<mi>h</mi>
<mi>t</mi>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>&rsqb;</mo>
<mo>&rsqb;</mo>
<mo>&rsqb;</mo>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
按下式确定所述积分方差表达式对ht的二阶导数解析式:
<mfenced open = "" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mfrac>
<mrow>
<msup>
<mi>d</mi>
<mn>2</mn>
</msup>
<mi>I</mi>
<mi>S</mi>
<mi>E</mi>
<mo>*</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>h</mi>
<mi>t</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msup>
<msub>
<mi>dh</mi>
<mi>t</mi>
</msub>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mrow>
<msqrt>
<mi>&pi;</mi>
</msqrt>
<msup>
<msub>
<mi>nh</mi>
<mi>t</mi>
</msub>
<mn>3</mn>
</msup>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<mn>1</mn>
<mrow>
<msqrt>
<mi>&pi;</mi>
</msqrt>
<msup>
<mi>n</mi>
<mn>2</mn>
</msup>
<msup>
<msub>
<mi>h</mi>
<mi>t</mi>
</msub>
<mn>3</mn>
</msup>
</mrow>
</mfrac>
<mo>&lsqb;</mo>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<munderover>
<munder>
<mi>&Sigma;</mi>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
</munder>
<mrow>
<mi>j</mi>
<mo>&NotEqual;</mo>
<mi>i</mi>
</mrow>
<mi>n</mi>
</munderover>
<mo>&lsqb;</mo>
<mi>exp</mi>
<mo>&lsqb;</mo>
<mo>-</mo>
<mfrac>
<mn>1</mn>
<mn>4</mn>
</mfrac>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<msub>
<mi>x</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<msub>
<mi>x</mi>
<mi>j</mi>
</msub>
</mrow>
<msub>
<mi>h</mi>
<mi>t</mi>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>&rsqb;</mo>
<mo>&rsqb;</mo>
<mo>-</mo>
<msup>
<mrow>
<mo>(</mo>
<msqrt>
<mn>2</mn>
</msqrt>
<mo>)</mo>
</mrow>
<mn>3</mn>
</msup>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<munderover>
<munder>
<mi>&Sigma;</mi>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
</munder>
<mrow>
<mi>j</mi>
<mo>&NotEqual;</mo>
<mi>i</mi>
</mrow>
<mi>n</mi>
</munderover>
<mo>&lsqb;</mo>
<mi>exp</mi>
<mo>&lsqb;</mo>
<mo>-</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<msub>
<mi>x</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<msub>
<mi>x</mi>
<mi>j</mi>
</msub>
</mrow>
<msub>
<mi>h</mi>
<mi>t</mi>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>&rsqb;</mo>
<mo>&rsqb;</mo>
<mo>&rsqb;</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<mfrac>
<mn>1</mn>
<mrow>
<mn>4</mn>
<msqrt>
<mi>&pi;</mi>
</msqrt>
<msup>
<mi>n</mi>
<mn>2</mn>
</msup>
<msup>
<msub>
<mi>h</mi>
<mi>t</mi>
</msub>
<mn>5</mn>
</msup>
</mrow>
</mfrac>
<mo>&lsqb;</mo>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<munderover>
<munder>
<mi>&Sigma;</mi>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
</munder>
<mrow>
<mi>j</mi>
<mo>&NotEqual;</mo>
<mi>i</mi>
</mrow>
<mi>n</mi>
</munderover>
<mo>&lsqb;</mo>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<msub>
<mi>x</mi>
<mi>j</mi>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mi>exp</mi>
<mo>&lsqb;</mo>
<mo>-</mo>
<mfrac>
<mn>1</mn>
<mn>4</mn>
</mfrac>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<msub>
<mi>x</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<msub>
<mi>x</mi>
<mi>j</mi>
</msub>
</mrow>
<msub>
<mi>h</mi>
<mi>t</mi>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>&rsqb;</mo>
<mo>&rsqb;</mo>
<mo>-</mo>
<mn>4</mn>
<msqrt>
<mn>2</mn>
</msqrt>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<munderover>
<munder>
<mi>&Sigma;</mi>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
</munder>
<mrow>
<mi>j</mi>
<mo>&NotEqual;</mo>
<mi>i</mi>
</mrow>
<mi>n</mi>
</munderover>
<mo>&lsqb;</mo>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<msub>
<mi>x</mi>
<mi>j</mi>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mi>exp</mi>
<mo>&lsqb;</mo>
<mo>-</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<msub>
<mi>x</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<msub>
<mi>x</mi>
<mi>j</mi>
</msub>
</mrow>
<msub>
<mi>h</mi>
<mi>t</mi>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>&rsqb;</mo>
<mo>&rsqb;</mo>
<mo>&rsqb;</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<mfrac>
<mn>1</mn>
<mrow>
<msqrt>
<mi>&pi;</mi>
</msqrt>
<msup>
<mi>n</mi>
<mn>2</mn>
</msup>
<msup>
<msub>
<mi>h</mi>
<mi>t</mi>
</msub>
<mn>5</mn>
</msup>
</mrow>
</mfrac>
<mo>&lsqb;</mo>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<munderover>
<munder>
<mi>&Sigma;</mi>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
</munder>
<mrow>
<mi>j</mi>
<mo>&NotEqual;</mo>
<mi>i</mi>
</mrow>
<mi>n</mi>
</munderover>
<mo>&lsqb;</mo>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<msub>
<mi>x</mi>
<mi>j</mi>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mi>exp</mi>
<mo>&lsqb;</mo>
<mo>-</mo>
<mfrac>
<mn>1</mn>
<mn>4</mn>
</mfrac>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<msub>
<mi>x</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<msub>
<mi>x</mi>
<mi>j</mi>
</msub>
</mrow>
<msub>
<mi>h</mi>
<mi>t</mi>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>&rsqb;</mo>
<mo>&rsqb;</mo>
<mo>-</mo>
<mn>4</mn>
<msqrt>
<mn>2</mn>
</msqrt>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<munderover>
<munder>
<mi>&Sigma;</mi>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
</munder>
<mrow>
<mi>j</mi>
<mo>&NotEqual;</mo>
<mi>i</mi>
</mrow>
<mi>n</mi>
</munderover>
<mo>&lsqb;</mo>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<msub>
<mi>x</mi>
<mi>j</mi>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mi>exp</mi>
<mo>&lsqb;</mo>
<mo>-</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<msub>
<mi>x</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<msub>
<mi>x</mi>
<mi>j</mi>
</msub>
</mrow>
<msub>
<mi>h</mi>
<mi>t</mi>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>&rsqb;</mo>
<mo>&rsqb;</mo>
<mo>&rsqb;</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>+</mo>
<mfrac>
<mn>1</mn>
<mrow>
<mn>8</mn>
<msqrt>
<mi>&pi;</mi>
</msqrt>
<msup>
<mi>n</mi>
<mn>2</mn>
</msup>
<msup>
<msub>
<mi>h</mi>
<mi>t</mi>
</msub>
<mn>7</mn>
</msup>
</mrow>
</mfrac>
<mo>&lsqb;</mo>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<munderover>
<munder>
<mi>&Sigma;</mi>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
</munder>
<mrow>
<mi>j</mi>
<mo>&NotEqual;</mo>
<mi>i</mi>
</mrow>
<mi>n</mi>
</munderover>
<mo>&lsqb;</mo>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<msub>
<mi>x</mi>
<mi>j</mi>
</msub>
<mo>)</mo>
</mrow>
<mn>4</mn>
</msup>
<mi>exp</mi>
<mo>&lsqb;</mo>
<mo>-</mo>
<mfrac>
<mn>1</mn>
<mn>4</mn>
</mfrac>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<msub>
<mi>x</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<msub>
<mi>x</mi>
<mi>j</mi>
</msub>
</mrow>
<msub>
<mi>h</mi>
<mi>t</mi>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>&rsqb;</mo>
<mo>&rsqb;</mo>
<mo>-</mo>
<mn>8</mn>
<msqrt>
<mn>2</mn>
</msqrt>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<munderover>
<munder>
<mi>&Sigma;</mi>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
</munder>
<mrow>
<mi>j</mi>
<mo>&NotEqual;</mo>
<mi>i</mi>
</mrow>
<mi>n</mi>
</munderover>
<mo>&lsqb;</mo>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<msub>
<mi>x</mi>
<mi>j</mi>
</msub>
<mo>)</mo>
</mrow>
<mn>4</mn>
</msup>
<mi>exp</mi>
<mo>&lsqb;</mo>
<mo>-</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<msub>
<mi>x</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<msub>
<mi>x</mi>
<mi>j</mi>
</msub>
</mrow>
<msub>
<mi>h</mi>
<mi>t</mi>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>&rsqb;</mo>
<mo>&rsqb;</mo>
<mo>&rsqb;</mo>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
上式中,xi为所述样本集合中第i个采样点对应的电能质量指标值,xj为所述样本集合中第i个采样点对应的电能质量指标值,n为采样点总数;
按下式利用所述积分方差表达式对ht的一阶导数解析式和二阶导数解析式确定迭代窗口宽度ht+1:
<mrow>
<msub>
<mi>h</mi>
<mrow>
<mi>t</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>h</mi>
<mi>t</mi>
</msub>
<mo>-</mo>
<mfrac>
<mrow>
<mfrac>
<mrow>
<mi>d</mi>
<mi>I</mi>
<mi>S</mi>
<mi>E</mi>
<mo>*</mo>
<mrow>
<mo>(</mo>
<mi>h</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>d</mi>
<mi>h</mi>
</mrow>
</mfrac>
<msub>
<mo>|</mo>
<mrow>
<mi>h</mi>
<mo>=</mo>
<msub>
<mi>h</mi>
<mi>t</mi>
</msub>
</mrow>
</msub>
</mrow>
<mrow>
<mfrac>
<mrow>
<msup>
<mi>d</mi>
<mn>2</mn>
</msup>
<mi>I</mi>
<mi>S</mi>
<mi>E</mi>
<mo>*</mo>
<mrow>
<mo>(</mo>
<mi>h</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msup>
<mi>dh</mi>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<msub>
<mo>|</mo>
<mrow>
<mi>h</mi>
<mo>=</mo>
<msub>
<mi>h</mi>
<mi>t</mi>
</msub>
</mrow>
</msub>
</mrow>
</mfrac>
<mo>.</mo>
</mrow>
9.如权利要求1所述的方法,其特征在于,所述利用所述电能质量指标的概率密度函数模型确定电能质量指标的统计特征值,包括:
按下式确定电能质量指标的均值E(X):
按下式确定电能质量指标的方差D(X):
按下式确定电能质量指标的合格概率Pok:
按下式确定电能质量指标的超出下限值的概率PLL:
按下式确定电能质量指标的超出上限值的概率PUL:
上式中,(a,b)为电能质量指标的变化区域,Xmin为电能质量指标的下限值,Xmax为电能质量指标的上限值,为所述电能质量指标的概率密度函数模型;
将所述电能质量指标的概率密度函数模型的变化特性曲线中峰值所在的点作为电能质量指标的众数。
10.一种分布式电源接入电网电能质量分析装置,其特征在于,所述装置包括:
获取单元,用于获取电能质量指标的样本集合;
优化单元,用于优化基于核密度函数的电能质量指标的概率密度函数模型中窗口宽度,其中,所述基于核密度函数的电能质量指标的概率密度函数模型根据所述样本集合确定;;
确定单元,用于利用所述电能质量指标的概率密度函数模型确定电能质量指标的统计特征值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710619336.XA CN107505519B (zh) | 2017-07-26 | 2017-07-26 | 一种分布式电源接入电网电能质量分析方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710619336.XA CN107505519B (zh) | 2017-07-26 | 2017-07-26 | 一种分布式电源接入电网电能质量分析方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107505519A true CN107505519A (zh) | 2017-12-22 |
CN107505519B CN107505519B (zh) | 2021-03-02 |
Family
ID=60690010
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710619336.XA Active CN107505519B (zh) | 2017-07-26 | 2017-07-26 | 一种分布式电源接入电网电能质量分析方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107505519B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108665376A (zh) * | 2018-04-06 | 2018-10-16 | 东北电力大学 | 一种空间负荷预测中确定元胞负荷最大值的核密度估计方法 |
CN109283416A (zh) * | 2018-11-20 | 2019-01-29 | 广东电网有限责任公司 | 一种敏感设备故障率的计算方法及装置 |
WO2020024210A1 (zh) * | 2018-08-02 | 2020-02-06 | 深圳大学 | 集成核密度估计器窗口参数优化方法、装置及终端设备 |
WO2020051789A1 (zh) * | 2018-09-12 | 2020-03-19 | 深圳大学 | 最小熵核密度估计器生成方法、装置和计算机可读存储介质 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20130158901A1 (en) * | 2011-12-19 | 2013-06-20 | Zafer Sahinoglu | Method and System for Detecting Unbalance in Power Grids |
CN104240140A (zh) * | 2014-06-04 | 2014-12-24 | 国家电网公司 | 一种基于状态估计的三相不平衡度概率评估方法 |
CN104504263A (zh) * | 2014-12-22 | 2015-04-08 | 国家电网公司 | 一种基于分布概率的光伏电站谐波水平评估方法 |
US20150278410A1 (en) * | 2014-04-01 | 2015-10-01 | Abb Technology Ag | Method for monitoring system variables of a distribution or transmission grid |
CN105403754A (zh) * | 2015-10-16 | 2016-03-16 | 云南电网有限责任公司曲靖供电局 | 一种基于核函数的谐波电流概率密度计算方法 |
CN105548809A (zh) * | 2015-12-18 | 2016-05-04 | 南瑞(武汉)电气设备与工程能效测评中心 | 一种电力系统谐波责任的确定方法 |
CN106651150A (zh) * | 2016-12-01 | 2017-05-10 | 全球能源互联网研究院 | 一种评估电网谐波影响的方法及装置 |
-
2017
- 2017-07-26 CN CN201710619336.XA patent/CN107505519B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20130158901A1 (en) * | 2011-12-19 | 2013-06-20 | Zafer Sahinoglu | Method and System for Detecting Unbalance in Power Grids |
US20150278410A1 (en) * | 2014-04-01 | 2015-10-01 | Abb Technology Ag | Method for monitoring system variables of a distribution or transmission grid |
CN104240140A (zh) * | 2014-06-04 | 2014-12-24 | 国家电网公司 | 一种基于状态估计的三相不平衡度概率评估方法 |
CN104504263A (zh) * | 2014-12-22 | 2015-04-08 | 国家电网公司 | 一种基于分布概率的光伏电站谐波水平评估方法 |
CN105403754A (zh) * | 2015-10-16 | 2016-03-16 | 云南电网有限责任公司曲靖供电局 | 一种基于核函数的谐波电流概率密度计算方法 |
CN105548809A (zh) * | 2015-12-18 | 2016-05-04 | 南瑞(武汉)电气设备与工程能效测评中心 | 一种电力系统谐波责任的确定方法 |
CN106651150A (zh) * | 2016-12-01 | 2017-05-10 | 全球能源互联网研究院 | 一种评估电网谐波影响的方法及装置 |
Non-Patent Citations (1)
Title |
---|
李国武等: "基于分布式电源接入的配电网电能质量分析与控制举措评估", 《电工技术》 * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108665376A (zh) * | 2018-04-06 | 2018-10-16 | 东北电力大学 | 一种空间负荷预测中确定元胞负荷最大值的核密度估计方法 |
CN108665376B (zh) * | 2018-04-06 | 2021-06-18 | 东北电力大学 | 一种空间负荷预测中确定元胞负荷最大值的核密度估计方法 |
WO2020024210A1 (zh) * | 2018-08-02 | 2020-02-06 | 深圳大学 | 集成核密度估计器窗口参数优化方法、装置及终端设备 |
WO2020051789A1 (zh) * | 2018-09-12 | 2020-03-19 | 深圳大学 | 最小熵核密度估计器生成方法、装置和计算机可读存储介质 |
CN109283416A (zh) * | 2018-11-20 | 2019-01-29 | 广东电网有限责任公司 | 一种敏感设备故障率的计算方法及装置 |
Also Published As
Publication number | Publication date |
---|---|
CN107505519B (zh) | 2021-03-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107505519B (zh) | 一种分布式电源接入电网电能质量分析方法及装置 | |
CN104156892B (zh) | 一种有源配电网电压跌落仿真与评估方法 | |
CN103825269B (zh) | 一种考虑电力系统功频静特性的快速概率潮流计算方法 | |
CN106709651B (zh) | 一种基于风险理论的电力系统安全性评估系统 | |
CN105790261B (zh) | 一种随机谐波潮流计算方法 | |
CN111797510A (zh) | 一种新能源场站短路比的计算方法及系统 | |
Qian et al. | N-1 static security assessment method for power grids with high penetration rate of renewable energy generation | |
CN113094887B (zh) | 移频电磁暂态仿真的优化方法、装置和电子设备 | |
CN110137948B (zh) | 一种双馈风电场次同步振荡风险的评估方法及装置 | |
WO2024104037A1 (zh) | 基于数学形态学和模式识别的直流电弧检测方法 | |
CN109830987A (zh) | 计及分布式光伏随机性的主动配电网概率稳定性分析方法 | |
CN113629880A (zh) | 一种基于数据中台的台区电压干扰装置检测方法及系统 | |
CN105514975A (zh) | 一种光伏发电系统的能效预测方法 | |
CN109861293B (zh) | 光伏不确定性对电力系统小信号稳定性影响的评估方法 | |
CN116050250A (zh) | 基于混合驱动模型的电网数字孪生模型建模方法及装置 | |
CN106022594A (zh) | 基于极值理论的电力系统在线运行安全风险评估方法 | |
CN113471972B (zh) | 一种用于电力系统的监测惯量的方法和系统 | |
CN113627041A (zh) | 基于多源数据融合的时变动态负荷在线建模方法及装置 | |
CN110896218B (zh) | 一种用于建立集合性居民负荷的谐波建模方法及系统 | |
CN111130098B (zh) | 一种含分布式电源配电网系统风险评估方法 | |
CN108038279B (zh) | 电网不对称故障下风机机电暂态恢复过程评价方法和装置 | |
Li et al. | Analytical Stochastic Power Flow for Power Systems with High Penetration of Renewable Energy | |
CN117494591B (zh) | 一种风电电能质量评估方法及装置 | |
CN111416350B (zh) | 区域电网灾变源荷匹配优化恢复方法及系统 | |
Chen et al. | Comprehensive evaluation of distributed power quality based on improved entropy method |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |