CN107133406A - 一种电力***静态电压稳定域边界的快速搜索方法 - Google Patents
一种电力***静态电压稳定域边界的快速搜索方法 Download PDFInfo
- Publication number
- CN107133406A CN107133406A CN201710327744.8A CN201710327744A CN107133406A CN 107133406 A CN107133406 A CN 107133406A CN 201710327744 A CN201710327744 A CN 201710327744A CN 107133406 A CN107133406 A CN 107133406A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- mtd
- voltage stability
- mtr
- 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
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
一种电力***静态电压稳定域边界的快速搜索方法,包括:确定静态电压稳定域边界通用数学模型;确定初始静态电压稳定域边界点;初始化静态电压稳定域边界的通用数学模型;预测相邻静态电压稳定域边界点校正预测点得到与初始静态电压稳定域边界点z0相邻的静态电压稳定域边界校正点将校正点作初值循环采用上述过程,搜索下一相邻的静态电压稳定域边界点,在循环过程中,将联立梯度矩阵和额外增加的方程中s替换为直至功率增长方向角β小于0和功率增长方向角β大于90°各出现一次结束。本发明构建的SVSR边界具有更高精度,可实现二维及高维有功注入空间中电力***静态电压稳定域边界的快速搜索。
Description
技术领域
本发明涉及一种静态电压稳定域边界搜索问题。特别是涉及一种电力***静态电压稳定域边界的快速搜索方法。
背景技术
区域电网互联规模不断扩大、高渗透率可再生能源大规模并网、电力电子设备广泛应用,在增强电力***输电能力、提高电力***运行灵活性和经济性的同时,也带来了功率注入波动性、潮流分布随机性等问题,加剧了电力***运行风险,对***电压稳定性带来了难以预估的隐患[1]-[5]。因此,研究高渗透率可再生能源并网和电力***电力电子化背景下的电压稳定性具有十分重要意义[6]-[7]。
目前,连续潮流(continuation power flow,CPF)法[8]仍是研究电力***静态电压稳定性的基本方法,该方法核心思想是根据确定的功率增长方向进行逐点迭代,计算***在满足电压稳定性前提下的最大负荷裕度,通过负荷裕度来评估***的电压稳定性。由于该算法实现容易,且具有较高的计算精度,因而在实际电力***中得到了广泛应用[9]-[11]。但随着可再生能源的大规模并网和需求侧响应技术的大规模应用,可再生能源和需求侧响应特有的间歇性与随机性增加了***功率增长方向的不确定性,仅仅依靠负荷预测和机组组合所确定的功率增长方向计算得到的负荷裕度已不能真实反映***的电压稳定性。若要评估***真实的电压稳定性,需针对可能出现的功率增长方向,重复调用CPF计算***的负荷裕度,这势必带来沉重的计算负担,不利于***静态电压稳定评估的效率。电力***的静态电压稳定域[12](static voltage stability region,SVSR)是描述确定网络拓扑结构和参数下,***具有静态电压稳定性的运行区域。通过所构建的SVSR,可准确评估出***在不确定性功率增长方向下的电压稳定性,避免了传统电压稳定分析方法在面临***负荷和电源出力不确定性时的不足,是一种分析和评估电力***电压稳定性的可靠方法[13]。
构建SVSR关键在于准确、快速搜索其稳定域边界,计算SVSR边界的方法主要包括拟合法[14]-[15]和近似法[16]-[18]两种。拟合法是由传统电压稳定分析方法发展而来的,根据SVSR 边界主要是由鞍节分岔(saddle node bifurcation,鞍节分岔点)点构成的特点,从基态潮流出发,重复调用CPF计算不同功率增长方向下的鞍节分岔点点,最终通过所搜索的鞍节分岔点点构成SVSR边界[14]-[15]。此类方法侧重于利用鞍节分岔点点来构建SVSR边界,具有较高的精度。但搜索鞍节分岔点的核心计算算法仍是CPF,因而在搜索整个SVSR边界时会带来大量的迭代运算,离线计算量庞大。近似法不再对SVSR边界点上的鞍节分岔点点进行逐点搜索,而是通过对SVSR边界的局部近似以提高SVSR构建速度。文献[16]-[18]通过对SVSR边界拓扑特性研究,发现SVSR局部边界呈现良好的线性特性,可采用超平面进行近似,在此基础上推导出能够近似表达SVSR边界的局部解析表达式。此类方法可根据实际需要,在满足计算精度前提下快速生成某一超平面,以近似替代SVSR的局部边界。该超平面可描述SVSR 边界的局部拓扑特性,但存在较大的保守性,不能保证SVSR边界的全局精度,且不同关键参数选取对SVSR边界拓扑结构影响较大,难以用统一解析式描述和获得通用性结论。目前,探究SVSR边界拓扑特性,在兼顾计算精度与计算效率基础上,构建SVSR仍是SVSR研究的难点。
发明内容
本发明所要解决的技术问题是,提供一种电力***静态电压稳定域边界的快速搜索方法。
本发明所采用的技术方案是:一种电力***静态电压稳定域边界的快速搜索方法,充分利用静态电压稳定域边界上相邻鞍节分岔点之间关联关系,提升搜索效率,具体包括如下步骤:
1)确定静态电压稳定域边界通用数学模型;
2)确定初始静态电压稳定域边界点;
3)初始化静态电压稳定域边界的通用数学模型;
4)预测相邻静态电压稳定域边界点
5)校正预测点得到与初始静态电压稳定域边界点z0相邻的静态电压稳定域边界校正点
6)将校正点作初值循环采用步骤4)~步骤5)搜索下一相邻的静态电压稳定域边界点,在循环过程中,将步骤4)第(3)步的联立梯度矩阵和额外增加的方程中s替换为直至功率增长方向角β小于0和功率增长方向角β大于90°各出现一次,结束。
步骤1)包括:
(1)将连续潮流方程与静态电压稳定域边界判据联立,构造描述给定功率增长方向d0所对应的静态电压稳定域边界点的数学模型如下:
式中,x为***静态状态变量的向量,g(x)=0为电力***常规潮流方程式;λ为***负荷裕度,d0为功率增长方向,L为潮流雅可比矩阵零特征值对应的左特征向量,JT(x)为潮流雅可比矩阵的转置矩阵;
(2)确定静态电压稳定域边界的通用数学模型如下:
式中,η和γ为新增功率参数,分别描述当功率增长方向d0在节点i、j二维有功注入空间内变化时,所对应的静态电压稳定域边界点在i、j节点上的有功功率变化量,ei和ej分别表示节点i和j有功注入功率单位变化的向量,ei和ej的长度与d0相同的向量,且分别仅有第i和第j个元素为1,其余元素皆为0;λ0d0为运行点在功率增长方向d0下的负荷裕度。
步骤2)包括:
(1)获取基础数据,包括***基态潮流状态变量、***拓扑结构;
(2)根据给定的功率增长方向d0,通过续潮流方程确定初始静态电压稳定域边界点对应的状态变量x0和负荷裕度λ0,从而获得初始静态电压稳定域边界点。
步骤3)包括:
(1)确定初始静态电压稳定域边界点处潮流雅克比矩阵J(x0)零特征值对应的左特征向量L0,设新增功率参数初值γ0=η0=0;
(2)确定式初始化静态电压稳定域边界的通用数学模型初值z0=[x0,L0,γ0,η0]T,
式中,x0为初始静态电压稳定域边界点对应的状态变量。
步骤4)包括:
(1)解静态电压稳定域边界的通用数学模型对各个变量偏导数的梯度矩阵如下:
(2)在梯度矩阵的基础上额外增加如下方程,用以指定切向量中的某一个分量的大小和方向,
定义功率增长方向角如下:
式中,为中对应于新增功率参数η的分量;e取±1;当e取+1时,指向功率增长方向角β增大的一侧,当e取-1则指向功率增长方向角β减小的一侧;J(x0)为初始静态电压稳定域边界点处潮流雅克比矩阵,JT(x)为潮流雅可比矩阵的转置矩阵,L0为J(x0)零特征值对应的左特征向量;
(3)联立梯度矩阵和额外增加的方程如下式,确定初始静态电压稳定域边界点的切向量如下式:
式中,s=[0 … 0 1],E=[0 … 0 e]T,二者长度与相等;
(4)设定步长σ,预测与初始静态电压稳定域边界点z0相邻的静态电压稳定域边界点如下式:
式中,为归一化后的切向量。
步骤5)包括:
(1)确定经过预测点且与切向量垂直的超平面,如下式:
式中z为自由变量;
(2)给出校正方程,以预测点为初值,采用牛顿法求解校正点如下:
本发明的一种电力***静态电压稳定域边界的快速搜索方法,可实现二维及高维有功注入空间中电力***静态电压稳定域边界的快速搜索,克服了近似性方法在构建SVSR边界存在保守性的不足,同时避免全局精确搜索SVSR边界时,采用CPF方法逐一迭代计算鞍节分岔点的计算负担。本发明的方法构建的SVSR边界具有更高精度,相对基于CPF的静态电压稳定域构建方法的计算效率,构建电力***静态电压稳定域计算时间大幅降低,显著提高了电力***电压稳定域的构建效率。
附图说明
图1是SVSR边界的快速搜索图示;
图2是本发明一种电力***静态电压稳定域边界的快速搜索方法的流程图;
图3是WECC-9测试***图;
图4是WECC-9***坐标轴为负荷节点的有功消耗的SVSR;
图5是SVSR边界的计算误差对比图。
具体实施方式
下面结合实施例和附图对本发明的一种电力***静态电压稳定域边界的快速搜索方法做出详细说明。
以下描述和附图充分地示出本发明的具体实施方案,以使本领域的技术人员能够实践它们。其他实施方案可以包括结构的、逻辑的、电气的、过程的以及其他的改变。实施例仅代表可能的变化。除非明确要求,否则单独的组件和功能是可选的,并且操作的顺序可以变化。一些实施方案的部分和特征可以被包括在或替换其他实施方案的部分和特征。本发明的实施方案的范围包括权利要求书的整个范围,以及权利要求书的所有可获得的等同物。在本文中,本发明的这些实施方案可以被单独地或总地用术语“发明”来表示,这仅仅是为了方便,并且如果事实上公开了超过一个的发明,不是要自动地限制该应用的范围为任何单个发明或发明构思。
电力***的静态电压稳定域是由***中所有鞍节分岔点点围成的多维空间,***运行在以鞍节分岔点点为边界构建的静态电压稳定域内是维持电压稳定的必要条件,因此SVSR边界的搜索可等效为鞍节分岔点点的搜索。
如图1、图2所示,本发明的一种电力***静态电压稳定域边界的快速搜索方法,其特征在于,充分利用静态电压稳定域边界上相邻鞍节分岔点之间关联关系,提升搜索效率,具体包括如下步骤:
1)确定静态电压稳定域边界通用数学模型;包括:
(1)将连续潮流方程与静态电压稳定域边界判据联立,构造描述给定功率增长方向d0所对应的静态电压稳定域边界点的数学模型如下:
式中,x为***静态状态变量的向量,g(x)=0为电力***常规潮流方程式;λ为***负荷裕度,d0为功率增长方向,L为潮流雅可比矩阵零特征值对应的左特征向量,JT(x)为潮流雅可比矩阵的转置矩阵;
(2)确定静态电压稳定域边界的通用数学模型如下式:
式中,η和γ为新增功率参数,分别描述当功率增长方向d0在节点i、j二维有功注入空间内变化时,所对应的静态电压稳定域边界点在i、j节点上的有功功率变化量,ei和ej分别表示节点i和j有功注入功率单位变化的向量,ei和ej的长度与d0相同的向量,且分别仅有第 i和第j个元素为1,其余元素皆为0;λ0d0为运行点在功率增长方向d0下的负荷裕度。
2)确定初始静态电压稳定域边界点;包括:
(1)获取基础数据,包括***基态潮流状态变量、***拓扑结构;
(2)根据给定的功率增长方向d0,通过续潮流方程确定初始静态电压稳定域边界点对应的状态变量x0和负荷裕度λ0,从而获得初始静态电压稳定域边界点。
3)初始化静态电压稳定域边界的通用数学模型;包括:
(1)确定初始静态电压稳定域边界点处潮流雅克比矩阵J(x0)零特征值对应的左特征向量L0,设新增功率参数初值γ0=η0=0;
(2)确定式初始化静态电压稳定域边界的通用数学模型初值z0=[x0,L0,γ0,η0]T,
式中,x0为初始静态电压稳定域边界点对应的状态变量。
4)预测相邻静态电压稳定域边界点包括:
(1)解静态电压稳定域边界的通用数学模型对各个变量偏导数的梯度矩阵如下:
(2)在梯度矩阵的基础上额外增加如下方程,用以指定切向量中的某一个分量的大小和方向,
定义功率增长方向角如下:
式中,Δη|z为中对应于新增功率参数η的分量;e取±1;当e取+1时,Δz|z指向功率增长方向角β增大的一侧,当e取-1则指向功率增长方向角β减小的一侧;J(x0)为初始静态电压稳定域边界点处潮流雅克比矩阵,JT(x)为潮流雅可比矩阵的转置矩阵,L0为J(x0)零特征值对应的左特征向量;
(3)联立梯度矩阵和额外增加的方程如下式,确定初始静态电压稳定域边界点的切向量如下式:
式中,s=[0 … 0 1],E=[0 … 0 e]T,二者长度与相等;
(4)设定步长σ,预测与初始静态电压稳定域边界点z0相邻的静态电压稳定域边界点如下式:
式中,为归一化后的切向量。
5)校正预测点得到与初始静态电压稳定域边界点z0相邻的静态电压稳定域边界校正点包括:
(1)确定经过预测点且与切向量垂直的超平面,如下式:
式中z为自由变量;
(2)给出校正方程,以预测点为初值,采用牛顿法求解校正点如下:
6)将校正点作初值循环采用步骤4)~步骤5)搜索下一相邻的静态电压稳定域边界点,在循环过程中,将步骤4)第(3)步的联立梯度矩阵和额外增加的方程中s替换为直至功率增长方向角β小于0和功率增长方向角β大于90°各出现一次,结束。
在采用校正点搜索下一相邻的静态电压稳定域边界点中,求校正点切向量时,将步骤4)第(3)步的联立梯度矩阵和额外增加的方程中s替换为即下一相邻的静态电压稳定域边界点预测环节中式(6)的增广梯度矩阵可直接用上一静态电压稳定域边界点校正环节中式(10)雅可比矩阵替代。以此类推,直至功率增长方向角β小于0,结束从初始静态电压稳定域边界点向β减小方向搜索静态电压稳定域边界点的过程。进一步,回到初始鞍节分岔点点,E中参数e取-1,再按上述步骤向β增大方向搜索静态电压稳定域边界点,直至β大于90°。
下面给出具体实例:
本实例是以搜索WECC3机9节点***SVSR边界为例,验证静态电压稳定域(SVSR)边界快速搜索方法的有效性,WECC3机9节点测试***如图3所示。
选择负荷节点5和7为电压稳定关键节点,在以节点5和7有功注入为坐标轴的二维有功负荷空间内采用本发明的方法搜索SVSR边界。如图4所示,将基态作为起始点,设定初始功率增长方向d0=[ΔS2,ΔS3,ΔS4,ΔS5,ΔS6,ΔS7,ΔS8,ΔS9]T=[0,0,0,0.9,0,1,0,0]T,功率增长引起的功率不平衡量由平衡节点所在发电机补偿,其中ΔSi为节点i的复功率,采用CPF沿d0追踪得图4中初始SNB点9,其坐标为(2.323,2.091)。根据初始SNB相关状态变量x0计算得L0,初始化z0=[x0,L0,γ0,η0]T,z0详细数值见表1,表1中θi和Vi分别为节点i的电压相角和幅值。
表1基于初始SNB采用预测-校正搜索下一SNB过程中相关参数计算结果
根据所选定关键节点5、7设置相应的单位向量e5、e7,选定预测步长σ=0.3、E中参数 e=1,计算得SNB点9的切向量Δz0=[Δx0,ΔL0,Δγ0,Δη0]T,其中Δγ0=-0.920 0、Δη0=1.000 0,对Δz0归一化后有Δγ0=-0.663 5、Δη0=0.720 3。得首个预测点其状态变量详细数值如表1所示。
将带入式(10)进行求解,得SNB点8处状态变量详细结果见表1。以作为初值采用本发明所提快速方法继续下一SNB点,可得SNB点7、6、 5、4、3、2、1、0,其坐标及γ和η在每次预测-校正过程的具体数值如表2所示。
表2β减小方向下的SNB搜索结果
编号 | ηpre | γpre | ηcor | γcor | 坐标 |
7 | 0.4423 | -0.4130 | 0.4188 | -0.4182 | (2.7421,1.6728) |
6 | 0.6206 | -0.6370 | 0.6167 | 0.6421 | (2.9401,1.4489) |
5 | 0.8103 | -0.8705 | 0.8062 | -0.8756 | (3.1295,1.2154) |
4 | 0.9907 | -1.1133 | 0.9861 | -1.1182 | (3.3095,0.9728) |
3 | 1.1609 | -1.3649 | 1.1559 | -1.3696 | (3.4792,0.7214) |
2 | 1.3198 | -1.6248 | 1.3144 | -1.6292 | (3.6377,0.4618) |
1 | 1.4667 | -1.8922 | 1.4608 | -18963 | (3.7841,0.1947) |
0 | 1.6006 | -2.1663 | 1.5942 | -2.1698 | (3.9176,-0.0788) |
其中点0位于图4中第一象限以外,有β≤0,标志着从初始SNB点开始在有功注入空间第 一象限内向β减小方向搜索SVSR临界点结束。此时需重新回到初始SNB点,在有功注入空 间第一象限内向β增大方向搜索SVSR临界点。因而重新读入初始SNB点信息,设置式(10)中参数e为-1,经上述预测-校正方法得首个预测和校正点和(对应图4中点10),其计算所得状态变量见表1,后续计算过程同上。当计算至SNB点19时,有β≥90°,结束计算, SNB点10至19的坐标及γ和η在每次预测-校正过程中的具体值如表3所示。将所有SNB点逐 一连接可得二维有功注入功率空间内的SVSR边界。
表3β增大方向下的SNB搜索结果
进一步图5对比了采用本发明的方法和CPF构建SVSR的计算精度。理论上,SVSR边界是有一系列SNB点构成的,对于任意一SVSR边界点,***潮流雅可比矩阵奇异,最小特征值为0,因此,可通过判断实际计算所得最小特征值与0之间的距离来判断本发明所提方法和CPF构建SVSR的精度。显然,除SNB点9外,采用本发明所提方法搜索所得SNB点的精度要远高于采用CPF追踪所得的SNB点。图中采用CPF追踪所得的SNB点3、5、6、 7、10、11、12、13、14和17的最小特征值均大于10-3,其最小值0.110×10-3出现在SNB点 16,最大值6.81×10-3出现在点SNB点18,所得平均最小特征值为2.63×10-3;而采用本发明所提方法计算所得的最小特征值平均值为7.25×10-4,且计算所得的所有最小特征值均小于 1.00×10-3,较CPF计算结果精度更高。
采用本发明所提方法和CPF搜索所得SVSR边界的计算时间(计算平台CPU IntelCore i5-4210H,主频2.9GHz,内存8GB)分别为13.9263s和1.6402s。结果表明:搜索SVSR边界上相同数量SNB点时,本发明的方法耗时约为CPF耗时的十分之一。
本发明的方法在同时兼顾了SVSR边界的精度与计算负担,高效率获得了高精度电力***静态电压稳定域边界。
以上实施例仅用以说明本发明的技术方案而非对其限制,尽管参照上述实施例对本发明进行了详细的说明,所属领域的普通技术人员依然可以对本发明的具体实施方式进行修改或者等同替换,这些未脱离本发明精神和范围的任何修改或者等同替换,均在申请待批的本发明的权利要求保护范围之内。
[参考文献]
[1]姜涛,李筱婧,李国庆,等,基于广域量测信息的负荷裕度灵敏度计算新方法[J].电工技术学报,2016,31(21):102-113.
[2]姜涛,贾宏杰,姜懿郎,等.跨区互联电网热稳定安全域边界近似方法[J].电工技术学报,2016,31(8):134-146.
[3]薛禹胜,雷兴,薛峰,等.关于风电不确定性对电力***影响的评述[J].中国电机工程学报,2014,34(29):5029-5040.
[4]竺炜,凌晨,张玉龙,等.基于电压水平约束的电网负荷安全域[J].中国电机工程学报,2016,36(23):6389-6400.
[5]陈厚合,李国庆,姜涛.计及静态电压稳定约束的交直流***可用输电能力[J].电网技术,2012,36(2):76-81.
[6]Yuan Haoyu,Jiang Tao,Jia Hongjie,et al.Real-time wide-area loadingmargin sensitivity(WALMS)in power systems[C].IEEE PES General Meeting,2015.
[7]Jiang Tao,Bai Lingquan,Jia Hongjie,et al.Identification of voltagestability critical injection region in bulk power systems based on therelative gain of voltage coupling[J].IET Generation Transmission&Distribution,2016,10(7):1495-1503.
[8]Chiang H D,Flueck AJ,Shah K S,et al.CPFLOW:a practical tool fortracing power system steady-state stationary behavior due to load andgeneration variations[J].IEEE Transactions on Power Systems,1995,10(2):623-634.
[9]Iba K,Suzuki H,Egawa M,et al.Calculation of critical loadingcondition with nose curve using homotopy continuation method[J].IEEETransactions on Power Systems,1991, 6(2):584-593.
[10]孙宏斌,李钦,张明晔,等.基于动态潮流方程的连续潮流模型与方法[J].中国电机工程学报,2011,31(7):77-82.
[11]赵晋泉,张伯明.改进连续潮流计算鲁棒性的策略研究[J].中国电机工程学报,2005, 25(22):7-11.
[12]Liu C C,Wu F F.Steady-state voltage stability regions of powersystems[C].IEEE Conference on Decision and Control,1984,6(1):23-31.
[13]余贻鑫,李鹏,孙强,等.电力***潮流可行域边界拓扑性质及边界算法[J].电力***自动化,2006,30(10):6-11.
[14]余贻鑫,宿吉峰,贾宏杰.电力大***电压稳定可行域可视化初探[J].电力***自动化,2001,25(22):1-5.
[15]苗伟威,贾宏杰,董泽寅.基于有功负荷注入空间静态电压稳定域的最小切负荷算法[J].中国电机工程学报,2012,32(16):44-41.
[16]王成山,许晓菲,余贻鑫,等.电力***电压稳定域的局部可视化描述及其应用[J].中国电机工程学报,2004,24(3):1-5.
[17]韩琪,余贻鑫,李慧玲,等.电力***注入空间静态电压稳定域边界的实用表达式 [J].中国电机工程学报,2005,25(5):8-14.
[18]韩琪,余贻鑫,贾宏杰,等.静态电压稳定域边界的非线性近似解析表达[J].电力***自动化,2005,29(11):10-14。
Claims (6)
1.一种电力***静态电压稳定域边界的快速搜索方法,其特征在于,充分利用静态电压稳定域边界上相邻鞍节分岔点之间关联关系,提升搜索效率,具体包括如下步骤:
1)确定静态电压稳定域边界通用数学模型
2)确定初始静态电压稳定域边界点;
3)初始化静态电压稳定域边界的通用数学模型;
4)预测相邻静态电压稳定域边界点
5)校正预测点得到与初始静态电压稳定域边界点z0相邻的静态电压稳定域边界校正点
6)将校正点作初值循环采用步骤4)~步骤5)搜索下一相邻的静态电压稳定域边界点,在循环过程中,将步骤4)第(3)步的联立梯度矩阵和额外增加的方程中s替换为直至功率增长方向角β小于0和功率增长方向角β大于90°各出现一次,结束。
2.根据权利要求1所述的一种电力***静态电压稳定域边界的快速搜索方法,其特征在于,步骤1)包括:
(1)将连续潮流方程与静态电压稳定域边界判据联立,构造描述给定功率增长方向d0所对应的静态电压稳定域边界点的数学模型如下:
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>g</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>&lambda;d</mi>
<mn>0</mn>
</msub>
<mo>=</mo>
<mn>0</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msup>
<mi>J</mi>
<mi>T</mi>
</msup>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
<mo>&CenterDot;</mo>
<mi>L</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msup>
<mi>L</mi>
<mi>T</mi>
</msup>
<mo>&CenterDot;</mo>
<mi>L</mi>
<mo>-</mo>
<mn>1</mn>
<mo>=</mo>
<mn>0</mn>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
式中,x为***静态状态变量的向量,g(x)=0为电力***常规潮流方程式;λ为***负荷裕度,d0为功率增长方向,L为潮流雅可比矩阵零特征值对应的左特征向量,JT(x)为潮流雅可比矩阵的转置矩阵;
(2)确定静态电压稳定域边界的通用数学模型如下:
<mrow>
<mi>&Phi;</mi>
<mrow>
<mo>(</mo>
<mi>z</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>g</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>&lambda;</mi>
<mn>0</mn>
</msub>
<msub>
<mi>d</mi>
<mn>0</mn>
</msub>
<mo>+</mo>
<msub>
<mi>&eta;e</mi>
<mi>i</mi>
</msub>
<mo>+</mo>
<msub>
<mi>&gamma;e</mi>
<mi>j</mi>
</msub>
<mo>=</mo>
<mn>0</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msup>
<mi>J</mi>
<mi>T</mi>
</msup>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
<mo>&CenterDot;</mo>
<mi>L</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msup>
<mi>L</mi>
<mi>T</mi>
</msup>
<mo>&CenterDot;</mo>
<mi>L</mi>
<mo>-</mo>
<mn>1</mn>
<mo>=</mo>
<mn>0</mn>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
式中,η和γ为新增功率参数,分别描述当功率增长方向d0在节点i、j二维有功注入空间内变化时,所对应的静态电压稳定域边界点在i、j节点上的有功功率变化量,ei和ej分别表示节点i和j有功注入功率单位变化的向量,ei和ej的长度与d0相同的向量,且分别仅有第i和第j个元素为1,其余元素皆为0;λ0d0为运行点在功率增长方向d0下的负荷裕度。
3.根据权利要求1所述的一种电力***静态电压稳定域边界的快速搜索方法,其特征在于,步骤2)包括:
(1)获取基础数据,包括***基态潮流状态变量、***拓扑结构;
(2)根据给定的功率增长方向d0,通过续潮流方程确定初始静态电压稳定域边界点对应的状态变量x0和负荷裕度λ0,从而获得初始静态电压稳定域边界点。
4.根据权利要求1所述的一种电力***静态电压稳定域边界的快速搜索方法,其特征在于,步骤3)包括:
(1)确定初始静态电压稳定域边界点处潮流雅克比矩阵J(x0)零特征值对应的左特征向量L0,设新增功率参数初值γ0=η0=0;
(2)确定式初始化静态电压稳定域边界的通用数学模型初值z0=[x0,L0,γ0,η0]T,
式中,x0为初始静态电压稳定域边界点对应的状态变量。
5.根据权利要求1所述的一种电力***静态电压稳定域边界的快速搜索方法,其特征在于,步骤4)包括:
(1)解静态电压稳定域边界的通用数学模型对各个变量偏导数的梯度矩阵如下:
<mrow>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>&Phi;</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>z</mi>
</mrow>
</mfrac>
<msub>
<mo>|</mo>
<msub>
<mi>z</mi>
<mn>0</mn>
</msub>
</msub>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>J</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<msub>
<mi>e</mi>
<mi>i</mi>
</msub>
</mtd>
<mtd>
<msub>
<mi>e</mi>
<mi>j</mi>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mfrac>
<mo>&part;</mo>
<mrow>
<mo>&part;</mo>
<mi>x</mi>
</mrow>
</mfrac>
<mo>&lsqb;</mo>
<msup>
<mi>J</mi>
<mi>T</mi>
</msup>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
<msub>
<mi>L</mi>
<mn>0</mn>
</msub>
<mo>&rsqb;</mo>
<msub>
<mo>|</mo>
<msub>
<mi>x</mi>
<mn>0</mn>
</msub>
</msub>
</mrow>
</mtd>
<mtd>
<mrow>
<msup>
<mi>J</mi>
<mi>T</mi>
</msup>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mrow>
<mn>2</mn>
<msubsup>
<mi>L</mi>
<mn>0</mn>
<mi>T</mi>
</msubsup>
</mrow>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>3</mn>
<mo>)</mo>
</mrow>
</mrow>
(2)在梯度矩阵的基础上额外增加如下方程,用以指定切向量中的某一个分量的大小和方向,
<mrow>
<mi>&Delta;</mi>
<mi>&eta;</mi>
<msub>
<mo>|</mo>
<msub>
<mi>z</mi>
<mn>0</mn>
</msub>
</msub>
<mo>=</mo>
<mi>e</mi>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>4</mn>
<mo>)</mo>
</mrow>
</mrow>
定义功率增长方向角如下:
<mrow>
<mi>&beta;</mi>
<mo>=</mo>
<mi>arctan</mi>
<mfrac>
<mrow>
<msub>
<mi>&lambda;</mi>
<mn>0</mn>
</msub>
<msub>
<mi>&Delta;P</mi>
<mi>j</mi>
</msub>
<mo>+</mo>
<mi>&gamma;</mi>
</mrow>
<mrow>
<msub>
<mi>&lambda;</mi>
<mn>0</mn>
</msub>
<msub>
<mi>&Delta;P</mi>
<mi>i</mi>
</msub>
<mo>+</mo>
<mi>&eta;</mi>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>5</mn>
<mo>)</mo>
</mrow>
</mrow>
式中,为中对应于新增功率参数η的分量;e取±1;当e取+1时,指向功率增长方向角β增大的一侧,当e取-1则指向功率增长方向角β减小的一侧;J(x0)为初始静态电压稳定域边界点处潮流雅克比矩阵,JT(x)为潮流雅可比矩阵的转置矩阵,L0为J(x0)零特征值对应的左特征向量;
(3)联立梯度矩阵和额外增加的方程如下式,确定初始静态电压稳定域边界点的切向量如下式:
<mrow>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>&Phi;</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>z</mi>
</mrow>
</mfrac>
<msub>
<mo>|</mo>
<msub>
<mi>z</mi>
<mn>0</mn>
</msub>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mi>s</mi>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>&CenterDot;</mo>
<mi>&Delta;</mi>
<mi>z</mi>
<msub>
<mo>|</mo>
<msub>
<mi>z</mi>
<mn>0</mn>
</msub>
</msub>
<mo>=</mo>
<mi>E</mi>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>6</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>&Delta;</mi>
<mi>z</mi>
<msub>
<mo>|</mo>
<msub>
<mi>z</mi>
<mn>0</mn>
</msub>
</msub>
<mo>=</mo>
<msup>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>&Phi;</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>z</mi>
</mrow>
</mfrac>
<msub>
<mo>|</mo>
<msub>
<mi>z</mi>
<mn>0</mn>
</msub>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mi>s</mi>
</mtd>
</mtr>
</mtable>
</mfenced>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mi>E</mi>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>7</mn>
<mo>)</mo>
</mrow>
</mrow>
式中,s=[0 … 0 1],E=[0 … 0 e]T,二者长度与相等;
(4)设定步长σ,预测与初始静态电压稳定域边界点z0相邻的静态电压稳定域边界点如下式:
<mrow>
<msubsup>
<mi>z</mi>
<mn>1</mn>
<mrow>
<mi>p</mi>
<mi>r</mi>
<mi>e</mi>
</mrow>
</msubsup>
<mo>=</mo>
<msub>
<mi>z</mi>
<mn>0</mn>
</msub>
<mo>+</mo>
<mi>&Delta;</mi>
<mover>
<mi>z</mi>
<mo>&OverBar;</mo>
</mover>
<msub>
<mo>|</mo>
<msub>
<mi>z</mi>
<mn>0</mn>
</msub>
</msub>
<mo>&CenterDot;</mo>
<mi>&sigma;</mi>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>8</mn>
<mo>)</mo>
</mrow>
</mrow>
式中,为归一化后的切向量。
6.根据权利要求1所述的一种电力***静态电压稳定域边界的快速搜索方法,其特征在于,步骤5)包括:
(1)确定经过预测点且与切向量垂直的超平面,如下式:
<mrow>
<mi>&Delta;</mi>
<msup>
<mover>
<mi>z</mi>
<mo>&OverBar;</mo>
</mover>
<mi>T</mi>
</msup>
<msub>
<mo>|</mo>
<msub>
<mi>z</mi>
<mn>0</mn>
</msub>
</msub>
<mo>&CenterDot;</mo>
<mrow>
<mo>(</mo>
<mi>z</mi>
<mo>-</mo>
<msub>
<mi>z</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mi>&sigma;</mi>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>9</mn>
<mo>)</mo>
</mrow>
</mrow>
式中z为自由变量;
(2)给出校正方程,以预测点为初值,采用牛顿法求解校正点如下:
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>&Phi;</mi>
<mrow>
<mo>(</mo>
<mi>z</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mn>0</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>&Delta;</mi>
<msup>
<mover>
<mi>z</mi>
<mo>&OverBar;</mo>
</mover>
<mi>T</mi>
</msup>
<msub>
<mo>|</mo>
<msub>
<mi>z</mi>
<mn>0</mn>
</msub>
</msub>
<mo>&CenterDot;</mo>
<mrow>
<mo>(</mo>
<mi>z</mi>
<mo>-</mo>
<msub>
<mi>z</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mi>&sigma;</mi>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>10</mn>
<mo>)</mo>
</mrow>
<mo>.</mo>
</mrow>
3
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710327744.8A CN107133406B (zh) | 2017-05-10 | 2017-05-10 | 一种电力***静态电压稳定域边界的快速搜索方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710327744.8A CN107133406B (zh) | 2017-05-10 | 2017-05-10 | 一种电力***静态电压稳定域边界的快速搜索方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107133406A true CN107133406A (zh) | 2017-09-05 |
CN107133406B CN107133406B (zh) | 2020-05-26 |
Family
ID=59731425
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710327744.8A Active CN107133406B (zh) | 2017-05-10 | 2017-05-10 | 一种电力***静态电压稳定域边界的快速搜索方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107133406B (zh) |
Cited By (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107529644A (zh) * | 2017-10-10 | 2018-01-02 | 东北电力大学 | 一种电力***静态电压稳定域边界线性近似方法 |
CN107947199A (zh) * | 2017-12-01 | 2018-04-20 | 东北电力大学 | 基于优化模型的电力***热稳定安全域边界的搜索方法 |
CN107968399A (zh) * | 2017-12-04 | 2018-04-27 | 国网辽宁省电力有限公司 | 一种快速搜索静态电压稳定域边界的方法 |
CN108376981A (zh) * | 2018-03-16 | 2018-08-07 | 东北电力大学 | 含柔性直流的电力***静态电压稳定域边界搜索方法 |
WO2019033820A1 (zh) * | 2017-08-17 | 2019-02-21 | 东北电力大学 | 一种电力***静态电压稳定域边界快速搜索的优化模型 |
CN110112725A (zh) * | 2019-04-25 | 2019-08-09 | 中国电力科学研究院有限公司 | 一种信息驱动的电网静态稳定域边界的构建方法及装置 |
CN111244963A (zh) * | 2020-03-13 | 2020-06-05 | 海仿(上海)科技有限公司 | 动态潮流计算方法、装置、设备和存储介质 |
CN111416348A (zh) * | 2020-04-21 | 2020-07-14 | 竺炜 | 全网静态电压安全及功角稳定的线路负荷安全域构建方法 |
CN111740425A (zh) * | 2020-06-30 | 2020-10-02 | 国网陕西省电力公司电力科学研究院 | 一种基于改进型连续潮流法的电力***静态电压稳定分析方法及*** |
CN112600201A (zh) * | 2020-12-07 | 2021-04-02 | 武汉大学 | 基于多维全纯嵌入法的高维静态电压稳定边界计算方法 |
CN113794204A (zh) * | 2021-08-19 | 2021-12-14 | 西安交通大学 | 电压稳定性分析方法、装置、电子设备及存储介质 |
CN114256868A (zh) * | 2021-11-16 | 2022-03-29 | 中国大唐集团科学技术研究院有限公司火力发电技术研究院 | 风电并网***宽频带振荡安全域构建方法 |
CN114352469A (zh) * | 2021-12-08 | 2022-04-15 | 国网江西省电力有限公司电力科学研究院 | 一种高阶水轮机调节***稳定域计算方法及*** |
CN116756198A (zh) * | 2023-05-08 | 2023-09-15 | 万合智慧(北京)能源科技有限公司 | 一种电力***静态电压稳定域边界的快速搜索方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2003134666A (ja) * | 2001-10-16 | 2003-05-09 | Takeo Kawamura | 非線形計画法による最適潮流計算システム |
CN101404412A (zh) * | 2008-11-05 | 2009-04-08 | 中国电力科学研究院 | 一种用于静态电压稳定性分析的方法 |
CN102055191A (zh) * | 2010-12-31 | 2011-05-11 | 重庆大学 | 一种考虑外网等值的静态电压稳定广域切负荷控制方法 |
CN104217090A (zh) * | 2013-06-04 | 2014-12-17 | 刘光晔 | 一种用泰勒级数解析计算电力***电压稳定临界点的方法 |
CN105514999A (zh) * | 2016-01-29 | 2016-04-20 | 江苏省电力公司电力经济技术研究院 | 一种电网静态电压稳定预防控制模型及其算法 |
CN105552960A (zh) * | 2015-12-31 | 2016-05-04 | 国网冀北电力有限公司电力科学研究院 | 风电场电力***的电压稳定分析方法及装置 |
-
2017
- 2017-05-10 CN CN201710327744.8A patent/CN107133406B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2003134666A (ja) * | 2001-10-16 | 2003-05-09 | Takeo Kawamura | 非線形計画法による最適潮流計算システム |
CN101404412A (zh) * | 2008-11-05 | 2009-04-08 | 中国电力科学研究院 | 一种用于静态电压稳定性分析的方法 |
CN102055191A (zh) * | 2010-12-31 | 2011-05-11 | 重庆大学 | 一种考虑外网等值的静态电压稳定广域切负荷控制方法 |
CN104217090A (zh) * | 2013-06-04 | 2014-12-17 | 刘光晔 | 一种用泰勒级数解析计算电力***电压稳定临界点的方法 |
CN105552960A (zh) * | 2015-12-31 | 2016-05-04 | 国网冀北电力有限公司电力科学研究院 | 风电场电力***的电压稳定分析方法及装置 |
CN105514999A (zh) * | 2016-01-29 | 2016-04-20 | 江苏省电力公司电力经济技术研究院 | 一种电网静态电压稳定预防控制模型及其算法 |
Non-Patent Citations (1)
Title |
---|
王刚 等: "静态电压稳定域边界的二次近似分析", 《中国电机工程学报》 * |
Cited By (23)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2019033820A1 (zh) * | 2017-08-17 | 2019-02-21 | 东北电力大学 | 一种电力***静态电压稳定域边界快速搜索的优化模型 |
US11050248B2 (en) | 2017-08-17 | 2021-06-29 | Northeast Electric Power University | Optimization model for quick track of SVSR boundary of power system |
CN107529644B (zh) * | 2017-10-10 | 2020-08-21 | 东北电力大学 | 一种电力***静态电压稳定域边界线性近似方法 |
CN107529644A (zh) * | 2017-10-10 | 2018-01-02 | 东北电力大学 | 一种电力***静态电压稳定域边界线性近似方法 |
CN107947199A (zh) * | 2017-12-01 | 2018-04-20 | 东北电力大学 | 基于优化模型的电力***热稳定安全域边界的搜索方法 |
CN107947199B (zh) * | 2017-12-01 | 2021-01-15 | 东北电力大学 | 基于优化模型的电力***热稳定安全域边界的搜索方法 |
CN107968399A (zh) * | 2017-12-04 | 2018-04-27 | 国网辽宁省电力有限公司 | 一种快速搜索静态电压稳定域边界的方法 |
CN108376981A (zh) * | 2018-03-16 | 2018-08-07 | 东北电力大学 | 含柔性直流的电力***静态电压稳定域边界搜索方法 |
CN110112725A (zh) * | 2019-04-25 | 2019-08-09 | 中国电力科学研究院有限公司 | 一种信息驱动的电网静态稳定域边界的构建方法及装置 |
CN110112725B (zh) * | 2019-04-25 | 2023-01-31 | 中国电力科学研究院有限公司 | 一种信息驱动的电网静态稳定域边界的构建方法及装置 |
CN111244963A (zh) * | 2020-03-13 | 2020-06-05 | 海仿(上海)科技有限公司 | 动态潮流计算方法、装置、设备和存储介质 |
CN111416348A (zh) * | 2020-04-21 | 2020-07-14 | 竺炜 | 全网静态电压安全及功角稳定的线路负荷安全域构建方法 |
CN111416348B (zh) * | 2020-04-21 | 2024-01-16 | 竺炜 | 全网静态电压安全及功角稳定的线路负荷安全域构建方法 |
CN111740425A (zh) * | 2020-06-30 | 2020-10-02 | 国网陕西省电力公司电力科学研究院 | 一种基于改进型连续潮流法的电力***静态电压稳定分析方法及*** |
CN112600201A (zh) * | 2020-12-07 | 2021-04-02 | 武汉大学 | 基于多维全纯嵌入法的高维静态电压稳定边界计算方法 |
CN112600201B (zh) * | 2020-12-07 | 2022-06-03 | 武汉大学 | 基于多维全纯嵌入法的高维静态电压稳定边界计算方法 |
CN113794204A (zh) * | 2021-08-19 | 2021-12-14 | 西安交通大学 | 电压稳定性分析方法、装置、电子设备及存储介质 |
CN113794204B (zh) * | 2021-08-19 | 2024-03-26 | 西安交通大学 | 电压稳定性分析方法、装置、电子设备及存储介质 |
CN114256868A (zh) * | 2021-11-16 | 2022-03-29 | 中国大唐集团科学技术研究院有限公司火力发电技术研究院 | 风电并网***宽频带振荡安全域构建方法 |
CN114256868B (zh) * | 2021-11-16 | 2023-01-24 | 中国大唐集团科学技术研究院有限公司火力发电技术研究院 | 风电并网***宽频带振荡安全域构建方法 |
CN114352469A (zh) * | 2021-12-08 | 2022-04-15 | 国网江西省电力有限公司电力科学研究院 | 一种高阶水轮机调节***稳定域计算方法及*** |
CN114352469B (zh) * | 2021-12-08 | 2024-02-09 | 国网江西省电力有限公司电力科学研究院 | 一种高阶水轮机调节***稳定域计算方法及*** |
CN116756198A (zh) * | 2023-05-08 | 2023-09-15 | 万合智慧(北京)能源科技有限公司 | 一种电力***静态电压稳定域边界的快速搜索方法 |
Also Published As
Publication number | Publication date |
---|---|
CN107133406B (zh) | 2020-05-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107133406A (zh) | 一种电力***静态电压稳定域边界的快速搜索方法 | |
Yuan et al. | Hybrid voltage stability assessment (VSA) for N− 1 contingency | |
CN104778367B (zh) | 基于单一状态断面的广域戴维南等值参数在线计算方法 | |
Chen et al. | Wide-area measurement-based voltage stability sensitivity and its application in voltage control | |
CN109494724B (zh) | 基于lu分解的大电网戴维南等值参数在线辨识方法 | |
CN113158573B (zh) | 一种基于深度学习的小水电群区域配网最优潮流获取方法 | |
CN108075480B (zh) | 一种交直流***的状态估计方法及*** | |
CN106655190A (zh) | 一种求解风电场概率最优潮流的方法 | |
CN105305439A (zh) | 一种考虑输入变量相关性的概率动态潮流计算方法及*** | |
CN104917174A (zh) | 一种针对主动配电网的静态电压稳定性判定方法 | |
Zang et al. | A robust state estimator for integrated electrical and heating networks | |
CN114935688B (zh) | 基于功率分段的电弧炉供电***谐波评估方法与*** | |
CN107749628A (zh) | 计及燃气机组无功调节及热电协调的多目标电压优化方法 | |
Mínguez et al. | State estimation sensitivity analysis | |
Gamm et al. | Solving several problems of power systems using spectral and singular analyses | |
CN106372440B (zh) | 一种并行计算的配电网自适应抗差状态估计方法及装置 | |
CN109698505B (zh) | 大电网静态电压稳定在线防控的调控量化映射计算方法 | |
CN107832959B (zh) | 一种考虑负荷特性和电源约束的电压稳定评估方法 | |
CN106356867B (zh) | 一种节点无功电压灵敏度的确定方法 | |
Hogan et al. | Towards effective clustering techniques for the analysis of electric power grids | |
CN109255541A (zh) | 一种基于最小二乘和一乘之和的配电网抗差状态估计方法 | |
Liu et al. | On expansion of estimated stability region: Theory, methodology, and application to power systems | |
CN108376981A (zh) | 含柔性直流的电力***静态电压稳定域边界搜索方法 | |
CN103390251A (zh) | 一种电力***状态估计的量测量权重设置方法 | |
Li et al. | Two-point estimate method for probabilistic optimal power flow computation including wind farms with correlated parameters |
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 |