CN103001214A - 一种基于中性点偏移的配电网三相潮流计算方法 - Google Patents
一种基于中性点偏移的配电网三相潮流计算方法 Download PDFInfo
- Publication number
- CN103001214A CN103001214A CN201210392893XA CN201210392893A CN103001214A CN 103001214 A CN103001214 A CN 103001214A CN 201210392893X A CN201210392893X A CN 201210392893XA CN 201210392893 A CN201210392893 A CN 201210392893A CN 103001214 A CN103001214 A CN 103001214A
- Authority
- CN
- China
- Prior art keywords
- centerdot
- phase
- node
- neutral point
- distribution network
- 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
Links
Images
Landscapes
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明涉及一种基于中性点偏移的配电网三相潮流计算方法,包括下述步骤:获取中压配电网网络拓扑,对计算馈线上各个节点和各条支路进行编号;引入各节点的中性点电压偏移量;采用所述中性点偏移量对所述节点的计算相电压值进行修正;采用所述修正后的计算相电压值计算节点处的负荷电流值;采用前推回代法计算各支路的三相支路电流及各节点的三相电压值;采用牛顿-拉夫逊法进行迭代,修正中性点电压偏移量;采用混合迭代法对三相电压进行求解;判断步骤G中的三相电压值是否满足迭代终止条件;计算三相潮流分布并输出结果。本发明引入中性点偏移量来处理三相电流之间的约束关系,采用前推回代与牛顿-拉夫逊结合的混合迭代法,利用二者的优点有效解决配电网三相潮流计算问题。
Description
技术领域
本发明涉及电力系统运行分析与控制技术领域,具体涉及一种基于中性点偏移的配电网三相潮流计算方法。
背景技术
潮流计算是电力系统分析中最基本、最重要的计算,是电力系统规划、运行分析、控制与优化的基础。从应用角度来看、电力系统的潮流计算可以分为两种,一种是基于电力系统的系统参数对称、负荷对称以及系统处于正常运行状态的单相潮流计算,另一种是基于系统网络参数三相不对称的三相潮流计算。长期以来,专家学者对输电网的潮流计算进行了大量深入的研究,输电系统一般是三相对称系统,通常采用单相潮流计算。随着我国国民经济的蓬勃发展及配电网在电力系统中的作用日益受到重视,配电网的潮流计算研究也开始受到关注。
配电网具有诸多不同于输电网络的特征,如开环运行(闭环设计、辐射状运行)、三相不平衡情况突出、线路电阻和电抗比值大、节点数目大等,使得传统的潮流计算方法不适用于配电网。为了提高配电网的自动化管理水平、保障配电网的安全可靠运行,必需对配电网进行及时、准确的潮流计算,因而有必要研究适于配电网的潮流计算方法。根据配电网的特征,不少专家学者提出了多种适于配电网潮流计算方法,如改进牛拉法,回路阻抗法,前推回代法等。前推回代法的收敛和迭代速度都比较快,计算量相对较小,适合辐射状配电网潮流计算,在单相潮流计算中得到了很好的验证;牛顿-拉夫逊迭代法具有二次收敛性,在病态情况下不易发散,适于大型潮流计算。
针对三相潮流问题,目前常用的方法是在考虑三相互阻抗后将三相线路等效成单相线路,再采用传统的牛顿法、前推回代法、回路阻抗法、近似解耦法等进行计算,未考虑三相电流之间的相互制约关系;或将中压侧各相输出的负荷大小等效成低压侧相应相的负荷量,再利用前推回代法进行计算,只适用于特定接线方式下的三相潮流计算。
发明内容
针对现有技术的不足,本发明提供一种基于中性点偏移的配电网三相潮流计算方法,本发明方法引入中性点偏移量来处理三相电流之间的约束关系,采用前推回代与牛顿-拉夫逊结合的混合迭代法,利用二者的优点有效解决配电网三相潮流计算问题。
本发明的目的是采用下述技术方案实现的:
一种基于中性点偏移的配电网三相潮流计算方法,其改进之处在于,所述方法包括下述步骤:
A、获取中压配电网网络拓扑,选取中压馈线作为计算单元,对馈线上各个节点和各条支路(馈线采用节点-支路模型,馈线上每一个需要计算电压电流值的点为一个节点,两节点之间的线路为支路)进行编号;
B、引入中性点电压偏移量;
C、采用所述中性点偏移量对所述节点的计算相电压值进行修正;
D、采用所述修正后的计算相电压值计算节点处的负荷电流值;
E、采用前推回代法计算各支路的三相电流值及各节点的三相电压值;
F、采用牛顿-拉夫逊法进行迭代,修正中性点电压偏移量;
G、采用混合迭代法对三相电压进行求解;
H、判断步骤G中的三相电压值是否满足迭代终止条件;
I、计算三相潮流分布并输出结果。
其中,所述步骤A中,所述中压配电网采用辐射状配电网或少环状配电网,在三相潮流计算中,选用中压馈线作为基本计算单元,馈线首端的中压母线作为三相对称的电压源;并对馈线上各支路和各节点进行编号。
其中,所述步骤C中,利用所述中性点电压偏移量用对相电压值进行修正,计算三相电压值用如下①式表示:
其中,所述步骤D中,节点i的A、B、C三相负荷分别为Sai、Sbi、Sci时,A、B、C三相上的负荷电流用下述②式表示:
得到三相负荷电流相互制约关系满足下式:
其中:M为与节点i直接相连的所有下层支路集;对于辐射状馈线,上下层是按电流流向划分的,对于节点i和节点j,若实际电源是从i流向j,则i为j的上层节点;若实际电源从j流向i,则i为j的下层节点;分别表示节点j在A、B、C三相支路电流。
其中,没有下层节点的节点为末端节点,与末端支路直接相连的支路为末端支路,对于末端支路有:
所述步骤E中,采用回代法计算节点的三相电压值,用下述⑦式表示:
对于首端节点的三相电压值,用下述⑧式表示:
其中:分别为首端节点在A、B、C三相电压值。
其中,所述步骤F中,采用牛顿-拉夫逊法,得中性点电压偏移量的迭代公式,用下述⑨式表示:
其中,所述步骤G中,利用所述中性点偏移量对三相相电压进行修正后,采用一次前推回代法得三相相电压值,利用牛顿-拉夫逊迭代法修正中性点偏移量,经过混合迭代即得三相相电压。
其中,所述步骤H中,所述迭代终止条件用下述式⑩表示:
其中,所述方法的负载接入方式为星形或三角形。
与现有技术比,本发明达到的有益效果是:
1、本发明提供的基于中性点偏移的配电网三相潮流计算方法,直接利用相电压电流进行计算,原理简单直观;
2、本发明提供的基于中性点偏移的配电网三相潮流计算方法,通过引入中性点电压偏移量对各相电压进行修正,合理考虑了三相电压电流之间的相互制约关系,使三相潮流计算模型更准确,同时也避免了对三相进行解耦;
3、本发明提供的基于中性点偏移的配电网三相潮流计算方法,采用混合迭代法进行求解,综合了前推回代与牛顿拉夫逊法的优化,计算三相潮流时迭代收敛性好,三相平衡时只需迭代几次,负荷重载或不平衡时也只需迭代几次到十几次,计算速度快;
4、本发明提供的基于中性点偏移的配电网三相潮流计算方法,既适应于三角形负载情况下,也适应于星形接线负载情况。
5、本发明提供的基于中性点偏移的配电网三相潮流计算方法,既适应于辐射状配电网,也适应于弱环配电网的潮流计算。
附图说明
图1是本发明提供的基于中性点偏移的各相相电压示意图;
图2是本发明提供的基于中性点偏移的三相潮流计算方法流程图。
具体实施方式
下面结合附图对本发明的具体实施方式作进一步的详细说明。
本发明提供的基于中性点偏移的配电网三相潮流计算方法,在已知A、B、C三相各相负荷的条件下,在逐次叠代潮流计算过程中,利用中性点电压偏移量,每次叠代中采用计算节点处的中性点电压偏移量修正该节点的各相计算相电压值,并采用修正后的计算相电压值计算该节点处的负荷电流,得到负荷电流后采用前推回代与牛拉法相结合的混合迭代法求解各相相电压,经数次迭代收敛最终得到潮流分布。本发明的原理及具体技术方案如下所示:
中压配电网络一般采用三相三线制,闭环设计,开环运行,为辐射状结构,基本单元为馈线,配电网三相潮流计算可以以馈线为基本单元。馈线首端为变电站的中压母线,在三相潮流计算中,作为三相对称的电压源节点,电源节点A、B、C三相相对中性点的电位分别为和电源节点三相电压和为0,电源中性点电位为0。对于三相线路,由于互阻抗相对自阻抗很小,本发明只考虑线路的自阻抗,支路i三相各相的阻抗分别为ZAi,ZBi、ZCi。
本发明提供的基于中性点偏移的三相潮流计算方法流程图如图2所示,包括下述步骤:
A、获取中压配电网网络拓扑,选取中压馈线作为计算单元,对馈线上各个节点和各条支路(馈线采用节点-支路模型,馈线上每一个需要计算电压电流值的点为一个节点,两节点之间的线路为支路)进行编号;
B、引入中性点电压偏移量:
C、采用所述中性点偏移量对所述节点的计算相电压值进行修正:利用中性点电压偏移量对各相相电压进行修正,各相计算相电压值为:
本发明提供的基于中性点偏移的各相相电压示意图如图1所示。
D、采用所述修正后的计算相电压值计算各节点处的负荷在该节点处三相上产生的负荷电流值;
配变低压侧各相负荷不平衡,每相负荷都会使高压侧三线电路上分别产生负荷电流,高压侧A相产生的负荷电流是低压侧三相负荷的叠加值,B相、C相上产生的负荷电流也一样,是叠加值。节点i三相负荷分别为Sai、Sbi、Sci时,各相上的负荷电流分别为:
即:
这一等式既适应于三角形接线,又适于星形接线,与配电变压器的接线方式无关。星形接线时,配电变压器三相绕组中性点即为电压的中性点,其电位与三相电压和的电压偏移量大小相等,而三角形接线方式无实际的中性点,可将三相电压和和电压偏移量等效成一虚拟中性点的电位,这样计算过程中只需考虑电压偏移量,无需考虑接线方式。
E、采用前推回代法计算节点的三相支路电流及三相电压值:
已知首端节点电压值和各节点负荷大小,首先假设各节点各相电压初值为首端节点电压值,中性点电压偏移量设为0,利用修正后的各相计算相电压值计算各节点各相的负荷电流,计算公式为:
其中M为与节点i直接相连的所有下层支路集。对于辐射状馈线,上下层是按电流流向划分的,对于节点i和节点j,若实际电源是从i流向j,则i为j的上层节点,若是从j流向i,则i为j的下层节点。
没有下层节点的节点为末端节点,与末端支路直接相连的支路为末端支路,对于末端支路有:
利用前推过程所得的各条支路上各相的电流值,采用回代法计算各节点各相电压:
其中:节点h为与节点直接相连的上层节点;分别为节点h在A、B、C三相电压值;
对于首端节点有:
F、采用牛顿-拉夫逊法进行迭代,修正中性点电压偏移量:
采用牛顿-拉夫逊法,可得中性点电压偏移量的迭代格式:
G、采用混合迭代法对三相电压进行求解:
利用中性点偏移量对各相相电压进行修正后,采用一次前推回代法可得各相相电压值,利用牛顿-拉夫逊迭代法可修正中性点偏移量,使其逼近真值,这样经过数次混合迭代即可得各相相电压。
H、判断步骤G中的三相电压值是否满足迭代终止条件:
迭代终止条件用下述式⑩表示:
其中:ε表示无穷小的数;分别表示第k次、第k-1次迭代过程中节点i处A、B、C三相的计算相电压值。
即前后两次迭代所得的各节点各相相电压值及中性点电压偏移量均小于一个足够小的数ε,即各相相电压值已非常接近真值。
计算流程如图2所示。前推回代法及牛顿-拉夫逊法均收敛,本方法也是收敛的,且收敛速度较快。
I、计算三相潮流分布并输出结果。
所述的三相潮流计算方法用于负载接入方式为星形或三角形的配电网三相潮流计算。
所述的三相潮流计算方法用于辐射状配电网或少环配电网三相潮流计算。
本发明提供的基于中性点偏移的配电网三相潮流计算方法,根据配电网三相不平衡时三相相电压和不为零的特征,在已知A、B、C三相各相负荷的条件下,在逐次叠代潮流计算过程中,利用中性点电压偏移量,每次叠代中采用计算节点处的中性点电压偏移量修正该节点的各相计算相电压值,并采用修正后的计算相电压值计算该节点处的负荷电流,得到负荷电流后采用前推回代与牛拉法相结合的混合迭代法求解各相相电压,经数次迭代收敛最终得到潮流分布。
最后应当说明的是:以上实施例仅用以说明本发明的技术方案而非对其限制,尽管参照上述实施例对本发明进行了详细的说明,所属领域的普通技术人员应当理解:依然可以对本发明的具体实施方式进行修改或者等同替换,而未脱离本发明精神和范围的任何修改或者等同替换,其均应涵盖在本发明的权利要求范围当中。
Claims (11)
1.一种基于中性点偏移的配电网三相潮流计算方法,其特征在于,所述方法包括下述步骤:
A、获取中压配电网网络拓扑,选取中压馈线作为计算单元,对馈线上各个节点和各条支路进行编号;
B、引入中性点电压偏移量;
C、采用所述中性点偏移量对所述节点的计算相电压值进行修正;
D、采用所述修正后的计算相电压值计算节点处的负荷电流值;
E、采用前推回代法计算各支路的三相电流值及各节点的三相电压值;
F、采用牛顿-拉夫逊法进行迭代,修正中性点电压偏移量;
G、采用混合迭代法对三相电压进行求解;
H、判断步骤G中的三相电压值是否满足迭代终止条件;
I、计算三相潮流分布并输出结果。
2.如权利要求1所述的配电网三相潮流计算方法,其特征在于,所述步骤A中,所述中压配电网采用辐射状配电网或少环状配电网,在三相潮流计算中,选用中压馈线作为基本计算单元,馈线首端的中压母线作为三相对称的电压源;并对馈线上各支路和各节点进行编号。
3.如权利要求1所述的配电网三相潮流计算方法,其特征在于,所述步骤B中,当配电网线路参数、负荷不平衡时,电源中性点电位不为零;所述中性点电压偏移量用表示。
9.如权利要求1所述的配电网三相潮流计算方法,其特征在于,所述步骤G中,利用所述中性点偏移量对三相相电压进行修正后,采用一次前推回代法得三相相电压值,利用牛顿-拉夫逊迭代法修正中性点偏移量,经过混合迭代即得三相相电压。
11.如权利要求1所述的配电网三相潮流计算方法,其特征在于,所述方法的负载接入方式为星形或三角形。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210392893.XA CN103001214B (zh) | 2012-10-16 | 2012-10-16 | 一种基于中性点偏移的配电网三相潮流计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210392893.XA CN103001214B (zh) | 2012-10-16 | 2012-10-16 | 一种基于中性点偏移的配电网三相潮流计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103001214A true CN103001214A (zh) | 2013-03-27 |
CN103001214B CN103001214B (zh) | 2015-08-05 |
Family
ID=47929484
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210392893.XA Active CN103001214B (zh) | 2012-10-16 | 2012-10-16 | 一种基于中性点偏移的配电网三相潮流计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103001214B (zh) |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103326354A (zh) * | 2013-05-23 | 2013-09-25 | 国家电网公司 | 基于负荷矩的闭环配电网潮流估计方法 |
CN103353918A (zh) * | 2013-04-24 | 2013-10-16 | 国家电网公司 | 一种配电网三相牛顿拉夫逊法潮流计算初值选取方法 |
CN103441518A (zh) * | 2013-08-06 | 2013-12-11 | 国家电网公司 | 单相负荷和相间负荷混合情况下的三相潮流分布确定方法 |
CN103825270A (zh) * | 2014-03-04 | 2014-05-28 | 国家电网公司 | 一种配电网三相状态估计雅可比矩阵常数化的处理方法 |
CN104269872A (zh) * | 2014-10-13 | 2015-01-07 | 国家电网公司 | 一种变压器三相节点导纳矩阵奇异的处理方法 |
CN104793509A (zh) * | 2015-01-30 | 2015-07-22 | 国家电网公司 | 一种时域仿真中pq型分布式电源与配电网接口处理方法 |
CN106339811A (zh) * | 2016-08-26 | 2017-01-18 | 中国电力科学研究院 | 一种低压配电网精确线损分析方法 |
CN107328981A (zh) * | 2017-08-11 | 2017-11-07 | 华北电力大学 | 一种变压器中性点故障电压的分析方法 |
CN107565515A (zh) * | 2017-08-17 | 2018-01-09 | 国家电网公司 | 一种采用角度补偿电压幅值的变压器过激磁保护方法 |
CN111478335A (zh) * | 2020-05-07 | 2020-07-31 | 国网山东省电力公司经济技术研究院 | 一种计及分布式光伏的配电网潮流计算方法及系统 |
CN114460409A (zh) * | 2022-04-13 | 2022-05-10 | 广东电网有限责任公司佛山供电局 | 基于中性点不平衡数据的电容器监测方法、装置及设备 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH0923578A (ja) * | 1995-07-03 | 1997-01-21 | Jinbatsu Shiyuukichi | 単相三線式配電需用家の構内配電改良方法 |
CN102636731A (zh) * | 2012-04-12 | 2012-08-15 | 上海交通大学 | 基于负荷测录系统的配电网单相断线判定及寻址方法 |
-
2012
- 2012-10-16 CN CN201210392893.XA patent/CN103001214B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH0923578A (ja) * | 1995-07-03 | 1997-01-21 | Jinbatsu Shiyuukichi | 単相三線式配電需用家の構内配電改良方法 |
CN102636731A (zh) * | 2012-04-12 | 2012-08-15 | 上海交通大学 | 基于负荷测录系统的配电网单相断线判定及寻址方法 |
Non-Patent Citations (1)
Title |
---|
张辉: "三相四线制配电网络潮流模型与算法研究", 《中国优秀硕士学位论文全文数据库》, no. 5, 15 November 2007 (2007-11-15) * |
Cited By (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103353918A (zh) * | 2013-04-24 | 2013-10-16 | 国家电网公司 | 一种配电网三相牛顿拉夫逊法潮流计算初值选取方法 |
CN103353918B (zh) * | 2013-04-24 | 2016-08-31 | 国家电网公司 | 一种配电网三相牛顿拉夫逊法潮流计算初值选取方法 |
CN103326354A (zh) * | 2013-05-23 | 2013-09-25 | 国家电网公司 | 基于负荷矩的闭环配电网潮流估计方法 |
CN103326354B (zh) * | 2013-05-23 | 2015-03-25 | 国家电网公司 | 基于负荷矩的闭环配电网潮流估计方法 |
CN103441518A (zh) * | 2013-08-06 | 2013-12-11 | 国家电网公司 | 单相负荷和相间负荷混合情况下的三相潮流分布确定方法 |
CN103441518B (zh) * | 2013-08-06 | 2016-04-20 | 国家电网公司 | 单相负荷和相间负荷混合情况下的三相潮流分布确定方法 |
CN103825270B (zh) * | 2014-03-04 | 2016-08-17 | 国家电网公司 | 一种配电网三相状态估计雅可比矩阵常数化的处理方法 |
CN103825270A (zh) * | 2014-03-04 | 2014-05-28 | 国家电网公司 | 一种配电网三相状态估计雅可比矩阵常数化的处理方法 |
CN104269872A (zh) * | 2014-10-13 | 2015-01-07 | 国家电网公司 | 一种变压器三相节点导纳矩阵奇异的处理方法 |
CN104793509A (zh) * | 2015-01-30 | 2015-07-22 | 国家电网公司 | 一种时域仿真中pq型分布式电源与配电网接口处理方法 |
CN106339811A (zh) * | 2016-08-26 | 2017-01-18 | 中国电力科学研究院 | 一种低压配电网精确线损分析方法 |
CN107328981A (zh) * | 2017-08-11 | 2017-11-07 | 华北电力大学 | 一种变压器中性点故障电压的分析方法 |
CN107328981B (zh) * | 2017-08-11 | 2020-01-10 | 华北电力大学 | 一种变压器中性点故障电压的分析方法 |
CN107565515A (zh) * | 2017-08-17 | 2018-01-09 | 国家电网公司 | 一种采用角度补偿电压幅值的变压器过激磁保护方法 |
CN111478335A (zh) * | 2020-05-07 | 2020-07-31 | 国网山东省电力公司经济技术研究院 | 一种计及分布式光伏的配电网潮流计算方法及系统 |
CN111478335B (zh) * | 2020-05-07 | 2021-07-13 | 国网山东省电力公司经济技术研究院 | 一种计及分布式光伏的配电网潮流计算方法及系统 |
CN114460409A (zh) * | 2022-04-13 | 2022-05-10 | 广东电网有限责任公司佛山供电局 | 基于中性点不平衡数据的电容器监测方法、装置及设备 |
Also Published As
Publication number | Publication date |
---|---|
CN103001214B (zh) | 2015-08-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103001214A (zh) | 一种基于中性点偏移的配电网三相潮流计算方法 | |
CN104362622B (zh) | 一种带小水电机组地区电网的负荷模型和建模方法 | |
CN103606921B (zh) | 基于实时广义戴维南等值的暂态电压稳定判断方法 | |
CN102104252B (zh) | 一种适用于电磁暂态分析的电力系统动态等值方法 | |
CN107064736B (zh) | 一种含多t接逆变型分布式电源配电网的故障定位方法 | |
CN102403716A (zh) | 一种多馈入交直流电网动态等值的方法 | |
CN103595040A (zh) | 含直流分布式电源的配电网及其综合负荷特性的分析方法 | |
CN102545252B (zh) | 基于三阶收敛牛顿法vsc-hvdc潮流计算方法 | |
CN103984822A (zh) | 一种基于全数字实时仿真装置的三相潮流实现方法 | |
CN103116704A (zh) | 一种基于局部几何参数化的连续潮流计算方法 | |
CN107666155A (zh) | 基于Markov模型的多能互补系统随机稳定性分析方法 | |
CN104917174A (zh) | 一种针对主动配电网的静态电压稳定性判定方法 | |
CN105762777A (zh) | 一种含多t接逆变型分布式电源配电网的纵联保护方法 | |
CN103199528A (zh) | 广域电力系统状态估计协调方法 | |
CN105071447A (zh) | 一种交直流混联多馈入直流受端电网强度获得方法 | |
CN103955594A (zh) | 一种电力系统动态等值方法 | |
CN106291046A (zh) | 混压同塔双回线单相跨单相的跨电压故障电流计算方法 | |
CN106877309A (zh) | 主从控制模式下直流微电网的稳定性判定方法 | |
CN103825270B (zh) | 一种配电网三相状态估计雅可比矩阵常数化的处理方法 | |
CN104268798B (zh) | 一种分布式电源及电动汽车对配电网影响的评估方法 | |
CN104283208B (zh) | 大型电网概率性可用输电能力的分解协调计算方法 | |
CN107196307A (zh) | 一种变压器故障后电网有功潮流快速估算的方法 | |
CN105514996B (zh) | 一种主配网协同潮流计算的方法和系统 | |
CN106340906A (zh) | 基于改进遗传算法的交直流系统低压减载优化方法 | |
CN103390890B (zh) | 基于电流分布系数的配电网潮流分析方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |