CN104037759A - 电力***注入转移分布因子的抗差最小二乘估计方法 - Google Patents
电力***注入转移分布因子的抗差最小二乘估计方法 Download PDFInfo
- Publication number
- CN104037759A CN104037759A CN201410261491.5A CN201410261491A CN104037759A CN 104037759 A CN104037759 A CN 104037759A CN 201410261491 A CN201410261491 A CN 201410261491A CN 104037759 A CN104037759 A CN 104037759A
- Authority
- CN
- China
- Prior art keywords
- distribution factor
- robust
- matrix
- square
- electric power
- 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
Landscapes
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明公开了一种电力***注入转移分布因子的抗差最小二乘估计方法,包括:根据注入转移分布因子定义式,推导普通最小二乘估计矩阵解式运算形式,在此基础上利用中位数法求取抗差最小二乘估计迭代初值,开始抗差估计迭代过程,以线路相对有功潮流残差作为收敛条件,满足收敛域的注入转移分布因子矩阵输出作为可用的合理值。本发明能够克服传统的基于线路电抗参数的直流估计方法由于线路参数不准确等原因造成的对注入转移分布因子估计的偏差,通过对量测好、坏数据的不同处理,在普通最小二乘估计基础上进一步提高了估计精度,可以为电力***实时调度监控线路有功潮流提供准确的注入转移分布因子。
Description
技术领域
本发明涉及一种电力***注入转移分布因子基于量测的抗差最小二乘估计方法。
背景技术
注入转移分布因子在电力***实时调度,静态安全分析,阻塞管理等领域具有重要应用。当前传统的注入转移分布因子的计算方法为基于直流潮流模型的直流估计方法,该计算方法计算前需要设定线路电抗参数、节点间线路连络关系、电网功率平衡节点。要想准确计算注入转移分布因子值,必须具备准确的线路电抗参数。但随着运行时间增加,由于设备维护不及时、气候、管理等各方面因素,线路参数将发生偏移,如果没有及时更新线路电抗参数,直流估计方法计算所得的注入转移分布因子将不准确。
当电网中发生元件退出运行或线路短路开断等故障时,节点间连络关系改变,即电网拓扑结构改变。同样,如果节点间连络关系没有及时更新,注入转移分布因子的计算将出现错误。
基于直流潮流模型的注入转移分布因子估计方法计算时设定了电网中功率平衡节点,很多情况下,该设置与电网实际运行中的功率平衡策略并不一致(如电网中含有多台调频机组时,应该有多个平衡节点)。因此平衡节点设置的问题也可能导致注入转移分布因子计算结果不准确。
实时调度应用中,如果使用传统直流估计方法计算得到的注入转移分布因子不准确,那么在线路发生传输功率越限等故障时,电力***运行人员将不能及时发现故障并执行相应策略排除故障,这将威胁到电力***的运行安全性可靠性。
因此,基于直流潮流模型的注入转移分布因子估计方法,由于其参数更新滞后、网络拓扑结构不能及时更新、平衡节点设置与实际不符等问题,已经逐渐不适用于当前多变、复杂的电力***。
另外,在交流潮流计算中,可以通过功率方程推导出电力***灵敏度(包括注入转移分布因子),但由于该方法同样依赖线路参数形成的节点导纳矩阵,且计算复杂,在注入转移分布因子的计算中很少采用。
注入转移分布因子在电力***中的应用,主要包括实时调度、静态安全分析、阻塞管理等方面,这些应用以注入转移分布因子为纽带,根据节点有功注入快速获得传输线路上的有功潮流,并用以做出控制决策。如果所用注入转移分布因子不准确,那么将对决策行为的准确性产生影响,进而对电力***造成损失甚至破坏,如,当注入转移分布因子不准确时,由它计算得到的线路潮流可能不能反映线路传输功率越限故障,电力***运行人员将不能及时发现该故障并予以切除。
发明内容
本发明的目的就是为了解决上述问题,提供了电力***注入转移分布因子的抗差最小二乘估计方法,该方法利用当前电力***日趋丰富的量测数据,结合抗差最小二乘法提高对注入转移分布因子的估计准确度。从而使电力***运行人员可以根据贴近实际的支路潮流情况发出准确的实时调度指令。
为了实现上述目的,本发明采用如下技术方案:
一种电力***注入转移分布因子的抗差最小二乘估计方法,包括:
步骤(1):根据注入转移分布因子定义,结合电力***SCADA量测数据,推导注入转移分布因子的普通最小二乘解的矩阵运算形式;
步骤(2):在普通最小二乘法的基础上,引入抗差最小二乘法:在注入转移分布因子的普通最小二乘估计解的形式中加入等价权矩阵,利用中位数法求取抗差最小二乘估计迭代初值,开始迭代过程,以线路有功潮流相对残差作为收敛条件,满足收敛域的注入转移分布因子矩阵作为可用的合理值输出。
所述步骤(1)中注入转移分布因子定义式为:
式中,为注入转移分布因子矩阵中的一个元素,表示***中某节点n对线路k有功功率的影响;注入转移分布因子根本上是电力***灵敏度的概念,因此由偏微分形式定义,当分子分母取较小变化量时,偏微分可近似等于两时段功率变化量之比,即时间t附近节点n有功注入变化引起的线路k有功潮流的变化量与节点n有功注入变化量ΔPn(t)之比。
线路k有功潮流总的变化量为***内所有N个节点注入变化引起的,即:
其中ΔPk(t)为线路k在t时刻与其相邻时刻的有功功率量测数据之差,ΔPn(t)为***中某节点n在t时刻及其相邻时刻有功注入量测数据之差,即
ΔPk(t)=Pk(ti+1)-Pk(ti),ΔPn(t)=Pn(ti+1)-Pn(ti);
为注入转移分布因子矩阵元素之一,表示***中某节点n对线路k有功功率的影响。
所述步骤(1)中注入转移分布因子的普通最小二乘解的矩阵运算形式:
其中,Pn为m×n维节点有功注入数据矩阵,Pk为m×1维线路k上的有功功率量测数据矩阵。
所述步骤(2)中利用抗差最小二乘法估计注入转移分布因子矩阵的迭代形式为:
其中,为注入转移分布因子矩阵Mk第(i+1)次迭代结果,Pn为m×n维节点有功注入数据矩阵,为第i次迭代后形成的m×m维等价权矩阵,Pk为m×1维线路k上的有功功率量测数据矩阵,εi+1为根据Mk第(i+1)次迭代结果产生的残差向量。
所述第i次迭代后形成的m×m维等价权矩阵为:
其中,pj为第i次迭代后形成的等价权矩阵的第j个对角线元素,为标准化残差ε'中的第j个元素的绝对值,j∈[1,2,…,m];标准化残差ε0为单位权中误差;k0,k1为根据电力***数据特点预先设定的两个常数,通常可以设定k0∈[1,1.5],k1∈[2.5,3]。
所述步骤(2)中的收敛条件具体为:
其中,是通过注入转移分布因子计算得到的线路k上的有功功率,Pk为线路k有功量测数据,A%为设定的相对残差,当相对残差小于A%时,认为注入转移分布因子矩阵Mk为合理值,迭代停止将Mk输出。
本发明的有益效果:
1.本发明首先根据注入转移分布因子定义推导得出注入转移分布因子的最小二乘估计解的形式,推导过程中,考虑电力***连续多时段SCADA量测数据特点,给出估计注入转移分布因子的矩阵运算形式;其次,为了在普通最小二乘估计的基础上进一步提高注入转移分布因子估计精度,本发明引入了抗差最小二乘法,抗差最小二乘法继承了普通最小二乘法求取注入转移分布因子的解的矩阵运算形式,保持了估计方法运算的快速性、简洁性;最后,抗差法的提高之处在于以迭代的方式求取注入转移分布因子,每次迭代过程中以等价权矩阵来实现对良好数据和坏数据的不同处理(良好数据权值视残差大小保持为1或者相应减小,坏数据权值赋零),迭代步骤以相对残差大小设置收敛域,满足收敛条件时输出的注入转移分布因子矩阵视为合理值,这一灵活性是抗差法提高注入转移分布因子估计精度的关键所在。
2.基于量测的抗差法估计得到的注入转移分布因子,可以通过丰富的量测数据反映电力***当前准确的运行状态,对于提高电力***控制的准确性有极大帮助。如果在实际应用中获取实时SCADA量测数据并结合抗差法估计注入转移分布因子,将不存在直流估计方法中线路参数变化、拓扑结构过时、平衡节点设置不准确等问题。
3.以所掌握河南电网某六节点九线路***SCADA量测数据,通过与直流估计法、普通最小二乘法对比线路有功功率残差指标,本发明提出的抗差法估计注入转移分布因子的有效性和实用性已经得到验证。
附图说明
图1为本发明整体流程图;
图2为本发明中所用等价权函数示意图;
图3为本发明实例验证所使用的***拓扑结构图;
图4为本发明方法与直流估计法相对残差时序分布对比图;
图5为本发明方法与直流估计法相对残差概率密度分布对比图;
图6为本发明方法与普通最小二乘估计法相对残差时序分布对比图;
图7为本发明方法与普通最小二乘估计法相对残差概率密度分布对比图。
具体实施方式
下面结合附图与实施例对本发明作进一步说明。
如图1所示,电力***注入转移分布因子的抗差最小二乘估计方法,主要包括如下步骤:
步骤(1):根据注入转移分布因子定义,结合电力***SCADA量测数据,推导注入转移分布因子的普通最小二乘解的矩阵运算形式;
步骤(2):在普通最小二乘法的基础上,引入抗差最小二乘法:在注入转移分布因子的普通最小二乘估计解的形式中加入等价权矩阵,利用中位数法求取抗差最小二乘估计迭代初值,开始迭代过程,以线路有功潮流相对残差作为收敛条件,满足收敛域的注入转移分布因子矩阵作为可用的合理值输出。
步骤(1)中注入转移分布因子定义式:
式中,为注入转移分布因子矩阵中一个元素,表示了节点n的有功注入对线路k有功潮流的影响。注入转移分布因子根本上是电力***灵敏度的概念,因此由偏微分形式定义,当分子分母取较小变化量时,偏微分可近似等于两时段功率变化量之比,即时间t附近节点n有功注入变化引起的线路k有功潮流的变化量与节点n有功注入变化量ΔPn(t)之比。
重写(1)式表达形式,
各节点注入数据可以直接由SCADA量测获得,通过相邻时段节点注入数据做差,可以得到ΔPn(t),
ΔPn(t)=Pn(ti+1)-Pn(ti) (3)
然而表示线路k有功潮流由于节点n有功注入引起的变化量(线路有功潮流变化量由***内所有节点注入变化引起),因此不能直接由量测数据获得。线路k有功潮流可以通过量测数据直接获得,通过连续时段线路有功数据做差,可以得到线路k有功总变化量ΔPk(t),
ΔPk(t)=Pk(ti+1)-Pk(ti) (4)
线路k有功总变化量ΔPk(t)由***内所有N个节点有功注入变化引起,即
将(2)式的形式代入(5),(5)式可以写做,
将(3)式形式代入(6)式等号右边,(4)式形式代入(6)式等号左边得,
(4)式可以由以下两个方程组成的方程组做差得到,
同样的,(8)方程组可以由更多按时间序列列写的方程构成,即
方程组(9)包含由t1-tm时段,共m个表示线路潮流与节点注入关系的方程组成,该关系由注入转移分布因子决定。
将(9)式写成矩阵形式,
为简记(10)式,令线路有功向量Pk=[Pk(t1),Pk(t2),…,Pk(tm)]T,令Pn代表(10)式中间m×N维表示节点有功注入数据的矩阵,令注入转移分布因子向量则(10)式可以简记为,
Pk=Pn·Mk (11)
其中,矩阵Pk,Pn,可以通过筛选出SCADA量测中的线路有功潮流、节点有功注入直接形成。注入转移分布因子矩阵Mk为待求量。为求得Mk的值,可以通过求解如下最小二乘问题获得,
其中,
ε=Pn·Mk-Pk (13)
该最小二乘问题解式为,
按照注入转移分布因子定义推导可得解式(14),从连续多时段SCADA量测数据中挑选出线路k有功潮流数据、***内所有N个节点有功注入数据,形成矩阵Pk、Pn,便可快速估计出注入转移分布因子矩阵Mk。
另外,如果希望同时估计出***内所有节点对于多条线路的注入转移分布因子,可以将矩阵Pk由一列扩展为多列,每一列是一条线路的连续多时段有功潮流数据;矩阵Pn不变,仍然为***内N个节点的有功注入数据;注入转移分布因子矩阵Mk按照矩阵Pk的扩展规则相应扩展便可。
以上便是根据注入转移分布因子定义,按照最小二乘思想推导注入转移分布因子解式的步骤。
步骤(2)在普通最小二乘基础上引入抗差最小二乘估计:
式(13)为回归结果与量测数据之差,即回归残差,
ε=Pn·Mk-Pk
其中Pn为m×N维节点有功注入数据,Mk为N×1维注入转移分布因子矩阵,Pk为m×1维线路有功矩阵。由此可以假设,残差ε=[ε1,ε2,…,εm]T。
抗差最小二乘法是一种特殊的加权最小二乘法。它的特殊之处在于,其权值矩阵在每一次迭代计算中,根据上一次迭代得到的残差及等价权函数不断更新。
首先给出抗差最小二乘法迭代格式,
其中为注入转移分布因子矩阵Mk第(i+1)次迭代结果,为第i次迭代后形成的m×m维等价权矩阵,εi+1为根据Mk第(i+1)次迭代结果产生的残差向量。
抗差最小二乘法可以灵活处理良好数据与坏数据,其灵活性主要通过等价权矩阵实现。等价权矩阵为对角线矩阵,设其对角线元素用表示,j∈[1,2,…,m],因此等价权矩阵
按照IGGⅢ方案,通过等价权函数可以给出等价权矩阵对角线元素与残差ε之间的映射关系,
其中为标准化残差ε'中的第j个元素的绝对值,j∈[1,2,…,m]。标准化残差ε0为单位权中误差。k0,k1为我们根据电力***数据特点预先设定的两个常数,通常可以设定k0∈[1,1.5],k1∈[2.5,3]。
通过等价权函数,等价权矩阵对角线元素可以由标准化残差ε'中对应元素获得。残差值越大,相应权值越小,坏数据权值赋零。
为了计算标准化残差进而得到等价权矩阵,需要有单位权中误差ε0的值,这里我们采用中位数法计算ε0。
如果我们有多组连续时段SCADA量测数据,可以取出其中q组,首先通过普通最小二乘法计算得出注入转移分布因子Mk的q组值,将所得Mk代入式(13),得到q组残差ε。对于每一个残差向量ε=[ε1,ε2,…,εm],通过如下取中位数运算获得一个εmed值,
q组残差向量ε将得到q个中位数εmed,取q个εmed中的最小值,记作εmed0,则单位权中误差估值为,
ε0=1.483εmed0 (18)
单位权中误差已通过中位数法成功求取,以ε/εmed0取绝对值后的值作为标准化残差初值,由(16)式获得等价权矩阵初值,代入抗差法迭代解式(15),即可进入上述迭代求解过程。
以上所述便是抗差最小二乘法估计注入转移分布因子矩阵Mk的迭代步骤。然而迭代停止需要设置收敛条件,残差表示回归方程与两侧值之差,因此我们可以使用相对残差值来作为判断收敛的条件,
其中是通过注入转移分布因子计算得到的线路k有功,Pk为线路k有功量测数据,等式左边即相对残差,例如当相对残差小于5%时,我们便可认为注入转移分布因子Mk为合理值,迭代停止将Mk输出。
以上便是为了进一步提高估计精度引入的注入转移分布因子抗差最小二乘估计法的计算步骤。
为直观展示抗差法中关键的等价权矩阵的形成原理,如图2所示为以标准化残差为自变量,权值为因变量的等价权函数示意图,由图可见,标准化残差元素越大,对应的量测数据被赋予的权值越小,绝对值超过k1的标准化残差对应的量测数据被视为坏数据,权值赋零,相当于将坏数据剔除。
抗差估计求取注入转移分布因子及其有效性验证,主要包括如下步骤:
1):注入转移分布因子的抗差估计应用在如图3所示的实际***之中,该***包含6个节点,9条线路,待求量注入转移分布因子矩阵M表示6个节点的有功注入(发电机注入为正,负荷视作负注入)对9条线路上的有功潮流的影像。我们已经掌握约1天的实际SCADA量测数据,数据文件按照每五分钟一次,时序排列。
2):明确了待求***的拓扑结构,具备了充足的量测数据,我们便可以按照抗差法估计注入转移分布因子的步骤来进行计算。SCADA量测数据为整个电网的多种量测数据,我们需要的仅是节点有功注入数据及线路有功潮流,因此首先可以使用MATLAB编写数据文件的读写程序,筛选出节点有功注入数据及线路有功潮流数据,形成mat数据文件,存储备用。
3):中位数法求出单位权中误差,为抗差法迭代过程求取等价权矩阵初值,从存储的mat数据文件中提取对应的节点有功注入与线路有功潮流数据,代入(15)式的迭代形式,开始迭代过程。
4):当求得的注入转移分布因子矩阵M满足诸如(19)式的收敛条件时,将该M矩阵输出,作为注入转移分布因子的合理值。
5):为了对比本发明提出的抗差估计方法与传统的直流估计法的准确性,同时证明在普通最小二乘的基础上引入抗差最小二乘估计的必要性,我们使用已知线路电抗参数,按照(3)式求出直流估计法得到注入转移分布因子矩阵M',使用量测数据,按照(14)式求出普通最小二乘法得到的注入转移分布因子矩阵M″。
6):将三种方法求得的注入转移分布因子矩阵M,M',M″代入(13)式,得到三种方法的三组线路有功潮流残差。根据三组残差绝对值,求得残差相对值,即相对残差,对比三组相对残差指标,便可以反映三种方法求得的注入转移分布因子准确度。分别做出直流估计法与抗差法残差的时序分布图(以线路3为例,见图4),概率密度分布对比图(图5);普通最小二乘与抗差最小二乘残差的时序分布图(以线路1为例,图6),概率密度分布对比图(图7)。
7):图4与图5为直流估计法与抗差估计相对残差对比图,由图4相对残差时序分布可见,直流估计得到的相对残差大于抗差估计得到的相对残差,由图5相对残差概率密度对比图可见,抗差估计得到的相对残差分布比直流估计的相对残差更加集中分布于0附近,图4、图5均表明本发明提出的抗差估计用于计算注入转移分布因子的方法优于现有的直流估计方法,本发明有效性得到验证;图6与图7为普通最小二乘估计与抗差估计相对残差对比图,同样,由图6相对残差时序分布对比图,图7相对残差概率密度分布对比图可见,抗差估计得到的注入转移分布因子更优,在普通最小二乘估计的基础上引入抗差估计的必要性也得到验证。
上述虽然结合附图对本发明的具体实施方式进行了描述,但并非对本发明保护范围的限制,所属领域技术人员应该明白,在本发明的技术方案的基础上,本领域技术人员不需要付出创造性劳动即可做出的各种修改或变形仍在本发明的保护范围以内。
Claims (7)
1.一种电力***注入转移分布因子的抗差最小二乘估计方法,其特征是,包括:
步骤(1):根据注入转移分布因子定义,结合电力***SCADA量测数据,推导注入转移分布因子的普通最小二乘解的矩阵运算形式;
步骤(2):在普通最小二乘法的基础上,引入抗差最小二乘法:在注入转移分布因子的普通最小二乘估计解的形式中加入等价权矩阵,利用中位数法求取抗差最小二乘估计迭代初值,开始迭代过程,以线路有功潮流相对残差作为收敛条件,将满足收敛域的注入转移分布因子矩阵作为可用的合理值输出。
2.如权利要求1所述的一种电力***注入转移分布因子的抗差最小二乘估计方法,其特征是,所述步骤(1)中注入转移分布因子定义式为:
式中,为注入转移分布因子矩阵中的一个元素,表示***中某节点n对线路k有功功率的影响;注入转移分布因子根本上是电力***灵敏度的概念,因此由偏微分形式定义,当分子分母取较小变化量时,偏微分可近似等于两时段功率变化量之比,即时间t附近节点n有功注入变化引起的线路k有功潮流的变化量与节点n有功注入变化量ΔPn(t)之比。
3.如权利要求2所述的一种电力***注入转移分布因子的抗差最小二乘估计方法,其特征是,线路k有功潮流总的变化量为***内所有N个节点注入变化引起的,即:
其中ΔPk(t)为线路k在t时刻与其相邻时刻的有功功率量测数据之差,ΔPn(t)为***中某节点n在t时刻及其相邻时刻有功注入量测数据之差,即
ΔPk(t)=Pk(ti+1)-Pk(ti),ΔPn(t)=Pn(ti+1)-Pn(ti);
为注入转移分布因子矩阵元素之一,表示***中某节点n对线路k有功功率的影响。
4.如权利要求1所述的一种电力***注入转移分布因子的抗差最小二乘估计方法,其特征是,所述步骤(1)中注入转移分布因子的普通最小二乘解的矩阵运算形式:
其中,Pn为m×n维节点有功注入数据矩阵,Pk为m×1维线路k上的有功功率量测数据矩阵。
5.如权利要求1所述的一种电力***注入转移分布因子的抗差最小二乘估计方法,其特征是,所述步骤(2)中利用抗差最小二乘法估计注入转移分布因子矩阵的迭代形式为:
其中,为注入转移分布因子矩阵Mk第(i+1)次迭代结果,Pn为m×n维节点有功注入数据矩阵,为第i次迭代后形成的m×m维等价权矩阵,Pk为m×1维线路k上的有功功率量测数据矩阵,εi+1为根据Mk第(i+1)次迭代结果产生的残差向量。
6.如权利要求5所述的一种电力***注入转移分布因子的抗差最小二乘估计方法,其特征是,所述第i次迭代后形成的m×m维等价权矩阵为:
其中,pj为第i次迭代后形成的等价权矩阵的第j个对角线元素,为标准化残差ε'中的第j个元素的绝对值,j∈[1,2,…,m];标准化残差ε0为单位权中误差;k0,k1为根据电力***数据特点预先设定的两个常数,通常可以设定k0∈[1,1.5],k1∈[2.5,3]。
7.如权利要求1所述的一种电力***注入转移分布因子的抗差最小二乘估计方法,其特征是,所述步骤(2)中的收敛条件具体为:
其中,是通过注入转移分布因子计算得到的线路k上的有功功率,Pk为线路k有功量测数据,A%为设定的相对残差,当相对残差小于A%时,认为注入转移分布因子矩阵Mk为合理值,迭代停止将Mk输出。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410261491.5A CN104037759B (zh) | 2014-06-12 | 2014-06-12 | 电力***注入转移分布因子的抗差最小二乘估计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410261491.5A CN104037759B (zh) | 2014-06-12 | 2014-06-12 | 电力***注入转移分布因子的抗差最小二乘估计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104037759A true CN104037759A (zh) | 2014-09-10 |
CN104037759B CN104037759B (zh) | 2016-01-20 |
Family
ID=51468423
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410261491.5A Active CN104037759B (zh) | 2014-06-12 | 2014-06-12 | 电力***注入转移分布因子的抗差最小二乘估计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104037759B (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104410108A (zh) * | 2014-11-26 | 2015-03-11 | 清华大学 | 用于联络线功率控制的发电转移分布因子快速计算方法 |
CN105119276A (zh) * | 2015-08-18 | 2015-12-02 | 河海大学 | 一种基于admm的电力***wlav抗差估计方法 |
CN107092775A (zh) * | 2017-03-17 | 2017-08-25 | 沈阳工业大学 | 基于激光测量***的锥螺纹数据处理算法 |
CN107359610A (zh) * | 2017-07-25 | 2017-11-17 | 国网山东省电力公司青岛供电公司 | 一种输电网注入转移分布因子估计方法及设备 |
CN109253718A (zh) * | 2018-11-22 | 2019-01-22 | 广东省重工建筑设计院有限公司 | 隧道变形自动化监测***和方法 |
CN112421636A (zh) * | 2020-10-31 | 2021-02-26 | 国网河南省电力公司漯河供电公司 | 一种计算电力***注入转移分布因子的方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20040181491A1 (en) * | 2003-03-12 | 2004-09-16 | Hitachi, Ltd. | Method, computer equipment and a program for planning of electric power generation and electric power trade |
CN103413043A (zh) * | 2013-08-09 | 2013-11-27 | 中国南方电网有限责任公司 | 一种电力***多维空间断面热稳定极限边界的求解方法 |
CN103514374A (zh) * | 2013-09-24 | 2014-01-15 | 清华大学 | 电力***在线滚动调度中不可行传输断面约束的辨识方法 |
-
2014
- 2014-06-12 CN CN201410261491.5A patent/CN104037759B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20040181491A1 (en) * | 2003-03-12 | 2004-09-16 | Hitachi, Ltd. | Method, computer equipment and a program for planning of electric power generation and electric power trade |
CN103413043A (zh) * | 2013-08-09 | 2013-11-27 | 中国南方电网有限责任公司 | 一种电力***多维空间断面热稳定极限边界的求解方法 |
CN103514374A (zh) * | 2013-09-24 | 2014-01-15 | 清华大学 | 电力***在线滚动调度中不可行传输断面约束的辨识方法 |
Non-Patent Citations (2)
Title |
---|
周登波: "准稳态灵敏度法在LMP计算中的应用", 《中国优秀硕士学位论文全文数据库(电子期刊)》, 26 July 2010 (2010-07-26) * |
沈瑜等: "发电联合转移分布因子及快速静态安全校核算法", 《电力***自动化》, vol. 27, no. 18, 25 September 2003 (2003-09-25) * |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104410108A (zh) * | 2014-11-26 | 2015-03-11 | 清华大学 | 用于联络线功率控制的发电转移分布因子快速计算方法 |
CN105119276A (zh) * | 2015-08-18 | 2015-12-02 | 河海大学 | 一种基于admm的电力***wlav抗差估计方法 |
CN107092775A (zh) * | 2017-03-17 | 2017-08-25 | 沈阳工业大学 | 基于激光测量***的锥螺纹数据处理算法 |
CN107092775B (zh) * | 2017-03-17 | 2020-06-05 | 沈阳工业大学 | 基于激光测量***的锥螺纹数据处理方法 |
CN107359610A (zh) * | 2017-07-25 | 2017-11-17 | 国网山东省电力公司青岛供电公司 | 一种输电网注入转移分布因子估计方法及设备 |
CN109253718A (zh) * | 2018-11-22 | 2019-01-22 | 广东省重工建筑设计院有限公司 | 隧道变形自动化监测***和方法 |
CN109253718B (zh) * | 2018-11-22 | 2020-11-17 | 广东省重工建筑设计院有限公司 | 隧道变形自动化监测***和方法 |
CN112421636A (zh) * | 2020-10-31 | 2021-02-26 | 国网河南省电力公司漯河供电公司 | 一种计算电力***注入转移分布因子的方法 |
Also Published As
Publication number | Publication date |
---|---|
CN104037759B (zh) | 2016-01-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104037759A (zh) | 电力***注入转移分布因子的抗差最小二乘估计方法 | |
Liu et al. | Modeling of local false data injection attacks with reduced network information | |
CN103840452B (zh) | 一种引入pmu量测信息的大电网状态估计方法 | |
CN102982379B (zh) | 基于集合论估计模型的电力***状态估计方法及*** | |
CN103944165B (zh) | 一种大电网参数辨识估计方法 | |
CN106707061A (zh) | 基于混合量测的配电网动态状态估计方法 | |
CN102508086B (zh) | 一种t形连接电力线路的参数估计方法 | |
CN106383296A (zh) | 一种基于相量分析的改进阻抗型有源配电网故障测距算法 | |
CN104009475B (zh) | 考虑等值元件全面性与参数物理约束的非拓扑静态等值法 | |
CN101666825B (zh) | 一种基于电压稳定性局部指标的电网电压稳定在线监测方法 | |
CN104362638B (zh) | 基于相量量测单元量测的电网聚合的关键节点电压调控方法 | |
CN105633956A (zh) | 一种基于Spiking神经网络伪量测建模的配电网三相状态估计方法 | |
CN101505061A (zh) | 一种基于时域仿真的戴维南等值参数跟踪的计算方法 | |
CN110289613A (zh) | 基于灵敏度矩阵的配电网拓扑识别与线路参数辨识方法 | |
CN102364490B (zh) | 基于层次分析模型的自动同调识别方法 | |
CN105375484A (zh) | 一种基于pmu的电力***分布式动态状态估计方法 | |
CN108376983A (zh) | 一种基于边界融合的电网状态估计***及方法 | |
CN102005758A (zh) | 一种基于内网实测信息的两端口外网静态等值方法 | |
CN105787604B (zh) | 一种输配协同的配网态势快速感知方法 | |
CN104617575B (zh) | 一种用于静态安全分析预想故障下潮流不可解评估方法 | |
CN105811404A (zh) | 一种输配协同的配网静态电压稳定态势监测方法 | |
CN102570457B (zh) | 一种基于内网量测数据的两端口外网静态等值方法 | |
CN105846437A (zh) | 一种计及互动相关性的cholesky分解半不变量潮流计算方法 | |
CN101350524A (zh) | 一种基于全微分方程的戴维南等值参数跟踪的计算方法 | |
CN107765179A (zh) | 一种适用于量测丢失的发电机动态状态估计方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |