CN113570117B - 一种基于二阶锥松弛转化方法的电-气综合能源微网优化调度方法 - Google Patents
一种基于二阶锥松弛转化方法的电-气综合能源微网优化调度方法 Download PDFInfo
- Publication number
- CN113570117B CN113570117B CN202110752294.3A CN202110752294A CN113570117B CN 113570117 B CN113570117 B CN 113570117B CN 202110752294 A CN202110752294 A CN 202110752294A CN 113570117 B CN113570117 B CN 113570117B
- Authority
- CN
- China
- Prior art keywords
- micro
- grid
- cost
- power
- order cone
- 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 63
- 238000006243 chemical reaction Methods 0.000 title claims abstract description 49
- 238000004146 energy storage Methods 0.000 claims abstract description 10
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 claims description 44
- CURLTUGMZLYLDI-UHFFFAOYSA-N Carbon dioxide Chemical compound O=C=O CURLTUGMZLYLDI-UHFFFAOYSA-N 0.000 claims description 24
- 239000007789 gas Substances 0.000 claims description 23
- 230000005611 electricity Effects 0.000 claims description 21
- 239000003345 natural gas Substances 0.000 claims description 19
- 239000000243 solution Substances 0.000 claims description 18
- 230000007613 environmental effect Effects 0.000 claims description 17
- 238000005457 optimization Methods 0.000 claims description 16
- UFHFLCQGNIYNRP-UHFFFAOYSA-N Hydrogen Chemical compound [H][H] UFHFLCQGNIYNRP-UHFFFAOYSA-N 0.000 claims description 12
- 229910002092 carbon dioxide Inorganic materials 0.000 claims description 12
- 230000008901 benefit Effects 0.000 claims description 11
- 230000009466 transformation Effects 0.000 claims description 11
- 239000002699 waste material Substances 0.000 claims description 11
- 239000001257 hydrogen Substances 0.000 claims description 9
- 229910052739 hydrogen Inorganic materials 0.000 claims description 9
- 230000008569 process Effects 0.000 claims description 9
- 238000002347 injection Methods 0.000 claims description 7
- 239000007924 injection Substances 0.000 claims description 7
- 239000001569 carbon dioxide Substances 0.000 claims description 6
- 238000004458 analytical method Methods 0.000 claims description 5
- 238000010438 heat treatment Methods 0.000 claims description 4
- 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 claims description 3
- 229910052799 carbon Inorganic materials 0.000 claims description 3
- 239000000446 fuel Substances 0.000 claims description 3
- 238000004519 manufacturing process Methods 0.000 abstract description 2
- 238000004422 calculation algorithm Methods 0.000 description 10
- 238000011160 research Methods 0.000 description 7
- 238000010586 diagram Methods 0.000 description 4
- 238000013473 artificial intelligence Methods 0.000 description 3
- 238000010248 power generation Methods 0.000 description 3
- 230000008878 coupling Effects 0.000 description 2
- 238000010168 coupling process Methods 0.000 description 2
- 238000005859 coupling reaction Methods 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 230000003044 adaptive effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 239000003795 chemical substances by application Substances 0.000 description 1
- 238000001816 cooling Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000003912 environmental pollution Methods 0.000 description 1
- 239000002803 fossil fuel Substances 0.000 description 1
- 238000012806 monitoring device Methods 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000011426 transformation method Methods 0.000 description 1
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
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/28—Arrangements for balancing of the load in a network by storage of energy
-
- 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
- H02J3/466—Scheduling the operation of the generators, e.g. connecting or disconnecting generators to meet a given 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
- H02J15/00—Systems for storing electric energy
- H02J15/008—Systems for storing electric energy using hydrogen as energy vector
-
- 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]
-
- 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
- H02J2300/00—Systems for supplying or distributing electric power characterised by decentralized, dispersed, or local generation
- H02J2300/20—The dispersed energy generation being of renewable origin
- H02J2300/28—The renewable source being wind energy
-
- 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
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E40/00—Technologies for an efficient electrical power generation, transmission or distribution
- Y02E40/70—Smart grids as climate change mitigation technology in the energy generation sector
-
- 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)
- Economics (AREA)
- Human Resources & Organizations (AREA)
- Power Engineering (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
技术领域
本发明属于微网优化调度策略领域，具体涉及一种基于二阶锥松弛转化方法的电-气综合能源微网优化调度方法。
背景技术
近些年，随着全球经济的迅速发展，对于能源的需求变得越来越大，这导致化石燃料迅速枯竭，全球环境污染日益加剧。微网作为一个包含分布式能源、储能系统、热电气负荷、保护与监控装置的小型系统，在未来的能源网络和结构中，将起到重要作用，它能充分发挥分布式电源的经济效益和环保效益，具有巨大的社会与经济意义，同时它对于维护系统稳定、保障某些重要负荷安全有着很大的优势，因而对于微网系统的研究日益重要。
目前对于微网的研究多为结构和算法方面。在微网结构方面，主要涉及微网中各个能源的类型、耦合关系等，相关的研究有基于电-冷-热结构，建立冷热电联供的系统模型；考虑沼-风-光之间的耦合关系，建立综合能源微网优化模型等等。此外储能系统因其能灵活调节的优势，在储能优化中起着重要的作用，所以包含储能系统的微网研究也很广泛，此外氢能作为一种清洁能源在微网系统中的应用正越来越广泛，考虑以氢能为主导的综合能源系统为未来综合能源系统的发展提供了方向。关于微网优化算法，已有一定的研究，常见的有以系统总成本最小为目标,采用混合整数规划方法求解的方法；为应对新能源的需要，自适应加权和算法相结合的近似动态规划算法被提出以较好地满足以新能源发电为主电源的独立微网日前调度需求；针对具体微网，粒子群算法和人工蜂群算法被引入以解决微网优化调度问题。此外，针对传统集中式优化调度方法难以全面反映综合能源微网内不同智能体的利益诉求,以及人工智能技术在综合能源调度方面的应用亟待进一步挖掘等问题，人工智能算法正成为一个新的研究方向。
微网系统进行建模时，由于微网潮流约束存在非线性约束边界条件，需进行线性转化。以往基于人工智能的启发式算法因其简单、较容易模拟复杂约束等优势，在求解非线性优化调度问题时得到了普遍应用。但智能算法在求解非线性模型时仍有明显不足，比如求解速度慢、容易陷入局部最优而无法保证全局最优，这些问题都促使研究人员探索研究更有效的求解算法，以满足全局最优、求解速度快等要求。
发明内容
为了克服现有技术的不足，本发明提出一种基于二阶锥松弛转化方法的电-气综合能源微网优化调度方法，对系统中各能源生产、转换及储能设备建模，建立电-气综合能源微网系统模型；运用二阶锥松弛方法对支路潮流非线性约束进行转化；在此基础上，基于所建立的微网模型，研究二阶锥松弛转化方法的应用条件，提出以微网总成本为目标函数的优化调度策略。本发明的优化调度策略可使微网总成本达到最优值，此外提出在目标函数中增加线路损耗的方法可以使微网系统在没有满足二阶锥规划应用条件的情况下正常运行，使总成本最优的同时，提高系统的适应性和灵活性。
为了实现上述目的，本发明的技术方案为：
一种基于二阶锥松弛转化方法的电-气综合能源微网优化调度方法，所述方法包括以下步骤：
S1：建立电-气综合能源微网系统及内部设备模型，设定设备运行相关参数，将包含运行成本和环境成本的综合运行成本设定为目标函数；
电-气综合能源微网系统涉及电力网络、氢储能系统、天然气网络三部分，微网系统运行的优化目标为综合运行成本最低，综合运行成本包括系统运行成本以及环境成本；
系统运行成本表示为：
f1＝f11+f12+f13 (1)
式(1)中f11为微网系统与主网能量交换的成本与收益，表示为：
其中Ce1为微网通过主网购电分时电价，而Ce2为微网向主网售电分时电价；表示t时段微网通过主网购电的有功功率，/>表示t时段微网向主网售电的有功功率；
式(1)中f12为微网向天然气网络购售天然气的成本与收益，表示为：
其中Cg1表示天然气购买的价格，Cg2表示天然气出售价格；表示t时段微网从气网购买的CH4功率，/>表示t时段微网出售的CH4功率；H表示天然气高热值；
式(1)中f13为微网系统运行能量损失惩罚成本，表示为：
其中KLoss为权重系数，表示系统对于能量损失的重视程度，0≤KLoss≤1；和/>分别表示微网系统中电解槽、燃料电池、氢转气装置和微型燃气轮机MT在能量转化时的能量损失；
为简化处理，在此微网系统中各种能量形式的转化均认为不会产生有害气体，因此环境成本由CO2排放量构成，环境成本表达式为：
其中是该区域电网单位供电平均CO2排放系数，/>是MT发电二氧化碳排放系数，/>是CH4的二氧化碳捕捉系数；
结合以上数学表达式，考虑经济性与环保性的综合运行成本目标函数为：
式中f1 max与分别为微网系统最大运行成本和最大环境成本，ω1与ω2为系统运行成本和环境成本各自的权重系数；通过标幺值方法，根据重视程度对系统运行成本和环境成本赋予不同的权重系数，将多目标优化问题转化为单目标问题；
S2：通过二阶锥松弛方法将非线性的潮流约束线性化，从而将非线性规划问题转化为SOCP规划问题进行求解，微电网潮流基本模型表示如下：
式中：p、q为节点注入有功和无功功率，P、Q为支路潮流有功和无功功率，V为节点电压，I为支路电流；
在式(8)-(13)中，支路阻抗rij、xij，节点对地导纳gj、bj均为常数，由于微网系统内部电压等级为400V，因此对地导纳gj、bj对微电网潮流的影响很小，为简便运算，将式(8)中的与/>进行忽略，节点电压V、支路电流I、支路潮流P、Q以及节点注入功率p、q为优化变量，显然它们之间为非线性关系，式(8)-(12)为一般电网潮流所共有的约束方程，式(13)为每个微电网系统所特有的约束条件，由微网自身的电力负荷、微源及其它微网内部结构决定；
通过二阶锥松弛方法将模型中的非线性约束线性化，从而将系统的非线性规划问题转化为二阶锥规划(second-order cone program,SOCP)，对式(10)进行转化处理，首先定义并进行二阶锥松弛转化，得到：
式(14)可以等价转化为2型范数表达的二阶锥松弛的标准形式，再将代入式(8)-(9)，即得到通过二阶锥松弛转化后的最优潮流约束方程：
S3：在微网系统没有满足二阶锥规划应用条件的情况下，提出通过增加线路损耗以解决问题的方法；
原问题的非凸可行域Coriginal经SOC松弛转化为凸可行域CSOC，如果由二阶锥松弛转化模型求解得到的最优解S还是Coriginal中的点，则二阶锥松弛转化是严格有效的，即最优解S也是原问题的最优解，其中所涉及的二阶锥实际上为旋转二阶锥，其标准形式为：
因为显然式(14)满足式(21)中x1,x2≥0的条件，且形式完全一致，微网系统结构为辐射状电力网络，对于辐射状电力网络，根据其建立的二阶锥松弛模型是严格准确的，即由二阶锥松弛模型计算得出的最优解也是原问题的最优解，同时，若微网运行的优化调度目标函数是凸函数，并且是支路电流的严格递增函数，则二阶锥松弛转化后的模型是准确收敛的；
对于微网系统，其系统优化成本目标是电网支路电流的严格递增函数，分析如下：
在满足式(14)的情况下，微电网中节点电压由于电压上下限的约束而基本不变(电压损耗为±5％)，而支路电流/>(实际为支路电流的平方)在满足电流上下限的约束下可以有大范围的取值，但当/>与支路电流实际值的差值越大，则微电网支路上的损耗/>增加也越大，从而使得电源需要提供更多的有功功率来满足式(16)中的有功功率平衡约束，导致系统成本目标增加；
此时微网系统成本目标是支路电流的严格递增函数，在追求系统成本目标最低的求解过程中，使得不断向式(14)的等号处进行收敛趋近，最终使得/>等于支路电流实际值，因而对电网潮流约束进行二阶锥松弛转化得到的系统模型是有效和准确的；
当微网中的风力发电机出现弃风时，由转化模型求解得出的微电网各支路电流经过验算与实际值有差异，并且弃风量也与实际值有差异，而在微网无弃风时，与实际值无差异，此时具有准确性；
在微网无弃风的情况下，上面关于目标函数是支路电流的严格递增函数的结论是正确的，然而在微网系统存在弃风的情况下，当/>大于实际值时，微电网支路上的损耗增加，进一步电源出力增加，但是此时风机作为电源，增加的支路损耗功率是用弃风功率来补偿的，因而这部分增加的损耗功率在求解成本目标函数最优的过程中不仅不会增加系统成本，反而还造成了一种微网系统弃风量减少的错误假象，从而使得求解二阶锥松弛转化模型得出的微网系统弃风功率减小，支路电流增大；
所以当微网系统出现弃风时，式(14)和式(16)的约束不足让支路电流收敛至实际值，因为在此时目标成本函数不再是/>的严格递增函数，为解决上述问题，需要在目标成本函数中另外增加关于/>的严格递增函数，而显然微电网支路有功损耗即为/>的严格递增函数，并且将支路损耗成本加入到目标函数也是有实际意义的。式(22)为增加的微电网支路有功损耗成本函数：
式中KLine表示支路损耗成本系数，该系数值很小，其作用是保证在任何时候目标成本函数都是微电网支路电流的严格递增函数，同时又几乎不会影响微网系统的优化调度；
在加入微电网支路损耗成本函数到目标成本函数中后，对微网系统的优化调度问题再次进行了求解计算，可以验证提出的在目标函数中增加线路损耗的方法可以使微网系统在没有满足二阶锥规划应用条件的情况下正常求解，使总成本最优的同时，提高系统的适应性和灵活性。
本发明的有益效果是：
1、本发明技术方案中，基于二阶锥松弛转化方法，建立电-气综合能源微网系统，提出电-气综合能源微网优化调度策略，在满足运行约束的情况下最优化微网系统的经济性。
2、引入氢储能系统，氢能作为一种清洁能源在微网系统中的应用正越来越广泛，考虑以氢能为主导的综合能源系统具有环保性和经济性。
3、提出使用二阶锥松弛转化方法对电网潮流中的非线性约束进行线性化，将非线性规划问题转化为SOCP规划问题进行求解，在微网系统没有满足二阶锥规划应用条件的情况下，提出通过增加线路损耗以解决问题的方法，使总成本最优的同时，提高系统的适应性和灵活性。
附图说明
图1是电-气综合能源微网结构图；
图2是微网系统结构图；
图3是微网负荷柱状图；
图4是风电预测出力曲线图；
图5是电力价格梯形图；
图6是无弃风时线路总损耗及弃风量曲线图；
图7是无弃风时支路L12和L36电流曲线图；
图8是有弃风时线路总损耗及弃风量曲线图；
图9是有弃风时支路L12和L36电流曲线图；
图10是改进后有弃风时线路总损耗及弃风量曲线图；
图11是改进后有弃风时支路L12和L36电流曲线图。
具体实施方式
下面结合附图对本发明作进一步说明。
参照图1～图11，一种基于二阶锥松弛转化方法的电-气综合能源微网优化调度方法，包括以下步骤：
S1：建立电-气综合能源微网系统及内部设备模型，设定设备运行相关参数，将包含运行成本和环境成本的综合运行成本设定为目标函数；
电-气综合能源微网系统涉及电力网络、氢储能系统、天然气网络三部分，微网系统运行的优化目标为综合运行成本最低，综合运行成本包括系统运行成本以及环境成本；
系统运行成本表示为：
f1＝f11+f12+f13 (1)
式(1)中f11为微网系统与主网能量交换的成本与收益，表示为：
其中Ce1为微网通过主网购电分时电价，而Ce2为微网向主网售电分时电价；表示t时段微网通过主网购电的有功功率，/>表示t时段微网向主网售电的有功功率；
式(1)中f12为微网向天然气网络购售天然气的成本与收益，表示为：
其中Cg1表示天然气购买的价格，Cg2表示天然气出售价格；表示t时段微网从气网购买的CH4功率，/>表示t时段微网出售的CH4功率；H表示天然气高热值；
式(1)中f13为微网系统运行能量损失惩罚成本，表示为：
其中KLoss为权重系数，表示系统对于能量损失的重视程度，0≤KLoss≤1；和/>分别表示微网系统中电解槽、燃料电池、氢转气装置和微型燃气轮机MT在能量转化时的能量损失；
为简化处理，在此微网系统中各种能量形式的转化均认为不会产生有害气体，因此环境成本由CO2排放量构成，环境成本表达式为：
其中是该区域电网单位供电平均CO2排放系数，/>是MT发电二氧化碳排放系数，/>是CH4的二氧化碳捕捉系数；
结合以上数学表达式，考虑经济性与环保性的综合运行成本目标函数为：
式中f1 max与分别为微网系统最大运行成本和最大环境成本，ω1与ω2为系统运行成本和环境成本各自的权重系数；通过标幺值方法，根据重视程度对系统运行成本和环境成本赋予不同的权重系数，将多目标优化问题转化为单目标问题；
S2：通过二阶锥松弛方法将非线性的潮流约束线性化，从而将非线性规划问题转化为SOCP规划问题进行求解，微电网潮流基本模型表示如下：
式中：p、q为节点注入有功和无功功率，P、Q为支路潮流有功和无功功率，V为节点电压，I为支路电流；
在式(8)-(13)中，支路阻抗rij、xij，节点对地导纳gj、bj均为常数，由于微网系统内部电压等级为400V，因此对地导纳gj、bj对微电网潮流的影响很小，为简便运算，可将式(8)中的与/>进行忽略。节点电压V、支路电流I、支路潮流P、Q以及节点注入功率p、q为优化变量，显然它们之间为非线性关系，式(8)-(12)为一般电网潮流所共有的约束方程，式(13)为每个微电网系统所特有的约束条件，由微网自身的电力负荷、微源及其它微网内部结构决定；
通过二阶锥松弛方法将模型中的非线性约束线性化，从而将系统的非线性规划问题转化为二阶锥规划(second-order cone program,SOCP)，对式(10)进行转化处理，首先定义并进行二阶锥松弛转化，得到：
式(14)可以等价转化为2型范数表达的二阶锥松弛的标准形式，再将代入式(8)-(9)，即可得到通过二阶锥松弛转化后的最优潮流约束方程：
S3：在微网系统没有满足二阶锥规划应用条件的情况下，提出通过增加线路损耗以解决问题的方法；
原问题的非凸可行域Coriginal经SOC松弛转化为凸可行域CSOC，如果由二阶锥松弛转化模型求解得到的最优解S还是Coriginal中的点，则二阶锥松弛转化是严格有效的，即最优解S也是原问题的最优解，其中所涉及的二阶锥实际上为旋转二阶锥，其标准形式为：
因为显然式(14)满足式(21)中x1,x2≥0的条件，且形式完全一致，微网系统结构为辐射状电力网络，对于辐射状电力网络，根据其建立的二阶锥松弛模型是严格准确的，即由二阶锥松弛模型计算得出的最优解也是原问题的最优解，同时，若微网运行的优化调度目标函数是凸函数，并且是支路电流的严格递增函数，则二阶锥松弛转化后的模型是准确收敛的；
对于微网系统，其系统优化成本目标是电网支路电流的严格递增函数，分析如下：
在满足式(14)的情况下，微电网中节点电压由于电压上下限的约束而基本不变(电压损耗为±5％)，而支路电流/>(实际为支路电流的平方)在满足电流上下限的约束下可以有大范围的取值，但当/>与支路电流实际值的差值越大，则微电网支路上的损耗/>增加也越大，从而使得电源需要提供更多的有功功率来满足式(16)中的有功功率平衡约束，导致系统成本目标增加；
此时微网系统成本目标是支路电流的严格递增函数，在追求系统成本目标最低的求解过程中，使得不断向式(14)的等号处进行收敛趋近，最终使得/>等于支路电流实际值，因而对电网潮流约束进行二阶锥松弛转化得到的系统模型是有效和准确的；
当微网中的风力发电机出现弃风时，由转化模型求解得出的微电网各支路电流经过验算与实际值有差异，并且弃风量也与实际值有差异，而在微网无弃风时，与实际值无差异，此时具有准确性；
在微网无弃风的情况下，上面关于目标函数是支路电流的严格递增函数的结论是正确的，然而在微网系统存在弃风的情况下，当/>大于实际值时，微电网支路上的损耗增加，进一步电源出力增加，但是此时风机作为电源，增加的支路损耗功率是用弃风功率来补偿的，因而这部分增加的损耗功率在求解成本目标函数最优的过程中不仅不会增加系统成本，反而还造成了一种微网系统弃风量减少的错误假象，从而使得求解二阶锥松弛转化模型得出的微网系统弃风功率减小，支路电流增大；
所以当微网系统出现弃风时，式(14)和式(16)的约束不足让支路电流收敛至实际值，因为在此时目标成本函数不再是/>的严格递增函数，为解决上述问题，需要在目标成本函数中另外增加关于/>的严格递增函数，而显然微电网支路有功损耗即为/>的严格递增函数，并且将支路损耗成本加入到目标函数也是有实际意义的。式(22)为增加的微电网支路有功损耗成本函数：
式中KLine表示支路损耗成本系数，该系数值很小，其作用是保证在任何时候目标成本函数都是微电网支路电流的严格递增函数，同时又几乎不会影响微网系统的优化调度；
在加入微电网支路损耗成本函数到目标成本函数中后，对微网系统的优化调度问题再次进行了求解计算，可以验证提出的在目标函数中增加线路损耗的方法可以使微网系统在没有满足二阶锥规划应用条件的情况下正常求解，使总成本最优的同时，提高系统的适应性和灵活性。
优化结果对比分析：
为了验证上述二阶锥松弛转化后的模型是有效性和准确性，本文对所建立的微网系统进行了仿真算例分析，时间步长Δt取为1h，调度时间T为24h。
图3为微网负荷柱状图，图4为风电预测出力曲线图，图5为电力价格梯形图。
图3中包含调度时间T内电负荷、气负荷和氢负荷的数值。风机在一天中的预测出力如图4所示，可以看到在白天风机出力较低，在夜晚风机出力较高。图5为调度时间T内的分时电价，包括微网通过主网的购电电价和微网向主网的售电电价，此外天然气价格给定为2.5/m3。
图6为无弃风时线路总损耗及弃风量曲线图，图7为无弃风时支路L12和L36电流曲线图。
由图6和图7可知，在微网无弃风的情况下，由二阶锥转化后的模型求解所得到的微电网线路总损耗、支路电流与实际值完全一致，此时二阶锥转化方法具有有效性和准确性，所以关于目标函数是支路电流的严格递增函数的结论也是完全正确的。
图8为有弃风时线路总损耗及弃风量曲线图，图9为有弃风时支路L12和L36电流曲线图。
当风机出力为正常情况下的1.5倍时，微网系统出现弃风。由图8和图9可知，在0:00-7:00和23:00-24:00时，微网出现弃风，由转化模型求解得出的微电网各支路电流经过验算与实际值有差异，并且弃风量也与实际值有差异。从图8可以看到由二阶锥松弛转化模型求得微网的弃风量小于实际弃风量，并且差值较大，而从图9看到由二阶锥松弛转化模型求得微网支路电流大于实际支路电流，并且差值也较大。分析可以得到的结果，在夜间时段，风机出力增大，微网出现弃风的时段线路总损耗、支路电流值与实际值有很大差异，而在白天时段，微网无弃风时，与实际值无差异，此时具有准确性。
图10为改进后有弃风时线路总损耗及弃风量曲线图，图11为改进后有弃风时支路L12和L36电流曲线图；
从图10和图11中可以看出在使用改进二阶锥松弛转化方法，在目标函数中加入线路损耗成本后，微网系统在0:00-7:00和23:00-24:00存在弃风时，二阶锥松弛转化模型的求解结果仍为实际有效值，验证了该方法的有效性。
在本说明书的描述中，参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本发明的至少一个实施例或示例中。在本说明书中，对上述术语的示意性表述不必须针对的是相同的实施例或示例。而且，描述的具体特征、结构、材料或者特点可以在任一个或多个实施例或示例中以合适的方式结合。此外，本领域的技术人员可以将本说明书中描述的不同实施例或示例进行结合和组合。
尽管上面已经展示出和描述了本发明的实施例，可以理解的是，本说明书实施例所述的内容仅仅是对发明构思的实现形式的列举，本发明的保护范围不应当被视为仅限于实施例所陈述的具体形式，本发明的保护范围也包括本领域技术人员根据本发明构思所能够想到的等同技术手段。
Claims (2)
1.一种基于二阶锥松弛转化方法的电-气综合能源微网优化调度方法，其特征在于，所述方法包括以下步骤：
S1：建立电-气综合能源微网系统及内部设备模型，设定设备运行相关参数，将包含运行成本和环境成本的综合运行成本设定为目标函数；
S2：通过二阶锥松弛方法将非线性的潮流约束线性化，从而将非线性规划问题转化为二阶锥规划问题进行求解；
S3：在微网系统没有满足二阶锥规划应用条件的情况下，提出通过增加线路损耗以解决问题的方法；
所述步骤S2中，微电网潮流基本模型表示如下：
式中：p、q为节点注入有功和无功功率，P、Q为支路潮流有功和无功功率，V为节点电压，V、分别表示节点电压的最大值和最小值，I为支路电流；
在式(8)-(13)中，支路阻抗rij、xij，节点对地导纳gj、bj均为常数，由于微网系统内部电压等级为400V，因此对地导纳gj、bj对微电网潮流的影响很小，为简便运算，将式(8)中的与/>进行忽略，节点电压V、支路电流I、支路潮流P、Q以及节点注入功率p、q为优化变量，它们之间为非线性关系，式(8)-(12)为一般电网潮流所共有的约束方程，式(13)为每个微电网系统所特有的约束条件，由微网自身的电力负荷、微源及其它微网内部结构决定；
通过二阶锥松弛方法将模型中的非线性约束线性化，从而将系统的非线性规划问题转化为二阶锥规划，对式(10)进行转化处理，首先定义并进行二阶锥松弛转化，得到：
式(14)等价转化为2型范数表达的二阶锥松弛的标准形式，再将 代入式(8)-(9)，即得到通过二阶锥松弛转化后的最优潮流约束方程：
所述步骤S3中，原问题的非凸可行域Coriginal经二阶锥松弛转化为凸可行域CSOC，如果由二阶锥松弛转化模型求解得到的最优解S还是Coriginal中的点，则二阶锥松弛转化是严格有效的，即最优解S也是原问题的最优解，其中所涉及的二阶锥为旋转二阶锥，其标准形式为：
因为显然式(14)满足式(21)中x1,x2≥0的条件，且形式完全一致，微网系统结构为辐射状电力网络，对于辐射状电力网络，根据其建立的二阶锥松弛模型是严格准确的，即由二阶锥松弛模型计算得出的最优解也是原问题的最优解，同时，若微网运行的优化调度目标函数是凸函数，并且是支路电流的严格递增函数，则二阶锥松弛转化后的模型是准确收敛的；
对于微网系统，其系统优化成本目标是电网支路电流的严格递增函数，分析如下：
在满足式(14)的情况下，微电网中节点电压由于电压上下限的约束而基本不变，而支路电流/>在满足电流上下限的约束下有大范围的取值，但当/>与支路电流实际值的差值越大，则微电网支路上的损耗/>增加也越大，从而使得电源需要提供更多的有功功率来满足式(16)中的有功功率平衡约束，导致系统成本目标增加；
此时微网系统成本目标是支路电流的严格递增函数，在追求系统成本目标最低的求解过程中，使得不断向式(14)的等号处进行收敛趋近，最终使得/>等于支路电流实际值；
当微网中的风力发电机出现弃风时，由转化模型求解得出的微电网各支路电流经过验算与实际值有差异，并且弃风量也与实际值有差异，而在微网无弃风时，与实际值无差异，此时具有准确性；
在微网无弃风的情况下，上面关于目标函数是支路电流的严格递增函数的结论是正确的，然而在微网系统存在弃风的情况下，当/>大于实际值时，微电网支路上的损耗/>增加，进一步电源出力增加，但是此时风机作为电源，增加的支路损耗功率是用弃风功率来补偿的，因而这部分增加的损耗功率在求解成本目标函数最优的过程中不仅不会增加系统成本，反而还造成了一种微网系统弃风量减少的错误假象，从而使得求解二阶锥松弛转化模型得出的微网系统弃风功率减小，支路电流增大；
所以当微网系统出现弃风时，式(14)和式(16)的约束不足让支路电流收敛至实际值，因为在此时目标成本函数不再是/>的严格递增函数，为解决上述问题，需要在目标成本函数中另外增加关于/>的严格递增函数，而显然微电网支路有功损耗即为/>的严格递增函数，并且将支路损耗成本加入到目标函数也是有实际意义的，式(22)为增加的微电网支路有功损耗成本函数：
式中KLine表示支路损耗成本系数，该系数值很小，其作用是保证在任何时候目标成本函数都是微电网支路电流的严格递增函数，同时不会影响微网系统的优化调度；
在加入微电网支路损耗成本函数到目标成本函数中后，对微网系统的优化调度问题再次进行了求解计算，验证提出的在目标函数中增加线路损耗的方法使微网系统在没有满足二阶锥规划应用条件的情况下正常求解，使总成本最优的同时，提高系统的适应性和灵活性。
2.如权利要求1所述的一种基于二阶锥松弛转化方法的电-气综合能源微网优化调度方法，其特征在于，所述步骤S1中，建立电-气综合能源微网系统及内部设备模型，该系统涉及电力网络、氢储能系统、天然气网络三部分，微网系统运行的优化目标为综合运行成本最低，综合运行成本包括系统运行成本以及环境成本；
系统运行成本表示为：
f1＝f11+f12+f13 (1)
式(1)中f11为微网系统与主网能量交换的成本与收益，表示为：
其中Ce1为微网通过主网购电分时电价，而Ce2为微网向主网售电分时电价；表示t时段微网通过主网购电的有功功率，/>表示t时段微网向主网售电的有功功率；
式(1)中f12为微网向天然气网络购售天然气的成本与收益，表示为：
其中Cg1表示天然气购买的价格，Cg2表示天然气出售的价格；表示t时段微网从气网购买的CH4功率，/>表示t时段微网出售的CH4功率；H表示天然气高热值；
式(1)中f13为微网系统运行能量损失惩罚成本，表示为：
其中KLoss为权重系数，表示系统对于能量损失的重视程度，0≤KLoss≤1；和/>分别表示微网系统中电解槽、燃料电池、氢转气装置和微型燃气轮机MT在能量转化时的能量损失；
为简化处理，在此微网系统中各种能量形式的转化均认为不会产生有害气体，因此环境成本由CO2排放量构成，环境成本表达式为：
其中是该区域电网单位供电平均CO2排放系数，/>是MT发电二氧化碳排放系数，是CH4的二氧化碳捕捉系数；
结合以上数学表达式，考虑经济性与环保性的综合运行成本目标函数为：
式中f1 max与分别为微网系统最大运行成本和最大环境成本，ω1与ω2为系统运行成本和环境成本各自的权重系数；通过标幺值方法，根据重视程度对系统运行成本和环境成本赋予不同的权重系数，将多目标优化问题转化为单目标问题。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110752294.3A CN113570117B (zh) | 2021-07-02 | 2021-07-02 | 一种基于二阶锥松弛转化方法的电-气综合能源微网优化调度方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110752294.3A CN113570117B (zh) | 2021-07-02 | 2021-07-02 | 一种基于二阶锥松弛转化方法的电-气综合能源微网优化调度方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113570117A CN113570117A (zh) | 2021-10-29 |
| CN113570117B true CN113570117B (zh) | 2024-02-09 |
Family
ID=78163652
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110752294.3A Active CN113570117B (zh) | 2021-07-02 | 2021-07-02 | 一种基于二阶锥松弛转化方法的电-气综合能源微网优化调度方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113570117B (zh) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115015474B (zh) * | 2022-04-15 | 2023-06-23 | 广东电网有限责任公司 | 一种电力用户的碳排放量检测方法及设备 |
| CN115713438A (zh) * | 2022-11-25 | 2023-02-24 | 国网江苏省电力有限公司常州供电分公司 | 综合能源系统的优化调度方法 |
| CN115859686B (zh) * | 2023-02-07 | 2023-05-09 | 山东科技大学 | 考虑扩展碳排放流的综合能源系统低碳调度方法及系统 |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109242366A (zh) * | 2018-11-06 | 2019-01-18 | 国电南瑞科技股份有限公司 | 一种电-气互联综合能源系统的多时段潮流优化方法 |
| CN109861232A (zh) * | 2019-02-22 | 2019-06-07 | 国家电网有限公司 | 一种基于二阶锥松弛方法的配电网动态重构降损方法 |
| CN110263387A (zh) * | 2019-05-29 | 2019-09-20 | 国网福建省电力有限公司 | 基于能源集线器的配电网、天然气网建立的综合能源系统优化运行模型及其线性处理方法 |
| CN110322051A (zh) * | 2019-06-06 | 2019-10-11 | 国网浙江省电力有限公司经济技术研究院 | 考虑n-1安全约束的综合能源系统优化配置方法 |
| CN111652441A (zh) * | 2020-06-04 | 2020-09-11 | 四川大学 | 考虑气电联合需求响应的气电综合能源系统配网优化方法 |
| WO2020237700A1 (zh) * | 2019-05-28 | 2020-12-03 | 国电南瑞科技股份有限公司 | 一种多能源系统运行调度方法 |
-
2021
- 2021-07-02 CN CN202110752294.3A patent/CN113570117B/zh active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109242366A (zh) * | 2018-11-06 | 2019-01-18 | 国电南瑞科技股份有限公司 | 一种电-气互联综合能源系统的多时段潮流优化方法 |
| CN109861232A (zh) * | 2019-02-22 | 2019-06-07 | 国家电网有限公司 | 一种基于二阶锥松弛方法的配电网动态重构降损方法 |
| WO2020237700A1 (zh) * | 2019-05-28 | 2020-12-03 | 国电南瑞科技股份有限公司 | 一种多能源系统运行调度方法 |
| CN110263387A (zh) * | 2019-05-29 | 2019-09-20 | 国网福建省电力有限公司 | 基于能源集线器的配电网、天然气网建立的综合能源系统优化运行模型及其线性处理方法 |
| CN110322051A (zh) * | 2019-06-06 | 2019-10-11 | 国网浙江省电力有限公司经济技术研究院 | 考虑n-1安全约束的综合能源系统优化配置方法 |
| CN111652441A (zh) * | 2020-06-04 | 2020-09-11 | 四川大学 | 考虑气电联合需求响应的气电综合能源系统配网优化方法 |
Non-Patent Citations (1)
| Title |
|---|
| 考虑氢能-天然气混合储能的电-气综合能源微网日前经济调度优化;刘继春;周春燕;高红均;郭焱林;朱雨薇;;电网技术(第01期);全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113570117A (zh) | 2021-10-29 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN113570117B (zh) | 一种基于二阶锥松弛转化方法的电-气综合能源微网优化调度方法 | |
| Hannan et al. | Binary particle swarm optimization for scheduling MG integrated virtual power plant toward energy saving | |
| Wang et al. | Day-ahead optimal economic dispatching of integrated port energy systems considering hydrogen | |
| CN105375479B (zh) | 一种基于模型预测控制的分布式能源能量管理方法 | |
| Zhou et al. | Integrated energy production unit: An innovative concept and design for energy transition toward low-carbon development | |
| Zhou et al. | Optimal sizing design and integrated cost-benefit assessment of stand-alone microgrid system with different energy storage employing chameleon swarm algorithm: A rural case in Northeast China | |
| Moazzami et al. | Application of multi-objective grey wolf algorithm on energy management of microgrids with techno-economic and environmental considerations | |
| Rouhani et al. | A comprehensive method for optimum sizing of hybrid energy systems using intelligence evolutionary algorithms | |
| Zhang et al. | Optimal scheduling strategy of virtual power plant with power-to-gas in dual energy markets | |
| Wang et al. | Study on the optimal configuration of a wind-solar-battery-fuel cell system based on a regional power supply | |
| Suo et al. | New energy wide area complementary planning method for multi-energy power system | |
| Zhu et al. | Multi-objective optimization of a hybrid energy system integrated with solar-wind-PEMFC and energy storage | |
| CN115408808A (zh) | 一种计及碳排放的能源配网系统多时间尺度优化调度方法 | |
| Gbadega et al. | JAYA algorithm-based energy management for a grid-connected micro-grid with PV-wind-microturbine-storage energy system | |
| Liu et al. | Low-carbon economic dispatch of integrated electricity-gas energy system considering carbon capture, utilization and storage | |
| Wang et al. | Flexible resource dynamic aggregation regulation method of virtual power plant to ensure more renewable energy generation | |
| CN109474015A (zh) | 基于两层优化的主动配电网储能系统优化配置方法及系统 | |
| Belkhier et al. | Novel design and adaptive coordinated energy management of hybrid fuel‐cells/tidal/wind/PV array energy systems with battery storage for microgrids | |
| Yasin et al. | Fuzzy logic based management of a stand-alone hybrid generator | |
| Zhang et al. | Optimal power flow algorithm based on second-order cone relaxation method for electricity-gas integrated energy microgrid | |
| Rouhani et al. | A teaching learning based optimization for optimal design of a hybrid energy system | |
| Phan-Van et al. | A comparison of different metaheuristic optimization algorithms on hydrogen storage-based microgrid sizing | |
| Li et al. | Optimal configuration for distributed generations in micro-grid system considering diesel as the main control source | |
| He et al. | Capacity Configuration in Integrated Energy Production Unit Considering Ladder-Type Trading Under CCER Quota | |
| Wei et al. | Optimal generation planning in a micro-grid for supplying electrical and thermal loads in order to reduce pollutant emissions |
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 |