CN117874562A - 一种基于改进混沌优化算法的分布式能源系统双层运行优化方法 - Google Patents
一种基于改进混沌优化算法的分布式能源系统双层运行优化方法 Download PDFInfo
- Publication number
- CN117874562A CN117874562A CN202410059667.2A CN202410059667A CN117874562A CN 117874562 A CN117874562 A CN 117874562A CN 202410059667 A CN202410059667 A CN 202410059667A CN 117874562 A CN117874562 A CN 117874562A
- Authority
- CN
- China
- Prior art keywords
- optimization model
- layer optimization
- algorithm
- follows
- optimal
- 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
- 238000005457 optimization Methods 0.000 title claims abstract description 225
- 238000000034 method Methods 0.000 title claims description 48
- 238000004146 energy storage Methods 0.000 claims abstract description 27
- 230000007613 environmental effect Effects 0.000 claims abstract description 19
- 238000012423 maintenance Methods 0.000 claims abstract description 9
- 238000010438 heat treatment Methods 0.000 claims abstract description 8
- 238000001816 cooling Methods 0.000 claims abstract description 7
- 239000007789 gas Substances 0.000 claims description 25
- 239000011159 matrix material Substances 0.000 claims description 24
- 230000000739 chaotic effect Effects 0.000 claims description 21
- 230000008569 process Effects 0.000 claims description 14
- 238000013486 operation strategy Methods 0.000 claims description 12
- 230000006870 function Effects 0.000 claims description 10
- 230000003247 decreasing effect Effects 0.000 claims description 9
- 230000005611 electricity Effects 0.000 claims description 8
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 claims description 8
- 230000009467 reduction Effects 0.000 claims description 7
- 230000007423 decrease Effects 0.000 claims description 6
- 230000003993 interaction Effects 0.000 claims description 6
- 238000012163 sequencing technique Methods 0.000 claims description 6
- 239000013598 vector Substances 0.000 claims description 6
- 238000010521 absorption reaction Methods 0.000 claims description 5
- 230000008859 change Effects 0.000 claims description 5
- 238000013461 design Methods 0.000 claims description 5
- 239000000446 fuel Substances 0.000 claims description 5
- 238000013178 mathematical model Methods 0.000 claims description 5
- 238000004364 calculation method Methods 0.000 claims description 4
- 239000003345 natural gas Substances 0.000 claims description 4
- 238000010606 normalization Methods 0.000 claims description 4
- 230000002776 aggregation Effects 0.000 claims description 3
- 238000004220 aggregation Methods 0.000 claims description 3
- 238000009826 distribution Methods 0.000 claims description 2
- 230000007246 mechanism Effects 0.000 claims description 2
- 238000012946 outsourcing Methods 0.000 claims description 2
- 238000005096 rolling process Methods 0.000 claims description 2
- 239000010410 layer Substances 0.000 abstract description 139
- 239000002356 single layer Substances 0.000 abstract description 18
- 230000006872 improvement Effects 0.000 abstract description 7
- CURLTUGMZLYLDI-UHFFFAOYSA-N Carbon dioxide Chemical compound O=C=O CURLTUGMZLYLDI-UHFFFAOYSA-N 0.000 abstract description 6
- 230000008901 benefit Effects 0.000 abstract description 4
- 229910002092 carbon dioxide Inorganic materials 0.000 abstract description 3
- 239000001569 carbon dioxide Substances 0.000 abstract description 2
- 238000011160 research Methods 0.000 abstract description 2
- 230000008878 coupling Effects 0.000 abstract 1
- 238000010168 coupling process Methods 0.000 abstract 1
- 238000005859 coupling reaction Methods 0.000 abstract 1
- 238000005265 energy consumption Methods 0.000 abstract 1
- 238000011161 development Methods 0.000 description 5
- 238000010586 diagram Methods 0.000 description 5
- 238000005057 refrigeration Methods 0.000 description 5
- 238000011156 evaluation Methods 0.000 description 4
- 238000005338 heat storage Methods 0.000 description 3
- 238000009434 installation Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 2
- 230000008092 positive effect Effects 0.000 description 2
- 229910052799 carbon Inorganic materials 0.000 description 1
- 230000000052 comparative effect Effects 0.000 description 1
- 238000003912 environmental pollution Methods 0.000 description 1
- 239000002803 fossil fuel Substances 0.000 description 1
- 238000003064 k means clustering Methods 0.000 description 1
- 230000035800 maturation Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 238000010248 power generation Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
- 239000002699 waste material Substances 0.000 description 1
Classifications
-
- 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
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明建立一个双层优化框架,以获得太阳能和热能储存耦合冷热电联供系统的最佳配置。外层以组件尺寸为决策变量,优化系统的年经济效益、年二氧化碳排放量和年一次能源消耗量;内层以设备输出为决策变量,优化运行和维护成本。以一家大型酒店为案例进行研究。优化结果表明,所提出的双层优化模型可以实现混合CCHP的最优配置方案和调度策略。与基于规则的多目标单层优化模型相比,所提出的双层优化模型在经济、能源和环境性能方面的最大改进率分别为30.63%、37.94%和38.16%。与基于规则的单目标单层优化模型相比,所提出的双层优化模型在经济、能源和环境性能方面的最大改进幅度分别为34.95%、35.81%和37.55%。此外,热能储存(TES)可大幅提高双层优化模型的最优性能。
Description
本发明属于分布式能源系统优化运行技术领域,尤其是一种基于改进混沌优化算法的分布式能源系统分层运行优化方法。
背景技术
大规模使用传统化石燃料带来的问题亟待紧急解决。充分利用可再生和分布式能源技术已成为解决能源危机和环境污染的有效途径。作为典型的分布式能源系统,冷热电三联供(CCHP)系统因其高能效、节能、环保和脱网运行等特性,已成为各国积极推广的技术热点。结合可再生能源建设CCHP系统有助于发挥多种形式可再生能源的优势,总体上实现丰富、清洁、可持续的发展。此外,CCHP系统的灵活能源供应模式能够很好地适应可再生能源的随机、波动和间歇性特性。建立清洁、低碳、安全、高效的能源供应模式,实现多能互补和能梯度利用,是未来能源系统发展的重要方向。随着新能源发电技术的迅猛发展和CCHP系统的成熟,综合能源系统研究已成为电力乃至整个能源领域改革的突破口,是未来能源领域的主要发展趋势。
发明内容
本发明针对考虑热能存储的混合冷热电联供系统,构建一个多目标双层CCHP优化模型,将k均值、非支配排序和相对决策度排序与两层优化模型结合,形成一个完整的两层CCHP优化模型,实现最优配置方案和运行策略;
本发明解决其技术问题是采取以下技术方案实现的:
一种基于改进混沌优化算法的分布式能源系统双层运行优化方法,包括以下步骤:
步骤1,建立包含外层优化和内层优化的双层模型;
步骤2,求解双层优化模型,将参考建筑的年负荷数据和历史气象数据输入矩阵;通过k-means算法对矩阵中的数据进行聚类,从而得到聚类后的情景及其相应的情景概率。将典型情景及其相应的情景概率输入外层优化模型;
步骤3,外层优化模型随机初始化一个设备容量场景群,内层优化模型读取初始化的设备容量群。内层优化模型采用混合整数线性规划(MILP)算法,以系统运行成本最小为目标,获得最佳运行策略,并将最佳运行方案传递给外层优化模型;
步骤4,外层优化模型以内层优化模型得到的运行方案为基础,利用多目标混沌优化(MOCGO)算法进行位置更新,得到新的设备装机容量群,并将新的装机容量群返回给内层优化模型;
步骤5,内层优化模型利用新的装机容量群再次求解最优运行方案,如此循环迭代。当达到最大迭代次数时,循环将停止。然后使用非支配排序和相对决策度排序对解决方案集进行排序,选择非支配解集中相对决策度最大的作为最优个体;
步骤1,构建的冷热电三联供系统的主要设备包括电制冷机(EC)、吸收式制冷机(AC)、热交换装置(HE)、热能储存罐(TES)、燃气锅炉(GB)、燃气轮机(GT)、太阳能集热器(ST)和光伏装置(PV)。双层优化模型包括外层优化和内层优化两部分,外层优化模型如下:
(1)决策变量:
X=[NEC,NAC,NTES,NGB,APV,AST,NMT] (1)
(2)在外层优化模型中,经济、能源和环境目标函数如下:
式中,fEco、fEne和fEnv分别为经济、能源和环境目标函数。IV是利率,lifem是设备寿命,Cm是第m个设备的投资成本,Nm是第m个设备的装机容量。fGas是燃料成本,fGrid是电网互动成本,fMain是设备维护成本。它们的定义如下:
式中,S为情景总数,Ps为情景概率,Cf为天然气价格,FGB和FGT分别为GB和GT的天然气消耗量,EOUT为售电量,EIN为外购电量,Ce为使用时间,CEC为EC的制冷量,QGB为GB的制热量,EGT为GT的发电量,ε为设备维护成本系数;
式中,FGrid是电网的燃料消耗量,λf和λe分别是燃气和电网的等效CO2排放系数。λf取值为0.968kg/kWh,λe取值为0.220kg/kWh;
(3)设备容量的约束条件如表1所示:
表1设计参数范围
(4)内层优化模型的决策变量主要是各设备的输出功率,其定义如下:
X=[EMT(1:24·S),QGB(1:24·S),…,CEC(1:24·S)] (8)
(5)内层优化模型的目标是最小化系统的年度运行和维护成本,其定义如下:
fInner=fMain+fGrid+fGas (9)
(6)约束条件如下:
ELoad+EEC+EGrid,in=EMT+EPV+EGrid,out (10)
QLoad+QAC+QTES,in=QHR+QTES,out+QGB (11)
CLoad=CAC+CEC (12)
步骤2,求解双层优化模型,将参考建筑的年负荷数据和历史气象数据输入矩阵,通过k-means算法对矩阵中的数据进行聚类,从而得到聚类后的情景及其相应的情景概率。将典型情景及其相应的情景概率输入外层优化模型;
(1)将一年的负荷和气象数据存储到一个矩阵中,并对该矩阵进行归一化。归一化的公式如下:
式中,Ai是归一化后的元素;Xi是原始数据在归一化前的元素;Xmax和Xmin分别是归一化前数据的最大值和最小值;
(2)使用k均值对归一化后的矩阵进行聚类。我们获取典型场景的数量k以及每个簇的中心日。与中心日对应的原始数据是该簇的典型场景,而每个簇中包含的天数即为这些典型场景的天数。每个场景的概率如下:
式中,Pk是每个场景的概率,Nk表示每个场景中包含的天数;
在使用k均值算法对数据进行聚类之前,需要指定一个确定的簇数。然而,如果聚类值的选择基于用户的主观意愿,将不可避免地使最终的聚类结果更加随机和不确定。而且,即使对于相同的样本集,当不同用户使用不同数量的簇来处理数据时,最终的聚类结果可能会有很大的差异。因此,采用适当的方法确定簇的数量并确保最终聚类结果的客观性是非常重要的。作为确定最优簇数的客观方法,肘部法则的核心思想包括以下几点,随着样本分类数的增加,每个簇中样本的聚集程度将逐渐增加,平方误差和(SSE)自然会逐渐减小。当簇的数量k小于最优簇数时,随着簇的数量k的增加,SSE的减小趋势会逐渐减缓。当簇的数量k等于最优簇数时,SSE的减小趋势会突然减小。此时,继续增加k的值,SSE将不会有明显的变化。因此,使得SSE减小幅度突然减小的簇的数量k,就是数据的实际簇数。肘部法则的核心度量指标是平方误差和(SSE),具体定义如下:
式中,k表示簇的数量;Ci表示第i个簇;p表示第i个簇的样本点;mi表示第i个簇的质心;
步骤3,外层优化模型随机初始化一个设备容量场景群,内层优化模型读取初始化的设备容量群。内层优化模型采用MILP算法,以系统运行成本最小为目标,获得最佳运行策略,并将最佳运行方案传递给外层优化模型;
步骤4,外层优化模型以内层优化模型得到的运行方案为基础,利用MOCGO进行位置更新,得到新的设备装机容量群,并将新的装机容量群返回给内层优化模型。所选多目标智能算法为多目标混沌优化算法,混沌博弈优化算法(CGO)是一种基于混沌理论原理提出的优化算法,它利用分形和混沌博弈的基本概念。其初始化公式如下:
(1)在该算法中,每个候选解Xi包含一些决策变量Xi,j,这些变量表示西尔平斯基三角形中这些合格种子的位置。在这个算法中,西尔平斯基三角形被视为候选解的搜索空间。其数学模型如下:
式中,n是西尔平斯基三角形中合格点的数量,d是这些点的维数。考虑的合格点的初始位置是在搜索空间中随机确定的,具体过程如下:
式中,是合格点的初始位置;/>和/>分别是第i个候选解的第j个决策变量的最小和最大允许值;rand.()是一个在区间[0,1]内的随机数;
(2)第一个种子位置更新的数学表示如下:
式中,βi和γi分别表示1或2的随机整数,用于模拟掷骰子的可能性。GB是当前的最优解,MGi是初始合格点的平均值,它们被视为第一个临时三角形的三个顶点。αi是一个随机生成的矩阵,Xi是第i个候选解;
(3)类似于生成第一个种子的过程,第二个种子可以沿着MGi和Xi之间的直线移动,并利用一些随机生成的因素来约束该移动。该过程的数学表示如下:
(3)第三个种子是基于MGi的位置随机生成的,其数学模型如下:
(4)第四个种子的位置更新是基于对随机选择的决策变量进行一些随机更改。第四个种子的数学建模如下:
式中,R是在[0,1]上均匀分布的随机数;k是一个在[1,d]范围内的随机整数向量。此外,为控制种子的移动限制,对于α提出四个不同的公式,可以表示如下:
式中,rand是在[0,1]范围内的随机向量;δ和ε是在[0,1]范围内的随机向量;
(5)对于多目标混沌优化算法,进行初始适应度。对于每个个体Xj,其强度值的定义如下:
S(Xj)=|{Xw|Xw∈P∪A∧Xj<Xw}| (23)
式中,Xj和Xw既是种群中的个体,也是存档个体。对于个体Xi来说,新的适应度值由两部分组成,密度信息D(Xi)和原始适应度值R(Xi)。R(Xi)的定义如下:
式中,R(Xi)个体在群体和外部档案中的主导地位决定Xi个体在群体和外部档案中的主导地位;
(6)关于基于移位的密度估计,如果某个个体在某个目标上的表现优于p,那么它就会被转移到与p在该目标上相同的位置;否则,它就会保持不变。个体p在种群D′(p,P)中的新密度可以表示为:
D′(p,P)=D(dist(p,q′1),dist(p,q′2),…dist(p,q′N-1)) (25)
式中,N为种群大小;dist(p,q′1)为两个个体相似度的程度;q′1为个体qi(qi≠pand qi∈P)的位移版本,其定义如下:
式中,m代表目标数,P(j)、qi(j)和分别代表个体p、q和q′的第j个目标值;
步骤5,内层优化模型利用新的装机容量群再次求解最优运行方案,如此循环迭代。当达到最大迭代次数时,循环将停止;
对于最优个体的选取,采用非支配排序和相对决策度排序获得帕累托前沿和最优个体,非支配排序可以筛选出不同个体的适应度优劣级别,相对决策度可以表达解集的分布情况。下面介绍基于相对决策度的非支配排序机制:
(1)对解集进行非支配排序,按照不同支配等级对解集进行划分。对于种群每次进化产生的帕累托前沿,采用外部存档的方法进行存储与更新;
(2)对于外部档案中划分好支配等级的解集,依据支配等级依次计算不同解的相对决策度并按降序排列,选择非支配解集中相对决策度最大的作为最优个体。相对决策度的计算公式如下:
式中,为第i个解在搜索空间中的决策度,Nq表示同一支配等级的解的数量,和/>为第i个解在第j个维度最接近的两个维度值,且满足/>和/>为同支配等级的解中第j个维度的最大值和最小值。RDMPi为第i个解在目标空间中的相对决策度,Nobj表示目标函数的数量;/>和/>为同支配等级的解中第j个目标函数的最大值和最小值;
上述的一种基于改进混沌优化算法的混合能源系统双层运行优化方法,所述的多目标优化算法,是为本技术领域技术人员所熟知的;
本发明的优点和积极效果是:
1、本发明构建一个多目标两层CCHP规划模型。我们将k均值和非支配排序和相对决策度排序方法与两层优化模型结合,形成一个完整的两层CCHP优化模型。提出一种新颖的多目标混沌优化算法(MOCGO)。通过对混沌博弈优化算法的改进,我们得到MOCGO;
2、本发明提出解决两层优化模型的新方法。我们将MOCGO与MILP结合起来解决构建的两层优化模型,并通过与其他算法的比较验证所提算法具有更卓越的优化结果。分析并比较在考虑和不考虑TES的情况下,两层优化模型下CCHP系统的设计性能;
3、本发明所提出的两层优化模型能够同时实现最优配置方案和运行策略。我们将所提出的两层优化模型与单层优化模型进行比较,并证明所提出的模型能够同时实现最优配置方案和调度策略;
附图说明
下面结合附图和实施例对本发明进一步说明;
图1是一种基于改进混沌优化算法的分布式能源系统双层运行优化方法总流程图;
图2是本实施例最终的聚类数图;
图3是本实施例具体的聚类结果图;
图4是改进混沌优化算法的工作流程图;
图5是本实施例年度逐日小时负荷负荷图;
图6是本实施例环境温度和辐照度图;
图7是本实施例分时电价和燃气价格图;
图8是本实施例混合CCHP系统的电能平衡图;
图9是本实施例混合CCHP系统的热能平衡图;
图10是本实施例混合CCHP系统的制冷能量平衡图;
图11是本实施例不考虑热能储存的情况下两层优化模型的最优性能因子图;
图1表明一种基于改进混沌优化算法的分布式能源系统双层运行优化方法的流程是:开始→建立包含外层优化和内层优化的双层模型→求解双层优化模型,将参考建筑的年负荷数据和历史气象数据输入矩阵;通过k-means算法对矩阵中的数据进行聚类,从而得到聚类后的情景及其相应的情景概率;将典型情景及其相应的情景概率输入外层优化模型→外层优化模型随机初始化一个设备容量场景群,内层优化模型读取初始化的设备容量群。内层优化模型采用MILP算法,以系统运行成本最小为目标,获得最佳运行策略,并将最佳运行方案传递给外层优化模型外层优化→模型以内层优化模型得到的运行方案为基础,利用MOCGO进行位置更新,得到新的设备装机容量群,并将新的装机容量群返回给内层优化模型→内层优化模型利用新的装机容量群再次求解最优运行方案,如此循环迭代。当达到最大迭代次数时,循环将停止。然后使用非支配排序和相对决策度排序对解决方案集进行排序,选择非支配解集中相对决策度最大的作为最优个体→结束;
实施例
本发明是一种基于改进混沌优化算法的分布式能源系统双层运行优化方法,双层运行优化方法的流程如图1所示;
步骤1,建立包含外层优化和内层优化的双层模型;
步骤2,求解双层优化模型,将参考建筑的年负荷数据和历史气象数据输入矩阵;通过k-means算法对矩阵中的数据进行聚类,从而得到聚类后的情景及其相应的情景概率;将典型情景及其相应的情景概率输入外层优化模型;
步骤3,外层优化模型随机初始化一个设备容量场景群,内层优化模型读取初始化的设备容量群。内层优化模型采用MILP算法,以系统运行成本最小为目标,获得最佳运行策略,并将最佳运行方案传递给外层优化模型;
步骤4,外层优化模型以内层优化模型得到的运行方案为基础,利用MOCGO进行位置更新,得到新的设备装机容量群,并将新的装机容量群返回给内层优化模型;
步骤5,内层优化模型利用新的装机容量群再次求解最优运行方案,如此循环迭代。当达到最大迭代次数时,循环将停止。然后使用非支配排序和相对决策度排序对解决方案集进行排序,选择非支配解集中相对决策度最大的作为最优个体;
步骤1,构建的冷热电三联供系统的主要设备包括电制冷机(EC)、吸收式制冷机(AC)、热交换装置(HE)、热能储存罐(TES)、燃气锅炉(GB)、燃气轮机(GT)、太阳能集热器(ST)和光伏装置(PV)。双层优化模型包括外层优化和内层优化两部分,对于外层优化模型,用式(1)建立决策变量,式(2)建立经济、能源和环境目标函数,用式(3)建立燃料成本,式(4)建立电网互动成本,式(5)、式(6)和式(7)建立设备维护成本,内层优化模型的决策变量主要是各设备的输出功率,用式(8)建立决策变量,内层优化模型的目标是最小化系统的年度运行和维护成本,用式(9)建立目标函数,用式(10)、(11)和(12)建立约束条件;
步骤2,求解双层优化模型,将参考建筑的年负荷数据和历史气象数据输入矩阵;通过k-means算法对矩阵中的数据进行聚类,从而得到聚类后的情景及其相应的情景概率;将典型情景及其相应的情景概率输入外层优化模型。将一年的负荷和气象数据存储到一个矩阵中,并对该矩阵进行归一化,用式(13)建立归一化。使用k均值对归一化后的矩阵进行聚类。我们获取典型场景的数量k以及每个簇的中心日。与中心日对应的原始数据是该簇的典型场景,而每个簇中包含的天数即为这些典型场景的天数,用式(14)建立每个场景的概率,在使用k均值算法对数据进行聚类之前,需要指定一个确定的簇数。然而,如果聚类值的选择基于用户的主观意愿,将不可避免地使最终的聚类结果更加随机和不确定。而且,即使对于相同的样本集,当不同用户使用不同数量的簇来处理数据时,最终的聚类结果可能会有很大的差异。因此,采用适当的方法确定簇的数量并确保最终聚类结果的客观性是非常重要的。作为确定最优簇数的客观方法,肘部法则的核心思想包括以下几点,随着样本分类数的增加,每个簇中样本的聚集程度将逐渐增加,平方误差和(SSE)自然会逐渐减小。当簇的数量k小于最优簇数时,随着簇的数量k的增加,SSE的减小趋势会逐渐减缓。当簇的数量k等于最优簇数时,SSE的减小趋势会突然减小。此时,继续增加k的值,SSE将不会有明显的变化。因此,使得SSE减小幅度突然减小的簇的数量k,就是数据的实际簇数。肘部法则的核心度量指标是平方误差和(SSE),用式(15)建立平方误差和,最终的聚类数如图2所示,使用K均值聚类算法对参考建筑的典型年度负荷和气象数据进行聚类,具体的聚类结果如图3所示;
步骤3,外层优化模型随机初始化一个设备容量场景群,内层优化模型读取初始化的设备容量群。内层优化模型采用MILP算法,以系统运行成本最小为目标,获得最佳运行策略,并将最佳运行方案传递给外层优化模型;
步骤4,外层优化模型以内层优化模型得到的运行方案为基础,利用MOCGO进行位置更新,得到新的设备装机容量群,并将新的装机容量群返回给内层优化模型。所选多目标智能算法为多目标混沌优化算法,混沌博弈优化算法(CGO)是一种基于混沌理论原理提出的优化算法,它利用分形和混沌博弈的基本概念,改进混沌优化算法的工作流程图如图4所示,在该算法中,每个候选解Xi包含一些决策变量Xi,j,这些变量表示西尔平斯基三角形中这些合格种子的位置。在这个算法中,西尔平斯基三角形被视为候选解的搜索空间,用式(16)建立数学模型,用式(17)确定初始位置,用式(18)建立第一个种子位置更新,类似于生成第一个种子的过程,第二个种子可以沿着MGi和Xi之间的直线移动,并利用一些随机生成的因素来约束该移动,第三个种子是基于MGi的位置随机生成的,用式(19)建立第二个种子位置更新,用式(20)建立第三个种子位置更新,第四个种子的位置更新是基于对随机选择的决策变量进行一些随机更改,用式(21)建立第四个种子的位置更新,此外,为控制种子的移动限制,用式(22)建立约束条件。对于多目标混沌优化算法,进行初始适应度,用式(23)定义强度值,用式(24)定义原始适应度值,关于基于移位的密度估计,如果某个个体在某个目标上的表现优于p,那么它就会被转移到与p在该目标上相同的位置;否则,它就会保持不变,用式(25)表示个体在种群中的新密度,用式(26)表示位移版本;
步骤5,内层优化模型利用新的装机容量群再次求解最优运行方案,如此循环迭代。当达到最大迭代次数时,循环将停止。然后使用非支配排序和相对决策度排序对解决方案集进行排序,选择非支配解集中相对决策度最大的作为最优个体,式(27)和(28)建立相对决策度的具体计算;
以美国能源部提供的一家大型酒店的能源数据为案例研究,本发明采用提出的优化方法解决混合CCHP两层优化模型,以实现系统的最优综合性能。图5展示年度逐日小时负荷,图6显示环境温度和辐照度,图7呈现分时电价和燃气价格;
图8展示混合CCHP系统的电能平衡。可以看到,酒店的电力消耗在整个年度内非常稳定,而燃气轮机几乎负责酒店的全部电负荷和电吸收制冷机的功耗。在构建的两层优化模型下,混合CCHP几乎可以独立于电网为客户供电,而不会产生多余的电力。此外,还发现在第五和第十二个典型场景中,电吸收制冷机的功耗最小。这种情况发生的原因是,在上述两个场景中,环境温度和辐射使得太阳能热产生大量的热能,因此空气制冷机几乎完全满足用户的制冷需求;
图9详细描述混合CCHP系统的热能平衡。可以观察到,混合CCHP系统在供热方面表现良好,避免能源浪费的问题。太阳能热在不同的典型场景中热能供应呈现较大变化,而燃气轮机的热能供应则非常稳定。考虑到设备的多次启停可能影响其使用寿命,本发明提出的两层优化模型对于延长燃气轮机和热回收设备的寿命具有积极作用。当环境条件不适宜太阳能热运行时,燃气锅炉会接管部分用户的热需求。此外,蓄热系统的储能工作时间、空调制冷工作时间与太阳能热的供热之间存在紧密联系,特别值得注意的是,在太阳能热产生的热能较大时,蓄热系统将储存部分热能,而空调将利用这部分热能进行制冷;而在太阳能热产生的热能较小时,蓄热系统则释放部分储存的热能,此时空调可能停止运行;
图10展示混合CCHP系统的制冷能量平衡。可以看到,混合CCHP系统可以可靠地满足客户的制冷负荷。此外,电吸收制冷机在满足客户制冷负荷方面起着更为关键的作用。另一方面,当太阳能热供热较多时,空调将几乎承担客户的全部制冷负荷,并在环境条件不适合太阳能热工作时,空调可能只会承担一小部分或甚至停机;
本部分设置4个案例用于验证所提方法的有效性;
案例1:采用本发明所提方法进行优化;
案例2:采用多目标单层优化模型进行优化;
案例3:采用单目标单层优化模型进行优化;
案例4:采用未考虑热能储存的双层优化模型进行优化;
在案例1中,为验证MOCGO在多能源系统规划运算中的适用性和有效性,将比较dMOPSO-MILP、MOEAD-MILP、MOPSOCD-MILP和MOCGO-MILP分别获得的最佳配置解在经济、能源和环境性能方面的不同表现;
表2不同算法得出的双层优化模型结果
表2展示解决两层优化模型的不同优化算法的最优结果。可以看出,所提出的MOCGO在经济成本、燃气消耗和环境排放方面均优于比较算法。这些对比结果说明,MOCGO在处理多研究的两层优化模型时更为高效,能够为用户提供更优的解决方案;
为证明本发明提出的两层优化模型能够同时实现混合CCHP的最优配置和运行策略,将该两层优化模型与基于规则的单层优化模型进行比较。对于单层优化模型,最大迭代次数为500。使用一次能源节约比(PESR)、年度成本节约率(ACSR)和二氧化碳排放减少率(CDERR)作为评价指标,比较单层优化模型和两层优化模型之间的性能差异;
在案例2中,为验证本研究提出的两层优化模型在提高系统经济、能源和环境性能方面的优势,将其与基于规则的单层多目标优化模型进行比较。基于规则的运行策略主要包括FTL、FEL、FTL-ECR、FEL-ECR和FHL;
表3多目标优化模型的评价指标值
表3呈现多目标优化模型的PESR、ACSR和CDERR数值。可以看到,在任何策略下,与单层优化模型相比,两层优化模型始终在提高系统的能源、经济和环境性能方面都是有效的。特别是在ACSR方面,两层优化模型实现最大经济性能的提升,达到30.63%。在PESR方面,两层优化模型实现高达37.94%的燃气消耗节约。而在CDERR方面,两层优化模型实现最大的CO2减排,为38.16%;
在光伏(PV)和太阳能热(ST)装机容量方面,发现在FEL、FTL、FEL-ECR、FTL-ECR和FHL策略下的多目标优化模型的最佳配置具有特定的PV和ST装机容量数值。然而,在两层优化模型中,其最佳PV装机容量为零,而最佳ST装机面积则为屋顶面积。上述情况可能是由于酒店全年消耗的热负荷较大且稳定。关于燃气轮机(GT)装机容量,发现两层和多目标优化模型之间只有很小的差异。至于上网电网交互限制,通过dMOPSO得到的FEL-ECR和FHL的最佳值为零,表明在这个时候混合型CCHP几乎可以实现脱网运行;
关于TES的装机容量,通过dMOPSO和MOPSOCD获得的FTL和FTL-ECR的最佳装机容量较小,因为系统几乎不产生多余的热能。就上网电网交互限制而言,通过dMOPSO得到的FEL-ECR和FHL的最佳值为0,表明在这个时候混合型CCHP几乎可以实现脱网运行。同样,两层优化模型也可以实现脱网运行。至于TES的装机容量,在FTL和FTL-ECR下,dMOPSO和MOPSOCD得到的最佳装机容量都相对较小,因为在FTL下系统几乎不产生多余的热能。然而,对于两层优化模型,TES的装机容量较大,这表明在两层优化模型下系统运行时产生更多的多余热能。对于GB的装机容量,通过dMOPSO得到的最佳值较小,因为系统的最佳GT装机容量较大。两层优化模型的最佳GB装机容量也较小,因为它具有较大的ST装机面积;
在案例3中,我们将两层优化模型与单一目标优化模型进行比较,以更深入地论证两层优化模型的卓越性;
表4单一目标优化模型的评估指标数值
表4显示单一目标的单层优化模型的评估指标数值。与之相比,两层优化模型显著提高能源、经济和环境性能。此外,两层优化模型在经济、能源和环境方面的最大改进分别达到34.95%、35.81%和37.55%;
两层优化模型和单一目标单层优化模型的最优配置解之间存在显著差异。单一目标单层优化模型具有特定的光伏(PV)和太阳能热(ST)安装面积,然而两层优化模型的ST占据全部屋顶面积。作为混合型CCHP的核心能源供应设备,燃气轮机(GT)在单一目标单层优化模型和两层优化模型之间的装机容量差异不大。单一目标单层优化模型具有较大的电网交互限制,这表明CCHP的独立运行性能较差。而两层优化模型使得CCHP能够实现脱网运行。对于GA的FTL和FTL-ECR、PSO的FTL以及DE的FTL-ECR,TES的装机容量较小。我们可以看到,在单一目标单层优化模型下,所有策略都具有较大的燃气锅炉(GB)装机容量,这表明GB在辅助供热方面发挥着重要作用;
图11展示在不考虑热能储存的情况下,两层优化模型的最优性能因子(PF)。可以观察到,即使在没有考虑热能储存的情况下,所提出的两层优化模型仍然提供一系列可行的解决方案,而非支配排序和相对决策度排序方法获得最优解。接下来,将比较考虑热能储存和不考虑热能储存的两层优化模型的最优解在性能和组件尺寸方面的特征,以及基于规则的运行策略;
表5不考虑TES的双层模型的最佳性能
表5描述在不考虑热能储存的情况下,两层优化模型的最优性能。与案例1中考虑热能储存的两层优化模型的最优性能相比,不考虑热能储存会显著降低CCHP的经济、能源和环境性能。此外,不考虑热能储存甚至会使CCHP的经济性能不如案例2中的多目标优化模型。这证明在两层优化模型中热能储存的重要性;
表6未考虑热能储存的两层模型的最优组件尺寸
表6提供在不考虑热能储存的情况下,两层优化模型的最优设计解。可以注意到,与考虑热能储存的两层优化模型相比(见第4.3节),不考虑热能储存会使CCHP的最优光伏(PV)装机面积显著增加,达到1214.55平方米.相比之下,太阳能热(ST)的装机面积减小至263.45平方米。这是因为此时太阳能热产生的多余热量没有热能储存器进行存储,导致最优太阳能热装机面积显著减小。燃气轮机(GT)的最优装机容量保持稳定,因为它是核心能源供应设备。太阳能热装机面积的显著减小导致最优燃气锅炉(GB)装机容量的显著增加,达到563.51千瓦。由于光伏装机面积的显著增加,光伏产生的多余电力不得不卖给电网,以保持系统处于稳定运行状态。最优电锅炉(EC)的装机容量减小至0千瓦,而最优空调(AC)的装机容量波动不大;
综上所述,所提出的两层优化模型可以显著提升混合型CCHP系统的经济、能源和环境性能。与基于规则的多目标单层优化模型相比,所提出的两层优化模型在经济、能源和环境性能方面最大可实现30.63%、37.94%和38.16%的改进。相较于基于规则的单目标单层优化模型,所提出的两层优化模型在经济、能源和环境性能方面最大可实现34.95%、35.81%和37.55%的改进。在两层优化模型下,热能储存(TES)对于混合型CCHP系统的发展至关重要。不考虑热能储存会使两层优化模型的最优系统经济、能源和环境性能达到328460.03美元、4385982.81千瓦、977244.11千克,这低于考虑热能储存的情况。此外,不考虑热能储存会导致CCHP的最优设计解与考虑热能储存时有很大不同;
需要强调的是,本发明所述的实施例是说明性的,而不是限定性的,因此本发明并不限于具体实施方式中所述的实施例,凡是由本领域技术人员根据本发明的技术方案得出的其他实施方式,同样属于本发明保护的范围。
Claims (3)
1.一种基于改进混沌优化算法的分布式能源系统双层运行优化方法,包括以下步骤:
步骤1,建立包含外层优化和内层优化的双层模型;
步骤2,求解双层优化模型,将参考建筑的年负荷数据和历史气象数据输入矩阵;通过k-means算法对矩阵中的数据进行聚类,从而得到聚类后的情景及其相应的情景概率。将典型情景及其相应的情景概率输入外层优化模型;
步骤3,外层优化模型随机初始化一个设备容量场景群,内层优化模型读取初始化的设备容量群。内层优化模型采用混合整数线性规划(MILP)算法,以系统运行成本最小为目标,获得最佳运行策略,并将最佳运行方案传递给外层优化模型;
步骤4,外层优化模型以内层优化模型得到的运行方案为基础,利用多目标混沌优化(MOCGO)算法进行位置更新,得到新的设备装机容量群,并将新的装机容量群返回给内层优化模型;
步骤5,内层优化模型利用新的装机容量群再次求解最优运行方案,如此循环迭代。当达到最大迭代次数时,循环将停止。然后使用非支配排序和相对决策度排序对解决方案集进行排序,选择非支配解集中相对决策度最大的作为最优个体;
步骤1,构建的冷热电三联供系统的主要设备包括电制冷机(EC)、吸收式制冷机(AC)、热交换装置(HE)、热能储存罐(TES)、燃气锅炉(GB)、燃气轮机(GT)、太阳能集热器(ST)和光伏装置(PV)。双层优化模型包括外层优化和内层优化两部分,外层优化模型如下:
(1)决策变量:
X=[NEC,NAC,NTES,NGB,APV,AST,NMT] (1)
(2)在外层优化模型中,经济、能源和环境目标函数如下:
式中,fEco、fEne和fEnv分别为经济、能源和环境目标函数。IV是利率,lifem是设备寿命,Cm是第m个设备的投资成本,Nm是第m个设备的装机容量。fGas是燃料成本,fGrid是电网互动成本,fMain是设备维护成本。它们的定义如下:
式中,S为情景总数,Ps为情景概率,Cf为天然气价格,FGB和FGT分别为GB和GT的天然气消耗量,EOUT为售电量,EIN为外购电量,Ce为使用时间,CEC为EC的制冷量,QGB为GB的制热量,EGT为GT的发电量,ε为设备维护成本系数;
式中,FGrid是电网的燃料消耗量,λf和λe分别是燃气和电网的等效CO2排放系数。λf取值为0.968kg/kWh,λe取值为0.220kg/kWh;
(3)设备容量的约束条件如表1所示:
表1设计参数范围
(4)内层优化模型的决策变量主要是各设备的输出功率,其定义如下:
X=[EMT(1:24·S),QGB(1:24·S),…,CEC(1:24·S)] (8)
(5)内层优化模型的目标是最小化系统的年度运行和维护成本,其定义如下:
fInner=fMain+fGrid+fGas (9)
(6)约束条件如下:
ELoad+EEC+EGrid,in=EMT+EPV+EGrid,out (10)
QLoad+QAC+QTES,in=QHR+QTES,out+QGB (11)
CLoad=CAC+CEC (12)
步骤2,求解双层优化模型,将参考建筑的年负荷数据和历史气象数据输入矩阵,通过k-means算法对矩阵中的数据进行聚类,从而得到聚类后的情景及其相应的情景概率。将典型情景及其相应的情景概率输入外层优化模型;
(1)将一年的负荷和气象数据存储到一个矩阵中,并对该矩阵进行归一化。归一化的公式如下:
式中,Ai是归一化后的元素;Xi是原始数据在归一化前的元素;Xmax和Xmin分别是归一化前数据的最大值和最小值;
(2)使用k均值对归一化后的矩阵进行聚类。我们获取典型场景的数量k以及每个簇的中心日。与中心日对应的原始数据是该簇的典型场景,而每个簇中包含的天数即为这些典型场景的天数。每个场景的概率如下:
式中,Pk是每个场景的概率,Nk表示每个场景中包含的天数;
在使用k均值算法对数据进行聚类之前,需要指定一个确定的簇数。然而,如果聚类值的选择基于用户的主观意愿,将不可避免地使最终的聚类结果更加随机和不确定。而且,即使对于相同的样本集,当不同用户使用不同数量的簇来处理数据时,最终的聚类结果可能会有很大的差异。因此,采用适当的方法确定簇的数量并确保最终聚类结果的客观性是非常重要的。作为确定最优簇数的客观方法,肘部法则的核心思想包括以下几点,随着样本分类数的增加,每个簇中样本的聚集程度将逐渐增加,平方误差和(SSE)自然会逐渐减小。当簇的数量k小于最优簇数时,随着簇的数量k的增加,SSE的减小趋势会逐渐减缓。当簇的数量k等于最优簇数时,SSE的减小趋势会突然减小。此时,继续增加k的值,SSE将不会有明显的变化。因此,使得SSE减小幅度突然减小的簇的数量k,就是数据的实际簇数。肘部法则的核心度量指标是平方误差和(SSE),具体定义如下:
式中,k表示簇的数量;Ci表示第i个簇;p表示第i个簇的样本点;mi表示第i个簇的质心;
步骤3,外层优化模型随机初始化一个设备容量场景群,内层优化模型读取初始化的设备容量群。内层优化模型采用MILP算法,以系统运行成本最小为目标,获得最佳运行策略,并将最佳运行方案传递给外层优化模型;
步骤4,外层优化模型以内层优化模型得到的运行方案为基础,利用MOCGO进行位置更新,得到新的设备装机容量群,并将新的装机容量群返回给内层优化模型。所选多目标智能算法为多目标混沌优化算法,混沌博弈优化算法(CGO)是一种基于混沌理论原理提出的优化算法,它利用分形和混沌博弈的基本概念。其初始化公式如下:
(1)在该算法中,每个候选解Xi包含一些决策变量Xi,j,这些变量表示西尔平斯基三角形中这些合格种子的位置。在这个算法中,西尔平斯基三角形被视为候选解的搜索空间。其数学模型如下:
式中,n是西尔平斯基三角形中合格点的数量,d是这些点的维数。考虑的合格点的初始位置是在搜索空间中随机确定的,具体过程如下:
式中,是合格点的初始位置;/>和/>分别是第i个候选解的第j个决策变量的最小和最大允许值;rand.()是一个在区间[0,1]内的随机数;
(2)第一个种子位置更新的数学表示如下:
式中,βi和γi分别表示1或2的随机整数,用于模拟掷骰子的可能性。GB是当前的最优解,MGi是初始合格点的平均值,它们被视为第一个临时三角形的三个顶点。αi是一个随机生成的矩阵,Xi是第i个候选解;
(3)类似于生成第一个种子的过程,第二个种子可以沿着MGi和Xi之间的直线移动,并利用一些随机生成的因素来约束该移动。该过程的数学表示如下:
(3)第三个种子是基于MGi的位置随机生成的,其数学模型如下:
(4)第四个种子的位置更新是基于对随机选择的决策变量进行一些随机更改。第四个种子的数学建模如下:
式中,R是在[0,1]上均匀分布的随机数;k是一个在[1,d]范围内的随机整数向量。此外,为控制种子的移动限制,对于α提出四个不同的公式,可以表示如下:
式中,rand是在[0,1]范围内的随机向量;δ和ε是在[0,1]范围内的随机向量;
(5)对于多目标混沌优化算法,进行初始适应度。对于每个个体Xj,其强度值的定义如下:
式中,Xj和Xw既是种群中的个体,也是存档个体。对于个体Xi来说,新的适应度值由两部分组成,密度信息D(Xi)和原始适应度值R(Xi)。R(Xi)的定义如下:
式中,R(Xi)个体在群体和外部档案中的主导地位决定Xi个体在群体和外部档案中的主导地位;
(6)关于基于移位的密度估计,如果某个个体在某个目标上的表现优于p,那么它就会被转移到与p在该目标上相同的位置;否则,它就会保持不变。个体p在种群D′(p,P)中的新密度可以表示为:
D′(p,P)=D(dist(p,q′1),dist(p,q′2),…dist(p,q′N-1)) (25)
式中,N为种群大小;dist(p,q′1)为两个个体相似度的程度;q′1为个体qi(qi≠p and qi∈P)的位移版本,其定义如下:
式中,m代表目标数,p(j)、qi(j)和分别代表个体p、q和q′的第j个目标值;
步骤5,内层优化模型利用新的装机容量群再次求解最优运行方案,如此循环迭代。当达到最大迭代次数时,循环将停止;
对于最优个体的选取,采用非支配排序和相对决策度排序获得帕累托前沿和最优个体,非支配排序可以筛选出不同个体的适应度优劣级别,相对决策度可以表达解集的分布情况。下面介绍基于相对决策度的非支配排序机制:
(1)对解集进行非支配排序,按照不同支配等级对解集进行划分。对于种群每次进化产生的帕累托前沿,采用外部存档的方法进行存储与更新;
(2)对于外部档案中划分好支配等级的解集,依据支配等级依次计算不同解的相对决策度并按降序排列,选择非支配解集中相对决策度最大的作为最优个体。
2.按照权利要求1所述一种相对决策度的计算,其特征在于:
相对决策度的计算公式如下:
式中,为第i个解在搜索空间中的决策度,Nq表示同一支配等级的解的数量,/>和为第i个解在第j个维度最接近的两个维度值,且满足/>和/>为同支配等级的解中第j个维度的最大值和最小值。RDMPi为第i个解在目标空间中的相对决策度,Nobj表示目标函数的数量;/>和/>为同支配等级的解中第j个目标函数的最大值和最小值。
3.按照权利要求1所述一种基于改进混沌优化算法的分布式能源系统双层运行优化方法,其特征在于:
所述多目标智能算法并非特指某种算法,可根据性能选择相适应的算法进行优化。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202410059667.2A CN117874562A (zh) | 2024-01-16 | 2024-01-16 | 一种基于改进混沌优化算法的分布式能源系统双层运行优化方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202410059667.2A CN117874562A (zh) | 2024-01-16 | 2024-01-16 | 一种基于改进混沌优化算法的分布式能源系统双层运行优化方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN117874562A true CN117874562A (zh) | 2024-04-12 |
Family
ID=90591751
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202410059667.2A Pending CN117874562A (zh) | 2024-01-16 | 2024-01-16 | 一种基于改进混沌优化算法的分布式能源系统双层运行优化方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN117874562A (zh) |
-
2024
- 2024-01-16 CN CN202410059667.2A patent/CN117874562A/zh active Pending
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Xu et al. | Smart energy systems: A critical review on design and operation optimization | |
CN109165788B (zh) | 一种冷热电联供系统的优化方法 | |
Yan et al. | Stochastic multi-scenario optimization for a hybrid combined cooling, heating and power system considering multi-criteria | |
Wang et al. | Synergistic planning of an integrated energy system containing hydrogen storage with the coupled use of electric-thermal energy | |
CN112036934A (zh) | 考虑热电协调运行的负荷聚合商参与需求响应的报价方法 | |
Chai et al. | Differential evolution-based system design optimization for net zero energy buildings under climate change | |
CN112418488B (zh) | 一种基于两阶段能量优化的综合能源系统调度方法及装置 | |
Liu et al. | Two-layer multiple scenario optimization framework for integrated energy system based on optimal energy contribution ratio strategy | |
CN116468215A (zh) | 计及源荷不确定性的综合能源系统调度方法及装置 | |
CN116187173A (zh) | 一种区域综合能源系统双层优化方法及系统 | |
Xu et al. | A data-based day-ahead scheduling optimization approach for regional integrated energy systems with varying operating conditions | |
CN114742276A (zh) | 一种考虑㶲效率的含orc的园区综合能源系统多目标优化调度方法 | |
CN113708418A (zh) | 微电网优化调度方法 | |
Li et al. | Multi-dimension day-ahead scheduling optimization of a community-scale solar-driven CCHP system with demand-side management | |
Zhou et al. | Multi-objective station-network synergy planning for regional integrated energy system considering energy cascade utilization and uncertainty | |
Ren et al. | Economic, energy and environmental analysis and evaluation of hybrid CCHP system considering different buildings: A two-level optimization model | |
Su et al. | Optimal placement and capacity sizing of energy storage systems via NSGA-II in active distribution network | |
Hongwei et al. | Robust stochastic optimal dispatching of integrated electricity-gas-heat systems with improved integrated demand response | |
Abdalla et al. | The impact of clustering strategies to site integrated community energy and harvesting systems on electrical demand and regional GHG reductions | |
Ren et al. | Design and analysis of solar hybrid combined cooling, heating and power system: A bi-level optimization model | |
CN116362421A (zh) | 一种基于能源综合统筹分析的供能分配预测系统及其方法 | |
CN117874562A (zh) | 一种基于改进混沌优化算法的分布式能源系统双层运行优化方法 | |
Hu et al. | Robust optimal scheduling of integrated energy systems considering multiple uncertainties | |
Altork et al. | Optimizing Hybrid Heating Systems: Identifying Ideal Stations and Conducting Economic Analysis Heating Houses in Jordan. | |
Chai et al. | Differential Evolution-based System for Net-zero Energy Buildings Under Climate Change |
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 |