CN109861280A - 基于虚拟同步发电机的微网变流器频率控制方法及系统 - Google Patents

基于虚拟同步发电机的微网变流器频率控制方法及系统 Download PDF

Info

Publication number
CN109861280A
CN109861280A CN201910086221.8A CN201910086221A CN109861280A CN 109861280 A CN109861280 A CN 109861280A CN 201910086221 A CN201910086221 A CN 201910086221A CN 109861280 A CN109861280 A CN 109861280A
Authority
CN
China
Prior art keywords
virtual synchronous
machine model
synchronous machine
current transformer
frequency
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.)
Granted
Application number
CN201910086221.8A
Other languages
English (en)
Other versions
CN109861280B (zh
Inventor
肖飞
王林
黄辉
魏亚龙
李长喜
唐启迪
祁招
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
State Grid Corp of China SGCC
Xuji Group Co Ltd
XJ Electric Co Ltd
Xian XJ Power Electronics Technology Co Ltd
Original Assignee
State Grid Corp of China SGCC
Xuji Group Co Ltd
XJ Electric Co Ltd
Xian XJ Power Electronics Technology Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by State Grid Corp of China SGCC, Xuji Group Co Ltd, XJ Electric Co Ltd, Xian XJ Power Electronics Technology Co Ltd filed Critical State Grid Corp of China SGCC
Priority to CN201910086221.8A priority Critical patent/CN109861280B/zh
Publication of CN109861280A publication Critical patent/CN109861280A/zh
Application granted granted Critical
Publication of CN109861280B publication Critical patent/CN109861280B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02BCLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO BUILDINGS, e.g. HOUSING, HOUSE APPLIANCES OR RELATED END-USER APPLICATIONS
    • Y02B70/00Technologies for an efficient end-user side electric power management and consumption
    • Y02B70/30Systems 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/3225Demand response systems, e.g. load shedding, peak shaving
    • YGENERAL 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
    • Y04INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
    • Y04SSYSTEMS 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/00Management or operation of end-user stationary applications or the last stages of power distribution; Controlling, monitoring or operating thereof
    • Y04S20/20End-user application control systems
    • Y04S20/222Demand response systems, e.g. load shedding, peak shaving

Landscapes

  • Control Of Eletrric Generators (AREA)

Abstract

本发明提供一种基于虚拟同步发电机的微网变流器频率控制方法及系统,属于微网变流器频率控制技术领域。该控制方法包括以下步骤:1)微网变流器采用虚拟同步机模型进行控制,通过对虚拟同步机模型的机械转矩进行控制,实现微网变流器频率的一次调节;2)在频率一次调节下,当微网变流器输出功率与负荷达到平衡时,以虚拟同步机模型的电磁功率作为动态负荷补偿量对虚拟同步机模型的机械功率进行补偿。本发明的有益效果是:当微网变流器输出功率与负荷达到平衡时,通过利用虚拟同步机模型的电磁功率作为动态负荷补偿量对虚拟同步机模型的机械功率进行补偿,实现了微网变流器频率的二次调节,使得微网变流器频率稳定在额定频率,实现简单。

Description

基于虚拟同步发电机的微网变流器频率控制方法及系统
技术领域
本发明涉及一种基于虚拟同步发电机的微网变流器频率控制方法及系统,属于微网变流器频率控制技术领域。
背景技术
对于交流微电网系统,频率稳定是系统可靠运行的必要条件。微网变流器频率控制最常见的是采用有功-频率下垂控制:微网变流器根据下垂曲线和检测的频率自适应调整输出功率和运行频率点,下垂控制为一次调节,频率偏移额定工作点。为实现频率的无差控制,通常会增加频率的二次调节,如采用PI等无差控制器进行调节。无差控制器的引入增大了控制参数设计难度,同时引入新的振荡模态,增加不稳定风险。
虚拟同步发电机(VSG)技术是近年来的研究热点,其将同步发电机本体及其控制器数学模型引入变流器的控制算法,通过调整内电势幅值和功角调节有功、无功输出。VSG技术用于微网变流器,可使变流器调频、调压具备惯量/阻尼特性,有效抑制并网点频率振荡,增强并网点电压强度。例如公布号为CN107887926A的发明专利申请文件,公开了一种园区型负荷虚拟同步协调控制系统,该控制系统根据频率变化情况,通过功频器控制机械功率改变实现一次调频,通过改变频率曲线使系统重新回到额定频率实现二次调频。该文件虽然公开了二次调频的内容,但是并没有给出具体手段。和常规变流器一样,基于VSG技术的微网变流器频率二次调节同样存在参数设计复杂和不稳定风险问题。
发明内容
本发明的目的是提供一种基于虚拟同步发电机的微网变流器频率控制方法及系统,用以解决目前采用VSG技术进行微电网变流器频率控制,存在的参数设计复杂、难以实现的问题。
为实现上述目的,本发明提供了一种基于虚拟同步发电机的微网变流器频率控制方法,该控制方法包括以下步骤:
1)微网变流器采用虚拟同步机模型进行控制,通过对虚拟同步机模型的机械转矩进行控制,实现微网变流器频率的一次调节;
2)在频率一次调节下,当微网变流器输出功率与负荷达到平衡时,以虚拟同步机模型的电磁功率作为动态负荷补偿量对虚拟同步机模型的机械功率进行补偿。
本发明还提供了一种基于虚拟同步发电机的微网变流器频率控制系统,包括处理器和存储器,所述处理器用于运行存储在所述存储器中的程序指令,以实现如下方法:
1)微网变流器采用虚拟同步机模型进行控制,通过对虚拟同步机模型的机械转矩进行控制,实现微网变流器频率的一次调节;
2)在频率一次调节下,当微网变流器输出功率与负荷达到平衡时,以虚拟同步机模型的电磁功率作为动态负荷补偿量对虚拟同步机模型的机械功率进行补偿。
本发明的基于虚拟同步发电机的微网变流器频率控制方法及系统的有益效果是:当微网变流器输出功率与负荷达到平衡时,通过利用虚拟同步机模型的电磁功率作为动态负荷补偿量对虚拟同步机模型的机械功率进行补偿,实现了微网变流器频率的二次调节,使得微网变流器频率稳定在额定频率,实现简单。
为了使得到的动态负荷补偿量更加精准,进一步提高频率调节的稳定性,作为对上述基于虚拟同步发电机的微网变流器频率控制方法及系统的一种改进,所述动态负荷补偿量通过将虚拟同步机模型的电磁功率进行低通滤波和求平均得到。
为了实现低通滤波,作为对上述基于虚拟同步发电机的微网变流器频率控制方法及系统的另一种改进,所述低通滤波采用的低通滤波器的传递函数G(s)为:
其中,G0为低通滤波器增益,ωn为低通滤波器截止频率,ξ为阻尼系数。
为了更精确地对虚拟同步机模型的机械转矩进行控制,作为对上述基于虚拟同步发电机的微网变流器频率控制方法及系统的又一种改进,所述步骤1)采用功频控制器对虚拟同步机模型的机械转矩进行控制,功频控制器的模型描述为:
Tm=[Pref-kf(f-f0)PN]/ω
其中,Tm为虚拟同步机模型的机械转矩,f、f0分别为虚拟同步机模型的机端电压频率、机端电压额定频率,PN为虚拟同步机模型的额定功率,Pref为有功功率指令值,kf为一次调频系数,ω为虚拟同步机模型的角速度。
为了更精确地对虚拟同步机模型的机械功率进行补偿,作为对上述基于虚拟同步发电机的微网变流器频率控制方法及系统的再一种改进,利用所述动态负荷补偿量对虚拟同步机模型的机械功率进行补偿,具体为:
Pm′=kf(fn-fr)+Pe_avg
其中,Pm′为补偿后的虚拟同步机模型的机械功率,kf为一次调频系数,fn为虚拟同步机模型的额定频率,fr为虚拟同步机模型的转子频率,Pe_avg为动态负荷补偿量。
为了更精确地对微网变流器频率进行调节,作为对上述基于虚拟同步发电机的微网变流器频率控制方法及系统的进一步改进,所述虚拟同步机模型为:
其中,ω、ω0分别为虚拟同步机模型的角速度、额定角速度,H为虚拟同步机模型的惯性时间常数,Tm、Te分别为虚拟同步机模型的机械转矩、电磁转矩,Kd为阻尼系数,R、L分别为虚拟同步机模型的定子电阻、定子电感,eabc、uabc、iabc分别为虚拟同步机模型的三相内电势、三相机端电压、三相定子电流,Pe为虚拟同步机模型的电磁功率,ia、ib、ic为单相定子电流,ea、eb、ec为单相内电势。
附图说明
图1是本发明的虚拟同步机模型图;
图2是本发明的机械功率补偿框图;
图3是本发明的频率二次调节原理图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及具体实施例对本发明进行进一步详细说明。
基于虚拟同步发电机的微网变流器频率控制方法实施例:
本发明的基于虚拟同步发电机的微网变流器频率控制方法,包括以下步骤:
1)微网变流器采用虚拟同步机模型进行控制,通过对虚拟同步机模型的机械转矩进行控制,实现微网变流器频率的一次调节。
其中,本发明采用功频控制器对虚拟同步机模型(也称VSG)的机械转矩进行控制,从而实现微网变流器频率的一次调节,其中VSG的机械转矩由功频控制器计算得到。
如图1所示,虚拟同步机模型为:
其中,ω、ω0分别为VSG的角速度、额定角速度,H为VSG的惯性时间常数,Tm、Te分别为VSG的机械转矩、电磁转矩,Pe为VSG的电磁功率,Kd为阻尼系数,R、L分别为VSG的定子电阻、定子电感,eabc、uabc、iabc分别为VSG的三相内电势、三相机端电压、三相定子电流,ia、ib、ic为单相定子电流,ea、eb、ec为单相内电势,θ为VSG的转子相角。
功频控制器的模型为:
Tm=[Pref-kf(f-f0)PN]/ω
其中,Tm为VSG的机械转矩,f、f0分别为VSG的机端电压频率、机端电压额定频率,PN为VSG的额定功率,Pref为有功功率指令值,kf为一次调频系数,ω为VSG的角速度。
实现微网变流器频率的一次调节的具体过程如下:
当微网变流器运行在稳态工作点时,虚拟同步机模型的机械功率与电磁功率达到平衡;
当负荷突增时,虚拟同步机模型输出电磁功率随之增大,其机械功率与电磁功率失衡,微网变流器频率降低。一次调节控制器(即功频控制器)根据频率偏差增大虚拟同步机模型的机械功率,使其机械功率与电磁功率重新达到稳态,但由于一次调节为有差调节,微网变流器新的稳态工作点频率向下偏移;
当负荷突减时,虚拟同步机模型输出电磁功率随之减小,其机械功率与电磁功率失衡,微网变流器频率升高。一次调节控制器根据频率偏差减小虚拟同步机模型的机械功率,使其机械功率与电磁功率重新达到稳态,但由于一次调节为有差调节,微网变流器新的稳态工作点频率向上偏移。
2)在频率一次调节下,当微网变流器输出功率与负荷达到平衡时,以虚拟同步机模型的电磁功率作为动态负荷补偿量对虚拟同步机模型的机械功率进行补偿。
其中,本发明通过将虚拟同步机模型的电磁功率进行低通滤波和求平均,得到动态负荷补偿量,并采用动态负荷补偿量对虚拟同步机模型的机械功率进行补偿,实现微网变流器频率的二次调节,从而使微网变流器频率稳定在额定频率。
动态负荷补偿量的计算过程如下:
由于稳态下,VSG的机械功率Pm、电磁功率Pe、系统负荷PL三者近似相等,因此,可将电磁功率Pe作为负荷动态估算值。同时,为防止补偿过程中的振荡失稳,需对电磁功率Pe进行低通滤波和求平均处理,最终得到动态负荷补偿量Pe_avg
低通滤波器的传递函数如下:
其中,G0为滤波器增益,ωn为滤波器截止频率,ξ为阻尼系数。工程上一般设置G0=1,ξ按“二阶最优系统”整定为0.707,ωn根据仿真或试验结果进行调整,优选ωn0/5。
如图2所示,利用动态负荷补偿量对虚拟同步机模型的机械功率进行补偿,具体为:
Pm′=kf(fn-fr)+Pe_avg
其中,Pm′为补偿后的虚拟同步机模型的机械功率,kf为一次调频系数,fn为虚拟同步机模型的额定频率,fr为虚拟同步机模型的转子频率,Pe_avg为动态负荷补偿量。
如图3所示,在功频控制器一次调节控制下,当微网变流器输出功率P与负荷PL达到平衡时,微网变流器运行于工作点A,微网变流器输出频率相对其额定频率fN的偏移量为-△f;加入频率二次调节后,一次调节频率-功率下垂曲线右移,稳态时,微网变流器运行于工作点A′,微网变流器输出额定频率fN,负荷仍为PL
基于虚拟同步发电机的微网变流器频率控制系统实施例:
本发明的基于虚拟同步发电机的微网变流器频率控制系统,包括处理器和存储器,处理器用于运行存储在存储器中的程序指令,以实现本发明的基于虚拟同步发电机的微网变流器频率控制方法,具体方法参见基于虚拟同步发电机的微网变流器频率控制方法实施例。

Claims (7)

1.一种基于虚拟同步发电机的微网变流器频率控制方法,其特征在于,该控制方法包括以下步骤:
1)微网变流器采用虚拟同步机模型进行控制,通过对虚拟同步机模型的机械转矩进行控制,实现微网变流器频率的一次调节;
2)在频率一次调节下,当微网变流器输出功率与负荷达到平衡时,以虚拟同步机模型的电磁功率作为动态负荷补偿量对虚拟同步机模型的机械功率进行补偿。
2.根据权利要求1所述的基于虚拟同步发电机的微网变流器频率控制方法,其特征在于,所述动态负荷补偿量通过将虚拟同步机模型的电磁功率进行低通滤波和求平均得到。
3.根据权利要求2所述的基于虚拟同步发电机的微网变流器频率控制方法,其特征在于,所述低通滤波采用的低通滤波器的传递函数G(s)为:
其中,G0为低通滤波器增益,ωn为低通滤波器截止频率,ξ为阻尼系数。
4.根据权利要求1所述的基于虚拟同步发电机的微网变流器频率控制方法,其特征在于,所述步骤1)采用功频控制器对虚拟同步机模型的机械转矩进行控制,功频控制器的模型描述为:
Tm=[Pref-kf(f-f0)PN]/ω
其中,Tm为虚拟同步机模型的机械转矩,f、f0分别为虚拟同步机模型的机端电压频率、机端电压额定频率,PN为虚拟同步机模型的额定功率,Pref为有功功率指令值,kf为一次调频系数,ω为虚拟同步机模型的角速度。
5.根据权利要求1所述的基于虚拟同步发电机的微网变流器频率控制方法,其特征在于,利用所述动态负荷补偿量对虚拟同步机模型的机械功率进行补偿,具体为:
P′m=kf(fn-fr)+Pe_avg
其中,P′m为补偿后的虚拟同步机模型的机械功率,kf为一次调频系数,fn为虚拟同步机模型的额定频率,fr为虚拟同步机模型的转子频率,Pe_avg为动态负荷补偿量。
6.根据权利要求1所述的基于虚拟同步发电机的微网变流器频率控制方法,其特征在于,所述虚拟同步机模型为:
其中,ω、ω0分别为虚拟同步机模型的角速度、额定角速度,H为虚拟同步机模型的惯性时间常数,Tm、Te分别为虚拟同步机模型的机械转矩、电磁转矩,Kd为阻尼系数,R、L分别为虚拟同步机模型的定子电阻、定子电感,eabc、uabc、iabc分别为虚拟同步机模型的三相内电势、三相机端电压、三相定子电流,Pe为虚拟同步机模型的电磁功率,ia、ib、ic为单相定子电流,ea、eb、ec为单相内电势。
7.一种基于虚拟同步发电机的微网变流器频率控制系统,其特征在于,包括处理器和存储器,所述处理器用于运行存储在所述存储器中的程序指令,以实现权利要求1-6任一项所述方法。
CN201910086221.8A 2019-01-29 2019-01-29 基于虚拟同步发电机的微网变流器频率控制方法及系统 Active CN109861280B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910086221.8A CN109861280B (zh) 2019-01-29 2019-01-29 基于虚拟同步发电机的微网变流器频率控制方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910086221.8A CN109861280B (zh) 2019-01-29 2019-01-29 基于虚拟同步发电机的微网变流器频率控制方法及系统

Publications (2)

Publication Number Publication Date
CN109861280A true CN109861280A (zh) 2019-06-07
CN109861280B CN109861280B (zh) 2022-07-22

Family

ID=66896705

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910086221.8A Active CN109861280B (zh) 2019-01-29 2019-01-29 基于虚拟同步发电机的微网变流器频率控制方法及系统

Country Status (1)

Country Link
CN (1) CN109861280B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113746114A (zh) * 2021-09-02 2021-12-03 中国电建集团江西省电力设计院有限公司 优化孤岛电网调频的控制方法及系统

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104638679A (zh) * 2015-02-06 2015-05-20 芜湖大学科技园发展有限公司 一种采用自适应调节的虚拟同步发电机频率控制方法
CN106026196A (zh) * 2016-07-28 2016-10-12 西安交通大学 一种基于逆变器的多虚拟同步发电机并联组网控制方法
CN106159999A (zh) * 2016-07-29 2016-11-23 合肥工业大学 基于虚拟同步发电机的柴储混合独立微网动态频率稳定策略
CN106684921A (zh) * 2017-03-20 2017-05-17 重庆大学 基于虚拟同步发电机的逆变器二次调频控制电路
CN108429431A (zh) * 2018-03-12 2018-08-21 许继集团有限公司 一种基于虚拟同步发电机的变流器及其控制方法
CN108429289A (zh) * 2018-03-19 2018-08-21 许继电气股份有限公司 一种基于虚拟同步发电机的控制方法及系统
US20180269819A1 (en) * 2017-03-14 2018-09-20 Abb Schweiz Ag Method And Control System For Controlling A Power Converter

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104638679A (zh) * 2015-02-06 2015-05-20 芜湖大学科技园发展有限公司 一种采用自适应调节的虚拟同步发电机频率控制方法
CN106026196A (zh) * 2016-07-28 2016-10-12 西安交通大学 一种基于逆变器的多虚拟同步发电机并联组网控制方法
CN106159999A (zh) * 2016-07-29 2016-11-23 合肥工业大学 基于虚拟同步发电机的柴储混合独立微网动态频率稳定策略
US20180269819A1 (en) * 2017-03-14 2018-09-20 Abb Schweiz Ag Method And Control System For Controlling A Power Converter
CN106684921A (zh) * 2017-03-20 2017-05-17 重庆大学 基于虚拟同步发电机的逆变器二次调频控制电路
CN108429431A (zh) * 2018-03-12 2018-08-21 许继集团有限公司 一种基于虚拟同步发电机的变流器及其控制方法
CN108429289A (zh) * 2018-03-19 2018-08-21 许继电气股份有限公司 一种基于虚拟同步发电机的控制方法及系统

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
张辉 等: "离网运行的多微源微网调频调压控制策略", 《高电压技术》 *
李斌 等: "基于改进虚拟同步发电机算法的微网逆变器二次调频方案", 《电网技术》 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113746114A (zh) * 2021-09-02 2021-12-03 中国电建集团江西省电力设计院有限公司 优化孤岛电网调频的控制方法及系统

Also Published As

Publication number Publication date
CN109861280B (zh) 2022-07-22

Similar Documents

Publication Publication Date Title
CN108429289B (zh) 一种基于虚拟同步发电机的控制方法及系统
KR20190127913A (ko) 전력 공급 네트워크에 전력을 공급하기 위한 방법
CN109659956A (zh) 多机并联运行的虚拟同步发电机惯性阻尼混合控制方法
CN110739721A (zh) 一种电压源型风电机组控制方法及系统
CN107370181B (zh) 并网控制方法及系统
CN112467784A (zh) 一种混合微网换流器自适应虚拟同步机控制方法
CN109888776B (zh) 针对直驱型风电场次同步谐振频率的预测方法及终端设备
CN110365051A (zh) 一种自适应指令滤波反演的虚拟同步电机控制方法
Khazaei et al. Impedance-model-based MIMO analysis of power synchronization control
CN111697595B (zh) 一种虚拟同步风机的频率动态优化控制方法及系统
CN109617082A (zh) 一种抑制微电网电压频率越限的暂态稳定控制方法
CN115882762A (zh) 构网型风电并网系统的频率优化控制方法
CN103762618B (zh) 一种具有致稳能力的发电系统及控制方法
CN113659611B (zh) 一种并网模式下虚拟同步发电机虚拟惯量的控制方法
CN116632866B (zh) 一种液流超容锂电池混合储能自适应惯量vsg控制方法
CN116667388B (zh) 一种液流超容锂电池混合储能抑制电力系统低频振荡方法
WO2024066537A1 (zh) 一种供电系统及构网控制方法
CN109861280A (zh) 基于虚拟同步发电机的微网变流器频率控制方法及系统
CN116683491A (zh) 一种新能源微电网的惯量控制方法
CN111146786A (zh) 一种附加阻尼调节的虚拟同步发电机功率控制方法和系统
CN113964858B (zh) 一种基于对偶同步原理的三相逆变器并网控制系统
CN110854883B (zh) 一种用于增强系统阻尼的储能控制方法及系统
CN109004680B (zh) 基于储能逆变器的风电场功率控制方法与系统
Ling et al. Adaptive control of grid connected inverter based on virtual synchronous machine
Xie et al. Improved active power dynamic response of voltage‐controlled doubly fed induction generator

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