CN115345350A - 一种交直流混合配电网规划方法、介质及系统 - Google Patents
一种交直流混合配电网规划方法、介质及系统 Download PDFInfo
- Publication number
- CN115345350A CN115345350A CN202210825774.2A CN202210825774A CN115345350A CN 115345350 A CN115345350 A CN 115345350A CN 202210825774 A CN202210825774 A CN 202210825774A CN 115345350 A CN115345350 A CN 115345350A
- Authority
- CN
- China
- Prior art keywords
- representing
- power
- energy storage
- node
- time
- 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
- 238000000034 method Methods 0.000 title claims abstract description 29
- 230000006870 function Effects 0.000 claims abstract description 40
- 238000004146 energy storage Methods 0.000 claims description 68
- 230000004044 response Effects 0.000 claims description 41
- 238000009434 installation Methods 0.000 claims description 16
- 239000002245 particle Substances 0.000 claims description 14
- 238000004422 calculation algorithm Methods 0.000 claims description 13
- 230000005540 biological transmission Effects 0.000 claims description 12
- 238000007599 discharging Methods 0.000 claims description 12
- 230000003828 downregulation Effects 0.000 claims description 11
- 230000003827 upregulation Effects 0.000 claims description 11
- 230000007812 deficiency Effects 0.000 claims description 8
- 238000004590 computer program Methods 0.000 claims description 6
- 230000008859 change Effects 0.000 claims description 4
- 230000005611 electricity Effects 0.000 claims description 4
- 230000008569 process Effects 0.000 claims description 4
- 239000000126 substance Substances 0.000 claims description 4
- 239000003344 environmental pollutant Substances 0.000 claims description 3
- 239000000446 fuel Substances 0.000 claims description 3
- 230000003993 interaction Effects 0.000 claims description 3
- 238000012423 maintenance Methods 0.000 claims description 3
- 231100000719 pollutant Toxicity 0.000 claims description 3
- 238000005067 remediation Methods 0.000 claims description 3
- 230000009194 climbing Effects 0.000 claims description 2
- 238000011002 quantification Methods 0.000 claims description 2
- 230000001105 regulatory effect Effects 0.000 claims 1
- 239000007789 gas Substances 0.000 description 25
- 238000005457 optimization Methods 0.000 description 6
- 238000010586 diagram Methods 0.000 description 4
- 238000011156 evaluation Methods 0.000 description 4
- 229910052739 hydrogen Inorganic materials 0.000 description 4
- 239000001257 hydrogen Substances 0.000 description 4
- UFHFLCQGNIYNRP-UHFFFAOYSA-N Hydrogen Chemical compound [H][H] UFHFLCQGNIYNRP-UHFFFAOYSA-N 0.000 description 2
- 230000033228 biological regulation Effects 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 150000002431 hydrogen Chemical class 0.000 description 2
- 238000006467 substitution reaction Methods 0.000 description 2
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 1
- 241001098636 Trichogramma alpha Species 0.000 description 1
- 230000002457 bidirectional effect Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000003064 k means clustering Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- MYVIATVLJGTBFV-UHFFFAOYSA-M thiamine(1+) chloride Chemical compound [Cl-].CC1=C(CCO)SC=[N+]1CC1=CN=C(C)N=C1N MYVIATVLJGTBFV-UHFFFAOYSA-M 0.000 description 1
- 238000012795 verification 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
- 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
- H02J4/00—Circuit arrangements for mains or distribution networks not specified as ac or dc
Landscapes
- Business, Economics & Management (AREA)
- Engineering & Computer Science (AREA)
- Economics (AREA)
- Human Resources & Organizations (AREA)
- Strategic Management (AREA)
- Tourism & Hospitality (AREA)
- Health & Medical Sciences (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Marketing (AREA)
- General Business, Economics & Management (AREA)
- Physics & Mathematics (AREA)
- Game Theory and Decision Science (AREA)
- Quality & Reliability (AREA)
- Operations Research (AREA)
- Entrepreneurship & Innovation (AREA)
- Development Economics (AREA)
- Power Engineering (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所示，直流配电网和交流配电网通过电压源换流器VSC(Voltage Source Converter)连接构成交直流混合配电网。直流配电网中接有各类直流负荷、光伏机组、储能装置等。交流配电网中接有交流负荷、风电机组、光伏机组等。通过高效调度各种泛在灵活性资源可实现交直流混合配电网的安全经济运行和高效决策。
含光伏和风电等分布式电源、交直流负荷以及储能在内的全要素交直流配电网络示意图如图2所示，其特点在于电力电子装置是其中的核心环节，因此交直流混合配电网具有高度电力电子化特征，可以大幅提升系统的稳定性、可控性和灵活性；另一方面，其兼容交直流两种电能形式，在为当前交流负荷供电的同时便于高比例可再生能源的接入，降低了系统复杂度，能够提高能量传输效率。
配电网的灵活性应具备应对不同时间尺度下单个或多个单元的波动性和随机不确定性，包括风机和光伏出力的不确定性以及电动汽车等新型负荷需求的不确定性。因此，配电网的灵活性具有多时间尺度、双向调节、多分布空间的多维特征，电力系统的供需平衡可以满足不同的时间尺度的变化。一方面，电力系统内的灵活性资源的可调度灵活性容量是随着系统的时序变化而发生变化，在系统运行的不同时刻，灵活性资源的参数也实时发生变化。另一方面，电力系统灵活性资源的响应时间和响应速度是不同的，和灵活性需求之间存在明显的差异。因此，在配电网规划中基于供需平衡考虑灵活性具有重要意义。
发明内容
本发明实施例提供一种交直流混合配电网规划方法、介质及系统，以解决现有技术进行交直流混合配电网规划未考虑灵活性资源的问题。
第一方面，提供一种交直流混合配电网规划方法，包括：
建立由下层运行层的目标函数和上层规划层的目标函数组成的交直流混合配电网规划的双层模型；
分别确定所述运行层和所述规划层的目标函数的约束条件；
分别在所述运行层和所述规划层的约束条件下，求解所述双层模型，输出最优的灵活性资源的位置和容量，使按照所述最优的灵活性资源的位置和容量规划所述交直流混合配电网；
其中，所述运行层的目标函数为：
minF1＝CF+COM+CENV+CDR+CG+Closs+CDRE；
F1表示系统年运行成本，CF表示燃料成本，COM表示灵活资源运维成本，CENV表示污染物治理成本，CDR表示需求响应负荷补偿成本，CG表示主网购电成本，Closs表示网损成本，CDRE表示弃风弃光惩罚费用；
其中，所述规划层的目标函数包括：
系统综合经济性目标函数：
minF2＝F1+Cinv；
F2表示系统年综合成本，Cinv表示灵活资源投资成本，
Cinv,MT表示微型燃气轮机单位容量投资成本，Cinv,ESS表示电力系统储能设备单位容量投资成本，Pi,MT表示节点i的微型燃气轮机安装容量，Ei,ESS表示节点i的电力系统储能设备安装容量，NMT表示微型燃气轮机待安装数量，NESS表示电力系统储能设备待安装数量，n表示灵活资源经济使用年限，r表示折现率；以及，
系统综合安全性目标函数：
CSCS表示系统综合安全性指标，CLSS,s,t,l表示场景s发生的总时间τ内第l条输电线路故障的负荷安全供电率指标，Ns表示场景数量，Nl表示输电线路数量，H表示安全分析时长；节点i为IEEE39节点中的节点，取值范围为1-39；
所述灵活性资源的位置和容量为节点i的微型燃气轮机安装容量和节点i的电力系统储能设备安装容量。
第二方面，提供一种计算机可读存储介质，所述计算机可读存储介质上存储有计算机程序指令；所述计算机程序指令被处理器执行时实现如上述第一方面实施例所述的交直流混合配电网规划方法。
第三方面，提供一种交直流混合配电网规划系统，包括：如上述第二方面实施例所述的计算机可读存储介质。
这样，本发明实施例，通过考虑系统灵活性供需关系，建立灵活资源运行-规划联合优化的双层模型。引入灵活性不足率和系统综合安全性指标进行评估，通过综合配置多种灵活资源能有效降低网损和系统弃风弃光率，提高系统经济性和灵活性资源的调节能力，从而提高交直流混合配电网的灵活性。
附图说明
为了更清楚地说明本发明实施例的技术方案，下面将对本发明实施例的描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。
图1是交直流混合配电网拓扑结构；
图2是全要素交直流配电网络示意图；
图3是本发明实施例的交直流混合配电网规划方法的流程图；
图4是本发明实施例的双层模型互动关系示意图；
图5是不同技术方案下网损与弃风弃光率关系示意图。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获取的所有其他实施例，都属于本发明保护的范围。
电力系统的核心是电力电量平衡，即根据规划年的负荷需求确定电力系统规划年的装机容量及需新增的容量，保证系统的可靠性及火电机组合理的利用小时数。传统规划仅仅考虑电力电量平衡，通过新建发电厂满足负荷增长及充裕度需求，没有考虑到灵活性的需求。然而随着可再生能源比例的增加，系统的需求不仅仅是单纯的负荷，灵活性需求也逐步凸显，并且已经远远超过传统电力系统的灵活性需求，有必要对灵活性的平衡与匹配进行专项分析。灵活性平衡已经成为系统安全稳定运行的关键和核心。
基于此，本发明实施例公开了一种交直流混合配电网规划方法。该方法考虑灵活性供需平衡。本发明实施例以微型燃气轮机(microturbine，MT)、电力系统储能设备(energy storage system，ESS)和需求响应(demand response，DR)负荷为灵活性资源，这三种可控资源都具备优良的灵活性调节能力。
如图3所示，该方法包括如下的步骤：
步骤S1：建立由下层运行层的目标函数和上层规划层的目标函数组成的交直流混合配电网规划的双层模型。
具体的，运行层的目标函数采用如下的式(1)表示：
(1)minF1＝CF+COM+CENV+CDR+CG+Closs+CDRE。
其中，F1表示系统年运行成本。CF表示燃料成本。COM表示灵活资源运维成本。CENV表示污染物治理成本。CDR表示需求响应负荷补偿成本。CG表示主网购电成本。Closs表示网损成本。CDRE表示弃风弃光惩罚费用。
具体的，规划层的目标函数包括：
一、系统综合经济性目标函数，具体采用如下的式(2)表示：
(2)minF2＝F1+Cinv。
其中，F2表示系统年综合成本。Cinv表示灵活资源投资成本。
具体的，Cinv采用如下的式(3)计算：
(3)
其中，Cinv,MT表示微型燃气轮机单位容量投资成本。Cinv,ESS表示电力系统储能设备单位容量投资成本。Pi,MT表示节点i的微型燃气轮机安装容量，Ei,ESS表示节点i的电力系统储能设备安装容量。NMT表示微型燃气轮机待安装数量。NESS表示电力系统储能设备待安装数量。n表示灵活资源经济使用年限。r表示折现率。
二、系统综合安全性目标函数，具体采用如下的式(4)表示：
(4)
其中，CSCS表示系统综合安全性指标。CLSS,s,t,l表示场景s发生的总时间τ内第l条输电线路故障的负荷安全供电率指标。Ns表示场景数量。Nl表示输电线路数量。H表示安全分析时长。在一具体的实施例中，H＝4h。
本发明实施例的规划层还引入了上述的系统综合安全性指标对系统安全分析各场景对应总时间τ内N-1故障下的负荷安全供电能力进行评判，以评估系统综合安全性。该指标描述了系统在发生支路故障下，由内部电源形成孤岛而维持系统负荷持续供电能力。
步骤S2：分别确定运行层和规划层的目标函数的约束条件。
具体的，运行层的目标函数的约束条件包括：
一、需求响应负荷变化量限值约束如式(6)所示：
(6)
其中，Pload,t表示t时刻用户负荷。αmin表示响应时段内移除负荷响应系数的下限。αmax表示响应时段内移除负荷响应系数的上限。βmin表示响应时段内增加负荷响应系数的下限。βmax表示响应时段内增加负荷响应系数的上限，PL,cut,t表示t时刻需求响应负荷移除的负荷。PL,add,t表示t时刻需求响应负荷增加的负荷。
在一具体的实施例中，αmin＝0，αmax＝0.5，βmin＝0，βmax＝0.5。
二、灵活运行指标约束如式(7)所示：
(7)floss<σ。
其中，floss表示系统灵活不足率。σ表示系统灵活性不足率上限。
具体的，系统灵活不足率floss的算式如式(8)所示：
(8)
其中，ξ表示权重系数。
表示系统上调灵活性不足率。
表示系统下调灵活性不足率。
表示所有场景内系统向上供给的和。
表示所有场景内系统向下供给的和。
若系统灵活性供应大于灵活性需求，则灵活性充足；否则灵活性不足。为了更直观量化灵活性供需关系，本发明实施例引入灵活性不足率作为系统灵活性评价指标。式(8)即为将系统上调灵活性不足率和下调灵活性不足率统一成系统灵活性不足率。具体的，系统上调灵活性不足率和系统下调灵活性不足率的算式如式(9)所示：
(9)
其中，ps表示场景s的发生概率。τ表示场景s发生的总时间。Ft D,down表示所有场景内向下灵活性需求之和。Ft D,up表示所有场景内向上灵活性需求之和。
表示t时刻微型燃气轮机上调灵活性供应。
表示t时刻微型燃气轮机下调灵活性供应。
表示t时刻电力系统储能设备上调灵活性供应。
表示t时刻电力系统储能设备下调灵活性供应。
表示t时刻需求响应负荷上调灵活性供应。
表示t时刻需求响应负荷下调灵活性供应。
具体的，微型燃气轮机上调灵活性供应和微型燃气轮机下调灵活性供应的算式如式(10)所示：
(10)
其中，Rup表示微型燃气轮机上调节功率速率。Rdown表示微型燃气轮机下调节功率速率。
表示微型燃气轮机最大输出功率。
表示微型燃气轮机最小输出功率。PMT,t表示t时刻微型燃气轮机的出力。ΔT表示时间间隔。
具体的，电力系统储能设备上调灵活性供应和电力系统储能设备下调灵活性供应的算式如式(11)所示：
(11)
其中，Pdis,max表示放电最大功率。Pch,max表示充电最大功率。EESS表示电力系统储能设备额定电量。ηdis表示电力系统储能设备放电效率。ηch表示电力系统储能设备充电效率。SOCmax表示电力系统储能设备荷电状态最大值。SOCmin表示电力系统储能设备荷电状态最小值。SOCt表示t时刻电力系统储能设备荷电状态。
具体的，需求响应负荷上调灵活性供应和需求响应负荷下调灵活性供应的算式如式(12)所示：
(12)
其中，FL,cut,t表示t时刻需求响应负荷移除的用电功率，FL,add,t表示t时刻需求响应负荷增加的用电功率。
三、灵活性资源出力约束
其中，灵活性资源出力约束包括电力系统储能设备和微型燃气轮机的出力量化交直流混合配电网灵活性供给能力的约束，如式(13)所示：
(13)
其中，I为电力系统储能设备和微型燃气轮机的节点集合，
分别表示t时刻上、下调灵活性供给；
分别为t时刻灵活性资源的上、下限出力，Pi t为t时刻灵活性资源实际出力，
分别为各灵活性资源的爬坡、滑坡速率，Δt表示时间间隔。具体的，规划层的目标函数的约束条件包括：
一、灵活资源容量约束如式(14)所示：
(14)
其中，
表示节点i的微型燃气轮机的最大安装容量。
表示节点i的电力系统储能设备最大安装容量。
二、系统综合安全性约束如式(15)所示：
(15)CSCS≥η。
其中，η表示系统综合安全性限值。
三、交流区域潮流约束
潮流约束是将交流区域原DisFlow潮流模型中的非线性项通过二阶锥松弛转换得到凸化线性潮流约束，如式(16)、(17)、(18)、(19)所示。
(16)
(17)
(18)
(19)
其中，rij、xij分别表示节点i→j之间电阻、电抗；u(j)、v(j)分别表示以节点j为末端和首端的支路节点集合；
分别表示在节点j处注入总有功、无功功率；
分别表示j节点处灵活性资源的有功、无功功率；
分别表示t时刻流过支路i→j的有功、无功功率；
分别表示t时刻主网流向配电网的有功、无功功率；
分别表示支路电流
和节点i电压
的平方；
表示节点j电压
的平方；
分别表示t时刻节点j处常规负荷、可中断负荷；
分别表示t时刻节点j处常规负荷、可中断负荷。
四、直流区域潮流约束
直流区域潮流方程经过线性化后与交流区域类似，如式(20)所示：
(20)
其中，
分别表示t时刻节点i处电力系统储能设备充、放电功率；
表示t时刻节点j处分布式光伏实际出力。
五、电力系统储能设备运行约束，如式(21)、(22)所示：
(21)
(22)
其中，
分别表示t时刻节点i处电力系统储能设备的最大充、放电功率；
分别表示t时刻节约束点i处电力系统储能设备充、放电状态，0表示充电，1表示放电；Ei,t、
分别表示t时刻节点i处电力系统储能设备电量与最大允许电量；Ei,t+1表示t+1时刻节点i处电力系统储能设备电量；ηi,ch、ηi,dis分别表示ESS充、放电效率，
分别表示t时刻节点i处电力系统储能设备的充、放电功率。
六、与主网交互功率约束，如式(23)所示：
(23)
其中，
分别表示变电站节点i注入有功功率上、下限；
分别表示变电站节点i注入无功功率上、下限，
分别表示t时刻主网流向配电网的节点i的有功、无功功率。
七、安全约束，如式(24)所示：
(24)
其中，Ui,max、Ui,min分别表示节点i电压安全运行上、下限；
分别表示t时刻交流支路i→j传输有功、无功功率；
表示t时刻直流支路i→j传输有功功率；
分别表示交、直流支路i→j最大传输容量；Ui表示节点i电压。
八、电力系统储能设备额定功率和额定容量上、下限约束，如式(25)所示：
(25)
其中，PESS,max、PESS,min分别表示电力系统储能设备额定功率上、下限；EESS,max、EESS,min分别表示电力系统储能设备的额定容量的上、下限；PESS,n表示电力系统储能设备功率，EESS,n表示电力系统储能设备的额定容量。
步骤S3：分别在运行层和规划层的约束条件下，求解双层模型，输出最优的灵活性资源的位置和容量，使按照最优的灵活性资源的位置和容量规划交直流混合配电网。
灵活性资源的位置和容量为节点i的微型燃气轮机安装容量和节点i的电力系统储能设备安装容量，即通过本发明实施例的方法最终得到在特定节点规划的微型燃气轮机安装容量和电力系统储能设备安装容量，从而可依此进行规划。
粒子群优化算法(particle swarm optimization，PSO)是一种利用微粒模拟飞鸟捕食行为，不断更新粒子位置和速度，寻找目标最优解的优化算法。该算法因收敛速度快，搜索能力强的特点而受到广泛应用。基于此，本发明实施例求解双层模型的算法为改进粒子群算法，其是一种采用惯性权重因子和学习因子调整的改进粒子群算法，可进一步优化粒子搜索能力，提高运算收敛性。本发明实施例采用的改进粒子群算法为如下的文献公开的改进粒子群算法：徐岩，张荟，孙易洲.基于变异粒子群算法的主动配电网故障恢复策略[J].电力自动化设备，2021，41(12):45-53。
采用改进粒子群算法求解的过程中，将灵活性资源的位置和容量作为改进例子群算法的粒子的位置和速度。
该步骤的具体过程如下：
(1)初始化灵活性资源的位置和容量。
(2)以灵活性资源出力为决策变量，在运行层的目标函数的约束条件下，求解运行层的目标函数，得到系统年运行成本。
其中，采用k-means聚类算法对风光荷年历史数据进行处理，将风光荷随机特性用不同季节不同气候下多个典型场景确定化描述，即得到算例分析中的初始参数，包括交流配电网电压等级、直流配电网电压等级、总有功负荷、总无功负荷、VSC容量、分布式光伏、分布式风电相关数据。这些数据输入算法中以进行相应的求解，得到系统年运行成本及灵活性资源出力。
(3)将系统年运行成本返回到规划层，以灵活性资源的位置和容量为决策变量，在规划层的目标函数的约束条件下，求解规划层的目标函数，得到更新的灵活性资源的位置和容量。
(4)重复上述的步骤(2)和步骤(3)，直到满足收敛条件，输出最优的灵活性资源的位置和容量。
其中，收敛条件为：前后迭代值之差达到迭代精度或迭代次数达到最大次数。迭代次数可根据经验确定。
其中运行层的决策变量为灵活性资源的出力，通过优化决策灵活性资源出力使目标函数达到最优，目标函数最优。
上述的计算过程，运行层以系统灵活性不足率为灵活运行评价指标，以系统年综合运行成本最优为目标函数，对灵活资源的出力进行优化，规划层还以系统年综合成本最优为目标，对灵活资源位置和容量进行全方位综合规划；运行层将优化各时段灵活性资源出力情况得到的年运行成本反馈到规划层，规划层和运行层进行上下互动，共同决策出灵活性资源最优配置方案，互动关系如图4所示。
为了验证本发明实施例的技术方案的微型燃气轮机、电力系统储能设备和需求响应负荷，这三种灵活资源对提升系统性能的有效性，设定以下三种方案进行验证：1)不考虑系统灵活性，即无灵活资源参与，不计σ和η约束；2)灵活资源仅为微型燃气轮机和需求响应负荷的双层模型，不计σ和η约束；3)灵活资源为微型燃气轮机、电力系统储能设备和需求响应负荷的双层模型，不计σ和η约束。对上述方案进行求解，得出三种方案下系统弃风弃光量与网损量关系见图5。
由图5可见，方案2)比方案1)的网损下降20.9％，弃风弃光率下降8.8％；方案3)又比方案2)降低网损23.2％，降低弃风弃光率5.9％。结果表明在含高比例分布式可再生能源(distributed renewable energy，DRE)配电网中配置灵活资源能有效提高DRE利用率，降低网损；且综合利用多种灵活资源时配电网能最大化减少网损和弃风弃光率，提高系统运行水平。
本发明实施例还公开了一种计算机可读存储介质，所述计算机可读存储介质上存储有计算机程序指令；所述计算机程序指令被处理器执行时实现如上述实施例所述的交直流混合配电网规划方法。
本发明实施例还公开了一种交直流混合配电网规划系统，包括：如上述实施例所述的计算机可读存储介质。
综上，本发明实施例，通过考虑系统灵活性供需关系，建立灵活资源运行-规划联合优化的双层模型。引入灵活性不足率和系统综合安全性指标进行评估，通过综合配置多种灵活资源能有效降低网损和系统弃风弃光率，提高系统经济性和灵活性资源的调节能力，从而提高交直流混合配电网的灵活性。
以上所述，仅为本发明的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，可轻易想到变化或替换，都应涵盖在本发明的保护范围之内。因此，本发明的保护范围应以权利要求的保护范围为准。
Claims (10)
1.一种交直流混合配电网规划方法，其特征在于，包括：
建立由下层运行层的目标函数和上层规划层的目标函数组成的交直流混合配电网规划的双层模型；
分别确定所述运行层和所述规划层的目标函数的约束条件；
分别在所述运行层和所述规划层的约束条件下，求解所述双层模型，输出最优的灵活性资源的位置和容量，使按照所述最优的灵活性资源的位置和容量规划所述交直流混合配电网；
其中，所述运行层的目标函数为：
min F1＝CF+COM+CENV+CDR+CG+Closs+CDRE；
F1表示系统年运行成本，CF表示燃料成本，COM表示灵活资源运维成本，CENV表示污染物治理成本，CDR表示需求响应负荷补偿成本，CG表示主网购电成本，Closs表示网损成本，CDRE表示弃风弃光惩罚费用；
其中，所述规划层的目标函数包括：
系统综合经济性目标函数：
min F2＝F1+Cinv；
F2表示系统年综合成本，Cinv表示灵活资源投资成本，
Cinv,MT表示微型燃气轮机单位容量投资成本，Cinv,ESS表示电力系统储能设备单位容量投资成本，Pi,MT表示节点i的微型燃气轮机安装容量，Ei,ESS表示节点i的电力系统储能设备安装容量，NMT表示微型燃气轮机待安装数量，NESS表示电力系统储能设备待安装数量，n表示灵活资源经济使用年限，r表示折现率；以及，
系统综合安全性目标函数：
CSCS表示系统综合安全性指标，CLSS,s,t,l表示场景s发生的总时间τ内第l条输电线路故障的负荷安全供电率指标，Ns表示场景数量，Nl表示输电线路数量，H表示安全分析时长；
所述灵活性资源的位置和容量为节点i的微型燃气轮机安装容量和节点i的电力系统储能设备安装容量。
2.根据权利要求1所述的交直流混合配电网规划方法，其特征在于，所述运行层的目标函数的约束条件包括：
需求响应负荷变化量限值约束：
其中，Pload,t表示t时刻用户负荷，αmin表示响应时段内移除负荷响应系数的下限，αmax表示响应时段内移除负荷响应系数的上限，βmin表示响应时段内增加负荷响应系数的下限，βmax表示响应时段内增加负荷响应系数的上限，PL,cut,t表示t时刻需求响应负荷移除的负荷，PL,add,t表示t时刻需求响应负荷增加的负荷；
灵活运行指标约束：floss<σ，其中，floss表示系统灵活不足率，σ表示系统灵活性不足率上限；
灵活性资源出力约束包括电力系统储能设备和微型燃气轮机的出力量化交直流混合配电网灵活性供给能力的约束：
其中，I为电力系统储能设备和微型燃气轮机的节点集合，
分别表示t时刻上、下调灵活性供给；
分别为t时刻灵活性资源的上、下限出力，Pi t为t时刻灵活性资源实际出力，
分别为各灵活性资源的爬坡、滑坡速率，Δt表示时间间隔。
3.根据权利要求2所述的交直流混合配电网规划方法，其特征在于，所述系统灵活不足率的算式为：
其中，ξ表示权重系数，
表示系统上调灵活性不足率，
表示系统下调灵活性不足率，
表示所有场景内系统向上供给的和，
表示所有场景内系统向下供给的和。
4.根据权利要求3所述的交直流混合配电网规划方法，其特征在于，所述系统上调灵活性不足率和所述系统下调灵活性不足率的算式为：
其中，ps表示场景s的发生概率，τ表示场景s发生的总时间，Ft D,down表示所有场景内向下灵活性需求之和，Ft D,up表示所有场景内向上灵活性需求之和，
表示t时刻微型燃气轮机上调灵活性供应，
表示t时刻微型燃气轮机下调灵活性供应，
表示t时刻电力系统储能设备上调灵活性供应，
表示t时刻电力系统储能设备下调灵活性供应，
表示t时刻需求响应负荷上调灵活性供应，
表示t时刻需求响应负荷下调灵活性供应。
5.根据权利要求4所述的交直流混合配电网规划方法，其特征在于，
所述微型燃气轮机上调灵活性供应和所述微型燃气轮机下调灵活性供应的算式包括：
其中，Rup表示微型燃气轮机上调节功率速率，Rdown表示微型燃气轮机下调节功率速率，
表示微型燃气轮机最大输出功率，
表示微型燃气轮机最小输出功率，PMT,t表示t时刻微型燃气轮机的出力，ΔT表示时间间隔；
所述电力系统储能设备上调灵活性供应和所述电力系统储能设备下调灵活性供应的算式包括：
其中，Pdis,max表示放电最大功率，Pch,max表示充电最大功率，EESS表示电力系统储能设备额定电量，ηdis表示电力系统储能设备放电效率，ηch表示电力系统储能设备充电效率，SOCmax表示电力系统储能设备荷电状态最大值，SOCmin表示电力系统储能设备荷电状态最小值，SOCt表示t时刻电力系统储能设备荷电状态；
所述需求响应负荷上调灵活性供应和需求响应负荷下调灵活性供应的算式包括：
其中，FL,cut,t表示t时刻需求响应负荷移除的用电功率，FL,add,t表示t时刻需求响应负荷增加的用电功率。
6.根据权利要求1所述的交直流混合配电网规划方法，其特征在于，所述规划层的目标函数的约束条件包括：
灵活资源容量约束：
其中，
表示节点i的微型燃气轮机的最大安装容量，
表示节点i的电力系统储能设备最大安装容量；
系统综合安全性约束：CSCS≥η，其中，η表示系统综合安全性限值；
交流区域潮流约束：
其中，rij、xij分别表示节点i→j之间电阻、电抗；u(j)、v(j)分别表示以节点j为末端和首端的支路节点集合；
分别表示在节点j处注入总有功、无功功率；
分别表示j节点处灵活性资源的有功、无功功率；
分别表示t时刻流过支路i→j的有功、无功功率；
分别表示t时刻主网流向配电网的有功、无功功率；
分别表示支路电流
和节点i电压
的平方；
表示节点j电压
的平方；
分别表示t时刻节点j处常规负荷、可中断负荷；
分别表示t时刻节点j处常规负荷、可中断负荷；
直流区域潮流约束：
其中，
分别表示t时刻节点i处电力系统储能设备充、放电功率；
表示t时刻节点j处分布式光伏实际出力；
电力系统储能设备运行约束：
其中，
分别表示t时刻节点i处电力系统储能设备的最大充、放电功率；
分别表示t时刻节约束点i处电力系统储能设备充、放电状态，0表示充电，1表示放电；Ei,t、
分别表示t时刻节点i处电力系统储能设备电量与最大允许电量；Ei,t+1表示t+1时刻节点i处电力系统储能设备电量；ηi,ch、ηi,dis分别表示电力系统储能设备充、放电效率，
分别表示t时刻节点i处电力系统储能设备的充、放电功率；
与主网交互功率约束：
其中，
分别表示变电站节点i注入有功功率上、下限；
分别表示变电站节点i注入无功功率上、下限，
分别表示t时刻主网流向配电网的节点i的有功、无功功率；
安全约束：
其中，Ui,max、Ui,min分别表示节点i电压安全运行上、下限；
分别表示t时刻交流支路i→j传输有功、无功功率；
表示t时刻直流支路i→j传输有功功率；
分别表示交、直流支路i→j最大传输容量；Ui表示节点i电压；
电力系统储能设备额定功率和额定容量上、下限约束：
其中，PESS,max、PESS,min分别表示电力系统储能设备额定功率上、下限；EESS,max、EESS,min分别表示电力系统储能设备的额定容量的上、下限；PESS,n表示电力系统储能设备功率，EESS,n表示电力系统储能设备的额定容量。
7.根据权利要求1所述的交直流混合配电网规划方法，其特征在于，所述求解所述双层模型的算法为改进粒子群算法，采用所述改进粒子群算法求解的过程中，将所述灵活性资源的位置和容量作为所述改进例子群算法的粒子的位置和速度。
8.根据权利要求7所述的交直流混合配电网规划方法，其特征在于，所述求解所述双层模型的步骤包括：
初始化灵活性资源的位置和容量；
以灵活性资源出力为决策变量，在运行层的目标函数的约束条件下，求解运行层的目标函数，得到系统年运行成本；
将系统年运行成本返回到规划层，以灵活性资源的位置和容量为决策变量，在规划层的目标函数的约束条件下，求解规划层的目标函数，得到更新的灵活性资源的位置和容量；
重复上述的求解运行层的目标函数的步骤和求解规划层的目标函数的步骤，直到满足收敛条件，输出最优的灵活性资源的位置和容量。
9.一种计算机可读存储介质，其特征在于：所述计算机可读存储介质上存储有计算机程序指令；所述计算机程序指令被处理器执行时实现如权利要求1～8中任一项所述的交直流混合配电网规划方法。
10.一种交直流混合配电网规划系统，其特征在于，包括：如权利要求9所述的计算机可读存储介质。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210825774.2A CN115345350A (zh) | 2022-07-14 | 2022-07-14 | 一种交直流混合配电网规划方法、介质及系统 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210825774.2A CN115345350A (zh) | 2022-07-14 | 2022-07-14 | 一种交直流混合配电网规划方法、介质及系统 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN115345350A true CN115345350A (zh) | 2022-11-15 |
Family
ID=83948575
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202210825774.2A Pending CN115345350A (zh) | 2022-07-14 | 2022-07-14 | 一种交直流混合配电网规划方法、介质及系统 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN115345350A (zh) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN118014164A (zh) * | 2024-04-08 | 2024-05-10 | 国网江西省电力有限公司经济技术研究院 | 一种考虑灵活性需求的储能容量配置双层优化方法及系统 |
-
2022
- 2022-07-14 CN CN202210825774.2A patent/CN115345350A/zh active Pending
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN118014164A (zh) * | 2024-04-08 | 2024-05-10 | 国网江西省电力有限公司经济技术研究院 | 一种考虑灵活性需求的储能容量配置双层优化方法及系统 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN103490410B (zh) | 一种基于多目标优化的微电网规划和容量配置方法 | |
| CN110119886B (zh) | 一种主动配网动态规划方法 | |
| CN110429649B (zh) | 考虑灵活性的高渗透率可再生能源集群划分方法 | |
| CN111293691A (zh) | 基于模型预测控制的微电网多时间尺度优化调度方法 | |
| CN107316125A (zh) | 一种基于经济运行域的主动配电网经济运行评价方法 | |
| CN104779611A (zh) | 基于集中式和分布式双层优化策略的微电网经济调度方法 | |
| CN110676849B (zh) | 一种孤岛微电网群能量调度模型的构建方法 | |
| CN113241803A (zh) | 一种基于新能源消纳的储能调度方法及计算机介质 | |
| CN108039741B (zh) | 计及微源余电上网的交直流混合微网优化运行方法 | |
| CN110867907B (zh) | 一种基于多类型发电资源同质化的电力系统调度方法 | |
| CN112383086A (zh) | 一种孤岛微电网日前能量-备用联合优化调度方法 | |
| CN115036914A (zh) | 考虑灵活性与新能源消纳的电网储能双层优化方法及系统 | |
| CN109149658B (zh) | 基于一致性理论的独立微电网分布式动态经济调度方法 | |
| CN117436773B (zh) | 一种含可中断负荷的独立微网资源规划方法和系统 | |
| CN113298407B (zh) | 一种工业园区电-气综合能源系统优化调度模型建立方法 | |
| CN115345350A (zh) | 一种交直流混合配电网规划方法、介质及系统 | |
| CN113361078B (zh) | 一种考虑光伏发电不确定性的配电网云边协同调控方法 | |
| CN113158547B (zh) | 计及经济性和可靠性的区域综合能源系统优化配置方法 | |
| CN113364043A (zh) | 一种基于条件风险价值的微电网群优化方法 | |
| CN114759616B (zh) | 一种考虑电力电子器件特性的微电网鲁棒优化调度方法 | |
| CN115693797A (zh) | 考虑v2g与需求响应的配电网调度方法、介质及系统 | |
| CN109787217A (zh) | 基于风电多状态模型及机会成本修正的备用出清方法 | |
| Fan et al. | Optimal allocation strategy of energy storage to improve the voltage quality of distribution network | |
| CN113902280B (zh) | 一种多微网联合经济调度方法 | |
| Wu | Economic benefit evaluation method for the micro-grid renewable energy system operation |
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 |