CN111525562B - 一种用于轨道交通供电的电压闪变快速评估方法及系统 - Google Patents

一种用于轨道交通供电的电压闪变快速评估方法及系统 Download PDF

Info

Publication number
CN111525562B
CN111525562B CN202010415827.4A CN202010415827A CN111525562B CN 111525562 B CN111525562 B CN 111525562B CN 202010415827 A CN202010415827 A CN 202010415827A CN 111525562 B CN111525562 B CN 111525562B
Authority
CN
China
Prior art keywords
voltage
carrying
probability
voltage flicker
model
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
Application number
CN202010415827.4A
Other languages
English (en)
Other versions
CN111525562A (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.)
Lanzhou Power Supply Co Of State Grid Gansu Electric Power Co
Original Assignee
Lanzhou Power Supply Co Of State Grid Gansu Electric Power Co
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 Lanzhou Power Supply Co Of State Grid Gansu Electric Power Co filed Critical Lanzhou Power Supply Co Of State Grid Gansu Electric Power Co
Priority to CN202010415827.4A priority Critical patent/CN111525562B/zh
Publication of CN111525562A publication Critical patent/CN111525562A/zh
Application granted granted Critical
Publication of CN111525562B publication Critical patent/CN111525562B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • H02J3/04Circuit arrangements for ac mains or ac distribution networks for connecting networks of the same frequency but supplied from different sources
    • H02J3/06Controlling transfer of power between connected networks; Controlling sharing of load between connected networks
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • H02J3/38Arrangements for parallely feeding a single network by two or more generators, converters or transformers
    • H02J3/46Controlling of the sharing of output between the generators, converters, or transformers
    • H02J3/48Controlling the sharing of the in-phase component
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2203/00Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
    • H02J2203/20Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2310/00The network for supplying or distributing electric power characterised by its spatial reach or by the load
    • H02J2310/40The network being an on-board power network, i.e. within a vehicle

Landscapes

  • Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Supply And Distribution Of Alternating Current (AREA)

Abstract

本发明公开了一种用于轨道交通供电的电压闪变快速评估方法及系统,涉及概率评估方法,其中的一种用于轨道交通供电的电压闪变快速评估方法,包括以下步骤:S1,建立地铁负荷的概率模型;S2,进行概率潮流计算;S3,根据电压的概率分布来进行电压闪变值的计算;S4,将所得的电压闪变值进行数据传输与共享,并根据所得数据进行电压闪变值的预测与预警。本发明用于轨道交通供电的电压闪变快速评估方法能确保轨道交通供电系统的安全性,提供良好的电能质量,并对可能发生闪变事故的线路提前做出预测。

Description

一种用于轨道交通供电的电压闪变快速评估方法及系统
技术领域
本发明涉及概率评估方法,具体涉及一种用于轨道交通供电的电压闪变快速评估方法及系统。
背景技术
目前,我国城市轨道交通的迅速发展,不仅促进了城市建设,而且也带来了一定的社会经济效益。但与此同时,地铁牵引负荷会对公共电网造成的影响也引发了人们的高度重视。因为地铁列车在运行过程中采用直流牵引,所以必然会有大量电子器件的使用,这就造成了地铁负荷具有冲击性强、且系统含有大量谐波等特点。此外,地铁一般建设于较为发达的城市,一旦对公共电网造成事故,可能对城市电力供应和稳定社会生产力造成重大影响。为了提前预防故障的发生,必须通过理论分析制定相应的措施,从而以减少电压波动与闪变对公共电网造成的危害,进而提高对城市供电的电能质量,确保对城市电网供电的可靠性和稳定性,从理论上消除了地铁负荷的接入对电网造成事故的潜在危害。
地铁牵引负荷已成为电力系统负荷中的重要一员,但又因为其具有的强随机性与波动性,本发明提出运用概率潮流的方法去处理电压闪变问题,而目前IEC标准下的闪变检测需要经过一系列的平方检测滤波、带通和加权滤波、平方器低通滤波和统计等流程,在实际数字闪变仪的设计中十分复杂,同时各滤波器的非理想幅频特性引起的误差会对测量结果产生很大的影响。为了在保证电压闪变测量的准确性同时简化检测步骤与模型,需要对闪变测量方法进行改进,用离散化方法处理电压闪变问题。
发明内容
本发明的主要目的在于提供一种用于轨道交通供电的电压闪变快速评估方法及系统,确保轨道交通供电系统的安全性,提供良好的电能质量,并对可能发生闪变事故的线路提前做出预测。
根据本发明的一个方面,提供了一种用于轨道交通供电的电压闪变快速评估方法,包括以下步骤:
S1,建立地铁负荷的概率模型;
S2,进行概率潮流计算;
S3,根据电压的概率分布来进行电压闪变值的计算;
S4,将所得的电压闪变值进行数据传输与共享,并根据所得数据进行电压闪变值的预测与预警。
进一步地,所述步骤S1包括地铁牵引负荷的建模,负荷模型的参数求解,参数求解又包括适应度函数的选取。
更进一步地,所述步骤S2包括采用拟蒙特卡洛采样生成样本,用Nataf变换进行负荷的相关性处理,用半不变量法进行潮流求解,得到供电系统电压的概率分布。
更进一步地,所述步骤S3包括抽取电压序列,对电压序列进行FFT变换,根据电压波动求取瞬时闪变视感度,最后求取电压闪变瞬时值。
更进一步地,所述步骤S4包括电压闪变值数据的处理、电压闪变值模型的建立及预测结果的输出与预警。
根据本发明的又一个方面,提供了一种用于轨道交通供电的电压闪变
快速评估系统,包括:
模型建立模块,用于建立地铁负荷的概率模型;
概率潮流计算模块,用于进行概率潮流计算;
电压闪变值计算模块,用于根据电压的概率分布来进行电压闪变值的
计算;
电压闪变值处理模块,用于将所得的电压闪变值进行数据传输与共享,并根据所得数据进行电压闪变值的预测与预警。
进一步地,所述模型建立模块用于地铁牵引负荷的建模,负荷模型的参数求解,参数求解又包括适应度函数的选取。
更进一步地,所述概率潮流计算模块用于采用拟蒙特卡洛采样生成样本,用Nataf变换进行负荷的相关性处理,用半不变量法进行潮流求解,得到供电系统电压的概率分布。
更进一步地,所述电压闪变值计算模块用于抽取电压序列,对电压序列进行FFT变换,根据电压波动求取瞬时闪变视感度,最后求取电压闪变瞬时值。
更进一步地,所述电压闪变值处理模块用于用于电压闪变值数据的处理、电压闪变值模型的建立及预测结果的输出与预警,电压闪变值模型采用灰色模型、线性模型和随机时间序列模型自适应加权和优选组合得到预测结果。
本发明的优点:
本发明用于轨道交通供电的电压闪变快速评估方法能确保轨道交通供电系统的安全性,提供良好的电能质量,并对可能发生闪变事故的线路提前做出预测;依据评估结果能有效分析地铁负荷的接入带来的电压闪变问题,确定整个轨道供电系统的安全性,掌握牵引负荷接入后未来时刻的电压闪变的变化状态,对改善轨道交通供电系统的结构、评估电压闪变问题以及根据预测模型,提前对电压闪变问题进行处理等方面有极大的进步意义。
除了上面所描述的目的、特征和优点之外,本发明还有其它的目的、特征和优点。下面将参照图,对本发明作进一步详细的说明。
附图说明
构成本申请的一部分的附图用来提供对本发明的进一步理解,本发明的示意性实施例及其说明用于解释本发明,并不构成对本发明的不当限定。
图1是本发明实施例的用于轨道交通供电的电压闪变快速评估方法整
体框图;
图2是本发明实施例的用于轨道交通供电的电压闪变快速评估方法流程图;
图3是本发明实施例的用于轨道交通供电的电压闪变快速评估系统结构框图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。
实施例1
参考图1和图2,如图1和图2所示,一种用于轨道交通供电的电压闪变快速评估方法,包括以下步骤:
S1,建立地铁负荷的概率模型;
S2,进行概率潮流计算;
S3,根据电压的概率分布来进行电压闪变值的计算;
S4,将所得的电压闪变值进行数据传输与共享,并根据所得数据进行电压闪变值的预测与预警。
本发明能确保轨道交通供电系统的安全性,提供良好的电能质量,并对可能发生闪变事故的线路提前做出预测;依据评估结果能有效分析地铁负荷的接入带来的电压闪变问题,确定整个轨道供电系统的安全性,掌握牵引负荷接入后未来时刻的电压闪变的变化状态,对改善轨道交通供电系统的结构、评估电压闪变问题以及根据预测模型,提前对电压闪变问题进行处理等方面有极大的进步意义。
本发明的方法具体包括以下步骤:
地铁牵引负荷的建模:
(1)主要对牵引负荷有功功率分布的影响有:①在某一固定时刻从牵引变电所获取功率的地铁数量;②每台地铁机车的有功功率。每辆地铁从牵引网获取功率的概率分布近似于正态分布,又由于正态分布具有可叠加性,所以可得出牵引负荷概率模型如下:
(2)模型参数的求解。为了避免粒子算法在后期陷入局部最优解,但又得确保在寻优前期能够较大范围地进行搜索,提出采用一种协同粒子群优(CPSO)算法进行参数辨识。
(3)适应度函数的选取。Fitness Function的选取对遗传算法的收敛性有着最显著影响。本文适应度函数选为:
概率潮流计算:
将地铁负荷接入供电系统,进行概率潮流计算。包括采用拟蒙特卡洛采样生成样本,用Nataf变换进行负荷的相关性处理,用半不变量法进行潮流求解,得到电压的概率分布。主要流程如下
(1) 用牛拉法进行交直流潮流的初始计算;
(2) 计算注入功率的累积分布函数,进行QMCS采样;
(3) 进而由可获得注入功率的半不变量;
(4) 计算得到入线口的节点电压的半不变量,由Gram-Charlier级数展开,得到节点电压;
(5) 测试进线口节点电压,计算是否越限,并记录。
根据电压的概率分布来进行电压闪变值的计算:
(1)首先运用拟蒙特卡洛采样法对牵引负荷进行抽样, 然后进行概率潮流计算得到电压值序列
(2)对序列进行FFT变换,得到其离散频谱数列/>
本发明采用FFT分解,该法要求采样点数N为2的某次幂,的频谱分辨率随采样点数的不同而不同,当采样点数取2048次时,即频谱分辨率为0.0244时,/>计算值基本接近于理论值。
(3)根据电压波动的定义可知,/>为/>频率下的电压波动值,其值为正弦电压方均根值曲线幅值的2倍。则,该频率下的/>通过下列公式计算:
对于冲击性负荷引起的电压波动与闪变问题,它们的闪变频率不一定在IEC标准给出的离散频率点上,因此需要对电压波动值的频率点其进行拟合。
(4)一信号均方根值的平方等于该信号的各次谐波均方根值的平方和,因此要计算某一电压序列的瞬时闪变视感度,只需将其经过傅里叶分解后得到各次谐波对应的/>值进行求和:
以上的离散化步骤,得到/>
(5)计算出所述电压信号瞬时闪变值
将所得的电压闪变值进行数据传输与共享,并根据所得数据进行电压闪变值的预测与预警。主要包括:
(1)数据的处理;
(2)模型的建立;
(3)预测结果的输出与预警。
实施例2
参考图3,如图3所示,一种用于轨道交通供电的电压闪变快速评估系
统,包括:
模型建立模块,用于建立地铁负荷的概率模型;
概率潮流计算模块,用于进行概率潮流计算;
电压闪变值计算模块,用于根据电压的概率分布来进行电压闪变值的
计算;
电压闪变值处理模块,用于将所得的电压闪变值进行数据传输与共享,并根据所得数据进行电压闪变值的预测与预警。
所述模型建立模块用于地铁牵引负荷的建模,负荷模型的参数求解,参数求解又包括适应度函数的选取。
所述概率潮流计算模块用于采用拟蒙特卡洛采样生成样本,用Nataf变换进行负荷的相关性处理,用半不变量法进行潮流求解,得到供电系统电压的概率分布。
所述电压闪变值计算模块用于抽取电压序列,对电压序列进行FFT变换,根据电压波动求取瞬时闪变视感度,最后求取电压闪变瞬时值。
所述电压闪变值处理模块用于用于电压闪变值数据的处理、电压闪变值模型的建立及预测结果的输出与预警,电压闪变值模型采用灰色模型、线性模型和随机时间序列模型自适应加权和优选组合得到预测结果。
以上所述仅为本发明的较佳实施例,并不用以限制本发明,凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。

Claims (2)

1.一种用于轨道交通供电的电压闪变快速评估方法,其特征在于,包括以下步骤:
S1,建立地铁负荷的概率模型;
S2,进行概率潮流计算;
S3,根据电压的概率分布来进行电压闪变值的计算;
S4,将所得的电压闪变值进行数据传输与共享,并根据所得数据进行电压闪变值的预测与预警;
所述步骤S1包括地铁牵引负荷的建模,负荷模型的参数求解,参数求解又包括适应度函数的选取;
地铁牵引负荷的建模:
牵引负荷概率模型如下:
模型参数的求解:采用一种协同粒子群优CPSO算法进行参数辨识;
适应度函数的选取:适应度函数选为:
所述步骤S2包括采用拟蒙特卡洛采样生成样本,用Nataf变换进行负荷的相关性处理,用半不变量法进行潮流求解,得到供电系统电压的概率分布;
概率潮流计算:
将地铁负荷接入供电系统,进行概率潮流计算,包括采用拟蒙特卡洛采样生成样本,用Nataf变换进行负荷的相关性处理,用半不变量法进行潮流求解,得到电压的概率分布;主要流程如下:
用牛拉法进行交直流潮流的初始计算;
计算注入功率的累积分布函数,进行QMCS采样;
进而由可获得注入功率的半不变量;
计算得到入线口的节点电压的半不变量,由Gram-Charlier级数展开,得到节点电压;
测试进线口节点电压,计算是否越限,并记录;
所述步骤S3包括抽取电压序列,对电压序列进行FFT变换,根据电压波动求取瞬时闪变视感度,最后求取电压闪变瞬时值;
根据电压的概率分布来进行电压闪变值的计算:
(1)首先运用拟蒙特卡洛采样法对牵引负荷进行抽样,然后进行概率潮流计算得到电压值序列U(N):
(2)对序列U(N)进行FFT变换,得到其离散频谱数列Uf(k)
Uf(k)=FFT[U(N)]
采用FFT分解,该法要求采样点数N为2的某次幂,U(N)的频谱分辨率随采样点数的不同而不同,当采样点数取2048次时,即频谱分辨率为0.0244时,S(t)计算值基本接近于理论值;
(3)根据电压波动的定义ΔU=Umax-Umin可知,d(i)为i频率下的电压波动值,其值为正弦电压方均根值曲线幅值的2倍;则,该频率下的Si(t)通过下列公式计算:
Si(t)=[d(i)/di]2
(4)计算某一电压序列U(N)的瞬时闪变视感度,将其经过傅里叶分解后得到各次谐波对应的Si(t)值进行求和:
以上是S(t)的离散化步骤,得到S(t);
(5)计算出电压信号瞬时闪变值
所述步骤S4包括电压闪变值数据的处理、电压闪变值模型的建立及预测结果的输出与预警,电压闪变值模型采用灰色模型、线性模型和随机时间序列模型自适应加权和优选组合得到预测结果。
2.一种如权利要求1所述的用于轨道交通供电的电压闪变快速评估方法的系统,其特征在于,包括:
模型建立模块,用于建立地铁负荷的概率模型;
概率潮流计算模块,用于进行概率潮流计算;
电压闪变值计算模块,用于根据电压的概率分布来进行电压闪变值的计算;
电压闪变值处理模块,用于将所得的电压闪变值进行数据传输与共享,并根据所得数据进行电压闪变值的预测与预警;
所述模型建立模块用于地铁牵引负荷的建模,负荷模型的参数求解,参数求解又包括适应度函数的选取;
所述概率潮流计算模块用于采用拟蒙特卡洛采样生成样本,用Nataf变换进行负荷的相关性处理,用半不变量法进行潮流求解,得到供电系统电压的概率分布;
所述电压闪变值计算模块用于抽取电压序列,对电压序列进行FFT变换,根据电压波动求取瞬时闪变视感度,最后求取电压闪变瞬时值;
所述电压闪变值处理模块用于用于电压闪变值数据的处理、电压闪变值模型的建立及预测结果的输出与预警。
CN202010415827.4A 2020-05-16 2020-05-16 一种用于轨道交通供电的电压闪变快速评估方法及系统 Active CN111525562B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010415827.4A CN111525562B (zh) 2020-05-16 2020-05-16 一种用于轨道交通供电的电压闪变快速评估方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010415827.4A CN111525562B (zh) 2020-05-16 2020-05-16 一种用于轨道交通供电的电压闪变快速评估方法及系统

Publications (2)

Publication Number Publication Date
CN111525562A CN111525562A (zh) 2020-08-11
CN111525562B true CN111525562B (zh) 2023-07-21

Family

ID=71912457

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010415827.4A Active CN111525562B (zh) 2020-05-16 2020-05-16 一种用于轨道交通供电的电压闪变快速评估方法及系统

Country Status (1)

Country Link
CN (1) CN111525562B (zh)

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103235981A (zh) * 2013-04-10 2013-08-07 东南大学 一种风电电能质量趋势预测方法
CN106846171A (zh) * 2016-12-27 2017-06-13 广州耐奇电气科技有限公司 一种城市轨道交通交直流系统电能管理系统
CN109038591A (zh) * 2018-08-16 2018-12-18 广东电网有限责任公司 一种半不变量法概率潮流计算方法
CN109066688A (zh) * 2018-09-06 2018-12-21 国网安徽省电力有限公司芜湖供电公司 基于可再生能源不确定性下的概率潮流数据获取方法
CN110940877A (zh) * 2019-11-29 2020-03-31 长园深瑞监测技术有限公司 一种基于加窗fft插值算法的电压闪变检测方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103235981A (zh) * 2013-04-10 2013-08-07 东南大学 一种风电电能质量趋势预测方法
CN106846171A (zh) * 2016-12-27 2017-06-13 广州耐奇电气科技有限公司 一种城市轨道交通交直流系统电能管理系统
CN109038591A (zh) * 2018-08-16 2018-12-18 广东电网有限责任公司 一种半不变量法概率潮流计算方法
CN109066688A (zh) * 2018-09-06 2018-12-21 国网安徽省电力有限公司芜湖供电公司 基于可再生能源不确定性下的概率潮流数据获取方法
CN110940877A (zh) * 2019-11-29 2020-03-31 长园深瑞监测技术有限公司 一种基于加窗fft插值算法的电压闪变检测方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
基于拟蒙特卡洛的半不变量法概率潮流计算;栗然等;《电力建设》;20171101(第11期);第145-147页 *
计及牵引负荷相关性的随机潮流计算;高锋阳;《铁道学报》;20181031;第40卷(第10期);第24-27页 *

Also Published As

Publication number Publication date
CN111525562A (zh) 2020-08-11

Similar Documents

Publication Publication Date Title
Rajesh et al. Optimally detecting and classifying the transmission line fault in power system using hybrid technique
CN105098763A (zh) 一种基于wams和ems的暂态电压稳定在线评估方法
Weckx et al. Parameter identification of unknown radial grids for theft detection
CN111401476B (zh) 基于边界区域重要性采样与核向量机的暂态安全评估方法
CN110703151A (zh) 基于振动盲源分离和贝叶斯模型的变压器故障诊断方法
Kaplan et al. A novel method based on Weibull distribution for short-term wind speed prediction
Samet et al. A wide nonlinear analysis of reactive power time series related to electric arc furnaces
CN108270239A (zh) 一种含分布式电源的配电网电能质量扰动源方向判定方法
CN115684776A (zh) 高压输电线路线损支路的确定方法及其装置、处理器
CN111525562B (zh) 一种用于轨道交通供电的电压闪变快速评估方法及系统
Jaramillo et al. Identification of distributed energy resources in low voltage distribution networks
Mahela et al. A protection scheme for distribution utility grid with wind energy penetration
Betancourt et al. A spatio-temporal processing Padé approach for visualizing harmonic distortion propagation on electrical networks
CN115343579B (zh) 一种电网故障分析方法、装置及电子设备
Lima et al. Hardware and software architecture for power quality analysis
Banerjee et al. Data-driven static load model parameter estimation with confidence factor
CN113177717B (zh) 一种基于影响增量灵敏度的输电系统韧性快速评估方法
CN105277852A (zh) 配电网线路工况分类识别方法
Yue et al. A power system harmonic problem based on the BP neural network learning algorithm
CN113640733A (zh) 一种台区负荷还原模拟仿真系统及其使用方法
Anshuman et al. A novel hybrid algorithm for event detection, localisation and classification
Li et al. Wide area synchronized measurements and inter-area oscillation study
Turunen et al. Recent developments in modal estimation of power system electromechanical oscillations
Samet et al. Quantizing branches' responsibility for voltage sags/swells and flicker events
Chen et al. A Spatial-temporal Frequency Nadir Prediction Method Based on ConvLSTM with Attention

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