CN112629526B - Tight combination navigation method for Beidou precise single-point positioning and inertial system - Google Patents

Tight combination navigation method for Beidou precise single-point positioning and inertial system Download PDF

Info

Publication number
CN112629526B
CN112629526B CN202011302395.2A CN202011302395A CN112629526B CN 112629526 B CN112629526 B CN 112629526B CN 202011302395 A CN202011302395 A CN 202011302395A CN 112629526 B CN112629526 B CN 112629526B
Authority
CN
China
Prior art keywords
carrier
updated
ins
ambiguity
equation
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
CN202011302395.2A
Other languages
Chinese (zh)
Other versions
CN112629526A (en
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.)
Information Engineering University of PLA Strategic Support Force
Original Assignee
Information Engineering University of PLA Strategic Support Force
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 Information Engineering University of PLA Strategic Support Force filed Critical Information Engineering University of PLA Strategic Support Force
Priority to CN202011302395.2A priority Critical patent/CN112629526B/en
Publication of CN112629526A publication Critical patent/CN112629526A/en
Application granted granted Critical
Publication of CN112629526B publication Critical patent/CN112629526B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01CMEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C21/00Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
    • G01C21/10Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by using measurements of speed or acceleration
    • G01C21/12Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by using measurements of speed or acceleration executed aboard the object being navigated; Dead reckoning
    • G01C21/16Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by using measurements of speed or acceleration executed aboard the object being navigated; Dead reckoning by integrating acceleration or speed, i.e. inertial navigation
    • G01C21/165Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by using measurements of speed or acceleration executed aboard the object being navigated; Dead reckoning by integrating acceleration or speed, i.e. inertial navigation combined with non-inertial navigation instruments
    • 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
    • G01S19/00Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/38Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
    • G01S19/39Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system the satellite radio beacon positioning system transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/42Determining position
    • G01S19/45Determining position by combining measurements of signals from the satellite radio beacon positioning system with a supplementary measurement
    • G01S19/47Determining position by combining measurements of signals from the satellite radio beacon positioning system with a supplementary measurement the supplementary measurement being an inertial measurement, e.g. tightly coupled inertial

Landscapes

  • Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Automation & Control Theory (AREA)
  • Navigation (AREA)
  • Position Fixing By Use Of Radio Waves (AREA)

Abstract

The invention relates to a tightly combined navigation method of a Beidou precise single-point positioning and inertial system, and belongs to the technical field of combined navigation positioning. When Kalman filtering is adopted to solve a state equation and an observation equation, only the state quantity and the covariance matrix related to the INS are updated when only time is updated, parameters related to the GNSS are not updated, and the ambiguity parameters and the ambiguity covariance of the last epoch are introduced for updating when measurement is updated. According to the method, the time updating frequency is faster than the measurement updating frequency, and only the state quantity and the covariance matrix related to the INS are updated when the time updating is controlled, so that the problems of poor calculation efficiency and poor stability caused by overlarge orders of the state transition matrix when a plurality of satellites are avoided, and the calculation efficiency and the instantaneity are improved.

Description

Tight combination navigation method for Beidou precise single-point positioning and inertial system
Technical Field
The invention relates to a tightly combined navigation method of a Beidou precise single-point positioning and inertial system, and belongs to the technical field of combined navigation positioning.
Background
And at the end of 7 months in 2020, the Beidou No. three global satellite navigation system formally opens services, and marks the globalization of satellite navigation opening in China. The Beidou satellite navigation system provides all-weather, all-day and high-precision positioning navigation and time service for Chinese and global users, and is a national important space-time infrastructure. At present, the Beidou system comprehensively serves various industries such as transportation, public safety, disaster relief and reduction, agriculture, forestry, animal husbandry, urban management and the like, and is integrated with national core infrastructure construction such as electric power, finance, communication and the like. However, the "vulnerability" of satellite signals results in limited application in complex environments. The signal is easy to be blocked, multipath is complex, and the signal is easy to be disturbed and deceptively produced in an electromagnetic environment, so that the signal is interrupted, and the continuous positioning capability is lost. The inertial navigation system outputs the angular velocity and acceleration of the motion carrier by means of the gyroscope and the accelerometer, and outputs the position, velocity and attitude information of the motion carrier in an integral mode, so that the inertial navigation system has the advantages of independence, no external interference and the like, but because the inertial navigation is a recursive navigation mode, the error can be gradually increased along with time.
The satellite/inertial combination system (GNSS/INS) can realize the complementary advantages of the two, and in order to obtain high-precision position information, the common combination means are DGNSS/INS combination and PPP/INS combination. The PPP/INS combination can obtain a high-precision navigation result by using only one user receiver, and has a wider application prospect and a wider application prospect. The extended Kalman filtering algorithm is generally used in the satellite/inertial combination system to perform system fusion solution and estimate relevant parameters of the satellite navigation system (GNSS) and the Inertial Navigation System (INS). Because the satellite ambiguity parameters need to be estimated, the order of a correlation matrix is increased during Kalman filtering, the complexity of a model is increased, the data processing efficiency is reduced, and the real-time calculation is not facilitated, and the data calculation by post-processing software is inconvenient.
Disclosure of Invention
The invention aims to provide a tightly combined navigation method of a Beidou precise single-point positioning and inertial system, which aims to solve the problems of complex data processing and low instantaneity in the existing tightly combined navigation process.
The invention provides a tightly combined navigation method of a Beidou precise single-point positioning and inertial system, which aims to solve the technical problems, and comprises the following steps:
1) Acquiring satellite observation data, ephemeris data, clock error data, correction file data and inertial data;
2) Establishing a state equation and an observation equation of a tight combination according to the acquired data;
3) And solving the state equation and the observation equation by adopting Kalman filtering, wherein in the solving process, only the state quantity and covariance related to the INS are updated when the time is updated, and the ambiguity parameter and the ambiguity covariance of the last epoch are introduced when the measurement is updated.
When Kalman filtering is adopted to solve a state equation and an observation equation, only the state quantity and the covariance matrix related to the INS are updated when only time is updated, parameters related to the GNSS are not updated, and the ambiguity parameters and the ambiguity covariance of the last epoch are introduced for updating when measurement is updated. According to the method, the time updating frequency is faster than the measurement updating frequency, and only the state quantity and the covariance matrix related to the INS are updated when the time updating is controlled, so that the problems of poor calculation efficiency and poor stability caused by overlarge orders of the state transition matrix when a plurality of satellites are avoided, and the calculation efficiency and the instantaneity are improved.
Further, the state equation is:
wherein the method comprises the steps ofA directional cosine matrix from the carrier b system to the geocentric earth system; />Representing accelerometer specific force output; />Representing a gyro angular rate output; δψ represents the attitude misalignment angle error; />Representing the rotation angle rate of the earth->A constructed oblique symmetrical array; epsilon a And epsilon g A system noise vector representing accelerometer and gyro noise; δr, δv, δφ, δb a 、δb g 、δT w 、δdt r And delta N IF The carrier position, carrier velocity, carrier attitude, carrier acceleration zero bias, carrier gyro zero bias, tropospheric wet delay, receiver clock bias, and error values of ambiguity parameters, respectively.
Further, the observation equation is:
wherein P is the pseudo-range observed quantity; l is the observed quantity of carrier phase; d is Doppler observed quantity; subscript IF is the ionosphere observation value and subscript INS is the inertial prediction value; e is a 1x3 matrix formed by the directional cosine of the receiver to the satellite; delta represents the error correction; t (T) w Represents tropospheric wet delay; dt (dt) r Representing receiver clock skew; n represents the ambiguity real number solution of the satellite; epsilon represents the observed noise of each observed quantity; δr, δv,δb a 、δb g 、δT w 、δdt r And delta N IF The carrier position, carrier velocity, carrier attitude, carrier acceleration zero bias, carrier gyro zero bias, tropospheric wet delay, receiver clock bias, and error values of ambiguity parameters, respectively.
Further, the time update is determined according to the output frequency of the INS, and the measurement update is determined according to the output frequency of the GNSS.
Drawings
FIG. 1 is a flow chart of a method of tightly integrated navigation of a Beidou precise single point positioning and inertial system of the present invention;
FIG. 2 is a flowchart of a specific implementation of a tightly combined model of the Beidou precise single-point positioning and inertial system of the present invention.
Detailed Description
The following describes the embodiments of the present invention further with reference to the drawings.
In the invention, only the state quantity related to INS and the covariance matrix thereof are updated in the time update of Kalman filtering, and the measurement equation is constructed by adopting PPP technology when the measurement is updated, and the related parameters of GNSS and INS are estimated. The implementation flow of the method is shown in fig. 1, and the specific implementation process is as follows.
1. And acquiring data related to the tightly combined navigation positioning, and carrying out INS mechanical arrangement.
The tightly combined navigation system formed by the GNSS and the INS comprises a Global Navigation Satellite System (GNSS) (the Beidou navigation system BDS can be adopted in the embodiment) and an Inertial Navigation System (INS), so that the acquired data comprise observation data, ephemeris data, clock error data, various correction file data and inertial measurement data of a satellite receiver.
The INS mechanics arrangement refers to a process of converting the speed increment and the angle increment output by inertial navigation into a navigation coordinate system and recursively calculating the initial position, the speed and the attitude to obtain the position, the speed and the attitude at the next moment. Because inertial navigation is essentially a gradual recursive process, it is also often referred to as INS mechanics orchestration.
2. And constructing a state equation and an observation equation, and solving by adopting Kalman filtering.
The integrated navigation system has the essential function of enabling a user to know the position, the speed and the gesture. Thus, the parameters to be solved are position r, velocity v, poseAt the same time, the inertial navigation system contains acceleration zero offset b a Zero offset b of gyro g The navigation accuracy of the inertial navigation part in the integrated navigation system is affected to a great extent, and therefore, the solution is also needed. Representing all the above parameters of interest by vector x, easily available +.>Wherein each parameter is a three-dimensional vector, and thus the state parameter x is a fifteen-dimensional vector.
The GNSS/INS tightly integrated navigation data fusion processing generally adopts an extended Kalman filtering algorithm, and in order to avoid the problem of running of an integrated navigation system caused by unstable observation coarse difference and state model, the extended Kalman filtering algorithm is provided with an robust self-adaptive algorithm, and the robustness of the integrated navigation is enhanced in a common way in engineering. A specific implementation flow of the tightly integrated navigation model is shown in fig. 2.
The satellite observation equation is linearized, and is unfolded near the initial state value, and the state parameter to be estimated becomes an error parameter of the state, namely, an error vector delta x= [ delta r, delta v, delta phi, delta b of the state parameter is used in the operation process a ,δb g ] T
To complete the extended Kalman filter solution, the state equation and the observation equation are first constructed, mainly the state equation of the form formula (1) and the observation equation of the form formula (2).
L=Hδx+e (2)
And then solving according to the standard solving step of the extended Kalman filtering to obtain a state parameter delta x.
For a general integrated navigation system, the state equation is determined by a dynamics model of an inertial navigation system, and can be generally expressed as formula (3).
In the method, in the process of the invention,a directional cosine matrix from the carrier b system to the geocentric earth system; />Representing accelerometer specific force output; />Representing a gyro angular rate output; δψ represents the attitude misalignment angle error; />Representing the rotation angle rate of the earth->A constructed oblique symmetrical array; epsilon a And epsilon g Representing a system noise vector composed of accelerometer and gyro noise.
δx=[δr,δv,δφ,δb a ,δb g ] T Is the state quantity related to INS, and the covariance isEssentially, a diagonal matrix is formed by the state quantities, and uncertainty of each quantity is represented. In the integrated navigation model, it is necessary to pass all the way to the next time. In the time update, only the covariance of the INS-related state quantity is updated, while the covariance of the GNSS-related state quantity remains unchanged, and is updated again at the time of measurement update.
(3) The formula is constructed according to an error equation of the INS, can be classified as a state quantity related to the INS, and can be abbreviated as:
in the PPP/INS tight combination, the receiver clock error, tropospheric wet delay and ambiguity parameters are also estimated, so the state equation needs to be extended, and the specific form can be expressed as:
after expanding the state quantity of the GNSS part, obtaining a shorthand form of a tightly combined state equation:
based on the above analysis, the observation equation for the PPP/INS tight combination is shown as follows:
also can be abbreviated as:
wherein P is the pseudo-range observed quantity; l is the observed quantity of carrier phase; d is Doppler observed quantity; subscript IF is the ionosphere observation value and subscript INS is the inertial prediction value; e is a 1x3 matrix formed by the directional cosine of the receiver to the satellite; delta represents the error correction; t (T) w Represents tropospheric wet delay; dt (dt) r Representing receiver clock skew; n represents the ambiguity real number solution of the satellite; epsilon represents the observed noise of each observed quantity.
The observation equation of PPP/INS tight combination is solved by Kalman filtering, in Kalman filtering, the time update is calculated according to the output frequency of INS, the general frequency is 100Hz, the measurement update is calculated according to the output frequency of GNSS, the general receiver frequency is 1Hz, namely, the measurement update is carried out for 100 times under the general condition. If the update is performed according to equation (6) in the time update, the state transition matrix dimension is very large, and the calculation efficiency of 100 times is relatively low. Therefore, in the time update, only the state quantity and covariance related to the INS are updated, and the ambiguity parameter and the ambiguity covariance of the last epoch are introduced during measurement update, so that the operation efficiency is improved. The specific implementation process is as follows:
assume that the measurement update of the k-1 epoch filter period yields the state parametersAnd its covariance matrix P k-1
Wherein:
performing Kalman filtering time updating: calculating state prediction vectors and covariance matrixes of the k-1 to k periods:
in the method, in the process of the invention,representing an INS-related system state transition matrix for periods k-1 through k; q (Q) k-1 A covariance matrix representing system noise, wherein:
and (3) performing time updating according to the above, and substituting the ambiguity state quantity and covariance of the k-1 epoch when the k epoch is measured and updated, so as to expand the state quantity and covariance of the k epoch:
update Kalmam filtered measurements:
wherein:
through the process, the method and the device can avoid the problems of poor calculation efficiency and stability caused by overlarge orders of the state transition matrix when a plurality of satellites are used, and can meet the requirement of instantaneity.

Claims (1)

1. The tightly combined navigation method of the Beidou precise single-point positioning and inertial system is characterized by comprising the following steps of:
1) Acquiring satellite observation data, ephemeris data, clock error data, correction file data and inertial data;
2) Establishing a state equation and an observation equation of a tight combination according to the acquired data; the state equation is as follows:
wherein the method comprises the steps ofA directional cosine matrix from the carrier b system to the geocentric earth system; />Representing accelerometer specific force output; />Representing a gyro angular rate output; δψ represents the attitude misalignment angle error; />Representing the rotation angle rate of the earth->A constructed oblique symmetrical array; epsilon a And epsilon g A system noise vector representing accelerometer and gyro noise; δr, δv, δφ, δb a 、δb g 、δT w 、δdt r And delta N IF The error values of carrier position, carrier speed, carrier posture, carrier acceleration zero bias, carrier gyro zero bias, tropospheric wet delay, receiver clock error and ambiguity parameters are respectively;
the observation equation is as follows:
wherein P is the pseudo-range observed quantity; l is the observed quantity of carrier phase; d is Doppler observed quantity; subscript IF is the ionosphere observation value and subscript INS is the inertial prediction value; e is a 1x3 matrix formed by the directional cosine of the receiver to the satellite; delta representsError correction; t (T) w Represents tropospheric wet delay; dt (dt) r Representing receiver clock skew; n represents the ambiguity real number solution of the satellite; epsilon represents the observed noise of each observed quantity; δr, δv,δb a 、δb g 、δT w 、δdt r And delta N IF The error values of carrier position, carrier speed, carrier posture, carrier acceleration zero bias, carrier gyro zero bias, tropospheric wet delay, receiver clock error and ambiguity parameters are respectively;
3) Solving the state equation and the observation equation by adopting Kalman filtering, wherein in the solving process, only the state quantity and covariance related to the INS are updated when the time is updated, and the ambiguity parameter and the ambiguity covariance of the last epoch are introduced when the measurement is updated;
the time update is determined according to the output frequency of the INS, and the measurement update is determined according to the output frequency of the GNSS.
CN202011302395.2A 2020-11-19 2020-11-19 Tight combination navigation method for Beidou precise single-point positioning and inertial system Active CN112629526B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011302395.2A CN112629526B (en) 2020-11-19 2020-11-19 Tight combination navigation method for Beidou precise single-point positioning and inertial system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011302395.2A CN112629526B (en) 2020-11-19 2020-11-19 Tight combination navigation method for Beidou precise single-point positioning and inertial system

Publications (2)

Publication Number Publication Date
CN112629526A CN112629526A (en) 2021-04-09
CN112629526B true CN112629526B (en) 2023-10-31

Family

ID=75304347

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011302395.2A Active CN112629526B (en) 2020-11-19 2020-11-19 Tight combination navigation method for Beidou precise single-point positioning and inertial system

Country Status (1)

Country Link
CN (1) CN112629526B (en)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113203414B (en) * 2021-05-21 2023-06-13 北京交通大学 Train positioning method based on GPS+BDS PPP/IMU tight combination
CN117784190A (en) * 2023-12-25 2024-03-29 开普勒卫星科技(武汉)有限公司 Root mean square information filtering method, system and related device based on time correlation

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102508277A (en) * 2011-10-27 2012-06-20 中国矿业大学 Precise point positioning and inertia measurement tightly-coupled navigation system and data processing method thereof
EP2722647A1 (en) * 2012-10-18 2014-04-23 Leica Geosystems AG Surveying System and Method
CN107422354A (en) * 2017-09-25 2017-12-01 武汉大学 A kind of PPP/SINS tight integration positioning and orientation methods that fuzziness is fixed
CN108226985A (en) * 2017-12-25 2018-06-29 北京交通大学 Train Combinated navigation method based on Static Precise Point Positioning
CN111694040A (en) * 2020-04-08 2020-09-22 中国人民解放军战略支援部队信息工程大学 Positioning method and device of satellite/inertial integrated navigation system

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8018377B2 (en) * 2009-01-23 2011-09-13 Her Majesty The Queen In Right Of Canada As Represented By The Minister Of Natural Resources Decoupled clock model with ambiguity datum fixing

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102508277A (en) * 2011-10-27 2012-06-20 中国矿业大学 Precise point positioning and inertia measurement tightly-coupled navigation system and data processing method thereof
EP2722647A1 (en) * 2012-10-18 2014-04-23 Leica Geosystems AG Surveying System and Method
CN107422354A (en) * 2017-09-25 2017-12-01 武汉大学 A kind of PPP/SINS tight integration positioning and orientation methods that fuzziness is fixed
CN108226985A (en) * 2017-12-25 2018-06-29 北京交通大学 Train Combinated navigation method based on Static Precise Point Positioning
CN111694040A (en) * 2020-04-08 2020-09-22 中国人民解放军战略支援部队信息工程大学 Positioning method and device of satellite/inertial integrated navigation system

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
三频差分GNSS/INS紧组合模型;肖凯等;《中国惯性技术学报》;20180430(第02期);全文 *
惯导辅助下的动态精密单点定位算法研究;韩舒文等;《全球定位系统》;20171231(第06期);全文 *
模糊度固定解PPP/INS紧组合模型;刘帅等;《系统工程与电子技术》;20161031(第10期);第2389-2394页 *

Also Published As

Publication number Publication date
CN112629526A (en) 2021-04-09

Similar Documents

Publication Publication Date Title
Han et al. A novel initial alignment scheme for low-cost INS aided by GPS for land vehicle applications
US10082583B2 (en) Method and apparatus for real-time positioning and navigation of a moving platform
US8164514B1 (en) Method and apparatus for fusing referenced and self-contained displacement measurements for positioning and navigation
CN103235328B (en) GNSS (global navigation satellite system) and MEMS (micro-electromechanical systems) integrated navigation method
Zhang et al. New optimal smoothing scheme for improving relative and absolute accuracy of tightly coupled GNSS/SINS integration
CN113203418B (en) GNSSINS visual fusion positioning method and system based on sequential Kalman filtering
CN107656300B (en) Satellite/inertia ultra-tight combination method based on Beidou/GPS dual-mode software receiver
JP2003043129A (en) Two-stage filter using multiple motion models for use in global positioning system
CN112629526B (en) Tight combination navigation method for Beidou precise single-point positioning and inertial system
AU2002315595A1 (en) An inertial/GPS navigation system
WO2002099453A2 (en) An inertial/gps navigation system
CN108344415A (en) A kind of integrated navigation information fusion method
Hide et al. GPS and low cost INS integration for positioning in the urban environment
Wendel et al. Time-differenced carrier phase measurements for tightly coupled GPS/INS integration
CN110646822B (en) Integer ambiguity Kalman filtering algorithm based on inertial navigation assistance
Anbu et al. Integration of inertial navigation system with global positioning system using extended kalman filter
Lyu et al. Optimal time difference-based TDCP-GPS/IMU navigation using graph optimization
Zou et al. A nonlinear transfer alignment of distributed POS based on adaptive second-order divided difference filter
CN114894181A (en) Real-time autonomous combined navigation positioning method and device
Yi On improving the accuracy and reliability of GPS/INS-based direct sensor georeferencing
Nagui et al. Improved GPS/IMU loosely coupled integration scheme using two kalman filter-based cascaded stages
CN115683094A (en) Vehicle-mounted double-antenna tight coupling positioning method and system in complex environment
CN114167472A (en) INS assisted GNSS PPP precise dynamic navigation positioning method and system
Abd Rabbou et al. Integration of multi-constellation GNSS precise point positioning and MEMS-based inertial systems using tightly coupled mechanization
TW448304B (en) Fully-coupled positioning process and system

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