CN107192990B - 外推法测雷达散射截面积 - Google Patents

外推法测雷达散射截面积 Download PDF

Info

Publication number
CN107192990B
CN107192990B CN201710470620.5A CN201710470620A CN107192990B CN 107192990 B CN107192990 B CN 107192990B CN 201710470620 A CN201710470620 A CN 201710470620A CN 107192990 B CN107192990 B CN 107192990B
Authority
CN
China
Prior art keywords
section
measured target
target body
radar cross
emission source
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
CN201710470620.5A
Other languages
English (en)
Other versions
CN107192990A (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.)
National Institute of Metrology
Original Assignee
National Institute of Metrology
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 National Institute of Metrology filed Critical National Institute of Metrology
Priority to CN201710470620.5A priority Critical patent/CN107192990B/zh
Publication of CN107192990A publication Critical patent/CN107192990A/zh
Application granted granted Critical
Publication of CN107192990B publication Critical patent/CN107192990B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/40Means for monitoring or calibrating
    • G01S7/4004Means for monitoring or calibrating of parts of a radar system

Landscapes

  • Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

本发明公开了一种外推法测雷达散射截面积的方法,包括:对发射源与被测目标体的位置进行对准,以使发射源正对被测目标体的被测位置;在发射源远离被测目标体的过程中,对应获取在多个测量位置下发射源的辐射功率或电磁波幅度,以及接收到的目标散射的空间辐射功率或电磁波幅度;根据多个测量位置下发射源的辐射功率或电磁波幅度,以及接收到的目标散射的空间辐射功率或电磁波幅度得到发射源和被测目标体随距离变化的雷达散射截面积测量曲线;根据雷达散射截面积测量曲线得到被测目标体的雷达散射截面积随距离变化的函数;对函数求距离无穷远极限,得到被测目标体的雷达散射截面积。本发明具有如下优点:测量准确度高、适用范围广。

Description

外推法测雷达散射截面积
技术领域
本发明涉及雷达技术领域,具体涉及一种外推法测雷达散射截面积。
背景技术
当前的雷达散射截面积(RCS)测量方法是通过球体对整个测试系统进行校准(或称“标定”,“定标”),通过先后将球体和被测目标在测试场中进行RCS测量值的比较,得到被测目标的RCS值,球体的RCS值是由理论计算得到(金属球:Pi*r^2,r为球体半径)。用球体标定RCS测量系统的缺点是,由于球体在所有观测方向上都是球体,360度转动RCS没有变化,从而缺少对RCS测量系统的方位角、俯仰角的校准信息。此外,球体对所接收到的电磁波反射小,从而使测量信号的信号噪声比(SNR)和信号杂波比(SCR)低,会造成标定后的RCS测量系统的测量灵敏度低,测量结果不确定度大。
如果采用带有角度RCS信息的其他反射体对RCS测量系统进行校准,必须首先知道该反射体的RCS量值,目前不能得到准确的复杂几何体的RCS量值。
发明内容
本发明旨在至少解决上述技术问题之一。
为此,本发明的目的在于提出一种外推法测雷达散射截面积的方法,可以提升RCS测量的准确性。
为了实现上述目的,本发明的实施例公开了一种外推法测雷达散射截面积的方法,包括以下步骤:S1:对发射源与被测目标体的位置进行对准,以使所述发射源正对所述被测目标体的被测位置,其中,所述发射源和所述被测目标体均设置在微波暗室中;S2:在所述发射源远离所述被测目标体的过程中,对应获取在多个测量位置下所述发射源的辐射功率或电磁波幅度,以及接收到的目标散射的空间辐射功率或电磁波幅度;S3:根据所述多个测量位置下所述发射源的辐射功率或电磁波幅度,以及接收到的目标散射的空间辐射功率或电磁波幅度,得到所述发射源和所述被测目标体随距离变化的雷达散射截面积测量曲线;S4:根据所述雷达散射截面积测量曲线得到所述被测目标体的雷达散射截面积随距离变化的函数;S5:对所述雷达散射截面积随距离变化的函数求距离无穷远极限,得到所述被测目标体的雷达散射截面积。
进一步地,步骤S2进一步包括:对所述每个测量位置进行多次测量,对测量数据进行相参积累得到所述每个测量位置的所述发射源的辐射功率或电磁波幅度,以及接收到的目标散射的空间辐射功率或电磁波幅度。
进一步地,步骤S4进一步包括:对所述雷达散射截面积测量曲线进行平滑处理,并采用多项式拟合的方式得到被测目标体的雷达散射截面积随距离变化的函数。
进一步地,当所述被测目标体为非球面体时,还包括:改变所述发射源和所述被测目标体之间的夹角并测量被测目标体的雷达散射截面积,以得到所述发射源和所述被测目标体之间的在不同夹角下所述被测目标体的雷达散射截面积。
进一步地,使用矢量网络分析仪将激光对准装置连接发射源的波导口面,并使用所述激光对准装置检测所述发射源与所述被测目标体之间的是否对准。
根据本发明实施例的外推法测雷达散射截面积的方法与现有技术相比具有如下优点:
测量准确,可以使目前国内RCS测量准确度提升5倍以上,本发明第二个优点是测量结果信息全面,可以准确测量得到被测目标在不同方位角,不同俯仰角下的RCS量值;
可以提升RCS测量准确度,具有更小的测量不确定度,更丰富的RCS测量信息,为隐身飞机、飞行器等国防领域及民用领域的飞机、飞行器的隐身效果、目标识别等方面提供精准测量;
可以通过对自行设计的RCS目标体进行精确测量,将该RCS目标体作为校准件发放给全国RCS测试场用于作为标定标准器,有利于建立我国RCS计量体系,统一全国RCS测量量值,达到统一、准确、可靠、有效。
本发明的附加方面和优点将在下面的描述中部分给出,部分将从下面的描述中变得明显,或通过本发明的实践了解到。
附图说明
本发明的上述和/或附加的方面和优点从结合下面附图对实施例的描述中将变得明显和容易理解,其中:
图1是本发明实施例的外推法测雷达散射截面积的方法的流程图;
图2是本发明一个实施例的外推法测雷达散射截面积采用相关设备的结构示意图。
具体实施方式
下面详细描述本发明的实施例,所述实施例的示例在附图中示出,其中自始至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。下面通过参考附图描述的实施例是示例性的,仅用于解释本发明,而不能理解为对本发明的限制。
在本发明的描述中,需要理解的是,术语“中心”、“纵向”、“横向”、“上”、“下”、“前”、“后”、“左”、“右”、“竖直”、“水平”、“顶”、“底”、“内”、“外”等指示的方位或位置关系为基于附图所示的方位或位置关系,仅是为了便于描述本发明和简化描述,而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作,因此不能理解为对本发明的限制。此外,术语“第一”、“第二”仅用于描述目的,而不能理解为指示或暗示相对重要性。
在本发明的描述中,需要说明的是,除非另有明确的规定和限定,术语“安装”、“相连”、“连接”应做广义理解,例如,可以是固定连接,也可以是可拆卸连接,或一体地连接;可以是机械连接,也可以是电连接;可以是直接相连,也可以通过中间媒介间接相连,可以是两个元件内部的连通。对于本领域的普通技术人员而言,可以具体情况理解上述术语在本发明中的具体含义。
参照下面的描述和附图,将清楚本发明的实施例的这些和其他方面。在这些描述和附图中,具体公开了本发明的实施例中的一些特定实施方式,来表示实施本发明的实施例的原理的一些方式,但是应当理解,本发明的实施例的范围不受此限制。相反,本发明的实施例包括落入所附加权利要求书的精神和内涵范围内的所有变化、修改和等同物。
以下结合附图描述本发明。
图1是本发明实施例的外推法测雷达散射截面积的方法的流程图。如图1所示,根据本发明实施例的外推法测雷达散射截面积的方法,包括以下步骤:
S1:对发射源与被测目标体的位置进行对准,以使所述发射源正对所述被测目标体的被测位置,其中,所述发射源和所述被测目标体均设置在微波暗室中;
S2:在所述发射源远离所述被测目标体的过程中,对应获取在多个测量位置下所述发射源的辐射功率或电磁波幅度,以及接收到的目标散射的空间辐射功率或电磁波幅度;
S3:根据所述多个测量位置下所述发射源的辐射功率或电磁波幅度,以及接收到的目标散射的空间辐射功率或电磁波幅度,得到所述发射源和所述被测目标体随距离变化的雷达散射截面积测量曲线;
S4:根据所述雷达散射截面积测量曲线得到所述被测目标体的雷达散射截面积随距离变化的函数;
S5:对所述雷达散射截面积随距离变化的函数求距离无穷远极限,得到所述被测目标体的雷达散射截面积。
图2是本发明一个实施例的外推法测雷达散射截面积采用相关设备的结构示意图。如图2所示,在本发明的一个实施例中,外推法测雷达散射截面积采用的相关设备包括:导轨引导系统、导轨覆盖系统、发射天线塔和驱动系统、上极化下方位天线定位器系统、对准设备等部分。
使用外推法测量RCS要求导轨引导系统精确对齐。导轨引导系统安装在外推法暗室的坑道内,用于支撑发射天线塔,并引导发射塔沿暗室的移动轴移动,同时保持发射塔在移动过程中的对准度在指定误差范围内。导轨引导系统被直接安装到混凝土地板上。导轨引导系统与外推法暗室内的地板和发射塔连接。
为了避免发射天线塔和接收天线塔之间导轨暴露所引起的反射,暗室内具备一个自动导轨覆盖系统。在本发明的一个示例中,导轨覆盖系统采用一个可移动地板,由沿暗室纵向排列在地面上的移动托架组成。导轨覆盖系统会自动对导轨覆盖吸波材料,同时为操作人员提供进出发射塔的通路。当轨道不覆盖吸波材料时,RCS为关状态;反之为开状态。
用外推法测量RCS时,需要测量S11(反射波/发射波)随天线间距离的变化。因此设定一个发射塔能够承载发射天线和其它相关设备,以及一个驱动系统来移动发射塔。发射天线塔用于承载微调对激光对准装置、测量仪器、发射天线和极化定位器。发射天线塔包括一个距离地面大约2m的工作台。其中,激光对准装置可以是激光对准仪,也可以是激光笔。
发射天线塔由长约9米的滚珠丝杠来驱动,采用带有绝对编码的直流无刷电机。电机和编码器由运动控制器控制,运动控制器本身又可以通过计算机远程控制,也可通过发射塔上的本地控制面板控制。运动控制器能够控制发射天线塔和导轨覆盖系统。根据发射天线塔的已知位置,运动控制器对导轨覆盖系统实行全自动控制。
在本发明的一个实施例中,被测目标体为金属方板,金属方板被固定在接收天线塔的天线架上,可以通过改变方板与发射天线口面之间的夹角,从而可以得到目标体在不同角度下的RCS值。
在本发明的一个实施例中,使用矢量网络分析仪将激光对准装置连接发射源的波导口面,并使用激光对准装置检测所述发射源与所述被测目标体之间的是否对准。
准备工作完成后,将微波暗室门关闭,实验操作人员在天线塔上开始进行实验操作,然后设定测量参数。参数设置时包括对时域门进行设置,用于滤除目标体周围物体的散射信号,尤其是目标体背后天线塔的散射信号。标准增益喇叭天线口面与标准散射体之间的距离为R,测量R取不同值时的S11
当待测体为金属球时,将金属球放置在泡沫支架上,将喇叭天线固定在发射塔上,并通过调节泡沫支架下方的升降台使喇叭天线口面中心与金属球球心对齐。喇叭天线口面到金属球球心的距离为R。
当待测体为铝合金方板时,将铝合金方板通过支架固定在天线架台上,将喇叭天线固定在发射塔上,喇叭天线口面到铝合金方板表面的距离为R。发射塔沿着导轨移动改变R的大小,并通过网络分析仪得到不同距离R下的S11值。
在测量过程中,每移动一定的距离,需要根据当前天线-目标之间的距离对时域门的center重新进行设置,从而得到更为精确的测量值。
在本发明的一个实施例中,在步骤S2中,对每个测量位置进行多次测量,对测量数据进行相参积累得到发射源的辐射功率或电磁波幅度,以及接收到的目标散射的空间辐射功率或电磁波幅度,从而提高信噪比和信杂比。
实验测得数据是不同距离下的S11值,将此值带入到以下公式中可以得到对应的RCS值:
其中,R为喇叭天线与被测目标体之间的距离,G为天线增益,λ为波长。
通过对RCS测量结果分析,发现RCS测量值随着距离的增大是波动上升的,并且最后上升趋势变缓,逐渐趋于定义值(无穷远RCS量值)。
在本发明的一个实施例中,步骤S4进一步包括:对所述雷达散射截面积测量曲线进行平滑处理,并采用多项式拟合的方式得到被测目标体的雷达散射截面积随距离变化的函数。
具体地,空间驻波形成的纹波可以通过滑动平均的方法滤掉。对于平滑后得到的外推曲线,采用多项式拟合的方式得到RCS随距离变化的函数f(R),对f(R)求取距离为无穷大时的极限值,得到被测目标的RCSDUT量值,即RCSDUT=limE→∞(f(R))。
本发明实施例的外推法测雷达散射截面积的方法与现有技术相比具有如下优点:
测量准确,可以使目前国内RCS测量准确度提升5倍以上,本发明第二个优点是测量结果信息全面,可以准确测量得到被测目标在不同方位角,不同俯仰角下的RCS量值;
可以提升我国RCS测量准确度,具有更小的测量不确定度,更丰富的RCS测量信息,为隐身飞机、飞行器等国防领域及民用领域的飞机、飞行器的隐身效果、目标识别等方面提供精准测量;
可以通过对自行设计的RCS目标体进行精确测量,将该RCS目标体作为校准件发放给全国RCS测试场用于作为标定标准器,有利于建立我国RCS计量体系,统一全国RCS测量量值,达到统一、准确、可靠、有效。
另外,本发明实施例的外推法测雷达散射截面积的方法的其它构成以及作用对于本领域的技术人员而言都是已知的,为了减少冗余,不做赘述。
在本说明书的描述中,参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本发明的至少一个实施例或示例中。在本说明书中,对上述术语的示意性表述不一定指的是相同的实施例或示例。而且,描述的具体特征、结构、材料或者特点可以在任何的一个或多个实施例或示例中以合适的方式结合。
尽管已经示出和描述了本发明的实施例,本领域的普通技术人员可以理解:在不脱离本发明的原理和宗旨的情况下可以对这些实施例进行多种变化、修改、替换和变型,本发明的范围由权利要求及其等同限定。

Claims (5)

1.一种外推法测雷达散射截面积的方法,其特征在于,包括以下步骤:
S1:对发射源与被测目标体的位置进行对准,以使所述发射源正对所述被测目标体的被测位置,其中,所述发射源和所述被测目标体均设置在微波暗室中;
S2:在所述发射源远离所述被测目标体的过程中,对应获取在多个测量位置下所述发射源的辐射功率或电磁波幅度,以及接收到的目标散射的空间辐射功率或电磁波幅度;
S3:根据所述多个测量位置下所述发射源的辐射功率或电磁波幅度,以及接收到的目标散射的空间辐射功率或电磁波幅度,得到所述被测目标体的雷达散射截面积随所述发射源和所述被测目标体之间的距离的变化而变化的雷达散射截面积测量曲线;
S4:根据所述雷达散射截面积测量曲线得到所述被测目标体的雷达散射截面积随距离变化的函数;
S5:对所述雷达散射截面积随距离变化的函数求距离无穷远极限,得到所述被测目标体的雷达散射截面积。
2.根据权利要求1所述的方法,其特征在于,步骤S2进一步包括:
对所述每个测量位置进行多次测量,对测量数据进行相参积累得到所述每个测量位置的所述发射源的辐射功率或电磁波幅度,以及接收到的目标散射的空间辐射功率或电磁波幅度。
3.根据权利要求1所述的方法,其特征在于,步骤S4进一步包括:
对所述雷达散射截面积测量曲线进行平滑处理,并采用多项式拟合的方式得到被测目标体的雷达散射截面积随距离变化的函数。
4.根据权利要求1所述的方法,其特征在于,当所述被测目标体为非球面体时,还包括:
改变所述发射源和所述被测目标体之间的夹角并测量被测目标体的雷达散射截面积,以得到所述发射源和所述被测目标体之间不同的夹角对应的所述被测目标体的雷达散射截面积。
5.根据权利要求1或4所述的方法,其特征在于,使用矢量网络分析仪将激光对准装置连接发射源的波导口面,并使用所述激光对准装置检测所述发射源与所述被测目标体之间是否对准。
CN201710470620.5A 2017-06-20 2017-06-20 外推法测雷达散射截面积 Active CN107192990B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710470620.5A CN107192990B (zh) 2017-06-20 2017-06-20 外推法测雷达散射截面积

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710470620.5A CN107192990B (zh) 2017-06-20 2017-06-20 外推法测雷达散射截面积

Publications (2)

Publication Number Publication Date
CN107192990A CN107192990A (zh) 2017-09-22
CN107192990B true CN107192990B (zh) 2019-09-17

Family

ID=59878622

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710470620.5A Active CN107192990B (zh) 2017-06-20 2017-06-20 外推法测雷达散射截面积

Country Status (1)

Country Link
CN (1) CN107192990B (zh)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109975778A (zh) * 2019-03-25 2019-07-05 中国计量科学研究院 一种天线雷达散射截面测量方法、装置及计算机设备
CN110261837B (zh) * 2019-06-27 2022-10-28 中国航空工业集团公司雷华电子技术研究所 一种基于航迹信息的复杂目标rcs计算方法
WO2021097636A1 (zh) * 2019-11-19 2021-05-27 深圳市大疆创新科技有限公司 标物雷达散射截面积确定方法、装置和存储介质
CN116047176B (zh) * 2022-12-05 2023-12-19 北京信凯达科技有限公司 一种带有自动避让装置的暗室电磁检测系统及其检测方法
CN116449327B (zh) * 2023-04-25 2023-10-13 中国计量科学研究院 一种基于脉冲压缩方法的外推法测量雷达散射截面积的方法和系统

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH1152044A (ja) * 1997-08-07 1999-02-26 Tech Res & Dev Inst Of Japan Def Agency レーダ反射断面積測定装置
CN102928829A (zh) * 2012-12-04 2013-02-13 上海无线电设备研究所 空间目标参数反演方法
CN102944872A (zh) * 2012-11-23 2013-02-27 北京航空航天大学 雷达散射截面近场到近场的变换方法
CN103777185A (zh) * 2014-01-15 2014-05-07 北京环境特性研究所 基于暗室扫频获取目标体雷达散射截面的方法及装置
CN105572652A (zh) * 2016-03-18 2016-05-11 西北工业大学 一种利用外推获得具有多次散射目标远场rcs的方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP5424588B2 (ja) * 2008-07-11 2014-02-26 三菱電機株式会社 レーダ断面積の測定装置及び方法並びにレーダ断面積の測定プログラム

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH1152044A (ja) * 1997-08-07 1999-02-26 Tech Res & Dev Inst Of Japan Def Agency レーダ反射断面積測定装置
CN102944872A (zh) * 2012-11-23 2013-02-27 北京航空航天大学 雷达散射截面近场到近场的变换方法
CN102928829A (zh) * 2012-12-04 2013-02-13 上海无线电设备研究所 空间目标参数反演方法
CN103777185A (zh) * 2014-01-15 2014-05-07 北京环境特性研究所 基于暗室扫频获取目标体雷达散射截面的方法及装置
CN105572652A (zh) * 2016-03-18 2016-05-11 西北工业大学 一种利用外推获得具有多次散射目标远场rcs的方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
"基于扫频时域法测量的RCS外推技术研究";李南京等;《微波学报》;20070831;第23卷(第4期);全文
"微波暗室目标RCS测试方法的研究与试验";郭静;《中国优秀硕士学位论文全文数据库 信息科技辑》;20090615;全文
"近场测量不同距离不同方位RCS 的一种方法";童广德等;《制导与引信》;20090630;第30卷(第2期);全文
"雷达目标RCS近远场变换";朱胜楠;《中国优秀硕士学位论文全文数据库 信息科技辑》;20110415;全文

Also Published As

Publication number Publication date
CN107192990A (zh) 2017-09-22

Similar Documents

Publication Publication Date Title
CN107192990B (zh) 外推法测雷达散射截面积
CN105548729B (zh) 一种阵列天线辐射特性的快速测量方法
CN103630761A (zh) 一种多探头球面近场通道校准装置及方法
CA2719193C (en) Method and device for measuring a radiation field
US8098194B2 (en) Determining characteristics of a radar cross section (RCS) test range
CN105866759B (zh) 用于目标rcs测量中多重定标和背景提取的装置设计及其信号处理方法
CN206369770U (zh) 雷达罩测试系统
CN105352978A (zh) 手持式吸波材料反射率测量装置
CN108562801A (zh) 一种阵列天线测试系统及其测试方法
CN115993584B (zh) 一种大俯仰角度雷达散射截面数据测量系统及其测量方法
US20120154122A1 (en) Method and apparatus of radio frequency testing
CN111025032A (zh) 一种基于升空平台的天线波束测量系统及方法
CN106019247B (zh) 地物波谱rcs测量系统运载平台圆周运动方位误差修正方法
CN109374990A (zh) 一种天线相位中心校准方法
CN209784528U (zh) 激光雷达及其激光垂直校准装置
CN106405256A (zh) 一种平面近场暗室散射的测试与补偿方法
CN208818765U (zh) 一种天线测试系统
US7498977B2 (en) Field probe form the angular response of a rigid body
Cote Automated swept-angle bistatic scattering measurements using continuous wave radar
CN114114171A (zh) 一种多功能内场散射成像测量系统、方法及应用
CN109782331A (zh) 一种束流位置探测器电中心标定装置及标定方法
CN212514997U (zh) 一种用于毫米波雷达的自动测试系统
CN108981922B (zh) 一种微波黑体发射率测量装置及测量方法
CN108489396B (zh) 一种二维转顶精度检测方法
KR100651741B1 (ko) 태그와 리더 사이의 인식거리 측정장치

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