CN107968409A - 一种考虑不平衡功率分配的概率潮流计算方法及*** - Google Patents
一种考虑不平衡功率分配的概率潮流计算方法及*** Download PDFInfo
- Publication number
- CN107968409A CN107968409A CN201711092861.7A CN201711092861A CN107968409A CN 107968409 A CN107968409 A CN 107968409A CN 201711092861 A CN201711092861 A CN 201711092861A CN 107968409 A CN107968409 A CN 107968409A
- Authority
- CN
- China
- Prior art keywords
- msub
- mrow
- cumulant
- power
- distribution
- 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
- 238000009826 distribution Methods 0.000 title claims abstract description 124
- 238000004364 calculation method Methods 0.000 title claims abstract description 32
- 238000000034 method Methods 0.000 claims abstract description 25
- 230000035945 sensitivity Effects 0.000 claims description 31
- 239000013598 vector Substances 0.000 claims description 22
- 230000004044 response Effects 0.000 claims description 13
- 238000005315 distribution function Methods 0.000 claims description 11
- 239000011159 matrix material Substances 0.000 claims description 10
- 238000005206 flow analysis Methods 0.000 claims description 6
- 230000008859 change Effects 0.000 claims description 4
- 239000000284 extract Substances 0.000 claims description 4
- 238000010606 normalization Methods 0.000 claims description 4
- 230000006870 function Effects 0.000 description 18
- 238000004458 analytical method Methods 0.000 description 9
- 238000004590 computer program Methods 0.000 description 7
- 238000010586 diagram Methods 0.000 description 6
- 238000012545 processing Methods 0.000 description 6
- 240000002853 Nelumbo nucifera Species 0.000 description 5
- 235000006508 Nelumbo nucifera Nutrition 0.000 description 5
- 235000006510 Nelumbo pentapetala Nutrition 0.000 description 5
- 239000000243 solution Substances 0.000 description 5
- 238000004422 calculation algorithm Methods 0.000 description 3
- 230000005611 electricity Effects 0.000 description 3
- 238000002347 injection Methods 0.000 description 3
- 239000007924 injection Substances 0.000 description 3
- 238000000342 Monte Carlo simulation Methods 0.000 description 2
- 238000013459 approach Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 238000010206 sensitivity analysis Methods 0.000 description 2
- 238000003860 storage Methods 0.000 description 2
- 238000006467 substitution reaction Methods 0.000 description 2
- 230000001052 transient effect Effects 0.000 description 2
- 230000009471 action Effects 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000005684 electric field Effects 0.000 description 1
- 239000004744 fabric Substances 0.000 description 1
- 238000009472 formulation Methods 0.000 description 1
- 238000012886 linear function Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 230000009897 systematic effect Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
Classifications
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/04—Circuit arrangements for ac mains or ac distribution networks for connecting networks of the same frequency but supplied from different sources
- H02J3/06—Controlling transfer of power between connected networks; Controlling sharing of load between connected networks
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0635—Risk analysis of enterprise or organisation activities
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
- H02J2203/20—Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
Landscapes
- Business, Economics & Management (AREA)
- Engineering & Computer Science (AREA)
- Human Resources & Organizations (AREA)
- Economics (AREA)
- Strategic Management (AREA)
- Health & Medical Sciences (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Physics & Mathematics (AREA)
- Marketing (AREA)
- General Business, Economics & Management (AREA)
- Entrepreneurship & Innovation (AREA)
- Tourism & Hospitality (AREA)
- Water Supply & Treatment (AREA)
- Primary Health Care (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Development Economics (AREA)
- Educational Administration (AREA)
- Power Engineering (AREA)
- Game Theory and Decision Science (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明提供了一种考虑不平衡功率分配的概率潮流计算方法及***,包括:获取基础数据;根据获取的基础数据计算母线负荷的半不变量和新能源机组出力的半不变量;根据母线负荷的半不变量和新能源机组出力的半不变量求取考虑注入功率波动、不平衡功率分配的电压和支路潮流的概率分布。本发明提供的技术方案解决了传统半不变量法概率潮流中不平衡功率由平衡节点承担的不合理问题,提高了半不变量法概率潮流的准确性和实用性。
Description
技术领域
本发明属于电力***自动化技术领域,具体涉及一种考虑不平衡功率分配的概率潮流计算方法及***。
背景技术
发电计划和停电计划的制定,需要通过未来态潮流和安全分析保证电网的安全稳定运行。目前,以未来态潮流和安全分析为核心的安全校核已经在我国各级调控中心得到广泛应用,并取得了良好的应用效果。但现有技术都是基于确定性的潮流模型,对于每一个分析断面,电源、网络以及负荷都是未来某一时刻的“期望值”,而未能考虑各种不确定因素引起的潮流分布变化以及电网安全风险。
概率潮流能够统筹考虑各种不确定性因素,便于发现电网中的薄弱环节,是电力***稳态运行分析的重要工具。概率潮流计算方法主要包括蒙特卡罗模拟法、点估计法和解析法,其中蒙特卡罗模拟法通常需要几千次以上的仿真计算,常作为其它方法结果的验证;点估计法算法简单,但获取的随机变量的高阶矩不够精确,不易处理随机变量的相关性;解析法包括卷积法和半不变量法,卷积法计算量大、耗时长,半不变量法将独立随机变量的卷积运算转化为简单的加减法,计算效率较高。
在传统的半不变量法概率潮流计算中将负荷及电源的随机波动产生的不平衡功率由平衡节点承担,不符合电网实际运行情况。因此急需一种针对交流潮流模型下电源、负荷的随机变动所引起的不平衡功率分配的概率潮流计算方法。
发明内容
本发明提出一种考虑不平衡功率分配的概率潮流计算方法,对于负荷波动等连续的随机变量,以准稳态灵敏度构建支路潮流与节点注入功率的线性关系,灵活的处理不平衡功率在各节点的分摊方式,从而得到更为符合电网实际运行的***潮流概率分布。
本发明提供的一种考虑不平衡功率分配的概率潮流计算方法,包括:
获取基础数据;
根据获取的基础数据计算母线负荷的半不变量和新能源机组出力的半不变量;
根据母线负荷的半不变量和新能源机组出力的半不变量求取考虑注入功率波动、不平衡功率分配的电压和支路潮流的概率分布;
所述基础数据包括:机组出力及机组出力概率分布描述、母线负荷期望及母线负荷概率分布描述。
所述计算母线负荷的半不变量包括:
计算母线负荷的各阶原点矩;
根据所述各阶原点矩,计算母线负荷的各阶半不变量。
所述计算新能源机组出力的半不变量包括:
计算新能源机组的各阶原点矩;
根据所述原点矩计算新能源机组出力的各阶半不变量。
所述求取考虑注入功率波动、不平衡功率分配的电压和支路潮流的概率分布包括:
基于准稳态灵敏度线性化潮流方程;
根据半不变量和Gram-Charlier级数展开理论,由注入功率的各阶半不变量计算节点电压和支路潮流的概率分布。
所述基于准稳态灵敏度线性化潮流方程的表达式如下所示:
式中,X为节点电压向量,Z为支路潮流向量,X0为节点电压向量的基准运行状态,Z0为支路潮流向量基准运行状态,ΔW0为注入功率的变化量,和分别为节点电压和支路潮流对注入功率变化的准稳态灵敏度,Fu为准稳态物理响应线性化后的系数矩阵,S0与T0分别为节点电压和支路潮流对注入功率变化的常规灵敏度。
所述准稳态物理响应线性化后的系数矩阵Fu如下所示:
其中,αn是n维列向量,n为***节点数,αi是第i节点的功率不平衡量分配系数,若αi=0,表示i节点不参与不平衡功率分担。
所述根据半不变量和Gram-Charlier级数展开理论,由注入功率的各阶半不变量计算节点电压和支路潮流的概率分布包括:
基于准稳态灵敏度线性化潮流方程,按照下式由注入功率各阶半不变量计算节点电压和支路潮流的各阶半不变量:
式中ΔWk是负荷各阶半不变量与新能源机组出力各阶半不变量相加得到的注入功率各阶半不变量,分别是中元素的k次幂矩阵;
节点电压和支路潮流随机变量规格化按下式表示:
其中,ξ表示节点电压和支路潮流随机变量,μ表示其均值,即一阶半不变量,σ为标准差,即二阶半不变量开平方;
节点电压和支路潮流的概率密度函数f(x)和概率分布函数F(x)按下式表示:
式中,和φ(x)分别是期望值为0,标准差为1的标准正态分布随机变量的概率密度函数和概率分布函数,cr为常系数,计算公式如下所示:
所述基础数据包括还电网模型和设备参数。
本发明提供一种考虑不平衡功率分配的概率潮流计算***,包括:
采集模块,用于获取基础数据;
第一计算模块,用于根据获取的基础数据计算母线负荷的半不变量和新能源机组出力的半不变量;
第二计算模块,用于根据母线负荷的半不变量和新能源机组出力的半不变量求取考虑注入功率波动、不平衡功率分配的电压和支路潮流的概率分布
所述基础数据包括:机组出力及机组出力概率分布描述、母线负荷期望及母线负荷概率分布描述。
所述第二计算模块包括:线性化子模块和计算子模块;
所述线性化子模块基于准稳态灵敏度线性化潮流方程;
所述计算子模块根据半不变量和Gram-Charlier级数展开理论,由注入功率的各阶半不变量计算节点电压和支路潮流的概率分布。
本发明提供一种考虑不平衡功率分配的概率潮流分析方法,包括:
根据权利要求1所述的考虑不平衡功率分配的概率潮流计算方法求取节点电压和支路潮流的概率分布;
根据所述节点电压和支路潮流的概率分布,确定是否存在潮流越限风险,若存在,则调整机组出力,否则输出节点电压和支路潮流的概率分布。
本发明提供一种考虑不平衡功率分配的概率潮流分析***,包括:
计算模块,用于根据权利要求1所述的考虑不平衡功率分配的概率潮流计算方法求取节点电压和支路潮流的概率分布;
判断模块,用于根据所述节点电压和支路潮流的概率分布,确定是否存在潮流越限风险,若存在,则调整机组出力,否则输出节点电压和支路潮流的概率分布。
与最接近的现有技术比,本发明提供的技术方案具有以下有益效果:
本发明提供的技术方案,采用半不变量法描述机组出力、负荷的随机波动特性,避免复杂的卷积运算;
本发明提供的技术方案,以准稳态灵敏度构建节点注入功率与支路潮流的线性关系,灵活合理地分配注入功率随机波动带来的不平衡功率;
本发明提供的技术方案,采用半不变量与级数展开相结合的方法,得到符合电网实际运行情况的概率潮流分布。
本发明提供的技术方案,解决了传统半不变量法概率潮流中不平衡功率由平衡节点承担的不合理问题,提高了半不变量法概率潮流的准确性和实用性。
附图说明
图1为本发明一种考虑不平衡功率分配的概率潮流计算方法流程图;
图2为本发明实施例基于考虑不平衡功率分配的概率潮流计算方法的概率潮流分析调整方法流程图。
具体实施方式
下面结合附图对本发明做进一步详细的说明:
如图1所示,本发明提供一种考虑不平衡功率分配的概率潮流计算方法,包括:
获取基础数据;
根据获取的基础数据计算母线负荷的半不变量和新能源机组出力的半不变量;
根据母线负荷的半不变量和新能源机组出力的半不变量求取考虑注入功率波动、不平衡功率分配的电压和支路潮流的概率分布;
所述基础数据包括:机组出力及机组出力概率分布描述、母线负荷期望及母线负荷概率分布描述。
所述基础数据包括电网模型、设备参数、机组出力及其概率分布描述、母线负荷期望及其概率分布描述。
所述计算母线负荷的半不变量包括:
计算母线负荷的各阶原点矩;
根据所述各阶原点矩,计算母线负荷的各阶半不变量。
所述计算新能源机组出力的半不变量包括:
计算新能源机组的各阶原点矩;
根据所述原点矩计算新能源机组出力的各阶半不变量。
所述求取考虑注入功率波动、不平衡功率分配的电压和支路潮流的概率分布包括:
基于准稳态灵敏度线性化潮流方程;
根据半不变量和Gram-Charlier级数展开理论,由注入功率的各阶半不变量计算节点电压和支路潮流的概率分布。
所述基于准稳态灵敏度线性化潮流方程的表达式如下所示:
式中,X为节点电压向量,Z为支路潮流向量,X0为节点电压向量的基准运行状态,Z0为支路潮流向量基准运行状态,ΔW0为注入功率的变化量,和分别为节点电压和支路潮流对注入功率变化的准稳态灵敏度,Fu为准稳态物理响应线性化后的系数矩阵,S0与T0分别为节点电压和支路潮流对注入功率变化的常规灵敏度。
所述准稳态物理响应线性化后的系数矩阵Fu如下所示:
其中,αn是n维列向量,n为***节点数,αi是第i节点的功率不平衡量分配系数,若αi=0,表示i节点不参与不平衡功率分担。
所述根据半不变量和Gram-Charlier级数展开理论,由注入功率的各阶半不变量计算节点电压和支路潮流的概率分布包括:
基于准稳态灵敏度线性化潮流方程,按照下式由注入功率各阶半不变量计算节点电压和支路潮流的各阶半不变量:
式中ΔWk是负荷各阶半不变量与新能源机组出力各阶半不变量相加得到的注入功率各阶半不变量,分别是中元素的k次幂矩阵;
节点电压和支路潮流随机变量规格化按下式表示:
其中,ξ表示节点电压和支路潮流随机变量,μ表示其均值,即一阶半不变量,σ为标准差,即二阶半不变量开平方;
节点电压和支路潮流的概率密度函数f(x)和概率分布函数F(x)按下式表示:
式中,和φ(x)分别是期望值为0,标准差为1的标准正态分布随机变量的概率密度函数和概率分布函数,cr为常系数,计算公式如下所示:
基于相同的构思本发明提供一种考虑不平衡功率分配的概率潮流计算***,可以包括:
采集模块,用于获取基础数据;
第一计算模块,用于根据获取的基础数据计算母线负荷的半不变量和新能源机组出力的半不变量;
第二计算模块,用于根据母线负荷的半不变量和新能源机组出力的半不变量求取考虑注入功率波动、不平衡功率分配的电压和支路潮流的概率分布;
所述基础数据包括:机组出力及机组出力概率分布描述、母线负荷期望及母线负荷概率分布描述。
所述第一计算模块计算母线负荷的半不变量包括:
第一计算子模块,用于计算母线负荷的各阶原点矩;
第二计算子模块,用于根据所述各阶原点矩,计算母线负荷的各阶半不变量。
所述第一计算模块计算新能源机组出力的半不变量包括:
第三计算子模块,用于计算新能源机组的各阶原点矩;
第四计算子模块,用于根据所述原点矩计算新能源机组出力的各阶半不变量。
所述第二计算模块包括:线性化子模块和计算子模块;
所述线性化子模块基于准稳态灵敏度线性化潮流方程;
所述计算子模块根据半不变量和Gram-Charlier级数展开理论,由注入功率的各阶半不变量计算节点电压和支路潮流的概率分布。
所述线性化子模块,基于准稳态灵敏度线性化潮流方程的表达式如下所示:
式中,X为节点电压向量,Z为支路潮流向量,X0为节点电压向量的基准运行状态,Z0为支路潮流向量基准运行状态,ΔW0为注入功率的变化量,和分别为节点电压和支路潮流对注入功率变化的准稳态灵敏度,Fu为准稳态物理响应线性化后的系数矩阵,S0与T0分别为节点电压和支路潮流对注入功率变化的常规灵敏度。
所述准稳态物理响应线性化后的系数矩阵Fu如下所示:
其中,αn是n维列向量,n为***节点数,αi是第i节点的功率不平衡量分配系数,若αi=0,表示i节点不参与不平衡功率分担。
所述计算子模块根据半不变量和Gram-Charlier级数展开理论,由注入功率的各阶半不变量计算节点电压和支路潮流的概率分布包括:
基于准稳态灵敏度线性化潮流方程,按照下式由注入功率各阶半不变量计算节点电压和支路潮流的各阶半不变量:
式中ΔWk是负荷各阶半不变量与新能源机组出力各阶半不变量相加得到的注入功率各阶半不变量,分别是中元素的k次幂矩阵;
节点电压和支路潮流随机变量规格化按下式表示:
其中,ξ表示节点电压和支路潮流随机变量,μ表示其均值,即一阶半不变量,σ为标准差,即二阶半不变量开平方;
节点电压和支路潮流的概率密度函数f(x)和概率分布函数F(x)按下式表示:
式中,和φ(x)分别是期望值为0,标准差为1的标准正态分布随机变量的概率密度函数和概率分布函数,cr为常系数,计算公式如下所示:
本发明还提供一种考虑不平衡功率分配的概率潮流分析方法,包括:
获取基础数据;
根据获取的基础数据计算母线负荷的半不变量;
根据获取的基础数据计算新能源机组出力的半不变量;
根据母线负荷的半不变量和新能源机组出力的半不变量求取考虑注入功率波动、不平衡功率分配的电压和支路潮流的概率分布;
根据所述节点电压和支路潮流的概率分布,确定是否存在潮流越限风险,若存在,则调整机组出力,否则输出节点电压和支路潮流的概率分布。
本发明提供一种考虑不平衡功率分配的概率潮流分析***,包括:
计算模块,用于根据考虑不平衡功率分配的概率潮流计算方法求取节点电压和支路潮流的概率分布;
判断模块,用于根据所述节点电压和支路潮流的概率分布,确定是否存在潮流越限风险,若存在,则调整机组出力,否则输出节点电压和支路潮流的概率分布。
实施例:
参照附图2:考虑不平衡功率分配的概率潮流计算方法的主要计算步骤及分析步骤流程图。图中,各个模块的工作功能如下:
基础数据准备,包括电网模型、设备参数、机组出力及其概率分布描述、母线负荷期望及其概率分布描述等。
半不变量法概率潮流计算
1)由母线负荷的概率分布特性,得到其原点矩。
对于连续分布的随机变量,设随机变量X的密度函数为f(x),则其r阶原点矩αr可由式(1)求得:
对于离散分布的随机变量,设离散随机变量X取值xi的概率为pi,则r阶原点矩由式(2)定义:
2)根据原点矩与半不变量的关系,计算母线负荷的各阶半不变量。
半不变量与原点矩的关系如下所示:
k1=α1 (3)
α1=k1 (5)
其中,r为阶数,αr为r阶原点矩,kr为r阶半不变量,为二项式系数。
3)计算新能源机组出力的半不变量,本发明以Weibull分布描述的风电为例:
风电场作为间歇性电源的典型代表之一,其在稳态模型研究方面已较为成熟,建立起以风速为变量、输出功率为因变量的近线性数学描述,可以基于此建立风电场出力的概率随机模型。本发明采用三参数Weibull模型,将模型的位置参数设定为风场所在地最小风速,则风速概率密度有:
式中:vw为风速,v0为位置参数,b(b>0)为尺度参数,反映风电场的平均风速;c(c>0)为形状参数。参数b,c可由平均风速μv和标准差σv计算的得到,其中Γ为Gamma函数:
根据机组功率特性和风速概率建立风力发电机组功率随机模型。一般可以认为在切入风速vci至额定风速vr区间成线性关系,则有:
式中Pw、Pr分别为机组实际输出有功功率和额定有功功率;k2=-k1vci。
由风电机组有功功率特性式(10)及式(7)所给出风速概率密度函数就可以构建出风电有功输出概率分布函数,由随机变量的函数分布可求出风力发电有功功率概率密度。通常只需考虑风速在切入风速和额定风速之间的情况,此时风电机组输出功率的概率分布函数和概率密度函数分别描述如下:
式中,a=k1v0+k2,β=k1b。
当已知随机变量的分布时,可求出它的各阶矩。由式(11)可得风电机组输出功率的r阶原点矩αr:
最后,根据式(3)-(6),原点矩与半不变量的关系,计算风电机组出力的各阶半不变量。
线性化交流潮流方程
极坐标形式下的***潮流方程可以表示为:
其中,W为节点注入功率向量,f为功率方程,x为节点电压构成的状态向量,Z为支路潮流向量,g为支路潮流方程。
将式(13)在基准运行点处利用泰勒级数展开,忽略二次及以上各阶项,可得:
式中,X为节点电压向量,下标0表示基准运行状态;ΔW为注入功率的变化量;S0与T0分别为节点电压和支路潮流对注入功率变化的灵敏度。
考虑注入功率随机波动引起的不平衡功率分配的概率潮流计算
(1)准稳态灵敏度计算
传统的解析法概率潮流计算中,负荷波动、新能源机组出力波动等产生的不平衡功率由平衡节点承担,并不符合电网的实际情况。本发明以准稳态灵敏度代替常规灵敏度,构建支路潮流与节点注入功率的线性关系,可以灵活的处理不平衡功率在各节点的分摊方式。
考虑电力***的准稳态物理响应的灵敏度分析方法称为准稳态灵敏度分析。这里“准稳态”的物理响应是指***在经受操作或扰动后,不计***暂态过程,但计及***扰动前后新旧稳态间的总变化。
由于电力***的物理响应,***在ΔW0的作用下,经过一段过渡过程后将达到一个新稳态。假设在新稳态下,真正作用在电力网络上的控制量为ΔW,则增量形式的准稳态物理响应可表示为:
ΔW=FuΔW0 (15)
其中,Fu为准稳态物理响应线性化后的系数矩阵。
式(15)结合式(14),可得:
准稳态灵敏度系数为:
不同的电力***有不同的准稳态物理响应,系数矩阵Fu可能千差万别。对于本发明中涉及的有功控制,考虑发电机调速器和负荷频率特性,当电力***达到一个新稳态时,***功率不平衡量将由所有的发电机和负荷来分担。设***的功率不平衡量按照分担系数αn在各节点母线之间进行分摊,则
其中,αn是n维列向量,n为***节点数,αi是第i节点的分配系数。若αi=0,表示i节点不参与不平衡功率分担。则ΔW由式(19)所示。
其中,Ιn是n维全1列向量。即,
采用半不变量与Gram-Charlie级数展开相结合的方法,求取支路潮流的概率分布。
①半不变量的性质:
性质1:若有n个相互独立的随机变量x(1),x(2),…,x(n),且各自有r阶半不变量存在,随机变量x为这n个随机变量之和,则x的r阶半不变量等于各独立随机变量的r阶半不变量之和,即:
性质2:随机变量y是随机变量x的线性函数,y=ax+b,为x的r阶半不变量,则随机变量y的各阶半不变量为:
②Gram-Charlier级数
对任一随机变量ξ,假定其均值为μ,标准差为σ。则其规格化随机变量ξ为
由Gram-Charlier级数展开理论,规格化随机变量的概率密度函数f(x)和概率分布函数F(x)分别为:
式中和φ(x)分别是期望值为0,标准差为1的标准正态分布随机变量的概率密度函数和概率分布函数,cr为常系数:
c0=1
c1=c2=0
…
③概率潮流分析
在式(16)所示的交流潮流模型中,采用准稳态灵敏度来描述节点功率变化量与节点电压与支路潮流变化量之间的关系,并由以下步骤求解考虑不平衡功率分配的概率潮流,以合理分配节点注入功率波动引起的功率不平衡量。
a.首先由母线负荷的概率分布特性,得到其原点矩。
b.根据原点矩与半不变量的关系,计算母线负荷的各阶半不变量。
c.由新能源机组的概率分布特性,得到其原点矩。
d.根据原点矩与半不变量的关系,计算新能源机组出力的各阶半不变量。
e.基于半不变量与Gram-Charlie级数展开相结合的方法,采用准稳态灵敏度来描述节点功率变化量与节点电压与支路潮流变化量之间的关系,求取节点电压和支路潮流的概率分布。
根据求取的节点电压和支路潮流的概率分布,判据是否存在电网潮流越限的安全风险,如果有,则可进一步调整机组出力(计划)以消除越限风险。
结果展示和输出,包括线路潮流及断面潮流的概率分布,越限概率等。
本领域内的技术人员应明白,本申请的实施例可提供为方法、***、或计算机程序产品。因此,本申请可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且,本申请可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本申请是参照根据本申请实施例的方法、设备(***)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器,使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品,该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上,使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理,从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
最后应当说明的是:以上实施例仅用于说明本发明的技术方案而非对其保护范围的限制,尽管参照上述实施例对本申请进行了详细的说明,所属领域的普通技术人员应当理解:本领域技术人员阅读本申请后依然可对申请的具体实施方式进行种种变更、修改或者等同替换,但这些变更、修改或者等同替换,均在申请待批的权利要求保护范围之内。
Claims (12)
1.一种考虑不平衡功率分配的概率潮流计算方法,其特征在于,包括:
获取基础数据;
根据获取的基础数据计算母线负荷的半不变量和新能源机组出力的半不变量;
根据母线负荷的半不变量和新能源机组出力的半不变量求取考虑注入功率波动、不平衡功率分配的电压和支路潮流的概率分布;
所述基础数据包括:机组出力及机组出力概率分布描述、母线负荷期望及母线负荷概率分布描述。
2.如权利要求1所述的考虑不平衡功率分配的概率潮流计算方法,其特征在于,所述计算母线负荷的半不变量包括:
计算母线负荷的各阶原点矩;
根据所述各阶原点矩,计算母线负荷的各阶半不变量。
3.如权利要求1所述的考虑不平衡功率分配的概率潮流计算方法,其特征在于,所述计算新能源机组出力的半不变量包括:
计算新能源机组的各阶原点矩;
根据所述原点矩计算新能源机组出力的各阶半不变量。
4.如权利要求1所述的考虑不平衡功率分配的概率潮流计算方法,其特征在于,所述求取考虑注入功率波动、不平衡功率分配的电压和支路潮流的概率分布包括:
基于准稳态灵敏度线性化潮流方程;
根据半不变量及Gram-Charlier级数展开理论,由注入功率的各阶半不变量计算节点电压和支路潮流的概率分布。
5.如权利要求4所述的考虑不平衡功率分配的概率潮流计算方法,其特征在于,所述基于准稳态灵敏度线性化潮流方程的表达式如下所示:
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>X</mi>
<mo>=</mo>
<msub>
<mi>X</mi>
<mn>0</mn>
</msub>
<mo>+</mo>
<mi>&Delta;</mi>
<mi>X</mi>
<mo>=</mo>
<msub>
<mi>X</mi>
<mn>0</mn>
</msub>
<mo>+</mo>
<msubsup>
<mi>S</mi>
<mn>0</mn>
<mi>R</mi>
</msubsup>
<msup>
<mi>&Delta;W</mi>
<mn>0</mn>
</msup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>Z</mi>
<mo>=</mo>
<msub>
<mi>Z</mi>
<mn>0</mn>
</msub>
<mo>+</mo>
<mi>&Delta;</mi>
<mi>Z</mi>
<mo>=</mo>
<msub>
<mi>Z</mi>
<mn>0</mn>
</msub>
<mo>+</mo>
<msubsup>
<mi>T</mi>
<mn>0</mn>
<mi>R</mi>
</msubsup>
<msup>
<mi>&Delta;W</mi>
<mn>0</mn>
</msup>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<msubsup>
<mi>S</mi>
<mn>0</mn>
<mi>R</mi>
</msubsup>
<mo>=</mo>
<msub>
<mi>S</mi>
<mn>0</mn>
</msub>
<msub>
<mi>F</mi>
<mi>u</mi>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msubsup>
<mi>T</mi>
<mn>0</mn>
<mi>R</mi>
</msubsup>
<mo>=</mo>
<msub>
<mi>T</mi>
<mn>0</mn>
</msub>
<msub>
<mi>F</mi>
<mi>u</mi>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
式中,X为节点电压向量,Z为支路潮流向量,X0为节点电压向量的基准运行状态,Z0为支路潮流向量基准运行状态,ΔW0为注入功率的变化量,和分别为节点电压和支路潮流对注入功率变化的准稳态灵敏度,Fu为准稳态物理响应线性化后的系数矩阵,S0与T0分别为节点电压和支路潮流对注入功率变化的常规灵敏度。
6.如权利要求5所述的考虑不平衡功率分配的概率潮流计算方法,其特征在于,所述准稳态物理响应线性化后的系数矩阵Fu如下所示:
<mrow>
<munder>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>&Element;</mo>
<mi>n</mi>
</mrow>
</munder>
<msub>
<mi>&alpha;</mi>
<mi>i</mi>
</msub>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<msub>
<mi>&alpha;</mi>
<mi>i</mi>
</msub>
<mo>&GreaterEqual;</mo>
<mn>0</mn>
</mrow>
其中,αn是n维列向量,n为***节点数,αi是第i节点的功率不平衡量分配系数,若αi=0,表示i节点不参与不平衡功率分配。
7.如权利要求4所述的考虑不平衡功率分配的概率潮流计算方法,其特征在于,所述根据半不变量和Gram-Charlier级数展开理论,由注入功率的各阶半不变量计算节点电压和支路潮流的概率分布包括:
基于准稳态灵敏度线性化潮流方程,按照下式由注入功率各阶半不变量计算节点电压和支路潮流的各阶半不变量:
<mrow>
<msub>
<mi>&Delta;X</mi>
<mi>k</mi>
</msub>
<mo>=</mo>
<msubsup>
<mi>S</mi>
<mn>0</mn>
<mrow>
<mi>R</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
</msubsup>
<msub>
<mi>&Delta;W</mi>
<mi>k</mi>
</msub>
</mrow>
<mrow>
<msub>
<mi>&Delta;Z</mi>
<mi>k</mi>
</msub>
<mo>=</mo>
<msubsup>
<mi>T</mi>
<mn>0</mn>
<mrow>
<mi>R</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
</msubsup>
<msub>
<mi>&Delta;W</mi>
<mi>k</mi>
</msub>
</mrow>
式中,ΔWk是负荷各阶半不变量与新能源机组出力各阶半不变量相加得到的注入功率各阶半不变量,分别是中元素的k次幂矩阵;
节点电压和支路潮流随机变量规格化按下式表示:
<mrow>
<mover>
<mi>&xi;</mi>
<mo>&OverBar;</mo>
</mover>
<mo>=</mo>
<mfrac>
<mrow>
<mi>&xi;</mi>
<mo>-</mo>
<mi>&mu;</mi>
</mrow>
<mi>&sigma;</mi>
</mfrac>
</mrow>
其中,ξ表示节点电压和支路潮流随机变量,μ表示其均值,即一阶半不变量,σ为标准差,即二阶半不变量开平方;
节点电压和支路潮流的概率密度函数f(x)和概率分布函数F(x)按下式表示:
式中,和φ(x)分别是期望值为0,标准差为1的标准正态分布随机变量的概率密度函数和概率分布函数,cr为常系数,计算公式如下所示:
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>c</mi>
<mn>0</mn>
</msub>
<mo>=</mo>
<mn>1</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>c</mi>
<mn>1</mn>
</msub>
<mo>=</mo>
<msub>
<mi>c</mi>
<mn>2</mn>
</msub>
<mo>=</mo>
<mn>0</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>c</mi>
<mn>3</mn>
</msub>
<mo>=</mo>
<mo>-</mo>
<mfrac>
<mrow>
<msub>
<mi>&alpha;</mi>
<mn>3</mn>
</msub>
<mo>-</mo>
<mn>3</mn>
<msub>
<mi>&alpha;</mi>
<mn>1</mn>
</msub>
<msub>
<mi>&alpha;</mi>
<mn>2</mn>
</msub>
<mo>+</mo>
<mn>2</mn>
<msubsup>
<mi>&alpha;</mi>
<mn>1</mn>
<mn>3</mn>
</msubsup>
</mrow>
<msup>
<mi>&sigma;</mi>
<mn>3</mn>
</msup>
</mfrac>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>c</mi>
<mn>4</mn>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>&alpha;</mi>
<mn>4</mn>
</msub>
<mo>-</mo>
<mn>4</mn>
<msub>
<mi>&alpha;</mi>
<mn>1</mn>
</msub>
<msub>
<mi>&alpha;</mi>
<mn>3</mn>
</msub>
<mo>+</mo>
<mn>6</mn>
<msubsup>
<mi>&alpha;</mi>
<mn>1</mn>
<mn>2</mn>
</msubsup>
<msub>
<mi>&alpha;</mi>
<mn>2</mn>
</msub>
<mo>-</mo>
<mn>3</mn>
<msubsup>
<mi>&alpha;</mi>
<mn>1</mn>
<mn>4</mn>
</msubsup>
</mrow>
<msup>
<mi>&sigma;</mi>
<mn>4</mn>
</msup>
</mfrac>
<mo>-</mo>
<mn>3</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>......</mn>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>.</mo>
</mrow>
8.如权利要求1所述的考虑不平衡功率分配的概率潮流计算方法,其特征在于,所述基础数据还包括电网模型和设备参数。
9.一种考虑不平衡功率分配的概率潮流计算***,其特征在于,包括:
采集模块,用于获取基础数据;
第一计算模块,用于根据获取的基础数据计算母线负荷的半不变量和新能源机组出力的半不变量;
第二计算模块,用于根据母线负荷的半不变量和新能源机组出力的半不变量求取考虑注入功率波动、不平衡功率分配的电压和支路潮流的概率分布;
所述基础数据包括:机组出力及机组出力概率分布描述、母线负荷期望及母线负荷概率分布描述。
10.如权利要求9所述的考虑不平衡功率分配的概率潮流计算***,其特征在于,所述第三计算模块包括:线性化子模块和计算子模块;
所述线性化子模块基于准稳态灵敏度线性化潮流方程;
所述计算子模块根据半不变量和Gram-Charlier级数展开理论,由注入功率的各阶半不变量计算节点电压和支路潮流的概率分布。
11.一种考虑不平衡功率分配的概率潮流分析方法,其特征在于,包括:
根据权利要求1所述的考虑不平衡功率分配的概率潮流计算方法求取节点电压和支路潮流的概率分布;
根据所述节点电压和支路潮流的概率分布,确定是否存在潮流越限风险,若存在,则调整机组出力,否则输出节点电压和支路潮流的概率分布。
12.一种考虑不平衡功率分配的概率潮流分析***,其特征在于,包括:
计算模块,用于根据权利要求1所述的考虑不平衡功率分配的概率潮流计算方法求取节点电压和支路潮流的概率分布;
判断模块,用于根据所述节点电压和支路潮流的概率分布,确定是否存在潮流越限风险,若存在,则调整机组出力,否则输出节点电压和支路潮流的概率分布。
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711092861.7A CN107968409A (zh) | 2017-11-08 | 2017-11-08 | 一种考虑不平衡功率分配的概率潮流计算方法及*** |
PCT/CN2018/114029 WO2019091356A1 (zh) | 2017-11-08 | 2018-11-05 | 一种概率潮流计算方法及*** |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711092861.7A CN107968409A (zh) | 2017-11-08 | 2017-11-08 | 一种考虑不平衡功率分配的概率潮流计算方法及*** |
Publications (1)
Publication Number | Publication Date |
---|---|
CN107968409A true CN107968409A (zh) | 2018-04-27 |
Family
ID=62000870
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711092861.7A Pending CN107968409A (zh) | 2017-11-08 | 2017-11-08 | 一种考虑不平衡功率分配的概率潮流计算方法及*** |
Country Status (2)
Country | Link |
---|---|
CN (1) | CN107968409A (zh) |
WO (1) | WO2019091356A1 (zh) |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109390946A (zh) * | 2018-10-08 | 2019-02-26 | 重庆大学 | 一种基于多参数规划理论的最优概率潮流快速计算方法 |
WO2019091356A1 (zh) * | 2017-11-08 | 2019-05-16 | 国家电网有限公司 | 一种概率潮流计算方法及*** |
CN109932587A (zh) * | 2018-12-04 | 2019-06-25 | 国网浙江省电力有限公司衢州供电公司 | 一种用户侧无功补偿设备均衡性监测方法 |
CN110048428A (zh) * | 2019-05-27 | 2019-07-23 | 河北工业大学 | 基于概率守恒原理的电力***概率潮流计算方法 |
CN110535136A (zh) * | 2019-08-16 | 2019-12-03 | 广东工业大学 | 一种交流潮流线性灵敏度分析方法 |
CN112069666A (zh) * | 2020-08-25 | 2020-12-11 | 贵州大学 | 一种基于概率潮流法的电网短期可靠性评估方法 |
CN113507128A (zh) * | 2021-07-14 | 2021-10-15 | 国网青海省电力公司 | 一种特高压直流换流站近场无功优化配置方法 |
CN113723821A (zh) * | 2021-08-31 | 2021-11-30 | 广东电网有限责任公司 | 一种基于潮流介数的电网故障预警方法及装置 |
CN113972659A (zh) * | 2021-11-04 | 2022-01-25 | 国网冀北电力有限公司经济技术研究院 | 一种考虑随机潮流的储能配置方法及*** |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114268102B (zh) * | 2021-12-23 | 2024-03-22 | 国网江苏省电力有限公司经济技术研究院 | 基于解析式概率潮流模型的电力***运行状态量化方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103995982A (zh) * | 2014-06-06 | 2014-08-20 | 国家电网公司 | 一种考虑机组随机故障的概率潮流计算方法 |
CN106099939A (zh) * | 2016-08-08 | 2016-11-09 | 国家电网公司 | 一种变电站无功设备对母线电压影响灵敏度的计算方法 |
CN106786595A (zh) * | 2016-11-29 | 2017-05-31 | 国电南瑞科技股份有限公司 | 一种考虑静态频率特性不确定性的概率潮流计算方法 |
CN107069721A (zh) * | 2017-06-21 | 2017-08-18 | 华北电力大学 | 一种基于随机集理论的电力***运行风险评估方法 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP4203993B2 (ja) * | 2003-08-01 | 2009-01-07 | 株式会社日立製作所 | 電力系統の解析支援方法 |
CN104410069B (zh) * | 2014-12-05 | 2017-02-08 | 国家电网公司 | 一种计及响应相关性的动态概率潮流计算方法 |
CN105811403B (zh) * | 2016-03-22 | 2018-05-22 | 江苏省电力试验研究院有限公司 | 基于半不变量和级数展开法的随机潮流算法 |
CN107968409A (zh) * | 2017-11-08 | 2018-04-27 | 中国电力科学研究院有限公司 | 一种考虑不平衡功率分配的概率潮流计算方法及*** |
-
2017
- 2017-11-08 CN CN201711092861.7A patent/CN107968409A/zh active Pending
-
2018
- 2018-11-05 WO PCT/CN2018/114029 patent/WO2019091356A1/zh active Application Filing
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103995982A (zh) * | 2014-06-06 | 2014-08-20 | 国家电网公司 | 一种考虑机组随机故障的概率潮流计算方法 |
CN106099939A (zh) * | 2016-08-08 | 2016-11-09 | 国家电网公司 | 一种变电站无功设备对母线电压影响灵敏度的计算方法 |
CN106786595A (zh) * | 2016-11-29 | 2017-05-31 | 国电南瑞科技股份有限公司 | 一种考虑静态频率特性不确定性的概率潮流计算方法 |
CN107069721A (zh) * | 2017-06-21 | 2017-08-18 | 华北电力大学 | 一种基于随机集理论的电力***运行风险评估方法 |
Non-Patent Citations (2)
Title |
---|
卫鹏等: "基于半不变量和Gram-Charlier级数展开法的随机潮流算法", 《电力工程技术》 * |
黄国栋等: "安全约束经济调度中有功潮流调整方法", 《电力***保护与控制》 * |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2019091356A1 (zh) * | 2017-11-08 | 2019-05-16 | 国家电网有限公司 | 一种概率潮流计算方法及*** |
CN109390946A (zh) * | 2018-10-08 | 2019-02-26 | 重庆大学 | 一种基于多参数规划理论的最优概率潮流快速计算方法 |
CN109932587B (zh) * | 2018-12-04 | 2021-01-08 | 国网浙江省电力有限公司衢州供电公司 | 一种用户侧无功补偿设备均衡性监测方法 |
CN109932587A (zh) * | 2018-12-04 | 2019-06-25 | 国网浙江省电力有限公司衢州供电公司 | 一种用户侧无功补偿设备均衡性监测方法 |
CN110048428A (zh) * | 2019-05-27 | 2019-07-23 | 河北工业大学 | 基于概率守恒原理的电力***概率潮流计算方法 |
CN110535136A (zh) * | 2019-08-16 | 2019-12-03 | 广东工业大学 | 一种交流潮流线性灵敏度分析方法 |
CN112069666A (zh) * | 2020-08-25 | 2020-12-11 | 贵州大学 | 一种基于概率潮流法的电网短期可靠性评估方法 |
CN112069666B (zh) * | 2020-08-25 | 2022-11-29 | 贵州大学 | 一种基于概率潮流法的电网短期可靠性评估方法 |
CN113507128A (zh) * | 2021-07-14 | 2021-10-15 | 国网青海省电力公司 | 一种特高压直流换流站近场无功优化配置方法 |
CN113723821A (zh) * | 2021-08-31 | 2021-11-30 | 广东电网有限责任公司 | 一种基于潮流介数的电网故障预警方法及装置 |
CN113723821B (zh) * | 2021-08-31 | 2024-04-26 | 广东电网有限责任公司 | 一种基于潮流介数的电网故障预警方法及装置 |
CN113972659A (zh) * | 2021-11-04 | 2022-01-25 | 国网冀北电力有限公司经济技术研究院 | 一种考虑随机潮流的储能配置方法及*** |
CN113972659B (zh) * | 2021-11-04 | 2024-02-02 | 国网冀北电力有限公司经济技术研究院 | 一种考虑随机潮流的储能配置方法及*** |
Also Published As
Publication number | Publication date |
---|---|
WO2019091356A1 (zh) | 2019-05-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107968409A (zh) | 一种考虑不平衡功率分配的概率潮流计算方法及*** | |
Xu et al. | Robust dispatch of high wind power-penetrated power systems against transient instability | |
Dong et al. | Probabilistic load flow analysis for power system containing wind farms | |
US20210367424A1 (en) | Multi-Objective Real-time Power Flow Control Method Using Soft Actor-Critic | |
CN103825269B (zh) | 一种考虑电力***功频静特性的快速概率潮流计算方法 | |
CN105811403A (zh) | 基于半不变量和级数展开法的随机潮流算法 | |
CN110429648B (zh) | 考虑风速随机波动的小干扰稳定裕度概率评估方法 | |
CN104156542B (zh) | 一种基于隐式投影的有源配电***稳定性仿真方法 | |
CN105634018B (zh) | 一种基于随机响应面法和内点法的含风电随机最优潮流求解方法 | |
CN115549187A (zh) | 一种风电并网的稳定性评估方法 | |
Gurung et al. | Probabilistic small-signal stability analysis of power system with solar farm integration | |
Hu et al. | Frequency prediction model combining ISFR model and LSTM network | |
Zadeh et al. | A novel probabilistic method for under frequency load shedding setting considering wind turbine response | |
Zhao et al. | Probabilistic voltage stability assessment considering stochastic load growth direction and renewable energy generation | |
Zaker et al. | PMU-based linear and nonlinear black-box modelling of power systems | |
CN106855909B (zh) | 一种适于有源配电网随机动态仿真的显隐混合积分方法 | |
Li et al. | Identifying critical nodes in power grids containing renewable energy based on electrical spreading probability | |
CN110048428B (zh) | 基于概率守恒原理的电力***概率潮流计算方法 | |
CN105808825B (zh) | 一种适于有源配电网动态仿真的滚动式随机投影积分方法 | |
Zhu et al. | Efficient identification of critical load model parameters affecting power system voltage stability | |
Chenxu et al. | An improved cumulant method for probabilistic load flow calculation | |
Philippe et al. | A copula-based uncertainty modeling of wind power generation for probabilistic power flow study | |
CN107947179A (zh) | 一种考虑风电调度策略的多点线性化随机潮流计算方法 | |
Bian et al. | Identification and improvement of probabilistic voltage instability modes of power system with wind power integration | |
Tajalli et al. | Maximizing social welfare considering the uncertainty of wind power plants using a distributed consensus-based algorithm |
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 | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20180427 |
|
RJ01 | Rejection of invention patent application after publication |