CN115099567A - 一种电气互联综合能源系统优化调度方法及终端 - Google Patents
一种电气互联综合能源系统优化调度方法及终端 Download PDFInfo
- Publication number
- CN115099567A CN115099567A CN202210554973.4A CN202210554973A CN115099567A CN 115099567 A CN115099567 A CN 115099567A CN 202210554973 A CN202210554973 A CN 202210554973A CN 115099567 A CN115099567 A CN 115099567A
- Authority
- CN
- China
- Prior art keywords
- constraint
- representing
- energy system
- establishing
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 35
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 claims abstract description 94
- 239000003345 natural gas Substances 0.000 claims abstract description 47
- 238000009826 distribution Methods 0.000 claims abstract description 30
- 238000005457 optimization Methods 0.000 claims abstract description 22
- 238000006243 chemical reaction Methods 0.000 claims abstract description 10
- 239000007789 gas Substances 0.000 claims description 57
- 230000006870 function Effects 0.000 claims description 20
- 230000005611 electricity Effects 0.000 claims description 17
- CURLTUGMZLYLDI-UHFFFAOYSA-N Carbon dioxide Chemical compound O=C=O CURLTUGMZLYLDI-UHFFFAOYSA-N 0.000 claims description 12
- 238000009434 installation Methods 0.000 claims description 12
- 238000012423 maintenance Methods 0.000 claims description 12
- 238000010248 power generation Methods 0.000 claims description 11
- 238000004590 computer program Methods 0.000 claims description 8
- 230000005540 biological transmission Effects 0.000 claims description 7
- 229910002092 carbon dioxide Inorganic materials 0.000 claims description 6
- 239000001569 carbon dioxide Substances 0.000 claims description 6
- 238000004519 manufacturing process Methods 0.000 claims description 5
- 238000005265 energy consumption Methods 0.000 claims description 4
- 238000005485 electric heating Methods 0.000 claims description 3
- 230000009467 reduction Effects 0.000 claims description 3
- 238000005067 remediation Methods 0.000 claims description 3
- 230000009466 transformation Effects 0.000 claims description 3
- 230000020169 heat generation Effects 0.000 claims description 2
- 238000004364 calculation method Methods 0.000 abstract description 8
- 230000008878 coupling Effects 0.000 abstract description 8
- 238000010168 coupling process Methods 0.000 abstract description 8
- 238000005859 coupling reaction Methods 0.000 abstract description 8
- 238000004088 simulation Methods 0.000 description 5
- 230000033228 biological regulation Effects 0.000 description 4
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 3
- 229910052799 carbon Inorganic materials 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 238000012545 processing Methods 0.000 description 2
- 238000004422 calculation algorithm Methods 0.000 description 1
- 230000000295 complement effect Effects 0.000 description 1
- 239000002131 composite material Substances 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000035699 permeability Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 230000002040 relaxant effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000004044 response Effects 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/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/06312—Adjustment or analysis of established resource schedule, e.g. resource or task levelling, or dynamic rescheduling
-
- 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
Landscapes
- Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- Engineering & Computer Science (AREA)
- Economics (AREA)
- Strategic Management (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Entrepreneurship & Innovation (AREA)
- Marketing (AREA)
- General Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- Tourism & Hospitality (AREA)
- Quality & Reliability (AREA)
- Game Theory and Decision Science (AREA)
- Operations Research (AREA)
- Development Economics (AREA)
- Health & Medical Sciences (AREA)
- Educational Administration (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为本发明实施例的电气互联综合能源系统扑图；
图4为本发明实施例的能源集线器内部结构示意图；
图5为本发明实施例二的有功出力情况；
图6为本发明实施例二的能源集线器的电能消费情况；
图7为本发明实施例二的变压器低压侧电压水平；
图8为本发明实施例二的最大负荷时刻和最小负荷时刻的系统电压水平；
标号说明：
1、一种电气互联综合能源系统优化调度终端；2、存储器；3、处理器。
具体实施方式
为详细说明本发明的技术内容、所实现目的及效果，以下结合实施方式并配合附图予以说明。
请参照图1，本发明实施例提供了一种电气互联综合能源系统优化调度方法，包括步骤：
以综合运行费用最小为目标建立电气互联综合能源系统优化调度模型的目标函数；
基于二阶锥和二次约束标准型转换建立电气互联综合能源系统优化调度模型的约束条件，所述约束条件包括能源集线器、配电网和天然气网的约束；
根据所述目标函数和约束条件建立电气互联综合能源系统优化调度模型，并采用求解器对所述调度模型进行求解，以得到综合运行费用最小的综合能源系统调度方案。
从上述描述可知，本发明的有益效果在于：以综合运行费用最小为目标建立电气互联综合能源系统优化调度模型的目标函数，基于二阶锥和二次约束标准型转换建立电气互联综合能源系统优化调度模型的约束条件，约束条件包括能源集线器、配电网和天然气网的约束，从而将电气互联综合能源系统优化调度问题转换为混合整数二阶锥-二次优化问题；根据目标函数和约束条件建立电气互联综合能源系统优化调度模型，并采用求解器对所述调度模型进行求解，以得到综合运行费用最小的综合能源系统调度方案，提升了计算效率，适合工程推广；充分考虑配电网、天然气网的深度耦合，更能适应多能深度耦合的综合能源系统。
进一步地，所述以综合运行费用最小为目标建立电气互联综合能源系统优化调度模型的目标函数包括：
以综合运行费用最小为目标建立目标函数，所述综合运行费用包括用能成本CEN、网损成本CL、弃风弃光惩罚CP和二氧化碳治理成本CT：
F1＝min(CEN+CL+CP+CT)。
由上述描述可知，使用综合运行费用最小为目标建立目标函数，便于后续基于目标函数建立电气互联综合能源系统优化调度模型。
进一步地，计算所述用能成本CEN的公式为：
其中，T表示时刻集合，ΩWTG和ΩPVG分别表示风机和光伏发电的安装集合，CE表示用电成本，CG表示用气成本，CDG表示分布式电源运维成本，fe表示度电费用，fg表示每立方米天然气费用，fWTG表示风机度电运维成本，fPVG表示光伏度电运维成本，Pgen表示上级电网电能输入，Vgas表示天然气输入，PWTG和PPVG分别表示风机和光伏发电的电力输入；
计算所述网损成本CL的公式为：
其中floss表示线损费用，Ploss表示系统总线损；
计算所述弃风弃光惩罚CP的公式为：
其中，fp.WTG和fp.PVG分别表示单位弃风和弃光的惩罚成本，PWTG.cut和PWTG.cut分别表示风机和光伏发电的有功削减量；
计算所述二氧化碳治理成本CT的公式为：
其中
和
分别表示电力和天然气的CO2排放因子。
由上述描述可知，通过计算用能成本、网损成本、弃风弃光惩罚和二氧化碳治理成本，便于计算综合运行费用的最小值并建立目标函数。
进一步地，所述建立电气互联综合能源系统优化调度模型的能源集线器约束包括：
建立能源集线器的能量平衡约束：
其中，ΩEH表示能源集线器的节点集合；t表示第t个调度时段
分别为能源集线器输入的电能和天然气；
分别为能源集线器节点的电负荷和热负荷；
分别为热电联产的产电和产热效率，ηGB为燃气锅炉的产热效率、ηT为变压器效率；
建立热电联产和燃气锅炉的电、热能量流边界约束：
其中
分别表示热电联产的电热功率，
表示燃气锅炉的热功率，
分别表示热电联产电功率的下限和上限，
分别表示热电联产热功率的下限和上限，
表示燃气锅炉热功率的下限和上限。
由上述描述可知，通过建立能源集线器的能量平衡约束以及热电联产和燃气锅炉的电、热能量流边界约束，能够得到能源集线器约束，以便于后续建立电气互联综合能源系统优化调度模型。
进一步地，基于二阶锥转换建立电气互联综合能源系统优化调度模型的配电网约束包括：
建立潮流约束：
其中，δ(j)表示以节点j为首节点的支路，π(j)表示以节点j为末节点的支路，Pj,t、Qj,t分别表示节点j的注入的有功和无功功率，Pij,t、Qij,t分别表示支路ij的有功和无功潮流，Rij、Xij分别表示支路ij的电阻和电抗，Iij,t表示支路电流；
对
进行松弛处理并改写为二阶锥形式：
由上述描述可知，通过将潮流约束进行松弛处理并改写为二阶锥形式，便于后续使用商业求解器进行模型求解。
进一步地，所述建立电气互联综合能源系统优化调度模型的配电网约束还包括：
建立分布式风机约束：
其中，ΩWTG表示风机安装节点，
分别表示风机投入数量和最大投入数量，
分别表示单台风机的有功无功和最大无功出力，sWTG表示风机的额定视在功率；
建立分布式光伏发电约束：
其中，ΩPVG表示风机安装节点，
分别表示风机投入数量和最大投入数量，
分别表示单台风机的有功和最大无功出力，sPVG表示光伏发电的额定视在功率；
建立变压器分接头约束：
其中Vseq表示变压器可调的电压序列，
表示0-1辅助变量。
由上述描述可知，在潮流约束进行松弛处理并改写为二阶锥形式之后，建立分布式风机、分布式光伏发电和变压器分接头约束，充分利用变压器、分布式电源等主动调节手段，深度挖掘电网潜力，满足负荷和分布式电源接入需求。
进一步地，基于二次约束标准型转换建立电气互联综合能源系统优化调度模型的天然气网约束包括：
将管道潮流约束转换为二次约束标准型：
其中fij.t表示管道ij内的气潮流，∏i,t、∏j,t分别表示气节点i和j气压的平方，φp表示天然气流传输系数；α1、α2均为二进制辅助变量；
建立气源约束：
Smin≤Sh≤Smax；
式中，Sh表示气源h的出气量，Smin、Smax分别表示气源出气量的下限和上限；
建立节点气压约束：
式中，∏min、∏max表示节点气压平方的下限和上限；
建立加压器约束：
式中，fc表示加压器c流过的气流，Γc表示加压器c的升压比例，
表示加压器c的传输容量上限；
建立节点气流平衡约束：
其中，σ(j)表示以j为首端节点的管道集合，即天然气支路jk、μ(j)分别表示以j为末端节点的管道集合，即天然气支路ij，fij表示管道ij的气流量，
Sj分别表示热负荷等效气流量和气源节点出气量，
表示能源集线器节点消耗的天然气负荷。
由上述描述可知，将管道潮流约束转换为二次约束标准型，以此方式，将电气互联综合能源系统优化调度问题转换为混合整数二阶锥-二次优化问题，便于后续采用现有的商业求解器进行求解，极大提升了计算效率，适合工程推广。
请参照图2，本发明另一实施例提供了一种电气互联综合能源系统优化调度终端，包括存储器、处理器及存储在所述存储器上并可在所述处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现上述一种电气互联综合能源系统优化调度方法。
本发明上述的一种电气互联综合能源系统优化调度方法及终端，适用于建立电气互联综合能源系统优化调度模型，实现综合能源系统优化调度，以下通过具体的实施方式进行说明：
实施例一
请参照图1，一种电气互联综合能源系统优化调度方法，包括步骤：
S1、以综合运行费用最小为目标建立电气互联综合能源系统优化调度模型的目标函数。
以综合运行费用最小为目标，综合运行费用包括用能成本CEN，网损成本CL，弃风弃光惩罚CP、二氧化碳治理成本CT：
F1＝min(CEN+CL+CP+CT)；
其中，T为时刻集合，ΩWTG和ΩPVG分别表示WTG(wind turbine generator，风机)和PVG(photovoltaic generator，光伏发电)的安装集合，CE表示用电成本，CG表示用气成本，CDG表示分布式电源运维成本，fe表示度电费用，fg表示每立方米天然气费用，fWTG表示WTG度电运维成本，fPVG表示PVG度电运维成本，Pgen表示上级电网电能输入，Vgas表示天然气输入，PWTG和PPVG分别表示WTG和PVG电力输入。
其中floss表示线损费用，Ploss表示系统总线损。
其中，fp.WTG和fp.PVG分别表示单位弃风弃光惩罚成本，PWTG.cut和PWTG.cut分别表示WTG和PVG的有功削减量。
其中
和
分别表示电力和天然气的CO2排放因子。
S2、基于二阶锥和二次约束标准型转换建立电气互联综合能源系统优化调度模型的约束条件，所述约束条件包括能源集线器、配电网和天然气网的约束。
S21、建立能源集线器约束。
能源集线器的能量平衡约束：
其中，ΩEH表示EH(energy hub，能源集线器)节点集合；t表示第t个调度时段
分别为EH输入的电能和天然气；
分别为EH节点的电负荷和热负荷；
分别为CHP的产电和产热效率，ηGB为GB的产热效率、ηT为变压器效率。
热电联产(CHP)和燃气锅炉(GB)的电、热能量流边界约束：
其中
分别表示CHP的电热功率，
表示GB的热功率。
表示CHP电功率的下限和上限；
表示CHP热功率的下限和上限；
表示GB热功率的下限和上限。
S22、建立配电网约束；
采用基于distflow支路潮流模型的最优潮流(branch flow model-optimalpower flow,BFM-OPF)模型，潮流约束具体为：
其中，δ(j)表示以节点j为首节点的支路，π(j)表示以节点j为末节点的支路，Pj,t、Qj,t分别表示节点j的注入的有功和无功功率，Pij,t、Qij,t分别表示支路ij的有功和无功潮流，Rij、Xij分别表示支路电阻和电抗，Iij,t表示支路电流。
配电网支路电流、节点电压约束如下：
其中Iij,min、Iij,max分别表示支路电流的下限和上限；Vj,min、Vj,max分别表示节点电压的下限和上限；
根据基尔霍夫定律，对于非EH的负荷节点，则有：
对于EH节点，由于EH不输出无功功率，其有功负荷改为EH的输入功率，即无功平衡为：
其中，Ω、ΩEH分别表示非EH节点集合和EH节点集合，
分别表示上级电网输入的有功、无功功率，
分别表示风机输出的有功和无功功率，
分别表示光伏输出的有功和无功功率。
配电网上级电源出力约束如下：
其中
分别表示电源有功出力的下限和上限；
分别表示电源无功出力的下限和上限；Ωgen表示电源节点集合。
OLTC(on-load tap charger，有载调压变压器)约束如下：
其中
表示OLTC电压调节下限，ΔVOLTC表示电压调节跨度，TRt表示分接头位置，TRmin、TRmax分别表示分接头位置的上下限。
令
则将潮流约束转换为：
松弛处理
得到，
将其改写为二阶锥形式：
因此，配电网支路电流、节点电压约束改写如下：
其中
分别表示支路电流平方的下限和上限；
分别表示节点电压平方的下限和上限；
其中的分布式WTG约束如下；
其中，ΩWTG表示风机安装节点，
分别表示风机投入数量和最大投入数量，
分别表示单台风机的有功无功和最大无功出力，sWTG表示WTG的额定视在功率；
同理可得分布式PVG约束如下；
其中，ΩPVG表示风机安装节点，
分别表示风机投入数量和最大投入数量，
分别表示单台风机的有功和最大无功出力，sPVG表示PVG的额定视在功率；
由于电压、电流经过变换导致OLTC约束变为二次约束，为解决OLTC约束的非线性问题，引入0-1辅助变量
其中Vseq表示OLTC可调的电压序列，其表达式为：
考虑到OLTC只可选择一个位置，因此可得：
S23、建立天然气网约束。
采用Weymouth稳态天然气管道潮流模型，其管道潮流约束为：
其中fij.t表示管道ij内的气潮流，πi.t、πj.t分别表示气节点i和j的气压，φp表示天然气流传输系数。
气源约束为：
Smin≤Sh≤Smax；
式中，Sh为气源h的出气量，Smin、Smax分别为气源出气量的下限和上限。
节点气压约束为：
加压器约束为：
式中，fc、πi,t和πj,t分别为加压器c流过的气流、进气口和出气口端的气压；Γc为加压器c的升压比例；
为加压器c的传输容量上限。
节点气流平衡约束为：
式中，σ(j)表示以j为首端节点的管道集合、μ(j)分别表示以j为末端节点的管道集合；fij为管道ij的气流量；
Sj分别为热负荷等效气流量和气源节点出气量；
表示EH节点消耗的天然气负荷；
令
则天然气潮流约束，并使用big_M法对管道潮流约束进行处理，转换为二次约束标准型；
其中fij.t表示管道ij内的气潮流，∏i,t、∏j,t分别表示气节点i和j气压的平方，φp表示天然气流传输系数；α1、α2均为二进制辅助变量；
节点气压约束为：
式中，∏min、∏max表示节点气压平方的下限和上限。
加压器约束为：
式中，fc为加压器c流过的气流，Γc为加压器c的升压比例，
为加压器c的传输容量上限。
S3、根据所述目标函数和约束条件建立电气互联综合能源系统优化调度模型，并采用商业求解器Gurobi对所述调度模型进行求解，以得到综合运行费用最小的综合能源系统调度方案。
具体的，求解得到的综合能源系统调度方案，其重点是配电网的优化调度，具体包括分布式电源的出力情况、OLTC分接头位置等参数。
与现有技术相比，本实施例具有以下有益效果：
1、充分利用OLTC、分布式电源等主动调节手段，深度挖掘电网潜力，满足负荷和分布式电源接入需求。
2、充分考虑配电网、天然气网的深度耦合，其优化调度计算方法更能适应多能深度耦合的综合能源系统。
3、将电-气互联综合能源系统优化调度问题转换为混合整数二阶锥-二次优化问题，并采用现有的商业求解器Gurobi进行求解，极大提升了计算效率，适合工程推广。
实施例二
本实施例采用如图3所示的电-气互联综合能源系统本发明所述计算方法进行仿真验证，能源集线器内部结构如图4所示。
在本实施例的仿真算例中，配电网部分是IEEE 33节点系统，天然气网是比利时20节点系统，通过4个能源集线器实现能量交换。能源集线器参数如表1所示。
表1能源集线器参数
本实施例中参数说明如下，配电网中，设节点2-18为居民负荷节点；节点26-33为商业负荷节点；节点19-25为工业负荷节点；WTG安装节点为节点7、17，安装数量为5、10组；PVG安装节点为节点12、26，安装数量分别为5、10组，单台WTG和PVG的额定容量均为0.2MVA，分布式电源总安装容量为6MVA，WTG额定功率因数为0.85，PVG额定功率因数为1，电网总负荷设置为基准负荷的1.5倍，即7.63MW+3.82Mvar，分布式电源渗透率约为70％，OLTC共有9个调节挡位，每个挡位的调节电压是0.0125p.u.，节点电压可调节范围是0.95～1.05p.u.；配电网电压允许偏差为±7％p.u.。
天然气网中，气源节点最大出气量1100m3/h，节点气压上下限分别为460kPa和380kPa，管道最大允许气流量为800m3/h，加压器安装于天然气网支路4-7和支路17-18，加压比分别为1.15和1.1。
能源集线器参数如下，CHP产热效率为0.57，产电效率为0.48；单位体积天然气热值为9.6kWh/m3。
系统网损成本为0.5元/kWh，弃风弃光惩罚费用为2元/kWh，WTG和PVG的运维成本均为0.2元/kWh，购电价格为0.42元/kWh，购气价格为5.5元/m3；
单位电能CO2排放为0.8647kg/kWh，单位体积天然气CO2排放为1.9kg/m3，碳交易价格为267元/t。
本实施例的年运行费用情况如表2所示，本实施例仿真所得的电-气互联综合能源系统总运行成本为3361.4万元，用能成本中，向上级电网购电成本为433.8万元，DG运维成本为840.6万元，这是由于虽然上级电网购电成本为0.42元/kWh，低于0.5元/kWh的DG运维成本，但上级电网产生的电能需要额外支付0.0231元/kWh的CO2治理成本，而DG则不需要额外支付CO2的治理费用，因此考虑CO2的治理成本后，DG使用成本更低，不会出现弃风弃光现象，弃风弃光成本为0。
表2本文方法的年运行费用
配电网有功出力情况请参照图5，CHP工作在冬季能够承担较多的电负荷供应，这是由于CHP工作在“以热定电”模式，冬季热负荷需求大，相应输出更多电能。高比例电-气互联系统采用风、光互补模式能够在较大程度上减少高峰电力需求，四季电负荷高峰均得到一定程度的削弱，但是由于分布式电源“夜晚无光、极热无风”的特点，且夏季无热负荷需求，CHP无法输出电能，导致夏季晚高峰问题仍较为突出，应进一步采取需求响应的方式，缓解夏季晚高峰供电压力。
图6为能源集线器的电能消耗，其中负值为能源集线器向电网反送电能，从图6可以看出，能源集线器向电网反送电能多发于春、秋、冬季的夜间，夜间电力负荷较低，且WTG出力较大，配电网电能消纳压力增大，若不采取对应调压措施，配电网系统极易发生过电压问题。
OLTC分接头位置请参照图7，最大、最小负荷时刻配电网电压情况请参照图8，其中最大负荷时刻为夏季12：00，最小负荷时刻为春季5：00。
请参照图7，考虑到目标函数中的线损成本，为降低线损成本，OLTC分接头会倾向于调高电网系统电压，在夏季，由于负荷较重，且综合能源系统无电力反送，电能消纳问题较小，OLTC低压侧电压多保持在较高水平，从图8可知，最大负荷时刻，OLTC低压侧电压为1.05p.u.。但在春、秋、冬季节，DG和CHP会向电网反送大量电能，电力消纳压力大，系统电压显著提升，过电压风险严重，从图7可以看出，OLTC倾向于调低系统电压，以满足电压约束，从图8可知，在最小负荷时刻，OLTC将低压侧电压调低至1.0125p.u.，避免系统出现过电压问题。
实施例三
请参照图2，一种电气互联综合能源系统优化调度终端1，包括存储器2、处理器3以及存储在所述存储器2上并可在处理器3上运行的计算机程序，所述处理器3执行所述计算机程序时实现实施例一或二的一种电气互联综合能源系统优化调度方法的各个步骤。
综上所述，本发明提供的一种电气互联综合能源系统优化调度方法及终端，基于电-气互联综合能源系统最优潮流算法，进一步引入配电网各主动管理元素，建立电-气互联综合能源系统优化调度模型，通过能源集线器耦合配电网和天然气气网，以经济性、碳排放为目标，优化综合能源系统各主动控制元素。引入辅助变量线性化OLTC电压调节约束，采用二阶锥松弛技术处理distflow配电网支路潮流，采用big_M法处理Weymouth天然气管道流量模型。通过数学转换，所提模型被转换为混合整数二阶锥-二次优化问题，采用Guribi求解器求解，以得到综合运行费用最小的综合能源系统调度方案。通过建立IEEE 33节点配电网系统和比利时20节点天然气网系统相耦合的电-气互联综合能源系统作为本发明实施例并进行仿真，仿真结果表明，本发明所提方法兼顾经济性和碳排放水平，具有较高的计算精度和较快的运算速度，可适用于电-气互联综合能源系统优化调度。
以上所述仅为本发明的实施例，并非因此限制本发明的专利范围，凡是利用本发明说明书及附图内容所作的等同变换，或直接或间接运用在相关的技术领域，均同理包括在本发明的专利保护范围内。
Claims (8)
1.一种电气互联综合能源系统优化调度方法，其特征在于，包括步骤：
以综合运行费用最小为目标建立电气互联综合能源系统优化调度模型的目标函数；
基于二阶锥和二次约束标准型转换建立电气互联综合能源系统优化调度模型的约束条件，所述约束条件包括能源集线器、配电网和天然气网的约束；
根据所述目标函数和约束条件建立电气互联综合能源系统优化调度模型，并采用求解器对所述调度模型进行求解，以得到综合运行费用最小的综合能源系统调度方案。
2.根据权利要求1所述的一种电气互联综合能源系统优化调度方法，其特征在于，所述以综合运行费用最小为目标建立电气互联综合能源系统优化调度模型的目标函数包括：
以综合运行费用最小为目标建立目标函数，所述综合运行费用包括用能成本CEN、网损成本CL、弃风弃光惩罚CP和二氧化碳治理成本CT：
F1＝min(CEN+CL+CP+CT)。
3.根据权利要求2所述的一种电气互联综合能源系统优化调度方法，其特征在于，计算所述用能成本CEN的公式为：
其中，T表示时刻集合，ΩWTG和ΩPVG分别表示风机和光伏发电的安装集合，CE表示用电成本，CG表示用气成本，CDG表示分布式电源运维成本，fe表示度电费用，fg表示每立方米天然气费用，fWTG表示风机度电运维成本，fPVG表示光伏度电运维成本，Pgen表示上级电网电能输入，Vgas表示天然气输入，PWTG和PPVG分别表示风机和光伏发电的电力输入；
计算所述网损成本CL的公式为：
其中floss表示线损费用，Ploss表示系统总线损；
计算所述弃风弃光惩罚CP的公式为：
其中，fp.WTG和fp.PVG分别表示单位弃风和弃光的惩罚成本，PWTG.cut和PWTG.cut分别表示风机和光伏发电的有功削减量；
计算所述二氧化碳治理成本CT的公式为：
其中
和
分别表示电力和天然气的CO2排放因子。
4.根据权利要求1所述的一种电气互联综合能源系统优化调度方法，其特征在于，所述建立电气互联综合能源系统优化调度模型的能源集线器约束包括：
建立能源集线器的能量平衡约束：
其中，ΩEH表示能源集线器的节点集合；t表示第t个调度时段
分别为能源集线器输入的电能和天然气；
分别为能源集线器节点的电负荷和热负荷；
分别为热电联产的产电和产热效率，ηGB为燃气锅炉的产热效率、ηT为变压器效率；
建立热电联产和燃气锅炉的电、热能量流边界约束：
其中
分别表示热电联产的电热功率，
表示燃气锅炉的热功率，
分别表示热电联产电功率的下限和上限，
分别表示热电联产热功率的下限和上限，
表示燃气锅炉热功率的下限和上限。
5.根据权利要求1所述的一种电气互联综合能源系统优化调度方法，其特征在于，基于二阶锥转换建立电气互联综合能源系统优化调度模型的配电网约束包括：
建立潮流约束：
其中，δ(j)表示以节点j为首节点的支路，π(j)表示以节点j为末节点的支路，Pj,t、Qj,t分别表示节点j的注入的有功和无功功率，Pij,t、Qij,t分别表示支路ij的有功和无功潮流，Rij、Xij分别表示支路ij的电阻和电抗，Iij,t表示支路电流；
对
进行松弛处理并改写为二阶锥形式：
6.根据权利要求5所述的一种电气互联综合能源系统优化调度方法，其特征在于，所述建立电气互联综合能源系统优化调度模型的配电网约束还包括：
建立分布式风机约束：
其中，ΩWTG表示风机安装节点，
分别表示风机投入数量和最大投入数量，
分别表示单台风机的有功无功和最大无功出力，sWTG表示风机的额定视在功率；
建立分布式光伏发电约束：
其中，ΩPVG表示风机安装节点，
分别表示风机投入数量和最大投入数量，
分别表示单台风机的有功和最大无功出力，sPVG表示光伏发电的额定视在功率；
建立变压器分接头约束：
其中Vseq表示变压器可调的电压序列，
表示0-1辅助变量。
7.根据权利要求1所述的一种电气互联综合能源系统优化调度方法，其特征在于，基于二次约束标准型转换建立电气互联综合能源系统优化调度模型的天然气网约束包括：
将管道潮流约束转换为二次约束标准型：
其中fij.t表示管道ij内的气潮流，∏i,t、∏j,t分别表示气节点i和j气压的平方，φp表示天然气流传输系数；α1、α2均为二进制辅助变量；
建立气源约束：
Smin≤Sh≤Smax；
式中，Sh表示气源h的出气量，Smin、Smax分别表示气源出气量的下限和上限；
建立节点气压约束：
式中，∏min、∏max表示节点气压平方的下限和上限；
建立加压器约束：
式中，fc表示加压器c流过的气流，Γc表示加压器c的升压比例，Fc max表示加压器c的传输容量上限；
建立节点气流平衡约束：
其中，σ(j)表示以j为首端节点的管道集合，即天然气支路jk、μ(j)分别表示以j为末端节点的管道集合，即天然气支路ij，fij表示管道ij的气流量，fj L、Sj分别表示热负荷等效气流量和气源节点出气量，
表示EH节点消耗的天然气负荷。
8.一种电气互联综合能源系统优化调度终端，其特征在于，包括存储器、处理器及存储在所述存储器上并可在所述处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现上述权利要求1-7任一项所述的一种电气互联综合能源系统优化调度方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210554973.4A CN115099567A (zh) | 2022-05-19 | 2022-05-19 | 一种电气互联综合能源系统优化调度方法及终端 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210554973.4A CN115099567A (zh) | 2022-05-19 | 2022-05-19 | 一种电气互联综合能源系统优化调度方法及终端 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN115099567A true CN115099567A (zh) | 2022-09-23 |
Family
ID=83289283
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202210554973.4A Pending CN115099567A (zh) | 2022-05-19 | 2022-05-19 | 一种电气互联综合能源系统优化调度方法及终端 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN115099567A (zh) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115345391A (zh) * | 2022-10-20 | 2022-11-15 | 广东电网有限责任公司 | 一种电-气能源系统的灾后恢复方法、装置及存储介质 |
| CN116384536A (zh) * | 2023-01-06 | 2023-07-04 | 国网江苏省电力有限公司扬州供电分公司 | 一种中大型能源用户综合能源协同规划方法及装置 |
-
2022
- 2022-05-19 CN CN202210554973.4A patent/CN115099567A/zh active Pending
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115345391A (zh) * | 2022-10-20 | 2022-11-15 | 广东电网有限责任公司 | 一种电-气能源系统的灾后恢复方法、装置及存储介质 |
| CN116384536A (zh) * | 2023-01-06 | 2023-07-04 | 国网江苏省电力有限公司扬州供电分公司 | 一种中大型能源用户综合能源协同规划方法及装置 |
| CN116384536B (zh) * | 2023-01-06 | 2024-05-07 | 国网江苏省电力有限公司扬州供电分公司 | 一种中大型能源用户综合能源协同规划方法 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108846507B (zh) | 基于混合整数二阶锥规划的电-气耦合系统日前经济调度方法 | |
| CN110245878B (zh) | 智能楼宇群体的分布式综合能源需求响应协同优化方法 | |
| CN106372742A (zh) | 考虑电转气多源储能型微网日前最优经济调度方法 | |
| CN110866627B (zh) | 一种考虑阶梯气价的多区域电-气耦合综合能源系统优化调度方法 | |
| CN115099567A (zh) | 一种电气互联综合能源系统优化调度方法及终端 | |
| CN107808216B (zh) | 电-气-热互联系统弃风弃光和电气热负荷削减综合最小优化模型构建方法 | |
| CN110955954B (zh) | 一种分层解耦的电气热综合能源系统最优负荷削减量方法 | |
| CN107358345B (zh) | 计及需求侧管理的分布式冷热电联供系统优化运行方法 | |
| CN112035984B (zh) | 一种电-气-储区域综合能源系统协同规划方法 | |
| CN110318877B (zh) | 基于双模态解耦的热电联产系统的风电消纳方法 | |
| CN111815068A (zh) | 二阶段约束边界紧缩算法求解城市综合能源网络优化方法 | |
| CN112836882B (zh) | 考虑设备负载率变化的区域综合能源系统运行优化方法 | |
| CN104617585A (zh) | 无功补偿配置方法 | |
| CN107404118A (zh) | 基于随机响应面法的电气互联系统概率最优潮流计算方法 | |
| CN114066204A (zh) | 一种综合能源系统的一体化优化规划与运行方法及装置 | |
| CN115186902A (zh) | 温室综合能源系统的调控方法、装置、终端及存储介质 | |
| CN108594658A (zh) | 一种电-气耦合系统概率最大负荷裕度多目标优化模型及其求解方法 | |
| CN113806952B (zh) | 一种考虑源-荷-储的冷热电综合能源系统及其优化运行方法 | |
| CN113131513B (zh) | 考虑碳排放的电、热、气转换优化运行配置方法和存储介质 | |
| CN112036652A (zh) | 基于机会约束规划的光伏-储能综合能源系统规划方法 | |
| CN104899659A (zh) | 一种智能园区多级能耗传递模型 | |
| CN114221346A (zh) | 一种综合能源系统潮流计算方法 | |
| Liu et al. | Low carbon economic dispatch for multi-energy distribution network with compressed air energy storage system as energy hub | |
| CN112701721A (zh) | 一种综合能源系统的协调规划方法 | |
| CN111724026A (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 |