CN115411735A - 一种配电网预期未供应电量计算方法、***及设备 - Google Patents
一种配电网预期未供应电量计算方法、***及设备 Download PDFInfo
- Publication number
- CN115411735A CN115411735A CN202211119436.3A CN202211119436A CN115411735A CN 115411735 A CN115411735 A CN 115411735A CN 202211119436 A CN202211119436 A CN 202211119436A CN 115411735 A CN115411735 A CN 115411735A
- Authority
- CN
- China
- Prior art keywords
- power
- load
- distribution network
- model
- 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.)
- Pending
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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/18—Network design, e.g. design based on topological or interconnect aspects of utility systems, piping, heating ventilation air conditioning [HVAC] or cabling
-
- 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/003—Load forecast, e.g. methods or systems for forecasting future load demand
-
- 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/388—Islanding, i.e. disconnection of local power supply from the network
-
- 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/46—Controlling of the sharing of output between the generators, converters, or transformers
- H02J3/466—Scheduling the operation of the generators, e.g. connecting or disconnecting generators to meet a given demand
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/04—Constraint-based CAD
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/04—Power grid distribution 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)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Physics (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Power Engineering (AREA)
- General Engineering & Computer Science (AREA)
- Geometry (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- Evolutionary Computation (AREA)
- Computing Systems (AREA)
- Computer Hardware Design (AREA)
- Computer Networks & Wireless Communication (AREA)
- Operations Research (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本申请涉及一种配电网预期未供应电量计算方法、***及设备,该方法包括获取含有分布式电源的配电网的网络拓扑模型并对其拓扑连接关系构成节点关联矩阵;获取网络拓扑模型中各个母线节点的负荷需求和各个分布式发电装置的安装容量并构建虚拟潮流平衡模型;根据故障条件和虚拟潮流平衡模型构建配电网的实际缺供电负荷模型以及计算实际缺供电负荷;获取配电网中电气元件的故障参数,根据实际缺供电负荷和故障参数构建预期未供应电量计算模型;根据虚拟潮流平衡模型、实际缺供电负荷模型和预期未供应电量计算模型构建线性规划优化模型并对其计算配电网的预期未供应电量。该方法计算配电网的预期未供应电量的速度快、计算准确,效率高。
Description
技术领域
本申请涉及配电网技术领域,尤其涉及一种配电网预期未供应电量计算方法、***及设备。
背景技术
配电网作为电力***中和用户直接相连的一部分,其供电可靠性直接影响着用户的用电体验。针对配电网的可靠性通常通过使用量化的可靠性指标进行定量评估。根据可靠性评估结果不仅能够找出配电网中的存在的隐患,并有针对性地排除隐患提高***可靠性,还能够对评估***进行可靠性的成本效益分析。配电网可靠性指标包括各负荷点、变电站以及配电网的缺电时间期望值、电量不足期望值、配电网的平均***断电频率、年均断电时间、用户平均断电频率、平均供电可用率和平均电量不足等。
现有对配电网的可靠性评估算法大多是像解析法或蒙特卡洛算法这类传统的可靠性评估算法,通常包含复杂的拓扑搜索和逻辑判断,难以运用数学解析式显式表达,因此求解效率不高,也难以与数学规划模型进行有效结合。
发明内容
本申请实施例提供了一种配电网预期未供应电量计算方法、***及设备,用于解决现有对配电网可靠性指标求解方式难以与数学规划模型结合,存在效率不高的技术问题。
为了实现上述目的,本申请实施例提供如下技术方案:
一种配电网预期未供应电量计算方法,包括以下步骤:
获取含有分布式电源的配电网的网络拓扑模型,根据所述网络拓扑模型的拓扑连接关系构成节点关联矩阵;
获取所述网络拓扑模型中各个母线节点的负荷需求和各个分布式发电装置的安装容量,根据所有的所述负荷需求和所述安装容量构建虚拟潮流平衡模型;
根据故障条件和所述虚拟潮流平衡模型,构建配电网发生线路元件故障时各支路下游的实际缺供电负荷模型,基于所述实际缺供电负荷模型获得配电网发生线路元件故障时各支路下游的实际缺供电负荷;
获取配电网中电气元件的故障参数,根据所述实际缺供电负荷和所述故障参数构建预期未供应电量计算模型;
根据所述虚拟潮流平衡模型、所述实际缺供电负荷模型和所述预期未供应电量计算模型构建线性规划优化模型,对所述线性规划优化模型进行优化求解,得到配电网的预期未供应电量。
优选地,对所述线性规划优化模型进行优化求解,得到配电网的预期未供应电量包括:采用Gurobi自带的分支定界法和目标条件对所述线性规划优化模型进行优化求解,得到配电网的预期未供应电量;其中,所述目标条件为选取所述实际缺供电负荷的最小值和Gurobi中松弛解与收敛解之间的距离为0.1%。
优选地,所述虚拟潮流平衡模型包括:负荷需求虚拟潮流平衡方程和容量虚拟潮流平衡方程,所述负荷需求虚拟潮流平衡方程为:
所述容量虚拟潮流平衡方程为:
式中,Ω为配电网的支路集合,Db为支路b下游的总最大负荷需求,为支路b在正方向上的虚拟负荷需求潮流,为支路b在负方向上的虚拟负荷需求潮流,MD为负荷需求虚拟潮流平衡方程取足够大的正整数,为第一辅助二值变量,Ai,b为节点关联矩阵中第i行第b列对应的元素数值,ψ为配电网的节点i集合,Di为负荷节点i的负荷需求;Ds为电源节点s的负荷需求;Gs为电源节点s下游的总负荷需求,CDGb为支路b下游的所有分布式发电装置计及换流效率的供电容量,为支路b在正方向上的虚拟负荷安装容量潮流,为支路b在负方向上的虚拟负荷安装容量潮流,为第二辅助二值变量,MDG为容量虚拟潮流平衡方程取足够大的正整数;为负荷节点i所安装分布式发电装置的交流安装容量,为负荷节点i所安装分布式发电装置的直流安装容量,为电源节点i所安装分布式发电装置的交流安装容量,为电源节点i所安装分布式发电装置的直流安装容量,ηivt为配电网中逆变器的换流效率,为电源节点s下游的计及换流效率分布式发电装置的安装总容量。
优选地,所述实际缺供电负荷模型采用损失负荷量方程表示,所述损失负荷量方程为:
式中,Ω为配电网的支路集合,Db为支路b下游的总最大负荷需求,CDGb为支路b下游的所有分布式发电装置计及换流效率的供电容量,Hb为配电网支路b故障而导致的下游负荷功率缺失量,MCP为损失负荷量方程取的足够大的正整数,为第三辅助二值变量。
优选地,所述预期未供应电量计算模型包括:用于计算配电网预期未供应电量的电量计算方程,所述电量计算方程为:
式中,EEENS为配电网的预期未供应电量,λb、λivt、λDCB分别为配电网支路b上线路、逆变器和直流断路器的故障率,τline、τivt、τDCB分别为配电网支路b上线路、逆变器和直流断路器的故障修复时间,Ω为配电网的支路集合,为负荷节点i所安装分布式发电装置的直流安装容量,Hb为配电网支路b故障而导致的下游负荷功率缺失量,ψ为配电网的节点i集合。
本申请还提供一种配电网预期未供应电量计算***,包括矩阵构建模块、潮流模型构建模块、供电负荷计算模块、电量模型构建模块和供应量计算模块;
所述矩阵构建模块,用于获取含有分布式电源的配电网的网络拓扑模型,根据所述网络拓扑模型的拓扑连接关系构成节点关联矩阵;
所述潮流模型构建模块,用于获取所述网络拓扑模型中各个母线节点的负荷需求和各个分布式发电装置的安装容量,根据所有的所述负荷需求和所述安装容量构建虚拟潮流平衡模型;
所述供电负荷计算模块,用于根据故障条件和所述虚拟潮流平衡模型,构建配电网发生线路元件故障时各支路下游的实际缺供电负荷模型,基于所述实际缺供电负荷模型获得配电网发生线路元件故障时各支路下游的实际缺供电负荷;
所述电量模型构建模块,用于获取配电网中电气元件的故障参数,根据所述实际缺供电负荷和所述故障参数构建预期未供应电量计算模型;
所述供应量计算模块,用于根据所述虚拟潮流平衡模型、所述实际缺供电负荷模型和所述预期未供应电量计算模型构建线性规划优化模型,对所述线性规划优化模型进行优化求解,得到配电网的预期未供应电量。
优选地,所述供应量计算模块还用于采用Gurobi自带的分支定界法和目标条件对所述线性规划优化模型进行优化求解,得到配电网的预期未供应电量;其中,所述目标条件为选取所述实际缺供电负荷的最小值和Gurobi中松弛解与收敛解之间的距离为0.1%。
优选地,所述虚拟潮流平衡模型包括:负荷需求虚拟潮流平衡方程和容量虚拟潮流平衡方程,所述负荷需求虚拟潮流平衡方程为:
所述容量虚拟潮流平衡方程为:
式中,Ω为配电网的支路集合,Db为支路b下游的总最大负荷需求,为支路b在正方向上的虚拟负荷需求潮流,为支路b在负方向上的虚拟负荷需求潮流,MD为负荷需求虚拟潮流平衡方程取足够大的正整数,为第一辅助二值变量,Ai,b为节点关联矩阵中第i行第b列对应的元素数值,ψ为配电网的节点i集合,Di为负荷节点i的负荷需求;Ds为电源节点s的负荷需求;Gs为电源节点s下游的总负荷需求,CDGb为支路b下游的所有分布式发电装置计及换流效率的供电容量,为支路b在正方向上的虚拟负荷安装容量潮流,为支路b在负方向上的虚拟负荷安装容量潮流,为第二辅助二值变量,MDG为容量虚拟潮流平衡方程取足够大的正整数;为负荷节点i所安装分布式发电装置的交流安装容量,为负荷节点i所安装分布式发电装置的直流安装容量,为电源节点i所安装分布式发电装置的交流安装容量,为电源节点i所安装分布式发电装置的直流安装容量,ηivt为配电网中逆变器的换流效率,为电源节点s下游的计及换流效率分布式发电装置的安装总容量。
优选地,所述实际缺供电负荷模型采用损失负荷量方程表示,所述损失负荷量方程为:
式中,Ω为配电网的支路集合,Db为支路b下游的总最大负荷需求,CDGb为支路b下游的所有分布式发电装置计及换流效率的供电容量,Hb为配电网支路b故障而导致的下游负荷功率缺失量,MCP为损失负荷量方程取的足够大的正整数,为第三辅助二值变量;
所述预期未供应电量计算模型包括:用于计算配电网预期未供应电量的电量计算方程,所述电量计算方程为:
式中,EEENS为配电网的预期未供应电量,λb、λivt、λDCB分别为配电网支路b上线路、逆变器和直流断路器的故障率,τline、τivt、τDCB分别为配电网支路b上线路、逆变器和直流断路器的故障修复时间,Ω为配电网的支路集合,为负荷节点i所安装分布式发电装置的直流安装容量,Hb为配电网支路b故障而导致的下游负荷功率缺失量,ψ为配电网的节点i集合。
本申请还提供一种终端设备,包括处理器和存储器;
所述存储器,用于存储程序代码,并将所述程序代码传输给所述处理器;
所述处理器,用于根据所述程序代码中的指令执行上述所述的配电网预期未供应电量计算方法。
从以上技术方案可以看出,本申请实施例具有以下优点:该配电网预期未供应电量计算方法、***及设备,该方法包括获取含有分布式电源的配电网的网络拓扑模型,根据网络拓扑模型的拓扑连接关系构成节点关联矩阵;获取网络拓扑模型中各个母线节点的负荷需求和各个分布式发电装置的安装容量,根据所有的负荷需求和安装容量构建虚拟潮流平衡模型;根据故障条件和虚拟潮流平衡模型,构建配电网发生线路元件故障时各支路下游的实际缺供电负荷模型,基于实际缺供电负荷模型获得配电网发生线路元件故障时各支路下游的实际缺供电负荷;获取配电网中电气元件的故障参数,根据实际缺供电负荷和故障参数构建预期未供应电量计算模型;根据虚拟潮流平衡模型、实际缺供电负荷模型和预期未供应电量计算模型构建线性规划优化模型,对线性规划优化模型进行优化求解,得到配电网的预期未供应电量。该配电网预期未供应电量计算方法构建的线性规划优化模型考虑了分布式电源上游发生故障孤岛供电情况且该线性规划优化模型是由数学表达式组成的,便于配电网的预期未供应电量的计算,让该配电网预期未供应电量计算方法实现与数学规划模块相融合,使得计算配电网的预期未供应电量的速度快、计算准确,效率高,解决了现有对配电网可靠性指标求解方式难以与数学规划模型结合,存在效率不高的技术问题。
附图说明
为了更清楚地说明本申请实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本申请的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其它的附图。
图1为本申请实施例所述的配电网预期未供应电量计算方法的步骤流程图;
图2为本申请实施例所述的配电网预期未供应电量计算***的框架图。
具体实施方式
为使得本申请的发明目的、特征、优点能够更加的明显和易懂,下面将结合本申请实施例中的附图,对本申请实施例中的技术方案进行清楚、完整地描述,显然,下面所描述的实施例仅仅是本申请一部分实施例,而非全部的实施例。基于本申请中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其它实施例,都属于本申请保护的范围。
本申请提出一种配电网预期未供应电量计算方法、***及设备,用于解决了现有对配电网可靠性指标求解方式难以与数学规划模型结合,存在效率不高的技术问题。
实施例一:
图1为本申请实施例所述的配电网预期未供应电量计算方法的步骤流程图。
如图1所示,本申请提供一种配电网预期未供应电量计算方法,包括以下步骤:
S10.获取含有分布式电源的配电网的网络拓扑模型,根据网络拓扑模型的拓扑连接关系构成节点关联矩阵。
需要说明的是,在步骤S10中,主要是根据配电网的网络拓扑模型中拓扑连接关系构建节点关联矩阵A。在本实施例中,节点关联矩阵A中第i行第b列对应的元素记为Ai,b,Ai,b为节点i与支路b的连接关系,Ai,b的数值取值规则为:支路远离节点时取1,支路朝向节点时取-1,支路与节点无连接关系时取0。
S20.获取网络拓扑模型中各个母线节点的负荷需求和各个分布式发电装置的安装容量,根据所有的负荷需求和安装容量构建虚拟潮流平衡模型。
需要说明的是,在步骤S20中,先获取配电网的网络拓扑模型中各个母线节点的负荷需求以及安装在各个节点上的分布式发电装置的安装容量;之后对获取的各个节点的负荷需求和各个安装容量分别构建虚拟潮流平衡方程。
进一步地,虚拟潮流平衡模型包括负荷需求虚拟潮流平衡方程和容量虚拟潮流平衡方程,负荷需求虚拟潮流平衡方程为:
容量虚拟潮流平衡方程为:
式中,Ω为配电网的支路集合,Db为支路b下游的总最大负荷需求,为支路b在正方向上的虚拟负荷需求潮流,为支路b在负方向上的虚拟负荷需求潮流,MD为负荷需求虚拟潮流平衡方程取足够大的正整数,为第一辅助二值变量,Ai,b为节点关联矩阵中第i行第b列对应的元素数值,ψ为配电网的节点i集合,Di为负荷节点i的负荷需求;Ds为电源节点s的负荷需求;Gs为电源节点s下游的总负荷需求,CDGb为支路b下游的所有分布式发电装置计及换流效率的供电容量,为支路b在正方向上的虚拟负荷安装容量潮流,为支路b在负方向上的虚拟负荷安装容量潮流,为第二辅助二值变量,MDG为容量虚拟潮流平衡方程取足够大的正整数;为负荷节点i所安装分布式发电装置的交流安装容量,为负荷节点i所安装分布式发电装置的直流安装容量,为电源节点i所安装分布式发电装置的交流安装容量,为电源节点i所安装分布式发电装置的直流安装容量,ηivt为配电网中逆变器的换流效率,为电源节点s下游的计及换流效率分布式发电装置的安装总容量。
需要说明的是,在负荷需求虚拟潮流平衡方程是基于不考虑分布式发电装置DG供电情况下配电网的负荷需求构建的,在负荷需求虚拟潮流平衡方程中第一辅助二值变量是用于约束和只能有一个取非0值。容量虚拟潮流平衡方程是基于仅考虑分布式发电装置DG供电情况下配电网的安装容量构建的。在容量虚拟潮流平衡方程的第二辅助二值变量是用于约束和只能有一个取非0值。
S30.根据故障条件和虚拟潮流平衡模型,构建配电网发生线路元件故障时各支路下游的实际缺供电负荷模型,基于实际缺供电负荷模型获得配电网发生线路元件故障时各支路下游的实际缺供电负荷。
需要说明的是,在步骤S30中,主要是构建配电网发生线路元件故障时各支路下游的实际缺供电负荷模型并通过该实际缺供电负荷模型计算配电网发生线路元件故障时各支路下游的实际缺供电负荷。在本实施例中,故障条件是配电网的上游线路故障后,分布式发电装置DG以孤岛方式实现配电网供电。即是计算得到的实际缺供电负荷是配电网考虑了故障后分布式发电装置DG为断电区域供电情境下的实际缺电的量。
进一步地,实际缺供电负荷模型采用损失负荷量方程表示,损失负荷量方程为:
式中,Ω为配电网的支路集合,Db为支路b下游的总最大负荷需求,CDGb为支路b下游的所有分布式发电装置计及换流效率的供电容量,Hb为配电网支路b故障而导致的下游负荷功率缺失量,MCP为损失负荷量方程取的足够大的正整数,为第三辅助二值变量。
需要说明的是,在损失负荷量方程中第三辅助二值变量是用于保证Hb非负,当Db大于CDGb时,只有取1才能满足负荷需求虚拟潮流平衡方程、容量虚拟潮流平衡方程和损失负荷量方程中所有不等式约束,此时Hb就被约束为Db-CDGb,即线路b故障而导致的下游负荷功率缺失量等于下游最大负荷需求总量与分布式发电装置DG供电量的差值;而当Db小于CDGb时,只有取0才能满足负荷需求虚拟潮流平衡方程、容量虚拟潮流平衡方程个损失负荷量方程所有不等式约束,此时Hb就被约束为0,而不会出现Hb取负数这种不符合实际的情况。
S40.获取配电网中电气元件的故障参数,根据实际缺供电负荷和故障参数构建预期未供应电量计算模型。
需要说明的是,在步骤S40中,根据获取故障参数和步骤S3得到的实际缺供电负荷Hb结合状态枚举法构建预期未供应电量计算模型,为后续计算配电网的预期未供应电量提供数据。其中,配电网中电气元件包括配电网支路上的线路、逆变器和直流断路器,故障参数包括故障率和故障修改时间,则配电网中电气元件的故障参数指的是线路故障率、逆变器故障率、直流断路器故障率、线路故障修复时间、逆变器故障修复时间和直流断路器故障修复时间。
进一步地,预期未供应电量计算模型包括:用于计算配电网预期未供应电量的电量计算方程,电量计算方程为:
式中,EEENS为配电网的预期未供应电量,λb、λivt、λDCB分别为配电网支路b上线路、逆变器和直流断路器的故障率,τline、τivt、τDCB分别为配电网支路b上线路、逆变器和直流断路器的故障修复时间,Ω为配电网的支路集合,为负荷节点i所安装分布式发电装置的直流安装容量,Hb为配电网支路b故障而导致的下游负荷功率缺失量,ψ为配电网的节点i集合。
S50.根据虚拟潮流平衡模型、实际缺供电负荷模型和预期未供应电量计算模型构建线性规划优化模型,对线性规划优化模型进行优化求解,得到配电网的预期未供应电量。
需要说明的是,在步骤S50中,主要是根据步骤S20至步骤S40中得到的负荷需求虚拟潮流平衡方程、容量虚拟潮流平衡方程、损失负荷量方程和电量计算方程构成线性规划优化模型,之后对线性规划优化模型进行优化求解,计算得到配电网的预期未供应电量。
进一步地,对线性规划优化模型进行优化求解,得到配电网的预期未供应电量包括:采用Gurobi自带的分支定界法和目标条件对线性规划优化模型进行优化求解,得到配电网的预期未供应电量;其中,目标条件为选取实际缺供电负荷的最小值和Gurobi中松弛解与收敛解之间的距离为0.1%。
需要说明的是,对线性规划优化模型进行优化求解也可以采用其他与Gurobi自带的分支定界法相同功能的商业求解器在目标条件下进行优化求解。
本申请提供的配电网预期未供应电量计算方法,包括获取含有分布式电源的配电网的网络拓扑模型,根据网络拓扑模型的拓扑连接关系构成节点关联矩阵;获取网络拓扑模型中各个母线节点的负荷需求和各个分布式发电装置的安装容量,根据所有的负荷需求和安装容量构建虚拟潮流平衡模型;根据故障条件和虚拟潮流平衡模型,构建配电网发生线路元件故障时各支路下游的实际缺供电负荷模型,基于实际缺供电负荷模型获得配电网发生线路元件故障时各支路下游的实际缺供电负荷;获取配电网中电气元件的故障参数,根据实际缺供电负荷和故障参数构建预期未供应电量计算模型;根据虚拟潮流平衡模型、实际缺供电负荷模型和预期未供应电量计算模型构建线性规划优化模型,对线性规划优化模型进行优化求解,得到配电网的预期未供应电量。该配电网预期未供应电量计算方法构建的线性规划优化模型考虑了分布式电源上游发生故障孤岛供电情况且该线性规划优化模型是由数学表达式组成的,便于配电网的预期未供应电量的计算,让该配电网预期未供应电量计算方法实现与数学规划模块相融合,使得计算配电网的预期未供应电量的速度快、计算准确,效率高,解决了现有对配电网可靠性指标求解方式难以与数学规划模型结合,存在效率不高的技术问题。
在本申请的实施例中,该配电网预期未供应电量计算方法也可以计算配电网***的平均停电频率指数以及配电网***的平均停电时间指数,也可针对用户数量构建虚拟潮流平衡方程,最终利用状态枚举的思路求解得到。通过该配电网预期未供应电量计算方法的步骤就可以快速计算得到辐射状配电网故障预期未供应电量,从而评估配电网的可靠性程度。其中,该配电网预期未供应电量计算方法也可以很方便地随着构建模型的精细化程度不断提高进行进一步地优化,例如考虑不同负荷重要程度时只需要对每个负荷进行加权计算也可以。
实施例二:
图2为本申请实施例所述的配电网预期未供应电量计算***的框架图。
如图2所示,本申请还提供一种配电网预期未供应电量计算***,包括矩阵构建模块10、潮流模型构建模块20、供电负荷计算模块30、电量模型构建模块40和供应量计算模块50;
矩阵构建模块10,用于获取含有分布式电源的配电网的网络拓扑模型,根据网络拓扑模型的拓扑连接关系构成节点关联矩阵;
潮流模型构建模块20,用于获取网络拓扑模型中各个母线节点的负荷需求和各个分布式发电装置的安装容量,根据所有的负荷需求和安装容量构建虚拟潮流平衡模型;
供电负荷计算模块30,用于根据故障条件和虚拟潮流平衡模型,构建配电网发生线路元件故障时各支路下游的实际缺供电负荷模型,基于实际缺供电负荷模型获得配电网发生线路元件故障时各支路下游的实际缺供电负荷;
电量模型构建模块40,用于获取配电网中电气元件的故障参数,根据实际缺供电负荷和故障参数构建预期未供应电量计算模型;
供应量计算模块50,用于根据虚拟潮流平衡模型、实际缺供电负荷模型和预期未供应电量计算模型构建线性规划优化模型,对线性规划优化模型进行优化求解,得到配电网的预期未供应电量。
在本申请实施例中,供应量计算模块50还用于采用Gurobi自带的分支定界法和目标条件对线性规划优化模型进行优化求解,得到配电网的预期未供应电量;其中,目标条件为选取实际缺供电负荷的最小值和Gurobi中松弛解与收敛解之间的距离为0.1%。
在本申请实施例中,虚拟潮流平衡模型包括:负荷需求虚拟潮流平衡方程和容量虚拟潮流平衡方程,负荷需求虚拟潮流平衡方程为:
容量虚拟潮流平衡方程为:
式中,Ω为配电网的支路集合,Db为支路b下游的总最大负荷需求,为支路b在正方向上的虚拟负荷需求潮流,为支路b在负方向上的虚拟负荷需求潮流,MD为负荷需求虚拟潮流平衡方程取足够大的正整数,为第一辅助二值变量,Ai,b为节点关联矩阵中第i行第b列对应的元素数值,ψ为配电网的节点i集合,Di为负荷节点i的负荷需求;Ds为电源节点s的负荷需求;Gs为电源节点s下游的总负荷需求,CDGb为支路b下游的所有分布式发电装置计及换流效率的供电容量,为支路b在正方向上的虚拟负荷安装容量潮流,为支路b在负方向上的虚拟负荷安装容量潮流,为第二辅助二值变量,MDG为容量虚拟潮流平衡方程取足够大的正整数;为负荷节点i所安装分布式发电装置的交流安装容量,为负荷节点i所安装分布式发电装置的直流安装容量,为电源节点i所安装分布式发电装置的交流安装容量,为电源节点i所安装分布式发电装置的直流安装容量,ηivt为配电网中逆变器的换流效率,为电源节点s下游的计及换流效率分布式发电装置的安装总容量。
在本申请实施例中,实际缺供电负荷模型采用损失负荷量方程表示,损失负荷量方程为:
式中,Ω为配电网的支路集合,Db为支路b下游的总最大负荷需求,CDGb为支路b下游的所有分布式发电装置计及换流效率的供电容量,Hb为配电网支路b故障而导致的下游负荷功率缺失量,MCP为损失负荷量方程取的足够大的正整数,为第三辅助二值变量;
预期未供应电量计算模型包括:用于计算配电网预期未供应电量的电量计算方程,电量计算方程为:
式中,EEENS为配电网的预期未供应电量,λb、λivt、λDCB分别为配电网支路b上线路、逆变器和直流断路器的故障率,τline、τivt、τDCB分别为配电网支路b上线路、逆变器和直流断路器的故障修复时间,Ω为配电网的支路集合,为负荷节点i所安装分布式发电装置的直流安装容量,Hb为配电网支路b故障而导致的下游负荷功率缺失量,ψ为配电网的节点i集合。
需要说明的是,实施例二中模块的内容对应于实施例一方法中的步骤,实施例一方法步骤的内容已经在实施例一中详细阐述了,在实施例二不再对***中模块的内容重复阐述。
实施例三:
本申请还提供一种终端设备,包括处理器和存储器;
存储器,用于存储程序代码,并将程序代码传输给处理器;
处理器,用于根据程序代码中的指令执行上述的配电网预期未供应电量计算方法。
所属领域的技术人员可以清楚地了解到,为描述的方便和简洁,上述描述的***,装置和单元的具体工作过程,可以参考前述方法实施例中的对应过程,在此不再赘述。
在本申请所提供的几个实施例中,应该理解到,所揭露的***,装置和方法,可以通过其它的方式实现。例如,以上所描述的装置实施例仅仅是示意性的,例如,所述单元的划分,仅仅为一种逻辑功能划分,实际实现时可以有另外的划分方式,例如多个单元或组件可以结合或者可以集成到另一个***,或一些特征可以忽略,或不执行。另一点,所显示或讨论的相互之间的耦合或直接耦合或通信连接可以是通过一些接口,装置或单元的间接耦合或通信连接,可以是电性,机械或其它的形式。
所述作为分离部件说明的单元可以是或者也可以不是物理上分开的,作为单元显示的部件可以是或者也可以不是物理单元,即可以位于一个地方,或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部单元来实现本实施例方案的目的。
另外,在本申请各个实施例中的各功能单元可以集成在一个处理单元中,也可以是各个单元单独物理存在,也可以两个或两个以上单元集成在一个单元中。上述集成的单元既可以采用硬件的形式实现,也可以采用软件功能单元的形式实现。
所述集成的单元如果以软件功能单元的形式实现并作为独立的产品销售或使用时,可以存储在一个计算机可读取存储介质中。基于这样的理解,本申请的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的全部或部分可以以软件产品的形式体现出来,该计算机软件产品存储在一个存储介质中,包括若干指令用以使得一台计算机设备(可以是个人计算机,服务器,或者网络设备等)执行本申请各个实施例所述方法的全部或部分步骤。而前述的存储介质包括:U盘、移动硬盘、只读存储器(ROM,Read-OnlyMemory)、随机存取存储器(RAM,Random Access Memory)、磁碟或者光盘等各种可以存储程序代码的介质。
以上所述,以上实施例仅用以说明本申请的技术方案,而非对其限制;尽管参照前述实施例对本申请进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本申请各实施例技术方案的精神和范围。
Claims (10)
1.一种配电网预期未供应电量计算方法,其特征在于,包括以下步骤:
获取含有分布式电源的配电网的网络拓扑模型,根据所述网络拓扑模型的拓扑连接关系构成节点关联矩阵;
获取所述网络拓扑模型中各个母线节点的负荷需求和各个分布式发电装置的安装容量,根据所有的所述负荷需求和所述安装容量构建虚拟潮流平衡模型;
根据故障条件和所述虚拟潮流平衡模型,构建配电网发生线路元件故障时各支路下游的实际缺供电负荷模型,基于所述实际缺供电负荷模型获得配电网发生线路元件故障时各支路下游的实际缺供电负荷;
获取配电网中电气元件的故障参数,根据所述实际缺供电负荷和所述故障参数构建预期未供应电量计算模型;
根据所述虚拟潮流平衡模型、所述实际缺供电负荷模型和所述预期未供应电量计算模型构建线性规划优化模型,对所述线性规划优化模型进行优化求解,得到配电网的预期未供应电量。
2.根据权利要求1所述的配电网预期未供应电量计算方法,其特征在于,对所述线性规划优化模型进行优化求解,得到配电网的预期未供应电量包括:采用Gurobi自带的分支定界法和目标条件对所述线性规划优化模型进行优化求解,得到配电网的预期未供应电量;其中,所述目标条件为选取所述实际缺供电负荷的最小值和Gurobi中松弛解与收敛解之间的距离为0.1%。
3.根据权利要求1所述的配电网预期未供应电量计算方法,其特征在于,所述虚拟潮流平衡模型包括:负荷需求虚拟潮流平衡方程和容量虚拟潮流平衡方程,所述负荷需求虚拟潮流平衡方程为:
所述容量虚拟潮流平衡方程为:
式中,Ω为配电网的支路集合,Db为支路b下游的总最大负荷需求,为支路b在正方向上的虚拟负荷需求潮流,为支路b在负方向上的虚拟负荷需求潮流,MD为负荷需求虚拟潮流平衡方程取足够大的正整数,为第一辅助二值变量,Ai,b为节点关联矩阵中第i行第b列对应的元素数值,ψ为配电网的节点i集合,Di为负荷节点i的负荷需求;Ds为电源节点s的负荷需求;Gs为电源节点s下游的总负荷需求,CDGb为支路b下游的所有分布式发电装置计及换流效率的供电容量,为支路b在正方向上的虚拟负荷安装容量潮流,为支路b在负方向上的虚拟负荷安装容量潮流,为第二辅助二值变量,MDG为容量虚拟潮流平衡方程取足够大的正整数;为负荷节点i所安装分布式发电装置的交流安装容量,为负荷节点i所安装分布式发电装置的直流安装容量,为电源节点i所安装分布式发电装置的交流安装容量,为电源节点i所安装分布式发电装置的直流安装容量,ηivt为配电网中逆变器的换流效率,为电源节点s下游的计及换流效率分布式发电装置的安装总容量。
6.一种配电网预期未供应电量计算***,其特征在于,包括矩阵构建模块、潮流模型构建模块、供电负荷计算模块、电量模型构建模块和供应量计算模块;
所述矩阵构建模块,用于获取含有分布式电源的配电网的网络拓扑模型,根据所述网络拓扑模型的拓扑连接关系构成节点关联矩阵;
所述潮流模型构建模块,用于获取所述网络拓扑模型中各个母线节点的负荷需求和各个分布式发电装置的安装容量,根据所有的所述负荷需求和所述安装容量构建虚拟潮流平衡模型;
所述供电负荷计算模块,用于根据故障条件和所述虚拟潮流平衡模型,构建配电网发生线路元件故障时各支路下游的实际缺供电负荷模型,基于所述实际缺供电负荷模型获得配电网发生线路元件故障时各支路下游的实际缺供电负荷;
所述电量模型构建模块,用于获取配电网中电气元件的故障参数,根据所述实际缺供电负荷和所述故障参数构建预期未供应电量计算模型;
所述供应量计算模块,用于根据所述虚拟潮流平衡模型、所述实际缺供电负荷模型和所述预期未供应电量计算模型构建线性规划优化模型,对所述线性规划优化模型进行优化求解,得到配电网的预期未供应电量。
7.根据权利要求6所述的配电网预期未供应电量计算***,其特征在于,所述供应量计算模块还用于采用Gurobi自带的分支定界法和目标条件对所述线性规划优化模型进行优化求解,得到配电网的预期未供应电量;其中,所述目标条件为选取所述实际缺供电负荷的最小值和Gurobi中松弛解与收敛解之间的距离为0.1%。
8.根据权利要求6所述的配电网预期未供应电量计算***,其特征在于,所述虚拟潮流平衡模型包括:负荷需求虚拟潮流平衡方程和容量虚拟潮流平衡方程,所述负荷需求虚拟潮流平衡方程为:
所述容量虚拟潮流平衡方程为:
式中,Ω为配电网的支路集合,Db为支路b下游的总最大负荷需求,为支路b在正方向上的虚拟负荷需求潮流,为支路b在负方向上的虚拟负荷需求潮流,MD为负荷需求虚拟潮流平衡方程取足够大的正整数,为第一辅助二值变量,Ai,b为节点关联矩阵中第i行第b列对应的元素数值,ψ为配电网的节点i集合,Di为负荷节点i的负荷需求;Ds为电源节点s的负荷需求;Gs为电源节点s下游的总负荷需求,CDGb为支路b下游的所有分布式发电装置计及换流效率的供电容量,为支路b在正方向上的虚拟负荷安装容量潮流,为支路b在负方向上的虚拟负荷安装容量潮流,为第二辅助二值变量,MDG为容量虚拟潮流平衡方程取足够大的正整数;为负荷节点i所安装分布式发电装置的交流安装容量,为负荷节点i所安装分布式发电装置的直流安装容量,为电源节点i所安装分布式发电装置的交流安装容量,为电源节点i所安装分布式发电装置的直流安装容量,ηivt为配电网中逆变器的换流效率,为电源节点s下游的计及换流效率分布式发电装置的安装总容量。
9.根据权利要求6所述的配电网预期未供应电量计算***,其特征在于,所述实际缺供电负荷模型采用损失负荷量方程表示,所述损失负荷量方程为:
式中,Ω为配电网的支路集合,Db为支路b下游的总最大负荷需求,CDGb为支路b下游的所有分布式发电装置计及换流效率的供电容量,Hb为配电网支路b故障而导致的下游负荷功率缺失量,MCP为损失负荷量方程取的足够大的正整数,为第三辅助二值变量;
所述预期未供应电量计算模型包括:用于计算配电网预期未供应电量的电量计算方程,所述电量计算方程为:
10.一种终端设备,其特征在于,包括处理器和存储器;
所述存储器,用于存储程序代码,并将所述程序代码传输给所述处理器;
所述处理器,用于根据所述程序代码中的指令执行如权利要求1-6任意一项所述的配电网预期未供应电量计算方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211119436.3A CN115411735A (zh) | 2022-09-14 | 2022-09-14 | 一种配电网预期未供应电量计算方法、***及设备 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211119436.3A CN115411735A (zh) | 2022-09-14 | 2022-09-14 | 一种配电网预期未供应电量计算方法、***及设备 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN115411735A true CN115411735A (zh) | 2022-11-29 |
Family
ID=84166651
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202211119436.3A Pending CN115411735A (zh) | 2022-09-14 | 2022-09-14 | 一种配电网预期未供应电量计算方法、***及设备 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115411735A (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116526480A (zh) * | 2023-07-05 | 2023-08-01 | 山西中控绿源科技有限公司 | 基于智慧能源管理平台的分布式供电***及方法 |
CN118137510A (zh) * | 2024-05-06 | 2024-06-04 | 国网天津市电力公司城西供电分公司 | 面向多类型负荷重组的配电网供电能力挖潜方法与设备 |
CN118137510B (zh) * | 2024-05-06 | 2024-07-09 | 国网天津市电力公司城西供电分公司 | 面向多类型负荷重组的配电网供电能力挖潜方法与设备 |
-
2022
- 2022-09-14 CN CN202211119436.3A patent/CN115411735A/zh active Pending
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116526480A (zh) * | 2023-07-05 | 2023-08-01 | 山西中控绿源科技有限公司 | 基于智慧能源管理平台的分布式供电***及方法 |
CN116526480B (zh) * | 2023-07-05 | 2023-10-13 | 北京中芯标准科技有限公司 | 基于智慧能源管理平台的分布式供电***及方法 |
CN118137510A (zh) * | 2024-05-06 | 2024-06-04 | 国网天津市电力公司城西供电分公司 | 面向多类型负荷重组的配电网供电能力挖潜方法与设备 |
CN118137510B (zh) * | 2024-05-06 | 2024-07-09 | 国网天津市电力公司城西供电分公司 | 面向多类型负荷重组的配电网供电能力挖潜方法与设备 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102214922B (zh) | 一种电网规划方案的评估*** | |
Gu et al. | Reliability-based distribution network pricing | |
CN103514570A (zh) | 一种含分布式电源的配电网扩展规划综合优化方法 | |
CN109728579B (zh) | 一种配电网运行效率的评估方法、评估装置及评估设备 | |
CN106709651B (zh) | 一种基于风险理论的电力***安全性评估*** | |
CN112993974B (zh) | 电力现货日前市场出清电价计算方法、装置、终端及介质 | |
CN112183855A (zh) | 区域电力现货市场出清方法、装置和电子设备 | |
Nayak et al. | Optimal allocation of energy storage system and its benefit analysis for unbalanced distribution network with wind generation | |
CN110676839A (zh) | 一种交直流混合配电网可靠性评估方法及*** | |
CN115411735A (zh) | 一种配电网预期未供应电量计算方法、***及设备 | |
Yin et al. | Operational reliability assessment of distribution network with energy storage systems | |
CN110929220A (zh) | 一种配电网指标权重计算方法及装置 | |
Shayesteh et al. | System reduction techniques for storage allocation in large power systems | |
CN107179688B (zh) | 考虑蒙特卡罗状态抽样截尾的电力***可靠性分析方法 | |
Ebrahimi et al. | A conditional value at risk based stochastic allocation of SOP in distribution networks | |
Shaaban et al. | A multi-objective approach for optimal DG allocation | |
CN111311433A (zh) | 一种日前电力现货市场运行备用出清方法、装置和设备 | |
Pan et al. | A flexible black‐start network partitioning strategy considering subsystem recovery time balance | |
CN115907163A (zh) | 电力***运行状态预测方法、装置和计算机设备 | |
Mahdavi et al. | Evaluation of GA performance in TNEP considering voltage level, network losses and number of bundle lines | |
CN112016838B (zh) | 配电网能效指标体系的贡献率计算方法、***及终端设备 | |
Zhu | A new methodology of analytical formula deduction and sensitivity analysis of EENS in bulk power system reliability assessment | |
CN111401729A (zh) | 一种新增黑启动机组布点确定方法及*** | |
CN111404154A (zh) | 一种配电网供电能力优化方法、设备及存储介质 | |
Khodayari et al. | Customer reliability based energy & SR scheduling at hierarchical level II |
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 |