CN108345039B - 一种消除地面核磁共振数据中邻频谐波干扰的方法 - Google Patents

一种消除地面核磁共振数据中邻频谐波干扰的方法 Download PDF

Info

Publication number
CN108345039B
CN108345039B CN201810029427.2A CN201810029427A CN108345039B CN 108345039 B CN108345039 B CN 108345039B CN 201810029427 A CN201810029427 A CN 201810029427A CN 108345039 B CN108345039 B CN 108345039B
Authority
CN
China
Prior art keywords
harmonic
frequency
adjacent
data
magnetic resonance
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
CN201810029427.2A
Other languages
English (en)
Other versions
CN108345039A (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.)
Jilin University
Original Assignee
Jilin 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 Jilin University filed Critical Jilin University
Priority to CN201810029427.2A priority Critical patent/CN108345039B/zh
Publication of CN108345039A publication Critical patent/CN108345039A/zh
Application granted granted Critical
Publication of CN108345039B publication Critical patent/CN108345039B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V3/00Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
    • G01V3/14Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation operating with electron or nuclear magnetic resonance
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00Arrangements or instruments for measuring magnetic variables
    • G01R33/20Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/46NMR spectroscopy
    • G01R33/4625Processing of acquired signals, e.g. elimination of phase errors, baseline fitting, chemometric analysis

Landscapes

  • Physics & Mathematics (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Remote Sensing (AREA)
  • Geology (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Environmental & Geological Engineering (AREA)
  • Geophysics (AREA)
  • Signal Processing (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Condensed Matter Physics & Semiconductors (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

本发明属于地面核磁共振数据处理领域,为一种消除地面核磁共振数据中邻频谐波的方法。该方法包括:录入地面核磁共振测量数据,测量实验地点地磁场,计算拉莫尔频率;搜索工频谐波基频,并对工频谐波进行建模,计算得到所有谐波系数;根据拉莫尔频率和每个谐波频率的频率差判断是否存在邻频谐波,若否,利用谐波基频和谐波系数计算得到谐波模型,在测量数据中减去谐波模型,若是,则去除领域谐波外的工频谐波后基于拟合方法去除领频谐波,本发明在不损失核磁共振信号的前提下,去除工频谐波中的邻频谐波对后续数据处理和反演的影响。尤其当实验地点的拉莫尔频率与工频谐波的高阶谐波相近时,可准确提取地面核磁共振信号特征,提高反演结果准确度。

Description

一种消除地面核磁共振数据中邻频谐波干扰的方法
技术领域
本发明属于地面核磁共振数据处理领域,具体地来讲为一种消除地面核磁共振数据中邻频谐波的方法。
背景技术
地面核磁共振是一种直接探测地下水赋存属性的地球物理方法,具有测量速度快,分辨率高和定量解释等优点,被广泛应用于水资源勘探、水文环境监测以及地下工程水害预警等领域。但是地面核磁共振方法获得的信号十分微弱,通常强度只有纳伏量级(nV,10-9V),易受到各种环境噪声的干扰,导致采集数据信噪比低,严重影响后续反演结果的准确性,制约了该方法在城镇、村庄以及地下工程等复杂环境中的应用。
地面核磁共振数据中主要包含三类噪声:尖峰噪声、工频及其谐波噪声和随机噪声,其中高压输送电网和电气设备等产生的工频谐波噪声最为强烈。对于工频谐波噪声,当存在其中一个谐波的频率与地面核磁共振信号的频率(称为拉莫尔频率)相近或相等时(小于等于3Hz),我们称之为邻频谐波。
专利CN103913778A公开了一种多个近端参考线圈的核磁共振信号实时噪声抵消装置。该方法首先消除参考线圈内的地面核磁共振信号,并利用参考通道与主通道的频率相关性,实现工频谐波的消除。但是当存在邻频谐波时,相关系数无法正确计算,且铺设参考线圈增加了野外实验的复杂度。
专利CN104898172A公开了一种互相关的核磁共振全波信号噪声滤除方法。该方法通过构造一个与地面核磁共振信号频率相同的正弦信号作为参考信号,并基于随机噪声和尖峰噪声与参考信号相关性较弱,而地面核磁共振信号与参考线圈相关性较强的原则,实现地面核磁共振数据的噪声滤除。但是当存在邻频谐波时,该方法无法有效区分谐波噪声和地面核磁共振信号。
专利CN106772646A公开了一种地面核磁共振信号提取方法。对于工频谐波噪声,提出了利用谐波建模消噪方法,该方法无需额外的参考通道,仅通过采集数据建立工频谐波模型。存在邻频谐波噪声时,该方法无法同样无法区分谐波噪声和地面核磁共振信号,导致建模结果产生错误。
发明内容
本发明所要解决的技术问题在于提供一种消除地面核磁共振数据中工频及邻频谐波干扰的方法,在不损失核磁共振信号的前提下,去除工频谐波中的邻频谐波对后续数据处理和反演的影响。尤其当实验地点的拉莫尔频率与工频谐波的高阶谐波相近时,可准确提取地面核磁共振信号特征,提高反演结果准确度。
本发明是这样实现的,
一种消除地面核磁共振数据中邻频谐波干扰的方法,该方法包括:
步骤1,录入地面核磁共振测量数据V=[V1,V2,L,VK],K为采样点数,测量实验地点地磁场B0,根据式fL=0.04258×B0计算拉莫尔频率fL
步骤2,搜索工频谐波基频f0,并对工频谐波进行建模,计算得到所有谐波系数X;
步骤3,根据拉莫尔频率fL和每个谐波频率fm=m·f0的频率差dfm=fL-fm判断是否存在邻频谐波,其中m=1,2,L,M,M是谐波个数。若所有的频率差|dfm|≥3Hz,说明不存在邻频谐波,利用步骤2得到的谐波基频fo和谐波系数X计算得到谐波模型Vharm,在测量数据V中减去谐波模型Vharm,实现工频及其谐波干扰的消除。
进一步地,步骤3中,若其中一个频率差|dfm|<3Hz,则存在邻频谐波,进行步骤4;
步骤4,将谐波系数X中邻频谐波对应的系数置零,进行谐波建模,得到除邻频谐波外的谐波模型V′harm,在测量数据V中减去谐波模型V′harm,首先完成除邻频谐波外所有工频谐波的消除,得到数据V′,基于拟合的方法,得到邻频谐波Vm的参数,在数据V′中减去邻频谐波Vm,得到去除谐波噪声后的地面核磁共振信号。
进一步地,步骤4中基于拟合的方法包括对数据V′以邻频谐波的频率fm进行希尔伯特变换,并经过低通滤波器后得到数据V′包络的实部VRe和虚部VIm
进一步地,对数据V′包络的实部VRe和虚部VIm进行非线性拟合获得6个未知参数x1~x6,其中四个参数对应于地面核磁共振信号参数,两个参数为邻频谐波的参数,获得邻频谐波的幅度和相位,计算得到邻频谐波模型Vm,在数据V′中减去邻频谐波模型Vm完成邻频谐波的消除。
经过低通滤波器后得到V′包络的实部VRe和虚部VIm分别为:
其中e0和θ分别为地面核磁共振信号的初始幅度,弛豫时间和相位,dfm=fL-fm是拉莫尔频率fL与邻波谐波频率fm的频率差。Am分别为邻频谐波的幅度和相位。εr和εi分别是随机噪声的实部和虚部分量。
进一步地,根据表达式(1)和式(2)对数据V′包络的实部VRe和虚部VIm进行非线性拟合获得6个未知参数x1~x6
其中,分别为经拟合得到的实部和虚部,x1~x4对应于地面核磁共振信号参数,x5和x6为邻频谐波的参数。
进一步地,步骤2包括:利用搜索得到的工频谐波基频f0,根据式(3)对工频及其谐波噪声建立模型Vharm
其中,Am分别是第m个谐波的幅度和相位,M是谐波个数,t为采样时间,t=[t1,t2,L,tK],对式(3)中的第m项进行线性化为:
其中,αm和βm为第m阶谐波对应系数,求解线性方程组,
得到谐波系数X=[α1,L,αM1L,βM]T,其中T是转置符号。
本发明与现有技术相比,有益效果在于:本发明公开的去除地面核磁共振测量数据中邻频谐波的方法,不需要参考线圈额外获取噪声数据。根据地面核磁共振信号和邻频谐波的幅度随时间变化的特征不同,利用非线性拟合方法同时估计两者的特征参数,进而实现对邻频谐波的建模和去除,且不损失有用信息。使得地面核磁共振方法可应用在拉莫尔频率和工频谐波的高阶谐波相近甚至重合的地点,拓宽了地面核磁共振方法的应用范围。
经试验,本发明方法对地面核磁共振数据中邻频谐波干扰的去除效率高,准确性好,提高了拉莫尔频率和工频谐波的高阶谐波相近的地点的探测结果的可信度。
附图说明
图1为地面核磁共振数据邻频谐波干扰去除方法的流程图;
图2地面核磁共振数据波形图和消除除邻频谐波外的其他谐波噪声的结果,(a)时域波形(b)频域波形;
图3希尔伯特变换和低通滤波后,地面核磁共振数据包络及拟合结果(a)实部(b)虚部;
图4地面核磁共振数据邻频谐波处理效果图(a)时域波形(b)频域波形。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。
以一组初始振幅为50nV,弛豫时间为0.3s,拉莫尔频率为2300Hz,初始相位为π/3rad的地面核磁共振信号;工频谐波基频为50.010Hz,100个谐波的幅度和相位呈随机分布;随机噪声为10nV的数据为例进行说明。
参见图1,步骤1,录入地面核磁共振数据V=[V1,V2,L,VK],K为采样点数,其时域和频域图见图2(a)和图2(b)黑色曲线所示。测量实验地点地磁场B0为54016nT,计算拉莫尔频率fL=0.04258×B0=2300Hz。
步骤2,利用搜索得到的工频谐波基频f0=50.010Hz,根据式(3)对工频及其谐波噪声建立模型:
其中,Am分别是第m个谐波的幅度和相位,M是谐波个数,t为采样时间,t=[t1,t2,L,tK]。对式(3)中的第m项进行线性化为:
其中,将f0=50.010Hz带入式(4),求解线性方程组
得到谐波系数X=[α1,L,αM1L,βM]T,其中T是转置符号。
步骤3,根据拉莫尔频率fL和每个谐波频率fm=m·f0的频率差dfm=fL-fm判断是否存在邻频谐波,其中第46阶谐波的频率为f46=n46×f0=2300.46Hz,|df46|=0.46Hz,小于3Hz,即存在邻频谐波。
步骤4,将谐波系数X中邻频谐波对应的系数置零,即α46=0,β46=0。带入式(3)得到V′harm,在测量数据V中减去V′harm得到数据V′,即实现除邻频外的工频及其谐波噪声消除,数据V′的时域和频域结果见图2中圆圈曲线。
对数据V′以邻频谐波频率f56=2300.46Hz进行希尔伯特变换,并经过低通滤波后,得到数据V′包络的实部和虚部(见图3(a)和图3(b)灰色曲线)分别为:
其中,e0和θ分别为地面核磁共振信号的初始幅度,弛豫时间和相位,df46=fL-f46=-0.46Hz是拉莫尔频率fL与邻波谐波频率f56的频率差。A56分别为邻频谐波的幅度和相位。εr和εi分别是随机噪声的实部和虚部分量。
根据表达式(1)和(2对测量数据的实部VRe和虚部VIm进行非线性拟合(拟合结果见图3中黑色圆圈曲线),
得到:x1=50.15(nV),x2=0.305(s),x3=-0.462(Hz),x4=1.048(rad),对应的地面核磁共振信号见图4(a)和图4(b)中黑色方框曲线,x5=43.35(nV)和x6=24.82(nV),对应的邻频谐波见图3中黑色虚线。得到邻频谐波的幅度和相位分别为:A56=49.95nV和带入式(3)得到邻频谐波模型V56,在数据V′中减去V56即完成了邻频谐波的消除,时域和频域结果见图4(a)和图4(b)中黑色曲线。
结果显示本发明方法对地面核磁共振数据中邻频谐波干扰的去除效率高,准确性好,提高了拉莫尔频率和工频谐波的高阶谐波相近的地点的探测结果的可信度。

Claims (4)

1.一种消除地面核磁共振数据中邻频谐波干扰的方法,其特征在于,该方法包括:
步骤1,录入地面核磁共振测量数据V=[V1,V2,…,VK],K为采样点数,测量实验地点地磁场B0,根据式fL=0.04258×B0计算拉莫尔频率fL
步骤2,搜索工频谐波基频f0,并对工频谐波进行建模,计算得到所有谐波系数X;
步骤3,根据拉莫尔频率fL和每个谐波频率fm=m·f0的频率差dfm=fL-fm判断是否存在邻频谐波,其中m=1,2,…,M,M是谐波个数,若所有的频率差|dfm|≥3Hz,说明不存在邻频谐波,利用步骤2得到的工频谐波基频f0 和谐波系数X计算得到谐波模型Vharm,在测量数据V中减去谐波模型Vharm,实现工频及其谐波干扰的消除;
步骤3中,若其中一个频率差|dfm|<3Hz,则存在邻频谐波,进行步骤4;
步骤4,将谐波系数X中邻频谐波对应的系数置零,进行谐波建模,得到除邻频谐波外的谐波模型V′harm,在测量数据V中减去谐波模型V′harm,首先完成除邻频谐波外所有工频谐波的消除,得到数据V′,基于拟合的方法,得到邻频谐波Vm的参数,在数据V′中减去邻频谐波Vm,得到去除谐波噪声后的地面核磁共振信号;
步骤4中基于拟合的方法包括对数据V′以邻频谐波的频率fm进行希尔伯特变换,并经过低通滤波器后得到数据V′包络的实部VRe和虚部VIm
对数据V′包络的实部VRe和虚部VIm进行非线性拟合获得6个未知参数x1~x6,其中四个参数对应于地面核磁共振信号参数,两个参数为邻频谐波的参数,获得邻频谐波的幅度和相位,计算得到邻频谐波Vm,在数据V′中减去邻频谐波Vm完成邻频谐波的消除。
2.按照权利要求1所述的方法,其特征在于,经过低通滤波器后得到V′包络的实部VRe和虚部VIm分别为:
其中e0和θ分别为地面核磁共振信号的初始幅度,弛豫时间和相位,dfm=fL-fm是拉莫尔频率fL与邻频谐波频率fm的频率差; Am分别为邻频谐波的幅度和相位,εr和εi分别是随机噪声的实部和虚部分量。
3.按照权利要求2所述的方法,其特征在于,根据表达式(1)和式(2)对数据V′包络的实部VRe和虚部VIm进行非线性拟合获得6个未知参数x1~x6
其中,分别为经拟合得到的实部和虚部,x1~x4对应于地面核磁共振信号参数,x5和x6为邻频谐波的参数。
4.按照权利要求1所述的方法,其特征在于,步骤2包括:利用搜索得到的工频谐波基频f0,根据式(3)对工频及其谐波噪声建立模型Vharm
其中,Am分别是第m个谐波的幅度和相位,M是谐波个数,t为采样时间,t=[t1,t2,…,tK],对式(3)中的第m项进行线性化为:
其中,αm和βm为第m阶谐波对应系数,求解线性方程组,
得到谐波系数X=[α1,…,αM1…,βM]T,其中T是转置符号。
CN201810029427.2A 2018-01-12 2018-01-12 一种消除地面核磁共振数据中邻频谐波干扰的方法 Active CN108345039B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810029427.2A CN108345039B (zh) 2018-01-12 2018-01-12 一种消除地面核磁共振数据中邻频谐波干扰的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810029427.2A CN108345039B (zh) 2018-01-12 2018-01-12 一种消除地面核磁共振数据中邻频谐波干扰的方法

Publications (2)

Publication Number Publication Date
CN108345039A CN108345039A (zh) 2018-07-31
CN108345039B true CN108345039B (zh) 2019-07-23

Family

ID=62961119

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810029427.2A Active CN108345039B (zh) 2018-01-12 2018-01-12 一种消除地面核磁共振数据中邻频谐波干扰的方法

Country Status (1)

Country Link
CN (1) CN108345039B (zh)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109782363B (zh) * 2019-02-15 2020-01-24 吉林大学 一种基于时域建模与频域对称的磁共振信号消噪方法
CN110018429B (zh) * 2019-03-29 2021-01-15 中国科学院电子学研究所 一种消除磁探测平台振动干扰磁场的方法和系统
CN110542925B (zh) * 2019-09-02 2020-07-28 吉林大学 一种基于峰值包络线的地震数据尖峰噪声识别及压制方法
CN111190233B (zh) * 2020-01-10 2021-08-27 吉林大学 一种基于展宽指数c的预极化场磁共振正反演方法
CN111239837B (zh) * 2020-02-20 2021-01-05 吉林大学 基于mcmc的地面磁共振信号参数提取方法
CN112235002B (zh) * 2020-12-14 2021-03-09 展讯通信(上海)有限公司 谐波消除方法及装置、存储介质、终端

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102508182A (zh) * 2011-11-30 2012-06-20 苏州安科医疗系统有限公司 一种用于磁共振成像的多阶谐波动态匀场方法
CN203502602U (zh) * 2013-10-05 2014-03-26 吉林大学 消除工频谐波干扰的核磁共振探测装置
CN104614778A (zh) * 2015-01-27 2015-05-13 吉林大学 基于ica的核磁共振地下水探测信号噪声消除方法
CN104796110A (zh) * 2015-05-05 2015-07-22 吉林大学 模拟工频梳状陷波器及其调整方法
CN104898172A (zh) * 2015-05-19 2015-09-09 吉林大学 一种基于互相关的核磁共振全波信号噪声滤除方法
CN105549097A (zh) * 2015-12-22 2016-05-04 吉林大学 一种瞬变电磁信号工频及其谐波干扰消除方法及装置
CN106772646A (zh) * 2016-12-19 2017-05-31 吉林大学 一种地面核磁共振信号提取方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9857443B2 (en) * 2013-03-15 2018-01-02 University Health Network Method and apparatus for the measurement, characterization and correction of geometric distortions in magnetic resonance imaging

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102508182A (zh) * 2011-11-30 2012-06-20 苏州安科医疗系统有限公司 一种用于磁共振成像的多阶谐波动态匀场方法
CN203502602U (zh) * 2013-10-05 2014-03-26 吉林大学 消除工频谐波干扰的核磁共振探测装置
CN104614778A (zh) * 2015-01-27 2015-05-13 吉林大学 基于ica的核磁共振地下水探测信号噪声消除方法
CN104796110A (zh) * 2015-05-05 2015-07-22 吉林大学 模拟工频梳状陷波器及其调整方法
CN104898172A (zh) * 2015-05-19 2015-09-09 吉林大学 一种基于互相关的核磁共振全波信号噪声滤除方法
CN105549097A (zh) * 2015-12-22 2016-05-04 吉林大学 一种瞬变电磁信号工频及其谐波干扰消除方法及装置
CN106772646A (zh) * 2016-12-19 2017-05-31 吉林大学 一种地面核磁共振信号提取方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
Separation of adjacent frequency components in Generalized Harmonic Analysis;Teruo Muraoka 等;《The 47th IEEE International Midwest Symposium on Circuits and Systems》;20041231;第II-157-160页
基于独立分量分析的同频、邻频及谐波干扰消除法;金贵斌 等;《电测与仪表》;20090630;第46卷(第522期);第1-4页
基于高阶累积量的MRS信号噪声滤除方法研究;王悦;《中国优秀硕士学位论文全文数据库 基础科学辑》;20160915(第09(2016)期);第A011-204页

Also Published As

Publication number Publication date
CN108345039A (zh) 2018-07-31

Similar Documents

Publication Publication Date Title
CN108345039B (zh) 一种消除地面核磁共振数据中邻频谐波干扰的方法
Markalous et al. Detection and location of partial discharges in power transformers using acoustic and electromagnetic signals
Yao et al. Fast S-transform for time-varying voltage flicker analysis
CN106772646B (zh) 一种地面核磁共振信号提取方法
CN103823244B (zh) 磁共振三分量消噪装置及消噪方法
Liu et al. Correlation analysis for spread-spectrum induced-polarization signal processing in electromagnetically noisy environments
Chen et al. An ameliorated synchroextracting transform based on upgraded local instantaneous frequency approximation
CN109581516B (zh) 曲波域统计量自适应阈值探地雷达数据去噪方法及系统
CN112083509B (zh) 时频电磁法中激发极化异常的检测方法
CN106646637A (zh) 一种去除核磁信号中尖峰噪声的方法
Jiang et al. Harmonic noise-elimination method based on the synchroextracting transform for magnetic-resonance sounding data
Cheng et al. Enhanced symplectic characteristics mode decomposition method and its application in fault diagnosis of rolling bearing
Kremer et al. Review of acquisition and signal processing methods for electromagnetic noise reduction and retrieval of surface nuclear magnetic resonance parameters
CN112487604A (zh) 海洋重力仪输出数据长时间非线性漂移补偿方法
CN104849590A (zh) 一种混合噪声干扰下微弱脉冲信号检测方法
CN110308332A (zh) 一种变电站接地网接地阻抗检测方法、系统及介质
Ma et al. An improved intrinsic time-scale decomposition method based on adaptive noise and its application in bearing fault feature extraction
Lan et al. Partial discharge location of power cables based on an improved phase difference method
Yan et al. Feature extraction by enhanced time–frequency analysis method based on Vold-Kalman filter
Xie et al. Suppressing the discrete spectral interference of the partial discharge signal based on bivariate empirical mode decomposition
CN108254794A (zh) 一种基于建模反恢复技术的磁共振消噪方法及装置
CN110456420B (zh) 基于近端参考线圈核磁共振地下水探测信号噪声消除方法
Xin et al. De-noising of GIS UHF Partial Discharge monitoring based on wavelet method
Li et al. Research on the detection method of MRS Signal initial amplitude based on chaotic detection system
Ma et al. Adaptive coherent noise suppression for three-axis fluxgate sensor

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