CN109002938B - 考虑n-1安全准则的交直流混合配电网双层规划方法 - Google Patents
考虑n-1安全准则的交直流混合配电网双层规划方法 Download PDFInfo
- Publication number
- CN109002938B CN109002938B CN201811085672.1A CN201811085672A CN109002938B CN 109002938 B CN109002938 B CN 109002938B CN 201811085672 A CN201811085672 A CN 201811085672A CN 109002938 B CN109002938 B CN 109002938B
- Authority
- CN
- China
- Prior art keywords
- direct current
- representing
- current
- node
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 39
- 238000005457 optimization Methods 0.000 claims abstract description 26
- 230000002068 genetic effect Effects 0.000 claims abstract description 14
- 238000010276 construction Methods 0.000 claims description 47
- 210000000349 chromosome Anatomy 0.000 claims description 12
- 230000008901 benefit Effects 0.000 claims description 8
- 230000005540 biological transmission Effects 0.000 claims description 6
- 230000005611 electricity Effects 0.000 claims description 6
- 238000011084 recovery Methods 0.000 claims description 6
- 230000001172 regenerating effect Effects 0.000 claims description 5
- 230000008878 coupling Effects 0.000 claims description 4
- 238000010168 coupling process Methods 0.000 claims description 4
- 238000005859 coupling reaction Methods 0.000 claims description 4
- 238000006243 chemical reaction Methods 0.000 claims description 3
- 239000013256 coordination polymer Substances 0.000 claims description 3
- 238000013461 design Methods 0.000 claims description 3
- 239000004576 sand Substances 0.000 claims description 3
- 238000004364 calculation method Methods 0.000 description 7
- 238000004891 communication Methods 0.000 description 4
- 230000006872 improvement Effects 0.000 description 2
- 230000035699 permeability Effects 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 230000015556 catabolic process Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 230000006855 networking Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
Images
Classifications
-
- 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
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/004—Artificial life, i.e. computing arrangements simulating life
- G06N3/006—Artificial life, i.e. computing arrangements simulating life based on simulated virtual individual or collective life forms, e.g. social simulations or particle swarm optimisation [PSO]
-
- 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/0631—Resource planning, allocation, distributing or scheduling for enterprises or organisations
- G06Q10/06315—Needs-based resource requirements planning or 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
Landscapes
- Business, Economics & Management (AREA)
- Engineering & Computer Science (AREA)
- Human Resources & Organizations (AREA)
- Economics (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Strategic Management (AREA)
- General Physics & Mathematics (AREA)
- Marketing (AREA)
- Tourism & Hospitality (AREA)
- General Business, Economics & Management (AREA)
- Entrepreneurship & Innovation (AREA)
- Health & Medical Sciences (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- General Health & Medical Sciences (AREA)
- Game Theory and Decision Science (AREA)
- Development Economics (AREA)
- Educational Administration (AREA)
- Artificial Intelligence (AREA)
- Biomedical Technology (AREA)
- Biophysics (AREA)
- Computational Linguistics (AREA)
- Data Mining & Analysis (AREA)
- Evolutionary Computation (AREA)
- Life Sciences & Earth Sciences (AREA)
- Molecular Biology (AREA)
- Computing Systems (AREA)
- General Engineering & Computer Science (AREA)
- Mathematical Physics (AREA)
- Software Systems (AREA)
- Public Health (AREA)
- Water Supply & Treatment (AREA)
- Primary Health Care (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明公开了一种考虑N‑1安全准则的交直流混合配电网双层规划方法,括以下步骤:1)初始化遗传算法,随机构造直流网架结构;2)构造针对直流配电网并考虑N‑1安全准则的下层优化模型;3)求解下层模型,检验是否存在弃风弃光;4)基于合格的直流网架结构,构造与直流网架结构相配套的上层交流系统规划模型;5)求解上层模型,判断是否满足遗传算法迭代上限。本发明在交直流混合配电网的规划模型基础上考虑了直流配电网N‑1安全准则,在保证交直流配电网规划经济性的同时,提高了其运行可靠性和灵活性,为出台交直流混合配电网的规划方案提供重要的指导和帮助。
Description
技术领域
本发明配电网规划领域,具体来说,涉及一种考虑N-1安全准则的交直流混合配电网双层规划方法。
背景技术
由于可再生能源大量且集中分布于偏远地区或沿海地区,离负荷中心较远,并且随着可再生能源渗透率的不断提高,需要建立全新的配电网规划方法。考虑可再生能源出力的波动性和间歇性,传统的交流配电网往往难以避免弃风弃光现象的发生。目前,交直流混合配电网得益于其更大的传输容量、更小的线路损耗以及更高的可再生能源渗透率,已成为解决可再生能源接入电力系统的有效方法。此外,N-1安全准则是电网安全可靠运行的重要指标之一。根据北美电力可靠性委员会(NERC)的标准,当任何一条线路由于故障被切除后,系统应当依然能够保持稳定运行。考虑到可再生能源集中分布的偏远地区往往地形复杂、气候恶劣,较易出现安全性问题。一旦出现N-1故障,往往导致大量弃风弃光现象的发生,甚至导致电网的局部崩溃,造成严重后果。因此在进行配电网规划的时候,往往需要考虑N-1安全准则。
事实上,N-1安全准则的交直流混合配电网的规划尚未得到充分研究。一方面,交直流规划的研究目前仅限于大电网传输线规划或楼宇电网规划,交直流混合配电网的规划模型仍处于起步阶段,特别是交直流电网的潮流分布及计算尚未有统一的方法。有学者在微网层面提出了交直流规划的模型,但是所提模型并没有考虑网架的规划方案,无法解决实际问题。另一方面,虽然N-1安全准则在机组组合优化或大电网规划方面有较为详细的研究,但在交直流混合规划领域尚未得到考虑。并且,对于考虑了N-1安全准则之后模型复杂度大大提高,以及双层规划问题难以直接求解的问题,国内外尚未有较为理想的求解方法。
发明内容
技术问题:本发明提供一种在最小化交直流混合配电网的投资和运行费用的同时,保证直流侧网架在发生N-1安全故障的条件下,不出现弃风弃光现象的考虑N-1安全准则的交直流混合配电网双层规划方法。
技术方案:本发明的考虑N-1安全准则的交直流混合配电网双层规划方法,包括以下步骤:
1)初始化遗传算法参数,基于0-1编码方式随机构造直流网架结构,网架数量等于种群数,令当前网架编号k=1;
2)基于鲁棒规划,针对第k个表征直流网架结构的染色体个体,构造针对直流配电网并考虑N-1安全准则的下层优化模型;
3)求解所述步骤2)构造的下层优化模型,检验直流网架结构在最恶劣N-1故障下是否存在弃风弃光,若存在弃风弃光,表明当前直流网架结构无法通过检验,则舍弃当前直流网架结构,基于0-1编码方式随机重新生成当前直流网架结构,并返回步骤2);如果不存在弃风弃光,则判断k是否小于种群数,如果是,则令k=k+1,返回步骤2),否则,进入步骤4)
4)针对当前种群中的每一个直流网架结构,构造与该直流网架结构相配套的上层交流系统规划模型;
5)求解所述步骤4)构造的上层交流系统规划模型,将得到的优化结果作为每一个直流网架结构的适应度值,判断当前迭代次数是否达到遗传算法迭代次数上限,若未达到迭代上限,则对直流网架结构染色体进行选择、交叉、变异操作后,令k=1并返回步骤2),否则,将当前种群中的直流网架结构所对应的适应度最小的交直流网架结构作为结果输出。
进一步的,本发明方法中,所述步骤1)中,设置遗传算法参数为:种群数量为25个,迭代次数为100次,交叉概率为0.8,变异概率为0.6。根据规划实际需求确定交直流配电网系统拓扑结构,构造交流线路、直流线路、变电站、VSC、风电和光伏规划候选集。0-1编码方式为0代表不建设,1代表建设,构建表征直流网架结构的染色体种群。
进一步的,本发明方法中,所述步骤2)中构造的下层优化模型包括:
3.1目标函数
其中和分别对应第i个节点在第v个负荷等级和第t个阶段下的弃风和弃光功率;CCW和CCP分别对应弃风和弃光的惩罚系数;WD为风电场候选节点集合,PV为光伏电站候选节点集合;bn,t为0-1变量,表示第n条线路是否发生故障;ΩCTG为N-1故障点集合。
3.2约束条件
其中,变量上方的横线“—”是直流系统变量的标识符,变量上方的尖号“^”是通过遗传算法给定的直流拓扑变量的标识符;下标中v代表第v个负荷等级,t代表第t个阶段;表示第n条线路上流过的功率,表示第i个节点流入VSC的功率,和表示第i个节点风电和光伏的功率,ηVSC为VSC的能量转换效率系数,表示与第i个节点相连的线路的集合,表示直流系统中VSC节点的集合;和分别表示第n条线路的首末端节点电压的平方,Rn表示第n条线路的电阻,表示第n条线路上流过的电流的平方,M为一个非常大的正整数,表示第n条线路的建设情况,为直流候选线路的集合;和SrWD分别对应于风电的实际功率和额定容量,表示第i个节点风电的建设情况;和SrPV分别对应于光伏的实际功率和额定容量,表示第i个节点光伏的建设情况;表示第i个节点电压的平方,Uupp和Ulow表示节点电压的上、下限,为直流节点的集合;为第n条线路的额定电流,SrL为线路的额定容量,为切负荷功率;ΩWD表示风电节点的集合,ΩPV表示光伏节点的集合;LDi,v,t表示直流负荷的大小。
进一步的,本发明方法中,所述步骤3)中,求解下层优化模型的方式为:遍历步骤2)中下层优化模型中所有可能的故障情况,对于每一种故障情况,采用CPLEX求解器在MATLAB环境中进行求解,在所有N-1故障中取出最大值,如果该最大值大于等于阈值0.01,表明存在弃风弃光现象,则舍弃当前直流网架结构,重新生成新的直流网架结构并返回步骤2);如果小于阈值,则保留当前直流网架结构并执行步骤4)。
进一步的,本发明方法中,所述步骤4)中构建的上层交流系统规划模型包括:
5.1目标函数
其中,变量上方的“—”和“~”是分别对应于直流和交流变量的标志符;下标中v代表第v个负荷等级,t代表第t个阶段;表示第n条线路的建设情况,表示第i个节点VSC的建设情况,表示第i个节点SVG的建设情况,表示第i个节点变电站的建设情况,且以上变量均为0-1变量;和分别对应于交流系统中线路、VSC、SVG和变电站的单位容量建设成本系数;和分别对应于交流系统中线路、VSC、风电和光伏的单位容量建设成本系数;表示候选交流线路集合,表示候选VSC节点集合,ΩD表示负荷节点的集合,ΩS表示候选变电站节点的集合;和为负荷节点的切负荷功率和切负荷惩罚系数;和分别表示变电站节点从上级电网购电的功率和购电成本系数,LDv表示第v个负荷等级的持续时间;和分别表示风电场的出力和运行费用,和分别表示光伏电站的出力和运行费用;和分别表示线路功率的流动方向,表示线路运行成本系数;RRL,RRVSC,RRSVG,RRS,RRWD和RRPV分别表示线路、VSC、SVG、变电站、风电场和光伏电站的资本回收率,设每一种电气设备的设计使用寿命为LTL,LTVSC,LTSVG,LTS和LTDG,则资本回收率的具体定义如下:
5.2建设约束
5.3交流DistFlow潮流约束
其中,表示由变电站注入节点i的有功功率,表示负荷阻抗角;表示线路n上流过的无功功率,表示由SVG注入节点i的无功功率,表示由VSC注入节点i的无功功率;Xn表示第n条线路的电抗。和分别表示以i节点为首末端节点的线路集合。
5.4直流DistFlow潮流约束
5.5辐射式运行约束
5.6交直流耦合约束
5.7上下限约束
进一步的,本发明方法中,所述步骤5)中求解上层交流系统规划模型的具体方法为:将步骤3)中确立的每一组直流网架结构代入步骤4)构造的上层交流系统规划模型中,采用CPLEX求解器在MATLAB环境中进行求解,将得到的全寿命周期最优成本效益作为对应于每一个直流网架结构染色体的适应度值。判断当前迭代次数是否达到上限。如果尚未达到迭代上限,则对当前种群进行选择、交叉和变异操作,得到一个新的种群,并返回步骤2)进行下一次迭代计算;如果达到迭代上限,则从当前种群中选择适应度最小的直流网架结构及其对应的交流网架配套建设方案,作为最优的交直流混合配电网规划方案输出。
为了交直流混合配电网规划的同时保证N-1安全准则,本发明方法将大量且集中分布在偏远或沿海地区的新能源以直流组网的形式接入直流网络,并通过线路将可再生能源输送到交流配电网当中。规划的目标是在最小化交直流混合配电网的投资和运行费用的同时,保证直流侧网架在发生N-1安全故障的条件下,不出现弃风弃光现象。
有益效果:本发明与现有技术相比,具有以下优点:
1)现有的交直流混合规划模型普遍局限于楼宇、微网等局部应用场景,并且不考虑N-1安全准则。本发明所提的交直流规划模型适用于配网层面,特别是可再生新能源集中分布的偏远地区或沿海地区,并且考虑了N-1安全准则,规划方案更为可靠。
2)现有的交直流潮流模型普遍为传统牛顿-拉夫逊潮流计算的改进,在相关公式中加入直流潮流项。其缺点在于潮流计算公式非凸,包含大量二次型变量,不能转化为可以快速求解的凸优化问题。本发明通过将交流配网潮流计算中的DistFlow模型通过一定简化,引入到直流配电网潮流计算中,保证了模型的凸性,使得配电网的交直流混合潮流可以通过商业软件快速求解。
3)现有的配电网规划模型普遍采用一种优化策略,即要么采用数值解法,要么基于智能算法,但这两种方法对于求解大规模优化问题均不理想。本发明集成了两者的优势,在求解模型的过程中既保证求解精度,又避免出现维数灾难,实现了大规模优化问题的可靠、快速求解。
附图说明:
图1是模型求解的算法流程图;
图2是用于模型测试的24节点系统拓扑图;
图3是以1、3、7和20号节点为直流接入交流的候选节点集的规划结果。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施案例对本发明进行深入地详细说明。应当理解,此处所描述的具体实施案例仅仅用以解释本发明,并不用于限定发明。
如图1所示,本发明提出一种考虑N-1安全准则的交直流混合配电网双层规划方法,包括以下步骤:
1)初始化遗传算法参数,基于0-1编码方式随机构造直流网架结构,直流网架结构数量等于种群数,令当前网架编号k=1;
2)基于鲁棒规划,针对第k个表征直流网架结构的染色体个体,构造针对直流配电网并考虑N-1安全准则的下层优化模型;
3)求解所述步骤2)构造的下层优化模型,检验直流网架结构在最恶劣N-1故障下是否存在弃风弃光,若存在弃风弃光,表明当前直流网架结构无法通过检验,则舍弃当前直流网架结构,基于步骤1)的网架生成方法重新生成当前直流网架结构,并返回步骤2);如果不存在弃风弃光,判断k是否小于种群数,如果是,k=k+1,返回步骤2);否则,进入步骤4)
4)针对当前种群中的每一个直流网架结构,构造与该直流网架结构相配套的上层交流系统规划模型;
5)求解所述步骤4)构造的上层交流系统规划模型,将优化结果作为每一个直流网架结构的适应度值。判断当前迭代次数是否达到遗传算法迭代次数上限,若未达到迭代上限,则对直流网架结构染色体进行选择、交叉、变异操作后,令k=1并返回步骤2),否则,将当前种群中的直流网架结构所对应的适应度最小的交直流网架结构作为结果输出。
步骤1)中,设置遗传算法参数为:种群数量为25个,迭代次数为100次,交叉概率为0.8,变异概率为0.6。根据规划实际需求确定交直流配电网系统拓扑结构,构造交流线路、直流线路、变电站、VSC、风电和光伏规划候选集。0-1编码方式为0代表不建设,1代表建设,构建表征直流网架结构的染色体种群。
步骤2)中构造的下层优化模型包括:
3.1目标函数
其中和分别对应第i个节点在第v个负荷等级和第t个阶段下的弃风和弃光功率;CCW和CCP分别对应弃风和弃光的惩罚系数;WD为风电场候选节点集合,PV为光伏电站候选节点集合;bn,t为0-1变量,表示第n条线路是否发生故障;ΩCTG为N-1故障点集合。
3.2约束条件
其中,变量上方的横线“—”是直流系统变量的标识符,变量上方的尖号“^”是通过遗传算法给定的直流拓扑变量的标识符;下标中v代表第v个负荷等级,t代表第t个阶段;表示第n条线路上流过的功率,表示第i个节点流入VSC的功率,和表示第i个节点风电和光伏的功率,ηVSC为VSC的能量转换效率系数,表示与第i个节点相连的线路的集合,表示直流系统中VSC节点的集合;和分别表示第n条线路的首末端节点电压的平方,Rn表示第n条线路的电阻,表示第n条线路上流过的电流的平方,M为一个非常大的正整数,表示第n条线路的建设情况,为直流候选线路的集合;和SrWD分别对应于风电的实际功率和额定容量,表示第i个节点风电的建设情况;和SrPV分别对应于光伏的实际功率和额定容量,表示第i个节点光伏的建设情况;表示第i个节点电压的平方,Uupp和Ulow表示节点电压的上、下限,为直流节点的集合;为第n条线路的额定电流,SrL为线路的额定容量,为切负荷功率;ΩWD表示风电节点的集合,ΩPV表示光伏节点的集合;LDi,v,t表示直流负荷的大小。
步骤3)中,求解下层优化模型的方式为:遍历步骤2)中下层优化模型中所有可能的故障情况,对于每一种故障情况,采用CPLEX求解器在MATLAB环境中进行求解,在所有N-1故障中取出最大值,如果该最大值大于等于阈值0.01,表明存在弃风弃光现象,则舍弃当前直流网架结构,重新生成新的直流网架结构并返回步骤2);如果小于阈值,则保留当前直流网架结构并执行步骤4)。
步骤4)中构建的上层交流系统规划模型包括:
5.1目标函数
其中,变量上方的“—”和“~”是分别对应于直流和交流变量的标志符;下标中v代表第v个负荷等级,t代表第t个阶段;表示第n条线路的建设情况,表示第i个节点VSC的建设情况,表示第i个节点SVG的建设情况,表示第i个节点变电站的建设情况,且以上变量均为0-1变量;和分别对应于交流系统中线路、VSC、SVG和变电站的单位容量建设成本系数;和分别对应于交流系统中线路、VSC、风电和光伏的单位容量建设成本系数;表示候选交流线路集合,表示候选VSC节点集合,ΩD表示负荷节点的集合,ΩS表示候选变电站节点的集合;和为负荷节点的切负荷功率和切负荷惩罚系数;和分别表示变电站节点从上级电网购电的功率和购电成本系数,LDv表示第v个负荷等级的持续时间;和分别表示风电场的出力和运行费用,和分别表示光伏电站的出力和运行费用;和分别表示线路功率的流动方向,表示线路运行成本系数;RRL,RRVSC,RRSVG,RRS,RRWD和RRPV分别表示线路、VSC、SVG、变电站、风电场和光伏电站的资本回收率,设每一种电气设备的设计使用寿命为LTL,LTVSC,LTSVG,LTS和LTDG,则资本回收率的具体定义如下:
5.2建设约束
5.3交流DistFlow潮流约束
其中,表示由变电站注入节点i的有功功率,表示负荷阻抗角;表示线路n上流过的无功功率,表示由SVG注入节点i的无功功率,表示由VSC注入节点i的无功功率;Xn表示第n条线路的电抗。和分别表示以i节点为首末端节点的线路集合。
5.4直流DistFlow潮流约束
5.5辐射式运行约束
5.6交直流耦合约束
5.7上下限约束
步骤5)中求解上层交流系统规划模型的具体方法为:将步骤3)中确立的每一组直流网架结构代入步骤4)构造的上层交流系统规划模型中,采用CPLEX求解器在MATLAB环境中进行求解,将得到的全寿命周期最优成本效益作为对应于每一个直流网架结构染色体的适应度值。判断当前迭代次数是否达到上限。如果尚未达到迭代上限,则对当前种群进行选择、交叉和变异操作,得到一个新的种群,并返回步骤2)进行下一次迭代计算;如果达到迭代上限,则从当前种群中选择适应度最小的直流网架结构及其对应的交流网架配套建设方案,作为最优的交直流混合配电网规划方案输出。
下面举一实例。以图2所示的系统为例,其系统参数如表1至表5所示。
表1线路参数
表2负荷参数(kW)
表3变电站参数
表4分布式能源参数
表5VSC参数
除表中参数外,交流线路的单位长度电阻为0.407Ω/km,电抗为0.380Ω/km,容量为6280kVA,投资成本为19140$/km;直流线路的单位长度电阻为0.247Ω/km,容量为5000kVA,投资成本为17500$/km。交直流线路的维护成本均为450$/年。负荷按照峰值负荷的70%,83%和100%分别划分年负荷持续曲线为2000小时/年,5760小时/年和1000小时/年。负荷阻抗角取0.9,从上级电网购电成本为0.07$/kVAh。SVG的投资成本为每台11000美元,容量为500kVA,上限为4台。此外,年利率取7.1%,所有设备的寿命均为20年。
基于所提发明步骤进行建模和优化,得到具体建设方案如图3所示,相关的建设和运行费用如表6所示。
表6成本统计表(百万美元)
在表6中,直流网架结构的建设成本占了总建设成本的绝大部分,这是因为考虑了多种新能源的建设费用以及交直流网架结构耦合处VSC的建设费用。而交流系统的运行成本则占了总运行成本的绝大多数,这主要是因为新能源的购电费用和运行费用几乎可以忽略不计,而交流系统则需要计及从上级公共电网购电的费用。采用了本文方法得到的规划方案总成本为1.5393亿美元。
图3给出了针对当前24节点测试系统的三阶段规划方案。可以看出,在不同规划方案下,直流网架结构始终满足N-1安全准则,并且交流系统的拓扑结构始终满足辐射式运行约束。在阶段1,位于26到29号节点的新能源电站首先建设,并通过两条直流线路27-1和28-3接入交流系统。交流系统中尚未建设新的变电站,但在1、6和7号节点处装设SVG以提供无功补偿和提高末端电压。在阶段2,位于25和31号节点的新能源电站投入运行,同时随着负荷功率的提高,原有变电站不能满足负荷的增长需求,因此新建位于23号节点的变电站以提供额外的功率,并且在17号节点新增了一台SVG。在阶段3,所有新能源电站均建成。考虑到新能源的总容量比阶段1和2都要大,原有的两条交直流传输线路不能满足N-1安全准则,因此新增了线路30-20以保障系统运行可靠性。同时,位于24号节点的变电站也在这一阶段建成,为18和20号节点供电。
Claims (7)
1.一种考虑N-1安全准则的交直流混合配电网双层规划方法,其特征在于,该方法包括以下步骤:
1)初始化遗传算法参数,基于0-1编码方式随机构造直流网架结构,网架结构的数量等于种群数,令当前网架编号k=1;
2)基于鲁棒规划,针对第k个表征直流网架结构的染色体个体,构造针对直流配电网并考虑N-1安全准则的下层优化模型;
3)求解所述步骤2)构造的下层优化模型,检验直流网架结构在最恶劣N-1故障下是否存在弃风弃光,若存在弃风弃光,表明当前直流网架结构无法通过检验,则舍弃当前直流网架结构,基于0-1编码方式随机重新生成当前直流网架结构,并返回步骤2);如果不存在弃风弃光,则判断k是否小于种群数,如果是,则令k=k+1,返回步骤2),否则,进入步骤4)
4)针对当前种群中的每一个直流网架结构,构造与该直流网架结构相配套的上层交流系统规划模型;
5)求解所述步骤4)构造的上层交流系统规划模型,将得到的优化结果作为每一个直流网架结构的适应度值,判断当前迭代次数是否达到遗传算法迭代次数上限,若未达到迭代上限,则对直流网架结构的染色体进行选择、交叉、变异操作后,令k=1并返回步骤2),否则,将当前种群中的直流网架结构所对应的适应度最小的交直流网架结构作为结果输出。
2.根据权利要求1所述的考虑N-1安全准则的交直流混合配电网双层规划方法,其特征在于,所述步骤1)中,初始化的遗传算法参数为:种群数量为25个,迭代次数为100次,交叉概率为0.8,变异概率为0.6,根据规划实际需求确定交直流配电网系统拓扑结构,构造交流线路、直流线路、变电站、VSC、风电和光伏规划候选集,0-1编码方式为0代表不建设,1代表建设,构建表征直流网架结构的染色体种群。
3.根据权利要求1所述的考虑N-1安全准则的交直流混合配电网双层规划方法,其特征在于,所述步骤2)中构造的下层优化模型包括:
3.1目标函数
其中,和分别对应第i个节点在第v个负荷等级和第t个阶段下的弃风和弃光功率;CCW和CCP分别对应弃风和弃光的惩罚系数;WD为风电场候选节点集合,PV为光伏电站候选节点集合;bn,t为0-1变量,表示第n条线路是否发生故障;ΩCTG为N-1故障点集合;
3.2约束条件
其中,变量上方的横线“—”是直流系统变量的标识符,变量上方的尖号“^”是通过遗传算法给定的直流拓扑变量的标识符;下标中v代表第v个负荷等级,t代表第t个阶段;表示第n条线路上流过的功率,表示第i个节点流入VSC的功率,和表示第i个节点风电和光伏的功率,ηVSC为VSC的能量转换效率系数,表示与第i个节点相连的线路的集合,表示直流系统中VSC节点的集合;和分别表示第n条线路的首末端节点电压的平方,Rn表示第n条线路的电阻,表示第n条线路上流过的电流的平方,M为一个非常大的正整数,表示第n条线路的建设情况,为直流候选线路的集合;和SrWD分别对应于风电的实际功率和额定容量,表示第i个节点风电的建设情况;和SrPV分别对应于光伏的实际功率和额定容量,表示第i个节点光伏的建设情况;表示第i个节点电压的平方,Uupp和Ulow表示节点电压的上、下限,为直流节点的集合;为第n条线路的额定电流,SrL为线路的额定容量,为切负荷功率;ΩWD表示风电节点的集合,ΩPV表示光伏节点的集合;LDi,v,t表示直流负荷的大小。
4.根据权利要求1、2或3所述的考虑N-1安全准则的交直流混合配电网双层规划方法,其特征在于,所述步骤3)中,求解下层优化模型的方式为:遍历步骤2)中下层优化模型中所有可能的故障情况,对于每一种故障情况,采用CPLEX求解器在MATLAB环境中进行求解,在所有N-1故障中取出最大值,如果该最大值大于或等于阈值0.01,则存在弃风弃光,否则不存在弃风弃光。
5.根据权利要求1、2或3所述的考虑N-1安全准则的交直流混合配电网双层规划方法,其特征在于,所述步骤4)中构建的上层交流系统规划模型包括:
5.1目标函数
其中,变量上方的“—”和“~”是分别对应于直流和交流变量的标志符;下标中v代表第v个负荷等级,t代表第t个阶段;表示第n条线路的建设情况,表示第i个节点VSC的建设情况,表示第i个节点SVG的建设情况,表示第i个节点变电站的建设情况,且以上变量均为0-1变量;和分别对应于交流系统中线路、VSC、SVG和变电站的单位容量建设成本系数;和分别对应于交流系统中线路、VSC、风电和光伏的单位容量建设成本系数;表示候选交流线路集合,表示候选VSC节点集合,ΩD表示负荷节点的集合,ΩS表示候选变电站节点的集合;和为负荷节点的切负荷功率和切负荷惩罚系数;和分别表示变电站节点从上级电网购电的功率和购电成本系数,LDv表示第v个负荷等级的持续时间;和分别表示风电场的出力和运行费用,和分别表示光伏电站的出力和运行费用;和分别表示线路功率的流动方向,表示线路运行成本系数;RRL,RRVSC,RRSVG,RRS,RRWD和RRPV分别表示线路、VSC、SVG、变电站、风电场和光伏电站的资本回收率,设每一种电气设备的设计使用寿命为LTL,LTVSC,LTSVG,LTS和LTDG,则资本回收率的具体定义如下:
5.2建设约束
5.3交流DistFlow潮流约束
其中,表示由变电站注入节点i的有功功率,表示负荷阻抗角;表示线路n上流过的无功功率,表示由SVG注入节点i的无功功率,表示由VSC注入节点i的无功功率;Xn表示第n条线路的电抗,和分别表示以i节点为首末端节点的线路集合;
5.4直流DistFlow潮流约束
5.5辐射式运行约束
5.6交直流耦合约束
5.7上下限约束
6.根据权利要求1、2或3所述的考虑N-1安全准则的交直流混合配电网双层规划方法,其特征在于,所述步骤5)中求解上层交流系统规划模型的具体方法为:将步骤3)中构造的下层优化模型中的每一组直流网架结构代入步骤4)构造的上层交流系统规划模型中,采用CPLEX求解器在MATLAB环境中进行求解,将得到的全寿命周期最优成本效益作为对应于每一个直流网架结构染色体的适应度值。
7.根据权利要求1、2或3所述的考虑N-1安全准则的交直流混合配电网双层规划方法,其特征在于,所述步骤5)中,最优的交直流网架结构为当前种群中适应度最小的直流网架结构及其对应的交流网架配套建设方案。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811085672.1A CN109002938B (zh) | 2018-09-17 | 2018-09-17 | 考虑n-1安全准则的交直流混合配电网双层规划方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811085672.1A CN109002938B (zh) | 2018-09-17 | 2018-09-17 | 考虑n-1安全准则的交直流混合配电网双层规划方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109002938A CN109002938A (zh) | 2018-12-14 |
CN109002938B true CN109002938B (zh) | 2021-06-01 |
Family
ID=64592582
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811085672.1A Active CN109002938B (zh) | 2018-09-17 | 2018-09-17 | 考虑n-1安全准则的交直流混合配电网双层规划方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109002938B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110350605B (zh) * | 2019-07-08 | 2020-12-29 | 华中科技大学 | 一种多阶段交直流混联配电网直流升级改造的规划方法 |
CN113036821B (zh) * | 2021-04-20 | 2023-08-01 | 国网四川省电力公司经济技术研究院 | 构建巨量直流送出的交直流协调骨干网架的方法及系统 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103049799A (zh) * | 2012-12-10 | 2013-04-17 | 河海大学 | 一种基于多目标优化的电网规划设计方法 |
CN104715294A (zh) * | 2015-04-02 | 2015-06-17 | 国家电网公司 | 一种配电网规划方案评价数据处理系统 |
CN105574632A (zh) * | 2014-10-14 | 2016-05-11 | 国家电网公司 | 一种交直流混合城市配电网综合效益评估方法 |
CN105740980A (zh) * | 2016-01-29 | 2016-07-06 | 国网福建省电力有限公司 | 一种含柔性直流输电系统的交直流混合电网优化规划方法 |
CN106815657A (zh) * | 2017-01-05 | 2017-06-09 | 国网福建省电力有限公司 | 一种考虑时序性和可靠性的配电网双层规划方法 |
CN107358337A (zh) * | 2017-06-08 | 2017-11-17 | 上海电力学院 | 一种基于网络重构的有源配电网规划方法 |
CN108122068A (zh) * | 2017-12-05 | 2018-06-05 | 中国电力科学研究院有限公司 | 一种配电网风险规划方法及系统 |
-
2018
- 2018-09-17 CN CN201811085672.1A patent/CN109002938B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103049799A (zh) * | 2012-12-10 | 2013-04-17 | 河海大学 | 一种基于多目标优化的电网规划设计方法 |
CN105574632A (zh) * | 2014-10-14 | 2016-05-11 | 国家电网公司 | 一种交直流混合城市配电网综合效益评估方法 |
CN104715294A (zh) * | 2015-04-02 | 2015-06-17 | 国家电网公司 | 一种配电网规划方案评价数据处理系统 |
CN105740980A (zh) * | 2016-01-29 | 2016-07-06 | 国网福建省电力有限公司 | 一种含柔性直流输电系统的交直流混合电网优化规划方法 |
CN106815657A (zh) * | 2017-01-05 | 2017-06-09 | 国网福建省电力有限公司 | 一种考虑时序性和可靠性的配电网双层规划方法 |
CN107358337A (zh) * | 2017-06-08 | 2017-11-17 | 上海电力学院 | 一种基于网络重构的有源配电网规划方法 |
CN108122068A (zh) * | 2017-12-05 | 2018-06-05 | 中国电力科学研究院有限公司 | 一种配电网风险规划方法及系统 |
Also Published As
Publication number | Publication date |
---|---|
CN109002938A (zh) | 2018-12-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106503839B (zh) | 一种海上风电场环形集电网络分层规划方法 | |
CN106655207A (zh) | 基于多数据分析的配电网无功优化系统及方法 | |
CN104537580B (zh) | 基于K‑means聚类的配电网网架的构建方法 | |
CN109598377B (zh) | 一种基于故障约束的交直流混合配电网鲁棒规划方法 | |
CN105512472B (zh) | 大型风电基地功率汇集系统拓扑结构分层优化设计方法 | |
CN108808715B (zh) | 考虑直流网络故障功率的多端柔直系统静态安全分析方法 | |
CN102930078B (zh) | 一种孤立微网分布式电源容量与布点优化方法 | |
CN111563691B (zh) | 一种接入新能源的交直流混合配电网性能评估方法 | |
CN105184669A (zh) | 一种基于节点集GN分裂算法的220kV城市环网分区方法 | |
CN105760971A (zh) | 一种基于可靠性比较分析的城市电网结构优化方法 | |
CN109002938B (zh) | 考虑n-1安全准则的交直流混合配电网双层规划方法 | |
CN114202186A (zh) | 一种海上风电场并网输送系统经济评估性方法 | |
CN110994697B (zh) | 含光储联合体的交直流配电网优化运行控制方法及系统 | |
CN105514987A (zh) | 10kV配电线路分布式光伏定容方法和系统 | |
CN107370182B (zh) | 一种主动配电网三元规划体系分布式电源接入规划方法 | |
CN111105089B (zh) | 一种考虑主配网协调的城市电网规划方法 | |
CN108988336B (zh) | 具有嵌套式微电网的充电桩系统的优化规划方法 | |
CN116681149A (zh) | 基于全生命周期的海上风电输电方式经济技术分析方法 | |
CN115495989A (zh) | 一种考虑风特性的海上风电场集电系统拓扑结构设计方法 | |
Wanchun et al. | Technical and economic analysis of power collection and transmission schemes adapted to large-scale offshore wind farm | |
CN105762832A (zh) | 配网分布式能源接入极限优化方法 | |
Guangyuan et al. | A Diversify Distribution Power Grid Planning Method Based on the Features of Areas in Smart Cities | |
CN105990829B (zh) | 一种lcc-vsc混合直流电网测试模型 | |
CN103473615A (zh) | 一种特高压交直流落点自动选择方法及装置 | |
CN113162101B (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 |