CN104535496B - 相位型表面等离子共振传感系统中的干涉光信号处理方法 - Google Patents

相位型表面等离子共振传感系统中的干涉光信号处理方法 Download PDF

Info

Publication number
CN104535496B
CN104535496B CN201410817251.9A CN201410817251A CN104535496B CN 104535496 B CN104535496 B CN 104535496B CN 201410817251 A CN201410817251 A CN 201410817251A CN 104535496 B CN104535496 B CN 104535496B
Authority
CN
China
Prior art keywords
signal
waveform
arcsin
value
phase
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
CN201410817251.9A
Other languages
English (en)
Other versions
CN104535496A (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.)
Dalian University of Technology
Original Assignee
Dalian University of 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 Dalian University of Technology filed Critical Dalian University of Technology
Priority to CN201410817251.9A priority Critical patent/CN104535496B/zh
Publication of CN104535496A publication Critical patent/CN104535496A/zh
Application granted granted Critical
Publication of CN104535496B publication Critical patent/CN104535496B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/17Systems in which incident light is modified in accordance with the properties of the material investigated
    • G01N21/55Specular reflectivity
    • G01N21/552Attenuated total reflection
    • G01N21/553Attenuated total reflection and using surface plasmons
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/17Systems in which incident light is modified in accordance with the properties of the material investigated
    • G01N21/41Refractivity; Phase-affecting properties, e.g. optical path length
    • G01N21/45Refractivity; Phase-affecting properties, e.g. optical path length using interferometric methods; using Schlieren methods

Landscapes

  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Investigating Or Analysing Materials By Optical Means (AREA)
  • Length Measuring Devices By Optical Means (AREA)

Abstract

本发明提供了一种相位型表面等离子共振传感系统中的干涉光信号处理方法,属于光信号处理技术领域。搭建基于马赫—泽德干涉仪结构的光路,He‑Ne激光器发出的激光经过表面等离子共振传感层后,激光携带着传感信号。在光路的输出端采用光电探测器将光信号接收,然后通过数据采集板将信号输入电脑。结合数学理论和信号处理知识,自主设计并编写算法将得到的信号精确地处理,进而解调出了高精度、高稳定性的相位信息。本发明的效果和益处是发明了一种新的检测和处理相位型表面等离子共振传感信号的方法,并实现了传感器的高精度测量,使得信号采集和信号处理实时进行,抑制了相位的漂移,提高了采集系统的角度分辨率。

Description

相位型表面等离子共振传感系统中的干涉光信号处理方法
技术领域
本发明属于光信号处理技术领域,涉及到光-电信号的转换,信号的提取、锁定、预处理、精确处理的方法,特别涉及到经典马赫—泽德干涉仪结构的光路装置。
背景技术
表面等离子体共振传感技术是一种高精度的光学折射率传感技术,以高灵敏度、所需待测样品少、响应速度快等优势,被广泛应用于化学和生物的传感应用领域中。传统的表面等离子体共振传感器主要基于对光的振幅信息的检测来实现折射率传感,包括强度型、角度型和波长型等几种类型。新型的相位型表面等离子体共振传感器基于相位信息检测,其传感分辨率优于传统的振幅检测的传感器2-3个数量级,因此在生物小分子相互作用和微量样品检测中有着巨大的应用潜力,其传感器构架和相位信号提取方式也是各种各样的。常见的几种提取方式为:光偏振测定法(检测激发光经过表面等离子体共振传感面发射前后的偏振状态改变来提取表面等离子体相位信号)、光外差法(通过光学拍频将电场振荡频率降低到几十kHz的低频振荡同时保留相位信息,然后用相位计读取低频振荡的相位信息)、光干涉法(检测表面等离子体共振中信号光的相位)。无论采取何种提取方式,其目的都是为了实现相位型传感器的高分辨率。
发明内容
本发明的目的是为了提高相位型表面等离子体共振传感器的分辨率。本发明要解决的技术问题是精确提取经过传感器两路干涉光信号之间的相位差信息。
本发明的技术方案是:
相位型表面等离子共振传感系统中的干涉光信号处理方法,
结合数学理论和信号处理知识,基于马赫—泽德干涉仪结构的光路,He-Ne激光器发出的激光经过表面等离子共振传感层后,会携带着传感信号。在光路的输出端采用光电探测器将光信号接收,然后通过数据采集板将信号输入电脑。自主设计并编写算法将得到的信号连续地经过提取、锁定、预处理、精确处理等步骤,进而解调出了高精度、高稳定性的相位信息。
相位型表面等离子共振传感系统中的干涉光信号处理方法,步骤如下:
1、提取信号
信号发生器中产生的三角波作为激励信号,将其采入系统,若第一个三角波波谷在波峰前面,则直接选取第一个波谷三角波到第一个三角波波峰之间对应的正弦波形。若第一个三角波波峰在第一个波谷三角波前面,则选取第一个三角波波谷和第二个三角波波峰之间对应的正弦波形,并选取三角波,采入系统中的第一个三角波完整上升沿所对应的正弦波形(干涉信号)作为步骤2处理的信号。循环采入三角波信号,并得到一系列干涉信号。
2、信号锁定
为了使每次采入系统中的干涉信号都在同一个相位处,确定锁定信号的方法使得每次开始采集干涉信号的位置随着随机直流分量的变化而变化,从而保证初始的相位是确定值。锁定信号的方法:找到步骤1中所选取的一列正弦波,并确定该列正弦波中接近中间位置的一个波峰。第一次循环结束后,记录上述波峰对应的初相值,以后每次循环均以第一次循环时记录的初相值为基准调整开始循环的位置,保证了每次循环的初始相位在较小的区间内变动,可得到多个正弦周期。例:本发明中,在步骤1中共选取了13个正弦波形,选取第6个正弦波形波峰的位置为基准位置,然后以这个波峰位置为基准,找到此波峰位 置向前40个点位置。
3、信号的预处理
步骤2得到的信号波形包含多个正弦周期,在同一上升沿,压电陶瓷依然有非线性和系统噪声的干扰,使得采集的各正弦周期的波峰和波谷值依然具有细微的偏差。记录得到的各正弦周期的波峰值和波谷值,然后对正弦周期进行采样和插值,得到信号波形的上下包络面(上包络面对应波峰值,下包络面对应波谷值),再根据上下包络面的值将波形归一化为标准正弦波形。
4、精确处理信号以及对相位的提取
将归一化后的标准正弦波形逐点进行反正弦运算,并将运算的结果存入数组中,将数组中的数据持续同时进行以下运算:
1)将两路标准正弦波形进行反正弦运算,并逐点做差并求绝对值,其中一路标准正弦波形反正弦运算后用arcsin p表示,另一路标准正弦波形反正弦运算后用arcsin s表示,表示计算得到的相位差值。
当arcsin p和arcsin s处于相同的单调变化趋势时
2)将两路标准正弦波形进行反正弦运算,并逐点求和,再用π减去和值,
当arcsin p和arcsin s处于不同的单调变化趋势时
3)将步骤1)和步骤2)求得的图像合并,并得到合并图像后的水平线。
4)将步骤3)水平线上的点求平均值作为最后的相位值,
步骤1)和步骤2)求得的图像在水平位置恰好互补,随着测量值的变化,参与计算的点数并不会发生变化,结果会保持很好的稳定性。对于两部分图像,应用不同的方法分别取其水平位置。对于arcsin p和arcsin s处于相同的单调变化趋势的曲线,选取水平部分的方法:先将整个波形中的峰值选出,再找到峰值中的最小值,将大于最小值的波形中的点组成水平部分;对于arcsin p和arcsin s处 于不同的单调变化趋势曲线,选取水平部分的方法:先将波形中的谷值选出,再找到谷值中的最大值,将小于最大值的波形中的点组成水平部分,最后将两路的水平部分合在一起处理。
本发明的效果和益处是发明了一种新的检测和处理相位型表面等离子共振传感信号的方法,并实现了传感器的高精度测量,使得信号采集和信号处理实时进行,抑制了相位的漂移,提高了采集系统的的角度分辨率。
附图说明
图1是本传感器光路图。
图2是初步提取到的两路光信号及其包络图。
图3是预处理两路光信号图。
图4是精确处理信号后将作和的和作差的叠加图。
图5是传感器实时测试性能时得到的时间—相位差图。
图6是传感器实时测试性能时得到的折射率—相位差图。
图7是精确处理信号后的高精度的相位信息图。
图中:1 He-Ne激光器;2 50×50的分光镜;3 三层结构的棱镜型传感层;
4 50×50的分光镜;5 沃拉斯顿棱镜;6 光电探测器;
7 USB-6210数据采集板;8 运算处理信号的PC机;
9 信号发生器驱动着压电陶瓷的反射镜;10 提取得到的p偏振光分量;
11 提取得到的s偏振光分量;12 两路光信号的包络;
13 精确处理信号以及对相位的提取步骤中两部分图像的水平位置。
具体实施方式
以下结合附图和技术方案,进一步说明本发明的具体实施方式。
He-Ne激光器发出线偏振光,经过分光镜将光分为两束,一束打到有三层结 构的传感层表面,其中P偏振光的相位对所加样品折射率敏感,而S偏振光的相位基本无变化;另一束光经粘贴着压电陶瓷的反射镜反射,压电陶瓷由有固定频率和电压的驱动器驱动。最后两路光信号在第二个分光镜处汇聚成为一束光,同时发生干涉。干涉光经渥拉斯顿棱镜将p偏振光和s偏振光分开,采用光电探测器将光信号接收,然后通过USB-6210数据采集板将电信号输入电脑,同时采集板也将驱动压电陶瓷的信号输入电脑。
将采集到的干涉光信号分为几个步骤进行处理:
1、提取信号
由于驱动产生干涉光信号的电器元件(压电陶瓷)的非线性和驱动电压直流分量的不稳定性导致采集到的干涉信号波形的疏密程度,波峰(谷)值是不一样的。因为压电陶瓷具有非线性,在激励信号三角波的上升沿压电陶瓷长度随电压变化的曲线与三角波下降沿压电陶瓷长度随电压的变化曲线并不重合且均不为直线。如果后续算法所用的波形同时包括有上升沿和下降沿,则会产生较大误差。信号发生器中产生的三角波作为激励信号,将其采入系统,(若第一个三角波波谷在波峰前面,则直接选取第一个波谷三角波到第一个三角波波峰之间对应的正弦波形。若第一个三角波波峰在第一个波谷三角波前面,则选取第一个三角波波谷和第二个三角波波峰之间对应的正弦波形)并选取三角波,采入系统中的第一个三角波完整上升沿所对应的正弦波形(干涉信号)作为步骤2处理的信号。循环采入三角波信号,并得到一系列干涉信号。通过此处理过程就可保证后续参与计算的干涉信号不受压电陶瓷上升沿和下降沿参数不同而导致的误差。
2、信号锁定
为了使每次采入系统中的干涉信号都在同一个相位处,确定锁定信号的方 法使得每次开始采集干涉信号的位置随着随机直流分量的变化而变化,从而保证初始的相位是确定值。锁定信号的方法:找到步骤1中所选取的一列正弦波,并确定该列正弦波中接近中间位置的一个波峰。第一次循环结束后,记录上述波峰对应的初相值,以后每次循环均以第一次循环时记录的初相值为基准调整开始循环的位置,保证了每次循环的初始相位在较小的区间内变动,可得到多个正弦周期。(例:本发明中,在步骤1中共选取了13个正弦波形,选取第6个正弦波形波峰的位置为基准位置),然后以这个波峰位置为基准,找到此波峰位置向前40个点位置)
3、信号的预处理
步骤2得到的信号波形包含多个正弦周期,在同一上升沿,压电陶瓷依然有非线性和系统噪声的干扰,使得采集的各正弦周期的波峰和波谷值依然具有细微的偏差。记录得到的各正弦周期的波峰值和波谷值,然后对正弦周期进行采样和插值,得到信号波形的上下包络面(上包络面对应波峰,下包络面对应波谷值),再根据上下包络面的值将波形归一化为标准正弦波形。
4、精确处理信号以及对相位的提取
将归一化后的标准正弦波形逐点进行反正弦运算,并将运算的结果存入数组中,将数组中的数据持续同时进行以下运算:
1)将两路标准正弦波形进行反正弦运算,并逐点做差并求绝对值,其中一路标准正弦波形反正弦运算后用arcsin p表示,另一路标准正弦波形反正弦运算后用arcsin s表示,表示计算的到的相位差值。
当arcsin p和arcsin s处于相同的单调变化趋势时
2)将两路标准正弦波形进行反正弦运算,并逐点求和,再用π减去和值,
当arcsin p和arcsin s处于不同的单调变化趋势时
3)将步骤1)和步骤2)求得的图像合并,并得到合并图像后的水平线。
4)将步骤3)水平线上的点求平均值作为最后的相位值,
步骤1)和步骤2)求得的图像在水平位置恰好互补,随着测量值的变化,参与计算的点数并不会发生变化,结果会保持很好的稳定性。对于两部分图像,应用不同的方法分别取其水平位置。对于arcsin p和arcsin s处于相同的单调变化趋势的曲线,选取水平部分的方法:先将整个波形中的峰值选出,再找到峰值中的最小值,将大于最小值的波形中的点组成水平部分;对于arcsin p和arcsin s处于不同的单调变化趋势曲线,选取水平部分的方法:先将波形中的谷值选出,再找到谷值中的最大值,将小于最大值的波形中的点组成水平部分,最后将两路的水平部分合在一起处理。从而利用该系统和信号处理算法将得到的光—电信号精确地处理,解调出了高精度、高稳定性的相位信息。

Claims (1)

1.一种相位型表面等离子共振传感系统中的干涉光信号处理方法,其特征在于,步骤如下:
基于马赫—泽德干涉仪结构的光路,He-Ne激光器发出线偏振光,经过分光镜将光分为两束,一束打到有三层结构的传感层表面,其中P偏振光的相位对所加样品折射率敏感,而S偏振光的相位基本无变化;另一束光经粘贴着压电陶瓷的反射镜反射,压电陶瓷由有固定频率和电压的驱动器驱动;最后两路光信号在第二个分光镜处汇聚成为一束光,同时发生干涉;干涉光经渥拉斯顿棱镜将p偏振光和s偏振光分开,采用光电探测器将光信号接收,然后通过USB-6210数据采集板将电信号输入电脑,同时采集板也将驱动压电陶瓷的信号输入电脑;
(1)提取信号
信号发生器中产生的三角波作为激励信号,将其采入系统,选取的原则:若第一个三角波波谷在波峰前面,直接选取第一个波谷三角波到第一个三角波波峰之间对应的正弦波形;若第一个三角波波峰在第一个波谷三角波前面,选取第一个三角波波谷和第二个三角波波峰之间对应的正弦波形,并选取三角波,采入系统中的第一个三角波完整上升沿所对应的正弦波形即干涉信号作为步骤(2)处理的信号;循环采入三角波,得到一系列干涉信号;
(2)信号锁定
确定步骤(1)中所选取的一列正弦波,并确定该列正弦波中接近中间位置的一个波峰;第一次循环结束后,记录上述波峰对应的初相值,每次循环均以第一次循环时记录的初相值为基准调整开始循环的位置,保证每次循环的初始相位在较小的区间内变动,得到多个正弦周期;
(3)信号预处理
记录步骤(2)得到的各正弦周期的波峰值和波谷值,对正弦周期进行采样和插值,得到信号波形的上下包络面,上包络面对应波峰值,下包络面对应波谷值,再根据上下包络面的值将波形归一化为标准正弦波形;
(4)精确处理信号以及对相位的提取
将归一化后的标准正弦波形逐点进行反正弦运算,并将运算的结果存入数组中,将数组中的数据持续同时进行以下运算:
1)将两路标准正弦波形进行反正弦运算,并逐点做差并求绝对值,其中一路标准正弦波形反正弦运算后用arcsin p表示,另一路标准正弦波形反正弦运算后用arcsin s表示,表示计算得到的相位差值;
当arcsin p和arcsin s处于相同的单调变化趋势时
2)将两路标准正弦波形进行反正弦运算,并逐点求和,再用π减去和值,当arcsin p和arcsin s处于不同的单调变化趋势时
3)将步骤1)和步骤2)求得的图像合并,并得到合并图像后的水平线;
4)将步骤3)水平线上的点求平均值作为最后的相位值,
对于arcsin p和arcsin s处于相同的单调变化趋势的曲线,选取水平部分的方法:先将整个波形中的峰值选出,再找到峰值中的最小值,将大于最小值的波形中的点组成水平部分;对于arcsin p和arcsin s处于不同的单调变化趋势曲线,选取水平部分的方法:先将波形中的谷值选出,再找到谷值中的最大值,将小于最大值的波形中的点组成水平部分,最后将两路的水平部分合在一起处理。
CN201410817251.9A 2014-12-24 2014-12-24 相位型表面等离子共振传感系统中的干涉光信号处理方法 Active CN104535496B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410817251.9A CN104535496B (zh) 2014-12-24 2014-12-24 相位型表面等离子共振传感系统中的干涉光信号处理方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410817251.9A CN104535496B (zh) 2014-12-24 2014-12-24 相位型表面等离子共振传感系统中的干涉光信号处理方法

Publications (2)

Publication Number Publication Date
CN104535496A CN104535496A (zh) 2015-04-22
CN104535496B true CN104535496B (zh) 2017-04-12

Family

ID=52851070

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410817251.9A Active CN104535496B (zh) 2014-12-24 2014-12-24 相位型表面等离子共振传感系统中的干涉光信号处理方法

Country Status (1)

Country Link
CN (1) CN104535496B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106770032A (zh) * 2016-12-07 2017-05-31 大连理工大学 相位探测型光纤spr生化传感系统
CN114285321B (zh) * 2021-12-25 2023-08-22 西安交通大学 一种基于共振驱动可多自由度钻行的大位移贯入式作动器
CN114614892B (zh) * 2022-05-11 2023-03-31 南京信息工程大学 一种基于等离子体光谱的光学信息传输系统

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2002001146A1 (en) * 2000-06-28 2002-01-03 The Charles Stark Draper Laboratory, Inc. Optical microcavity resonator sensor
CN1588064A (zh) * 2004-08-27 2005-03-02 清华大学 空间相位调制干涉阵列检测生物芯片的方法及系统
CN101398378A (zh) * 2008-01-28 2009-04-01 国家纳米科学中心 一种表面等离子共振的相位测量方法及其测量系统
CN101449135A (zh) * 2006-04-17 2009-06-03 比奥普蒂克斯有限责任公司 基于偏振的干涉测量检测器
CN101477044A (zh) * 2009-01-19 2009-07-08 浙江大学 一种表面等离子共振传感器
TW201445128A (zh) * 2013-05-30 2014-12-01 Univ Nat Cheng Kung 具不對稱粒子形狀之定域化表面電漿共振檢測系統

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7027676B2 (en) * 2004-03-08 2006-04-11 Agilent Technologies, Inc. Optical phase measurement of target
TWI426260B (zh) * 2010-05-04 2014-02-11 Univ Nat Central Surface plasmon resonance measurement device
US8792102B2 (en) * 2010-10-28 2014-07-29 General Electric Company Interferometric spectral imaging of a two-dimensional array of samples using surface plasmon resonance

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2002001146A1 (en) * 2000-06-28 2002-01-03 The Charles Stark Draper Laboratory, Inc. Optical microcavity resonator sensor
CN1588064A (zh) * 2004-08-27 2005-03-02 清华大学 空间相位调制干涉阵列检测生物芯片的方法及系统
CN101449135A (zh) * 2006-04-17 2009-06-03 比奥普蒂克斯有限责任公司 基于偏振的干涉测量检测器
CN101398378A (zh) * 2008-01-28 2009-04-01 国家纳米科学中心 一种表面等离子共振的相位测量方法及其测量系统
CN101477044A (zh) * 2009-01-19 2009-07-08 浙江大学 一种表面等离子共振传感器
TW201445128A (zh) * 2013-05-30 2014-12-01 Univ Nat Cheng Kung 具不對稱粒子形狀之定域化表面電漿共振檢測系統

Also Published As

Publication number Publication date
CN104535496A (zh) 2015-04-22

Similar Documents

Publication Publication Date Title
CN110779464B (zh) 一种时域频域联合分析宽光谱相干测量方法及系统
US7407817B2 (en) Surface plasmon resonance sensors and method for detecting samples using the same
CN102095538B (zh) 保偏光纤应力传感的数据解调方法
CN106842095B (zh) 基于里德堡原子量子相干效应的射频电场校对方法及装置
CN110487313A (zh) 光频域反射技术中光源扫频非线性自校正方法
CN104677398B (zh) 一种基于频漂动态补偿的φ‑otdr增敏方法及装置
CN104535496B (zh) 相位型表面等离子共振传感系统中的干涉光信号处理方法
CN103604774A (zh) 基于非线性调谐提高激光气体分析灵敏度的方法和装置
CN104748671B (zh) 角位移单频激光干涉仪非线性误差修正方法及装置
CN109297581A (zh) 一种用于补偿相位敏感光时域反射计中频率漂移的二次相位差分测量方法
CN110702210A (zh) 光纤迈克尔逊干涉测振仪非线性误差修正方法
CN112946611A (zh) 基于相似三角插值采样的扫频非线性矫正测距方法
CN110048765A (zh) 一种基于整体最小二乘拟合的φ-otdr定量测量的方法
CN104061998A (zh) 漫反射式零差正交激光测振仪
CN110701998A (zh) 光纤迈克尔逊干涉仪非线性误差修正方法
CN105911605B (zh) 一种光学干涉式重力仪中的闭环信号采集方法
CN106124032A (zh) 一种光频调制器调制延迟的数字化测量方法及装置
CN108489959A (zh) 一种相干反斯托克斯拉曼光谱扫描装置和方法
Svitlov et al. Comparison of three digital fringe signal processing methods in a ballistic free-fall absolute gravimeter
CN109883350A (zh) 一种异形曲面结构内部形貌的高精度测量系统和测量方法
CN205581298U (zh) 一种基于f-p标准具的高精度调频连续波激光测距系统
Jiang et al. Improved transition detection algorithm for a self-mixing displacement sensor
CN104655029B (zh) 一种位相增强型薄膜厚度测量方法和系统
CN116610908A (zh) 一种基于多层感知机的偏振调制测距系统及方法
WO2019218166A1 (zh) 一种激光器的偏置电流控制方法及装置

Legal Events

Date Code Title Description
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant