CN112799054A - Method for acquiring multi-period underground three-dimensional form of dynamic crack based on ground penetrating radar - Google Patents

Method for acquiring multi-period underground three-dimensional form of dynamic crack based on ground penetrating radar Download PDF

Info

Publication number
CN112799054A
CN112799054A CN202011566874.5A CN202011566874A CN112799054A CN 112799054 A CN112799054 A CN 112799054A CN 202011566874 A CN202011566874 A CN 202011566874A CN 112799054 A CN112799054 A CN 112799054A
Authority
CN
China
Prior art keywords
crack
underground
film
ground penetrating
penetrating radar
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
Application number
CN202011566874.5A
Other languages
Chinese (zh)
Other versions
CN112799054B (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.)
Shaanxi Shanbei Mining Hanjiawan Coal Co ltd
China University of Mining and Technology Beijing CUMTB
Original Assignee
Shaanxi Shanbei Mining Hanjiawan Coal Co ltd
China University of Mining and Technology Beijing CUMTB
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 Shaanxi Shanbei Mining Hanjiawan Coal Co ltd, China University of Mining and Technology Beijing CUMTB filed Critical Shaanxi Shanbei Mining Hanjiawan Coal Co ltd
Priority to CN202011566874.5A priority Critical patent/CN112799054B/en
Publication of CN112799054A publication Critical patent/CN112799054A/en
Application granted granted Critical
Publication of CN112799054B publication Critical patent/CN112799054B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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
    • G01S13/00Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
    • G01S13/88Radar or analogous systems specially adapted for specific applications
    • G01S13/885Radar or analogous systems specially adapted for specific applications for ground probing
    • 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/36Means for anti-jamming, e.g. ECCM, i.e. electronic counter-counter measures

Landscapes

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

Abstract

The invention discloses a method for acquiring a multi-period underground three-dimensional form of a dynamic fracture based on a ground penetrating radar, which comprises the following steps: laying a film on the surface of the area where the crack is located, and ensuring that the surface of the crack is completely covered by the film; pouring liquid on the film; aligning a ground penetrating radar to the crack at the position close to the ground surface of the crack, and scanning at intervals along the vertical direction of the trend of the line of the crack as a path to obtain electromagnetic wave data of a plurality of sections of the underground form of the crack at the present stage; pumping out the liquid for storage and removing the film; carrying out data processing on all the acquired electromagnetic wave data in the current period, and drawing and obtaining the contour of the crack at each section in the current period by a visual interpretation method; constructing a three-dimensional model of the current underground form of the crack by utilizing the contour of the crack to obtain the current underground three-dimensional form of the crack and related data thereof; and continuously repeating the steps to obtain the multi-stage underground three-dimensional form and the related data thereof in the dynamic change process of the crack development. The method can nondestructively acquire the multi-stage underground form of the dynamic fracture.

Description

Method for acquiring multi-period underground three-dimensional form of dynamic crack based on ground penetrating radar
Technical Field
The invention relates to the field of land reclamation of coal mining subsidence lands, in particular to a method for acquiring multi-period underground three-dimensional forms of dynamic cracks based on a ground penetrating radar.
Background
In western regions of China, particularly western windage sand areas and loess areas, ground cracks are one of geological environment problems caused by coal mining, which cause the problems of building deformation, underground pipeline damage, farmland damage, accelerated soil moisture evaporation, vegetation damage, water and soil loss and the like, bring great difficulty to mining area management workers, and are also an important link for land reclamation in mining areas. In order to study the influence of the ground fissure on the ecological environment, the underground morphology of the ground fissure is necessary for evaluating the risk and studying the development rule of the ground fissure. At present, a student injects gypsum slurry into a crack, and obtains the underground form of the crack through a ground penetrating radar technology after the gypsum is solidified, but the method can damage the structure of the crack and influence the subsequent development process of the crack, and the gypsum slurry is not uniformly poured, so that the filling inside the crack is insufficient, and the method cannot be used for carrying out the underground form research of dynamic cracks. Meanwhile, a student excavates a working section on one side of a ground crack, and then acquires the underground form of the crack through a ground penetrating radar technology transverse scanner.
Disclosure of Invention
In order to solve the problems, the invention provides a method for acquiring the multi-stage underground three-dimensional form of the dynamic fracture based on the ground penetrating radar, which has the characteristics of simple process, easy realization, no influence on the surface environment, extremely small required engineering amount, time saving and labor saving, does not influence the form of the fracture and the subsequent development process thereof, can acquire the multi-stage underground three-dimensional form in the dynamic fracture development process in a nondestructive manner, and has important significance for researching the development process and rule of the dynamic fracture underground form along with the change of time.
In order to achieve the purpose, the invention adopts the technical scheme that:
the invention provides a method for acquiring a multi-period underground three-dimensional form of a dynamic fracture based on a ground penetrating radar, which is characterized by comprising the following steps of:
step 1) laying a film on the surface of an area where a crack is located, and ensuring that the surface of the crack is completely covered by the film; pouring liquid on the film in the crack to enable the liquid to fill an area A formed by the film in the whole crack in an enclosing mode, and isolating the liquid from the crack through the film; aligning a ground penetrating radar to the crack at the near-surface position of the area A, and scanning at intervals from one end to the other end along the vertical direction of the trend of the row of cracks to obtain electromagnetic wave data of a plurality of sections of the underground form of the crack at the present stage; pumping out the liquid in the area A for storage and removing the film;
step 2) carrying out data processing on all the electromagnetic wave data acquired in the step 1) at the current period, and describing by a visual interpretation method to obtain the contour of the crack at each section at the current period;
step 3) drawing the contour of the fracture at each section in the current period obtained in the step 2) in three-dimensional image processing software, and constructing a three-dimensional model of the underground form of the fracture in the current period in an interpolation mode so as to obtain the underground three-dimensional form of the fracture in the current period and related data thereof;
and 4) continuously repeating the steps 1) to 3) according to the set frequency to obtain the multi-period underground three-dimensional form and related data thereof in the dynamic change process of the crack development.
Further, the liquid in step 1) is selected to avoid interference on radar wave energy.
Further, in the step 1), the ground penetrating radar is not in contact with liquid. Optionally, the ground penetrating radar is kept out of contact with the liquid by laying a thin plate on the film.
Further, the thin plate is selected to avoid interference with radar wave energy.
Further, in the step 1), when scanning is carried out at intervals from one end to the other end along the direction of the column seam, a marking line is made at intervals of 10cm in the direction perpendicular to the direction of the seam, and the ground penetrating radar is used for scanning from a position of 30cm on one side of the seam to a position of 30cm on the other side of the seam along the marking line.
Further, in step 2), the data processing includes background drying, one-dimensional filtering, wavelet transformation and gain operation.
Further, in step 3), the correlation data includes an average width, an average depth, a surface area, and a volume of the fracture.
The invention has the characteristics and beneficial effects that:
the invention provides a method for acquiring a multi-period underground three-dimensional form of a dynamic crack based on a ground penetrating radar, which solves the defect that the underground form of the crack can be damaged by injecting gypsum slurry into the crack, simultaneously avoids the defects of extremely large engineering amount and time and labor consumption when a ground penetrating radar technology is adopted for excavating a working section at one side of the crack, ensures that the crack can not permeate into soil by injecting liquid which can be extracted into the crack and can not damage the underground form of the crack and using a film, ensures that the ground penetrating radar can not contact with the liquid so that the ground penetrating radar can scan, has the characteristics of simple process, easy realization, no influence on the ground surface environment, extremely small required engineering amount, time and labor saving, simultaneously does not influence the form of the crack and the subsequent development process thereof, can nondestructively model the underground form of the dynamic crack, the method has the advantages that morphological parameters such as the average width, the average depth, the surface area and the volume of the crack are obtained, and the method has important significance for researching the development process and rule of the underground three-dimensional shape of the dynamic crack along with the change of time.
Drawings
FIGS. 1 (a) - (d) are respectively a three-dimensional view, a front view, a side view and a top view of a method of applying a film and filling a liquid in a ground fracture in step 1) of an embodiment of the method of the present invention;
fig. 2 (a) - (d) are respectively a three-dimensional view, a front view, a side view and a top view of the scanning mode of the ground penetrating radar in step 2) of the embodiment of the method of the present invention;
FIG. 3 is an electromagnetic wave image scanned by the ground penetrating radar in the embodiment of the present invention;
FIG. 4 is a cross-sectional view of a fracture subterranean formation depicted using visual interpretation in an embodiment of the present invention;
FIG. 5 is a cross-sectional view of a fracture subterranean formation of different scan profiles acquired in an embodiment of the present invention;
FIG. 6 is a schematic diagram of a multi-stage three-dimensional model during dynamic fracture development constructed according to an embodiment of the present invention.
Detailed Description
In order to make the objects, technical solutions and advantages of the present invention more apparent, the present invention is further described in detail below with reference to the accompanying drawings and embodiments. It should be understood that the detailed description and specific examples, while indicating the scope of the invention, are intended for purposes of illustration only and are not intended to limit the scope of the invention.
For better understanding of the present invention, an application example of the method for acquiring the multi-phase underground three-dimensional morphology of the dynamic fracture based on the ground penetrating radar is described in detail below.
The method for acquiring the multi-period underground three-dimensional form of the dynamic fracture based on the ground penetrating radar comprises the following steps:
step 1) scanning and acquiring electromagnetic wave data of crack area by using ground penetrating radar
Referring to fig. 1, a film which is large enough is laid on the surface of the area where the crack is located, and the surface of the crack is ensured to be completely covered by the film; pouring liquid on the film in the crack to make the liquid fill the area A formed by the film in the crack, and isolating the liquid from the crack through the film; the liquid should be selected from liquid with small interference degree on radar wave energy, such as water, oil, carbohydrate and the like, the thin film is solvent-resistant, friction-resistant, tear-resistant and liquid-sealed, the ground penetrating radar is aligned to the crack at the position close to the earth surface of the area A, scanning is carried out at intervals from one end to the other end along the trend of the line of the crack, and the vertical direction of the trend of the crack is used as a scanning path to obtain electromagnetic wave data of a plurality of sections of the underground form of the crack at the present time; then, the liquid in the area A is pumped out for storage and the film is removed.
In this embodiment, the specific scanning process is as follows:
referring to fig. 2, in order to ensure that the liquid in the crack does not contact the ground penetrating radar so as to avoid damaging the ground penetrating radar, a thin plate (except for isolating the liquid from the ground penetrating radar by the thin plate, other modes for isolating the liquid from the ground penetrating radar can be used) is laid on the plastic film in the scanning range, the thin plate can be a plastic plate or other polymeric thin plates which do not interfere with the scanning result of the ground penetrating radar, a marking line is made at every 10cm in the vertical direction of the crack trend, and the high-frequency ground penetrating radar is used for scanning from the 30cm position on one side of the crack to the 30cm position on the other side of the crack along the marking line, so that electromagnetic wave data of the underground form of the crack in the current period inclined to.
In the embodiment, the ground penetrating radar adopts a 1.6GHz high-frequency ground penetrating radar of China university of mining (Beijing) coal resources and a key laboratory of the safe mining country.
Step 2) processing all the electromagnetic wave data acquired in the step 1) at the current period, and drawing by a visual interpretation method to obtain the contour of the crack at each section at the current period
And (3) carrying out data processing on the image data of the region where the crack is located, which is obtained in the step 1), by using data processing software. In this embodiment, the data processing software adopts a GR radar processing and analyzing system developed by the university of mining china (beijing). And (2) respectively carrying out conventional background drying removal, one-dimensional filtering processing, wavelet transformation, gain and other operations on each image data acquired in the step 1) to obtain an electromagnetic wave image with an interference region as shown in fig. 3.
Referring to fig. 4, the contour of the fracture of the cross section is judged according to the interference area in the electromagnetic wave image, and the contour of the fracture at each cross section at present is depicted by a visual interpretation method.
Step 3) modeling based on SolidWorks three-dimensional morphology
Referring to fig. 5, the contour of the crack at each section in the current period obtained in the step 2) is drawn in a three-dimensional image processing software SolidWorks software according to the proportion of 1:1, and a three-dimensional model of the underground form in the current period of the crack is constructed in an interpolation mode, so that the underground three-dimensional form in the current period of the crack and related data thereof are obtained.
And 4) continuously repeating the steps 1) to 3) according to the set frequency to obtain the multi-period underground three-dimensional form and related data thereof in the dynamic change process of the crack development.
Referring to fig. 6, the present example obtained a three-phase subsurface three-dimensional morphology with the morphological parameters shown in table 1, with the volume accuracy being the fracture volume and the volume of fluid injected into the fracture.
TABLE 1 fracture morphology parameters
Figure BDA0002861242410000041
From table 1, it can be seen that the method of the present invention can obtain the volume of the fracture three-dimensional morphological modeling without damage, the error is about 20%, the main causes of the error are the error of the ground penetrating radar scanning image and the error in the three-dimensional modeling process, if more accurate depth information is required, the number of scanning layers can be encrypted or more accurate ground penetrating radar equipment and image processing software can be used.
The present invention and its embodiments have been described above schematically, without limitation, and what is shown in the drawings is only one of the embodiments of the present invention and is not actually limited thereto. Therefore, if the person skilled in the art receives the teaching, it is within the scope of the present invention to design the similar manner and embodiments without departing from the spirit of the invention.

Claims (8)

1. A method for acquiring a dynamic fracture multi-period underground three-dimensional form based on a ground penetrating radar is characterized by comprising the following steps:
step 1) laying a film on the surface of an area where a crack is located, and ensuring that the surface of the crack is completely covered by the film; pouring liquid on the film in the crack to enable the liquid to fill an area A formed by the film in the whole crack in an enclosing mode, and isolating the liquid from the crack through the film; aligning a ground penetrating radar to the crack at the near-surface position of the area A, and scanning at intervals from one end to the other end along the direction of the line of the crack to obtain electromagnetic wave data of a plurality of sections of the underground form of the crack at the present stage; pumping out the liquid in the area A for storage and removing the film;
step 2) carrying out data processing on all the electromagnetic wave data acquired in the step 1) at the current period, and describing by a visual interpretation method to obtain the contour of the crack at each section at the current period;
step 3) drawing the contour of the fracture at each section in the current period obtained in the step 2) in three-dimensional image processing software, and constructing a three-dimensional model of the underground form of the fracture in the current period in an interpolation mode so as to obtain the underground three-dimensional form of the fracture in the current period and related data thereof;
and 4) continuously repeating the steps 1) to 3) according to the set frequency to obtain the multi-period underground three-dimensional form and related data thereof in the dynamic change process of the crack development.
2. The method according to claim 1, wherein the liquid in step 1) is selected to avoid interference with radar wave energy.
3. The method according to claim 1, wherein in step 1), the ground penetrating radar is not in contact with the liquid.
4. A method according to claim 3, wherein the ground penetrating radar is kept out of contact with the liquid by laying a thin plate on the film.
5. The method of claim 4, wherein the sheet is selected to avoid interference with radar wave energy.
6. The method as claimed in claim 1, wherein in step 1), when scanning along the line of cracks from one end to the other end at intervals, a mark line is made at every 10cm interval in the direction perpendicular to the direction of the crack, and the ground penetrating radar is used for scanning from a position 30cm on one side of the crack to a position 30cm on the other side of the crack along the mark line.
7. The method according to claim 1, wherein in step 2), the data processing comprises background dessication, one-dimensional filtering processing, wavelet transformation and gain operation.
8. The method according to any one of claims 1 to 7, wherein in step 3) the correlation data comprises the average width, average depth, surface area and volume of the fracture.
CN202011566874.5A 2020-12-25 2020-12-25 Method for acquiring multi-period underground three-dimensional morphology of dynamic crack based on ground penetrating radar Active CN112799054B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011566874.5A CN112799054B (en) 2020-12-25 2020-12-25 Method for acquiring multi-period underground three-dimensional morphology of dynamic crack based on ground penetrating radar

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011566874.5A CN112799054B (en) 2020-12-25 2020-12-25 Method for acquiring multi-period underground three-dimensional morphology of dynamic crack based on ground penetrating radar

Publications (2)

Publication Number Publication Date
CN112799054A true CN112799054A (en) 2021-05-14
CN112799054B CN112799054B (en) 2024-03-22

Family

ID=75804948

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011566874.5A Active CN112799054B (en) 2020-12-25 2020-12-25 Method for acquiring multi-period underground three-dimensional morphology of dynamic crack based on ground penetrating radar

Country Status (1)

Country Link
CN (1) CN112799054B (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113433547A (en) * 2021-07-05 2021-09-24 长沙理工大学 Ground penetrating radar hidden crack offset imaging method, system, terminal and medium

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102590874A (en) * 2012-01-16 2012-07-18 中国矿业大学(北京) Method for detecting ground surface crack of upland coal-mining subsidence paddy field
CN102944195A (en) * 2012-11-28 2013-02-27 水利部交通运输部国家能源局南京水利科学研究院 Method for detecting depth of crack
CN104267449A (en) * 2014-09-15 2015-01-07 中国神华能源股份有限公司 Method for detecting coal mining earth surface crack patterns
CN108519596A (en) * 2018-03-05 2018-09-11 长沙理工大学 A kind of section of jurisdiction Hidden Fracture recognition methods based on match tracing and wavelet transformation
CN109696674A (en) * 2019-02-14 2019-04-30 中国矿业大学(北京) A method of ground fissure form is obtained based on Ground Penetrating Radar
US20190242246A1 (en) * 2018-02-08 2019-08-08 Saudi Arabian Oil Company Mapping fracture length using downhole ground penetrating radar

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102590874A (en) * 2012-01-16 2012-07-18 中国矿业大学(北京) Method for detecting ground surface crack of upland coal-mining subsidence paddy field
CN102944195A (en) * 2012-11-28 2013-02-27 水利部交通运输部国家能源局南京水利科学研究院 Method for detecting depth of crack
CN104267449A (en) * 2014-09-15 2015-01-07 中国神华能源股份有限公司 Method for detecting coal mining earth surface crack patterns
US20190242246A1 (en) * 2018-02-08 2019-08-08 Saudi Arabian Oil Company Mapping fracture length using downhole ground penetrating radar
CN108519596A (en) * 2018-03-05 2018-09-11 长沙理工大学 A kind of section of jurisdiction Hidden Fracture recognition methods based on match tracing and wavelet transformation
CN109696674A (en) * 2019-02-14 2019-04-30 中国矿业大学(北京) A method of ground fissure form is obtained based on Ground Penetrating Radar

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
兰天庆;胡泊洲;董文楠;张昕;: "砂煤岩互层水力裂缝扩展规律的数值模拟研究", 能源与环保, no. 10, 2 November 2018 (2018-11-02) *
张健;毕银丽;彭苏萍;: "采动地表裂缝三维形态探测方法及精度评价研究", 煤炭科学技术, vol. 48, no. 09, 15 September 2020 (2020-09-15) *
戴华阳;王祥;李军;郭俊廷;阎跃观;: "形态参数下地表采动裂缝三维建模及可视化方法", 测绘通报, no. 02, 25 February 2018 (2018-02-25) *
程远方;杨柳;吴百烈;董丙响;: "定向井压裂裂缝三维扩展形态的可视化仿真", 计算机仿真, no. 12, 15 December 2012 (2012-12-15) *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113433547A (en) * 2021-07-05 2021-09-24 长沙理工大学 Ground penetrating radar hidden crack offset imaging method, system, terminal and medium

Also Published As

Publication number Publication date
CN112799054B (en) 2024-03-22

Similar Documents

Publication Publication Date Title
WO2020248817A1 (en) Water supply pipeline leakage detection method based on three-dimensional image attribute analysis using ground penetrating radar
Acikgoz et al. Evaluation of the response of a vaulted masonry structure to differential settlements using point cloud data and limit analyses
US20220383475A1 (en) Method for identifying potential landslide hazard of reservoir bank based on rock mass degradation feature
CN102912780B (en) Sandy soil water-immersion testing method for loess collapsible deformation
CN105239611B (en) Determine the method that water stopping curtain seepage influences on surrounding enviroment below excavation of foundation pit face
WO2021008630A1 (en) Method, apparatus, and system for quantitative analysis of sealing strength of trap fault
Barraca et al. A case study of the use of GPR for rehabilitation of a classified Art Deco building: The InovaDomus house
CN103290836B (en) Vacuum decompression compacting and consolidating construction method for soft foundation treatment
Wu et al. The pixel crack reconstruction method: From fracture image to crack geological model for fracture evolution simulation
CN103091471B (en) A kind of model accelerating rinse fill consolidation and test method thereof
CN112799054B (en) Method for acquiring multi-period underground three-dimensional morphology of dynamic crack based on ground penetrating radar
CN110987765A (en) Rock mass microscopic fracture testing method based on three-dimensional digital speckle
CN109610525A (en) A kind of method of quick judgement water-stop curtain construction quality
CN109696674B (en) Method for acquiring ground fracture form based on ground penetrating radar
CN110363855B (en) Rock-fill dam transparentization modeling method
CN105975990B (en) A kind of fine expression of Reclaimed Soil section
CN105297753B (en) A kind of determination method of the hydrogeological parameter under water-stop curtain effect
CN106338257A (en) Method for quantitatively measuring embedding depth of proppant
Wang et al. Multi-dimensional size effects and representative elements for non-persistent fractured rock masses: A perspective of geometric parameter distribution
CN117313343A (en) Digital model acquisition method based on three-dimensional laser scanner blasting area
CN112554217A (en) Construction method of concrete impervious wall
Gillespie et al. A 3d model of complex A, la Venta, Mexico
CN113722663B (en) Method for evaluating underwater dredging and mud throwing of limited space
CN115577616A (en) Carbonatite fracture-cave earthquake depicting method and device based on deep learning
CN103628510A (en) Underground diaphragm wall model for test of underground diaphragm wall defect detection method

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