CN108992061B - 运动训练负荷监测方法 - Google Patents

运动训练负荷监测方法 Download PDF

Info

Publication number
CN108992061B
CN108992061B CN201810616253.XA CN201810616253A CN108992061B CN 108992061 B CN108992061 B CN 108992061B CN 201810616253 A CN201810616253 A CN 201810616253A CN 108992061 B CN108992061 B CN 108992061B
Authority
CN
China
Prior art keywords
electrocardiogram
exercise
data
ring
training load
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
CN201810616253.XA
Other languages
English (en)
Other versions
CN108992061A (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.)
Hangzhou Xinchang Technology Co ltd
Original Assignee
Hangzhou Xinchang 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 Hangzhou Xinchang Technology Co ltd filed Critical Hangzhou Xinchang Technology Co ltd
Priority to CN201810616253.XA priority Critical patent/CN108992061B/zh
Publication of CN108992061A publication Critical patent/CN108992061A/zh
Application granted granted Critical
Publication of CN108992061B publication Critical patent/CN108992061B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/24Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
    • A61B5/316Modalities, i.e. specific diagnostic methods
    • A61B5/318Heart-related electrical modalities, e.g. electrocardiography [ECG]
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/72Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7235Details of waveform analysis
    • A61B5/7246Details of waveform analysis using correlation, e.g. template matching or determination of similarity

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Molecular Biology (AREA)
  • Animal Behavior & Ethology (AREA)
  • Pathology (AREA)
  • Physics & Mathematics (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Medical Informatics (AREA)
  • Veterinary Medicine (AREA)
  • Surgery (AREA)
  • Biophysics (AREA)
  • General Health & Medical Sciences (AREA)
  • Public Health (AREA)
  • Cardiology (AREA)
  • Artificial Intelligence (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Physiology (AREA)
  • Psychiatry (AREA)
  • Signal Processing (AREA)
  • Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)

Abstract

本发明公开了一种运动训练负荷监测方法,包括以下步骤:1)对监测对象在某种程度下的运动负荷训练后的特定时间,采集该监测对象多导联心电图数据;2)提取所述心电图数据中心电信号的ST‑T段数据,对心电图数据中的ST‑T段或T波数据进行心电动力学建模,得到心电动力学图;3)计算所述心电动力学图的时间异质度和空间异质度,得到所述检测对象的运动训练负荷监测指标。本发明将心电信号看作一类动态模式,使用动态模式建模方法对心电信号进行分析,提取其内部深层次的动力学信息,突破了已有提取心电信号时频特征、统计特征等范畴,并将提取到的动力学信息用于运动训练负荷的评估。

Description

运动训练负荷监测方法
技术领域
本发明涉及运动训练技术领域,具体地说是一种运动训练负荷监测方法。
背景技术
加强运动员训练和全民健身的科学化水平,一个关键问题是如何根据个体的生理特点和运动应激反应选择适宜的运动训练负荷。理想的训练计划应当遵循运动训练中循序渐进的超负荷原则,即逐渐增加运动负荷,逐步提高运动能力,待身体适应该运动负荷后,再继续增加运动负荷,使运动能力不断提高。只有准确、及时地发现运动员在训练负荷下的身体机能变化状况,才能合理地安排和调整运动训练计划,避免过度疲劳,减少运动损伤,提高运动训练的有效性和安全性。
虽然肌酸激酶、肌钙蛋白等指标已经被证明为监测运动负荷的有效指标,但是,这些传统的运动负荷监测方法的共同缺点是:1)需要采血,有创测试,不受运动员和普通运动健身者欢迎,也不可能长期频繁使用;2)时间滞后,一般需要几个小时到几天出结果,造成的运动损伤已经很难挽回;3)测试程序比较复杂,测试价格比较昂贵,许多专业运动队都很难频繁使用,更不适用于普通运动健身者的科学健身指导。
心电图运动负荷试验又称心电运动试验,是通过一定量的运动增加心脏负荷,观察心电图变化,对已知或怀疑患有心血管疾病,尤其是冠状动脉粥样硬化性心脏病进行临床评估的方法。在运动心电图中,可以观察的心电参数包括心率、QRS波、ST段、T波和心律失常。心电运动试验在临床上的应用包括:诊断冠心病;评定已确诊冠心病的严重程度;预测心血管事件和心源性死亡;评定运动能力和耐量;评定运动相关的症状;评估心率变时性功能、心律失常和对植入性器械治疗的反应;评估对医学介入治疗的反应。
综上,现有技术无法快速、便捷地实现运动训练负荷监测,而且有创测试存在着固有缺陷。现有技术心电运动试验还不能实现运动训练负荷监测,鉴于心电运动试验的无创特性和便捷性,若能将其用于运动训练负荷监测,可能会有较大优势。
发明内容
有鉴于此,本发明针对上述与运动训练负荷监测相关的现有技术存在的有创、测试程序复杂、测试价格昂贵等技术问题,提供一种无创、快速、便捷的运动训练负荷监测方法,并降低了监测成本。
本发明的技术解决方案是,提供一种运动训练负荷监测方法,包括以下步骤:
1)对监测对象在某种程度下的运动负荷训练后的特定时间,采集该监测对象多导联心电图数据;
2)提取所述心电图数据中心电信号的ST-T段数据,对心电图数据中的ST-T段或T波数据进行心电动力学建模,得到心电动力学图;
3)计算所述心电动力学图的时间异质度和空间异质度,得到所述检测对象的运动训练负荷监测指标。
作为可选,所述对心电图数据中的ST-T段或T波数据进行动力学建模由以下方法实现:将提取的ST-T段或T波数据转换为三维数据,得到ST-T环或T环;
采用神经网络辨识器,对ST-T环或T环的内在系统动态利用确定学习算法进行局部神经网络逼近,获得关于心电图ST-T段或T环数据内在的动力学特征;
将利用神经网络得到的关于心电图ST-T段或T环数据内在的动力学特征沿着ST-T环或T环轨迹进行三维可视化显示,得到心电动力学图;所述心电动力学图包含心电信号ST-T段的状态信息和沿着ST-T段或T环状态轨迹的动力学特征。
作为可选,通过计算ST-T环或T环的空间异质度SI和时间异质度TI,得到运动训练负荷监测指标E,其中
Figure BDA0001696917700000021
E=a×TI-b×SI+c;其中a、b和c为可变系数,根据运动训练负荷监测指标E与血检结果的一致性确定,运动训练负荷监测指标E>0时,表征运动训练负荷过大,E<0时,表征运动训练负荷正常。
作为可选,运动训练负荷监测指标E=-0.0018TI-SI+0.4,E>0时,表征运动训练负荷过大,E<0时,表征运动训练负荷正常。
作为可选,通过采集20-60秒12导联心电图,提取其中18-20个心动周期的ST-T段或T波数据进行动力学建模分析,从而得到运动训练负荷监测指标。
作为可选,通过采集20秒12导联心电图,提取其中20个心动周期的ST-T段或T波数据进行动力学建模分析,从而得到运动训练负荷监测指标。
作为可选,分别在无运动负荷时、低运动负荷训练、中运动负荷训练,和高运动负荷训练三种运动训练后即刻和24小时后,采集样本人群的静脉血,分别测量20个心动周期的12导联心电图相关数据,并填写主观感觉疲劳问卷;通过静脉血标本,采集受试对象的肌酸激酶-MB型同工酶(CK-MB)、高敏肌钙蛋白I(hsTnI),采用RPE(ratingofperceivedexertion)量表测量样本人群的主观疲劳程度,以采血的检测结果作为对心电图测试结果的验证和校正的标准,并对心电图测试结果进行准确度评价。
采用以上运动训练负荷监测方法,本发明与现有技术相比,具有以下优点:(1)本发明的心电动力学图采用动态模式建模与识别的新方法,对心电信号进行动力学建模,提取逐拍心电信号ST段或T波中的动力学特征,突破已有的仅提取心电信号时频特征、统计特征的范畴,对心电信号中微小的T波电交替进行更为敏感的检测,以达到运动训练负荷监测的目的,通过建立心电动力学图进一步提出通过心电信号ST-T段或T波时间异质度和空间异质度这两个指标来检测微伏级心电改变,从而达到评估运动负荷水平的目的;(2)本发明将心电信号看作一类动态模式,使用动态模式建模方法对心电信号进行分析,提取其内部深层次的动力学信息,突破了已有提取心电信号时频特征、统计特征等范畴,并将提取到的动力学信息用于运动训练负荷的评估。
附图说明
图1为健康成年男性的标准12导联心电图;
图2为健康成年男性的心电动力学图;
图3为运动负荷过大引起心肌缺血的运动员的标准12导联心电图;
图4为运动负荷过大引起心肌缺血的运动员的心电动力学图;
具体实施方式
下面结合附图和具体实施例对本发明作进一步说明,但本发明并不仅仅限于这些实施例。
本发明涵盖任何在本发明的精髓和范围上做的替代、修改、等效方法以及方案。为了使公众对本发明有彻底的了解,在以下本发明优选实施例中详细说明了具体的细节,而对本领域技术人员来说没有这些细节的描述也可以完全理解本发明。
本发明提供一种运动训练负荷监测方法,包括以下步骤:
1)对监测对象在某种程度下的运动负荷训练后的特定时间,采集该监测对象多导联心电图数据;
2)提取所述心电图数据中心电信号的ST-T段数据,对心电图数据中的ST-T段或T波数据进行心电动力学建模,得到心电动力学图;
3)计算所述心电动力学图的时间异质度和空间异质度,得到所述检测对象的运动训练负荷监测指标。
所述对心电图数据中的ST-T段或T波数据进行动力学建模由以下方法实现:
将提取的ST-T段或T波数据转换为三维数据,得到ST-T环或T环;
采用神经网络辨识器,对ST-T环或T环的内在系统动态利用确定学习算法进行局部神经网络逼近,获得关于心电图ST-T段或T环数据内在的动力学特征;
将利用神经网络得到的关于心电图ST-T段或T环数据内在的动力学特征沿着ST-T环或T环轨迹进行三维可视化显示,得到心电动力学图;所述心电动力学图包含心电信号ST-T段的状态信息和沿着ST-T段或T环状态轨迹的动力学特征。
更具体地,所述神经网络辨识器为动态径向基函数(radial basis function,RBF)神经网络辨识器,对ST-T环的内在系统动态利用确定学习算法进行局部准确RBF神经网络逼近,获得关于心电图ST-T段信号内在的动力学特征,这种系统动力学特征,与静态心电信号特征相比,是对心电信号内部更深层次更本质的一种特征描述;
把RBF神经网络的建模结果沿着ST-T环轨迹进行三维可视化显示,得到心电动力学图。心电动力学图既包含了心电信号ST-T段的状态信息,也包含了沿着ST-T段状态轨迹的动力学特征,是一种全息特征,在空间和时间上刻画出心脏复极过程电活动的动力学信息。
通过计算ST-T环或T环的空间异质度SI和时间异质度TI,得到运动训练负荷监测指标E,其中
Figure BDA0001696917700000041
E=a×TI-b×SI+c;其中a、b和c为可变系数,根据运动训练负荷监测指标E与血检结果的一致性确定,运动训练负荷监测指标E>0时,表征运动训练负荷过大,E<0时,表征运动训练负荷正常。以a=-0.0018,b=1,c=0.4为例,E=-0.0018TI-SI+0.4。
如图2所示,健康成年男性的心电动力学图比较规整。而运动训练负荷过量引起心肌缺血时,其心电动力学图明显比较紊乱,如图4所示。
本发明的另一个实施例,对心电信号进行动力学建模包括以下步骤:
1)截取心电信号的T波数据,将12导联T波数据转换为三维数据,得到T环;
2)采用动态径向基函数(radial basis function,RBF)神经网络辨识器,对T环的内在系统动态利用确定学习算法进行局部准确RBF神经网络逼近,获得关于心电图T波信号内在的动力学特征,这种系统动力学特征,与静态心电信号特征相比,是对心电信号内部更深层次更本质的一种特征描述;
3)把RBF神经网络的建模结果沿着T环轨迹进行三维可视化显示,得到心电动力学图。
通过计算T环的空间异质度SI和时间异质度TI,得到运动训练负荷监测指标E,其中
Figure BDA0001696917700000051
E=-0.0018TI-SI+0.4。运动训练负荷监测指标E>0时,运动训练负荷过大,E<0时,运动训练负荷正常。
通过采集20-60秒12导联心电图,提取其中18-20个心动周期的ST-T段或T波数据进行动力学建模分析,从而得到运动训练负荷监测指标。
通过采集20秒12导联心电图,提取其中20个心动周期的ST-T段或T波数据进行动力学建模分析,从而得到运动训练负荷监测指标。
得到带有时间戳的运动训练负荷监测指标后,将其传输至云端个人健康档案,通过手持终端访问云端数据,并设置访问权限,只有本人和经过本人授权的其他人员才能具有访问权限,不可篡改运动训练负荷监测指标数据。手持终端上设有北斗定位系统,在通过手持终端访问数据时,手持终端强制定位,当定位结果发现与被授权人的常住地不符合时,则手持终端关闭数据访问权限,需要本人生物识别特征验证后才能重新获得权限,以防止个人运动信息以及指标数据的泄露。大量的特定人群的运动训练负荷监测指标在云端形成大数据,通过大数据形成数据分析报告,包括年龄、职业、性别、身体质量指数、高血压、糖尿病等信息与运动训练负荷监测指标的关系。
本发明的方法还可利用得到的心电动力学图,即运动负荷心电图,和12导联动态心电图作为一个综合的指标,精准度更高。在做12导联心电图的时候同时测量心电图∑T/R,利用血清心肌肌钙蛋白,血浆心纳素,心电图∑T/R三个指标综合,作为验证和校正的标准,并对心电图测试结果进行准确度评价。这些技术特征均可与上述两个实施例进行结合。
以300名高水平运动员和资深运动健身人群作为研究对象,分别在无运动负荷时、低运动负荷训练、中运动负荷训练,和高运动负荷训练三种运动训练后即刻和24小时后,分别采集静脉血,分别测量20个心动周期的12导联心电图相关数据,并填写主观感觉疲劳问卷。通过静脉血标本,采集受试对象的肌酸激酶-MB型同工酶(CK-MB)、高敏肌钙蛋白I(hsTnI),采用RPE(rating ofperceived exertion)量表测量受试对象的主观疲劳程度。通过对20个心动周期的12导联心电图建立心电动力学图,计算运动训练负荷监测指标E。以CK-MB测试结果为金标准,CDG测试的准确率、灵敏度和特异性均为80%。以hsTnI测试结果为金标准,CDG测试的准确率为80.4%,灵敏度为71.5%,特异度为81.8%。
以上仅就本发明较佳的实施例作了说明,但不能理解为是对权利要求的限制。本发明不仅局限于以上实施例,其具体结构允许有变化。总之,凡在本发明独立权利要求的保护范围内所作的各种变化均在本发明的保护范围内。

Claims (5)

1.一种运动训练负荷监测方法,包括以下步骤:
1)对监测对象在某种程度下的运动负荷训练后的特定时间,采集该监测对象多导联心电图数据;
2)提取所述心电图数据中心电信号的ST-T段数据,对心电图数据中的ST-T段或T波数据进行心电动力学建模,得到心电动力学图;
3)计算所述心电动力学图的时间异质度和空间异质度,得到所述监 测对象的运动训练负荷监测指标;
所述对心电图数据中的ST-T段或T波数据进行动力学建模由以下方法实现:将提取的ST-T段或T波数据转换为三维数据,得到ST-T环或T环;
采用神经网络辨识器,对ST-T环或T环的内在系统动态利用确定学习算法进行局部神经网络逼近,获得关于心电图ST-T段或T环数据内在的动力学特征;
将利用神经网络得到的关于心电图ST-T段或T环数据内在的动力学特征沿着ST-T环或T环轨迹进行三维可视化显示,得到心电动力学图;所述心电动力学图包含心电信号ST-T段的状态信息和沿着ST-T段或T环状态轨迹的动力学特征;
通过计算ST-T环或T环的空间异质度SI和时间异质度TI,得到运动训练负荷监测指标E,其中
Figure FDA0002639257590000011
E=a×TI-b×SI+c;其中a、b和c为可变系数,根据运动训练负荷监测指标E与血检结果的一致性确定,运动训练负荷监测指标E>0时,表征运动训练负荷过大,E<0时,表征运动训练负荷正常。
2.如权利要求1所述的运动训练负荷监测方法,其特征在于:运动训练负荷监测指标E=-0.0018TI-SI+0.4,E>0时,表征运动训练负荷过大,E<0时,表征运动训练负荷正常。
3.如权利要求2所述的运动训练负荷监测方法,其特征在于:通过采集20-60秒12导联心电图,提取其中18-20个心动周期的ST-T段或T波数据进行动力学建模分析,从而得到运动训练负荷监测指标。
4.如权利要求3所述的运动训练负荷监测方法,其特征在于:通过采集20秒12导联心电图,提取其中20个心动周期的ST-T段或T波数据进行动力学建模分析,从而得到运动训练负荷监测指标。
5.如权利要求2所述的运动训练负荷监测方法,其特征在于:分别在无运动负荷时、低运动负荷训练、中运动负荷训练,和高运动负荷训练三种运动训练后即刻和24小时后,采集样本人群的静脉血,分别测量20个心动周期的12导联心电图相关数据,并填写主观感觉疲劳问卷;通过静脉血标本,采集受试对象的肌酸激酶-MB型同工酶(CK-MB)、高敏肌钙蛋白I(hsTnI),采用RPE(rating ofperceived exertion)量表测量样本人群的主观疲劳程度,以采血的检测结果作为对心电图测试结果的验证和校正的标准,并对心电图测试结果进行准确度评价。
CN201810616253.XA 2018-06-14 2018-06-14 运动训练负荷监测方法 Active CN108992061B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810616253.XA CN108992061B (zh) 2018-06-14 2018-06-14 运动训练负荷监测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810616253.XA CN108992061B (zh) 2018-06-14 2018-06-14 运动训练负荷监测方法

Publications (2)

Publication Number Publication Date
CN108992061A CN108992061A (zh) 2018-12-14
CN108992061B true CN108992061B (zh) 2020-09-25

Family

ID=64601215

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810616253.XA Active CN108992061B (zh) 2018-06-14 2018-06-14 运动训练负荷监测方法

Country Status (1)

Country Link
CN (1) CN108992061B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN118266874A (zh) * 2024-06-04 2024-07-02 山东大学 基于确定学习的单通道气流睡眠呼吸暂停实时检测系统

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1291904A (zh) * 1998-02-26 2001-04-18 欧姆龙株式会社 运动器械、体力评价方法以及搏动计
CN102397067A (zh) * 2011-12-07 2012-04-04 秦皇岛市康泰医学系统有限公司 一种运动负荷心电波形叠加分析方法及其系统
CN106921729A (zh) * 2017-01-16 2017-07-04 中州大学 一种农业墒情监测系统

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20080057590A1 (en) * 2006-06-07 2008-03-06 Mickey Urdea Markers associated with arteriovascular events and methods of use thereof
US10022060B2 (en) * 2012-09-21 2018-07-17 Beth Israel Deaconess Medical Center, Inc. High throughput arrhythmia risk assessment using multilead residua signals

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1291904A (zh) * 1998-02-26 2001-04-18 欧姆龙株式会社 运动器械、体力评价方法以及搏动计
CN102397067A (zh) * 2011-12-07 2012-04-04 秦皇岛市康泰医学系统有限公司 一种运动负荷心电波形叠加分析方法及其系统
CN106921729A (zh) * 2017-01-16 2017-07-04 中州大学 一种农业墒情监测系统

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
A new method for early detection of myocardial ischemia: cardiodynamicsgram(CDG);Cong Wang, et al;《Science China:Information Sciences》;20160131;全文 *
基于心电动力学图的远程心肌缺血检测系统的设计与实现;雷鹏;《中国优秀硕士论文数据库》;20171231;全文 *

Also Published As

Publication number Publication date
CN108992061A (zh) 2018-12-14

Similar Documents

Publication Publication Date Title
Ge et al. Cardiac arrhythmia classification using autoregressive modeling
US6390986B1 (en) Classification of heart rate variability patterns in diabetics using cepstral analysis
Brage et al. Reliability and validity of the combined heart rate and movement sensor Actiheart
Cohn et al. Use of treadmill score to quantify ischemic response and predict extent of coronary disease.
KR102008196B1 (ko) 심전도 데이터를 이용한 혈중 칼륨농도 예측모델 생성장치 및 그 방법
WO2007149985A2 (en) Assessing dementia and dementia-type disorders
US20020107455A1 (en) Method and apparatus for estimating degree of neuronal impairment in brain cortex
WO2001097066A1 (en) Methods and apparatus for providing an indicator of autonomic nervous system function
Tiwari et al. A comparative study of stress and anxiety estimation in ecological settings using a smart-shirt and a smart-bracelet
US20020156392A1 (en) Method and apparatus for inspecting biological rhythms
Cataldo et al. Resting sympatho-vagal balance is related to 10 km running performance in master endurance athletes
CN108992061B (zh) 运动训练负荷监测方法
CN108786074B (zh) 运动训练负荷监测设备
Ausland et al. Heart rhythm assessment in elite endurance athletes: a better method?
Garatachea et al. Estimation of energy expenditure in healthy adults from the YMCA submaximal cycle ergometer test
Guidi et al. Heart rate variability analysis during muscle fatigue due to prolonged isometric contraction
Hounaida et al. ST-based Deep Learning analysis of COVID-19 patients
Chaikovsky et al. ECG Universal score system: new instrument for electrocardiogram analysis
Manukova et al. An Approach to Evaluation of Clinically Healthy People by Preventive Cardio Control
Dandu et al. Determining physiological significance of inertial gait features in multiple sclerosis
Hikmah et al. Sleep Quality Assessment from Robust Heart and Muscle Fatigue Estimation Using Supervised Machine Learning.
Caesaria et al. Effect of Muscle Fatigue on EMG Signal and Maximum Heart Rate for Pre and Post Physical Activity
Moghadam et al. Quantification of Parkinson tremor intensity based on EMG signal analysis using fast orthogonal search algorithm
Moody et al. Distinguishing ischemic from non-ischemic ST changes: the Physionet/Computers in Cardiology Challenge 2003
Moatamed et al. Sport analytics platform for athletic readiness assessment

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