CN113572163B - 基于改进仿射算法的分布式电源主动配电网优化调度方法 - Google Patents
基于改进仿射算法的分布式电源主动配电网优化调度方法 Download PDFInfo
- Publication number
- CN113572163B CN113572163B CN202110859394.6A CN202110859394A CN113572163B CN 113572163 B CN113572163 B CN 113572163B CN 202110859394 A CN202110859394 A CN 202110859394A CN 113572163 B CN113572163 B CN 113572163B
- Authority
- CN
- China
- Prior art keywords
- value
- interval
- node
- formula
- voltage
- 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
- PXFBZOLANLWPMH-UHFFFAOYSA-N 16-Epiaffinine Natural products C1C(C2=CC=CC=C2N2)=C2C(=O)CC2C(=CC)CN(C)C1C2CO PXFBZOLANLWPMH-UHFFFAOYSA-N 0.000 title claims abstract description 64
- 238000000034 method Methods 0.000 title claims abstract description 46
- 238000004422 calculation algorithm Methods 0.000 title claims abstract description 36
- 238000004364 calculation method Methods 0.000 claims abstract description 33
- 230000006870 function Effects 0.000 claims abstract description 25
- 238000005457 optimization Methods 0.000 claims abstract description 17
- 238000011156 evaluation Methods 0.000 claims abstract description 7
- 238000004146 energy storage Methods 0.000 claims description 26
- 239000000243 solution Substances 0.000 claims description 20
- 230000005611 electricity Effects 0.000 claims description 15
- 238000012937 correction Methods 0.000 claims description 14
- 230000007704 transition Effects 0.000 claims description 12
- 230000009194 climbing Effects 0.000 claims description 9
- 239000011159 matrix material Substances 0.000 claims description 6
- 238000011144 upstream manufacturing Methods 0.000 claims description 6
- 238000011161 development Methods 0.000 claims description 5
- 230000005540 biological transmission Effects 0.000 claims description 3
- 238000006243 chemical reaction Methods 0.000 claims description 3
- 238000002347 injection Methods 0.000 claims description 3
- 239000007924 injection Substances 0.000 claims description 3
- 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 2
- 229910052739 hydrogen Inorganic materials 0.000 claims 1
- 239000001257 hydrogen Substances 0.000 claims 1
- 125000004435 hydrogen atom Chemical class [H]* 0.000 claims 1
- 101000712674 Homo sapiens TGF-beta receptor type-1 Proteins 0.000 description 7
- 102100033456 TGF-beta receptor type-1 Human genes 0.000 description 7
- 238000010586 diagram Methods 0.000 description 5
- 238000009825 accumulation Methods 0.000 description 3
- 230000008901 benefit Effects 0.000 description 3
- 230000006872 improvement Effects 0.000 description 3
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 2
- 229910052799 carbon Inorganic materials 0.000 description 2
- 238000005315 distribution function Methods 0.000 description 2
- 230000007246 mechanism Effects 0.000 description 2
- 238000010248 power generation Methods 0.000 description 2
- 101000813777 Homo sapiens Splicing factor ESS-2 homolog Proteins 0.000 description 1
- 102100039575 Splicing factor ESS-2 homolog Human genes 0.000 description 1
- 238000010521 absorption reaction Methods 0.000 description 1
- 230000033228 biological regulation Effects 0.000 description 1
- 230000000903 blocking effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 238000007599 discharging Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000003912 environmental pollution Methods 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 238000006386 neutralization reaction Methods 0.000 description 1
- 230000035699 permeability Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
- 230000002087 whitening effect Effects 0.000 description 1
Images
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
- H02J3/004—Generation forecast, e.g. methods or systems for forecasting future energy generation
-
- 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
- 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
- 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/04—Circuit arrangements for ac mains or ac distribution networks for connecting networks of the same frequency but supplied from different sources
- H02J3/06—Controlling transfer of power between connected networks; Controlling sharing of load between connected networks
-
- 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
- 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]
-
- 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/22—The renewable source being solar energy
- H02J2300/24—The renewable source being solar energy of photovoltaic origin
-
- 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
- Y02E10/00—Energy generation through renewable energy sources
- Y02E10/50—Photovoltaic [PV] energy
- Y02E10/56—Power conversion systems, e.g. maximum power point trackers
-
- 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
- Y02E10/00—Energy generation through renewable energy sources
- Y02E10/70—Wind energy
- Y02E10/76—Power conversion electric or electronic aspects
Landscapes
- Business, Economics & Management (AREA)
- Engineering & Computer Science (AREA)
- Human Resources & Organizations (AREA)
- Economics (AREA)
- Strategic Management (AREA)
- Tourism & Hospitality (AREA)
- Physics & Mathematics (AREA)
- Power Engineering (AREA)
- Theoretical Computer Science (AREA)
- Marketing (AREA)
- Entrepreneurship & Innovation (AREA)
- General Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- Operations Research (AREA)
- Development Economics (AREA)
- Quality & Reliability (AREA)
- Health & Medical Sciences (AREA)
- Game Theory and Decision Science (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
基于改进仿射算法的分布式电源主动配电网优化调度方法,包括以下步骤:1)基于灰色马尔科夫模型预测出未来一天光伏及风电的输出功率区间值;2)构建以运行费用最小、负荷波动最小的多目标优化模型,并采用TOPSIS法将多目标转换为单目标;3)建立辅助模型求解各决策变量及目标函数最优值的区间解的上、下界值;4)将优化求解得到的区间调度结果及各功率区间值转换为仿射形式;5)提出一种改进的仿射算法对配电网进行潮流计算,输出区间潮流结果;6)输出优化后的区间调度结果,并计算评价指标。本发明能精准预测风光出力区间,在降低潮流区间结果保守性的同时提高调度计算效率,并为调度人员提供数据参考。
Description
技术领域
本发明涉及主动配电网优化调度技术领域,具体涉及一种基于改进仿射算法的分布式电源主动配电网优化调度方法。
背景技术
全球化石能源枯竭、环境污染日益严重的现实挑战下,大力推动可再生能源发展、加快传统能源结构的转变,实现“碳达峰”和“碳中和”已成为我国重要战略目标。由于可再生能源的不确定性和间歇性,多样式的分布式可再生能源和储能系统接入配电网,随之引起配电网电能质量、继电保护、负荷波动等诸多问题。传统配电网逐渐转变为主动配电网,主动配电网可以通过调控手段自主调节分布式发电(DG)、储能装置(ESS)、柔性负荷(FL)等灵活性资源,并根据配电网的实际运行状态调节其网络结构、发电单元及负荷,来达到配电网的安全和经济运行等目标。因此,为了减轻对配电网运行安全的影响,有必要研究含高渗透率分布式可再生能源的主动配电网优化调度问题,而其关键是构建科学合理的调度模型和提出切实可行的求解计算方法。
中国专利“基于风、光随机性的主动配电网多目标日前优化调度方法”(申请号202110017757.1)公开了一种基于风、光随机性的主动配电网多目标日前优化调度方法,采用概率分布函数表征风光出力,并通过随机数模拟风光出力的随机性;该方法存在概率分布函数获取较难、不够准确等问题。而区间分析法可以表示所有风光出力的可能值,具有建模的简单性和结果的完备性等优点。仿射算术作为对区间算术的改进,优化了区间计算中区间扩张的问题,并且仿射算术能够考虑不同的不确定性变量间的关联性和耦合性,在电力系统潮流计算中越来越受到人们的关注。因此,本发明提出一种基于改进仿射算法的分布式电源主动配电网优化调度方法。
发明内容
为解决上述技术问题,本发明提供一种基于改进仿射算法的分布式电源主动配电网优化调度方法,能精准预测风光出力区间,在降低潮流区间结果保守性的同时提高调度计算效率,并为调度人员提供数据参考。
本发明采取的技术方案为:
基于改进仿射算法的分布式电源主动配电网优化调度方法,包括以下步骤:
步骤一:基于灰色马尔科夫模型,得到风机和光伏的区间预测值;
步骤二:构建以运行费用最小、负荷波动最小的多目标优化模型,并采用TOPSIS法将多目标转换为单目标;
步骤三:建立辅助模型,求解各决策变量及目标函数最优值的区间解的上、下界值;
步骤四:将优化求解得到的区间调度结果、以及各功率区间值转换为仿射形式;
步骤五:基于改进仿射算法对配电网进行潮流计算,输出区间潮流结果;
步骤六:输出优化后的区间调度结果,并计算评价指标。
通过上述步骤完成主动配电网优化调度。
所述步骤一包括以下步骤:
S1.1:采用灰色GM(1,1)模型进行预测:
式(9)中:n为样本数量,a、b分别为发展系数和灰作用量,x(0)(1)为原始数据序列,x(0)(k+1) 为灰色GM(1,1)模型预测值序列。
S1.2:引入误差修正模型,对预测值进行修正,修正公式如下所示:
式(10)中:X(0)为历史数据序列,X(2)为预测值序列,E(0)为误差修正序列,n为样本数量;
S1.3:对误差修正序列E(0)进行预测,得到E(1);再对预测得到的序列X(2)进行二次修正,得到预测修正后的序列XΔ(2)。
S1.4:采用马尔科夫模型,得到如下各状态的转移概率、以及状态转移概率矩阵:
式(11)、(12)中:qij(k)为状态Ei经过k步转移到状态Ej的概率,Mi为处于状态Ei的样本数,Mij(k)为状态Ei经过k步转移到状态Ej的样本数,i、j=1,…,n;q(k)为状态转移概率矩阵;
S1.5:根据q(k)计算各步的自相关系数和权重系数,最终确定出历史数据对应预测值的每种状态的概率qi,并结合预测修正后的序列XΔ(2),得到最终预测值X:
式(13)中:XΔ(2)为预测修正后的序列,qi为每种状态的转移概率,Ei,mid为第i个状态区间的中点值,X为最终预测值。
所述步骤二包括以下步骤:
S2.1:以配电网运行总费用最小,建立的第一目标函数为:
式(14)中:T为总时段数,Ness、Nw、Npv和Nint分别为储能、风机、光伏和可控负荷的数量;Pt g和Pt loss分别为t时段ADN的购电价格、购电量和有功网损;和分别为储能、风机和光伏的平准化度电成本;和分别为储能、风机、光伏和可控负荷在t时段的功率,且均为仿射值。
S2.2:以配电网负荷波动最小,建立的第二目标函数为:
式(15)中:Pt Load和Pav分别为ADN在t时段的总负荷和整个调度周期内负荷均值。
S2.3:采用TOPSIS法,将多目标转换为如下单目标:
所述步骤二还包括确定调度问题的基本约束条件:
1)、有功功率平衡约束:
式(17)中:Pt g和Pt loss分别是t时段AND的购电量和有功网损;和分别为风机、光伏、储能和可控负荷在t时段的功率,且均为仿射值;为t时段节点i 负荷的有功功率;Nw、Npv、Ness和Nint分别为风机、光伏、储能和可控负荷的数量;I为AND 的节点总数。
2)、区间潮流等式约束:
3)、DG出力与爬坡率约束:
式(19)、(20)中:和Ppvmin,np/Ppvmax,np分别为风机和光伏出力的最小/大值。-Rwmax,nw/Rwmax,nw、-Rpvmax,np/Rpvmax,np分别为风机、光伏爬坡率的最小/大值;分别为风机、光伏在t时段的功率;Rwt,nw、Rpvt,np分别为风机、光伏在t时段的爬坡率。
4)、储能系统约束:
式(21)、(22)中:-Pessmax,ne/Pessmax,ne为ESS充/放电功率最大值;SOCessmin,ne和SOCessmax,ne分别为ESS荷电状态的最小值和最大值;为ESS在t时刻的负荷;SOCesst,n为ESS在t时刻的荷电状态。
5)、配变关口功率约束:
6)、节点电压约束:
7)、支路有功功率约束:
8)、净负荷波动率约束:
所述步骤三中:建立第一辅助模型,求解决策变量的上、下限:
第一辅助模型:
式(29)、(30)中:Pt,n为第n个决策变量在t时刻的取值,h(P)=[hl,hu]为等式约束条件,g(P)≤[gl,gu]为不等式约束条件。
其中,将第一辅助模型求解的结果作为Pt,nl和Pt,nu,决策变量Pt,n满足约束条件的取值区间[Pt,nl,Pt,nu]。
建立第二辅助模型,求解目标函数最优值的区间解的上、下界值:
式(31)、(32)中:W(Pt,n)为决策变量Pt,n的目标函数,h(P)=[hl,hu]为等式约束条件, g(P)≤[gl,gu]为不等式约束条件。
其中,将第二辅助模型求解的结果,分别作为目标函数最优值的下界和上界值,并由相对应的各决策变量Pt,n的取值作为各自区间解的边界值,得出在目标函数取最优值时各可调度机组的出力区间。
所述步骤五包括如下步骤:
S5.1:将配电网所有节点进行编号,输入ADN各节点参数和各支路参数;
S5.3:从末节点开始计算每节点的注入电流:
其中,式(34)具体的计算方法如下:
式(36)中:x0、z0为仿射中心值,xi、zj为仿射系数值,εn+1为新增噪声元。
S5.3.3:根据KCL计算上游节点电流,具体计算公式如下:
式(37)中:Ni为i节点及下游节点的集合,j为上游节点。
S5.4:前推计算每个节点的电压:
S5.6:计算配电网有功网损:
S5.7:输出区间潮流结果。
所述步骤六,计算评价指标包括:
①:电压稳定性指标FVSI:
式(40)中:FVSIij为支路ij电压稳定性指标;Vi为节点i的电压;Vj为节点j的电压。
②:线路容量裕度指标FML:
③:算法改进精度指标CV:
本发明一种基于改进仿射算法的分布式电源主动配电网优化调度方法,技术效果如下: 1.本发明提出的改进后的灰色马尔科夫模型能够较为精准的预测出各类可再生能源的出力区间;
2.本发明提出的通过区间泰勒公式改进后的仿射除法形式简洁,由于固定了噪声元个数,区间计算结果的保守性更小,缩短了计算时间,增大了计算效率;
3.本发明将改进的仿射算法应用于含高渗透率分布式电源主动配电网区间优化调度中,区间优化调度结果可为调度人员提供更加直观的主动配电网状态量上下界信息,相较于传统放射算法,本发明所提方法得到的结果更精确、保守性更好。
4.本发明提出了一种带有误差修正机制的灰色马尔科夫模型,预测了风机和光伏的功率区间值,并建立了综合考虑主动配电网运行约束和灵活性指标、以综合运行费用最低和净负荷波动最小为目标的多目标区间优化调度数学模型。在此基础上,提出改进仿射算法,将区间泰勒公式应用于仿射计算的非线性运算中。通过算例仿真可得,本发明所提方法计算效率更高、区间结果的保守性更好,有效改善了区间结果的可靠性和有效性。
附图说明
下面结合附图和实施例,对本发明的技术方案作进一步说明:
图1是本发明提供的基于改进仿射算法的含高渗透率分布式电源主动配电网优化调度方法流程图;
图2是本发明提供的改进的IEEE33节点拓扑图;
图3是本发明提供的通过灰色马尔科夫模型得到的光伏及风机功率区间预测结果图;
图4(1)是本发明提供的区间调度结果的区间中点图;
图4(2)是本发明提供的区间调度结果的区间半径图;
图5是本发明提供的储能电池1的荷电状态区间轨迹和边界曲线;
图6是本发明提供的IEEE33节点电压区间图;
图7是本发明提供的电压稳定性指标结果图;
图8是本发明提供的线路容量裕度指标结果图;
图9是本发明提供的算法改进精度示意图;
图10是本发明提供的20:00时节点电压区间计算结果对比图。
具体实施方式
图1是所提基于改进仿射算法的含高渗透率分布式电源主动配电网优化调度方法流程图。图1中,本发明提供的基于改进仿射算法的含高渗透率分布式电源主动配电网优化调度方法包括以下步骤:
步骤一:基于灰色马尔科夫模型,预测出未来一天光伏及风电的输出功率区间值;
步骤二:构建以运行费用最小、负荷波动最小的多目标优化模型,并采用TOPSIS法将多目标转换为单目标;
步骤三:建立辅助模型,求解各决策变量及目标函数最优值的区间解的上、下界值;
步骤四:将优化求解得到的区间调度结果及各功率区间值转换为仿射形式;
步骤五:提出一种改进的仿射算法,对配电网进行潮流计算,输出区间潮流结果;
步骤六:输出优化后的区间调度结果,并计算评价指标。
所述步骤一,采用带有误差修正机制的灰色马尔科夫模型得到风机和光伏的区间预测值,包括以下步骤:
生成灰色GM(1,1)模型序列如下:
X(0)={x(0)(1),x(0)(2),…,x(0)(n)} (1)
X(1)={x(1)(1),x(1)(2),…,x(1)(n)} (2)
Z(1)={z(1)(2),z(1)(3),…,z(1)(n)} (3)
式中:X(0)为原始数据序列、X(1)为累加生成序列、Z(1)为临近生成序列,n为样本数量。
建立白化微分方程如下:
微分方程的解如下:
式中:n为样本数量,a、b分别为发展系数和灰作用量,x(0)(1)为原始数据序列,x(0)(k+1) 为灰色GM(1,1)模型预测值序列。
将得到的累加生成序列的预测值进行一次累减,得到第k+1个数据的预测值:
式中:n为历史数据的总个数,x(0)(k+1)为灰色GM(1,1)模型预测得到的数据。
引入误差修正模型对预测值进行修正,修正公式如下所示:
式中:X(0)为历史数据序列,X(2)为预测得到的序列,E(0)为误差修正序列,n为样本数量。
对误差修正序列E(0)进行预测得到E(1),再对预测序列X(2)进行二次修正得到预测修正后的序列XΔ(2)。
采用马尔科夫模型,得到如下各状态的转移概率及状态转移概率矩阵:
式(11)、(12)中:q(k)为状态转移概率矩阵,qij(k)为状态Ei经过k步转移到状态Ej的概率,Mi为处于状态Ei的样本数,Mij(k)为状态Ei经过k步转移到状态Ej的样本数,i、 j=1,…,n。
根据q(k)计算各步的自相关系数和权重系数,最终确定出历史数据对应预测值的每种状态的概率qi,并结合XΔ(2)得到最终预测值X:
式中:XΔ(2)为预测修正后的序列,qi为每种状态的转移概率,Ei,mid为第i个状态区间的中点值,X为最终预测值。
所述步骤二中,以分布式电源,储能系统为调度对象,以运行费用最小、负荷波动最小为多目标进行优化建模,采用TOPSIS法将多目标转换为单目标。包括以下步骤:
1)以配电网运行总费用最小建立的第一目标函数为:
式中:T为总时段数,Ness、Nw、Npv和Nint分别为储能、风机、光伏和可控负荷的数量; 和分别为t时段ADN的购电价格、购电量和有功网损;Lwnw和分别为储能、风机和光伏的平准化度电成本;和分别为储能、风机、光伏和可控负荷在t时段的功率,且均为仿射值。
2)以配电网负荷波动最小建立的第二目标函数为:
3)采用TOPSIS法将多目标转换为如下单目标:
所述步骤二中,确定调度问题的基本约束条件包括:
1)有功功率平衡约束:
式(17)中:和分别是t时段AND的购电量和有功网损;和分别为风机、光伏、储能和可控负荷在t时段的功率,且均为仿射值;为t时段节点i负荷的有功功率;Nw、Npv、Ness和Nint分别为风机、光伏、储能和可控负荷的数量;I为AND的节点总数。
2)区间潮流等式约束:
3)DG出力与爬坡率约束:
式(19)、(20)中:和Ppvmin,np/Ppvmax,np分别为风机和光伏出力的最小/大值。-Rwmax,nw/Rwmax,nw、-Rpvmax,np/Rpvmax,np分别为风机、光伏爬坡率的最小/大值;分别为风机、光伏在t时段的功率;Rwt,nw、Rpvt,np分别为风机、光伏在t时段的爬坡率。
4)储能系统约束:
式(21)、(22)中:-Pessmax,ne/Pessmax,ne为ESS充/放电功率最大值;SOCessmin,ne和SOCessmax,ne分别为ESS荷电状态的最小值和最大值;为ESS在t时刻的负荷;SOCesst,n为ESS在t时刻的荷电状态。
5)配变关口功率约束:
6)节点电压约束:
7)支路有功功率约束:
8)净负荷波动率约束:
所述步骤三包括:
1)建立第一辅助模型,求解决策变量的上、下限:
第一辅助模型:
式(29)、(30)中:Pt,n为第n个决策变量在t时刻的取值,h(P)=[hl,hu]为等式约束条件, g(P)≤[gl,gu]为不等式约束条件。
2)建立第二辅助模型,求解目标函数最优值的区间解的上、下界值:
式(31)、(32)中:W(Pt,n)为决策变量Pt,n的目标函数,h(P)=[hl,hu]为等式约束条件,g(P)≤[gl,gu]为不等式约束条件。
其中,将第二辅助模型,求解的结果分别作为目标函数最优值的下界和上界值,并由相对应的各决策变量Pt,n的取值作为各自区间解的边界值,得出在目标函数取最优值时各可调度机组的出力区间。
所述步骤五中,提出一种改进的仿射算法对配电网进行潮流计算,包括如下步骤:
S5.1:将配电网所有节点进行编号,输入ADN各节点参数和各支路参数;
S5.3:从末节点开始计算每节点的注入电流:
其中,式(34)具体的计算方法如下:
式(36)中:x0、z0为仿射中心值,xi、zj为仿射系数值,εn+1为新增噪声元。
S5.3.3:根据KCL计算上游节点电流,具体计算公式如下:
式(37)中:Ni为i节点及下游节点的集合,j为上游节点。
S5.4:前推计算每个节点的电压:
S5.6:计算配电网有功网损:
S5.7:输出区间潮流结果。
所述步骤六中,计算评价指标包括:
①、电压稳定性指标FVSI
式(40)中:FVSIij为支路ij电压稳定性指标;Vi为节点i的电压;Vj为节点j的电压。
②、线路容量裕度指标FML
③、算法改进精度指标CV
图2是本发明所提供的改进的IEEE33节点系统拓扑图,以此为例,本发明以该改进的IEEE33节点系统为例进行仿真分析。该主动配电网包含3个500kW的风力发电和2个600kW的光伏,配置200kW/800kWh的储能系统;3个单独储能系统为200kW/800kWh。结点1至5最大允许负载电流为500A,其余为300A。主动配电网、储能系统和分布式可再生能源的最大爬坡率分别为1000kW/h、200kW/h和200kW/h,为0.3/0.9,可控负荷变化的区间为[100,250]kW。
通过实际的光伏和风机的历史数据,采用上述灰色马尔科夫预测模型分别对光伏1、光伏2、风机1、风机2、风机3的功率区间进行预测,预测得到的区间结果如图3所示;如图可以看出光伏出力区间受时间影响较大,在00:00时刻减小为零,正午12:00时刻增大达到最大值;而风电出力区间受时间影响较小;通过将灰色马尔科夫预测得到的出力区间值改写成仿射形式即可带入优化调度模型中进行计算。
通过采用上述区间优化调度模型,对该主动配电网中分布式发电及储能系统进行区间优化调度,为了更加直观的展示,将区间调度结果分成区间中点和区间半径,区间调度结果的区间中点值如图4(1)所示;区间调度结果的区间半径值如图4(2)所示。如图可以看到,在1-5时段和22-24时段:ESS4、ESS5和ESS6通过风机向ADN供电;又因为是低电价时段,ESS1、ESS2和ESS3的充电,但在此时段内购电区间相对较小,表明配电网的波动较小。在10-16时段:光伏出力处于高峰,并通过ESS7、ESS8向ADN输入实现消纳,为了满足功率平衡,ADN需向主网大量购电;尤其是在14时段的购电区间相对较大但购电量最低,因为此时光伏输出功率和不确定性较大,导致ESS7与ESS8注入的不确定量大增,使得配电网购电电量的区间相应增大;ESS4在该时段的输出功率区间基本为0,输出较为稳定。根据区间调度结果得到运行费用和负荷波动区间分别为[73280.04,78051.7]元和[21710,33490]kW。
以ESS1为例,呈现其经调度后的荷电状态区间轨迹和边界曲线如图5所示;由图5可见,在13-17时段和21-24时段:ESS1的充放电功率区间较大,说明其受到ADN不确定性的影响的较大。7-9时段与23时段以后有达到SOC最小值的可能性,存在过放电情形,结合图4(1)和图4(2)可知,此时段内ESS1的调度区间基本为0,此时ESS1不会继续放电,因此该时段ESS1的SOC区间大小保持不变;同理在15-17时段及22时段以后存在充电过饱和情形,根据图4(1)和图4(2)可知,ESS1的调度区间同样为0,因此该时段的区间大小不变。
区间优化调度后,对主动配电网进行仿射潮流计算,得到系统的节点电压区间如图6 所示;可见:10-18和26-33等末节点的电压区间较宽,由于这些节点附近接入的DRESs存在不确定性,即节点注入功率的区间较宽。特别在0-8和20-24时段,28-33节点电压区间波动较大,存在电压过低的情况,需要引起相关人员注意。
通过计算系统电压稳定性指标及线路容量裕度指标,从而判断系统稳定性,其中电压稳定性指标计算结果如图7所示;图7直观体现了ADN的电压稳定性,可见:28-33节点的电压稳定性波动较大,若FVSI警戒值降低,则有电压崩溃的可能,其他节点的FVSI 区间较窄,电压崩溃的风险较小,而1和20节点附近的电压最稳定。
线路容量裕度指标如图8所示;其中,线路1-4的容量裕度区间有小于0的部分,存在线路过载阻塞即灵活性不足的情况,并不能很好地适应负荷的波动。其他线路的容量裕度区间均大于0,说明线路的灵活性较好,无过载现象发生,这些线路能够承受较大的负荷波动。
通过对比改进后的仿射潮流算法以及传统潮流算法,得到算法改进精度示意图如图9 所示;可以看出,所有CV值均大于0,说明与AA相比,IAA计算出的区间解更窄,能够有效的减小区间结果的保守性,使得区间结果更具有参考意义,特别是在节点20附近区间的保守性减小得最为明显。其中,算法效率对比如表1所示;
表1是本发明提供的计算效率对比表
IAA | AA | |
平均计算时间(s) | 15.5204 | 17.0458 |
由表1可以看出,IAA计算效率较AA更好,因为噪声元越多计算效率越差,若给定的区间仿射形式所含噪声元个数增多,那么IAA计算效率也越为明显。选取其中20:00时节点电压计算结果对比如图10所示;可见,IAA的结果均包含在AA内部,并且越靠近线路末端,其区间差异越大,IAA的保守性更小。
本发明采用改进后的灰色马尔科夫模型能够较为精准的预测出各类可再生能源的出力区间;通过区间泰勒公式改进后的仿射除法形式较为简单,由于固定了噪声元个数,缩短了计算时间,增大了计算效率;根据算法改进精度可知,通过仿射计算出的区间结果包含了改进仿射的区间结果,即改进仿射所受的误差影响较小,区间结果的保守性更小。
Claims (6)
1.基于改进仿射算法的分布式电源主动配电网优化调度方法,其特征在于包括以下步骤:
步骤一:基于灰色马尔科夫模型,得到风机和光伏的区间预测值;
步骤二:构建以运行费用最小、负荷波动最小的多目标优化模型,并采用TOPSIS法将多目标转换为单目标;
步骤三:建立辅助模型,求解各决策变量及目标函数最优值的区间解的上、下界值;
所述步骤三中:
建立第一辅助模型,求解决策变量的上、下限:
第一辅助模型:
Min Pt,n
s.t. h(P)=[hl,hu] (29)
g(P)≤[gl,gu]
t=1,2,…,T;n=1,2,…,N
Max Pt,n
s.t. h(P)=[hl,hu] (30)
g(P)≤[gl,gu]
t=1,2,…,T;n=1,2,…,N
式(29)、(30)中:Pt,n为第n个决策变量在t时刻的取值,h(P)=[hl,hu]为等式约束条件,g(P)≤[gl,gu]为不等式约束条件;
其中,将第一辅助模型求解的结果作为Pt,nl和Pt,nu,决策变量Pt,n满足约束条件的取值区间[Pt,nl,Pt,nu];
建立第二辅助模型,求解目标函数最优值的区间解的上、下界值:
Min MinW(Pt,n)
s.t.h(P)=[hl,hu](31)
g(P)≤[gl,gu]
Max MinW(Pt,n)
s.t.h(P)=[hl,hu](32)
g(P)≤[gl,gu]
式(31)、(32)中:W(Pt,n)为决策变量Pt,n的目标函数,h(P)=[hl,hu]为等式约束条件,g(P)≤[gl,gu]为不等式约束条件;
其中,将第二辅助模型求解的结果,分别作为目标函数最优值的下界和上界值,并由相对应的各决策变量Pt,n的取值作为各自区间解的边界值,得出在目标函数取最优值时各可调度机组的出力区间;
步骤四:将优化求解得到的区间调度结果、以及各功率区间值转换为仿射形式;
步骤五:基于改进仿射算法对配电网进行潮流计算,输出区间潮流结果;
所述步骤五包括如下步骤:
S5.1:将配电网所有节点进行编号,输入ADN各节点参数和各支路参数;
S5.3:从末节点开始计算每节点的注入电流:
其中,式(34)具体的计算方法如下:
式(36)中:x0、z0为仿射中心值,xi、zj为仿射系数值,εn+1为新增噪声元;
S5.3.3:根据KCL计算上游节点电流,具体计算公式如下:
式(37)中:Ni为i节点及下游节点的集合,j为上游节点;
S5.4:前推计算每个节点的电压:
S5.6:计算配电网有功网损:
S5.7:输出区间潮流结果;
步骤六:输出优化后的区间调度结果,并计算评价指标;
通过上述步骤完成主动配电网优化调度。
2.根据权利要求1所述基于改进仿射算法的分布式电源主动配电网优化调度方法,其特征在于:所述步骤一包括以下步骤:
S1.1:采用灰色GM(1,1)模型进行预测:
式(9)中:n为样本数量,a、b分别为发展系数和灰作用量,x(0)(1)为原始数据序列,x(0)(k+1)为灰色GM(1,1)模型预测值序列;
S1.2:引入误差修正模型,对预测值进行修正,修正公式如下所示:
式(10)中:X(0)为历史数据序列,X(2)为预测值序列,E(0)为误差修正序列,n为样本数量;
S1.3:对误差修正序列E(0)进行预测,得到E(1);再对预测得到的序列X(2)进行二次修正,得到预测修正后的序列XΔ(2);
S1.4:采用马尔科夫模型,得到如下各状态的转移概率、以及状态转移概率矩阵:
式(11)、(12)中:qij(k)为状态Ei经过k步转移到状态Ej的概率,Mi为处于状态Ei的样本数,Mij(k)为状态Ei经过k步转移到状态Ej的样本数,i、j=1,…,n;q(k)为状态转移概率矩阵;
S1.5:根据q(k)计算各步的自相关系数和权重系数,最终确定出历史数据对应预测值的每种状态的概率qi,并结合预测修正后的序列XΔ(2),得到最终预测值X:
式(13)中:XΔ(2)为预测修正后的序列,qi为每种状态的转移概率,Ei,mid为第i个状态区间的中点值,X为最终预测值。
3.根据权利要求1所述基于改进仿射算法的分布式电源主动配电网优化调度方法,其特征在于:所述步骤二包括以下步骤:
S2.1:以配电网运行总费用最小,建立的第一目标函数为:
式(14)中:T为总时段数,Ness、Nw、Npv和Nint分别为储能、风机、光伏和可控负荷的数量;和分别为t时段ADN的购电价格、购电量和有功网损;和分别为储能、风机和光伏的平准化度电成本;和分别为储能、风机、光伏和可控负荷在t时段的功率,且均为仿射值;
S2.2:以配电网负荷波动最小,建立的第二目标函数为:
式(15)中:Pt Load和Pav分别为ADN在t时段的总负荷和整个调度周期内负荷均值;
S2.3:采用TOPSIS法,将多目标转换为如下单目标:
4.根据权利要求3所述基于改进仿射算法的分布式电源主动配电网优化调度方法,其特征在于:所述步骤二还包括确定调度问题的基本约束条件:
1)、有功功率平衡约束:
式(17)中:和分别是t时段AND的购电量和有功网损;和分别为风机、光伏、储能和可控负荷在t时段的功率,且均为仿射值;为t时段节点i负荷的有功功率;Nw、Npv、Ness和Nint分别为风机、光伏、储能和可控负荷的数量;I为AND的节点总数;
2)、区间潮流等式约束:
3)、DG出力与爬坡率约束:
4)、储能系统约束:
5)、配变关口功率约束:
6)、节点电压约束:
7)、支路有功功率约束:
8)、净负荷波动率约束:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110859394.6A CN113572163B (zh) | 2021-07-28 | 2021-07-28 | 基于改进仿射算法的分布式电源主动配电网优化调度方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110859394.6A CN113572163B (zh) | 2021-07-28 | 2021-07-28 | 基于改进仿射算法的分布式电源主动配电网优化调度方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113572163A CN113572163A (zh) | 2021-10-29 |
CN113572163B true CN113572163B (zh) | 2023-04-07 |
Family
ID=78168659
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110859394.6A Active CN113572163B (zh) | 2021-07-28 | 2021-07-28 | 基于改进仿射算法的分布式电源主动配电网优化调度方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113572163B (zh) |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114221344A (zh) * | 2021-11-05 | 2022-03-22 | 三峡大学 | 基于改进仿射算法的前推回代潮流计算方法 |
CN114204548A (zh) * | 2021-11-19 | 2022-03-18 | 国网辽宁省电力有限公司鞍山供电公司 | 一种配电线路带电作业多维数据熵权法状态识别方法 |
CN114169627B (zh) * | 2021-12-14 | 2022-11-18 | 湖南工商大学 | 一种深度强化学习分布式光伏发电激励的系统 |
CN114336686B (zh) * | 2022-01-10 | 2023-11-07 | 武汉大学 | 一种基于自适应预测控制系统的电网区域振荡抑制方法 |
CN114254765A (zh) * | 2022-03-01 | 2022-03-29 | 之江实验室 | 一种面向模拟推演的主动式序列决策方法、装置和介质 |
CN114662756A (zh) * | 2022-03-23 | 2022-06-24 | 福州大学 | 考虑源-荷不确定性的热电微能网双层仿射优化调度方法 |
CN116522584B (zh) * | 2023-03-07 | 2023-10-27 | 北京智中能源科技发展有限公司 | 一种配电网供电爬坡能力最大化计算的优化方法 |
CN118117625B (zh) * | 2024-04-30 | 2024-08-23 | 宁德时代新能源科技股份有限公司 | 储能系统功率调度分配方法、装置、存储介质及能源管理系统 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108564206A (zh) * | 2018-03-27 | 2018-09-21 | 中国农业大学 | 一种基于分布式优化和空间相关性的风电功率预测方法 |
Family Cites Families (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2009016113A1 (en) * | 2007-07-27 | 2009-02-05 | Abb Research Ltd | Controlling a power system |
US8566266B2 (en) * | 2010-08-27 | 2013-10-22 | Mitsubishi Electric Research Laboratories, Inc. | Method for scheduling the operation of power generators using factored Markov decision process |
CN106786570B (zh) * | 2017-03-10 | 2019-05-10 | 国网山东省电力公司经济技术研究院 | 一种含风电电力系统区间经济调度方法 |
CN106849190B (zh) * | 2017-03-21 | 2019-05-17 | 国网天津市电力公司 | 一种基于Rollout算法的多能互补微网实时调度方法 |
CN110518569A (zh) * | 2019-06-12 | 2019-11-29 | 广西电网有限责任公司电力科学研究院 | 一种计及电网运行收益与风险的概率优化调度模型 |
CN112217202B (zh) * | 2020-09-29 | 2022-06-14 | 东南大学 | 考虑灵活性投资的分布式新能源、储能与配电网规划方法 |
CN112736926A (zh) * | 2020-12-25 | 2021-04-30 | 华北电力大学 | 一种分布式新能源接入配电网区间仿射潮流动态优化方法 |
-
2021
- 2021-07-28 CN CN202110859394.6A patent/CN113572163B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108564206A (zh) * | 2018-03-27 | 2018-09-21 | 中国农业大学 | 一种基于分布式优化和空间相关性的风电功率预测方法 |
Also Published As
Publication number | Publication date |
---|---|
CN113572163A (zh) | 2021-10-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113572163B (zh) | 基于改进仿射算法的分布式电源主动配电网优化调度方法 | |
Hemmati et al. | Market based transmission expansion and reactive power planning with consideration of wind and load uncertainties | |
CN107591844A (zh) | 考虑节点注入功率不确定性的主动配电网鲁棒重构方法 | |
CN109508499A (zh) | 多时段多场景分布式电源最优接入位置与容量研究方法 | |
CN107317361B (zh) | 一种考虑区域自治能力的主动配电网全局优化调度方法 | |
CN107968439A (zh) | 基于混合整数线性规划的主动配电网联合优化算法 | |
CN112736926A (zh) | 一种分布式新能源接入配电网区间仿射潮流动态优化方法 | |
CN107706932A (zh) | 一种基于动态自适应模糊逻辑控制器的能量调度优化方法 | |
CN115640963A (zh) | 一种考虑投资运营模式的海上风电接入系统鲁棒规划方法 | |
CN103715688A (zh) | 并网型微网优化配置方法及装置 | |
CN107834540B (zh) | 一种基于概率约束的分布式光伏接入容量的确定方法 | |
CN115423282A (zh) | 一种基于奖惩阶梯型碳交易的电-氢-储综合能源网多目标优化调度模型 | |
Xiao et al. | A multi‐energy complementary coordinated dispatch method for integrated system of wind‐photovoltaic‐hydro‐thermal‐energy storage | |
Khoubseresht et al. | An analytical optimum method for simultaneous integration of PV, wind turbine and BESS to maximize technical benefits | |
CN114884136A (zh) | 一种计及风电相关性的主动配电网鲁棒优化调度方法 | |
CN105633950B (zh) | 一种考虑风电注入不确定性的多目标随机模糊动态最优潮流求解方法 | |
CN109840621A (zh) | 考虑储能系统影响的并网型微电网日前多目标优化方法 | |
An et al. | Real-time optimal operation control of micro energy grid coupling with electricity-thermal-gas considering prosumer characteristics | |
CN116995740A (zh) | 一种配电网分布式风电及储能优化配置方法和系统 | |
CN112260274A (zh) | 一种基于全景理论的虚拟电厂组建方法 | |
Luo et al. | Two-stage optimal scheduling of virtual power plant considering demand response and forecast errors | |
CN114759616B (zh) | 一种考虑电力电子器件特性的微电网鲁棒优化调度方法 | |
Hamdi et al. | Optimum configuration of a dispatchable hybrid renewable energy plant using artificial neural networks: Case study of Ras Ghareb, Egypt. | |
CN111404162B (zh) | 计及电压稳定约束的区域电网调度域划分方法 | |
Krim et al. | Joint optimisation of sizing and fuzzy logic power management of a hybrid storage system considering economic reliability indices |
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 |