CN112861292B - 一种电-气综合能源系统恢复改善方法 - Google Patents
一种电-气综合能源系统恢复改善方法 Download PDFInfo
- Publication number
- CN112861292B CN112861292B CN202110036776.9A CN202110036776A CN112861292B CN 112861292 B CN112861292 B CN 112861292B CN 202110036776 A CN202110036776 A CN 202110036776A CN 112861292 B CN112861292 B CN 112861292B
- Authority
- CN
- China
- Prior art keywords
- gas
- natural gas
- recovery
- stage
- node
- 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
- 238000011084 recovery Methods 0.000 title claims abstract description 160
- 238000000034 method Methods 0.000 title claims abstract description 40
- 230000006872 improvement Effects 0.000 title claims abstract description 17
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 claims abstract description 422
- 239000007789 gas Substances 0.000 claims abstract description 231
- 239000003345 natural gas Substances 0.000 claims abstract description 226
- 238000005457 optimization Methods 0.000 claims abstract description 17
- 238000012545 processing Methods 0.000 claims abstract description 12
- 238000012423 maintenance Methods 0.000 claims description 23
- 230000005540 biological transmission Effects 0.000 claims description 14
- 238000006243 chemical reaction Methods 0.000 claims description 13
- 238000004364 calculation method Methods 0.000 claims description 12
- 238000012546 transfer Methods 0.000 claims description 11
- 230000001419 dependent effect Effects 0.000 claims description 7
- 150000001875 compounds Chemical class 0.000 claims description 6
- 239000004243 E-number Substances 0.000 claims description 3
- 230000008439 repair process Effects 0.000 abstract description 2
- 230000008878 coupling Effects 0.000 description 3
- 238000010168 coupling process Methods 0.000 description 3
- 238000005859 coupling reaction Methods 0.000 description 3
- 239000000446 fuel Substances 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 238000007796 conventional method Methods 0.000 description 2
- 230000005611 electricity Effects 0.000 description 2
- 238000010248 power generation Methods 0.000 description 2
- 230000008719 thickening Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000003990 capacitor Substances 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 239000001257 hydrogen Substances 0.000 description 1
- 229910052739 hydrogen Inorganic materials 0.000 description 1
- 125000004435 hydrogen atom Chemical class [H]* 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/18—Network design, e.g. design based on topological or interconnect aspects of utility systems, piping, heating ventilation air conditioning [HVAC] or cabling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- 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
-
- 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
-
- 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/48—Controlling the sharing of the in-phase component
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/04—Constraint-based 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
- 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]
-
- 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
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A30/00—Adapting or protecting infrastructure or their operation
- Y02A30/60—Planning or developing urban green infrastructure
Landscapes
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Economics (AREA)
- Strategic Management (AREA)
- Human Resources & Organizations (AREA)
- Tourism & Hospitality (AREA)
- Geometry (AREA)
- General Business, Economics & Management (AREA)
- Marketing (AREA)
- Entrepreneurship & Innovation (AREA)
- Development Economics (AREA)
- Operations Research (AREA)
- Computer Hardware Design (AREA)
- Health & Medical Sciences (AREA)
- Game Theory and Decision Science (AREA)
- Power Engineering (AREA)
- Evolutionary Computation (AREA)
- Quality & Reliability (AREA)
- General Engineering & Computer Science (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- Computational Mathematics (AREA)
- Educational Administration (AREA)
- Mathematical Analysis (AREA)
- Public Health (AREA)
- Water Supply & Treatment (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种电‑气综合能源系统恢复改善方法。建立电‑气综合能源系统网络结构拓扑模型;在电‑气综合能源系统网络结构拓扑模型下,确定初始系统运行状态;建立电‑气综合能源系统恢复序列优化模型;根据恢复序列处理得到电‑气综合能源系统恢复力参数,获得电‑气综合能源系统恢复情况,并根据恢复情况进行判断实施改善。本发明能够得到电‑气综合能源系统大范围故障后电力线路和天然气管道的恢复序列,指导系统施工人员按照所得到的恢复序列修复电力线路和天然气管道,快速恢复系统中的电力负荷和天然气负荷,从而有效改善电‑气综合能源系统整体的恢复能力。
Description
技术领域
本发明属于电力系统、综合能源系统技术领域的一种电-气综合能源系统恢复优化方法,涉及实现电力系统与天然气系统深度耦合情形下的电-气综合能源系统恢复力评估与改善方法。
背景技术
近年来,燃气机组装机容量不断增大,燃气机组的发电量在电力系统总发电量中占据了很大的比重,导致电力系统的运行状况更加依赖于天然气系统。与此同时,随着电转气设备的引入,天然气系统的产气量也逐渐依赖于电力系统的电力负荷。然而,电力系统和天然气系统的这种相互依赖特性,使得两个系统的恢复过程也更加复杂。一方面,电力系统中燃气机组的恢复不仅取决于电力系统自身的恢复情况,还取决于来自天然气系统的天然气燃料供应情况;另一方面,电转气设备的恢复也取决于电力系统中的电力负荷是否供给。
目前,复杂工程系统的恢复方法主要是针对单一工程系统的元件恢复序列进行优化,如基于复杂网络理论或电力系统运行分析理论研究电力系统的恢复方式,然后提出改善系统恢复能力的措施。但是这些方法仅考虑电力系统独立运行的场景,未考虑与电力系统紧密耦合的其他工程系统的影响,在恢复方式优化方面的研究成果无法有效应用于综合能源系统。
因此,亟需综合考虑电力系统和天然气系统的相互依赖特性,兼顾电力系统、天然气系统自身内部的恢复特性以及两个系统之间相互依赖、相互影响的特性,提出一种电-气综合能源系统恢复方式优化方法,实现电-气综合能源系统整体的恢复能力优化,通过增加系统维修设备数量、加粗电力线路、加厚天然气管道等措施,有效提升电-气综合能源系统的恢复能力。
发明内容
为了解决背景技术中存在的问题,本发明的目的是提供一种电-气综合能源系统恢复改善方法。
本发明能够快速计算获得电-气综合能源系统大范围故障后电力线路和天然气管道的恢复序列,系统维修设备通过按照所得到的恢复序列修复电力线路和天然气管道,可以快速恢复电-气综合能源系统中的电力负荷和天然气负荷,有效提升电-气综合能源系统的恢复能力。
为实现上述目的,本发明采用的具体技术方案是包括以下步骤:
步骤1、建立电-气综合能源系统网络结构拓扑模型;
步骤2、在电-气综合能源系统网络结构拓扑模型下,确定初始系统运行状态,包括:初始正常电力线路和天然气管道集合、待恢复输电线路和天然气管道集合、电力系统和天然气系统维修设备数量、电-气综合能源系统所需的恢复阶段数;
步骤3、建立电-气综合能源系统恢复序列优化模型,具体包括:恢复目标函数、恢复资源约束、电力系统运行约束、天然气系统运行约束、电力-天然气相互依赖约束,采用反向贪心算法得到电力线路和天然气管道的恢复序列;
步骤4、根据恢复序列处理得到电-气综合能源系统恢复力参数,获得电-气综合能源系统恢复情况,并根据恢复情况进行判断实施改善。
所述的电-气综合能源系统是由电力系统、天然气系统、电力-天然气相互依赖环节三个部分组成。
电力系统包含电节点和电力线路。
电节点指电力系统中的发电机组、变电站、电负荷设备,其中,发电机组包含燃气机组和非燃气机组(如燃煤机组、核电机组、水电机组等)两种,电负荷设备消耗的电力负荷包含常规电负荷和转气负荷两种,转气负荷用于电转气设备的正常工作,常规电负荷用于与天然气系统无关的电负荷设备的正常工作。电力线路指电力系统中的输电线路、变压器支路。电节点之间通过电力线路连接,每条电力线路首尾两端的电节点分别称为首端节点和末端节点。
天然气系统包含气节点和天然气管道。气节点包含天然气系统中的气源设备、压缩机设备、气负荷设备,其中,气源设备包括常规气源和电转气设备两种,气负荷设备消耗的天然气负荷包含常规气负荷和转电负荷两种,转电负荷用于燃气机组发电,常规气负荷用于与电力系统无关的气负荷设备的正常工作。天然气管道指天然气系统中的输气管道和压缩机支路。气节点之间通过天然气管道连接,每条天然气管道首尾两端的气节点分别称为首端节点和末端节点。
电力-天然气相互依赖环节包括燃气机组依赖链接和电转气设备依赖链接,燃气机组依赖链接是指燃气机组在电力系统中所在的电节点与在天然气系统中所在的气节点之间的链接;电转气设备依赖链接是指电转气设备在电力系统中所在的电节点与在天然气系统中所在的气节点之间的链接。燃气机组发电所消耗的天然气燃料依赖于与天然气系统的气节点;电转气设备维持正常工作所需要的电力负荷依赖于电力系统中的电节点。
所述的步骤1,具体如下:
所述的电-气综合能源系统网络结构拓扑模型包含了电力系统网络结构拓扑模型、天然气系统网络结构拓扑模型和电力-天然气相互依赖环节网络结构拓扑模型;
步骤1.1、建立电力系统网络结构拓扑模型:
Ge=(I,L)(1)
式中,GE表示电力系统网络结构拓扑模型,I表示电节点集合,L表示电力线路集合;
步骤1.2、建立天然气系统网络结构拓扑模型:
Gg=(J,P)(2)
式中,GG表示天然气网络结构拓扑模型,J表示气节点集合,P表示天然气管道集合;
步骤1.3、建立电力-天然气相互依赖环节网络结构拓扑模型:
Ξ=(Ξgfu,Ξptg)(3)
式中,Ξ表示电力-天然气相互依赖环节网络结构拓扑模型;Ξgfu表示燃气机组依赖链接集合;Ξptg表示电转气设备依赖链接集合。
所述步骤2,具体如下:
步骤2.1、确定初始正常电力线路集合Nl和初始正常天然气管道集合Np;
步骤2.2、确定待恢复电力线路集合Nl和天然气管道集合Np;
步骤2.3、确定电力系统维修设备数量Re和天然气系统维修设备数量Rg;
步骤2.4、确定电力系统所需的恢复阶段数Te、天然气系统所需的恢复阶段数Tg,得到电-气综合能源系统所需的总恢复阶段数:
T=max{Te,Tg}(4)
式中,T表示电-气综合能源系统所需的总恢复阶段数;Te表示电力系统所需的恢复阶段数;Tg表示天然气系统所需的恢复阶段数;max{Te,Tg}表示取Te和Tg两者中的最大值。
所述步骤3,具体如下:
步骤3.1、以电力负荷和天然气负荷恢复量的加权之和最大为目标Max r(t),建立第t阶段内电-气综合能源系统恢复序列优化的目标函数:
式中,r(t)表示-气综合能源系统恢复序列优化的目标函数;t表示恢复阶段的序数;υ表示电力系统恢复权重系数;Ne表示电力系统中电节点的总数;Ng表示天然气系统中气节点的总数;i表示电节点的序数;j表示气节点的序数;表示第t阶段电节点i上的电力负荷;表示系统发生故障前电节点i上的电力负荷;表示第t阶段气节点j上的天然气负荷;表示系统发生故障前气节点j上的天然气负荷;
所述的电力系统恢复权重系数表示电-气综合能源系统运行人员对电力负荷恢复的重视程度,取值范围在0到1之间。电力系统恢复权重系数越大,表示电-气综合能源系统运行人员越重视对电力负荷的恢复,而对天然气负荷的恢复越不重视。
步骤3.2、建立第t阶段内电-气综合能源系统恢复序列优化的约束条件:
步骤3.2.1、建立第t阶段内的恢复资源约束:
式中,表示第t阶段电力线路l的恢复决策变量,为0-1整型变量,当第t阶段电力线路l被修复时,否则l表示电力线路l的序数;Nl表示电力系统中电力线路的总数;表示第t阶段天然气管道p的恢复决策变量,为0-1整型变量,当第t阶段天然气管道被修复时,否则p表示天然气管道的序数;Np表示天然气系统中天然气管道的总数;
式(6)为可恢复的电力线路数量不能超过电力系统维修设备数量Re约束,式(7)为可恢复的天然气管道数量不能超过天然气系统维修设备数量Rg约束,式(8)为电力线路和天然气管道恢复决策变量满足0-1整型变量约束。
步骤3.2.2、建立第t阶段内的电力系统运行约束:
式中,表示第t阶段电力线路l的状态变量,为0-1整型变量,当第t阶段电力线路l为正常状态时,否则m表示电力系统恢复阶段的序数;表示第t阶段电节点i的状态变量,为0-1整型变量,当第t阶段电节点i为正常状态时,否则l∈i表示与电节点i相连的电力线路集合;表示第t阶段电节点i上的非燃气机组的出力;表示第t阶段电节点i上的燃气机组的出力;表示第t阶段电力线路l流过的潮流;l∈Ξ1(i)表示所有首端节点与电节点i相连的电力线路的集合;l∈Ξ2(i)表示所有末端节点与电节点i相连的电力线路的集合;表示第t阶段电节点i上的常规电负荷;表示第t阶段电节点i上的转气负荷;和分别表示第t阶段电力线路l首端节点和末端节点的电压相角;xbranch,l表示电力线路l的电抗;M表示一个很大的数;和分别表示电节点i上非燃气机组的最小技术出力和最大技术出力;和分别表示电节点i上燃气机组的最小技术出力和最大技术出力;表示电力线路l的传输容量;
式(9)为电力线路状态变量与恢复决策变量关系约束,式(10)为电力线路状态变量与电节点状态变量关系约束,式(11)为电力线路和电节点状态变量满足0-1整型变量约束,式(12)为电节点潮流平衡约束,式(13)为电节点负荷功率约束,式(14)为电力线路潮流方程约束,式(15)-(16)分别为非燃气机组和燃气机组出力上下限约束,式(17)为电力线路传输容量约束,式(18)为电节点可恢复电力负荷约束。
步骤3.2.3、建立第t阶段内的天然气系统运行约束:
式中,表示第t阶段天然气管道p的状态变量,为0-1整型变量,当第t阶段天然气管道p为正常状态时,否则n表示天然气系统恢复阶段的序数;表示第t阶段气节点j的状态变量,为0-1整型变量,当第t阶段气节点j为正常状态时,否则p∈j表示与气节点j相连的天然气管道集合;表示第t阶段气节点j上的常规气源的出力;表示第t阶段气节点j上的电转气设备的出力;表示第t阶段天然气管道p流过的气流;p∈Ξ1(j)表示所有首端节点与气节点j相连的天然气管道的集合;p∈Ξ2(j)表示所有末端节点与气节点j相连的天然气管道的集合;表示第t阶段气节点j上的常规气负荷;表示第t阶段气节点j上的转电负荷;和分别表示第t阶段天然气管道p首端节点和末端节点的气压;Cp表示天然气管道p的传输系数;和分别表示气节点j上的常规气源的最小技术出力和最大技术出力;和分别表示气节点j上的电转气设备的最小技术出力和最大技术出力;和分别表示气节点j上的气压下限值和上限值;表示第t阶段气节点j上的气压;
式(19)为天然气管道状态变量与恢复决策变量关系约束,式(20)为天然气管道状态变量与气节点状态变量关系约束,式(21)为天然气管道和气节点状态变量满足0-1整型变量约束,式(22)为气节点气流平衡约束,式(23)为气节点负荷功率约束,式(24)为天然气管道气流方程约束,式(25)-(26)分别为常规气源和电转气设备出力上下限约束,式(27)为天然气管道传输容量约束,式(28)为气节点气压上下限约束,式(29)为气节点可恢复天然气负荷约束。
步骤3.2.4、建立第t阶段内的电力-天然气相互依赖约束:
式中,ηgfu表示燃气机组的工作效率;ηptg表示电转气设备的工作效率;j∈Ξgfu(i)表示电节点i上的转气负荷向气节点j上的电转气设备供电;i∈Ξptg(j)表示气节点j上的转电负荷向电节点i上的燃气机组供气;
式(30)为燃气机组气流消耗量与电功率产出量的转化约束;式(31)为电转气设备电功率消耗量与气流产出量的转化约束。
步骤3.3、基于反向贪心算法,对每一阶段的恢复决策进行计算。从当前待恢复电力线路、天然气管道集合中,找出使步骤3.1中优化目标值最小且满足步骤3.2中约束条件的电力线路、天然气管道作为本阶段的恢复对象;
步骤3.4、将步骤3.3中得到的恢复对象从当前待恢复电力线路、天然气管道集合中移除,判断当前待恢复电力线路、天然气管道集合是否为空集:
如果是,则说明所有待恢复电力线路、天然气管道已经被全部恢复,基于各阶段的恢复决策得到电力线路、天然气管道的恢复序列,进行步骤4;
否则,返回步骤3.3,更新当前待恢复电力线路、天然气管道集合,进行下一阶段的计算。
所述步骤4中,具体如下:
步骤4.1、根据步骤3中得到的恢复序列,计算各阶段恢复的电力负荷,按照以下公式处理获得电力系统恢复力参数:
步骤4.2、根据步骤3中得到的恢复序列,计算各阶段恢复的天然气负荷,按照以下公式处理获得天然气系统恢复力参数:
步骤4.3、根据步骤4.1、步骤4.2中得到的电力系统恢复力参数和天然气系统恢复力参数,按照以下公式处理获得电-气综合能源系统恢复力参数:
Rres=υRres,e+(1-υ)Rres,g(34)
步骤4.4、根据计算得到电-气综合能源系统恢复力结果实施改善:
若电-气综合能源系统的恢复力参数值小于预设的恢复力阈值,则电-气综合能源系统恢复力水平低,按照以下措施进行改善:
对于电力系统,增加电力系统维修设备数量以增加每一阶段可恢复的电力线路数量,加粗电力线路以避免电力线路故障;
对于天然气系统,增加天然气系统维修设备数量以增加每一阶段可恢复的天然气管道数量,加厚天然气管道以避免天然气管道故障。
否则电-气综合能源系统恢复力水平高,不进行操作。
通过上述方法完整实施,加快电-气综合能源系统中的电力负荷和天然气负荷的恢复速度,有效改善电-气综合能源系统整体的恢复能力。
在上述电-气综合能源系统恢复改善方法中,步骤3.2.2中:所述式(9),即电力线路状态变量与恢复决策变量关系约束,具体为:只有电力线路的恢复决策变量取值为1时,该电力线路的状态变量才可能取值为1;如果电力线路的恢复决策变量取值为0时,该电力线路的状态变量也取值为0。
在上述电-气综合能源系统恢复改善方法中,步骤3.2.2中:所述式(10),即电力线路状态变量与电节点状态变量关系约束,具体为:如果电节点所相连的电力线路中,至少有一条电力线路的状态变量为1,则该电节点的状态变量为1;如果电节点所相连的电力线路的状态变量均为0,则该电节点的状态变量为0。
在上述电-气综合能源系统恢复改善方法中,步骤3.2.3中:所述式(19),即天然气管道状态变量与恢复决策变量关系约束,具体为:只有天然气管道的恢复决策变量取值为1时,该天然气管道的状态变量才可能取值为1;如果天然气管道的恢复决策变量取值为0时,该天然气管道的状态变量也取值为0。
在上述电-气综合能源系统恢复改善方法中,步骤3.2.3中:所述式(20),即天然气管道状态变量与气节点状态变量关系约束,具体为:如果气节点所相连的天然气管道中,至少有一条天然气管道的状态变量为1,则该气节点的状态变量为1;如果气节点所相连的天然气管道的状态变量均为0,则该气节点的状态变量为0。
在上述电-气综合能源系统恢复改善方法中,步骤3.3中:所述的反向贪心算法,具体为:从第T阶段到第1阶段反向进行,对每一阶段的恢复决策进行计算,得到各阶段的恢复决策,基于各阶段的恢复决策得到恢复序列,然后从第1阶段到第T阶段正向进行恢复。
在上述电-气综合能源系统恢复改善方法中,步骤4.3中:所述的电-气综合能源系统恢复力参数,反映了电-气综合能源系统的恢复能力。电-气综合能源系统恢复力参数越大,表明电力负荷和天然气负荷的恢复量越大,电-气综合能源系统的恢复能力越强。
本发明具有的有益效果如下:
本发明的方法为一种在电力系统与天然气系统深度耦合情形下,实现电-气综合能源系统大范围故障后的恢复改善方法,克服以往恢复方法仅考虑单一能源系统独立运行下的恢复力进行分析和改善的不足;本发明方法可以兼顾电力系统、天然气系统自身内部的恢复特性以及两个系统之间相互依赖、相互影响的特性,实现电-气综合能源系统整体的恢复方式优化,有效提升电-气综合能源系统的恢复能力。
本发明能够满足未来电力系统和天然气系统深度耦合情形下的工程应用需求,快速对大范围故障后电-气综合能源系统中的电力负荷和天然气负荷进行恢复。对于电力系统,增加电力系统维修设备数量以增加每一阶段可恢复的电力线路数量,加粗电力线路以避免电力线路故障;对于天然气系统,增加天然气系统维修设备数量以增加每一阶段可恢复的天然气管道数量,加厚天然气管道以避免天然气管道故障。通过上述措施,可以加快电-气综合能源系统中的电力负荷和天然气负荷的恢复速度,有效改善电-气综合能源系统整体的恢复能力,最大程度地减小大范围故障造成的损失。
附图说明
图1是本发明的电-气综合能源系统的实施例结构示意图;
图2是本发明方法流程图。
图3是本发明的电-气综合能源系统中电力负荷恢复计算结果图。
图4是本发明的电-气综合能源系统中天然气负荷恢复计算结果图。
具体实施方式
下面结合附图和实施例对本发明作进一步说明。
按照本发明内容完整方法实施的具体实施例如下:
以IEEE30节点电力系统和比利时20节点天然气系统组成的电-气综合能源系统为例,结合技术方案和附图详细叙述本发明的具体实施。
如图1所示,电-气综合能源系统包含电力系统、天然气系统、电力-天然气相互依赖环节。其中,电力系统包含30个电节点和41条电力线路,天然气系统包含20个气节点和19条天然气管道。令电节点5、7、8上的发电机组为燃气机组,需要的天然气燃料分别由气节点6、7、10上的转电负荷提供,其余发电机组为非燃气机组。令气节点1、2、5、8上的气源为电转气设备,需要的电力负荷分别由电节点2、7、5、14上的转气负荷提供,其余气源为常规气源。
另外,为了方便计算,对比利时20节点系统作如下修改:将双回天然气管道合并为单回;以初始运行状态下天然气管道流过的天然气流量的1.5倍作为天然气管道的传输容量,不满8×103m3·h-1按8×103m3·h-1计算。对IEEE30节点系统作如下修改:将所有线路均简化为无向有权边,不计并联电容支路;以初始运行状态下电力线路流过潮流的1.5倍作为电力线路的传输容量,不满30MW按30MW计算。
本发明实施流程如图2所示,具体步骤如下:
步骤1、建立电-气综合能源系统网络结构拓扑模型,具体如下:
步骤1.1、建立电力系统网络结构拓扑模型;
步骤1.2、建立天然气系统网络结构拓扑模型;
步骤1.3、建立电力-天然气相互依赖环节网络结构拓扑模型;
步骤2、在电-气综合能源系统网络结构拓扑模型下,确定初始系统运行状态,包括:初始正常电力线路和天然气管道集合、待恢复电力线路和天然气管道集合、电力系统和天然气系统维修设备数量、电-气综合能源系统所需的恢复阶段数,具体如下:
步骤2.1、确定初始正常电力线路集合Nl和初始正常天然气管道集合Np;
步骤2.2、确定待恢复电力线路集合Nl和待恢复天然气管道集合Np;
其中,初始正常电力线路集合Nl为图1电力系统区域中方框包含的4条电力线路,其余电力线路组成了待恢复电力线路集合Nl;初始正常天然气管道集合Np为图1天然气系统区域中方框包含的1条天然气管道,其余天然气管道组成了待恢复天然气管道集合Np。
步骤2.3、确定电力系统维修设备数量Re和天然气系统维修设备数量Rg;
其中,电力系统维修设备数量Re取为1,天然气系统维修设备数量Rg取为1。
步骤2.4、确定电力系统所需的恢复阶段数Te、天然气系统所需的恢复阶段数Tg,得到电-气综合能源系统所需的总恢复阶段数T;
其中,电力系统所需的恢复阶段数Te为37,天然气系统所需的恢复阶段数Tg为18,电-气综合能源系统所需的总恢复阶段数T为37。
所述步骤3,建立电-气综合能源系统恢复序列优化模型,具体包括:恢复目标、恢复资源约束、电力系统运行约束、天然气系统运行约束、相互依赖约束,基于反向贪心算法,得到电力线路和天然气管道的恢复序列,具体如下:
步骤3.1、以电力负荷和天然气负荷恢复量的加权之和最大为目标Max r(t),建立第t阶段内电-气综合能源系统恢复序列优化的目标函数:
式中,r(t)表示-气综合能源系统恢复序列优化的目标函数;t表示恢复阶段的序数;υ表示电力系统恢复权重系数;Ne表示电力系统中电节点的总数;Ng表示天然气系统中气节点的总数;i表示电节点的序数;j表示气节点的序数;表示第t阶段电节点i上的电力负荷;表示系统发生故障前电节点i上的电力负荷;表示第t阶段气节点j上的天然气负荷;表示系统发生故障前气节点j上的天然气负荷;
其中,电力系统恢复权重系数υ取为0.5。
步骤3.2、建立第t阶段内电-气综合能源系统恢复序列优化的约束条件:
步骤3.2.1、建立第t阶段内的恢复资源约束;
步骤3.2.2、建立第t阶段内的电力系统运行约束;
步骤3.2.3、建立第t阶段内的天然气系统运行约束;
步骤3.2.4、建立第t阶段内的电力-天然气相互依赖约束;
步骤3.3、基于反向贪心算法,对每一阶段的恢复决策进行计算。从当前待恢复电力线路、天然气管道集合中,找出使步骤3.1中优化目标值最小且满足步骤3.2中约束条件的电力线路、天然气管道作为本阶段的恢复对象;
步骤3.4、将步骤3.3中得到的恢复对象从当前待恢复电力线路、天然气管道集合中移除。判断当前待恢复电力线路、天然气管道集合是否为空集,如果是,则说明所有待恢复电力线路、天然气管道已经被全部恢复,基于各阶段的恢复决策得到电力线路、天然气管道的恢复序列,进行步骤4;否则,返回步骤3.3,更新当前待恢复电力线路、天然气管道集合,进行下一阶段的计算;
步骤4中,建立电-气综合能源系统恢复力参数,计算得到电-气综合能源系统恢复力结果并实施改善,具体如下:
步骤4.1、根据步骤3中得到的恢复序列,计算各阶段恢复的电力负荷,如图3中的黑线所示,按照以下公式处理获得电力系统恢复力参数,结果如表1所示;
为了说明本发明方法的优点,利用常规方法(每阶段选择传输容量最大的电力线路进行恢复)计算得到的各阶段恢复的电力负荷如图3中的灰线所示。比较图3中的黑线和灰线,可以看出本发明方法得到的电力线路恢复序列能够更快速地恢复电力负荷,提高电力系统的恢复能力。
步骤4.2、根据步骤3中得到的恢复序列,计算各阶段恢复的天然气负荷,如图4的黑线所示,按照以下公式处理获得天然气系统恢复力参数,结果如表1所示;
为了说明本发明方法的优点,利用常规方法(每阶段选择传输容量最大的天然气管道进行恢复)计算得到的各阶段恢复的天然气负荷如图4中的灰线所示。比较图3中的黑线和灰线,可以看出本发明方法得到的天然气管道恢复序列能够更快速地恢复天然气负荷,提高天然气系统的恢复能力。
步骤4.3、根据步骤4.1、步骤4.2中得到的电力系统恢复力参数和天然气系统恢复力参数,按照以下公式处理获得电-气综合能源系统恢复力参数,结果如表1所示;
表1恢复力参数计算结果
从表1可以看出,电力系统的恢复力结果大于天然气系统的恢复力计算结果,表明电力系统的恢复能力由于天然气系统的恢复能力,这是因为电力系统的网络拓扑为环状,能够通过网络重构实现电力负荷的快速恢复,而天然气系统的网络拓扑为辐射状,无法过网络重构实现天然气负荷的快速恢复。
步骤4.4、根据计算得到电-气综合能源系统恢复力结果实施改善:若电-气综合能源系统的恢复力参数值过小,表明系统恢复力水平较低。
设电-气综合能源系统恢复力目标值为0.9,通过表1可以看出电-气综合能源系统恢复力计算结果为0.785,低于设定的目标值,需要通过实施改善措施来提高电-气综合能源系统恢复力。
对于电力系统,增加电力系统维修设备数量以增加每一阶段可恢复的电力线路数量,将步骤2.3中的电力系统维修设备数量Re改为2,增加天然气系统维修设备数量以增加每一阶段可恢复的天然气管道数量,将步骤2.3中的天然气系统维修设备数量Rg改为2。同时,加粗部分电力线路以避免电力线路故障,加厚部分天然气管道以避免天然气管道故障。重新计算电-气综合能源系统恢复力参数,结果如表2所示。
表2实施改进措施后的恢复力参数计算结果
比较表1和表2可以看出,通过上述措施,提高了电-气综合能源系统恢复力结果,高于设定的目标值,有效改善了电-气综合能源系统整体的恢复能力。
Claims (3)
1.一种电-气综合能源系统恢复改善方法,其特征在于:
步骤1、建立电-气综合能源系统网络结构拓扑模型;
步骤2、在电-气综合能源系统网络结构拓扑模型下,确定初始系统运行状态,包括:初始正常电力线路和天然气管道集合、待恢复输电线路和天然气管道集合、电力系统和天然气系统维修设备数量、电-气综合能源系统所需的恢复阶段数;
步骤3、建立电-气综合能源系统恢复序列优化模型,具体包括:恢复目标函数、恢复资源约束、电力系统运行约束、天然气系统运行约束、电力-天然气相互依赖约束,采用反向贪心算法得到电力线路和天然气管道的恢复序列;
步骤4、根据恢复序列处理得到电-气综合能源系统恢复力参数,获得电-气综合能源系统恢复情况,并根据恢复情况进行判断实施改善;
所述的步骤1,具体如下:
所述的电-气综合能源系统网络结构拓扑模型包含了电力系统网络结构拓扑模型、天然气系统网络结构拓扑模型和电力-天然气相互依赖环节网络结构拓扑模型;
步骤1.1、建立电力系统网络结构拓扑模型:
Ge=(I,L) (1)
式中,GE表示电力系统网络结构拓扑模型,I表示电节点集合,L表示电力线路集合;
步骤1.2、建立天然气系统网络结构拓扑模型:
Gg=(J,P) (2)
式中,GG表示天然气网络结构拓扑模型,J表示气节点集合,P表示天然气管道集合;
步骤1.3、建立电力-天然气相互依赖环节网络结构拓扑模型:
Ξ=(Ξgfu,Ξptg) (3)
式中,Ξ表示电力-天然气相互依赖环节网络结构拓扑模型;Ξgfu表示燃气机组依赖链接集合;Ξptg表示电转气设备依赖链接集合;
所述步骤3,具体如下:
步骤3.1、以电力负荷和天然气负荷恢复量的加权之和最大为目标Max r(t),建立第t阶段内电-气综合能源系统恢复序列优化的目标函数:
式中,r(t)表示-气综合能源系统恢复序列优化的目标函数;t表示恢复阶段的序数;υ表示电力系统恢复权重系数;Ne表示电力系统中电节点的总数;Ng表示天然气系统中气节点的总数;i表示电节点的序数;j表示气节点的序数;表示第t阶段电节点i上的电力负荷;表示系统发生故障前电节点i上的电力负荷;表示第t阶段气节点j上的天然气负荷;表示系统发生故障前气节点j上的天然气负荷;
步骤3.2、建立第t阶段内电-气综合能源系统恢复序列优化的约束条件:
步骤3.2.1、建立第t阶段内的恢复资源约束:
式中,表示第t阶段电力线路l的恢复决策变量,当第t阶段电力线路l被修复时,否则l表示电力线路l的序数;Nl表示电力系统中电力线路的总数;表示第t阶段天然气管道p的恢复决策变量,为0-1整型变量,当第t阶段天然气管道被修复时,否则p表示天然气管道的序数;Np表示天然气系统中天然气管道的总数;Re表示电力系统维修设备数量,Rg表示天然气系统维修设备数量;
步骤3.2.2、建立第t阶段内的电力系统运行约束:
式中,表示第t阶段电力线路l的状态变量,为0-1整型变量,当第t阶段电力线路l为正常状态时,否则m表示电力系统恢复阶段的序数;表示第t阶段电节点i的状态变量,为0-1整型变量,当第t阶段电节点i为正常状态时,否则l∈i表示与电节点i相连的电力线路集合;表示第t阶段电节点i上的非燃气机组的出力;表示第t阶段电节点i上的燃气机组的出力;表示第t阶段电力线路l流过的潮流;l∈Ξ1(i)表示所有首端节点与电节点i相连的电力线路的集合;l∈Ξ2(i)表示所有末端节点与电节点i相连的电力线路的集合;表示第t阶段电节点i上的常规电负荷;表示第t阶段电节点i上的转气负荷;和分别表示第t阶段电力线路l首端节点和末端节点的电压相角;xbranch,l表示电力线路l的电抗;M表示一个很大的数;和分别表示电节点i上非燃气机组的最小技术出力和最大技术出力;和分别表示电节点i上燃气机组的最小技术出力和最大技术出力;表示电力线路l的传输容量;
步骤3.2.3、建立第t阶段内的天然气系统运行约束:
式中,表示第t阶段天然气管道p的状态变量,为0-1整型变量,当第t阶段天然气管道p为正常状态时,否则n表示天然气系统恢复阶段的序数;表示第t阶段气节点j的状态变量,为0-1整型变量,当第t阶段气节点j为正常状态时,否则p∈j表示与气节点j相连的天然气管道集合;表示第t阶段气节点j上的常规气源的出力;表示第t阶段气节点j上的电转气设备的出力;表示第t阶段天然气管道p流过的气流;p∈Ξ1(j)表示所有首端节点与气节点j相连的天然气管道的集合;p∈Ξ2(j)表示所有末端节点与气节点j相连的天然气管道的集合;表示第t阶段气节点j上的常规气负荷;表示第t阶段气节点j上的转电负荷;和分别表示第t阶段天然气管道p首端节点和末端节点的气压;Cp表示天然气管道p的传输系数;和分别表示气节点j上的常规气源的最小技术出力和最大技术出力;和分别表示气节点j上的电转气设备的最小技术出力和最大技术出力;和分别表示气节点j上的气压下限值和上限值;表示第t阶段气节点j上的气压;
步骤3.2.4、建立第t阶段内的电力-天然气相互依赖约束:
式中,ηgfu表示燃气机组的工作效率;ηptg表示电转气设备的工作效率;j∈Ξgfu(i)表示电节点i上的转气负荷向气节点j上的电转气设备供电;i∈Ξptg(j)表示气节点j上的转电负荷向电节点i上的燃气机组供气;
步骤3.3、基于反向贪心算法,对每一阶段的恢复决策进行计算;从当前待恢复电力线路、天然气管道集合中,找出使步骤3.1中优化目标值最小且满足步骤3.2中约束条件的电力线路、天然气管道作为本阶段的恢复对象;
步骤3.4、将步骤3.3中得到的恢复对象从当前待恢复电力线路、天然气管道集合中移除,判断当前待恢复电力线路、天然气管道集合是否为空集:
如果是,则说明所有待恢复电力线路、天然气管道已经被全部恢复,基于各阶段的恢复决策得到电力线路、天然气管道的恢复序列,进行步骤4;
否则,返回步骤3.3,更新当前待恢复电力线路、天然气管道集合,进行下一阶段的计算。
2.根据权利要求1所述的一种电-气综合能源系统恢复改善方法,其特征在于:所述步骤2,具体如下:
步骤2.2、确定待恢复电力线路集合Nl和天然气管道集合Np;
步骤2.3、确定电力系统维修设备数量Re和天然气系统维修设备数量Rg;
步骤2.4、确定电力系统所需的恢复阶段数Te、天然气系统所需的恢复阶段数Tg,得到电-气综合能源系统所需的总恢复阶段数:
T=max{Te,Tg} (4)
式中,T表示电-气综合能源系统所需的总恢复阶段数;Te表示电力系统所需的恢复阶段数;Tg表示天然气系统所需的恢复阶段数;max{Te,Tg}表示取Te和Tg两者中的最大值。
3.根据权利要求2所述的一种电-气综合能源系统恢复改善方法,其特征在于:所述步骤4中,具体如下:
步骤4.1、根据步骤3中得到的恢复序列,计算各阶段恢复的电力负荷,按照以下公式处理获得电力系统恢复力参数:
步骤4.2、根据步骤3中得到的恢复序列,计算各阶段恢复的天然气负荷,按照以下公式处理获得天然气系统恢复力参数:
步骤4.3、根据步骤4.1、步骤4.2中得到的电力系统恢复力参数和天然气系统恢复力参数,按照以下公式处理获得电-气综合能源系统恢复力参数:
Rres=υRres,e+(1-υ)Rres,g (34)
步骤4.4、根据计算得到电-气综合能源系统恢复力结果实施改善:
若电-气综合能源系统的恢复力参数值小于预设的恢复力阈值,则电-气综合能源系统恢复力水平低,按照以下措施进行改善:
对于电力系统,增加电力系统维修设备数量以增加每一阶段可恢复的电力线路数量,加粗电力线路以避免电力线路故障;
对于天然气系统,增加天然气系统维修设备数量以增加每一阶段可恢复的天然气管道数量,加厚天然气管道以避免天然气管道故障。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110036776.9A CN112861292B (zh) | 2021-01-12 | 2021-01-12 | 一种电-气综合能源系统恢复改善方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110036776.9A CN112861292B (zh) | 2021-01-12 | 2021-01-12 | 一种电-气综合能源系统恢复改善方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112861292A CN112861292A (zh) | 2021-05-28 |
CN112861292B true CN112861292B (zh) | 2022-08-05 |
Family
ID=76002886
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110036776.9A Active CN112861292B (zh) | 2021-01-12 | 2021-01-12 | 一种电-气综合能源系统恢复改善方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112861292B (zh) |
Family Cites Families (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108667007B (zh) * | 2018-04-16 | 2019-12-13 | 清华大学 | 计及电-气耦合系统约束的电压稳定裕度计算方法 |
CN110263387B (zh) * | 2019-05-29 | 2023-04-07 | 国网福建省电力有限公司 | 基于配电网、天然气网的能源系统运行装置 |
CN110263435B (zh) * | 2019-06-20 | 2021-02-09 | 燕山大学 | 基于电-气耦合综合能源系统的双层优化故障恢复方法 |
CN111444593B (zh) * | 2020-03-02 | 2022-05-03 | 浙江大学 | 一种电-气综合能源系统元件脆弱性改善方法 |
CN111222721B (zh) * | 2020-03-25 | 2022-08-12 | 广东电网有限责任公司广州供电局 | 一种提升电力系统恢复力的电-气网扩展规划方法 |
CN111276976B (zh) * | 2020-03-26 | 2021-05-11 | 华北电力大学 | 一种用于复杂能源互联系统极端事件的应急响应方法 |
CN111652441B (zh) * | 2020-06-04 | 2021-02-12 | 四川大学 | 考虑气电联合需求响应的气电综合能源系统配网优化方法 |
CN111987700B (zh) * | 2020-08-19 | 2022-07-08 | 国网综合能源服务集团有限公司 | 基于两阶段优化的综合能源系统型配电网故障恢复方法 |
CN112200695A (zh) * | 2020-10-30 | 2021-01-08 | 东南大学 | 一种城市级综合能源系统优化调度方法 |
-
2021
- 2021-01-12 CN CN202110036776.9A patent/CN112861292B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN112861292A (zh) | 2021-05-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107332234B (zh) | 一种考虑可再生能源间歇性的主动配电网多故障修复方法 | |
CN110503250B (zh) | 考虑电-热转移负荷转移量不确定性的综合能源规划方法 | |
CN111444593B (zh) | 一种电-气综合能源系统元件脆弱性改善方法 | |
CN109615248B (zh) | 一种考虑天然气系统影响的电力系统连锁故障估计方法 | |
CN112260261B (zh) | 一种电-气综合能源系统脆弱性评估与改善方法 | |
CN105529701B (zh) | 一种基于人工蜂群算法的含直流换流站电源恢复路径优化方法 | |
CN107658895B (zh) | 一种直流换相失败预测环节门槛值整定方法 | |
CN113742917B (zh) | 一种考虑多阶段恢复过程的综合能源系统韧性提升方法 | |
CN111931331B (zh) | 一种适用于配电网机组组合的分布鲁棒优化方法 | |
CN111371123B (zh) | 多类型电源汇集直流外送系统切机策略优化方法和装置 | |
CN115995790A (zh) | 一种配电网故障恢复方法、系统、设备及介质 | |
CN112928753A (zh) | 一种多能协同的配电网主动解列控制方法 | |
CN102545709B (zh) | 一种适用于太阳能发电装置的能量优化装置 | |
CN112861292B (zh) | 一种电-气综合能源系统恢复改善方法 | |
CN111581856B (zh) | 一种基于孤立电网的黑启动流程建模与仿真模拟方法 | |
CN105356508A (zh) | 基于psd-bpa的电网风电接入评价系统及其方法 | |
Saadeh et al. | Comparative performance analysis of HVDC and HVAC transmission systems in the presence of PV generation: a case study using the IEEE-5-bus network | |
WO2023103455A1 (zh) | 一种大规模海上风电场集电系统拓扑结构优化方法及系统 | |
CN112448404B (zh) | 一种电-气-热互联背景下的配电网可靠性增效计算方法 | |
CN113094854B (zh) | 一种电-气耦合系统骨干网架重构方法 | |
CN111786382A (zh) | 一种考虑加权潮流熵的配电网负荷恢复量计算方法 | |
CN115510618A (zh) | 计及元件重要度辨识的电气耦合系统两阶段韧性恢复方法 | |
CN115995815B (zh) | 一种基于多模块嵌套迭代的负荷故障恢复方法 | |
CN112329231B (zh) | 一种适用于交直流混联受端电网的电压弹性力评估方法 | |
CN113006895B (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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |