CN112186744B - 一种适用于含分布式电源配电网的供电恢复方法及应用 - Google Patents
一种适用于含分布式电源配电网的供电恢复方法及应用 Download PDFInfo
- Publication number
- CN112186744B CN112186744B CN202010975468.8A CN202010975468A CN112186744B CN 112186744 B CN112186744 B CN 112186744B CN 202010975468 A CN202010975468 A CN 202010975468A CN 112186744 B CN112186744 B CN 112186744B
- Authority
- CN
- China
- Prior art keywords
- power
- period
- node
- distribution network
- power supply
- 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
- 238000009826 distribution Methods 0.000 title claims abstract description 129
- 238000000034 method Methods 0.000 title claims abstract description 98
- 238000011084 recovery Methods 0.000 title claims abstract description 64
- 238000005457 optimization Methods 0.000 claims abstract description 19
- 238000004458 analytical method Methods 0.000 claims abstract description 7
- 238000004146 energy storage Methods 0.000 claims description 31
- 230000008569 process Effects 0.000 claims description 18
- 238000004519 manufacturing process Methods 0.000 claims description 17
- 238000007599 discharging Methods 0.000 claims description 9
- 238000005094 computer simulation Methods 0.000 claims description 6
- 230000000694 effects Effects 0.000 claims description 4
- 230000008859 change Effects 0.000 claims description 3
- 238000010276 construction Methods 0.000 claims description 3
- 238000010438 heat treatment Methods 0.000 claims description 3
- 238000005065 mining Methods 0.000 claims description 3
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 3
- 238000004088 simulation Methods 0.000 description 8
- 238000010586 diagram Methods 0.000 description 5
- 238000010248 power generation Methods 0.000 description 4
- 238000011156 evaluation Methods 0.000 description 3
- 238000004836 empirical method Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 1
- 230000000052 comparative effect Effects 0.000 description 1
- 230000002708 enhancing effect Effects 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 230000008092 positive effect Effects 0.000 description 1
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
-
- 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/12—Circuit arrangements for ac mains or ac distribution networks for adjusting voltage in ac networks by changing a characteristic of the network load
- H02J3/14—Circuit arrangements for ac mains or ac distribution networks for adjusting voltage in ac networks by changing a characteristic of the network load by switching loads on to, or off from, network, e.g. progressively balanced loading
- H02J3/144—Demand-response operation of the power transmission or distribution network
-
- 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/388—Islanding, i.e. disconnection of local power supply from the network
-
- 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
-
- 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/50—Controlling the sharing of the out-of-phase component
-
- 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/10—The dispersed energy generation being of fossil origin, e.g. diesel 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
- 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
- Y02B—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO BUILDINGS, e.g. HOUSING, HOUSE APPLIANCES OR RELATED END-USER APPLICATIONS
- Y02B70/00—Technologies for an efficient end-user side electric power management and consumption
- Y02B70/30—Systems integrating technologies related to power network operation and communication or information technologies for improving the carbon footprint of the management of residential or tertiary loads, i.e. smart grids as climate change mitigation technology in the buildings sector, including also the last stages of power distribution and the control, monitoring or operating management systems at local level
- Y02B70/3225—Demand response systems, e.g. load shedding, peak shaving
-
- 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
- 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
- Y04S20/00—Management or operation of end-user stationary applications or the last stages of power distribution; Controlling, monitoring or operating thereof
- Y04S20/20—End-user application control systems
- Y04S20/222—Demand response systems, e.g. load shedding, peak shaving
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明公开一种适用于含分布式电源配电网的供电恢复方法。所述配电网供电恢复方法包括:根据产业类型,对配电网内的负荷用户进行分类;通过间接分析法计算每个种类负荷用户的单小时用户停电损失值;建立以总用户停电损失值最小为目标的多时段配电网供电恢复最优化问题模型;求解所述多时段配电网供电恢复最优化问题,得到多时段含分布式电源配电网负荷供电恢复方案。本发明能够解决变电站故障仅分布式电源可用时的配电网供电恢复问题,协调控制分布式电源的输出,能达到最小化减少用户损失、提高配网系统的恢复能力和韧性的目的。
Description
技术领域
本发明涉及配电网供电恢复领域,特别是涉及一种适用于含分布式电源配电网的供电恢复方法。
背景技术
近些年,随着全球气候的变化,极端天气事件发生次数增加,这使得提高配电网的韧性变得越来越重要。为此,在由极端天气事件导致的配电网停电事故发生后,需要采用适当的配电网供电恢复方法来快速恢复重要负载,提升配电网韧性。分布式电源(Distributed Energy Source,DER)和微电网在配电网中的引入也提高了系统的韧性,并为更新颖和有效的配电网供电恢复方法的出现提供了机遇。目前,已有不少文献对含DER的配电网供电恢复方法进行了研究,然而现有的研究很少将用户停电损失(CustomerInterruption Cost,CIC)与供电恢复方法直接联系起来,且这些研究多采用经验法直接选择成本系数来简单描述用户的负荷断电成本,对用户停电损失的评估缺乏系统性和准确性。
发明内容
鉴于以上所述现有技术的缺点,本发明的目的在于提供一种适用于含分布式电源配电网的供电恢复方法,用于解决现有技术中不能将用户停电经济损失与供电恢复方法准确联系,无法得到使用户停电损失最小的配电网的供电恢复方案的问题。
为实现上述目的及其他相关目的,本发明提供一种适用于含分布式电源配电网的供电恢复方法,所述配电网供电恢复方法包括如下步骤:
根据产业类型,对含分布式电源配电网内的负荷用户进行分类;
通过间接分析法获取每个种类负荷用户的单小时用户停电损失值;
建立以总用户停电损失值最小为目标的多时段含分布式电源配电网供电恢复模型;
对所述多时段含分布式电源配电网供电恢复模型进行求解,得到多时段含分布式电源配电网负荷供电恢复方案。
进一步的,根据产业类型,对含分布式电源配电网内的负荷用户进行分类的方法如下:
根据是否从事生产活动将配电网负荷用户分为住宅用户和产业用户两类;
根据产业用户从事的产业类型,将产业用户进一步分为第一产业用户、第二产业用户、第三产业用户;其中,所述第一产业用户为从事农业、林业、牧业、渔业的用户;所述第二产业用户为从事矿业、制造业、建筑业以及电力、热力、燃气及水生产和供应业的用户;所述第三产业用户为从事服务业,即第一产业与第二产业之外行业的用户。
进一步的,所述通过间接分析法获取每个种类负荷用户的单小时用户停电损失值,具体如下:
住宅用户的单小时停电损失值可通过下式计算:
其中,CICresi表示住宅用户的单小时用户停电损失值;w表示该地区平均时薪;Pmax表示住宅用户的用电功率峰值;
产业用户的停电损失可通过下式获取:
其中,CICindu表示产业用户的单小时用户停电损失值;VA表示产业用户产值的平均年增加值;Pfp表示产业用户全产能生产时的平均用电需求;TP表示产业用户一年内平均用于生产的时间;cva为比例系数,表示产业用户平均总停电损失与平均断电停产损失之比。
进一步的,所述建立以总用户停电损失值最小为目标的多时段含分布式电源配电网供电恢复模型的步骤包括:
建立以总用户停电损失值最小为目标的多时段含分布式电源配电网供电恢复目标函数;
确定包含系统建模约束和系统运行约束的多时段含分布式电源配电网供电恢复最优化问题约束条件;
对所述多时段含分布式电源配电网供电恢复最优化问题约束条件中的非线性部分进行线性化处理。
进一步的,所述以总用户停电损失值最小为目标的多时段含分布式电源配电网供电恢复目标函数为:
式中,CICi,t表示t时段节点i处负荷的单小时用户停电损失;Pi D表示t时段节点i处负荷的有功功率需求;ri,t为决策变量,表示t时段节点i处负荷的恢复状态,ri,t为1表示此时节点i处负荷已恢复供电,ri,t为0表示此时节点i处负荷未恢复供电;Tint表示每个断电恢复时段的时长;为配电网所有节点的集合;/>为所有停电恢复时段的集合;
所述系统建模约束包括分布式电源输出约束和辐射状拓扑约束;
其中,所述分布式电源输出约束如下:
式中,表示t时段节点i处所连所有分布式电源输出有功功率之和,/>表示t时段节点i处所连所有分布式电源输出无功功率之和;/>表示t时段节点i处所连可控分布式电源输出有功功率,/>表示t时段节点i处所连可控分布式电源输出无功功率;/>表示t时段节点i处所连不可控分布式电源输出有功功率,/>表示t时段节点i处所连不可控分布式电源输出无功功率;/>和/>分别表示t时段节点i处所连储能系统的放电功率和充电功率;Sd为配电网连接可控分布式电源的节点集合;Snd为配电网连接不可控分布式电源的节点集合;Sε为配电网连接储能系统的节点集合;
所述节点i处所连可控分布式电源,包括但不限于微型燃气轮机、柴油发电机等,均满足下列输出约束:
式中,Pi dDGmin和Pi dDGmax分别表示t时段节点i处所连可控分布式电源输出有功功率下限和上限;和/>分别表示t时段节点i处所连可控分布式电源输出无功功率下限和上限;/>表示节点i处所连可控分布式电源输出有功功率的变化率的上限。
所述节点i处所连不可控分布式电源,包括但不限于风力发电、光伏发电等,均满足下列输出约束:
式中,表示t时段节点i处所连不可控分布式电源输出有功功率的预测值,γi表示节点i处所连不可控分布式电源的最大功率因数角;
所述节点i处所连储能系统,满足下列输出约束:
式中,和/>为决策变量,分别表示t时段节点i处所连储能系统的充电状态和放电状态,若/>为1,表示t时段节点i处所连ESS正在充电;若/>为0,表示t时段节点i处所连ESS未在充电;若/>为1,表示t时段节点i处所连ESS正在放电;若/>为0,表示t时段节点i处所连ESS未在放电;Pi dischmax和Pi chmax分别表示节点i处所连储能系统的最大放电功率和最大充电功率;/>表示节点i处所连储能系统的荷电状态的初始值,/>和/>分别表示节点i处所连储能系统的荷电状态的下限和上限;Crated,i表示节点i处所连储能系统的额定能量容量;/>和/>分别表示节点i处所连储能系统的充电效率和放电效率;
所述辐射状拓扑约束如下:
式中,bij,t为辅助变量,表示t时段支路ij两端节点的父子关系,若节点i是节点j的父节点,bij,t=0;反之,bij,t=1;αij,t为决策变量,表示t时段支路ij的开关状态,αij,t为0表示t时段支路ij断开;αij,t为1表示t时段支路ij连通;Fij,t表示t时段流经支路ij的虚拟潮流;Dj表示t时段节点j处的虚拟负荷需求;M表示一个无穷大的正数;ε表示配电网所有支路的集合;表示配电网所有根节点的集合;
所述系统运行约束包括配电网潮流约束、安全运行约束以及系统热备用容量约束;
其中,所述配电网潮流约束如下所示:
式中,Pij,t和Qij,t分别表示t时段流经支路ij的有功功率潮流和无功功率潮流;表示t时段节点j处负荷的无功功率需求;Rij和Xij分别表示配电网支路ij的电阻和电抗;表示t时段节点i处电压幅值的平方;/>表示t时段流经支路ij的电流幅值的平方;
所述安全运行约束如下所示
式中,表示t时段流经支路ij的电流幅值的上限;/>分别表示t时段节点i处电压幅值的下限、上限;/>表示t时段流经支路ij的有功功率潮流的上限;
所述系统热备用容量约束如下所示:
式中,μ为系统热备用率;表示t时段节点i处所连所有分布式电源的输出上限。
进一步的,对所述多时段含分布式电源配电网供电恢复最优化问题约束条件中的非线性部分进行线性化处理的步骤包括:
对所述多时段含分布式电源配电网供电恢复最优化问题约束条件中的功率潮流平方项进行分段近似线性化;具体如下:
所述多时段配电网供电恢复最优化问题约束条件中的功率潮流平方项包括有功功率潮流平方项和无功功率潮流平方项/>所述有功功率潮流平方项/>的分段线性化近似过程如下:
式中,表示所述有功功率潮流平方项/>的近似值;/>和/>为线性化过程的辅助积分变量;a为计数变量;A为分段近似的段数;
无功功率潮流平方项的近似值/>可通过相同过程计算得到。
对所述多时段含分布式电源配电网供电恢复最优化问题约束条件中的电压幅值平方项以常数标准值近似线性化;如下所示:
式中,Unsqr表示电压幅值平方项的常数标准值。
进一步的,应用基于MATLAB平台的Gurobi求解器提供的混合整数线性规划(mixed-integer linear programming,MILP)求解方法对所述多时段含分布式电源配电网供电恢复模型进行求解,得到的多时段含分布式电源配电网负荷供电恢复方案为以各个恢复时段负荷恢复状态、馈线开关状态和分布式电源输出为决策变量的动态孤岛划分方案和分布式电源输出分配方案。
一种适用于含分布式电源配电网的供电恢复方法的应用,所述适用于含分布式电源配电网的供电恢复方法能应用在极端天气下,变电站无法正常工作,配电网仅以分布式电源作为供电恢复电能来源的情况下。
本发明的优点和积极效果是:
本发明的方法与现有的基于加权恢复负荷量和基于加权恢复供电时间的配电网供电恢复技术相比,本发明的方法能将用户停电损失与供电恢复方法直接联系起来,得到使用户停电损失最小的供电恢复结果,降低用户因意外停电产生的经济损失;针对现有基于停电成本的配电网供电恢复技术中多采用经验法直接选择成本系数来简单描述用户的负荷断电成本,对用户停电损失的评估缺乏系统性和准确性的问题,本发明的方法利用权威统计机构公开的区域生产经济数据,采用宏观经济学方法对用户停电损失进行评估,具有较高的真实性和可信度,可解决用户停电损失的评估不准确的问题。
附图说明
图1为本申请实施例提供的一种适用于含分布式电源配电网的供电恢复方法的工作流程图;
图2为本申请实施例提供的含分布式电源的33节点配电网系统图;
图3为本申请实施例提供的风机和光伏24小时出力折线图;
图4为本申请实施例提供的使用本发明供电恢复方法时各种分布式电源在各个恢复时刻的有功出力仿真结果柱状图;
图5为本申请实施例提供的本发明供电恢复方法与基于时间模型的供电恢复方法对比的每时刻恢复供电负荷功率仿真结果柱状图。
具体实施方式
以下通过特定的具体实例说明本发明的实施方式,本领域技术人员可由本说明书所揭露的内容轻易地了解本发明的其他优点与功效。本发明还可以通过另外不同的具体实施方式加以实施或应用,本说明书中的各项细节也可以基于不同观点与应用,在没有背离本发明的精神下进行各种修饰或改变。需说明的是,在不冲突的情况下,以下实施例及实施例中的特征可以相互组合。
需要说明的是,以下实施例中所提供的图示仅以示意方式说明本发明的基本构想,遂图示中仅显示与本发明中有关的组件而非按照实际实施时的组件数目、形状及尺寸绘制,其实际实施时各组件的型态、数量及比例可为一种随意的改变,且其组件布局型态也可能更为复杂。
图1为本申请实施例提供的一种适用于含分布式电源配电网的供电恢复方法的工作流程图;本发明提供一种适用于含分布式电源配电网的供电恢复方法,所述配电网供电恢复方法包括如下步骤:
S1、根据产业类型,对含分布式电源配电网内的负荷用户进行分类;
S2、通过间接分析法获取每个种类负荷用户的单小时用户停电损失值;
S3、建立以总用户停电损失值最小为目标的多时段含分布式电源配电网供电恢复模型;
S4、对所述多时段含分布式电源配电网供电恢复模型进行求解,得到多时段含分布式电源配电网负荷供电恢复方案。
具体的,根据产业类型,对含分布式电源配电网内的负荷用户进行分类的方法如下:
根据是否从事生产活动将配电网负荷用户分为住宅用户和产业用户两类;
根据产业用户从事的产业类型,将产业用户进一步分为第一产业用户、第二产业用户、第三产业用户;其中,所述第一产业用户为从事农业、林业、牧业、渔业的用户;所述第二产业用户为从事矿业、制造业、建筑业以及电力、热力、燃气及水生产和供应业的用户;所述第三产业用户为从事服务业,即第一产业与第二产业之外行业的用户。
其中,所述通过间接分析法获取每个种类负荷用户的单小时用户停电损失值,具体如下:
住宅用户的单小时停电损失值可通过下式计算:
其中,CICresi表示住宅用户的单小时用户停电损失值;w表示该地区平均时薪;Pmax表示住宅用户的用电功率峰值;
产业用户的停电损失可通过下式获取:
其中,CICindu表示产业用户的单小时用户停电损失值;VA表示产业用户产值的平均年增加值;Pfp表示产业用户全产能生产时的平均用电需求;TP表示产业用户一年内平均用于生产的时间;cva为比例系数,表示产业用户平均总停电损失与平均断电停产损失之比,其中,产业用户平均总停电损失、平均断电停产损失均可通过当地的产业相关数据获取。
此外,所述建立以总用户停电损失值最小为目标的多时段含分布式电源配电网供电恢复模型的步骤包括:
建立以总用户停电损失值最小为目标的多时段含分布式电源配电网供电恢复目标函数;
确定包含系统建模约束和系统运行约束的多时段含分布式电源配电网供电恢复最优化问题约束条件;
对所述多时段含分布式电源配电网供电恢复最优化问题约束条件中的非线性部分进行线性化处理。
进一步的,所述以总用户停电损失值最小为目标的多时段含分布式电源配电网供电恢复目标函数为:
式中,CICi,t表示t时段节点i处负荷的单小时用户停电损失;Pi D表示t时段节点i处负荷的有功功率需求;ri,t为决策变量,表示t时段节点i处负荷的恢复状态,ri,t为1表示此时节点i处负荷已恢复供电,ri,t为0表示此时节点i处负荷未恢复供电;Tint表示每个断电恢复时段的时长;为配电网所有节点的集合;/>为所有停电恢复时段的集合;
所述系统建模约束包括分布式电源输出约束和辐射状拓扑约束;
其中,所述分布式电源输出约束如下:
式中,表示t时段节点i处所连所有分布式电源输出有功功率之和,/>表示t时段节点i处所连所有分布式电源输出无功功率之和;/>表示t时段节点i处所连可控分布式电源输出有功功率,/>表示t时段节点i处所连可控分布式电源输出无功功率;/>表示t时段节点i处所连不可控分布式电源输出有功功率,/>表示t时段节点i处所连不可控分布式电源输出无功功率;/>和/>分别表示t时段节点i处所连储能系统的放电功率和充电功率;Sd为配电网连接可控分布式电源的节点集合;Snd为配电网连接不可控分布式电源的节点集合;Sε为配电网连接储能系统的节点集合;
所述节点i处所连可控分布式电源,包括但不限于微型燃气轮机、柴油发电机等,均满足下列输出约束:
式中,Pi dDGmin和Pi dDGmax分别表示t时段节点i处所连可控分布式电源输出有功功率下限和上限;和/>分别表示t时段节点i处所连可控分布式电源输出无功功率下限和上限;/>表示节点i处所连可控分布式电源输出有功功率的变化率的上限。
所述节点i处所连不可控分布式电源,包括但不限于风力发电、光伏发电等,均满足下列输出约束:
式中,表示t时段节点i处所连不可控分布式电源输出有功功率的预测值,γi表示节点i处所连不可控分布式电源的最大功率因数角;
所述节点i处所连储能系统,满足下列输出约束:
式中,和/>为决策变量,分别表示t时段节点i处所连储能系统的充电状态和放电状态,若/>为1,表示t时段节点i处所连ESS正在充电;若/>为0,表示t时段节点i处所连ESS未在充电;若/>为1,表示t时段节点i处所连ESS正在放电;若/>为0,表示t时段节点i处所连ESS未在放电;Pi dischmax和Pi chmax分别表示节点i处所连储能系统的最大放电功率和最大充电功率;/>表示节点i处所连储能系统的荷电状态的初始值,/>和/>分别表示节点i处所连储能系统的荷电状态的下限和上限;Crated,i表示节点i处所连储能系统的额定能量容量;/>和/>分别表示节点i处所连储能系统的充电效率和放电效率;
所述辐射状拓扑约束如下:
式中,bij,t为辅助变量,表示t时段支路ij两端节点的父子关系,若节点i是节点j的父节点,bij,t=0;反之,bij,t=1;αij,t为决策变量,表示t时段支路ij的开关状态,αij,t为0表示t时段支路ij断开;αij,t为1表示t时段支路ij连通;Fij,t表示t时段流经支路ij的虚拟潮流;Dj表示t时段节点j处的虚拟负荷需求;M表示一个无穷大的正数;ε表示配电网所有支路的集合;表示配电网所有根节点的集合;
所述系统运行约束包括配电网潮流约束、安全运行约束以及系统热备用容量约束;
其中,所述配电网潮流约束如下所示:
式中,Pij,t和Qij,t分别表示t时段流经支路ij的有功功率潮流和无功功率潮流;表示t时段节点j处负荷的无功功率需求;Rij和Xij分别表示配电网支路ij的电阻和电抗;表示t时段节点i处电压幅值的平方;/>表示t时段流经支路ij的电流幅值的平方;
所述安全运行约束如下所示
式中,表示t时段流经支路ij的电流幅值的上限;/>分别表示t时段节点i处电压幅值的下限、上限;/>表示t时段流经支路ij的有功功率潮流的上限;
所述系统热备用容量约束如下所示:
式中,μ为系统热备用率;表示t时段节点i处所连所有分布式电源的输出上限。
进一步的,对所述多时段含分布式电源配电网供电恢复最优化问题约束条件中的非线性部分进行线性化处理的步骤包括:
对所述多时段含分布式电源配电网供电恢复最优化问题约束条件中的功率潮流平方项进行分段近似线性化;具体如下:
所述多时段配电网供电恢复最优化问题约束条件中的功率潮流平方项包括有功功率潮流平方项和无功功率潮流平方项/>所述有功功率潮流平方项/>的分段线性化近似过程如下:
式中,表示所述有功功率潮流平方项/>的近似值;/>和/>为线性化过程的辅助积分变量;a为计数变量;A为分段近似的段数;
无功功率潮流平方项的近似值/>可通过相同过程计算得到。
对所述多时段含分布式电源配电网供电恢复最优化问题约束条件中的电压幅值平方项以常数标准值近似线性化;如下所示:
式中,Unsqr表示电压幅值平方项的常数标准值。
最后,应用基于MATLAB平台的Gurobi求解器提供的混合整数线性规划(mixed-integer linear programming,MILP)求解方法对所述多时段含分布式电源配电网供电恢复模型进行求解,得到的多时段含分布式电源配电网负荷供电恢复方案为以各个恢复时段负荷恢复状态、馈线开关状态和分布式电源输出为决策变量的动态孤岛划分方案和分布式电源输出分配方案。
一种适用于含分布式电源配电网的供电恢复方法的应用,所述适用于含分布式电源配电网的供电恢复方法能应用在极端天气下,变电站无法正常工作,配电网仅以分布式电源作为供电恢复电能来源的情况下。
作为举例,在本实施例中,为了验证本发明方案的有效性,进行如下仿真实验:
请参阅图2、图3,图2为本申请实施例提供的含分布式电源的33节点配电网系统图,图3为本申请实施例提供的风机和光伏24小时出力折线图;本发明使用了含分布式电源的33节点配电网系统,其中包含4个分布式电源。每个分布式电源包含一个微型燃气轮机、一个风机、一个光伏机组和一个储能系统。每个分布式电源中相同种类的电源具有相同的参数。微型燃气轮机的峰值输出功率为600kw,风机和光伏的峰值输出功率分别为500kW和300kw,且风机和光伏的最小允许功率因数均为0.9,其各自的24小时出力如图3所示。储能系统的参数如表一所示。
表一储能系统参数
根据中国某市统计局公布的地区工资和生产数据,通过计算得出每种用户的单小时用户停电损失值,如表二所示。为了更贴合实际,设置工作时间为9:00-17:00,这段时间以外的用户停电损失值更小。所有用户负载都设定为可中断负载。配电网供电恢复时间段的开始时间是上午11点,总恢复时间是10小时,每个时间段时长是1小时。
表二每种用户的单小时用户停电损失值(单位:元)
用户类型 | 工作时间 | 非工作时间 |
住宅用户 | 3.427 | 2.742 |
第一产业用户 | 15.94 | 12.75 |
第二产业用户 | 26.11 | 20.89 |
第三产业用户 | 84.37 | 67.50 |
请参阅图4。图4为本申请实施例提供的使用本发明供电恢复方法时各种分布式电源在各个恢复时刻的有功出力仿真结果柱状图。从图4可以看出,本发明提出的配电网供电恢复方法能够协调微型燃气轮机和储能系统的输出功率,以应对风机和光伏机组的不确定性,保持用于负载恢复的输出功率稳定。
为了进一步验证本发明所提配电网供电恢复方法的有效性,进行了与目标函数为最大化加权供电时间的基于时间模型的配电网供电恢复方法的对比仿真,并与本发明中的配电网供电恢复方法在相同约束条件下进行求解。为了便于比较,设置住宅、第一产业用户、第二产业用户和第三产业用户的权重分别为0.1、1、10、100。两种供电恢复方法的部分仿真结果比较如表三所示。
表三两种供电恢复方法的部分仿真结果
由表三可以明显看出,与基于时间模型的相比供电恢复方法,基于用户停电损失的供电恢复方法具有更大的总恢复负荷电能,更小的总总用户停电损失值,更高的可再生能源利用效率,在增强配电系统的弹性方面更具优势。
请参阅图5,图5为本申请实施例提供的本发明供电恢复方法与基于时间模型的供电恢复方法对比的每时刻恢复供电负荷功率仿真结果柱状图。从图,5中可以看出,尽管两种供电恢复方法都能保证产业用户的供电,但本发明中的供电恢复方法能恢复更多的住宅用户负荷,并一直保持稳定较高的电能供应水平。
上述实施例仅例示性说明本发明的原理及其功效,而非用于限制本发明。任何熟悉此技术的人士皆可在不违背本发明的精神及范畴下,对上述实施例进行修饰或改变。因此,举凡所属技术领域中具有通常知识者在未脱离本发明所揭示的精神与技术思想下所完成的一切等效修饰或改变,仍应由本发明的权利要求所涵盖。
Claims (6)
1.一种适用于含分布式电源配电网的供电恢复方法,其特征在于,所述配电网供电恢复方法包括如下步骤:
根据产业类型,对含分布式电源配电网内的负荷用户进行分类;
通过间接分析法获取每个种类负荷用户的单小时用户停电损失值;
建立以总用户停电损失值最小为目标的多时段含分布式电源配电网供电恢复模型;
对所述多时段含分布式电源配电网供电恢复模型进行求解,得到多时段含分布式电源配电网负荷供电恢复方案;
所述建立以总用户停电损失值最小为目标的多时段含分布式电源配电网供电恢复模型的步骤包括:
建立以总用户停电损失值最小为目标的多时段含分布式电源配电网供电恢复目标函数;
确定包含系统建模约束和系统运行约束的多时段含分布式电源配电网供电恢复最优化问题约束条件;
对所述多时段含分布式电源配电网供电恢复最优化问题约束条件中的非线性部分进行线性化处理;
所述以总用户停电损失值最小为目标的多时段含分布式电源配电网供电恢复目标函数为:
式中,CICi,t表示t时段节点i处负荷的单小时用户停电损失;Pi D表示t时段节点i处负荷的有功功率需求;ri,t为决策变量,表示t时段节点i处负荷的恢复状态,ri,t为1表示此时节点i处负荷已恢复供电,ri,t为0表示此时节点i处负荷未恢复供电;Tint表示每个断电恢复时段的时长;为配电网所有节点的集合;/>为所有停电恢复时段的集合;
所述系统建模约束包括分布式电源输出约束和辐射状拓扑约束;
其中,所述分布式电源输出约束如下:
式中,表示t时段节点i处所连所有分布式电源输出有功功率之和,/>表示t时段节点i处所连所有分布式电源输出无功功率之和;/>表示t时段节点i处所连可控分布式电源输出有功功率,/>表示t时段节点i处所连可控分布式电源输出无功功率;/>表示t时段节点i处所连不可控分布式电源输出有功功率,/>表示t时段节点i处所连不可控分布式电源输出无功功率;/>和/>分别表示t时段节点i处所连储能系统的放电功率和充电功率;/>为配电网连接可控分布式电源的节点集合;/>为配电网连接不可控分布式电源的节点集合;/>为配电网连接储能系统的节点集合;
所述节点i处所连可控分布式电源,均满足下列输出约束:
式中,Pi dDGmin和Pi dDGmax分别表示t时段节点i处所连可控分布式电源输出有功功率下限和上限;Qi dDGmin和Qi dDGmax分别表示t时段节点i处所连可控分布式电源输出无功功率下限和上限;Ri dDGmax表示节点i处所连可控分布式电源输出有功功率的变化率的上限;
所述节点i处所连不可控分布式电源,均满足下列输出约束:
式中,表示t时段节点i处所连不可控分布式电源输出有功功率的预测值,γi表示节点i处所连不可控分布式电源的最大功率因数角;
所述节点i处所连储能系统,满足下列输出约束:
式中,和/>为决策变量,分别表示t时段节点i处所连储能系统的充电状态和放电状态,若/>为1,表示t时段节点i处所连ESS正在充电;若/>为0,表示t时段节点i处所连ESS未在充电;若/>为1,表示t时段节点i处所连ESS正在放电;若/>为0,表示t时段节点i处所连ESS未在放电;Pi dischmax和Pi chmax分别表示节点i处所连储能系统的最大放电功率和最大充电功率;/>表示节点i处所连储能系统的荷电状态的初始值,/>和/>分别表示节点i处所连储能系统的荷电状态的下限和上限;Crated,i表示节点i处所连储能系统的额定能量容量;/>和/>分别表示节点i处所连储能系统的充电效率和放电效率;
所述辐射状拓扑约束如下:
式中,bij,t为辅助变量,表示t时段支路ij两端节点的父子关系,若节点i是节点j的父节点,bij,t=0;反之,bij,t=1;αij,t为决策变量,表示t时段支路ij的开关状态,αij,t为0表示t时段支路ij断开;αij,t为1表示t时段支路ij连通;Fij,t表示t时段流经支路ij的虚拟潮流;Dj表示t时段节点j处的虚拟负荷需求;M表示一个无穷大的正数;表示配电网所有支路的集合;/>表示配电网所有根节点的集合;
所述系统运行约束包括配电网潮流约束、安全运行约束以及系统热备用容量约束;
其中,所述配电网潮流约束如下所示:
式中,Pij,t和Qij,t分别表示t时段流经支路ij的有功功率潮流和无功功率潮流;表示t时段节点j处负荷的无功功率需求;Rij和Xij分别表示配电网支路ij的电阻和电抗;/>表示t时段节点i处电压幅值的平方;/>表示t时段流经支路ij的电流幅值的平方;
所述安全运行约束如下所示
式中,表示t时段流经支路ij的电流幅值的上限;/>分别表示t时段节点i处电压幅值的下限、上限;/>表示t时段流经支路ij的有功功率潮流的上限;
所述系统热备用容量约束如下所示:
式中,μ为系统热备用率;表示t时段节点i处所连所有分布式电源的输出上限。
2.根据权利要求1所述的一种适用于含分布式电源配电网的供电恢复方法,其特征在于,根据产业类型,对含分布式电源配电网内的负荷用户进行分类的方法如下:
根据是否从事生产活动将配电网负荷用户分为住宅用户和产业用户两类;
根据产业用户从事的产业类型,将产业用户进一步分为第一产业用户、第二产业用户、第三产业用户;其中,所述第一产业用户为从事农业、林业、牧业、渔业的用户;所述第二产业用户为从事矿业、制造业、建筑业以及电力、热力、燃气及水生产和供应业的用户;所述第三产业用户为从事服务业,即第一产业与第二产业之外行业的用户。
3.根据权利要求1所述的一种适用于含分布式电源配电网的供电恢复方法,其特征在于,所述通过间接分析法获取每个种类负荷用户的单小时用户停电损失值,具体如下:
住宅用户的单小时停电损失值可通过下式计算:
其中,CICresi表示住宅用户的单小时用户停电损失值;w表示该住宅地区平均时薪;Pmax表示住宅用户的用电功率峰值;
产业用户的停电损失可通过下式获取:
其中,CICindu表示产业用户的单小时用户停电损失值;VA表示产业用户产值的平均年增加值;Pfp表示产业用户全产能生产时的平均用电需求;TP表示产业用户一年内平均用于生产的时间;cva为比例系数,表示产业用户平均总停电损失与平均断电停产损失之比。
4.根据权利要求1所述的一种适用于含分布式电源配电网的供电恢复方法,其特征在于,对所述多时段含分布式电源配电网供电恢复最优化问题约束条件中的非线性部分进行线性化处理的步骤包括:
对所述多时段含分布式电源配电网供电恢复最优化问题约束条件中的功率潮流平方项进行分段近似线性化;
对所述多时段含分布式电源配电网供电恢复最优化问题约束条件中的电压幅值平方项以常数标准值近似线性化。
5.根据权利要求1所述的一种适用于含分布式电源配电网的供电恢复方法,其特征在于,对所述多时段含分布式电源配电网供电恢复模型进行求解,得到的多时段含分布式电源配电网负荷供电恢复方案为以各个恢复时段负荷恢复状态、馈线开关状态和分布式电源输出为决策变量的动态孤岛划分方案和分布式电源输出分配方案。
6.权利要求1-5任一项所述的一种适用于含分布式电源配电网的供电恢复方法的应用,其特征在于:
所述适用于含分布式电源配电网的供电恢复方法能应用在极端天气下,变电站无法正常工作,配电网仅以分布式电源作为供电恢复电能来源的情况下。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010975468.8A CN112186744B (zh) | 2020-09-16 | 2020-09-16 | 一种适用于含分布式电源配电网的供电恢复方法及应用 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010975468.8A CN112186744B (zh) | 2020-09-16 | 2020-09-16 | 一种适用于含分布式电源配电网的供电恢复方法及应用 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112186744A CN112186744A (zh) | 2021-01-05 |
CN112186744B true CN112186744B (zh) | 2024-03-08 |
Family
ID=73921486
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010975468.8A Active CN112186744B (zh) | 2020-09-16 | 2020-09-16 | 一种适用于含分布式电源配电网的供电恢复方法及应用 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112186744B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113131468B (zh) * | 2021-04-16 | 2023-04-07 | 东南大学 | 一种考虑用户停电损失的多目标配电网供电恢复方法 |
CN113705964B (zh) * | 2021-07-16 | 2024-01-02 | 国网上海市电力公司 | 一种配电网韧性恢复灾前预案制定方法及装置 |
CN114069618B (zh) * | 2021-11-15 | 2024-08-16 | 国网江苏省电力有限公司常州供电分公司 | 基于最小总停电损失的配电网供电恢复方法 |
CN115619080A (zh) * | 2022-11-11 | 2023-01-17 | 国网吉林省电力有限公司长春供电公司 | 一种考虑用户停电损失的配电网灾后应急抢修方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104809524A (zh) * | 2015-05-04 | 2015-07-29 | 国家电网公司 | 基于拟最小路法的智能配电网供电分区恢复与优化方法 |
CN109802387A (zh) * | 2019-03-11 | 2019-05-24 | 燕山大学 | 一种含微网的弹性配电网多级供电恢复方法 |
CN110350508A (zh) * | 2019-05-16 | 2019-10-18 | 东南大学 | 一种同时考虑重构与孤岛划分的主动配电网故障恢复统一模型的方法 |
CN111293683A (zh) * | 2020-02-13 | 2020-06-16 | 东方电子股份有限公司 | 一种兼顾安全性与经济性的配网自愈优化方法 |
CN111478358A (zh) * | 2019-11-06 | 2020-07-31 | 南瑞集团有限公司 | 一种计及分布式电源不确定性的配电网鲁棒恢复决策方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105140917B (zh) * | 2015-09-06 | 2017-05-10 | 清华大学 | 适用于不确定性环境下的主动配电网鲁棒恢复控制方法 |
-
2020
- 2020-09-16 CN CN202010975468.8A patent/CN112186744B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104809524A (zh) * | 2015-05-04 | 2015-07-29 | 国家电网公司 | 基于拟最小路法的智能配电网供电分区恢复与优化方法 |
CN109802387A (zh) * | 2019-03-11 | 2019-05-24 | 燕山大学 | 一种含微网的弹性配电网多级供电恢复方法 |
CN110350508A (zh) * | 2019-05-16 | 2019-10-18 | 东南大学 | 一种同时考虑重构与孤岛划分的主动配电网故障恢复统一模型的方法 |
CN111478358A (zh) * | 2019-11-06 | 2020-07-31 | 南瑞集团有限公司 | 一种计及分布式电源不确定性的配电网鲁棒恢复决策方法 |
CN111293683A (zh) * | 2020-02-13 | 2020-06-16 | 东方电子股份有限公司 | 一种兼顾安全性与经济性的配网自愈优化方法 |
Non-Patent Citations (2)
Title |
---|
多源协同的配电网多时段负荷恢复优化决策方法;许寅;王颖;和敬涵;李晨;;电力系统自动化(第02期);全文 * |
考虑韧性增强策略的配电网负荷恢复优化;吴文仙;韩冬;孙伟卿;奚培锋;;电力科学与工程(第04期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN112186744A (zh) | 2021-01-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112186744B (zh) | 一种适用于含分布式电源配电网的供电恢复方法及应用 | |
CN109858759B (zh) | 一种工业园区综合能源平衡调度方法 | |
Lupangu et al. | A review of technical issues on the development of solar photovoltaic systems | |
CN113131468B (zh) | 一种考虑用户停电损失的多目标配电网供电恢复方法 | |
CN109767078B (zh) | 一种基于混合整数规划的多类型电源检修安排方法 | |
CN109888835A (zh) | 一种基于改进粒子群的分布式光伏配电网规划方法 | |
CN108599206A (zh) | 高比例不确定电源场景下的配电网混合储能配置方法 | |
CN101710702A (zh) | 实现电力系统动态节能调度的方法 | |
CN105006844A (zh) | 一种间歇式发电并网条件下的电力系统日前鲁棒调度系统 | |
CN110661301B (zh) | 一种水光蓄多能互补发电系统的容量配置优化方法 | |
CN103577891A (zh) | 一种含分布式电源的多孤岛微网优化合作运行方法 | |
CN111049165A (zh) | 一种用于新能源电力系统储能配置的方法及系统 | |
CN110994606A (zh) | 一种基于复杂适应系统理论的多能源电源容量配置方法 | |
CN113128844A (zh) | 一种基于供电设备容量限制下分布式电源规划方法 | |
Zhou et al. | Optimal sizing of pv system and bess for smart household under stepwise power tariff | |
CN106340901B (zh) | 一种基于惯性矩的最优分位点配电网可控资源配置方法 | |
Le et al. | Design, sizing and operation of a hybrid renewable energy system for farming | |
Yang et al. | Over-voltage mitigation control strategies for distribution system with high PV penetration | |
Laribi et al. | Optimal dimensioning and operation of a grid-supporting energy storage system | |
CN115549138A (zh) | 一种多种互补外送系统中的储能容量优化配置方法及系统 | |
Li et al. | Proportion optimization of wind and solar power capacity for regional power network | |
Julianto et al. | Confronting the Duck Curve Problem Using Dynamic Economic Emission Dispatch with CAES. | |
Laribi et al. | Combined grid-supporting and market-based operation strategy for battery storage systems | |
CN107959293B (zh) | 一种基于最优潮流和分段激励的居民负荷削减调度方法 | |
CN106300443B (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 |