CN117372639B - Three-dimensional geological-geophysical modeling method combining sparse diaphysis section with auxiliary surface at any position - Google Patents
Three-dimensional geological-geophysical modeling method combining sparse diaphysis section with auxiliary surface at any position Download PDFInfo
- Publication number
- CN117372639B CN117372639B CN202311306202.4A CN202311306202A CN117372639B CN 117372639 B CN117372639 B CN 117372639B CN 202311306202 A CN202311306202 A CN 202311306202A CN 117372639 B CN117372639 B CN 117372639B
- Authority
- CN
- China
- Prior art keywords
- geological
- model
- data
- inversion
- dimensional
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 37
- 210000003275 diaphysis Anatomy 0.000 title claims description 3
- 238000012800 visualization Methods 0.000 claims abstract description 15
- 238000010276 construction Methods 0.000 claims abstract description 14
- 238000005553 drilling Methods 0.000 claims abstract description 11
- 238000011160 research Methods 0.000 claims abstract description 9
- 230000010365 information processing Effects 0.000 claims abstract description 6
- 230000033558 biomineral tissue development Effects 0.000 claims description 12
- 238000005065 mining Methods 0.000 claims description 11
- 230000008569 process Effects 0.000 claims description 11
- 239000011435 rock Substances 0.000 claims description 8
- 238000005516 engineering process Methods 0.000 claims description 7
- 230000000704 physical effect Effects 0.000 claims description 7
- 230000019771 cognition Effects 0.000 claims description 4
- 238000012545 processing Methods 0.000 claims description 4
- 238000004458 analytical method Methods 0.000 claims description 3
- 238000004364 calculation method Methods 0.000 claims description 3
- 230000008859 change Effects 0.000 claims description 3
- 238000012937 correction Methods 0.000 claims description 3
- 238000013461 design Methods 0.000 claims description 3
- 238000003708 edge detection Methods 0.000 claims description 3
- 238000000605 extraction Methods 0.000 claims description 3
- 238000001914 filtration Methods 0.000 claims description 3
- 230000005484 gravity Effects 0.000 claims description 3
- 238000000926 separation method Methods 0.000 claims description 3
- 230000006872 improvement Effects 0.000 claims description 2
- 229910052500 inorganic mineral Inorganic materials 0.000 abstract description 11
- 239000011707 mineral Substances 0.000 abstract description 11
- 230000001788 irregular Effects 0.000 abstract 1
- 230000015572 biosynthetic process Effects 0.000 description 3
- 230000007480 spreading Effects 0.000 description 3
- 230000000694 effects Effects 0.000 description 2
- 230000010354 integration Effects 0.000 description 2
- 238000011835 investigation Methods 0.000 description 2
- 230000002093 peripheral effect Effects 0.000 description 2
- 230000002159 abnormal effect Effects 0.000 description 1
- 230000005856 abnormality Effects 0.000 description 1
- 230000006978 adaptation Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000000452 restraining effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
- G06T17/05—Geographic models
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V3/00—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
- G01V3/38—Processing data, e.g. for analysis, for interpretation, for correction
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T19/00—Manipulating 3D models or images for computer graphics
Abstract
The invention relates to the technical field of mineral exploration, in particular to a three-dimensional geological-geophysical modeling method combining a sparse backbone profile with an auxiliary surface at any position, which mainly comprises the following steps: data arrangement, information processing interpretation, backbone geological section construction, three-dimensional geological model construction, three-dimensional constraint inversion and model visualization. The invention firstly utilizes comprehensive data to establish a backbone profile, controls a whole construction grid of a research area by using a sparse backbone cross profile, then uses the backbone profile as a basis to construct a three-dimensional geological model by an interpolation method, then converts the geological model into a physical model, performs three-dimensional grid inversion by taking geophysics, drilling information and geological data as constraints, and further modifies the model by gradually adding auxiliary surfaces (which can be irregular profiles) at any positions until the inversion result meets the set precision. The auxiliary surface can be a plane or a section with any inclination angle and trend.
Description
Technical Field
The invention relates to the technical field of mineral exploration, in particular to a three-dimensional geological-geophysical modeling method by combining a sparse backbone profile with an auxiliary surface at any position.
Background
With the continuous progress of mineral exploitation, china faces the crisis of resource exhaustion, so how to expand resource reserves in deep and peripheral areas of a working area (particularly a plurality of potential mineral collection areas) is one of the challenges faced inevitably, and the discovery of new mineral deposits in the deep and peripheral areas mainly depends on the knowledge of geological structures of a research area, including the shape and scale of the space of a geological body, the spatial spread of rock mass, a mineral control stratum and fracture structure, and the like, and in addition, the deep extension of a mineral formation zone is tracked, and a new mineral formation target area is searched in the deep layer, so that the method is also an important direction for realizing future deep exploration breakthrough. Therefore, the transparency of the ore collection area is the basis of further resource investigation, and the three-dimensional geological-geophysical modeling technology can trace back to the deep spreading of the ore-related geologic body, give out the control ore structure form or judge the important ore formation abnormal characteristics, and is an effective means for realizing the transparency of the research area.
The current three-dimensional geology-geophysical modeling method is more suitable for areas with relatively rich geological data and better stratum continuity, and for areas with sparse data and complex geological structures, actual geological conditions and deep geological structures are difficult to reflect, for example, areas with large amounts of coverage on the earth surface, too concentrated important data such as drilling holes and the like, or areas with extremely complex geological conditions are not reasonably and effectively solved.
Based on the reasons, a three-dimensional geological-geophysical modeling method combining a sparse backbone profile with an auxiliary surface at any position is designed, and the application of the method can indicate the direction for further prospecting work of a mining area and a mining collecting area, and provides technical support for mineral resource exploration, especially for exploration work of the mining collecting area and coverage area with sparse data.
Disclosure of Invention
The invention aims to overcome the defects of the prior art, provides a three-dimensional geological-geophysical modeling method combining a sparse backbone profile with an auxiliary surface at any position, indicates a direction for further prospecting work of a mining area/a mining collecting area, and provides technical support for mineral resource exploration, especially mineral collecting area, coverage area and deep prospecting work with sparse data.
The basic idea of the method is to build a sparse backbone profile by utilizing comprehensive data (geology, geophysics, geochemistry and the like) to control a grid of the whole structure of a research area, build a three-dimensional geological model by taking the backbone profile as a basis through an interpolation method, then convert the geological model into a physical model to perform heavy magnetic three-dimensional inversion, and modify the model by gradually adding auxiliary surfaces at any positions until the inversion result is satisfactory. The auxiliary surface can be a plane or a section with any inclination angle and trend.
The modeling flow mainly comprises 6 steps: data arrangement, information processing interpretation, backbone geological section construction, three-dimensional geological model construction, three-dimensional constraint inversion and model visualization.
S1, data arrangement:
the related data mainly comprise geological, drilling, physical, geophysical and geochemical data; the collected or collected information is divided into two types according to the requirements, one type is called constraint information, namely information capable of indirectly constraining the model form, the constraint information needs to be processed and converted into information capable of directly participating in inversion or indicating the model form, and the constraint information mainly comprises geophysical data, geochemical data and drilling data; the other type is called basic information, which is indispensable basic information used as a basic condition to participate in the inversion modeling process all the time and mainly comprises geological data and physical data;
s2, information processing interpretation:
the method mainly aims at constraint information to perform mining, enhancement and extraction of effective information. For geophysical data, performing edge detection, grid inversion and other processes on the basis of conventional gridding, filtering and bit field separation processes, and extracting effective information about fracture, rock mass and stratum geological units;
s3, constructing a backbone geological profile:
the section is deduced and drawn by using constraint information interpretation results with the overall construction grid of the investigation region as a target for planning the section position, so as to reflect knowledge of the spatial distribution of stratum, fracture, rock mass and ore body of the section passing through the region. After the preliminary construction is completed, vectorizing the profile, importing inversion software, correcting the profile by adopting a mature 2.5D gravity magnetic inversion modeling technology, and finally determining the backbone geological profile;
s4, constructing a three-dimensional geological model:
determining a modeling space range, adding terrain data, importing inversion corrected backbone geologic profile, actually measured and interpreted structural information and a surface geologic map into a three-dimensional space, and establishing a three-dimensional geologic model by utilizing a visualization and interpolation technology, wherein the model mainly characterizes an integral geologic structure grid of a research area, has a gap from an actual geologic condition in detail, and performs multi-constraint inversion work in the later stage, so that the process is essentially continuous correction and perfection of the geologic model;
s5, three-dimensional constraint inversion:
the method comprises the steps of taking the gravity-magnetic data separated by a potential field as inversion basic data, converting a three-dimensional geologic model into a physical model according to physical property analysis results, carrying out constraint inversion, wherein geologic body boundaries related to earth surface, drilling and backbone sections are not easy to change as known conditions, preferentially adjusting physical properties of geologic units in the inversion process, gradually adding auxiliary surfaces to modify the geologic body boundaries to fit measured data, and finally obtaining the three-dimensional model which accords with geological cognition and fits the measured gravity-magnetic data;
s6, model visualization:
importing the built model to a visualization platform; deep analyzing the spatial structure of the geologic body, extracting geological information, carrying out deep mineralization prediction, or analyzing the spatial relationship between geologic bodies related to mineralization according to the deep mineralization prediction, and establishing an mineralization mode; and carrying out mine design, reserve calculation and deep-side ore body prediction.
The visualization platform in S6 is Encom PA, voxel or a platform supporting formats such as vox/dxf/ts/stl/wrl.
The beneficial technical effects of the invention are as follows:
the modeling method and the modeling flow have good application effects on the geological model with wide construction area, relatively sparse data and complex geological structure area, and various existing data (geology, geophysics, geochemistry, physical properties and the like) are fully fused in a multi-source data integration mode and are used for restraining the three-dimensional model so as to establish a reliable complex geological structure model;
the form of gradually increasing the auxiliary surface can more reasonably and efficiently modify the model to fit the actual measurement abnormality, so that the establishment work of unnecessary sections is effectively reduced, the modeling time is shortened, and the working efficiency is improved;
the auxiliary surface with any direction, inclination angle, depth and length can make more full use of the known information, and can simultaneously give consideration to the changes of geological units in different directions so as to more accurately reflect the actual geological conditions;
the obtained modeling result can enable researchers to know the spreading forms of deep geological structures, fractures, rock masses and the like of the region more clearly, improves the cognition of basic geological conditions, can know the spatial spreading forms of underground structures, invaded bodies, strata and faults with depth of 3000m in a shallow range, indicates the direction for further prospecting work, and simultaneously provides technical support for deep and covered mining work.
Description of the drawings:
FIG. 1 is a schematic flow chart of the method for modeling of the present invention.
Detailed Description
Referring to fig. 1, the invention provides a three-dimensional geological-geophysical modeling method combining a sparse diaphyseal profile with an auxiliary surface at any position.
The method comprises the following steps:
s1, data arrangement:
the related data mainly comprises geological, drilling, physical, geophysical and geochemical data. The collected or collected information is divided into two types according to the requirements, one type is called constraint information, namely information capable of indirectly constraining the model form, the constraint information needs to be processed and converted into information capable of directly participating in inversion or indicating the model form, and the constraint information mainly comprises geophysical data, geochemical data and drilling data; the other type is called basic information, which is indispensable basic information used as a basic condition to participate in the inversion modeling process all the time and mainly comprises geological data and physical data;
s2, information processing interpretation:
the method mainly aims at constraint information to perform mining, enhancement and extraction of effective information. For geophysical data, performing edge detection and grid inversion processing on the basis of conventional gridding, filtering and bit field separation processing, and extracting effective information about fracture, rock mass and stratum geological units;
s3, constructing a backbone geological profile:
the overall construction grid of the research area is controlled to serve as a target planning section position, and the section is deduced and drawn by utilizing constraint information interpretation results so as to reflect the knowledge of the spatial distribution of stratum, fracture, rock mass and ore body of the section passing through the area; after the preliminary construction is completed, vectorizing the profile, importing inversion software, correcting the profile by adopting a mature 2.5D gravity magnetic inversion modeling technology, and finally determining the backbone geological profile;
s4, constructing a three-dimensional geological model:
determining a modeling space range, adding terrain data, importing inversion corrected backbone geological section, actually measured and interpreted structural information and a surface geological map into a three-dimensional space, and establishing a three-dimensional geological model by utilizing a visualization and interpolation technology; at the moment, the model mainly characterizes the overall geological structure grid of the research area, the difference is still remained between the details and the actual geological condition, and the constraint inversion work is carried out for multiple times in the later period, so that the continuous correction and improvement process of the geological model is essential;
s5, three-dimensional constraint inversion:
the method comprises the steps of taking the gravity-magnetic data separated by a potential field as inversion basic data, converting a three-dimensional geologic model into a physical model according to physical property analysis results, carrying out constraint inversion, wherein geologic body boundaries related to earth surface, drilling and backbone sections are not easy to change as known conditions, preferentially adjusting physical properties of geologic units in the inversion process, gradually adding auxiliary surfaces to modify the geologic body boundaries to fit measured data, and finally obtaining the three-dimensional model which accords with geological cognition and fits the measured gravity-magnetic data;
s6, model visualization:
importing the built model to a visualization platform; deep analyzing the spatial structure of the geologic body, extracting geological information, carrying out deep mineralization prediction, or analyzing the spatial relationship between geologic bodies related to mineralization according to the deep mineralization prediction, and establishing an mineralization mode; and carrying out mine design, reserve calculation and deep-side ore body prediction.
The visualization platform in S6 is Encom PA, voxel or a platform supporting formats such as vox/dxf/ts/stl/wrl.
The above is only a preferred embodiment of the present invention, only for helping to understand the method and the core idea of the present application, and the protection scope of the present invention is not limited to the above examples, and all technical solutions belonging to the concept of the present invention belong to the protection scope of the present invention. It should be noted that modifications and adaptations to the present invention may occur to one skilled in the art without departing from the principles of the present invention and are intended to be within the scope of the present invention.
The invention fundamentally solves the problem that the three-dimensional geology-geophysical modeling work of areas with sparse and uneven geological data, complex geological conditions and the like in the prior art has no reasonable and effective solution, especially the mining work of a coverage area or a deep part, fully utilizes various existing data through a multi-source data integration mode, and efficiently and reliably establishes a complex geological structure model, thereby more accurately simulating the spatial structure form of the geological body and reflecting the geological problems, and is particularly suitable for mining areas with larger area, relatively sparse geological data and complex geological conditions. The method can be used for knowing the spatial distribution form of underground structures, invaded bodies, strata and faults with depth of 3000m in a shallow range, indicating the direction for further prospecting work and providing technical support for deep and covered area prospecting work.
Claims (2)
1. A three-dimensional geology-geophysical modeling method combining a sparse diaphysis section with an auxiliary surface at any position is characterized in that the modeling process mainly comprises 6 steps: data arrangement, information processing interpretation, backbone geological section construction, three-dimensional geological model construction, three-dimensional constraint inversion and model visualization;
s1, data arrangement:
the related data mainly comprise geological, drilling, physical, geophysical and geochemical data; the collected or collected information is divided into two types according to the requirements, one type is called constraint information, namely information capable of indirectly constraining the model form, the constraint information needs to be processed and converted into information capable of directly participating in inversion or indicating the model form, and the constraint information mainly comprises geophysical data, geochemical data and drilling data;
the other type is called basic information, which is indispensable basic information used as a basic condition to participate in the inversion modeling process all the time and mainly comprises geological data and physical data;
s2, information processing interpretation:
the method mainly aims at constraint information to perform mining, enhancement and extraction of effective information; for geophysical data, performing edge detection and grid inversion processing on the basis of conventional gridding, filtering and bit field separation processing, and extracting effective information about fracture, rock mass and stratum geological units;
s3, constructing a backbone geological profile:
the overall construction grid of the research area is controlled to serve as a target planning section position, and the section is deduced and drawn by utilizing constraint information interpretation results so as to reflect the knowledge of the spatial distribution of stratum, fracture, rock mass and ore body of the section passing through the area; after the preliminary construction is completed, vectorizing the profile, importing inversion software, correcting the profile by adopting a mature 2.5D gravity magnetic inversion modeling technology, and finally determining the backbone geological profile;
s4, constructing a three-dimensional geological model:
determining a modeling space range, adding terrain data, importing inversion corrected backbone geological section, actually measured and interpreted structural information and a surface geological map into a three-dimensional space, and establishing a three-dimensional geological model by utilizing a visualization and interpolation technology; at the moment, the model mainly characterizes the overall geological structure grid of the research area, the difference is still remained between the details and the actual geological condition, and the constraint inversion work is carried out for multiple times in the later period, so that the continuous correction and improvement process of the geological model is essential;
s5, three-dimensional constraint inversion:
the method comprises the steps of taking the gravity-magnetic data separated by a potential field as inversion basic data, converting a three-dimensional geologic model into a physical model according to physical property analysis results, carrying out constraint inversion, wherein geologic body boundaries related to earth surface, drilling and backbone sections are not easy to change as known conditions, preferentially adjusting physical properties of geologic units in the inversion process, gradually adding auxiliary surfaces to modify the geologic body boundaries to fit measured data, and finally obtaining the three-dimensional model which accords with geological cognition and fits the measured gravity-magnetic data;
s6, model visualization:
importing the built model to a visualization platform; deep analyzing the spatial structure of the geologic body, extracting geological information, carrying out deep mineralization prediction, or analyzing the spatial relationship between geologic bodies related to mineralization according to the deep mineralization prediction, and establishing an mineralization mode; and carrying out mine design, reserve calculation and deep-side ore body prediction.
2. The three-dimensional geologic-geophysical modeling method of a sparse diaphyseal cross-section in combination with arbitrary position aid according to claim 1, wherein the visualization platform in S6 is Encom PA, voxel or a platform supporting the format of vox/dxf/ts/stl/wrl.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202311306202.4A CN117372639B (en) | 2023-10-10 | 2023-10-10 | Three-dimensional geological-geophysical modeling method combining sparse diaphysis section with auxiliary surface at any position |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202311306202.4A CN117372639B (en) | 2023-10-10 | 2023-10-10 | Three-dimensional geological-geophysical modeling method combining sparse diaphysis section with auxiliary surface at any position |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN117372639A CN117372639A (en) | 2024-01-09 |
| CN117372639B true CN117372639B (en) | 2024-03-26 |
Family
ID=89397537
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202311306202.4A Active CN117372639B (en) | 2023-10-10 | 2023-10-10 | Three-dimensional geological-geophysical modeling method combining sparse diaphysis section with auxiliary surface at any position |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN117372639B (en) |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104866653A (en) * | 2015-04-29 | 2015-08-26 | 中国地质科学院矿产资源研究所 | Method for acquiring underground three-dimensional density structure |
| CN115082634A (en) * | 2022-05-05 | 2022-09-20 | 中国矿业大学(北京) | Three-dimensional geological structure method for weathered layer granite rare earth ore |
| WO2022257930A1 (en) * | 2021-06-07 | 2022-12-15 | 中国地质大学（北京） | Method and apparatus for simulating mineralization in mineralization research area |
| CN115932988A (en) * | 2022-12-20 | 2023-04-07 | 中国地质科学院 | Heavy magnetic data three-dimensional inversion modeling method and system constrained by prior information |
-
2023
- 2023-10-10 CN CN202311306202.4A patent/CN117372639B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104866653A (en) * | 2015-04-29 | 2015-08-26 | 中国地质科学院矿产资源研究所 | Method for acquiring underground three-dimensional density structure |
| WO2022257930A1 (en) * | 2021-06-07 | 2022-12-15 | 中国地质大学（北京） | Method and apparatus for simulating mineralization in mineralization research area |
| CN115082634A (en) * | 2022-05-05 | 2022-09-20 | 中国矿业大学(北京) | Three-dimensional geological structure method for weathered layer granite rare earth ore |
| CN115932988A (en) * | 2022-12-20 | 2023-04-07 | 中国地质科学院 | Heavy magnetic data three-dimensional inversion modeling method and system constrained by prior information |
Also Published As
| Publication number | Publication date |
|---|---|
| CN117372639A (en) | 2024-01-09 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN111750822B (en) | Coal mining induced overlying strata and surface subsidence collaborative dynamic prediction method | |
| CN103837908A (en) | Rapid prospecting positioning method applicable to hidden sandstone-type uranium mine | |
| CN105785471A (en) | Impact danger evaluation method of mine pre-exploiting coal seam | |
| CN107065019A (en) | Applied to road disaster and the 3 D electromagnetic imaging device and application method that collapse detection | |
| CN102877828A (en) | CT (Computed Tomography) imaging method of three-dimensional multi-well combined well land | |
| CN112255693A (en) | Goaf filling space detection method under coal mine fully mechanized caving coal mining process condition | |
| CN113536414B (en) | Rock slope stability analysis method, system and medium based on three-dimensional modeling | |
| CN109285221A (en) | Deposit three-dimensional Geological Modeling | |
| CN117372639B (en) | Three-dimensional geological-geophysical modeling method combining sparse diaphysis section with auxiliary surface at any position | |
| CN111505733A (en) | Multi-source heterogeneous information fusion analysis and dynamic prediction method for occurrence of special geologic body of underground engineering fractured rock mass | |
| CN114814978B (en) | Granite area tungsten-tin ore exploration method based on multiple depth scales | |
| CN110376659A (en) | A kind of tunnel geological forecast device and its system | |
| CN115932988A (en) | Heavy magnetic data three-dimensional inversion modeling method and system constrained by prior information | |
| CN115185015A (en) | Deep lithium beryllium ore exploration method | |
| Wang et al. | Multi-source data integration and multi-scale modeling framework for progressive prediction of complex geological interfaces in tunneling | |
| CN113887046A (en) | Coal mine tunnel modeling method based on three-dimensional geologic body | |
| CN112734076A (en) | Large-scale positioning prediction method for hard rock type uranium ore resources | |
| Ji et al. | An automated method to build 3D multi-scale geological models for engineering sedimentary layers with stratum lenses | |
| CN112696181A (en) | Fracturing prediction method and device based on coupling relation between ground stress and natural fractures | |
| CN105824042A (en) | Maximum longitudinal distance design method based on optimal lighting energy | |
| CN113267827B (en) | Earthquake prediction method and device using petroleum earthquake and drilling data | |
| CN114814982B (en) | Method for predicting favorable ore-forming part of granite uranium ore | |
| CN115828636B (en) | Anti-channeling construction parameter optimization method for shale gas well group fracturing | |
| CN114594518B (en) | Fine stratum contrast method for complex fault blocks in later development period based on well-seismic alternation | |
| Liu* et al. | Application of microseismic monitoring technology in low production and low efficiency reservoir reconstruction |
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 |