CN113742827A - Method for constructing highway slope monitoring network system based on finite difference analysis - Google Patents
Method for constructing highway slope monitoring network system based on finite difference analysis Download PDFInfo
- Publication number
- CN113742827A CN113742827A CN202111031241.9A CN202111031241A CN113742827A CN 113742827 A CN113742827 A CN 113742827A CN 202111031241 A CN202111031241 A CN 202111031241A CN 113742827 A CN113742827 A CN 113742827A
- Authority
- CN
- China
- Prior art keywords
- slope
- monitoring
- displacement
- road
- monitoring network
- 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.)
- Granted
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/13—Architectural design, e.g. computer-aided architectural design [CAAD] related to design of buildings, bridges, landscapes, production plants or roads
-
- E—FIXED CONSTRUCTIONS
- E02—HYDRAULIC ENGINEERING; FOUNDATIONS; SOIL SHIFTING
- E02D—FOUNDATIONS; EXCAVATIONS; EMBANKMENTS; UNDERGROUND OR UNDERWATER STRUCTURES
- E02D1/00—Investigation of foundation soil in situ
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
Abstract
The invention provides a method for constructing a road slope monitoring network system based on finite difference analysis, which comprises the following steps: selecting a monitoring section according to the on-site survey result of the road slope to obtain a first arrangement scheme of monitoring points; carrying out image acquisition on the road slope through unmanned aerial vehicle oblique photography, and establishing a road slope geological model; extracting three-dimensional point cloud data of a road slope geological model, and performing grid refinement to obtain a grid file; performing numerical simulation by using an intensity reduction method to obtain a slope safety factor, a deformation displacement cloud picture, a displacement contour map and a shear strain increment cloud picture; finding out the position of the slope instability, and obtaining a second arrangement scheme of the monitoring points according to the position of the slope instability; and obtaining the highway slope monitoring network according to the first arrangement scheme and the second arrangement scheme of the monitoring points. The method can solve the technical problems that monitoring point positions of the highway side slope are not reasonably arranged, and the main deformation of the landslide and the change condition of the stress concentration part cannot be accurately reflected.
Description
Technical Field
The invention relates to the technical field of highway slope monitoring, in particular to a method for constructing a highway slope monitoring network system based on finite difference analysis.
Background
More sufficient geological data and real-time slope dynamics can be obtained through road slope monitoring, so that an unstable area of a suspicious slope is judged, a sliding mode, a sliding direction and a sliding speed of the unstable slope are determined, a slope development change rule is mastered, and an important basis is provided for taking necessary protection measures.
For monitoring the road side slope, the side slope geological condition and the side slope potential failure mode need to be fully considered, the construction safety monitoring and the long-term monitoring are mainly considered, the overall situation is considered, and the key points are highlighted. In the prior art, for monitoring point locations for monitoring road side slopes, survey lines are mainly laid according to terrains and typical sections, corresponding monitoring points are laid according to the survey lines, and finally a monitoring network is formed, wherein the plurality of laid monitoring point locations are generally in a grid shape at equal intervals.
However, for slope monitoring, although the above-mentioned grid-like arrangement of monitoring points can monitor a certain section of area of a highway slope, the arrangement of monitoring points is evenly distributed and lacks pertinence, and the arrangement of monitoring points of the highway slope is not reasonable enough, and cannot accurately reflect the main deformation of the landslide and the change condition of the stress concentration portion.
Disclosure of Invention
Aiming at the defects in the prior art, the invention provides a method for constructing a road slope monitoring network system based on finite difference analysis, and aims to solve the technical problems that the monitoring point position arrangement of a road slope is not reasonable enough, and the main deformation of the landslide and the change condition of a stress concentration part cannot be accurately reflected in the prior art.
The invention adopts the technical scheme that a method for constructing a road slope monitoring network system based on finite difference analysis comprises the following steps:
carrying out on-site investigation on the road side slope;
selecting a monitoring section according to the on-site investigation result, and obtaining a first arrangement scheme of monitoring points according to the monitoring section;
acquiring images of the highway side slope through unmanned aerial vehicle oblique photography, and establishing a highway side slope geological model by using three-dimensional modeling software according to the acquired images;
extracting three-dimensional point cloud data of a road slope geological model, and performing grid refinement on the three-dimensional point cloud data to obtain a grid file;
carrying out numerical simulation by using an intensity reduction method according to the grid file to obtain a slope safety factor, a deformation displacement cloud picture, a displacement contour map and a strain increment cloud picture;
finding out the position of the side slope with instability according to the safety coefficient of the side slope, the deformation displacement cloud chart, the displacement contour map and the shearing strain increment cloud chart, and obtaining a second arrangement scheme of monitoring points according to the position of the side slope with instability;
and obtaining the highway slope monitoring network according to the first arrangement scheme and the second arrangement scheme of the monitoring points.
Furthermore, when the road side slope is surveyed on the spot, the surveyed objects comprise geological phenomena with poor landform, stratum, geological structure and rock-soil property.
Furthermore, all monitoring points are distributed at equal intervals in the first distribution scheme of the monitoring points.
Further, extracting dense point cloud of the road slope geological model as three-dimensional point cloud data.
Further, the mesh refinement processing is carried out on the three-dimensional point cloud data, and the mesh refinement processing comprises the following steps: and generating a mesh file by using 3D modeling software through a curved surface solid modeling method.
Further, numerical simulations were performed using FLAC3D software.
Further, the method for finding out the position of the instability of the side slope according to the safety coefficient of the side slope, the deformation displacement cloud picture, the displacement contour map and the shearing strain increment cloud picture comprises the following steps:
and extracting surface layer displacement and deep layer displacement of the side slope according to the deformation displacement cloud chart and the displacement contour map by using a fish language compiled by FLAC3D software, and combining the surface layer displacement and the deep layer displacement of the side slope with the safety coefficient of the side slope and the shearing strain increment cloud chart to form a displacement curve-reduction coefficient relation curve, wherein the maximum displacement is the position of instability of the side slope.
Further, highway slope monitoring network comprises the field data acquisition instrument on a plurality of monitoring points, includes: rainfall monitor, crack monitor, ground water level monitor, inclinometer.
According to the technical scheme, the beneficial technical effects of the invention are as follows:
through the technical scheme of the embodiment, a more reasonable road slope monitoring network with monitoring points can be obtained, and the main deformation of the landslide and the change condition of the stress concentration part can be accurately reflected.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings used in the description of the embodiments or the prior art will be briefly described below. Throughout the drawings, like elements or portions are generally identified by like reference numerals. In the drawings, elements or portions are not necessarily drawn to scale.
FIG. 1 is a flow chart of a method for constructing a monitoring grid architecture according to an embodiment of the present invention;
fig. 2 is a schematic view of a bedding rock slope according to an embodiment of the present invention.
Detailed Description
Embodiments of the present invention will be described in detail below with reference to the accompanying drawings. The following examples are only for illustrating the technical solutions of the present invention more clearly, and therefore are only examples, and the protection scope of the present invention is not limited thereby.
It is to be noted that, unless otherwise specified, technical or scientific terms used herein shall have the ordinary meaning as understood by those skilled in the art to which the invention pertains.
Examples
The embodiment provides a method for constructing a road slope monitoring network system based on finite difference analysis, as shown in fig. 1, the method comprises the following steps:
s1, carrying out on-site investigation on road side slope
When the road side slope is subjected to field exploration, the objects to be surveyed comprise the geological phenomena of poor landform, stratum, geological structure and rock-soil property.
S2, selecting a monitoring section according to the on-site investigation result, and planning a first arrangement scheme of monitoring points according to the monitoring section;
in a specific embodiment, the monitoring section can be selected from a potentially damaged part such as a geological condition poor, a deformation large, a fault, a crack, a dangerous rock mass and the like, or a part such as a steep part and poor stability. After a monitoring section is selected, firstly, a first arrangement scheme of monitoring points is planned on the monitoring section; in the first arrangement scheme of the monitoring points, all the monitoring points are arranged at equal intervals.
S3, carrying out image acquisition on the road slope through unmanned aerial vehicle oblique photography, and establishing a road slope geological model by using three-dimensional modeling software according to the acquired image
In a specific embodiment, the image acquisition is performed on the road slope in a manner that can be realized by any one of the existing unmanned aerial vehicle oblique photography technologies.
According to the collected images, the existing three-dimensional modeling software is used for establishing a road slope geological model, such as PhotoSacan three-dimensional modeling software.
S4, extracting three-dimensional point cloud data of the road slope geological model, and carrying out grid refinement processing on the three-dimensional point cloud data to obtain a grid file
In a specific embodiment, a three-dimensional point cloud processing software (such as Meshlab) is used for processing the road slope geological model, and dense point cloud of the road slope geological model is extracted as three-dimensional point cloud data. The extracted three-dimensional point cloud data is gridded data.
And then, carrying out mesh refinement processing on the meshed three-dimensional point cloud data to obtain a mesh file. In a specific embodiment, a mesh file is generated by a face solid modeling method using 3D modeling software (such as a Rhino plug-in), and the format of the mesh file is optimized to be f3 grid. The grid file obtained by processing according to the method in the step is basically consistent with the actual situation of the road side slope geology, and the side slope catastrophe process can be accurately reflected.
S5, performing numerical simulation by using an intensity reduction method according to the grid file to obtain a slope safety factor, a deformation displacement cloud picture, a displacement contour map and a shear strain increment cloud picture
The strength reduction method is a method for calculating a finite element of slope stability, wherein the safety coefficient of a slope is continuously reduced in the calculation, and the reduced parameters are continuously substituted into a model for repeated calculation until the model reaches the limit and is damaged, and the value before the damage is the safety coefficient of the slope. In a specific implementation mode, the FLAC3D software is used for carrying out numerical simulation, and a deformation displacement cloud picture, a displacement contour map and a shear strain increment cloud picture can be obtained by using a calculation process of intensity reduction on a grid file.
FLAC3D software is fast Lagrange finite difference analysis software, can perform three-dimensional structure stress characteristic simulation and plastic flow analysis of soil, rocks and other materials, and fits the actual structure by adjusting polyhedral units in a three-dimensional grid; the unit material can adopt a linear or nonlinear constitutive model, and when the material generates yield flow under the action of external force, the grid can correspondingly deform and move. The FLAC3D employs an explicit Lagrangian algorithm and a hybrid-discrete partitioning technique, which can very accurately simulate plastic failure and flow of materials.
S6, finding out the position of the side slope with instability according to the safety coefficient of the side slope, the deformation displacement cloud chart, the displacement contour map and the shearing strain increment cloud chart, and obtaining a second arrangement scheme of monitoring points according to the position of the side slope with instability
In a specific implementation mode, a fish language compiled by FLAC3D software is used, slope surface layer displacement and deep layer displacement are extracted according to a deformation displacement cloud chart and a displacement contour map, the slope surface layer displacement and the deep layer displacement are combined with a slope safety coefficient and a shear strain increment cloud chart to form a displacement curve-reduction coefficient relation curve, and the maximum displacement position is a position where the slope is unstable, namely the position where the slope is most likely to be damaged. And the position of the slope with instability is taken as a monitoring point, and the arrangement positions of the monitoring points form a second arrangement scheme of the monitoring points.
The bedding rock slope is used as an example to illustrate the instability of the slope:
the schematic diagram of the bedding rock slope is shown in fig. 2, the slope body of the bedding rock slope is in the same direction with the rock stratum surface, the inclined included angle is less than 30 degrees, the rock mass is mostly in interbedded and interbedded dislocation zone, and is often a penetrating structural surface, and the rock mass slope is controlled by the combination relationship of the rock stratum inclination angle and the slope angle, the structural surface development condition and the strength. The destabilization pattern is shown in table 1:
TABLE 1 layered slope destabilization mode
Because different slopes may have different instability modes, the positions of the slopes can be found by combining the results (stress, strain, displacement and other changes) obtained by numerical simulation with different slope failure modes. When a curve of the relationship between the displacement or the maximum displacement of a certain part and the reduction coefficient is established in the calculation process, the inflection point on the curve can be used as a critical point of the slope in damage; due to the property of the rock-soil mass, the occurrence of the damage of the rock-soil mass is closely related to the expansion and distribution of the plastic region, and the plastic strain region is communicated when the slope is damaged.
S7, obtaining the highway side slope monitoring network according to the first arrangement scheme and the second arrangement scheme of the monitoring points
The second arrangement scheme of the monitoring points is that theoretical arrangement positions of the monitoring points are calculated according to the model, and the theoretical arrangement positions of the monitoring points cannot be actually constructed due to field geological conditions during actual arrangement, so that side slope sections are selected to be arranged according to field investigation results, such as geological conditions, hydrological conditions, cracks and the like of a road side slope field, and different monitoring sections form a side slope monitoring network.
In a specific implementation mode, the constructed road slope monitoring network is composed of various field data acquisition instruments on a plurality of monitoring points, and comprises the following steps: rainfall monitor, crack monitor, ground water level monitor, inclinometer.
Through the technical scheme of the embodiment, a more reasonable road slope monitoring network with monitoring points can be obtained, and the main deformation of the landslide and the change condition of the stress concentration part can be accurately reflected.
Finally, it should be noted that: the above examples are only intended to illustrate the technical solution of the present invention, but not to limit it; while the invention has been described in detail and with reference to the foregoing embodiments, it will be understood by those skilled in the art that: the technical solutions described in the foregoing embodiments may still be modified, or some or all of the technical features may be equivalently replaced; such modifications and substitutions do not depart from the spirit and scope of the present invention, and they should be construed as being included in the following claims and description.
Claims (8)
1. A method for constructing a road slope monitoring network system based on finite difference analysis is characterized by comprising the following steps:
carrying out on-site investigation on the road side slope;
selecting a monitoring section according to the on-site investigation result, and obtaining a first arrangement scheme of monitoring points according to the monitoring section;
acquiring images of the highway side slope through unmanned aerial vehicle oblique photography, and establishing a highway side slope geological model by using three-dimensional modeling software according to the acquired images;
extracting three-dimensional point cloud data of a road slope geological model, and performing grid refinement on the three-dimensional point cloud data to obtain a grid file;
carrying out numerical simulation by using an intensity reduction method according to the grid file to obtain a slope safety factor, a deformation displacement cloud picture, a displacement contour map and a strain increment cloud picture;
finding out the position of the side slope with instability according to the safety coefficient of the side slope, the deformation displacement cloud chart, the displacement contour map and the shearing strain increment cloud chart, and obtaining a second arrangement scheme of monitoring points according to the position of the side slope with instability;
and obtaining the highway slope monitoring network according to the first arrangement scheme and the second arrangement scheme of the monitoring points.
2. The method for constructing a road slope monitoring network system based on finite difference analysis according to claim 1, wherein when the road slope is surveyed on the spot, the surveyed objects include topographic features, strata, geological structures and geological phenomena with poor geotechnical properties.
3. The finite difference analysis-based road slope monitoring network system construction method according to claim 1, wherein the monitoring points are distributed at equal intervals in the first distribution scheme of the monitoring points.
4. The finite difference analysis-based road slope monitoring network system construction method according to claim 1, characterized in that dense point cloud of a road slope geological model is extracted as the three-dimensional point cloud data.
5. The finite difference analysis-based road slope monitoring network system construction method according to claim 1, wherein the mesh refinement processing is performed on the three-dimensional point cloud data, and comprises the following steps: and generating a mesh file by using 3D modeling software through a curved surface solid modeling method.
6. The finite difference analysis-based road slope monitoring network architecture construction method according to claim 1, wherein the numerical simulation is performed by using FLAC3D software.
7. The method for constructing the road slope monitoring network system based on the finite difference analysis according to claim 1, wherein the position of the slope instability is found according to the slope safety factor, the deformation displacement cloud chart, the displacement contour chart and the shear strain increment cloud chart, and the method comprises the following steps:
and extracting surface layer displacement and deep layer displacement of the side slope according to the deformation displacement cloud chart and the displacement contour map by using a fish language compiled by FLAC3D software, and combining the surface layer displacement and the deep layer displacement of the side slope with the safety coefficient of the side slope and the shearing strain increment cloud chart to form a displacement curve-reduction coefficient relation curve, wherein the maximum displacement is the position of instability of the side slope.
8. The method for constructing a road slope monitoring network system based on finite difference analysis according to claim 1, wherein the road slope monitoring network is composed of field data acquisition instruments on a plurality of monitoring points, and comprises: rainfall monitor, crack monitor, ground water level monitor, inclinometer.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111031241.9A CN113742827B (en) | 2021-09-03 | 2021-09-03 | Highway slope monitoring network system construction method based on finite difference analysis |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111031241.9A CN113742827B (en) | 2021-09-03 | 2021-09-03 | Highway slope monitoring network system construction method based on finite difference analysis |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113742827A true CN113742827A (en) | 2021-12-03 |
CN113742827B CN113742827B (en) | 2023-05-19 |
Family
ID=78735429
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111031241.9A Active CN113742827B (en) | 2021-09-03 | 2021-09-03 | Highway slope monitoring network system construction method based on finite difference analysis |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113742827B (en) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114722560A (en) * | 2021-12-23 | 2022-07-08 | 西南交通大学 | Safety coefficient determination method based on soil slope strength reduction method |
CN116229371A (en) * | 2023-03-13 | 2023-06-06 | 浙江省交通运输科学研究院 | Sky-ground integrated highway side slope safety prevention and control method |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106405675A (en) * | 2016-08-25 | 2017-02-15 | 山东科技大学 | Dynamic monitoring system and method for early warning against slope slide of tailing pond of strip mining pit |
CN110162932A (en) * | 2019-06-21 | 2019-08-23 | 招商局重庆交通科研设计院有限公司 | The calculation method of automatic monitoring slope stability based on finite element lattice |
CN110378056A (en) * | 2019-07-25 | 2019-10-25 | 四川大学 | It is a kind of for the slope stability measuring method of slope geological mechanical model and application |
CN111784829A (en) * | 2020-06-01 | 2020-10-16 | 重庆地质矿产研究院 | Three-dimensional geological model-based landslide monitoring key point arrangement method |
CN113192084A (en) * | 2021-05-07 | 2021-07-30 | 中国公路工程咨询集团有限公司 | Machine vision-based highway slope micro-displacement deformation monitoring method |
-
2021
- 2021-09-03 CN CN202111031241.9A patent/CN113742827B/en active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106405675A (en) * | 2016-08-25 | 2017-02-15 | 山东科技大学 | Dynamic monitoring system and method for early warning against slope slide of tailing pond of strip mining pit |
CN110162932A (en) * | 2019-06-21 | 2019-08-23 | 招商局重庆交通科研设计院有限公司 | The calculation method of automatic monitoring slope stability based on finite element lattice |
CN110378056A (en) * | 2019-07-25 | 2019-10-25 | 四川大学 | It is a kind of for the slope stability measuring method of slope geological mechanical model and application |
CN111784829A (en) * | 2020-06-01 | 2020-10-16 | 重庆地质矿产研究院 | Three-dimensional geological model-based landslide monitoring key point arrangement method |
CN113192084A (en) * | 2021-05-07 | 2021-07-30 | 中国公路工程咨询集团有限公司 | Machine vision-based highway slope micro-displacement deformation monitoring method |
Non-Patent Citations (3)
Title |
---|
刘中帅: "抗滑键加固渝黔高速公路顺层岩质边坡模拟研究", 《重庆大学学报(自然科学版)》 * |
张颖 等: "FLAC3D在边坡数值分析及位移检测中的应用", 《探矿工程》 * |
杨真子: "基于强降雨的膨胀土边坡稳定性分析", 《公路交通技术》 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114722560A (en) * | 2021-12-23 | 2022-07-08 | 西南交通大学 | Safety coefficient determination method based on soil slope strength reduction method |
CN116229371A (en) * | 2023-03-13 | 2023-06-06 | 浙江省交通运输科学研究院 | Sky-ground integrated highway side slope safety prevention and control method |
CN116229371B (en) * | 2023-03-13 | 2023-10-13 | 浙江省交通运输科学研究院 | Sky-ground integrated highway side slope safety prevention and control method |
Also Published As
Publication number | Publication date |
---|---|
CN113742827B (en) | 2023-05-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Li et al. | Quantitative assessment for the rockfall hazard in a post-earthquake high rock slope using terrestrial laser scanning | |
CN114663627B (en) | Mine digital model building method based on three-dimensional point cloud database | |
Wang et al. | Cliff face rock slope stability analysis based on unmanned arial vehicle (UAV) photogrammetry | |
CN113742827A (en) | Method for constructing highway slope monitoring network system based on finite difference analysis | |
CN115690354B (en) | Dynamic control method for shallow tunnel construction based on three-dimensional live-action numerical analysis | |
CN108842821B (en) | Calculation method for reasonable buried depth of submarine tunnel constructed by drilling and blasting method | |
Xie et al. | Three-dimensional stability evaluation of landslides and a sliding process simulation using a new geographic information systems component | |
US20210109248A1 (en) | System, method, and device for real-time sinkhole detection | |
CN115062368A (en) | Digital twin body and fine modeling system and method for long tunnel shed | |
CN110827405A (en) | Digital remote sensing geological mapping method and system | |
CN113536414B (en) | Rock slope stability analysis method, system and medium based on three-dimensional modeling | |
CN116486025A (en) | Urban geological data processing platform based on big data cloud computing technology | |
CN116776444A (en) | Digital twinning-based immersed tube tunnel construction method and system | |
Dong et al. | Combined methodology for three-dimensional slope stability analysis coupled with time effect: a case study in Germany | |
VERMA | Micro-Zonation of Landslide Hazards Between Aizawl City and Lengpui Airport, Mizoram, India, Using Geo-informatics | |
Wang et al. | Automatic identification of rock discontinuity and stability analysis of tunnel rock blocks using terrestrial laser scanning | |
Fan et al. | Advanced stability analysis method for the tunnel face in jointed rock mass based on DFN-DEM | |
Yan et al. | Detection and numerical simulation of potential hazard in oil pipeline areas based on UAV surveys | |
Gao et al. | 3D visualization monitoring and early warning of surface deformation in subsidence area based on GIS | |
CN115713599A (en) | Method for accurately monitoring and preventing underground coal mining mine | |
CN114608661B (en) | Certainty index evaluation method for landslide rock mass structure in mountain area | |
Stead et al. | Rock slope characterization and geomechanical modelling | |
CN113269886B (en) | Slope three-dimensional digital twin model building method based on multi-source data fusion | |
CN116797755B (en) | Modeling method for multi-time-space three-dimensional geological structure of mixed rock zone | |
Saylan | Three dimensional modeling and monitoring of a man made slope using images taken from unmanned aerial vehicle |
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 |