CN112383045A - 新能源发电不确定性的暂态稳定越限概率计算方法和装置 - Google Patents
新能源发电不确定性的暂态稳定越限概率计算方法和装置 Download PDFInfo
- Publication number
- CN112383045A CN112383045A CN202011053005.2A CN202011053005A CN112383045A CN 112383045 A CN112383045 A CN 112383045A CN 202011053005 A CN202011053005 A CN 202011053005A CN 112383045 A CN112383045 A CN 112383045A
- Authority
- CN
- China
- Prior art keywords
- new energy
- transient stability
- energy station
- active
- expected
- 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
- 230000001052 transient effect Effects 0.000 title claims abstract description 199
- 238000010248 power generation Methods 0.000 title claims abstract description 42
- 238000004364 calculation method Methods 0.000 title claims abstract description 36
- 238000005070 sampling Methods 0.000 claims abstract description 38
- 238000000034 method Methods 0.000 claims abstract description 28
- 102100036790 Tubulin beta-3 chain Human genes 0.000 claims description 6
- 238000011156 evaluation Methods 0.000 claims description 5
- 102100036788 Tubulin beta-4A chain Human genes 0.000 claims description 4
- 238000011002 quantification Methods 0.000 claims description 3
- 238000011158 quantitative evaluation Methods 0.000 abstract description 4
- 230000006870 function Effects 0.000 description 12
- 238000010586 diagram Methods 0.000 description 8
- 238000012545 processing Methods 0.000 description 6
- 238000004458 analytical method Methods 0.000 description 5
- 238000004590 computer program Methods 0.000 description 5
- 238000012163 sequencing technique Methods 0.000 description 3
- 230000009467 reduction Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000001186 cumulative effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000004880 explosion Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000013139 quantization Methods 0.000 description 1
- 230000006641 stabilisation Effects 0.000 description 1
- 238000011105 stabilization Methods 0.000 description 1
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
-
- 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/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; 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
- 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
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/38—Arrangements for parallely feeding a single network by two or more generators, converters or transformers
- H02J3/381—Dispersed generators
-
- 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/10—Power transmission or distribution systems management focussing at grid-level, e.g. load flow analysis, node profile computation, meshed network optimisation, active network management or spinning reserve management
-
- 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]
-
- 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
- H02J2300/00—Systems for supplying or distributing electric power characterised by decentralized, dispersed, or local generation
- H02J2300/20—The dispersed energy generation being of renewable origin
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
- Y04S—SYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
- Y04S10/00—Systems supporting electrical power generation, transmission or distribution
- Y04S10/50—Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications
Landscapes
- Business, Economics & Management (AREA)
- Engineering & Computer Science (AREA)
- Human Resources & Organizations (AREA)
- Economics (AREA)
- Power Engineering (AREA)
- Strategic Management (AREA)
- Entrepreneurship & Innovation (AREA)
- Theoretical Computer Science (AREA)
- Marketing (AREA)
- Health & Medical Sciences (AREA)
- Educational Administration (AREA)
- Tourism & Hospitality (AREA)
- Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- General Physics & Mathematics (AREA)
- Development Economics (AREA)
- General Health & Medical Sciences (AREA)
- Water Supply & Treatment (AREA)
- Public Health (AREA)
- Game Theory and Decision Science (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Primary Health Care (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明公开一种考虑新能源发电不确定性的暂态稳定越限概率计算方法和装置,该方法包括:根据暂态稳定关键模式的暂态稳定裕度和新能源场站参与因子,计算新能源场站对电网暂态稳定的综合影响因子,确定主导新能源场站;根据新能源场站有功概率分布密度函数、有功实时计划值及其对电网暂态稳定的综合影响因子,确定主导新能源场站的有功采样区间,生成各主导新能源场站有功分档区间组合;计算各组合对应的电网运行状态下预想故障的暂态稳定裕度,得到暂态稳定裕度越限组合;根据预想故障发生概率和暂态稳定裕度越限组合的概率,计算新能源发电不确定性的暂态稳定越限概率。本发明实现了电网暂态稳定越限风险的量化评估。
Description
技术领域
本发明涉及一种新能源发电不确定性的暂态稳定越限概率计算 方法和装置,属于电力***运行与控制技术领域。
背景技术
电网的安全稳定性随电网运行状态变化而变化,新能源发电不确 定性会导致电网运行状态的不确定,新能源并网下电网的安全稳定性 呈现出概率特性。随着新能源发电占比逐渐增大,电网运行状态的不 确定性会逐渐加剧,电网安全稳定的概率特性已不容忽略,电网安全 稳定越限概率成为调度运行控制人员监控电网安全风险的一项重要 指标。
对于新能源发电不确定性的处理,通常采用场景法。通过新能源 场站发电离散化,生成新能源场站发电组合场景,采用针对确定的电 网运行状态的安全稳定分析方法对组合场景进行安全稳定评估,实现 新能源发电不确定性的电网安全稳定概率特性分析。在新能源场站数 极少的情况下,现有方法的计算速度能够满足在线安全稳定分析的实 时性要求。但是,随着新能源场站数目增加,新能源场站发电组合场 景数急剧增加,存在组合***问题。现有方法不能满足大量新能源场 站接入的实际电网在线安全稳定分析对计算速度的要求。
发明内容
针对现有技术的不足,本发明的目的在于提供一种考虑新能源发 电不确定性的暂态稳定越限概率计算方法和装置,以解决现有方法不 能满足大量新能源场站接入的实际电网在线安全稳定分析计算速度 要求的问题。
本发明具体采用如下技术方案:新能源发电不确定性的暂态稳定 越限概率计算方法,所述方法包括:
步骤SS1:针对设定时刻t的电网预想运行状态S,进行预想故障 集F中各预想故障下暂态稳定量化评估,得到预想故障集F中各预想 故障的暂态稳定关键模式集M及暂态稳定关键模式集M中各关键模 式对应的暂态稳定裕度,和新能源场站集N中各新能源场站的暂态 稳定参与因子;
步骤SS2:针对预想故障集F中各预想故障,分别计算预想故障 下的新能源场站集N中各新能源场站对电网暂态稳定的综合影响因 子,生成预想故障关联的主导新能源场站集D;
步骤SS3:针对预想故障集F中各预想故障,分别根据t时刻主 导新能源场站集D中各新能源场站有功概率分布密度函数、有功实 时计划值及各新能源场站对电网暂态稳定的综合影响因子,确定预想 故障下D中各新能源场站的有功采样区间,生成主导新能源场站集D 中各新能源场站有功分档区间组合集C;
步骤SS4:针对预想故障集F中各预想故障,分别计算有功分档 区间组合集C中各组合对应的电网运行状态下预想故障的暂态稳定 裕度,得到预想故障下暂态稳定裕度越限组合集U;
步骤SS5:根据t时刻预想故障集F中各预想故障发生概率和预 想故障下暂态稳定裕度越限组合集U中各组合的概率,计算考虑新 能源发电不确定性的暂态稳定越限概率。
作为一种较佳的实施例,所述步骤SS1具体包括:预想故障集F 中暂态稳定关键模式集M为暂态稳定裕度与所有暂态稳定模式中暂 态稳定裕度最小值之差小于设定值的暂态稳定模式组成的集合。
作为一种较佳的实施例,所述步骤SS2具体包括:
预想故障下的新能源场站集N中各新能源场站对电网暂态稳定 综合影响因子的计算公式为:
λn=∑m∈M[λm.n(1-ηm)α],n∈N (1)
式中,λn为预想故障下的新能源场站集N中新能源场站n对电 网暂态稳定的综合影响因子,λm.n为该预想故障的暂态稳定关键模式 集M中关键模式m下新能源场站集N中新能源场站n的暂态稳定参 与因子,ηm为该预想故障的暂态稳定关键模式集M中关键模式m的暂态稳定裕度,取值范围为[-1,1],α为设定参数,大于0。
作为一种较佳的实施例,所述步骤SS2还包括:
生成预想故障关联的主导新能源场站集D的方法为:将预想故 障下不满足如下公式(2)的新能源场站从新能源场站集N中剔除后, 作为与该预想故障关联的主导新能源场站集D;
式中,β为设定参数,大于0,小于1。
作为一种较佳的实施例,所述步骤SS2还包括:
将满足如下公式(3)的主导新能源场站集D中新能源场站分别 用一个主导新能源等值场站替代,更新主导新能源场站集D,将主导 新能源等值场站代表的主导新能源场站对电网暂态稳定综合影响因 子的平均值作为其对电网暂态稳定的综合影响因子,
式中,ε为设定参数,大于0;λi为预想故障下的主导新能源场站 集D中主导新能源场站i对电网暂态稳定的综合影响因子;λj为预想 故障下的主导新能源场站集D中主导新能源场站j对电网暂态稳定的 综合影响因子;λd为预想故障下的主导新能源场站集D中主导新能 源场站d对电网暂态稳定的综合影响因子。
作为一种较佳的实施例,所述步骤SS3中,确定预想故障下主导 新能源场站集D中各新能源场站的有功采样区间,生成主导新能源 场站集D中各新能源场站有功分档区间组合集C,包括以下步骤:
根据设定的置信水平参数和t时刻主导新能源场站集D中各新能 源场站有功概率分布密度函数,确定t时刻D中各新能源场站的有功 置信区间;
针对主导新能源场站集D中对电网暂态稳定的综合影响因子大 于0的新能源场站,将t时刻有功实时计划值作为其有功采样区间下 限、t时刻有功置信区间上限作为其有功采样区间上限;针对主导新 能源场站集D中对电网暂态稳定的综合影响因子小于0的新能源场 站,将t时刻有功实时计划值作为其有功采样区间上限、t时刻有功置 信区间下限作为其有功采样区间下限;
通过公式(4)计算主导新能源场站集D中各新能源场站有功采 样区间的分档数,分别按其分档数将有功采样区间划分为区间长度相 等的有功分档区间;
式中,kd为t时刻主导新能源场站集D中新能源场站d有功采样 区间的分档数,δp为设定的功率分档参数,Pd.u、Pd.d分别为t时刻主 导新能源场站集D中新能源场站d有功采样区间上限和下限,λd为 主导新能源场站集D中新能源场站d对电网暂态稳定的综合影响因子;
采用枚举法分别从主导新能源场站集D中各新能源场站的有功 分档区间中,任选一个有功分档区间进行穷尽式组合,生成的组合组 成的集合作为新能源场站有功分档区间组合集C。
作为一种较佳的实施例,所述步骤SS4中,计算新能源场站有功 分档区间组合集C中各组合对应的电网运行状态下预想故障的暂态 稳定裕度,得到预想故障下暂态稳定裕度越限组合集U,具体包括如 下步骤:
采用如下公式(5)计算新能源场站有功分档区间组合集C中各 组合的排序指标,按照新能源场站有功分档区间组合集C中各组合 排序指标由大到小的顺序,针对对应的电网运行状态进行预想故障下 的暂态稳定裕度计算,直至序号连续的多个组合下暂态稳定裕度都不 小于暂态稳定越限的裕度定值,且这些组合中排序指标最大值与最小 值之比大于设定参数,所述设定参数大于1,中止余下组合对应的电 网运行状态生成和暂态稳定裕度计算,将暂态稳定裕度小于暂态稳定 越限的裕度定值的新能源场站有功分档区间组合集C中组合组成的 集合作为U;
sc=∑d∈D[λd(Pd.c-Pd.0)],c∈C (5)
式中,sc为新能源场站有功分档区间组合集C中组合c的排序指 标,Pd.c为C中组合c对应的主导新能源场站集D中新能源场站d有 功分档区间的有功期望值,Pd.0为t时刻主导新能源场站集D中新能 源场站d有功实时计划值。
作为一种较佳的实施例,所述步骤SS4具体还包括:基于电网预 想运行状态S,以组合中各新能源场站有功分档区间的有功期望值作 为新能源场站的有功,通过潮流计算得到新能源场站有功分档区间组 合集C中各组合对应的电网运行状态。
作为一种较佳的实施例,所述步骤SS5中,新能源发电不确定性 的暂态稳定越限概率的计算公式为:
式中,pol为考虑t时刻新能源发电不确定性的暂态稳定越限概率, pf为预想故障集F中预想故障f在t时刻发生的概率,Uf为预想故障 集F中预想故障f下暂态稳定裕度越限组合集,为Uf中组合uf对 应的组合概率。
本发明还提出新能源发电不确定性的暂态稳定越限概率计算装 置,包括:
暂态稳定量化评估模块,用于执行:针对设定时刻t的电网预想 运行状态S,进行预想故障集F中各预想故障下暂态稳定量化评估, 得到预想故障集F中各预想故障的暂态稳定关键模式集M及暂态稳 定关键模式集M中各关键模式对应的暂态稳定裕度,和新能源场站 集N中各新能源场站的暂态稳定参与因子;
主导新能源场站确定模块,用于执行:针对预想故障集F中各预 想故障,分别计算预想故障下的新能源场站集N中各新能源场站对 电网暂态稳定的综合影响因子,生成预想故障关联的主导新能源场站 集D;
新能源场站有功分档区间组合模块,用于执行:针对预想故障集 F中各预想故障,分别根据t时刻主导新能源场站集D中各新能源场 站有功概率分布密度函数、有功实时计划值及各新能源场站对电网暂 态稳定的综合影响因子,确定预想故障下D中各新能源场站的有功 采样区间,生成主导新能源场站集D中各新能源场站有功分档区间 组合集C;
暂态稳定越限组合确定模块,用于执行:针对预想故障集F中各 预想故障,分别计算有功分档区间组合集C中各组合对应的电网运 行状态下预想故障的暂态稳定裕度,得到预想故障下暂态稳定裕度越 限组合集U;
暂态稳定越限概率计算模块,用于执行:根据t时刻预想故障集 F中各预想故障发生概率和预想故障下暂态稳定裕度越限组合集U中 各组合的概率,计算考虑新能源发电不确定性的暂态稳定越限概率。
本发明所达到的有益效果:第一,本发明提出新能源发电不确定 性的暂态稳定越限概率计算方法,通过进行预想故障集F中各预想故 障下暂态稳定量化评估,得到预想故障集F中各预想故障的暂态稳定 关键模式集M及暂态稳定关键模式集M中各关键模式对应的暂态稳 定裕度,和新能源场站集N中各新能源场站的暂态稳定参与因子, 进而计算预想故障下的新能源场站集N中各新能源场站对电网暂态 稳定的综合影响因子,生成预想故障关联的主导新能源场站集D,确 定预想故障下D中各新能源场站的有功采样区间,生成主导新能源 场站集D中各新能源场站有功分档区间组合集C,再计算有功分档区 间组合集C中各组合对应的电网运行状态下预想故障的暂态稳定裕 度,得到预想故障下暂态稳定裕度越限组合集U,大幅削减参与计算 的新能源场站数,根据t时刻预想故障集F中各预想故障发生概率和 预想故障下暂态稳定裕度越限组合集U中各组合的概率,计算考虑 新能源发电不确定性的暂态稳定越限概率,实现了计算量的大幅缩减; 第二,本发明提出新能源发电不确定性的暂态稳定越限概率计算装置, 通过执行暂态稳定量化评估模块、主导新能源场站确定模块、新能源 场站有功分档区间组合模块、暂态稳定越限组合确定模块、暂态稳定 越限概率计算模块,根据新能源场站有功置信区间、有功实时计划值 和综合影响因子符号,确定用于计算的新能源场站有功采样区间,缩 小新能源场站有功取值范围,结合综合影响因子对采样区间进行分档, 有效降低新能源场站有功分档区间组合数,实现了计算场景的大幅缩 减;第三,本发明可满足考虑新能源发电不确定性的暂态稳定越限概 率在线计算的实时性要求。
附图说明
图1是本发明的新能源发电不确定性的暂态稳定越限概率计算 方法的流程图。
具体实施方式
下面结合附图对本发明作进一步描述。以下实施例仅用于更加清 楚地说明本发明的技术方案,而不能以此来限制本发明的保护范围。
实施例1:如图1所示,本发明提出新能源发电不确定性的暂态 稳定越限概率计算方法,包括以下步骤:
步骤1:针对设定时刻t电网预想运行状态S进行预想故障集F 中各预想故障下暂态稳定量化评估,得到F中各预想故障的暂态稳定 关键模式集M及M中各关键模式对应的暂态稳定裕度和新能源场站 集N中各新能源场站的暂态稳定参与因子;
其中,M为暂态稳定裕度与所有暂态稳定模式中暂态稳定裕度最 小值之差小于设定值(通常设置为0.2)的暂态稳定模式组成的集合。
参与因子大于0,表示新能源场站有功增加,暂态稳定裕度下降; 参与因子小于0,表示新能源场站有功减少,暂态稳定裕度下降。
步骤2:针对F中各预想故障,分别计算预想故障下N中各新能 源场站对电网暂态稳定的综合影响因子,生成预想故障关联的主导新 能源场站集D;
其中,预想故障下N中各新能源场站对电网暂态稳定综合影响 因子的计算公式为:
λn=∑m∈M[λm.n(1-ηm)α],n∈N (1)
式中,λn为预想故障下N中新能源场站n对电网暂态稳定的综 合影响因子,λm.n为该预想故障的M中关键模式m下N中新能源场 站n的暂态稳定参与因子,ηm为该预想故障的M中关键模式m的暂 态稳定裕度,取值范围为[-1,1],α为设定参数,大于0,通常设置为 2。
生成预想故障关联的主导新能源场站集D,方法为:
将预想故障下不满足公式(2)的新能源场站从新能源场站集N 中剔除后,作为与该预想故障关联的D;
式中,β为设定参数,大于0,小于1,通常设置为0.5。
将满足公式(3)的D中新能源场站分别用一个主导新能源等值 场站替代,更新D,将主导新能源等值场站代表的主导新能源场站对 电网暂态稳定综合影响因子的平均值作为其对电网暂态稳定的综合 影响因子,
式中,ε为设定参数,大于0,通常设置为0.05。
步骤3:针对F中各预想故障,分别根据t时刻D中各新能源场 站有功概率分布密度函数、有功实时计划值及其对电网暂态稳定的综 合影响因子,确定预想故障下D中各新能源场站的有功采样区间, 生成D中各新能源场站有功分档区间组合集C;
对于D中的新能源等值场站,根据新能源等值场站代表的各新 能源场站有功概率分布密度函数,生成新能源等值场站的有功概率分 布密度函数,将新能源等值场站代表的各新能源场站有功实时计划值 之和,作为新能源等值场站的有功实时计划值;
其中,确定预想故障下D中各新能源场站的有功采样区间,生 成D中各新能源场站有功分档区间组合集C,包括以下步骤:
根据设定的置信水平参数(通常设置为0.95)和t时刻D中各新 能源场站有功概率分布密度函数,确定t时刻D中各新能源场站的有 功置信区间;
针对D中对电网暂态稳定的综合影响因子大于0的新能源场站, 将t时刻有功实时计划值作为其有功采样区间下限、t时刻有功置信区 间上限作为其有功采样区间上限;针对D中对电网暂态稳定的综合 影响因子小于0的新能源场站,将t时刻有功实时计划值作为其有功 采样区间上限、t时刻有功置信区间下限作为其有功采样区间下限;
通过公式(4)计算D中各新能源场站有功采样区间的分档数, 分别按其分档数将有功采样区间划分为区间长度相等的有功分档区 间;
式中,kd为t时刻D中新能源场站d有功采样区间的分档数,δp为 设定的功率分档参数,Pd.u、Pd.d分别为t时刻D中新能源场站d有功 采样区间上限和下限,λd为D中新能源场站d对电网暂态稳定的综 合影响因子。
采用枚举法分别从D中各新能源场站的有功分档区间中任选一 个有功分档区间进行穷尽式组合生成的组合组成的集合,作为C。
步骤4:针对F中各预想故障,分别计算C中各组合对应的电网 运行状态下预想故障的暂态稳定裕度,得到预想故障下暂态稳定裕度 越限组合集U;
其中,计算C中各组合对应的电网运行状态下预想故障的暂态 稳定裕度,得到预想故障下暂态稳定裕度越限组合集U,方法为:
采用公式(5)计算C中各组合的排序指标,按照C中各组合排 序指标由大到小的顺序,针对其对应的电网运行状态进行预想故障下 的暂态稳定裕度计算,直至序号连续的多个组合下暂态稳定裕度都不 小于暂态稳定越限的裕度定值且这些组合中排序指标最大值与最小 值之比大于设定参数,该参数大于1(通常设置为1.2),中止余下组 合对应的电网运行状态生成和暂态稳定裕度计算,将暂态稳定裕度小 于暂态稳定越限的裕度定值的C中组合组成的集合作为U;
sc=∑d∈D[λd(Pd.c-Pd.0)],c∈C (5)
式中,sc为C中组合c的排序指标,Pd.c为C中组合c对应的D 中新能源场站d有功分档区间的有功期望值,Pd.0为t时刻D中新能 源场站d有功实时计划值。
基于S,以组合中各新能源场站有功分档区间的有功期望值作为 新能源场站的有功,通过潮流计算得到C中各组合对应的电网运行 状态。
在潮流计算中考虑t时刻电网中所有参与一次调频的新能源场站、 非新能源发电厂频率控制模型参数,以及负荷频率特性模型参数和直 流***频率调制策略;
其中,按照C中各组合排序指标由大到小的顺序,针对其对应 的电网运行状态进行预想故障下的暂态稳定裕度计算,可以按组合的 排序,针对多个序号连续的组合,采用多组合并行处理的方式提高计 算速度。
步骤5:根据t时刻F中各预想故障发生概率和预想故障下U中 各组合的概率,计算考虑新能源发电不确定性的暂态稳定越限概率。
其中,新能源发电不确定性的暂态稳定越限概率的计算公式为:
组合概率为组合内D中各新能源场站有功分档区间所对应累积 概率的乘积。
实施例2:本发明提出新能源发电不确定性的暂态稳定越限概率 计算装置,包括:
暂态稳定量化评估模块,用于执行:针对设定时刻t电网预想运 行状态S进行预想故障集F中各预想故障下暂态稳定量化评估,得到 F中各预想故障的暂态稳定关键模式集M及M中各关键模式对应的 暂态稳定裕度和新能源场站集N中各新能源场站的暂态稳定参与因 子;
主导新能源场站确定模块,用于执行:针对F中各预想故障,分 别计算预想故障下N中各新能源场站对电网暂态稳定的综合影响因 子,生成预想故障关联的主导新能源场站集D;
新能源场站有功分档区间组合模块,用于执行:针对F中各预想 故障,分别根据t时刻D中各新能源场站有功概率分布密度函数、有 功实时计划值及其对电网暂态稳定的综合影响因子,确定预想故障下 D中各新能源场站的有功采样区间,生成D中各新能源场站有功分 档区间组合集C;
暂态稳定越限组合确定模块,用于执行:针对F中各预想故障, 分别计算C中各组合对应的电网运行状态下预想故障的暂态稳定裕 度,得到预想故障下暂态稳定裕度越限组合集U;
暂态稳定越限概率计算模块,用于执行:根据t时刻F中各预想 故障发生概率和预想故障下U中各组合的概率,计算考虑新能源发 电不确定性的暂态稳定越限概率。
本申请是参照根据本申请实施例的方法、设备(***)和计算机 程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指 令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和 /或方框图中的流程和/或方框的结合。可提供这些计算机程序指令 到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设 备的处理器以产生一个机器,使得通过计算机或其他可编程数据处理 设备的处理器执行的指令产生用于实现在流程图一个流程或多个流 程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数 据处理设备以特定方式工作的计算机可读存储器中,使得存储在该计 算机可读存储器中的指令产生包括指令装置的制造品,该指令装置实 现在流程图一个流程或多个流程和/或方框图一个方框或多个方框 中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理 设备上,使得在计算机或其他可编程设备上执行一系列操作步骤以产 生计算机实现的处理,从而在计算机或其他可编程设备上执行的指令 提供用于实现在流程图一个流程或多个流程和/或方框图一个方框 或多个方框中指定的功能的步骤。
以上仅为本发明的实施例而已,并不用于限制本发明,凡在本发 明的精神和原则之内,所做的任何修改、等同替换、改进等,均包含 在申请待批的本发明的权利要求范围之内。
Claims (10)
1.新能源发电不确定性的暂态稳定越限概率计算方法,其特征在于,所述方法包括:
步骤SS1:针对设定时刻t的电网预想运行状态S,进行预想故障集F中各预想故障下暂态稳定量化评估,得到预想故障集F中各预想故障的暂态稳定关键模式集M及暂态稳定关键模式集M中各关键模式对应的暂态稳定裕度,和新能源场站集N中各新能源场站的暂态稳定参与因子;
步骤SS2:针对预想故障集F中各预想故障,分别计算预想故障下的新能源场站集N中各新能源场站对电网暂态稳定的综合影响因子,生成预想故障关联的主导新能源场站集D;
步骤SS3:针对预想故障集F中各预想故障,分别根据t时刻主导新能源场站集D中各新能源场站有功概率分布密度函数、有功实时计划值及各新能源场站对电网暂态稳定的综合影响因子,确定预想故障下D中各新能源场站的有功采样区间,生成主导新能源场站集D中各新能源场站有功分档区间组合集C;
步骤SS4:针对预想故障集F中各预想故障,分别计算有功分档区间组合集C中各组合对应的电网运行状态下预想故障的暂态稳定裕度,得到预想故障下暂态稳定裕度越限组合集U;
步骤SS5:根据t时刻预想故障集F中各预想故障发生概率和预想故障下暂态稳定裕度越限组合集U中各组合的概率,计算考虑新能源发电不确定性的暂态稳定越限概率。
2.根据权利要求1所述的新能源发电不确定性的暂态稳定越限概率计算方法,其特征在于,所述步骤SS1具体包括:预想故障集F中暂态稳定关键模式集M为暂态稳定裕度与所有暂态稳定模式中暂态稳定裕度最小值之差小于设定值的暂态稳定模式组成的集合。
3.根据权利要求1所述的新能源发电不确定性的暂态稳定越限概率计算方法,其特征在于,所述步骤SS2具体包括:
预想故障下的新能源场站集N中各新能源场站对电网暂态稳定综合影响因子的计算公式为:
λn=∑m∈M[λm.n(1-ηm)α],n∈N (1)
式中,λn为预想故障下的新能源场站集N中新能源场站n对电网暂态稳定的综合影响因子,λm.n为该预想故障的暂态稳定关键模式集M中关键模式m下新能源场站集N中新能源场站n的暂态稳定参与因子,ηm为该预想故障的暂态稳定关键模式集M中关键模式m的暂态稳定裕度,取值范围为[-1,1],α为设定参数,大于0。
5.根据权利要求4所述的新能源发电不确定性的暂态稳定越限概率计算方法,其特征在于,所述步骤SS2还包括:
将满足如下公式(3)的主导新能源场站集D中新能源场站分别用一个主导新能源等值场站替代,更新主导新能源场站集D,将主导新能源等值场站代表的主导新能源场站对电网暂态稳定综合影响因子的平均值作为其对电网暂态稳定的综合影响因子,
式中,ε为设定参数,大于0;λi为预想故障下的主导新能源场站集D中主导新能源场站i对电网暂态稳定的综合影响因子;λj为预想故障下的主导新能源场站集D中主导新能源场站j对电网暂态稳定的综合影响因子;λd为预想故障下的主导新能源场站集D中主导新能源场站d对电网暂态稳定的综合影响因子。
6.根据权利要求1所述的新能源发电不确定性的暂态稳定越限概率计算方法,其特征在于,所述步骤SS3中,确定预想故障下主导新能源场站集D中各新能源场站的有功采样区间,生成主导新能源场站集D中各新能源场站有功分档区间组合集C,包括以下步骤:
根据设定的置信水平参数和t时刻主导新能源场站集D中各新能源场站有功概率分布密度函数,确定t时刻D中各新能源场站的有功置信区间;
针对主导新能源场站集D中对电网暂态稳定的综合影响因子大于0的新能源场站,将t时刻有功实时计划值作为其有功采样区间下限、t时刻有功置信区间上限作为其有功采样区间上限;针对主导新能源场站集D中对电网暂态稳定的综合影响因子小于0的新能源场站,将t时刻有功实时计划值作为其有功采样区间上限、t时刻有功置信区间下限作为其有功采样区间下限;
通过公式(4)计算主导新能源场站集D中各新能源场站有功采样区间的分档数,分别按其分档数将有功采样区间划分为区间长度相等的有功分档区间;
式中,kd为t时刻主导新能源场站集D中新能源场站d有功采样区间的分档数,δp为设定的功率分档参数,Pd.u、Pd.d分别为t时刻主导新能源场站集D中新能源场站d有功采样区间上限和下限,λd为主导新能源场站集D中新能源场站d对电网暂态稳定的综合影响因子;
采用枚举法分别从主导新能源场站集D中各新能源场站的有功分档区间中,任选一个有功分档区间进行穷尽式组合,生成的组合组成的集合作为新能源场站有功分档区间组合集C。
7.根据权利要求1所述的新能源发电不确定性的暂态稳定越限概率计算方法,其特征在于,所述步骤SS4中,计算新能源场站有功分档区间组合集C中各组合对应的电网运行状态下预想故障的暂态稳定裕度,得到预想故障下暂态稳定裕度越限组合集U,具体包括如下步骤:
采用如下公式(5)计算新能源场站有功分档区间组合集C中各组合的排序指标,按照新能源场站有功分档区间组合集C中各组合排序指标由大到小的顺序,针对对应的电网运行状态进行预想故障下的暂态稳定裕度计算,直至序号连续的多个组合下暂态稳定裕度都不小于暂态稳定越限的裕度定值,且这些组合中排序指标最大值与最小值之比大于设定参数,所述设定参数大于1,中止余下组合对应的电网运行状态生成和暂态稳定裕度计算,将暂态稳定裕度小于暂态稳定越限的裕度定值的新能源场站有功分档区间组合集C中组合组成的集合作为U;
sc=∑d∈D[λd(Pd.c-Pd.0)],c∈C (5)
式中,sc为新能源场站有功分档区间组合集C中组合c的排序指标,Pd.c为C中组合c对应的主导新能源场站集D中新能源场站d有功分档区间的有功期望值,Pd.0为t时刻主导新能源场站集D中新能源场站d有功实时计划值。
8.根据权利要求7所述的新能源发电不确定性的暂态稳定越限概率计算方法,其特征在于,所述步骤SS4具体还包括:基于电网预想运行状态S,以组合中各新能源场站有功分档区间的有功期望值作为新能源场站的有功,通过潮流计算得到新能源场站有功分档区间组合集C中各组合对应的电网运行状态。
10.新能源发电不确定性的暂态稳定越限概率计算装置,其特征在于,包括:
暂态稳定量化评估模块,用于执行:针对设定时刻t的电网预想运行状态S,进行预想故障集F中各预想故障下暂态稳定量化评估,得到预想故障集F中各预想故障的暂态稳定关键模式集M及暂态稳定关键模式集M中各关键模式对应的暂态稳定裕度,和新能源场站集N中各新能源场站的暂态稳定参与因子;
主导新能源场站确定模块,用于执行:针对预想故障集F中各预想故障,分别计算预想故障下的新能源场站集N中各新能源场站对电网暂态稳定的综合影响因子,生成预想故障关联的主导新能源场站集D;
新能源场站有功分档区间组合模块,用于执行:针对预想故障集F中各预想故障,分别根据t时刻主导新能源场站集D中各新能源场站有功概率分布密度函数、有功实时计划值及各新能源场站对电网暂态稳定的综合影响因子,确定预想故障下D中各新能源场站的有功采样区间,生成主导新能源场站集D中各新能源场站有功分档区间组合集C;
暂态稳定越限组合确定模块,用于执行:针对预想故障集F中各预想故障,分别计算有功分档区间组合集C中各组合对应的电网运行状态下预想故障的暂态稳定裕度,得到预想故障下暂态稳定裕度越限组合集U;
暂态稳定越限概率计算模块,用于执行:根据t时刻预想故障集F中各预想故障发生概率和预想故障下暂态稳定裕度越限组合集U中各组合的概率,计算考虑新能源发电不确定性的暂态稳定越限概率。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011053005.2A CN112383045B (zh) | 2020-09-29 | 2020-09-29 | 新能源发电不确定性的暂态稳定越限概率计算方法和装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011053005.2A CN112383045B (zh) | 2020-09-29 | 2020-09-29 | 新能源发电不确定性的暂态稳定越限概率计算方法和装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112383045A true CN112383045A (zh) | 2021-02-19 |
CN112383045B CN112383045B (zh) | 2022-08-30 |
Family
ID=74580884
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011053005.2A Active CN112383045B (zh) | 2020-09-29 | 2020-09-29 | 新能源发电不确定性的暂态稳定越限概率计算方法和装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112383045B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113469410A (zh) * | 2021-05-26 | 2021-10-01 | 国电南瑞科技股份有限公司 | 考虑新能源不确定性的电网最大消纳能力在线计算方法及装置 |
CN113765153A (zh) * | 2021-08-24 | 2021-12-07 | 国网电力科学研究院有限公司 | 一种基于风险评估的新能源极限接入能力计算方法及*** |
CN115498641A (zh) * | 2022-11-18 | 2022-12-20 | 中国电力科学研究院有限公司 | 一种基于新能源多场站短路比的决策控制方法及*** |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110556843A (zh) * | 2019-09-17 | 2019-12-10 | 国电南瑞科技股份有限公司 | 一种新能源场站的暂态功角稳定性能指标设计方法及*** |
CN111130100A (zh) * | 2020-01-03 | 2020-05-08 | 国电南瑞科技股份有限公司 | 考虑新能源不确定性的日内计划安全稳定校核方法及*** |
-
2020
- 2020-09-29 CN CN202011053005.2A patent/CN112383045B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110556843A (zh) * | 2019-09-17 | 2019-12-10 | 国电南瑞科技股份有限公司 | 一种新能源场站的暂态功角稳定性能指标设计方法及*** |
CN111130100A (zh) * | 2020-01-03 | 2020-05-08 | 国电南瑞科技股份有限公司 | 考虑新能源不确定性的日内计划安全稳定校核方法及*** |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113469410A (zh) * | 2021-05-26 | 2021-10-01 | 国电南瑞科技股份有限公司 | 考虑新能源不确定性的电网最大消纳能力在线计算方法及装置 |
CN113469410B (zh) * | 2021-05-26 | 2023-10-31 | 国电南瑞科技股份有限公司 | 考虑新能源不确定性的电网最大消纳能力在线计算方法及装置 |
CN113765153A (zh) * | 2021-08-24 | 2021-12-07 | 国网电力科学研究院有限公司 | 一种基于风险评估的新能源极限接入能力计算方法及*** |
CN113765153B (zh) * | 2021-08-24 | 2023-09-12 | 国网电力科学研究院有限公司 | 一种基于风险评估的新能源极限接入能力计算方法及*** |
CN115498641A (zh) * | 2022-11-18 | 2022-12-20 | 中国电力科学研究院有限公司 | 一种基于新能源多场站短路比的决策控制方法及*** |
Also Published As
Publication number | Publication date |
---|---|
CN112383045B (zh) | 2022-08-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112383045B (zh) | 新能源发电不确定性的暂态稳定越限概率计算方法和装置 | |
WO2019184286A1 (zh) | 一种机组恢复在线动态决策方法及*** | |
CN107256449B (zh) | 一种智能变电站继电保护装置状态评价与评估方法 | |
CN113077075B (zh) | 新能源不确定性电力***安全风险预防控制方法及装置 | |
CN102117450A (zh) | 基于经验反馈的核电设备预防性更换周期的优化方法 | |
CN109993665B (zh) | 电力***在线安全稳定评估方法、装置及*** | |
CN111680712B (zh) | 基于日内相似时刻的变压器油温预测方法、装置及*** | |
CN110556823B (zh) | 基于模型降维的安全约束机组组合快速计算方法及*** | |
CN106056305B (zh) | 一种基于状态聚类的发电***可靠性快速评估方法 | |
CN112734141B (zh) | 一种多元化负荷区间预测方法及装置 | |
CN117909864A (zh) | 一种电力故障预测***及方法 | |
CN117557047A (zh) | 一种基于深度强化学习的配电设备运维优化方法及*** | |
CN106022594B (zh) | 基于极值理论的电力***在线运行安全风险评估方法 | |
JPH11346438A (ja) | 電力需要量予測方法 | |
CN115829122A (zh) | 一种电力客服工单投诉预警方法及*** | |
CN110797865A (zh) | 一种数据驱动的供电薄弱环节识别方法及*** | |
CN112467806B (zh) | 一种基于dmn电网调控规则确定电网运行状态的方法及*** | |
CN107844909B (zh) | 核电厂日常生产风险管理的滚动裕度风险控制方法及*** | |
CN114266427A (zh) | 一种基于拓扑分析的新设备投运启动方案生成方法 | |
CN111697562B (zh) | 一种基于svm的暂态稳定预防控制策略的生成方法及装置 | |
CN113112296A (zh) | 一种集中式现货市场的典型场景库构建方法及装置 | |
CN112241587A (zh) | 配电线路风险评估模型构建方法、风险评估方法及*** | |
CN111428756A (zh) | 基于时间序列信息熵的规划数据融合实时态方法及装置 | |
CN111952959A (zh) | 压缩电网过程模拟时间的方法、装置及存储介质 | |
CN113822470B (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 |