CN110048428B - 基于概率守恒原理的电力***概率潮流计算方法 - Google Patents
基于概率守恒原理的电力***概率潮流计算方法 Download PDFInfo
- Publication number
- CN110048428B CN110048428B CN201910444918.8A CN201910444918A CN110048428B CN 110048428 B CN110048428 B CN 110048428B CN 201910444918 A CN201910444918 A CN 201910444918A CN 110048428 B CN110048428 B CN 110048428B
- Authority
- CN
- China
- Prior art keywords
- power
- probability
- probability density
- power system
- branch
- 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.)
- Expired - Fee Related
Links
Images
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
-
- 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
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明为基于概率守恒原理的电力***概率潮流计算方法,该方法包括:获取含有新能源的电力***网络的各项数据;对常规发电机出力、新能源出力和负荷功率历史数据进行统计分析,并生成电力***所有节点注入功率的概率模型;建立电力***潮流计算方程的线性化模型;按照线性化方程的关系,基于概率论中的概率守恒原理的原理建立不同随机变量的概率密度函数之间的关系式;利用δ函数的定义和性质,对不同随机变量概率密度函数之间的关系式进行处理生成电力***内节点电压、支路潮流和支路有功损耗的概率密度函数。该方法能够得到准确的概率潮流计算的结果,合理反映了***内节点电压、支路潮流和支路有功损耗的概率特性,且计算效率高。
Description
技术领域
本发明是一种基于概率守恒原理的电力***概率潮流计算方法,以概率论中概率守恒原理为理论基础,对电力***中节点电压、节点注入功率、支路潮流和有功损耗之间概率密度函数之间的关系进行推导,并得到概率潮流计算的结果,属于电力***概率潮流计算研究领域。
背景技术
电力***作为典型的多变量、高维、复杂、非线性***,不可避免地存在着不确定性和风险,伴随新能源电力渗透率不断提高,电力***运行的不确定性更为突出。由于新能源出力的波动性、间歇性与随机性,其发电出力的不确定性使得其接入配电网后对电网的无功潮流分布、有功潮流分布、电压稳定性等也会产生很大的冲击,这直接影响到了配电网的安全稳定性,关系到用户的用电效率和用电质量。而且,由于分布式电源接入配电网的位置和功率大小的不同,会引起整条线路上面潮流流向的改变,出现功率倒送的情况,对传统的配电网线路和保护***产生很大的影响。
潮流是电力***稳态运行分析的基本计算内容,是电力***规划、调度运行相关各项分析的基础。在含有新能源发电的电力***运行中,由于自然条件导致的新能源出力的随机性和波动性是一个导致电力***不确定性的主要因素。此时,常规的确定性潮流计算已不能准确地分析随机因素对***运行的影响,不确定性潮流分析方法应运而生,其中概率潮流计算是不确定性潮流计算的主要方法之一。
概率潮流在1974年由Borkwska提出,已发展出模拟法、近似法和解析法。模拟法以蒙特卡洛法(Monter Carlo,MC)为代表,其精确性需要大量统计取样及统计模拟,导致计算量大、时间长。近似法以点估计法为代表,根据已知输入随机变量的概率分布,求输出随机变量各阶矩的概率统计,计算速度较快,但输出随机变量的高阶矩误差较大,***分析中需要兼顾精确性与快速性。解析法以半不变量法为代表,通过求解随机变量的半不变量矩阵,进行矩阵运算,最后通过级数展开式得到潮流响应的概率密度函数或概率分布函数和统计特征,但级数展开式容易出现不收敛的情况,导致概率密度函数出现负值的情况。
发明内容
针对现有技术的不足,本发明拟解决的技术问题是,提供了一种基于概率守恒原理的电力***概率潮流计算方法。该方法能够得到准确的含有新能源发电的电力***的概率潮流计算的结果,合理反映了***内节点电压、支路潮流和支路有功损耗的概率特性,且计算效率高,为电网人员对电力***不确定分析的工作奠定基础。
本发明解决所述技术问题采用的技术方案是:
一种基于概率守恒原理的电力***概率潮流计算方法,所述方法包括以下步骤:
步骤1.获取含有新能源的电力***网络的各项数据;
步骤2.依据步骤1中的数据建立用于分析电力***潮流的电力***网络模型,形成n个节点,m个支路,并对步骤1中常规发电机数据、新能源出力和负荷功率历史数据进行统计分析,生成电力***每个节点注入功率的联合概率密度函数,该联合概率密度函数包括常规发电机出力的概率密度函数、新能源出力的概率密度函数和负荷功率的概率密度函数,然后再对所有节点上的注入功率的联合概率密度函数进行累乘,得到整个***所有节点注入功率的概率密度函数pY(y);
步骤3.以节点注入功率Y、支路潮流Z、支路有功损耗U、节点电压X作为随机变量,建立电力***潮流计算方程的线性化模型;
步骤4.按照步骤3所述的线性化模型,基于概率论中的概率守恒原理的原理,建立不同随机变量的概率密度函数之间的关系式;
步骤5.利用δ函数的定义和性质,对步骤4所述的不同随机变量的概率密度函数之间的关系式进行处理,用节点注入功率表示支路潮流Z、支路有功损耗U、节点电压X的概率密度函数,即生成电力***内节点电压、支路潮流和支路有功损耗的概率密度函数。
步骤6.输出概率潮流计算结果。
与现有技术相比,本发明的有益效果在于:
(1)本发明该方法包括如下步骤:获取含有新能源的电力***网络的各项数据;对常规发电机出力、新能源出力和负荷功率历史数据进行统计分析,并生成电力***所有节点注入功率的概率模型;建立电力***潮流计算方程的线性化模型;按照线性化方程的关系,基于概率论中的概率守恒原理的原理建立不同随机变量的概率密度函数之间的关系式;利用δ函数的定义和性质,对不同随机变量概率密度函数之间的关系式进行处理生成电力***内节点电压、支路潮流和支路有功损耗的概率密度函数;输出概率潮流计算结果。
通过电力***潮流计算方程的线性化模型,基于概率论中的概率守恒原理,推导出不同随机变量概率密度函数之间关系的推导,进而引入δ函数后,直接生成了潮流响应(即节点电压、支路潮流、支路有功损耗)的概率密度函数,不会产生上述的半不变量法级数展开过程中不收敛导致的概率密度函数存在负值的情况。
(2)本发明方法不需要运用类似蒙特卡洛法这样的随机抽样方式来对潮流计算的过程进行大量的循环模拟和统计,计算时间短,效率高。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明基于概率守恒原理的电力***概率潮流计算方法的流程图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
如图1所示为本发明一种基于概率守恒原理的电力***概率潮流计算方法的流程图,其步骤包括:
步骤1.获取含有新能源的电力***网络的各项数据;
步骤2.依据步骤1中的数据建立用于分析电力***潮流的电力***网络模型,形成n个节点,m个支路,并对步骤1中常规(火电、水电)发电机数据、新能源出力和负荷功率历史数据进行统计分析,生成电力***每个节点注入功率的联合概率密度函数,该联合概率密度函数包括发电机出力的概率密度函数、新能源出力的概率密度函数和负荷功率的概率密度函数,然后再对所有节点上的注入功率的联合概率密度函数进行累乘,得到整个***所有节点注入功率的概率密度函数pY(y);
步骤3.以节点注入功率Y、支路潮流Z、支路有功损耗U、节点电压X作为随机变量,建立电力***潮流计算方程的线性化模型;
步骤4.按照步骤3所述的线性化模型,基于概率论中的概率守恒原理的原理,建立不同随机变量的概率密度函数之间的关系式;
步骤5.利用δ函数的定义和性质,对步骤4所述的不同随机变量的概率密度函数之间的关系式进行处理,用节点注入功率表示支路潮流Z、支路有功损耗U、节点电压X的概率密度函数,即生成电力***内节点电压、支路潮流和支路有功损耗的概率密度函数。
步骤6.输出概率潮流计算结果。
进一步的,所述的步骤1中,需要获取的含有新能源的电力***网络的各项数据包括:电力***网络结构、线路电压等级、线路阻抗值、常规发电机的功率因数和故障率、新能源历史出力、负荷数据等。
进一步的,所述的步骤2包括以下步骤:
步骤201:按照步骤1中获取的含有新能源的电力***网络的各项数据建立电力***网络模型;
步骤202:规定常规发电机仅有即正常运行与发生故障两种运行状态,单台发电机的故障率为f,常规发电机总台数为t,有i台发电机发生故障的概率pi可以用二项分布来描述,如式(1):
式中,pi表示有i台发电机发生故障的概率;i为发生故障的发电机台数;C表示组合运算符号。
式中,yG为节点上常规发电机的注入功率即发电机的出力;yGi为有i台发电机发生故障的发电机实际出力,δ(·)为狄拉克函数。
步骤203:对一定时间周期内新能源出力的历史数据进行统计,采用Weibull分布和Beta分布分别拟合风电出力和光伏出力的概率密度函数,构成新能源出力概率密度函数;
步骤204:随着电力***实际运行,测量记录的节点上的负荷功率,对一定时间周期内的负荷数据进行统计,计算负荷数据的均值和方差,采用高斯分布拟合负荷功率的概率密度函数;
步骤205:将对应的节点上的常规发电机出力概率密度函数、新能源出力概率密度函数和负荷功率的概率密度函数相乘,生成某个节点注入功率的联合概率密度函数。
再对所有节点上得注入功率的联合概率密度函数进行累乘,得到整个***所有节点注入功率的概率密度函数pY(y)
式中,n为电力***节点个数。
进一步的,所述的步骤3包括:
以节点注入功率Y、支路潮流Z、支路有功损耗U、节点电压X作为随机变量,电力***潮流计算方程简化为:
式中,n为电力***节点个数;m为支路条数;X为节点电压(包括节点电压的幅值Vm和相角Va),其中的元素用x表示;Y为节点注入功率(包括节点注入的有功功率P和无功功率Q),其中的元素用y表示;Z为支路潮流(包括支路上流过的有功功率Pl和无功功率Ql),其中的元素用z表示;U为支路有功损耗Ploss,其中的元素用u表示;f、g、h为随机变量之间的映射关系;
对式(3)进行线性化,得到电力***潮流计算方程的线性化模型:
式中,J0、G0、S0分别为节点注入功率、支路潮流和支路有功损耗对节点电压求偏导生成的矩阵。
所述的步骤4包括以下步骤:
步骤401:基于概率论中的概率守恒原理,建立不同随机变量的概率密度函数的关系为式(6):
式中,p(·)为不同随机变量的概率密度函数;Ω为每个随机变量中所有元素组成的区域;Y代表节点注入功率,y是Y其中的元素;X为节点电压,其中的元素用x表示;Z为支路潮流,其中的元素用z表示;U为支路有功损耗Ploss,其中的元素用u表示;pY(y)为所有节点注入功率的概率密度函数;
步骤402:通过式(5)的线性化模型对式(6)进行变换,推导出不同随机变量的概率密度函数的关系如下式:
式中,表示矩阵的行列式运算符号。
所述的步骤5包括:
利用δ函数的定义和性质,对步骤4所述的不同随机变量的概率密度函数之间的关系式(7)进行变换,生成电力***内节点电压、支路潮流和支路有功损耗的概率密度函数:
式中,δ(·)为狄拉克函数;Dy为随机变量节点注入功率Y中每个元素的定义域。
进而,对于第i个节点的节点电压,第j条支路的支路潮流和支路有功损耗的概率密度函数为:
进一步的,所述的步骤6包括:
以表格数据的形式输出步骤5得到的随机变量(支路潮流Z、支路有功损耗U、节点电压X)的概率密度函数,生成对应的曲线,分析电力***节点电压、支路潮流和支路有功损耗的概率分布。
仿真实验表明,在相同实施例下,本发明实验结果仿真时间可以控制在10000次蒙特卡洛(段玉兵,龚宇雷,谭兴国,et al.基于蒙特卡罗模拟的微电网随机潮流计算方法[J].电工技术学报,2011,26(1增).)仿真时间的30%内,误差可以控制在10000次蒙特卡洛仿真结果的0.01%数量级之内。
本发明中应用了具体实施例对本发明的原理及实施方式进行了阐述,以上实施例的说明只是用于帮助理解本发明的方法及其核心思想;同时,对于本领域的一般技术人员,依据本发明的思想,在具体实施方式及应用范围上均会有改变之处,综上所述,本说明书内容不应理解为对本发明的限制。
本发明未述及之处适用于现有技术。
Claims (5)
1.一种基于概率守恒原理的电力***概率潮流计算方法,所述方法包括以下步骤:
步骤1.获取含有新能源的电力***网络的各项数据;
步骤2.依据步骤1中的数据建立用于分析电力***潮流的电力***网络模型,形成n个节点,m个支路,并对步骤1中的电力***网络的常规发电机数据、新能源出力和负荷功率历史数据进行统计分析,生成电力***每个节点注入功率的联合概率密度函数,该联合概率密度函数包括常规发电机出力的概率密度函数、新能源出力的概率密度函数和负荷功率的概率密度函数,然后再对所有节点上的注入功率的联合概率密度函数进行累乘,得到整个***所有节点注入功率的概率密度函数;
步骤3.以节点注入功率、支路潮流、支路有功损耗、节点电压作为随机变量,建立电力***潮流计算方程的线性化模型;
步骤4.按照步骤3所述的线性化模型,基于概率论中的概率守恒原理,建立不同随机变量的概率密度函数之间的关系式;
步骤5.利用δ函数的定义和性质,对步骤4所述的不同随机变量的概率密度函数之间的关系式进行处理,用节点注入功率表示支路潮流、支路有功损耗、节点电压的概率密度函数,即生成电力***内节点电压、支路潮流和支路有功损耗的概率密度函数;
步骤6.输出概率潮流计算结果。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910444918.8A CN110048428B (zh) | 2019-05-27 | 2019-05-27 | 基于概率守恒原理的电力***概率潮流计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910444918.8A CN110048428B (zh) | 2019-05-27 | 2019-05-27 | 基于概率守恒原理的电力***概率潮流计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110048428A CN110048428A (zh) | 2019-07-23 |
CN110048428B true CN110048428B (zh) | 2020-06-30 |
Family
ID=67283687
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910444918.8A Expired - Fee Related CN110048428B (zh) | 2019-05-27 | 2019-05-27 | 基于概率守恒原理的电力***概率潮流计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110048428B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113723821B (zh) * | 2021-08-31 | 2024-04-26 | 广东电网有限责任公司 | 一种基于潮流介数的电网故障预警方法及装置 |
CN113884820A (zh) * | 2021-11-19 | 2022-01-04 | 安徽南瑞中天电力电子有限公司 | 基于电能表的线路阻抗测量方法及线路故障类型分析方法 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103208798B (zh) * | 2013-03-26 | 2014-10-08 | 河海大学 | 一种含风电场电力***概率潮流的计算方法 |
CN104410069B (zh) * | 2014-12-05 | 2017-02-08 | 国家电网公司 | 一种计及响应相关性的动态概率潮流计算方法 |
CN107611979B (zh) * | 2017-09-26 | 2019-08-30 | 华中科技大学 | 一种基于概率潮流计算的电力***运行校正控制方法 |
CN107968409A (zh) * | 2017-11-08 | 2018-04-27 | 中国电力科学研究院有限公司 | 一种考虑不平衡功率分配的概率潮流计算方法及*** |
-
2019
- 2019-05-27 CN CN201910444918.8A patent/CN110048428B/zh not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
CN110048428A (zh) | 2019-07-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Williams et al. | Probabilistic load flow modeling comparing maximum entropy and Gram-Charlier probability density function reconstructions | |
CN105429129B (zh) | 一种考虑网络约束的间歇性能源发电容量置信度评估方法 | |
CN102709908B (zh) | 大规模风电接入电网后的网损预测方法 | |
Li et al. | Stochastic event-triggered cubature Kalman filter for power system dynamic state estimation | |
CN110429648B (zh) | 考虑风速随机波动的小干扰稳定裕度概率评估方法 | |
CN109066688A (zh) | 基于可再生能源不确定性下的概率潮流数据获取方法 | |
Aien et al. | Probabilistic load flow considering wind generation uncertainty | |
CN107968409A (zh) | 一种考虑不平衡功率分配的概率潮流计算方法及*** | |
CN104156542B (zh) | 一种基于隐式投影的有源配电***稳定性仿真方法 | |
CN107730111A (zh) | 一种考虑用户负荷和新能源接入的配网电压风险评估模型 | |
CN110048428B (zh) | 基于概率守恒原理的电力***概率潮流计算方法 | |
CN104992017A (zh) | 计及风电功率随机波动的广义负荷建模方法与装置 | |
Liu et al. | A coordinate-descent algorithm for tracking solutions in time-varying optimal power flows | |
CN107221933B (zh) | 一种概率潮流计算方法 | |
Yang et al. | Multi-objective optimal control approach for static voltage stability of power system considering interval uncertainty of the wind farm output | |
CN109830987A (zh) | 计及分布式光伏随机性的主动配电网概率稳定性分析方法 | |
Rahmouni et al. | Transient stability analysis of the IEEE 39-bus power system using gear and block methods | |
Yong et al. | Investigating the influence of storage on Renewable Energy Capacity Credit | |
Han et al. | Reliability assessment of power system containing wind farm based on steady-state power flow | |
Olamaei et al. | An efficient method for load flow analysis of distribution networks including PV nodes | |
Kayalvizhi et al. | Stochastic optimal power flow in presence of wind generations using harmony search algorithm | |
Peng et al. | Continuous power flow for hybrid AC/DC microgrid considering uncertainty of intermittent DG output | |
CN113113908A (zh) | 适用于现代大电网频率响应的时域解析方法与*** | |
Tang et al. | Reliability evaluation of EGIES with highly nonlinear modeling of gas system components | |
Kan et al. | Calculation of probabilistic harmonic power flow based on improved three-point estimation method and maximum entropy as distributed generators access to distribution network |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20200630 |
|
CF01 | Termination of patent right due to non-payment of annual fee |