CN112858979B - 基于三维大气电场测量的雷暴云点电荷定位海拔校正方法 - Google Patents

基于三维大气电场测量的雷暴云点电荷定位海拔校正方法 Download PDF

Info

Publication number
CN112858979B
CN112858979B CN202110036291.XA CN202110036291A CN112858979B CN 112858979 B CN112858979 B CN 112858979B CN 202110036291 A CN202110036291 A CN 202110036291A CN 112858979 B CN112858979 B CN 112858979B
Authority
CN
China
Prior art keywords
electric field
correction
point
altitude
charge
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
CN202110036291.XA
Other languages
English (en)
Other versions
CN112858979A (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.)
Nanjing University of Information Science and Technology
Original Assignee
Nanjing University of Information Science and Technology
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 Nanjing University of Information Science and Technology filed Critical Nanjing University of Information Science and Technology
Priority to CN202110036291.XA priority Critical patent/CN112858979B/zh
Publication of CN112858979A publication Critical patent/CN112858979A/zh
Application granted granted Critical
Publication of CN112858979B publication Critical patent/CN112858979B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R35/00Testing or calibrating of apparatus covered by the other groups of this subclass
    • G01R35/005Calibrating; Standards or reference devices, e.g. voltage or resistance standards, "golden" references
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R29/00Arrangements for measuring or indicating electric quantities not covered by groups G01R19/00 - G01R27/00
    • G01R29/08Measuring electromagnetic field characteristics
    • G01R29/0807Measuring electromagnetic field characteristics characterised by the application
    • G01R29/0814Field measurements related to measuring influence on or from apparatus, components or humans, e.g. in ESD, EMI, EMC, EMP testing, measuring radiation leakage; detecting presence of micro- or radiowave emitters; dosimetry; testing shielding; measurements related to lightning
    • G01R29/0842Measurements related to lightning, e.g. measuring electric disturbances, warning systems
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R29/00Arrangements for measuring or indicating electric quantities not covered by groups G01R19/00 - G01R27/00
    • G01R29/08Measuring electromagnetic field characteristics
    • G01R29/0864Measuring electromagnetic field characteristics characterised by constructional or functional features
    • G01R29/0892Details related to signal analysis or treatment; presenting results, e.g. displays; measuring specific signal features other than field strength, e.g. polarisation, field modes, phase, envelope, maximum value
    • 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
    • Y02ATECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/10Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Electromagnetism (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

本发明公开了一种基于三维大气电场测量的雷暴云点电荷定位海拔校正方法,包括以下步骤:首先建立雷暴云点电荷定位海拔校正模型,定义校正参数;然后利用三维大气电场仪,测得三维大气电场分量;最后推导出海拔校正后的点电荷坐标计算公式。基于间接测量误差理论,研究大气电场仪所在位置的海拔高度对雷暴云点电荷定位性能的影响,对校正前后的定位误差进行比较分析。本发明在雷暴监测领域可以显著降低了在单个大气电场仪下获取雷暴云方位数据的误差,通过数据证明,相比于校正前,该方法在实际应用中具有更好的性能,减少了海拔对雷暴云点电荷定位性能的负面影响,为雷暴云预报预警研究提供了一种新的误差补偿手段。

Description

基于三维大气电场测量的雷暴云点电荷定位海拔校正方法
技术领域
本发明涉及一种雷电探测方法,尤其涉及一种基于三维大气电场测量的雷暴云点电荷定位海拔校正方法。
背景技术
在日常生活中,避雷针可以有效减少了雷电灾害带来的经济损失,但这些属于被动防雷范畴,无法准确预报雷电发生的时间和位置。近些年雷电预警技术得到了快速发展,通过分析晴天和雷暴天气时三维大气电场的时域特征,为雷暴发生阶段追踪雷暴云运动轨迹提供了一种思路,然而,现有的研究在促进基于三维大气电场测量的雷暴云探测工作的同时,忽略了大气电场仪所在位置的海拔高度对雷暴云点电荷定位精度的负面影响,在预测精度上面难有显著的突破。
发明内容
发明目的:本发明的目的在于提供一种显著降低雷暴云方位数据误差的基于三维大气电场测量的雷暴云点电荷定位海拔校正方法。
技术方案:本发明的基于三维大气电场测量的雷暴云点电荷定位海拔校正方法,包括以下步骤:
步骤一、建立雷暴云点电荷定位海拔校正模型,以N点为坐标原点,建立三维直角坐标系,其中,S(x,y,z)为海拔校正前雷暴云点电荷所在位置;N(0,0,0)为三维大气电场仪观测点所在位置;海拔校正点M(0,0,-h)为电场仪所在Z轴与海平面的交点;h表示大气电场仪自身高度与其所处位置的海拔之和,海拔校正前的雷暴云点电荷水平偏角和仰角分别为α,β;r为点电荷S到电场仪N的距离;
步骤二、大气电场仪在N处,测得S的电场强度为E,由大气电场仪测得x、y、z方向上校准前的电场分量分别为Ex、Ey、Ez
步骤三、在测点M处,得到S的电场强度为E';对电场分量Ex、Ey、Ez进行校准,校准后的电场分量分别为Ex'、Ey'、Ez',校正后的点电荷水平偏角和仰角分别为α',β',r'为S到M的距离;
步骤四、海拔校正后的点电荷S的球坐标为S'(r',α',β'),根据三维直角坐标系的矢量关系得到校正后的点电荷直角坐标S'(x',y',z');
步骤五、通过仿真对比分析海拔校正前后的测距误差。
进一步地,步骤一中,所述海拔校正前雷暴云点电荷所在位置S(x,y,z)与α、β、r的关系分别为:
Figure GDA0003695801090000021
Figure GDA0003695801090000022
其中,
Figure GDA0003695801090000023
进一步地,步骤三中,所述校正后的三维电场分量Ex'、Ey'、Ez'分别为:
Figure GDA0003695801090000024
由式(3)进一步得到:
Figure GDA0003695801090000025
其中,
Figure GDA0003695801090000026
进一步地,步骤四中,海拔校正后的点电荷球坐标S'(r',α',β')为:
Figure GDA0003695801090000031
其中,
Figure GDA0003695801090000032
进一步地,步骤四中,根据矢量关系,校正后的点电荷直角坐标S'(x',y',z')为:
Figure GDA0003695801090000033
有益效果:与现有技术相比,本发明具有如下显著优点:
本发明提出了一种基于三维大气电场测量的点电荷定位海拔校正方法,在雷暴监测领域可以显著降低了在单个大气电场仪下获取雷暴云方位数据的误差,通过数据证明,相比于校正前,该方法在实际应用中具有更好的性能,减少了海拔对雷暴云点电荷定位性能的负面影响,为雷暴云预报预警研究提供了一种新的误差补偿手段。
附图说明
图1为本发明雷暴云点电荷定位海拔校正模型图;
图2为本发明雷暴云点电荷到校正点距离、校正高度与测距误差的关系图;
图3为本发明海拔校正前后雷暴云点电荷定位测距误差对比图;
图4为本发明雷暴云点电荷到校正点距离、校正高度与水平偏角测量误差的关系图;
图5为本发明海拔校正前后雷暴云点电荷定位水平偏角测量误差对比图;
图6为本发明雷暴云点电荷到校正点距离、校正高度与仰角测量误差的关系图;
图7为本发明海拔校正前后雷暴云点电荷定位仰角测量误差对比图。
具体实施方式
下面结合附图对本发明的技术方案作进一步说明。
本发明提供一种基于三维大气电场测量的雷暴云点电荷定位海拔校正方法,包括如下步骤:
一、建立雷暴云点电荷定位海拔校正模型,其具体为:
如图1所示,建立的雷暴云点电荷定位海拔校正模型中,以N点为坐标原点,建立三维直角坐标系,其中:S(x,y,z)为海拔校正前雷暴云点电荷所在位置;N(0,0,0)为三维大气电场仪观测点所在位置;海拔校正点M(0,0,-h)为电场仪所在Z轴与海平面的交点;h表示大气电场仪自身高度与其所处位置的海拔之和,记为校正高度h。海拔校正前的雷暴云点电荷水平偏角和仰角分别为α,β;r为点电荷S到电场仪N的距离;在大气电场仪N处,测得S的电场强度为E,由大气电场仪测得x、y、z方向上校准前的电场分量分别为Ex、Ey、Ez,校准后的电场分量分别为Ex'、Ey'、Ez'。校正后的点电荷水平偏角和仰角分别为α',β';r'为S到M的距离;在测点M处,得到S的电场强度为E';海拔校正后的点电荷S的球坐标为S'(r',α',β')。
二、提出一种基于三维大气电场测量的雷暴云点电荷定位海拔校正方法,其具体为:
在建立雷暴云电场测量模型的基础上,利用镜像法得到了未进行海拔校正的三维大气电场分量与雷暴云点电荷球坐标(r,α,β)、直角坐标(x,y,z)的关系分别为:
Figure GDA0003695801090000041
Figure GDA0003695801090000042
其中,
Figure GDA0003695801090000043
基于图1所示模型,得到校正后的三维电场分量Ex'、Ey'、Ez'分别为:
Figure GDA0003695801090000044
由式(3)进一步得到:
Figure GDA0003695801090000051
其中,
Figure GDA0003695801090000052
利用式(4),得到海拔校正后的点电荷球坐标S'(r',α',β')为:
Figure GDA0003695801090000053
根据图1所示矢量关系,校正后的点电荷直角坐标S'(x',y',z')为:
Figure GDA0003695801090000054
根据上述方法,以下对基于三维大气电场测量的雷暴云点电荷定位海拔校正方法从测距测向方面进行性能分析:
基于间接测量误差理论,分析大气电场仪所在位置的海拔高度对雷暴云点电荷定位测距测向精度的影响,进一步研究点电荷定位性能。另外,设大气电场分量测量标准偏差均为
Figure GDA0003695801090000055
利用雷暴云点电荷S与校正点M的距离r',校正后的雷暴云水平偏角α'、仰角β'的计算公式,由电场分量测量误差
Figure GDA0003695801090000056
引起距离r',水平偏角α',仰角β'的测量误差σr,σα,σβ为:
Figure GDA0003695801090000057
由式(7)可知,测距误差σr与距离r'、高度h、仰角β'和电场测量误差
Figure GDA0003695801090000058
均有关,仅与水平偏角α'无关;同样地,测向误差σα,σβ也是如此,仅与水平偏角α'无关;雷暴云点电荷S到校正点M的距离r'、校正高度h和电场测量误差
Figure GDA0003695801090000061
越小,点电荷定位测距测向精度越高。对式(7)的分析,可初步判断,高度h与测距测向误差存在相关性,未考虑高度h进行定位性能分析是不准确的。特别地,结合式(7),较大的高度h会对定位性能产生负面影响。
1、雷暴云点电荷定位测距性能分析
利用式(7),研究距离,校正高度与测距误差的关系,仿真结果如图2所示。图2中,测距误差σr均随着距离r'和校正高度h的增大而增大;当距离r'在0到0.6km时,误差σr受高度h的影响较小,σr几乎为零;当距离r'大于0.6km时,误差σr随着高度h的增大而缓慢增大,最大达到0.061km。无论高度h取值如何,误差σr均会随着距离r'的增加而增大。
为研究高度对点电荷定位测距性能的影响,通过仿真对比分析海拔校正前后的测距误差,结果如图3所示。从图3看出,当校正高度一定时,测距误差σr随距离r'的增大而增大,特别是海拔校正前的测距误差变化地更为剧烈,最大误差达到0.053km。能够发现,虽然校正后的测距误差σr仍存在,但是其随距离r'的变化较为缓和,且相比于校正前,由间接测量带来的测距误差更小。特别是,当校正高度越小时,测距误差σr也越小,最大高度h为0.1km时的误差仅为0.038km。因此,在实际实验中,将三维大气电场仪置于较空旷的低海拔处进行雷暴云探测,会有更好的点电荷定位测距效果。
2、雷暴云点电荷定位测向性能分析
研究距离r',校正高度h与水平偏角测量误差σα的关系,仿真结果如图4所示。图4中,随着距离r'和校正高度h的增大,水平偏角测量误差σα也随之增大。当距离r'在0到0.6km时,误差σα受高度h影响较小,误差σα小于0.005度;当距离r'大于0.6km时,误差σα随高度h的增加而缓慢增加,最大达到0.061度。当高度h处于0到0.05km时,误差σα几乎不受距离r'影响,σα小于0.038度;而当距离r'和高度h均较大时,误差σα随着两者的增大而骤升至0.061度。
为进一步研究高度h对水平偏角测量的影响,通过仿真对比分析海拔校正前后的水平偏角误差,结果如图5所示。图5中,海拔校正前后的水平偏角测量误差σα均随距离r'的增大而增大,特别是当距离r'大于1km时,海拔校正前的误差急剧增大,最大骤升至0.093度。相反地,随着距离r'的增大,校正后的误差增加的较缓,校正高度h为0.1km时的最大误差仅为0.038度。此外,不难发现,校正高度越小,海拔校正方法对测向性能的积极影响越大。
同样地,研究距离r',校正高度h与仰角测量误差σβ的关系,仿真结果如图6所示。图6中,仰角测量误差σβ均随距离r'和校正高度h的增大而增大;当距离r'在0到1km时,误差σβ受高度h影响较小,误差σβ小于0.001度;而当距离r'大于1km时,随着高度h的增大,误差σβ爬升速度随h的增大而变快,特别当h取到0.5km时,σβ达到最大值0.0087度。此外,当高度h小于0.1km时,误差σβ几乎不受距离r'影响,σβ小于0.0013度,而当h大于0.1km时,随着距离r'的增大,误差σβ也会增大,特别是h越大导致该变化越剧烈。
进一步地,通过仿真对比分析海拔校正前后的仰角误差,结果如图7所示。图7中,随着距离r'的不断增大,校正前的仰角测量误差σβ增加的愈加剧烈,当r'增加至2km时,误差σβ爬升至0.0425度。相比之下,海拔校正方法表现出了很好的校正效果。校正后的误差σβ几乎不随r'的变化而变化。当r'大于1.2km时,相比于0.02km的校正高度,进一步增大高度h会引起误差σβ的小幅增大,但最大误差仅为0.001度,仍处于较小的误差水平。
综上所述,海拔校正前后的水平偏角测量误差和仰角测量误差均存在,但校正后的误差得到了有效降低,特别是对减小仰角误差具有更好的效果。结合图1可知,海拔校正与否主要影响雷暴云点电荷到校正点距离和仰角的测量,引入海拔校正方法后,能较好地提升雷暴云点电荷定位测距测向性能。

Claims (1)

1.一种基于三维大气电场测量的雷暴云点电荷定位海拔校正方法,其特征在于,包括以下步骤:
步骤一、建立雷暴云点电荷定位海拔校正模型,以N点为坐标原点,建立三维直角坐标系,其中,S(x,y,z)为海拔校正前雷暴云点电荷所在位置;N(0,0,0)为三维大气电场仪观测点所在位置;海拔校正点M(0,0,-h)为大气电场仪所在Z轴与海平面的交点;h表示大气电场仪自身高度与其所处位置的海拔之和,海拔校正前的雷暴云点电荷水平偏角和仰角分别为α,β;r为点电荷S到大气电场仪N的距离;
步骤二、大气电场仪在N处,测得S的电场强度为E,由大气电场仪测得x、y、z方向上校正前的电场分量分别为Ex、Ey、Ez
步骤三、在海拔校正点M处,得到S的电场强度为E',对电场分量Ex、Ey、Ez进行校正,校正后的电场分量分别为Ex'、Ey'、Ez',校正后的点电荷水平偏角和仰角分别为α'、β',r'为S到M的距离;
步骤四、海拔校正后的点电荷在S处的球坐标为S'(r',α',β'),根据三维直角坐标系的矢量关系得到校正后的点电荷直角坐标S'(x',y',z');
步骤五、通过仿真对比分析海拔校正前后的测距误差;
步骤一中,所述海拔校正前雷暴云点电荷所在位置S(x,y,z)与α、β、r的关系分别为:
Figure FDA0003695801080000011
Figure FDA0003695801080000012
其中,
Figure FDA0003695801080000013
步骤三中,所述校正后的三维电场分量Ex'、Ey'、Ez'分别为:
Figure FDA0003695801080000021
由式(3)进一步得到:
Figure FDA0003695801080000022
其中,
Figure FDA0003695801080000023
步骤四中,海拔校正后的点电荷球坐标S'(r',α',β')为:
Figure FDA0003695801080000024
其中,
Figure FDA0003695801080000025
根据矢量关系,校正后的点电荷直角坐标S'(x',y',z')为:
Figure FDA0003695801080000026
CN202110036291.XA 2021-01-12 2021-01-12 基于三维大气电场测量的雷暴云点电荷定位海拔校正方法 Active CN112858979B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110036291.XA CN112858979B (zh) 2021-01-12 2021-01-12 基于三维大气电场测量的雷暴云点电荷定位海拔校正方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110036291.XA CN112858979B (zh) 2021-01-12 2021-01-12 基于三维大气电场测量的雷暴云点电荷定位海拔校正方法

Publications (2)

Publication Number Publication Date
CN112858979A CN112858979A (zh) 2021-05-28
CN112858979B true CN112858979B (zh) 2022-08-23

Family

ID=76002846

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110036291.XA Active CN112858979B (zh) 2021-01-12 2021-01-12 基于三维大气电场测量的雷暴云点电荷定位海拔校正方法

Country Status (1)

Country Link
CN (1) CN112858979B (zh)

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105372724A (zh) * 2015-12-10 2016-03-02 国网电力科学研究院武汉南瑞有限责任公司 一种针对山区水域障碍物大气电场仪观测结果的修订方法
JP2017009512A (ja) * 2015-06-25 2017-01-12 Kddi株式会社 気圧値を用いて高度及び誤差を推定する高度推定装置、プログラム及び方法
CN109142896A (zh) * 2018-07-25 2019-01-04 南京信息工程大学 基于三维大气电场和memd的雷电预警方法
CN110174557A (zh) * 2019-01-11 2019-08-27 南京信息工程大学 一种基于三维大气电场仪观测视角的雷暴云定位校准方法

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9097792B2 (en) * 2009-08-12 2015-08-04 The Johns Hopkins University System and method for atmospheric correction of information
CN106597417A (zh) * 2017-01-10 2017-04-26 北京航天计量测试技术研究所 一种远距离扫描激光雷达测量误差的修正方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2017009512A (ja) * 2015-06-25 2017-01-12 Kddi株式会社 気圧値を用いて高度及び誤差を推定する高度推定装置、プログラム及び方法
CN105372724A (zh) * 2015-12-10 2016-03-02 国网电力科学研究院武汉南瑞有限责任公司 一种针对山区水域障碍物大气电场仪观测结果的修订方法
CN109142896A (zh) * 2018-07-25 2019-01-04 南京信息工程大学 基于三维大气电场和memd的雷电预警方法
CN110174557A (zh) * 2019-01-11 2019-08-27 南京信息工程大学 一种基于三维大气电场仪观测视角的雷暴云定位校准方法

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
Hong-yan Xing等.Analysis on Electric Field Based on Three Dimensional Atmospheric Electric Field Apparatus.《J Electr Eng Technol》.2018,第13卷(第4期), *
Hongyan Xing等.Thunderstorm Cloud Localization Algorithm and Performance Analysis of a Three‑Dimensional Atmospheric Electric Field Apparatus.《Journal of Electrical Engineering & Technolog》.2019, *
Wei Xu等.Inversion of a Thunderstorm Cloud Charging Model Based on a 3D Atmospheric Electric Field.《Applied sciences》.2018, *
周俊驰等.大气电场仪观测结果的修订.《南京信息工程大学学报·自然科学版》.2011,第3卷(第3期), *
行鸿彦等.大气电场仪观测数据的海拔校正及联网.《解放军理工大学学报(自然科学版)》.2014,第15卷(第6期), *

Also Published As

Publication number Publication date
CN112858979A (zh) 2021-05-28

Similar Documents

Publication Publication Date Title
CN109001711B (zh) 多线激光雷达标定方法
CN112859051B (zh) 激光雷达点云运动畸变的矫正方法
KR101914550B1 (ko) 레이더의 표적 위치 추적 방법
CN110490932B (zh) 单目红外共面光标迭代优化测量行车吊臂空间位姿方法
CN110081905B (zh) 一种基于单站光电经纬仪的光波大气折射误差计算方法
CN107843895B (zh) 一种双多普勒雷达三维风场反演方法
CN109883320B (zh) 一种土地面积测量方法及系统
CN109917151B (zh) 一种基于动态权值分配的船舶相对风融合算法
CN113129377B (zh) 一种三维激光雷达快速鲁棒slam方法和装置
CN112270320A (zh) 一种基于卫星影像校正的输电线路杆塔坐标校准方法
CN110927765B (zh) 激光雷达与卫星导航融合的目标在线定位方法
JP7529337B1 (ja) 非等方性を考慮した水平方向対流圏遅延分類方法
CN115683170A (zh) 基于雷达点云数据融合误差的校准方法
CN115824223A (zh) 一种基于多源融合的室内外无缝定位方法
CN118010000A (zh) 一种基于激光点云的高压电塔垂直度检测方法
CN112858979B (zh) 基于三维大气电场测量的雷暴云点电荷定位海拔校正方法
CN108562870B (zh) 一种声源定位校准方法
CN112766127B (zh) 基于互补集合模态分解和sg滤波的雷云点电荷定位方法
CN111679261B (zh) 一种基于反光板的激光雷达定位方法及系统
CN108562269A (zh) 一种相对高度测量方法及装置
CN116973895A (zh) 一种激光点云匹配的实时位姿校正方法
CN104075710A (zh) 一种基于航迹预测的机动扩展目标轴向姿态实时估计方法
CN108508444A (zh) 一种采动影响区输电线路弧垂计算方法及系统
CN111609850B (zh) 一种用于绳索悬挂式机器人定位补偿方法
CN115201823A (zh) 一种利用BDS-InSAR数据融合的地表形变监测方法

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