CN105891851B - 一种基于导航卫星漏信号进行定位在轨试验验证方法 - Google Patents

一种基于导航卫星漏信号进行定位在轨试验验证方法 Download PDF

Info

Publication number
CN105891851B
CN105891851B CN201510036026.6A CN201510036026A CN105891851B CN 105891851 B CN105891851 B CN 105891851B CN 201510036026 A CN201510036026 A CN 201510036026A CN 105891851 B CN105891851 B CN 105891851B
Authority
CN
China
Prior art keywords
earth
detector
navigation
satellite
aeronautical
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
CN201510036026.6A
Other languages
English (en)
Other versions
CN105891851A (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.)
Beijing Institute of Spacecraft System Engineering
Original Assignee
Beijing Institute of Spacecraft System Engineering
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 Beijing Institute of Spacecraft System Engineering filed Critical Beijing Institute of Spacecraft System Engineering
Priority to CN201510036026.6A priority Critical patent/CN105891851B/zh
Publication of CN105891851A publication Critical patent/CN105891851A/zh
Application granted granted Critical
Publication of CN105891851B publication Critical patent/CN105891851B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Position Fixing By Use Of Radio Waves (AREA)

Abstract

本发明方法提出一种基于导航卫星漏信号进行定位的在轨试验验证方法,首先确定在轨测试飞行高度,计算探测器、导航卫星和地球的相对位置关系,根据三者位置关系计算任意位置相对导航卫星的地球矢量与探测器矢量的夹角,判断导航卫星的可见性;根据导航卫星可见性判定结果,确定星载导航接收机开机工作的高度范围和开机时间,根据导航可见性的判定结果,分析在轨飞行姿态,确定GNSS导航接收设备接收天线的波束角范围;在轨测试记录导航漏信号的接收情况及导航定位结果;比较分析GNSS接收机导航漏信号的接收能力及其对导航性能的影响,本发明方法解决了深空探测和高轨道地球卫星的导航卫星漏信号接收设备的全方位功能性能测试和利用导航卫星漏信号进行定位的系统设计问题。

Description

一种基于导航卫星漏信号进行定位在轨试验验证方法
技术领域
本发明属于深空探测技术领域,具体涉及一种基于导航卫星漏信号的导航定位在轨试验验证方法。
背景技术
空间导航卫星漏信号即地球对面的导航卫星信号和导航卫星的旁瓣信号,也称作弱导航信号,通常信号强度在-180dBW~-188dBW的范围内。随着探测器距离地球距离增加,接收设备接收到的导航信号呈指数级衰减,对接收设备的性能要求越苛刻。
目前常用的GNSS导航定位试验方法主要通过地面模拟试验验证设备性能,或在轨利用静止轨道卫星或其他高轨道卫星进行利用高增益天线验证弱导航信号的接收情况,尚未在深空探测任务中使用。地面试验方法是根据安装导航接收设备的卫星轨道数据,利用导航卫星信号模拟源来验证设备性能,高轨卫星在轨利用GNSS导航信号接收设备接收GNSS导航卫星信号,验证接收设备对弱导航信号的接收能力,根据测试结果验证的设备的可用性及在轨性能,上述两种GNSS导航信号接收设备的性能测试方法主要存在如下缺点:
(1)利用地面模拟源不能完全准确模拟设备在轨工作时导航信号的状态,不能准确反映设备的在轨工作性能;
(2)设备在地面模拟试验过程中的工作环境与在轨飞行过程中的空间环境有明显差别,地面试验不能充分验证设备对空间环境的适应性;
(3)利用静止轨道卫星进行飞行试验时,GNSS导航接收设备飞行高度和飞行速度为固定值,不能充分考核GNSS导航接收设备在不同飞行高度和飞行速度下的工作性能,验证不充分。
发明内容
本发明解决的技术问题是提出了一种导航卫星漏信号接收设备在轨试验验证方法,以克服已有技术的缺陷,解决深空探测和高轨道地球卫星的导航卫星漏信号接收设备的全方位功能性能测试和利用导航卫星漏信号进行定位的系统设计问题。
本发明提出一种基于导航卫星漏信号进行定位在轨试验验证方法,包括以下步骤:
1)确定在轨测试飞行高度,通过对GNSS系统的覆盖特性分析,高于地面3000km的范围内为弱导航信号的覆盖区域,确定GNSS导航卫星接收设备在3000km以上进行在轨测试;
2)根据探测器飞行轨道动力学,获取地球、任意导航卫星和探测器在地心惯性系Oi-XiYiZi下的位置,分别计算地球→导航卫星、地球→探测器、探测器→导航卫星的距离关系;根据所述的距离关系得到探测器的不同位置,通过计算探测器与导航卫星连线和导航卫星与地球连线的夹角,判断出导航卫星的可见性;
3)根据探测器运行轨道,当探测器距离高度大于3000km时,每隔500km对探测器和所有GNSS导航卫星的可见性进行取点分析,当有4颗导航卫星满足导航卫星可见性要求时,将对应的飞行高度确定为GNSS导航接收设备在轨测试开始时刻;
4)对探测器在轨飞行姿态进行分析,通过计算导航卫星接收天线轴向与探测器-地心连线矢量的夹角以及探测器对GNSS导航卫星的半张角,确定GNSS导航接收天线的波束角范围;
5)GNSS导航接收设备在轨测试过程中,记录地球到其他天体转移过程中正常接收导航信号GNSS导航卫星小于4颗时的飞行高度,该飞行高度作为返回地球过程中GNSS导航接收设备开始在轨测试的飞行高度,从而确定返回地球过程中GNSS导航接收设备的开机工作时间;
6)记录地球到其他天体转移段探测器远离地球的过程中接收导航卫星漏信号的情况和其他天体到地球返回段探测器接近地球的过程中接收到的由弱渐强的导航漏信号的情况,比较在不同高度范围内,对于可见导航卫星信号的捕获的情况,从而对导航接收机信号捕获性能进行评估,验证分析GNSS导航接收设备信号接收的性能;
7)根据地球到其他天体转移和其他天体到地球转移过程中的GNSS导航接收设备的在轨数据,完成探测器在轨飞行过程中轨道确定,同时根据在轨飞行结果,完成GNSS导航接收设备的性能评估,确认GNSS导航设备的在轨适用范围。
本方法充分利用自由返回任务飞行轨道的特点,在不同的飞行高度和飞行速度范围内完成高轨GNSS导航接收设备的在轨飞行试验,能够充分考核GNSS导航接收设备的性能和任务适用性,与以往试验方法相比,具有如下优点:
(1)本方法利用实际的在轨飞行完成高轨GNSS导航接收设备的试验验证,能够真实模拟GNSS导航卫星的弱导航信号状态,更加准确的反映高轨导航接收设备的性能。
(2)高轨导航接收设备在轨飞行过程中经历空间辐射环境的考核,能够验证高轨导航接收设备对空间环境的适应性。
(3)利用深空探测任务进行飞行试验,能够考核高轨GNSS导航接收设备在数千公里至数万公里(最远可达40万公里)的由远及近和由近及远飞行过程中的性能,飞行速度最大可达到11km/s,不仅能够验证高轨GNSS导航接收设备对地球静止轨道卫星的适用性,同时能够充分考核对深空探测任务的适用性,在很大程度上拓展了高轨GNSS导航接收设备的工程适用范围。
附图说明
图1为基于导航卫星漏信号进行定位的在轨试验验证流程图;
图2为探测器对导航卫星可见性分析示意图;
图3为GNSS导航接收天线波速角分析示意图;
图4为月地高速再入返回飞行器飞行高度变化情况示意图;
图5为月地高速再入返回飞行器飞行速度变化情况示意图。
具体实施方式
下面结合附图及具体实施方式对本发明进行详细说明。
高轨卫星及深空探测器的全球卫星导航定位系统(GNSS)导航信号接收设备通过接收来自地球对面导航卫星信号和导航卫星的旁瓣信号等漏信号实现定位,能够解决地面测控站覆盖不足,导致测定轨精度不能满足飞行任务要求的问题,为高轨道卫星和深空探测器提供新的测定轨手段。利用自由返回轨道的特点,能够在较大的高度和速度范围内充分验证导航卫星信号接收设备在距离由远及近及由近及远过程中的性能,明确设备的适用范围。本方法适用于月球探测及其它深空探测任务中利用GNSS导航卫星信号实现导航定位的在轨试验和方案设计。
本发明方法是对基于深空探测任务GNSS导航接收系统进行系统设计的一种方法,要点是充分利用深空探测任务的特点,充分考虑地球出发转移和返回过程中高度和速度的变化,设计在轨工作的时间,并考虑飞行过程中探测器的飞行姿态,对GNSS导航接收天线的波束角进行设计,实施过程如图1所示,具体如下:
步骤一、对GNSS系统的覆盖特性进行分析,GNSS系统能够覆盖地面至距离地面3000km高度范围内的导航定位要求,在高于地面3000km的范围内为弱导航信号的覆盖区域,因此GNSS导航接收设备应该在3000km以上进行在轨测试,充分考核设备接收和处理弱导航信号的能力。
步骤二、根据探测器飞行轨道动力学,根据探测器飞行轨道动力学,获取地球、任意导航卫星和探测器在地心J2000惯性坐标系Oi-XiYiZi下位置的直角坐标:(Xei,Yei,Zei),(XGi,YGi,ZGi),(XHi,YHi,ZHi)。分别计算地球→导航卫星、地球→探测器、探测器→导航卫星的距离关系ReG、ReH、RHG
根据式①计算探测器在不同位置时,探测器与导航卫星连线和导航卫星与地球连线的夹角α,即计算导航卫星的可见性。其中:
已知a=Re+a/ReG,Re+a=地球半径+大气层高度,b为导航卫星天线波束最大幅角(含旁瓣),均可视为常值。则探测器对导航卫星可见性判据为:
a≤α≤b ③
当探测器满足公式③规定的条件时,则对该导航卫星可见,否则为不可见。如图2所示。
步骤三、根据探测器运行轨道,当探测器距离高度大于3000km时,每隔500km对探测器和所有GNSS导航卫星的可见性进行取点分析,当有4颗导航卫星满足公式③的要求时,将对应的飞行高度确定为GNSS导航接收设备在轨测试开始时刻。
步骤四、对探测器在轨飞行姿态进行分析,确定GNSS导航接收天线的波束角范围。如图3所示,S为导航接收天线轴向,ReH为探测器与地心连线距离,R为GNSS导航卫星轨道半径,卫星导航接收天线轴向与探测器-地心连线矢量ReH的夹角为σ,γ为探测器对GNSS导航卫星的半张角。
当R<ReH时,
γ=arcsin(R/ReH) ④
当R>ReH时,
γ=180° ⑤
GNSS导航卫星接收天线的半波束角应在满足接收要求的情况下尽可能小(波束角越小对弱信号的接收能力越强),GNSS导航卫星接收天线的半波束角β为:
当γ+σ<180°时,
β=γ+σ ⑥
当γ+σ>180°时,
β=180° ⑦
步骤五、GNSS导航接收设备在轨测试过程中,记录地球-其他天体转移过程(如地球到月球的转移过程)中能够正常接收导航信号GNSS导航卫星小于4颗时的飞行高度h,作为返回地球过程中GNSS导航接收设备开始在轨测试的飞行高度,从而确定返回地球过程中GNSS导航接收设备的开机工作时间。
步骤六、记录地球到其他天体(如地球-月球)转移段探测器远离地球的过程中接收导航卫星漏信号的情况和其他天体到地球(如月球-地球)返回段探测器接近地球的过程中接收到的由弱渐强的导航漏信号的情况,包括可见相对距离、可见导航卫星数量、旁瓣信号范围等,比较在不同高度范围内,对于可见导航卫星信号的捕获的情况,从而对导航接收机信号捕获性能进行评估,验证GNSS导航接收设备的设计方案,分析GNSS导航接收设备信号接收的性能。
步骤七、根据地球到其他天体转移(如地月转移)和其他天体到地球转移(如月地转移)过程中的GNSS导航接收设备的在轨数据,完成探测器在轨飞行过程中轨道确定,同时可以根据在轨飞行结果,完成GNSS导航接收设备的性能评估,确认GNSS导航设备的在轨适用范围。
自此,完成了/实现了高轨GNSS导航接收设备在轨试验的系统设计和在轨实施。
具体以月球到地球探测器高轨GNSS导航接收设备在轨试验为例。
月地高速再入返回飞行器采用地月自由返回轨道飞行,全过程飞行时间196h,在地月转移和月地转移过程中其飞行高度和飞行速度变化范围较大,如图4和图5所示。利用地月自由返回轨道特点能够充分考核高轨GNSS导航接收设备对不同高度和飞行速度的适应能力。
月地高度再入返回飞行器高轨GNSS导航接收设备在轨试验的系统设计的主要步骤如下:
(一)确定在轨试验开始的飞行高度
(1)根据月地高速再入返回飞行器距离地面不同飞行高度时在地球赤道惯性坐标系中的位置ReH,(本例对应的飞行高度为[3000km,3500km,4000km,4500km,5000km,5500km,6000km…])根据公式③确定月地高速再入返回飞行器对GNSS导航卫星的可见性。
(2)根据分析结果,确认对至少4颗GNSS导航卫星可见时的最低飞行高度h,作为GNSS导航接收设备开机测试的飞行高度。本例h=5000km。
(二)确定GNSS导航接收天线的波束角
(1)分析全飞行过程中GNSS导航接收天线轴向与探测器-地心连线的夹角σ。本例中,根据月地高速再入返回飞行器的飞行姿态和GNSS导航接收天线的安装方向分析可知,GNSS导航接收天线轴向与探测器-地心连线的夹角σ变化范围为0°~180°之间。
(2)利用公式④和⑤计算探测器对GNSS导航卫星轨道的半张角。本例中,月地高速再入返回飞行器与地心的距离ReH在5000km至380000km之间变化,GNSS导航卫星的轨道半径R=42164km,月地高速再入返回飞行器对GNS导航卫星轨道的半张角γ最大值为180°。
(3)根据公式⑥和⑦计算得到天线的波束角β=180°。
(三)确定月地返回过程的开机时刻
根据地月转移过程中高轨GNSS导航接收设备可以接收4颗导航卫星信号时的飞行高度,确定月地返回过程中的开机高度h。本例中,高轨GNSS导航接收设备在地月转移过程中距离地面约70000km高度处能够接收的导航卫星信号少于4颗,因此月地返回过程中在距离地面高度h=70000km处进行导航接收设备在轨测试。
(四)在轨试验结果评估
月地高速再入返回飞行器在轨完成了高轨GNSS导航接收设备的在轨试验,导航设备在地月转移段探测器距地6000km~63000km完成定位,月地转移段50000km~5000km完成定位,定位精度优于100m;在地月转移段导航漏信号由强渐弱,开机定位时间不超过5min,距离增加对定位精度影响较小;在月地转移段导航漏信号由弱渐强,GNSS导航接收机开机时处于盲捕状态,开机定位时间约30min,在70000km~50000km能接收到漏信号,但最终定位高度约为50000km,随着距地距离减小,漏信号强度增加,导航定位精度提高,进一步充分验证了高轨GNSS导航接收设备对深空探测任务的适用性。
本发明中未说明部分属于本领域的公知技术。

Claims (3)

1.一种基于导航卫星漏信号进行定位在轨试验验证方法,其特征在于,包括以下步骤:
1)确定在轨测试飞行高度,通过对GNSS系统的覆盖特性分析,高于地面3000km的范围内为弱导航信号的覆盖区域,确定GNSS导航卫星接收设备在3000km以上进行在轨测试;
2)根据探测器飞行轨道动力学,获取地球、任意导航卫星和探测器在地心惯性系Oi-XiYiZi下的位置,分别计算地球→导航卫星、地球→探测器、探测器→导航卫星的距离关系;根据所述的距离关系得到探测器的不同位置,通过计算探测器与导航卫星连线和导航卫星与地球连线的夹角,判断出导航卫星的可见性;
3)根据探测器运行轨道,当探测器距离高度大于3000km时,每隔500km对探测器和所有GNSS导航卫星的可见性进行取点分析,当有4颗导航卫星满足导航卫星可见性要求时,将对应的飞行高度确定为GNSS导航接收设备在轨测试开始时刻;
4)对探测器在轨飞行姿态进行分析,通过计算导航卫星接收天线轴向与探测器-地心连线矢量的夹角以及探测器对GNSS导航卫星的半张角,确定GNSS导航接收天线的波束角范围;
5)GNSS导航接收设备在轨测试过程中,记录地球到其他天体转移过程中正常接收导航信号GNSS导航卫星小于4颗时的飞行高度,该飞行高度作为返回地球过程中GNSS导航接收设备开始在轨测试的飞行高度,从而确定返回地球过程中GNSS导航接收设备的开机工作时间;
6)记录地球到其他天体转移段探测器远离地球的过程中接收导航卫星漏信号的情况和其他天体到地球返回段探测器接近地球的过程中接收到的由弱渐强的导航漏信号的情况,比较在不同高度范围内,对于可见导航卫星信号的捕获的情况,从而对导航接收机信号捕获性能进行评估,验证分析GNSS导航接收设备信号接收的性能;
7)根据地球到其他天体转移和其他天体到地球转移过程中的GNSS导航接收设备的在轨数据,完成探测器在轨飞行过程中轨道确定,同时根据在轨飞行结果,完成GNSS导航接收设备的性能评估,确认GNSS导航设备的在轨适用范围。
2.根据权利要求1所述的基于导航卫星漏信号进行定位在轨试验验证方法,其特征在于,所述的步骤2)中所述的地球、任意导航卫星和探测器在地心惯性系Oi-XiYiZi下的位置为(Xei,Yei,Zei),(XGi,YGi,ZGi),(XHi,YHi,ZHi),所述的地球→导航卫星、地球→探测器、探测器→导航卫星的距离关系为ReG、ReH、RHG
所述的探测器与导航卫星连线和导航卫星与地球连线的夹角α:
已知a=Re+a /ReG,Re+a=地球半径+大气层高度,b为导航卫星天线波束最大幅角,均可视为常值,探测器对导航卫星可见性判据为:
a≤α≤b ③
当探测器满足公式③规定的条件时,对该导航卫星可见,否则为不可见。
3.根据权利要求1或2所述的基于导航卫星漏信号进行定位在轨试验验证方法,其特征在于,所述的步骤4)所述的导航卫星接收天线轴向与探测器-地心连线矢量ReH的夹角为σ,探测器对GNSS导航卫星的半张角为γ,导航卫星接收天线轴向为S,探测器与地心连线距离为ReH,GNSS导航卫星轨道半径为R:
当R<ReH时,γ=arcsin(R/ReH); ④
当R>ReH时,γ=180°; ⑤
所述的GNSS导航卫星接收天线的半波束角β为:
当γ+σ<180°时,β=γ+σ; ⑥
当γ+σ>180°时,β=180° ⑦。
CN201510036026.6A 2015-01-23 2015-01-23 一种基于导航卫星漏信号进行定位在轨试验验证方法 Active CN105891851B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510036026.6A CN105891851B (zh) 2015-01-23 2015-01-23 一种基于导航卫星漏信号进行定位在轨试验验证方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510036026.6A CN105891851B (zh) 2015-01-23 2015-01-23 一种基于导航卫星漏信号进行定位在轨试验验证方法

Publications (2)

Publication Number Publication Date
CN105891851A CN105891851A (zh) 2016-08-24
CN105891851B true CN105891851B (zh) 2018-06-08

Family

ID=56999580

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510036026.6A Active CN105891851B (zh) 2015-01-23 2015-01-23 一种基于导航卫星漏信号进行定位在轨试验验证方法

Country Status (1)

Country Link
CN (1) CN105891851B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108408087B (zh) * 2018-02-12 2019-07-16 北京空间技术研制试验中心 低轨长寿命载人航天器的在轨测试方法
CN111381261B (zh) * 2018-12-29 2022-05-27 广州市泰斗电子科技有限公司 一种定位解算方法、装置及卫星导航接收机
CN111273687A (zh) * 2020-02-17 2020-06-12 上海交通大学 基于gnss观测量和机间测距的多无人机协同相对导航方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0777128B1 (fr) * 1995-12-01 2001-11-07 Astrium SAS Procédé et dispositif de mesure d'attitude de satellite
CN102116641A (zh) * 2009-12-31 2011-07-06 北京控制工程研究所 一种深空自主导航星敏感器的半物理仿真试验系统及方法
CN103076618A (zh) * 2012-12-26 2013-05-01 北京空间飞行器总体设计部 一种星载测量型gnss接收机地面验证和性能评测方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0777128B1 (fr) * 1995-12-01 2001-11-07 Astrium SAS Procédé et dispositif de mesure d'attitude de satellite
CN102116641A (zh) * 2009-12-31 2011-07-06 北京控制工程研究所 一种深空自主导航星敏感器的半物理仿真试验系统及方法
CN103076618A (zh) * 2012-12-26 2013-05-01 北京空间飞行器总体设计部 一种星载测量型gnss接收机地面验证和性能评测方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
Attitude Navigation for Deep Space Explorer based on Mission Planning;Fanyu Zhao;《2014 26th Chinese Control and Decision Conference》;20140602;第2088-2093页 *
深空探测器自主技术发展现状与趋势;崔平远 等;《航空学报》;20130722;第35卷(第1期);第13-28页 *

Also Published As

Publication number Publication date
CN105891851A (zh) 2016-08-24

Similar Documents

Publication Publication Date Title
US6735523B1 (en) Process and system of coupled real-time GPS/IMU simulation with differential GPS
US7817087B2 (en) Method and apparatus for relative navigation using reflected GPS signals
CN104406605A (zh) 机载多导航源综合导航仿真系统
CN102736091B (zh) 星球表面广域探测的卫星导航方法及系统
CN108594272A (zh) 一种基于鲁棒卡尔曼滤波的抗欺骗干扰组合导航方法
CN103017772B (zh) 一种基于可观性分析的光学和脉冲星融合自主导航方法
Gaylor et al. GPS/INS Kalman filter design for spacecraft operating in the proximity of International Space Station
CN103968836A (zh) 一种基于时序伪距差分推算移动目标位置的方法及装置
CN105891851B (zh) 一种基于导航卫星漏信号进行定位在轨试验验证方法
EP3924748A1 (en) Techniques for determining geolocations
Sabatini et al. Differential global positioning system (DGPS) for flight testing
Specht Experimental studies on the relationship between HDOP and position error in the GPS system
Khanafseh et al. Autonomous airborne refueling of unmanned air vehicles using the global positioning system
CN111679242A (zh) 适用于指向在轨航天器的地面天线导引方法
US6782320B1 (en) Method and system of single-antenna determination of position, time, and attitude of a moving object by satellite navigation
CN112394381B (zh) 基于球卫星的全自主月面导航和数据通信方法
Bamford et al. A GPS receiver for lunar missions
CN111427000A (zh) 适用于地面测站天线对卫星指向的目标视向量确定方法
Qiao et al. A multiple GNSS-based orbit determination algorithm for geostationary satellites
Dougherty et al. GPS modeling for designing aerospace vehicle navigation systems
CN106019346B (zh) 基于双向通信协同定位辅助的单星定位方法
CN111366126A (zh) 地面测站天线对卫星指向的视向量计算系统
Kanamori et al. A Fingerprint-like Positioning Method Using Low Earth Orbit Satellites
Katona On mean revisit frequency of non-repeating satellite orbits with finite sensor range
CN112803985B (zh) 长基线动态星间通信和星间测量验证方法、系统

Legal Events

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