CN114326358B - 一种多链并行分割高精度fpga时间数字转换方法 - Google Patents

一种多链并行分割高精度fpga时间数字转换方法 Download PDF

Info

Publication number
CN114326358B
CN114326358B CN202111559443.0A CN202111559443A CN114326358B CN 114326358 B CN114326358 B CN 114326358B CN 202111559443 A CN202111559443 A CN 202111559443A CN 114326358 B CN114326358 B CN 114326358B
Authority
CN
China
Prior art keywords
delay
chain
equivalent
delay chain
units
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
CN202111559443.0A
Other languages
English (en)
Other versions
CN114326358A (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.)
Shanghai Institute of Optics and Fine Mechanics of CAS
Original Assignee
Shanghai Institute of Optics and Fine Mechanics of CAS
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 Shanghai Institute of Optics and Fine Mechanics of CAS filed Critical Shanghai Institute of Optics and Fine Mechanics of CAS
Priority to CN202111559443.0A priority Critical patent/CN114326358B/zh
Publication of CN114326358A publication Critical patent/CN114326358A/zh
Application granted granted Critical
Publication of CN114326358B publication Critical patent/CN114326358B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Measurement Of Unknown Time Intervals (AREA)

Abstract

一种多链并行分割高精度FPGA时间数字转换方法,将待测脉冲输入N条并行的延迟链,设待测脉冲在第i个延迟链中从输入到被捕捉经过τi个延迟单元;将每条延迟链得到的延迟单元τi相加,作为待测脉冲在等效延迟链中经过的延迟单元τ,即将位置不同、各延迟单元不同的并行延迟链相互分割,最终合并得到一条等效延迟链。本发明合并得到的等效延迟链改善了单条延迟链延迟单元不均匀、有较大延迟单元的问题,使TDC的测量分辨率和测量精度得到明显提升,并且简化了电路,有效降低了环境温度变化带来的影响。

Description

一种多链并行分割高精度FPGA时间数字转换方法
技术领域
本发明涉及高精度时间测量,尤其涉及一种应用在现场可编程逻辑器件FPGA的TDC高精度实现方法。可应用于时间传输、基础研究等领域。
背景技术
时间测量是现代社会各种重要技术的基础,常用于时间传输、时钟比对、定时同步、精密测量、光通信等领域。高精度的时间间隔测量技术,尤其是皮秒量级的技术更为重要。目前,应用最广泛的高精度时间间隔测量技术是时间数字转换技术(time to digitalconverter,TDC),相比于ASIC的TDC,FPGA-TDC有着开发周期短、设计灵活等优点。在基于FPGA的TDC设计中,使用加法进位链构成的延迟链对时钟周期进行内插的方法已被广泛使用。然而基于加法进位链构成的延迟链有以下两个特点:延迟单元的不一致性和存在较大延迟单元。其中延迟单元的不一致性是由FPGA制造工艺决定的,较大延迟单元是由多余走线和跨时钟区域造成的。上述特点是限制TDC测量精度的主要因素。另外,环境温度和工作电压的改变也会影响内插单元的延时。
为了提高FPGA-TDC的测量精度,人们提出了多种解决方案。在先技术一:Wu,J.,and Z.Shi."The 10-Ps Wave Union Tdc:Improving Fpga Tdc Resolution Beyond ItsCell Delay."Paper presented at the IEEE Nuclear Science Symposium ConferenceRecord,2009.提出了wave union的方法来提高测量精度。该方法通过将待测信号转换为相应的周期信号,从而在单条延迟链上实现多次测量。由于每次测量的起始位置不同,所以每次落在的延迟单元也不同。这样通过在单条延迟链上进行多次测量就实现了不同延迟单元间的切割,从而将较大延迟单元切割成小的延迟单元。因此突破了延迟单元的固有物理时延,使TDC的分辨率和测量精度都得到了有效的提高。但是该方法会引入较大的死时间,因为每次测量至少花费一个时钟周期,那么N次测量就需要至少N个时钟周期。该结构的TDC在一些要求高吞吐量的场景下无法使用,并且其译码电路比较复杂。
在先技术二:Qin,X.,L.Wang,D.Liu,Y.Zhao,X.Rong,and J.Du."A1.15-Ps BinSize and 3.5-Ps Single-Shot Precision Time-to-Digital Converter with on-BoardOffset Correction in an Fpga."IEEE Transactions on Nuclear Science 64,no.12(2017):2951-57.提出了多链并行测量求平均的方法。该方法将待测信号同时输入N条延迟链,把得到的N个时间戳求平均作为最终的细测量结果,理论上分辨率可以提高倍。该方法可以在一个时钟周期内完成N次测量,在改善TDC分辨率的同时兼顾了死时间。但是该结构TDC的多条延迟链带来大的资源开销,除此之外还需要复杂的校准模块和计算模块。校准模块不仅要校准、存储每条延迟链的真实时延,还需要校准不同链间的走线时延,最后通过计算模块来得到平均值。
发明内容
为克服上述现有技术的不足,本发明提出了一种多链并行分割高精度FPGA时间数字转换方法,将多条并行延迟链分割合并为一条等效延迟链,使用该等效延迟链对一个时钟周期进行内插,从而实现高精度的TDC。
本发明原理如下:
使用原语将FPGA内部相邻专用进位单元的cout与cin相连,由此级联形成延迟链。每个cout作为D触发器的输入,当系统时钟上升沿到来时,D触发器采样得到此时cout的数值,由此得到一串温度计码111..000,其中1的个数即为待测脉冲上升沿从进入延迟链到下一个系统时钟上升沿这段时间经过的延迟单元个数。
本发明的技术解决方案如下:
一种多链并行分割高精度FPGA时间数字转换方法,其特点在于,
将待测脉冲输入N条并行的延迟链,设待测脉冲在第i个延迟链中从输入到被捕捉经过τi个延迟单元;
将每条延迟链得到的延迟单元τi相加,作为待测脉冲在等效延迟链中经过的延迟单元τ,即
将位置不同、各延迟单元不同的并行延迟链相互分割,最终合并得到一条等效延迟链。使用该等效延迟链对一个时钟周期内插,从而提高TDC的测量精度。等效延迟链的延迟单元约为单条延迟链的N倍,因此平均分辨率相较于单链TDC提高了N倍。
所述延迟链由FPGA的专用进位单元级联构成。
将相邻的专用进位单元cin与cout相连,由此级联形成所述的延迟链,且每个专用进位单元cout作为D触发器的输入,当系统时钟上升沿到来时,D触发器对cout数据完成采样。
所述的等效延迟链经过校准阶段后将得到延迟单元信息构建为查找表存储在RAM中;在测量阶段得到待测脉冲在等效延迟链中经过τ个延迟单元后,查表即可得到细测量时间。
本发明的特点和优点是:
1)本发明将N条延迟链分割合并成一条等效延迟链,不仅对较大延迟单元完成了分割,而且等效延迟单元增加了N倍。这样通过等效延迟链对一个时钟周期进行内插,使其平均分辨率增强了N倍,测量精度也得到有效提升。
2)本发明相对于在先技术一,在打破延迟单元物理时延,提高分辨率和测量精度的同时,没有引入多余的死时间。相对于在先技术二,只需校准、存储等效延迟链的真实时延,而且无需后续的求平均计算。大大简化了电路,节约了资源使用。
3)本发明得益于多链间的相互分割,等效延迟链有了颗粒度更小且数量更多的延迟单元,使得TDC分辨率对环境温度变化不敏感。
附图说明
图1是本发明的多链分割合并为等效延迟链的原理图;
图2是本发明的电路功能框图。
图3是本发明中粗细测量时序图
具体实施方式
下面结合实施例和附图对本发明作进一步说明,但不应以此限制本发明的保护范围。
先请参阅图1,图1是本发明中多链分割合并为等效延迟链的原理图,由图可见,不同的矩形代表不同宽度的延迟单元。延迟链的延迟单元是不均匀的,两条延迟链之间也不相同。另一方面,不同的延迟链物理位置不同,则输入信号到每条延迟链之间的走线时延也是不同的。所以同一个输入信号在不同的延迟链上经过的延迟单元数是不同的。如图1所示,输入信号上升沿落在第一条链的第3个延迟单元,而落在第二条链的第2个延迟单元。将两条链的延迟单元数相加作为等效延迟链的延迟单元数,由此得到一条等效延迟链,则输入信号上升沿落在等效延迟链的第5个延迟单元。可以看到,两条链互相分割合并后的等效延迟链,其平均分辨率提高了两倍,且延迟链1中的第三个较大延迟单元也得到了有效分割。由此可知,将N条延迟链分割合并为一条等效延迟链后,其延迟单元为单链的N倍,则平均分辨率提高了N倍。同时单链里存在的较大延迟单元也被分割为正常大小,测量精度得到有效提升。链间分割同样改善了由环境温度改变对单链延迟单元大小造成的影响。
本发明的实施例如下:
利用本发明提出的多链并行分割法在赛灵思FPGA上实现双通道TDC,TDC框图如图2所示,将待测脉冲1和待测脉冲2分别输入通道1和通道2,通道1和通道2的结构相同。以通道1来进行说明,使用原语调用carry4单元(如果是intel FPGA,则为LE的进位链单元),将前后carry4单元的cout和cin相连,实现将多个carry4单元级联为延迟链。延迟链的长度由系统时钟决定。延迟链的每个抽头接到D触发器的输入端,构成了触发器阵列。当系统时钟上升沿到来,触发器阵列对延迟链抽头进行采样,得到形如111…000的温度计码,再将温度计码输入译码电路,把1的个数转换为二进制值,也即是传输经过延迟单元的个数。使用同样的方法构造N条这样的延迟链,使其在FPGA芯片上并行排列。校准模块得到等效延迟链,把等效延迟链延迟单元大小的查找表存放在RAM中。粗计数模块由计数器构成,计数整周期的个数,实现大范围的时延测量。将粗测量结果和细测量结果都存储在RAM中,当测量结束后,将粗测量结果和细测量结果组合得到最终的时延测量结果。
所述的校准模块将译码电路得到的N个值相加作为等效延迟链的值。在正式测量前,使用码密度法校准得到等效延迟链每个延迟单元的大小。由独立晶振产生的随机脉冲信号在延迟链上传输。对于一个校准信号上升沿,在不同延迟链上被不同延迟单元捕捉,将落在每条延迟链上的延迟单元数相加作为落在等效延迟链上的延迟单元数。经过大量随机脉冲信号校准后,通过统计方法得到等效延迟链每个延迟单元的大小,即其中n为校准脉冲信号的次数,k为落在每个延迟单元bin上的次数。由此便可得到等效延迟链每一个延迟单元的大小,将每个延迟单元对应的延迟大小作为查找表存放在RAM中,完成对等效延迟链的校准。
校准完成后即可进行时延测量,将待测脉冲1和待测脉冲2分别接入TDC的通道1和通道2。当通道1细测量有效信号置高时,将细测量值和此时粗计数值n1存入ram。当通道2细测量有效信号置高时,将细测量值和此时粗计数值n2存入RAM。细测量模块则记录了脉冲上升沿到下一个时钟上升沿之间经过了多少个延迟单元,由存储在RAM中等效延迟链延迟单元大小的查找表即可得知细测量时间t1和t2。图3是粗细测量时序图,如图所示,两个待测脉冲间的时间间隔即为TI=nTclk+t1-t2,其中n=n2-n1
本发明将多条加法进位链采用链间分割合并成一条等效延迟链,实现对一个时钟周期的内插,完成细时间测量。合并得到的等效延迟链改善了单条延迟链延迟单元不均匀、有较大延迟单元的问题,使TDC的测量分辨率和测量精度得到明显提升,并且简化了电路,有效降低了环境温度变化带来的影响。本发明属于时间测量技术领域,可应用于时间传输、时钟比对、定时同步、精密测量、光通信等领域。

Claims (3)

1.一种多链并行分割高精度FPGA时间数字转换方法,其特征在于,
将待测脉冲输入N条并行的延迟链,设待测脉冲在第i个延迟链中从输入到被捕捉经过个延迟单元;
将每条延迟链得到的延迟单元相加,作为待测脉冲在等效延迟链中经过的延迟单元,即/>
将位置不同、各延迟单元不同的并行延迟链相互分割,最终合并得到一条等效延迟链;所述的等效延迟链经过校准阶段后将得到延迟单元信息构建为查找表存储在RAM中;在测量阶段得到待测脉冲在等效延迟链中经过个延迟单元后,查表即可得到细测量时间;
正式测量前,使用码密度法校准得到等效延迟链每个延迟单元的大小;由独立晶振产生的随机脉冲信号在延迟链上传输;对于一个校准信号上升沿,在不同延迟链上被不同延迟单元捕捉,将落在每条延迟链上的延迟单元数相加作为落在等效延迟链上的延迟单元数;经过大量随机脉冲信号校准后,通过统计方法得到等效延迟链每个延迟单元的大小,即,其中n为校准脉冲信号的次数,k为落在每个延迟单元bin上的次数;由此便可得到等效延迟链每一个延迟单元的大小,将每个延迟单元对应的延迟大小作为查找表存放在RAM中,完成对等效延迟链的校准。
2.根据权利要求1所述的多链并行分割高精度FPGA时间数字转换方法,其特征在于,所述延迟链由FPGA的专用进位单元级联构成。
3.根据权利要求2所述的多链并行分割高精度FPGA时间数字转换方法,其特征在于,将相邻的专用进位单元cin与cout相连,由此级联形成所述的延迟链,且每个专用进位单元cout作为D触发器的输入,当系统时钟上升沿到来时,D触发器对cout数据完成采样。
CN202111559443.0A 2021-12-20 2021-12-20 一种多链并行分割高精度fpga时间数字转换方法 Active CN114326358B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111559443.0A CN114326358B (zh) 2021-12-20 2021-12-20 一种多链并行分割高精度fpga时间数字转换方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111559443.0A CN114326358B (zh) 2021-12-20 2021-12-20 一种多链并行分割高精度fpga时间数字转换方法

Publications (2)

Publication Number Publication Date
CN114326358A CN114326358A (zh) 2022-04-12
CN114326358B true CN114326358B (zh) 2024-05-17

Family

ID=81052587

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111559443.0A Active CN114326358B (zh) 2021-12-20 2021-12-20 一种多链并行分割高精度fpga时间数字转换方法

Country Status (1)

Country Link
CN (1) CN114326358B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117368568A (zh) * 2023-12-07 2024-01-09 国仪量子技术(合肥)股份有限公司 一种信号频率的测量方法、装置和设备

Citations (21)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101467067A (zh) * 2006-06-15 2009-06-24 皇家飞利浦电子股份有限公司 用于飞行时间pet的集成多通道时间-数字转换器
CN101520640A (zh) * 2008-11-08 2009-09-02 中国工程物理研究院流体物理研究所 基于fpga的时间间隔测量仪
CN201331680Y (zh) * 2008-11-08 2009-10-21 中国工程物理研究院流体物理研究所 基于fpga的时间间隔测量仪
RU2455672C1 (ru) * 2011-02-14 2012-07-10 Российская Федерация, от имени которой выступает Государственная корпорация по атомной энергии "Росатом" - Госкорпорация "Росатом" Многоканальный измеритель временных интервалов
CN103676622A (zh) * 2013-10-28 2014-03-26 中国电子科技集团公司第四十一研究所 一种高精度的正负时间间隔测量方法及装置
CN104216279A (zh) * 2014-09-23 2014-12-17 西安宏泰时频技术有限公司 一种基于fpga的时间间隔测量装置
CN104597748A (zh) * 2015-02-12 2015-05-06 中国科学技术大学 一种基于fpga的时间数字变换器
CN104614976A (zh) * 2015-02-12 2015-05-13 中国科学技术大学 一种基于fpga的时间数字变换器
CN105867102A (zh) * 2009-03-30 2016-08-17 高通股份有限公司 具有改进的分辨率的时间-数字转换器(tdc)
WO2016127357A1 (zh) * 2015-02-12 2016-08-18 中国科学技术大学 一种基于fpga的时间数字变换器
CN106019923A (zh) * 2016-05-18 2016-10-12 中国科学技术大学 一种基于fpga的时间数字变换器
CN108401445A (zh) * 2017-06-30 2018-08-14 深圳市大疆创新科技有限公司 用于测量时间的电路、方法及相关芯片、系统和设备
CN109286463A (zh) * 2018-12-05 2019-01-29 北京中创为南京量子通信技术有限公司 一种基于fpga的高精度时间测量方法
CN109407500A (zh) * 2018-11-22 2019-03-01 深圳天眼激光科技有限公司 一种基于fpga的时间间隔测量方法
CN110262209A (zh) * 2019-06-03 2019-09-20 中国科学技术大学 基于fpga的时间数字变换器
JP2019176357A (ja) * 2018-03-29 2019-10-10 株式会社メガチップス 時間デジタル変換回路及び時間デジタル変換方法
CN110673463A (zh) * 2018-07-02 2020-01-10 陈昊昌 一种高线性度多通道抽头延时线时间数字转换器
CN111033312A (zh) * 2017-08-31 2020-04-17 深圳市大疆创新科技有限公司 光学距离测量设备的延迟时间校准及相关联的系统和方法
CN111723539A (zh) * 2020-05-14 2020-09-29 天津大学 一种双测时模式tdc芯片设计制造
CN112506031A (zh) * 2020-11-30 2021-03-16 中国计量科学研究院 一种激光干涉条纹信号的高精度时间间隔测量系统
CN113253597A (zh) * 2021-05-11 2021-08-13 深圳力策科技有限公司 时间数字转换装置和光学测距传感器

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE102005020903B3 (de) * 2005-05-07 2006-11-09 Infineon Technologies Ag Steuerbare Verzögerungseinrichtung
JP4666409B2 (ja) * 2006-02-17 2011-04-06 ヴェリジー(シンガポール) プライベート リミテッド 較正パルス注入による時間−デジタル変換
KR101503732B1 (ko) * 2013-06-14 2015-03-20 연세대학교 산학협력단 시간-디지털 변환기
US11664813B2 (en) * 2019-09-30 2023-05-30 Seiko Epson Corporation Delay circuit, time to digital converter, and A/D conversion circuit

Patent Citations (21)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101467067A (zh) * 2006-06-15 2009-06-24 皇家飞利浦电子股份有限公司 用于飞行时间pet的集成多通道时间-数字转换器
CN101520640A (zh) * 2008-11-08 2009-09-02 中国工程物理研究院流体物理研究所 基于fpga的时间间隔测量仪
CN201331680Y (zh) * 2008-11-08 2009-10-21 中国工程物理研究院流体物理研究所 基于fpga的时间间隔测量仪
CN105867102A (zh) * 2009-03-30 2016-08-17 高通股份有限公司 具有改进的分辨率的时间-数字转换器(tdc)
RU2455672C1 (ru) * 2011-02-14 2012-07-10 Российская Федерация, от имени которой выступает Государственная корпорация по атомной энергии "Росатом" - Госкорпорация "Росатом" Многоканальный измеритель временных интервалов
CN103676622A (zh) * 2013-10-28 2014-03-26 中国电子科技集团公司第四十一研究所 一种高精度的正负时间间隔测量方法及装置
CN104216279A (zh) * 2014-09-23 2014-12-17 西安宏泰时频技术有限公司 一种基于fpga的时间间隔测量装置
CN104597748A (zh) * 2015-02-12 2015-05-06 中国科学技术大学 一种基于fpga的时间数字变换器
CN104614976A (zh) * 2015-02-12 2015-05-13 中国科学技术大学 一种基于fpga的时间数字变换器
WO2016127357A1 (zh) * 2015-02-12 2016-08-18 中国科学技术大学 一种基于fpga的时间数字变换器
CN106019923A (zh) * 2016-05-18 2016-10-12 中国科学技术大学 一种基于fpga的时间数字变换器
CN108401445A (zh) * 2017-06-30 2018-08-14 深圳市大疆创新科技有限公司 用于测量时间的电路、方法及相关芯片、系统和设备
CN111033312A (zh) * 2017-08-31 2020-04-17 深圳市大疆创新科技有限公司 光学距离测量设备的延迟时间校准及相关联的系统和方法
JP2019176357A (ja) * 2018-03-29 2019-10-10 株式会社メガチップス 時間デジタル変換回路及び時間デジタル変換方法
CN110673463A (zh) * 2018-07-02 2020-01-10 陈昊昌 一种高线性度多通道抽头延时线时间数字转换器
CN109407500A (zh) * 2018-11-22 2019-03-01 深圳天眼激光科技有限公司 一种基于fpga的时间间隔测量方法
CN109286463A (zh) * 2018-12-05 2019-01-29 北京中创为南京量子通信技术有限公司 一种基于fpga的高精度时间测量方法
CN110262209A (zh) * 2019-06-03 2019-09-20 中国科学技术大学 基于fpga的时间数字变换器
CN111723539A (zh) * 2020-05-14 2020-09-29 天津大学 一种双测时模式tdc芯片设计制造
CN112506031A (zh) * 2020-11-30 2021-03-16 中国计量科学研究院 一种激光干涉条纹信号的高精度时间间隔测量系统
CN113253597A (zh) * 2021-05-11 2021-08-13 深圳力策科技有限公司 时间数字转换装置和光学测距传感器

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
基于FPGA的高精度多通道时间数字转换器设计;王巍;董永孟;李捷;熊拼搏;周浩;杨正琳;王冠宇;袁军;周玉涛;;微电子学;20151220(第06期);第698-705页 *

Also Published As

Publication number Publication date
CN114326358A (zh) 2022-04-12

Similar Documents

Publication Publication Date Title
US5166959A (en) Picosecond event timer
CN107643674B (zh) 一种基于FPGA进位链的Vernier型TDC电路
Kwiatkowski et al. Efficient implementation of multiple time coding lines-based TDC in an FPGA device
CN109143832B (zh) 一种高精度多通道的时间数字转换器
WO2016127357A1 (zh) 一种基于fpga的时间数字变换器
CN109634089B (zh) 一种应用于工艺非受控检测的两级tdc电路
CN108061848B (zh) 基于fpga的加法进位链延时的测量方法及系统
US11782393B2 (en) Time to digital conversion
CN114326358B (zh) 一种多链并行分割高精度fpga时间数字转换方法
CN110515292B (zh) 基于双向运行环形进位链的tdc电路及测量方法
CN112486008B (zh) 基于tdc的分辨率可调时间测量统计系统及方法
CN113253597A (zh) 时间数字转换装置和光学测距传感器
CN110703583A (zh) 基于soc的多通道高精度大量程时间数字转换器
CN106681127A (zh) 移位寄存器电路、相位差计算方法及时间数字转换器
KR101503732B1 (ko) 시간-디지털 변환기
CN112578661A (zh) 一种用于fpga型时间数字转换器的延迟线校准电路
Aloisio et al. FPGA implementation of a high-resolution time-to-digital converter
CN113917830B (zh) 循环游标延时链电路、时间数字转换器和信号选取方法
Sesta et al. A novel sub-10 ps resolution TDC for CMOS SPAD array
CN114637182B (zh) 基于fpga进位链的tdc细时间测量系统及方法
CN103675383B (zh) 一种量测波形的电路
CN114967411B (zh) 一种具备自动复位机制的多级时间数字转换器
CN109143833B (zh) 一种应用于高分辨率时间数字转换器的小数部分测量电路
CN215769392U (zh) 一种时间数字转换装置和光学测距传感器
CN116318140A (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