[go: up one dir, main page]

CN111753397B - 一种计及分布式电源相关性的配电系统可靠性评估方法 - Google Patents

一种计及分布式电源相关性的配电系统可靠性评估方法 Download PDF

Info

Publication number
CN111753397B
CN111753397B CN202010447877.0A CN202010447877A CN111753397B CN 111753397 B CN111753397 B CN 111753397B CN 202010447877 A CN202010447877 A CN 202010447877A CN 111753397 B CN111753397 B CN 111753397B
Authority
CN
China
Prior art keywords
load
reliability
matrix
island
distributed 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.)
Active
Application number
CN202010447877.0A
Other languages
English (en)
Other versions
CN111753397A (zh
Inventor
宣菊琴
张林垚
郑洁云
吴桂联
何华琴
郑维明
黄东明
黄超
王震
胡志坚
王珏莹
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
State Grid Fujian Electric Power Co Ltd
Economic and Technological Research Institute of State Grid Fujian Electric Power Co Ltd
Original Assignee
XIAMEN GREAT POWER GEO INFORMATION TECHNOLOGY CO LTD
Wuhan University WHU
State Grid Fujian Electric Power Co Ltd
Economic and Technological Research Institute of State Grid Fujian Electric Power Co Ltd
Quanzhou Power Supply Co of State Grid Fujian Electric Power Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by XIAMEN GREAT POWER GEO INFORMATION TECHNOLOGY CO LTD, Wuhan University WHU, State Grid Fujian Electric Power Co Ltd, Economic and Technological Research Institute of State Grid Fujian Electric Power Co Ltd, Quanzhou Power Supply Co of State Grid Fujian Electric Power Co Ltd filed Critical XIAMEN GREAT POWER GEO INFORMATION TECHNOLOGY CO LTD
Priority to CN202010447877.0A priority Critical patent/CN111753397B/zh
Publication of CN111753397A publication Critical patent/CN111753397A/zh
Application granted granted Critical
Publication of CN111753397B publication Critical patent/CN111753397B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING OR CALCULATING; COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G—PHYSICS
    • G06—COMPUTING OR CALCULATING; 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/0639—Performance analysis of employees; Performance analysis of enterprise or organisation operations
    • G06Q10/06393—Score-carding, benchmarking or key performance indicator [KPI] analysis
    • G—PHYSICS
    • G06—COMPUTING OR CALCULATING; 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
    • G—PHYSICS
    • G06—COMPUTING OR CALCULATING; COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/02—Reliability analysis or reliability optimisation; Failure analysis, e.g. worst case scenario performance, failure mode and effects analysis [FMEA]

Landscapes

  • Business, Economics & Management (AREA)
  • Engineering & Computer Science (AREA)
  • Human Resources & Organizations (AREA)
  • Theoretical Computer Science (AREA)
  • Economics (AREA)
  • Physics & Mathematics (AREA)
  • Strategic Management (AREA)
  • General Physics & Mathematics (AREA)
  • Tourism & Hospitality (AREA)
  • Health & Medical Sciences (AREA)
  • Marketing (AREA)
  • Development Economics (AREA)
  • Educational Administration (AREA)
  • Entrepreneurship & Innovation (AREA)
  • General Business, Economics & Management (AREA)
  • Game Theory and Decision Science (AREA)
  • Quality & Reliability (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • Operations Research (AREA)
  • Public Health (AREA)
  • Water Supply & Treatment (AREA)
  • General Health & Medical Sciences (AREA)
  • Primary Health Care (AREA)
  • Supply And Distribution Of Alternating Current (AREA)

Abstract

本发明涉及一种计及分布式电源相关性的配电系统可靠性评估方法，包括如下步骤：步骤S1：考虑分布式电源的不确定性和风力、光照强度之间的相关性，分析风光发电出力规律，基于拉丁超立方抽样，结合Spearman秩相关系数和Cholesky分解获得指定相关系数的分布式电源出力样本；步骤S2：基于IEEE可靠性测试系统的峰值负荷百分比模型,采用负荷聚类技术进行聚类获得各负荷点的负荷情况；步骤S3：基于计及开关故障的双向层级结构的配电系统可靠性评估算法，对系统可靠性进行评估；步骤S4：考虑分布式电源接入后配电网的孤岛运行模式，以负荷切除量最大提出改进的启发式负荷削减策略，并对孤岛内负荷的可靠性指标进行修正。本发明提高了系统的可靠性计算效率。

Description

一种计及分布式电源相关性的配电系统可靠性评估方法
技术领域
本发明涉及分布式电源领域，具体涉及一种计及分布式电源相关性的配电系统可靠性评估方法。
背景技术
随着现代社会的快速发展，用户对电力系统供电可靠性的要求不断提高，如何快速准确地对配电系统可靠性进行评估并可靠地进行供电至关重要。分布式电源(Distributed Generation，DG)发电方式灵活、绿色环保，大量地应用于配电网络。DG的接入改变了配电网的结构和运行方式，与此同时也存在着输出功率具有波动性、不确定性等问题。并且同一地区的风速和光照强度之间具有相关性，如果没有对这一因素加以考虑，会影响配电网的可靠性评估，对配电网的可靠运行带来影响。
目前含DG的配电系统可靠性研究很少将相关性因素考虑进去，导致评估结果不切合工程实际。随着配电网复杂程度的增加，系统元件数量增大，传统可靠性评估算法计算效率低下，并且大多数情况下忽略开关故障的影响。考虑开关故障的双向层级结构可靠性评估算法能够对配电网可靠性进行快速、精确评估，在考虑DG相关性的同时，在孤岛内DG出力不足时采用改进的启发式负荷削减策略，确保切除最大负荷，提高了供电可靠性。
发明内容
有鉴于此，本发明的目的在于提供一种计及分布式电源相关性的配电系统可靠性评估方法，考虑DG的同时，还考虑了DG之间的相关性对配电网可靠性的影响，适应了新形势下电网发展需要，能够更加切合工程实际地对配电网进行可靠性评估；提出考虑开关故障的双向层级结构可靠性评估算法能够对可靠性指标进行快速评估，同时考虑多种因素的影响，提高了评估的精确性；采用的改进的启发式负荷削减策略在DG出力不足时能够改善配电网的可靠性。
为实现上述目的，本发明采用如下技术方案：
一种计及分布式电源相关性的配电系统可靠性评估方法，包括如下步骤：
步骤S1:考虑分布式电源的不确定性和风力、光照强度之间的相关性，分析风光发电出力规律，基于拉丁超立方抽样，结合Spearman秩相关系数和Cholesky分解获得指定相关系数的分布式电源出力样本；
步骤S2:基于IEEE可靠性测试系统的峰值负荷百分比模型,采用负荷聚类技术进行聚类获得各负荷点的负荷情况；
步骤S3：基于计及开关故障的双向层级结构的配电系统可靠性评估算法，对系统可靠性进行评估；
步骤S4：考虑分布式电源接入后配电网的孤岛运行模式，以负荷切除量最大提出改进的启发式负荷削减策略，并对孤岛内负荷的可靠性指标进行修正。
进一步的，所述结合Spearman秩相关系数、拉丁超立方抽样、Cholesky分解生成相关性样本具体为：
分布式电源由风力发电机和光伏发电机组成，二者出力皆具有波动性和不确定性。风力发电输出功率受风速大小影响，大量实测研究数据表明，风速的概率分布符合两参数威布尔分布：
Figure BDA0002506458100000021
Figure BDA0002506458100000022
式中，V是风力发电机组轮毂处的风速；k＞0是形状参数，c＞0是尺度参数；f(V)为风速的概率密度函数；F(V)为风速的累计概率分布函数。
光伏发电输出功率与温度和光照强度均有关，并主要取决于光照强度。大量研究表明光照强度服从Beta分布，其概率分布函数可表示如下：
Figure BDA0002506458100000023
式中，I为光照强度；Imax为光强最大值；α和β为形状参数；Γ(·)为伽玛函数。
步骤S11:将威布尔分布或Beta分布的累计概率分布函数均匀地分为N个区间，对每个区间随机抽取任意值，进行逆变换得到该区间的采样值；依次对M个变量进行上述操作可得M×N阶初始样本矩阵P；
步骤S12:根据变量的历史数据所得到的秩相关性矩阵ρobj，对其进行Cholesky分解，其中H为下三角矩阵；
ρobj＝HHT (4)
步骤S13:随机生成M×N阶矩阵R，矩阵R每行由不大于N且不重复的正整数组成，计算其秩相关系数矩阵ρR后进行Cholesky分解，其中L为分解后的下三角矩阵，并得到中间矩阵G；
ρR＝LLT (5)
G＝L-1R (6)
步骤S14:根据矩阵ρobj的秩相关性将矩阵G更新为Gn，再根据Gn中各行元素的排列顺序将矩阵R更新为Rn；
Gn＝HL-1R (7)
步骤S15:根据Rn的排序顺序将矩阵P更新为Pn，获得与矩阵ρobj相同相关性的多变量样本空间。
进一步的，所述Spearman秩相关系数具体为：
设随机变量(X,Y)的K对样本(xi,yi)，则样本的秩相关系数ρ可表示为：
Figure BDA0002506458100000031
Figure BDA0002506458100000032
Figure BDA0002506458100000033
式中，rxi、ryi分别为随机变量排列后xi、yi所对应的秩；
其中ρ＞0表示Y随X增加而增加；ρ＜0表示Y随X增加而减少；ρ＝0表示X、Y之间没有相关性；并且，|ρ|越大说明随机变量(X,Y)单调变化的关系越强。
进一步的，所述步骤S2具体为：
步骤S21:设负荷集合为Q＝{q1,q2,…,qn}，采用快速排序法对Q进行排序，
并将负荷水平相同的进行合并,不同的负荷水平从小到大排序形成数组L，其长度为Length，并用数组Num存储各个负荷水平合并后包含的负荷点个数，每个负荷水平即为初始聚类负荷中心；
步骤S22:计算L中各个相邻负荷中心的欧式距离h，并设定负荷簇内的距离阈值T：
h(i)＝|L(i+1)-L(i)| (11)
T＝5[max(L)-min(L)]/Length (12)
步骤S23:将距离最近的一对负荷值归为第一类CL1，更新聚类后的聚类中心及负荷的数目；
Figure BDA0002506458100000041
Figure BDA0002506458100000042
式中：Ici为负荷簇i中的负荷点集合；
步骤S24:继续计算剩余负荷与负荷簇CL1聚类中心的欧式距离，将距离小于T的负荷归并至CL1中，并将已经归并的负荷从负荷集合中删去；
步骤S25:在剩余的负荷找到距离最近的一对负荷，设最小距离为hm，若hm＜T，说明这对负荷属于同一个负荷簇，将其归并至CLi中，并按照步骤S24所述的方法将负荷距离小于T的负荷归并至CLi中，否则转向步骤S26；
步骤S26:这对负荷分别归并至不同的负荷簇，继续步骤S25；
步骤S27:如果最后只剩下一个负荷，则将该负荷归并为一类。
进一步的，所述步骤S3具体为：
步骤S31:以断路器或隔离开关为界进行分块；
步骤S32:计算元件块等效故障率和故障持续时间；
步骤S33:分支线路的首端断路器融合在父块中，其它的开关故障融合在下游元件块中，顺向传递可靠性参数；
步骤S34:开关故障累加至逆向相邻下游元件块；
步骤S35:开关故障融合在逆向下游相邻元件块中，并随着隔离开关块逆向传递可靠性参数；
步骤S36:计算可靠性指标。
进一步的，所述可靠性指标包括负荷点可靠性指标和系统可靠性指标。
进一步的，所述负荷点可靠性指标包括三个指标：年平均停电频率λi、年平均停电时间Ui和平均停电持续时间ri。
进一步的，系统可靠性指标包含如下：
(1)系统平均停电频率指标
Figure BDA0002506458100000051
式中，n表示负荷点总数，Ni表示第i个负荷点包含的用户数；
(2)系统平均停电持续时间指标
Figure BDA0002506458100000052
(3)用户平均停电持续时间指标
Figure BDA0002506458100000053
(4)平均供电可用率指标
Figure BDA0002506458100000054
(5)系统期望缺供电量指标
Figure BDA0002506458100000055
式中，Lai表示第i个负荷点的年平均负荷功率。
进一步的，所述步骤S4：
步骤S41:考虑分布式电源接入后配电网的孤岛运行模式，以负荷切除量最大提出改进的启发式负荷削减策略；
步骤S42:当分布式电源出力不足时，优先削减负荷总量最大的负荷点，如果分布式电源出力仍然不足，则继续削减负荷总量最大的负荷点，直至满足要求为止；
步骤S43:对孤岛内负荷的可靠性指标进行修正：
Figure BDA0002506458100000061
Figure BDA0002506458100000062
式中：PCi为负荷削减策略下负荷i恢复供电的概率，λ3、U3分别为负荷削减策略后的负荷故障率和年平均故障时间。
进一步的，所述步骤S41具体为：分布式电源的出力具有不确定性的特点，只有当孤岛内分布式电源总出力大于孤岛内总负荷时，孤岛才能成功形成，孤岛成功形成的概率PIS：
Figure BDA0002506458100000063
式中，N为采样空间，∑PDGi、∑PLi分别为孤岛内第i个场景DG总出力及总负荷；
设LPi为孤岛内的任意负荷，λ1、λ2分别表示修正前后LPi的年平均停电频率，U1、U2分别表示修正前后LPi的年平均故障时间，DG接入后LPi的可靠性指标修正为：
Figure BDA0002506458100000064
Figure BDA0002506458100000065
式中，m、λki和Uki分别表示孤岛上游主馈线元件个数、各元件故障率和年平均故障时间。
本发明与现有技术相比具有以下有益效果：
1、本发明在进行配电网可靠性评估时，考虑了DG相关性的影响，更加切合工程实际；
2、本发明采用了计及开关故障的双向层级结构的配电系统可靠性评估算法，评估速度更快，并且考虑了开关故障的影响，评估结果更加精确；
3、本发明采用了改进的启发式负荷削减策略，在DG出力不足时能够改善配电网的可靠性。
附图说明
图1为本发明计及DG相关性的可靠性评估流程。
图2为本发明计及开关故障的双向层级结构可靠性评估算法流程。
图3为本发明IEEERBTS BUS6系统图。
图4为本发明系统部分负荷点的年平均停电频率变化情况。
图5为本发明不同场景的系统缺供电量随相关性的变化情况。
具体实施方式
下面结合附图及实施例对本发明做进一步说明。
如附图1和附图2所示，本发明的一种计及DG相关性的配电系统可靠性评估方法，具体实现如下：
本发明所述技术方案的步骤S1的具体实现如下：
首先假设DG主要由风力发电机(WindTurbine Generator，WTG)以及光伏发电机(Photovoltaic Generator，PVG)组成，二者出力皆具有波动性以及不确定性，可以分别用威布尔分布和Beta分布进行表征。
Spearman秩相关系数可以表征DG间出力的相关性。设随机变量(X,Y)的K对样本(xi,yi)，则样本的秩相关系数ρ可表示为：
Figure BDA0002506458100000071
Figure BDA0002506458100000072
Figure BDA0002506458100000073
式中，rxi、ryi分别为随机变量排列后xi、yi所对应的秩。
其中ρ＞0表示Y随X增加而增加；ρ＜0表示Y随X增加而减少；ρ＝0表示X、Y之间没有相关性。并且，|ρ|越大说明随机变量(X,Y)单调变化的关系越强。
结合Spearman秩相关系数、拉丁超立方抽样、Cholesky分解生成相关性样本的步骤如下：
(1)将威布尔分布或Beta分布的累计概率分布函数均匀地分为N个区间，对每个区间随机抽取任意值，进行逆变换得到该区间的采样值；依次对M个变量进行上述操作可得M×N阶初始样本矩阵P；
(2)根据变量的历史数据所得到的秩相关性矩阵ρobj，对其进行Cholesky分解，其中H为下三角矩阵；
ρobj＝HHT (4)
(3)随机生成M×N阶矩阵R，矩阵R每行由不大于N且不重复的正整数组成，计算其秩相关系数矩阵ρR后进行Cholesky分解，其中L为分解后的下三角矩阵，并得到中间矩阵G；
ρR＝LLT (5)
G＝L-1R (6)
(4)根据矩阵ρobj的秩相关性将矩阵G更新为Gn，再根据Gn中各行元素的排列顺序将矩阵R更新为Rn；
Gn＝HL-1R (7)
(5)根据Rn的排序顺序将矩阵P更新为Pn，最终获得与矩阵ρobj相同相关性的多变量样本空间。
在本实施例中，所述步骤S2的具体实现如下：
年负荷曲线所包含的负荷等级众多，为简化计算，可以采用负荷聚类技术进行归类。设负荷集合为Q＝{q1,q2,…,qn}，负荷聚类的具体流程如下：
(1)采用快速排序法对Q进行排序，并将负荷水平相同的进行合并。不同的负荷水平从小到大排序形成数组L，其长度为Length，并用数组Num存储各个负荷水平合并后包含的负荷点个数，每个负荷水平即为初始聚类负荷中心；
(2)计算L中各个相邻负荷中心的欧式距离h，并设定负荷簇内的距离阈值T：
h(i)＝|L(i+1)-L(i)| (8)
T＝5[max(L)-min(L)]/Length (9)
(3)将距离最近的一对负荷值归为第一类CL1，更新聚类后的聚类中心以及负荷的数目；
Figure BDA0002506458100000091
Figure BDA0002506458100000092
式中：Ici为负荷簇i中的负荷点集合。
(4)继续计算剩余负荷与负荷簇CL1聚类中心的欧式距离，将距离小于T的负荷归并至CL1中，并将已经归并的负荷从负荷集合中删去；
(5)在剩余的负荷找到距离最近的一对负荷，设最小距离为hm，若hm＜T，说明这对负荷属于同一个负荷簇，将其归并至CLi中，并按照(4)所述的方法将负荷距离小于T的负荷归并至CLi中，否则转向(6)；
(6)这对负荷分别归并至不同的负荷簇，继续(5)；
(7)如果最后只剩下一个负荷，则将该负荷归并为一类。
在本实施例中，所述步骤S3的具体实现如下：
考虑开关故障的双向层级结构可靠性评估算法将开关作为边界对系统进行分块，计算效率高，并且考虑了开关故障的影响，能够对系统进行精确的可靠性评估。
配电系统的可靠性指标包括负荷点可靠性指标和系统可靠性指标。负荷点可靠性指标包括三个指标：年平均停电频率λi、年平均停电时间Ui和平均停电持续时间ri。系统可靠性指标包括：
(1)系统平均停电频率指标(system average Interruption Frequency Index，SAIFI)。
Figure BDA0002506458100000101
式中，n表示负荷点总数，Ni表示第i个负荷点包含的用户数。
(2)系统平均停电持续时间指标(system average interruption durationindex，SAIDI)。
Figure BDA0002506458100000102
(3)用户平均停电持续时间指标(customer average interruption durationindex，CAIDI)。
Figure BDA0002506458100000103
(4)平均供电可用率指标(average service availability index，ASAI)。
Figure BDA0002506458100000104
(5)系统期望缺供电量指标(energy not supplied index，ENSI)。
Figure BDA0002506458100000105
式中，Lai表示第i个负荷点的年平均负荷功率。
有了上述评估系统的可靠性指标后，就可以进行可靠性评估。考虑开关故障的双向层级结构可靠性评估算法的具体步骤如下：
(1)以断路器或隔离开关为界进行分块；
(2)计算元件块等效故障率和故障持续时间；
(3)分支线路的首端断路器融合在父块中，其它的开关故障融合在下游元件块中，顺向传递可靠性参数；
(4)开关故障累加至逆向相邻下游元件块；
(5)开关故障融合在逆向下游相邻元件块中，并随着隔离开关块逆向传递可靠性参数；
(6)计算可靠性指标。
在本实施例中，所述步骤S4的具体实现如下：
(1)孤岛运行
配电网发生故障时，DG只对孤岛内负荷继续供电，因此DG只影响孤岛内负荷的可靠性。
考虑到DG的出力具有不确定性的特点，只有当孤岛内DG总出力大于孤岛内总负荷时，孤岛才能成功形成，孤岛成功形成的概率PIS：
Figure BDA0002506458100000111
式中，N为采样空间，∑PDGi、∑PLi分别为孤岛内第i个场景DG总出力及总负荷。
只有当孤岛上游主馈线元件发生故障时，才形成孤岛运行，只需对孤岛上游主馈线元件发生故障情况进行修正即可。设LPi为孤岛内任意负荷，λ1、λ2分别表示修正前后LPi的年平均停电频率，U1、U2分别表示修正前后LPi的年平均故障时间，DG接入后LPi的可靠性指标修正为：
Figure BDA0002506458100000112
Figure BDA0002506458100000113
式中，m、λki和Uki分别表示孤岛上游主馈线元件个数、各元件故障率和年平均故障时间。
(2)负荷削减策略
由于DG的出力具有不确定性，发生故障时并不能保证任意场景下都能给孤岛内所有负荷进行供电，因此就需要切除部分负荷。
本实施例采用改进的启发式负荷削减策略，具体为：为了使DG能最大化地为负荷进行供电，当DG出力不足时，优先削减负荷总量最大的负荷点，如果DG出力仍然不足，则继续削减负荷总量最大的负荷点，直至满足要求为止。考虑负荷削减策略后，需要对孤岛内负荷的可靠性指标进行修正：
Figure BDA0002506458100000121
Figure BDA0002506458100000122
式中：PCi为负荷削减策略下负荷i恢复供电的概率，λ3、U3分别为负荷削减策略后的负荷故障率和年平均故障时间。
实施例1
如附图3所示，本实施例采用改造的IEEE RBTS BUS6系统。该系统包括1条母线，23个熔断器，23台配电变压器，23个负荷，30条线路，4台断路器和1个隔离开关。线路长度、各元件可靠性数据、负荷等具体信息分别如表1-表6所示。
表1线路长度
Figure BDA0002506458100000123
表2各元件可靠性参数
Figure BDA0002506458100000124
表3各负荷点负荷参数
Figure BDA0002506458100000125
Figure BDA0002506458100000131
表4各周负荷峰值与年负荷峰值之比
Figure BDA0002506458100000132
表5各日负荷峰值与周负荷峰值之比
Figure BDA0002506458100000133
表6各小时负荷峰值与日负荷峰值之比
Figure BDA0002506458100000134
Figure BDA0002506458100000141
设置威布尔分布的形状参数为3.97，尺度参数10.7，风机切入风速3m/s、额定风速14m/s、切出风速25m/s；Beta分布的形状参数α＝2.0、β＝0.8，PVG的最大光照强度为600W/m2，在馈线53、59处均安装了DG，形成计划孤岛；拉丁超立方抽样规模为1000。
在馈线53、59处均安装1台WTG、1台PVG，WTG与PVG的额定容量相同。风速之间的相关系数为0.8，光照强度之间的相关系数为0.8，风速与光照强度之间的相关系数为-0.6。在不同的DG额定容量情况下，系统部分负荷点的年平均停电频率变化情况如附图4所示。由附图4可知：
(1)随着DG额定容量的增加，孤岛内各负荷点的年平均停电频率逐渐减小。可知DG的接入可以有效提高孤岛内负荷的可靠性指标；
(2)随着DG额定容量的增加，LP16负荷点的年平均停电频率不变。这是因为DG的接入只会提升孤岛范围内负荷点的可靠性；
(3)由于LP8负荷点的负荷功率较小，LP23负荷点的负荷功率较大，采用的改进的启发式负荷削减策略在DG出力不足时优先削减负荷功率较大的负荷点，因此随着DG额定容量的增大，LP8负荷点相对于LP23负荷点的年平均停电频率改善效果更为明显；
(4)在DG额定容量的增大的初始阶段，孤岛内负荷点年平均停电频率迅速减小，而随着DG额定容量不断增大，孤岛内负荷点年平均停电频率减小缓慢，直至饱和。
为了更好地对DG相关性进行研究，设置以下三种场景。设置每台DG的额定容量均为1.2MW。
(1)场景一：所有DG均为WTG，风速之间的相关性依次为0.1、0.2、0.3、0.4、0.5、0.6、0.7、0.8。
(2)场景二：所有DG均为PVG，光照强度之间的相关性依次为0.1、0.2、0.3、0.4、0.5、0.6、0.7、0.8。
(3)场景三：每个孤岛均接入一台WTG、一台PVG，风速之间的相关性为0.8，光照强度之间的相关性为0.8，风速与光照强度之间的相关性依次为-0.1、-0.2、-0.3、-0.4、-0.5、-0.6、-0.7、-0.8。
三种场景的ENSI指标随相关性的变化如附图5所示，并且由附图5可知：
(1)根据场景一可知，随着风速与风速的正相关性增加，ENSI逐渐增大。这是因为风速与风速的正相关性越大，风机输出功率同时增大或减小的可能性增大，导致系统的可靠性降低；
(2)根据场景二可知，随着光照强度与光照强度的正相关性增加，ENSI逐渐增大。这是因为光照强度与光照强度的正相关性越大，PVG输出功率同时增大或减小的可能性增大，导致系统的可靠性降低；
(3)根据场景三可知，随着风速和光照强度的负相关性增加，ENSI逐渐减小。这是因为风速和光照强度的负相关性越大，WTG和PVG的输出功率同时增大或减小的可能性减小，所以系统可靠性提升。
将本文采用的改进的启发式负荷削减策略(策略1)与考虑负荷重要系数的削减策略(策略2)进行可靠性对比分析，结果如表7所示。
表7不同负荷削减策略的系统可靠性指标
Figure BDA0002506458100000161
由表7可知，策略1与策略2相比，采用策略1所产生的系统可靠性更好。这是因为每次进行负荷削减时，策略1每次只削减负荷最大的负荷点，有时候能够避免削减多处负荷，因此可靠性指标更好。
以上所述仅为本发明的较佳实施例，凡依本发明申请专利范围所做的均等变化与修饰，皆应属本发明的涵盖范围。

Claims (8)

1.一种计及分布式电源相关性的配电系统可靠性评估方法，其特征在于，包括如下步骤：
步骤S1:考虑分布式电源的不确定性和风力、光照强度之间的相关性，分析风光发电出力规律，基于拉丁超立方抽样，结合Spearman秩相关系数和Cholesky分解获得指定相关系数的分布式电源出力样本；
步骤S2:基于IEEE可靠性测试系统的峰值负荷百分比模型，采用负荷聚类技术进行聚类获得各负荷点的负荷情况；
步骤S3：基于计及开关故障的双向层级结构的配电系统可靠性评估算法，对系统可靠性进行评估；
步骤S4：考虑分布式电源接入后配电网的孤岛运行模式，以负荷切除量最大提出改进的启发式负荷削减策略，并对孤岛内负荷的可靠性指标进行修正；
所述步骤S4：
步骤S41:考虑分布式电源接入后配电网的孤岛运行模式，以负荷切除量最大提出改进的启发式负荷削减策略；
步骤S42:当分布式电源出力不足时，优先削减负荷总量最大的负荷点，如果分布式电源出力仍然不足，则继续削减负荷总量最大的负荷点，直至满足要求为止；
步骤S43:对孤岛内负荷的可靠性指标进行修正：
Figure FDA0003934061080000011
Figure FDA0003934061080000012
式中：PCi为负荷削减策略下负荷i恢复供电的概率，λ3、U3分别为负荷削减策略后的负荷故障率和年平均故障时间；
所述步骤S41具体为：分布式电源的出力具有不确定性的特点，只有当孤岛内分布式电源总出力大于孤岛内总负荷时，孤岛才能成功形成，孤岛成功形成的概率PIS：
Figure FDA0003934061080000021
式中，N为采样空间，∑PDGi、∑PLi分别为孤岛内第i个场景DG总出力及总负荷；
设LPi为孤岛内的任意负荷，λ1、λ2分别表示修正前后LPi的年平均停电频率，U1、U2分别表示修正前后LPi的年平均故障时间，DG接入后LPi的可靠性指标修正为：
Figure FDA0003934061080000022
Figure FDA0003934061080000023
式中，m、λki和Uki分别表示孤岛上游主馈线元件个数、各元件故障率和年平均故障时间。
2.根据权利要求1所述的一种计及分布式电源相关性的配电系统可靠性评估方法，其特征在于，所述结合Spearman秩相关系数、拉丁超立方抽样、Cholesky分解生成相关性样本具体为：
步骤S11:将威布尔分布或Beta分布的累计概率分布函数均匀地分为N个区间，对每个区间随机抽取任意值，进行逆变换得到该区间的采样值；依次对M个变量进行上述操作可得M×N阶初始样本矩阵P；
步骤S12:根据变量的历史数据所得到的秩相关性矩阵ρobj，对其进行Cholesky分解，其中H为下三角矩阵；
ρobj＝HHT (4)
步骤S13:随机生成M×N阶矩阵R，矩阵R每行由不大于N且不重复的正整数组成，计算其秩相关系数矩阵ρR后进行Cholesky分解，其中L为分解后的下三角矩阵，并得到中间矩阵G；
ρR＝LLT (5)
G＝L-1R (6)
步骤S14:根据矩阵ρobj的秩相关性将矩阵G更新为Gn，再根据Gn中各行元素的排列顺序将矩阵R更新为Rn；
Gn＝HL-1R (7)
步骤S15:根据Rn的排序顺序将矩阵P更新为Pn，获得与矩阵ρobj相同相关性的多变量样本空间。
3.根据权利要求1所述的一种计及分布式电源相关性的配电系统可靠性评估方法，其特征在于，所述Spearman秩相关系数具体为：
设随机变量(X,Y)的K对样本(xi,yi)，则样本的秩相关系数ρ可表示为：
Figure FDA0003934061080000031
Figure FDA0003934061080000032
Figure FDA0003934061080000033
式中，rxi、ryi分别为随机变量排列后xi、yi所对应的秩；
其中ρ＞0表示Y随X增加而增加；ρ＜0表示Y随X增加而减少；ρ＝0表示X、Y之间没有相关性；并且，|ρ|越大说明随机变量(X,Y)单调变化的关系越强。
4.根据权利要求1所述的一种计及分布式电源相关性的配电系统可靠性评估方法，其特征在于，所述步骤S2具体为：
步骤S21:设负荷集合为Q＝{q1,q2,…,qn}，采用快速排序法对Q进行排序，并将负荷水平相同的进行合并,不同的负荷水平从小到大排序形成数组L，其长度为Length，并用数组Num存储各个负荷水平合并后包含的负荷点个数，每个负荷水平即为初始聚类负荷中心；
步骤S22:计算L中各个相邻负荷中心的欧式距离h(i)，并设定负荷簇内的距离阈值T：
h(i)＝|L(i+1)-L(i)| (11)
T＝5[max(L)-min(L)]/Length (12)
步骤S23:将距离最近的一对负荷值归为第一类CL1，更新聚类后的聚类中心及负荷的数目；
Figure FDA0003934061080000041
Figure FDA0003934061080000042
式中：Ici为负荷簇i中的负荷点集合；
步骤S24:继续计算剩余负荷与负荷簇CL1聚类中心的欧式距离，将距离小于T的负荷归并至CL1中，并将已经归并的负荷从负荷集合中删去；
步骤S25:在剩余的负荷找到距离最近的一对负荷，设最小距离为hm，若hm＜T，说明这对负荷属于同一个负荷簇，将其归并至CLi中，并按照步骤S24所述的方法将负荷距离小于T的负荷归并至CLi中，否则转向步骤S26；
步骤S26:这对负荷分别归并至不同的负荷簇，继续步骤S25；
步骤S27:如果最后只剩下一个负荷，则将该负荷归并为一类。
5.根据权利要求1所述的一种计及分布式电源相关性的配电系统可靠性评估方法，其特征在于，所述步骤S3具体为：
步骤S31:以断路器或隔离开关为界进行分块；
步骤S32:计算元件块等效故障率和故障持续时间；
步骤S33:分支线路的首端断路器融合在父块中，其它的开关故障融合在下游元件块中，顺向传递可靠性参数；
步骤S34:开关故障累加至逆向相邻下游元件块；
步骤S35:开关故障融合在逆向下游相邻元件块中，并随着隔离开关块逆向传递可靠性参数；
步骤S36:计算可靠性指标。
6.根据权利要求5所述的一种计及分布式电源相关性的配电系统可靠性评估方法，其特征在于，所述可靠性指标包括负荷点可靠性指标和系统可靠性指标。
7.根据权利要求6所述的一种计及分布式电源相关性的配电系统可靠性评估方法，其特征在于，所述负荷点可靠性指标包括三个指标：年平均停电频率λi、年平均停电时间Ui和平均停电持续时间ri。
8.根据权利要求6所述的一种计及分布式电源相关性的配电系统可靠性评估方法，其特征在于，系统可靠性指标包含如下：
(1)系统平均停电频率指标
Figure FDA0003934061080000051
式中，n表示负荷点总数，Ni表示第i个负荷点包含的用户数；
(2)系统平均停电持续时间指标
Figure FDA0003934061080000052
(3)用户平均停电持续时间指标
Figure FDA0003934061080000053
(4)平均供电可用率指标
Figure FDA0003934061080000054
(5)系统期望缺供电量指标
Figure FDA0003934061080000055
式中，Lai表示第i个负荷点的年平均负荷功率。
CN202010447877.0A 2020-05-25 2020-05-25 一种计及分布式电源相关性的配电系统可靠性评估方法 Active CN111753397B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010447877.0A CN111753397B (zh) 2020-05-25 2020-05-25 一种计及分布式电源相关性的配电系统可靠性评估方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010447877.0A CN111753397B (zh) 2020-05-25 2020-05-25 一种计及分布式电源相关性的配电系统可靠性评估方法

Publications (2)

Publication Number Publication Date
CN111753397A CN111753397A (zh) 2020-10-09
CN111753397B true CN111753397B (zh) 2023-05-05

Family

ID=72673866

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010447877.0A Active CN111753397B (zh) 2020-05-25 2020-05-25 一种计及分布式电源相关性的配电系统可靠性评估方法

Country Status (1)

Country Link
CN (1) CN111753397B (zh)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112821451B (zh) * 2021-01-11 2023-05-09 国网福建省电力有限公司泉州供电公司 基于需求侧管理与储能的城镇配电网光伏接入应对方法
CN112685963A (zh) * 2021-01-27 2021-04-20 重庆大学 基于数据驱动的电力系统可靠性逆问题建模与计算方法
CN113128071B (zh) * 2021-05-08 2024-02-09 南京工程学院 一种含光伏发电的发电系统可靠性评估方法
CN113205266B (zh) * 2021-05-12 2022-06-14 浙江浙能技术研究院有限公司 分布式区域能源布局用规划与评估系统

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5798939A (en) * 1995-03-31 1998-08-25 Abb Power T&D Company, Inc. System for optimizing power network design reliability
CN101059555A (zh) * 2006-04-17 2007-10-24 盈正豫顺电子股份有限公司 分布式电源系统孤岛运转的检测装置及其方法
CN103855707A (zh) * 2014-02-20 2014-06-11 深圳供电局有限公司 一种含分布式电源配电网的供电可靠性评估方法
CN105406470A (zh) * 2015-12-21 2016-03-16 国家电网公司 一种基于开关边界分区的有源配电网可靠性评估方法
CN109103874A (zh) * 2017-11-27 2018-12-28 广东电网有限责任公司佛山供电局 考虑负荷特性与分布式电源接入的配电网可靠性评估方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5798939A (en) * 1995-03-31 1998-08-25 Abb Power T&D Company, Inc. System for optimizing power network design reliability
CN101059555A (zh) * 2006-04-17 2007-10-24 盈正豫顺电子股份有限公司 分布式电源系统孤岛运转的检测装置及其方法
CN103855707A (zh) * 2014-02-20 2014-06-11 深圳供电局有限公司 一种含分布式电源配电网的供电可靠性评估方法
CN105406470A (zh) * 2015-12-21 2016-03-16 国家电网公司 一种基于开关边界分区的有源配电网可靠性评估方法
CN109103874A (zh) * 2017-11-27 2018-12-28 广东电网有限责任公司佛山供电局 考虑负荷特性与分布式电源接入的配电网可靠性评估方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
基于动态孤岛混合整数线性规划模型的主动配电网可靠性分析;孙雨潇等;《电力建设》;20190501(第05期);全文 *

Also Published As

Publication number Publication date
CN111753397A (zh) 2020-10-09

Similar Documents

Publication Publication Date Title
CN111753397B (zh) 一种计及分布式电源相关性的配电系统可靠性评估方法
CN104376504B (zh) 一种基于解析法的配电系统概率可靠性评估方法
CN109103874B (zh) 考虑负荷特性与分布式电源接入的配电网可靠性评估方法
CN111861029A (zh) 一种计及孤岛划分和网络重构的供电可靠性评估方法
CN109146124B (zh) 一种基于时变失效率的配电终端改造决策方法
CN106354985B (zh) 一种考虑分布式电源的配电系统可靠性评估方法
US20170271880A1 (en) Resiliency controller for frequency regulation in islanded microgrid
CN116345565A (zh) 一种新能源与储能容量联合优化方法、系统、设备和介质
CN112865089A (zh) 一种改进的主动配电网大规模场景分析方法
CN119362409A (zh) 基于贝叶斯网络与序贯蒙特卡洛法的微网配电系统可靠性评估方法
CN104636993B (zh) 配电系统可靠性算法
CN117937620A (zh) 计及新能源出力不确定性和相继故障风险的电力调度方法
CN117741456B (zh) 一种动态可重构电池网络故障诊断方法、系统及电子设备
CN116937695B (zh) 一种光伏电源的数据分析方法、装置、设备及存储介质
CN120511764A (zh) 一种非预期性电力系统日前-日内联合优化调度方法
CN109586309B (zh) 基于大数据自由熵理论和场景匹配的配电网无功优化方法
CN118735073A (zh) 配网电储能能量时间尺度平移能力指数预测方法
CN110707712B (zh) 一种考虑需求响应的微电网供电可靠性分析方法
CN116307838A (zh) 一种电网可靠性评估中配电网精确等值方法及系统、介质
CN115986836A (zh) 计及场景概率及相关性的输电网络灵活性评价方法及系统
CN117977655B (zh) 海上风电储能配置方法及装置、存储介质、终端
CN119518968B (zh) 煤炭场景不确定性影响下风光储发电系统出力优化方法
CN116976060B (zh) 分布式新能源集群-负荷-网络动态重构的灵活组网方法
CN110854928B (zh) 面向分布式电源和电动汽车的大型配电网风险控制优化方法
CN120875805A (zh) 一种配电网气候适应性重构方法

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
TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20230818

Address after: 350003 No. 257, 54 Road, Gulou District, Fuzhou, Fujian

Patentee after: STATE GRID FUJIAN ELECTRIC POWER Co.,Ltd.

Patentee after: STATE GRID FUJIAN ECONOMIC Research Institute

Address before: 350003 No. 257, 54 Road, Gulou District, Fuzhou, Fujian

Patentee before: STATE GRID FUJIAN ELECTRIC POWER Co.,Ltd.

Patentee before: STATE GRID FUJIAN ECONOMIC Research Institute

Patentee before: Quanzhou Power Supply Company of State Grid Fujian Electric Power Company

Patentee before: XIAMEN GREAT POWER GEO INFORMATION TECHNOLOGY Co.,Ltd.

Patentee before: WUHAN University