CN111065106B - 移动通信网络中基于异常检测和核密度估计kde的指标突变小区检测方法 - Google Patents

移动通信网络中基于异常检测和核密度估计kde的指标突变小区检测方法 Download PDF

Info

Publication number
CN111065106B
CN111065106B CN201911402828.9A CN201911402828A CN111065106B CN 111065106 B CN111065106 B CN 111065106B CN 201911402828 A CN201911402828 A CN 201911402828A CN 111065106 B CN111065106 B CN 111065106B
Authority
CN
China
Prior art keywords
kpi
abnormal
cell
sequence
kde
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
CN201911402828.9A
Other languages
English (en)
Other versions
CN111065106A (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.)
Chongqing University
Original Assignee
Chongqing University
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 Chongqing University filed Critical Chongqing University
Priority to CN201911402828.9A priority Critical patent/CN111065106B/zh
Publication of CN111065106A publication Critical patent/CN111065106A/zh
Application granted granted Critical
Publication of CN111065106B publication Critical patent/CN111065106B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04WWIRELESS COMMUNICATION NETWORKS
    • H04W16/00Network planning, e.g. coverage or traffic planning tools; Network deployment, e.g. resource partitioning or cells structures
    • H04W16/22Traffic simulation tools or models
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04WWIRELESS COMMUNICATION NETWORKS
    • H04W24/00Supervisory, monitoring or testing arrangements
    • H04W24/02Arrangements for optimising operational condition
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04WWIRELESS COMMUNICATION NETWORKS
    • H04W24/00Supervisory, monitoring or testing arrangements
    • H04W24/04Arrangements for maintaining operational condition

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Mobile Radio Communication Systems (AREA)
  • Information Retrieval, Db Structures And Fs Structures Therefor (AREA)

Abstract

本发明公开了移动通信网络中基于异常检测和核密度估计KDE的指标突变小区检测方法,主要步骤为:1)建立历史数据序列和测试数据序列。2)筛选掉全部KPI序列均正常的小区,保留异常小区。3)利用基于高斯核的KDE计算异常小区差分前的异常KPI序列的异常分数,并建立异常分数矩阵A。4)判断异常分数矩阵A中是否存在元素a>α,若是,则记元素a对应的KPI为突变指标,对应小区记为指标突变小区。本发明可精确定位指标突变小区,灵活性高,鲁棒性强,性能稳定。

Description

移动通信网络中基于异常检测和核密度估计KDE的指标突变 小区检测方法
技术领域
本发明涉及移动通信网络的网络优化技术,具体是移动通信网络中基于异常检测和核密度估计KDE的指标突变小区检测方法。
背景技术
在移动通信网络中,指标突变小区的检测一直是一个重要的研究内容。指标突变小区的检测应用范围十分广泛,尤其是在快速定位网络恶化、指定指标优化、指标突变小区检测等方面起着重要作用。
在长期演进中,某些小区的一些关键性能指标(KPI,Key PerformanceIndicators)数据相对长时间的历史数据发生了较大的变化,被认为相对之前的数据范围产生了较大的异常。此外,在这些异常的小区中,某些小区的KPI相对近期的KPI急剧变化。这些因为KPI急剧变化导致指标异常的小区即为指标突变小区。
目前的指标突变小区的检测主要分为两类:一类是基于长期网络优化的经验判断,即基于长期演进的经验来划定正常小区指标的门限,根据这个门限和实时的网络数据进行比较,如果某小区的指标数据超出了这个门限,则将该小区标记为指标突变小区;一类是使用实时数据进行估计,它用密度估计的方法,计算每个小区的每个KPI的突变程度,如果某个指标突变程度过大,该指标即为突变指标,该小区即为指标突变小区。
在现有的指标突变小区检测的方法中,基于经验门限的方法设置了一个固定的门限值来筛选出异常小区,在实际生产中,使用经验门限的方法要求在检测之前给定一个固定的门限范围,这个范围通常来自于长时间网络优化的经验。然而,这个门限值不能随着检测经验而快速更新,也不具有根据检测数据本身更新的特性,从而影响了异常检测的效果。此外,应当指出的是,这种基于设定好的门限检测得到的异常小区并不全是指标突变小区,还包含了很多其他造成异常的原因,因此仅仅用该检测方法是不完善的。另一方面,基于密度函数检测的方法使用每个数据点的临近点数据来估计每个数据点的突变程度。密度函数的检测需要一定量的训练数据,在现实生产中,该方法的训练数据往往来源于实时数据周边的少量数据样本点。由于这些样本点的数据并不一定都是正常数据,因此可能会对训练本身造成干扰。同时,由于这些样本点本身数量较少,直接使用这些数据作为训练数据存在训练量不足的弊端。
发明内容
本发明的目的是解决现有技术中存在的问题。
为实现本发明目的而采用的技术方案是这样的,移动通信网络中基于异常检测和核密度估计KDE的指标突变小区检测方法,主要包括以下步骤:
1)获取移动通信网络所有小区关键性能指标KPI序列,并对每个小区的KPI序列进行数据清洗,主要步骤如下:
1.1)基于设定的粒度g,筛选掉粒度不合规的小区。
1.2)基于设定的数据长度S=T/g,筛选掉时间长度不合规的小区。T为指标突变小区检测周期。
1.3)以0填充剩余小区缺失值。
1.4)基于业务告警信息,筛选掉业务告警小区。
2)基于清洗后的KPI序列,建立历史数据序列和测试数据序列。
进一步,历史数据序列为t0-t1时间段的KPI序列。测试数据序列为t2-t3时间段的KPI序列。t0为KPI序列起始时间。t1∈(t0,t3)为历史数据序列终结时间。t2∈(t0,t3)为测试数据序列起始时间。t3为KPI序列终结时间。
3)基于历史数据序列和测试数据序列,筛选掉全部KPI序列均正常的小区,保留异常小区,并记录每个小区的所有异常KPI的第一个发生异常的时间点t1st
确定异常小区的主要步骤如下:
3.1)利用季节差分法对历史数据序列和测试数据序列进行差分。将差分后的历史数据序列和测试数据序列输入到加法模型中,得到历史趋势数据集和测试趋势数据集。
加法模型Y[t]=T[t]+S[t]+e[t]。其中,T[t]、S[t]、e[t]为KPI差分后得到的三个分解因子。Y[t]表示加法模型输出的历史趋势数据或测试趋势数据。
3.2)基于历史趋势数据集,分别计算不同小区中不同KPI序列的均值μ和标准差σ,并确定每个KPI的正常区间I。
KPI的正常区间I=(μ-ω×σ,μ+ω×σ)。ω为常数。
3.3)判断测试趋势数据集中每个元素是否位于对应KPI序列的正常区间I内,若否,则将对应KPI标记为异常指标,并将所述异常指标所在小区标记为异常小区。
4)利用基于高斯核的KDE计算异常小区差分前的异常KPI序列的异常分数,并建立异常分数矩阵A。
建立异常分数矩阵A的主要步骤如下:
4.1)以t0-t1st时间段的异常KPI序列作为KDE的训练数据集,对KDE进行训练。
4.2)以t1st-t3时间段的异常KPI序列作为KDE的测试数据集,并将测试数据集输入到KDE中。
4.3)利用KDE对测试数据集拟合,计算出t1st-t3时间段的每个时间点对应的异常分数值,并取对数,得到每个时间点对应的异常分数。将异常分数写入异常分数矩阵A中。
5)判断异常分数矩阵A中是否存在元素a>α,若是,则记元素a对应的KPI值v为突变指标,KPI值v对应小区记为指标突变小区。α为预设的阈值。
值得说明的是,为了发挥经验门限检测方法和密度函数检测方法的优势并弥补两者的不足,本发明使用每个小区每个KPI的历史数据计算出该指标的门限值,基于这个门限值找到异常小区和异常KPI。之后,对于这些异常小区的异常KPI,本发明使用基于高斯核心的KDE检测方法来判断该指标是否为突变指标,接着判断该小区是否为指标突变小区。这样可以尽可能的避免两种方法的弊端的同时,以尽可能高的准确率定位指标突变小区,帮助指导定点优化。
本发明的技术效果是毋庸置疑的。本发明解决了现有的经验门限检测方法门限无法实时更新和不完善的缺陷和密度函数存在非正常训练数据干扰和训练数据不足的问题。本发明使用历史数据来动态生成门限以检测得到异常小区和异常的KPI,并使用异常KPI的正常的历史数据作为训练数据来拟合高斯核心的KDE,输出每个数据点的指标突变程度,以得到指标突变小区。该检测方法可精确定位指标突变小区,灵活性高,鲁棒性强,性能稳定。
附图说明
图1为网络中检测指标突变小区的过程示意图。
图2为指标突变检测系统的具体数据流程图。
具体实施方式
下面结合实施例对本发明作进一步说明,但不应该理解为本发明上述主题范围仅限于下述实施例。在不脱离本发明上述技术思想的情况下,根据本领域普通技术知识和惯用手段,做出各种替换和变更,均应包括在本发明的保护范围内。
实施例1:
参见图1至图2,移动通信网络中基于异常检测和核密度估计KDE的指标突变小区检测方法,主要包括以下步骤:
1)获取移动通信网络所有小区KPI(关键性能指标,Key Performance Indicator)数据序列,并对每个小区的关键性能指标KPI序列进行数据清洗,主要步骤如下:
1.1)基于设定的粒度g,筛选掉粒度不合规的小区。
1.2)基于设定的数据长度S=T/g,筛选掉时间长度不合规的小区。T为指标突变小区检测周期。
1.3)以0填充剩余小区缺失值。
1.4)基于业务告警信息,筛选掉业务告警小区。
2)基于清洗后的KPI序列,建立历史数据序列和测试数据序列。
进一步,历史数据序列为t0-t1时间段的KPI序列。测试数据序列为t2-t3时间段的KPI序列。t0为KPI序列起始时间。t1∈(t0,t3)为历史数据序列终结时间。t2∈(t0,t3)为测试数据序列起始时间。t3为KPI序列终结时间。
3)基于历史数据序列和测试数据序列,筛选掉全部KPI序列均正常的小区,保留异常小区,并记录每个小区的所有异常KPI的第一个发生异常的时间点t1st
确定异常小区的主要步骤如下:
3.1)利用季节差分法对历史数据序列和测试数据序列进行差分。将差分后的历史数据序列和测试数据序列输入到加法模型中,得到历史趋势数据集和测试趋势数据集。
加法模型Y[t]=T[t]+S[t]+e[t]。其中,T[t]、S[t]、e[t]为KPI差分后得到的三个分解因子。Y[t]表示加法模型输出的历史趋势数据或测试趋势数据。
3.2)基于历史趋势数据集,分别计算不同小区中不同KPI序列的均值μ和标准差σ,并确定每个KPI的正常区间I。
KPI的正常区间I=(μ-ω×σ,μ+ω×σ)。ω为常数。
3.3)判断测试趋势数据集中每个元素是否位于对应KPI序列的正常区间I内,若否,则将对应KPI标记为异常指标,并将所述异常指标所在小区标记为异常小区。
4)利用基于高斯核的KDE(核密度估计,Kernel Density Estimation)计算异常小区差分前的异常KPI序列的异常分数,并建立异常分数矩阵A。
建立异常分数矩阵A的主要步骤如下:
4.1)以t0-t1st时间段的异常KPI序列作为KDE的训练数据集,对KDE进行训练。
4.2)以t1st-t3时间段的异常KPI序列作为KDE的测试数据集,并将测试数据集输入到KDE中。
4.3)利用KDE对测试数据集拟合,计算出t1st-t3时间段的每个时间点对应的异常分数值,并取对数,得到每个时间点对应的异常分数。将异常分数写入异常分数矩阵A中。
5)判断异常分数矩阵A中是否存在元素a>α,若是,则记元素a对应的v关键性能指标KPI为突变指标,v对应小区记为指标突变小区。α为预设的阈值。
实施例2:
移动通信网络中基于异常检测和核密度估计KDE的指标突变小区检测方法,主要包括以下步骤:
1)获取移动通信网络所有小区KPI序列,并对每个小区的KPI序列进行数据清洗。
2)基于清洗后的KPI序列,建立历史数据序列和测试数据序列。
3)基于历史数据序列和测试数据序列,筛选掉全部KPI序列均正常的小区,保留异常小区,并记录每个小区的所有异常KPI的第一个发生异常的时间点t1st
4)利用基于高斯核的KDE(Kernel Density Estimation)计算异常小区差分前的异常KPI序列的异常分数,并建立异常分数矩阵A。
5)判断异常分数矩阵A中是否存在元素a>α,若是,则记元素a对应的KPI为突变指标,对应小区记为指标突变小区。α为预设的阈值。
实施例3:
移动通信网络中基于异常检测和核密度估计(KDE)指标突变小区检测方法,主要步骤见实施例2,其中,确定异常小区的主要步骤如下:
1)利用季节差分法对历史数据序列和测试数据序列进行差分。将差分后的历史数据序列和测试数据序列输入到加法模型中,得到历史趋势数据集和测试趋势数据集。
加法模型Y[t]=T[t]+S[t]+e[t]。其中,T[t]、S[t]、e[t]为KPI差分后得到的三个分解因子。Y[t]表示加法模型输出的历史趋势数据或测试趋势数据。
2)基于历史趋势数据集,分别计算不同小区中不同KPI序列的均值μ和标准差σ,并确定每个KPI的正常区间I。
KPI的正常区间I=(μ-ω×σ,μ+ω×σ)。ω为常数。
3)判断测试趋势数据集中每个元素是否位于对应KPI序列的正常区间I内,若否,则将对应KPI标记为异常指标,并将所述异常指标所在小区标记为异常小区。
实施例4:
移动通信网络中基于异常检测和核密度估计KDE的指标突变小区检测方法,主要步骤见实施例2,其中,建立异常分数矩阵A的主要步骤如下:
1)以t0-t1st时间段的异常KPI序列作为KDE的训练数据集,对KDE进行训练。
2)以t1st-t3时间段的异常KPI序列作为KDE的测试数据集,并将测试数据集输入到KDE中。
3)利用KDE对测试数据集拟合,计算出t1st-t3时间段的每个时间点对应的异常分数值,并取对数,得到每个时间点对应的异常分数。将异常分数写入异常分数矩阵A中。

Claims (6)

1.移动通信网络中基于异常检测和核密度估计KDE的指标突变小区检测方法,其特征在于,主要包括以下步骤:
1)获取移动通信网络所有小区KPI序列,并对每个小区的KPI序列进行数据清洗;
2)基于清洗后的KPI序列,建立历史数据序列和测试数据序列;
3)基于历史数据序列和测试数据序列,筛选掉全部KPI序列均正常的小区,保留异常小区,并记录每个小区的所有异常KPI的第一个发生异常的时间点t1st
4)利用基于高斯核的KDE计算异常小区差分前的异常KPI序列的异常分数,并建立所述异常分数矩阵A;
建立异常分数矩阵A的主要步骤如下:
4.1)以t0-t1st时间段的异常KPI序列作为KDE的训练数据集,对KDE进行训练;t0为KPI序列起始时间;
4.2)以t1st-t3时间段的异常KPI序列作为KDE的测试数据集,并将测试数据集输入到KDE中;t3为KPI序列终结时间;
4.3)利用KDE对测试数据集拟合,计算出t1st-t3时间段的每个时间点对应的异常分数值,并取对数,得到每个时间点对应的异常分数;将异常分数写入异常分数矩阵A中;
5)判断异常分数矩阵A中是否存在元素a>α,若是,则记元素a对应的KPI值v为突变指标,KPI值v对应小区记为指标突变小区;α为预设的阈值。
2.根据权利要求1所述的移动通信网络中基于异常检测和核密度估计KDE的指标突变小区检测方法,其特征在于,对每个小区的KPI序列进行数据清洗的主要步骤如下:
1)基于设定的粒度g,筛选掉粒度不合规的小区;
2)基于设定的数据长度S=T/g,筛选掉时间长度不合规的小区;T为指标突变小区检测周期;
3)以0填充剩余小区缺失值;
4)基于业务告警信息,筛选掉业务告警小区。
3.根据权利要求1或2所述的移动通信网络中基于异常检测和核密度估计KDE的指标突变小区检测方法,其特征在于,历史数据序列为t0-t1时间段的KPI序列;测试数据序列为t2-t3时间段的KPI序列;t0为KPI序列起始时间;t1∈(t0,t3)为历史数据序列终结时间;t2∈(t0,t3)为测试数据序列起始时间;t3为KPI序列终结时间。
4.根据权利要求1或2所述的移动通信网络中基于异常检测和核密度估计KDE的指标突变小区检测方法,其特征在于,确定异常小区的主要步骤如下:
1)利用季节差分法对历史数据序列和测试数据序列进行差分;将差分后的历史数据序列和测试数据序列输入到加法模型中,得到历史趋势数据集和测试趋势数据集;
2)基于历史趋势数据集,分别计算不同小区中不同KPI序列的均值μ和标准差σ,并确定每个KPI的正常区间I;
3)判断测试趋势数据集中每个元素是否位于对应KPI序列的正常区间I内,若否,则将对应KPI标记为异常指标,并将所述异常指标所在小区标记为异常小区。
5.根据权利要求4所述的移动通信网络中基于异常检测和核密度估计KDE的指标突变小区检测方法,其特征在于,加法模型Y[t]=T[t]+S[t]+e[t];其中,T[t]、S[t]、e[t]为KPI差分后得到的三个分解因子;Y[t]表示加法模型输出的历史趋势数据或测试趋势数据。
6.根据权利要求4所述的移动通信网络中基于异常检测和核密度估计KDE的指标突变小区检测方法,其特征在于,KPI的正常区间I=(μ-ω×σ,μ+ω×σ);ω为常数。
CN201911402828.9A 2019-12-31 2019-12-31 移动通信网络中基于异常检测和核密度估计kde的指标突变小区检测方法 Active CN111065106B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911402828.9A CN111065106B (zh) 2019-12-31 2019-12-31 移动通信网络中基于异常检测和核密度估计kde的指标突变小区检测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911402828.9A CN111065106B (zh) 2019-12-31 2019-12-31 移动通信网络中基于异常检测和核密度估计kde的指标突变小区检测方法

Publications (2)

Publication Number Publication Date
CN111065106A CN111065106A (zh) 2020-04-24
CN111065106B true CN111065106B (zh) 2022-04-26

Family

ID=70305097

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911402828.9A Active CN111065106B (zh) 2019-12-31 2019-12-31 移动通信网络中基于异常检测和核密度估计kde的指标突变小区检测方法

Country Status (1)

Country Link
CN (1) CN111065106B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111813645B (zh) * 2020-07-16 2022-04-12 莫毓昌 一种基于决策图的在线指标异常检测方法
CN111897695B (zh) * 2020-07-31 2022-06-17 平安科技(深圳)有限公司 获取kpi异常数据样本的方法、装置和计算机设备

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103945442A (zh) * 2014-05-07 2014-07-23 东南大学 移动通信系统中基于线性预测原理的系统异常检测方法
CN105230063A (zh) * 2013-05-02 2016-01-06 华为技术有限公司 网络优化的方法、网络优化的装置和网络优化的设备
CN107092582A (zh) * 2017-03-31 2017-08-25 江苏方天电力技术有限公司 一种基于残差后验的异常值在线检测及置信度评估方法
CN110278121A (zh) * 2018-03-15 2019-09-24 中兴通讯股份有限公司 一种检测网络性能异常的方法、装置、设备及存储介质

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20160088502A1 (en) * 2013-05-14 2016-03-24 Nokia Solutions And Networks Oy Method and network device for cell anomaly detection
US9544812B2 (en) * 2015-01-30 2017-01-10 Alcatel Lucent System and method for mitigating network congestion using fast congestion detection in a wireless radio access network (RAN)

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105230063A (zh) * 2013-05-02 2016-01-06 华为技术有限公司 网络优化的方法、网络优化的装置和网络优化的设备
CN103945442A (zh) * 2014-05-07 2014-07-23 东南大学 移动通信系统中基于线性预测原理的系统异常检测方法
CN107092582A (zh) * 2017-03-31 2017-08-25 江苏方天电力技术有限公司 一种基于残差后验的异常值在线检测及置信度评估方法
CN110278121A (zh) * 2018-03-15 2019-09-24 中兴通讯股份有限公司 一种检测网络性能异常的方法、装置、设备及存储介质

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
"Anomalous Cell Detection with Kernel Density-Based Local Outlier Factor";Miao Dandan等;《中国通信》;20150915;第12卷(第09期);73-84 *
"Anomaly Detection Algorithms for the Sleeping Cell Detection in LTE Networks";Sergey Chernov;《2015 IEEE 81st Vehicular Technology Conference (VTC Spring)》;20150702;全文 *
"基于大数据的时间序列预测研究与应用";程艳云等;《计算机技术与发展》;20160322;第26卷(第06期);181-184 *

Also Published As

Publication number Publication date
CN111065106A (zh) 2020-04-24

Similar Documents

Publication Publication Date Title
CN111065106B (zh) 移动通信网络中基于异常检测和核密度估计kde的指标突变小区检测方法
CN108304348B (zh) 一种基于二元维纳过程的轴承剩余寿命预测方法
KR101677008B1 (ko) Tsd 기반 데이터 예측 방법
CN111148142B (zh) 基于异常检测和集成学习的休眠小区检测方法
CN107133635A (zh) 一种多源异构飞行事故航迹数据融合方法
CN110458791B (zh) 质量缺陷检测方法和检测设备
CN110738346A (zh) 一种基于威布尔分布的批次电能表可靠度预测方法
KR102158100B1 (ko) 이상 감지를 이용한 모니터링 자동화 방법 및 장치
CN115705413A (zh) 异常日志的确定方法及装置
CN107330587B (zh) 一种基于ais的近岸航道水深利用率计算方法
CN101789091A (zh) 视频清晰度自动识别系统及其方法
CN117312769A (zh) 一种基于BiLSTM的物联网时序数据异常检测方法
CN111090585A (zh) 一种基于众测过程的众测任务关闭时间自动预测方法
KR101945131B1 (ko) 비정규분포 공정에서의 극소불량률 관리 방법 및 장치
CN116545111A (zh) 一种基于变电站内外网的信息交互系统
CN115546108A (zh) 基于边云协同和ar的汽车轮胎外观质量智能检测方法
CN111245531B (zh) 通过分析设备模块功率进行参数自适应调节的方法
CN108260148B (zh) 一种故障检测方法及装置
CN115442271B (zh) 一种网络性能指标时序数据异常检测方法
CN110569277A (zh) 一种配置数据信息自动识别与归类方法及系统
CN117406161B (zh) 一种计量装置量值偏离度预警方法、系统、设备及介质
CN115936266B (zh) 轨道交通设备的可靠度预测方法、系统、设备和介质
CN113822258B (zh) 一种线缆长度计量方法、装置及存储介质
CN117455911B (zh) 一种网线水晶头生产测试方法及系统
CN117870779B (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