CN116415708B - 考虑置信水平的电网鲁棒规划方法 - Google Patents
考虑置信水平的电网鲁棒规划方法 Download PDFInfo
- Publication number
- CN116415708B CN116415708B CN202211732410.6A CN202211732410A CN116415708B CN 116415708 B CN116415708 B CN 116415708B CN 202211732410 A CN202211732410 A CN 202211732410A CN 116415708 B CN116415708 B CN 116415708B
- Authority
- CN
- China
- Prior art keywords
- planning
- power grid
- cost
- power
- node
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 55
- 238000011156 evaluation Methods 0.000 claims abstract description 26
- 230000005540 biological transmission Effects 0.000 claims description 53
- 239000011159 matrix material Substances 0.000 claims description 15
- 238000010248 power generation Methods 0.000 claims description 15
- 238000004364 calculation method Methods 0.000 claims description 13
- 238000010276 construction Methods 0.000 claims description 12
- 230000014509 gene expression Effects 0.000 claims description 10
- 238000012545 processing Methods 0.000 claims description 10
- 230000007774 longterm Effects 0.000 claims description 8
- 230000007613 environmental effect Effects 0.000 claims description 7
- 239000003344 environmental pollutant Substances 0.000 claims description 6
- 231100000719 pollutant Toxicity 0.000 claims description 6
- 238000012423 maintenance Methods 0.000 claims description 5
- 101150067055 minC gene Proteins 0.000 claims description 3
- 238000012163 sequencing technique Methods 0.000 claims description 3
- 238000004458 analytical method Methods 0.000 claims description 2
- 238000000354 decomposition reaction Methods 0.000 claims description 2
- 230000000694 effects Effects 0.000 claims description 2
- 238000012360 testing method Methods 0.000 claims description 2
- 238000012795 verification Methods 0.000 claims description 2
- 230000005611 electricity Effects 0.000 claims 1
- 238000013178 mathematical model Methods 0.000 claims 1
- 238000005457 optimization Methods 0.000 description 11
- 238000005516 engineering process Methods 0.000 description 5
- 238000011160 research Methods 0.000 description 4
- 238000007726 management method Methods 0.000 description 3
- 238000010606 normalization Methods 0.000 description 3
- 238000004088 simulation Methods 0.000 description 3
- 230000003044 adaptive effect Effects 0.000 description 2
- 230000007423 decrease Effects 0.000 description 2
- 238000011084 recovery Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 229910003460 diamond Inorganic materials 0.000 description 1
- 239000010432 diamond Substances 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
- 238000009472 formulation Methods 0.000 description 1
- 239000000446 fuel Substances 0.000 description 1
- 230000005283 ground state Effects 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 230000035515 penetration Effects 0.000 description 1
- 238000003786 synthesis reaction Methods 0.000 description 1
- 230000033772 system development Effects 0.000 description 1
- 238000010977 unit operation Methods 0.000 description 1
- 230000000007 visual effect Effects 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
-
- 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/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- 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/008—Circuit arrangements for ac mains or ac distribution networks involving trading of energy or energy transmission rights
-
- 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
- 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
-
- 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]
-
- 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
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Power Engineering (AREA)
- Economics (AREA)
- Human Resources & Organizations (AREA)
- Strategic Management (AREA)
- Health & Medical Sciences (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Marketing (AREA)
- General Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- Tourism & Hospitality (AREA)
- Game Theory and Decision Science (AREA)
- Quality & Reliability (AREA)
- Operations Research (AREA)
- Entrepreneurship & Innovation (AREA)
- Development Economics (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
本发明公开了考虑置信水平的电网鲁棒规划方法,包括:收集规划区电网参数、负荷分布情况、新能源出力的历史数据;构建不同置信水平的不确定集合,生成不同组极限场景;建立电网扩展规划模型;确定电网拓扑状态变量为第一阶段变量,将系统运行变量作为第二阶段变量,获得两阶段电网鲁棒扩展规划模型;对两阶段电网鲁棒规划模型进行求解,获得不同置信水平下的电网规划方案;建立电网扩展规划评价体系,对不同规划方案进行评估,确定电网最终规划方案。本发明考虑了不同置信水平的不确定集合对电网规划的影响,降低了规划方案的不确定性;本发明得到的规划方案综合效益更好,有效平衡了规划方案的鲁棒性与成本,可为电网规划人员提供有效参考。
Description
技术领域
本发明属于电力系统规划技术领域,具体涉及一种考虑置信水平的输电网鲁棒扩展规划方法。
背景技术
高比例可再生能源并网是未来电力系统发展的必然趋势,可再生能源接入带来的强不确定性使电力系统的基本形态和运行特性发生根本性改变。目前,国内外学者对随机变量的不确定性已有大量研究,其中鲁棒优化因具有求解效率高、无需随机变量的精确概率分布、求解结果能够适应最恶劣场景等优势,被广泛应用于电力系统调度运行。现有的鲁棒优化研究多简单考虑新能源出力的上下限,极少对新能源场出力的时空相关性进行研究,虽有少量鲁棒优化方法考虑了新能源场出力的时空相关性,但此类方法尚未应用于电力系统规划中;且现有研究一般采用不确定集合描述随机变量的波动范围,并未考虑不同置信水平下的不确定集合对调度运行或系统规划的影响。基于不确定集合决策得到电力系统的规划方案或调度计划,其决策结果的保守性与不确定集合的置信水平密切相关。随着置信水平的降低,决策结果鲁棒性降低、成本控制增强。
因此,研究将考虑随机变量时空相关性的鲁棒优化技术与电力系统规划技术有机结合,在此基础上改变不确定集合的置信水平,得到对应的不同组规划方案,通过评价机制决策出最终规划方案。此类新型规划技术方法更加契合高比例可再生能源并网带来的新特点,在保证系统安全可靠运行的同时,有效降低系统规划阶段的各项成本。
发明内容
本发明的目的是针对上述问题,提供一种考虑置信水平的输电网鲁棒扩展规划方法:通过缩放不确定集合的体积,改变考虑的可再生能源出力的历史场景数量,得到不同置信度的不确定集合,提取不同组极限场景;建立极限场景法的两阶段输电网规划模型,求解可得不同组电网规划方案;建立高比例新能源并网下电网扩展规划评价体系,对不同规划方案进行评估,根据综合评估分值,确定电网最终规划方案,有效平衡了规划方案的鲁棒性与成本,为规划人员提供参考。
本发明的技术方案是考虑置信水平的输电网鲁棒规划方法,包括以下步骤:
步骤1:收集规划区电网参数信息、负荷分布情况、可再生能源出力的历史数据;
步骤2:将收集的可再生能源出力的历史数据,以每一组出力的历史数据作为一个历史场景;获取历史数据的统计信息,考虑可再生能源场站出力的随机波动性和时空相关性,构建不确定集合;在不确定集合的基础上,改变不确定集合的置信水平,确定可再生能源出力的极限值,生成不同组极限场景;
步骤3:以规划水平年总成本最低为优化目标,综合考虑电网拓扑重构约束、投资预算约束、系统潮流约束、节点功率平衡约束、常规机组出力约束、弃风约束和切负荷约束,建立考虑N-1预想故障的输电网中长期扩展规划模型;
步骤4:对步骤3的模型进行线性处理,确定电网拓扑状态变量为第一阶段变量,将系统运行变量作为第二阶段变量,得到两阶段输电网中长期鲁棒扩展规划模型;
步骤5:对步骤4的两阶段输电网中长期鲁棒扩展规划模型进行求解,获得不同置信水平下的电网规划方案;
步骤6:建立高比例新能源并网下输电网扩展规划评价体系,对不同规划方案进行评估;
步骤7:根据综合评估分值,确定电网最终规划方案。
进一步地,所述电网参数信息包括电网各输电走廊已存线路数、各输电走廊可扩建的回数、各走廊输电线长度、各走廊输电线电抗、各走廊输电线传输容量、各条线路的编号以及首尾节点编号。
优选地,所述利用可再生能源出力的历史数据构建不确定集合,具体形式为广义凸包不确定集合。
广义凸包不确定集合从椭球不确定集合发展而来,不仅计及随机变量之间的相关性,还保留了盒式集合的线性表示方法;当有两个随机变量时,生成的广义凸包在二维空间呈现菱形形状,当有n个随机变量时,生成的广义凸包在n维空间呈现多面凸包体形状,n≥3。
n个可再生场站出力构建的广义凸包不确定集合的计算式如下:
式中Ucor表示广义凸包不确定性集合;χ为n个可再生能源出力的历史场景;为对应节点的可再生能源出力的第i个极限场景;Ne为历史场景数量。
进一步地,步骤2中,所述改变不确定集合的置信水平,置信水平为凸包内所考虑的场景的概率数,表达式为:
t=Pr{χ∈(ωP-1Ue′+c)}ω∈[0,kmax] (2)
式中:ω为置信度t对应的常量;Pr{·}为概率式;P为对椭球集合的对称正定系数矩阵Q进行正交分解的正交矩阵;Ue′为轴向初始极限场景集合;c为椭球集合的中心;kmax为放缩系数。
进一步地,步骤2,所述极限场景为任意一个可再生能源出力取极限值时的场景,在广义凸包中取凸包的顶点进行描述。
基于广义凸包不确定集合提取的极限场景集的计算式为:
式中,为极限场景集合。
步骤3中,考虑N-1预想故障的输电网扩展规划模型的目标函数为线路投资成本与常规机组运行成本、弃风弃负荷风险成本之和最小:
minCTotal=CInv+mT(CGen+CRis) (4)
式中,CTotal为规划总成本;CInv为线路投入成本;CGen为常规机组的发电成本;CRis为系统弃可再生能源成本和切负荷风险成本;m为系统运行成本与投资成本之间的等效因子;T为全年运行时长。
考虑N-1预想故障的输电网扩展规划模型的约束条件包括:
1)电网拓扑重构约束:
式中:为节点i、j之间第k回线路的0-1决策变量;/>为输电走廊ij待建线路集合;/>分别为输电走廊ij所允许建设线路条数的最小值和最大值;其中,式(5)为待建线路状态约束,当/>值为1时表示节点i、j之间第k回线路投建,若为0则表示不投建;式(6)为架线顺序约束,对于端点相同的输电线路,应满足建设序列回路约束。即在同一建线走廊下,若第k回线路投建,则第k+1回线路可能投建,若第k回线路不投建,则第k+1回及以上的所有线路均不投建;式(7)架线总数上下限约束,每条线路走廊通道的新建线路数应小于所允许的最大可建线路数。
2)投资预算约束:
式中,Ω为全部输电线路集;为走廊ij内第k条线路的投资费用;ψ为电网待建线路投资成本总预算;
3)系统潮流约束:
其中,nij为输电走廊ij已建线路数式(9)、(10)为已存线路潮流约束,式(11)、(12)为新建线路潮流约束;和/>分别为线路i、j间第k回线路的有功输送量与线路电纳;θi、θj分别为节点i、j的电压相角;/>为线路ij的最大传输容量;
4)节点功率平衡约束:
θmin≤θi≤θmax i∈N (14)
θref=0 (15)
其中,式(13)为节点功率平衡约束,式(14)为节点电压相角约束,式(15)为参考节点电压相角约束;式中,NG、NW分别为火电机组和风电机组接入的节点集合;Nstart、Nend分别为以节点i为始节点和末节点的节点集合;Pi G、Pi W和Pi D分别为节点i处的火电机组有功出力、风电机组有功出力和有功负荷;ΔPi W、ΔPi D分别为节点i处的弃风和切负荷量;θmax、θmin分别为节点i电压相角的最大、最小值;θref为系统中的平衡节点处电压相角值;
5)常规机组出力约束:
Pi G,min≤Pi G≤Pi G,max i∈NG (16)
式中,Pi G,max和Pi G,min分别为节点i处火力发电机允许的最大和最小有功出力。火力发电机组有功出力水平应在规定范围内。
6)弃风约束
0≤ΔPi W≤ri WPi W i∈NW (17)
式中,ri W为节点i所允许的最大弃风比例。
7)切负荷约束
0≤ΔPi D≤ri DPi D i∈ND (18)
式中,ri D为节点i所允许的最大切负荷比例。
进一步地,步骤4对考虑N-1预想故障的输电网扩展规划模型进行线性处理,经过逐一线性化处理后,原混合整数非线性规划模型转化为混合整数线性规划。
根据变量是否随可再生能源出力变化,确定电网拓扑状态变量为第一阶段变量,将系统运行变量作为第二阶段变量;电网规划第一阶段变量在运行模拟之前已确定,运行过程中不随可再生能源随机出力波动而变化,第二阶段变量在第一阶段决策变量确定后,与可再生能源出力随机场景相适应的可调节变量;采用极限场景约束形式描述可再生能源的随机出力,获得采用极限场景法的两阶段电网鲁棒规划模型:
式中,x为决策变量,即第一阶段变量;y为运行变量,即第二阶段变量;χ1,…,χV分别为随机变量的极限场景,V为极限场景个数;h(·)、q(·)分别代表等式、不等式约束条件;C(x,y)表示极限场景下电网规划模型的目标函数。
步骤6中,所述建立高比例新能源并网下电网扩展规划评价体系,考虑了高比例新能源并网对电网规划的影响,采用三层评估体系,其中目标层为“高比例新能源并网下电网规划评价指标体系”;准层次有“可靠性”、“成本控制”、“环保性”;指标层有“设备满载率”、“投资成本”、“运行成本”、“风险成本”、“清洁能源发电占比”、“污染物排放量”。根据指标层各项指标的数学表达式,计算不同组规划方案的数值。
进一步地,采用层次分析法对不同方案进行评估,具体步骤包括:根据专家经验将各层指标按重要程度排序,采用AHP九级标度法确定相邻指标的重要性标度值,构建相应的判断矩阵;设置准则层权重系数并进行一致性检验;设置综合权重系数并进行一致性校验;根据专家打分和综合权重值,得到最终规划方案。
相比现有技术,本发明的有益效果包括:
1)本发明采用广义凸包不确定集合处理可再生能源随机出力,相比传统的不确定集合,在考虑可再生能源出力具有随机波动性的基础上同时考虑了时空相关性,可有效减少不确定集合的空间体积。
2)本发明在广义凸包不确定集合的基础上提取极限场景,利用极限场景约束替代所有由随机变量确定的场景约束,相比传统鲁棒规划技术,较为直观易懂,无需复杂的模型转换处理,在实际工程中易于推广应用。
3)本发明所提考虑置信水平的输电网鲁棒扩展规划方法,通过改变考虑的历史场景的数量,获得不同置信水平下的广义凸包,其对应不同组电网规划方案,有效改善了电网拓扑结构,可适应规划期内电源、负荷的动态波动,保证电网安全可靠运行。
4)本发明所提考虑置信水平的输电网鲁棒扩展规划方法,通过建立高比例新能源并网下电网扩展规划评价体系,并采用层次分析法进行评估,综合考虑了电网的可靠性、成本和环保性,经评价后所得最终方案有效平衡了鲁棒性与成本,可为规划人员提供有效参考。
附图说明
下面结合附图和实施例对本发明作进一步说明。
图1是本发明实施例的考虑置信水平的输电网鲁棒扩展规划方法的流程示意图。
图2是本发明实施例的Garver-6节点系统拓扑结构。
图3是本发明实施例的广义凸包不确定集合构建过程的流程示意图。
图4是本发明实施例的广义凸包不确定集合鲁棒性调节示意图。
图5是本发明实施例的高比例风电并网下输电网规划评价指标体系。
图6是本发明实施例的五种规划方案的各项评估分值。
具体实施方式
如图1所示,考虑置信水平的输电网鲁棒规划方法,包括以下步骤:
步骤1:收集规划区域的电网各项信息、规划水平年的电源选址和负荷分布情况、可再生分布式电源出力的历史数据;
实施例的网架拓扑结构如图2所示,各项电网技术参数如表1所示。系统基准电压为220kV,基准功率为100MVA,选取节点2为参考节点。系统基本配置情况:节点1、2、6处配置火电机组,具体参数如表2所示;节点4、6处分别配置1台150MW的风电机组;系统总负荷为600MW,负荷分布情况如表3所示。线路投资费用为93168美元/km,线路维护管理费用为4658.4美元/km,资金折现率为10%,线路回收周期为10年。弃风惩罚成本系数为50美元/(MW·h),切负荷惩罚成本系数为500美元/(MW·h)。中国城市电网的可靠性为99.948%,考虑全年N-1故障为5h。
表1网架支路参数信息表
表2火电机组信息表
表3节点信息表
步骤2:将收集的新能源出力的历史数据,以每一组出力的历史数据作为一个历史场景,n个风电场的随机出力的历史场景集可表示为:
式中:Uh表示风电随机出力的历史场景集;χh,i为第i个历史场景;Nh为历史场景数量。
根据历史数据,获取新能源出力的预测值和预测误差,考虑可再生能源场站出力的时空相关性,参考公开号为CN107944638A的中国发明申请“一种基于时空相关性的新能源不确定集合建模方法”公开的时空相关性表示方法,构建广义凸包不确定集合,构建过程如图3所示;构建结果如图4所示。
改变凸包不确定集合的置信水平,确定新能源出力的极限值,生成不同组极限场景;取广义凸包不确定集合的顶点作为极限场景。
步骤3:以规划总成本最低为优化目标,综合考虑投资约束和系统安全运行约束,建立考虑N-1预想故障的输电网扩展规划模型,
步骤4:分析模型数学形式并进行线性处理;确定电网拓扑状态变量为第一阶段变量,将系统运行变量作为第二阶段变量,得到两阶段输电网中长期鲁棒扩展规划模型;
步骤4.1:分析N-1预想故障的输电网中长期扩展规划模型:
目标函数:以规划方案的总成本最低为优化目标,包括线路投资成本和系统在N-1预想故障下的运行成本。
其中,规划水平年线路投入成本包括新建线路的建线成本和规划期内所有线路的维护管理费用,具体表示为:
式中:Ω为全部输电线路集;分别为输电走廊ij待建线路和可存线路集合;分别为走廊ij内第k条线路的投资费用和维修管理费用;/>为节点i、j之间第k回线路的0-1决策变量;x为线路资金回收系数:
x的计算式如下
式中:r为每年的折扣率,y为输电网寿命年限。
机组消纳燃料的发电成本为:
CGen=ai(Pi G)2+biPi G+ci,i∈NG (23)
式中:NG为发电机节点集;ai、bi、ci分别为节点i处的发电机组的发电成本二次项系数、一次项系数和常数项系数;Pi G为节点i处发电机组的有功出力。
系统风险成本包括弃风风险成本和切负荷风险成本,计算式如下:
式中:分别为节点的弃风惩罚系数和弃负荷惩罚系数;NW表示风电机组接入的节点集合。
约束条件包括电网拓扑重构约束、投资预算约束、系统潮流约束、节点功率平衡约束、常规机组出力约束、弃风约束、切负荷约束。
综上,建立的考虑N-1预想故障的输电网中长期扩展规划模型的目标函数如下:
minCTotal=CInv+μT(CGen+CRis) (25)
步骤4.2:模型线性化;
步骤4.1的优化模型中,有二进制变量连续型变量:Pi G、ΔPi W、ΔPi D、/>θi等。式(23)即CGen的计算式中包含运行变量的二次项,且系统潮流约束中(θi-θj)与/>的乘积为双变量乘积。综上,此模型为一个复杂的混合整数非线性规划问题,求解效率不高,且难以获得最优解。因此,需对该模型进行线性化处理。
针对式(23)包含运行变量的二次项问题,采用分段线性化处理技术,参照2006年第21期《IEEE Transactions on Power Systems》刊登的Carrion M等的论文“AComputationally Efficient Mixed-Integer Linear Formulation for the ThermalUnit Commitment Problem”公开的分段线性化方法,该方法将曲线按因变量均分为S段,通过求出每段的函数值,进而求出每条线段的斜率li,s和截距δi,s。最终式(23)即CGen的计算式可用近似表示如下:
其中,各段的截距δi,s为新引入的连续变量,通过求解δi,s进而求解各火电机组出力Pi G。
针对系统潮流约束中双变量乘积问题,采用大M处理技术,参照2020年第35期《IEEE Transactions on Power Systems》刊登的Zhuo Z等的论文“IncorporatingMassive Scenarios in Transmission Expansion Planning With High RenewableEnergy Penetration”公开的大M法,处理过程如下所示:
式中,M为足够大正数,M>>0。
式(29)中当待选线路投建即时,M不起作用,式(29)等价转化为已建线路潮流方程式,即式(9),当待选线路不投建即/>时,由式(29)与式(12)共同约束可得/>即线路传输功率被限定为零。
综上,经过逐一线性化处理后,原混合整数非线性规划模型转化为如下所示的混合整数线性规划模型,
步骤4.3:采用极限场景法构建含N-1预想故障的两阶段输电网中长期鲁棒扩展规划模型;
N-1预想故障下输电网优化模型可描述为步骤4.2的混合整数线性规划模型,该优化问题可进一步抽象为如下优化模型:
式中:上标B和F分别表示系统基态与N-1预想故障态;x为决策变量即代表某条建线是否投建,当风电功率变化时,x在运行模拟之前已确定,运行过程中不随风电随机波动而变化,属于输电网规划第一阶段变量;y为运行变量,包括线路传输功率、节点电压相角、火电机组出力等,其值在第一阶段决策变量确定后,与风电出力随机场景χ相适应的可调节变量,属于输电网规划的第二阶段变量;χ为随机变量,即风电随机出力场景;h(·)、q(·)分别代表等式、不等式约束。从数学形式上看,该模型为随机混合整数线性规划模型。极限场景法能够有效处理该模型中的随机变量,该方法首先根据随机变量的取值空间确定极限场景,极限场景为随机变量都处于波动极限值的场景,一般取封闭凸包的顶点。然后用极限场景约束替代所有由随机变量确定的场景约束。可以证明,若线路投建状态变量x能适应极限场景χV,则对取值空间内的所有场景χi也具有适应性。
最终采用极限场景的两阶段输电网鲁棒规划模型可表示为:
式中,χ1,…,χV为随机变量取值空间的极限场景,V为极限场景个数。
步骤5:将步骤2中生成的不同组极限场景分别代入步骤4的规划模型中,对规划模型进行求解,获得不同置信水平下的电网规划方案;
采用MATLAB R2018b平台,利用YALMIP建模,调用CPLEX软件进行求解。本实施例中,选取中国西北某风场出力的实际数据,提取风电功率的均值、方差等统计信息,据此产生相关性系数为0.46的场景集,每组包含8755个场景,进行仿真计算。以此为输入数据,构建广义凸包不确定集合,在此基础上,通过调整放缩系数使集合处于不同置信水平,获取不同置信水平下不确定集合的风电极限场景,采用极限场景法得到不同组建线方案。鲁棒调节过程和运行结果如表4所示,其中NA表示弃风弃负荷约束条件越线,无可行解。
表4不同置信水平所得规划方案的各项成本
由表4可知,降低不确定集合的置信水平,所决策的规划方案成本控制增强,当置信水平进一步降低时(τ=74.09%),受规划模型中弃风和弃负荷约束限制,无可行解。因此,共有五种不同的线路规划方案,具体结果如表5所示。
表5五种规划方案
结合表4、表5可知,置信水平由100%将至98.72%时,此时不确定集合剔除了出现概率极低的极恶劣场景,线路投资成本大幅度降低;置信水平降低至96.78%时,为使总成本最低,在系统运行状态不变时,需增加新建线路以减少弃风和弃负荷量,故线路投资成本呈现小幅度回增;置信水平进一步降低时,为保证不出现严重的弃风和弃负荷现象,系统调整运行状态,火电机组出力增加,新建线路数减少;当置信水平降至74.09%时,采用此置信水平下的极限场景进行规划,系统出现较为严重的弃风或弃负荷现象,所得规划结果难以满足安全运行要求。
步骤6:建立高比例新能源并网下电网扩展规划评价体系,对不同规划方案进行评估;
广义凸包不确定集合涵盖了所有的可能场景,去除少量出现概率低的场景可缩小不确定集合,但鲁棒性会降低。为平衡规划方案的鲁棒性与成本,本发明提出考虑置信水平的输电网鲁棒规划方法,该方法经过步骤1至步骤2后,产生多个不同置信水平的广义凸包不确定集合,然后根据步骤3至步骤5,得到多个规划方案。为全面、直观地说明步骤5产生的五组不同规划方案的优劣,本发明充分考虑高比例可再生能源并网的特点,建立三层评估指标体系,采用层次分析法进行评估,具体过程包括:
1)建立高比例风电并网下输电网规划评价指标体系,体系结构如图5所示。
2)根据各项指标计算式,计算得各方案对应的指标值。本评价体系的各项指标表达式为:
A.可靠性指标:
设备满载率A11,设备若长期处于满载,可能出现过载或者故障情况,可说明电网结构设置的不合理性。
设备满载率A11的具体计算式为:
式中:N1为处于满载状态的设备个数,Nall为系统总设备个数。
B.成本控制指标:
设备投资成本B11,设备投资成本包含线路的新建成本和维护成本,与CInv的计算式一致,
B11=CInv (34)
系统运行成本B12的表达式为:
B12=CGen (35)
系统风险成本B13的表达式为:
B13=CRis (36)
C.环保性指标:
可再生能源发电占比C11,新能源发电可替代部分传统的火力发电,因此高比例风电的电网较传统电网环保效应上有显著提升。
可再生能源发电占比C11的具体计算公式为:
式中:为此刻系统风电出力总量;/>为对应时刻发电机出力总量。
火力发电的污染物排放量如表6所示,其中CO2含量最高,应给予较高重视。
CO2排放量C12的计算式如下:
式中:K(CO2)为CO2排放系数。
其他污染物排放量C13的计算式如下:
式中:K(X)分别为NOX、SO2、CO排放系数。
表6火力发电污染物排放系数
3)根据专家经验将各层指标按重要程度排序,采用AHP九级标度法确定相邻指标的重要性标度值,构建相应的判断矩阵,并归一化处理。
A.准则层判断矩阵以及归一化:
准则层判断矩阵为对应的最大特征值为λM,max=3.0291;则一致性指标为CI=0.01455;随机一致性比例为CR=0.0280<0.1,通过一致性校验。/>
B.指标层判断矩阵以及归一化:
指标层判断矩阵对应的最大特征值为λB,max=3.0858;一致性指标为CI=0.0429;随机一致性比例为CR=0.0825<0.1,通过一致性校验。
指标层判断矩阵对应的最大特征值为λC,max=3.0940;一致性指标为CI=0.0470;随机一致性比例为CR=0.0904<0.1,通过一致性校验。
4)计算各指标综合权重。
A.准则层各指标权重:
判断矩阵M的特征向量为ωM=[0.9640 0.2286 0.1355]T;则各指标的权重依次为:0.7259、0.1722、0.1019;
判断矩阵B的特征向量为ωB=[0.9048 0.1352 0.4038]T;则各项指标的权重依次为:0.6267、0.0937、0.2796;
判断矩阵C的特征向量为ωC=[0.9584 0.2601 0.1177]T;则各项指标的权重依次为:0.7173、0.1947、0.0880;
B.将权重折算到指标层:
WB13=0.0482;WC11=0.0731、WC12=0.0198、WC13=0.0090。
5)计算四组不同方案下各指标的数值,进行专家打分。根据各指标计算式,由规划和运维人员组成专家组评估得出各指标分值。专家打分和综合权重值,如表7所示。
表7五种规划方案各项指标的分值
步骤7:根据综合评估分值,确定电网最终规划方案。由步骤6的专家打分和综合权重值,计算方案的综合分值,结果如图6所示。
由图6可知,经上述体系评价,方案四评估分值最高,综合考虑了电网的可靠性、成本和环保性,可为规划人员提供有效参考。
Claims (4)
1.考虑置信水平的电网鲁棒规划方法,其特征在于,包括以下步骤:
步骤1:收集规划区电网参数信息、负荷分布情况、可再生能源出力的历史数据;
步骤2:根据可再生能源出力的历史数据获取统计信息,构建不同置信水平下的不确定集合,确定可再生能源出力的极限值,生成不同组极限场景;
步骤3:建立考虑N-1预想故障的输电网中长期扩展规划模型;
步骤4:对步骤3的模型进行线性处理,确定电网拓扑状态变量为第一阶段变量,将系统运行变量作为第二阶段变量,得到两阶段输电网中长期鲁棒扩展规划模型;
步骤5:对步骤4的两阶段输电网中长期鲁棒扩展规划模型进行求解,获得不同置信水平下的电网规划方案;
步骤6:建立高比例新能源并网下电网扩展规划评价体系,对不同规划方案进行评估;
步骤7:根据综合评估分值,确定电网最终规划方案;
步骤2中,所述置信水平为构建的不确定集合的置信水平,不确定集合的置信水平τ的表达式为:
τ=Pr{χ∈Rn×1|U(χ)≤ω} (1)
式中:ω为置信度τ对应的常量;Pr{·}表示概率;χ为n个新能源出力的历史场景;U为不确定集合;
步骤2中,n个新能源场站出力构建的不确定集合的计算式如下:
式中,Ucor表示考虑相关性的不确定性集合;为对应节点的新能源出力的第i个极限场景;Ne为历史场景数量;ui为极限场景/>的系数;
步骤2中,所述新能源出力的极限场景,计算式为:
式中,为极限场景集合;c为椭球不确定集合的中心;kmax为放缩系数;P为对椭球集合的对称正定系数矩阵进行正交分解的正交矩阵;U′e为轴向初始极限场景集合;
步骤3中,所述考虑N-1预想故障的输电网中长期扩展规划模型,其优化目标为投资成本与运行成本之和最小:
minCTotal=CInv+μT(CGen+CRis) (4)
式中CTotal为规划总成本;CInv为线路投入成本;CGen为常规机组的发电费用;CRis为系统弃可再生能源成本和切负荷风险成本;μ为系统运行成本与投资成本之间的等效因子,T为全年运行时长;
所述考虑N-1预想故障的输电网扩展规划模型的约束条件包括:
1)电网拓扑重构约束:
式中:为节点i、j之间第k回线路的0-1决策变量;/>为输电走廊ij待建线路集合;分别为输电走廊ij所允许建设线路条数的最小值和最大值;其中,式(5)为待建线路状态约束,当/>值为1时表示节点i、j之间第k回线路投建,若为0则表示不投建;式(6)为架线顺序约束,对于端点相同的输电线路,应满足建设序列回路约束,即在同一建线走廊下,若第k回线路投建,则第k+1回线路可能投建,若第k回线路不投建,则第k+1回及以上的所有线路均不投建;式(7)架线总数上下限约束,每条线路走廊通道的新建线路数应小于所允许的最大可建线路数;
2)投资预算约束:
式中,Ω为全部输电线路集;为走廊ij内第k条线路的投资费用;ψ为电网待建线路投资成本总预算;
3)系统潮流约束:
其中,nij为输电走廊ij已建线路数式(9)、(10)为已存线路潮流约束,式(11)、(12)为新建线路潮流约束;和/>分别为线路i、j间第k回线路的有功输送量与线路电纳;θi、θj分别为节点i、j的电压相角;/>为线路ij的最大传输容量;
4)节点功率平衡约束:
θmin≤θi≤θmax i∈N (14)
θref=0 (15)
其中,式(13)为节点功率平衡约束,式(14)为节点电压相角约束,式(15)为参考节点电压相角约束;式中,NG、NW分别为火电机组和风电机组接入的节点集合;Nstart、Nend分别为以节点i为始节点和末节点的节点集合;Pi G、Pi W和Pi D分别为节点i处的火电机组有功出力、风电机组有功出力和有功负荷;ΔPi W、ΔPi D分别为节点i处的弃风和切负荷量;θmax、θmin分别为节点i电压相角的最大、最小值;θref为系统中的平衡节点处电压相角值;
5)常规机组出力约束:
Pi G,min≤Pi G≤Pi G,max i∈NG (16)
式中,Pi G,max和Pi G,min分别为节点i处火力发电机允许的最大和最小有功出力;
6)弃风约束
0≤ΔPi W≤ri WPi W i∈NW (17)
式中,ri W为节点i所允许的最大弃风比例;
7)切负荷约束
0≤ΔPi D≤ri DPi D i∈ND (18)
式中,ri D为节点i所允许的最大切负荷比例。
2.根据权利要求1所述的考虑置信水平的电网鲁棒规划方法,其特征在于,步骤4中,所述两阶段含N-1预想故障的输电网鲁棒扩展规划模型的数学模型如下:
式中,x为决策变量,即第一阶段变量;y为运行变量,即第二阶段变量;χ1,…,χV为随机变量的极限场景,V为极限场景个数;h(·)、q(·)分别代表等式、不等式约束条件;C(x,y)表示极限场景下电网规划模型的目标函数。
3.根据权利要求2所述的考虑置信水平的电网鲁棒规划方法,其特征在于,所述评价体系为三层评估体系,具体包括:
1)可靠性指标:
设备满载率A11,设备若长期处于满载,可能出现过载或者故障情况,可说明电网结构设置的不合理性;
设备满载率A11的具体计算式为:
式中:N1为处于满载状态的设备个数,Nall为系统总设备个数;
2)成本控制指标:
设备投资成本B11,设备投资成本包含线路的新建成本和维护成本,与线路投入成本CInv表达式一致,即:
B11=CInv (21)
其中CInv表示线路投入成本;
系统运行成本B12的表达式为:
B12=CGen (22)
其中CGen表示常规机组的发电成本;
系统风险成本B13的表达式为:
B13=CRis (23)
其中CRis表示系统弃可再生能源成本和切负荷风险成本;
3)环保性指标:
可再生能源发电占比C11,新能源发电可替代部分传统的火力发电,因此高比例风电的电网较传统电网环保效应上有显著提升;
可再生能源发电占比C11的具体计算公式为:
式中:为此刻系统风电出力总量;/>为对应时刻发电机出力总量;
CO2排放量C12,计算式如下:
式中:K(CO2)为CO2排放系数;
其他污染物排放量C13,计算式如下:
式中:K(X)分别为NOX、SO2、CO污染物的排放系数。
4.根据权利要求3所述的考虑置信水平的电网鲁棒规划方法,其特征在于,步骤6中采用层次分析法对规划方案进行评估,具体过程包括:
根据专家经验将各层指标按重要程度排序,采用AHP九级标度法确定相邻指标的重要性标度值,构建相应的判断矩阵;
设置准则层权重系数并进行一致性检验;设置综合权重系数并进行一致性校验;
根据专家打分和综合权重值,得到最终方案。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211732410.6A CN116415708B (zh) | 2022-12-30 | 2022-12-30 | 考虑置信水平的电网鲁棒规划方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211732410.6A CN116415708B (zh) | 2022-12-30 | 2022-12-30 | 考虑置信水平的电网鲁棒规划方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN116415708A CN116415708A (zh) | 2023-07-11 |
CN116415708B true CN116415708B (zh) | 2024-03-29 |
Family
ID=87057134
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202211732410.6A Active CN116415708B (zh) | 2022-12-30 | 2022-12-30 | 考虑置信水平的电网鲁棒规划方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116415708B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117595267A (zh) * | 2023-09-20 | 2024-02-23 | 南京工业大学 | 一种交流潮流配电网综合能源分布式合作博弈方法及系统 |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107528314A (zh) * | 2017-08-15 | 2017-12-29 | 东南大学 | 一种不确定环境下考虑安全校验的电网鲁棒规划方法 |
CN108493998A (zh) * | 2018-04-14 | 2018-09-04 | 华中科技大学 | 考虑需求响应与n-1预想故障的鲁棒输电网规划方法 |
CN109242177A (zh) * | 2018-08-30 | 2019-01-18 | 国网江西省电力有限公司经济技术研究院 | 有源配电网规划方法 |
CN110472830A (zh) * | 2019-07-12 | 2019-11-19 | 中国电力科学研究院有限公司 | 一种考虑新能源接入的输电网规划方法和系统 |
KR20200057942A (ko) * | 2018-11-19 | 2020-05-27 | 한국전기연구원 | 신재생에너지 공급확대를 반영한 계통 순부하의 변동성 산정 방법 및 시스템 |
CN112671047A (zh) * | 2021-01-09 | 2021-04-16 | 三峡大学 | 考虑极限场景的主动配电网重构与无功联合鲁棒优化方法 |
CN113421004A (zh) * | 2021-06-30 | 2021-09-21 | 国网山东省电力公司潍坊供电公司 | 输配协同的主动配电网分布式鲁棒扩展规划系统及方法 |
CN114336663A (zh) * | 2022-01-07 | 2022-04-12 | 华北电力科学研究院有限责任公司 | 一种新型电力系统源网协同规划方法及装置 |
CN115471006A (zh) * | 2022-09-30 | 2022-12-13 | 国网能源研究院有限公司 | 一种考虑风电出力不确定性的电源规划方法及系统 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US11966989B2 (en) * | 2020-03-23 | 2024-04-23 | Xendee Corporation | Robust and fast design of microgrids, der systems, and other energy systems using staged hybrid investment planning |
-
2022
- 2022-12-30 CN CN202211732410.6A patent/CN116415708B/zh active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107528314A (zh) * | 2017-08-15 | 2017-12-29 | 东南大学 | 一种不确定环境下考虑安全校验的电网鲁棒规划方法 |
CN108493998A (zh) * | 2018-04-14 | 2018-09-04 | 华中科技大学 | 考虑需求响应与n-1预想故障的鲁棒输电网规划方法 |
CN109242177A (zh) * | 2018-08-30 | 2019-01-18 | 国网江西省电力有限公司经济技术研究院 | 有源配电网规划方法 |
KR20200057942A (ko) * | 2018-11-19 | 2020-05-27 | 한국전기연구원 | 신재생에너지 공급확대를 반영한 계통 순부하의 변동성 산정 방법 및 시스템 |
CN110472830A (zh) * | 2019-07-12 | 2019-11-19 | 中国电力科学研究院有限公司 | 一种考虑新能源接入的输电网规划方法和系统 |
CN112671047A (zh) * | 2021-01-09 | 2021-04-16 | 三峡大学 | 考虑极限场景的主动配电网重构与无功联合鲁棒优化方法 |
CN113421004A (zh) * | 2021-06-30 | 2021-09-21 | 国网山东省电力公司潍坊供电公司 | 输配协同的主动配电网分布式鲁棒扩展规划系统及方法 |
CN114336663A (zh) * | 2022-01-07 | 2022-04-12 | 华北电力科学研究院有限责任公司 | 一种新型电力系统源网协同规划方法及装置 |
CN115471006A (zh) * | 2022-09-30 | 2022-12-13 | 国网能源研究院有限公司 | 一种考虑风电出力不确定性的电源规划方法及系统 |
Non-Patent Citations (2)
Title |
---|
《一种考虑极限场景的配电网鲁棒扩展规划方法》;黄家祺;《电力建设》;第41卷(第7期);1-8 * |
《考虑直接负荷控制与风电不确定性的输电网扩展规划》;蒋霖;《电力系统保护与控制》;第48卷(第3期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN116415708A (zh) | 2023-07-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107506854B (zh) | 一种考虑差异化场景的220kV电网网架结构规划方法 | |
Zhao et al. | Distributed risk-limiting load restoration for wind power penetrated bulk system | |
CN107591841B (zh) | 适用于新能源大规模接入下的电网演化模拟方法 | |
CN115640963A (zh) | 一种考虑投资运营模式的海上风电接入系统鲁棒规划方法 | |
CN110543720A (zh) | 基于sdae-elm伪量测模型的状态估计方法 | |
CN103326348A (zh) | 一种提高地区电网供电能力分析及全过程在线监测的系统 | |
CN116415708B (zh) | 考虑置信水平的电网鲁棒规划方法 | |
CN115310378A (zh) | 一种极端台风灾害下电网韧性评估及差异化规划方法 | |
CN109672219A (zh) | 一种求解含风电场无功优化模型的方法及系统 | |
CN115017854A (zh) | 一种基于多维评估指标体系的配电网dg最大准入容量计算方法 | |
CN115000968A (zh) | 基于博弈和风电出力相关性的海上风电接入系统规划方法 | |
CN116187165A (zh) | 基于改进粒子群算法的电网弹性提升方法 | |
CN115051388A (zh) | 一种基于分布鲁棒的“源-网-荷-储”两阶段调度优化方法 | |
CN113344283A (zh) | 基于边缘智能的能源互联网新能源消纳能力评估方法 | |
CN116316611B (zh) | 基于低压台区下的供电方法及系统 | |
CN116722536A (zh) | 基于数据驱动的柔性配电网分布鲁棒扩展规划方法及系统 | |
CN115313422B (zh) | 主动支撑主网频率的电-气综合能源系统运行控制方法 | |
CN116231749A (zh) | 基于数字孪生的新能源电力系统调度方法 | |
CN114006410B (zh) | 一种基于机会约束规划的大规模海上风电接入点优化方法 | |
CN115954868A (zh) | 一种电网安全分析方法、装置及介质 | |
CN114925916A (zh) | 一种基于数字孪生的重要用户保电应急辅助决策方法 | |
CN104242354A (zh) | 计及风光出力相关性的新能源集中外送运行特性评估方法 | |
CN112651540A (zh) | 一种配电网规划项目投资优化方法 | |
Wang et al. | A bi-level programming model for distribution network expansion planning with distributed generations and energy storage systems | |
CN113505468B (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 |